Preponderant clonal evolution of Trypanosoma cruzi I from Argentinean Chaco revealed by Multilocus Sequence Typing (MLST)

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.

Microevolution and subspecific taxonomy of Trypanosoma cruzi

Trypanosoma cruzi, the agent of Chagas disease, is a highly polymorphic species, subdivided into 6 main evolutionary lineages or near-clades (formerly discrete typing units or DTUs). An additional near-clade (TC-bat) has recently been evidenced. This pattern is considered to be the result of predominant clonal evolution (PCE). PCE is compatible with occasional mating/hybridization, which do not break the prevalent pattern of clonal evolution, the main trait of it being the presence of Multigene Bifurcating Trees (MGBTs) at all evolutionary levels ("clonal frame"). The development of highly resolutive genetic (microsatellites*) and genomic (sequencing and multi-single nucleotide polymorphism {SNP}* typing) markers shows that PCE also operates at a microevolutionary* level within each of the near-clades ("Russian doll pattern"), in spite of occasional meiosis and hybridization events. Within each near-clade, one can evidence widespread clonal multilocus genotypes*, linkage disequilibrium*, Multigene Bifurcating Trees and lesser near-clades. The within near-clade population structure is like a miniature picture of that of the whole species, suggesting gradual rather than saltatory evolution. Additional data are required to evaluate the stability of these lesser near-clades in the long run and to evaluate the need for an adequate nomenclature for this microevolutionary level.

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.

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.

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.

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.

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.

Evidence and importance of genetic exchange among field populations of Trypanosoma cruzi

Many eukaryotic pathogenic microorganisms that were previously assumed to propagate clonally have retained cryptic sexual cycles. The principal reproductive mode of Trypanosoma cruzi, the aetiological agent of Chagas disease, remains a controversial topic. Despite the existence of two recent natural hybrid lineages, a pervasive view is that recombination has been restrained at an evolutionary scale and is of little epidemiological relevance to contemporary parasite populations. This article reviews the growing number of field studies which indicate that natural hybridization in T. cruzi may be frequent, non-obligatory and idiosyncratic; potentially involving independent exchange of kinetoplast and nuclear genetic material as well as canonical meiotic mechanisms. Together these observations now challenge the traditional paradigm of preponderate clonal evolution in T. cruzi and highlight the need for additional, intensive and appropriately sampled field surveys, complemented by high resolution, combined nuclear and mitochondrial population genetics analyses.

The population genetics of Trypanosoma cruzi revisited in the light of the predominant clonal evolution model

Comparing the population structure of Trypanosoma cruzi with that of other pathogens, including parasitic protozoa, fungi, bacteria and viruses, shows that the agent of Chagas disease shares typical traits with many other species, related to a predominant clonal evolution (PCE) pattern: statistically significant linkage disequilibrium, overrepresented multilocus genotypes, near-clades (genetic subdivisions somewhat blurred by occasional genetic exchange/hybridization) and "Russian doll" patterns (PCE is observed, not only at the level of the whole species, but also, within the near-clades). Moreover, T. cruzi population structure exhibits linkage with the diversity of several strongly selected genes, with gene expression profiles, and with some major phenotypic traits. We discuss the evolutionary significance of these results, and their implications in terms of applied research (molecular epidemiology/strain typing, analysis of genes of interest, vaccine and drug design, immunological diagnosis) and of experimental evolution. Lastly, we revisit the long-term debate of describing new species within the T. cruzi taxon.