Motif conservation, stability, and host gene expression are the main drivers of snoRNA expression across vertebrates.

Small nucleolar RNAs (snoRNAs) are structured noncoding RNAs present in multiple copies within eukaryotic genomes. snoRNAs guide chemical modifications on their target RNA and regulate processes like ribosome assembly and splicing. Most human snoRNAs are embedded within host gene introns, the remainder being independently expressed from intergenic regions. We recently characterized the abundance of snoRNAs and their host gene across several healthy human tissues and found that the level of most snoRNAs does not correlate with that of their host gene, with the observation that snoRNAs embedded within the same host gene often differ drastically in abundance. To better understand the determinants of snoRNA expression, we trained machine learning models to predict whether snoRNAs are expressed or not in human tissues based on more than 30 collected features related to snoRNAs and their genomic context. By interpreting the models' predictions, we find that snoRNAs rely on conserved motifs, a stable global structure and terminal stem, and a transcribed locus to be expressed. We observe that these features explain well the varying abundance of snoRNAs embedded within the same host gene. By predicting the expression status of snoRNAs across several vertebrates, we notice that only one-third of all annotated snoRNAs are expressed per genome, as in humans. Our results suggest that ancestral snoRNAs disseminated within vertebrate genomes, sometimes leading to the development of new functions and a probable gain in fitness and thereby conserving features favorable to the expression of these few snoRNAs, the large remainder often degenerating into pseudogenes.

snoDB 2.0: an enhanced interactive database, specializing in human snoRNAs.

snoDB is an interactive database of human small nucleolar RNAs (snoRNAs) that includes up-to-date information on snoRNA features, genomic location, conservation, host gene, snoRNA-RNA targets and snoRNA abundance and provides links to other resources. In the second edition of this database (snoDB 2.0), we added an entirely new section on ribosomal RNA (rRNA) chemical modifications guided by snoRNAs with easy navigation between the different rRNA versions used in the literature and experimentally measured levels of modification. We also included new layers of information, including snoRNA motifs, secondary structure prediction, snoRNA-protein interactions, copy annotations and low structure bias expression data in a wide panel of tissues and cell lines to bolster functional probing of snoRNA biology. Version 2.0 features updated identifiers, more links to external resources and duplicate entry resolution. As a result, snoDB 2.0, which is freely available at https://bioinfo-scottgroup.med.usherbrooke.ca/snoDB/, represents a one-stop shop for snoRNA features, rRNA modification targets, functional impact and potential regulators.

The regulatory roles of small nucleolar RNAs within their host locus.

Small nucleolar RNAs (snoRNAs) are a class of conserved noncoding RNAs forming complexes with proteins to catalyse site-specific modifications on ribosomal RNA. Besides this canonical role, several snoRNAs are now known to regulate diverse levels of gene expression. While these functions are carried out in trans by mature snoRNAs, evidence has also been emerging of regulatory roles of snoRNAs in cis, either within their genomic locus or as longer transcription intermediates during their maturation. Herein, we review recent findings that snoRNAs can interact in cis with their intron to regulate the expression of their host gene. We also explore the ever-growing diversity of longer host-derived snoRNA extensions and their functional impact across the transcriptome. Finally, we discuss the role of snoRNA duplications into forging these new layers of snoRNA-mediated regulation, as well as their involvement in the genomic imprinting of their host locus.

snoDB: an interactive database of human snoRNA sequences, abundance and interactions.

Small nucleolar RNAs (snoRNAs) are an abundant type of non-coding RNA with conserved functions in all known eukaryotes. Classified into two main families, the box C/D and H/ACA snoRNAs, they enact their most well characterized role of guiding site specific modifications in ribosomal RNA, through the formation of specific ribonucleoprotein complexes, with fundamental implications in ribosome biogenesis. However, it is becoming increasingly clear that the landscape of snoRNA cellular functionality is much broader than it once seemed with novel members, non-uniform expression patterns, new and diverse targets as well as several emerging non-canonical functions ranging from the modulation of alternative splicing to the regulation of chromatin architecture. In order to facilitate the further characterization of human snoRNAs in a holistic manner, we introduce an online interactive database tool: snoDB. Its purpose is to consolidate information on human snoRNAs from different sources such as sequence databases, target information, both canonical and non-canonical from the literature and from high-throughput RNA-RNA interaction datasets, as well as high-throughput sequencing data that can be visualized interactively.

