The cellular landscape of mid-size noncoding RNA

Transcribing the enigma: the B chromosome as a territory of uncharted RNAs

B chromosomes are supernumerary elements found in several groups of eukaryotes, including fungi, plants, and animals. Typically, these chromosomes either originate from their hosts through errors in meiosis or interspecifically through horizontal transfer. While many B chromosomes are primarily heterochromatic and possess a low number of coding genes, these additional elements are still capable of transcribing sequences and exerting influence on the expression of host genes. How B chromosomes escape elimination and which impacts can be promoted in the cell always intrigued the cytogeneticists. In pursuit of understanding the behavior and functional impacts of these extra elements, cytogenetic studies meet the advances of molecular biology, incorporating various techniques into investigating B chromosomes from a functional perspective. In this review, we present a timeline of studies investigating B chromosomes and RNAs, highlighting the advances and key findings throughout their history. Additionally, we identified which RNA classes are reported in the B chromosomes and emphasized the necessity for further investigation into new perspectives on the B chromosome functions. In this context, we present a phylogenetic tree that illustrates which branches either report B chromosome presence or have functional RNA studies related to B chromosomes. We propose investigating other unexplored RNA classes and conducting functional analysis in conjunction with cytogenetic studies to enhance our understanding of the B chromosome from an RNA perspective.

7SK small nuclear RNA (Rn7SK) induces apoptosis through intrinsic and extrinsic pathways in human embryonic kidney cell line

BACKGROUND

Rn7SK, a highly conserved small nuclear non-coding RNA, controls Polymerase II transcription machinery by activating of the Positive Transcriptional Elongation Factor b (P-TEFb). Apart from its role in transcriptional regulation, the potential functions of Rn7SK in cell apoptosis are poorly understood. In a previous study, we demonstrated that overexpression of 7SK induces apoptosis in HEK cells. However, it remains unclear whether 7SK-mediated apoptosis induction is exerted through the intrinsic or extrinsic pathways.

METHODS AND RESULTS

Rn7SK was overexpressed in HEK 293T cell line using Lipofectamine 2000 reagent to investigate its potential apoptotic functions. The overexpression of Rn7SK resulted in reduced cell viability through the induction of apoptosis, as evidenced by MTT assay and Annexin V/PI staining. Concurrently, alterations in the expression levels of key apoptosis-related genes were observed, as determined by quantitative RT-PCR. Furthermore, Rn7SK overexpression led to a decrease in cell proliferation, as assessed by colony formation assay and growth curve analysis. This reduction was associated with downregulated expression of key proliferative-related genes. Additionally, the migration and invasion capabilities of cells were significantly inhibited upon upregulation of Rn7SK, as demonstrated by transwell assays.

CONCLUSIONS

This study suggests the apoptotic role of 7SK through both intrinsic and extrinsic pathways, necessitating further investigation into its underlying mechanisms.

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.

SNORD11B-mediated 2'-O-methylation of primary let-7a in colorectal carcinogenesis

Evidence indicates that small nucleolar RNAs (snoRNAs) participate in tumorigenesis and development and could be promising biomarkers for colorectal cancer (CRC). Here, we examine the profile of snoRNAs in CRC and find that expression of SNORD11B is increased in CRC tumor tissues and cell lines, with a significant positive correlation between SNORD11B expression and that of its host gene NOP58. SNORD11B promotes CRC cell proliferation and invasion and inhibits apoptosis. Mechanistically, SNORD11B promotes the processing and maturation of 18 S ribosomal RNA (rRNA) by mediating 2'-O-methylated (Nm) modification on the G509 site of 18 S rRNA. Intriguingly, SNORD11B mediates Nm modification on the G225 site of MIRLET7A1HG (pri-let-7a) with a canonical motif, resulting in degradation of pri-let-7a, inhibition of DGCR8 binding, reduction in mature tumor suppressor gene let-7a-5p expression, and upregulation of downstream oncogene translation. SNORD11B performs comparably to CEA and CA199 in diagnosing CRC. High expression of SNORD11B is significantly correlated with a more advanced TNM stage and lymph node metastasis, which indicates poor prognosis.

