Aedes albopictus (Skuse, 1884) (Diptera: Culicidae), a mosquito native to Asia, has recently invaded all five continents. In Central Africa it was first reported in the early 2000s, and has since been implicated in the emergence of arboviruses such as dengue and chikungunya in this region. Recent genetic studies of invasive species have shown that multiple introductions are a key factor for successful expansion in new areas. As a result, phenotypic characters such as vector competence and insecticide susceptibility may vary within invasive pest species, potentially affecting vector efficiency and pest management. Here we assessed the genetic variability and population genetics of Ae. albopictus isolates in Cameroon (Central Africa), thereby deducing their likely geographic origin.
Methods and Results
Mosquitoes were sampled in 2007 in 12 localities in southern Cameroon and analyzed for polymorphism at six microsatellite loci and in two mitochondrial DNA regions (ND5 and COI). All the microsatellite markers were successfully amplified and were polymorphic, showing moderate genetic structureamong geographic populations (FST = 0.068, P<0.0001). Analysis of mtDNA sequences revealed four haplotypes each for the COI and ND5 genes, with a dominant haplotype shared by all Cameroonian samples. The weak genetic variation estimated from the mtDNA genes is consistent with the recent arrival of Ae. albopictus in Cameroon. Phylogeographic analysis based on COI polymorphism indicated that Ae. albopictus populations from Cameroon are related to tropical rather than temperate or subtropical outgroups.
Citation: Kamgang B, Brengues C, Fontenille D, Njiokou F, Simard F, Paupy C (2011) Genetic Structure of the Tiger Mosquito, Aedes albopictus, in Cameroon (Central Africa). PLoS ONE 6(5): e20257. https://doi.org/10.1371/journal.pone.0020257
Editor: Pedro Lagerblad Oliveira, Universidade Federal do Rio de Janeiro, Brazil
Received: March 16, 2011; Accepted: April 22, 2011; Published: May 24, 2011
Copyright: © 2011 Kamgang 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 research was supported by the French government (Agence Nationale pour la Recherche) through the Epidengue project (ANR 05 SEST 010 -01). Basile Kamgang was supported by an Institut de Recherche pour le Développement/Département Soutien et Formation PhD fellowship. 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.
The Asian “tiger” mosquito, Aedes (Stegomyia) albopictus (Skuse, 1894) (Diptera: Culicidae), is a major vector of dengue (DENV, Flaviriridae, flavivirus) and chikungunya (CHIKV, Togaviridae, alphavirus) viruses in tropical areas –. This mosquito species originated in southeast Asia , and has rapidly invaded numerous countries in Africa, the Middle East, Europe and the Americas over the last three decades . This rapid global spread was favored by international trade, especially in used tires . In Africa, Aedes albopictus was first reported in South Africa in 1989, where it was immediately controlled , and two years later in Nigeria, where it pullulated . In Central Africa, Ae. albopictus was first recorded in Cameroon in 2000 , followed by rapid expansion in most of the south . Aedes albopictus is now the dominant species at the periphery of urban centers such as Yaoundé and Douala , where it has gradually replaced the resident species, Aedes aegypti (L., 1762).
Coinciding with the spread of Ae. albopictus, numerous DENV and CHIKV epidemics occurred in Cameroon  and Gabon  during the 2000s. Ae. albopictus was incriminated as the main vector of both viruses during a 2007 dual outbreak in Libreville, Gabon , . In addition, Ae. albopictus isolates in Cameroon proved to be orally susceptible to DEN-2 virus and highly competent for CHIKV .
Aedes albopictus exhibits strong physiological and ecological plasticity, allowing it to thrive in a wide range of climates and habitats. Known as a sylvatic species in Asia , Ae. albopictus has now adapted to human environments and is preferentially found in suburban areas, although it has also been recorded in densely populated urban areas such as Rome, Italy . The range of larval development sites suitable for Ae. albopictus is extremely broad and includes natural water collections (e.g., bamboo stumps, bromeliads and tree holes) and artificial containers (e.g., water storage containers, abandoned car parts and all other discarded items where rainwater can accumulate) . The mosquito preferentially feeds on mammals , but females can feed on most groups of cold- to warm-blooded vertebrates, including reptiles, birds and amphibians . The feeding preferences of Ae. albopictus vary with the geographic origin of local mosquito populations . This ability of Ae. albopictus to feed on different animal species enhances not only its biological traits (e.g., fecundity and survival) but also its capacity to propagate zoonotic pathogens across animal species and to humans . Finally, the broad host spectrum would facilitate invasion of diverse environments ranging from forests to urban areas.
Several phylogeographic studies have been undertaken to determine the origin of invasive Ae. albopictus populations. An early study based on isoenzymatic markers suggested that North and South American populations (USA and Brazil) were closely related to Japanese populations . This was confirmed by a study suggesting that “temperate populations” in the USA, Japan and Italy form a genetic cluster distinct from tropical Asian populations . More recently, however, studies of mitochondrial markers challenged the common origin of North and South American populations , and suggested that Brazilian populations were related to Southeast Asian populations rather than to temperate Asian populations . Differences in biological traits, such as the winter diapause exhibited by USA populations but not by Brazilian populations , support this hypothesis. It has also been shown that the vector competence of Ae. albopictus varies according to its geographic origin –. As biological traits of Ae. albopictus are genetically determined, and as these traits may influence virus transmission in newly colonized areas, it is important to determine the geographic origin of invading populations.
As little is known of the geographic origin or genetics of Ae. albopictus in Africa, we undertook a genetic study of Ae. albopictus in Cameroon.
Materials and Methods
An institutional clearance for this study, including the sampling of mosquito, was granted by Health Ministry of Cameroon.
Aedes albopictus mosquitoes were sampled as larvae or pupae in 12 localities in South Cameroon during April 2007. The geographic locations and main characteristics of the sampling sites are summarized in Figure 1 and Table 1. Aquatic stages (field generation, F0) were collected from domestic containers (tanks) and peri-domestic water collections (used tires, discarded containers, etc.) and forwarded in plastic boxes to insectaries. Each sample from a given locality consisted of pooled larvae from 4–6 breeding places, to avoid inbreeding. In the insectaries, mosquito larvae and pupae were transferred to pans containing water for rearing to adulthood. Adults were examined morphologically to confirm they were Ae. albopictus and were then stored in individual tubes containing a desiccant at −20°C until molecular analysis.
Total DNA was extracted from whole mosquito bodies, as described by Morlais et al. . DNA pellets were resuspended in sterile water and stored at −20°C. Genetic polymorphism was assessed at six microsatellite markers: AealbA9, AealbB6, AealbB51, AealbB52, AealbD2 and AealbF3 .
DNA was amplified in a GeneAmp 9600 thermal cycler (PerkinElmer, USA) in 25-µl reaction mixes containing 2 µl of 1/5-diluted total DNA, 2.5 µl of 10× reaction buffer (Qiagen, USA), 0.5 mM MgCl2 (for AealbB6 and AealbD2), 0.2 mM each dNTP (Eurogentec, Belgium), 10 pmol of each primer, and 0.5 U of Taq polymerase (Qiagen, USA). The 5′ end of the forward primer was labeled with a fluorescent dye (FAM, NED, PET and VIC). Cycling conditions were as follows: 5 min at 94°C, followed by 35 cycles of 45 s at 94°C, 30 s at the annealing temperature (TA), 45 s at 72°C, and 30 min at 72°C. PCR products were diluted in water (1/50 to 1/200, according to the sensitivity of each marker) and pooled with other compatible products according to allele size range and fluorescent dye. Pools were prepared by adding 0.4 µl of GS 500 Liz internal size standard™ (Applied Biosystems, USA) and HD formamide to 1 µl of each diluted amplification product, for a total volume of 20 µl. The mixture was heated at 95°C for 3 minutes before migration in an ABI Prism™ 3130 XL automatic sequencer (Applied Biosystems). Microsatellite alleles were scored with GeneMapper software (Applied Biosystems).
Genetic diversity per locus and per sample was characterized by estimating allele richness, Rs , and unbiased expected heterozygosity He , using FSTAT22.214.171.124  and GENETIX4.05.2 . The frequency of null alleles at each locus within each population, linkage disequilibrium between pairs of loci, deviations from Hardy-Weinberg equilibrium (HWE), and genetic differentiation indices were estimated with GENEPOP4.07 . FIS and FST estimates were calculated as described in , and were tested for statistical significance with exact tests available in GENEPOP4.07. The overall significance of multiple tests was estimated with Fisher's combined probability test. Critical significance levels for multiple testing were corrected with the sequential Bonferroni procedure . The frequency of null alleles at each locus within each population was determined with GENEPOP V4.07, and the allele and genotype frequencies were then adjusted accordingly in MICROCHECKER V2.2.3 . The null allele adjusted dataset was compared to the original dataset to explore the effect of null alleles on estimates of genetic differentiation.
Isolation by distance (IBD) was estimated with Mantel's test in GENEPOP4.0.7, using the correlation between genetic and geographic distances by the regression of pairwise FST/(1−FST) on the natural logarithm (Ln) of straight-line geographic distance . UPGMA trees were computed with POWERMARKER 3.25 , based on microsatellite Nei's genetic distance . The robustness of population trees was evaluated by bootstrapping over loci. One thousand bootstrapped trees were summarized to obtain a consensus tree, using the CONSENS module in the Phylip package , followed by editing with TREEVIEW .
A Bayesian approach available in STRUCTURE 2.2 was used to infer the number of clusters (K) in the data set, without prior information of the sampling locations . A model in which allele frequencies 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. We used 20 independent runs for each value of K (K = 1 to 12), with a burn-in period of 100 000 iterations and 400 000 replications. The method of Evanno et al.  was used to determine the most likely number of clusters. This approach uses an ad hoc quantity, ΔK, based on the second-order rate of change of the likelihood function between successive values of K.
Mitochondrial DNA analysis
At least ten individuals per locality, selected among those previously analyzed for microsatellite markers, were used to explore sequence polymorphism in the mitochondrial genes ND5 and COI. DNA samples were used as templates to amplify a 400-bp fragment of ND5 and a 550-bp fragment of COI. Two sets of primers were used: 5′TCCTTAGAATAAAATCCCGC-3′(FOR-ND5) and 5′-GTTTCTGCTTTAGTTCATTCTTC-3′ (REV-ND5) for ND5 , and 5′-GGAGGATTTGGAAATTGATTAGTTC-3′ (FOR-COI) and 5′-CCCGGTAAAATTAAAATATAAACTTC-3′ (REV-COI) for COI . PCR amplification used a GeneAmp 9600 thermal cycler (PerkinElmer) with 21-µl reaction mixes containing 4 µl of 1/50-diluted template DNA, 2.5 µl of 10× reaction buffer (Qiagen), 0.5 mM MgCl2 (for ND5), 0.2 mM dNTP, 10 pmol of each primer and 0.5 U of Taq polymerase (Qiagen). The thermal cycling conditions consisted of 3 min at 94°C, 35 cycles of 30 s at 94°C, 45 s at TA (52°C for ND5 and 54°C for COI), and 45 s at 72°C, followed by 10 min at 72°C. PCR products were detected by agarose gel electrophoresis in Tris-borate-EDTA buffer (TBE), stained with 0.5 µg/ml ethidium bromide in TBE, and visualized over UV light. PCR products were purified with the Ampure kit (Beckman Coulter, France, Paris) and sequenced on both strands using ABI Prism BigDye Terminator version 3.1 (Applied Biosystems). Each 10-µl sequencing reaction mix contained 1 µl of Ready Reaction mix (Applied Biosystems), 1 µl of 5× sequencing buffer, 3.2 pmol/µl primer, and 0.5–1.0 µl of purified PCR product. The reaction consisted of an initial denaturation step of 1 min at 96°C, followed by 25 cycles of 10 s at 96°C, 5 s at 50°C, and 4 min at 60°C, with a final elongation step of 3 min at 72°C. Sequencing reaction products were purified with the Agencourt CleanSeq kit (Beckman Coulter) and analyzed on an ABI 3130XL automatic sequencer (Applied Biosystem).
Sequences were analyzed with SEQSCAPE software 2.5 (Applied Biosytems) and aligned with Clustal W . ND5 and COI sequences were numbered according to the reference sequences (GeneBank ID AJ971026.1 and AJ971015.1, respectively). Basic sequence statistics, including the number of haplotypes per sample (Hp), the number of segregating sites (S), haplotype diversity (HpD), nucleotide diversity (π) and the average number of nucleotide differences were computed with DnaSP 4.10.9 . Neutrality tests used the statistics of Tajima , Fu and Li  and Fu  in DnaSP 4.10.9. The phylogenetic relationships between COI and ND5 haplotypes recorded in Cameroon and previously published sequences of Ae. albopictus from Asia, the Americas, the Indian Ocean and Europe were explored with Bayesian inference analysis. MrModeltest v2.2  was first used to select the model best fitting the COI and ND5 nucleotide sequence data (under Akaike's information criterion). Analyses were performed with MrBayes 3.1.2  and four Markov chains were run for 2 000 000 generations (sampling every 100 generations) to allow adequate time for convergence. The first 2500 trees were discarded as burn-in, and the remaining 17 500 sampled trees were used to estimate the 50% majority rule consensus tree and the Bayesian posterior probabilities. All MCMC runs were repeated twice to confirm consistent approximation of the posterior parameter distributions.
Genetic variability and Hardy Weinberg expectations.
A total of 452 Ae. albopictus specimens collected in 12 localities in Cameroon were successfully genotyped at all 6 microsatellite loci (Table 2). All loci were polymorphic, showing a number of distinct alleles ranging from 2 at locus AaelbB51 to 17 at loci AaelbA9 and AaelbB6. However, at AaelbB51, the dominant and shared allele across all samples was fixed in Garoua-Boulai, Abong-Mbang, Buea, Bamenda and Bafoussam. Over all the loci, the diversity indexes (Rs and He) were lowest in Bafoussam and highest in Buéa (Table 2). Departures from Hardy-Weinberg equilibrium associated with significant heterozygote frequency deficits were detected in 7 of 67 tests that could be conducted (one in AaebB51 and six in AaelbB6, see Table 2). All these deviations were attributed to the presence of null alleles, with frequencies of unobserved alleles ranging from 0.013 to 0.274 in loci with significant heterozygote deficiency (see Table S1). However, null alleles did not significantly bias our interpretation, as re-analysis of adjusted datasets gave similar results. Significant linkage disequilibrium was revealed in three occasions (AaelbA9-AaelbB6 and AaelbA9-AaelbD2 in Bamenda, and AaelbB6-AaelbD2 in Bertoua).
Overall genetic differentiation between samples was moderate and statistically significant (FST = 0.068, P<0.0001). Pairwise analysis of populations indicated significant genetic differentiation for most sample pairs (61/66), associated with FST estimates ranging from 0.006 (Bertoua-Abong-Mbang) to 0.176 (Buea-Yaounde) (Table 3). The UPGMA unrooted tree constructed with pairwise genetic distances and bootstrap values supporting node clearly indicated that the Buea sample was genetically distant from all the other samples (Figure 2). The topology of the tree shown in Figure 2 suggested that the pattern of genetic differentiation was not shaped by the geographic distance between sampling sites, in agreement with Mantel's test (P = 0.33). Bayesian cluster analysis identified 2 distinct genetic clusters in the dataset, although most of geographic populations appeared as a mixture of specimens from the two clusters except Buea and Garoua-Boulaï that appear made up of a single dominant cluster (see Figure S1).
The genetic relationship among 12 Ae. albopictus populations sampled in Cameroon is shown.
Mitochondrial DNA analysis
Nucleotide sequences of mtDNA ND5 (400 bp) and COI (500 bp) were obtained for 153 Ae. albopictus specimens sampled at the 12 Cameroonian locations. Overall, the polymorphism in the ND5 region was low, with only three segregating sites defining four closely related haplotypes, resulting in low haplotype and nucleotide diversity indexes (Table 4 and 5). The most frequent haplotype, H1 (86.3%), was detected in all geographic samples. Similarly, low levels of polymorphism were detected in the COI gene (Table 4 and 5). The most frequent haplotype, H1 (86.9%), was found in all geographic samples. Tajima's D and Fu's Fs statistics were negative for both genes (D = −0.8950 and Fs = −1.827 for ND5; D = −0.347 and Fs = −0.558 for COI), but not statistically significant. Negative values of these statistics might reflect either a selective sweep or recent demographic expansion .
In order to determine the geographic origin of Ae. albopictus populations invading Cameroon, phylogenetic relationships between COI and ND5 sequences recorded in Cameroon and previously published sequences were assessed by using Bayesian inference. COI sequences segregated into two lineages (Figure 3). The first lineage encompassed specimens from temperate and subtropical areas (France, Greece, USA, Reunion Island and Madagascar) and the second encompassed specimens from tropical areas (Brazil, India, Cambodia, Thailand and Vietnam), including all Cameroonian specimens. Interestingly, a SNP at position 322 (Table 4) distinguished all tropical populations (G at 322) from all temperate/subtropical populations (A at 322). All the sequences were monophyletic at ND5. The sequence of the most frequent haplotype, H1, perfectly matched a cosmopolitan haplotype previously described in South-East Asia (Thailand, Cambodia, Vietnam), the Indian Ocean (Reunion Island and Madagascar), Europe (France), Hawaï and the USA , . Haplotypes H2 and H4 have previously been described in Cameroon , whereas H3 is new.
The phylogeny was constructed using MrBayes 3.1.2, ngen = 2 000 000. Best-fitting models selected using MRmodeltest (AIC) were HKY for COI and HKY+I+G for the ND5 nucleotide dataset. Branch support is indicated by the posterior probability values. Accession numbers of COI and ND5 outgroup sequences are given in supporting information file Table S2.
This is the first study of the genetic variation and differentiation of Ae. albopictus populations invading continental Africa. Twelve mosquito population samples from Cameroon were analyzed using a set of six microsatellite DNA markers and two mtDNA genes (COI and ND5).
The set of microsatellite makers used here had already been tested on a single Ae. albopictus population from Indonesia . The number of alleles detected in our study of a larger number of field populations was significantly higher at all loci, except for AaelB51. Overall, the observed level of polymorphism and the small number of deviations from Hardy-Weinberg equilibrium across markers and samples suggest that this set of microsatellite markers is suitable for population genetic studies of Ae. albopictus in Africa.
MtDNA markers have been extensively used to assess the genetic diversity of Ae. albopictus populations across most of its geographic range. The level of polymorphism found here within COI and ND5 sequences was low (four haplotypes for each gene) and consistent with previous studies of populations sampled in newly invaded areas , , , , in which the number of haplotypes per country never exceeded five, whatever the mtDNA marker used (ND5, COI or Cytb). The discrepancy in the levels of polymorphism obtained with the microsatellite and mtDNA markers was partly attributable to differences in mutation rates between maker classes. Microsatellites exhibit extreme intraspecific variability, Mendelian inheritance and high mutation rates, whereas mtDNA coding regions are relatively stable, haploidic and maternally inherited . In addition, the mosquito samples we used for sequencing were small (10 individuals per sample) compared to those used for microsatellite genotyping. Even with larger numbers of sequenced specimens (27 and 40, in Garoua-Boulai and Buea respectively), the haplotype diversity indexes remained similar to those observed in other populations. This suggests that the low mtDNA diversity observed in Ae. albopictus in Cameroon was not related to the sample size. This low overall level of genetic diversity is consistent with recent introduction of Ae. albopictus in Cameroon.
The pattern of pairwise genetic differentiation estimated from the microsatellite data indicated a significant level of genetic structure among Ae. albopictus populations in Cameroon. No obvious relation between genetic and geographic distances was found, suggesting that the genetic structure has been shaped by additional biotic or abiotic factors. The Bayesian cluster analysis (i.e. existence of 2 genetic clusters co-occuring in most of cameroonian samples) supports a multiple introduction hypothesis. Recent studies have shown the important role of multiple introductions in successful colonization by invasive species  and in their genetic structure . The introduction phase is usually subject to severe bottlenecks that reduce the genetic variability of founding populations to levels incompatible with an expansion phase. It is now thought that multiple introductions may help to reach sufficient additive genetic variability for efficient expansion and adaptation to novel environments. The Buea population was highly differentiated from other geographic samples, possibly reflecting distinct introduction events with different geographic origins. Indeed, Buea is geographically and culturally close to Nigeria, where the first waves of Ae. albopictus colonization in Africa were recorded in the 1990s . Intensive trade between Buea and Nigeria may have favored passive dispersion of Ae. albopictus through accidental transportation of eggs, larvae or adults. Further analyses including samples from Nigeria would help to elucidate the relationships between Nigerian and Cameroonian populations of Ae. albopictus. Besides the effect of multiple introductions on the genetic structure of Ae. albopictus in Cameroon, the role of genetic drift must be considered. Indeed, the genetic variation observed in local populations of Ae. albopictus in Brazil  and in Reunion Island  has been attributed to genetic drift associated with restricted gene flow in recently established populations.
Analyses of COI sequences revealed that Cameroonian Ae. albopictus were related to specimens originating from tropical rather than temperate or subtropical areas. Unfortunately, we were unable to determine the geographic origin more precisely. Several studies have shown that Ae. albopictus populations that invaded Brazil ,  were genetically related to tropical rather than temperate Asian populations. These results, together with our findings, suggest that tropical lineages possess genetic and physiological traits that allow them to colonize new tropical areas more efficiently than temperate lineages. Among these putative physiological differences between temperate and tropical populations, Hawley et al.  has reported distinct winter diapause behavior in laboratory conditions.
Despite the low genetic variation across the sampled Ae. albopictus populations, this study suggests that Ae. albopictus invasion of Cameroon involved multiple introductions from tropical sources. Further studies with more extensive sampling and a larger number of genetic markers, including other nuclear loci, are needed to confirm this hypothesis and to determine the precise geographic origin of Cameroonian Ae. albopictus.
Microsatellite based Bayesian cluster analysis using STRUCTURE . (A). Estimates of Δ K, based on the second order rate of change of the likelihood function with respect to K , to determine the most likely number of clusters (K) in the data set. (B). Graphical representation of the data set for the most likely K = 2, where each color 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) Abong-Mbang, 2) Bertoua, 3) Garoua-Boulai, 4) Ayos, 5) Bafia, 6) Yaoundé, 7) Buea, 8) Douala, 9) Pouma, 10) Bamenda, 11) Banganté and 12) Bafoussam.
Estimates of null allele frequencies per locus and per sample of Ae. albopictus. Bolded: loci showing significant (P<0.05) heterozygote deficiency after correction for multiple testing.
We thank Elysée Nchoutpouen and Philippe Bousses for their assistance in mosquito sampling.
Conceived and designed the experiments: BK CP. Performed the experiments: BK CB. Analyzed the data: BK FS CP. Contributed reagents/materials/analysis tools: BK CB CP. Wrote the paper: BK FN FS DF CP.
- 1. Kow CY, Koon LL, Pang FY (2001) Detection of Dengue Viruses in Field Caught Male Aedes aegypti and Aedes albopictus (Diptera: Culicidae) in Singapore by Type-Specific PCR. J Med Entomol 38: 475–479.
- 2. Reiter P, Fontenille D, Paupy C (2006) Aedes albopictus as an epidemic vector of chikungunya virus: another emerging problem? Lancet Infect Dis 6: 463–464.
- 3. Paupy C, Ollomo B, Kamgang B, Moutailler S, Rousset D, et al. (2010) Role of Ae. albopictus versus Ae. aegypti in dengue and chikungunya emergences in Central Africa: Laboratory and field evidence. Vector Borne Zoonotic Dis 10: 259–266.
- 4. Smith CEG (1956) The history of dengue in tropical Asia and its probable relationship to the mosquitoes Aedes aegypti. Journal of Tropical Medicine and Hygiene 59: 3–11.
- 5. Gratz NG (2004) Critical review of the vector status of Aedes albopictus. Med Vet Entomol 18: 215–227.
- 6. Reiter P (1998) Aedes albopictus and the world trade in used tires, 1988–1995: the shape of things to come? J Am Mosq Control Assoc 14: 83–94.
- 7. Cornel AJ, Hunt RH (1991) Aedes albopictus in Africa? First records of live specimens in imported tires in Cape Town. J Am Mosq Control Assoc 7: 107–108.
- 8. Savage HM, Ezike VI, Nwankwo ACN, Spiegel R, Miller BR (1992) First record of breeding populations of Aedes albopictus in continental Africa: implications for arboviral transmission. J Am Mosq Control Assoc 8: 101–103.
- 9. Fontenille D, Toto JC (2001) Aedes (Stegomyia) albopictus (Skuse), a potential new dengue vector in Southern Cameroon. Emerg Infect Dis 7: 1066–1067.
- 10. Simard F, Nchoutpouen E, Toto JC, Fontenille D (2005) Geographic distribution and breeding site preference of Aedes albopictus and Aedes aegypti (Diptera: Culicidae) in Cameroon, Central Africa. J Med Entomol 42: 726–731.
- 11. Kamgang B, Happi JY, Boisier P, Njiokou F, Hervé JP, et al. (2010) Geographic and ecological distribution of the dengue and chikungunya virus vectors Aedes aegypti and Aedes albopictus in three major Cameroonian towns. Med Vet Entomol 24: 132–141.
- 12. Peyrefitte CN, Rousset D, Pastorino BA, Pouillot R, Bessaud M, et al. (2007) Chikungunya virus, Cameroon, 2006. Emerging Infectious Diseases 13: 768–771.
- 13. Leroy EM, Nkoghe D, Ollomo B, Nze-Nkogue C, Becquart P, et al. (2009) Concurrent chikungunya and dengue virus infections during simultaneous outbreaks, Gabon. Emerg Infect Dis 15: 591–593.
- 14. Dalla Pozza G, Majori G (1992) First record of Aedes albopictus establishment in Italy. J Am Mosq Control Assoc 8: 318–320.
- 15. Hawley AH (1988) The biology of Aedes albopictus. J Am Mosq Control Assoc 4: 2–39.
- 16. Scholte EJ, Schaffner F (2007) Waiting the tiger: establishment and spread of the Aedes albopictus mosquito in Europe. In: Takken W, Knols BGJ, editors. Emerging Pests and Vector-Borne Diseases in Europe. Wageningen Academic Publishers, Wageningen, The Netherlands.
- 17. Delatte H, Desvars A, Bouétard A, Bord S, Gimonneau G, et al. (2010) Blood-feeding behavior of Aedes albopictus, a vector of Chikungunya on La Réunion. Vector Borne Zoonotic Dis 10: 249–258.
- 18. Gubler DJ (2003) Aedes albopictus in Africa. The Lancet Infect Dis 3: 751–752.
- 19. Kambhampati S, Black WC, Rai KS (1991) Geographic origin of the US and Brazilian Aedes albopictus inferred from allozyme analysis. Heredity 67: 85–93.
- 20. Urbanelli S, Bellini R, Carrieri M, Sallicandro P, Celli G (2000) Population structure of Aedes albopictus (Skuse): the mosquito which is colonizing Mediterranean countries. Heredity 84: 331–337.
- 21. Birungi J, Munstermann LE (2002) Genetic structure of Aedes albopictus (Diptera: Culicidae) populations based on mitochondrial ND5 sequences: evidence for an independent invasion into Brazil and United States. Ann Entomol Soc Am 95: 125–132.
- 22. Mousson L, Dauga C, Garrigues T, Schaffner F, Vazeille M, et al. (2005) Phylogeography of Aedes (Stegomyia) aegypti (L.) and Aedes (Stegomyia) albopictus (Skuse) (Diptera: Culicidae) based on mitochondrial DNA variations. Genet Res 85: 1–11.
- 23. Gubler DJ, Rosen L (1976) Variation among geographic strains of Aedes albopictus in susceptibility to infection with dengue viruses. Am J Trop Med Hyg 25: 318–325.
- 24. Tesh RB, Gubler DJ, Rosen L (1976) Variation among geographic strains of Ae. albopictus in susceptibility to infection with chikungunya virus. Am J Trop Med Hyg 33: 176–181.
- 25. Lourenço-de-Oliveira R, Vazeille M, Bispo-de-Filippis AM, Failloux AB (2003) Large genetic differentiation and low variation in vector competence for dengue and yellow fever viruses of Aedes albopictus from Brazil, the United States, and the Cayman islands. Am J Trop Med Hyg 69: 105–114.
- 26. Morlais I, Ponçon N, Simard F, Cohuet A, Fontenille D (2004) Intraspecific nucleotide variation in An. gambiae promises new insights into the biology of malaria vectors. Am J Trop Med Hyg 71: 795–802.
- 27. Porreta D, Gargani M, Bellini R, Calviti M, Urbanelli S (2006) Isolation of microsatellite markers in the tiger mosquito Aedes albopictus (Skuse). Mol Ecol Notes 6: 880–881.
- 28. El Mousadik A, Petit RJ (1996) High level of genetic differentiation for allelic richness among populations of the argan tree (Argania spinosa L. Skeels) endemic of Morocco. Theor Appl Genet 92: 832–839.
- 29. Nei M (1987) Molecular evolutionary genetics. Columbia University Press, New York.
- 30. Goudet J (1995) FSTAT (version 1.2): a computer software to calculate F-statistics. J Heredity 86: 485–486.
- 31. Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F (2004) GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Laboratoire Génome, Populations, Interactions, CNRS UMR 5000, Université de Montpellier II, Montpellier (France), http://www.genetix.univmontp2.fr/genetix/genetix.htm.
- 32. Raymond M, Rousset F (1995) GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered 86: 248–249.
- 33. Weir BS, Cockerham CC (1984) Estimating F-statistics for the analysis of population structure. Evolution 38: 1358–1370.
- 34. Holm S (1979) A simple sequentially rejective multiple test procedure. Scand J Stat 6: 65–70.
- 35. Van Oosterhout C, Hutchinson WF, Wills DP, Shipley P (2004) MICROCHECKER: Software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol 4: 535–538.
- 36. Rousset F (1997) Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics 145: 1219–1228.
- 37. Liu K, Muse SV (2005) POWERMARKER: Integrated analysis environment for genetic marker data. Bioinformatics 21: 2128–2129.
- 38. Nei M (1983) Genetic polymorphism and the role of mutation in evolution. In: Nei M, Koehn R, editors. Evolution of genes and proteins. Sunderland, MA: Sinauer Associates.
- 39. Felsenstein J (1989) PHYLIP-phylogeny inference package (version 3.2). Cladistics 5: 164–166.
- 40. Page RDM (1996) TREEVIEW: An application to display phylogenetic trees on personal computers. Comput Appl Biosci 12: 357–358.
- 41. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959.
- 42. 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.
- 43. Bonacum J, DeSalle R, O'Grady P, Olivera D, Wintermute J, et al. (2001) New nuclear and mitochondrial primers for systematics and comparative genomics in Drosophilidae. Dros Inf Serv 84: 201–204.
- 44. Thompson JD, Higgins D, Gibson TJ (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22: 4673–4680.
- 45. Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R (2003) DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 19: 2496–2497.
- 46. Tajima F (1989) The effect of change in population size on DNA polymorphism. Genetics 123: 597–601.
- 47. Fu YX, Li WH (1993) Statistical tests of neutrality mutations. Genetics 133: 693–709.
- 48. Fu YX (1997) Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147: 915–925.
- 49. Nylander JAA (2004) MrModeltest v2. Evolutionary BiologyCentre, Uppsala University. Program distributed by the author.
- 50. Ronquist F, Huelsenbeck JP (2003) MRBAYES v3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574.
- 51. Usmani-Brown S, Cohnstaedt L, Munstermann EL (2009) Population Genetics of Aedes albopictus (Diptera: Culicidae) Invading Populations, Using Mitochondrial nicotinamide Adenine Dinucleotide Dehydrogenase Subunit 5 Sequences. Ann Entomol Soc Am 102: 144–150.
- 52. Maia RT, Scarpassa VM, Maciel-Litaiff LH, Tadei WP (2009) Reduced levels of genetic variation in Aedes albopictus (Diptera: Culicidae) from Manaus, Amazonas State, Brazil, based on analysis of the mitochondrial DNA ND5 gene. Gen Mol Res 8: 998–1007.
- 53. Mas-Coma S, Bargues MD (2009) Populations, hybrids and the systematic concepts of species and subspecies in Chagas disease triatomine vectors inferred from nuclear ribosomal and mitochondrial DNA. Acta Tropica 110: 112–136.
- 54. Kolbe JJ, Glor RE, Schettino LR, Lara AC, Larson A, et al. (2004) Genetic variation increases during biological invasion by a Cuban lizard. Nature 431: 177–181.
- 55. Zalewski A, Michalska-Parda A, Bartoszewicz M, Kozakiewicz M, Brzeziński M (2010) Multiple introductions determine the genetic structure of an invasive species population: American mink Neovison vison in Poland. Biol Conserv 1436: 1355–1363.
- 56. Ayres CF, Romao TP, Melo-Santos MA, Furtado AF (2002) Genetic diversity in Brazilian populations of Aedes albopictus. Mem Inst Oswaldo Cruz 97: 871–875.
- 57. Paupy C, Girod R, Salvan M, Rodhain F, Failloux AB (2001) Population structure of Aedes albopictus from La Réunion Island (Indian Ocean) with respect to susceptibility to a dengue virus. Heredity 87: 273–283.
- 58. Hawley WA, Reiter P, Copeland RS, Pumpuni CB, Craig GB Jr (1987) Aedes albopictus in North America: probable introduction in used tires from northern Asia. Science 236: 1114–1116.