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

MicroRNA (miRNA) 5′-isoforms, or 5′-isomiRs, are small-RNA species that originate from the same genomic loci as the major miRNAs with their 5′ ends shifted from the 5′ ends of the miRNAs by a few nucleotides. Although 5′-isomiRs have been reported, their origins, properties and potential functions remain to be examined. We systematically studied 5′-isomiRs in human, mouse, fruitfly and worm by analysing a large collection of small non-coding RNA and mRNA profiling data. The results revealed a broad existence of 5′-isomiRs in the four species, many of which were conserved and could arise from genomic loci of canonical and non-canonical miRNAs. The well-conserved 5′-isomiRs have several features, including a preference of the 3p over the 5p arms of hairpins of conserved mammalian miRNAs, altered 5′-isomiRs across species and across tissues, and association with structural variations of miRNA hairpins. Importantly, 5′-isomiRs and their major miRNAs may have different mRNA targets and thus potentially play distinct roles of gene regulation, as shown by an integrative analysis combining miRNA and mRNA profiling data from psoriatic and normal human skin and from murine miRNA knockout assays. Indeed, 18 5′-isomiRs had aberrant expression in psoriatic human skin, suggesting their potential function in psoriasis pathogenesis. The results of the current study deepened our understanding of the diversity and conservation of miRNAs, their plasticity in gene regulation and potential broad function in complex diseases.     

Large-Scale Screens of miRNA-mRNA Interactions Unveiled That the 3′UTR of a Gene Is Targeted by Multiple miRNAs

Animal microRNA (miRNA) target prediction is still a challenge, although many prediction programs have been exploited. MiRNAs exert their function through partially binding the messenger RNAs (mRNAs; likely at 3′ untranslated regions [3′UTRs]), which makes it possible to detect the miRNA-mRNA interactions in vitro by co-transfection of miRNA and a luciferase reporter gene containing the target mRNA fragment into mammalian cells under a dual-luciferase assay system. Here, we constructed a human miRNA expression library and used a dual-luciferase assay system to perform large-scale screens of interactions between miRNAs and the 3′UTRs of seven genes, which included more than 3,000 interactions with triplicate experiments for each interaction. The screening results showed that the 3′UTR of one gene can be targeted by multiple miRNAs. Among the prediction algorithms, a Bayesian phylogenetic miRNA target identification algorithm and a support vector machine (SVM) presented a relatively better performance (27% for EIMMo and 24.7% for miRDB) against the average precision (17.3%) of the nine prediction programs used here. Additionally, we noticed that a relatively high conservation level was shown at the miRNA 3′ end targeted regions, as well as the 5′ end (seed region) binding sites.     

Functional and Evolutionary Significance of Human MicroRNA Seed Region Mutations

MicroRNAs have emerged in recent years as important regulators of cell function in both normal and diseased cells. MiRNAs coordinately regulate large suites of target genes by mRNA degradation and/or translational inhibition. The mRNA target specificities of miRNAs in animals are primarily encoded within a 7 nt “seed region” mapping to positions 2–8 at the molecule''s 5′ end. We here combine computational analyses with experimental studies to explore the functional significance of sequence variation within the seed region of human miRNAs. The results indicate that a substitution of even a single nucleotide within the seed region changes the spectrum of mRNA targets by >50%. The high functional cost of even single nucleotide changes within seed regions is consistent with their high sequence conservation among miRNA families both within and between species and suggests processes that may underlie the evolution of miRNA regulatory control.     

Comprehensive thermodynamic analysis of 3' double-nucleotide overhangs neighboring Watson-Crick terminal base pairs

