Selected isomiR expression profiles via arm switching?

The intricacies of isomiRs: from classification to clinical relevance

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

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.

Age-associated changes in microRNAs affect the differentiation potential of human mesenchymal stem cells: Novel role of miR-29b-1-5p expression

Age-associated osteoporosis is widely accepted as involving the disruption of osteogenic stem cell populations and their functioning. Maintenance of the local bone marrow (BM) microenvironment is critical for regulating proliferation and differentiation of the multipotent BM mesenchymal stromal/stem cell (BMSC) population with age. The potential role of microRNAs (miRNAs) in modulating BMSCs and the BM microenvironment has recently gained attention. However, miRNAs expressed in rapidly isolated BMSCs that are naïve to the non-physiologic standard tissue culture conditions and reflect a more accurate in vivo profile have not yet been reported. Here we directly isolated CD271 positive (+) BMSCs within hours from human surgical BM aspirates without culturing and performed microarray analysis to identify the age-associated changes in BMSC miRNA expression. One hundred and two miRNAs showed differential expression with aging. Target prediction and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed that the up-regulated miRNAs targeting genes in bone development pathways were considerably enriched. Among the differentially up-regulated miRNAs the novel passenger strand miR-29b-1-5p was abundantly expressed as a mature functional miRNA with aging. This suggests a critical arm-switching mechanism regulates the expression of the miR-29b-1-5p/3p pair shifting the normally degraded arm, miR-29b-1-5p, to be the dominantly expressed miRNA of the pair in aging. The normal guide strand miR-29b-1-3p is known to act as a pro-osteogenic miRNA. On the other hand, overexpression of the passenger strand miR-29b-1-5p in culture-expanded CD271+ BMSCs significantly down-regulated the expression of stromal cell-derived factor 1 (CXCL12)/ C-X-C chemokine receptor type 4 (SDF-1(CXCL12)/CXCR4) axis and other osteogenic genes including bone morphogenetic protein-2 (BMP-2) and runt-related transcription factor 2 (RUNX2). In contrast, blocking of miR-29b-1-5p function using an antagomir inhibitor up-regulated expression of BMP-2 and RUNX2 genes. Functional assays confirmed that miR-29b-1-5p negatively regulates BMSC osteogenesis in vitro. These novel findings provide evidence of a pathogenic anti-osteogenic role for miR-29b-1-5p and other miRNAs in age-related defects in osteogenesis and bone regeneration.

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.

Overview of host miRNA properties and their association with epigenetics, long non-coding RNAs, and Xeno-infectious factors

MicroRNA-derived structures play impressive roles in various biological processes. So dysregulation of miRNAs can lead to different human diseases. Recent studies have extended our comprehension of the control of miRNA function and features. Here, we overview some remarkable miRNA properties that have potential implications for the miRNA functions, including different variants of a miRNA called isomiRs, miRNA arm selection/arm switching, and the effect of these factors on miRNA target selection. Besides, we review some aspects of miRNA interactions such as the interaction between epigenetics and miRNA (different miRNAs and their related processing enzymes are epigenetically regulated by multiple DNA methylation enzymes. moreover, DNA methylation could be controlled by diverse mechanisms related to miRNAs), direct and indirect crosstalk between miRNA and lnc (Long Non-Coding) RNAs as a further approach to conduct intercellular regulation called "competing endogenous RNA" (ceRNA) that is involved in the pathogenesis of different diseases, and the interaction of miRNA activities and some Xeno-infectious (virus/bacteria/parasite) factors, which result in modulation of the pathogenesis of infections. This review provides some related studies to a better understanding of miRNA involvement mechanisms and overcoming the complexity of related diseases that may be applicable and useful to prognostic, diagnostic, therapeutic purposes and personalized medicine in the future.

miRSwitch: detecting microRNA arm shift and switch events

