sRNAflow: A Tool for the Analysis of Small RNA-Seq Data

The analysis of small RNA sequencing data across a range of biofluids is a significant research area, given the diversity of RNA types that hold potential diagnostic, prognostic, and predictive value. The intricate task of segregating the complex mixture of small RNAs from both human and other species, including bacteria, fungi, and viruses, poses one of the most formidable challenges in the analysis of small RNA sequencing data, currently lacking satisfactory solutions. This study introduces sRNAflow, a user-friendly bioinformatic tool with a web interface designed for the analysis of small RNAs obtained from biological fluids. Tailored to the unique requirements of such samples, the proposed pipeline addresses various challenges, including filtering potential RNAs from reagents and environment, classifying small RNA types, managing small RNA annotation overlap, conducting differential expression assays, analysing isomiRs, and presenting an approach to identify the sources of small RNAs within samples. sRNAflow also encompasses an alternative alignment-free analysis of RNA-seq data, featuring clustering and initial RNA source identification using BLAST. This comprehensive approach facilitates meaningful comparisons of results between different analytical methods.

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.

Genetic Regulation of Human isomiR Biogenesis

MicroRNAs play a critical role in regulating gene expression post-transcriptionally. Variations in mature microRNA sequences, known as isomiRs, arise from imprecise cleavage and nucleotide substitution or addition. These isomiRs can target different mRNAs or compete with their canonical counterparts, thereby expanding the scope of miRNA post-transcriptional regulation. Our study investigated the relationship between cis-acting single-nucleotide polymorphisms (SNPs) in precursor miRNA regions and isomiR composition, represented by the ratio of a specific 5'-isomiR subtype to all isomiRs identified for a particular mature miRNA. Significant associations between 95 SNP-isomiR pairs were identified. Of note, rs6505162 was significantly associated with both the 5'-extension of hsa-miR-423-3p and the 5'-trimming of hsa-miR-423-5p. Comparison of breast cancer and normal samples revealed that the expression of both isomiRs was significantly higher in tumors than in normal tissues. This study sheds light on the genetic regulation of isomiR maturation and advances our understanding of post-transcriptional regulation by microRNAs.

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.

A comprehensive pan-cancer analysis reveals cancer-associated robust isomiR expression landscapes in miRNA arm switching

MicroRNAs (miRNAs), important regulators of gene expression, play critical roles in various biological processes and tumorigenesis. To reveal the potential relationships between multiple isomiRs and arm switching, we performed a comprehensive pan-cancer analysis to discuss their roles in tumorigenesis and cancer prognosis. Our results showed that many miR-#-5p and miR-#-3p pairs from the two arms of pre-miRNA may have abundant expression levels, and they are often involved in distinct functional regulatory networks by targeting different mRNAs, although they may also interact with common targets. The two arms may show diverse isomiR expression landscapes, and their expression ratio might vary, mainly depending on tissue type. Dominantly expressed isomiRs can be used to determine distinct cancer subtypes that are associated with clinical outcome, indicating that they may be potential prognostic biomarkers. Our findings indicate robust and flexible isomiR expression landscapes that will enrich the study of miRNAs/isomiRs and aid in revealing the potential roles of multiple isomiRs yielded by arm switching in tumorigenesis.

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.

Structural basis of microRNA biogenesis by Dicer-1 and its partner protein Loqs-PB

In animals and plants, Dicer enzymes collaborate with double-stranded RNA-binding domain (dsRBD) proteins to convert precursor-microRNAs (pre-miRNAs) into miRNA duplexes. We report six cryo-EM structures of Drosophila Dicer-1 that show how Dicer-1 and its partner Loqs‑PB cooperate (1) before binding pre-miRNA, (2) after binding and in a catalytically competent state, (3) after nicking one arm of the pre-miRNA, and (4) following complete dicing and initial product release. Our reconstructions suggest that pre-miRNA binds a rare, open conformation of the Dicer‑1⋅Loqs‑PB heterodimer. The Dicer-1 dsRBD and three Loqs‑PB dsRBDs form a tight belt around the pre-miRNA, distorting the RNA helix to place the scissile phosphodiester bonds in the RNase III active sites. Pre-miRNA cleavage shifts the dsRBDs and partially closes Dicer-1, which may promote product release. Our data suggest a model for how the Dicer‑1⋅Loqs‑PB complex affects a complete cycle of pre-miRNA recognition, stepwise endonuclease cleavage, and product release.

miRNome Profiling and Functional Analysis Reveal Involvement of hsa-miR-1246 in Colon Adenoma-Carcinoma Transition by Targeting AXIN2 and CFTR

Regulatory changes occurring early in colorectal cancer development remain poorly investigated. Since the majority of cases develop from polyps in the adenoma-carcinoma transition, a search of early molecular features, such as aberrations in miRNA expression occurring prior to cancer development, would enable identification of potentially causal, rather than consequential, candidates in the progression of polyp to cancer. In the current study, by employing small RNA-seq profiling of colon biopsy samples, we described differentially expressed miRNAs and their isoforms in the adenoma-carcinoma transition. Analysis of healthy-adenoma-carcinoma sequence in an independent validation group enabled us to identify early deregulated miRNAs including hsa-miR-1246 and hsa-miR-215-5p, the expressions of which are, respectively, gradually increasing and decreasing. Loss-of-function experiments revealed that inhibition of hsa-miR-1246 lead to reduced cell viability, colony formation, and migration rate, thereby indicating an oncogenic effect of this miRNA in vitro. Subsequent western blot and luciferase reporter assay provided evidence of hsa-miR-1246 being involved in the regulation of target AXIN2 and CFTR genes’ expression. To conclude, the present study revealed possible involvement of hsa-miR-1246 in early colorectal cancer development and regulation of tumor suppressors AXIN2 and CFTR.

Global MicroRNA Profiling of Vascular Endothelial Cells

