Phenomic assessment of genetic buffering by kinetic analysis of cell arrays

High-resolution yeast quiescence profiling in human-like media reveals complex influences of auxotrophy and nutrient availability

Yeast cells survive in stationary phase culture by entering quiescence, which is measured by colony-forming capacity upon nutrient re-exposure. Yeast chronological lifespan (CLS) studies, employing the comprehensive collection of gene knockout strains, have correlated weakly between independent laboratories, which is hypothesized to reflect differential interaction between the deleted genes, auxotrophy, media composition, and other assay conditions influencing quiescence. This hypothesis was investigated by high-throughput quiescence profiling of the parental prototrophic strain, from which the gene deletion strain libraries were constructed, and all possible auxotrophic allele combinations in that background. Defined media resembling human cell culture media promoted long-term quiescence and was used to assess effects of glucose, ammonium sulfate, auxotrophic nutrient availability, target of rapamycin signaling, and replication stress. Frequent, high-replicate measurements of colony-forming capacity from cultures aged past 60 days provided profiles of quiescence phenomena such as gasping and hormesis. Media acidification was assayed in parallel to assess correlation. Influences of leucine, methionine, glucose, and ammonium sulfate metabolism were clarified, and a role for lysine metabolism newly characterized, while histidine and uracil perturbations had less impact. Interactions occurred between glucose, ammonium sulfate, auxotrophy, auxotrophic nutrient limitation, aeration, TOR signaling, and/or replication stress. Weak correlation existed between media acidification and maintenance of quiescence. In summary, experimental factors, uncontrolled across previous genome-wide yeast CLS studies, influence quiescence and interact extensively, revealing quiescence as a complex metabolic and developmental process that should be studied in a prototrophic context, omitting ammonium sulfate from defined media, and employing highly replicable protocols.

A yeast phenomic model for the influence of Warburg metabolism on genetic buffering of doxorubicin

BACKGROUND

The influence of the Warburg phenomenon on chemotherapy response is unknown. Saccharomyces cerevisiae mimics the Warburg effect, repressing respiration in the presence of adequate glucose. Yeast phenomic experiments were conducted to assess potential influences of Warburg metabolism on gene-drug interaction underlying the cellular response to doxorubicin. Homologous genes from yeast phenomic and cancer pharmacogenomics data were analyzed to infer evolutionary conservation of gene-drug interaction and predict therapeutic relevance.

METHODS

Cell proliferation phenotypes (CPPs) of the yeast gene knockout/knockdown library were measured by quantitative high-throughput cell array phenotyping (Q-HTCP), treating with escalating doxorubicin concentrations under conditions of respiratory or glycolytic metabolism. Doxorubicin-gene interaction was quantified by departure of CPPs observed for the doxorubicin-treated mutant strain from that expected based on an interaction model. Recursive expectation-maximization clustering (REMc) and Gene Ontology (GO)-based analyses of interactions identified functional biological modules that differentially buffer or promote doxorubicin cytotoxicity with respect to Warburg metabolism. Yeast phenomic and cancer pharmacogenomics data were integrated to predict differential gene expression causally influencing doxorubicin anti-tumor efficacy.

RESULTS

Yeast compromised for genes functioning in chromatin organization, and several other cellular processes are more resistant to doxorubicin under glycolytic conditions. Thus, the Warburg transition appears to alleviate requirements for cellular functions that buffer doxorubicin cytotoxicity in a respiratory context. We analyzed human homologs of yeast genes exhibiting gene-doxorubicin interaction in cancer pharmacogenomics data to predict causality for differential gene expression associated with doxorubicin cytotoxicity in cancer cells. This analysis suggested conserved cellular responses to doxorubicin due to influences of homologous recombination, sphingolipid homeostasis, telomere tethering at nuclear periphery, actin cortical patch localization, and other gene functions.

CONCLUSIONS

