Evolutionary simulations to detect functional lineage-specific genes

Motivation: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. Results: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (∼1400 bp), restricted to hominoids, can be predicted to be functional in ∼94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. Availability: ReEVOLVER is available at . Contact: [email protected] Supplementary Data: Supplementary Data are available at Bioinformatics onlin

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

The Urea Carboxylase and Allophanate Hydrolase Activities of Urea Amidolyase Are Functionally Independent

Urea amidolyase (UAL) is a multifunctional biotin-dependent enzyme that contributes to both bacterial and fungal pathogenicity by catalyzing the ATP-dependent cleavage of urea into ammonia and CO2. UAL is comprised of two enzymatic components: urea carboxylase (UC) and allophanate hydrolase (AH). These enzyme activities are encoded on separate but proximally related genes in prokaryotes while, in most fungi, they are encoded by a single gene that produces a fusion enzyme on a single polypeptide chain. It is unclear whether the UC and AH activities are connected through substrate channeling or other forms of direct communication. Here, we use multiple biochemical approaches to demonstrate that there is no substrate channeling or interdomain/intersubunit communication between UC and AH. Neither stable nor transient interactions can be detected between prokaryotic UC and AH and the catalytic efficiencies of UC and AH are independent of one another. Furthermore, an artificial fusion of UC and AH does not significantly alter the AH enzyme activity or catalytic efficiency. These results support the surprising functional independence of AH from UC in both the prokaryotic and fungal UAL enzymes and serve as an important reminder that the evolution of multifunctional enzymes through gene fusion events does not always correlate with enhanced catalytic function

RNA from a simple-tandem repeat is required for sperm maturation and male fertility in Drosophila melanogaster.

Tandemly-repeated DNAs, or satellites, are enriched in heterochromatic regions of eukaryotic genomes and contribute to nuclear structure and function. Some satellites are transcribed, but we lack direct evidence that specific satellite RNAs are required for normal organismal functions. Here, we show satellite RNAs derived from AAGAG tandem repeats are transcribed in many cells throughout Drosophila melanogaster development, enriched in neurons and testes, often localized within heterochromatic regions, and important for viability. Strikingly, we find AAGAG transcripts are necessary for male fertility, and that AAGAG RNA depletion results in defective histone-protamine exchange, sperm maturation and chromatin organization. Since these events happen late in spermatogenesis when the transcripts are not detected, we speculate that AAGAG RNA in primary spermatocytes 'primes' post-meiosis steps for sperm maturation. In addition to demonstrating essential functions for AAGAG RNAs, comparisons between closely related Drosophila species suggest that satellites and their transcription evolve quickly to generate new functions

Convergent origination of a Drosophila-like dosage compensation mechanism in a reptile lineage

Sex chromosomes differentiated from different ancestral autosomes in various vertebrate lineages. Here, we trace the functional evolution of the XY Chromosomes of the green anole lizard (Anolis carolinensis), on the basis of extensive high-throughput genome, transcriptome and histone modification sequencing data and revisit dosage compensation evolution in representative mammals and birds with substantial new expression data. Our analyses show that Anolis sex chromosomes represent an ancient XY system that originated at least ≈160 million years ago in the ancestor of Iguania lizards, shortly after the separation from the snake lineage. The age of this system approximately coincides with the ages of the avian and two mammalian sex chromosomes systems. To compensate for the almost complete Y Chromosome degeneration, X-linked genes have become twofold up-regulated, restoring ancestral expression levels. The highly efficient dosage compensation mechanism of Anolis represents the only vertebrate case identified so far to fully support Ohno’s original dosage compensation hypothesis. Further analyses reveal that X up-regulation occurs only in males and is mediated by a male-specific chromatin machinery that leads to global hyperacetylation of histone H4 at lysine 16 specifically on the X Chromosome. The green anole dosage compensation mechanism is highly reminiscent of that of the fruit fly, Drosophila melanogaster. Altogether, our work unveils the convergent emergence of a Drosophila-like dosage compensation mechanism in an ancient reptilian sex chromosome system and highlights that the evolutionary pressures imposed by sex chromosome dosage reductions in different amniotes were resolved in fundamentally different ways

Exon-phase symmetry and intrinsic structural disorder promote modular evolution in the human genome

