Optimized multilocus sequence typing (MLST) scheme for Trypanosoma cruzi

New insights into Trypanosoma cruzi genetic diversity, and its influence on parasite biology and clinical outcomes

Chagas disease, caused by Trypanosoma cruzi, remains a serious public health problem worldwide. The parasite was subdivided into six distinct genetic groups, called "discrete typing units" (DTUs), from TcI to TcVI. Several studies have indicated that the heterogeneity of T. cruzi species directly affects the diversity of clinical manifestations of Chagas disease, control, diagnosis performance, and susceptibility to treatment. Thus, this review aims to describe how T. cruzi genetic diversity influences the biology of the parasite and/or clinical parameters in humans. Regarding the geographic dispersion of T. cruzi, evident differences were observed in the distribution of DTUs in distinct areas. For example, TcII is the main DTU detected in Brazilian patients from the central and southeastern regions, where there are also registers of TcVI as a secondary T. cruzi DTU. An important aspect observed in previous studies is that the genetic variability of T. cruzi can impact parasite infectivity, reproduction, and differentiation in the vectors. It has been proposed that T. cruzi DTU influences the host immune response and affects disease progression. Genetic aspects of the parasite play an important role in determining which host tissues will be infected, thus heavily influencing Chagas disease's pathogenesis. Several teams have investigated the correlation between T. cruzi DTU and the reactivation of Chagas disease. In agreement with these data, it is reasonable to suppose that the immunological condition of the patient, whether or not associated with the reactivation of the T. cruzi infection and the parasite strain, may have an important role in the pathogenesis of Chagas disease. In this context, understanding the genetics of T. cruzi and its biological and clinical implications will provide new knowledge that may contribute to additional strategies in the diagnosis and clinical outcome follow-up of patients with Chagas disease, in addition to the reactivation of immunocompromised patients infected with T. cruzi.

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.

Guide RNA Repertoires in the Main Lineages of Trypanosoma cruzi: High Diversity and Variable Redundancy Among Strains

Trypanosoma cruzi, as other kinetoplastids, has a complex mechanism of editing of mitochondrial mRNAs that requires guide RNAs (gRNAs) coded in DNA minicircles in the kinetoplast. There are many variations on this mechanism among species. mRNA editing and gRNA repertoires are almost unknown in T. cruzi. Here, gRNAs were inferred based on deep-sequenced minicircle hypervariable regions (mHVRs) and editing cascades were rebuilt in strains belonging to the six main T. cruzi lineages. Inferred gRNAs were clustered according to their sequence similarity to constitute gRNA classes. Extreme diversity of gRNA classes was observed, which implied highly divergent gRNA repertoires among different lineages, even within some lineages. In addition, a variable gRNA class redundancy (i.e., different gRNA classes editing the same mRNA region) was detected among strains. Some strains had upon four times more gRNA classes than others. Such variations in redundancy affected gRNA classes of all mRNAs in a concerted way, i.e., there are correlated variations in the number of gRNAs classes editing each mRNA. Interestingly, cascades were incomplete for components of the respiratory complex I in several strains. Finally, gRNA classes of different strains may potentially edit mitochondrial mRNAs from other lineages in the same way as they edit their own mitochondrial mRNAs, which is a prerequisite for biparental inheritance of minicircle in hybrids. We propose that genetic exchange and biparental inheritance of minicircles combined with minicircle drift due to (partial) random segregation of minicircles during kDNA replication is a suitable hypothesis to explain the divergences among strains and the high levels of gRNA redundancy in some strains. In addition, our results support that the complex I may not be required in some stages in the life cycle as previously shown and that linkage (in the same minicircle) of gRNAs that edit different mRNAs may prevent gRNA class lost in such stage.

PhyloQuant approach provides insights into Trypanosoma cruzi evolution using a systems-wide mass spectrometry-based quantitative protein profile

The etiological agent of Chagas disease, Trypanosoma cruzi, is a complex of seven genetic subdivisions termed discrete typing units (DTUs), TcI-TcVI and Tcbat. The relevance of T. cruzi genetic diversity to the variable clinical course of the disease, virulence, pathogenicity, drug resistance, transmission cycles and ecological distribution requires understanding the parasite origin and population structure. In this study, we introduce the PhyloQuant approach to infer the evolutionary relationships between organisms based on differential mass spectrometry-based quantitative features. In particular, large scale quantitative bottom-up proteomics features (MS1, iBAQ and LFQ) were analyzed using maximum parsimony, showing a correlation between T. cruzi DTUs and closely related trypanosomes' protein expression and sequence-based clustering. Character mapping enabled the identification of synapomorphies, herein the proteins and their respective expression profiles that differentiate T. cruzi DTUs and trypanosome species. The distance matrices based on phylogenetics and PhyloQuant clustering showed statistically significant correlation highlighting the complementarity between the two strategies. Moreover, PhyloQuant allows the identification of differentially regulated and strain/DTU/species-specific proteins, and has potential application in the identification of specific biomarkers and candidate therapeutic targets.

Differential gene expression of virulence factors modulates infectivity of TcI Trypanosoma cruzi strains

