Mouse let-7 miRNA populations exhibit RNA editing that is constrained in the 5'-seed/ cleavage/anchor regions and stabilize predicted mmu-let-7a:mRNA duplexes

Genome-wide investigation and functional analysis of RNA editing sites in wheat

Wheat is an important cereal and half of the world population consumed it. Wheat faces environmental stresses and different techniques (CRISPR, gene silencing, GWAS, etc.) were used to enhance its production but RNA editing (RESs) is not fully explored in wheat. RNA editing has a special role in controlling environmental stresses. The genome-wide identification and functional characterization of RESs in different types of wheat genotypes was done. We employed six wheat genotypes by RNA-seq analyses to achieve RESs. The findings revealed that RNA editing events occurred on all chromosomes equally. RNA editing sites were distributed randomly and 10–12 types of RESs were detected in wheat genotypes. Higher number of RESs were detected in drought-tolerant genotypes. A-to-I RNA editing (2952, 2977, 1916, 2576, 3422, and 3459) sites were also identified in six wheat genotypes. Most of the genes were found to be engaged in molecular processes after a Gene Ontology analysis. PPR (pentatricopeptide repeat), OZ1 (organelle zinc-finger), and MORF/RIP gene expression levels in wheat were also examined. Normal growth conditions diverge gene expression of these three different gene families, implying that normal growth conditions for various genotypes can modify RNA editing events and have an impact on gene expression levels. While the expression of PPR genes was not change. We used Variant Effect Predictor (VEP) to annotate RNA editing sites, and Local White had the highest RESs in the CDS region of the protein. These findings will be useful for prediction of RESs in other crops and will be helpful in drought tolerance development in wheat.

Detecting and Characterizing A-To-I microRNA Editing in Cancer

Adenosine to inosine (A-to-I) editing consists of an RNA modification where single adenosines along the RNA sequence are converted into inosines. Such a biochemical transformation is catalyzed by enzymes belonging to the family of adenosine deaminases acting on RNA (ADARs) and occurs either co- or post-transcriptionally. The employment of powerful, high-throughput detection methods has recently revealed that A-to-I editing widely occurs in non-coding RNAs, including microRNAs (miRNAs). MiRNAs are a class of small regulatory non-coding RNAs (ncRNAs) acting as translation inhibitors, known to exert relevant roles in controlling cell cycle, proliferation, and cancer development. Indeed, a growing number of recent researches have evidenced the importance of miRNA editing in cancer biology by exploiting various detection and validation methods. Herein, we briefly overview early and currently available A-to-I miRNA editing detection and validation methods and discuss the significance of A-to-I miRNA editing in human cancer.

Presynaptic dysfunction in CASK-related neurodevelopmental disorders

CASK-related disorders are genetically defined neurodevelopmental syndromes. There is limited information about the effects of CASK mutations in human neurons. Therefore, we sought to delineate CASK-mutation consequences and neuronal effects using induced pluripotent stem cell-derived neurons from two mutation carriers. One male case with autism spectrum disorder carried a novel splice-site mutation and a female case with intellectual disability carried an intragenic tandem duplication. We show reduction of CASK protein in maturing neurons from the mutation carriers, which leads to significant downregulation of genes involved in presynaptic development and of CASK protein interactors. Furthermore, CASK-deficient neurons showed decreased inhibitory presynapse size as indicated by VGAT staining, which may alter the excitatory-inhibitory (E/I) balance in developing neural circuitries. Using in vivo magnetic resonance spectroscopy quantification of GABA in the male mutation carrier, we further highlight the possibility to validate in vitro cellular data in the brain. Our data show that future pharmacological and clinical studies on targeting presynapses and E/I imbalance could lead to specific treatments for CASK-related disorders.

microRNA-mRNA analysis in pituitary and hypothalamus of sterile Japanese flounder

Double haploidy is an advantageous situation for genetic mapping and genome sequencing studies. In the present study, the hypothalamus and pituitary gland from sterile and fertile double-haploid (DH) Japanese flounders (aged 5 years) were used as experimental materials for studying the expression of genes in individuals with reproductive disorders, using high-throughput sequencing technology. The results revealed abnormal levels of some hormones in sterile DHs during the breeding season. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that the significantly different microRNAs and messenger RNAs were related to metabolism, signal transduction, and melanogenesis; those related to steroid hormone synthesis and secretion related pathways were not detected. Our results suggest that the key to sterility in DHs was the arrested ovary development. However, the reason for arrested ovary development was mainly related to the lower levels of expression of genes involved in steroid biosynthesis in gonads, and was not related to the pituitary. For maintaining homeostasis, the hypothalamus and pituitary would have large differences in several processes, including signal transduction, metabolism, and immune response. The present study provides primary data for further studies on sterility in fish, and even in other animals.

Integrated mRNA and miRNA expression profile analyses reveal the potential roles of sex-biased miRNA-mRNA pairs in gonad tissues of the Chinese concave-eared torrent frog (Odorrana tormota)

The Chinese concave-eared torrent frog (Odorrana tormota) is typically sexually dimorphic. Females are significantly less common than males in the wild. Until now, the molecular mechanisms of reproduction and sex differentiation of frogs remain unclear. Here, we integrated mRNA and microRNA (miRNA) expression profiles to reveal the molecular mechanisms of reproduction and sex differentiation in O. tormota. We identified 234 differentially expressed miRNAs (DEMs) and 18,551 differentially expressed transcripts. Of these, 12,053 mRNAs and 64 miRNAs were upregulated in testes, and 6,498 mRNAs and 170 miRNAs were upregulated in ovaries. Integrated analysis of the miRNA and mRNA expression profiles predicted 75,602 potential miRNA-mRNA interaction sites, with 42,065 negative miRNA-mRNA interactions. We found 36 differentially expressed genes (DEGs) related to reproduction and sex differentiation, of which 15 DEGs formed 92 negative miRNA-mRNA interactions with 34 known DEMs. Thus, miRNAs may play other important roles in O. tormota. Furthermore, Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed reproductive-related processes, such as the gonadotropinreleasing hormone signaling pathway and ovarian steroidogenesis. Based on functional annotation and the literature, the retinoic acid signaling pathway, the SOX9-AMH pathway, and the process of spermatogenesis may be involved in the molecular mechanisms of reproduction and sex differentiation in O. tormota, and may be regulated by miRNAs. The miRNA-mRNA pairs described may provide further understanding of the regulatory mechanisms associated with reproduction and sex differentiation, and the molecular mechanism of reproduction in O. tormota.

H19 contributes to poor clinical features in NSCLC patients and leads to enhanced invasion in A549 cells through regulating miRNA-203-mediated epithelial-mesenchymal transition

Recent studies have demonstrated that the overexpression of H19 may contribute towards development of tumorigenesis in various types of cancer. To investigate the role of H19 in the development of non-small cell lung cancer (NSCLC), 76 NSCLC tissues samples and their adjacent normal tissue samples were collected. Expression level of H19, and its association with clinicopathological features and overall survival was analyzed. It was found that compared with normal adjacent tissues, H19 expression was elevated in NSCLC tissues along with a decreased miR-203 expression level. It was also found that patients who were in advanced clinical stages had a higher H19 and a lower miR-203 expression compared to normal tissues. The overall survival time of patients with higher H19 expression was shorter compared with the lower H19 expression group. Upregulation of A549 enhanced cell proliferation and promoted invasion. Overexpression of H19 stimulated the epithelial-mesenchymal transition (EMT) process in lung cancer cells and demonstrated typical morphological characteristics of EMT. The level of mesenchymal marker protein, such as Vimentin and SNAI1 increased; while CDH1 protein level decreased. Also, H19 negatively regulated miR-203. Inhibition of H19 attenuated miR-203 induced EMT process. Upregulation of H19 contributes to poor clinical features in patients with NSCLC, induces occurrence of EMT, promotes proliferation and stimulates cell invasion in NSCLC cell line through regulating miRNA-203 mediated EMT.

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

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

