The evolution of evolvability in microRNA target sites in vertebrates

Regulation of soldier caste differentiation by microRNAs in Formosan subterranean termite (Coptotermes formosanus Shiraki)

The soldier caste is one of the most distinguished castes inside the termite colony. The mechanism of soldier caste differentiation has mainly been studied at the transcriptional level, but the function of microRNAs (miRNAs) in soldier caste differentiation is seldom studied. In this study, the workers of Coptotermes formosanus Shiraki were treated with methoprene, a juvenile hormone analog which can induce workers to transform into soldiers. The miRNomes of the methoprene-treated workers and the controls were sequenced. Then, the differentially expressed miRNAs (DEmiRs) were corrected with the differentially expressed genes DEGs to construct the DEmiR-DEG regulatory network. Afterwards, the DEmiR-regulated DEGs were subjected to GO enrichment and KEGG enrichment analysis. A total of 1,324 miRNAs were identified, among which 116 miRNAs were screened as DEmiRs between the methoprene-treated group and the control group. A total of 4,433 DEmiR-DEG pairs were obtained. No GO term was recognized as significant in the cellular component, molecular function, or biological process categories. The KEGG enrichment analysis of the DEmiR-regulated DEGs showed that the ribosome biogenesis in eukaryotes and circadian rhythm-fly pathways were enriched. This study demonstrates that DEmiRs and DEGs form a complex network regulating soldier caste differentiation in termites.

Diploid and triploid Chinook salmon (Oncorhynchus tshawytscha) have altered microRNA responses in immune tissues after infection with Vibrio anguillarum

Production of sterile fishes through artificial retention of a third set of chromosomes (triploidy) is a sustainable alternative for aquaculture since it reduces escapee pressure on wild populations. However, these fishes have reduced survival in stressful conditions and in response to infection. In this study, the impact of Vibrio anguillarum infection on diploid and triploid Chinook salmon (Oncorhynchus tshawytscha) was investigated to identify if there was any significant immune regulation by microRNAs (miRNA). Small RNAs from hindgut, head kidney, and spleen were sequenced to determine if miRNA transcript abundance was altered due to ploidy and infection in nine-month old full-sibling diploids and triploids. All three tissues had differentially expressed miRNA prior to infection, indicating subtle changes in epigenetic regulation due to increased ploidy. Additionally, miRNA were altered by infection, but there was only a difference in spleen miRNA expression between diploids and triploids at three days of infection. Furthermore, one miRNA (ssa-miR-2188-3p) was confirmed as having an altered response to infection in triploids compared to diploids, implicating potential immune dysregulation due to increased ploidy. The miRNAs identified in this study are predicted to target immune pathways, providing evidence for their importance in regulating responses to pathogens. This study is the first to investigate how increased ploidy alters miRNA expression in response to infection. Additionally, it provides evidence for epigenetic dysregulation in triploid fishes, which may contribute to their poor performance in response to stress.

Canalization of Phenotypes-When the Transcriptome is Constantly but Weakly Perturbed

Recent studies have increasingly pointed to microRNAs (miRNAs) as the agent of gene regulatory network (GRN) stabilization as well as developmental canalization against constant but small environmental perturbations. To analyze mild perturbations, we construct a Dicer-1 knockdown line (dcr-1 KD) in Drosophila that modestly reduces all miRNAs by, on average, ∼20%. The defining characteristic of stabilizers is that, when their capacity is compromised, GRNs do not change their short-term behaviors. Indeed, even with such broad reductions across all miRNAs, the changes in the transcriptome are very modest during development in stable environment. By comparison, broad knockdowns of other regulatory genes (esp. transcription factors) by the same method should lead to drastic changes in the GRNs. The consequence of destabilization may thus be in long-term development as postulated by the theory of canalization. Flies with modest miRNA reductions may gradually deviate from the developmental norm, resulting in late-stage failures such as shortened longevity. In the optimal culture condition, the survival to adulthood is indeed normal in the dcr-1 KD line but, importantly, adult longevity is reduced by ∼90%. When flies are stressed by high temperature, dcr-1 KD induces lethality earlier in late pupation and, as the perturbations are shifted earlier, the affected stages are shifted correspondingly. Hence, in late stages of development with deviations piling up, GRN would be increasingly in need of stabilization. In conclusion, miRNAs appear to be a solution to weak but constant environmental perturbations.

Challenging Cellular Homeostasis: Spatial and Temporal Regulation of miRNAs

Mature microRNAs (miRNAs) are single-stranded non-coding RNA (ncRNA) molecules that act in post-transcriptional regulation in animals and plants. A mature miRNA is the end product of consecutive, highly regulated processing steps of the primary miRNA transcript. Following base-paring of the mature miRNA with its mRNA target, translation is inhibited, and the targeted mRNA is degraded. There are hundreds of miRNAs in each cell that work together to regulate cellular key processes, including development, differentiation, cell cycle, apoptosis, inflammation, viral infection, and more. In this review, we present an overlooked layer of cellular regulation that addresses cell dynamics affecting miRNA accessibility. We discuss the regulation of miRNA local storage and translocation among cell compartments. The local amounts of the miRNAs and their targets dictate their actual availability, which determines the ability to fine-tune cell responses to abrupt or chronic changes. We emphasize that changes in miRNA storage and compactization occur under induced stress and changing conditions. Furthermore, we demonstrate shared principles on cell physiology, governed by miRNA under oxidative stress, tumorigenesis, viral infection, or synaptic plasticity. The evidence presented in this review article highlights the importance of spatial and temporal miRNA regulation for cell physiology. We argue that limiting the research to mature miRNAs within the cytosol undermines our understanding of the efficacy of miRNAs to regulate cell fate under stress conditions.

Systematic characterization of seed overlap microRNA cotargeting associated with lupus pathogenesis

BACKGROUND

Combinatorial gene regulation by multiple microRNAs (miRNAs) is widespread and closely spaced target sites often act cooperatively to achieve stronger repression ("neighborhood" miRNA cotargeting). While miRNA cotarget sites are suggested to be more conserved and implicated in developmental control, the pathological significance of miRNA cotargeting remains elusive.

RESULTS

Here, we report the pathogenic impacts of combinatorial miRNA regulation on inflammation in systemic lupus erythematosus (SLE). In the SLE mouse model, we identified the downregulation of two miRNAs, miR-128 and miR-148a, by TLR7 stimulation in plasmacytoid dendritic cells. Functional analyses using human cell lines demonstrated that miR-128 and miR-148a additively target KLF4 via extensively overlapping target sites ("seed overlap" miRNA cotargeting) and suppress the inflammatory responses. At the transcriptome level, "seed overlap" miRNA cotargeting increases susceptibility to downregulation by two miRNAs, consistent with additive but not cooperative recruitment of two miRNAs. Systematic characterization further revealed that extensive "seed overlap" is a prevalent feature among broadly conserved miRNAs. Highly conserved target sites of broadly conserved miRNAs are largely divided into two classes-those conserved among eutherian mammals and from human to Coelacanth, and the latter, including KLF4-cotargeting sites, has a stronger association with both "seed overlap" and "neighborhood" miRNA cotargeting. Furthermore, a deeply conserved miRNA target class has a higher probability of haplo-insufficient genes.

