IsomiRage: From Functional Classification to Differential Expression of miRNA Isoforms

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.

Computational tools supporting known miRNA identification

The study of small RNAs is a field that is expanding quickly. Other functional short RNA molecules other than microRNAs, and gene expression regulators, have been found in animals and plants. MicroRNAs play a significant role in host-microbe interactions, and parasite microRNAs may affect the host's innate immunity. Furthermore, short RNAs are intriguing non-invasive biomarker possibilities because they can be found in physiological fluids. These trends suggest that for many researchers, quick and simple techniques for expression profiling and subsequent downstream analysis of miRNA-seq data are crucial. We selected sRNAtoolbox to make integrated sRNA research easier. Each tool can be used separately or to explore and analyze sRNAbench results in further depth. A special focus was placed on the tools' usability. We review available miRNA research tools to have an overview of the evaluation of the tools. Mainly we evaluate the tool sRNAtoolbox.

Phospho-RNA-Seq Highlights Specific Small RNA Profiles in Plasma Extracellular Vesicles

Small RNAs (sRNAs) are bioactive molecules that can be detected in biofluids, reflecting physiological and pathological states. In plasma, sRNAs are found within extracellular vesicles (EVs) and in extravesicular compartments, offering potential sources of highly sensitive biomarkers. Deep sequencing strategies to profile sRNAs favor the detection of microRNAs (miRNAs), the best-known class of sRNAs. Phospho-RNA-seq, through the enzymatic treatment of sRNAs with T4 polynucleotide kinase (T4-PNK), has been recently developed to increase the detection of thousands of previously inaccessible RNAs. In this study, we investigated the value of phospho-RNA-seq on both the EVs and extravesicular plasma subfractions. Phospho-RNA-seq increased the proportion of sRNAs used for alignment and highlighted the diversity of the sRNA transcriptome. Unsupervised clustering analysis using sRNA counts matrices correctly classified the EVs and extravesicular samples only in the T4-PNK treated samples, indicating that phospho-RNA-seq stresses the features of sRNAs in each plasma subfraction. Furthermore, T4-PNK treatment emphasized specific miRNA variants differing in the 5'-end (5'-isomiRs) and certain types of tRNA fragments in each plasma fraction. Phospho-RNA-seq increased the number of tissue-specific messenger RNA (mRNA) fragments in the EVs compared with the extravesicular fraction, suggesting that phospho-RNA-seq favors the discovery of tissue-specific sRNAs in EVs. Overall, the present data emphasizes the value of phospho-RNA-seq in uncovering RNA-based biomarkers in EVs.

miRNAs and isomiRs: Serum-Based Biomarkers for the Development of Intellectual Disability and Autism Spectrum Disorder in Tuberous Sclerosis Complex

Tuberous sclerosis complex (TSC) is a rare multi-system genetic disorder characterized by a high incidence of epilepsy and neuropsychiatric manifestations known as tuberous-sclerosis-associated neuropsychiatric disorders (TANDs), including autism spectrum disorder (ASD) and intellectual disability (ID). MicroRNAs (miRNAs) are small regulatory non-coding RNAs that regulate the expression of more than 60% of all protein-coding genes in humans and have been reported to be dysregulated in several diseases, including TSC. In the current study, RNA sequencing analysis was performed to define the miRNA and isoform (isomiR) expression patterns in serum. A Receiver Operating Characteristic (ROC) curve analysis was used to identify circulating molecular biomarkers, miRNAs, and isomiRs, able to discriminate the development of neuropsychiatric comorbidity, either ASD, ID, or ASD + ID, in patients with TSC. Part of our bioinformatics predictions was verified with RT-qPCR performed on RNA isolated from patients’ serum. Our results support the notion that circulating miRNAs and isomiRs have the potential to aid standard clinical testing in the early risk assessment of ASD and ID development in TSC patients.

OUP accepted manuscript

It is currently unknown how many RNA transcripts are able to induce degradation of microRNAs (miRNA) via the mechanism known as target-directed miRNA degradation (TDMD). We developed TDMDfinder, a computational pipeline that identifies ‘high confidence’ TDMD interactions in the Human and Mouse transcriptomes by combining sequence alignment and feature selection approaches. Our predictions suggested that TDMD is widespread, with potentially every miRNA controlled by endogenous targets. We experimentally tested 37 TDMDfinder predictions, of which 17 showed TDMD effects as measured by RT-qPCR and small RNA sequencing, linking the miR-17, miR-19, miR-30, miR-221, miR-26 and miR-23 families to novel endogenous TDMDs. In some cases, TDMD was found to affect different members of the same miRNA family selectively. Features like complementarity to the miRNA 3′ region, bulge size and hybridization energy appeared to be the main factors determining sensitivity. Computational analyses performed using the multiomic TCGA platform substantiated the involvement of many TDMD transcripts in human cancer and highlighted 36 highly significant interactions, suggesting TDMD as a new potential oncogenic mechanism. In conclusion, TDMDfinder provides the first inventory of bona fide human and mouse TDMDs. Available as a free webtool, TDMDfinder allows users to search for any TDMD interaction of interest by customizing its selection criteria.

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.

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

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

MiRNA post-transcriptional modification dynamics in T cell activation

T cell activation leads to extensive changes in the miRNA repertoire. Although overall miRNA expression decreases within a few hours of T cell activation, some individual miRNAs are specifically upregulated. Using next-generation sequencing, we assessed miRNA expression and post-transcriptional modification kinetics in human primary CD4+ T cells upon T cell receptor (TCR) or type I interferon stimulation. This analysis identified differential expression of multiple miRNAs not previously linked to T cell activation. Remarkably, upregulated miRNAs showed a higher frequency of 3' adenylation. TCR stimulation was followed by increased expression of RNA modifying enzymes and the RNA degrading enzymes Dis3L2 and Eri1. In the midst of this adverse environment, 3' adenylation may serve a protective function that could be exploited to improve miRNA stability for T cell-targeted therapy.

An integrated workflow for biomarker development using microRNAs in extracellular vesicles for cancer precision medicine

EV-miRNAs are microRNA (miRNA) molecules encapsulated in extracellular vesicles (EVs), which play crucial roles in tumor pathogenesis, progression, and metastasis. Recent studies about EV-miRNAs have gained novel insights into cancer biology and have demonstrated a great potential to develop novel liquid biopsy assays for various applications. Notably, compared to conventional liquid biomarkers, EV-miRNAs are more advantageous in representing host-cell molecular architecture and exhibiting higher stability and specificity. Despite various available techniques for EV-miRNA separation, concentration, profiling, and data analysis, a standardized approach for EV-miRNA biomarker development is yet lacking. In this review, we performed a substantial literature review and distilled an integrated workflow encompassing important steps for EV-miRNA biomarker development, including sample collection and EV isolation, EV-miRNA extraction and quantification, high-throughput data preprocessing, biomarker prioritization and model construction, functional analysis, as well as validation. With the rapid growth of "big data", we highlight the importance of efficient mining of high-throughput data for the discovery of EV-miRNA biomarkers and integrating multiple independent datasets for in silico and experimental validations to increase the robustness and reproducibility. Furthermore, as an efficient strategy in systems biology, network inference provides insights into the regulatory mechanisms and can be used to select functionally important EV-miRNAs to refine the biomarker candidates. Despite the encouraging development in the field, a number of challenges still hinder the clinical translation. We finally summarize several common challenges in various biomarker studies and discuss potential opportunities emerging in the related fields.

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.

