Evolutionary conservation of microRNA regulatory circuits: an examination of microRNA gene complexity and conserved microRNA-target interactions through metazoan phylogeny

During the last decade, a variety of critical biological processes, including early embryo development, cell proliferation, differentiation, apoptosis, and metabolic regularity, have been shown to be genetically regulated by a large gene family encoding a class of tiny RNA molecules termed microRNAs (miRNAs). All miRNAs share a common biosynthetic pathway and reaction mechanisms. The sequence of many miRNAs is found to be conserved, in their mature form, among different organisms. In addition, the evolutionary appearance of multicellular organisms appears to correlate with the appearance of the miRNA pathway for regulating gene expression. The miRNA pathway has the potential to regulate vast networks of gene products in a coordinate manner. Recent evidence has not only implicated the miRNA pathway in regulating a vast array of basic cellular processes but also specialized processes that are required for cellular identity and tissue specificity. A survey of the literature shows that some miRNA pathways are conserved virtually intact throughout phylogeny while miRNA diversity also correlates with speciation. The number of miRNA genes, the expression of miRNAs, and target diversities of miRNAs tend to be positively correlated with morphological complexities observed in animals. Thus, organismal complexity can be estimated by the complexity of the miRNA circuitry. The complexity of the miRNA gene families establishes a link between genotypic complexity and phenotypic complexity in animal evolution. In this paper, we start with the discussion of miRNA conservation. Then we interpret the trends in miRNA conservation to deduce miRNA evolutionary trends in metazoans. Based on these conservation patterns observed in each component of the miRNA regulatory system, we attempt to propose a global insight on the probable consistency between morphological evolution in animals and the molecular evolution of miRNA gene activity in the cell.     

High-Throughput miRNA and mRNA Sequencing of Paired Colorectal Normal,Tumor and Metastasis Tissues and Bioinformatic Modeling of miRNA-1 Therapeutic Applications

MiRNAs are discussed as diagnostic and therapeutic molecules. However, effective miRNA drug treatments with miRNAs are, so far, hampered by the complexity of the miRNA networks. To identify potential miRNA drugs in colorectal cancer, we profiled miRNA and mRNA expression in matching normal, tumor and metastasis tissues of eight patients by Illumina sequencing. We validated six miRNAs in a large tissue screen containing 16 additional tumor entities and identified miRNA-1, miRNA-129, miRNA-497 and miRNA-215 as constantly de-regulated within the majority of cancers. Of these, we investigated miRNA-1 as representative in a systems-biology simulation of cellular cancer models implemented in PyBioS and assessed the effects of depletion as well as overexpression in terms of miRNA-1 as a potential treatment option. In this system, miRNA-1 treatment reverted the disease phenotype with different effectiveness among the patients. Scoring the gene expression changes obtained through mRNA-Seq from the same patients we show that the combination of deep sequencing and systems biological modeling can help to identify patient-specific responses to miRNA treatments. We present this data as guideline for future pre-clinical assessments of new and personalized therapeutic options.     

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

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

Evolution of microRNA genes in Oryza sativa and Arabidopsis thaliana: an update of the inverted duplication model

The origin and evolution of microRNA (miRNA) genes, which are of significance in tuning and buffering gene expressions in a number of critical cellular processes, have long attracted evolutionary biologists. However, genome-wide perspectives on their origins, potential mechanisms of their de novo generation and subsequent evolution remain largely unsolved in flowering plants. Here, genome-wide analyses of Oryza sativa and Arabidopsis thaliana revealed apparently divergent patterns of miRNA gene origins. A large proportion of miRNA genes in O. sativa were TE-related and MITE-related miRNAs in particular, whereas the fraction of these miRNA genes much decreased in A. thaliana. Our results show that the majority of TE-related and pseudogene-related miRNA genes have originated through inverted duplication instead of segmental or tandem duplication events. Based on the presented findings, we hypothesize and illustrate the four likely molecular mechanisms to de novo generate novel miRNA genes from TEs and pseudogenes. Our rice genome analysis demonstrates that non-MITEs and MITEs mediated inverted duplications have played different roles in de novo generating miRNA genes. It is confirmed that the previously proposed inverted duplication model may give explanations for non-MITEs mediated duplication events. However, many other miRNA genes, known from the earlier proposed model, were rather arisen from MITE transpositions into target genes to yield binding sites. We further investigated evolutionary processes spawned from de novo generated to maturely-formed miRNA genes and their regulatory systems. We found that miRNAs increase the tunability of some gene regulatory systems with low gene copy numbers. The results also suggest that gene balance effects may have largely contributed to the evolution of miRNA regulatory systems.     

