Genome-wide searching for pseudouridylation guide snoRNAs: analysis of the Saccharomyces cerevisiae genome

Peak analysis of cell-free RNA finds recurrently protected narrow regions with clinical potential

BACKGROUND

Cell-free RNAs (cfRNAs) can be detected in biofluids and have emerged as valuable disease biomarkers. Accurate identification of the fragmented cfRNA signals, especially those originating from pathological cells, is crucial for understanding their biological functions and clinical value. However, many challenges still need to be addressed for their application, including developing specific analysis methods and translating cfRNA fragments with biological support into clinical applications.

RESULTS

We present cfPeak, a novel method combining statistics and machine learning models to detect the fragmented cfRNA signals effectively. When test in real and artificial cfRNA sequencing (cfRNA-seq) data, cfPeak shows an improved performance compared with other applicable methods. We reveal that narrow cfRNA peaks preferentially overlap with protein binding sites, vesicle-sorting sites, structural sites, and novel small non-coding RNAs (sncRNAs). When applied in clinical cohorts, cfPeak identified cfRNA peaks in patients' plasma that enable cancer detection and are informative of cancer types and metastasis.

CONCLUSIONS

Our study fills the gap in the current small cfRNA-seq analysis at fragment-scale and builds a bridge to the scientific discovery in cfRNA fragmentomics. We demonstrate the significance of finding low abundant tissue-derived signals in small cfRNA and prove the feasibility for application in liquid biopsy.

Epitranscriptomic rRNA fingerprinting reveals tissue-of-origin and tumor-specific signatures

Mammalian ribosomal RNA (rRNA) molecules are highly abundant RNAs, decorated with over 220 rRNA modifications. Previous works have shown that some rRNA modification types can be dynamically regulated; however, how and when the mammalian rRNA modification landscape is remodeled remains largely unexplored. Here, we employ direct RNA sequencing to chart the human and mouse rRNA epitranscriptome across tissues, developmental stages, cell types, and disease. Our analyses reveal multiple rRNA sites that are differentially modified in a tissue- and/or developmental stage-specific manner, including previously unannotated modified sites. We demonstrate that rRNA modification patterns can be used for tissue and cell-type identification, which we hereby term "epitranscriptomic fingerprinting." We then explore rRNA modification patterns in normal-tumor matched samples from lung cancer patients, finding that epitranscriptomic fingerprinting accurately classifies clinical samples into normal and tumor groups from only 250 reads per sample, demonstrating the potential of rRNA modifications as diagnostic biomarkers.

Landscape of RNA pseudouridylation in archaeon Sulfolobus islandicus

Pseudouridine, one of the most abundant RNA modifications, is synthesized by stand-alone or RNA-guided pseudouridine synthases. Here, we comprehensively mapped pseudouridines in rRNAs, tRNAs and small RNAs in the archaeon Sulfolobus islandicus and identified Cbf5-associated H/ACA RNAs. Through genetic deletion and in vitro modification assays, we determined the responsible enzymes for these modifications. The pseudouridylation machinery in S. islandicus consists of the stand-alone enzymes aPus7 and aPus10, and six H/ACA RNA-guided enzymes that account for all identified pseudouridines. These H/ACA RNAs guide the modification of all eleven sites in rRNAs, two sites in tRNAs, and two sites in CRISPR RNAs. One H/ACA RNA shows exceptional versatility by targeting eight different sites. aPus7 and aPus10 are responsible for modifying positions 13, 54 and 55 in tRNAs. We identified four atypical H/ACA RNAs that lack the lower stem and the ACA motif and confirmed their function both in vivo and in vitro. Intriguingly, atypical H/ACA RNAs can be modified by Cbf5 in a guide-independent manner. Our data provide the first global view of pseudouridylation in archaea and reveal unexpected structures, substrates, and activities of archaeal H/ACA RNPs.

Current Views on Pathophysiology and Potential Therapeutic Targets in Sjögren's Syndrome: A Review from the Perspective of Viral Infections, Toll-like Receptors, and Long-Noncoding RNAs

Sjögren's syndrome (SS) is a rheumatic disease characterized by sicca and extraglandular symptoms, such as interstitial lung disease and renal tubular acidosis. SS potentially affects the prognosis of patients, especially in cases of complicated extraglandular symptoms; however, only symptomatic therapies against xerophthalmia and xerostomia are currently included in the practice guidelines as recommended therapies for SS. Considering that SS is presumed to be a multifactorial entity caused by genetic and environmental factors, a multidisciplinary approach is necessary to clarify the whole picture of its pathogenesis and to develop disease-specific therapies for SS. This review discusses past achievements and future prospects for pursuing the pathophysiology and therapeutic targets for SS, especially from the perspectives of viral infections, toll-like receptors (TLRs), long-noncoding RNAs (lncRNAs), and related signals. Based on the emerging roles of viral infections, TLRs, long-noncoding RNAs and related signals, antiviral therapy, hydroxychloroquine, and vitamin D may lower the risk of or mitigate SS. Janus-kinase (JAK) inhibitors are also potential novel therapeutic options for several rheumatic diseases involving the JAK-signal transducer and activator of transcription pathways, which are yet to be ascertained in a randomized controlled study targeting SS.

Web Services for RNA-RNA Interaction Prediction

Non-coding RNAs have various biological functions such as translational regulation, and RNA-RNA interactions play essential roles in the mechanisms of action of these RNAs. Therefore, RNA-RNA interaction prediction is an important problem in bioinformatics, and many tools have been developed for the computational prediction of RNA-RNA interactions. In addition to the development of novel algorithms with high accuracy, the development and maintenance of web services is essential for enhancing usability by experimental biologists. In this review, we survey web services for RNA-RNA interaction predictions and introduce how to use primary web services. We present various prediction tools, including general interaction prediction tools, prediction tools for specific RNA classes, and RNA-RNA interaction-based RNA design tools. Additionally, we discuss the future perspectives of the development of RNA-RNA interaction prediction tools and the sustainability of web services.

Cms1 coordinates stepwise local 90S pre-ribosome assembly with timely snR83 release

Ribosome synthesis begins in the nucleolus with 90S pre-ribosome construction, but little is known about how the many different snoRNAs that modify the pre-rRNA are timely guided to their target sites. Here, we report a role for Cms1 in such a process. Initially, we discovered CMS1 as a null suppressor of a nop14 mutant impaired in Rrp12-Enp1 factor recruitment to the 90S. Further investigations detected Cms1 at the 18S rRNA 3' major domain of an early 90S that carried H/ACA snR83, which is known to guide pseudouridylation at two target sites within the same subdomain. Cms1 co-precipitates with many 90S factors, but Rrp12-Enp1 encircling the 3' major domain in the mature 90S is decreased. We suggest that Cms1 associates with the 3' major domain during early 90S biogenesis to restrict premature Rrp12-Enp1 binding but allows snR83 to timely perform its modification role before the next 90S assembly steps coupled with Cms1 release take place.

