The impact of age, biogenesis, and genomic clustering on Drosophila microRNA evolution

Parameters of clustered suboptimal miRNA biogenesis

The nuclear cleavage of a suboptimal primary miRNA hairpin by the Drosha/DGCR8 complex ("Microprocessor") can be enhanced by an optimal miRNA neighbor, a phenomenon termed cluster assistance. Several features and biological impacts of this new layer of miRNA regulation are not fully known. Here, we elucidate the parameters of cluster assistance of a suboptimal miRNA and also reveal competitive interactions amongst optimal miRNAs within a cluster. We exploit cluster assistance as a functional assay for suboptimal processing and use this to invalidate putative suboptimal substrates, as well as identify a "solo" suboptimal miRNA. Finally, we report complexity in how specific mutations might affect the biogenesis of clustered miRNAs in disease contexts. This includes how an operon context can buffer the effect of a deleterious processing variant, but reciprocally how a point mutation can have a nonautonomous effect to impair the biogenesis of a clustered, suboptimal, neighbor. These data expand our knowledge regarding regulated miRNA biogenesis in humans and represent a functional assay for empirical definition of suboptimal Microprocessor substrates.

Accurate microRNA annotation of animal genomes using trained covariance models of curated microRNA complements in MirMachine

The annotation of microRNAs depends on the availability of transcriptomics data and expert knowledge. This has led to a gap between the availability of novel genomes and high-quality microRNA complements. Using >16,000 microRNAs from the manually curated microRNA gene database MirGeneDB, we generated trained covariance models for all conserved microRNA families. These models are available in our tool MirMachine, which annotates conserved microRNAs within genomes. We successfully applied MirMachine to a range of animal species, including those with large genomes and genome duplications and extinct species, where small RNA sequencing is hard to achieve. We further describe a microRNA score of expected microRNAs that can be used to assess the completeness of genome assemblies. MirMachine closes a long-persisting gap in the microRNA field by facilitating automated genome annotation pipelines and deeper studies into the evolution of genome regulation, even in extinct organisms.

Splice and Dice: Intronic microRNAs, Splicing and Cancer

Introns span only a quarter of the human genome, yet they host around 60% of all known microRNAs. Emerging evidence indicates the adaptive advantage of microRNAs residing within introns is attributed to their complex co-regulation with transcription and alternative splicing of their host genes. Intronic microRNAs are often co-expressed with their host genes, thereby providing functional synergism or antagonism that is exploited or decoupled in cancer. Additionally, intronic microRNA biogenesis and the alternative splicing of host transcript are co-regulated and intertwined. The importance of intronic microRNAs is under-recognized in relation to the pathogenesis of cancer.

Biogenesis, characterization, and functions of mirtrons

MicroRNAs (miRNAs) are major post-transcriptional regulators of gene expression. They base pair with the complementary target mRNA at the 3'UTR and modulate cellular processes by repressing the mRNA translation or degrading the mRNA. There are well-documented mechanisms of biogenesis of miRNA; however, a sizeable number of miRNAs are also produced by non-canonical pathways. Mirtrons represent a predominant class of non-canonical miRNAs. Mirtrons originate from intronic regions and are produced in a splicing-dependent and Drosha-independent manner. Mirtrons constitute about 15% of all miRNAs produced in a human body and have caught attention of researchers worldwide due to their unconventional origin, sequence characteristics, evolutionary dynamics, ability to regulate variety of cellular processes and their immense potential in disease therapeutics. In this comprehensive review we collate the research done in the past decade including biogenesis, sequence characteristics, regulation, and emerging therapeutic roles of mirtrons. This article is categorized under: RNA Processing > Processing of Small RNAs Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs Regulatory RNAs/RNAi/Riboswitches > RNAi: Mechanisms of Action.

Run or Die in the Evolution of New MicroRNAs-Testing the Red Queen Hypothesis on De Novo New Genes

