- Open Access
Molecular markers for analyses of intraspecific genetic diversity in the Asian Tiger mosquito, Aedes albopictus
Parasites & Vectorsvolume 8, Article number: 188 (2015)
The dramatic worldwide expansion of Aedes albopictus (the Asian tiger mosquito) and its vector competence for numerous arboviruses represent a growing threat to public health security. Molecular markers are crucially needed for tracking the rapid spread of this mosquito and to obtain a deeper knowledge of population structure. This is a fundamental requirement for the development of strict monitoring protocols and for the improvement of sustainable control measures.
Wild population samples from putative source areas and from newly colonised regions were analysed for variability at the ribosomal DNA internal transcribed spacer 2 (ITS2). Moreover, a new set of 23 microsatellite markers (SSR) was developed. Sixteen of these SSRs were tested in an ancestral (Thailand) and two adventive Italian populations.
Seventy-six ITS2 sequences representing 52 unique haplotypes were identified, and AMOVA indicated that most of their variation occurred within individuals (74.36%), while only about 8% was detected among populations. Spatial analyses of molecular variance revealed that haplotype genetic similarity was not related to the geographic proximity of populations and the haplotype phylogeny clearly indicated that highly related sequences were distributed across populations from different geographical regions. The SSR markers displayed a high level of polymorphism both in the ancestral and in adventive populations, and F ST estimates suggested the absence of great differentiation. The ancestral nature of the Thai population was corroborated by its higher level of variability.
The two types of genetic markers here implemented revealed the distribution of genetic diversity within and between populations and provide clues on the dispersion dynamics of this species. It appears that the diffusion of this mosquito does not conform to a progressive expansion from the native Asian source area, but to a relatively recent and chaotic propagule distribution mediated by human activities. Under this scenario, multiple introductions and admixture events probably play an important role in maintaining the genetic diversity and in avoiding bottleneck effects. The polymorphic SSR markers here implemented will provide an important tool for reconstructing the routes of invasion followed by this mosquito.
The Asian tiger mosquito, Aedes (Stegomyia) albopictus is a vector of numerous arboviruses including Dengue, Chikungunya, and a variety of epizootic viruses [1-3]. It is considered to be the most invasive mosquito species in the world , and it is listed among the top hundred most dangerous invasive species (http://www.issg.org/database/species/search.asp?st=100ss).
Originally a native of the tropical forests of Southeast Asia, it has spread during the last 30–40 years to both tropical and temperate regions worldwide [4,5]. Global warming can result in wetter and warmer conditions that favour the spread of this mosquito and its aquatic larvae [6,7]. The dramatic expansion of this aggressive daytime biting mosquito has increased public health concerns for arbovirus-related disease outbreaks [8-11]. Some biological traits such as feeding preferences , winter diapause , and vector competence  of Ae. albopictus vary according to its geographic origin. Given that these traits may be genetically determined and may influence pathogen transmission, it is of paramount importance to define the geographic origin of invading populations.
Furthermore, as the threat represented by this mosquito is growing due to the lack of sustainable control measures and of the progressive spread of insecticide-resistance [15,16], efficient surveillance methods may depend on the detailed knowledge of the population structure and dynamics of this vector [10,17]. Moreover, the potential presence of cryptic subgroups, that may be unaffected by control methods, should not be under-estimated. Indeed, great variation in terms of genome size has been previously identified in strains from different regions [18-21]. Variations in the rDNA cistron have also been reported within and among populations of this mosquito . In this context, the Internal Transcribed Spacer 2 (ITS2) of nuclear ribosomal DNA has been shown to be a useful marker, because of its power in detecting differentiation, not only among cryptic species [23-25], but also among conspecific populations [26-28]. As ITS2 is one of the less functionally restricted regions of the rDNA, it can accumulate mutations within isolated populations relatively quickly and can thus be an indicator of genetic discontinuity between populations [29,30]. The genetic diversity and the phylogeographic relationships among populations have been extensively analysed using mitochondrial DNA (mtDNA) markers [31-35]. The data from populations collected in various regions of the world were consistent in displaying the presence of weak sequence variation at mtDNA loci, probably due to the recent spread of the species, but also due to the low mutation rate of these markers. Indeed, when a set of 13 microsatellite markers were used to analyse the invasion of the Australian region , they appeared to be more informative, compared to mtDNA markers, in resolving the population structure and the invasion dynamics. Moreover, considering that Ae. albopictus populations are infected by the symbiont Wolbachia , whose presence can impact mtDNA diversity within host populations, the utility of mtDNA as a neutral marker is reduced .
In the present study we used the ITS2 region to address the presence of intraspecific differentiation (genetic barriers) among groups of populations from different ancestral and recently colonized geographic areas. Moreover, we developed additional novel microsatellite markers for this mosquito species. We provide evidence of their usefulness for the analysis of the genetic composition of populations.
Nine wild population samples of Ae. albopictus were collected across the species range (Table 1). Two samples were from Thailand (Ban Rai, Lampang Province in the North and Phato, Chumphon Province in the South); two from Réunion (St. Denis in the North and St. Pierre in the South-West); and five from Northern Italy (Brescia, Castellanza, Pavia, Modena, Cesena). The Thai samples are from a region that pertains to the putative home range of the species, while Réunion and the Italian samples are representative of old (17th and 18th century)  and recently colonized (1990)  areas, respectively. For each of the samples, eggs were collected using ovitraps and the emerging adults were reared in the insectary at 27°C, with 60-80% relative humidity and a 12:12 h (L:D) photoperiod. The identities of all the samples were confirmed using the morphological keys of Rueda . The nine population samples were screened for sequence variability of the ITS2 marker.
The Rimini laboratory strain (F35) was used for the construction of the microsatellite-enriched libraries. This strain was established from mosquitoes collected in Rimini (Emilia Romagna) at the Centro Agricoltura Ambiente ‘G. Nicoli’ (Crevalcore, Italy). Three of the nine wild samples previously described, Ban Rai, Cesena and Brescia, were also screened to assess the variability of the newly developed microsatellites markers.
ITS2 marker characterization
DNA was extracted according to Baruffi et al.  from individual females of the population samples. The extracted DNA was resuspended in TE (10 mM Tris–HCl, pH 8, 1 mM EDTA) and its concentration was quantified using a Nanodrop ND-1000 spectrophotometer (Nano-drop Technologies Inc., Wilmington, DE, USA). The ribosomal DNA ITS2 region was amplified from the DNA samples using the universal primers 5.8Sf 5′-atcactcggctcgtggatcg-3′ and 28Sr 5′-atgcttaaatttagggggtagtc-3′ . These primers anneal to highly conserved sequences in the 5.8S and 28S rDNA genes flanking the ITS2 region. Reactions were performed in a volume of 25 μl with approximately 3 ng DNA, 1.5 mM MgCl2, Reaction buffer (10 mM Tris–HCl pH 9, 50 mM KCl), 200 µM dNTP, 10 pmol of each primer and 1 unit Taq DNA polymerase (Invitrogen). Amplifications were performed using an Eppendorf MasterCycler with the following cycle conditions: an initial denaturing step at 94°C for 2 min; twenty-five cycles of denaturation at 94°C for 20 s, annealing at 54°C for 15 s, extension at 72°C for 25 s, followed by a final extension at 72°C for 5 min.
The PCR products were separated on 1.5% agarose gels, stained with ethidium bromide and photographed under ultraviolet light. Amplification products were ligated into the TOPO 2.1 vector and transformed into competent One Shot cells using the TOPO TA cloning kit (Invitrogen). DNA inserts were sequenced on both strands.
ITS2 data analysis
The boundaries of the ITS2 region were determined by comparison with the available Ae. albopictus sequences in GenBank. The ITS2 sequences were aligned using MUSCLE 3.8 . Nucleotide diversity and gene diversity parameters were determined using Arlequin 22.214.171.124 . Geographic structure in the ITS2 data set was investigated by spatial analyses of molecular variance implemented in SAMOVA 1.0 . This approach identifies groups of populations that are genetically and geographically homogeneous and maximally differentiated from each other. The method requires the a priori definition of the number of groups (K) of populations that exist and generates F-statistics (F SC, F ST and F CT) using an AMOVA approach. Different numbers of groups (K) were tested, and a simulated annealing procedure permitted the identification of the composition of each of the K groups that maximizes the F CT index (proportion of total genetic variance due to differences between groups). The program was run for two to eight groups (K = 2 to K = 8) each time with the simulated annealing process repeated 100 times, starting each time with a different partition of the population samples into the K groups. A phylogenetic analysis was also performed including the unique ITS2 sequences or haplotypes identified in the nine population samples by Arlequin 126.96.36.199 . The haplotype sequences together with an Ae. flavopictus ITS2 outgroup sequence (GenBank accession number: AF353559.1) were aligned using MUSCLE 3.8 . The Kimura-2-parameter and Γ distribution model indicated by the Akaike Information Criterion (AIC = 3064) and the Bayesian Information Criterion (BIC = 3899), was used to perform Maximum Likelihood analyses in MEGA 5.2.2 . The reliability of the resulting haplotype tree was determined by 1000 bootstrap replications. The mid-point rooted tree was drawn using FigTree v1.4 (http://tree.bio.ed.ac.uk/software/figtree/).
Microsatellite isolation and characterization
Total genomic DNA was extracted individually from 64 mosquitoes (40 males and 24 females) from the Rimini strain using the method reported by Baruffi et al. . The DNA was subsequently pooled and submitted to Genetic Identification Services, GIS (Chatsworth, CA, USA, http://www.genetic-id-services.com) for library construction. The library was enriched for four different di- and tri-nucleotide microsatellite motifs, namely (CA)n, (GA)n, (AAC)n, and (ATG)n. Insert DNA from individual clones was amplified by PCR following GIS guidelines. The PCR products were purified using QIAquick columns (Qiagen) and sequenced using the DYEnamic ET Terminator Cycle Sequencing Kit (Amersham Bioscience) with the universal M13 primers (M13f 5′-aggaaacagctatgaccatg-3′ and M13r 5′-acgacgttgtaaaacgacgg-3′). Clones with inserts less than 200 bp were not considered for sequencing. Sequencing products were resolved on the Applied Biosystems model 377 DNA Sequencer (Applied Biosystems). The resulting sequences were screened for the presence of microsatellite motifs using WebSat . All the clone sequences were checked for redundancy using BLASTn .
Microsatellite sequences that shared one or both the 5′/3′ flanking regions in common were discarded to exclude potential multi-locus microsatellites families. The sequences that were found to be unique within the library were subjected to BLASTn analysis against the nr database to identify eventual hits with previously described sequences. Primer pairs were designed for the flanking regions of the identified SSR sequences using Primer 3 . Preliminary PCR screenings were performed on 20 individuals (10 males and 10 females) from the Rimini strain to validate the primer efficiencies in single SSR locus amplifications. DNA amplifications were performed on a Mastercycler gradient (Eppendorf). The PCR reaction mixture consisted of 50 ng genomic DNA, 1x PCR buffer, 1.5 mM MgCl2, 270 µM dNTPs (Invitrogen), 1 U Taq polymerase (Invitrogen), and 10 μM of each primer, one of which was 5′ labelled with a fluorescent dye, in a final volume of 15 μl. PCR cycling conditions were 94°C for 3 min, 30 cycles of 94°C for 30 s, 57-60°C for 30 s, and 72°C for 30 s, followed by a final extension step of 72°C for 10 minutes. Aliquots of PCR products were separated by electrophoresis on 2% agarose gels stained with ethidium bromide, and visualized under UV light. Each PCR product was then diluted 1:10 in ddH2O water and 1 μl of this dilution was added to 9 μl of a mixture of deionized formamide and GeneScan-500 ROX size standard (Applied Biosystems). After denaturation for 4 min at 94°C the fragments were resolved on an ABI PRISM 310 Genetic Analyzer (Applied Biosystems) with the GENESCAN software package (Applied Biosystems). To avoid genotyping errors, the program TANDEM  was used to automate binning of the microsatellite allele length variants.
Microsatellite data analysis
The parameters of intra- and inter-population genetic diversity, such as observed allele size, number of actual alleles (na), number of effective alleles (ne), number of private alleles (np) and AMOVA were computed using GenALEx 6.5 . Observed and expected heterozygosity (HO and HE, respectively), gene diversity (HS) and F ST values were computed using MICROSATELLITE ANALYSER (MSA) V.4.05 . Linkage disequilibrium, the frequency of null alleles (An) and the deviations from Hardy-Weinberg expectations were computed using GENEPOP version 4.2 [54,55]. For loci with fewer than five alleles, an exact test of Hardy-Weinberg proportions was performed. For loci with five alleles or more, an unbiased estimate of the exact probability was obtained using the Markov chain method of Guo & Thompson  for each combination of locus and population. The allelic Polymorphic Information Content (PIC) was derived using CERVUS .
ITS2 sequence variability
ITS2 sequences were obtained from an average of three individuals from each of the nine population samples collected in Thailand, Réunion and Italy (Table 2). In addition, to determine whether multiple copies of the ITS2 region varied within individuals, up to four clones were sequenced from each of these individuals. The primers amplified a product of 502–588 bp in length that included part of 5.8S and 28S sequences. The ITS2 sequence itself ranged from 322 to 408 bp in length with an average GC composition of 51%. The 76 ITS2 sequences, obtained from a total of 26 individuals, represent 52 unique sequence-types or haplotypes (Table 2). The gene diversity estimates were very high across the analysed sequences and ranged from 0.893 ± 0.111 in St. Denis (Réunion) and Castellanza (Italy) to 1.000 ± 0.052 in Brescia (Italy). Nucleotide diversity (π), was very low, ranging from 0.005 ± 0.004 in St. Pierre (Réunion) to 0.0191 ± 0.012 in Cesena (Italy). These low nucleotide diversity estimates are expected, as in the 415 bp MUSCLE ITS2 alignment only 44 sites (~10%) contained substitutions.
An estimate of variability distribution in and among the nine tested population samples is reported in Table 3. AMOVA indicated that most of the variation occurs within individuals (74.36%) while only about 8% of total variation is detected among populations. The inter-individual variation within samples accounted for 17.51% (P < 0.05) of the total variance. When the three main geographical population groups (Thailand, Réunion and Italy) are accounted for (Table 4), similar and relatively low percentages of variation were observed among groups (8.12%) and among populations within groups (7.09%), indicating that the genetic similarity is not related to the geographic proximity of populations. Indeed, according to the SAMOVA analysis (Table 5) the best partitioning of genetic diversity was obtained when the nine samples were grouped into seven groups: F CT =0.186, P = 0.011. Here the genetic diversity is not only among geographical regions, as even the two ancestral Thai populations are separated into two groups, while the North-Italian populations are separated into three groups and the two Réunion populations into two groups.
ITS2 haplotype phylogeny and population distribution
The phylogenetic analysis based on the 52 ITS2 haplotype sequences identified in the nine geographic populations resulted in a Maximum Likelihood tree with a log likelihood of −1428.3 (Figure 1). The ITS2 sequences form a compact cluster with very few nodes supported by bootstrap values greater than 50%. Highly related sequences are distributed across populations from different geographical regions; some evidence of clustering of tightly related sequences is present in Ban Rai, Phato (Thailand) and to a lesser degree in St. Pierre (Réunion).
The majority (47of 52) of the haplotypes are present in only one population. Only three haplotypes are shared by two populations (haplotype 18 in Pavia and St. Denis, haplotype 21 in Pavia and St. Pierre, and haplotype 39 in Cesena and Modena). Haplotype 23 is shared by three Italian populations (Pavia, Brescia and Modena), and haplotype 13 is present in four populations (Phato, Modena, St. Denis and St. Pierre).
SSR: Isolation and characterization
Out of a total of 109 sequenced clones obtained from the Ae. albopictus library, 92 contained tandem DNA repeats with at least 5 repeat units. BLASTn analyses showed that one or both flanking regions were shared in 42 sequences and they were thus excluded from the subsequent characterization. Within the remaining 50 sequenced clones, five harboured more than one tandem DNA repeat, resulting in the identification of 55 microsatellite motifs. Among these, trinucleotide motifs were the most abundant (28), followed by dinucleotide (22), mononucleotide (3), heptanucleotide (1) and hendecanucleotide motifs (1). Among the trinucleotides motifs, CAA was highly abundant (17) followed by TGA (10). Among the dinucleotide motifs, TG was predominant (14) followed by CT (8). Thirty-seven of these microsatellite motifs contained perfect repeat arrays. Twelve microsatellites contained imperfect repeats and the other six microsatellites contained compound repeats.
Forty-five microsatellite sequences were considered for further validation, as the others lacked sufficient flanking regions to design primers. After the PCR screening on 20 individuals from the Rimini laboratory strain, 23 SSR sequences were validated as loci (Table 6). The nucleotide sequences of these validated SSR loci were submitted to GenBank [Genbank: KP859591-KP859613].
Evaluation of microsatellite polymorphism in ancestral and adventive populations
On the basis of a preliminary screening performed on the Rimini strain, 16 (Aealbmic1 to Aealbmic16) of the 23 validated SSR loci (Table 6) were used to analyse the degree of polymorphism in 79 mosquitoes from three populations: Ban Rai (Thailand), Brescia and Cesena (Italy) (Table 7). The average Polymorphic Information Content (PIC) estimated across the 16 loci was 0.56 ± 0.29 in Ban Rai, 0.50 ± 0.26 in Cesena and 0.47 ± 0.21 in Brescia, suggesting that these loci are informative for population analyses. Heterozygous individuals were present in both sexes for all 16 loci, excluding any condition of sex linkage. Tests for Hardy-Weinberg equilibrium (HWE), after sequential Bonferroni correction , revealed that five loci did not meet the expectations of this model (Table 7): three (Aealbmic4, 11, 13) in the Ban Rai sample, one in Cesena (Aealbmic3) and two in Brescia (Aealbmic13, 16). However, the locus/population combinations that were not in HWE were not concentrated at any locus or in any population. As no evidence of linkage disequilibrium between loci was assessed, these 16 loci can be considered to be independent. In the three populations most of the loci display relatively high levels of gene diversity (Hs). The frequency of null alleles across the loci was generally relatively low apart from Aealbmic15 with an estimate of 0.98 in Cesena. Two alleles were detected in Cesena at the Aealbmic15 locus, but one was very rare. In the other populations this locus was monomorphic for the common allele. Aealbmic1 is monomorphic in Ban Rai and Cesena, while in Brescia two alleles were present. Across the loci a higher level of polymorphism was observed in Ban Rai, with a mean number of alleles per locus of 6.19 (Table 7). Private alleles (np) were detected in all populations, with the highest average value in the Ban Rai sample (2.25). Analysis of molecular variance (AMOVA) across the three populations indicates that the greatest portion of variance, 76.62%, is found within individuals, while 16.80% is among individuals and only the remaining 6.57% is among populations. The FST values between each pair of populations are statistically significant and indicate that the degree of differentiation between the two Northern Italian samples, Cesena and Brescia (FST = 0.061), is not very different from those between these two Italian samples and the Ban Rai sample (FST = 0.072 in both comparisons).
In the absence of a fully annotated genome, molecular markers for Ae. albopictus are crucially needed to track the rapid spread of this mosquito that is causing increasing concern due to its status as a vector of several arboviruses of public health importance [3,9]. The two types of genetic markers here implemented give a picture of the distribution of genetic diversity within and between populations providing clues on the dispersion dynamics of this species. ITS2 sequence variation has been analysed in mosquitoes from the supposed native range, Thailand, from an old colonised area, Réunion, and from a newly invaded area, Italy. It appears that substantial variability, both in terms of haplotype number and gene diversity, is present in the native as well in the two areas of introduction. The haplotype variability is mostly concentrated within individuals (74.36%) or within populations (84.78%), while only a relatively small proportion into geographical differences, both at single population (8.13%) and at geographic area level (8.12%). This implies that, according to SAMOVA, the haplotype variation, besides being high, has a distribution that is independent from geography. Different haplotypes are present in each individual of the nine samples, and the tree clearly demonstrates that they are similar and dispersed across native and adventive populations (Figure 1). Substantially concordant conclusions can be obtained with the newly developed SSR markers, which display a high level of polymorphism both in the ancestral (Ban Rai, Thailand) and adventive populations (North Italian samples), previously scored for ITS2 variability. Although the chromosomal location of these markers remains unknown, the assessed linkage equilibrium between them suggests that they are statistically independent and their variability patterns might reflect genome wide patterns across populations. Based on the population samples here analysed, the SSR markers were, as expected, more informative than ITS2 in revealing the slight genetic diversity between native and derived populations both in terms of variability and differentiation. About 7% of total variability (AMOVA) is represented by the differences between the three populations, while the highest variability was found within populations. This is an agreement with a previous study on mosquitoes collected in Réunion . The FST estimates indicate that the differentiation between the two Italian samples is slightly smaller (FST = 0.061) than that between the Italian and Thai samples (FST = 0.072), suggesting the absence of great differentiation between ancestral and derived populations. The ancestral nature of the Thai population is also corroborated by its high level of variability.
Taken together these genetic data raise questions on the way this mosquito is spreading and becoming established worldwide. Clearly, our data indicate that the diffusion of this mosquito does not conform to a progressive expansion away from the native Asian source area, but supports a relatively recent and chaotic propagule distribution mediated by human activities. Under this scenario, multiple introductions and admixture events probably play an important role in maintaining the genetic diversity and in avoiding bottleneck effects.
In this study we have exploited nuclear molecular markers that have allowed us to highlight the presence of high intraspecific variability and a low/moderate level of differentiation of Ae. albopictus populations collected in different eco-geographic areas. Most of the genetic variation was detected at the individual level and this contributed to the high genetic variation observed within populations. This genetic feature, also revealed using other markers [7,34,36,59] may have facilitated the characteristic ability of this invasive species to expand and adapt to novel environments.
The highly polymorphic SSR markers here implemented, together with the set previously described [36,59,60], represent an important tool for reconstructing the routes of invasion and for the identification of the origins of mosquito outbreaks.
Gratz NG. Critical review of the vector status of Aedes albopictus. Med Vet Entomol. 2004;18(3):215–27.
Paupy C, Delatte H, Bagny L, Corbel V, Fontenille D. Aedes albopictus, an arbovirus vector: from the darkness to the light. Microbes Infect Inst Pasteur. 2009;11(14–15):1177–85.
Lambrechts L, Scott TW, Gubler DJ. Consequences of the expanding global distribution of Aedes albopictus for dengue virus transmission. PLoS Negl Trop Dis. 2010;4(5):e646.
Medlock JM, Hansford KM, Schaffner F, Versteirt V, Hendrickx G, Zeller H, et al. A Review of the invasive mosquitoes in Europe: ecology, public health risks, and control options. Vector Borne Zoonotic Dis. 2012;12(6):435–47.
Enserink M. A mosquito goes global. Science. 2008;320(5878):864–6.
Erickson RA, Hayhoe K, Presley SM, Allen LJS, Long KR, Cox SB. Potential impacts of climate change on the ecology of dengue and its mosquito vector the Asian tiger mosquito (Aedes albopictus). Environ Res Lett. 2012;7:034003.
Medley KA, Jenkins DG, Hoffman EA. Human-aided and natural dispersal drive gene flow across the range of an invasive mosquito. Mol Ecol. 2015;24(2):284–95.
Tatem AJ, Hay SI, Rogers DJ. Global traffic and disease vector dispersal. Proc Natl Acad Sci. 2006;103(16):6242–7.
Benedict MQ, Levine RS, Hawley WA, Lounibos LP. Spread of the tiger: global risk of invasion by the mosquito Aedes albopictus. Vector Borne Zoonotic Dis. 2007;7(1):76–85.
Gasperi G, Bellini R, Malacrida AR, Crisanti A, Dottori M, Aksoy S. A new threat looming over the Mediterranean basin: emergence of viral diseases transmitted by Aedes albopictus mosquitoes. PLoS Negl Trop Dis. 2012;6(9):e1836.
Schaffner F, Medlock JM, Van Bortel W. Public health significance of invasive mosquitoes in Europe. Clin Microbiol Infect. 2013;19(8):685–92.
Delatte H, Desvars A, Bouétard A, Bord S, Gimonneau G, Vourc’h G, et al. Blood-feeding behavior of Aedes albopictus, a vector of Chikungunya on La Réunion. Vector Borne Zoonotic Dis. 2010;10(3):249–58.
Poelchau MF, Reynolds JA, Elsik CG, Denlinger DL, Armbruster PA. Deep sequencing reveals complex mechanisms of diapause preparation in the invasive mosquito, Aedes albopictus. Proc Biol Sci. 2013;280(1759):20130143.
Lourenço de Oliveira R, Vazeille M, de Filippis AM, Failloux AB. 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. 2003;69(1):105–14.
Tantely ML, Tortosa P, Alout H, Berticat C, Berthomieu A, Rutee A, et al. Insecticide resistance in Culex pipiens quinquefasciatus and Aedes albopictus mosquitoes from La Réunion Island. Insect Biochem Mol Biol. 2010;40(4):317–24.
Kawada H, Maekawa Y, Abe M, Ohashi K, Ohba S, Takagi M. Spatial distribution and pyrethroid susceptibility of mosquito larvae collected from catch basins in parks in Nagasaki city, Nagasaki, Japan. Jpn J Infect Dis. 2010;63(1):19–24.
Bonizzoni M, Gasperi G, Chen X, James AA. The invasive mosquito species Aedes albopictus: current knowledge and future perspectives. Trends Parasitol. 2013;29(9):460–8.
Rao PN, Rai K. Inter and intraspecific variation in nuclear DNA content in Aedes mosquitoes. Heredity. 1987;59(2):253–8.
Kumar A, Rai KS. Intraspecific variation in nuclear DNA content among world populations of a mosquito, Aedes albopictus (Skuse). Theor Appl Genet. 1990;79(6):748–52.
McLain DK, Rai KS, Fraser MJ. Intraspecific and interspecific variation in the sequence and abundance of highly repeated DNA among mosquitoes of the Aedes albopictus subgroup. Heredity. 1987;58(3):373–81.
Black WC, Rai KS. Genome evolution in mosquitoes: intraspecific and interspecific variation in repetitive DNA amounts and organization. Genet Res. 1988;51(3):185–96.
Black WC, McLain DK, Rai KS. Patterns of variation in the rDNA cistron within and among World populations of a mosquito, Aedes albopictus (Skuse). Genetics. 1989;121(3):539–50.
Gomulski LM, Meiswinkel R, Delécolle J-C, Goffredo M, Gasperi G. Phylogeny of the subgenus Culicoides and related species in Italy, inferred from internal transcribed spacer 2 ribosomal DNA sequences. Med Vet Entomol. 2006;20(2):229–38.
Young I, Coleman AW. The advantages of the ITS2 region of the nuclear rDNA cistron for analysis of phylogenetic relationships of insects: a Drosophila example. Mol Phylogenet Evol. 2004;30(1):236–42.
Zomuanpuii R, Ringngheti L, Brindha S, Gurusubramanian G, Senthil KN. ITS2 characterization and Anopheles species identification of the subgenus Cellia. Acta Trop. 2013;125(3):309–19.
Wörheide G, Hooper JN, Degnan BM. Phylogeography of western Pacific Leucetta “chagosensis” (Porifera: Calcarea) from ribosomal DNA sequences: implications for population history and conservation of the Great Barrier Reef World Heritage Area (Australia). Mol Ecol. 2002;11(9):1753–68.
Behbahani A, Dutton TJ, Davies N, Townson H, Sinkins SP. Population differentiation and Wolbachia phylogeny in mosquitoes of the Aedes scutellaris group. Med Vet Entomol. 2005;19(1):66–71.
Giessler S, Wolinska J. Capturing the population structure of microparasites: using ITS-sequence data and a pooled DNA approach. Mol Ecol Resour. 2013;13(5):918–28.
Müller T, Philippi N, Dandekar T, Schultz J, Wolf M. Distinguishing species. RNA. 2007;13(9):1469–72.
Coleman AW. Is there a molecular key to the level of “biological species” in eukaryotes? A DNA guide. Mol Phylogenet Evol. 2009;50(1):197–203.
Birungi J, Munstermann LE. 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. 2002;95(1):125–32.
Mousson L, Dauga C, Garrigues T, Schaffner F, Vazeille M, Failloux AB. Phylogeography of Aedes (Stegomyia) aegypti (L.) and Aedes (Stegomyia) albopictus (Skuse) (Diptera: Culicidae) based on mitochondrial DNA variations. Genet Res. 2005;86(1):1–11.
Maia RT, Scarpassa VM, Maciel-Litaiff LH, Tadei WP. Reduced levels of genetic variation in Aedes albopictus (Diptera: Culicidae) from Manaus, Amazonas State, Brazil, based on analysis of the mitochondrial DNA ND5 gene. Genet Mol Res. 2009;8(3):998–1007.
Ashfaq M, Hebert PD, Mirza JH, Khan AM, Zafar Y, Mirza MS. Analyzing mosquito (Diptera: culicidae) diversity in Pakistan by DNA barcoding. PLoS One. 2014;9(5):e97268.
Kamgang B, Brengues C, Fontenille D, Njiokou F, Simard F, Paupy C. Genetic structure of the tiger mosquito, Aedes albopictus, in Cameroon (Central Africa). PLoS One. 2011;6(5):e20257.
Beebe NW, Ambrose L, Hill LA, Davis JB, Hapgood G, Cooper RD, et al. Tracing the tiger: population genetics provides valuable insights into the Aedes (Stegomyia) albopictus invasion of the Australasian Region. PLoS Negl Trop Dis. 2013;7(8):e2361.
Minard G, Tran FH, Dubost A, Tran-Van V, Mavingui P, Moro CV. Pyrosequencing 16S rRNA genes of bacteria associated with wild tiger mosquito Aedes albopictus: a pilot study. Front Cell Infect Microbiol. 2014;4:59.
Graham RI, Wilson K. Male-killing Wolbachia and mitochondrial selective sweep in a migratory African insect. BMC Evol Biol. 2012;12:204.
Mattingly PF. The sub-genus Stegomyia (Diptera: Culicidae) in the Ethiopian Region II. Distribution of species confined to the East and South African sub-region. Bull Br Mus (Natural History) Entomol. 1953;3(1):1–65.
Sabatini A, Raineri V, Trovato G, Coluzzi M. Aedes albopictus in Italy and possible diffusion of the species into the Mediterranean area. Parassitologia. 1990;32(3):301–4.
Rueda LM. Pictorial keys for the identification of mosquitoes (Diptera: Culicidae) associated with Dengue virus transmission. Zootaxa. 2004;589:1–60.
Baruffi L, Damiani G, Guglielmino CR, Bandi C, Malacrida AR, Gasperi G. Polymorphism within and between populations of Ceratitis capitata: comparison between RAPD and multilocus enzyme electrophoresis data. Heredity. 1995;74(4):425–37.
De Jong L, Moreau X, Dalia J, Coustau C, Thiéry A. Molecular characterization of the invasive Asian tiger mosquito, Aedes (Stegomyia) albopictus (Diptera: Culicidae) in Corsica. Acta Trop. 2009;112(3):266–9.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.
Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10(3):564–7.
Dupanloup I, Schneider S, Excoffier L. A simulated annealing approach to define the genetic structure of populations. Mol Ecol. 2002;11(12):2571–81.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.
Martins WS, Lucas DC, Neves KF, Bertioli DJ. WebSat - A web software for microsatellite marker development. Bioinformation. 2009;3(6):282–3.
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.
Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000;132:365–86.
Matschiner M, Salzburger W. TANDEM: integrating automated allele binning into genetics and genomics workflows. Bioinformatics. 2009;25(15):1982–3.
Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research - an update. Bioinformatics. 2012;28(19):2537–9.
Dieringer D, Schlötterer C. Microsatellite analyser (MSA): a platform independent analysis tool for large microsatellite data sets. Mol Ecol Notes. 2003;3(1):167–9.
Rousset F. genepop’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour. 2008;8(1):103–6.
Raymond M, Rousset F. GENEPOP (Version 1.2): Population genetics software for exact tests and ecumenicism. J Hered. 1995;86(3):248–9.
Guo SW, Thompson EA. Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics. 1992;48(2):361–72.
Kalinowski ST, Taper ML, Marshall TC. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007;16(5):1099–106.
Rice WR. Analyzing Tables of Statistical Tests. Evolution. 1989;43(1):223–5.
Delatte H, Toty C, Boyer S, Bouetard A, Bastien F, Fontenille D. Evidence of habitat structuring Aedes albopictus populations in Réunion Island. PLoS Negl Trop Dis. 2013;7(3):e2111.
Porretta D, Gargani M, Bellini R, Calvitti M, Urbanelli S. Isolation of microsatellite markers in the tiger mosquito Aedes albopictus (Skuse). Mol Ecol Notes. 2006;6(3):880–1.
This research is part of the EU-FP7 Research Infrastructures project (INFRAVEC) Grant N° 228421 and was partially supported by the Italian Ministry of Health (RF-2010-2318965) and the Fondazione Banca del Monte di Lombardia (Pavia, Italy). The study was performed within the FAO/IAEA research CRP programme “Exploring Genetic, Molecular, Mechanical and Behavioural Methods of Sex Separation in Mosquitoes (FS)”. Sincere thanks go to Dr. Romeo Bellini (Centro Agricoltura Ambiente “G. Nicoli”, Crevalcore, Italy), for providing Aedes albopictus specimens.
The authors declare that they have no competing interests.
ARM, GG, CRG, MM and LMG contributed to the study design. PS, GT, LMG and FS contributed to field mosquito egg collection. MM, LMG, NA, FS and GT performed the experimental work. MM, LMG, CRG and ARM analysed the data. ARM, CRG, MM, LMG and GG wrote the manuscript. All authors read and approved the final manuscript.
About this article
- Aedes albopictus
- Simple sequence repeats (SSRs)
- Internal transcribed spacer 2 (ITS2)
- Virus vectors
- Population genetics