The evolution of plant microRNAs: insights from a basal eudicot sacred lotus

microRNAs (miRNAs) are important noncoding small RNAs that regulate mRNAs in eukaryotes. However, under which circumstances different miRNAs/miRNA families exhibit different evolutionary trajectories in plants remains unclear. In this study, we sequenced the small RNAs and degradome from a basal eudicot, sacred lotus (Nelumbo nucifera or lotus), to identify miRNAs and their targets. Combining with public miRNAs, we predicted 57 pre‐eudicot miRNA families from different evolutionary stages. We found that miRNA families featuring older age, higher copy and target number tend to show lower propensity for miRNA family loss (PGL) and stronger signature of purifying selection during divergence of temperate and tropical lotus. Further analyses of lotus genome revealed that there is an association between loss of miRNA families in descendent plants and in duplicated genomes. Gene dosage balance is crucial in maintaining those preferentially retained MIRNA duplicates by imposing stronger purifying selection. However, these factors and selection influencing miRNA family evolution are not applicable to the putative MIRNA‐likes. Additionally, the MIRNAs participating in lotus pollen–pistil interaction, a conserved process in angiosperms, also have a strong signature of purifying selection. Functionally, sequence divergence in MIRNAs escalates expression divergence of their target genes between temperate and tropical lotus during rhizome and leaf growth. Overall, our study unravels several important factors and selection that determine the miRNA family distribution in plants and duplicated genomes, and provides evidence for functional impact of MIRNA sequence evolution.     

MicroRNA genes derived from repetitive elements and expanded by segmental duplication events in mammalian genomes

MicroRNAs (miRNAs) are a class of small noncoding RNAs that regulate gene expression by targeting mRNAs for translation repression or mRNA degradation. Many miRNAs are being discovered and studied, but in most cases their origin, evolution and function remain unclear. Here, we characterized miRNAs derived from repetitive elements and miRNA families expanded by segmental duplication events in the human, rhesus and mouse genomes. We applied a comparative genomics approach combined with identifying miRNA paralogs in segmental duplication pair data in a genome-wide study to identify new homologs of human miRNAs in the rhesus and mouse genomes. Interestingly, using segmental duplication pair data, we provided credible computational evidence that two miRNA genes are located in the pseudoautosomal region of the human Y chromosome. We characterized all the miRNAs whether they were derived from repetitive elements or not and identified significant differences between the repeat-related miRNAs (RrmiRs) and non-repeat-derived miRNAs in (1) their location in protein-coding and intergenic regions in genomes, (2) the minimum free energy of their hairpin structures, and (3) their conservation in vertebrate genomes. We found some lineage-specific RrmiR families and three lineage-specific expansion families, and provided evidence indicating that some RrmiR families formed and expanded during evolutionary segmental duplication events. We also provided computational and experimental evidence for the functions of the conservative RrmiR families in the three species. Together, our results indicate that repetitive elements contribute to the origin of miRNAs, and large segmental duplication events could prompt the expansion of some miRNA families, including RrmiR families. Our study is a valuable contribution to the knowledge of evolution and function of non-coding region in genome.     

Computational identification of microRNA gene loci and precursor microRNA sequences in CHO cell lines

MicroRNAs (miRNAs) have recently entered Chinese hamster ovary (CHO) cell culture technology, due to their severe impact on the regulation of cellular phenotypes. Applications of miRNAs that are envisioned range from biomarkers of favorable phenotypes to cell engineering targets. These applications, however, require a profound knowledge of miRNA sequences and their genomic organization, which exceeds the currently available information of ~400 conserved mature CHO miRNA sequences. Based on these recently published sequences and two independent CHO-K1 genome assemblies, this publication describes the computational identification of CHO miRNA genomic loci. Using BLAST alignment, 415 previously reported CHO miRNAs were mapped to the reference genomes, and subsequently assigned to a distinct genomic miRNA locus. Sequences of the respective precursor-miRNAs were extracted from both reference genomes, folded in silico to verify correct structures and cross-compared. In the end, 212 genomic loci and pre-miRNA sequences representing 319 expressed mature miRNAs (approximately 50% of miRNAs represented matching pairs of 5' and 3' miRNAs) were submitted to the miRBase miRNA repository. As a proof-of-principle for the usability of the published genomic loci, four likely polycistronic miRNA cluster were chosen for PCR amplification using CHO-K1 and DHFR (-) genomic DNA. Overall, these data on the genomic context of miRNA expression in CHO will simplify the development of tools employing stable overexpression or deletion of miRNAs, allow the identification of miRNA promoters and improve detection methods such as microarrays.     