The Red Queen hypothesis depicts evolution as the continual struggle to adapt. According to this hypothesis, new genes, especially those originating from nongenic sequences (i.e., de novo genes), are eliminated unless they evolve continually in adaptation to a changing environment. Here, we analyze two Drosophila de novo miRNAs that are expressed in a testis-specific manner with very high rates of evolution in their DNA sequence. We knocked out these miRNAs in two sibling species and investigated their contributions to different fitness components. We observed that the fitness contributions of miR-975 in Drosophila simulans seem positive, in contrast to its neutral contributions in D. melanogaster, whereas miR-983 appears to have negative contributions in both species, as the fitness of the knockout mutant increases. As predicted by the Red Queen hypothesis, the fitness difference of these de novo miRNAs indicates their different fates.

Elucidation of the microRNA Transcriptome in Western Corn Rootworm Reveals Its Dynamic and Evolutionary Complexity

Diabrotica virgifera virgifera (western corn rootworm, WCR) is one of the most destructive agricultural insect pests in North America. It is highly adaptive to environmental stimuli and crop protection technologies. However, little is known about the underlying genetic basis of WCR behavior and adaptation. More specifically, the involvement of small RNAs (sRNAs), especially microRNAs (miRNAs), a class of endogenous small non-coding RNAs that regulate various biological processes, has not been examined, and the datasets of putative sRNA sequences have not previously been generated for WCR. To achieve a comprehensive collection of sRNA transcriptomes in WCR, we constructed, sequenced, and analyzed sRNA libraries from different life stages of WCR and northern corn rootworm (NCR), and identified 101 conserved precursor miRNAs (pre-miRNAs) in WCR and other Arthropoda. We also identified 277 corn rootworm specific pre-miRNAs. Systematic analyses of sRNA populations in WCR revealed that its sRNA transcriptome, which includes PIWI-interacting RNAs (piRNAs) and miRNAs, undergoes a dynamic change throughout insect development. Phylogenetic analysis of miRNA datasets from model species reveals that a large pool of species-specific miRNAs exists in corn rootworm; these are potentially evolutionarily transient. Comparisons of WCR miRNA clusters to other insect species highlight conserved miRNA-regulated processes that are common to insects. Parallel Analysis of RNA Ends (PARE) also uncovered potential miRNA-guided cleavage sites in WCR. Overall, this study provides a new resource for studying the sRNA transcriptome and miRNA-mediated gene regulation in WCR and other Coleopteran insects.

Cloning and characterization of Bubaline mammary miRNAs: An in silico approach

MicroRNAs (miRNAs) are small non-coding RNAs ~ 19-25 nucleotides long that are involved in the regulation of gene expression. They negatively regulate the gene expression via inhibition or complete degradation of mRNAs by binding the complementary target sequences in 3' untranslated region. The present investigation was aimed at profiling of miRNAs expressed in the Bubaline mammary tissue at dry stage of lactation cycle. Small RNAs were isolated from freshly collected mammary tissues and T4 RNA ligase was used to ligate the enriched miRNAs with 3' and 5' linker sequences in two separate reactions. cDNA copies were synthesized from linkered small RNAs followed by the PCR amplification. The PCR products were resolved on 15% non-denaturing polyacrylamide gel electrophoresis by gelstar staining. The PCR products were cloned using pGEM®-T easy vector system and the desired clones (with linkered small RNA sequences) were confirmed using restriction digestion of plasmids with EcoRI. Out of 15 Bubaline small RNA sequences, eight sequences (Seq. ID I-VIII) matched the size range of miRNA molecules i.e., 18-26 nucleotides. The Bubaline small RNA sequences II and III showed partial alignment with various mammalian and non-mammalian miRNAs. The small RNA sequences obtained in the present study did not show any perfect match with already reported mRNA, rRNA or tRNA sequences in different databases. Hence, only the Bubaline small RNA sequences that showed partial homology with miRNAs were considered as putative Bubaline miRNAs. The present study established the basic repertoire of miRNAs expressed at dry stage of lactation in Bubaline mammary gland.

