Exploring the yeast acetylome using functional genomics.

Lysine acetylation is a dynamic posttranslational modification with a well-defined role in regulating histones. The impact of acetylation on other cellular functions remains relatively uncharacterized. We explored the budding yeast acetylome with a functional genomics approach, assessing the effects of gene overexpression in the absence of lysine deacetylases (KDACs). We generated a network of 463 synthetic dosage lethal (SDL) interactions involving class I and II KDACs, revealing many cellular pathways regulated by different KDACs. A biochemical survey of genes interacting with the KDAC RPD3 identified 72 proteins acetylated in vivo. In-depth analysis of one of these proteins, Swi4, revealed a role for acetylation in G1-specific gene expression. Acetylation of Swi4 regulates interaction with its partner Swi6, both components of the SBF transcription factor. This study expands our view of the yeast acetylome, demonstrates the utility of functional genomic screens for exploring enzymatic pathways, and provides functional information that can be mined for future studies.

VID22 counteracts G-quadruplex-induced genome instability.

Genome instability is a condition characterized by the accumulation of genetic alterations and is a hallmark of cancer cells. To uncover new genes and cellular pathways affecting endogenous DNA damage and genome integrity, we exploited a Synthetic Genetic Array (SGA)-based screen in yeast. Among the positive genes, we identified VID22, reported to be involved in DNA double-strand break repair. vid22Δ cells exhibit increased levels of endogenous DNA damage, chronic DNA damage response activation and accumulate DNA aberrations in sequences displaying high probabilities of forming G-quadruplexes (G4-DNA). If not resolved, these DNA secondary structures can block the progression of both DNA and RNA polymerases and correlate with chromosome fragile sites. Vid22 binds to and protects DNA at G4-containing regions both in vitro and in vivo. Loss of VID22 causes an increase in gross chromosomal rearrangement (GCR) events dependent on G-quadruplex forming sequences. Moreover, the absence of Vid22 causes defects in the correct maintenance of G4-DNA rich elements, such as telomeres and mtDNA, and hypersensitivity to the G4-stabilizing ligand TMPyP4. We thus propose that Vid22 is directly involved in genome integrity maintenance as a novel regulator of G4 metabolism.

Skp, Cullin, F-box (SCF)-Met30 and SCF-Cdc4-Mediated Proteolysis of CENP-A Prevents Mislocalization of CENP-A for Chromosomal Stability in Budding Yeast.

Restricting the localization of the histone H3 variant CENP-A (Cse4 in yeast, CID in flies) to centromeres is essential for faithful chromosome segregation. Mislocalization of CENP-A leads to chromosomal instability (CIN) in yeast, fly and human cells. Overexpression and mislocalization of CENP-A has been observed in many cancers and this correlates with increased invasiveness and poor prognosis. Yet genes that regulate CENP-A levels and localization under physiological conditions have not been defined. In this study we used a genome-wide genetic screen to identify essential genes required for Cse4 homeostasis to prevent its mislocalization for chromosomal stability. We show that two Skp, Cullin, F-box (SCF) ubiquitin ligases with the evolutionarily conserved F-box proteins Met30 and Cdc4 interact and cooperatively regulate proteolysis of endogenous Cse4 and prevent its mislocalization for faithful chromosome segregation under physiological conditions. The interaction of Met30 with Cdc4 is independent of the D domain, which is essential for their homodimerization and ubiquitination of other substrates. The requirement for both Cdc4 and Met30 for ubiquitination is specifc for Cse4; and a common substrate for Cdc4 and Met30 has not previously been described. Met30 is necessary for the interaction between Cdc4 and Cse4, and defects in this interaction lead to stabilization and mislocalization of Cse4, which in turn contributes to CIN. We provide the first direct link between Cse4 mislocalization to defects in kinetochore structure and show that SCF-mediated proteolysis of Cse4 is a major mechanism that prevents stable maintenance of Cse4 at non-centromeric regions, thus ensuring faithful chromosome segregation. In summary, we have identified essential pathways that regulate cellular levels of endogenous Cse4 and shown that proteolysis of Cse4 by SCF-Met30/Cdc4 prevents mislocalization and CIN in unperturbed cells.

A genome-scale yeast library with inducible expression of individual genes.

The ability to switch a gene from off to on and monitor dynamic changes provides a powerful approach for probing gene function and elucidating causal regulatory relationships. Here, we developed and characterized YETI (Yeast Estradiol strains with Titratable Induction), a collection in which > 5,600 yeast genes are engineered for transcriptional inducibility with single-gene precision at their native loci and without plasmids. Each strain contains SGA screening markers and a unique barcode, enabling high-throughput genetics. We characterized YETI using growth phenotyping and BAR-seq screens, and we used a YETI allele to identify the regulon of Rof1, showing that it acts to repress transcription. We observed that strains with inducible essential genes that have low native expression can often grow without inducer. Analysis of data from eukaryotic and prokaryotic systems shows that native expression is a variable that can bias promoter-perturbing screens, including CRISPRi. We engineered a second expression system, Z3 EB42, that gives lower expression than Z3 EV, a feature enabling conditional activation and repression of lowly expressed essential genes that grow without inducer in the YETI library.