IsomiR_Window: a system for analyzing small-RNA-seq data in an integrative and user-friendly manner

BACKGROUND

IsomiRs are miRNA variants that vary in length and/or sequence when compared to their canonical forms. These variants display differences in length and/or sequence, including additions or deletions of one or more nucleotides (nts) at the 5' and/or 3' end, internal editings or untemplated 3' end additions. Most available tools for small RNA-seq data analysis do not allow the identification of isomiRs and often require advanced knowledge of bioinformatics. To overcome this, we have developed IsomiR Window, a platform that supports the systematic identification, quantification and functional exploration of isomiR expression in small RNA-seq datasets, accessible to users with no computational skills.

METHODS

IsomiR Window enables the discovery of isomiRs and identification of all annotated non-coding RNAs in RNA-seq datasets from animals and plants. It comprises two main components: the IsomiR Window pipeline for data processing; and the IsomiR Window Browser interface. It integrates over ten third-party softwares for the analysis of small-RNA-seq data and holds a new algorithm that allows the detection of all possible types of isomiRs. These include 3' and 5'end isomiRs, 3' end tailings, isomiRs with single nucleotide polymorphisms (SNPs) or potential RNA editings, as well as all possible fuzzy combinations. IsomiR Window includes all required databases for analysis and annotation, and is freely distributed as a Linux virtual machine, including all required software.

RESULTS

IsomiR Window processes several datasets in an automated manner, without restrictions of input file size. It generates high quality interactive figures and tables which can be exported into different formats. The performance of isomiR detection and quantification was assessed using simulated small-RNA-seq data. For correctly mapped reads, it identified different types of isomiRs with high confidence and 100% accuracy. The analysis of a small RNA-seq data from Basal Cell Carcinomas (BCCs) using isomiR Window confirmed that miR-183-5p is up-regulated in Nodular BCCs, but revealed that this effect was predominantly due to a novel 5'end variant. This variant displays a different seed region motif and 1756 isoform-exclusive mRNA targets that are significantly associated with disease pathways, underscoring the biological relevance of isomiR-focused analysis. IsomiR Window is available at https://isomir.fc.ul.pt/ .

A longitudinal study defined circulating microRNAs as reliable biomarkers for disease prognosis and progression in ALS human patients

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease associated with motor neuron degeneration, muscle atrophy and paralysis. To date, multiple panels of biomarkers have been described in ALS patients and murine models. Nevertheless, none of them has sufficient specificity and thus the molecular signature for ALS prognosis and progression remains to be elucidated. Here we overcome this limitation through a longitudinal study, analyzing serum levels of circulating miRNAs, stable molecules that are recently used as promising biomarkers for many types of human disorders, in ALS patients during the progression of the pathology. We performed next-generation sequencing (NGS) analysis and absolute RT quantification of serum samples of ALS patients and healthy controls. The expression levels of five selected miRNAs were quantitatively analyzed during disease progression in each patient and we demonstrated that high levels of miR-206, miR-133a and miR-151a-5p can predict a slower clinical decline of patient functionality. In particular, we found that miR-206 and miR-151a-5p serum levels were significantly up-regulated at the mild stage of ALS pathology, to decrease in the following moderate and severe stages, whereas the expression levels of miR-133a and miR-199a-5p remained low throughout the course of the disease, showing a diagnostic significance in moderate and severe stages for miR-133a and in mild and terminal ones for miR-199a-5p. Moreover, we found that miR-423–3p and 151a-5p were significantly downregulated respectively in mild and terminal stages of the disease. These data suggest that these miRNAs represent potential prognostic markers for ALS disease.

Discovery and functional understanding of MiRNAs in molluscs: a genome-wide profiling approach

Small non-coding RNAs play a pivotal role in gene regulation, repression of transposable element and viral activity in various organisms. Among the various categories of these small non-coding RNAs, microRNAs (miRNAs) guide post-translational gene regulation in cellular development, proliferation, apoptosis, oncogenesis, and differentiation. Here, we performed a genome-wide computational prediction of miRNAs to improve the understanding of miRNA observation and function in molluscs. As an initial step, hundreds of conserved miRNAs were predicted in 35 species of molluscs through genome scanning. Afterwards, the miRNAs' population, isoforms, organization, and function were characterized in detail. Furthermore, the key miRNA biogenesis factors, including AGO2, DGCR8, DICER, DROSHA, TRABP2, RAN, and XPO5, were elucidated based on homologue sequence searching. We also summarized the miRNAs' function in biomineralization, immune and stress response, as well as growth and development in molluscs. Because miRNAs play a vital role in various lifeforms, this study will provide insight into miRNA biogenesis and function in molluscs, as well as other invertebrates.

MicroRNAs Regulating Autophagy in Neurodegeneration

Social and economic impacts of neurodegenerative diseases (NDs) become more prominent in our constantly aging population. Currently, due to the lack of knowledge about the aetiology of most NDs, only symptomatic treatment is available for patients. Hence, researchers and clinicians are in need of solid studies on pathological mechanisms of NDs. Autophagy promotes degradation of pathogenic proteins in NDs, while microRNAs post-transcriptionally regulate multiple signalling networks including autophagy. This chapter will critically discuss current research advancements in the area of microRNAs regulating autophagy in NDs. Moreover, we will introduce basic strategies and techniques used in microRNA research. Delineation of the mechanisms contributing to NDs will result in development of better approaches for their early diagnosis and effective treatment.

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

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

Computational methods for annotation of plant regulatory non-coding RNAs using RNA-seq

Plant transcriptome encompasses numerous endogenous, regulatory non-coding RNAs (ncRNAs) that play a major biological role in regulating key physiological mechanisms. While studies have shown that ncRNAs are extremely diverse and ubiquitous, the functions of the vast majority of ncRNAs are still unknown. With ever-increasing ncRNAs under study, it is essential to identify, categorize and annotate these ncRNAs on a genome-wide scale. The use of high-throughput RNA sequencing (RNA-seq) technologies provides a broader picture of the non-coding component of transcriptome, enabling the comprehensive identification and annotation of all major ncRNAs across samples. However, the detection of known and emerging class of ncRNAs from RNA-seq data demands complex computational methods owing to their unique as well as similar characteristics. Here, we discuss major plant endogenous, regulatory ncRNAs in an RNA sample followed by computational strategies applied to discover each class of ncRNAs using RNA-seq. We also provide a collection of relevant software packages and databases to present a comprehensive bioinformatics toolbox for plant ncRNA researchers. We assume that the discussions in this review will provide a rationale for the discovery of all major categories of plant ncRNAs.

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.