CONCLUSIONS

Our study collectively suggests the complexity of distinct modes of miRNA cotargeting and the importance of their perturbations in human diseases.

A microRNA Cluster Controls Fat Cell Differentiation and Adipose Tissue Expansion By Regulating SNCG

The H19X-encoded miR-424(322)/503 cluster regulates multiple cellular functions. Here, it is reported for the first time that it is also a critical linchpin of fat mass expansion. Deletion of this miRNA cluster in mice results in obesity, while increasing the pool of early adipocyte progenitors and hypertrophied adipocytes. Complementary loss and gain of function experiments and RNA sequencing demonstrate that miR-424(322)/503 regulates a conserved genetic program involved in the differentiation and commitment of white adipocytes. Mechanistically, it is demonstrated that miR-424(322)/503 targets γ-Synuclein (SNCG), a factor that mediates this program rearrangement by controlling metabolic functions in fat cells, allowing adipocyte differentiation and adipose tissue enlargement. Accordingly, diminished miR-424(322) in mice and obese humans co-segregate with increased SNCG in fat and peripheral blood as mutually exclusive features of obesity, being normalized upon weight loss. The data unveil a previously unknown regulatory mechanism of fat mass expansion tightly controlled by the miR-424(322)/503 through SNCG.

MicroRNA ssa-mir-196a-4 deceases lgr8 expression in testis development of Chinese tongue sole (Cynoglossus semilaevis)

MicroRNAs (miRNAs) contribute to gonadal development in animals. However, there is little information about miRNA regulation function involved in gonadal development in fish. Our group previously identified sex-related miRNAs of Chinese tongue sole (Cynoglossus semilaevis) during sex determination and differentiation by small RNA sequencing. In the present study, we characterized ssa-mir-196a-4 and its expression in testis and verified its interaction with lgr8. miRNA ssa-mir-196a-4 precursor was predicted to have a typical hairpin structure and highly conserved among various fish species. Fluorescence in situ hybridization (FISH) of ssa-mir-196a-4 in the testis of Chinese tongue sole showed that it is mainly expressed in the cytoplasm of Sertoli cells. We determined that ssa-mir-196a-4 interacted with lgr8 by bioinformatics analysis using miRanda software. According to the dual-luciferase gene reporter assay, lgr8 is a direct target of ssa-mir-196a-4. Overexpression of ssa-mir-196a-4 in the cells of the testis cell line of Chinese tongue sole decreased the expression levels of lgr8 messenger RNA (mRNA) and protein by targeting its coding sequence (CDS) region. These results suggest that ssa-mir-196a-4 acts as a post-transcriptional regulator of lgr8 and plays an important role in developing testes of Chinese tongue sole.

Variation of microRNA expression in the human placenta driven by population identity and sex of the newborn

BACKGROUND

Analysis of lymphocyte cell lines revealed substantial differences in the expression of mRNA and microRNA (miRNA) among human populations. The extent of such population-associated differences in actual human tissues remains largely unexplored. The placenta is one of the few solid human tissues that can be collected in substantial numbers in a controlled manner, enabling quantitative analysis of transient biomolecules such as RNA transcripts. Here, we analyzed microRNA (miRNA) expression in human placental samples derived from 36 individuals representing four genetically distinct human populations: African Americans, European Americans, South Asians, and East Asians. All samples were collected at the same hospital following a unified protocol, thus minimizing potential biases that might influence the results.

RESULTS

Sequence analysis of the miRNA fraction yielded 938 annotated and 70 novel miRNA transcripts expressed in the placenta. Of them, 82 (9%) of annotated and 11 (16%) of novel miRNAs displayed quantitative expression differences among populations, generally reflecting reported genetic and mRNA-expression-based distances. Several co-expressed miRNA clusters stood out from the rest of the population-associated differences in terms of miRNA evolutionary age, tissue-specificity, and disease-association characteristics. Among three non-environmental influenced demographic parameters, the second largest contributor to miRNA expression variation after population was the sex of the newborn, with 32 miRNAs (3% of detected) exhibiting significant expression differences depending on whether the newborn was male or female. Male-associated miRNAs were evolutionarily younger and correlated inversely with the expression of target mRNA involved in neuron-related functions. In contrast, both male and female-associated miRNAs appeared to mediate different types of hormonal responses. Demographic factors further affected reported imprinted expression of 66 placental miRNAs: the imprinting strength correlated with the mother's weight, but not height.

CONCLUSIONS

Our results showed that among 12 assessed demographic variables, population affiliation and fetal sex had a substantial influence on miRNA expression variation among human placental samples. The effect of newborn-sex-associated miRNA differences further led to expression inhibition of the target genes clustering in specific functional pathways. By contrast, population-driven miRNA differences might mainly represent neutral changes with minimal functional impacts.

Comparative epigenetics in animal physiology: An emerging frontier

The unprecedented access to annotated genomes now facilitates the investigation of the molecular basis of epigenetic phenomena in phenotypically diverse animals. In this critical review, we describe the roles of molecular epigenetic mechanisms in regulating mitotically and meiotically stable spatiotemporal gene expression, phenomena that provide the molecular foundation for the intra-, inter-, and trans-generational emergence of physiological phenotypes. By focusing principally on emerging comparative epigenetic roles of DNA-level and transcriptome-level epigenetic mark dynamics in the emergence of phenotypes, we highlight the relationship between evolutionary conservation and innovation of specific epigenetic pathways, and their interplay as a priority for future study. This comparative approach is expected to significantly advance our understanding of epigenetic phenomena, as animals show a diverse array of strategies to epigenetically modify physiological responses. Additionally, we review recent technological advances in the field of molecular epigenetics (single-cell epigenomics and transcriptomics and editing of epigenetic marks) in order to (1) investigate environmental and endogenous factor dependent epigenetic mark dynamics in an integrative manner; (2) functionally test the contribution of specific epigenetic marks for animal phenotypes via genome and transcript-editing tools. Finally, we describe advantages and limitations of emerging animal models, which under the Krogh principle, may be particularly useful in the advancement of comparative epigenomics and its potential translational applications in animal science, ecotoxicology, ecophysiology, climate change science and wild-life conservation, as well as organismal health.

Dynamic landscape and evolution of m6A methylation in human

m⁶A is a prevalent internal modification in mRNAs and has been linked to the diverse effects on mRNA fate. To explore the landscape and evolution of human m⁶A, we generated 27 m⁶A methylomes across major adult tissues. These data reveal dynamic m⁶A methylation across tissue types, uncover both broadly or tissue-specifically methylated sites, and identify an unexpected enrichment of m⁶A methylation at non-canonical cleavage sites. A comparison of fetal and adult m⁶A methylomes reveals that m⁶A preferentially occupies CDS regions in fetal tissues. Moreover, the m⁶A sub-motifs vary between fetal and adult tissues or across tissue types. From the evolutionary perspective, we uncover that the selection pressure on m⁶A sites varies and depends on their genic locations. Unexpectedly, we found that ∼40% of the 3′UTR m⁶A sites are under negative selection, which is higher than the evolutionary constraint on miRNA binding sites, and much higher than that on A-to-I RNA modification. Moreover, the recently gained m⁶A sites in human populations are clearly under positive selection and associated with traits or diseases. Our work provides a resource of human m⁶A profile for future studies of m⁶A functions, and suggests a role of m⁶A modification in human evolutionary adaptation and disease susceptibility.