Warburg status alters the genetic network required for yeast to buffer doxorubicin toxicity. Integration of yeast phenomic and cancer pharmacogenomics data suggests evolutionary conservation of gene-drug interaction networks and provides a new experimental approach to model their influence on chemotherapy response. Thus, yeast phenomic models could aid the development of precision oncology algorithms to predict efficacious cytotoxic drugs for cancer, based on genetic and metabolic profiles of individual tumors.

A Humanized Yeast Phenomic Model of Deoxycytidine Kinase to Predict Genetic Buffering of Nucleoside Analog Cytotoxicity

Knowledge about synthetic lethality can be applied to enhance the efficacy of anticancer therapies in individual patients harboring genetic alterations in their cancer that specifically render it vulnerable. We investigated the potential for high-resolution phenomic analysis in yeast to predict such genetic vulnerabilities by systematic, comprehensive, and quantitative assessment of drug–gene interaction for gemcitabine and cytarabine, substrates of deoxycytidine kinase that have similar molecular structures yet distinct antitumor efficacy. Human deoxycytidine kinase (dCK) was conditionally expressed in the Saccharomyces cerevisiae genomic library of knockout and knockdown (YKO/KD) strains, to globally and quantitatively characterize differential drug–gene interaction for gemcitabine and cytarabine. Pathway enrichment analysis revealed that autophagy, histone modification, chromatin remodeling, and apoptosis-related processes influence gemcitabine specifically, while drug–gene interaction specific to cytarabine was less enriched in gene ontology. Processes having influence over both drugs were DNA repair and integrity checkpoints and vesicle transport and fusion. Non-gene ontology (GO)-enriched genes were also informative. Yeast phenomic and cancer cell line pharmacogenomics data were integrated to identify yeast–human homologs with correlated differential gene expression and drug efficacy, thus providing a unique resource to predict whether differential gene expression observed in cancer genetic profiles are causal in tumor-specific responses to cytotoxic agents.

Gene-nutrient interaction markedly influences yeast chronological lifespan

PURPOSE

Research into the genetic mechanisms of aging has expanded rapidly over the past two decades. This has in part been the result of the use of model organisms (particularly yeast, worms and flies) and high-throughput technologies, combined with a growing interest in aging research. Despite this progress, widespread consensus regarding the pathways that are fundamental to the modulation of cellular aging and lifespan for all organisms has been limited due to discrepancies between different studies. We have compared results from published genome-wide, chronological lifespan (CLS) screens of individual gene deletion strains in Saccharomyces cerevisiae in order to identify gene deletion strains with consistent influences on longevity as possible indicators of fundamental aging processes from this single-celled, eukaryotic model organism.

METHODS

Three previous reports have described genetic modifiers of chronological aging in the budding yeast (S. cerevisiae) using the yeast gene deletion strain collection. We performed a comparison among the data sets using correlation and decile distribution analysis to describe concordance between screens and identify strains that consistently increased or decreased CLS. We used gene enrichment analysis in an effort to understand the biology underlying genes identified in multiple studies. We attempted to replicate the different experimental conditions employed by the screens to identify potential sources of variability in CLS worth further investigating.

RESULTS

Among 3209 strains present in all three screens, nine deletions strains were in common in the longest-lived decile (2.80%) and thirteen were in common in the shortest-lived decile (4.05%) of all three screens. Similarly, pairwise overlap between screens was low. When the same comparison was extended to three deciles to include more mutants studied in common between the three screens, enrichment of cellular processes based on gene ontology analysis in the long-lived strains remained very limited. To test the hypothesis that different parental strain auxotrophic requirements or media formulations employed by the respective genome-wide screens might contribute to the lack of concordance, different CLS assay conditions were assessed in combination with strains having different ploidy and auxotrophic requirements (all relevant to differences in the way the three genome-wide CLS screens were performed). This limited but systematic analysis of CLS with respect to auxotrophy, ploidy, and media revealed several instances of gene-nutrient interaction.

