Asymmetric functional divergence of duplicate genes in yeast

Most duplicate genes are eliminated from a genome shortly after duplication, but those that remain are an important source of biochemical diversity. Here, I present evidence from genome-scale protein-protein interaction data, microarray expression data, and large-scale gene knockout data that this diversification is often asymmetrical: one duplicate usually shows significantly more molecular or genetic interactions than the other. I propose a model that can explain this divergence pattern if asymmetrically diverging duplicate gene pairs show increased robustness to deleterious mutations.     

The evolutionary demography of duplicate genes

Although gene duplication has generally been viewed as a necessary source of material for the origin of evolutionary novelties, the rates of origin, loss, and preservation of gene duplicates are not well understood. Applying steady-state demographic techniques to the age distributions of duplicate genes censused in seven completely sequenced genomes, we estimate the average rate of duplication of a eukaryotic gene to be on the order of 0.01/gene/million years, which is of the same order of magnitude as the mutation rate per nucleotide site. However, the average half-life of duplicate genes is relatively small, on the order of 4.0 million years. Significant interspecific variation in these rates appears to be responsible for differences in species-specific genome sizes that arise as a consequence of a quasi-equilibrium birth-death process. Most duplicated genes experience a brief period of relaxed selection early in their history and a minority exhibit the signature of directional selection, but those that survive more than a few million years eventually experience strong purifying selection. Thus, although most theoretical work on the gene-duplication process has focused on issues related to adaptive evolution, the origin of a new function appears to be a very rare fate for a duplicate gene. A more significant role of the duplication process may be the generation of microchromosomal rearrangements through reciprocal silencing of alternative copies, which can lead to the passive origin of post-zygotic reproductive barriers in descendant lineages of incipient species.     

External factors accelerate expression divergence between duplicate genes

We examined the evolution of expression of duplicate genes in Arabidopsis thaliana, by analyzing 512 data sets of gene expression microarrays and 2022 recent duplicate gene pairs. Expression divergence between gene duplicates is significantly greater in response to environmental stress than to developmental processes. A slow rate of expression divergence during development might offer dosage-dependent selective advantage, whereas rapid expression divergence in response to external changes might accelerate adaptation.     

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.     

Adaptive functional divergence among triplicated alpha-globin genes in rodents

The functional divergence of duplicated genes is thought to play an important role in the evolution of new developmental and physiological pathways, but the role of positive selection in driving this process remains controversial. The objective of this study was to test whether amino acid differences among triplicated alpha-globin paralogs of the Norway rat (Rattus norvegicus) and the deer mouse (Peromyscus maniculatus) are attributable to a relaxation of purifying selection or to a history of positive selection that has adapted the gene products to new or modified physiological tasks. In each rodent species, the two paralogs at the 5'-end of the alpha-globin gene cluster (HBA-T1 and HBA-T2) are evolving in concert and are therefore identical or nearly identical in sequence. However, in each case, the HBA-T1 and HBA-T2 paralogs are distinguished from the third paralog at the 3'-end of the gene cluster (HBA-T3) by multiple amino acid substitutions. An analysis of genomic sequence data from several rodent species revealed that the HBA-T3 genes of Rattus and Peromyscus originated via independent, lineage-specific duplication events. In the independently derived HBA-T3 genes of both species, a likelihood analysis based on a codon-substitution model revealed that accelerated rates of amino acid substitution are attributable to positive directional selection, not to a relaxation of purifying selection. As a result of functional divergence among the triplicated alpha-globin genes in Rattus and Peromyscus, the red blood cells of both rodent species contain a mixture of functionally distinct alpha-chain hemoglobin isoforms that are predicted to have different oxygen-binding affinities. In P. maniculatus, a species that is able to sustain physiological function under conditions of chronic hypoxia at high altitude, the coexpression of distinct hemoglobin isoforms with graded oxygen affinities is expected to broaden the permissible range of arterial oxygen tensions for pulmonary/tissue oxygen transport.     

Expression diversity and evolutionary dynamics of rice duplicate genes

