A computational screen for methylation guide snoRNAs in yeast

Computational limitations and future needs to unravel the full potential of 2'-O-Methylation and C/D box snoRNAs

This review evaluates the current state of C/D snoRNA databases and prediction tools in relation to 2'-O-methylation (2'-O-Me). It highlights the limitations of existing resources in accurately annotating and predicting guide snoRNAs, particularly for newly identified 2"-O-Me sites. We emphasize the need for advanced computational approaches specifically tailored to 2"-O-Me to enable the discovery and functional analysis of snoRNAs. Given the growing importance of 2'-O-Me in areas such as cancer epitranscriptomics, ribosome biogenesis, and heterogeneity, existing tools remain inadequate. As 2'-O-Me gains recognition as a potential biomarker and therapeutic target, more sophisticated methods are urgently needed to improve snoRNA annotation and prediction, facilitating biomedical advancements.

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.

Mapping snoRNA-target RNA interactions in an RNA-binding protein-dependent manner with chimeric eCLIP

BACKGROUND

Small nucleolar RNAs (snoRNAs) are non-coding RNAs that function in ribosome and spliceosome biogenesis, primarily by guiding modifying enzymes to specific sites on ribosomal RNA (rRNA) and spliceosomal RNA (snRNA). However, many orphan snoRNAs remain uncharacterized, with unidentified or unvalidated targets, and studies on additional snoRNA-associated proteins are limited.

RESULTS

We adapted an enhanced chimeric eCLIP approach to comprehensively profile snoRNA-target RNA interactions using both core and accessory snoRNA-binding proteins as baits. Using core snoRNA-binding proteins, we confirmed most annotated snoRNA-rRNA and snoRNA-snRNA interactions in mouse and human cell lines and called novel, high-confidence interactions for orphan snoRNAs. While some of these interactions result in chemical modification, others may have modification-independent functions. We showed that snoRNA ribonucleoprotein complexes containing certain accessory proteins, like WDR43 and NOLC1, enriched for specific subsets of snoRNA-target RNA interactions with distinct roles in ribosome and spliceosome biogenesis. Notably, we discovered that SNORD89 guides 2'-O-methylation at two neighboring sites in U2 snRNA that fine-tune splice site recognition.

CONCLUSIONS

Chimeric eCLIP of snoRNA-associating proteins enables a comprehensive framework for studying snoRNA-target interactions in an RNA-binding protein-dependent manner, revealing novel interactions and regulatory roles in RNA biogenesis.

IMA GENOME - F20 A draft genome assembly of Agroatheliarolfsii, Ceratobasidiumpapillatum, Pyrenopezizabrassicae, Neopestalotiopsismacadamiae, Sphaerellopsisfilum and genomic resources for Colletotrichumspaethianum and Colletotrichumfructicola

This is a genome announacment there is no abstract

SNORD113-114 cluster maintains haematopoietic stem cell self-renewal via orchestrating the translation machinery

