# Cryptic Genetic Diversity within the Anopheles nili group of Malaria Vectors in the Equatorial Forest Area of Cameroon (Central Africa)

Cyrille Ndo, Frédéric Simard, Pierre Kengne, Parfait H Awono-Ambéné,
Isabelle Morlais, Igor Sharakhov, Didier Fontenille, Christophe
Antonio-Nkondjio

► To cite this version:

Cyrille Ndo, Frédéric Simard, Pierre Kengne, Parfait H Awono-Ambéné, Isabelle Morlais, et al.. Cryptic Genetic Diversity within the Anopheles nili group of Malaria Vectors in the Equatorial Forest Area of Cameroon (Central Africa). PLoS ONE, 2013, 8 (3), pp.e58862. 10.1371/journal.pone.0058862. ird-02897002

HAL Id: ird-02897002

https://ird.hal.science/ird-02897002

Submitted on 25 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Cryptic Genetic Diversity within the *Anopheles nili* group of Malaria Vectors in the Equatorial Forest Area of Cameroon (Central Africa)

Cyrille Ndo¹,²,³*, Frédéric Simard², Pierre Kengne¹,², Parfait Awono-Ambene¹, Isabelle Morlais¹,², Igor Sharakhov⁴, Didier Fontenille², Christophe Antonio-Nkondjio¹,⁵

¹Laboratoire de Recherche sur le Paludisme, Organisation de Coordination pour la lutte Contre les Endémies en Afrique Centrale, Yaoundé, Cameroon, ²Unité de Recherche Mixte Maladies Infectieuses et Vecteurs : Ecologie, Génétique, Evolution et Contrôle, Institut de Recherche pour le Développement, Montpellier, France, ³Faculty of Medicine and Pharmaceutical Sciences, University of Douala, Douala, Cameroon, ⁴Department of Entomology, Virginia Tech, Blacksburg, Virginia, United States of America, ⁵Faculty of Health Sciences, University of Bamenda, Bambili, Cameroon

## Abstract

**Background:** The *Anopheles nili* group of mosquitoes includes important vectors of human malaria in equatorial forest and humid savannah regions of sub-Saharan Africa. However, it remains largely understudied, and data on its populations' bionomics and genetic structure are crucially lacking. Here, we used a combination of nuclear (i.e. microsatellite and ribosomal DNA) and mitochondrial DNA markers to explore and compare the level of genetic polymorphism and divergence among populations and species of the group in the savannah and forested areas of Cameroon, Central Africa.

**Principal Findings:** All the markers provided support for the current classification within the *An. nili* group. However, they revealed high genetic heterogeneity within *An. nili s.s.* in deep equatorial forest environment. Nuclear markers showed the species to be composed of five highly divergent genetic lineages that differed by 1.8 to 12.9% of their Internal Transcribed Spacer 2 (ITS2) sequences, implying approximate divergence time of 0.82 to 5.86 million years. However, mitochondrial data only detected three major subdivisions, suggesting different evolutionary histories of the markers.

**Conclusions/Significance:** This study enlightened additional cryptic genetic diversity within *An. nili s.s.* in the deep equatorial forest environment of South Cameroon, reflecting a complex demographic history for this major vector of malaria in this environment. These preliminary results should be complemented by further studies which will shed light on the distribution, epidemiological importance and evolutionary history of this species group in the African rainforest, providing opportunities for in-depth comparative studies of local adaptation and speciation in major African malaria vectors.

**Citation:** Ndo C, Simard F, Kengne P, Awono-Ambene P, Morlais I, et al. (2013) Cryptic Genetic Diversity within the *Anopheles nili* group of Malaria Vectors in the Equatorial Forest Area of Cameroon (Central Africa). PLoS ONE 8(3): e58862. doi:10.1371/journal.pone.0058862

**Editor:** Roman Ganta, Kansas State University, United States of America

**Received** October 31, 2012; **Accepted** February 7, 2013; **Published** March 14, 2013

**Copyright:** © 2013 Ndo et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This study was funded by a UNICEF/United Nations Development Programme/World Bank/World Health Organization/TDR/Molecular Entomology grant A60347 to CAN and a National Institutes of Health grant R21 AI079350 to IS. CN was supported by a student fellowship grant from the French Institut de Recherche pour le Développement (IRD). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: cyrndo@yahoo.fr

## Introduction

In the equatorial forest regions of Africa, malaria remains endemic despite significant efforts in treatment of the disease and vector control [1]. Part of explanation for this situation resides in the fact that a number of genetically, behaviorally and ecologically distinct species and populations of *Anopheles* mosquitoes can and do transmit *Plasmodium*, the malignant agent of the disease, simultaneously or replace each other seasonally insuring year-round transmission [2]. Unfortunately, many of these vectors might escape traditional vector control operations by Indoor Residual Spraying of insecticides (IRS) and Insecticide-Treated bed Nets (ITNs), which are directed against indoor biting and/or resting mosquitoes [3,4]. This is exemplified in the *Anopheles nili* group of malaria vectors.

The *An. nili* group includes four closely related species namely *An. nili sensu stricto* (thereafter, *An. nili s.s.*), *An. somalicus*, *An. carnevalei* and *An. ovengensis* which can be identified through slight morphological diagnostic characters present at the larval and/or adult stages [5–7], and by a species-specific PCR assay [8]. Their larvae typically develop in river networks with a marked preference for the edges of fast running streams and rivers exposed to light and containing vegetation and/or floating debris [9]. *Anopheles nili s.s.* is widespread and is found throughout sub-Saharan Africa. It is highly anthropophilic, biting man and resting indoors as well as outdoors [5,10–16]. The other members of the group are more localized within the equatorial forest region of Central Africa, although their true distribution range remains largely unknown [9,17,18]. *Anopheles somalicus* is not involved in malaria transmission because of its preference for non-human blood and exophilic habits [2,19]. A comprehensive study in
Cameroon confirmed that *An. nili* s.s. is the major malaria vector of the group and emphasized the exophagic behavior of *An. ovengensis* and *An. carnevalei* [20]. In Equatorial Guinea, sporozoite rates in *An. ovengensis* can reach 4.1% (n = 74), which is higher than that of *An. gambiae* in the same area (3.3%, n = 603, [17]), confirming its major role in the epidemiology of malaria. Moreover, recent studies found the circulation of various *Plasmodium* species and lineages related to *P. falciparum* in great apes and infection of gorillas and guenons by *P. falciparum* in the central forest region of Cameroon and Gabon [21–25]. However, the potential mosquito vectors have not been identified yet. These discoveries further bring urgency to the study of the taxonomic status of members of the *An. nili* group and their possible role in transmission of various *Plasmodium* species.

A recent study using 11 microsatellite markers demonstrated a significant genetic differentiation of the forest *An. nili* s.s. population of Kenge in the Democratic Republic of Congo (DRC) when compared to other populations from humid savannah of Central and West Africa [18]. Both local adaptation and geographic isolation could cause this differentiation. Extensive allele sharing between populations and homogeneity across loci suggested that enhanced genetic drift rather than selection was responsible for the observed pattern. Fluorescence *in situ* hybridization (FISH) mapping of the microsatellite loci on a recently developed cytogenetic map for *An. nili* [26,27] demonstrated physical independence of the markers, strengthening the view that enhanced genetic drift, rather than selection, was responsible for reduced variability and increased differentiation of the Kenge population. These data strongly suggest existence of ecological and geographical barriers that affect gene flow among *An. nili* natural populations in equatorial Africa.

Here, we used nuclear (i.e. microsatellites and rDNA) and mitochondrial DNA markers to further explore the genetic structure of members of the *An. nili* group in South Cameroon, where all four members of the group are known to occur [7,8,20]. We document further cryptic genetic substructure within *An. nili* s.s., reflecting a complex demographic history. We advocate for further evaluation of the impact of this genetic heterogeneity and its relevance for malaria epidemiology and control.

## Materials and Methods

### Ethics Statement

This study received the approval of the National Ethics Committee of Cameroon (Ref: FWA IRB0001954). No other specific permits were required. No location was privately-owned or protected and field studies did not involve endangered or protected species. Consent was obtained from all volunteers participating in the study, and they were monitored for signs of fever and treated for all diagnosed malaria cases.

### Mosquito Collection and Identification

