Asymmetric Functional Divergence of Duplicate Genes in Yeast

Functional similarities and differences among subunits of the nascent polypeptide-associated complex (NAC) of Saccharomyces cerevisiae

Protein factors bind ribosomes near the tunnel exit, facilitating protein trafficking and folding. In eukaryotes, the heterodimeric nascent polypeptide-associated complex (NAC) is the most abundant - equimolar to ribosomes. Saccharomyces cerevisiae has a minor β-type subunit (Nacβ2) in addition to abundant Nacβ1, and therefore two NAC heterodimers, α/β1 and α/β12. The additional beta NAC gene arose at the time of the whole genome duplication that occurred in the S. cerevisiae lineage. Nacβ2 has been implicated in regulating the fate of mRNA encoding ribosomal protein Rpl4 during translation via its interaction with the Caf130 subunit of the regulatory CCR4-Not complex. We found that Nacβ2 residues just C-terminal to the globular domain are required for its interaction with Caf130 and its negative effect on growth of cells lacking Acl4, the specialized chaperone for Rpl4. Substitution of these Nacβ2 residues at homologous positions in Nacβ1 results in a chimeric protein that interacts with Caf130 and slows the growth of ∆acl4 cells lacking Nacβ2. Furthermore, alteration of residues in the N-terminus of Nacβ2 or chimeric Nacβ1 previously shown to affect ribosome binding overcomes the growth defect of ∆acl4. Our results are consistent with a model in which Nacβ2's ribosome association per se, or its precise positioning, is necessary for productive recruitment of CCR4-Not via its interaction with the Caf130 subunit to drive Rpl4 mRNA degradation.

Models for the retention of duplicate genes and their biological underpinnings

Gene content in genomes changes through several different processes, with gene duplication being an important contributor to such changes. Gene duplication occurs over a range of scales from individual genes to whole genomes, and the dynamics of this process can be context dependent. Still, there are rules by which genes are retained or lost from genomes after duplication, and probabilistic modeling has enabled characterization of these rules, including their context-dependence. Here, we describe the biology and corresponding mathematical models that are used to understand duplicate gene retention and its contribution to the set of biochemical functions encoded in a genome.

How gene duplication diversifies the landscape of protein oligomeric state and function

Oligomeric proteins are central to cellular life and the duplication and divergence of their genes is a key driver of evolutionary innovations. The duplication of a gene coding for an oligomeric protein has numerous possible outcomes, which motivates questions on the relationship between structural and functional divergence. How do protein oligomeric states diversify after gene duplication? In the simple case of duplication of a homo-oligomeric protein gene, what properties can influence the fate of descendant paralogs toward forming independent homomers or maintaining their interaction as a complex? Furthermore, how are functional innovations associated with the diversification of oligomeric states? Here, we review recent literature and present specific examples in an attempt to illustrate and answer these questions.

Evidence from Drosophila Supports Higher Duplicability of Faster Evolving Genes

The faster rate of evolution of duplicated genes relative to singletons has been well documented in multiple lineages. This observation has generally been attributed to a presumed release from constraint following creation of a redundant, duplicate copy. However, it is not obvious that the relationship operates in this direction. An alternative possibility—that the faster rate of evolution predates the duplication event and the observed differences result from a higher propensity to duplicate in fast-evolving genes—has been tested in primates and in insects. However, these studies arrived at different conclusions and clarity is needed on whether these contrasting results relate to differences in methodology or legitimate biological differences between the lineages selected. Here, we test whether duplicable genes are faster evolving independent of duplication in the Drosophila lineage and find that our results support the conclusion that faster evolving genes are more likely to duplicate, in agreement with previous work in primates. Our findings indicate that this characteristic of gene duplication is not restricted to a single lineage and has broad implications for the interpretation of the impact of gene duplication. We identify a subset of “singletons” which defy the general trends and appear to be faster evolving. Further investigation implicates homology detection failure and suggests that these may be duplicable genes with unidentifiable paralogs.

Novel Antarctic yeast adapts to cold by switching energy metabolism and increasing small RNA synthesis

The novel extremophilic yeast Rhodotorula frigidialcoholis, formerly R. JG1b, was isolated from ice-cemented permafrost in University Valley (Antarctic), one of coldest and driest environments on Earth. Phenotypic and phylogenetic analyses classified R. frigidialcoholis as a novel species. To characterize its cold-adaptive strategies, we performed mRNA and sRNA transcriptomic analyses, phenotypic profiling, and assessed ethanol production at 0 and 23 °C. Downregulation of the ETC and citrate cycle genes, overexpression of fermentation and pentose phosphate pathways genes, growth without reduction of tetrazolium dye, and our discovery of ethanol production at 0 °C indicate that R. frigidialcoholis induces a metabolic switch from respiration to ethanol fermentation as adaptation in Antarctic permafrost. This is the first report of microbial ethanol fermentation utilized as the major energy pathway in response to cold and the coldest temperature reported for natural ethanol production. R. frigidialcoholis increased its diversity and abundance of sRNAs when grown at 0 versus 23 °C. This was consistent with increase in transcription of Dicer, a key protein for sRNA processing. Our results strongly imply that post-transcriptional regulation of gene expression and mRNA silencing may be a novel evolutionary fungal adaptation in the cryosphere.

Retention of duplicated genes in evolution

Gene duplication is a prevalent phenomenon across the tree of life. The processes that lead to the retention of duplicated genes are not well understood. Functional genomics approaches in model organisms, such as yeast, provide useful tools to test the mechanisms underlying retention with functional redundancy and divergence of duplicated genes, including fates associated with neofunctionalization, subfunctionalization, back-up compensation, and dosage amplification. Duplicated genes may also be retained as a consequence of structural and functional entanglement. Advances in human gene editing have enabled the interrogation of duplicated genes in the human genome, providing new tools to evaluate the relative contributions of each of these factors to duplicate gene retention and the evolution of genome structure.

Horizontally acquired cysteine synthase genes undergo functional divergence in lepidopteran herbivores

Horizontal gene transfer (HGT) plays an important role in evolutionary processes as organisms adapt to their environments, and now cases of gene duplication after HGT in eukaryotes are emerging at an increasing rate. However, the fate and roles of the duplicated genes over time in eukaryotes remain unclear. Here we conducted a comprehensive analysis of the evolution of cysteine synthase (CYS) in lepidopteran insects. Our results indicate that HGT-derived CYS genes are widespread and have undergone duplication following horizontal transfer in many lepidopteran insects. Moreover, lepidopteran CYS proteins not only have β-cyanoalanine synthase activity but also possess cysteine synthase activity that is involved in sulfur amino acid biosynthesis. Duplicated CYS genes show marked divergence in gene expression patterns and enzymatic properties, suggesting that they probably have undergone subfunctionalization and/or neofunctionalization in Lepidoptera. The gene transfer of CYS genes and subsequent duplication appears to have facilitated the adaptation of lepidopteran insects to different diets and promoted their ecological diversification. Overall, this study provides valuable insights into the ecological and evolutionary contributions of CYS in lepidopteran insects.

The divergence of alternative splicing between ohnologs in teleost fishes

BACKGROUND

Gene duplication and alternative splicing (AS) are two distinct mechanisms generating new materials for genetic innovations. The evolutionary link between gene duplication and AS is still controversial, due to utilizing duplicates from inconsistent ages of duplication events in earlier studies. With the aid of RNA-seq data, we explored evolutionary scenario of AS divergence between duplicates with ohnologs that resulted from the teleost genome duplication event in zebrafish, medaka, and stickleback.

RESULTS

Ohnologs in zebrafish have fewer AS forms compared to their singleton orthologs, supporting the function-sharing model of AS divergence between duplicates. Ohnologs in stickleback have more AS forms compared to their singleton orthologs, which supports the accelerated model of AS divergence between duplicates. The evolution of AS in ohnologs in medaka supports a combined scenario of the function-sharing and the accelerated model of AS divergence between duplicates. We also found a small number of ohnolog pairs in each of the three teleosts showed significantly asymmetric AS divergence. For example, the well-known ovary-factor gene cyp19a1a has no AS form but its ohnolog cyp19a1b has multiple AS forms in medaka, suggesting that functional divergence between duplicates might have result from AS divergence.