A key signature of module exchange in the genome is phase symmetry of exons, suggestive of exon shuffling events that occurred without disrupting translation reading frame. At the protein level, intrinsic structural disorder may be another key element because disordered regions often serve as functional elements that can be effectively integrated into a protein structure. Therefore, we asked whether exon-phase symmetry in the human genome and structural disorder in the human proteome are connected, signalling such evolutionary mechanisms in the assembly of multi-exon genes. We found an elevated level of structural disorder of regions encoded by symmetric exons and a preferred symmetry of exons encoding for mostly disordered regions (>70% predicted disorder). Alternatively spliced symmetric exons tend to correspond to the most disordered regions. The genes of mostly disordered proteins (>70% predicted disorder) tend to be assembled from symmetric exons, which often arise by internal tandem duplications. Preponderance of certain types of short motifs (e.g. SH3-binding motif) and domains (e.g. high-mobility group domains) suggests that certain disordered modules have been particularly effective in exon-shuffling events. Our observations suggest that structural disorder has facilitated modular assembly of complex genes in evolution of the human genome. © 2013 The Author(s)

The Pseudogenes of Barley

Pseudogenes have a reputation of being ‘evolutionary relics’ or ‘junk DNA’. While they are well characterized in mammals, studies in more complex plant genomes were so far hampered by the absence of reference genome sequences. Barley is one of the economically most important cereals and has a genome size of 5.1 Gb. With the first high-quality genome reference assembly available for a Triticeae crop, we conducted a whole genome assessment of pseudogenes on the barley genome. We identified, characterized, and classified 89,440 gene fragments and pseudogenes, scattered along the chromosomes with occasional hotspots and higher densities at the chromosome ends. Full-length pseudogenes (11,015) have preferentially retained their exon-intron structure. Retrotransposition of processed mRNAs only plays a marginal role in their creation. However, the distribution of retroposed pseudogenes reflects the Rabl configuration of barley chromosomes and thus hints towards founding mechanisms. While defense-response related parent genes were found under-represented in cultivated barley, we detected several defense related pseudogenes in wild barley accessions. 7.2% of the pseudogenes are transcriptionally active and may potentially adopt new regulatory roles.The barley genome is rich in pseudogenes and small gene fragments mainly located towards chromosome tips or as tandemly repeated units. Our results indicate non-random duplication and pseudogenization preferences and improve our understanding of gene birth and death dynamics in large plant genomes and the mechanisms that lead to evolutionary innovations201

The cellular origins and evolution of the maternal-fetal interface in mammals : Final report DFG Project number: 433034324

The placenta is a fascinating evolutionary innovation—an entirely new organ that arose in mammals to enable pregnancy and mediate physiological exchange between mother and fetus. Yet, how this organ originated, diversified, and is regulated at the cellular level across species remains poorly understood. In this project, we set out to explore the evolution and development of the placenta through a deep, comparative lens, using single-cell transcriptomic and epigenomic profiling across a diverse set of mammals. We generated nearly 400,000 single-nucleus transcriptomes from the maternal and fetal components of the placenta in nine species: human, marmoset, mouse, rat, guinea pig, rabbit, sheep, horse, and opossum. These datasets span key stages of pregnancy and enabled the generation of high-resolution, cross-species cell atlases of the maternal-fetal interface. In parallel, we produced high-quality epigenomic data for the mouse placenta, to understand the gene regulatory architecture driving cell differentiation. Our analyses revealed striking evolutionary innovations in murid rodents, such as the emergence of novel trophoblast cell (sub)types—including the split of syncytiotrophoblast into two layer-specific subtypes, sinusoidal giant cells, and spongiotrophoblasts—that are absent in humans, primates, and other mammals. Using gene expression similarity and phylogenetic mapping across six core species, we reconstructed the evolutionary history of trophoblast cell types, showing that new cell types emerged in the murid lineage ~27 million years ago, fundamentally reshaping the placental interface in rodents. To investigate the regulatory mechanisms behind these changes, we developed novel methods to integrate the transcriptomic and epigenomic data, combining autoencoder-based neural networks, tailored statistical models, and custom strategies to overcome challenges posed by large developmental time gaps. This enabled us to trace the regulatory divergence between syncytiotrophoblast subtypes, highlighting transcription factors such as CREB5 and Jun-AP1 as key drivers of subtype specification. These TFs target genes involved in cytoskeletal remodeling, suggesting a functional link between gene regulation and the morphological adaptation of trophoblasts in fetal versus maternal regions of the placenta. Our findings position the placenta as a unique outlier among mammalian organs—defined not by the conservation of ancestral cell types, but by remarkable cell-type innovation. This exceptional evolutionary plasticity makes it an ideal system for addressing one of the most fundamental questions in biology: how do new cell types evolve

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

Functional diversification of duplicate genes through subcellular adaptation of encoded proteins

Analysis of the subcellular localization patterns of duplicate genes revealed that protein subcellular adaptation represents a common mechanism for the functional diversification of duplicate genes