Comparative ribosome profiling reveals extensive translational complexity in different Trypanosoma brucei life cycle stages

Transcriptomic analysis of N-terminal mutated Trypanosoma cruzi UBP1 knockdown underlines the importance of this RNA-binding protein in parasite development

BACKGROUND

During its life cycle, the human pathogen Trypanosoma cruzi must quickly adapt to different environments, in which the variation in the gene expression of the regulatory U-rich RNA-binding protein 1 (TcUBP1) plays a crucial role. We have previously demonstrated that the overexpression of TcUBP1 in insect-dwelling epimastigotes orchestrates an RNA regulon to promote differentiation to infective forms.

METHODS

In an attempt to generate TcUBP1 knockout parasites by using CRISPR-Cas9 technology, in the present study, we obtained a variant transcript that encodes a protein with 95% overall identity and a modified N-terminal sequence. The expression of this mutant protein, named TcUBP1mut, was notably reduced compared to that of the endogenous form found in normal cells. TcUBP1mut-knockdown epimastigotes exhibited normal growth and differentiation into infective metacyclic trypomastigotes and were capable of infecting mammalian cells.

RESULTS

We analyzed the RNA-Seq expression profiles of these parasites and identified 276 up- and 426 downregulated genes with respect to the wildtype control sample. RNA-Seq comparison across distinct developmental stages revealed that the transcriptomic profile of these TcUBP1mut-knockdown epimastigotes significantly differs not only from that of epimastigotes in the stationary phase but also from the gene expression landscape characteristic of infective forms. This is both contrary to and consistent with the results of our recent study involving TcUBP1-overexpressing cells.

CONCLUSION

Together, our findings demonstrate that the genes exhibiting opposite changes under overexpression and knockdown conditions unveil key mRNA targets regulated by TcUBP1. These mostly encompass transcripts that encode for trypomastigote-specific surface glycoproteins and ribosomal proteins, supporting a role for TcUBP1 in determining the molecular characteristics of the infective stage.

Translational control by Trypanosoma brucei DRBD18 contributes to the maintenance of the procyclic state

Trypanosoma brucei occupies distinct niches throughout its life cycle, within both the mammalian and tsetse fly hosts. The immunological and biochemical complexity and variability of each of these environments require a reshaping of the protein landscape of the parasite both to evade surveillance and face changing metabolic demands. In kinetoplastid protozoa, including T. brucei, posttranscriptional control mechanisms are the primary means of gene regulation, and these are often mediated by RNA-binding proteins. DRBD18 is a T. brucei RNA-binding protein that reportedly interacts with ribosomal proteins and translation factors. Here, we tested a role for DRBD18 in translational control. We validate the DRBD18 interaction with translating ribosomes and the translation initiation factor, eIF3a. We further show that DRBD18 depletion by RNA interference leads to altered polysomal profiles with a specific depletion of heavy polysomes. Ribosome profiling analysis reveals that 101 transcripts change in translational efficiency (TE) upon DRBD18 depletion: 41 exhibit decreased TE and 60 exhibit increased TE. A further 66 transcripts are buffered, that is, changes in transcript abundance are compensated by changes in TE such that the total translational output is expected not to change. In DRBD18-depleted cells, a set of transcripts that codes for procyclic form-specific proteins is translationally repressed while, conversely, transcripts that code for bloodstream form- and metacyclic form-specific proteins are translationally enhanced. RNA immunoprecipitation/qRT-PCR indicates that DRBD18 associates with members of both repressed and enhanced cohorts. These data suggest that DRBD18 contributes to the maintenance of the procyclic state through both positive and negative translational control of specific mRNAs.

Synthesis and characterization of I-BET151 derivatives for use in identifying protein targets in the African trypanosome

Trypanosoma brucei, the causative agent of Human African Trypanosomiasis (HAT) and animal trypanosomiases, cycles between a bloodstream form in mammals and a procyclic form in the gut of its insect vector. We previously discovered that the human bromodomain inhibitor I-BET151 causes transcriptome changes that resemble the transition from the bloodstream to the procyclic form. In particular, I-BET151 induces replacement of variant surface glycoprotein (VSG) with procyclin protein. While modest binding of I-BET151 to TbBdf2 and TbBdf3 has been demonstrated, it is unknown whether I-BET151 binds to other identified T. brucei bromodomain proteins and/or other targets. To identify target(s) in T. brucei, we have synthesized I-BET151 derivatives maintaining the key pharmacophoric elements with functionality useful for chemoproteomic approaches. We identified compounds that are potent in inducing expression of procyclin, delineating a strategy towards the design of drugs against HAT and other trypanosomiases. Furthermore, these derivatives represent useful chemical probes to elucidate the molecular mechanism underlying I-BET151-induced differentiation.

Refinement of Leishmania donovani Genome Annotations in the Light of Ribosome-Protected mRNAs Fragments (Ribo-Seq Data)

Advances in next-generation sequencing methodologies have facilitated the assembly of an ever-increasing number of genomes. Gene annotations are typically conducted via specialized software, but the most accurate results require additional manual curation that incorporates insights derived from functional and bioinformatic analyses (e.g., transcriptomics, proteomics, and phylogenetics). In this study, we improved the annotation of the Leishmania donovani (strain HU3) genome using publicly available data from the deep sequencing of ribosome-protected mRNA fragments (Ribo-Seq). As a result of this analysis, we uncovered 70 previously non-annotated protein-coding genes and improved the annotation of around 600 genes. Additionally, we present evidence for small upstream open reading frames (uORFs) in a significant number of transcripts, indicating their potential role in the translational regulation of gene expression. The bioinformatics pipelines developed for these analyses can be used to improve the genome annotations of other organisms for which Ribo-Seq data are available. The improvements provided by these studies will bring us closer to the ultimate goal of a complete and accurately annotated L. donovani genome and will enhance future transcriptomics, proteomics, and genetics studies.

Subcellular protein localisation of Trypanosoma brucei bloodstream form-upregulated proteins maps stage-specific adaptations

Background: Genome-wide subcellular protein localisation in Trypanosoma brucei, through our TrypTag project, has comprehensively dissected the molecular organisation of this important pathogen. Powerful as this resource is , T. brucei has multiple developmental forms and we previously only analysed the procyclic form. This is an insect life cycle stage, leaving the mammalian bloodstream form unanalysed. The expectation is that between life stages protein localisation would not change dramatically (completely unchanged or shifting to analogous stage-specific structures). However, this has not been specifically tested. Similarly, which organelles tend to contain proteins with stage-specific expression can be predicted from known stage specific adaptations but has not been comprehensively tested. Methods: We used endogenous tagging with mNG to determine the sub-cellular localisation of the majority of proteins encoded by transcripts significantly upregulated in the bloodstream form, and performed comparison to the existing localisation data in procyclic forms. Results: We have confirmed the localisation of known stage-specific proteins and identified the localisation of novel stage-specific proteins. This gave a map of which organelles tend to contain stage specific proteins: the mitochondrion for the procyclic form, and the endoplasmic reticulum, endocytic system and cell surface in the bloodstream form. Conclusions: This represents the first genome-wide map of life cycle stage-specific adaptation of organelle molecular machinery in T. brucei.

Profiling the bloodstream form and procyclic form Trypanosoma brucei cell cycle using single-cell transcriptomics

African trypanosomes proliferate as bloodstream forms (BSFs) and procyclic forms in the mammal and tsetse fly midgut, respectively. This allows them to colonise the host environment upon infection and ensure life cycle progression. Yet, understanding of the mechanisms that regulate and drive the cell replication cycle of these forms is limited. Using single-cell transcriptomics on unsynchronised cell populations, we have obtained high resolution cell cycle regulated (CCR) transcriptomes of both procyclic and slender BSF Trypanosoma brucei without prior cell sorting or synchronisation. Additionally, we describe an efficient freeze-thawing protocol that allows single-cell transcriptomic analysis of cryopreserved T. brucei. Computational reconstruction of the cell cycle using periodic pseudotime inference allowed the dynamic expression patterns of cycling genes to be profiled for both life cycle forms. Comparative analyses identify a core cycling transcriptome highly conserved between forms, as well as several genes where transcript levels dynamics are form specific. Comparing transcript expression patterns with protein abundance revealed that the majority of genes with periodic cycling transcript and protein levels exhibit a relative delay between peak transcript and protein expression. This work reveals novel detail of the CCR transcriptomes of both forms, which are available for further interrogation via an interactive webtool.

Translational reprogramming as a driver of antimony-drug resistance in Leishmania

