SLaP mapper: a webserver for identifying and quantifying spliced-leader addition and polyadenylation site usage in kinetoplastid genomes

Leishmania infantum (JPCM5) Transcriptome, Gene Models and Resources for an Active Curation of Gene Annotations

Leishmania infantum is one of the causative agents of visceral leishmaniases, the most severe form of leishmaniasis. An improved assembly for the L. infantum genome was published five years ago, yet delineation of its transcriptome remained to be accomplished. In this work, the transcriptome annotation was attained by a combination of both short and long RNA-seq reads. The good agreement between the results derived from both methodologies confirmed that transcript assembly based on Illumina RNA-seq and further delimitation according to the positions of spliced leader (SAS) and poly-A (PAS) addition sites is an adequate strategy to annotate the transcriptomes of Leishmania, a procedure previously used for transcriptome annotation in other Leishmania species and related trypanosomatids. These analyses also confirmed that the Leishmania transcripts boundaries are relatively slippery, showing extensive heterogeneity at the 5′- and 3′-ends. However, the use of RNA-seq reads derived from the PacBio technology (referred to as Iso-Seq) allowed the authors to uncover some complex transcription patterns occurring at particular loci that would be unnoticed by the use of short RNA-seq reads alone. Thus, Iso-Seq analysis provided evidence that transcript processing at particular loci would be more dynamic than expected. Another noticeable finding was the observation of a case of allelic heterozygosity based on the existence of chimeric Iso-Seq reads that might be generated by an event of intrachromosomal recombination. In addition, we are providing the L. infantum gene models, including both UTRs and CDS regions, that would be helpful for undertaking whole-genome expression studies. Moreover, we have built the foundations of a communal database for the active curation of both gene/transcript models and functional annotations for genes and proteins.

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.

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.

UTRme: A Scoring-Based Tool to Annotate Untranslated Regions in Trypanosomatid Genomes

Most signals involved in post-transcriptional regulatory networks are located in the untranslated regions (UTRs) of the mRNAs. Therefore, to deepen our understanding of gene expression regulation, delimitation of these regions with high accuracy is needed. The trypanosomatid lineage includes a variety of parasitic protozoans causing a significant worldwide burden on human health. Given their peculiar mechanisms of gene expression, these organisms depend on post-transcriptional regulation as the main level of gene expression control. In this context, the definition of the UTR regions becomes of key importance. We have developed UTR-mini-exon (UTRme), a graphical user interface (GUI) stand-alone application to identify and annotate 5′ and 3′ UTR regions in a highly accurate way. UTRme implements a multiple scoring system tailored to address the issue of false positive UTR assignment that frequently arise because of the characteristics of the intergenic regions. Even though it was developed for trypanosomatids, the tool can be used to predict 3′ sites in any eukaryote and 5′ UTRs in any organism where trans-splicing occurs (such as the model organism C. elegans). UTRme offers a way for non-bioinformaticians to precisely determine UTRs from transcriptomic data. The tool is freely available via the conda and github repositories.

An Alternative Strategy for Trypanosome Survival in the Mammalian Bloodstream Revealed through Genome and Transcriptome Analysis of the Ubiquitous Bovine Parasite Trypanosoma (Megatrypanum) theileri

There are hundreds of Trypanosoma species that live in the blood and tissue spaces of their vertebrate hosts. The vast majority of these do not have the ornate system of antigenic variation that has evolved in the small number of African trypanosome species, but can still maintain long-term infections in the face of the vertebrate adaptive immune system. Trypanosoma theileri is a typical example, has a restricted host range of cattle and other Bovinae, and is only occasionally reported to cause patent disease although no systematic survey of the effect of infection on agricultural productivity has been performed. Here, a detailed genome sequence and a transcriptome analysis of gene expression in bloodstream form T. theileri have been performed. Analysis of the genome sequence and expression showed that T. theileri has a typical kinetoplastid genome structure and allowed a prediction that it is capable of meiotic exchange, gene silencing via RNA interference and, potentially, density-dependent growth control. In particular, the transcriptome analysis has allowed a comparison of two distinct trypanosome cell surfaces, T. brucei and T. theileri, that have each evolved to enable the maintenance of a long-term extracellular infection in cattle. The T. theileri cell surface can be modeled to contain a mixture of proteins encoded by four novel large and divergent gene families and by members of a major surface protease gene family. This surface composition is distinct from the uniform variant surface glycoprotein coat on African trypanosomes providing an insight into a second mechanism used by trypanosome species that proliferate in an extracellular milieu in vertebrate hosts to avoid the adaptive immune response.

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.

Genome sequencing reveals metabolic and cellular interdependence in an amoeba-kinetoplastid symbiosis

Endosymbiotic relationships between eukaryotic and prokaryotic cells are common in nature. Endosymbioses between two eukaryotes are also known; cyanobacterium-derived plastids have spread horizontally when one eukaryote assimilated another. A unique instance of a non-photosynthetic, eukaryotic endosymbiont involves members of the genus Paramoeba, amoebozoans that infect marine animals such as farmed fish and sea urchins. Paramoeba species harbor endosymbionts belonging to the Kinetoplastea, a diverse group of flagellate protists including some that cause devastating diseases. To elucidate the nature of this eukaryote-eukaryote association, we sequenced the genomes and transcriptomes of Paramoeba pemaquidensis and its endosymbiont Perkinsela sp. The endosymbiont nuclear genome is ~9.5 Mbp in size, the smallest of a kinetoplastid thus far discovered. Genomic analyses show that Perkinsela sp. has lost the ability to make a flagellum but retains hallmark features of kinetoplastid biology, including polycistronic transcription, trans-splicing, and a glycosome-like organelle. Mosaic biochemical pathways suggest extensive ‘cross-talk’ between the two organisms, and electron microscopy shows that the endosymbiont ingests amoeba cytoplasm, a novel form of endosymbiont-host communication. Our data reveal the cell biological and biochemical basis of the obligate relationship between Perkinsela sp. and its amoeba host, and provide a foundation for understanding pathogenicity determinants in economically important Paramoeba.

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.

