Spliced leader trapping reveals widespread alternative splicing patterns in the highly dynamic transcriptome of Trypanosoma brucei

Unique features of conventional and nonconventional introns in Euglena gracilis

BACKGROUND

Nuclear introns in Euglenida have been understudied. This study aimed to investigate nuclear introns in Euglenida by identifying a large number of introns in Euglena gracilis (E. gracilis), including cis-spliced conventional and nonconventional introns, as well as trans-spliced outrons. We also examined the sequence characteristics of these introns.

RESULTS

A total of 28,337 introns and 11,921 outrons were identified. Conventional and nonconventional introns have distinct splice site features; the former harbour canonical GT/C-AG splice sites, whereas the latter are capable of forming structured motifs with their terminal sequences. We observed that short introns had a preference for canonical GT-AG introns. Notably, conventional introns and outrons in E. gracilis exhibited a distinct cytidine-rich polypyrimidine tract, in contrast to the thymidine-rich tracts observed in other organisms. Furthermore, the SL-RNAs in E. gracilis, as well as in other trans-splicing species, can form a recently discovered motif called the extended U6/5' ss duplex with the respective U6s. We also describe a novel type of alternative splicing pattern in E. gracilis. The tandem repeat sequences of introns in this protist were determined, and their contents were comparable to those in humans.

CONCLUSIONS

Our findings highlight the unique features of E. gracilis introns and provide insights into the splicing mechanism of these introns, as well as the genomics and evolution of Euglenida.

Molecular insights and antibody response to Dr20/22 in dogs naturally infected with Dirofilaria repens

Subcutaneous dirofilariasis, caused by the parasitic nematode Dirofilaria repens, is a growing concern in Europe, affecting both dogs and humans. This study focused on D. repens Dr20/22, a protein encoded by an alt (abundant larval transcript) gene family. While well-documented in L3 larvae of other filariae species, this gene family had not been explored in dirofilariasis. The research involved cloning Dr20/22 cDNA, molecular characterization, and evaluating its potential application in the diagnosis of dirofilariasis. Although Real-Time analysis revealed mRNA expression in both adult worms and microfilariae, the native protein remained undetected in lysates from both developmental stages. This suggests the protein's specificity for L3 larvae and may be related to a process called SLTS (spliced leader trans-splicing), contributing to stage-specific gene expression. The specificity of the antigen for invasive larvae positions it as a promising early marker for dirofilariasis. However, ELISA tests using sera from infected and uninfected dogs indicated limited diagnostic utility. While further research is required, our findings contribute to a deeper understanding of the molecular and immunological aspects of host-parasite interactions and could offer insights into the parasite's strategies for evading the immune system.

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.

Chemical Inhibition of Bromodomain Proteins in Insect-Stage African Trypanosomes Perturbs Silencing of the Variant Surface Glycoprotein Repertoire and Results in Widespread Changes in the Transcriptome

The eukaryotic protozoan parasite Trypanosoma brucei is transmitted by the tsetse fly to both humans and animals, where it causes a fatal disease called African trypanosomiasis. While the parasite lacks canonical DNA sequence-specific transcription factors, it does possess histones, histone modifications, and proteins that write, erase, and read histone marks. Chemical inhibition of chromatin-interacting bromodomain proteins has previously been shown to perturb bloodstream specific trypanosome processes, including silencing of the variant surface glycoprotein (VSG) genes and immune evasion. Transcriptomic changes that occur in bromodomain-inhibited bloodstream parasites mirror many of the changes that occur as parasites developmentally progress from the bloodstream to the insect stage. We performed transcriptome sequencing (RNA-seq) time courses to determine the effects of chemical bromodomain inhibition in insect-stage parasites using the compound I-BET151. We found that treatment with I-BET151 causes large changes in the transcriptome of insect-stage parasites and also perturbs silencing of VSG genes. The transcriptomes of bromodomain-inhibited parasites share some features with early metacyclic-stage parasites in the fly salivary gland, implicating bromodomain proteins as important for regulating transcript levels for developmentally relevant genes. However, the downregulation of surface procyclin protein that typically accompanies developmental progression is absent in bromodomain-inhibited insect-stage parasites. We conclude that chemical modulation of bromodomain proteins causes widespread transcriptomic changes in multiple trypanosome life cycle stages. Understanding the gene-regulatory processes that facilitate transcriptome remodeling in this highly diverged eukaryote may shed light on how these mechanisms evolved. IMPORTANCE The disease African trypanosomiasis imposes a severe human and economic burden for communities in sub-Saharan Africa. The parasite that causes the disease is transmitted to the bloodstream of a human or ungulate via the tsetse fly. Because the environments of the fly and the bloodstream differ, the parasite modulates the expression of its genes to accommodate two different lifestyles in these disparate niches. Perturbation of bromodomain proteins that interact with histone proteins around which DNA is wrapped (chromatin) causes profound changes in gene expression in bloodstream-stage parasites. This paper reports that gene expression is also affected by chemical bromodomain inhibition in insect-stage parasites but that the genes affected differ depending on life cycle stage. Because trypanosomes diverged early from model eukaryotes, an understanding of how trypanosomes regulate gene expression may lend insight into how gene-regulatory mechanisms evolved. This could also be leveraged to generate new therapeutic strategies.

The RNA export factor TbMex67 connects transcription and RNA export in Trypanosoma brucei and sets boundaries for RNA polymerase I

TbMex67 is the major mRNA export factor known to date in trypanosomes, forming part of the docking platform within the nuclear pore. To explore its role in co-transcriptional mRNA export, recently reported in Trypanosoma brucei, pulse labelling of nascent RNAs with 5-ethynyl uridine (5-EU) was performed with cells depleted of TbMex67 and complemented with a dominant-negative mutant (TbMex67-DN). RNA polymerase (Pol) II transcription was unaffected, but the procyclin loci, which encode mRNAs transcribed by Pol I from internal sites on chromosomes 6 and 10, showed increased levels of 5-EU incorporation. This was due to Pol I readthrough transcription, which proceeded beyond the procyclin and procyclin-associated genes up to the Pol II transcription start site on the opposite strand. Complementation by TbMex67-DN also increased Pol I-dependent formation of R-loops and γ-histone 2A foci. The DN mutant exhibited reduced nuclear localisation and binding to chromatin compared to wild-type TbMex67. Together with its interaction with chromatin remodelling factor TbRRM1 and Pol II, and transcription-dependent association of Pol II with nucleoporins, our findings support a role for TbMex67 in connecting transcription and export in T. brucei. In addition, TbMex67 stalls readthrough by Pol I in specific contexts, thereby limiting R-loop formation and replication stress.

Aminoacyl tRNA Synthetases: Implications of Structural Biology in Drug Development against Trypanosomatid Parasites

The ensemble of aminoacyl tRNA synthetases is regarded as a key component of the protein translation machinery. With the progressive increase in structure-based studies on tRNA synthetase-ligand complexes, the detailed picture of these enzymes is becoming clear. Having known their critical role in deciphering the genetic code in a living system, they have always been chosen as one of the important targets for development of antimicrobial drugs. Later on, the role of aminoacyl tRNA synthetases (aaRSs) on the survivability of trypanosomatids has also been validated. It became evident through several gene knockout studies that targeting even one of these enzymes affected parasitic growth drastically. Such successful studies have inspired researchers to search for inhibitors that could specifically target trypanosomal aaRSs, and their never-ending efforts have provided fruitful results. Taking all such studies into consideration, these macromolecules of prime importance deserve further investigation for the development of drugs that cure spectrum of infections caused by trypanosomatids. In this review, we have compiled advancements of over a decade that have taken place in the pursuit of devising drugs by using trypanosomatid aaRSs as a major target of interest. Several of these inhibitors work on an exemplary low concentration range without posing any threat to the mammalian cells which is a very critical aspect of the drug discovery process. Advancements have been made in terms of using structural biology as an important tool to analyze the architecture of the trypanosomatids aaRSs and concoction of inhibitors with augmented specificities toward their targets. Some of the inhibitors that have been tested on other parasites successfully but their efficacy has so far not been validated against these trypanosomatids have also been appended.

