Large-scale analysis of microRNA evolution

The biogenesis and regulation of animal microRNAs

MicroRNAs (miRNAs) are small, yet profoundly influential, non-coding RNAs that base-pair with mRNAs to induce RNA silencing. Although the basic principles of miRNA biogenesis and function have been established, recent breakthroughs have yielded important new insights into the molecular mechanisms of miRNA biogenesis. In this Review, we discuss the metazoan miRNA biogenesis pathway step-by-step, focusing on the key biogenesis machinery, including the Drosha-DGCR8 complex (Microprocessor), exportin-5, Dicer and Argonaute. We also highlight newly identified cis-acting elements and their impact on miRNA maturation, informed by advanced high-throughput and structural studies, and discuss recently discovered mechanisms of clustered miRNA processing, target recognition and target-directed miRNA decay (TDMD). Lastly, we explore multiple regulatory layers of miRNA biogenesis, mediated by RNA-protein interactions, miRNA tailing (uridylation or adenylation) and RNA modifications.

Small and long non-coding RNAs: Past, present, and future

Since the introduction of the central dogma of molecular biology in 1958, various RNA species have been discovered. Messenger RNAs transmit genetic instructions from DNA to make proteins, a process facilitated by housekeeping non-coding RNAs (ncRNAs) such as small nuclear RNAs (snRNAs), ribosomal RNAs (rRNAs), and transfer RNAs (tRNAs). Over the past four decades, a wide array of regulatory ncRNAs have emerged as crucial players in gene regulation. In celebration of Cell's 50th anniversary, this Review explores our current understanding of the most extensively studied regulatory ncRNAs-small RNAs and long non-coding RNAs (lncRNAs)-which have profoundly shaped the field of RNA biology and beyond. While small RNA pathways have been well documented with clearly defined mechanisms, lncRNAs exhibit a greater diversity of mechanisms, many of which remain unknown. This Review covers pivotal events in their discovery, biogenesis pathways, evolutionary traits, action mechanisms, functions, and crosstalks among ncRNAs. We also highlight their roles in pathophysiological contexts and propose future research directions to decipher the unknowns of lncRNAs by leveraging lessons from small RNAs.

A time course analysis through diapause reveals dynamic temporal patterns of microRNAs associated with endocrine regulation in the butterfly Pieris napi

Organisms inhabiting highly seasonal environments must cope with a wide range of environmentally induced challenges. Many seasonal challenges require extensive physiological modification to survive. In winter, to survive extreme cold and limited resources, insects commonly enter diapause, which is an endogenously derived dormant state associated with minimized cellular processes and low energetic expenditure. Due to the high degree of complexity involved in diapause, substantial cellular regulation is required, of which our understanding primarily derives from the transcriptome via messenger RNA expression dynamics. Here we aim to advance our understanding of diapause by investigating microRNA (miRNA) expression in diapausing and direct developing pupae of the butterfly Pieris napi. We identified coordinated patterns of miRNA expression throughout diapause in both head and abdomen tissues of pupae, and via miRNA target identification, found several expression patterns to be enriched for relevant diapause-related physiological processes. We also identified two candidate miRNAs, miR-14-5p and miR-2a-3p, that are likely involved in diapause progression through their activity in the ecdysone pathway, a critical regulator of diapause termination. miR-14-5p targets phantom, a gene in the ecdysone synthesis pathway, and is upregulated early in diapause. miR-2a-3p has been found to be expressed in response to ecdysone, and is upregulated during diapause termination. Together, the expression patterns of these two miRNAs match our current understanding of the timing of hormonal regulation of diapause in P. napi and provide interesting candidates to further explore the mechanistic role of microRNAs in diapause regulation.

A comprehensive coding and microRNA transcriptome of vertebral bone in postlarvae and juveniles of Senegalese sole (Solea senegalensis)

Understanding vertebral bone development is essential to prevent skeletal malformations in farmed fish related to genetic and environmental factors. This is an important issue in Solea senegalensis, with special impact of spinal anomalies in postlarval and juvenile stages. Vertebral bone transcriptomics in farmed fish mainly comes from coding genes, and barely on miRNA expression. Here, we used RNA-seq of spinal samples to obtain the first comprehensive coding and miRNA transcriptomic repertoire for postlarval and juvenile vertebral bone, covering different vertebral phenotypes and egg-incubation temperatures related to skeleton health in S. senegalensis. Coding genes, miRNA and pathways regulating bone development and growth were identified. Differential transcriptomic profiles and suggestive mRNA-miRNA interactions were found between postlarvae and juveniles. Bone-related genes and functions were associated with the extracellular matrix, development and regulatory processes, calcium binding, retinol and lipid metabolism or response to stimulus, including those revealed by the miRNA targets related to signaling, cellular and metabolic processes, growth, cell proliferation and biological adhesion. Pathway enrichment associated with fish skeleton were identified when comparing postlarvae and juveniles: growth and bone development functions in postlarvae, while actin cytoskeleton, focal adhesion and proteasome related to bone remodeling in juveniles. The transcriptome data disclosed candidate coding and miRNA gene markers related to bone cell processes, references for functional studies of the anosteocytic bone of S. senegalensis. This study establishes a broad transcriptomic foundation to study healthy and anomalous spines under early thermal conditions across life-stages in S. senegalensis, and for comparative analysis of skeleton homeostasis and pathology in fish and vertebrates.

