Small subunit ribosomal metabarcoding reveals extraordinary trypanosomatid diversity in Brazilian bats

Background: Bats are a highly successful, globally dispersed order of mammals that occupy a wide array of ecological niches. They are also intensely parasitized and implicated in multiple viral, bacterial and parasitic zoonoses. Trypanosomes are thought to be especially abundant and diverse in bats. In this study, we used 18S ribosomal RNA metabarcoding to probe bat trypanosome diversity in unprecedented detail. Methodology/Principal Findings: Total DNA was extracted from the blood of 90 bat individuals (17 species) captured along Atlantic Forest fragments of Espírito Santo state, southeast Brazil. 18S ribosomal RNA was amplified by standard and/or nested PCR, then deep sequenced to recover and identify Operational Taxonomic Units (OTUs) for phylogenetic analysis. Blood samples from 34 bat individuals (13 species) tested positive for infection by 18S rRNA amplification. Amplicon sequences clustered to 14 OTUs, of which five were identified as Trypanosoma cruzi I, T. cruzi III/V, Trypanosoma cruzi marinkellei, Trypanosoma rangeli, and Trypanosoma dionisii, and seven were identified as novel genotypes monophyletic to basal T. cruzi clade types of the New World. Another OTU was identified as a trypanosome like those found in reptiles. Surprisingly, the remaining OTU was identified as Bodo saltans–closest non-parasitic relative of the trypanosomatid order. While three blood samples featured just one OTU (T. dionisii), all others resolved as mixed infections of up to eight OTUs. Conclusions/Significance: This study demonstrates the utility of next-generation barcoding methods to screen parasite diversity in mammalian reservoir hosts. We exposed high rates of local bat parasitism by multiple trypanosome species, some known to cause fatal human disease, others non-pathogenic, novel or yet little understood. Our results highlight bats as a long-standing nexus among host-parasite interactions of multiple niches, sustained in part by opportunistic and incidental infections of consequence to evolutionary theory as much as to public health. Author summary: Bats make up a mega-diverse, intensely parasitized order of volant mammals whose unique behavioural and physiological adaptations promote infection by a vast array of microorganisms. Trypanosomes stand out as ancient protozoan parasites of bats. As cryptic morphology, low parasitaemia and selective growth in culture have recurrently biased survey, we used 18S ribosomal RNA metabarcoding to resolve bat trypanosomatid diversity in Atlantic Forest fragments of southeast Brazil. Next to several unknown species, our deep sequence-based detection and assignment protocol recognized multiple known human-pathogenic trypanosomes, another linked to reptile hosts as well as a non-parasitic kinetoplastid in the blood of various phyllostomid bats. The striking permissivity exposed here, in a region where bat trypanosomes recently featured in a fatal case of Chagas disease, compels further research on bats’ role in the dispersal and spill-over of various microorganisms among humans and wildlife

Leishmania donovani populations in Eastern Sudan: temporal structuring and a link between human and canine transmission.

BACKGROUND: Visceral leishmaniasis (VL), caused by the members of the Leishmania donovani complex, has been responsible for devastating VL epidemics in the Sudan. Multilocus microsatellite and sequence typing studies can provide valuable insights into the molecular epidemiology of leishmaniasis, when applied at local scales. Here we present population genetic data for a large panel of strains and clones collected in endemic Sudan between 1993 and 2001. METHODS: Genetic diversity was evaluated at fourteen microsatellite markers and eleven nuclear sequence loci across 124 strains and clones. RESULTS: Microsatellite data defined six genetic subpopulations with which the nuclear sequence data were broadly congruent. Pairwise estimates of FST (microsatellite) and KST (sequence) indicated small but significant shifts among the allelic repertoires of circulating strains year on year. Furthermore, we noted the co-occurrence of human and canine L. donovani strains in three of the six clusters defined. Finally, we identified widespread deficit in heterozygosity in all four years tested but strong deviation from inter-locus linkage equilibrium in two years. CONCLUSIONS: Significant genetic diversity is present among L. donovani in Sudan, and minor population structuring between years is characteristic of entrenched, endemic disease transmission. Seasonality in vector abundance and transmission may, to an extent, explain the shallow temporal clines in allelic frequency that we observed. Genetically similar canine and human strains highlight the role of dogs as important local reservoirs of visceral leishmaniasis

Ecological host fitting of Trypanosoma cruzi TcI in Bolivia: mosaic population structure, hybridization and a role for humans in Andean parasite dispersal.

