A computational investigation of kinetoplastid trans-splicing

Structural analysis of PpSP15 and PsSP9 sand fly salivary proteins designed with a self-cleavable linker as a live vaccine candidate against cutaneous leishmaniasis

Background

Leishmania parasites are deposited in the host through sand fly bites along with sand fly saliva. Therefore, salivary proteins are promising vaccine candidates for controlling leishmaniasis. Herein, two immunogenic salivary proteins, PpSP15 from Phlebotomus papatasi and PsSP9 from Phlebotomus sergenti, were selected as vaccine candidates to be delivered by live Leishmania tarentolae as vector. The stepwise in silico protocol advantaged in this study for multi-protein design in L. tarentolae is then described in detail.

Methods

All possible combinations of two salivary proteins, PpSP15 and PsSP9, with or without T2A peptide were designed at the mRNA and protein levels. Then, the best combination for the vaccine candidate was selected based on mRNA and protein stability along with peptide analysis.

Results

At the mRNA level, the most favored secondary structure was PpSP15-T2A-PsSP9. At the protein level, the refined three-dimensional models of all combinations were structurally valid; however, local quality estimation showed that the PpSp15-T2A-PsSP9 fusion had higher stability for each amino acid position, with low root-mean-square deviation (RMSD), compared with the original proteins. In silico evaluation confirmed the PpSP15-T2A-PsSP9 combination as a good Th1-polarizing candidate in terms of high IFN-γ production and low IL-10/TGF-β ratio in response to three consecutive immunizations. Potential protein expression was then confirmed by Western blotting.

Conclusions

The approach presented herein is among the first studies to have privileged protein homology modeling along with mRNA analysis for logical live vaccine design-coding multi-proteins.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05437-x.

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.

Evolution of Alternative Splicing in Eudicots

Alternative pre-mRNA splicing (AS) is prevalent in plants and is involved in many interactions between plants and environmental stresses. However, the patterns and underlying mechanisms of AS evolution in plants remain unclear. By analyzing the transcriptomes of four eudicot species, we revealed that the divergence of AS is largely due to the gains and losses of AS events among orthologous genes. Furthermore, based on a subset of AS, in which AS can be directly associated with specific transcripts, we found that AS that generates transcripts containing premature termination codons (PTC), are likely more conserved than those that generate non-PTC containing transcripts. This suggests that AS coupled with nonsense-mediated decay (NMD) might play an important role in affecting mRNA levels post-transcriptionally. To understand the mechanisms underlying the divergence of AS, we analyzed the key determinants of AS using a machine learning approach. We found that the presence/absence of alternative splice site (SS) within the junction, the distance between the authentic SS and the nearest alternative SS, the size of exon-exon junctions were the major determinants for both alternative 5' donor site and 3' acceptor site among the studied species, suggesting a relatively conserved AS mechanism. The comparative analysis further demonstrated that variations of the identified AS determinants significantly contributed to the AS divergence among closely related species in both Solanaceae and Brassicaceae taxa. Together, these results provide detailed insights into the evolution of AS in plants.

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.

Programmed Genome Rearrangements in the Ciliate Oxytricha

The ciliate Oxytricha is a microbial eukaryote with two genomes, one of which experiences extensive genome remodeling during development. Each round of conjugation initiates a cascade of events that construct a transcriptionally active somatic genome from a scrambled germline genome, with considerable help from both long and small noncoding RNAs. This process of genome remodeling entails massive DNA deletion and reshuffling of remaining DNA segments to form functional genes from their interrupted and scrambled germline precursors. The use of Oxytricha as a model system provides an opportunity to study an exaggerated form of programmed genome rearrangement. Furthermore, studying the mechanisms that maintain nuclear dimorphism and mediate genome rearrangement has demonstrated a surprising plasticity and diversity of noncoding RNA pathways, with new roles that go beyond conventional gene silencing. Another aspect of ciliate genetics is their unorthodox patterns of RNA-mediated, epigenetic inheritance that rival Mendelian inheritance. This review takes the reader through the key experiments in a model eukaryote that led to fundamental discoveries in RNA biology and pushes the biological limits of DNA processing.

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

