Many ways to generate microRNA-like small RNAs: non-canonical pathways for microRNA production

Large-scale identification of mirtrons in Arabidopsis and rice

A new catalog of microRNA (miRNA) species called mirtrons has been discovered in animals recently, which originate from spliced introns of the gene transcripts. However, only one putative mirtron, osa-MIR1429, has been identified in rice (Oryza sativa). We employed a high-throughput sequencing (HTS) data- and structure-based approach to do a genome-wide search for the mirtron candidate in both Arabidopsis (Arabidopsis thaliana) and rice. Five and eighteen candidates were discovered in the two plants respectively. To investigate their biological roles, the targets of these mirtrons were predicted and validated based on degradome sequencing data. The result indicates that the mirtrons could guide target cleavages to exert their regulatory roles post-transcriptionally, which needs further experimental validation.

Complexity of murine cardiomyocyte miRNA biogenesis, sequence variant expression and function

microRNAs (miRNAs) are critical to heart development and disease. Emerging research indicates that regulated precursor processing can give rise to an unexpected diversity of miRNA variants. We subjected small RNA from murine HL-1 cardiomyocyte cells to next generation sequencing to investigate the relevance of such diversity to cardiac biology. ∼40 million tags were mapped to known miRNA hairpin sequences as deposited in miRBase version 16, calling 403 generic miRNAs as appreciably expressed. Hairpin arm bias broadly agreed with miRBase annotation, although 44 miR* were unexpectedly abundant (>20% of tags); conversely, 33 -5p/-3p annotated hairpins were asymmetrically expressed. Overall, variability was infrequent at the 5′ start but common at the 3′ end of miRNAs (5.2% and 52.3% of tags, respectively). Nevertheless, 105 miRNAs showed marked 5′ isomiR expression (>20% of tags). Among these was miR-133a, a miRNA with important cardiac functions, and we demonstrated differential mRNA targeting by two of its prevalent 5′ isomiRs. Analyses of miRNA termini and base-pairing patterns around Drosha and Dicer cleavage regions confirmed the known bias towards uridine at the 5′ most position of miRNAs, as well as supporting the thermodynamic asymmetry rule for miRNA strand selection and a role for local structural distortions in fine tuning miRNA processing. We further recorded appreciable expression of 5 novel miR*, 38 extreme variants and 8 antisense miRNAs. Analysis of genome-mapped tags revealed 147 novel candidate miRNAs. In summary, we revealed pronounced sequence diversity among cardiomyocyte miRNAs, knowledge of which will underpin future research into the mechanisms involved in miRNA biogenesis and, importantly, cardiac function, disease and therapy.

Identification of conserved and novel microRNAs that are responsive to heat stress in Brassica rapa

The species Brassica rapa includes various vegetable crops. Production of these vegetable crops is usually impaired by heat stress. Some microRNAs (miRNAs) in Arabidopsis have been considered to mediate gene silencing in plant response to abiotic stress. However, it remains unknown whether or what miRNAs play a role in heat resistance of B. rapa. To identify genomewide conserved and novel miRNAs that are responsive to heat stress in B. rapa, we defined temperature thresholds of non-heading Chinese cabbage (B. rapa ssp. chinensis) and constructed small RNA libraries from the seedlings that had been exposed to high temperature (46 °C) for 1 h. By deep sequencing and data analysis, we selected a series of conserved and novel miRNAs that responded to heat stress. In total, Chinese cabbage shares at least 35 conserved miRNA families with Arabidopsis thaliana. Among them, five miRNA families were responsive to heat stress. Northern hybridization and real-time PCR showed that the conserved miRNAs bra-miR398a and bra-miR398b were heat-inhibitive and guided heat response of their target gene, BracCSD1; and bra-miR156h and bra-miR156g were heat-induced and its putative target BracSPL2 was down-regulated. According to the criteria of miRNA and miRNA* that form a duplex, 21 novel miRNAs belonging to 19 miRNA families were predicted. Of these, four were identified to be heat-responsive by Northern blotting and/or expression analysis of the putative targets. The two novel miRNAs bra-miR1885b.3 and bra-miR5718 negatively regulated their putative target genes. 5'-Rapid amplification of cDNA ends PCR indicated that three novel miRNAs cleaved the transcripts of their target genes where their precursors may have evolved from. These results broaden our perspective on the important role of miRNA in plant responses to heat.