Genetic interaction profiles of regulatory kinases differ between environmental conditions and cellular states.

Cell growth and quiescence in eukaryotic cells is controlled by an evolutionarily conserved network of signaling pathways. Signal transduction networks operate to modulate a wide range of cellular processes and physiological properties when cells exit proliferative growth and initiate a quiescent state. How signaling networks function to respond to diverse signals that result in cell cycle exit and establishment of a quiescent state is poorly understood. Here, we studied the function of signaling pathways in quiescent cells using global genetic interaction mapping in the model eukaryotic cell, Saccharomyces cerevisiae (budding yeast). We performed pooled analysis of genotypes using molecular barcode sequencing (Bar-seq) to test the role of ~4,000 gene deletion mutants and ~12,000 pairwise interactions between all non-essential genes and the protein kinase genes TOR1, RIM15, and PHO85 in three different nutrient-restricted conditions in both proliferative and quiescent cells. We detect up to 10-fold more genetic interactions in quiescent cells than proliferative cells. We find that both individual gene effects and genetic interaction profiles vary depending on the specific pro-quiescence signal. The master regulator of quiescence, RIM15, shows distinct genetic interaction profiles in response to different starvation signals. However, vacuole-related functions show consistent genetic interactions with RIM15 in response to different starvation signals, suggesting that RIM15 integrates diverse signals to maintain protein homeostasis in quiescent cells. Our study expands genome-wide genetic interaction profiling to additional conditions, and phenotypes, and highlights the conditional dependence of epistasis.

Functional annotation of chemical libraries across diverse biological processes.

Chemical-genetic approaches offer the potential for unbiased functional annotation of chemical libraries. Mutations can alter the response of cells in the presence of a compound, revealing chemical-genetic interactions that can elucidate a compound's mode of action. We developed a highly parallel, unbiased yeast chemical-genetic screening system involving three key components. First, in a drug-sensitive genetic background, we constructed an optimized diagnostic mutant collection that is predictive for all major yeast biological processes. Second, we implemented a multiplexed (768-plex) barcode-sequencing protocol, enabling the assembly of thousands of chemical-genetic profiles. Finally, based on comparison of the chemical-genetic profiles with a compendium of genome-wide genetic interaction profiles, we predicted compound functionality. Applying this high-throughput approach, we screened seven different compound libraries and annotated their functional diversity. We further validated biological process predictions, prioritized a diverse set of compounds, and identified compounds that appear to have dual modes of action.

Errata: Functional annotation of chemical libraries across diverse biological processes.

This corrects the article DOI: 10.1038/nchembio.2436.

Errata: Functional annotation of chemical libraries across diverse biological processes.

This corrects the article DOI: 10.1038/nchembio.2436.

Systematic mapping of genetic interaction networks.

Genetic interactions influencing a phenotype of interest can be identified systematically using libraries of genetic tools that perturb biological systems in a defined manner. Systematic screens conducted in the yeast Saccharomyces cerevisiae have identified thousands of genetic interactions and provided insight into the global structure of biological networks. Techniques enabling systematic genetic interaction mapping have been extended to other single-celled organisms, the bacteria Escherichia coli and the yeast Schizosaccharomyces pombe, opening the way to comparative investigations of interaction networks. Genetic interaction screens in Caenorhabditis elegans, Drosophila melanogaster, and mammalian models are helping to improve our understanding of metazoan-specific signaling pathways. Together, our emerging knowledge of the genetic wiring diagrams of eukaryotic and prokaryotic cells is providing a new understanding of the relationship between genotype and phenotype.

Genetic interaction networks: toward an understanding of heritability.

Understanding the relationship between the genotypes and phenotypes of individuals is key for identifying genetic variants responsible for disease and developing successful therapeutic strategies. Mapping the phenotypic effects of individual genetic variants and their combinations in human populations presents numerous practical and statistical challenges. However, model organisms, such as the budding yeast Saccharomyces cerevisiae, provide an incredible set of molecular tools and advanced technologies that should be able to efficiently perform this task. In particular, large-scale genetic interaction screens in yeast and other model systems have revealed common properties of genetic interaction networks, many of which appear to be maintained over extensive evolutionary distances. Indeed, despite relatively low conservation of individual genes and their pairwise interactions, the overall topology of genetic interaction networks and the connections between broad biological processes may be similar in most organisms. Taking advantage of these general principles should provide a fundamental basis for mapping and predicting genetic interaction networks in humans.