CONCLUSIONS

We found that a combined scenario of function-sharing and accelerated models for AS evolution in ohnologs in teleosts and rule out the independent model that assumes a lack of correlation between gene duplication and AS. Our study thus provided insights into the link between gene duplication and AS in general and ohnolog divergence in teleosts from AS perspective in particular.

Out of the testis, into the ovary: biased outcomes of gene duplication and deletion in Drosophila

Gene turnover is a key source of adaptive variation. Yet most evolutionary studies have focused on gene duplication, dismissing gene deletion as a mechanism that simply eradicates redundancy. Here, I use genome-scale sequence and multi-tissue expression data from Drosophila melanogaster and Drosophila pseudoobscura to simultaneously assess the evolutionary outcomes of gene duplication and deletion in Drosophila. I find that gene duplication is more frequent than gene deletion in both species, indicating that it may play a more important role in Drosophila evolution. However, examination of several genic properties reveals that genes likely possess distinct functions after duplication that diverge further before deletion, suggesting that loss of redundancy cannot explain a majority of gene deletion events in Drosophila. Moreover, in addition to providing support for the well-known "out of the testis" origin of young duplicate genes, analyses of gene expression profiles uncover a preferential bias against deletion of old ovary-expressed genes. Therefore, I propose a novel "into the ovary" hypothesis for gene deletion in Drosophila, in which gene deletion may promote adaptation by salvaging genes that contribute to the evolution of female reproductive phenotypes. Under this combined "out of the testis, into the ovary" evolutionary model, gene duplication and deletion work in concert to generate and maintain a balanced repertoire of genes that promote sex-specific adaptation in Drosophila.

Duplication and diversification of trehalase confers evolutionary advantages on lepidopteran insects

Gene duplication provides a major source of new genes for evolutionary novelty and ecological adaptation. However, the maintenance of duplicated genes and their relevance to adaptive evolution has long been debated. Insect trehalase (Treh) plays key roles in energy metabolism, growth, and stress recovery. Here, we show that the duplication of Treh in Lepidoptera (butterflies and moths) is linked with their adaptation to various environmental stresses. Generally, two Treh genes are present in insects: Treh1 and Treh2. We report three distinct forms of Treh in lepidopteran insects, where Treh1 was duplicated into two gene clusters (Treh1a and Treh1b). These gene clusters differ in gene expression patterns, enzymatic properties, and subcellular localizations, suggesting that the enzymes probably underwent sub- and/or neofunctionalization in the lepidopteran insects. Interestingly, selective pressure analysis provided significant evidence of positive selection on duplicate Treh1b gene in lepidopteran insect lineages. Most positively selected sites were located in the alpha-helical region, and several sites were close to the trehalose binding and catalytic sites. Subcellular adaptation of duplicate Treh1b driven by positive selection appears to have occurred as a result of selected changes in specific sequences, allowing for rapid reprogramming of duplicated Treh during evolution. Our results suggest that gene duplication of Treh and subsequent functional diversification could increase the survival rate of lepidopteran insects through various regulations of intracellular trehalose levels, facilitating their adaptation to diverse habitats. This study provides evidence regarding the mechanism by which gene family expansion can contribute to species adaptation through gene duplication and subsequent functional diversification.

Genome-wide identification and characterization of TALE superfamily genes in cotton reveals their functions in regulating secondary cell wall biosynthesis

BACKGROUND

Cotton fiber length and strength are both key traits of fiber quality, and fiber strength (FS) is tightly correlated with secondary cell wall (SCW) biosynthesis. The three-amino-acid-loop-extension (TALE) superclass homeoproteins are involved in regulating diverse biological processes in plants, and some TALE members has been identified to play a key role in regulating SCW formation. However, little is known about the functions of TALE members in cotton (Gossypium spp.).

RESULTS

In the present study, based on gene homology, 46, 47, 88 and 94 TALE superfamily genes were identified in G. arboreum, G. raimondii, G. barbadense and G. hirsutum, respectively. Phylogenetic and evolutionary analysis showed the evolutionary conservation of two cotton TALE families (including BEL1-like and KNOX families). Gene structure analysis also indicated the conservation of GhTALE members under selection. The analysis of promoter cis-elements and expression patterns suggested potential transcriptional regulation functions in fiber SCW biosynthesis and responses to some phytohormones for GhTALE proteins. Genome-wide analysis of colocalization of TALE transcription factors with SCW-related QTLs revealed that some BEL1-like genes and KNAT7 homologs may participate in the regulation of cotton fiber strength formation. Overexpression of GhKNAT7-A03 and GhBLH6-A13 significantly inhibited the synthesis of lignocellulose in interfascicular fibers of Arabidopsis. Yeast two-hybrid (Y2H) experiments showed extensive heteromeric interactions between GhKNAT7 homologs and some GhBEL1-like proteins. Yeast one-hybrid (Y1H) experiments identified the upstream GhMYB46 binding sites in the promoter region of GhTALE members and defined the downstream genes that can be directly bound and regulated by GhTALE heterodimers.

CONCLUSION

We comprehensively identified TALE superfamily genes in cotton. Some GhTALE members are predominantly expressed during the cotton fiber SCW thicking stage, and may genetically correlated with the formation of FS. Class II KNOX member GhKNAT7 can interact with some GhBEL1-like members to form the heterodimers to regulate the downstream targets, and this regulatory relationship is partially conserved with Arabidopsis. In summary, this study provides important clues for further elucidating the functions of TALE genes in regulating cotton growth and development, especially in the fiber SCW biosynthesis network, and it also contributes genetic resources to the improvement of cotton fiber quality.

Ribosomal flavours: an acquired taste for specific mRNAs?

The regulation of translation is critical in almost every aspect of gene expression. Nonetheless, the ribosome is historically viewed as a passive player in this process. However, evidence is accumulating to suggest that variations in the ribosome can have an important influence on which mRNAs are translated. Scope for variation is provided via multiple avenues, including heterogeneity at the level of both ribosomal proteins and ribosomal RNAs and their covalent modifications. Together, these variations provide the potential for hundreds, if not thousands, of flavours of ribosome, each of which could have idiosyncratic preferences for the translation of certain messenger RNAs. Indeed, perturbations to this heterogeneity appear to affect specific subsets of transcripts and manifest as cell-type-specific diseases. This review provides a historical perspective of the ribosomal code hypothesis, before outlining the various sources of heterogeneity, their regulation and functional consequences for the cell.

Pattern of expression and interaction specificity of multiple G-protein beta (Gβ) subunit isoforms with their potential target proteins reveal functional dominance of BjuGβ1 in the allotetraploid Brassica juncea

Heterotrimeric G-protein, consisting Gα, Gβ and Gγ subunits, interacts with various upstream and downstream effector (target) proteins to regulate a large array of conserved and species-specific biological functions. The targets of G-protein components are recently reported in model plant Arabidopsis thaliana; however limited information is available from crop species. In this study, we utilized yeast two-hybrid (Y2H) assay to screen the diversity of interacting partners of multiple Gβ subunit isoforms from allotetraploid Brassica juncea, a globally important oilseed and vegetable crop. The three BjuGβ genes (BjuGβ1-3), resulted from whole genome triplication event in Brassica lineage, showed distinct expression profile during plant developmental stages with maximal transcript abundance during reproductive stages. Protein-protein interaction of three BjuGβ proteins (bait) against the Y2H cDNA library (prey) identified a total of 14 and 1 non-redundant targets for BjuGβ1 and BjuGβ2, whereas BjuGβ3 screening surprisingly did not yield any genuine target, thereby suggesting functional dominance of BjuGβ1. The triplicated BjuGβ isoforms showed a high degree of interaction strength and specificity with the identified target proteins, which are known to be involved in diverse biological functions in plants. qRT-PCR analysis further indicated that the expression of BjuGβ-target genes was developmentally regulated under various tissue types studied and showed a high degree of co-expression pattern with the BjuGβ genes, particularly during flower and silique development in B. juncea. Taken together, our data provides novel insights on pattern of expression and interaction specificity governing functional divergence of multiple Gβ subunit proteins in polyploid B. juncea.

