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.     

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.     

Duplicated genes evolve slower than singletons despite the initial rate increase

Background

Gene duplication is an important mechanism that can lead to the emergence of new functions during evolution. The impact of duplication on the mode of gene evolution has been the subject of several theoretical and empirical comparative-genomic studies. It has been shown that, shortly after the duplication, genes seem to experience a considerable relaxation of purifying selection.

Results

Here we demonstrate two opposite effects of gene duplication on evolutionary rates. Sequence comparisons between paralogs show that, in accord with previous observations, a substantial acceleration in the evolution of paralogs occurs after duplication, presumably due to relaxation of purifying selection. The effect of gene duplication on evolutionary rate was also assessed by sequence comparison between orthologs that have paralogs (duplicates) and those that do not (singletons). It is shown that, in eukaryotes, duplicates, on average, evolve significantly slower than singletons. Eukaryotic ortholog evolutionary rates for duplicates are also negatively correlated with the number of paralogs per gene and the strength of selection between paralogs. A tally of annotated gene functions shows that duplicates tend to be enriched for proteins with known functions, particularly those involved in signaling and related cellular processes; by contrast, singletons include an over-abundance of poorly characterized proteins.

Conclusions

These results suggest that whether or not a gene duplicate is retained by selection depends critically on the pre-existing functional utility of the protein encoded by the ancestral singleton. Duplicates of genes of a higher biological import, which are subject to strong functional constraints on the sequence, are retained relatively more often. Thus, the evolutionary trajectory of duplicated genes appears to be determined by two opposing trends, namely, the post-duplication rate acceleration and the generally slow evolutionary rate owing to the high level of functional constraints.

     

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.     

Correlated asymmetry of sequence and functional divergence between duplicate proteins of Saccharomyces cerevisiae

The role of sequence divergence in functional divergence of duplicate genes is a topic of great interest. In this study, we compare the numbers of amino acid substitutions in each sequence since two yeast duplicates diverged, using a preduplication ancestral outgroup. Using this strategy, we explored the relationship between sequence divergence and functional divergence between duplicate partners. We show that the degree of relative functional asymmetry between duplicate proteins is proportional to the relative sequence divergence between them. Furthermore, of the two duplicates, the copy closer to their ancestral sequence (fewer number of amino acid substitutions) interacts with more proteins and affects fitness more severely when deleted. Therefore, asymmetric sequence divergence between duplicates is correlated with asymmetric functional divergence and may underlie the duplicate's role in genetic robustness against mutations. Among the functional traits considered, protein abundance appears to have the strongest correlation with the nonsynonymous divergence between duplicates. Taken together with the results from whole-genome analyses, our results indicate that within-species duplicates are subject to the same evolutionary force that acts on interspecific sequence and functional divergence. In particular, we detect signs of purifying selection on the more slowly evolving duplicate.     

Evolutionary mechanisms underlying the functional divergence of duplicate genes involved in vertebrates' circadian rhythm pathway

Circadian rhythms, that are governed physiologically and behaviorally by endogenous clock, have been described in many species. Living organisms use this endogenous circadian clock to anticipate environmental transitions, perform activities at biologically advantageous times during the day, and undergo characteristic seasonal responses. Gene duplication is one of the most important mechanisms in the evolution of gene diversity. After duplication, one or both of duplicates can accumulate amino acid changes, thereby promoting functional divergence through the action of natural selection. The circadian system, like many other multigene families, has undergone this genetic revolution, and so circadian genes that are found in single copies in insects are duplicated in vertebrates. We analyzed six groups of genes involved in vertebrates' circadian rhythm pathway to find signatures of molecular evolutionary processes such as gene duplication, natural selection, recombination, and functional divergence. The obtained results, then, were used to determine what evolutionary forces have influenced the fates of duplicated genes of each group. We showed in this research that recombination has not been widespread during the evolution of circadian genes and that purifying selection has been the prominent natural pressure operating on circadian genes. We also showed that the evolution of circadian genes has been depended on gene duplication and functional divergence. Finally, we put forward models best describing the evolutionary fates of circadian duplicates.     

The Ghost of Selection Past: Rates of Evolution and Functional Divergence of Anciently Duplicated Genes

The duplication of genes and even complete genomes may be a prerequisite for major evolutionary transitions and the origin of evolutionary novelties. However, the evolutionary mechanisms of gene evolution and the origin of novel gene functions after gene duplication have been a subject of many debates. Recently, we compiled 26 groups of orthologous genes, which included one gene from human, mouse, and chicken, one or two genes from the tetraploid Xenopus and two genes from zebrafish. Comparative analysis and mapping data showed that these pairs of zebrafish genes were probably produced during a fish-specific genome duplication that occurred between 300 and 450 Mya, before the teleost radiation (Taylor et al. 2001). As discussed here, many of these retained duplicated genes code for DNA binding proteins. Different models have been developed to explain the retention of duplicated genes and in particular the subfunctionalization model of Force et al. (1999) could explain why so many developmental control genes have been retained. Other models are harder to reconcile with this particular set of duplicated genes. Most genes seem to have been subjected to strong purifying selection, keeping properties such as charge and polarity the same in both duplicates, although some evidence was found for positive Darwinian selection, in particular for Hox genes. However, since only the cumulative pattern of nucleotide substitutions can be studied, clear indications of positive Darwinian selection or neutrality may be hard to find for such anciently duplicated genes. Nevertheless, an increase in evolutionary rate in about half of the duplicated genes seems to suggest that either positive Darwinian selection has occurred or that functional constraints have been relaxed at one point in time during functional divergence. Received: 4 January 2001 / Accepted: 29 March 2001     