Arm selection, the preferential expression of a 3′ or 5′ mature microRNA (miRNA), is a highly dynamic and tissue-specific process. Time-dependent expression shifts or switches between the arms are also relevant for human diseases. We present miRSwitch, a web server to facilitate the analysis and interpretation of arm selection events. Our species-independent tool evaluates pre-processed small non-coding RNA sequencing (sncRNA-seq) data, i.e. expression matrices or output files from miRNA quantification tools (miRDeep2, miRMaster, sRNAbench). miRSwitch highlights potential changes in the distribution of mature miRNAs from the same precursor. Group comparisons from one or several user-provided annotations (e.g. disease states) are possible. Results can be dynamically adjusted by choosing from a continuous range of highly specific to very sensitive parameters. Users can compare potential arm shifts in the provided data to a human reference map of pre-computed arm shift frequencies. We created this map from 46 tissues and 30 521 samples. As case studies we present novel arm shift information in a Alzheimer’s disease biomarker data set and from a comparison of tissues in Homo sapiens and Mus musculus. In summary, miRSwitch offers a broad range of customized arm switch analyses along with comprehensive visualizations, and is freely available at: https://www.ccb.uni-saarland.de/mirswitch/.

Expressing a Target Mimic of miR156fhl-3p Enhances Rice Blast Disease Resistance Without Yield Penalty by Improving SPL14 Expression

MicroRNAs (miRNAs) play essential roles in the regulation of plant growth and defense responses. More and more, miRNA-3ps are reported to act in plant development and immunity. miR156 is a conserved miRNA, and most previous studies focus on its roles in plant growth, development, and yield determinacy. Here, we show that expressing a target mimic of miR156fhl-3p led to enhanced rice blast disease resistance without a yield penalty. miR156fhl-3p was differentially responsive to Magnaporthe oryzae in susceptible and resistant accessions. Transgenic lines expressing a target mimic of miR156fhl-3p (MIM156-3p) exhibited enhanced rice blast disease resistance and increased expression of defense-related genes. MIM156-3p also enhanced the mRNA abundance of SPL14 and WRKY45 by down-regulating miR156-5p and pre-miR156. Moreover, MIM156-3p lines displayed a decreased number of second rachis branches per panicle but enlarged grains, leading to unchanged yield per plant. Consistently, overexpressing miR156h (OX156) led to enhanced susceptibility to M. oryzae and decreased the expression of SPL14 and WRKY45. Our results indicate that miR156fhl-3p mounts a regulatory role on miR156-5p, which subsequently regulates the expression of SPL14 and WRKY45 to improve rice blast disease resistance.

Transient up-regulation of miR-155-3p by lipopolysaccharide in primary human monocyte-derived macrophages results in RISC incorporation but does not alter TNF expression

Background: The innate immune response is a tightly regulated process that reacts rapidly in response to pathogen-associated molecular patterns (PAMPs) such as lipopolysaccharide (LPS). Evidence is accumulating that microRNAs contribute to this, although few studies have examined the early events that constitute the “primary” response.

Methods: LPS-dependent changes to miRNA expression were studied in primary human monocyte-derived macrophages (1°MDMs). An unbiased screen by microarray was validated by qPCR and a method for the absolute quantitation of miRNAs was also developed, utilising 5’ phosphorylated RNA oligonucleotide templates. RNA immunoprecipitation was performed to explore incorporation of miRNAs into the RNA-induced silencing complex (RISC). The effect of miRNA functional inhibition on TNF expression (mRNA and secretion) was investigated.

Results: Of the 197 miRNAs expressed in 1°MDMs, only five were induced >1.5-fold. The most strongly induced was miR-155-3p, the partner strand to miR-155-5p, which are both derived from the MIR155HG/BIC gene (pri-miR-155). The abundance of miR-155-3p was induced transiently ~250-fold at 2-4hrs and then returned towards baseline, mirroring pri-miR-155. Other PAMPs, IL-1β, and TNF caused similar responses. IL-10, NF-κB, and JNK inhibition reduced these responses, unlike cytokine-suppressing mycolactone. Absolute quantitation revealed that miRNA abundance varies widely from donor-to-donor, and showed that miR-155-3p abundance is substantially less than miR-155-5p in unstimulated cells. However, at its peak there were 446-1,113 copies/cell, and miR-155-3p was incorporated into the RISC with an efficiency similar to miR-16-5p and miR-155-5p. Inhibition of neither miRNA affected TNF secretion after 2hrs in 1°MDMs, but technical challenges here are noted.

Conclusions: Dynamic regulation of miRNAs during the primary response is rare, with the exception of miR-155-3p. Further work is required to establish whether its low abundance, even at the transient peak, is sufficient for biological activity and to determine whether there are specific mechanisms determining its biogenesis from miR-155 precursors