Sequence variation of MicroRNAs and their binding sites in Arabidopsis

Major differences exist between plants and animals both in the extent of microRNA (miRNA)-based gene regulation and the sequence complementarity requirements for miRNA-messenger RNA pairing. Whether these differences affect how these sites evolve at the molecular level is unknown. To determine the extent of sequence variation at miRNAs and their targets in a plant species, we resequenced 16 miRNA families (66 miRNAs in total) and all 52 of the characterized binding sites for these miRNAs in the plant model Arabidopsis (Arabidopsis thaliana), accounting for around 50% of the known miRNAs and binding sites in this species. As has been shown previously in humans, we find that both miRNAs and their target binding sites have very low nucleotide variation and divergence compared to their flanking sequences in Arabidopsis, indicating strong purifying selection on these sites in this species. Sequence data flanking the mature miRNAs, however, exhibit normal levels of polymorphism for the accessions in this study and, in some cases, nonneutral evolution or subtle effects on predicted pre-miRNA secondary structure, suggesting that there is raw material for the differential function of miRNA alleles. Overall, our results show that despite differences in the architecture of miRNA-based regulation, miRNAs and their targets are similarly constrained in both plants and animals.     

Conservation and divergence in plant microRNAs

MicroRNAs (miRNAs) are a class of small, non-coding RNAs that regulate gene expression in eukaryotic cells. The past decade has seen an explosion in our understanding of the sets of miRNA genes encoded in the genomes in different species of plants and the mechanisms by which miRNAs interact with target RNAs. A subset of miRNA families (and their binding sites in target RNAs) are conserved between angiosperms and basal plants, suggesting they predate the divergence of existing lineages of plants. However, the majority of miRNA families expressed by any given plant species have a narrow phylogenetic distribution. As a group, these "young" miRNAs genes appear to be evolutionarily fluid and lack clearly understood biological function. The goal of this review is to summarize our understanding of the sets of miRNA genes and miRNA targets that exist in various plant species and to discuss hypotheses that explain the patterns of conservation and divergence observed among microRNAs in plants.     

Identification of microRNAs of the herpesvirus family

Epstein-Barr virus (EBV or HHV4), a member of the human herpesvirus (HHV) family, has recently been shown to encode microRNAs (miRNAs). In contrast to most eukaryotic miRNAs, these viral miRNAs do not have close homologs in other viral genomes or in the genome of the human host. To identify other miRNA genes in pathogenic viruses, we combined a new miRNA gene prediction method with small-RNA cloning from several virus-infected cell types. We cloned ten miRNAs in the Kaposi sarcoma-associated virus (KSHV or HHV8), nine miRNAs in the mouse gammaherpesvirus 68 (MHV68) and nine miRNAs in the human cytomegalovirus (HCMV or HHV5). These miRNA genes are expressed individually or in clusters from either polymerase (pol) II or pol III promoters, and share no substantial sequence homology with one another or with the known human miRNAs. Generally, we predicted miRNAs in several large DNA viruses, and we could neither predict nor experimentally identify miRNAs in the genomes of small RNA viruses or retroviruses.     

Young intragenic miRNAs are less coexpressed with host genes than old ones: implications of miRNA-host gene coevolution

