Competitive genomic screens of barcoded yeast libraries

By virtue of advances in next generation sequencing technologies, we have access to new genome sequences almost daily. The tempo of these advances is accelerating, promising greater depth and breadth. In light of these extraordinary advances, the need for fast, parallel methods to define gene function becomes ever more important. Collections of genome-wide deletion mutants in yeasts and E. coli have served as workhorses for functional characterization of gene function, but this approach is not scalable, current gene-deletion approaches require each of the thousands of genes that comprise a genome to be deleted and verified. Only after this work is complete can we pursue high-throughput phenotyping. Over the past decade, our laboratory has refined a portfolio of competitive, miniaturized, high-throughput genome-wide assays that can be performed in parallel. This parallelization is possible because of the inclusion of DNA 'tags', or 'barcodes,' into each mutant, with the barcode serving as a proxy for the mutation and one can measure the barcode abundance to assess mutant fitness. In this study, we seek to fill the gap between DNA sequence and barcoded mutant collections. To accomplish this we introduce a combined transposon disruption-barcoding approach that opens up parallel barcode assays to newly sequenced, but poorly characterized microbes. To illustrate this approach we present a new Candida albicans barcoded disruption collection and describe how both microarray-based and next generation sequencing-based platforms can be used to collect 10,000-1,000,000 gene-gene and drug-gene interactions in a single experiment.

Design, synthesis, and characterization of a highly effective Hog1 inhibitor: a powerful tool for analyzing MAP kinase signaling in yeast

The Saccharomyces cerevisiae High-Osmolarity Glycerol (HOG) pathway is a conserved mitogen-activated protein kinase (MAPK) signal transduction system that often serves as a model to analyze systems level properties of MAPK signaling. Hog1, the MAPK of the HOG-pathway, can be activated by various environmental cues and it controls transcription, translation, transport, and cell cycle adaptations in response to stress conditions. A powerful means to study signaling in living cells is to use kinase inhibitors; however, no inhibitor targeting wild-type Hog1 exists to date. Herein, we describe the design, synthesis, and biological application of small molecule inhibitors that are cell-permeable, fast-acting, and highly efficient against wild-type Hog1. These compounds are potent inhibitors of Hog1 kinase activity both in vitro and in vivo. Next, we use these novel inhibitors to pinpoint the time of Hog1 action during recovery from G(1) checkpoint arrest, providing further evidence for a specific role of Hog1 in regulating cell cycle resumption during arsenite stress. Hence, we describe a novel tool for chemical genetic analysis of MAPK signaling and provide novel insights into Hog1 action.

Dosage suppression genetic interaction networks enhance functional wiring diagrams of the cell

Dosage suppression is a genetic interaction in which overproduction of one gene rescues a mutant phenotype of another gene. Although dosage suppression is known to map functional connections among genes, the extent to which it might illuminate global cellular functions is unclear. Here we analyze a network of interactions linking dosage suppressors to 437 essential genes in yeast. For 424 genes, we curated interactions from the literature. Analyses revealed that many dosage suppression interactions occur between functionally related genes and that the majority do not overlap with other types of genetic or physical interactions. To confirm the generality of these network properties, we experimentally identified dosage suppressors for 29 genes from pooled populations of temperature-sensitive mutant cells transformed with a high-copy molecular-barcoded open reading frame library, MoBY-ORF 2.0. We classified 87% of the 1,640 total interactions into four general types of suppression mechanisms, which provided insight into their relative frequencies. This work suggests that integrating the results of dosage suppression studies with other interaction networks could generate insights into the functional wiring diagram of a cell.

A comprehensive platform for highly multiplexed mammalian functional genetic screens

Background

Genome-wide screening in human and mouse cells using RNA interference and open reading frame over-expression libraries is rapidly becoming a viable experimental approach for many research labs. There are a variety of gene expression modulation libraries commercially available, however, detailed and validated protocols as well as the reagents necessary for deconvolving genome-scale gene screens using these libraries are lacking. As a solution, we designed a comprehensive platform for highly multiplexed functional genetic screens in human, mouse and yeast cells using popular, commercially available gene modulation libraries. The Gene Modulation Array Platform (GMAP) is a single microarray-based detection solution for deconvolution of loss and gain-of-function pooled screens.

