snoSeeker: an advanced computational package for screening of guide and orphan snoRNA genes in the human genome

TGF-β2 Induces Ribosome Activity, Alters Ribosome Composition and Inhibits IRES-Mediated Translation in Chondrocytes

Alterations in cell fate are often attributed to (epigenetic) regulation of gene expression. An emerging paradigm focuses on specialized ribosomes within a cell. However, little evidence exists for the dynamic regulation of ribosome composition and function. Here, we stimulated a chondrocytic cell line with transforming growth factor beta (TGF-β2) and mapped changes in ribosome function, composition and ribosomal RNA (rRNA) epitranscriptomics. 35S Met/Cys incorporation was used to evaluate ribosome activity. Dual luciferase reporter assays were used to assess ribosomal modus. Ribosomal RNA expression and processing were determined by RT-qPCR, while RiboMethSeq and HydraPsiSeq were used to determine rRNA modification profiles. Label-free protein quantification of total cell lysates, isolated ribosomes and secreted proteins was done by LC-MS/MS. A three-day TGF-β2 stimulation induced total protein synthesis in SW1353 chondrocytic cells and human articular chondrocytes. Specifically, TGF-β2 induced cap-mediated protein synthesis, while IRES-mediated translation was not (P53 IRES) or little affected (CrPv IGR and HCV IRES). Three rRNA post-transcriptional modifications (PTMs) were affected by TGF-β2 stimulation (18S-Gm1447 downregulated, 18S-ψ1177 and 28S-ψ4598 upregulated). Proteomic analysis of isolated ribosomes revealed increased interaction with eIF2 and tRNA ligases and decreased association of eIF4A3 and heterogeneous nuclear ribonucleoprotein (HNRNP)s. In addition, thirteen core ribosomal proteins were more present in ribosomes from TGF-β2 stimulated cells, albeit with a modest fold change. A prolonged stimulation of chondrocytic cells with TGF-β2 induced ribosome activity and changed the mode of translation. These functional changes could be coupled to alterations in accessory proteins in the ribosomal proteome.

RMBase v3.0: decode the landscape, mechanisms and functions of RNA modifications

