The Trypanosoma cruzi Sylvio X10 strain maxicircle sequence: the third musketeer

Mitochondrial DNA Structure in Trypanosoma cruzi

Kinetoplastids display a single, large mitochondrion per cell, with their mitochondrial DNA referred to as the kinetoplast. This kinetoplast is a network of concatenated circular molecules comprising a maxicircle (20-64 kb) and up to thousands of minicircles varying in size depending on the species (0.5-10 kb). In Trypanosoma cruzi, maxicircles contain typical mitochondrial genes found in other eukaryotes. They consist of coding and divergent/variable regions, complicating their assembly due to repetitive elements. However, next-generation sequencing (NGS) methods have resolved these issues, enabling the complete sequencing of maxicircles from different strains. Furthermore, several insertions and deletions in the maxicircle sequences have been identified among strains, affecting specific genes. Unique to kinetoplastids, minicircles play a crucial role in a particular U-insertion/deletion RNA editing system by encoding guide RNAs (gRNAs). These gRNAs are essential for editing and maturing maxicircle mRNAs. In Trypanosoma cruzi, although only a few studies have utilized NGS methods to date, the structure of these molecules suggests a classification into four main groups of minicircles. This classification is based on their size and the number of highly conserved regions (mHCRs) and hypervariable regions (mHVRs).

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.

In silico investigation of cytochrome bc1 molecular inhibition mechanism against Trypanosoma cruzi

Chagas' disease is a neglected tropical disease caused by the kinetoplastid protozoan Trypanosoma cruzi. The only therapies are the nitroheterocyclic chemicals nifurtimox and benznidazole that cause various adverse effects. The need to create safe and effective medications to improve medical care remains critical. The lack of verified T. cruzi therapeutic targets hinders medication research for Chagas' disease. In this respect, cytochrome bc1 has been identified as a promising therapeutic target candidate for antibacterial medicines of medical and agricultural interest. Cytochrome bc1 belongs to the mitochondrial electron transport chain and transfers electrons from ubiquinol to cytochrome c1 by the action of two catalytic sites named Qi and Qo. The two binding sites are highly selective, and specific inhibitors exist for each site. Recent studies identified the Qi site of the cytochrome bc1 as a promising drug target against T. cruzi. However, a lack of knowledge of the drug mechanism of action unfortunately hinders the development of new therapies. In this context, knowing the cause of binding site selectivity and the mechanism of action of inhibitors and substrates is crucial for drug discovery and optimization processes. In this paper, we provide a detailed computational investigation of the Qi site of T. cruzi cytochrome b to shed light on the molecular mechanism of action of known inhibitors and substrates. Our study emphasizes the action of inhibitors at the Qi site on a highly unstructured portion of cytochrome b that could be related to the biological function of the electron transport chain complex.

Trypanosoma cruzi strain and starvation-driven mitochondrial RNA editing and transcriptome variability

Trypanosoma cruzi is a unicellular protistan parasitic species that is comprised of strains and isolates exhibiting high levels of genetic and metabolic variability. In the insect vector, it is known to be highly responsive to starvation, a signal for progression to a life stage in which it can infect mammalian cells. Most mRNAs encoded in its mitochondrion require the targeted insertion and deletion of uridines to become translatable transcripts. This study defined differences in uridine-insertion/deletion RNA editing among three strains and established the mechanism whereby abundances of edited (and, thus, translatable) mitochondrial gene products increase during starvation. Our approach utilized our custom T-Aligner toolkit to describe transcriptome-wide editing events and reconstruct editing products from high-throughput sequencing data. We found that the relative abundance of mitochondrial transcripts and the proportion of mRNAs that are edited varies greatly between analyzed strains, a characteristic that could potentially impact metabolic capacity. Starvation typically led to an increase in overall editing activity rather than affecting a specific step in the process. We also determined that transcripts CR3, CR4, and ND3 produce multiple open reading frames that, if translated, would generate different proteins. Finally, we quantitated the inherent flexibility of editing in T. cruzi and found it to be higher relative to that in a related trypanosomatid lineage. Over time, new editing domains or patterns could prove advantageous to the organism and become more widespread within individual transcriptomes or among strains.

Maxicircle architecture and evolutionary insights into Trypanosoma cruzi complex

