In silico, biologically-inspired modelling of genomic variation generation in surface proteins of Trypanosoma cruzi

The Elusive Trypanosoma cruzi Disperse Gene Protein Family (DGF-1)

Chagas disease, caused by Trypanosoma cruzi infections, is included in the group of neglected diseases, and efforts to develop new therapeutic or immunoprevention approaches have not been successful. After the publication of the T. cruzi genome, the number of molecular and biochemical studies on this parasite has increased considerably, many of which are focused on families of variant surface proteins, especially trans-sialidases, mucins, and mucin-associated proteins. The disperse gene protein 1 family (DGF-1) is one of the most abundant families in the T. cruzi genome; however, the large gene size, high copy numbers, and low antibody titers detected in infected humans make it an unattractive study target. However, here we argue that given the ubiquitous presence in all T. cruzi species, and physicochemical characteristics, the DGF-1 gene family may play and important role in host-parasite interactions.

Genomic Organization and Generation of Genetic Variability in the RHS (Retrotransposon Hot Spot) Protein Multigene Family in Trypanosoma cruzi

Retrotransposon Hot Spot (RHS) is the most abundant gene family in Trypanosoma cruzi, with unknown function in this parasite. The aim of this work was to shed light on the organization and expression of RHS in T. cruzi. The diversity of the RHS protein family in T. cruzi was demonstrated by phylogenetic and recombination analyses. Transcribed sequences carrying the RHS domain were classified into ten distinct groups of monophyletic origin. We identified numerous recombination events among the RHS and traced the origins of the donors and target sequences. The transcribed RHS genes have a mosaic structure that may contain fragments of different RHS inserted in the target sequence. About 30% of RHS sequences are located in the subtelomere, a region very susceptible to recombination. The evolution of the RHS family has been marked by many events, including gene duplication by unequal mitotic crossing-over, homologous, as well as ectopic recombination, and gene conversion. The expression of RHS was analyzed by immunofluorescence and immunoblotting using anti-RHS antibodies. RHS proteins are evenly distributed in the nuclear region of T. cruzi replicative forms (amastigote and epimastigote), suggesting that they could be involved in the control of the chromatin structure and gene expression, as has been proposed for T. brucei.

Comparative Analysis of the Secretome and Interactome of Trypanosoma cruzi and Trypanosoma rangeli Reveals Species Specific Immune Response Modulating Proteins

Chagas disease, a zoonosis caused by the flagellate protozoan Trypanosoma cruzi, is a chronic and systemic parasitic infection that affects ~5–7 million people worldwide, mainly in Latin America. Chagas disease is an emerging public health problem due to the lack of vaccines and effective treatments. According to recent studies, several T. cruzi secreted proteins interact with the human host during cell invasion. Moreover, some comparative studies with T. rangeli, which is non-pathogenic in humans, have been performed to identify proteins directly involved in the pathogenesis of the disease. In this study, we present an integrated analysis of canonical putative secreted proteins (PSPs) from both species. Additionally, we propose an interactome with human host and gene family clusters, and a phylogenetic inference of a selected protein. In total, we identified 322 exclusively PSPs in T. cruzi and 202 in T. rangeli. Among the PSPs identified in T. cruzi, we found several trans-sialidases, mucins, MASPs, proteins with phospholipase 2 domains (PLA2-like), and proteins with Hsp70 domains (Hsp70-like) which have been previously characterized and demonstrated to be related to T. cruzi virulence. PSPs found in T. rangeli were related to protozoan metabolism, specifically carboxylases and phosphatases. Furthermore, we also identified PSPs that may interact with the human immune system, including heat shock and MASP proteins, but in a lower number compared to T. cruzi. Interestingly, we describe a hypothetical hybrid interactome of PSPs which reveals that T. cruzi secreted molecules may be down-regulating IL-17 whilst T. rangeli may enhance the production of IL-15. These results will pave the way for a better understanding of the pathophysiology of Chagas disease and may ultimately lead to the identification of molecular targets, such as key PSPs, that could be used to minimize the health outcomes of Chagas disease by modulating the immune response triggered by T. cruzi infection.

Trypanosoma cruzi Genome 15 Years Later: What Has Been Accomplished?

