The evolution of duplicate gene expression in mammalian organs.

Gene duplications generate genomic raw material that allows the emergence of novel functions, likely facilitating adaptive evolutionary innovations. However, global assessments of the functional and evolutionary relevance of duplicate genes in mammals were until recently limited by the lack of appropriate comparative data. Here, we report a large-scale study of the expression evolution of DNA-based functional gene duplicates in three major mammalian lineages (placental mammals, marsupials, egg-laying monotremes) and birds, on the basis of RNA sequencing (RNA-seq) data from nine species and eight organs. We observe dynamic changes in tissue expression preference of paralogs with different duplication ages, suggesting differential contribution of paralogs to specific organ functions during vertebrate evolution. Specifically, we show that paralogs that emerged in the common ancestor of bony vertebrates are enriched for genes with brain-specific expression and provide evidence for differential forces underlying the preferential emergence of young testis- and liver-specific expressed genes. Further analyses uncovered that the overall spatial expression profiles of gene families tend to be conserved, with several exceptions of pronounced tissue specificity shifts among lineage-specific gene family expansions. Finally, we trace new lineage-specific genes that may have contributed to the specific biology of mammalian organs, including the little-studied placenta. Overall, our study provides novel and taxonomically broad evidence for the differential contribution of duplicate genes to tissue-specific transcriptomes and for their importance for the phenotypic evolution of vertebrates.

Sex-biased microRNA expression in mammals and birds reveals underlying regulatory mechanisms and a role in dosage compensation.

Sexual dimorphism depends on sex-biased gene expression, but the contributions of microRNAs (miRNAs) have not been globally assessed. We therefore produced an extensive small RNA sequencing data set to analyze male and female miRNA expression profiles in mouse, opossum, and chicken. Our analyses uncovered numerous cases of somatic sex-biased miRNA expression, with the largest proportion found in the mouse heart and liver. Sex-biased expression is explained by miRNA-specific regulation, including sex-biased chromatin accessibility at promoters, rather than piggybacking of intronic miRNAs on sex-biased protein-coding genes. In mouse, but not opossum and chicken, sex bias is coordinated across tissues such that autosomal testis-biased miRNAs tend to be somatically male-biased, whereas autosomal ovary-biased miRNAs are female-biased, possibly due to broad hormonal control. In chicken, which has a Z/W sex chromosome system, expression output of genes on the Z Chromosome is expected to be male-biased, since there is no global dosage compensation mechanism that restores expression in ZW females after almost all genes on the W Chromosome decayed. Nevertheless, we found that the dominant liver miRNA, miR-122-5p, is Z-linked but expressed in an unbiased manner, due to the unusual retention of a W-linked copy. Another Z-linked miRNA, the male-biased miR-2954-3p, shows conserved preference for dosage-sensitive genes on the Z Chromosome, based on computational and experimental data from chicken and zebra finch, and acts to equalize male-to-female expression ratios of its targets. Unexpectedly, our findings thus establish miRNA regulation as a novel gene-specific dosage compensation mechanism.

The evolution of lncRNA repertoires and expression patterns in tetrapods.

Only a very small fraction of long noncoding RNAs (lncRNAs) are well characterized. The evolutionary history of lncRNAs can provide insights into their functionality, but the absence of lncRNA annotations in non-model organisms has precluded comparative analyses. Here we present a large-scale evolutionary study of lncRNA repertoires and expression patterns, in 11 tetrapod species. We identify approximately 11,000 primate-specific lncRNAs and 2,500 highly conserved lncRNAs, including approximately 400 genes that are likely to have originated more than 300 million years ago. We find that lncRNAs, in particular ancient ones, are in general actively regulated and may function predominantly in embryonic development. Most lncRNAs evolve rapidly in terms of sequence and expression levels, but tissue specificities are often conserved. We compared expression patterns of homologous lncRNA and protein-coding families across tetrapods to reconstruct an evolutionarily conserved co-expression network. This network suggests potential functions for lncRNAs in fundamental processes such as spermatogenesis and synaptic transmission, but also in more specific mechanisms such as placenta development through microRNA production.

The life history of retrocopies illuminates the evolution of new mammalian genes.