We sequenced maxicircles from T. cruzi strains representative of the species evolutionary diversity by using long-read sequencing, which allowed us to uncollapse their repetitive regions, finding that their real lengths range from 35 to 50 kb. T. cruzi maxicircles have a common architecture composed of four regions: coding region (CR), AT-rich region, short (SR) and long repeats (LR). Distribution of genes, both in order and in strand orientation are conserved, being the main differences the presence of deletions affecting genes coding for NADH dehydrogenase subunits, reinforcing biochemical findings that indicate that complex I is not functional in T. cruzi. Moreover, the presence of complete minicircles into maxicircles of some strains lead us to think about the origin of minicircles. Finally, a careful phylogenetic analysis was conducted using coding regions of maxicircles from up to 29 strains, and 1108 single copy nuclear genes from all of the DTUs, clearly establishing that taxonomically T. cruzi is a complex of species composed by group 1 that contains clades A (TcI), B (TcIII) and D (TcIV), and group 2 (1 and 2 do not coincide with groups I and II described decades ago) containing clade C (TcII), being all hybrid strains of the BC type. Three variants of maxicircles exist in T. cruzi: a, b and c, in correspondence with clades A, B, and C from mitochondrial phylogenies. While A and C carry maxicircles a and c respectively, both clades B and D carry b maxicircle variant; hybrid strains also carry the b- variant. We then propose a new nomenclature that is self-descriptive and makes use of both the phylogenetic relationships and the maxicircle variants present in T. cruzi.

The Complete Mitochondrial DNA of Trypanosoma cruzi: Maxicircles and Minicircles

The mitochondrial DNA of Trypanosomatids, known as the kinetoplast DNA or kDNA or mtDNA, consists of a few maxicircles and thousands of minicircles concatenated together into a huge complex network. These structures present species-specific sizes, from 20 to 40 Kb in maxicircles and from 0.5 to 10 Kb in minicircles. Maxicircles are equivalent to other eukaryotic mitochondrial DNAs, while minicircles contain coding guide RNAs involved in U-insertion/deletion editing processes exclusive of Trypanosomatids that produce the maturation of the maxicircle-encoded transcripts. The knowledge about this mitochondrial genome is especially relevant since the expression of nuclear and mitochondrial genes involved in oxidative phosphorylation must be coordinated. In Trypanosoma cruzi (T. cruzi), the mtDNA has a dual relevance; the production of energy, and its use as a phylogenetic marker due to its high conservation among strains. Therefore, this study aimed to assemble, annotate, and analyze the complete repertoire of maxicircle and minicircle sequences of different T. cruzi strains by using DNA sequencing. We assembled and annotated the complete maxicircle sequence of the Y and Bug2148 strains. For Bug2148, our results confirm that the maxicircle sequence is the longest assembled to date, and is composed of 21 genes, most of them conserved among Trypanosomatid species. In agreement with previous results, T. cruzi minicircles show a conserved structure around 1.4 Kb, with four highly conserved regions and other four hypervariable regions interspersed between them. However, our results suggest that the parasite minicircles display several sizes and numbers of conserved and hypervariable regions, contrary to those previous studies. Besides, this heterogeneity is also reflected in the three conserved sequence blocks of the conserved regions that play a key role in the minicircle replication. Our results using sequencing technologies of second and third-generation indicate that the different consensus sequences of the maxicircles and minicircles seem to be more complex than previously described indicating at least four different groups in T. cruzi minicircles.

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.

Common Structural Patterns in the Maxicircle Divergent Region of Trypanosomatidae

Maxicircles of all kinetoplastid flagellates are functional analogs of mitochondrial genome of other eukaryotes. They consist of two distinct parts, called the coding region and the divergent region (DR). The DR is composed of highly repetitive sequences and, as such, remains the least explored segment of a trypanosomatid genome. It is extremely difficult to sequence and assemble, that is why very few full length maxicircle sequences were available until now. Using PacBio data, we assembled 17 complete maxicircles from different species of trypanosomatids. Here we present their large-scale comparative analysis and describe common patterns of DR organization in trypanosomatids.

A link between mitochondrial gene expression and life stage morphologies in Trypanosoma cruzi