Mosquitoes were collected from 2007 to 2009 in eight riverside villages situated in different ecological settings in Cameroon including Afan-Essokyé, Ako, Ekelemba, Kentzou, Moloundou, Mbébé, Nkolbisson, and Nyabessan (Figure 1, Table 1). Field-collected larvae and adult *An. nili* s.l. mosquitoes were first sorted using morphological diagnostic characters [5–7,28]. Adult females were stored individually in numbered tubes containing a desiccant whereas larvae were preserved in tubes containing 95% ethanol. Then, genomic DNA was extracted using a standard protocol [29] and morphological identification was confirmed by PCR [8].

## Microsatellite Genotyping and Analysis

Eleven microsatellite loci developed in *An. nili* s.s. [27,30] were used to amplify *An. carnevalei*, *An. ovengensis* and *An. somalicus* DNA. Genotypes were determined for individual mosquitoes as previously described [18]. Genetic diversity at each microsatellite locus within each sample was characterized by estimates of allelic richness (Rs) [31] and unbiased expected heterozygosity (He) under Hardy-Weinberg Equilibrium (HWE), using FSAT V2.9.3.2 [32]. This program was also used to test for conformance of genotype frequencies to HWE for each locus and each sample, and to explore Linkage Disequilibrium (LD) between microsatellite loci. The frequency of null alleles at each locus within each population was determined using GENEPOP V4.0 [33,34], and the allele and genotype frequencies were then adjusted accordingly in MICROCHECKER V2.2.3 [35]. The null allele adjusted dataset was compared to the original dataset to explore the effect of null alleles on estimates of genetic differentiation.

Heterozygosity tests were implemented in BOTTLENECK V1.0.02 [36] in order to assess deviation from Mutation-Drift Equilibrium (MDE). The tests were performed using the Stepwise Mutation Model (SMM) and the Two Phase Model (TPM), with 10–30% indels larger than one repeat unit. Statistical significance of the tests was assessed for each sample across all loci by the Wilcoxon signed-rank test available in BOTTLENECK.

Bayesian clustering analysis was implemented in STRUCTURE V2.2 [37] to estimate the most probable number of discrete populations or genetic clusters (K) in the dataset, while minimizing Hardy-Weinberg and gametic phase disequilibrium between loci within clusters, with no a priori assumptions of population structure. The admixture model, which assumes that a fraction of the genome of each individual can come from one of each population, was used. Allele frequencies were allowed to be correlated assuming that the populations diverged from a single ancestral population at some point in the past, and the differences in their allele frequencies are the result of drift that has occurred since their divergence [37]. The method of Evanno *et al.* [38] was used to determine the most likely number of clusters. Differentiation among clusters was examined with Wright F-statistics [39] computed according to Weir and Cockerham [40] using FSTAT. For the calculation of *Fst*, monomorphic loci were excluded from the analysis and significance tests were conducted using a randomization approach as implemented in FSTAT.

The number of putative migrants (*Nm*) per generation between populations was estimated from pairwise *Fst* [41]. The contribution of geographical distance to the pattern of differentiation observed was investigated by the regression of *Fst*/(1-*Fst*) on geographic distance, assuming migration in one direction (i.e. along hydrographic networks) [42]. The significance of this regression was tested using a mantel test available in GENEPOP [43]. The Bonferroni correction [44] was used to adjust the nominal significance level whenever multiple tests were performed.

## DNA Sequencing and Phylogenetic Analysis

Sequencing was performed using the primers set and protocol described in Ndo *et al.* [18]. Ten specimens per geographical sample (but twenty from Moloundou) among those used for microsatellite analysis were sequenced for the Internal Transcribed Spacer 2 (ITS2) and domain-3 (D3) of nuclear 28S ribosomal DNA, and for *Cytochrome Oxydase subunit II (COII)* and *NADH dehydrogenase subunit IV (ND4)* of mitochondrial DNA.

Basic sequence polymorphism statistics and various population parameters were computed using DnaSP V5.10.01 [45] and MEGA V5.05 [46]. Genetic diversity was measured by assessing haplotype diversity (hd) and nucleotide diversity (π). MEGA was
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Map of Cameroon showing collection sites and the 8 microsatellites clusters detected by Bayesian clustering analysis using STRUCTURE. Each colored rectangle represents a distinct genetic cluster within which each vertical bar corresponds to a single individual. Note that Ako and Nkolbisson specimens belong to the same genetic cluster (see text).
doi:10.1371/journal.pone.0058862.g001</figcaption></figure>
<table><caption>**Table 1.** Characteristics of sampling sites and species collection methods and dates.</caption><thead><tr><th>Sites</th><th>Geographical coordinates</th><th>Hydrographic network</th><th>Ecological domain<sup>b</sup></th><th>Species<sup>a</sup></th><th>Stage</th><th>Method of Collection<sup>c</sup></th><th>Collection date</th></tr></thead><tbody><tr><td>Ako</td><td>6°49'N, 10°43'E</td><td>Atoro and Akon</td><td>HS</td><td>An. nili s.s.</td><td>Adult</td><td>HLC</td><td>03/2007</td></tr><tr><td>Nkolbisson</td><td>3°52'N; 11°25'E</td><td>Mefou</td><td>DF</td><td>An. nili s.s.</td><td>Adult</td><td>PSC</td><td>04/2008</td></tr><tr><td>Kentzou</td><td>4°09'N, 15°02'E</td><td>Kadei</td><td>DF</td><td>An. nili s.s.</td><td>Larva</td><td>LC</td><td>02/2009</td></tr><tr><td>Moloundou</td><td>2°02'N, 15°12'E</td><td>Boumba and Ngoko</td><td>DEF</td><td>An. nili s.s.</td><td>Larva+Adult</td><td>HLC/LC</td><td>10/2008 & 02/2009</td></tr><tr><td>Ekelemba</td><td>1°49'N, 15°54'E</td><td>Baad</td><td>DEF</td><td>An. nili s.s.</td><td>Larva</td><td>LC</td><td>02/2009</td></tr><tr><td>Mbébé</td><td>4°10'N; 11°04'E</td><td>Sanaga</td><td>DF</td><td>An. somalicus</td><td>Larva</td><td>LC</td><td>03/2009</td></tr><tr><td>Nyabessan</td><td>2°10'N; 10°40'E</td><td>Ntem and Njoh</td><td>DEF</td><td>An. ovengensis</td><td>Adult</td><td>HLC</td><td>12/2007</td></tr><tr><td>Afan-Essokyé</td><td>2°22'N; 9°59'E</td><td>Bitandé</td><td>DEF</td><td>An. carnevalei</td><td>Adult</td><td>HLC</td><td>12/2007</td></tr></tbody></table>

also used to compute genetic distances between haplotypes and/or populations [47]. The rate of the ITS2 clock of 2.2% per million years (myr), as assessed within the *Drosophila melanogaster* and *D. virilis* groups [48], was used to estimate divergence time between sequences. For mtDNA sequences, the conventional mtDNA rate of 2% per myr was used [49].

The Neighbour-Joining (NJ) method [50] performed in MEGA was used for the construction of phylogenetic trees, which were rooted using *An. moucheti* sequences ([63]; Ndo, unpublished) as an outgroup.

## Results

### Samples Composition

All mosquitoes from Ako, Ekelemba, Kentzou, Moloundou and Nkolbisson were morphologically and molecularly identified as *An. nili s.s.*, although some specimens from Moloundou and Ekelemba failed to amplify even after repeated attempts, possibly due to poor DNA quality. All individuals from Nyabessan, Afan-Essokyé and Mbébé were *An. ovengensis*, *An. carnevalei* and *An. somalicus* respectively, but *An. nili s.s.* is also known to be present in the two later localities [9]. We never collected *An. somalicus* on human baits owing to its zoophilic behavior. Similarly, attempts to collect *An. nili s.s.* on human baits in Kentzou were unsuccessful suggesting possible zoophilic behavior of the species in this locality.

### Microsatellite Analysis

**Genetic variability.** A total of 363 individual mosquitoes were genotyped and there were marked differences in amplification profiles and polymorphism of loci among species, as well as between *An. nili s.s.* populations (Table 2). All 11 loci readily amplified and were polymorphic in *An. nili s.s.* samples from Ako and Nkolbisson, whereas only 10 (90.9%) successfully amplified in Kentzou (9 polymorphic) and Moloundou (9 polymorphic), falling down to 8 (73%) in Ekelemba (with only 5 loci being polymorphic). Similarly, only 9 (82%) loci successfully amplified *An. carnevalei* and *An. ovengensis* DNA, among which 7 and 6 respectively were polymorphic. In *An. somalicus*, 8 (73%) loci were scored successfully, with only 5 showing more than one allele (Table 2). These differences in locus amplification and polymorphism resulted in high genetic variability within the *An. nili s.s.* populations from Ako and Nkolbisson (Rs: 8.33–8.93, He: 0.833–0.857), whereas polymorphism was much lower in the other samples (Rs: 1.77–6.46, He: 0.231–0.659).