Time- and sex-dependent delayed effects of acute radiation exposure manifest via miRNA dysregulation

The detrimental effects of high-dose ionizing radiation on human health are well-known, but the influence of sex differences on the delayed effects of acute radiation exposure (DEARE) remains unclear. Here, we conducted six-month animal experiments using escalating radiation doses (7-9 Gy) on male and female C57BL/6 mice. The results show that female mice exhibited greater resistance to radiation, showing increased survival at six months post-total body irradiation. LD50/30 (lethal dose expected to cause 50% lethality in 30 days) for female mice is 8.08 Gy, while for male mice it is 7.76 Gy. DEARE causes time- and sex-dependent dysregulation of microRNA expression, processing enzymes, and the HOTAIR regulatory pathway. Differential regulation of molecular patterns associated with growth, development, apoptosis, and cancer is also observed in male and female mice. These findings shed light on the molecular basis of age and sex differences in DEARE response and emphasize the importance of personalized medicine for mitigating radiation-induced injuries and diseases.

Mitochondrial-related microRNAs and their roles in cellular senescence

Aging is a natural aspect of mammalian life. Although cellular mortality is inevitable, various diseases can hasten the aging process, resulting in abnormal or premature senescence. As cells age, they experience distinctive morphological and biochemical shifts, compromising their functions. Research has illuminated that cellular senescence coincides with significant alterations in the microRNA (miRNA) expression profile. Notably, a subset of aging-associated miRNAs, originally encoded by nuclear DNA, relocate to mitochondria, manifesting a mitochondria-specific presence. Additionally, mitochondria themselves house miRNAs encoded by mitochondrial DNA (mtDNA). These mitochondria-residing miRNAs, collectively referred to as mitochondrial miRNAs (mitomiRs), have been shown to influence mtDNA transcription and protein synthesis, thereby impacting mitochondrial functionality and cellular behavior. Recent studies suggest that mitomiRs serve as critical sensors for cellular senescence, exerting control over mitochondrial homeostasis and influencing metabolic reprogramming, redox equilibrium, apoptosis, mitophagy, and calcium homeostasis-all processes intimately connected to senescence. This review synthesizes current findings on mitomiRs, their mitochondrial targets, and functions, while also exploring their involvement in cellular aging. Our goal is to shed light on the potential molecular mechanisms by which mitomiRs contribute to the aging process.

ncOrtho: efficient and reliable identification of miRNA orthologs

MicroRNAs (miRNAs) are post-transcriptional regulators that finetune gene expression via translational repression or degradation of their target mRNAs. Despite their functional relevance, frameworks for the scalable and accurate detection of miRNA orthologs are missing. Consequently, there is still no comprehensive picture of how miRNAs and their associated regulatory networks have evolved. Here we present ncOrtho, a synteny informed pipeline for the targeted search of miRNA orthologs in unannotated genome sequences. ncOrtho matches miRNA annotations from multi-tissue transcriptomes in precision, while scaling to the analysis of hundreds of custom-selected species. The presence-absence pattern of orthologs to 266 human miRNA families across 402 vertebrate species reveals four bursts of miRNA acquisition, of which the most recent event occurred in the last common ancestor of higher primates. miRNA families are rarely modified or lost, but notable exceptions for both events exist. miRNA co-ortholog numbers faithfully indicate lineage-specific whole genome duplications, and miRNAs are powerful markers for phylogenomic analyses. Their exceptionally low genetic diversity makes them suitable to resolve clades where the phylogenetic signal is blurred by incomplete lineage sorting of ancestral alleles. In summary, ncOrtho allows to routinely consider miRNAs in evolutionary analyses that were thus far reserved to protein-coding genes.

Role of the proline-rich disordered domain of DROSHA in intronic microRNA processing

DROSHA serves as a gatekeeper of the microRNA (miRNA) pathway by processing primary transcripts (pri-miRNAs). While the functions of structured domains of DROSHA have been well documented, the contribution of N-terminal proline-rich disordered domain (PRD) remains elusive. Here we show that the PRD promotes the processing of miRNA hairpins located within introns. We identified a DROSHA isoform (p140) lacking the PRD, which is produced by proteolytic cleavage. Small RNA sequencing revealed that p140 is significantly impaired in the maturation of intronic miRNAs. Consistently, our minigene constructs demonstrated that PRD enhances the processing of intronic hairpins, but not those in exons. Splice site mutations did not affect the PRD's enhancing effect on intronic constructs, suggesting that the PRD acts independently of splicing reaction by interacting with sequences residing within introns. The N-terminal regions from zebrafish and Xenopus DROSHA can replace the human counterpart, indicating functional conservation despite poor sequence alignment. Moreover, we found that rapidly evolving intronic miRNAs are generally more dependent on PRD than conserved ones, suggesting a role of PRD in miRNA evolution. Our study reveals a new layer of miRNA regulation mediated by a low-complexity disordered domain that senses the genomic contexts of miRNA loci.

Critical Considerations for Investigating MicroRNAs during Tumorigenesis: A Case Study in Conceptual and Contextual Nuances of miR-211-5p in Melanoma

