A Population-Genetic Lens into the Process of Gene Loss Following Whole-Genome Duplication

Whole-genome duplications (WGDs) have occurred in many eukaryotic lineages. However, the underlying evolutionary forces and molecular mechanisms responsible for the long-term retention of gene duplicates created by WGDs are not well understood. We employ a population-genomic approach to understand the selective forces acting on paralogs and investigate ongoing duplicate-gene loss in multiple species of Paramecium that share an ancient WGD. We show that mutations that abolish protein function are more likely to be segregating in retained WGD paralogs than in single-copy genes, most likely because of ongoing nonfunctionalization post-WGD. This relaxation of purifying selection occurs in only one WGD paralog, accompanied by the gradual fixation of nonsynonymous mutations and reduction in levels of expression, and occurs over a long period of evolutionary time, “marking” one locus for future loss. Concordantly, the fitness effects of new nonsynonymous mutations and frameshift-causing indels are significantly more deleterious in the highly expressed copy compared with their paralogs with lower expression. Our results provide a novel mechanistic model of gene duplicate loss following WGDs, wherein selection acts on the sum of functional activity of both duplicate genes, allowing the two to wander in expression and functional space, until one duplicate locus eventually degenerates enough in functional efficiency or expression that its contribution to total activity is too insignificant to be retained by purifying selection. Retention of duplicates by such mechanisms predicts long times to duplicate-gene loss, which should not be falsely attributed to retention due to gain/change in function.     

Preferential retention of the slowly evolving gene in pairs of duplicates in angiosperm genomes

Gene duplication provides raw material for functional innovation, but gene duplicability varies considerably. Previous studies have found widespread asymmetrical sequence evolution between paralogs. However, it remains unknown whether the rate of evolution among paralogs affects their propensity of being retained after another round of whole-genome duplication (WGD). In this study, we investigated gene groups that have experienced two successive WGDs to determine which of two older duplicates with different evolutionary rates was more likely to retain both younger duplicates. To uncouple the measurement of evolutionary rates from any assignment of duplicate or singleton status, we measured the evolutionary rates of singleton genes in out-lineages but classified these singleton genes according to whether they are retained or not in a crown group of species. We found that genes that retained younger duplicates in the crown group of genomes were more constrained prior to the younger duplication event than those that failed to leave duplicates. In addition, we also found that the retained clades have more genes in out-lineages. Subsequent analyses showed that genes in the retained clades were expressed more broadly and highly than genes in the singleton clades. We concluded that the set of repeatedly retained genes after two WGDs is biased toward slowly evolving genes in angiosperms, suggesting that the potential of genes for both functional conservation and divergence likely affects their propensity of being retained after WGD in angiosperms.     

New genes from non-coding sequence: the role of de novo protein-coding genes in eukaryotic evolutionary innovation

The origin of novel protein-coding genes de novo was once considered so improbable as to be impossible. In less than a decade, and especially in the last five years, this view has been overturned by extensive evidence from diverse eukaryotic lineages. There is now evidence that this mechanism has contributed a significant number of genes to genomes of organisms as diverse as Saccharomyces, Drosophila, Plasmodium, Arabidopisis and human. From simple beginnings, these genes have in some instances acquired complex structure, regulated expression and important functional roles. New genes are often thought of as dispensable late additions; however, some recent de novo genes in human can play a role in disease. Rather than an extremely rare occurrence, it is now evident that there is a relatively constant trickle of proto-genes released into the testing ground of natural selection. It is currently unknown whether de novo genes arise primarily through an ‘RNA-first’ or ‘ORF-first’ pathway. Either way, evolutionary tinkering with this pool of genetic potential may have been a significant player in the origins of lineage-specific traits and adaptations.     

Drosophila duplicate genes evolve new functions on the fly

Gene duplication is thought to play a key role in phenotypic innovation. While several processes have been hypothesized to drive the retention and functional evolution of duplicate genes, their genomic contributions have never been determined. We recently developed the first genome-wide method to classify these processes by comparing distances between expression profiles of duplicate genes and their ancestral single-copy orthologs. Application of our approach to spatial gene expression profiles in two Drosophila species revealed that a majority of young duplicate genes possess new functions, and that new functions are acquired rapidly—often within a few million years. Surprisingly, new functions tend to arise in younger copies of duplicate gene pairs. Moreover, we found that young duplicates are often specifically expressed in testes, whereas old duplicates are broadly expressed across several tissues, providing strong support for the hypothetical “out-of-testes” origin of new genes. In this Extra View, I discuss our findings in the context of theoretical predictions about gene duplication, with a particular emphasis on the importance of natural selection in the evolution of novel phenotypes.     

