Computational identification of four spliceosomal snRNAs from the deep-branching eukaryote Giardia intestinalis

Spliceosomal introns in the diplomonad parasite Giardia duodenalis revisited

Complete reference genomes, including correct feature annotations, are a fundamental aspect of genomic biology. In the case of protozoan species such as Giardia duodenalis, a major human and animal parasite worldwide, accurate genome annotation can deepen our understanding of the evolution of parasitism and pathogenicity by identifying genes underlying key traits and clinically relevant cellular mechanisms, and by extension, the development of improved prevention strategies and treatments. This study used bioinformatics analyses of Giardia mRNA libraries to characterize known introns and identify new intron candidates, working towards completion of the G. duodenalis assemblage A strain 'WB' genome and further elucidating Giardia's gene expression. By using a set of experimentally validated positive control loci to calibrate our intron detection pipeline, we were able to detect evidence of previously missed candidate splice junctions directly from expressed transcript data. These intron candidates were further studied in silico using NMDS (non-metric multidimensional scaling) clustering to determine shared characteristics and their relative importance such as secondary structure, splicing efficiency and motif conservation, and thus to refine intron models. Results from this study identified 34 new intron candidates, with several potential introns showing evidence that secondary structure of the mRNA molecule might play a more significant role in splicing than previously reported eukaryotic splicing activity mediated by a reduced spliceosome present in G. duodenalis.

An Unusual U2AF2 Inhibits Splicing and Attenuates the Virulence of the Human Protozoan Parasite Entamoeba histolytica

E. histolytica is the etiological agent of intestinal amebiasis and liver abscesses, which still poses public health threat globally. Metronidazole is the drug of choice against amebiasis. However, metronidazole-resistant amoebic clinical isolates and strains have been reported recently, challenging the efforts for amebiasis eradication. In search of alternative treatments, E. histolytica transcriptomes have shown the association of genes involved in RNA metabolism with the virulence of the parasite. Among the upregulated genes in amoebic liver abscesses are the splicing factors EhU2AF2 and a paralog of EhSF3B1. For this reason and because EhU2AF2 contains unusual KH-QUA2 (84KQ) motifs in its lengthened C-terminus domain, here we investigated how the role of EhU2AF2 in pre-mRNA processing impacts the virulence of the parasite. We found that 84KQ is involved in splicing inhibition/intron retention of several virulence and non-virulence-related genes. The 84KQ domain interacts with the same domain of the constitutive splicing factor SF1 (SF1KQ), both in solution and when SF1KQ is bound to branchpoint signal RNA probes. The 84KQ–SF1KQ interaction prevents splicing complex E to A transition, thus inhibiting splicing. Surprisingly, the deletion of the 84KQ domain in EhU2AF2 amoeba transformants increased splicing and enhanced the in vitro and in vivo virulence phenotypes. We conclude that the interaction of the 84KQ and SF1KQ domains, probably involving additional factors, tunes down Entamoeba virulence by favoring intron retention.

The Integrator complex regulates differential snRNA processing and fate of adult stem cells in the highly regenerative planarian Schmidtea mediterranea

In multicellular organisms, cell type diversity and fate depend on specific sets of transcript isoforms generated by post-transcriptional RNA processing. Here, we used Schmidtea mediterranea, a flatworm with extraordinary regenerative abilities and a large pool of adult stem cells, as an in vivo model to study the role of Uridyl-rich small nuclear RNAs (UsnRNAs), which participate in multiple RNA processing reactions including splicing, in stem cell regulation. We characterized the planarian UsnRNA repertoire, identified stem cell-enriched variants and obtained strong evidence for an increased rate of UsnRNA 3'-processing in stem cells compared to their differentiated counterparts. Consistently, components of the Integrator complex showed stem cell-enriched expression and their depletion by RNAi disrupted UsnRNA processing resulting in global changes of splicing patterns and reduced processing of histone mRNAs. Interestingly, loss of Integrator complex function disrupted both stem cell maintenance and regeneration of tissues. Our data show that the function of the Integrator complex in UsnRNA 3'-processing is conserved in planarians and essential for maintaining their stem cell pool. We propose that cell type-specific modulation of UsnRNA composition and maturation contributes to in vivo cell fate choices, such as stem cell self-renewal in planarians.

Sense-antisense gene overlap is probably a cause for retaining the few introns in Giardia genome and the implications

Background