### Hardy-Weinberg Equilibrium and Linkage Disequilibrium

Only polymorphic loci were considered for HWE and LD analyses. A total of 21 (33.33%) tests out of 63 were found out of HWE (P<0.001) after Bonferroni correction for multiple testing (Table 2). Within *An. nili s.s.*, the sample from Moloundou had the highest number of deviations. In fact, all 9 polymorphic loci were out of HWE, associated with positive *Fis* values suggesting population substructure, although null alleles could also account for some deviation. By contrast, no more than 2 significant deviations were observed in all other *An. nili s.s.* samples. The number of significant deviations from HWE detected in the other members of the group varied from 1 (out of 5 polymorphic loci) in *An. somalicus* to 3 (out of 7 polymorphic loci) in *An. carnevalei* (Table 2). The frequency of null alleles at those loci showing significant deviation from HWE ranged from 8.73 to 41.65% (Table S1), suggesting mutations in the primer-annealing sequences were probably involved. Nevertheless, null alleles did not significantly bias our subsequent interpretation of populations' differentiation.

Exact tests for LD within each *An. nili s.s.* population resulted in 5 (2.6%) significant values (P<0.05) out of 192 comparisons after correction by the Bonferroni procedure. All these deviations were detected in Moloundou where all polymorphic loci were found out of HWE, and involved loci 1D80 and F41, 2Ateta and F41, B115 and F41, 2Ateta and 1A27, and B115 and F56. Because all these loci are located on different chromosomal arms [27], this result further supports HWE analysis in pointing towards splitting of the gene pool in the Moloundou population (i.e. the Wahlund's effect). No significant LD between loci was detected in *An. carnevalei*, *An. ovengensis* and *An. somalicus* populations.

### Genetic Differentiation and Gene Flow

The Bayesian clustering analysis using the software STRUCTURE was implemented in the expectation that the genotypes would segregate into four main clusters corresponding to the four previously recognized species within the *An. nili* group. The optimal number of clusters detected based on the method of Evano *et al.* [38] was K=8 (Figure 1). *Anopheles carnevalei*, *An. somalicus* and *An. ovengensis* samples each formed a distinct and homogeneous genetic cluster, while five clusters were defined within *An. nili s.s.* Specimens sampled in savannah (Ako) and degraded forest (Nkolbisson) areas were merged together into a single genetic cluster, whereas four clusters were detected among populations sampled within the evergreen forest block: one in

<footnote><sup>a</sup>morphological and molecular identification;</footnote>

<footnote><sup>b</sup>HS, Humid Savannah; DF, Degraded forest; DEF, Deep evergreen forest;</footnote>

<footnote><sup>c</sup>HLC, Human Landing Collection; PSC, Pyrethrum Spray Catch; LC, Larval Collection through the dipping technique [67].</footnote>

<footnote>doi:10.1371/journal.pone.0058862.t001</footnote>
<table><caption>**Table 2.** Genetic variability at eleven microsatellite loci in 8 *An. nili* s.l. populations from Cameroon.</caption><thead><tr><th rowspan="2">Locus</th><th colspan="5">An. nili s.s.</th><th>An. carnevalei</th><th>An. ovengensis</th><th>An. somalicus</th></tr><tr><th>Ako</th><th>Nkolbisson</th><th>Kentzou</th><th>Moloundou</th><th>Ekelemba</th><th>Afan- Essokyé</th><th>Nyabessan</th><th>Mbébé</th></tr></thead><tbody><tr><td>1D80</td><td>41</td><td>38</td><td>53</td><td>71</td><td>32</td><td>35</td><td>48</td><td>46</td></tr><tr><td>Rs</td><td>10.81</td><td>8.27</td><td>5.99</td><td>10.50</td><td>3.50</td><td>6.04</td><td>7.82</td><td>5.11</td></tr><tr><td>He</td><td>0.917</td><td>0.869</td><td>0.713</td><td>0.913</td><td>0.310</td><td>0.775</td><td>0.829</td><td>0.765</td></tr><tr><td>Fis</td><td>0.151**</td><td>0.074</td><td>0.093</td><td><strong>0.232***</strong></td><td>-0.115</td><td>0.052</td><td><strong>0.731***</strong></td><td>0.148</td></tr><tr><td>1A27</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rs</td><td>13.03</td><td>11.75</td><td>11.63</td><td>12.61</td><td>6.93</td><td>-</td><td>5.07</td><td>-</td></tr><tr><td>He</td><td>0.946</td><td>0.937</td><td>0.927</td><td>0.948</td><td>0.739</td><td>-</td><td>0.571</td><td>-</td></tr><tr><td>Fis</td><td>0.103*</td><td><strong>0.269***</strong></td><td>0.091*</td><td><strong>0.239***</strong></td><td>-0.358***</td><td>-</td><td>0.203*</td><td>-</td></tr><tr><td>2Ateta</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rs</td><td>9.22</td><td>8.92</td><td>7.54</td><td>5.11</td><td>1.31</td><td>7.36</td><td>6.73</td><td>3.82</td></tr><tr><td>He</td><td>0.861</td><td>0.886</td><td>0.849</td><td>0.696</td><td>0.030</td><td>0.817</td><td>0.784</td><td>0.610</td></tr><tr><td>Fis</td><td>0.011</td><td>0.112</td><td>0.095</td><td><strong>0.585***</strong></td><td>0.000</td><td>0.044</td><td>0.036</td><td>0.144</td></tr><tr><td>A14</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rs</td><td>8.60</td><td>6.06</td><td>1.19</td><td>1.30</td><td>1.00</td><td>1.00</td><td>1.00</td><td>5.83</td></tr><tr><td>He</td><td>0.858</td><td>0.802</td><td>0.019</td><td>0.083</td><td>0.000</td><td>0.000</td><td>0.000</td><td>0.787</td></tr><tr><td>Fis</td><td>0.112</td><td><strong>0.255**</strong></td><td>0.000</td><td><strong>0.830***</strong></td><td>-</td><td>-</td><td>-</td><td><strong>0.501***</strong></td></tr><tr><td>A154</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rs</td><td>9.31</td><td>7.17</td><td>1.45</td><td>8.65</td><td>3.00</td><td>6.11</td><td>6.16</td><td>-</td></tr><tr><td>He</td><td>0.825</td><td>0.680</td><td>0.051</td><td>0.887</td><td>0.660</td><td>0.793</td><td>0.769</td><td>-</td></tr><tr><td>Fis</td><td>0.114*</td><td>0.002</td><td>1.000</td><td><strong>0.821***</strong></td><td>-0.086</td><td><strong>0.678***</strong></td><td><strong>0.501***</strong></td><td>-</td></tr><tr><td>2C157</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rs</td><td>5.92</td><td>4.84</td><td>1.00</td><td>1.00</td><td>1.00</td><td>2.00</td><td>1.00</td><td>1.00</td></tr><tr><td>He</td><td>0.766</td><td>0.670</td><td>0.000</td><td>0.000</td><td>0.000</td><td>0.507</td><td>0.000</td><td>0.000</td></tr><tr><td>Fis</td><td>-0.018</td><td>0.063</td><td>-</td><td>-</td><td>-</td><td><strong>1.000***</strong></td><td>-</td><td>-</td></tr><tr><td>F56</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rs</td><td>10.34</td><td>9.30</td><td>4.82</td><td>8.45</td><td>-</td><td>4.99</td><td>5.57</td><td>3.29</td></tr><tr><td>He</td><td>0.907</td><td>0.891</td><td>0.665</td><td>0.901</td><td>-</td><td>0.844</td><td>0.617</td><td>0.830</td></tr><tr><td>Fis</td><td><strong>0.360***</strong></td><td>0.089</td><td>0.250**</td><td><strong>0.232***</strong></td><td>-</td><td><strong>0.556***</strong></td><td>0.027</td><td>0.199**</td></tr><tr><td>B115</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rs</td><td>6.20</td><td>6.22</td><td>4.58</td><td>7.19</td><td>1.62</td><td>1.00</td><td>6.99</td><td>1.00</td></tr><tr><td>He</td><td>0.802</td><td>0.814</td><td>0.677</td><td>0.791</td><td>0.088</td><td>0.000</td><td>0.838</td><td>0.000</td></tr><tr><td>Fis</td><td>-0.143</td><td>-0.027</td><td>-0.117</td><td><strong>0.499***</strong></td><td>-0.033</td><td>-</td><td>0.036</td><td>-</td></tr><tr><td>F41</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rs</td><td>11.76</td><td>14.43</td><td>5.45</td><td>8.93</td><td>1.00</td><td>6.63</td><td>1.00</td><td>11.34</td></tr><tr><td>He</td><td>0.935</td><td>0.960</td><td>0.733</td><td>0.864</td><td>0.000</td><td>0.785</td><td>0.000</td><td>0.919</td></tr><tr><td>Fis</td><td>0.080</td><td>-0.005</td><td>-0.159*</td><td><strong>0.323***</strong></td><td>-</td><td>0.187</td><td>-</td><td>0.064</td></tr><tr><td>1F43</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rs</td><td>6.57</td><td>6.41</td><td>4.53</td><td>7.32</td><td>-</td><td>-</td><td>-</td><td>1.00</td></tr><tr><td>He</td><td>0.817</td><td>0.804</td><td>0.400</td><td>0.508</td><td>-</td><td>-</td><td>-</td><td>0.000</td></tr><tr><td>Fis</td><td><strong>0.320***</strong></td><td>-0.048</td><td>0.116</td><td><strong>0.641***</strong></td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1G13</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rs</td><td>6.51</td><td>8.30</td><td>-</td><td>-</td><td>-</td><td>6.84</td><td>-</td><td>-</td></tr><tr><td>He</td><td>0.783</td><td>0.851</td><td>-</td><td>-</td><td>-</td><td>0.692</td><td>-</td><td>-</td></tr><tr><td>Fis</td><td>0.009</td><td><strong>0.254***</strong></td><td>-</td><td>-</td><td>-</td><td>0.281**</td><td>-</td><td>-</td></tr><tr><td>Mean over loci</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rs</td><td>8.93</td><td>8.33</td><td>4.38</td><td>6.46</td><td>1.77</td><td>3.81</td><td>3.76</td><td>3.29</td></tr><tr><td>He</td><td>0.857</td><td>0.833</td><td>0.503</td><td>0.659</td><td>0.231</td><td>0.579</td><td>0.490</td><td>0.489</td></tr></tbody></table>
<table><caption>**Table 2.** Cont.</caption><thead><tr><th rowspan="2">Locus</th><th colspan="5">An. nili s.s.</th><th>An. carnevalei</th><th>An. ovengensis</th><th>An. somalicus</th></tr><tr><th>Ako</th><th>Nkolbisson</th><th>Kentzou</th><th>Moloundou</th><th>Ekelemba</th><th>Afan- Essokyé</th><th>Nyabessan</th><th>Mbébé</th></tr></thead><tbody><tr><td>Fis</td><td>0.104</td><td><strong>0.098***</strong></td><td>0.060*</td><td><strong>0.437***</strong></td><td>0.060</td><td><strong>0.374***</strong></td><td><strong>0.269***</strong></td><td><strong>0.210***</strong></td></tr></tbody></table>

Kentzou, one in Ekelemba and two in Moloundou, namely 'Moloundou A' and 'Moloundou B' in subsequent analyses.

High levels of genetic differentiation were recorded between species. The overall *Fst* estimate across all loci was 0.373 and was highly significant (P<0.0001). Mean pairwise *Fst* over all loci varied from 0.201 to 0.659 (P<0.0001) with corresponding *Nn* values in the range of 0.14 to 1.16, suggesting very restricted gene flow between species (Table 3).

Within *An. nili s.s.*, high levels of genetic differentiation (*Fst*: 0.16–0.45; P<0.0001) and correspondingly low migration index (*Nn*: 0.28–1.07) were estimated between populations, except between Ako and Nkolbisson which belong to the same genetic cluster (Table 3).

No significant correlation was detected between genetic and geographic distances (R² = 0.0178; P = 0.304) and there was no indication of recent population contraction in any of the populations tested (data not shown).

## DNA Sequencing Analysis

**rDNA polymorphism and divergence.** The D3 domain of the 28S rDNA was amplified in 90 specimens (10 from each geographic population, and 20 from Moloundou) and the alignment length, including insertions/deletions (indels), was 382 bp (Figure S1). No sequence polymorphism (π = 0) was detected within each geographical population, except in Moloundou where two highly divergent sequence variants were observed (π = 0.034), corresponding to the Moloundou A and Moloundou B microsatellites clusters. Haplotypes detected in *An. carnevalei*, *An. ovengensis* and *An. somalicus* populations, and in *An. nili s.s.* from Ako

and Nkolbisson were identical to those published by Kengne *et al.* [8] and Ndo *et al.* [18] (EMBL accession numbers: AJ429053 for *An. nili s.s.*, AJ429052 for *An. ovengensis*, AJ438690 for *An. somalicus* and AJ429051 for *An. carnevalei*). The four remaining *An. nili s.s.* haplotypes were new and found exclusively in deep forest populations (GenBank accession numbers: KC189970–KC189973). They differed from one another, and from the deposited *An. nili s.s.* haplotype by 13 to 16 fixed nucleotide substitutions and many indels. Genetic distances between species' haplotypes ranged from 0.3 to 5%. Similar genetic distances were also estimated between *An. nili s.s.* haplotypes, ranging from 0.5 to 5.9% (Table S2).

Likewise, there was no intra-population variation for the ITS2 region and the alignment length was 480 bp for *An. carnevalei*, 503 bp for *An. ovengensis*, 513 bp for *An. somalicus*, 451 bp for *An. nili s.s.* from Ako and Nkolbisson, 445 bp for *An. nili s.s.* from Kentzou and Moloundou A and 480 bp for sequences from Ekelemba and Moloundou B (Figure S2). ITS2 sequences also segregated into eight haplotypes, in straight correspondence with the four microsatellite clusters and the D3 sequence analysis. Four of these haplotypes were identical to those described earlier (EMBL accession numbers: AJ429048 for *An. nili s.s.*, AJ429049 for *An. ovengensis*, AJ438689 for *An. somalicus* and AJ429050 for *An. carnevalei*) [8,18]. The four new haplotypes were also observed in the deep forest *An. nili s.s.* populations (GenBank accession numbers: KC189966–KC189969). They differed from one another, and from the deposited *An. nili s.s.* haplotype by 30 to 95 nucleotide substitutions and many indels, suggesting high genetic divergence (Figure S2). Genetic distances between species

<table><caption>**Table 3.** Pairwise microsatellite *Fst* (below) and corresponding *Nm* (above) values between *An. nili s.l.* populations from Cameroon.</caption><thead><tr><th rowspan="2"></th><th colspan="5">An. nili s.s.</th><th>An. carnevalei</th><th>An. ovengensis</th><th>An. somalicus</th></tr><tr><th>Ako</th><th>Nkolbisson</th><th>Kentzou</th><th>Moloundou A</th><th>Moloundou B</th><th>Ekelemba</th><th>Afan-Essokyé</th><th>Nyabessan</th><th>Mbébé</th></tr></thead><tbody><tr><td>Ako</td><td>-</td><td>23.97</td><td>0.56</td><td>0.88</td><td>0.50</td><td>0.35</td><td>0.74</td><td>0.55</td><td>0.54</td></tr><tr><td>Nkolbisson</td><td>0.011***</td><td>-</td><td>0.51</td><td>0.78</td><td>0.48</td><td>0.33</td><td>0.67</td><td>0.49</td><td>0.51</td></tr><tr><td>Kentzou</td><td><strong>0.303***</strong></td><td><strong>0.313***</strong></td><td>-</td><td>1.07</td><td>0.40</td><td>0.28</td><td>0.28</td><td>0.55</td><td>0.26</td></tr><tr><td>Moloundou A</td><td><strong>0.221***</strong></td><td><strong>0.236***</strong></td><td><strong>0.160***</strong></td><td>-</td><td>0.62</td><td>0.40</td><td>0.40</td><td>1.16</td><td>0.30</td></tr><tr><td>Moloundou B</td><td><strong>0.350***</strong></td><td><strong>0.367***</strong></td><td><strong>0.343***</strong></td><td><strong>0.296***</strong></td><td>-</td><td>0.40</td><td>0.23</td><td>0.48</td><td>0.20</td></tr><tr><td>Ekelemba</td><td><strong>0.418***</strong></td><td><strong>0.435***</strong></td><td><strong>0.450***</strong></td><td><strong>0.397***</strong></td><td><strong>0.416***</strong></td><td>-</td><td>0.14</td><td>0.26</td><td>0.13</td></tr><tr><td>Afan-Essokyé</td><td><strong>0.274***</strong></td><td><strong>0.295***</strong></td><td><strong>0.483***</strong></td><td><strong>0.423***</strong></td><td><strong>0.580***</strong></td><td><strong>0.646***</strong></td><td>-</td><td>0.28</td><td>0.36</td></tr><tr><td>Nyabessan</td><td><strong>0.315***</strong></td><td><strong>0.341***</strong></td><td><strong>0.311***</strong></td><td><strong>0.201***</strong></td><td><strong>0.380***</strong></td><td><strong>0.531***</strong></td><td><strong>0.495***</strong></td><td>-</td><td>0.26</td></tr><tr><td>Mbébé</td><td><strong>0.313***</strong></td><td><strong>0.319***</strong></td><td><strong>0.476***</strong></td><td><strong>0.435***</strong></td><td><strong>0.575***</strong></td><td><strong>0.659***</strong></td><td><strong>0.452***</strong></td><td><strong>0.491***</strong></td><td>-</td></tr></tbody></table>

<footnote>*P<0.05;</footnote>

<footnote>**P<0.01;</footnote>

<footnote>***P<0.001; Bold characters denote a significant heterozygote deficiency after correction for multiple testing by the sequential Bonferroni procedure;</footnote>

<footnote>He: expected heterozygosity; Rs: allelic richness calculated based on minimum size of 10 individuals per population [31] - no amplification or no polymorphism detected; Fis are estimated following Weir & Cockerham [40].</footnote>

<footnote>doi:10.1371/journal.pone.0058862.t002</footnote>

<footnote>***P<0.001. Bold characters denote a significant heterozygote deficiency after correction for multiple testing by the sequential Bonferroni procedure.</footnote>

<footnote>doi:10.1371/journal.pone.0058862.t003</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 2. Neighbor-Joining (NJ) trees of *An. nili* s.l. D3 haplotypes from Cameroon.** H1...Hn: Haplotype 1 to Haplotype n.
doi:10.1371/journal.pone.0058862.g002</figcaption></figure>

were high, ranging from 5.8 to 18.6% which correspond to
divergence times between 2.63 and 8.45 myr assuming a
divergence rate of 2.2% per myr (above). Within *An. nili s.s.*,
genetic distances between haplotypes ranged from 1.8 to 12.9%,
implying approximate divergence time of 0.82 to 5.86 myr (Table
S2).

The topologies of the D3 and ITS2 phylogenetic trees were
congruent, showing *An. nili* group to be monophyletic although
relationship between species changed according to the marker. All
the trees revealed three major groups among *Anopheles nili s.s.*
specimens analyzed: 1) *An. nili s.s.* from savannah/degraded forest
(Ako and Nkolbisson), 2) *An. nili s.s.* from deep forest group 1
(Kentzou and Moloundou A) and 3) *An. nili s.s.* from deep forest
group 2 (Ekelemba and Moloundou B) (Figures 2 and 3).
Interestingly, the two forest clusters appeared more closely related
to *An. ovengensis* and other members of the group, than to *An. nili s.s.*
from savannah/degraded forest areas.

