Next Article in Journal
Effects of Prey Distribution and Heterospecific Interactions on the Functional Response of Harmonia axyridis and Aphidius gifuensis to Myzus persicae
Next Article in Special Issue
Wing Geometric Morphometrics as a Tool for the Identification of Culex Subgenus Mosquitoes of Culex (Diptera: Culicidae)
Previous Article in Journal
Transcriptome Analysis of the Oriental Fruit Fly Bactrocera dorsalis Early Embryos
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phylogeny of Anopheles (Kerteszia) (Diptera: Culicidae) Using Mitochondrial Genes

by
Karin Kirchgatter
1,*,
Lilian de Oliveira Guimarães
1,
Henrry Hugo Yañez Trujillano
2,
Fernando Rafael Arias
3,
Abraham Germán Cáceres
4,5,
Ana Maria Ribeiro de Castro Duarte
1,
Rosely dos Santos Malafronte
6,
Rosa Maria Tubaki
7 and
Maria Anice Mureb Sallum
8
1
Laboratório de Bioquímica e Biologia Molecular, Superintendência de Controle de Endemias/Instituto de Medicina Tropical de São Paulo, Universidade de São Paulo, SP 01027-000, Brazil
2
Area de Vigilancia y Control Vectorial, Dirección de Salud Ambiental, Dirección Ejecutiva de Salud Ambiental, Dirección Regional de Salud Cusco, Cusco 08003, Peru
3
Laboratorio de Referencia Regional de Salud Pública, Dirección Regional de Salud Puno, Puno 21002, Peru
4
Seccion de Entomología, Instituto de Medicina Tropical “Daniel A. Carrión” y Departamento Académico de Microbiología Médica, Facultad de Medicina, Universidad Nacional Mayor de San Marcos, Lima 15081, Peru
5
Laboratorio de Entomologia, Instituto Nacional de Salud, Lima 15064, Peru
6
Laboratório de Protozoologia, Instituto de Medicina Tropical de São Paulo, Universidade de São Paulo, São Paulo, SP 05403-000, Brazil
7
Laboratório de Entomologia Médica, Superintendência de Controle de Endemias, São Paulo, SP 01027-000, Brazil
8
Departamento de Epidemiologia, Faculdade de Saúde Pública, Universidade de São Paulo, São Paulo, SP 01246-904, Brazil
*
Author to whom correspondence should be addressed.
Insects 2020, 11(5), 324; https://doi.org/10.3390/insects11050324
Submission received: 3 May 2020 / Revised: 19 May 2020 / Accepted: 21 May 2020 / Published: 24 May 2020
(This article belongs to the Special Issue Vector-Borne Diseases in a Changing World)

Abstract

:
Identification of mosquito species is necessary for determining the entomological components of malaria transmission, but it can be difficult in morphologically similar species. DNA sequences are largely used as an additional tool for species recognition, including those that belong to species complexes. Kerteszia mosquitoes are vectors of human and simian malaria in the Neotropical Region, but there are few DNA sequences of Kerteszia species in public databases. In order to provide relevant information about diversity and improve knowledge in taxonomy of Kerteszia species in Peru, we sequenced part of the mitochondrial genome, including the cytochrome c oxidase I (COI) barcode region. Phylogenetic analyses structured all species of mosquitoes collected in Peru into a single clade, separate from the Brazilian species. The Peruvian clade was composed of two lineages, encompassing sequences from Anopheles (Kerteszia) boliviensis and Anopheles (Kerteszia) pholidotus. An. pholidotus sequences were recorded for the first time in Peru, whereas An. boliviensis sequences were for the first time published in the GenBank database. Sequences generated from specimens morphologically identified as Anopheles (Kerteszia) cruzii clustered into three separate clades according to the collection localities of Serra do Mar, Serra da Mantiqueira, and Serra da Cantareira, confirming An. cruzii as a species complex, composed of at least three putative species.

1. Introduction