Leishmania is a unicellular protozoan that has a limited transcriptional control and mostly uses post-transcriptional regulation of gene expression, although the molecular mechanisms of the process are still poorly understood. Treatments of leishmaniasis, pathologies associated with Leishmania infections, are limited due to drug resistance. Here, we report dramatic differences in mRNA translation in antimony drug-resistant and sensitive strains at the full translatome level. The major differences (2431 differentially translated transcripts) were demonstrated in the absence of the drug pressure supporting that complex preemptive adaptations are needed to efficiently compensate for the loss of biological fitness once they are exposed to the antimony. In contrast, drug-resistant parasites exposed to antimony activated a highly selective translation of only 156 transcripts. This selective mRNA translation is associated with surface protein rearrangement, optimized energy metabolism, amastins upregulation, and improved antioxidant response. We propose a novel model that establishes translational control as a major driver of antimony-resistant phenotypes in Leishmania.

Ribosome Specialization in Protozoa Parasites

Ribosomes, in general, are viewed as constitutive macromolecular machines where protein synthesis takes place; however, this view has been recently challenged, supporting the hypothesis of ribosome specialization and opening a completely new field of research. Recent studies have demonstrated that ribosomes are heterogenous in their nature and can provide another layer of gene expression control by regulating translation. Heterogeneities in ribosomal RNA and ribosomal proteins that compose them favor the selective translation of different sub-pools of mRNAs and functional specialization. In recent years, the heterogeneity and specialization of ribosomes have been widely reported in different eukaryotic study models; however, few reports on this topic have been made on protozoa and even less on protozoa parasites of medical importance. This review analyzes heterogeneities of ribosomes in protozoa parasites highlighting the specialization in their functions and their importance in parasitism, in the transition between stages in their life cycle, in the change of host and in response to environmental conditions.

RNA-seq reveals that overexpression of TcUBP1 switches the gene expression pattern towards that of the infective form of Trypanosoma cruzi

Trypanosomes regulate gene expression mainly by using post-transcriptional mechanisms. Key factors responsible for carrying out this regulation are RNA-binding proteins (RBPs), affecting subcellular localization, translation, and/or transcript stability. Trypanosoma cruzi U-rich RBP 1 (TcUBP1) is a small protein that modulates the expression of several surface glycoproteins of the trypomastigote infective stage of the parasite. Its mRNA targets are known but the impact of its overexpression at the transcriptome level in the insect-dwelling epimastigote cells has not yet been investigated. Thus, in the present study, by using a tetracycline-inducible system, we generated a population of TcUBP1-overexpressing parasites and analyzed its effect by RNA-seq methodology. This allowed us to identify 793 up- and 371 down-regulated genes with respect to the wild-type control sample. Among the up-regulated genes, it was possible to identify members coding for the TcS superfamily, MASP, MUCI/II, and protein kinases, whereas among the down-regulated transcripts, we found mainly genes coding for ribosomal, mitochondrial, and synthetic pathway proteins. RNA-seq comparison with two previously published datasets revealed that the expression profile of this TcUBP1-overexpressing replicative epimastigote form resembles the transition to the infective metacyclic trypomastigote stage. We identified novel cis-regulatory elements in the 3'-untranslated region of the affected transcripts and confirmed that UBP1m -a signature TcUBP1 binding element previously characterized in our lab- is enriched in the list of stabilized genes. We can conclude that the overall effect of TcUBP1 overexpression on the epimastigote transcriptome is mainly the stabilization of mRNAs coding for proteins that are important for parasite infection.

Cell-to-flagellum attachment and surface architecture in kinetoplastids

A key morphological feature of kinetoplastid parasites is the position and length of flagellum attachment to the cell body. This lateral attachment is mediated by the flagellum attachment zone (FAZ), a large complex cytoskeletal structure, which is essential for parasite morphogenesis and pathogenicity. Despite the complexity of the FAZ only two transmembrane proteins, FLA1 and FLA1BP, are known to interact and connect the flagellum to the cell body. Across the different kinetoplastid species, each only has a single FLA/FLABP pair, except in Trypanosoma brucei and Trypanosoma congolense where there has been an expansion of these genes. Here, we focus on the selection pressure behind the evolution of the FLA/FLABP proteins and the likely impact this will have on host-parasite interactions.

Translational control by Trypanosoma brucei DRBD18 contributes to the maintenance of the procyclic state

Trypanosoma brucei occupies distinct niches throughout its life cycle, within both the mammalian and tsetse fly hosts. The immunological and biochemical complexity and variability of each of these environments require a reshaping of the protein landscape of the parasite both to evade surveillance and face changing metabolic demands. Whereas most well-studied organisms rely on transcriptional control as the main regulator of gene expression, post-transcriptional control mechanisms are particularly important in T. brucei , and these are often mediated by RNA-binding proteins. DRBD18 is a T. brucei RNA-binding protein that interacts with ribosomal proteins and translation factors. Here, we tested a role for DRBD18 in translational control. We show that DRBD18 depletion by RNA interference leads to altered polysomal profiles with a specific depletion of heavy polysomes. Ribosome profiling analysis reveals that 101 transcripts change in translational efficiency (TE) upon DRBD18 depletion: 41 exhibit decreased TE and 60 exhibit increased TE. A further 66 transcripts are buffered, i.e . changes in transcript abundance are compensated by changes in TE such that the total translational output is expected not to change. Proteomic analysis validates these data. In DRBD18-depleted cells, a cohort of transcripts that codes for procyclic form-specific proteins is translationally repressed while, conversely, transcripts that code for bloodstream form- and metacyclic form-specific proteins are translationally enhanced. These data suggest that DRBD18 contributes to the maintenance of the procyclic state through both positive and negative translational control of specific mRNAs.

Subcellular protein localisation of Trypanosoma brucei bloodstream form-upregulated proteins maps stage-specific adaptations

Background: Genome-wide subcellular protein localisation in Trypanosoma brucei, through our TrypTag project, has comprehensively dissected the molecular organisation of this important pathogen. Powerful as this resource is, T. brucei has multiple developmental forms and we previously only analysed the procyclic form. This is an insect life cycle stage, leaving the mammalian bloodstream form unanalysed. The expectation is that between life stages protein localisation would not change dramatically (completely unchanged or shifting to analogous stage-specific structures). However, this has not been specifically tested. Similarly, which organelles tend to contain proteins with stage-specific expression can be predicted from known stage specific adaptations but has not been comprehensively tested. Methods: We used endogenous tagging with mNG to determine the sub-cellular localisation of the majority of proteins encoded by transcripts significantly upregulated in the bloodstream form, and performed comparison to the existing localisation data in procyclic forms. Results: We have confirmed the localisation of known and identified the localisation of novel stage-specific proteins. This gave a map of which organelles tend to contain stage specific proteins: the mitochondrion for the procyclic form, and the endoplasmic reticulum, endocytic system and cell surface in the bloodstream form. Conclusions: This represents the first genome-wide map of life cycle stage-specific adaptation of organelle molecular machinery in T. brucei.

Trypanosoma brucei Acyl-Protein Thioesterase-like (TbAPT-L) Is a Lipase with Esterase Activity for Short and Medium-Chain Fatty Acids but Has No Depalmitoylation Activity

Dynamic post-translational modifications allow the rapid, specific, and tunable regulation of protein functions in eukaryotic cells. S-acylation is the only reversible lipid modification of proteins, in which a fatty acid, usually palmitate, is covalently attached to a cysteine residue of a protein by a zDHHC palmitoyl acyltransferase enzyme. Depalmitoylation is required for acylation homeostasis and is catalyzed by an enzyme from the alpha/beta hydrolase family of proteins usually acyl-protein thioesterase (APT1). The enzyme responsible for depalmitoylation in Trypanosoma brucei parasites is currently unknown. We demonstrate depalmitoylation activity in live bloodstream and procyclic form trypanosomes sensitive to dose-dependent inhibition with the depalmitoylation inhibitor, palmostatin B. We identified a homologue of human APT1 in Trypanosoma brucei which we named TbAPT-like (TbAPT-L). Epitope-tagging of TbAPT-L at N- and C- termini indicated a cytoplasmic localization. Knockdown or over-expression of TbAPT-L in bloodstream forms led to robust changes in TbAPT-L mRNA and protein expression but had no effect on parasite growth in vitro, or cellular depalmitoylation activity. Esterase activity in cell lysates was also unchanged when TbAPT-L was modulated. Unexpectedly, recombinant TbAPT-L possesses esterase activity with specificity for short- and medium-chain fatty acid substrates, leading to the conclusion, TbAPT-L is a lipase, not a depalmitoylase.

The Trypanosoma brucei RNA-binding protein DRBD18 ensures correct mRNA trans splicing and polyadenylation patterns