On 15 July 2020 was the 15th anniversary of the Science Magazine issue that reported three trypanosomatid genomes, namely Leishmania major, Trypanosoma brucei, and Trypanosoma cruzi. That publication was a milestone for the research community working with trypanosomatids, even more so, when considering that the first draft of the human genome was published only four years earlier after 15 years of research. Although nowadays, genome sequencing has become commonplace, the work done by researchers before that publication represented a huge challenge and a good example of international cooperation. Research in neglected diseases often faces obstacles, not only because of the unique characteristics of each biological model but also due to the lower funds the research projects receive. In the case of Trypanosoma cruzi the etiologic agent of Chagas disease, the first genome draft published in 2005 was not complete, and even after the implementation of more advanced sequencing strategies, to this date no final chromosomal map is available. However, the first genome draft enabled researchers to pick genes a la carte, produce proteins in vitro for immunological studies, and predict drug targets for the treatment of the disease or to be used in PCR diagnostic protocols. Besides, the analysis of the T. cruzi genome is revealing unique features about its organization and dynamics. In this work, I briefly summarize the actions of Latin American researchers that contributed to the first publication of the T. cruzi genome and discuss some features of the genome that may help to understand the parasite's robustness and adaptive capabilities.

An Evolutionary View of Trypanosoma Cruzi Telomeres

Like in most eukaryotes, the linear chromosomes of Trypanosoma cruzi end in a nucleoprotein structure called the telomere, which is preceded by regions of variable length called subtelomeres. Together telomeres and subtelomeres are dynamic sites where DNA sequence rearrangements can occur without compromising essential interstitial genes or chromosomal synteny. Good examples of subtelomeres involvement are the expansion of human olfactory receptors genes, variant surface antigens in Trypanosoma brucei, and Saccharomyces cerevisiae mating types. T. cruzi telomeres are made of long stretches of the hexameric repeat 5′-TTAGGG-OH-3′, and its subtelomeres are enriched in genes and pseudogenes from the large gene families RHS, TS and DGF1, DEAD/H-RNA helicase and N-acetyltransferase, intermingled with sequences of retrotransposons elements. In particular, members of the Trans-sialidase type II family appear to have played a role in shaping the current T. cruzi telomere structure. Although the structure and function of T. cruzi telomeric and subtelomeric regions have been documented, recent experiments are providing new insights into T. cruzi's telomere-subtelomere dynamics. In this review, I discuss the co-evolution of telomere, subtelomeres and the TS gene family, and the role that these regions may have played in shaping T. cruzi's genome.

Assessment of genetic mutation frequency induced by oxidative stress in Trypanosoma cruzi

Trypanosoma cruzi is the etiological agent of Chagas disease, a public health challenge due to its morbidity and mortality rates, which affects around 6-7 million people worldwide. Symptoms, response to chemotherapy, and the course of Chagas disease are greatly influenced by T. cruzi's intra-specific variability. Thus, DNA mutations in this parasite possibly play a key role in the wide range of clinical manifestations and in drug sensitivity. Indeed, the environmental conditions of oxidative stress faced by T. cruzi during its life cycle can generate genetic mutations. However, the lack of an established experimental design to assess mutation rates in T. cruzi precludes the study of conditions and mechanisms that potentially produce genomic variability in this parasite. We developed an assay that employs a reporter gene that, once mutated in specific positions, convert G418-sensitive into G418-insenstitive T. cruzi. We were able to determine the frequency of DNA mutations in T. cruzi exposed and non-exposed to oxidative insults assessing the number of colony-forming units in solid selective media after plating a defined number of cells. We verified that T. cruzi's spontaneous mutation frequency was comparable to those found in other eukaryotes, and that exposure to hydrogen peroxide promoted a two-fold increase in T. cruzi's mutation frequency. We hypothesize that genetic mutations in T. cruzi can arise from oxidative insults faced by this parasite during its life cycle.

Multigene families in Trypanosoma cruzi and their role in infectivity

The Trypanosoma cruzi genome contains the most widely expanded content (∼12,000 genes) of the trypanosomatids sequenced to date. This expansion is reflected in the high number of repetitive sequences and particularly in the large quantity of genes that make up its multigene families. Recently it was discovered that the contents of these families vary between phylogenetically unrelated strains. We review the basic characteristics of trans-sialidases and mucins as part of the mechanisms of immune evasion of T. cruzi and as ligands and factors involved in the cross talk between the host cell and the parasite. We also show recently published data describing two new multigene families, DGF-1 and MASP, that form an important part of the scenario representing the complex biology of T. cruzi.

Pathogenesis of Chagas disease: time to move on

Trypanosoma cruzi is the etiologic agent of Chagas disease. The contributions of parasite and immune system for disease pathogenesis remain unresolved and controversial. The possibility that Chagas disease was an autoimmune progression triggered by T. cruzi infection led some to question the benefit of treating chronically T. cruzi-infected persons with drugs. Furthermore, it provided the rationale for not investing in research aimed at a vaccine which might carry a risk of inducing autoimmunity or exacerbating inflammation. This viewpoint was adopted by cash-strapped health systems in the developing economies where the disease is endemic and has been repeatedly challenged by researchers and clinicians in recent years and there is now a considerable body of evidence and broad consensus that parasite persistence is requisite for pathogenesis and that antiparasitic immunity can be protective against T. cruzi pathogenesis without eliciting autoimmune pathology. Thus, treatment of chronically infected patients is likely to yield positive outcomes and efforts to understand immunity and vaccine development should be recognized as a priority area of research for Chagas disease.