Sorting duplicated loci disentangles complexities of polyploid genomes masked by genotyping by sequencing

Many plants and animals of polyploid origin are currently enjoying a genomics explosion enabled by modern sequencing and genotyping technologies. However, routine filtering of duplicated loci in most studies using genotyping by sequencing introduces an unacceptable, but often overlooked, bias when detecting selection. Retained duplicates from ancient whole‐genome duplications (WGDs) may be found throughout genomes, whereas retained duplicates from recent WGDs are concentrated at distal ends of some chromosome arms. Additionally, segmental duplicates can be found at distal ends or nearly anywhere in a genome. Evidence shows that these duplications facilitate adaptation through one of two pathways: neo‐functionalization or increased gene expression. Filtering duplicates removes distal ends of some chromosomes, and distal ends are especially known to harbour adaptively important genes. Thus, filtering of duplicated loci impoverishes the interpretation of genomic data as signals from contiguous duplicated genes are ignored. We review existing strategies to genotype and map duplicated loci; we focus in detail on an overlooked strategy of using gynogenetic haploids (1N) as a part of new genotyping by sequencing studies. We provide guidelines on how to use this haploid strategy for studies on polyploid‐origin vertebrates including how it can be used to screen duplicated loci in natural populations. We conclude by discussing areas of research that will benefit from better inclusion of polyploid loci; we particularly stress the sometimes overlooked fact that basing genomic studies on dense maps provides value added in the form of locating and annotating outlier loci or colocating outliers into islands of divergence.     

Degree of Functional Divergence in Duplicates Is Associated with Distinct Roles in Plant Evolution

Gene duplication is a major mechanism to create new genes. After gene duplication, some duplicated genes undergo functionalization, whereas others largely maintain redundant functions. Duplicated genes comprise various degrees of functional diversification in plants. However, the evolutionary fate of high and low diversified duplicates is unclear at genomic scale. To infer high and low diversified duplicates in Arabidopsis thaliana genome, we generated a prediction method for predicting whether a pair of duplicate genes was subjected to high or low diversification based on the phenotypes of knock-out mutants. Among 4,017 pairs of recently duplicated A. thaliana genes, 1,052 and 600 are high and low diversified duplicate pairs, respectively. The predictions were validated based on the phenotypes of generated knock-down transgenic plants. We determined that the high diversified duplicates resulting from tandem duplications tend to have lineage-specific functions, whereas the low diversified duplicates produced by whole-genome duplications are related to essential signaling pathways. To assess the evolutionary impact of high and low diversified duplicates in closely related species, we compared the retention rates and selection pressures on the orthologs of A. thaliana duplicates in two closely related species. Interestingly, high diversified duplicates resulting from tandem duplications tend to be retained in multiple lineages under positive selection. Low diversified duplicates by whole-genome duplications tend to be retained in multiple lineages under purifying selection. Taken together, the functional diversities determined by different duplication mechanisms had distinct effects on plant evolution.     

Functional divergence of gene duplicates through ectopic recombination

Gene duplication stimulates evolutionary innovation as the resulting paralogs acquire mutations that lead to sub‐ or neofunctionalization. A comprehensive in silico analysis of paralogs in Saccharomyces cerevisiae reveals that duplicates of cell‐surface and subtelomeric genes also undergo ectopic recombination, which leads to new chimaeric alleles. Mimicking such intergenic recombination events in the FLO (flocculation) family of cell‐surface genes shows that chimaeric FLO alleles confer different adhesion phenotypes than the parental genes. Our results indicate that intergenic recombination between paralogs can generate a large set of new alleles, thereby providing the raw material for evolutionary adaptation and innovation.     

De novo origin of new genes with introns in Plasmodium vivax

The origin of new genes is critical for organisms adapting to new niches. Here, we present evidence for a recent de novo origin of at least 13 protein-coding genes in the genome of Plasmodium vivax. Although recently de novo originated genes have often been suggested to be initially intronless, five of the genes identified in our analysis contain introns in their coding regions. Further investigations revealed that these introns likely evolved from previously intergenic regions together with the coding sequences. We discuss the potential mechanisms for intron formation in these genes and propose that intronization be considered in the formation of de novo originated genes.     