Results

Experiments with specially constructed lentiviral-based plasmid pools containing ~78,000 shRNAs demonstrated that the GMAP is capable of deconvolving genome-wide shRNA "dropout" screens. Further experiments with a larger, ~90,000 shRNA pool demonstrate that equivalent results are obtained from plasmid pools and from genomic DNA derived from lentivirus infected cells. Parallel testing of large shRNA pools using GMAP and next-generation sequencing methods revealed that the two methods provide valid and complementary approaches to deconvolution of genome-wide shRNA screens. Additional experiments demonstrated that GMAP is equivalent to similar microarray-based products when used for deconvolution of open reading frame over-expression screens.

Conclusion

Herein, we demonstrate four major applications for the GMAP resource, including deconvolution of pooled RNAi screens in cells with at least 90,000 distinct shRNAs. We also provide detailed methodologies for pooled shRNA screen readout using GMAP and compare next-generation sequencing to GMAP (i.e. microarray) based deconvolution methods.

The synthetic genetic interaction network reveals small molecules that target specific pathways in Sacchromyces cerevisiae

High-throughput elucidation of synthetic genetic interactions (SGIs) has contributed to a systems-level understanding of genetic robustness and fault-tolerance encoded in the genome. Pathway targets of various compounds have been predicted by comparing chemical-genetic synthetic interactions to a network of SGIs. We demonstrate that the SGI network can also be used in a powerful reverse pathway-to-drug approach for identifying compounds that target specific pathways of interest. Using the SGI network, the method identifies an indicator gene that may serve as a good candidate for screening a library of compounds. The indicator gene is selected so that compounds found to produce sensitivity in mutants deleted for the indicator gene are likely to abrogate the target pathway. We tested the utility of the SGI network for pathway-to-drug discovery using the DNA damage checkpoint as the target pathway. An analysis of the compendium of synthetic lethal interactions in yeast showed that superoxide dismutase 1 (SOD1) has significant SGI connectivity with a large subset of DNA damage checkpoint and repair (DDCR) genes in Saccharomyces cerevisiae, and minimal SGIs with non-DDCR genes. We screened a sod1Δ strain against three National Cancer Institute (NCI) compound libraries using a soft agar high-throughput halo assay. Fifteen compounds out of ∼3100 screened showed selective toxicity toward sod1Δ relative to the isogenic wild type (wt) strain. One of these, 1A08, caused a transient increase in growth in the presence of sublethal doses of DNA damaging agents, suggesting that 1A08 inhibits DDCR signaling in yeast. Genome-wide screening of 1A08 against the library of viable homozygous deletion mutants further supported DDCR as the relevant targeted pathway of 1A08. When assayed in human HCT-116 colorectal cancer cells, 1A08 caused DNA-damage resistant DNA synthesis and blocked the DNA-damage checkpoint selectively in S-phase.

Systematic exploration of essential yeast gene function with temperature-sensitive mutants

Conditional temperature-sensitive (ts) mutations are valuable reagents for studying essential genes in the yeast Saccharomyces cerevisiae. We constructed 787 ts strains, covering 497 (∼45%) of the 1,101 essential yeast genes, with ∼30% of the genes represented by multiple alleles. All of the alleles are integrated into their native genomic locus in the S288C common reference strain and are linked to a kanMX selectable marker, allowing further genetic manipulation by synthetic genetic array (SGA)-based, high-throughput methods. We show two such manipulations: barcoding of 440 strains, which enables chemical-genetic suppression analysis, and the construction of arrays of strains carrying different fluorescent markers of subcellular structure, which enables quantitative analysis of phenotypes using high-content screening. Quantitative analysis of a GFP-tubulin marker identified roles for cohesin and condensin genes in spindle disassembly. This mutant collection should facilitate a wide range of systematic studies aimed at understanding the functions of essential genes.

Displaying chemical information on a biological network using Cytoscape