An improved understanding of how a parasite species exploits its genetic repertoire to colonize novel hosts and environmental niches is crucial to establish the epidemiological risk associated with emergent pathogenic genotypes. Trypanosoma cruzi, a genetically heterogeneous, multi-host zoonosis, provides an ideal system to examine the sylvatic diversification of parasitic protozoa. In Bolivia, T. cruzi I, the oldest and most widespread genetic lineage, is pervasive across a range of ecological clines. High-resolution nuclear (26 loci) and mitochondrial (10 loci) genotyping of 199 contemporaneous sylvatic TcI clones was undertaken to provide insights into the biogeographical basis of T. cruzi evolution. Three distinct sylvatic parasite transmission cycles were identified: one highland population among terrestrial rodent and triatomine species, composed of genetically homogenous strains (Ar = 2.95; PA/L = 0.61; DAS = 0.151), and two highly diverse, parasite assemblages circulating among predominantly arboreal mammals and vectors in the lowlands (Ar = 3.40 and 3.93; PA/L = 1.12 and 0.60; DAS = 0.425 and 0.311, respectively). Very limited gene flow between neighbouring terrestrial highland and arboreal lowland areas (distance ~220 km; FST = 0.42 and 0.35) but strong connectivity between ecologically similar but geographically disparate terrestrial highland ecotopes (distance >465 km; FST = 0.016-0.084) strongly supports ecological host fitting as the predominant mechanism of parasite diversification. Dissimilar heterozygosity estimates (excess in highlands, deficit in lowlands) and mitochondrial introgression among lowland strains may indicate fundamental differences in mating strategies between populations. Finally, accelerated parasite dissemination between densely populated, highland areas, compared to uninhabited lowland foci, likely reflects passive, long-range anthroponotic dispersal. The impact of humans on the risk of epizootic Chagas disease transmission in Bolivia is discussed

Wild Trypanosoma cruzi I genetic diversity in Brazil suggests admixture and disturbance in parasite populations from the Atlantic Forest region.

BACKGROUND: Trypanosoma cruzi (Kinetoplastida, Trypanosomatidae) infection is an ancient and widespread zoonosis distributed throughout the Americas. Ecologically, Brazil comprises several distinct biomes: Amazonia, Cerrado, Caatinga, Pantanal and the Atlantic Forest. Sylvatic T. cruzi transmission is known to occur throughout these biomes, with multiple hosts and vectors involved. Parasite species-level genetic diversity can be a useful marker for ecosystem health. Our aims were to: investigate sylvatic T. cruzi genetic diversity across different biomes, detect instances of genetic exchange, and explore the possible impact of ecological disturbance on parasite diversity at an intra-species level. METHODS: We characterised 107 isolates of T. cruzi I (TcI; discrete typing unit, DTU I) from different major Brazilian biomes with twenty-seven nuclear microsatellite loci. A representative subset of biologically cloned isolates was further characterised using ten mitochondrial gene loci. We compared these data generated from Brazilian TcI isolates from around America. RESULTS: Genetic diversity was remarkably high, including one divergent cluster that branched outside the known genetic diversity of TcI in the Americas. We detected evidence for mitochondrial introgression and genetic exchange between the eastern Amazon and Caatinga. Finally, we found strong signatures of admixture among isolates from the Atlantic Forest region by comparison to parasites from other study sites. CONCLUSIONS: Atlantic Forest sylvatic TcI populations are highly fragmented and admixed by comparison to others around Brazil. We speculate on: the possible causes of Atlantic Forest admixture; the role of T. cruzi as a sentinel for ecosystem health, and the impact disrupted sylvatic transmission cycles might have on accurate source attribution in oral outbreaks

Trypanosoma cruzi IIc: phylogenetic and phylogeographic insights from sequence and microsatellite analysis and potential impact on emergent Chagas disease.