Testis-specific miRNA-469 up-regulated in gonadotropin-regulated testicular RNA helicase (GRTH/DDX25)-null mice silences transition protein 2 and protamine 2 messages at sites within coding region: implications of its role in germ cell development

Gonadotropin-regulated testicular RNA helicase (GRTH/DDX25), a testis-specific member of the DEAD-box family, is an essential post-transcriptional regulator of spermatogenesis. Failure of expression of Transition protein 2 (TP2) and Protamine 2 (Prm2) proteins (chromatin remodelers, essential for spermatid elongation and completion of spermatogenesis) with preservation of their mRNA expression was observed in GRTH-null mice (azoospermic due to failure of spermatids to elongate). These were identified as target genes for the testis-specific miR-469, which is increased in the GRTH-null mice. Further analysis demonstrated that miR-469 repressed TP2 and Prm2 protein expression at the translation level with minor effect on mRNA degradation, through binding to the coding regions of TP2 and Prm2 mRNAs. The corresponding primary-microRNAs and the expression levels of Drosha and DGCR8 (both mRNA and protein) were increased significantly in the GRTH-null mice. miR-469 silencing of TP2 and Prm2 mRNA in pachytene spermatocytes and round spermatids is essential for their timely translation at later times of spermiogenesis, which is critical to attain mature sperm. Collectively, these studies indicate that GRTH, a multifunctional RNA helicase, acts as a negative regulator of miRNA-469 biogenesis and consequently their function during spermatogenesis.

Prognostic values of the miR-17-92 cluster and its paralogs in colon cancer

BACKGROUND

MicroRNAs have been shown to offer great potential in both the diagnosis and prognosis of cancer. Despite the well-established role of the miR-17-92 in cancer formation and progression, the contribution of each individual miRNA remains to be characterized. Thus, we investigated whether deregulation of the miR-17-92 associated with colon cancer prognosis.

METHODS

Expression levels of the miR-17-92 cluster and its paralogs were determined in 48 colon tumor and 48 paired normal tissues by real-time qRT-PCR. Associations with miRNA expression, age, sex, TNM staging, and survival prognosis were evaluated.

RESULTS

MiR-17-92 cluster and its paralogs were significantly overexpressed in colon tumor. No significant associations were found between the deregulation of certain miRNAs and the clinical and pathologic characteristics observed in patients. Kaplan-Meier curves demonstrated significantly reduced overall survival in patients expressing high levels of miR-17. In multivariate Cox models, miR-17 overexpression (HR 2.67; P = 0.007) and TNM staging (HR 8.87; P = 0.002) were significantly associated with a risk of death.

CONCLUSIONS

The miR-17-92 cluster and its paralogs were significantly elevated in patients with colon cancer, and heightened expression of miR-17 was associated with poor survival. Moreover, miR-17 and TNM staging were both identified as significant, but independent, prognostic biomarkers in colon cancer.

deepBlockAlign: a tool for aligning RNA-seq profiles of read block patterns

MOTIVATION

High-throughput sequencing methods allow whole transcriptomes to be sequenced fast and cost-effectively. Short RNA sequencing provides not only quantitative expression data but also an opportunity to identify novel coding and non-coding RNAs. Many long transcripts undergo post-transcriptional processing that generates short RNA sequence fragments. Mapped back to a reference genome, they form distinctive patterns that convey information on both the structure of the parent transcript and the modalities of its processing. The miR-miR* pattern from microRNA precursors is the best-known, but by no means singular, example.