Human dyskerin binds to cytoplasmic H/ACA-box-containing transcripts affecting nuclear hormone receptor dependence

Background

Dyskerin is a nuclear protein involved in H/ACA box snoRNA-guided uridine modification of RNA. In humans, its defective function is associated with cancer development and induces specific post-transcriptional alterations of gene expression. In this study, we seek to unbiasedly identify mRNAs regulated by dyskerin in human breast cancer-derived cells.

Results

We find that dyskerin depletion affects the expression and the association with polysomes of selected mRNA isoforms characterized by the retention of H/ACA box snoRNA-containing introns. These snoRNA retaining transcripts (snoRTs) are bound by dyskerin in the cytoplasm in the form of shorter 3′ snoRT fragments. We then characterize the whole cytoplasmic dyskerin RNA interactome and find both H/ACA box snoRTs and protein-coding transcripts which may be targeted by the snoRTs’ guide properties. Since a fraction of these protein-coding transcripts is involved in the nuclear hormone receptor binding, we test to see if this specific activity is affected by dyskerin. Obtained results indicate that dyskerin dysregulation may alter the dependence on nuclear hormone receptor ligands in breast cancer cells. These results are paralleled by consistent observations on the outcome of primary breast cancer patients stratified according to their tumor hormonal status. Accordingly, experiments in nude mice show that the reduction of dyskerin levels in estrogen-dependent cells favors xenograft development in the absence of estrogen supplementation.

Conclusions

Our work suggests a cytoplasmic function for dyskerin which could affect mRNA post-transcriptional networks relevant for nuclear hormone receptor functions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02746-3.

Concerted modification of nucleotides at functional centers of the ribosome revealed by single-molecule RNA modification profiling

Nucleotides in RNA and DNA are chemically modified by numerous enzymes that alter their function. Eukaryotic ribosomal RNA (rRNA) is modified at more than 100 locations, particularly at highly conserved and functionally important nucleotides. During ribosome biogenesis, modifications are added at various stages of assembly. The existence of differently modified classes of ribosomes in normal cells is unknown because no method exists to simultaneously evaluate the modification status at all sites within a single rRNA molecule. Using a combination of yeast genetics and nanopore direct RNA sequencing, we developed a reliable method to track the modification status of single rRNA molecules at 37 sites in 18 S rRNA and 73 sites in 25 S rRNA. We use our method to characterize patterns of modification heterogeneity and identify concerted modification of nucleotides found near functional centers of the ribosome. Distinct, undermodified subpopulations of rRNAs accumulate upon loss of Dbp3 or Prp43 RNA helicases, suggesting overlapping roles in ribosome biogenesis. Modification profiles are surprisingly resistant to change in response to many genetic and acute environmental conditions that affect translation, ribosome biogenesis, and pre-mRNA splicing. The ability to capture single-molecule RNA modification profiles provides new insights into the roles of nucleotide modifications in RNA function.

Emergence of the primordial pre-60S from the 90S pre-ribosome

Synthesis of ribosomes begins in the nucleolus with formation of the 90S pre-ribosome, during which the pre-40S and pre-60S pathways diverge by pre-rRNA cleavage. However, it remains unclear how, after this uncoupling, the earliest pre-60S subunit continues to develop. Here, we reveal a large-subunit intermediate at the beginning of its construction when still linked to the 90S, the precursor to the 40S subunit. This primordial pre-60S is characterized by the SPOUT domain methyltransferase Upa1-Upa2, large α-solenoid scaffolds, Mak5, one of several RNA helicases, and two small nucleolar RNA (snoRNAs), C/D box snR190 and H/ACA box snR37. The emerging pre-60S does not efficiently disconnect from the 90S pre-ribosome in a dominant mak5 helicase mutant, allowing a 70-nm 90S-pre-60S bipartite particle to be visualized by electron microscopy. Our study provides insight into the assembly pathway when the still-connected nascent 40S and 60S subunits are beginning to separate.

Guide RNA acrobatics: positioning consecutive uridines for pseudouridylation by H/ACA pseudouridylation loops with dual guide capacity

Site-specific pseudouridylation of human ribosomal and spliceosomal RNAs is directed by H/ACA guide RNAs composed of two hairpins carrying internal pseudouridylation guide loops. The distal "antisense" sequences of the pseudouridylation loop base-pair with the target RNA to position two unpaired target nucleotides 5'-UN-3', including the 5' substrate U, under the base of the distal stem topping the guide loop. Therefore, each pseudouridylation loop is expected to direct synthesis of a single pseudouridine (Ψ) in the target sequence. However, in this study, genetic depletion and restoration and RNA mutational analyses demonstrate that at least four human H/ACA RNAs (SNORA53, SNORA57, SCARNA8, and SCARNA1) carry pseudouridylation loops supporting efficient and specific synthesis of two consecutive pseudouridines (ΨΨ or ΨNΨ) in the 28S (Ψ3747/Ψ3749), 18S (Ψ1045/Ψ1046), and U2 (Ψ43/Ψ44 and Ψ89/Ψ91) RNAs, respectively. In order to position two substrate Us for pseudouridylation, the dual guide loops form alternative base-pairing interactions with their target RNAs. This remarkable structural flexibility of dual pseudouridylation loops provides an unexpected versatility for RNA-directed pseudouridylation without compromising its efficiency and accuracy. Besides supporting synthesis of at least 6% of human ribosomal and spliceosomal Ψs, evidence indicates that dual pseudouridylation loops also participate in pseudouridylation of yeast and archaeal rRNAs.

Nhp2 is a reader of H2AQ105me and part of a network integrating metabolism with rRNA synthesis

Ribosome biogenesis is an essential cellular process that requires integration of extracellular cues, such as metabolic state, with intracellular signalling, transcriptional regulation and chromatin accessibility at the ribosomal DNA. Here, we demonstrate that the recently identified histone modification, methylation of H2AQ105 (H2AQ105me), is an integral part of a dynamic chromatin network at the rDNA locus. Its deposition depends on a functional mTor signalling pathway and acetylation of histone H3 at position K56, thus integrating metabolic and proliferative signals. Furthermore, we identify a first epigenetic reader of this modification, the ribonucleoprotein Nhp2, which specifically recognizes H2AQ105me. Based on functional and proteomic data, we suggest that Nhp2 functions as an adapter to bridge rDNA chromatin with components of the small subunit processome to efficiently coordinate transcription of rRNA with its post‐transcriptional processing. We support this by showing that an H2AQ105A mutant has a mild defect in early processing of rRNA.

CRISPR mediated genome editing, a tool to dissect RNA modification processes

Though over 100 distinct RNA modifications have been identified, the roles for many of these modifications in vivo remain unknown. Genome editing is one tool investigators are using to better understand the roles these modifications play and the consequences of their absence. In this chapter, we describe how CRISPR mediated genome editing can be used to interrogate the process of RNA modification in C. albicans. Furthermore, we discuss how the protocols described can be altered to meet experimental demands. The underlying theory on which these protocols are based are applicable to a variety of model systems. The protocols described utilize the widely used S. pyogenes Cas9, but the field of genome editing is quickly evolving. We discuss the recent developments of more flexible CRISPR systems that can target a greater number of sites in the genome. These and other advancements make CRISPR mediated genome editing a practical methodology to investigate RNA modification.