Factors Involved in miRNA Processing Change Its Expression Level during Imitation of Hypoxia in BeWo b30 Cells

One of the main complications of pregnancy and causes of maternal and perinatal mortality is preeclampsia. The pathophysiology of preeclampsia is associated with the development of placenta and fetal hypoxia and secretion of a number of effective molecules. The human choriocarcinoma cell line BeWo b30 is often used as a model of the placental barrier. It was shown that oxyquinoline derivatives can mimic hypoxia by suppressing HIF-prolyl hydroxylases and the accumulation of HIF-1α. This effect also leads to a change in the expression of microRNAs and their target genes. However, with hypoxia in cells, not only the level of individual miRNAs but also the ratio of miRNA isoforms (isomiRs) can change, presumably due to inaccuracies in the work of the Drosha and Dicer enzymes. In this work, we showed a change in the expression of the factors involved in the maturation of miRNAs when simulating hypoxia in BeWo b30 cells with an oxyquinoline derivative, which may be one of the causes for the change in the ratio of miRNA isoforms.

Computational tools for plant small RNA detection and categorization

Small RNAs (sRNAs) are important short-length molecules with regulatory functions essential for plant development and plasticity. High-throughput sequencing of total sRNA populations has revealed that the largest share of sRNA remains uncategorized. To better understand the role of sRNA-mediated cellular regulation, it is necessary to create accurate and comprehensive catalogues of sRNA and their sequence features, a task that currently relies on nontrivial bioinformatic approaches. Although a large number of computational tools have been developed to predict features of sRNA sequences, these tools are mostly dedicated to microRNAs and none integrates the functionalities necessary to describe units from all sRNA pathways thus far discovered in plants. Here, we review the different classes of sRNA found in plants and describe available bioinformatics tools that can help in their detection and categorization.

Single nucleotide alterations in MicroRNAs and human cancer-A not fully explored field

MicroRNAs are ~20 nt long small noncoding RNAs that are processed from stem-looped precursors and function mainly as posttranscriptional regulators of protein coding genes through binding to 3'-untranslated regions of messenger RNAs to inhibit the translation or cause RNA degradation. It is predicted microRNAs could regulate up to half of all human genes and are proved to play important roles in human diseases including cancer. They bind to target mRNAs based on complementary binding which is dominated by the so-called "seed" region which are the 5' 2-8 bases of the microRNA. Due to the small size in nature, even a single nucleotide variation in the precursor region especially those located in the seed regions could show big influence. Here, I summarized and reviewed the current knowledge of these single nucleotide alterations in microRNAs in human cancer including (i) common SNPs in the precursor region, (ii) isomiRs, (iii) somatic mutations of microRNAs. Briefly, this is an underexploited field and clearly, warrants further studies to reveal their biological and clinical significances. I believe they will be key to advancing personalized medicine.

Deciphering miRNAs' Action through miRNA Editing

MicroRNAs (miRNAs) are small non-coding RNAs with the capability of modulating gene expression at the post-transcriptional level either by inhibiting messenger RNA (mRNA) translation or by promoting mRNA degradation. The outcome of a myriad of physiological processes and pathologies, including cancer, cardiovascular and metabolic diseases, relies highly on miRNAs. However, deciphering the precise roles of specific miRNAs in these pathophysiological contexts is challenging due to the high levels of complexity of their actions. Indeed, regulation of mRNA expression by miRNAs is frequently cell/organ specific; highly dependent on the stress and metabolic status of the organism; and often poorly correlated with miRNA expression levels. Such biological features of miRNAs suggest that various regulatory mechanisms control not only their expression, but also their activity and/or bioavailability. Several mechanisms have been described to modulate miRNA action, including genetic polymorphisms, methylation of miRNA promoters, asymmetric miRNA strand selection, interactions with RNA-binding proteins (RBPs) or other coding/non-coding RNAs. Moreover, nucleotide modifications (A-to-I or C-to-U) within the miRNA sequences at different stages of their maturation are also critical for their functionality. This regulatory mechanism called "RNA editing" involves specific enzymes of the adenosine/cytidine deaminase family, which trigger single nucleotide changes in primary miRNAs. These nucleotide modifications greatly influence a miRNA's stability, maturation and activity by changing its specificity towards target mRNAs. Understanding how editing events impact miRNA's ability to regulate stress responses in cells and organs, or the development of specific pathologies, e.g., metabolic diseases or cancer, should not only deepen our knowledge of molecular mechanisms underlying complex diseases, but can also facilitate the design of new therapeutic approaches based on miRNA targeting. Herein, we will discuss the current knowledge on miRNA editing and how this mechanism regulates miRNA biogenesis and activity.

An Approach to Identify Individual Functional Single Nucleotide Polymorphisms and Isoform MicroRNAs

MicroRNAs (miRNAs) and single nucleotide polymorphisms (SNPs) play important roles in disease risk and development, especially cancer. Importantly, when SNPs are located in pre-miRNAs, they affect their splicing mechanism and change the function of miRNAs. To improve disease risk assessment, we propose an approach and developed a software tool, IsomiR_Find, to identify disease/phenotype-related SNPs and isomiRs in individuals. Our approach is based on the individual's samples, with SNP information extracted from the 1000 Genomes Project. SNPs were mapped to pre-miRNAs based on whole-genome coordinates and then SNP-pre-miRNA sequences were constructed. Moreover, we developed matpred2, a software tool to identify the four splicing sites of mature miRNAs. Using matpred2, we identified isomiRs and then verified them by searching within individual miRNA sequencing data. Our approach yielded biomarkers for biological experiments, mined functions of miRNAs and SNPs, improved disease risk assessment, and provided a way to achieve individualized precision medicine.

Sleeping beauty genetic screen identifies miR-23b::BTBD7 gene interaction as crucial for colorectal cancer metastasis

Background

Metastatic colorectal cancer (CRC) remains a deadly disease. Identifying locally advanced CRC patients with high risk of developing metastasis and improving outcome of metastatic CRC patients require discovering master regulators of metastasis. In this context, the non-coding part of the human genome is still largely unexplored.

Methods

To interrogate the non-coding part of the human genome and disclose regulators of CRC metastasis, we combined a transposon-based forward genetic screen with a novel in vitro assay, which forces cells to grow deprived of cell-substrate and cell-cell contacts (i.e. forced single cell suspension assay - fSCS).

Findings