Trypanosoma cruzi is the etiological agent of Chagas disease, whose clinical outcome ranges from asymptomatic individuals to chronic fatal megasyndromes. Despite being central to pathogenesis, the regulation of parasite virulence factors' expression remains largely unknown. In this work, the relative expression of several parasite virulence factors between two TcI strains (Ninoa, low virulence and Qro, high virulence) was assessed by qRT-PCR of total and of polysome-associated mRNA, as well as by western blots. Trypomastigotes were also incubated with specific anti-sense morpholino oligonucleotides to block the translation of a selected virulence factor, calreticulin, in both strains. Ninoa trypomastigotes showed significantly lower levels of trypomastigote-decay acceleration factor, complement regulatory protein, complement C2 receptor inhibitor trispanning, and glycoproteins 82 and 90 mRNAs compared with Qro. There was a significantly lower recruitment of complement regulatory protein and complement C2 receptor inhibitor trispanning mRNAs to polysomes and higher recruitment of MASP mRNA to monosomes in Ninoa strain. Calreticulin mRNA displayed both a higher total mRNA level and recruitment to translationally active polysomes in the Ninoa strain (low virulence) than in the Qro strain (high virulence). When calreticulin was downregulated by ≈ 50% by anti-sense morpholino oligonucleotides, a significant decrease of parasite invasion in mammalian cells was found in both strains. Calreticulin downregulation, however, only increased significantly the activation of the complement system by Ninoa trypomastigotes. These results suggest a role for the regulation of virulence factors' gene expression in the differential virulence among T. cruzi strains. Furthermore, a possible function of calreticulin in parasite invasion not related to its binding to complement factors is shown.

Multilocus sequence typing analysis of Leishmania clinical isolates from cutaneous leishmaniasis patients of Iran

Cutaneous leishmaniasis (CL) is mainly caused by L. major and L. tropica in Old World and might be represented as typical skin lesion(s) or sometimes as a spectrum of atypical manifestations. We applied multilocus sequence typing (MLST) to explore genetic variations of Leishmania strains isolated from atypical vs. typical CL patients from Iran. A PCR-sequencing was performed for seven housekeeping genes (g6pd, mpi, asat, icd, 6pgd, fh, and trys) and genetic diversity indices and phylogenetic relationships were analyzed. A total of 41 isolates of L. major (28/41) and L. tropica (13/41) from 21 (51.2%) atypical CL and 20 (48.8%) typical CL cases were included. A set of additional sequences of 41 strains of 17 species of Leishmania were retrieved from databases. Different SNP variations were detected and the highest rate of heterozygous sites was found in g6pd and 6pgd genes (6 sites) for L. tropica and in asat and 6pgd genes (7 sites) for L. major strains. All strains were clustered into 58 unique sequence types (STs) including 17 STs related to 41 strains of Leishmania of this study. Concatenated tree clustered all strains in 6 main clades (A to F) including L. major (clade D) and L. tropica (clade B) strains. Two strains of L. major (codes 28 and 42) with highest nucleotide variations were more close to L. tropica and were grouped in Clade B. All of the STs were related in clonal complexes by using eBURST with the prediction of founder genotypes. A high rate of genetic variations and heterozygocity was evident in L. tropica and L. major strains; nevertheless, there was no significant difference in the diversity of Leishmania strains between typical CL and atypical CL groups. This study represents the first successful application of MLST approach to L. tropica and L. major strains in Iran.

Development of a Multilocus Sequence Typing Scheme for Giardia intestinalis

Giardia intestinalis is an intestinal protozoan most commonly found in humans. It has been grouped into 8 assemblages (A-H). Markers such as the glutamate dehydrogenase gene, triose phosphate isomerase and beta-giardin (β-giardin) have been widely used for genotyping. In addition, different genetic targets have been proposed as a valuable alternative to assess diversity and genetics of this microorganism. Thus, our objective was to evaluate new markers for the study of the diversity and intra-taxa genetic structure of G. intestinalis in silico and in DNA obtained from stool samples. We analysed nine constitutive genes in 80 complete genome sequences and in a group of 24 stool samples from Colombia. Allelic diversity was evaluated by locus and for the concatenated sequence of nine loci that could discriminate up to 53 alleles. Phylogenetic reconstructions allowed us to identify AI, AII and B assemblages. We found evidence of intra- and inter-assemblage recombination events. Population structure analysis showed genetic differentiation among the assemblages analysed.

Evidence of hybridization, mitochondrial introgression and biparental inheritance of the kDNA minicircles in Trypanosoma cruzi I

Background

Genetic exchange in Trypanosoma cruzi is controversial not only in relation to its frequency, but also to its mechanism. Parasexual genetic exchange has been proposed based on laboratory hybrids, but population genomics strongly suggests meiosis in T. cruzi. In addition, mitochondrial introgression has been reported several times in natural isolates although its mechanism is not fully understood yet. Moreover, hybrid T. cruzi DTUs (TcV and TcVI) have inherited at least part of the kinetoplastic DNA (kDNA = mitochondrial DNA) from both parents.

Methodology/Principal findings

