A meta-analysis revealed insights into the sources, conservation and impact of microRNA 5'-isoforms in four model species

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.

Genome-wide analysis of Panonychus citri microRNAs with a focus on potential insecticidal activity of 4 microRNAs to eggs and nymphs

Panonychus citri McGregor (Acari: Tetranychidae), a destructive citrus pest, causes considerable annual economic losses due to its short lifespan and rapid resistance development. MicroRNA (miRNA)-induced RNA interference is a promising approach for pest control because of endogenous regulation of pest growth and development. To search for miRNAs with potential insecticidal activity in P. citri, genome-wide analysis of miRNAs at different developmental stages was conducted, resulting in the identification of 136 miRNAs, including 73 known and 63 novel miRNAs. A total of 17 isomiRNAs and 12 duplicated miRNAs were characterized. MiR-1 and miR-252-5p were identified as reference miRNAs for P. citri and Tetranychus urticae. Based on differential expression analysis, treatments with miR-let-7a and miR-315 mimics and the miR-let-7a antagomir significantly reduced the egg hatch rate and resulted in abnormal egg development. Overexpression or downregulation of miR-34-5p and miR-305-5p through feeding significantly decreased the adult eclosion rate and caused molting defects. The 4 miRNAs, miR-let-7a, miR-315, miR-34-5p, and miR-305-5p, had important regulatory functions and insecticidal properties in egg hatching and adult eclosion. In general, these data advance our understanding of miRNAs in mite biology, which can assist future studies on insect-specific miRNA-based green pest control technology.

Adipose Tissue Exosome circ_sxc Mediates the Modulatory of Adiposomes on Brain Aging by Inhibiting Brain dme-miR-87-3p

Aging of the brain usually leads to the decline of neurological processes and is a major risk factor for various neurodegenerative diseases, including sleep disturbances and cognitive decline. Adipose tissue exosomes, as adipocyte-derived vesicles, may mediate the regulatory processes of adipose tissue on other organs, including the brain; however, the regulatory mechanisms remain unclear. We analyzed the sleep-wake behavior of young (10 days) and old (40 days) Drosophila and found that older Drosophila showed increased sleep fragmentation, which is similar to mammalian aging characteristics. To investigate the cross-tissue regulatory mechanisms of adiposity on brain aging, we extracted 10-day and 40-day Drosophila adipose tissue exosomes and identified circRNAs with age-dependent expression differences by RNA-seq and differential analysis. Furthermore, by combining data from 3 datasets of the GEO database (GSE130158, GSE24992, and GSE184559), circ_sxc that was significantly downregulated with age was finally screened out. Moreover, dme-miR-87-3p, a conserved target of circ_sxc, accumulates in the brain with age and exhibits inhibitory effects in predicted binding relationships with neuroreceptor ligand genes. In summary, the current study showed that the Drosophila brain could obtain circ_sxc by uptake of adipose tissue exosomes which crossed the blood-brain barrier. And circ_sxc suppressed brain miR-87-3p expression through sponge adsorption, which in turn regulated the expression of neurological receptor ligand proteins (5-HT1B, GABA-B-R1, Rdl, Rh7, qvr, NaCP60E) and ensured brain neuronal synaptic signaling normal function of synaptic signaling. However, with aging, this regulatory mechanism is dysregulated by the downregulation of the adipose exosome circ_sxc, which contributes to the brain exhibiting sleep disturbances and other "aging" features.

Pathogenesis of psoriasis in the "omic" era. Part II. Genetic, genomic and epigenetic changes in psoriasis

Psoriasis is a multifactorial disease in which genetic, environmental and epigenetic factors regulating gene expression play a key role. In the “genomic era”, genome-wide association studies together with target genotyping platforms performed in different ethnic populations have found more than 50 genetic susceptible markers associated with the risk of psoriasis which have been identified so far. Up till now, the strongest association with the risk of the disease has been proved for HLA-C*06 gene. The majority of other psoriasis risk SNPs are situated near the genes encoding molecules involved in adaptive and innate immunity, and skin barrier function. Many contemporary studies indicate that the epigenetic changes: histone modification, promoter methylations, long non-coding and micro-RNA hyperexpression are considered as factors contributing to psoriasis pathogenesis as they regulate abnormal keratinocyte differentiation and proliferation, aberrant keratinocytes – inflammatory cells communication, neoangiogenesis and chronic inflammation. The circulating miRNAs detected in the blood may become specific markers in the diagnosis, prognosis and response to the treatment of the disease. The inhibition of expression in selected miRNAs may be a new promising therapy option for patients with psoriasis.

MicroRNAs in Daphnia magna identified and characterized by deep sequencing, genome mapping and manual curation