New genes contribute substantially to adaptive evolutionary innovation, but the functional evolution of new mammalian genes has been little explored at a broad scale. Previous work established mRNA-derived gene duplicates, known as retrocopies, as models for the study of new gene origination. Here we combine mammalian transcriptomic and epigenomic data to unveil the processes underlying the evolution of stripped-down retrocopies into complex new genes. We show that although some robustly expressed retrocopies are transcribed from preexisting promoters, most evolved new promoters from scratch or recruited proto-promoters in their genomic vicinity. In particular, many retrocopy promoters emerged from ancestral enhancers (or bivalent regulatory elements) or are located in CpG islands not associated with other genes. We detected 88-280 selectively preserved retrocopies per mammalian species, illustrating that these mechanisms facilitated the birth of many functional retrogenes during mammalian evolution. The regulatory evolution of originally monoexonic retrocopies was frequently accompanied by exon gain, which facilitated co-option of distant promoters and allowed expression of alternative isoforms. While young retrogenes are often initially expressed in the testis, increased regulatory and structural complexities allowed retrogenes to functionally diversify and evolve somatic organ functions, sometimes as complex as those of their parents. Thus, some retrogenes evolved the capacity to temporarily substitute for their parents during the process of male meiotic X inactivation, while others rendered parental functions superfluous, allowing for parental gene loss. Overall, our reconstruction of the "life history" of mammalian retrogenes highlights retroposition as a general model for understanding new gene birth and functional evolution.

Combinatorial Gene Regulatory Functions Underlie Ultraconserved Elements in Drosophila.

Ultraconserved elements (UCEs) are discrete genomic elements conserved across large evolutionary distances. Although UCEs have been linked to multiple facets of mammalian gene regulation their extreme evolutionary conservation remains largely unexplained. Here, we apply a computational approach to investigate this question in Drosophila, exploring the molecular functions of more than 1,500 UCEs shared across the genomes of 12 Drosophila species. Our data indicate that Drosophila UCEs are hubs for gene regulatory functions and suggest that UCE sequence invariance originates from their combinatorial roles in gene control. We also note that the gene regulatory roles of intronic and intergenic UCEs (iUCEs) are distinct from those found in exonic UCEs (eUCEs). In iUCEs, transcription factor (TF) and epigenetic factor binding data strongly support iUCE roles in transcriptional and epigenetic regulation. In contrast, analyses of eUCEs indicate that they are two orders of magnitude more likely than the expected to simultaneously include protein-coding sequence, TF-binding sites, splice sites, and RNA editing sites but have reduced roles in transcriptional or epigenetic regulation. Furthermore, we use a Drosophila cell culture system and transgenic Drosophila embryos to validate the notion of UCE combinatorial regulatory roles using an eUCE within the Hox gene Ultrabithorax and show that its protein-coding region also contains alternative splicing regulatory information. Taken together our experiments indicate that UCEs emerge as a result of combinatorial gene regulatory roles and highlight common features in mammalian and insect UCEs implying that similar processes might underlie ultraconservation in diverse animal taxa.

Evolution of Hox post-transcriptional regulation by alternative polyadenylation and microRNA modulation within 12 Drosophila genomes.

Hox genes encode a family of transcriptional regulators that operate differential developmental programs along the anteroposterior axis of bilateral animals. Regulatory changes affecting Hox gene expression are believed to have been crucial for the evolution of animal body plans. In Drosophila melanogaster, Hox expression is post-transcriptionally regulated by microRNAs (miRNAs) acting on target sites located in the 3' untranslated regions (3'UTRs) of Hox mRNAs. Notably, recent work has shown that during D. melanogaster development Hox genes produce mRNAs with variable 3'UTRs (short and long forms) in different sets of tissues as a result of alternative polyadenylation; importantly, Hox short and long 3'UTRs contain very different target sites for miRNAs. Here, we use a computational approach to explore the evolution of Hox 3'UTRs treated with especial regard to miRNA regulation. Our work is focused on the 12 Drosophila species for which genomic sequences are available and shows, first, that alternative polyadenylation of Hox transcripts is a feature shared by all drosophilids tested in the study. Second, that the regulatory impact of miRNAs is evolving very fast within the Drosophila group. Third, that in contrast to the low degree of primary sequence conservation, Hox 3'UTR regions within the group show very similar RNA topology indicating that RNA structure is under strong selective pressure. Finally, we also demonstrate that Hox alternative polyadenylation can remodel the control regions seen by miRNAs by at least two mechanisms: via adding new cis-regulatory sequences-in the form of miRNA target sites-to short 3'UTR forms as well as by modifying the regulatory impact of miRNA target sites in short 3'UTR forms through changes in RNA secondary structure caused by the use of distal polyadenylation signals.

