Remarkable genetic diversity of Trypanosoma cruzi and Trypanosoma rangeli in two localities of southern Ecuador identified via deep sequencing of mini-exon gene amplicons

Whole-genome assembly of a hybrid Trypanosoma cruzi strain assembled with Nanopore sequencing alone

Trypanosoma cruzi is the causative agent of Chagas disease, which causes 10,000 deaths per year. Despite the high mortality associated with Chagas, relatively few parasite genomes have been assembled to date, with genome assemblies unavailable even for some commonly used laboratory strains. This is at least partially due to T. cruzi's highly complex and highly repetitive genome, which defies investigation using traditional short-read sequencing methods. In this study, we have generated a high-quality whole-genome assembly of the hybrid Tulahuen strain, a commercially available type VI strain, using long-read Nanopore sequencing without short-read scaffolding. The assembled genome contains 25% repeat regions, 17% variable multigene family members, and 27% transposable elements (TEs) and is of comparable quality with T. cruzi genome assemblies that utilized both long- and short-read data. Notably, we find that regions with TEs are significantly enriched for multicopy surface proteins, and that surface proteins are, on average, closer to TEs than to other coding regions. This finding suggests that mobile genetic elements such as transposons may drive recombination within surface protein gene families. This work demonstrates the feasibility of Nanopore sequencing to resolve complex regions of T. cruzi genomes, and with these resolved regions, provides support for a possible mechanism for genomic diversification.

Wide reference databases for typing Trypanosoma cruzi based on amplicon sequencing of the minicircle hypervariable region

BACKGROUND

Trypanosoma cruzi, the etiological agent of Chagas Disease, exhibits remarkable genetic diversity and is classified into different Discrete Typing Units (DTUs). Strain typing techniques are crucial for studying T. cruzi, because their DTUs have significant biological differences from one another. However, there is currently no methodological strategy for the direct typing of biological materials that has sufficient sensitivity, specificity, and reproducibility. The high diversity and copy number of the minicircle hypervariable regions (mHVRs) makes it a viable target for typing.

METHODOLOGY/PRINCIPAL FINDINGS

Approximately 24 million reads obtained by amplicon sequencing of the mHVR were analyzed for 62 strains belonging to the six main T. cruzi DTUs. To build reference databases of mHVR diversity for each DTU and to evaluate this target as a typing tool. Strains of the same DTU shared more mHVR clusters than strains of different DTUs, and clustered together. Different identity thresholds were used to build the reference sets of the mHVR sequences (85% and 95%, respectively). The 95% set had a higher specificity and was more suited for detecting co-infections, whereas the 85% set was excellent for identifying the primary DTU of a sample. The workflow's capacity for typing samples obtained from cultures, a set of whole-genome data, under various simulated PCR settings, in the presence of co-infecting lineages and for blood samples was also assessed.

CONCLUSIONS/SIGNIFICANCE

