- Open Access
Altitudinal population structure and microevolution of the malaria vector Anopheles cruzii (Diptera: Culicidae)
© Lorenz et al.; licensee BioMed Central. 2014
- Received: 26 June 2014
- Accepted: 29 November 2014
- Published: 16 December 2014
In Brazil, the autochthonous transmission of extra-Amazonian malaria occurs mainly in areas of the southeastern coastal Atlantic Forest, where Anopheles cruzii is the primary vector. In these locations, the population density of the mosquito varies with altitude (5–263 m above sea level), prompting us to hypothesise that gene flow is also unevenly distributed. Describing the micro-geographical and temporal biological variability of this species may be a key to understanding the dispersion of malaria in the region. We explored the homogeneity of the An. cruzii population across its altitudinal range of distribution using wing shape and mtDNA gene analysis. We also assessed the stability of wing geometry over time.
Larvae were sampled from lowland (5–20 m) and hilltop (81–263 m) areas in a primary Atlantic Forest region, in the municipality of Cananéia (State of São Paulo, Brazil). The right wings of males and females were analysed by standard geometric morphometrics. Eighteen landmarks were digitised for each individual and a discriminant analysis was used to compare samples from the hilltop and lowland. A 400-bp DNA fragment of the mitochondrial cytochrome oxidase gene subunit I (CO-I) was PCR-amplified and sequenced.
Wing shapes were distinct between lowland and hilltop population samples. Results of cross-validated tests based on Mahalanobis distances showed that the individuals from both micro-environments were correctly reclassified in a range of 54–96%. The wings of hilltop individuals were larger. The CO-I gene was highly polymorphic (haplotypic diversity = 0.98) and altitudinally structured (Фst = 0.085 and Jaccard = 0.033). We found 60 different haplotypes but only two were shared by the lowland and hilltop populations. Wing shape changed over the brief study period (2009–2013).
Wing geometry and CO-I gene analysis indicated that An. cruzii is vertically structured. Wing shape varied rapidly, but altitude structure was maintained. Future investigations should identify the biotic/abiotic causes of these patterns and their implications in the local epidemiology of malaria.
- Geometric morphometrics
- Mitochondrial gene
- Atlantic Forest
Anopheles (Kerteszia) cruzii is a neotropical mosquito that employs the bromeliad phytotelmata as a larval habitat. Females blood-feed on humans as well as on other mammals and birds . This mosquito occurs in natural forest ecosystems , and forest fragments . In spite of its occurrence in forest ecosystems impacted by human activities, An. cruzii is mainly a sylvatic species with a low synanthropy index . This mosquito can feed on blood during the day and at night; however, its activity peaks in twilight periods -. It is a primary vector of Plasmodium spp. parasites in areas within the Atlantic Forest biome in Brazil, especially in the States of São Paulo, Paraná, and Santa Catarina ,,.
Individuals of An. cruzii were found naturally infected with Plasmodium spp. oocysts in the intestine and sporozoites in the salivary glands -. In rural areas of the Juquitiba and São Vicente municipalities, eastern Sao Paulo state, An. cruzii was infected with Plasmodium vivax (0.149%), Plasmodium vivax strain VK247 (0.086%), and either Plasmodium brasilianum or Plasmodium malariae. The occurrence of these species was positively associated with maintenance of transmission of human Plasmodium in the Atlantic Forest in the states of Sao Paulo and Rio de Janeiro  in cycles that involve humans, Alouatta, and Cebus primates . The infective biting rate was low in areas in which infected An. cruzii were found; however, the abundance of the mosquito and voracious blood-feeding behaviour maintain the endemic circulation of human Plasmodium in some areas within the Atlantic Forest domain. Furthermore, An. cruzii is a vector of Plasmodium simium and Plasmodium brasilianum, which cause simian malaria ,. The coastal region of the Atlantic Forest shelters primate species that can be infected by P. simium and P. brasilianum; these parasites can be occasionally transmitted to humans by infective Kerteszia bites ,.
In an ecological study conducted in the municipality of Cananéia as biological markers and sampled individuals of An. cruzii from in Sao Paulo state, Marques et al. found that An. cruzii is heterogeneously distributed across different altitudes in the Atlantic Forest. The authors assessed the distribution of An. cruzii in three microenvironments, which were grouped based on altitude (lowland: 5–20 m, hill slope: 33–54 m, and hilltop: 81–263 m altitude). Follow-up was monthly for one year; results showed that this species occurred more densely in the hilltop, presumably due to ecological preferences. We then hypothesised that gene flow is not homogeneously distributed across the altitude range and sought to determine the transience of this scenario. The aim of this study was to determine 1) if individuals of Anopheles cruzii from lowland and hilltop are morphogenetically similar and 2) if the biological variability of An. cruzii is stable over time. To address these questions, we used wing geometry and CO-I mitochondrial DNA sequences as biological markers and sampled individuals of An. cruzii from Cananéia municipality from 2009–2013.
Immature An. cruzii individuals were taken from water accumulated inside the tanks of terrestrial and epiphytic bromeliads located up to three meters from the ground. Larvae, pupae, and eggs were collected from the water with a manual pump  and kept in the laboratory until emergence of the adults. Species identifications were based on characteristics of the fourth-instar larvae and pupae using the identification key proposed by Forattini et al.  and characteristics suggested by Sallum et al. .
Collection sites were separated into two microenvironments, lowland (5–20 m altitude) and hilltop (81–263 m altitude) landscape categories, both inside a well preserved forest ecosystem. The lowland areas were characterized by high humidity, high tree density, and little sunlight at ground level. In contrast, hilltop areas were abundant in rocky outcrops, tree density and humidity were low, and there was more sunlight at ground level. Co-ordinates of collection sites are described in Additional file 1. A ground distance of 470 m separates the closest sites between lowland and hilltop; the maximum distance was 1880 m.
Field collections were carried out in January 2012, July 2012, and January 2013. Representatives from the January 2009 population were those reported by Marques et al. . The water from each bromeliad tank was kept separate in 500-mL plastic containers in the laboratory under controlled temperature (25 ± 1°C) and humidity (80 ± 10%). Food and larval density in the container were similar for all collections. After emerging, the adults were euthanized and stored in 95% ethanol at -80°C.
Number of individuals used for geometric morphometrics and date of capture of Anopheles cruzii
Extraction, amplification, and purification of mtDNA
Mosquitoes were homogenized as described by Jowett  for extraction of genomic DNA. A total of 96 specimens of An. cruzii were used for genetic analysis. We used all the individuals collected in January 2012 (Table 1); the same sample was used for geometric morphometrics analyses. The CO-I gene fragment was PCR amplified in a reaction mixture of 1× buffer (20 mM Tris-HCl, pH = 8.4), 0.4 mM dNTP, 2.5 mM MgCl2, 0.5 mM each primer, 0.5 μl Taq polymerase, 2.5 μl genomic DNA, and sterile water to a final volume of 20 μl. A 407-bp fragment of the CO-I gene was amplified with forward primer UEA-7 and reverse primer UEA-10 , and did not include the barcode region of mitochondrial DNA. Cycling conditions were as follows: 95°C for 3 min followed by 40 cycles of 94°C for 40 s, 52°C for 40 s, and 72°C for 1 min, and a final extension of 72°C for 10 min. The GenBank accession numbers for the sequences obtained are: KC992738–KC992770. The PCR products were cleaned with a PureLink™ PCR Purification kit (Invitrogen Corporation, Melle, Germany) according to manufacturer instructions. The fragment was sequenced in the forward and reverse directions with ABI PRISM dGTP BigDye® Terminator v3 (Lincoln Centre Drive, Foster City, CA).
Sequencing mtDNA and data analysis
The purified PCR products and an aliquot of the oligonucleotides specific for the CO-I gene (40 pmol/L) were sent to bio-molecular company Genomic and Molecular Engineering for sequencing on an Applied Biosystems model 3130xl sequencer. These sequences were aligned and edited with the program MEGA 5.0 . The hypothesis of strict neutrality was examined with statistics D  and F  and was tested with the program DnaSP v5 . Analysis of haplotype diversity and the number of polymorphic sites were calculated with MEGA 5.0  using Kimura 2-parameter distance. The Φst value, which estimates the genetic differentiation by molecular variance (AMOVA), was generated by the software ARLEQUIN 188.8.131.52 . Genetic analysis of population differentiation and nucleotide diversity were also calculated with Arlequin software and a haplotype network was constructed with TCS 1.12 . To test the similarity between hilltop and lowland haplotypes, we used Jaccard’s coefficient according to Real .
Altitudinal comparisons of chronological samples
Data of sexual dimorphism of An. cruzii using geometric morphometrics
#10 and #17
#2 and #9
#2 and #9
#10 and #17
#2 and #9
#9 and #10
#9 and #10
#9 and #10
Altitudinal comparisons between hilltop and lowland using geometric morphometrics
Centroid size (mm)
Hilltop X Lowland
#1 and #2
#13 and #15
#1 and #2
#3 and #15
#9 and #18
#1 and #13
#1 and #2
#2 and #13
CO-I mitochondrial DNA
Summary statistics for genetic data of Anopheles cruzii : haplotypes, nucleotide diversity and neutrality’s tests
Hilltop + Lowland
The An. cruzii hilltop and lowland populations showed a pattern of wing divergence between all sampled sets. The minimum ground distance between these micro-regions is approximately 470 m, so the mosquitoes could fly between them . Apparently, there are other ecological reasons that prevent them from freely circulating between altitudes. The preference for specific biotopes and bromeliads within each microenvironment may differ in each population. The micro-environmental context may influence behavior, genetic drift, natural selection and indirectly drive genotype and phenotype. This hypothesis was also discussed in a study case conducted in the same locality: Marques et al. found that the volume of water in lowland bromeliads is significantly higher than in hilltop plants; this could be a determining factor in the preferences of each population. The landscapes of hilltop and lowland also differ in other respects: the levels of sunlight, vegetation cover, relative humidity, and resident plant species may influence the preference of anopheline populations and, thus, their differentiation.
Veloso  observed that the species is indifferent to special microclimates, being found in large quantity in various types of biotopes. If the individuals were under selective pressure, as has already been observed in other species of culicids, there would be a tendency to reduce the intraspecific variability of this characteristic ; this was not observed in our study. However, the separation of groups by wing shape indicates a divergence within the species An. cruzii. Micro-environmental variation might affect wing shape, and this interference can be related to flight performance in Culicidae . In Drosophila, numerous genes control natural variation in wing shape , and the development of shape requires a cascade of genes that act throughout development; this may explain why changes in wing shape are sensitive to a variety of environmental stressors  and that minor perturbations during development can lead to large changes in shape .
Although wing shape varied over time, altitudinal structure and sexual wing dimorphism were maintained. Despite its dynamic and rich variability, wings are evolutionarily informative and appear to be canalized. The morphological variation differential for each sex can be the result of different selective pressures that can shape the wings as part of their adaptations. Models by Lande  and Lande and Arnold  assume that the genetic basis for sexual dimorphism is polygenic. Males and females have different ecological roles in the environment and use their wings differently; this could motivate sex-specific natural selection . The comparison of wing shape consensus after Procrustes superimposition revealed that the most influential landmark in the majority of populations was landmark #9, localized in the wing border. It seems that this anatomical mark is important in sexual dimorphism of this species. According to Dujardin , the landmarks are differentially affected by the same displacement. The posterior border of the wing moved in opposite ways, depending on the external stimulus . In terms of allometry, the contribution of size to the shape variation was statistically significant for most populations; however, the elimination of allometry before CVA allowed us to say that the observed variation probably is not due to plasticity.
In this study, we found a very high number of haplotypes for each microenvironment sampled, only two of which were shared between lowland and hilltop insects. The nucleotide diversity was also high, indicating that the analysed fragment of CO-I is extremely polymorphic in An. cruzii. The large number of haplotypes found differs from other studies of the anopheline CO-I gene ,. The high frequency of transitions in the third codon position indicates that the analysed groups diverged recently . There have been no published population-based studies of mtDNA in An. cruzii. The large quantity of haplotypes found in micro-environmental comparisons may reflect the proximity to the centre of origin of this species. The distribution of An. cruzii is restricted to South America ; thus, its centre of origin is also located in this region. In general, older populations have a higher diversity than younger populations ,.
Comparing the values of genetic diversity (ΦST) with those of morphological diversity (Qst), we noted that Qst of shape and size, were higher. This reveals a greater degree of evolution in the wing than in the CO-I gene; therefore, these markers evolve at different rates. If the markers were neutral they should converge to the same value , which is not the case here. In An. cruzii, the CO-I gene was so polymorphic that might not be the most appropriate marker to demonstrate a clear pattern of divergence between populations, as in Anopheles darlingi, or Aedes aegypti. Although Tajima  and Fu and Li’s  F neutrality tests did not observe significant negative population values (Table 4), the results of mismatch distribution analyses indicated that the population probably did not have a constant size, but remained in growth-decline (Figure 7).
Studies of Culex and Aedes have demonstrated that altitudinal stratification can be revealed by individual wing morphology; however those studies have compared distant populations, with altitudes ranging from 800 to 2130 m. Our study revealed variation between closely located populations, where the geographical distance between microenvironments was not substantial (see Additional file 1). An incidental hypothesis of our work is that the lowland/hilltop populations are undergoing an incipient speciation process. The presence of cryptic species under the nominal species An. cruzii has been reported in other localities, but the occurrence of this phenomenon in Cananéia is still controversial ,. To test this hypothesis, new studies should be performed using the other taxonomical markers traditionally employed to diagnose cryptic species in An. cruzii: chromosomes ,,.
The data for wing shape and gene CO-I were concordant, which probably indicates the vertical structure of An. cruzii in the Cananéia region. Despite microevolution over time, the altitudinal structure and sexual wing dimorphism were maintained in this species. The hilltop and lowland populations differ in several aspects, which might be also reflected in their vectorial capacity. Our proposed scenario relies on theoretical conditions conducive to development of parapatric speciation, in which there is no specific extrinsic barrier to gene flow. The population is continuous, but the individuals do not cross randomly . They are more prone to cross with their geographical neighbours than with individuals of a different group. There is continuity between the landscape microenvironments, because the bromeliaceous breeding sites are present throughout the forest and there are no physical barriers between lowland and hilltop. However, the morphological and molecular evidence suggests that the anophelines do not pass randomly through the vertical stratum. Knowledge of these phenomena will support our understanding of the transmission dynamics of the Plasmodium sp. It is also possible that the gene flow patterns differ throughout the year or that population makes demographic substitution. The gene variability in populations is unknown; however, if this scenario changes, it can be used to understand the epidemiology of malaria in the region. It opens a demand for studies on sexual dimorphism and genetic basis of altitudinal preference.
A number of people provided field and laboratory assistance or mosquito samples used in the study. We are especially grateful to Luiz Carlos de Oliveira, Décio Wach, Jair Donizete da Silva, Francisco Alves dos Santos, Rui de Lima and Solange Schiessl Lorenz. Financial support: CNPq no. 140964/2013-4; CAPES Research Grant #23038.005274/2011-24; FAPESP no. 2011/30397-7; FAPESP no. 2013/05521-9.
- Consoli RAGB, Lourenço-de-Oliveira R: Principais mosquitos de importância sanitária no Brasil. 1994, Fiocruz Press, Rio de JaneiroGoogle Scholar
- Branquinho MS, Marrelli MT, Curado I, Natal D, Barata JMS, Tubaki R, Carréri-Bruno GC, Menezes RT, Kloetzel JK: Infecção do Anopheles (Kerteszia) cruzii por Plasmodium vivax e Plasmodium vivax variante VK247 nos Municípios de São Vicente e Juquitiba, São Paulo. Rev Panam Salud Publica. 1997, 2: 122-134. 10.1590/S1020-49891997000900004.View ArticleGoogle Scholar
- Forattini OP: Culicidologia médica. 2002, EDUSP Press, São PauloGoogle Scholar
- Forattini OP, Gomes AC, Natal D, Santos JLF: Observações sobre a atividade de mosquitos Culicidae em matas primitivas da planície e perfis epidemiológicos de vários ambientes do Vale do Ribeira, São Paulo, Brasil. Rev Saude Publ. 1986, 20: 178-203.Google Scholar
- Forattini OP, Kakitani I, Massad E, Marucci D: Studies on mosquitoes (Diptera: Culicidae) and anthropic environment: biting activity and blood-seeking parity of Anopheles (Kerteszia) in South-Eastern Brazil. Rev Saude Publ. 1996, 30: 107-114. 10.1590/S0034-89101996000200001.View ArticleGoogle Scholar
- Forattini OP, Gomes ADC, Natal D, Kakitani I, Marucci D: Freqüência domiciliar e endofilia de mosquitos Culicidae no Vale do Ribeira, São Paulo, Brasil. Rev Saude Publ. 1987, 21: 188-192. 10.1590/S0034-89101987000300003.View ArticleGoogle Scholar
- Gadelha P: From “forest malaria” to “bromeliad malaria”: a case-study of scientific controversy and malaria control. Parassitologia. 1994, 36: 175-195.PubMedGoogle Scholar
- Tubaki RM, Carreri-Bruno GC, Glasser CM, Ciaravolo RMC: Biting activity of Anopheles (Kerteszia) cruzii (Diptera, Culicidae) in domiciliary habitats in the Southern Atlantic Forest, Peruibe, State of São Paulo, Brazil. Rev Bras Entomol. 1993, 37: 569-575.Google Scholar
- Azevedo AL: Aspectos da epidemiologia da malária e da biologia de Anopheles (Kerteszia) cruzii Dyar & Knab em vales montanhosos do sistema de Mata Atlântica. 1997, Instituto Oswaldo Cruz Press, Rio de JaneiroGoogle Scholar
- Carvalho MED, Glasser CM, Ciaravolo RM, Etzel A, Dos Santos LA, Ferreira CS: Sorologia de malária vivax no foco Aldeia dos Índios, município de Peruíbe, Estado de São Paulo, 1984 a 1986. Cadernos de Saúde Pública. 1988, 4: 276-292. 10.1590/S0102-311X1988000300003.View ArticleGoogle Scholar
- Curado I, Duarte AMR, Lal AA, Oliveira SG, Kloetzel JK: Antibodies anti-bloodstream and circumsporozoite antigens (Plasmodium vivax and Plasmodium malariae/P. brasilianum) in areas of very low malaria endemicity in Brazil. Mem Inst Oswaldo Cruz. 1997, 92: 235-243. 10.1590/S0074-02761997000200017.View ArticlePubMedGoogle Scholar
- Amaral J: Infecção natural do Nyssorhynchus (Kerteszia) cruzii e bellator (Diptera Culicidae). Folha med. 1942, 15: 1-9.Google Scholar
- Coutinho JO, Rachou RG, Ricciardi I: Contribuição para o conhecimento dos transmissores de malária no Brasil - Anopheles (K) cruzii e bellator. Arq Hig. 1943, 13: 115-127.Google Scholar
- Duarte AMR, Pereira DM, de Paula MB, Fernandes A, Urbinatti PR, Ribeiro AF, Malafronte RS: Natural infection in anopheline species and its implications for autochthonous malaria in the Atlantic forest in Brazil. Parasit Vectors. 2013, 6: 58-10.1186/1756-3305-6-58.PubMed CentralView ArticlePubMedGoogle Scholar
- Rachou RG, Ferreira MO:Kerteszia as the group most responsible for malaria in the State of Santa Catarina. Rev Bras Malar Doen Trop. 1966, 18: 691-Google Scholar
- Deane LM, Ferreira Neto JA, Deane MP, Silveira IP:Anopheles (Kerteszia) cruzii, a natural vector of the monkey malaria parasites, Plasmodium simium and Plasmodium brasilianum. Trans Roy Soc Trop Med Hyg. 1970, 4: 64-Google Scholar
- Deane LM: Simian malaria in Brazil. Mem Inst Oswaldo Cruz. 1992, 87: 1-20. 10.1590/S0074-02761992000700001.View ArticlePubMedGoogle Scholar
- Deane LM, Ferreira-Neto JA, MOURA LM: 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. 10.1590/S0074-02761984000400011.View ArticlePubMedGoogle Scholar
- Marques TC, Bourke BP, Laporta GZ, Sallum MAM: Mosquito (Diptera: Culicidae) assemblages associated with Nidularium and Vriesia bromeliads in Serra do mar, Atlantic Forest. Brazil Parasit Vectors. 2012, 5: 41-10.1186/1756-3305-5-41.View ArticlePubMedGoogle Scholar
- Prefeitura Municipal da Estância de Cananéia. São Paulo. Brazil [http://www.cananeia.sp.gov.br/index.php?exibir=secoes&ID=50]
- Centro de Pesquisas Meteorológicas e Climáticas aplicadas à Agricultura (CEPAGRI). Meteorologia Universidade Estadual de Campinas – UNICAMP; [http://www.cepagri.unicamp.br/]
- IBGE (Instituto Brasileiro de Geografia e Estatística). [http://www.ibge.gov.br/home/]
- Belkin JN, Schick RX, Heinemann SJ: Mosquito studies (Diptera: Culicidae). VI. Mosquitoes originally described from North America. Contrib Am Entomol Inst. 1966, 1: 1-39.Google Scholar
- Sallum MAM, Dos Santos CL, Wilkerson RC: Studies on Anopheles (Kerteszia) homunculus Komp (Diptera: Culicidae). Walter Reed Army Inst of Res Silver Spring MD Dept Entomol. 2009, 3: 40-Google Scholar
- Lorenz C, Marques TC, Sallum MAM, Suesdek L: Morphometrical diagnosis of the malaria vectors Anopheles cruzii, An. homunculus and An. bellator. Parasit Vectors. 2012, 5: 1-7. 10.1186/1756-3305-5-257.View ArticleGoogle Scholar
- Vidal PO, Rocha LS, Peruzin MCJ: Wing diagnostic characters for mosquitões Culex quinquefasciatus and Culex nigripalpus (Diptera; Culicidae). Rev Bras Entomol. 2011, 55: 132-142. 10.1590/S0085-56262011000100022.View ArticleGoogle Scholar
- Rohlf FJ: TpsDig, digitize landmarks and outlines [Computer program]. Version 2.05. Dept Ecol.& Evol., St. Univ. New York: Stony Brook; 2005.Google Scholar
- Rohlf FJ: TpsRelw, relative warps analysis [Computer program]. Version 1.36. Dept Ecol.& Evol., St. Univ. New York: Stony Brook; 2003.Google Scholar
- StatSoft, Inc: Statistica for Windows [Computer program manual]. Tulsa: OK: StatSoft, Inc; 2004. Available at: [http://www.statsoft.com]
- Dujardin JP: Morphometrics applied to medical entomology. Infect Genet Evol. 2008, 8: 875-890. 10.1016/j.meegid.2008.07.011.View ArticlePubMedGoogle Scholar
- Spitze K: Population structure in Daphnia obtusa: quantitative genetic and allozymic variation. Genetics. 1993, 135: 367-374.PubMed CentralPubMedGoogle Scholar
- Dujardin JP: COO, MOG and COV for Windows [Computer program manual]. Marseille: Institut de Recherches pour le Développement; 2010.Google Scholar
- MOME (Morphometrics in Medical Entomology). [http://mome-clic.com/]
- Klingenberg CP: MorphoJ: an integrated software package for geometric morphometrics. Mol Ecol Res. 2011, 11: 353-357. 10.1111/j.1755-0998.2010.02924.x.View ArticleGoogle Scholar
- TPS (Morphometrics at SUNY Stony Brook). [http://life.bio.sunysb.edu/morph/]
- Jowett T: Preparation of Nucleic Acids. Drosophila, a Practical Approach. 1986, Oxford, U.K, PressGoogle Scholar
- Zhang DX, Hewitt E: Assessment of the universality and utility of a set of conserved mitochondrial COI primers in insects. Insect Mol Biol. 1996, 6: 143-150. 10.1111/j.1365-2583.1997.tb00082.x.View ArticleGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA 5.0: Molecular Evolutionary Genetics Analysis (MEGA) software version 5.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123: 585-595.PubMed CentralPubMedGoogle Scholar
- Fu YX, Li WH: Statistical tests of neutrality of mutations. Genetics. 1993, 133: 693-709.PubMed CentralPubMedGoogle Scholar
- Rozas J, Sánchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphisms analyses by the coalescente and other methods. Bioinformatics Application Note. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.View ArticleGoogle Scholar
- Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evol Bioinformatics Online. 2005, 1: 47-50.Google Scholar
- Clement M, Posada DCKA, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-1659. 10.1046/j.1365-294x.2000.01020.x.View ArticlePubMedGoogle Scholar
- Real R: Tables of significant values of Jaccard's index of similarity. Miscellània Zoològica. 1999, 22: 29-40.Google Scholar
- Bernasconi MV, Valsangiacomo C, Piffaretti JC, Ward PI: Phylogenetic relationships among Muscoidea (Diptera: Calyptratae) based on mitochondrial DNA sequences. Insect Mol Biol. 2000, 9: 67-74. 10.1046/j.1365-2583.2000.00158.x.View ArticlePubMedGoogle Scholar
- Mirabello L, Conn JE: Molecular population genetics of the malaria vector Anopheles darlingi in Central and South America. Heredity. 2006, 96: 311-321. 10.1038/sj.hdy.6800805.View ArticlePubMedGoogle Scholar
- Ferreira E, Castro-Filho J, Tomich A, Sá FT: Estudo sobre o raio de vôo do Anopheles (Kerteszia) cruzii e do Anopheles (Kerteszia) bellator em Guaratuba, litoral do Estado do Paraná. Brasil Rev Bras Malariol Doen Trop. 1969, 21: 819-822.Google Scholar
- Veloso HP, Moura JV, Klein RM: Delimitação Ecológica dos anofelinos do subgênero Kerteszia na região costeira do sul do Brasil. Mem Inst Oswaldo Cruz. 1956, 54: 517-520.PubMedGoogle Scholar
- Dujardin JP, Costa J, Bustamante D, Jaramillo N, Catalá S: Deciphering morphology in Triatominae: the evolutionary signals. Acta Trop. 2009, 110: 101-111. 10.1016/j.actatropica.2008.09.026.View ArticlePubMedGoogle Scholar
- Hoffmann AA, Shirriffs J: Geographic variation for wing shape in Drosophila serrata. Evolution. 2002, 56: 1068-1073. 10.1111/j.0014-3820.2002.tb01418.x.View ArticlePubMedGoogle Scholar
- Weber K, Eisman R, Morey L, Patty A, Sparks J, Tausek M, Zeng ZB: An analysis of polygenes affecting wing shape on chromosome 3 in Drosophila melanogaster. Genetics. 1999, 153: 773-786.PubMed CentralPubMedGoogle Scholar
- Garcia-Bellido A, De Celis JF: Developmental genetics of the venation pattern of Drosophila. Annu Rev Genet. 1992, 26: 275-302. 10.1146/annurev.genet.26.1.275.View ArticleGoogle Scholar
- Emlen DJ, Nijhout HF: The development and evolution of exaggerated morphologies in insects. Annu Rev Entomol. 2000, 45: 661-708. 10.1146/annurev.ento.45.1.661.View ArticlePubMedGoogle Scholar
- Lande R: Sexual dimorphism, sexual selection, and adaptation in polygenic characters. Evolution. 1980, 4: 292-305. 10.2307/2407393.View ArticleGoogle Scholar
- Lande R, Arnold SJ: Evolution of mating preference and sexual dimorphism. J Theor Biol. 1985, 117: 651-664. 10.1016/S0022-5193(85)80245-9.View ArticlePubMedGoogle Scholar
- Wekesa JW, Brogdon WG, Hawley WA, Besansky NJ: Flight tone of field collected populations of Anopheles gambie and An. arabiensis (Diptera: Culicidae). Physiol Entomol. 1998, 23: 289-294. 10.1046/j.1365-3032.1998.233087.x.View ArticleGoogle Scholar
- Oshagui MA, Shemshad KH, Yaghobi-Ershadi MR, Pedram M, Vatandoost H, Abaie MR, Akbarzadeh K, Mohtarami F: Genetic structure of the malaria vector Anopheles superpictus in Iran using mitochondrial cytochrome oxidase (COI and COII) and morphologic markers: a new species complex?. Acta Trop. 2006, 101: 241-248. 10.1016/j.actatropica.2007.02.006.View ArticleGoogle Scholar
- Patsoula E, Samanidou-Voyadjoglou A, Spanakos G, Kremastinou J, Nasioulas G, Vakalis C: Molecular characterization of the Anopheles maculipennis complex during surveillance for the, Olympic Games in Athens. Med Vet Entomol. 2004, 2007 (21): 36-43.Google Scholar
- Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P: Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann Entomol Soc Am. 1994, 87: 651-701. 10.1093/aesa/87.6.651.View ArticleGoogle Scholar
- Loaiza JR, Bermingham E, Sanjur OI, Scott ME, Bickersmith SA, Conn JE: Review of genetic diversity in malaria vectors (Culicidae: Anophelinae). Infec Gen Evol. 2012, 12: 1-12. 10.1016/j.meegid.2011.08.004.View ArticleGoogle Scholar
- Molina-Cruz A, De Merida AMP, Mills K, Rodriguez F, Schoua C, Yurrita MM, Black WC: Gene flow among Anopheles albimanus populations in Central America, South America, and the Caribbean assessed by microsatellites and mitochondrial DNA. Am J Trop Med Hyg. 2004, 71: 350-359.PubMedGoogle Scholar
- Kambhampati S, Rai KS: Mitochondrial DNA variation within and among populations of the mosquito Aedes albopictus. Genome. 1991, 34: 288-292. 10.1139/g91-046.View ArticlePubMedGoogle Scholar
- Motoki MT, Suesdek L, Bergo ES, Sallum MAM: Wing geometry of Anopheles darlingi Root (Diptera: Culicidae) in five major Brazilian ecoregions. Infect Gen Evol. 2012, 12: 1246-1252. 10.1016/j.meegid.2012.04.002.View ArticleGoogle Scholar
- Moutinho PR, Gil LH, Cruz RB, Ribolla PE: Population dynamics, structure and behavior of Anopheles darlingi in a rural settlement in the Amazon rainforest of Acre, Brazil. Malaria J. 2011, 10: 174-10.1186/1475-2875-10-174.View ArticleGoogle Scholar
- Scarpassa VM, Cardoza TB, Cardoso RP: Population genetics and phylogeography of Aedes aegypti (Diptera: Culicidae) from Brazil. Am J Trop Med Hyg. 2008, 78: 895-903.PubMedGoogle Scholar
- Demirci B, Lee Y, Lanzaro GC, Alten B: Altitudinal genetic and morphometric variation among populations of Culex theileri Theobald (Diptera: Culicidae) from Northeastern Turkey. J Vector Ecol. 2011, 37: 197-209. 10.1111/j.1948-7134.2012.00217.x.View ArticleGoogle Scholar
- Kuclu O, Aldemir A, Demirci B: Altitudinal variation in the morphometric characteristics of Aedes vexans Meigen from northeastern Turkey. J Vector Ecol. 2011, 36: 30-41. 10.1111/j.1948-7134.2011.00138.x.View ArticlePubMedGoogle Scholar
- Carvalho-Pinto CJD, Lourenço-de-Oliveira R: Isoenzimatic analysis of four Anopheles (Kerteszia) cruzii (Diptera: Culicidae) populations of Brazil. Mem Inst Oswaldo Cruz. 2004, 99: 471-475. 10.1590/S0074-02762004000500002.View ArticlePubMedGoogle Scholar
- Malafronte RS, Marrelli MT, Ramirez CCL, Nassar MN, Marinotti O: Intraspecific variation of second internal transcribed spacer of nuclear ribosomal DNA among populations of Anopheles (Kerteszia) cruzii (Diptera: Culicidae). J Med Entomol. 2007, 44: 538-542. 10.1603/0022-2585(2007)44[538:IVOSIT]2.0.CO;2.View ArticleGoogle Scholar
- Ramirez CCL, Dessen EMB: Cytogenetic analysis of a natural population of Anopheles cruzii. Rev Bras Gen. 1994, 17: 41-46.Google Scholar
- Ramírez CCL, Dessen EMB: Chromosomal evidence for sibling species of the malaria vector Anopheles cruzii. Genome. 2000, 43: 143-151. 10.1139/g99-103.View ArticlePubMedGoogle Scholar
- Templeton AR: Genética de populações e teoria microevolutiva. 2011, Holos Press, São PauloGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.