We proved that fSCS selects CRC cells with mesenchymal and pro-metastatic traits. Moreover, we found that the transposon insertions conferred CRC cells resistance to fSCS and thus metastatic advantage. Among the retrieved transposon insertions, we demonstrated that the one located in the 3′UTR of BTBD7 disrupts miR-23b::BTBD7 interaction and contributes to pro-metastatic traits. In addition, miR-23b and BTBD7 correlate with CRC metastasis both in preclinical experiments and in clinical samples.

Interpretation

fSCS is a simple and scalable in vitro assay to investigate pro-metastatic traits and transposon-based genetic screens can interrogate the non-coding part of the human genome (e.g. miRNA::target interactions). Finally, both Btbd7 and miR-23b represent promising prognostic biomarkers and therapeutic targets in CRC.

Fund

This work was supported by Marie Curie Actions (CIG n. 303877) and Friuli Venezia Giulia region (Grant Agreement n°245574), Italian Association for Cancer Research (AIRC, MFAG n°13589), Italian Ministry of Health (GR-2010-2319387 and PE-2016-02361040) and 5x1000 to CRO Aviano.

Computational Resources for Prediction and Analysis of Functional miRNA and Their Targetome

microRNAs are evolutionarily conserved, endogenously produced, noncoding RNAs (ncRNAs) of approximately 19-24 nucleotides (nts) in length known to exhibit gene silencing of complementary target sequence. Their deregulated expression is reported in various disease conditions and thus has therapeutic implications. In the last decade, various computational resources are published in this field. In this chapter, we have reviewed bioinformatics resources, i.e., miRNA-centered databases, algorithms, and tools to predict miRNA targets. First section has enlisted more than 75 databases, which mainly covers information regarding miRNA registries, targets, disease associations, differential expression, interactions with other noncoding RNAs, and all-in-one resources. In the algorithms section, we have compiled about 140 algorithms from eight subcategories, viz. for the prediction of precursor (pre-) and mature miRNAs. These algorithms are developed on various sequence, structure, and thermodynamic based features incorporated into different machine learning techniques (MLTs). In addition, computational identification of miRNAs from high-throughput next generation sequencing (NGS) data and their variants, viz. isomiRs, differential expression, miR-SNPs, and functional annotation, are discussed. Prediction and analysis of miRNAs and their associated targets are also evaluated under miR-targets section providing knowledge regarding novel miRNA targets and complex host-pathogen interactions. In conclusion, we have provided comprehensive review of in silico resources published in miRNA research to help scientific community be updated and choose the appropriate tool according to their needs.

Diff isomiRs: Large-scale detection of differential isomiRs for understanding non-coding regulated stress omics in plants

Plants have an amazing ability to cope with wide variety of stresses by regulating the expression of genes and thus by altering the physiological status. In the past few years, canonical microRNA variants (isomiRs) have been shown to play pivotal roles by acting as regulators of the transcriptional machinery. In the present research, we present Diff isomiRs, a web-based exploratory repository of differential isomiRs across 16 sequenced plant species representing a total of 433 datasets across 21 different stresses and 158 experimental states. Diff isomiRs provides the high-throughput detection of differential isomiRs using mapping-based and model-based differential analysis revealing a total of 16,157 and 2,028 differential isomiRs, respectively. Easy-to-use and web-based exploration of differential isomiRs provides several features such as browsing of the differential isomiRs according to stress or species, as well as association of the differential isomiRs to targets and plant endogenous target mimics (PeTMs). Diff isomiRs also provides the relationship between the canonical miRNAs, isomiRs and the miRNA-target interactions. This is the first web-based large-scale repository for browsing differential isomiRs and will facilitate better understanding of the regulatory role of the isomiRs with respect to the canonical microRNAs. Diff isomiRs can be accessed at: www.mcr.org.in/diffisomirs.

An Introduction to Methods for Discovery and Functional Analysis of MicroRNAs in Plants

MicroRNAs play important roles in posttranscriptional regulation of plant development, metabolism, and abiotic stress responses. The recent generation of massive amounts of small RNA sequence data, along with development of bioinformatic tools to identify miRNAs and their mRNA targets, has led to an explosion of newly identified putative miRNAs in plants. Genome editing techniques like CRISPR-Cas9 will allow us to study the biological role of these potential novel miRNAs by efficiently targeting both the miRNA and its mRNA target. In this chapter, we review bioinformatic tools and experimental methods for the identification and functional characterization of miRNAs and their target mRNAs in plants.

An integrative piRNA analysis of mouse gametes and zygotes reveals new potential origins and gene regulatory roles

Piwi-interacting RNAs (piRNAs) are a subclass of the small non-coding RNAs (sncRNAs). Their main reported function was to exert control over transposable elements (TEs) in mammalian germline. In this study undertaking a deeper bioinformatics analysis of piRNAs present in mouse oocytes, sperm cells and zygotes, we first elaborated a new piRNA database based on sequences identified as piRNAs by immunoprecipitation with PIWI proteins. Our bioinformatics analysis revealed that, at least in gametes and zygotes, piRNAs could encompass multifunctional cell-dependent regulatory molecules. Indeed, genome analysis of the piRNA mapping density (reads/kb) evidenced in all samples an enrichment of intron-derived piRNAs. Further, piRNA population was classified into sequences not associated to TEs or repeats (NRapiRNAs) and associated to repetitive genome elements (RapiRNAs). In oocytes most of the NRapiRNAs mapped to the 5′UTRs of coding mRNAs, while higher proportion of NRapiRNAs was detected in sperm cells associated to the 3′UTRs of mRNAs. This piRNA complementarity to mRNA UTRs suggests key post-transcriptional regulatory roles over mRNAs such as those encoding MHC genes. In addition, a striking association of RapiRNA with long non-coding RNAs (lncRNAs) was identified. piRNAs associated with relevant lncRNAs such as: Rab26os and GAS5 and key mRNAs, were particularly assessed.

Identifying and characterizing functional 3' nucleotide addition in the miRNA pathway

Over the past decade, modifications to microRNAs (miRNAs) via 3' end nucleotide addition have gone from a deep-sequencing curiosity to experimentally confirmed drivers of a range of regulatory activities. Here we overview the methods that have been deployed by researchers seeking to untangle these diverse functional roles and include characterizing not only the nucleotidyl transferases catalyzing the additions but also the nucleotides being added, and the timing of their addition during the miRNA pathway. These methods and their further development are key to clarifying the diverse and sometimes contradictory functional findings presently attributed to these nucleotide additions.

Endogenous transcripts control miRNA levels and activity in mammalian cells by target-directed miRNA degradation

