Remarkable diversity of intron-1 of the para voltage-gated sodium channel gene in an Anopheles gambiae/Anopheles coluzzii hybrid zone
© Santolamazza et al.; licensee Biomed Central. 2015
Received: 6 August 2014
Accepted: 20 December 2014
Published: 21 January 2015
Genomic differentiation between Anopheles gambiae and Anopheles coluzzii - the major malaria vectors in sub-Saharan Africa - is localized into large “islands” toward the centromeres of chromosome-X and the two autosomes. Linkage disequilibrium between these genomic islands was first detected between species-specific polymorphisms within ribosomal DNA genes (IGS-rDNA) on the X-chromosome and a single variant at position 702 of intron 1 (Int-1702) of the para Voltage-Gated Sodium Channel (VGSC) gene on chromosome arm 2 L. Intron-1 sequence data from West and Central Africa revealed two clearly distinct and species-specific haplogroups, each characterized by very low polymorphism, which has been attributed to a selective sweep. The aim of this study was to analyse Int-1 sequence diversity in A. gambiae and A. coluzzii populations from the Far-West of their range, in order to assess whether this selective-sweep signature could persist in a zone of high interspecific hybridization.
A 531 bp region of VGSC Int-1 was sequenced in 21 A. coluzzii, 31 A. gambiae, and 12 hybrids from The Gambia and Guinea Bissau, located within the Far-West geographical region, and in 53 A. gambiae s.l. samples from the rest of the range.
Far-West samples exhibit dramatic Int-1 polymorphism, far higher within each country than observed throughout the rest of the species range. Moreover, patterning of haplotypes within A. coluzzii confirms previous evidence of a macro-geographic subdivision into a West and a Central African genetic cluster, and reveals a possible genetic distinction of A. coluzzii populations from the Far-West.
The results suggest a relaxation of selective pressures acting across the VGSC gene region in the hybrid zone. Genetic differentiation in the Far-West could be attributable to a founder effect within A. coluzzii, with subsequent extensive gene flow with secondarily-colonizing A. gambiae, potentially yielding a novel insight on the dynamic processes impacting genetic divergence of these key malaria vectors.
KeywordsMosquito Malaria Hybridization
Mosquito species belonging to the Afro-tropical Anopheles gambiae complex represent a valuable model for studies of ecological speciation (‘speciation with gene flow’) [1,2]. In fact, owing to the major role of some of these species in malaria transmission, their genetic divergence has been studied extensively for more than half century, revealing repeated events of ‘ecotypic speciation’ [3,4]. These studies have revealed the existence of morphologically indistinguishable, but chromosomally/genetically distinct species with very different roles as malaria vectors (A. gambiae sensu stricto, Anopheles arabiensis, Anopheles quadriannulatus, Anopheles amharicus, Anopheles melas, Anopheles merus, Anopheles bwambae) [5,6]. In addition to these sibling species, which are isolated by both pre-and post-mating mechanisms, more recent studies have highlighted a further subdivision within the nominal species A. gambiae s.s., the most synanthropic and efficient malaria vector of the complex [7,8]. These studies have revealed the existence of two taxonomic units initially named the ‘M’ and ‘S’ molecular forms and now formally raised to species as Anopheles coluzzii and A. gambiae s.s. (hereafter referred as A. gambiae) , respectively. The two species are isolated only by (partially understood) pre-mating mechanisms and show limited genomic differentiation, which is localized most prominently in low-recombination peri-centromeric regions of chromosome-X and of the two autosomes [10-13]. These regions have been postulated to represent ‘genomic islands of speciation’, expanding in size by selection across linked loci connected to reproductive isolation, as predicted by ecological speciation with gene flow models [1,2]. However, their large size is probably enhanced by locally reduced recombination [10,14,15], leading to suggestions that high centromeric differentiation is primarily a result of recurrent background selection and hitchhiking unrelated to speciation [2,16]. Interestingly, the species-specific linkage disequilibrium among the three physically-unlinked centromeric regions is widespread , but has been lost at the western extreme of their range (i.e. the ‘Far-West’, from The Gambia to Guinea Bissau). In this putative secondary contact zone between the two species, high frequencies of hybrids are found and pronounced inter-specific differentiation is maintained only on the chromosome-X centromere [17-21].
The first evidence of genetic linkage between the chromosome-X and chromosome-2 genomic island regions came from pre-genomic studies showing that the SNPs within the ribosomal intergenic spacer region (IGS- rDNA), which define the two species , are in strong linkage disequilibrium (LD) with a SNP at position 702 of Intron-1 (hereafter named Int-1702) of the para Voltage-Gated Sodium Channel (VGSC) gene [22,23]. The VGSC gene has been studied extensively because of the presence of two mutations in the exon immediately downstream of Int-1, causing a change from Leucine to Serine (L1014S) or to Phenylalanine (L1014F). Each mutation can confer knock-down resistance (kdr) to DDT and pyrethroid insecticides [22,23]. In West and Central Africa the Int-1702 SNP, which exhibits diagnostic nucleotides for each species (i.e. A. coluzzii = Int-1C; A. gambiae = Int-1T), defines two clearly distinct and species-specific haplotype groups each characterized by very low polymorphism. This led to the hypothesis of a selective sweep, pre-dating kdr mutations, centred on favourable variants in nearby genes that might contribute to the segregation of the two species .
The aim of this study was to investigate sequence diversity of the Int-1 region of the VGSC in A. gambiae and A. coluzzii populations from the Far-West to assess whether the selective-sweep observed in the rest of the species range is maintained in the face of high interspecific gene flow.
Anopheles gambiae sensu lato sample identification and kdr genotyping
Anopheles gambiae s.l. specimens were identified as A. arabiensis, A. quadriannulatus, A. melas or A. merus using the method of Scott et al. , and as A. coluzzii or A. gambiae s.s. by using both PCR-RFLP of the intergenic spacer (IGS) rDNA region  and the SINE-PCR  methods. The latter method is based on an A. coluzzii-specific and irreversible single-locus insertion of a SINE200 retrotransposon in the X-chromosome centromeric region, about 1.5 Mb from IGS region containing A. coluzzii vs A. gambiae species-specific SNPs. The kdr 1014 locus was genotyped by either allele-specific-PCR (AS-PCR) or Hot Oligonucleotide Ligation Assay (HOLA) methods [27-29] in A. coluzzii and A. gambiae. Genotypes of kdr 1014 position in A. arabiensis, A. quadriannulatus, A. melas and A. merus were obtained by direct sequencing of a 531 bp fragment .
Sequencing of Int-1 of VGSC gene
Amplification of a 531 bp region of Int-1 of the VGSC gene (AGAP004707 in genome build AgamP4.1) was carried out using the “Ganest” primer (5′- CAT ACA TTG CTT AAA GCT CTA ATT ATC -3′), located upstream at positions 388–414 in the Int-1 region, coupled with the “Montrev” primer (5′- CAC AAG GCA CAC GAT ACG -3′), located downstream at positions 995–1013 at the end of intron-2 (nucleotide positions as in ). The PCR mixture contained: 10× PCR Buffer (Bioline), 2 mM MgCl2, 200 μM dNTPs equimolar mix, 1 U Taq DNA polymerase (Bioline) and 0.25 μM of each primer, in a total volume of 25 μl. Cycling conditions were as follows: 94°C for 5 min., 35 cycles each with 94°C for 30 sec./ 50°C for 35 sec./ 72°C for 1 min., followed by a final extension step of 10 min at 72°C. Direct sequencing was performed according to  and .
Chromatograms were edited using the Staden Package ver. 2003.1.6  and sequences were aligned using MAFFT ver. 5 . Haplotype inference was performed with the PHASE algorithm  as implemented in DnaSP v5.10.01 . Haplotype sequences reported in Gentile et al.  were also included for comparison. Sequences were deposited in Genbank under accession n° KP300645-KP300752.
DnaSP v. 5.10.01  was used to produce estimates of Int-1 polymorphism and perform Tajima’s D and Fu and Li’s D* and F* [34,35] neutrality tests. Genealogical relationships among Int-1 haplotypes were reconstructed by computing parsimony networks using TCS 1.21 . Haplotypes names were retained for those previously reported in , while novel M1- and S1-related haplotypes were either named with consecutive numbers or with specific codes when exclusive of the Far-West region (i.e. GU for Guinea Bissau, GA for The Gambia, and GUGA for both countries). Frequency and distribution of Int-1 haplotypes in the Far West, West and Central African geographic regions in A. coluzzii, A. gambiae and putative hybrids (by merging original and previous data from ) were also computed. Inter- and intra-specific Fst estimates of genetic differentiation were calculated using Arlequin 3.11 .
A 531 bp region of the VGSC gene Int-1 was sequenced in 21 A. coluzzii, 31 A. gambiae, 12 A. gambiae x A. coluzzii hybrids from the Far-West region. These sequences were aligned with those obtained by Gentile et al.  and with 53 additional original sequences from A. coluzzii and A. gambiae populations from the rest of the distribution range, as well as from other species of the complex (Additional file 1: Table S1), resulting in a total of 490 sequences. Note that A. coluzzii and A. gambiae were identified by both the IGS marker and SINE insertion, which provided completely consistent identification in all samples, with the notable exceptions of those from The Gambia (55% inconsistent) and Guinea Bissau (51% inconsistent). Inconsistently identified individuals were classified based on the SINE marker, as this is not biased toward the hybrid genotype by intra-chromosome recombination known to occur in the IGS multicopy rDNA region . Specimens were genotyped for kdr-mutations: all carried the wild type/insecticide kds susceptible allele of VGSC gene (i.e. TTA) with the exception of two out of eight individuals from Rwanda carrying the TCA allele (hereafter L1014S) in homozygosis and three carrying the TTT allele (hereafter L1014F) in heterozygosis (Additional file 2: Table S2).
Intron-1 of VGSC gene polymorphism and summary statistics in Anopheles coluzzii and A. gambiae samples
Fu & Li’s D
Fu & Li’s F
Data on Int-1 of the VGSC gene highlight dramatic genetic differences between A. coluzzii and A. gambiae populations from West and Central Africa and those from the western extreme of the species range, where inter-specific gene-flow is elevated and inter-specific differentiation reduced.
In West and Central Africa, Int-1 genetic differentiation between A. coluzzii and A. gambiae is high and comparable to that observed in their sibling A. arabiensis (Figure 4). Across this wide geographical region, as previously shown , the two species exhibit low nucleotide/haplotype diversity and are strongly segregated based on two main Int-1 haplotypes (A. coluzzii-M1 and A. gambiae-S1, separated by a single C-T mutational step at site Int-1702) and species-specific rare variants stemming from these. To explain this low polymorphism, a selective sweep centered on a favourable variant in a nearby gene was suggested [22,23]. Since the VGSC gene includes mutations conferring kdr resistance to insecticides, a possible hypothesis is that Int-1 haplotypes are in linkage with these strongly selected alleles. Introgression of kdr alleles and adjacent genomic regions has repeatedly been reported from different geographic areas [19,39]. Indeed, whole genome sequence data have recently shown that over 3 Mb has introgressed from A. gambiae to A. coluzzii in Ghana along with the kdr L1014F mutation . It is possible that selection on kdr-associated Int-1 haplotypes may have had a role in reducing diversity, as kdr-resistance has been reported in some of the A. gambiae analysed populations particularly from West Africa [41,42] [but not in the Far-West region, Pinto et al., unpublished observations; see discussion below]. However, the close physical proximity of kdr and Int-1 would make recombination between resistant and susceptible kdr alleles and their linked Int-1 polymorphisms highly unlikely over a short timescale. Moreover, in A. coluzzii, as the samples analyzed came from populations where kdr-alleles were either absent or present at moderate to low frequencies (i.e. populations from Benin, 40.0%, Nigeria, 19.5% and Cameroon, 6.3% ), and are thus likely to have introgressed very recently [19,22,39]. Recent genomic studies have given additional hints for understanding the reduction in genetic variation at Int-1 and its linkage disequilibrium with markers on the physically unlinked X-centromeric region defining the two species. In fact, as already mentioned, the VGSC gene is located within the chromosome-2 “genomic island” of highest divergence between A. coluzzii and A. gambiae. Under the “speciation island”  scenario, it can be hypothesized that a hitchhiking effect on Int-1 has occurred due to diversifying selection on a chromosome-2 “island” gene participating in the building-up of pre-mating barriers, or conferring differential ecological adaptation and niche segregation between A. gambiae and A. coluzzii. Alternatively, under the “incidental island” scenario [2,13,16], substitutions at Int-1702 may have become fixed after species splitting and accumulated little genetic differentiation due to reduced recombination in the chromosome-2 centromeric region. Interestingly, however, the association between chromosome-X and −2 “islands” is neither observed in Rwanda (Additional file 2: Table S2) nor in Tanzania . In these East African sites, both M1 and S1 haplotypes were found segregating in A. gambiae populations. Thus, if Int-1C represents the ancestral allele in the A. gambiae complex (Figure 1), then Int-1C/T may be considered an ancestral polymorphism retained in A. gambiae populations from East Africa (where A. coluzzii is absent), which became fixed in westward sympatric areas after the splitting of the two species.
In the Far-West region, a strong reduction of inter-specific genetic divergence between A. gambiae and A. coluzzii is found - as indicated by the lower Fst observed in this region (0.14) as opposed to the rest of the range (0.51-0.69) (Figure 4) - and a preferential introgression of M1-related Int-1 haplotypes (“typical” of A. coluzzii) into A. gambiae is observed. These data are consistent with previous studies showing weak association between chromosome-X and −2 centromeric regions and occurrence of asymmetric introgression from A. coluzzii into A. gambiae in the westernmost extreme of their range [17-21]. Furthermore, 15 exclusive Far-West haplotypes were inferred through PHASE  and found interspersed and connected to M1 and S1 geographically widespread variants and to the Far-West-specific M3 (Figure 2). The presence of such private haplotypes might indicate that selective pressures on the chromosome-2 centromere observed in Central African populations (Table 1) are relaxed in the Far-West. Note that, although recombination along the centromeric 500-bp Int-1 fragment analyzed would normally be considered minimal, some reduction in accuracy might occur when reconstructing haplotypes using the PHASE algorithm in the Far-West region, where LD along the 2 L-centromere is known to be lower than in the rest of the species range . However, the PHASE results are supported by summary statistics (Table 1) also indicating extreme Int-1 diversity and recent introgression events in Far-West populations of both species.
There are contrasting possible explanations for the remarkable Int-1 polymorphism observed in the Far-West region. Under the ‘speciation island’ hypothesis  relaxation of diversifying selection on a key isolating trait on chromosome-2 centromeric “island” (of which Int-1 is a part) may have contributed to weaken pre-mating barriers between A. gambiae and A. coluzzii and promote a higher rate of gene-flow in the Far-West region. This hypothesis, however, is in contrast with data from other West and Central African areas, where introgression from A. gambiae to A. coluzzii of a kdr-related genomic portion in linkage with Int-1 does not produce an increase in hybridization rates . Alternatively, the genomic region linked to Int-1 may be not related to speciation [2,16] and the observed pattern in the Far-West region could be attributed to a relaxation of purifying selection operating separately within each species on adaptive genetic traits not directly (or only weakly) involved in reproductive isolation. Hence, following this hypothesis, increased Int-1 polymorphism in the “Far-West” region might be the consequence of an increased recombination rate within the 2 L-centromeric “island” (and Int-1) following disruption of linked (background) selection. Resolution of these competing hypotheses requires assessment of the role of hybridization on the extent of linkage disruption throughout the 2 L-“island” and understanding of whether and how this might affect association with traits critical to speciation.
Finally, the frequency and distribution of Int-1 haplotypes within A. coluzzii across its range provides some hints on further intra-specific geographical patterns (Figure 3). In fact, populations from the West and Far-West regions are characterized by the exclusive presence of haplotype M5, not observed in those from Central Africa. This is consistent with results obtained by other nuclear markers (e.g. microsatellites) showing a macro-geographic subdivision into two distinct West and Central African genetic clusters, corresponding to the forest-savannah biome transition, which may have acted as an ecological barrier to gene flow [26,43,44]. Moreover, the high frequency of the Far-West exclusive M3-haplotype (separated from the major and widespread M1-haplotype by 5 mutational steps) allows speculation that a founder effect (followed by either selection or drift) affected A. coluzzii populations colonizing this region in the past. This last point merits further investigation through a multi-locus approach at a wider genome scale to shed light on the genetic characteristics of source populations originating the A. gambiae/A. coluzzii hybrid zone.
We are grateful to all scientists and entomology teams who contributed to mosquito samples, with particular reference to L. Iyikirenga (National Malaria Control Programme–TRAC Plus, Ministry of Health, Kigali, Rwanda), T.G.T. Jaenson (Medical Entomology Unit, Department of Systematic Biology, Uppsala University, Uppsala, Sweden) and K. Palsson (KTH Royal Institute of Technology, Stockholm, Sweden). This work was supported by European Union’s INFRAVEC project (grant agreement no. 228421 under FP7 program) to AdT and JP, Ricerca Scientifica 2012 grant by Università di Roma SAPIENZA to AdT and MIUR-FIRB “Futuro in Ricerca 2010” grant to BC (Grant N° RBFR106NTE) and by CIRM ISS - Italian Malaria Network.
- Feder JL, Egan SP, Nosil P. The genomics of speciation-with-gene-flow. Trends Genet. 2012;28:342–50.PubMedView ArticleGoogle Scholar
- Turner TL, Hahn MW. Genomic islands of speciation or genomic islands and speciation? Mol Ecol. 2010;19:848–50.PubMedView ArticleGoogle Scholar
- Coluzzi M. Spatial distribution of chromosomal inversions and speciation in anopheline mosquitoes. In: Barigozzi C, editor. Mechanisms of Speciation. New York: Alan R. Liss Inc; 1982. p. 143–53.Google Scholar
- Manoukis NC, Powell JR, Touré MB, Sacko A, Edillo FE, Coulibaly MB, et al. A test of the chromosomal theory of ecotypic speciation in Anopheles gambiae. Proc Natl Acad Sci USA. 2008;105:2940–5.PubMed CentralPubMedView ArticleGoogle Scholar
- Coluzzi M, Sabatini A, della Torre A, Di Deco MA, Petrarca V. A polytene chromosome analysis of the Anopheles gambiae species complex. Science. 2002;298:1415–8.PubMedView ArticleGoogle Scholar
- Ayala FJ, Coluzzi M. Chromosome speciation: humans, Drosophila, and mosquitoes. Proc Natl Acad Sci USA. 2005;102:6535–42.PubMed CentralPubMedView ArticleGoogle Scholar
- della Torre A, Fanello C, Akogbeto M, Dossou-yovo J, Favia G, Petrarca V, et al. Molecular evidence of incipient speciation within Anopheles gambiae s.s. in West Africa. Insect Mol Biol. 2001;10:9–18.PubMedView ArticleGoogle Scholar
- Della Torre A, Costantini C, Besansky NJ, Caccone A, Petrarca V, Powell JR, et al. Speciation within Anopheles gambiae: the glass is half full. Science. 2002;298:115–7.PubMedView ArticleGoogle Scholar
- Coetzee M, Hunt RH, Wilkerson R, della Torre A, Coulibaly MB, Besansky NJ. Anopheles coluzzii and Anopheles amharicus, new members of the Anopheles gambiae complex. Zootaxa. 2013;3619:246–74.View ArticleGoogle Scholar
- Turner TL, Hahn MW, Nuzhdin SV. Genomic islands of speciation in Anopheles gambiae. PLoS Biol. 2005;3:e285.PubMed CentralPubMedView ArticleGoogle Scholar
- Lawniczak MK, Emrich SJ, Holloway AK, Regier AP, Olson M, White B, et al. Widespread divergence between incipient Anopheles gambiae species revealed by whole genome sequences. Science. 2010;330:512–4.PubMed CentralPubMedView ArticleGoogle Scholar
- Neafsey DE, Lawniczak MK, Park DJ, Redmond SN, Coulibaly MB, Traoré SF, et al. SNP genotyping defines complex gene-flow boundaries among African malaria vector mosquitoes. Science. 2010;330:514–7.PubMedView ArticleGoogle Scholar
- White BJ, Cheng C, Simard F, Costantini C, Besansky NJ. Genetic association of physically unlinked islands of genomic divergence in incipient species of Anopheles gambiae. Mol Ecol. 2010;19:925–39.PubMed CentralPubMedView ArticleGoogle Scholar
- Via S, West J. The genetic mosaic suggests a new role for hitchhiking in ecological speciation. Mol Ecol. 2008;17:4334–45.PubMedView ArticleGoogle Scholar
- Feder JL, Nosil P. The efficacy of divergence hitchhiking in generating genomic islands during ecological speciation. Evolution. 2010;64:1729–47.PubMedView ArticleGoogle Scholar
- Cruickshank TE, Hahn MW. Re-analysis suggests that genomic islands of speciation are due to reduced diversity, not reduced gene flow. Mol Ecol. 2014;23:3133–57.PubMedView ArticleGoogle Scholar
- Caputo B, Santolamazza F, Vicente JL, Nwakanma DC, Jawara M, Palsson K, et al. The “far-west” of Anopheles gambiae molecular forms. PLoS One. 2011;6:e16415.PubMed CentralPubMedView ArticleGoogle Scholar
- Nieman CC, Sanford MR, Dinis J, Martins C, Rodrigues A, Cornel AJ, et al. Asymmetric introgression between the M and S forms of the malaria vector, Anopheles gambiae, maintains divergence despite extensive hybridization. Mol Ecol. 2011;20:4983–94.PubMed CentralPubMedView ArticleGoogle Scholar
- Weetman D, Wilding CS, Steen K, Pinto J, Donnelly MJ. Gene flow-dependent genomic divergence between Anopheles gambiae M and S forms. Mol Biol Evol. 2012;29:279–91.PubMed CentralPubMedView ArticleGoogle Scholar
- Nwakanma DC, Neafsey DE, Jawara M, Adiamoh M, Lund E, Rodrigues A, et al. Breakdown in the process of incipient speciation in Anopheles gambiae. Genetics. 2013;193:1221–31.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee Y, Marsden CD, Norris LC, Collier TC, Main BJ, Fofana A, et al. Spatiotemporal dynamics of gene flow and hybrid fitness between the M and S forms of the malaria mosquito, Anopheles gambiae. Proc Natl Acad Sci USA. 2013;110:19854–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Weill M, Chandre F, Brengues C, Manguin S, Akogbeto M, Pasteur N, et al. The kdr mutation occurs in the Mopti form of Anopheles gambiae s.s. through introgression. Insect Mol Biol. 2000;9:451–5.PubMedView ArticleGoogle Scholar
- Gentile G, Santolamazza F, Fanello C, Petrarca V, Caccone A, della Torre A. Variation in an intron sequence of the voltage-gated sodium channel gene correlates with genetic differentiation between Anopheles gambiae s.s. molecular forms. Insect Mol Biol. 2004;13:371–7.PubMedView ArticleGoogle Scholar
- Scott JA, Brogdon WG, Collins FH. Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg. 1993;49:520–9.PubMedGoogle Scholar
- Fanello C, Santolamazza F, della Torre A. Simultaneous identification of species and molecular forms of the Anopheles gambiae complex by PCR-RFLP. Med Vet Entomol. 2002;16:461–4.PubMedView ArticleGoogle Scholar
- Santolamazza F, Mancini E, Simard F, Qi Y, Tu Z, della Torre A. Insertion polymorphisms of SINE200 retrotransposons within speciation islands of Anopheles gambiae molecular forms. Malar J. 2008;7:163.PubMed CentralPubMedView ArticleGoogle Scholar
- Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Berge JB, Devonshire AL, et al. Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998;7:179–84.PubMedView ArticleGoogle Scholar
- Ranson H, Jensen B, Vulule JM, Wang X, Hemingway J, Collins FH. Identification of a point mutation in the voltage-gated sodium channel gene of Kenyan Anopheles gambiae associated with resistance to DDT and pyrethroids. Insect Mol Biol. 2000;9:491–7.PubMedView ArticleGoogle Scholar
- Lynd A, Ranson H, McCall PJ, Randle NP, Black WC, Walker ED, et al. A simplified high-throughput method for pyrethroid knock-down resistance (kdr) detection in Anopheles gambiae. Malar J. 2005;4:16.PubMed CentralPubMedView ArticleGoogle Scholar
- Staden R, Beal KF, Bonfield JK. The Staden package. Methods Mol Biol. 2000;132:115–30.PubMedGoogle Scholar
- Katoh K, Kuma K, Toh H, Miyata T. MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005;33:511–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Stephens M, Smith N, Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001;68:978–89.PubMed CentralPubMedView ArticleGoogle Scholar
- Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.PubMedView ArticleGoogle Scholar
- Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.PubMed CentralPubMedGoogle Scholar
- Fu YX, Li WH. Statistical tests of neutrality of mutations. Genetics. 1993;133:693–709.PubMed CentralPubMedGoogle Scholar
- Clement M, Posada D, Crandall KA. TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000;9:1657–9.PubMedView ArticleGoogle Scholar
- Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform Online. 2005;1:47–50.PubMed CentralGoogle Scholar
- Santolamazza F, Caputo B, Calzetta M, Vicente JL, Mancini E, Petrarca V, et al. Comparative analyses reveal discrepancies among results of commonly used methods for Anopheles gambiae molecular form identification. Malar J. 2011;10:215.PubMed CentralPubMedView ArticleGoogle Scholar
- Etang J, Vicente JL, Nwane P, Chouaibou M, Morlais I, Do Rosario VE, et al. Polymorphism of intron-1 in the voltage-gated sodium channel gene of Anopheles gambiae s.s. populations from Cameroon with emphasis on insecticide knockdown resistance mutations. Mol Ecol. 2009;18:3076–86.PubMedView ArticleGoogle Scholar
- Clarkson CS, Weetman D, Essandoh J, Yawson AE, Maslen G, Manske M, et al. Adaptive introgression between Anopheles sibling species eliminates a major genomic island but not reproductive isolation. Nature Comm. 2014;5:4248.View ArticleGoogle Scholar
- Pinto J, Lynd A, Vicente JL, Santolamazza F, Randle NP, Moreno G, et al. Multiple origins of knockdown resistance mutations in the Afrotropical mosquito vector Anopheles gambiae. PLoS One. 2007;2:e1243.PubMed CentralPubMedView ArticleGoogle Scholar
- Santolamazza F, Calzetta M, Etang J, Barrese E, Dia I, Caccone A, et al. Distribution of knock-down resistance mutations in Anopheles gambiae (Diptera: Culicidae) molecular forms in west and west-central Africa. Malar J. 2008;7:74.PubMed CentralPubMedView ArticleGoogle Scholar
- Slotman MA, Tripet F, Cornel AJ, Meneses CR, Lee Y, Reimer LJ, et al. Evidence for subdivision within the M molecular form of Anopheles gambiae. Mol Ecol. 2007;16:639–49.PubMedView ArticleGoogle Scholar
- Pinto J, Egyir-Yawson A, Vicente J, Gomes B, Santolamazza F, Moreno M, et al. Geographic population structure of the African malaria vector Anopheles gambiae suggests a role for the forest-savannah biome transition as a barrier to gene flow. Evol Appl. 2013;6:910–24.PubMed CentralPubMedView ArticleGoogle Scholar
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.