Regulation of Large Number of Weak Targets-New Insights from Twin-microRNAs

Each animal microRNA (miRNA) targets many genes for repression. Down-regulation of most of these targets is weak and has no detectable individual phenotypic effect. Whether this extensive weak repression is biologically relevant is a central issue in the debate on miRNA functionality. In the “small (target) pool” view, weak repression is nonfunctional and should be gradually removed during evolution. However, since the selective advantage of removing individual targets is small, testing this hypothesis is a challenge. We propose a novel approach by using miRNAs we call twin-miRs, which produce two mature products from the hairpin of the same miRNA precursor. Loss of the minor miR partner would affect all its targets and thus could be visible to selection. Since the minor miRs repress all their targets weakly, the “small pool” hypothesis would predict the elimination of twin-miRs over time. Surveying and sequencing 45 small RNA libraries in Drosophila, we found that nearly 40% of miRNAs produce twin-miRs. The minor forms are expressed in nontrivial abundance and repress their targets weakly. Interestingly, twin-miRs are often evolutionarily old, highly conserved, and comparable to solo-miRs in expression. Since there is no measurable trend toward reduction in target pool size, we conclude that at least some of the weak repression interactions are functional. A companion study using the May–Wigner theory of network stability suggests that distributed weak repression cumulatively contributes to stability of gene regulatory networks.

Predicting microRNA targeting efficacy in Drosophila

Background

MicroRNAs (miRNAs) are short regulatory RNAs that derive from hairpin precursors. Important for understanding the functional roles of miRNAs is the ability to predict the messenger RNA (mRNA) targets most responsive to each miRNA. Progress towards developing quantitative models of miRNA targeting in Drosophila and other invertebrate species has lagged behind that of mammals due to the paucity of datasets measuring the effects of miRNAs on mRNA levels.

Results

We acquired datasets suitable for the quantitative study of miRNA targeting in Drosophila. Analyses of these data expanded the types of regulatory sites known to be effective in flies, expanded the mRNA regions with detectable targeting to include 5′ untranslated regions, and identified features of site context that correlate with targeting efficacy in fly cells. Updated evolutionary analyses evaluated the probability of conserved targeting for each predicted site and indicated that more than a third of the Drosophila genes are preferentially conserved targets of miRNAs. Based on these results, a quantitative model was developed to predict targeting efficacy in insects. This model performed better than existing models, and it drives the most recent version, v7, of TargetScanFly.

Conclusions

Our evolutionary and functional analyses expand the known scope of miRNA targeting in flies and other insects. The existence of a quantitative model that has been developed and trained using Drosophila data will provide a valuable resource for placing miRNAs into gene regulatory networks of this important experimental organism.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1504-3) contains supplementary material, which is available to authorized users.

Importance of miRNA stability and alternative primary miRNA isoforms in gene regulation during Drosophila development

Mature microRNAs (miRNAs) are processed from primary transcripts (pri-miRNAs), and their expression is controlled at transcriptional and post-transcriptional levels. However, how regulation at multiple levels achieves precise control remains elusive. Using published and new datasets, we profile a time course of mature and pri-miRNAs in Drosophila embryos and reveal the dynamics of miRNA production and degradation as well as dynamic changes in pri-miRNA isoform selection. We found that 5' nucleotides influence stability of mature miRNAs. Furthermore, distinct half-lives of miRNAs from the mir-309 cluster shape their temporal expression patterns, and the importance of rapid degradation of the miRNAs in gene regulation is detected as distinct evolutionary signatures at the target sites in the transcriptome. Finally, we show that rapid degradation of miR-3/-309 may be important for regulation of the planar cell polarity pathway component Vang. Altogether, the results suggest that complex mechanisms regulate miRNA expression to support normal development.

MicroRNA duplication accelerates the recruitment of new targets during vertebrate evolution