MicroRNA对多细胞动物复杂性进化的影响

MicroRNA(miRNA)是一种长度约为22个碱基的非编码单链小分子RNA。作为一类重要的转录后基因表达调控因子,miRNA参与了广泛的生物学过程,如发育时程调控、细胞分化、凋亡、肿瘤以及病毒抵抗等。然而,除了在个体发生过程中的重要功能外,越来越多的研究表明,miRNA在系统发生中也扮演着关键的角色。基因表达模式的不同被广泛地认为是物种内和物种间表型差异的根源,动物物种间miRNA的保守性和多样性研究提示miRNA对物种间表型差异以及动物进化起着重要的作用。文章介绍了miRNA产生过程和作用机制,重点探讨了miRNA在动物进化过程中的作用,从miRNA的进化速度、miRNA表达的时空特异性、miRNA作用靶位点变异以及miRNA基因的扩增与丢失4个方面论述miRNA介导的基因调控网络对多细胞动物发育复杂性进化的影响,推测miRNA在多细胞动物进化过程中驱动了复杂性的增加。     

miR-155-3p: processing by-product or rising star in immunity and cancer?

MicroRNAs (miRNAs) are key players in gene regulation that target specific mRNAs for degradation or translational repression. Each miRNA is synthesized as a miRNA duplex comprising two strands (5p and 3p). However, only one of the two strands becomes active and is selectively incorporated into the RNA-induced silencing complex in a process known as miRNA strand selection. Recently, significant progress has been made in understanding the factors and processes involved in strand selection. Here, we explore the selection and functionality of the miRNA star strand (either 5p or 3p), which is generally present in the cell at low levels compared to its partner strand and, historically, has been thought to possess no biological activity. We also highlight the concepts of miRNA arm switching and miRNA isomerism. Finally, we offer insights into the impact of aberrant strand selection on immunity and cancer. Leading us through this journey is miR-155, a well-established regulator of immunity and cancer, and the increasing evidence that its 3p strand plays a role in these arenas. Interestingly, the miR-155-5p/-3p ratio appears to vary dependent on the timing of the immune response, and the 3p strand seems to play a regulatory role upon its partner 5p strand.     

Selection on Synonymous Sites for Increased Accessibility around miRNA Binding Sites in Plants

Synonymous codons are widely selected for various biological mechanisms in both prokaryotes and eukaryotes. Recent evidence suggests that microRNA (miRNA) function may affect synonymous codon choices near miRNA target sites. To better understand this, we perform genome-wide analysis on synonymous codon usage around miRNA target sites in four plant genomes. We observed a general trend of increased site accessibility around miRNA target sites in plants. Guanine-cytosine (GC)-poor codons are preferred in the flank region of miRNA target sites. Within-genome analyses show significant variation among miRNA targets in species. GC content of the target gene can partly explain the variation of site accessibility among miRNA targets. miRNA targets in GC-rich genes show stronger selection signals than those in GC-poor genes. Gene's codon usage bias and the conservation level of miRNA and its target also have some effects on site accessibility, but the expression level of miRNA or its target and the mechanism of miRNA activity do not contribute to site accessibility differences among miRNA targets. We suggest that synonymous codons near miRNA targets are selected for efficient miRNA binding and proper miRNA function. Our results present a new dimension of natural selection on synonymous codons near miRNA target sites in plants, which will have important implications of coding sequence evolution.     

The role of pleiotropy vs signaller-receiver gene epistasis in life history trade-offs: dissecting the genomic architecture of organismal design in social systems

