Beyond the one-locus-one-miRNA paradigm: microRNA isoforms enable deeper insights into breast cancer heterogeneity

MiR-30a-5p isoform -1|1 promotes the progression of gastric cancer by inhibiting TMEM66 and reducing intratumoral cytotoxic T cells

Gastric cancer is histologically classified into the intestinal subtype, which forms tubular structures, and the aggressive diffuse subtype, characterized by rapid invasion and poor prognosis. The variety and quantity of miRNA isoforms between different histological subtypes of gastric cancer were unknown. Through systematic filtering, we found that more diverse miR-30a-5p isoforms was present in the diffuse subtype of gastric cancer, and was associated with patients' worse survival independent of tumor stage based on the TCGA miRNA-seq data. Among all nine isoforms of miR-30a-5p, miR-30a-5p -1|1 was more abundant than the archetype of miR-30a-5p. Higher expression of miR-30a-5p -1|1 was observed in patients with advanced tumor stage and poor survival. Furthermore, miR-30a-5p -1|1 could promote the metastasis of gastric cancer cells both in vitro and in vivo by down-regulating TMEM66. In clinical samples, decreased expression of TMEM66 was characteristic of gastric cancer, and the low level of TMEM66 correlated with deceased CD8 positive cells in the tumor microenvironment probably due to decreased cytokines production. In conclusion, the variety of miR-30a-5p isoforms correlates with worse survival in gastric cancer patients. Moreover, miR-30a-5p -1|1 could promote gastric cancer metastasis by inhibiting TMEM66 and the infiltration of intratumoral CD8 positive cells.

5'-Isoforms of miR-1246 Have Distinct Targets and Stronger Functional Impact Compared with Canonical miR-1246 in Colorectal Cancer Cells In Vitro

Colorectal cancer (CRC) is a multifactorial disease involving genetic and epigenetic factors, such as miRNAs. Sequencing-based studies have revealed that miRNAs have many isoforms (isomiRs) with modifications at the 3'- and 5'-ends or in the middle, resulting in distinct targetomes and, consequently, functions. In the present study, we aimed to evaluate the putative targets and functional role of miR-1246 and its two 5'-isoforms (ISO-miR-1246_a and ISO-miR-1246_G) in vitro. Commercial Caco-2 cells of CRC origin were analyzed for the expression of WT-miR-1246 and its 5'-isoforms using small RNA sequencing data, and the overabundance of the two miR-1246 isoforms was determined in cells. The transcriptome analysis of Caco-2 cells transfected with WT-miR-1246, ISO-miR-1246_G, and ISO-miR-1246_a indicated the minor overlap of the targetomes between the studied miRNA isoforms. Consequently, an enrichment analysis showed the involvement of the potential targets of the miR-1246 isoforms in distinct signaling pathways. Cancer-related pathways were predominantly more enriched in dysregulated genes in ISO-miR-1246_G and ISO-miR-1246_a, whereas cell cycle pathways were more enriched in WT-miR-1246. The functional analysis of WT-miR-1246 and its two 5'-isoforms revealed that the inhibition of any of these molecules had a tumor-suppressive role (reduced cell viability and migration and promotion of early cell apoptosis) in CRC cells. However, the 5'-isoforms had a stronger effect on viability compared with WT-miR-1246. To conclude, this research shows that WT-miR-1246 and its two 5'-isoforms have different targetomes and are involved in distinct signaling pathways but collectively play an important role in CRC pathogenesis.

Circulating isomiRs May Be Superior Biomarkers Compared to Their Corresponding miRNAs: A Pilot Biomarker Study of Using isomiR-Ome to Detect Coronary Calcium-Based Cardiovascular Risk in Patients with NAFLD

Circulating miRNAs are increasingly being considered as biomarkers in various medical contexts, but the value of analyzing isomiRs (isoforms of canonical miRNA sequences) has not frequently been assessed. Here we hypothesize that an in-depth analysis of the full circulating miRNA landscape could identify specific isomiRs that are stronger biomarkers, compared to their corresponding miRNA, for identifying increased CV risk in patients with non-alcoholic fatty liver disease (NAFLD)-a clinical unmet need. Plasma miRNAs were sequenced with next-generation sequencing (NGS). Liver fat content was measured with magnetic-resonance spectrometry (MRS); CV risk was determined, beyond using traditional biomarkers, by a CT-based measurement of coronary artery calcium (CAC) score and the calculation of a CAC score-based CV-risk percentile (CAC-CV%). This pilot study included n = 13 patients, age > 45 years, with an MRS-measured liver fat content of ≥5% (wt/wt), and free of overt CVD. NGS identified 1103 miRNAs and 404,022 different isomiRs, of which 280 (25%) and 1418 (0.35%), respectively, passed an abundance threshold. Eighteen (sixteen/two) circulating miRNAs correlated positively/negatively, respectively, with CAC-CV%, nine of which also significantly discriminated between high/low CV risk through ROC-AUC analysis. IsomiR-ome analyses uncovered 67 isomiRs highly correlated (R ≥ 0.55) with CAC-CV%. Specific isomiRs of miRNAs 101-3p, 144-3p, 421, and 484 exhibited stronger associations with CAC-CV% compared to their corresponding miRNA. Additionally, while miRNAs 140-3p, 223-3p, 30e-5p, and 342-3p did not correlate with CAC-CV%, specific isomiRs with altered seed sequences exhibited a strong correlation with coronary atherosclerosis burden. Their predicted isomiRs-specific targets were uniquely enriched (compared to their canonical miRNA sequence) in CV Disease (CVD)-related pathways. Two of the isomiRs exhibited discriminative ROC-AUC, and another two showed a correlation with reverse cholesterol transport from cholesterol-loaded macrophages to ApoB-depleted plasma. In summary, we propose a pipeline for exploring circulating isomiR-ome as an approach to uncover novel and strong CVD biomarkers.

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/.

Skin treatment with non-thermal plasma modulates the immune system through miR-223-3p and its target genes

Non-thermal plasma, a partially ionized gas, holds significant potential for clinical applications, including wound-healing support, oral therapies, and anti-tumour treatments. While its applications showed promising outcomes, the underlying molecular mechanisms remain incompletely understood. We thus apply non-thermal plasma to mouse auricular skin and conducted non-coding RNA sequencing, as well as single-cell blood sequencing. In a time-series analysis (five timepoints spanning 2 hours), we compare the expression of microRNAs in the plasma-treated left ears to the unexposed right ears of the same mice as well as to the ears of unexposed control mice. Our findings indicate specific effects in the treated ears for a set of five miRNAs: mmu-miR-144-5p, mmu-miR-144-3p, mmu-miR-142a-5p, mmu-miR-223-3p, and mmu-miR-451a. Interestingly, mmu-miR-223-3p also exhibits an increase over time in the right non-treated ear of the exposed mice, suggesting systemic effects. Notably, this miRNA, along with mmu-miR-142a-5p and mmu-miR-144-3p, regulates genes and pathways associated with wound healing and tissue regeneration (namely ErbB, FoxO, Hippo, and PI3K-Akt signalling). This co-regulation is particularly remarkable considering the significant seed dissimilarities among the miRNAs. Finally, single-cell sequencing of PBMCs reveals the downregulation of 12 from 15 target genes in B-cells, Cd4+ and Cd8+ T-cells. Collectively, our data provide evidence for a systemic effect of non-thermal plasma.

Benchmarking of bioinformatics tools for NGS-based microRNA profiling with RT-qPCR method

MicroRNAs are vital gene expression regulators, extensively studied worldwide. The large-scale characterization of miRNAomes is possible using next-generation sequencing (NGS). This technology offers great opportunities, but these cannot be fully exploited without proper and comprehensive bioinformatics analysis. This may be achieved by the use of reliable dedicated software; however, different programs may generate divergent results, leading to additional discrepancies. Thus, the aim of this study was to compare three bioinformatic algorithms dedicated to NGS-based microRNA profiling and validate them using an alternative method, namely RT-qPCR. The comparison analysis revealed differences in the number and sets of identified miRNAs. The qPCR confirmed the expression of the investigated microRNAs. The correlation analysis of NGS and qPCR measurements showed strong and significant coefficients for a subset of the tested miRNAs, including those detected by all three algorithms. Single miRNA variants (isomiRs) showed different levels of correlation with the qPCR data. The obtained results revealed the good performance of all tested programs, despite the observed differences. Moreover, they implied that some specific miRNAs may be differentially estimated using NGS technology and the qPCR method, regardless of the used bioinformatics software. These discrepancies may stem from many factors, including the composition of the isomiR profile, their abundance, length, and investigated species. In conclusion, in this study, we shed light on the bioinformatics aspects of miRNAome profiling, elucidating its complexity and pinpointing potential features influencing validation. Thus, qPCR validation results should be open to interpretation when not fully concordant with NGS results until further, additional analyses are conducted.