Expression and molecular regulation of non-coding RNAs in HPV-positive head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) is the sixth most prevalent malignancy worldwide. Accumulating evidence suggests that persistent HPV infection is closely related to a subset of HNSCC types, and the incidence of human papillomavirus (HPV)-positive HNSCC has been annually increasing in recent decades. Although the carcinogenesis of HPV-positive HNSCC has not been completely elucidated, it has been well confirmed that E6 and E7, the main viral oncoproteins are responsible for the maintenance of malignant transformation, promotion of cell proliferation, and increase in tumor invasion. Moreover, compared with HPV-negative HNSCC, HPV-positive HNSCC shows some special clinical-pathological features, which are possibly related to HPV infection and their specific regulatory mechanisms. Non-coding RNA (ncRNA) is a class of RNA lacking the protein-coding function and playing a critical regulatory role via multiple complex molecular mechanisms. NcRNA is an important regulatory pattern of epigenetic modification, which can exert significant effects on HPV-induced tumorigenesis and progression by deregulating downstream genes. However, the knowledge of ncRNAs is still limited, hence, a better understanding of ncRNAs could provide some insights for exploring the carcinogenesis mechanism and identifying valuable biomarkers in HPV-positive HNSCC. Therefore, in this review, we mainly focused on the expression profile of ncRNAs (including lncRNA, miRNA, and circRNA) and explored their regulatory role in HPV-positive HNSCC, aiming to clarify the regulatory mechanism of ncRNAs and identify valuable biomarkers for HPV-positive HNSCC.

Circulatory long noncoding RNAs (circulatory-LNC-RNAs) as novel biomarkers and therapeutic targets in cardiovascular diseases: Implications for cardiovascular diseases complications

Cardiovascular diseases (CVDs) are a group of disorders with major global health consequences. The prevalence of CVDs continues to grow due to population-aging and lifestyle modifications. Non-coding RNAs (ncRNAs) as key regulators of cell signaling pathways have gained attention in the occurrence and development of CVDs. Exosomal-lncRNAs (exos-lncRNAs) are emerging biomarkers due to their high sensitivity and specificity, stability, accuracy and accessibility in the biological fluids. Recently, circulatory and exos-based-lncRNAs are emerging and novel bio-tools in various pathogenic conditions. It is worth mentioning that dysregulation of these molecules has been found in different types of CVDs. In this regard, we aimed to discuss the knowledge gaps and suggest research priorities regarding circulatory and exos-lncRNAs as novel bio-tools and therapeutic targets for CVDs.

Research progress on plant noncoding RNAs in response to low-temperature stress

Low temperature (LT) is an important factor limiting plant growth and distribution. Plants have evolved sophisticated adaptive mechanisms to cope with hypothermia. RNA silencing is the orchestrator of these cellular responses. RNA silencing, which modifies gene expression through noncoding RNAs (ncRNAs), is a strategy used by plants to combat environmental stress. ncRNAs, which have very little protein-coding capacity, work by binding reverse complementary endogenous transcripts. In plants, ncRNAs include small non-coding RNAs (sncRNAs), medium-sized non-coding RNAs (mncRNAs), and long non-coding RNAs (lncRNAs). Apart from describing the biogenesis of different ncRNAs (miRNAs, siRNAs, and lncRNAs), we thoroughly discuss the functions of these ncRNAs during cold acclimation. Two major classes of sncRNAs, microRNAs and siRNAs, play essential regulatory roles in cold response processes through the posttranscriptional gene silencing (PTGS) pathway or transcriptional gene silencing (TGS) pathway. Microarray or transcriptome sequencing analysis can reveal a large number of cold-responsive miRNAs in plants. In this review, the cold-response patterns of miRNAs verified by Northern blotting or quantitative PCR in Arabidopsis thaliana, rice, and many other important crops are discussed. The detailed molecular mechanisms of several miRNAs in Arabidopsis (miR397, miR408, miR402, and miR394) and rice (Osa-miR156, Osa-miR319, and Osa-miR528) that regulate plant cold resistance are elucidated. In addition, the regulatory mechanism of the lncRNA SVALKA in the cold signaling pathway is explained in detail. Finally, we present the challenges for understanding the roles of small ncRNAs in cold signal transduction.

Midsize noncoding RNAs in cancers: a new division that clarifies the world of noncoding RNA or an unnecessary chaos?