MicroRNA sequencing (miRNA-seq) enables the detection and characterization of the cell miRNome, including miRNA isoforms (isomiRs) and novel miRNA species. In roughly half of the cases, the most abundant isomiR in the cells is not the reference miRNA given in miRBase, which highlights the importance of isomiR-specific analysis. Here, we describe a gel-free protocol for global miRNA profiling in vascular endothelial cells and the main steps of the subsequent data analysis with two alternative analysis methods. In addition to endothelial cells, the protocol is suitable for other cell and tissue types and has been successfully used to obtain miRNA-seq data from human cardiac tissue, plasma, pericardial fluid, and biofluid exosomes.

Protein profiling reveals potential isomiR-associated cross-talks among RNAs in cholangiocarcinoma

•

Protein profiling identified crucial genes that were used to screen related ncRNAs.

•

Both miRNAs and isomiRs were involved in coding-non-coding RNA regulatory network.

•

IsomiR-associated ceRNA networks implied the complex interactions among RNAs.

Cholangiocarcinomas (CCAs) are tumors that arise from the cholangiocytes. Although some genes have been shown with important roles in pathological process, interactions or cross-talks among different RNAs are important to understand the detailed molecular mechanisms in cancer development, especially discussing cross-talks among isomiRs and other RNAs. Herein, to characterize crucial genes in CCA, the protein expression profile was performed to survey potential crucial mRNAs and related non-coding RNAs (ncRNAs) in mRNA-ncRNA network, mainly including miRNAs/isomiRs and lncRNAs. Deregulated mRNAs were firstly obtained if consistent expression patterns were found at protein and mRNA levels, and related miRNAs/isomiRs were screened according to regulatory relationships. Diverse isomiRs from a given miRNA locus also contributed to interactions between the small RNAs and target mRNAs, and miRNAs were further used to survey related lncRNAs to expand the interactions. Thus, several groups of RNAs were constructed as candidate competitive endogenous RNA (ceRNA) networks. Finally, we found that RAB11FIP1:miR-101-3p:MIR3142HG may be a potential ceRNA network, and the interactions among them may be more complex due to variety of isomiRs. Simultaneously, RAB11FIP1 and miR-194-5p were also detected other related lncRNAs (FBXL19-AS1, SNHG1 and PVT1) that may be crucial in coding-non-coding RNA regulatory network. Our results show that diverse isomiRs with sequence and expression heterogeneities contribute to ceRNA regulatory network that may have crucial roles in CCA, which will expand our understanding of interactions among diverse RNAs and their contributions in cancer development.

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.

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.

Comprehensive analysis of the isomiRome in the vegetative organs of the conifer Pinus pinaster under contrasting water availability

An increasing number of microRNAs (miRNAs) and miRNA-related sequences produced during miRNA biogenesis, comprising the isomiRome, have been recently highlighted in different species as critical mediators of environmental stress responses. Conifers have some of the largest known genomes but an extensive characterization of the isomiRome from any conifer species has been lacking. We provide here a comprehensive overview of the Pinus pinaster isomiRome expressed in roots, stem and needles under well-watered and drought conditions. From the 13,441 unique small RNA sequences identified, 2,980 were annotated as canonical miRNAs or miRNA* and the remaining were classified as isomiRNA or miRNA-like sequences. A survey of their expression patterns highlighted roots as the most responsive organ under drought, where specific sequences of which a 24-nt novel miRNA stood out, were strongly down-regulated. Given the putative roles of the miRNA-targeted transcripts validated specifically in root tissues, some of the miRNAs, conserved and novel, are shortlisted as potential regulators of drought response. These results provide a valuable resource for comparative studies between gymnosperms and angiosperms. Furthermore, it evidences high transferability of the isomiRome between pine species being a useful basis for further molecular regulation and physiological studies, and especially those focused on adaptation to drought conditions.

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.

miRNA repertoires of cystic fibrosis ex vivo models highlight miR-181a and miR-101 that regulate WISP1 expression

Cystic fibrosis (CF), a genetic disorder, is characterized by chronic lung disease. Small non-coding RNAs are key regulators of gene expression and participate in various processes, which are dysregulated in CF; however, they remain poorly studied. Here, we determined the complete microRNAs (miRNAs) expression pattern in three CF ex vivo models. The miRNA profiles of air-liquid interface cultures of airway epithelia (bronchi, nasal cells, and nasal polyps) samples from patients with CF and non-CF controls were obtained by deep sequencing. Compared with non-CF controls, several miRNAs were deregulated in CF samples; for instance, miR-181a-5p and the miR-449 family were upregulated. Moreover, mature miRNAs often showed variations (i.e. isomiRs) relative to their reference sequence, such as miR-101, suggesting that miRNAs consist of heterogeneous repertoires of multiple isoforms with different effects on gene expression. Analysis of miR-181a-5p and miR-101-3p roles indicated that they regulate the expression of WISP1, a key component of cell proliferation/migration programs. We showed that miR-101 and miR-181a-5p participated in aberrant recapitulation of wound healing programs by controlling WISP1 mRNA and protein level. Our miRNA expression data bring new insights into CF physiopathology and define new potential therapeutic targets in CF. © 2020 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

microRNA-seq of cartilage reveals an overabundance of miR-140-3p which contains functional isomiRs