•

A web-server for identification of spliced-leader and polyadenylation addition sites.

•

Fully automated site quantification and gene assignment.

•

Multiple species within the Kinetoplastida.

The Kinetoplastida are a diverse and globally distributed class of free-living and parasitic single-celled eukaryotes that collectively cause a significant burden on human health and welfare. In kinetoplastids individual genes do not have promoters, but rather all genes are arranged downstream of a small number of RNA polymerase II transcription initiation sites and are thus transcribed in polycistronic gene clusters. Production of individual mRNAs from this continuous transcript occurs co-transcriptionally by trans-splicing of a ∼39 nucleotide capped RNA and subsequent polyadenylation of the upstream mRNA. SLaP mapper (Spliced-Leader and Polyadenylation mapper) is a fully automated web-service for identification, quantitation and gene-assignment of both spliced-leader and polyadenylation addition sites in Kinetoplastid genomes. SLaP mapper only requires raw read data from paired-end Illumina RNAseq and performs all read processing, mapping, quality control, quantification, and analysis in a fully automated pipeline. To provide usage examples and estimates of the quantity of sequence data required we use RNAseq obtained from two different library preparations from both Trypanosoma brucei and Leishmania mexicana to show the number of expected reads that are obtained from each preparation type. SLaP mapper is an easy to use, platform independent webserver that is freely available for use at http://www.stevekellylab.com/software/slap. Example files are provided on the website.

The calmodulin intergenic spacer as molecular target for characterization of Leishmania species

Background

Human leishmaniasis is a neglected disease caused by parasites of the genus Leishmania. Clinical aspects of this disease can vary significantly, reflecting the wide range of parasites in the genus Leishmania. Knowing accurately the Leishmania species infecting humans is important for clinical case management and evaluation of epidemiological risk. Calmodulin is an essential gene in trypanosomatids that modulates the calcium metabolism in various cellular activities. Despite its strong conservation in trypanosomatids, it has been recently observed that its untranslated regions (UTR) diverge among species.

Methods

In this study we analyzed the sequences and the absolute dinucleotide frequency of the intergenic spacer of the calmodulin gene (containing both, 3′ and 5′UTR) in nine reference Leishmania species and ten clinical isolates obtained from patients with cutaneous leishmaniasis.

Results

We show that the short calmodulin intergenic spacers exhibit features that make them interesting for applications in molecular characterization and phylogenetic studies of Leishmania. Dendrograms based on sequence alignments and on the dinucleotide frequency indicate that this particular region of calmodulin gene might be useful for species typing between the Leishmania and Viannia subgenera.

Conclusions

Mutations and composition of the calmodulin intergenic spacer from Leishmania species might have taxonomic value as parameters to define if an isolate is identical to a certain species or belongs to one of the two current subgenera.

"Eco-omics": a review of the application of genomics, transcriptomics, and proteomics for the study of the ecology of harmful algae

The implementation of molecular techniques has been widely adopted throughout the life sciences except in the marine sciences. The latter trend is quickly being reversed as even more cutting-edge molecular platforms, referred to collectively as 'omics-related technologies, are being used in a number of laboratories that study various aspects of life in the marine environment. This review provides a brief overview of just a few representative studies that have used genomics, transcriptomics, or proteomics approaches to deepen our understanding, specifically, about the underlying molecular biology of harmful algae. The examples of the studies described here are particularly relevant in showing how the information gleaned from these technologies can uncover the genetic capacity of harmful algal bloom-forming species, can generate new hypotheses about mechanistic relationships that bridge gene-environment interactions, and can impinge on our understanding surrounding the ecology of these organisms.

