IsoMiRmap: fast, deterministic and exhaustive mining of isomiRs from short RNA-seq datasets

Comprehensive profiling of small RNAs and their changes and linkages to mRNAs in schizophrenia and bipolar disorder

We investigated small non-coding RNAs (sncRNAs) from the prefrontal cortex of 93 individuals diagnosed with schizophrenia (SCZ) or bipolar disorder (BD) and 77 controls. We uncovered recurring complex sncRNA profiles, with 98% of all sncRNAs being accounted for by miRNA isoforms (60.6%), tRNA-derived fragments (17.8%), rRNA-derived fragments (11.4%), and Y RNA-derived fragments (8.3%). In SCZ, 15% of all sncRNAs exhibit statistically significant changes in their abundance. In BD, the fold changes (FCs) are highly correlated with those in SCZ but less acute. Non-templated nucleotide additions to the 3´-ends of many miRNA isoforms determine their FC independently of miRNA identity or genomic locus of origin. In both SCZ and BD, disease- and age-associated sncRNAs and mRNAs reveal accelerated aging. Co-expression modules between sncRNAs and mRNAs align with the polarities of SCZ changes and implicate sncRNAs in critical processes, including synaptic signaling, neurogenesis, memory, behavior, and cognition.

The 2024 Nobel Prize in Physiology or Medicine: microRNA Takes Center Stage

The Non-coding Journal Editorial Board Members would like to congratulate Victor Ambros and Gary Ruvkun, who were jointly awarded the 2024 Nobel Prize in Physiology or Medicine for their groundbreaking discovery of microRNAs and the role of microRNAs in post-transcriptional gene regulation, uncovering a previously unknown layer of gene control in eukaryotes [...].

athisomiRDB: A comprehensive database of Arabidopsis isomiRs

Several pieces of evidence challenge the traditional view of miRNAs as static molecules, revealing dynamic isomiRs originating from each miRNA precursor arm. In plants, isomiRs, which result from imprecise cleavage during pre-miRNA processing and post-transcriptional alterations, serve as crucial regulators of target microRNAs (miRNAs). Despite numerous studies on Arabidopsis miRNAs, the systematic identification and annotation of isomiRs across various tissues and conditions remain limited. Due to the lack of systematically collected isomiR information, we introduce the athisomiRDB database, which houses 20 764 isomiRs from Arabidopsis small RNA-sequencing (sRNA-seq) libraries. It comprises >2700 diverse samples and allows exploration at the sample, miRNA, or isomiR levels, offering insights into the presence or absence of isomiRs. The athisomiRDB includes exclusive and ambiguous isomiRs, each with features such as transcriptional origin, variant-containing isomiRs, and identifiers for frequent single-nucleotide polymorphisms from the 1001 Genomes Project. It also provides 3' nontemplated post-transcriptional additions, isomiR-target interactions, and trait associations for each isomiR. We anticipate that athisomiRDB will play a pivotal role in unraveling the regulatory nature of the Arabidopsis miRNAome and enhancing sRNA research by leveraging isomiR profiles from extensive sRNA-seq datasets. Database URL: https://www.nipgr.ac.in/athisomiRDB.

The subcellular distribution of miRNA isoforms, tRNA-derived fragments, and rRNA-derived fragments depends on nucleotide sequence and cell type

BACKGROUND

MicroRNA isoforms (isomiRs), tRNA-derived fragments (tRFs), and rRNA-derived fragments (rRFs) represent most of the small non-coding RNAs (sncRNAs) found in cells. Members of these three classes modulate messenger RNA (mRNA) and protein abundance and are dysregulated in diseases. Experimental studies to date have assumed that the subcellular distribution of these molecules is well-understood, independent of cell type, and the same for all isoforms of a sncRNA.

RESULTS

