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

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

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

3'IsomiR Species Composition Affects Reliable Quantification of miRNA/isomiR Variants by Poly(A) RT-qPCR: Impact on Small RNA-Seq Profiling Validation

Small RNA-sequencing (small RNA-seq) has revealed the presence of small RNA-naturally occurring variants such as microRNA (miRNA) isoforms or isomiRs. Due to their small size and the sequence similarity among miRNA isoforms, their validation by RT-qPCR is challenging. We previously identified two miR-31-5p isomiRs-the canonical and a 3'isomiR variant (3' G addition)-which were differentially expressed between individuals with azoospermia of different origin. Here, we sought to determine the discriminatory capacity between these two closely-related miRNA isoforms of three alternative poly(A) based-RT-qPCR strategies in both synthetic and real biological context. We found that these poly(A) RT-qPCR strategies exhibit a significant cross-reactivity between these miR-31-5p isomiRs which differ by a single nucleotide, compromising the reliable quantification of individual miRNA isoforms. Fortunately, in the biological context, given that the two miRNA variants show changes in the same direction, RT-qPCR results were consistent with the findings of small RNA-seq study. We suggest that miRNA selection for RT-qPCR validation should be performed with care, prioritizing those canonical miRNAs that, in small RNA-seq, show parallel/homogeneous expression behavior with their most prevalent isomiRs, to avoid confounding RT-qPCR-based results. This is suggested as the current best strategy for robust biomarker selection to develop clinically useful tests.

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.

MicroRNA 3' ends shorten during adolescent brain maturation

MicroRNA (miRNA) dysregulation is well-documented in psychiatric disease, but miRNA dynamics remain poorly understood during adolescent and early adult brain maturation, when symptoms often first appear. Here, we use RNA sequencing to examine miRNAs and their mRNA targets in cortex and hippocampus from early-, mid-, and late-adolescent and adult mice. Furthermore, we use quantitative proteomics by tandem mass tag mass spectrometry (TMT-MS) to examine protein dynamics in cortex from the same subjects. We found that ~25% of miRNAs' 3' ends shorten with age due to increased 3' trimming and decreased U tailing. Particularly, shorter but functionally competent isoforms (isomiRs) of miR-338-3p increase up to 10-fold during adolescence and only in brain. MiRNAs that undergo 3' shortening exhibit stronger negative correlations with targets that decrease with age and stronger positive correlations with targets that increase with age, than miRNAs with stable 3' ends. Increased 3' shortening with age was also observed in available mouse and human miRNA-seq data sets, and stronger correlations between miRNAs that undergo shortening and their mRNA targets were observed in two of the three available data sets. We conclude that age-associated miRNA 3' shortening is a well-conserved feature of postnatal brain maturation.

IsomiR-eQTL: A Cancer-Specific Expression Quantitative Trait Loci Database of miRNAs and Their Isoforms

The identification of expression quantitative trait loci (eQTL) is an important component in efforts to understand how genetic variants influence disease risk. MicroRNAs (miRNAs) are short noncoding RNA molecules capable of regulating the expression of several genes simultaneously. Recently, several novel isomers of miRNAs (isomiRs) that differ slightly in length and sequence composition compared to their canonical miRNAs have been reported. Here we present isomiR-eQTL, a user-friendly database designed to help researchers find single nucleotide polymorphisms (SNPs) that can impact miRNA (miR-eQTL) and isomiR expression (isomiR-eQTL) in 30 cancer types. The isomiR-eQTL includes a total of 152,671 miR-eQTLs and 2,390,805 isomiR-eQTLs at a false discovery rate (FDR) of 0.05. It also includes 65,733 miR-eQTLs overlapping known cancer-associated loci identified through genome-wide association studies (GWAS). To the best of our knowledge, this is the first study investigating the impact of SNPs on isomiR expression at the genome-wide level. This database may pave the way for researchers toward finding a model for personalised medicine in which miRNAs, isomiRs, and genotypes are utilised.

Quantification of microRNA editing using two-tailed RT-qPCR for improved biomarker discovery

