Maxicircle architecture and evolutionary insights into Trypanosoma cruzi complex

High-throughput analysis of the Trypanosoma cruzi minicirculome (mcDNA) unveils structural variation and functional diversity

Trypanosoma cruzi causes Chagas disease and has a unique extranuclear genome enclosed in a structure called the kinetoplast, which contains circular genomes known as maxi- and minicircles. While the structure and function of maxicircles are well-understood, many aspects of minicircles remain to be discovered. Here, we performed a high-throughput analysis of the minicirculome (mcDNA) in 50 clones isolated from Colombia's diverse T. cruzi I populations. Results indicate that mcDNA comprises four diverse subpopulations with different structures, lengths, and numbers of interspersed semi-conserved (previously termed ultra-conserved regions mHCV) and hypervariable (mHVPs) regions. Analysis of mcDNA ancestry and inter-clone differentiation indicates the interbreeding of minicircle sequence classes is placed along diverse strains and hosts. These results support evidence of the multiclonal dynamics and random bi-parental segregation. Finally, we disclosed the guide RNA repertoire encoded by mcDNA at a clonal scale, and several attributes of its abundance and function are discussed.

A phased genome assembly of a Colombian Trypanosoma cruzi TcI strain and the evolution of gene families

Chagas is an endemic disease in tropical regions of Latin America, caused by the parasite Trypanosoma cruzi. High intraspecies variability and genome complexity have been challenges to assemble high quality genomes needed for studies in evolution, population genomics, diagnosis and drug development. Here we present a chromosome-level phased assembly of a TcI T. cruzi strain (Dm25). While 29 chromosomes show a large collinearity with the assembly of the Brazil A4 strain, three chromosomes show both large heterozygosity and large divergence, compared to previous assemblies of TcI T. cruzi strains. Nucleotide and protein evolution statistics indicate that T. cruzi Marinkellei separated before the diversification of T. cruzi in the known DTUs. Interchromosomal paralogs of dispersed gene families and histones appeared before but at the same time have a more strict purifying selection, compared to other repeat families. Previously unreported large tandem arrays of protein kinases and histones were identified in this assembly. Over one million variants obtained from Illumina reads aligned to the primary assembly clearly separate the main DTUs. We expect that this new assembly will be a valuable resource for further studies on evolution and functional genomics of Trypanosomatids.

Identification of a unique conserved region from a kinetoplastid genome of Leishmania orientalis (formerly named Leishmania siamensis) strain PCM2 in Thailand

Mitochondrial DNAs (mtDNAs) appear in almost all eukaryotic species and are useful molecular markers for phylogenetic studies and species identification. Kinetoplast DNAs (kDNAs) are structurally complex circular mtDNA networks in kinetoplastids, divided into maxicircles and minicircles. Despite several kDNAs of many Leishmania species being examined, the kDNAs of the new species, Leishmania orientalis (formerly named Leishmania siamensis) strain PCM2, have not been explored. This study aimed to investigate the maxicircle and minicircle DNAs of L. orientalis strain PCM2 using hybrid genome sequencing technologies and bioinformatic analyses. The kDNA sequences were isolated and assembled using the SPAdes hybrid assembler from the Illumina short-read and PacBio long-read data. Circular contigs of the maxicircle and minicircle DNAs were reconstructed and confirmed by BLASTn and rKOMICs programs. The kDNA genome was annotated by BLASTn before the genome comparison and phylogenetic analysis by progressiveMauve, MAFFT, and MEGA programs. The maxicircle of L. orientalis strain PCM2 (18,215 bp) showed 99.92% similarity and gene arrangement to Leishmania enriettii strain LEM3045 maxicircle with variation in the 12s rRNA gene and divergent region. Phylogenetics of the whole sequence, coding regions, divergent regions, and 12s rRNA gene also confirmed this relationship and subgenera separation. The identified 105 classes of minicircles (402-1177 bp) were clustered monophyletically and related to the Leishmania donovani minicircles. The kinetoplast maxicircle and minicircle DNAs of L. orientalis strain PCM2 contained a unique conserved region potentially useful for specific diagnosis of L. orientalis and further exploration of this parasite population genetics in Thailand and related regions.

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.

Comparative microRNA profiling of Trypanosoma cruzi infected human cells

Introduction

Trypanosoma cruzi, the causative agent of Chagas disease, can infect almost any nucleated cell in the mammalian host. Although previous studies have described the transcriptomic changes that occur in host cells during parasite infection, the understanding of the role of post-transcriptional regulation in this process is limited. MicroRNAs, a class of short non-coding RNAs, are key players in regulating gene expression at the post-transcriptional level, and their involvement in the host-T. cruzi interplay is a growing area of research. However, to our knowledge, there are no comparative studies on the microRNA changes that occur in different cell types in response to T. cruzi infection.

Methods and results

Here we investigated microRNA changes in epithelial cells, cardiomyocytes and macrophages infected with T. cruzi for 24 hours, using small RNA sequencing followed by careful bioinformatics analysis. We show that, although microRNAs are highly cell type-specific, a signature of three microRNAs -miR-146a, miR-708 and miR-1246, emerges as consistently responsive to T. cruzi infection across representative human cell types. T. cruzi lacks canonical microRNA-induced silencing mechanisms and we confirm that it does not produce any small RNA that mimics known host microRNAs. We found that macrophages show a broad response to parasite infection, while microRNA changes in epithelial and cardiomyocytes are modest. Complementary data indicated that cardiomyocyte response may be greater at early time points of infection.