We tested these assumptions by investigating the subcellular distribution of isomiRs, tRFs, and rRFs in biological replicates from three cell lines from the same tissue and same-sex donors that model the same cancer subtype. In each cell line, we profiled the isomiRs, tRFs, and rRFs in the nucleus, cytoplasm, whole mitochondrion (MT), mitoplast (MP), and whole cell. Using a rigorous mathematical model we developed, we accounted for cross-fraction contamination and technical errors and adjusted the measured abundances accordingly. Analyses of the adjusted abundances show that isomiRs, tRFs, and rRFs exhibit complex patterns of subcellular distributions. These patterns depend on each sncRNA's exact sequence and the cell type. Even in the same cell line, isoforms of the same sncRNA whose sequences differ by a few nucleotides (nts) can have different subcellular distributions.

CONCLUSIONS

SncRNAs with similar sequences have different subcellular distributions within and across cell lines, suggesting that each isoform could have a different function. Future computational and experimental studies of isomiRs, tRFs, and rRFs will need to distinguish among each molecule's various isoforms and account for differences in each isoform's subcellular distribution in the cell line at hand. While the findings add to a growing body of evidence that isomiRs, tRFs, rRFs, tRNAs, and rRNAs follow complex intracellular trafficking rules, further investigation is needed to exclude alternative explanations for the observed subcellular distribution of sncRNAs.

Incautious design of shRNAs for stable overexpression of miRNAs could result in generation of undesired isomiRs

shRNA-mediated strategy of miRNA overexpression based on RNA Polymerase III (Pol III) expression cassettes is widely used for miRNA functional studies. For some miRNAs, e.g., encoded in the genome as a part of a polycistronic miRNA cluster, it is most likely the only way for their individual stable overexpression. Here we have revealed that expression of miRNAs longer than 19 nt (e.g. 23 nt in length hsa-miR-93-5p) using such approach could be accompanied by undesired predominant generation of 5' end miRNA isoforms (5'-isomiRs). Extra U residues (up to five) added by Pol III at the 3' end of the transcribed shRNA during transcription termination could cause a shift in the Dicer cleavage position of the shRNA. This results in the formation of 5'-isomiRs, which have a significantly altered seed region compared to the initially encoded canonical hsa-miR-93-5p. We demonstrated that the commonly used qPCR method is insensitive to the formation of 5'-isomiRs and cannot be used to confirm miRNA overexpression. However, the predominant expression of 5'-isomiRs without three or four first nucleotides instead of the canonical isoform could be disclosed based on miRNA-Seq analysis. Moreover, mRNA sequencing data showed that the 5'-isomiRs of hsa-miR-93-5p presumably regulate their own mRNA targets. Thus, omitting miRNA-Seq analysis may lead to erroneous conclusions regarding revealed mRNA targets and possible molecular mechanisms in which studied miRNA is involved. Overall, the presented results show that structures of shRNAs for stable overexpression of miRNAs requires careful design to avoid generation of undesired 5'-isomiRs.

The intricacies of isomiRs: from classification to clinical relevance

MicroRNAs (miRNAs) and isoforms of their archetype, called isomiRs, regulate gene expression via complementary base-pair binding to messenger RNAs (mRNAs). The partially evolutionarily conserved isomiR sequence variations are differentially expressed among tissues, populations, and genders, and between healthy and diseased states. Aiming towards the clinical use of isomiRs as diagnostic biomarkers and for therapeutic purposes, several challenges need to be addressed, including (i) clarification of isomiR definition, (ii) improved annotation in databases with new standardization (such as the mirGFF3 format), and (iii) improved methods of isomiR detection, functional verification, and in silico analysis. In this review we discuss the respective challenges, and highlight the opportunities for clinical use of isomiRs, especially in the light of increasing amounts of next-generation sequencing (NGS) data.

MINRbase: a comprehensive database of nuclear- and mitochondrial-ribosomal-RNA-derived fragments (rRFs)