Evolutionary dynamics of microRNA target sites across vertebrate evolution

MicroRNAs (miRNAs) control the abundance of the majority of the vertebrate transcriptome. The recognition sequences, or target sites, for bilaterian miRNAs are found predominantly in the 3' untranslated regions (3'UTRs) of mRNAs, and are amongst the most highly conserved motifs within 3'UTRs. However, little is known regarding the evolutionary pressures that lead to loss and gain of such target sites. Here, we quantify the selective pressures that act upon miRNA target sites. Notably, selective pressure extends beyond deeply conserved binding sites to those that have undergone recent substitutions. Our approach reveals that even amongst ancient animal miRNAs, which exert the strongest selective pressures on 3'UTR sequences, there are striking differences in patterns of target site evolution between miRNAs. Considering only ancient animal miRNAs, we find three distinct miRNA groups, each exhibiting characteristic rates of target site gain and loss during mammalian evolution. The first group both loses and gains sites rarely. The second group shows selection only against site loss, with site gains occurring at a neutral rate, whereas the third loses and gains sites at neutral or above expected rates. Furthermore, mutations that alter the strength of existing target sites are disfavored. Applying our approach to individual transcripts reveals variation in the distribution of selective pressure across the transcriptome and between miRNAs, ranging from strong selection acting on a small subset of targets of some miRNAs, to weak selection on many targets for other miRNAs. miR-20 and miR-30, and many other miRNAs, exhibit broad, deeply conserved targeting, while several other comparably ancient miRNAs show a lack of selective constraint, and a small number, including mir-146, exhibit evidence of rapidly evolving target sites. Our approach adds valuable perspective on the evolution of miRNAs and their targets, and can also be applied to characterize other 3'UTR regulatory motifs.

Chronic Voluntary Ethanol Drinking in Cynomolgus Macaques Elicits Gene Expression Changes in Prefrontal Cortical Area 46

BACKGROUND

Genome-wide profiling to examine brain transcriptional features associated with excessive ethanol (EtOH) consumption has been applied to a variety of species including rodents, nonhuman primates (NHPs), and humans. However, these data were obtained from cross-sectional samples which are particularly vulnerable to individual variation when obtained from small outbred populations typical of human and NHP studies. In the current study, a novel within-subject design was used to examine the effects of voluntary EtOH consumption on prefrontal cortex (PFC) gene expression in a NHP model.

METHODS

Two cohorts of cynomolgus macaques (n = 23) underwent a schedule-induced polydipsia procedure to establish EtOH self-administration followed by 6 months of daily open access to EtOH (4% w/v) and water. Individual daily EtOH intakes ranged from an average of 0.7 to 3.7 g/kg/d. Dorsal lateral PFC area 46 (A46) brain biopsies were collected in EtOH-naïve and control monkeys; contralateral A46 biopsies were collected from the same monkeys following the 6 months of fluid consumption. Gene expression changes were assessed using RNA-Seq paired analysis, which allowed for correction of individual baseline differences in gene expression.

RESULTS

A total of 675 genes were significantly down-regulated following EtOH consumption; these were functionally enriched for immune response, cell adhesion, plasma membrane, and extracellular matrix. A total of 567 genes that were up-regulated following EtOH consumption were enriched in microRNA target sites and included target sites associated with Toll-like receptor pathways. The differentially expressed genes were also significantly enriched in transcription factor binding sites.

CONCLUSIONS

The data presented here are the first to use a longitudinal biopsy strategy to examine how chronic EtOH consumption affects gene expression in the primate PFC. Prominent effects were seen in both cell adhesion and neuroimmune pathways; the latter contained both pro- and antiinflammatory genes. The data also indicate that changes in miRNAs and transcription factors may be important epigenetic regulators of EtOH consumption.

Gene regulatory network stabilized by pervasive weak repressions: microRNA functions revealed by the May-Wigner theory

Food web and gene regulatory networks (GRNs) are large biological networks, both of which can be analyzed using the May-Wigner theory. According to the theory, networks as large as mammalian GRNs would require dedicated gene products for stabilization. We propose that microRNAs (miRNAs) are those products. More than 30% of genes are repressed by miRNAs, but most repressions are too weak to have a phenotypic consequence. The theory shows that (i) weak repressions cumulatively enhance the stability of GRNs, and (ii) broad and weak repressions confer greater stability than a few strong ones. Hence, the diffuse actions of miRNAs in mammalian cells appear to function mainly in stabilizing GRNs. The postulated link between mRNA repression and GRN stability can be seen in a different light in yeast, which do not have miRNAs. Yeast cells rely on non-specific RNA nucleases to strongly degrade mRNAs for GRN stability. The strategy is suited to GRNs of small and rapidly dividing yeast cells, but not the larger mammalian cells. In conclusion, the May-Wigner theory, supplanting the analysis of small motifs, provides a mathematical solution to GRN stability, thus linking miRNAs explicitly to 'developmental canalization'.

miRNAs responsive to the diabetic microenvironment in the human beta cell line EndoC-βH1 may target genes in the FOXO, HIPPO and Lysine degradation pathways

Altered expression of miRNAs is evident in the islets of diabetic human donors, but the effects of specific aspects of the diabetic microenvironment and identity of gene ontology pathways demonstrating target gene enrichment in response to each is understudied. We assessed changes in the miRNA milieu in response to high/low glucose, hypoxia, dyslipidaemia and inflammatory factors in a humanised EndoC-βH1 beta cell culture system and performed miRPath analysis for each treatment individually. The 10 miRNAs demonstrating the greatest dysregulation across treatments were then independently validated and Gene Set Enrichment Analysis to confirm targeted pathways undertaken. 171 of 392 miRNAs displayed altered expression in response to one or more cellular stressors. miRNA changes were treatment specific, but their target genes were enriched in conserved pathways. 5 miRNAs (miR-136-5p, miR299-5p, miR-454-5p, miR-152 and miR-185) were dysregulated in response to multiple stressors and survived validation in independent samples (p = 0.008, 0.002, 0.012, 0.005 and 0.024 respectively). Target genes of dysregulated miRNAs were clustered into FOXO1, HIPPO and Lysine degradation pathways (p = 0.02, p = 5.84 × 10^-5 and p = 3.00 × 10^-3 respectively). We provide evidence that the diabetic microenvironment may induce changes to the expression of miRNAs targeting genes enriched in pathways involved in cell stress response and cell survival.

miR-15a/-16 Inhibit Angiogenesis by Targeting the Tie2 Coding Sequence: Therapeutic Potential of a miR-15a/16 Decoy System in Limb Ischemia