The repertoire of miRNAs has considerably expanded during metazoan evolution, and duplication is an important mechanism for generating new functional miRNAs. However, relatively little is known about the functional divergence between paralogous miRNAs and the possible coevolution between duplicated miRNAs and the genomic contexts. By systematically examining small RNA expression profiles across various human tissues and interrogating the publicly available miRNA:mRNA pairing chimeras, we found that changes in expression patterns and targeting preferences are widespread for duplicated miRNAs in vertebrates. Both the empirical interactions and target predictions suggest that evolutionarily conserved homo-seed duplicated miRNAs pair with significantly higher numbers of target sites compared to the single-copy miRNAs. Our birth-and-death evolutionary analysis revealed that the new target sites of miRNAs experienced frequent gains and losses during function development. Our results suggest that a newly emerged target site has a higher probability to be functional and maintained by natural selection if it is paired to a seed shared by multiple paralogous miRNAs rather than being paired to a single-copy miRNA. We experimentally verified the divergence in target repression between two paralogous miRNAs by transfecting let-7a and let-7b mimics into kidney-derived cell lines of four mammalian species and measuring the resulting transcriptome alterations by extensive high-throughput sequencing. Our results also suggest that the gains and losses of let-7 target sites might be associated with the evolution of repressiveness of let-7 across mammalian species.

On the possibility of death of new genes - evidence from the deletion of de novo microRNAs

Background

New genes are constantly formed, sometimes from non-genic sequences, creating what is referred to as de novo genes. Since the total number of genes remains relatively steady, gene deaths likely balance out new births. In metazoan genomes, microRNAs (miRs) genes, small and non-coding, account for the bulk of functional de novo genes and are particularly suited to the investigation of gene death.

Results

In this study, we discover a Drosophila-specific de novo miRNA (mir-977) that may be facing impending death. Strikingly, after this testis-specific gene is deleted from D. melanogaster, most components of male fitness increase, rather than decrease as had been expected. These components include male viability, fertility and males’ ability to repress female re-mating. Given that mir-977 has a negative fitness effect in D. melanogaster, this de novo gene with an adaptive history for over 60 Myrs may be facing elimination. In some other species where mir-977 is not found, gene death may have already happened.

Conclusion

The surprising result suggests that de novo genes, constantly rising and falling during evolution, may often be transiently adaptive and then purged from the genome.

Electronic supplementary material

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

Deep experimental profiling of microRNA diversity, deployment, and evolution across the Drosophila genus

To assess miRNA evolution across the Drosophila genus, we analyzed several billion small RNA reads across 12 fruit fly species. These data permit comprehensive curation of species- and clade-specific variation in miRNA identity, abundance, and processing. Among well-conserved miRNAs, we observed unexpected cases of clade-specific variation in 5' end precision, occasional antisense loci, and putatively noncanonical loci. We also used strict criteria to identify a large set (649) of novel, evolutionarily restricted miRNAs. Within the bulk collection of species-restricted miRNAs, two notable subpopulations are splicing-derived mirtrons and testes-restricted, recently evolved, clustered (TRC) canonical miRNAs. We quantified miRNA birth and death using our annotation and a phylogenetic model for estimating rates of miRNA turnover. We observed striking differences in birth and death rates across miRNA classes defined by biogenesis pathway, genomic clustering, and tissue restriction, and even identified flux heterogeneity among Drosophila clades. In particular, distinct molecular rationales underlie the distinct evolutionary behavior of different miRNA classes. Mirtrons are associated with high rates of 3' untemplated addition, a mechanism that impedes their biogenesis, whereas TRC miRNAs appear to evolve under positive selection. Altogether, these data reveal miRNA diversity among Drosophila species and principles underlying their emergence and evolution.

Comparative genomic analysis of upstream miRNA regulatory motifs in Caenorhabditis