Cytoscape is an open-source software package that is widely used to integrate and visualize diverse data sets in biology. This chapter explains how to use Cytoscape to integrate open-source chemical information with a biological network. By visualizing information about known compound-target interactions in the context of a biological network of interest, one can rapidly identify novel avenues to perturb the system with compounds and, for example, potentially identify therapeutically relevant targets. Herein, two different protocols are explained in detail, with no prior knowledge of Cytoscape assumed, which demonstrate how to incorporate data from the ChEMBL database with either a gene-gene or a protein-protein interaction network. ChEMBL is a very large, open-source repository of compound-target information available from the European Molecular Biology Laboratory.

The automated cell: compound and environment screening system (ACCESS) for chemogenomic screening

The automated cell, compound and environment screening system (ACCESS) was developed as an automated platform for chemogenomic research. In the yeast Saccharomyces cerevisiae, a number of genomic screens rely on the modulation of gene dose to determine the mode of action of bioactive compounds or the effects of environmental/compound perturbations. These and other phenotypic experiments have been shown to benefit from high-resolution growth curves and a highly automated controlled environment system that enables a wide range of multi-well assays that can be run over many days without any manual intervention. Furthermore, precise control of drug dosing, timing of drug exposure, and precise timing of cell harvesting at specific generation times are important for optimal results. Some of these benefits include the ability to derive fine distinctions between growth rates of mutant strains (1) and the discovery of novel compounds and drug targets (2). The automation has also enabled large-scale screening projects with over 100,000 unique compounds screened to date including a thousand genome-wide screens (3). The ACCESS system also has a diverse set of software tools to enable users to set up, run, annotate, and evaluate complex screens with minimal training.

Onge RP, Deutschbauer A, Nislow C. Gene annotation and drug target discovery in Candida albicans with a tagged transposon mutant collection

Candida albicans is the most common human fungal pathogen, causing infections that can be lethal in immunocompromised patients. Although Saccharomyces cerevisiae has been used as a model for C. albicans, it lacks C. albicans' diverse morphogenic forms and is primarily non-pathogenic. Comprehensive genetic analyses that have been instrumental for determining gene function in S. cerevisiae are hampered in C. albicans, due in part to limited resources to systematically assay phenotypes of loss-of-function alleles. Here, we constructed and screened a library of 3633 tagged heterozygous transposon disruption mutants, using them in a competitive growth assay to examine nutrient- and drug-dependent haploinsufficiency. We identified 269 genes that were haploinsufficient in four growth conditions, the majority of which were condition-specific. These screens identified two new genes necessary for filamentous growth as well as ten genes that function in essential processes. We also screened 57 chemically diverse compounds that more potently inhibited growth of C. albicans versus S. cerevisiae. For four of these compounds, we examined the genetic basis of this differential inhibition. Notably, Sec7p was identified as the target of brefeldin A in C. albicans screens, while S. cerevisiae screens with this compound failed to identify this target. We also uncovered a new C. albicans-specific target, Tfp1p, for the synthetic compound 0136-0228. These results highlight the value of haploinsufficiency screens directly in this pathogen for gene annotation and drug target identification.

Candida albicans is a normal inhabitant in our bodies, yet it can become pathogenic and cause infections that range from the superficial in healthy individuals to deadly in the immunocompromised. Comprehensive genetic analysis of C. albicans to identify mechanisms of virulence and new treatment strategies has been hampered by limited, publically accessible genomic resources. By combining the principles of Saccharomyces cerevisiae strain tagging with transposon mutagenesis to generate individually tagged mutants, we created the first entirely public resource that allows simultaneous measurement of strain fitness of ∼60% of the genome in a wide range of experimental treatments. By identifying genes that confer a fitness or growth defect when reduced in copy number, we uncovered genes whose protein products represent potential antifungal targets. Moreover, screening this strain collection with chemical compounds allowed us to identify anticandidal chemicals while concurrently gaining insight into their cellular mechanism of action. This resource, combined with straightforward screening methodology, provides powerful tools to generate hypotheses for functional annotation of the genome, and our results highlight the value of direct versus model-based pathogen studies.