Discovery of 17 conserved structural RNAs in fungi

Many non-coding RNAs with known functions are structurally conserved: their intramolecular secondary and tertiary interactions are maintained across evolutionary time. Consequently, the presence of conserved structure in multiple sequence alignments can be used to identify candidate functional non-coding RNAs. Here, we present a bioinformatics method that couples iterative homology search with covariation analysis to assess whether a genomic region has evidence of conserved RNA structure. We used this method to examine all unannotated regions of five well-studied fungal genomes (Saccharomyces cerevisiae, Candida albicans, Neurospora crassa, Aspergillus fumigatus, and Schizosaccharomyces pombe). We identified 17 novel structurally conserved non-coding RNA candidates, which include four H/ACA box small nucleolar RNAs, four intergenic RNAs and nine RNA structures located within the introns and untranslated regions (UTRs) of mRNAs. For the two structures in the 3' UTRs of the metabolic genes GLY1 and MET13, we performed experiments that provide evidence against them being eukaryotic riboswitches.

H/ACA Small Ribonucleoproteins: Structural and Functional Comparison Between Archaea and Eukaryotes

During ribosome synthesis, ribosomal RNA is modified through the formation of many pseudouridines and methylations which contribute to ribosome function across all domains of life. In archaea and eukaryotes, pseudouridylation of rRNA is catalyzed by H/ACA small ribonucleoproteins (sRNPs) utilizing different H/ACA guide RNAs to identify target uridines for modification. H/ACA sRNPs are conserved in archaea and eukaryotes, as they share a common general architecture and function, but there are also several notable differences between archaeal and eukaryotic H/ACA sRNPs. Due to the higher protein stability in archaea, we have more information on the structure of archaeal H/ACA sRNPs compared to eukaryotic counterparts. However, based on the long history of yeast genetic and other cellular studies, the biological role of H/ACA sRNPs during ribosome biogenesis is better understood in eukaryotes than archaea. Therefore, this review provides an overview of the current knowledge on H/ACA sRNPs from archaea, in particular their structure and function, and relates it to our understanding of the roles of eukaryotic H/ACA sRNP during eukaryotic ribosome synthesis and beyond. Based on this comparison of our current insights into archaeal and eukaryotic H/ACA sRNPs, we discuss what role archaeal H/ACA sRNPs may play in the formation of ribosomes.

Identification of protein-protected mRNA fragments and structured excised intron RNAs in human plasma by TGIRT-seq peak calling

Human plasma contains > 40,000 different coding and non-coding RNAs that are potential biomarkers for human diseases. Here, we used thermostable group II intron reverse transcriptase sequencing (TGIRT-seq) combined with peak calling to simultaneously profile all RNA biotypes in apheresis-prepared human plasma pooled from healthy individuals. Extending previous TGIRT-seq analysis, we found that human plasma contains largely fragmented mRNAs from > 19,000 protein-coding genes, abundant full-length, mature tRNAs and other structured small non-coding RNAs, and less abundant tRNA fragments and mature and pre-miRNAs. Many of the mRNA fragments identified by peak calling correspond to annotated protein-binding sites and/or have stable predicted secondary structures that could afford protection from plasma nucleases. Peak calling also identified novel repeat RNAs, miRNA-sized RNAs, and putatively structured intron RNAs of potential biological, evolutionary, and biomarker significance, including a family of full-length excised intron RNAs, subsets of which correspond to mirtron pre-miRNAs or agotrons.

A Census and Categorization Method of Epitranscriptomic Marks

In the past few years, thorough investigation of chemical modifications operated in the cells on ribonucleic acid (RNA) molecules is gaining momentum. This new field of research has been dubbed “epitranscriptomics”, in analogy to best-known epigenomics, to stress the potential of ensembles of RNA modifications to constitute a post-transcriptional regulatory layer of gene expression orchestrated by writer, reader, and eraser RNA-binding proteins (RBPs). In fact, epitranscriptomics aims at identifying and characterizing all functionally relevant changes involving both non-substitutional chemical modifications and editing events made to the transcriptome. Indeed, several types of RNA modifications that impact gene expression have been reported so far in different species of cellular RNAs, including ribosomal RNAs, transfer RNAs, small nuclear RNAs, messenger RNAs, and long non-coding RNAs. Supporting functional relevance of this largely unknown regulatory mechanism, several human diseases have been associated directly to RNA modifications or to RBPs that may play as effectors of epitranscriptomic marks. However, an exhaustive epitranscriptome’s characterization, aimed to systematically classify all RNA modifications and clarify rules, actors, and outcomes of this promising regulatory code, is currently not available, mainly hampered by lack of suitable detecting technologies. This is an unfortunate limitation that, thanks to an unprecedented pace of technological advancements especially in the sequencing technology field, is likely to be overcome soon. Here, we review the current knowledge on epitranscriptomic marks and propose a categorization method based on the reference ribonucleotide and its rounds of modifications (“stages”) until reaching the given modified form. We believe that this classification scheme can be useful to coherently organize the expanding number of discovered RNA modifications.

Reducing the structure bias of RNA-Seq reveals a large number of non-annotated non-coding RNA

The study of RNA expression is the fastest growing area of genomic research. However, despite the dramatic increase in the number of sequenced transcriptomes, we still do not have accurate estimates of the number and expression levels of non-coding RNA genes. Non-coding transcripts are often overlooked due to incomplete genome annotation. In this study, we use annotation-independent detection of RNA reads generated using a reverse transcriptase with low structure bias to identify non-coding RNA. Transcripts between 20 and 500 nucleotides were filtered and crosschecked with non-coding RNA annotations revealing 111 non-annotated non-coding RNAs expressed in different cell lines and tissues. Inspecting the sequence and structural features of these transcripts indicated that 60% of these transcripts correspond to new snoRNA and tRNA-like genes. The identified genes exhibited features of their respective families in terms of structure, expression, conservation and response to depletion of interacting proteins. Together, our data reveal a new group of RNA that are difficult to detect using standard gene prediction and RNA sequencing techniques, suggesting that reliance on actual gene annotation and sequencing techniques distorts the perceived architecture of the human transcriptome.

5'-End sequencing in Saccharomyces cerevisiae offers new insights into 5'-ends of tRNAH is and snoRNAs

