On the Regulatory Evolution of New Genes Throughout Their Life History

The lawc gene emerged de novo from conserved genomic elements and acquired a broad expression pattern in Drosophila

It has recently become evident that the de novo emergence of genes is widespread and documented for a variety of organisms. De novo genes frequently emerge in proximity to existing genes, forming gene overlaps. Here, we present an analysis of the evolutionary history of a putative de novo gene, lawc, which overlaps with the conserved Trf2 gene, which encodes a general transcription factor in Drosophila melanogaster. We demonstrate that lawc emerged approximately 68 million years ago in the 5'-untranslated region of Trf2 and displays an extensive spatiotemporal expression pattern. One of the most remarkable features of the lawc evolutionary history is that its emergence was facilitated by the engagement of Drosophilidae-specific short highly conserved regions located in Trf2 introns. This represents a unique example of putative de novo gene birth involving conserved DNA regions localized in introns of conserved gene. The observed lawc expression pattern may be due to overlap of lawc with the 5'-untranslated region of Trf2. This study not only enriches our understanding of gene evolution but also highlights the complex interplay between genetic conservation and innovation.

Transposable element exonization generates a reservoir of evolving and functional protein isoforms

Alternative splicing enhances protein diversity in different ways, including through exonization of transposable elements (TEs). Recent transcriptomic analyses identified thousands of unannotated spliced transcripts with exonizing TEs, but their contribution to the proteome and biological relevance remains unclear. Here, we use transcriptome assembly, ribosome profiling, and proteomics to describe a population of 1,227 unannotated TE exonizing isoforms generated by mRNA splicing and recurrent in human populations. Despite being shorter and lowly expressed, these isoforms are shared between individuals and efficiently translated. Functional analyses show stable expression, specific cellular localization, and, in some cases, modified functions. Exonized TEs are rich in ancient genes, whereas the involved splice sites are recent and can be evolutionarily conserved. In addition, exonized TEs contribute to the secondary structure of the emerging isoforms, supporting their functional relevance. We conclude that TE-spliced isoforms represent a diversity reservoir of functional proteins on which natural selection can act.

The three-dimensional genome drives the evolution of asymmetric gene duplicates via enhancer capture-divergence

Previous evolutionary models of duplicate gene evolution have overlooked the pivotal role of genome architecture. Here, we show that proximity-based regulatory recruitment by distally duplicated genes is an efficient mechanism for modulating tissue-specific production of preexisting proteins. By leveraging genomic asymmetries, we performed a coexpression analysis on Drosophila melanogaster tissue data to show the generality of enhancer capture-divergence (ECD) as a significant evolutionary driver of asymmetric, distally duplicated genes. We use the recently evolved gene HP6/Umbrea as an example of the ECD process. By assaying genome-wide chromosomal conformations in multiple Drosophila species, we show that HP6/Umbrea was inserted near a preexisting, long-distance three-dimensional genomic interaction. We then use this data to identify a newly found enhancer (FLEE1), buried within the coding region of the highly conserved, essential gene MFS18, that likely neofunctionalized HP6/Umbrea. Last, we demonstrate ancestral transcriptional coregulation of HP6/Umbrea's future insertion site, illustrating how enhancer capture provides a highly evolvable, one-step solution to Ohno's dilemma.

De novo gene integration into regulatory networks via interaction with conserved genes in peach

De novo genes can evolve "from scratch" from noncoding sequences, acquiring novel functions in organisms and integrating into regulatory networks during evolution to drive innovations in important phenotypes and traits. However, identifying de novo genes is challenging, as it requires high-quality genomes from closely related species. According to the comparison with nine closely related Prunus genomes, we determined at least 178 de novo genes in P. persica "baifeng". The distinct differences were observed between de novo and conserved genes in gene characteristics and expression patterns. Gene ontology enrichment analysis suggested that Type I de novo genes originated from sequences related to plastid modification functions, while Type II genes were inferred to have derived from sequences related to reproductive functions. Finally, transcriptome sequencing across different tissues and developmental stages suggested that de novo genes have been evolutionarily recruited into existing regulatory networks, playing important roles in plant growth and development, which was also supported by WGCNA analysis and quantitative trait loci data. This study lays the groundwork for future research on the origins and functions of genes in Prunus and related taxa.

New genes driven by segmental duplications share a testis-specific expression pattern in the chromosome-level genome assembly of tree sparrow

Based on a chromosome-level genome assembly, a burst of new genes with different structures but a similar testis-specific expression pattern was detected in tree sparrow.

