Guihong Tan

A global genetic interaction network maps a wiring diagram of cellular function

INTRODUCTION Genetic interactions occur when mutations in two or more genes combine to generate an unexpected phenotype. An extreme negative or synthetic lethal genetic interaction occurs when two mutations, neither lethal individually, combine to cause cell death. Conversely, positive genetic interactions occur when two mutations produce a phenotype that is less severe than expected. Genetic interactions identify functional relationships between genes and can be harnessed for biological discovery and therapeutic target identification. They may also explain a considerable component of the undiscovered genetics associated with human diseases. Here, we describe construction and analysis of a comprehensive genetic interaction network for a eukaryotic cell. RATIONALE Genome sequencing projects are providing an unprecedented view of genetic variation. However, our ability to interpret genetic information to predict inherited phenotypes remains limited, in large part due to the extensive buffering of genomes, making most individual eukaryotic genes dispensable for life. To explore the extent to which genetic interactions reveal cellular function and contribute to complex phenotypes, and to discover the general principles of genetic networks, we used automated yeast genetics to construct a global genetic interaction network. RESULTS We tested most of the ~6000 genes in the yeast Saccharomyces cerevisiae for all possible pairwise genetic interactions, identifying nearly 1 million interactions, including ~550,000 negative and ~350,000 positive interactions, spanning ~90% of all yeast genes. Essential genes were network hubs, displaying five times as many interactions as nonessential genes. The set of genetic interactions or the genetic interaction profile for a gene provides a quantitative measure of function, and a global network based on genetic interaction profile similarity revealed a hierarchy of modules reflecting the functional architecture of a cell. Negative interactions connected functionally related genes, mapped core bioprocesses, and identified pleiotropic genes, whereas positive interactions often mapped general regulatory connections associated with defects in cell cycle progression or cellular proteostasis. Importantly, the global network illustrates how coherent sets of negative or positive genetic interactions connect protein complex and pathways to map a functional wiring diagram of the cell. CONCLUSION A global genetic interaction network highlights the functional organization of a cell and provides a resource for predicting gene and pathway function. This network emphasizes the prevalence of genetic interactions and their potential to compound phenotypes associated with single mutations. Negative genetic interactions tend to connect functionally related genes and thus may be predicted using alternative functional information. Although less functionally informative, positive interactions may provide insights into general mechanisms of genetic suppression or resiliency. We anticipate that the ordered topology of the global genetic network, in which genetic interactions connect coherently within and between protein complexes and pathways, may be exploited to decipher genotype-to-phenotype relationships. A global network of genetic interaction profile similarities. (Left) Genes with similar genetic interaction profiles are connected in a global network, such that genes exhibiting more similar profiles are located closer to each other, whereas genes with less similar profiles are positioned farther apart. (Right) Spatial analysis of functional enrichment was used to identify and color network regions enriched for similar Gene Ontology bioprocess terms. We generated a global genetic interaction network for Saccharomyces cerevisiae, constructing more than 23 million double mutants, identifying about 550,000 negative and about 350,000 positive genetic interactions. This comprehensive network maps genetic interactions for essential gene pairs, highlighting essential genes as densely connected hubs. Genetic interaction profiles enabled assembly of a hierarchical model of cell function, including modules corresponding to protein complexes and pathways, biological processes, and cellular compartments. Negative interactions connected functionally related genes, mapped core bioprocesses, and identified pleiotropic genes, whereas positive interactions often mapped general regulatory connections among gene pairs, rather than shared functionality. The global network illustrates how coherent sets of genetic interactions connect protein complex and pathway modules to map a functional wiring diagram of the cell.

An extended set of yeast-based functional assays accurately identifies human disease mutations

We can now routinely identify coding variants within individual human genomes. A pressing challenge is to determine which variants disrupt the function of disease-associated genes. Both experimental and computational methods exist to predict pathogenicity of human genetic variation. However, a systematic performance comparison between them has been lacking. Therefore, we developed and exploited a panel of 26 yeast-based functional complementation assays to measure the impact of 179 variants (101 disease- and 78 non-disease-associated variants) from 22 human disease genes. Using the resulting reference standard, we show that experimental functional assays in a 1-billion-year diverged model organism can identify pathogenic alleles with significantly higher precision and specificity than current computational methods.

Systematic analysis of complex genetic interactions