RESULTS

deepBlockAlign introduces a two-step approach to align RNA-seq read patterns with the aim of quickly identifying RNAs that share similar processing footprints. Overlapping mapped reads are first merged to blocks and then closely spaced blocks are combined to block groups, each representing a locus of expression. In order to compare block groups, the constituent blocks are first compared using a modified sequence alignment algorithm to determine similarity scores for pairs of blocks. In the second stage, block patterns are compared by means of a modified Sankoff algorithm that takes both block similarities and similarities of pattern of distances within the block groups into account. Hierarchical clustering of block groups clearly separates most miRNA and tRNA, and also identifies about a dozen tRNAs clustering together with miRNA. Most of these putative Dicer-processed tRNAs, including eight cases reported to generate products with miRNA-like features in literature, exhibit read blocks distinguished by precise start position of reads.

AVAILABILITY

The program deepBlockAlign is available as source code from http://rth.dk/resources/dba/.

CONTACT

gorodkin@rth.dk; studla@bioinf.uni-leipzig.de

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Posttranslational modification of Argonautes and their role in small RNA-mediated gene regulation

Shortly after their discovery, repertoires of miRNA were identified, together with proteins involved in their biogenesis and action. It is now obvious that miRNA-mediated gene regulation itself is regulated at multiple levels. Identifying the regulatory mechanisms that underpin small RNA homeostasis by modulation of their biogenesis and action has become a key issue, which can be partly resolved by identifying mediators of Argonautes turnover. An emerging theme in the control of Argonaute stability and activity is through posttranslational modifications, which are the focus of this review.

Gatekeepers for Piwi-piRNA complexes to enter the nucleus

RNA silencing pathways are now recognized to participate in essential cellular functions ranging from the regulation of mRNA turnover to the suppression of the activity of potentially deleterious transposable elements (TEs). Piwi-interacting RNAs (piRNAs) are germline-specific, small silencing RNAs that suppress TE activity and maintain genome integrity during germline development. In Drosophila ovarian somatic cells, piRNAs are processed from long single-stranded RNAs by a Dicer-independent pathway and are loaded onto Piwi in the cytoplasm. The Piwi-piRNA complexes are then transported into the nucleus to exert TE silencing. This mechanism involves gatekeepers for a functional Piwi-piRNA complex to be imported, which parallels with the Tetrahymena Twi1p-scan RNA pathway used to carry out the programmed DNA elimination.

Are the Ro RNP-associated Y RNAs concealing microRNAs? Y RNA-derived miRNAs may be involved in autoimmunity

Here we discuss the hypothesis that the RNA components of the Ro ribonucleoproteins (RNPs), the Y RNAs, can be processed into microRNAs (miRNAs). Although Ro RNPs, whose main protein components Ro60 and La are targeted by the immune system in several autoimmune diseases, were discovered many years ago, their function is still poorly understood. Indeed, recent data show that miRNA-sized small RNAs can be generated from Y RNAs. This hypothesis leads also to a model in which Ro60 acts as a modulator in the Y RNA-derived miRNA biogenesis pathway. The implications of these Y RNA-derived miRNAs, which may be specifically produced under pathological circumstances such as in autoimmunity or during viral infections, for the enigmatic function of Ro RNPs are discussed.

RNA in pieces

Eukaryotic genomes accommodate numerous types of information within diverse DNA and RNA sequence elements. At many loci, these elements overlap and the same sequence is read multiple times during the production, processing, localization, function and turnover of a single transcript. Moreover, two or more transcripts from the same locus might use a common sequence in different ways, to perform distinct biological roles. Recent results show that many transcripts also undergo post-transcriptional cleavage to release specific fragments, which can then function independently. This phenomenon appears remarkably widespread, with even well-documented transcript classes such as messenger RNAs yielding fragments. RNA fragmentation significantly expands the already extraordinary spectrum of transcripts present within eukaryotic cells, and also calls into question how the 'gene' should be defined.