MicroRNAs (miRNAs) are small non-coding RNAs that function in RNA silencing and post-transcriptional regulation of gene expression in most organisms. The water flea, Daphnia magna is a key model to study phenotypic, physiological and genomic responses to environmental cues and miRNAs can potentially mediate these responses. By using deep sequencing, genome mapping and manual curations, we have characterised the miRNAome of D. magna. We identified 66 conserved miRNAs and 13 novel miRNAs; all of these were found in the three studied life stages of D. magna (juveniles, subadults, adults), but with variation in expression levels between stages. Forty-one of the miRNAs were clustered into 13 genome clusters also present in the D. pulex genome. Most miRNAs contained sequence variants (isomiRs). The highest expressed isomiRs were 3′ template variants with one nucleotide deletion or 3′ non-template variants with addition of A or U at the 3′ end. We also identified offset RNAs (moRs) and loop RNAs (loRs). Our work extends the base for further work on all species (miRNA, isomiRs, moRNAs, loRNAs) of the miRNAome of Daphnia as biomarkers in response to chemical substances and environment cues, and underline age dependency.

miRBaseMiner, a tool for investigating miRBase content

microRNAs are small non-coding RNA molecules playing a central role in gene regulation. miRBase is the standard reference source for analysis and interpretation of experimental studies. However, the richness and complexity of the annotation is often underappreciated by users. Moreover, even for experienced users, the size of the resource can make it difficult to explore annotation to determine features such as species coverage, the impact of specific characteristics and changes between successive releases. A further consideration is that each new miRBase release contains entries that have had limited review and which may subsequently be removed in a future release to ensure the quality of annotation. To aid the miRBase user, we developed a software tool, miRBaseMiner, for investigating miRBase annotation and generating custom annotation sets. We apply the tool to characterize each release from v9.2 to v22 to examine how annotation has changed across releases and highlight some of the annotation features that users should keep in mind when using for miRBase for data analysis. These include: (1) entries with identical or very similar sequences; (2) entries with multiple annotated genome locations; (3) hairpin precursor entries with extremely low-estimated minimum free energy; (4) entries possessing reverse complementary; (5) entries with 3ʹ poly(A) ends. As each of these factors can impact the identification of dysregulated features and subsequent clinical or biological conclusions, miRBaseMiner is a valuable resource for any user using miRBase as a reference source.

MicroRNA-like RNAs from the same miRNA precursors play a role in cassava chilling responses

MicroRNAs (miRNAs) are known to play important roles in various cellular processes and stress responses. MiRNAs can be identified by analyzing reads from high-throughput deep sequencing. The reads realigned to miRNA precursors besides canonical miRNAs were initially considered as sequencing noise and ignored from further analysis. Here we reported a small-RNA species of phased and half-phased miRNA-like RNAs different from canonical miRNAs from cassava miRNA precursors detected under four distinct chilling conditions. They can form abundant multiple small RNAs arranged along precursors in a tandem and phased or half-phased fashion. Some of these miRNA-like RNAs were experimentally confirmed by re-amplification and re-sequencing, and have a similar qRT-PCR detection ratio as their cognate canonical miRNAs. The target genes of those phased and half-phased miRNA-like RNAs function in process of cell growth metabolism and play roles in protein kinase. Half-phased miR171d.3 was confirmed to have cleavage activities on its target gene P-glycoprotein 11, a broad substrate efflux pump across cellular membranes, which is thought to provide protection for tropical cassava during sharp temperature decease. Our results showed that the RNAs from miRNA precursors are miRNA-like small RNAs that are viable negative gene regulators and may have potential functions in cassava chilling responses.

MiRNA profiling of gastrointestinal stromal tumors by next-generation sequencing

Deregulation of miRNAs has been observed virtually in all major types of cancer, whereas the miRNA signature in GIST is not well characterized yet. In this study the first high-throughput miRNA profiling of 15 paired GIST and adjacent normal tissue samples was performed using small RNA-seq approach and differentially expressed miRNAs as well as isomiRNAs were defined. Highly significantly deregulated miRNAs were selected for validation by Taq-Man low-density array in replication group of 40 paired samples. Validated miRNAs were further subjected to enrichment analysis, which revealed significantly enriched KEGG pathways in the main GIST associated pathways. Further, we used an integrated analysis of miRNA-mRNA correlations for KIT and PDGFRA target genes and found a significant correlation between all of the enriched miRNAs and their target gene KIT. Results of the phenotype analysis showed miR-509-3p to be up-regulated in epithelioid and mixed cell types compared to spindle type, whereas miR-215-5p showed negative correlation with risk grade of GIST. These data reveal a detailed miRNA profile of GIST and highlight new candidates that may be important in the development of malignant disease.