MicroRNA-15a (miR-15a) and miR-16, which are transcribed from the miR-15a/miR-16-1 cluster, inhibit post-ischemic angiogenesis. MicroRNA (miRNA) binding to mRNA coding sequences (CDSs) is a newly emerging mechanism of gene expression regulation. We aimed to (1) identify new mediators of the anti-angiogenic action of miR-15a and -16, (2) develop an adenovirus (Ad)-based miR-15a/16 decoy system carrying a luciferase reporter (Luc) to both sense and inhibit miR-15a/16 activity, and (3) investigate Ad.Luc-Decoy-15a/16 therapeutic potential in a mouse limb ischemia (LI) model. LI increased miR-15a and -16 expression in mouse muscular endothelial cells (ECs). The miRNAs also increased in cultured human umbilical vein ECs (HUVECs) exposed to serum starvation, but not hypoxia. Using bioinformatic tools and luciferase activity assays, we characterized miR-15a and -16 binding to Tie2 CDS. In HUVECs, miR-15a or -16 overexpression reduced Tie2 at the protein, but not the mRNA, level. Conversely, miR-15a or -16 inhibition improved angiogenesis in a Tie2-dependent manner. Local Ad.Luc-Decoy-15a/16 delivery increased Tie2 levels in ischemic skeletal muscle and improved post-LI angiogenesis and perfusion recovery, with reduced toe necrosis. Bioluminescent imaging (in vivo imaging system [IVIS]) provided evidence that the Ad.Luc-Decoy-15a/16 system responds to miR-15a/16 increases. In conclusion, we have provided novel mechanistic evidence of the therapeutic potential of local miR-15a/16 inhibition in LI.

Serum miR-182 is a predictive biomarker for dichotomization of risk of hepatocellular carcinoma in rats

Exploration of animal models leads to discoveries that can reveal candidate biomarkers for translation to human populations. Herein, a model of hepatocarcinogenesis and protection was used in which rats treated with aflatoxin (AFB₁ ) daily for 28 days (200 µg/kg BW) developed tumors compared with rats completely protected from tumors by concurrent administration of the chemoprotective agent, 1-[2-cyano-3-,12-dioxooleana-1,9(11)-dien-28-oyl]imidazole (CDDO-Im). Differential expression of miRNAs in tumors (AFB₁ ) and nontumor (AFB₁  + CDDO-Im) bearing livers and their levels in sera over the life-course of the animals was determined. miRNA transcriptome analysis identified 17 miRNAs significantly upregulated at greater than five-fold in the tumors. The ten most dysregulated miRNAs judged by fold-change and biological significance were selected for further study, including liver-specific miR-122-5p. Validation of sequencing results by real-time PCR confirmed the upregulation of the majority of these miRNAs in tumors, including miR-182, as well as miR-224-5p as the most dysregulated of these miRNAs (over 400-fold). The longitudinal analysis of levels of miR-182 in sera demonstrated significant and persistent increases (5.13-fold; 95% CI: 4.59-5.70). The increase in miR-182 was detected months before any clinical symptoms were present in the animals. By the terminal time point of the study, in addition to elevated levels of serum miR-182, serum miR-122-5p was also found to be increased (>1.5-fold) in animals that developed hepatocarcinomas. Thus, using the data from an unbiased discovery approach of the tissue findings, serum miR-182 was found to track across the complex, multistage process of hepatocarcinogenesis opening an opportunity for translation to human populations.

Identification and expression of microRNAs in european eels Anguilla anguilla from two natural sites with different pollution levels

MicroRNAs (miRNAs) are a class of small non-coding RNA that control multiple biological processes through negative post-transcriptional regulation of gene expression. Recently a role of miRNAs in the response of aquatic organisms to environmental toxicants emerged. Toxicant-induced changes in miRNA expression might then represent novel biomarkers to evaluate the health status of these organisms. In this study, we aimed to identify the miRNA repertoire in the liver of the European eel Anguilla anguilla and to compare their differential expression between a polluted site located in the Gironde Estuary and a pristine site in Arcachon Bay (France). A total of 299 mature miRNAs were identified. In polluted water, 19 miRNAs were up-regulated and 22 were down-regulated. We predicted that these differentially expressed miRNAs could target 490 genes that were involved in ribosome biogenesis, response to hormones, response to chemical and chromatin modification. Moreover, we observed only few examples (29) of negative correlation between the expression levels of miRNAs and their targets suggesting that, in the system studied, miRNAs might not only regulate gene expression directly by degrading mRNA but also by inhibiting protein translation or by regulating other epigenetic processes. This study is the first example of in situ investigation of the role of miRNAs in the response of a fish species to water quality. Our findings provide new insights into the involvement of epigenetic mechanisms in the response of animals chronically exposed to pollution and pave the way for the utilization of miRNAs in aquatic ecotoxicology.

Social status regulates the hepatic miRNAome in rainbow trout: Implications for posttranscriptional regulation of metabolic pathways

Juvenile rainbow trout develop social hierarchies when held in dyads, and the development of socially subordinate (SS) and social dominance (SD) phenotypes in this context has been linked to specific changes in the hepatic energy metabolism of all major macronutrients. Following our recently reported finding that transcript abundance of drosha, a key component of the microRNA (miRNA) biogenesis pathway, is increased in paired juvenile rainbow trout irrespective of social status compared to socially isolated (SI) controls, we here determined global changes of the hepatic miRNA pathway genes in detail at the transcript and protein level. Both SD and SS rainbow trout exhibited increased Ago2 protein abundance compared to SI rainbow trout, suggesting that hepatic miRNA function is increased in rainbow trout maintained in dyads. Given the well-described differences in hepatic intermediary metabolism between SD and SS rainbow trout, and the important role of miRNAs in the posttranscriptional regulation of metabolic pathways, we also identified changes in hepatic miRNA abundance between SS and SD rainbow trout using small RNA next generation sequencing. We identified a total of 24 differentially regulated miRNAs, with 15 miRNAs that exhibited increased expression, and 9 miRNAs that exhibited decreased expression in the liver of SS trout compared to SD trout. To identify potential miRNA-dependent posttranscriptional regulatory pathways important for social status-dependent regulation of hepatic metabolism in rainbow trout, we used an in silico miRNA target prediction and pathway enrichment approach. We identified enrichment for pathways related to metabolism of carbohydrates, lipids and proteins in addition to organelle-specific processes involved in energy metabolism, especially mitochondrial fusion and fission. Select predicted miRNA-mRNA target pairs within these categories were quantitatively analyzed by real-time RT-PCR to validate candidates for future studies that will probe the functional metabolic roles of specific hepatic miRNAs in the development of SD and SS metabolic phenotypes.

The evolutionary dynamics of microRNAs in domestic mammals

