Genome-Wide Mapping of Binding Sites Reveals Multiple Biological Functions of the Transcription Factor Cst6p in Saccharomyces cerevisiae

Analysis of Pneumocystis Transcription Factor Evolution and Implications for Biology and Lifestyle

Pneumocystis jirovecii kills hundreds of thousands of immunocompromised patients each year. Yet many aspects of the biology of this obligate pathogen remain obscure because it is not possible to culture the fungus in vitro independently of its host. Consequently, our understanding of Pneumocystis pathobiology is heavily reliant upon bioinformatic inferences. We have exploited a powerful combination of genomic and phylogenetic approaches to examine the evolution of transcription factors in Pneumocystis species. We selected protein families (Pfam families) that correspond to transcriptional regulators and used bioinformatic approaches to compare these families in the seven Pneumocystis species that have been sequenced to date with those from other yeasts, other human and plant pathogens, and other obligate parasites. Some Pfam families of transcription factors have undergone significant reduction during their evolution in the Pneumocystis genus, and other Pfam families have been lost or appear to be in the process of being lost. Meanwhile, other transcription factor families have been retained in Pneumocystis species, and some even appear to have undergone expansion. On this basis, Pneumocystis species seem to have retained transcriptional regulators that control chromosome maintenance, ribosomal gene regulation, RNA processing and modification, and respiration. Meanwhile, regulators that promote the assimilation of alternative carbon sources, amino acid, lipid, and sterol biosynthesis, and iron sensing and homeostasis appear to have been lost. Our analyses of transcription factor retention, loss, and gain provide important insights into the biology and lifestyle of Pneumocystis. IMPORTANCE Pneumocystis jirovecii is a major fungal pathogen of humans that infects healthy individuals, colonizing the lungs of infants. In immunocompromised and transplant patients, this fungus causes life-threatening pneumonia, and these Pneumocystis infections remain among the most common and serious infections in HIV/AIDS patients. Yet we remain remarkably ignorant about the biology and epidemiology of Pneumocystis due to the inability to culture this fungus in vitro. Our analyses of transcription factor retentions, losses, and gains in sequenced Pneumocystis species provide valuable new views of their specialized biology, suggesting the retention of many metabolic and stress regulators and the loss of others that are essential in free-living fungi. Given the lack of in vitro culture methods for Pneumocystis, this powerful bioinformatic approach has advanced our understanding of the lifestyle of P. jirovecii and the nature of its dependence on the host for survival.

Association of ScV-LA Virus with Host Protein Metabolism Determined by Proteomics Analysis and Cognate RNA Sequencing

Saccharomyces yeasts are highly dispersed in the environment and microbiota of higher organisms. The yeast killing phenotype, encoded by the viral system, was discovered to be a significant property for host survival. Minor alterations in transcription patterns underpin the reciprocal relationship between LA and M viruses and their hosts, suggesting the fine-tuning of the transcriptional landscape. To uncover the principal targets of both viruses, we performed proteomics analysis of virus-enriched subsets of host proteins in virus type-specific manner. The essential pathways of protein metabolism-from biosynthesis and folding to degradation-were found substantially enriched in virus-linked subsets. The fractionation of viruses allowed separation of virus-linked host RNAs, investigated by high-content RNA sequencing. Ribosomal RNA was found to be inherently associated with LA-lus virus, along with other RNAs essential for ribosome biogenesis. This study provides a unique portrayal of yeast virions through the characterization of the associated proteome and cognate RNAs, and offers a background for understanding ScV-LA viral infection persistency.

Mediator Engineering of Saccharomyces cerevisiae To Improve Multidimensional Stress Tolerance