Profiling and Bioinformatic Analyses Indicate Differential circRNA and miRNA/isomiR Expression and Interactions

As a novel class of noncoding RNAs, circular RNAs (circRNAs) have been reported to play a role in various biological processes. Some circRNAs may serve as microRNA (miRNA) sponges, regulating transcription or splicing. Herein, we investigated the expression profiles and interactions of miRNAs/isomiRs and circRNAs in male patients with esophageal cancer. We found that some miRNA genes generated two deregulated miRNA products (miR-#-5p and miR-#-3p), and these products were consistently abnormally expressed. Some circRNAs were predicted to be miRNA sponges for specific miRNAs. Some of these typically showed opposing expression patterns in cancer tissues: one upregulated and the other downregulated. Although fewer miRNAs were predicted to interact with circRNAs, the number of predicted interactions would be substantially increased if detailed isomiRs were involved. High sequence similarity across multiple isomiRs suggested that they might interact with circRNAs, similar to the interaction of homologous miRNAs with circRNAs. At the isomiR level, due to the characteristics of the sequences and expression patterns involved, the cross-talk between different ncRNAs is complicated despite simplification of the isomiRs involved through clustering. We expect that our results may provide methods for further study of the cross-talk among ncRNAs and elucidate their biological roles in human diseases.

High-throughput sequencing of pituitary and hypothalamic microRNA transcriptome associated with high rate of egg production

Background

MicroRNAs exist widely in viruses, plants and animals. As endogenous small non-coding RNAs, miRNAs regulate a variety of biological processes. Tissue miRNA expression studies have discovered numerous functions for miRNAs in various tissues of chicken, but the regulation of miRNAs in chicken pituitary and hypothalamic development related to high and low egg-laying performance has remained unclear.

Results

In this study, using high-throughput sequencing technology, we sequenced two tissues (pituitary and hypothalamus) in 3 high- and 3 low-rate egg production Luhua chickens at the age of 300 days. By comparing low- and high-rate egg production chickens, 46 known miRNAs and 27 novel miRNAs were identified as differentially expressed (P < 0.05). Six differentially expressed known miRNAs, which are expressed in both tissues, were used in RT-qPCR validation and SNP detection. Among them, seven SNPs in two miRNA precursors (gga-miR-1684a and gga-miR-1434) were found that might enhance or reduce the production of the mature miRNAs. In addition, 124 and 30 reciprocally expressed miRNA-target pairs were identified by RNA-seq in pituitary and hypothalamic tissues, respectively and randomly selected candidate miRNA and miRNA-target pairs were validated by RT-qPCR in Jiuyuan black fowl. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotation illustrated that a large number of egg laying-related pathways were enriched in the high-rate egg production chickens, including ovarian steroidogenesis and steroid hormone biosynthesis.

Conclusions

These differentially expressed miRNAs and their predicted target genes, especially identified reciprocally expressed miRNA-target pairs, advance the study of miRNA function and egg production associated miRNA identification. The analysis of the miRNA-related SNPs and their effects provided insights into the effects of SNPs on miRNA biogenesis and function. The data generated in this study will further our understanding of miRNA regulation mechanisms in the chicken egg-laying process.

Electronic supplementary material

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

RNA binding to APOBEC deaminases; Not simply a substrate for C to U editing

Apolipoprotein B mRNA Editing Catalytic Polypeptide-like 1 or APOBEC1 was discovered in 1993 as the zinc-dependent cytidine deaminase responsible for the production of an in frame stop codon in apoB mRNA through modification of cytidine at nucleotide position 6666 to uridine. At the time of this discovery there was much speculation concerning the mechanism of base modification RNA editing which has been rekindled by the discovery of multiple C to U RNA editing events in the 3′ UTRs of mRNAs and the finding that other members of the APOBEC family while able to bind RNA, have the biological function of being DNA mutating enzymes. Current research is addressing the mechanism for these nucleotide modification events that appear not to adhere to the mooring sequence-dependent model for APOBEC1 involving the assembly of a multi protein containing editosome. This review will summarize our current understanding of the structure and function of APOBEC proteins and examine how RNA binding to them may be a regulatory mechanism.

Expressed microRNA associated with high rate of egg production in chicken ovarian follicles

MicroRNA (miRNA) is a highly conserved class of small noncoding RNA about 19-24 nucleotides in length that function in a specific manner to post-transcriptionally regulate gene expression in organisms. Tissue miRNA expression studies have discovered a myriad of functions for miRNAs in various aspects, but a role for miRNAs in chicken ovarian tissue at 300 days of age has not hitherto been reported. In this study, we performed the first miRNA analysis of ovarian tissues in chickens with low and high rates of egg production using high-throughput sequencing. By comparing low rate of egg production chickens with high rate of egg production chickens, 17 significantly differentially expressed miRNAs were found (P < 0.05), including 11 known and six novel miRNAs. We found that all 11 known miRNAs were involved mainly in pathways of reproduction regulation, such as steroid hormone biosynthesis and dopaminergic synapse. Additionally, expression profiling of six randomly selected differentially regulated miRNAs were validated by quantitative real-time polymerase chain reaction (RT-qPCR). Some miRNAs, such as gga-miR-34b, gga-miR-34c and gga-miR-216b, were reported to regulate processes such as proliferation, cell cycle, apoptosis and metastasis and were expressed differentially in ovaries of chickens with high rates of egg production, suggesting that these miRNAs have an important role in ovary development and reproductive management of chicken. Furthermore, we uncovered that a significantly up-regulated miRNA-gga-miR-200a-3p-is ubiquitous in reproduction-regulation-related pathways. This miRNA may play a special central role in the reproductive management of chicken, and needs to be further studied for confirmation.

Accurate detection for a wide range of mutation and editing sites of microRNAs from small RNA high-throughput sequencing profiles

Various types of mutation and editing (M/E) events in microRNAs (miRNAs) can change the stabilities of pre-miRNAs and/or complementarities between miRNAs and their targets. Small RNA (sRNA) high-throughput sequencing (HTS) profiles can contain many mutated and edited miRNAs. Systematic detection of miRNA mutation and editing sites from the huge volume of sRNA HTS profiles is computationally difficult, as high sensitivity and low false positive rate (FPR) are both required. We propose a novel method (named MiRME) for an accurate and fast detection of miRNA M/E sites using a progressive sequence alignment approach which refines sensitivity and improves FPR step-by-step. From 70 sRNA HTS profiles with over 1.3 billion reads, MiRME has detected thousands of statistically significant M/E sites, including 3'-editing sites, 57 A-to-I editing sites (of which 32 are novel), as well as some putative non-canonical editing sites. We demonstrated that a few non-canonical editing sites were not resulted from mutations in genome by integrating the analysis of genome HTS profiles of two human cell lines, suggesting the existence of new editing types to further diversify the functions of miRNAs. Compared with six existing studies or methods, MiRME has shown much superior performance for the identification and visualization of the M/E sites of miRNAs from the ever-increasing sRNA HTS profiles.

Comprehensive analysis of small RNA-seq data reveals that combination of miRNA with its isomiRs increase the accuracy of target prediction in Arabidopsis thaliana

Along with the canonical miRNA, distinct miRNA-like sequences called sibling miRNAs (sib-miRs) are generated from the same pre-miRNA. Among them, isomeric sequences featuring slight variations at the terminals, relative to the canonical miRNA, constitute a pool of isomeric sibling miRNAs (isomiRs). Despite the high prevalence of isomiRs in eukaryotes, their features and relevance remain elusive. In this study, we performed a comprehensive analysis of mature precursor miRNA (pre-miRNA) sequences from Arabidopsis to understand their features and regulatory targets. The influence of isomiR terminal heterogeneity in target binding was examined comprehensively. Our comprehensive analyses suggested a novel computational strategy that utilizes miRNA and its isomiRs to enhance the accuracy of their regulatory target prediction in Arabidopsis. A few targets are shared by several members of isomiRs; however, this phenomenon was not typical. Gene Ontology (GO) enrichment analysis showed that commonly targeted mRNAs were enriched for certain GO terms. Moreover, comparison of these commonly targeted genes with validated targets from published data demonstrated that the validated targets are bound by most isomiRs and not only the canonical miRNA. Furthermore, the biological role of isomiRs in target cleavage was supported by degradome data. Incorporating this finding, we predicted potential target genes of several miRNAs and confirmed them by experimental assays. This study proposes a novel strategy to improve the accuracy of predicting miRNA targets through combined use of miRNA with its isomiRs.

Identification of miRNAs during mouse postnatal ovarian development and superovulation

BACKGROUND

MicroRNAs are small noncoding RNAs that play critical roles in regulation of gene expression in wide array of tissues including the ovary through sequence complementarity at post-transcriptional level. Tight regulation of multitude of genes involved in ovarian development and folliculogenesis could be regulated at transcription level by these miRNAs. Therefore, tissue specific miRNAs identification is considered a key step towards understanding the role of miRNAs in biological processes.

METHODS

To investigate the role of microRNAs during ovarian development and folliculogenesis we sequenced eight different libraries using Illumina deep sequencing technology. Different developmental stages were selected to explore miRNAs expression pattern at different stages of gonadal maturation with/without treatment of PMSG/hCG for superovulation.

RESULTS

From massive sequencing reads, clean reads of 16-26 bp were selected for further analysis of differential expression analysis and novel microRNA annotation. Expression analysis of all miRNAs at different developmental stages showed that some miRNAs were present ubiquitously while others were differentially expressed at different stages. Among differentially expressed miRNAs we reported 61 miRNAs with a fold change of more than 2 at different developmental stages among all libraries. Among the up-regulated miRNAs, mmu-mir-1298 had the highest fold change with 4.025 while mmu-mir-150 was down-regulated more than 3 fold. Furthermore, we found 2659 target genes for 20 differentially expressed microRNAs using seven different target predictions programs (DIANA-mT, miRanda, miRDB, miRWalk, RNAhybrid, PICTAR5, TargetScan). Analysis of the predicted targets showed certain ovary specific genes targeted by single or multiple microRNAs. Furthermore, pathway annotation and Gene ontology showed involvement of these microRNAs in basic cellular process.

CONCLUSIONS

These results suggest the presence of different miRNAs at different stages of ovarian development and superovulation. Potential role of these microRNAs was elucidated using bioinformatics tools in regulation of different pathways, biological functions and cellular components underlying ovarian development and superovulation. These results provide a framework for extended analysis of miRNAs and their roles during ovarian development and superovulation. Furthermore, this study provides a base for characterization of individual miRNAs to discover their role in ovarian development and female fertility.

An exploration of evolution, maturation, expression and function relationships in mir-23 ∼ 27 ∼ 24 cluster

The study aims to explore the potential relationships of evolution, maturation, expression and function between homologous/clustered miRNAs. mir-23∼27∼24 gene cluster, including the two gene clusters (mir-23a and mir-23b) and the three miRNA gene families (mir-23, mir-27 and mir-24), was typically selected as an example. These related miRNAs show similar evolutionary patterns and various expression patterns. Most of them show consistent isomiR expression pattern, and the “switching” phenomenon can be found between different abundant isomiR species. These findings suggest that these sequence or location related miRNAs show the similar miRNA processing and maturation processes, and the robust selection of the most dominant isomiR exists in specific tissues. Functional analysis show that these miRNAs show similar distributions of enriched gene categories, suggesting the close functional prelateships via direct or indirect coordinate regulation in biological processes. The study reveals the close evolutionary, expression and functional relationships between related homologous/clustered miRNAs, which will further enrich miRNA studies and understand direct or indirect interactions between miRNAs.

A genome-wide screen for non-template nucleotides and isomiR repertoires in miRNAs indicates dynamic and versatile microRNAome

miRNA variants (termed isomiRs) have been reported as potential functional molecules that may affect miRNA stability or target selection. The aims of the present study were to comprehensively survey and characterize non-template nucleotides (NTNs) and isomiR repertoires in miRNAs. Over 50 % of the NTNs were located in the 3' ends (also termed 3' additions), followed by the 5' ends and adjacent positions to 5' and 3' ends. The similar distributions of NTNs and isomiR repertoires might be detected between homologous or clustered miRNAs. miRNA might be stably expressed based on the typical analysis, but its isomiRs might be strongly up- or down-regulated. IsomiRs with novel seed sequences were mainly derived from "seed shifting" events in 5' isomiRs, NTNs in 5' ends or in seed sequences. IsomiRs from a miRNA locus or homologous miRNA loci maybe have the same seed sequences, but they would have various enrichment levels and 3' ends. Interestingly, isomiR species with novel seed sequences via NTNs in the seed region were always stably expressed. These novel seed sequences could lead to novel functional roles through driving the potential novel target mRNAs. Integrated predicted target mRNAs and further microarray validation showed that these isomiRs have versatile biological roles. Collectively, multiple isomiR products and miRNA maturation processes provide opportunities to perform versatile roles in the regulatory network, which further enriches and complicates the regulation of biological processes.