Thermodynamic parameters are reported for duplex formation of 48 self-complementary RNA duplexes containing Watson–Crick terminal base pairs (GC, AU and UA) with all 16 possible 3′ double-nucleotide overhangs; mimicking the structures of short interfering RNAs (siRNA) and microRNAs (miRNA). Based on nearest-neighbor analysis, the addition of a second dangling nucleotide to a single 3′ dangling nucleotide increases stability of duplex formation up to 0.8 kcal/mol in a sequence dependent manner. Results from this study in conjunction with data from a previous study [A. S. O'Toole, S. Miller and M. J. Serra (2005) RNA, 11, 512.] allows for the development of a refined nearest-neighbor model to predict the influence of 3′ double-nucleotide overhangs on the stability of duplex formation. The model improves the prediction of free energy and melting temperature when tested against five oligomers with various core duplex sequences. Phylogenetic analysis of naturally occurring miRNAs was performed to support our results. Selection of the effector miR strand of the mature miRNA duplex appears to be dependent upon the identity of the 3′ double-nucleotide overhang. Thermodynamic parameters for 3′ single terminal overhangs adjacent to a UA pair are also presented.     

Evolution after Whole-Genome Duplication: Teleost MicroRNAs

MicroRNAs (miRNAs) are important gene expression regulators implicated in many biological processes, but we lack a global understanding of how miRNA genes evolve and contribute to developmental canalization and phenotypic diversification. Whole-genome duplication events likely provide a substrate for species divergence and phenotypic change by increasing gene numbers and relaxing evolutionary pressures. To understand the consequences of genome duplication on miRNA evolution, we studied miRNA genes following the teleost genome duplication (TGD). Analysis of miRNA genes in four teleosts and in spotted gar, whose lineage diverged before the TGD, revealed that miRNA genes were retained in ohnologous pairs more frequently than protein-coding genes, and that gene losses occurred rapidly after the TGD. Genomic context influenced retention rates, with clustered miRNA genes retained more often than nonclustered miRNA genes and intergenic miRNA genes retained more frequently than intragenic miRNA genes, which often shared the evolutionary fate of their protein-coding host. Expression analyses revealed both conserved and divergent expression patterns across species in line with miRNA functions in phenotypic canalization and diversification, respectively. Finally, major strands of miRNA genes experienced stronger purifying selection, especially in their seeds and 3′-complementary regions, compared with minor strands, which nonetheless also displayed evolutionary features compatible with constrained function. This study provides the first genome-wide, multispecies analysis of the mechanisms influencing metazoan miRNA evolution after whole-genome duplication.     

The Fate of miRNA* Strand through Evolutionary Analysis: Implication for Degradation As Merely Carrier Strand or Potential Regulatory Molecule?

Background

During typical microRNA (miRNA) biogenesis, one strand of a ∼22 nt RNA duplex is preferentially selected for entry into a silencing complex, whereas the other strand, known as the passenger strand or miRNA* strand, is degraded. Recently, some miRNA* sequences were reported as guide miRNAs with abundant expression. Here, we intended to discover evolutionary implication of the fate of miRNA* strand by analyzing miRNA/miRNA* sequences across vertebrates.

Principal Findings

Mature miRNAs based on gene families were well conserved especially for their seed sequences across vertebrates, while their passenger strands always showed various divergence patterns. The divergence mainly resulted from divergence of different animal species, homologous miRNA genes and multicopy miRNA hairpin precursors. Some miRNA* sequences were phylogenetically conserved in seed and anchor sequences similar to mature miRNAs, while others revealed high levels of nucleotide divergence despite some of their partners were highly conserved. Most of those miRNA precursors that could generate abundant miRNAs from both strands always were well conserved in sequences of miR-#-5p and miR-#-3p, especially for their seed sequences.

Conclusions

The final fate of miRNA* strand, either degraded as merely carrier strand or expressed abundantly as potential functional guide miRNA, may be destined across evolution. Well-conserved miRNA* strands, particularly conservation in seed sequences, maybe afford potential opportunities for contributing to regulation network. The study will broaden our understanding of potential functional miRNA* species.     

HEN1 recognizes 21-24 nt small RNA duplexes and deposits a methyl group onto the 2' OH of the 3' terminal nucleotide

microRNAs (miRNAs) and small interfering RNAs (siRNAs) in plants bear a methyl group on the ribose of the 3′ terminal nucleotide. We showed previously that the methylation of miRNAs and siRNAs requires the protein HEN1 in vivo and that purified HEN1 protein methylates miRNA/miRNA* duplexes in vitro. In this study, we show that HEN1 methylates both miRNA/miRNA* and siRNA/siRNA* duplexes in vitro with a preference for 21–24 nt RNA duplexes with 2 nt overhangs. We also demonstrate that HEN1 deposits the methyl group on to the 2′ OH of the 3′ terminal nucleotide. Among various modifications that can occur on the ribose of the terminal nucleotide, such as 2′-deoxy, 3′-deoxy, 2′-O-methyl and 3′-O-methyl, only 2′-O-methyl on a small RNA inhibits the activity of yeast poly(A) polymerase (PAP). These findings indicate that HEN1 specifically methylates miRNAs and siRNAs and implicate the importance of the 2′-O-methyl group in the biology of RNA silencing.     

Repertoire of Bovine miRNA and miRNA-Like Small Regulatory RNAs Expressed upon Viral Infection

MicroRNA (miRNA) and other types of small regulatory RNAs play a crucial role in the regulation of gene expression in eukaryotes. Several distinct classes of small regulatory RNAs have been discovered in recent years. To extend the repertoire of small RNAs characterized in mammals and to examine relationship between host miRNA expression and viral infection we used Illumina''s ultrahigh throughput sequencing approach. We sequenced three small RNA libraries prepared from cell line derived from the adult bovine kidney under normal conditions and upon infection of the cell line with Bovine herpesvirus 1. We used a bioinformatics approach to distinguish authentic mature miRNA sequences from other classes of small RNAs and short RNA fragments represented in the sequencing data. Using this approach we detected 219 out of 356 known bovine miRNAs and 115 respective miRNA* sequences. In addition we identified five new bovine orthologs of known mammalian miRNAs and discovered 268 new cow miRNAs many of which are not identifiable in other mammalian genomes and thus might be specific to the ruminant lineage. In addition we found seven new bovine mirtron candidates. We also discovered 10 small nucleolar RNA (snoRNA) loci that give rise to small RNA with possible miRNA-like function. Results presented in this study extend our knowledge of the biology and evolution of small regulatory RNAs in mammals and illuminate mechanisms of small RNA biogenesis and function. New miRNA sequences and the original sequencing data have been submitted to miRNA repository (miRBase) and NCBI GEO archive respectively. We envisage that these resources will facilitate functional annotation of the bovine genome and promote further functional and comparative genomics studies of small regulatory RNA in mammals.     

Characterization of the mammalian miRNA turnover landscape

Steady state cellular microRNA (miRNA) levels represent the balance between miRNA biogenesis and turnover. The kinetics and sequence determinants of mammalian miRNA turnover during and after miRNA maturation are not fully understood. Through a large-scale study on mammalian miRNA turnover, we report the co-existence of multiple cellular miRNA pools with distinct turnover kinetics and biogenesis properties and reveal previously unrecognized sequence features for fast turnover miRNAs. We measured miRNA turnover rates in eight mammalian cell types with a combination of expression profiling and deep sequencing. While most miRNAs are stable, a subset of miRNAs, mostly miRNA*s, turnovers quickly, many of which display a two-step turnover kinetics. Moreover, different sequence isoforms of the same miRNA can possess vastly different turnover rates. Fast turnover miRNA isoforms are enriched for 5′ nucleotide bias against Argonaute-(AGO)-loading, but also additional 3′ and central sequence features. Modeling based on two fast turnover miRNA*s miR-222-5p and miR-125b-1-3p, we unexpectedly found that while both miRNA*s are associated with AGO, they strongly differ in HSP90 association and sensitivity to HSP90 inhibition. Our data characterize the landscape of genome-wide miRNA turnover in cultured mammalian cells and reveal differential HSP90 requirements for different miRNA*s. Our findings also implicate rules for designing stable small RNAs, such as siRNAs.