Interrogating the Role of miR-125b and Its 3'isomiRs in Protection against Hypoxia

MiR-125b has therapeutic potential in the amelioration of myocardial ischemic injury. MicroRNA isomiRs, with either 5' or 3' addition or deletion of nucleotide(s), have been reported from next-generation sequencing data (NGS). However, due to technical challenges, validation and functional studies of isomiRs are few. In this study, we discovered using NGS, four 3'isomiRs of miR-125b, i.e., addition of A (adenosine), along with deletions of A, AG (guanosine) and AGU (uridine) from rat and sheep heart. These findings were validated using RT-qPCR. Comprehensive functional studies were carried out in the H9C2 hypoxia model. After miR-125b, isomiRs of Plus A, Trim A, AG and AGU mimic transfection, the H9C2 cells were subjected to hypoxic challenge. As assessed using cell viability, apoptosis, CCK-8 and LDH release, miR-125b and isomiRs were all protective against hypoxia. However, Plus A and Trim A were more effective than miR-125b, whilst Trim AG and Trim AGU had far weaker effects than miR-125b. Interestingly, both the gene regulation profile and apoptotic gene validation indicated a major overlap among miR-125b, Plus A and Trim A, whilst Trims AG and AGU revealed a different profile compared to miR-125b. Conclusions: miR-125b and its 3' isomiRs are expressed stably in the heart. miR-125b and isomiRs with addition or deletion of A might function concurrently and concordantly under specific physiological and pathophysiological conditions. In-depth understanding of isomiRs' metabolism and function will contribute to better miRNA therapeutic drug design.

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.

miR-183-5p overexpression orchestrates collective invasion in salivary adenoid cystic carcinoma through the FAT1/YAP1 signaling pathway

The invasion of cancer cells into interstitial tissues in a cohesive unit is termed collective invasion, and it is important for the invasion and metastasis of salivary adenoid cystic carcinoma (SACC). However, the underlying mechanisms regulating SACC collective invasion are still poorly understood. Here, we found that SACC tissues exhibited remarkable FAT1 downregulation and YAP1 upregulation at the invasive front, which was closely associated with collective invasion and distant metastasis. Decreasing FAT1 expression significantly activated the YAP1 signaling pathway and promoted collective invasion. Moreover, miR-183-5p was identified as the candidate regulator of FAT1 by bioinformatic analysis and an online database algorithm. A dual luciferase reporter experiment further confirmed that miR-183-5p directly targeted the FAT1 3'-UTR to reduce FAT1 expression. Increasing or decreasing miR-183-5p expression promoted or attenuated collective invasion, which was reversed by YAP1 siRNA or FAT1 siRNA, respectively. In addition, knocking down miR-183-5p reduced tumor burden and attenuated collective invasion in vivo. Together, these results suggest that the miR-183-5p/FAT1/YAP1 signaling pathway is a critical driver of SACC collective invasion, revealing a novel therapeutic target.

One locus, several functional RNAs-emerging roles of the mechanisms responsible for the sequence variability of microRNAs

With the development of modern molecular genetics, the original "one gene-one enzyme" hypothesis has been outdated. For protein coding genes, the discovery of alternative splicing and RNA editing provided the biochemical background for the RNA repertoire of a single locus, which also serves as an important pillar for the enormous protein variability of the genomes. Non-protein coding RNA genes were also revealed to produce several RNA species with distinct functions. The loci of microRNAs (miRNAs), encoding for small endogenous regulatory RNAs, were also found to produce a population of small RNAs, rather than a single defined product. This review aims to present the mechanisms contributing to the astonishing variability of miRNAs revealed by the new sequencing technologies. One important source is the careful balance of arm selection, producing sequentially different 5p- or 3p-miRNAs from the same pre-miRNA, thereby broadening the number of regulated target RNAs and the phenotypic response. In addition, the formation of 5', 3' and polymorphic isomiRs, with variable end and internal sequences also leads to a higher number of targeted sequences, and increases the regulatory output. These miRNA maturation processes, together with other known mechanisms such as RNA editing, further increase the potential outcome of this small RNA pathway. By discussing the subtle mechanisms behind the sequence diversity of miRNAs, this review intends to reveal this engaging aspect of the inherited "RNA world", how it contributes to the almost infinite molecular variability among living organisms, and how this variability can be exploited to treat human diseases.

Genetics and RNA Regulation of Uveal Melanoma

Uveal melanoma (UM) is the most common intraocular malignant tumor and the most frequent melanoma not affecting the skin. While the rate of UM occurrence is relatively low, about 50% of patients develop metastasis, primarily to the liver, with lethal outcome despite medical treatment. Notwithstanding that UM etiopathogenesis is still under investigation, a set of known mutations and chromosomal aberrations are associated with its pathogenesis and have a relevant prognostic value. The most frequently mutated genes are BAP1, EIF1AX, GNA11, GNAQ, and SF3B1, with mutually exclusive mutations occurring in GNAQ and GNA11, and almost mutually exclusive ones in BAP1 and SF3B1, and BAP1 and EIF1AX. Among chromosomal aberrations, monosomy of chromosome 3 is the most frequent, followed by gain of chromosome 8q, and full or partial loss of chromosomes 1 and 6. In addition, epigenetic mechanisms regulated by non-coding RNAs (ncRNA), namely microRNAs and long non-coding RNAs, have also been investigated. Several papers investigating the role of ncRNAs in UM have reported that their dysregulated expression affects cancer-related processes in both in vitro and in vivo models. This review will summarize current findings about genetic mutations, chromosomal aberrations, and ncRNA dysregulation establishing UM biology.

Dynamic and static circulating cancer microRNA biomarkers - a validation study

For cancers and other pathologies, early diagnosis remains the most promising path to survival. Profiling of longitudinal cohorts facilitates insights into trajectories of biomarkers. We measured microRNA expression in 240 serum samples from patients with colon, lung, and breast cancer and from cancer-free controls. Each patient provided at least two serum samples, one prior to diagnosis and one following diagnosis. The median time interval between the samples was 11.6 years. Using computational models, we evaluated the circulating profiles of 21 microRNAs. The analysis yielded two sets of biomarkers, static ones that show an absolute difference between certain cancer types and controls and dynamic ones where the level over time provided higher diagnostic information content. In the first group, miR-99a-5p stands out for all three cancer types. In the second group, miR-155-5p allows to predict lung cancers and colon cancers. Classification in samples from cancer and non-cancer patients using gradient boosted trees reached an average accuracy of 79.9%. The results suggest that individual change over time or an absolute value at one time point may predict a disease with high specificity and sensitivity.

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.

5’isomiR-183-5p|+2 elicits tumor suppressor activity in a negative feedback loop with E2F1

Background

MicroRNAs (miRNAs) and isomiRs play important roles in tumorigenesis as essential regulators of gene expression. 5’isomiRs exhibit a shifted seed sequence compared to the canonical miRNA, resulting in different target spectra and thereby extending the phenotypic impact of the respective common pre-miRNA. However, for most miRNAs, expression and function of 5’isomiRs have not been studied in detail yet. Therefore, this study aims to investigate the functions of miRNAs and their 5’isomiRs.

Methods

The expression of 5’isomiRs was assessed in The Cancer Genome Atlas (TCGA) breast cancer patient dataset. Phenotypic effects of miR-183 overexpression in triple-negative breast cancer (TNBC) cell lines were investigated in vitro and in vivo by quantifying migration, proliferation, tumor growth and metastasis. Direct targeting of E2F1 by miR-183-5p|+2 was validated with a 3’UTR luciferase assay and linked to the phenotypes of isomiR overexpression.