Widespread anti-sense transcription in apple is correlated with siRNA production and indicates a large potential for transcriptional and/or post-transcriptional control

Characterizing the transcriptome of eukaryotic organisms is essential for studying gene regulation and its impact on phenotype. The realization that anti-sense (AS) and noncoding RNA transcription is pervasive in many genomes has emphasized our limited understanding of gene transcription and post-transcriptional regulation. Numerous mechanisms including convergent transcription, anti-correlated expression of sense and AS transcripts, and RNAi remain ill-defined. Here, we have combined microarray analysis and high-throughput sequencing of small RNAs (sRNAs) to unravel the complexity of transcriptional and potential post-transcriptional regulation in eight organs of apple (Malus × domestica). The percentage of AS transcript expression is higher than that identified in annual plants such as rice and Arabidopsis thaliana. Furthermore, we show that a majority of AS transcripts are transcribed beyond 3'UTR regions, and may cover a significant portion of the predicted sense transcripts. Finally we demonstrate at a genome-wide scale that anti-sense transcript expression is correlated with the presence of both short (21-23 nt) and long (> 30 nt) siRNAs, and that the sRNA coverage depth varies with the level of AS transcript expression. Our study provides a new insight on the functional role of anti-sense transcripts at the genome-wide level, and a new basis for the understanding of sRNA biogenesis in plants.

A challenge for miRNA: multiple isomiRs in miRNAomics

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

Next-Generation Sequencing RNA-Seq Library Construction