The protozoan Trypanosoma cruzi has a complicated dual-host life cycle, and starvation can trigger transition from the replicating insect stage to the mammalian-infectious nonreplicating insect stage (epimastigote to trypomastigote differentiation). Abundance of some mature RNAs derived from its mitochondrial genome increase during culture starvation of T. cruzi for unknown reasons. Here, we examine T. cruzi mitochondrial gene expression in the mammalian intracellular replicating life stage (amastigote), and uncover implications of starvation-induced changes in gene expression. Mitochondrial RNA levels in general were found to be lowest in actively replicating amastigotes. We discovered that mitochondrial respiration decreases during starvation in insect stage cells, despite the previously observed increases in mitochondrial mRNAs encoding electron transport chain (ETC) components. Surprisingly, T. cruzi epimastigotes in replete medium grow at normal rates when we genetically compromised their ability to perform insertion/deletion editing and thereby generate mature forms of some mitochondrial mRNAs. However, these cells, when starved, were impeded in the epimastigote to trypomastigote transition. Further, they experience a short-flagella phenotype that may also be linked to differentiation. We hypothesize a scenario where levels of mature RNA species or editing in the single T. cruzi mitochondrion are linked to differentiation by a yet-unknown signaling mechanism.

Trypanosomatid mitochondrial RNA editing: dramatically complex transcript repertoires revealed with a dedicated mapping tool

RNA editing by targeted insertion and deletion of uridine is crucial to generate translatable mRNAs from the cryptogenes of the mitochondrial genome of kinetoplastids. This type of editing consists of a stepwise cascade of reactions generally proceeding from 3′ to 5′ on a transcript, resulting in a population of partially edited as well as pre-edited and completely edited molecules for each mitochondrial cryptogene of these protozoans. Often, the number of uridines inserted and deleted exceed the number of nucleotides that are genome-encoded. Thus, analysis of kinetoplastid mitochondrial transcriptomes has proven frustratingly complex. Here we present our analysis of Leptomonas pyrrhocoris mitochondrial cDNA deep sequencing reads using T-Aligner, our new tool which allows comprehensive characterization of RNA editing, not relying on targeted transcript amplification and on prior knowledge of final edited products. T-Aligner implements a pipeline of read mapping, visualization of all editing states and their coverage, and assembly of canonical and alternative translatable mRNAs. We also assess T-Aligner functionality on a more challenging deep sequencing read input from Trypanosoma cruzi. The analysis reveals that transcripts of cryptogenes of both species undergo very complex editing that includes the formation of alternative open reading frames and whole categories of truncated editing products.

Mitochondrial genomics of human pathogenic parasite Leishmania (Viannia) panamensis

Background

The human parasite Leishmania (V.) panamensis is one of the pathogenic species responsible for cutaneous leishmaniasis in Central and South America. Despite its importance in molecular parasitology, its mitochondrial genome, divided into minicircles and maxicircles, haven’t been described so far.

Methods

Using NGS-based sequencing (454 and ILLUMINA), and combining de novo genome assembly and mapping strategies, we report the maxicircle kDNA annotated genome of L. (V.) panamensis, the first reference of this molecule for the subgenus Viannia. A comparative genomics approach is performed against other Leishmania and Trypanosoma species.

Results

The results show synteny of mitochondrial genes of L. (V.) panamensis with other kinetoplastids. It was also possible to identify nucleotide variants within the coding regions of the maxicircle, shared among some of them and others specific to each strain. Furthermore, we compared the minicircles kDNA sequences of two strains and the results show that the conserved and divergent regions of the minicircles exhibit strain-specific associations.

High throughput sequencing revolution reveals conserved fundamentals of U-indel editing