We present reference databases of mHVR sequences and an optimized typing workflow for T. cruzi including a simple online tool for deep amplicon sequencing analysis (https://ntomasini.github.io/cruzityping/). The results show that the workflow displays an equivalent resolution to that of the other typing methods. Owing to its specificity, sensitivity, relatively low cost, and simplicity, the proposed workflow could be an alternative for screening different types of samples.

Revisiting gene typing and phylogeny of Trypanosoma cruzi reference strains: Comparison of the relevance of mitochondrial DNA, single-copy nuclear DNA, and the intergenic region of mini-exon gene

Chagas disease is a widespread neglected disease in Latin America. Trypanosoma cruzi, the causative agent of the disease, is currently subdivided into six DTUs (discrete typing units) named TcI-TcVI, and although no clear association has been found between parasite genetics and different clinical outcomes of the disease or different transmission cycles, genetic characterization of T. cruzi strains remains crucial for integrated epidemiological studies. Numerous markers have been used for this purpose, although without consensus. These include mitochondrial genes, single or multiple-copy nuclear genes, ribosomal RNA genes, and the intergenic region of the repeated mini-exon gene. To increase our knowledge of these gene sequences and their usefulness for strain typing, we sequenced fragments of three mitochondrial genes, nine single-copy nuclear genes, and the repeated intergenic part of the mini-exon gene by Next Generation Sequencing (NGS) on a sample constituted of 16 strains representative of T. cruzi genetic diversity, to which we added the corresponding genetic data of the 38 T. cruzi genomes fully sequenced until 2022. Our results show that single-copy nuclear genes remain the gold standard for characterizing T. cruzi strains; the phylogenetic tree from concatenated genes (3959 bp) confirms the six DTUs previously recognized and provides additional information about the alleles present in the hybrid strains. In the tree built from the three mitochondrial concatenated genes (1274 bp), three main clusters are identified, including one with TcIII, TcIV, TcV, and TcVI DTUs which are not separated. Nevertheless, mitochondrial markers remain necessary for detecting introgression and heteroplasmy. The phylogenetic tree built from the sequence alignment of the repeated mini-exon gene fragment (327 bp) displayed six clusters, but only TcI was associated with a single cluster. The sequences obtained from strains belonging to the other DTUs were scattered into different clusters. Therefore, while the mini-exon marker may bring, for some biological samples, some advantages in terms of sensibility due to its repeated nature, mini-exon sequences must be used with caution and, when possible, avoided for T. cruzi typing and phylogenetic studies.

Amplicon Sequencing Reveals Complex Infection in Infants Congenitally Infected With Trypanosoma Cruzi and Informs the Dynamics of Parasite Transmission

Congenital transmission of Trypanosoma cruzi is an important source of new Chagas infections worldwide. The mechanisms of congenital transmission remain poorly understood, but there is evidence that parasite factors are involved. Investigating changes in parasite strain diversity during transmission could provide insight into the parasite factors that influence the process. Here we use amplicon sequencing of a single copy T. cruzi gene to evaluate the diversity of infection in clinical samples from Chagas positive mothers and their infected infants. Several infants and mothers were infected with multiple parasite strains, mostly of the same TcV lineage, and parasite strain diversity was higher in infants than mothers. Two parasite haplotypes were detected exclusively in infant samples, while one haplotype was never found in infants. Together, these data suggest multiple parasites initiate a congenital infection and that parasite factors influence the probability of vertical transmission.

Whole Genome Assembly of a Hybrid Trypanosoma cruzi Strain Assembled with Nanopore Sequencing Alone

Trypanosoma cruzi is the causative agent of Chagas disease, which causes 10,000 deaths per year. Despite the high mortality caused by the pathogen, relatively few parasite genomes have been assembled to date; even some commonly used laboratory strains do not have publicly available genome assemblies. This is at least partially due to T. cruzi's highly complex and highly repetitive genome: while describing the variation in genome content and structure is critical to better understanding T. cruzi biology and the mechanisms that underlie Chagas disease, the complexity of the genome defies investigation using traditional short read sequencing methods. Here, we have generated a high-quality whole genome assembly of the hybrid Tulahuen strain, a commercially available Type VI strain, using long read Nanopore sequencing without short read scaffolding. Using automated tools and manual curation for annotation, we report a genome with 25% repeat regions, 17% variable multigene family members, and 27% transposable elements. Notably, we find that regions with transposable elements are significantly enriched for surface proteins, and that on average surface proteins are closer to transposable elements compared to other coding regions. This finding supports a possible mechanism for diversification of surface proteins in which mobile genetic elements such as transposons facilitate recombination within the gene family. This work demonstrates the feasibility of nanopore sequencing to resolve complex regions of T. cruzi genomes, and with these resolved regions, provides support for a possible mechanism for genomic diversification.

An overview of the trypanosomatid (Kinetoplastida: Trypanosomatidae) parasites infecting several mammal species in Colombia

BACKGROUND

Trypanosomatids are among the most critical parasites for public health due to their impact on human, animal, and plant health. Diseases associated with these pathogens manifest mainly in poor and vulnerable populations, where social, environmental, and biological factors modulate the case incidence and geographical distribution.

METHODS

We used Sanger and amplicon-based next-generation sequencing (NGS) in samples from different mammals to identify trypanosomatid infections in several departments in Colombia. A total of 174 DNA samples (18 humans, 83 dogs, and 73 wild mammals) were analyzed by conventional PCR using a fragment of the heat shock protein 70 (Hsp70) gene and Sanger sequenced the positive samples. Twenty-seven samples were sent for amplicon-based NGS using the same gene fragment. Data obtained were used to perform diversity analyses.

RESULTS

One hundred and thirteen samples were positive for PCR by Hsp70 fragment; these corresponded to 22.1% Leishmania spp., 18.6% L. amazonensis, 9.7% L. braziliensis, 14.2% L. infantum, 8% L. panamensis, and 27.4% Trypanosoma cruzi. Comparison of the identified species by the two sequencing technologies used resulted in 97% concordance. Alpha and beta diversity indices were significant, mainly for dogs; there was an interesting index of coinfection events in the analyzed samples: different Leishmania species and the simultaneous presence of T. cruzi and even T. rangeli in one of the samples analyzed. Moreover, a low presence of L. braziliensis was observed in samples from wild mammals. Interestingly, to our knowledge, this is the first report of Leishmania detection in Hydrochaeris hydrochaeris (capybara) in Colombia.

CONCLUSIONS

The Hsp70 fragment used in this study is an optimal molecular marker for trypanosomatid identification in many hosts and allows the identification of different species in the same sample when amplicon-based sequencing is used. However, the use of this fragment for molecular diagnosis through conventional PCR should be carefully interpreted because of this same capacity to identify several parasites. This point is of pivotal importance in highly endemic countries across South America because of the co-circulation of different genera from the Trypanosomatidae family. The findings show an interesting starting point for One Health approaches in which coevolution and vector-host interactions can be studied.

Molecular detection of Trypanosoma (Trypanosomatidae) in bats from Thailand, with their phylogenetic relationships

The vast majority of trypanosome species is vector-borne parasites, with some of them being medically and veterinary important (such as Trypanosoma cruzi and Trypanosoma brucei) and capable of causing serious illness in vertebrate hosts. The discovery of trypanosomes in bats emphasizes the importance of bats as an important reservoir. Interestingly, there is a hypothesis that bats are ancestral hosts of T. cruzi. Trypanosome diversity has never been investigated in bats in Thailand, despite being in a biodiversity hot spot. To gain a better understanding of the diversity and evolutionary relationship of trypanosomes, polymerase chain reaction-based surveys were carried out from 2018 to 2020 in 17 sites. A total of 576 bats were captured, representing 23 species. A total of 38 (6.6%) positive samples was detected in ten bat species. Trypanosoma dionisii and Trypanosoma noyesi were identified from Myotis siligorensis and Megaderma spasma, respectively. The remaining 18S rRNA sequences of trypanosomes were related to other trypanosomes previously reported elsewhere. The sequences in the current study showed nucleotide identity as low as 90.74% compared to those of trypanosomes in the GenBank database, indicating the possibility of new species. All bat trypanosomes identified in the current study fall within the T. cruzi clade. The current study adds to evidence linking T. noyesi to a bat trypanosome and further supports the bat host origin of the T. cruzi clade. To the best of authors' knowledge, this is the first study on bat trypanosomes in Thailand and their phylogenetic relationships with global isolates.

Development of an Amplicon-Based Next-Generation Sequencing Protocol to Identify Leishmania Species and Other Trypanosomatids in Leishmaniasis Endemic Areas

Trypanosomatid infections are an important public health threat affecting many low-income countries across the tropics, particularly in the Americas. Trypanosomatids can infect many vertebrate, invertebrate, and plant species and play an important role as human pathogens. Among these clinically relevant pathogens are species from the genera Leishmania and Trypanosoma. Mixed trypanosomatid infections remain a largely unexplored phenomenon. Herein, we describe the application of an amplicon-based next-generation sequencing (NGS) assay to detect and identify trypanosomatid species in mammalian reservoirs, human patients, and sand fly vectors throughout regions of Leishmania endemicity. Sixty-five samples from different departments of Colombia, including two samples from Venezuela, were analyzed: 49 samples from cutaneous leishmaniasis (CL) patients, 8 from sand flies, 2 from domestic reservoirs (Canis familiaris), and 6 from wild reservoirs (Phyllostomus hastatus). DNA from each sample served to identify the presence of trypanosomatids through conventional PCR using heat shock protein 70 (HSP70) gene as the target. PCR products underwent sequencing by Sanger sequencing and NGS, and trypanosomatid species were identified by using BLASTn against a reference database built from trypanosomatid-derived HSP70 sequences. The alpha and beta diversity indexes of amplicon sequence variants were calculated for each group. The results revealed the presence of mixed infections with more than two Leishmania species in 34% of CL samples analyzed. Trypanosoma cruzi was identified in samples from wild reservoirs, as well as in sand fly vectors. Coinfection events with three different Leishmania species were identified in domestic reservoirs. These findings depose the traditional paradigm of leishmaniasis as being a single-species-driven infection and redraw the choreography of host-pathogen interaction in the context of multiparasitism. Further research is needed to decipher how coinfections may influence disease progression. This knowledge is key to developing an integrated approach for diagnosis and treatment.

IMPORTANCE Traditionally, there has been a frequent, yet incorrect assumption that phlebotomine vectors, animal reservoirs, and human hosts are susceptible to Leishmania infection by a single parasite species. However, current evidence supports that these new vector-parasite-reservoir associations lend vectors and reservoirs greater permissiveness to certain Leishmania species, thus promoting the appearance of coinfection events, particularly in disease-endemic regions. The present study describes the application of an amplicon-based next-generation sequencing (NGS) assay to detect and identify trypanosomatid species in mammalian reservoirs, human patients, and sand fly vectors from regions of endemicity for leishmaniasis. This changes our understanding of the clinical course of leishmaniasis in areas of endemicity.

Identification of Multiple Blastocystis Subtypes in Domestic Animals From Colombia Using Amplicon-Based Next Generation Sequencing

Blastocystis is frequently reported in fecal samples from animals and humans worldwide, and a variety of subtypes (STs) have been observed in wild and domestic animals. In Colombia, few studies have focused on the transmission dynamics and epidemiological importance of Blastocystis in animals. In this study, we characterized the frequency and subtypes of Blastocystis in fecal samples of domestic animals including pigs, minipigs, cows, dogs, horses, goats, sheep, and llama from three departments of Colombia. Of the 118 fecal samples included in this study 81.4% (n = 96) were positive for Blastocystis using a PCR that amplifies a fragment of the small subunit ribosomal RNA (SSU rRNA) gene. PCR positive samples were sequenced by next generation amplicon sequencing (NGS) to determine subtypes. Eleven subtypes were detected, ten previously reported, ST5 (50.7%), ST10 (47.8%), ST25 (34.3%), ST26 (29.8%), ST21 (22.4%), ST23 (22.4%), ST1 (17.9%), ST14 (16.4%), ST24 (14.9%), ST3 (7.5%), and a novel subtype, named ST32 (3.0%). Mixed infection and/or intra -subtype variations were identified in most of the samples. Novel ST32 was observed in two samples from a goat and a cow. To support novel subtype designation, a MinION based sequencing strategy was used to generate the full-length of the SSU rRNA gene. Comparison of full-length nucleotide sequences with those from current valid subtypes supported the designation of ST32. This is the first study in Colombia using NGS to molecularly characterize subtypes of Blastocystis in farm animals. A great diversity of subtypes was observed in domestic animals including subtypes previously identified in humans. Additionally, subtype overlap between the different hosts examined in this study were observed. These findings highlight the presence of Blastocystis subtypes with zoonotic potential in farm animals indicating that farm animals could play a role in transmission to humans.