We describe the Mitochondrial and Nuclear rRNA fragment database (MINRbase), a knowledge repository aimed at facilitating the study of ribosomal RNA-derived fragments (rRFs). MINRbase provides interactive access to the profiles of 130 238 expressed rRFs arising from the four human nuclear rRNAs (18S, 5.8S, 28S, 5S), two mitochondrial rRNAs (12S, 16S) or four spacers of 45S pre-rRNA. We compiled these profiles by analyzing 11 632 datasets, including the GEUVADIS and The Cancer Genome Atlas (TCGA) repositories. MINRbase offers a user-friendly interface that lets researchers issue complex queries based on one or more criteria, such as parental rRNA identity, nucleotide sequence, rRF minimum abundance and metadata keywords (e.g. tissue type, disease). A 'summary' page for each rRF provides a granular breakdown of its expression by tissue type, disease, sex, ancestry and other variables; it also allows users to create publication-ready plots at the click of a button. MINRbase has already allowed us to generate support for three novel observations: the internal spacers of 45S are prolific producers of abundant rRFs; many abundant rRFs straddle the known boundaries of rRNAs; rRF production is regimented and depends on 'personal attributes' (sex, ancestry) and 'context' (tissue type, tissue state, disease). MINRbase is available at https://cm.jefferson.edu/MINRbase/.

The Typical tRNA Co-Expresses Multiple 5' tRNA Halves Whose Sequences and Abundances Depend on Isodecoder and Isoacceptor and Change with Tissue Type, Cell Type, and Disease

Transfer RNA-derived fragments (tRFs) are noncoding RNAs that arise from either mature transfer RNAs (tRNAs) or their precursors. One important category of tRFs comprises the tRNA halves, which are generated through cleavage at the anticodon. A given tRNA typically gives rise to several co-expressed 5'-tRNA halves (5'-tRHs) that differ in the location of their 3' ends. These 5'-tRHs, even though distinct, have traditionally been treated as indistinguishable from one another due to their near-identical sequences and lengths. We focused on co-expressed 5'-tRHs that arise from the same tRNA and systematically examined their exact sequences and abundances across 10 different human tissues. To this end, we manually curated and analyzed several hundred human RNA-seq datasets from NCBI's Sequence Run Archive (SRA). We grouped datasets from the same tissue into their own collection and examined each group separately. We found that a given tRNA produces different groups of co-expressed 5'-tRHs in different tissues, different cell lines, and different diseases. Importantly, the co-expressed 5'-tRHs differ in their sequences, absolute abundances, and relative abundances, even among tRNAs with near-identical sequences from the same isodecoder or isoacceptor group. The findings suggest that co-expressed 5'-tRHs that are produced from the same tRNA or closely related tRNAs have distinct, context-dependent roles. Moreover, our analyses show that cell lines modeling the same tissue type and disease may not be interchangeable when it comes to experimenting with tRFs.

Characterization of all small RNAs in and comparisons across cultured megakaryocytes and platelets of healthy individuals and COVID-19 patients

BACKGROUND

The small noncoding RNAs (sncRNAs) in megakaryocytes (MKs) and platelets are not well characterized. Neither is the impact of SARS-CoV-2 infection on the sncRNAs of platelets.

OBJECTIVES

To investigate the sorting of MK sncRNAs into platelets, and the differences in the platelet sncRNAomes of healthy donors (HDs) and COVID-19 patients.

METHODS

We comprehensively profiled sncRNAs from MKs cultured from cord blood-derived CD34+ cells, platelets from HDs, and platelets from patients with moderate and severe SARS-CoV-2 infection. We also comprehensively profiled Argonaute (AGO)-bound sncRNAs from the cultured MKs.

RESULTS

We characterized the sncRNAs in MKs and platelets and can account for ∼95% of all sequenced reads. We found that MKs primarily comprise microRNA isoforms (isomiRs), tRNA-derived fragments (tRFs), rRNA-derived fragments (rRFs), and Y RNA-derived fragments (yRFs) in comparable abundances. The platelets of HDs showed a skewed distribution by comparison: 56.7% of all sncRNAs are yRFs, 34.4% are isomiRs, and <2.0% are tRFs and rRFs. Most isomiRs in MKs and platelets are either noncanonical, nontemplated, or both. When comparing MKs and platelets from HDs, we found numerous isomiRs, tRFs, rRFs, and yRFs showing opposite enrichments or depletions, including molecules from the same parental miRNA arm, tRNA, rRNA, or Y RNA. The sncRNAome of platelets from patients with COVID-19 is skewed compared to that of HDs with only 19.8% of all sncRNAs now being yRFs, isomiRs increasing to 63.6%, and tRFs and rRFs more than tripling their presence to 6.1%.