Transient up-regulation of miR-155-3p by lipopolysaccharide in primary human monocyte-derived macrophages results in RISC incorporation but does not alter TNF expression

Background: The innate immune response is a tightly regulated process that reacts rapidly in response to pathogen-associated molecular patterns (PAMPs) such as lipopolysaccharide (LPS). Evidence is accumulating that microRNAs contribute to this, although few studies have examined the early events that constitute the “primary” response. Methods: LPS-dependent changes to miRNA expression were studied in primary human monocyte-derived macrophages (1°MDMs). An unbiased screen by microarray was validated by qPCR and a method for the absolute quantitation of miRNAs was also developed, utilising 5’ phosphorylated RNA oligonucleotide templates. RNA immunoprecipitation was performed to explore incorporation of miRNAs into the RNA-induced silencing complex (RISC). The effect of miRNA functional inhibition on TNF expression (mRNA and secretion) was investigated. Results: Of the 197 miRNAs expressed in 1°MDMs, only five were induced >1.5-fold. The most strongly induced was miR-155-3p, the partner strand to miR-155-5p, which are both derived from the BIC gene (B cell integration cluster, MIR155HG). The abundance of miR-155-3p was induced transiently ~250-fold at 2-4hrs and then returned towards baseline, mirroring the BIC mRNA. Other PAMPs, IL-1β, and TNF caused similar responses. IL-10, NF-κB, and JNK inhibition suppressed these responses, unlike cytokine-suppressing mycolactone. Absolute quantitation showed that miRNA abundance varies widely from donor-to-donor, and showed that miR-155-3p abundance is substantially less than miR-155-5p in unstimulated cells. However, at its peak there were 446-1,113 copies/cell, and miR-155-3p was incorporated into the RISC with an efficiency similar to miR-16-5p and miR-155-5p. Inhibition of neither miRNA affected TNF expression in 1°MDMs, but technical challenges here are noted. Conclusions: Dynamic regulation of miRNAs during the primary response is rare, with the exception of miR-155-3p, which transiently achieves levels that might have a biological effect. Further work on this candidate would need to overcome the technical challenges of the broad-ranging effects of liposomes on 1°MDMs.

Substantial Dysregulation of miRNA Passenger Strands Underlies the Vascular Response to Injury

Vascular smooth muscle cell (VSMC) dedifferentiation is a common feature of vascular disorders leading to pro-migratory and proliferative phenotypes, a process induced through growth factor and cytokine signaling cascades. Recently, many studies have demonstrated that small non-coding RNAs (miRNAs) can induce phenotypic effects on VSMCs in response to vessel injury. However, most studies have focused on the contribution of individual miRNAs. Our study aimed to conduct a detailed and unbiased analysis of both guide and passenger miRNA expression in vascular cells in vitro and disease models in vivo. We analyzed 100 miRNA stem loops by TaqMan Low Density Array (TLDA) from primary VSMCs in vitro. Intriguingly, we found that a larger proportion of the passenger strands was significantly dysregulated compared to the guide strands after exposure to pathological stimuli, such as platelet-derived growth factor (PDGF) and IL-1α. Similar findings were observed in response to injury in porcine vein grafts and stent models in vivo. In these studies, we reveal that the miRNA passenger strands are predominantly dysregulated in response to vascular injury.

miRNA arm switching identifies novel tumour biomarkers

Background

microRNAs have been reported to play critical roles in cancer and to have potential as diagnostic biomarkers. During miRNA biogenesis, one strand of the miRNA hairpin precursor is preferentially selected as a functionally mature miRNA, while the other strand is typically degraded. Arm switching occurs when the strand preference is changed. This preference can be different and can change dynamically depending upon the species, tissue types, or development stages. Due to recent advances in next-generation sequencing methods, arm switching has been observed in a variety of cancers.

Methods

A tumour miRNA-Seq dataset was collected from The Cancer Genome Atlas (TCGA). The support vector machine (SVM) method combined with 5-fold cross validation was applied to select the best combination of arm-switched miRNA tumour markers. Survival analysis was also applied to identify patient survival associated miRNA markers.

Findings

We observed 51 arm-switched miRNAs and of these, 7 were associated with patient survival. Twenty-three 1-combination arm switching miRNAs with excellent diagnostic value were identified. Interestingly, ovarian cancer showed a significant difference in arm switching pattern compared with 32 other cancers.