Trigenic interactions in yeast link bioprocesses To dissect the genotype-phenotype landscape of a cell, it is necessary to understand interactions between genes. Building on the digenic protein-protein interaction network, Kuzmin et al. created a trigenic landscape of yeast by using a synthetic genetic array (see the Perspective by Walhout). Triple-mutant analyses indicated that the majority of genes with trigenic associations functioned within the same biological processes. These converged on networks identified in the digenic interaction landscape. Although the overall effects were weaker for trigenic than for digenic interactions, trigenic interactions were more likely to bridge biological processes in the cell. Science, this issue p. eaao1729; see also p. 269 Trigenic interactions in yeast link bioprocesses are explored. INTRODUCTION Genetic interactions occur when mutations in different genes combine to result in a phenotype that is different from expectation based on those of the individual mutations. Negative genetic interactions occur when a combination of mutations leads to a fitness defect that is more exacerbated than expected. For example, synthetic lethality occurs when two mutations, neither of which is lethal on its own, generate an inviable double mutant. Alternatively, positive genetic interactions occur when genetic perturbations combine to generate a double mutant with a greater fitness than expected. Global digenic interaction studies have been useful for understanding the functional wiring diagram of the cell and may also provide insight into the genotype-to-phenotype relationship, which is important for tracking the missing heritability of human health and disease. Here we describe a network of higher-order trigenic interactions and explore its implications. RATIONALE Variation in phenotypic outcomes in different individuals is caused by genetic determinants that act as modifiers. Modifier loci are prevalent in human populations, but knowledge regarding how variants interact to modulate phenotype in different individuals is lacking. Similarly, in yeast, traits including conditional essentiality—in which certain genes are essential in one genetic background but nonessential in another—often result from an interplay of multiple modifier loci. Because complex modifiers may underlie the genetic basis of physiological states found in natural populations, it is critical to understand the landscape of higher-order genetic interactions. RESULTS To survey trigenic interactions, we designed query strains that sampled key features of the global digenic interaction network: (i) digenic interaction strength, (ii) average number of digenic interactions, and (iii) digenic interaction profile similarity. In total, we tested ~400,000 double and ~200,000 triple mutants for fitness defects and identified ~9500 digenic and ~3200 trigenic negative interactions. Although trigenic interactions tend to be weaker than digenic interactions, they were both enriched for functional relationships. About one-third of trigenic interactions identified “novel” connections that were not observed in our digenic control network, whereas the remaining approximately two-thirds of trigenic interactions “modified” a digenic interaction, suggesting that the global digenic interaction network is important for understanding the trigenic interaction network. Despite their functional enrichment, trigenic interactions also bridged distant bioprocesses. We estimate that the global trigenic interaction network is ~100 times as large as the global digenic network, highlighting the potential for complex genetic interactions to affect the biology of inheritance. CONCLUSION The extensive network of trigenic interactions and their ability to generate functionally diverse phenotypes suggest that higher-order genetic interactions may play a key role in the genotype-to-phenotype relationship, genome size, and speciation. Systematic analysis of trigenic interactions. We surveyed for trigenic interactions and found that they are ~100 times as prevalent as digenic interactions, often modify a digenic interaction, and connect functionally related genes as well as genes in more diverse bioprocesses (multicolored nodes). PPI, protein-protein interaction. To systematically explore complex genetic interactions, we constructed ~200,000 yeast triple mutants and scored negative trigenic interactions. We selected double-mutant query genes across a broad spectrum of biological processes, spanning a range of quantitative features of the global digenic interaction network and tested for a genetic interaction with a third mutation. Trigenic interactions often occurred among functionally related genes, and essential genes were hubs on the trigenic network. Despite their functional enrichment, trigenic interactions tended to link genes in distant bioprocesses and displayed a weaker magnitude than digenic interactions. We estimate that the global trigenic interaction network is ~100 times as large as the global digenic network, highlighting the potential for complex genetic interactions to affect the biology of inheritance, including the genotype-to-phenotype relationship.

A framework for exhaustively mapping functional missense variants

Although we now routinely sequence human genomes, we can confidently identify only a fraction of the sequence variants that have a functional impact. Here, we developed a deep mutational scanning framework that produces exhaustive maps for human missense variants by combining random codon mutagenesis and multiplexed functional variation assays with computational imputation and refinement. We applied this framework to four proteins corresponding to six human genes: UBE2I (encoding SUMO E2 conjugase), SUMO1 (small ubiquitin‐like modifier), TPK1 (thiamin pyrophosphokinase), and CALM1/2/3 (three genes encoding the protein calmodulin). The resulting maps recapitulate known protein features and confidently identify pathogenic variation. Assays potentially amenable to deep mutational scanning are already available for 57% of human disease genes, suggesting that DMS could ultimately map functional variation for all human disease genes.