Little is known about miRNA decay. A target-directed miRNA degradation mechanism (TDMD) has been suggested, but further investigation on endogenous targets is necessary. Here, we identify hundreds of targets eligible for TDMD and show that an endogenous RNA (Serpine1) controls the degradation of two miRNAs (miR-30b-5p and miR-30c-5p) in mouse fibroblasts. In our study, TDMD occurs when the target is expressed at relatively low levels, similar in range to those of its miRNAs (100-200 copies per cell), and becomes more effective at high target:miRNA ratios (>10:1). We employ CRISPR/Cas9 to delete the miR-30 responsive element within Serpine1 3'UTR and interfere with TDMD. TDMD suppression increases miR-30b/c levels and boosts their activity towards other targets, modulating gene expression and cellular phenotypes (i.e., cell cycle re-entry and apoptosis). In conclusion, a sophisticated regulatory layer of miRNA and gene expression mediated by specific endogenous targets exists in mammalian cells.

3' RNA Uridylation in Epitranscriptomics, Gene Regulation, and Disease

Emerging evidence implicates a wide range of post-transcriptional RNA modifications that play crucial roles in fundamental biological processes including regulating gene expression. Collectively, they are known as epitranscriptomics. Recent studies implicate 3' RNA uridylation, the non-templated addition of uridine(s) to the terminal end of RNA, as a key player in epitranscriptomics. In this review, we describe the functional roles and significance of 3' terminal RNA uridylation that has diverse functions in regulating both mRNAs and non-coding RNAs. In mammals, three Terminal Uridylyl Transferases (TUTases) are primarily responsible for 3' RNA uridylation. These enzymes are also referred to as polyU polymerases. TUTase 1 (TUT1) is implicated in U6 snRNA maturation via uridylation. The TUTases TUT4 and/or TUT7 are the predominant mediators of all other cellular uridylation. Terminal uridylation promotes turnover for many polyadenylated mRNAs, replication-dependent histone mRNAs that lack polyA-tails, and aberrant structured noncoding RNAs. In addition, uridylation regulates biogenesis of a subset of microRNAs and generates isomiRs, sequent variant microRNAs that have altered function in specific cases. For example, the RNA binding protein and proto-oncogene LIN28A and TUT4 work together to polyuridylate pre-let-7, thereby blocking biogenesis and function of the tumor suppressor let-7 microRNA family. In contrast, monouridylation of Group II pre-miRNAs creates an optimal 3' overhang that promotes recognition and subsequent cleavage by the Dicer-TRBP complex that then yields the mature microRNA. Also, uridylation may play a role in non-canonical microRNA biogenesis. The overall significance of 3' RNA uridylation is discussed with an emphasis on mammalian development, gene regulation, and disease, including cancer and Perlman syndrome. We also introduce recent changes to the HUGO-approved gene names for multiple terminal nucleotidyl transferases that affects in part TUTase nomenclature (TUT1/TENT1, TENT2/PAPD4/GLD2, TUT4/ZCCHC11/TENT3A, TUT7/ZCCHC6/TENT3B, TENT4A/PAPD7, TENT4B/PAPD5, TENT5A/FAM46A, TENT5B/FAM46B, TENT5C/FAM46C, TENT5D/FAM46D, MTPAP/TENT6/PAPD1).

A Comprehensive Approach to Sequence-oriented IsomiR annotation (CASMIR): demonstration with IsomiR profiling in colorectal neoplasia

Background

MicroRNA (miRNA) profiling is an important step in studying biological associations and identifying marker candidates. miRNA exists in isoforms, called isomiRs, which may exhibit distinct properties. With conventional profiling methods, limitations in assay and analysis platforms may compromise isomiR interrogation.

Results

We introduce a comprehensive approach to sequence-oriented isomiR annotation (CASMIR) to allow unbiased identification of global isomiRs from small RNA sequencing data. In this approach, small RNA reads are maintained as independent sequences instead of being summarized under miRNA names. IsomiR features are identified through step-wise local alignment against canonical forms and precursor sequences. Through customizing the reference database, CASMIR is applicable to isomiR annotation across species. To demonstrate its application, we investigated isomiR profiles in normal and neoplastic human colorectal epithelia. We also ran miRDeep2, a popular miRNA analysis algorithm to validate isomiRs annotated by CASMIR. With CASMIR, specific and biologically relevant isomiR patterns could be identified. We note that specific isomiRs are often more abundant than their canonical forms. We identify isomiRs that are commonly up-regulated in both colorectal cancer and advanced adenoma, and illustrate advantages in targeting isomiRs as potential biomarkers over canonical forms.

Conclusions

Studying miRNAs at the isomiR level could reveal new insight into miRNA biology and inform assay design for specific isomiRs. CASMIR facilitates comprehensive annotation of isomiR features in small RNA sequencing data for isomiR profiling and differential expression analysis.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4794-7) contains supplementary material, which is available to authorized users.

MicroRNA profiling in MDA-MB-231 human breast cancer cell exposed to the Phaleria macrocarpa (Boerl.) fruit ethyl acetate fraction (PMEAF) through IIlumina Hi-Seq technologies and various in silico bioinformatics tools

ETHNOPHARMACOLOGICAL RELEVANCE

Phaleria macrocarpa (Scheff) Boerl, is a famous traditional medicinal plant which exhibited cytotoxicity against various cancerous cells. Traditionally, P. macrocarpa has been used to control cancer, impotency, hemorrhoids, diabetes mellitus, allergies, liver and heart disease, kidney disorders, blood diseases, acne, stroke, migraine, and various skin diseases.

AIM OF THE STUDY

Recent studies have demonstrated a potent anticancer potential of P. macrocarpa, especially against HeLa cell. The objective of this study was to investigate the regulation of miRNAs on MDA-MB-231 treated with P. macrocarpa ethyl acetate fraction (PMEAF).

MATERIALS AND METHODS

The regulation of miRNAs on MDA-MB-231 cells treated with PMEAF was studied through IIlumina, Hi-Seq. 2000 platform of Next Generation Sequencing (NGS) and various in silico bioinformatics tools.

RESULTS

The PMEAF treatment against MDA-MB-231 cells identified 10 upregulated and 10 downregulated miRNAs. A set of 606 target genes of 10 upregulated miRNAs and 517 target genes of 10 downregulated miRNAs were predicted based on computational and validated databases by using miRGate DB Query. Meanwhile, results from DAVID Bioinformatics Resources 6.8 specified the functional annotation of the upregulated miRNAs involvement in cancer pathway by suppressing the oncogenes and downregulating miRNAs by expressing the tumour suppressor genes in the regulation of apoptosis pathway.

CONCLUSION

In conclusion, the results of this study proved that PMEAF is a promising anticancer agent with high cytotoxicity against MDA-MB-231 breast cancer cells and it induced apoptotic cell death mechanism through the regulation of miRNAs. PMEAF might be the best candidate for developing more potent anticancer drugs or chemo preventive supplements.