This unit presents protocols for construction of next-generation sequencing (NGS) directional RNA sequencing libraries for the Illumina HiSeq and MiSeq from a wide variety of input RNA sources. The protocols are based on the New England Biolabs (NEB) small RNA library preparation set for Illumina, although similar kits exist from different vendors. The protocol preserves the orientation of the original RNA in the final sequencing library, enabling strand-specific analysis of the resulting data. These libraries have been used for differential gene expression analysis and small RNA discovery and are currently being tested for de novo transcriptome assembly. The protocol is robust and applicable to a broad range of RNA input types and RNA quality, making it ideal for high-throughput laboratories.

Imperfect centered miRNA binding sites are common and can mediate repression of target mRNAs

Background

MicroRNAs (miRNAs) bind to mRNAs and target them for translational inhibition or transcriptional degradation. It is thought that most miRNA-mRNA interactions involve the seed region at the 5′ end of the miRNA. The importance of seed sites is supported by experimental evidence, although there is growing interest in interactions mediated by the central region of the miRNA, termed centered sites. To investigate the prevalence of these interactions, we apply a biotin pull-down method to determine the direct targets of ten human miRNAs, including four isomiRs that share centered sites, but not seeds, with their canonical partner miRNAs.

Results

We confirm that miRNAs and their isomiRs can interact with hundreds of mRNAs, and that imperfect centered sites are common mediators of miRNA-mRNA interactions. We experimentally demonstrate that these sites can repress mRNA activity, typically through translational repression, and are enriched in regions of the transcriptome bound by AGO. Finally, we show that the identification of imperfect centered sites is unlikely to be an artifact of our protocol caused by the biotinylation of the miRNA. However, the fact that there was a slight bias against seed sites in our protocol may have inflated the apparent prevalence of centered site-mediated interactions.

Conclusions

Our results suggest that centered site-mediated interactions are much more frequent than previously thought. This may explain the evolutionary conservation of the central region of miRNAs, and has significant implications for decoding miRNA-regulated genetic networks, and for predicting the functional effect of variants that do not alter protein sequence.

MicroRNA expression profiles in human breast cancer cells after multifraction and single-dose radiation treatment

MicroRNAs (miRNAs) are small non-coding RNAs that contribute to modulating signaling pathways after radiation exposure and have emerged as a potential therapeutic target or biomarker in the radiation response of cancer. Exposing breast cancer cells to single-dose (SD) or multifractionated (MF) radiation may affect the cells differently. However, the roles of miRNAs in breast cancer cells after the response to SD or MF is not thoroughly understood. Therefore, the purpose of the present study was to comprehensively investigate the response of miRNAs in MDA-MB-361 by using various radiation exposing protocols. Our results revealed that only a small fraction of miRNAs exhibiting differential expressions (>1.5‑fold) was identified after MDA-MB-361 cells were exposed to SD (10 Gy) or MF radiation (2 Gy x 5 MF). In addition, we observed that several miRNAs in the MDA-MB-361 cells frequently exhibited differential responses to various types of radiation treatment. Among these miRNAs, the expression levels of an oncogenic miR-17-92 cluster increased following SD radiation treatment. Conversely, miR-19a-3p, miR-20a-5p, and miR-19b-3p expressions were inhibited by >1.5-fold in the following MF treatment. Further analysis of the miR-17-92 cluster expression levels revealed that miR-17, miR-18a, miR-19a/b and miR-20a were significantly overexpressed and miR-92a was downregulated in breast cancer. Functional annotation demonstrated that target genes of the miR-17-92 cluster were predominantly involved in the regulation of radiation-associated signal pathways such as mitogen-activated protein kinase (MAPK), ErbB, p53, Wnt, transforming growth factor-β (TGF-β), mTOR signaling pathways and cell cycles with an FDR <0.05. Overall, the results of the present study revealed distinct differences in the response of miRNAs to SD and MF radiation exposure, and these radiation-associated miRNAs may contribute to radiosensitivity and can be used as biomarkers for radiotherapy.

Comprehensive microRNA profiling of prostate cancer cells after ionizing radiation treatment

MicroRNAs (miRNAs) are small, non-coding RNAs that negatively regulate gene expression and have emerged as potential biomarkers in radiation response to human cancer. Only a few miRNAs have been identified in radiation response to prostate cancer and the involvement of the radiation-associated miRNA machinery in the response of prostate cancer cells to radiation is not thoroughly understood. Therefore, the purpose of the present study was to comprehensively investigate the expression levels, arm selection preference and isomiRs of radiation-response miRNAs in radiation-treated PC3 cells using a next-generation sequencing (NGS) approach. Our data revealed that the arm selection preference and 3′ modification of miRNAs may be altered in prostate cancer after radiation exposure. In addition, the proportion of AA dinucleotide modifications at the end of the read gradually increased in a time-dependent manner after PC3 radiation treatment. We also identified 6 miRNAs whose expression increased and 16 miRNAs whose expression decreased after exposure to 10 Gy of radiation. A pathway enrichment analysis revealed that the target genes of these radiation-induced miRNAs significantly co-modulated the radiation response pathway, including the mitogen-activated protein kinase (MAPK), Wnt, transforming growth factor-β (TGF-β) and ErbB signaling pathways. Furthermore, analysis of The Cancer Genome Atlas (TCGA) database revealed that the expression of these radiation-induced miRNAs was frequently dysregulated in prostate cancer. Our study identified radiation-induced miRNA candidates which may contribute to radiosensitivity and can be used as biomarkers for radiotherapy.

MicroRNAs in platelet physiology and pathology

MicroRNAs (miRNAs), highly conserved, short (approx. 22 nucleotides) non-coding RNAs, exhibit a fine-tune control over gene expression by complementary sequence binding and translational repression of protein coding mRNA transcripts. Recently, the role of miRNAs has been increasingly investigated in various physiological or pathophysiological events. Circulating platelets are crucial for coagulation physiology to maintain haemostatic balance and are involved in various pathologies such as atherosclerosis and thrombosis. Anucleate platelets lack genomic DNA but inherit diverse array of functional coding or non-coding RNAs and translational machinery from their parent cells - megakaryocytes enabling activated platelets to synthesize proteins which suggests the possibility of post transcriptional gene regulation in platelets. Expression of functionally active miRNAs in platelets changes during platelet activation indicating a role in platelet biology. Here, we present a review on recently identified platelet miRNAs and their role in platelet physiology that is essential for maintaining haemostasis.

MicroRNAs: new candidates for the regulation of the human cumulus-oocyte complex

STUDY QUESTION

What is the expression pattern of microRNAs (miRNAs) in human cumulus-oocyte complexes (COCs)?

SUMMARY ANSWER

Several miRNAs are enriched in cumulus cells (CCs) or oocytes, and are predicted to target genes involved in biological functions of the COC.

WHAT IS KNOWN ALREADY

The transcriptional profiles of human MII oocytes and the surrounding CCs are known. However, very limited data are available about post-transcriptional regulators, such as miRNAs. This is the first study focussing on the identification and quantification of small RNAs, including miRNAs, in human oocytes and CCs using a deep-sequencing approach.

STUDY DESIGN, SIZE, DURATION

MII oocytes and CCs were collected from women who underwent IVF.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Using the Illumina/deep-sequencing technology, we analyzed the small RNAome of pooled MII oocytes (n = 24) and CC samples (n = 20). The mRNA targets of CC and MII oocyte miRNAs were identified using in silico prediction algorithms. Using oligonucleotide microarrays, genome-wide gene expression was studied in oocytes (10 pools of 19 ± 3 oocytes/each) and 10 individual CC samples. TaqMan miRNA assays were used to confirm the sequencing results in independent pools of MII oocytes (3 pools of 8 ± 3 oocytes/each) and CC samples (3 pools of 7 ± 3 CCs/each). The functional role of one miRNA, MIR23a, was assessed in primary cultures of human CCs.