MicroRNAs are non-coding RNAs fundamental to metazoan development and disease. Although the aberrant regulation of microRNAs during mammalian tumorigenesis is well established, investigations into the contributions of individual microRNAs are wrought with conflicting observations. The underlying cause of these inconsistencies is often attributed to context-specific functions of microRNAs. We propose that consideration of both context-specific factors, as well as underappreciated fundamental concepts of microRNA biology, will permit a more harmonious interpretation of ostensibly diverging data. We discuss the theory that the biological function of microRNAs is to confer robustness to specific cell states. Through this lens, we then consider the role of miR-211-5p in melanoma progression. Using literature review and meta-analyses, we demonstrate how a deep understating of domain-specific contexts is critical for moving toward a concordant understanding of miR-211-5p and other microRNAs in cancer biology.

Clustering pattern and evolution characteristic of microRNAs in grass carp (Ctenopharyngodon idella)

BACKGROUND

A considerable fraction of microRNAs (miRNAs) are highly conserved, and certain miRNAs correspond to genomic clusters. The clustering of miRNAs can be advantageous, possibly by allowing coordinated expression. However, little is known about the evolutionary forces responsible for the loss and acquisition of miRNA and miRNA clusters.

RESULTS

The results demonstrated that several novel miRNAs arose throughout grass carp evolution. Duplication and de novo production were critical strategies for miRNA cluster formation. Duplicates accounted for a smaller fraction of the expansion in the grass carp miRNA than de novo creation. Clustered miRNAs are more conserved and change slower, whereas unique miRNAs usually have high evolution rates and low expression levels. The expression level of miRNA expression in clusters is strongly correlated.

CONCLUSIONS

This study examines the genomic distribution, evolutionary background, and expression regulation of grass carp miRNAs. Our findings provide novel insights into the genesis and development of miRNA clusters in teleost.

Evolution and Phylogeny of MicroRNAs - Protocols, Pitfalls, and Problems

MicroRNAs are important regulators in many eukaryotic lineages. Typical miRNAs have a length of about 22nt and are processed from precursors that form a characteristic hairpin structure. Once they appear in a genome, miRNAs are among the best-conserved elements in both animal and plant genomes. Functionally, they play an important role in particular in development. In contrast to protein-coding genes, miRNAs frequently emerge de novo. The genomes of animals and plants harbor hundreds of mutually unrelated families of homologous miRNAs that tend to be persistent throughout evolution. The evolution of their genomic miRNA complement closely correlates with important morphological innovation. In addition, miRNAs have been used as valuable characters in phylogenetic studies. An accurate and comprehensive annotation of miRNAs is required as a basis to understand their impact on phenotypic evolution. Since experimental data on miRNA expression are limited to relatively few species and are subject to unavoidable ascertainment biases, it is inevitable to complement miRNA sequencing by homology based annotation methods. This chapter reviews the state of the art workflows for homology based miRNA annotation, with an emphasis on their limitations and open problems.

Evolution after Whole-Genome Duplication: Teleost MicroRNAs

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

A prescient evolutionary model for genesis, duplication and differentiation of MIR160 homologs in Brassicaceae

MicroRNA160 is a class of nitrogen-starvation responsive genes which governs establishment of root system architecture by down-regulating AUXIN RESPONSE FACTOR genes (ARF10, ARF16 and ARF17) in plants. The high copy number of MIR160 variants discovered by us from land plants, especially polyploid crop Brassicas, posed questions regarding genesis, duplication, evolution and function. Absence of studies on impact of whole genome and segmental duplication on retention and evolution of MIR160 homologs in descendent plant lineages prompted us to undertake the current study. Herein, we describe ancestry and fate of MIR160 homologs in Brassicaceae in context of polyploidy driven genome re-organization, copy number and differentiation. Paralogy amongst Brassicaceae MIR160a, MIR160b and MIR160c was inferred using phylogenetic analysis of 468 MIR160 homologs from land plants. The evolutionarily distinct MIR160a was found to represent ancestral form and progenitor of MIR160b and MIR160c. Chronology of evolutionary events resulting in origin and diversification of genomic loci containing MIR160 homologs was delineated using derivatives of comparative synteny. A prescient model for causality of segmental duplications in establishment of paralogy in Brassicaceae MIR160, with whole genome duplication accentuating the copy number increase, is being posited in which post-segmental duplication events viz. differential gene fractionation, gene duplications and inversions are shown to drive divergence of chromosome segments. While mutations caused the diversification of MIR160a, MIR160b and MIR160c, duplicated segments containing these diversified genes suffered gene rearrangements via gene loss, duplications and inversions. Yet the topology of phylogenetic and phenetic trees were found congruent suggesting similar evolutionary trajectory. Over 80% of Brassicaceae genomes and subgenomes showed a preferential retention of single copy each of MIR160a, MIR160b and MIR160c suggesting functional relevance. Thus, our study provides a blue-print for reconstructing ancestry and phylogeny of MIRNA gene families at genomics level and analyzing the impact of polyploidy on organismal complexity. Such studies are critical for understanding the molecular basis of agronomic traits and deploying appropriate candidates for crop improvement.

Large-Scale Annotation and Evolution Analysis of MiRNA in Insects