In order to address such topics, we sequenced and analyzed fourteen nuclear DNA fragments and three kDNA maxicircle genes in three TcI stocks which are natural clones potentially involved in events of genetic exchange. We also deep-sequenced (a total of 6,146,686 paired-end reads) the minicircle hypervariable region (mHVR) of the kDNA in such three strains. In addition, we analyzed the DNA content by flow cytometry to address cell ploidy. We observed that most polymorphic sites in nuclear loci showed a hybrid pattern in one cloned strain and the other two cloned strains were compatible as parental strains (or nearly related to the true parents). The three clones had almost the same ploidy and the DNA content was similar to the reference strain Sylvio (a nearly diploid strain). Despite maxicircle genes evolve faster than nuclear housekeeping ones, we detected no polymorphisms in the sequence of three maxicircle genes showing mito-nuclear discordance. Lastly, the hybrid stock shared 66% of its mHVR clusters with one putative parent and 47% with the other one; in contrast, the putative parental stocks shared less than 30% of the mHVR clusters between them.

Conclusions/significance

The results suggest a reductive division, a natural hybridization, biparental inheritance of the minicircles in the hybrid and maxicircle introgression. The models including such phenomena and explaining the relationships between these three clones are discussed.

Chagas disease, an important public health problem in Latin America, is caused by the parasite Trypanosoma cruzi. Despite being a widely studied parasite, several questions on the biology of genetic exchange remain unanswered. Population genomic studies have inferred meiosis in T. cruzi, but this cellular division mechanism has not been observed in laboratory yet. In addition, previous results suggest that mitochondrial DNA (called kDNA) may be inherited from both parents in hybrids. Here, we analyzed a hybrid strain and its potential parents to address the mechanisms of genetic exchange at nuclear and mitochondrial levels. We observed that the hybrid strain had heterozygous patterns and DNA content compatible with a meiosis event. Also, we observed that the evolutionary histories of nuclear DNA and kDNA maxicircles were discordant and that the three strains shared identical DNA sequences. Mitochondrial introgression of maxicircle DNA from one genotype to another may explain this observation. In addition, we demonstrated that the hybrid strain shared kDNA minicircles with both parental strains. Our results suggest that hybridization implied meiosis and biparental inheritance of the kDNA. Further research is required to address such phenomena in detail.

Development of a Multilocus sequence typing (MLST) scheme for Pan-Leishmania

Since the description of the Leishmania genus, its identification and organization have been a challenge. A high number of molecular markers have been developed to resolve phylogenetic differences at the species level and for addressing key epidemiological and population genetics questions. Based on Multilocus enzyme electrophoresis (MLEE), Multilocus sequence typing (MLST) schemes have been developed using different gene candidates. From 38 original gene targets proposed by other authors, 27 of them were chosen. In silico selection was made by analyzing free access genomic sequence data of 33 Leishmania species, one Paraleishmania representative, and one outgroup, in order to select the best 15 loci. De novo amplifications and primers redesign of these 15 genes were analyzed over a panel of 20 reference strains and isolates. Phylogenetic analysis was made at every step. Two MLST schemes were selected. The first one was based on the analysis of three-gene fragments, and it is suitable for species assignment as well as basic phylogenetic studies. By the addition of seven-genes, an approach based on the analysis of ten-gene fragments was also proposed. This is the first work that two optimized MLST schemes have been suggested, validated against a phylogenetically diverse panel of Leishmania isolates. MLST is potentially a powerful phylogenetic approach, and most probably the new gold standard for Leishmania spp. characterization.

Chagas Disease in the United States: a Public Health Approach

Trypanosoma cruzi is the etiological agent of Chagas disease, usually transmitted by triatomine vectors. An estimated 20 to 30% of infected individuals develop potentially lethal cardiac or gastrointestinal disease. Sylvatic transmission cycles exist in the southern United States, involving 11 triatomine vector species and infected mammals such as rodents, opossums, and dogs. Nevertheless, imported chronic T. cruzi infections in migrants from Latin America vastly outnumber locally acquired human cases. Benznidazole is now FDA approved, and clinical and public health efforts are under way by researchers and health departments in a number of states. Making progress will require efforts to improve awareness among providers and patients, data on diagnostic test performance and expanded availability of confirmatory testing, and evidence-based strategies to improve access to appropriate management of Chagas disease in the United States.

Multilocus Sequence Typing (MLST) and Random Polymorphic DNA (RAPD) Comparisons of Geographic Isolates of Neoparamoeba perurans, the Causative Agent of Amoebic Gill Disease

Neoparamoba perurans, is the aetiological agent of amoebic gill disease (AGD), a disease that affects farmed Atlantic salmon worldwide. Multilocus sequence typing (MLST) and Random Amplified Polymorphic DNA (RAPD) are PCR-based typing methods that allow for the highly reproducible genetic analysis of population structure within microbial species. To the best of our knowledge, this study represents the first use of these typing methods applied to N. perurans with the objective of distinguishing geographical isolates. These analyses were applied to a total of 16 isolates from Australia, Canada, Ireland, Scotland, Norway, and the USA. All the samples from Australia came from farm sites on the island state of Tasmania. Genetic polymorphism among isolates was more evident from the RAPD analysis compared to the MLST that used conserved housekeeping genes. Both techniques consistently identified that isolates of N. perurans from Tasmania, Australia were more similar to each other than to the isolates from other countries. While genetic differences were identified between geographical isolates, a BURST analysis provided no evidence of a founder genotype. This suggests that emerging outbreaks of AGD are not due to rapid translocation of this important salmonid pathogen from the same area.

