Pervasive and persistent redundancy among duplicated genes in yeast

Transcriptomic signatures of schizophrenia revealed by dopamine perturbation in an ex vivo model

The dopaminergic hypothesis of schizophrenia (SZ) postulates that dopaminergic over activity causes psychosis, a central feature of SZ, based on the observation that blocking dopamine (DA) improves psychotic symptoms. DA is known to have both receptor- and non-receptor-mediated effects, including oxidative mechanisms that lead to apoptosis. The role of DA-mediated oxidative processes in SZ has been little studied. Here, we have used a cell perturbation approach and measured transcriptomic profiles by RNAseq to study the effect of DA exposure on transcription in B-cell transformed lymphoblastoid cell lines (LCLs) from 514 SZ cases and 690 controls. We found that DA had widespread effects on both cell growth and gene expression in LCLs. Overall, 1455 genes showed statistically significant differential DA response in SZ cases and controls. This set of differentially expressed genes is enriched for brain expression and for functions related to immune processes and apoptosis, suggesting that DA may play a role in SZ pathogenesis through modulating those systems. Moreover, we observed a non-significant enrichment of genes near genome-wide significant SZ loci and with genes spanned by SZ-associated copy number variants (CNVs), which suggests convergent pathogenic mechanisms detected by both genetic association and gene expression. The study suggests a novel role of DA in the biological processes of immune and apoptosis that may be relevant to SZ pathogenesis. Furthermore, our results show the utility of pathophysiologically relevant perturbation experiments to investigate the biology of complex mental disorders.

Genetic Interaction Network as an Important Determinant of Gene Order in Genome Evolution

Although it is generally accepted that eukaryotic gene order is not random, the basic principles of gene arrangement on a chromosome remain poorly understood. Here, we extended existing population genetics theories that were based on two-locus models and proposed a hypothesis that genetic interaction networks drive the evolution of eukaryotic gene order. We predicted that genes with positive epistasis would move toward each other in evolution, during which a negative correlation between epistasis and gene distance formed. We tested and confirmed our prediction with computational simulations and empirical data analyses. Importantly, we demonstrated that gene order in the budding yeast could be successfully predicted from the genetic interaction network. Taken together, our study reveals the role of the genetic interaction network in the evolution of gene order, extends our understanding of the encoding principles in genomes, and potentially offers new strategies to improve synthetic biology.

INCURVATA11 and CUPULIFORMIS2 Are Redundant Genes That Encode Epigenetic Machinery Components in Arabidopsis

All critical developmental and physiological events in a plant's life cycle depend on the proper activation and repression of specific gene sets, and this often involves epigenetic mechanisms. Some Arabidopsis thaliana mutants with disorders of the epigenetic machinery exhibit pleiotropic defects, including incurved leaves and early flowering, due to the ectopic and heterochronic derepression of developmental regulators. Here, we studied one such mutant class, the incurvata11 (icu11) loss-of-function mutants. We have identified ICU11 as the founding member of a small gene family that we have named CUPULIFORMIS (CP). This family is part of the 2-oxoglutarate/Fe(II)-dependent dioxygenase superfamily. ICU11 and its closest paralog, CP2, have unequally redundant functions: although cp2 mutants are phenotypically wild type, icu11 cp2 double mutants skip vegetative development and flower upon germination. This phenotype is reminiscent of loss-of-function mutants of the Polycomb-group genes EMBRYONIC FLOWER1 (EMF1) and EMF2 Double mutants harboring icu11 alleles and loss-of-function alleles of genes encoding components of the epigenetic machinery exhibit synergistic, severe phenotypes, and some are similar to those of emf mutants. Hundreds of genes are misexpressed in icu11 plants, including SEPALLATA3 (SEP3), and derepression of SEP3 causes the leaf phenotype of icu11 ICU11 and CP2 are nucleoplasmic proteins that act as epigenetic repressors through an unknown mechanism involving histone modification, but not DNA methylation.

Whole-Genome Analysis of Three Yeast Strains Used for Production of Sherry-Like Wines Revealed Genetic Traits Specific to Flor Yeasts

Flor yeast strains represent a specialized group of Saccharomyces cerevisiae yeasts used for biological wine aging. We have sequenced the genomes of three flor strains originated from different geographic regions and used for production of sherry-like wines in Russia. According to the obtained phylogeny of 118 yeast strains, flor strains form very tight cluster adjacent to the main wine clade. SNP analysis versus available genomes of wine and flor strains revealed 2,270 genetic variants in 1,337 loci specific to flor strains. Gene ontology analysis in combination with gene content evaluation revealed a complex landscape of possibly adaptive genetic changes in flor yeast, related to genes associated with cell morphology, mitotic cell cycle, ion homeostasis, DNA repair, carbohydrate metabolism, lipid metabolism, and cell wall biogenesis. Pangenomic analysis discovered the presence of several well-known "non-reference" loci of potential industrial importance. Events of gene loss included deletions of asparaginase genes, maltose utilization locus, and FRE-FIT locus involved in iron transport. The latter in combination with a flor-yeast-specific mutation in the Aft1 transcription factor gene is likely to be responsible for the discovered phenotype of increased iron sensitivity and improved iron uptake of analyzed strains. Expansion of the coding region of the FLO11 flocullin gene and alteration of the balance between members of the FLO gene family are likely to positively affect the well-known propensity of flor strains for velum formation. Our study provides new insights in the nature of genetic variation in flor yeast strains and demonstrates that different adaptive properties of flor yeast strains could have evolved through different mechanisms of genetic variation.

IGF-Binding Proteins: Why Do They Exist and Why Are There So Many?

Insulin-like growth factors (IGFs) are key growth-promoting peptides that act as both endocrine hormones and autocrine/paracrine growth factors. In the bloodstream and in local tissues, most IGF molecules are bound by one of the members of the IGF-binding protein (IGFBP) family, of which six distinct types exist. These proteins bind to IGF with an equal or greater affinity than the IGF1 receptor and are thus in a key position to regulate IGF signaling globally and locally. Binding to an IGFBP increases the half-life of IGF in the circulation and blocks its potential binding to the insulin receptor. In addition to these classical roles, IGFBPs have been shown to modulate IGF signaling locally under various conditions. Although members of the IGFBP family share significant sequence homology, they each have unique structural features and play distinct roles. These IGFBP genes also have different modes of regulation and distinct expression patterns. Some IGFBPs have been found to bind to their own receptors or to translocate into the interior compartments of cells where they may execute IGF-independent actions. In spite of this functional and regulatory diversity, it has been puzzling that loss-of-function studies have yielded relatively little information about the physiological functions of IGFBPs. In this review, we suggest that evolution has tended to retain an array of IGFBPs in order to facilitate fine-tuning of IGF signaling. We explore the emerging explanation that many IGFBP functions have evolved to allow the targeted adjustment of IGF signaling under stressful or irregular conditions, which would likely not be revealed in a standard laboratory setting.

Genome-wide identification and expression profile analysis of CCH gene family in Populus

