Ancestral aneuploidy and stable chromosomal duplication resulting in differential genome structure and gene expression control in trypanosomatid parasites.

Aneuploidy is widely observed in both unicellular and multicellular eukaryotes, usually associated with adaptation to stress conditions. Chromosomal duplication stability is a tradeoff between the fitness cost of having unbalanced gene copies and the potential fitness gained from increased dosage of specific advantageous genes. Trypanosomatids, a family of protozoans that include species that cause neglected tropical diseases, are a relevant group to study aneuploidies. Their life cycle has several stressors that could select for different patterns of chromosomal duplications and/or losses, and their nearly universal use of polycistronic transcription increases their reliance on gene expansion/contraction, as well as post-transcriptional control as mechanisms for gene expression regulation. By evaluating the data from 866 isolates covering seven trypanosomatid genera, we have revealed that aneuploidy tolerance is an ancestral characteristic of trypanosomatids but has a reduced occurrence in a specific monophyletic clade that has undergone large genomic reorganization and chromosomal fusions. We have also identified an ancient chromosomal duplication that was maintained across these parasite's speciation, named collectively as the trypanosomatid ancestral supernumerary chromosome (TASC). TASC has most genes in the same coding strand, is expressed as a disomic chromosome (even having four copies), and has increased potential for functional variation, but it purges highly deleterious mutations more efficiently than other chromosomes. The evidence of stringent control over gene expression in this chromosome suggests that these parasites have adapted to mitigate the fitness cost associated with this ancient chromosomal duplication.

dgfr: an R package to assess sequence diversity of gene families.

BACKGROUND: Gene families are groups of homologous genes that often have similar biological functions. These families are formed by gene duplication events throughout evolution, resulting in multiple copies of an ancestral gene. Over time, these copies can acquire mutations and structural variations, resulting in members that may vary in size, motif ordering and sequence. Multigene families have been described in a broad range of organisms, from single-celled bacteria to complex multicellular organisms, and have been linked to an array of phenomena, such as host-pathogen interactions, immune evasion and embryonic development. Despite the importance of gene families, few approaches have been developed for estimating and graphically visualizing their diversity patterns and expression profiles in genome-wide studies. RESULTS: Here, we introduce an R package named dgfr, which estimates and enables the visualization of sequence divergence within gene families, as well as the visualization of secondary data such as gene expression. The package takes as input a multi-fasta file containing the coding sequences (CDS) or amino acid sequences from a multigene family, performs a pairwise alignment among all sequences, and estimates their distance, which is subjected to dimension reduction, optimal cluster determination, and gene assignment to each cluster. The result is a dataset that allows for the visualization of sequence divergence and expression within the gene family, an approximation of the number of clusters present in the family. CONCLUSIONS: dgfr provides a way to estimate and study the diversity of gene families, as well as visualize the dispersion and secondary profile of the sequences. The dgfr package is available at https://github.com/lailaviana/dgfr under the GPL-3 license.

Leishmania amazonensis from distinct clinical forms/hosts has polymorphisms in Lipophosphoglycans, displays variations in immunomodulatory properties and, susceptibility to antileishmanial drugs.

Lipophosphoglycan (LPG), the major Leishmania glycoconjugate, induces pro-inflammatory/immunosuppressive innate immune responses. Here, we evaluated functional/biochemical LPG properties from six Leishmania amazonensis strains from different hosts/clinical forms. LPGs from three strains (GV02, BA276, and LV79) had higher pro-inflammatory profiles for most of the mediators, including tumor necrosis factor alpha and interleukin 6. For this reason, glycoconjugates from all strains were biochemically characterized and had polymorphisms in their repeat units. They consisted of three types: type I, repeat units devoid of side chains; type II, containing galactosylated side chains; and type III, containing glucosylated side chains. No relationship was observed between LPG type and the pro-inflammatory properties. Finally, to evaluate the susceptibility against antileishmanial agents, two strains with high (GV02, BA276) and one with low (BA336) pro-inflammatory activity were selected for chemotherapeutic tests in THP-1 cells. All analyzed strains were susceptible to amphotericin B (AmB) but displayed various responses against miltefosine (MIL) and glucantime (GLU). The GV02 strain (canine visceral leishmaniasis) had the highest IC50 for MIL (3.34 µM), whereas diffuse leishmaniasis strains (BA276 and BA336) had a higher IC50 for GLU (6.87-12.19 mM). The highest IC50 against MIL shown by the GV02 strain has an impact on clinical management. Miltefosine is the only drug approved for dog treatment in Brazil. Further studies into drug susceptibility of L. amazonensis strains are warranted, especially in areas where dog infection by this species overlaps with those caused by Leishmania infantum.

Trypanosoma cruzi genetic diversity: impact on transmission cycles and Chagas disease.

