<title>Background</title> <p>Malaria remains one of the world's major health problems claiming at least one million deaths each year in Africa [##UREF##0##1##]. In the forested areas of equatorial Africa, where malaria transmission occurs all year long, <italic>Anopheles moucheti </italic>mosquitoes can sustain malaria transmission intensities as high as 100–300 infected bites per man per year in villages located at the vicinity of large rivers and slow-moving streams where its larvae develop [##UREF##1##2##, ####UREF##2##3##, ##REF##15780341##4##, ##REF##17162956##5####17162956##5##]. However, despite playing such an important epidemiological role in malaria transmission, this group of mosquitoes remains poorly known and insufficiently studied. Data on its bionomics and genetic structure are currently lacking although such data are of paramount importance for a comprehensive implementation and monitoring of malaria vector control in Central Africa [##REF##15225985##6##].</p> <p><italic>Anopheles moucheti </italic>is a group of three morphological forms: <italic>A. moucheti moucheti, A. moucheti nigeriensis </italic>and <italic>A. moucheti bervoetsi</italic>, which can be distinguished by minor variations in the size and distribution of pale fringe spots and pale vein spots on the wings at the adult stage, and at the larval stage by the number of branches of the saddle hair (&gt;5 branches for <italic>A. m. bervoetsi </italic>and &lt;5 branches for the two others) [##UREF##3##7##,##UREF##4##8##]. However, population genetics studies using allozyme markers revealed that these morphological variations were not segregating between the different taxonomic units that build up the <italic>A. moucheti </italic>group in Central Africa and were therefore of poor diagnostic value [##REF##12798002##9##]. More recently, DNA sequence differences were detected in the mitochondrial gene encoding the cytochrome B (CytB) and the ribosomal DNA Internal Transcribed Spacers (ITS) 1 and 2 between specimens of <italic>A. m. moucheti</italic>, <italic>A. m. nigeriensis </italic>and <italic>A. m. bervoetsi </italic>and a diagnostic PCR assay was subsequently developed allowing straightforward identification of the three taxonomic units within the <italic>A. moucheti </italic>group [##REF##17550437##10##]. This study further suggested that <italic>A. m. moucheti </italic>is widespread throughout the forested areas of Central Africa, whereas <italic>A. m. nigeriensis </italic>and <italic>A. m. bervoetsi </italic>were found only in their type localities in Nigeria and the Democratic Republic of Congo (DRC), respectively. Microsatellite DNA markers have been isolated from <italic>A. moucheti </italic>[##UREF##5##11##] and these were demonstrated to be suitable tools for population genetics studies within this group of mosquitoes [##REF##17227668##12##]. Very low levels of genetic differentiation (Fst&lt;0.0275) were detected between <italic>A. moucheti </italic>populations situated 65–400 km apart in Cameroon, suggesting high levels of gene flow at this geographical scale [##REF##17227668##12##].</p> <p>Previous findings from Cameroon [##REF##17227668##12##] are expanded through the inclusion of mosquitoes sampled in DRC and Uganda, to further explore the level of genetic structuring between populations of the <italic>A. moucheti </italic>group and to precise the taxonomic status of <italic>An. m. bervoetsi</italic>. Analytical methods, based on various aspects of the data, are used to provide insights into the role and relative importance of geographic distance, demographic parameters (<italic>eg </italic>effective population size and demographic instability) and natural barriers to gene flow such as habitat discontinuities and speciation in shaping the observed population structure.</p>
<title>Methods</title> <title>Mosquitoes sampling and collection sites</title> <p>The mosquito samples obtained from four villages in Cameroon that were used in this study were described in details previously [##REF##17227668##12##]. Additional adult mosquitoes were collected by pyrethrum spray catches and/or bednet traps from two villages in DRC including the type locality of <italic>A. m. bervoetsi</italic>, Tsakalakuku (5°51'S; 17°23'E) and Kenge (5°19'S; 19°58'E); and from the island of Bufumira (0°19'S; 32°22'E) on Lake Victoria in Uganda (Figure ##FIG##0##1##). Collections were attempted in Nigeria, in and around the village of Akaka (6°27'N; 3°24'E) in the Lagos area where <italic>A. m. nigeriensis </italic>was originally described [##UREF##3##7##], but were unsuccessful.</p> <p>Tsakalakuku is situated in the tropical wet savannas region of Africa. The area is characterized by a succession of hills covered with grass fields and valleys occupied by forest galleries along streams. The climate comprises a dry and a rainy season of six months each. All other collection sites are located within the Congo-Guinean phytogeographic zone, characterized by a typical equatorial climate with two rainy seasons extending from March to June and from September to November (total rainfall around 1,500 mm per year).</p> <p>Mosquito collections were conducted from July 2003 to April 2004 in Cameroon [##REF##17227668##12##], in December 2003 in DRC and in October 2004 in Uganda. <italic>Anopheles moucheti </italic>specimens were visually sorted from other anophelines according to morphological identification keys [##UREF##3##7##,##UREF##6##13##]. All specimens were stored individually and kept at -20°C until further analysis.</p> <title>DNA extraction and genotype scoring</title> <p>Genomic DNA was extracted from wings or legs of each individual mosquito as described earlier [##REF##17227668##12##]. Morphological identification was confirmed through the recently developed PCR based assay [##REF##17550437##10##]. Genotypes at 10 microsatellite loci were determined for the DRC (N = 64 in Tsakalakuku and N = 11 in Kenge) and Uganda (N = 57) samples as previously described [##REF##17227668##12##].</p> <title>Data analysis</title> <p>Genetic diversity within samples and overall was measured at each locus by estimating allele richness <italic>Rs</italic>, an unbiased estimator of the number of alleles in each sample accounting for differences in sample sizes [##UREF##7##14##], and <italic>He </italic>[##UREF##8##15##], the unbiased expected heterozygosity under Hardy-Weinberg equilibrium (HWE), using the software FSTAT V2.9.3.2 [##UREF##9##16##]. Genotypic frequencies were tested against HWE for each locus in the pooled population and in each sample. Statistical significance was assessed by the exact probability test available in GENEPOP V3.2 [##UREF##10##17##]. Linkage disequilibrium between loci was tested by exact tests on contingency tables, also available in GENEPOP.</p> <p>Genetic differentiation between populations was assessed by estimating Wright's F-statistics [##REF##748322##18##], calculated according to Weir &amp; Cockerham [##UREF##11##19##]. Statistical significance of Fst was assessed using G-based exact tests for genotypic differentiation [##REF##8978076##20##], available in GENEPOP. The correlation between genetic and geographic distances, assuming isolation by distance, was assessed by the regression of pairwise Fst estimates on the logarithm (ln) of geographic distances between sampling sites [##REF##9093870##21##], and tested by the Mantel test available in GENEPOP. A Bayesian approach was further implemented to infer the number of genetic clusters (<italic>K</italic>) in the dataset without prior information on the sampling locations, using STRUCTURE 2.2 [##REF##10835412##22##]. A model where the allele frequencies were correlated within populations was assumed (λ was set at 1, the default value). The software was run with the option of admixture, allowing for some mixed ancestry within individuals, and α was allowed to vary. Twenty independent runs were done for each value of <italic>K </italic>(K = 1 to 8), with a burn-in period of 100,000 iterations and 100,000 replications. The method of Evanno <italic>et al </italic>[##REF##15969739##23##] was used to determine the most likely number of clusters. This approach uses an <italic>ad hoc </italic>quantity, Δ<italic>K</italic>, based on the second order rate of change of the likelihood function between successive values of <italic>K</italic>.</p> <p>Because demographic instability such as recent population bottleneck and/or expansion might bias genetic differentiation estimates to a significant extent [##REF##10589511##24##,##REF##11420373##25##], heterozygosity tests were used to test for Mutation-Drift Equilibrium (MDE) within each sample, as implemented in BOTTLENECK 1.2.02 [##REF##8978083##26##]. At selectively neutral loci, the expected heterozygosities calculated from allele frequencies data (<italic>He</italic>) and from the number of alleles and sample sizes (<italic>Heq</italic>) are expected not to be significantly different in a population at MDE. Comparing <italic>He </italic>to <italic>Heq </italic>across loci, therefore, provides the basis for testing this hypothesis. If the population recently experienced a bottleneck, rare alleles will be rapidly lost and therefore <italic>Heq </italic>will decrease faster than <italic>He</italic>. Thus, observing a significant number of loci with <italic>He&gt;Heq </italic>suggests the focal population recently experienced a bottleneck while the reverse trend (i.e., <italic>He</italic>&lt;<italic>Heq</italic>) may suggest population expansion. Estimates of expected heterozygosity under MDE were calculated assuming a Stepwise Mutation Model (SMM) and a Two Phase Model (TPM) with 10–30% indels larger than the repeat unit. Statistical significance of the deviation from MDE was assessed for each sample across all loci by the Wilcoxon signed ranks tests and sign tests available in BOTTLENECK.</p> <p>Differences in effective population size (Ne) between samples might further increase estimates of genetic differentiation between populations because differences in Ne violates assumptions of the island model of population migration, assumed to hold true when devising F-statistics [##UREF##11##19##]. Estimates of "long-term" effective population size [##UREF##8##15##] were calculated for each sample based on the expected heterozygosity at each microsatellite locus assuming a SMM using the formula Neμ = {[1/(1-He)]<sup>2</sup>-1}/8 [##UREF##8##15##,##REF##4777279##27##], where He is the expected heterozygosity under HWE and μ is the microsatellite mutation rate.</p> <p>Because the average mutation rate does not vary much even between well separated species such as pigs (7 × 10<sup>-5</sup>, [##REF##7626894##28##]) and mice (4.5 × 10<sup>-5</sup>, [##REF##1353738##29##]), the value of 10<sup>-4 </sup>proposed by Lehmann <italic>et al </italic>[##REF##9501493##30##] for <italic>Anopheles gambiae </italic>was conservatively adopted for estimating <italic>A. moucheti </italic>long-term Ne. Nevertheless, inferences were drawn on a relative scale, using the product of Neμ as a proxy of long-term Ne for each population, therefore alleviating any bias due to incorrect estimation of the mutation rate.</p> <p>In all instances where multiple tests were conducted simultaneously, the sequential Bonferroni procedure [##UREF##12##31##] was applied to adjust the nominal significance level.</p>
<title>Results</title> <title>Genetic variability within populations</title> <p>A total of 355 mosquitoes of the <italic>A. moucheti </italic>group were analysed in this study, including 223 mosquitoes from Cameroon that were previously genotyped [##REF##17227668##12##]. Of these, 64 specimens collected in Tsakalakuku (DRC) were <italic>A. m. bervoetsi </italic>and the rest were <italic>A. m. moucheti</italic>. Genotypes at 10 microsatellites were determined. All loci amplified successfully in all <italic>A. m. moucheti </italic>populations and were highly polymorphic, showing between 11 (AM13) and 17 (AM5 and AM15) distinct alleles. By contrast, only six loci could be amplified in the <italic>A. m. bervoetsi </italic>sample, among which AM13 was not polymorphic (Table ##TAB##0##1##). Consequently, this population showed the lowest allele richness (even when only five polymorphic loci were considered), and lowest expected heterozygosity (Table ##TAB##0##1##), followed by the <italic>A. m. moucheti </italic>population collected from the island of Bufurima (Uganda). All mainland populations of <italic>A. m. moucheti </italic>showed similar average allele richness (range 5.53–6.66) and expected heterozygosity (range 0.771–0.833).</p> <p>Hardy-Weinberg expectations were significantly rejected (P &lt; 0.001) for seven out of 10 loci when considering the pooled samples as belonging to a single panmictic population, with heterozygote deficits being evidenced at all loci, as expected when different gene pools are mixed. At the population level, 22 out of 65 tests did not conform to Hardy-Weinberg expectations after the multi-test analysis was taken into account. Significant deviation from HWE varied across loci in a population-dependent manner. The Uganda population from Bufumira island had the highest number of loci in departure from HWE (6 of 10) while the Kenge population had the fewest (1 of 10).</p> <p>Exact tests for linkage disequilibrium within each of the seven populations resulted in three significant values out of 276 comparisons after correction by the Bonferroni procedure (two in Mouloundou (AM2-AM6, AM9-AM20) and one in Simbock (AM9-AM10)). No pair of loci appeared in linkage disequilibrium in more than one population, suggesting genetic independence between loci. When the test was performed in the pooled populations, two pairs of loci (AM2-AM16 and AM2-AM20) out of 45 possible combinations showed highly significant P values (&lt;10<sup>-6</sup>).</p> <title>Genetic differentiation between populations</title> <p>Table ##TAB##1##2## shows Fst estimates for all pairwise population comparisons. Genotypic frequencies were highly significantly different among samples (G-test, P &lt; 0.001). Low to moderate levels of genetic differentiation were measured among mainland populations of <italic>A. m. moucheti </italic>from Cameroon and DRC, with mean Fst estimates ranging 0.009–0.049 (P &lt; 0.001). The island population of Bufurima (Uganda) showed higher levels of differentiation with this core group, with pairwise Fst estimates in the range 0.167–0.223 (P &lt; 0.001). The highest levels of genetic differentiation were observed in all comparisons involving the <italic>A. m. bervoetsi </italic>sample collected in Tsakalakuku (DRC), with Fst estimates ranging 0.343–0.448 (P &lt; 0.001).</p> <p>In agreement with results based on Fst, the Bayesian cluster analysis showed that the most likely K value identified was K = 3. This corresponded to three distinct genetic clusters: (1) mainland <italic>A. m. moucheti </italic>from Cameroon and DRC, (2) <italic>A. m. moucheti </italic>from Bufumira island in Uganda and (3) <italic>A. m. bervoetsi </italic>(Figure ##FIG##1##2##).</p> <p>Because isolation by distance is likely to play a major role in shaping the distribution of genetic diversity across continuous habitats [##REF##9093870##21##,##UREF##13##32##], only the mainland populations of <italic>A. m. moucheti </italic>from Cameroon and DRC (<italic>i.e</italic>. excluding the samples from the island of Bufurima and the <italic>A. m. bervoetsi </italic>sample from Tsakalakuku and focusing on "cluster 1" described above) were used for the Mantel test. Positive and highly significant correlation (P &lt; 0.008, Mantel test) was found between genetic (Fst) and geographic distances. Using the equation of the regression line of Fst on the logarithm of distance between sampling sites (Figure ##FIG##2##3##), the expected level of genetic differentiation between all mainland samples of <italic>A. m. moucheti </italic>and the <italic>A. m. moucheti </italic>sample from Bufurima island and <italic>A. m. bervoetsi</italic>, respectively were predicted under the hypothesis that geographic distance between populations was the main determinant of genetic differentiation. As can be graphically seen on Figure ##FIG##2##3##, the predicted Fst estimates were three to nine folds lower than the observed value for the Bufurima sample, and up to 8–25 folds lower than the observed value for the <italic>A. m. bervoetsi </italic>population. As such, distance alone contributed to less than 30% of the observed level of differentiation between the island <italic>A. m. moucheti </italic>sample from Bufurima and all other <italic>A. m. moucheti </italic>samples, while it explained less than 15% of the differentiation with the <italic>A. m. bervoetsi </italic>sample.</p> <title>Effective population size and demographic stability</title> <p>Estimates of long-term Ne were similar for all mainland <italic>A. m. moucheti </italic>populations (Table ##TAB##2##3##). These were significantly lower for the Bufurima island population and the <italic>A. m. bervoetsi </italic>sample. Calculation of the relative ratio of Neμ compared to the <italic>A. m. moucheti </italic>population with the smallest effective population size (Bufumira island), showed that the estimates were five to nine folds higher for mainland populations of <italic>A. m. moucheti</italic>, whereas they were at least five folds lower for the <italic>A. m. bervoetsi </italic>sample. This demonstrates significant heterogeneity in effective population size between the different genetic clusters identified above.</p> <p>Estimates of genetic differentiation and effective population size however are based on the assumption of MDE. Results of the heterozygosity tests (Table ##TAB##3##4##) did not reveal any evidence for departure from MDE in any of the mainland populations of <italic>A. m. moucheti</italic>, nor in the <italic>A. m. bervoetsi </italic>population. However, a consistent trend for lower-than-expected heterozygosity (i.e., He&lt;Heq) was evidenced for the Bufurima island population, suggesting recent demographic expansion.</p>
<title>Discussion</title> <p>In this study, six <italic>A. m. moucheti </italic>populations from different geographic locations and one <italic>A. m. bervoetsi </italic>population were compared for variation in polymorphism and allele distribution at 10 microsatellite loci. Successful amplification at each microsatellite locus was obtained for all <italic>A. m. moucheti </italic>specimens while only six loci could be amplified in the <italic>A. m. bervoetsi </italic>sample, one of which did not show any polymorphism, all specimens investigated showing the same single allele at a homozygous state. This result provides further support for speciation within the <italic>A. moucheti </italic>group of malaria vectors in Central Africa and reflects, for the first time, genome-wide differentiation between <italic>A. m. moucheti </italic>and <italic>A. m. bervoetsi</italic>. Indeed, although the exact cytological location of the microsatellite markers is not known yet, linkage disequilibrium analysis revealed no evidence for genetic linkage between loci, suggesting they provided independent replicates for genome-wide estimation of genetic differentiation between samples. Successful amplification of microsatellite alleles was demonstrated among closely related species such as humans and great apes [##REF##7510853##33##] as well as between sibling species of wasps [##REF##9561790##34##] and members of anophelines species complexes [##REF##9675871##35##,##REF##10403334##36##]. However, the proportion of loci developed for one species that can amplify in another decreases rapidly with increasing evolutionary distance [##REF##1901092##37##,##REF##8688957##38##]. These results are, therefore, in straight agreement with previous studies based on morphological [##UREF##3##7##,##UREF##4##8##] and molecular data (mtDNA CytB and rDNA ITS; [##REF##17550437##10##]), prompting for elevation of <italic>A. m. bervoetsi </italic>to full specific rank as a closely related sibling species of <italic>A. m. moucheti</italic>.</p> <p>Linkage disequilibrium analysis further ruled out the hypothesis that the fairly high number of loci which were found out of HWE in several collections was indicative of inbreeding and/or population subdivision (within samples). If this was the case, genome-wide signatures of departure from HWE and high linkage disequilibrium between loci should be evidenced, because members of the different sub-populations would have different probabilities to carry certain combinations of alleles [##REF##9501493##30##]. Such trends were not observed in the dataset, suggesting null alleles, rather than population subdivision may be responsible for the deviations observed. Null alleles are a common finding in anophelines' population genetics studies [##REF##10583542##39##, ####REF##12721225##40##, ##REF##15638948##41####15638948##41##]. Because the frequencies of such null alleles might differ between sub-populations, they contribute to the overall genetic differentiation between populations. Fst estimates between populations were therefore calculated using all the information available from all loci and all samples.</p> <p>Fst estimates recorded between <italic>A. m. bervoetsi </italic>and each of <italic>A. m. moucheti </italic>populations were very high and statistically significant (Fst&gt;0.34, P &lt; 0.001), falling in the upper range of values reported between well separated anophelines sibling species using various molecular markers [##REF##9409838##42##, ####REF##17038689##43##, ##REF##17257119##44####17257119##44##]. This result, as well as results from the Bayesian analysis clearly identified <italic>A. m. bervoetsi </italic>as a genetically distinct entity within the <italic>A. moucheti </italic>group. Accordingly, it seems reasonable to consider this taxon as a full, independently evolving species within the <italic>A. moucheti </italic>group and, henceforth, to refer to this species as <italic>Anopheles bervoetsi</italic>. However, considering that <italic>A. bervoetsi </italic>has never been reported to occur in sympatry with <italic>A. moucheti s.s</italic>., nor outside of its type locality, this assertion should be validated through traditional crossing experiments, which are yet impossible to implement because members of the <italic>A. moucheti </italic>group have never been maintained successfully under insectary conditions. Preliminary analysis of 237 field-collected <italic>A. bervoetsi </italic>specimens after ELISA detected three females infected by <italic>Plasmodium falciparum </italic>(Antonio-Nkondjio C, Ndo C, Awono-Ambene HP and Simard F, unpublished). Although incrimination of this species as a malaria vector through dissection of its salivary gland still has to be processed, this points to a possible and previously unrecognized role of this mosquito in malaria transmission in Central Africa.</p> <p>Significant isolation by distance was revealed between mainland <italic>A. m. moucheti </italic>populations from Cameroon and DRC, separated by distances &gt;1,000 km (Figure ##FIG##0##1##), suggesting continuous habitat suitability for <italic>A. moucheti </italic>in these forested environments. Extrapolating the level of differentiation expected under the sole influence of geographical distance between mainland populations and the <italic>A. m. moucheti </italic>population collected on the island of Bufurima on Lake Victoria (Uganda) showed that the observed level of differentiation was three to nine folds higher than expected. Such high Fst estimates probably reflect the contribution of large water bodies separating this island population from mainland ones, acting as a barrier to gene flow by restricting opportunities for migration between populations, as was demonstrated for <italic>A. gambiae </italic>in this area [##REF##15638948##41##] and elsewhere [##REF##10403334##36##,##REF##16313568##45##,##REF##17937805##46##]. Moreover, significant differences in effective population sizes (Ne) were demonstrated, the island population of Bufurima showing significantly lower Ne than its mainland counterparts, and these might further increase Fst estimates [##REF##6614868##47##]. Although Ne estimates based on He are criticizable because they rely on a number of assumptions including correct estimation of microsatellite mutation rate and mutation model, populations at MDE and selective neutrality of the loci, the comparison of Neμ on a relative scale allowed relieving some of these assumptions. Lower effective population size on the islands of Lake Victoria was indeed demonstrated for <italic>A. gambiae </italic>compared to neighbouring mainland populations [##REF##15638948##41##] and the results presented here suggest the same probably applies within the <italic>A. moucheti </italic>group. As mentioned above however, estimates of Ne derived from He are sensitive to deviation from MDE. No significant deviation from MDE was evidenced within <italic>A. m. moucheti </italic>populations, although, to some extent, trends for recent population expansion were revealed for the island population of Bufurima. Colonization of Lake Victoria islands by anthropophilic malaria vectors probably followed initial human settlements on these islands in the early 1900s [##REF##15638948##41##,##UREF##14##48##]. Initial founder effect might have occurred at that time but experimental studies and simulations have shown such events are only detectable for a small number of generations before a new equilibrium is reached, especially when considering molecular markers with high mutation rates such as microsatellite loci [##REF##8978083##26##,##REF##11050547##49##]. Inferences suggesting population expansion are generally more robust but still need to be ascertained through in-depth investigations [##REF##11420373##25##].</p> <p>Clearly, a more comprehensive picture of the genetic structure and distribution of genetic diversity within and among natural populations of members of the <italic>A. moucheti </italic>group of malaria vectors would have been obtained with the inclusion in this study of specimens of <italic>A. m. nigeriensis</italic>. However, as mentioned above, collections conducted in and around the type locality of this species were unsuccessful. Earlier investigations allowed collection of a few representative specimens [##REF##17550437##10##] but sample sizes were far too small to allow reliable microsatellite allelic frequencies assessment. The Lagos area has recently undergone significant levels of anthropogenic environmental reshaping and urban expansion [##UREF##15##50##] and this might have led to a significant drop in <italic>A. m. nigeriensis </italic>populations, as was observed for <italic>A. m. moucheti </italic>in areas of southern Cameroon [##REF##15780341##4##]. As formerly highlighted [##REF##15225985##6##], the availability of PCR-based diagnostic tools and other molecular markers and their increased use in routine entomological surveys might allow more refined assessments of the diversity, geographic distribution ranges and relative epidemiological importance of the distinct anopheline species that constitute the extraordinary diverse and fluctuating malaria vector system in Africa. Such knowledge is of paramount importance for a comprehensive, efficient and sustainable implementation of vector control as a means to alleviate the malaria burden in Africa.</p>
<title>Conclusion</title> <p>In conclusion, this study provides strong support for considering <italic>A. m. bervoetsi </italic>as a full-rank, genetically independent, species within the <italic>A. moucheti </italic>group of malaria vectors. The species should henceforth be named <italic>Anopheles bervoetsi</italic>. However, its epidemiological role as a vector of human malaria parasites still deserves further investigation because roughly nothing is known to date on its biology and behaviour. Isolation by distance seems to be the major factor shaping <italic>A. moucheti s.s</italic>. populations' genetic structure throughout its distribution range across forested areas of Central Africa but significant geographical barriers to gene flow exist, as evidenced from reduced effective population size and high levels of genetic differentiation observed in a population collected from an island on Lake Victoria. Such genetically isolated populations in a geographically confined environment might be of considerable interest for a safe assessment of new and innovative vector control strategies aiming at population suppression and/or replacement, such as those based on the release of sterile or otherwise genetically altered mosquitoes.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p><italic>Anopheles moucheti </italic>is a major malaria vector in forested areas of Africa. However, despite its important epidemiological role, it remains poorly known and insufficiently studied. Here, levels of genetic differentiation were estimated between different <italic>A. moucheti </italic>populations sampled throughout its distribution range in Central Africa.</p> <title>Methods</title> <p>Polymorphism at ten microsatellite markers was compared in mosquitoes sampled in Cameroon, the Democratic Republic of Congo and an island on Lake Victoria in Uganda. Microsatellite data were used to estimate genetic diversity within populations, their relative long-term effective population size, and the level of genetic differentiation between them.</p> <title>Results</title> <p>All specimens collected in Tsakalakuku (Democratic Republic of Congo) were identified as <italic>A. m. bervoetsi </italic>while other samples consisted of <italic>A. m. moucheti</italic>. Successful amplification was obtained at all microsatellite loci within all <italic>A. m. moucheti </italic>samples while only six loci amplified in <italic>A. m. bervoetsi</italic>. Allelic richness and heterozygosity were high for all populations except the island population of Uganda and <italic>A. m. bervoetsi</italic>. High levels of genetic differentiation were recorded between <italic>A. m. bervoetsi </italic>and each <italic>A. m. moucheti </italic>sample as well as between the island population of <italic>A. m. moucheti </italic>and mainland populations. Significant isolation by distance was evidenced between mainland populations.</p> <title>Conclusion</title> <p>High levels of genetic differentiation supports complete speciation of <italic>A. m. bervoetsi </italic>which should henceforth be recognized as a full species and named <italic>A. bervoetsi</italic>. Isolation by distance is the main force driving differentiation between mainland populations of <italic>A. m. moucheti</italic>. Genetically and geographically isolated populations exist on Lake Victoria islands, which might serve as relevant field sites for evaluation of innovative vector control strategies.</p>
<title>Authors' contributions</title> <p>CAN was involved in the study design and implementation, conducted field sampling, microsatellite genotyping and data analysis,  prepared and drafted the manuscript. FS supervised the study conception and design, contributed to data analysis and drafted the manuscript. CN participated to field sampling in Cameroon and microsatellite genotyping. LM organized and greatly contributed to sampling in Uganda. PK was involved in molecular analysis and helped with markers selection resources. PAA and DF were involved in the conception of the study and revised the manuscript.</p>
<title>Acknowledgements</title> <p>We are indebted to E. Mazambi and F. Watsenga of the Institut National de Recherche Biomédicale (INRB), Kinshasa, DRC; F. Senfuka of the Uganda Virus Research Institute (UVRI) and G. Ollivier of the French Foreign Office in DRC for valuable help during mosquito collections. This work was supported by grants no. A20727 and A60347 from the UNICEF/UNDP/World Bank/WHO Special Programme for Research and Training in Tropical Diseases (TDR) to C.A.N and the French Institut de Recherche pour le Développement (IRD/UR016).</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p>A schematic map of Africa showing sampling sites for <italic>A. moucheti </italic>in Cameroon, Democratic Republic of Congo (DRC) and Uganda.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Bayesian cluster analysis using STRUCTURE </bold>[##REF##10835412##22##]. Graphical representation of the data set for the most likely <italic>K </italic>(<italic>K </italic>= 3), where each colour corresponds to a suggested cluster and each individual is represented by a vertical bar. The numbers in the X-axis correspond to a specific sample: 1-Simbock, 2-Olama, 3-Nyabessan, 4-Mouloundou, 5-Kenge, 6-Uganda, 7-Tsakalakuku. The Y-axis represents the probability of assignment of an individual to each cluster.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Correlation between average Fst estimates over 10 microsatellite loci and logarithm of geographic distance between collection sites for pairwise comparisons of seven <italic>A. moucheti </italic>populations from Cameroon, DRC and Uganda.</bold> The name of each sampling site refers to comparisons involving these populations.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Genetic diversity at 10 microsatellite loci in <italic>Anopheles moucheti </italic>from Cameroon<sup>a</sup>, DRC and Uganda.</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td/><td align="center" colspan="4">Cameroon<sup>a</sup></td><td align="center" colspan="2">DRC</td><td align="right">Uganda</td><td/></tr><tr><td/><td/><td colspan="4"><hr/></td><td colspan="2"><hr/></td><td colspan="1"><hr/></td><td/></tr><tr><td align="center">Locus</td><td/><td align="right">Simbock (2n = 118)</td><td align="right">Olama (2n = 112)</td><td align="right">Nyabessan (2n = 108)</td><td align="right">Mouloundou (2n = 108)</td><td align="right">Kenge (2n = 22)</td><td align="right">Tsakalakuku (2n = 128)</td><td align="right">Bufumira (2n = 114)</td><td align="right">All (2n = 710)</td></tr></thead><tbody><tr><td align="center">AM1</td><td align="right"><italic>Rs</italic></td><td align="right">5.99</td><td align="right">6.30</td><td align="right">5.35</td><td align="right">6.92</td><td align="right">5.87</td><td align="right"><italic>NA</italic></td><td align="right">4.41</td><td align="right">6.46</td></tr><tr><td/><td align="right">He</td><td align="right">0.805</td><td align="right">0.823</td><td align="right">0.805</td><td align="right">0.854</td><td align="right">0.853</td><td align="right">-</td><td align="right">0.662</td><td align="right">0.834</td></tr><tr><td/><td align="right">Fis</td><td align="right">+0.139</td><td align="right"><bold>+0.264</bold></td><td align="right"><bold>+0.347</bold></td><td align="right"><bold>+0.338</bold></td><td align="right">+0.372</td><td align="right">-</td><td align="right"><bold>-0.037</bold></td><td align="right"><bold>+0.287</bold></td></tr><tr><td align="center">AM2</td><td align="right"><italic>Rs</italic></td><td align="right">6.06</td><td align="right">6.87</td><td align="right">6.0</td><td align="right">6.67</td><td align="right">7.75</td><td align="right">2</td><td align="right">4.58</td><td align="right">6.75</td></tr><tr><td/><td align="right">He</td><td align="right">0.807</td><td align="right">0.847</td><td align="right">0.812</td><td align="right">0.836</td><td align="right">0.895</td><td align="right">0.503</td><td align="right">0.719</td><td align="right">0.845</td></tr><tr><td/><td align="right">Fis</td><td align="right">+0.119</td><td align="right">+0.186</td><td align="right"><bold>+0.342</bold></td><td align="right"><bold>+0.174</bold></td><td align="right">+0.267</td><td align="right">-1</td><td align="right"><bold>+0.321</bold></td><td align="right"><bold>+0.188</bold></td></tr><tr><td align="center">AM5</td><td align="right"><italic>Rs</italic></td><td align="right">8.06</td><td align="right">7.22</td><td align="right">8.63</td><td align="right">7.37</td><td align="right">5.70</td><td align="right">3.39</td><td align="right">5.76</td><td align="right">8.87</td></tr><tr><td/><td align="right">He</td><td align="right">0.881</td><td align="right">0.862</td><td align="right">0.898</td><td align="right">0.870</td><td align="right">0.723</td><td align="right">0.639</td><td align="right">0.772</td><td align="right">0.906</td></tr><tr><td/><td align="right">Fis</td><td align="right"><bold>+0.150</bold></td><td align="right">+0.018</td><td align="right">+0.069</td><td align="right">+0.168</td><td align="right"><bold>+0.758</bold></td><td align="right">-0.027</td><td align="right"><bold>+0.355</bold></td><td align="right">+0.21</td></tr><tr><td align="center">AM6</td><td align="right"><italic>Rs</italic></td><td align="right">7.24</td><td align="right">6.40</td><td align="right">7.33</td><td align="right">7.13</td><td align="right">3.72</td><td align="right"><italic>NA</italic></td><td align="right">2.69</td><td align="right">7.00</td></tr><tr><td/><td align="right">He</td><td align="right">0.872</td><td align="right">0.826</td><td align="right">0.870</td><td align="right">0.850</td><td align="right">0.671</td><td align="right">-</td><td align="right">0.482</td><td align="right">0.847</td></tr><tr><td/><td align="right">Fis</td><td align="right"><bold>+0.283</bold></td><td align="right"><bold>+0.298</bold></td><td align="right"><bold>+0.263</bold></td><td align="right"><bold>+0.202</bold></td><td align="right">+0.605</td><td align="right">-</td><td align="right"><bold>+0.700</bold></td><td align="right"><bold>+0.382</bold></td></tr><tr><td align="center">AM9</td><td align="right"><italic>Rs</italic></td><td align="right">6.54</td><td align="right">6.58</td><td align="right">7.80</td><td align="right">5.85</td><td align="right">5.38</td><td align="right"><italic>NA</italic></td><td align="right">2.31</td><td align="right">6.63</td></tr><tr><td/><td align="right">He</td><td align="right">0.762</td><td align="right">0.800</td><td align="right">0.875</td><td align="right">0.785</td><td align="right">0.745</td><td align="right">-</td><td align="right">0.201</td><td align="right">0.785</td></tr><tr><td/><td align="right">Fis</td><td align="right">-0.026</td><td align="right">-0.011</td><td align="right"><bold>+0.042</bold></td><td align="right">+0.039</td><td align="right">+0.152</td><td align="right">-</td><td align="right">-0.074</td><td align="right"><bold>+0.137</bold></td></tr><tr><td align="center">AM10</td><td align="right"><italic>Rs</italic></td><td align="right">4.28</td><td align="right">4.14</td><td align="right">5.49</td><td align="right">4.73</td><td align="right">5.00</td><td align="right">1.94</td><td align="right">4.31</td><td align="right">4.85</td></tr><tr><td/><td align="right">He</td><td align="right">0.704</td><td align="right">0.698</td><td align="right">0.790</td><td align="right">0.764</td><td align="right">0.775</td><td align="right">0.260</td><td align="right">0.636</td><td align="right">0.738</td></tr><tr><td/><td align="right">Fis</td><td align="right">-0.065</td><td align="right">+0.028</td><td align="right">+0.110</td><td align="right">+0.019</td><td align="right">+0.205</td><td align="right">+0.107</td><td align="right">-0.159</td><td align="right">+0.139</td></tr><tr><td align="center">AM13</td><td align="right"><italic>Rs</italic></td><td align="right">5.93</td><td align="right">6.12</td><td align="right">5.49</td><td align="right">6.16</td><td align="right">5.37</td><td align="right">1</td><td align="right">4.78</td><td align="right">5.92</td></tr><tr><td/><td align="right">He</td><td align="right">0.823</td><td align="right">0.824</td><td align="right">0.802</td><td align="right">0.830</td><td align="right">0.737</td><td align="right">-</td><td align="right">0.669</td><td align="right">0.848</td></tr><tr><td/><td align="right">Fis</td><td align="right">+0.008</td><td align="right">+0.042</td><td align="right">+0.177</td><td align="right">+0.077</td><td align="right">-0.091</td><td align="right">-</td><td align="right">-0.042</td><td align="right">+0.292</td></tr><tr><td align="center">AM15</td><td align="right"><italic>Rs</italic></td><td align="right">7.10</td><td align="right">6.85</td><td align="right">8.21</td><td align="right">7.17</td><td align="right">6.00</td><td align="right">6.13</td><td align="right">1.60</td><td align="right">7.85</td></tr><tr><td/><td align="right">He</td><td align="right">0.820</td><td align="right">0.822</td><td align="right">0.869</td><td align="right">0.818</td><td align="right">0.834</td><td align="right">0.825</td><td align="right">0.083</td><td align="right">0.844</td></tr><tr><td/><td align="right">Fis</td><td align="right"><bold>+0.188</bold></td><td align="right">+0.094</td><td align="right">+0.049</td><td align="right">+0.119</td><td align="right">+0.263</td><td align="right"><bold>+0.193</bold></td><td align="right">-0.025</td><td align="right"><bold>+0.264</bold></td></tr><tr><td align="center">AM16</td><td align="right"><italic>Rs</italic></td><td align="right">5.98</td><td align="right">6.69</td><td align="right">6.20</td><td align="right">6.37</td><td align="right">6.59</td><td align="right">1.45</td><td align="right">3.20</td><td align="right">6.39</td></tr><tr><td/><td align="right">He</td><td align="right">0.823</td><td align="right">0.843</td><td align="right">0.822</td><td align="right">0.829</td><td align="right">0.868</td><td align="right">0.058</td><td align="right">0.605</td><td align="right">0.836</td></tr><tr><td/><td align="right">Fis</td><td align="right">+0.051</td><td align="right">+0.101</td><td align="right">+0.076</td><td align="right">+0.076</td><td align="right">+0.203</td><td align="right">+0.327</td><td align="right"><bold>+0.192</bold></td><td align="right"><bold>+0.273</bold></td></tr><tr><td align="center">AM20</td><td align="right"><italic>Rs</italic></td><td align="right">7.17</td><td align="right">7.22</td><td align="right">6.06</td><td align="right">7.16</td><td align="right">3.88</td><td align="right"><italic>NA</italic></td><td align="right">3.41</td><td align="right">6.94</td></tr><tr><td/><td align="right">He</td><td align="right">0.844</td><td align="right">0.832</td><td align="right">0.788</td><td align="right">0.848</td><td align="right">0.634</td><td align="right">-</td><td align="right">0.532</td><td align="right">0.839</td></tr><tr><td/><td align="right">Fis</td><td align="right">+0.218</td><td align="right">+0.066</td><td align="right"><bold>+0.203</bold></td><td align="right">+0.083</td><td align="right">+0.130</td><td align="right">-</td><td align="right"><bold>+0.271</bold></td><td align="right"><bold>+0.249</bold></td></tr><tr><td colspan="10"><hr/></td></tr><tr><td align="center">Mean across all loci</td><td align="right"><italic>Rs</italic></td><td align="right">6.44</td><td align="right">6.44</td><td align="right">6.66</td><td align="right">6.55</td><td align="right">5.53</td><td align="right">2.65</td><td align="right">3.71</td><td align="right">6.77</td></tr><tr><td/><td align="right">He</td><td align="right">0.814</td><td align="right">0.818</td><td align="right">0.833</td><td align="right">0.830</td><td align="right">0.771</td><td align="right">0.228</td><td align="right">0.534</td><td align="right">0.689</td></tr><tr><td/><td align="right">Fis</td><td align="right">+0.112</td><td align="right">+0.114</td><td align="right">+0.166</td><td align="right">+0.180</td><td align="right">+0.291</td><td align="right">-0.136</td><td align="right">+0.175</td><td align="right"><italic>NC</italic></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Pairwise Fst estimates between <italic>A. moucheti </italic>populations from Cameroon<sup>a</sup>, DRC and Uganda.</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td/><td align="center" colspan="4">Cameroon<sup>a</sup></td><td align="center" colspan="2">DRC</td></tr><tr><td/><td/><td colspan="4"><hr/></td><td colspan="2"><hr/></td></tr><tr><td/><td/><td align="center">Simbock</td><td align="center">Olama</td><td align="center">Nyabessan</td><td align="center">Mouloundou</td><td align="center">Kenge</td><td align="center">Tsakalakuku</td></tr></thead><tbody><tr><td align="center">Cameroon<sup>a</sup></td><td align="right">Simbock</td><td/><td/><td/><td/><td/><td/></tr><tr><td/><td align="right">Olama</td><td align="center">0.011*</td><td/><td/><td/><td/><td/></tr><tr><td/><td align="right">Nyabessan</td><td align="center">0.014*</td><td align="center">0.009*</td><td/><td/><td/><td/></tr><tr><td/><td align="right">Mouloundou</td><td align="center">0.028*</td><td align="center">0.023*</td><td align="center">0.017*</td><td/><td/><td/></tr><tr><td align="center">DRC</td><td align="right">Kenge</td><td align="center">0.049*</td><td align="center">0.037*</td><td align="center">0.040*</td><td align="center">0.032*</td><td/><td/></tr><tr><td/><td align="right">Tsakalakuku</td><td align="center">0.378*</td><td align="center">0.372*</td><td align="center">0.343*</td><td align="center">0.346*</td><td align="center">0.422*</td><td/></tr><tr><td align="center">Uganda</td><td align="right">Bufumira</td><td align="center">0.167*</td><td align="center">0.172*</td><td align="center">0.187*</td><td align="center">0.168*</td><td align="center">0.223*</td><td align="center">0.448*</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Long-term <italic>Ne </italic>estimates based on genetic diversity (expected heterozygosity) in each collection site, assuming microsatellite loci follow an SMM (see text).</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center" colspan="2">Collection site</td><td align="center">He</td><td align="center">Ne (± SD)</td><td align="center">Neμ</td><td align="center">RR</td></tr></thead><tbody><tr><td align="center">Cameroon</td><td align="right">Simbock</td><td align="center">0.814</td><td align="center">42,396 ± 10,934</td><td align="center">3.488</td><td align="right">7.7:1</td></tr><tr><td/><td align="right">Olama</td><td align="center">0.818</td><td align="center">40,670 ± 10,651</td><td align="center">3.649</td><td align="right">8.1:1</td></tr><tr><td/><td align="right">Nyabessan</td><td align="center">0.833</td><td align="center">53,046 ± 14,269</td><td align="center">4.357</td><td align="right">9.3:1</td></tr><tr><td/><td align="right">Mouloundou</td><td align="center">0.830</td><td align="center">45,606 ± 12,092</td><td align="center">4.20</td><td align="right">9.3:1</td></tr><tr><td align="center">DRC</td><td align="right">Kenge</td><td align="center">0.771</td><td align="center">36,901 ± 24,361</td><td align="center">2.259</td><td align="right">5:1</td></tr><tr><td/><td align="right">Tsakalakuku</td><td align="center">0.228</td><td align="center">5,230 ± 1,602</td><td align="center">0.085</td><td align="right">0.188:1</td></tr><tr><td align="center">Uganda</td><td align="right">Bufumira</td><td align="center">0.534</td><td align="center">7,983 ± 2,316</td><td align="center">0.451</td><td align="right">1</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>Estimates of P-value for the heterozygosity tests for each population of the <italic>A. moucheti </italic>group.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Collection site</td><td/><td/><td align="center" colspan="3">TPM</td><td align="left">SMM</td></tr><tr><td/><td/><td/><td colspan="3"><hr/></td><td/></tr><tr><td/><td/><td/><td align="left">70%<sup>a</sup></td><td align="left">80%<sup>a</sup></td><td align="left">90%<sup>a</sup></td><td/></tr></thead><tbody><tr><td align="left">Cameroon</td><td align="right">Simbock</td><td align="right"><italic>He&gt;Heq</italic></td><td align="left">7</td><td align="left">6</td><td align="left">4</td><td align="left">3</td></tr><tr><td/><td/><td align="right"><italic>Sign test</italic></td><td align="left">0.360</td><td align="left">0.619</td><td align="left">0.192</td><td align="left">0.066</td></tr><tr><td/><td/><td align="right"><italic>Wilcoxon test</italic></td><td align="left">0.322</td><td align="left">0.695</td><td align="left">0.625</td><td align="left">0.024*</td></tr><tr><td/><td align="right">Olama</td><td align="right"><italic>He&gt;Heq</italic></td><td align="left">7</td><td align="left">6</td><td align="left">5</td><td align="left">3</td></tr><tr><td/><td/><td align="right"><italic>Sign test</italic></td><td align="left">0.377</td><td align="left">0.609</td><td align="left">0.396</td><td align="left">0.067</td></tr><tr><td/><td/><td align="right"><italic>Wilcoxon test</italic></td><td align="left">0.275</td><td align="left">1</td><td align="left">0.492</td><td align="left">0.084</td></tr><tr><td/><td align="right">Nyabessan</td><td align="right"><italic>He&gt;Heq</italic></td><td align="left">8</td><td align="left">8</td><td align="left">6</td><td align="left">4</td></tr><tr><td/><td/><td align="right"><italic>Sign test</italic></td><td align="left">0.158</td><td align="left">0.149</td><td align="left">0.608</td><td align="left">0.184</td></tr><tr><td/><td/><td align="right"><italic>Wilcoxon test</italic></td><td align="left">0.010*</td><td align="left">0.016*</td><td align="left">0.625</td><td align="left">0.432</td></tr><tr><td/><td align="right">Mouloundou</td><td align="right"><italic>He&gt;Heq</italic></td><td align="left">7</td><td align="left">5</td><td align="left">4</td><td align="left">2</td></tr><tr><td/><td/><td align="right"><italic>Sign test</italic></td><td align="left">0.362</td><td align="left">0.376</td><td align="left">0.183</td><td align="left">0.015*</td></tr><tr><td/><td/><td align="right"><italic>Wilcoxon test</italic></td><td align="left">0.557</td><td align="left">0.846</td><td align="left">0.275</td><td align="left">0.019*</td></tr><tr><td align="left">DRC</td><td align="right">Kenge</td><td align="right"><italic>He&gt;Heq</italic></td><td align="left">7</td><td align="left">7</td><td align="left">7</td><td align="left">7</td></tr><tr><td/><td/><td align="right"><italic>Sign test</italic></td><td align="left">0.355</td><td align="left">0.395</td><td align="left">0.351</td><td align="left">0.376</td></tr><tr><td/><td/><td align="right"><italic>Wilcoxon test</italic></td><td align="left">0.019*</td><td align="left">0.105</td><td align="left">0.432</td><td align="left">0.557</td></tr><tr><td/><td align="right">Tsakalakuku</td><td align="right"><italic>He&gt;Heq</italic></td><td align="left">4</td><td align="left">4</td><td align="left">4</td><td align="left">3</td></tr><tr><td/><td/><td align="right"><italic>Sign test</italic></td><td align="left">0.212</td><td align="left">0.212</td><td align="left">0.225</td><td align="left">0.575</td></tr><tr><td/><td/><td align="right"><italic>Wilcoxon test</italic></td><td align="left">0.625</td><td align="left">0.625</td><td align="left">0.625</td><td align="left">1</td></tr><tr><td align="left">Uganda</td><td align="right">Bufurima</td><td align="right"><italic>He&gt;Heq</italic></td><td align="left">2</td><td align="left">2</td><td align="left">1</td><td align="left">1</td></tr><tr><td/><td/><td align="right"><italic>Sign test</italic></td><td align="left">0.015*</td><td align="left">0.015*</td><td align="left">0.002**</td><td align="left">0.002**</td></tr><tr><td/><td/><td align="right"><italic>Wilcoxon test</italic></td><td align="left">0.032</td><td align="left">0.010*</td><td align="left">0.003**</td><td align="left">0.002**</td></tr></tbody></table></table-wrap>
<table-wrap-foot><p>All: refers to populations pooled. 2n, number of chromosomes scored; <italic>Rs</italic>, allele richness [##UREF##7##14##]; He, expected heterozygosity under Hardy-Weinberg equilibrium [##REF##17248844##51##]. Fis was calculated according to Weir &amp; Cockerham [##UREF##11##19##] and goodness of fit to Hardy-Weinberg equilibrium was estimated by the exact test available in Genepop 3.2 [##UREF##10##17##]. Bolded values: P &lt; 0.05 after taking into account multiple tests [##UREF##12##31##]. <italic>NA</italic>, No PCR product could be detected; -, irrelevant because no polymorphism was detected;<italic>NC</italic>, not computed. <sup>a </sup>Data adapted from [##REF##17227668##12##].</p></table-wrap-foot> <table-wrap-foot><p>*P &lt; 0.001; DRC, Democratic Republic of Congo; <sup>a </sup>pairwise Fst estimates between Cameroon populations were already published in [##REF##17227668##12##].</p></table-wrap-foot> <table-wrap-foot><p>He, unbiased heterozygosity [##UREF##8##15##]; Ne, mean effective population size calculated across all loci (± standard deviation); μ, mutation rate; RR, relative ratio of Neμ compared to the <italic>A. m. moucheti </italic>population with the smallest effective population size (Bufurima island).</p></table-wrap-foot> <table-wrap-foot><p>TPM, two-phase mutation model with <sup>a </sup>% single step mutation; SMM, stepwise mutation model.</p><p><italic>He&gt;Heq</italic>, number of loci with He&gt;Heq (out of 10 loci tested in each sample, except the <italic>A. m. bervoetsi </italic>sample from Tsakalakuku where only 5 polymorphic loci were considered). * P &lt; 0.05 and ** P &lt; 0.01 (two tails P-values for deviation from MDE) after correction for multiple testing.</p></table-wrap-foot>
<graphic xlink:href="1475-2875-7-120-1"/> <graphic xlink:href="1475-2875-7-120-2"/> <graphic xlink:href="1475-2875-7-120-3"/>
[{"collab": ["WHO/UNICEF/RBM"], "article-title": ["World Malaria Report 2005"], "source": ["WHO/HTM/MAL/20051102"], "year": ["2005"], "publisher-name": ["World Health Organization, Geneva"]}, {"surname": ["Njan Nloga", "Robert", "Toto", "Carnevale"], "given-names": ["A", "V", "JC", "P"], "italic": ["Anopheles moucheti"], "source": ["Bulletin de Liaison et de Documentation de l'OCEAC"], "year": ["1993"], "volume": ["26"], "fpage": ["63"], "lpage": ["67"]}, {"surname": ["Antonio-Nkondjio", "Awono-Ambene", "Toto", "Meunier", "Zebaze-Kemleu", "Nyambam", "Wondji", "Tchuinkam", "Fontenille"], "given-names": ["C", "P", "JC", "JY", "S", "R", "CS", "T", "D"], "article-title": ["High malaria transmission intensity in a village close to Yaound\u00e9, the capital city of Cameroon"], "source": ["J Med Entomol"], "year": ["2002"], "volume": ["39"], "fpage": ["330"], "lpage": ["355"]}, {"surname": ["Gillies", "De Meillon"], "given-names": ["MT", "B"], "source": ["The Anophelinae of Africa South of the Sahara (Ethiopian zoogeographical region)"], "year": ["1968"], "edition": ["2"], "publisher-name": ["Publications of the South African Institute for Medical Research, Johannesburg"]}, {"surname": ["Brunhes", "Le Goff", "Manga", "Geoffroy"], "given-names": ["J", "G", "L", "B"], "article-title": ["Anoph\u00e8les afro-tropicaux. IV \u2013 Mise au point sur les esp\u00e8ces et sous-esp\u00e8ces du groupe "], "italic": ["Anopheles ", "Cellia", "moucheti", "d'An", "multicinctus ", "d'An", "Cellia", "garnhami "], "source": ["Annales de la Soci\u00e9t\u00e9 Entomologique de France"], "year": ["1998"], "volume": ["34"], "fpage": ["397"], "lpage": ["405"]}, {"surname": ["Annan", "Kengne", "Berthomieu", "Antonio-Nkondjio", "Rousset", "Fontenille", "Weill"], "given-names": ["Z", "P", "A", "C", "F", "D", "M"], "article-title": ["Isolation and characterisation of polymorphic microsatellite markers from the mosquito "], "italic": ["Anopheles moucheti"], "source": ["Mol Ecol Notes"], "year": ["2003"], "volume": ["3"], "fpage": ["56"], "lpage": ["57"], "pub-id": ["10.1046/j.1471-8286.2003.00347.x"]}, {"surname": ["Gillies", "Coetzee"], "given-names": ["MT", "M"], "source": ["A supplement to the Anophelinae of Africa south of the Sahara"], "year": ["1987"], "publisher-name": ["Publications of the South African Institute for Medical Research, Johannesburg"]}, {"surname": ["El Mousadik", "Petit"], "given-names": ["A", "R"], "article-title": ["High level of genetic differentiation for allelic richness among populations of the argan tree ("], "italic": ["Argania spinosa "], "source": ["Theoret Appl Genet"], "year": ["1996"], "volume": ["92"], "fpage": ["832"], "lpage": ["839"], "pub-id": ["10.1007/BF00221895"]}, {"surname": ["Nei"], "given-names": ["M"], "source": ["Molecular evolutionary genetics"], "year": ["1987"], "publisher-name": ["Colombia University press, New York"]}, {"surname": ["Goudet"], "given-names": ["J"], "article-title": ["FSTAT version 2.9.3.2. A computer software to calculate F-statistics"], "source": ["J Hered"], "year": ["1995"], "volume": ["86"], "fpage": ["485"], "lpage": ["486"]}, {"surname": ["Raymond", "Rousset"], "given-names": ["M", "F"], "article-title": ["GENEPOP, Version 1.2. A population genetics software for exact tests and ecumenicism"], "source": ["J Hered"], "year": ["1995"], "volume": ["86"], "fpage": ["248"], "lpage": ["249"]}, {"surname": ["Weir", "Cockerham"], "given-names": ["BS", "CC"], "article-title": ["Estimating F-statistics for the analysis of population structure"], "source": ["Evolution"], "year": ["1984"], "volume": ["38"], "fpage": ["1358"], "lpage": ["1370"], "pub-id": ["10.2307/2408641"]}, {"surname": ["Holm"], "given-names": ["S"], "article-title": ["A simple sequentially rejective multiple test procedure"], "source": ["Scand J Stat"], "year": ["1979"], "volume": ["6"], "fpage": ["65"], "lpage": ["70"]}, {"surname": ["Slatkin", "Real LA"], "given-names": ["M"], "article-title": ["Gene flow and population structure"], "source": ["Ecological genetics"], "year": ["1994"], "fpage": ["3"], "lpage": ["17"]}, {}, {"surname": ["Cohen"], "given-names": ["B"], "article-title": ["Urbanization in developing countries: Current trends, future projections, and key challenges for sustainability"], "source": ["Technology in Society"], "year": ["2006"], "volume": ["28"], "fpage": ["63"], "lpage": ["80"], "pub-id": ["10.1016/j.techsoc.2005.10.005"]}]
Malar J. 2008 Jul 4; 7:120