Trans-splicing in the cestode Hymenolepis microstoma is constitutive across the life cycle and depends on gene structure and composition

Spliced leader (SL) trans-splicing is a key process during mRNA maturation of many eukaryotes, in which a short sequence (SL) is transferred from a precursor SL-RNA into the 5' region of an immature mRNA. This mechanism is present in flatworms, in which it is known to participate in the resolution of polycistronic transcripts. However, most trans-spliced transcripts are not part of operons, and it is not clear if this process may participate in additional regulatory mechanisms in this group. In this work, we present a comprehensive analysis of SL trans-splicing in the model cestode Hymenolepis microstoma. We identified four different SL-RNAs which are indiscriminately trans-spliced to 622 gene models. SL trans-splicing is enriched in constitutively expressed genes and does not appear to be regulated throughout the life cycle. Operons represented at least 20% of all detected trans-spliced gene models, showed conservation to those of the cestode Echinococcus multilocularis, and included complex loci such as an alternative operon (processed as either a single gene through cis-splicing or as two genes of a polycistron). Most insertion sites were identified in the 5' untranslated region (UTR) of monocistronic genes. These genes frequently contained introns in the 5' UTR, in which trans-splicing used the same acceptor sites as cis-splicing. These results suggest that, unlike other eukaryotes, trans-splicing is associated with internal intronic promoters in the 5' UTR, resulting in transcripts with strong splicing acceptor sites without competing cis-donor sites, pointing towards a simple mechanism driving the evolution of novel SL insertion sites.

slORFfinder: a tool to detect open reading frames resulting from trans-splicing of spliced leader sequences

Trans-splicing of a spliced leader (SL) to the 5' ends of mRNAs is used to produce mature mRNAs in several phyla of great importance to human health and the marine ecosystem. One of the consequences of the addition of SL sequences is the change or disruption of the open reading frames (ORFs) in the recipient transcripts. Given that most SL sequences have one or more of the trinucleotide NUG, including AUG in flatworms, trans-splicing of SL sequences can potentially supply a start codon to create new ORFs, which we refer to as slORFs, in the recipient mRNAs. Due to the lack of a tool to precisely detect them, slORFs were usually neglected in previous studies. In this work, we present the tool slORFfinder, which automatically links the SL sequences to the recipient mRNAs at the trans-splicing sites identified from SL-containing reads of RNA-Seq and predicts slORFs according to the distribution of ribosome-protected footprints (RPFs) on the trans-spliced transcripts. By applying this tool to the analyses of nematodes, ascidians and euglena, whose RPFs are publicly available, we find wide existence of slORFs in these taxa. Furthermore, we find that slORFs are generally translated at higher levels than the annotated ORFs in the genomes, suggesting they might have important functions. Overall, this study provides a tool, slORFfinder (https://github.com/songbo446/slORFfinder), to identify slORFs, which can enhance our understanding of ORFs in taxa with SL machinery.

Sequences and proteins that influence mRNA processing in Trypanosoma brucei: Evolutionary conservation of SR-domain and PTB protein functions

BACKGROUND

Spliced leader trans splicing is the addition of a short, capped sequence to the 5' end of mRNAs. It is widespread in eukaryotic evolution, but factors that influence trans splicing acceptor site choice have been little investigated. In Kinetoplastids, all protein-coding mRNAs are 5' trans spliced. A polypyrimidine tract is usually found upstream of the AG splice acceptor, but there is no branch point consensus; moreover, splicing dictates polyadenylation of the preceding mRNA, which is a validated drug target.

METHODOLOGY AND PRINCIPAL FINDINGS

We here describe a trans splicing reporter system that can be used for studies and screens concerning the roles of sequences and proteins in processing site choice and efficiency. Splicing was poor with poly(U) tracts less than 9 nt long, and was influenced by an intergenic region secondary structure. A screen for signals resulted in selection of sequences that were on average 45% U and 35% C. Tethering of either the splicing factor SF1, or the cleavage and polyadenylation factor CPSF3 within the intron stimulated processing in the correct positions, while tethering of two possible homologues of Opisthokont PTB inhibited processing. In contrast, tethering of SR-domain proteins RBSR1, RBSR2, or TSR1 or its interaction partner TSR1IP, promoted use of alternative signals upstream of the tethering sites. RBSR1 interacts predominantly with proteins implicated in splicing, whereas the interactome of RBSR2 is more diverse.

CONCLUSIONS

Our selectable constructs are suitable for screens of both sequences, and proteins that affect mRNA processing in T. brucei. Our results suggest that the functions of PTB and SR-domain proteins in splice site definition may already have been present in the last eukaryotic common ancestor.

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.

A diversified and segregated mRNA spliced-leader system in the parasitic Perkinsozoa

Spliced-leader trans-splicing (SLTS) has been described in distantly related eukaryotes and acts to mark mRNAs with a short 5′ exon, giving different mRNAs identical 5′ sequence-signatures. The function of these systems is obscure. Perkinsozoa encompasses a diversity of parasitic protists that infect bivalves, toxic-tide dinoflagellates, fish and frog tadpoles. Here, we report considerable sequence variation in the SLTS-system across the Perkinsozoa and find that multiple variant SLTS-systems are encoded in parallel in the ecologically important Perkinsozoa parasite Parvilucifera sinerae. These results demonstrate that the transcriptome of P. sinerae is segregated based on the addition of different spliced-leader (SL) exons. This segregation marks different gene categories, suggesting that SL-segregation relates to functional differentiation of the transcriptome. By contrast, both sets of gene categories are present in the single SL-transcript type sampled from Maranthos, implying that the SL-segregation of the Parvilucifera transcriptome is a recent evolutionary innovation. Furthermore, we show that the SLTS-system marks a subsection of the transcriptome with increased mRNA abundance and includes genes that encode the spliceosome system necessary for SLTS-function. Collectively, these data provide a picture of how the SLTS-systems can vary within a major evolutionary group and identify how additional transcriptional-complexity can be achieved through SL-segregation.

Paving the Way: Contributions of Big Data to Apicomplexan and Kinetoplastid Research

In the age of big data an important question is how to ensure we make the most out of the resources we generate. In this review, we discuss the major methods used in Apicomplexan and Kinetoplastid research to produce big datasets and advance our understanding of Plasmodium, Toxoplasma, Cryptosporidium, Trypanosoma and Leishmania biology. We debate the benefits and limitations of the current technologies, and propose future advancements that may be key to improving our use of these techniques. Finally, we consider the difficulties the field faces when trying to make the most of the abundance of data that has already been, and will continue to be, generated.

Genomic Occupancy of the Bromodomain Protein Bdf3 Is Dynamic during Differentiation of African Trypanosomes from Bloodstream to Procyclic Forms

Trypanosoma brucei, the causative agent of human and animal African trypanosomiasis, cycles between a mammalian host and a tsetse fly vector. The parasite undergoes huge changes in morphology and metabolism during adaptation to each host environment. These changes are reflected in the different transcriptomes of parasites living in each host. However, it remains unclear whether chromatin-interacting proteins help mediate these changes. Bromodomain proteins localize to transcription start sites in bloodstream parasites, but whether the localization of bromodomain proteins changes as parasites differentiate from bloodstream to insect stages remains unknown. To address this question, we performed cleavage under target and release using nuclease (CUT&RUN) against bromodomain protein 3 (Bdf3) in parasites differentiating from bloodstream to insect forms. We found that Bdf3 occupancy at most loci increased at 3 h following onset of differentiation and decreased thereafter. A number of sites with increased bromodomain protein occupancy lie proximal to genes with altered transcript levels during differentiation, such as procyclins, procyclin-associated genes, and invariant surface glycoproteins. Most Bdf3-occupied sites are observed throughout differentiation. However, one site appears de novo during differentiation and lies proximal to the procyclin gene locus housing genes essential for remodeling surface proteins following transition to the insect stage. These studies indicate that occupancy of chromatin-interacting proteins is dynamic during life cycle stage transitions and provide the groundwork for future studies on the effects of changes in bromodomain protein occupancy. Additionally, the adaptation of CUT&RUN for Trypanosoma brucei provides other researchers with an alternative to chromatin immunoprecipitation (ChIP).