Genomic background predicts the fate of duplicated genes: evidence from the yeast genome

Gene duplication with subsequent divergence plays a central role in the acquisition of genes with novel function and complexity during the course of evolution. With reduced functional constraints or through positive selection, these duplicated genes may experience accelerated evolution. Under the model of subfunctionalization, loss of subfunctions leads to complementary acceleration at sites with two copies, and the difference in average rate between the sequences may not be obvious. On the other hand, the classical model of neofunctionalization predicts that the evolutionary rate in one of the two duplicates is accelerated. However, the classical model does not tell which of the duplicates experiences the acceleration in evolutionary rate. Here, we present evidence from the Saccharomyces cerevisiae genome that a duplicate located in a genomic region with a low-recombination rate is likely to evolve faster than a duplicate in an area of high recombination. This observation is consistent with population genetics theory that predicts that purifying selection is less effective in genomic regions of low recombination (Hill-Robertson effect). Together with previous studies, our results suggest the genomic background (e.g., local recombination rate) as a potential force to drive the divergence between nontandemly duplicated genes. This implies the importance of structure and complexity of genomes in the diversification of organisms via gene duplications.     

Importance of lineage-specific expansion of plant tandem duplicates in the adaptive response to environmental stimuli

Plants have substantially higher gene duplication rates compared with most other eukaryotes. These plant gene duplicates are mostly derived from whole genome and/or tandem duplications. Earlier studies have shown that a large number of duplicate genes are retained over a long evolutionary time, and there is a clear functional bias in retention. However, the influence of duplication mechanism, particularly tandem duplication, on duplicate retention has not been thoroughly investigated. We have defined orthologous groups (OGs) between Arabidopsis (Arabidopsis thaliana) and three other land plants to examine the functional bias of retained duplicate genes during vascular plant evolution. Based on analysis of Gene Ontology categories, it is clear that genes in OGs that expanded via tandem duplication tend to be involved in responses to environmental stimuli, while those that expanded via nontandem mechanisms tend to have intracellular regulatory roles. Using Arabidopsis stress expression data, we further demonstrated that tandem duplicates in expanded OGs are significantly enriched in genes that are up-regulated by biotic stress conditions. In addition, tandem duplication of genes in an OG tends to be highly asymmetric. That is, expansion of OGs with tandem genes in one organismal lineage tends to be coupled with losses in the other. This is consistent with the notion that these tandem genes have experienced lineage-specific selection. In contrast, OGs with genes duplicated via nontandem mechanisms tend to experience convergent expansion, in which similar numbers of genes are gained in parallel. Our study demonstrates that the expansion of gene families and the retention of duplicates in plants exhibit substantial functional biases that are strongly influenced by the mechanism of duplication. In particular, genes involved in stress responses have an elevated probability of retention in a single-lineage fashion following tandem duplication, suggesting that these tandem duplicates are likely important for adaptive evolution to rapidly changing 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.     

Positive Selection and Functional Divergence After Melanopsin Gene Duplication

A newly discovered melanopsin gene (Opn4) encodes a member of the opsins, melanopsin. Two melanopsin genes, mammalian-like Opn4m and Xenopus-like Opn4x, have been described in nonmammalian vertebrates, but the underlying evolutionary mechanisms behind the duplication of melanopsin genes remain unclear. We conducted a comprehensive evolutionary analysis within a phylogenetic framework. In our phylogenetic tree, the duplication of Opn4m and Opn4x probably occurred prior to the emergence of vertebrates, and subsequently Opn4x disappeared in the lineages leading to mammalian species. Evolutionary analyses show strong purifying selection during melanopsin evolution. We also provide evidence that Opn4x underwent positive selection after the early gene duplication events. It has been indicated that functional divergence and altered functional constraints occurred between Opn4m and Opn4x duplicates with the identification of positively selected amino acids. Our findings highlight the evolutionary malleability in vertebrate melanopsin genes and provide a genetic basis for comparative studies of functional properties of these two melanopsins.     

Mouse duplicate genes are as essential as singletons

Duplicate genes in mouse are widely thought to have functional redundancy, and to be less essential than singleton genes. We analyzed nearly 3900 individually knocked out mouse genes and discovered that the proportion of essential genes is approximately 55% in both singletons and duplicates. This suggests that mammalian duplicates rarely compensate for each other, and that the absence of phenotypes in mice deficient for a duplicate gene should not be automatically attributed to paralogous compensation.     

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.     