MicroRNAs (miRNAs) have emerged as key regulators of gene expression. Intragenic miRNAs account for ～50% of mammalian miRNAs. Classic studies reported that they are usually coexpressed with host genes. Here, using genome-wide miRNA and gene expression profiles from five sample sets, we show that evolutionarily conserved ('old') intragenic miRNAs tend to be coexpressed with host genes, but non-conserved ('young') ones rarely do so. This result is robust: in all sample sets, the coexpression rate of young miRNAs is significantly lower than that of conserved ones even after controlling for abundance. As a result, although young miRNAs dominate in human genome, the majority of intragenic miRNAs that show coexpression with host genes are phylogenetically old ones. For younger miRNAs, extrapolation of their expression profiles from those of their host genes should be treated with caution. We propose a model to explain this phenomenon in which the majority of young miRNAs are unlikely to be coexpressed with host genes; however, for some fraction of young miRNAs coexpression with their host genes, initially imbued by chromatin level effects, is advantageous and these are the ones likely to embed into the system and evolve ever higher levels of coexpression, possibly by evolving piggybacking mechanisms.     

microRNA complements in deuterostomes: origin and evolution of microRNAs

SUMMARY Although numerous studies have emphasized the role of microRNAs (miRNAs) in the control of many different cellular processes, they might also exert a profound effect on the macroevolution of animal body plans. It has been hypothesized that, because miRNAs increase genic precision and are continuously being added to metazoan genomes through geologic time, miRNAs might be instrumental for canalization of development and morphological evolution. Nonetheless, an outstanding question remains: how are new miRNAs constantly evolving? To address this question, we assessed the miRNA complements of four deuterostome species, chosen because of their sequenced genomes and well‐resolved phylogeny. Our comparative analysis shows that each of these four species is characterized by a unique repertoire of miRNAs, with few instances of miRNA loss. Moreover, we find that almost half of the miRNAs identified in this study are located in intronic regions of protein coding genes, suggesting that new miRNAs might arise from intronic regions in a process we term intronic exaptation. We also show that miRNAs often occur within cotranscribed clusters, and describe the biological function of one of these conserved clusters, the miR‐1/miR‐133 cluster. Taken together, our work shows that miRNAs can easily emerge within already transcribed regions of DNA, whether it be introns or preexisting clusters of miRNAs and/or miRNAs and protein coding genes, and because of their regulatory roles, these novel players change the structure of gene regulatory networks, with potential macroevolutionary results.     

Sequence Fingerprints of MicroRNA Conservation

It is known that the conservation of protein-coding genes is associated with their sequences both various species, such as animals and plants. However, the association between microRNA (miRNA) conservation and their sequences in various species remains unexplored. Here we report the association of miRNA conservation with its sequence features, such as base content and cleavage sites, suggesting that miRNA sequences contain the fingerprints for miRNA conservation. More interestingly, different species show different and even opposite patterns between miRNA conservation and sequence features. For example, mammalian miRNAs show a positive/negative correlation between conservation and AU/GC content, whereas plant miRNAs show a negative/positive correlation between conservation and AU/GC content. Further analysis puts forward the hypothesis that the introns of protein-coding genes may be a main driving force for the origin and evolution of mammalian miRNAs. At the 5′ end, conserved miRNAs have a preference for base U, while less-conserved miRNAs have a preference for a non-U base in mammals. This difference does not exist in insects and plants, in which both conserved miRNAs and less-conserved miRNAs have a preference for base U at the 5′ end. We further revealed that the non-U preference at the 5′ end of less-conserved mammalian miRNAs is associated with miRNA function diversity, which may have evolved from the pressure of a highly sophisticated environmental stimulus the mammals encountered during evolution. These results indicated that miRNA sequences contain the fingerprints for conservation, and these fingerprints vary according to species. More importantly, the results suggest that although species share common mechanisms by which miRNAs originate and evolve, mammals may develop a novel mechanism for miRNA origin and evolution. In addition, the fingerprint found in this study can be predictor of miRNA conservation, and the findings are helpful in achieving a clearer understanding of miRNA function and evolution.     

Global analysis of microRNA target gene expression reveals that miRNA targets are lower expressed in mature mouse and Drosophila tissues than in the embryos