Trypanosoma cruzi, the agent of Chagas disease (ChD), exhibits remarkable biological and genetic diversity, along with eco-epidemiological complexity. In order to facilitate communication among researchers aiming at the characterisation of biological and epidemiological aspects of T. cruzi, parasite isolates and strains were partitioned into seven discrete typing units (DTUs), TcI-TcVI and TcBat, identifiable by reproducible genotyping protocols. Here we present the potential origin of the genetic diversity of T. cruzi and summarise knowledge about eco-epidemiological associations of DTUs with mammalian reservoirs and vectors. Circumstantial evidence of a connection between T. cruzi genotype and ChD manifestations is also discussed emphasising the role of the host's immune response in clinical ChD progression. We describe genomic aspects of DTUs focusing on polymorphisms in multigene families encoding surface antigens that play essential functions for parasite survival both in the insect vector and the mammalian host. Such antigens most probably contributed to the parasite success in establishing infections in different hosts and exploring several niches. Gaps in the current knowledge and challenges for future research are pointed out.

BepFAMN: A Method for Linear B-Cell Epitope Predictions Based on Fuzzy-ARTMAP Artificial Neural Network.

The public health system is extremely dependent on the use of vaccines to immunize the population from a series of infectious and dangerous diseases, preventing the system from collapsing and millions of people dying every year. However, to develop these vaccines and effectively monitor these diseases, it is necessary to use accurate diagnostic methods capable of identifying highly immunogenic regions within a given pathogenic protein. Existing experimental methods are expensive, time-consuming, and require arduous laboratory work, as they require the screening of a large number of potential candidate epitopes, making the methods extremely laborious, especially for application to larger microorganisms. In the last decades, researchers have developed in silico prediction methods, based on machine learning, to identify these markers, to drastically reduce the list of potential candidate epitopes for experimental tests, and, consequently, to reduce the laborious task associated with their mapping. Despite these efforts, the tools and methods still have low accuracy, slow diagnosis, and offline training. Thus, we develop a method to predict B-cell linear epitopes which are based on a Fuzzy-ARTMAP neural network architecture, called BepFAMN (B Epitope Prediction Fuzzy ARTMAP Artificial Neural Network). This was trained using a linear averaging scheme on 15 properties that include an amino acid ratio scale and a set of 14 physicochemical scales. The database used was obtained from the IEDB website, from which the amino acid sequences with the annotations of their positive and negative epitopes were taken. To train and validate the knowledge models, five-fold cross-validation and competition techniques were used. The BepiPred-2.0 database, an independent database, was used for the tests. In our experiment, the validation dataset reached sensitivity = 91.50%, specificity = 91.49%, accuracy = 91.49%, MCC = 0.83, and an area under the curve (AUC) ROC of approximately 0.9289. The result in the testing dataset achieves a significant improvement, with sensitivity = 81.87%, specificity = 74.75%, accuracy = 78.27%, MCC = 0.56, and AOC = 0.7831. These achieved values demonstrate that BepFAMN outperforms all other linear B-cell epitope prediction tools currently used. In addition, the architecture provides mechanisms for online training, which allow the user to find a new B-cell linear epitope, and to improve the model without need to re-train itself with the whole dataset. This fact contributes to a considerable reduction in the number of potential linear epitopes to be experimentally validated, reducing laboratory time and accelerating the development of diagnostic tests, vaccines, and immunotherapeutic approaches.

Genomics and functional genomics in Leishmania and Trypanosoma cruzi: statuses, challenges and perspectives.

The availability of Trypanosomatid genomic data in public databases has opened myriad experimental possibilities that have contributed to a more comprehensive understanding of the biology of these parasites and their interactions with hosts. In this review, after brief remarks on the history of the Trypanosoma cruzi and Leishmania genome initiatives, we present an overview of the relevant contributions of genomics, transcriptomics and functional genomics, discussing the primary obstacles, challenges, relevant achievements and future perspectives of these technologies.

Innate immune recognition of an AT-rich stem-loop DNA motif in the Plasmodium falciparum genome.

Although Toll-like receptor 9 (TLR9) has been implicated in cytokine and type I interferon (IFN) production during malaria in humans and mice, the high AT content of the Plasmodium falciparum genome prompted us to examine the possibility that malarial DNA triggered TLR9-independent pathways. Over 6000 ATTTTTAC ("AT-rich") motifs are present in the genome of P. falciparum, which we show here potently induce type I IFNs. Parasite DNA, parasitized erythrocytes and oligonucleotides containing the AT-rich motif induce type I IFNs via a pathway that did not involve the previously described sensors TLR9, DAI, RNA polymerase-III or IFI16/p204. Rather, AT-rich DNA sensing involved an unknown receptor that coupled to the STING, TBK1 and IRF3-IRF7 signaling pathway. Mice lacking IRF3, IRF7, the kinase TBK1 or the type I IFN receptor were resistant to otherwise lethal cerebral malaria. Collectively, these observations implicate AT-rich DNA sensing via STING, TBK1 and IRF3-IRF7 in P. falciparum malaria.

Conditional genome engineering reveals canonical and divergent roles for the Hus1 component of the 9-1-1 complex in the maintenance of the plastic genome of Leishmania.