Trypanosoma cruzi, the etiological agent of Chagas disease, is highly genetically diverse. Numerous lines of evidence point to the existence of six stable genetic lineages or DTUs: TcI, TcIIa, TcIIb, TcIIc, TcIId, and TcIIe. Molecular dating suggests that T. cruzi is likely to have been an endemic infection of neotropical mammalian fauna for many millions of years. Here we have applied a panel of 49 polymorphic microsatellite markers developed from the online T. cruzi genome to document genetic diversity among 53 isolates belonging to TcIIc, a lineage so far recorded almost exclusively in silvatic transmission cycles but increasingly a potential source of human infection. These data are complemented by parallel analysis of sequence variation in a fragment of the glucose-6-phosphate isomerase gene. New isolates confirm that TcIIc is associated with terrestrial transmission cycles and armadillo reservoir hosts, and demonstrate that TcIIc is far more widespread than previously thought, with a distribution at least from Western Venezuela to the Argentine Chaco. We show that TcIIc is truly a discrete T. cruzi lineage, that it could have an ancient origin and that diversity occurs within the terrestrial niche independently of the host species. We also show that spatial structure among TcIIc isolates from its principal host, the armadillo Dasypus novemcinctus, is greater than that among TcI from Didelphis spp. opossums and link this observation to differences in ecology of their respective niches. Homozygosity in TcIIc populations and some linkage indices indicate the possibility of recombination but cannot yet be effectively discriminated from a high genome-wide frequency of gene conversion. Finally, we suggest that the derived TcIIc population genetic data have a vital role in determining the origin of the epidemiologically important hybrid lineages TcIId and TcIIe

Importation of Hybrid Human-Associated Trypanosoma cruzi Strains of Southern South American Origin, Colombia.

We report the characterization of Trypanosoma cruzi of southern South American origin among humans, domestic vectors, and peridomestic hosts in Colombia using high-resolution nuclear and mitochondrial genotyping. Expanding our understanding of the geographic range of lineage TcVI, which is associated with severe Chagas disease, will help clarify risk of human infection for improved disease control

2b-RAD genotyping for population genomic studies of Chagas disease vectors: Rhodnius ecuadoriensis in Ecuador

Background: Rhodnius ecuadoriensis is the main triatomine vector of Chagas disease, American trypanosomiasis, in Southern Ecuador and Northern Peru. Genomic approaches and next generation sequencing technologies have become powerful tools for investigating population diversity and structure which is a key consideration for vector control. Here we assess the effectiveness of three different 2b restriction site-associated DNA (2b-RAD) genotyping strategies in R. ecuadoriensis to provide sufficient genomic resolution to tease apart microevolutionary processes and undertake some pilot population genomic analyses. Methodology/Principal findings: The 2b-RAD protocol was carried out in-house at a non-specialized laboratory using 20 R. ecuadoriensis adults collected from the central coast and southern Andean region of Ecuador, from June 2006 to July 2013. 2b-RAD sequencing data was performed on an Illumina MiSeq instrument and analyzed with the STACKS de novo pipeline for loci assembly and Single Nucleotide Polymorphism (SNP) discovery. Preliminary population genomic analyses (global AMOVA and Bayesian clustering) were implemented. Our results showed that the 2b-RAD genotyping protocol is effective for R. ecuadoriensis and likely for other triatomine species. However, only BcgI and CspCI restriction enzymes provided a number of markers suitable for population genomic analysis at the read depth we generated. Our preliminary genomic analyses detected a signal of genetic structuring across the study area. Conclusions/Significance: Our findings suggest that 2b-RAD genotyping is both a cost effective and methodologically simple approach for generating high resolution genomic data for Chagas disease vectors with the power to distinguish between different vector populations at epidemiologically relevant scales. As such, 2b-RAD represents a powerful tool in the hands of medical entomologists with limited access to specialized molecular biological equipment. Author summary: Understanding Chagas disease vector (triatomine) population dispersal is key for the design of control measures tailored for the epidemiological situation of a particular region. In Ecuador, Rhodnius ecuadoriensis is a cause of concern for Chagas disease transmission, since it is widely distributed from the central coast to southern Ecuador. Here, a genome-wide sequencing (2b-RAD) approach was performed in 20 specimens from four communities from Manabí (central coast) and Loja (southern) provinces of Ecuador, and the effectiveness of three type IIB restriction enzymes was assessed. The findings of this study show that this genotyping methodology is cost effective in R. ecuadoriensis and likely in other triatomine species. In addition, preliminary population genomic analysis results detected a signal of population structure among geographically distinct communities and genetic variability within communities. As such, 2b-RAD shows significant promise as a relatively low-tech solution for determination of vector population genomics, dynamics, and spread

The potential role of gut microbiota in shaping host energetics and metabolic rate