Mitochondrial DNA Polymorphism and Divergence

Mitochondrial DNA sequences were determined for a mini-
mum of 81 individual mosquitoes and summary statistics for
genetic polymorphism are given in Table 4. Alignment length was
603 bp and 319 bp for the *COII* and *ND4*, respectively. Variations
among sequences only included base substitutions. Absence of

indels and heterozygous sites therefore suggests that there was no
nuclear transposed copy.

In contrast to rDNA data, there was no strong support for 8
distinct genetic clusters, although the existence of distinct genetic
lineages within *An. nili s.s.* populations was confirmed. All the
phylogenetic trees were congruent, showing *An. nili* group as
monophyletic, with six clades: 1-*An. carnevalei*, 2-*An. ovengensis*, 3-*An. somalicus*, 4-*An. nili s.s.* from Ako and Nkolbisson 5-*An. nili s.s.* from Kentzou and Moloundou A (deep forest group 1), and 6-*An. nili s.s.* from Ekelemba and Moloundou B (deep forest group 2) (Figures 4 and 5). However, the relationship between these clades was not clear as their position changed according to the gene tree. As observed above, *An. nili s.s.* specimens sampled in the deep forest appeared more closely related to *An. ovengensis* than to savannah/degraded forest populations of *An. nili s.s.*

