Experimental evolution of S. cerevisiae for caffeine tolerance alters multidrug resistance and TOR signaling pathways.

Caffeine is a natural compound that inhibits the major cellular signaling regulator TOR, leading to widespread effects including growth inhibition. S. cerevisiae yeast can adapt to tolerate high concentrations of caffeine in coffee and cacao fermentations and in experimental systems. While many factors affecting caffeine tolerance and TOR signaling have been identified, further characterization of their interactions and regulation remain to be studied. We used experimental evolution of S. cerevisiae to study the genetic contributions to caffeine tolerance in yeast, through a collaboration between high school students evolving yeast populations coupled with further research exploration in university labs. We identified multiple evolved yeast populations with mutations in PDR1 and PDR5, which contribute to multidrug resistance, and showed that gain-of-function mutations in multidrug resistance family transcription factors PDR1, PDR3, and YRR1 differentially contribute to caffeine tolerance. We also identified loss-of-function mutations in TOR effectors SIT4, SKY1, and TIP41, and show that these mutations contribute to caffeine tolerance. These findings support the importance of both the multidrug resistance family and TOR signaling in caffeine tolerance, and can inform future exploration of networks affected by caffeine and other TOR inhibitors in model systems and industrial applications.

Functional Analysis of G6PD Variants Associated With Low G6PD Activity in the All of Us Research Program.

Glucose-6-phosphate dehydrogenase (G6PD) protects red blood cells against oxidative damage through regeneration of NADPH. Individuals with G6PD polymorphisms (variants) that produce an impaired G6PD enzyme are usually asymptomatic, but at risk of hemolytic anemia from oxidative stressors, including certain drugs and foods. Prevention of G6PD deficiency-related hemolytic anemia is achievable through G6PD genetic testing or whole-genome sequencing (WGS) to identify affected individuals who should avoid hemolytic triggers. However, accurately predicting the clinical consequence of G6PD variants is limited by over 800 G6PD variants which remain of uncertain significance. There also remains significant variability in which deficiency-causing variants are included in pharmacogenomic testing arrays across institutions: many panels only include c.202G>A, even though dozens of other variants can also cause G6PD deficiency. Here, we seek to improve G6PD genotype interpretation using data available in the All of Us Research Program and using a yeast functional assay. We confirm that G6PD coding variants are the main contributor to decreased G6PD activity, and that 13% of individuals in the All of Us data with deficiency-causing variants would be missed if only the c.202G>A variant were tested for. We expand clinical interpretation for G6PD variants of uncertain significance; reporting that c.595A>G, known as G6PD Dagua or G6PD Açores, and the newly identified variant c.430C>G, reduce activity sufficiently to lead to G6PD deficiency. We also provide evidence that five missense variants of uncertain significance are unlikely to lead to G6PD deficiency, since they were seen in hemi- or homozygous individuals without a reduction in G6PD activity. We also applied the new WHO guidelines and were able to classify two synonymous variants as WHO class C. We anticipate these results will improve the accuracy, and prompt increased use, of G6PD genetic tests through a more complete clinical interpretation of G6PD variants. As the All of Us data increases from 245,000 to 1 million participants, and additional functional assays are carried out, we expect this research to serve as a template to enable complete characterization of G6PD deficiency genotypes. With an increased number of interpreted variants, genetic testing of G6PD will be more informative for preemptively identifying individuals at risk for drug- or food-induced hemolytic anemia.

yEvo: A modular eukaryotic genetics and evolution research experience for high school students.