Saccharomyces cerevisiae is a well-performing workhorse in chemical production, which encounters complex environmental stresses during industrial processes. We constructed a multiple stress tolerance mutant, Med15^V76R/R84K, that was obtained by engineering the KIX domain of Mediator tail subunit Med15. Med15^V76R/R84K interacted with transcription factor Hap5 to improve ARV1 expression for sterol homeostasis for decreasing membrane fluidity and thereby enhancing acid tolerance. Med15^V76R/R84K interacted with transcription factor Mga2 to improve GIT1 expression for phospholipid biosynthesis for increasing membrane integrity and thereby improving oxidative tolerance. Med15^V76R/R84K interacted with transcription factor Aft1 to improve NFT1 expression for inorganic ion transport for reducing membrane permeability and thereby enhancing osmotic tolerance. Based on this Med15 mutation, Med15^V76R/R84K, the engineered S. cerevisiae strain, showed a 28.1% increase in pyruvate production in a 1.0-L bioreactor compared to that of S. cerevisiae with its native Med15. These results indicated that Mediator engineering provides a potential alternative for improving multidimensional stress tolerance in S. cerevisiae. IMPORTANCE This study identified the role of the KIX domain of Mediator tail subunit Med15 in response to acetic acid, H₂O₂, and NaCl in S. cerevisiae. Engineered KIX domain by protein engineering, the mutant strain Med15^V76R/R84K, increased multidimensional stress tolerance and pyruvate production compared with that of S. cerevisiae with its native Med15. The Med15^V76R/R84K could increase membrane related genes expression possibly by enhancing interaction with transcription factor to improve membrane physiological functions under stress conditions.

Genome-scale modeling of yeast metabolism: retrospectives and perspectives

Yeasts have been widely used for production of bread, beer and wine, as well as for production of bioethanol, but they have also been designed as cell factories to produce various chemicals, advanced biofuels and recombinant proteins. To systematically understand and rationally engineer yeast metabolism, genome-scale metabolic models (GEMs) have been reconstructed for the model yeast Saccharomyces cerevisiae and nonconventional yeasts. Here, we review the historical development of yeast GEMs together with their recent applications, including metabolic flux prediction, cell factory design, culture condition optimization and multi-yeast comparative analysis. Furthermore, we present an emerging effort, namely the integration of proteome constraints into yeast GEMs, resulting in models with improved performance. At last, we discuss challenges and perspectives on the development of yeast GEMs and the integration of proteome constraints.

Loss-of-Function ROX1 Mutations Suppress the Fluconazole Susceptibility of upc2AΔ Mutation in Candida glabrata, Implicating Additional Positive Regulators of Ergosterol Biosynthesis

Two of the major classes of antifungal drugs in clinical use target ergosterol biosynthesis. Despite its importance, our understanding of the transcriptional regulation of ergosterol biosynthesis genes in pathogenic fungi is essentially limited to the role of hypoxia and sterol-stress-induced transcription factors such as Upc2 and Upc2A as well as homologs of sterol response element binding (SREB) factors. To identify additional regulators of ergosterol biosynthesis in Candida glabrata, an important human fungal pathogen with reduced susceptibility to ergosterol biosynthesis inhibitors relative to other Candida spp., we used a serial passaging strategy to isolate suppressors of the fluconazole hypersusceptibility of a upc2AΔ deletion mutant. This led to the identification of loss-of-function mutations in two genes: ROX1, the homolog of a hypoxia gene transcriptional suppressor in Saccharomyces cerevisiae, and CST6, a transcription factor that is involved in the regulation of carbon dioxide response in C. glabrata. Here, we describe a detailed analysis of the genetic interaction of ROX1 and UPC2A. In the presence of fluconazole, loss of Rox1 function restores ERG11 expression to the upc2AΔ mutant and inhibits the expression of ERG3 and ERG6, leading to increased levels of ergosterol and decreased levels of the toxic sterol 14α methyl-ergosta-8,24(28)-dien-3β, 6α-diol, relative to the upc2AΔ mutant. Our observations establish that Rox1 is a negative regulator of ERG gene biosynthesis and indicate that a least one additional positive transcriptional regulator of ERG gene biosynthesis must be present in C. glabrata.

IMPORTANCECandida glabrata is one of the most important human fungal pathogens and has reduced susceptibility to azole-class inhibitors of ergosterol biosynthesis. Although ergosterol is the target of two of the three classes of antifungal drugs, relatively little is known about the regulation of this critical cellular pathway. Sterols are both essential components of the eukaryotic plasma membrane and potential toxins; therefore, sterol homeostasis is critical for cell function. Here, we identified two new negative regulators in C. glabrata of ergosterol (ERG) biosynthesis gene expression. Our results also indicate that in addition to Upc2A, the only known activator of ERG genes, additional positive regulators of this pathway must exist.