Most of the human genome is made out of noncoding RNAs (ncRNAs). These ncRNAs do not code for proteins but carry a vast number of important functions in human cells such as: modification and processing other RNAs (tRNAs, rRNAs, snRNAs, snoRNAs, miRNAs), help in the synthesis of ribosome proteins, initiation of DNA replication, regulation of transcription, processing of pre-messenger mRNA during its maturation and much more. The ncRNAs also have a significant impact on many events that occur during carcinogenesis in cancer cells, such as: regulation of cell survival, cellular signaling, apoptosis, proliferation or even influencing the metastasis process. The ncRNAs may be divided based on their length, into short and long, where 200 nucleotides is the "magic" border. However, a new division was proposed, suggesting the creation of the additional group called midsize noncoding RNAs, with the length ranging from 50-400 nucleotides. This new group may include: transfer RNA (tRNA), small nuclear RNAs (snRNAs) with 7SK and 7SL, small nucleolar RNAs (snoRNAs), small Cajal body-specific RNAs (scaRNAs) and YRNAs. In this review their structure, biogenesis, function and influence on carcinogenesis process will be evaluated. What is more, a question will be answered of whether this new division is a necessity that clears current knowledge or just creates an additional misunderstanding in the ncRNA world?

CODA: a combo-Seq data analysis workflow

The analysis of the combined mRNA and miRNA content of a biological sample can be of interest for answering several research questions, like biomarkers discovery, or mRNA-miRNA interactions. However, the process is costly and time-consuming, separate libraries need to be prepared and sequenced on different flowcells. Combo-Seq is a library prep kit that allows us to prepare combined mRNA-miRNA libraries starting from very low total RNA. To date, no dedicated bioinformatics method exists for the processing of Combo-Seq data. In this paper, we describe CODA (Combo-seq Data Analysis), a workflow specifically developed for the processing of Combo-Seq data that employs existing free-to-use tools. We compare CODA with exceRpt, the pipeline suggested by the kit manufacturer for this purpose. We also evaluate how Combo-Seq libraries analysed with CODA perform compared with conventional poly(A) and small RNA libraries prepared from the same samples. We show that using CODA more successfully trimmed reads are recovered compared with exceRpt, and the difference is more dramatic with short sequencing reads. We demonstrate how Combo-Seq identifies as many genes and fewer miRNAs compared to the standard libraries, and how miRNA validation favours conventional small RNA libraries over Combo-Seq. The CODA code is available at https://github.com/marta-nazzari/CODA.

Quantum computing algorithms: getting closer to critical problems in computational biology

The recent biotechnological progress has allowed life scientists and physicians to access an unprecedented, massive amount of data at all levels (molecular, supramolecular, cellular and so on) of biological complexity. So far, mostly classical computational efforts have been dedicated to the simulation, prediction or de novo design of biomolecules, in order to improve the understanding of their function or to develop novel therapeutics. At a higher level of complexity, the progress of omics disciplines (genomics, transcriptomics, proteomics and metabolomics) has prompted researchers to develop informatics means to describe and annotate new biomolecules identified with a resolution down to the single cell, but also with a high-throughput speed. Machine learning approaches have been implemented to both the modelling studies and the handling of biomedical data. Quantum computing (QC) approaches hold the promise to resolve, speed up or refine the analysis of a wide range of these computational problems. Here, we review and comment on recently developed QC algorithms for biocomputing, with a particular focus on multi-scale modelling and genomic analyses. Indeed, differently from other computational approaches such as protein structure prediction, these problems have been shown to be adequately mapped onto quantum architectures, the main limit for their immediate use being the number of qubits and decoherence effects in the available quantum machines. Possible advantages over the classical counterparts are highlighted, along with a description of some hybrid classical/quantum approaches, which could be the closest to be realistically applied in biocomputation.

Discovery of RUF6 ncRNA-interacting proteins involved in P. falciparum immune evasion

