Skip to main content

High genetic diversity but no geographical structure of Aedes albopictus populations in Réunion Island

Abstract

Background

In recent years, the Asian tiger mosquito Aedes albopictus has emerged as a species of major medical concern following its global expansion and involvement in many arbovirus outbreaks. On Réunion Island, Ae. albopictus was responsible for a large chikungunya outbreak in 2005–2006 and more recently an epidemic of dengue which began at the end of 2017 and is still ongoing at the time of writing. This dengue epidemic has seen a high number of human cases in south and west coastal regions, while few cases have been reported in the north and east of the island. To better understand the role of mosquito populations in such spatial patterns of dengue virus transmission in Réunion Island, we examined the genetic diversity and population structure of Ae. albopictus sampled across the island.

Results

Between November 2016 and March 2017, a total of 564 mosquitoes were collected from 19 locations in three main climatic regions (West, East and Center) of Réunion Island and were genotyped using 16 microsatellite loci. A high genetic diversity was observed with 2–15 alleles per locus and the average number of alleles per population varying between 4.70–5.90. Almost all FIS values were significantly positive and correlated to individual relatedness within populations using a hierarchical clustering approach based on principal components analyses (HCPC). However, the largest part of genetic variance was among individuals within populations (97%) while only 3% of genetic variance was observed among populations within regions. Therefore, no distinguishable population structure or isolation by distance was evidenced, suggesting high rates of gene flow at the island scale.

Conclusions

Our results show high genetic diversity but no genetic structure of Ae. albopictus populations in Réunion Island thus reflecting frequent movements of mosquitoes between populations probably due to human activity. These data should help in the understanding of Ae. albopictus vector capacity and the design of effective mosquito control strategies.

Background

The Asian tiger mosquito, Aedes (Stegomyia) albopictus [1], originated from Southeast Asia [2] is currently among the most invasive mosquito species in the world [3, 4]. Since the second half of the 20th century, this mosquito species has spread worldwide except for Antarctica [5], probably mainly disseminating through the transportation of eggs via the international trade of used tires [3, 6] and exotic plants [7]. Its ecological plasticity in different traits including egg diapause, the ability to use natural or urban larval breeding sites [8, 9] and opportunistic feeding behavior [10, 11] may have aided dispersal and adaptation to newly colonized environments across many contrasting climatic conditions; ranging from tropical to temperate [4]. In addition to its aggressive daytime human-biting behavior, Ae. albopictus has been shown in laboratory tests to be a proven competent vector for over 20 arboviruses including yellow fever virus, West Nile virus [8] and Zika virus [12]. The continued activity and dispersal of this mosquito has significant importance to public health, for example, due to its involvement in recent arbovirus outbreaks such as chikungunya outbreaks in Réunion Island in 2005–2006 [13] and in Italy in 2007 [14] and 2017 [15].

Réunion Island, an overseas French territory, is a volcanic island of 2512 km2 located about 700 km east of Madagascar. It hosts about 850,000 inhabitants, mainly living in coastal regions. Population density and climate vary greatly across the island: the west coast is dry and densely populated, the east coast is wet with a lower human population density and the center of the island is a humid mountainous region. Twelve mosquito species (from four genera: Aedes, Anopheles, Culex and Orthopodomyia) including Ae. albopictus and Ae. aegypti are currently known to inhabit the island [16, 17]. However, Ae. albopictus is one of the most abundant mosquito species in Réunion Island and is common all over the island in urban, periurban and rural areas sometimes reaching up 1200 m [17, 18]. This mosquito species was first recorded in Réunion Island in 1913 [19] while the presence of Ae. aegypti has been reported since 1902 [20]. Aedes albopictus gradually replaced Ae. aegypti whose distributional range was considerably reduced [17, 21, 22] probably due to dichlorodiphenyltrichloroethane (DDT) treatments for mosquito-control campaigns against malaria vectors in the 1950s [16].

The Asian tiger mosquito has become a major concern for public health in Réunion Island because it has been recognized as the main vector in several arbovirus outbreaks. In 2005–2006, this island has experienced a chikungunya virus outbreak that was responsible of 265,000 human cases (34% of the population) with 237 associated deaths [23]. More recently, an epidemic caused by a dengue virus type-2 (DENV-2) started at the end of 2017 and is still ongoing on the island [24, 25]. The French public health agency has reported more than 49,000 cases of dengue since the beginning of the epidemic, with most of the cases reported in the south and west coasts of the island [24]. As vector capacity for mosquito-borne viruses depends largely on vector and virus genetics [26,27,28,29], we asked whether genetic structure of Ae. albopictus populations may explain the observed geographical distribution of human cases of dengue in Réunion Island. Previous studies of Ae. albopictus population structure from Réunion Island showed contradictory results. Paupy et al. [30] and Delatte et al. [9], using allozymes and microsatellite loci, respectively, found that populations from the west and east coasts were genetically differentiated while, more recently, Sherpa et al. [31] did not find evidence for any genetic structure using 1561 single nucleotide polymorphisms (SNPs).

The present study aims at a better understanding the population structure of Ae. albopictus in Réunion Island. For this purpose, 564 mosquitoes from 19 populations sampled in different locations across the island (West, East and Center), were genotyped using 16 microsatellite loci that have previously shown their effectiveness in describing population structure of Ae. albopictus [32,33,34,35]. Genetic diversity as well as gene flow among populations were then assessed. This study could help understanding patterns of dengue transmission in Réunion Island in order to implement effective strategies to protect human populations.

Methods

Mosquito samples