MicroRNAs (miRNAs) comprise a class of short noncoding RNA molecules that play diverse developmental and physiological roles by controlling mRNA abundance and protein output of the vast majority of transcripts. Despite the importance of miRNAs in regulating gene function, we still lack a complete understanding of how miRNAs themselves are transcriptionally regulated. To fill this gap, we predicted regulatory sequences by searching for abundant short motifs located upstream of miRNAs in eight species of Caenorhabditis nematodes. We identified three conserved motifs across the Caenorhabditis phylogeny that show clear signatures of purifying selection from comparative genomics, patterns of nucleotide changes in motifs of orthologous miRNAs, and correlation between motif incidence and miRNA expression. We then validated our predictions with transgenic green fluorescent protein reporters and site-directed mutagenesis for a subset of motifs located in an enhancer region upstream of let-7 We demonstrate that a CT-dinucleotide motif is sufficient for proper expression of GFP in the seam cells of adult C. elegans, and that two other motifs play incremental roles in combination with the CT-rich motif. Thus, functional tests of sequence motifs identified through analysis of molecular evolutionary signatures provide a powerful path for efficiently characterizing the transcriptional regulation of miRNA genes.

The Drosophila Dicer-1 Partner Loquacious Enhances miRNA Processing from Hairpins with Unstable Structures at the Dicing Site

In Drosophila, Dicer-1 binds Loquacious-PB (Loqs-PB) as its major co-factor. Previous analyses indicated that loqs mutants only partially impede miRNA processing, but the activity of minor isoforms or maternally deposited Loqs was not eliminated in these studies. We addressed this by generating a cell line from loqs-null embryos and found that only ∼40% of miRNAs showed clear Loqs dependence. Genome-wide comparison of the hairpin structure and Loqs dependence suggested that Loqs substrates are influenced by base-pairing status at the dicing site. Artificial alteration of base-pairing stability at this position in model miRNA hairpins resulted in predicted changes in Loqs dependence, providing evidence for this hypothesis. Finally, we found that evolutionarily young miRNA genes tended to be Loqs dependent. We propose that Loqs may have roles in assisting the de novo emergence of miRNA genes by facilitating dicing of suboptimal hairpin substrates.

Splicing-dependent expression of microRNAs of mirtron origin in human digestive and excretory system cancer cells

Background

An abundant class of intronic microRNAs (miRNAs) undergoes atypical Drosha-independent biogenesis in which the spliceosome governs the excision of hairpin miRNA precursors, called mirtrons. Although nearly 500 splicing-dependent miRNA candidates have been recently predicted via bioinformatic analysis of human RNA-Seq datasets, only a few of them have been experimentally validated. The detailed mechanism of miRNA processing by the splicing machinery and the roles of mirtronic miRNAs in cancer are yet to be uncovered.

Methods

We experimentally examined whether biogenesis of certain miRNAs is under a splicing control by analyzing their expression levels in response to alterations in the 5′- and 3′-splice sites of a series of intron-containing minigenes carrying appropriate miRNAs. The expression levels of the miRNAs processed from mirtrons were determined by quantitative real-time PCR in five digestive tract (pancreas PANC-1, SU.86.86, T3M4, stomach KATOIII, colon HCT116) and two excretory system (kidney CaKi-1, 786-O) carcinoma cell lines as well as in pancreatic, stomach, and colorectal tumors. Transiently expressed SRSF1 and SRSF2 splicing factors were quantified by western blotting in the nuclear fractions of HCT116 cells.

Results

We found that biogenesis of the human hsa-miR-1227-3p, hsa-miR-1229-3p, and hsa-miR-1236-3p is splicing-dependent; therefore, these miRNAs can be assigned to the class of miRNAs processed by a non-canonical mirtron pathway. The expression analysis revealed a differential regulation of human mirtronic miRNAs in various cancer cell lines and tumors. In particular, hsa-miR-1229-3p is selectively upregulated in the pancreatic and stomach cancer cell lines derived from metastatic sites. Compared with the healthy controls, the expression of hsa-miR-1226-3p was significantly higher in stomach tumors but extensively downregulated in colorectal tumors. Furthermore, we provided evidence that overexpression of SRSF1 or SRSF2 can upregulate the processing of individual mirtronic miRNAs in HCT116 cells.