Elucidating diversity in the class composition of the minicircle hypervariable region of Trypanosoma cruzi: New perspectives on typing and kDNA inheritance

Background

Trypanosoma cruzi, the protozoan causative of Chagas disease, is classified into six main Discrete Typing Units (DTUs): TcI-TcVI. This parasite has around 10⁵ copies of the minicircle hypervariable region (mHVR) in their kinetoplastic DNA (kDNA). The genetic diversity of the mHVR is virtually unknown. However, cross-hybridization assays using mHVRs showed hybridization only between isolates belonging to the same genetic group. Nowadays there is no methodologic approach with a good sensibility, specificity and reproducibility for direct typing on biological samples. Due to its high copy number and apparently high diversity, mHVR becomes a good target for typing.

Methodology/Principal findings

Around 22 million reads, obtained by amplicon sequencing of the mHVR, were analyzed for nine strains belonging to six T. cruzi DTUs. The number and diversity of mHVR clusters was variable among DTUs and even within a DTU. However, strains of the same DTU shared more mHVR clusters than strains of different DTUs and clustered together. In addition, hybrid DTUs (TcV and TcVI) shared similar percentages (1.9–3.4%) of mHVR clusters with their parentals (TcII and TcIII). Conversely, just 0.2% of clusters were shared between TcII and TcIII suggesting biparental inheritance of the kDNA in hybrids. Sequencing at low depth (20,000–40,000 reads) also revealed 95% of the mHVR clusters for each of the analyzed strains. Finally, the method revealed good correlation in cluster identity and abundance between different replications of the experiment (r = 0.999).

Conclusions/Significance

Our work sheds light on the sequence diversity of mHVRs at intra and inter-DTU level. The mHVR amplicon sequencing workflow described here is a reproducible technique, that allows multiplexed analysis of hundreds of strains and results promissory for direct typing on biological samples in a future. In addition, such approach may help to gain knowledge on the mechanisms of the minicircle evolution and phylogenetic relationships among strains.

Chagas disease is an important public health problem in Latin America showing a wide diversity of clinical manifestations and epidemiological patterns. It is caused by the parasite Trypanosoma cruzi. This parasite is genetically diverse and classified into six main lineages. However, the relationship between intra-specific genetic diversity and clinical or epidemiological features is not clear, mainly because low sensitivity for direct typing on biological samples. For this reason, genetic markers with high copy number are required to achieve sensitivity. Here, we deep sequenced and analyzed a DNA region present in the large mitochondria of the parasite (named as mHVR, 10⁵ copies per parasite) from strains belonging to the six main lineages in order to analyze mHVR diversity and to evaluate its usefulness for typing. Despite the high sequence diversity, strains of the same lineage shared more sequences than strains of different lineages. Curiously, hybrid lineages shared mHVR sequences with both parents suggesting that mHVR (and DNA minicircles from the mitochondria) are inherited from both parentals. The mHVR amplicon sequencing workflow proposed here is reproducible and, potentially, it would be useful for typing hundreds of biological samples at time. It also provides a valuable approach to perform evolutionary and functional studies.

Evaluation of the multispecies coalescent method to explore intra-Trypanosoma cruzi I relationships and genetic diversity

Chagas Disease is a zoonosis caused by the parasite Trypanosoma cruzi. Several high-resolution markers have subdivided T. cruzi taxon into at least seven lineages or Discrete Typing Units (DTUs) (TcI-TcVI and TcBat). Trypanosoma cruzi I is the most diverse and geographically widespread DTU. Recently a TcI genotype related to domestic cycles was proposed and named as TcIDOM. Herein, we combined traditional markers and housekeeping genes and applied a Multispecies Coalescent method to explore intra-TcI relationships, lineage boundaries and genetic diversity in a random set of isolates and DNA sequences retrieved from Genbank from different countries in the Americas. We found further evidence supporting TcIDOM as an independent and emerging genotype of TcI at least in Colombia and Venezuela. We also found evidence of high phylogenetic incongruence between parasite's gene trees (including introgression) and embedded species trees, and a lack of genetic structure among geography and hosts, illustrating the complex dynamics and epidemiology of TcI across the Americas. These findings provide novel insights into T. cruzi systematics and epidemiology and support the need to assess parasite diversity and lineage boundaries through hypothesis testing using different approaches to those traditionally employed, including the Bayesian Multispecies coalescent method.

Dissecting the phyloepidemiology of Trypanosoma cruzi I (TcI) in Brazil by the use of high resolution genetic markers

BACKGROUND

Trypanosoma cruzi, the causal agent of Chagas disease, is monophyletic but genetically heterogeneous. It is currently represented by six genetic lineages (Discrete Typing Units, DTUs) designated TcI-TcVI. TcI is the most geographically widespread and genetically heterogeneous lineage, this as is evidenced by a wide range of genetic markers applied to isolates spanning a vast geographic range in Latin America.

METHODOLOGY/PRINCIPAL FINDINGS