Aedes albopictus eggs were sampled from 19 sites across Réunion Island (Fig. 1) in November and December 2016 and in February and March 2017. Most of the samples were collected in periurban areas with most buildings being family homes with gardens or allotment gardens. Six ovitraps were randomly placed at each site and spaced 5–10 m apart. Eggs from the six ovitraps (200–1500 eggs) were pooled together according to the site and transported to an insect laboratory where they were reared to adulthood. After morphological identification of adults using Theobald’s criteria [36], individuals were stored at − 20 °C until molecular analyses. Although egg sampling can show some limitation for population genetics due to sampling of siblings, this method was best suited for the collection of the large number of Ae. albopictus specimens required and at different locations across the island. Each of the 19 sites was characterized according to the following environmental variables: altitude, temperature, rainfall, vegetation type and land cover (Additional file 1: Table S1). We used mean temperatures normalized over 30 years (1981–2010) from the French meteorological weather service (Météo-France, https://donneespubliques.meteofrance.fr/). Rainfall data were provided by climatic stations installed by the French Agricultural Research Centre for International Development (CIRAD). The major type of vegetation around the sampled sites was determined from a mapping established by the Conservatoire Botanique National de Mascarin. Finally, the land-use was assessed using a mapping developed by the CIRAD from satellite data obtained in 2016–2017. Multiple Factor Analysis (MFA) was conducted using the FactoMineR package v.1.41 [37] in R software [38] to gather the 19 sites into regions with similar environmental characteristics.

Fig. 1
figure 1

Map of sampling sites where Aedes albopictus mosquitoes were collected in Réunion Island. For each site, GPS coordinates and environmental characteristics are indicated in Additional file 1: Table S1. Solid lines: national roads; dotted lines: secondary roads (not all secondary roads are shown); thick grey lines: outline of cirques. The sampling sites are colored according to the climatic regions: blue, West; orange, East; green, Center (see Fig. 2 for more details on climatic regions). Population codes are given in parentheses

Microsatellite processing and genotyping

DNA was extracted from individual mosquitoes using a cetyltrimethylammonium bromide (CTAB) protocol [39]. A total of 564 individuals were genotyped (24–30 individuals per site). In each site, half of the analyzed individuals were males and the other half females. A set of 16 microsatellite loci previously developed for Ae. albopictus was used [32, 34] (Additional file 2: Table S2). Microsatellite locus primers were pooled in three PCR-mixes based on non-overlapping size ranges and amplified by multiplex PCRs using fluorescent primers (Additional file 2: Table S2). All PCR reactions (10 µl) contained 0.5 U of Taq DNA polymerase, 6 µM of dNTPs, 37.5 µM of MgCl2, 1.2–7 µM of each primer (Additional file 2: Table S2) and 10 ng of mosquito DNA. The following program was used to amplify all 16 loci: 10 min at 95 °C followed by 35 cycles at 95 °C for 30 s, 57 °C for 30 s and 72 °C for 1 min and a final elongation step of 72 °C for 10 min. PCR products were run on a 3730XL DNA Genetic Analyzer (Applied Biosystems, California, USA) with a GeneScan 500LIZ size standard. Microsatellite alleles were read using the software GeneMapper v. 4.0 (Applied Biosystems).

Statistical analyses

Genetic diversity

The presence of null alleles was checked for each of the 16 microsatellite loci at the population level using Micro-Checker v.2.2.3 [40]. Multi-locus genotypes (MLGs) were assessed using the allelematch R package v.2.5 [38, 41]. The inbreeding coefficient (FIS) was assessed using Arlequin v.3.5 [42]. The average and effective number of alleles (Na and Ne respectively), the number of private alleles (NP) and the observed and expected heterozygosities (HO and HE respectively) were estimated using GenAlEx v.6.5 [43]. Tests of Hardy-Weinberg equilibrium (HWE) and linkage disequilibrium (LD) were performed using GENEPOP v.4.2 on the web [44, 45] and test_LD functions in the genepop R package v.1.05 [44], respectively. Sequential Bonferroni corrections were applied for all multiple comparisons [46].

A large departure from panmixia and strong linkage disequilibrium were detected for all microsatellite loci that might reflect the existence of substructure in the populations. Such a cryptic structure could be due to the sampling method that can favor the sampling of close relatives (i.e. half-sibs or full-sibs). To detect possible close relatives in our samples, individuals of each population were assigned to different genetic clusters, potentially corresponding to families of related individuals, using a hierarchical clustering approach based on principal components analyses (HCPC function, the Facto MineR R package v.1.41 [37]). To limit the potential bias induced by the sampling of related individuals, 100 subsampled datasets were generated by randomly sampling a maximum of two individuals per genetic cluster in each population. Hardy-Weinberg equilibrium (HWE), linkage disequilibrium (LD), the effective number of alleles (Ne) and the expected heterozygosity (HE) were estimated for the 100 subsampled datasets, as described above. Although the subsampling procedure does not solve the problem of sampling relatives, it allows assessment of the potential effect of relatedness in the dataset in population genetic analyses.

Population structure

Unsupervised clustering

Aedes albopictus population structure was investigated using principle component analysis and the Bayesian clustering algorithm implemented in STRUCTURE v.2.3.4 [47], with the admixture model considering correlated allele frequencies. The number of groups (K) was varied from 1–19. A ‛burn-in’ period of 104 iterations was followed by 10 runs of 105 iterations leading to the estimation of the membership coefficients. Runs were pooled using the CLUMPAK server [48]. The choice of the best number of clusters (K) was based on the Evanno’s method [49].

Population differentiation

For each pair of populations, the fixation index (Fst) was computed following Weir & Cokerham method [50]. For each microsatellite locus, the allelic frequencies were computed either (i) without considering individual relatedness using GenAlEx v.6.5 [43] or (ii) by taking the relationships among individuals into account using Colony v.2.0.6.4, likelihood method [51]. FST values were computed on the basis of the allelic frequencies estimated with both methods to detect the potential effect of relatedness on the estimation of genetic differentiation.

The distribution of genetic variation among individuals, populations and geographical regions was assessed in hierarchical analyses of molecular variance (AMOVA; [52, 53]) using GenAlEx v.6.5 [43]. Three geographical regions (West, East and Center) were considered according to environmental characterization of the sampling sites.

Isolation by distance was analyzed using R software [38], computing the linear relationship between pairwise estimates of FST and the geographical distance by road. Sequential Bonferroni corrections were applied for all multiple comparisons [46].

Results

To characterize the genetic differentiation of Ae. albopictus populations, mosquitoes from 19 sites across Réunion Island were genotyped (Fig. 1). To cluster the sites according to their environmental characteristics, multiple factor analysis (MFA) was performed based on five factors: two climatic variables (temperature and rainfall), altitude, land use and the type of vegetation (Fig. 2). The five factors contributed equivalently to the first axis of the MFA. The land-use and the vegetation type contributed more to the second axis (37% for each of the two factors) than the climatic variables (21% for both temperature and rainfall) and the altitude (5%). The first axis of the MFA divided the sites according to the temperature and the altitude, while the second axis divided the sites according to the rainfall. The 19 sites were thus clustered into three regions (Fig. 2): the West coast region with nine sites (W-PRO, W-LPO, W-ERM, W-P3B, W-ESL, W-PLA, W-LDP, W-PGB and W-SJO); the East coast region with six sites (E-PDA, E-PBS, E-PBSB, E-PNDL, E-PCD and E-PCP) and the Center region with four sites (C-PTC, C-PHY, C-PDP and C-PSA). The west coast was characterized by altitude ranging from 9–437 m a.s.l., temperatures from 19–25 °C and rainfall from 490–1540 mm (Additional file 1: Table S1). The east coast was wetter than the west coast and the center, with altitude ranging from 20–314 m a.s.l., temperatures from 20–24 °C and rainfall from 2130–3700 mm. The center of the island was mainly characterized by sites with high altitude, ranging from 449–1161 m a.s.l., temperatures from 14–19 °C and rainfall from 1150–3830 mm.

Fig. 2
figure 2

A two-dimensional plot from a multiple factor analysis (MFA) performed on the 19 Aedes albopictus sampling sites. The analysis was based on five environmental variables: two climatic variables (temperature and rainfall), altitude, land use and type of vegetation

High genetic diversity of Ae. albopictus in Réunion Island

Overall 564 mosquitoes from the 19 sites (24–30 mosquitoes per site) were genotyped using 16 microsatellite loci [32, 34]. For six microsatellite loci (Aealbmic4, Aealbmic5, Aealbmic12, Aealbmic13, Albdi6, Albtri3), null alleles were observed in more than 70% of the populations (Additional file 3: Table S3). These six microsatellite loci were therefore removed from downstream population genetics analyses. We also removed 32 individuals with two or more missing alleles in their multilocus genotypes (MLGs). Among the remaining 532 individuals, 530 different MLGs (i.e. differentiated by at least one allele) were identified; 528 were unique (i.e. identified in only one individual), while two MLGs were shared, each by two mosquitoes from the same population. Only unique MLGs were kept for population genetics analyses. Thus, our final dataset consisted of 530 individuals genotyped with 10 microsatellite loci (Aealbmic2, Aealbmic3, Aealbmic6, Aealbmic7, Aealbmic8, Aealbmic9, Aealbmic10, Aealbmic11, Aealbmic16 and Albtri45).

Allelic richness of the 10 microsatellite loci was high with the number of alleles per locus varying from 2 (Aealbmic2) to 15 (Aealbmic3) and the average number of alleles per population (Na) varying from 4.70 (in W-P3B) to 5.90 (in W-LDP and W-PRO; Additional file 4: Table S4).

Cryptic intrapopulation structure due to individual relatedness

In each population, the average number of effective alleles (Ne) was lower than the observed average number of alleles (Na) thus reflecting imbalanced allelic frequencies in populations. Indeed, for most of the loci, one to three alleles were overrepresented within populations (Fig. 3). Besides, almost all alleles were shared by mosquitoes from different populations, very few private alleles were observed: two in W-ERM (Aealbmic11 and Aealbmic3) and in W-PRO (Aealbmic16 and Aealbmic9) and one in W-LPO, C-PHY (Aealbmic8) and in C-PSA (Aealbmic11; Additional file 4: Table S4). The observed heterozygosity (HO) varied from 0.47–0.62 while the expected heterozygosity (HE) varied from 0.55–0.66. In all populations, HE was higher than HO and the inbreeding coefficients (FIS) were significantly positive, indicating a deficit in heterozygotes except for E-PNDL and C-PHY for which it was not significantly different from zero (Additional file 4: Table S4). Moreover, significant linkage disequilibrium was observed between all pairs of loci (Chi-square test; 46 < χ2 < 127; df = 38; P < 0.013). However, no deviation from HWE and no evidence of linkage disequilibrium were observed after sequential Bonferroni correction for multiple testing.

Fig. 3
figure 3

Allelic frequencies for 10 microsatellite loci in the 19 mosquito populations. Only alleles with a frequency of at least 25% in one or more populations are colored. Numbers correspond to the size of alleles for each microsatellite locus

We therefore examined whether there were clusters of individuals within each population. For this purpose, we used a hierarchical clustering approach based on principal components analyses (HCPC) that grouped individuals according to their genetic proximity (Additional file 5: Figure S1). Two (E-PBS, E-PCP and C-PDP) to six (E-PBSB) clusters (with 1–24 individuals per cluster) were found per population (Table 1). The presence of clusters of individuals within each population suggests the presence of relatives. A high number of relatives within each genetic cluster may explain the observed deviations from HWE and the significant linkage disequilibrium between microsatellite markers. To break down this clustering effect due to individualsʼ relatedness (e.g. half-sib structure), we generated subsampled datasets by randomly sampling only two individuals per genetic cluster within each population (except for clusters with only one individual). Overall 100 subsampled datasets (with 4–10 individuals per population and 118 individuals per dataset) were generated. Each subsampled dataset was further used to perform population genetics analyses. Low sample size per population was used (4–10 individuals, i.e. no more than 2 individuals per cluster) because our preliminary analyses showed that including a higher number of individuals from the same cluster increased the probability to observe a departure from panmixia. Most of the populations appeared to be at HWE (Table 1). Indeed, all populations were at HWE in more than 60 subsampled datasets over the 100 subsampled datasets tested for HWE tests, except W-LDP. In addition, no linkage disequilibrium was observed between all pairs of loci (Additional file 6: Figure S2).

Table 1 Summary statistics in each Aedes albopictus population estimated using the 100 subsampled datasets

No genetic structure of Ae. albopictus populations

The fixation index, FST, was then used to investigate population differentiation. For each pair of populations, FST was computed by using the full dataset directly or the likelihood method (see “Methods” section) that allows relatedness among individuals to be taken into account when estimating the allelic frequencies. Both methods showed similar results (Pearson’s product moment correlation, R2 = 0.79; P < 0.001; n = 171; Additional file 7: Figure S3). Therefore, only the results obtained with the full dataset (i.e. without taking relatedness among individuals into account) were considered. The FST values ranged from 0.011–0.057 and were significant in most pairwise comparisons (Table 2). However, no pairwise comparison remained significantly different from 0 after the sequential Bonferroni correction. The average FST among populations of the West coast region (FST = 0.02 ± 0.01) was lower than that measured among populations of the East coast or Center (FST = 0.03 ± 0.004 and 0.04 ± 0.01, respectively). A linear regression model between matrices of FST and geographical distances was not significant (R2 = 1.36 × 10−4, P = 0.88; n = 171; Additional file 8: Figure S4).

Table 2 Pairwise FST between the 19 Aedes albopictus populations calculated using the Weir and Cockerham’s method

Analysis of molecular variance (AMOVA; Table 3) across the 19 Ae. albopictus populations considering the three regions (West, East and Center) revealed that 97% of the variance was found among individuals within populations (P = 0.001), while only 3% was found between populations (P = 0.001) and 0% between regions (P = 0.80; Table 3). Finally, neither principal components analysis (Fig. 4a) nor Bayesian clustering approach (STRUCTURE v.2.3.4; Fig. 4b) revealed any population structure.

Table 3 Analyses of molecular variance (AMOVA) for the 19 Aedes albopictus populations
Fig. 4
figure 4

Genetic structure of Aedes albopictus in Réunion Island. a A two-dimensional plot from a principal component analysis (PCA) based on individual genotypes for 10 microsatellite loci. b Bayesian analysis of population structure for the 19 mosquito populations produced using STRUCTURE. Each vertical bar represents an individual, each color represents a cluster (K) and the color of the bar indicates the probability of assignment to each cluster. Only the results for K = 2 to K = 4 are presented

Discussion

We assessed the genetic diversity and the genetic structure of 19 Ae. albopictus populations from Réunion Island using 10 microsatellite loci. We highlighted a high genetic diversity but found no evidence of genetic differentiation between populations despite the contrasted environments between the three regions (West, East and Center) where mosquitoes were collected.

A high genetic diversity reflecting long-term established mosquito populations

The genetic diversity observed in the 19 Ae. albopictus populations was high: on average more than five alleles were observed per population and 530 MLGs were identified among the 532 mosquitoes examined. Using microsatellite loci, Manni et al. [35] observed a higher genetic diversity in Ae. albopictus populations from Réunion Island compared to populations from Asia (e.g. from Japan and China), Europe (from Greece, Albania, Italy) and America (from Virginia, USA). A recent study based on 1561 genome-wide SNPs also described a higher genetic diversity in Ae. albopictus populations from Réunion Island compared to European populations [31]. Two main hypotheses may explain the highest genetic diversity of Ae. albopictus populations from Réunion Island. The first hypothesis is that Ae. albopictus populations in Réunion Island result from an admixture of several populations from different geographical origins. According to Manni et al. [35] investigations of population genetics and Approximate Bayesian Computation (ABC) analyses, Ae. albopictus populations from Réunion Island might result from admixture between populations from Thailand and Japan. By comparing mtDNA polymorphism between ancient (collected in 1956) and more recent samples (collected in 2007), Delatte et al. [54] suggested that two waves of Ae. albopictus invasion may have occurred in Réunion Island, the first one several centuries ago and the second one in the 1990s, leading to a replacement of populations. The assumption of population replacement from 1990s seems an unlikely explanation for our results because the genetic diversity in Réunion Island is higher than observed in Albania [31, 35] where Ae. albopictus has been recorded since 1979 [55] or in Italy where this mosquito species has been established since 1990 [56]. However, it is likely that several introduction events of Ae. albopictus individuals from different geographical origins might have occurred in Réunion Island. The second hypothesis is that long-term established mosquito populations would have favored the accumulation of genetic variability within populations, countering the effects of any initial demographic bottleneck. This would be supported by the fact that Ae. albopictus was first recorded in Réunion Island in 1913 [19] and it became one of the most abundant species in the island.

No genetic structure revealing the impact of humans in the distribution of Ae. albopictus

The analysis of the 19 mosquito populations with the 10 genotyped microsatellite loci revealed significant linkage disequilibrium between all pairs of loci. This could be explained by the imbalanced allelic frequencies observed in all populations and for all loci: the most frequent alleles are more often associated, regardless of their position on chromosomes. The over-representation of some alleles may be due to a founder effect (only few individuals have contributed to the establishment of the populations) but may also reflect the fact that each of the sampled populations was probably represented by the progenies from few females. Indeed, we measured significant positive FIS in almost all populations reflecting a deficit in heterozygotes. Ignoring closely related individuals can lead to bias in the estimation of allelic frequencies and thus in the resulting population genetics analyses such as tests of HWE and linkage disequilibrium [57]. In this study, the use of the full dataset for population genetics analyses led to significant deviations from HWE and linkage disequilibrium. However, when genetically close related individuals were removed from the dataset, almost all populations were in panmixia and no linkage disequilibrium was observed between pairs of microsatellite loci. Previous studies based on allozymes also reported heterozygote deficits leading to significant deviation from HWE [30, 58,59,60,61,62]. Significant positive FIS values were also observed in several studies based on microsatellites (FIS ranging from 0.10–0.20; [9, 33, 34]) and SNPs (mean FIS of 0.19 for populations from Réunion Island; [31]) in populations from both native and non-native areas of the Asian tiger mosquito. These positive FIS values could be explained by inbreeding resulting from preferential mating among relatives or by sampling methods [63]. Indeed, we used samples collected as eggs in ovitraps (six ovitraps for each site) and eggs collected in the same ovitrap could belong to the same progeny. Moreover, Ae. albopictus has low dispersal capacities (lower than 500 m) [64, 65] and females can lay eggs in several containers [66]. Hence, it is likely that very few females laid eggs in the six ovitraps used in each sampled site. Therefore, although we genotyped a high number of mosquitoes in each population (24–30 mosquitoes per population), it is likely that this does not represent the sampling site diversity. In future investigations, it would be interesting to use adults (which is not straightforward with Ae. albopictus) or eggs sampled at different times to increase the number of sampled progenies. Anyway, we recommend testing for relatedness between individuals when significantly positive FIS are detected in populations. This is all the more relevant given the rapid spread of Ae. albopictus and its involvement in arbovirus outbreaks that leads to an increasing number of studies to understand its genetics and adaptation.

Both the Bayesian clustering algorithm implemented in STRUCTURE and the analysis of molecular variance (AMOVA) revealed an absence of population structure despite a high genetic diversity and a great variability of bioclimatic conditions (climate, vegetation type or anthropogenic activities). This suggests frequent exchange between individuals from different locations across the island. The absence of isolation by distance strongly suggested that human-assisted long-distance gene flows are frequent and play a key role in the dispersal of the Asian tiger mosquito in Réunion Island as has already been shown elsewhere at both large [67, 68] and local scales [69]. Our results are consistent with those of Sherpa et al. [31] who did not detect differentiation among populations from Réunion Island using thousands of genome-wide SNPs. However, previous investigations in Réunion Island using allozymes [30] or microsatellite loci (designed for Ae. aegypti [9]) highlighted genetic differentiation between Ae. albopictus populations from the east and west coasts related to ecological factors [9, 30]. Apart from molecular markers, the difference between these previous studies and our results could be explained by sampling methods. Indeed, we used only eggs from ovitraps (which are artificial sites) that were laid in periurban areas whereas different stages (larvae, pupae and adults) collected in both natural (gully, rock or tree holes) and artificial (tires or flower pots) sites were used in these previous investigations [9, 30].

Conclusions

In this study, we observed a high genetic diversity in Ae. albopictus populations from Réunion Island, probably due to their long-term establishment. We also observed a lack of genetic structure reflecting a strong connectivity among populations. The absence of population structure and isolation by distance represent key information for the management of vector control strategies against Ae. albopictus in Réunion Island. Vector control strategies could be applied at the island scale by taking into account environmental characteristics and mosquito population densities through time because gene flow between mosquito populations may have a limited impact in their effectiveness. Additionally, these results increase our understanding of the evolutionary history of the Asian tiger mosquito in Réunion Island and will allow deciphering its vector capacity. Further studies on vector competence using the epidemic DENV-2 strain and Ae. albopictus populations sampled across Réunion Island should help understanding the transmission patterns of dengue virus currently causing an outbreak in the island. As vector capacity can also be affected by extrinsic or abiotic factors such as temperature, rainfall, vector density, vector survival and the probability a vector feeds on a host [29], these factors should also be taken into account to better understand the transmission patterns of dengue virus in Réunion Island.

Availability of data and materials

Data supporting the conclusions of this article are included within the article and its additional files. Raw data are available on Zenodo (https://doi.org/10.5281/zenodo.3475984).

Abbreviations

HCPC:

hierarchical clustering approach based on principal components analyses

DENV–2:

dengue virus type-2

SNPs:

single nucleotide polymorphisms

MFA:

multiple factor analysis

MLGs:

multilocus genotypes

HWE:

Hardy-Weinberg equilibrium

LD:

linkage disequilibrium

References

  1. Skuse FA. The banded mosquito of Bengal. Indian Museum Notes. 1894;3:1–20.

    Google Scholar 

  2. Smith CE. The history of dengue in tropical Asia and its probable relationship to the mosquito Aedes aegypti. J Trop Med Hyg. 1956;59:243–51.

    CAS  PubMed  Google Scholar 

  3. 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:76–85.

    Article  PubMed  Google Scholar 

  4. Bonizzoni M, Gasperi G, Chen X, James AA. The invasive mosquito species Aedes albopictus: current knowledge and future perspectives. Trends Parasitol. 2013;29:460–8.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Kraemer MU, Sinka ME, Duda KA, Mylne A, Shearer FM, Barker CM, et al. The global distribution of the arbovirus vectors Aedes aegypti and Ae. albopictus. Elife. 2015;4:e08347.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Reiter P, Sprenger D. The used tire trade: a mechanism for the worldwide dispersal of container breeding mosquitoes. J Am Mosq Control Assoc. 1987;3:494–501.

    CAS  PubMed  Google Scholar 

  7. Enserink M. A mosquito goes global. Science. 2008;320:864–6.

    Article  CAS  PubMed  Google Scholar 

  8. Paupy C, Delatte H, Bagny L, Corbel V, Fontenille D. Aedes albopictus, an arbovirus vector: from the darkness to the light. Microbes Infect. 2009;11:1177–85.

    Article  CAS  PubMed  Google Scholar 

  9. 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:e2111.

    Article  PubMed  PubMed Central  Google Scholar 

  10. 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:249–58.

    Article  PubMed  Google Scholar 

  11. Kamgang B, Nchoutpouen E, Simard F, Paupy C. Notes on the blood-feeding behavior of Aedes albopictus (Diptera: Culicidae) in Cameroon. Parasit Vectors. 2012;5:57.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Wong PS, Li MZ, Chong CS, Ng LC, Tan CH. Aedes (Stegomyia) albopictus (Skuse): a potential vector of Zika virus in Singapore. PLoS Negl Trop Dis. 2013;7:e2348.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Josseran L, Paquet C, Zehgnoun A, Caillere N, Le Tertre A, Solet JL, et al. Chikungunya disease outbreak, Réunion Island. Emerg Infect Dis. 2006;12:1994–5.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Rezza G, Nicoletti L, Angelini R, Romi R, Finarelli A, Panning M, et al. Infection with chikungunya virus in Italy: an outbreak in a temperate region. Lancet. 2007;370:1840–6.

    Article  CAS  PubMed  Google Scholar 

  15. Venturi G, Di Luca M, Fortuna C, Remoli ME, Riccardo F, Severini F, et al. Detection of a chikungunya outbreak in central Italy, August to September 2017. EuroSurveillance. 2017;22:17–00646.

    Article  PubMed Central  Google Scholar 

  16. Hamon J. Etude biologique et systématique des Culicidae de l’ile de La Réunion. Mem Inst Sci Madagascar Série E. 1953;4:521–41.

    Google Scholar 

  17. Boussès P, Dehecq JS, Brengues C, Fontenille D. Inventaire actualisé des moustiques (Diptera: Culicidae) de l’île de La Réunion, océan Indien. Bull Soc Pathol Exot. 2013;106:113–25.

    Article  PubMed  Google Scholar 

  18. Delatte H, Paupy C, Dehecq JS, Thiria J, Failloux AB, Fontenille D. Aedes albopictus, vector of chikungunya and dengue viruses in Réunion Island: biology and control. Parasite. 2008;15:3–13.

    Article  CAS  PubMed  Google Scholar 

  19. Edwards FW. Notes on the mosquitos of Madagascar, Mauritius and Réunion. Bull Entomol Res. 1920;11:133–8.

    Article  Google Scholar 

  20. Vassal JJ. Le paludisme à l’île de La Réunion. Per Gli Stud Della Maria. 1907;8:18–27.

    Google Scholar 

  21. Salvan M, Mouchet J. Aedes albopictus et Aedes aegypti à l’île de la Réunion. Ann Soc Belg Med Trop. 1994;74:323–6.

    CAS  PubMed  Google Scholar 

  22. Bagny L, Delatte H, Quilici S, Fontenille D. Progressive decrease in Aedes aegypti distribution in Réunion Island since the 1900s. J Med Entomol. 2009;46:1541–5.

    Article  PubMed  Google Scholar 

  23. Renault P, Solet JL, Sissoko D, Balleydier E, Larrieu S, Filleul L, et al. A major epidemic of chikungunya virus infection on Réunion Island, France, 2005–2006. Am J Trop Med Hyg. 2007;77:727–31.

    Article  PubMed  Google Scholar 

  24. Santé Publique France. Le point épidémio / Dengue à La Réunion / Point de situation no 19 -19/03/2019. https://www.santepubliquefrance.fr/regions/ocean-indien/documents/bulletin-regional/2019/surveillance-de-la-dengue-a-la-reunion.-point-epidemiologique-au-19-mars-2019.

  25. Pascalis H, Turpin J, Roche M, Krejbich P, Gadea G, Atyame CM, et al. The epidemic of dengue virus type-2 cosmopolitan genotype on Réunion Island relates to its active circulation in the southwestern Indian Ocean neighboring islands. Heliyon. 2019;5:e01455.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Lambrechts L, Halbert J, Durand P, Gouagna LC, Koella JC. Host genotype by parasite genotype interactions underlying the resistance of anopheline mosquitoes to Plasmodium falciparum. Malar J. 2005;4:3.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Lambrechts L, Chevillon C, Albright RG, Thaisomboonsuk B, Richardson JH, Jarman RG, et al. Genetic specificity and potential for local adaptation between dengue viruses and mosquito vectors. BMC Evol Biol. 2009;9:160.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Fansiri T, Fontaine A, Diancourt L, Caro V, Thaisomboonsuk B, Richardson JH, et al. Genetic mapping of specific interactions between Aedes aegypti mosquitoes and dengue viruses. PLoS Genet. 2013;9:e1003621.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Kramer LD, Ciota AT. Dissecting vectorial capacity for mosquito-borne viruses. Curr Opin Virol. 2015;15:112–8.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Paupy C, Girod R, Salvan M, Rodhain F, Failloux AB. Population structure of Aedes albopictus from La Réunion Island (Indian Ocean) with respect to susceptibility to a dengue virus. Heredity. 2001;87:273–83.

    Article  CAS  PubMed  Google Scholar 

  31. Sherpa S, Rioux D, Pougnet-Lagarde C, Després L. Genetic diversity and distribution differ between long-established and recently introduced populations in the invasive mosquito Aedes albopictus. Infect Genet Evol. 2018;58:145–56.

    Article  PubMed  Google Scholar 

  32. 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:e2361.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Minard G, Tran FH, Van VT, Goubert C, Bellet C, Lambert G, et al. French invasive Asian tiger mosquito populations harbor reduced bacterial microbiota and genetic diversity compared to Vietnamese autochthonous relatives. Front Microbiol. 2015;6:970.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Manni M, Gomulski LM, Aketarawong N, Tait G, Scolari F, Somboon P, et al. Molecular markers for analyses of intraspecific genetic diversity in the Asian tiger mosquito, Aedes albopictus. Parasit Vectors. 2015;8:188.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Manni M, Guglielmino CR, Scolari F, Vega-Rúa A, Failloux AB, Somboon P, et al. Genetic evidence for a worldwide chaotic dispersion pattern of the arbovirus vector, Aedes albopictus. PLoS Negl Trop Dis. 2017;11:e0005332.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Theobald FV. A monograph of the Culicidae or mosquitoes. London: The British Museum (Natural History); 1910.

    Google Scholar 

  37. Lê S, Josse J, Husson F. FactoMineR: an R package for multivariate analysis. J Stat Softw. 2008;25:1–18.

    Article  Google Scholar 

  38. R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018.

    Google Scholar 

  39. Rogers SO, Bendich AJ. Extraction of DNA from plant tissues. In: Gelvin S, Schilperoort RA, editors. Plant molecular biology manual. Boston: Kluwer Academic Publishers; 1988. p. 1–10.

    Google Scholar 

  40. Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535–8.

    Article  CAS  Google Scholar 

  41. Galpern P, Manseau M, Hettinga P, Smith K, Wilson P. Allelematch: an R package for identifying unique multilocus genotypes where genotyping error and missing data may be present. Mol Ecol Resour. 2012;12:771–8.

    Article  PubMed  Google Scholar 

  42. Excoffier L, Lischer HE. Arlequin ver. 3.5: a new series of programs to preform population genetics analysis under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.

    Article  PubMed  Google Scholar 

  43. Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics. 2012;28:2537–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Rousset F. Genepop’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour. 2008;8:103–6.

    Article  PubMed  Google Scholar 

  45. Raymond M, Rousset F. GENEPOP (Version 1.2): population genetics software for exact tests and ecumenicism. J Hered. 1995;86:248–9.

    Article  Google Scholar 

  46. Holm S. A simple sequentially rejective multiple test procedure. Scand J Stat. 1979;6:65–70.

    Google Scholar 

  47. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. Kopelman NM, Mayzel J, Jakobsson M, Rosenberg NA, Mayrose I. Clumpak: a program for identifying clustering modes and packaging population structure inferences across K. Mol Ecol Resour. 2015;15:1179–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14:2611–20.

    Article  CAS  PubMed  Google Scholar 

  50. Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38:1358–70.

    CAS  PubMed  Google Scholar 

  51. Jones OR, Wang J. COLONY: A program for parentage and sibship inference from multilocus genotype data. Mol Ecol Resour. 2010;10:551–5.

    Article  PubMed  Google Scholar 

  52. Excoffier L, Smouse PE, Quattro JM. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992;131:479–91.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. Michalakis Y, Excoffier L. A generic estimation of population subdivision using distances between alleles with special reference for microsatellite loci. Genetics. 1996;142:1061–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  54. Delatte H, Bagny L, Brengue C, Bouetard A, Paupy C, Fontenille D. The invaders: phylogeography of dengue and chikungunya viruses Aedes vectors, on the South West islands of the Indian Ocean. Infect Genet Evol. 2011;11:1769–81.

    Article  CAS  PubMed  Google Scholar 

  55. Adhami J, Murati N. The presence of the mosquito Aedes albopictus in Albania. Rev Mjekësore. 1987;1:13–6.

    Google Scholar 

  56. 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:301–4.

    CAS  PubMed  Google Scholar 

  57. Wang J. Effects of sampling close relatives on some elementary population genetics analyses. Mol Ecol Resour. 2018;18:41–54.

    Article  PubMed  Google Scholar 

  58. Black WC, Ferrari JA, Rai KS, Sprenger D. Breeding structure of a colonising species: Aedes albopictus (Skuse) in the United States. Heredity. 1988;60:173–81.

    Article  PubMed  Google Scholar 

  59. Kambhampati S, Black WC, Rai KS. Geographic origin of the US and Brazilian Aedes albopictus inferred from allozyme analysis. Heredity. 1991;67:85–93.

    Article  PubMed  Google Scholar 

  60. Urbanelli S, Bellini R, Carrieri M, Sallicandro P, Celli G. Population structure of Aedes albopictus (Skuse): the mosquito which is colonizing Mediterranean countries. Heredity. 2000;84:331–7.

    Article  CAS  PubMed  Google Scholar 

  61. Vazeille M, Mousson L, Rakatoarivony I, Villeret R, Rodhain F, Duchemin JB, et al. Population genetic structure and competence as a vector for dengue type 2 virus of Aedes aegypti and Aedes albopictus from Madagascar. Am J Trop Med Hyg. 2001;65:491–7.

    Article  CAS  PubMed  Google Scholar 

  62. de Lourenço 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:105–14.

    Article  Google Scholar 

  63. Goubert C, Minard G, Vieira C, Boulesteix M. Population genetics of the Asian tiger mosquito Aedes albopictus, an invasive vector of human diseases. Heredity. 2016;117:125–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Marini F, Caputo B, Pombi M, Tarsitani G, della Torre A. Study of Aedes albopictus dispersal in Rome, Italy, using sticky traps in mark-release-recapture experiments. Med Vet Entomol. 2010;24:361–8.

    Article  CAS  PubMed  Google Scholar 

  65. Le Goff G, Damiens D, Ruttee AH, Payet L, Lebon C, Dehecq JS, et al. Field evaluation of seasonal trends in relative population sizes and dispersal pattern of Aedes albopictus males in support of the design of a sterile male release strategy. Parasit Vectors. 2019;12:81.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Hawley WA. The biology of Aedes albopictus. J Am Mosq Control Assoc Suppl. 1988;1:1–39.

    CAS  PubMed  Google Scholar 

  67. 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:284–95.

    Article  PubMed  Google Scholar 

  68. Egizi A, Kiser J, Abadam C, Fonseca DM. The hitchhiker’s guide to becoming invasive: exotic mosquitoes spread across a US state by human transport not autonomous flight. Mol Ecol. 2016;25:3033–47.

    Article  PubMed  Google Scholar 

  69. Schmidt TL, Rašić G, Zhang D, Zheng X, Xi Z, Hoffmann AA. Genome-wide SNPs reveal the drivers of gene flow in an urban population of the Asian tiger mosquito, Aedes albopictus. PLoS Negl Trop Dis. 2017;11:e0006009.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are very grateful to Cyrille Lebon for his help during the sampling of mosquitoes in the field; Guillaume Minard for his advice in the choice of microsatellite loci; Annelise Tran for providing us databases on mean temperatures and rainfall; David Wilkinson for the improvement of English language in the manuscript.

Funding

This study was supported by the FEDER INTERREG V OCEAN-INDIEN 2014-2020 under convention no. GURDTI/2017-0788-6877 and the consortium ZIKAlliance project (European Union-Horizon 2020 program under grant agreement no. 735548). The work of ACL was supported by the University of Réunion Island as A.T.E.R.

Author information

Authors and Affiliations

Authors

Contributions

CMA and PMa conceived and designed the experiments. ACL and CMA performed the experiments. ACL and PMi carried out the statistical analyses. ACL, CMA, PMi, HM and PMa wrote the paper. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Patrick Mavingui or Célestine M. Atyame.

Ethics declarations

Ethics approval and consent to participate

Not applicable. None of the samples used in this study were collected in protected areas and the Ae. albopictus mosquito is not considered as an endangered or protected species. Therefore, no specific permission was required for mosquito sampling.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1

. Environmental characteristics of the 19 sampling sites of Aedes albopictus in Réunion Island.

Additional file 2: Table S2

. Information on the microsatellite loci used for the genotyping of Aedes albopictus populations. Abbreviations: F, forward; R, reverse.

Additional file 3: Table S3

. Presence/absence of null alleles for each microsatellite locus in the 19 mosquito populations.

Additional file 4: Table S4

. Microsatellite variation in the 19 Aedes albopictus populations. Abbreviations: n, number of individuals analyzed; Na, number of alleles; Ne, effective number of alleles; Np, number of private alleles; HO, observed heterozygosities; HE, expected heterozygosities; FIS, inbreeding coefficient (*P < 0.05, **P < 0.01, ***P < 0.001); SE, standard error.

Additional file 5: Figure S1

. Hierarchical clustering of Aedes albopictus individuals from the same population. Example showing the case of the population C-PSA (n = 25 individuals). a Dendrogram generated by the hierarchical clustering. b Principal component analysis. Abbreviation: Ind, individual.

Additional file 6: Figure S2

. Distribution of linkage disequilibrium Chi-square statistic estimated for each couple of loci. In grey: over the 100 subsampled datasets; in black: for the full dataset.

Additional file 7: Figure S3.

Relationship between FST values for Aedes albopictus pairs of populations. a Comparisons of FST estimated using R (x-axis) and GenAlEx v.6.5 [43] (y-axis) without taking into account relatives. b Comparisons of FST estimated using R without taking relatives into account (x-axis) and by taking relatives into account (y-axis). Full line: regression line; dotted line: y = x.

Additional file 8: Figure S4.

Relationship between the FST and the distances by road calculated between pairs of mosquito populations. y = − 3.58 × 10−6 x + 2.61 × 10−2 (Pearsonʼs test: n = 171, t = − 0.15, R2 = 1.36 × 10−4, P = 0.88).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Latreille, A.C., Milesi, P., Magalon, H. et al. High genetic diversity but no geographical structure of Aedes albopictus populations in Réunion Island. Parasites Vectors 12, 597 (2019). https://doi.org/10.1186/s13071-019-3840-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-019-3840-x

Keywords