The parasite Trypanosoma brucei grows as bloodstream forms in mammals, and as procyclic forms in tsetse flies. Transcription is polycistronic, all mRNAs are trans spliced, and polyadenylation sites are defined by downstream splicing signals. Expression regulation therefore depends heavily on post-transcriptional mechanisms. The RNA-binding protein DRBD18 was previously implicated in the export of some mRNAs from the nucleus in procyclic forms. It copurifies the outer ring of the nuclear pore, mRNA export factors and exon-junction-complex proteins. We show that for more than 200 mRNAs, DRBD18 depletion caused preferential accumulation of versions with shortened 3'-untranslated regions, arising from use of polyadenylation sites that were either undetectable or rarely seen in nondepleted cells. The shortened mRNAs were often, but not always, more abundant in depleted cells than the corresponding longer versions in normal cells. Their appearance was linked to the appearance of trans-spliced, polyadenylated RNAs containing only downstream 3'-untranslated region-derived sequences. Experiments with one mRNA suggested that nuclear retention alone, through depletion of MEX67, did not affect mRNA length, suggesting a specific effect of DRBD18 on processing. DRBD18-bound mRNAs were enriched in polypyrimidine tract motifs, and DRBD18 was found in both the nucleus and the cytoplasm. We therefore suggest that in the nucleus, DRBD18 might bind to polypyrimidine tracts in 3'-UTRs of mRNA precursors. Such binding might both prevent recognition of mRNA-internal polypyrimidine tracts by splicing factors, and promote export of the processed bound mRNAs to the cytosol.

Stage-specific transcription activator ESB1 regulates monoallelic antigen expression in Trypanosoma brucei

Variant surface glycoprotein (VSG) coats bloodstream form Trypanosoma brucei parasites, and monoallelic VSG expression underpins the antigenic variation necessary for pathogenicity. One of thousands of VSG genes is transcribed by RNA polymerase I in a singular nuclear structure called the expression site body (ESB), but how monoallelic VSG transcription is achieved remains unclear. Using a localization screen of 153 proteins we found one, ESB-specific protein 1 (ESB1), that localized only to the ESB and is expressed only in VSG-expressing life cycle stages. ESB1 associates with DNA near the active VSG promoter and is necessary for VSG expression, with overexpression activating inactive VSG promoters. Mechanistically, ESB1 is necessary for recruitment of a subset of ESB components, including RNA polymerase I, revealing that the ESB has separately assembled subdomains. Because many trypanosomatid parasites have divergent ESB1 orthologues yet do not undergo antigenic variation, ESB1 probably represents an important class of transcription regulators.

Gene expression study to elucidate the anti-trypanosomal activity of quinapyramine methyl sulphate (QPS)

The kinetoplastid protozoan parasite, Trypanosoma evansi causes a fatal disease condition known as Surra in equines throughout the globe. Disease condition being acute in nature, entrust a huge economic and health impact on the equine industry. Till date, quinapyramine methyl sulphate (QPS) is the first line of treatment and a panacea for the T. evansi infection in equines. Still after the >70 years of its discovery, there is no clue about the mode of action of QPS in T. evansi. The establishment of in vitro cultivation of T. evansi in HMI-9 media has provided opportunity to study the alteration in mRNA expression of parasite on exposure to the drug. With this research gap, the present study aimed to investigate the relative mRNA expression of 13 important drug target genes to elucidate the anti-trypanosomal activity of QPS against T. evansi. The IC₅₀ of QPS against a pony isolate of T. evansi was determined as 276.4 nM(147.21 ng/ mL) in the growth inhibitory assay. The in vitro cultured T. evansi population were further exposed to IC₅₀ of QPS and their relative mRNA expression was studied at 12 h, 24 h and 48 h interval.The mRNA expression of several genes such as hexokinase, trypanothione reductase, aurora kinase, oligopeptidase B and ribonucleotide reductase II were found refractory (non-significant, p > 0.1234) to the exposure of QPS. Significant up-regulation of trans-sialidase (p < 0.0001), ESAG8 (p < 0.0021), ribonucleotide reductase I (p < 0.0001), ornithine decarboxylase (p < 0.0001), topoisomerase II (p < 0.0021) and casein kinase I (p < 0.0021) were recorded after exposure with QPS. The arginine kinase 1 and calcium ATPase I showed highly significant (p < 0.0001) down-regulation in the drug kinetics. Therefore, the arginine kinase 1 and calcium ATPase I can be explored further to elucidate the trypanocidal activity of QPS. The preliminary data generated provide the potential of arginine kinase 1 and calcium ATPase I mRNA mediated pathway of trypanocidal action of QPS. Further, transcriptomics approach is required to investigate the possible mechanism of action of drugs at molecular level against the targeted organism.

A long noncoding RNA promotes parasite differentiation in African trypanosomes

The parasite Trypanosoma brucei causes African sleeping sickness that is fatal to patients if untreated. Parasite differentiation from a replicative slender form into a quiescent stumpy form promotes host survival and parasite transmission. Long noncoding RNAs (lncRNAs) are known to regulate cell differentiation in other eukaryotes. To determine whether lncRNAs are also involved in parasite differentiation, we used RNA sequencing to survey the T. brucei genome, identifying 1428 previously uncharacterized lncRNA genes. We find that grumpy lncRNA is a key regulator that promotes parasite differentiation into the quiescent stumpy form. This function is promoted by a small nucleolar RNA encoded within the grumpy lncRNA. snoGRUMPY binds to messenger RNAs of at least two stumpy regulatory genes, promoting their expression. grumpy overexpression reduces parasitemia in infected mice. Our analyses suggest that T. brucei lncRNAs modulate parasite-host interactions and provide a mechanism by which grumpy regulates cell differentiation in trypanosomes.

Open chromatin analysis in Trypanosoma cruzi life forms highlights critical differences in genomic compartments and developmental regulation at tDNA loci

BACKGROUND

Genomic organization and gene expression regulation in trypanosomes are remarkable because protein-coding genes are organized into codirectional gene clusters with unrelated functions. Moreover, there is no dedicated promoter for each gene, resulting in polycistronic gene transcription, with posttranscriptional control playing a major role. Nonetheless, these parasites harbor epigenetic modifications at critical regulatory genome features that dynamically change among parasite stages, which are not fully understood.

RESULTS

Here, we investigated the impact of chromatin changes in a scenario commanded by posttranscriptional control exploring the parasite Trypanosoma cruzi and its differentiation program using FAIRE-seq approach supported by transmission electron microscopy. We identified differences in T. cruzi genome compartments, putative transcriptional start regions, and virulence factors. In addition, we also detected a developmental chromatin regulation at tRNA loci (tDNA), which could be linked to the intense chromatin remodeling and/or the translation regulatory mechanism required for parasite differentiation. We further integrated the open chromatin profile with public transcriptomic and MNase-seq datasets. Strikingly, a positive correlation was observed between active chromatin and steady-state transcription levels.

CONCLUSION

Taken together, our results indicate that chromatin changes reflect the unusual gene expression regulation of trypanosomes and the differences among parasite developmental stages, even in the context of a lack of canonical transcriptional control of protein-coding genes.

Low Dose Gamma Irradiation of Trypanosoma evansi Parasites Identifies Molecular Changes That Occur to Repair Radiation Damage and Gene Transcripts That May Be Involved in Establishing Disease in Mice Post-Irradiation

The protozoan parasite Trypanosoma evansi is responsible for causing surra in a variety of mammalian hosts and is spread by many vectors over a wide geographical area making it an ideal target for irradiation as a tool to study the initial events that occur during infection. Parasites irradiated at the representative doses 100Gy, 140Gy, and 200Gy were used to inoculate BALB/c mice revealing that parasites irradiated at 200Gy were unable to establish disease in all mice. Cytokine analysis of mice inoculated with 200Gy of irradiated parasites showed significantly lower levels of interleukins when compared to mice inoculated with non-irradiated and 100Gy irradiated parasites. Irradiation also differentially affected the abundance of gene transcripts in a dose-dependent trend measured at 6- and 20-hours post-irradiation with 234, 325, and 484 gene transcripts affected 6 hours post-irradiation for 100Gy-, 140Gy- and 200Gy-irradiated parasites, respectively. At 20 hours post-irradiation, 422, 381, and 457 gene transcripts were affected by irradiation at 100Gy, 140Gy, and 200Gy, respectively. A gene ontology (GO) term analysis was carried out for the three representative doses at 6 hours and 20 hours post-irradiation revealing different processes occurring at 20 hours when compared to 6 hours for 100Gy irradiation. The top ten most significant processes had a negative Z score. These processes fall in significance at 140Gy and even further at 200Gy, revealing that they were least likely to occur at 200Gy, and thus may have been responsible for infection in mice by 100Gy and 140Gy irradiated parasites. When looking at 100Gy irradiated parasites 20 hours post-irradiation processes with a positive Z score, we identified genes that were involved in multiple processes and compared their fold change values at 6 hours and 20 hours. We present these genes as possibly necessary for repair from irradiation damage at 6 hours and suggestive of being involved in the establishment of disease in mice at 20 hours post-irradiation. A potential strategy using this information to develop a whole parasite vaccine is also postulated.