It is widely accepted that the last eukaryotic common ancestor and early eukaryotes were intron-rich and intron loss dominated subsequent evolution, thus the presence of only very few introns in some modern eukaryotes must be the consequence of massive loss. But it is striking that few eukaryotes were found to have completely lost introns. Despite extensive research, the causes of massive intron losses remain elusive. Actually the reverse question -- how the few introns can be retained under the evolutionary selection pressure of intron loss -- is equally significant but was rarely studied, except that it was conjectured that the essential functions of some introns prevent their loss. The situation that extremely few (eight) spliceosome-mediated cis-spliced introns present in the relatively simple genome of Giardia lamblia provides an excellent opportunity to explore this question.

Results

Our investigation found three types of distribution patterns of the few introns in the intron-containing genes: ancient intron in ancient gene, later-evolved intron in ancient gene, and later-evolved intron in later-evolved gene, which can reflect to some extent the dynamic evolution of introns in Giardia. Without finding any special features or functional importance of these introns responsible for their retention, we noticed and experimentally verified that some intron-containing genes form sense-antisense gene pairs with transcribable genes on their complementary strands, and that the introns just reside in the overlapping regions.

Conclusions

In Giardia’s evolution, despite constant evolutionary selection pressure of intron loss, intron gain can still occur in both ancient and later-evolved genes, but only a few introns are retained; at least the evolutionary retention of some of the introns might not be due to the functional constraint of the introns themselves but the causes outside of introns, such as the constraints imposed by other genomic functional elements overlapping with the introns. These findings can not only provide some clues to find new genomic functional elements -- in the areas overlapping with introns, but suggest that “functional constraint” of introns may not be necessarily directly associated with intron loss and gain, and that the real functions are probably still outside of our current knowledge.

Reviewers

This article was reviewed by Mikhail Gelfand, Michael Gray, and Igor Rogozin.

Electronic supplementary material

The online version of this article (10.1186/s13062-018-0226-5) contains supplementary material, which is available to authorized users.

Cardinal Epigenetic Role of non-coding Regulatory RNAs in Circadian Rhythm

Circadian rhythm which governs basic physiological activities like sleeping, feeding and energy consumption is regulated by light-controlled central clock genes in the pacemaker neuron. The timekeeping machinery with unique transcriptional and post-transcriptional feedback loops is controlled by different small regulatory RNAs in the brain. Roles of the multiple neuronal genes, especially post-transcriptional regulation, splicing, polyadenylation, mature mRNA editing, and stability of translation products, are controlled by epigenetic activities orchestrated via small RNAs. Collectively, these mechanisms regulate clock and light-controlled genes for effecting pacemaker activity and entrainment. Regulatory small RNAs of the circadian circuit, timekeeping mechanism, synchronization of regular entrainment, oscillation, and rhythmicity are regulated by diversified RNA molecules. Regulatory small RNAs operate critical roles in brain activities including the neuronal clock activity. In this report, we propose the emergence of the earlier unexpected small RNAs for a historic perspective of epigenetic regulation of the brain clock system.

Reconstruction of Sugar Metabolic Pathways of Giardia lamblia

Giardia lamblia is an “important” pathogen of humans, but as a diplomonad excavate it is evolutionarily distant from other eukaryotes and relatively little is known about its core metabolic pathways. KEGG, the widely referenced site for providing information of metabolism, does not yet include many enzymes from Giardia species. Here we identify Giardia's core sugar metabolism using standard bioinformatic approaches. By comparing Giardia proteomes with known enzymes from other species, we have identified enzymes in the glycolysis pathway, as well as some enzymes involved in the TCA cycle and oxidative phosphorylation. However, the majority of enzymes from the latter two pathways were not identifiable, indicating the likely absence of these functionalities. We have also found enzymes from the Giardia glycolysis pathway that appear more similar to those from bacteria. Because these enzymes are different from those found in mammals, the host organisms for Giardia, we raise the possibility that these bacteria-like enzymes could be novel drug targets for treating Giardia infections.

Evolutionarily divergent spliceosomal snRNAs and a conserved non-coding RNA processing motif in Giardia lamblia