Malaria is a mosquito-borne tropical disease that remains an important public health problem in some tropical and subtropical countries. In South America, the majority of cases occur in areas of the Amazon rainforest, and are mainly concentrated in Venezuela, Brazil, and Peru, which registered 30%, 24%, and 19% of the American malaria cases in 2015, respectively [1].
In Peru, malaria transmission occurs throughout the year in rural areas below 2000 m altitude, with the highest risk in the department of Loreto in the Amazon Region, which registered 96% of the 55,218 cases reported in the country in 2017 [2]. In 2017, the department of Cusco, where mosquitoes from this study were collected, registered an increase of 41% in the malaria cases in relation to 2016, being the fifth department in Peru in number of cases [2]. In Lima and areas at high altitude (including the Inca Trail, Machu Picchu, and Lake Titicaca) the malaria risk is low. In 2019, there was a decrease in cases of malaria, with 56% fewer reported cases compared to the same period in 2018. However, an outbreak of 34 cases of malaria due to Plasmodium vivax was reported in Tumbes Department, in the northern part of the country, which included 22 autochthonous cases and 12 cases imported from Venezuela; this outbreak highlights the risk of re-introduction of malaria in a historically endemic area where transmission had been interrupted in recent years [3].
Plasmodium vivax is responsible for 80% of the Peruvian malaria cases [4]. However, as in other endemic countries, the incidence of Plasmodium malariae may be underestimated due to difficulties in microscopic differentiation from P. vivax, and its presence in coinfections. Plasmodium malariae has been identified in the vast majority of blood samples from Amerindian ethnic groups from Peru, including coinfection with P. vivax [5]. In areas of the Peruvian Amazon, the primary vector is Anopheles darlingi. This species accounts for >85% of the anopheline mosquitoes that were observed to be blood-feeding on humans, and it is responsible for malaria transmission, particularly in frontier areas [6,7,8,9]. In Iquitos, Peruvian Amazon, a high proportion of An. darlingi was found blood-feeding on humans and avians, showing the vector a generalist mosquito [10].
In Brazil, 99% of malaria cases registered are from the Amazon Region, where An. darlingi is the primary vector, and other species can participate as secondary vectors in specific areas [11]. Outside the Amazon region, malaria endemic areas are mainly related to the Atlantic forest biome where Kerteszia species are associated with phytotelma habitats, such as those in the Bromeliaceae family, where they use the axils as a larval habitat [12]. In the Atlantic Forest, Plasmodium infection is usually asymptomatic, caused by sub-patent levels of P. vivax or P. malariae [13], and Anopheles (Kerteszia) cruzii and Anopheles (Kerteszia) bellator are their vectors, with An. cruzii involved in the transmission of both human and simian malaria [14,15,16,17].
Anopheles species that use bromeliads phytotelmata as larval habitats have been recognized as a distinct group by Knab (1913). The relationship between the subgenera Kerteszia and Nyssorhynchus has been hypothesized by Root (1922). Following this, Christophers (1924) and Edwards (1932) included Kerteszia in the subgenus Nyssorhynchus, but Dyar (1928) and Komp (1937) treat them as a distinct subgenera [18]. Recently, Foster et al. [19] proposed Kertezia as a genus on the basis of the robust results of mitogenome phylogeny of the subfamily Anophelinae. There are 12 valid Kerteszia species: Anopheles (Kerteszia) auyantepuiensis, Anopheles (Kerteszia) bambusicolus, Anopheles (Kerteszia) bellator, Anopheles (Kerteszia) boliviensis, Anopheles (Kerteszia) cruzii, Anopheles (Kerteszia) gonzalezrinconesi, Anopheles (Kerteszia) homunculus, Anopheles (Kerteszia) laneanus, Anopheles (Kerteszia) lepidotus, Anopheles (Kerteszia) neivai, Anopheles (Kerteszia) pholidotus, and Anopheles (Kerteszia) rollai [20].
Kerteszia mosquitoes are widely distributed from Mexico to southern Brazil [18,21,22]. In Peru, there are records of An. bambusicolus, An. boliviensis, An. cruzii, An. homunculus, An. laneanus, An. lepidotus, An. neivai, and An. pholidotus (ref. in [17]). In Brazil, the occurrences of An. bambusicolus (in areas of the interior of the Atlantic Forest), An. bellator, An. cruzii and An. homunculus (in coastal and mountainous areas of the Serra do Mar in the Atlantic Forest Region), An. laneanus (in areas of Serra da Mantiqueira), and An. neivai (in the Amazon Region) have been registered (ref. in [17]). An. boliviensis was registered in Brazil [23] but the absence of recent records makes its identification questionable [20].
Anopheles homunculus, An. cruzii, An. bellator, and An. laneanus form a well-defined species group, which is always easily separated from the other Kerteszia species [18]. Anopheles homunculus and An. cruzii are morphologically similar on the basis of the external characteristics of the females. However, they can be easily identified by characteristics of male genitalia, as well as fourth-instar larva and pupa [24,25]. Anopheles cruzii has been demonstrated to be a species complex on the basis of cytological differences on the banding patterns of polythene chromosomes [26]. Recently, Oliveira et al. [27], using the mitogenome, showed that the specimens from Serra da Cantareira, São Paulo, Brazil, may belong to distinct lineage that is sister to An. cruzii.
Differences in the ecology and behavior of species can affect the transmission of malaria and the success of the methods used in the control, such as insecticides and insecticide-impregnated bed nets [28]. Therefore, the characterization of species composition of a local population for the design of vector control projects is of vital importance. DNA-based methods have been widely used for species identification. Mitochondrial genes have been broadly employed in studies focusing on taxonomy and phylogeny of the Culicidae and Anophelinae subfamilies. The cytochrome c oxidase subunit I (COI) gene of the mitogenome has been largely used to identify taxonomic groups at species level [29,30]. There are COI sequences available in the GenBank and BOLD databases, including sequences generated for the Mosquito Barcoding Initiative (MBI) (http://barcodeoflife.ning.com/group/mosquitobarcoding).
In the current study, we initially used in the phylogenetic analysis a fragment of mitochondrial cytochrome c oxidase I (COI) that represents the barcode region [31,32]. The aim was to analyze the haplotype diversity and verify the genetic variability of specimens collected in Peru. This mitochondrial region is the main source of reference sequences, and has been successfully used to identify mosquitoes [25,29,33,34] and to evaluate intraspecific variability, especially in Neotropical Anopheles [35,36,37]. This analysis included mosquitoes from Peru, which were previously identified as An. cruzii, An. boliviensis, An. Laneanus, and An. lepidotus using morphological characteristics. Lastly, we generated two more partial sequences (another COI region and ND4, NADH dehydrogenase subunit 4 gene), totalizing 10% of the whole mtDNA (mitochondrial DNA) (≈1500 bp), and compared with reference sequences from a previously published complete mitochondrial DNA genome, to verify whether An. cruzii is a species complex.

2. Materials and Methods

2.1. Mosquito Sampling and Handling

Females were collected in different localities and dates (Table 1 and Table 2). In Brazil, mosquitoes were collected using CDC (Centers for Disease Control) light traps with CO2 (dry ice) for 12 h, beginning at twilight, or with Shannon trap for 3 h, from 6:00 p.m. to 9:00 p.m. Peruvian mosquitoes were collected with Shannon traps, from 6:00 p.m. to 9:00 p.m. Mosquitoes were killed with chloroform vapor and transported to the laboratory for identification, using morphological taxonomy key by Zavortink [18], Consoli and Lourenço-de-Oliveira [38], and Forattini [12]. Mosquitoes were individually stored at −20 °C in 1.5 mL plastic tubes sealed with parafilm.

2.2. Genomic DNA Extraction and PCR Amplification of Mitochondrial Gene Fragments

Mosquitoes were crushed with disposable tips in lysis buffer, and genomic DNA (gDNA) was obtained using PureLink Genomic DNA Purification Kit (Invitrogen), according to the manufacturer’s instructions.
A fragment of 710 base pairs (bp) of the barcode region of the mitochondrial COI gene was generated by PCR amplification using the primer pair LCO1490 and HCO2198 [39], following the protocol proposed by Ruiz et al. [40]. A fragment of ≈600 bp of COI gene (gene region used in phylogeny) was obtained as described [41] with the primers COIF (5′-GGA TTA TTA GGA TTT ATT GT-3′) and COIR (5′-GCA AAT AAT GAA ATT GTT CT-3′). A fragment of ≈400 bp of the mitochondrial ND4 gene was generated using the primers ND4+ (5′-GTD YAT TTA TGA TTR CCT AA-3′) and ND4-(5′-CTT CGD CTT CCW ADW CGT TC-3′) and the PCR program, as described [42].

2.3. Sequencing, Alignment, and Phylogenetic Analysis

Amplified fragments were sequenced directly in both directions by BigDye Terminator v3.0 Cycle Sequencing Kit in ABI Genetic Analyzer (Applied Biosystems®, Foster City, CA, USA), using the corresponding flanking primers. The COI and ND4 sequences obtained from newly sequenced mosquitoes were deposited in the GenBank database (COI: MH589352-MH589421; ND4: MH560305-MH560342). Sequences were aligned with the reference sequences (Table 3) using Clustal X version 1.81 [43]. Alignments were then inspected and edited in MEGA version 5.1 [44].
Two methods of phylogenetic reconstruction were implemented using the Bayesian approach. In the first round of analysis, 624 bp of the COI barcoding region was employed, due the higher availability of sequences of this COI gene region for the Anopheles diversity in the GenBank database, from which all the sequence data of reference Kerteszia species were downloaded. In the second round of analysis, phylogenetic relationship was investigated using ND4 and COI genes. COI and ND4 datasets were concatenated, after trimming the primers off from the sequences—including 338 bp of the ND4 gene, 511 bp of the COI, and 624 bp of the barcode region of the COI, generating a total of 1473 bp. In this analysis only reference sequences from complete mtDNA were used, avoiding an incomplete final matrix and missing data. All sequence accession numbers used in both analyses are provided in Table 1, Table 2 and Table 3 and added into phylogenetic trees.
The phylogenetic reconstructions were performed using the Bayesian approach implemented in MrBayes v3.2.0 [45]. Bayesian inferences were run with two Markov Chain Monte Carlo searches of 3 million generations, each with sampling of 1 in 300 trees. After a burn-in of 25%, the remaining 15,002 trees were used to generate a 50% majority-rule consensus tree. The results of analyses were visualized using FigTree version 1.4.0 [46]. The topologies were rooted using sequences obtained for species of the subgenera Anopheles or Nyssorhynchus (Table 2 and Table 3).

3. Results

The results of the analyses using only DNA sequences of the COI barcoding region are shown in Figure 1. This approach was used with the aim of increasing the number of sequences in the phylogenetic tree, as there are no complete mitochondrial sequences available for some Kerteszia species such as An. neivai, An. pholidotus, and An. lepidotus. The sequence alignment composed of 85 sequences with 624 bp length—9 from Peru, 26 from Brazil, and 50 references from GenBank—had 168 variable sites. The Bayesian phylogenetic tree (Figure 1) shows that the sequences from nine specimens from Peru grouped together in a strongly supported clade. The Peruvian species clade contains two major groups. One group encompasses sequences from specimens identified as An. boliviensis, An. cruzii, An. laneanus, and An. lepidotus. The second group contains sequences from An. lepidotus that clustered with An. pholidotus from Venezuela.
For the second round of Bayesian analysis, approximately 1.5 kb (10% of the whole mtDNA) was employed, including COI and ND4 genes of 39 specimens of An. (Kerteszia) species from Peru and Brazil. Twelve new sequences were blasted with sequences available in GenBank for other Kerteszia species to verify identification and potential presence of stop codon in each sequence. The alignment generated 1473 bp with 371 variable sites. Sequences of Anopheles (Anopheles) species were used as outgroup. All clades recovered in the Bayesian phylogenetic analysis were strongly supported (Bayesian Posterior Probability—BPP = 1). All the Anopheles (Kerteszia) species were grouped in a clade that was separated into two clades: one containing all Peruvian species, and the second clade being formed of species collected in Brazil (Figure 2). The Peruvian species clade was split in two minor groups. Species from Brazil were separated in two clades: one clade composed of An. bellator and An. homunculus, and a second clade composed of An. laneanus and An. cruzii, with four mosquitoes morphologically identified as An. bellator (one of which was collected in 1946). Sequences of An. cruzii from Brazil were separated into clades in accordance with the geographical localities of Serra do Mar, Serra da Mantiqueira and Serra da Cantareira (Figure 2 and Figure 3).

4. Discussion

Recently, the number of publications using mitochondrial sequences for inferring phylogenetic relationships within the Kerteszia has increased [47,48,49,50,51,52,53]. However, there are few studies focusing on populations of An. cruzii [19,27,54], with no publication regarding Kerteszia species that occur in Peru. In order to provide relevant information about diversity and improve knowledge in taxonomy of Kerteszia species from Peru, we sequenced two mitochondrial genes (COI and ND4), including the DNA barcode region.
The analysis of genetic differentiation of COI and ND4 genes indicated that there are at least two Kerteszia species of that occur in Peru. Preliminary identification of specimens of An. cruzii and An. lepidotus using morphological characters of the female was not confirmed by phylogenetic analyses of mitochondrial genes. In fact, specimens from Peru that were identified as An. cruzii had scales on the abdominal terga, indicating that they belonged to other species. In the Bayesian topology, those specimens of An. cruzii with scales on the abdomen terga grouped with most of the mosquitoes from Peru, which were identified as An. boliviensis. In addition, one specimen morphologically identified as An. lepidotus clustered with An. pholidotus, with an identity of 98% in the COI barcoding region. The clade encompassing An. pholidotus sequences was recovered and was well separated from An. lepidotus. Moreover, a Peruvian specimen preliminary identified as An. laneanus clustered with An. boliviensis. For accurate identification of An. boliviensis, An. pholidotus, and An. lepidotus using morphology, it is important to examine characters of the male and fourth-instar larva in addition to female characters [18]. Furthermore, until 2012 [25], the only morphological feature that had been proposed to differentiate females of An. pholidotus from An. lepidotus was the size of the scales on proximal tergites, being moderately wide or broad for An. lepidotus and predominantly narrow to moderately wide for An. pholidotus [18]. This characteristic was shown vague and difficult to interpret. In fact, in other countries of South America also, such as Colombia, molecular analyses showed that the species previously believed to be An. lepidotus corresponded to An. pholidotus [47].
Chromosomal banding pattern of the ovarian polytene chromosomes had already suggested that there are three putative species under the name An. cruzii, termed as A, B, and C [26,55]. More recently, studies employing DNA sequences from nuclear genes have corroborated this hypothesis [56,57,58,59]. Our analysis, using DNA sequences from mitochondrial genes, also supported the hypotheses of An. cruzii representing a complex of morphologically similar species. The sequences employed in this study that were generated from specimens identified as An. cruzii were separated in accordance to their geographical origin. Sequence with the GenBank code KU551284 was recovered in a clade that is sister to An. laneanus and yet separate from other specimens of An. cruzii. This KU551284 sequence was obtained from a mosquito collected in the Parque da Serra da Cantareira, situated in the mountain range north of the municipality of São Paulo [27]. The KU551285 sequence was also sorted separately from all the other cruzii sequences and was obtained from a mosquito collected in Itatiaia National Park, in the Mantiqueira mountain range, between the states of Rio de Janeiro and Minas Gerais [27]. All the other An. cruzii sequences were grouped in a branch composed by sequences from mosquitoes collected only in the Serra do Mar region, a mountainous chain that extends by approximately 1500 km along the coast, from Rio de Janeiro state to the north of Santa Catarina state. Indeed, at least two genetically distinct An. cruzii lineages have been shown to occur in the mountains covered by the Atlantic Forest in South-East Brazil [60]. Here, through another methodology, we confirmed these two populations and showed that a third, coming from another mountain range, the Serra da Cantareira, also separates from the other two in the phylogenetic trees. These three populations of An. cruzii were consistently separated regardless of the methodology used (analysis of mitochondrial or nuclear genes). Moreover, in our phylogenetic trees, the addition of another COI gene region as well as another mitochondrial gene (ND4) has not changed the topology of the tree, showing that the most important information is contained in the barcoding region, confirming that the COI gene is more useful for identifying sibling and cryptic species, as has been demonstrated [51].
We also demonstrated that mosquitoes morphologically identified as An. bellator (CA1, B1, B2, and B3) are likely An. cruzii. In fact, adult females of these species are morphologically similar and identification can be problematic when the adults were damaged during field collections [18]. In addition to the placement of these mosquitoes on the phylogenetic trees, their grouping based on the geographical region of the field collections confirms that An. cruzii is a species complex.

5. Conclusions

The analysis of genetic differentiation of COI and ND4 genes indicated that there are at least two Kerteszia species occurring in Peru: An. boliviensis and An. pholidotus. The occurrence of these species in Peru is herein registered for the first-time, using COI DNA sequences. Moreover, we confirmed the occurrence of two subpopulations of An. cruzii in the Atlantic rainforest and showed a third population in the mountain range of the Parque da Serra da Cantareira. This population is likely a distinct species that clustered separately from the remaining two lineages of An. cruzii in the phylogenetic trees. One lineage encompasses specimens from Serra do Mar on the Atlantic coast, and the second lineage is formed by An. cruzii from Serra da Mantiqueira. Besides their placement in the phylogenetic tree, the grouping of An. cruzii sequences reflects the geographical origin of the adults (Serra da Cantareira, Serra do Mar, and Serra da Mantiqueira). Thus, these findings add additional support for the An. cruzii species complex.

Author Contributions

Conceptualization, K.K.; formal analysis, K.K. and M.A.M.S.; investigation, L.d.O.G.; resources, K.K., H.H.Y.T., F.R.A., A.G.C., A.M.R.d.C.D., R.d.S.M., and R.M.T.; writing—original draft, K.K. and L.d.O.G.; writing—review and editing, K.K. and M.A.M.S. All authors have read and agreed to the published version of the manuscript.

Funding

We thank the Vicerrectorado for Research of the Universidad Nacional Mayor de San Marcos for their sponsorship (research code: 100101171). M.A.M.S. is financially supported by Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), grant number 2014/26229-7, and Conselho Nacional de Pesquisa (CNPq), grant number 301877/2016-5. A.M.R.d.C.D. is supported by grant FAPESP number 2008/52016-0. L.O.G. is supported by a postdoctoral fellowship (FAPESP 2018/16232-1). K.K. is a CNPq research fellow, grant number 308678/2018-4.