Environment-driven shifts in inter-individual variation and phenotypic integration within subnetworks of the mussel transcriptome and proteome

The environment can alter the magnitude of phenotypic variation among individuals, potentially influencing evolutionary trajectories. However, environmental influences on variation are complex and remain understudied. Populations in heterogeneous environments might exhibit more variation, the amount of variation could differ between benign and stressful conditions, and/or variation might manifest in different ways among stages of the gene‐to‐protein expression cascade or among physiological functions. Here, we explore these three issues by quantifying patterns of inter‐individual variation in both transcript and protein expression levels among California mussels, Mytilus californianus Conrad. Mussels were exposed to five ecologically relevant treatments that varied in the mean and interindividual heterogeneity of body temperature. To target a diverse set of physiological functions, we assessed variation within 19 expression subnetworks, including canonical stress‐response pathways and empirically derived coexpression clusters that represent a diffuse set of cellular processes. Variation in expression was particularly pronounced in the treatments with high mean and heterogeneous body temperatures. However, with few exceptions, environment‐dependent shifts of variation in the transcriptome were not reflected in the proteome. A metric of phenotypic integration provided evidence for a greater degree of constraint on relative expression levels (i.e., stronger correlation) within expression subnetworks in benign, homogeneous environments. Our results suggest that environments that are more stressful on average – and which also tend to be more heterogeneous – can relax these expression constraints and reduce phenotypic integration within biochemical subnetworks. Context‐dependent “unmasking” of functional variation may contribute to interindividual differences in physiological phenotype and performance in stressful environments.

Several different sequences are implicated in bloodstream-form-specific gene expression in Trypanosoma brucei

The parasite Trypanosoma brucei grows as bloodstream forms in mammalian hosts, and as procyclic forms in tsetse flies. In trypanosomes, gene expression regulation depends heavily on post-transcriptional mechanisms. Both the RNA-binding protein RBP10 and glycosomal phosphoglycerate kinase PGKC are expressed only in mammalian-infective forms. RBP10 targets procyclic-specific mRNAs for destruction, while PGKC is required for bloodstream-form glycolysis. Developmental regulation of both is essential: expression of either RBP10 or PGKC in procyclic forms inhibits their proliferation. We show that the 3’-untranslated region of the RBP10 mRNA is extraordinarily long—7.3kb—and were able to identify six different sequences, scattered across the untranslated region, which can independently cause bloodstream-form-specific expression. The 3’-untranslated region of the PGKC mRNA, although much shorter, still contains two different regions, of 125 and 153nt, that independently gave developmental regulation. No short consensus sequences were identified that were enriched either within these regulatory regions, or when compared with other mRNAs with similar regulation, suggesting that more than one regulatory RNA-binding protein is important for repression of mRNAs in procyclic forms. We also identified regions, including an AU repeat, that increased expression in bloodstream forms, or suppressed it in both forms. Trypanosome mRNAs that encode RNA-binding proteins often have extremely extended 3’-untranslated regions. We suggest that one function of this might be to act as a fail-safe mechanism to ensure correct regulation even if mRNA processing or expression of trans regulators is defective.

The parasite Trypanosoma brucei causes sleeping sickness in humans, and nagana in cattle, and is transmitted by Tsetse flies. It grows in the bloodstream and tissue fluids of mammalian hosts, as "bloodstream forms", and as "procyclic forms" in the midgut of tsetse flies. Several hundred proteins are expressed in a stage-specific fashion, and this is essential for parasite survival in the different environments. RBP10 is an RNA-binding protein that is expressed only in bloodstream forms. It binds to procyclic-specific mRNAs, and causes their destruction. PGKC is an enzyme that is also specifically expressed in bloodstream forms. Developmental regulation of both is essential: expression of either RBP10 or PGKC in procyclic forms prevents their growth. The mRNAs encoding both proteins are very unstable in procyclic forms, and the sequences responsible are in an "untranslated region" of the mRNA—sequences that follow the part that codes for protein. We here show that the mRNA encoding PGKC has two regions that independently cause developmental regulation, and that the very long untranslated region of the RBP10 mRNA has no fewer than six regulatory regions, but there were no obvious similarities between them. We suggest that the presence of several different regulatory sequences in trypanosome mRNAs might be a fail-safe mechanism to ensure correct regulation.

RiboReport - benchmarking tools for ribosome profiling-based identification of open reading frames in bacteria

Small proteins encoded by short open reading frames (ORFs) with 50 codons or fewer are emerging as an important class of cellular macromolecules in diverse organisms. However, they often evade detection by proteomics or in silico methods. Ribosome profiling (Ribo-seq) has revealed widespread translation in genomic regions previously thought to be non-coding, driving the development of ORF detection tools using Ribo-seq data. However, only a handful of tools have been designed for bacteria, and these have not yet been systematically compared. Here, we aimed to identify tools that use Ribo-seq data to correctly determine the translational status of annotated bacterial ORFs and also discover novel translated regions with high sensitivity. To this end, we generated a large set of annotated ORFs from four diverse bacterial organisms, manually labeled for their translation status based on Ribo-seq data, which are available for future benchmarking studies. This set was used to investigate the predictive performance of seven Ribo-seq-based ORF detection tools (REPARATION_blast, DeepRibo, Ribo-TISH, PRICE, smORFer, ribotricer and SPECtre), as well as IRSOM, which uses coding potential and RNA-seq coverage only. DeepRibo and REPARATION_blast robustly predicted translated ORFs, including sORFs, with no significant difference for ORFs in close proximity to other genes versus stand-alone genes. However, no tool predicted a set of novel, experimentally verified sORFs with high sensitivity. Start codon predictions with smORFer show the value of initiation site profiling data to further improve the sensitivity of ORF prediction tools in bacteria. Overall, we find that bacterial tools perform well for sORF detection, although there is potential for improving their performance, applicability, usability and reproducibility.

Polysomal mRNA Association and Gene Expression in Trypanosoma brucei

Background: The contrasting physiological environments of  Trypanosoma brucei procyclic (insect vector) and bloodstream (mammalian host) forms necessitates deployment of different molecular processes and, therefore, changes in protein expression. Transcriptional regulation is unusual in  T. brucei because the arrangement of genes is polycistronic; however, genes which are transcribed together are subsequently cleaved into separate mRNAs by  trans-splicing. Following pre-mRNA processing, the regulation of mature mRNA stability is a tightly controlled cellular process. While many stage-specific transcripts have been identified, previous studies using RNA-seq suggest that changes in overall transcript level do not necessarily reflect the abundance of the corresponding protein.

Methods: To better understand the regulation of gene expression in  T. brucei, we performed a bioinformatic analysis of RNA-seq on total, sub-polysomal, and polysomal mRNA samples. We further cross-referenced our dataset with a previously published proteomics dataset to identify new protein coding sequences.

Results: Our analyses showed that several long non-coding RNAs are more abundant in the sub-polysome samples, which possibly implicates them in regulating cellular differentiation in  T. brucei. We also improved the annotation of the  T.brucei genome by identifying new putative protein coding transcripts that were confirmed by mass spectrometry data.

Conclusions: Several long non-coding RNAs are more abundant in the sub-polysome cellular fractions and might pay a role in the regulation of gene expression. We hope that these data will be of wide general interest, as well as being of specific value to researchers studying gene regulation expression and life stage transitions in  T. brucei.

Translation complex stabilization on messenger RNA and footprint profiling to study the RNA responses and dynamics of protein biosynthesis in the cells

During protein biosynthesis, ribosomes bind to messenger (m)RNA, locate its protein-coding information, and translate the nucleotide triplets sequentially as codons into the corresponding sequence of amino acids, forming proteins. Non-coding mRNA features, such as 5' and 3' untranslated regions (UTRs), start sites or stop codons of different efficiency, stretches of slower or faster code and nascent polypeptide interactions can alter the translation rates transcript-wise. Most of the homeostatic and signal response pathways of the cells converge on individual mRNA control, as well as alter the global translation output. Among the multitude of approaches to study translational control, one of the most powerful is to infer the locations of translational complexes on mRNA based on the mRNA fragments protected by these complexes from endonucleolytic hydrolysis, or footprints. Translation complex profiling by high-throughput sequencing of the footprints allows to quantify the transcript-wise, as well as global, alterations of translation, and uncover the underlying control mechanisms by attributing footprint locations and sizes to different configurations of the translational complexes. The accuracy of all footprint profiling approaches critically depends on the fidelity of footprint generation and many methods have emerged to preserve certain or multiple configurations of the translational complexes, often in challenging biological material. In this review, a systematic summary of approaches to stabilize translational complexes on mRNA for footprinting is presented and major findings are discussed. Future directions of translation footprint profiling are outlined, focusing on the fidelity and accuracy of inference of the native in vivo translation complex distribution on mRNA.

Extensive Translational Regulation through the Proliferative Transition of Trypanosoma cruzi Revealed by Multi-Omics