Non-coding RNAs (ncRNAs) are emerging regulators of immune evasion and transmission of Plasmodium falciparum RUF6 is an ncRNA gene family that is transcribed by RNA polymerase III but actively regulates the Pol II-transcribed var virulence gene family. Understanding how RUF6 ncRNA connects to downstream effectors is lacking. We developed an RNA-directed proteomic discovery (ChIRP-MS) protocol to identify in vivo RUF6 ncRNA-protein interactions. The RUF6 ncRNA interactome was purified with biotinylated antisense oligonucleotides. Quantitative label-free mass spectrometry identified several unique proteins linked to gene transcription including RNA Pol II subunits, nucleosome assembly proteins, and a homologue of DEAD box helicase 5 (DDX5). Affinity purification of Pf-DDX5 identified proteins originally found by our RUF6-ChIRP protocol, validating the technique's robustness for identifying ncRNA interactomes in P. falciparum Inducible displacement of nuclear Pf-DDX5 resulted in significant down-regulation of the active var gene. Our work identifies a RUF6 ncRNA-protein complex that interacts with RNA Pol II to sustain the var gene expression, including a helicase that may resolve G-quadruplex secondary structures in var genes to facilitate transcriptional activation and progression.

Current Status of Regulatory Non-Coding RNAs Research in the Tritryp

Trypanosomatids are protozoan parasites that cause devastating vector-borne human diseases. Gene expression regulation of these organisms depends on post-transcriptional control in responding to diverse environments while going through multiple developmental stages of their complex life cycles. In this scenario, non-coding RNAs (ncRNAs) are excellent candidates for a very efficient, quick, and economic strategy to regulate gene expression. The advent of high throughput RNA sequencing technologies show the presence and deregulation of small RNA fragments derived from canonical ncRNAs. This review seeks to depict the ncRNA landscape in trypanosomatids, focusing on the small RNA fragments derived from functional RNA molecules observed in RNA sequencing studies. Small RNA fragments derived from canonical ncRNAs (tsRNAs, snsRNAs, sdRNAs, and sdrRNAs) were identified in trypanosomatids. Some of these RNAs display changes in their levels associated with different environments and developmental stages, demanding further studies to determine their functional characterization and potential roles. Nevertheless, a comprehensive and detailed ncRNA annotation for most trypanosomatid genomes is still needed, allowing better and more extensive comparative and functional studies.

Extracellular vesicles carry distinct proteo-transcriptomic signatures that are different from their cancer cell of origin

Circulating extracellular vesicles (EVs) contain molecular footprints-lipids, proteins, RNA, and DNA-from their cell of origin. Consequently, EV-associated RNA and proteins have gained widespread interest as liquid-biopsy biomarkers. Yet, an integrative proteo-transcriptomic landscape of EVs and comparison with their cell of origin remains obscure. Here, we report that EVs enrich distinct proteo-transcriptome that does not linearly correlate with their cell of origin. We show that EVs enrich endosomal and extracellular proteins, small RNA (∼13-200 nucleotides) associated with cell differentiation, development, and Wnt signaling. EVs cargo specific RNAs (RNY3, vtRNA, and MIRLET-7) and their complementary proteins (YBX1, IGF2BP2, and SRSF1/2). To ensure an unbiased and independent analyses, we studied 12 cancer cell lines, matching EVs (inhouse and exRNA database), and serum EVs of patients with prostate cancer. Together, we show that EV-RNA-protein complexes may constitute a functional interaction network to protect and regulate molecular access until a function is achieved.

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.

The tRNA-Derived Fragment tRF-24-V29K9UV3IU Functions as a miRNA-like RNA to Prevent Gastric Cancer Progression by Inhibiting GPR78 Expression