Leveraging DNA damage response signaling to identify yeast genes controlling genome stability.

Oncogenesis frequently is accompanied by rampant genome instability, which fuels genetic heterogeneity and resistance to targeted cancer therapy. We have developed an approach that allows precise, quantitative measurement of genome instability in high-throughput format in the Saccharomyces cerevisiae model system. Our approach takes advantage of the strongly DNA damage-inducible gene RNR3, in conjunction with the reporter synthetic genetic array methodology, to infer mutants exhibiting genome instability by assaying for increased Rnr3 abundance. We screen for genome instability across a set of ~1000 essential and ~4200 nonessential mutant yeast alleles in untreated conditions and in the presence of the DNA-damaging agent methylmethane sulfonate. Our results provide broad insights into the cellular processes and pathways required for genome maintenance. Through comparison with existing genome instability screens, we isolated 130 genes that had not previously been linked to genome maintenance, 51% of which have human homologs. Several of these homologs are associated with a genome instability phenotype in human cells or are causally mutated in cancer. A comprehensive understanding of the processes required to prevent genome instability will facilitate a better understanding of its sources in oncogenesis.

Rapid and Efficient Plasmid Construction by Homologous Recombination in Yeast

The cloning of DNA fragments is a fundamental aspect of molecular biology. Traditional DNA cloning techniques rely on the ligation of an insert and a linearized plasmid that have been digested with restriction enzymes and the subsequent introduction of the ligated DNA into Escherichia coli for propagation. However, this method is limited by the availability of restriction sites, which often becomes problematic when cloning multiple or large DNA fragments. Furthermore, using traditional methods to clone multiple DNA fragments requires experience and multiple laborious steps. In this protocol, we describe a simple and efficient cloning method that relies on homologous recombination in the yeast Saccharomyces cerevisiae to assemble multiple DNA fragments, with 30-bp homology regions between the fragments, into one sophisticated construct. This method can easily be extended to clone plasmids for other organisms, such as bacteria, plants, and mammalian cells.

Expanding the Atlas of Functional Missense Variation for Human Genes

Although we now routinely sequence human genomes, we can confidently identify only a fraction of the sequence variants that have a functional impact. Here we developed a deep mutational scanning framework that produces exhaustive maps for human missense variants by combining random codon-mutagenesis and multiplexed functional variation assays with computational imputation and refinement. We applied this framework to four proteins corresponding to six human genes: UBE2I (encoding SUMO E2 conjugase), SUMO1 (small ubiquitin-like modifier), TPK1 (thiamin pyrophosphokinase), and CALM1/2/3 (three genes encoding the protein calmodulin). The resulting maps recapitulate known protein features, and confidently identify pathogenic variation. Assays potentially amenable to deep mutational scanning are already available for 57% of human disease genes, suggesting that DMS could ultimately map functional variation for all human disease genes.

Identifying pathogenicity of human variants via paralog-based yeast complementation

To better understand the health implications of personal genomes, we now face a largely unmet challenge to identify functional variants within disease-associated genes. Functional variants can be identified by trans-species complementation, e.g., by failure to rescue a yeast strain bearing a mutation in an orthologous human gene. Although orthologous complementation assays are powerful predictors of pathogenic variation, they are available for only a few percent of human disease genes. Here we systematically examine the question of whether complementation assays based on paralogy relationships can expand the number of human disease genes with functional variant detection assays. We tested over 1,000 paralogous human-yeast gene pairs for complementation, yielding 34 complementation relationships, of which 33 (97%) were novel. We found that paralog-based assays identified disease variants with success on par with that of orthology-based assays. Combining all homology-based assay results, we found that complementation can often identify pathogenic variants outside the homologous sequence region, presumably because of global effects on protein folding or stability. Within our search space, paralogy-based complementation more than doubled the number of human disease genes with a yeast-based complementation assay for disease variation.

Construction of Multifragment Plasmids by Homologous Recombination in Yeast

Over the past decade, the focus of cloning has shifted from constructing plasmids that express a single gene of interest to creating multigenic constructs that contain entire pathways or even whole genomes. Traditional cloning methods that rely on restriction digestion and ligation are limited by the number and size of fragments that can efficiently be combined. Here, we focus on the use of homologous-recombination-based DNA manipulation in the yeast Saccharomyces cerevisiae for the construction of plasmids from multiple DNA fragments. Owing to its simplicity and high efficiency, cloning by homologous recombination in yeast is very accessible and can be applied to high-throughput construction procedures. Its applications extend beyond yeast-centered purposes and include the cloning of large mammalian DNA sequences and entire bacterial genomes.