MAIN RESULTS AND THE ROLE OF CHANCE

Deep sequencing of small RNAs yielded more than 1 million raw reads. By mapping reads with a single location to the human genome, known miRNAs that were abundant in MII oocytes (MIR184, MIR100 and MIR10A) or CCs (MIR29a, MIR30d, MIR21, MIR93, MIR320a, MIR125a and the LET7 family) were identified. Predicted target genes of the oocyte miRNAs were associated with the regulation of transcription and cell cycle, whereas genes targeted by CC miRNAs were involved in extracellular matrix and apoptosis. Comparison of the predicted miRNA target genes and mRNA microarray data resulted in a list of 224 target genes that were differentially expressed in MII oocytes and CCs, including PTGS2, CTGF and BMPR1B that are important for cumulus-oocyte communication. Functional analysis using primary CC cultures revealed that BCL2 and CYP19A1 mRNA levels were decreased upon MIR23a overexpression.

LIMITATIONS, REASONS FOR CAUTION

Only known miRNAs were investigated in the present study on COCs. Moreover, the source of the material is MII oocytes that failed to fertilize.

WIDER IMPLICATIONS OF THE FINDINGS

The present findings suggest that miRNA could play a role in the regulation of the oocyte and CC crosstalk.

STUDY FUNDING/COMPETING INTEREST(S)

This work was partially supported by a grant from Ferring Pharmaceuticals. The authors of the study have no conflict of interest to report.

TRIAL REGISTRATION NUMBER

Not applicable.

Identification of miRNAs associated with sexual maturity in chicken ovary by Illumina small RNA deep sequencing

Background

MicroRNAs have been suggested to play important roles in the regulation of gene expression in various biological processes. To investigate the function of miRNAs in chicken ovarian development and folliculogenesis, two small RNA libraries constructed from sexually mature (162-day old) and immature (42-day old) ovary tissues of Single Comb White Leghorn chicken were sequenced using Illumina small RNA deep sequencing.

Results

In the present study, 14,545,100 and 14,774,864 clean reads were obtained from sexually mature (162-d) and sexually immature (42-d) ovaries, respectively. In total, 202 known miRNAs were identified, and 93 of them were found to be significantly differentially expressed: 42 miRNAs were up-regulated and 51 miRNAs were down-regulated in the mature ovary compared to the immature ovary. Among the up-regulated miRNAs, gga-miR-1a has the largest fold-change (6.405-fold), while gga-miR-375 has the largest fold-change (11.345-fold) among the down-regulated miRNAs. The three most abundant miRNAs in the chicken ovary are gga-miR-10a, gga-let-7 and gga-miR-21. Five differentially expressed miRNAs (gga-miR-1a, 21, 26a, 137 and 375) were validated by real-time quantitative RT-PCR (qRT-PCR). Furthermore, the expression patterns of the five miRNAs were analyzed in different developmental stages of chicken ovary and follicles of various sizes.

Conclusion

The present study provides the first miRNA profile in sexually immature and mature chicken ovaries. Some miRNAs such as gga-miR-1a and gga-miR-21are expressed differentially in immature and mature chicken ovaries as well as among different sized follicles, suggesting an important role in the follicular growth or ovulation mechanism in the chicken.

Complexity of microRNA function and the role of isomiRs in lipid homeostasis

MicroRNAs (miRNAs) are key posttranscriptional regulators of biological pathways that govern lipid metabolic phenotypes. Recent advances in high-throughput small RNA sequencing technology have revealed the complex and dynamic repertoire of miRNAs. Specifically, it has been demonstrated that a single genomic locus can give rise to multiple, functionally distinct miRNA isoforms (isomiR). There are several mechanisms by which isomiRs can be generated, including processing heterogeneity and posttranscriptional modifications, such as RNA editing, exonuclease-mediated nucleotide trimming, and/or nontemplated nucleotide addition (NTA). NTAs are dominant at the 3'-end of a miRNA, are most commonly uridylation or adenlyation events, and are catalyzed by one or more of several nucleotidyl transferase enzymes. 3' NTAs can affect miRNA stability and/or activity and are physiologically regulated, whereas modifications to the 5'-ends of miRNAs likely alter miRNA targeting activity. Recent evidence also suggests that the biogenesis of specific miRNAs, or small RNAs that act as miRNAs, can occur through unconventional mechanisms that circumvent key canonical miRNA processing steps. The unveiling of miRNA diversity has significantly added to our view of the complexity of miRNA function. In this review we present the current understanding of the biological relevance of isomiRs and their potential role in regulating lipid metabolism.

Characterization and differential expression of microRNAs in the ovaries of pregnant and non-pregnant goats (Capra hircus)

Background

Ovarian follicular development and hormone secretion are complex and coordinated biological processes which will usually be altered during pregnancy. Ovarian function is tightly regulated by a multitude of genes, and also by some specific miRNAs. It is necessary to identify the differentially expressed miRNAs in the ovaries of pregnant and non-pregnant mammals, in order to further understand the role of miRNA-mediated post-transcriptional regulation in mammalian reproduction. Here, we performed a comprehensive search for hircine miRNAs using two small RNA sequencing libraries prepared from the ovaries of pregnant and non-pregnant goats.

Results

617 conserved and 7 putative novel miRNAs were identified in the hircine ovaries. A total of 471 conserved miRNAs (76.34%) were co-expressed in both pregnant and non-pregnant libraries, and 90 pregnancy-specific and 56 non-pregnancy-specific conserved miRNAs were identified. Additionally, 407 unique miRNAs (65.96%) were significantly differentially expressed in the pregnant and non-pregnant libraries, of which 294 were upregulated and 113 were downregulated in the pregnant library compared to the non-pregnant library. Further analysis showed that miR-143 was predicted to bind to the target sequences of Frizzled-6 and -3 receptor genes in the Wnt/beta-catenin signaling pathway, and let-7b may target the Activin receptor I and Smad 2/3 genes in the TGF-beta signaling pathway. The expression level of 5 randomly selected miRNAs were analyzed by quantitative real-time PCR (q-PCR), and the results demonstrated that the expression patterns were consistent with the Solexa sequencing results.

Conclusions

The identification and characterization of differentially expressed miRNAs in the ovaries of pregnant and non-pregnant goats provides important information on the role of miRNA in the regulation of the ovarian development and function. This data will be helpful to facilitate studies on the regulation of miRNAs during mammalian reproduction.

Comprehensive analysis of microRNAs in breast cancer

Background

MicroRNAs (miRNAs) are short noncoding RNAs (approximately 22 nucleotides in length) that play important roles in breast cancer progression by downregulating gene expression. The detailed mechanisms and biological functions of miRNA molecules in breast carcinogenesis have yet to be fully elucidated. This study used bioinformatics and experimental approaches to conduct detailed analysis of the dysregulated miRNAs, arm selection preferences, 3' end modifications, and position shifts in isoforms of miRNAs (isomiRs) in breast cancer.

Methods

Next-generation sequencing (NGS) data on breast cancer was obtained from the NCBI Sequence Read Archive (SRA). The miRNA expression profiles and isomiRs in normal breast and breast tumor tissues were determined by mapping the clean reads back to human miRNAs. Differences in miRNA expression and pre-miRNA 5p/3p arm usage between normal and breast tumor tissues were further investigated using stem-loop reverse transcription and real-time polymerase chain reaction.

Results

