- Open Access
Genetic structure and gene flow of the flea Xenopsylla cheopis in Madagascar and Mayotte
Parasites & Vectors volume 10, Article number: 347 (2017)
The flea Xenopsylla cheopis (Siphonaptera: Pulicidae) is a vector of plague. Despite this insect’s medical importance, especially in Madagascar where plague is endemic, little is known about the organization of its natural populations. We undertook population genetic analyses (i) to determine the spatial genetic structure of X. cheopis in Madagascar and (ii) to determine the potential risk of plague introduction in the neighboring island of Mayotte.
We genotyped 205 fleas from 12 sites using nine microsatellite markers. Madagascan populations of X. cheopis differed, with the mean number of alleles per locus per population ranging from 1.78 to 4.44 and with moderate to high levels of genetic differentiation between populations. Three distinct genetic clusters were identified, with different geographical distributions but with some apparent gene flow between both islands and within Malagasy regions. The approximate Bayesian computation (ABC) used to test the predominant direction of flea dispersal implied a recent population introduction from Mayotte to Madagascar, which was estimated to have occurred between 1993 and 2012. The impact of this flea introduction in terms of plague transmission in Madagascar is unclear, but the low level of flea exchange between the two islands seems to keep Mayotte free of plague for now.
This study highlights the occurrence of genetic structure among populations of the flea vector of plague, X. cheopis, in Madagascar and suggests that a flea population from Mayotte has been introduced to Madagascar recently. As plague has not been reported in Mayotte, this introduction is unlikely to present a major concern for plague transmission. Nonetheless, evidence of connectivity among flea populations in the two islands indicates a possibility for dispersal by fleas in the opposite direction and thus a risk of plague introduction to Mayotte.
The oriental rat flea, Xenopsylla cheopis (Siphonaptera: Pulicidae) is a holometabolous insect ectoparasite and was first described in Egypt , which is believed to represent its origin . This species of flea is now cosmopolitan because of widespread dispersal (principally on ships) by its preferred rodent host, the black rat Rattus rattus (Rodentia: Muridae) . Xenopsylla cheopis is also a frequent parasite on the brown rat R. norvegicus and can parasitize other small mammals [4, 5]. This flea draws particular attention because of its role as a vector of pathogens responsible for human diseases such as plague and murine typhus [6, 7]. Indeed, it is thought to be the most efficient vector of the plague bacterium, Yersinia pestis  and can transmit the plague both between rodent hosts and to humans.
Plague is a re-emerging disease occurring in many regions of the World  but, more than 90% of worldwide cases in 2014–2015 were reported in Africa ; Madagascar being the most affected country . Plague was introduced to the eastern coastal region of Madagascar in 1898 , apparently by ships from India  and subsequently, spread to other ports before reaching the Central Highlands in 1921, where it became endemic  inside the “plague focus” (altitude >800 m) . Nevertheless, some areas outside this focus have epidemics, such as the District of Ikongo (altitude ~750 m) in 1998 [15, 16] and the District of Ambilobe (altitude <500 m) in 2011 . Although apparently absent from coastal areas since the 1930s, plague re-emerged in Mahajanga, a port in the north-west in 1991 with annual outbreaks between 1995 and 1999 . Eighty-six per cent of suspected cases reported in Madagascar between 2007 and 2011 were classified as bubonic plague , reflecting the important role played by flea vectors in the transmission cycle.
In Madagascar, Rattus rattus, the principal plague reservoir is abundant and found in diverse habitats . Two flea species (Pulicidae) are reported to be the main vectors of Y. pestis: Xenopsylla cheopis , and an endemic species Synopsyllus fonquerniei . Xenopsylla cheopis can be found throughout most of Madagascar, independent of altitude , while S. fonquerniei is largely absent below 800 m . Unlike S. fonquerniei, X. cheopis is predominantly found on R. rattus captured inside houses [15, 18, 22] and is therefore thought to be an important vector for human cases. Population genetic and phylogeographic studies have provided important insights into the invasion history and population ecology of R. rattus and Y. pestis in Madagascar. A large nationwide study of R. rattus indicated the existence of only two clearly defined genetic groups corresponding to two separate introduction events, with one largely confined to the far north of the island and one reflecting a large spatial expansion of an introduction that occurred in the south . In contrast, Y. pestis exhibited significant geographical separation among 15 identified subclades, implying largely local epidemiological cycles with limited gene flow . However, the same study also found evidence of long distance transfers, probably human-mediated .
Despite the importance of X. cheopis for plague transmission in Madagascar, no studies of its population genetic structure have been conducted. Most studies of X. cheopis involve laboratory studies of the flea-bacteria relationship , the host-flea relationship  or insecticide treatment effects [26,27,28]. Field studies on population structure and dispersal in Madagascar and surrounding areas are required to improve our understanding of vector dynamics and the associated epidemiological risks.
For example, studies of natural populations of the plague flea vector Oropsylla hirsuta, which transmits plague bacterium to the black-tailed prairie dog, highlighted the lack of isolation by distance and spatial genetic structure of the flea, and demonstrated that re-colonization of fleas from plague-free to plague zones occurred and caused a flea population expansion after epizootics . The implication is that this flea has reasonable dispersal ability at this spatial scale, with the high estimated rate of gene flow exhibited by O. hirsuta potentially associated with dispersal by its hosts (prairie dogs and other mammals) .
Although the population genetic structure of parasites is often linked to the population genetic structure and dispersal of their hosts, the extent of this congruence will depend on the intimacy of the parasite-host association , and factors such as host specificity and time spent on host. In many cases, parasites are expected to show stronger structuring than their hosts due to limited dispersal abilities and lower effective population sizes . A comparative study of the fur flea Listropsylla agrippinae and the nest flea Chiastopsylla rossi revealed different phylogeographic patterns, with the nest flea showing higher genetic divergence between sampling localities, presumably due to more restricted dispersal as a consequence of less time spent on hosts . In other cases, a lack of concordance in the phylogeographic and population genetic structure between some flea species and their hosts suggests the potentially important role of dispersal by other sympatric hosts [30, 34]. As a fur flea, adult X. cheopis are predicted to have relatively frequent opportunities to be dispersed by their hosts. However, in Madagascar, unlike its principal host, X. cheopis is strongly geographically restricted to houses  and is therefore spatially restricted to a small subset of R. rattus populations. Other mammals that are less abundant than Rattus rattus, such as Mus musculus and Suncus murinus, can be infested with X. cheopis . However, unlike R. rattus, these are peridomestic and largely restricted to around houses. Thus, we predict that X. cheopis populations will show much greater genetic structure than R. rattus populations due to a combination of limited dispersal and lower effective population sizes.
The dispersal of rodents and their fleas from the plague focus in Madagascar poses a serious potential health threat to other areas of Madagascar and to neighboring islands, such as Mayotte. Plague has not been reported in Mayotte, however, there is an important maritime trade route between Mayotte and north-west Madagascar. Our objectives were therefore: (i) to determine the genetic diversity and spatial genetic structure of X. cheopis populations in Madagascar; (ii) to determine the extent and pattern of any gene flow between X. cheopis populations in Madagascar and Mayotte, which may constitute an indication of plague introduction risk.
Samples of X. cheopis were collected in two areas of Mahajanga, Madagascar (the Port of Mahajanga and Marolaka) and the village of Longoni in Mayotte (Fig. 1). Additional specimens of X. cheopis were available from a further nine sites in Madagascar. All fleas were collected from small mammal hosts that had been trapped following published protocols . Some specimens in Madagascar were collected as part of human plague outbreak responses and/or rodent-vector monitoring (Table 1). Fleas were collected alive on their hosts and preserved in 70% ethanol. Two hundred and five fleas were collected (Table 1) and species identification was performed by an expert taxonomist. When available, more than 10 specimens per site were analysed genetically.
Isolation of microsatellite loci
Genomic DNA was extracted using a DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany) from 25 X. cheopis fleas collected in Madagascar. Partial genomic libraries enriched for microsatellite repeats were constructed using the protocol described by [36, 37] and summarized in Additional file 1. After an initial screening of 44 candidate microsatellite loci for suitability for genotyping (Additional file 1), a panel of 12 loci were deemed suitable for genotyping based on ease of scoring (e.g. no spurious PCR products, few stutter bands) and apparent polymorphism (Table 2).
DNA extraction, PCR and genotyping
Genomic DNA for genotyping was extracted using the Instagene™ Matrix kit (Bio-Rad Laboratories Inc., California, USA) following the manufacturer’s instructions. Multiplex PCR were performed by pooling microsatellite primers into three pools of 4 pairs of primers per pool (Table 2). Multiplex PCRs contained 7.5 μl of Type-it microsatellite PCR kit (2×) (Qiagen, Hilden, Germany), 0.4 μl of each forward and reverse primers of the pool (10 μM), 1 μl of DNA sample and sufficient water for a 15 μl final reaction volume. Thermal amplification conditions were 95 °C for 5 min, followed by 35 cycles of 95 °C for 30 s, 58 °C for 45 s, 72 °C for 55 s and a final extension step of 72 °C for 10 min. PCR products were sized by capillary electrophoresis on an ABI prism 3130 (Applied Biosystems, California, USA) with GeneScan™ 500 LIZ® Size Standard (Applied Biosystems, California, USA) and GENEMAPPER software (Applied Biosystems, California, USA).
Genotypic linkage disequilibrium between pairs of loci was tested with the exact test implemented in the GENEPOP v.4.3 . The same software was used to estimate basic measures of genetic variability: mean number of alleles per locus (Na), observed (Ho) and expected (He) heterozygosities and FIS  for each population. Significance test for FIS values were performed using Hardy-Weinberg exact test implemented in GENEPOP. Allelic data were checked for null alleles, allelic dropout and stutter bands using MICROCHECKER v.2.2.3 with the Oosterhout algorithm . Population differentiation (FST) was estimated using ARLEQUIN v.220.127.116.11 ; the ENA method was used also to calculate an unbiased FST (FST ENA) using FREENA .
Identification of genetic clusters
STRUCTURE simultaneously identifies potential populations (clusters) and probabilistically assigns individuals to each of the K populations based on the sample genotypes. STRUCTURE runs were performed using the admixture model and correlated allele frequencies, with one million iterations of the Markov Monte Carlo Chain (MCMC) used as ‘burn-in’ that were followed by ten million MCMC iterations; the probability to observe the data [Ln P(D)] was calculated for values of K ranging from 1 to 12, with five iterations for each K-value. The best estimate of K was taken to be the maximum value observed before the plateau of the curve Ln P(D) against K . STRUCTURE HARVESTER  was also used to identify the most pronounced level of population structure using the method of Evanno et al. . CLUMPP v.1.1.2  was used to find the optimal alignment from replicate STRUCTURE runs, with the summary of results generated using DISTRUCT v.1.1 .
DAPC is a multivariate statistical method which uses the k-means clustering approach . DAPC first transforms data using a principal components analysis (PCA) and subsequently identifies clusters using discriminant analysis (DA). This method thus defines genetic clusters, assigns individuals to clusters and allows visual assessment between-population differentiation. DAPC was implemented using the adegenet package  of the R software . The identification of the number of genetic clusters was done using the function find.clusters with the prior that the maximum number of clusters is equal to 12 (K = 12).
Inference of population introduction
An approximate Bayesian computation (ABC) analysis was conducted to infer the history of population introduction of X. cheopis between Mayotte and Madagascar. DIYABC v.2.0  was used to test seven different scenarios of introduction history (Fig. 2). Three populations were chosen to test these scenarios: one from Mayotte (the port of Longoni) and two from Madagascar (Marolaka located in the coastal region of Mahajanga and Amparafaravola located in the Central Highlands; Fig. 1). Amparafaravola was chosen as the site to represent the Central Highlands as this sample was the most representative of the existing genetic variability. For each scenario and each population, the following demographic parameters were estimated: dates of founding of the different populations (as number of generations) (ti), current effective population size (as number of diploid individuals, Ni), number of founders in the introduced populations (Nbi) and the duration of the initial bottleneck (dbi), which may be considered as a latency phase after each introduction event.
For all demographic parameters, prior distribution ranges were implemented according to current knowledge of X. cheopis. The generation time of X. cheopis varies between 18 days to 20 months according to biotic and abiotic conditions . Laboratory conditions of 25 ± 2 °C temperature and 70–80% relative humidity stimulate an average 12 generations per year. However, we note that natural populations of X. cheopis show seasonal dynamics in both the Central Highlands and Mahajanga . Larval and pupal development times and survival can be affected by temperature and humidity  and so, generation times may be longer under field conditions. As these conditions are similar to those in our study sites, time (expressed in generations before sampling) was translated into years assuming 12 generations per year. One million simulations were conducted under each scenario (a total of 7 million generations). Posterior probabilities of each scenario were computed by performing direct approach and a polychotomous weighted logistic regression on the 1% of simulated datasets closest to the observed dataset [54, 55] after linear discriminant analysis on summary statistics . Principal components analysis (PCA) was performed using summary statistics of simulated datasets and of observed dataset to check the suitability of the model (scenarios and prior parameters). Confidence in scenario choice was further tested using additional simulations. This included estimation of the probability of a type-II error (the probability of selecting the chosen scenario when it is not correct). After scenario choice, we proceeded to parameter inference estimated from the modes and 95% confidence intervals (CI) of their posterior distributions.
To determine directional relative migration between the three chosen populations, we used the divMigrate function from the R-package diveRsity v.1.9.89  using Jost’s D as measure of genetic distance  and with a bootstrap value of 1000. DivMigrate uses the method described in Sundqvist et al. .
Nine microsatellite markers produced PCR amplicons for all samples (XC009, XC037, XC021, SF009, XC013, XC023, XC044, XC016 and XC002; cf. Table 2) and only these nine loci were used in genetic analyses. The number of alleles at these 9 loci ranged from 4 to 11 (Table 2) and the mean number of alleles per locus and per populations (Na) ranged from 1.78 to 4.44 (Table 1). The expected heterozygosity (He) ranged from 0.25 to 0.50 while the observed heterozygosity (Ho) varied from 0.07 up to 0.41. All populations except those from Tsiroanomandidy, Toliara and Farafangana, showed a heterozygote deficit (i.e. He > Ho and thus deviation from expected Hardy-Weinberg equilibrium conditions), and with 7 of the 12 sites having significantly positive FIS values (Hardy-Weinberg exact test; P < 0.05) (Table 1). Only one of the 36 locus-pair combinations indicated significant linkage disequilibrium (XC002-XC009; P < 0.05). Estimated frequencies of null alleles per locus per population ranged from 0.16 to 0.45 (5 cases ≤ 0.20; 9 cases between 0.20 and 0.30; 2 cases between 0.30 and 0.40; and 1 case > 0.40) (Additional file 2: Table S1).
Both Bayesian clustering and multivariate methods identified three main genetic clusters: changes in Ln P(D) approached a plateau at K = 4 (Additional file 3: Figure S1), thus the best K value (the maximum K value observed before the plateau, ) was K = 3, the curve of DeltaK  showed a slope with a break occurring at K = 3 and the graphical output yielded by DAPC also supported three genetic clusters (Fig. 3).
Genetic clusters had some underlying geographical structure, whereby all individuals from Mayotte were assigned to one cluster and individuals from Madagascar were assigned to three clusters including the cluster that contained samples from Mayotte (Fig. 4). All fleas from Mayotte had a membership coefficient > 90% to cluster 1. Fleas from some regions in Madagascar, namely Marolaka, the port of Mahajanga, Mandritsara and Antananarivo, also had substantial membership to cluster 1 as well as a second cluster (indicated in red in Fig. 4) that was only identified in Madagascar. Most individuals from Toliara and Farafangana had high membership coefficients to cluster 2, while fleas from Ankasina and Ankazobe, had membership coefficients that indicated they belonged to cluster 2, but with more apparent genetic influence of the third cluster (indicated in yellow in Fig. 4). Fleas from Tsiroanomandidy and Ambondromisotra presented a membership coefficient greater than 90% to cluster 3; in Amparafaravola, cluster 3 was still dominant though some individuals appeared to belong to cluster 2 or cluster 1 (Fig. 4).
Division of populations into genetic groups was also reflected by the differentiation among these populations suggested by values of FST (and FST ENA; Additional file 4: Tables S2, S3). The FST values varied from 0 (Tsiroanomandidy vs Ambondromisotra) to 0.51 (Antananarivo vs Ambondromisotra), with most values quite high (> 0.25) and almost all values of FST significantly different from 0. The unbiased FST ENA values ranged from 0.01 (Port of Mahajanga vs Marolaka) to 0.45 (Antananarivo vs Ambondromisotra, Tsiroanomandidy vs Longoni and Toliara vs Tsiroanomandidy) (Additional file 4: Tables S2, S3). Comparable with estimates of FST, most pairwise estimates of FST ENA were relatively high (> 0.25).
Using the direct approach, the highest posterior probability (p) was obtained for scenario 6 (p = 0.50), with scenarios 7 and 5 having substantially lower probabilities (p = 0.19 and 0.18 respectively; Additional file 5: Table S4); the logistic approach altered slightly the ranking of demographic scenarios, marginally favoring Scenario 6 (P = 0.48) and with Scenario 5 favored over 7 (P = 0.43 and 0.06, respectively). All other scenarios for both approaches had low (< 0.1) probabilities (Additional file 5: Table S4 and Additional file 6: Figure S2). The type II error rate (Additional file 7: Table S5) associated with Scenario 6 (probability to select Scenario 6 though it is not correct) was relatively low (P = 0.08) providing confidence in the scenario choice. The favored Scenario 6 assumes that the population from Marolaka was a result of an admixture of a population from Longoni and from Amparafaravola (Fig. 2). DivMigrate analysis supported Scenario 6 and suggested that asymmetric migrations occurred between flea populations of Mayotte and Madagascar. Flea population sources from Longoni (Mayotte) and Amparafaravola (Madagascar) migrated and constituted the sink population of Marolaka (Madagascar). High significant relative migration rates were obtained: 0.87 (Longoni to Marolaka) and 1 (Amparafaravola to Marolaka) (Additional file 8: Table S6).
Estimates for the different parameters inferred using scenario 6 are given in Table 3. It was estimated that about 47 generations had occurred (95% CI: 16.4–241) between the time of sampling (t0) and the time of population admixture (t1) (Fig. 2; Scenario 6). Assuming 12 generations per year for X. cheopis and that the samples were collected in 2014, the population from Marolaka is estimated to have been affected by an introduction that occurred between 1993 and 2012 (t1: 2010 [95% CI: 1993–2012]) at an admixture rate of r2 = 55.8% (95% CI: 21–82%) (Table 3). If fewer generations per year occur under field conditions, a date towards the beginning of this time period is more likely.
Plague is a health problem in Madagascar and there is a need to understand the natural organization and dynamics of the key vector, Xenopsylla cheopis. This first population genetic analysis of X. cheopis from Madagascar and Mayotte makes important advances in our understanding: Xenopsylla cheopis populations are genetically and geographically structured in Madagascar, with interesting differences compared to previous studies of the principal host R. rattus and, a recent population migration of this flea occurred from Mayotte to Madagascar.
Spatial genetic structure of fleas in Madagascar
Within Madagascar, although there was fairly high genetic differentiation between populations, there was evidence of some gene flow between populations. As predicted, X. cheopis appears to show stronger population genetic structure than rats, with three clusters present in Madagascar compared to two for rats , and higher FST values (most values >0.25), compared with a range of 0.01–0.21 for rats . Although these patterns may suggest more limited dispersal in the parasite compared to the host, other factors such as lower effective population size may also play a role [31, 32]. Moreover, as several populations showed shared membership of different genetic clusters and the best-supported scenarios from the ABC analysis (Scenarios 6, 7 and 5) assume population admixture, whilst flea populations are genetically structured, there is also clear evidence of gene flow among populations.
Thus, comparisons of our results for X. cheopis with the results for R. rattus reveal some similarities and some interesting differences. It is notable that the R. rattus genetic structure appears to broadly follow a north/south divide . The northern genetic type of rats was largely restricted to the far north, with STRUCTURE analysis revealing only three sites out of 35 with evidence of mixed ancestry for some individuals. These three sites were geographically located between the northerly and southerly sites and included Mahajanga and a site close to Mandritsara, both sites in our study of X. cheopis which showed shared membership of different clusters. Although our results therefore suggest some congruence with the phylogeographic structure of hosts, as suggested for another fur flea, Listropsylla agrippinae, in South Africa , in general gene flow in X. cheopis did not appear to be strongly associated with geographical distance and R. rattus genetic structure. Specifically, some sites further south also showed shared membership of different clusters (Antananarivo, Amparafaravola), whilst, unlike for R. rattus , different flea populations within the Central Highlands could differ substantially in their proportion of membership to different clusters (e.g. Ankazobe vs Tsiroanomandidy and Ambondromisotra), indicating that in some cases little gene flow occurs between populations separated by relatively short distances. These apparently contradictory results could be explained by differences in the relative importance of human-mediated dispersal between R. rattus and X. cheopis and between sites.
A number of factors may lead to more frequent successful dispersal of fleas than black rats. Individual rats may carry a large number of fleas (in a study of rural highland villages in Madagascar, the average number of fleas on rats in houses was >2 ). Flea eggs, larvae and pupae are found in dust and debris from activities such as rice pounding , suggesting stages other than adults may also be potentially dispersed by human activity. Moreover, the competitive advantage of residents relative to migrants may be more significant for rats than fleas. Critically though, in urban sites that are linked by major roads and a high frequency of human movements, and are therefore more likely to experience human-mediated dispersal, R. rattus have largely been replaced by R. norvegicus (e.g. Mahajanga and Antananarivo ). Thus, at least in the last 30–40 years , human-mediated dispersal of X. cheopis may be much higher than for R. rattus. In contrast, dispersal of R. rattus and fleas in more remote, rural sites may be dominated by non-human mediated short-distance dispersal events. As R. rattus is abundant and widespread, occurring in diverse habitats , this would favor cumulative high levels of gene flow across rural landscapes. Indeed, in a population genetic study of rats at a more local landscape scale, rat populations showed only weak genetic differentiation between adjacent villages, with genetic structure at least partly related to topographic relief . In contrast, the restriction of X. cheopis to rats and other peridomestic mammals living inside houses is likely to mean that successful, non-human mediated dispersal between villages is relatively rare.
Whilst fleas did show more genetic structuring than rats, given the patchy distribution of X. cheopis in rural landscapes and the likely effects on dispersal rates and effective population sizes, the number of clusters detected in Madagascar (three) is perhaps surprisingly low. This is especially true when compared to the more significant, local clustering of Y. pestis , despite the more widespread distribution of Y. pestis in the landscape (transmitted between rats in habitats other than houses by the endemic flea S. fonquerniei) and presumably increased opportunities for dispersal. Although our results do suggest significant levels of gene flow, likely via human-mediated dispersal, they may also reflect the restricted nature of our study as sampling too few individuals per site can lead to underestimation of the number of clusters  and incorrect cluster assignment . However, simulation studies have indicated that reasonably limited sampling (6–10 individuals per site) can detect cryptic population structure and that any effects of low sample size are less with hierarchical population structure . In our study, some sites in the south had low sample size from few hosts (Ambondromisotra, Toliara and Farafangana). Despite this, whilst Ambondromisotra had high membership coefficient to cluster 3 and Farafangana to cluster 2, Toliara showed a mix of individuals with high membership coefficient to clusters 2 and 3.
Other factors may also influence the observed patterns of genetic variation within and between populations. Nine among twelve populations of X. cheopis showed heterozygote deficit which is significant in almost all cases. Evolutionary processes such as genetic drift, natural selection, inbreeding, mutation, population bottleneck or gene flow may influence allele frequencies and cause heterozygote deficit. Several of these may have contributed to our results. In addition to the gene admixture and gene flow between populations discussed above, fleas obtained from individual rats or rat families may be related, whilst population bottlenecks could occur within the spatially restricted populations of X. cheopis or due to seasonal dynamics in highland populations . Increasing the number of sites at both the national scale and landscape scale, and standardizing the sample size within sites, would undoubtedly clarify the population structure of X. cheopis within Madagascar and elucidate the role of rat dispersal and human-mediated dispersal. However, our rather limited study already yields important insights, possibly due to the strong hierarchical structure inherent within parasite populations.
Flea exchange between Mayotte and Madagascar
Our results indicate that limited flea exchange does occur between Mayotte and north Madagascar, with both the ABC analysis and the divMigrate analyses suggesting asymmetric gene flow from Mayotte to Madagascar. This exchange is likely to be linked to commerce and shipping routes in the Indian Ocean and reflects similar patterns in previous studies of other species. The mitochondrial haplotype group associated with the R. rattus introduction to north Madagascar  has previously been found in East Africa and Grande Comore, whilst the mitochondrial haplotypes associated with the southern R. rattus introduction was found in Mayotte . This led to the conclusion that Mayotte was colonized by R. rattus from Madagascar . Interestingly, although we find evidence of a relatively recent successful introduction of X. cheopis from Mayotte, only one R. rattus individual in Madagascar has been shown to carry a mitochondrial haplotype related to European human colonization , again possibly reflecting a greater success of flea migrants compared to rat migrants as discussed above. Successful introduction is conditioned by the introduction of some rats carrying fleas that are able to successfully reproduce.
Relevance of findings for plague introduction and dynamics
Our findings are important for assessing the epidemiological risk of plague introduction in plague free areas such as Mayotte. Mayotte already supports populations of both a plague-competent vector species (X. cheopis) and plague-competent mammal host species such as R. rattus and an insectivore found to be infected with Y. pestis in Mahajanga (Suncus murinus) . Although port cities such as Mahajanga or Longoni are hubs allowing plague introduction and intercontinental spread , the infrequent movement of fleas (and presumably their hosts) between the two ports indicated by this study suggest that the risk is real but limited. As there have been no plague outbreaks in humans in Mahajanga after the outbreaks in 1990s  and genetic studies of Y. pestis indicate that Mahajanga outbreaks were triggered by dispersal of infected rodents or fleas from the Central Highlands , the risk may have been further reduced because plague is no longer circulating in this region. However, studies of the peridomestic mammal and flea communities in Mahajanga are needed to assess this.
Even if plague did arrive in Mayotte due to the movement of infected fleas or hosts, many other factors would influence plague establishment and the risk of epidemics in the human population. One such factor is vector competence. Although the genetic differentiation among populations of X. cheopis is not associated with morphological differences, genetic differences might affect vector competence or resistance to pathogens  and therefore impact on disease transmission. Correlation between genetic structure of the insect vector and heterogeneity of vector competence has been reported for other disease-vector systems, for example Aedes albopictus (insect vector of dengue and yellow fever viruses) [65,66,67]. Interestingly, Tsiroanomandidy and Ambondromisotra, where cluster 3 occurs, are ‘plague focus’ areas where outbreaks occur every year. More work would need to be conducted to determine if different natural populations of X. cheopis differ in their vector competence.
Clearly, to further understand the risk of plague introduction and establishment in neighboring countries such as Mayotte as well as plague outbreaks in Madagascar outside the plague focus, there is a need for further, more extensive studies of flea and rat dynamics and dispersal, including for the endemic vector species, S. fonquerniei.
This study shows strong spatial structure among populations of the flea vector of plague, X. cheopis, from Madagascar and the nearby island of Mayotte. Gene flow occurred between Madagascar and Mayotte, but with evidence for a flea population having been introduced recently from Mayotte to Madagascar. As X. cheopis is the main vector of plague in Madagascar, the introduction of individuals to Mayotte may present a risk of plague introduction to this island.
Approximate Bayesian computation
Discriminant analysis of principal components
Excluding Null Allele
Markov Monte Carlo Chain
Principal components analysis
Rothschild NC. New species of Siphonaptera from Egypt and the Soudan. Entomol Mon Mag. 1903;39:83–7.
Gratz N. Rodent reservoirs and flea vectors of natural foci of plague. In: WHO Division of communicable diseases, editors. Plague manual: epidemiology, distribution, surveillance and control. Geneva, Switzerland: World Health Organization; 1999. p. 63–96.
Smit FGAM. The fleas of New Zealand (Siphonaptera). J Roy Soc New Zeal. 1979;9(2):143–232.
Beaucournu J-C, Fontenille D. Contribution à un catalogue des puces de Madagascar (Insecta, Siphonaptera). Arch Inst Pasteur Madagascar. 1991;59:57–98.
Klein J, Uilenberg G. Données faunistiques et écologiques sur les puces de Madagascar (Siphonaptera). Cah ORSTOM Sér Entomol Méd Parasitol. 1966;4(8):31–60.
Bitam I, Dittmar K, Parola P, Whiting MF, Raoult D. Fleas and flea-borne diseases. Int J of Infect Dis. 2010;14:e667–e76.
Eisen RJ, Gage KL. Transmission of flea-borne zoonotic agents. Annu Rev Entomol. 2012;57:61–82.
Eisen RJ, Wilder AP, Bearden SW, Montenieri JA, Gage KL. Early-phase transmission of Yersinia pestis by unblocked Xenopsylla cheopis (Siphonaptera: Pulicidae) is as efficient as transmission by blocked fleas. J Med Entomol. 2007;44:678–82.
Bertherat EG. Plague in Madagascar: overview of the 2014–2015 epidemic season. Week Epidemiol Rec. 2015;90(20):250–2.
Bertherat EG. Plague around the world, 2010–2015. Week Epidemiol Rec. 2016;91(8):89–93.
Brygoo ER. Epidemiologie de la peste à Madagascar. Arch Inst Pasteur Madagascar. 1966;35:9–147.
Morelli G, Song Y, Mazzoni CJ, Eppinger M, Roumagnac P, Wagner DM, et al. Yersinia pestis genome sequencing identifies patterns of global phylogenetic diversity. Nat Genet. 2010;42:1140–3.
Chanteau S, Ratsitorahina M, Rahalison L, Rasoamanana B, Chan F, Boisier P, et al. Current epidemiology of human plague in Madagascar. Microbes Infect. 2000;2(1):25–31.
Andrianaivoarimanana V, Kreppel K, Elissa N, Duplantier JM, Carniel E, Rajerison M, Jambou R. Understanding the persistence of plague foci in Madagascar. PLoS Negl Trop Dis. 2013;7:1–8.
Duplantier J-M, Duchemin JB, Ratsitorahina M, Rahalison L, Chanteau S. Résurgence de la peste dans le district d'Ikongo à Madagascar en 1998. 2. Réservoirs et vecteurs impliqués. Bull Soc Pathol Exot. 2001;94(2):119–22.
Migliani R, Ratsitorahina M, Rahalison L, Rakotoarivony I, Duchemin JB, Duplantier JM, et al. Résurgence de la peste dans le District d'Ikongo à Madagascar en 1998. 1. Aspects épidémiologiques dans la population humaine. Bull Soc Pathol Exot. 2001;94(2):115–8.
Richard V, Riehm JM, Herindrainy P, Soanandrasana R, Ratsitorahina M, Rakotomanana F, et al. Pneumonic plague outbreak, northern Madagascar, 2011. Emerg Infect Dis. 2015;21:8–15.
Rahelinirina S, Duplantier JM, Ratovonjato J, Ramilijaona O, Ratsimba M, Rahalison L. Study on the movement of Rattus rattus and evaluation of the plague dispersion in Madagascar. Vector Borne Zoonotic Dis. 2010;10(1):77–84.
Beaucournu J-C. Diversité des puces vectrices en fonction des foyers pesteux. Bull Soc Pathol Exot. 1999;92(5 Bis):419–21.
Brygoo ER, Dodin A. Transmission expérimentale de la peste par Synopsyllus fonquerniei de Madagascar. C R hebdomadaires Sci Acad Sci Série D: Sci Nat. 1967;264(4):681–2.
Duplantier J-M, Rakotondravony D. The rodent problem in Madagascar: agricultural pest and threat to human health. In: Singleton GR, Hinds LA, Leirs H, Zhang Z, editors. Ecologically-based rodent management. Canberra: ACIAR publications; 1999. p. 441–59.
Duplantier J-M, Duchemin J-B, Chanteau S, Carniel E. From the recent lessons of the Malagasy foci towards a global understanding of the factors involved in plague reemergence. Vet Res. 2005;36:437–53.
Brouat C, Tollenaere C, Estoup A, Loiseau A, Sommer S, Soanandrasana R, et al. Invasion genetics of a human commensal rodent: the black rat Rattus rattus in Madagascar. Mol Ecol. 2014;23:4153–67.
Vogler AJ, Chan F, Wagner DM, Roumagnac P, Lee J, Nera R, et al. Phylogeography and molecular epidemiology of Yersinia pestis in Madagascar. PLoS Negl Trop Dis. 2011;5(9):e1319.
Mears S, Clark F, Greenwood M, Larsen K. Host location, survival and fecundity of the oriental rat flea Xenopsylla cheopis (Siphonaptera: Pulicidae) in relation to black rat Rattus rattus (Rodentia: Muridae) host. Bull Entomol Res. 2002;92:375–84.
Boyer S, Miarinjara A, Elissa N. Xenopsylla cheopis (Siphonaptera: Pulicidae) susceptibility to deltamethrin in Madagascar. PLoS One. 2014;9(11):e111998.
Ratovonjato J, Duchemin JB, Duplantier J-M, Chanteau S. Xenopsylla cheopis (Siphonaptera: Xenopsyllinae), puces des foyers ruraux de peste des Hautes Terres malgaches: niveau de sensibilité au DDT, aux pyréthrinoïdes et aux carbamates après 50 années de lutte chimique. Arch Inst Pasteur Madagascar. 2000;66:9–12.
Zhang Y-C, Qi Y-M. Histological effects of fenoxycarb on the oriental rat flea, Xenopsylla cheopis (Siphonaptera: Pulicidae). Acta Entomol Sin. 2008;51(5):504–8.
Brinkerhoff RJ, Martin AP, Jones RT, Collinge SK. Population genetic structure of the prairie dog flea and plague vector, Oropsylla hirsuta. Parasitology. 2011;138(1):71–9.
Jones PH, Britten HB. The absence of concordant population genetic structure in the black-tailed prairie dog and the flea, Oropsylla hirsuta, with implications for the spread of Yersinia pestis. Mol Ecol. 2010;19:2038–49.
Huyse T, Poulin R, Theron A. Speciation in parasites: a population genetics approach. Trends Parasitol. 2005;21(10):469–75.
du Toit N, van Vuuren BJ, Matthee S, Matthee CA. Biogeography and host related factors trump parasite life history: limited congruence among the genetic structures of specific ectoparasitic lice and their rodent hosts. Mol Ecol. 2013;22(20):5185–204.
van der Mescht L, Matthee S, Matthee C. Comparative phylogeography between two generalist flea species reveal a complex interaction between parasite life history and host vicariance: parasite-host association matters. BMC Evol Biol. 2015;15:105.
Lin G, Zhao F, Chen H, Deng X, Su J, Zhang T. Comparative phylogeography of the plateau zokor (Eospalax baileyi) and its host-associated flea (Neopsylla paranoma) in the Qinghai-Tibet Plateau. BMC Evol Biol. 2014;14(180):2–9.
Rahelinirina S, Léon A, Harstskeerl RA, Sertour N, Ahmed A, Raharimanana C, et al. First isolation and direct evidence for the existence of large small-mammal reservoirs of Leptospira sp. in Madagascar. PloS One. 2010;5(11):e14111.
Gardner MG, Cooper SJB, Bull CM, Grant WN. Isolation of microsatellite loci from a social lizard, Egernia stokesii, using a modified enrichment procedure. J Hered. 1999;90(2):301–4.
Watts PC, Wu JH, Westgarth C, Thompson DJ, Kemp SJ. A panel of microsatellite loci for the southern damselfly, Coenagrion mercuriale (Odonata: Coenagrionidae). Conserv Genet. 2004;5(1):117–9.
Rousset F. GENEPOP'007: a complete re-implementation of the GENEPOP software for windows and Linux. Mol Ecol Res. 2008;8:103–6.
Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38(6):1358–70.
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.
Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and windows. Mol Ecol Res. 2010;10(3):564–7.
Chapuis M-P, Estoup A. Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol. 2007;24(3):621–31.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59.
Jombart T, Devillard S, Balloux F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 2010;11(1):94.
Earl DA, vonHoldt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Res. 2012;4(2):359–61.
Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14(8):2611–20.
Jakobsson M, Rosenberg NA. CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007;23(14):1801–6.
Rosenberg NA. DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes. 2004;4(1):137–8.
Jombart T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008;24(11):1403–5.
R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. 2008. https://www.r-project.org/. Accessed 2 April 2014.
Cornuet JM, Pudlo P, Veyssier J, Dehne-Garcia A, Gautier M, Leblois R, et al. DIYABC v2.0: a software to make approximate Bayesian computation inferences about population history using single nucleotide polymorphism, DNA sequence and microsatellite data. Bioinformatics. 2014;30:1187–9.
Rodhain F, Perez C. Précis d'entomologie médicale et vétérinaire. Paris: Maloine; 1985.
Kreppel KS, Telfer S, Rajerison M, Morse A, Baylis M. Effect of temperature and relative humidity on the development times and survival of Synopsyllus fonquerniei and Xenopsylla cheopis, the flea vectors of plague in Madagascar. Parasit Vectors. 2016;9(82):2–10.
Cornuet J-M, Ravigné V, Estoup A. Inference on population history and model checking using DNA sequence and microsatellite data with the software DIYABC (v1. 0). BMC Bioinform. 2010;11(1):401.
Cornuet J-M, Santos F, Beaumont MA, Robert CP, Marin J-M, Balding DJ, et al. Inferring population history with DIY ABC: a user-friendly approach to approximate Bayesian computation. Bioinformatics. 2008;24(23):2713–9.
Keenan K, McGinnity P, Cross TF, Crozier WW, Prodöhl PA. diveRsity: an R package for the estimation and exploration of population genetics parameters and their associated errors. Met Ecol Evol. 2013;4(8):782–8.
Jost L. GST and its relatives do not measure differentiation. Mol Ecol. 2008;17(18):4015–26.
Sundqvist L, Keenan K, Zackrisson M, Prodöhl P, Kleinhans D. Directional genetic differentiation and relative migration. Ecol Evol. 2016;6(11):3461–75.
Brouat C, Rahelinirina S, Loiseau A, Rahalison L, Rajerison M, Laffly D, et al. Plague circulation and population genetics of the reservoir Rattus rattus: the influence of topographic relief on the distribution of the disease within the Madagascan focus. PLoS Negl Trop Dis. 2013;7(6):e2266.
Fogelqvist J, Niittyvuopio A, Agren J, Savolainen O, Lascoux M. Cryptic population genetic structure: the number of inferred clusters depends on sample size. Mol Ecol Res. 2010;10:314–23.
Rosenberg NA, Burke T, Elo K, Feldman MW, Freidlin PJ, Groenen MAM, et al. Empirical evaluation of genetic clustering methods using multilocus genotypes from 20 chicken breeds. Genetics. 2001;159(2):699–713.
Tollenaere C, Brouat C, Duplantier J-M, Rahalison L, Rahelinirina S, Pascal M, et al. Phylogeography of the introduced species Rattus rattus in the western Indian Ocean, with special emphasis on the colonization history of Madagascar. J Biogeogr. 2010;37:398–410.
Vogler AJ, Chan F, Nottingham R, Andersen G, Drees K, Beckstrom-Sternberg SM, et al. A decade of plague in Mahajanga, Madagascar: insights into the global maritime spread of pandemic plague. MBio. 2013;4(1):e00623.
McCoy K. The population genetic structure of vectors and our understanding of disease epidemiology. Parasite. 2008;15:444–8.
Herrera F, Urdaneta L, Rivero J, Zoghbi N, Ruiz J, Carrasquel G, et al. Population genetic structure of the dengue mosquito Aedes aegypti in Venezuela. Mem Inst Oswaldo Cruz. 2006;101(6):625–33.
Paupy C, Chantha N, Vazeille M, Reynes J-M, Rodhain F, Failloux A-B. Variation over space and time of Aedes aegypti in Phnom Penh (Cambodia): genetic structure and oral susceptibility to a dengue virus. Genet Res. 2003;82(3):171–82.
Tien TK, Vazeille-Falcoz M, Mousson L, Hoang TH, Rodhain F, Huong NT, Failloux A-B. Aedes aegypti in ho chi Minh City (Viet Nam): susceptibility to dengue 2 virus and genetic differentiation. Trans R Soc Trop Med Hyg. 1999;93(6):581–6.
Thanks to the GIS team of the Institut Pasteur of Madagascar (IPM) for supplying base maps. Also, thanks to the IBIZA platform of the Mixed Research Unit (UMR PVBMT) of CIRAD in La Réunion for hosting genotyping lab works. Many thanks for Christophe Simiand for his help during genotyping lab-works.
Epidemiological surveys and laboratory works were funded by the Agence Régional de la Santé - Océan Indien (ARS-OI) (n°2/DSP/Etudes et Statistiques/2013) and the Institut Pasteur of Madagascar. Microsatellite primer development was supported by the Wellcome Trust (a part of Research Career Development Fellowship to S. Telfer #081705 and Senior Fellowship to S. Telfer #095171).
Availability of data and materials
The microsatellite datasets supporting the conclusions of this article are available in the GenBank repository (see Table 1 for accession numbers).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Protocol describing microsatellite primers development (partial genomic libraries enrichment). (DOCX 19 kb)
Estimated frequencies of null alleles per locus per population using Van Oosterhout’s method. (XLSX 10 kb)
Pairwise F ST and values between all populations. Table S3. Pairwise F ST ENA values among populations. (XLSX 12 kb)
Probability for each of the seven tested scenarios of population introduction using direct and logistic approaches. (XLSX 9 kb)
Raw results obtained for the seven scenarios tested using DIYABC software. (a) Posterior probabilities of scenarios obtained through a logistic regression computed every 10% (between 10 and 100%) of the number of selected datasets. (b) PCA plot allowing to visualize how close datasets simulated (each small dot) under each scenario (different colors) are from the observed dataset (large yellow dot). The most relevant scenario chosen was scenario 6. The probabilities for each scenario using direct and logistic approaches are given. (PDF 786 kb)
Confidence in scenario choice (using posterior based error computations). Pseudo-observed datasets (pods) were drawn from 500 simulated datasets closest to the observed dataset (s = 500). (DOCX 17 kb)
Bidirectional estimate of relative migration rates. The data give relative migration values by pairwise comparisons of the three populations from Madagascar and Mayotte. Analysis was performed using the function divMigrate from the package diveRsity . Jost’s D statistic was used to measure the genetic distance  and the bootstrap replicates were set to 1000. (XLSX 9 kb)
About this article
Cite this article
Harimalala, M., Telfer, S., Delatte, H. et al. Genetic structure and gene flow of the flea Xenopsylla cheopis in Madagascar and Mayotte. Parasites Vectors 10, 347 (2017). https://doi.org/10.1186/s13071-017-2290-6
- Xenopsylla cheopis
- Genetic structure
- Gene flow