Trypanosoma cruzi is the etiological agent for Chagas disease, a neglected parasitic disease in Latin America. Gene transcription control governs the eukaryotic cell replication but is absent in trypanosomatids; thus, it must be replaced by posttranscriptional regulatory events. We investigated the entrance into the T. cruzi replicative cycle using ribosome profiling and proteomics on G₁/S epimastigote cultures synchronized with hydroxyurea. We identified 1,784 translationally regulated genes (change > 2, false-discovery rate [FDR] < 0.05) and 653 differentially expressed proteins (change > 1.5, FDR < 0.05), respectively. A major translational remodeling accompanied by an extensive proteome change is found, while the transcriptome remains largely unperturbed at the replicative entrance of the cell cycle. The differentially expressed genes comprise specific cell cycle processes, confirming previous findings while revealing candidate cell cycle regulators that undergo previously unnoticed translational regulation. Clusters of genes showing a coordinated regulation at translation and protein abundance share related biological functions such as cytoskeleton organization and mitochondrial metabolism; thus, they may represent posttranscriptional regulons. The translatome and proteome of the coregulated clusters change in both coupled and uncoupled directions, suggesting that complex cross talk between the two processes is required to achieve adequate protein levels of different regulons. This is the first simultaneous assessment of the transcriptome, translatome, and proteome of trypanosomatids, which represent a paradigm for the absence of transcriptional control. The findings suggest that gene expression chronology along the T. cruzi cell cycle is controlled mainly by translatome and proteome changes coordinated using different mechanisms for specific gene groups. IMPORTANCE Trypanosoma cruzi is an ancient eukaryotic unicellular parasite causing Chagas disease, a potentially life-threatening illness that affects 6 to 7 million people, mostly in Latin America. The antiparasitic treatments for the disease have incomplete efficacy and adverse reactions; thus, improved drugs are needed. We study the mechanisms governing the replication of the parasite, aiming to find differences with the human host, valuable for the development of parasite-specific antiproliferative drugs. Transcriptional regulation is essential for replication in most eukaryotes, but in trypanosomatids, it must be replaced by subsequent gene regulation steps since they lack transcription initiation control. We identified the genome-wide remodeling of mRNA translation and protein abundance during the entrance to the replicative phase of the cell cycle. We found that translation is strongly regulated, causing variation in protein levels of specific cell cycle processes, representing the first simultaneous study of the translatome and proteome in trypanosomatids.

Allele-specific assembly of a eukaryotic genome corrects apparent frameshifts and reveals a lack of nonsense-mediated mRNA decay

To date, most reference genomes represent a mosaic consensus sequence in which the homologous chromosomes are collapsed into one sequence. This approach produces sequence artefacts and impedes analyses of allele-specific mechanisms. Here, we report an allele-specific genome assembly of the diploid parasite Trypanosoma brucei and reveal allelic variants affecting gene expression. Using long-read sequencing and chromosome conformation capture data, we could assign 99.5% of all heterozygote variants to a specific homologous chromosome and build a 66 Mb long allele-specific genome assembly. The phasing of haplotypes allowed us to resolve hundreds of artefacts present in the previous mosaic consensus assembly. In addition, it revealed allelic recombination events, visible as regions of low allelic heterozygosity, enabling the lineage tracing of T. brucei isolates. Interestingly, analyses of transcriptome and translatome data of genes with allele-specific premature termination codons point to the absence of a nonsense-mediated decay mechanism in trypanosomes. Taken together, this study delivers a reference quality allele-specific genome assembly of T. brucei and demonstrates the importance of such assemblies for the study of gene expression control. We expect the new genome assembly will increase the awareness of allele-specific phenomena and provide a platform to investigate them.

Polysomal mRNA Association and Gene Expression in Trypanosoma brucei

Background: The contrasting physiological environments of  Trypanosoma brucei procyclic (insect vector) and bloodstream (mammalian host) forms necessitates deployment of different molecular processes and, therefore, changes in protein expression. Transcriptional regulation is unusual in  T. brucei because the arrangement of genes is polycistronic; however, genes which are transcribed together are subsequently cleaved into separate mRNAs by  trans-splicing. Following pre-mRNA processing, the regulation of mature mRNA stability is a tightly controlled cellular process. While many stage-specific transcripts have been identified, previous studies using RNA-seq suggest that changes in overall transcript level do not necessarily reflect the abundance of the corresponding protein. Methods: To better understand the regulation of gene expression in  T. brucei, we performed a bioinformatic analysis of RNA-seq on total, sub-polysomal, and polysomal mRNA samples. We further cross-referenced our dataset with a previously published proteomics dataset to identify new protein coding sequences. Results: Our analyses showed that several long non-coding RNAs are more abundant in the sub-polysome samples, which possibly implicates them in regulating cellular differentiation in  T. brucei. We also improved the annotation of the  T.brucei genome by identifying new putative protein coding transcripts that were confirmed by mass spectrometry data. Conclusions: Several long non-coding RNAs are more abundant in the sub-polysome cellular fractions and might pay a role in the regulation of gene expression. We hope that these data will be of wide general interest, as well as being of specific value to researchers studying gene regulation expression and life stage transitions in  T. brucei.

Polysomal mRNA Association and Gene Expression in Trypanosoma brucei

Background: The contrasting physiological environments of  Trypanosoma brucei procyclic (insect vector) and bloodstream (mammalian host) forms necessitates deployment of different molecular processes and, therefore, changes in protein expression. Transcriptional regulation is unusual in  T. brucei because the arrangement of genes is polycistronic; however, genes which are transcribed together are subsequently cleaved into separate mRNAs by  trans-splicing. Following pre-mRNA processing, the regulation of mature mRNA stability is a tightly controlled cellular process. While many stage-specific transcripts have been identified, previous studies using RNA-seq suggest that changes in overall transcript level do not necessarily reflect the abundance of the corresponding protein. Methods: To better understand the regulation of gene expression in  T. brucei, we performed a bioinformatic analysis of RNA-seq on total, sub-polysomal, and polysomal mRNA samples. We further cross-referenced our dataset with a previously published proteomics dataset to identify new protein coding sequences. Results: Our analyses showed that several long non-coding RNAs are more abundant in the sub-polysome samples, which possibly implicates them in regulating cellular differentiation in  T. brucei. We also improved the annotation of the  T.brucei genome by identifying new putative protein coding transcripts that were confirmed by mass spectrometry data. Conclusions: Several long non-coding RNAs are more abundant in the sub-polysome cellular fractions and might pay a role in the regulation of gene expression. We hope that these data will be of wide general interest, as well as being of specific value to researchers studying gene regulation expression and life stage transitions in  T. brucei.

Epitranscriptome machinery in Trypanosomatids: New players on the table?

Trypanosoma and Leishmania parasites cause devastating tropical diseases resulting in serious global health consequences. These organisms have complex life cycles with mammalian hosts and insect vectors. The parasites must, therefore, survive in different environments, demanding rapid physiological and metabolic changes. These responses depend upon regulation of gene expression, which primarily occurs posttranscriptionally. Altering the composition or conformation of RNA through nucleotide modifications is one posttranscriptional mechanism of regulating RNA fate and function, and modifications including N6-methyladenosine (m6A), N1-methyladenosine (m1A), N5-methylcytidine (m5C), N4-acetylcytidine (ac4C), and pseudouridine (Ψ), dynamically regulate RNA stability and translation in diverse organisms. Little is known about RNA modifications and their machinery in Trypanosomatids, but we hypothesize that they regulate parasite gene expression and are vital for survival. Here, we identified Trypanosomatid homologs for writers of m1A, m5C, ac4C, and Ψ and analyze their evolutionary relationships. We systematically review the evidence for their functions and assess their potential use as therapeutic targets. This work provides new insights into the roles of these proteins in Trypanosomatid parasite biology and treatment of the diseases they cause and illustrates that Trypanosomatids provide an excellent model system to study RNA modifications, their molecular, cellular, and biological consequences, and their regulation and interplay.

Nutrient availability regulates proline/alanine transporters in Trypanosoma brucei

Trypanosoma brucei is a species of unicellular parasite that can cause severe diseases in livestock and humans, including African trypanosomiasis and Chagas disease. Adaptation to diverse environments and changes in nutritional conditions is essential for T. brucei to establish an infection when changing hosts or during invasion of different host tissues. One such adaptation is the ability of T. brucei to rapidly switch its energy metabolism from glucose metabolism in the mammalian blood to proline catabolism in the insect stages and vice versa. However, the mechanisms that support the parasite's response to nutrient availability remain unclear. Using RNAseq and qRT-PCR, we investigated the response of T. brucei to amino acid or glucose starvation and found increased mRNA levels of several amino acid transporters, including all genes of the amino acid transporter AAT7-B subgroup. Functional characterization revealed that AAT7-B members are plasma membrane-localized in T. brucei and when expressed in Saccharomyces cerevisiae supported the uptake of proline, alanine, and cysteine, while other amino acids were poorly recognized. All AAT7-B members showed a preference for proline, which is transported with high or low affinity. RNAi-mediated AAT7-B downregulation resulted in a reduction of intracellular proline concentrations and growth arrest under low proline availability in cultured procyclic form parasites. Taken together, these results suggest a role of AAT7-B transporters in the response of T. brucei to proline starvation and proline catabolism.