Copper plays key roles in plant physiological activities. To maintain copper cellular homeostasis, copper chaperones have important functions in binding and transporting copper to target proteins. Detailed characterization and function analysis of a copper chaperone, CCH, is presently limited to Arabidopsis. This study reports the identification of 21 genes encoding putative CCH proteins in Populus trichocarpa. Besides sharing the conserved metal-binding motif MXCXXC and forming a βαββαβ secondary structure at the N-terminal, all the PtCCHs possessed the plant-exclusive extended C-terminal. Based on their gene structure, conserved motifs, and phylogenetic analysis, the PtCCHs were divided into three subgroups. Our analysis indicated that whole-genome duplication and tandem duplication events likely contributed to expansion of the CCH gene family in Populus. Tissue-specific data from PlantGenIE revealed that PtCCH genes had broad expression patterns in different tissues. Quantitative real-time RT-PCR (qRT-PCR) analysis revealed that PnCCH genes of P. simonii × P. nigra also had different tissue-specific expression traits, as well as different inducible-expression patterns in response to copper stresses (excessive and deficiency). In summary, our study of CCH genes in the Populus genome provides a comprehensive analysis of this gene family, and lays an important foundation for further investigation of their roles in copper homeostasis of poplar.

Accessing the Phenotype Gap: Enabling Systematic Investigation of Paralog Functional Complexity with CRISPR

Single-gene knockout experiments can fail to reveal function in the context of redundancy, which is frequently observed among duplicated genes (paralogs) with overlapping functions. We discuss the complexity associated with studying paralogs and outline how recent advances in CRISPR will help address the "phenotype gap" and impact biomedical research.

Growth condition dependency is the major cause of non-responsiveness upon genetic perturbation

Investigating the role and interplay between individual proteins in biological processes is often performed by assessing the functional consequences of gene inactivation or removal. Depending on the sensitivity of the assay used for determining phenotype, between 66% (growth) and 53% (gene expression) of Saccharomyces cerevisiae gene deletion strains show no defect when analyzed under a single condition. Although it is well known that this non-responsive behavior is caused by different types of redundancy mechanisms or by growth condition/cell type dependency, it is not known what the relative contribution of these different causes is. Understanding the underlying causes of and their relative contribution to non-responsive behavior upon genetic perturbation is extremely important for designing efficient strategies aimed at elucidating gene function and unraveling complex cellular systems. Here, we provide a systematic classification of the underlying causes of and their relative contribution to non-responsive behavior upon gene deletion. The overall contribution of redundancy to non-responsive behavior is estimated at 29%, of which approximately 17% is due to homology-based redundancy and 12% is due to pathway-based redundancy. The major determinant of non-responsiveness is condition dependency (71%). For approximately 14% of protein complexes, just-in-time assembly can be put forward as a potential mechanistic explanation for how proteins can be regulated in a condition dependent manner. Taken together, the results underscore the large contribution of growth condition requirement to non-responsive behavior, which needs to be taken into account for strategies aimed at determining gene function. The classification provided here, can also be further harnessed in systematic analyses of complex cellular systems.

Paralog-Specific Functions of RPL7A and RPL7B Mediated by Ribosomal Protein or snoRNA Dosage in Saccharomyces cerevisiae

Most ribosomal proteins in Saccharomyces cerevisiae are encoded by two paralogs that additively produce the optimal protein level for cell growth. Nonetheless, deleting one paralog of most ribosomal protein gene pairs results in a variety of phenotypes not observed when the other paralog is deleted. To determine whether paralog-specific phenotypes associated with deleting RPL7A or RPL7B stem from distinct functions or different levels of the encoded isoforms, the coding region and introns of one paralog, including an intron-embedded snoRNA (small nucleolar RNA) gene, were exchanged with that of the other paralog. Among mutants harboring a single native or chimeric RPL7 allele, expression from the RPL7A locus exceeded that from the RPL7B locus, and more Rpl7a was expressed from either locus than Rpl7b. Phenotypic differences in tunicamycin sensitivity, ASH1 mRNA localization, and mobility of the Ty1 retrotransposon were strongly correlated with Rpl7 and ribosome levels, but not with the Rpl7 or snoRNA isoform expressed. Although Ty1 RNA is cotranslationally localized, depletion of Rpl7 minimally affected synthesis of Ty1 Gag protein, but strongly influenced Ty1 RNA localization. Unlike the other processes studied, Ty1 cDNA accumulation was influenced by both the level and isoform of Rpl7 or snoRNA expressed. These cellular processes had different minimal threshold values for Rpl7 and ribosome levels, but all were functional when isoforms of either paralog were expressed from the RPL7A locus or both RPL7 loci. This study illustrates the broad range of phenotypes that can result from depleting ribosomes to different levels.

Increased rates of protein evolution and asymmetric deceleration after the whole-genome duplication in yeasts

Background

Whole-genome duplication (WGD) events have shaped the genomes of eukaryotic organisms. Relaxed selection after duplication along with inherent functional constraints are thought to determine the fate of the paralogs and, ultimately, the evolution of gene function. Here, we investigated the rate of protein evolution (as measured by dN/dS ratios) before and after the WGD in the hemiascomycete yeasts, and the way in which changes in such rates relate to molecular and biological function.

Results

For most groups of orthologous genes (81%) we observed a change in the rates of evolution after genome duplication. Genes with atypically-low dN/dS ratio before the WGD were prone to increase their rates of evolution after duplication. Importantly, the paralogs were often different in their rates of evolution after the WGD (50% cases), however, this was more consistent with an asymmetric deceleration in the protein-evolution rates, rather than an asymmetric increase of the initial rates. Functional-category analysis showed that regulatory proteins such as protein kinases and transcription factors were enriched in genes that increase their rates of evolution after the WGD. While changes in the rate of protein-sequence evolution were associated to protein abundance, content of disordered regions, and contribution to fitness, these features were an attribute of specific functional classes.

Conclusions

Our results indicate that strong purifying selection in ancestral pre-duplication sequences is a strong predictor of increased rates after the duplication in yeasts and that asymmetry in evolution rate is established during the deceleration phase. In addition, changes in the rates at which paralogous sequences evolve before and after WGD are different for specific protein functions; increased rates of protein evolution after duplication occur preferentially in specific protein functions.

Electronic supplementary material

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

The Phenotypic Plasticity of Duplicated Genes in Saccharomyces cerevisiae and the Origin of Adaptations