Results

TCGA breast cancer patient data indicated that three variants of miR-183-5p are highly expressed and upregulated, namely miR-183-5p|0, miR-183-5p|+1 and miR-183-5p|+2. However, TNBC cell lines displayed reduced proliferation and invasion upon overexpression of pre-miR-183. While invasion was reduced individually by all three isomiRs, proliferation and cell cycle progression were specifically inhibited by overexpression of miR-183-5p|+2. Proteomic analysis revealed reduced expression of E2F target genes upon overexpression of this isomiR, which could be attributed to direct targeting of E2F1, specifically by miR-183-5p|+2. Knockdown of E2F1 partially phenocopied the effect of miR-183-5p|+2 overexpression on cell proliferation and cell cycle. Gene set enrichment analysis of TCGA and METABRIC patient data indicated that the activity of E2F strongly correlated with the expression of miR-183-5p, suggesting transcriptional regulation of the miRNA by a factor of the E2F family. Indeed, in vitro, expression of miR-183-5p was regulated by E2F1. Hence, miR-183-5p|+2 directly targeting E2F1 appears to be part of a negative feedback loop potentially fine-tuning its activity.

Conclusions

This study demonstrates that 5’isomiRs originating from the same arm of the same pre-miRNA (i.e. pre-miR-183-5p) may exhibit different functions and thereby collectively contribute to the same phenotype. Here, one of three isomiRs was shown to counteract expression of the pre-miRNA by negatively regulating a transcriptional activator (i.e. E2F1). We speculate that this might be part of a regulatory mechanism to prevent uncontrolled cell proliferation, which is disabled during cancer progression.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-022-02380-8.

Updated review of advances in microRNAs and complex diseases: experimental results, databases, webservers and data fusion

MicroRNAs (miRNAs) are gene regulators involved in the pathogenesis of complex diseases such as cancers, and thus serve as potential diagnostic markers and therapeutic targets. The prerequisite for designing effective miRNA therapies is accurate discovery of miRNA-disease associations (MDAs), which has attracted substantial research interests during the last 15 years, as reflected by more than 55 000 related entries available on PubMed. Abundant experimental data gathered from the wealth of literature could effectively support the development of computational models for predicting novel associations. In 2017, Chen et al. published the first-ever comprehensive review on MDA prediction, presenting various relevant databases, 20 representative computational models, and suggestions for building more powerful ones. In the current review, as the continuation of the previous study, we revisit miRNA biogenesis, detection techniques and functions; summarize recent experimental findings related to common miRNA-associated diseases; introduce recent updates of miRNA-relevant databases and novel database releases since 2017, present mainstream webservers and new webserver releases since 2017 and finally elaborate on how fusion of diverse data sources has contributed to accurate MDA prediction.

isomiRTar: a comprehensive portal of pan-cancer 5'-isomiR targeting

Inaccurate cleavage of pri- and pre-miRNA hairpins by Drosha and Dicer results in the generation of miRNA isoforms known as isomiRs. isomiRs with 5'-end variations (5'-isomiRs) create a new dimension in miRNA research since they have different seed regions and distinct targetomes. We developed isomiRTar (https://isomirtar.hse.ru)-a comprehensive portal that allows one to analyze expression profiles and targeting activity of 5'-isomiRs in cancer. Using the Cancer Genome Atlas sequencing data, we compiled the list of 1022 5'-isomiRs expressed in 9282 tumor samples across 31 cancer types. Sequences of these isomiRs were used to predict target genes with miRDB and TargetScan. The putative interactions were then subjected to the co-expression analysis in each cancer type to identify isomiR-target pairs supported by significant negative correlations. Downstream analysis of the data deposited in isomiRTar revealed both cancer-specific and cancer-conserved 5'-isomiR expression landscapes. Pairs of isomiRs differing in one nucleotide shift from 5'-end had poorly overlapping targetomes with the median Jaccard index of 0.06. The analysis of colorectal cancer 5'-isomiR-mediated regulatory networks revealed promising candidate tumor suppressor isomiRs: hsa-miR-203a-3p-+1, hsa-miR-192-5p-+1 and hsa-miR-148a-3p-0. In summary, we believe that isomiRTar will help researchers find novel mechanisms of isomiR-mediated gene silencing in different types of cancer.

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.

The placenta epigenome-brain axis: placental epigenomic and transcriptomic responses that preprogram cognitive impairment

Aim: The placenta-brain axis reflects a developmental linkage where disrupted placental function is associated with impaired neurodevelopment later in life. Placental gene expression and the expression of epigenetic modifiers such as miRNAs may be tied to these impairments and are understudied. Materials & methods: The expression levels of mRNAs (n = 37,268) and their targeting miRNAs (n = 2083) were assessed within placentas collected from the ELGAN study cohort (n = 386). The ELGAN adolescents were assessed for neurocognitive function at age 10 and the association with placental mRNA/miRNAs was determined. Results: Placental mRNAs related to inflammatory and apoptotic processes are under miRNA control and associated with cognitive impairment at age 10. Conclusion: Findings highlight key placenta epigenome-brain relationships that support the developmental origins of health and disease hypothesis.

From Molecular Mechanisms to Therapeutics: Understanding MicroRNA-21 in Cancer

MicroRNAs (miRNAs) are small noncoding RNAs that play an important role in regulating gene expression at a posttranscriptional level. As one of the first discovered oncogenic miRNAs, microRNA-21 (miR-21) has been highlighted for its critical role in cancers, such as glioblastoma, pancreatic adenocarcinoma, non-small cell lung cancer, and many others. MiR-21 targets many vital components in a wide range of cancers and acts on various cellular processes ranging from cancer stemness to cell death. Expression of miR-21 is elevated within cancer tissues and circulating miR-21 is readily detectable in biofluids, making it valuable as a cancer biomarker with significant potential for use in diagnosis and prognosis. Advances in RNA-based therapeutics have revealed additional avenues by which miR-21 can be utilized as a promising target in cancer. The purpose of this review is to outline the roles of miR-21 as a key modulator in various cancers and its potential as a therapeutic target.

A deep learning method for miRNA/isomiR target detection

Accurate identification of microRNA (miRNA) targets at base-pair resolution has been an open problem for over a decade. The recent discovery of miRNA isoforms (isomiRs) adds more complexity to this problem. Despite the existence of many methods, none considers isomiRs, and their performance is still suboptimal. We hypothesize that by taking the isomiR–mRNA interactions into account and applying a deep learning model to study miRNA–mRNA interaction features, we may improve the accuracy of miRNA target predictions. We developed a deep learning tool called DMISO to capture the intricate features of miRNA/isomiR–mRNA interactions. Based on tenfold cross-validation, DMISO showed high precision (95%) and recall (90%). Evaluated on three independent datasets, DMISO had superior performance to five tools, including three popular conventional tools and two recently developed deep learning-based tools. By applying two popular feature interpretation strategies, we demonstrated the importance of the miRNA regions other than their seeds and the potential contribution of the RNA-binding motifs within miRNAs/isomiRs and mRNAs to the miRNA/isomiR–mRNA interactions.

EpisomiR, a New Family of miRNAs, and Its Possible Roles in Human Diseases

MicroRNAs (miRNAs) are synthesized through a canonical pathway and play a role in human diseases, such as cancers and cardiovascular, neurodegenerative, psychiatric, and chronic inflammatory diseases. The development of sequencing technologies has enabled the identification of variations in noncoding miRNAs. These miRNA variants, called isomiRs, are generated through a non-canonical pathway, by several enzymes that alter the length and sequence of miRNAs. The isomiR family is, now, expanding further to include episomiRs, which are miRNAs with different modifications. Since recent findings have shown that isomiRs reflect the cell-specific biological function of miRNAs, knowledge about episomiRs and isomiRs can, possibly, contribute to the optimization of diagnosis and therapeutic technology for precision medicine.

Isoforms of miR-148a and miR-203a are putative suppressors of colorectal cancer

MicroRNAs are short non-coding molecules which regulate translation in a gene-specific manner. MicroRNA isoforms that differ by few extra or missing nucleotides at the 5'-terminus (5'-isomiR) show strikingly different target specificity. This study aimed to identify functional roles of 5′-isomiR in colorectal cancers. Transcriptomic targets of microRNA isoforms were predicted using bioinformatics tools miRDB and TargetScan. The sets of putative targets identified for 5′-isomiR were integrated with mRNA and microRNA sequencing data for primary colorectal tumors retrieved from The Cancer Genome Atlas Colon Adenocarcinoma (TCGA-COAD) database. The network of interactions among miRNA, their targets and transcription factors was built using the miRGTF-net algorithm. The results indicate that microRNA isoforms highly expressed in colorectal cancer and differing by a single nucleotide position at the 5'-terminus have ≤ 30% common targets. The regulatory network of interactions enables identification of the most engaged microRNA isoforms. Anti-correlated expression levels of canonical microRNA hsa-miR-148a-3p and its putative targets including CSF1, ETS1, FLT1, ITGA5, MEIS1, MITF and RUNX2 proliferation regulators suggest an anti-tumor role for this molecule. The canonical microRNA hsa-miR-203a-3p|0 and its 5′-isoform bind different sets of anti-correlated putative targets, although both of them interact with genes involved in the epithelial-mesenchymal transition: SNAI2 and TNC.

