Abstract
Fleas (Insecta, Siphonaptera) are the known vectors of serious bacterial pathogens, such as Yersinia pestis and Rickettsia typhi. The microbiota of fleas has been poorly investigated although it has a known influence on vector competence. Here, we report and analyse the microbiota of 577 flea specimens from Madagascar, a hotspot of plague transmission. Importantly, endemic Malagasy fleas show low host specificity, allowing addressing the importance of vertebrate host species in microbiota composition, as well as that of abiotic variables. We describe through Illumina sequencing of 2 hypervariable regions of 16 S rDNA the bacterial composition of 577 flea specimens of Madagascar. We address the importance of biotic (mammalian host and flea species) and abiotic (season and sampling site) variables on bacterial community composition. Bacterial composition appears driven by flea species and season, but interestingly not by the flea’s vertebrate hosts. These results suggest that the flea microbiota is at least in part acquired before they become adult, possibly in the immature off-host stages. Taken together, our results suggest that the microbiota of sampled fleas are composed of bacterial taxa with vertical transmission, such as Wolbachia which are prevalent in the present dataset, together with several bacterial taxa for which the occurrence is driven by environment factors, especially season and habitat. Given the importance of the microbiota in vector competence, we discuss the epidemiological consequences of environmentally-driven acquisition of microbiota in fleas on plague transmission in Madagascar.
Introduction
The holobiont concept has been proposed for invertebrate vectors, in which the host genome should not be considered as the only compartment at play in vector transmission1. It has been shown on several occasions that microparasites infecting vectors have major consequences on invertebrate biology and should thus considered as an intrinsic component of the vector. These microparasites include fungi, viruses and symbiotic bacteria, which all may drastically affect several life history traits of their invertebrate host, including their vectorial competence. For mostly technical reasons, the majority of investigations to date have explored the diversity of bacterial communities in arthropod vectors, as well as insights into community composition origin2,3. Such communities include core and non-core microbes, roughly equivalent to essential and non-essential microbes, respectively. The exploration of these aspects using bacterial communities has revealed that microbiota composition can be environmentally driven, and depend on a number of variables such as sampling site, season or vertebrate host species4. By contrast, some obligate symbionts are vertically transmitted and may hence display co-evolutionary patterns with their hosts. Indeed, obligate symbionts (also referred as mutualists) provide essential nutrients and vitamins to their arthropod hosts with sometimes extremely restricted diets such as bat flies5, ticks6 or aphids7. Endosymbionts such as vertically transmitted Wolbachia are well known sex manipulators that favour the fitness of infected females over that of males and uninfected females. These symbionts may in turn provide their hosts with protection against predators as established for aphids8 or viral infection9,10, which can be used to control the transmission of arboviruses by major mosquito vectors such as Aedes albopictus or Aedes aegypti11,12.
Fleas (Insecta, Siphonaptera) are a diversified group of hematophagous ectoparasites represented worldwide by over 2500 described species and subspecies13. These arthropods occupy a wide range of habitats and hosts, and are distributed unevenly across the globe14. Although the greatest diversity of fleas is found in the Palaearctic region15, fleas are found on all continents, including in extreme environments such as Antarctica where Glaciopsyllus antarcticus has been identified in seabirds16,17. Adult fleas feed parasitically on the blood of their vertebrate hosts, specifically mammals and birds. While fleas are rarely specific to their vertebrate host, some taxa associate with a particular host group, as exemplified by Xenopsylla spp., which are found specifically on rodents18. Fleas are known vectors of pathogens of significant medical importance such as the causative agents of plague (Yersinia pestis)19,20, murine typhus (Rickettsia typhi)21 or cat scratch disease (Bartonella henselae)22.
These vectors are of major medical importance on southwestern Indian Ocean islands, especially Madagascar, which reports among the highest incidence rates of plague worldwide23,24,25. In Madagascar, 49 flea species have been documented, of which 42 are endemic26. Malagasy fleas provide an interesting biological context for understanding the drivers of bacterial community composition in these disease vectors for two important reasons: (1) the unique endemism of the Malagasy fleas fauna allows exploring long-term evolution of their symbionts and (2) given the lack of reported host-specificity for endemic fleas from the Central Highlands, this area of the island provides the means to understand if vertebrate hosts are a strong driver of microbiota composition in fleas.
Using a sample composed of 12 species endemic to Madagascar sampled in different locations on the island, we addressed how biotic (i.e. flea and vertebrate host species) and abiotic (sampling site and season) variables influence the composition of the microbiota of fleas. The resulting data provide new insights into mechanisms driving the composition of bacterial communities in adult fleas, which may in turn have important consequences on the epidemiology of certain zoonotic diseases, especially plague.
Materials and methods
Flea specimens included in the analysis
Five hundred seventy-seven (221 from Ambohitantely, 312 from Ankazomivady and 44 from Lakato) adult fleas sampled either in Madagascar were included in the study. This sample is composed of 12 flea species endemic to Madagascar and representing 4 genera: Centetipsylla (C. madagascarensis), Paractenopsyllus (P. vauceli, P. duplantieri, P. grandidieri, P. petiti, P. rouxi and P. raxworthyi), Synopsyllus (S. robici, S. estradei and S. fonquerniei) and Dynopsyllus (D. brachypeten and D. flacourti). These specimens were previously sampled in the context of different research programs aiming at investigating the ecology of vector-borne microbes with zoonotic potential on islands from the southwestern Indian Ocean27,28. Malagasy fleas were collected on both introduced (M. musculus, R. rattus, S. murinus) and endemic Malagasy terrestrial mammals. The endemic Malagasy terrestrial mammals include 4 endemic rodent species of the subfamily Nesomyinae (Eliurus majori, E. minor, E. tanala and Nesomys rufus) and 10 endemic tenrecs (Hemicentetes nigriceps, Microgale cowani, M. dobsoni, M. longicaudata, M. gymnorhyncha, M. parvula, M. soricoides, M. taiva, Oryzorictes hova and Setifer setosus)29. A subsample from Ambohitantely was selected for exploring the impact of season on microbial composition as this site is located in a plague focus. This site was sampled both in winter (May 2014) and at the beginning of summer (October 2014), corresponding to the annual human plague season in Madagascar. The climatic aspects of seasonality at Ambohitantely include annual rainfall of about 1460 mm of which 86% falls between November and April, a cold season extending from May to August, and a warm season from December to February30.
DNA extraction procedures
Fleas included in the present study were obtained from previous studies27,28. As a result, no extraction controls could be included before Illumina sequencing. DNA extraction from Malagasy fleas was detailed previously28 and followed a non-destructive procedure maintaining the integrity of insect cuticles allowing morphological diagnosis and the production of molecular barcodes. Briefly, the individuals were placed on filter paper to remove the alcohol. Then, using a binocular, an incision was made in the abdomen with a syringe needle. Each individual was then placed in 100 µl of Instagate solution and PBS solution from the Instagene TM Matrix kit (Bio-Rad Laboratories, Hercules, California)28. DNA was pooled per location, flea species and host species. Sample size in each pool was composed of 2 to 52 DNA samples (see supplementary S1 Table for the composition of each pool).
Description of flea associated bacteriomes through Miseq sequencing of v3-v4 and v4-v5 regions of 16 S
In order to address the diversity of flea bacteriomes, both V3-V4 and V4-V5 regions of the bacterial 16SrRNA were PCR amplified and sequenced on an Illumina Miseq sequencing platform at Genoscreen (Lille, France). We used both the V3-V4 and V4-V5 regions of the bacterial 16SrDNA in order to limit PCR biases and maximize the detection of most bacterial diversity. The 3’ ends of V3-forward (TACGGRAGGCAGCAG) and V4-reverse (GGACTACCAGGGTATCTAAT) or V4-forward (GTGYCAGCMGCCGCGGTAA) and V5-reverse (CCGYCAATTYMTTTRAGTTT) bacterium-specific primers were associated at the 5’ end with multiplex identifier (MID) tags, a GsFLX key, and GsFLX adapters. Each pool was independently amplified twice with distinct MID tags, allowing the individual identification of each pool, as well as each of the two PCR duplicates. Sequences underwent a first quality control and analysis process before being demultiplexed, assembled in paired-end reads and depleted of primers sequences and low-quality reads. Both V3-V4 and V4-V5 raw datasets were first analysed using the SILVA database for the identification of OTUs. Then, V3-V4 and V4-V5 datasets were processed separately using a pipeline implemented in Mothur31. V3-V4 and V4-V5 datasets were depleted of sequences < 350 bp and < 300 bp, respectively. Sequences with ambiguous calls and homopolymers longer than 10 bp were removed. Chimeras were then removed and OTUs were defined by clustering at 97%. OTUs were then aligned on a customized SILVA database including the V3-V5 region of 16 S rRNA of bacteria.
Biodiversity analyses
Both V3-V4 and V4-V5 datasets, each consisting of 49 pools composed of fleas from a single species sampled on the same vertebrate host species, at the same study site and during the same sampling season were analysed independently using R software (R-4.2.2) in Rstudio (2022.07.2). The α-diversity was assessed using the vegan package, which allowed the evaluation of the Chao1 index, measuring the total number of unique taxa and the inverse of Simpson’s index estimating the diversity and evenness of communities. The structure of β-diversity was compared among flea genera and different sampling sites. For this, permutational MANOVA (PERMANOVA) tests with 999 permutations were performed and non-metric multidimensional scaling (NMDS) ordinations were conducted on Bray–Curtis dissimilarities, calculated from rarefied sequence counts, after square root transformation and Wisconsin standardization, to produce plots. Higher values for these indices indicate greater diversity (species richness and abundance). Measures of flea microbial diversity were examined in relation to flea species, host species, season, and collection site. The ANOVA test was used to examine significant differences in diversity measures and the test was considered significant when p < 0.05. Of the 49 pools, 23 corresponded to fleas collected in Réserve Spéciale d’Ambohitantely in the Central Highlands. This subsample was composed of fleas from five endemic flea species (Paractenopsyllus duplantieri, Paractenopsyllus grandidieri, Paractenopsyllus petiti, Synopsyllus estradei and Synopsyllus fonquerniei) sampled on 7 endemic (Microgale cowani, Microgale dobsoni, Microgale gymnorhyncha, Microgale longicaudata Microgale majori, Setifer setosus, Eliurus majori) and 1 introduced (Rattus rattus) mammal species. The two (V3-V4 and V4-V5) datasets from these 23 pools were eventually analysed independently as they allowed addressing the diversity of the flea microbiota at the same site and during two different seasons.
Results
Composition of fleas microbiota
Using the normalised data, a total of 1,380,057 and 1,588,722 reads was obtained for the V3-V4 and V4-V5 regions, respectively. The phyla Proteobacteria (V3-V4 = 88.2%; V4-V5 = 92.4%), Firmicutes (V3-V4 = 7.1%; V4-V5 = 5.1%) and Bacteroidetes (V3-V4 = 3.0%; V4-V5 = 2.5%) accounted for more than 96.0% of the sequences for each of the studied regions (see Supplementary Fig. 1). The analyses of V3-V4 reads revealed 223 bacterial genera of which 5 genera represented slightly more than 53% of all reads, namely Wolbachia (19.2%), Thorsellia (12.4%), Sphingomonas (8.5%), Rhizobium (6.8%) and Mesorhizobium (6.6%) (Fig. 1a). V4-V5 dataset revealed 196 bacterial genera of which 5 accounted for more than 66% of the reads namely Wolbachia (35.8%), Rickettsia (10.2%), Rhizobium (7.8%), Sphingomonas (6.5%) and Thorsellia (6.2%) (Fig. 1b).
Season but not host species influence bacterial communities of fleas
We investigated how flea species, vertebrate host species and sampling sites influenced bacterial community composition using the Chao1 and Simpson diversity indices (see supplementary S2 Table for all the indices). The Simpson diversity index was different according to flea species but this was significant for V4-V5 dataset only (Table 1). Importantly, for both the V3-V4 and V4-V5 data sets, no significant difference was detected in the bacterial community diversity when comparing a specific flea species sampled on different hosts (Table 1).
Table 1 Anova test P-values of Simpson Diversity and Chao1 Indices as a function of flea species, host species, being either endemic or introduced, collection site and season, for V3-V4 and V4-V5 datasets.
Full size table
We further addressed these patterns using a subsample composed of fleas from Ambohitantely, which were sampled at the same site during two successive seasons. Using this subset, an impact of the season was observed using both indices (Chao 1 and Simpson diversity indices) and, with bacterial diversity being much higher during the wet than the dry season, regardless of the 16 S sequenced region (Table 1).
Fig. 1
figure 1
Heatmap showing detected genera for which an abundance > 1% was detected in the entire dataset. V3-V4 (a) and V4-V5 (b) according to flea species. Abbreviations of flea genera include C = Centetipsylla, D = Dynopsyllus, P = Paractenopsyllus and S = Synopsyllus.
Full size image
Fig. 2
figure 2
Structure of bacterial communities depending on the flea genus using V3-V4 (a) and V4-V5 (b) data sets, and the sampling site using V3-V4 (c) and V4-V5 (d) data sets. NMDS plots were generated using Bray–Curtis distance matrices. The P-value of the PERMANOVA is indicated in the upper right corner.
Full size image
Importantly, this subsample confirmed that vertebrate hosts did not control bacterial community composition, using both indices and datasets (see Table 1, Ambo V3-V4 and Ambo V4-V5 columns). By contrast, the difference in bacterial diversity depending on flea species was detected using the Simpson index and V4-V5 dataset only (Table 1).
Lastly, using Permanova analyses on the whole sample, the bacterial community composition was significantly different between sampling sites for both V3-V4 (χ2 = 3.0653, df = 2, P < 0.001; Fig. 2c) and V4-V5 (χ2 = 3.417, df = 2, P < 0.01; Fig. 2d) datasets. No significant difference was detected for the bacterial community composition between flea genera and flea species for V3-V4 and V4-V5 datasets.
Discussion
The present study aimed at describing the bacterial community composition of a large diversity of endemic fleas from Madagascar and at investigating the biotic and abiotic variables driving the composition of these communities. We first explored the taxonomic composition of bacteria in fleas and then determined the effect of a number of biotic (flea genus and species, vertebrate host species) and abiotic (season, sampling site) variables on the composition of flea microbiota.
For both datasets, three phyla, namely Proteobacteria, Firmicutes and Bacteroides represented over 96.0% of the obtained sequences. These results are comparable to previous studies conducted on the flea microbiota in North America and Uganda32,33. At the genus level, both datasets showed that Wolbachia, Thorsellia, Sphingomonas, Rickettsia and environmental bacteria (Mesorhizobium and Rhizobium) were dominant and represented over 53% of reads.
It is important to note that the study was carried out using DNA of fleas that were extracted for other purposes27,28,34, and for which no extraction control was originally included. Therefore, presented data likely includes contaminant reads resulting from DNA extraction that we could not remove from raw sequence data. However, all fleas were extracted in the same way and in the same batch28, and contaminant DNA is thus likely homogenous throughout the samples. Therefore, although presented community composition includes possible extraction contaminants, we consider that differences in community composition depending on biotic/abiotic variable result from each analysed variable.
Our analyses showed that α-diversity of the microbiota is influenced by flea species but only for V4-V5 dataset. By contrast, vertebrate host species and sampling sites did not significantly influence the α-diversity of flea microbiota (see Table 1). In addition, Permanova analysis revealed that β-diversity was influenced by sampling site but not by flea species or genus. We then analysed independently fleas from Ambohitantely site in order to control for the effect of sampling site, address the importance of season in community composition and further test the importance of vertebrate host in driving bacterial community composition. The diversity of the flea microbiota at this locality appeared higher during the wet than during the dry season. Importantly, the bacterial community composition was not significantly affected by vertebrate host (see Table 1), confirming the pattern detected with the whole dataset. This strongly supports that bacterial community composition is driven by environmental factors as exemplified by the importance of season at this sampling site.
One of the most striking results of the present analyses is that the microbiota of fleas in general appeared influenced by non-biotic variables such as season but not by vertebrate host species, which is different from what is reported for other host-parasite systems, such as seabird ticks4. Indeed, ticks feed on their hosts at all development stages and the composition of their microbiota is largely influenced by their host species35,36. This pattern strongly suggests that most of the flea bacterial community is acquired by these invertebrates before parasitizing their vertebrate host. Fleas are characterized by a detritivore larval off-host stage37,38. Larvae feed on organic matter present in their environment (e.g. excrements of arthropod or vertebrate hosts). We propose that a fraction of the microbial diversity associated with this environment, including those in their diet, could be acquired by larval fleas and then trans-stadially transmitted to older stages.
The microbiota of arthropod vectors is increasingly recognized as controlling the vector capacity of their hosts, and is even used to block, for example, arboviruses transmitted by the mosquito Aedes aegypti. Given the medical importance of plague in Madagascar, it is of upmost importance to address whether some habitats could favour bacterial community structures conducive to the transmission of Yersinia pestis by these vectors. Bubonic plague transmission in Madagascar is largely restricted to the Central Highlands, with the notable exception of Mahajanga, a harbour town located along the northwest coast where Y. pestis maintenance is unstable, transmitted by fleas of the introduced shrew Suncus murinus, but recurrently introduced from the Central Highlands39. The distribution of bubonic plague in the Central Highlands is thought to be controlled by the distribution of the main vectors, i.e. introduced Xenopsylla cheopis and endemic Synopsylla fonquernei. Data presented herein strongly suggest that the environment shapes the microbial community of fleas, which might in turn have consequences on vector competence, in line with experimental infection carried out on fleas raised on different substrates32. Our data together with previously published investigations32,33 support that fleas microbiota is mostly structured by the environment, which has in turn consequences on Y. pestis transmission. It is thus important to address the role of the environment in plague epidemiology in Madagascar.
Data availability
Illumina data generated during the current study is accessible on Genbank, Sequence Read Archive (SRA) to PRJNA931254 (https://www.ncbi.nlm.nih.gov/sra/?term=PRJNA931254).
References
Rosenberg, E. & Zilber-Rosenberg, I. Symbiosis and development: The hologenome concept. Birth Defects Res. Part C Embryo. Today Rev. 93, 56–66 (2011).
ArticleMATHGoogle Scholar
Muturi, E. J., Ramirez, J. L., Rooney, A. P. & Kim, C.-H. Comparative analysis of gut microbiota of mosquito communities in central Illinois. PLoS Negl. Trop. Dis. 11, e0005377 (2017).
ArticlePubMedPubMed CentralGoogle Scholar
Gao, H., Cui, C., Wang, L., Jacobs-Lorena, M. & Wang, S. Mosquito microbiota and implications for disease control. Trends Parasitol. 36, 98–111 (2020).
ArticlePubMedMATHGoogle Scholar
Gomard, Y. et al. Changes in bacterial diversity, composition and interactions during the development of the seabird tick Ornithodoros maritimus (Argasidae). Microb. Ecol. 81, 770–783 (2021).
ArticleADSPubMedGoogle Scholar
Duron, O. et al. Origin, acquisition and diversification of heritable bacterial endosymbionts in louse flies and bat flies. Mol. Ecol. 23, 2105–2117 (2014).
ArticlePubMedGoogle Scholar
Duron, O. et al. Tick-bacteria mutualism depends on B vitamin synthesis pathways. Curr. Biol. 28, 1896-1902.e5 (2018).
ArticlePubMedGoogle Scholar
Wilson, A. C. C. et al. Genomic insight into the amino acid relations of the pea aphid, Acyrthosiphon pisum, with its symbiotic bacterium Buchnera aphidicola. Insect. Mol. Biol. 19(Suppl 2), 249–258 (2010).
ArticlePubMedGoogle Scholar
Tsuchida, T. et al. Symbiotic bacterium modifies aphid body color. Science 330, 1102–1104 (2010).
ArticleADSPubMedGoogle Scholar
Teixeira, L., Ferreira, A. & Ashburner, M. The bacterial symbiont Wolbachia induces resistance to RNA viral infections in Drosophila melanogaster. PLoS Biol. 6, e2 (2008).
ArticlePubMedMATHGoogle Scholar
Johnson, K. N. The impact of Wolbachia on virus infection in mosquitoes. Viruses 7, 5705–5717 (2015).
ArticlePubMedPubMed CentralMATHGoogle Scholar
Moreira, L. A. et al. A Wolbachia symbiont in Aedes aegypti limits infection with dengue, chikungunya, and plasmodium. Cell 139, 1268–1278 (2009).
ArticlePubMedGoogle Scholar
Dodson, B. L. et al. Wolbachia enhances West Nile Virus (WNV) infection in the mosquito Culex tarsalis. PLoS Negl. Trop. Dis. 8, e2965 (2014).
ArticlePubMedPubMed CentralGoogle Scholar
Whiting, M. F., Whiting, A. S., Hastriter, M. W. & Dittmar, K. A molecular phylogeny of fleas (Insecta: Siphonaptera): Origins and host associations. Cladistics 24, 677–707 (2008).
ArticleMATHGoogle Scholar
Medvedev, S. Specific features of the distribution and host associations of fleas (Siphonaptera). Entomol. Rev. 82, 1165–1177 (2002).
MATHGoogle Scholar
Krasnov, B. R. Functional and evolutionary ecology of fleas: A model for ecological parasitology (Cambridge University Press, 2008).
BookMATHGoogle Scholar
Bell, P. J., Burton, H. R. & van Franeker, J. A. Aspects of the biology of Glaciopsyllus antarcticus (Siphonaptera: Ceratophyllidae) during the breeding season of a host (Fulmarus glacialoides). Polar Biol. 8, 403–410 (1988).
ArticleGoogle Scholar
Steele, W., Pilgrim, R. & Palma, R. Occurrence the flea Glaciopsyllus antarcticus and avian lice in central Dronning Maud Land. Polar Biol. 18, 292–294 (1997).
ArticleGoogle Scholar
Bitam, I., Dittmar, K., Parola, P., Whiting, M. F. & Raoult, D. Fleas and flea-borne diseases. Int. J. Infect. Dis. 14, e667–e676 (2010).
ArticlePubMedGoogle Scholar
Hinnebusch, B. J. The evolution of flea-borne transmission in Yersinia pestis. Curr. Issues Mol. Biol. 7, 197–212 (2005).
PubMedMATHGoogle Scholar
Dennis, D. T. & Mead, P. S. Chapter 42 - Plague. In Trop Infect Dis 2nd edn (eds Guerrant, R. L. et al.) 471–481 (Churchill Livingstone, 2006).
ChapterGoogle Scholar
Azad, A. F. Epidemiology of murine typhus. Annu. Rev. Entomol. 35, 553–570 (1990).
ArticlePubMedMATHGoogle Scholar
Biancardi, A. L. & Curi, A. L. L. Cat-scratch disease. Ocul. Immunol. Inflamm. 22, 148–154 (2014).
ArticlePubMedGoogle Scholar
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. 36, 437–453 (2005).
ArticlePubMedMATHGoogle Scholar
Andrianaivoarimanana, V. et al. Understanding the persistence of plague foci in Madagascar. PLoS Negl. Trop. Dis. 7, e2382 (2013).
ArticlePubMedPubMed CentralGoogle Scholar
Alderson, J., Quastel, M., Wilson, E. & Bellamy, D. Factors influencing the re-emergence of plague in Madagascar. Emerg. Top. Life Sci. 4, 423–433 (2020).
ArticlePubMed CentralGoogle Scholar
Duchemin, J. B., Harimalala, M. & Duplantier, J.-M. Siphonaptera, fleas, parasy. New Nat Hist Madag. (Princeton: Princeton University Press; 2022). 1054–80.
Goodman, S. M., Randrenjarison Andriniaina, H. R., Soarimalala, V. & Beaucournu, J.-C. The fleas of endemic and introduced small mammals in central highland forests of Madagascar: Faunistics, species diversity, and absence of host specificity. J. Med. Entomol. 52, 1135–1143 (2015).
ArticlePubMedGoogle Scholar
Harimalala, M., Adélaïde, M., Duchemin, J.-B., Ramihangihajason, T. & Boyer, S. Species diversity and phylogeny of fleas of small terrestrial mammals in the forests of the Central Highlands of Madagascar. Zootaxa 4399, 181 (2018).
ArticlePubMedGoogle Scholar
Goodman, S. M. & Soarimalala, V. Introduction to mammals. New Nat Hist Madag. Steven M. Goodman. (Princeton: Princeton University Press; 2022). 1737–69.
Goodman, S. M., Raherilalao, M. J. & Wohlhauser, S., editors. Les aires protégées terrestres de Madagascar: Leur histoire, description et biota, tome 1: Introduction. (Association Vahatra in Antananarivo; 2021).
Schloss, P. D. et al. Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541 (2009).
ArticleADSPubMedPubMed CentralGoogle Scholar
Jones, R. T., Vetter, S. M. & Gage, K. L. Exposing laboratory-reared fleas to soil and wild flea feces increases transmission of Yersinia pestis. Am. J. Trop. Med. Hyg. 89, 784–787 (2013).
ArticlePubMedPubMed CentralGoogle Scholar
Jones, R. T., Borchert, J., Eisen, R., MacMillan, K., Boegler, K. & Gage, K. L. Flea-associated bacterial communities across an environmental transect in a plague-endemic region of Uganda. PLoS One. 10, (2015).
Guernier, V. et al. Fleas of small mammals on Reunion Island: Diversity, distribution and epidemiological consequences. PLoS Negl. Trop. Dis. 8, e3129 (2014).
ArticlePubMedPubMed CentralGoogle Scholar
Nakao, R. et al. A novel approach, based on BLSOMs (Batch Learning Self-Organizing Maps), to the microbiome analysis of ticks. ISME J. 7, 1003–1015 (2013).
ArticlePubMedPubMed CentralGoogle Scholar
Kueneman, J. G., Esser, H. J., Weiss, S. J., Jansen, P. A. & Foley, J. E. Tick microbiomes in neotropical forest fragments are best explained by tick-associated and environmental factors rather than host blood source. Appl. Environ. Microbiol. 87, e02668-e2720 (2021).
ADSPubMedPubMed CentralGoogle Scholar
Silverman, J. & Appel, A. G. Adult cat flea (Siphonaptera: Pulicidae) excretion of host blood proteins in relation to larval nutrition. J. Med. Entomol. 31, 265–271 (1994).
ArticlePubMedGoogle Scholar
O’Donnell, M. & Elston, D. M. What’s eating you? Human flea (Pulex irritans). Cutis. 106, (2020).
Vogler, A. J. et al. Temporal phylogeography of Yersinia pestis in Madagascar: Insights into the long-term maintenance of plague. PLoS Negl Trop Dis. 11, e0005887 (2017).
ArticlePubMedPubMed CentralGoogle Scholar
Download references
Acknowledgements
Authors acknowledge Muriel Dietrich for participating in fleas sampling in Madagascar and bringing advice on statistical analyses, and Sébastien Boyer for providing DNA samples from fleas collected in Madagascar by SMG and VS.
Funding
Funding was provided by a grant from the Grainger Bioinformatics Center at the Field Museum of Natural History and by the Programme d’Equipement Prioritaire et de Recherche Maladies Infectieuses Emergentes (PEPR MIE) Deciphering the Environmental, Biological and Societal factors driving the risk of Plague reemergence in France (DEBS-Plague).
Author information
Authors and Affiliations
Université de La Réunion, Unité Mixte de Recherche Processus Infectieux en Milieu Insulaire Tropical (UMR PIMIT), CNRS 9192, INSERM 1187, IRD 249, Plateforme Technologique CYROI, Sainte Clotilde, La Réunion, France
Victoria Carcauzon, Colette Cordonin & Pablo Tortosa
Université de Lille, INSERM 1019, CNRS 9017, CHU Lille, Institut Pasteur de Lille, Center for Infection and Immunity of Lille (CIIL), 59000, Lille, France
Florent Sebbane
Association Vahatra, BP 3972, 101, Antananarivo, Madagascar
Voahangy Soarimalala & Steven M. Goodman
Field Museum of Natural History, Chicago, IL, 60605-2496, USA
Steven M. Goodman
Authors
Victoria Carcauzon
View author publications
You can also search for this author in PubMedGoogle Scholar
2. Colette Cordonin
View author publications
You can also search for this author in PubMedGoogle Scholar
3. Florent Sebbane
View author publications
You can also search for this author in PubMedGoogle Scholar
4. Voahangy Soarimalala
View author publications
You can also search for this author in PubMedGoogle Scholar
5. Steven M. Goodman
View author publications
You can also search for this author in PubMedGoogle Scholar
6. Pablo Tortosa
View author publications
You can also search for this author in PubMedGoogle Scholar
Contributions
The project was initiated by CC, PT and SMG. All biological material from Madagascar was col-lected by MD, SMG and VS. VC, CC and PT analysed the data and wrote the initial draft of the article. FS participated to draft writing and made substantial contribution to the interpretation of the data.
Corresponding author
Correspondence to Victoria Carcauzon.
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Supplementary Material 1
Supplementary Material 2
Supplementary Material 3
Supplementary Material 4
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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-nc-nd/4.0/.
Reprints and permissions
About this article
Check for updates. Verify currency and authenticity via CrossMark
Cite this article
Carcauzon, V., Cordonin, C., Sebbane, F. et al. Malagasy flea microbiota results from a combination of vertically transmitted and environmentally acquired microbes. Sci Rep 15, 8461 (2025). https://doi.org/10.1038/s41598-025-90670-4
Download citation
Received:28 October 2024
Accepted:14 February 2025
Published:12 March 2025
DOI:https://doi.org/10.1038/s41598-025-90670-4
Share this article
Anyone you share the following link with will be able to read this content:
Get shareable link
Sorry, a shareable link is not currently available for this article.
Copy to clipboard
Provided by the Springer Nature SharedIt content-sharing initiative
Keywords
Plague
Fleas
Microbiota
Vector competence
Madagascar