Insects are among the most diverse and successful groups of animals and exhibit great morphological diversity and complexity. The innovation of wings and metamorphosis are some examples of the fascinating biological evolution of insects. Most microRNAs (miRNAs) contribute to canalization by conferring robustness to gene networks and thus increase the heritability of important phenotypes. Though previous studies have demonstrated how miRNAs regulate important phenotypes, little is still known about miRNA evolution in insects. Here, we used both small RNA-seq data and homology searching methods to annotate the miRNA repertoires of 152 arthropod species, including 135 insects and 17 noninsect arthropods. We identified 16,212 miRNA genes, and classified them into highly conserved (62), insect-conserved (90), and lineage-specific (354) miRNA families. The phylogenetic relationship of miRNA binary presence/absence dynamics implies that homoplastic loss of conserved miRNA families tends to occur in far-related morphologically simplified taxa, including scale insects (Coccoidea) and twisted-wing insects (Strepsiptera), leading to inconsistent phylogenetic tree reconstruction. The common ancestor of Insecta shares 62 conserved miRNA families, of which five were rapidly gained in the early winged-insects (Pterygota). We also detected extensive miRNA losses in Paraneoptera that are correlated with morphological reduction, and miRNA gains in early Endopterygota around the time holometabolous metamorphosis appeared. This was followed by abundant miRNA gains in Hymenoptera and Lepidoptera. In summary, we provide a comprehensive data set and a detailed evolutionary analysis of miRNAs in insects. These data will be important for future studies on miRNA functions associated with insect morphological innovation and trait biodiversity.

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.

Bioinformatics Analysis of Evolution and Human Disease Related Transposable Element-Derived microRNAs

Transposable element (TE) has the ability to insert into certain parts of the genome, and due to this event, it is possible for TEs to generate new factors and one of these factors are microRNAs (miRNA). miRNAs are non-coding RNAs made up of 19 to 24 nucleotides and numerous miRNAs are derived from TE. In this study, to support general knowledge on TE and miRNAs derived from TE, several bioinformatics tools and databases were used to analyze miRNAs derived from TE in two aspects: evolution and human disease. The distribution of TEs in diverse species presents that almost half of the genome is covered with TE in mammalians and less than a half in other vertebrates and invertebrates. Based on selected evolution-related miRNAs studies, a total of 51 miRNAs derived from TE were found and analyzed. For the human disease-related miRNAs, total of 34 miRNAs derived from TE were organized from the previous studies. In summary, abundant miRNAs derived from TE are found, however, the function of miRNAs derived from TE is not informed either. Therefore, this study provides theoretical understanding of miRNAs derived from TE by using various bioinformatics tools.

MicroRNAs: From Mechanism to Organism

MicroRNAs (miRNAs) are short, regulatory RNAs that act as post-transcriptional repressors of gene expression in diverse biological contexts. The emergence of small RNA-mediated gene silencing preceded the onset of multicellularity and was followed by a drastic expansion of the miRNA repertoire in conjunction with the evolution of complexity in the plant and animal kingdoms. Along this process, miRNAs became an essential feature of animal development, as no higher metazoan lineage tolerated loss of miRNAs or their associated protein machinery. In fact, ablation of the miRNA biogenesis machinery or the effector silencing factors results in severe embryogenesis defects in every animal studied. In this review, we summarize recent mechanistic insight into miRNA biogenesis and function, while emphasizing features that have enabled multicellular organisms to harness the potential of this broad class of repressors. We first discuss how different mechanisms of regulation of miRNA biogenesis are used, not only to generate spatio-temporal specificity of miRNA production within an animal, but also to achieve the necessary levels and dynamics of expression. We then explore how evolution of the mechanism for small RNA-mediated repression resulted in a diversity of silencing complexes that cause different molecular effects on their targets. Multicellular organisms have taken advantage of this variability in the outcome of miRNA-mediated repression, with differential use in particular cell types or even distinct subcellular compartments. Finally, we present an overview of how the animal miRNA repertoire has evolved and diversified, emphasizing the emergence of miRNA families and the biological implications of miRNA sequence diversification. Overall, focusing on selected animal models and through the lens of evolution, we highlight canonical mechanisms in miRNA biology and their variations, providing updated insight that will ultimately help us understand the contribution of miRNAs to the development and physiology of multicellular organisms.

Small RNAs in eucaryotes: new clues for amplifying microRNA benefits

miRNAs, the smallest nucleotide molecules able to regulate gene expression at post transcriptional level, are found in both animals and plants being involved in fundamental processes for growth and development of living organisms. The number of miRNAs has been hypothesized to increase when some organisms specialized the process of mastication and grinding of food. Further to the vertical transmission, miRNAs can undergo horizontal transmission among different species, in particular between plants and animals. In the last years, an increasing number of studies reported that miRNA passage occurs through feeding, and that in animals, plant miRNAs can survive the gastro intestinal digestion and transferred by blood into host cells, where they can exert their functions modulating gene expression. The present review reports studies on miRNAs during evolution, with particular focus on biogenesis and mechanisms regulating their stability in plants and animals. The different biogenesis and post biogenesis modifications allow to discriminate miRNAs of plant origin from those of animal origin, and make it possible to better clarify the controversial question on whether a possible cross-kingdom miRNA transfer through food does exist. The majority of human medicines and supplements derive from plants and a regular consumption of plant food is suggested for their beneficial effects in the prevention of metabolic diseases, cancers, and dietary related disorders. So far, these beneficial effects have been generally attributed to the content of secondary metabolites, whereas mechanisms regarding other components remain unclear. Therefore, in light of the above reported studies miRNAs could result another component for the medical properties of plants. miRNAs have been mainly studied in mammals characterizing their sequences and molecular targets as available in public databases. The herein presented studies provide evidences that miRNA situation is much more complex than the static situation reported in databases. Indeed, miRNAs may have redundant activities, variable sequences, different methods of biogenesis, and may be differently influenced by external and environmental factors. In-depth knowledge of mechanisms of synthesis, regulation and transfer of plant miRNAs to other species can open new frontiers in the therapy of many human diseases, including cancer.