Although over 170 chemical modifications have been identified, their prevalence, mechanism and function remain largely unknown. To enable integrated analysis of diverse RNA modification profiles, we have developed RMBase v3.0 (http://bioinformaticsscience.cn/rmbase/), a comprehensive platform consisting of eight modules. These modules facilitate the exploration of transcriptome-wide landscape, biogenesis, interactome and functions of RNA modifications. By mining thousands of epitranscriptome datasets with novel pipelines, the 'RNA Modifications' module reveals the map of 73 RNA modifications of 62 species. the 'Genes' module allows to retrieve RNA modification profiles and clusters by gene and transcript. The 'Mechanisms' module explores 23 382 enzyme-catalyzed or snoRNA-guided modified sites to elucidate their biogenesis mechanisms. The 'Co-localization' module systematically formulates potential correlations between 14 histone modifications and 6 RNA modifications in various cell-lines. The 'RMP' module investigates the differential expression profiles of 146 RNA-modifying proteins (RMPs) in 18 types of cancers. The 'Interactome' integrates the interactional relationships between 73 RNA modifications with RBP binding events, miRNA targets and SNPs. The 'Motif' illuminates the enriched motifs for 11 types of RNA modifications identified from epitranscriptome datasets. The 'Tools' introduces a novel web-based 'modGeneTool' for annotating modifications. Overall, RMBase v3.0 provides various resources and tools for studying RNA modifications.

snoRNAs in hematopoiesis and blood malignancies: A comprehensive review

Small nucleolar RNAs (snoRNAs) are noncoding RNA molecules of highly variable size, usually ranging from 60 to 150 nucleotides. They are classified into H/ACA box snoRNAs, C/D box snoRNAs, and scaRNAs. Their functional profile includes biogenesis of ribosomes, processing of rRNAs, 2'-O-methylation and pseudouridylation of RNAs, alternative splicing and processing of mRNAs and the generation of small RNA molecules like miRNA. The snoRNAs have been observed to have an important role in hematopoiesis and malignant hematopoietic conditions including leukemia, lymphoma, and multiple myeloma. Blood malignancies arise in immune system cells or the bone marrow due to chromosome abnormalities. It has been estimated that annually over 1.25 million cases of blood cancer occur worldwide. The snoRNAs often show a differential expression profile in blood malignancies. Recent reports associate the abnormal expression of snoRNAs with the inhibition of apoptosis, uncontrolled cell proliferation, angiogenesis, and metastasis. This implies that targeting snoRNAs could be a potential way to treat hematologic malignancies. In this review, we describe the various functions of snoRNAs, their role in hematopoiesis, and the consequences of their dysregulation in blood malignancies. We also evaluate the potential of the dysregulated snoRNAs as biomarkers and therapeutic targets for blood malignancies.

A Dynamic rRNA Ribomethylome Drives Stemness in Acute Myeloid Leukemia

The development and regulation of malignant self-renewal remain unresolved issues. Here, we provide biochemical, genetic, and functional evidence that dynamics in ribosomal RNA (rRNA) 2'-O-methylation regulate leukemia stem cell (LSC) activity in vivo. A comprehensive analysis of the rRNA 2'-O-methylation landscape of 94 patients with acute myeloid leukemia (AML) revealed dynamic 2'-O-methylation specifically at exterior sites of ribosomes. The rRNA 2'-O-methylation pattern is closely associated with AML development stage and LSC gene expression signature. Forced expression of the 2'-O-methyltransferase fibrillarin (FBL) induced an AML stem cell phenotype and enabled engraftment of non-LSC leukemia cells in NSG mice. Enhanced 2'-O-methylation redirected the ribosome translation program toward amino acid transporter mRNAs enriched in optimal codons and subsequently increased intracellular amino acid levels. Methylation at the single site 18S-guanosine 1447 was instrumental for LSC activity. Collectively, our work demonstrates that dynamic 2'-O-methylation at specific sites on rRNAs shifts translational preferences and controls AML LSC self-renewal.

SIGNIFICANCE

We establish the complete rRNA 2'-O-methylation landscape in human AML. Plasticity of rRNA 2'-O-methylation shifts protein translation toward an LSC phenotype. This dynamic process constitutes a novel concept of how cancers reprogram cell fate and function. This article is highlighted in the In This Issue feature, p. 247.

Non-coding RNAs in human health and disease: potential function as biomarkers and therapeutic targets

Human diseases have been a critical threat from the beginning of human history. Knowing the origin, course of action and treatment of any disease state is essential. A microscopic approach to the molecular field is a more coherent and accurate way to explore the mechanism, progression, and therapy with the introduction and evolution of technology than a macroscopic approach. Non-coding RNAs (ncRNAs) play increasingly important roles in detecting, developing, and treating all abnormalities related to physiology, pathology, genetics, epigenetics, cancer, and developmental diseases. Noncoding RNAs are becoming increasingly crucial as powerful, multipurpose regulators of all biological processes. Parallel to this, a rising amount of scientific information has revealed links between abnormal noncoding RNA expression and human disorders. Numerous non-coding transcripts with unknown functions have been found in addition to advancements in RNA-sequencing methods. Non-coding linear RNAs come in a variety of forms, including circular RNAs with a continuous closed loop (circRNA), long non-coding RNAs (lncRNA), and microRNAs (miRNA). This comprises specific information on their biogenesis, mode of action, physiological function, and significance concerning disease (such as cancer or cardiovascular diseases and others). This study review focuses on non-coding RNA as specific biomarkers and novel therapeutic targets.

Single-base resolution mapping of 2′-O-methylation sites by an exoribonuclease-enriched chemical method

2'-O-methylation (Nm) is one of the most abundant RNA epigenetic modifications and plays a vital role in the post-transcriptional regulation of gene expression. Current Nm mapping approaches are normally limited to highly abundant RNAs and have significant technical hurdles in mRNAs or relatively rare non-coding RNAs (ncRNAs). Here, we developed a new method for enriching Nm sites by using RNA exoribonuclease and periodate oxidation reactivity to eliminate 2'-hydroxylated (2'-OH) nucleosides, coupled with sequencing (Nm-REP-seq). We revealed several novel classes of Nm-containing ncRNAs as well as mRNAs in humans, mice, and drosophila. We found that some novel Nm sites are present at fixed positions in different tRNAs and are potential substrates of fibrillarin (FBL) methyltransferase mediated by snoRNAs. Importantly, we discovered, for the first time, that Nm located at the 3'-end of various types of ncRNAs and fragments derived from them. Our approach precisely redefines the genome-wide distribution of Nm and provides new technologies for functional studies of Nm-mediated gene regulation.

Small Nucleolar RNAs and Their Comprehensive Biological Functions in Hepatocellular Carcinoma

Small nucleolar RNAs (snoRNAs) are a class of highly conserved, stable non-coding RNAs involved in both post-transcriptional modification of RNA and in ribosome biogenesis. Recent research shows that the dysfunction of snoRNAs plays a pivotal role in hepatocellular carcinoma (HCC) and related etiologies, such as hepatitis B virus (HBV), hepatitis C virus (HCV), and non-alcoholic fatty liver disease (NAFLD). Growing evidence suggests that snoRNAs act as oncogenes or tumor suppressors in hepatocellular carcinoma (HCC) through multiple mechanisms. Furthermore, snoRNAs are characterized by their stability in body fluids and their clinical relevance and represent promising tools as diagnostic and prognostic biomarkers. SnoRNAs represent an emerging area of cancer research. In this review, we summarize the classification, biogenesis, activity, and functions of snoRNAs, as well as highlight the mechanism and roles of snoRNAs in HCC and related diseases. Our findings will aid in the understanding of complex processes of tumor occurrence and development, as well as suggest potential diagnostic markers and treatment targets. Furthermore, we discuss several limitations and suggest future research and application directions.

Genome-wide analysis and functional annotation of chromatin-enriched noncoding RNAs in rice during somatic cell regeneration

BACKGROUND

Plants have the remarkable ability to generate callus, a pluripotent cell mass that acquires competence for subsequent tissue regeneration. Global chromatin remodeling is required for this cell fate transition, but how the process is regulated is not fully understood. Chromatin-enriched noncoding RNAs (cheRNAs) are thought to play important roles in maintaining chromatin state. However, whether cheRNAs participate in somatic cell regeneration in plants has not yet been clarified.

RESULTS

To uncover the characteristics and functions of cheRNAs during somatic cell reprogramming in plants, we systematically investigate cheRNAs during callus induction, proliferation and regeneration in rice. We identify 2284 cheRNAs, most of which are novel long non-coding RNAs or small nucleolar RNAs. These cheRNAs, which are highly conserved across plant species, shuttle between chromatin and the nucleoplasm during somatic cell regeneration. They positively regulate the expression of neighboring genes via specific RNA motifs, which may interact with DNA motifs around cheRNA loci. Large-scale mutant analysis shows that cheRNAs are associated with plant size and seed morphology. Further detailed functional investigation of two che-lncRNAs demonstrates that their loss of function impairs cell dedifferentiation and plant regeneration, highlighting the functions of cheRNAs in regulating the expression of neighboring genes via specific motifs. These findings support cis- regulatory roles of cheRNAs in influencing a variety of rice traits.

CONCLUSIONS

cheRNAs are a distinct subclass of regulatory non-coding RNAs that are required for somatic cell regeneration and regulate rice traits. Targeting cheRNAs has great potential for crop trait improvement and breeding in future.

The role of OncoSnoRNAs and Ribosomal RNA 2'-O-methylation in Cancer

Ribosomes are essential nanomachines responsible for all protein production in cells. Ribosome biogenesis and function are energy costly processes, they are tightly regulated to match cellular needs. In cancer, major pathways that control ribosome biogenesis and function are often deregulated to ensure cell survival and to accommodate the continuous proliferation of tumour cells. Ribosomal RNAs (rRNAs) are abundantly modified with 2'-O-methylation (Nm, ribomethylation) being one of the most common modifications. In eukaryotic ribosomes, ribomethylation is performed by the methyltransferase Fibrillarin guided by box C/D small nucleolar RNAs (snoRNAs). Accumulating evidences indicate that snoRNA expression and ribosome methylation profiles are altered in cancer. Here we review our current knowledge on differential snoRNA expression and rRNA 2ʹ-O methylation in the context of human malignancies, and discuss the consequences and opportunities for cancer diagnostics, prognostics, and therapeutics.

Genome-Wide Identification and Analysis of Small Nucleolar RNAs and Their Roles in Regulating Latex Regeneration in the Rubber Tree (Hevea brasiliensis)

Small nucleolar RNAs (snoRNAs) are a class of conserved nuclear RNAs that play important roles in the modification of ribosomal RNAs (rRNAs) in plants. In rubber trees, rRNAs are run off with latex flow during tapping and need to be regenerated for maintaining the functions of the laticifer cells. SnoRNAs are expected to play essential roles in the regeneration of rRNAs. However, snoRNAs in the rubber tree have not been sufficiently characterized thus far. In this study, we performed nuclear RNA sequencing (RNA-seq) to identify snoRNAs globally and investigate their roles in latex regeneration. We identified a total of 3,626 snoRNAs by computational prediction with nuclear RNA-seq data. Among these snoRNAs, 50 were highly expressed in latex; furthermore, the results of reverse transcription polymerase chain reaction (RT-PCR) showed the abundant expression of 31 of these snoRNAs in latex. The correlation between snoRNA expression and adjusted total solid content (TSC/C) identified 13 positively yield-correlated snoRNAs. To improve the understanding of latex regeneration in rubber trees, we developed a novel insulated tapping system (ITS), which only measures the latex regenerated in specific laticifers. Using this system, a laticifer-abundant snoRNA, HbsnoR28, was found to be highly correlated with latex regeneration. To the best of our knowledge, this is the first report to globally identify snoRNAs that might be involved in latex regeneration regulation and provide new clues for unraveling the mechanisms underlying the regulation of latex regeneration.

Small nucleolar RNA is potential as a novel player in leukemogenesis and clinical application

Lack of clarity of the mechanisms that underlie leukemogenesis obstructs the diagnosis, prognosis, and treatment of leukemia. Research has found that small nuclear RNA (snoRNA) plays an essential role in leukemia. These small non-coding RNAs are involved in ribosome biogenesis, including the 2′-O-methylation and pseudouridylation of precursor ribosomal RNA (pre-rRNA), and pre-rRNA splicing. Recently, many snoRNAs were found to be orphans that have no predictable RNA modification targets, but these RNAs have always been found to be located in different subcellular organelles, and they play diverse roles. Using high-throughput technology, snoRNA expression profiles have been revealed in leukemia, and some of the deregulated snoRNAs may regulate the cell cycle, differentiation, proliferation, and apoptosis in leukemic cells and confer drug resistance during leukemia treatment. In this review, we discuss the expression profiles and functions of snoRNAs, particularly orphan snoRNAs, in leukemia. It is possible that the dysregulated snoRNAs are promising diagnosis and prognosis markers for leukemia, which may serve as potential therapeutic targets in leukemia treatment.

Emerging Classes of Small Non-Coding RNAs With Potential Implications in Diabetes and Associated Metabolic Disorders

Most of the sequences in the human genome do not code for proteins but generate thousands of non-coding RNAs (ncRNAs) with regulatory functions. High-throughput sequencing technologies and bioinformatic tools significantly expanded our knowledge about ncRNAs, highlighting their key role in gene regulatory networks, through their capacity to interact with coding and non-coding RNAs, DNAs and proteins. NcRNAs comprise diverse RNA species, including amongst others PIWI-interacting RNAs (piRNAs), involved in transposon silencing, and small nucleolar RNAs (snoRNAs), which participate in the modification of other RNAs such as ribosomal RNAs and transfer RNAs. Recently, a novel class of small ncRNAs generated from the cleavage of tRNAs or pre-tRNAs, called tRNA-derived small RNAs (tRFs) has been identified. tRFs have been suggested to regulate protein translation, RNA silencing and cell survival. While for other ncRNAs an implication in several pathologies is now well established, the potential involvement of piRNAs, snoRNAs and tRFs in human diseases, including diabetes, is only beginning to emerge. In this review, we summarize fundamental aspects of piRNAs, snoRNAs and tRFs biology. We discuss their biogenesis while emphasizing on novel sequencing technologies that allow ncRNA discovery and annotation. Moreover, we give an overview of genomic approaches to decrypt their mechanisms of action and to study their functional relevance. The review will provide a comprehensive landscape of the regulatory roles of these three types of ncRNAs in metabolic disorders by reporting their differential expression in endocrine pancreatic tissue as well as their contribution to diabetes incidence and diabetes-underlying conditions such as inflammation. Based on these discoveries we discuss the potential use of piRNAs, snoRNAs and tRFs as promising therapeutic targets in metabolic disorders.

A small nucleolar RNA, SNORD126, promotes adipogenesis in cells and rats by activating the PI3K-AKT pathway

Small nucleolar RNA (snoRNA) plays important role in various histogenesis. Whether snoRNA plays a role in adipogenesis is unknown. SNORD126 is a C/D box snoRNA. We previously demonstrated that SNORD126 promoted hepatocellular carcinoma cell growth by activating the phosphoinositide 3-kinase-protein kinase B (Akt) pathway through upregulating fibroblast growth factor receptor 2 expression. In the present study, we found that the expression of SNORD126 was downregulated in the obesity-related tissues in high-fat diet-fed rats. Overexpression of SNORD126 in 3T3-L1 cells promoted adipocytes differentiation. SNORD126 significantly increased the expression of CCAAT/enhancer-binding protein α, fatty acid-binding protein 4, peroxisome proliferative-activated receptor-γ, and the phosphorylation of Akt and p70S6K. Overexpression of SNORD126 in human adipose-derived stem cells stimulated adipogenesis and increased phosphorylation of Akt. Meanwhile, SNORD126 increased the messenger RNA and protein levels of cyclin D1 and cyclin-dependent kinase 2, which promoted mitotic clonal expansion progression during the early stage of 3T3-L1 cell differentiation. We further found that SNORD126 accelerated the growth of the groin fat pad and increased phosphorylation of Akt and p70S6K in rats. Overall, our results suggested that SNORD126 promoted adipocyte differentiation through increasing phosphorylation of Akt and p70S6K both in vitro and in vivo.

Conservation of Archaeal C/D Box sRNA-Guided RNA Modifications

Post-transcriptional modifications fulfill many important roles during ribosomal RNA maturation in all three domains of life. Ribose 2'-O-methylations constitute the most abundant chemical rRNA modification and are, for example, involved in RNA folding and stabilization. In archaea, these modification sites are determined by variable sets of C/D box sRNAs that guide the activity of the rRNA 2'-O-methyltransferase fibrillarin. Each C/D box sRNA contains two guide sequences that can act in coordination to bridge rRNA sequences. Here, we will review the landscape of archaeal C/D box sRNA genes and their target sites. One focus is placed on the apparent accelerated evolution of guide sequences and the varied pairing of the two individual guides, which results in different rRNA modification patterns and RNA chaperone activities.

Small nucleolar RNAs: continuing identification of novel members and increasing diversity of their molecular mechanisms of action

Identified five decades ago amongst the most abundant cellular RNAs, small nucleolar RNAs (snoRNAs) were initially described as serving as guides for the methylation and pseudouridylation of ribosomal RNA through direct base pairing. In recent years, however, increasingly powerful high-throughput genomic approaches and strategies have led to the discovery of many new members of the family and surprising diversity in snoRNA functionality and mechanisms of action. SnoRNAs are now known to target RNAs of many biotypes for a wider range of modifications, interact with diverse binding partners, compete with other binders for functional interactions, recruit diverse players to targets and affect protein function and accessibility through direct interaction. This mini-review presents the continuing characterization of the snoRNome through the identification of new snoRNA members and the discovery of their mechanisms of action, revealing a highly versatile noncoding family playing central regulatory roles and connecting the main cellular processes.

SNORD126 Promotes Hepatitis C Virus Infection by Upregulating Claudin-1 via Activation of PI3K-AKT Signaling Pathway

Hepatitis C virus (HCV) infection involves a variety of viral and host factors, some of which promote the infection process. A small nucleolar RNA, C/D box 126 (SNORD126), was previously shown to be associated with hepatocellular carcinoma (HCC). However, the role of SNORD126 in HCV infection, which is one of the primary reasons for HCC development, has not been elucidated. In the present study, using small nucleolar RNA profiling, we observed that SNORD126 was significantly downregulated during HCV infection in both Huh7 and Huh7.5.1 cells. In addition, overexpression of SNORD126 enhanced HCV entry into host cells, whereas SNORD126 knockdown showed the opposite effect, suggesting that SNORD126 promotes HCV infection, especially through viral entry. Further functional analysis revealed that SNORD126 could enhance the expression level of claudin-1 (CLDN1), a key HCV entry factor, by increasing the levels of phosphorylated AKT. Additionally, the function of SNORD126 in HCV infection was associated with ribonucleoprotein (RNP) complexes. In summary, our findings demonstrate that oncogenic SNORD126 levels are decreased during HCV infection probably due to the host defense reaction, and SNORD126 may be important to promote viral entry by increasing CLDN1 expression through activation of the PI3K-AKT pathway, the mechanism of which is partly associated with SNORD126-mediated snoRNA RNP (snoRNP) function. Our work here provides initial evidence that endogenous snoRNA takes part in HCV infection and shows potential as a diagnostic or antiviral agent.

Functional diversity of small nucleolar RNAs

Small nucleolar RNAs (snoRNAs) are short non-protein-coding RNAs with a long-recognized role in tuning ribosomal and spliceosomal function by guiding ribose methylation and pseudouridylation at targeted nucleotide residues of ribosomal and small nuclear RNAs, respectively. SnoRNAs are increasingly being implicated in regulation of new types of post-transcriptional processes, for example rRNA acetylation, modulation of splicing patterns, control of mRNA abundance and translational efficiency, or they themselves are processed to shorter stable RNA species that seem to be the principal or alternative bioactive isoform. Intriguingly, some display unusual cellular localization under exogenous stimuli, or tissue-specific distribution. Here, we discuss the new and unforeseen roles attributed to snoRNAs, focusing on the presumed mechanisms of action. Furthermore, we review the experimental approaches to study snoRNA function, including high resolution RNA:protein and RNA:RNA interaction mapping, techniques for analyzing modifications on targeted RNAs, and cellular and animal models used in snoRNA biology research.

ESDA: An Improved Approach to Accurately Identify Human snoRNAs for Precision Cancer Therapy

Background:

SnoRNAs (Small nucleolar RNAs) are small RNA molecules with approximately 60-300 nucleotides in sequence length. They have been proved to play important roles in cancer occurrence and progression. It is of great clinical importance to identify new snoRNAs as fast and accurately as possible.

Objective:

A novel algorithm, ESDA (Elastically Sparse Partial Least Squares Discriminant Analysis), was proposed to improve the speed and the performance of recognizing snoRNAs from other RNAs in human genomes.

Methods:

In ESDA algorithm, to optimize the extracted information, kernel features were selected from the variables extracted from both primary sequences and secondary structures. Then they were used by SPLSDA (sparse partial least squares discriminant analysis) algorithm as input variables for the final classification model training to distinguish snoRNA sequences from other Human RNAs. Due to the fact that no prior biological knowledge is request to optimize the classification model, ESDA is a very practical method especially for completely new sequences.

Results:

89 H/ACA snoRNAs and 269 C/D snoRNAs of human were used as positive samples and 3403 non-snoRNAs as negative samples to test the identification performance of the proposed ESDA. For the H/ACA snoRNAs identification, the sensitivity and specificity were respectively as high as 99.6% and 98.8%. For C/D snoRNAs, they were respectively 96.1% and 98.3%. Furthermore, we compared ESDA with other widely used algorithms and classifiers: SnoReport, RF (Random Forest), DWD (Distance Weighted Discrimination) and SVM (Support Vector Machine). The highest improvement of accuracy obtained by ESDA was 25.1%.

Conclusion:

Strongly proved the superiority performance of ESDA and make it promising for identifying SnoRNAs for further development of the precision medicine for cancers.

"Lost and Found": snoRNA Annotation in the Xenopus Genome and Implications for Evolutionary Studies

Small nucleolar RNAs (snoRNAs) function primarily as guide RNAs for posttranscriptional modification of rRNAs and spliceosomal snRNAs, both of which are functionally important and evolutionarily conserved molecules. It is commonly believed that snoRNAs and the modifications they mediate are highly conserved across species. However, most relevant data on snoRNA annotation and RNA modification are limited to studies on human and yeast. Here, we used RNA-sequencing data from the giant oocyte nucleus of the frog Xenopus tropicalis to annotate a nearly complete set of snoRNAs. We compared the frog data with snoRNA sets from human and other vertebrate genomes, including mammals, birds, reptiles, and fish. We identified many Xenopus-specific (or nonhuman) snoRNAs and Xenopus-specific domains in snoRNAs from conserved RNA families. We predicted that some of these nonhuman snoRNAs and domains mediate modifications at unexpected positions in rRNAs and snRNAs. These modifications were mapped as predicted when RNA modification assays were applied to RNA from nine vertebrate species: frogs X. tropicalis and X. laevis, newt Notophthalmus viridescens, axolotl Ambystoma mexicanum, whiptail lizard Aspidoscelis neomexicana, zebrafish Danio rerio, chicken, mouse, and human. This analysis revealed that only a subset of RNA modifications is evolutionarily conserved and that modification patterns may vary even between closely related species. We speculate that each functional domain in snoRNAs (half of an snoRNA) may evolve independently and shuffle between different snoRNAs.

Developmental changes of rRNA ribose methylations in the mouse

A sequencing-based profiling method (RiboMeth-seq) for ribose methylations was used to study methylation patterns in mouse adult tissues and during development. In contrast to previous reports based on studies of human cancer cell lines, almost all methylation sites were close to fully methylated in adult tissues. A subset of sites was differentially modified in developing tissues compared to their adult counterparts and showed clear developmental dynamics. This provides the first evidence for ribosome heterogeneity at the level of rRNA modifications during mouse development. In a prominent example, the expression levels of SNORD78 during development appeared to be regulated by alternative splicing of the Gas5 host-gene and to correlate with the methylation level of its target site at LSU-G4593. The results are discussed in the context of the specialized ribosome hypothesis.

Analysis of Expression Pattern of snoRNAs in Different Cancer Types with Machine Learning Algorithms

Small nucleolar RNAs (snoRNAs) are a new type of functional small RNAs involved in the chemical modifications of rRNAs, tRNAs, and small nuclear RNAs. It is reported that they play important roles in tumorigenesis via various regulatory modes. snoRNAs can both participate in the regulation of methylation and pseudouridylation and regulate the expression pattern of their host genes. This research investigated the expression pattern of snoRNAs in eight major cancer types in TCGA via several machine learning algorithms. The expression levels of snoRNAs were first analyzed by a powerful feature selection method, Monte Carlo feature selection (MCFS). A feature list and some informative features were accessed. Then, the incremental feature selection (IFS) was applied to the feature list to extract optimal features/snoRNAs, which can make the support vector machine (SVM) yield best performance. The discriminative snoRNAs included HBII-52-14, HBII-336, SNORD123, HBII-85-29, HBII-420, U3, HBI-43, SNORD116, SNORA73B, SCARNA4, HBII-85-20, etc., on which the SVM can provide a Matthew’s correlation coefficient (MCC) of 0.881 for predicting these eight cancer types. On the other hand, the informative features were fed into the Johnson reducer and repeated incremental pruning to produce error reduction (RIPPER) algorithms to generate classification rules, which can clearly show different snoRNAs expression patterns in different cancer types. The analysis results indicated that extracted discriminative snoRNAs can be important for identifying cancer samples in different types and the expression pattern of snoRNAs in different cancer types can be partly uncovered by quantitative recognition rules.

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.

Diverse RNA interference strategies in early-branching metazoans

BACKGROUND

Micro RNAs (miRNAs) and piwi interacting RNAs (piRNAs), along with the more ancient eukaryotic endogenous small interfering RNAs (endo-siRNAs) constitute the principal components of the RNA interference (RNAi) repertoire of most animals. RNAi in non-bilaterians - sponges, ctenophores, placozoans and cnidarians - appears to be more diverse than that of bilaterians, and includes structurally variable miRNAs in sponges, an enormous number of piRNAs in cnidarians and the absence of miRNAs in ctenophores and placozoans.

RESULTS

Here we identify thousands of endo-siRNAs and piRNAs from the sponge Amphimedon queenslandica, the ctenophore Mnemiopsis leidyi and the cnidarian Nematostella vectensis using a computational approach that clusters mapped small RNA sequences and annotates each cluster based on the read length and relative abundance of the constituent reads. This approach was validated on 11 small RNA libraries in Drosophila melanogaster, demonstrating the successful annotation of RNAi-associated loci with properties consistent with previous reports. In the non-bilaterians we uncover seven new miRNAs from Amphimedon and four from Nematostella as well as sub-populations of candidate cis-natural antisense transcript (cis-NAT) endo-siRNAs. We confirmed the absence of miRNAs in Mnemiopsis but detected an abundance of endo-siRNAs in this ctenophore. Analysis of putative piRNA structure suggests that conserved localised secondary structures in primary transcripts may be important for the production of mature piRNAs in Amphimedon and Nematostella, as is also the case for endo-siRNAs.

CONCLUSION

Together, these findings suggest that the last common ancestor of extant animals did not have the entrained RNAi system that typifies bilaterians. Instead it appears that bilaterians, cnidarians, ctenophores and sponges express unique repertoires and combinations of miRNAs, piRNAs and endo-siRNAs.

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.

Nuclear fate of yeast snoRNA is determined by co-transcriptional Rnt1 cleavage

Small nucleolar RNA (snoRNA) are conserved and essential non-coding RNA that are transcribed by RNA Polymerase II (Pol II). Two snoRNA classes, formerly distinguished by their structure and ribonucleoprotein composition, act as guide RNA to target RNA such as ribosomal RNA, and thereby introduce specific modifications. We have studied the 5'end processing of individually transcribed snoRNA in S. cerevisiae to define their role in snoRNA biogenesis and functionality. Here we show that pre-snoRNA processing by the endonuclease Rnt1 occurs co-transcriptionally with removal of the m7G cap facilitating the formation of box C/D snoRNA. Failure of this process causes aberrant 3'end processing and mislocalization of snoRNA to the cytoplasm. Consequently, Rnt1-dependent 5'end processing of box C/D snoRNA is critical for snoRNA-dependent methylation of ribosomal RNA. Our results reveal that the 5'end processing of box C/D snoRNA defines their distinct pathway of maturation.

SNORD126 promotes HCC and CRC cell growth by activating the PI3K-AKT pathway through FGFR2

Small nucleolar RNA (snoRNA) dysfunctions have been associated with cancer development. SNORD126 is an orphan C/D box snoRNA that is encoded within introns 5-6 of its host gene, cyclin B1-interacting protein 1 (CCNB1IP1). The cancer-associated molecular mechanisms triggered by SNORD126 are not fully understood. Here, we demonstrate that SNORD126 is highly expressed in hepatocellular carcinoma (HCC) and colorectal cancer (CRC) patient samples. SNORD126 increased Huh-7 and SW480 cell growth and tumorigenicity in nude mice. Knockdown of SNORD126 inhibited HepG2 and LS174T cell growth. We verified that SNORD126 was not processed into small RNAs with miRNA activity. Moreover, SNORD126 did not show a significant expression correlation with CCNB1IP1 in HCC samples or regulate CCNB1IP1 expression. Our gene expression profile analysis indicated that SNORD126-upregulated genes frequently mapped to the PI3K-AKT pathway. SNORD126 overexpression increased the levels of phosphorylated AKT, GSK-3β, and p70S6K and elevated fibroblast growth factor receptor 2 (FGFR2) expression. siRNA-mediated knockdown or AZD4547-mediated inactivation of FGFR2 in SNORD126-overexpressing Huh-7 cells inhibited AKT phosphorylation and suppressed cell growth. These findings indicate an oncogenic role for SNORD126 in cancer and suggest its potential as a therapeutic target.

Clinical significance and biological roles of small nucleolar RNAs in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is among the most common and fatal cancers. It is a multistage and multifactorial carcinoma, in which a number of factors serve roles in its initiation and progression. Small nucleolar RNAs (snoRNAs), considered to serve a role in various cancers, have recently been identified to have significant contributions to HCC tumorigenesis. Recent studies suggest that snoRNAs have a critical role in the pathogenesis of HCC. Moreover, detailed studies have demonstrated that various snoRNAs are involved in a range of biological processes associated with HCC, including initiation, proliferation, tumor growth, the cell cycle, apoptosis and metastasis. In the present review, an overview of recent studies to date has been provided, focusing on the association of snoRNAs with HCC. Based on the findings, further studies focusing on the association of snoRNAs with HCC are required to verify the diagnostic and therapeutic capacities of snoRNAs in HCC.

Using Baseline Transcriptional Connectomes in Rat to Identify Genetic Pathways Associated with Predisposition to Complex Traits

Although rat is a critical model organism in preclinical medications development, its use in systems genetics studies remains sparse. The PhenoGen database and website contain detailed information on the qualitative and quantitative aspects of the rat brain, liver, heart, and brown adipose transcriptome. This database has been generated using the HXB/BXH recombinant inbred panel and is being expanded to a hybrid rat diversity panel that includes many common inbred strains as well. By using such a panel, the PhenoGen project has created a renewable and cumulative resource for the rat genomics community. The database has been used to reconstruct the brain transcriptome identifying both annotated and unannotated transcribed elements that range in size from 20 nucleotides to over 30,000 nucleotides and elements that have a wide variety of roles in the cell including generation of proteins and regulation of the transcription and translation processes. In all 4 tissues, baseline transcriptional connectomes have been generated to model the relationships among transcripts. These connectomes can be used to identify genetic pathways associated with complex traits and to gain insight into biological function of individual transcripts. The PhenoGen website contains tools that allow the user to explore qualitative features of individual genes and to see how the gene relates to other genes within a tissue. The PhenoGen database and website continue to grow and to make use of the latest statistical methods for systems genetics creating a national resource for the rat genomics community.

Fungal Genome Annotation

The term "genome annotation" includes identification of protein-coding and noncoding sequences (e.g., repeats, rDNA, and ncRNA) in genome assemblies and attaching functional information (metadata) to these annotated features. Here, we describe the basic outline of fungal nuclear and mitochondrial genome annotation as performed at the US Department of Energy Joint Genome Institute (JGI).

Regulatory RNPs: a novel class of ribonucleoproteins that potentially contribute to ribosome heterogeneity

Many ribonucleoproteins (RNPs), which are comprised of noncoding RNA and associated proteins, are involved in essential cellular processes such as translation and pre-mRNA splicing. One class of RNP is the small Cajal body-specific RNP (scaRNP), which contributes to the biogenesis of small nuclear RNPs (snRNPs) that are central components of the spliceosome. Three scaRNAs are internally processed, generating stable nucleolus-enriched RNAs of unknown function. Here, we provide data that show that these RNAs become part of RNPs we term regulatory RNPs (regRNPs). Most modifications within rRNA (predominantly pseudouridylation and ribose 2′-O-methylation) are conducted by small nucleolar RNPs (snoRNPs), and we provide evidence that the activity of at least some of these snoRNPs is under the control of regRNPs. Because modifications within rRNA can vary in different physiological or pathological situations, rRNA modifications are thought to be the major source of ribosome heterogeneity. Our identification of regRNPs thus provides a potential mechanism for how ribosome heterogeneity may be accomplished. This work also provides additional functional connections between the Cajal body and the nucleolus.

Summary: Processed scaRNAs give rise to a novel regulatory RNP, which regulates the modification of ribosomal RNA. These findings provide insight into the mechanisms governing ribosome heterogeneity.

Protein coding genes as hosts for noncoding RNA expression

With the emergence of high-throughput sequence characterization methods and the subsequent improvements in gene annotations, it is becoming increasingly clear that a large proportion of eukaryotic protein-coding genes (as many as 50% in human) serve as host genes for non-coding RNA genes. Amongst the most extensively characterized embedded non-coding RNA genes, small nucleolar RNAs and microRNAs represent abundant families. Encoded individually or clustered, in sense or antisense orientation with respect to their host and independently expressed or dependent on host expression, the genomic characteristics of embedded genes determine their biogenesis and the extent of their relationship with their host gene. Not only can host genes and the embedded genes they harbour be co-regulated and mutually modulate each other, many are functionally coupled playing a role in the same cellular pathways. And while host-non-coding RNA relationships can be highly conserved, mechanisms have been identified, and in particular an association with transposable elements, allowing the appearance of copies of non-coding genes nested in host genes, or the migration of embedded genes from one host gene to another. The study of embedded non-coding genes and their relationship with their host genes increases the complexity of cellular networks and provides important new regulatory links that are essential to properly understand cell function.

Identification of urinary exosomal noncoding RNAs as novel biomarkers in chronic kidney disease

In chronic kidney disease (CKD), the decline in the glomerular filtration rate is associated with increased morbidity and mortality and thus poses a major challenge for healthcare systems. While the contribution of tissue-derived miRNAs and mRNAs to CKD progression has been extensively studied, little is known about the role of urinary exosomes and their association with CKD. Exosomes are small, membrane-derived endocytic vesicles that contribute to cell-to-cell communication and are present in various body fluids, such as blood or urine. Next-generation sequencing approaches have revealed that exosomes are enriched in noncoding RNAs and thus exhibit great potential for sensitive nucleic acid biomarkers in various human diseases. Therefore, in this study we aimed to identify urinary exosomal ncRNAs as novel biomarkers for diagnosis of CKD. Since up to now most approaches have focused on the class of miRNAs, we extended our analysis to several other noncoding RNA classes, such as tRNAs, tRNA fragments (tRFs), mitochondrial tRNAs, or lincRNAs. For their computational identification from RNA-seq data, we developed a novel computational pipeline, designated as ncRNASeqScan. By these analyses, in CKD patients we identified 30 differentially expressed ncRNAs, derived from urinary exosomes, as suitable biomarkers for early diagnosis. Thereby, miRNA-181a appeared as the most robust and stable potential biomarker, being significantly decreased by about 200-fold in exosomes of CKD patients compared to healthy controls. Using a cell culture system for CKD indicated that urinary exosomes might indeed originate from renal proximal tubular epithelial cells.

SnoReport 2.0: new features and a refined Support Vector Machine to improve snoRNA identification

Background

snoReport uses RNA secondary structure prediction combined with machine learning as the basis to identify the two main classes of small nucleolar RNAs, the box H/ACA snoRNAs and the box C/D snoRNAs. Here, we present snoReport 2.0, which substantially improves and extends in the original method by: extracting new features for both box C/D and H/ACA box snoRNAs; developing a more sophisticated technique in the SVM training phase with recent data from vertebrate organisms and a careful choice of the SVM parameters C and γ; and using updated versions of tools and databases used for the construction of the original version of snoReport. To validate the new version and to demonstrate its improved performance, we tested snoReport 2.0 in different organisms.

Results

Results of the training and test phases of boxes H/ACA and C/D snoRNAs, in both versions of snoReport, are discussed. Validation on real data was performed to evaluate the predictions of snoReport 2.0. Our program was applied to a set of previously annotated sequences, some of them experimentally confirmed, of humans, nematodes, drosophilids, platypus, chickens and leishmania. We significantly improved the predictions for vertebrates, since the training phase used information of these organisms, but H/ACA box snoRNAs identification was improved for the other ones.

Conclusion

We presented snoReport 2.0, to predict H/ACA box and C/D box snoRNAs, an efficient method to find true positives and avoid false positives in vertebrate organisms. H/ACA box snoRNA classifier showed an F-score of 93 % (an improvement of 10 % regarding the previous version), while C/D box snoRNA classifier, an F-Score of 94 % (improvement of 14 %). Besides, both classifiers exhibited performance measures above 90 %. These results show that snoReport 2.0 avoid false positives and false negatives, allowing to predict snoRNAs with high quality. In the validation phase, snoReport 2.0 predicted 67.43 % of vertebrate organisms for both classes. For Nematodes and Drosophilids, 69 % and 76.67 %, for H/ACA box snoRNAs were predicted, respectively, showing that snoReport 2.0 is good to identify snoRNAs in vertebrates and also H/ACA box snoRNAs in invertebrates organisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1345-6) contains supplementary material, which is available to authorized users.

An updated human snoRNAome

Small nucleolar RNAs (snoRNAs) are a class of non-coding RNAs that guide the post-transcriptional processing of other non-coding RNAs (mostly ribosomal RNAs), but have also been implicated in processes ranging from microRNA-dependent gene silencing to alternative splicing. In order to construct an up-to-date catalog of human snoRNAs we have combined data from various databases, de novo prediction and extensive literature review. In total, we list more than 750 curated genomic loci that give rise to snoRNA and snoRNA-like genes. Utilizing small RNA-seq data from the ENCODE project, our study characterizes the plasticity of snoRNA expression identifying both constitutively as well as cell type specific expressed snoRNAs. Especially, the comparison of malignant to non-malignant tissues and cell types shows a dramatic perturbation of the snoRNA expression profile. Finally, we developed a high-throughput variant of the reverse-transcriptase-based method for identifying 2'-O-methyl modifications in RNAs termed RimSeq. Using the data from this and other high-throughput protocols together with previously reported modification sites and state-of-the-art target prediction methods we re-estimate the snoRNA target RNA interaction network. Our current results assign a reliable modification site to 83% of the canonical snoRNAs, leaving only 76 snoRNA sequences as orphan.

RMBase: a resource for decoding the landscape of RNA modifications from high-throughput sequencing data

Although more than 100 different types of RNA modifications have been characterized across all living organisms, surprisingly little is known about the modified positions and their functions. Recently, various high-throughput modification sequencing methods have been developed to identify diverse post-transcriptional modifications of RNA molecules. In this study, we developed a novel resource, RMBase (RNA Modification Base, http://mirlab.sysu.edu.cn/rmbase/), to decode the genome-wide landscape of RNA modifications identified from high-throughput modification data generated by 18 independent studies. The current release of RMBase includes ∼ 9500 pseudouridine (Ψ) modifications generated from Pseudo-seq and CeU-seq sequencing data, ∼ 1000 5-methylcytosines (m(5)C) predicted from Aza-IP data, ∼ 124 200 N6-Methyladenosine (m(6)A) modifications discovered from m(6)A-seq and ∼ 1210 2'-O-methylations (2'-O-Me) identified from RiboMeth-seq data and public resources. Moreover, RMBase provides a comprehensive listing of other experimentally supported types of RNA modifications by integrating various resources. It provides web interfaces to show thousands of relationships between RNA modification sites and microRNA target sites. It can also be used to illustrate the disease-related SNPs residing in the modification sites/regions. RMBase provides a genome browser and a web-based modTool to query, annotate and visualize various RNA modifications. This database will help expand our understanding of potential functions of RNA modifications.

Functional and Structural Characterization of FAU Gene/Protein from Marine Sponge Suberites domuncula

Finkel-Biskis-Reilly murine sarcoma virus (FBR-MuSV) ubiquitously expressed (FAU) gene is down-regulated in human prostate, breast and ovarian cancers. Moreover, its dysregulation is associated with poor prognosis in breast cancer. Sponges (Porifera) are animals without tissues which branched off first from the common ancestor of all metazoans. A large majority of genes implicated in human cancers have their homologues in the sponge genome. Our study suggests that FAU gene from the sponge Suberites domuncula reflects characteristics of the FAU gene from the metazoan ancestor, which have changed only slightly during the course of animal evolution. We found pro-apoptotic activity of sponge FAU protein. The same as its human homologue, sponge FAU increases apoptosis in human HEK293T cells. This indicates that the biological functions of FAU, usually associated with "higher" metazoans, particularly in cancer etiology, possess a biochemical background established early in metazoan evolution. The ancestor of all animals possibly possessed FAU protein with the structure and function similar to evolutionarily more recent versions of the protein, even before the appearance of true tissues and the origin of tumors and metastasis. It provides an opportunity to use pre-bilaterian animals as a simpler model for studying complex interactions in human cancerogenesis.

The emerging landscape of small nucleolar RNAs in cell biology

Small nucleolar RNAs (snoRNAs) are a large class of small noncoding RNAs present in all eukaryotes sequenced thus far. As a family, they have been well characterized as playing a central role in ribosome biogenesis, guiding either the sequence-specific chemical modification of pre-rRNA (ribosomal RNA) or its processing. However, in higher eukaryotes, numerous orphan snoRNAs were described over a decade ago, with no known target or ascribed function, suggesting the possibility of alternative cellular functionality. In recent years, thanks in great part to advances in sequencing methodologies, we have seen many examples of the diversity that exists in the snoRNA family on multiple levels. In this review, we discuss the identification of novel snoRNA members, of unexpected binding partners, as well as the clarification and extension of the snoRNA target space and the characterization of diverse new noncanonical functions, painting a new and extended picture of the snoRNA landscape. Under the deluge of novel features and functions that have recently come to light, snoRNAs emerge as a central, dynamic, and highly versatile group of small regulatory RNAs. WIREs RNA 2015, 6:381–397. doi: 10.1002/wrna.1284

Genome-wide small nucleolar RNA expression analysis of lung cancer by next-generation deep sequencing

Emerging evidence indicates that small nucleolar RNAs (snoRNAs), a class of small noncoding RNAs, may play important function in tumorigenesis. Nonsmall-cell lung cancer (NSCLC) is the number one cancer killer for men and women. Systematically characterizing snoRNAs in NSCLC will develop biomarkers for its early detection and prognostication. We used next-generation deep sequencing to comprehensively characterize snoRNA profiles in 12 NSCLC tissues. We used quantitative reverse transcription polymerase chain reaction (qRT-PCR) to verify the findings in 40 surgical Stage I NSCLC specimens and 126 frozen NSCLC tissues of different stages. The 126 NSCLC tissues were divided into a training set and a testing set. Deep sequencing identified 458 snoRNAs, of which, 29 had a ≥3.0-fold expression level change in Stage I NSCLC tissues versus normal tissues. qRT-PCR analysis showed that 16 of 29 snoRNAs exhibited consistent changes with deep sequencing data. The 16 snoRNAs exhibited 0.75-0.94 area under receiver-operator characteristic curve values in distinguishing lung tumor from normal lung tissues (all ≤0.0001) with 70.0-95.0% sensitivity and 70.0-95.0% specificity. Six genes (snoRA47, snoRA68, snoRA78, snoRA21, snoRD28 and snoRD66) were identified whose expressions were associated with overall survival of the NSCLC patients. A prediction model consisting of three genes (snoRA47, snoRA68 and snoRA78) was developed in the training set of 77 cases, which could significantly predict overall survival of the NSCLC patients (p < 0.0001). The prognostic performance of the prediction model was confirmed in the testing set of 49 NSCLC patients. The identified snoRNA signatures may provide potential biomarkers for the early detection and prognostication of NSCLC.

Rapid birth-and-death evolution of imprinted snoRNAs in the Prader-Willi syndrome locus: implications for neural development in Euarchontoglires

Imprinted small nucleolar RNAs (snoRNAs) are only found in eutherian genomes and closely related to brain functions. A complex human neurological disease, Prader-Willi syndrome (PWS), is primarily attributed to the deletion of imprinted snoRNAs in chromosome 15q11-q13. Here we investigated the snoRNA repertoires in the PWS locus of 12 mammalian genomes and their evolution processes. A total of 613 imprinted snoRNAs were identified in the PWS homologous loci and the gene number was highly variable across lineages, with a peak in Euarchontoglires. Lineage-specific gene gain and loss events account for most extant genes of the HBII-52 (SNORD115) and the HBII-85 (SNORD116) gene family, and remarkable high gene-birth rates were observed in the primates and the rodents. Meanwhile, rapid sequence substitution occurred only in imprinted snoRNA genes, rather than their flanking sequences or the protein-coding genes located in the same imprinted locus. Strong selective constraints on the functional elements of these imprinted snoRNAs further suggest that they are subjected to birth-and-death evolution. Our data suggest that the regulatory role of HBII-52 on 5-HT_2CR pre-mRNA might originate in the Euarchontoglires through adaptive process. We propose that the rapid evolution of PWS-related imprinted snoRNAs has contributed to the neural development of Euarchontoglires.

Small nucleolar RNAs and RNA-guided post-transcriptional modification

snoRNAs (small nucleolar RNAs) constitute one of the largest and best-studied classes of non-coding RNAs that confer enzymatic specificity. With associated proteins, these snoRNAs form ribonucleoprotein complexes that can direct 2'-O-methylation or pseudouridylation of target non-coding RNAs. Aided by computational methods and high-throughput sequencing, new studies have expanded the diversity of known snoRNA functions. Complexes incorporating snoRNAs have dynamic specificity, and include diverse roles in RNA silencing, telomerase maintenance and regulation of alternative splicing. Evidence that dysregulation of snoRNAs can cause human disease, including cancer, indicates that the full scope of snoRNA roles remains an unfinished story. The diversity in structure, genomic origin and function between snoRNAs found in different complexes and among different phyla illustrates the surprising plasticity of snoRNAs in evolution. The ability of snoRNAs to direct highly specific interactions with other RNAs is a consistent thread in their newly discovered functions. Because they are ubiquitous throughout Eukarya and Archaea, it is likely they were a feature of the last common ancestor of these two domains, placing their origin over two billion years ago. In the present chapter, we focus on recent advances in our understanding of these ancient, but functionally dynamic RNA-processing machines.

Improved annotation of 3' untranslated regions and complex loci by combination of strand-specific direct RNA sequencing, RNA-Seq and ESTs

The reference annotations made for a genome sequence provide the framework for all subsequent analyses of the genome. Correct and complete annotation in addition to the underlying genomic sequence is particularly important when interpreting the results of RNA-seq experiments where short sequence reads are mapped against the genome and assigned to genes according to the annotation. Inconsistencies in annotations between the reference and the experimental system can lead to incorrect interpretation of the effect on RNA expression of an experimental treatment or mutation in the system under study. Until recently, the genome-wide annotation of 3′ untranslated regions received less attention than coding regions and the delineation of intron/exon boundaries. In this paper, data produced for samples in Human, Chicken and A. thaliana by the novel single-molecule, strand-specific, Direct RNA Sequencing technology from Helicos Biosciences which locates 3′ polyadenylation sites to within +/− 2 nt, were combined with archival EST and RNA-Seq data. Nine examples are illustrated where this combination of data allowed: (1) gene and 3′ UTR re-annotation (including extension of one 3′ UTR by 5.9 kb); (2) disentangling of gene expression in complex regions; (3) clearer interpretation of small RNA expression and (4) identification of novel genes. While the specific examples displayed here may become obsolete as genome sequences and their annotations are refined, the principles laid out in this paper will be of general use both to those annotating genomes and those seeking to interpret existing publically available annotations in the context of their own experimental data.

Computational approaches in detecting non- coding RNA

The important role of non coding RNAs (ncRNAs) in the cell has made their identification a critical issue in the biological research. However, traditional approaches such as PT-PCR and Northern Blot are costly. With recent progress in bioinformatics and computational prediction technology, the discovery of ncRNAs has become realistically possible. This paper aims to introduce major computational approaches in the identification of ncRNAs, including homologous search, de novo prediction and mining in deep sequencing data. Furthermore, related software tools have been compared and reviewed along with a discussion on future improvements.

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.

Combined miRNA and mRNA signature identifies key molecular players and pathways involved in chikungunya virus infection in human cells

Since its discovery, Chikungunya fever caused by a virus (CHIKV) has ravaged most of Africa and Southeast Asia. Despite there being more than a million reported cases in India alone and the seriousness of the disease in the chronic phase, a clear understanding of the disease pathogenesis and host response remains elusive. Here, we use microarray technology and quantitative PCR method to establish the complete miRNA, snoRNA and mRNA signature of host response upon CHIKV infection in human cell line infection model, HEK293T. The results were further validated in human primary cells (dermal fibroblasts). miRNA expression profiling revealed regulation of 152 miRNAs post CHIKV infection. An interesting overlap in miRNA signature was seen majorly with HCV, HPV and HIV1 virus. The microarray data further validated by qRT-PCR revealed induction of miR-744, miR-638, miR-503 and others among the top upregulated miRNAs. Notably, we found induction of snoRNAs belonging to C/D cluster including close paralogs of U3, U44, U76 and U78 snoRNAs. Genes were found to be differentially expressed along 3 major pathways; TGF-β, endocytosis and the cell cycle pathways. qRT-PCR data confirmed strong induction of TGF-β (SMAD6, JUN, SKIL) and endocytosis pathway (CXCR4, HSPA8, ADRB1) genes while downregulation of cell cycle genes (CDC27 and CDC23). Interestingly, use of TGF-β inhibitor, SB-431542, increased CHIKV mediated cell death. Overall, this study aims at providing the first complete transcriptome signature of host response upon CHIKV infection to aid identification of possible biomarkers and therapeutic targets.

snoStrip: a snoRNA annotation pipeline

MOTIVATION

Although small nucleolar RNAs form an important class of non-coding RNAs, no comprehensive annotation efforts have been undertaken, presumably because the task is complicated by both the large number of distinct small nucleolar RNA families and their relatively rapid pace of sequence evolution.

RESULTS

With snoStrip we present an automatic annotation pipeline developed specifically for comparative genomics of small nucleolar RNAs. It makes use of sequence conservation, canonical box motifs as well as secondary structure and predicts putative targets.

AVAILABILITY AND IMPLEMENTATION

The snoStrip web service and the download version is available at http://snostrip.bioinf.uni-leipzig.de/

Insights into snoRNA biogenesis and processing from PAR-CLIP of snoRNA core proteins and small RNA sequencing

BACKGROUND

In recent years, a variety of small RNAs derived from other RNAs with well-known functions such as tRNAs and snoRNAs, have been identified. The functional relevance of these RNAs is largely unknown. To gain insight into the complexity of snoRNA processing and the functional relevance of snoRNA-derived small RNAs, we sequence long and short RNAs, small RNAs that co-precipitate with the Argonaute 2 protein and RNA fragments obtained in photoreactive nucleotide-enhanced crosslinking and immunoprecipitation (PAR-CLIP) of core snoRNA-associated proteins.

RESULTS

Analysis of these data sets reveals that many loci in the human genome reproducibly give rise to C/D box-like snoRNAs, whose expression and evolutionary conservation are typically less pronounced relative to the snoRNAs that are currently cataloged. We further find that virtually all C/D box snoRNAs are specifically processed inside the regions of terminal complementarity, retaining in the mature form only 4-5 nucleotides upstream of the C box and 2-5 nucleotides downstream of the D box. Sequencing of the total and Argonaute 2-associated populations of small RNAs reveals that despite their cellular abundance, C/D box-derived small RNAs are not efficiently incorporated into the Ago2 protein.

CONCLUSIONS

We conclude that the human genome encodes a large number of snoRNAs that are processed along the canonical pathway and expressed at relatively low levels. Generation of snoRNA-derived processing products with alternative, particularly miRNA-like, functions appears to be uncommon.

Bridging epigenomics and complex disease: the basics

The DNA sequence largely defines gene expression and phenotype. However, it is becoming increasingly clear that an additional chromatin-based regulatory network imparts both stability and plasticity to genome output, modifying phenotype independently of the genetic blueprint. Indeed, alterations in this "epigenetic" control layer underlie, at least in part, the reason for monozygotic twins being discordant for disease. Functionally, this regulatory layer comprises post-translational modifications of DNA and histones, as well as small and large noncoding RNAs. Together these regulate gene expression by changing chromatin organization and DNA accessibility. Successive technological advances over the past decade have enabled researchers to map the chromatin state with increasing accuracy and comprehensiveness, catapulting genetic research into a genome-wide era. Here, aiming particularly at the genomics/epigenomics newcomer, we review the epigenetic basis that has helped drive the technological shift and how this progress is shaping our understanding of complex disease.

miRDeep*: an integrated application tool for miRNA identification from RNA sequencing data

miRDeep and its varieties are widely used to quantify known and novel micro RNA (miRNA) from small RNA sequencing (RNAseq). This article describes miRDeep*, our integrated miRNA identification tool, which is modeled off miRDeep, but the precision of detecting novel miRNAs is improved by introducing new strategies to identify precursor miRNAs. miRDeep* has a user-friendly graphic interface and accepts raw data in FastQ and Sequence Alignment Map (SAM) or the binary equivalent (BAM) format. Known and novel miRNA expression levels, as measured by the number of reads, are displayed in an interface, which shows each RNAseq read relative to the pre-miRNA hairpin. The secondary pre-miRNA structure and read locations for each predicted miRNA are shown and kept in a separate figure file. Moreover, the target genes of known and novel miRNAs are predicted using the TargetScan algorithm, and the targets are ranked according to the confidence score. miRDeep* is an integrated standalone application where sequence alignment, pre-miRNA secondary structure calculation and graphical display are purely Java coded. This application tool can be executed using a normal personal computer with 1.5 GB of memory. Further, we show that miRDeep* outperformed existing miRNA prediction tools using our LNCaP and other small RNAseq datasets. miRDeep* is freely available online at http://www.australianprostatecentre.org/research/software/mirdeep-star.