Interpretation

These results suggest that arm switching miRNAs could be used as potential biomarkers for various cancers.

Fund

This work was partially supported by the National Natural Science Foundation of China (no. 61472158, 61572227), and University of Macau Faculty of Health Sciences (MYRG2016-00101-FHS).

Interplay between RNA interference and heat shock response systems in Drosophila melanogaster

The genome expression pattern is strongly modified during the heat shock response (HSR) to form an adaptive state. This may be partly achieved by modulating microRNA levels that control the expression of a great number of genes that are embedded within the gene circuitry. Here, we investigated the cross-talk between two highly conserved and universal house-keeping systems, the HSR and microRNA machinery, in Drosophila melanogaster We demonstrated that pronounced interstrain differences in the microRNA levels are alleviated after heat shock (HS) to form a uniform microRNA pattern. However, individual strains exhibit different patterns of microRNA expression during the course of recovery. Importantly, HS-regulated microRNAs may target functionally similar HS-responsive genes involved in the HSR. Despite the observed general downregulation of primary microRNA precursor expression as well as core microRNA pathway genes after HS, the levels of many mature microRNAs are upregulated. This indicates that the regulation of miRNA expression after HS occurs at transcriptional and post-transcriptional levels. It was also shown that deletion of all hsp70 genes had no significant effect on microRNA biogenesis but might influence the dynamics of microRNA expression during the HSR.

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.

IsomiR expression patterns in canonical and Dicer‑independent microRNAs

Multiple microRNA (miRNA) variants, known as isomiRs, are extensively distributed in miRNA loci and predominantly derive from the alternative cleavage of Drosha/Dicer and 3′addition events. The present study aimed to investigate the expression patterns of multiple isomiRs in typical miRNA and Dicer-independent miRNA loci by conducting evolutionary and expression analysis using public datasets. Although different miRNA maturation processes exist, multiple isomiRs can be detected by similar expression distributions. However, isomiR expression in Dicer-independent miRNA loci tends to be at a moderate level, particularly for random distribution in the ends that are split by Dicer in the typical miRNA loci. Compared with the mature miRNA locus (dominant miRNA locus), the non-dominant miRNA locus indicates an expression distribution similar to that of the Dicer-independent miRNA locus. These results increase the understanding of multiple isomiRs in the progression of diseases.

Loss of miR-26a-5p promotes proliferation, migration, and invasion in prostate cancer through negatively regulating SERBP1

The biological role of miR-26a involved in the carcinogenesis of prostate cancer (PC) has been controversial. Besides, the underlying mechanism by which miR-26a plays a role in PC has been unclear. To investigate the role of miR-26a-5p in the PC, miR-26a-5p was detected and statistically analyzed in clinical PC tissues and a panel of PC cell lines. Using bioinformatics analysis, we found that serpine1 messenger RNA (mRNA) binding protein 1 (SERBP1) was a potential downstream target of miR-26a-5p. Using luciferase reporter and western blot, we identified that miR-26a-5p negatively regulated SERBP1 on the PC cell line level. It was confirmed that miR-26a-5p was markedly downregulated in PC tissues compared with normal controls whose reduced expression was significantly associated with metastasis and poor overall prognosis and found that miR-26a-5p was able to prevent proliferation and motility of PC cells in vitro. Additionally, SERBP1 was identified as a downstream target of miR-26a-5p. Moreover, it was observed that SERBP1 was markedly upregulated in prostate cancer tissues and was significantly associated with tissue metastasis and Gleason score. Taken together, our results for the first time demonstrate that the loss of miR-26a-5p promotes proliferation, migration, and invasion through targeting SERBP1 in PC, supporting the tumor-suppressing role of miR-26a-5p in PC.

Dengue virus infection alters post-transcriptional modification of microRNAs in the mosquito vector Aedes aegypti

Recent discoveries regarding the importance of isomiRs have increased our understanding of the regulatory complexities of the miRNAome. Observed changes in the miRNA profiles in mosquitoes infected with flaviviruses have implicated small RNAs in the interactions between viruses and their vectors. Here we analysed the isomiR profiles of both uninfected and infected Aedes aegypti mosquitoes with the major human pathogen dengue virus (DENV). We found that several specific isomiRs were significantly altered in their abundance patterns in response to DENV infection potentially affecting their target repertoire. Notable among these were isomiR variants which displayed arm-switching. We also demonstrate that modifications to the 3p end of miRNAs are vastly more prevalent than those at the 5p ends. We also observed that in only 45% of Ae. aegypti miRNAs the most abundant read matches the exact sequence reported in miRBase. Further, we found positive correlations between the number of mature miRNA reads, pre-miRNA length, GC content and secondary structure minimum free energy with the number of isomiRs. The findings presented here provide some evidence that isomiR production is not a random phenomenon and may be important in DENV replication in its vector.