Traditional life history theory ignores trade-offs due to social interactions, yet social systems expand the set of possible trade-offs affecting a species evolution--by introducing asymmetric interactions between the sexes, age classes and invasion of alternative strategies. We outline principles for understanding gene epistasis due to signaller-receiver dynamics, gene interactions between individuals, and impacts on life history trade-offs. Signaller-receiver epistases create trade-offs among multiple correlated traits that affect fitness, and generate multiple fitness optima conditional on frequency of alternative strategies. In such cases, fitness epistasis generated by selection can maintain linkage disequilibrium, even among physically unlinked loci. In reviewing genetic methods for studying life history trade-offs, we conclude that current artificial selection or gene manipulation experiments focus on pleiotropy. Multi-trait selection experiments, multi-gene engineering methods or multiple endocrine manipulations can test for epistasis and circumvent these limitations. In nature, gene mapping in field pedigrees is required to study social gene epistases and associated trade-offs. Moreover, analyses of correlational selection and frequency-dependent selection are necessary to study epistatic social system trade-offs, which can be achieved with group-structured versions of Price's (1970) equation.     

Genome-wide screen for aberrantly expressed miRNAs reveals miRNA profile signature in breast cancer

Dysregulation in the expression of miRNAs contributes to the occurrence and development of many human cancers. We herein attempted to obtain the potential association between miRNA expression profile and breast cancer by applying high-throughput sequencing technology. Small RNAs from seven paired tumor and adjacent normal tissue samples were sequenced. To determine the miRNA expression profiles in tissues and sera, another five equally pooled serum samples from 20 patients and 30 normal women were sequenced. Despite a similar number in abundantly expressed miRNAs across samples, we detected varying miRNA expression profiles. Some miRNAs showed inconsistent or opposite dysregulation trends across different tumor tissues, including some abundantly expressed miRNA gene clusters and gene families. Wilcoxon sign-rank test for paired samples analysis revealed that abnormal miRNAs showed a higher level of variation across the seven tumor samples. We also completely surveyed abnormal miRNAs expressed in tumor and serum tissues in the mixed datasets based on the relative expression levels. Most of these miRNAs were significantly down-regulated in tumor samples, but nine abnormal miRNAs (miR-18a, 19a, 20a, 30a, 103b, 126, 126*, 192, 1287) were consistently expressed in tumor tissues and serum samples. Based on experimentally validated target mRNAs, functional enrichment analysis indicated that these abnormal miRNAs and miRNA groups (miRNA gene clusters and gene families) have important roles in multiple biological processes. Dynamic miRNA expression profiles, various abnormal miRNA profiles and complexity of the miRNA regulatory network reveal that the miRNA expression profile is a potential biomarker for classifying or detecting human disease.     

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.     

Two molecular features contribute to the Argonaute specificity for the microRNA and RNAi pathways in C. elegans

In Caenorhabditis elegans, specific Argonaute proteins are dedicated to the RNAi and microRNA pathways. To uncover how the precise Argonaute selection occurs, we designed dsRNA triggers containing both miRNA and siRNA sequences. While dsRNA carrying nucleotides mismatches can only enter the miRNA pathway, a fully complementary dsRNA successfully rescues let-7 miRNA function and initiates silencing by RNAi. We demonstrated that RDE-1 is essential for RNAi induced by the perfectly paired trigger, yet is not required for silencing by the let-7 miRNA. In contrast, ALG-1/ALG-2 are required for the miRNA function, but not for the siRNA-directed gene silencing. Finally, a dsRNA containing a bulged miRNA and a perfectly paired siRNA can enter both pathways suggesting that the sorting of small RNAs occurs after that the dsRNA trigger has been processed by Dicer. Thus, our data suggest that the selection of Argonaute proteins is affected by two molecular features: (1) the structure of the small RNA duplex; and (2) the Argonautes specific characteristics.     

Natural variation in biogenesis efficiency of individual Arabidopsis thaliana microRNAs