Sequence divergence between consensus haplotypes from each
population ranged from 0.4 to 6.2% and from 0.5 to 11.2% for
*COII* and *ND4* genes, respectively (Table S3). Because the
conventional mtDNA rate of molecular evolution of 2% per myr
is used exclusively for lineages separated by less than 8% sequence
divergence, only *COII* divergence time was estimated and found to
vary approximately between 0.2 and 3 myr. All new mtDNA
sequences have been deposited to GenBank (Accession numbers:
KC189974 - KC190023).

<figure><img src="image_4.png" /><figcaption>**Figure 3. Neighbor-Joining (NJ) trees of *An. nili* s.l. ITS2 haplotypes from Cameroon.** H1...Hn: Haplotype 1 to Haplotype n.
doi:10.1371/journal.pone.0058862.g003</figcaption></figure>
<table><caption>**Table 4.** Summary statistics for *COII* and *ND4* polymorphism in *An. nili* s.l. populations from Cameroon.</caption><thead><tr><th rowspan="2"></th><th colspan="6">An. nili s.s.</th><th>An. carnevalei</th><th>An. ovengensis</th><th>An. somalicus</th></tr><tr><th>Ako</th><th>Nkolbisson</th><th>Kentzou</th><th>Moloundou A</th><th>Moloundou B</th><th>Ekelemba</th><th>Afan-Essokyé</th><th>Nyabessan</th><th>Mbébé</th></tr></thead><tbody><tr><td colspan="10"><strong>COII</strong></td></tr><tr><td>n</td><td>9</td><td>9</td><td>10</td><td>10</td><td>9</td><td>10</td><td>7</td><td>8</td><td>10</td></tr><tr><td>ns</td><td>603</td><td>603</td><td>603</td><td>603</td><td>603</td><td>603</td><td>603</td><td>603</td><td>603</td></tr><tr><td>s</td><td>0</td><td>0</td><td>6</td><td>3</td><td>6</td><td>1</td><td>12</td><td>3</td><td>7</td></tr><tr><td>h</td><td>1</td><td>1</td><td>5</td><td>3</td><td>2</td><td>2</td><td>5</td><td>4</td><td>6</td></tr><tr><td>dh</td><td>0</td><td>0</td><td>0.667</td><td>0.378</td><td>0.500</td><td>0.533</td><td>0.857</td><td>0.786</td><td>0.867</td></tr><tr><td>π</td><td>0</td><td>0</td><td>0.002</td><td>0.001</td><td>0.005</td><td>0.001</td><td>0.011</td><td>0.002</td><td>0.003</td></tr><tr><td>k</td><td>0</td><td>0</td><td>1.200</td><td>0.600</td><td>3.000</td><td>0.533</td><td>6.000</td><td>1.000</td><td>1.911</td></tr><tr><td>D<sub>T</sub></td><td>nc</td><td>nc</td><td>-1.796*</td><td>-1.562</td><td>1.566</td><td>1.302</td><td>1.225</td><td>-0.431</td><td>-0.964</td></tr><tr><td>D</td><td>nc</td><td>nc</td><td>-2.081*</td><td>-1.784</td><td>1.372</td><td>0.804</td><td>0.747</td><td>-0.622</td><td>-1.131</td></tr><tr><td>F</td><td>nc</td><td>nc</td><td>-2.256*</td><td>-1.934</td><td>1.580</td><td>1.026</td><td>0.938</td><td>-0.642</td><td>-1.224</td></tr><tr><td>F<sub>s</sub></td><td>nc</td><td>nc</td><td>-1.803*</td><td>-0.459</td><td>5.000</td><td>1.029</td><td>0.648</td><td>-1.164</td><td>-1.953*</td></tr><tr><td colspan="10"><strong>ND4</strong></td></tr><tr><td>n</td><td>8</td><td>8</td><td>10</td><td>10</td><td>8</td><td>10</td><td>7</td><td>8</td><td>10</td></tr><tr><td>ns</td><td>319</td><td>319</td><td>319</td><td>319</td><td>319</td><td>319</td><td>319</td><td>319</td><td>319</td></tr><tr><td>s</td><td>0</td><td>0</td><td>5</td><td>3</td><td>9</td><td>1</td><td>2</td><td>17</td><td>3</td></tr><tr><td>h</td><td>1</td><td>1</td><td>5</td><td>4</td><td>5</td><td>2</td><td>3</td><td>4</td><td>5</td></tr><tr><td>hd</td><td>0</td><td>0</td><td>0.822</td><td>0.533</td><td>0.857</td><td>0.200</td><td>0.524</td><td>0.750</td><td>0.822</td></tr><tr><td>π</td><td>0</td><td>0</td><td>0.005</td><td>0.002</td><td>0.016</td><td>0.001</td><td>0.002</td><td>0.020</td><td>0.003</td></tr><tr><td>k</td><td>0</td><td>0</td><td>1.489</td><td>0.600</td><td>5.035</td><td>0.200</td><td>0.571</td><td>7.929</td><td>1.022</td></tr><tr><td>D<sub>T</sub></td><td>nc</td><td>nc</td><td>-0.631</td><td>-1.562</td><td>1.513</td><td>-1.112</td><td>-1.237</td><td>1.079</td><td>-1.063</td></tr><tr><td>D</td><td>nc</td><td>nc</td><td>-0.686</td><td>-1.784</td><td>1.161</td><td>-1.243</td><td>-1.296</td><td>1.571*</td><td>-0.338</td></tr><tr><td>F</td><td>nc</td><td>nc</td><td>-0.754</td><td>-1.933</td><td>1.378</td><td>-1.347</td><td>-1.374</td><td>1.619*</td><td>-0.577</td></tr><tr><td>F<sub>s</sub></td><td>nc</td><td>nc</td><td>-1.284</td><td>-1.964</td><td>0.807</td><td>-0.339</td><td>-0.922</td><td>3.474</td><td>-2.20**</td></tr></tbody></table>

