Microevolution of Trypanosoma cruzi reveals hybridization and clonal mechanisms driving rapid genome diversification

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.

Differential expression of meiosis and homologous recombination-related genes in the life cycle of Trypanosoma cruzi

Trypanosoma cruzi has a complex life cycle consisting of four morphological and distinct biological stages. Although some authors suggest that T. cruzi primarily follows clonal reproduction, recent genomic and transcriptomic studies indicate an unorthodox capacity for recombination. We aimed to estimate the differential gene expression of 10 meiosis/homologous recombination-related genes during the T. cruzi life cycle, including epimastigotes, under two different types of stress (oxidative stress and pH changes). We performed RT-qPCR tests using novel-designed primers to estimate the differential gene expression (∆Ct and ∆∆Ct) of nine genes (SPO11, HAP2, RAD50, MRN complex, BRCA2, DMC1, MND1, and RPA1) and RAD51, which was previously reported. Our results show basal expression of all genes during the life cycle, indicating their hypothetical role in several cellular processes but with specific signatures of differential gene expression during the life cycle (HAP2, RPA, RAD50, BRCA2, MND1, and DMC1) and oxidative stress (RPA, MRE11, NBS1, BRCA2, MND1, and RAD51). Additionally, we found that the MRN complex has an independent level of expression in T. cruzi, with profiles of MRE11 and NBS1 upregulated in some stages. Recent studies on other trypanosomatids have highlighted the influence of HAP2 and RPA in recombination and hybridization. If T. cruzi uses the same repertoire of genes, our findings could suggest that metacyclogenesis may be the putative step that the parasite uses to undergo recombination. Likewise, our study reveals the differential profiles of genes expressed in response to oxidative and pH stress. Further studies are necessary to confirm our findings and understand the recombination mechanism in T. cruzi.

Life in plastic, it's fantastic! How Leishmania exploit genome instability to shape gene expression

Leishmania are kinetoplastid pathogens that cause leishmaniasis, a debilitating and potentially life-threatening infection if untreated. Unusually, Leishmania regulate their gene expression largely post-transcriptionally due to the arrangement of their coding genes into polycistronic transcription units that may contain 100s of functionally unrelated genes. Yet, Leishmania are capable of rapid and responsive changes in gene expression to challenging environments, often instead correlating with dynamic changes in their genome composition, ranging from chromosome and gene copy number variations to the generation of extrachromosomal DNA and the accumulation of point mutations. Typically, such events indicate genome instability in other eukaryotes, coinciding with genetic abnormalities, but for Leishmania, exploiting these products of genome instability can provide selectable substrates to catalyse necessary gene expression changes by modifying gene copy number. Unorthodox DNA replication, DNA repair, replication stress factors and DNA repeats are recognised in Leishmania as contributors to this intrinsic instability, but how Leishmania regulate genome plasticity to enhance fitness whilst limiting toxic under- or over-expression of co-amplified and co-transcribed genes is unclear. Herein, we focus on fresh, and detailed insights that improve our understanding of genome plasticity in Leishmania. Furthermore, we discuss emerging models and factors that potentially circumvent regulatory issues arising from polycistronic transcription. Lastly, we highlight key gaps in our understanding of Leishmania genome plasticity and discuss future studies to define, in higher resolution, these complex regulatory interactions.

Accessing the Variability of Multicopy Genes in Complex Genomes using Unassembled Next-Generation Sequencing Reads: The Case of Trypanosoma cruzi Multigene Families

Repetitive elements cause assembly fragmentation in complex eukaryotic genomes, limiting the study of their variability. The genome of Trypanosoma cruzi, the parasite that causes Chagas disease, has a high repetitive content, including multigene families. Although many T. cruzi multigene families encode surface proteins that play pivotal roles in host-parasite interactions, their variability is currently underestimated, as their high repetitive content results in collapsed gene variants. To estimate sequence variability and copy number variation of multigene families, we developed a read-based approach that is independent of gene-specific read mapping and de novo assembly. This methodology was used to estimate the copy number and variability of MASP, TcMUC, and Trans-Sialidase (TS), the three largest T. cruzi multigene families, in 36 strains, including members of all six parasite discrete typing units (DTUs). We found that these three families present a specific pattern of variability and copy number among the distinct parasite DTUs. Inter-DTU hybrid strains presented a higher variability of these families, suggesting that maintaining a larger content of their members could be advantageous. In addition, in a chronic murine model and chronic Chagasic human patients, the immune response was focused on TS antigens, suggesting that targeting TS conserved sequences could be a potential avenue to improve diagnosis and vaccine design against Chagas disease. Finally, the proposed approach can be applied to study multicopy genes in any organism, opening new avenues to access sequence variability in complex genomes. IMPORTANCE Sequences that have several copies in a genome, such as multicopy-gene families, mobile elements, and microsatellites, are among the most challenging genomic segments to study. They are frequently underestimated in genome assemblies, hampering the correct assessment of these important players in genome evolution and adaptation. Here, we developed a new methodology to estimate variability and copy numbers of repetitive genomic regions and employed it to characterize the T. cruzi multigene families MASP, TcMUC, and transsialidase (TS), which are important virulence factors in this parasite. We showed that multigene families vary in sequence and content among the parasite's lineages, whereas hybrid strains have a higher sequence variability that could be advantageous to the parasite's survivability. By identifying conserved sequences within multigene families, we showed that the mammalian host immune response toward these multigene families is usually focused on the TS multigene family. These TS conserved and immunogenic peptides can be explored in future works as diagnostic targets or vaccine candidates for Chagas disease. Finally, this methodology can be easily applied to any organism of interest, which will aid in our understanding of complex genomic regions.

Genetic diversity of Trypanosoma cruzi strains isolated from chronic chagasic patients and non-human hosts in the state of São Paulo, Brazil

BACKGROUND

Trypanosoma cruzi shows an exuberant genetic diversity. Currently, seven phylogenetic lineages, called discrete typing units (DTUs), are recognised: TcI-TcVI and Tcbat. Despite advances in studies on T. cruzi and its populations, there is no consensus regarding its heterogeneity.

OBJECTIVES

This study aimed to perform molecular characterisation of T. cruzi strains, isolated in the state of São Paulo, to identify the DTUs involved and evaluate their genetic diversity.

METHODS

T. cruzi strains were isolated from biological samples of chronic chagasic patients, marsupials and triatomines through culture techniques and subjected to molecular characterisation using the fluorescent fragment length barcoding (FFLB) technique. Subsequently, the results were correlated with complementary information to enable better discrimination between the identified DTUs.

FINDINGS

It was possible to identify TcI in two humans and two triatomines; TcII/VI in 19 humans, two marsupials and one triatomine; and TcIII in one human host, an individual that also presented a result for TcI, which indicated the possibility of a mixed infection. Regarding the strains characterised by the TcII/VI profile, the correlation with complementary information allowed to suggest that, in general, these parasite populations indeed correspond to the TcII genotype.

MAIN CONCLUSIONS

The TcII/VI profile, associated with domestic cycles and patients with chronic Chagas disease, was the most prevalent among the identified DTUs. Furthermore, the correlation of the study results with complementary information made it possible to suggest that TcII is the predominant lineage of this work.