Nuclear Compartmentalization Contributes to Stage-Specific Gene Expression Control in Trypanosoma cruzi

In the protozoan parasite Trypanosoma cruzi, as in other trypanosomatids, transcription of protein coding genes occurs in a constitutive fashion, producing large polycistronic transcription units. These units are composed of non-functionally related genes which are pervasively processed to yield each mRNA. Therefore, post-transcriptional processes are crucial to regulate gene expression. Considering that nuclear compartmentalization could contribute to gene expression regulation, we comparatively studied the nuclear, cytoplasmic and whole cell transcriptomes of the non-infective epimastigote stage of T. cruzi, using RNA-Seq. We found that the cytoplasmic transcriptome tightly correlates with the whole cell transcriptome and both equally correlate with the proteome. Nonetheless, 1,200 transcripts showed differential abundance between the nuclear and cytoplasmic fractions. For the genes with transcript content augmented in the nucleus, significant structural and compositional differences were found. The analysis of the reported epimastigote translatome and proteome, revealed scarce ribosome footprints and encoded proteins for them. Ontology analyses unveiled that many of these genes are distinctive of other parasite life-cycle stages. Finally, the relocalization of transcript abundance in the metacyclic trypomastigote infective stage was confirmed for specific genes. While gene expression is strongly dependent on transcript steady-state level, we here highlight the importance of the distribution of transcripts abundance between compartments in T. cruzi. Particularly, we show that nuclear compartmentation is playing an active role in the developmental stage determination preventing off-stage expression.

Genome Sequence of Phytomonas françai, a Cassava (Manihot esculenta) Latex Parasite

Here, we report the genome sequence of the cassava (Manihot esculenta) latex parasite Phytomonas françai P. françai infection is linked with the yield-loss disease "chochamento de raizes" (empty roots) in the Unha variety of cassava, a disease characterized by poor root development and chlorosis of the leaves.

Morpholino antisense oligo inhibits trans-splicing of pre-inositol 1,4,5-trisphosphate receptor mRNA of Trypanosoma cruzi and suppresses parasite growth and infectivity

Morpholino antisense oligos (MAOs) are used to investigate physiological gene function by inhibiting gene translation or construction of specific alternative splicing variants by blocking cis-splicing. MAOs are attractive drug candidates for viral- and bacterial-infectious disease therapy because of properties such as in vivo stability and specificity to target genes. Recently, we showed that phosphorothioate antisense oligos against Trypanosoma cruzi inositol 1,4,5-trisphosphate receptor (TcIP(3)R) mRNA inhibit the parasite host cell infection. In the present study, we identified the spliced leader (SL) acceptor of pre-TcIP(3)R mRNA and synthesized MAO, which inhibited trans-splicing of the transcript (MAO-1). MAO-1 was found to inhibit the addition of SL-RNA to pre-TcIP(3)R mRNA by real-time RT-PCR analysis. Treatment of the parasites with MAO-1 significantly impaired the growth and infectivity into host cells. These results indicate that MAO-1 is a potential novel drug for Chagas disease and that MAOs inhibiting trans-splicing can be used to investigate the physiology of trypanosomal genes leading to the development of novel drugs.

Comparative Life Cycle Transcriptomics Revises Leishmania mexicana Genome Annotation and Links a Chromosome Duplication with Parasitism of Vertebrates

Leishmania spp. are protozoan parasites that have two principal life cycle stages: the motile promastigote forms that live in the alimentary tract of the sandfly and the amastigote forms, which are adapted to survive and replicate in the harsh conditions of the phagolysosome of mammalian macrophages. Here, we used Illumina sequencing of poly-A selected RNA to characterise and compare the transcriptomes of L. mexicana promastigotes, axenic amastigotes and intracellular amastigotes. These data allowed the production of the first transcriptome evidence-based annotation of gene models for this species, including genome-wide mapping of trans-splice sites and poly-A addition sites. The revised genome annotation encompassed 9,169 protein-coding genes including 936 novel genes as well as modifications to previously existing gene models. Comparative analysis of gene expression across promastigote and amastigote forms revealed that 3,832 genes are differentially expressed between promastigotes and intracellular amastigotes. A large proportion of genes that were downregulated during differentiation to amastigotes were associated with the function of the motile flagellum. In contrast, those genes that were upregulated included cell surface proteins, transporters, peptidases and many uncharacterized genes, including 293 of the 936 novel genes. Genome-wide distribution analysis of the differentially expressed genes revealed that the tetraploid chromosome 30 is highly enriched for genes that were upregulated in amastigotes, providing the first evidence of a link between this whole chromosome duplication event and adaptation to the vertebrate host in this group. Peptide evidence for 42 proteins encoded by novel transcripts supports the idea of an as yet uncharacterised set of small proteins in Leishmania spp. with possible implications for host-pathogen interactions.