MiRNAs are crucial regulators of gene expression found across both the plant and animal kingdoms. While the number of annotated miRNAs deposited in miRBase has greatly increased in recent years, few studies provided comparative analyses across sets of related species, or investigated the role of miRNAs in the evolution of gene regulation. We generated small RNA libraries across 5 mammalian species (cow, dog, horse, pig and rabbit) from 4 different tissues (brain, heart, kidney and testis). We identified 1676 miRBase and 413 novel miRNAs by manually curating the set of computational predictions obtained from miRCat and miRDeep2. Our dataset spanning five species has enabled us to investigate the molecular mechanisms and selective pressures driving the evolution of miRNAs in mammals. We highlight the important contributions of intronic sequences (366 orthogroups), duplication events (135 orthogroups) and repetitive elements (37 orthogroups) in the emergence of new miRNA loci. We use this framework to estimate the patterns of gains and losses across the phylogeny, and observe high levels of miRNA turnover. Additionally, the identification of lineage-specific losses enables the characterisation of the selective constraints acting on the associated target sites. Compared to the miRBase subset, novel miRNAs tend to be more tissue specific. 20 percent of novel orthogroups are restricted to the brain, and their target repertoires appear to be enriched for neuron activity and differentiation processes. These findings may reflect an important role for young miRNAs in the evolution of brain expression plasticity. Many seed sequences appear to be specific to either the cow or the dog. Analyses on the associated targets highlight the presence of several genes under artificial positive selection, suggesting an involvement of these miRNAs in the domestication process. Altogether, we provide an overview on the evolutionary mechanisms responsible for miRNA turnover in 5 domestic species, and their possible contribution to the evolution of gene regulation.

miRWalk: An online resource for prediction of microRNA binding sites

miRWalk is an open-source platform providing an intuitive interface that generates predicted and validated miRNA-binding sites of known genes of human, mouse, rat, dog and cow. The core of miRWalk is the miRNA target site prediction with the random-forest-based approach software TarPmiR searching the complete transcript sequence including the 5’-UTR, CDS and 3’-UTR. Moreover, it integrates results other databases with predicted and validated miRNA-target interactions. The focus is set on a modular design and extensibility as well as a fast update cycle. The database is available using Python, MySQL and HTML/Javascript Database URL: http://mirwalk.umm.uni-heidelberg.de.

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.

Direct measurement of pervasive weak repression by microRNAs and their role at the network level

BACKGROUND

A gene regulatory network (GRN) comprises many weak links that are often regulated by microRNAs. Since miRNAs rarely repress their target genes by more than 30%, doubts have been expressed about the biological relevance of such weak effects. These doubts raise the possibility of under-estimation as miRNA repression is usually estimated indirectly from equilibrium expression levels.

RESULTS

To measure miRNA repression directly, we inhibited transcript synthesis in Drosophila larvae and collected time-course data on mRNA abundance, the decline of which reflects transcript degradation. The rate of target degradation in the absence of miR310s, a moderately expressed miRNA family, was found to decrease by 5 to 15%. A conventional analysis that does not remove transcript synthesis yields an estimate of 6.5%, within the range of the new estimates. These data permit further examinations of the repression mechanisms by miRNAs including seed matching types, APA (alternative polyadenylation) sites, effects of other highly-expressed miRNAs and the length of 3'UTR. Our direct measurements suggest the latter two factors have a measurable effect on decay rate.

CONCLUSION

The direct measurement confirms pervasive weak repression by miRNAs, supporting the conclusions based on indirect assays. The confirmation suggests that this weak repression may indeed be miRNAs' main function. In this context, we discuss the recent proposal that weak repression is "cumulatively powerful" in stabilizing GRNs.

Evolution of New miRNAs and Cerebro-Cortical Development

The noncoding portion of the genome, including microRNAs, has been fertile evolutionary soil for cortical development in primates. A major contribution to cortical expansion in primates is the generation of novel precursor cell populations. Because miRNA expression profiles track closely with cell identity, it is likely that numerous novel microRNAs have contributed to cellular diversity in the brain. The tools to determine the genomic context within which novel microRNAs emerge and how they become integrated into molecular circuitry are now in hand.

Epigenetics in teleost fish: From molecular mechanisms to physiological phenotypes

While the field of epigenetics is increasingly recognized to contribute to the emergence of phenotypes in mammalian research models across different developmental and generational timescales, the comparative biology of epigenetics in the large and physiologically diverse vertebrate infraclass of teleost fish remains comparatively understudied. The cypriniform zebrafish and the salmoniform rainbow trout and Atlantic salmon represent two especially important teleost orders, because they offer the unique possibility to comparatively investigate the role of epigenetic regulation in 3R and 4R duplicated genomes. In addition to their sequenced genomes, these teleost species are well-characterized model species for development and physiology, and therefore allow for an investigation of the role of epigenetic modifications in the emergence of physiological phenotypes during an organism's lifespan and in subsequent generations. This review aims firstly to describe the evolution of the repertoire of genes involved in key molecular epigenetic pathways including histone modifications, DNA methylation and microRNAs in zebrafish, rainbow trout, and Atlantic salmon, and secondly, to discuss recent advances in research highlighting a role for molecular epigenetics in shaping physiological phenotypes in these and other teleost models. Finally, by discussing themes and current limitations of the emerging field of teleost epigenetics from both theoretical and technical points of view, we will highlight future research needs and discuss how epigenetics will not only help address basic research questions in comparative teleost physiology, but also inform translational research including aquaculture, aquatic toxicology, and human disease.

Rainbow trout exposed to benzo[a]pyrene yields conserved microRNA binding sites in DNA methyltransferases across 500 million years of evolution

The objective of this study was to examine the regulation of DNA methylation following acute (24 h) and prolonged (14 d) exposure to low (1 ng/L) and high (10 ng/L) benzo[a]pyrene. However, with the recent release of the rainbow trout genome, we were able to conduct a more detailed analysis regarding the regulation of the enzymes involved in DNA methylation; DNA methyltransferases (DNMTs). Bioinformatic approaches were used to identify candidate microRNA (miRNA) that potentially bind to the DNMT1 and DNMT3a 3′UTR. Results indicated a significant decrease in global methylation in both liver and muscle, with an associated decrease in DNA methyltransferase activity and DNMT3a transcript abundance. There was a significant increase in one specific candidate miRNA (miR29a) that was predicted to bind to DNMT3a. Taking a comparative genomics approach, the binding sites of miR29a to the DNMT3a 3′UTR was compared across species, spanning fish to mammals, and revealed a highly conserved binding motif that has been maintained since the vertebrate ancestor, approximately 500 million years ago. This research establishes that miRNA act as an essential mediator between the environment and DNA methylation patterns via DNMTs, which is further confirmed by a genomic regulatory mechanism that has been deeply conserved throughout evolution.

Argonaute CLIP Defines a Deregulated miR-122-Bound Transcriptome that Correlates with Patient Survival in Human Liver Cancer

MicroRNA-122, an abundant and conserved liver-specific miRNA, regulates hepatic metabolism and functions as a tumor suppressor, yet systematic and direct biochemical elucidation of the miR-122 target network remains incomplete. To this end, we performed Argonaute crosslinking immunoprecipitation (Argonaute [Ago]-CLIP) sequencing in miR-122 knockout and control mouse livers, as well as in matched human hepatocellular carcinoma (HCC) and benign liver tissue to identify miRNA target sites transcriptome-wide in two species. We observed a majority of miR-122 binding on 3' UTRs and coding exons followed by extensive binding to other genic and non-genic sites. Motif analysis of miR-122-dependent binding revealed a G-bulged motif in addition to canonical motifs. A large number of miR-122 targets were found to be species specific. Upregulation of several common mouse and human targets, most notably BCL9, predicted survival in HCC patients. These results broadly define the molecular consequences of miR-122 downregulation in hepatocellular carcinoma.