Genomic analyses, gene expression and antigenic profile of the trans-sialidase superfamily of Trypanosoma cruzi reveal an undetected level of complexity

The protozoan parasite Trypanosoma cruzi is the etiologic agent of Chagas disease, a highly debilitating human pathology that affects millions of people in the Americas. The sequencing of this parasite's genome reveals that trans-sialidase/trans-sialidase-like (TcS), a polymorphic protein family known to be involved in several aspects of T. cruzi biology, is the largest T. cruzi gene family, encoding more than 1,400 genes. Despite the fact that four TcS groups are well characterized and only one of the groups contains active trans-sialidases, all members of the family are annotated in the T. cruzi genome database as trans-sialidase. After performing sequence clustering analysis with all TcS complete genes, we identified four additional groups, demonstrating that the TcS family is even more heterogeneous than previously thought. Interestingly, members of distinct TcS groups show distinctive patterns of chromosome localization. Members of the TcSgroupII, which harbor proteins involved in host cell attachment/invasion, are preferentially located in subtelomeric regions, whereas members of the largest and new TcSgroupV have internal chromosomal locations. Real-time RT-PCR confirms the expression of genes derived from new groups and shows that the pattern of expression is not similar within and between groups. We also performed B-cell epitope prediction on the family and constructed a TcS specific peptide array, which was screened with sera from T. cruzi-infected mice. We demonstrated that all seven groups represented in the array are antigenic. A highly reactive peptide occurs in sixty TcS proteins including members of two new groups and may contribute to the known cross-reactivity of T. cruzi epitopes during infection. Taken together, our results contribute to a better understanding of the real complexity of the TcS family and open new avenues for investigating novel roles of this family during T. cruzi infection.

Localization and developmental regulation of a dispersed gene family 1 protein in Trypanosoma cruzi

The dispersed gene family 1 (DGF-1) is the fifth largest gene family in the Trypanosoma cruzi genome, with over 500 members (11). Many of the predicted DGF-1 protein products have several transmembrane domains and N-glycosylation and phosphorylation sites and were thought to localize in the plasma membrane. Here, we report that affinity-purified antibodies against a region of one of these proteins (DGF-1.2) localized it intracellularly in different stages of the parasite. DGF-1.2 is more abundant in the amastigote stage than in trypomastigotes and epimastigotes, as detected by immunofluorescence and Western blot analyses. The protein changed localization during intracellular or extracellular differentiation from the trypomastigote to the amastigote stage, where it finally localized to small bodies in close contact with the inner side of the amastigote plasma membrane. DGF-1.2 did not colocalize with markers of other subcellular organelles, such as acidocalcisomes, glycosomes, reservosomes, lipid droplets, or endocytic vesicles. During extracellular differentiation, the protein was detected in the culture medium from 0 to 22 h, peaking at 14 h. The presence of DGF-1.2 in the differentiation culture medium was confirmed by mass spectrometry analysis. Finally, when epimastigotes were subjected to starvation, there was a decrease in the labeling of the cells and, in Western blots, the appearance of bands of lower molecular mass, suggesting its cleavage. These results represent the first report of direct immunodetection and developmental expression and secretion of a DGF-1 protein.

Machine learning techniques for the automated classification of adhesin-like proteins in the human protozoan parasite Trypanosoma cruzi

This paper reports on the evaluation of different machine learning techniques for the automated classification of coding gene sequences obtained from several organisms in terms of their functional role as adhesins. Diverse, biologically-meaningful, sequence-based features were extracted from the sequences and used as inputs to the in silico prediction models. Another contribution of this work is the generation of potentially novel and testable predictions about the surface protein DGF-1 family in Trypanosoma cruzi. Finally, these techniques are potentially useful for the automated annotation of known adhesin-like proteins from the trans-sialidase surface protein family in T. cruzi, the etiological agent of Chagas disease.

Molecular analysis of surface glycoprotein multigene family TrGP expressed on the plasma membrane of Trypanosoma rangeli epimastigotes forms

Trypanosoma rangeli, a non-pathogenic hemoflagelate that in Central and South America infects humans, shares with Trypanosoma cruzi reservoirs and triatomine vectors, as well as geographical distribution. Recently, we have described in T. rangeli a truncated gene copy belonging to the group II of the trans-sialidase superfamily (TrGP). This superfamily, collectively known in T. cruzi as gp85/TS, includes members that are involved in host cell invasion and infectivity. To confirm the presence of this superfamily in the genome of T. rangeli and obtain a better knowledge of its characteristics, we designed a PCR and RT-PCR cloning strategy to allow sequence analysis of both genomic and transcribed copies. We identified two full-length copies of TrGP, some pseudogenes, and N- and C-terminal sequences of several genes. We also analyzed the expression and cellular localization of these proteins in epimastigote forms of a Venezuelan T. rangeli isolate using polyclonal antibodies made against a recombinant peptide from the N-terminal region of a TrGP member. We confirmed that TrGP is a multigenic family that shares many features with T. cruzi gp85/TS, including the telomeric location of some of its members, and by immunofluorescence analysis that its location is at the surface of the parasite.