Gene and genome duplication are the major sources of biological innovations in plants and animals. Functional and transcriptional divergence between the copies after gene duplication has been considered the main driver of innovations . However, here we show that increased phenotypic plasticity after duplication plays a more major role than thought before in the origin of adaptations. We perform an exhaustive analysis of the transcriptional alterations of duplicated genes in the unicellular eukaryote Saccharomyces cerevisiae when challenged with five different environmental stresses. Analysis of the transcriptomes of yeast shows that gene duplication increases the transcriptional response to environmental changes, with duplicated genes exhibiting signatures of adaptive transcriptional patterns in response to stress. The mechanism of duplication matters, with whole-genome duplicates being more transcriptionally altered than small-scale duplicates. The predominant transcriptional pattern follows the classic theory of evolution by gene duplication; with one gene copy remaining unaltered under stress, while its sister copy presents large transcriptional plasticity and a prominent role in adaptation. Moreover, we find additional transcriptional profiles that are suggestive of neo- and subfunctionalization of duplicate gene copies. These patterns are strongly correlated with the functional dependencies and sequence divergence profiles of gene copies. We show that, unlike singletons, duplicates respond more specifically to stress, supporting the role of natural selection in the transcriptional plasticity of duplicates. Our results reveal the underlying transcriptional complexity of duplicated genes and its role in the origin of adaptations.

Computational Methods for Annotation Transfers from Sequence

Surveys of public sequence resources show that experimentally supported functional information is still completely missing for a considerable fraction of known proteins and is clearly incomplete for an even larger portion. Bioinformatics methods have long made use of very diverse data sources alone or in combination to predict protein function, with the understanding that different data types help elucidate complementary biological roles. This chapter focuses on methods accepting amino acid sequences as input and producing GO term assignments directly as outputs; the relevant biological and computational concepts are presented along with the advantages and limitations of individual approaches.

The RICE MINUTE-LIKE1 (RML1) gene, encoding a ribosomal large subunit protein L3B, regulates leaf morphology and plant architecture in rice

Mutations of ribosomal proteins (RPs) are known to cause developmental abnormalities in yeast, mammals, and dicotyledonous plants; however, their effects have not been studied in rice. Here, we identifiy a ribosomal biogenesis mutant, rice minute-like1 (rml1) that displays a minute phenotype as evidenced by retarded growth and defects in the vascular system. We determine that RML1 encodes a ribosome large subunit protein 3B (RPL3B) in rice by means of map-based cloning and genetic complementation. RPL3B is abundantly expressed in all the tissues, whereas RPL3A, another RPL3 gene family member, is expressed at low levels. Notably, the expression level of RPL3A in the rml1 mutant is similar to that in the wild-type, suggesting that RPL3A provides no functional compensation for RPL3B in rml1 plants. Ribosomal profiles show that mutation of RPL3B leads to a significant reduction in free 60S ribosomal subunits and polysomes, indicating a ribosomal insufficiency in the rml1 mutant. Our results demonstrate that the ribosomal protein gene RPL3B is required for maintaining normal leaf morphology and plant architecture in rice through its regulation of ribosome biogenesis.

Gene Duplicability of Core Genes Is Highly Consistent across All Angiosperms

Gene duplication is an important mechanism for adding to genomic novelty. Hence, which genes undergo duplication and are preserved following duplication is an important question. It has been observed that gene duplicability, or the ability of genes to be retained following duplication, is a nonrandom process, with certain genes being more amenable to survive duplication events than others. Primarily, gene essentiality and the type of duplication (small-scale versus large-scale) have been shown in different species to influence the (long-term) survival of novel genes. However, an overarching view of "gene duplicability" is lacking, mainly due to the fact that previous studies usually focused on individual species and did not account for the influence of genomic context and the time of duplication. Here, we present a large-scale study in which we investigated duplicate retention for 9178 gene families shared between 37 flowering plant species, referred to as angiosperm core gene families. For most gene families, we observe a strikingly consistent pattern of gene duplicability across species, with gene families being either primarily single-copy or multicopy in all species. An intermediate class contains gene families that are often retained in duplicate for periods extending to tens of millions of years after whole-genome duplication, but ultimately appear to be largely restored to singleton status, suggesting that these genes may be dosage balance sensitive. The distinction between single-copy and multicopy gene families is reflected in their functional annotation, with single-copy genes being mainly involved in the maintenance of genome stability and organelle function and multicopy genes in signaling, transport, and metabolism. The intermediate class was overrepresented in regulatory genes, further suggesting that these represent putative dosage-balance-sensitive genes.

Global analysis of human duplicated genes reveals the relative importance of whole-genome duplicates originated in the early vertebrate evolution

Background

Gene duplication is a genetic mutation that creates functionally redundant gene copies that are initially relieved from selective pressures and may adapt themselves to new functions with time. The levels of gene duplication may vary from small-scale duplication (SSD) to whole genome duplication (WGD). Studies with yeast revealed ample differences between these duplicates: Yeast WGD pairs were functionally more similar, less divergent in subcellular localization and contained a lesser proportion of essential genes. In this study, we explored the differences in evolutionary genomic properties of human SSD and WGD genes, with the identifiable human duplicates coming from the two rounds of whole genome duplication occurred early in vertebrate evolution.

Results

We observed that these two groups of duplicates were also dissimilar in terms of their evolutionary and genomic properties. But interestingly, this is not like the same observed in yeast. The human WGDs were found to be functionally less similar, diverge more in subcellular level and contain a higher proportion of essential genes than the SSDs, all of which are opposite from yeast. Additionally, we explored that human WGDs were more divergent in their gene expression profile, have higher multifunctionality and are more often associated with disease, and are evolutionarily more conserved than human SSDs.

Conclusions

Our study suggests that human WGD duplicates are more divergent and entails the adaptation of WGDs to novel and important functions that consequently lead to their evolutionary conservation in the course of evolution.

Electronic supplementary material

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

On the Relationship of Protein and mRNA Dynamics in Vertebrate Embryonic Development

A biochemical explanation of development from the fertilized egg to the adult requires an understanding of the proteins and RNAs expressed over time during embryogenesis. We present a comprehensive characterization of protein and mRNA dynamics across early development in Xenopus. Surprisingly, we find that most protein levels change little and duplicated genes are expressed similarly. While the correlation between protein and mRNA levels is poor, a mass action kinetics model parameterized using protein synthesis and degradation rates regresses protein dynamics to RNA dynamics, corrected for initial protein concentration. This study provides detailed data for absolute levels of ∼10,000 proteins and ∼28,000 transcripts via a convenient web portal, a rich resource for developmental biologists. It underscores the lasting impact of maternal dowry, finds surprisingly few cases where degradation alone drives a change in protein level, and highlights the importance of transcription in shaping the dynamics of the embryonic proteome.

Iron-Sulfur Cluster Biogenesis Chaperones: Evidence for Emergence of Mutational Robustness of a Highly Specific Protein-Protein Interaction