Duplicate genes are believed to be a major source of new gene functions over evolutionary time. In order to evaluate the evolutionary dynamics of rice duplicate genes, formed principally by paleoployploidization prior to the speciation of the Poaceae family, we have employed a public microarray dataset including 155 gene expression omnibus sample plates and bioinformatics tools. At least 57.4% of old ~70 million years ago (MYA) duplicate gene pairs exhibit divergences in expression over the given experimental set, whereas at least 50.9% of young ~7.7-MYA duplicate gene pairs were shown to be divergent. When grouping the rice duplicate genes according to functional categories, we noted a striking and significant enrichment of divergent duplicate metabolism-associated genes, as compared to that observed in non-divergent duplicate genes. While both non-synonymous substitution (Ka) and synonymous substitution (Ks) values between non- and divergent duplicate gene pairs evidenced significant differences, the Ka/Ks values between them exhibited no significant differences. Interestingly, the average numbers of conserved motifs of the duplicate gene pairs revealed a pattern of decline along with an increase in expression diversity, partially supporting the subfunctionalization model with degenerative complementation in regulatory motifs. Duplicate gene pairs with high local similarity (HLS) segments, which might be formed via conversion between rice paleologs, evidenced higher expression correlations than were observed in the gene pairs without the HLS segments; this probably resulted in an increased likelihood of gene conversion in promoters of the gene pairs harboring HLS segments. More than 60% of the rice gene families exhibited similar high expression diversity between members as compared to that of randomly selected gene pairs. These findings are likely reflective of the evolutionary dynamics of rice duplicate genes for gene retention. An erratum to this article can be found at     

Polymorphism and divergence at three duplicate genes in Brassica nigra

The CONSTANS-like gene family has been shown to evolve exceptionally fast in Brassicaceae. In the present study we analyzed sequence polymorphism and divergence of three genes from this family: COL1 (CONSTANS-LIKE 1) and two copies of CO (CONSTANS), COa and COb, in B. nigra. There was a significant fourfold difference in overall nucleotide diversity among the three genes, with BniCOb having twice as much variation as BniCOL1, which in turn was twice as variable as BniCOa. The ratio of nonsynonymous-to-synonymous substitutions (dN/dS) was high for all three genes, confirming previous studies. While we did not detect evidence of selection at BniCOa and BniCOb, there was a significant excess of polymorphic synonymous mutations in a McDonald-Kreitman test comparing COL1 in B. nigra and A. thaliana. This is apparently the result of an increase in selective constraint on COL1 in B. nigra combined with a decrease in A. thaliana. In conclusion, a complex scenario involving both demography and selection seems to have shaped the pattern of polymorphism at the three genes.     

The Zn-peptidase superfamily: functional convergence after evolutionary divergence.

Zn-dependent carboxypeptidases (ZnCP) cleave off the C-terminal amino acid residues from proteins and peptides. Here we describe a superfamily that unites classical ZnCP with other enzymes, most of which are known (or likely) to participate in metal-dependent peptide bond cleavage, but not necessarily in polypeptide substrates. It is demonstrated that aspartoacylase (ASP gene) and succinylglutamate desuccinylase (ASTE gene) are members of the ZnCP family. The Zn-binding site along with the structural core of the protein is shown to be conserved between ZnCP and another large family of hydrolases that includes mostly aminopeptidases (ZnAP). Both families (ZnCP and ZnAP) include not only proteases but also enzymes that perform N-deacylation, and enzymes that catalyze N-desuccinylation of amino acids. This is a result of functional convergence that apparently occurred after the divergence of the two families.     

Rapid divergence in expression between duplicate genes inferred from microarray data

For more than 30 years, expression divergence has been considered as a major reason for retaining duplicated genes in a genome, but how often and how fast duplicate genes diverge in expression has not been studied at the genomic level. Using yeast microarray data, we show that expression divergence between duplicate genes is significantly correlated with their synonymous divergence (K_S) and also with their nonsynonymous divergence (K_A) if K_A ≤ 0.3. Thus, expression divergence increases with evolutionary time, and K_A is initially coupled with expression divergence. More interestingly, a large proportion of duplicate genes have diverged quickly in expression and the vast majority of gene pairs eventually become divergent in expression. Indeed, more than 40% of gene pairs show expression divergence even when K_S is ≤ 0.10, and this proportion becomes >80% for K_S > 1.5. Only a small fraction of ancient gene pairs do not show expression divergence.     

Predicted functional RNAs within coding regions constrain evolutionary rates of yeast proteins