tRNA^{H is} guanylyltransferase (Thg1) specifies eukaryotic tRNA^{H is} identity by catalysing a 3'-5' non-Watson-Crick (WC) addition of guanosine to the 5'-end of tRNA^{H is} . Thg1 family enzymes in Archaea and Bacteria, called Thg1-like proteins (TLPs), catalyse a similar but distinct 3'-5' addition in an exclusively WC-dependent manner. Here, a genetic system in Saccharomyces cerevisiae was employed to further assess the biochemical differences between Thg1 and TLPs. Utilizing a novel 5'-end sequencing pipeline, we find that a Bacillus thuringiensis TLP sustains the growth of a thg1Δ strain by maintaining a WC-dependent addition of U_-1 across from A₇₃ . Additionally, we observe 5'-end heterogeneity in S. cerevisiae small nucleolar RNAs (snoRNAs), an observation that may inform methods of annotation and mechanisms of snoRNA processing.

Nucleomorph Small RNAs in Cryptophyte and Chlorarachniophyte Algae

The regulation of gene expression and RNA maturation underlies fundamental processes such as cell homeostasis, development, and stress acclimation. The biogenesis and modification of RNA is tightly controlled by an array of regulatory RNAs and nucleic acid-binding proteins. While the role of small RNAs (sRNAs) in gene expression has been studied in-depth in select model organisms, little is known about sRNA biology across the eukaryotic tree of life. We used deep sequencing to explore the repertoires of sRNAs encoded by the miniaturized, endosymbiotically derived “nucleomorph” genomes of two single-celled algae, the cryptophyte Guillardia theta and the chlorarachniophyte Bigelowiella natans. A total of 32.3 and 35.3 million reads were generated from G. theta and B. natans, respectively. In G. theta, we identified nucleomorph U1, U2, and U4 spliceosomal small nuclear RNAs (snRNAs) as well as 11 C/D box small nucleolar RNAs (snoRNAs), five of which have potential plant and animal homologs. The snoRNAs are predicted to perform 2′-O methylation of rRNA (but not snRNA). In B. natans, we found the previously undetected 5S rRNA as well as six orphan sRNAs. Analysis of chlorarachniophyte snRNAs shed light on the removal of the miniature 18–21 nt introns found in B. natans nucleomorph genes. Neither of the nucleomorph genomes appears to encode RNA pseudouridylation machinery, and U5 snRNA cannot be found in the cryptophyte G. theta. Considering the central roles of U5 snRNA and RNA modifications in other organisms, cytoplasm-to-nucleomorph RNA shuttling in cryptophyte algae is a distinct possibility.

Post-transcriptional pseudouridylation in mRNA as well as in some major types of noncoding RNAs

Pseudouridylation is a post-transcriptional isomerization reaction that converts a uridine to a pseudouridine (Ψ) within an RNA chain. Ψ has chemical properties that are distinct from that of uridine and any other known nucleotides. Experimental data accumulated thus far have indicated that Ψ is present in many different types of RNAs, including coding and noncoding RNAs. Ψ is particularly concentrated in rRNA and spliceosomal snRNAs, and plays an important role in protein translation and pre-mRNA splicing, respectively. Ψ has also been found in mRNA, but its function there remains essentially unknown. In this review, we discuss the mechanisms and functions of RNA pseudouridylation, focusing on rRNA, snRNA and mRNA. We also discuss the methods, which have been developed to detect Ψs in RNAs. This article is part of a Special Issue entitled: mRNA modifications in gene expression control edited by Dr. Soller Matthias and Dr. Fray Rupert.

RNA-Seq employing a novel rRNA depletion strategy reveals a rich repertoire of snoRNAs in Euglena gracilis including box C/D and Ψ-guide RNAs targeting the modification of rRNA extremities

Previous mRNA transcriptome studies of Euglena gracilis have shown that this organism possesses a large and diverse complement of protein coding genes; however, the study of non-coding RNA classes has been limited. The natural extensive fragmentation of the E. gracilis large subunit ribosomal RNA presents additional barriers to the identification of non-coding RNAs as size-selected small RNA libraries will be dominated by rRNA sequences. In this study we have developed a strategy to significantly reduce rRNA amplification prior to RNA-Seq analysis thereby producing a ncRNA library allowing for the identification of many new E. gracilis small RNAs. Library analysis reveals 113 unique new small nucleolar (sno) RNAs and a large collection of snoRNA isoforms, as well as the first significant collection of nuclear tRNAs in this organism. A 3' end AGAUGN consensus motif and conserved structural features can now be defined for E. gracilis pseudouridine guide RNAs. snoRNAs of both classes were identified that target modification of the 3' extremities of rRNAs utilizing predicted base-pairing interactions with internally transcribed spacers (ITS), providing insight into the timing of steps in rRNA maturation. Cumulatively, this represents the most comprehensive analysis of small ncRNAs in Euglena gracilis to date.

The fungal snoRNAome

Small nucleolar RNAs (snoRNAs) are essential players in the rRNA biogenesis due to their involvement in the nucleolytic processing of the precursor and the subsequent guidance of nucleoside modifications. Within the kingdom Fungi, merely a few species-specific surveys have explored their snoRNA repertoire. However, the wide range of the snoRNA landscape spanning all major fungal lineages has not been mapped so far, mainly because of missing tools for automatized snoRNA detection and functional analysis. For the first time, we report here a comprehensive inventory of fungal snoRNAs together with a functional analysis and an in-depth investigation of their evolutionary history including innovations, deletions, and target switches. This large-scale analysis, incorporating more than 120 snoRNA families with more than 7700 individual snoRNA sequences, catalogs and clarifies the landscape of fungal snoRNA families, assigns functions to previously orphan snoRNAs, and increases the number of sequences by 450%. We also show that the snoRNAome is subject to ongoing rearrangements and adaptations, e.g., through lineage-specific targets and redundant guiding functions.

The Glaciozyma antarctica genome reveals an array of systems that provide sustained responses towards temperature variations in a persistently cold habitat

Extremely low temperatures present various challenges to life that include ice formation and effects on metabolic capacity. Psyhcrophilic microorganisms typically have an array of mechanisms to enable survival in cold temperatures. In this study, we sequenced and analysed the genome of a psychrophilic yeast isolated in the Antarctic region, Glaciozyma antarctica. The genome annotation identified 7857 protein coding sequences. From the genome sequence analysis we were able to identify genes that encoded for proteins known to be associated with cold survival, in addition to annotating genes that are unique to G. antarctica. For genes that are known to be involved in cold adaptation such as anti-freeze proteins (AFPs), our gene expression analysis revealed that they were differentially transcribed over time and in response to different temperatures. This indicated the presence of an array of adaptation systems that can respond to a changing but persistent cold environment. We were also able to validate the activity of all the AFPs annotated where the recombinant AFPs demonstrated anti-freeze capacity. This work is an important foundation for further collective exploration into psychrophilic microbiology where among other potential, the genes unique to this species may represent a pool of novel mechanisms for cold survival.

Identification of 15 candidate structured noncoding RNA motifs in fungi by comparative genomics

BACKGROUND

With the development of rapid and inexpensive DNA sequencing, the genome sequences of more than 100 fungal species have been made available. This dataset provides an excellent resource for comparative genomics analyses, which can be used to discover genetic elements, including noncoding RNAs (ncRNAs). Bioinformatics tools similar to those used to uncover novel ncRNAs in bacteria, likewise, should be useful for searching fungal genomic sequences, and the relative ease of genetic experiments with some model fungal species could facilitate experimental validation studies.