Uncovering Divergence in Gene Expression Regulation in the Adaptation of Yeast to Nitrogen Scarcity

Saccharomyces cerevisiae rewires its transcriptional output to survive stressful environments, such as nitrogen scarcity under fermentative conditions. Although divergence in nitrogen metabolism among natural yeast populations has been reported, the impact of regulatory genetic variants modulating gene expression and nitrogen consumption remains to be investigated. Here, we employed an F1 hybrid from two contrasting S. cerevisiae strains, providing a controlled genetic environment to map cis factors involved in the divergence of gene expression regulation in response to nitrogen scarcity. We used a dual approach to obtain genome-wide allele-specific profiles of chromatin accessibility, transcription factor binding, and gene expression through ATAC-seq (assay for transposase accessible chromatin) and RNA-seq (transcriptome sequencing). We observed large variability in allele-specific expression and accessibility between the two genetic backgrounds, with a third of these differences specific to a deficient nitrogen environment. Furthermore, we discovered events of allelic bias in gene expression correlating with allelic bias in transcription factor binding solely under nitrogen scarcity, where the majority of these transcription factors orchestrates the nitrogen catabolite repression regulatory pathway and demonstrates a cis × environment-specific response. Our approach allowed us to find cis variants modulating gene expression, chromatin accessibility, and allelic differences in transcription factor binding in response to low nitrogen culture conditions. IMPORTANCE Historically, coding variants were prioritized when searching for causal mechanisms driving adaptation of natural populations to stressful environments. However, the recent focus on noncoding variants demonstrated their ubiquitous role in adaptation. Here, we performed genome-wide regulatory variation profiles between two divergent yeast strains when facing nitrogen nutritional stress. The open chromatin availability of several regulatory regions changes in response to nitrogen scarcity. Importantly, we describe regulatory events that deviate between strains. Our results demonstrate a widespread variation in gene expression regulation between naturally occurring populations in response to stressful environments.

Acute ethanol stress induces sumoylation of conserved chromatin structural proteins in Saccharomyces cerevisiae

Stress is ubiquitous to life and can irreparably damage essential biomolecules and organelles in cells. To survive, organisms must sense and adapt to stressful conditions. One highly conserved adaptive stress response is through the posttranslational modification of proteins by the small ubiquitin-like modifier (SUMO). Here, we examine the effects of acute ethanol stress on protein sumoylation in the budding yeast Saccharomyces cerevisiae. We found that cells exhibit a transient sumoylation response after acute exposure to ≤7.5% vol/vol ethanol. By contrast, the sumoylation response becomes chronic at 10% ethanol exposure. Mass spectrometry analyses identified 18 proteins that are sumoylated after acute ethanol exposure, with 15 known to associate with chromatin. Upon further analysis, we found that the chromatin structural proteins Smc5 and Smc6 undergo ethanol-induced sumoylation that depends on the activity of the E3 SUMO ligase Mms21. Using cell-cycle arrest assays, we observed that Smc5 and Smc6 ethanol-induced sumoylation occurs during G1 and G2/M phases but not S phase. Acute ethanol exposure also resulted in the formation of Rad52 foci at levels comparable to Rad52 foci formation after exposure to the DNA alkylating agent methyl methanesulfonate (MMS). MMS exposure is known to induce the intra-S-phase DNA damage checkpoint via Rad53 phosphorylation, but ethanol exposure did not induce Rad53 phosphorylation. Ethanol abrogated the effect of MMS on Rad53 phosphorylation when added simultaneously. From these studies, we propose that acute ethanol exposure induces a change in chromatin leading to sumoylation of specific chromatin structural proteins.

A proteome-integrated, carbon source dependent genetic regulatory network in Saccharomyces cerevisiae