Genome-wide screen in Saccharomyces cerevisiae identifies vacuolar protein sorting, autophagy, biosynthetic, and tRNA methylation genes involved in life span regulation

The study of the chronological life span of Saccharomyces cerevisiae, which measures the survival of populations of non-dividing yeast, has resulted in the identification of homologous genes and pathways that promote aging in organisms ranging from yeast to mammals. Using a competitive genome-wide approach, we performed a screen of a complete set of approximately 4,800 viable deletion mutants to identify genes that either increase or decrease chronological life span. Half of the putative short-/long-lived mutants retested from the primary screen were confirmed, demonstrating the utility of our approach. Deletion of genes involved in vacuolar protein sorting, autophagy, and mitochondrial function shortened life span, confirming that respiration and degradation processes are essential for long-term survival. Among the genes whose deletion significantly extended life span are ACB1, CKA2, and TRM9, implicated in fatty acid transport and biosynthesis, cell signaling, and tRNA methylation, respectively. Deletion of these genes conferred heat-shock resistance, supporting the link between life span extension and cellular protection observed in several model organisms. The high degree of conservation of these novel yeast longevity determinants in other species raises the possibility that their role in senescence might be conserved.

A survey of yeast genomic assays for drug and target discovery

Over the past decade, the development and application of chemical genomic assays using the model organism Saccharomyces cerevisiae has provided powerful methods to identify the mechanism of action of known drugs and novel small molecules in vivo. These assays identify drug target candidates, genes involved in buffering drug target pathways and also help to define the general cellular response to small molecules. In this review, we examine current yeast chemical genomic assays and summarize the potential applications of each approach.

A universal TagModule collection for parallel genetic analysis of microorganisms

Systems-level analyses of non-model microorganisms are limited by the existence of numerous uncharacterized genes and a corresponding over-reliance on automated computational annotations. One solution to this challenge is to disrupt gene function using DNA tag technology, which has been highly successful in parallelizing reverse genetics in Saccharomyces cerevisiae and has led to discoveries in gene function, genetic interactions and drug mechanism of action. To extend the yeast DNA tag methodology to a wide variety of microorganisms and applications, we have created a universal, sequence-verified TagModule collection. A hallmark of the 4280 TagModules is that they are cloned into a Gateway entry vector, thus facilitating rapid transfer to any compatible genetic system. Here, we describe the application of the TagModules to rapidly generate tagged mutants by transposon mutagenesis in the metal-reducing bacterium Shewanella oneidensis MR-1 and the pathogenic yeast Candida albicans. Our results demonstrate the optimal hybridization properties of the TagModule collection, the flexibility in applying the strategy to diverse microorganisms and the biological insights that can be gained from fitness profiling tagged mutant collections. The publicly available TagModule collection is a platform-independent resource for the functional genomics of a wide range of microbial systems in the post-genome era.

Highly-multiplexed barcode sequencing: an efficient method for parallel analysis of pooled samples

Next-generation sequencing has proven an extremely effective technology for molecular counting applications where the number of sequence reads provides a digital readout for RNA-seq, ChIP-seq, Tn-seq and other applications. The extremely large number of sequence reads that can be obtained per run permits the analysis of increasingly complex samples. For lower complexity samples, however, a point of diminishing returns is reached when the number of counts per sequence results in oversampling with no increase in data quality. A solution to making next-generation sequencing as efficient and affordable as possible involves assaying multiple samples in a single run. Here, we report the successful 96-plexing of complex pools of DNA barcoded yeast mutants and show that such 'Bar-seq' assessment of these samples is comparable with data provided by barcode microarrays, the current benchmark for this application. The cost reduction and increased throughput permitted by highly multiplexed sequencing will greatly expand the scope of chemogenomics assays and, equally importantly, the approach is suitable for other sequence counting applications that could benefit from massive parallelization.

Genotype to phenotype: a complex problem