Leishmania: Responding to environmental signals and challenges without regulated transcription

Here we describe the non-canonical control of gene expression in Leishmania, a single-cell parasite that is responsible for one of the major neglected tropical diseases. We discuss the lack of regulated RNA synthesis, the post-transcriptional gene regulation including RNA stability and regulated translation. We also show that genetic adaptations such as mosaic aneuploidy, gene copy number variations and DNA sequence polymorphisms are important means for overcoming drug challenge and environmental diversity. These mechanisms are discussed in the context of the unique flow of genetic information found in Leishmania and related protists.

Large-Scale Phylogenetic Analysis of Trypanosomatid Adenylate Cyclases Reveals Associations with Extracellular Lifestyle and Host-Pathogen Interplay

Receptor adenylate cyclases (RACs) on the surface of trypanosomatids are important players in the host–parasite interface. They detect still unidentified environmental signals that affect the parasites’ responses to host immune challenge, coordination of social motility, and regulation of cell division. A lesser known class of oxygen-sensing adenylate cyclases (OACs) related to RACs has been lost in trypanosomes and expanded mostly in Leishmania species and related insect-dwelling trypanosomatids. In this work, we have undertaken a large-scale phylogenetic analysis of both classes of adenylate cyclases (ACs) in trypanosomatids and the free-living Bodo saltans. We observe that the expanded RAC repertoire in trypanosomatids with a two-host life cycle is not only associated with an extracellular lifestyle within the vertebrate host, but also with a complex path through the insect vector involving several life cycle stages. In Trypanosoma brucei, RACs are split into two major clades, which significantly differ in their expression profiles in the mammalian host and the insect vector. RACs of the closely related Trypanosoma congolense are intermingled within these two clades, supporting early RAC diversification. Subspecies of T. brucei that have lost the capacity to infect insects exhibit high numbers of pseudogenized RACs, suggesting many of these proteins have become redundant upon the acquisition of a single-host life cycle. OACs appear to be an innovation occurring after the expansion of RACs in trypanosomatids. Endosymbiont-harboring trypanosomatids exhibit a diversification of OACs, whereas these proteins are pseudogenized in Leishmania subgenus Viannia. This analysis sheds light on how ACs have evolved to allow diverse trypanosomatids to occupy multifarious niches and assume various lifestyles.

The computational approaches of lncRNA identification based on coding potential: Status quo and challenges

Long noncoding RNAs (lncRNAs) make up a large proportion of transcriptome in eukaryotes, and have been revealed with many regulatory functions in various biological processes. When studying lncRNAs, the first step is to accurately and specifically distinguish them from the colossal transcriptome data with complicated composition, which contains mRNAs, lncRNAs, small RNAs and their primary transcripts. In the face of such a huge and progressively expanding transcriptome data, the in-silico approaches provide a practicable scheme for effectively and rapidly filtering out lncRNA targets, using machine learning and probability statistics. In this review, we mainly discussed the characteristics of algorithms and features on currently developed approaches. We also outlined the traits of some state-of-the-art tools for ease of operation. Finally, we pointed out the underlying challenges in lncRNA identification with the advent of new experimental data.

Molecular characterization of 60S ribosomal protein L12 of E. tenella

This study analyzed the large-subunit (60S) ribosomal protein L12 of Eimeria tenella (Et60s-RPL12). A full-length cDNA was cloned, and the recombinant protein was expressed in E. coli BL21 and inoculated in rabbits to produce the polyclonal antibody. Quantitative real-time polymerase chain reaction and western blotting were used to analyze the transcription levels of Et60s-RPL12 and translation levels in different developmental stages of E. tenella. The results showed that the mRNA transcription level of Et60s-RPL12 was highest in second-generation merozoites, whereas the translation level was highest in unsporulated oocysts. Indirect immunofluorescence showed that Et60s-RPL12 was localized to the anterior region and surface of sporozoites, except for the two refractile bodies. As the invasion of DF-1 cells progressed, fluorescence intensity was increased, and Et60s-RPL12 was localized to the parasitophorous vacuole membrane (PVM). The secretion assay results using staurosporine indicated that this protein was secreted, but not from micronemes. The role of Et60s-RPL12 in invasion was evaluated in vitro. The results of the invasion assay showed that polyclonal antibody inhibited host cell invasion by the parasite, which reached about 12%. However, the rate of invasion was not correlated with the concentration of IgG.

The RNA-associated proteins MKT1 and MKT1L form alternative PBP1-containing complexes in Trypanosoma brucei

Control of gene expression in kinetoplastids such as trypanosomes depends heavily on RNA-binding proteins that influence mRNA decay and translation. We previously showed that the trypanosome protein MKT1 forms a multicomponent protein complex: MKT1 interacts with PBP1, which in turn recruits LSM12 and poly(A)-binding protein. MKT1 is recruited to mRNAs by sequence-specific RNA-binding proteins, resulting in stabilization of the bound mRNA. We here show that PBP1, LSM12, and a 117-residue protein, XAC1 (Tb927.7.2780), are present in complexes that contain either MKT1 or an MKT1-like protein, MKT1L (Tb927.10.1490). All five proteins are present predominantly in the complexes, and we found evidence for a minor subset of complexes containing both MKT1 and MKT1L. XAC1-containing complexes reproducibly contained RNA-binding proteins that were previously found associated with MKT1. Moreover, XAC1- or MKT1-containing complexes specifically recruited one of the two poly(A)-binding proteins, PABP2, and one of the six cap-binding translation initiation complexes, EIF4E6-EIF4G5. Yeast two-hybrid assay results indicated that MKT1 directly interacts with EIF4G5. MKT1-PBP1 complexes can therefore interact with mRNAs via their poly(A) tails and caps, as well as through sequence-specific RNA-binding proteins. Correspondingly, MKT1 is associated with many mRNAs, although not with those encoding ribosomal proteins. Meanwhile, MKT1L resembles MKT1 at the C terminus but additionally features an N-terminal extension with low-complexity regions. Although MKT1L depletion inhibited cell proliferation, we found no evidence that it specifically interacts with RNA-binding proteins or mRNA. We speculate that MKT1L may compete with MKT1 for PBP1 binding and thereby modulate the function of MKT1-containing complexes.

Transcriptomic changes across the life cycle of Trypanosoma cruzi II

Trypanosoma cruzi is a flagellated protozoan that causes Chagas disease; it presents a complex life cycle comprising four morphological stages: epimastigote (EP), metacyclic trypomastigote (MT), cell-derived trypomastigote (CDT) and amastigote (AM). Previous transcriptomic studies on three stages (EPs, CDTs and AMs) have demonstrated differences in gene expressions among them; however, to the best of our knowledge, no studies have reported on gene expressions in MTs. Therefore, the present study compared differentially expressed genes (DEGs), and signaling pathway reconstruction in EPs, MTs, AMs and CDTs. The results revealed differences in gene expressions in the stages evaluated; these differences were greater between MTs and AMs-PTs. The signaling pathway that presented the highest number of DEGs in all the stages was associated with ribosomes protein profiles, whereas the other related pathways activated were processes related to energy metabolism from glucose, amino acid metabolism, or RNA regulation. However, the role of autophagy in the entire life cycle of T. cruzi and the presence of processes such as meiosis and homologous recombination in MTs (where the expressions of SPO11 and Rad51 plays a role) are crucial. These findings represent an important step towards the full understanding of the molecular basis during the life cycle of T. cruzi.

The mitochondrial phosphate carrier TbMCP11 is essential for mitochondrial function in the procyclic form of Trypanosoma brucei

Conserved amongst all eukaryotes is a family of mitochondrial carrier proteins (SLC25A) responsible for the import of various solutes across the inner mitochondrial membrane. We previously reported that the human parasite Trypanosoma brucei possesses 26 SLC25A proteins (TbMCPs) amongst which two, TbMCP11 and TbMCP8, were predicted to function as phosphate importers. The transport of inorganic phosphate into the mitochondrion is a prerequisite to drive ATP synthesis by substrate level and oxidative phosphorylation and thus crucial for cell viability. In this paper we describe the functional characterization of TbMCP11. In procyclic form T. brucei, the RNAi of TbMCP11 blocked ATP synthesis on mitochondrial substrates, caused a drop of the mitochondrial oxygen consumption and drastically reduced cell viability. The functional complementation in yeast and mitochondrial swelling experiments suggested a role for TbMCP11 as inorganic phosphate carrier. Interestingly, procyclic form T. brucei cells in which TbMCP11 was depleted displayed an inability to either replicate or divide the kinetoplast DNA, which resulted in a severe cytokinesis defect.