SCPP Genes and Their Relatives in Gar: Rapid Expansion of Mineralization Genes in Osteichthyans

Gar is an actinopterygian that has bone, dentin, enameloid, and ganoin (enamel) in teeth and/or scales. Mineralization of these tissues involves genes encoding various secretory calcium-binding phosphoproteins (SCPPs) in osteichthyans, but no SCPP genes have been identified in chondrichthyans to date. In the gar genome, we identified 38 SCPP genes, seven of which encode "acidic-residue-rich" proteins and 31 encode "Pro/Gln (P/Q) rich" proteins. These gar SCPP genes constitute the largest known repertoire, including many newly identified P/Q-rich genes expressed in teeth and/or scales. Among gar SCPP genes, six acidic and three P/Q-rich genes were identified as orthologs of sarcopterygian genes. The sarcopterygian orthologs of most of these acidic genes are involved in bone and/or dentin formation, and sarcopterygian orthologs of all three P/Q-rich genes participate in enamel formation. The finding of these genes in gar suggests that an elaborate SCPP gene-based genetic system for tissue mineralization was already present in stem osteichthyans. While SCPP genes have been thought to originate from ancient SPARCL1, SPARCL1L1 appears to be more closely related to these genes, because it established a structure similar to acidic SCPP genes probably in stem gnathostomes, perhaps at about the same time with the origin of tissue mineralization. Assuming enamel evolved in stem osteichthyans, all P/Q-rich SCPP genes likely arose within the osteichthyan lineage. Furthermore, the absence of acidic SCPP genes in chondrichthyans might be explained by the secondary loss of earliest acidic genes. It appears that many SCPP genes expanded rapidly in stem osteichthyans and in basal actinopterygians.

Binding interface change and cryptic variation in the evolution of protein-protein interactions

Background

Physical interactions between proteins are essential for almost all biological functions and systems. To understand the evolution of function it is therefore important to understand the evolution of molecular interactions. Of key importance is the evolution of binding specificity, the set of interactions made by a protein, since change in specificity can lead to “rewiring” of interaction networks. Unfortunately, the interfaces through which proteins interact are complex, typically containing many amino-acid residues that collectively must contribute to binding specificity as well as binding affinity, structural integrity of the interface and solubility in the unbound state.

Results

In order to study the relationship between interface composition and binding specificity, we make use of paralogous pairs of yeast proteins. Immediately after duplication these paralogues will have identical sequences and protein products that make an identical set of interactions. As the sequences diverge, we can correlate amino-acid change in the interface with any change in the specificity of binding. We show that change in interface regions correlates only weakly with change in specificity, and many variants in interfaces are functionally equivalent. We show that many of the residue replacements within interfaces are silent with respect to their contribution to binding specificity.

Conclusions

We conclude that such functionally-equivalent change has the potential to contribute to evolutionary plasticity in interfaces by creating cryptic variation, which in turn may provide the raw material for functional innovation and coevolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0608-1) contains supplementary material, which is available to authorized users.

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.

Mean field theory for biology inspired duplication-divergence network model

The duplication-divergence network model is generally thought to incorporate key ingredients underlying the growth and evolution of protein-protein interaction networks. Properties of the model have been elucidated through numerous simulation studies. However, a comprehensive theoretical study of the model is lacking. Here, we derived analytic expressions for quantities describing key characteristics of the network-the average degree, the degree distribution, the clustering coefficient, and the neighbor connectivity-in the mean-field, large-N limit of an extended version of the model, duplication-divergence complemented with heterodimerization and addition. We carried out extensive simulations and verified excellent agreement between simulation and theory except for one partial case. All four quantities obeyed power-laws even at moderate network size ( N∼10(4)), except the degree distribution, which had an additional exponential factor observed to obey power-law. It is shown that our network model can lead to the emergence of scale-free property and hierarchical modularity simultaneously, reproducing the important topological properties of real protein-protein interaction networks.

Rapid divergence and diversification of mammalian duplicate gene functions

BACKGROUND

Gene duplication provides raw material for the evolution of functional innovation. We recently developed a phylogenetic method that classifies evolutionary processes driving the retention of duplicate genes by quantifying divergence between their spatial gene expression profiles and that of their single-copy orthologous gene in a closely related sister species.

RESULTS

Here, we apply our classification method to pairs of duplicate genes in eight mammalian genomes, using data from 11 tissues to construct spatial gene expression profiles. We find that young mammalian duplicates are often functionally conserved, and that expression divergence rapidly increases over evolutionary time. Moreover, expression divergence results in increased tissue specificity, with an overrepresentation of expression in male kidney, underrepresentation of expression in female liver, and strong underrepresentation of expression in testis. Thus, duplicate genes acquire a diversity of new tissue-specific functions outside of the testis, possibly contributing to the origin of a multitude of complex phenotypes during mammalian evolution.

CONCLUSIONS

Our findings reveal that mammalian duplicate genes are initially functionally conserved, and then undergo rapid functional divergence over evolutionary time, acquiring diverse tissue-specific biological roles. These observations are in stark contrast to the much faster expression divergence and acquisition of broad housekeeping roles we previously observed in Drosophila duplicate genes. Due to the smaller effective population sizes of mammals relative to Drosophila, these analyses implicate natural selection in the functional evolution of duplicate genes.

Multispecies Analysis of Expression Pattern Diversification in the Recently Expanded Insect Ly6 Gene Family

Gene families often consist of members with diverse expression domains reflecting their functions in a wide variety of tissues. However, how the expression of individual members, and thus their tissue-specific functions, diversified during the course of gene family expansion is not well understood. In this study, we approached this question through the analysis of the duplication history and transcriptional evolution of a rapidly expanding subfamily of insect Ly6 genes. We analyzed different insect genomes and identified seven Ly6 genes that have originated from a single ancestor through sequential duplication within the higher Diptera. We then determined how the original embryonic expression pattern of the founding gene diversified by characterizing its tissue-specific expression in the beetle Tribolium castaneum, the butterfly Bicyclus anynana, and the mosquito Anopheles stephensi and those of its duplicates in three higher dipteran species, representing various stages of the duplication history (Megaselia abdita, Ceratitis capitata, and Drosophila melanogaster). Our results revealed that frequent neofunctionalization episodes contributed to the increased expression breadth of this subfamily and that these events occurred after duplication and speciation events at comparable frequencies. In addition, at each duplication node, we consistently found asymmetric expression divergence. One paralog inherited most of the tissue-specificities of the founder gene, whereas the other paralog evolved drastically reduced expression domains. Our approach attests to the power of combining a well-established duplication history with a comprehensive coverage of representative species in acquiring unequivocal information about the dynamics of gene expression evolution in gene families.

Systematic lipidomic analysis of yeast protein kinase and phosphatase mutants reveals novel insights into regulation of lipid homeostasis

The regulatory pathways required to maintain eukaryotic lipid homeostasis are largely unknown. We developed a systematic approach to uncover new players in the regulation of lipid homeostasis. Through an unbiased mass spectrometry-based lipidomic screening, we quantified hundreds of lipid species, including glycerophospholipids, sphingolipids, and sterols, from a collection of 129 mutants in protein kinase and phosphatase genes of Saccharomyces cerevisiae. Our approach successfully identified known kinases involved in lipid homeostasis and uncovered new ones. By clustering analysis, we found connections between nutrient-sensing pathways and regulation of glycerophospholipids. Deletion of members of glucose- and nitrogen-sensing pathways showed reciprocal changes in glycerophospholipid acyl chain lengths. We also found several new candidates for the regulation of sphingolipid homeostasis, including a connection between inositol pyrophosphate metabolism and complex sphingolipid homeostasis through transcriptional regulation of AUR1 and SUR1. This robust, systematic lipidomic approach constitutes a rich, new source of biological information and can be used to identify novel gene associations and function.