Conclusions

An interplay of different splicing factors, such as SRSF1 or SRSF2, may alter the levels of miRNAs of mirtron origin in a cell. Our findings underline the specific expression profiles of mirtronic miRNAs in colorectal, stomach, and pancreatic cancer.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-016-0200-y) contains supplementary material, which is available to authorized users.

The significance and scope of evolutionary developmental biology: a vision for the 21st century

Evolutionary developmental biology (evo-devo) has undergone dramatic transformations since its emergence as a distinct discipline. This paper aims to highlight the scope, power, and future promise of evo-devo to transform and unify diverse aspects of biology. We articulate key questions at the core of eleven biological disciplines-from Evolution, Development, Paleontology, and Neurobiology to Cellular and Molecular Biology, Quantitative Genetics, Human Diseases, Ecology, Agriculture and Science Education, and lastly, Evolutionary Developmental Biology itself-and discuss why evo-devo is uniquely situated to substantially improve our ability to find meaningful answers to these fundamental questions. We posit that the tools, concepts, and ways of thinking developed by evo-devo have profound potential to advance, integrate, and unify biological sciences as well as inform policy decisions and illuminate science education. We look to the next generation of evolutionary developmental biologists to help shape this process as we confront the scientific challenges of the 21st century.

Selective Suppression of the Splicing-Mediated MicroRNA Pathway by the Terminal Uridyltransferase Tailor

Several terminal uridyltransferases (TUTases) are known to modulate small RNA biogenesis and/or function via diverse mechanisms. Here, we demonstrate that Drosophila splicing-derived pre-miRNAs (mirtrons) are efficiently modified by the previously uncharacterized TUTase, Tailor. Tailor is necessary and sufficient for mirtron hairpin uridylation, and this modification inhibits mirtron biogenesis. Genome-wide analyses demonstrate that mirtrons are dominant Tailor substrates, and three features contribute to substrate specificity. First, reprogramming experiments show Tailor preferentially identifies splicing-derived miRNAs. Second, in vitro tests indicate Tailor prefers substrate hairpins over mature miRNAs. Third, Tailor exhibits sequence preference for 3'-terminal AG, a defining mirtron characteristic. Our work supports the notion that Tailor preferentially suppresses biogenesis of mirtrons, an evolutionarily adventitious pre-miRNA substrate class. Moreover, we detect preferential activity of Tailor on 3'-G canonical pre-miRNAs, and specific depletion of such loci from the pool of conserved miRNAs. Thus, Tailor activity may have had collateral impact on shaping populations of canonical miRNAs.

Evolutionary Comparison of Two Combinatorial Regulators of SBP-Box Genes, MiR156 and MiR529, in Plants

A complete picture of the evolution of miRNA combinatorial regulation requires the synthesis of information on all miRNAs and their targets. MiR156 and miR529 are two combinatorial regulators of squamosa promoter binding protein-like (SBP-box) genes. Previous studies have clarified the evolutionary dynamics of their targets; however, there have been no reports on the evolutionary patterns of two miRNA regulators themselves to date. In this study, we investigated the evolutionary differences between these two miRNA families in extant land plants. Our work found that miR529 precursor, especially of its mature miRNA sequence, has a higher evolutionary rate. Such accelerating evolution of miR529 has significantly effects on its structural stability, and sequence conservation against existence of itself. By contrast, miR156 evolves more rapidly in loop region of the stable secondary structure, which may contribute to its functional diversity. Moreover, miR156 and miR529 genes have distinct rates of loss after identical duplication events. MiR529 genes have a higher average loss rate and asymmetric loss rate in duplicated gene pairs, indicating preferred miR529 gene losses become another predominant mode of inactivation, that are implicated in the contraction of this family. On the contrary, duplicated miR156 genes have a low loss rate, and could serve as another new source for functional diversity. Taken together, these results provide better insight into understanding the evolutionary divergence of miR156 and miR529 family in miRNA combinational regulation network.