Acknowledgments

We thank Bruno Mathias for his help in editing the figures.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization. World Malaria Report; WHO: Geneva, Switzerland, 2016. [Google Scholar]
  2. Centro Nacional de Epidemiología, Prevención y Control de Enfermedades—MINSA, Perú. Available online: http://www.dge.gob.pe/portal/docs/vigilancia/sala/2017/SE51/malaria.pdf (accessed on 22 May 2019).
  3. Pan American Health Organization; World Health Organization. Epidemiological Update: Malaria in the Americas; PAHO: Washington, DC, USA; WHO: Geneva, Switzerland, 2019. [Google Scholar]
  4. Recht, J.; Siqueira, A.M.; Monteiro, W.M.; Herrera, S.M.; Herrera, S.; Lacerda, M.V.G. Malaria in Brazil, Colombia, Peru and Venezuela: Current challenges in malaria control and elimination. Malar. J. 2017, 16, 273. [Google Scholar] [CrossRef]
  5. Sulzer, A.J.; Cantella, R.; Colichon, A.; Gleason, N.N.; Walls, K.W. A focus of hyperendemic Plasmodium malariaeP. vivax with no P. falciparum in a primitive population in the Peruvian Amazon jungle. Bull. World Health Organ. 1975, 52, 273–278. [Google Scholar]
  6. Aramburu, G.J.; Asayag, C.R.; Witzig, R. Malaria reemergence in the peruvian amazon region. Emerg. Infect. Dis. 1999, 5, 209–215. [Google Scholar]
  7. Manguin, S.R.; Wilkerson, R.C.; Conn, J.E.; Rubio-Palis, Y.; Danoff-Burg, J.A.; Roberts, D.R. Population structure of the primary malaria vector in South America, Anopheles darlingi, using isoenzyme, ITS2, RAPD, and morphological markers. Am. J. Trop. Med. Hyg. 1999, 60, 364–376. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Flores, C.; Fernández, R.; Escobedo, K.S.; Vela, Q.; Schoeler, G.B. Natural Plasmodium infections in Anopheles darlingi and Anopheles benarrochi (Diptera: Culicidae) from Eastern Peru. J. Med. Entomol. 2004, 41, 489–494. [Google Scholar] [CrossRef] [PubMed]
  9. Reinbold-Wasson, D.D.; Sardelis, M.R.; Jones, J.W.; Watts, D.M.; Fernandez, R.; Carbajal, F.; Pecor, J.E.; Calampa, C.; Klein, T.A.; Turell, M.J. Determinants of Anopheles seasonal distribution patterns across a forest to periurban gradient near Iquitos, Peru. Am. J. Trop. Med. Hyg. 2012, 86, 459–463. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Moreno, M.; Saavedra, M.P.; Bickersmith, S.A.; Prussing, C.; Michalski, A.; Tong Rios, C.; Vinetz, J.M.; Conn, J.E. Intensive trapping of blood-fed Anopheles darlingi in Amazonian Peru reveals unexpectedly high proportions of avian blood-meals. PLoS Negl. Trop. Dis. 2017, 11, e0005337. [Google Scholar] [CrossRef] [PubMed]
  11. Oliveira-Ferreira, J.; Lacerda, M.V.; Brasil, P.; Ladislau, J.L.; Tauil, P.L.; Daniel-Ribeiro, C.T. Malaria in Brazil: An overview. Malar. J. 2010, 9, 115. [Google Scholar] [CrossRef] [Green Version]
  12. Forattini, O. P Culicidologia Médica, vol. 2: Identificação, Biologia, Epidemiologia; Edusp: São Paulo, Brazil, 2002. [Google Scholar]
  13. Carlos, B.C.; Rona, L.D.P.; Christophides, G.K.; Souza-Neto, J.A. A comprehensive analysis of malaria transmission in Brazil. Pathog. Glob. Health 2019, 113, 1–13. [Google Scholar] [CrossRef] [Green Version]
  14. Deane, L.M.; Ferreira-Neto, J.A.; Deane, S.P.; Silveira, I.P. Anopheles (Kerteszia) cruzii, a natural vector of the monkey malaria parasites, Plasmodium simium and Plasmodium Brasilianum. Trans. R. Soc. Trop. Med. Hyg. 1970, 64, 647. [Google Scholar] [CrossRef]
  15. Deane, L.M.; Ferreira-Neto, J.A.; Lima, M.M. The vertical dispersion of Anopheles (Kerteszia) cruzii in a forest in southern Brazil suggests that human cases of malaria of simian origin might be expected. Mem. Inst. Oswaldo Cruz 1984, 79, 461–463. [Google Scholar] [CrossRef] [PubMed]
  16. Demari-Silva, B.; Laporta, G.Z.; Oliveira, T.; Sallum, M. Plasmodium infection in Kerteszia cruzii (Diptera: Culicidae) in the Atlantic tropical rain forest, southeastern Brazil. Infect. Genet. Evol. 2020, 78, 104061. [Google Scholar] [CrossRef] [PubMed]
  17. Multini, L.C.; Wilke, A.B.B.; Marrelli, M.T. Neotropical Anopheles (Kerteszia) mosquitoes associated with bromeliad-malaria transmission in a changing world. Acta Trop. 2020, 205, 105413. [Google Scholar] [CrossRef] [PubMed]
  18. Zavortink, T.J. Mosquito studies XXIX. A review of the subgenus Kerteszia of Anopheles. Contr. Amer. Entomol. Inst. 1973, 9, 1–54. [Google Scholar]
  19. Foster, P.G.; de Oliveira, T.M.P.; Bergo, E.S.; Conn, J.E.; Sant’Ana, D.C.; Nagaki, S.S.; Nihei, S.; Lamas, C.E.; González, C.; Moreira, C.C.; et al. Phylogeny of Anophelinae using mitochondrial protein coding genes. R. Soc. Open Sci. 2017, 4, 170758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Collucci, E.; Sallum, M.A.M. Phylogenetic analysis of the subgenus Kerteszia of Anopheles (Diptera: Culicidae: Anophelinae) based on morphological characters. Insect Syst. Evol. 2003, 34, 361–372. [Google Scholar] [CrossRef]
  21. Marrelli, M.T.; Malafronte, R.S.; Sallum, M.A.; Natal, D. Kerteszia subgenus of Anopheles associated with the Brazilian Atlantic rainforest: Current knowledge and future challenges. Malar. J. 2007, 6, 127. [Google Scholar] [CrossRef]
  22. Walter Reed Biosystematics Unit (WRBU). Systematic Catalog of Culicidae. 2019. Available online: http://www.mosquitocatalog.org/taxon_descr.aspx?ID=16082 (accessed on 31 August 2019).
  23. Guimarães, J.H. Systematic Database of Diptera of the Americas South of the United States (Family Culicidae); Plêiade/FAPESP: São Paulo, Brazil, 1997. [Google Scholar]
  24. Sallum, M.A.M.; Dos Santos, C.L.S.; Wilkerson, R.C. Studies on Anopheles (Kerteszia) homunculus Komp (Diptera: Culicidae). Zootaxa 2009, 2299, 1–18. [Google Scholar] [CrossRef]
  25. Harrison, B.A.; Ruiz-Lopez, F.; Falero, G.C.; Savage, H.M.; Pecor, J.E.; Wilkerson, R.C. Anopheles (Kerteszia) lepidotus (Diptera: Culicidae), not the malaria vector we thought it was: Revised male and female morphology; larva, pupa, and male genitalia characters; and molecular verification. Zootaxa 2012, 3218, 1–17. [Google Scholar] [CrossRef]
  26. Ramirez, C.C.; Dessen, E.M. Chromosome differentiated populations of Anopheles cruzii: Evidence for a third sibling species. Genetica 2000, 108, 73–80. [Google Scholar] [CrossRef]
  27. Oliveira, T.M.; Foster, P.G.; Bergo, E.S.; Nagaki, S.S.; Sanabani, S.S.; Marinotti, O.; Marinotti, P.N.; Sallum, M.A. Mitochondrial Genomes of Anopheles (Kerteszia) (Diptera: Culicidae) From the Atlantic Forest, Brazil. J. Med. Entomol. 2016, 53, 790–797. [Google Scholar] [CrossRef] [PubMed]
  28. Multini, L.C.; Marrelli, M.T.; Beier, J.C.; Wilke, A.B.B. Increasing Complexity Threatens the Elimination of Extra-Amazonian Malaria in Brazil. Trends Parasitol. 2019, 35, 383–387. [Google Scholar] [CrossRef] [PubMed]
  29. Ruiz-Lopez, F.; Wilkerson, R.C.; Conn, J.E.; McKeon, S.N.; Levin, D.M.; Quiñones, M.L.; Póvoa, M.M.; Linton, Y.M. DNA barcoding reveals both known and novel taxa in the albitarsis group (Anopheles: Nyssorhynchus) of Neotropical malaria vectors. Parasit. Vectors 2012, 5, 44. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Bourke, B.P.; Oliveira, T.P.; Suesdek, L.; Bergo, E.S.; Sallum, M.A. A multi-locus approach to barcoding in the Anopheles strodei subgroup (Diptera: Culicidae). Parasit. Vectors 2013, 6, 111. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Hebert, P.D.N.; Cywinska, A.; Ball, S.L.; deWaard, J.R. Biological identifications through DNA barcodes. Proc. Biol. Sci. 2003, 270, 313–321. [Google Scholar] [CrossRef] [Green Version]
  32. Hebert, P.D.N.; Ratnasingham, S.; deWaard, J.R. Barcoding animal life: Cytochrome c oxidase subunit 1 divergences among closely related species. Proc. Biol. Sci. 2003, 270, S96–S99. [Google Scholar] [CrossRef] [Green Version]
  33. Cywinska, A.; Hunter, F.F.; Hebert, P.D.N. Identifying Canadian mosquito species through DNA barcodes. Med. Vet. Entomol. 2006, 20, 413–424. [Google Scholar] [CrossRef]
  34. Kumar, N.P.; Rajavel, A.R.; Natarajan, R.; Jambulingam, P. DNA barcodes can distinguish species of Indian mosquitoes (Diptera: Culicidae). J. Med. Entomol. 2007, 44, 1–7. [Google Scholar] [CrossRef]
  35. Mirabello, L.; Conn, J.E. Molecular population genetics of the malaria vector Anopheles darlingi in Central and South America. Heredity 2006, 96, 311–321. [Google Scholar] [CrossRef]
  36. Gutiérrez, L.; Orrego, L.M.; Gómez, G.F.; López, A.; Luckhart, S.; Conn, J.E.; Correa, M.M. A new mtDNA COI gene lineage closely related to Anopheles janconnae of the Albitarsis complex in the Caribbean region of Colombia. Mem. Inst. Oswaldo Cruz 2010, 105, 1019–1025. [Google Scholar] [CrossRef]
  37. Jaramillo, L.M.; Gutiérrez, L.A.; Luckhart, S.; Conn, J.E.; Correa, M.M. Molecular evidence for a single taxon, Anopheles nuneztovari s.l., from two endemic malaria regions in Colombia. Mem. Inst. Oswaldo Cruz 2011, 106, 1017–1023. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Consoli, R.A.G.B.; Lourenço-de-Oliveira, R. Principais Mosquitos de Importância Sanitária No Brasil; Editora Fiocruz: Rio de Janeiro, Brazil, 1994. [Google Scholar]
  39. Folmer, O.; Black, M.; Hoeh, W.; Lutz, R.; Vrijenhoek, R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotechnol. 1994, 3, 294–299. [Google Scholar] [PubMed]
  40. Ruiz, F.; Linton, Y.M.; Ponsonby, D.J.; Conn, J.E.; Herrera, M.; Quiñones, M.L.; Vélez, I.D.; Wilkerson, R.C. Molecular comparison of topotypic specimens confirms Anopheles (Nyssorhynchus) dunhami Causey (Diptera: Culicidae) in the Colombian Amazon. Mem. Inst. Oswaldo Cruz 2010, 105, 899–903. [Google Scholar] [CrossRef] [PubMed]
  41. Sallum, M.A.M.; Schultz, T.R.; Foster, P.G.; Aronstein, K.; Wirtz, R.A.; Wilkerson, R.C. Phylogeny of Anophelinae (Diptera: Culicidae) based on nuclear ribosomal and mitochondrial DNA sequences. Syst. Entomol. 2002, 27, 361–382. [Google Scholar] [CrossRef] [Green Version]
  42. Gorrochotegui-Escalante, N.; Munoz, M.L.; Fernandez-Salas, I.; Beaty, B.J.; Black, W.C. 4th. Genetic isolation by distance among Aedes aegypti populations along the northeastern coast of Mexico. Am. J. Trop. Med. Hyg. 2000, 62, 200–209. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Thompson, J.D.; Gibson, T.J.; Plewniak, F.; Jeanmougin, F.; Higgins, D.G. The CLUSTAL_X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25, 4876–4882. [Google Scholar] [CrossRef] [Green Version]
  44. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [Green Version]
  45. Huelsenbeck, J.P.; Ronquist, F. MRBAYES: Bayesian inference of phylogeny. Bioinformatics 2001, 17, 754–755. [Google Scholar] [CrossRef] [Green Version]
  46. Rambaut, A. FigTree: Tree Figure Drawing Tool Version 1.4.0, Institute of Evolutionary Biology, University of Edinburgh. 2006–2012. Available online: http://tree.bio.ed.ac.uk/software/figtree (accessed on 22 May 2019).
  47. Escovar, J.E.; González, R.; Quiñones, M.L.; Wilkerson, R.C.; Ruiz, F.; Harrison, B.A. Morphology of the larvae, male genitalia and DNA sequences of Anopheles (Kerteszia) pholidotus (Diptera: Culicidae) from Colombia. Mem. Inst. Oswaldo Cruz 2014, 109, 473–479. [Google Scholar] [CrossRef] [Green Version]
  48. Freitas, L.A.; Russo, C.A.; Voloch, C.M.; Mutaquiha, O.C.; Marques, L.P.; Schrago, C.G. Diversification of the Genus Anopheles and a Neotropical Clade from the Late Cretaceous. PLoS ONE 2015, 10, e0134462. [Google Scholar] [CrossRef] [Green Version]
  49. Hua, Y.Q.; Ding, Y.R.; Yan, Z.T.; Si, F.L.; Luo, Q.C.; Chen, B. The complete mitochondrial genome of Anopheles minimus (Diptera: Culicidae) and the phylogenetics of known Anopheles mitogenomes. Insect Sci. 2016, 23, 353–365. [Google Scholar] [CrossRef] [PubMed]
  50. López-Rubio, A.; Suaza-Vasco, J.; Marcet, P.L.; Ruíz-Molina, N.; Cáceres, L.; Porter, C.; Uribe, S. Use of DNA barcoding to distinguish the malaria vector Anopheles neivai in Colombia. Zootaxa 2016, 4175, 377–389. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Wang, G.; Li, C.; Zheng, W.; Song, F.; Guo, X.; Wu, Z.; Luo, P.; Yang, Y.; He, L.; Zhao, T. An evaluation of the suitability of COI and COII gene variation for reconstructing the phylogeny of, and identifying cryptic species in, anopheline mosquitoes (Diptera Culicidae). Mitochondrial Dna A Dna Mapp. Seq. Anal. 2017, 28, 769–777. [Google Scholar] [CrossRef] [PubMed]
  52. Chen, K.; Wang, Y.; Li, X.Y.; Peng, H.; Ma, Y.J. Sequencing and analysis of the complete mitochondrial genome in Anopheles sinensis (Diptera: Culicidae). Infect. Dis. Poverty 2017, 6, 149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. López-Rubio, A.; Suaza, J.D.; Porter, C.; Uribe, S.; Bedoya, G.; Vélez, I.D. Phylogenetic signal at the Cytb-SertRNA-IG1-ND1 mitochondrial region in Anopheles (Kerteszia) neivai Howard, Dyar & Knab, 1913. Biomedica 2017, 37, 143–154. [Google Scholar] [PubMed] [Green Version]
  54. Lorenz, C.; Patané, J.S.; Suesdek, L. Morphogenetic characterisation, date of divergence, and evolutionary relationships of malaria vectors Anopheles cruzii and Anopheles homunculus. Infect. Genet. Evol. 2015, 35, 144–152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Ramirez, C.C.L.; Dessen, E.M.B. The polytene chromosomes of the mosquito Anopheles bellator compared with those of Anopheles cruzii. Braz. J. Genet. 1996, 19, 555–558. [Google Scholar] [CrossRef] [Green Version]
  56. Rona, L.D.; Carvalho-Pinto, C.J.; Gentile, C.; Grisard, E.C.; Peixoto, A.A. Assessing the molecular divergence between Anopheles (Kerteszia) cruzii populations from Brazil using the timeless gene: Further evidence of a species complex. Malar. J. 2009, 8, 60. [Google Scholar] [CrossRef] [Green Version]
  57. Rona, L.D.; Carvalho-Pinto, C.J.; Mazzoni, C.J.; Peixoto, A.A. Estimation of divergence time between two sibling species of the Anopheles (Kerteszia) cruzii complex using a multilocus approach. BMC Evol. Biol. 2010, 10, 91. [Google Scholar] [CrossRef] [Green Version]
  58. Rona, L.D.; Carvalho-Pinto, C.J.; Peixoto, A.A. Evidence for the occurrence of two sympatric sibling species within the Anopheles (Kerteszia) cruzii complex in southeast Brazil and the detection of asymmetric introgression between them using a multilocus analysis. BMC Evol. Biol. 2013, 13, 207. [Google Scholar] [CrossRef] [Green Version]
  59. Rona, L.D.; Carvalho-Pinto, C.J.; Peixoto, A.A. Molecular evidence for the occurrence of a new sibling species within the Anopheles (Kerteszia) cruzii complex in south-east Brazil. Malar. J. 2010, 9, 33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. de Rezende Dias, G.; Fujii, T.T.S.; Fogel, B.F.; Lourenço-de-Oliveira, R.; Silva-do-Nascimento, T.F.; Pitaluga, A.N.; Carvalho-Pinto, C.J.; Carvalho, A.B.; Peixoto, A.A.; Rona, L.D.P. Cryptic diversity in an Atlantic Forest malaria vector from the mountains of South-East Brazil. Parasit. Vectors 2018, 11, 36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Phylogeny of the Kerteszia species using cytochrome c oxidase I (COI) barcoding region by Bayesian inference tree. Posterior probability values are shown to each clade. Yellow, An. cruzii clade; grey, An. laneanus clade; white, An. bellator; blue, An. lepidotus; dark green, An. homunculus; orange, An. neivai; pink, An. pholidotus; green, An. boliviensis; dark blue, An. (Anopheles) species; dark grey, An. (Nyssorhynchus) species. Specimens reclassified are highlighted in bold (see Table 1 and Table 2).
