Hybrid sterility in crosses between two Brazilian sibling species of the Anopheles albitarsis complex
Parasites & Vectorsvolume 7, Article number: 559 (2014)
Complexes of cryptic species are common in several taxa and this is also the case in the Anopheles genus, a group including all known human malaria vectors. The Anopheles albitarsis complex comprises at least nine cryptic species, some of which are implicated as vectors of human malaria. Several different types of data have been generated for this species complex such as cytogenetics, alloenzymes, morphological and feeding behavioral, hybridization experiments, RAPD-PCR and RFLP and mitochondrial and nuclear markers. Studies focused on its postzygotic isolation are still somewhat rare in the literature despite their importance to understand the speciation process and the level of gene flow potentially occurring among the different sibling species.
Hybridization experiments between Anopheles albitarsis s.s. and Anopheles marajoara, as well as backcrosses between hybrids and Anopheles albitarsis s.s., were performed using the induced mating technique. Results were compared to intraspecific crosses. Larva-to-adult viability and sex ratio were also assessed.
Male hybrids show very low insemination rates and nearly complete sterility, apparently due to abnormalities in their reproductive organs. Evidence of partial sterility among the hybrid females was also observed.
Our data indicated that Anopheles albitarsis s.s. and Anopheles marajoara show a high level of postzygotic isolation with a strong hybrid male sterility. This result is consistent with the Haldane’s rule which states that in interspecific crosses the heterogametic sex is the first to be affected. However, the fact that the females are not completely sterile raises the possibility of introgression between these two siblings species.
Complexes of cryptic species are common in several taxa , and this is the case in the Anopheles genus, a group including all known human malaria vectors ,. A better understanding of cryptic species can be epidemiologically relevant since it could offer clues about differential vector capacity (reviewed in ).
Anopheles albitarsis s.l. Lynch-Arribálzaga is widely distributed from northern Guatemala to northern Argentina  and is a complex of cryptic species that includes some of the most important regional vectors in America -. The first indication of An. albitarsis s.l. occurring as a species complex was described by Kreutzer et al. . Since then, several approaches were applied to investigate and characterize the species members: cytogenetics, alloenzymes, morphological and feeding behavioral analysis, hybridization experiments, RAPD-PCR and RFLP and DNA sequencing of mitochondrial and nuclear markers ,-.
Current knowledge suggests that the An. albitarsis complex comprises at least nine species: (i) An. albitarsis s.s. Lynch-Arribálzaga or An. albitarsis A, found in Argentina, Paraguay and Brazil (States of Bahia, Paraná, Santa Catarina, and São Paulo); (ii) Anopheles oryzalimnetes Wilkerson & Motoki or An. albitarsis B, widespread in Paraguay and probably the species with the wider geographical distribution in Brazil (States of Bahia, Ceará, Espírito Santo, Mato Grosso, Pará, Paraná, Rio de Janeiro and São Paulo); (iii) Anopheles marajoara Galvão & Damasceno or An. albitarsis C, found in Paraguay, Venezuela and in Brazil (States of Amapá, Amazonas, Mato Grosso, Pará, Rondônia and São Paulo); (iv) Anopheles deaneorum Rosa-Freitas or An. albitarsis D, distributed in northern Argentina and in Brazil (States of Acre, Mato Grosso, Paraná and Rondônia); (v) Anopheles janconnae Wilkerson & Sallum or An. albitarsis E, found in northern Brazil (States of Amapá, Pará and Roraima) and in Venezuela; (vi) An. albitarsis F, in Colombia, Venezuela and Trinidad; (vii) An. albitarsis G, exclusive to Brazil (States of Amazonas, Bahia and Pará); (viii) An. albitarsis H, also restricted to Brazil (States of Mato Grosso and Rondônia); and (ix) An. albitarsis I, found in Colombia and Venezuela ,,-. Moreover, the number of species within the Anopheles albitarsis complex can be even higher as An. deaneorum and An. marajoara might include more than one sibling species ,.
In an evolutionary context, analysis of the reproductive isolation between cryptic species of insect vectors allows the identification of potential gene flow among siblings. In addition, laboratory crossing experiments are likely to reveal the exchange of genes potentially involved in vectorial capacity or insecticide resistance, being therefore important for vector control programs ,-. In the Albi tarsis species complex only two hybridization studies between siblings have been published so far ,. This is probably due to the difficulty of keeping these species in captivity, since neotropical Anopheles albitarsis species rarely perform free mating crosses under laboratory conditions. In 1991, Klein et al.  observed male hybrid sterility in crosses between An. deaneorum and a species that they called “B”, later defined as An. marajoara (Richard C. Wilkerson, personal communication). In the second study, Lima et al.  observed that hybrid males from crosses between An. albitarsis s.s. and An. deaneorum showed very low insemination rates and suggested that the abnormalities in their reproductive organs could explain this. In the present work, we performed crossing experiments between two Brazilian sibling species of the An. albitarsis complex, An. albitarsis s.s. and An. marajoara. We used these two species in order to analyze the degree of reproductive isolation between them and compare with previous ones. The present work represents an interesting study from an evolutionary point of view of An. albitarsis complex and is epidemiologically relevant focusing a Neotropical malaria vector such as An. marajoara.
We analyzed specimens of An. albitarsis s.s. from a colony maintained at the Instituto de Biologia do Exército (Rio de Janeiro State, Brazil) following the protocol described by Horosko et al. with slight modifications. The colony was established in 1993 with mosquitoes collected from Massaranduba municipality (26° 36′ 39″ S, 49° 0′ 33″ W) (Santa Catarina State, Brazil) . This species was treated as a control in crosses since the specimens used were obtained from a stable colony . In addition, An. marajoara females were collected in Mazagão municipality (0° 6′ 58″ S, 51° 17′ 10″ W) (Amapá State, Brazil). Specimens were identified as Anopheles albitarsis s.l. according to Faran & Linthicum  and were considered as An. marajoara based on Conn et al.. The identification was then confirmed by DNA barcoding (see below). Laboratory rearing of An. marajoara was performed according to Lima et al.. The An. marajoara specimens belonging to F1-F3 generations were used in crossing experiments. All procedures were carried out under controlled conditions: 27°C for the immature stages (larvae to pupae) and 25°C and 70% relative humidity for the adult mosquitoes.
Identification by DNA barcoding
DNA barcodes (658 bp of the mtDNA Cytochrome c Oxidase – COI) were generated for 11 An. albitarsis s.s. and 12 An. marajoara specimens (Accession numbers: KM391793 to KM391815). DNA extraction was carried out as in Jowett  and PCR amplification using the LCO1490 and HCO2198 primers of Folmer et al. . Sequencing reactions were carried out in both directions with the Big Dye Terminator Kit and ABI 3730 automated sequencer (PE Applied Biosystems). Sequences were edited in Bioedit 7.2.3  and the delineation of species within Albitarsis group by phylogenetic analysis was carried out using specimens from each distinctive COI lineage found by NJ-K2P analysis. The sequences added were: An. albitarsis H: GenBank: DQ076222, DQ076223, DQ076224; An. deaneorum: GenBank: DQ076226, DQ076227, DQ076229, DQ076230; An. albitarsis G: GenBank: DQ076221, DQ076225; An. oryzalimnetes: GenBank: DQ076210, DQ076211, DQ076213, DQ0762105. We used the Kimura-two-parameter distance model (K2P), Neighbor-joining analysis (NJ) and 1,000 replicates to produce an unrooted consensus tree.
Methods for crossing experiments were performed as described in Lima et al.. Males and females were separated every twelve hours to guarantee the collection of virgin females . Females were fed according to the protocol approved by the Fiocruz Ethical Committee on Use of Animals (CEUA) and subjected to the artificial-mating technique . Four days after mating, the females were placed in a plastic device (8.5 cm-diameter × 4.5 cm-depth) closed with a nylon mesh and submerged in a bowl (18 cm-diameter × 7.5 cm-depth) filled with dechlorinated water. The internal surface of the plastic device was covered with moist filter paper to keep eggs from drying out. After two days, the spermathecae (from alive and dead females) were dissected and examined by light microscopy; mating success was scored based on the presence of sperm.
Eggs were counted after forced egg laying. Two days after hatching, larvae were counted and transferred to water-filled plastic bowls. Every two days numbers of larvae and living pupae were counted. The immature stages were reared following the same protocol already mentioned for the colony maintenance. Adults were used in the subsequent crosses. Intra and interspecific mating experiments were performed and the hybrid offspring were used in backcrosses according to Lima et al. . The parameters evaluated for each crossing experiment were: the number of inseminated females, the total number of eggs, hatching eggs, larvae, pupae and adults, and also the male and female ratio.
Evaluation of the male genital apparatus
Male genital structures of the offspring were dissected and examined using stereoscopy microscopy. External and internal genital structures were mounted in slides with saline solution, and in some cases, were stained using the buffer reference standard pH = 4.0 (Sigma) for a better visualization. Chi-square statistical analysis was estimated by the program GraphPad Prism 5.0 .
Figure 1 shows a phylogenetic tree with six clusters clearly defined: 1) An. albitarsis H; 2) An. deaneorum; 3) An. albitarsis G; 4) An. oryzalimnetes; 5) An. albitarisis s.s. and 6) An. marajoara. The latter group included all individual sequences of An. marajoara showing that we used a pure lineage on our crossing experiments.
Table 1 shows larvae-to-adult viability of the offspring from intraspecific and interspecific crosses between An. albitarsis s.s. and An. marajoara, and also backcrosses of these hybrids with An. albitarsis s.s.. In each case, eight different crosses were conducted. The hatching number from the different groups of crosses was highly variable. The two intraspecific and the interspecific crosses between male An. marajoara × female An albitarsis s.s. showed very similar hatching percentages (ranged between 38% and 45%). Nevertheless, a remarkable reduction on hatching rates was observed in all backcrosses. Furthermore, crosses between hybrid males and An. albitarsis s.s. females resulted in null or near to zero rates. In general, the pupation was a more homogeneous developmental event than hatching. A significant reduction of the larvae viability was observed in crosses that involved females of An. marajoara, in conspecific crosses of An. marajoara (X2 = 26.34, df = 1, P < 0.001) and male An. albitarsis s.s. × female An. marajoara (X2 = 4.48, df = 1, P = 0.034). Differing from the dramatically decreased hatching in backcrosses, pupation and adult emergence seemed unaffected when hybrid females were used. Intraspecific and interspecific crosses also resulted in similar and high percentages of emergence. Finally, we observed that in all cases the sexual ratio was not statistically different from 1:1 indicating similar mortality between genders.
Evaluation of insemination
After oviposition, the eggs were counted and the spermathecae dissected and classified as positive or negative depending on the presence or lack of sperm, respectively. Females with positive spermathecae were considered inseminated. Table 2 shows the number of mated and positive females for each performed cross. As expected, a higher percentage of inseminated females was observed in crosses involving An. albitarsis s.s. males since they result from induced mating, which introduces a level of selection . The low An. marajoara male effectiveness can be attributed to the use of F1 and F2 mosquito generations not adapted to laboratory conditions. The crosses with An. marajoara males were really challenging, even compared with other Anopheline species, such as An. deaneorum. The number of inseminated females was significantly lower in crosses with hybrid males than with An. albitarsis s.s. In some cases, females mated with interspecific males showed immotile sperm or even agglomerated sperm with dark pigmentation (data not shown).
Figure 2 shows comparisons between different crosses using the normalized insemination rates estimated from Table 2 data. Normalization was carried out using results from the intraspecific cross of An. albitarsis s.s.. Figure 2A illustrates the normalized insemination rates involving females of different genotypes crossed to An. albitarsis s.s. males. In all cases, the difference between female insemination rates was not significant, regardless of their genotype. The results observed for the two types of hybrid females were expected since they have a similar genotype. Figure 2B shows the normalized insemination data of An. albitarsis s.s. females mated to males of different genotypes. Comparing to the intraspecific cross, a significant reduction of insemination rates was observed in all other crosses. Anopheles marajoara males inefficiently inseminate An. albitarsis s.s. females (P < 0.01) and hybrid males performed even less successfully (P < 0.001). However, insemination rates of crosses using both types of hybrid males did not differ significantly. Although some females inseminated by hybrid males laid eggs, their viability was null (see Table 1).
Figure 3 shows the male reproductive organs of An. albitarsis s.s. and An. marajoara (Figure 3A and B, respectively) and hybrids resulted from interspecific crosses (Figure 3C-E). Invariably, abnormalities were detected only in male hybrids. Malformations consisted of reduced testis with long and thin vas deferens (Figure 3C), abnormal testis with a lobular structure and a malformed and irregular vas deferens (Figure 3D) or fused testis with short vas deferens (Figure 3E).
Variable patterns of reproductive isolation are observed between pairs of siblings of Anopheles species complexes. For example in the An. gambiae complex, two molecularly distinguished forms, known as M and S, do not show reduction in hybrid males and females’ viability or fertility. This result suggests that other mechanisms, apart from postzygotic developmental barriers , are acting in this speciation process. In this sense, it has been claimed that differences in wing-beat frequencies recognition plays a relevant role . Hybridization between Anopheles minimus A and C, from the An. minimus complex, results in fertile males and females with normal ovaries . However, crosses between An. minimus A and E produce sterile males and normal females .
Concerning the postzygotic isolation mechanisms, the hybrid sterility seems to be one of the first outcomes observed in anophelines ,,-. Crosses between An. albitarsis s.s. and An. marajoara produced viable hybrid offspring from both genders in normal rates suggesting the absence of prezygotic isolation mechanisms. However, our results point to a high degree of postzygotic reproductive isolation between these sibling species. The observed pattern of crosses fits the Haldane’s rule, which states that in interspecific crosses the heterogametic sex shows sterility or viability problems before the homogametic one . Eggs resulted from crosses involving hybrid males do not hatch at all. This indicates that the reproductive isolation is a consequence of males’ sterility, possibly due to abnormalities in their reproductive organs. Although hybrid females exhibit a similar degree of insemination when compared to An. albitarsis s.s. control, the hatching percentage of their eggs also suggests postzygotic isolation. Moreover, this isolation phenomenon seems to have a lower effect in females than in males. Interestingly, the development rate of the hybrid females’ offspring is comparable to the intraspecific An. albitarsis s.s. control. This indicates some degree of hybrid females’ infertility, a potential further point of control of interspecific gene flow, at the embryonic stage. In the An. albitarsis complex, previous analysis of reproductive isolation between An. albitarsis s.s. and An. deaneorum, also revealed postzygotic isolation due to sterility of hybrid males . Differing from the present observations, those hybrid males exhibited different levels of sterility. In that case, hybrid males carrying an An. deaneorum X chromosome were the most affected. Moreover, An. deaneorum males were less successful in inseminating An. albitarsis s.s females than their conspecific females, whereas no significant difference was observed in the reciprocal cross. This result suggests that in areas where both species occur in sympatry, asymmetrical introgression might happen. This asymmetry in the hybrid sterility is frequent in interspecific crosses of closely related species . In spite of that, our results suggest an established and effective barrier to this interspecific introgression considering the almost complete failure of both An. albitaris s.s. and An. marajoara males. Further crossing experiments and studies of natural populations using molecular markers might help to determine the genetic relationships and the gene flow among siblings of the Albitarsis complex that includes some of the malaria vectors of America.
Our data confirm that An. albitarsis s.s. and An. marajoara are distinct species of the Albitarsis complex with a high degree of postzygotic reproductive isolation between them. Hybrid males show sterility probably caused by abnormalities in their reproductive organs and a subtle effect was observed in the hybrid females. This result is consistent with the Haldane’s rule which states that in interspecific crosses the heterogametic sex is the first to be affected. The fact that the females are not completely sterile raises the possibility of introgression between these two siblings species.
Bickford D, Lohman DJ, Sodhi NS, Ng PK, Meier R, Winker K, Ingram KK, Das I: Cryptic species as a window on diversity and conservation. Trends of Ecol Evol. 2007, 22: 148-155. 10.1016/j.tree.2006.11.004.
Pfenninger M, Schwenk K: Cryptic animal species are homogeneously distributed among taxa and biogeographical regions. BMC Evol Biol. 2007, 7: 121-10.1186/1471-2148-7-121.
Sinka ME, Bangs MJ, Manguin A, Rubio-Palis Y, Chareonviriyaphap T, Coetzee M, Mbogo CM, Hemingway J, Patil AP, Temperley WH, Gething PW, Kabaria CW, Burkot TR, Harbach RE, HaY SI: A global map of dominant malaria vectors. Parasit Vectors. 2012, 5: 69-10.1186/1756-3305-5-69.
Otarigho B, Falade MO: Molecular evolution and phylogenomics of the Anopheles gambiae complex. J Mosq Res. 2013, 3 (9): 65-70.
Ayala FJ, Coluzzi M: Chromosome speciation: humans, Drosophila, and mosquitoes. Proc Natl Acad Sci U S A. 2005, 102: 6535-6542. 10.1073/pnas.0501847102.
Rosa-Freitas MG, Deane LM, Momen H: A morphological, isoenzymatic and behavioral study of ten populations of Anopheles (Nyssorhynchus) albitarsis Lynch-Arribalzaga, 1878 (Diptera: Culicidae) including from the Type-locality-Baradero, Argentina. Mem Inst Oswaldo Cruz. 1990, 85: 275-289. 10.1590/S0074-02761990000300004.
Conn JE, Wilkerson RC, Segura MNO, De-Sousa RTL, Schlichting CD, Wirtz A, Póvoa MM: Emergence of a new neotropical malaria vector facilitated by human migration and changes in land use. Am J Trop Med Hyg. 2002, 66 (1): 18-22.
Povoa MM, de Souza RTL, Lacerda RNL, Rosa ES, Galiza D, de Souza JR, Wirtz RA, Schlichting CD, Conn JE: The importance of Anopheles albitarsis E and Anopheles darlingi in humam malária transmission in Boa Vista state of Roraima, Brazil. Mem Inst Oswaldo Cruz. 2006, 101: 163-168. 10.1590/S0074-02762006000200008.
Mckeon SN, Lehr MA, Wilkerson RC, Ruiz JF, Sallum MA, Lima JBP, Povoa MM, Conn JE: Lineage divergence detected in the malaria vector Anopheles marajoara (Diptera: Culicidae) in Amazonian Brazil. Malar J. 2010, 9: 271-10.1186/1475-2875-9-271.
Kreutzer RD, Kitzmiller JB, Rabbani MG: Cytogenetically distinguishable sympatric and allopatric populations of the mosquito Anopheles albitarsis. Acta Amazon. 1976, 6 (4): 473-481.
Rosa-Freitas MG, Deane LM: Anopheles (Nyssorhynchus) deaneorum: a new species in the albitarsis complex (Diptera: Culicidae). Mem Inst Oswaldo Cruz. 1989, 84: 289-302. 10.1590/S0074-02761989000300002.
Narang SK, Klein TA, Perera OP, Lima JB, Tang AT: Genetic evidence for the existence of cryptic species in the Anopheles albitarsis complex in Brazil: allozymes and mitochondrial DNA restriction fragment length polymorphisms. Biochem Genet. 1993, 31: 97-112. 10.1007/BF02399823.
Wilkerson RC, Gaffigan TV, Lima JB: Identification of species related to Anopheles (Nyssorhynchus) albitarsis by random amplified polymorphic DNA-polymerase chain reaction (Diptera: Culicidae). Mem Inst Oswaldo Cruz. 1995, 90 (6): 721-732. 10.1590/S0074-02761995000600013.
Wilkerson RC, Parsons TJ, Klein TA, Gaffigan TV, Bergo E, Consolim J: Diagnosis by random amplified polymorphic DNA polymerase chain reaction of four cryptic species related to Anopheles (Nyssorhynchus) albitarsis (Diptera: Culicidae) from Paraguay, Argentina, and Brazil. J Med Entomol. 1995, 32: 697-704. 10.1093/jmedent/32.5.697.
Wilkerson RC, Foster PG, Li C, Sallum MAM: Molecular phylogeny of neotropical Anopheles (Nyssorhynchus) albitarsis species complex (Diptera: Culicidae). An Entomol Soc Amer. 2005, 98: 918-925. 10.1603/0013-8746(2005)098[0918:MPONAN]2.0.CO;2.
Lima JBP, Valle D, Peixoto AA: Analysis of reproductive isolation between sibling species Anopheles albitarsis s.s. and Anopheles deaneorum, two malaria vectors belonging to the Albitarsis Complex (Diptera: Culicidae). J Med Entomol. 2004, 41 (5): 888-893. 10.1603/0022-2585-41.5.888.
Li C, Wilkerson RC: Identification of Anopheles (Nyssorhynchus) albitarsis complex species (Diptera: Culicidae) using rDNA internal transcribed sapacer 2-based polymerase chain reaction primers. Mem Inst Oswaldo Cruz. 2005, 100 (5): 495-500. 10.1590/S0074-02762005000500009.
Li C, Wilkerson RC: Intragenomic rDNA ITS2 variation in the Neotropical Anopheles (Nyssorhynchus) albitarsis complex (Diptera: Culicidae). J Heredity. 2006, 98: 51-59. 10.1093/jhered/esl037.
Merritt TJA, Young CR, Wilkerson RC, Quattro JM: Intron retention identifies a malaria vector within the Anopheles (Nysshorhynchus) albitarsis complex (Diptera: Culicidae). Mol Phyl and Evol. 2005, 35: 719-724. 10.1016/j.ympev.2005.03.009.
Brochero HHL, Li C, Wilkerson RC: A newly recognized species in the Anopheles (Nyssorhynchus) albitarsis complex (Diptera: Culicidae) from Puerto Carreno. Colômbia Am J Trop Med Hyg. 2007, 76 (6): 1113-1117.
Lehr MA, Kilpatrick CW, Wilkerson RC, Conn JE: Cryptic species in the Anopheles (Nyssorhynchus) albitarsis (Diptera: Culicidae) complex: incongruence between random amplified polymorphic DNA-Polymerase Chain Reaction identification and analysis of mitochondrial DNA COI gene sequences. Ann Entomol Soc Am. 2005, 98 (6): 908-917. 10.1603/0013-8746(2005)098[0908:CSITAN]2.0.CO;2.
Motoki MT, Wilkerson RC, Sallum MAM: The Anopheles albitarsis complex with the recognition of Anopheles oryzalimnetes Wilkerson and Motoki, n. sp. and Anopheles janconnae Wilkerson and Sallum, n. sp. (Diptera: Culicidae). Mem Inst Oswaldo Cruz. 2009, 104 (6): 823-850. 10.1590/S0074-02762009000600004.
Sinka ME, Rubio-Palis Y, Manguin S, Patil AP, Temperley WH, Gething PW, Van Boeckel T, Kabaria CW, Harbach RE, Hay SI: The dominat Anopheles vectors of human malaria in the Americas: occurrence data, distribution maps and bionomic précis. Parasit Vectors. 2010, 3: 72-10.1186/1756-3305-3-72.
Ruiz-Lopez F, Wilkerson RC, Conn JE, McKeon SN, Levin DM, Quiñones ML, Póvoa MM, Linton YM: DNA barcoding reveals both known and novel taxa in the Albitarsis Group (Anopheles: Nyssorhynchus) of Neotropical malaria vectors. Parasit Vectors. 2012, 5: 44-10.1186/1756-3305-5-44.
Bourke BP, Foster PG, Bergo ES, Calado DC, Sallum MA: Phylogenetic relationships among species of Anopheles (Nyssorhynchus) (Diptera: Culicidae) based on nuclear and mitochondrial gene sequences. Acta Trop. 2010, 114: 88-96. 10.1016/j.actatropica.2010.01.009.
Besansky NJ, Krzywinski J, Lehmann T, Simard F, Kern M, Mukabayire O, Fontenille D, Touré Y, Sagnon N: Semipermeable species boundaries between Anopheles gambiae and Anopheles arabiensis: Evidence from multilocus DNA sequence variation. Proc Natl AcadSci U S A. 2003, 100: 10818-10823. 10.1073/pnas.1434337100.
Slotman MA, Della Torre A, Calzetta M, Powell JR: Differential introgression of chromosomal regions between Anopheles gambiae and An. arabiensis. Am J Trop Med Hyg. 2005, 73: 326-335.
Djogbénou L, Chandre F, Berthomieu A, Dabiré R, Koffi A, Alout H, Weill M: Evidence of introgression of the ace-1R mutation and of the ace-1 duplication in West African Anopheles gambiae s. s. PLoS One. 2008, 3: e2172-10.1371/journal.pone.0002172.
Klein TA, Lima JBP, Toda-Tang A: Hybridization evidence supporting separate especies status for Anopheles albitarsis and Anopheles deaneorum (Diptera: Culicidae) in Brazil. J Am Mosq Control Assoc. 1991, 7: 301-303.
Horosko S, Lima JBP, Brandolini MB: Establishment of a free mating colony of Anopheles albitarsis from Brazil. J Am Mosq Control Assoc. 1997, 13: 95-96.
Faran ME, Linthicum KJ: A handbook of the Amazon species of Anopheles (Nyssorhynchus) (Diptera: Culicidae). Mosq System. 1981, 13: 1-81.
Jowett T: Preparation of Nucleic Acids. Drosophila: A Practical Approach. Edited by: Roberts DB. 1998, IRL Press, Oxford, 347-371.
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.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.
Clements AN: The Biology of Mosquitoes. Development, Nutrition and Reproduction. Volume I. 1992, Chapmann & Hall, London
Ow-Yang CF, Sta Maria FL, Wharton RH: Maintenance of a laboratory colony of Anopheles maculatus Theobald by artificial mating. Mosq News. 1963, 23: 34-35.
Motulsky HJ: Prism 5 Statistics Guide. GraphPad Software. 2007, Inc, San Diego CA
Lima JBP, Valle D, Peixoto AA: Adaptation of a South American malaria vector to laboratory colonization suggests faster-male evolution for mating ability. BMC Evol Biol. 2004, 4: 12-10.1186/1471-2148-4-12.
Diabaté A, Dabire RK, Milongo N, Lehmann T: Evaluating the effect of postzygotic isolation between molecular forms of Anopheles gambiae (Diptera: Culicidae). J Med Entomol. 2007, 44 (1): 60-64. 10.1603/0022-2585(2007)44[60:ETEOPI]2.0.CO;2.
Pennetier C, Warren B, Dabiré KR, Russell IJ, Gibson G: “Singing on the Wing” as a mechanism for species recognition in the malaria mosquito Anopheles gambiae. Curr Biol. 2010, 20: 131-136. 10.1016/j.cub.2009.11.040.
Choochote W, Rongsriyam Y, Leemingsawat S, Jitpakdi A, Komalamisra N, Surathin K, Somboon P, Chen B, Wongkamchai S, Jariyapan N, Tippawangkosol P, Pitasawat B, Riyong D: Intraspecific hybridization of Anopheles minimus (Diptera: Culicidae) species A and C in Thailand. Southeast Asian J Trop Med Public Health. 2002, 33 (3): 23-28.
Somboon P, Thongwat D, Somwang P, Teng HJ, Tsuda Y, Takagi M: The specific status of Anopheles minimus s.l. collected from Taiwan. Southeast Asian J Trop Med Public Health. 2005, 36 (3): 605-608.
Kitzmiller JB, Frizzi G, Baker RH: Evolution and speciation within the maculipennis complex of the genus Anopheles. Genet Insect Vectors Dis. 1967, 5: 151-210. Elsevier Publishing Company
Kitzmiller JB: Genetics, cytogenetics and evolution of mosquitoes. Adv Genet. 1976, 18: 316-433.
Coluzzi M: Sibling species in Anopheles and their importance in malariology. Misc Publ Entomol Sot Am. 1970, 7: 63-72.
Orr HA: The population genetics of speciation: the evolution of hybrid incompabilities. Genetics. 1995, 139: 1805-1813.
Orr HA: Haldane’s rule. Ann Rev Ecol Syst. 1997, 28: 195-2018. 10.1146/annurev.ecolsys.28.1.195.
Coyne JA, Orr HA: Patterns of speciation revisited. Evolution. 1997, 51: 295-303. 10.2307/2410984.
We would like thank Robson Costa da Silva for technical assistance and Denise Valle for helpful advice and suggestions. This work was supported by grants from CNPq, CAPES and FIOCRUZ. We dedicate this manuscript to Alexandre Afranio Peixoto in memorian.
The authors declare that they have no competing interests.
NGF and ASA performed crossing experiments and analysis, and drafted the manuscript. RVDMA participated in the data analysis and drafting the manuscript. AKRG was involved in the field collection. AAP and JBPL participated in the design of the study, data analysis, drafting the manuscript, general supervision of the research group and funding acquisition. All authors read and approved the final document.