Biogenesis of iron-sulfur clusters (FeS) is a highly conserved process involving Hsp70 and J-protein chaperones. However, Hsp70 specialization differs among species. In most eukaryotes, including Schizosaccharomyces pombe, FeS biogenesis involves interaction between the J-protein Jac1 and the multifunctional Hsp70 Ssc1. But, in Saccharomyces cerevisiae and closely related species, Jac1 interacts with the specialized Hsp70 Ssq1, which emerged through duplication of SSC1. As little is known about how gene duplicates affect the robustness of their protein interaction partners, we analyzed the functional and evolutionary consequences of Ssq1 specialization on the ubiquitous J-protein cochaperone Jac1, by comparing S. cerevisiae and S. pombe. Although deletion of JAC1 is lethal in both species, alanine substitutions within the conserved His-Pro-Asp (HPD) motif, which is critical for Jac1:Hsp70 interaction, have species-specific effects. They are lethal in S. pombe, but not in S. cerevisiae. These in vivo differences correlated with in vitro biochemical measurements. Charged residues present in the J-domain of S. cerevisiae Jac1, but absent in S. pombe Jac1, are important for tolerance of S. cerevisiae Jac1 to HPD alterations. Moreover, Jac1 orthologs from species that encode Ssq1 have a higher sequence divergence. The simplest interpretation of our results is that Ssq1's coevolution with Jac1 resulted in expansion of their binding interface, thus increasing the efficiency of their interaction. Such an expansion could in turn compensate for negative effects of HPD substitutions. Thus, our results support the idea that the robustness of Jac1 emerged as consequence of its highly efficient and specific interaction with Ssq1.

Distinct Regulatory Changes Underlying Differential Expression of TEOSINTE BRANCHED1-CYCLOIDEA-PROLIFERATING CELL FACTOR Genes Associated with Petal Variations in Zygomorphic Flowers of Petrocosmea spp. of the Family Gesneriaceae

CYCLOIDEA (CYC)-like genes, belonging to the plant-specific TCP transcription factor family that is named after TEOSINTE BRANCHED1 (TB1) from maize (Zea mays), CYC from Antirrhinum majus, and the PROLIFERATING CELL FACTORS (PCF) from rice (Oryza sativa), have conserved dorsal identity function in patterning floral zygomorphy mainly through specific expression in dorsal petals of a flower. Their expression changes are usually related to morphological diversity of zygomorphic flowers. However, it is still a challenge to elucidate the molecular mechanism underlying their expression differentiation. It is also unknown whether CINCINNATA (CIN)-like TCP genes, locally controlling cell growth and proliferation, are involved in the evolution of floral zygomorphy. To address these questions, we selected two closely related species, i.e. Petrocosmea glabristoma and Petrocosmea sinensis, with distinct petal morphology to conduct expression, hybridization, mutant, and allele-specific expression analyses. The results show that the size change of the dorsal petals between the two species is mainly mediated by the expression differentiation of CYC1C and CYC1D, while the shape variation of all petals is related to the expression change of CIN1. In reciprocal F1 hybrids, the expression of CYC1C, CYC1D, and CIN1 conforms to an additive inheritance mode, consistent with the petal phenotypes of hybrids. Through allele-specific expression analyses, we find that the expression differentiation of these TCP genes is underlain by distinctly different types of regulatory changes. We suggest that highly redundant paralogs with identical expression patterns and interspecific expression differentiation may be controlled by remarkably different regulatory pathways because natural selection may favor different regulatory modifications rather than coding sequence changes of key developmental genes in generating morphological diversity.

Complete biosynthesis of opioids in yeast

Opioids are the primary drugs used in Western medicine for pain management and palliative care. Farming of opium poppies remains the sole source of these essential medicines, despite diverse market demands and uncertainty in crop yields due to weather, climate change, and pests. We engineered yeast to produce the selected opioid compounds thebaine and hydrocodone starting from sugar. All work was conducted in a laboratory that is permitted and secured for work with controlled substances. We combined enzyme discovery, enzyme engineering, and pathway and strain optimization to realize full opiate biosynthesis in yeast. The resulting opioid biosynthesis strains required the expression of 21 (thebaine) and 23 (hydrocodone) enzyme activities from plants, mammals, bacteria, and yeast itself. This is a proof of principle, and major hurdles remain before optimization and scale-up could be achieved. Open discussions of options for governing this technology are also needed in order to responsibly realize alternative supplies for these medically relevant compounds.

A Minimal Set of Glycolytic Genes Reveals Strong Redundancies in Saccharomyces cerevisiae Central Metabolism

As a result of ancestral whole-genome and small-scale duplication events, the genomes of Saccharomyces cerevisiae and many eukaryotes still contain a substantial fraction of duplicated genes. In all investigated organisms, metabolic pathways, and more particularly glycolysis, are specifically enriched for functionally redundant paralogs. In ancestors of the Saccharomyces lineage, the duplication of glycolytic genes is purported to have played an important role leading to S. cerevisiae's current lifestyle favoring fermentative metabolism even in the presence of oxygen and characterized by a high glycolytic capacity. In modern S. cerevisiae strains, the 12 glycolytic reactions leading to the biochemical conversion from glucose to ethanol are encoded by 27 paralogs. In order to experimentally explore the physiological role of this genetic redundancy, a yeast strain with a minimal set of 14 paralogs was constructed (the "minimal glycolysis" [MG] strain). Remarkably, a combination of a quantitative systems approach and semiquantitative analysis in a wide array of growth environments revealed the absence of a phenotypic response to the cumulative deletion of 13 glycolytic paralogs. This observation indicates that duplication of glycolytic genes is not a prerequisite for achieving the high glycolytic fluxes and fermentative capacities that are characteristic of S. cerevisiae and essential for many of its industrial applications and argues against gene dosage effects as a means of fixing minor glycolytic paralogs in the yeast genome. The MG strain was carefully designed and constructed to provide a robust prototrophic platform for quantitative studies and has been made available to the scientific community.

Gene duplication of type-B ARR transcription factors systematically extends transcriptional regulatory structures in Arabidopsis

Many of duplicated genes are enriched in signaling pathways. Recently, gene duplication of kinases has been shown to provide genetic buffering and functional diversification in cellular signaling. Transcription factors (TFs) are also often duplicated. However, how duplication of TFs affects their regulatory structures and functions of target genes has not been explored at the systems level. Here, we examined regulatory and functional roles of duplication of three major ARR TFs (ARR1, 10, and 12) in Arabidopsis cytokinin signaling using wild-type and single, double, and triple deletion mutants of the TFs. Comparative analysis of gene expression profiles obtained from Arabidopsis roots in wild-type and these mutants showed that duplication of ARR TFs systematically extended their transcriptional regulatory structures, leading to enhanced robustness and diversification in functions of target genes, as well as in regulation of cellular networks of target genes. Therefore, our results suggest that duplication of TFs contributes to robustness and diversification in functions of target genes by extending transcriptional regulatory structures.

Extensive heterosis in growth of yeast hybrids is explained by a combination of genetic models