nc: not calculated. n: number of sequences; s: number of polymorphic sites; h: number of haplotypes; hd: haplotype diversity; π: nucleotide diversity. DT: Tajima's D; D: Fu and Li's D; F: Fu and Li's F; Fs: Fu's Fs statistics. *P<0.05, **P<0.01 and ***P<0.001.

doi:10.1371/journal.pone.0058862.t004

Discussion

Genotypes at 11 microsatellite loci and sequence variation in two rDNA regions (D3, ITS2) and two mtDNA genes (*COII*, *ND4*) were used to assess the level of genetic variability and differentiation among species of the *An. nili* group, and between *An. nili s.s.* populations from different ecological settings in South Cameroon. All the markers provided support for the current classification in the group and revealed further cryptic genetic diversity within the major vector *An. nili s.s.*, despite morphological similarities.

Heterologous microsatellite loci are commonly used across many insect species, particularly in *Anopheles* [51–54]. In fact, the use of primers originally designed to amplify DNA from a species in a different one can save time and resources, especially for microsatellite loci. However, depending on the genetic distance between species, microsatellite loci could either readily amplify, fail to amplify or show low polymorphism due to the presence of mismatches in priming sites within flanking sequences or absence of nucleotide repeat units variability [52,55–57]. Therefore, it is advocated to screen heterologous loci in few samples to detect suitable markers, before their use in large population genetic studies.

In the present study, 11 *An. nili s.s.* loci were screened for use in *An. carnevalei*, *An. somalicus* and *An. ovengensis* and up to 70% cross-amplified in all species. This provides a clear indication that the

four species are closely related, since the probability of cross-amplification of microsatellite loci is inversely related with genetic distance between taxa [55,56]. However, the large differences in microsatellite allele frequencies distribution and sequence profiles, resulting in high pairwise genetic differentiation indices and low migration rates fully support speciation between *An. nili s.s.*, *An. carnevalei*, *An. somalicus* and *An. ovengensis*, in agreement with previous studies using allozyme and rDNA markers [8,58]. Besides, although a high proportion of *An. nili s.s.* loci amplified in the sibling species of the group, only 6 in both *An. carnevalei* and *An. ovengensis*, and 5 in *An. somalicus* showed the qualities that make them suitable genetic markers (e.g., lack of null alleles, high genetic variability, conformity to HWE, statistical independence) and could be used for population genetic studies in these species.

Previously unrecognized cryptic genetic diversity not correlated to geographical distance between sampling sites was evidenced within the nominal species of the group *An. nili s.s.*, by detecting highly divergent lineages in the deep equatorial forest, despite morphological similarities of all specimens analyzed. This genetic substructure was detected by a range of molecular markers with different evolutionary dynamics, providing strong support for phylogenetic inferences. Indeed, rDNA haplotypes found in these populations differed from one another, and from the *An. nili s.s.* haplotypes from savannah/degraded forest areas by a large number of fixed mutations and indels leading to genetic distance
<figure><img src="image_5.png" /><figcaption>**Figure 4. Neighbor-Joining (NJ) trees of An. nili s.l. COII haplotypes from Cameroon.** H1...Hn: Haplotype 1 to Haplotype n.
doi:10.1371/journal.pone.0058862.g004</figcaption></figure>

estimates 4 to 8 fold higher than those commonly reported among cryptic *Anopheles* species [8,59–64]. According to accepted rates of molecular evolution of 2.2% per myr for ITS2 [48] and 2% per myr for mtDNA [49], we estimated the divergence time of the new molecular forms (lineages) at about 0.8 to 6 Myr or 0.2 to 3 Myr, respectively. This high genetic divergence within the *An. nili* group suggests that these taxa might actually belong to distinct species groups. However, these results generated using a small sample size and on a limited sampling range should be considered as preliminary. Extending the study area throughout the equatorial forest block of Central Africa will allow improving knowledge on the taxonomic status and the geographic distribution of forest *An*.

*nili* s.l. populations, their ecology and behaviour, and their epidemiological role in malaria parasite transmission.

The evergreen forest block of Central Africa is known to be of low overall quality for the development of *An. nili s.s.* which is replaced in this area by the other members of the group, namely *An. carnevalei* and *An. ovengensis* [9,18,65]. *Anopheles nili s.s.* is rather more frequent at the savannah/degraded forest ecotone, where it is characterized by shallow population structure with a weak effect of geographical distance on genetic differentiation [18]. Violation of MDE in deep forest populations suggests unstable demographic history for the species in the evergreen forest block of Central Africa. The low dispersal of *An. nili s.s.* in this environment probably restricts opportunities for genetic exchange between
<figure><img src="image_6.png" /><figcaption>**Figure 5. Neighbor-Joining (NJ) trees of *An. nili* s.l. ND4 haplotypes from Cameroon.** H1...Hn: Haplotype 1 to Haplotype n.
doi:10.1371/journal.pone.0058862.g005</figcaption></figure>

geographically isolated populations, prompting for genetic diver-
gence under the action of increased genetic drift [18,66]. This is in
agreement with the low genetic variability observed in forest
populations suggesting small effective population size, coupled
with high levels of genetic differentiation corresponding to weak
migration index estimated between populations. Moreover, the
fact that many of the microsatellite markers that were polymorphic
in samples from Ako and Nkolbisson did not amplify at all or were
not polymorphic (i.e. only one allele could be detected) in deep
forest populations is also consistent with the action of increased
genetic drift in the later populations. Finally, by demonstrating
physical independence between markers, recent mapping of

microsatellite loci on the chromosomal arms of *An. nili* provided
further support for a major role of genetic drift in shaping the
population genetic structure of *An. nili* in Central Africa, reflected
in genome-wide differentiation patterns [27]. The development of
additional molecular markers is now needed to allow more fine-
grained analysis of the level of genetic polymorphism and its
distribution along the genome of *An. nili* in search for signatures of
selection in ecologically-relevant genes or genomic regions.
Indeed, original information gained on the genetic structure, gene
flow and species diversity within this major group of malaria
vectors can further be used to explore global processes such as the
genetics of ecological adaptation, speciation and susceptibility to
*Plasmodium*, within a comparative framework that will use and complement information available for other major human malaria vectors such as mosquitoes from the *An. gambiae* complex and the *An. funestus* group.