Haplotype identification within Trypanosoma cruzi I in Colombian isolates from several reservoirs, vectors and humans

Genetic variability in the Trypanosoma cruzi I group has recently been revealed in Colombian isolates from humans, reservoirs and vectors. Genomic rearrangements and the polymorphic regions in taxonomic markers, such as the miniexon gene, have led to the development of molecular tools to identify phylogenetic haplotypes in T. cruzi isolates. From genetic polymorphisms found in T. cruzi I isolates, they have been classified into four haplotypes according to their epidemiologic transmission cycles. Haplotype Ia is associated with domestic isolates, from Rhodnius prolixus; haplotype Ib, with the domestic and peridomestic cycle, mainly associated with Triatoma dimidiata; haplotype Ic is a poorly characterized group, which has been associated with the peridomestic cycle; and haplotype Id has been related to the sylvatic cycle. In order to demonstrate that the circulating T. cruzi I isolates in Colombia can be classified in the four proposed haplotypes, specific primers were designed on polymorphic regions of the miniexon gene's intergenic sequences. This set of primers allowed the molecular characterization of 33 Colombian isolates, classifying them into three of the four proposed haplotypes (Ia, Ib and Id). Results obtained from maximum parsimony and maximum-likelihood-based phylogenetic analyses correlated with the molecular classification of the isolates and their transmission cycles. This study brings insights into the Chagas disease epidemiology and the parasite's transmission dynamics.

Homology, paralogy and function of DGF-1, a highly dispersed Trypanosoma cruzi specific gene family and its implications for information entropy of its encoded proteins

Surface adhesion proteins are essential for Trypanosoma cruzi invasion of mammalian cells. Here we show that Dispersed Gene Family-1 (DGF-1) members, previously identified as nuclear repeated sequences present in several chromosomes and comprising the third largest T. cruzi specific gene family, have conserved adhesin motifs including four segments with significant similarity to human beta 7 integrin. Flow cytometry and biotinylation assays with anti-DGF-1 antibodies indicated that, as expected, DGF-1 members are expressed on the trypomastigote surface. The DGF-1 genealogy, inferred using T. cruzi Genome Project data and network phylogeny algorithms, suggests that this gene family is separated in at least three groups with differential distribution of functional domains. To identify which members of this gene family are expressed we used a combined approach of RT-PCR and codon usage profiles, showing that expressed members have a very biased codon usage favoring GC, whereas non-expressed members have a homogeneous distribution. Shannon information entropy was used to measure sequence variability and revealed four major high entropy segments in the extracellular domain of DGF-1 overlapping with important putative functional modules of the predicted proteins. Testing for natural selection, however, indicated that these high entropy segments were not under positive selection, which contradicts the notion that positive selection is the cause of high variability in specific domains of a protein relative to other less variable regions in the same molecule. We conjectured that members of the DGF-1 family might be associated with the ability of T. cruzi to bind extracellular matrix proteins, such as fibronectin and laminin, and speculated on mechanisms that would be generating the localized diversity in these molecules in the absence of selection.

An exploration of the genetic robustness landscape of surface protein families in the human protozoan parasite Trypanosoma cruzi

The ability of genes to be robust to mutations at the codon level has been suggested as a key factor for understanding adaptation features. It has been proposed that genes relevant to host-parasite interactions will tend to exhibit high volatility or "antirobust" patterns, which may be related to the capacity of the parasite to evade the host immune system. We compared two superfamilies of surface proteins, trans-sialidase (TS)-like proteins and putative surface protein dispersed gene family-1 (DGF-1), in the parasite Trypanosoma cruzi in terms of a measure of gene volatility. We proposed alternative codon robustness indicators based on cross entropy and impurity of amino acids encoded by point-mutations, which were compared to a volatility estimator previously published. This allowed us to present a more detailed description of the differences between families. A significant difference was observed in terms of these scores, with the TS-MVar1 and the DGF-1 families showing the highest and lowest gene volatility values respectively. The cross entropy and impurity estimators suggest that the MVar1 levels of volatility are linearly correlated with their capacity to generate diverse sets of amino acids as a consequence of potential mutations. This study indicates the feasibility of applying different measures of genetic robustness to detect variations between potential drug targets at the protein level.