A comprehensive evolutionary scenario for the origin and neofunctionalization of the Drosophila speciation gene Odysseus (OdsH)

Odysseus (OdsH) was the first speciation gene described in Drosophila related to hybrid sterility in offspring of mating between Drosophila mauritiana and Drosophila simulans. Its origin is attributed to the duplication of the gene unc-4 in the subgenus Sophophora. By using a much larger sample of Drosophilidae species, we showed that contrary to what has been previously proposed, OdsH origin occurred 62 MYA. Evolutionary rates, expression, and transcription factor-binding sites of OdsH evidence that it may have rapidly experienced neofunctionalization in male sexual functions. Furthermore, the analysis of the OdsH peptide allowed the identification of mutations of D. mauritiana that could result in incompatibility in hybrids. In order to find if OdsH could be related to hybrid sterility, beyond Sophophora, we explored the expression of OdsH in Drosophila arizonae and Drosophila mojavensis, a pair of sister species with incomplete reproductive isolation. Our data indicated that OdsH expression is not atypical in their male-sterile hybrids. In conclusion, we have proposed that the origin of OdsH occurred earlier than previously proposed, followed by neofunctionalization. Our results also suggested that its role as a speciation gene might be restricted to D. mauritiana and D. simulans.

Four classic "de novo" genes all have plausible homologs and likely evolved from retro-duplicated or pseudogenic sequences

Despite being previously regarded as extremely unlikely, the idea that entirely novel protein-coding genes can emerge from non-coding sequences has gradually become accepted over the past two decades. Examples of "de novo origination", resulting in lineage-specific "orphan" genes, lacking coding orthologs, are now produced every year. However, many are likely cases of duplicates that are difficult to recognize. Here, I re-examine the claims and show that four very well-known examples of genes alleged to have emerged completely "from scratch"- FLJ33706 in humans, Goddard in fruit flies, BSC4 in baker's yeast and AFGP2 in codfish-may have plausible evolutionary ancestors in pre-existing genes. The first two are likely highly diverged retrogenes coding for regulatory proteins that have been misidentified as orphans. The antifreeze glycoprotein, moreover, may not have evolved from repetitive non-genic sequences but, as in several other related cases, from an apolipoprotein that could have become pseudogenized before later being reactivated. These findings detract from various claims made about de novo gene birth and show there has been a tendency not to invest the necessary effort in searching for homologs outside of a very limited syntenic or phylostratigraphic methodology. A robust approach is used for improving detection that draws upon similarities, not just in terms of statistical sequence analysis, but also relating to biochemistry and function, to obviate notable failures to identify homologs.

Integrating massive RNA-seq data to elucidate transcriptome dynamics in Drosophila melanogaster

The volume of ribonucleic acid (RNA)-seq data has increased exponentially, providing numerous new insights into various biological processes. However, due to significant practical challenges, such as data heterogeneity, it is still difficult to ensure the quality of these data when integrated. Although some quality control methods have been developed, sample consistency is rarely considered and these methods are susceptible to artificial factors. Here, we developed MassiveQC, an unsupervised machine learning-based approach, to automatically download and filter large-scale high-throughput data. In addition to the read quality used in other tools, MassiveQC also uses the alignment and expression quality as model features. Meanwhile, it is user-friendly since the cutoff is generated from self-reporting and is applicable to multimodal data. To explore its value, we applied MassiveQC to Drosophila RNA-seq data and generated a comprehensive transcriptome atlas across 28 tissues from embryogenesis to adulthood. We systematically characterized fly gene expression dynamics and found that genes with high expression dynamics were likely to be evolutionarily young and expressed at late developmental stages, exhibiting high nonsynonymous substitution rates and low phenotypic severity, and they were involved in simple regulatory programs. We also discovered that human and Drosophila had strong positive correlations in gene expression in orthologous organs, revealing the great potential of the Drosophila system for studying human development and disease.

Promoter evolution of mammalian gene duplicates

BACKGROUND

Gene duplication is thought to be a central process in evolution to gain new functions. The factors that dictate gene retention following duplication as well paralog gene divergence in sequence, expression and function have been extensively studied. However, relatively little is known about the evolution of promoter regions of gene duplicates and how they influence gene duplicate divergence. Here, we focus on promoters of paralog genes, comparing their similarity in sequence, in the sets of transcription factors (TFs) that bind them, and in their overall promoter architecture.

RESULTS