IMPORTANCE The parasite Trypanosoma brucei is the causative agent of human and animal African trypanosomiasis (sleeping sickness). Trypanosomiasis, which affects humans and cattle, is fatal if untreated. Existing drugs have significant side effects. Thus, these parasites impose a significant human and economic burden in sub-Saharan Africa, where trypanosomiasis is endemic. T. brucei cycles between the mammalian host and a tsetse fly vector, and parasites undergo huge changes in morphology and metabolism to adapt to different hosts. Here, we show that DNA-interacting bromodomain protein 3 (Bdf3) shows changes in occupancy at its binding sites as parasites transition from the bloodstream to the insect stage. Additionally, a new binding site appears near the locus responsible for remodeling of parasite surface proteins during transition to the insect stage. Understanding the mechanisms behind host adaptation is important for understanding the life cycle of the parasite.

N6-methyladenosine in poly(A) tails stabilize VSG transcripts

RNA modifications are important regulators of gene expression¹. In Trypanosoma brucei, transcription is polycistronic and thus most regulation happens post-transcriptionally². N⁶-methyladenosine (m⁶A) has been detected in this parasite, but its function remains unknown³. Here we found that m⁶A is enriched in 342 transcripts using RNA immunoprecipitation, with an enrichment in transcripts encoding variant surface glycoproteins (VSGs). Approximately 50% of the m⁶A is located in the poly(A) tail of the actively expressed VSG transcripts. m⁶A residues are removed from the VSG poly(A) tail before deadenylation and mRNA degradation. Computational analysis revealed an association between m⁶A in the poly(A) tail and a 16-mer motif in the 3' untranslated region of VSG genes. Using genetic tools, we show that the 16-mer motif acts as a cis-acting motif that is required for inclusion of m⁶A in the poly(A) tail. Removal of this motif from the 3' untranslated region of VSG genes results in poly(A) tails lacking m⁶A, rapid deadenylation and mRNA degradation. To our knowledge, this is the first identification of an RNA modification in the poly(A) tail of any eukaryote, uncovering a post-transcriptional mechanism of gene regulation.

Monodopsis and Vischeria Genomes Shed New Light on the Biology of Eustigmatophyte Algae

Members of eustigmatophyte algae, especially Nannochloropsis and Microchloropsis, have been tapped for biofuel production owing to their exceptionally high lipid content. Although extensive genomic, transcriptomic, and synthetic biology toolkits have been made available for Nannochloropsis and Microchloropsis, very little is known about other eustigmatophytes. Here we present three near-chromosomal and gapless genome assemblies of Monodopsis strains C73 and C141 (60 Mb) and Vischeria strain C74 (106 Mb), which are the sister groups to Nannochloropsis and Microchloropsis in the order Eustigmatales. These genomes contain unusually high percentages of simple repeats, ranging from 12% to 21% of the total assembly size. Unlike Nannochloropsis and Microchloropsis, long interspersed nuclear element repeats are abundant in Monodopsis and Vischeria and might constitute the centromeric regions. We found that both mevalonate and nonmevalonate pathways for terpenoid biosynthesis are present in Monodopsis and Vischeria, which is different from Nannochloropsis and Microchloropsis that have only the latter. Our analysis further revealed extensive spliced leader trans-splicing in Monodopsis and Vischeria at 36-61% of genes. Altogether, the high-quality genomes of Monodopsis and Vischeria not only serve as the much-needed outgroups to advance Nannochloropsis and Microchloropsis research, but also shed new light on the biology and evolution of eustigmatophyte algae.

Single-cell transcriptomic analysis of bloodstream Trypanosoma brucei reconstructs cell cycle progression and developmental quorum sensing

Developmental steps in the trypanosome life-cycle involve transition between replicative and non-replicative forms specialised for survival in, and transmission between, mammalian and tsetse fly hosts. Here, using oligopeptide-induced differentiation in vitro, we model the progressive development of replicative 'slender' to transmissible 'stumpy' bloodstream form Trypanosoma brucei and capture the transcriptomes of 8,599 parasites using single cell transcriptomics (scRNA-seq). Using this framework, we detail the relative order of biological events during asynchronous development, profile dynamic gene expression patterns and identify putative regulators. We additionally map the cell cycle of proliferating parasites and position stumpy cell-cycle exit at early G1 before progression to a distinct G0 state. A null mutant for one transiently elevated developmental regulator, ZC3H20 is further analysed by scRNA-seq, identifying its point of failure in the developmental atlas. This approach provides a paradigm for the dissection of differentiation events in parasites, relevant to diverse transitions in pathogen biology.

Peculiarities of aminoacyl-tRNA synthetases from trypanosomatids

Trypanosomatid parasites are responsible for various human diseases, such as sleeping sickness, animal trypanosomiasis, or cutaneous and visceral leishmaniases. The few available drugs to fight related parasitic infections are often toxic and present poor efficiency and specificity, and thus, finding new molecular targets is imperative. Aminoacyl-tRNA synthetases (aaRSs) are essential components of the translational machinery as they catalyze the specific attachment of an amino acid onto cognate tRNA(s). In trypanosomatids, one gene encodes both cytosolic- and mitochondrial-targeted aaRSs, with only three exceptions. We identify here a unique specific feature of aaRSs from trypanosomatids, which is that most of them harbor distinct insertion and/or extension sequences. Among the 26 identified aaRSs in the trypanosome Leishmania tarentolae, 14 contain an additional domain or a terminal extension, confirmed in mature mRNAs by direct cDNA nanopore sequencing. Moreover, these RNA-Seq data led us to address the question of aaRS dual localization and to determine splice-site locations and the 5'-UTR lengths for each mature aaRS-encoding mRNA. Altogether, our results provided evidence for at least one specific mechanism responsible for mitochondrial addressing of some L. tarentolae aaRSs. We propose that these newly identified features of trypanosomatid aaRSs could be developed as relevant drug targets to combat the diseases caused by these parasites.

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.

Evaluation of various RNA-seq approaches for identification of gene outrons in the flatworm Opisthorchis felineus

The parasitic flatworm Opisthorchis felineus is one of the causative agents of opisthorchiasis in humans. Recently, we assembled the O. felineus genome, but the correct genome annotation by means of standard methods was hampered by the presence of spliced leader trans-splicing (SLTS). As a result of SLTS, the original 5’-end (outron) of the transcripts is replaced by a short spliced leader sequence donated from a specialized SL RNA. SLTS is involved in the RNA processing of more than half of O. felineus genes, making it hard to determine the structure of outrons and bona fide transcription start sites of the corresponding genes and operons, being based solely on mRNA-seq data. In the current study, we tested various experimental approaches for identifying the sequences of outrons in O. felineus using massive parallel sequencing. Two of them were developed by us for targeted sequencing of already processed branched outrons. One was based on sequence-specific reverse transcription from the SL intron toward the 5’-end of the Y-branched outron. The other used outron hybridization with an immobilized single-stranded DNA probe complementary to the SL intron. Additionally, two approaches to the sequencing of rRNA-depleted total RNA were used, allowing the identification of a wider range of transcripts compared to mRNAseq. One is based on the enzymatic elimination of overrepresented cDNAs, the other utilizes exonucleolytic degradation of uncapped RNA by Terminator enzyme. By using the outron-targeting methods, we were not able to obtain the enrichment of RNA preparations by processed outrons, which is most likely indicative of a rapid turnover of these trans-splicing intermediate products. Of the two rRNA depletion methods, a method based on the enzymatic normalization of cDNA (Zymo-Seq RiboFree) showed high efficiency. Compared to mRNA-seq, it provides an approximately twofold increase in the fraction of reads originating from outrons and introns. The results suggest that unprocessed nascent transcripts are the main source of outron sequences in the RNA pool of O. felineus.