Heterosis, also known as hybrid vigor, is the superior performance of a heterozygous hybrid relative to its homozygous parents. Despite the scientific curiosity of this phenotypic phenomenon and its significance for food production in agriculture, its genetic basis is insufficiently understood. Studying heterosis in yeast can potentially yield insights into its genetic basis, can allow one to test the different hypotheses that have been proposed to explain the phenomenon and allows better understanding of how to take advantage of this phenomenon to enhance food production. We therefore crossed 16 parental yeast strains to form 120 yeast hybrids, and measured their growth rates under five environmental conditions. A considerable amount of dominant genetic variation was found in growth performance, and heterosis was measured in 35% of the hybrid-condition combinations. Despite previous reports of correlations between heterosis and measures of sequence divergence between parents, we detected no such relationship. We used several analyses to examine which genetic model might explain heterosis. We found that dominance complementation of recessive alleles, overdominant interactions within loci and epistatic interactions among loci each contribute to heterosis. We concluded that in yeast heterosis is a complex phenotype created by the combined contribution of different genetic interactions.

Preservation of genetic and regulatory robustness in ancient gene duplicates of Saccharomyces cerevisiae

Biological systems remain robust against certain genetic and environmental challenges. Robustness allows the exploration of ecological adaptations. It is unclear what factors contribute to increasing robustness. Gene duplication has been considered to increase genetic robustness through functional redundancy, accelerating the evolution of novel functions. However, recent findings have questioned the link between duplication and robustness. In particular, it remains elusive whether ancient duplicates still bear potential for innovation through preserved redundancy and robustness. Here we have investigated this question by evolving the yeast Saccharomyces cerevisiae for 2200 generations under conditions allowing the accumulation of deleterious mutations, and we put mechanisms of mutational robustness to a test. S. cerevisiae declined in fitness along the evolution experiment, but this decline decelerated in later passages, suggesting functional compensation of mutated genes. We resequenced 28 genomes from experimentally evolved S. cerevisiae lines and found more mutations in duplicates—mainly small-scale duplicates—than in singletons. Genetically interacting duplicates evolved similarly and fixed more amino acid–replacing mutations than expected. Regulatory robustness of the duplicates was supported by a larger enrichment for mutations at the promoters of duplicates than at those of singletons. Analyses of yeast gene expression conditions showed a larger variation in the duplicates’ expression than that of singletons under a range of stress conditions, sparking the idea that regulatory robustness allowed a wider range of phenotypic responses to environmental stresses, hence faster adaptations. Our data support the persistence of genetic and regulatory robustness in ancient duplicates and its role in adaptations to stresses.

Effect of duplicate genes on mouse genetic robustness: an update

In contrast to S. cerevisiae and C. elegans, analyses based on the current knockout (KO) mouse phenotypes led to the conclusion that duplicate genes had almost no role in mouse genetic robustness. It has been suggested that the bias of mouse KO database toward ancient duplicates may possibly cause this knockout duplicate puzzle, that is, a very similar proportion of essential genes (P_E) between duplicate genes and singletons. In this paper, we conducted an extensive and careful analysis for the mouse KO phenotype data and corroborated a strong effect of duplicate genes on mouse genetics robustness. Moreover, the effect of duplicate genes on mouse genetic robustness is duplication-age dependent, which holds after ruling out the potential confounding effect from coding-sequence conservation, protein-protein connectivity, functional bias, or the bias of duplicates generated by whole genome duplication (WGD). Our findings suggest that two factors, the sampling bias toward ancient duplicates and very ancient duplicates with a proportion of essential genes higher than that of singletons, have caused the mouse knockout duplicate puzzle; meanwhile, the effect of genetic buffering may be correlated with sequence conservation as well as protein-protein interactivity.

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.

Genomic evidence for adaptation by gene duplication

Gene duplication is widely believed to facilitate adaptation, but unambiguous evidence for this hypothesis has been found in only a small number of cases. Although gene duplication may increase the fitness of the involved organisms by doubling gene dosage or neofunctionalization, it may also result in a simple division of ancestral functions into daughter genes, which need not promote adaptation. Hence, the general validity of the adaptation by gene duplication hypothesis remains uncertain. Indeed, a genome-scale experiment found similar fitness effects of deleting pairs of duplicate genes and deleting individual singleton genes from the yeast genome, leading to the conclusion that duplication rarely results in adaptation. Here we contend that the above comparison is unfair because of a known duplication bias among genes with different fitness contributions. To rectify this problem, we compare homologous genes from the budding yeast Saccharomyces cerevisiae and the fission yeast Schizosaccharomyces pombe. We discover that simultaneously deleting a duplicate gene pair in S. cerevisiae reduces fitness significantly more than deleting their singleton counterpart in S. pombe, revealing post-duplication adaptation. The duplicates-singleton difference in fitness effect is not attributable to a potential increase in gene dose after duplication, suggesting that the adaptation is owing to neofunctionalization, which we find to be explicable by acquisitions of binary protein-protein interactions rather than gene expression changes. These results provide genomic evidence for the role of gene duplication in organismal adaptation and are important for understanding the genetic mechanisms of evolutionary innovation.

Expression evolution facilitated the convergent neofunctionalization of a sodium channel gene

Ion channels have played a substantial role in the evolution of novel traits across all of the domains of life. A fascinating example of a novel adaptation is the convergent evolution of electric organs in the Mormyroid and Gymnotiform electric fishes. The regulated currents that flow through ion channels directly generate the electrical signals which have evolved in these fish. Here, we investigated how the expression evolution of two sodium channel paralogs (Scn4aa and Scn4ab) influenced their convergent molecular evolution following the teleost-specific whole-genome duplication. We developed a reliable assay to accurately measure the expression stoichiometry of these genes and used this technique to analyze relative expression of the duplicate genes in a phylogenetic context. We found that before a major shift in expression from skeletal muscle and neofunctionalization in the muscle-derived electric organ, Scn4aa was first downregulated in the ancestors of both electric lineages. This indicates that underlying the convergent evolution of this gene, there was a greater propensity toward neofunctionalization due to its decreased expression relative to its paralog Scn4ab. We investigated another derived muscle tissue, the sonic organ of Porichthys notatus, and show that, as in the electric fishes, Scn4aa again shows a radical shift in expression away from the ancestral muscle cells into the evolutionarily novel muscle-derived tissue. This study presents evidence that expression downregulation facilitates neofunctionalization after gene duplication, a pattern that may often set the stage for novel trait evolution after gene duplication.

Contribution of nonohnologous duplicated genes to high habitat variability in mammals

The mechanism by which genetic systems affect environmental adaptation is a focus of considerable attention in the fields of ecology, evolution, and conservation. However, the genomic characteristics that constrain adaptive evolution have remained unknown. A recent study showed that the proportion of duplicated genes in whole Drosophila genomes correlated with environmental variability within habitat, but it remains unclear whether the correlation is observed even in vertebrates whose genomes including a large number of duplicated genes generated by whole-genome duplication (WGD). Here, we focus on fully sequenced mammalian genomes that experienced WGD in early vertebrate lineages and show that the proportion of small-scale duplication (SSD) genes in the genome, but not that of WGD genes, is significantly correlated with habitat variability. Moreover, species with low habitat variability have a higher proportion of lost duplicated genes, particularly SSD genes, than those with high habitat variability. These results indicate that species that inhabit variable environments may maintain more SSD genes in their genomes and suggest that SSD genes are important for adapting to novel environments and surviving environmental changes. These insights may be applied to predicting invasive and endangered species.