In total, 78 TcI isolated from hosts and vectors distributed in 5 different biomes of Brazil, were analyzed using 6 nuclear housekeeping genes, 25 microsatellite loci and one mitochondrial marker. Nuclear markers reveal substantial genetic diversity, significant gene flow between biomes, incongruence in phylogenies, and haplotypic analysis indicative of intra-DTU genetic exchange. Phylogenetic reconstructions based on mitochondrial and nuclear loci were incongruent, and consistent with introgression. Structure analysis of microsatellite data reveals that, amongst biomes, the Amazon is the most genetically diverse and experiences the lowest level of gene flow. Investigation of population structure based on the host species/genus, indicated that Didelphis marsupialis might play a role as the main disperser of TcI.

CONCLUSIONS/SIGNIFICANCE

The present work considers a large TcI sample from different hosts and vectors spanning multiple ecologically diverse biomes in Brazil. Importantly, we combine fast and slow evolving markers to contribute to the epizootiological understanding of TcI in five distinct Brazilian biomes. This constitutes the first instance in which MLST analysis was combined with the use of MLMT and maxicircle markers to evaluate the genetic diversity of TcI isolates in Brazil. Our results demonstrate the existence of substantial genetic diversity and the occurrence of introgression events. We provide evidence of genetic exchange in TcI isolates from Brazil and of the relative isolation of TcI in the Amazon biome. We observe the absence of strict associations with TcI genotypes to geographic areas and/or host species.

Multilocus sequence typing: genetic diversity in Trypanosoma cruzi I (TcI) isolates from Brazilian didelphids

BACKGROUND

Trypanosoma cruzi is a protozoan parasite characterized by extensive genetic heterogeneity. There are currently six recognised, genetically distinct, monophyletic clades designated discrete typing units (DTUs). TcI has the broadest geographical range and most genetic diversity evidenced by a wide range of genetic markers applied to isolates spanning a vast geographical range across Latin America. However, little is known of the diversity of TcI that exists within sylvatic mammals across the geographical expanse of Brazil.

RESULTS

Twenty-nine sylvatic TcI isolates spanning multiple ecologically diverse biomes across Brazil were analyzed by the application of multilocus sequence typing (MLST) using four nuclear housekeeping genes. Results revealed extensive genetic diversity and also incongruence among individual gene trees. There was no association of intralineage genotype with geography or with any particular biome, with the exception of isolates from Caatinga that formed a single cluster. However, haplotypic analyses of METIII and LYT1 constitutive markers provided evidence of recombination events in two isolates derived from Didelphis marsupialis and D. albiventris, respectively. For diversity studies all possible combinations of markers were assessed with the objective of selecting the combination of gene targets that are most resolutive using the minimum number of genes. A panel of just three gene fragments (DHFR-TS, LYT1 and METIII) discriminated 26 out of 35 genotypes.

CONCLUSIONS

These findings showed geographical association of genotypes clustering in Caatinga but more characteristically TcI genotypes widely distributed without specific association to geographical areas or biomes. Importantly, we detected the signature of recombination events at the nuclear level evidenced by haplotypic analysis and incongruence.

Introgression of the Kinetoplast DNA: An Unusual Evolutionary Journey in Trypanosoma cruzi

INTRODUCTION

Phylogenetic relationships between different lineages of Trypanosoma cruzi, the agent of Chagas disease, have been controversial for several years. However, recent phylogenetic and phylogenomic analyses clarified the nuclear relationships among such lineages. However, incongruence between nuclear and kinetoplast DNA phylogenies has emerged as a new challenge. This incongruence implies several events of mitochondrial introgression at evolutionary level. However, the mechanism that gave origin to introgressed lineages is unknown. Here, I will review and discuss how maxicircles of the kinetoplast were horizontally and vertically transferred between different lineages of T. cruzi.

CONCLUSION

Finally, I will discuss what we know - and what we don't - about the kDNA transference and inheritance in the context of sexual reproduction in this parasite.

Cytochrome c oxidase subunit 1 gene as a DNA barcode for discriminating Trypanosoma cruzi DTUs and closely related species

BACKGROUND

The DNA barcoding system using the cytochrome c oxidase subunit 1 mitochondrial gene (cox1 or COI) is highly efficient for discriminating vertebrate and invertebrate species. In the present study, we examined the suitability of cox1 as a marker for Trypanosoma cruzi identification from other closely related species. Additionally, we combined the sequences of cox1 and the nuclear gene glucose-6-phosphate isomerase (GPI) to evaluate the occurrence of mitochondrial introgression and the presence of hybrid genotypes.

METHODS

Sixty-two isolates of Trypanosoma spp. obtained from five of the six Brazilian biomes (Amazon Forest, Atlantic Forest, Caatinga, Cerrado and Pantanal) were sequenced for cox1 and GPI gene fragments. Phylogenetic trees were reconstructed using neighbor-joining, maximum likelihood, parsimony and Bayesian inference methods. Molecular species delimitation was evaluated through pairwise intraspecific and interspecific distances, Automatic Barcode Gap Discovery, single-rate Poisson Tree Processes and multi-rate Poisson Tree Processes.

RESULTS