The large repertoire of 2'-O-methylation guided by C/D snoRNAs on Trypanosoma brucei rRNA

The parasite Trypanosoma brucei cycles between insect and mammalian hosts, and is the causative agent of sleeping sickness. Here, we performed genome-wide mapping of 2'-O-methylations (Nms) on trypanosome rRNA using three high-throughput sequencing methods; RibOxi-seq, RiboMeth-seq and 2'-OMe-seq. This is the first study using three genome-wide mapping approaches on rRNA from the same species showing the discrepancy among the methods. RibOxi-seq detects all the sites, but RiboMeth-seq is the only method to evaluate the level of a single Nm site. The sequencing revealed at least ninety-nine Nms guided by eighty-five snoRNAs among these thirty-eight Nms are trypanosome specific sites. We present the sequence and target of the C/D snoRNAs guiding on rRNA. This is the highest number of Nms detected to date on rRNA of a single cell parasite. Based on RiboMeth-seq, several Nm sites were found to be differentially regulated at the two stages of the parasite life cycle, the insect procyclic form (PCF) versus the bloodstream form (BSF) in the mammalian host.

Upstream ORFs Influence Translation Efficiency in the Parasite Trypanosoma cruzi

It is generally accepted that the presence of ORFs in the 5' untranslated region of eukaryotic transcripts modulates the production of proteins by controlling the translation initiation rate of the main CDS. In trypanosomatid parasites, which almost exclusively depend on post-transcriptional mechanisms to regulate gene expression, translation has been identified as a key step. However, the mechanisms of control of translation are not fully understood. In the present work, we have annotated the 5'UTRs of the Trypanosoma cruzi genome both in epimastigotes and metacyclic trypomastigotes and, using a stringent classification approach, we identified putative regulatory uORFs in about 9% of the analyzed 5'UTRs. The translation efficiency (TE) and translational levels of transcripts containing putative repressive uORFs were found to be significantly reduced. These findings are supported by the fact that proteomic methods only identify a low number of proteins coded by transcripts containing repressive uORF. We additionally show that AUG is the main translation initiator codon of repressive uORFs in T. cruzi. Interestingly, the decrease in TE is more pronounced when the uORFs overlaps the main CDS. In conclusion, we show that the presence of the uORF and features such as initiation codon and/or location of the uORFs may be acting to fine tune translation levels in these parasites.

Polysome Profiling and Metabolic Labeling Methods to Measure Translation in Trypanosoma brucei

The amount of a protein that is made in a cell is determined not only by the corresponding mRNA level but also by the efficiency with which the mRNA is translated. Very powerful transcriptome-wide methods are available to analyze both the density of ribosomes on each mRNA and the rate at which polypeptides are elongated. However, for many research questions, simpler, less expensive methods are more suitable. Here we describe two methods to assess the general translation status of cells: polysome profiling by sucrose density gradient centrifugation and metabolic labeling using radioactive amino acids. Both methods can also be used to examine translation of individual mRNAs.

Organising the cell cycle in the absence of transcriptional control: Dynamic phosphorylation co-ordinates the Trypanosoma brucei cell cycle post-transcriptionally

The cell division cycle of the unicellular eukaryote Trypanosome brucei is tightly regulated despite the paucity of transcriptional control that results from the arrangement of genes in polycistronic units and lack of dynamically regulated transcription factors. To identify the contribution of dynamic phosphorylation to T. brucei cell cycle control we have combined cell cycle synchronisation by centrifugal elutriation with quantitative phosphoproteomic analysis. Cell cycle regulated changes in phosphorylation site abundance (917 sites, average 5-fold change) were more widespread and of a larger magnitude than changes in protein abundance (443 proteins, average 2-fold change) and were mostly independent of each other. Hierarchical clustering of co-regulated phosphorylation sites according to their cell cycle profile revealed that a bulk increase in phosphorylation occurs across the cell cycle, with a significant enrichment of known cell cycle regulators and RNA binding proteins (RBPs) within the largest clusters. Cell cycle regulated changes in essential cell cycle kinases are temporally co-ordinated with differential phosphorylation of components of the kinetochore and eukaryotic initiation factors, along with many RBPs not previously linked to the cell cycle such as eight PSP1-C terminal domain containing proteins. The temporal profiles demonstrate the importance of dynamic phosphorylation in co-ordinating progression through the cell cycle, and provide evidence that RBPs play a central role in post-transcriptional regulation of the T. brucei cell cycle. Data are available via ProteomeXchange with identifier PXD013488.

Mapping benznidazole resistance in trypanosomatids and exploring evolutionary histories of nitroreductases and ABCG transporter protein sequences

The nitro-heterocyclic compound benznidazole (BZ) is the first-line drug for the treatment of Chagas disease, caused by the protozoan Trypanosoma cruzi. However, therapeutic failures are common for reasons that include the influences of parasite and host genetics, the effects of toxicity on adherence to treatment, and difficulties in demonstrating parasitological cure. To obtain information on the origin of the resistance to BZ and eliminate from the scenery the participation of the host, initially we mapped the susceptibility to the drug in thirteen species of seven genera of the family Trypanosomatidae. We verified that all Trypanosoma species are sensitive to low concentrations of the drug (IC₅₀ 2.7 to 25 µM) while Non-Trypanosoma species are highly resistant to these concentrations. The two groups of parasites correspond to the major phylogenetic lineages of trypanosomatids. Next, we searched in the trypanosomatid genome databases homologs of two type-I nitroreductases (NTR-1 and OYE) and an ABC transporter (ABCG1) that have been associated with BZ resistance in T. cruzi. The predicted proteins were characterized regarding domains and used for phylogenetic analyses. Homologous NTR-1 genes were found in all trypanosomatids investigated and the structural characteristics of the enzyme suggest that it may be functional. OYE genes were absent in BZ-sensitive African trypanosomes, which excludes the participation of this enzyme in BZ bio-activation. Two copies of ABCG1 genes were observed in most BZ resistant species, while Trypanosoma species exhibit only one copy per haploid genome. Functional studies are required to verify the involvement of these genes in BZ resistance. In addition, since multiple mechanisms can contribute to BZ susceptibility, our study poses a range of organisms highly resistant to BZ in which these aspects can be investigated. Preliminary studies on BZ uptake indicate marked differences between BZ-sensitive and BZ-resistant species.

On the Function of Trans-Splicing: No Evidence for Widespread Proteome Diversification in Trypanosomes

A long-standing mystery of genomic/transcriptomic structure involves spliced leader trans-splicing (SLTS), in which short RNA “tags” transcribed from a distinct genomic locus is added near the 5′ end of RNA transcripts by the spliceosome. SLTS has been observed in diverse eukaryotes in a phylogenetic pattern implying recurrent independent evolution. This striking convergence suggests important functions for SLTS, however no general novel function is known. Recent findings of frequent alternative SLTS (ALT-TS) suggest that ALT-TS could impart widespread functionality. Here, we tested the hypothesis that ALT-TS diversifies proteomes by comparing splicing patterns in orthologous genes between two deeply diverged trypanosome parasites. We also tested proteome diversification functions of ALT-TS by utilizing ribosome profiling sequence data. Finally, we investigated ALT-TS as a mechanism to regulate the expression of unproductive transcripts. Although our results indicate the functional importance of some cases of trans-splicing, we find no evidence for the hypothesis that proteome diversification is a general function of trans-splicing.

GWIPS-viz: 2018 update