Emerging studies have proved that tRNA-derived fragments (tRFs) play vital roles in tumor metastasis; however, the function of tRFs in gastric cancer (GC) remains largely unclear. We investigated the role of tRF-24-V29K9UV3IU in growth and metastasis of GC using a xenograft mouse model. Differential gene expression downstream of tRF-24-V29K9UV3IU was identified by transcriptome sequencing, and interaction was then verified by a dual luciferase reporter and RNA immunoprecipitation. MKN-45 cells were also used to explore the biological functions of tRF-24-V29K9UV3IU in vitro. Here, knockdown of tRF-24-V29K9UV3IU promoted tumor growth and metastasis of GC in vivo. The expression of tRF-24-V29K9UV3IU and E-cadherin (epithelial cell marker) was down-regulated in tumors of mice following tRF-24-V29K9UV3IU knockdown, whereas the mesenchymal cell markers N-cadherin and vimentin displayed an opposite trend. Transcriptome sequencing identified 87 differentially expressed genes (DEGs) down-regulated in the tRF-24-V29K9UV3IU-overexpressed groups compared with the control group. Among them, G-protein–coupled receptor 78 (GPR78), the most significantly down-regulated DEG, was also predicted to be a target of tRF-24-V29K9UV3IU. Moreover, tRF-24-V29K9UV3IU could function as a miRNA-like fragment and bind to AGO2 and directly silence GPR78 expression by complementing with the 3′-untranslated region of the GPR78 mRNA. Functionally, overexpression of tRF-24-V29K9UV3IU significantly suppressed proliferation, migration, and invasion and promoted apoptosis of MKN-45 cells, whereas GPR78 attenuated these effects. Therefore, our data suggest that tRF-24-V29K9UV3IU functions as a miRNA-like fragment to suppress GPR78 expression and thus inhibit GC progression. These observations suggest that the tRF-24-V29K9UV3IU/GPR78 axis serves as a potential therapeutic target in GC.

Potential Roles of Long Noncoding RNAs as Therapeutic Targets in Organ Transplantation

Organ transplantation is the most preferred treatment option for end-stage organ diseases; however, allograft rejection is the major hurdle in successful long-term transplant survival. In spite of developing better HLA matching and more effective immunosuppressive regimen, one-year graft survival has been increased by nearly 90% and the incidence of acute rejection by one-year post-transplantation has been decreased by 12.2% in the last decades, chronic allograft rejection has remained as one of the major obstacles to the long-lasting survival of the transplanted allograft. Therefore, seemingly preventing the allograft rejection and inducing immunological tolerance against transplanted allografts is one of the primary goals in transplantation research to enable long-lasting graft survival. Various mechanisms such as long noncoding RNAs (lncRNAs) have been proposed that induce immune tolerance by modulating the gene expression and regulating innate and adaptive immune responses during transplantation. Besides, because of involvement in regulating epigenetic, transcriptional, and post-translational mechanisms, lncRNAs could affect allograft status. Therefore, these molecules could be considered as the potential targets for prediction, prognosis, diagnosis, and treatment of graft rejection. It is suggested that the noninvasive predictive biomarkers hold promise to overcome the current limitations of conventional tissue biopsy in the diagnosis of rejection. Hence, this review aims to provide a comprehensive overview of lncRNAs and their function to facilitate diagnosis, prognosis, and prediction of the risk of graft rejection, and the suggestive therapeutic choices after transplantation.

High-grade ovarian cancer associated H/ACA snoRNAs promote cancer cell proliferation and survival

Small nucleolar RNAs (snoRNAs) are an omnipresent class of non-coding RNAs involved in the modification and processing of ribosomal RNA (rRNA). As snoRNAs are required for ribosome production, the increase of which is a hallmark of cancer development, their expression would be expected to increase in proliferating cancer cells. However, assessing the nature and extent of snoRNAs' contribution to cancer biology has been largely limited by difficulties in detecting highly structured RNA. In this study, we used a dedicated midsize non-coding RNA (mncRNA) sensitive sequencing technique to accurately survey the snoRNA abundance in independently verified high-grade serous ovarian carcinoma (HGSC) and serous borderline tumour (SBT) tissues. The results identified SNORA81, SNORA19 and SNORA56 as an H/ACA snoRNA signature capable of discriminating between independent sets of HGSC, SBT and normal tissues. The expression of the signature SNORA81 correlates with the level of ribosomal RNA (rRNA) modification and its knockdown inhibits 28S rRNA pseudouridylation and accumulation leading to reduced cell proliferation and migration. Together our data indicate that specific subsets of H/ACA snoRNAs may promote tumour aggressiveness by inducing rRNA modification and synthesis.

The Emerging Landscape of Long Non-Coding RNAs in Wilms Tumor