SnoRNA copy regulation affects family size, genomic location and family abundance levels.

BACKGROUND: Small nucleolar RNAs (snoRNAs) are an abundant class of noncoding RNAs present in all eukaryotes and best known for their involvement in ribosome biogenesis. In mammalian genomes, many snoRNAs exist in multiple copies, resulting from recombination and retrotransposition from an ancestral snoRNA. To gain insight into snoRNA copy regulation, we used Rfam classification and normal human tissue expression datasets generated using low structure bias RNA-seq to characterize snoRNA families. RESULTS: We found that although box H/ACA families are on average larger than box C/D families, the number of expressed members is similar for both types. Family members can cover a wide range of average abundance values, but importantly, expression variability of individual members of a family is preferred over the total variability of the family, especially for box H/ACA snoRNAs, suggesting that while members are likely differentially regulated, mechanisms exist to ensure uniformity of the total family abundance across tissues. Box C/D snoRNA family members are mostly embedded in the same host gene while box H/ACA family members tend to be encoded in more than one different host, supporting a model in which box C/D snoRNA duplication occurred mostly by cis recombination while box H/ACA snoRNA families have gained copy members through retrotransposition. And unexpectedly, snoRNAs encoded in the same host gene can be regulated independently, as some snoRNAs within the same family vary in abundance in a divergent way between tissues. CONCLUSIONS: SnoRNA copy regulation affects family sizes, genomic location of the members and controls simultaneously member and total family abundance to respond to the needs of individual tissues.

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.

Intronic small nucleolar RNAs regulate host gene splicing through base pairing with their adjacent intronic sequences.

BACKGROUND: Small nucleolar RNAs (snoRNAs) are abundant noncoding RNAs best known for their involvement in ribosomal RNA maturation. In mammals, most expressed snoRNAs are embedded in introns of longer genes and produced through transcription and splicing of their host. Intronic snoRNAs were long viewed as inert passengers with little effect on host expression. However, a recent study reported a snoRNA influencing the splicing and ultimate output of its host gene. Overall, the general contribution of intronic snoRNAs to host expression remains unclear. RESULTS: Computational analysis of large-scale human RNA-RNA interaction datasets indicates that 30% of detected snoRNAs interact with their host transcripts. Many snoRNA-host duplexes are located near alternatively spliced exons and display high sequence conservation suggesting a possible role in splicing regulation. The study of the model SNORD2-EIF4A2 duplex indicates that the snoRNA interaction with the host intronic sequence conceals the branch point leading to decreased inclusion of the adjacent alternative exon. Extended SNORD2 sequence containing the interacting intronic region accumulates in sequencing datasets in a cell-type-specific manner. Antisense oligonucleotides and mutations that disrupt the formation of the snoRNA-intron structure promote the splicing of the alternative exon, shifting the EIF4A2 transcript ratio away from nonsense-mediated decay. CONCLUSIONS: Many snoRNAs form RNA duplexes near alternative exons of their host transcripts, placing them in optimal positions to control host output as shown for the SNORD2-EIF4A2 model system. Overall, our study supports a more widespread role for intronic snoRNAs in the regulation of their host transcript maturation.

Altered tRNA processing is linked to a distinct and unusual La protein in Tetrahymena thermophila.