CONCLUSION

The sncRNAomes of MKs and platelets are very rich and more complex than it has been believed. The evidence suggests complex mechanisms that sort MK sncRNAs into platelets. SARS-CoV-2 infection acutely alters the contents of platelets by changing the relative proportions of their sncRNAs.

A novel approach for a joint analysis of isomiR and mRNA expression data reveals features of isomiR targeting in breast cancer

A widely used procedure for selecting significant miRNA-mRNA or isomiR-mRNA pairs out of predicted interactions involves calculating the correlation between expression levels of miRNAs/isomiRs and mRNAs in a series of samples. In this manuscript, we aimed to assess the validity of this procedure by comparing isomiR-mRNA correlation profiles in sets of sequence-based predicted target mRNAs and non-target mRNAs (negative controls). Target prediction was carried out using RNA22 and TargetScan algorithms. Spearman's correlation analysis was conducted using miRNA and mRNA sequencing data of The Cancer Genome Atlas Breast Invasive Carcinoma (TCGA-BRCA) project. Luminal A, luminal B, basal-like breast cancer subtypes, and adjacent normal tissue samples were analyzed separately. Using the sets of putative targets and non-targets, we introduced adjusted isomiR targeting activity (ITA)-the number of negatively correlated potential isomiR targets adjusted by the background (estimated using non-target mRNAs). We found that for most isomiRs a significant negative correlation between isomiR-mRNA expression levels appeared more often in a set of predicted targets compared to the non-targets. This trend was detected for both classical seed region binding types (8mer, 7mer-m8, 7mer-A1, 6mer) predicted by TargetScan and the non-classical ones (G:U wobbles and up to one mismatch or unpaired nucleotide within seed sequence) predicted by RNA22. Adjusted ITA distributions were similar for target sites located in 3'-UTRs and coding mRNA sequences, while 5'-UTRs had much lower scores. Finally, we observed strong cancer subtype-specific patterns of isomiR activity, highlighting the differences between breast cancer molecular subtypes and normal tissues. Surprisingly, our target prediction- and correlation-based estimates of isomiR activities were practically non-correlated with the average isomiR expression levels neither in cancerous nor in normal samples.

Pan-Cancer Analysis of Canonical and Modified miRNAs Enhances the Resolution of the Functional miRNAome in Cancer

Epitranscriptomic studies of miRNAs have added a new layer of complexity to the cancer field. Although there is fast-growing interest in adenosine-to-inosine (A-to-I) miRNA editing and alternative cleavage that shifts miRNA isoforms, simultaneous evaluation of both modifications in cancer is still missing. Here, we concurrently profiled multiple miRNA modification types, including A-to-I miRNA editing and shifted miRNA isoforms, in >13,000 adult and pediatric tumor samples across 38 distinct cancer cohorts from The Cancer Genome Atlas and The Therapeutically Applicable Research to Generate Effective Treatments data sets. The differences between canonical miRNAs and the wider miRNAome in terms of expression, clustering, dysregulation, and prognostic standpoint were investigated. The combination of canonical miRNAs and modified miRNAs boosted the quality of clustering results, outlining unique clinicopathologic features among cohorts. Certain modified miRNAs showed opposite expression from their canonical counterparts in cancer, potentially impacting their targets and function. Finally, a shifted and edited miRNA isoform was experimentally validated to directly bind and suppress a unique target. These findings outline the importance of going beyond the well-established paradigm of one mature miRNA per miRNA arm to elucidate novel mechanisms related to cancer progression.

SIGNIFICANCE

Modified miRNAs may act as cancer biomarkers and function as allies or antagonists of their canonical counterparts in gene regulation, suggesting the concurrent consideration of canonical and modified miRNAs can boost patient stratification.

Circulating microRNA signatures associated with disease severity and outcome in COVID-19 patients

Background

SARS-CoV-2 induces a spectrum of clinical conditions ranging from asymptomatic infection to life threatening severe disease. Host microRNAs have been involved in the cytokine storm driven by SARS-CoV-2 infection and proposed as candidate biomarkers for COVID-19.