Hundreds of Out-of-Frame Remodeled Gene Families in the Escherichia coli Pangenome

All genomes include gene families with very limited taxonomic distributions that potentially represent new genes and innovations in protein-coding sequence, raising questions on the origins of such genes. Some of these genes are hypothesized to have formed de novo, from noncoding sequences, and recent work has begun to elucidate the processes by which de novo gene formation can occur. A special case of de novo gene formation, overprinting, describes the origin of new genes from noncoding alternative reading frames of existing open reading frames (ORFs). We argue that additionally, out-of-frame gene fission/fusion events of alternative reading frames of ORFs and out-of-frame lateral gene transfers could contribute to the origin of new gene families. To demonstrate this, we developed an original pattern-search in sequence similarity networks, enhancing the use of these graphs, commonly used to detect in-frame remodeled genes. We applied this approach to gene families in 524 complete genomes of Escherichia coli. We identified 767 gene families whose evolutionary history likely included at least one out-of-frame remodeling event. These genes with out-of-frame components represent ∼2.5% of all genes in the E. coli pangenome, suggesting that alternative reading frames of existing ORFs can contribute to a significant proportion of de novo genes in bacteria.     

Evolution of Stress-Regulated Gene Expression in Duplicate Genes of Arabidopsis thaliana

Due to the selection pressure imposed by highly variable environmental conditions, stress sensing and regulatory response mechanisms in plants are expected to evolve rapidly. One potential source of innovation in plant stress response mechanisms is gene duplication. In this study, we examined the evolution of stress-regulated gene expression among duplicated genes in the model plant Arabidopsis thaliana. Key to this analysis was reconstructing the putative ancestral stress regulation pattern. By comparing the expression patterns of duplicated genes with the patterns of their ancestors, duplicated genes likely lost and gained stress responses at a rapid rate initially, but the rate is close to zero when the synonymous substitution rate (a proxy for time) is >~0.8. When considering duplicated gene pairs, we found that partitioning of putative ancestral stress responses occurred more frequently compared to cases of parallel retention and loss. Furthermore, the pattern of stress response partitioning was extremely asymmetric. An analysis of putative cis-acting DNA regulatory elements in the promoters of the duplicated stress-regulated genes indicated that the asymmetric partitioning of ancestral stress responses are likely due, at least in part, to differential loss of DNA regulatory elements; the duplicated genes losing most of their stress responses were those that had lost more of the putative cis-acting elements. Finally, duplicate genes that lost most or all of the ancestral responses are more likely to have gained responses to other stresses. Therefore, the retention of duplicates that inherit few or no functions seems to be coupled to neofunctionalization. Taken together, our findings provide new insight into the patterns of evolutionary changes in gene stress responses after duplication and lay the foundation for testing the adaptive significance of stress regulatory changes under highly variable biotic and abiotic environments.     

Asymmetric evolution of duplicate genes encoding the CCAAT-binding factor NF-Y in plant genomes

NF-Y is a ubiquitous CCAAT-binding factor composed of NF-YA, NF-YB and NF-YC. Multiple genes encoding NF-Y subunits have been identified in plant genomes. It remains unclear whether the duplicate genes underwent different evolutionary patterns. Likelihood-ratio tests were used to examine whether the amino acid substitution rates are the same between duplicate genes. The influences of selection on evolution were evaluated by comparing the conservative and radical amino acid substitution rates, as well as maximum-likelihood analysis. Some NF-YB and NF-YC duplicates showed significant evidence of asymmetric evolution but not the NF-YA duplicates. Most amino acid replacements in the NF-YB and NF-YC duplicates result in changes in hydropathy, polar requirement and polarity. The physicochemical changes in the sequences of NF-YB seem to be coupled to asymmetric divergence in gene function. Plant NF-Y genes have evolved in different patterns. Relaxed selective constraints following gene duplication are most likely responsible for the unequal evolutionary rates and distinct divergence patterns of duplicate NF-Y genes. Positive selection may have promoted amino acid hydropathy changes in the NF-YC duplicates.     

The Roles of Whole-Genome and Small-Scale Duplications in the Functional Specialization of Saccharomyces cerevisiae Genes