Small RNA-Sequencing for Analysis of Circulating miRNAs: Benchmark Study

Small RNA-sequencing (RNA-Seq) is being increasingly used for profiling of circulating microRNAs (miRNAs), a new group of promising biomarkers. Unfortunately, small RNA-Seq protocols are prone to biases limiting quantification accuracy, which motivated development of several novel methods. Here, we present comparison of all small RNA-Seq library preparation approaches that are commercially available for quantification of miRNAs in biofluids. Using synthetic and human plasma samples, we compared performance of traditional two-adaptor ligation protocols (Lexogen, Norgen), as well as methods using randomized adaptors (NEXTflex), polyadenylation (SMARTer), circularization (RealSeq), capture probes (EdgeSeq), or unique molecular identifiers (QIAseq). There was no single protocol outperforming others across all metrics. Limited overlap of measured miRNA profiles was documented between methods largely owing to protocol-specific biases. Methods designed to minimize bias largely differ in their performance, and contributing factors were identified. Usage of unique molecular identifiers has rather negligible effect and, if designed incorrectly, can even introduce spurious results. Together, these results identify strengths and weaknesses of all current methods and provide guidelines for applications of small RNA-Seq in biomarker research.

ExhauFS: exhaustive search-based feature selection for classification and survival regression

Feature selection is one of the main techniques used to prevent overfitting in machine learning applications. The most straightforward approach for feature selection is an exhaustive search: one can go over all possible feature combinations and pick up the model with the highest accuracy. This method together with its optimizations were actively used in biomedical research, however, publicly available implementation is missing. We present ExhauFS-the user-friendly command-line implementation of the exhaustive search approach for classification and survival regression. Aside from tool description, we included three application examples in the manuscript to comprehensively review the implemented functionality. First, we executed ExhauFS on a toy cervical cancer dataset to illustrate basic concepts. Then, multi-cohort microarray breast cancer datasets were used to construct gene signatures for 5-year recurrence classification. The vast majority of signatures constructed by ExhauFS passed 0.65 threshold of sensitivity and specificity on all datasets, including the validation one. Moreover, a number of gene signatures demonstrated reliable performance on independent RNA-seq dataset without any coefficient re-tuning, i.e., turned out to be cross-platform. Finally, Cox survival regression models were used to fit isomiR signatures for overall survival prediction for patients with colorectal cancer. Similarly to the previous example, the major part of models passed the pre-defined concordance index threshold 0.65 on all datasets. In both real-world scenarios (breast and colorectal cancer datasets), ExhauFS was benchmarked against state-of-the-art feature selection models, including L₁-regularized sparse models. In case of breast cancer, we were unable to construct reliable cross-platform classifiers using alternative feature selection approaches. In case of colorectal cancer not a single model passed the same 0.65 threshold. Source codes and documentation of ExhauFS are available on GitHub: https://github.com/s-a-nersisyan/ExhauFS.

Non-canonical features of microRNAs: paradigms emerging from cardiovascular disease

Research showing that microRNAs (miRNAs) are versatile regulators of gene expression has instigated tremendous interest in cardiovascular research. The overwhelming majority of studies are predicated on the dogmatic notion that miRNAs regulate the expression of specific target mRNAs by inhibiting mRNA translation or promoting mRNA decay in the RNA-induced silencing complex (RISC). These efforts mostly identified and dissected contributions of multiple regulatory networks of miRNA-target mRNAs to cardiovascular pathogenesis. However, evidence from studies in the past decade indicates that miRNAs also operate beyond this canonical paradigm, featuring non-conventional regulatory functions and cellular localizations that have a pathophysiological role in cardiovascular disease. In this Review, we highlight the functional relevance of atypical miRNA biogenesis and localization as well as RISC heterogeneity. Moreover, we delineate remarkable non-canonical examples of miRNA functionality, including direct interactions with proteins beyond the Argonaute family and their role in transcriptional regulation in the nucleus and in mitochondria. We scrutinize the relevance of non-conventional biogenesis and non-canonical functions of miRNAs in cardiovascular homeostasis and pathology, and contextualize how uncovering these non-conventional properties can expand the scope of translational research in the cardiovascular field and beyond.

A Novel Tissue Atlas and Online Tool for the Interrogation of Small RNA Expression in Human Tissues and Biofluids

One promising goal for utilizing the molecular information circulating in biofluids is the discovery of clinically useful biomarkers. Extracellular RNAs (exRNAs) are one of the most diverse classes of molecular cargo, easily assayed by sequencing and with expressions that rapidly change in response to subject status. Despite diverse exRNA cargo, most evaluations from biofluids have focused on small RNA sequencing and analysis, specifically on microRNAs (miRNAs). Another goal of characterizing circulating molecular information, is to correlate expression to injuries associated with specific tissues of origin. Biomarker candidates are often described as being specific, enriched in a particular tissue or associated with a disease process. Likewise, miRNA data is often reported to be specific, enriched for a tissue, without rigorous testing to support the claim. Here we provide a tissue atlas of small RNAs from 30 different tissues and three different blood cell types. We analyzed the tissues for enrichment of small RNA sequences and assessed their expression in biofluids: plasma, cerebrospinal fluid, urine, and saliva. We employed published data sets representing physiological (resting vs. acute exercise) and pathologic states (early- vs. late-stage liver fibrosis, and differential subtypes of stroke) to determine differential tissue-enriched small RNAs. We also developed an online tool that provides information about exRNA sequences found in different biofluids and tissues. The data can be used to better understand the various types of small RNA sequences in different tissues as well as their potential release into biofluids, which should help in the validation or design of biomarker studies.

Non-coding RNAs as liquid biopsy biomarkers in cancer

Although non-coding RNAs have long been considered as non-functional "junk" RNAs, accumulating evidence in the past decade indicates that they play a critical role in pathogenesis of various cancers. In addition to their biological significance, the recognition that their expression levels are frequently dysregulated in multiple cancers have fueled the interest for exploiting their clinical potential as cancer biomarkers. In particular, microRNAs (miRNAs), a subclass of small non-coding RNAs that epigenetically modulate gene-transcription, have become one of the most well-studied substrates for the development of liquid biopsy biomarkers for cancer patients. The emergence of high-throughput sequencing technologies has enabled comprehensive molecular characterisation of various non-coding RNA expression profiles in multiple cancers. Furthermore, technological advances for quantifying lowly expressed RNAs in the circulation have facilitated robust identification of previously unrecognised and undetectable biomarkers in cancer patients. Here we summarise the latest progress on the utilisation of non-coding RNAs as non-invasive cancer biomarkers. We evaluated the suitability of multiple non-coding RNA types as blood-based cancer biomarkers and examined the impact of recent technological breakthroughs on the development of non-invasive molecular biomarkers in cancer.

Experimental MicroRNA Detection Methods