We generated a high-resolution whole-genome sequence and individually deleted 5100 genes in Sigma1278b, a Saccharomyces cerevisiae strain closely related to reference strain S288c. Similar to the variation between human individuals, Sigma1278b and S288c average 3.2 single-nucleotide polymorphisms per kilobase. A genome-wide comparison of deletion mutant phenotypes identified a subset of genes that were conditionally essential by strain, including 44 essential genes unique to Sigma1278b and 13 unique to S288c. Genetic analysis indicates the conditional phenotype was most often governed by complex genetic interactions, depending on multiple background-specific modifiers. Our comprehensive analysis suggests that the presence of a complex set of modifiers will often underlie the phenotypic differences between individuals.

Systematic analysis of genome-wide fitness data in yeast reveals novel gene function and drug action

The relationship between co-fitness and co-inhibition of genes in chemicogenomic yeast screens provides insights into gene function and drug target prediction.

We systematically analyzed the relationships between gene fitness profiles (co-fitness) and drug inhibition profiles (co-inhibition) from several hundred chemogenomic screens in yeast. Co-fitness predicted gene functions distinct from those derived from other assays and identified conditionally dependent protein complexes. Co-inhibitory compounds were weakly correlated by structure and therapeutic class. We developed an algorithm predicting protein targets of chemical compounds and verified its accuracy with experimental testing. Fitness data provide a novel, systems-level perspective on the cell.

Onge RP, VanderSluis B, Makhnevych T, Vizeacoumar FJ, Alizadeh S, Bahr S, Brost RL, Chen Y, Cokol M, Deshpande R, Li Z, Lin ZY, Liang W, Marback M, Paw J, San Luis BJ, Shuteriqi E, Hin Yan Tong A, van Dyk N, Wallace IM, Whitney JA, Weirauch MT, Zhong G, Zhu H, Houry WA, Brudno M, Ragibizadeh S, Papp B, Pál C, Roth FP, Giaever G, Nislow C, Troyanskaya OG, Bussey H, Bader GD, Gingras AC, Morris QD, Kim PM, Kaiser CA, Myers CL, Andrews BJ, Boone C. The genetic landscape of a cell

A genome-scale genetic interaction map was constructed by examining 5.4 million gene-gene pairs for synthetic genetic interactions, generating quantitative genetic interaction profiles for approximately 75% of all genes in the budding yeast, Saccharomyces cerevisiae. A network based on genetic interaction profiles reveals a functional map of the cell in which genes of similar biological processes cluster together in coherent subsets, and highly correlated profiles delineate specific pathways to define gene function. The global network identifies functional cross-connections between all bioprocesses, mapping a cellular wiring diagram of pleiotropy. Genetic interaction degree correlated with a number of different gene attributes, which may be informative about genetic network hubs in other organisms. We also demonstrate that extensive and unbiased mapping of the genetic landscape provides a key for interpretation of chemical-genetic interactions and drug target identification.

Quantitative phenotyping via deep barcode sequencing

Next-generation DNA sequencing technologies have revolutionized diverse genomics applications, including de novo genome sequencing, SNP detection, chromatin immunoprecipitation, and transcriptome analysis. Here we apply deep sequencing to genome-scale fitness profiling to evaluate yeast strain collections in parallel. This method, Barcode analysis by Sequencing, or "Bar-seq," outperforms the current benchmark barcode microarray assay in terms of both dynamic range and throughput. When applied to a complex chemogenomic assay, Bar-seq quantitatively identifies drug targets, with performance superior to the benchmark microarray assay. We also show that Bar-seq is well-suited for a multiplex format. We completely re-sequenced and re-annotated the yeast deletion collection using deep sequencing, found that approximately 20% of the barcodes and common priming sequences varied from expectation, and used this revised list of barcode sequences to improve data quality. Together, this new assay and analysis routine provide a deep-sequencing-based toolkit for identifying gene-environment interactions on a genome-wide scale.

Novel insights into iron metabolism by integrating deletome and transcriptome analysis in an iron deficiency model of the yeast Saccharomyces cerevisiae

Background

Iron-deficiency anemia is the most prevalent form of anemia world-wide. The yeast Saccharomyces cerevisiae has been used as a model of cellular iron deficiency, in part because many of its cellular pathways are conserved. To better understand how cells respond to changes in iron availability, we profiled the yeast genome with a parallel analysis of homozygous deletion mutants to identify essential components and cellular processes required for optimal growth under iron-limited conditions. To complement this analysis, we compared those genes identified as important for fitness to those that were differentially-expressed in the same conditions. The resulting analysis provides a global perspective on the cellular processes involved in iron metabolism.

