Assessing the molecular divergence between Anopheles (Kerteszia) cruzii populations from Brazil using the timeless gene: further evidence of a species complex
© Rona et al; licensee BioMed Central Ltd. 2009
Received: 24 November 2008
Accepted: 09 April 2009
Published: 09 April 2009
Anopheles (Kerteszia) cruzii was the most important vector of human malaria in southern Brazil between 1930–1960. Nowadays it is still considered an important Plasmodium spp. vector in southern and south-eastern Brazil, incriminated for oligosymptomatic malaria. Previous studies based on the analysis of X chromosome banding patterns and inversion frequencies in An. cruzii populations from these areas have suggested the occurrence of three sibling species. In contrast, two genetically distinct groups among An. cruzii populations from south/south-east and north-east Brazil have been revealed by isoenzyme analysis. Therefore, An. cruzii remains unclear.
In this study, a partial sequence of the timeless gene (~400 bp), a locus involved in the control of circadian rhythms, was used as a molecular marker to assess the genetic differentiation between An. cruzii populations from six geographically distinct areas of Brazil.
The timeless gene revealed that An. cruzii from Itaparica Island, Bahia State (north-east Brazil), constitutes a highly differentiated group compared with the other five populations from south and south-east Brazil. In addition, significant genetic differences were also observed among some of the latter populations.
Analysis of the genetic differentiation in the timeless gene among An. cruzii populations from different areas of Brazil indicated that this malaria vector is a complex of at least two cryptic species. The data also suggest that further work might support the occurrence of other siblings within this complex in Brazil.
Anopheles cruzii is one of the few mosquito species belonging to the subgenus Kerteszia. Immature stages of this species are found associated with water trapped in the interfoliar space of plants from the Bromeliaceae family, which are abundant in the Brazilian Atlantic forest [1–3]. Accordingly, the distribution of these bromeliad-breeding mosquitoes is restricted to the Atlantic forest, which stretches from the coast of Rio Grande do Sul State (southern Brazil) to Sergipe State (north-eastern Brazil) [4, 5].
The adults are found in a variety of habitats, from sea level in coastal areas to the mountains. Females are strongly anthropophilic and preferably bite during the evening [2, 6, 7], perhaps biting more than one host to complete egg maturation, which is epidemiologically relevant for malaria transmission [8–10].
Between 1930 and 1960, An. cruzii together with Anopheles bellator and Anopheles homunculus, which also belong to Kerteszia, were considered the main vectors of malaria when the disease was endemic in southern Brazil. Vector control measures have significantly reduced or even interrupted malaria transmission in some areas, but eradication of the pathogen was not achieved and An. cruzii is still responsible for several oligosymptomatic malaria cases in southern and south-eastern Brazil.
The Amazon region is highly endemic for human malaria, caused by Plasmodium vivax and Plasmodium falciparum, and imported cases are frequently reported in different states due to emigration from this region [11, 12]. However, several autochthonous cases were reported in a study in Santa Catarina State, southern Brazil . In the states of São Paulo and Rio de Janeiro, as well as in the state of Bahia, where An. cruzii and Anopheles (Nyssorhynchus) spp. are considered the main vectors of the disease, respectively [3, 7, 13, 14], several imported and autochthonous cases of malaria are reported every year in the Atlantic forest region . Reinforcing the epidemiological importance of An. cruzii as a malaria vector in south-east Brazil, another recent study in Espírito Santo State, including the locality of Santa Teresa, suggested that this species is the potential vector of recent autochthonous cases of malaria in this state .
Anopheles cruzii is also a natural vector of simian malaria in Rio de Janeiro and São Paulo States . Studies on seasonal and vertical distribution of An. cruzii in coastal São Paulo State demonstrated high vertical mobility from ground level to tree tops, with significantly more activity in the uppermost branch layer of the forest . This behaviour could be responsible for human infection by simian Plasmodium species [19, 20].
Epidemiological surveillance and the use of control measures are required to avoid the expansion or introduction of malaria in areas where vector species are abundant and susceptible humans are present. Thus, assessment of the epidemiological status of such localities as well as knowledge concerning the biology, behaviour and the genetic characteristics of the vector species are relevant to prevent the occurrence of outbreaks or to lead control strategies, especially in formerly endemic areas.
Despite its epidemiological importance, there are only a few population genetic studies of An. cruzii [18, 21], and its taxonomic status is unclear. Anopheles cruzii is polymorphic for chromosome rearrangements. Differences in inversions frequencies, and X chromosome banding patterns from south-eastern and southern Brazil, have suggested the existence of three sibling species [21–24]. On the other hand, isoenzymes indicated two genetically isolated groups, one from Bahia State (north-eastern Brazil), and the other from south-eastern and southern Brazil (Rio de Janeiro, São Paulo and Santa Catarina States) . Finally, in a recent study based on sequence analysis of the second Internal Transcribed Spacer of the nuclear ribosomal DNA (ITS2), the authors found no conclusive evidence for sibling species among samples of An. cruzii from south-eastern and southern Brazilian localities .
The activity and feeding rhythms of insect vectors are very important to disease transmission. These patterns are controlled by endogenous circadian clocks, which are under genetic control . Furthermore, clock genes are also involved in the control of mating rhythms that are potentially important in maintaining sexual isolation between closely related species [28, 29].
The circadian rhythms of malaria vectors belonging to the subgenus Kerteszia were formerly studied by Pittendrigh  and, recently, these rhythms were also studied in An. cruzii . The timeless gene is involved in the control of activity rhythms in Drosophila , and controls differences in mating rhythms between closely related Drosophila species .
In the present study, a fragment of ~400 bp of the An. cruzii timeless gene was used as a molecular marker to assess intraspecific variability and genetic divergence among six populations of An. cruzii captured in different locations within the geographic distribution range of this species in Brazil.
Isolation of the An. cruzii timeless gene sequence
Sequence of primers used to amplify the timeless gene fragments
Sequence of primers at 5' → 3'
Interpopulational analysis of the An. cruzii timeless gene
Females were processed individually and genomic DNA was extracted as above . PCR amplification was carried out for 35 cycles at 94°C for 30 s, 62°C for 60 s and 72°C for 90 s using the proofreading Pfu DNA polymerase (Biotools) and primers 5'acbatim02a or 5'cruziitim02 and 3'cruziitim03 (Table 1). Negative controls (no DNA added) were included in all amplification reactions and pre- and post-PCR procedures did not share equipment or reagents. After cloning the fragments obtained as above, at least eight clones of each mosquito were sequenced and two consensus sequences representing both alleles were generated. When only one haplotype was observed among the eight sequences the mosquito was considered a homozygote. The probability that a heterozygote will be mistakenly classified as a homozygote with this procedure is less than 1%. Five mosquitoes were classified as homozygotes in Itatiaia, none in Florianópolis and one in each of the other four populations. The sequences obtained in homozygote mosquitoes were duplicated prior to analysis. However, the population genetics analysis was also carried out without duplicating the homozygote sequences and the results were very similar.
DNA sequence analysis
The timeless gene fragments were aligned using the GCG package (Wisconsin Package Version 10.2, Genetics Computer Group) and ClustalX software . Analyses of the polymorphism and differentiation between populations were performed using DNASP4.0  and PROSEQ programs . F ST was calculated as described by Hudson et al  and significance was evaluated by 1,000 random permutations. Phylogenetic analysis was carried out using MEGA 4.0  using the default parameters.
Isolation of An. cruzii timeless gene fragment
Different PCR schemes were tested to amplify a fragment of the An. cruzii timeless gene (see Methods). Figure 2 shows an alignment of the predicted amino acid sequence encoded by this fragment obtained from An. cruzii compared to the TIMELESS protein of other insect species (D. melanogaster, Aedes aegypti and An. gambiae). A fairly high degree of inter-specific similarity is observed, but the putative protein encoded by 5' end of this fragment is variable, presenting some amino acid changes among the species compared. Figure 2 also shows the approximate positions of the two introns that occur in this region of the gene, as well as the location of the primers used to amplify the fragment from An. cruzii used for the population genetics analysis described below.
Molecular variation and divergence among An. cruzii populations
The geographic distribution of the six Brazilian populations of An. cruzii used in this study is shown at Figure 1. Initially, using the primers 5'cruziitim02 and 3'cruziitim03 (see Figure 2), a ~420 bp fragment of the timeless gene was amplified in all An. cruzii populations analyzed, with the exception of samples from Bahia State (Itaparica Island), which revealed a ~400 bp amplification product, indicating the existence of length variation among the studied populations. The sample from Bahia, however, displayed lower amplification in some cases using these primers, and so a new internal forward primer named 5'acbatim02a (Table 1) was designed based on the initial sequences obtained. Using this new primer in conjunction with 3'cruziitim03, a ~410 bp fragment of timeless gene was obtained for all An. cruzii populations from south and south-east Brazil and a ~390 bp from Bahia.
A total of 24 sequences from Florianópolis, 24 from Cananéia, 22 from Juquitiba, 24 from Itatiaia, 12 from Santa Teresa and 28 from Itaparica (Bahia State) populations were obtained. The sequences were submitted to GenBank (accession numbers: FJ408732 – FJ408865). A full alignment of all sequences is shown in Additional file 1. Most of the base substitutions were silent or occurred within the two introns, which show a number of indels. A few non-synonymous changes were also observed, causing seven amino acid differences among the sequences.
Polymorphisms of all An. cruzii populations
Genetic differentiation between all An. cruzii populations
1. Florianópolis × Cananéia
2. Juquitiba × Santa Teresa
3. Florianópolis × Juquitiba
4. Cananéia × Juquitiba
5. Florianópolis × Itatiaia
6. Florianópolis × Santa Teresa
7. Itatiaia × Santa Teresa
8. Juquitiba × Itatiaia
9. Cananéia × Santa Teresa
10. Cananéia × Itatiaia
11. Florianópolis × Bahia
12. Juquitiba × Bahia
13. Cananéia × Bahia
14. Santa Teresa × Bahia
15. Itatiaia × Bahia
16. *An. cruzii × Bahia
Non-synonymous changes on the timeless gene fragment
05 (first codon base)
Individuals from all populations analysed
06 (second codon base)
Individuals from all populations analysed
Juq66a; Juq66b; Can06b; Can12b
08 (first codon base)
Individuals from south and south-east populations
All individuals from Bahia population and Flo37a; Can02b
18 (second codon base)
All individuals from south and south-east populations and Bahia19a; Bahia33a; Bahia20b
Gl utamin e
Individuals from Bahia population
11 (first codon base)
All individuals from Florianópolis, Cananéia, Juquitiba, Itatiaia and Santa Teresa populations
All individuals from Bahia population
188 (first codon base)
All individuals from Florianópolis, Cananéia, Juquitiba, Itatiaia and Santa Teresa populations
All individuals from Bahia population
275 (first codon base)
All individuals from Florianópolis, Cananéia, Juquitiba, Itatiaia and Santa Teresa populations
All individuals from Bahia population
The values using only coding regions (shown in parentheses in Table 3) show some differences compared with those obtained with the whole sequence. Yet even using the more conserved coding regions, the values of differentiation between the population from Bahia and all others revealed a high number of fixed differences and only a few shared polymorphisms. Among the southern and south-eastern populations, there were shared polymorphisms and no fixed differences, suggesting they belong to the same or to very closely related species.
Divergence time between An. cruzii populations
The estimate of the time of divergence between An. cruzii populations from Bahia and the others were calculated using the Da value based on the third codon positions. This estimate assumed that substitutions rates observed between An. cruzii from Bahia State and the other populations originally from southern regions of Brazil are similar to the estimated rates in the same fragment of the timeless gene between closely related Drosophila persimilis and Drosophila pseudoobscura, species that diverged around 0.85 millions of years ago (MYA) (FlyBase Accession Numbers FBtr0185090 and FBtr0282161, respectively) . The divergence observed for the timeless gene between these two Drosophila species based on the third codon positions is 0.03030. Based on the Da value (0.05426), the estimated time of divergence between An. cruzii populations from south and south-east Brazil and that from Bahia State, is approximately 1.5 MYA.
Genealogy of the An. cruzii timeless sequences
Zavortink  pointed out morphological differences in the larval stage of populations of An. cruzii from Rio de Janeiro and Santa Catarina States, suggesting that An. cruzii could represent more than a single species. A moderately high F ST value between Florianópolis (Santa Catarina State) and Itatiaia (Rio de Janeiro State) populations was reported here. In addition, comparison of Itatiaia with the other populations (excluding Bahia) revealed even higher F ST values, perhaps suggesting that this population is indeed in a process of differentiation and incipient speciation. Moreover, sequences from Itatiaia showed some clustering in the Neighbour-joining tree (Figure 3). Itatiaia was also the least polymorphic population of south and south-east Brazil and showed the highest number of homozygotes suggesting some inbreeding. It is possible that this reflects a smaller effective size and the relative isolation of this population, since its location in a valley between two mountain chains (Serra do Mar and Serra da Mantiqueira – Figure 1) might reduce gene flow with other populations.
In a recent review, Ayala and Coluzzi  argue that many siblings are outcomes of recent speciation processes associated with paracentric inversions, mostly involving the X chromosome. Ramirez and Dessen [23, 24], studying the X chromosome banding patterns and inversion frequencies of distinct populations of An. cruzii from south and south-east Brazil, showed that there are three X chromosomal forms (A, B and C), suggesting a process of incipient speciation acting on An. cruzii populations. Among the localities analysed in this study, only Juquitiba and Cananéia were also investigated by Ramirez and Dessen [23, 24]. They observed that in Juquitiba the majority of mosquitoes had form A and the remainder had form C, while in Cananéia form B predominated with the remainder having form A [23, 24]. Although the differentiation in the timeless gene between these two populations is not high, the F ST value is significant and does not contradict the results of the chromosomal analysis. The relatively low differentiation in timeless among most populations from south and south-east Brazil might reflect introgression at this locus. It would be interesting to analyse the same populations with an X-linked molecular marker to see whether a higher level of differentiation is found.
Recently, Malafronte et al  compared sequences of ITS2 (Internal Spacer Region 2) from several An. cruzii populations from south and south-east Brazil. Although, they found some differences between sequences from different localities, including Juquitiba and Cananéia, they considered premature to conclude based on their results that there are distinct sibling species in the areas they investigated. Similar results were observed by Calado et al  using PCR-RAPD and PCR-RFLP of the ITS2 region.
Very strong evidence was presented here that confirms the existence of a different species in Bahia State, a finding that supports a previous isoenzyme study . The extremely high F ST values detected between this population and the other five populations studied, as well as the higher number of fixed differences observed, show that Bahia represents a different species. This population also presented lower levels of variability than those from south and south-east Brazil, possibly indicating a smaller population size or past founder effects. However, although the isoenzyme heterozygosity reported for Bahia is lower than Cananéia it is similar to that observed in Florianópolis .
A very rough estimate suggests that the divergence between the Bahia population and the more southern populations of An. cruzii possibly occurred around 1.5 MYA, during the Pleistocene. Climate changes during this period such as an intense precipitation decrease and more arid conditions fragmented the Brazilian Atlantic forest  creating refugia that played an important role in the differentiation among populations of a number of forest species, such as marmosets , tree frogs and many others . Forest fragmentation has also been proposed to explain differentiation among populations of the Atlantic forest mosquito Sabethes albiprivus . Since An. cruzii is also a forest-obligate species, it is possible that the Bahia and southern populations of this species complex suffered fragmentation due a constriction of the forest. Although Tajima's D and Fu & Li's D and F statistics were non-significant, they were negative in most cases and that is consistent with population expansion following the forest recovery after the Pleistocene. Analysis of a number of other molecular markers will allow more precise estimates of the divergence time between the Bahia population and those of south and south-east Brazil. It may also help in determining whether further An. cruzii siblings exist in the latter area.
Finally, although malaria cases are reported annually in Bahia State, the main vector implicated in Plasmodium spp. transmission in this area is An. darlingi and not An. cruzii, the most important vector in the southern states. This suggests that the differentiation observed within the An. cruzii complex might also explain aspects of the vectorial capacity of these mosquitoes, however further studies are needed to confirm or reject this hypothesis.
Analysis of the molecular polymorphism and genetic differentiation of the timeless gene among Brazilian populations of An. cruzii indicates that this malaria vector is a complex of at least two cryptic species, one occurring in the north-east (Bahia State) and another in south and south-east Brazil. In addition, the data also suggest that populations of the latter region might also constitute different incipient species and that further work might support the occurrence of other siblings within this complex in Brazil.
The authors are indebted to Dr Rosely dos Santos Malafronte (Instituto de Medicina Tropical de São Paulo), Dr Monique de Albuquerque Motta (FIOCRUZ – Rio de Janeiro) and Claudiney Biral dos Santos (Unidade de Medicina Tropical – Universidade Federal do Espirito Santo) for providing mosquitoes from Cananéia, Juquitiba, Itatiaia and Santa Teresa; to Paulo Amoretty and Robson Costa da Silva for their technical assistance, to Dr André Nóbrega Pitaluga for preparing Figure 1, to Dr Julian Gray for reading the manuscript, and to PDTIS-FIOCRUZ for use of its DNA sequencing facility. This work was supported by grants from the Howard Hughes Medical Institute, FIOCRUZ, Faperj and CNPq.
- Pittendrigh CS: The ectopic specialization of Anopheles homunculus, and its relation to competition with An. bellator. Evolution. 1949, 4: 64-78. 10.2307/2405534.View ArticleGoogle Scholar
- Veloso HP, De Moura JV, Klein RM: Ecological limitation of Anopheles of the Subgenus Kerteszia in the coastal region of Southern Brazil. Mem Inst Oswaldo Cruz. 1956, 54: 517-548.PubMedGoogle Scholar
- Rachou RG: Anofelinos do Brasil: Comportamento dasespécies vetoras de malária. Rev Bras Malariol Doencas Trop. 1958, 10: 145-181.Google Scholar
- Zavortink TJ: A review of the subgenus Kerteszia of Anopheles. Cont Am Entomol Inst. 1973, 9: 1-54.Google Scholar
- Consoli RAGB, Lourenço-de-Oliveira R: Principais mosquitos de importância sanitária no Brasil. 1994, Rio de Janeiro: Ed. FiocruzGoogle Scholar
- Corrêa RR, Forattini OP, Guarita OF, Rabello EX: Observations on the flight of Anopheles (Kerteszia) cruzii and of A. (K.) bellator, vectors of malaria (Diptera, Culicidae). Arq Hig Saude Publica. 1961, 26: 333-342.PubMedGoogle Scholar
- Aragão MB: Geographic distribution and abundance of Anopheles species (Kerteszia) (Diptera, Culicidae). Rev Bras Malariol Doencas Trop. 1964, 16: 73-109.PubMedGoogle Scholar
- Wilkerson RC, Peyton EL: The Brazilian malaria vector Anopheles (Kerteszia) cruzii: Life stages and biology (Diptera: Culicidae). Mosq Syst. 1991, 23: 110-122.Google Scholar
- Forattini OP, Kakitani I, Massad E, Gomes Ade C: Studies on mosquitoes (Diptera: Culicidae) and anthropic environment. 1 – Parity of blood seeking Anopheles (Kerteszia) in south-eastern Brazil. Rev Saude Publica. 1993, 27: 1-8.View ArticlePubMedGoogle Scholar
- Bona AC, Navarro-Silva MA: Anopheles cruzii parity in dense rain forest in Southern Brazil. Rev Saude Publica. 2006, 40: 1118-1123.View ArticlePubMedGoogle Scholar
- MS – Ministério da Saúde, Brazil (Brazilian Ministry of Health). 2006, [http://portal.saude.gov.br/portal/arquivos/pdf/folder_malaria_2006_web.pdf]
- Machado RL, D' Almeida Couto AA, Cavasini CE, Calvosa VS: Malaria outside the Brazilian Amazonian region: the situation in Santa Catarina State. Rev Soc Bras Med Trop. 2003, 36: 581-586. 10.1590/S0037-86822003000500007.View ArticlePubMedGoogle Scholar
- SESAB – Secretaria da Saúde do Estado da Bahia (Bahia State Health Department). [http://www.saude.ba.gov.br]
- Davis NC, Kumm HHW: Further incrimination of Anopheles darlingi Root as a transmitter of malaria. Am J Trop Med. 1932, 12: 93-95.Google Scholar
- SINAN – Sistema de Informação de Agravos de Notificação, Ministério da Saúde. [http://dtr2004.saude.gov.br/sinanweb/novo/]
- Rezende HR, Cerutti C, Santos CB: Aspectos atuais da distribuição geográfica de Anopheles (Kerteszia) cruzii Dyar & Knab, 1908 no Estado do Espírito Santo, Brasil. Entomol Vect. 2005, 12: 123-126.View ArticleGoogle Scholar
- Deane LM, Ferreira-Neto JA, Deane SP, Silveira IP: Anopheles (Kerteszia) cruzii, a natural vector of the monkey malaria parasites, Plasmodium simium and Plamodium brasilianum. Trans R Soc Trop Med Hyg. 1970, 64: 647-10.1016/0035-9203(70)90088-X.View ArticlePubMedGoogle Scholar
- Marrelli MT, Malafronte RS, Sallum MA, Natal D: Kerteszia subgenus of Anopheles associated with the Brazilian Atlantic rainforest: current knowledge and future challenges. Malar J. 2007, 6: 127-134. 10.1186/1475-2875-6-127.PubMed CentralView ArticlePubMedGoogle Scholar
- Deane LM, Ferreira-Neto JA, Lima MM: The vertical dispersion of Anopheles (Kerteszia) cruzii in a forest in southern Brazil suggests that human cases of simian origin be expect. Mem Inst Oswaldo Cruz. 1984, 79: 461-463.View ArticlePubMedGoogle Scholar
- Ueno HM, Forattini OP, Kakitani I: Vertical and seasonal distribution of Anopheles (Kerteszia) in Ilha Comprida, Southeastern Brazil. Rev Saude Publica. 2007, 41: 269-275.View ArticlePubMedGoogle Scholar
- Ramirez CC, Dessen EM: Cytogenetics analysis of a natural population of Anopheles cruzii. Rev Bras Genet. 1994, 17: 41-46.Google Scholar
- Ramirez CC, Dessen EM, Otto PA: Inversion polymorphism in a natural population of Anopheles cruzii. Caryologia. 1994, 47: 121-130.View ArticleGoogle Scholar
- Ramirez CC, Dessen EM: Chromosomal evidence for sibling species of the malaria vector Anopheles cruzii. Genome. 2000, 43: 143-151. 10.1139/gen-43-1-143.View ArticlePubMedGoogle Scholar
- Ramirez CC, Dessen EM: Chromosome differentiated populations of Anopheles cruzii: evidence for a third sibling species. Genetica. 2000, 108: 73-80. 10.1023/A:1004020904877.View ArticlePubMedGoogle Scholar
- Carvalho-Pinto CJ, Lourenço-de-Oliveira R: Isoenzymatic analysis of four Anopheles (Kerteszia) cruzii (Díptera: Culicidae) populations of Brazil. Mem Inst Oswaldo Cruz. 2004, 99: 471-475. 10.1590/S0074-02762004000500002.View ArticlePubMedGoogle Scholar
- Malafronte Rdos S, Marrelli MT, Ramirez CC, 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 ArticlePubMedGoogle Scholar
- Hardin PE: The Circadian Timekeeping System of Drosophila. Curr Biol. 2005, 15: 714-722. 10.1016/j.cub.2005.08.019.View ArticleGoogle Scholar
- Sakai T, Ishida N: Circadian rhythms of female mating activity governed by clock genes in Drosophila. Proc Natl Acad Sci. 2001, 98: 9221-9225. 10.1073/pnas.151443298.PubMed CentralView ArticlePubMedGoogle Scholar
- Tauber E, Roe H, Costa R, Hennessy JM, Kyriacou CP: Temporal mating isolation driven by a behavioral gene in Drosophila. Curr Biol. 2003, 13: 140-145. 10.1016/S0960-9822(03)00004-6.View ArticlePubMedGoogle Scholar
- Pittendrigh CS: The quantitative evaluation of Kerteszia breeding grounds. Am J Trop Med Hyg. 1950, 30: 457-468.PubMedGoogle Scholar
- Chahad-Ehlers S, Lozovei AL, Marques MD: Reproductive and post-embryonic daily rhythm patterns of the malaria vector Anopheles (Kerteszia) cruzii: aspects of the life cycle. Chronobiol Int. 2007, 24: 289-304. 10.1080/07420520701282174.View ArticlePubMedGoogle Scholar
- Jowett T: Preparation of nucleic acids. Drosophila, A Practical Approach. 1998, IRL press, Oxford: Roberts DB, 347-371.Google Scholar
- Gentile C, Meireles-Filho AC, Britto C, Lima JB, Valle D, Peixoto AA: Cloning and daily expression of the timeless gene in Aedes aegypti (Diptera:Culicidae). Insect Biochem Mol Biol. 2006, 36: 878-884. 10.1016/j.ibmb.2006.08.008.View ArticlePubMedGoogle Scholar
- GenBank database. [http://www.ncbi.nlm.nih.gov/BLAST/]
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralView ArticlePubMedGoogle Scholar
- Rozas J, Sánchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.View ArticlePubMedGoogle Scholar
- Filatov DA, Charlesworth D: DNA polimorphism, haplotype structure and balancing selection in the Leavenworthia PgiC locus. Genetics. 1999, 153: 1423-1434.PubMed CentralPubMedGoogle Scholar
- Hudson RR, Slatkin M, Maddison WP: Estimation of levels of gene flow from DNA sequence data. Genetics. 1992, 132: 583-589.PubMed CentralPubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.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
- Tamura K, Subramanian S, Kumar S: Temporal Patterns of Fruit Fly (Drosophila) Evolution Revealed by Mutation Clocks. Mol Biol Evol. 2004, 21 (1): 36-44. 10.1093/molbev/msg236.View ArticlePubMedGoogle Scholar
- Ayala FJ, Coluzzi M: Chromosome speciation: Humans, Drosophila, and mosquitoes. Proc Natl Acad Sci. 2005, 102: 6535-6542. 10.1073/pnas.0501847102.PubMed CentralView ArticlePubMedGoogle Scholar
- Calado DC, Navarro-Silva MA, Sallum MAM: PCR-RAPD and PCR-RFLP polymorphism detected in Anopheles cruzii (Diptera, Culicidae). Rev Bras Entomol. 2006, 50: 423-430. 10.1590/S0085-56262006000300014.View ArticleGoogle Scholar
- Vasconcelos PM, Becker TA, Renne PR, Brimhall GH: Age and duration of weathering by 40K-40Ar and 40Ar/39Ar analysis of potassium-manganese oxides. Science. 1992, 258: 451-455. 10.1126/science.258.5081.451.View ArticlePubMedGoogle Scholar
- Marroig G, Cropp S, Cheverud JM: Systematics and evolution of the Jacchus group of marmosets (Platyrrhini). Am J Phys Anthropol. 2004, 123: 11-22. 10.1002/ajpa.10146.View ArticlePubMedGoogle Scholar
- Carnaval AC, Hickerson MJ, Haddad CF, Rodrigues MT, Moritz C: Stability predicts genetic diversity in the Brazilian Atlantic forest hotspot. Science. 2009, 323: 785-789. 10.1126/science.1166955.View ArticlePubMedGoogle Scholar
- Pedro PM, Sallum MA, Butlin RK: Forest-obligate Sabethes mosquitoes suggest palaeoecological perturbations. Heredity. 2008, 101: 186-95. 10.1038/hdy.2008.45.View ArticlePubMedGoogle Scholar
- Nei M, Kumar S: Molecular Evolution and Phylogenetics. 2000, New York: Oxford University PressGoogle 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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.