We observe that promoters of recent duplications display higher sequence similarity between them and that sequence similarity rapidly declines between promoters of more ancient paralogs. In contrast, similarity in cis-regulation, as measured by the set of TFs that bind promoters of both paralogs, does not simply decrease with time from duplication and is instead related to promoter architecture-paralogs with CpG Islands (CGIs) in their promoters share a greater fraction of TFs, while CGI-less paralogs are more divergent in their TF binding set. Focusing on recent duplication events and partitioning them by their duplication mechanism enables us to uncover promoter properties associated with gene retention, as well as to characterize the evolution of promoters of newly born genes: In recent retrotransposition-mediated duplications, we observe asymmetry in cis-regulation of paralog pairs: Retrocopy genes are lowly expressed and their promoters are bound by fewer TFs and are depleted of CGIs, in comparison with the original gene copy. Furthermore, looking at recent segmental duplication regions in primates enable us to compare successful retentions versus loss of duplicates, showing that duplicate retention is associated with fewer TFs and with CGI-less promoter architecture.

CONCLUSIONS

In this work, we profiled promoters of gene duplicates and their inter-paralog divergence. We also studied how their characteristics are associated with duplication time and duplication mechanism, as well as with the fate of these duplicates. These results underline the importance of cis-regulatory mechanisms in shaping the evolution of new genes and their fate following duplication.

Decoding the Spermatogenesis Program: New Insights from Transcriptomic Analyses

Spermatogenesis is a complex differentiation process coordinated spatiotemporally across and along seminiferous tubules. Cellular heterogeneity has made it challenging to obtain stage-specific molecular profiles of germ and somatic cells using bulk transcriptomic analyses. This has limited our ability to understand regulation of spermatogenesis and to integrate knowledge from model organisms to humans. The recent advancement of single-cell RNA-sequencing (scRNA-seq) technologies provides insights into the cell type diversity and molecular signatures in the testis. Fine-grained cell atlases of the testis contain both known and novel cell types and define the functional states along the germ cell developmental trajectory in many species. These atlases provide a reference system for integrated interspecies comparisons to discover mechanistic parallels and to enable future studies. Despite recent advances, we currently lack high-resolution data to probe germ cell-somatic cell interactions in the tissue environment, but the use of highly multiplexed spatial analysis technologies has begun to resolve this problem. Taken together, recent single-cell studies provide an improvedunderstanding of gametogenesis to examine underlying causes of infertility and enable the development of new therapeutic interventions.

Spatiotemporal Regulation of a Single Adaptively Evolving Trans-Regulatory Element Contributes to Spermatogenetic Expression Divergence in Drosophila

Due to extensive pleiotropy, trans-acting elements are often thought to be evolutionarily constrained. While the impact of trans-acting elements on gene expression evolution has been extensively studied, relatively little is understood about the contribution of a single trans regulator to interspecific expression and phenotypic divergence. Here, we disentangle the effects of genomic context and miR-983, an adaptively evolving young microRNA, on expression divergence between Drosophila melanogaster and D. simulans. We show miR-983 effects promote interspecific expression divergence in testis despite its antagonism with the often-predominant context effects. Single-cyst RNA-seq reveals that distinct sets of genes gain and lose miR-983 influence under disruptive or diversifying selection at different stages of spermatogenesis, potentially helping minimize antagonistic pleiotropy. At the round spermatid stage, the effects of miR-983 are weak and distributed, coincident with the transcriptome undergoing drastic expression changes. Knocking out miR-983 causes reduced sperm length with increased within-individual variation in D. melanogaster but not in D. simulans, and the D. melanogaster knockout also exhibits compromised sperm defense ability. Our results provide empirical evidence for the resolution of antagonistic pleiotropy and also have broad implications for the function and evolution of new trans regulators.

Transcriptional and post-transcriptional regulation of young genes in plants

Background

New genes continuously emerge from non-coding DNA or by diverging from existing genes, but most of them are rapidly lost and only a few become fixed within the population. We hypothesized that young genes are subject to transcriptional and post-transcriptional regulation to limit their expression and minimize their exposure to purifying selection.

Results

We performed a protein-based homology search across the tree of life to determine the evolutionary age of protein-coding genes present in the rice genome. We found that young genes in rice have relatively low expression levels, which can be attributed to distal enhancers, and closed chromatin conformation at their transcription start sites (TSS). The chromatin in TSS regions can be re-modeled in response to abiotic stress, indicating conditional expression of young genes. Furthermore, transcripts of young genes in Arabidopsis tend to be targeted by nonsense-mediated RNA decay, presenting another layer of regulation limiting their expression.