MicroRNAs as regulatory elements in psoriasis

Psoriasis is a chronic, autoimmune, and complex genetic disorder that affects 23% of the European population. The symptoms of Psoriatic skin are inflammation, raised and scaly lesions. microRNA, which is short, nonprotein-coding, regulatory RNAs, plays critical roles in psoriasis. microRNA participates in nearly all biological processes, such as cell differentiation, development and metabolism. Recent researches reveal that multitudinous novel microRNAs have been identified in skin. Some of these substantial novel microRNAs play as a class of posttranscriptional gene regulator in skin disease, such as psoriasis. In order to insight into microRNAs biological functions and verify microRNAs biomarker, we review diverse references about characterization, profiling and subtype of microRNAs. Here we will share our opinions about how and which microRNAs are as regulatory in psoriasis.

Identification and characterization of novel and conserved microRNAs in several tissues of the Chinese rare minnow (Gobiocypris rarus) based on illumina deep sequencing technology

BACKGROUND

MicroRNAs (miRNAs), which comprise a large family of endogenous small non-coding RNA molecules, play important roles in the regulation of gene expression in various biological processes. The Chinese rare minnow (Gobiocypris rarus) is a Chinese native fish species and is used extensively as an experimental fish in China; however, relevant biological data, especially miRNA transcriptome data, have not been well documented. To discover conserved and potential novel miRNAs in Chinese rare minnows, a pool of equal amounts of RNA obtained from 6 different adult rare minnow tissues (brain, eye, gill, liver, muscle and heart) was sequenced using illumina deep sequencing technology.

RESULTS

In the present study, 26,930,553 raw reads, representing 2,118,439 unique high-quality reads, were obtained from the pooled small RNA library. Using bioinformatics analysis, 352 conserved and 112 novel Chinese rare minnow miRNAs were first discovered and characterized in this study. Moreover, we found extensive sequence variations (isomiRs) in rare minnow miRNAs, including internal miRNA isomiRs and terminal isomiRs at both the 5' and 3' ends and nucleotide variants. Six conserved and 4 novel miRNAs were selected and validated in 6 different adult rare minnow tissues using quantitative real-time PCR (qPCR). The results showed that miR-30a, miR-30b, and Novel-37 are ubiquitously expressed in a variety of tissues. miR-16a, miR-9, miR-125b, miR-34a, and Novel-69 were predominantly expressed in the brain. Novel-115 and Novel-7 were highly expressed in gills, but were relatively weakly expressed in other tissues. These results provided the expression patterns of miRNA genes in Chinese rare minnow. Finally, based on bioinformatics predictions, we mainly found that Novel-94 and Novel-1b-5p were simultaneously targeted to the 3'UTR of Dmrt1, which controls sex determination and/or sexual differentiation in a variety of metazoans at different sites. Novel-29b targeted the 3'UTR of Foxl2, which is involved in the maintenance of ovarian function and the transcriptional regulation of gonadal differentiation-related genes. Novel-62 and Novel-53 targeted the 3'UTR of ERbeta1 and ERbeta2 (which regulate the transcription of target genes), respectively.

CONCLUSIONS

Rare minnow is a widely used model for assessing the risk of environmental pollution in China. Identifying and characterizing rare minnow miRNA genes is necessary to discover the biological function of miRNAs and to screen for new molecule biomarkers to assess the risk of environmental pollution in the future.

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.

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

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

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

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

Deep sequencing, profiling and detailed annotation of microRNAs in Takifugu rubripes

Background

microRNAs (miRNAs) in fish have not been as extensively studied as those in mammals. The fish species Takifugu rubripes is an intensively studied model organism whose genome has been sequenced. The T. rubripes genome is approximately eight times smaller than the human genome, but has a similar repertoire of protein-coding genes. Therefore, it is useful for identifying non-coding genes, including miRNA genes. To identify miRNA expression patterns in different organs of T. rubripes and give fundamental information to aid understanding of miRNA populations in this species, we extracted small RNAs from tissues and performed deep sequencing analysis to profile T. rubripes miRNAs. These data will be of assistance in functional studies of miRNAs in T. rubripes.

Results