Nascent pre-tRNAs are transcribed by RNA polymerase III and immediately bound by La proteins on the UUU-3'OH sequence, using a tandem arrangement of the La motif and an adjacent RNA recognition motif-1 (RRM1), resulting in protection from 3'-exonucleases and promotion of pre-tRNA folding. The Tetrahymena thermophila protein Mlp1 has been previously classified as a genuine La protein, despite the predicted absence of the RRM1. We find that Mlp1 functions as a La protein through binding of pre-tRNAs, and affects pre-tRNA processing in Tetrahymena thermophila and when expressed in fission yeast. However, unlike in other examined eukaryotes, depletion of Mlp1 results in 3'-trailer stabilization. The 3'-trailers in Tetrahymena thermophila are uniquely short relative to other examined eukaryotes, and 5'-leaders have evolved to disfavour pre-tRNA leader/trailer pairing. Our data indicate that this variant Mlp1 architecture is linked to an altered, novel mechanism of tRNA processing in Tetrahymena thermophila.

Annotation of snoRNA abundance across human tissues reveals complex snoRNA-host gene relationships.

BACKGROUND: Small nucleolar RNAs (snoRNAs) are mid-size non-coding RNAs required for ribosomal RNA modification, implying a ubiquitous tissue distribution linked to ribosome synthesis. However, increasing numbers of studies identify extra-ribosomal roles of snoRNAs in modulating gene expression, suggesting more complex snoRNA abundance patterns. Therefore, there is a great need for mapping the snoRNome in different human tissues as the blueprint for snoRNA functions. RESULTS: We used a low structure bias RNA-Seq approach to accurately quantify snoRNAs and compare them to the entire transcriptome in seven healthy human tissues (breast, ovary, prostate, testis, skeletal muscle, liver, and brain). We identify 475 expressed snoRNAs categorized in two abundance classes that differ significantly in their function, conservation level, and correlation with their host gene: 390 snoRNAs are uniformly expressed and 85 are enriched in the brain or reproductive tissues. Most tissue-enriched snoRNAs are embedded in lncRNAs and display strong correlation of abundance with them, whereas uniformly expressed snoRNAs are mostly embedded in protein-coding host genes and are mainly non- or anticorrelated with them. Fifty-nine percent of the non-correlated or anticorrelated protein-coding host gene/snoRNA pairs feature dual-initiation promoters, compared to only 16% of the correlated non-coding host gene/snoRNA pairs. CONCLUSIONS: Our results demonstrate that snoRNAs are not a single homogeneous group of housekeeping genes but include highly regulated tissue-enriched RNAs. Indeed, our work indicates that the architecture of snoRNA host genes varies to uncouple the host and snoRNA expressions in order to meet the different snoRNA abundance levels and functional needs of human tissues.

Endothelial deletion of PKCδ prevents VEGF inhibition and restores blood flow reperfusion in diabetic ischemic limb.

AIMS: Peripheral artery disease is a complication of diabetes leading to critical hindlimb ischemia. Diabetes-induced inhibition of VEGF actions is associated with the activation of protein kinase Cδ (PKCδ). We aim to specifically investigate the role of PKCδ in endothelial cell (EC) function and VEGF signaling. METHODS: Nondiabetic and diabetic mice, with (ec-Prkcd-/-) or without (ec-Prkcdf/f) endothelial deletion of PKCδ, underwent femoral artery ligation. Blood flow reperfusion was assessed up to 4 weeks post-surgery. Capillary density, EC apoptosis and VEGF signaling were evaluated in the ischemic muscle. Src homology region 2 domain-containing phosphatase-1 (SHP-1) phosphatase activity was assessed in vitro using primary ECs. RESULTS: Ischemic muscle of diabetic ec-Prkcdf/f mice exhibited reduced blood flow reperfusion and capillary density while apoptosis increased as compared to nondiabetic ec-Prkcdf/f mice. In contrast, blood flow reperfusion and capillary density were significantly improved in diabetic ec-Prkcd-/- mice. VEGF signaling pathway was restored in diabetic ec-Prkcd-/- mice. The deletion of PKCδ in ECs prevented diabetes-induced VEGF unresponsiveness through a reduction of SHP-1 phosphatase activity. CONCLUSIONS: Our data provide new highlights in mechanisms by which PKCδ activation in EC contributed to poor collateral vessel formation, thus, offering novel therapeutic targets to improve angiogenesis in the diabetic limb.