Transposable elements: insertion pattern and impact on gene expression evolution in hominids.

Transposable elements (TEs) can affect the regulation of nearby genes through several mechanisms. Here, we examine to what extent recent TE insertions have contributed to the evolution of gene expression in hominids. We compare expression levels of human and chimpanzee orthologs and detect a weak increase in expression divergence (ED) for genes with species-specific TE insertions compared with unaffected genes. However, we show that genes with TE insertions predating the human-chimpanzee split also exhibit a similar increase in ED and therefore conclude that the increase is not due to the transcriptional influence of the TEs. These results are further confirmed by lineage-specific analysis of ED, using rhesus macaque as an outgroup: Human-chimpanzee ortholog pairs, where one ortholog has suffered TE insertion but not the other, do not show increased ED along the lineage where the insertion occurred, relative to the other lineage. We also show that genes with recent TE insertions tend to produce more alternative transcripts but find no evidence that the TEs themselves promote transcript diversity. Finally, we observe that TEs are enriched upstream relative to downstream of genes and show that this is due to insertional bias, rather than selection, because this bias is only observed in genes expressed in the germ line. This provides an alternative neutral explanation for the accumulation of TEs in upstream sequences.

Repurposing of promoters and enhancers during mammalian evolution.

Promoters and enhancers-key controllers of gene expression-have long been distinguished from each other based on their function. However, recent work suggested that common architectural and functional features might have facilitated the conversion of one type of element into the other during evolution. Here, based on cross-mammalian analyses of epigenome and transcriptome data, we provide support for this hypothesis by detecting 445 regulatory elements with signatures of activity turnover (termed P/E elements). Most events represent transformations of putative ancestral enhancers into promoters, leading to the emergence of species-specific transcribed loci or 5' exons. Distinct GC sequence compositions and stabilizing 5' splicing (U1) regulatory motif patterns may have predisposed P/E elements to regulatory repurposing, and changes in the U1 and polyadenylation signal densities and distributions likely drove the evolutionary activity switches. Our work suggests that regulatory repurposing facilitated regulatory innovation and the origination of new genes and exons during evolution.

Conserved microRNA editing in mammalian evolution, development and disease.

BACKGROUND: Mammalian microRNAs (miRNAs) are sometimes subject to adenosine-to-inosine RNA editing, which can lead to dramatic changes in miRNA target specificity or expression levels. However, although a few miRNAs are known to be edited at identical positions in human and mouse, the evolution of miRNA editing has not been investigated in detail. In this study, we identify conserved miRNA editing events in a range of mammalian and non-mammalian species. RESULTS: We demonstrate deep conservation of several site-specific miRNA editing events, including two that date back to the common ancestor of mammals and bony fishes some 450 million years ago. We also find evidence of a recent expansion of an edited miRNA family in placental mammals and show that editing of these miRNAs is associated with changes in target mRNA expression during primate development and aging. While global patterns of miRNA editing tend to be conserved across species, we observe substantial variation in editing frequencies depending on tissue, age and disease state: editing is more frequent in neural tissues compared to heart, kidney and testis; in older compared to younger individuals; and in samples from healthy tissues compared to tumors, which together suggests that miRNA editing might be associated with a reduced rate of cell proliferation. CONCLUSIONS: Our results show that site-specific miRNA editing is an evolutionarily conserved mechanism, which increases the functional diversity of mammalian miRNA transcriptomes. Furthermore, we find that although miRNA editing is rare compared to editing of long RNAs, miRNAs are greatly overrepresented among conserved editing targets.

Germ cell-specific targeting of DICER or DGCR8 reveals a novel role for endo-siRNAs in the progression of mammalian spermatogenesis and male fertility.