Mapping the microRNA Expression Profiles in Glyoxalase Over-expressing Salinity Tolerant Rice

In the recent years, glyoxalase pathway has been an active area of research in both human and plants. This pathway is reported to confer stress tolerance in plants, by modulating the glutathione homeostasis to achieve detoxification of a potent cytotoxic and mutagenic compound, methylglyoxal. The microRNAs (miRNAs) are also reported to play significant role in stress tolerance for plants. However, the cross-talk of miRNAs with the metabolism regulated by glyoxalase in the salinity-tolerance is unexplored. We therefore investigated whether expression profiles of miRNAs are altered in response to glyoxalase overexpression, and if any of these are also responsible for modulating the stress responses of plants. In this study, the Next Generation Sequencing (NGS) was employed to profile miRNA expression levels from glyoxalase overexpressing transgenic lines. The associated targets of differentially expressed miRNAs were predicted and their functional annotation was carried out using Gene Ontology (GO) and KEGG Orthology (KO), which showed their involvement in several crucial biological pathways. The analysis of NGS datasets also identified other isoforms or isomiRs of selected miRNAs, which may have an active role in developing tolerance against salt stress. Different aspects of miRNA modifications were also studied in glyoxalase overexpressing lines.

Uncovering the Stability of Mature miRNAs by 4-Thio-Uridine Metabolic Labeling

MicroRNAs (miRNAs) are an evolutionary conserved class of short, single-stranded noncoding RNAs (<18-22 nt in length) that act in posttranscriptional regulation of gene expression in higher eukaryotes. The abundance of a miRNA is a key feature in control of its activity and, therefore, a number of mechanisms finely regulate miRNA levels, acting at both transcriptional and posttranscriptional level. Recent evidences, including our research, highlighted the role of miRNA decay as a mechanism controlling the miRNA pool. We describe in this chapter an optimized methodology to determine miRNA degradation rates in mammalian cells. Our approach is based on metabolic pulse labeling with 4-thiouridine (4sU), a uridine analog that is incorporated in nascent RNA and allows thiol-specific biotinylation and selective pull-down of labeled RNA. In particular, given the long average half-life and the complex biogenetic process of miRNAs, we developed a "pulse-chase" protocol where 4sU is removed from the medium after a long labeling period (2-3 h pulse), and labeled RNA is purified at different time points to measure the decay of labeled molecules. By combining the 4sU-based "pulse-chase" approach with high-throughput small RNA sequencing (sRNAseq), it is possible to quantify at genome-wide level miRNA degradation rates.

Plant IsomiR Atlas: Large Scale Detection, Profiling, and Target Repertoire of IsomiRs in Plants

microRNAs (miRNAs) play an important role as key regulators controlling the post-transcriptional events in plants across development, abiotic and biotic stress, tissue polarity and also in defining the evolutionary basis of the origin of the post-transcriptional machinery. Identifying patterns of regulated and co-regulated small RNAs, in particular miRNAs and their sequence variants with the availability of next generation sequencing approaches has widely demonstrated the role of miRNAs and their temporal regulation in maintaining plant development and their response to stress conditions. Although the role of canonical miRNAs has been widely explored and functional diversity is revealed, those works for isomiRs are still limited and urgent to be carried out across plants. This relative lack of information with respect to isomiRs might be attributed to the non-availability of large-scale detection of isomiRs across wide plant species. In the present research, we addressed this by developing Plant isomiR Atlas, which provides large-scale detection of isomiRs across 23 plant species utilizing 677 smallRNAs datasets and reveals a total of 98,374 templated and non-templated isomiRs from 6,167 precursors. Plant isomiR Atlas provides several visualization features such as species specific isomiRs, isomiRs and canonical miRNAs overlap, terminal modification classifications, target identification using psRNATarget and TargetFinder and also canonical miRNAs:target interactions. Plant isomiR Atlas will play a key role in understanding the regulatory nature of miRNAome and will accelerate to understand the functional role of isomiRs. Plant isomiR Atlas is available at www.mcr.org.in/isomir. One Sentence Summary  Plant isomiR Atlas will play a key role in understanding the regulatory nature of miRNAome and will accelerate the understanding and diversity of functional targets of plants isomiRs.

Evaluation of high-throughput isomiR identification tools: illuminating the early isomiRome of Tribolium castaneum

Background

MicroRNAs carry out post-transcriptional gene regulation in animals by binding to the 3' untranslated regions of mRNAs, causing their degradation or translational repression. MicroRNAs influence many biological functions, and dysregulation can therefore disrupt development or even cause death. High-throughput sequencing and the mining of animal small RNA data has shown that microRNA genes can yield differentially expressed isoforms, known as isomiRs. Such isoforms are particularly relevant during early development, and the extension or truncation of the 5' end can change the profile of mRNA targets compared to the original mature sequence. We used the publicly available small RNA dataset of the model beetle Tribolium castaneum to create the first comparative isomiRome of early developmental stages in this species. Standard microRNA analysis software does not specifically account for isomiRs. We therefore carried out the first comparative evaluation of the specialized tools isomiRID, isomiR-SEA and miraligner, which can be downloaded for local use and can handle next generation sequencing data.

Results

We compared the performance of isomiRID, isomiR-SEA and miraligner using simulated Illumina HiSeq2000 and MiSeq data to test the impact of technical errors. We also created artificial microRNA isoforms to determine the effect of biological variants on the performance of each algorithm. We found that isomiRID achieved the best true positive rate among the three algorithms, but only accounted for one mutation at a time. In contrast, miraligner reported all variations simultaneously but with 78% sensitivity, yielding isomiRs with 3' or 5' deletions. Finally, isomiR-SEA achieved a sensitivity of 25–33% when the seed region was mutated or partly deleted, but was the only tool that could accommodate more than one mismatch. Using the best tool, we performed a complete isomiRome analysis of the early developmental stages of T. castaneum.

Conclusions

Our findings will help researchers to select the most suitable isomiR analysis tools for their experiments. We confirmed the dynamic expression of 3′ non-template isomiRs and expanded the isomiRome by all known isomiR modifications during the early development of T. castaneum.

Electronic supplementary material

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

Chronic low-dose exposure to a mixture of environmental endocrine disruptors induces microRNAs/isomiRs deregulation in mouse concomitant with intratesticular estradiol reduction