Results

Using functional profiling, we identified several genes known to be involved in high affinity iron uptake, in addition to novel genes that may play a role in iron metabolism. Our results provide support for the primary involvement in iron homeostasis of vacuolar and endosomal compartments, as well as vesicular transport to and from these compartments. We also observed an unexpected importance of the peroxisome for growth in iron-limited media. Although these components were essential for growth in low-iron conditions, most of them were not differentially-expressed. Genes with altered expression in iron deficiency were mainly associated with iron uptake and transport mechanisms, with little overlap with those that were functionally required. To better understand this relationship, we used expression-profiling of selected mutants that exhibited slow growth in iron-deficient conditions, and as a result, obtained additional insight into the roles of CTI6, DAP1, MRS4 and YHR045W in iron metabolism.

Conclusion

Comparison between functional and gene expression data in iron deficiency highlighted the complementary utility of these two approaches to identify important functional components. This should be taken into consideration when designing and analyzing data from these type of studies. We used this and other published data to develop a molecular interaction network of iron metabolism in yeast.

Computationally driven, quantitative experiments discover genes required for mitochondrial biogenesis

Mitochondria are central to many cellular processes including respiration, ion homeostasis, and apoptosis. Using computational predictions combined with traditional quantitative experiments, we have identified 100 proteins whose deficiency alters mitochondrial biogenesis and inheritance in Saccharomyces cerevisiae. In addition, we used computational predictions to perform targeted double-mutant analysis detecting another nine genes with synthetic defects in mitochondrial biogenesis. This represents an increase of about 25% over previously known participants. Nearly half of these newly characterized proteins are conserved in mammals, including several orthologs known to be involved in human disease. Mutations in many of these genes demonstrate statistically significant mitochondrial transmission phenotypes more subtle than could be detected by traditional genetic screens or high-throughput techniques, and 47 have not been previously localized to mitochondria. We further characterized a subset of these genes using growth profiling and dual immunofluorescence, which identified genes specifically required for aerobic respiration and an uncharacterized cytoplasmic protein required for normal mitochondrial motility. Our results demonstrate that by leveraging computational analysis to direct quantitative experimental assays, we have characterized mutants with subtle mitochondrial defects whose phenotypes were undetected by high-throughput methods.

Mitochondria are the proverbial powerhouses of the cell, running the fundamental biochemical processes that produce energy from nutrients using oxygen. These processes are conserved in all eukaryotes, from humans to model organisms such as baker's yeast. In humans, mitochondrial dysfunction plays a role in a variety of diseases, including diabetes, neuromuscular disorders, and aging. In order to better understand fundamental mitochondrial biology, we studied genes involved in mitochondrial biogenesis in the yeast S. cerevisiae, discovering over 100 proteins with novel roles in this process. These experiments assigned function to 5% of the genes whose function was not known. In order to achieve this rapid rate of discovery, we developed a system incorporating highly quantitative experimental assays and an integrated, iterative process of computational protein function prediction. Beginning from relatively little prior knowledge, we found that computational predictions achieved about 60% accuracy and rapidly guided our laboratory work towards hundreds of promising candidate genes. Thus, in addition to providing a more thorough understanding of mitochondrial biology, this study establishes a framework for successfully integrating computation and experimentation to drive biological discovery. A companion manuscript, published in PLoS Computational Biology (doi:10.1371/journal.pcbi.1000322), discusses observations and conclusions important for the computational community.

Combining chemical genomics screens in yeast to reveal spectrum of effects of chemical inhibition of sphingolipid biosynthesis

BACKGROUND

Single genome-wide screens for the effect of altered gene dosage on drug sensitivity in the model organism Saccharomyces cerevisiae provide only a partial picture of the mechanism of action of a drug.

RESULTS