RESULTS

We have adapted a bioinformatics pipeline for discovering bacterial ncRNAs to systematically analyze many fungal genomes. This comparative genomics pipeline integrates information on conserved RNA sequence and structural features with alternative splicing information to reveal fungal RNA motifs that are candidate regulatory domains, or that might have other possible functions. A total of 15 prominent classes of structured ncRNA candidates were identified, including variant HDV self-cleaving ribozyme representatives, atypical snoRNA candidates, and possible structured antisense RNA motifs. Candidate regulatory motifs were also found associated with genes for ribosomal proteins, S-adenosylmethionine decarboxylase (SDC), amidase, and HexA protein involved in Woronin body formation. We experimentally confirm that the variant HDV ribozymes undergo rapid self-cleavage, and we demonstrate that the SDC RNA motif reduces the expression of SAM decarboxylase by translational repression. Furthermore, we provide evidence that several other motifs discovered in this study are likely to be functional ncRNA elements.

CONCLUSIONS

Systematic screening of fungal genomes using a computational discovery pipeline has revealed the existence of a variety of novel structured ncRNAs. Genome contexts and similarities to known ncRNA motifs provide strong evidence for the biological and biochemical functions of some newly found ncRNA motifs. Although initial examinations of several motifs provide evidence for their likely functions, other motifs will require more in-depth analysis to reveal their functions.

The yeast noncoding RNA interaction network

This article describes the creation of the first expert manually curated noncoding RNA interaction networks for S. cerevisiae The RNA-RNA and RNA-protein interaction networks have been carefully extracted from the experimental literature and made available through the IntAct database (www.ebi.ac.uk/intact). We provide an initial network analysis and compare their properties to the much larger protein-protein interaction network. We find that the proteins that bind to ncRNAs in the network contain only a small proportion of classical RNA binding domains. We also see an enrichment of WD40 domains suggesting their direct involvement in ncRNA interactions. We discuss the challenges in collecting noncoding RNA interaction data and the opportunities for worldwide collaboration to fill the unmet need for this data.

Synthesis, Function, and Heterogeneity of snoRNA-Guided Posttranscriptional Nucleoside Modifications in Eukaryotic Ribosomal RNAs

Ribosomal RNAs contain numerous 2'-O-methylated nucleosides and pseudouridines. Methylation of the 2' oxygen of ribose moieties and isomerization of uridines into pseudouridines are catalyzed by C/D and H/ACA small nucleolar ribonucleoprotein particles, respectively. We review the composition, structure, and mode of action of archaeal and eukaryotic C/D and H/ACA particles. Most rRNA modifications cluster in functionally crucial regions of the rRNAs, suggesting they play important roles in translation. Some of these modifications promote global translation efficiency or modulate translation fidelity. Strikingly, recent quantitative nucleoside modification profiling methods have revealed that a subset of modification sites is not always fully modified. The finding of such ribosome heterogeneity is in line with the concept of specialized ribosomes that could preferentially translate specific mRNAs. This emerging concept is supported by findings that some human diseases are caused by defects in the rRNA modification machinery correlated with a significant alteration of IRES-dependent translation.

Transcriptional regulation of snRNAs and its significance for plant development

Small nuclear RNA (snRNA) represents a distinct class of non-coding RNA molecules. As these molecules have fundamental roles in RNA metabolism, including pre-mRNA splicing and ribosomal RNA processing, it is essential that their transcription be tightly regulated in eukaryotic cells. The genome of each organism contains hundreds of snRNA genes. Although the structures of these genes are highly diverse among organisms, the trans-acting factors that regulate snRNA transcription are evolutionarily conserved. Recent studies of the Arabidopsis thaliana srd2-1 mutant, which is defective in the snRNA transcription factor, provide insight into the physiological significance of snRNA regulation in plant development. Here, I review the current understanding of the molecular mechanisms underlying snRNA transcription.

Tuning the ribosome: The influence of rRNA modification on eukaryotic ribosome biogenesis and function

rRNAs are extensively modified during their transcription and subsequent maturation in the nucleolus, nucleus and cytoplasm. RNA modifications, which are installed either by snoRNA-guided or by stand-alone enzymes, generally stabilize the structure of the ribosome. However, they also cluster at functionally important sites of the ribosome, such as the peptidyltransferase center and the decoding site, where they facilitate efficient and accurate protein synthesis. The recent identification of sites of substoichiometric 2'-O-methylation and pseudouridylation has overturned the notion that all rRNA modifications are constitutively present on ribosomes, highlighting nucleotide modifications as an important source of ribosomal heterogeneity. While the mechanisms regulating partial modification and the functions of specialized ribosomes are largely unknown, changes in the rRNA modification pattern have been observed in response to environmental changes, during development, and in disease. This suggests that rRNA modifications may contribute to the translational control of gene expression.

Identification of Non-Coding RNAs in the Candida parapsilosis Species Group

The Candida CTG clade is a monophyletic group of fungal species that translates CTG as serine, and includes the pathogens Candida albicans and Candida parapsilosis. Research has typically focused on identifying protein-coding genes in these species. Here, we use bioinformatic and experimental approaches to annotate known classes of non-coding RNAs in three CTG-clade species, Candida parapsilosis, Candida orthopsilosis and Lodderomyces elongisporus. We also update the annotation of ncRNAs in the C. albicans genome. The majority of ncRNAs identified were snoRNAs. Approximately 50% of snoRNAs (including most of the C/D box class) are encoded in introns. Most are within mono- and polycistronic transcripts with no protein coding potential. Five polycistronic clusters of snoRNAs are highly conserved in fungi. In polycistronic regions, splicing occurs via the classical pathway, as well as by nested and recursive splicing. We identified spliceosomal small nuclear RNAs, the telomerase RNA component, signal recognition particle, RNase P RNA component and the related RNase MRP RNA component in all three genomes. Stem loop IV of the U2 spliceosomal RNA and the associated binding proteins were lost from the ancestor of C. parapsilosis and C. orthopsilosis, following the divergence from L. elongisporus. The RNA component of the MRP is longer in C. parapsilosis, C. orthopsilosis and L. elongisporus than in S. cerevisiae, but is substantially shorter than in C. albicans.

Identification and Characterization of Small Noncoding RNAs in Genome Sequences of the Edible Fungus Pleurotus ostreatus