Serum microRNAs in buffalo cows: Potential biomarkers of pregnancy

MiRNAs (microRNA) constitute a large family of single-stranded, non-coding small RNAs. Although the functions and target genes of most miRNAs are still unknown, it has been well established that they are involved in embryogenesis, organogenesis and neonatal birth. In recent decades, interest in buffalo breeding has largely increased worldwide thus focusing the attention on this species as a dairy purpose animal. Problems related to long calving interval, late puberty and seasonal anestrus hamper reproductive efficiency in this species. Early pregnancy diagnosis is important to shorten the calving interval and increase lifetime production on dairy animals. MicroRNAs have recently emerged as key molecules in fertilization of several species even though in buffalo, few previous studies have investigated miRNAs. The aim of this research was to identify the best miRNA reference in serum among miR-191, miR-25-3p, SNORD44, and SNORD48. Consequently, assess the expression levels of miR-103, miR-200b, miR-301a, miR-423-5p, miR-375, miR-451 and miR-452 involved in buffalo progesterone-maturation oocyte and pregnancy. Interestingly, we found that all the miRNAs analyzed at 40days after artificial insemination were differentially expressed among pregnant and non-pregnant buffalo. Our research may be a first step for gain further insight in the biological function of circulating miRNAs in this species.

Unveiling the Impact of the Genomic Architecture on the Evolution of Vertebrate microRNAs

Eukaryotic genomes frequently exhibit interdependency between transcriptional units, as evidenced by regions of high gene density. It is well recognized that vertebrate microRNAs (miRNAs) are usually embedded in those regions. Recent work has shown that the genomic context is of utmost importance to determine miRNA expression in time and space, thus affecting their evolutionary fates over long and short terms. Consequently, understanding the inter- and intraspecific changes on miRNA genomic architecture may bring novel insights on the basic cellular processes regulated by miRNAs, as well as phenotypic evolution and disease-related mechanisms.

Prediction of human miRNA target genes using computationally reconstructed ancestral mammalian sequences

MicroRNAs (miRNA) are short single-stranded RNA molecules derived from hairpin-forming precursors that play a crucial role as post-transcriptional regulators in eukaryotes and viruses. In the past years, many microRNA target genes (MTGs) have been identified experimentally. However, because of the high costs of experimental approaches, target genes databases remain incomplete. Although several target prediction programs have been developed in the recent years to identify MTGs in silico, their specificity and sensitivity remain low. Here, we propose a new approach called MirAncesTar, which uses ancestral genome reconstruction to boost the accuracy of existing MTGs prediction tools for human miRNAs. For each miRNA and each putative human target UTR, our algorithm makes uses of existing prediction tools to identify putative target sites in the human UTR, as well as in its mammalian orthologs and inferred ancestral sequences. It then evaluates evidence in support of selective pressure to maintain target site counts (rather than sequences), accounting for the possibility of target site turnover. It finally integrates this measure with several simpler ones using a logistic regression predictor. MirAncesTar improves the accuracy of existing MTG predictors by 26% to 157%. Source code and prediction results for human miRNAs, as well as supporting evolutionary data are available at http://cs.mcgill.ca/∼blanchem/mirancestar.

Discovering MicroRNA-Regulatory Modules in Multi-Dimensional Cancer Genomic Data: A Survey of Computational Methods

MicroRNAs (miRs) are small single-stranded noncoding RNA that function in RNA silencing and post-transcriptional regulation of gene expression. An increasing number of studies have shown that miRs play an important role in tumorigenesis, and understanding the regulatory mechanism of miRs in this gene regulatory network will help elucidate the complex biological processes at play during malignancy. Despite advances, determination of miR–target interactions (MTIs) and identification of functional modules composed of miRs and their specific targets remain a challenge. A large amount of data generated by high-throughput methods from various sources are available to investigate MTIs. The development of data-driven tools to harness these multi-dimensional data has resulted in significant progress over the past decade. In parallel, large-scale cancer genomic projects are allowing new insights into the commonalities and disparities of miR–target regulation across cancers. In the first half of this review, we explore methods for identification of pairwise MTIs, and in the second half, we explore computational tools for discovery of miR-regulatory modules in a cancer-specific and pan-cancer context. We highlight strengths and limitations of each of these tools as a practical guide for the computational biologists.

LimiTT: link miRNAs to targets

Background

MicroRNAs (miRNAs) impact various biological processes within animals and plants. They complementarily bind target mRNAs, effecting a post-transcriptional negative regulation on mRNA level. The investigation of miRNA target interactions (MTIs) by high throughput screenings is challenging, as frequently used in silico target prediction tools are prone to emit false positives. This issue is aggravated for niche model organisms, where validated miRNAs and MTIs both have to be transferred from well described model organisms. Even though DBs exist that contain experimentally validated MTIs, they are limited in their search options and they utilize different miRNA and target identifiers.

Results

The implemented pipeline LimiTT integrates four existing DBs containing experimentally validated MTIs. In contrast to other cumulative databases (DBs), LimiTT includes MTI data of 26 species. Additionally, the pipeline enables the identification and enrichment analysis of MTIs with and without species specificity based on dynamic quality criteria. Multiple tabular and graphical outputs are generated to permit the detailed assessment of results.

Conclusion