MicroRNAs (miRNAs) are considerably small yet highly important riboregulators involved in nearly all cellular processes. Due to their critical roles in posttranscriptional regulation of gene expression, they have the potential to be used as biomarkers in addition to their use as drug targets. Although computational approaches speed up the initial genomewide identification of putative miRNAs, experimental approaches are essential for further validation and functional analyses of differentially expressed miRNAs. Therefore, sensitive, specific, and cost-effective microRNA detection methods are imperative for both individual and multiplex analysis of miRNA expression in different tissues and during different developmental stages. There are a number of well-established miRNA detection methods that can be exploited depending on the comprehensiveness of the study (individual miRNA versus multiplex analysis), the availability of the sample and the location and intracellular concentration of miRNAs. This review aims to highlight not only traditional but also novel strategies that are widely used in experimental identification and quantification of microRNAs.

Coupling miR/isomiR and mRNA Expression Signatures Unveils New Molecular Layers of Endometrial Receptivity

Embryo implantation depends on endometrial receptivity (ER). To achieve ER, the preparation of the uterine lining requires controlled priming by ovarian hormones and the expression of numerous genes in the endometrial tissue. microRNAs (miRs) have emerged as critical genetic regulators of ER in fertility and of the diseases that are associated with infertility. With the rapid development of next-generation sequencing technologies, it has become clear that miR genes can produce canonical miRs and variants—isomiRs. Here, we describe miR/isomiR expression dynamics across the four time points of natural chorionic gonadotropin (hCG)-administered cycles. Sequencing of the small RNAs (sRNA-seq) revealed that the most significant expression changes during the transition from the pre-receptive to the receptive phase occurred in the isomiR families of miR-125a, miR-125b, miR-10a, miR-10b, miR-449c, miR-92a, miR-92b, and miR-99a. Pairing the analysis of the differentially expressed (DE) miRs/isomiRs and their predicted DE mRNA targets uncovered 280 negatively correlating pairs. In the receptive endometrium, the 5′3′-isomiRs of miR-449c, which were among the most highly up-regulated isomiRs, showed a negative correlation with their target, transcription factor (TF) MYCN, which was down-regulated. Joint analysis of the miR/isomiR and TF expression identified several regulatory interactions. Based on these data, a regulatory TF-miR/isomiR gene-target circuit including let7g-5p and miR-345; the isomiR families of miR-10a, miR-10b, miR-92a, and miR-449c; and MYCN and TWIST1 was proposed to play a key role in the establishment of ER. Our work uncovers the complexity and dynamics of the endometrial isomiRs that can act cooperatively with miRs to control the functionally important genes that are critical to ER. Further studies of miR/isomiR expression patterns that are paired with those of their target mRNAs may provide a more in-depth picture of the endometrial pathologies that are associated with implantation failure.

ECM-Receptor Regulatory Network and Its Prognostic Role in Colorectal Cancer