The GWIPS-viz browser (http://gwips.ucc.ie/) is an on-line genome browser which is tailored for exploring ribosome profiling (Ribo-seq) data. Since its publication in 2014, GWIPS-viz provides Ribo-seq data for an additional 14 genomes bringing the current total to 23. The integration of new Ribo-seq data has been automated thereby increasing the number of available tracks to 1792, a 10-fold increase in the last three years. The increase is particularly substantial for data derived from human sources. Following user requests, we added the functionality to download these tracks in bigWig format. We also incorporated new types of data (e.g. TCP-seq) as well as auxiliary tracks from other sources that help with the interpretation of Ribo-seq data. Improvements in the visualization of the data have been carried out particularly for bacterial genomes where the Ribo-seq data are now shown in a strand specific manner. For higher eukaryotic datasets, we provide characteristics of individual datasets using the RUST program which includes the triplet periodicity, sequencing biases and relative inferred A-site dwell times. This information can be used for assessing the quality of Ribo-seq datasets. To improve the power of the signal, we aggregate Ribo-seq data from several studies into Global aggregate tracks for each genome.

Ribosome Profiling in Trypanosomatids

Ribosomes are the machinery responsible for reading mRNAs and translating them into proteins. The ribosome profiling approach is based on high-throughput sequencing of ribosome-protected mRNAs. RNAs not harboring ribosomes are removed by nuclease digestion leaving the so-called ribosome "footprints." The purified "footprint" RNA molecules are processed into DNA libraries and their individual abundance is determined by deep sequencing. Ribosome profiling reveals the portion of transcripts which are actually protein-coding and can be used for differential gene expression analysis addressing rates of protein synthesis, and translational control and efficiency.

Regulation of gene expression in trypanosomatids: living with polycistronic transcription

In trypanosomes, RNA polymerase II transcription is polycistronic and individual mRNAs are excised by trans-splicing and polyadenylation. The lack of individual gene transcription control is compensated by control of mRNA processing, translation and degradation. Although the basic mechanisms of mRNA decay and translation are evolutionarily conserved, there are also unique aspects, such as the existence of six cap-binding translation initiation factor homologues, a novel decapping enzyme and an mRNA stabilizing complex that is recruited by RNA-binding proteins. High-throughput analyses have identified nearly a hundred regulatory mRNA-binding proteins, making trypanosomes valuable as a model system to investigate post-transcriptional regulation.

Simultaneous Ribosome Profiling of Human Host Cells Infected with Toxoplasma gondii

Toxoplasma gondii is a single-celled parasite that has infected up to one-third of the world’s population. Significant overhauls in gene expression in both the parasite and the host cell accompany parasite invasion, and a better understanding of these changes may lead to the development of new therapeutic agents. In this study, we employed ribosome profiling to determine the changes that occur at the levels of transcription and translation in both the parasite and the infected host cell at the same time. We discovered features of Toxoplasma mRNAs that suggest a means for controlling parasite gene expression under stressful conditions. We also show that differences in host gene expression occur depending on whether they are confluent or not. Our findings demonstrate the feasibility of using ribosomal profiling to interrogate the host-parasite dynamic under a variety of conditions.

Toxoplasma gondii is a ubiquitous obligate intracellular parasite that infects the nucleated cells of warm-blooded animals. From within the parasitophorous vacuole in which they reside, Toxoplasma tachyzoites secrete an arsenal of effector proteins that can reprogram host gene expression to facilitate parasite survival and replication. Gaining a better understanding of how host gene expression is altered upon infection is central for understanding parasite strategies for host invasion and for developing new parasite therapies. Here, we applied ribosome profiling coupled with mRNA measurements to concurrently study gene expression in the parasite and in host human foreskin fibroblasts. By examining the parasite transcriptome and translatome, we identified potential upstream open reading frames that may permit the stress-induced preferential translation of parasite mRNAs. We also determined that tachyzoites reduce host death-associated pathways and increase survival, proliferation, and motility in both quiescent and proliferative host cell models of infection. Additionally, proliferative cells alter their gene expression in ways that are consistent with massive transcriptional rewiring, while quiescent cells were best characterized by reentry into the cell cycle. We also identified a translational control regimen consistent with mechanistic target of rapamycin (mTOR) activation in quiescent cells and, to a lesser degree, in proliferative cells. This study illustrates the utility of the method for dissection of gene expression programs simultaneously in the parasite and host.

IMPORTANCEToxoplasma gondii is a single-celled parasite that has infected up to one-third of the world’s population. Significant overhauls in gene expression in both the parasite and the host cell accompany parasite invasion, and a better understanding of these changes may lead to the development of new therapeutic agents. In this study, we employed ribosome profiling to determine the changes that occur at the levels of transcription and translation in both the parasite and the infected host cell at the same time. We discovered features of Toxoplasma mRNAs that suggest a means for controlling parasite gene expression under stressful conditions. We also show that differences in host gene expression occur depending on whether they are confluent or not. Our findings demonstrate the feasibility of using ribosomal profiling to interrogate the host-parasite dynamic under a variety of conditions.

Comparative molecular cell biology of phototrophic euglenids and parasitic trypanosomatids sheds light on the ancestor of Euglenozoa

Parasitic trypanosomatids and phototrophic euglenids are among the most extensively studied euglenozoans. The phototrophic euglenid lineage arose relatively recently through secondary endosymbiosis between a phagotrophic euglenid and a prasinophyte green alga that evolved into the euglenid secondary chloroplast. The parasitic trypanosomatids (i.e. Trypanosoma spp. and Leishmania spp.) and the freshwater phototrophic euglenids (i.e. Euglena gracilis) are the most evolutionary distant lineages in the Euglenozoa phylogenetic tree. The molecular and cell biological traits they share can thus be considered as ancestral traits originating in the common euglenozoan ancestor. These euglenozoan ancestral traits include common mitochondrial presequence motifs, respiratory chain complexes containing various unique subunits, a unique ATP synthase structure, the absence of mitochondria-encoded transfer RNAs (tRNAs), a nucleus with a centrally positioned nucleolus, closed mitosis without dissolution of the nuclear membrane and nucleoli, a nuclear genome containing the unusual 'J' base (β-D-glucosyl-hydroxymethyluracil), processing of nucleus-encoded precursor messenger RNAs (pre-mRNAs) via spliced-leader RNA (SL-RNA) trans-splicing, post-transcriptional gene silencing by the RNA interference (RNAi) pathway and the absence of transcriptional regulation of nuclear gene expression. Mitochondrial uridine insertion/deletion RNA editing directed by guide RNAs (gRNAs) evolved in the ancestor of the kinetoplastid lineage. The evolutionary origin of other molecular features known to be present only in either kinetoplastids (i.e. polycistronic transcripts, compaction of nuclear genomes) or euglenids (i.e. monocistronic transcripts, huge genomes, many nuclear cis-spliced introns, polyproteins) is unclear.

Functional genomics in sand fly-derived Leishmania promastigotes

Background

Leishmania development in the sand fly gut leads to highly infective forms called metacyclic promastigotes. This process can be routinely mimicked in culture. Gene expression–profiling studies by transcriptome analysis have been performed with the aim of studying promastigote forms in the sand fly gut, as well as differences between sand fly–and culture-derived promastigotes.

Findings

Transcriptome analysis has revealed the crucial role of the microenvironment in parasite development within the sand fly gut because substantial differences and moderate correlation between the transcriptomes of cultured and sand fly–derived promastigotes have been found. Sand fly–derived metacyclics are more infective than metacyclics in culture. Therefore, some caution should be exercised when using cultured promastigotes, depending on the experimental design. The most remarkable examples are the hydrophilic acidic surface protein/small endoplasmic reticulum protein (HASP/SHERP) cluster, the glycoprotein 63 (gp63), and autophagy genes, which are up-regulated in sand fly–derived promastigotes compared with cultured promastigotes. Because HASP/SHERP genes are up-regulated in nectomonad and metacyclic promastigotes in the sand fly, the encoded proteins are not metacyclic specific. Metacyclic promastigotes are distinguished by morphology and high infectivity. Isolating them from the sand fly gut is not exempt from technical difficulty, because other promastigote forms remain in the gut even 15 days after infection. Leishmania major procyclic promastigotes within the sand fly gut up-regulate genes involved in cell cycle regulation and glucose catabolism, whereas metacyclics increase transcript levels of fatty acid biosynthesis and ATP-coupled proton transport genes. Most parasite's signal transduction pathways remain uncharacterized. Future elucidation may improve understanding of parasite development, particularly signaling molecule-encoding genes in sand fly versus culture and between promastigote forms in the sand fly gut.

Conclusions

Transcriptome analysis has been demonstrated to be technically efficacious to study differential gene expression in sand fly gut promastigote forms. Transcript and protein levels are not well correlated in these organisms (approximately 25% quantitative coincidences), especially under stress situations and at differentiation processes. However, transcript and protein levels behave similarly in approximately 60% of cases from a qualitative point of view (increase, decrease, or no variation). Changes in translational efficiency observed in other trypanosomatids strongly suggest that the differences are due to translational regulation and regulation of the steady-state protein levels. The lack of low-input sample strategies does not allow translatome and proteome analysis of sand fly–derived promastigotes so far.

Ribosome Profiling: Global Views of Translation

The translation of messenger RNA (mRNA) into protein and the folding of the resulting protein into an active form are prerequisites for virtually every cellular process and represent the single largest investment of energy by cells. Ribosome profiling-based approaches have revolutionized our ability to monitor every step of protein synthesis in vivo, allowing one to measure the rate of protein synthesis across the proteome, annotate the protein coding capacity of genomes, monitor localized protein synthesis, and explore cotranslational folding and targeting. The rich and quantitative nature of ribosome profiling data provides an unprecedented opportunity to explore and model complex cellular processes. New analytical techniques and improved experimental protocols will provide a deeper understanding of the factors controlling translation speed and its impact on protein function and cell physiology as well as the role of ribosomal RNA and mRNA modifications in regulating translation.