Conclusions

These data suggest that transcriptional and post-transcriptional mechanisms contribute to the conditional expression of young genes, which may alleviate purging selection while providing an opportunity for phenotypic exposure and functionalization.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01339-7.

Orphan gene in Littorina: An unexpected role of symbionts in the host evolution

Mechanisms of reproductive isolation between closely related sympatric species are of high evolutionary significance as they may function as initial drivers of speciation and protect species integrity afterwards. Proteins involved in the establishment of reproductive barriers often evolve fast and may be key players in cessation of gene flow between the incipient species. The five Atlantic Littorina (Neritrema) species represent a notable example of recent radiation. The geographic ranges of these young species largely overlap and the mechanisms of reproductive isolation are poorly understood. In this study, we performed a detailed analysis of the reproductive protein LOSP, previously identified in Littorina. We showed that this protein is evolutionary young and taxonomically restricted to the genus Littorina. It has high sequence variation both within and between Littorina species, which is compatible with its presumable role in the reproductive isolation. The strongest differences in the LOSP structure were detected between Littorina subgenera with distinctive repetitive motifs present exclusively in the Neritrema species, but not in L. littorea. Moreover, the sequence of these repetitive structural elements demonstrates a high homology with genetic elements of bacteria, identified as components of Littorina associated microbiomes. We suggest that these elements were acquired from a symbiotic bacterial donor via horizontal genetic transfer (HGT), which is indirectly confirmed by the presence of multiple transposable elements in the LOSP flanking and intronic regions. Furthermore, we hypothesize that this HGT-driven evolutionary innovation promoted LOSP function in reproductive isolation, which might be one of the factors determining the intensive cladogenesis in the Littorina (Neritrema) lineage in contrast to the anagenesis in the L. littorea clade.

On the Origin and Evolution of Drosophila New Genes during Spermatogenesis

The origin of functional new genes is a basic biological process that has significant contribution to organismal diversity. Previous studies in both Drosophila and mammals showed that new genes tend to be expressed in testes and avoid the X chromosome, presumably because of meiotic sex chromosome inactivation (MSCI). Here, we analyze the published single-cell transcriptome data of Drosophila adult testis and find an enrichment of male germline mitotic genes, but an underrepresentation of meiotic genes on the X chromosome. This can be attributed to an excess of autosomal meiotic genes that were derived from their X-linked mitotic progenitors, which provides direct cell-level evidence for MSCI in Drosophila. We reveal that new genes, particularly those produced by retrotransposition, tend to exhibit an expression shift toward late spermatogenesis compared with their parental copies, probably due to the more intensive sperm competition or sexual conflict. Our results dissect the complex factors including age, the origination mechanisms and the chromosomal locations that influence the new gene origination and evolution in testes, and identify new gene cases that show divergent cell-level expression patterns from their progenitors for future functional studies.

The Physics-Biology continuum challenges darwinism: Evolution is directed by the homeostasis-dependent bidirectional relation between genome and phenotype

The physics-biology continuum relies on the fact that life emerged from prebiotic molecules. Here, I argue that life emerged from the coupling between nucleic acid and protein synthesis during which proteins (or proto-phenotypes) maintained the physicochemical parameter equilibria (or proto-homeostasis) in the proximity of their encoding nucleic acids (or proto-genomes). This protected the proto-genome physicochemical integrity (i.e., atomic composition) from environmental physicochemical constraints, and therefore increased the probability of reproducing the proto-genome without variation. From there, genomes evolved depending on the biological activities they generated in response to environmental fluctuations. Thus, a genome maintaining homeostasis (i.e., internal physicochemical parameter equilibria), despite and in response to environmental fluctuations, maintains its physicochemical integrity and has therefore a higher probability to be reproduced without variation. Consequently, descendants have a higher probability to share the same phenotype than their parents. Otherwise, the genome is modified during replication as a consequence of the imbalance of the internal physicochemical parameters it generates, until new mutation-deriving biological activities maintain homeostasis in offspring. In summary, evolution depends on feedforward and feedback loops between genome and phenotype, as the internal physicochemical conditions that a genome generates ─ through its derived phenotype in response to environmental fluctuations ─ in turn either guarantee its stability or direct its variation. Evolution may not be explained by the Darwinism-derived, unidirectional principle (random mutations-phenotypes-natural selection) but rather by the bidirectional relationship between genome and phenotype, in which the phenotype in interaction with the environment directs the evolution of the genome it derives from.