After analyzing a total of 139 million reads, we found miRNA species in nine tissues (fast and slow muscles, heart, eye, brain, intestine, liver, ovaries, and testes). We identified 1420 known miRNAs, many of which were strongly expressed in certain tissues with expression patterns similar to those described for other animals in previous reports. Most miRNAs were expressed in tissues other than the ovaries or testes. However, some miRNA families were highly abundant in the gonads, but expressed only at low levels in somatic tissue, suggesting specific function in germ cells. The most abundant isomiRs (miRNA variants) of many miRNAs had identical sequences in the 5′ region. However, isomiRs of some miRNAs, including fru-miR-462-5p, varied in the 5′ region in some tissues, suggesting that they may target different mRNA transcripts. Longer small RNAs (26–31 nt), which were abundant in the gonads, may be putative piRNAs because of their length and their origin from repetitive elements. Additionally, our data include possible novel classes of small RNAs.

Conclusions

We elucidated miRNA expression patterns in various organs of T. rubripes. Most miRNA sequences are conserved in vertebrates, indicating that the basic functions of vertebrate miRNAs share a common evolution. Some miRNA species exhibit different distributions of isomiRs between tissues, suggesting that they have a broad range of functions.

Electronic supplementary material

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

Variability of miRNA expression during the differentiation of human embryonic stem cells into retinal pigment epithelial cells

Embryonic stem cells (ESCs) and induced pluripotent stem cells can be induced to differentiate into retinal pigment epithelium (RPE). MiRNAs have been characterized and found playing important roles in the differentiation process of ESCs, but their length and sequence heterogeneity (isomiRs), and their non-canonical forms of miRNAs are underestimated or ignored. In this report, we found some non-canonical miRNAs (dominant isomiRs) in all differentiation stages, and 27 statistically significant editing sites were identified in 24 different miRNAs. Moreover, we found marked major-to-minor arm-switching events in 14 pre-miRNAs during the hESC to RPE cell differentiation phases. Our study for the first time reports exploring the variability of miRNA expression during the differentiation of hESCs into RPE cells and the results show that miRNA variability is a ubiquitous phenomenon in the ESC differentiation.

MIR137 variants identified in psychiatric patients affect synaptogenesis and neuronal transmission gene sets

Sequence analysis of 13 microRNA (miRNA) genes expressed in the human brain and located in genomic regions associated with schizophrenia and/or bipolar disorder, in a northern Swedish patient/control population, resulted in the discovery of two functional variants in the MIR137 gene. On the basis of their location and the allele frequency differences between patients and controls, we explored the hypothesis that the discovered variants impact the expression of the mature miRNA and consequently influence global mRNA expression affecting normal brain functioning. Using neuronal-like SH-SY5Y cells, we demonstrated significantly reduced mature miR-137 levels in the cells expressing the variant miRNA gene. Subsequent transcriptome analysis showed that the reduction in miR-137 expression led to the deregulation of gene sets involved in synaptogenesis and neuronal transmission, all implicated in psychiatric disorders. Our functional findings add to the growing data, which implicate that miR-137 has an important role in the etiology of psychiatric disorders and emphasizes its involvement in nervous system development and proper synaptic function.

YM500v2: a small RNA sequencing (smRNA-seq) database for human cancer miRNome research

We previously presented YM500, which is an integrated database for miRNA quantification, isomiR identification, arm switching discovery and novel miRNA prediction from 468 human smRNA-seq datasets. Here in this updated YM500v2 database (http://ngs.ym.edu.tw/ym500/), we focus on the cancer miRNome to make the database more disease-orientated. New miRNA-related algorithms developed after YM500 were included in YM500v2, and, more significantly, more than 8000 cancer-related smRNA-seq datasets (including those of primary tumors, paired normal tissues, PBMC, recurrent tumors, and metastatic tumors) were incorporated into YM500v2. Novel miRNAs (miRNAs not included in the miRBase R21) were not only predicted by three independent algorithms but also cleaned by a new in silico filtration strategy and validated by wetlab data such as Cross-Linked ImmunoPrecipitation sequencing (CLIP-seq) to reduce the false-positive rate. A new function 'Meta-analysis' is additionally provided for allowing users to identify real-time differentially expressed miRNAs and arm-switching events according to customer-defined sample groups and dozens of clinical criteria tidying up by proficient clinicians. Cancer miRNAs identified hold the potential for both basic research and biotech applications.

MicroRNAs in normal and psoriatic skin

Psoriasis is a chronic and common human skin disorder currently with no cure. Psoriatic skin displays inflammatory, raised, and scaly lesions with widely aberrant gene expression. Recent studies have revealed critical roles that microRNAs play as a class of posttranscriptional gene regulator in skin development and skin diseases. A substantial number of novel microRNAs have been identified in skin, and much has been learned about the dysregulated expression and functional roles of microRNAs in psoriasis, as well as the robustness and plasticity of microRNA-mediated gene expression regulation. Here we review recent progresses in discovery, profiling, and characterization of microRNAs in human psoriatic skin, discuss insights to their biological functions, and share our view on remaining challenges to be addressed.