An insulin-to-insulin regulatory network orchestrates phenotypic specificity in development and physiology

Insulin-like peptides (ILPs) play highly conserved roles in development and physiology. Most animal genomes encode multiple ILPs. Here we identify mechanisms for how the forty Caenorhabditis elegans ILPs coordinate diverse processes, including development, reproduction, longevity and several specific stress responses. Our systematic studies identify an ILP-based combinatorial code for these phenotypes characterized by substantial functional specificity and diversity rather than global redundancy. Notably, we show that ILPs regulate each other transcriptionally, uncovering an ILP-to-ILP regulatory network that underlies the combinatorial phenotypic coding by the ILP family. Extensive analyses of genetic interactions among ILPs reveal how their signals are integrated. A combined analysis of these functional and regulatory ILP interactions identifies local genetic circuits that act in parallel and interact by crosstalk, feedback and compensation. This organization provides emergent mechanisms for phenotypic specificity and graded regulation for the combinatorial phenotypic coding we observe. Our findings also provide insights into how large hormonal networks regulate diverse traits.

Human gut microbes use multiple transporters to distinguish vitamin B₁₂ analogs and compete in the gut

Genomic and metagenomic sequencing efforts, including human microbiome projects, reveal that microbes often encode multiple systems that appear to accomplish the same task. Whether these predictions reflect actual functional redundancies is unclear. We report that the prominent human gut symbiont Bacteroides thetaiotaomicron employs three functional, homologous vitamin B₁₂ transporters that in at least two cases confer a competitive advantage in the presence of distinct B₁₂ analogs (corrinoids). In the mammalian gut, microbial fitness can be determined by the presence or absence of a single transporter. The total number of distinct corrinoid transporter families in the human gut microbiome likely exceeds those observed in B. thetaiotaomicron by an order of magnitude. These results demonstrate that human gut microbes use elaborate mechanisms to capture and differentiate corrinoids in vivo and that apparent redundancies observed in these genomes can instead reflect hidden specificities that determine whether a microbe will colonize its host.

All or nothing: protein complexes flip essentiality between distantly related eukaryotes

In the budding yeast Saccharomyces cerevisiae, the subunits of any given protein complex are either mostly essential or mostly nonessential, suggesting that essentiality is a property of molecular machines rather than individual components. There are exceptions to this rule, however, that is, nonessential genes in largely essential complexes and essential genes in largely nonessential complexes. Here, we provide explanations for these exceptions, showing that redundancy within complexes, as revealed by genetic interactions, can explain many of the former cases, whereas “moonlighting,” as revealed by membership of multiple complexes, can explain the latter. Surprisingly, we find that redundancy within complexes cannot usually be explained by gene duplication, suggesting alternate buffering mechanisms. In the distantly related Schizosaccharomyces pombe, we observe the same phenomenon of modular essentiality, suggesting that it may be a general feature of eukaryotes. Furthermore, we show that complexes flip essentiality in a cohesive fashion between the two species, that is, they tend to change from mostly essential to mostly nonessential, or vice versa, but not to mixed patterns. We show that these flips in essentiality can be explained by differing lifestyles of the two yeasts. Collectively, our results support a previously proposed model where proteins are essential because of their involvement in essential functional modules rather than because of specific topological features such as degree or centrality.

Molecular mechanisms of paralogous compensation and the robustness of cellular networks

Robustness is the ability of a system to maintain its function despite environmental or genetic perturbation. Genetic robustness is a key emerging property of living systems and is achieved notably by the presence of partially redundant parts that result from gene duplication. Functional overlap between paralogs allows them to compensate for each other's loss, as commonly revealed by aggravating genetic interactions. However, the molecular mechanisms linking the genotype (loss of function of a gene) to the phenotype (genetic buffering by a paralog) are still poorly understood and the molecular aspects of this compensation are rarely addressed in studies of gene duplicates. Here, we review molecular mechanisms of functional compensation between paralogous genes, many of which from studies that were not meant to study this phenomenon. We propose a standardized terminology and, depending on whether or not the molecular behavior of the intact gene is modified in response to the deletion of its paralog, we classify mechanisms of compensation into passive and active events. We further describe three non-exclusive mechanisms of active paralogous compensation for which there is evidence in the literature: changes in abundance, in localization, and in protein interactions. This review will serve as a framework for the genetic and molecular analysis of paralogous compensation, one of the universal features of genetic systems.

Systems biology and the integration of mechanistic explanation and mathematical explanation

The paper discusses how systems biology is working toward complex accounts that integrate explanation in terms of mechanisms and explanation by mathematical models-which some philosophers have viewed as rival models of explanation. Systems biology is an integrative approach, and it strongly relies on mathematical modeling. Philosophical accounts of mechanisms capture integrative in the sense of multilevel and multifield explanations, yet accounts of mechanistic explanation (as the analysis of a whole in terms of its structural parts and their qualitative interactions) have failed to address how a mathematical model could contribute to such explanations. I discuss how mathematical equations can be explanatorily relevant. Several cases from systems biology are discussed to illustrate the interplay between mechanistic research and mathematical modeling, and I point to questions about qualitative phenomena (rather than the explanation of quantitative details), where quantitative models are still indispensable to the explanation. Systems biology shows that a broader philosophical conception of mechanisms is needed, which takes into account functional-dynamical aspects, interaction in complex networks with feedback loops, system-wide functional properties such as distributed functionality and robustness, and a mechanism's ability to respond to perturbations (beyond its actual operation). I offer general conclusions for philosophical accounts of explanation.

Contemporary, yeast-based approaches to understanding human genetic variation

Determining how genetic variation contributes to human health and disease is a critical challenge. As one of the most genetically tractable model organisms, yeast has played a central role in meeting this challenge. The advent of new technologies, including high-throughput DNA sequencing and synthesis, proteomics, and computational methods, has vastly increased the power of yeast-based approaches to determine the consequences of human genetic variation. Recent successes include systematic exploration of the effects of gene dosage, large-scale analysis of the effect of coding variation on gene function, and the use of humanized yeast to model disease. By virtue of its manipulability, small genome size, and genetic tractability, yeast is poised to help us understand human genetic variation.

Evolution after whole-genome duplication: a network perspective

Gene duplication plays an important role in the evolution of genomes and interactomes. Elucidating how evolution after gene duplication interplays at the sequence and network level is of great interest. In this work, we analyze a data set of gene pairs that arose through whole-genome duplication (WGD) in yeast. All these pairs have the same duplication time, making them ideal for evolutionary investigation. We investigated the interplay between evolution after WGD at the sequence and network levels and correlated these two levels of divergence with gene expression and fitness data. We find that molecular interactions involving WGD genes evolve at rates that are three orders of magnitude slower than the rates of evolution of the corresponding sequences. Furthermore, we find that divergence of WGD pairs correlates strongly with gene expression and fitness data. Because of the role of gene duplication in determining redundancy in biological systems and particularly at the network level, we investigated the role of interaction networks in elucidating the evolutionary fate of duplicated genes. We find that gene neighborhoods in interaction networks provide a mechanism for inferring these fates, and we developed an algorithm for achieving this task. Further epistasis analysis of WGD pairs categorized by their inferred evolutionary fates demonstrated the utility of these techniques. Finally, we find that WGD pairs and other pairs of paralogous genes of small-scale duplication origin share similar properties, giving good support for generalizing our results from WGD pairs to evolution after gene duplication in general.