Figure 1. Phylogeny of the Kerteszia species using cytochrome c oxidase I (COI) barcoding region by Bayesian inference tree. Posterior probability values are shown to each clade. Yellow, An. cruzii clade; grey, An. laneanus clade; white, An. bellator; blue, An. lepidotus; dark green, An. homunculus; orange, An. neivai; pink, An. pholidotus; green, An. boliviensis; dark blue, An. (Anopheles) species; dark grey, An. (Nyssorhynchus) species. Specimens reclassified are highlighted in bold (see Table 1 and Table 2).
Insects 11 00324 g001
Figure 2. Phylogeny of the Kerteszia species using COI and NADH dehydrogenase subunit 4 (ND4) genes by Bayesian inference tree. Posterior probability values are shown to each clade. An. cruzii from is Serra da Cantareira represented in red. An. cruzii from Serra da Mantiqueira is represented in blue. An. cruzii from Serra do Mar is represented as a yellow box. Specimens reclassified are highlighted in bold (see Table 1 and Table 2).
Figure 2. Phylogeny of the Kerteszia species using COI and NADH dehydrogenase subunit 4 (ND4) genes by Bayesian inference tree. Posterior probability values are shown to each clade. An. cruzii from is Serra da Cantareira represented in red. An. cruzii from Serra da Mantiqueira is represented in blue. An. cruzii from Serra do Mar is represented as a yellow box. Specimens reclassified are highlighted in bold (see Table 1 and Table 2).
Insects 11 00324 g002
Figure 3. (a) Origin localities map for An. cruzii sequences from Peru and Brazil used in this study. (b) Localities of Itatiaia/RJ (Serra da Mantiqueira, blue dot); Parque Estadual da Cantareira/SP (Serra da Cantareira, red dot); Cananéia/SP, Caraguatatuba/SP, Itanhaém/SP, and São Paulo/SP (all in Serra do Mar, yellow dots) are shown.
Figure 3. (a) Origin localities map for An. cruzii sequences from Peru and Brazil used in this study. (b) Localities of Itatiaia/RJ (Serra da Mantiqueira, blue dot); Parque Estadual da Cantareira/SP (Serra da Cantareira, red dot); Cananéia/SP, Caraguatatuba/SP, Itanhaém/SP, and São Paulo/SP (all in Serra do Mar, yellow dots) are shown.
Insects 11 00324 g003
Table 1. Main features of the specimens from Peru used in this study and GenBank accession numbers.
Table 1. Main features of the specimens from Peru used in this study and GenBank accession numbers.
Mosquito IDOriginCollection YearSpecies (According to Morphological Identification)Species (According to Phylogeny Analysis)COI * (COIF/COIR Primers)ND4 *COI * (LCO/HCO Primers)
1AMiraflores, Yanatile, Calca, Cusco, Peru2007Anopheles (Kerteszia) boliviensisAnopheles (Kerteszia) boliviensisMH589387MH560305--
1BMiraflores, Yanatile, Calca, Cusco, Peru2007Anopheles (Kerteszia) boliviensisAnopheles (Kerteszia) boliviensisMH589388MH560306MH589377
2ACorimayo, Yanatile, Calca, Cusco, Peru2008Anopheles (Kerteszia) boliviensisAnopheles (Kerteszia) boliviensisMH589389MH560307--
2BCorimayo, Yanatile, Calca, Cusco, Peru2008Anopheles (Kerteszia) boliviensisAnopheles (Kerteszia) boliviensisMH589390MH560308MH589378
3AColca, Yanatile, Calca, Cusco, Peru2008Anopheles (Kerteszia) boliviensisAnopheles (Kerteszia) boliviensisMH589391MH560309MH589379
3BColca, Yanatile Calca, Cusco, Peru2008Anopheles (Kerteszia) boliviensisAnopheles (Kerteszia) boliviensisMH589392MH560310MH589380
4AMiraflores, Yanatile, Calca, Cusco, Peru2008Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) boliviensisMH589393MH560311MH589381
4BMiraflores, Yanatile, Calca, Cusco, Peru2008Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) boliviensisMH589394MH563012MH589382
5AColca, Yanatile, Calca, Cusco, Peru2008Anopheles (Kerteszia) lepidotusAnopheles (Kerteszia) pholidotusMH589395MH563013--
5BColca, Yanatile, Calca, Cusco, Peru2008Anopheles (Kerteszia) lepidotusAnopheles (Kerteszia) pholidotusMH589396MH560314--
6AMiraflores, Yanatile, Calca, Cusco, Peru2008Anopheles (Kerteszia) lepidotusAnopheles (Kerteszia) boliviensisMH589397MH560315MH589383
6BMiraflores, Yanatile, Calca, Cusco, Peru2008Anopheles (Kerteszia) lepidotusAnopheles (Kerteszia) pholidotusMH589398MH563016MH589384
LANSelva Alegre, Alto Inambari, Sandia, Puno, Peru2000Anopheles (Kerteszia) laneanusAnopheles (Kerteszia) boliviensisMH589409---MH589367
* COI, cytochrome c oxidase subunit I gene; ND4, NADH dehydrogenase subunit 4 gene.
Table 2. Main features of the specimens from Brazil used in this study and GenBank accession numbers.
Table 2. Main features of the specimens from Brazil used in this study and GenBank accession numbers.
Mosquito IDOriginCollection YearSpecies (According Morphological Identification)Species (According Phylogeny Analysis)COI * (COIF/COIR Primers)ND4 *COI * (LCO/HCO Primers)
CA1Caraguatatuba, SP, Brazil1946Anopheles (Kerteszia) bellatorAnopheles (Kerteszia) cruziiMH589406MH560328MH589361
B1Marsilac, São Paulo, SP, Brazil2009Anopheles (Kerteszia) bellatorAnopheles (Kerteszia) cruzii--MH560317MH589362
B2Marsilac, São Paulo, SP, Brazil2009Anopheles (Kerteszia) bellatorAnopheles (Kerteszia) cruziiMH589399MH560318MH589363
B3Marsilac, São Paulo, SP, Brazil2009Anopheles (Kerteszia) bellatorAnopheles (Kerteszia) cruzii--MH560319MH589364
PAC1Parelheiros, São Paulo, SP, Brazil2002Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589412MH560333MH589352
PAC2Parelheiros, São Paulo, SP, Brazil2002Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589414MH560335MH589353
PAC3Parelheiros, São Paulo, SP, Brazil2002Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589415MH560336MH589354
PAC4Parelheiros, São Paulo, SP, Brazil2002Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589416MH560337MH589355
PAC5Parelheiros, São Paulo, SP, Brazil2002Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589417MH560338MH589356
PAC6Parelheiros, São Paulo, SP, Brazil2002Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589418--MH589357
PAC7Parelheiros, São Paulo, SP, Brazil2002Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589419MH560339MH589358
PAC8Parelheiros, São Paulo, SP, Brazil2002Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589420MH560340MH589359
PAC9Parelheiros, São Paulo, SP, Brazil2002Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589421MH560341MH589360
PAC10XParelheiros, São Paulo, SP, Brazil2002Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589413MH560334--
C1Marsilac, São Paulo, SP, Brazil2009Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruzii--MH560320MH589369
C2Marsilac, São Paulo, SP, Brazil2009Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589400MH560321MH589370
C3Marsilac, São Paulo, SP, Brazil2009Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruzii--MH560322MH589371
C4Marsilac, São Paulo, SP, Brazil2009Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589401MH560323MH589372
C5Marsilac, São Paulo, SP, Brazil2009Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589402MH560324MH589373
C6Marsilac, São Paulo, SP, Brazil2009Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589403MH560325MH589374
C7Marsilac, São Paulo, SP, Brazil2009Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589404MH560326MH589375
C8Marsilac, São Paulo, SP, Brazil2009Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589405MH560327MH589376
L1Campos do Jordão, SP, Brazil 1997Anopheles (Kerteszia) laneanusAnopheles (Kerteszia) laneanusMH589407MH560329MH589365
L3Campos do Jordão, SP, Brazil 1997Anopheles (Kerteszia) laneanusAnopheles (Kerteszia) laneanusMH589408MH560330MH589366
M95Itanhaém, SP, Brazil2010Anopheles (Kerteszia) cruziiAnopheles (Kerteszia) cruziiMH589411MH560332MH589386
M83Itanhaém, SP, Brazil2010Anopheles (Anopheles) maculipes/pseudomaculipesAnopheles (Anopheles) intermediusMH589410MH560331MH589385
PAN1Parelheiros, São Paulo, SP, Brazil2002Anopheles (Nyssorhynchus) sp.Anopheles (Nysssorhynchus) strodei--MH560342MH589368
* COI, cytochrome c oxidase subunit I gene; ND4, NADH dehydrogenase subunit 4 gene.
Table 3. GenBank accession numbers of the reference mosquitoes used in the study.
Table 3. GenBank accession numbers of the reference mosquitoes used in the study.
SpeciesMosquito IDOrigin# GenBank
An. lepidotus51917879Napo/EcuadorJQ041282 *
An. lepidotus51917894Napo/EcuadorJQ041284 *
An. pholidotus51917595Tachira/VenezuelaJQ041287 *
An. pholidotus51917607Tachira/VenezuelaJQ041288 *
An. homunculusBA22_33 Camacan, BA/BrazilJQ291234 *
An. homunculusES10_2 Santa Teresa, ES/BrazilJQ291235 *
An. homunculusRS20_21 Maquiné, RS/BrazilJQ291236 *
An. homunculusRS20_33 Maquiné, RS/BrazilJQ291237 *
An. homunculusRS30_56 Maquiné, RS/BrazilJQ291238 *
An. homunculusSP23_1 Cananéia, SP/BrazilJQ291240 *
An. homunculusST26Cananéia, SP/BrazilJQ291242 *
An. lepidotusCCDB-10764-A9Orellana/EcuadorKF670991 *
An. cruzii???Cananéia/SPKJ701506
An. cruzii104Tapiraí, SP/BrazilKM507293 *
An. cruzii110Tapiraí, SP/BrazilKM507294 *
An. cruzii112Tapiraí, SP/BrazilKM507295 *
An. cruzii127Tapiraí, SP/BrazilKM507296 *
An. cruzii192Tapiraí, SP/BrazilKM507297 *
An. cruzii198Tapiraí, SP/BrazilKM507298 *
An. cruzii199Tapiraí, SP/BrazilKM507299 *
An. cruzii200Tapiraí, SP/BrazilKM507300 *
An. cruzii217Tapiraí, SP/BrazilKM507301 *
An. cruzii240Tapiraí, SP/BrazilKM507302 *
An. cruzii241Tapiraí, SP/BrazilKM507303 *
An. cruzii242Tapiraí, SP/BrazilKM507304 *
An. cruzii255Tapiraí, SP/BrazilKM507305 *
An. cruzii257Tapiraí, SP/BrazilKM507306 *
An. cruzii258Tapiraí, SP/BrazilKM507307 *
An. cruzii259Tapiraí, SP/BrazilKM507308 *
An. cruzii260Tapiraí, SP/BrazilKM507309 *
An. cruzii261Tapiraí, SP/BrazilKM507310 *
An. cruzii262Tapiraí, SP/BrazilKM507311 *
An. cruzii267Tapiraí, SP/BrazilKM507312 *
An. cruzii269Tapiraí, SP/BrazilKM507314 *
An. homunculusBA22_32Camacan/BAKU551283
An. cruziiPEC_2_7São Paulo/SPKU551284
An. cruziiRJ03_2Itatiaia/RJKU551285
An. cruziiRS32_11_7Maquiné/RSKU551286
An. bellatorSP24_3_1Cananéia/SPKU551287
An. laneanusSP52-103Pindamonhangaba/SPKU551288
An. cruziiST16Cananéia/SPKU551289
An. neivaiMB10165Petit-Saut/French GuyanaMF172273 *
An. neivaiMB10252Petit-Saut/French GuyanaMF172274 *
An. neivaiMB10253Petit-Saut/French GuyanaMF172275 *
An. neivaiMB10254Petit-Saut/French GuyanaMF172276 *
An. homunculusBA22_31Camacan/BAMF381605
An. laneanusSP52-17Pindamonhangaba/SPMF381613
An. homunculusRS30_158 Maquiné/RSMF381697
An. (An.) intermediusSP02_8_1Pariquera-Açú/SPMF381700
An. (Nyr.) strodeiSP104_18_1Pindamonhangaba/SPMF381724 *
An. homunculusST19 Cananéia/SPMF381729
* Only for COI barcoding region analysis. # GenBank accession number (https://www.ncbi.nlm.nih.gov/genbank/).