Relationship between gene duplicability and diversifiability in the topology of biochemical networks

Background

Selective gene duplicability, the extensive expansion of a small number of gene families, is universal. Quantitatively, the number of genes (P_(K)) with K duplicates in a genome decreases precipitously as K increases, and often follows a power law (P_(k)∝k^-α). Functional diversification, either neo- or sub-functionalization, is a major evolution route for duplicate genes.

Results

Using three lines of genomic datasets, we studied the relationship between gene duplicability and diversifiability in the topology of biochemical networks. First, we explored scenario where two pathways in the biochemical networks antagonize each other. Synthetic knockout of respective genes for the two pathways rescues the phenotypic defects of each individual knockout. We identified duplicate gene pairs with sufficient divergences that represent this antagonism relationship in the yeast S. cerevisiae. Such pairs overwhelmingly belong to large gene families, thus tend to have high duplicability. Second, we used distances between proteins of duplicate genes in the protein interaction network as a metric of their diversification. The higher a gene’s duplicate count, the further the proteins of this gene and its duplicates drift away from one another in the networks, which is especially true for genetically antagonizing duplicate genes. Third, we computed a sequence-homology-based clustering coefficient to quantify sequence diversifiability among duplicate genes – the lower the coefficient, the more the sequences have diverged. Duplicate count (K) of a gene is negatively correlated to the clustering coefficient of its duplicates, suggesting that gene duplicability is related to the extent of sequence divergence within the duplicate gene family.

Conclusion

Thus, a positive correlation exists between gene diversifiability and duplicability in the context of biochemical networks – an improvement of our understanding of gene duplicability.

Stochastic galactokinase expression underlies GAL gene induction in a GAL3 mutant of Saccharomyces cerevisiae

GAL1 and GAL3 are paralogous signal transducers that functionally inactivate Gal80p to activate the Gal4p-dependent transcriptional activation of GAL genes in Saccharomyces cerevisiae in response to galactose. Unlike a wild-type strain, the gal3∆ strain shows delayed growth kinetics as a result of the signaling function of GAL1. The mechanism ensuring that GAL1 is eventually expressed to turn on the GAL switch in the gal3∆ strain remains a paradox. Using galactose and histidine growth complementation assays, we demonstrate that 0.3% of the gal3∆ cell population responds to galactose. This is corroborated by flow cytometry and microscopic analysis. The galactose responders and nonresponders isolated from the galactose-adapted population attain the original bimodal state and this phenotype is found to be as hard wired as a genetic trait. Computational analysis suggests that the log-normal distribution in GAL4 synthesis can lead to bimodal expression of GAL80, resulting in the bimodal expression of GAL genes. Heterozygosity at the GAL80 but not at the GAL1, GAL2 or GAL4 locus alters the extent of bimodality of the gal3∆ cell population. We suggest that the asymmetric expression pattern between GAL1 and GAL3 results in the ability of S. cerevisiae to activate the GAL pathway by conferring nongenetic heterogeneity.

The evolution and functional divergence of the beta-carotene oxygenase gene family in teleost fish--exemplified by Atlantic salmon

In mammals, two carotenoid cleaving oxygenases are known; beta-carotene 15,15'-monooxygenase (BCMO1) and beta-carotene 9',10'-oxygenase (BCO2). BCMO1 is a key enzyme in vitamin A synthesis by symmetrically cleaving beta-carotene into 2 molecules of all-trans-retinal, while BCO2 is responsible for asymmetric cleavage of a broader range of carotenoids. Here, we show that the Atlantic salmon beta-carotene oxygenase (bco) gene family contains 5 members, three bco2 and two bcmo1 paralogs. Using public sequence databases, multiple bco genes were also found in several additional teleost species. Phylogenetic analysis indicates that bco2a and bco2b originate from the teleost fish specific genome duplication (FSGD or 3R), while the third and more distant paralog, bco2 like, might stem from a prior duplication event in the teleost lineage. The two bcmo1 paralogs (bcmo1 and bcmo1 like) appear to be the result of an ancient duplication event that took place before the divergence of ray-finned (Actinopterygii) and lobe-finned fish (Sarcopterygii), with subsequent nonfunctionalization and loss of one Sarcopterygii paralog. Gene expression analysis of the bcmo1 and bco2 paralogs in Atlantic salmon reveals regulatory divergence with tissue specific expression profiles, suggesting that the beta-carotene oxygenase subtypes have evolved functional divergences. We suggest that teleost fish have evolved and maintained an extended repertoire of beta-carotene oxygenases compared to the investigated Sarcopterygii species, and hypothesize that the main driver behind this functional divergence is the exposure to a diverse set of carotenoids in the aquatic environment.

Genetic robustness and functional evolution of gene duplicates

Gene duplications are a major source of evolutionary innovations. Understanding the functional divergence of duplicates and their role in genetic robustness is an important challenge in biology. Previously, analyses of genetic robustness were primarily focused on duplicates essentiality and epistasis in several laboratory conditions. In this study, we use several quantitative data sets to understand compensatory interactions between Saccharomyces cerevisiae duplicates that are likely to be relevant in natural biological populations. We find that, owing to their high functional load, close duplicates are unlikely to provide substantial backup in the context of large natural populations. Interestingly, as duplicates diverge from each other, their overall functional load is reduced. At intermediate divergence distances the quantitative decrease in fitness due to removal of one duplicate becomes smaller. At these distances, yeast duplicates display more balanced functional loads and their transcriptional control becomes significantly more complex. As yeast duplicates diverge beyond 70% sequence identity, their ability to compensate for each other becomes similar to that of random pairs of singletons.

Neofunctionalization of young duplicate genes in Drosophila

Gene duplication is a key source of genetic innovation that plays a role in the evolution of phenotypic complexity. Although several evolutionary processes can result in the long-term retention of duplicate genes, their relative contributions in nature are unknown. Here we develop a phylogenetic approach for comparing genome-wide expression profiles of closely related species to quantify the roles of conservation, neofunctionalization, subfunctionalization, and specialization in the preservation of duplicate genes. Application of our method to pairs of young duplicates in Drosophila shows that neofunctionalization, the gain of a novel function in one copy, accounts for the retention of almost two-thirds of duplicate genes. Surprisingly, novel functions nearly always originate in younger (child) copies, whereas older (parent) copies possess functions similar to those of ancestral genes. Further examination of such pairs reveals a strong bias toward RNA-mediated duplication events, implicating asymmetric duplication and positive selection in the evolution of new functions. Moreover, we show that young duplicate genes are expressed primarily in testes and that their expression breadth increases over evolutionary time. This finding supports the "out-of-testes" hypothesis, which posits that testes are a catalyst for the emergence of new genes that ultimately evolve functions in other tissues. Thus, our study highlights the importance of neofunctionalization and positive selection in the retention of young duplicates in Drosophila and illustrates how duplicates become incorporated into novel functional networks over evolutionary time.

Tissue-specific profiling reveals transcriptome alterations in Arabidopsis mutants lacking morphological phenotypes

Traditional genetic analysis relies on mutants with observable phenotypes. Mutants lacking visible abnormalities may nevertheless exhibit molecular differences useful for defining gene function. To examine this, we analyzed tissue-specific transcript profiles from Arabidopsis thaliana transcription factor gene mutants with known roles in root epidermis development, but lacking a single-gene mutant phenotype due to genetic redundancy. We discovered substantial transcriptional changes in each mutant, preferentially affecting root epidermal genes in a manner consistent with the known double mutant effects. Furthermore, comparing transcript profiles of single and double mutants, we observed remarkable variation in the sensitivity of target genes to the loss of one or both paralogous genes, including preferential effects on specific branches of the epidermal gene network, likely reflecting the pathways of paralog subfunctionalization during evolution. In addition, we analyzed the root epidermal transcriptome of the transparent testa glabra2 mutant to clarify its role in the network. These findings provide insight into the molecular basis of genetic redundancy and duplicate gene diversification at the level of a specific gene regulatory network, and they demonstrate the usefulness of tissue-specific transcript profiling to define gene function in mutants lacking informative visible changes in phenotype.

