Evolution of human gene expression

During evolution, biological differences between species can arise not only due to structural differences between genes, but also following changes in how, where and when genes are active. However, we know much less about this second aspect, because large-scale comparative transcriptomics only became feasible relatively recently. In this thesis, I will therefore investigate several aspects of gene expression evolution, with emphasis on our own species. A first step to understanding regulatory evolution is to determine how variation in gene expression is created. Transposable elements (TEs) are genomic parasites that can affect their host genome in a number of ways, including gene expression. In Chapter 2, I investigate to what extent transposable elements (TEs) have contributed to expression differences between humans and chimpanzees. Once expression variation has been established, a combination of selection and drift will decide which variants are passed on to future generations. It is of particular interest to identify changes that were established through positive selection, as these are adaptive. In Chapter 3, I describe a new method to detect positive selection acting on gene expression and apply it to data from humans and chimpanzees. Human gene expression is regulated through several mechanisms associated with transcription and post-transcriptional processing. In Chapter 4, I consider the long-term evolution of the human genome and investigate whether genes have reached their maximum capacity in terms of regulatory complexity. Finally, in Chapter 5, I explore the relationship between gene regulation and sequence conservation by identifying and analysing extremely conserved elements in the genome of the fruit fly Drosophila melanogaster

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

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

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

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 life cycle of Drosophila orphan genes

Orphans are genes restricted to a single phylogenetic lineage and emerge at high rates. While this predicts an accumulation of genes, the gene number has remained remarkably constant through evolution. This paradox has not yet been resolved. Because orphan genes have been mainly analyzed over long evolutionary time scales, orphan loss has remained unexplored. Here we study the patterns of orphan turnover among close relatives in the Drosophila obscura group. We show that orphans are not only emerging at a high rate, but that they are also rapidly lost. Interestingly, recently emerged orphans are more likely to be lost than older ones. Furthermore, highly expressed orphans with a strong male-bias are more likely to be retained. Since both lost and retained orphans show similar evolutionary signatures of functional conservation, we propose that orphan loss is not driven by high rates of sequence evolution, but reflects lineage specific functional requirements.Comment: 47 pages, 19 figure

ProteinHistorian: Tools for the Comparative Analysis of Eukaryote Protein Origin

The evolutionary history of a protein reflects the functional history of its ancestors. Recent phylogenetic studies identified distinct evolutionary signatures that characterize proteins involved in cancer, Mendelian disease, and different ontogenic stages. Despite the potential to yield insight into the cellular functions and interactions of proteins, such comparative phylogenetic analyses are rarely performed, because they require custom algorithms. We developed ProteinHistorian to make tools for performing analyses of protein origins widely available. Given a list of proteins of interest, ProteinHistorian estimates the phylogenetic age of each protein, quantifies enrichment for proteins of specific ages, and compares variation in protein age with other protein attributes. ProteinHistorian allows flexibility in the definition of protein age by including several algorithms for estimating ages from different databases of evolutionary relationships. We illustrate the use of ProteinHistorian with three example analyses. First, we demonstrate that proteins with high expression in human, compared to chimpanzee and rhesus macaque, are significantly younger than those with human-specific low expression. Next, we show that human proteins with annotated regulatory functions are significantly younger than proteins with catalytic functions. Finally, we compare protein length and age in many eukaryotic species and, as expected from previous studies, find a positive, though often weak, correlation between protein age and length. ProteinHistorian is available through a web server with an intuitive interface and as a set of command line tools; this allows biologists and bioinformaticians alike to integrate these approaches into their analysis pipelines. ProteinHistorian's modular, extensible design facilitates the integration of new datasets and algorithms. The ProteinHistorian web server, source code, and pre-computed ages for 32 eukaryotic genomes are freely available under the GNU public license at http://lighthouse.ucsf.edu/ProteinHistorian/

Adaptive sequence evolution is driven by biotic stress in a pair of orchid species (Dactylorhiza) with distinct ecological optima

This work was funded by an Austrian Science Fund (FWF) project (Y661-B16) awarded to OP and a Marie Curie IEF fellowship (PERG-GA-2011-299608-TRANSADAPTATION) to FB.The orchid family is the largest in the angiosperms, but little is known about the molecular basis of the significant variation they exhibit. We investigate here the transcriptomic divergence between two European terrestrial orchids, Dactylorhiza incarnata and D. fuchsii, and integrate these results in the context of their distinct ecologies that we also document. Clear signals of lineage-specific adaptive evolution of protein-coding sequences are identified, notably targeting elements of biotic defence, including both physical and chemical adaptations in the context of divergent pools of pathogens and herbivores. In turn, a substantial regulatory divergence between the two species appears linked to adaptation/acclimation to abiotic conditions. Several of the pathways affected by differential expression are also targeted by deviating post-transcriptional regulation via sRNAs. Finally, Dactylorhiza incarnata appears to suffer from insufficient sRNA control over the activity of RNA-dependent DNA polymerase, resulting in increased activity of class I transposable elements and, over time, in larger genome size than that of D. fuchsii. The extensive molecular divergence between the two species suggests significant genomic and transcriptomic shock in their hybrids and offers insights into the difficulty of coexistence at the homoploid level. Altogether, biological response to selection, accumulated during the history of these orchids, appears governed by their microenvironmental context, in which biotic and abiotic pressures act synergistically to shape transcriptome structure, expression and regulation.Publisher PDFPeer reviewe

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