MicroRNAs (miRNAs) are non-coding small RNAs of ~22 nt that regulate the gene expression by base pairing with target mRNAs, leading to mRNA cleavage or translational repression. It is currently estimated that miRNAs account for ~1% of predicted genes in higher eukaryotic genomes and that up to 30% of genes might be regulated by miRNAs. However, only very few miRNAs have been functionally characterized and the general functions of miRNAs are not globally studied. In this study, we systematically analyzed the expression patterns of miRNA targets using several public microarray profiles. We found that the expression levels of miRNA targets are lower in all mouse and Drosophila tissues than in the embryos. We also found miRNAs more preferentially target ubiquitously expressed genes than tissue-specifically expressed genes. These results support the current suggestion that miRNAs are likely to be largely involved in embryo development and maintaining of tissue identity.     

Accuracy of MicroRNA Discovery Pipelines in Non-Model Organisms Using Closely Related Species Genomes

Mapping small reads to genome reference is an essential and more common approach to identify microRNAs (miRNAs) in an organism. Using closely related species genomes as proxy references can facilitate miRNA expression studies in non-model species that their genomes are not available. However, the level of error this introduces is mostly unknown, as this is the result of evolutionary distance between the proxy reference and the species of interest. To evaluate the accuracy of miRNA discovery pipelines in non-model organisms, small RNA library data from a mosquito, Aedes aegypti, were mapped to three well annotated insect genomes as proxy references using miRanalyzer with two strict and loose mapping criteria. In addition, another web-based miRNA discovery pipeline (DSAP) was used as a control for program performance. Using miRanalyzer, more than 80% reduction was observed in the number of mapped reads using strict criterion when proxy genome references were used; however, only 20% reduction was recorded for mapped reads to other species known mature miRNA datasets. Except a few changes in ranking, mapping criteria did not make any significant differences in the profile of the most abundant miRNAs in A. aegypti when its original or a proxy genome was used as reference. However, more variation was observed in miRNA ranking profile when DSAP was used as analysing tool. Overall, the results also suggested that using a proxy reference did not change the most abundant miRNAs’ differential expression profiles when infected or non-infected libraries were compared. However, usage of a proxy reference could provide about 67% of the original outcome from more extremely up- or down-regulated miRNA profiles. Although using closely related species genome incurred some losses in the number of miRNAs, the most abundant miRNAs along with their differential expression profile would be acceptable based on the sensitivity level of each project.     

Identification of tissue-specific microRNAs from mouse

MicroRNAs (miRNAs) are a new class of noncoding RNAs, which are encoded as short inverted repeats in the genomes of invertebrates and vertebrates. It is believed that miRNAs are modulators of target mRNA translation and stability, although most target mRNAs remain to be identified. Here we describe the identification of 34 novel miRNAs by tissue-specific cloning of approximately 21-nucleotide RNAs from mouse. Almost all identified miRNAs are conserved in the human genome and are also frequently found in nonmammalian vertebrate genomes, such as pufferfish. In heart, liver, or brain, it is found that a single, tissue-specifically expressed miRNA dominates the population of expressed miRNAs and suggests a role for these miRNAs in tissue specification or cell lineage decisions. Finally, a miRNA was identified that appears to be the fruitfly and mammalian ortholog of C. elegans lin-4 stRNA.     

Mechanisms of microRNA-mediated gene regulation

microRNAs (miRNAs) are identified as a class of non-protein regulators and a new source for broad control of gene expression in eukaryotes. The past years have witnessed substantial progress in understanding miRNA functions and mechanisms, although a few controversies remain. Various hypotheses and models have been suggested for the mechanisms of miRNA repression, including translational inhibition at the level of initiation or elongation, rapid degradation of the nascent peptide, mRNA degradation, and mRNA sequestration into P bodies (processing bodies) and SGs (stress granules) for degradation or/and storage. Recently, some noncanonical miRNA regulation, such as miRNA activation and de-repression of miRNA inhibition, have been uncovered. This review discusses some recent advances about how miRNAs regulate their targets and various modes of miRNA function.