Conservation and novelty in the microRNA genomic landscape of hyperdiverse cichlid fishes

MicroRNAs (miRNAs) play crucial roles in the post-transcriptional control of messenger RNA (mRNA). These miRNA-mRNA regulatory networks are present in nearly all organisms and contribute to development, phenotypic divergence, and speciation. To examine the miRNA landscape of cichlid fishes, one of the most species-rich families of vertebrates, we profiled the expression of both miRNA and mRNA in a diverse set of cichlid lineages. Among these, we found that conserved miRNAs differ from recently arisen miRNAs (i.e. lineage specific) in average expression levels, number of target sites, sequence variability, and physical clustering patterns in the genome. Furthermore, conserved miRNA target sites tend to be enriched at the 5' end of protein-coding gene 3' UTRs. Consistent with the presumed regulatory role of miRNAs, we detected more negative correlations between the expression of miRNA-mRNA functional pairs than in random pairings. Finally, we provide evidence that novel miRNA targets sites are enriched in genes involved in protein synthesis pathways. Our results show how conserved and evolutionarily novel miRNAs differ in their contribution to the genomic landscape and highlight their particular evolutionary roles in the adaptive diversification of cichlids.

MicroRNA Gene Regulation in Extremely Young and Parallel Adaptive Radiations of Crater Lake Cichlid Fish

Cichlid fishes provide textbook examples of explosive phenotypic diversification and sympatric speciation, thereby making them ideal systems for studying the molecular mechanisms underlying rapid lineage divergence. Despite the fact that gene regulation provides a critical link between diversification in gene function and speciation, many genomic regulatory mechanisms such as microRNAs (miRNAs) have received little attention in these rapidly diversifying groups. Therefore, we investigated the posttranscriptional regulatory role of miRNAs in the repeated sympatric divergence of Midas cichlids (Amphilophus spp.) from Nicaraguan crater lakes. Using miRNA and mRNA sequencing of embryos from five Midas species, we first identified miRNA binding sites in mRNAs and highlighted the presences of a surprising number of novel miRNAs in these adaptively radiating species. Then, through analyses of expression levels, we identified putative miRNA/gene target pairs with negatively correlated expression level that were consistent with the role of miRNA in downregulating mRNA. Furthermore, we determined that several miRNA/gene pairs show convergent expression patterns associated with the repeated benthic/limnetic sympatric species divergence implicating these miRNAs as potential molecular mechanisms underlying replicated sympatric divergence. Finally, as these candidate miRNA/gene pairs may play a central role in phenotypic diversification in these cichlids, we characterized the expression domains of selected miRNAs and their target genes via in situ hybridization, providing further evidence that miRNA regulation likely plays a role in the Midas cichlid adaptive radiation. These results provide support for the hypothesis that extremely quickly evolving miRNA regulation can contribute to rapid evolutionary divergence even in the presence of gene flow.

Import of Non-Coding RNAs into Human Mitochondria: A Critical Review and Emerging Approaches

Mitochondria harbor their own genetic system, yet critically depend on the import of a number of nuclear-encoded macromolecules to ensure their expression. In all eukaryotes, selected non-coding RNAs produced from the nuclear genome are partially redirected into the mitochondria, where they participate in gene expression. Therefore, the mitochondrial RNome represents an intricate mixture of the intrinsic transcriptome and the extrinsic RNA importome. In this review, we summarize and critically analyze data on the nuclear-encoded transcripts detected in human mitochondria and outline the proposed molecular mechanisms of their mitochondrial import. Special attention is given to the various experimental approaches used to study the mitochondrial RNome, including some recently developed genome-wide and in situ techniques.

Well-Annotated microRNAomes Do Not Evidence Pervasive miRNA Loss

microRNAs are conserved noncoding regulatory factors implicated in diverse physiological and developmental processes in multicellular organisms, as causal macroevolutionary agents and for phylogeny inference. However, the conservation and phylogenetic utility of microRNAs has been questioned on evidence of pervasive loss. Here, we show that apparent widespread losses are, largely, an artefact of poorly sampled and annotated microRNAomes. Using a curated data set of animal microRNAomes, we reject the view that miRNA families are never lost, but they are rarely lost (92% are never lost). A small number of families account for a majority of losses (1.7% of families account for >45% losses), and losses are associated with lineages exhibiting phenotypic simplification. Phylogenetic analyses based on the presence/absence of microRNA families among animal lineages, and based on microRNA sequences among Osteichthyes, demonstrate the power of these small data sets in phylogenetic inference. Perceptions of widespread evolutionary loss of microRNA families are due to the uncritical use of public archives corrupted by spurious microRNA annotations, and failure to discriminate false absences that occur because of incomplete microRNAome annotation.

In silico identification of microRNAs predicted to regulate N-myristoyltransferase and Methionine Aminopeptidase 2 functions in cancer and infectious diseases