Ab initio identification of novel regulatory elements in the genome of Trypanosoma brucei by Bayesian inference on sequence segmentation

BACKGROUND

The rapid increase in the availability of genome information has created considerable demand for both comparative and ab initio predictive bioinformatic analyses. The biology laid bare in the genomes of many organisms is often novel, presenting new challenges for bioinformatic interrogation. A paradigm for this is the collected genomes of the kinetoplastid parasites, a group which includes Trypanosoma brucei the causative agent of human African trypanosomiasis. These genomes, though outwardly simple in organisation and gene content, have historically challenged many theories for gene expression regulation in eukaryotes.

METHODOLOGY/PRINCIPLE FINDINGS

Here we utilise a Bayesian approach to identify local changes in nucleotide composition in the genome of T. brucei. We show that there are several elements which are found at the starts and ends of multicopy gene arrays and that there are compositional elements that are common to all intergenic regions. We also show that there is a composition-inversion element that occurs at the position of the trans-splice site.

CONCLUSIONS/SIGNIFICANCE

The nature of the elements discovered reinforces the hypothesis that context dependant RNA secondary structure has an important influence on gene expression regulation in Trypanosoma brucei.

Genomic organization of leishmania species

Leishmania is a protozoan parasite belonging to the family Trypanosomatidae, which is found among 88 different countries. The parasite lives as an amastigote in vertebrate macrophages and as a promastigote in the digestive tract of sand fly. It can be cultured in the laboratory using appropriate culture media. Although the sexual cycle of Leishmania has not been observed during the promastigote and amastigote stages, it has been reported by some researchers. Leishmania has eukaryotic cell organization. Cell culture is convenient and cost effective, and because posttranslational modifications are common processes in the cultured cells, the cells are used as hosts for preparing eukaryotic recombinant proteins for research. Several transcripts of rDNA in the Leishmania genome are suitable regions for conducting gene transfer. Old World Leishmania spp. has 36 chromosomes, while New World Leishmania spp. has 34 or 35 chromosomes. The genomic organization and parasitic characteristics have been investigated. Leishmania spp. has a unique genomic organization among eukaryotes; the genes do not have introns, and the chromosomes are smaller with larger numbers of genes confined to a smaller space within the nucleus. Leishmania spp. genes are organized on one or both DNA strands and are transcribed as polycistronic (prokaryotic-like) transcripts from undefined promoters. Regulation of gene expression in the members of Trypanosomatidae differs from that in other eukaryotes. The trans-splicing phenomenon is a necessary step for mRNA processing in lower eukaryotes and is observed in Leishmania spp. Another particular feature of RNA editing in Leishmania spp. is that mitochondrial genes encoding respiratory enzymes are edited and transcribed. This review will discuss the chromosomal and mitochondrial (kinetoplast) genomes of Leishmania spp. as well as the phenomenon of RNA editing in the kinetoplast genome.

The transcriptome of the human pathogen Trypanosoma brucei at single-nucleotide resolution

The genome of Trypanosoma brucei, the causative agent of African trypanosomiasis, was published five years ago, yet identification of all genes and their transcripts remains to be accomplished. Annotation is challenged by the organization of genes transcribed by RNA polymerase II (Pol II) into long unidirectional gene clusters with no knowledge of how transcription is initiated. Here we report a single-nucleotide resolution genomic map of the T. brucei transcriptome, adding 1,114 new transcripts, including 103 non-coding RNAs, confirming and correcting many of the annotated features and revealing an extensive heterogeneity of 5′ and 3′ ends. Some of the new transcripts encode polypeptides that are either conserved in T. cruzi and Leishmania major or were previously detected in mass spectrometry analyses. High-throughput RNA sequencing (RNA-Seq) was sensitive enough to detect transcripts at putative Pol II transcription initiation sites. Our results, as well as recent data from the literature, indicate that transcription initiation is not solely restricted to regions at the beginning of gene clusters, but may occur at internal sites. We also provide evidence that transcription at all putative initiation sites in T. brucei is bidirectional, a recently recognized fundamental property of eukaryotic promoters. Our results have implications for gene expression patterns in other important human pathogens with similar genome organization (Trypanosoma cruzi, Leishmania sp.) and revealed heterogeneity in pre-mRNA processing that could potentially contribute to the survival and success of the parasite population in the insect vector and the mammalian host.