Different patterns of gene structure divergence following gene duplication in Arabidopsis

BACKGROUND

Divergence in gene structure following gene duplication is not well understood. Gene duplication can occur via whole-genome duplication (WGD) and single-gene duplications including tandem, proximal and transposed duplications. Different modes of gene duplication may be associated with different types, levels, and patterns of structural divergence.

RESULTS

In Arabidopsis thaliana, we denote levels of structural divergence between duplicated genes by differences in coding-region lengths and average exon lengths, and the number of insertions/deletions (indels) and maximum indel length in their protein sequence alignment. Among recent duplicates of different modes, transposed duplicates diverge most dramatically in gene structure. In transposed duplications, parental loci tend to have longer coding-regions and exons, and smaller numbers of indels and maximum indel lengths than transposed loci, reflecting biased structural changes in transposed duplications. Structural divergence increases with evolutionary time for WGDs, but not transposed duplications, possibly because of biased gene losses following transposed duplications. Structural divergence has heterogeneous relationships with nucleotide substitution rates, but is consistently positively correlated with gene expression divergence. The NBS-LRR gene family shows higher-than-average levels of structural divergence.

CONCLUSIONS

Our study suggests that structural divergence between duplicated genes is greatly affected by the mechanisms of gene duplication and may be not proportional to evolutionary time, and that certain gene families are under selection on rapid evolution of gene structure.

The ribosomal protein Rpl22 controls ribosome composition by directly repressing expression of its own paralog, Rpl22l1

Most yeast ribosomal protein genes are duplicated and their characterization has led to hypotheses regarding the existence of specialized ribosomes with different subunit composition or specifically-tailored functions. In yeast, ribosomal protein genes are generally duplicated and evidence has emerged that paralogs might have specific roles. Unlike yeast, most mammalian ribosomal proteins are thought to be encoded by a single gene copy, raising the possibility that heterogenous populations of ribosomes are unique to yeast. Here, we examine the roles of the mammalian Rpl22, finding that Rpl22^−/− mice have only subtle phenotypes with no significant translation defects. We find that in the Rpl22^−/− mouse there is a compensatory increase in Rpl22-like1 (Rpl22l1) expression and incorporation into ribosomes. Consistent with the hypothesis that either ribosomal protein can support translation, knockdown of Rpl22l1 impairs growth of cells lacking Rpl22. Mechanistically, Rpl22 regulates Rpl22l1 directly by binding to an internal hairpin structure and repressing its expression. We propose that ribosome specificity may exist in mammals, providing evidence that one ribosomal protein can influence composition of the ribosome by regulating its own paralog.

Translation is the process by which proteins are made within a cell. Ribosomes are the main macromolecular complexes involved in this process. Ribosomes are composed of ribosomal RNA and ribosomal proteins. Ribosomal proteins are generally thought to be structural components of the ribosome but recent findings have suggested that they might have a regulatory function as well. A growing number of human diseases have been linked to mutations in genes encoding factors involved in ribosome biogenesis and translation. These include developmental malformations, inherited bone marrow failure syndromes and cancer in a variety of organisms. Here, we describe the role of one ribosomal protein regulating another. We provide evidence that ribosomal proteins can influence the composition of the ribosome, which we hypothesize, may impact the function of the ribosome.

Specialized yeast ribosomes: a customized tool for selective mRNA translation

Evidence is now accumulating that sub-populations of ribosomes - so-called specialized ribosomes - can favour the translation of subsets of mRNAs. Here we use a large collection of diploid yeast strains, each deficient in one or other copy of the set of ribosomal protein (RP) genes, to generate eukaryotic cells carrying distinct populations of altered ‘specialized’ ribosomes. We show by comparative protein synthesis assays that different heterologous mRNA reporters based on luciferase are preferentially translated by distinct populations of specialized ribosomes. These mRNAs include reporters carrying premature termination codons (PTC) thus allowing us to identify specialized ribosomes that alter the efficiency of translation termination leading to enhanced synthesis of the wild-type protein. This finding suggests that these strains can be used to identify novel therapeutic targets in the ribosome. To explore this further we examined the translation of the mRNA encoding the extracellular matrix protein laminin β3 (LAMB3) since a LAMB3-PTC mutant is implicated in the blistering skin disease Epidermolysis bullosa (EB). This screen identified specialized ribosomes with reduced levels of RP L35B as showing enhanced synthesis of full-length LAMB3 in cells expressing the LAMB3-PTC mutant. Importantly, the RP L35B sub-population of specialized ribosomes leave both translation of a reporter luciferase carrying a different PTC and bulk mRNA translation largely unaltered.

A systematic approach for the genetic dissection of protein complexes in living cells

Cells contain many important protein complexes involved in performing and regulating structural, metabolic, and signaling functions. One major challenge in cell biology is to elucidate the organization and mechanisms of robustness of these complexes in vivo. We developed a systematic approach to study structural dependencies within complexes in living cells by deleting subunits and measuring pairwise interactions among other components. We used our methodology to perturb two conserved eukaryotic complexes: the retromer and the nuclear pore complex. Our results identify subunits that are critical for the assembly of these complexes, reveal their structural architecture, and uncover mechanisms by which protein interactions are modulated. Our results also show that paralogous proteins play a key role in the robustness of protein complexes and shape their assembly landscape. Our approach paves the way for studying the response of protein interactomes to mutations and enhances our understanding of genotype-phenotype maps.

Human monogenic disease genes have frequently functionally redundant paralogs

Mendelian disorders are often caused by mutations in genes that are not lethal but induce functional distortions leading to diseases. Here we study the extent of gene duplicates that might compensate genes causing monogenic diseases. We provide evidence for pervasive functional redundancy of human monogenic disease genes (MDs) by duplicates by manifesting 1) genes involved in human genetic disorders are enriched in duplicates and 2) duplicated disease genes tend to have higher functional similarities with their closest paralogs in contrast to duplicated non-disease genes of similar age. We propose that functional compensation by duplication of genes masks the phenotypic effects of deleterious mutations and reduces the probability of purging the defective genes from the human population; this functional compensation could be further enhanced by higher purification selection between disease genes and their duplicates as well as their orthologous counterpart compared to non-disease genes. However, due to the intrinsic expression stochasticity among individuals, the deleterious mutations could still be present as genetic diseases in some subpopulations where the duplicate copies are expressed at low abundances. Consequently the defective genes are linked to genetic disorders while they continue propagating within the population. Our results provide insight into the molecular basis underlying the spreading of duplicated disease genes.