Conclusions

Our findings emphasize the significance of considering microRNA changes at the cellular level and complement previous studies conducted at higher organizational levels, such as heart samples. While miR-146a has been previously implicated in T. cruzi infection, similarly to its involvement in many other immunological responses, miR-1246 and miR-708 are demonstrated here for the first time. Given their expression in multiple cell types, we anticipate our work as a starting point for future investigations into their role in the post-transcriptional regulation of T. cruzi infected cells and their potential as biomarkers for Chagas disease.

Trypanosoma cruzi Isolates Naturally Adapted to Congenital Transmission Display a Unique Strategy of Transplacental Passage

Chagas disease is mainly transmitted by vertical transmission (VT) in nonendemic areas and in endemic areas where vector control programs have been successful. For the present study, we isolated natural Trypanosoma cruzi strains vertically transmitted through three generations and proceeded to study their molecular mechanism of VT using mice. No parasitemia was detected in immunocompetent mice, but the parasites were able to induce an immune response and colonize different organs. VT experiments revealed that infection with different strains did not affect mating, pregnancy, or resorption, but despite low parasitemia, VT strains reached the placenta and resulted in higher vertical transmission rates than strains of either moderate or high virulence. While the virulent strain modulated more than 2,500 placental genes, VT strains modulated 150, and only 29 genes are shared between them. VT strains downregulated genes associated with cell division and replication and upregulated immunomodulatory genes, leading to anti-inflammatory responses and tolerance. The virulent strain stimulated a strong proinflammatory immune response, and this molecular footprint correlated with histopathological analyses. We describe a unique placental response regarding the passage of T. cruzi VT isolates across the maternal-fetal interphase, challenging the current knowledge derived mainly from studies of laboratory-adapted or highly virulent strains. IMPORTANCE The main findings of this study are that we determined that there are Trypanosoma cruzi strains adapted to transplacental transmission and completely different from the commonly used laboratory reference strains. This implies a specific strategy for the vertical transmission of Chagas disease. It is impressive that the strains specialized for vertical transmission modify the gene expression of the placenta in a totally different way than the reference strains. In addition, we describe isolates of T. cruzi that cannot be transmitted transplacentally. Taken together, these results open up new insights into the molecular mechanisms of this insect vector-independent transmission form.

Novel insertions in the mitochondrial maxicircle of Trypanosoma musculi, a mouse trypanosome

Trypanosoma musculi is a, globally distributed, mouse-specific haemoflagellate, of the family Trypanosomatidae, which shares similar characteristics in morphology with Trypanosoma lewisi. The kinetoplast (mitochondrial) DNA of Trypanosomatidae flagellates is comprised of catenated maxicircles and minicircles. However, genetic information on the T. musculi kinetoplast remains largely unknown. In this study, the T. musculi maxicircle genome was completely assembled, with PacBio and Illumina sequencing, and the size was confirmed at 34 606 bp. It consisted of 2 distinct parts: the coding region and the divergent regions (DRs, DRI and II). In comparison with other trypanosome maxicircles (Trypanosoma brucei, Trypanosoma cruzi and T. lewisi), the T. musculi maxicircle has a syntenic distribution of genes and shares 73.9, 78.0 and 92.7% sequence identity, respectively, over the whole coding region. Moreover, novel insertions in MURF2 (630 bp) and in ND5 (1278 bp) were found, respectively, which are homologous to minicircles. These findings support an evolutionary scenario similar to the one proposed for insertions in Trypanosoma cruzi, the pathogen of American trypanosomiasis. These novel insertions, together with a deletion (281 bp) in ND4, question the role of Complex I in T. musculi. A detailed analysis of DRII indicated that it contains numerous repeat motifs and palindromes, the latter of which are highly conservative and contain A5C elements. The comprehensively annotated kinetoplast maxicircle of T. musculi reveals a high degree of similarity between this parasite and the maxicircle of T. lewisi and suggests that the DRII could be a valuable marker for distinguishing these evolutionarily related species.

Genetic Diversity of Trypanosoma cruzi in Panama Inferred by Multi-locus Sequence Typing of Mitochondrial Genes

The objective of this study was to provide information on Trypanosoma cruzi genetic diversity among isolates obtained from different biological sources circulating in endemic areas of Panama. Initial discrete typing units (DTUs) assignment was performed evaluating three single locus molecular markers (mini-exon, heat shock protein 60 and glucose-6-phosphate isomerase genes). Further diversity within TcI lineages was explored using a multi-locus sequence typing approach with six maxicircle genes. Haplotype network analysis and evolutionary divergency estimations were conducted to investigate the genetic relatedness between Panamanian TcI isolates and isolates from different endemic regions in the Americas. Our molecular approach validated that TcI is the predominant DTU circulating in Panama across different hosts and vector species, but also confirmed the presence of TcIII and TcVI circulating in the country. The phylogenetic tree topography for most Panamanian TcI isolates displayed a high level of genetic homogeneity between them. The haplotype network analysis inferred a higher genetic diversity within Panamanian TcI isolates, displaying eight different haplotypes circulating in endemic regions of the country, and revealed geographical structuring among TcI from different endemic regions in the Americas. This study adds novelty on the genetic diversity of T. cruzi circulating in Panama and complements regional phylogeographic studies regarding intra-TcI variations.