Long noncoding RNAs (LncRNAs) are transcripts of nucleic acid sequences with a length of more than 200 bp, which have only partial coding capabilities. Recent studies have shown that lncRNAs located in the nucleus or cytoplasm can be used as gene expression regulatory elements due to their important regulatory effects in a variety of biological processes. Wilms tumor (WT) is a common abdominal tumor in children whose pathogenesis remains unclear. In recent years, many specifically expressed lncRNAs have been found in WT, which affect the occurrence and development of WT. At the same time, lncRNAs may have the capacity to become novel biomarkers for the diagnosis and prognosis of WT. This article reviews related research progress on the relationship between lncRNAs and WT, to provide a new direction for clinical diagnosis and treatment of WT.

MGcount: a total RNA-seq quantification tool to address multi-mapping and multi-overlapping alignments ambiguity in non-coding transcripts

BACKGROUND

Total-RNA sequencing (total-RNA-seq) allows the simultaneous study of both the coding and the non-coding transcriptome. Yet, computational pipelines have traditionally focused on particular biotypes, making assumptions that are not fullfilled by total-RNA-seq datasets. Transcripts from distinct RNA biotypes vary in length, biogenesis, and function, can overlap in a genomic region, and may be present in the genome with a high copy number. Consequently, reads from total-RNA-seq libraries may cause ambiguous genomic alignments, demanding for flexible quantification approaches.

RESULTS

Here we present Multi-Graph count (MGcount), a total-RNA-seq quantification tool combining two strategies for handling ambiguous alignments. First, MGcount assigns reads hierarchically to small-RNA and long-RNA features to account for length disparity when transcripts overlap in the same genomic position. Next, MGcount aggregates RNA products with similar sequences where reads systematically multi-map using a graph-based approach. MGcount outputs a transcriptomic count matrix compatible with RNA-sequencing downstream analysis pipelines, with both bulk and single-cell resolution, and the graphs that model repeated transcript structures for different biotypes. The software can be used as a python module or as a single-file executable program.

CONCLUSIONS

MGcount is a flexible total-RNA-seq quantification tool that successfully integrates reads that align to multiple genomic locations or that overlap with multiple gene features. Its approach is suitable for the simultaneous estimation of protein-coding, long non-coding and small non-coding transcript concentration, in both precursor and processed forms. Both source code and compiled software are available at https://github.com/hitaandrea/MGcount .

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07757-1.

Emerging Data on the Diversity of Molecular Mechanisms Involving C/D snoRNAs

Box C/D small nucleolar RNAs (C/D snoRNAs) represent an ancient family of small non-coding RNAs that are classically viewed as housekeeping guides for the 2′-O-methylation of ribosomal RNA in Archaea and Eukaryotes. However, an extensive set of studies now argues that they are involved in mechanisms that go well beyond this function. Here, we present these pieces of evidence in light of the current comprehension of the molecular mechanisms that control C/D snoRNA expression and function. From this inventory emerges that an accurate description of these activities at a molecular level is required to let the snoRNA field enter in a second age of maturity.

Pan-cancer chromatin analysis of the human vtRNA genes uncovers their association with cancer biology

Background: The vault RNAs (vtRNAs) are a class of 84-141-nt eukaryotic non-coding RNAs transcribed by RNA polymerase III, associated to the ribonucleoprotein complex known as vault particle. Of the four human vtRNA genes, vtRNA1-1, vtRNA1-2 and vtRNA1-3, clustered at locus 1, are integral components of the vault particle, while vtRNA2-1 is a more divergent homologue located in a second locus. Gene expression studies of vtRNAs in large cohorts have been hindered by their unsuccessful sequencing using conventional transcriptomic approaches.

Methods: VtRNA expression in The Cancer Genome Atlas (TCGA) Pan-Cancer cohort was estimated using the genome-wide DNA methylation and chromatin accessibility data (ATAC-seq) of their genes as surrogate variables. The association between vtRNA expression and patient clinical outcome, immune subtypes and transcriptionally co-regulated gene programs was analyzed in the dataset.