Integrated regulatory networks can be powerful tools to examine and test properties of cellular systems, such as modelling environmental effects on the molecular bioeconomy, where protein levels are altered in response to changes in growth conditions. Although extensive regulatory pathways and protein interaction data sets exist which represent such networks, few have formally considered quantitative proteomics data to validate and extend them. We generate and consider such data here using a label-free proteomics strategy to quantify alterations in protein abundance for S. cerevisiae when grown on minimal media using glucose, galactose, maltose and trehalose as sole carbon sources. Using a high quality-controlled subset of proteins observed to be differentially abundant, we constructed a proteome-informed network, comprising 1850 transcription factor interactions and 37 chaperone interactions, which defines the major changes in the cellular proteome when growing under different carbon sources. Analysis of the differentially abundant proteins involved in the regulatory network pointed to their significant roles in specific metabolic pathways and function, including glucose homeostasis, amino acid biosynthesis, and carbohydrate metabolic process. We noted strong statistical enrichment in the differentially abundant proteome of targets of known transcription factors associated with stress responses and altered carbon metabolism. This shows how such integrated analysis can lend further experimental support to annotated regulatory interactions, since the proteomic changes capture both magnitude and direction of gene expression change at the level of the affected proteins. Overall this study highlights the power of quantitative proteomics to help define regulatory systems pertinent to environmental conditions.

A Transcriptional Regulatory Map of Iron Homeostasis Reveals a New Control Circuit for Capsule Formation in Cryptococcus neoformans

Iron is essential for the growth of the human fungal pathogen Cryptococcus neoformans within the vertebrate host, and iron sensing contributes to the elaboration of key virulence factors, including the formation of the polysaccharide capsule. C. neoformans employs sophisticated iron acquisition and utilization systems governed by the transcription factors Cir1 and HapX. However, the details of the transcriptional regulatory networks that are governed by these transcription factors and connections to virulence remain to be defined. Here, we used chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) and transcriptome analysis (RNA-seq) to identify genes directly regulated by Cir1 and/or HapX in response to iron availability. Overall, 40 and 100 genes were directly regulated by Cir1, and 171 and 12 genes were directly regulated by HapX, under iron-limited and replete conditions, respectively. More specifically, we found that Cir1 directly controls the expression of genes required for iron acquisition and metabolism, and indirectly governs capsule formation by regulating specific protein kinases, a regulatory connection not previously revealed. HapX regulates the genes responsible for iron-dependent pathways, particularly under iron-depleted conditions. By analyzing target genes directly bound by Cir1 and HapX, we predicted the binding motifs for the transcription factors and verified that the purified proteins bind these motifs in vitro Furthermore, several direct target genes were coordinately and reciprocally regulated by Cir1 and HapX, suggesting that these transcription factors play conserved roles in the response to iron availability. In addition, biochemical analyses revealed that Cir1 and HapX are iron-containing proteins, implying that the regulatory networks of Cir1 and HapX may be influenced by the incorporation of iron into these proteins. Taken together, our identification of the genome-wide transcriptional networks provides a detailed understanding of the iron-related regulatory landscape, establishes a new connection between Cir1 and kinases that regulate capsule, and underpins genetic and biochemical analyses that reveal iron-sensing mechanisms for Cir1 and HapX in C. neoformans.

ChIP-exo analysis highlights Fkh1 and Fkh2 transcription factors as hubs that integrate multi-scale networks in budding yeast

The understanding of the multi-scale nature of molecular networks represents a major challenge. For example, regulation of a timely cell cycle must be coordinated with growth, during which changes in metabolism occur, and integrate information from the extracellular environment, e.g. signal transduction. Forkhead transcription factors are evolutionarily conserved among eukaryotes, and coordinate a timely cell cycle progression in budding yeast. Specifically, Fkh1 and Fkh2 are expressed during a lengthy window of the cell cycle, thus are potentially able to function as hubs in the multi-scale cellular environment that interlocks various biochemical networks. Here we report on a novel ChIP-exo dataset for Fkh1 and Fkh2 in both logarithmic and stationary phases, which is analyzed by novel and existing software tools. Our analysis confirms known Forkhead targets from available ChIP-chip studies and highlights novel ones involved in the cell cycle, metabolism and signal transduction. Target genes are analyzed with respect to their function, temporal expression during the cell cycle, correlation with Fkh1 and Fkh2 as well as signaling and metabolic pathways they occur in. Furthermore, differences in targets between Fkh1 and Fkh2 are presented. Our work highlights Forkhead transcription factors as hubs that integrate multi-scale networks to achieve proper timing of cell division in budding yeast.

Predictive models of eukaryotic transcriptional regulation reveals changes in transcription factor roles and promoter usage between metabolic conditions