miR-140 is selectively expressed in cartilage. Deletion of the entire Mir140 locus in mice results in growth retardation and early-onset osteoarthritis-like pathology; however, the relative contribution of miR-140-5p or miR-140-3p to the phenotype remains to be determined. An unbiased small RNA sequencing approach identified miR-140-3p as significantly more abundant (>10-fold) than miR-140-5p in human cartilage. Analysis of these data identified multiple miR-140-3p isomiRs differing from the miRBase annotation at both the 5' and 3' end, with >99% having one of two seed sequences (5' bases 2-8). Canonical (miR-140-3p.2) and shifted (miR-140-3p.1) seed isomiRs were overexpressed in chondrocytes and transcriptomics performed to identify targets. miR-140-3p.1 and miR-140-3p.2 significantly down-regulated 694 and 238 genes, respectively, of which only 162 genes were commonly down-regulated. IsomiR targets were validated using 3'UTR luciferase assays. miR-140-3p.1 targets were enriched within up-regulated genes in rib chondrocytes of Mir140-null mice and within down-regulated genes during human chondrogenesis. Finally, through imputing the expression of miR-140 from the expression of the host gene WWP2 in 124 previously published data sets, an inverse correlation with miR-140-3p.1 predicted targets was identified. Together these data suggest the novel seed containing isomiR miR-140-3p.1 is more functional than original consensus miR-140-3p seed containing isomiR.

Differences in the Drosha and Dicer Cleavage Profiles in Colorectal Cancer and Normal Colon Tissue Samples

Human colorectal adenocarcinoma cell line Caco-2 is often used as a model of healthy intestinal epithelium, in particular, in miRNA studies. The work of the enzymes Drosha and Dicer is an integral part of the process of miRNA formation. Inaccuracies in the work of these enzymes lead to a change in the nucleotide sequences of miRNAs with the formation of new isoforms, which, in turn, can change intracellular regulatory mechanisms. In the framework of this study, it was shown that the quantitative estimates of inaccuracies in Drosha and Dicer activity significantly differ between the specimens of normal colon tissue and malignant colorectal tumors.

A comprehensive overview of bull sperm-borne small non-coding RNAs and their diversity across breeds

Background

Mature sperm carry thousands of RNAs, including mRNAs, lncRNAs, tRNAs, rRNAs and sncRNAs, though their functional significance is still a matter of debate. Growing evidence suggests that sperm RNAs, especially sncRNAs, are selectively retained during spermiogenesis or specifically transferred during epididymis maturation, and are thus delivered to the oocyte at fertilization, providing resources for embryo development. However , a deep characterization of the sncRNA content of bull sperm and its expression profile across breeds is currently lacking. To fill this gap, we optimized a guanidinium–Trizol total RNA extraction protocol to prepare high-quality RNA from frozen bull sperm collected from 40 representative bulls from six breeds. Deep sequencing was performed (40 M single 50-bp reads per sample) to establish a comprehensive repertoire of cattle sperm sncRNA.

Results

Our study showed that it comprises mostly piRNAs (26%), rRNA fragments (25%), miRNAs (20%) and tRNA fragments (tsRNA, 14%). We identified 5p-halves as the predominant tsRNA subgroup in bull sperm, originating mostly from Gly and Glu isoacceptors. Our study also increased by ~ 50% the sperm repertoire of known miRNAs and identified 2022 predicted miRNAs. About 20% of sperm miRNAs were located within genomic clusters, expanding the list of known polycistronic pri-miRNA clusters and defining several networks of co-expressed miRNAs. Strikingly, our study highlighted the great diversity of isomiRs, resulting mainly from deletions and non-templated additions (A and U) at the 3p end. Substitutions within miRNA sequence accounted for 40% of isomiRs, with G>A, U>C and C>U substitutions being the most frequent variations. In addition, many sncRNAs were found to be differentially expressed across breeds.

Conclusions

Our study provides a comprehensive overview of cattle sperm sncRNA, and these findings will pave the way for future work on the role of sncRNAs in embryo development and their relevance as biomarkers of semen fertility.

MicroRNA-411 and Its 5'-IsomiR Have Distinct Targets and Functions and Are Differentially Regulated in the Vasculature under Ischemia

MicroRNAs are posttranscriptional regulators of gene expression. As microRNAs can target many genes simultaneously, microRNAs can regulate complex multifactorial processes, including post-ischemic neovascularization, a major recovery pathway in cardiovascular disease. MicroRNAs select their target mRNAs via full complementary binding with their seed sequence, i.e., nucleotides 2-8 from the 5' end of a microRNA. The exact sequence of a mature microRNA, and thus of its 5' and 3' ends, is determined by two sequential cleavage steps of microRNA precursors, Drosha/DGCR8 and Dicer. When these cleavage steps result in nucleotide switches at the 5' end, forming a so-called 5'-isomiR, this results in a shift in the mature microRNA's seed sequence. The role of 5'-isomiRs in cardiovascular diseases is still unknown. Here, we characterize the expression and function of the 5'-isomiR of miR-411 (ISO-miR-411). ISO-miR-411 is abundantly expressed in human primary vascular cells. ISO-miR-411 has a different "targetome" from WT-miR-411, with only minor overlap. The ISO-miR-411/WT-miR-411 ratio is downregulated under acute ischemia, both in cells and a murine ischemia model, but is upregulated instead in chronically ischemic human blood vessels. ISO-miR-411 negatively influences vascular cell migration, whereas WT-miR-411 does not. Our data demonstrate that isomiR formation is a functional pathway that is actively regulated during ischemia.

Silver, titanium dioxide, and zinc oxide nanoparticles trigger miRNA/isomiR expression changes in THP-1 cells that are proportional to their health hazard potential

After over a decade of nanosafety research, it is indisputable that the vast majority of nano-sized particles induce a plethora of adverse cellular responses - the severity of which is linked to the material's physicochemical properties. Differentiated THP-1 cells were previously exposed for 6 h and 24 h to silver, titanium dioxide, and zinc oxide nanoparticles at the maximum molar concentration at which no more than 15% cellular cytotoxicity was observed. All three nanoparticles differed in extent of induction of biological pathways corresponding to immune response signaling and metal ion homeostasis. In this study, we integrated gene and miRNA expression profiles from the same cells to propose miRNA biomarkers of adverse exposure to metal-based nanoparticles. We employed RNA sequencing together with a quantitative strategy that also enables analysis of the overlooked repertoire of length and sequence miRNA variants called isomiRs. Whilst only modest changes in expression were observed within the first 6 h of exposure, the miRNA/isomiR (miR) profiles of each nanoparticle were unique. Via canonical correlation and pathway enrichment analyses, we identified a co-regulated miR-mRNA cluster, predicted to be highly relevant for cellular response to metal ion homeostasis. These miRs were annotated to be canonical or variant isoforms of hsa-miR-142-5p, -342-3p, -5100, -6087, -6894-3p, and -7704. Hsa-miR-5100 was differentially expressed in response to each nanoparticle in both the 6 h and 24 h exposures. Taken together, this co-regulated miR-mRNA cluster could represent potential biomarkers of sub-toxic metal-based nanoparticle exposure.

IsomiRs and tRNA-derived fragments are associated with metastasis and patient survival in uveal melanoma

Uveal melanoma (UVM) is the most common primary intraocular malignancy in adults. With over 50% of patients developing metastatic disease, there is an unmet need for improved diagnostic and therapeutic options. Efforts to understand the molecular biology of the disease have revealed several markers that correlate with patient prognosis, including the copy number of chromosome 3, genetic alterations in the BAP1, EIF1AX and SF3B1 genes, and other transcriptional features. Here, we expand upon previous reports by comprehensively characterizing the short RNA-ome in 80 primary UVM tumor samples. In particular, we describe a previously unseen complex network involving numerous regulatory molecules that comprise microRNA (miRNAs), novel UVM-specific miRNA loci, miRNA isoforms (isomiRs), and tRNA-derived fragments (tRFs). Importantly, we show that the abundance profiles of isomiRs and tRFs associate with various molecular phenotypes, metastatic disease, and patient survival. Our findings suggest deep involvement of isomiRs and tRFs in the disease etiology of UVM. We posit that further study and characterization of these novel molecules will improve understanding of the mechanisms underlying UVM, and lead to the development of new diagnostic and therapeutic approaches.

Poly(A)-specific ribonuclease sculpts the 3' ends of microRNAs

The 3' ends of metazoan microRNAs (miRNAs) are initially defined by the RNase III enzymes during maturation, but subsequently experience extensive modifications by several enzymatic activities. For example, terminal nucleotidyltransferases (TENTs) elongate miRNAs by adding one or a few nucleotides to their 3' ends, which occasionally leads to differential regulation of miRNA stability or function. However, the catalytic entities that shorten miRNAs and the molecular consequences of such shortening are less well understood, especially in vertebrates. Here, we report that poly(A)-specific ribonuclease (PARN) sculpts the 3' ends of miRNAs in human cells. By generating PARN knockout cells and characterizing their miRNAome, we demonstrate that PARN digests the 3' extensions of miRNAs that are derived from the genome or attached by TENTs, thereby effectively reducing the length of miRNAs. Surprisingly, PARN-mediated shortening has little impact on miRNA stability, suggesting that this process likely operates to finalize miRNA maturation, rather than to initiate miRNA decay. PARN-mediated shortening is pervasive across most miRNAs and appears to be a conserved mechanism contributing to the 3' end formation of vertebrate miRNAs. Our findings add miRNAs to the expanding list of noncoding RNAs whose 3' end formation depends on PARN.

Small RNA expression from viruses, bacteria and human miRNAs in colon cancer tissue and its association with microsatellite instability and tumor location

BACKGROUND

MicroRNAs (miRNA) and other small RNAs are frequently dysregulated in cancer and are promising biomarkers for colon cancer. Here we profile human, virus and bacteria small RNAs in normal and tumor tissue from early stage colon cancer and correlate the expression with clinical parameters.

METHODS

Small RNAs from colon cancer tissue and adjacent normal mucosa of 48 patients were sequenced using Illumina high-throughput sequencing. Clinical parameters were correlated with the small RNA expression data using linear models. We performed a meta-analysis by comparing publicly available small RNA sequencing datasets with our original sequencing data to confirm the main findings.

RESULTS

We identified 331 differentially expressed miRNAs between tumor and normal samples. We found that the major changes in miRNA expression between left and right colon are due to miRNAs located within the Hox-developmental genes, including miR-10b, miR-196b and miR-615. Further, we identified new miRNAs associated with microsatellite instability (MSI), including miR-335, miR-26 and miR-625. We performed a meta-analysis on all publicly available miRNA-seq datasets and identified 117 common miRNAs that were differentially expressed between tumor and normal tissue. The miRNAs miR-135b and miR-31 were the most significant upregulated miRNA in tumor across all datasets. The miRNA miR-133a was the most strongly downregulated miRNA in our dataset and also showed consistent downregulation in the other datasets. The miRNAs associated with MSI and tumor location in our data showed similar changes in the other datasets. Finally, we show that small RNAs from Epstein-Barr virus and Fusobacterium nucleatum are differentially expressed between tumor and normal adjacent tissue.

CONCLUSIONS

Small RNA profiling in colon cancer tissue revealed novel RNAs associated with MSI and tumor location. We show that Fusobacterium nucleatum are detectable at the RNA-level in colon tissue, and that both Fusobacterium nucleatum and Epstein-Barr virus separate tumor and normal tissue.

Tumor classification and biomarker discovery based on the 5'isomiR expression level

Background

The miRNA isoforms (isomiRs) have been suggested to regulate the same pathways as the canonical miRNA and play an important biological role in miRNA-mediated gene regulation. Recently, a study has demonstrated that the presence or absence of all isomiRs could efficiently discriminate amongst 32 TCGA cancer types. Besides, an effective reduction of distinguishing isomiR features for multiclass tumor discrimination must have a major impact on our understanding of the disease and treatment of cancer.

Methods

In this study, we have constructed a combination of the genetic algorithms (GA) with Random Forest (RF) algorithms to detect reliable sets of cancer-associated 5’isomiRs from TCGA isomiR expression data for multiclass tumor classification.

Results

We obtained 100 sets of the optimal predictive features, each of which comprised of 50–5’isomiRs that could effectively classify with an average sensitivity of 92% samples from 32 different tumor types. We calculated the frequency with which a 5’isomiR found in these sets as measuring its importance for tumor classification. Many highly frequent 5’isomiRs with different 5′ loci from canonical miRNAs were detected in these sets, supporting that the isomiRs play a significant role in the multiclass tumor classification. The further functional enrichment analysis showed that the target genes of the 10 most frequently appearing 5’isomiRs were involved in the activity of transcription activator and protein kinase and cell-cell adhesion.

Conclusions

The findings of the present study indicated that the 5’isomiRs might be employed for multiclass tumor classification and the suggested that GA/RF model could perform effective tumor classification by a series of largely independent optimal predictor 5′ isomiR sets.

miRNA length variation during macrophage stimulation confounds the interpretation of results: implications for miRNA quantification by RT-qPCR

Most microRNAs (miRNAs) are expressed as a mix of length isoforms (referred to as isomiRs). IsomiR stoichiometry can be differentially impacted upon cell stimulation, as recently evidenced by our group in the context of immune responses induced by type-I interferon (IFN). Here, we revisit published RNA-seq data sets of human and mouse macrophages stimulated with bacterial products at the isomiR level. We demonstrate that for several miRNAs, macrophage stimulation induces changes in isomiR stoichiometry. Critically, we find that changes in miRNA expression can be misinterpreted when miRNAs are quantified by RT-qPCR, as primers directed against canonical miRNA sequences may not equally target the different isomiRs that are regulated endogenously. Beyond the case of phagocyte stimulation, our analyses reinforce the concept that analysis of miRNA expression at the isoform level should become standard practice.

Global microRNA and isomiR expression associated with liver metabolism is induced by organophosphorus flame retardant exposure in male Chinese rare minnow (Gobiocypris rarus)

To reveal the adverse effects of organophosphorus flame retardants (OPFRs) on aquatic organisms at the epigenetic level, male Chinese rare minnows were exposed to 0.24 mg/L tris(2‑butoxyethyl) phosphate (TBOEP), 0.04 mg/L tris(1,3‑dichloro‑2‑propyl) phosphate (TDCIPP), or 0.012 mg/L triphenyl phosphate (TPHP) for 14 days. The effects of sub-acute OPFR exposure on liver miRNA and the 3' isomiR expression profiles of Chinese rare minnows were investigated. Through small RNA sequencing and bioinformatics analysis, a total of 32, 84, and 19 differentially expressed miRNAs were detected for TBOEP, TDCIPP, and TPHP exposure, respectively (p < 0.05). Target prediction of the differentially expressed miRNAs and pathway enrichment analysis indicated that predicted altered mRNAs for all three OPFRs were associated with metabolic pathways, whereas base excision repair was only predicted to be perturbed by the TPHP treatment. In addition, 3' isomiR-Us were unexpectedly abundant in all groups (e.g., miR-143), and TDCIPP strongly increased the ratio of 3' isomiR-U expression. Finally, histological examination and metabolic enzyme activity analyses werein agreement with the predicted metabolic pathways. As such, our study indicates that the investigation of epigenetics changes in miRNA gene transcription is a considerable method for the assessment of aquatic toxicity.

miRge 2.0 for comprehensive analysis of microRNA sequencing data

Background

miRNAs play important roles in the regulation of gene expression. The rapidly developing field of microRNA sequencing (miRNA-seq; small RNA-seq) needs comprehensive, robust, user-friendly and standardized bioinformatics tools to analyze these large datasets. We present miRge 2.0, in which multiple enhancements were made towards these goals.

Results

miRge 2.0 has become more comprehensive with increased functionality including a novel miRNA detection method, A-to-I editing analysis, integrated standardized GFF3 isomiR reporting, and improved alignment to miRNAs. The novel miRNA detection method uniquely uses both miRNA hairpin sequence structure and composition of isomiRs resulting in higher specificity for potential miRNA identification. Using known miRNA data, our support vector machine (SVM) model predicted miRNAs with an average Matthews correlation coefficient (MCC) of 0.939 over 32 human cell datasets and outperformed miRDeep2 and miRAnalyzer regarding phylogenetic conservation. The A-to-I editing detection strongly correlated with a reference dataset with adjusted R² = 0.96. miRge 2.0 is the most up-to-date aligner with custom libraries to both miRBase v22 and MirGeneDB v2.0 for 6 species: human, mouse, rat, fruit fly, nematode and zebrafish; and has a tool to create custom libraries. For user-friendliness, miRge 2.0 is incorporated into bcbio-nextgen and implementable through Bioconda.

Conclusions

miRge 2.0 is a redesigned, leading miRNA RNA-seq aligner with several improvements and novel utilities. miRge 2.0 is freely available at: https://github.com/mhalushka/miRge.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2287-y) contains supplementary material, which is available to authorized users.

Dumbbell-PCR for Discriminative Quantification of a Small RNA Variant

Cellular RNAs are often expressed as multiple isoforms of complex heterogeneity in both length and terminal sequences. IsomiRs, the isoforms of microRNAs, are such an example. Distinct quantification of each RNA variant is necessary to unravel the biogenesis mechanism and biological significance of heterogenetic RNA expression. Here we describe Dumbbell-PCR (Db-PCR), a TaqMan RT-PCR-based method that distinctively quantifies a specific small RNA variant with single-nucleotide resolution at terminal sequences. Db-PCR enables the quantitative analysis of RNA terminal heterogeneity without performing Next-Generation Sequencing.

Strong conservation of inbred mouse strain microRNA loci but broad variation in brain microRNAs due to RNA editing and isomiR expression

Diversity in the structure and expression of microRNAs, important regulators of gene expression, arises from SNPs, duplications followed by divergence, production of isomiRs, and RNA editing. Inbred mouse strains and crosses using them are important reference populations for genetic mapping, and as models of human disease. We determined the nature and extent of interstrain miRNA variation by (i) identifying miRNA SNPs in whole-genome sequence data from 36 strains, and (ii) examining miRNA editing and expression in hippocampus (Hpc) and frontal cortex (FCx) of six strains, to facilitate the study of miRNAs in neurobehavioral phenotypes. miRNA loci were strongly conserved among the 36 strains, but even the highly conserved seed region contained 16 SNPs. In contrast, we identified RNA editing in 58.9% of miRNAs, including 11 consistent editing events in the seed region. We confirmed the functional significance of three conserved edits in the miR-379/410 cluster, demonstrating that edited miRNAs gained novel target mRNAs not recognized by the unedited miRNAs. We found significant interstrain differences in miRNA and isomiR expression: Of 779 miRNAs expressed in Hpc and 719 in FCx, 262 were differentially expressed (190 in Hpc, 126 in FCx, 54 in both). We also identified 32 novel miRNA candidates using miRNA prediction tools. Our studies provide the first comprehensive analysis of SNP, isomiR, and RNA editing variation in miRNA loci across inbred mouse strains, and a detailed catalog of expressed miRNAs in Hpc and FCx in six commonly used strains. These findings will facilitate the molecular analysis of neurological and behavioral phenotypes in this model organism.

miR-222 isoforms are differentially regulated by type-I interferon

Endogenous microRNAs (miRNAs) often exist as multiple isoforms (known as "isomiRs") with predominant variation around their 3'-end. Increasing evidence suggests that different isomiRs of the same family can have diverse functional roles, as recently demonstrated with the example of miR-222-3p 3'-end variants. While isomiR levels from a same miRNA family can vary between tissues and cell types, change of templated isomiR stoichiometry to stimulation has not been reported to date. Relying on small RNA-sequencing analyses, we demonstrate here that miR-222-3p 3'-end variants >23 nt are specifically decreased upon interferon (IFN) β stimulation of human fibroblasts, while shorter isoforms are spared. This length-dependent dynamic regulation of long miR-222-3p 3'-isoforms and >40 other miRNA families was confirmed in human monocyte-derived dendritic cells following infection with Salmonella Typhimurium, underlining the breadth of 3'-length regulation by infection, beyond the example of miR-222-3p. We further show that stem-loop miRNA Taqman RT-qPCR exhibits selectivity between 3'-isoforms, according to their length, and that this can lead to misinterpretation of results when these isoforms are differentially regulated. Collectively, and to our knowledge, this work constitutes the first demonstration that the stoichiometry of highly abundant templated 3'-isoforms of a same miRNA family can be dynamically regulated by a stimulus. Given that such 3'-isomiRs can have different functions, our study underlines the need to consider isomiRs when investigating miRNA-based regulation.

MicroRNA dynamics at the onset of primordial germ and somatic cell sex differentiation during mouse embryonic gonad development

In mammals, commitment and specification of germ cell lines involves complex programs that include sex differentiation, control of proliferation, and meiotic initiation. Regulation of these processes is genetically controlled by fine-tuned mechanisms of gene regulation in which microRNAs (miRNAs) are involved. We have characterized, by small-RNA-seq and bioinformatics analyses, the miRNA expression patterns of male and female mouse primordial germ cells (PGCs) and gonadal somatic cells at embryonic stages E11.5, E12.5, and E13.5. Differential expression analyses revealed differences in the regulation of key miRNA clusters such as miR-199-214, miR-182-183-96, and miR-34c-5p, whose targets have defined roles during gonadal sexual determination in both germ and somatic cells. Extensive analyses of miRNA sequences revealed an increase in noncanonical isoforms on PGCs at E12.5 and dramatic changes of 3' isomiR expression and 3' nontemplate nucleotide additions in female PGCs at E13.5. Additionally, RT-qPCR analyses of genes encoding proteins involved in miRNA biogenesis and 3' nucleotide addition uncovered sexually and developmentally specific expression, characterized by the decay of Drosha, Dgcr8, and Xpo5 expression along gonadal development. These results demonstrate that miRNAs, their isomiRs, and miRNA machinery are differentially regulated and participate actively in gonadal sexual differentiation in both PGCs and gonadal somatic cells.

Small Molecules in Rare Tumors: Emerging Role of MicroRNAs in GIST

Gastrointestinal stromal tumors (GISTs) are the most common mesenchymal tumors of gastrointestinal tract. GISTs have very different clinical phenotypes and underlying molecular characteristics that are not yet completely understood. microRNAs (miRNAs) have been shown to participate in carcinogenesis pathways through post-transcriptional regulation of gene expression in different tumors. Over the last years emerging evidence has highlighted the role of miRNAs in GISTs. This review provides an overview of original research papers that analyze miRNA deregulation patterns, functional role, diagnostic, therapeutic and prognostic implications in GIST as well as provides directions for further research in the field.

Modified Polyadenylation-Based RT-qPCR Increases Selectivity of Amplification of 3'-MicroRNA Isoforms

MicroRNA (miRNA) detection by reverse transcription (RT) quantitative real-time PCR (RT-qPCR) is the most popular method currently used to measure miRNA expression. Although the majority of miRNA families are constituted of several 3′-end length variants (“isomiRs”), little attention has been paid to their differential detection by RT-qPCR. However, recent evidence indicates that 3′-end miRNA isoforms can exhibit 3′-length specific regulatory functions, underlining the need to develop strategies to differentiate 3′-isomiRs by RT-qPCR approaches. We demonstrate here that polyadenylation-based RT-qPCR strategies targeted to 20–21 nt isoforms amplify entire miRNA families, but that primers targeted to >22 nt isoforms were specific to >21 nt isoforms. Based on this observation, we developed a simple method to increase selectivity of polyadenylation-based RT-qPCR assays toward shorter isoforms, and demonstrate its capacity to help distinguish short RNAs from longer ones, using synthetic RNAs and biological samples with altered isomiR stoichiometry. Our approach can be adapted to many polyadenylation-based RT-qPCR technologies already exiting, providing a convenient way to distinguish long and short 3′-isomiRs.

The dynamic recruitment of TRBP to neuronal membranes mediates dendritogenesis during development

MicroRNAs are important regulators of local protein synthesis during neuronal development. We investigated the dynamic regulation of microRNA production and found that the majority of the microRNA‐generating complex, consisting of Dicer, TRBP, and PACT, specifically associates with intracellular membranes in developing neurons. Stimulation with brain‐derived neurotrophic factor (BDNF), which promotes dendritogenesis, caused the redistribution of TRBP from the endoplasmic reticulum into the cytoplasm, and its dissociation from Dicer, in a Ca²⁺‐dependent manner. As a result, the processing of a subset of neuronal precursor microRNAs, among them the dendritically localized pre‐miR16, was impaired. Decreased production of miR‐16‐5p, which targeted the BDNF mRNA itself, was rescued by expression of a membrane‐targeted TRBP. Moreover, miR‐16‐5p or membrane‐targeted TRBP expression blocked BDNF‐induced dendritogenesis, demonstrating the importance of neuronal TRBP dynamics for activity‐dependent neuronal development. We propose that neurons employ specialized mechanisms to modulate local gene expression in dendrites, via the dynamic regulation of microRNA biogenesis factors at intracellular membranes of the endoplasmic reticulum, which in turn is crucial for neuronal dendrite complexity and therefore neuronal circuit formation and function.

Altered microRNA profiles during early colon adenoma progression in a porcine model of familial adenomatous polyposis

MicroRNAs are dysregulated in various cancers including colorectal cancer, and are potential useful biomarkers of disease development. We used next generation sequencing to investigate miRNA expression profiles in low- and high-grade intraepithelial dysplastic polyps from pigs carrying a mutation in the adenomatous polyposis coli tumour suppressor (APC¹³¹¹ , orthologous to human APC¹³⁰⁹ ) that model an inherited predisposition to colorectal cancer, familial adenomatous polyposis. We identified several miRNAs and their isomiRs significantly (P < 0.05) differentially expressed between low and high-grade intraepithelial dysplastic polyps. Of these, ssc-let-7e, ssc-miR-98, ssc-miR-146a-5p, ssc-miR-146b, ssc-miR-183 and ssc-miR-196a were expressed at higher level and ssc-miR-126-3p at lower level in high-grade intraepithelial dysplastic polyps. Functional miRNA target analysis revealed significant (P < 0.001) over-representation of cancer-related pathways, including 'microRNAs in cancer', 'proteoglycans in cancer', 'pathways in cancer' and 'colorectal cancer'. This is the first study to reveal miRNAs associated with premalignant transformation of colon polyps.

Identification of metastasis-associated microRNAs in serum from rectal cancer patients

MicroRNAs (miRNAs) are promising prognostic and diagnostic biomarkers due to their high stability in blood. Here we investigate the expression of miRNAs and other noncoding (nc) RNAs in serum of rectal cancer patients. Serum from 96 rectal cancer patients was profiled using small RNA sequencing and expression of small RNAs was correlated with the clinicopathological characteristics of the patients. Multiple classes of RNAs were detected, including miRNAs and fragments of tRNAs, snoRNAs, long ncRNAs, and other classes of RNAs. Several miRNAs, miRNA variants (isomiRs) and other ncRNAs were differentially expressed between Stage IV and Stage I-III rectal cancer patients, including several members of the miR-320 family. Furthermore, we show that high expression of miR-320d as well as one tRNA fragment is associated with poor survival. We also show that several miRNAs and isomiRs are differentially expressed between patients receiving preoperative chemoradiotherapy and patients who did not receive any treatment before serum collection. In summary, our study shows that the expression of miRNAs and other small ncRNAs in serum may be used to predict distant metastasis and survival in rectal cancer.

Assessment of isomiR Discrimination Using Commercial qPCR Methods

We sought to determine whether commercial quantitative polymerase chain reaction (qPCR) methods are capable of distinguishing isomiRs: variants of mature microRNAs (miRNAs) with sequence endpoint differences. We used two commercially available miRNA qPCR methods to quantify miR-21-5p in both synthetic and real cell contexts. We find that although these miRNA qPCR methods possess high sensitivity for specific sequences, they also pick up background signals from closely related isomiRs, which influences the reliable quantification of individual isomiRs. We conclude that these methods do not possess the requisite specificity for reliable isomiR quantification.

Characterizing isomiR variants within the microRNA-34/449 family

miR‐34/449 microRNAs are conserved regulators of multiciliated cell differentiation. Here, we evidence and characterize expression of two isomiR variant sequences from the miR‐34/449 family in human airway epithelial cells. These isomiRs differ from their canonical counterparts miR‐34b and miR‐449c by one supplemental uridine at their 5′‐end, leading to a one‐base shift in their seed region. Overexpression of canonical miR‐34/449 or 5′‐isomiR‐34/449 induces distinct gene expression profiles and biological effects. However, some target transcripts and functional activities are shared by both canonical microRNAs and isomiRs. Indeed, both repress important targets that result in cell cycle blockage and Notch pathway inhibition. Our findings suggest that 5′‐isomiR‐34/449 may represent additional mechanisms by which miR‐34/449 family finely controls several pathways to drive multiciliogenesis.

Characterization of miRNAs involved in response to poly(I:C) in porcine airway epithelial cells

MicroRNAs (miRNA) have been implicated in a variety of pathological conditions including infectious diseases. Knowledge of the miRNAs affected by poly(I:C), a synthetic analog of viral double-stranded RNA, in porcine airway epithelial cells (PAECs) contributes to understanding the mechanisms of swine viral respiratory diseases, which bring enormous economic loss worldwide every year. In this study, we used high throughput sequencing to profile miRNA expression in PAECs treated with poly(I:C) as compared to the untreated control. This approach revealed 23 differentially expressed miRNAs (DEMs), five of which have not been implicated in viral infection before. Nineteen of the 23 miRNAs were down-regulated including members of the miR-17-92 cluster, a well-known polycistronic oncomir and extensively involved in viral infection in humans. Target genes of DEMs, predicted using bioinformatic methods and validated by luciferase reporter analysis on two representative DEMs, were significantly enriched in several pathways including transforming growth factor-β signaling. A large quantity of sequence variations (isomiRs) were found including a substitution at position 5, which was verified to redirect miRNAs to a new spectrum of targets by luciferase reporter assay together with bioinformatics analysis. Twelve novel porcine miRNAs conserved in other species were identified by homology analysis together with cloning verification. Furthermore, the expression analysis revealed the potential importance of three novel miRNAs in porcine immune response to viruses. Overall, our data contribute to clarifying the mechanisms underlying the host immune response against respiratory viruses in pigs, and enriches the repertoire of porcine miRNAs.

isomiR-SEA: an RNA-Seq analysis tool for miRNAs/isomiRs expression level profiling and miRNA-mRNA interaction sites evaluation

Background

Massive parallel sequencing of transcriptomes, revealed the presence of many miRNAs and miRNAs variants named isomiRs with a potential role in several cellular processes through their interaction with a target mRNA. Many methods and tools have been recently devised to detect and quantify miRNAs from sequencing data. However, all of them are implemented on top of general purpose alignment methods, thus providing poorly accurate results and no information concerning isomiRs and conserved miRNA-mRNA interaction sites.

Results

To overcome these limitations we present a novel algorithm named isomiR-SEA, that is able to provide users with very accurate miRNAs expression levels and both isomiRs and miRNA-mRNA interaction sites precise classifications. Tags are mapped on the known miRNAs sequences thanks to a specialized alignment algorithm developed on top of biological evidence concerning miRNAs structure. Specifically, isomiR-SEA checks for miRNA seed presence in the input tags and evaluates, during all the alignment phases, the positions of the encountered mismatches, thus allowing to distinguish among the different isomiRs and conserved miRNA-mRNA interaction sites.

Conclusions

isomiR-SEA performances have been assessed on two public RNA-Seq datasets proving that the implemented algorithm is able to account for more reliable and accurate miRNAs expression levels with respect to those provided by two compared state of the art tools. Moreover, differently from the few methods currently available to perform isomiRs detection, the proposed algorithm implements the evaluation of isomiRs and conserved miRNA-mRNA interaction sites already in the first alignment phases, thus avoiding any additional filtering stages potentially responsible for the loss of useful information.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-0958-0) contains supplementary material, which is available to authorized users.