Identifying genes essential for survival in the host is fundamental to unraveling the biology of human pathogens and understanding mechanisms of pathogenesis. The protozoan parasite Trypanosoma brucei causes devastating diseases in humans and animals in sub-Saharan Africa, and the publication in 2005 of the genome sequence provided the first glance at the coding potential of this organism. Although at present there is a catalogue of predicted protein coding genes, the challenge remains to identify all authentic genes, including their boundaries. We used next generation RNA sequencing (RNA-Seq) to map transcribed regions and RNA polymerase II transcription initiation sites on a genome-wide scale. This approach allowed us to improve and correct the current annotation, to reveal a widespread heterogeneity of RNA processing sites (trans-splicing and polyadenylation) and to estimate that most genes are expressed at levels corresponding to 1 to 10 mRNAs per cell. Our data indicate that different transcript forms representing the same gene are present stochastically within the mRNA population. This unanticipated scenario may contribute to determining gene expression landscapes to adapt to different environments in the parasite life cycle.

Transgenic, fluorescent Leishmania mexicana allow direct analysis of the proteome of intracellular amastigotes

Investigating the proteome of intracellular pathogens is often hampered by inadequate methodologies to purify the pathogen free of host cell material. This has also precluded direct proteome analysis of the intracellular, amastigote form of Leishmania spp., protozoan parasites that cause a spectrum of diseases that affect some 12 million patients worldwide. Here a method is presented that combines classic, isopycnic density centrifugation with fluorescent particle sorting for purification by exploiting transgenic, fluorescent parasites to allow direct proteome analysis of the purified organisms. By this approach the proteome of intracellular Leishmania mexicana amastigotes was compared with that of extracellular promastigotes that are transmitted by insect vectors. In total, 509 different proteins were identified by mass spectrometry and database search. This number corresponds to approximately 6% of gene products predicted from the reference genome of Leishmania major. Intracellular amastigotes synthesized significantly more proteins with basic pI and showed a greater abundance of enzymes of fatty acid catabolism, which may reflect their living in acidic habitats and metabolic adaptation to nutrient availability, respectively. Bioinformatics analyses of the genes corresponding to the protein data sets produced clear evidence for skewed codon usage and translational bias in these organisms. Moreover analysis of the subset of genes whose products were more abundant in amastigotes revealed characteristic sequence motifs in 3'-untranslated regions that have been linked to translational control elements. This suggests that proteome data sets may be used to identify regulatory elements in mRNAs. Last but not least, at 6% coverage the proteome identified all vaccine antigens tested to date. Thus, the present data set provides a valuable resource for selection of candidate vaccine antigens.

Improving the prediction of mRNA extremities in the parasitic protozoan Leishmania

Background

Leishmania and other members of the Trypanosomatidae family diverged early on in eukaryotic evolution and consequently display unique cellular properties. Their apparent lack of transcriptional regulation is compensated by complex post-transcriptional control mechanisms, including the processing of polycistronic transcripts by means of coupled trans-splicing and polyadenylation. Trans-splicing signals are often U-rich polypyrimidine (poly(Y)) tracts, which precede AG splice acceptor sites. However, as opposed to higher eukaryotes there is no consensus polyadenylation signal in trypanosomatid mRNAs.

Results