Transcription factors (TF) are central to transcriptional regulation, but they are often studied in relative isolation and without close control of the metabolic state of the cell. Here, we describe genome-wide binding (by ChIP-exo) of 15 yeast TFs in four chemostat conditions that cover a range of metabolic states. We integrate this data with transcriptomics and six additional recently mapped TFs to identify predictive models describing how TFs control gene expression in different metabolic conditions. Contributions by TFs to gene regulation are predicted to be mostly activating, additive and well approximated by assuming linear effects from TF binding signal. Notably, using TF binding peaks from peak finding algorithms gave distinctly worse predictions than simply summing the low-noise and high-resolution TF ChIP-exo reads on promoters. Finally, we discover indications of a novel functional role for three TFs; Gcn4, Ert1 and Sut1 during nitrogen limited aerobic fermentation. In only this condition, the three TFs have correlated binding to a large number of genes (enriched for glycolytic and translation processes) and a negative correlation to target gene transcript levels.

Sumoylation of DNA-bound transcription factor Sko1 prevents its association with nontarget promoters

Sequence-specific transcription factors (TFs) represent one of the largest groups of proteins that is targeted for SUMO post-translational modification, in both yeast and humans. SUMO modification can have diverse effects, but recent studies showed that sumoylation reduces the interaction of multiple TFs with DNA in living cells. Whether this relates to a general role for sumoylation in TF binding site selection, however, has not been fully explored because few genome-wide studies aimed at studying such a role have been reported. To address this, we used genome-wide analysis to examine how sumoylation regulates Sko1, a yeast bZIP TF with hundreds of known binding sites. We find that Sko1 is sumoylated at Lys 567 and, although many of its targets are osmoresponse genes, the level of Sko1 sumoylation is not stress-regulated and the modification does not depend or impinge on its phosphorylation by the osmostress kinase Hog1. We show that Sko1 mutants that cannot bind DNA are not sumoylated, but attaching a heterologous DNA binding domain restores the modification, implicating DNA binding as a major determinant for Sko1 sumoylation. Genome-wide chromatin immunoprecipitation (ChIP-seq) analysis shows that a sumoylation-deficient Sko1 mutant displays increased occupancy levels at its numerous binding sites, which inhibits the recruitment of the Hog1 kinase to some induced osmostress genes. This strongly supports a general role for sumoylation in reducing the association of TFs with chromatin. Extending this result, remarkably, sumoylation-deficient Sko1 binds numerous additional promoters that are not normally regulated by Sko1 but contain sequences that resemble the Sko1 binding motif. Our study points to an important role for sumoylation in modulating the interaction of a DNA-bound TF with chromatin to increase the specificity of TF-DNA interactions.

Reconstruction of a Global Transcriptional Regulatory Network for Control of Lipid Metabolism in Yeast by Using Chromatin Immunoprecipitation with Lambda Exonuclease Digestion

To build transcription regulatory networks, transcription factor binding must be analyzed in cells grown under different conditions because their responses and targets differ depending on environmental conditions. We performed whole-genome analysis of the DNA binding of five Saccharomyces cerevisiae transcription factors involved in lipid metabolism, Ino2, Ino4, Hap1, Oaf1, and Pip2, in response to four different environmental conditions in chemostat cultures, which allowed us to keep the specific growth rate constant. Chromatin immunoprecipitation with lambda exonuclease digestion (ChIP-exo) enabled the detection of binding events at a high resolution. We discovered a large number of unidentified targets and thus expanded functions for each transcription factor (e.g., glutamate biosynthesis as a target of Oaf1 and Pip2). Moreover, condition-dependent binding of transcription factors in response to cell metabolic state (e.g., differential binding of Ino2 between fermentative and respiratory metabolic conditions) was clearly suggested. Combining the new binding data with previously published data from transcription factor deletion studies revealed the high complexity of the transcriptional regulatory network for lipid metabolism in yeast, which involves the combinatorial and complementary regulation by multiple transcription factors. We anticipate that our work will provide insights into transcription factor binding dynamics that will prove useful for the understanding of transcription regulatory networks. IMPORTANCE Transcription factors play a crucial role in the regulation of gene expression and adaptation to different environments. To better understand the underlying roles of these adaptations, we performed experiments that give us high-resolution binding of transcription factors to their targets. We investigated five transcription factors involved in lipid metabolism in yeast, and we discovered multiple novel targets and condition-specific responses that allow us to draw a better regulatory map of the lipid metabolism.