RNA structures in alternative splicing and back-splicing

Alternative splicing greatly expands the transcriptomic and proteomic diversities related to physiological and developmental processes in higher eukaryotes. Splicing of long noncoding RNAs, and back- and trans- splicing further expanded the regulatory repertoire of alternative splicing. RNA structures were shown to play an important role in regulating alternative splicing and back-splicing. Application of novel sequencing technologies made it possible to identify genome-wide RNA structures and interaction networks, which might provide new insights into RNA splicing regulation in vitro to in vivo. The emerging transcription-folding-splicing paradigm is changing our understanding of RNA alternative splicing regulation. Here, we review the insights into the roles and mechanisms of RNA structures in alternative splicing and back-splicing, as well as how disruption of these structures affects alternative splicing and then leads to human diseases. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

SLFinder, a pipeline for the novel identification of splice-leader sequences: a good enough solution for a complex problem

BACKGROUND

Spliced Leader trans-splicing is an important mechanism for the maturation of mRNAs in several lineages of eukaryotes, including several groups of parasites of great medical and economic importance. Nevertheless, its study across the tree of life is severely hindered by the problem of identifying the SL sequences that are being trans-spliced.

RESULTS

In this paper we present SLFinder, a four-step pipeline meant to identify de novo candidate SL sequences making very few assumptions regarding the SL sequence properties. The pipeline takes transcriptomic de novo assemblies and a reference genome as input and allows the user intervention on several points to account for unexpected features of the dataset. The strategy and its implementation were tested on real RNAseq data from species with and without SL Trans-Splicing.

CONCLUSIONS

SLFinder is capable to identify SL candidates with good precision in a reasonable amount of time. It is especially suitable for species with unknown SL sequences, generating candidate sequences for further refining and experimental validation.

Identification of clinically approved small molecules that inhibit growth and affect transcript levels of developmentally regulated genes in the African trypanosome

Trypanosoma brucei are unicellular parasites endemic to Sub-Saharan Africa that cause fatal disease in humans and animals. Infection with these parasites is caused by the bite of the tsetse fly vector, and parasites living extracellularly in the blood of infected animals evade the host immune system through antigenic variation. Existing drugs for Human and Animal African Trypanosomiasis are difficult to administer and can have serious side effects. Resistance to some drugs is also increasing, creating an urgent need for alternative trypanosomiasis therapeutics. We screened a library of 1,585 U.S. or foreign-approved drugs and identified 154 compounds that inhibit trypanosome growth. As all of these compounds have already undergone testing for human toxicity, they represent good candidates for repurposing as trypanosome therapeutics. In addition to identifying drugs that inhibit trypanosome growth, we wished to identify small molecules that can induce bloodstream form parasites to differentiate into forms adapted for the insect vector. These insect stage parasites lack the immune evasion mechanisms prevalent in bloodstream forms, making them vulnerable to the host immune system. To identify drugs that increase transcript levels of an invariant, insect-stage specific surface protein called procyclin, we engineered bloodstream reporter parasites that express Green Fluorescent Protein (GFP) following induction or stabilization of the procyclin transcript. Using these bloodstream reporter strains in combination with automated flow cytometry, we identified eflornithine, spironolactone, and phenothiazine as small molecules that increase abundance of procyclin transcript. Both eflornithine and spironolactone also affect transcript levels for a subset of differentiation associated genes. While we failed to identify compounds that increase levels of procyclin protein on the cell surface, this study is proof of principle that these fluorescent reporter parasites represent a useful tool for future small molecule or genetic screens aimed at identifying molecules or processes that initiate remodeling of the parasite surface during life cycle stage transitions.

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.

Unexpected diversity in eukaryotic transcription revealed by the retrotransposon hotspot family of Trypanosoma brucei

The path from DNA to RNA to protein in eukaryotes is guided by a series of factors linking transcription, mRNA export and translation. Many of these are conserved from yeast to humans. Trypanosomatids, which diverged early in the eukaryotic lineage, exhibit unusual features such as polycistronic transcription and trans-splicing of all messenger RNAs. They possess basal transcription factors, but lack recognisable orthologues of many factors required for transcription elongation and mRNA export. We show that retrotransposon hotspot (RHS) proteins fulfil some of these functions and that their depletion globally impairs nascent RNA synthesis by RNA polymerase II. Three sub-families are part of a coordinated process in which RHS6 is most closely associated with chromatin, RHS4 is part of the Pol II complex and RHS2 connects transcription with the translation machinery. In summary, our results show that the components of eukaryotic transcription are far from being universal, and reveal unsuspected plasticity in the course of evolution.

Fugacium Spliced Leader Genes Identified from Stranded RNA-Seq Datasets

Trans-splicing mechanisms have been documented in many lineages that are widely distributed phylogenetically, including dinoflagellates. The spliced leader (SL) sequence itself is conserved in dinoflagellates, although its gene sequences and arrangements have diversified within or across different species. In this study, we present 18 Fugacium kawagutii SL genes identified from stranded RNA-seq reads. These genes typically have a single SL but can contain several partial SLs with lengths ranging from 103 to 292 bp. Unexpectedly, we find the SL gene transcripts contain sequences upstream of the canonical SL, suggesting that generation of mature transcripts will require additional modifications following trans-splicing. We have also identified 13 SL-like genes whose expression levels and length are comparable to Dino-SL genes. Lastly, introns in these genes were identified and a new site for Sm-protein binding was proposed. Overall, this study provides a strategy for fast identification of SL genes and identifies new sequences of F. kawagutii SL genes to supplement our understanding of trans-splicing.

The Trypanosoma brucei RNA-Binding Protein TbRRM1 is Involved in the Transcription of a Subset of RNA Pol II-Dependent Genes

It has been long thought that RNA Polymerase (Pol) II transcriptional regulation does not operate in trypanosomes. However, recent reports have suggested that these organisms could regulate RNA Pol II transcription by epigenetic mechanisms. In this paper, we investigated the role of TbRRM1 in transcriptional regulation of RNA Pol II-dependent genes by focusing both in genes located in a particular polycistronic transcription unit (PTU) and in the monocistronic units of the SL-RNA genes. We showed that TbRRM1 is recruited throughout the PTU, with a higher presence on genes than intergenic regions. However, its depletion leads both to the decrease of nascent RNA and to chromatin compaction only of regions located distal to the main transcription start site. These findings suggest that TbRRM1 facilitates the RNA Pol II transcriptional elongation step by collaborating to maintain an open chromatin state in particular regions of the genome. Interestingly, the SL-RNA genes do not recruit TbRRM1 and, after TbRRM1 knockdown, nascent SL-RNAs accumulate while the chromatin state of these regions remains unchanged. Although it was previously suggested that TbRRM1 could regulate RNA Pol II-driven genes, we provide here the first experimental evidence which involves TbRRM1 to transcriptional regulation.

Heterodera glycines utilizes promiscuous spliced leaders and demonstrates a unique preference for a species-specific spliced leader over C. elegans SL1

Spliced leader trans-splicing (SLTS) plays a part in the maturation of pre-mRNAs in select species across multiple phyla but is particularly prevalent in Nematoda. The role of spliced leaders (SL) within the cell is unclear and an accurate assessment of SL occurrence within an organism is possible only after extensive sequencing data are available, which is not currently the case for many nematode species. SL discovery is further complicated by an absence of SL sequences from high-throughput sequencing results due to incomplete sequencing of the 5'-ends of transcripts during RNA-seq library preparation, known as 5'-bias. Existing datasets and novel methodology were used to identify both conserved SLs and unique hypervariable SLs within Heterodera glycines, the soybean cyst nematode. In H. glycines, twenty-one distinct SL sequences were found on 2,532 unique H. glycines transcripts. The SL sequences identified on the H. glycines transcripts demonstrated a high level of promiscuity, meaning that some transcripts produced as many as nine different individual SL-transcript combinations. Most uniquely, transcriptome analysis revealed that H. glycines is the first nematode to demonstrate a higher SL trans-splicing rate using a species-specific SL over well-conserved Caenorhabditis elegans SL-like sequences.