Humans are environmentally exposed not only to single endocrine-disrupting chemicals (EDCs) but to mixtures that affect their reproductive health. In reproductive tissues, microRNAs (miRNAs) are emerging as key targets of EDCs. Here, we analysed changes in the testis "miRNome" (and their biogenesis mechanism) in chronically exposed adult mice to a cocktail of five EDCs containing 0.3 mg/kg-body weight (BW)/day of each phthalate (DEHP, DBP, BBP) and 0.05 mg/kg-BW/day of each alkylphenol (NP, OP), from conception to adulthood. The testis "miRNome" was characterised using next-generation sequencing (NGS). Expression levels of genes involved in miRNA biogenesis were measured by RT-qPCR, as well as several physiological and cytological parameters. We found two up-regulated, and eight down-regulated miRNAs and thirty-six differentially expressed isomiRs along with an over-expression of Drosha, Adar and Zcchc11. A significant decrease of intratesticular estradiol but not testosterone was detected. Functional analysis showed altered spermatogenesis, germ cell apoptosis and negative correlation of miR-18a-5p with Nr1h2 involved in the deregulation of the steroidogenesis pathway. Here, we present the first association between miRNA/isomiRs deregulation, their mechanisms of biogenesis and histopathological and hormonal alterations in testes of adult mice exposed to a mixture of low-dose EDCs, which can play a role in male infertility.

Serum microRNA signatures as "liquid biopsies" for interrogating hepatotoxic mechanisms and liver pathogenesis in human

MicroRNAs (miRNAs) released into the peripheral circulation upon cellular injury have shown a promise as a new class of tissue-specific biomarkers. We were first to demonstrate that next-generation sequencing analysis of serum from human subjects with acetaminophen-induced liver injury revealed a specific signature of circulating miRNAs. We consequently hypothesized that different types of hepatic liver impairments might feature distinct signatures of circulating miRNAs and that this approach might be useful as minimally invasive diagnostic "liquid biopsies" enabling the interrogation of underlying molecular mechanisms of injury in distant tissues. Therefore we examined serum circulating miRNAs in a total of 72 serum samples from a group of 53 subjects that included patients with accidental acetaminophen overdose, hepatitis B infection, liver cirrhosis and type 2 diabetes as well as gender- and age-matched healthy subjects with no evidence of liver disease. The miRNA signatures were identified using next-generation sequencing that provided analysis for the whole miRNome, including miRNA isoforms. Compared to the healthy subjects, a total of 179 miRNAs showed altered serum levels across the diseased subjects. Although many subjects have elevated alanine aminotransferase suggesting liver impairments, we identified distinct miRNA signatures for different impairments with minimum overlap. Furthermore, the bioinformatics analysis of miRNA signatures revealed relevant molecular pathways associated with the mechanisms of toxicity and or pathogenesis of disease. Interestingly, the high proportion of miRNA isoforms present in the respective signatures indicated a new level of complexity in cellular response to stress or disease. Our study demonstrates for the first time that signatures of circulating miRNAs show specificity for liver injury phenotypes and, once validated, might become useful for diagnosis of organ pathologies as "liquid biopsies".

MicroRNA-independent functions of DGCR8 are essential for neocortical development and TBR1 expression

Recent evidence indicates that the miRNA biogenesis factors DROSHA, DGCR8, and DICER exert non-overlapping functions, and have also roles in miRNA-independent regulatory mechanisms. However, it is currently unknown whether miRNA-independent functions of DGCR8 play any role in the maintenance of neuronal progenitors and during corticogenesis. Here, by phenotypic comparison of cortices from conditional Dgcr8 and Dicer knockout mice, we show that Dgcr8 deletion, in contrast to Dicer depletion, leads to premature differentiation of neural progenitor cells and overproduction of TBR1-positive neurons. Remarkably, depletion of miRNAs upon DCGR8 loss is reduced compared to DICER loss, indicating that these phenotypic differences are mediated by miRNA-independent functions of DGCR8. We show that Dgcr8 mutations induce an earlier and stronger phenotype in the developing nervous system compared to Dicer mutants and that miRNA-independent functions of DGCR8 are critical for corticogenesis. Finally, our data also suggest that the Microprocessor complex, with DROSHA and DGCR8 as core components, directly regulates the Tbr1 transcript, containing evolutionarily conserved hairpins that resemble miRNA precursors, independently of miRNAs.

DeAnnIso: a tool for online detection and annotation of isomiRs from small RNA sequencing data

Small RNA (sRNA) Sequencing technology has revealed that microRNAs (miRNAs) are capable of exhibiting frequent variations from their canonical sequences, generating multiple variants: the isoforms of miRNAs (isomiRs). However, integrated tool to precisely detect and systematically annotate isomiRs from sRNA sequencing data is still in great demand. Here, we present an online tool, DeAnnIso (Detection and Annotation of IsomiRs from sRNA sequencing data). DeAnnIso can detect all the isomiRs in an uploaded sample, and can extract the differentially expressing isomiRs from paired or multiple samples. Once the isomiRs detection is accomplished, detailed annotation information, including isomiRs expression, isomiRs classification, SNPs in miRNAs and tissue specific isomiR expression are provided to users. Furthermore, DeAnnIso provides a comprehensive module of target analysis and enrichment analysis for the selected isomiRs. Taken together, DeAnnIso is convenient for users to screen for isomiRs of their interest and useful for further functional studies. The server is implemented in PHP + Perl + R and available to all users for free at: http://mcg.ustc.edu.cn/bsc/deanniso/ and http://mcg2.ustc.edu.cn/bsc/deanniso/.

miR-isomiRExp: a web-server for the analysis of expression of miRNA at the miRNA/isomiR levels

MicroRNA (miRNA) locus has been found that can generate a series of varied isomiR sequences. Most studies always focus on determining miRNA level, however, the canonical miRNA sequence is only a specific member in the multiple isomiRs. Some studies have shown that isomiR sequences play versatile roles in biological progress, and the analysis and research should be simultaneously performed at the miRNA/isomiR levels. Based on the biological characteristics of miRNA and isomiR, we developed miR-isomiRExp to analyze expression pattern of miRNA at the miRNA/isomiR levels, provide insights into tracking miRNA/isomiR maturation and processing mechanisms, and reveal functional characteristics of miRNA/isomiR. Simultaneously, we also performed expression analysis of specific human diseases using public small RNA sequencing datasets based on the analysis platform, which may help in surveying the potential deregulated miRNA/isomiR expression profiles, especially sequence and function-related isomiRs for further interaction analysis and study. The miR-isomiRExp platform provides miRNA/isomiR expression patterns and more information to study deregulated miRNA loci and detailed isomiR sequences. This comprehensive analysis will enrich experimental miRNA studies. miR-isomiRExp is available at http://server.malab.cn/miRisomiRExp/.

Degradation dynamics of microRNAs revealed by a novel pulse-chase approach

The regulation of miRNAs is critical to the definition of cell identity and behavior in normal physiology and disease. To date, the dynamics of miRNA degradation and the mechanisms involved in remain largely obscure, in particular, in higher organisms. Here, we developed a pulse-chase approach based on metabolic RNA labeling to calculate miRNA decay rates at genome-wide scale in mammalian cells. Our analysis revealed heterogeneous miRNA half-lives, with many species behaving as stable molecules (T_1/2 > 24 h), while others, including passenger miRNAs and a number (25/129) of guide miRNAs, are quickly turned over (T_1/2 = 4–14 h). Decay rates were coupled with other features, including genomic organization, transcription rates, structural heterogeneity (isomiRs), and target abundance, measured through quantitative experimental approaches. This comprehensive analysis highlighted functional mechanisms that mediate miRNA degradation, as well as the importance of decay dynamics in the regulation of the miRNA pool under both steady-state conditions and during cell transitions.