Interplay of Energetics and ER Stress Exacerbates Alzheimer's Amyloid-β (Aβ) Toxicity in Yeast

Alzheimer's disease (AD) is a progressive neurodegeneration. Oligomers of amyloid-β peptides (Aβ) are thought to play a pivotal role in AD pathogenesis, yet the mechanisms involved remain unclear. Two major isoforms of Aβ associated with AD are Aβ40 and Aβ42, the latter being more toxic and prone to form oligomers. Here, we took a systems biology approach to study two humanized yeast AD models which expressed either Aβ40 or Aβ42 in bioreactor cultures. Strict control of oxygen availability and culture pH, strongly affected chronological lifespan and reduced variations during cell growth. Reduced growth rates and biomass yields were observed upon Aβ42 expression, indicating a redirection of energy from growth to maintenance. Quantitative physiology analyses furthermore revealed reduced mitochondrial functionality and ATP generation in Aβ42 expressing cells, which matched with observed aberrant mitochondrial structures. Genome-wide expression level analysis showed that Aβ42 expression triggered strong ER stress and unfolded protein responses. Equivalent expression of Aβ40, however, induced only mild ER stress, which resulted in hardly affected physiology. Using AD yeast models in well-controlled cultures strengthened our understanding on how cells translate different Aβ toxicity signals into particular cell fate programs, and further enhance their potential as a discovery platform to identify possible therapies.

CO2 sensing in fungi: at the heart of metabolic signaling

Adaptation to the changing environmental CO₂ levels is essential for all living cells. In particular, microorganisms colonizing and infecting the human body are exposed to highly variable concentrations, ranging from atmospheric 0.04 to 5% and more in blood and specific host niches. Carbonic anhydrases are highly conserved metalloenzymes that enable fixation of CO₂ by its conversion into bicarbonate. This process is not only crucial to ensure the supply of adequate carbon amounts for cellular metabolism, but also contributes to several signaling processes in fungi, including morphology and communication. The fungal specific carbonic anhydrase gene NCE103 is transcribed in response to CO₂ availability. As recently shown, this regulation relies on the ATF/CREB transcription factor Cst6 and the AGC family protein kinase Sch9. Here, we review the regulatory mechanisms which control NCE103 expression in the model organism Saccharomyces cerevisiae and the pathogenic yeasts Candida albicans and Candida glabrata and discuss which additional factors might contribute in this novel CO₂ sensing cascade.

Systems Biology of Metabolism

Metabolism is highly complex and involves thousands of different connected reactions; it is therefore necessary to use mathematical models for holistic studies. The use of mathematical models in biology is referred to as systems biology. In this review, the principles of systems biology are described, and two different types of mathematical models used for studying metabolism are discussed: kinetic models and genome-scale metabolic models. The use of different omics technologies, including transcriptomics, proteomics, metabolomics, and fluxomics, for studying metabolism is presented. Finally, the application of systems biology for analyzing global regulatory structures, engineering the metabolism of cell factories, and analyzing human diseases is discussed.

Engineering of synthetic, stress-responsive yeast promoters

Advances in synthetic biology and our understanding of the rules of promoter architecture have led to the development of diverse synthetic constitutive and inducible promoters in eukaryotes and prokaryotes. However, the design of promoters inducible by specific endogenous or environmental conditions is still rarely undertaken. In this study, we engineered and characterized a set of strong, synthetic promoters for budding yeast Saccharomyces cerevisiae that are inducible under acidic conditions (pH ≤ 3). Using available expression and transcription factor binding data, literature on transcriptional regulation, and known rules of promoter architecture we improved the low-pH performance of the YGP1 promoter by modifying transcription factor binding sites in its upstream activation sequence. The engineering strategy outlined for the YGP1 promoter was subsequently applied to create a response to low pH in the unrelated CCW14 promoter. We applied our best promoter variants to low-pH fermentations, enabling ten-fold increased production of lactic acid compared to titres obtained with the commonly used, native TEF1 promoter. Our findings outline and validate a general strategy to iteratively design and engineer synthetic yeast promoters inducible to environmental conditions or stresses of interest.