Transcription Factors Drive Opposite Relationships between Gene Age and Tissue Specificity in Male and Female Drosophila Gonads

Evolutionarily young genes are usually preferentially expressed in the testis across species. Although it is known that older genes are generally more broadly expressed than younger genes, the properties that shaped this pattern are unknown. Older genes may gain expression across other tissues uniformly, or faster in certain tissues than others. Using Drosophila gene expression data, we confirmed previous findings that younger genes are disproportionately testis biased and older genes are disproportionately ovary biased. We found that the relationship between gene age and expression is stronger in the ovary than any other tissue and weakest in testis. We performed ATAC-seq on Drosophila testis and found that although genes of all ages are more likely to have open promoter chromatin in testis than in ovary, promoter chromatin alone does not explain the ovary bias of older genes. Instead, we found that upstream transcription factor (TF) expression is highly predictive of gene expression in ovary but not in testis. In the ovary, TF expression is more predictive of gene expression than open promoter chromatin, whereas testis gene expression is similarly influenced by both TF expression and open promoter chromatin. We propose that the testis is uniquely able to express younger genes controlled by relatively few TFs, whereas older genes with more TF partners are broadly expressed with peak expression most likely in the ovary. The testis allows widespread baseline expression that is relatively unresponsive to regulatory changes, whereas the ovary transcriptome is more responsive to trans-regulation and has a higher ceiling for gene expression.

Transcriptome Analysis in Yeast Reveals the Externality of Position Effects

The activity of a gene newly integrated into a chromosome depends on the genomic context of the integration site. This “position effect” has been widely reported, although the other side of the coin, that is, how integration affects the local chromosomal environment, has remained largely unexplored, as have the mechanism and phenotypic consequences of this “externality” of the position effect. Here, we examined the transcriptome profiles of approximately 250 Saccharomyces cerevisiae strains, each with GFP integrated into a different locus of the wild-type strain. We found that in genomic regions enriched in essential genes, GFP expression tended to be lower, and the genes near the integration site tended to show greater expression reduction. Further joint analysis with public genome-wide histone modification profiles indicated that this effect was associated with H3K4me2. More importantly, we found that changes in the expression of neighboring genes, but not GFP expression, significantly altered the cellular growth rate. As a result, genomic loci that showed high GFP expression immediately after integration were associated with growth disadvantages caused by elevated expression of neighboring genes, ultimately leading to a low total yield of GFP in the long run. Our results were consistent with competition for transcriptional resources among neighboring genes and revealed a previously unappreciated facet of position effects. This study highlights the impact of position effects on the fate of exogenous gene integration and has significant implications for biological engineering and the pathology of viral integration into the host genome.

Eutherian-Specific Gene TRIML2 Attenuates Inflammation in the Evolution of Placentation

Evolution of highly invasive placentation in the stem lineage of eutherians and subsequent extension of pregnancy set eutherians apart from other mammals, that is, marsupials with short-lived placentas, and oviparous monotremes. Recent studies suggest that eutherian implantation evolved from marsupial attachment reaction, an inflammatory process induced by the direct contact of fetal placenta with maternal endometrium after the breakdown of the shell coat, and shortly before the onset of parturition. Unique to eutherians, a dramatic downregulation of inflammation after implantation prevents the onset of premature parturition, and is critical for the maintenance of gestation. This downregulation likely involved evolutionary changes on maternal as well as fetal/placental side. Tripartite-motif family-like2 (TRIML2) only exists in eutherian genomes and shows preferential expression in preimplantation embryos, and trophoblast-derived structures, such as chorion and placental disc. Comparative genomic evidence supports that TRIML2 originated from a gene duplication event in the stem lineage of Eutheria that also gave rise to eutherian TRIML1. Compared with TRIML1, TRIML2 lost the catalytic RING domain of E3 ligase. However, only TRIML2 is induced in human choriocarcinoma cell line JEG3 with poly(I:C) treatment to simulate inflammation during viral infection. Its knockdown increases the production of proinflammatory cytokines and reduces trophoblast survival during poly(I:C) stimulation, while its overexpression reduces proinflammatory cytokine production, supporting TRIML2’s role as a regulatory inhibitor of the inflammatory pathways in trophoblasts. TRIML2’s potential virus-interacting PRY/SPRY domain shows significant signature of selection, suggesting its contribution to the evolution of eutherian-specific inflammation regulation during placentation.