Like protein-coding genes, loci that produce microRNAs (miRNAs) are generally considered to be under purifying selection, consistent with miRNA polymorphisms being able to cause disease. Nevertheless, it has been hypothesized that variation in miRNA genes may contribute to phenotypic diversity. Here we demonstrate that a naturally occurring polymorphism in the MIR164A gene affects leaf shape and shoot architecture in Arabidopsis thaliana, with the effects being modified by additional loci in the genome. A single base pair substitution in the miRNA complementary sequence alters the predicted stability of the miRNA:miRNA(?) duplex. It thereby greatly reduces miRNA accumulation, probably because it interferes with precursor processing. We demonstrate that this is not a rare exception and that natural strains of Arabidopsis thaliana harbor dozens of similar polymorphisms that affect processing of a wide range of miRNA precursors. Our results suggest that natural variation in miRNA biogenesis resulting from cis mutations is a common contributor to phenotypic variation in plants.     

Characterization and 454 pyrosequencing of Major Histocompatibility Complex class I genes in the great tit reveal complexity in a passerine system

ABSTRACT: BACKGROUND: The critical role of Major Histocompatibility Complex (Mhc) genes in disease resistance and their highly polymorphic nature make them exceptional candidates for studies investigating genetic effects on survival, mate choice and conservation. Species that harbor many Mhc loci and high allelic diversity are particularly intriguing as they are potentially under strong selection and studies of such species provide valuable information as to the mechanisms maintaining Mhc diversity. However comprehensive genotyping of complex multilocus systems has been a major challenge to date with the result that little is known about the consequences of this complexity in terms of fitness effects and disease resistance. RESULTS: In this study, we genotyped the Mhc class I exon 3 of the great tit (Parus major) from two nest-box breeding populations near Oxford, UK that have been monitored for decades. Characterization of Mhc class I exon 3 was adopted and bidirectional sequencing was carried using the 454 sequencing platform. Full analysis of sequences through a stepwise variant validation procedure allowed reliable typing of more than 800 great tits based on 214,357 reads; from duplicates we estimated the repeatability of typing as 0.94. A total of 862 alleles were detected, and the presence of at least 16 functional loci was shown - the highest number characterized in a wild bird species. Finally, the functional alleles were grouped into 17 supertypes based on their antigen binding affinities. CONCLUSIONS: We found extreme complexity at the Mhc class I of the great tit both in terms of allelic diversity and gene number. The presence of many functional loci was shown, together with a pseudogene family and putatively non-functional alleles; there was clear evidence that functional alleles were under strong balancing selection. This study is the first step towards an in-depth analysis of this gene complex in this species, which will help understanding how parasite-mediated and sexual selection shape and maintain host genetic variation in nature. We believe that study systems like ours can make important contributions to the field of evolutionary biology and emphasize the necessity of integrating long-term field-based studies with detailed genetic analysis to unravel complex evolutionary processes.     

Interactive Web-based Annotation of Plant MicroRNAs with iwa-miRNA

MicroRNAs (miRNAs) are important regulators of gene expression. The large-scale detection and profiling of miRNAs have been accelerated with the development of high-throughput small RNA sequencing (sRNA-Seq) techniques and bioinformatics tools. However, generating high-quality comprehensive miRNA annotations remains challenging due to the intrinsic complexity of sRNA-Seq data and inherent limitations of existing miRNA prediction tools. Here, we present iwa-miRNA, a Galaxy-based framework that can facilitate miRNA annotation in plant species by combining computational analysis and manual curation. iwa-miRNA is specifically designed to generate a comprehensive list of miRNA candidates, bridging the gap between already annotated miRNAs provided by public miRNA databases and new predictions from sRNA-Seq datasets. It can also assist users in selecting promising miRNA candidates in an interactive mode, contributing to the accessibility and reproducibility of genome-wide miRNA annotation. iwa-miRNA is user-friendly and can be easily deployed as a web application for researchers without programming experience. With flexible, interactive, and easy-to-use features, iwa-miRNA is a valuable tool for the annotation of miRNAs in plant species with reference genomes. We also illustrate the application of iwa-miRNA for miRNA annotation using data from plant species with varying genomic complexity. The source codes and web server of iwa-miRNA are freely accessible at http://iwa-miRNA.omicstudio.cloud/.