It is increasingly recognised that symbiotic microbiota (especially those present in the gut) have important influences on the functioning of their host. Here we review the interplay between this microbial community and the growth, metabolic rate and nutritional energy harvest of the host. We show how recent developments in experimental and analytical methods have allowed much easier characterisation of the nature, and increasingly the functioning, of the gut microbiota. Manipulation studies that remove or augment gut microorganisms or transfer them between hosts have allowed unprecedented insights into their impact. While much of the information to date has come from studies of laboratory model organisms, recent studies have used a more diverse range of host species, including those living in natural conditions, revealing their ecological relevance. The gut microbiota can provide the host with dietary nutrients that would be otherwise unobtainable, as well as allow the host flexibility in its capacity to cope with changing environments. The composition of the gut microbial community of a species can vary seasonally or when the host moves between environments (e.g. fresh and sea water in the case of migratory fish). It can also change with host diet choice, metabolic rate (or demands) and life stage. These changes in gut microbial community composition enable the host to live within different environments, adapt to seasonal changes in diet and maintain performance throughout its entire life history, highlighting the ecological relevance of the gut microbiota. While it is evident that gut microbes can underpin host metabolic plasticity, the causal nature of associations between particular microorganisms and host performance is not always clear unless a manipulative approach has been used. Many studies have focussed on a correlative approach by characterising microbial community composition, but there is now a need for more experimental studies in both wild and laboratory‐based environments, to reveal the true role of gut microbiota in influencing the functioning of their hosts, including its capacity to tolerate environmental change. We highlight areas where these would be particularly fruitful in the context of ecological energetics

temporal structuring and a link between human and canine transmission

Background Visceral leishmaniasis (VL), caused by the members of the Leishmania donovani complex, has been responsible for devastating VL epidemics in the Sudan. Multilocus microsatellite and sequence typing studies can provide valuable insights into the molecular epidemiology of leishmaniasis, when applied at local scales. Here we present population genetic data for a large panel of strains and clones collected in endemic Sudan between 1993 and 2001. Methods Genetic diversity was evaluated at fourteen microsatellite markers and eleven nuclear sequence loci across 124 strains and clones. Results Microsatellite data defined six genetic subpopulations with which the nuclear sequence data were broadly congruent. Pairwise estimates of FST (microsatellite) and KST (sequence) indicated small but significant shifts among the allelic repertoires of circulating strains year on year. Furthermore, we noted the co-occurrence of human and canine L. donovani strains in three of the six clusters defined. Finally, we identified widespread deficit in heterozygosity in all four years tested but strong deviation from inter-locus linkage equilibrium in two years. Conclusions Significant genetic diversity is present among L. donovani in Sudan, and minor population structuring between years is characteristic of entrenched, endemic disease transmission. Seasonality in vector abundance and transmission may, to an extent, explain the shallow temporal clines in allelic frequency that we observed. Genetically similar canine and human strains highlight the role of dogs as important local reservoirs of visceral leishmaniasis

High Trypanosoma spp. diversity is maintained by bats and triatomines in Espírito Santo state, Brazil

The aim of this study was to reevaluate the ecology of an area in the Atlantic Forest, southeast Brazil, where Chagas disease (CD) has been found to occur. In a previous study, immediately after the occurrence of a CD case, we did not observe any sylvatic small mammals or dogs with Trypanosoma cruzi cruzi infections, but Triatoma vitticeps presented high T. c. cruzi infection rates. In this study, we investigated bats together with non-volant mammals, dogs, and triatomines to explore other possible T. c. cruzi reservoirs/hosts in the area. Seventy-three non-volant mammals and 186 bats were captured at three sites within the Guarapari municipality, Espírito Santo state. Rio da Prata and Amarelos sites exhibited greater richness in terms of non-volant mammals and bats species, respectively. The marsupial Metachirus nudicaudatus, the rodent Trinomys paratus, and the bats Artibeus lituratus and Carollia perspicillata were the most frequently captured species. As determined by positive hemocultures, only two non-volant mammals were found to be infected by Trypanosoma species: Monodelphis americana, which was infected by T. cascavelli, T. dionisii and Trypanosoma sp., and Callithrix geoffroyi, which was infected by T. minasense. Bats presented T. c. cruzi TcI and TcIII/V, T. c. marinkellei, T. dionisii, T. rangeli B and D, and Trypanosoma sp. infections. Seven dogs were infected with T. cruzi based only on serological exams. The triatomines T. vitticeps and Panstrongylus geniculatus were found to be infected by trypanosomes via microscopy. According to molecular characterization, T. vitticeps specimens were infected with T. c. cruzi TcI, TcII, TcIII/V, and TcIV, T. c. marinkellei and T. dionisii. We observed high trypanosome diversity in a small and fragmented region of the Atlantic Forest. This diversity was primarily maintained by bats and T. vitticeps. Our findings show that the host specificity of the Trypanosoma genus should be thoroughly reviewed. In addition, our data show that CD cases can occur without an enzootic cycle near residential areas