Moreover, the recent findings of circulation of *P. falciparum* along with other *Plasmodium* species in great apes in Central African forest region [24,25] raise concerns about pathogen transfer between humans and primates and highlight the need to improve knowledge of forest malaria vectors. Although it is not yet known whether the gorilla carriers of *P. falciparum* acquired the parasite from humans (or from other primates), nor the pathogenicity of these infections in gorillas, the continuously increasing contact between humans and primate populations, mostly due to logging and deforestation, increases opportunities for transmission of new pathogens from primates to humans and vice versa [25]. The possibility that forest mosquitoes such as members of the *An. nili* group might contribute to the emergence of zoonosis in humans needs to be assessed.

## Supporting Information

**Figure S1** Alignment of D3 sequences of *An. nili* s.l. haplotypes. AN: *An. nili* s.s.; AO: *An. ovengensis*; AC: *An. carnevalei*; AS: *An. somalicus*. AK: Ako; NK: Nkolbisson; KT: Kentzou; MOA: Moloundou A; MOB: Moloundou B; EK: Ekelemba; AE: Afan-Essokyé; NY: Nyabessan; MB: Mbébé. (PDF)

**Figure S2** Alignment of ITS2 sequences of *An. nili* s.l. haplotypes. AN: *An. nili* s.s.; AO: *An. ovengensis*; AC: *An. carnevalei*;

AS: *An. somalicus*. AK: Ako; NK: Nkolbisson; KT: Kentzou; MOA: Moloundou A; MOB: Moloundou B; EK: Ekelemba; AE: Afan-Essokyé; NY: Nyabessan; MB: Mbébé. (PDF)

**Table S1** Estimates of null alleles frequencies per locus and per geographical population of *An. nili* s.l. from Cameroon. (DOCX)

**Table S2** Pairwise genetic distance estimates between *An. nili* s.l. D3 (above the diagonal) and ITS2 (below the diagonal) haplotypes from Cameroon. (DOCX)

**Table S3** Pairwise genetic distance estimates between *An. nili* s.l. ND4 (above the diagonal) and COII (below the diagonal) haplotypes from Cameroon. (DOCX)

## Acknowledgments

The authors are grateful to Billy Tene for help during field collections, Anna Cohuet and Diego Ayala for support during genetic data analysis, and Daniel Couret for map conception.

## Author Contributions

Critically revised the manuscript: FS IM IS DF CAN. Conceived and designed the experiments: CN FS IS DF CAN. Performed the experiments: CN PK PAA CAN. Analyzed the data: CN FS CAN. Contributed reagents/materials/analysis tools: IM IS DF CAN. Wrote the paper: CN.

## References

1. World Health Organization (WHO). World malaria report 2011 (2011) Available on http://www.who.int/malaria/world_malaria_report_2011.

2. Fontenille D, Simard F (2004) Unravelling complexities in human malaria transmission dynamics in Africa through a comprehensive knowledge of vector populations. Comp Immun Microbiol Infect Dis 27: 357-375.

3. Pates H, Curtis C (2005) Mosquito behavior and vector control. Ann Rev Entomol 50: 53-70.

4. Greenwood B (2009). Can malaria be eliminated? Trans Roy Soc Trop Med Hyg: 103: S2-S5.

5. Gillies MT, De Meillon B (1968) The Anophelinae of Africa South of the Sahara (Ethiopian Zoogeographical Region), 2nd ed. The South African Institute for Medical Research, Johannesburg.

6. Brunhes G, Le Goff G, Geoffroy B (1999) Afro-tropical anopheline mosquitoes. III. Description of three new species: Anopheles carnevalei sp. Nov., and A. hervyi sp. Nov., and A. dualaensis sp. Nov. and resurrection of A. rageaui Mattingly and Adam. J Am Mosq Control Assoc 15 (4): 552-558.

7. Awono-Ambene HP, Kengue P, Simard F, Antonio-Nkondjio C, Fontenille D (2004) Description and biometrics of Anopheles (Cellia) ovengensis (Diptera: Culicidae): a new malaria vector species of the Anopheles nili group from south Cameroon. J Med Entomol 41: 561-568.

8. Kengne P, Awono-Ambene HP, Antonio-Nkondjio C, Simard F, Fontenille D (2003) Molecular identification of the Anopheles nili group African malaria vectors. Med Vet Entomol 17: 67-74.

9. Antonio-Nkondjio C, Ndo C, Costantini C, Awono-Ambene P, Fontenille D, Simard F (2009) Distribution and larval habitat characterization of Anopheles nili and An. moucheti along river networks in south Cameroon. Acta Trop 112: 270-276.

10. Hamon J, Mouchet J (1961) Les vecteurs secondaires du paludisme humain en Afrique. Med Trop 21: 643-60.

11. Carnevale P, Le Goff G, Toto JC, Robert V (1992) Anopheles nili as the main vector of human malaria in villages of southern Cameroon. Med Vet Entomol 6: 135-8.

12. Antonio-Nkondjio, C, Awono-Ambene P, Toto JC, Meunier JY, Zebaze-Kemleu S, et al. (2002a) High malaria transmission intensity in a village close to Yaoundé, the capital city of Cameroon. J Med Entomol 39: 350-355.

13. Dia I, Diop T, Rakotoarivony I, Kengne P, Fontenille D (2003) Bionomics of Anopheles gambiae Giles, A. arabiensis Patton, A. funestus Giles and A. nili (Theobald) (Diptera: Culicidae) and transmission of Plasmodium falciparum in a Sudano-Guinean zone (Ngari, Senegal). J Med Entomol 40: 279-83.

14. Adja AM, N'goran KE, Kengne P, Koudou GB, Toure M, et al. (2006) Transmission vectorielle du paludisme en savane arborée à Gansé en Côte d'Ivoire. Med Trop 66: 449-455.

15. Adja AM, N'goran EK, Koudou GB, Dia I, Kengne P, et al. (2011) Contribution of Anopheles funestus, An. gambiae and An. nili (Diptera: Culicidae) to the perennial malaria transmission in the southern and western forest areas of Cote d'Ivoire. Ann Trop Med Parasitol 105(1): 13-24.

16. Awono-Ambene P, Antonio-Nkondjio C, Toto JC, Ndo C, Etang J, et al. (2009) Epidemiological importance of the Anopheles nili group of malaria vectors in equatorial villages of Cameroon, Central Africa. Sci Med Afr 1: 13-20.

17. Ridl FC, Bass C, Torrez M, Govender D, Ramdeen V, et al. (2008) A pre-intervention study of malaria vector abundance in Rio Muni, Equatorial Guinea: their role in malaria transmission and the incidence of insecticide resistance alleles. Malar J 7: 194.

18. Ndo C, Antonio-Nkondjio C, Cohuet A, Ayala D, Kengne P, et al. (2010) Population genetic structure of the malaria vector Anopheles nili in sub-Saharan Africa. Malar J 9: 161.

19. Rivola E, Holstein MH (1957) Note sur une variété d'Anopheles nili Theo. Bull Soc Path Exo 50: 382-387.

20. Antonio-Nkondjio C, Kerah Hinzoumbe C, Simard F, Awono-Ambene P, Tchuinkam T, et al. (2006) Complexity of the malaria vectorial system in Cameroon: Contribution of secondary vectors to malaria transmission. J Med Entomol 43: 1215-1221.

21. Ollomo B, Durand P, Prugnolle F, Douzery E, Arnathau C, et al. (2009) A new malaria agent in African hominids. PLoS Pathog 5: e1000446.

22. Duval L, Fourment M, Nerrienet E, Rousset D, Sadeuh SA, et al. (2010) African apes as reservoirs of Plasmodium falciparum and the origin and diversification of the Laverania subgenus. Proc Natl Acad Sci USA 107: 10561-10566.

23. Liu W, Li Y, Learn GH, Rudicell RS, Robertson JD, et al. (2010) Origin of the human malaria parasite Plasmodium falciparum in gorillas. Nature 467: 420-5.