Among Euglenozoans, mitochondrial RNA editing occurs in the diplonemids and in the kinetoplastids that include parasitic trypanosomes. Yet U-indel editing, in which open reading frames (ORFs) on mRNAs are generated by insertion and deletion of uridylates in locations dictated by guide RNAs, appears confined to kinetoplastids. The nature of guide RNA and edited mRNA populations has been cursorily explored in a surprisingly extensive number of species over the years, although complete sets of fully edited mRNAs for most kinetoplast genomes are largely missing. Now, however, high throughput sequencing technologies have had an enormous impact on what we know and will learn about the mechanisms, benefits, and final edited products of U-indel editing. Tools including PARERS, TREAT, and T-Aligner function to organize and make sense of U-indel mRNA transcriptomes, which are comprised of mRNAs harboring uridylate indels both consistent and inconsistent with translatable products. From high throughput sequencing data come arguments that partially edited mRNAs containing "junction regions" of noncanonical editing are editing intermediates, and conversely, arguments that they are dead-end products. These data have also revealed that the percent of a given transcript population that is fully or partially edited varies dramatically between transcripts and organisms. Outstanding questions that are being addressed include the prevalence of sequences that apparently encode alternative ORFs, diversity of editing events in ORF termini and 5' and 3' untranslated regions, and the differences that exist in this byzantine process between species. High throughput sequencing technologies will also undoubtedly be harnessed to probe U-indel editing's evolutionary origins. This article is categorized under: RNA Processing > RNA Editing and Modification RNA Evolution and Genomics > Computational Analyses of RNA.

Gene expression to mitochondrial metabolism: Variability among cultured Trypanosoma cruzi strains

The insect-transmitted protozoan parasite Trypanosoma cruzi experiences changes in nutrient availability and rate of flux through different metabolic pathways across its life cycle. The species encompasses much genetic diversity of both the nuclear and mitochondrial genomes among isolated strains. The genetic or expression variation of both genomes are likely to impact metabolic responses to environmental stimuli, and even steady state metabolic function, among strains. To begin formal characterization these differences, we compared aspects of metabolism between genetically similar strains CL Brener and Tulahuen with less similar Esmeraldo and Sylvio X10 strains in a culture environment. Epimastigotes of all strains took up glucose at similar rates. However, the degree of medium acidification that could be observed when glucose was absent from the medium varied by strain, indicating potential differences in excreted metabolic byproducts. Our main focus was differences related to electron transport chain function. We observed differences in ATP-coupled respiration and maximal respiratory capacity, mitochondrial membrane potential, and mitochondrial morphology between strains, despite the fact that abundances of two nuclear-encoded proteins of the electron transport chain are similar between strains. RNA sequencing reveals strain-specific differences in abundances of mRNAs encoding proteins of the respiratory chain but also other metabolic processes. From these differences in metabolism and mitochondrial phenotypes we have generated tentative models for the differential metabolic fluxes or differences in gene expression that may underlie these results.

The kinetoplast DNA of the Australian trypanosome, Trypanosoma copemani, shares features with Trypanosoma cruzi and Trypanosoma lewisi

Kinetoplast DNA (kDNA) is the mitochondrial genome of trypanosomatids. It consists of a few dozen maxicircles and several thousand minicircles, all catenated topologically to form a two-dimensional DNA network. Minicircles are heterogeneous in size and sequence among species. They present one or several conserved regions that contain three highly conserved sequence blocks. CSB-1 (10 bp sequence) and CSB-2 (8 bp sequence) present lower interspecies homology, while CSB-3 (12 bp sequence) or the Universal Minicircle Sequence is conserved within most trypanosomatids. The Universal Minicircle Sequence is located at the replication origin of the minicircles, and is the binding site for the UMS binding protein, a protein involved in trypanosomatid survival and virulence. Here, we describe the structure and organisation of the kDNA of Trypanosoma copemani, a parasite that has been shown to infect mammalian cells and has been associated with the drastic decline of the endangered Australian marsupial, the woylie (Bettongia penicillata). Deep genomic sequencing showed that T. copemani presents two classes of minicircles that share sequence identity and organisation in the conserved sequence blocks with those of Trypanosoma cruzi and Trypanosoma lewisi. A 19,257 bp partial region of the maxicircle of T. copemani that contained the entire coding region was obtained. Comparative analysis of the T. copemani entire maxicircle coding region with the coding regions of T. cruzi and T. lewisi showed they share 71.05% and 71.28% identity, respectively. The shared features in the maxicircle/minicircle organisation and sequence between T. copemani and T. cruzi/T. lewisi suggest similarities in their process of kDNA replication, and are of significance in understanding the evolution of Australian trypanosomes.

Analysis of the mitochondrial maxicircle of Trypanosoma lewisi, a neglected human pathogen

BACKGROUND