In-depth analysis of the genome of Trypanosoma evansi, an etiologic agent of surra

Trypanosoma evansi is the causative agent of the animal trypanosomiasis surra, a disease with serious economic burden worldwide. The availability of the genome of its closely related parasite Trypanosoma brucei allows us to compare their genetic and evolutionarily shared and distinct biological features. The complete genomic sequence of the T. evansi YNB strain was obtained using a combination of genomic and transcriptomic sequencing, de novo assembly, and bioinformatic analysis. The genome size of the T. evansi YNB strain was 35.2 Mb, showing 96.59% similarity in sequence and 88.97% in scaffold alignment with T. brucei. A total of 8,617 protein-coding genes, accounting for 31% of the genome, were predicted. Approximately 1,641 alternative splicing events of 820 genes were identified, with a majority mediated by intron retention, which represented a major difference in post-transcriptional regulation between T. evansi and T. brucei. Disparities in gene copy number of the variant surface glycoprotein, expression site-associated genes, microRNAs, and RNA-binding protein were clearly observed between the two parasites. The results revealed the genomic determinants of T. evansi, which encoded specific biological characteristics that distinguished them from other related trypanosome species.

Nuclear Phosphatidylinositol 5-Phosphatase Is Essential for Allelic Exclusion of Variant Surface Glycoprotein Genes in Trypanosomes

Allelic exclusion of variant surface glycoprotein (VSG) genes is essential for African trypanosomes to evade the host antibody response by antigenic variation. The mechanisms by which this parasite expresses only one of its ∼2,000 VSG genes at a time are unknown. We show that nuclear phosphatidylinositol 5-phosphatase (PIP5Pase) interacts with repressor activator protein 1 (RAP1) in a multiprotein complex and functions in the control of VSG allelic exclusion. RAP1 binds PIP5Pase substrate phosphatidylinositol 3,4,5-trisphosphate [PI(3,4,5)P3], and catalytic mutation of PIP5Pase that inhibits PI(3,4,5)P3 dephosphorylation results in simultaneous transcription of VSGs from all telomeric expression sites (ESs) and from silent subtelomeric VSG arrays. PIP5Pase and RAP1 bind to telomeric ESs, especially at 70-bp repeats and telomeres, and their binding is altered by PIP5Pase inactivation or knockdown, implying changes in ES chromatin organization. Our data suggest a model whereby PIP5Pase controls PI(3,4,5)P3 binding by RAP1 and, thus, RAP1 silencing of telomeric and subtelomeric VSG genes. Hence, allelic exclusion of VSG genes may entail control of nuclear phosphoinositides.

Quantitative RNA Analysis Using RNA-Seq

High-throughput sequencing of cDNA copies of mRNA (RNA-seq) provides a digital readout of mRNA levels over several orders of magnitude, as well as mapping the transcripts to the nucleotide level. Here we describe two different RNA-seq approaches, including one that exploits the 39-nucleotide mini-exon or spliced leader (SL) sequence found at the 5' end of all Leishmania (and other trypanosomatid) mRNAs.

Transcriptome Studies in Trypanosoma cruzi Using RNA-seq

RNA-seq is a powerful method to study in detail transcriptome changes in defined conditions, providing enormous amount of information on RNA stability and gene regulation. In this chapter, we describe a directional and a nondirectional library preparation protocol for RNA-seq in Trypanosoma cruzi, as well as a pipeline for bioinformatic analysis, which includes read trimming, alignment to a reference genome, and differential expression.

Glucose Signaling Is Important for Nutrient Adaptation during Differentiation of Pleomorphic African Trypanosomes

The African trypanosome has evolved mechanisms to adapt to changes in nutrient availability that occur during its life cycle. During transition from mammalian blood to insect vector gut, parasites experience a rapid reduction in environmental glucose. Here we describe how pleomorphic parasites respond to glucose depletion with a focus on parasite changes in energy metabolism and growth. Long slender bloodstream form parasites were rapidly killed as glucose concentrations fell, while short stumpy bloodstream form parasites persisted to differentiate into the insect-stage procyclic form parasite. The rate of differentiation was lower than that triggered by other cues but reached physiological rates when combined with cold shock. Both differentiation and growth of resulting procyclic form parasites were inhibited by glucose and nonmetabolizable glucose analogs, and these parasites were found to have upregulated amino acid metabolic pathway component gene expression. In summary, glucose transitions from the primary metabolite of the blood-stage infection to a negative regulator of cell development and growth in the insect vector, suggesting that the hexose is not only a key metabolic agent but also an important signaling molecule.IMPORTANCE As the African trypanosome Trypanosoma brucei completes its life cycle, it encounters many different environments. Adaptation to these environments includes modulation of metabolic pathways to parallel the availability of nutrients. Here, we describe how the blood-dwelling life cycle stages of the African trypanosome, which consume glucose to meet their nutritional needs, respond differently to culture in the near absence of glucose. The proliferative long slender parasites rapidly die, while the nondividing short stumpy parasite remains viable and undergoes differentiation to the next life cycle stage, the procyclic form parasite. Interestingly, a sugar analog that cannot be used as an energy source inhibited the process. Furthermore, the growth of procyclic form parasite that resulted from the event was inhibited by glucose, a behavior that is similar to that of parasites isolated from tsetse flies. Our findings suggest that glucose sensing serves as an important modulator of nutrient adaptation in the parasite.

Life cycle adapted upstream open reading frames (uORFs) in Trypanosoma congolense: A post-transcriptional approach to accurate gene regulation

The presented work explores the regulatory influence of upstream open reading frames (uORFs) on gene expression in Trypanosoma congolense. More than 31,000 uORFs in total were identified and characterized here. We found evidence for the uORFs’ appearance in the transcriptome to be correlated with proteomic expression data, clearly indicating their repressive potential in T. congolense, which has to rely on post-transcriptional gene expression regulation due to its unique genomic organization. Our data show that uORF’s translation repressive potential does not only correlate with elemental sequence features such as length, position and quantity, but involves more subtle components, in particular the codon and amino acid profiles. This corresponds with the popular mechanistic model of a ribosome shedding initiation factors during the translation of a uORF, which can prevent reinitiation at the downstream start codon of the actual protein-coding sequence, due to the former extensive consumption of crucial translation components. We suggest that uORFs with uncommon codon and amino acid usage can slow down the translation elongation process in T. congolense, systematically deplete the limited factors, and restrict downstream reinitiation, setting up a bottleneck for subsequent translation of the protein-coding sequence. Additionally we conclude that uORFs dynamically influence the T. congolense life cycle. We found evidence that transition to epimastigote form could be supported by gain of uORFs due to alternative trans-splicing, which down-regulate housekeeping genes’ expression and render the trypanosome in a metabolically reduced state of endurance.

Inhibition of maternal embryonic leucine zipper kinase with OTSSP167 displays potent anti-leukemic effects in chronic lymphocytic leukemia