Both cox1 and GPI genes recognized and differentiated T. cruzi, Trypanosoma cruzi marinkellei, Trypanosoma dionisii and Trypanosoma rangeli. Cox1 discriminated Tcbat, TcI, TcII, TcIII and TcIV. Additionally, TcV and TcVI were identified as a single group. Cox1 also demonstrated diversity in the discrete typing units (DTUs) TcI, TcII and TcIII and in T. c. marinkellei and T. rangeli. Cox1 and GPI demonstrated TcI and TcII as the most genetically distant branches, and the position of the other T. cruzi DTUs differed according to the molecular marker. The tree reconstructed with concatenated cox1 and GPI sequences confirmed the separation of the subgenus Trypanosoma (Schizotrypanum) sp. and the T. cruzi DTUs TcI, TcII, TcIII and TcIV. The evaluation of single nucleotide polymorphisms (SNPs) was informative for DTU differentiation using both genes. In the cox1 analysis, one SNP differentiated heterozygous hybrids from TcIV sequences. In the GPI analysis one SNP discriminated Tcbat from TcI, while another SNP distinguished TcI from TcIII.

CONCLUSIONS

DNA barcoding using the cox1 gene is a reliable tool to distinguish T. cruzi from T. c. marinkellei, T. dionisii and T. rangeli and identify the main T. cruzi genotypes.

Real-time PCR strategy for the identification of Trypanosoma cruzi discrete typing units directly in chronically infected human blood

The protozoan Trypanosoma cruzi is the causative agent of Chagas disease, a major public health problem in Latin America. This parasite has a complex population structure comprised by six or seven major evolutionary lineages (discrete typing units or DTUs) TcI-TcVI and TcBat, some of which have apparently resulted from ancient hybridization events. Because of the existence of significant biological differences between these lineages, strain characterization methods have been essential to study T. cruzi in its different vectors and hosts. However, available methods can be laborious and costly, limited in resolution or sensitivity. In this study, a new genotyping strategy by real-time PCR to identify each of the six DTUs in clinical blood samples have been developed and evaluated. Two nuclear (SL-IR and 18S rDNA) and two mitochondrial genes (COII and ND1) were selected to develop original primers. The method was evaluated with eight genomic DNA of T. cruzi populations belonging to the six DTUs, one genomic DNA of Trypanosoma rangeli, and 53 blood samples from individuals with chronic Chagas disease. The assays had an analytical sensitivity of 1-25fg of DNA per reaction tube depending on the DTU analyzed. The selectivity of trials with 20fg/μL of genomic DNA identified each DTU, excluding non-targets DTUs in every test. The method was able to characterize 67.9% of the chronically infected clinical samples with high detection of TcII followed by TcI. With the proposed original genotyping methodology, each DTU was established with high sensitivity after a single real-time PCR assay. This novel protocol reduces carryover contamination, enables detection of each DTU independently and in the future, the quantification of each DTU in clinical blood samples.

Identification of Trypanosoma cruzi Discrete Typing Units (DTUs) in Latin-American migrants in Barcelona (Spain)

Trypanosoma cruzi, the causative agent of Chagas disease, is divided into six Discrete Typing Units (DTUs): TcI-TcVI. We aimed to identify T. cruzi DTUs in Latin-American migrants in the Barcelona area (Spain) and to assess different molecular typing approaches for the characterization of T. cruzi genotypes. Seventy-five peripheral blood samples were analyzed by two real-time PCR methods (qPCR) based on satellite DNA (SatDNA) and kinetoplastid DNA (kDNA). The 20 samples testing positive in both methods, all belonging to Bolivian individuals, were submitted to DTU characterization using two PCR-based flowcharts: multiplex qPCR using TaqMan probes (MTq-PCR), and conventional PCR. These samples were also studied by sequencing the SatDNA and classified as type I (TcI/III), type II (TcII/IV) and type I/II hybrid (TcV/VI). Ten out of the 20 samples gave positive results in the flowcharts: TcV (5 samples), TcII/V/VI (3) and mixed infections by TcV plus TcII (1) and TcV plus TcII/VI (1). By SatDNA sequencing, we classified the 20 samples, 19 as type I/II and one as type I. The most frequent DTU identified by both flowcharts, and suggested by SatDNA sequencing in the remaining samples with low parasitic loads, TcV, is common in Bolivia and predominant in peripheral blood. The mixed infection by TcV-TcII was detected for the first time simultaneously in Bolivian migrants. PCR-based flowcharts are very useful to characterize DTUs during acute infection. SatDNA sequence analysis cannot discriminate T. cruzi populations at the level of a single DTU but it enabled us to increase the number of characterized cases in chronically infected patients.

Is Predominant Clonal Evolution a Common Evolutionary Adaptation to Parasitism in Pathogenic Parasitic Protozoa, Fungi, Bacteria, and Viruses?

We propose that predominant clonal evolution (PCE) in microbial pathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure. The main features of PCE are (1) strong linkage disequilibrium, (2) the widespread occurrence of stable genetic clusters blurred by occasional bouts of genetic exchange ('near-clades'), (3) the existence of a "clonality threshold", beyond which recombination is efficiently countered by PCE, and near-clades irreversibly diverge. We hypothesize that the PCE features are not mainly due to natural selection but also chiefly originate from in-built genetic properties of pathogens. We show that the PCE model obtains even in microbes that have been considered as 'highly recombining', such as Neisseria meningitidis, and that some clonality features are observed even in Plasmodium, which has been long described as panmictic. Lastly, we provide evidence that PCE features are also observed in viruses, taking into account their extremely fast genetic turnover. The PCE model provides a convenient population genetic framework for any kind of micropathogen. It makes it possible to describe convenient units of analysis (clones and near-clades) for all applied studies. Due to PCE features, these units of analysis are stable in space and time, and clearly delimited. The PCE model opens up the possibility of revisiting the problem of species definition in these organisms. We hypothesize that PCE constitutes a major evolutionary strategy for protozoa, fungi, bacteria, and viruses to adapt to parasitism.