The haemoflagellate Trypanosoma lewisi is a kinetoplastid parasite which, as it has been recently reported to cause human disease, deserves increased attention. Characteristic features of all kinetoplastid flagellates are a uniquely structured mitochondrial DNA or kinetoplast, comprised of a network of catenated DNA circles, and RNA editing of mitochondrial transcripts. The aim of this study was to describe the kinetoplast DNA of T. lewisi.

METHODS/RESULTS

In this study, purified kinetoplast DNA from T. lewisi was sequenced using high-throughput sequencing in combination with sequencing of PCR amplicons. This allowed the assembly of the T. lewisi kinetoplast maxicircle DNA, which is a homologue of the mitochondrial genome in other eukaryotes. The assembly of 23,745 bp comprises the non-coding and coding regions. Comparative analysis of the maxicircle sequence of T. lewisi with Trypanosoma cruzi, Trypanosoma rangeli, Trypanosoma brucei and Leishmania tarentolae revealed that it shares 78%, 77%, 74% and 66% sequence identity with these parasites, respectively. The high GC content in at least 9 maxicircle genes of T. lewisi (ATPase6; NADH dehydrogenase subunits ND3, ND7, ND8 and ND9; G-rich regions GR3 and GR4; cytochrome oxidase subunit COIII and ribosomal protein RPS12) implies that their products may be extensively edited. A detailed analysis of the non-coding region revealed that it contains numerous repeat motifs and palindromes.

CONCLUSIONS

We have sequenced and comprehensively annotated the kinetoplast maxicircle of T. lewisi. Our analysis reveals that T. lewisi is closely related to T. cruzi and T. brucei, and may share similar RNA editing patterns with them rather than with L. tarentolae. These findings provide novel insight into the biological features of this emerging human pathogen.

Comparative genomic analysis of human infective Trypanosoma cruzi lineages with the bat-restricted subspecies T. cruzi marinkellei

Background

Trypanosoma cruzi marinkellei is a bat-associated parasite of the subgenus Schizotrypanum and it is regarded as a T. cruzi subspecies. Here we report a draft genome sequence of T. c. marinkellei and comparison with T. c. cruzi. Our aims were to identify unique sequences and genomic features, which may relate to their distinct niches.

Results

The T. c. marinkellei genome was found to be ~11% smaller than that of the human-derived parasite T. c. cruzi Sylvio X10. The genome size difference was attributed to copy number variation of coding and non-coding sequences. The sequence divergence in coding regions was ~7.5% between T. c. marinkellei and T. c. cruzi Sylvio X10. A unique acetyltransferase gene was identified in T. c. marinkellei, representing an example of a horizontal gene transfer from eukaryote to eukaryote. Six of eight examined gene families were expanded in T. c. cruzi Sylvio X10. The DGF gene family was expanded in T. c. marinkellei. T. c. cruzi Sylvio X10 contained ~1.5 fold more sequences related to VIPER and L1Tc elements. Experimental infections of mammalian cell lines indicated that T. c. marinkellei has the capacity to invade non-bat cells and undergo intracellular replication.

Conclusions

Several unique sequences were identified in the comparison, including a potential subspecies-specific gene acquisition in T. c. marinkellei. The identified differences reflect the distinct evolutionary trajectories of these parasites and represent targets for functional investigation.

Multiple mitochondrial introgression events and heteroplasmy in trypanosoma cruzi revealed by maxicircle MLST and next generation sequencing

Background

Mitochondrial DNA is a valuable taxonomic marker due to its relatively fast rate of evolution. In Trypanosoma cruzi, the causative agent of Chagas disease, the mitochondrial genome has a unique structural organization consisting of 20–50 maxicircles (∼20 kb) and thousands of minicircles (0.5–10 kb). T. cruzi is an early diverging protist displaying remarkable genetic heterogeneity and is recognized as a complex of six discrete typing units (DTUs). The majority of infected humans are asymptomatic for life while 30–35% develop potentially fatal cardiac and/or digestive syndromes. However, the relationship between specific clinical outcomes and T. cruzi genotype remains elusive. The availability of whole genome sequences has driven advances in high resolution genotyping techniques and re-invigorated interest in exploring the diversity present within the various DTUs.

Methodology/Principal Findings