Leishmania species are protozoan parasites whose remarkably plastic genome limits the establishment of effective genetic manipulation and leishmaniasis treatment. The strategies used by Leishmania to maintain its genome while allowing variability are not fully understood. Here, we used DiCre-mediated conditional gene deletion to show that HUS1, a component of the 9-1-1 (RAD9-RAD1-HUS1) complex, is essential and is required for a G2/M checkpoint. By analyzing genome-wide instability in HUS1 ablated cells, HUS1 is shown to have a conserved role, by which it preserves genome stability and also a divergent role, by which it promotes genome variability. These roles of HUS1 are related to distinct patterns of formation and resolution of single-stranded DNA and γH2A, throughout the cell cycle. Our findings suggest that Leishmania 9-1-1 subunits have evolved to co-opt canonical genomic maintenance and genomic variation functions. Hence, this study reveals a pivotal function of HUS1 in balancing genome stability and transmission in Leishmania. These findings may be relevant to understanding the evolution of genome maintenance and plasticity in other pathogens and eukaryotes.

Comparative transcriptome profiling of virulent and non-virulent Trypanosoma cruzi underlines the role of surface proteins during infection.

Trypanosoma cruzi, the protozoan that causes Chagas disease, has a complex life cycle involving several morphologically and biochemically distinct stages that establish intricate interactions with various insect and mammalian hosts. It has also a heterogeneous population structure comprising strains with distinct properties such as virulence, sensitivity to drugs, antigenic profile and tissue tropism. We present a comparative transcriptome analysis of two cloned T. cruzi strains that display contrasting virulence phenotypes in animal models of infection: CL Brener is a virulent clone and CL-14 is a clone that is neither infective nor pathogenic in in vivo models of infection. Gene expression analysis of trypomastigotes and intracellular amastigotes harvested at 60 and 96 hours post-infection (hpi) of human fibroblasts revealed large differences that reflect the parasite's adaptation to distinct environments during the infection of mammalian cells, including changes in energy sources, oxidative stress responses, cell cycle control and cell surface components. While extensive transcriptome remodeling was observed when trypomastigotes of both strains were compared to 60 hpi amastigotes, differences in gene expression were much less pronounced when 96 hpi amastigotes and trypomastigotes of CL Brener were compared. In contrast, the differentiation of the avirulent CL-14 from 96 hpi amastigotes to extracellular trypomastigotes was associated with considerable changes in gene expression, particularly in gene families encoding surface proteins such as trans-sialidases, mucins and the mucin associated surface proteins (MASPs). Thus, our comparative transcriptome analysis indicates that the avirulent phenotype of CL-14 may be due, at least in part, to a reduced or delayed expression of genes encoding surface proteins that are associated with the transition of amastigotes to trypomastigotes, an essential step in the establishment of the infection in the mammalian host. Confirming the role of members of the trans-sialidase family of surface proteins for parasite differentiation, transfected CL-14 constitutively expressing a trans-sialidase gene displayed faster kinetics of trypomastigote release in the supernatant of infected cells compared to wild type CL-14.

Whole genome sequencing of Trypanosoma cruzi field isolates reveals extensive genomic variability and complex aneuploidy patterns within TcII DTU.

BACKGROUND: Trypanosoma cruzi, the etiologic agent of Chagas disease, is currently divided into six discrete typing units (DTUs), named TcI-TcVI. TcII is among the major DTUs enrolled in human infections in South America southern cone, where it is associated with severe cardiac and digestive symptoms. Despite the importance of TcII in Chagas disease epidemiology and pathology, so far, no genome-wide comparisons of the mitochondrial and nuclear genomes of TcII field isolates have been performed to track the variability and evolution of this DTU in endemic regions. RESULTS: In the present work, we have sequenced and compared the whole nuclear and mitochondrial genomes of seven TcII strains isolated from chagasic patients from the central and northeastern regions of Minas Gerais, Brazil, revealing an extensive genetic variability within this DTU. A comparison of the phylogeny based on the nuclear or mitochondrial genomes revealed that the majority of branches were shared by both sequences. The subtle divergences in the branches are probably consequence of mitochondrial introgression events between TcII strains. Two T. cruzi strains isolated from patients living in the central region of Minas Gerais, S15 and S162a, were clustered in the nuclear and mitochondrial phylogeny analysis. These two strains were isolated from the other five by the Espinhaço Mountains, a geographic barrier that could have restricted the traffic of insect vectors during T. cruzi evolution in the Minas Gerais state. Finally, the presence of aneuploidies was evaluated, revealing that all seven TcII strains have a different pattern of chromosomal duplication/loss. CONCLUSIONS: Analysis of genomic variability and aneuploidies suggests that there is significant genomic variability within Minas Gerais TcII strains, which could be exploited by the parasite to allow rapid selection of favorable phenotypes. Also, the aneuploidy patterns vary among T. cruzi strains and does not correlate with the nuclear phylogeny, suggesting that chromosomal duplication/loss are recent and frequent events in the parasite evolution.