Noncoding RNAs (ncRNAs) have been identified in many fungi. However, no genome-scale identification of ncRNAs has been inventoried for basidiomycetes. In this research, we detected 254 small noncoding RNAs (sncRNAs) in a genome assembly of an isolate (CCEF00389) of Pleurotus ostreatus, which is a widely cultivated edible basidiomycetous fungus worldwide. The identified sncRNAs include snRNAs, snoRNAs, tRNAs, and miRNAs. SnRNA U1 was not found in CCEF00389 genome assembly and some other basidiomycetous genomes by BLASTn. This implies that if snRNA U1 of basidiomycetes exists, it has a sequence that varies significantly from other organisms. By analyzing the distribution of sncRNA loci, we found that snRNAs and most tRNAs (88.6%) were located in pseudo-UTR regions, while miRNAs are commonly found in introns. To analyze the evolutionary conservation of the sncRNAs in P. ostreatus, we aligned all 254 sncRNAs to the genome assemblies of some other Agaricomycotina fungi. The results suggest that most sncRNAs (77.56%) were highly conserved in P. ostreatus, and 20% were conserved in Agaricomycotina fungi. These findings indicate that most sncRNAs of P. ostreatus were not conserved across Agaricomycotina fungi.

Accurate placement of substrate RNA by Gar1 in H/ACA RNA-guided pseudouridylation

H/ACA RNA-guided ribonucleoprotein particle (RNP), the most complicated RNA pseudouridylase so far known, uses H/ACA guide RNA for substrate capture and four proteins (Cbf5, Nop10, L7Ae and Gar1) for pseudouridylation. Although it was shown that Gar1 not only facilitates the product release, but also enhances the catalytic activity, the chemical role that Gar1 plays in this complicated machinery is largely unknown. Kinetics measurement on Pyrococcus furiosus RNPs at different temperatures making use of fluorescence anisotropy showed that Gar1 reduces the catalytic barrier through affecting the activation entropy instead of enthalpy. Site-directed mutagenesis combined with molecular dynamics simulations demonstrated that V149 in the thumb loop of Cbf5 is critical in placing the target uridine to the right position toward catalytic D85 of Cbf5. The enzyme elegantly aligns the position of uridine in the catalytic site with the help of Gar1. In addition, conversion of uridine to pseudouridine results in a rigid syn configuration of the target nucleotide in the active site and causes Gar1 to pull out the thumb. Both factors guarantee the efficient release of the product.

Enhanced leavening ability of baker’s yeast by overexpression of SNR84 with PGM2 deletion

Dough-leavening ability is one of the main aspects considered when selecting a baker’s yeast strain for baking industry. Generally, modification of maltose metabolic pathway and known regulatory networks of maltose metabolism were used to increase maltose metabolism to improve leavening ability in lean dough. In this study, we focus on the effects of PGM2 (encoding for the phosphoglucomutase) and SNR84 (encoding for the H/ACA snoRNA) that are not directly related to both the maltose metabolic pathway and known regulatory networks of maltose metabolism on the leavening ability of baker’s yeast in lean dough. The results show that the modifications on PGM2 and/or SNR84 are effective ways in improving leavening ability of baker’s yeast in lean dough. Deletion of PGM2 decreased cellular glucose-1-phosphate and overexpression of SNR84 increased the maltose permease activity. These changes resulted in 11, 19 and 21 % increases of the leavening ability for PGM2 deletion, SNR84 overexpression and SNR84 overexpression combining deleted PGM2, respectively.

Genome-wide analysis of small nucleolar RNAs of Leishmania major reveals a rich repertoire of RNAs involved in modification and processing of rRNA

Trypanosomatids are protozoan parasites and the causative agent of infamous infectious diseases. These organisms regulate their gene expression mainly at the post-transcriptional level and possess characteristic RNA processing mechanisms. In this study, we analyzed the complete repertoire of Leishmania major small nucleolar (snoRNA) RNAs by performing RNA-seq analysis on RNAs that were affinity-purified using the C/D snoRNA core protein, SNU13, and the H/ACA core protein, NHP2. This study revealed a large collection of C/D and H/ACA snoRNAs, organized in gene clusters generally containing both snoRNA types. Abundant snoRNAs were identified and predicted to guide trypanosome-specific rRNA cleavages. The repertoire of snoRNAs was compared to that of the closely related Trypanosoma brucei, and 80% of both C/D and H/ACA molecules were found to have functional homologues. The comparative analyses elucidated how snoRNAs evolved to generate molecules with analogous functions in both species. Interestingly, H/ACA RNAs have great flexibility in their ability to guide modifications, and several of the RNA species can guide more than one modification, compensating for the presence of single hairpin H/ACA snoRNA in these organisms. Placing the predicted modifications on the rRNA secondary structure revealed hypermodification regions mostly in domains which are modified in other eukaryotes, in addition to trypanosome-specific modifications.

Improvement of stress tolerance and leavening ability under multiple baking-associated stress conditions by overexpression of the SNR84 gene in baker's yeast

During the bread-making process, industrial baker's yeast cells are exposed to multiple baking-associated stresses, such as elevated high-temperature, high-sucrose and freeze-thaw stresses. There is a high demand for baker's yeast strains that could withstand these stresses with high leavening ability. The SNR84 gene encodes H/ACA snoRNA (small nucleolar RNA), which is known to be involved in pseudouridylation of the large subunit rRNA. However, the function of the SNR84 gene in baker's yeast coping with baking-associated stresses remains unclear. In this study, we explored the effect of SNR84 overexpression on baker's yeast which was exposed to high-temperature, high-sucrose and freeze-thaw stresses. These results suggest that overexpression of the SNR84 gene conferred tolerance of baker's yeast cells to high-temperature, high-sucrose and freeze-thaw stresses and enhanced their leavening ability in high-sucrose and freeze-thaw dough. These findings could provide a valuable insight for breeding of novel stress-resistant baker's yeast strains that are useful for baking.

Non-coding RNA gene families in the genomes of anopheline mosquitoes

Background

Only a small fraction of the mosquito species of the genus Anopheles are able to transmit malaria, one of the biggest killer diseases of poverty, which is mostly prevalent in the tropics. This diversity has genetic, yet unknown, causes. In a further attempt to contribute to the elucidation of these variances, the international “Anopheles Genomes Cluster Consortium” project (a.k.a. “16 Anopheles genomes project”) was established, aiming at a comprehensive genomic analysis of several anopheline species, most of which are malaria vectors. In the frame of the international consortium carrying out this project our team studied the genes encoding families of non-coding RNAs (ncRNAs), concentrating on four classes: microRNA (miRNA), ribosomal RNA (rRNA), small nuclear RNA (snRNA), and in particular small nucleolar RNA (snoRNA) and, finally, transfer RNA (tRNA).

Results

Our analysis was carried out using, exclusively, computational approaches, and evaluating both the primary NGS reads as well as the respective genome assemblies produced by the consortium and stored in VectorBase; moreover, the results of RNAseq surveys in cases in which these were available and meaningful were also accessed in order to obtain supplementary data, as were “pre-genomic era” sequence data stored in nucleic acid databases. The investigation included the identification and analysis, in most species studied, of ncRNA genes belonging to several families, as well as the analysis of the evolutionary relations of some of those genes in cross-comparisons to other members of the genus Anopheles.

Conclusions

Our study led to the identification of members of these gene families in the majority of twenty different anopheline taxa. A set of tools for the study of the evolution and molecular biology of important disease vectors has, thus, been obtained.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1038) contains supplementary material, which is available to authorized users.