To describe intra-DTU diversity, we developed a highly resolutive maxicircle multilocus sequence typing (mtMLST) scheme based on ten gene fragments. A panel of 32 TcI isolates was genotyped using the mtMLST scheme, GPI, mini-exon and 25 microsatellite loci. Comparison of nuclear and mitochondrial data revealed clearly incongruent phylogenetic histories among different geographical populations as well as major DTUs. In parallel, we exploited read depth data, generated by Illumina sequencing of the maxicircle genome from the TcI reference strain Sylvio X10/1, to provide the first evidence of mitochondrial heteroplasmy (heterogeneous mitochondrial genomes in an individual cell) in T. cruzi.

Conclusions/Significance

mtMLST provides a powerful approach to genotyping at the sub-DTU level. This strategy will facilitate attempts to resolve phenotypic variation in T. cruzi and to address epidemiologically important hypotheses in conjunction with intensive spatio-temporal sampling. The observations of both general and specific incidences of nuclear-mitochondrial phylogenetic incongruence indicate that genetic recombination is geographically widespread and continues to influence the natural population structure of TcI, a conclusion which challenges the traditional paradigm of clonality in T. cruzi.

Chagas disease, caused by the protozoan parasite Trypanosoma cruzi, is an important public health problem in Latin America. While molecular techniques can differentiate the major T. cruzi genetic lineages, few have sufficient resolution to describe diversity among closely related strains. The online availability of three mitochondrial genomes allowed us to design a multilocus sequence typing (mtMLST) scheme to exploit these rapidly evolving markers. We compared mtMLST with current nuclear typing tools using isolates belonging to the oldest and most widely occurring lineage TcI. T. cruzi is generally believed to reproduce clonally. However, in this study, distinct branching patterns between mitochondrial and nuclear phylogenetic trees revealed multiple incidences of genetic exchange within different geographical populations and major lineages. We also examined Illumina sequencing data from the TcI genome strain which revealed multiple different mitochondrial genomes within an individual parasite (heteroplasmy) that were, however, not sufficiently divergent to represent a major source of typing error. We strongly recommend this combined nuclear and mitochondrial genotyping methodology to reveal cryptic diversity and genetic exchange in T. cruzi. The level of resolution that this mtMLST provides should greatly assist attempts to elucidate the complex interactions between parasite genotype, clinical outcome and disease distribution.

Trypanosomatid comparative genomics: Contributions to the study of parasite biology and different parasitic diseases

In 2005, draft sequences of the genomes of Trypanosoma brucei, Trypanosoma cruzi and Leishmania major, also known as the Tri-Tryp genomes, were published. These protozoan parasites are the causative agents of three distinct insect-borne diseases, namely sleeping sickness, Chagas disease and leishmaniasis, all with a worldwide distribution. Despite the large estimated evolutionary distance among them, a conserved core of ~6,200 trypanosomatid genes was found among the Tri-Tryp genomes. Extensive analysis of these genomic sequences has greatly increased our understanding of the biology of these parasites and their host-parasite interactions. In this article, we review the recent advances in the comparative genomics of these three species. This analysis also includes data on additional sequences derived from other trypanosmatid species, as well as recent data on gene expression and functional genomics. In addition to facilitating the identification of key parasite molecules that may provide a better understanding of these complex diseases, genome studies offer a rich source of new information that can be used to define potential new drug targets and vaccine candidates for controlling these parasitic infections.

The revised Trypanosoma cruzi subspecific nomenclature: rationale, epidemiological relevance and research applications

The protozoan Trypanosoma cruzi, its mammalian reservoirs, and vectors have existed in nature for millions of years. The human infection, named Chagas disease, is a major public health problem for Latin America. T. cruzi is genetically highly diverse and the understanding of the population structure of this parasite is critical because of the links to transmission cycles and disease. At present, T. cruzi is partitioned into six discrete typing units (DTUs), TcI-TcVI. Here we focus on the current status of taxonomy-related areas such as population structure, phylogeographical and eco-epidemiological features, and the correlation of DTU with natural and experimental infection. We also summarize methods for DTU genotyping, available for widespread use in endemic areas. For the immediate future multilocus sequence typing is likely to be the gold standard for population studies. We conclude that greater advances in our knowledge on pathogenic and epidemiological features of these parasites are expected in the coming decade through the comparative analysis of the genomes from isolates of various DTUs.