Considerations for optimization of microRNA PCR assays for molecular diagnosis

The remarkable stability of microRNAs in biofluids underlies their potential as biomarkers, but their small size presents challenges for detection by RT-qPCR. The heterogeneity of microRNAs, with each one comprising a series of variants or 'isomiRs', adds additional complexity. Presented here are the key considerations for use of RT-qPCR to measure microRNAs and their isomiRs, with a focus on plasma. Modified nucleotides can be incorporated into primer sequences to enhance affinity and provide increased specificity and sensitivity for RT-qPCR assays. Approaches based upon polyA tailing and use of a common oligo(dT)-based reverse transcription oligonucleotide will detect most isomiRs. Conversely, stem-loop RT oligonucleotides and sequence specific probes can enable detection of specific isomiRs of interest. Next generation sequencing of all the products of a microRNA RT-PCR reaction is a promising new approach for both microRNA quantification and characterization.

miRNA Nomenclature: A View Incorporating Genetic Origins, Biosynthetic Pathways, and Sequence Variants

High-throughput sequencing of miRNAs has revealed the diversity and variability of mature and functional short noncoding RNAs, including their genomic origins, biogenesis pathways, sequence variability, and newly identified products such as miRNA-offset RNAs (moRs). Here we review known cases of alternative mature miRNA-like RNA fragments and propose a revised definition of miRNAs to encompass this diversity. We then review nomenclature guidelines for miRNAs and propose to extend nomenclature conventions to align with those for protein-coding genes established by international consortia. Finally, we suggest a system to encompass the full complexity of sequence variations (i.e., isomiRs) in the analysis of small RNA sequencing experiments.