Methods

To discover signatures of circulating miRNAs associated with COVID-19, disease severity and mortality, small RNA-sequencing was performed on serum samples collected from 89 COVID-19 patients (34 severe, 29 moderate, 26 mild) at hospital admission and from 45 healthy controls (HC). To search for possible sources of miRNAs, investigation of differentially expressed (DE) miRNAs in relevant human cell types in vitro.

Results

COVID-19 patients showed upregulation of miRNAs associated with lung disease, vascular damage and inflammation and downregulation of miRNAs that inhibit pro-inflammatory cytokines and chemokines, angiogenesis, and stress response. Compared with mild/moderate disease, patients with severe COVID-19 had a miRNA signature indicating a profound impairment of innate and adaptive immune responses, inflammation, lung fibrosis and heart failure. A subset of the DE miRNAs predicted mortality. In particular, a combination of high serum miR-22-3p and miR-21-5p, which target antiviral response genes, and low miR-224-5p and miR-155-5p, targeting pro-inflammatory factors, discriminated severe from mild/moderate COVID-19 (AUROC 0.88, 95% CI 0.80-0.95, p<0.0001), while high leukocyte count and low levels of miR-1-3p, miR-23b-3p, miR-141-3p, miR-155-5p and miR-4433b-5p predicted mortality with high sensitivity and specificity (AUROC 0.95, 95% CI 0.89-1.00, p<0.0001). In vitro experiments showed that some of the DE miRNAs were modulated directly by SARS-CoV-2 infection in permissive lung epithelial cells.

Conclusions

We discovered circulating miRNAs associated with COVID-19 severity and mortality. The identified DE miRNAs provided clues on COVID-19 pathogenesis, highlighting signatures of impaired interferon and antiviral responses, inflammation, organ damage and cardiovascular failure as associated with severe disease and death.

Epigenomic Modifications in Modern and Ancient Genomes

Epigenetic changes have been identified as a major driver of fundamental metabolic pathways. More specifically, the importance of epigenetic regulatory mechanisms for biological processes like speciation and embryogenesis has been well documented and revealed the direct link between epigenetic modifications and various diseases. In this review, we focus on epigenetic changes in animals with special attention on human DNA methylation utilizing ancient and modern genomes. Acknowledging the latest developments in ancient DNA research, we further discuss paleoepigenomic approaches as the only means to infer epigenetic changes in the past. Investigating genome-wide methylation patterns of ancient humans may ultimately yield in a more comprehensive understanding of how our ancestors have adapted to the changing environment, and modified their lifestyles accordingly. We discuss the difficulties of working with ancient DNA in particular utilizing paleoepigenomic approaches, and assess new paleoepigenomic data, which might be helpful in future studies.

MicroRNA Isoforms Contribution to Melanoma Pathogenesis

Cutaneous melanoma (CM) is the most lethal tumor among skin cancers, and its incidence is constantly increasing. A deeper understanding of the molecular processes guiding melanoma pathogenesis could improve diagnosis, treatment and prognosis. MicroRNAs play a key role in melanoma biology. Recently, next generation sequencing (NGS) experiments, designed to assess small-RNA expression, revealed the existence of microRNA variants with different length and sequence. These microRNA isoforms are known as isomiRs and provide an additional layer to the complex non-coding RNA world. Here, we collected data from NGS experiments to provide a comprehensive characterization of miRNA and isomiR dysregulation in benign nevi (BN) and early-stage melanomas. We observed that melanoma and BN express different and specific isomiRs and have a different isomiR abundance distribution. Moreover, isomiRs from the same microRNA can have opposite expression trends between groups. Using The Cancer Genome Atlas (TCGA) dataset of skin cancers, we analyzed isomiR expression in primary melanoma and melanoma metastasis and tested their association with NF1, BRAF and NRAS mutations. IsomiRs differentially expressed were identified and catalogued with reference to the canonical form. The reported non-random dysregulation of specific isomiRs contributes to the understanding of the complex melanoma pathogenesis and serves as the basis for further functional studies.