Non-coding RNAs (ncRNAs) have diverse essential biological functions in all organisms, and in eukaryotes, two such classes of ncRNAs are the small nucleolar (sno) and small nuclear (sn) RNAs. In this study, we have identified and characterized a collection of sno and snRNAs in Giardia lamblia, by exploiting our discovery of a conserved 12 nt RNA processing sequence motif found in the 3' end regions of a large number of G. lamblia ncRNA genes. RNA end mapping and other experiments indicate the motif serves to mediate ncRNA 3' end formation from mono- and di-cistronic RNA precursor transcripts. Remarkably, we find the motif is also utilized in the processing pathway of all four previously identified trans-spliced G. lamblia introns, revealing a common RNA processing pathway for ncRNAs and trans-spliced introns in this organism. Motif sequence conservation then allowed for the bioinformatic and experimental identification of additional G. lamblia ncRNAs, including new U1 and U6 spliceosomal snRNA candidates. The U6 snRNA candidate was then used as a tool to identity novel U2 and U4 snRNAs, based on predicted phylogenetically conserved snRNA-snRNA base-pairing interactions, from a set of previously identified G. lamblia ncRNAs without assigned function. The Giardia snRNAs retain the core features of spliceosomal snRNAs but are sufficiently evolutionarily divergent to explain the difficulties in their identification. Most intriguingly, all of these snRNAs show structural features diagnostic of U2-dependent/major and U12-dependent/minor spliceosomal snRNAs.

Characterizing ncRNAs in Human Pathogenic Protists Using High-Throughput Sequencing Technology

ncRNAs are key genes in many human diseases including cancer and viral infection, as well as providing critical functions in pathogenic organisms such as fungi, bacteria, viruses, and protists. Until now the identification and characterization of ncRNAs associated with disease has been slow or inaccurate requiring many years of testing to understand complicated RNA and protein gene relationships. High-throughput sequencing now offers the opportunity to characterize miRNAs, siRNAs, small nucleolar RNAs (snoRNAs), and long ncRNAs on a genomic scale, making it faster and easier to clarify how these ncRNAs contribute to the disease state. However, this technology is still relatively new, and ncRNA discovery is not an application of high priority for streamlined bioinformatics. Here we summarize background concepts and practical approaches for ncRNA analysis using high-throughput sequencing, and how it relates to understanding human disease. As a case study, we focus on the parasitic protists Giardia lamblia and Trichomonas vaginalis, where large evolutionary distance has meant difficulties in comparing ncRNAs with those from model eukaryotes. A combination of biological, computational, and sequencing approaches has enabled easier classification of ncRNA classes such as snoRNAs, but has also aided the identification of novel classes. It is hoped that a higher level of understanding of ncRNA expression and interaction may aid in the development of less harsh treatment for protist-based diseases.

TEMPERATURE-DEPENDENT STRUCTURAL VARIABILITY OF RNAs: SPLICED LEADER RNAs AND THEIR EVOLUTIONARY HISTORY

The structures attained by RNA molecules depend not only on their sequence but also on environmental parameters such as their temperature. So far, this effect has been largely neglected in bioinformatics studies. Here, we show that structural comparisons can be facilitated and more coherent structural models can be obtained when differences in environmental parameters are taken into account. We re-evaluate the secondary structures of the spliced leader (SL) RNAs from the seven eukaryotic phyla in which SL RNA trans-splicing has been described. Adjusting structure prediction to the natural growth temperatures and considering energetically similar secondary structures, we observe striking similarities among Euglenida, Kinetoplastida, Dinophyceae, Cnidaria, Rotifera, Nematoda, Platyhelminthes, and Tunicata that cannot be explained easily by the independent innovation of SL RNAs in each of these phyla. Supplementary Table is available at .

Characterization of RNase MRP RNA and novel snoRNAs from Giardia intestinalis and Trichomonas vaginalis

Background

Eukaryotic cells possess a complex network of RNA machineries which function in RNA-processing and cellular regulation which includes transcription, translation, silencing, editing and epigenetic control. Studies of model organisms have shown that many ncRNAs of the RNA-infrastructure are highly conserved, but little is known from non-model protists. In this study we have conducted a genome-scale survey of medium-length ncRNAs from the protozoan parasites Giardia intestinalis and Trichomonas vaginalis.

Results

We have identified the previously 'missing' Giardia RNase MRP RNA, which is a key ribozyme involved in pre-rRNA processing. We have also uncovered 18 new H/ACA box snoRNAs, expanding our knowledge of the H/ACA family of snoRNAs.

Conclusions

Results indicate that Giardia intestinalis and Trichomonas vaginalis, like their distant multicellular relatives, contain a rich infrastructure of RNA-based processing. From here we can investigate the evolution of RNA processing networks in eukaryotes.

The short non-coding transcriptome of the protozoan parasite Trypanosoma cruzi