Evolution of cis-regulatory elements in duplicated genes of yeast

An increasing number of studies report that functional divergence in duplicated genes is accompanied by gene expression changes, although the evolutionary mechanism behind this process remains unclear. Our genomic analysis on the yeast Saccharomyces cerevisiae shows that the number of shared regulatory motifs in the duplicates decreases with evolutionary time, whereas the total number of regulatory motifs remains unchanged. Moreover, genes with numerous paralogs in the yeast genome do not have especially low number of regulatory motifs. These findings indicate that degenerative complementation is not the sole mechanism behind expression divergence in yeast. Moreover, we found some evidence for the action of positive selection on cis-regulatory motifs after gene duplication. These results suggest that the evolution of functional novelty has a substantial role in yeast duplicate gene evolution.     

Relative evolutionary rates of NBS-encoding genes revealed by soybean segmental duplication

It is well known that nucleotide binding site (NBS)-encoding genes are duplicate-rich and fast-evolving genes. However, there is little information on the relative importance of tandem and segmental NBS duplicates and their exact evolutionary rates. The two rounds of large-scale duplication that have occurred in soybean provide a unique opportunity to investigate these issues. Comparison of NBS and non-NBS genes on segments of syntenic homoeologs shows that NBS-encoding genes evolve at least 1.5-fold faster (~1.5-fold higher synonymous and approximately 2.3-fold higher nonsynonymous substitution rates) and lose their genes approximately twofold faster than the flanking non-NBS genes. Compared with segmental duplicates, tandem NBS duplicates are more abundant in soybean, suggesting that tandem duplication is the major driving force in the expansion of NBS genes. Notably, significant sequence exchanges along with significantly positive selection were detected in most tandem-duplicated NBS gene families. The results suggest that the rapid evolution of NBS genes may be due to the combined effects of diversifying selection and frequent sequence exchanges. Interestingly, TIR-NBS-LRR genes (TNLs) have a higher nucleotide substitution rate than non-TNLs, indicating that these types of NBS genes may have a rather different evolutionary pattern. It is important to determine the exact relative evolutionary rates of TNL, non-TNL, and non-NBS genes in order to understand how fast the host plant can adjust its response to rapidly evolving pathogens in a coevolutionary context.     

Comparative analysis of teleost fish genomes reveals preservation of different ancient clock duplicates in different fishes

Clock (Circadian locomotor output cycle kaput) was the first vertebrate circadian clock gene identified in a mouse forward genetics mutagenesis screen. It encodes a bHLH-PAS protein that is highly conserved throughout evolution. Tetrapods also have the second Clock gene, Clock2 or Npas2 (Neuronal PAS domain protein 2). Conversely, the fruit fly, an invertebrate, has only one clock gene. Interrogation of the five teleost fish genome databases revealed that the zebrafish and the Japanese pufferfish (fugu) each have three clock genes, whereas the green spotted pufferfish (tetraodon), the Japanese medaka fish and the three-spine stickleback each have two clock genes. Phylogenetic and splice site analyses indicated that zebrafish and fugu each have two clock1 genes, clock1a and clock1b and one clock2; tetraodon also have clock1a and clock1b but do not have clock2; and medaka and stickleback each have clock1b and one clock2. Genome neighborhood analysis further showed that clock1a/clock1b in zebrafish, fugu and tetraodon is an ancient duplicate. While the dN/dS ratios of these three fish clock duplicates are all <1, indicating that purifying selection has acted upon them; the Tajima relative rate test showed that all three fish clock duplicates have asymmetric evolutionary rates, implicating that one of these duplicates have been under positive selection or relaxed functional constraint. These results support the view that teleost fish clock genes were generated from an ancient genome-wide duplication, and differential gene loss after the duplication resulted in retention of different ancient duplicates in different teleost fishes, which could have contributed to the evolution of the distinct fish circadian clock mechanisms.     

Functional analysis of gene duplications in Saccharomyces cerevisiae

Gene duplication can occur on two scales: whole-genome duplications (WGD) and smaller-scale duplications (SSD) involving individual genes or genomic segments. Duplication may result in functionally redundant genes or diverge in function through neofunctionalization or subfunctionalization. The effect of duplication scale on functional evolution has not yet been explored, probably due to the lack of global knowledge of protein function and different times of duplication events. To address this question, we used integrated Bayesian analysis of diverse functional genomic data to accurately evaluate the extent of functional similarity and divergence between paralogs on a global scale. We found that paralogs resulting from the whole-genome duplication are more likely to share interaction partners and biological functions than smaller-scale duplicates, independent of sequence similarity. In addition, WGD paralogs show lower frequency of essential genes and higher synthetic lethality rate, but instead diverge more in expression pattern and upstream regulatory region. Thus, our analysis demonstrates that WGD paralogs generally have similar compensatory functions but diverging expression patterns, suggesting a potential of distinct evolutionary scenarios for paralogs that arose through different duplication mechanisms. Furthermore, by identifying these functional disparities between the two types of duplicates, we reconcile previous disputes on the relationship between sequence divergence and expression divergence or essentiality.