Insight into the mechanisms and functions of spliceosomal snRNA pseudouridylation

Pseudouridines (Ψs) are the most abundant and highly conserved modified nucleotides found in various stable RNAs of all organisms. Most Ψs are clustered in regions that are functionally important for pre-mRNA splicing. Ψ has an extra hydrogen bond donor that endows RNA molecules with distinct properties that contribute significantly to RNA-mediated cellular processes. Experimental data indicate that spliceosomal snRNA pseudouridylation can be catalyzed by both RNA-dependent and RNA-independent mechanisms. Recent work has also demonstrated that pseudouridylation can be induced at novel positions under stress conditions, suggesting a regulatory role for Ψ.

Deep profiling of the novel intermediate-size noncoding RNAs in intraerythrocytic Plasmodium falciparum

Intermediate-size noncoding RNAs (is-ncRNAs) have been shown to play important regulatory roles in the development of several eukaryotic organisms. However, they have not been thoroughly explored in Plasmodium falciparum, which is the most virulent malaria parasite infecting human being. By using Illumina/Solexa paired-end sequencing of an is-ncRNA-specific library, we performed a systematic identification of novel is-ncRNAs in intraerythrocytic P. falciparum, strain 3D7. A total of 1,198 novel is-ncRNA candidates, including antisense, intergenic, and intronic is-ncRNAs, were identified. Bioinformatics analyses showed that the intergenic is-ncRNAs were the least conserved among different Plasmodium species, and antisense is-ncRNAs were more conserved than their sense counterparts. Twenty-two novel snoRNAs were identified, and eight potential novel classes of P. falciparum is-ncRNAs were revealed by clustering analysis. The expression of randomly selected novel is-ncRNAs was confirmed by RT-PCR and northern blotting assays. An obvious different expressional profile of the novel is-ncRNA between the early and late intraerythrocytic developmental stages of the parasite was observed. The expression levels of the antisense RNAs correlated with those of their cis-encoded sense RNA counterparts, suggesting that these is-ncRNAs are involved in the regulation of gene expression of the parasite. In conclusion, we accomplished a deep profiling analysis of novel is-ncRNAs in P. falciparum, analysed the conservation and structural features of these novel is-ncRNAs, and revealed their differential expression patterns during the development of the parasite. These findings provide important information for further functional characterisation of novel is-ncRNAs during the development of P. falciparum.

The ribosomal protein rpl26 promoter is required for its 3' sense terminus ncRNA transcription in Schizosaccharomyces pombe, implicating a new transcriptional mechanism for ncRNAs

Transcriptome studies have revealed that many non-coding RNAs (ncRNAs) are located near the 3' sense terminus of protein-coding genes. However, the transcription and function of these RNAs remain elusive. Here, we identify a 3' sense termini-associated sRNA (TASR) downstream of rpl26 in Schizosaccharomyces pombe (S. pombe). Structure and function assays indicate that the TASR is an H/ACA box snoRNA required for 18S rRNA pseudouridylation at U121 and U305 sites and is therefore a cognate of snR49 from the budding yeast. Transcriptional studies show that pre-snR49 overlaps most of the coding sequence (CDS) of rpl26. Using scanning deletion analysis within promoter region, we show that the rpl26 promoter is required for the 3' TASR transcription. Interestingly, chromosomal synteny of rpl26-snR49 is found in the Schizosaccharomyces groups. Taken together, we have revealed a new transcriptional mechanism for 3' sense TASRs, which are transcribed by the same promoter as their upstream protein genes. These results further suggest that the origin and function of 3' sense ncRNAs are associated with upstream genes in higher eukaryotes.

Annotating functional RNAs in genomes using Infernal

Many different types of functional non-coding RNAs participate in a wide range of important cellular functions but the large majority of these RNAs are not routinely annotated in published genomes. Several programs have been developed for identifying RNAs, including specific tools tailored to a particular RNA family as well as more general ones designed to work for any family. Many of these tools utilize covariance models (CMs), statistical models of the conserved sequence, and structure of an RNA family. In this chapter, as an illustrative example, the Infernal software package and CMs from the Rfam database are used to identify RNAs in the genome of the archaeon Methanobrevibacter ruminantium, uncovering some additional RNAs not present in the genome's initial annotation. Analysis of the results and comparison with family-specific methods demonstrate some important strengths and weaknesses of this general approach.

Class-specific prediction of ncRNAs

Many RNA families, i.e., groups of homologous RNA genes, belong to RNA classes, such as tRNAs, snoRNAs, or microRNAs, that are characterized by common sequence motifs and/or common secondary structure features. The detection of new members of RNA classes, as well as the comprehensive annotation of genomes with members of RNA classes is a challenging task that goes beyond simple homology search. Computational methods addressing this problem typically use a three-tiered approach: In the first step an efficient and sensitive filter is employed. In the second step the candidate set is narrowed down using computationally expensive methods geared towards specificity. In the final step the hits are annotated with class-specific features and scored. Here we review the tools that are currently available for a diverse set of RNA classes.

RNA-guided isomerization of uridine to pseudouridine--pseudouridylation

Box H/ACA ribonucleoproteins (RNPs), each consisting of one unique guide RNA and 4 common core proteins, constitute a family of complex enzymes that catalyze, in an RNA-guided manner, the isomerization of uridines to pseudouridines (Ψs) in RNAs, a reaction known as pseudouridylation. Over the years, box H/ACA RNPs have been extensively studied revealing many important aspects of these RNA modifying machines. In this review, we focus on the composition, structure, and biogenesis of H/ACA RNPs. We explain the mechanism of how this enzyme family recognizes and specifies its target uridine in a substrate RNA. We discuss the substrates of box H/ACA RNPs, focusing on rRNA (rRNA) and spliceosomal small nuclear RNA (snRNA). We describe the modification product Ψ and its contribution to RNA function. Finally, we consider possible mechanisms of the bone marrow failure syndrome dyskeratosis congenita and of prostate and other cancers linked to mutations in H/ACA RNPs.

Identification and characterisation of non-coding small RNAs in the pathogenic filamentous fungus Trichophyton rubrum

BACKGROUND

Accumulating evidence demonstrates that non-coding RNAs (ncRNAs) are indispensable components of many organisms and play important roles in cellular events, regulation, and development.

RESULTS

Here, we analysed the small non-coding RNA (ncRNA) transcriptome of Trichophyton rubrum by constructing and sequencing a cDNA library from conidia and mycelia. We identified 352 ncRNAs and their corresponding genomic loci. These ncRNA candidates included 198 entirely novel ncRNAs and 154 known ncRNAs classified as snRNAs, snoRNAs and other known ncRNAs. Further bioinformatic analysis detected 96 snoRNAs, including 56 snoRNAs that had been annotated in other organisms and 40 novel snoRNAs. All snoRNAs belonged to two major classes--C/D box snoRNAs and H/ACA snoRNAs--and their potential target sites in rRNAs and snRNAs were predicted. To analyse the evolutionary conservation of the ncRNAs in T. rubrum, we aligned all 352 ncRNAs to the genomes of six dermatophytes and to the NCBI non-redundant nucleotide database (NT). The results showed that most of the identified snRNAs were conserved in dermatophytes. Of the 352 ncRNAs, 102 also had genomic loci in other dermatophytes, and 27 were dermatophyte-specific.