Our freely available web-based pipeline LimiTT (https://bioinformatics.mpi-bn.mpg.de/) is optimized to determine MTIs with and without species specification. It links miRNAs and/or putative targets with high granularity. The integrated mapping to homologous target identifiers enables the identification of MTIs not only for standard models, but for niche model organisms as well.

Electronic supplementary material

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

Host gene constraints and genomic context impact the expression and evolution of human microRNAs

Increasing evidence has shown that recent miRNAs tend to emerge within coding genes. Here we conjecture that human miRNA evolution is tightly influenced by the genomic context, especially by host genes. Our findings show a preferential emergence of intragenic miRNAs within old genes. We found that miRNAs within old host genes are significantly more broadly expressed than those within young ones. Young miRNAs within old genes are more broadly expressed than their intergenic counterparts, suggesting that young miRNAs have an initial advantage by residing in old genes, and benefit from their hosts' expression control and from the exposure to diverse cellular contexts and target genes. Our results demonstrate that host genes may provide stronger expression constraints to intragenic miRNAs in the long run. We also report associated functional implications, highlighting the genomic context and host genes as driving factors for the expression and evolution of human miRNAs.

Recent miRNAs tend to emerge within coding genes. Here, by analysing miRNA expression data from six species and comparing genomes from 13 species, the authors report that host genes may provide stronger expression constraints to intragenic miRNAs in the long run.

General hallmarks of microRNAs in brain evolution and development

MicroRNAs (miRNAs) are endogenous, small non-coding RNA molecules that mediate post-transcriptional gene suppression by incomplete matches with their host mRNAs. In the central nervous system, miRNAs that functionally interact with their target genes constitute a flexible, robust and buffered regulatory network, exerting diverse roles in brain evolution and development. However, distinct variation either in hub miRNA expression levels or patterns may initiate and/or progress various adult-onset nerve-related diseases. In this review, we will summarize the current knowledge about the general hallmarks of brain miRNAs that act as vital determinants in increasingly complicated neural activities. We endeavor to provide a constructive insight into the neuroscience research in the quest to comprehend molecular underpinnings of physiological functions and pathological disorders in central nervous system.

Compartmentalization of functions and predicted miRNA regulation among contiguous regions of the nematode intestine

The intestine of parasitic nematodes has proven an important target for therapies aimed at prevention and treatment of diseases caused by these pathogens in humans, animals and plants. We have developed a unique research model with the intestine of Ascaris suum, the large round worm of swine and humans, that will enhance biological research on this tissue. To expand utility of this model, we quantitatively compared expression of 15,382 coding RNAs and 277 noncoding, micro RNAs (miRNAs) among 3 contiguous regions of the adult A. suum intestine. Differentially expressed transcripts were identified among regions, with the largest number expressed at significantly higher levels in the anterior region, identifying this region as the most functionally unique compared to middle and posterior regions. We further identified 64 exon splice variants (from 47 genes) that are differentially expressed among these regions. A total of 2,063 intestinal mRNA transcripts were predicted to be targeted by intestinal miRNA, and negative correlation coefficients for miRNA:mRNA abundances predicted 22 likely influential miRNAs and 503 likely associated miRNA:mRNA pairs. A. suum intestinal miRNAs were identified that are conserved with intestinal miRNAs from C. elegans (10 mature sequences and 13 seed sequences conserved), and prospective intestinal miRNAs from the murine gastrointestinal nematode, Heligmosomoides polygyrus (5 mature and 11 seeds). Most of the conserved intestinal miRNAs were also high abundance miRNAs. The data provide the most comprehensive compilation of constitutively and differentially expressed genes along the length of the intestine for any nematode species. The information will guide prospective development of many hypotheses on nematode intestinal functions encoded by mRNAs, miRNAs and interactions between these RNA populations.

Micromanaging metabolism-a role for miRNAs in teleost energy metabolism

MicroRNAs (miRNAs) are small, non-protein coding RNA sequences, which are found in most eukaryotes. Since their initial discovery, miRNAs have emerged as important regulators of many biological processes. One of the most important processes profoundly regulated by miRNAs is energy metabolism. Traditionally, metabolic functions of miRNAs have been studied in genome-sequenced mammalian organisms, especially the mouse model. However, partially driven by commercial interest in aquaculture, increasingly feasible large-scale molecular techniques have resulted in the characterization of miRNA repertoires, and importantly, several genome sequences of several (commercially important) teleost species, which also hold important roles as research models in the comparative physiology of energy metabolism. This review aims to introduce the recent advances in miRNA research in teleost fish and to describe the current knowledge of miRNA function in teleost energy metabolism. The most pressing research needs and questions to determine metabolic roles of miRNAs in teleost models are presented, as well as applicable technical approaches and current bottlenecks. Rainbow trout, which possess the advantages of newly available molecular tools and a long history as comparative research model in teleost energy metabolism, are discussed as a promising research model to address these questions.

Integrated analyses of zebrafish miRNA and mRNA expression profiles identify miR-29b and miR-223 as potential regulators of optic nerve regeneration

Background

Unlike mammals, zebrafish have the ability to regenerate damaged parts of their central nervous system (CNS) and regain functionality of the affected area. A better understanding of the molecular mechanisms involved in zebrafish regeneration may therefore provide insight into how CNS repair might be induced in mammals. Although many studies have described differences in gene expression in zebrafish during CNS regeneration, the regulatory mechanisms underpinning the differential expression of these genes have not been examined.

Results

We used microarrays to analyse and integrate the mRNA and microRNA (miRNA) expression profiles of zebrafish retina after optic nerve crush to identify potential regulatory mechanisms that underpin central nerve regeneration. Bioinformatic analysis identified 3 miRNAs and 657 mRNAs that were differentially expressed after injury. We then combined inverse correlations between our miRNA expression and mRNA expression, and integrated these findings with target predictions from TargetScan Fish to identify putative miRNA-gene target pairs. We focused on two over-expressed miRNAs (miR-29b and miR-223), and functionally validated seven of their predicted gene targets using RT-qPCR and luciferase assays to confirm miRNA-mRNA binding. Gene ontology analysis placed the miRNA-regulated genes (eva1a, layna, nefmb, ina, si:ch211-51a6.2, smoc1, sb:cb252) in key biological processes that included cell survival/apoptosis, ECM-cytoskeleton signaling, and heparan sulfate proteoglycan binding,

Conclusion

Our results suggest a key role for miR-29b and miR-223 in zebrafish regeneration. The identification of miRNA regulation in a zebrafish injury model provides a framework for future studies in which to investigate not only the cellular processes required for CNS regeneration, but also how these mechanisms might be regulated to promote successful repair and return of function in the injured mammalian brain.

Electronic supplementary material

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

Two Virus-Induced MicroRNAs Known Only from Teleost Fishes Are Orthologues of MicroRNAs Involved in Cell Cycle Control in Humans

MicroRNAs (miRNAs) are ~22 base pair-long non-coding RNAs which regulate gene expression in the cytoplasm of eukaryotic cells by binding to specific target regions in mRNAs to mediate transcriptional blocking or mRNA cleavage. Through their fundamental roles in cellular pathways, gene regulation mediated by miRNAs has been shown to be involved in almost all biological phenomena, including development, metabolism, cell cycle, tumor formation, and host-pathogen interactions. To address the latter in a primitive vertebrate host, we here used an array platform to analyze the miRNA response in rainbow trout (Oncorhynchus mykiss) following inoculation with the virulent fish rhabdovirus Viral hemorrhagic septicaemia virus. Two clustered miRNAs, miR-462 and miR-731 (herein referred to as miR-462 cluster), described only in teleost fishes, were found to be strongly upregulated, indicating their involvement in fish-virus interactions. We searched for homologues of the two teleost miRNAs in other vertebrate species and investigated whether findings related to ours have been reported for these homologues. Gene synteny analysis along with gene sequence conservation suggested that the teleost fish miR-462 and miR-731 had evolved from the ancestral miR-191 and miR-425 (herein called miR-191 cluster), respectively. Whereas the miR-462 cluster locus is found between two protein-coding genes (intergenic) in teleost fish genomes, the miR-191 cluster locus is found within an intron of a protein-coding gene (intragenic) in the human genome. Interferon (IFN)-inducible and immune-related promoter elements found upstream of the teleost miR-462 cluster locus suggested roles in immune responses to viral pathogens in fish, while in humans, the miR-191 cluster functionally associated with cell cycle regulation. Stimulation of fish cell cultures with the IFN inducer poly I:C accordingly upregulated the expression of miR-462 and miR-731, while no stimulatory effect on miR-191 and miR-425 expression was observed in human cell lines. Despite high sequence conservation, evolution has thus resulted in different regulation and presumably also different functional roles of these orthologous miRNA clusters in different vertebrate lineages.

Functional conservation of both CDS- and 3'-UTR-located microRNA binding sites between species

MicroRNAs (miRNAs) mediate gene regulation posttranscriptionally through pairing of their seed (2-7 nt) to 3'-untranslated regions (3'-UTRs) or coding regions (coding sequences [CDSs]) of their target genes. CDS target sites generally show weaker repression effects than 3'-UTR sites. However, little is known about the conservation of the function, that is, repression effect, for these two groups of target sites. In addition, no systematic analysis of the evolutionary constraint on CDS sites exists to date. To address these questions, we performed RNA-sequencing to quantify the regulatory effect of miR-15a/miR-16 and miR-92a on their CDS and 3'-UTR targets in human and macaque cells. These miRs were knocked down transiently so the repression effect could be tracked immediately. Although on average CDS targets are less derepressed than 3'-UTR targets in both species, both the 3'-UTR targets and the CDS targets are functionally conserved. The evolutionary analysis of miRNA target sites shows that CDS sites are more conserved than nontarget control, albeit to a lesser extent than 3'-UTR sites. In conclusion, CDS target sites are functional, even though they are subject to less functional constraint than 3'-UTR target sites.