Functional RNAs (fRNAs) are being recognized as an important regulatory component in biological processes. Interestingly, recent computational studies suggest that the number and biological significance of functional RNAs within coding regions (coding fRNAs) may have been underestimated. We hypothesized that such coding fRNAs will impose additional constraint on sequence evolution because the DNA primary sequence has to simultaneously code for functional RNA secondary structures on the messenger RNA in addition to the amino acid codons for the protein sequence. To test this prediction, we first utilized computational methods to predict conserved fRNA secondary structures within multiple species alignments of Saccharomyces sensu strico genomes. We predict that as much as 5% of the genes in the yeast genome contain at least one functional RNA secondary structure within their protein-coding region. We then analyzed the impact of coding fRNAs on the evolutionary rate of protein-coding genes because a decrease in evolutionary rate implies constraint due to biological functionality. We found that our predicted coding fRNAs have a significant influence on evolutionary rates (especially at synonymous sites), independent of other functional measures. Thus, coding fRNA may play a role on sequence evolution. Given that coding regions of humans and flies contain many more predicted coding fRNAs than yeast, the impact of coding fRNAs on sequence evolution may be substantial in genomes of higher eukaryotes.     

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.     

Different evolutionary patterns between young duplicate genes in the human genome

Background  

Following gene duplication, two duplicate genes may experience relaxed functional constraints or acquire different mutations, and may also diverge in function. Whether the two copies will evolve in different patterns remains unclear, however, because previous studies have reached conflicting conclusions. In order to resolve this issue, by providing a general picture, we studied 250 independent pairs of young duplicate genes from the whole human genome.     

Rapid and asymmetric divergence of duplicate genes in the human gene coexpression network

Background  

While gene duplication is known to be one of the most common mechanisms of genome evolution, the fates of genes after duplication are still being debated. In particular, it is presently unknown whether most duplicate genes preserve (or subdivide) the functions of the parental gene or acquire new functions. One aspect of gene function, that is the expression profile in gene coexpression network, has been largely unexplored for duplicate genes.     

Functional evolutionary divergence of proteolytic enzymes and their inhibitors

The amino acid residues of protease inhibitors that interact directly with proteases are hypervariable, suggesting recent positive evolutionary selection for functional differences between them. A similar hypervariability is apparent in some proteases, which may be co-evolving with the inhibitors.     

The role of cis-regulatory motifs and genetical control of expression in the divergence of yeast duplicate genes

Expression divergence of duplicate genes is widely believedto be important for their retention and evolution of new function,although the mechanism that determines their expression divergenceremains unclear. We use a genetical genomics approach to exploredivergence in genetical control of yeast duplicate genes createdby a whole-genome duplication that occurred about 100 MYA andthose with a younger duplication age. The analysis reveals thatduplicate genes have a significantly higher probability of sharingcommon genetic control than pairs of singleton genes. The expressionquantitative trait loci (eQTLs) have diverged completely fora high proportion of duplicate pairs, whereas a substantiallylarger proportion of duplicates share common regulatory motifsafter 100 Myr of divergent evolution. The similarity in bothgenetical control and cis motif structure for a duplicate pairis a reflection of its evolutionary age. This study revealsthat up to 20% of variation in expression between ancient duplicategene pairs in the yeast genome can be explained by both cismotif divergence (8%) and by trans eQTL divergence (10%). Initially,divergence in all 3 aspects of cis motif structure, trans-geneticalcontrol, and expression evolves coordinately with the codingsequence divergence of both young and old duplicate pairs. Thesefindings highlight the importance of divergence in both cismotif structure and trans-genetical control in the diverse setof mechanisms underlying the expression divergence of yeastduplicate genes.     

An improved likelihood ratio test for detecting site-specific functional divergence among clades of protein-coding genes

Maximum likelihood codon substitution models have proven useful for studying when and how protein function evolves, but they have recently been criticized on a number of fronts. The strengths and weaknesses of such methods must therefore be identified and improved upon. Here, using simulations, we show that the Clade model C versus M1a test for functional divergence among clades is prone to false positives under simple evolutionary conditions. We then propose a new null model (M2a_rel) that better accounts for among-site variation in selective constraint. We show that the revised test has an improved false-positive rate and good power. Applying this test to previously analyzed data sets of primate ribonucleases and mammalian rhodopsins reveals that some conclusions may have been misled by the original method. The improved test should prove useful for identifying patterns of divergence in selective constraint among paralogous gene families and among orthologs from ecologically divergent species.     

Duplications and functional divergence of ADP-glucose pyrophosphorylase genes in plants

Background  

ADP-glucose pyrophosphorylase (AGPase), which catalyses a rate limiting step in starch synthesis, is a heterotetramer comprised of two identical large and two identical small subunits in plants. Although the large and small subunits are equally sensitive to activity-altering amino acid changes when expressed in a bacterial system, the overall rate of non-synonymous evolution is ~2.7-fold greater for the large subunit than for the small subunit. Herein, we examine the basis for their different rates of evolution, the number of duplications in both large and small subunit genes and document changes in the patterns of AGPase evolution over time.