Interactions of the extracellular matrix (ECM) and cellular receptors constitute one of the crucial pathways involved in colorectal cancer progression and metastasis. With the use of bioinformatics analysis, we comprehensively evaluated the prognostic information concentrated in the genes from this pathway. First, we constructed a ECM-receptor regulatory network by integrating the transcription factor (TF) and 5'-isomiR interaction databases with mRNA/miRNA-seq data from The Cancer Genome Atlas Colon Adenocarcinoma (TCGA-COAD). Notably, one-third of interactions mediated by 5'-isomiRs was represented by noncanonical isomiRs (isomiRs, whose 5'-end sequence did not match with the canonical miRBase version). Then, exhaustive search-based feature selection was used to fit prognostic signatures composed of nodes from the network for overall survival prediction. Two reliable prognostic signatures were identified and validated on the independent The Cancer Genome Atlas Rectum Adenocarcinoma (TCGA-READ) cohort. The first signature was made up by six genes, directly involved in ECM-receptor interaction: AGRN, DAG1, FN1, ITGA5, THBS3, and TNC (concordance index 0.61, logrank test p = 0.0164, 3-years ROC AUC = 0.68). The second hybrid signature was composed of three regulators: hsa-miR-32-5p, NR1H2, and SNAI1 (concordance index 0.64, logrank test p = 0.0229, 3-years ROC AUC = 0.71). While hsa-miR-32-5p exclusively regulated ECM-related genes (COL1A2 and ITGA5), NR1H2 and SNAI1 also targeted other pathways (adhesion, cell cycle, and cell division). Concordant distributions of the respective risk scores across four stages of colorectal cancer and adjacent normal mucosa additionally confirmed reliability of the models.

Sperm-borne sncRNAs: potential biomarkers for semen fertility?

Semen infertility or sub-fertility, whether in humans or livestock species, remains a major concern for clinicians and technicians involved in reproduction. Indeed, they can cause tragedies in human relationships or have a dramatic overall negative impact on the sustainability of livestock breeding. Understanding and predicting semen fertility issues is therefore crucial and quality control procedures as well as biomarkers have been proposed to ensure sperm fertility. However, their predictive values appeared to be too limited and additional relevant biomarkers are still required to diagnose sub-fertility efficiently. During the last decade, the study of molecular mechanisms involved in spermatogenesis and sperm maturation highlighted the regulatory role of a variety of small non-coding RNAs (sncRNAs) and led to the discovery that sperm sncRNAs comprise both remnants from spermatogenesis and post-testicular sncRNAs acquired through interactions with extracellular vesicles along epididymis. This has led to the hypothesis that sncRNAs may be a source of relevant biomarkers, associated either with sperm functionality or embryo development. This review aims at providing a synthetic overview of the current state of knowledge regarding implication of sncRNA in spermatogenesis defects and their putative roles in sperm maturation and embryo development, as well as exploring their use as fertility biomarkers.

Hsa-miR-21-3p associates with breast cancer patient survival and targets genes in tumor suppressive pathways

Breast cancer is the cancer most often diagnosed in women. MicroRNAs (MIRs) are short RNA molecules that bind mRNA resulting in their downregulation. MIR21 has been shown to be an oncomiR in most cancer types, including breast cancer. Most of the effects of miR-21 have been attributed to hsa-miR-21-5p that is transcribed from the leading strand of MIR21, but hsa-miR-21-3p (miR-21-3p), transcribed from the lagging strand, is much less studied. The aim of the study is to analyze whether expression of miR-21-3p is prognostic for breast cancer. MiR-21-3p association with survival, clinical and pathological characteristics was analyzed in a large breast cancer cohort and validated in three separate cohorts, including TCGA and METABRIC. Analytical tools were also used to infer miR-21-3p function and to identify potential target genes and functional pathways. The results showed that in the exploration cohort, high miR-21-3p levels associated with shorter survival and lymph node positivity. In the three validation cohorts, high miR-21-3p levels associated with pathological characteristics that predict worse prognosis. Specifically, in the largest validation cohort, METABRIC (n = 1174), high miR-21-3p levels associated with large tumors, a high grade, lymph node and HER2 positivity, and shorter breast-cancer-specific survival (HR = 1.38, CI 1.13–1.68). This association remained significant after adjusting for confounding factors. The genes with expression levels that correlated with miR-21-3p were enriched in particular pathways, including the epithelial-to-mesenchymal transition and proliferation. Among the most significantly downregulated targets were MAT2A and the tumor suppressive genes STARD13 and ZNF132. The results from this study emphasize that both 3p- and 5p-arms from a MIR warrant independent study. The data show that miR-21-3p overexpression in breast tumors is a marker of worse breast cancer progression and it affects genes in pathways that drive breast cancer by down-regulating tumor suppressor genes. The results suggest miR-21-3p as a potential biomarker.

Targeted designing functional markers revealed the role of retrotransposon derived miRNAs as mobile epigenetic regulators in adaptation responses of pistachio

We developed novel miRNA-based markers based on salt responsive miRNA sequences to detect polymorphisms in miRNA sequences and locations. The validation of 76 combined miRNA + miRNA and miRNA + ISSR markers in the three extreme pistachio populations led to the identification of three selected markers that could link salt tolerance phenotype to genotype and divided pistachio genotypes and Pistacia species into three clusters. This novel functional marker system, in addition to more efficient performance, has higher polymorphisms than previous miRNA-based marker systems. The functional importance of the target gene of five miRNAs in the structure of the three selected markers in regulation of different genes such as ECA2, ALA10, PFK, PHT1;4, PTR3, KUP2, GRAS, TCP, bHLH, PHD finger, PLATZ and genes involved in developmental, signaling and biosynthetic processes shows that the polymorphism associated with these selected miRNAs can make a significant phenotypic difference between salt sensitive and tolerant pistachio genotypes. The sequencing results of selected bands showed the presence of conserved miRNAs in the structure of the mitochondrial genome. Further notable findings of this study are that the sequences of PCR products of two selected markers were annotated as Gypsy and Copia retrotransposable elements. The transposition of retrotransposons with related miRNAs by increasing the number of miRNA copies and changing their location between nuclear and organellar genomes can affect the regulatory activity of these molecules. These findings show the crucial role of retrotransposon-derived miRNAs as mobile epigenetic regulators between intracellular genomes in regulating salt stress responses as well as creating new and tolerant phenotypes for adaptation to environmental conditions.

Identification and Elucidation of the Protective isomiRs in Lung Cancer Patient Prognosis

MicroRNAs (miRNAs) are approximately 20–22 nucleotides in length, which are well known to participate in the post-transcriptional modification. The mature miRNAs were observed to be varied on 5′ or 3′ that raise another term—the isoforms of mature miRNAs (isomiRs), which have been proven not the artifacts and discussed widely recently. In our research, we focused on studying the 5′ isomiRs in lung adenocarcinoma (LUAD) in The Cancer Genome Atlas (TCGA). We characterized 75 isomiRs significantly associated with better prognosis and 43 isomiRs with poor prognosis. The 75 protective isomiRs can successfully distinguish tumors from normal samples and are expressed differently between patients of early and late stages. We also found that most of the protective isomiRs tend to be with downstream shift and upregulated compared with those with upstream shift, implying that a possible selection occurs during cancer development. Among these protective isomiRs, we observed a highly positive and significant correlation, as well as in harmful isomiRs, suggesting cooperation within the group. However, between protective and harmful, there is no such a concordance but conversely more negative correlation, suggesting the possible antagonistic effect between protective and harmful isomiRs. We also identified that two isomiRs miR-181a-3p|-3 and miR-181a-3p|2, respectively, belong to the harmful and protective groups, suggesting a bidirectional regulation of their originated archetype—miR-181a-3p. Additionally, we found that the protective isomiRs of miR-21-5p, which is an oncomiR, may be evolved as the tumor suppressors through producing isomiRs to hinder metastasis. In summary, these results displayed the characteristics of the protective isomiRs and their potential for developing the treatment of lung cancer.

Emerging Role of isomiRs in Cancer: State of the Art and Recent Advances

The advent of Next Generation Sequencing technologies brought with it the discovery of several microRNA (miRNA) variants of heterogeneous lengths and/or sequences. Initially ascribed to sequencing errors/artifacts, these isoforms, named isomiRs, are now considered non-canonical variants that originate from physiological processes affecting the canonical miRNA biogenesis. To date, accurate IsomiRs abundance, biological activity, and functions are not completely understood; however, the study of isomiR biology is an area of great interest due to their high frequency in the human miRNome, their putative functions in cooperating with the canonical miRNAs, and potential for exhibiting novel functional roles. The discovery of isomiRs highlighted the complexity of the small RNA transcriptional landscape in several diseases, including cancer. In this field, the study of isomiRs could provide further insights into the miRNA biology and its implication in oncogenesis, possibly providing putative new cancer diagnostic, prognostic, and predictive biomarkers as well. In this review, a comprehensive overview of the state of research on isomiRs in different cancer types, including the most common tumors such as breast cancer, colorectal cancer, melanoma, and prostate cancer, as well as in the less frequent tumors, as for example brain tumors and hematological malignancies, will be summarized and discussed.

miR-21 and miR-146a: The microRNAs of inflammaging and age-related diseases

The first paper on "inflammaging" published in 2001 paved the way for a unifying theory on how and why aging turns out to be the main risk factor for the development of the most common age-related diseases (ARDs). The most exciting challenge on this topic was explaining how systemic inflammation steeps up with age and why it shows different rates among individuals of the same chronological age. The "epigenetic revolution" in the past twenty years conveyed that the assessment of the individual genetic make-up is not enough to depict the trajectories of age-related inflammation. Accordingly, others and we have been focusing on the role of non-coding RNA, i.e. microRNAs (miRNAs), in inflammaging. The results obtained in the latest 10 years underpinned the key role of a miRNA subset that we have called inflammamiRs, owing to their ability to master (NF-κB)-driven inflammatory pathways. In this review, we will focus on two inflammamiRs, i.e. miR-21-5p and miR-146a-5p, which target a variety of molecules belonging to the NF-κB/NLRP3 pathways. The interplay between miR-146a-5p and IL-6 in the context of aging and ARDs will also be highlighted. We will also provide the most relevant evidence suggesting that circulating inflammamiRs, along with IL-6, can measure the degree of inflammaging.

MicroRNA Interactome Multiomics Characterization for Cancer Research and Personalized Medicine: An Expert Review

MicroRNAs (miRNAs) that are mutually modulated by their interacting partners (interactome) are being increasingly noted for their significant role in pathogenesis and treatment of various human cancers. Recently, miRNA interactome dissected with multiomics approaches has been the subject of focus since individual tools or methods failed to provide the necessary comprehensive clues on the complete interactome. Even though single-omics technologies such as proteomics can uncover part of the interactome, the biological and clinical understanding still remain incomplete. In this study, we present an expert review of studies involving multiomics approaches to identification of miRNA interactome and its application in mechanistic characterization, classification, and therapeutic target identification in a variety of cancers, and with a focus on proteomics. We also discuss individual or multiple miRNA-based interactome identification in various pathological conditions of relevance to clinical medicine. Various new single-omics methods that can be integrated into multiomics cancer research and the computational approaches to analyze and predict miRNA interactome are also highlighted in this review. In all, we contextulize the power of multiomics approaches and the importance of the miRNA interactome to achieve the vision and practice of predictive, preventive, and personalized medicine in cancer research and clinical oncology.

Hairpin sequence and structure is associated with features of isomiR biogenesis

MiRNA isoforms (isomiRs) are single stranded small RNAs originating from the same pri-miRNA hairpin as a result of cleavage by Drosha and Dicer enzymes. Variations at the 5'-end of a miRNA alter the seed region of the molecule, thus affecting the targetome of the miRNA. In this manuscript, we analysed the distribution of miRNA cleavage positions across 31 different cancers using miRNA sequencing data of TCGA project. As a result, we found that the processing positions are not tissue specific and that all miRNAs could be correctly classified as ones exhibiting homogeneous or heterogeneous cleavage at one of the four cleavage sites. In 42% of cases (42 out of 100 miRNAs), we observed imprecise 5'-end Dicer cleavage, while this fraction was only 14% for Drosha (14 out of 99). To the contrary, almost all cleavage sites of 3'-ends (either Drosha or Dicer) were heterogeneous. With the use of only four nucleotides surrounding a 5'-end Dicer cleavage position we built a model which allowed us to distinguish between homogeneous and heterogeneous cleavage with the reliable quality (ROC AUC = 0.68). Finally, we showed the possible applications of the study by the analysis of two 5'-end isoforms originating from the same exogeneous shRNA hairpin. It turned out that the less expressed shRNA variant was functionally active, which led to the increased off-targeting. Thus, the obtained results could be applied to the design of shRNAs whose processing will result in a single 5'-variant.

Small RNA-Sequencing: Approaches and Considerations for miRNA Analysis

MicroRNAs (miRNAs) are a class of small RNA molecules that have an important regulatory role in multiple physiological and pathological processes. Their disease-specific profiles and presence in biofluids are properties that enable miRNAs to be employed as non-invasive biomarkers. In the past decades, several methods have been developed for miRNA analysis, including small RNA sequencing (RNA-seq). Small RNA-seq enables genome-wide profiling and analysis of known, as well as novel, miRNA variants. Moreover, its high sensitivity allows for profiling of low input samples such as liquid biopsies, which have now found applications in diagnostics and prognostics. Still, due to technical bias and the limited ability to capture the true miRNA representation, its potential remains unfulfilled. The introduction of many new small RNA-seq approaches that tried to minimize this bias, has led to the existence of the many small RNA-seq protocols seen today. Here, we review all current approaches to cDNA library construction used during the small RNA-seq workflow, with particular focus on their implementation in commercially available protocols. We provide an overview of each protocol and discuss their applicability. We also review recent benchmarking studies comparing each protocol's performance and summarize the major conclusions that can be gathered from their usage. The result documents variable performance of the protocols and highlights their different applications in miRNA research. Taken together, our review provides a comprehensive overview of all the current small RNA-seq approaches, summarizes their strengths and weaknesses, and provides guidelines for their applications in miRNA research.

Novel, non-invasive markers for detecting therapy induced neuroendocrine differentiation in castration-resistant prostate cancer patients

Neuroendocrine prostate cancer (NEPC), a highly aggressive variant of castration-resistant prostate cancer (CRPC), often emerges upon treatment with androgen pathway inhibitors, via neuroendocrine differentiation. Currently, NEPC diagnosis is challenging as available markers are not sufficiently specific. Our objective was to identify novel, extracellular vesicles (EV)-based biomarkers for diagnosing NEPC. Towards this, we performed small RNA next generation sequencing in serum EVs isolated from a cohort of CRPC patients with adenocarcinoma characteristics (CRPC-Adeno) vs CRPC-NE and identified significant dysregulation of 182 known and 4 novel miRNAs. We employed machine learning algorithms to develop an 'EV-miRNA classifier' that could robustly stratify 'CRPC-NE' from 'CRPC-Adeno'. Examination of protein repertoire of exosomes from NEPC cellular models by mass spectrometry identified thrombospondin 1 (TSP1) as a specific biomarker. In view of our results, we propose that a miRNA panel and TSP1 can be used as novel, non-invasive tools to identify NEPC and guide treatment decisions. In conclusion, our study identifies for the first time, novel non-invasive exosomal/extracellular vesicle based biomarkers for detecting neuroendocrine differentiation in advanced castration resistant prostate cancer patients with important translational implications in clinical management of these patients that is currently extremely challenging.

Tumor IsomiR Encyclopedia (TIE): a pan-cancer database of miRNA isoforms

SUMMARY

MicroRNAs (miRNAs) are master regulators of gene expression in cancers. Their sequence variants or isoforms (isomiRs) are highly abundant and possess unique functions. Given their short sequence length and high heterogeneity, mapping isomiRs can be challenging; without adequate depth and data aggregation, low frequency events are often disregarded. To address these challenges, we present the Tumor IsomiR Encyclopedia (TIE): a dynamic database of isomiRs from over 10 000 adult and pediatric tumor samples in The Cancer Genome Atlas (TCGA) and The Therapeutically Applicable Research to Generate Effective Treatments (TARGET) projects. A key novelty of TIE is its ability to annotate heterogeneous isomiR sequences and aggregate the variants obtained across all datasets. Results can be browsed online or downloaded as spreadsheets. Here, we show analysis of isomiRs of miR-21 and miR-30a to demonstrate the utility of TIE.

AVAILABILITY AND IMPLEMENTATION

TIE search engine and data are freely available to use at https://isomir.ccr.cancer.gov/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

On the impact of batch effect correction in TCGA isomiR expression data

MicroRNAs (miRNAs) are small non-coding RNAs with diverse functions in post-transcriptional regulation of gene expression. Sequence and length variants of miRNAs are called isomiRs and can exert different functions compared to their canonical counterparts. The Cancer Genome Atlas (TCGA) provides isomiR-level expression data for patients of various cancer entities collected in a multi-center approach over several years. However, the impact of batch effects within individual cohorts has not been systematically investigated and corrected for before. Therefore, the aim of this study was to identify relevant cohort-specific batch variables and generate batch-corrected isomiR expression data for 16 TCGA cohorts. The main batch variables included sequencing platform, plate, sample purity and sequencing depth. Platform bias was related to certain length and sequence features of individual recurrently affected isomiRs. Furthermore, significant downregulation of reported tumor suppressive isomiRs in lung tumor tissue compared to normal samples was only observed after batch correction, highlighting the importance of working with corrected data. Batch-corrected datasets for all cohorts including quality control are provided as supplement. In summary, this study reveals that batch effects present in the TCGA dataset might mask biologically relevant effects and provides a valuable resource for research on isomiRs in cancer (accessible through GEO: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE164767).

The Role of miRNAs, miRNA Clusters, and isomiRs in Development of Cancer Stem Cell Populations in Colorectal Cancer

MicroRNAs (miRNAs or miRs) have a critical role in regulating stem cells (SCs) during development and altered expression can cause developmental defects and/or disease. Indeed, aberrant miRNA expression leads to wide-spread transcriptional dysregulation which has been linked to many cancers. Mounting evidence also indicates a role for miRNAs in the development of the cancer SC (CSC) phenotype. Our goal herein is to provide a review of: (i) current research on miRNAs and their targets in colorectal cancer (CRC), and (ii) miRNAs that are differentially expressed in colon CSCs. MicroRNAs can work in clusters or alone when targeting different SC genes to influence CSC phenotype. Accordingly, we discuss the specific miRNA cluster classifications and isomiRs that are predicted to target the ALDH1, CD166, BMI1, LRIG1, and LGR5 SC genes. miR-23b and miR-92A are of particular interest because our previously reported studies on miRNA expression in isolated normal versus malignant human colonic SCs showed that miR-23b and miR-92a are regulators of the LGR5 and LRIG1 SC genes, respectively. We also identify additional miRNAs whose expression inversely correlated with mRNA levels of their target genes and associated with CRC patient survival. Altogether, our deliberation on miRNAs, their clusters, and isomiRs in regulation of SC genes could provide insight into how dysregulation of miRNAs leads to the emergence of different CSC populations and SC overpopulation in CRC.

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

MOTIVATION

MicroRNA (miRNA) precursor arms give rise to multiple isoforms simultaneously called "isomiRs." IsomiRs from the same arm typically differ by a few nucleotides at either their 5´ or 3´ termini, or both. In humans, the identities and abundances of isomiRs depend on a person's sex, population of origin, race/ethnicity, and on tissue type, tissue state, and disease type/subtype. Moreover, nearly half of the time the most abundant isomiR differs from the miRNA sequence found in public databases. Accurate mining of isomiRs from deep sequencing data is thus important.

RESULTS

We developed isoMiRmap, a fast, standalone, user-friendly mining tool that identifies and quantifies all isomiRs by directly processing short RNA-seq datasets. IsoMiRmap is a portable "plug-and-play" tool, requires minimal setup, has modest computing and storage requirements, and can process an RNA-seq dataset with 50 million reads in just a few minutes on an average laptop. IsoMiRmap deterministically and exhaustively reports all isomiRs in a given deep sequencing dataset and quantifies them accurately (no double-counting). IsoMiRmap comprehensively reports all miRNA precursor locations from which an isomiR may be transcribed, tags as 'ambiguous' isomiRs whose sequences exist both inside and outside of the space of known miRNA sequences and reports the public identifiers of common single-nucleotide polymorphisms and documented somatic mutations that may be present in an isomiR. IsoMiRmap also identifies isomiRs with 3´ non-templated post-transcriptional additions. Compared to similar tools, isoMiRmap is the fastest, reports more bona fide isomiRs, and provides the most comprehensive information related to an isomiR's transcriptional origin.

AVAILABILITY

The codes for isoMiRmap are freely available at https://cm.jefferson.edu/isoMiRmap/ and https://github.com/TJU-CMC-Org/isoMiRmap/. IsomiR profiles for the datasets of the 1000 Genomes Project, spanning five population groups, and The Cancer Genome Atlas (TCGA), spanning 33 cancer studies, are also available at https://cm.jefferson.edu/isoMiRmap/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Rola miRNA w rozwoju wybranych nowotworów – potencjalne zastosowanie w diagnostyce*

MikroRNA (miRNA) są małymi cząsteczkami kwasu rybonukleinowego, które mimo że nie podlegają procesowi translacji, pełnią ważną funkcję regulacyjną w komórkach eukariotycznych. Ich fizjologiczną funkcją jest utrzymywanie homeostazy komórek. Zaburzona ekspresja miRNA może spowodować rozwój wielu chorób, w tym chorób nowotworowych. Działanie miRNA polega na hamowaniu tworzenia się białek, w tym białek o właściwościach onkogennych i antyonkogennych. Mutacje w miejscach kodowania miRNA mogą prowadzić do nadmiernego lub zmniejszonego wytwarzania wspomnianych białek. Odkrycie miRNA i poznanie ich roli w komórce otworzyło nowe możliwości dla diagnostyki chorób nowotworowych. Zmiany poziomu odpowiednich miRNA, w krwiobiegu lub innych płynach ustrojowych, mogą być markerem diagnostycznym chorób. Diagnostyka onkologiczna mogłaby przebiegać na podstawie badań profilu miRNA pacjenta i porównania go z opracowanymi wcześniej profilami zmian miRNA powiązanymi z występowaniem danego rodzaju choroby nowotworowej. Informacja o zmianach profilu miRNA podstawowych w regulacji ekspresji genów związanych z procesami nowotworzenia, mogłaby się przyczynić do opracowania terapii eksperymentalnych opartych na przywróceniu pierwotnego poziomu miRNA w komórkach, a tym samym, na przywróceniu prawidłowej regulacji ekspresji genów. Coraz nowsze metody wyciszania i włączania ekspresji miRNA mogą w przyszłości zaowocować skutecznymi rozwiązaniami terapeutycznymi.

isomiRs-Hidden Soldiers in the miRNA Regulatory Army, and How to Find Them?

Numerous studies on microRNAs (miRNA) in cancer and other diseases have been accompanied by diverse computational approaches and experimental methods to predict and validate miRNA biological and clinical significance as easily accessible disease biomarkers. In recent years, the application of the next-generation deep sequencing for the analysis and discovery of novel RNA biomarkers has clearly shown an expanding repertoire of diverse sequence variants of mature miRNAs, or isomiRs, resulting from alternative post-transcriptional processing events, and affected by (patho)physiological changes, population origin, individual's gender, and age. Here, we provide an in-depth overview of currently available bioinformatics approaches for the detection and visualization of both mature miRNA and cognate isomiR sequences. An attempt has been made to present in a systematic way the advantages and downsides of in silico approaches in terms of their sensitivity and accuracy performance, as well as used methods, workflows, and processing steps, and end output dataset overlapping issues. The focus is given to the challenges and pitfalls of isomiR expression analysis. Specifically, we address the availability of tools enabling research without extensive bioinformatics background to explore this fascinating corner of the small RNAome universe that may facilitate the discovery of new and more reliable disease biomarkers.

Encyclopedia of tools for the analysis of miRNA isoforms

RNA sequencing data sets rapidly increase in quantity. For microRNAs (miRNAs), frequently dozens to hundreds of billion reads are generated per study. The quantification of annotated miRNAs and the prediction of new miRNAs are leading computational tasks. Now, the increased depth of coverage allows to gain deeper insights into the variability of miRNAs. The analysis of isoforms of miRNAs (isomiRs) is a trending topic, and a range of computational tools for the analysis of isomiRs has been developed. We provide an overview on 27 available computational solutions for the analysis of isomiRs. These include both stand-alone programs (17 tools) and web-based solutions (10 tools) and span a publication time range from 2010 to 2020. Seven of the tools were published in 2019 and 2020, confirming the rising importance of the topic. While most of the analyzed tools work for a broad range of organisms or are completely independent of a reference organism, several tools have been tailored for the analysis of human miRNA data or for plants. While 14 of the tools are general analysis tools of miRNAs, and isomiR analysis is one of their features, the remaining 13 tools have specifically been developed for isomiR analysis. A direct comparison on 20 deep sequencing data sets for selected tools provides insights into the heterogeneity of results. With our work, we provide users a comprehensive overview on the landscape of isomiR analysis tools and in that support the selection of the most appropriate tool for their respective research task.

Long-term exposure of human endothelial cells to metformin modulates miRNAs and isomiRs

Increasing evidence suggest that the glucose-lowering drug metformin exerts a valuable anti-senescence role. The ability of metformin to affect the biogenesis of selected microRNAs (miRNAs) was recently suggested. MicroRNA isoforms (isomiRs) are distinct variations of miRNA sequences, harboring addition or deletion of one or more nucleotides at the 5′ and/or 3′ ends of the canonical miRNA sequence. We performed a comprehensive analysis of miRNA and isomiR profile in human endothelial cells undergoing replicative senescence in presence of metformin. Metformin treatment was associated with the differential expression of 27 miRNAs (including miR-100-5p, -125b-5p, -654-3p, -217 and -216a-3p/5p). IsomiR analysis revealed that almost 40% of the total miRNA pool was composed by non-canonical sequences. Metformin significantly affects the relative abundance of 133 isomiRs, including the non-canonical forms of the aforementioned miRNAs. Pathway enrichment analysis suggested that pathways associated with proliferation and nutrient sensing are modulated by metformin-regulated miRNAs and that some of the regulated isomiRs (e.g. the 5′ miR-217 isomiR) are endowed with alternative seed sequences and share less than half of the predicted targets with the canonical form. Our results show that metformin reshapes the senescence-associated miRNA/isomiR patterns of endothelial cells, thus expanding our insight into the cell senescence molecular machinery.

The Role of Non-Coding RNAs in Uveal Melanoma

Simple Summary

The development of uveal melanoma is a multifactorial and multi-step process, in which abnormal gene expression plays a key role. Recently, several studies have highlighted the role of non-coding RNAs in the progression of uveal melanoma by affecting different signaling pathways. As important agents in the regulation of genes, non-coding RNAs have enormous potential to open up therapeutic pathways, predict response to treatment, and anticipate patient outcome for uveal melanoma. This review aims to provide a comprehensive view of what we know about ncRNAs in uveal melanoma currently.

Abstract

Uveal melanoma (UM) is the most common primary intraocular tumor in adulthood. Approximately 50% of patients develop metastatic disease, which typically affects the liver and is usually fatal within one year. This type of cancer is heterogeneous in nature and is divided into two broad groups of tumors according to their susceptibility to develop metastasis. In the last decade, chromosomal abnormalities and the aberrant expression of several signaling pathways and oncogenes in uveal melanomas have been described. Recently, importance has been given to the association of the mentioned deregulation with the expression of non-coding RNAs (ncRNAs). Here, we review the different classes of ncRNAs—such as long non-coding RNAs (lncRNAs) and microRNAs (miRNAs)—and their contribution to the development of UM. Special attention is given to miRNAs and their regulatory role in physiopathology and their potential as biomarkers. As important agents in gene regulation, ncRNAs have a huge potential for opening up therapeutic pathways, predicting response to treatment, and anticipating patient outcome for UM.

Unification of miRNA and isomiR research: the mirGFF3 format and the mirtop API

MOTIVATION

MicroRNAs (miRNAs) are small RNA molecules (∼22 nucleotide long) involved in post-transcriptional gene regulation. Advances in high-throughput sequencing technologies led to the discovery of isomiRs, which are miRNA sequence variants. While many miRNA-seq analysis tools exist, the diversity of output formats hinders accurate comparisons between tools and precludes data sharing and the development of common downstream analysis methods.

RESULTS

To overcome this situation, we present here a community-based project, miRNA Transcriptomic Open Project (miRTOP) working towards the optimization of miRNA analyses. The aim of miRTOP is to promote the development of downstream isomiR analysis tools that are compatible with existing detection and quantification tools. Based on the existing GFF3 format, we first created a new standard format, mirGFF3, for the output of miRNA/isomiR detection and quantification results from small RNA-seq data. Additionally, we developed a command line Python tool, mirtop, to create and manage the mirGFF3 format. Currently, mirtop can convert into mirGFF3 the outputs of commonly used pipelines, such as seqbuster, isomiR-SEA, sRNAbench, Prost! as well as BAM files. Some tools have also incorporated the mirGFF3 format directly into their code, such as, miRge2.0, IsoMIRmap and OptimiR. Its open architecture enables any tool or pipeline to output or convert results into mirGFF3. Collectively, this isomiR categorization system, along with the accompanying mirGFF3 and mirtop API, provide a comprehensive solution for the standardization of miRNA and isomiR annotation, enabling data sharing, reporting, comparative analyses and benchmarking, while promoting the development of common miRNA methods focusing on downstream steps of miRNA detection, annotation and quantification.

AVAILABILITY AND IMPLEMENTATION

https://github.com/miRTop/mirGFF3/ and https://github.com/miRTop/mirtop.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.