Even though microRNAs have been viewed as promising biomarkers for years, their clinical implementation is still lagging far behind. This is in part due to the lack of RT-qPCR technologies that can differentiate between microRNA isoforms. For example, A-to-I editing of microRNAs through adenosine deaminase acting on RNA (ADAR) enzymes can affect their expression levels and functional roles, but editing isoform-specific assays are not commercially available. Here, we describe RT-qPCR assays that are specific for editing isoforms, using microRNA-379 (miR-379) as a model. The assays are based on two-tailed RT-qPCR, and we show them to be compatible both with SYBR Green and hydrolysis-based chemistries, as well as with both qPCR and digital PCR. The assays could readily detect different miR-379 editing isoforms in various human tissues as well as changes of editing levels in ADAR-overexpressing cell lines. We found that the miR-379 editing frequency was higher in prostate cancer samples compared to benign prostatic hyperplasia samples. Furthermore, decreased expression of unedited miR-379, but not edited miR-379, was associated with treatment resistance, metastasis, and shorter overall survival. Taken together, this study presents the first RT-qPCR assays that were demonstrated to distinguish A-to-I-edited microRNAs, and shows that they can be useful in the identification of biomarkers that previously have been masked by other isoforms.

mi-IsoNet: systems-scale microRNA landscape reveals rampant isoform-mediated gain of target interaction diversity and signaling specificity

MicroRNA (miRNA) is not a single sequence, but a series of multiple variants (also termed isomiRs) with sequence and expression heterogeneity. Whether and how these isoforms contribute to functional variation and complexity at the systems and network levels remain largely unknown. To explore this question systematically, we comprehensively analyzed the expression of small RNAs and their target sites to interrogate functional variations between novel isomiRs and their canonical miRNA sequences. Our analyses of the pan-cancer landscape of miRNA expression indicate that multiple isomiRs generated from the same miRNA locus often exhibit remarkable variation in their sequence, expression and function. We interrogated abundant and differentially expressed 5' isomiRs with novel seed sequences via seed shifting and identified many potential novel targets of these 5' isomiRs that would expand interaction capabilities between small RNAs and mRNAs, rewiring regulatory networks and increasing signaling circuit complexity. Further analyses revealed that some miRNA loci might generate diverse dominant isomiRs that often involved isomiRs with varied seeds and arm-switching, suggesting a selective advantage of multiple isomiRs in regulating gene expression. Finally, experimental validation indicated that isomiRs with shifted seed sequences could regulate novel target mRNAs and therefore contribute to regulatory network rewiring. Our analysis uncovers a widespread expansion of isomiR and mRNA interaction networks compared with those seen in canonical small RNA analysis; this expansion suggests global gene regulation network perturbations by alternative small RNA variants or isoforms. Taken together, the variations in isomiRs that occur during miRNA processing and maturation are likely to play a far more complex and plastic role in gene regulation than previously anticipated.

Genistein alleviates chronic vascular inflammatory response via the miR‑21/NF‑κB p65 axis in lipopolysaccharide‑treated mice

Chronic vascular inflammatory response is an important pathological basis of cardiovascular disease. Genistein (GEN), a natural compound, exhibits anti‑inflammatory effects. The aim of the present study was to investigate the effects of GEN on lipopolysaccharide (LPS)‑induced chronic vascular inflammatory response in mice and explore the underlying anti‑inflammatory mechanisms. C57BL/6 mice were fed with a high‑fat diet combined with intraperitoneal injection of LPS to induce chronic vascular inflammation. The expression levels of TNF‑α, IL‑6 and microRNA (miR)‑21 in the vasculature were detected via reverse transcription‑quantitative (RT‑q)PCR. The protein levels of inducible nitric oxide synthase (iNOS) and NF‑κB p65 were detected via western blotting. NF‑κB p65 was also analyzed via immunohistochemistry and immunofluorescence (IF). In addition, after transfection with miR‑21 mimic or inhibitor for 24 h, vascular endothelial cells (VECs) were treated with GEN and LPS. RT‑qPCR and western blot analyses were performed to detect the expression of TNF‑α, IL‑6, miR‑21 and iNOS, and the protein levels of iNOS and NF‑κB p65, respectively. IF was used to measure NF‑κB p65 nuclear translocation. The results revealed that GEN significantly decreased the expression of inflammation‑associated vascular factors in LPS‑treated C57BL/6 mice, including TNF‑α, IL‑6, iNOS, NF‑κB p65 and miR‑21. Furthermore, miR‑21 antagomir enhanced the anti‑inflammatory effects of GEN. In LPS‑induced VECs, miR‑21 mimic increased inflammation‑associated factor expression and attenuated the anti‑inflammatory effects of GEN, whereas miR‑21 inhibitor induced opposing effects. Therefore, the results of the present study suggested that GEN inhibited chronic vascular inflammatory response in mice, which may be associated with the inhibition of VEC inflammatory injury via the miR‑21/NF‑κB p65 pathway.