Protein myristoylation is a key protein modification carried out by N-Myristoyltransferase (NMT) after Methionine aminopeptidase 2 (MetAP2) removes methionine from the amino-terminus of the target protein. Protein myristoylation by NMT augments several signaling pathways involved in a myriad of cellular processes, including developmental pathways and pathways that when dysregulated lead to cancer or immune dysfunction. The emerging evidence pointing to NMT-mediated myristoylation as a major cellular regulator underscores the importance of understanding the framework of this type of signaling event. Various studies have investigated the role that myristoylation plays in signaling dysfunction by examining differential gene or protein expression between normal and diseased states, such as cancers or following HIV-1 infection, however no study exists that addresses the role of microRNAs (miRNAs) in the regulation of myristoylation. By performing a large scale bioinformatics and functional analysis of the miRNAs that target key genes involved in myristoylation (NMT1, NMT2, MetAP2), we have narrowed down a list of promising candidates for further analysis. Our condensed panel of miRNAs identifies 35 miRNAs linked to cancer, 21 miRNAs linked to developmental and immune signaling pathways, and 14 miRNAs linked to infectious disease (primarily HIV). The miRNAs panel that was analyzed revealed several NMT-targeting mRNAs (messenger RNA) that are implicated in diseases associated with NMT signaling alteration, providing a link between the realms of miRNA and myristoylation signaling. These findings verify miRNA as an additional facet of myristoylation signaling that must be considered to gain a full perspective. This study provides the groundwork for future studies concerning NMT-transcript-binding miRNAs, and will potentially lead to the development of new diagnostic/prognostic biomarkers and therapeutic targets for several important diseases.

Fishing Into the MicroRNA Transcriptome

In the last decade, several studies have been focused on revealing the microRNA (miRNA) repertoire and determining their functions in farm animals such as poultry, pigs, cattle, and fish. These small non-protein coding RNA molecules (18-25 nucleotides) are capable of controlling gene expression by binding to messenger RNA (mRNA) targets, thus interfering in the final protein output. MiRNAs have been recognized as the main regulators of biological features of economic interest, including body growth, muscle development, fat deposition, and immunology, among other highly valuable traits, in aquatic livestock. Currently, the miRNA repertoire of some farmed fish species has been identified and characterized, bringing insights about miRNA functions, and novel perspectives for improving health and productivity. In this review, we summarize the current advances in miRNA research by examining available data on Neotropical and other key species exploited by fisheries and in aquaculture worldwide and discuss how future studies on Neotropical fish could benefit from this knowledge. We also make a horizontal comparison of major results and discuss forefront strategies for miRNA manipulation in aquaculture focusing on forward-looking ideas for forthcoming research.

Biological function in the twilight zone of sequence conservation

Strong DNA conservation among divergent species is an indicator of enduring functionality. With weaker sequence conservation we enter a vast ‘twilight zone’ in which sequence subject to transient or lower constraint cannot be distinguished easily from neutrally evolving, non-functional sequence. Twilight zone functional sequence is illuminated instead by principles of selective constraint and positive selection using genomic data acquired from within a species’ population. Application of these principles reveals that despite being biochemically active, most twilight zone sequence is not functional.

MicroRNAs as biomarkers for psychiatric disorders with a focus on autism spectrum disorder: Current progress in genetic association studies, expression profiling, and translational research

MicroRNAs (miRNAs) are a group of small noncoding RNA molecules, 18-25 nucleotides in length, which can negatively regulate gene expression at the post-transcriptional level by binding to messenger RNAs. About half of all identified miRNAs in humans are expressed in the brain and display regulatory functions important for many biological processes related to the development of the central nervous system (CNS). Disruptions in miRNA biogenesis and miRNA-target interaction have been related to CNS diseases, including psychiatric disorders. In this review, we focus on the role of miRNAs in autism spectrum disorder (ASD) and summarize recent findings about ASD-associated genetic variants in miRNA genes, in miRNA biogenesis genes, and miRNA targets. We discuss deregulation of miRNA expression in ASD and functional validation of ASD-related miRNAs in animal models. Including miRNAs in studies of ASD will contribute to our understanding of its etiology and pathogenesis and facilitate the discrimination between different disease subgroups. Autism Res 2017. © 2017 International Society for Autism Research, Wiley Periodicals, Inc. Autism Res 2017, 10: 1184-1203. © 2017 International Society for Autism Research, Wiley Periodicals, Inc.

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

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

Pervasive microRNA Duplication in Chelicerates: Insights from the Embryonic microRNA Repertoire of the Spider Parasteatoda tepidariorum

MicroRNAs are small (∼22 nt) noncoding RNAs that repress translation and therefore regulate the production of proteins from specific target mRNAs. microRNAs have been found to function in diverse aspects of gene regulation within animal development and many other processes. Among invertebrates, both conserved and novel, lineage specific, microRNAs have been extensively studied predominantly in holometabolous insects such as Drosophila melanogaster However little is known about microRNA repertoires in other arthropod lineages such as the chelicerates. To understand the evolution of microRNAs in this poorly sampled subphylum, we characterized the microRNA repertoire expressed during embryogenesis of the common house spider Parasteatoda tepidariorum We identified a total of 148 microRNAs in P. tepidariorum representing 66 families. Approximately half of these microRNA families are conserved in other metazoans, while the remainder are specific to this spider. Of the 35 conserved microRNAs families 15 had at least two copies in the P. tepidariorum genome. A BLAST-based approach revealed a similar pattern of duplication in other spiders and a scorpion, but not among other chelicerates and arthropods, with the exception of a horseshoe crab. Among the duplicated microRNAs we found examples of lineage-specific tandem duplications, and the duplication of entire microRNA clusters in three spiders, a scorpion, and in a horseshoe crab. Furthermore, we found that paralogs of many P. tepidariorum microRNA families exhibit arm switching, which suggests that duplication was often followed by sub- or neofunctionalization. Our work shows that understanding the evolution of microRNAs in the chelicerates has great potential to provide insights into the process of microRNA duplication and divergence and the evolution of animal development.