Results: VtRNAs promoters are enriched in transcription factors related to viral infection. VtRNA2-1 is likely the most independently regulated homologue. VtRNA1-1 has the most accessible chromatin, followed by vtRNA1-2, vtRNA2-1 and vtRNA1-3. VtRNA1-1 and vtRNA1-3 chromatin status does not significantly change in cancer tissues. Meanwhile, vtRNA2-1 and vtRNA1-2 expression is widely deregulated in neoplastic tissues and its alteration is compatible with a broad oncogenic role for vtRNA1-2, and both tumor suppressor and oncogenic functions for vtRNA2-1. Yet, vtRNA1-1, vtRNA1-2 and vtRNA2-1 promoter DNA methylation predicts a shorter patient overall survival cancer-wide. In addition, gene ontology analyses of vtRNAs co-regulated genes identify a chromosome regulatory domain, epithelial differentiation, immune and thyroid cancer gene sets for specific vtRNAs. Furthermore, vtRNA expression patterns are associated with cancer immune subtypes and vtRNA1-2 expression is positively associated with cell proliferation and wound healing.

Conclusions: Our study presents the landscape of vtRNA chromatin status cancer-wide, identifying co-regulated gene networks and ontological pathways associated with the different vtRNA genes that may account for their diverse roles in cancer.

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.

The medium-size noncoding RNA transcriptome of Ostreococcus tauri, the smallest living eukaryote, reveals a large family of small nucleolar RNAs displaying multiple genomic expression strategies

The small nucleolar RNAs (snoRNAs), essential for ribosome biogenesis, constitute a major family of medium-size noncoding RNAs (mncRNAs) in all eukaryotes. We present here, for the first time in a marine unicellular alga, the characterization of the snoRNAs family in Ostreococcus tauri, the smallest photosynthetic eukaryote. Using a transcriptomic approach, we identified 131 O. tauri snoRNAs (Ot–snoRNA) distributed in three classes: the C/D snoRNAs, the H/ACA snoRNAs and the MRP RNA. Their genomic organization revealed a unique combination of both the intronic organization of animals and the polycistronic organization of plants. Remarkably, clustered genes produced Ot–snoRNAs with unusual structures never previously described in plants. Their abundances, based on quantification of reads and northern blots, showed extreme differences in Ot–snoRNA accumulation, mainly determined by their differential stability. Most of these Ot–snoRNAs were predicted to target rRNAs or snRNAs. Seventeen others were orphan Ot–snoRNAs that would not target rRNA. These were specific to O. tauri or Mamiellophyceae and could have functions unrelated to ribosome biogenesis. Overall, these data reveal an ‘evolutionary response’ adapted to the extreme compactness of the O. tauri genome that accommodates the essential Ot–snoRNAs, developing multiple strategies to optimize their coordinated expression with a minimal cost on regulatory circuits.

One locus with two roles: microRNA-independent functions of microRNA-host-gene locus-encoded long noncoding RNAs

Long noncoding RNAs (lncRNAs) are RNA transcripts longer than 200 nucleotides that do not code for proteins. LncRNAs play crucial regulatory roles in several biological processes via diverse mechanisms and their aberrant expression is associated with various diseases. LncRNA genes are further subcategorized based on their relative organization in the genome. MicroRNA (miRNA)-host-gene-derived lncRNAs (lnc-MIRHGs) refer to lncRNAs whose genes also harbor miRNAs. There exists crosstalk between the processing of lnc-MIRHGs and the biogenesis of the encoded miRNAs. Although the functions of the encoded miRNAs are usually well understood, whether those lnc-MIRHGs play independent functions are not fully elucidated. Here, we review our current understanding of lnc-MIRHGs, including their biogenesis, function, and mechanism of action, with a focus on discussing the miRNA-independent functions of lnc-MIRHGs, including their involvement in cancer. Our current understanding of lnc-MIRHGs strongly indicates that this class of lncRNAs could play important roles in basic cellular events as well as in diseases. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

In-Depth Analysis Reveals Production of Circular RNAs from Non-Coding Sequences

The sequencing of total RNA depleted for ribosomal sequences remains the method of choice for the study of circRNAs. Our objective was to characterize non-canonical circRNAs, namely not originating from back splicing and circRNA produced by non-coding genes. To this end, we analyzed a dataset from porcine testis known to contain about 100 intron-derived circRNAs. Labelling reads containing a circular junction and originating from back splicing provided information on the very small contribution of long non-coding genes to the production of canonical circRNAs. Analyses of the other reads revealed two origins for non-canonical circRNAs: (1) Intronic sequences for lariat-derived intronic circRNAs and intron circles, (2) Mono-exonic genes (mostly non-coding) for either a new type of circRNA (including only part of the exon: sub-exonic circRNAs) or, even more rarely, mono-exonic canonical circRNAs. The most complex set of sub-exonic circRNAs was produced by RNase_MRP (ribozyme RNA). We specifically investigated the intronic circRNA of ATXN2L, which is probably an independently transcribed sisRNA (stable intronic sequence RNA). We may be witnessing the emergence of a new non-coding gene in the porcine genome. Our results are evidence that most non-canonical circRNAs originate from non-coding sequences.