The pathway for RNA interference is widespread in metazoans and participates in numerous cellular tasks, from gene silencing to chromatin remodeling and protection against retrotransposition. The unicellular eukaryote Trypanosoma cruzi is missing the canonical RNAi pathway and is unable to induce RNAi-related processes. To further understand alternative RNA pathways operating in this organism, we have performed deep sequencing and genome-wide analyses of a size-fractioned cDNA library (16–61 nt) from the epimastigote life stage. Deep sequencing generated 582,243 short sequences of which 91% could be aligned with the genome sequence. About 95–98% of the aligned data (depending on the haplotype) corresponded to small RNAs derived from tRNAs, rRNAs, snRNAs and snoRNAs. The largest class consisted of tRNA-derived small RNAs which primarily originated from the 3′ end of tRNAs, followed by small RNAs derived from rRNA. The remaining sequences revealed the presence of 92 novel transcribed loci, of which 79 did not show homology to known RNA classes.

Chagas' disease is a major health problem in Latin America and is caused by the protozoan parasite Trypanosoma cruzi. T. cruzi lacks the pathway for RNA interference, which is widespread among eukaryotes, and is therefore unable to induce RNAi-related processes. In many organisms, small RNAs play an important role in regulating gene expression and other cellular processes. In order to understand if other small RNA pathways are operating in this organism, we performed high throughput sequencing and genome-wide analyses of the short transcriptome. We identified an abundance of small RNAs derived from non-coding RNA genes, including transfer RNAs, ribosomal RNAs as well as small nucleolar RNAs and small nuclear RNAs. Certain tRNA types were overrepresented as precursors for small RNAs. Further, we identified 79 novel small non-coding RNAs, not previously reported. We did not identify canonical small RNAs, like microRNAs and small interfering RNAs, and concluded that these do not exist in T. cruzi. This study has provided insights into the short transcriptome of a major human pathogen and provided starting points for further functional investigation of small RNAs and their biological roles.

Numerous fragmented spliceosomal introns, AT-AC splicing, and an unusual dynein gene expression pathway in Giardia lamblia

Spliceosomal introns are hallmarks of eukaryotic genomes, dividing coding regions into separate exons, which are joined during mRNA intron removal catalyzed by the spliceosome. With few known exceptions, spliceosomal introns are cis-spliced, that is, removed from one contiguous pre-mRNA transcript. The protistan intestinal parasite Giardia lamblia exhibits one of the most reduced eukaryotic genomes known, with short intergenic regions and only four known spliceosomal introns. Our genome-wide search for additional introns revealed four unusual cases of spliceosomal intron fragmentation, with consecutive exons of conserved protein-coding genes being dispersed to distant genomic sites. Independent transcripts are trans-spliced to yield contiguous mature mRNAs. Most strikingly, a dynein heavy chain subunit is both interrupted by two fragmented introns and also predicted to be assembled as two separately translated polypeptides, a remarkably complex expression pathway for a nuclear-encoded sequence. For each case, we observe extensive base-pairing potential between intron halves. This base pairing provides both a rationale for the in vivo association of independently transcribed mRNAs transcripts and the apparent specificity of splicing. Similar base-pairing potential in two cis-spliced G. lamblia introns suggests an evolutionary pathway whereby intron fragmentation of cis-spliced introns is permissible and a preliminary evolutionary step to complete gene fission. These results reveal remarkably complex genome dynamics in a severely genomically reduced parasite.

Homology-based annotation of non-coding RNAs in the genomes of Schistosoma mansoni and Schistosoma japonicum

BACKGROUND

Schistosomes are trematode parasites of the phylum Platyhelminthes. They are considered the most important of the human helminth parasites in terms of morbidity and mortality. Draft genome sequences are now available for Schistosoma mansoni and Schistosoma japonicum. Non-coding RNA (ncRNA) plays a crucial role in gene expression regulation, cellular function and defense, homeostasis, and pathogenesis. The genome-wide annotation of ncRNAs is a non-trivial task unless well-annotated genomes of closely related species are already available.

RESULTS

A homology search for structured ncRNA in the genome of S. mansoni resulted in 23 types of ncRNAs with conserved primary and secondary structure. Among these, we identified rRNA, snRNA, SL RNA, SRP, tRNAs and RNase P, and also possibly MRP and 7SK RNAs. In addition, we confirmed five miRNAs that have recently been reported in S. japonicum and found two additional homologs of known miRNAs. The tRNA complement of S. mansoni is comparable to that of the free-living planarian Schmidtea mediterranea, although for some amino acids differences of more than a factor of two are observed: Leu, Ser, and His are overrepresented, while Cys, Meth, and Ile are underrepresented in S. mansoni. On the other hand, the number of tRNAs in the genome of S. japonicum is reduced by more than a factor of four. Both schistosomes have a complete set of minor spliceosomal snRNAs. Several ncRNAs that are expected to exist in the S. mansoni genome were not found, among them the telomerase RNA, vault RNAs, and Y RNAs.