Extensive rewiring and complex evolutionary dynamics in a C. elegans multiparameter transcription factor network

Gene duplication results in two identical paralogs that diverge through mutation, leading to loss or gain of interactions with other biomolecules. Here, we comprehensively characterize such network rewiring for C. elegans transcription factors (TFs) within and across four newly delineated molecular networks. Remarkably, we find that even highly similar TFs often have different interaction degrees and partners. In addition, we find that most TF families have a member that is highly connected in multiple networks. Further, different TF families have opposing correlations between network connectivity and phylogenetic age, suggesting that they are subject to different evolutionary pressures. Finally, TFs that have similar partners in one network generally do not in another, indicating a lack of pressure to retain cross-network similarity. Our multiparameter analyses provide unique insights into the evolutionary dynamics that shaped TF networks.

Function relaxation followed by diversifying selection after whole-genome duplication in flowering plants

Episodes of whole-genome duplication (WGD) followed by gene loss dominate the evolutionary history of flowering plants. Despite the importance of understanding gene evolution following WGD, little is known about the evolutionary dynamics of this process. In this study, we analyzed duplicated genes from three WGD events in the Arabidopsis (Arabidopsis thaliana) lineage using multiple data types. Most duplicated genes that have survived from the most recent WGD (α) are under purifying selection in modern Arabidopsis populations. Using the number of identified protein-protein interactions as a proxy for functional divergence, approximately 92.7% of α-duplicated genes were diverged in function from one another in modern Arabidopsis populations, indicating that their preservation is no longer explicable by dosage balance. Dosage-balanced retention declines with antiquity of duplication: 24.1% of α-duplicated gene pairs in Arabidopsis remain in dosage balance with interacting partners, versus 12.9% and 9.4% for the earlier β-duplication and γ-triplication. GO-slim (a cut-down version of gene ontologies) terms reinforce evidence from protein-protein interactions, showing that the putatively diverged gene pairs are adapted to different cellular components. We identified a group of α-duplicated genes that show higher than average single-nucleotide polymorphism density, indicating that a period of positive selection, potentially driving functional divergence, may have preceded the current phase of purifying selection. We propose three possible paths for the evolution of duplicated genes following WGD.

To the beat of a different drum: determinants implicated in the asymmetric sequence divergence of Caenorhabditis elegans paralogs

Background

Gene duplicates often exhibit asymmetric rates of molecular evolution in their early evolutionary existence. This asymmetry in rates is thought to signify the maintenance of the ancestral function by one copy and the removal of functional constraint on the other copy, enabling it to embark on a novel evolutionary trajectory. Here I focused on a large population of evolutionarily young gene duplicates (K_S ≤ 0.14) in the Caenorhabditis elegans genome in order to conduct the first combined analysis of four predictors (evolutionary age, chromosomal location, structural resemblance between duplicates, and duplication span) which may be implicated in the asymmetric sequence divergence of paralogs at the nucleotide and amino acid level. In addition, I investigate if either paralog is equally likely to embark on a trajectory of accelerated sequence evolution or whether the derived paralog is more likely to exhibit faster sequence evolution.

Results

Three predictors (evolutionary age of duplicates, chromosomal location and duplication span) serve as major determinants of sequence asymmetry between C. elegans paralogs. Paralogs diverge asymmetrically in sequence with increasing evolutionary age, the relocation of one copy to a different chromosome and attenuated duplication spans that likely fail to capture the entire ancestral repertoire of coding sequence and regulatory elements. Furthermore, for paralogs residing on the same chromosome, opposite transcriptional orientation and increased genomic distance do not increase sequence asymmetry between paralogs. For a subset of duplicate pairs wherein the ancestral versus derived paralog could be distinguished, the derived paralogs are more likely to evolve at accelerated rates.

Conclusions

This genome-wide study of evolutionarily young duplicates stemming primarily from DNA-mediated small-scale duplication events demonstrates that genomic relocation to a new chromosome has important consequences for asymmetric divergence of paralogs, akin to paralogs arising from RNA-mediated duplication events. Additionally, the duplication span is negatively correlated with sequence rate asymmetry among paralogs, suggesting that attenuated duplication spans stemming from incomplete duplication of the ORF and/or ancestral regulatory elements further accelerate sequence divergence between paralogs. Cumulatively, derived copies exhibit accelerated rates of sequence evolution suggesting that they are primed for a divergent evolutionary trajectory by changes in structure and genomic context at inception.

A kinetic model of the evolution of a protein interaction network

Background

Known protein interaction networks have very particular properties. Old proteins tend to have more interactions than new ones. One of the best statistical representatives of this property is the node degree distribution (distribution of proteins having a given number of interactions). It has previously been shown that this distribution is very close to the sum of two distinct exponential components. In this paper, we asked: What are the possible mechanisms of evolution for such types of networks? To answer this question, we tested a kinetic model for simplified evolution of a protein interactome. Our proposed model considers the emergence of new genes and interactions and the loss of old ones. We assumed that there are generally two coexisting classes of proteins. Proteins constituting the first class are essential only for ecological adaptations and are easily lost when ecological conditions change. Proteins of the second class are essential for basic life processes and, hence, are always effectively protected against deletion. All proteins can transit between the above classes in both directions. We also assumed that the phenomenon of gene duplication is always related to ecological adaptation and that a new copy of a duplicated gene is not essential. According to this model, all proteins gain new interactions with a rate that preferentially increases with the number of interactions (the rich get richer). Proteins can also gain interactions because of duplication. Proteins lose their interactions both with and without the loss of partner genes.

Results

The proposed model reproduces the main properties of protein-protein interaction networks very well. The connectivity of the oldest part of the interaction network is densest, and the node degree distribution follows the sum of two shifted power-law functions, which is a theoretical generalization of the previous finding. The above distribution covers the wide range of values of node degrees very well, much better than a power law or generalized power law supplemented with an exponential cut-off. The presented model also relates the total number of interactome links to the total number of interacting proteins. The theoretical results were for the interactomes of A. thaliana, B. taurus, C. elegans, D. melanogaster, E. coli, H. pylori, H. sapiens, M. musculus, R. norvegicus and S. cerevisiae.

Conclusions

Using these approaches, the kinetic parameters could be estimated. Finally, the model revealed the evolutionary kinetics of proteome formation, the phenomenon of protein differentiation and the process of gaining new interactions.

Asymmetric functional divergence of young, dispersed gene duplicates in Arabidopsis thaliana

One prediction of the classic Ohno model of gene duplication predicts that new genes form from the asymmetric functional divergence of a newly arisen, redundant duplicate locus. In order to understand the mechanisms which give rise to functional divergence of newly formed dispersed duplicates, we assessed the expression and molecular evolutionary divergence of a suite of 19 highly similar dispersed duplicates in Arabidopsis thaliana. These duplicates have a K sil equal to or less than 5 % and are specific to the A. thaliana lineage; thus, they predictably represent some of the youngest duplicates in the A. thaliana genome. We found that the majority of young duplicate loci exhibit asymmetric expression patterns, with the daughter locus exhibiting reduced expression across all tissues analyzed relative to the progenitor locus or simply not expressed. Furthermore, daughter loci, on the whole, have significantly more nonsynonymous substitutions than the progenitor loci. We also identified four pairs of loci which exhibit significant (P < 0.05) evolutionary rate asymmetry, three of which exhibit elevated dN/dS in the duplicate copy. We suggest, based on these data, that functional diversification initially takes the form of asymmetric regulatory divergence that can be a direct consequence of the mode of duplication. The reduced and/or absence of expression in the daughter copy relaxes functional constraint on its protein coding sequence leading to the asymmetric accumulation of nonsynonymous mutations. Thus, our data both affirm Ohno's prediction while explaining the mechanism by which functional divergence initially occurs following duplication for dispersed gene 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.