Molecular Diversity of Trypanosoma cruzi Detected in the Vector Triatoma protracta from California, USA

BACKGROUND

Trypanosoma cruzi, causative agent of Chagas disease in humans and dogs, is a vector-borne zoonotic protozoan parasite that can cause fatal cardiac disease. While recognized as the most economically important parasitic infection in Latin America, the incidence of Chagas disease in the United States of America (US) may be underreported and even increasing. The extensive genetic diversity of T. cruzi in Latin America is well-documented and likely influences disease progression, severity and treatment efficacy; however, little is known regarding T. cruzi strains endemic to the US. It is therefore important to expand our knowledge on US T. cruzi strains, to improve upon the recognition of and response to locally acquired infections.

METHODOLOGY/PRINCIPLE FINDINGS

We conducted a study of T. cruzi molecular diversity in California, augmenting sparse genetic data from southern California and for the first time investigating genetic sequences from northern California. The vector Triatoma protracta was collected from southern (Escondido and Los Angeles) and northern (Vallecito) California regions. Samples were initially screened via sensitive nuclear repetitive DNA and kinetoplast minicircle DNA PCR assays, yielding an overall prevalence of approximately 28% and 55% for southern and northern California regions, respectively. Positive samples were further processed to identify discrete typing units (DTUs), revealing both TcI and TcIV lineages in southern California, but only TcI in northern California. Phylogenetic analyses (targeting COII-ND1, TR and RB19 genes) were performed on a subset of positive samples to compare Californian T. cruzi samples to strains from other US regions and Latin America. Results indicated that within the TcI DTU, California sequences were similar to those from the southeastern US, as well as to several isolates from Latin America responsible for causing Chagas disease in humans.

CONCLUSIONS/SIGNIFICANCE

Triatoma protracta populations in California are frequently infected with T. cruzi. Our data extend the northern limits of the range of TcI and identify a novel genetic exchange event between TcI and TcIV. High similarity between sequences from California and specific Latin American strains indicates US strains may be equally capable of causing human disease. Additional genetic characterization of Californian and other US T. cruzi strains is recommended.

The diversity and expansion of the trans-sialidase gene family is a common feature in Trypanosoma cruzi clade members

Trans-sialidase (TS) is a polymorphic protein superfamily described in members of the protozoan genus Trypanosoma. Of the eight TS groups recently described, TS group I proteins (some of which have catalytic activity) are present in the distantly related Trypanosoma brucei and Trypanosoma cruzi phylogenetic clades, whereas other TS groups have only been described in some species belonging to the T. cruzi clade. In the present study we analyzed the repertoire, distribution and phylogenetic relationships of TS genes among species of the T. cruzi clade based on sequence similarity, multiple sequence alignment and tree-reconstruction approaches using TS sequences obtained with the aid of PCR-based strategies or retrieved from genome databases. We included the following representative isolates of the T. cruzi clade from South America: T. cruzi, T. cruzi Tcbat, Trypanosoma cruzi marinkellei, Trypanosoma dionisii, Trypanosoma rangeli and Trypanosoma conorhini. The cloned sequences encoded conserved TS protein motifs Asp-box and VTVxNVxLYNR but lacked the FRIP motif (conserved in TS group I). The T. conorhini sequences were the most divergent. The hybridization patterns of TS probes with chromosomal bands confirmed the abundance of these sequences in species in the T. cruzi clade. Divergence and relationship analysis placed most of the TS sequences in the groups defined in T. cruzi. Further examination of members of TS group II, which includes T. cruzi surface glycoproteins implicated in host cell attachment and invasion, showed that sequences of T. cruzi Tcbat grouped with those of T. cruzi genotype TcI. Our analysis indicates that different members of the T. cruzi clade, with different vertebrate hosts, vectors and pathogenicity, share the extensive expansion and sequence diversification of the TS gene family. Altogether, our results are congruent with the evolutionary history of the T. cruzi clade and represent a contribution to the understanding of the molecular evolution and role of TS proteins in trypanosomes.

Genetic diversity of Trypanosoma cruzi in bats, and multilocus phylogenetic and phylogeographical analyses supporting Tcbat as an independent DTU (discrete typing unit)