Functional wiring of the yeast kinome revealed by global analysis of genetic network motifs.

A combinatorial genetic perturbation strategy was applied to interrogate the yeast kinome on a genome-wide scale. We assessed the global effects of gene overexpression or gene deletion to map an integrated genetic interaction network of synthetic dosage lethal (SDL) and loss-of-function genetic interactions (GIs) for 92 kinases, producing a meta-network of 8700 GIs enriched for pathways known to be regulated by cognate kinases. Kinases most sensitive to dosage perturbations had constitutive cell cycle or cell polarity functions under standard growth conditions. Condition-specific screens confirmed that the spectrum of kinase dosage interactions can be expanded substantially in activating conditions. An integrated network composed of systematic SDL, negative and positive loss-of-function GIs, and literature-curated kinase-substrate interactions revealed kinase-dependent regulatory motifs predictive of novel gene-specific phenotypes. Our study provides a valuable resource to unravel novel functional relationships and pathways regulated by kinases and outlines a general strategy for deciphering mutant phenotypes from large-scale GI networks.

SGAtools: one-stop analysis and visualization of array-based genetic interaction screens.

Screening genome-wide sets of mutants for fitness defects provides a simple but powerful approach for exploring gene function, mapping genetic networks and probing mechanisms of drug action. For yeast and other microorganisms with global mutant collections, genetic or chemical-genetic interactions can be effectively quantified by growing an ordered array of strains on agar plates as individual colonies, and then scoring the colony size changes in response to a genetic or environmental perturbation. To do so, requires efficient tools for the extraction and analysis of quantitative data. Here, we describe SGAtools (http://sgatools.ccbr.utoronto.ca), a web-based analysis system for designer genetic screens. SGAtools outlines a series of guided steps that allow the user to quantify colony sizes from images of agar plates, correct for systematic biases in the observations and calculate a fitness score relative to a control experiment. The data can also be visualized online to explore the colony sizes on individual plates, view the distribution of resulting scores, highlight genes with the strongest signal and perform Gene Ontology enrichment analysis.

A negative genetic interaction map in isogenic cancer cell lines reveals cancer cell vulnerabilities.

Improved efforts are necessary to define the functional product of cancer mutations currently being revealed through large-scale sequencing efforts. Using genome-scale pooled shRNA screening technology, we mapped negative genetic interactions across a set of isogenic cancer cell lines and confirmed hundreds of these interactions in orthogonal co-culture competition assays to generate a high-confidence genetic interaction network of differentially essential or differential essentiality (DiE) genes. The network uncovered examples of conserved genetic interactions, densely connected functional modules derived from comparative genomics with model systems data, functions for uncharacterized genes in the human genome and targetable vulnerabilities. Finally, we demonstrate a general applicability of DiE gene signatures in determining genetic dependencies of other non-isogenic cancer cell lines. For example, the PTEN(-/-) DiE genes reveal a signature that can preferentially classify PTEN-dependent genotypes across a series of non-isogenic cell lines derived from the breast, pancreas and ovarian cancers. Our reference network suggests that many cancer vulnerabilities remain to be discovered through systematic derivation of a network of differentially essential genes in an isogenic cancer cell model.

A systems biology approach reveals the role of a novel methyltransferase in response to chemical stress and lipid homeostasis.

Using small molecule probes to understand gene function is an attractive approach that allows functional characterization of genes that are dispensable in standard laboratory conditions and provides insight into the mode of action of these compounds. Using chemogenomic assays we previously identified yeast Crg1, an uncharacterized SAM-dependent methyltransferase, as a novel interactor of the protein phosphatase inhibitor cantharidin. In this study we used a combinatorial approach that exploits contemporary high-throughput techniques available in Saccharomyces cerevisiae combined with rigorous biological follow-up to characterize the interaction of Crg1 with cantharidin. Biochemical analysis of this enzyme followed by a systematic analysis of the interactome and lipidome of CRG1 mutants revealed that Crg1, a stress-responsive SAM-dependent methyltransferase, methylates cantharidin in vitro. Chemogenomic assays uncovered that lipid-related processes are essential for cantharidin resistance in cells sensitized by deletion of the CRG1 gene. Lipidome-wide analysis of mutants further showed that cantharidin induces alterations in glycerophospholipid and sphingolipid abundance in a Crg1-dependent manner. We propose that Crg1 is a small molecule methyltransferase important for maintaining lipid homeostasis in response to drug perturbation. This approach demonstrates the value of combining chemical genomics with other systems-based methods for characterizing proteins and elucidating previously unknown mechanisms of action of small molecule inhibitors.

Quantitative analysis of fitness and genetic interactions in yeast on a genome scale.

Global quantitative analysis of genetic interactions is a powerful approach for deciphering the roles of genes and mapping functional relationships among pathways. Using colony size as a proxy for fitness, we developed a method for measuring fitness-based genetic interactions from high-density arrays of yeast double mutants generated by synthetic genetic array (SGA) analysis. We identified several experimental sources of systematic variation and developed normalization strategies to obtain accurate single- and double-mutant fitness measurements, which rival the accuracy of other high-resolution studies. We applied the SGA score to examine the relationship between physical and genetic interaction networks, and we found that positive genetic interactions connect across functionally distinct protein complexes revealing a network of genetic suppression among loss-of-function alleles.

Global analysis of the yeast knockout phenome.

Genome-wide phenotypic screens in the budding yeast Saccharomyces cerevisiae, enabled by its knockout collection, have produced the largest, richest, and most systematic phenotypic description of any organism. However, integrative analyses of this rich data source have been virtually impossible because of the lack of a central data repository and consistent metadata annotations. Here, we describe the aggregation, harmonization, and analysis of ~14,500 yeast knockout screens, which we call Yeast Phenome. Using this unique dataset, we characterized two unknown genes (YHR045W and YGL117W) and showed that tryptophan starvation is a by-product of many chemical treatments. Furthermore, we uncovered an exponential relationship between phenotypic similarity and intergenic distance, which suggests that gene positions in both yeast and human genomes are optimized for function.

Systematic exploration of synergistic drug pairs.

Drug synergy allows a therapeutic effect to be achieved with lower doses of component drugs. Drug synergy can result when drugs target the products of genes that act in parallel pathways ('specific synergy'). Such cases of drug synergy should tend to correspond to synergistic genetic interaction between the corresponding target genes. Alternatively, 'promiscuous synergy' can arise when one drug non-specifically increases the effects of many other drugs, for example, by increased bioavailability. To assess the relative abundance of these drug synergy types, we examined 200 pairs of antifungal drugs in S. cerevisiae. We found 38 antifungal synergies, 37 of which were novel. While 14 cases of drug synergy corresponded to genetic interaction, 92% of the synergies we discovered involved only six frequently synergistic drugs. Although promiscuity of four drugs can be explained under the bioavailability model, the promiscuity of Tacrolimus and Pentamidine was completely unexpected. While many drug synergies correspond to genetic interactions, the majority of drug synergies appear to result from non-specific promiscuous synergy.

Protein complexes are central in the yeast genetic landscape.

If perturbing two genes together has a stronger or weaker effect than expected, they are said to genetically interact. Genetic interactions are important because they help map gene function, and functionally related genes have similar genetic interaction patterns. Mapping quantitative (positive and negative) genetic interactions on a global scale has recently become possible. This data clearly shows groups of genes connected by predominantly positive or negative interactions, termed monochromatic groups. These groups often correspond to functional modules, like biological processes or complexes, or connections between modules. However it is not yet known how these patterns globally relate to known functional modules. Here we systematically study the monochromatic nature of known biological processes using the largest quantitative genetic interaction data set available, which includes fitness measurements for ∼5.4 million gene pairs in the yeast Saccharomyces cerevisiae. We find that only 10% of biological processes, as defined by Gene Ontology annotations, and less than 1% of inter-process connections are monochromatic. Further, we show that protein complexes are responsible for a surprisingly large fraction of these patterns. This suggests that complexes play a central role in shaping the monochromatic landscape of biological processes. Altogether this work shows that both positive and negative monochromatic patterns are found in known biological processes and in their connections and that protein complexes play an important role in these patterns. The monochromatic processes, complexes and connections we find chart a hierarchical and modular map of sensitive and redundant biological systems in the yeast cell that will be useful for gene function prediction and comparison across phenotypes and organisms. Furthermore the analysis methods we develop are applicable to other species for which genetic interactions will progressively become more available.

DRYGIN: a database of quantitative genetic interaction networks in yeast.

Genetic interactions are highly informative for deciphering the underlying functional principles that govern how genes control cell processes. Recent developments in Synthetic Genetic Array (SGA) analysis enable the mapping of quantitative genetic interactions on a genome-wide scale. To facilitate access to this resource, which will ultimately represent a complete genetic interaction network for a eukaryotic cell, we developed DRYGIN (Data Repository of Yeast Genetic Interactions)-a web database system that aims at providing a central platform for yeast genetic network analysis and visualization. In addition to providing an interface for searching the SGA genetic interactions, DRYGIN also integrates other data sources, in order to associate the genetic interactions with pathway information, protein complexes, other binary genetic and physical interactions, and Gene Ontology functional annotation. DRYGIN version 1.0 currently holds more than 5.4 million measurements of genetic interacting pairs involving approximately 4500 genes, and is available at http://drygin.ccbr.utoronto.ca.

Genetic interactions reveal the evolutionary trajectories of duplicate genes.

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.