Transcriptional divergence plays a role in the rewiring of protein interaction networks after gene duplication

Gene duplication plays a key role in the evolution of protein-protein interaction (PPI) networks. After a gene duplication event, paralogous proteins may diverge through the gain and loss of PPIs. This divergence can be explained by two non-exclusive mechanisms. First, mutations may accumulate in the coding sequences of the paralogs and affect their protein sequences, which can modify, for instance, their binding interfaces and thus their interaction specificity. Second, mutations may accumulate in the non-coding region of the genes and affect their regulatory sequences. The resulting changes in expression profiles can lead to paralogous proteins being differentially expressed and occurring in the cell with different sets of potential interaction partners. These changes could also alter splicing regulation and lead to the inclusion or exclusion of alternative exons. The evolutionary role of these regulatory mechanisms remains largely unexplored. We use bioinformatics analyses of existing PPI data and proteome-wide PPI screening to show that the divergence of transcriptional regulation between paralogs plays a significant role in determining their PPI specificity. Because many gene duplication events are followed by rapid changes in transcriptional regulation, our results suggest that PPI networks may be rewired by gene duplication, without the need for protein to diverge in their binding specificities. This article is part of a Special Issue entitled: From protein structures to clinical applications.

Gene duplications contribute to the overrepresentation of interactions between proteins of a similar age

BACKGROUND

The study of biological networks and how they have evolved is fundamental to our understanding of the cell. By investigating how proteins of different ages are connected in the protein interaction network, one can infer how that network has expanded in evolution, without the need for explicit reconstruction of ancestral networks. Studies that implement this approach show that proteins are often connected to proteins of a similar age, suggesting a simultaneous emergence of interacting proteins. There are several theories explaining this phenomenon, but despite the importance of gene duplication in genome evolution, none consider protein family dynamics as a contributing factor.

RESULTS

In an S. cerevisiae protein interaction network we investigate to what extent edges that arise from duplication events contribute to the observed tendency to interact with proteins of a similar age. We find that part of this tendency is explained by interactions between paralogs. Age is usually defined on the level of protein families, rather than individual proteins, hence paralogs have the same age. The major contribution however, is from interaction partners that are shared between paralogs. These interactions have most likely been conserved after a duplication event. To investigate to what extent a nearly neutral process of network growth can explain these results, we adjust a well-studied network growth model to incorporate protein families. Our model shows that the number of edges between paralogs can be amplified by subsequent duplication events, thus explaining the overrepresentation of interparalog edges in the data. The fact that interaction partners shared by paralogs are often of the same age as the paralogs does not arise naturally from our model and needs further investigation.

CONCLUSION

We amend previous theories that explain why proteins of a similar age prefer to interact by demonstrating that this observation can be partially explained by gene duplication events. There is an ongoing debate on whether the protein interaction network is predominantly shaped by duplication and subfunctionalization or whether network rewiring is most important. Our analyses of S. cerevisiae protein interaction networks demonstrate that duplications have influenced at least one property of the protein interaction network: how proteins of different ages are connected.

Pervasive indels and their evolutionary dynamics after the fish-specific genome duplication

Insertions and deletions (indels) in protein-coding genes are important sources of genetic variation. Their role in creating new proteins may be especially important after gene duplication. However, little is known about how indels affect the divergence of duplicate genes. We here study thousands of duplicate genes in five fish (teleost) species with completely sequenced genomes. The ancestor of these species has been subject to a fish-specific genome duplication (FSGD) event that occurred approximately 350 Ma. We find that duplicate genes contain at least 25% more indels than single-copy genes. These indels accumulated preferentially in the first 40 my after the FSGD. A lack of widespread asymmetric indel accumulation indicates that both members of a duplicate gene pair typically experience relaxed selection. Strikingly, we observe a 30-80% excess of deletions over insertions that is consistent for indels of various lengths and across the five genomes. We also find that indels preferentially accumulate inside loop regions of protein secondary structure and in regions where amino acids are exposed to solvent. We show that duplicate genes with high indel density also show high DNA sequence divergence. Indel density, but not amino acid divergence, can explain a large proportion of the tertiary structure divergence between proteins encoded by duplicate genes. Our observations are consistent across all five fish species. Taken together, they suggest a general pattern of duplicate gene evolution in which indels are important driving forces of evolutionary change.

Redundancy and rewiring of genetic networks following genome-wide duplication events

Polyploidy or whole-genome duplication is a frequent phenomenon within the plant kingdom and has been associated with the occurrence of evolutionary novelty and increase in biological complexity. Because genome-wide duplication events duplicate whole molecular networks it is of interest to investigate how these networks evolve subsequent to such events. Although genome duplications are generally followed by massive gene loss, at least part of the network is usually retained in duplicate and can rewire to execute novel functions. Alternatively, the network can remain largely redundant and as such confer robustness against mutations. The increasing availability of high-throughput data makes it possible to study evolution following whole genome duplication events at the network level. Here we discuss how the use of 'omics' data in network analysis can provide novel insights on network redundancy and rewiring and conclude with some directions for future research.

Characterization and evolutionary analysis of protein-protein interaction networks

While researchers have known the importance of the protein-protein interaction for decades, recent innovations in large-scale screening techniques have caused a shift in the paradigm of protein function analysis. Where the focus was once on the individual protein, attention is now directed to the surrounding network of protein associations. As protein interaction networks can provide useful insights into the potential function of and phenotypes associated with proteins, the increasing availability of large-scale protein interaction data suggests that molecular biologists can extract more meaningful hypotheses through examination of these large networks. Further, increasing availability of high-quality protein interaction data in multiple species has allowed interpretation of the properties of networks (i.e., the presence of hubs and modularity) from an evolutionary perspective. In this chapter, we discuss major previous findings derived from analyses of large-scale protein interaction data, focusing on approaches taken by landmark assays in evaluating the structure and evolution of these networks. We then outline basic techniques for protein interaction network analysis with the goal of pointing out the benefits and potential limitations of these approaches. As the majority of large-scale protein interaction data has been generated in budding yeast, literature described here focuses on this important model organism with references to other species included where possible.

Cloning, characterization and gene expression of a metallothionein isoform in the edible cockle Cerastoderma edule after cadmium or mercury exposure

Metallothionein (MT) genes encode crucial metal-binding proteins ubiquitously expressed in living organisms and which play important roles in homeostasis of essential metals and detoxification processes. Here, the molecular organization of the first metallothionein gene of the edible cockle Cerastoderma edule and its expression after cadmium (Cd) or mercury (Hg) exposures were determined. The resulting sequence (Cemt1) exhibits unusual features. The full length cDNA encodes a protein of 73 amino acids with nine classical Cys-X((1-3))-Cys motifs, but also one Cys-Cys not generally found in molluscan MT. Moreover, characterization of the molecular organization of the Cemt1 gene revealed two different alleles (A1 and A2) with length differences due to large deletion events in their intronic sequences involving direct Short Interspersed repeated Elements (SINE), while their exonic sequences were identical. To our knowledge, such large excision mechanisms have never before been reported in a bivalve gene sequence. After 10 days of Cd exposure at environmentally relevant doses, quantitative real-time PCR revealed a strong induction of Cemt1 in gills of C. edule. Surprisingly, neither induction of the Cemt1 gene nor of MT protein was shown after Hg exposure, despite the fact that this organism is able to bioaccumulate a high amount of this trace metal which is theoretically one of the most powerful inducers of MT biosynthesis.

Selection for higher gene copy number after different types of plant gene duplications