We refined a previously reported method to target 5' splice junctions by incorporating the pyrimidine content of query sequences into a scoring function. We also investigated a novel approach for predicting polyadenylation (poly(A)) sites in-silico, by comparing query sequences to polyadenylated expressed sequence tags (ESTs) using position-specific scanning matrices (PSSMs). An additional analysis of the distribution of putative splice junction to poly(A) distances helped to increase prediction rates by limiting the scanning range. These methods were able to simplify splice junction prediction without loss of precision and to increase polyadenylation site prediction from 22% to 47% within 100 nucleotides.

Conclusion

We propose a simplified trans-splicing prediction tool and a novel poly(A) prediction tool based on comparative sequence analysis. We discuss the impact of certain regions surrounding the poly(A) sites on prediction rates and contemplate correlating biological mechanisms. This work aims to sharpen the identification of potentially functional untranslated regions (UTRs) in a large-scale, comparative genomics framework.

Developmental regulation of gene expression in trypanosomatid parasitic protozoa

Kinetoplastids branched early from the eukaryotic lineage and include several parasitic protozoan species. Up to several hundred kinetoplastid genes are co-transcribed into polycistronic RNAs and individual mRNAs are resolved by coupled co-transcriptional trans-splicing of a universal splice-leader RNA (SL-RNA) and 3'-end maturation processes. Protein-coding genes lack RNA polymerase II promoters. Consequently, most of gene regulation in these organisms occurs post-transcriptionally. Over the last few years, many more genes that are regulated at the mRNA stability level and a few at the translation level have been reported. Almost all major trypanosome homologues of yeast/mammalian mRNA degradation enzymes have been functionally characterized and major pathways identified. Novel paradigms have also recently emerged: regulated post-transcriptional processing of cytoplasmic RNAs, SL-RNA transcriptional silencing-mediated global stress response, and Leishmania-specific large-scale modulation of post-transcriptional gene expression via inactive degenerated retroelements. Several of these developments have greatly benefited from the recently completed genomic sequences and functional genomic studies.

Codon usage suggests that translational selection has a major impact on protein expression in trypanosomatids

Background

Different proteins are required in widely different quantities to build a living cell. In most organisms, transcription control makes a major contribution to differential expression. This is not the case in trypanosomatids where most genes are transcribed at an equivalent rate within large polycistronic clusters. Thus, trypanosomatids must use post-transcriptional control mechanisms to balance gene expression requirements.

Results

Here, the evidence for translational selection, the enrichment of 'favoured' codons in more highly expressed genes, is explored. A set of highly expressed, tandem-repeated genes display codon bias in Trypanosoma cruzi, Trypanosoma brucei and Leishmania major. The tRNA complement reveals forty-five of the sixty-one possible anticodons indicating widespread use of 'wobble' tRNAs. Consistent with translational selection, cognate tRNA genes for favoured codons are over-represented. Importantly, codon usage (Codon Adaptation Index) correlates with predicted and observed expression level. In addition, relative codon bias is broadly conserved among syntenic genes from different trypanosomatids.

Conclusion

Synonymous codon bias is correlated with tRNA gene copy number and with protein expression level in trypanosomatids. Taken together, the results suggest that translational selection is the dominant mechanism underlying the control of differential protein expression in these organisms. The findings reveal how trypanosomatids may compensate for a paucity of canonical Pol II promoters and subsequent widespread constitutive RNA polymerase II transcription.

The bursicon gene in mosquitoes: an unusual example of mRNA trans-splicing

The bursicon gene in Anopheles gambiae is encoded by two loci. Burs124 on chromosome arm 2L contains exons 1, 2, and 4, while burs3 on arm 2R contains exon 3. Exon 3 is efficiently spliced into position in the mature transcript. This unusual gene arrangement is ancient within mosquitoes, being shared by Aedes aegypti and Culex pipiens.

The untranslated regions of genes from Trypanosoma cruzi: perspectives for functional characterization of strains and isolates