Direct Determination of Pseudouridine in RNA by Mass Spectrometry Coupled with Stable Isotope Labeling

Pseudouridine (Ψ) is the only "mass-silent" nucleoside produced by post-transcriptional RNA modification. We developed a mass spectrometry (MS)-based technique coupled with in vivo deuterium (D) labeling of uridines for direct determination of Ψs in cellular RNA and applied it to the comprehensive analysis of post-transcriptional modifications in human ribosomal RNAs. The method utilizes human TK6/mouse FM3A cells deficient in uridine monophosphate synthase using a CRISPR-Cas9 technique to turn off de novo uridine synthesis and fully labels uridines with D at uracil positions 5 and 6 by cultivating the cells in a medium containing uridine-5,6-D₂. The pseudouridylation reaction in those cells results in the exchange of the D at the C5 of uracil with hydrogen from solvent, which produces a -1 Da mass shift, thus allowing MS-based determination of RNA Ψs. We present here the experimental details of this method and show that it allows the identification of all Ψs in human major nuclear and nucleolar RNAs, including several previously unknown Ψs. Because the method allows direct determination of Ψs at the femtomole level of RNA, it will serve as a useful tool for structure/function studies of a wide variety of noncoding RNAs.

Handling multi-mapped reads in RNA-seq

Many eukaryotic genomes harbour large numbers of duplicated sequences, of diverse biotypes, resulting from several mechanisms including recombination, whole genome duplication and retro-transposition. Such repeated sequences complicate gene/transcript quantification during RNA-seq analysis due to reads mapping to more than one locus, sometimes involving genes embedded in other genes. Genes of different biotypes have dissimilar levels of sequence duplication, with long-noncoding RNAs and messenger RNAs sharing less sequence similarity to other genes than biotypes encoding shorter RNAs. Many strategies have been elaborated to handle these multi-mapped reads, resulting in increased accuracy in gene/transcript quantification, although separate tools are typically used to estimate the abundance of short and long genes due to their dissimilar characteristics. This review discusses the mechanisms leading to sequence duplication, the biotypes affected, the computational strategies employed to deal with multi-mapped reads and the challenges that still remain to be overcome.

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.

Introns: Good Day Junk Is Bad Day Treasure

Introns are ubiquitous in eukaryotic transcripts. They are often viewed as junk RNA but the huge energetic burden of transcribing, removing, and degrading them suggests a significant evolutionary advantage. Ostensibly, an intron functions within the host pre-mRNA to regulate its splicing, transport, and degradation. However, recent studies have revealed an entirely new class of trans-acting functions where the presence of intronic RNA in the cell impacts the expression of other genes in trans. Here, we review possible new mechanisms of intron functions, with a focus on the role of yeast introns in regulating the cell growth response to starvation.

The cellular landscape of mid-size noncoding RNA

Noncoding RNA plays an important role in all aspects of the cellular life cycle, from the very basic process of protein synthesis to specialized roles in cell development and differentiation. However, many noncoding RNAs remain uncharacterized and the function of most of them remains unknown. Mid-size noncoding RNAs (mncRNAs), which range in length from 50 to 400 nucleotides, have diverse regulatory functions but share many fundamental characteristics. Most mncRNAs are produced from independent promoters although others are produced from the introns of other genes. Many are found in multiple copies in genomes. mncRNAs are highly structured and carry many posttranscriptional modifications. Both of these facets dictate their RNA-binding protein partners and ultimately their function. mncRNAs have already been implicated in translation, catalysis, as guides for RNA modification, as spliceosome components and regulatory RNA. However, recent studies are adding new mncRNA functions including regulation of gene expression and alternative splicing. In this review, we describe the different classes, characteristics and emerging functions of mncRNAs and their relative expression patterns. Finally, we provide a portrait of the challenges facing their detection and annotation in databases. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.