The evolutionary origins of the multitude of duplicate genes in the plant genomes are still incompletely understood. To gain an appreciation of the potential selective forces acting on these duplicates, we phylogenetically inferred the set of metabolic gene families from 10 flowering plant (angiosperm) genomes. We then compared the metabolic fluxes for these families, predicted using the Arabidopsis thaliana and Sorghum bicolor metabolic networks, with the families' duplication propensities. For duplications produced by both small scale (small-scale duplications) and genome duplication (whole-genome duplications), there is a significant association between the flux and the tendency to duplicate. Following this global analysis, we made a more fine-scale study of the selective constraints observed on plant sodium and phosphate transporters. We find that the different duplication mechanisms give rise to differing selective constraints. However, the exact nature of this pattern varies between the gene families, and we argue that the duplication mechanism alone does not define a duplicated gene's subsequent evolutionary trajectory. Collectively, our results argue for the interplay of history, function, and selection in shaping the duplicate gene evolution in plants.

Local synteny and codon usage contribute to asymmetric sequence divergence of Saccharomyces cerevisiae gene duplicates

Background

Duplicated genes frequently experience asymmetric rates of sequence evolution. Relaxed selective constraints and positive selection have both been invoked to explain the observation that one paralog within a gene-duplicate pair exhibits an accelerated rate of sequence evolution. In the majority of studies where asymmetric divergence has been established, there is no indication as to which gene copy, ancestral or derived, is evolving more rapidly. In this study we investigated the effect of local synteny (gene-neighborhood conservation) and codon usage on the sequence evolution of gene duplicates in the S. cerevisiae genome. We further distinguish the gene duplicates into those that originated from a whole-genome duplication (WGD) event (ohnologs) versus small-scale duplications (SSD) to determine if there exist any differences in their patterns of sequence evolution.

Results

For SSD pairs, the derived copy evolves faster than the ancestral copy. However, there is no relationship between rate asymmetry and synteny conservation (ancestral-like versus derived-like) in ohnologs. mRNA abundance and optimal codon usage as measured by the CAI is lower in the derived SSD copies relative to ancestral paralogs. Moreover, in the case of ohnologs, the faster-evolving copy has lower CAI and lowered expression.

Conclusions

Together, these results suggest that relaxation of selection for codon usage and gene expression contribute to rate asymmetry in the evolution of duplicated genes and that in SSD pairs, the relaxation of selection stems from the loss of ancestral regulatory information in the derived copy.

Differences in the number of intrinsically disordered regions between yeast duplicated proteins, and their relationship with functional divergence

Background

Intrinsically disordered regions are enriched in short interaction motifs that play a critical role in many protein-protein interactions. Since new short interaction motifs may easily evolve, they have the potential to rapidly change protein interactions and cellular signaling. In this work we examined the dynamics of gain and loss of intrinsically disordered regions in duplicated proteins to inspect if changes after genome duplication can create functional divergence. For this purpose we used Saccharomyces cerevisiae and the outgroup species Lachancea kluyveri.

Principal Findings

We find that genes duplicated as part of a genome duplication (ohnologs) are significantly more intrinsically disordered than singletons (p<2.2^e-16, Wilcoxon), reflecting a preference for retaining intrinsically disordered proteins in duplicate. In addition, there have been marked changes in the extent of intrinsic disorder following duplication. A large number of duplicated genes have more intrinsic disorder than their L. kluyveri ortholog (29% for duplicates versus 25% for singletons) and an even greater number have less intrinsic disorder than the L. kluyveri ortholog (37% for duplicates versus 25% for singletons). Finally, we show that the number of physical interactions is significantly greater in the more intrinsically disordered ohnolog of a pair (p = 0.003, Wilcoxon).

Conclusion

This work shows that intrinsic disorder gain and loss in a protein is a mechanism by which a genome can also diverge and innovate. The higher number of interactors for proteins that have gained intrinsic disorder compared with their duplicates may reflect the acquisition of new interaction partners or new functional roles.

Evidence for short-time divergence and long-time conservation of tissue-specific expression after gene duplication

Gene duplication is one of the main mechanisms by which genomes can acquire novel functions. It has been proposed that the retention of gene duplicates can be associated to processes of tissue expression divergence. These models predict that acquisition of divergent expression patterns should be acquired shortly after the duplication, and that larger divergence in tissue expression would be expected for paralogs, as compared to orthologs of a similar age. Many studies have shown that gene duplicates tend to have divergent expression patterns and that gene family expansions are associated with high levels of tissue specificity. However, the timeframe in which these processes occur have rarely been investigated in detail, particularly in vertebrates, and most analyses do not include direct comparisons of orthologs as a baseline for the expected levels of tissue specificity in absence of duplications. To assess the specific contribution of duplications to expression divergence, we combine here phylogenetic analyses and expression data from human and mouse. In particular, we study differences in spatial expression among human-mouse paralogs, specifically duplicated after the radiation of mammals, and compare them to pairs of orthologs in the same species. Our results show that gene duplication leads to increased levels of tissue specificity and that this tends to occur promptly after the duplication event.

Genetic interactions reveal the evolutionary trajectories of duplicate genes

Duplicate genes show significantly fewer interactions than singleton genes, and functionally similar duplicates can exhibit dissimilar profiles because common interactions are ‘hidden' due to buffering.

Genetic interaction profiles provide insights into evolutionary mechanisms of duplicate retention by distinguishing duplicates under dosage selection from those retained because of some divergence in function.

The genetic interactions of duplicate genes evolve in an extremely asymmetric way and the directionality of this asymmetry correlates well with other evolutionary properties of duplicate genes.

Genetic interaction profiles can be used to elucidate the divergent function of specific duplicate pairs.

Gene duplication and divergence serves as a primary source for new genes and new functions, and as such has broad implications on the evolutionary process. Duplicate genes within S. cerevisiae have been shown to retain a high degree of similarity with regard to many of their functional properties (Papp et al, 2004; Guan et al, 2007; Wapinski et al, 2007; Musso et al, 2008), and perturbation of duplicate genes has been shown to result in smaller fitness defects than singleton genes (Gu et al, 2003; DeLuna et al, 2008; Dean et al, 2008; Musso et al, 2008). Individual genetic interactions between pairs of genes and profiles of such interactions across the entire genome provide a new context in which to examine the properties of duplicate compensation.

In this study we use the most recent and comprehensive set of genetic interactions in yeast produced to date (Costanzo et al, 2010) to address questions of duplicate retention and redundancy. We show that the ability for duplicate genes to buffer the deletion of a partner has three main consequences. First it agrees with previous work demonstrating that a high proportion of duplicate pairs are synthetic lethal, a classic indication of the ability to buffer one another functionally (DeLuna et al, 2008; Dean et al, 2008; Musso et al, 2008). Second, it reduces the number of genetic interactions observed between duplicate genes and the rest of the genome by masking interactions relating to common function from experimental detection. Third, this buffering of common interactions serves to reduce profile similarity in spite of common function (Figure 1). The compensatory ability of functionally similar duplicates buffers genetic interactions related to their common function (reducing the number of genetic interactions overall), while allowing the measurement of interactions related to any divergent function. Thus, even functionally similar duplicates may have dissimilar genetic interaction profiles. As previously surmised (Ihmels et al, 2007), duplicate genes under selection for dosage amplification have differing profile characteristics. We show that dosage-mediated duplicates have much higher genetic interaction profile similarity than do other duplicate pairs. Furthermore, we show in a comparison with local neighbors on a protein–protein interaction (PPI) network, that although dosage-mediated duplicates more often have higher similarity to each other than they do to their neighbors, the reverse is true for duplicates in general. That is, slightly divergent duplicate genes more often exhibit a higher similarity with a common neighbor on the PPI network than they do with each other, and that observation is consistent with the idea that common interactions are buffered while interactions corresponding to divergent functions are observed.