Researchers have long been enthralled with the idea that gene duplication can generate novel functions, crediting this process with great evolutionary importance. Empirical data shows that whole-genome duplications (WGDs) are more likely to be retained than small-scale duplications (SSDs), though their relative contribution to the functional fate of duplicates remains unexplored. Using the map of genetic interactions and the re-sequencing of 27 Saccharomyces cerevisiae genomes evolving for 2,200 generations we show that SSD-duplicates lead to neo-functionalization while WGD-duplicates partition ancestral functions. This conclusion is supported by: (a) SSD-duplicates establish more genetic interactions than singletons and WGD-duplicates; (b) SSD-duplicates copies share more interaction-partners than WGD-duplicates copies; (c) WGD-duplicates interaction partners are more functionally related than SSD-duplicates partners; (d) SSD-duplicates gene copies are more functionally divergent from one another, while keeping more overlapping functions, and diverge in their sub-cellular locations more than WGD-duplicates copies; and (e) SSD-duplicates complement their functions to a greater extent than WGD–duplicates. We propose a novel model that uncovers the complexity of evolution after gene duplication.     

DNA methylation profiling of pseudogene-parental gene pairs and two gene families

A substantial proportion of human genes contain tissue-specifically DNA-methylated regions (TDMRs). However, little is known about the evolutionary conservation of differentially methylated loci, how they evolve, and the signals that regulate them. We have studied TDMR conservation in the PLG and TBX gene families and in 32 pseudogene–parental gene pairs. Among the members of the recently evolved PLG gene family, 5′-UTR methylation is conserved and inversely correlated with the cognate gene expression, indicating as well a conserved regulatory role of DNA methylation. Conversely, many genes of the much older TBX family display complementary tissue-specific methylation, suggesting an epigenetic complementation in the evolution of this gene family. Similar to gene families, unprocessed pseudogenes arose from gene duplications and we found TDMR conservation in some pseudogene–parental gene pairs displaying short evolutionary distances. However, for the majority of unprocessed pseudogenes and for all processed pseudogenes examined, we found that tissue-specific methylation arose de novo after gene duplication.     

An ancient frame-shifting event in the highly conserved KPNA gene family has undergone extensive compensation by natural selection in vertebrates

One of the prevailing arguments in evolutionary theory is that the duplicates of genes can acquire novel functionality. This is because only one of the paralogs need maintain the ancestral function, leaving room for natural experimentation due to a respite in purifying selection. Although many duplicates can subsequently become disabled by nullifying mutations, a few may also go on to diverge along a novel evolutionary trajectory. Here, evidence is provided that demonstrates how this scenario may not always be true. Rather, in the case of the highly conserved KPNA importin family, an initial relaxation in selection induced a frameshift that was later suppressed and heavily compensated for as part of a reparative and optimizing process. Despite a resulting divergence, there remains a distinct preservation of both sequence and functionality among the paralogs. This would indicate that duplicates can be retained by selection for reasons related to their redundant functionality. It also shows that, even when positive selection is inferred in duplicate genes, this may be of a compensatory nature rather than one representing any biochemical innovation. Generally, this development would perhaps be a more common outcome for gene duplication than is currently maintained.     

On the origins of Mendelian disease genes in man: the impact of gene duplication

Over 3,000 human diseases are known to be linked to heritable genetic variation, mapping to over 1,700 unique genes. Dating of the evolutionary age of these disease-associated genes has suggested that they have a tendency to be ancient, specifically coming into existence with early metazoa. The approach taken by past studies, however, assumes that the age of a disease is the same as the age of its common ancestor, ignoring the fundamental contribution of duplication events in the evolution of new genes and function. Here, we date both the common ancestor and the duplication history of known human disease-associated genes. We find that the majority of disease genes (80%) are genes that have been duplicated in their evolutionary history. Periods for which there are more disease-associated genes, for example, at the origins of bony vertebrates, are explained by the emergence of more genes at that time, and the majority of these are duplicates inferred to have arisen by whole-genome duplication. These relationships are similar for different disease types and the disease-associated gene's cellular function. This indicates that the emergence of duplication-associated diseases has been ongoing and approximately constant (relative to the retention of duplicate genes) throughout the evolution of life. This continued until approximately 390 Ma from which time relatively fewer novel genes came into existence on the human lineage, let alone disease genes. For single-copy genes associated with disease, we find that the numbers of disease genes decreases with recency. For the majority of duplicates, the disease-associated mutation is associated with just one of the duplicate copies. A universal explanation for heritable disease is, thus, that it is merely a by-product of the evolutionary process; the evolution of new genes (de novo or by duplication) results in the potential for new diseases to emerge.