CONCLUSION

The ncRNA sequences and structures presented here represent the most complete dataset of ncRNA from any lophotrochozoan reported so far. This data set provides an important reference for further analysis of the genomes of schistosomes and indeed eukaryotic genomes at large.

Draft genome sequencing of giardia intestinalis assemblage B isolate GS: is human giardiasis caused by two different species?

Giardia intestinalis is a major cause of diarrheal disease worldwide and two major Giardia genotypes, assemblages A and B, infect humans. The genome of assemblage A parasite WB was recently sequenced, and the structurally compact 11.7 Mbp genome contains simplified basic cellular machineries and metabolism. We here performed 454 sequencing to 16× coverage of the assemblage B isolate GS, the only Giardia isolate successfully used to experimentally infect animals and humans. The two genomes show 77% nucleotide and 78% amino-acid identity in protein coding regions. Comparative analysis identified 28 unique GS and 3 unique WB protein coding genes, and the variable surface protein (VSP) repertoires of the two isolates are completely different. The promoters of several enzymes involved in the synthesis of the cyst-wall lack binding sites for encystation-specific transcription factors in GS. Several synteny-breaks were detected and verified. The tetraploid GS genome shows higher levels of overall allelic sequence polymorphism (0.5 versus <0.01% in WB). The genomic differences between WB and GS may explain some of the observed biological and clinical differences between the two isolates, and it suggests that assemblage A and B Giardia can be two different species.

Giardia intestinalis is a major contributor to the enormous burden of diarrheal diseases with 250 million symptomatic infections per year, and it is part of the WHO neglected disease initiative. Nonetheless, there is poor insight into how Giardia causes disease; it is not invasive, secretes no known toxin and both the duration and symptoms of giardiasis are highly variable. Currently, there are seven defined variants (assemblages) of G. intestinalis, with only assemblages A and B being known to infect humans. Although assemblage B is the most prevalent worldwide, it is inconclusive whether the various genotypes are associated with different disease outcomes. We have used the 454 sequencing technology to sequence the first assemblage B isolate, and the genome was compared to the earlier sequenced assemblage A isolate. Large genetic differences were detected in genes involved in survival of the parasite during infections. The genomic differences between assemblage A and B can explain some of the observed biological and clinical differences between the two assemblages. Our data suggest that assemblage A and B Giardia can be two different species. The identification of genomic differences between assemblages is indeed very important for further studies of the disease and in the development of new methods for diagnosis and treatment of giardiasis.

High throughput genome-wide survey of small RNAs from the parasitic protists Giardia intestinalis and Trichomonas vaginalis

RNA interference (RNAi) is a set of mechanisms which regulate gene expression in eukaryotes. Key elements of RNAi are small sense and antisense RNAs from 19 to 26 nt generated from double-stranded RNAs. MicroRNAs (miRNAs) are a major type of RNAi-associated small RNAs and are found in most eukaryotes studied to date. To investigate whether small RNAs associated with RNAi appear to be present in all eukaryotic lineages, and therefore present in the ancestral eukaryote, we studied two deep-branching protozoan parasites, Giardia intestinalis and Trichomonas vaginalis. Little is known about endogenous small RNAs involved in RNAi of these organisms. Using Illumina Solexa sequencing and genome-wide analysis of small RNAs from these distantly related deep-branching eukaryotes, we identified 10 strong miRNA candidates from Giardia and 11 from Trichomonas. We also found evidence of Giardia short-interfering RNAs potentially involved in the expression of variant-specific surface proteins. In addition, eight new small nucleolar RNAs from Trichomonas are identified. Our results indicate that miRNAs are likely to be general in ancestral eukaryotes and therefore are likely to be a universal feature of eukaryotes.

The RNA infrastructure: dark matter of the eukaryotic cell?

Eukaryotes express many functional non-protein-coding RNAs (ncRNAs) that participate in the processing and regulation of other RNA molecules. By focusing on connections between RNA-based processes, common patterns emerge that form a network-like RNA infrastructure. Owing to the intracellular movement of RNA during its processing (both between nuclear compartments and between the nucleus and cytoplasm), the RNA infrastructure contains both spatial and temporal connections. As research moves away from being protein-centric and focuses more on genomics, it is timely to explore these often 'hidden' aspects of the eukaryotic cell. The general and ancestral nature of most basic RNA-processing steps places a new focus on the generality of the spatial and temporal steps in RNA processing.