mirPRo-a novel standalone program for differential expression and variation analysis of miRNAs

Being involved in many important biological processes, miRNAs can regulate gene expression by targeting mRNAs to facilitate their degradation or translational inhibition. Many miRNA sequencing studies reveal that miRNA variations such as isomiRs and “arm switching” are biologically relevant. However, existing standalone tools usually do not provide comprehensive, detailed information on miRNA variations. To deepen our understanding of miRNA variability, we developed a new standalone tool called “mirPRo” to quantify known miRNAs and predict novel miRNAs. Compared with the most widely used standalone program, miRDeep2, mirPRo offers several new functions including read cataloging based on genome annotation, optional seed region check, miRNA family expression quantification, isomiR identification and categorization, and “arm switching” detection. Our comparative data analyses using three datasets from mouse, human and chicken demonstrate that mirPRo is more accurate than miRDeep2 by avoiding over-counting of sequence reads and by implementing different approaches in adapter trimming, mapping and quantification. mirPRo is an open-source standalone program (https://sourceforge.net/projects/mirpro/).

A Systematic Evaluation of Feature Selection and Classification Algorithms Using Simulated and Real miRNA Sequencing Data

Sequencing is widely used to discover associations between microRNAs (miRNAs) and diseases. However, the negative binomial distribution (NB) and high dimensionality of data obtained using sequencing can lead to low-power results and low reproducibility. Several statistical learning algorithms have been proposed to address sequencing data, and although evaluation of these methods is essential, such studies are relatively rare. The performance of seven feature selection (FS) algorithms, including baySeq, DESeq, edgeR, the rank sum test, lasso, particle swarm optimistic decision tree, and random forest (RF), was compared by simulation under different conditions based on the difference of the mean, the dispersion parameter of the NB, and the signal to noise ratio. Real data were used to evaluate the performance of RF, logistic regression, and support vector machine. Based on the simulation and real data, we discuss the behaviour of the FS and classification algorithms. The Apriori algorithm identified frequent item sets (mir-133a, mir-133b, mir-183, mir-937, and mir-96) from among the deregulated miRNAs of six datasets from The Cancer Genomics Atlas. Taking these findings altogether and considering computational memory requirements, we propose a strategy that combines edgeR and DESeq for large sample sizes.

Bioinformatic Interrogation of 5p-arm and 3p-arm Specific miRNA Expression Using TCGA Datasets

MicroRNAs (miRNAs) play important roles in cellular functions and developmental processes. They are also implicated in oncogenesis mechanisms and could serve as potential cancer biomarkers. Using high-throughput miRNA sequencing information, expression of both the 5p-arm and 3p-arm mature miRNAs were demonstrated and generated from the single miRNA hairpin precursor. However, current miRNA annotations lack comprehensive 5p-arm/3p-arm feature annotations. Among known human mature miRNAs, only half of them are annotated with arm features. This generated ambiguous results in many miRNA-Sequencing (miRNA-Seq) studies. In this report, we have interrogated the TCGA (the Cancer Genome Atlas) miRNA expression datasets with an improved, fully annotated human 5p-arm and 3p-arm miRNA reference list. By utilizing this comprehensive miRNA arm-feature annotations, enhanced determinations and clear annotations were achieved for the miRNA isoforms (isomiRs) recognized from the sequencing reads. In the gastric cancer (STAD) dataset, as an example, 32 5p-arm/3p-arm OPEN ACCESS J. Clin. Med. 2015, 4 1799 specific miRNAs were found to be down-regulated and 24 5p-arm/3p-arm specific miRNAs were found to be up-regulated. We have further extended miRNA biomarker discoveries to additional TCGA miRNA-Seq datasets and provided extensive expression information on 5p-arm/3p-arm miRNAs across multiple cancer types. Our results identified several miRNAs that could be potential common biomarkers for human cancers.

Similarity computation strategies in the microRNA-disease network: a survey

Various microRNAs have been demonstrated to play roles in a number of human diseases. Several microRNA-disease network reconstruction methods have been used to describe the association from a systems biology perspective. The key problem for the network is the similarity computation model. In this article, we reviewed the main similarity computation methods and discussed these methods and future works. This survey may prompt and guide systems biology and bioinformatics researchers to build more perfect microRNA-disease associations and may make the network relationship clear for medical researchers.

Evolutionary and expression analysis of miR-#-5p and miR-#-3p at the miRNAs/isomiRs levels

We mainly discussed miR-#-5p and miR-#-3p under three aspects: (1) primary evolutionary analysis of human miRNAs; (2) evolutionary analysis of miRNAs from different arms across the typical 10 vertebrates; (3) expression pattern analysis of miRNAs at the miRNA/isomiR levels using public small RNA sequencing datasets. We found that no bias can be detected between the numbers of 5p-miRNA and 3p-miRNA, while miRNAs from miR-#-5p and miR-#-3p show variable nucleotide compositions. IsomiR expression profiles from the two arms are always stable, but isomiR expressions in diseased samples are prone to show larger degree of dispersion. miR-#-5p and miR-#-3p have relative independent evolution/expression patterns and datasets of target mRNAs, which might also contribute to the phenomena of arm selection and/or arm switching. Simultaneously, miRNA/isomiR expression profiles may be regulated via arm selection and/or arm switching, and the dynamic miRNAome and isomiRome will adapt to functional and/or evolutionary pressures. A comprehensive analysis and further experimental study at the miRNA/isomiR levels are quite necessary for miRNA study.

Entropy-based model for miRNA isoform analysis

MiRNAs have been widely studied due to their important post-transcriptional regulatory roles in gene expression. Many reports have demonstrated the evidence of miRNA isoform products (isomiRs) in high-throughput small RNA sequencing data. However, the biological function involved in these molecules is still not well investigated. Here, we developed a Shannon entropy-based model to estimate isomiR expression profiles of high-throughput small RNA sequencing data extracted from miRBase webserver. By using the Kolmogorov-Smirnov statistical test (KS test), we demonstrated that the 5p and 3p miRNAs present more variants than the single arm miRNAs. We also found that the isomiR variant, except the 3’ isomiR variant, is strongly correlated with Minimum Free Energy (MFE) of pre-miRNA, suggesting the intrinsic feature of pre-miRNA should be one of the important factors for the miRNA regulation. The functional enrichment analysis showed that the miRNAs with high variation, particularly the 5’ end variation, are enriched in a set of critical functions, supporting these molecules should not be randomly produced. Our results provide a probabilistic framework for miRNA isoforms analysis, and give functional insights into pre-miRNA processing.

Comparisons of isomiR patterns and classification performance using the rank-based MANOVA and 10-fold cross-validation

Next generation sequencing technology has identified a series of miRNA variants (named "isomiRs"), which might be associated with cancer progression. We provide a new strategy to reanalyze the miR-seq datasets through a view of the isomiR spectrum. Firstly, differentially expressed (DE) isomiRs were detected with the DESeq algorithm based on negative binomial distribution. Secondly, the rank-based MANOVA was adopted to compare the isomiR patterns between normal and tumor tissues. Moreover, a comprehensive survey on classification performance of three features was conducted, including the logistic regression, k-nearest neighbors and Random Forest. Finally, functional enrichment analysis was performed with the putative targets of specific isomiRs to elucidate their biological functions. Furthermore, the methods were applied to the downloaded miR-seq datasets of breast invasive carcinoma from TCGA. We found that the expression levels of multiple isomiRs derived from the same miRNA locus showed significant inconsistency between normal and tumor samples. In most cases, logistic regression with multiple DE isomiRs was superior to the others, with highest AUC and lowest AIC. Similarly, DE isomiRs performed best in the average accuracy of standard classifiers. Integrated targets were significantly enriched in some cancer-related pathways, including MAPK signaling pathway, and focal adhesion. Collectively, we could recommend the rank-based MANOVA for comparing different isomiR patterns, and further investigation on isomiRs needs to be considered in miRNA sequencing research.

miRSeq: a user-friendly standalone toolkit for sequencing quality evaluation and miRNA profiling

MicroRNAs (miRNAs) present diverse regulatory functions in a wide range of biological activities. Studies on miRNA functions generally depend on determining miRNA expression profiles between libraries by using a next-generation sequencing (NGS) platform. Currently, several online web services are developed to provide small RNA NGS data analysis. However, the submission of large amounts of NGS data, conversion of data format, and limited availability of species bring problems. In this study, we developed miRSeq to provide alternatives. To test the performance, we had small RNA NGS data from four species, including human, rat, fly, and nematode, analyzed with miRSeq. The alignments results indicate that miRSeq can precisely evaluate the sequencing quality of samples regarding percentage of self-ligation read, read length distribution, and read category. miRSeq is a user-friendly standalone toolkit featuring a graphical user interface (GUI). After a simple installation, users can easily operate miRSeq on a PC or laptop by using a mouse. Within minutes, miRSeq yields useful miRNA data, including miRNA expression profiles, 3′ end modification patterns, and isomiR forms. Moreover, miRSeq supports the analysis of up to 105 animal species, providing higher flexibility.

Virus versus host plant microRNAs: who determines the outcome of the interaction?

Considering the importance of microRNAs (miRNAs) in the regulation of essential processes in plant pathogen interactions, it is not surprising that, while plant miRNA sequences counteract viral attack via antiviral RNA silencing, viruses in turn have developed antihost defense mechanisms blocking these RNA silencing pathways and establish a counter-defense. In the current study, computational and stem-loop Reverse Transcription – Polymerase Chain Reaction (RT-PCR) approaches were employed to a) predict and validate virus encoded mature miRNAs (miRs) in 39 DNA-A sequences of the bipartite genomes of African cassava mosaic virus (ACMV) and East African cassava mosaic virus-Uganda (EACMV-UG) isolates, b) determine whether virus encoded miRs/miRs* generated from the 5′/3′ harpin arms have the capacity to bind to genomic sequences of the host plants Jatropha or cassava and c) investigate whether plant encoded miR/miR* sequences have the potential to bind to the viral genomes. Different viral pre-miRNA hairpin sequences and viral miR/miR* length variants occurring as isomiRs were predicted in both viruses. These miRNAs were located in three Open Reading Frames (ORFs) and in the Intergenic Region (IR). Moreover, various target genes for miRNAs from both viruses were predicted and annotated in the host plant genomes indicating that they are involved in biotic response, metabolic pathways and transcription factors. Plant miRs/miRs* from conserved and highly expressed families were identified, which were shown to have potential targets in the genome of both begomoviruses, representing potential plant miRNAs mediating antiviral defense. This is the first assessment of predicted viral miRs/miRs* of ACMV and EACMV-UG and host plant miRNAs, providing a reference point for miRNA identification in pathogens and their hosts. These findings will improve the understanding of host- pathogen interaction pathways and the function of viral miRNAs in Euphorbiaceous crop plants.