The analysis identified and confirmed the aberrant expression of 22 miRNAs in breast cancer. Results from pathway enrichment analysis further indicated that the aberrantly expressed miRNAs play important roles in breast carcinogenesis by regulating the mitogen-activated protein kinase (MAPK) signaling pathway. Data also indicated that the position shifts in isomiRs and 3' end modifications were consistent in breast tumor and adjacent normal tissues, and that 5p/3p arm usage of some miRNAs displayed significant preferences in breast cancer.

Conclusions

Expression pattern and arm selection of miRNAs are significantly varied in breast cancers through analyzing NGS data and experimental approach. These miRNA candidates have high potential to play critical roles in the progression of breast cancer and could potentially provide as targets for future therapy.

MicroRNA 3' end nucleotide modification patterns and arm selection preference in liver tissues

Background

The expression of microRNA (miRNA) genes undergoes several maturation steps. Recent studies brought new insights into the maturation process, but also raised debates on the maturation mechanism. To understand the mechanism better, we downloaded small RNA sequence reads from NCBI SRA and quantified the expression profiles of miRNAs in normal and tumor liver tissues.

Results

From these miRNA expression profiles, we studied several issues related to miRNA biogenesis. First of all, the 3' ends of mature miRNAs usually carried modified nucleotides, generated from nucleotide addition or RNA editing. We found that adenine accounted for more than 50% of all miRNA 3' end modification events in all libraries. However, uracil dominated over adenine in several miRNA types. Moreover, the miRNA reads in the HBV-associated libraries have much lower rates of nucleotide modification. These results indicate that miRNA 3' end modifications are miRNA specific and may differ between normal and tumor tissues. Secondly, according to the hydrogen-bonding theory, the expression ratio of 5p arm to 3p arm miRNAs, derived from the same pre-miRNA, should be constant over tissues. However, a comparison of the expression profiles of the 5p arm and 3p arm miRNAs showed that one arm is preferred in the normal liver tissue whereas the other is preferred in the tumor liver tissue. In other words, different liver tissues have their own preferences on selecting either arm to be mature miRNAs.

Conclusions

The results suggest that besides the traditional miRNA biogenesis theory, another mechanism may also participate in the miRNA biogenesis pathways.

The repertoire and features of human platelet microRNAs

Playing a central role in the maintenance of hemostasis as well as in thrombotic disorders, platelets contain a relatively diverse messenger RNA (mRNA) transcriptome as well as functional mRNA-regulatory microRNAs, suggesting that platelet mRNAs may be regulated by microRNAs. Here, we elucidated the complete repertoire and features of human platelet microRNAs by high-throughput sequencing. More than 492 different mature microRNAs were detected in human platelets, whereas the list of known human microRNAs was expanded further by the discovery of 40 novel microRNA sequences. As in nucleated cells, platelet microRNAs bear signs of post-transcriptional modifications, mainly terminal adenylation and uridylation. In vitro enzymatic assays demonstrated the ability of human platelets to uridylate microRNAs, which correlated with the presence of the uridyltransferase enzyme TUT4. We also detected numerous microRNA isoforms (isomiRs) resulting from imprecise Drosha and/or Dicer processing, in some cases more frequently than the reference microRNA sequence, including 5′ shifted isomiRs with redirected mRNA targeting abilities. This study unveils the existence of a relatively diverse and complex microRNA repertoire in human platelets, and represents a mandatory step towards elucidating the intraplatelet and extraplatelet role, function and importance of platelet microRNAs.

YM500: a small RNA sequencing (smRNA-seq) database for microRNA research