24. Prugnolle F, Durand P, Neel C, Ollomo B, Ayala FJ, et al. (2010) African great apes are natural hosts of multiple related malaria species, including Plasmodium falciparum. Proc Natl Acad Sci USA 107: 1458-63.

25. Prugnolle F, Ollomo B, Durand P, Yalcindag E, Arnathau C, et al. (2011) African monkeys are infected by Plasmodium falciparum nonhuman primate-specific strains. PNAS 108 (29): 11948-11953.

26. Sharakhova MV, Antonio-Nkondjio C, Xia A, Ndo C, Awono-Ambene P, et al. (2011) Cytogenetic map for Anopheles nili: Application for population genetics and comparative physical mapping. Infect Genet Evol 11: 746-754.

27. Peery A, Sharakhova M, Antonio-Nkondjio C, Ndo C, Weill M, et al. (2011) Improving the population genetics toolbox for the study of the African malaria vector Anopheles nili: microsatellite mapping chromosomes. Para Vect 4: 202.

28. Gillies MT, Coetzee M (1987) A Supplement to the Anophelinae of Africa south of the Sahara. South African Institute of Medical Research, Johannesburg, p. 143.
29. Morlais I, Ponçon N, Simard F, Cohuet A, Fontenille D (2004) Intraspecific nucleotide variation in Anopheles gambiae: new insights into the biology of malaria vectors. *Am J Trop Med Hyg* 71(6): 795-802.

30. Berthomieu A, Kengne P, Awono-Ambene HP, Raymond M, Fontenille D, et al. (2003) Isolation and characterization of microsatellite DNA markers in the malaria vector Anopheles nili. *Mol Ecol* Notes 3: 392-393.

31. El Mousadik A, Petit R (1996) High level of genetic differentiation for allelic richness among populations of the argan tree (Argania spinosa L.) Skeels endemic to Morocco. *Theo Appl Gen* (92): 832-839.

32. Goudet J (1995) FSTAT version 2.9.3.2. A computer software to calculate F-statistics. *J Hered* 86: 485-486.

33. Raymond M, Rousset F (1995) GENEPOP, Version 1.2. A population genetics software for exact tests and ecumenicism. *J Hered* 86: 248-249.

34. Rousset F (2007) *Genepop'007* a complete re-implementation of the *Genepop* software for windows and linux. *Mol Ecol* Notes 8: 103-106.

35. Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P (2004) MICRO-CHECKER: Software for identifying and correcting genotyping errors in microsatellite data. *Mol Ecol* 4: 535-538.

36. Cornuet JM, Luikart G (1996) Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. *Genetics* 144: 2001-2014.

37. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. *Genetics* 155: 945-959.

38. Evanno G, Regnaut S, Goudet J (2005) Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. *Mol Ecol* 14: 2611-2620.

39. Wright S (1978) *Evolution and Genetics of populations*. Variability within and among natural populations. Vol 4, Chicago, University of Chicago press.

40. Weir BS, Cockerham CC (1984) Estimating F-statistics for the analysis of population structure. *Evolution* 38: 1358-1370.

41. Slatkin M (1995) A measure of population subdivision based on microsatellite allele frequencies. *Genetics* 139: 457-462.

42. Rousset F (1997) Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. *Genetics* 145(4): 1219-1228.

43. Rousset F (2007) *Genepop'007* a complete re-implementation of the *Genepop* software for windows and linux. *Mol Ecol* Notes 8: 103-106.

44. Holm S (1979) A simple sequentially rejective multiple test procedure. *Scand J Stat* 6: 65-70.

45. Librado P, Rozas J (2009) DnaSP v5: a software for comprehensive analysis of DNA data. *Bioinformatics* 25: 1451-1452.

46. Kumar S, Tamura K, Nei M (2004) MEGA 3: integrate software for Molecular Evolutionary Genetics Analysis and sequence alignment. *Brief Bioinformatics* 5: 150-163.

47. Tajima F, Nei M (1984) Estimation of evolutionary distance between nucleotide sequences. *Mol Biol Evol* 1(3): 269-85.

48. Schlötterer C, Hausser MT, Von Haeseler A, Tautz D (1994) Comparative evolutionary analysis of rDNA ITS regions in Drosophila. *Mol Biol Evol* 11 (3): 513-522.

49. Desalle R, Templeton A, Mori I, Pletscher S, Johnston JS (1987) Temporal and spatial heterogeneity of mtDNA polymorphism in natural populations of Drosophila mercatorum. *Genetics* 116: 215-223.

50. Saitou N, Nei M (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees. *Mol Biol Evol* 4: 406-425.

51. Walton C, Thelwell TJ, Priestman A, Butlin RK (1998) The use of microsatellites to study gene flow in natural populations of Anopheles malaria vectors in Africa: potential and pitfalls. *J Am Mosq. Control Assoc* 14: 266-272.

52. Kamau L, Lehmann T, Hawley WA, Orago ASS, Ke Z, et al. (1998) Use of short tandem repeats for analysis of genetic variability in sympatric populations of An. gambiae and An. arabiensis. *Heredity* 80: 675-682.

53. Simard F, Fontenille D, Lehmann T, Girod R, Brutus L (1999). High amounts of genetic differentiation between populations of the malaria vector Anopheles arabiensis from West Africa and eastern outer islands. *Am J Trop Med Hyg* 60: 1000-1009.

54. Antonio-Nkondjio C, Ndo C, Kengne P, Mukwaya L, Awono-Ambene HP, et al. (2008) Population structure of the malaria vector Anopheles moucheti in the equatorial forest region of Africa. *Malar J* 7: 120.

55. Weber JL, Kwitek AE, May PE, Killyam AM (1990) Dinucleotide repeat polymorphism at the, D4s 174-locus. *Nucl Acids Res* 18: 4636-4638.

56. Primmer CR, Moller AP, Ellegren H (1996) A wide-range survey of cross-species microsatellite amplification in birds. *Mol Ecol* 5: 365-378.

57. Deitz KC, Reddy VP, Reddy MR., Satyanarayanan N, Lindsey, et al. (2012) Limited usefulness of microsatellite markers from the malaria vector Anopheles gambiae when applied to the closely related species Anopheles melas. *J Hered* 103 (4): 585-593.

58. Awono-Ambene HP, Simard F, Antonio-Nkondjio C, Cohuet A, Kengne P, et al. (2006). Multilocus enzyme electrophoresis supports speciation within the Anopheles nili group of Malaria vectors in Cameroon. *Am J Trop Med Hyg* 75: 656-658.

59. Paskewitz SM, Ng K, Coetzee M, Hunt RH (1993) Evaluation of the polymerase chain reaction method for identifying members of the Anopheles gambiae (Diptera: Culicidae) complex in Southern Africa. *J Med Entomol* 30: 953-957.

60. Scott JA, Brogdon WG, Collins FH (1993) Identification of single specimens of the Anopheles gambaie complex by the polymerase chain reaction. *Am J Trop Med Hyg* 49 (4): 520-529.

61. Collins FH, Paskewitz SM (1996) A review of the use of ribosomal DNA (rDNA) to differentiate among cryptic Anopheles species. *Ins Mol Biol* 5 (1): 1-9.

62. Cohuet A, Simard F, Toto JC, Kengne P, Coetzee M, et al. (2003) Species identification within the Anopheles funestus group of malaria vectors in Cameroon and evidence for a new species. *Am J Trop Med Hyg* 69 (2): 200-205.

63. Kengne P, Antonio-Nkondjio C, Awono-Ambene HP, Simard F, Awolola TS, et al. (2007) Molecular differentiation of three closely related members of the mosquito species complex, Anopheles moucheti, by mitochondrial and ribosomal DNA polymorphism. *Med Vet Entomol* 21: 177-182.

64. Manonmani AM, Sadanandane C, Sahu SS, Mathivanan A, Jambulingam P (2007) rDNA-ITS2-PCR assay for grouping the cryptic species of Anopheles culicifacies complex (Diptera: Culicidae). *Acta Trop* 104 (1): 72-77.

65. Ayala D, Costantini C, Ose K, Kamdem GC, Antonio-Nkondjio C, et al. (2009) Habitat suitability and ecological niche profile of major malaria vectors in Cameroon. *Malar J* 8: 307.

66. Le Goff G, Carnevale P, Robert V (1997) Low dispersion of anopheline malaria vectors in the african equatorial forest. *Parasite 2*: 187-189.

67. Silver JB (2008) Mosquito ecology: field sampling methods. Springer, 3rd Edition.