CONCLUSIONS

There is surprisingly little overlap between the results of three independently performed genome-wide screens of CLS in S. cerevisiae. However, differences in strain genetic background (ploidy and specific auxotrophic requirements) were present, as well as different media and experimental conditions (e.g., aeration and pooled vs. individual culturing), which, along with stochastic effects such as genetic drift or selection of secondary mutations that suppress the loss of function from gene deletion, could in theory account for some of the lack of consensus between results. Considering the lack of overlap in CLS phenotypes among the set of genes reported by all three screens, and the results of a CLS experiment that systematically tested (incorporating extensive controls) for interactions between variables existing between the screens, we propose that discrepancies can be reconciled through deeper understanding of the influence of cell intrinsic factors such as auxotrophic requirements ploidy status, extrinsic factors such as media composition and aeration, as well as interactions that may occur between them, for example as a result of different pooling vs. individually aging cultures. Such factors may have a more significant impact on CLS outcomes than previously realized. Future studies that systematically account for these contextual factors, and can thus clarify the interactions between genetic and nutrient factors that alter CLS phenotypes, should aid more complete understanding of the underlying biology so that genetic principles of CLS in yeast can be extrapolated to differential cellular aging observed in animal models.

Long-range coupling between the extracellular gates and the intracellular ATP binding domains of multidrug resistance protein pumps and cystic fibrosis transmembrane conductance regulator channels

The ABCC transporter subfamily includes pumps, the long and short multidrug resistance proteins (MRPs), and an ATP-gated anion channel, the cystic fibrosis transmembrane conductance regulator (CFTR). We show that despite their thermodynamic differences, these ABCC transporter subtypes use broadly similar mechanisms to couple their extracellular gates to the ATP occupancies of their cytosolic nucleotide binding domains. A conserved extracellular phenylalanine at this gate was a prime location for producing gain of function (GOF) mutants of a long MRP in yeast (Ycf1p cadmium transporter), a short yeast MRP (Yor1p oligomycin exporter), and human CFTR channels. Extracellular gate mutations rescued ATP binding mutants of the yeast MRPs and CFTR by increasing ATP sensitivity. Control ATPase-defective MRP mutants could not be rescued by this mechanism. A CFTR double mutant with an extracellular gate mutation plus a cytosolic GOF mutation was highly active (single-channel open probability >0.3) in the absence of ATP and protein kinase A, each normally required for CFTR activity. We conclude that all 3 ABCC transporter subtypes use similar mechanisms to couple their extracellular gates to ATP occupancy, and highly active CFTR channels that bypass defects in ATP binding or phosphorylation can be produced.

Yeast Phenomics: An Experimental Approach for Modeling Gene Interaction Networks that Buffer Disease

The genome project increased appreciation of genetic complexity underlying disease phenotypes: many genes contribute each phenotype and each gene contributes multiple phenotypes. The aspiration of predicting common disease in individuals has evolved from seeking primary loci to marginal risk assignments based on many genes. Genetic interaction, defined as contributions to a phenotype that are dependent upon particular digenic allele combinations, could improve prediction of phenotype from complex genotype, but it is difficult to study in human populations. High throughput, systematic analysis of S. cerevisiae gene knockouts or knockdowns in the context of disease-relevant phenotypic perturbations provides a tractable experimental approach to derive gene interaction networks, in order to deduce by cross-species gene homology how phenotype is buffered against disease-risk genotypes. Yeast gene interaction network analysis to date has revealed biology more complex than previously imagined. This has motivated the development of more powerful yeast cell array phenotyping methods to globally model the role of gene interaction networks in modulating phenotypes (which we call yeast phenomic analysis). The article illustrates yeast phenomic technology, which is applied here to quantify gene X media interaction at higher resolution and supports use of a human-like media for future applications of yeast phenomics for modeling human disease.