The sequencing of Trypanosoma cruzi genome has been completed and a great deal of information is now available. However, the organization of protozoa genomes is somewhat elusive and much effort must be applied to reveal all the information coded in the nucleotide sequences. Among the DNA segments that needs further investigation are the untranslated regions of genes. Many of the T. cruzi genes that were revealed by the genome sequencing lack information about the untranslated regions. In this paper, some features of these untranslated segments as well as their applications in T. cruzi populations are discussed.

Spliced leader RNA trans-splicing in dinoflagellates

Through the analysis of hundreds of full-length cDNAs from fifteen species representing all major orders of dinoflagellates, we demonstrate that nuclear-encoded mRNAs in all species, from ancestral to derived lineages, are trans-spliced with the addition of the 22-nt conserved spliced leader (SL), DCCGUAGCCAUUUUGGCUCAAG (D = U, A, or G), to the 5' end. SL trans-splicing has been documented in a limited but diverse number of eukaryotes, in which this process makes it possible to translate polycistronically transcribed nuclear genes. In SL trans-splicing, SL-donor transcripts (SL RNAs) contain two functional domains: an exon that provides the SL for mRNA and an intron that contains a spliceosomal (Sm) binding site. In dinoflagellates, SL RNAs are unusually short at 50-60 nt, with a conserved Sm binding motif (AUUUUGG) located in the SL (exon) rather than the intron. The initiation nucleotide is predominantly U or A, an unusual feature that may affect capping, and hence the translation and stability of the recipient mRNA. The core SL element was found in mRNAs coding for a diverse array of proteins. Among the transcripts characterized were three homologs of Sm-complex subunits, indicating that the role of the Sm binding site is conserved, even if the location on the SL is not. Because association with an Sm-complex often signals nuclear import for U-rich small nuclear RNAs, it is unclear how this Sm binding site remains on mature mRNAs without impeding cytosolic localization or translation of the latter. The sequences reported in this paper have been deposited in the GenBank database (accession nos. AF 512889, DQ 864761-DQ 864971, DQ 867053-DQ 867070, DQ 884413-DQ 884451, EF 133854-EF 133905, EF 133961-EF 134003, EF 134083-EF 134402, EF 141835, and EF 143070-EF 143105).

Mapping the antigenicity of the parasites in Leishmania donovani infection by proteome serology

Background

Leishmaniasis defines a cluster of protozoal diseases with diverse clinical manifestations. The visceral form caused by Leishmania donovani is the most severe. So far, no vaccines exist for visceral leishmaniasis despite indications of naturally developing immunity, and sensitive immunodiagnostics are still at early stages of development.

Methodology/Principle Findings

Establishing a proteome-serological methodology, we mapped the antigenicity of the parasites and the specificities of the immune responses in human leishmaniasis. Using 2-dimensional Western blot analyses with sera and parasites isolated from patients in India, we detected immune responses with widely divergent specificities for up to 330 different leishmanial antigens. 68 antigens were assigned to proteins in silver- and fluorochrome-stained gels. The antigenicity of these proteins did not correlate with the expression levels of the proteins. Although some antigens are shared among different parasite isolates, there are extensive differences and no immunodominant antigens, but indications of antigenic drift in the parasites. Six antigens were identified by mass spectrometry.

Conclusions/Significance

Proteomics-based dissection of the serospecificities of leishmaniasis patients provides a comprehensive inventory of the complexity and interindividual heterogeneity of the host-responses to and variations in the antigenicity of the Leishmania parasites. This information can be instrumental in the development of vaccines and new immune monitoring and diagnostic devices.

Correction: genetic algorithm learning as a robust approach to RNA editing site site prediction

After the publication of [1], we were alerted to an error in our data. The error was an one-off miscalculation in the extraction of position information for our set of true negatives. Our data set should have used randomly selected non-edited cytosines (C) as true negatives, but the data generation phase resulted in a set of nucleotides that were each one nucleotide downstream of known, unedited cytosines. The consequences of this error are reflected in changes to our results, although the general conclusions presented in our original publication remain largely unchanged.