Synteny and comparative analysis of miRNA retention, conservation, and structure across Brassicaceae reveals lineage- and sub-genome-specific changes

The recent availability of genome sequences together with syntenic block information for Brassicaceae offers an opportunity to study microRNA (miRNA) evolution across this family. We employed a synteny-based comparative genomics strategy to unambiguously identify miRNA homologs from the genome sequence of members of Brassicaceae. Such an analysis of miRNA across Brassicaceae allowed us to classify miRNAs as conserved, lineage-, karyotype- and sub-genome-specific. The differential loss of miRNA from sub-genomes in polyploid genomes of Brassica rapa and Brassica oleracea shows that miRNA also follows the rules of gene fractionation as observed in the case of protein-coding genes. The study of mature and miR* region of precursors revealed instances of in-dels and SNPs which reflect the evolutionary history of the genomes. High level of conservation in miR* regions in some cases points to their functional relevance which needs to be further investigated. We further show that sequence and length variability in precursor sequences can affect the free energy and foldback structure of miRNA which may ultimately affect their biogenesis and expression in the biological system.

A Uniform System for the Annotation of Vertebrate microRNA Genes and the Evolution of the Human microRNAome

Although microRNAs (miRNAs) are among the most intensively studied molecules of the past 20 years, determining what is and what is not a miRNA has not been straightforward. Here, we present a uniform system for the annotation and nomenclature of miRNA genes. We show that less than a third of the 1,881 human miRBase entries, and only approximately 16% of the 7,095 metazoan miRBase entries, are robustly supported as miRNA genes. Furthermore, we show that the human repertoire of miRNAs has been shaped by periods of intense miRNA innovation and that mature gene products show a very different tempo and mode of sequence evolution than star products. We establish a new open access database--MirGeneDB ( http://mirgenedb.org )--to catalog this set of miRNAs, which complements the efforts of miRBase but differs from it by annotating the mature versus star products and by imposing an evolutionary hierarchy upon this curated and consistently named repertoire.

Evidence for the expression of abundant microRNAs in the locust genome

Substantial accumulation of neutral sequences accounts for genome size expansion in animal genomes. Numerous novel microRNAs (miRNAs), which evolve in a birth and death manner, are considered evolutionary neutral sequences. The migratory locust is an ideal model to determine whether large genomes contain abundant neutral miRNAs because of its large genome size. A total of 833 miRNAs were discovered, and several miRNAs were randomly chosen for validation by Northern blot and RIP-qPCR. Three additional verification methods, namely, processing-dependent methods of miRNA biogenesis using RNAi, evolutionary comparison with closely related species, and evidence supported by tissue-specific expression, were applied to provide compelling results that support the authenticity of locust miRNAs. We observed that abundant local duplication events of miRNAs, which were unique in locusts compared with those in other insects with small genome sizes, may be responsible for the substantial acquisition of miRNAs in locusts. Together, multiple evidence showed that the locust genome experienced a burst of miRNA acquisition, suggesting that genome size expansion may have considerable influences of miRNA innovation. These results provide new insight into the genomic dynamics of miRNA repertoires under genome size evolution.

The mitochondrial lncRNA ASncmtRNA-2 is induced in aging and replicative senescence in Endothelial Cells

Age-associated cardiovascular diseases are at least partially ascribable to vascular cell senescence. Replicative senescence (RS) and stress-induced premature senescence (SIPS) are provoked respectively by endogenous (telomere erosion) and exogenous (H2O2, UV) stimuli resulting in cell cycle arrest in G1 and G2 phases. In both scenarios, mitochondria-derived ROS are important players in senescence initiation. We aimed to define whether a mtDNA-transcribed long-non-coding-RNA (lncRNA), ASncmtRNA-2, has a role in vascular aging and senescence. Aortas of old mice, characterized by increased senescence, showed an increment in ASncmtRNA-2 expression. In vitro analysis of Endothelial Cells (EC) and Vascular Smooth Muscle Cells (VSMC) established that ASncmtRNA-2 is induced in EC, but not in VSMC, during RS. Surprisingly, ASncmtRNA-2 is not upregulated in two different EC SIPS scenarios, treated with H2O2 and UV. The p16 gene displayed similar ASncmtRNA-2 expression patterns, suggesting a possible co-regulation of the two genes. Interestingly, the expression of two miRNAs, hsa-miR-4485 and hsa-miR-1973, with perfect homology to the double strand region of ASncmtRNA-2 and originating at least in part from a mitochondrial transcript, was induced in RS, opening to the possibility that this lncRNA functions as a non-canonical precursor of these miRNAs. Cell cycle analysis of EC transiently over-expressing ASncmtRNA-2 revealed an accumulation of EC in the G2/M phase, but not in the G1 phase. We propose that ASncmtRNA-2 in EC might be involved in the RS establishment by participating in the cell cycle arrest in G2/M phase, possibly through the production of hsa-miR-4485 and hsa-miR-1973. This article is part of a Special Issue entitled: Mitochondria.

A critical appraisal of the use of microRNA data in phylogenetics

Recent progress in resolving the tree of life continues to expose relationships that resist resolution, which drives the search for novel sources of information to solve these difficult phylogenetic problems. A recent example, the presence and absence of microRNA families, has been vigorously promoted as an ideal source of phylogenetic data and has been applied to several perennial phylogenetic problems. The utility of such data for phylogenetic inference hinges critically both on developing stochastic models that provide a reasonable description of the process that give rise to these data, and also on the careful validation of those models in real inference scenarios. Remarkably, however, the statistical behavior and phylogenetic utility of microRNA data have not yet been rigorously characterized. Here we explore the behavior and performance of microRNA presence/absence data under a variety of evolutionary models and reexamine datasets from several previous studies. We find that highly heterogeneous rates of microRNA gain and loss, pervasive secondary loss, and sampling error collectively render microRNA-based inference of phylogeny difficult. Moreover, our reanalyses fundamentally alter the conclusions for four of the five studies that we reexamined. Our results indicate that the capacity of miRNA data to resolve the tree of life has been overstated, and we urge caution in their application and interpretation.

MicroRNA-target binding structures mimic microRNA duplex structures in humans

Traditionally, researchers match a microRNA guide strand to mRNA sequences using sequence comparisons to predict its potential target genes. However, many of the predictions can be false positives due to limitations in sequence comparison alone. In this work, we consider the association of two related RNA structures that share a common guide strand: the microRNA duplex and the microRNA-target binding structure. We have analyzed thousands of such structure pairs and found many of them share high structural similarity. Therefore, we conclude that when predicting microRNA target genes, considering just the microRNA guide strand matches to gene sequences may not be sufficient--the microRNA duplex structure formed by the guide strand and its companion passenger strand must also be considered. We have developed software to translate RNA binding structure into encoded representations, and we have also created novel automatic comparison methods utilizing such encoded representations to determine RNA structure similarity. Our software and methods can be utilized in the other RNA secondary structure comparisons as well.

Construction of microRNA functional families by a mixture model of position weight matrices

MicroRNAs (miRNAs) are small regulatory molecules that repress the translational processes of their target genes by binding to their 3′ untranslated regions (3′ UTRs). Because the target genes are predominantly determined by their sequence complementarity to the miRNA seed regions (nucleotides 2–7) which are evolutionarily conserved, it is inferred that the target relationships and functions of the miRNA family members are conserved across many species. Therefore, detecting the relevant miRNA families with confidence would help to clarify the conserved miRNA functions, and elucidate miRNA-mediated biological processes. We present a mixture model of position weight matrices for constructing miRNA functional families. This model systematically finds not only evolutionarily conserved miRNA family members but also functionally related miRNAs, as it simultaneously generates position weight matrices representing the conserved sequences. Using mammalian miRNA sequences, in our experiments, we identified potential miRNA groups characterized by similar sequence patterns that have common functions. We validated our results using score measures and by the analysis of the conserved targets. Our method would provide a way to comprehensively identify conserved miRNA functions.

miRNAs: small genes with big potential in metazoan phylogenetics

microRNAs (miRNAs) are a key component of gene regulatory networks and have been implicated in the regulation of virtually every biological process found in multicellular eukaryotes. What makes them interesting from a phylogenetic perspective is the high conservation of primary sequence between taxa, their accrual in metazoan genomes through evolutionary time, and the rarity of secondary loss in most metazoan taxa. Despite these properties, the use of miRNAs as phylogenetic markers has not yet been discussed within a clear conceptual framework. Here we highlight five properties of miRNAs that underlie their utility in phylogenetics: 1) The processes of miRNA biogenesis enable the identification of novel miRNAs without prior knowledge of sequence; 2) The continuous addition of miRNA families to metazoan genomes through evolutionary time; 3) The low level of secondary gene loss in most metazoan taxa; 4) The low substitution rate in the mature miRNA sequence; and 5) The small probability of convergent evolution of two miRNAs. Phylogenetic analyses using both Bayesian and parsimony methods on a eumetazoan miRNA data set highlight the potential of miRNAs to become an invaluable new tool, especially when used as an additional line of evidence, to resolve previously intractable nodes within the tree of life.

Detecting and comparing non-coding RNAs in the high-throughput era

In recent years there has been a growing interest in the field of non-coding RNA. This surge is a direct consequence of the discovery of a huge number of new non-coding genes and of the finding that many of these transcripts are involved in key cellular functions. In this context, accurately detecting and comparing RNA sequences has become important. Aligning nucleotide sequences is a key requisite when searching for homologous genes. Accurate alignments reveal evolutionary relationships, conserved regions and more generally any biologically relevant pattern. Comparing RNA molecules is, however, a challenging task. The nucleotide alphabet is simpler and therefore less informative than that of amino-acids. Moreover for many non-coding RNAs, evolution is likely to be mostly constrained at the structural level and not at the sequence level. This results in very poor sequence conservation impeding comparison of these molecules. These difficulties define a context where new methods are urgently needed in order to exploit experimental results to their full potential. This review focuses on the comparative genomics of non-coding RNAs in the context of new sequencing technologies and especially dealing with two extremely important and timely research aspects: the development of new methods to align RNAs and the analysis of high-throughput data.