Divergent target recognition by coexpressed 5'-isomiRs of miR-142-3p and selective viral mimicry

Sequence heterogeneity at the ends of mature microRNAs (miRNAs) is well documented, but its effects on miRNA function are largely unexplored. Here we studied the impact of miRNA 5'-heterogeneity, which affects the seed region critical for target recognition. Using the example of miR-142-3p, an emerging regulator of the hematopoietic lineage in vertebrates, we show that naturally coexpressed 5'-variants (5'-isomiRs) can recognize largely distinct sets of binding sites. Despite this, both miR-142-3p isomiRs regulate exclusive and shared targets involved in actin dynamics. Thus, 5'-heterogeneity can substantially broaden and enhance regulation of one pathway. Other 5'-isomiRs, in contrast, recognize largely overlapping sets of binding sites. This is exemplified by two herpesviral 5'-isomiRs that selectively mimic one of the miR-142-3p 5'-isomiRs. We hypothesize that other cellular and viral 5'-isomiRs can similarly be grouped into those with divergent or convergent target repertoires, based on 5'-sequence features. Taken together, our results provide a detailed characterization of target recognition by miR-142-3p and its 5'-isomiR-specific viral mimic. We furthermore demonstrate that miRNA 5'-end variation leads to differential targeting and can thus broaden the target range of miRNAs.

A comprehensive characterization of the nuclear microRNA repertoire of post-mitotic neurons

MicroRNAs (miRNAs) are small non-coding RNAs with important functions in the development and plasticity of post-mitotic neurons. In addition to the well-described cytoplasmic function of miRNAs in post-transcriptional gene regulation, recent studies suggested that miRNAs could also be involved in transcriptional and post-transcriptional regulatory processes in the nuclei of proliferating cells. However, whether miRNAs localize to and function within the nucleus of post-mitotic neurons is unknown. Using a combination of microarray hybridization and small RNA deep sequencing, we identified a specific subset of miRNAs which are enriched in the nuclei of neurons. Nuclear enrichment of specific candidate miRNAs (miR-25 and miR-92a) could be independently validated by Northern blot, quantitative real-time PCR (qRT-PCR) and fluorescence in situ hybridization (FISH). By cross-comparison to published reports, we found that nuclear accumulation of miRNAs might be linked to a down-regulation of miRNA expression during in vitro development of cortical neurons. Importantly, by generating a comprehensive isomiR profile of the nuclear and cytoplasmic compartments, we found a significant overrepresentation of guanine nucleotides (nt) at the 3′-terminus of nuclear-enriched isomiRs, suggesting the presence of neuron-specific mechanisms involved in miRNA nuclear localization. In conclusion, our results provide a starting point for future studies addressing the nuclear function of specific miRNAs and the detailed mechanisms underlying subcellular localization of miRNAs in neurons and possibly other polarized cell types.