TP53 pathway defects contributed to therapy resistance and adverse clinical outcome in chronic lymphocytic leukemia (CLL), which represents an unmet clinical need with few therapeutic options. Maternal embryonic leucine zipper kinase (MELK) is a novel oncogene, which plays crucial roles in mitotic progression and stem cell maintenance. OTSSP167, an orally administrated inhibitor targeting MELK, is currently in a phase I/II clinical trial in patients with advanced breast cancer and acute myeloid leukemia. Yet, no investigation has been elucidated to date regarding the oncogenic role of MELK and effects of OTSSP167 in chronic lymphocytic leukemia (CLL). Previous studies confirmed MELK inhibition abrogated cancer cell survival via p53 signaling pathway. Thus, we aimed to determine the biological function of MELK and therapeutic potential of OTSSP167 in CLL. Herein, MELK over-expression was observed in CLL cells, and correlated with higher WBC count, advanced stage, elevated LDH, increased β2-MG level, unmutated IGHV, positive ZAP-70, deletion of 17p13 and inferior prognosis of CLL patients. In accordance with functional enrichment analyses in gene expression profiling, CLL cells with depletion or inhibition of MELK exhibited impaired cell proliferation, enhanced fast-onset apoptosis, induced G2/M arrest, attenuated cell chemotaxis and promoted sensitivity to fludarabine and ibrutinib. However, gain-of-function assay showed increased cell proliferation and cell chemotaxis. In addition, OTSSP167 treatment reduced phosphorylation of AKT and ERK1/2. It decreased FoxM1 phosphorylation, expression of FoxM1, cyclin B1 and CDK1, while up-regulating p53 and p21 expression. Taken together, MELK served as a candidate of therapeutic target in CLL. OTSSP167 exhibits potent anti-tumor activities in CLL cells, highlighting a novel molecule-based strategy for leukemic interventions.

Identification of a novel fused gene family implicates convergent evolution in eukaryotic calcium signaling

Background

Both calcium signals and protein phosphorylation responses are universal signals in eukaryotic cell signaling. Currently three pathways have been characterized in different eukaryotes converting the Ca²⁺ signals to the protein phosphorylation responses. All these pathways have based mostly on studies in plants and animals.

Results

Based on the exploration of genomes and transcriptomes from all the six eukaryotic supergroups, we report here in Metakinetoplastina protists a novel gene family. This family, with a proposed name SCAMK, comprises SnRK3 fused calmodulin-like III kinase genes and was likely evolved through the insertion of a calmodulin-like3 gene into an SnRK3 gene by unequal crossover of homologous chromosomes in meiosis cell. Its origin dated back to the time intersection at least 450 million-year-ago when Excavata parasites, Vertebrata hosts, and Insecta vectors evolved. We also analyzed SCAMK’s unique expression pattern and structure, and proposed it as one of the leading calcium signal conversion pathways in Excavata parasite. These characters made SCAMK gene as a potential drug target for treating human African trypanosomiasis.

Conclusions

This report identified a novel gene fusion and dated its precise fusion time in Metakinetoplastina protists. This potential fourth eukaryotic calcium signal conversion pathway complements our current knowledge that convergent evolution occurs in eukaryotic calcium signaling.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4685-y) contains supplementary material, which is available to authorized users.

Landscape of the spliced leader trans-splicing mechanism in Schistosoma mansoni

Spliced leader dependent trans-splicing (SLTS) has been described as an important RNA regulatory process that occurs in different organisms, including the trematode Schistosoma mansoni. We identified more than seven thousand putative SLTS sites in the parasite, comprising genes with a wide spectrum of functional classes, which underlines the SLTS as a ubiquitous mechanism in the parasite. Also, SLTS gene expression levels span several orders of magnitude, showing that SLTS frequency is not determined by the expression level of the target gene, but by the presence of particular gene features facilitating or hindering the trans-splicing mechanism. Our in-depth investigation of SLTS events demonstrates widespread alternative trans-splicing (ATS) acceptor sites occurring in different regions along the entire gene body, highlighting another important role of SLTS generating alternative RNA isoforms in the parasite, besides the polycistron resolution. Particularly for introns where SLTS directly competes for the same acceptor substrate with cis-splicing, we identified for the first time additional and important features that might determine the type of splicing. Our study substantially extends the current knowledge of RNA processing by SLTS in S. mansoni, and provide basis for future studies on the trans-splicing mechanism in other eukaryotes.

The Pumilio-domain protein PUF6 contributes to SIDER2 retroposon-mediated mRNA decay in Leishmania

Leishmania and other trypanosomatid protozoa lack control at the level of transcription initiation and regulate gene expression exclusively post-transcriptionally. We have reported previously that Leishmania harbors a unique class of short interspersed degenerate retroposons (SIDERs) that are predominantly located within 3'UTRs and play a major role in post-transcriptional control. We have shown that members of the SIDER2 subfamily initiate mRNA decay through endonucleolytic cleavage within the second conserved 79-nt signature sequence of SIDER2 retroposons. Here, we have developed an optimized MS2 coat protein tethering system to capture trans-acting factor(s) regulating SIDER2-mediated mRNA decay. Tethering of the MS2 coat protein to a reporter RNA harboring two MS2 stem-loop aptamers and the cognate SIDER2-containing 3'UTR in combination with immunoprecipitation and mass spectrometry analysis led to the identification of RNA-binding proteins with known functions in mRNA decay. Among the candidate SIDER2-interacting proteins that were individually tethered to a SIDER2 reporter RNA, the Pumilio-domain protein PUF6 was shown to enhance degradation and reduce transcript half-life. Furthermore, we showed that PUF6 binds to SIDER2 sequences that include the regulatory 79-nt signature motif, hence contributing to the mRNA decay process. Consistent with a role of PUF6 in SIDER2-mediated decay, genetic inactivation of PUF6 resulted in increased accumulation and higher stability of endogenous SIDER2-bearing transcripts. Overall, these studies provide new insights into regulated mRNA decay pathways in Leishmania controlled by SIDER2 retroposons and propose a broader role for PUF proteins in mRNA decay within the eukaryotic kingdom.

Transporters of Trypanosoma brucei-phylogeny, physiology, pharmacology

Trypanosoma brucei comprise the causative agents of sleeping sickness, T. b. gambiense and T. b. rhodesiense, as well as the livestock-pathogenic T. b. brucei. The parasites are transmitted by the tsetse fly and occur exclusively in sub-Saharan Africa. T. brucei are not only lethal pathogens but have also become model organisms for molecular parasitology. We focus here on membrane transport proteins of T. brucei, their contribution to homeostasis and metabolism in the context of a parasitic lifestyle, and their pharmacological role as potential drug targets or routes of drug entry. Transporters and channels in the plasma membrane are attractive drug targets as they are accessible from the outside. Alternatively, they can be exploited to selectively deliver harmful substances into the trypanosome's interior. Both approaches require the targeted transporter to be essential: in the first case to kill the trypanosome, in the second case to prevent drug resistance due to loss of the transporter. By combining functional and phylogenetic analyses, we were mining the T. brucei predicted proteome for transporters of pharmacological significance. Here, we review recent progress in the identification of transporters of lipid precursors, amino acid permeases and ion channels in T. brucei.

Genomic and Transcriptomic Analysis Reveals Spliced Leader Trans-Splicing in Cryptomonads

Spliced leader trans-splicing (SLTS) is a poorly understood mechanism that is found in a diversity of eukaryotic lineages. In SLTS, a short RNA sequence is added near the 5′ ends of the transcripts of protein-coding genes by a modified spliceosomal reaction. Available data suggest that SLTS has evolved many times, and might be more likely to evolve in animals. That SLTS might be more likely to evolve in the context of the generally complex transcriptomes characteristic of animals suggests the possibility that SLTS functions in gene regulation or transcriptome diversification, however no general novel function for SLTS is known. Here, I report SLTS in a lineage of cellularly complex unicellular eukaryotes. Cryptomonads are a group of eukaryotic algae that acquired photosynthetic capacity by secondary endosymbiosis of a red alga, and that retain a reduced copy of the nucleus of the engulfed alga. I estimate that at least one-fifth of genes in the model cryptomonad Guillardia theta and its relative Hanusia phi undergo SLTS. I show that hundreds of genes in G. theta generate alternative transcripts by SLTS at alternative sites, however I find little evidence for alternative protein production by alternative SLTS splicing. Interestingly, I find no evidence for substantial operon structure in the G. theta genome, in contrast to previous findings in other lineages with SLTS. These results extend SLTS to another major group of eukaryotes, and heighten the mystery of the evolution of SLTS and its association with cellular and transcriptomic complexity.