Small non-coding RNAs act as critical regulators of gene expression and are essential for male germ cell development and spermatogenesis. Previously, we showed that germ cell-specific inactivation of Dicer1, an endonuclease essential for the biogenesis of micro-RNAs (miRNAs) and endogenous small interfering RNAs (endo-siRNAs), led to complete male infertility due to alterations in meiotic progression, increased spermatocyte apoptosis and defects in the maturation of spermatozoa. To dissect the distinct physiological roles of miRNAs and endo-siRNAs in spermatogenesis, we compared the testicular phenotype of mice with Dicer1 or Dgcr8 depletion in male germ cells. Dgcr8 mutant mice, which have a defective miRNA pathway while retaining an intact endo-siRNA pathway, were also infertile and displayed similar defects, although less severe, to Dicer1 mutant mice. These included cumulative defects in meiotic and haploid phases of spermatogenesis, resulting in oligo-, terato-, and azoospermia. In addition, we found by RNA sequencing of purified spermatocytes that inactivation of Dicer1 and the resulting absence of miRNAs affected the fine tuning of protein-coding gene expression by increasing low level gene expression. Overall, these results emphasize the essential role of miRNAs in the progression of spermatogenesis, but also indicate a role for endo-siRNAs in this process.

Evolution of the correlation between expression divergence and protein divergence in mammals.

Divergence of protein sequences and gene expression patterns are two fundamental mechanisms that generate organismal diversity. Here, we have used genome and transcriptome data from eight mammals and one bird to study the positive correlation of these two processes throughout mammalian evolution. We demonstrate that the correlation is stable over time and most pronounced in neural tissues, which indicates that it is the result of strong negative selection. The correlation is not driven by genes with specific functions and may instead best be viewed as an evolutionary default state, which can nevertheless be evaded by certain gene types. In particular, genes with developmental and neural functions are skewed toward changes in gene expression, consistent with selection against pleiotropic effects associated with changes in protein sequences. Surprisingly, we find that the correlation between expression divergence and protein divergence is not explained by between-gene variation in expression level, tissue specificity, protein connectivity, or other investigated gene characteristics, suggesting that it arises independently of these gene traits. The selective constraints on protein sequences and gene expression patterns also fluctuate in a coordinate manner across phylogenetic branches: We find that gene-specific changes in the rate of protein evolution in a specific mammalian lineage tend to be accompanied by similar changes in the rate of expression evolution. Taken together, our findings highlight many new aspects of the correlation between protein divergence and expression divergence, and attest to its role as a fundamental property of mammalian genome evolution.

A selection index for gene expression evolution and its application to the divergence between humans and chimpanzees.

The importance of gene regulation in animal evolution is a matter of long-standing interest, but measuring the impact of selection on gene expression has proven a challenge. Here, we propose a selection index of gene expression as a straightforward method for assessing the mode and strength of selection operating on gene expression levels. The index is based on the widely used McDonald-Kreitman test and requires the estimation of four quantities: the within-species and between-species expression variances as well as the sequence heterozygosity and divergence of neutrally evolving sequences. We apply the method to data from human and chimpanzee lymphoblastoid cell lines and show that gene expression is in general under strong stabilizing selection. We also demonstrate how the same framework can be used to estimate the proportion of adaptive gene expression evolution.

The accumulation of gene regulation through time.

Gene expression is governed by an intricate combination of transcription factors (TFs), microRNAs (miRNAs), splicing factors, and other regulators. Genes cannot support infinitely complex regulation due to sequence constraints and the increased likelihood of harmful errors. However, the upper limit of regulatory complexity in the genome is not known. Here, we provide evidence that human genes are currently not operating at their maximum capacity in terms of gene regulation. We analyze genes spanning the full spectrum of eukaryote evolution, from primate-specific genes to genes present in the eukaryote ancestor, and show that older genes tend to be bound by more TFs, have more conserved upstream sequences, generate more alternative isoforms, house more miRNA targets, and are more likely to be affected by nonsense-mediated decay and RNA editing. These results cannot be explained by overrepresentation of certain functional categories among younger or older genes. Furthermore, the increase in complexity is continuous over evolutionary time, without signs of saturation, leading to the conclusion that most genes, at least in the human genome, have the capacity to evolve even more complex gene regulation in the future.