The resources for carrying out and analyzing microbial evolution experiments have become more accessible, making it possible to expand these studies beyond the research laboratory and into the classroom. We developed five connected, standards-aligned yeast evolution laboratory modules, called "yEvo," for high school students. The modules enable students to take agency in answering open-ended research questions. In Module 1, students evolve baker's yeast to tolerate an antifungal drug, and in subsequent modules, investigate how evolved yeasts adapted to this stressful condition at both the phenotype and genotype levels. We used pre- and post-surveys from 72 students at two different schools and post-interviews with students and teachers to assess our program goals and guide module improvement over 3 years. We measured changes in student conceptions, confidence in scientific practices, and interest in STEM careers. Students who participated in yEvo showed improvements in understanding of activity-specific concepts and reported increased confidence in designing a valid biology experiment. Student experimental data replicated literature findings and has led to new insights into antifungal resistance. The modules and provided materials, alongside "proof of concept" evaluation metrics, will serve as a model for other university researchers and K - 16 classrooms interested in engaging in open-ended research questions using yeast as a model system.

Systematic profiling of ale yeast protein dynamics across fermentation and repitching.

Studying the genetic and molecular characteristics of brewing yeast strains is crucial for understanding their domestication history and adaptations accumulated over time in fermentation environments, and for guiding optimizations to the brewing process itself. Saccharomyces cerevisiae (brewing yeast) is among the most profiled organisms on the planet, yet the temporal molecular changes that underlie industrial fermentation and beer brewing remain understudied. Here, we characterized the genomic makeup of a Saccharomyces cerevisiae ale yeast widely used in the production of Hefeweizen beers, and applied shotgun mass spectrometry to systematically measure the proteomic changes throughout 2 fermentation cycles which were separated by 14 rounds of serial repitching. The resulting brewing yeast proteomics resource includes 64,740 protein abundance measurements. We found that this strain possesses typical genetic characteristics of Saccharomyces cerevisiae ale strains and displayed progressive shifts in molecular processes during fermentation based on protein abundance changes. We observed protein abundance differences between early fermentation batches compared to those separated by 14 rounds of serial repitching. The observed abundance differences occurred mainly in proteins involved in the metabolism of ergosterol and isobutyraldehyde. Our systematic profiling serves as a starting point for deeper characterization of how the yeast proteome changes during commercial fermentations and additionally serves as a resource to guide fermentation protocols, strain handling, and engineering practices in commercial brewing and fermentation environments. Finally, we created a web interface (https://brewing-yeast-proteomics.ccbb.utexas.edu/) to serve as a valuable resource for yeast geneticists, brewers, and biochemists to provide insights into the global trends underlying commercial beer production.

Systematic Profiling of Ale Yeast Protein Dynamics across Fermentation and Repitching.

Studying the genetic and molecular characteristics of brewing yeast strains is crucial for understanding their domestication history and adaptations accumulated over time in fermentation environments, and for guiding optimizations to the brewing process itself. Saccharomyces cerevisiae (brewing yeast) is amongst the most profiled organisms on the planet, yet the temporal molecular changes that underlie industrial fermentation and beer brewing remain understudied. Here, we characterized the genomic makeup of a Saccharomyces cerevisiae ale yeast widely used in the production of Hefeweizen beers, and applied shotgun mass spectrometry to systematically measure the proteomic changes throughout two fermentation cycles which were separated by 14 rounds of serial repitching. The resulting brewing yeast proteomics resource includes 64,740 protein abundance measurements. We found that this strain possesses typical genetic characteristics of Saccharomyces cerevisiae ale strains and displayed progressive shifts in molecular processes during fermentation based on protein abundance changes. We observed protein abundance differences between early fermentation batches compared to those separated by 14 rounds of serial repitching. The observed abundance differences occurred mainly in proteins involved in the metabolism of ergosterol and isobutyraldehyde. Our systematic profiling serves as a starting point for deeper characterization of how the yeast proteome changes during commercial fermentations and additionally serves as a resource to guide fermentation protocols, strain handling, and engineering practices in commercial brewing and fermentation environments. Finally, we created a web interface (https://brewing-yeast-proteomics.ccbb.utexas.edu/) to serve as a valuable resource for yeast geneticists, brewers, and biochemists to provide insights into the global trends underlying commercial beer production.

Caffeine-tolerant mutations selected through an at-home yeast experimental evolution teaching lab.

yEvo is a curriculum for high school students centered around evolution experiments in S. cerevisiae . To adapt the curriculum for remote instruction, we created a new protocol to evolve non-engineered yeast in the presence of caffeine. Evolved strains had increased caffeine tolerance and distinct colony morphologies. Many possessed copy number variations, transposon insertions, and mutations affecting genes with known relationships to caffeine and TOR signaling - which is inhibited by caffeine - and in other genes not previously connected with caffeine. This demonstrates that our accessible, at-home protocol is sufficient to permit novel insights into caffeine tolerance.

Caffeine-tolerant mutations selected through an at-home yeast experimental evolution teaching lab.

yEvo is a curriculum for high school students centered around evolution experiments in S. cerevisiae . To adapt the curriculum for remote instruction, we created a new protocol to evolve non-GMO yeast in the presence of caffeine. Evolved strains had increased caffeine tolerance and distinct colony morphologies. Many possessed copy number variations, transposon insertions, and mutations affecting genes with known relationships to caffeine and TOR signaling - which is inhibited by caffeine - and in other genes not previously connected with caffeine. This demonstrates that our accessible, at-home protocol is sufficient to permit novel insights into caffeine tolerance.

Functional interpretation, cataloging, and analysis of 1,341 glucose-6-phosphate dehydrogenase variants.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency affects over 500 million individuals who can experience anemia in response to oxidative stressors such as certain foods and drugs. Recently, the World Health Organization (WHO) called for revisiting G6PD variant classification as a priority to implement genetic medicine in low- and middle-income countries. Toward this goal, we sought to collect reports of G6PD variants and provide interpretations. We identified 1,341 G6PD variants in population and clinical databases. Using the ACMG standards and guidelines for the interpretation of sequence variants, we provided interpretations for 268 variants, including 186 variants that were not reported or of uncertain significance in ClinVar, bringing the total number of variants with non-conflicting interpretations to 400. For 414 variants with functional or clinical data, we analyzed associations between activity, stability, and current classification systems, including the new 2022 WHO classification. We corroborated known challenges with classification systems, including phenotypic variation, emphasizing the importance of comparing variant effects across individuals and studies. Biobank data made available by All of Us illustrate the benefit of large-scale sequencing and phenotyping by adding additional support connecting variants to G6PD-deficient anemia. By leveraging available data and interpretation guidelines, we created a repository for information on G6PD variants and nearly doubled the number of variants with clinical interpretations. These tools enable better interpretation of G6PD variants for the implementation of genetic medicine.

JAK-STAT Signaling in Inflammatory Breast Cancer Enables Chemotherapy-Resistant Cell States.

Inflammatory breast cancer (IBC) is a difficult-to-treat disease with poor clinical outcomes due to high risk of metastasis and resistance to treatment. In breast cancer, CD44+CD24- cells possess stem cell-like features and contribute to disease progression, and we previously described a CD44+CD24-pSTAT3+ breast cancer cell subpopulation that is dependent on JAK2/STAT3 signaling. Here we report that CD44+CD24- cells are the most frequent cell type in IBC and are commonly pSTAT3+. Combination of JAK2/STAT3 inhibition with paclitaxel decreased IBC xenograft growth more than either agent alone. IBC cell lines resistant to paclitaxel and doxorubicin were developed and characterized to mimic therapeutic resistance in patients. Multi-omic profiling of parental and resistant cells revealed enrichment of genes associated with lineage identity and inflammation in chemotherapy-resistant derivatives. Integrated pSTAT3 chromatin immunoprecipitation sequencing and RNA sequencing (RNA-seq) analyses showed pSTAT3 regulates genes related to inflammation and epithelial-to-mesenchymal transition (EMT) in resistant cells, as well as PDE4A, a cAMP-specific phosphodiesterase. Metabolomic characterization identified elevated cAMP signaling and CREB as a candidate therapeutic target in IBC. Investigation of cellular dynamics and heterogeneity at the single cell level during chemotherapy and acquired resistance by CyTOF and single cell RNA-seq identified mechanisms of resistance including a shift from luminal to basal/mesenchymal cell states through selection for rare preexisting subpopulations or an acquired change. Finally, combination treatment with paclitaxel and JAK2/STAT3 inhibition prevented the emergence of the mesenchymal chemo-resistant subpopulation. These results provide mechanistic rational for combination of chemotherapy with inhibition of JAK2/STAT3 signaling as a more effective therapeutic strategy in IBC. SIGNIFICANCE: Chemotherapy resistance in inflammatory breast cancer is driven by the JAK2/STAT3 pathway, in part via cAMP/PKA signaling and a cell state switch, which can be overcome using paclitaxel combined with JAK2 inhibitors.

yEvo: experimental evolution in high school classrooms selects for novel mutations that impact clotrimazole resistance in Saccharomyces cerevisiae.

Antifungal resistance in pathogenic fungi is a growing global health concern. Nonpathogenic laboratory strains of Saccharomyces cerevisiae are an important model for studying mechanisms of antifungal resistance that are relevant to understanding the same processes in pathogenic fungi. We have developed a series of laboratory modules in which high school students used experimental evolution to study antifungal resistance by isolating azole-resistant S. cerevisiae mutants and examining the genetic basis of resistance. We have sequenced 99 clones from these experiments and found that all possessed mutations previously shown to impact azole resistance, validating our approach. We additionally found recurrent mutations in an mRNA degradation pathway and an uncharacterized mitochondrial protein (Csf1) that have possible mechanistic connections to azole resistance. The scale of replication in this initiative allowed us to identify candidate epistatic interactions, as evidenced by pairs of mutations that occur in the same clone more frequently than expected by chance (positive epistasis) or less frequently (negative epistasis). We validated one of these pairs, a negative epistatic interaction between gain-of-function mutations in the multidrug resistance transcription factors Pdr1 and Pdr3. This high school-university collaboration can serve as a model for involving members of the broader public in the scientific process to make meaningful discoveries in biomedical research.

PI3K drives the de novo synthesis of coenzyme A from vitamin B5.

In response to hormones and growth factors, the class I phosphoinositide-3-kinase (PI3K) signalling network functions as a major regulator of metabolism and growth, governing cellular nutrient uptake, energy generation, reducing cofactor production and macromolecule biosynthesis1. Many of the driver mutations in cancer with the highest recurrence, including in receptor tyrosine kinases, Ras, PTEN and PI3K, pathologically activate PI3K signalling2,3. However, our understanding of the core metabolic program controlled by PI3K is almost certainly incomplete. Here, using mass-spectrometry-based metabolomics and isotope tracing, we show that PI3K signalling stimulates the de novo synthesis of one of the most pivotal metabolic cofactors: coenzyme A (CoA). CoA is the major carrier of activated acyl groups in cells4,5 and is synthesized from cysteine, ATP and the essential nutrient vitamin B5 (also known as pantothenate)6,7. We identify pantothenate kinase 2 (PANK2) and PANK4 as substrates of the PI3K effector kinase AKT8. Although PANK2 is known to catalyse the rate-determining first step of CoA synthesis, we find that the minimally characterized but highly conserved PANK49 is a rate-limiting suppressor of CoA synthesis through its metabolite phosphatase activity. Phosphorylation of PANK4 by AKT relieves this suppression. Ultimately, the PI3K-PANK4 axis regulates the abundance of acetyl-CoA and other acyl-CoAs, CoA-dependent processes such as lipid metabolism and proliferation. We propose that these regulatory mechanisms coordinate cellular CoA supplies with the demands of hormone/growth-factor-driven or oncogene-driven metabolism and growth.

Measuring Pharmacogene Variant Function at Scale Using Multiplexed Assays.

As costs of next-generation sequencing decrease, identification of genetic variants has far outpaced our ability to understand their functional consequences. This lack of understanding is a central challenge to a key promise of pharmacogenomics: using genetic information to guide drug selection and dosing. Recently developed multiplexed assays of variant effect enable experimental measurement of the function of thousands of variants simultaneously. Here, we describe multiplexed assays that have been performed on nearly 25,000 variants in eight key pharmacogenes (ADRB2, CYP2C9, CYP2C19, NUDT15, SLCO1B1, TMPT, VKORC1, and the LDLR promoter), discuss advances in experimental design, and explore key challenges that must be overcome to maximize the utility of multiplexed functional data.

FGFR-inhibitor-mediated dismissal of SWI/SNF complexes from YAP-dependent enhancers induces adaptive therapeutic resistance.

How cancer cells adapt to evade the therapeutic effects of drugs targeting oncogenic drivers is poorly understood. Here we report an epigenetic mechanism leading to the adaptive resistance of triple-negative breast cancer (TNBC) to fibroblast growth factor receptor (FGFR) inhibitors. Prolonged FGFR inhibition suppresses the function of BRG1-dependent chromatin remodelling, leading to an epigenetic state that derepresses YAP-associated enhancers. These chromatin changes induce the expression of several amino acid transporters, resulting in increased intracellular levels of specific amino acids that reactivate mTORC1. Consistent with this mechanism, addition of mTORC1 or YAP inhibitors to FGFR blockade synergistically attenuated the growth of TNBC patient-derived xenograft models. Collectively, these findings reveal a feedback loop involving an epigenetic state transition and metabolic reprogramming that leads to adaptive therapeutic resistance and provides potential therapeutic strategies to overcome this mechanism of resistance.

Metabolic pathway alterations in microvascular endothelial cells in response to hypoxia.

The vasculature within a tumor is highly disordered both structurally and functionally. Endothelial cells that comprise the vasculature are poorly connected causing vessel leakage and exposing the endothelium to a hypoxic microenvironment. Therefore, most anti-angiogenic therapies are generally inefficient and result in acquired resistance to increased hypoxia due to elimination of the vasculature. Recent studies have explored the efficacy of targeting metabolic pathways in tumor cells in combination with anti-angiogenic therapy. However, the metabolic alterations of endothelial cells in response to hypoxia have been relatively unexplored. Here, we measured polar metabolite levels in microvascular endothelial cells exposed to short- and long-term hypoxia with the goal of identifying metabolic vulnerabilities that can be targeted to normalize tumor vasculature and improve drug delivery. We found that many amino acid-related metabolites were altered by hypoxia exposure, especially within alanine-aspartate-glutamate, serine-threonine, and cysteine-methionine metabolism. Additionally, there were significant changes in de novo pyrimidine synthesis as well as glutathione and taurine metabolism. These results provide key insights into the metabolic alterations that occur in endothelial cells in response to hypoxia, which serve as a foundation for future studies to develop therapies that lead to vessel normalization and more efficient drug delivery.

The INPP4B Tumor Suppressor Modulates EGFR Trafficking and Promotes Triple-Negative Breast Cancer.

Inactivation of the tumor suppressor lipid phosphatase INPP4B is common in triple-negative breast cancer (TNBC). We generated a genetically engineered TNBC mouse model deficient in INPP4B. We found a dose-dependent increase in tumor incidence in INPP4B homozygous and heterozygous knockout mice compared with wild-type (WT), supporting a role for INPP4B as a tumor suppressor in TNBC. Tumors derived from INPP4B knockout mice are enriched for AKT and MEK gene signatures. Consequently, mice with INPP4B deficiency are more sensitive to PI3K or MEK inhibitors compared with WT mice. Mechanistically, we found that INPP4B deficiency increases PI(3,4)P2 levels in endocytic vesicles but not at the plasma membrane. Moreover, INPP4B loss delays degradation of EGFR and MET, while promoting recycling of receptor tyrosine kinases (RTK), thus enhancing the duration and amplitude of signaling output upon growth factor stimulation. Therefore, INPP4B inactivation in TNBC promotes tumorigenesis by modulating RTK recycling and signaling duration. SIGNIFICANCE: Inactivation of the lipid phosphatase INPP4B is frequent in TNBC. Using a genetically engineered mouse model, we show that INPP4B functions as a tumor suppressor in TNBC. INPP4B regulates RTK trafficking and degradation, such that loss of INPP4B prolongs both PI3K and ERK activation.This article is highlighted in the In This Issue feature, p. 1079.

Inhibition of the polyamine synthesis enzyme ornithine decarboxylase sensitizes triple-negative breast cancer cells to cytotoxic chemotherapy.

Treatment of patients with triple-negative breast cancer (TNBC) is limited by a lack of effective molecular therapies targeting this disease. Recent studies have identified metabolic alterations in cancer cells that can be targeted to improve responses to standard-of-care chemotherapy regimens. Using MDA-MB-468 and SUM-159PT TNBC cells, along with LC-MS/MS and HPLC metabolomics profiling, we found here that exposure of TNBC cells to the cytotoxic chemotherapy drugs cisplatin and doxorubicin alter arginine and polyamine metabolites. This alteration was because of a reduction in the levels and activity of a rate-limiting polyamine biosynthetic enzyme, ornithine decarboxylase (ODC). Using gene silencing and inhibitor treatments, we determined that the reduction in ODC was mediated by its negative regulator antizyme, targeting ODC to the proteasome for degradation. Treatment with the ODC inhibitor difluoromethylornithine (DFMO) sensitized TNBC cells to chemotherapy, but this was not observed in receptor-positive breast cancer cells. Moreover, TNBC cell lines had greater sensitivity to single-agent DFMO, and ODC levels were elevated in TNBC patient samples. The alterations in polyamine metabolism in response to chemotherapy, as well as DFMO-induced preferential sensitization of TNBC cells to chemotherapy, reported here suggest that ODC may be a targetable metabolic vulnerability in TNBC.

Oncogenic PI3K promotes methionine dependency in breast cancer cells through the cystine-glutamate antiporter xCT.

The precursor homocysteine is metabolized either through the methionine cycle to produce methionine or through the transsulfuration pathway to synthesize cysteine. Alternatively, cysteine can be obtained through uptake of its oxidized form, cystine. Many cancer cells exhibit methionine dependency such that their proliferation is impaired in growth media in which methionine is replaced by homocysteine. We showed that oncogenic PIK3CA and decreased expression of SLC7A11, a gene that encodes a cystine transporter also known as xCT, correlated with increased methionine dependency in breast cancer cells. Oncogenic PIK3CA was sufficient to confer methionine dependency to mammary epithelial cells, partly by decreasing cystine uptake through the transcriptional and posttranslational inhibition of xCT. Manipulation of xCT activity altered the proliferation of breast cancer cells in methionine-deficient, homocysteine-containing media, suggesting that it functionally contributed to methionine dependency. We propose that concurrent with decreased cystine uptake through xCT, PIK3CA mutant cells use homocysteine through the transsulfuration pathway to synthesize cysteine. Consequently, less homocysteine is available to produce methionine, contributing to methionine dependency. These results indicate that oncogenic PIK3CA alters methionine and cysteine utilization, partly by inhibiting xCT to contribute to the methionine dependency phenotype in breast cancer cells.

Nonessential amino acid metabolism in breast cancer.

Interest in studying cancer metabolism has risen in recent years, as it has become evident that the relationship between cancer and metabolic pathways could reveal novel biomarkers and therapeutic targets. Metabolic starvation therapy is particularly promising due to its low toxicity. Nonessential amino acids are promising metabolites for such therapy because they become essential in many tumor cells, including breast cancer cells. This review will focus on four nonessential amino acid metabolism pathways: glutamine-glutamate, serine-glycine, cysteine, and arginine-proline metabolism. Recent studies of these amino acids have revealed metabolic enzymes that have the potential to be effective as cancer therapy targets or biomarkers for response to metabolic starvation therapy. The review will also discuss features of nonessential amino acid metabolism that merit further investigation to determine their relevancy to breast cancer treatment.