Trypanosoma cruzi specific mRNA amplification by in vitro transcription improves parasite transcriptomics in host-parasite RNA mixtures

Background

Trypanosomatids are a group of protozoan parasites that includes the etiologic agents of important human illnesses as Chagas disease, sleeping sickness and leishmaniasis. These parasites have a significant distinction from other eukaryotes concerning mRNA structure, since all mature mRNAs have an identical species-specific sequence of 39 nucleotides at the 5′ extremity, named spliced leader (SL). Considering this peculiar aspect of trypanosomatid mRNA, the aim of the present work was to develop a Trypanosoma cruzi specific in vitro transcription (IVT) linear mRNA amplification method in order to improve parasite transcriptomics analyses.

Methods

We designed an oligonucleotide complementary to the last 21 bases of T. cruzi SL sequence, bearing an upstream T7 promoter (T7SL primer), which was used to direct the synthesis of second-strand cDNA. Original mRNA was then amplified by IVT using T7 RNA polymerase. T7SL-amplified RNA from two distinct T. cruzi stages (epimastigotes and trypomastigotes) were deep sequenced in SOLiD platform. Usual poly(A) + RNA and and T7-oligo(dT) amplified RNA (Eberwine method) were also sequenced. RNA-Seq reads were aligned to our new and improved T. cruzi Dm28c genome assembly (PacBio technology) and resulting transcriptome pattern from these three RNA preparation methods were compared, mainly concerning the conservation of mRNA transcritional levels and DEGs detection between epimastigotes and trypomastigotes.

Results

T7SL IVT method detected more potential differentially expressed genes in comparison to either poly(A) + RNA or T7dT IVT, and was also able to produce reliable quantifications of the parasite transcriptome down to 3 ng of total RNA. Furthermore, amplification of parasite mRNA in HeLa/epimastigote RNA mixtures showed that T7SL IVT generates transcriptome quantification with similar detection of differentially expressed genes when parasite RNA mass was only 0.1% of the total mixture (R = 0.78 when compared to poly(A) + RNA).

Conclusions

The T7SL IVT amplification method presented here allows the detection of more potential parasite differentially expressed genes (in comparison to poly(A) + RNA) in host-parasite mixtures or samples with low amount of RNA. This method is especially useful for trypanosomatid transcriptomics because it produces less bias than PCR-based mRNA amplification. Additionally, by simply changing the complementary region of the T7SL primer, the present method can be applied to any trypanosomatid species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4163-y) contains supplementary material, which is available to authorized users.

Variation in Basal Body Localisation and Targeting of Trypanosome RP2 and FOR20 Proteins

TOF-LisH-PLL motifs define FOP family proteins; some members are involved in flagellum assembly. The critical role of FOP family protein FOR20 is poorly understood. Here, we report relative localisations of the four FOP family proteins in parasitic Trypanosoma brucei: TbRP2, TbOFD1 and TbFOP/FOP1-like are mature basal body proteins whereas TbFOR20 is present on pro- and mature basal bodies - on the latter it localises distal to TbRP2. We discuss how the data, together with published work for another protist Giardia intestinalis, informs on likely FOR20 function. Moreover, our localisation study provides convincing evidence that the antigen recognised by monoclonal antibody YL1/2 at trypanosome mature basal bodies is FOP family protein TbRP2, not tyrosinated α-tubulin as widely stated in the literature. Curiously, FOR20 proteins from T. brucei and closely related African trypanosomes possess short, negatively-charged N-terminal extensions absent from FOR20 in other trypanosomatids and other eukaryotes. The extension is necessary for protein targeting, but insufficient to re-direct TbRP2 to probasal bodies. Yet, FOR20 from the American trypanosome T. cruzi, which lacks any extension, localises to pro- and mature basal bodies when expressed in T. brucei. This identifies unexpected variation in FOR20 architecture that is presently unique to one clade of trypanosomatids.

Expression of the RNA-binding protein RBP10 promotes the bloodstream-form differentiation state in Trypanosoma brucei

In nearly all eukaryotes, cellular differentiation is governed by changes in transcription, and stabilized by chromatin and DNA modification. Gene expression control in the pathogen Trypanosoma brucei, in contrast, relies almost exclusively on post-transcriptional mechanisms, so RNA binding proteins must assume the burden that is usually borne by transcription factors. T. brucei multiply in the blood of mammals as bloodstream forms, and in the midgut of Tsetse flies as procyclic forms. We show here that a single RNA-binding protein, RBP10, promotes the bloodstream-form trypanosome differentiation state. Depletion of RBP10 from bloodstream-form trypanosomes gives cells that can grow only as procyclic forms; conversely, expression of RBP10 in procyclic forms converts them to bloodstream forms. RBP10 binds to procyclic-specific mRNAs containing an UAUUUUUU motif, targeting them for translation repression and destruction. Products of RBP10 target mRNAs include not only the major procyclic surface protein and enzymes of energy metabolism, but also protein kinases and stage-specific RNA-binding proteins: this suggests that alterations in RBP10 trigger a regulatory cascade.

The proteome and transcriptome of the infectious metacyclic form of Trypanosoma brucei define quiescent cells primed for mammalian invasion

The infectious metacyclic forms of Trypanosoma brucei result from a complex development in the tsetse fly vector. When they infect mammals, they cause African sleeping sickness in humans. Due to scarcity of biological material and difficulties of the tsetse fly as an experimental system, very limited information is available concerning the gene expression profile of metacyclic forms. We used an in vitro system based on expressing the RNA binding protein 6 to obtain infectious metacyclics and determined their protein and mRNA repertoires by mass-spectrometry (MS) based proteomics and mRNA sequencing (RNA-Seq) in comparison to non-infectious procyclic trypanosomes. We showed that metacyclics are quiescent cells, and propose this influences the choice of a monocistronic variant surface glycoprotein expression site. Metacyclics have a largely bloodstream-form type transcriptome, and thus are programmed to translate a bloodstream-form type proteome upon entry into the mammalian host and resumption of cell division. Genes encoding cell surface components showed the largest changes between procyclics and metacyclics, observed at both the transcript and protein levels. Genes encoding metabolic enzymes exhibited expression in metacyclics with features of both procyclic and bloodstream forms, suggesting that this intermediate-type metabolism is dictated by the availability of nutrients in the tsetse fly vector.

GT-rich promoters can drive RNA pol II transcription and deposition of H2A.Z in African trypanosomes

Genome‐wide transcription studies are revealing an increasing number of “dispersed promoters” that, unlike “focused promoters”, lack well‐conserved sequence motifs and tight regulation. Dispersed promoters are nevertheless marked by well‐defined chromatin structures, suggesting that specific sequence elements must exist in these unregulated promoters. Here, we have analyzed regions of transcription initiation in the eukaryotic parasite Trypanosoma brucei, in which RNA polymerase II transcription initiation occurs over broad regions without distinct promoter motifs and lacks regulation. Using a combination of site‐specific and genome‐wide assays, we identified GT‐rich promoters that can drive transcription and promote the targeted deposition of the histone variant H2A.Z in a genomic context‐dependent manner. In addition, upon mapping nucleosome occupancy at high resolution, we find nucleosome positioning to correlate with RNA pol II enrichment and gene expression, pointing to a role in RNA maturation. Nucleosome positioning may thus represent a previously unrecognized layer of gene regulation in trypanosomes. Our findings show that even highly dispersed, unregulated promoters contain specific DNA elements that are able to induce transcription and changes in chromatin structure.

The Cytological Events and Molecular Control of Life Cycle Development of Trypanosoma brucei in the Mammalian Bloodstream