A challenge for miRNA: multiple isomiRs in miRNAomics

Accumulating evidence suggests that a single microRNA (miRNA) locus can generate a series of sequences during miRNA maturation process. These multiple sequences, called miRNA variants, or isomiRs, have different lengths and different 5' and 3' ends. Some of these isomiRs are detected as varied nucleotides and 3' additional non-template nucleotides. As physiological miRNA isoforms, they have drawn attention for possible regulatory biological roles. The present work mainly reviews miRNA/isomiR biogenesis, isomiR expression patterns, and functional and evolutionary implications, especially between isomiRs from homologous and clustered miRNA loci. The phenomenon of multiple isomiRs and their biological roles indicates that analysis performed at the miRNA and isomiR levels should be included in miRNA studies. This may enrich and complicate miRNA biogenesis and coding-non-coding RNA regulatory networks.

Flexible microRNA arm selection in rice

MicroRNAs act at the post-transcriptional level and guide Argonaute proteins to cleave their corresponding target transcripts. However, little attention has been paid to arm selection in miRNA precursors. In this study, small RNA high-throughput sequencing data from 29 different rice libraries were pooled to investigate tissue- and abiotic stress-specific dynamic expression of miRNAs. We found that more than half of pre-miRNAs showed changes in arm selection in different tissues and/or under different abiotic stresses. Our findings suggest that miRNA selection is remarkably prevalent in plants, providing new insights into the role of miRNAs in plant growth and development.