MicroRNAs (miRNAs) are small RNAs ∼22 nt in length that are involved in the regulation of a variety of physiological and pathological processes. Advances in high-throughput small RNA sequencing (smRNA-seq), one of the next-generation sequencing applications, have reshaped the miRNA research landscape. In this study, we established an integrative database, the YM500 (http://ngs.ym.edu.tw/ym500/), containing analysis pipelines and analysis results for 609 human and mice smRNA-seq results, including public data from the Gene Expression Omnibus (GEO) and some private sources. YM500 collects analysis results for miRNA quantification, for isomiR identification (incl. RNA editing), for arm switching discovery, and, more importantly, for novel miRNA predictions. Wetlab validation on >100 miRNAs confirmed high correlation between miRNA profiling and RT-qPCR results (R = 0.84). This database allows researchers to search these four different types of analysis results via our interactive web interface. YM500 allows researchers to define the criteria of isomiRs, and also integrates the information of dbSNP to help researchers distinguish isomiRs from SNPs. A user-friendly interface is provided to integrate miRNA-related information and existing evidence from hundreds of sequencing datasets. The identified novel miRNAs and isomiRs hold the potential for both basic research and biotech applications.

miRNA-miRNA interaction implicates for potential mutual regulatory pattern

BACKGROUND

Natural or endogenous sense/antisense miRNAs, located on sense and antisense strands in the same genomic region, respectively, are detected recently. However, little is known about these miRNA pairs, especially for their distributions in different animal species. We herein present systematic analysis of them in human, mouse and rat miRNAs, and their expression patterns based on deep sequencing datasets.

METHODS AND RESULTS

The phenomenon of miRNA-miRNA interaction could be detected in different animal species. The common miRNAs pairs were found across species. These miRNA pairs could form miRNA:miRNA duplex with complete complementary structure, and were prone to be located on specific chromosomes. They might be homologous miRNA genes (especially in human), or clustered in a gene cluster (especially in rat), or simultaneously detected in different genomic regions due to multicopy pre-miRNAs. Remarkably, some miRNA pairs, located in different genomic regions, also showed complementarity as well as endogenous sense/antisense miRNAs. Based on published deep sequencing datasets, one member of miRNA pairs always was abundantly expressed, whereas another was quite rare. Rare common target mRNAs of these miRNA pairs were predicted.

CONCLUSIONS

Interaction between miRNAs and significant expression divergence implied complex potential mutual regulatory pattern in the miRNA world. The study would enrich miRNA regulatory network.

[Application of next generation sequencing in microRNA detection]

MicroRNAs (miRNAs) are a class of ~22nt long non-coding RNAs. They are evolutionarily conserved and play essential roles in the regulation of post-transcriptional gene expression. The rapidly developing next generation sequencing (NGS) has important applications in miRNA detection. This review is focused on the mechanism of three NGS platforms and their applications in miRNA detection. In contrast to traditional methods, NGS has major advantages: high throughput, precise, accurate, and repeatable. Its application includes new miRNAs exploration, detection of miRNA*, miRNA editing, and isomiR and target mRNA detection. As NGS develops, the cost of sequencing is declining which makes it possible for NGS to be widely used in the coming years. Next generation sequencing will greatly promote researches of miRNAs.

Transcriptome analysis of microRNAs in developing cerebral cortex of rat

BACKGROUND

The morphogenesis of the cerebral cortex depends on the precise control of gene expression during development. Small non-coding RNAs, including microRNAs and other groups of small RNAs, play profound roles in various physiological and pathological processes via their regulation of gene expression. A systematic analysis of the expression profile of small non-coding RNAs in developing cortical tissues is important for clarifying the gene regulation networks mediating key developmental events during cortical morphogenesis.

RESULTS

Global profiling of the small RNA transcriptome was carried out in rat cerebral cortex from E10 till P28 using next-generation sequencing technique. We found an extraordinary degree of developmental stage-specific expression of a large group of microRNAs. A group of novel microRNAs with functional hints were identified, and brain-enriched expression and Dicer-dependent production of high-abundant novel microRNAs were validated. Profound editing of known microRNAs at "seed" sequence and flanking sequence was observed, with much higher editing events detected at late postnatal stages than embryonic stages, suggesting the necessity of microRNA editing for the fine tuning of gene expression during the formation of complicated synaptic connections at postnatal stages.

CONCLUSION

Our analysis reveals extensive regulation of microRNAs during cortical development. The dataset described here will be a valuable resource for clarifying new regulatory mechanisms for cortical development and diseases and will greatly contribute to our understanding of the divergence, modification, and function of microRNAs.

The microRNA Transcriptome of Human Cytomegalovirus (HCMV)

The purpose of the present study was to characterize the microRNA transcriptome (miRNAome) of the human cytomegalovirus (HCMV or HHV5). We used deep sequencing and real time PCR (qPCR) together with bioinformatics to analyze the pattern of small RNA expression in cells infected with low-passage isolates of HCMV as well as in plasma and amniotic fluid. We report here on the discovery of four new precursors and ten new miRNAs as well as eleven microRNA-offset-RNAs (moRs) that are all encoded by HCMV. About eighty percent of the total HCMV reads were perfectly mapped to HCMV miRNAs, strongly suggestive of their important biological role that in large remains still to be defined and characterized. Taken altogether, the results of this study demonstrate the power and usefulness of the combined bioinformatics/biological approach in discovering additional important members of HCMV- encoded small RNAs and can be applied to the study of other viruses as well.

An antagomir to microRNA Let7f promotes neuroprotection in an ischemic stroke model

We previously showed that middle-aged female rats sustain a larger infarct following experimental stroke as compared to younger female rats, and paradoxically, estrogen treatment to the older group is neurotoxic. Plasma and brain insulin-like growth factor-1 (IGF-1) levels decrease with age. However, IGF-1 infusion following stroke, prevents estrogen neurotoxicity in middle-aged female rats. IGF1 is neuroprotective and well tolerated, but also has potentially undesirable side effects. We hypothesized that microRNAs (miRNAs) that target the IGF-1 signaling family for translation repression could be alternatively suppressed to promote IGF-1-like neuroprotection. Here, we report that two conserved IGF pathway regulatory microRNAs, Let7f and miR1, can be inhibited to mimic and even extend the neuroprotection afforded by IGF-1. Anti-mir1 treatment, as late as 4 hours following ischemia, significantly reduced cortical infarct volume in adult female rats, while anti-Let7 robustly reduced both cortical and striatal infarcts, and preserved sensorimotor function and interhemispheric neural integration. No neuroprotection was observed in animals treated with a brain specific miRNA unrelated to IGF-1 (anti-miR124). Remarkably, anti-Let7f was only effective in intact females but not males or ovariectomized females indicating that the gonadal steroid environment critically modifies miRNA action. Let7f is preferentially expressed in microglia in the ischemic hemisphere and confirmed in ex vivo cultures of microglia obtained from the cortex. While IGF-1 was undetectable in microglia harvested from the non-ischemic hemisphere, IGF-1 was expressed by microglia obtained from the ischemic cortex and was further elevated by anti-Let7f treatment. Collectively these data support a novel miRNA-based therapeutic strategy for neuroprotection following stroke.

Comprehensive analysis of RNA-Seq data reveals extensive RNA editing in a human transcriptome

RNA editing is a post-transcriptional event that recodes hereditary information. Here we describe a comprehensive profile of the RNA editome of a male Han Chinese individual based on analysis of ∼767 million sequencing reads from poly(A)(+), poly(A)(-) and small RNA samples. We developed a computational pipeline that carefully controls for false positives while calling RNA editing events from genome and whole-transcriptome data of the same individual. We identified 22,688 RNA editing events in noncoding genes and introns, untranslated regions and coding sequences of protein-coding genes. Most changes (∼93%) converted A to I(G), consistent with known editing mechanisms based on adenosine deaminase acting on RNA (ADAR). We also found evidence of other types of nucleotide changes; however, these were validated at lower rates. We found 44 editing sites in microRNAs (miRNAs), suggesting a potential link between RNA editing and miRNA-mediated regulation. Our approach facilitates large-scale studies to profile and compare editomes across a wide range of samples.

miRNA arm selection and isomiR distribution in gastric cancer

BACKGROUND

MicroRNAs (miRNAs) are small non-protein-coding RNAs. miRNA genes need several biogenesis steps to form function miRNAs. However, the precise mechanism and biology involved in the mature miRNA molecules are not clearly investigated. In this study, we conducted in-depth analyses to examine the arm selection and isomiRs using NGS platform.

METHODS

We sequenced small RNAs from one pair of normal and gastric tumor tissues with Solexa platform. By analyzing the NGS data, we quantified the expression profiles of miRNAs and isomiRs in gastric tissues. Then, we measured the expression ratios of 5p arm to 3p arm of the same pre-miRNAs. And, we used Kolmogorov-Smirnov (KS) test to examine isomiR pattern difference between tissues.

RESULTS

Our result showed the 5p arm and 3p arm miRNA derived from the same pre-miRNAs have different tissue expression preference, one preferred normal tissue and the other preferred tumor tissue, which strongly implied that there could be other mechanism controlling mature miRNA selection in addition to the known hydrogen-bonding selection rule. Furthermore, by using the KS test, we demonstrated that some isomiR types preferentially occur in normal gastric tissue but other types prefer tumor gastric tissue.

CONCLUSIONS

Arm selections and isomiR patterns are significantly varied in human cancers by using deep sequencing NGS data. Our results provided a novel research topic in miRNA regulation study. With advanced bioinformatics and molecular biology studies, more robust conclusions and insight into miRNA regulation can be achieved in the near future.

MicroRNAs and their isomiRs function cooperatively to target common biological pathways

Background

Variants of microRNAs (miRNAs), called isomiRs, are commonly reported in deep-sequencing studies; however, the functional significance of these variants remains controversial. Observational studies show that isomiR patterns are non-random, hinting that these molecules could be regulated and therefore functional, although no conclusive biological role has been demonstrated for these molecules.

Results

To assess the biological relevance of isomiRs, we have performed ultra-deep miRNA-seq on ten adult human tissues, and created an analysis pipeline called miRNA-MATE to align, annotate, and analyze miRNAs and their isomiRs. We find that isomiRs share sequence and expression characteristics with canonical miRNAs, and are generally strongly correlated with canonical miRNA expression. A large proportion of isomiRs potentially derive from AGO2 cleavage independent of Dicer. We isolated polyribosome-associated mRNA, captured the mRNA-bound miRNAs, and found that isomiRs and canonical miRNAs are equally associated with translational machinery. Finally, we transfected cells with biotinylated RNA duplexes encoding isomiRs or their canonical counterparts and directly assayed their mRNA targets. These studies allow us to experimentally determine genome-wide mRNA targets, and these experiments showed substantial overlap in functional mRNA networks suppressed by both canonical miRNAs and their isomiRs.

Conclusions

Together, these results find isomiRs to be biologically relevant and functionally cooperative partners of canonical miRNAs that act coordinately to target pathways of functionally related genes. This work exposes the complexity of the miRNA-transcriptome, and helps explain a major miRNA paradox: how specific regulation of biological processes can occur when the specificity of miRNA targeting is mediated by only 6 to 11 nucleotides.

Accurate identification of A-to-I RNA editing in human by transcriptome sequencing

RNA editing enhances the diversity of gene products at the post-transcriptional level. Approaches for genome-wide identification of RNA editing face two main challenges: separating true editing sites from false discoveries and accurate estimation of editing levels. We developed an approach to analyze transcriptome sequencing data (RNA-seq) for global identification of RNA editing in cells for which whole-genome sequencing data are available. We applied the method to analyze RNA-seq data of a human glioblastoma cell line, U87MG. Around 10,000 DNA-RNA differences were identified, the majority being putative A-to-I editing sites. These predicted A-to-I events were associated with a low false-discovery rate (∼5%). Moreover, the estimated editing levels from RNA-seq correlated well with those based on traditional clonal sequencing. Our results further facilitated unbiased characterization of the sequence and evolutionary features flanking predicted A-to-I editing sites and discovery of a conserved RNA structural motif that may be functionally relevant to editing. Genes with predicted A-to-I editing were significantly enriched with those known to be involved in cancer, supporting the potential importance of cancer-specific RNA editing. A similar profile of DNA-RNA differences as in U87MG was predicted for another RNA-seq data set obtained from primary breast cancer samples. Remarkably, significant overlap exists between the putative editing sites of the two transcriptomes despite their difference in cell type, cancer type, and genomic backgrounds. Our approach enabled de novo identification of the RNA editome, which sets the stage for further mechanistic studies of this important step of post-transcriptional regulation.

MicroRNA identity and abundance in developing swine adipose tissue as determined by Solexa sequencing

MicroRNAs (miRNAs) are small ∼22-nt regulatory RNAs that regulate the stability and translation of cognate mRNAs. MiRNAs participate in the regulation of adipogenesis, and identification of the full repertoire of miRNAs expressed in adipose tissue is likely to significantly increase our understanding of adipose tissue growth and development. Here, we adopted a deep sequencing approach to determine the identity and abundance of miRNAs in developing swine adipose tissue. Via this approach, we identified the sequences and relative expression levels of 227 conserved miRNAs (of which 59 were novel) and 66 potential porcine miRNAs. The expression levels displayed a large range, as reflected by the number of sequence reads, which varied from several counts for rare miRNAs to several million reads for the most abundant miRNAs. The abundant miRNAs principally belonged to 32 miRNA gene families, including miR-143, miR-103, let-7, and miR-148. Of the conserved miRNAs, 93 miRNAs were up-regulated and 33 miRNAs were down-regulated in the adult pig adipose tissue. Moreover, we observed sequence variants and seed edits of the miRNAs. KEGG pathway analysis and GO term enrichment suggested that highly expressed miRNAs are involved in adipose tissue development, signal transduction, cell-cell and cell-extracellular matrix communication, neural development and function, and lipid metabolism including carboxylic acid, oxoacid, fatty acid, steroid, glycerolipid, alcohol and phospholipid metabolism. Our results expand the number of known porcine miRNAs and provide a thorough account of the miRNA transcriptome in porcine adipose tissue.

MicroRNA: Implications for Alzheimer Disease and other Human CNS Disorders

Understanding complex diseases such as sporadic Alzheimer disease (AD) has been a major challenge. Unlike the familial forms of AD, the genetic and environmental risks factors identified for sporadic AD are extensive. MicroRNAs are one of the major noncoding RNAs that function as negative regulators to silence or suppress gene expression via translational inhibition or message degradation. Their discovery has evoked great excitement in biomedical research for their promise as potential disease biomarkers and therapeutic targets. Key microRNAs have been identified as essential for a variety of cellular events including cell lineage determination, proliferation, apoptosis, DNA repair, and cytoskeletal organization; most, if not all, acting to fine-tune gene expression at the post-transcriptional level in a host of cellular signaling networks. Dysfunctional microRNA-mediated regulation has been implicated in the pathogenesis of many disease states. Here, the current understanding of the role of miRNAs in the central nervous system is reviewed with emphasis on their impact on the etiopathogenesis of sporadic AD.

A comprehensive survey of miRNA repertoire and 3' addition events in the placentas of patients with pre-eclampsia from high-throughput sequencing

Background

To gain insight into potential roles of isomiR spectrum and isomiRs with 3′ additions in pre-eclampsia, we performed a comprehensive survey of miRNA repertoire and 3′ addition events from placental samples with different degrees of pre-eclampsia by applying SOLiD sequencing platform.

Principal Findings

Over 30% isomiRs were detected with 3′ non-template additional nucleotides, especially for additional nucleotide of adenosine. However, these modified isomiRs showed a lower percentage of total miRNA expression (<15%). Generally, 1-3 abundant isomiRs from a given miRNA locus were identified, but none of them was detected with 3′ additions. Different miRNAs indicated various isomiR spectrums and expression patterns. The most abundant isomiR spectrum, isomiR profile and expression pattern always were stability, but herein we found several exceptions across samples, especially between normal and diseased samples. At isomiR level, we detected a distinct subset of differentially expressed modified isomiRs between normal and diseased samples or between mild and severe samples. Gene Ontology analysis of their experimentally validated target genes revealed enrichment for specific biological process categories.

Conclusions

The phenomenon of multiple isomiRs, especially for isomiRs with 3′ additions, is not a random event during pre-miRNA processing. Varieties of isomiRs and expression patterns reveal potential functional implication and should be taken into account. The study enriches association of miRNAs and human disease, including potential roles of various miRNA variants and 3′ addition events.

Song exposure regulates known and novel microRNAs in the zebra finch auditory forebrain

BACKGROUND

In an important model for neuroscience, songbirds learn to discriminate songs they hear during tape-recorded playbacks, as demonstrated by song-specific habituation of both behavioral and neurogenomic responses in the auditory forebrain. We hypothesized that microRNAs (miRNAs or miRs) may participate in the changing pattern of gene expression induced by song exposure. To test this, we used massively parallel Illumina sequencing to analyse small RNAs from auditory forebrain of adult zebra finches exposed to tape-recorded birdsong or silence.

RESULTS

In the auditory forebrain, we identified 121 known miRNAs conserved in other vertebrates. We also identified 34 novel miRNAs that do not align to human or chicken genomes. Five conserved miRNAs showed significant and consistent changes in copy number after song exposure across three biological replications of the song-silence comparison, with two increasing (tgu-miR-25, tgu-miR-192) and three decreasing (tgu-miR-92, tgu-miR-124, tgu-miR-129-5p). We also detected a locus on the Z sex chromosome that produces three different novel miRNAs, with supporting evidence from Northern blot and TaqMan qPCR assays for differential expression in males and females and in response to song playbacks. One of these, tgu-miR-2954-3p, is predicted (by TargetScan) to regulate eight song-responsive mRNAs that all have functions in cellular proliferation and neuronal differentiation.

CONCLUSIONS

The experience of hearing another bird singing alters the profile of miRNAs in the auditory forebrain of zebra finches. The response involves both known conserved miRNAs and novel miRNAs described so far only in the zebra finch, including a novel sex-linked, song-responsive miRNA. These results indicate that miRNAs are likely to contribute to the unique behavioural biology of learned song communication in songbirds.

Genome-wide approaches in the study of microRNA biology

MicroRNAs (miRNAs), a class of ∼21-23 nucleotide long non-coding RNAs (ncRNAs), have critical roles in diverse biological processes that encompass development, proliferation, apoptosis, stress response, and fat metabolism. miRNAs recognize their target mRNA transcripts by partial sequence complementarity and collectively have been estimated to regulate the majority of human genes. Consequently, misregulation of miRNAs or disruption of their target sites in genes has been implicated in a variety of human diseases ranging from cancer metastasis to neurological disorders. With the development and availability of genomic technologies and computational approaches, the field of miRNA biology has advanced tremendously over the last decade. Here we review the genome-wide approaches that have allowed for the discovery of new miRNAs, the characterization of their targets, and a systems-level view of their impact.