Mechanisms and role of microRNA deregulation in cancer onset and progression

MicroRNAs are key regulators of various fundamental biological processes and, although representing only a small portion of the genome, they regulate a much larger population of target genes. Mature microRNAs (miRNAs) are single-stranded RNA molecules of 20-23 nucleotide (nt) length that control gene expression in many cellular processes. These molecules typically reduce the stability of mRNAs, including those of genes that mediate processes in tumorigenesis, such as inflammation, cell cycle regulation, stress response, differentiation, apoptosis and invasion. MicroRNA targeting is mostly achieved through specific base-pairing interactions between the 5' end ('seed' region) of the miRNA and sites within coding and untranslated regions (UTRs) of mRNAs; target sites in the 3' UTR diminish mRNA stability. Since miRNAs frequently target hundreds of mRNAs, miRNA regulatory pathways are complex. Calin and Croce were the first to demonstrate a connection between microRNAs and increased risk of developing cancer, and meanwhile the role of microRNAs in carcinogenesis has definitively been evidenced. It needs to be considered that the complex mechanism of gene regulation by microRNAs is profoundly influenced by variation in gene sequence (polymorphisms) of the target sites. Thus, individual variability could cause patients to present differential risks regarding several diseases. Aiming to provide a critical overview of miRNA dysregulation in cancer, this article reviews the growing number of studies that have shown the importance of these small molecules and how these microRNAs can affect or be affected by genetic and epigenetic mechanisms.

Mirtrons: microRNA biogenesis via splicing

A well-defined mechanism governs the maturation of most microRNAs (miRNAs) in animals, via stepwise cleavage of precursor hairpin transcripts by the Drosha and Dicer RNase III enzymes. Recently, several alternative miRNA biogenesis pathways were elucidated, the most prominent of which substitutes Drosha cleavage with splicing. Such short hairpin introns are known as mirtrons, and their study has uncovered related pathways that combine splicing with other ribonucleolytic machinery to yield Dicer substrates for miRNA biogenesis. In this review, we consider the mechanisms of splicing-mediated miRNA biogenesis, computational strategies for mirtron discovery, and the evolutionary implications of the existence of multiple miRNA biogenesis pathways. Altogether, the features of mirtron pathways illustrate unexpected flexibility in combining RNA processing pathways, and highlight how multiple functions can be encoded by individual transcripts.

A novel class of heat-responsive small RNAs derived from the chloroplast genome of Chinese cabbage (Brassica rapa)

BACKGROUND

Non-coding small RNAs play critical roles in various cellular processes in a wide spectrum of eukaryotic organisms. Their responses to abiotic stress have become a popular topic of economic and scientific importance in biological research. Several studies in recent years have reported a small number of non-coding small RNAs that map to chloroplast genomes. However, it remains uncertain whether small RNAs are generated from chloroplast genome and how they respond to environmental stress, such as high temperature. Chinese cabbage is an important vegetable crop, and heat stress usually causes great losses in yields and quality. Under heat stress, the leaves become etiolated due to the disruption and disassembly of chloroplasts. In an attempt to determine the heat-responsive small RNAs in chloroplast genome of Chinese cabbage, we carried out deep sequencing, using heat-treated samples, and analysed the proportion of small RNAs that were matched to chloroplast genome.

RESULTS

Deep sequencing provided evidence that a novel subset of small RNAs were derived from the chloroplast genome of Chinese cabbage. The chloroplast small RNAs (csRNAs) include those derived from mRNA, rRNA, tRNA and intergenic RNA. The rRNA-derived csRNAs were preferentially located at the 3'-ends of the rRNAs, while the tRNA-derived csRNAs were mainly located at 5'-termini of the tRNAs. After heat treatment, the abundance of csRNAs decreased in seedlings, except those of 24 nt in length. The novel heat-responsive csRNAs and their locations in the chloroplast were verified by Northern blotting. The regulation of some csRNAs to the putative target genes were identified by real-time PCR. Our results reveal that high temperature suppresses the production of some csRNAs, which have potential roles in transcriptional or post-transcriptional regulation.