Diverse modes of evolutionary emergence and flux of conserved microRNA clusters

Many animal miRNA loci reside in genomic clusters that generate multicistronic primary-miRNA transcripts. While clusters that contain copies of the same miRNA hairpin are clearly products of local duplications, the evolutionary provenance of clusters with disparate members is less clear. Recently, it was proposed that essentially all such clusters in Drosophila derived from de novo formation of miRNA-like hairpins within existing miRNA transcripts, and that the maintenance of multiple miRNAs in such clusters was due to evolutionary hitchhiking on a major cluster member. However, this model seems at odds with the fact that many such miRNA clusters are composed of well-conserved miRNAs. In an effort to trace the birth and expansion of miRNA clusters that are presently well-conserved across Drosophilids, we analyzed a broad swath of metazoan species, with particular emphasis on arthropod evolution. Beyond duplication and de novo birth, we highlight a diversity of modes that contribute to miRNA evolution, including neofunctionalization of miRNA copies, fissioning of locally duplicated miRNA clusters, miRNA deletion, and miRNA cluster expansion via the acquisition and/or neofunctionalization of miRNA copies from elsewhere in the genome. In particular, we suggest that miRNA clustering by acquisition represents an expedient strategy to bring cohorts of target genes under coordinate control by miRNAs that had already been individually selected for regulatory impact on the transcriptome.

Microevolution of nematode miRNAs reveals diverse modes of selection

Micro-RNA (miRNA) genes encode abundant small regulatory RNAs that play key roles during development and in homeostasis by fine tuning and buffering gene expression. This layer of regulatory control over transcriptional networks is preserved by selection across deep evolutionary time, yet selection pressures on individual miRNA genes in contemporary populations remain poorly characterized in any organism. Here, we quantify nucleotide variability for 129 miRNAs in the genome of the nematode Caenorhabditis remanei to understand the microevolution of this important class of regulatory genes. Our analysis of three population samples and C. remanei's sister species revealed ongoing natural selection that constrains evolution of all sequence domains within miRNA hairpins. We also show that new miRNAs evolve faster than older miRNAs but that selection nevertheless favors their persistence. Despite the ongoing importance of purging of new mutations, we discover a trove of >400 natural miRNA sequence variants that include single nucleotide polymorphisms in seed motifs, indels that ablate miRNA functional domains, and origination of new miRNAs by duplication. Moreover, we demonstrate substantial nucleotide divergence of pre-miRNA hairpin alleles between populations and sister species. These findings from the first global survey of miRNA microevolution in Caenorhabditis support the idea that changes in gene expression, mediated through divergence in miRNA regulation, can contribute to phenotypic novelty and adaptation to specific environments in the present day as well as the distant past.

Adaptive evolution of testis-specific, recently evolved, clustered miRNAs in Drosophila

The propensity of animal miRNAs to regulate targets bearing modest complementarity, most notably via pairing with miRNA positions ∼2-8 (the "seed"), is believed to drive major aspects of miRNA evolution. First, minimal targeting requirements have allowed most conserved miRNAs to acquire large target cohorts, thus imposing strong selection on miRNAs to maintain their seed sequences. Second, the modest pairing needed for repression suggests that evolutionarily nascent miRNAs may generally induce net detrimental, rather than beneficial, regulatory effects. Hence, levels and activities of newly emerged miRNAs are expected to be limited to preserve the status quo of gene expression. In this study, we unexpectedly show that Drosophila testes specifically express a substantial miRNA population that contravenes these tenets. We find that multiple genomic clusters of testis-restricted miRNAs harbor recently evolved miRNAs, whose experimentally verified orthologs exhibit divergent sequences, even within seed regions. Moreover, this class of miRNAs exhibits higher expression and greater phenotypic capacities in transgenic misexpression assays than do non-testis-restricted miRNAs of similar evolutionary age. These observations suggest that these testis-restricted miRNAs may be evolving adaptively, and several methods of evolutionary analysis provide strong support for this notion. Consistent with this, proof-of-principle tests show that orthologous miRNAs with divergent seeds can distinguish target sensors in a species-cognate manner. Finally, we observe that testis-restricted miRNA clusters exhibit extraordinary dynamics of miRNA gene flux in other Drosophila species. Altogether, our findings reveal a surprising tissue-directed influence of miRNA evolution, involving a distinct mode of miRNA function connected to adaptive gene regulation in the testis.