African trypanosomes cause devastating disease in sub-Saharan Africa in humans and livestock. The parasite lives extracellularly within the bloodstream of mammalian hosts and is transmitted by blood-feeding tsetse flies. In the blood, trypanosomes exhibit two developmental forms: the slender form and the stumpy form. The slender form proliferates in the bloodstream, establishes the parasite numbers and avoids host immunity through antigenic variation. The stumpy form, in contrast, is non-proliferative and is adapted for transmission. Here, we overview the features of slender and stumpy form parasites in terms of their cytological and molecular characteristics and discuss how these contribute to their distinct biological functions. Thereafter, we describe the technical developments that have enabled recent discoveries that uncover how the slender to stumpy transition is enacted in molecular terms. Finally, we highlight new understanding of how control of the balance between slender and stumpy form parasites interfaces with other components of the infection dynamic of trypanosomes in their mammalian hosts. This interplay between the host environment and the parasite’s developmental biology may expose new vulnerabilities to therapeutic attack or reveal where drug control may be thwarted by the biological complexity of the parasite’s lifestyle.

Multiplexed Spliced-Leader Sequencing: A high-throughput, selective method for RNA-seq in Trypanosomatids

High throughput sequencing techniques are poorly adapted for in vivo studies of parasites, which require prior in vitro culturing and purification. Trypanosomatids, a group of kinetoplastid protozoans, possess a distinctive feature in their transcriptional mechanism whereby a specific Spliced Leader (SL) sequence is added to the 5'end of each mRNA by trans-splicing. This allows to discriminate Trypansomatid RNA from mammalian RNA and forms the basis of our new multiplexed protocol for high-throughput, selective RNA-sequencing called SL-seq. We provided a proof-of-concept of SL-seq in Leishmania donovani, the main causative agent of visceral leishmaniasis in humans, and successfully applied the method to sequence Leishmania mRNA directly from infected macrophages and from highly diluted mixes with human RNA. mRNA profiles obtained with SL-seq corresponded largely to those obtained from conventional poly-A tail purification methods, indicating both enumerate the same mRNA pool. However, SL-seq offers additional advantages, including lower sequencing depth requirements, fast and simple library prep and high resolution splice site detection. SL-seq is therefore ideal for fast and massive parallel sequencing of parasite transcriptomes directly from host tissues. Since SLs are also present in Nematodes, Cnidaria and primitive chordates, this method could also have high potential for transcriptomics studies in other organisms.

Trypanosoma brucei TbIF1 inhibits the essential F1-ATPase in the infectious form of the parasite

The mitochondrial (mt) F_oF₁-ATP synthase of the digenetic parasite, Trypanosoma brucei, generates ATP during the insect procyclic form (PF), but becomes a perpetual consumer of ATP in the mammalian bloodstream form (BF), which lacks a canonical respiratory chain. This unconventional dependence on F_oF₁-ATPase is required to maintain the essential mt membrane potential (Δψm). Normally, ATP hydrolysis by this rotary molecular motor is restricted to when eukaryotic cells experience sporadic hypoxic conditions, during which this compulsory function quickly depletes the cellular ATP pool. To protect against this cellular treason, the highly conserved inhibitory factor 1 (IF1) binds the enzyme in a manner that solely inhibits the hydrolytic activity. Intriguingly, we were able to identify the IF1 homolog in T. brucei (TbIF1), but determined that its expression in the mitochondrion is tightly regulated throughout the life cycle as it is only detected in PF cells. TbIF1 appears to primarily function as an emergency brake in PF cells, where it prevented the restoration of the Δψm by F_oF₁-ATPase when respiration was chemically inhibited. In vitro, TbIF1 overexpression specifically inhibits the hydrolytic activity but not the synthetic capability of the F_oF₁-ATP synthase in PF mitochondria. Furthermore, low μM amounts of recombinant TbIF1 achieve the same inhibition of total mt ATPase activity as the F_oF₁-ATPase specific inhibitors, azide and oligomycin. Therefore, even minimal ectopic expression of TbIF1 in BF cells proved lethal as the indispensable Δψm collapsed due to inhibited F_oF₁-ATPase. In summary, we provide evidence that T. brucei harbors a natural and potent unidirectional inhibitor of the vital F_oF₁-ATPase activity that can be exploited for future structure-based drug design.

Enzymes are catalysts that drive both a forward and reverse chemical reaction depending on the thermodynamic properties. F_oF₁-ATP synthase is a multiprotein enzyme that under normal physiological conditions generates ATP. However, when respiration is impeded, this rotary molecular machine reverses and hydrolyzes ATP to pump protons and maintain the essential mitochondrial membrane potential. While this activity is exceptional in most eukaryotic cells, the unique composition of the Trypanosoma brucei mitochondrion dictates that the infectious stage of this human parasite is utterly dependent on the hydrolytic activity of F_oF₁-ATPase. While searching for better chemotherapeutics against Human African Trypanosomiasis, several trypanocidal compounds were determined to interact with this enzyme, but they indiscriminately inhibit both the ATP hydrolytic and synthetic activities. A more promising approach involves the conserved eukaryotic protein IF1, a unidirectional inhibitor that prevents just ATP hydrolysis. Auspiciously, we identified this protein homolog in T. brucei (TbIF1) and its expression is tightly regulated between life stages of the parasite. Importantly, the introduction of exogenous TbIF1 protein specifically inhibits F_oF₁-ATPase and is lethal for the infectious stage of T. brucei. Therefore, we have identified a natural inhibitor of an essential and druggable enzyme that can be exploited for future structure-based drug design.

The role of the PI(3,5)P2 kinase TbFab1 in endo/lysosomal trafficking in Trypanosoma brucei

Protein trafficking through endo/lysosomal compartments is critically important to the biology of the protozoan parasite Trypanosoma brucei, but the routes material may take to the lysosome, as well as the molecular factors regulating those routes, remain incompletely understood. Phosphoinositides are signaling phospholipids that regulate many trafficking events by recruiting specific effector proteins to discrete membrane subdomains. In this study, we investigate the role of one phosphoinositide, PI(3,5)P₂ in T. brucei. We find a low steady state level of PI(3,5)P₂ in bloodstream form parasites comparable to that of other organisms. RNAi knockdown of the putative PI(3)P-5 kinase TbFab1 decreases the PI(3,5)P₂ pool leading to rapid cell death. TbFab1 and PI(3,5)P₂ both localize strongly to late endo/lysosomes. While most trafficking functions were intact in TbFab1 deficient cells, including both endocytic and biosynthetic trafficking to the lysosome, lysosomal turnover of an endogenous ubiquitinylated membrane protein, ISG65, was completely blocked suggesting that TbFab1 plays a role in the ESCRT-mediated late endosomal/multivesicular body degradative pathways. Knockdown of a second component of PI(3,5)P₂ metabolism, the PI(3,5)P₂ phosphatase TbFig4, also resulted in delayed turnover of ISG65. Together, these results demonstrate an essential role for PI(3,5)P₂ in the turnover of ubiquitinylated membrane proteins and in trypanosome endomembrane biology.

Untranslated regions of mRNA and their role in regulation of gene expression in protozoan parasites

Protozoan parasites are one of the oldest living entities in this world that throughout their existence have shown excellent resilience to the odds of survival and have adapted beautifully to ever changing rigors of the environment. In view of the dynamic environment encountered by them throughout their life cycle, and in establishing pathogenesis, it is unsurprising that modulation of gene expression plays a fundamental role in their survival. In higher eukaryotes, untranslated regions (UTRs) of transcripts are one of the crucial regulators of gene expression (influencing mRNA stability and translation efficiency). Parasitic protozoan genome studies have led to the characterization (in silico, in vitro and in vivo) of a large number of their genes. Comparison of higher eukaryotic UTRs with parasitic protozoan UTRs reveals the existence of several similar and dissimilar facets of the UTRs. This review focuses on the elements of UTRs of medically important protozoan parasites and their regulatory role in gene expression. Such information may be useful to researchers in designing gene targeting strategies linked with perturbation of host-parasite relationships leading to control of specific parasites.