The role of RNA structure at 5' untranslated region in microRNA-mediated gene regulation

Recent studies have suggested that the secondary structure of the 5' untranslated region (5' UTR) of messenger RNA (mRNA) is important for microRNA (miRNA)-mediated gene regulation in humans. mRNAs that are targeted by miRNA tend to have a higher degree of local secondary structure in their 5' UTR; however, the general role of the 5' UTR in miRNA-mediated gene regulation remains unknown. We systematically surveyed the secondary structure of 5' UTRs in both plant and animal species and found a universal trend of increased mRNA stability near the 5' cap in mRNAs that are regulated by miRNA in animals, but not in plants. Intra-genome comparison showed that gene expression level, GC content of the 5' UTR, number of miRNA target sites, and 5' UTR length may influence mRNA structure near the 5' cap. Our results suggest that the 5' UTR secondary structure performs multiple functions in regulating post-transcriptional processes. Although the local structure immediately upstream of the start codon is involved in translation initiation, RNA structure near the 5' cap site, rather than the structure of the full-length 5' UTR sequences, plays an important role in miRNA-mediated gene regulation.

3'LIFE: a functional assay to detect miRNA targets in high-throughput

MicroRNAs (miRNAs) are short non-coding RNAs that regulate gene output at the post-transcriptional level by targeting degenerate elements primarily in 3'untranslated regions (3'UTRs) of mRNAs. Individual miRNAs can regulate networks of hundreds of genes, yet for the majority of miRNAs few, if any, targets are known. Misexpression of miRNAs is also a major contributor to cancer progression, thus there is a critical need to validate miRNA targets in high-throughput to understand miRNAs' contribution to tumorigenesis. Here we introduce a novel high-throughput assay to detect miRNA targets in 3'UTRs, called Luminescent Identification of Functional Elements in 3'UTRs (3'LIFE). We demonstrate the feasibility of 3'LIFE using a data set of 275 human 3'UTRs and two cancer-relevant miRNAs, let-7c and miR-10b, and compare our results to alternative methods to detect miRNA targets throughout the genome. We identify a large number of novel gene targets for these miRNAs, with only 32% of hits being bioinformatically predicted and 27% directed by non-canonical interactions. Functional analysis of target genes reveals consistent roles for each miRNA as either a tumor suppressor (let-7c) or oncogenic miRNA (miR-10b), and preferentially target multiple genes within regulatory networks, suggesting 3'LIFE is a rapid and sensitive method to detect miRNA targets in high-throughput.

A tri-component conservation strategy reveals highly confident microRNA-mRNA interactions and evolution of microRNA regulatory networks

MicroRNAs are small non-coding RNAs that can regulate expressions of their target genes at the post-transcriptional level. In this study, we propose a tri-component strategy that combines the conservation of microRNAs, homology of mRNA coding regions, and conserved microRNA binding sites in the 3′ untranslated regions to discover conserved microRNA-mRNA interactions. To validate the performance of our conservation strategy, we collected the experimentally validated microRNA-mRNA interactions from three databases as the golden standard. We found that the proposed strategy can improve the performance of existing target prediction algorithms by approximately 2–4 fold. In addition, we demonstrated that the proposed strategy could efficiently retain highly confident interactions from the intersection results of the existing algorithms and filter out the possible false positive predictions in the union one. Furthermore, this strategy can facilitate our ability to trace the homologues in different species that are targeted by the same miRNA family because it combines these three features to identify the conserved miRNA-mRNA interactions during evolution. Through an extensive application of the proposed conservation strategy to a study of the miR-1/206 regulatory network, we demonstrate that the target mRNA recruiting process could be associated with expansion of miRNA family during its evolution. We also uncovered the functional evolution of the miR-1/206 regulatory network. In this network, the early targeted genes tend to participate in more general and development-related functions. In summary, the conservation strategy is capable of helping to highlight the highly confident miRNA-mRNA interactions and can be further applied to reveal the evolutionary features of miRNA regulatory network and functions.

No gene in the genome makes sense except in the light of evolution

Evolutionary conservation has been an accurate predictor of functional elements across the first decade of metazoan genomics. More recently, there has been a move to define functional elements instead from biochemical annotations. Evolutionary methods are, however, more comprehensive than biochemical approaches can be and can assess quantitatively, especially for subtle effects, how biologically important--how injurious after mutation--different types of elements are. Evolutionary methods are thus critical for understanding the large fraction (up to 10%) of the human genome that does not encode proteins and yet might convey function. These methods can also capture the ephemeral nature of much noncoding functional sequence, with large numbers of functional elements having been gained and lost rapidly along each mammalian lineage. Here, we review how different strengths of purifying selection have impacted on protein-coding and non-protein-coding loci and on transcription factor binding sites in mammalian and fruit fly genomes.

Inferring the evolutionary history of primate microRNA binding sites: overcoming motif counting biases

The first microRNAs (miRNAs) were identified as essential, conserved regulators of gene expression, targeting the same genes across nearly all bilaterians. However, there are also prominent examples of conserved miRNAs whose functions appear to have shifted dramatically, sometimes over very brief periods of evolutionary time. To determine whether the functions of conserved miRNAs are stable or dynamic over evolutionary time scales, we have here defined the neutral turnover rates of short sequence motifs in predicted primate 3'-UTRs. We find that commonly used approaches to quantify motif turnover rates, which use a presence/absence scoring in extant lineages to infer ancestral states, are inherently biased to infer the accumulation of new motifs, leading to the false inference of continually increasing regulatory complexity over time. Using a maximum likelihood approach to reconstruct individual ancestral nucleotides, we observe that binding sites of conserved miRNAs in fact have roughly equal numbers of gain and loss events relative to ancestral states and turnover extremely slowly relative to nearly identical permutations of the same motif. Contrary to case studies showing examples of functional turnover, our systematic study of miRNA binding sites suggests that in primates, the regulatory roles of conserved miRNAs are strongly conserved. Our revised methodology may be used to quantify the mechanism by which regulatory networks evolve.