Identification and consequences of miRNA-target interactions--beyond repression of gene expression

Comparative genomics analyses and high-throughput experimental studies indicate that a microRNA (miRNA) binds to hundreds of sites across the transcriptome. Although the knockout of components of the miRNA biogenesis pathway has profound phenotypic consequences, most predicted miRNA targets undergo small changes at the mRNA and protein levels when the expression of the miRNA is perturbed. Alternatively, miRNAs can establish thresholds in and increase the coherence of the expression of their target genes, as well as reduce the cell-to-cell variability in target gene expression. Here, we review the recent progress in identifying miRNA targets and the emerging paradigms of how miRNAs shape the dynamics of target gene expression.

Sex-biased expression of microRNAs in Drosophila melanogaster

Most animals have separate sexes. The differential expression of gene products, in particular that of gene regulators, is underlying sexual dimorphism. Analyses of sex-biased expression have focused mostly on protein-coding genes. Several lines of evidence indicate that microRNAs, a class of major gene regulators, are likely to have a significant role in sexual dimorphism. This role has not been systematically explored so far. Here, I study the sex-biased expression pattern of microRNAs in the model species Drosophila melanogaster. As with protein-coding genes, sex-biased microRNAs are associated with the reproductive function. Strikingly, contrary to protein-coding genes, male-biased microRNAs are enriched in the X chromosome, whereas female microRNAs are mostly autosomal. I propose that the chromosomal distribution is a consequence of high rates of de novo emergence, and a preference for new microRNAs to be expressed in the testis. I also suggest that demasculinization of the X chromosome may not affect microRNAs. Interestingly, female-biased microRNAs are often encoded within protein-coding genes that are also expressed in females. MicroRNAs with sex-biased expression do not preferentially target sex-biased gene transcripts. These results strongly suggest that the sex-biased expression of microRNAs is mainly a consequence of high rates of microRNA emergence in the X chromosome (male bias) or hitchhiked expression by host genes (female bias).

Human-specific microRNA regulation of FOXO1: implications for microRNA recognition element evolution

MicroRNAs (miRNAs) have been established as important negative post-transcriptional regulators for gene expression. Within the past decade, miRNAs targeting transcription factors (TFs) has emerged as an important mechanism for gene expression regulation. Here, we tested the hypothesis that in TF 3'UTRs, human-specific single nucleotide change(s) that create novel miRNA recognition elements (MREs) contribute to species-specific differences in TF expression. From several potential human-specific TF MREs, one candidate, a member of the Forkhead Box O (FOXO) subclass in the Forkhead family known as Forkhead Box O1 (FOXO1; FKHR; NM_002015) was tested further. Human FOXO1 contains two sites predicted to confer miR-183-mediated post-transcriptional regulation: one specific to humans and the other conserved. Utilizing dual luciferase expression reporters, we show that only the human FOXO1 3'UTR contains a functional miR-183 site, not found in chimpanzee or mouse 3'untranslated regions (UTRs). Site-directed mutagenesis supports functionality of the human-specific miR-183 site, but not the conserved miR-183 site. Via overexpression and target site protection assays, we show that human FOXO1 is regulated by miR-183, but mouse FOXO1 is not. Finally, FOXO1-regulated cellular phenotypes, including cell invasion and proliferation, are impacted by miR-183 targeting only in human cells. These results provide strong evidence for human-specific gain of TF MREs, a process that may underlie evolutionary differences between phylogenic groups.