Duplication and retention biases of essential and non-essential genes revealed by systematic knockdown analyses

When a duplicate gene has no apparent loss-of-function phenotype, it is commonly considered that the phenotype has been masked as a result of functional redundancy with the remaining paralog. This is supported by indirect evidence showing that multi-copy genes show loss-of-function phenotypes less often than single-copy genes and by direct tests of phenotype masking using select gene sets. Here we take a systematic genome-wide RNA interference approach to assess phenotype masking in paralog pairs in the Caenorhabditis elegans genome. Remarkably, in contrast to expectations, we find that phenotype masking makes only a minor contribution to the low knockdown phenotype rate for duplicate genes. Instead, we find that non-essential genes are highly over-represented among duplicates, leading to a low observed loss-of-function phenotype rate. We further find that duplicate pairs derived from essential and non-essential genes have contrasting evolutionary dynamics: whereas non-essential genes are both more often successfully duplicated (fixed) and lost, essential genes are less often duplicated but upon successful duplication are maintained over longer periods. We expect the fundamental evolutionary duplication dynamics presented here to be broadly applicable.

Duplicate genes occur in all organisms. It has been found that mutations in duplicate genes cause defects much less often than when single copy genes are mutated. It is widely believed that this is due to functional redundancy—that is, the two genes can carry out similar functions so that the non-mutated duplicate gene can cover for or “mask” the phenotype of the mutation in the first duplicate. To determine whether this hypothesis is true, it is necessary to test systematically whether defects indeed occur in the organism when both duplicate genes are inhibited. We have for the first time carried out such an analysis in a multicellular organism, the nematode Caenorhabditis elegans. In contrast to expectations, we observed that when both copies of duplicate genes are inhibited deleterious effects are very rare. We show that this is because duplicate genes are much more often non-essential compared to genes where there is only a single copy. Non-essential genes are also lost from the genome much more often than essential genes. However, when essential genes are duplicated, they remain present in the genome over longer periods. Our results give a framework to explain the evolutionary dynamics of duplications in the genome.

The role of protein interactions in mediating essentiality and synthetic lethality

Genes are characterized as essential if their knockout is associated with a lethal phenotype, and these "essential genes" play a central role in biological function. In addition, some genes are only essential when deleted in pairs, a phenomenon known as synthetic lethality. Here we consider genes displaying synthetic lethality as "essential pairs" of genes, and analyze the properties of yeast essential genes and synthetic lethal pairs together. As gene duplication initially produces an identical pair or sets of genes, it is often invoked as an explanation for synthetic lethality. However, we find that duplication explains only a minority of cases of synthetic lethality. Similarly, disruption of metabolic pathways leads to relatively few examples of synthetic lethality. By contrast, the vast majority of synthetic lethal gene pairs code for proteins with related functions that share interaction partners. We also find that essential genes and synthetic lethal pairs cluster in the protein-protein interaction network. These results suggest that synthetic lethality is strongly dependent on the formation of protein-protein interactions. Compensation by duplicates does not usually occur mainly because the genes involved are recent duplicates, but is more commonly due to functional similarity that permits preservation of essential protein complexes. This unified view, combining genes that are individually essential with those that form essential pairs, suggests that essentiality is a feature of physical interactions between proteins protein-protein interactions, rather than being inherent in gene and protein products themselves.

The protein chaperone HSP90 can facilitate the divergence of gene duplicates

The heat-shock protein 90 (HSP90) acts as a chaperone by ensuring proper maturation and folding of its client proteins. The HSP90 capacitor hypothesis holds that interactions with HSP90 allow proteins to accumulate mutations while maintaining function. Following this logic, HSP90 clients would be predicted to show relaxed selection compared with nonclients. In this study, we identify a new HSP90 client in the plant steroid hormone pathway: the transcription factor BES1. Its closest paralog, BZR1, is not an HSP90 client. This difference in HSP90 client status in two highly similar proteins enabled a direct test of the capacitor hypothesis. We find that BES1 shows relaxed selection compared to BZR1, hallmarks of neo- and subfunctionalization, and dynamic HSP90 client status across independent evolutionary paths. These results suggested that HSP90's influence on gene evolution may be detectable if we compare gene duplicates because duplicates share most other properties influencing evolutionary rate that might otherwise conceal the chaperone's effect. We test this hypothesis using systematically identified HSP90 clients in yeast and observe a significant trend of HSP90 clients evolving faster than their nonclient paralogs. This trend was not detected when yeast clients and nonclients were compared without considering paralog status. Our data provide evidence that HSP90 influences selection on genes encoding its clients and facilitates divergence between 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.

Gene functional trade-offs and the evolution of pleiotropy

Pleiotropy is the property of genes affecting multiple functions or characters of an organism. Genes vary widely in their degree of pleiotropy, but this variation is often considered a by-product of their evolutionary history. We present a functional theory of how pleiotropy may itself evolve. We consider genes that contribute to two functions, where contributing more to one function detracts from allocation to the second function. We show that whether genes become pleiotropic or specialize on a single function depends on the nature of trade-offs as gene activities contribute to different traits and on how the functionality of these traits affects fitness. In general, when a gene product can perform well at two functions, it evolves to do so, but not when pleiotropy would greatly disrupt each function. Consequently, reduced pleiotropy should often evolve, with genes specializing on the trait that is currently more important to fitness. Even when pleiotropy does evolve, not all genes are expected to become equally pleiotropic; genes with higher levels of expression are more likely to evolve greater pleiotropy. For the case of gene duplicates, we find that perfect subfunctionalization evolves only under stringent conditions. More often, duplicates are expected to maintain a certain degree of functional redundancy, with the gene contributing more to trait functionality evolving the highest degree of pleiotropy. Gene product interactions can facilitate subfunctionalization, but whether they do so depends on the curvature of the fitness surface. Finally, we find that stochastic gene expression favors pleiotropy by selecting for robustness in fitness components.

Feature Identification of Compensatory Gene Pairs without Sequence Homology in Yeast

Genetic robustness refers to a compensatory mechanism for buffering deleterious mutations or environmental variations. Gene duplication has been shown to provide such functional backups. However, the overall contribution of duplication-based buffering for genetic robustness is rather small. In this study, we investigated whether transcriptional compensation also exists among genes that share similar functions without sequence homology. A set of nonhomologous synthetic-lethal gene pairs was assessed by using a coexpression network, protein-protein interactions, and other types of genetic interactions in yeast. Our results are notably different from those of previous studies on buffering paralogs. The low expression similarity and the conditional coexpression alone do not play roles in identifying the functionally compensatory genes. Additional properties such as synthetic-lethal interaction, the ratio of shared common interacting partners, and the degree of coregulation were, at least in part, necessary to extract functional compensatory genes. Our network-based approach is applicable to select several well-documented cases of compensatory gene pairs and a set of new pairs. The results suggest that transcriptional reprogramming plays a limited role in functional compensation among nonhomologous genes. Our study aids in understanding the mechanism and features of functional compensation more in detail.