We then asked whether duplicates' genetic interactions that are not buffered appear in a symmetric or an asymmetric fashion. Previous work has established asymmetric patterns with regard to PPI degree (Wagner, 2002; He and Zhang, 2005), sequence divergence (Conant and Wagner, 2003; Zhang et al, 2003; Kellis et al, 2004; Scannell and Wolfe, 2008) and expression patterns (Gu et al, 2002b; Tirosh and Barkai, 2007). Although genetic interactions are further removed from mechanism than protein–protein interactions, for example, they do offer a more direct measurement of functional consequence and, thus, may give a better indication of the functional differences between a duplicate pair. We found that duplicates exhibit a strikingly asymmetric pattern of genetic interactions, with the ratio of interactions between sisters commonly exceeding 7:1 (Figure 4A). The observations differ significantly from random simulations in which genetic interactions were redistributed between sisters with equal probability (Figure 4A). Moreover, the directionality of this interaction asymmetry agrees with other physiological properties of duplicate pairs. For example, the sister with more genetic interactions also tends to have more protein–protein interactions and also tends to evolve at a slower rate (Figure 4B).

Genetic interaction degree and profiles can be used to understand the functional divergence of particular duplicates pairs. As a case example, we consider the whole-genome-duplication pair CIK1–VIK1. Each of these genes encode proteins that form distinct heterodimeric complexes with the microtubule motor protein Kar3 (Manning et al, 1999). Although each of these proteins depend on a direct physical interaction with Kar3, Cik1 has a much higher profile similarity to Kar3 than does Vik1 (r=0.5 and r=0.3, respectively). Consistent with its higher similarity, Δcik1 and Δkar3 exhibit several similar phenotypes, including abnormally short spindles, chromosome loss and delayed cell cycle progression (Page et al, 1994; Manning et al, 1999). In contrast, a Δvik1 mutant strain exhibits no overt phenotype (Manning et al, 1999).

The characterization of functional redundancy and divergence between duplicate genes is an important step in understanding the evolution of genetic systems. Large-scale genetic network analysis in Saccharomyces cerevisiae provides a powerful perspective for addressing these questions through quantitative measurements of genetic interactions between pairs of duplicated genes, and more generally, through the study of genome-wide genetic interaction profiles associated with duplicated genes. We show that duplicate genes exhibit fewer genetic interactions than other genes because they tend to buffer one another functionally, whereas observed interactions are non-overlapping and reflect their divergent roles. We also show that duplicate gene pairs are highly imbalanced in their number of genetic interactions with other genes, a pattern that appears to result from asymmetric evolution, such that one duplicate evolves or degrades faster than the other and often becomes functionally or conditionally specialized. The differences in genetic interactions are predictive of differences in several other evolutionary and physiological properties of duplicate pairs.

Difference in gene duplicability may explain the difference in overall structure of protein-protein interaction networks among eukaryotes

Background

A protein-protein interaction network (PIN) was suggested to be a disassortative network, in which interactions between high- and low-degree nodes are favored while hub-hub interactions are suppressed. It was postulated that a disassortative structure minimizes unfavorable cross-talks between different hub-centric functional modules and was positively selected in evolution. However, by re-examining yeast PIN data, several researchers reported that the disassortative structure observed in a PIN might be an experimental artifact. Therefore, the existence of a disassortative structure and its possible evolutionary mechanism remains unclear.

Results

In this study, we investigated PINs from the yeast, worm, fly, human, and malaria parasite including four different yeast PIN datasets. The analyses showed that the yeast, worm, fly, and human PINs are disassortative while the malaria parasite PIN is not. By conducting simulation studies on the basis of a duplication-divergence model, we demonstrated that a preferential duplication of low- and high-degree nodes can generate disassortative and non-disassortative networks, respectively. From this observation, we hypothesized that the difference in degree dependence on gene duplications accounts for the difference in assortativity of PINs among species. Comparison of 55 proteomes in eukaryotes revealed that genes with lower degrees showed higher gene duplicabilities in the yeast, worm, and fly, while high-degree genes tend to have high duplicabilities in the malaria parasite, supporting the above hypothesis.

Conclusions

These results suggest that disassortative structures observed in PINs are merely a byproduct of preferential duplications of low-degree genes, which might be caused by an organism's living environment.

Improving evolutionary models of protein interaction networks

MOTIVATION

Theoretical models of biological networks are valuable tools in evolutionary inference. Theoretical models based on gene duplication and divergence provide biologically plausible evolutionary mechanics. Similarities found between empirical networks and their theoretically generated counterpart are considered evidence of the role modeled mechanics play in biological evolution. However, the method by which these models are parameterized can lead to questions about the validity of the inferences. Selecting parameter values in order to produce a particular topological value obfuscates the possibility that the model may produce a similar topology for a large range of parameter values. Alternately, a model may produce a large range of topologies, allowing (incorrect) parameter values to produce a valid topology from an otherwise flawed model. In order to lend biological credence to the modeled evolutionary mechanics, parameter values should be derived from the empirical data. Furthermore, recent work indicates that the timing and fate of gene duplications are critical to proper derivation of these parameters.

RESULTS

We present a methodology for deriving evolutionary rates from empirical data that is used to parameterize duplication and divergence models of protein interaction network evolution. Our method avoids shortcomings of previous methods, which failed to consider the effect of subsequent duplications. From our parameter values, we find that concurrent and existing existing duplication and divergence models are insufficient for modeling protein interaction network evolution. We introduce a model enhancement based on heritable interaction sites on the surface of a protein and find that it more closely reflects the high clustering found in the empirical network.

Targeted inactivation of Mapk4 in mice reveals specific nonredundant functions of Erk3/Erk4 subfamily mitogen-activated protein kinases

Erk4 and Erk3 are atypical members of the mitogen-activated protein (MAP) kinase family. The high sequence identity of Erk4 and Erk3 proteins and the similar organization of their genes imply that the two protein kinases are paralogs. Recently, we have shown that Erk3 function is essential for neonatal survival and critical for the establishment of fetal growth potential and pulmonary function. To investigate the specific functions of Erk4, we have generated mice with a targeted disruption of the Mapk4 gene. We show that Erk4-deficient mice are viable and fertile and exhibit no gross morphological or physiological anomalies. Loss of Erk4 is not compensated by changes in Erk3 expression or activity during embryogenesis or in adult tissues. We further demonstrate that additional loss of Erk4 does not exacerbate the fetal growth restriction and pulmonary immaturity phenotypes of Erk3(-/-) mice and does not compromise the viability of Erk3(+/-) neonates. Interestingly, behavioral phenotyping revealed that Erk4-deficient mice manifest depression-like behavior in the forced-swimming test. Our analysis indicates that the MAP kinase Erk4 is dispensable for mouse embryonic development and reveals that Erk3 and Erk4 have acquired specialized functions through evolutionary diversification.

Asymmetric and non-uniform evolution of recently duplicated human genes

BACKGROUND

Gene duplications are a source of new genes and protein functions. The innovative role of duplication events makes families of paralogous genes an interesting target for studies in evolutionary biology. Here we study global trends in the evolution of human genes that resulted from recent duplications.

RESULTS

The pressure of negative selection is weaker during a short time immediately after a duplication event. Roughly one fifth of genes in paralogous gene families are evolving asymmetrically: one of the proteins encoded by two closest paralogs accumulates amino acid substitutions significantly faster than its partner. This asymmetry cannot be explained by differences in gene expression levels. In asymmetric gene pairs the number of deleterious mutations is increased in one copy, while decreased in the other copy as compared to genes constituting non-asymmetrically evolving pairs. The asymmetry in the rate of synonymous substitutions is much weaker and not significant.

CONCLUSIONS

The increase of negative selection pressure over time after a duplication event seems to be a major trend in the evolution of human paralogous gene families. The observed asymmetry in the evolution of paralogous genes shows that in many cases one of two gene copies remains practically unchanged, while the other accumulates functional mutations. This supports the hypothesis that slowly evolving gene copies preserve their original functions, while fast evolving copies obtain new specificities or functions.