CONCLUSIONS

In addition to nucleus, the chloroplast is another important organelle that generates a number of small RNAs. Many members of csRNA families are highly sensitive to heat stress. Some csRNAs respond to heat stress by silencing target genes. We suggest that proper temperature is important for production of chloroplast small RNAs, which are associated with plant resistance to abiotic stress.

From snoRNA to miRNA: Dual function regulatory non-coding RNAs

Small nucleolar RNAs (snoRNAs) are an ancient class of small non-coding RNAs present in all eukaryotes and a subset of archaea that carry out a fundamental role in the modification and processing of ribosomal RNA. In recent years, however, a large proportion of snoRNAs have been found to be further processed into smaller molecules, some of which display different functionality. In parallel, several studies have uncovered extensive similarities between snoRNAs and other types of small non-coding RNAs, and in particular microRNAs. Here, we explore the extent of the relationship between these types of non-coding RNA and the possible underlying evolutionary forces that shaped this subset of the current non-coding RNA landscape.

Precursor miR-886, a novel noncoding RNA repressed in cancer, associates with PKR and modulates its activity

Noncoding RNAs have drawn significant attention in biology recently. Whereas the current research is highly inclined to microRNAs, research on other noncoding RNAs has lagged behind. Here, we investigated a novel noncoding RNA that has been known as precursor microRNA miR-886 (pre-miR-886). Pre-miR-886 has been proposed also as a vault RNA, a component of the vault complex implicated in cancer drug resistance. We identified pre-miR-886 as a 102-nucleotide-long, abundant cytoplasmic RNA that is neither a genuine pre-microRNA nor a vault RNA. Pre-miR-886 is physically associated with PKR (Protein Kinase RNA-activated), an interferon-inducible and double-stranded RNA dependent kinase. The suppression of pre-miR-886 activates PKR and its downstream pathways, eIF2α phosphorylation and the NF-κB pathway, leading to impaired cell proliferation. We also found that pre-miR-886 is suppressed in a wide-range of cancer cell lines and in clinical specimens. This study is the first intense characterization of pre-miR-886 as well as the initial report on its function as a PKR regulator, which suggests a critical role in tumorigenesis.

The role of the precursor structure in the biogenesis of microRNA

The human genome contains more than 1,000 microRNA (miRNA) genes, which are transcribed mainly by RNA polymerase II. The canonical pathway of miRNA biogenesis includes the nuclear processing of primary transcripts (pri-miRNAs) by the ribonuclease Drosha and further cytoplasmic processing of pre-miRNAs by the ribonuclease Dicer. This review discusses the issue of miRNA end heterogeneity generated primarily by Drosha and Dicer cleavage and focuses on the structural aspects of the Dicer step of miRNA biogenesis. We examine the structures of miRNA precursors, both predicted and experimentally determined, as well as the influence of various motifs that disturb the regularity of pre-miRNA structure on Dicer cleavage specificity. We evaluate the structural determinants of the length diversity of miRNA generated by Dicer from different precursors and highlight the importance of asymmetrical motifs. Finally, we discuss the impact of Dicer protein partners on cleavage efficiency and specificity and propose the contribution of pre-miRNA structural plasticity to the dynamics of the dicing complex.

Vive la différence: biogenesis and evolution of microRNAs in plants and animals

MicroRNAs are pervasive in both plants and animals, but many aspects of their biogenesis, function and evolution differ. We reveal how these differences contribute to characteristic features of microRNA evolution in the two kingdoms.

Vive la différence: biogenesis and evolution of microRNAs in plants and animals