Trypanosoma cruzi is a complex of phenotypically and genetically diverse isolates distributed in six discrete typing units (DTUs) designated as TcI-TcVI. Five years ago, T. cruzi isolates from Brazilian bats showing unique patterns of traditional ribosomal and spliced leader PCRs not clustering into any of the six DTUs were designated as the Tcbat genotype. In the present study, phylogenies inferred using SSU rRNA (small subunit of ribosomal rRNA), gGAPDH (glycosomal glyceraldehyde 3-phosphate dehydrogenase) and Cytb (cytochrome b) genes strongly supported Tcbat as a monophyletic lineage prevalent in Brazil, Panama and Colombia. Providing strong support for Tcbat, sequences from 37 of 47 nuclear and 12 mitochondrial genes (retrieved from a draft genome of Tcbat) and reference strains of all DTUs available in databanks corroborated Tcbat as an independent DTU. Consistent with previous studies, multilocus analysis of most nuclear genes corroborated the evolution of T. cruzi from bat trypanosomes its divergence into two main phylogenetic lineages: the basal TcII; and the lineage clustering TcIV, the clade comprising TcIII and the sister groups TcI-Tcbat. Most likely, the common ancestor of Tcbat and TcI was a bat trypanosome. However, the results of the present analysis did not support Tcbat as the ancestor of all DTUs. Despite the insights provided by reports of TcIII, TcIV and TcII in bats, including Amazonian bats harbouring TcII, further studies are necessary to understand the roles played by bats in the diversification of all DTUs. We also demonstrated that in addition to value as molecular markers for DTU assignment, Cytb, ITS rDNA and the spliced leader (SL) polymorphic sequences suggest spatially structured populations of Tcbat. Phylogenetic and phylogeographical analyses, multiple molecular markers specific to Tcbat, and the degrees of sequence divergence between Tcbat and the accepted DTUs strongly support the definitive classification of Tcbat as a new DTU.

Between a bug and a hard place: Trypanosoma cruzi genetic diversity and the clinical outcomes of Chagas disease

Over the last 30 years, concomitant with successful transnational disease control programs across Latin America, Chagas disease has expanded from a neglected, endemic parasitic infection of the rural poor to an urbanized chronic disease, and now a potentially emergent global health problem. Trypanosoma cruzi infection has a highly variable clinical course, ranging from complete absence of symptoms to severe and often fatal cardiovascular and/or gastrointestinal manifestations. To date, few correlates of clinical disease progression have been identified. Elucidating a putative role for T. cruzi strain diversity in Chagas disease pathogenesis is complicated by the scarcity of parasites in clinical specimens and the limitations of our contemporary genotyping techniques. This article systematically reviews the historical literature, given our current understanding of parasite genetic diversity, to evaluate the evidence for any association between T. cruzi genotype and chronic clinical outcome, risk of congenital transmission or reactivation and orally transmitted outbreaks.

Evolution of Trypanosoma cruzi: clarifying hybridisations, mitochondrial introgressions and phylogenetic relationships between major lineages

Several different models of Trypanosoma cruzi evolution have been proposed. These models suggest that scarce events of genetic exchange occurred during the evolutionary history of this parasite. In addition, the debate has focused on the existence of one or two hybridisation events during the evolution of T. cruzi lineages. Here, we reviewed the literature and analysed available sequence data to clarify the phylogenetic relationships among these different lineages. We observed that TcI, TcIII and TcIV form a monophyletic group and that TcIII and TcIV are not, as previously suggested, TcI-TcII hybrids. Particularly, TcI and TcIII are sister groups that diverged around the same time that a widely distributed TcIV split into two clades (TcIV_S and TcIV_N). In addition, we collected evidence that TcIII received TcIV_S kDNA by introgression on several occasions. Different demographic hypotheses (surfing and asymmetrical introgression) may explain the origin and expansion of the TcIII group. Considering these hypotheses, genetic exchange should have been relatively frequent between TcIII and TcIV_S in the geographic area in which their distributions overlapped. In addition, our results support the hypothesis that two independent hybridisation events gave rise to TcV and TcVI. Consequently, TcIV_S kDNA was first transferred to TcIII and later to TcV and TcVI in TcII/TcIII hybridisation events.

Molecular genotyping of Trypanosoma cruzi for lineage assignment and population genetics

Trypanosoma cruzi, the etiological agent of Chagas disease, remains a major public health problem in Latin America. Infection with T. cruzi is lifelong and can lead to a spectrum of pathological sequelae ranging from subclinical to lethal cardiac and/or gastrointestinal complications. Isolates of T. cruzi can be assigned to six genetic lineages or discrete typing units (DTUs), which are broadly associated with disparate ecologies, transmission cycles, and geographical distributions. This extensive genetic diversity is also believed to contribute to the clinical variation observed among chagasic patients. Unravelling the population structure of T. cruzi is fundamental to understanding Chagas disease epidemiology, developing control strategies, and resolving the relationship between parasite genotype and clinical prognosis. To date, no single, widely validated, genetic target allows unequivocal resolution to DTU-level. In this chapter we present standardized methods for strain DTU assignment using PCR-restriction fragment length polymorphism analysis (PCR-RFLP) and nuclear multilocus sequence typing (MLST). PCR-RFLPs have the advantages of simplicity and reproducibility, requiring limited expertise and few laboratory consumables. MLST data are more laborious to generate but more informative; DNA sequences are readily transferable between research groups and amenable to recombination detection and intra-lineage analyses. We also recommend a mitochondrial (maxicircle) MLST scheme and a panel of 28 microsatellite loci for higher resolution population genetics studies. Due to the scarcity of T. cruzi in blood and tissue, all of these genotyping techniques have limited sensitivity when applied directly to clinical or biological specimens, particularly when targets are single (MLST) or low copy number (PCR-RFLPs). We therefore describe essential protocols to isolate parasites, derive biological clones, and extract T. cruzi genomic DNA from field and clinical samples.