Dynamic mRNP Remodeling in Response to Internal and External Stimuli

Signal transduction and the regulation of gene expression are fundamental processes in every cell. RNA-binding proteins (RBPs) play a key role in the post-transcriptional modulation of gene expression in response to both internal and external stimuli. However, how signaling pathways regulate the assembly of RBPs with mRNAs remains largely unknown. Here, we summarize observations showing that the formation and composition of messenger ribonucleoprotein particles (mRNPs) is dynamically remodeled in space and time by specific signaling cascades and the resulting post-translational modifications. The integration of signaling events with gene expression is key to the rapid adaptation of cells to environmental changes and stress. Only a combined approach analyzing the signal transduction pathways and the changes in post-transcriptional gene expression they cause will unravel the mechanisms coordinating these important cellular processes.

An estimate of the total number of true human miRNAs

While the number of human miRNA candidates continuously increases, only a few of them are completely characterized and experimentally validated. Toward determining the total number of true miRNAs, we employed a combined in silico high- and experimental low-throughput validation strategy. We collected 28 866 human small RNA sequencing data sets containing 363.7 billion sequencing reads and excluded falsely annotated and low quality data. Our high-throughput analysis identified 65% of 24 127 mature miRNA candidates as likely false-positives. Using northern blotting, we experimentally validated miRBase entries and novel miRNA candidates. By exogenous overexpression of 108 precursors that encode 205 mature miRNAs, we confirmed 68.5% of the miRBase entries with the confirmation rate going up to 94.4% for the high-confidence entries and 18.3% of the novel miRNA candidates. Analyzing endogenous miRNAs, we verified the expression of 8 miRNAs in 12 different human cell lines. In total, we extrapolated 2300 true human mature miRNAs, 1115 of which are currently annotated in miRBase V22. The experimentally validated miRNAs will contribute to revising targetomes hypothesized by utilizing falsely annotated miRNAs.

qPCR-based methods for expression analysis of miRNAs

Several approaches for miRNA expression analysis have been developed in recent years. In this article, we provide an updated and comprehensive review of available qPCR-based methods for miRNA expression analysis and discuss their advantages and disadvantages. Existing techniques involve the use of stem-loop reverse transcriptase-PCR, polyadenylation of RNAs, ligation of adapters or RT with complex primers, using universal or miRNA-specific qPCR primers and/or probes. Many of these methods are oriented towards the expression analysis of mature miRNAs and few are designed for the study of pre-miRNAs and pri-miRNAs. We also discuss findings from articles that compare results from existing methods. Finally, we suggest key points for the improvement of available techniques and for the future development of additional methods.

miR-625-3p and lncRNA GAS5 in Liquid Biopsies for Predicting the Outcome of Malignant Pleural Mesothelioma Patients Treated with Neo-Adjuvant Chemotherapy and Surgery

Combining neo-adjuvant chemotherapy and surgery is part of multimodality treatment of malignant pleural mesothelioma (MPM), but not all patients benefit from this approach. In this exploratory analysis, we investigated the prognostic value of circulating miR-625-3p and lncRNA GAS5 after neo-adjuvant chemotherapy. 36 MPM patients from the SAKK 17/04 trial (NCT00334594), whose blood was available before and after chemotherapy were investigated. RNA was isolated from plasma and reverse transcribed into cDNA. miR-16-5p and β-actin were used as a reference gene for miR-625-3p and GAS5, respectively. After exclusion of samples due to hemolysis or RNA degradation, paired plasma samples from 32 patients before and after chemotherapy were further analyzed. Quantification of miR-625-3p levels in all 64 samples revealed a bimodal distribution and cloning and sequencing of miR-625-3p qPCR product revealed the presence of miR-625-3p isomiRs. Relative change of the circulating miR-625-3p and GAS5 levels after chemotherapy showed that increased circulating miR-625-3p and decreased GAS5 was significantly associated with disease progression (Fisher's test, p = 0.0393). In addition, decreased levels of circulating GAS5 were significantly associated with shorter overall and progression-free survival. Our exploratory analysis revealed a potential value of circulating non-coding RNA for selection of patients likely to benefit from surgery after platinum-based adjuvant chemotherapy.

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.