CONCLUSIONS

Our systematic analysis may provide important clues to the function and evolution of ncRNAs in T. rubrum. These results also provide important information to complement the current annotation of the T. rubrum genome, which primarily comprises protein-coding genes.

RNA pseudouridylation: new insights into an old modification

Pseudouridine is the most abundant post-transcriptionally modified nucleotide in various stable RNAs of all organisms. Pseudouridine is derived from uridine via base-specific isomerization, resulting in an extra hydrogen-bond donor that distinguishes it from other nucleotides. In eukaryotes, uridine-to-pseudouridine isomerization is catalyzed primarily by box H/ACA RNPs, ribonucleoproteins that act as pseudouridylases. When introduced into RNA, pseudouridine contributes significantly to RNA-mediated cellular processes. It was recently discovered that pseudouridylation can be induced by stress, suggesting a regulatory role for pseudouridine. It has also been reported that pseudouridine can be artificially introduced into mRNA by box H/ACA RNPs and that such introduction can mediate nonsense-to-sense codon conversion, thus demonstrating a new means of generating coding or protein diversity.

Role of Ess1 in growth, morphogenetic switching, and RNA polymerase II transcription in Candida albicans

Candida albicans is a fungal pathogen that causes potentially fatal infections among immune-compromised individuals. The emergence of drug resistant C. albicans strains makes it important to identify new antifungal drug targets. Among potential targets are enzymes known as peptidyl-prolyl cis/trans isomerases (PPIases) that catalyze isomerization of peptide bonds preceding proline. We are investigating a PPIase called Ess1, which is conserved in all major human pathogenic fungi. Previously, we reported that C. albicans Ess1 is essential for growth and morphogenetic switching. In the present study, we re-evaluated these findings using more rigorous genetic analyses, including the use of additional CaESS1 mutant alleles, distinct marker genes, and the engineering of suitably-matched isogenic control strains. The results confirm that CaEss1 is essential for growth in C. albicans, but show that reduction of CaESS1 gene dosage by half (δ/+) does not interfere with morphogenetic switching. However, further reduction of CaEss1 levels using a conditional allele does reduce morphogenetic switching. We also examine the role of the linker α-helix that distinguishes C. albicans Ess1 from the human Pin1 enzyme, and present results of a genome-wide transcriptome analysis. The latter analysis indicates that CaEss1 has a conserved role in regulation of RNA polymerase II function, and is required for efficient termination of small nucleolar RNAs and repression of cryptic transcription in C. albicans.

Computational analysis of noncoding RNAs

Noncoding RNAs have emerged as important key players in the cell. Understanding their surprisingly diverse range of functions is challenging for experimental and computational biology. Here, we review computational methods to analyze noncoding RNAs. The topics covered include basic and advanced techniques to predict RNA structures, annotation of noncoding RNAs in genomic data, mining RNA-seq data for novel transcripts and prediction of transcript structures, computational aspects of microRNAs, and database resources.

Comparative genomics of eukaryotic small nucleolar RNAs reveals deep evolutionary ancestry amidst ongoing intragenomic mobility

BACKGROUND

Small nucleolar (sno)RNAs are required for posttranscriptional processing and modification of ribosomal, spliceosomal and messenger RNAs. Their presence in both eukaryotes and archaea indicates that snoRNAs are evolutionarily ancient. The location of some snoRNAs within the introns of ribosomal protein genes has been suggested to belie an RNA world origin, with the exons of the earliest protein-coding genes having evolved around snoRNAs after the advent of templated protein synthesis. Alternatively, this intronic location may reflect more recent selection for coexpression of snoRNAs and ribosomal components, ensuring rRNA modification by snoRNAs during ribosome synthesis. To gain insight into the evolutionary origins of this genetic organization, we examined the antiquity of snoRNA families and the stability of their genomic location across 44 eukaryote genomes.

RESULTS

We report that dozens of snoRNA families are traceable to the Last Eukaryotic Common Ancestor (LECA), but find only weak similarities between the oldest eukaryotic snoRNAs and archaeal snoRNA-like genes. Moreover, many of these LECA snoRNAs are located within the introns of host genes independently traceable to the LECA. Comparative genomic analyses reveal the intronic location of LECA snoRNAs is not ancestral however, suggesting the pattern we observe is the result of ongoing intragenomic mobility. Analysis of human transcriptome data indicates that the primary requirement for hosting intronic snoRNAs is a broad expression profile. Consistent with ongoing mobility across broadly-expressed genes, we report a case of recent migration of a non-LECA snoRNA from the intron of a ubiquitously expressed non-LECA host gene into the introns of two LECA genes during the evolution of primates.

CONCLUSIONS

Our analyses show that snoRNAs were a well-established family of RNAs at the time when eukaryotes began to diversify. While many are intronic, this association is not evolutionarily stable across the eukaryote tree; ongoing intragenomic mobility has erased signal of their ancestral gene organization, and neither introns-first nor evolved co-expression adequately explain our results. We therefore present a third model - constrained drift - whereby individual snoRNAs are intragenomically mobile and may occupy any genomic location from which expression satisfies phenotype.

A differential sequencing-based analysis of the C. elegans noncoding transcriptome

Noncoding RNAs are increasingly being recognized as important players in eukaryote biology. However, despite major efforts in mapping the Caenorhabditis elegans transcriptome over the last couple of years, nonpolyadenylated and intermediate-size noncoding RNAs (is-ncRNAs) are still incompletely explored. We have combined an enzymatic approach with full-length RNA-Seq of is-ncRNAs in C. elegans. A total of 473 novel is-ncRNAs has been identified, of which a substantial fraction was associated with transcription factor binding sites and developmentally regulated expression patterns. Analysis of sequence and secondary structure permitted classification of more than 200 is-ncRNAs into several known RNA classes, while another 33 is-ncRNAs were identified as belonging to two previously uncharacterized groups of is-ncRNAs. Three of the unclassified is-ncRNAs contain the 5' Alu domain common to SRP RNAs and specifically bound with the SRP9/14 heterodimer in vitro. One of these (inc394) showed 65% sequence identity with the human, neuron-specific BC200 RNA. Structure-based clustering analysis and in vitro binding experiments supported the notion that the nematode stem-bulge RNAs (sbRNAs) are homologs (or functional analogs) of the Y RNAs. Moreover, analysis of the differential libraries showed that some mature snoRNAs undergo secondary 5' cap modification after processing of the primary transcript, thus suggesting the existence of a wider range of functional RNAs arising from processed and modified fragments of primary transcripts.