MicroRNAs are pervasive in both plants and animals, but many aspects of their biogenesis, function and evolution differ. We reveal how these differences contribute to characteristic features of microRNA evolution in the two kingdoms.

Characterization of human platelet microRNA by quantitative PCR coupled with an annotation network for predicted target genes

Platelets are anucleate blood cells that play a crucial role in thrombosis and hemostasis. Despite their lack of nuclear DNA, platelets contain significant amounts of microRNA (miRNA) that may have vital functions in post-transcriptional gene regulation. Here, we combined comprehensive miRNA expression profiling by quantitative PCR with target prediction analysis for the most abundant miRNAs in human platelets. A network composed of predicted platelet miRNA target genes was then constructed, using annotations available in Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. In addition, we evaluated possible differences in miRNA levels between resting and thrombin-stimulated platelets. We identified 281 transcripts, including 228 mature miRNAs and 53 minor miRNAs (or miR*), of which six miRNAs (miR-15 a, miR-339-3 p, miR-365, miR-495, miR-98, and miR-361-3 p) were up- or down-regulated in activated human platelets (P ≤ 0.001). A redundancy-reduced network was established that encompassed 246 genes in five statistically significant functional clusters representing platelet miRNA regulating pathways. Comparison of the 246 network genes with the platelet mRNA expression data available at ArrayExpress database confirmed that most of these genes (89%) are expressed in human platelets. In conclusion, this work affirms a recent microarray study reporting a wide-spread existence of miRNAs in human platelets. Further, we observed that thrombin stimulation was associated with altered levels of some miRNAs in platelets. The proposed functional network, combining computational prediction analysis with annotations from experimental observations, may in addition provide some information about probable miRNA target pathways in human platelets.

Identification of human miRNA precursors that resemble box C/D snoRNAs

There are two main classes of small nucleolar RNAs (snoRNAs): the box C/D snoRNAs and the box H/ACA snoRNAs that function as guide RNAs to direct sequence-specific modification of rRNA precursors and other nucleolar RNA targets. A previous computational and biochemical analysis revealed a possible evolutionary relationship between miRNA precursors and some box H/ACA snoRNAs. Here, we investigate a similar evolutionary relationship between a subset of miRNA precursors and box C/D snoRNAs. Computational analyses identified 84 intronic miRNAs that are encoded within either box C/D snoRNAs, or in precursors showing similarity to box C/D snoRNAs. Predictions of the folded structures of these box C/D snoRNA-like miRNA precursors resemble the structures of known box C/D snoRNAs, with the boxes C and D often in close proximity in the folded molecule. All five box C/D snoRNA-like miRNA precursors tested (miR-27b, miR-16-1, mir-28, miR-31 and let-7g) bind to fibrillarin, a specific protein component of functional box C/D snoRNP complexes. The data suggest that a subset of small regulatory RNAs may have evolved from box C/D snoRNAs.

Dicer-independent, Ago2-mediated microRNA biogenesis in vertebrates

A canonical biogenesis pathway involving sequential cleavage by the Drosha and Dicer RNAse III enzymes governs the maturation of most animal microRNAs. However, there exist a variety of alternative miRNA biogenesis pathways, most of which bypass Drosha processing. Recently, three groups described for the first time a vertebrate microRNA pathway that bypasses Dicer cleavage. This mechanism was characterized with respect to the highly conserved vertebrate gene mir-451, for which Drosha processing yields a short (42 nucleotide) hairpin that is directly loaded into Ago2, the sole vertebrate "Slicer" Argonaute. Ago2-mediated cleavage of this hairpin yields a 30 nucleotide intermediate, whose 3' end is resected to generate the dominantly cloned ~23 nucleotide mature miR-451. Knowledge of this pathway provides an unprecedented tool with which to express microRNAs and small interfering RNAs in Dicer mutant cells. More generally, the mir-451 backbone constitutes a new platform for gene silencing that complements existing shRNA technology.