Using the example of the tumor cell invasion inhibitor dihydromotuporamine C, we show that a more complete picture of drug action can be obtained by combining different chemical genomics approaches--analysis of the sensitivity of rho0 cells lacking mitochondrial DNA, drug-induced haploinsufficiency, suppression of drug sensitivity by gene overexpression and chemical-genetic synthetic lethality screening using strains deleted of nonessential genes. Killing of yeast by this chemical requires a functional mitochondrial electron-transport chain and cytochrome c heme lyase function. However, we find that it does not require genes associated with programmed cell death in yeast. The chemical also inhibits endocytosis and intracellular vesicle trafficking and interferes with vacuolar acidification in yeast and in human cancer cells. These effects can all be ascribed to inhibition of sphingolipid biosynthesis by dihydromotuporamine C.

CONCLUSION

Despite their similar conceptual basis, namely altering drug sensitivity by modifying gene dosage, each of the screening approaches provided a distinct set of information that, when integrated, revealed a more complete picture of the mechanism of action of a drug on cells.

Combination chemical genetics

Predicting the behavior of living organisms is an enormous challenge given their vast complexity. Efforts to model biological systems require large datasets generated by physical binding experiments and perturbation studies. Genetic perturbations have proven important and are greatly facilitated by the advent of comprehensive mutant libraries in model organisms. Small-molecule chemical perturbagens provide a complementary approach, especially for systems that lack mutant libraries, and can easily probe the function of essential genes. Though single chemical or genetic perturbations provide crucial information associating individual components (for example, genes, proteins or small molecules) with pathways or phenotypes, functional relationships between pathways and modules of components are most effectively obtained from combined perturbation experiments. Here we review the current state of and discuss some future directions for 'combination chemical genetics', the systematic application of multiple chemical or mixed chemical and genetic perturbations, both to gain insight into biological systems and to facilitate medical discoveries.

Off-target effects of psychoactive drugs revealed by genome-wide assays in yeast

To better understand off-target effects of widely prescribed psychoactive drugs, we performed a comprehensive series of chemogenomic screens using the budding yeast Saccharomyces cerevisiae as a model system. Because the known human targets of these drugs do not exist in yeast, we could employ the yeast gene deletion collections and parallel fitness profiling to explore potential off-target effects in a genome-wide manner. Among 214 tested, documented psychoactive drugs, we identified 81 compounds that inhibited wild-type yeast growth and were thus selected for genome-wide fitness profiling. Many of these drugs had a propensity to affect multiple cellular functions. The sensitivity profiles of half of the analyzed drugs were enriched for core cellular processes such as secretion, protein folding, RNA processing, and chromatin structure. Interestingly, fluoxetine (Prozac) interfered with establishment of cell polarity, cyproheptadine (Periactin) targeted essential genes with chromatin-remodeling roles, while paroxetine (Paxil) interfered with essential RNA metabolism genes, suggesting potential secondary drug targets. We also found that the more recently developed atypical antipsychotic clozapine (Clozaril) had no fewer off-target effects in yeast than the typical antipsychotics haloperidol (Haldol) and pimozide (Orap). Our results suggest that model organism pharmacogenetic studies provide a rational foundation for understanding the off-target effects of clinically important psychoactive agents and suggest a rational means both for devising compound derivatives with fewer side effects and for tailoring drug treatment to individual patient genotypes.

An integrated platform of genomic assays reveals small-molecule bioactivities

Bioactive compounds are widely used to modulate protein function and can serve as important leads for drug development. Identifying the in vivo targets of these compounds remains a challenge. Using yeast, we integrated three genome-wide gene-dosage assays to measure the effect of small molecules in vivo. A single TAG microarray was used to resolve the fitness of strains derived from pools of (i) homozygous deletion mutants, (ii) heterozygous deletion mutants and (iii) genomic library transformants. We demonstrated, with eight diverse reference compounds, that integration of these three chemogenomic profiles improves the sensitivity and specificity of small-molecule target identification. We further dissected the mechanism of action of two protein phosphatase inhibitors and in the process developed a framework for the rational design of multidrug combinations to sensitize cells with specific genotypes more effectively. Finally, we applied this platform to 188 novel synthetic chemical compounds and identified both potential targets and structure-activity relationships.