Haematopoietic stem cells (HSCs) self-renew and differentiate to replenish the pool of blood cells, which require a low but finely tuned protein synthesis rate. Nonetheless, the translatome landscape in HSCs and how the translation machinery orchestrates HSC self-renewal remain largely elusive. Here we perform ultra-low-input Ribo-seq in HSCs, progenitor and lineage cells, and reveal HSC-specific translated genes involved in rRNA processing. We systematically profile small nucleolar RNAs (snoRNAs) and uncover an indispensable role of the SNORD113-114 cluster in regulating HSC self-renewal. Maternal knockout (Mat-KO) of this cluster substantially impairs HSC self-renewal, whereas loss of the paternal allele shows no obvious phenotype. Mechanistically, Mat-KO results in dysregulation of translation machinery (rRNA 2'-O-Me modifications, pre-rRNA processing, 60S ribosome assembly and translation) and induces nucleolar stress in HSCs, which exempts p53 from Mdm2-mediated proteasomal degradation and leads to apoptosis. Collectively, our study provides a promising facet to our understanding of snoRNA-mediated regulation in HSC homeostasis.

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.

Mapping snoRNA-target RNA interactions in an RNA binding protein-dependent manner with chimeric eCLIP

Small nucleolar RNAs (snoRNAs) are non-coding RNAs that function in ribosome and spliceosome biogenesis, primarily by guiding modifying enzymes to specific sites on ribosomal RNA (rRNA) and spliceosomal RNA (snRNA). However, many orphan snoRNAs remain uncharacterized, with unidentified or unvalidated targets, and studies on additional snoRNA-associated proteins are limited. We adapted an enhanced chimeric eCLIP approach to comprehensively profile snoRNA-target RNA interactions using both core and accessory snoRNA binding proteins as baits. Using core snoRNA binding proteins, we confirmed most annotated snoRNA-rRNA and snoRNA-snRNA interactions in mouse and human cell lines and called novel, high-confidence interactions for orphan snoRNAs. While some of these interactions result in chemical modification, others may have modification-independent functions. We then showed that snoRNA ribonucleoprotein complexes containing certain accessory proteins, like WDR43 and NOLC1, enriched for specific subsets of snoRNA-target RNA interactions with distinct roles in ribosome and spliceosome biogenesis. Notably, we discovered that SNORD89 guides 2'-O-methylation at two neighboring sites in U2 snRNA that are important for activating splicing, but also appear to ensure imperfect splicing for a subset of near-constitutive exons. Thus, chimeric eCLIP of snoRNA-associating proteins enables a comprehensive framework for studying snoRNA-target interactions in an RNA binding protein-dependent manner, revealing novel interactions and regulatory roles in RNA biogenesis.

Chromosome-level genome assembly and annotation of the black sea urchin Arbacia lixula (Linnaeus, 1758)

The black sea urchin (Arbacia lixula) is a keystone species inhabiting the coastal shallow waters of the Mediterranean Sea, which is a key driver of littoral communities' structure. Here, we present the first genome assembly and annotation of this species, standing as the first Arbacioida genome, including both nuclear and mitochondrial genomes. To obtain a chromosome-level assembly, we used a combination of PacBio high fidelity (HiFi) reads and chromatin capture reads (Omni-C). In addition, we generated a high-quality nuclear annotation of both coding and non-coding genes, by using published RNA-Seq data from several individuals of A. lixula and gene models from closely related species. The nuclear genome assembly has a total span of 607.91 Mb, being consistent with its experimentally estimated genome size. The assembly contains 22 chromosome-scale scaffolds (96.52% of the total length), which coincides with its known karyotype. A total of 72,767 transcripts were predicted from the nuclear genome, 24,171 coding, and 48,596 non-coding that included lncRNA, snoRNA, and tRNAs. The circularized mitochondrial genome had 15,740 bp comprising 13 protein-coding genes, 2 rRNA, and 22 tRNA. This reference genome will enhance ongoing A. lixula studies and benefit the wider sea urchin scientific community.

Growth-regulated co-occupancy of Mediator and Lsm3 at intronic ribosomal protein genes

Mediator is a well-known transcriptional co-regulator and serves as an adaptor between gene-specific regulatory proteins and RNA polymerase II. Studies on the chromatin-bound form of Mediator revealed interactions with additional protein complexes involved in various transcription-related processes, such as the Lsm2-8 complex that is part of the spliceosomal U6 small nuclear ribonucleoprotein complex. Here, we employ Chromatin Immunoprecipitation sequencing (ChIP-seq) of chromatin associated with the Lsm3 protein and the Med1 or Med15 Mediator subunits. We identify 86 genes co-occupied by both Lsm3 and Mediator, of which 73 were intron-containing ribosomal protein genes. In logarithmically growing cells, Mediator primarily binds to their promoter regions but also shows a second, less pronounced occupancy at their 3'-exons. During the late exponential phase, we observe a near-complete transition of Mediator from these promoters to a position in their 3'-ends, overlapping the Lsm3 binding sites ∼250 bp downstream of their last intron-exon boundaries. Using an unbiased RNA sequencing approach, we show that transition of Mediator from promoters to the last exon of these genes correlates to reduction of both their messenger RNA levels and splicing ratios, indicating that the Mediator and Lsm complexes cooperate to control growth-regulated expression of intron-containing ribosomal protein genes at the levels of transcription and splicing.

Loss of Conserved rRNA Modifications in the Peptidyl Transferase Center Leads to Diminished Protein Synthesis and Cell Growth in Budding Yeast

Ribosomal RNAs (rRNAs) are extensively modified during the transcription and subsequent maturation. Three types of modifications, 2'-O-methylation of ribose moiety, pseudouridylation, and base modifications, are introduced either by a snoRNA-driven mechanism or by stand-alone enzymes. Modified nucleotides are clustered at the functionally important sites, including peptidyl transferase center (PTC). Therefore, it has been hypothesised that the modified nucleotides play an important role in ensuring the functionality of the ribosome. In this study, we demonstrate that seven 25S rRNA modifications, including four evolutionarily conserved modifications, in the proximity of PTC can be simultaneously depleted without loss of cell viability. Yeast mutants lacking three snoRNA genes (snR34, snR52, and snR65) and/or expressing enzymatically inactive variants of spb1(D52A/E679K) and nop2(C424A/C478A) were constructed. The results show that rRNA modifications in PTC contribute collectively to efficient translation in eukaryotic cells. The deficiency of seven modified nucleotides in 25S rRNA resulted in reduced cell growth, cold sensitivity, decreased translation levels, and hyperaccurate translation, as indicated by the reduced missense and nonsense suppression. The modification m5C2870 is crucial in the absence of the other six modified nucleotides. Thus, the pattern of rRNA-modified nucleotides around the PTC is essential for optimal ribosomal translational activity and translational fidelity.

Comprehensive map of ribosomal 2'-O-methylation and C/D box snoRNAs in Drosophila melanogaster

During their maturation, ribosomal RNAs (rRNAs) are decorated by hundreds of chemical modifications that participate in proper folding of rRNA secondary structures and therefore in ribosomal function. Along with pseudouridine, methylation of the 2'-hydroxyl ribose moiety (Nm) is the most abundant modification of rRNAs. The majority of Nm modifications in eukaryotes are placed by Fibrillarin, a conserved methyltransferase belonging to a ribonucleoprotein complex guided by C/D box small nucleolar RNAs (C/D box snoRNAs). These modifications impact interactions between rRNAs, tRNAs and mRNAs, and some are known to fine tune translation rates and efficiency. In this study, we built the first comprehensive map of Nm sites in Drosophila melanogaster rRNAs using two complementary approaches (RiboMethSeq and Nanopore direct RNA sequencing) and identified their corresponding C/D box snoRNAs by whole-transcriptome sequencing. We de novo identified 61 Nm sites, from which 55 are supported by both sequencing methods, we validated the expression of 106 C/D box snoRNAs and we predicted new or alternative rRNA Nm targets for 31 of them. Comparison of methylation level upon different stresses show only slight but specific variations, indicating that this modification is relatively stable in D. melanogaster. This study paves the way to investigate the impact of snoRNA-mediated 2'-O-methylation on translation and proteostasis in a whole organism.

Complicated target recognition by archaeal box C/D guide RNAs

Box C/D RNAs guide the site-specific formation of 2'-O-methylated nucleotides (Nm) of RNAs in eukaryotes and archaea. Although C/D RNAs have been profiled in several archaea, their targets have not been experimentally determined. Here, we mapped Nm in rRNAs, tRNAs, and abundant small RNAs (sRNAs) and profiled C/D RNAs in the crenarchaeon Sulfolobus islandicus. The targets of C/D RNAs were assigned by analysis of base-pairing interactions, in vitro modification assays, and gene deletion experiments, revealing a complicated landscape of C/D RNA-target interactions. C/D RNAs widely use dual antisense elements to target adjacent sites in rRNAs, enhancing modification at weakly bound sites. Two consecutive sites can be guided with the same antisense element upstream of box D or D', a phenomenon known as double-specificity that is exclusive to internal box D' in eukaryotic C/D RNAs. Several C/D RNAs guide modification at a single non-canonical site. This study reveals the global landscape of RNA-guided 2'-O-methylation in an archaeon and unexpected targeting rules employed by C/D RNA.

LinearCoFold and LinearCoPartition: linear-time algorithms for secondary structure prediction of interacting RNA molecules

Many RNAs function through RNA-RNA interactions. Fast and reliable RNA structure prediction with consideration of RNA-RNA interaction is useful, however, existing tools are either too simplistic or too slow. To address this issue, we present LinearCoFold, which approximates the complete minimum free energy structure of two strands in linear time, and LinearCoPartition, which approximates the cofolding partition function and base pairing probabilities in linear time. LinearCoFold and LinearCoPartition are orders of magnitude faster than RNAcofold. For example, on a sequence pair with combined length of 26,190 nt, LinearCoFold is 86.8× faster than RNAcofold MFE mode, and LinearCoPartition is 642.3× faster than RNAcofold partition function mode. Surprisingly, LinearCoFold and LinearCoPartition's predictions have higher PPV and sensitivity of intermolecular base pairs. Furthermore, we apply LinearCoFold to predict the RNA-RNA interaction between SARS-CoV-2 genomic RNA (gRNA) and human U4 small nuclear RNA (snRNA), which has been experimentally studied, and observe that LinearCoFold's prediction correlates better with the wet lab results than RNAcofold's.

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.

The fission yeast methyl phosphate capping enzyme Bmc1 guides 2'-O-methylation of the U6 snRNA

Splicing requires the tight coordination of dynamic spliceosomal RNAs and proteins. U6 is the only spliceosomal RNA transcribed by RNA Polymerase III and undergoes an extensive maturation process. In humans and fission yeast, this includes addition of a 5' γ-monomethyl phosphate cap by members of the Bin3/MePCE family as well as snoRNA guided 2'-O-methylation. Previously, we have shown that the Bin3/MePCE homolog Bmc1 is recruited to the S. pombe telomerase holoenzyme by the LARP7 family protein Pof8, where it acts in a catalytic-independent manner to protect the telomerase RNA and facilitate holoenzyme assembly. Here, we show that Bmc1 and Pof8 are required for the formation of a distinct U6 snRNP that promotes 2'-O-methylation of U6, and identify a non-canonical snoRNA that guides this methylation. We also show that the 5' γ-monomethyl phosphate capping activity of Bmc1 is not required for its role in promoting snoRNA guided 2'-O-methylation, and that this role relies on different regions of Pof8 from those required for Pof8 function in telomerase. Our results are consistent with a novel role for Bmc1/MePCE family members in stimulating 2'-O-methylation and a more general role for Bmc1 and Pof8 in guiding noncoding RNP assembly beyond the telomerase RNP.

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.

Computational Methods and Deep Learning for Elucidating Protein Interaction Networks

Protein interactions play a critical role in all biological processes, but experimental identification of protein interactions is a time- and resource-intensive process. The advances in next-generation sequencing and multi-omics technologies have greatly benefited large-scale predictions of protein interactions using machine learning methods. A wide range of tools have been developed to predict protein-protein, protein-nucleic acid, and protein-drug interactions. Here, we discuss the applications, methods, and challenges faced when employing the various prediction methods. We also briefly describe ways to overcome the challenges and prospective future developments in the field of protein interaction biology.

iSnoDi-LSGT: identifying snoRNA-disease associations based on local similarity constraints and global topological constraints

Growing evidence proves that small nucleolar RNAs (snoRNAs) have important functions in various biological processes, the malfunction of which leads to the emergence and development of complex diseases. However, identifying snoRNA-disease associations is an ongoing challenging task due to the considerable time- and money-consuming biological experiments. Therefore, it is urgent to design efficient and economical methods for the identification of snoRNA-disease associations. In this regard, we propose a computational method named iSnoDi-LSGT, which utilizes snoRNA sequence similarity and disease similarity as local similarity constraints. The iSnoDi-LSGT predictor further employs network embedding technology to extract topological features of snoRNAs and diseases, based on which snoRNA topological similarity and disease topological similarity are calculated as global topological constraints. To the best of our knowledge, the iSnoDi-LSGT is the first computational method for snoRNA-disease association identification. The experimental results indicate that the iSnoDi-LSGT predictor can effectively predict unknown snoRNA-disease associations. The web server of the iSnoDi-LSGT predictor is freely available at http://bliulab.net/iSnoDi-LSGT.

The small nucleolar RNA SnoR28 regulates plant growth and development by directing rRNA maturation

Small nucleolar RNAs (snoRNAs) are noncoding RNAs (ncRNAs) that guide chemical modifications of structural RNAs, which are essential for ribosome assembly and function in eukaryotes. Although numerous snoRNAs have been identified in plants by high-throughput sequencing, the biological functions of most of these snoRNAs remain unclear. Here, we identified box C/D SnoR28.1s as important regulators of plant growth and development by screening a CRISPR/Cas9-generated ncRNA deletion mutant library in Arabidopsis thaliana. Deletion of the SnoR28.1 locus, which contains a cluster of three genes producing SnoR28.1s, resulted in defects in root and shoot growth. SnoR28.1s guide 2'-O-ribose methylation of 25S rRNA at G2396. SnoR28.1s facilitate proper and efficient pre-rRNA processing, as the SnoR28.1 deletion mutants also showed impaired ribosome assembly and function, which may account for the growth defects. SnoR28 contains a 7-bp antisense box, which is required for 2'-O-ribose methylation of 25S rRNA at G2396, and an 8-bp extra box that is complementary to a nearby rRNA methylation site and is partially responsible for methylation of G2396. Both of these motifs are required for proper and efficient pre-rRNA processing. Finally, we show that SnoR28.1s genetically interact with HIDDEN TREASURE2 and NUCLEOLIN1. Our results advance our understanding of the roles of snoRNAs in Arabidopsis.

The snoGloBe interaction predictor reveals a broad spectrum of C/D snoRNA RNA targets

Box C/D small nucleolar RNAs (snoRNAs) are a conserved class of RNA known for their role in guiding ribosomal RNA 2'-O-ribose methylation. Recently, C/D snoRNAs were also implicated in regulating the expression of non-ribosomal genes through different modes of binding. Large scale RNA-RNA interaction datasets detect many snoRNAs binding messenger RNA, but are limited by specific experimental conditions. To enable a more comprehensive study of C/D snoRNA interactions, we created snoGloBe, a human C/D snoRNA interaction predictor based on a gradient boosting classifier. SnoGloBe considers the target type, position and sequence of the interactions, enabling it to outperform existing predictors. Interestingly, for specific snoRNAs, snoGloBe identifies strong enrichment of interactions near gene expression regulatory elements including splice sites. Abundance and splicing of predicted targets were altered upon the knockdown of their associated snoRNA. Strikingly, the predicted snoRNA interactions often overlap with the binding sites of functionally related RNA binding proteins, reinforcing their role in gene expression regulation. SnoGloBe is also an excellent tool for discovering viral RNA targets, as shown by its capacity to identify snoRNAs targeting the heavily methylated SARS-CoV-2 RNA. Overall, snoGloBe is capable of identifying experimentally validated binding sites and predicting novel sites with shared regulatory function.

Small and intermediate size structural RNAs in the unicellular parasite Cryptosporidium parvum as revealed by sRNA-seq and comparative genomics

Small and intermediate-size noncoding RNAs (sRNAs and 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 Cryptosporidium parvum, an obligate zoonotic protist parasite responsible for the diarrhoeal disease cryptosporidiosis. Using Illumina sequencing of a small RNA library, a systematic identification of novel small and is-ncRNAs was performed in C. parvum excysted sporozoites. A total of 79 novel is-ncRNA candidates, including antisense, intergenic and intronic is-ncRNAs, were identified, including 7 new small nucleolar RNAs (snoRNAs). Expression of select novel is-ncRNAs was confirmed by RT-PCR. Phylogenetic conservation was analysed using covariance models (CMs) in related Cryptosporidium and apicomplexan parasite genome sequences. A potential new type of small ncRNA derived from tRNA fragments was observed. Overall, a deep profiling analysis of novel is-ncRNAs in C. parvum and related species revealed structural features and conservation of these novel is-ncRNAs. Covariance models can be used to detect is-ncRNA genes in other closely related parasites. These findings provide important new sequences for additional functional characterization of novel is-ncRNAs in the protist pathogen C. parvum.

The Chromosome-Scale Reference Genome of Macadamia tetraphylla Provides Insights Into Fatty Acid Biosynthesis

Macadamia is an evergreen tree belonging to the Proteaceae family. The two commercial macadamia species, Macadamia integrifolia and M. tetraphylla, are highly prized for their edible kernels. The M. integrifolia genome was recently sequenced, but the genome of M. tetraphylla has to date not been published, which limits the study of biological research and breeding in this species. This study reports a high-quality genome sequence of M. tetraphylla based on the Oxford Nanopore Technologies technology and high-throughput chromosome conformation capture techniques (Hi-C). An assembly of 750.87 Mb with 51.11 Mb N50 length was generated, close to the 740 and 758 Mb size estimates by flow cytometry and k-mer analysis, respectively. Genome annotation indicated that 61.42% of the genome is composed of repetitive sequences and 34.95% is composed of long terminal repeat retrotransposons. Up to 31,571 protein-coding genes were predicted, of which 92.59% were functionally annotated. The average gene length was 6,055 bp. Comparative genome analysis revealed that the gene families associated with defense response, lipid transport, steroid biosynthesis, triglyceride lipase activity, and fatty acid metabolism are expanded in the M. tetraphylla genome. The distribution of fourfold synonymous third-codon transversion showed a recent whole-genome duplication event in M. tetraphylla. Genomic and transcriptomic analysis identified 187 genes encoding 33 crucial oil biosynthesis enzymes, depicting a comprehensive map of macadamia lipid biosynthesis. Besides, the 55 identified WRKY genes exhibited preferential expression in root as compared to that in other tissues. The genome sequence of M. tetraphylla provides novel insights for breeding novel varieties and genetic improvement of agronomic traits.

The RNA-binding protein RBP45D of Arabidopsis promotes transgene silencing and flowering time

RNA-directed DNA methylation (RdDM) helps to defend plants against invasive nucleic acids. In the canonical form of RdDM, 24-nt small interfering RNAs (siRNAs) are produced by DICER-LIKE 3 (DCL3). The siRNAs are loaded onto ARGONAUTE (AGO) proteins leading ultimately to de novo DNA methylation. Here, we introduce the Arabidopsis thaliana prors1 (LUC) transgenic system, in which 24-nt siRNAs are generated to silence the promoter-LUC construct. A forward genetic screen performed with this system identified, besides known components of RdDM (NRPD2A, RDR2, AGO4 and AGO6), the RNA-binding protein RBP45D. RBP45D is involved in CHH (where H is A, C or T) DNA methylation, and maintains siRNA production originating from the LUC transgene. RBP45D is localized to the nucleus, where it is associated with small nuclear RNAs (snRNAs) and small nucleolar RNAs (snoRNAs). RNA-Seq analysis showed that in CRISPR/Cas-mediated rbp-ko lines FLOWERING LOCUS C (FLC) mRNA levels are upregulated and several loci differentially spliced, among them FLM. In consequence, loss of RBP45D delays flowering, presumably mediated by the release of FLC levels and/or alternative splicing of FLM. Moreover, because levels and processing of transcripts of known RdDM genes are not altered in rbp-ko lines, RBP45D should have a more direct function in transgene silencing, probably independent of the canonical RdDM pathway. We suggest that RBP45D facilitates siRNA production by stabilizing either the precursor RNA or the slicer protein. Alternatively, RBP45D could be involved in chromatin modifications, participate in retention of Pol IV transcripts and/or in Pol V-dependent lncRNA retention in chromatin to enable their scaffold function.

Fractional 2'-O-methylation in the ribosomal RNA of Dictyostelium discoideum supports ribosome heterogeneity in Amoebozoa

A hallmark of ribosomal RNA (rRNA) are 2′-O-methyl groups that are introduced sequence specifically by box C/D small nucleolar RNAs (snoRNAs) in ribonucleoprotein particles. Most data on this chemical modification and its impact on RNA folding and stability are derived from organisms of the Opisthokonta supergroup. Using bioinformatics and RNA-seq data, we identify 30 novel box C/D snoRNAs in Dictyostelium discoideum, many of which are differentially expressed during the multicellular development of the amoeba. By applying RiboMeth-seq, we find 49 positions in the 17S and 26S rRNA 2′-O-methylated. Several of these nucleotides are substoichiometrically modified, with one displaying dynamic modification levels during development. Using homology-based models for the D. discoideum rRNA secondary structures, we localize many modified nucleotides in the vicinity of the ribosomal A, P and E sites. For most modified positions, a guiding box C/D snoRNA could be identified, allowing to determine idiosyncratic features of the snoRNA/rRNA interactions in the amoeba. Our data from D. discoideum represents the first evidence for ribosome heterogeneity in the Amoebozoa supergroup, allowing to suggest that it is a common feature of all eukaryotes.

Systematic mapping of rRNA 2'-O methylation during frog development and involvement of the methyltransferase Fibrillarin in eye and craniofacial development in Xenopus laevis

Ribosomes are essential nanomachines responsible for protein production. Although ribosomes are present in every living cell, ribosome biogenesis dysfunction diseases, called ribosomopathies, impact particular tissues specifically. Here, we evaluate the importance of the box C/D snoRNA-associated ribosomal RNA methyltransferase fibrillarin (Fbl) in the early embryonic development of Xenopus laevis. We report that in developing embryos, the neural plate, neural crest cells (NCCs), and NCC derivatives are rich in fbl transcripts. Fbl knockdown leads to striking morphological defects affecting the eyes and craniofacial skeleton, due to lack of NCC survival caused by massive p53-dependent apoptosis. Fbl is required for efficient pre-rRNA processing and 18S rRNA production, which explains the early developmental defects. Using RiboMethSeq, we systematically reinvestigated ribosomal RNA 2’-O methylation in X. laevis, confirming all 89 previously mapped sites and identifying 15 novel putative positions in 18S and 28S rRNA. Twenty-three positions, including 10 of the new ones, were validated orthogonally by low dNTP primer extension. Bioinformatic screening of the X. laevis transcriptome revealed candidate box C/D snoRNAs for all methylated positions. Mapping of 2’-O methylation at six developmental stages in individual embryos indicated a trend towards reduced methylation at specific positions during development. We conclude that fibrillarin knockdown in early Xenopus embryos causes reduced production of functional ribosomal subunits, thus impairing NCC formation and migration.

Ribosomes are essential nanomachines responsible for protein production in all cells. Ribosomopathies are diseases caused by improper ribosome formation due to mutations in ribosomal proteins or ribosome assembly factors. Such diseases primarily affect the brain and blood, and it is unclear how malfunctioning of a process as general as ribosome formation can lead to tissue-specific diseases. Here we have examined how fibrillarin, an enzyme which modifies ribosomal RNA by adding methyl groups at specific sites, affects early embryonic development in the frog Xenopus laevis. We have revealed its importance in the maturation of cells forming an embryonic structure called the neural crest. Fibrillarin depletion leads to reduced eye size and abnormal head shape, reminiscent of other conditions such as Treacher Collins syndrome. Molecularly, the observed phenotypes are explainable by increased p53-dependent programmed cell death triggered by inhibition of certain pre-rRNA processing steps. Our systematic investigation of the ribosomal RNA 2’-O methylation repertoire across development has further revealed hypomodification at a late stage of development, which might play a role in late developmental transitions involving differential translation by compositionally different ribosomes.

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.

Finding differentially expressed sRNA-Seq regions with srnadiff

Small RNAs (sRNAs) encompass a great variety of molecules of different kinds, such as microRNAs, small interfering RNAs, Piwi-associated RNA, among others. These sRNAs have a wide range of activities, which include gene regulation, protection against virus, transposable element silencing, and have been identified as a key actor in determining the development of the cell. Small RNA sequencing is thus routinely used to assess the expression of the diversity of sRNAs, usually in the context of differentially expression, where two conditions are compared. Tools that detect differentially expressed microRNAs are numerous, because microRNAs are well documented, and the associated genes are well defined. However, tools are lacking to detect other types of sRNAs, which are less studied, and whose precursor RNA is not well characterized. We present here a new method, called srnadiff, which finds all kinds of differentially expressed sRNAs. To the extent of our knowledge, srnadiff is the first tool that detects differentially expressed sRNAs without the use of external information, such as genomic annotation or additional sequences of sRNAs.

Translational control by DHX36 binding to 5'UTR G-quadruplex is essential for muscle stem-cell regenerative functions

Skeletal muscle has a remarkable ability to regenerate owing to its resident stem cells (also called satellite cells, SCs). SCs are normally quiescent; when stimulated by damage, they activate and expand to form new fibers. The mechanisms underlying SC proliferative progression remain poorly understood. Here we show that DHX36, a helicase that unwinds RNA G-quadruplex (rG4) structures, is essential for muscle regeneration by regulating SC expansion. DHX36 (initially named RHAU) is barely expressed at quiescence but is highly induced during SC activation and proliferation. Inducible deletion of Dhx36 in adult SCs causes defective proliferation and muscle regeneration after damage. System-wide mapping in proliferating SCs reveals DHX36 binding predominantly to rG4 structures at various regions of mRNAs, while integrated polysome profiling shows that DHX36 promotes mRNA translation via 5′-untranslated region (UTR) rG4 binding. Furthermore, we demonstrate that DHX36 specifically regulates the translation of Gnai2 mRNA by unwinding its 5′ UTR rG4 structures and identify GNAI2 as a downstream effector of DHX36 for SC expansion. Altogether, our findings uncover DHX36 as an indispensable post-transcriptional regulator of SC function and muscle regeneration acting through binding and unwinding rG4 structures at 5′ UTR of target mRNAs.

Skeletal muscle stem cells (or satellite cells, SCs) are normally quiescent but activate and expand in response to injury. Here the authors show that induction of DHX36 helicase during SC activation promotes mRNA translation by binding to 5′UTR mRNA G-quadruplexes (rG4) in targets including Gnai2 and unwinding them.

Saccharomyces cerevisiae, a Powerful Model for Studying rRNA Modifications and Their Effects on Translation Fidelity

Ribosomal RNA is a major component of the ribosome. This RNA plays a crucial role in ribosome functioning by ensuring the formation of the peptide bond between amino acids and the accurate decoding of the genetic code. The rRNA carries many chemical modifications that participate in its maturation, the formation of the ribosome and its functioning. In this review, we present the different modifications and how they are deposited on the rRNA. We also describe the most recent results showing that the modified positions are not 100% modified, which creates a heterogeneous population of ribosomes. This gave rise to the concept of specialized ribosomes that we discuss. The knowledge accumulated in the yeast Saccharomyces cerevisiae is very helpful to better understand the role of rRNA modifications in humans, especially in ribosomopathies.

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.

Profiling of RNA ribose methylation in Arabidopsis thaliana

Eukaryotic rRNAs and snRNAs are decorated with abundant 2′-O-methylated nucleotides (Nm) that are predominantly synthesized by box C/D snoRNA-guided enzymes. In the model plant Arabidopsis thaliana, C/D snoRNAs have been well categorized, but there is a lack of systematic mapping of Nm. Here, we applied RiboMeth-seq to profile Nm in cytoplasmic, chloroplast and mitochondrial rRNAs and snRNAs. We identified 111 Nm in cytoplasmic rRNAs and 19 Nm in snRNAs and assigned guide for majority of the detected sites using an updated snoRNA list. At least four sites are directed by guides with multiple specificities as shown in yeast. We found that C/D snoRNAs frequently form extra pairs with nearby sequences of methylation sites, potentially facilitating the substrate binding. Chloroplast and mitochondrial rRNAs contain five almost identical methylation sites, including two novel sites mediating ribosomal subunit joining. Deletion of FIB1 or FIB2 gene reduced the accumulation of C/D snoRNA and rRNA methylation with FIB1 playing a bigger role in methylation. Our data reveal the comprehensive 2′-O-methylation maps for Arabidopsis rRNAs and snRNAs and would facilitate study of their function and biosynthesis.

Analysis of fungal genomes reveals commonalities of intron gain or loss and functions in intron-poor species

Previous evolutionary reconstructions have concluded that early eukaryotic ancestors including both the last common ancestor of eukaryotes and of all fungi had intron-rich genomes. By contrast, some extant eukaryotes have few introns, underscoring the complex histories of intron–exon structures, and raising the question as to why these few introns are retained. Here, we have used recently available fungal genomes to address a variety of questions related to intron evolution. Evolutionary reconstruction of intron presence and absence using 263 diverse fungal species supports the idea that massive intron reduction through intron loss has occurred in multiple clades. The intron densities estimated in various fungal ancestors differ from zero to 7.6 introns per 1 kb of protein-coding sequence. Massive intron loss has occurred not only in microsporidian parasites and saccharomycetous yeasts, but also in diverse smuts and allies. To investigate the roles of the remaining introns in highly-reduced species, we have searched for their special characteristics in eight intron-poor fungi. Notably, the introns of ribosome-associated genes RPL7 and NOG2 have conserved positions; both intron-containing genes encoding snoRNAs. Furthermore, both the proteins and snoRNAs are involved in ribosome biogenesis, suggesting that the expression of the protein-coding genes and noncoding snoRNAs may be functionally coordinated. Indeed, these introns are also conserved in three-quarters of fungi species. Our study shows that fungal introns have a complex evolutionary history and underappreciated roles in gene expression.

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.

Draft genome of the herbaceous bamboo Raddia distichophylla

Bamboos are important nontimber forest plants widely distributed in the tropical and subtropical regions of Asia, Africa, America, and Pacific islands. They comprise the Bambusoideae in the grass family (Poaceae), including approximately 1700 described species in 127 genera. In spite of the widespread uses of bamboo for food, construction, and bioenergy, the gene repertoire of bamboo still remains largely unexplored. Raddia distichophylla (Schrad. ex Nees) Chase, belonging to the tribe Olyreae (Bambusoideae, Poaceae), is a diploid herbaceous bamboo with only slightly lignified stems. In this study, we report a draft genome assembly of the ∼589 Mb whole-genome sequence of R. distichophylla with a contig N50 length of 86.36 Kb. Repeat sequences account for ∼49.08% of the genome assembly, of which LTR retrotransposons occupy ∼35.99% of the whole genome. A total of 30,763 protein-coding genes were annotated in the R. distichophylla genome with an average transcript size of 2887 bp. Access to this herbaceous bamboo genome sequence will provide novel insights into biochemistry, molecular marker-assisted breeding programs, and germplasm conservation for bamboo species worldwide.

RMDisease: a database of genetic variants that affect RNA modifications, with implications for epitranscriptome pathogenesis

Deciphering the biological impacts of millions of single nucleotide variants remains a major challenge. Recent studies suggest that RNA modifications play versatile roles in essential biological mechanisms, and are closely related to the progression of various diseases including multiple cancers. To comprehensively unveil the association between disease-associated variants and their epitranscriptome disturbance, we built RMDisease, a database of genetic variants that can affect RNA modifications. By integrating the prediction results of 18 different RNA modification prediction tools and also 303,426 experimentally-validated RNA modification sites, RMDisease identified a total of 202,307 human SNPs that may affect (add or remove) sites of eight types of RNA modifications (m6A, m5C, m1A, m5U, Ψ, m6Am, m7G and Nm). These include 4,289 disease-associated variants that may imply disease pathogenesis functioning at the epitranscriptome layer. These SNPs were further annotated with essential information such as post-transcriptional regulations (sites for miRNA binding, interaction with RNA-binding proteins and alternative splicing) revealing putative regulatory circuits. A convenient graphical user interface was constructed to support the query, exploration and download of the relevant information. RMDisease should make a useful resource for studying the epitranscriptome impact of genetic variants via multiple RNA modifications with emphasis on their potential disease relevance. RMDisease is freely accessible at: www.xjtlu.edu.cn/biologicalsciences/rmd.

Epitranscriptomics of Mammalian Mitochondrial Ribosomal RNA

Modified nucleotides are present in all ribosomal RNA molecules. Mitochondrial ribosomes are unique to have a set of methylated residues that includes universally conserved ones, those that could be found either in bacterial or in archaeal/eukaryotic cytosolic ribosomes and those that are present exclusively in mitochondria. A single pseudouridine within the mt-rRNA is located in the peptidyltransferase center at a position similar to that in bacteria. After recent completion of the list of enzymes responsible for the modification of mammalian mitochondrial rRNA it became possible to summarize an evolutionary history, functional role of mt-rRNA modification enzymes and an interplay of the mt-rRNA modification and mitoribosome assembly process, which is a goal of this review.

The Existence and Localization of Nuclear snoRNAs in Arabidopsis thaliana Revisited

Ribosome biogenesis is one cell function-defining process. It depends on efficient transcription of rDNAs in the nucleolus as well as on the cytosolic synthesis of ribosomal proteins. For newly transcribed rRNA modification and ribosomal protein assembly, so-called small nucleolar RNAs (snoRNAs) and ribosome biogenesis factors (RBFs) are required. For both, an inventory was established for model systems like yeast and humans. For plants, many assignments are based on predictions. Here, RNA deep sequencing after nuclei enrichment was combined with single molecule species detection by northern blot and in vivo fluorescence in situ hybridization (FISH)-based localization studies. In addition, the occurrence and abundance of selected snoRNAs in different tissues were determined. These approaches confirm the presence of most of the database-deposited snoRNAs in cell cultures, but some of them are localized in the cytosol rather than in the nucleus. Further, for the explored snoRNA examples, differences in their abundance in different tissues were observed, suggesting a tissue-specific function of some snoRNAs. Thus, based on prediction and experimental confirmation, many plant snoRNAs can be proposed, while it cannot be excluded that some of the proposed snoRNAs perform alternative functions than are involved in rRNA modification.

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.

Draft genomes of two outcrossing wild rice, Oryza rufipogon and O. longistaminata, reveal genomic features associated with mating-system evolution

Oryza rufipogon and O. longistaminata are important wild relatives of cultivated rice, harboring a promising source of novel genes for rice breeding programs. Here, we present de novo assembled draft genomes and annotation of O. rufipogon and O. longistaminata. Our analysis reveals a considerable number of lineage-specific gene families associated with the self-incompatibility (SI) and formation of reproductive separation. We show how lineage-specific expansion or contraction of gene families with functional enrichment of the recognition of pollen, thus enlightening their reproductive diversification. We documented a large number of lineage-specific gene families enriched in salt stress, antifungal response, and disease resistance. Our comparative analysis further shows a genome-wide expansion of genes encoding NBS-LRR proteins in these two outcrossing wild species in contrast to six other selfing rice species. Conserved noncoding sequences (CNSs) in the two wild rice genomes rapidly evolve relative to selfing rice species, resulting in the reduction of genomic variation owing to shifts of mating systems. We find that numerous genes related to these rapidly evolving CNSs are enriched in reproductive structure development, flower development, and postembryonic development, which may associate with SI in O. rufipogon and O. longistaminata.

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.

The draft genome sequence of an upland wild rice species, Oryza granulata

Exploiting novel gene sources from wild relatives has proven to be an efficient approach to advance crop genetic breeding efforts. Oryza granulata, with the GG genome type, occupies the basal position of the Oryza phylogeny and has the second largest genome (~882 Mb). As an upland wild rice species, it possesses renowned traits that distinguish it from other Oryza species, such as tolerance to shade and drought, immunity to bacterial blight and resistance to the brown planthopper. Here, we generated a 736.66-Mb genome assembly of O. granulata with 40,131 predicted protein-coding genes. With Hi-C data, for the first time, we anchored ~98.2% of the genome assembly to the twelve pseudo-chromosomes. This chromosome-length genome assembly of O. granulata will provide novel insights into rice genome evolution, enhance our efforts to search for new genes for future rice breeding programmes and facilitate the conservation of germplasm of this endangered wild rice species.

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.

Dual-initiation promoters with intertwined canonical and TCT/TOP transcription start sites diversify transcript processing

Variations in transcription start site (TSS) selection reflect diversity of preinitiation complexes and can impact on post-transcriptional RNA fates. Most metazoan polymerase II-transcribed genes carry canonical initiation with pyrimidine/purine (YR) dinucleotide, while translation machinery-associated genes carry polypyrimidine initiator (5’-TOP or TCT). By addressing the developmental regulation of TSS selection in zebrafish we uncovered a class of dual-initiation promoters in thousands of genes, including snoRNA host genes. 5’-TOP/TCT initiation is intertwined with canonical initiation and used divergently in hundreds of dual-initiation promoters during maternal to zygotic transition. Dual-initiation in snoRNA host genes selectively generates host and snoRNA with often different spatio-temporal expression. Dual-initiation promoters are pervasive in human and fruit fly, reflecting evolutionary conservation. We propose that dual-initiation on shared promoters represents a composite promoter architecture, which can function both coordinately and divergently to diversify RNAs.

The functional significance of start site choice in promoter architectures is little understood. Here the authors identify in zebrafish development and mammalian cells a class of dual-initiation promoters, in which non-canonical YC dinucleotides reflecting 5’-TOP/TCT initiation are intertwined with canonical YR-initiation.

A single H/ACA small nucleolar RNA mediates tumor suppression downstream of oncogenic RAS

Small nucleolar RNAs (snoRNAs) are a diverse group of non-coding RNAs that direct chemical modifications at specific residues on other RNA molecules, primarily on ribosomal RNA (rRNA). SnoRNAs are altered in several cancers; however, their role in cell homeostasis as well as in cellular transformation remains poorly explored. Here, we show that specific subsets of snoRNAs are differentially regulated during the earliest cellular response to oncogenic RAS^G12V expression. We describe a novel function for one H/ACA snoRNA, SNORA24, which guides two pseudouridine modifications within the small ribosomal subunit, in RAS-induced senescence in vivo. We find that in mouse models, loss of Snora24 cooperates with RAS^G12V to promote the development of liver cancer that closely resembles human steatohepatitic hepatocellular carcinoma (HCC). From a clinical perspective, we further show that human HCCs with low SNORA24 expression display increased lipid content and are associated with poor patient survival. We next asked whether ribosomes lacking SNORA24-guided pseudouridine modifications on 18S rRNA have alterations in their biophysical properties. Single-molecule Fluorescence Resonance Energy Transfer (FRET) analyses revealed that these ribosomes exhibit perturbations in aminoacyl-transfer RNA (aa-tRNA) selection and altered pre-translocation ribosome complex dynamics. Furthermore, we find that HCC cells lacking SNORA24-guided pseudouridine modifications have increased translational miscoding and stop codon readthrough frequencies. These findings highlight a role for specific snoRNAs in safeguarding against oncogenic insult and demonstrate a functional link between H/ACA snoRNAs regulated by RAS and the biophysical properties of ribosomes in cancer.

Methylation guide RNA evolution in archaea: structure, function and genomic organization of 110 C/D box sRNA families across six Pyrobaculum species

Archaeal homologs of eukaryotic C/D box small nucleolar RNAs (C/D box sRNAs) guide precise 2′-O-methyl modification of ribosomal and transfer RNAs. Although C/D box sRNA genes constitute one of the largest RNA gene families in archaeal thermophiles, most genomes have incomplete sRNA gene annotation because reliable, fully automated detection methods are not available. We expanded and curated a comprehensive gene set across six species of the crenarchaeal genus Pyrobaculum, particularly rich in C/D box sRNA genes. Using high-throughput small RNA sequencing, specialized computational searches and comparative genomics, we analyzed 526 Pyrobaculum C/D box sRNAs, organizing them into 110 families based on synteny and conservation of guide sequences which determine methylation targets. We examined gene duplications and rearrangements, including one family that has expanded in a pattern similar to retrotransposed repetitive elements in eukaryotes. New training data and inclusion of kink-turn secondary structural features enabled creation of an improved search model. Our analyses provide the most comprehensive, dynamic view of C/D box sRNA evolutionary history within a genus, in terms of modification function, feature plasticity, and gene mobility.

Landscape of the complete RNA chemical modifications in the human 80S ribosome

During ribosome biogenesis, ribosomal RNAs acquire various chemical modifications that ensure the fidelity of translation, and dysregulation of the modification processes can cause proteome changes as observed in cancer and inherited human disorders. Here, we report the complete chemical modifications of all RNAs of the human 80S ribosome as determined with quantitative mass spectrometry. We assigned 228 sites with 14 different post-transcriptional modifications, most of which are located in functional regions of the ribosome. All modifications detected are typical of eukaryotic ribosomal RNAs, and no human-specific modifications were observed, in contrast to a recently reported cryo-electron microscopy analysis. While human ribosomal RNAs appeared to have little polymorphism regarding the post-transcriptional modifications, we found that pseudouridylation at two specific sites in 28S ribosomal RNA are significantly reduced in ribosomes of patients with familial dyskeratosis congenita, a genetic disease caused by a point mutation in the pseudouridine synthase gene DKC1. The landscape of the entire epitranscriptomic ribosomal RNA modifications provides a firm basis for understanding ribosome function and dysfunction associated with human disease.

Structural and evolutionary insights into ribosomal RNA methylation

Methylation of nucleotides in ribosomal RNAs (rRNAs) is a ubiquitous feature that occurs in all living organisms. Identification of all enzymes responsible for rRNA methylation, as well as mapping of all modified rRNA residues, is now complete for a number of model species, such as Escherichia coli and Saccharomyces cerevisiae. Recent high-resolution structures of bacterial ribosomes provided the first direct visualization of methylated nucleotides. The structures of ribosomes from various organisms and organelles have also lately become available, enabling comparative structure-based analysis of rRNA methylation sites in various taxonomic groups. In addition to the conserved core of modified residues in ribosomes from the majority of studied organisms, structural analysis points to the functional roles of some of the rRNA methylations, which are discussed in this Review in an evolutionary context.