miRNAs regulated by estrogens, tamoxifen, and endocrine disruptors and their downstream gene targets

MicroRNAs (miRNAs) are short (22 nucleotides), single-stranded, non-coding RNAs that form complimentary base-pairs with the 3' untranslated region of target mRNAs within the RNA-induced silencing complex (RISC) and block translation and/or stimulate mRNA transcript degradation. The non-coding miRBase (release 21, June 2014) reports that human genome contains ∼ 2588 mature miRNAs which regulate ∼ 60% of human protein-coding mRNAs. Dysregulation of miRNA expression has been implicated in estrogen-related diseases including breast cancer and endometrial cancer. The mechanism for estrogen regulation of miRNA expression and the role of estrogen-regulated miRNAs in normal homeostasis, reproduction, lactation, and in cancer is an area of great research and clinical interest. Estrogens regulate miRNA transcription through estrogen receptors α and β in a tissue-specific and cell-dependent manner. This review focuses primarily on the regulation of miRNA expression by ligand-activated ERs and their bona fide gene targets and includes miRNA regulation by tamoxifen and endocrine disrupting chemicals (EDCs) in breast cancer and cell lines.

New global analysis of the microRNA transcriptome of primary tumors and lymph node metastases of papillary thyroid cancer

Background

Papillary Thyroid Cancer (PTC) is the most prevalent type of endocrine cancer. Its incidence has rapidly increased in recent decades but little is known regarding its complete microRNA transcriptome (miRNome). In addition, there is a need for molecular biomarkers allowing improved PTC diagnosis.

Methods

We performed small RNA deep-sequencing of 3 PTC, their matching normal tissues and lymph node metastases (LNM). We designed a new bioinformatics framework to handle each aspect of the miRNome: whole expression profiles, isomiRs distribution, non-templated additions distributions, RNA-editing or mutation. Results were validated experimentally by qRT-PCR on normal samples, tumors and LNM from 14 independent patients and in silico using the dataset from The Cancer Genome Atlas (small RNA deepsequencing of 59 normal samples, 495 PTC, and 8 LNM).

Results

We performed small RNA deep-sequencing of 3 PTC, their matching normal tissues and lymph node metastases (LNM). We designed a new bioinformatics framework to handle each aspect of the miRNome: whole expression profiles, isomiRs distribution, non-templated additions distributions, RNA-editing or mutation. Results were validated experimentally by qRT-PCR on normal samples, tumors and LNM from 14 independent patients and in silico using the dataset from The Cancer Genome Atlas (small RNA deep-sequencing of 59 normal samples, 495 PTC, and 8 LNM). We confirmed already described up-regulations of microRNAs in PTC, such as miR-146b-5p or miR-222-3p, but we also identified down-regulated microRNAs, such as miR-7-5p or miR-30c-2-3p. We showed that these down-regulations are linked to the tumorigenesis process of thyrocytes. We selected the 14 most down-regulated microRNAs in PTC and we showed that they are potential biomarkers of PTC samples. Nevertheless, they can distinguish histological classical variants and follicular variants of PTC in the TCGA dataset. In addition, 12 of the 14 down-regulated microRNAs are significantly less expressed in aggressive PTC compared to non-aggressive PTC. We showed that the associated aggressive expression profile is mainly due to the presence of the BRAF V600E mutation. In general, primary tumors and LNM presented similar microRNA expression profiles but specific variations like the down-regulation of miR-7-2-3p and miR-30c-2-3p in LNM were observed. Investigations of the 5p-to-3p arm expression ratios, non-templated additions or isomiRs distributions revealed no major implication in PTC tumorigenesis process or LNM appearance.

Conclusions

Our results showed that down-regulated microRNAs can be used as new potential common biomarkers of PTC and to distinguish main subtypes of PTC. MicroRNA expressions can be linked to the development of LNM of PTC. The bioinformatics framework that we have developed can be used as a starting point for the global analysis of any microRNA deep-sequencing data in an unbiased way.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2082-3) contains supplementary material, which is available to authorized users.

miRMOD: a tool for identification and analysis of 5' and 3' miRNA modifications in Next Generation Sequencing small RNA data

In the past decade, the microRNAs (miRNAs) have emerged to be important regulators of gene expression across various species. Several studies have confirmed different types of post-transcriptional modifications at terminal ends of miRNAs. The reports indicate that miRNA modifications are conserved and functionally significant as it may affect miRNA stability and ability to bind mRNA targets, hence affecting target gene repression. Next Generation Sequencing (NGS) of the small RNA (sRNA) provides an efficient and reliable method to explore miRNA modifications. The need for dedicated software, especially for users with little knowledge of computers, to determine and analyze miRNA modifications in sRNA NGS data, motivated us to develop miRMOD. miRMOD is a user-friendly, Microsoft Windows and Graphical User Interface (GUI) based tool for identification and analysis of 5′ and 3′ miRNA modifications (non-templated nucleotide additions and trimming) in sRNA NGS data. In addition to identification of miRNA modifications, the tool also predicts and compares the targets of query and modified miRNAs. In order to compare binding affinities for the same target, miRMOD utilizes minimum free energies of the miRNA:target and modified-miRNA:target interactions. Comparisons of the binding energies may guide experimental exploration of miRNA post-transcriptional modifications. The tool is available as a stand-alone package to overcome large data transfer problems commonly faced in web-based high-throughput (HT) sequencing data analysis tools. miRMOD package is freely available at http://bioinfo.icgeb.res.in/miRMOD.

A Uniform System for the Annotation of Vertebrate microRNA Genes and the Evolution of the Human microRNAome

Although microRNAs (miRNAs) are among the most intensively studied molecules of the past 20 years, determining what is and what is not a miRNA has not been straightforward. Here, we present a uniform system for the annotation and nomenclature of miRNA genes. We show that less than a third of the 1,881 human miRBase entries, and only approximately 16% of the 7,095 metazoan miRBase entries, are robustly supported as miRNA genes. Furthermore, we show that the human repertoire of miRNAs has been shaped by periods of intense miRNA innovation and that mature gene products show a very different tempo and mode of sequence evolution than star products. We establish a new open access database--MirGeneDB ( http://mirgenedb.org )--to catalog this set of miRNAs, which complements the efforts of miRBase but differs from it by annotating the mature versus star products and by imposing an evolutionary hierarchy upon this curated and consistently named repertoire.

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