Share and Cite

MDPI and ACS Style

Kirchgatter, K.; de Oliveira Guimarães, L.; Hugo Yañez Trujillano, H.; Rafael Arias, F.; Cáceres, A.G.; de Castro Duarte, A.M.R.; dos Santos Malafronte, R.; Tubaki, R.M.; Mureb Sallum, M.A. Phylogeny of Anopheles (Kerteszia) (Diptera: Culicidae) Using Mitochondrial Genes. Insects 2020, 11, 324. https://doi.org/10.3390/insects11050324

AMA Style

Kirchgatter K, de Oliveira Guimarães L, Hugo Yañez Trujillano H, Rafael Arias F, Cáceres AG, de Castro Duarte AMR, dos Santos Malafronte R, Tubaki RM, Mureb Sallum MA. Phylogeny of Anopheles (Kerteszia) (Diptera: Culicidae) Using Mitochondrial Genes. Insects. 2020; 11(5):324. https://doi.org/10.3390/insects11050324

Chicago/Turabian Style

Kirchgatter, Karin, Lilian de Oliveira Guimarães, Henrry Hugo Yañez Trujillano, Fernando Rafael Arias, Abraham Germán Cáceres, Ana Maria Ribeiro de Castro Duarte, Rosely dos Santos Malafronte, Rosa Maria Tubaki, and Maria Anice Mureb Sallum. 2020. "Phylogeny of Anopheles (Kerteszia) (Diptera: Culicidae) Using Mitochondrial Genes" Insects 11, no. 5: 324. https://doi.org/10.3390/insects11050324

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop