Structure and function of a transcriptional network activated by the MAPK Hog1

Shared and redundant proteins coordinate signal cross-talk between MAPK pathways in yeast

All cells must detect, interpret, and adapt to multiple and concurrent stimuli. While signaling pathways are highly specialized, different pathways often share components or have components with overlapping functions. In the yeast Saccharomyces cerevisiae, the high osmolarity glycerol (HOG) pathway has two seemingly redundant branches, mediated by Sln1 and Sho1. Both branches are activated by osmotic pressure, leading to phosphorylation of the MAPKs Hog1 and Kss1. The mating pathway is activated by pheromone, leading to phosphorylation of the MAPKs Fus3 and Kss1. Given that Kss1 is shared by the two pathways, we investigated its role in signal coordination. We activated both pathways with a combination of salt and pheromone, in cells lacking the shared MAPK and in cells lacking either of the redundant branches of the HOG pathway. By systematically evaluating MAPK activation, translocation, and transcription programs, we determined that Sho1 mediates cross talk between the HOG and mating pathways and does so through Kss1. Further, we show that Kss1 initiates a transcriptional program that is distinct from that induced by Hog1 and Fus3. Our findings reveal how redundant and shared components coordinate concurrent signals and thereby adapt to sudden environmental changes.

Unraveling the interplay between a small RNA and RNase E in bacteria

Small RNAs (sRNAs) are major regulators of gene expression in bacteria, exerting their regulation primarily via base pairing with their target transcripts and modulating translation. Accumulating evidence suggest that sRNAs can also affect the stability of their target transcripts by altering their accessibility to endoribonucleases. Yet, the effects of sRNAs on transcript stability and the mechanisms underlying them have not been studied in wide scale. Here we employ large-scale RNA-seq-based methodologies in the model bacterium Escherichia coli to quantitatively study the functional interaction between a sRNA and an endoribonuclease in regulating gene expression, using the well-established sRNA, GcvB, and the major endoribonuclease, RNase E. Studying single and double mutants of gcvB and rne and analysing their RNA-seq results by the Double Mutant Cycle approach, we infer distinct modes of the interplay between GcvB and RNase E. Transcriptome-wide mapping of RNase E cleavage sites provides further support to the results of the RNA-seq analysis, identifying cleavage sites in targets in which the functional interaction between GcvB and RNase E is evident. Together, our results indicate that the most dominant mode of GcvB-RNase E functional interaction is GcvB enhancement of RNase E cleavage, which varies in its magnitude between different targets.

Mechanism of benzoxazinoids affecting the growth and development of Fusarium oxysporum f. sp. fabae

Continuous cropping of faba bean (Vicia faba L.) has led to a high incidence of wilt disease. The implementation of an intercropping system involving wheat and faba bean can effectively control the propagation of faba bean wilt disease. To investigate the mechanisms of wheat in mitigating faba bean wilt disease in a wheat-faba bean intercropping system. A comprehensive investigation was conducted to assess the temporal variations in Fusarium oxysporum f. sp. fabae (FOF) on the chemotaxis of benzoxazinoids (BXs) and wheat root through indoor culture tests. The effects of BXs on FOF mycelial growth, spore germination, spore production, and electrical conductivity were examined. The influence of BXs on the ultrastructure of FOF was investigated through transmission electron microscopy. Eukaryotic mRNA sequencing was utilized to analyze the differentially expressed genes in FOF upon treatment with BXs. FOF exhibited a significant positive chemotactic effect on BXs in wheat roots and root secretions. BXs possessed the potential to exert significant allelopathic effects on the mycelial growth, spore germination, and sporulation of FOF. In addition, BXs demonstrated a remarkable ability to disrupt the structural integrity and stability of the membrane and cell wall of the FOF mycelia. BXs possessed the capability of posing threats to the integrity and stability of the cell membrane and cell wall. This ultimately resulted in physiological dysfunction, effectively inhibiting the regular growth and developmental processes of the FOF.

Implication of polymerase recycling for nascent transcript quantification by live cell imaging

Transcription enables the production of RNA from a DNA template. Due to the highly dynamic nature of transcription, live-cell imaging methods play a crucial role in measuring the kinetics of this process. For instance, transcriptional bursts have been visualized using fluorescent phage-coat proteins that associate tightly with messenger RNA (mRNA) stem loops formed on nascent transcripts. To convert the signal emanating from a transcription site into meaningful estimates of transcription dynamics, the influence of various parameters on the measured signal must be evaluated. Here, the effect of gene length on the intensity of the transcription site focus was analyzed. Intuitively, a longer gene can support a larger number of transcribing polymerases, thus leading to an increase in the measured signal. However, measurements of transcription induced by hyper-osmotic stress responsive promoters display independence from gene length. A mathematical model of the stress-induced transcription process suggests that the formation of gene loops that favor the recycling of polymerase from the terminator to the promoter can explain the observed behavior. One experimentally validated prediction from this model is that the amount of mRNA produced from a short gene should be higher than for a long one as the density of active polymerase on the short gene will be increased by polymerase recycling. Our data suggest that this recycling contributes significantly to the expression output from a gene and that polymerase recycling is modulated by the promoter identity and the cellular state.

The cAMP-PKA signalling crosstalks with CWI and HOG-MAPK pathways in yeast cell response to osmotic and thermal stress

The yeast Saccharomyces cerevisiae is widely used in food and non-food industries. During industrial fermentation yeast strains are exposed to fluctuations in oxygen concentration, osmotic pressure, pH, ethanol concentration, nutrient availability and temperature. Fermentation performance depends on the ability of the yeast strains to adapt to these changes. Suboptimal conditions trigger responses to the external stimuli to allow homeostasis to be maintained. Stress-specific signalling pathways are activated to coordinate changes in transcription, translation, protein function, and metabolic fluxes while a transient arrest of growth and cell cycle progression occur. cAMP-PKA, HOG-MAPK and CWI signalling pathways are turned on during stress response. Comprehension of the mechanisms involved in the responses and in the adaptation to these stresses during fermentation is key to improving this industrial process. The scope of this review is to outline the advancement of knowledge about the cAMP-PKA signalling and the crosstalk of this pathway with the CWI and HOG-MAPK cascades in response to the environmental challenges heat and hyperosmotic stress.

Deletion of the transcription factors Hsf1, Msn2 and Msn4 in yeast uncovers transcriptional reprogramming in response to proteotoxic stress

The response to proteotoxic stresses such as heat shock allows organisms to maintain protein homeostasis under changing environmental conditions. We asked what happens if an organism can no longer react to cytosolic proteotoxic stress. To test this, we deleted or depleted, either individually or in combination, the stress-responsive transcription factors Msn2, Msn4, and Hsf1 in Saccharomyces cerevisiae. Our study reveals a combination of survival strategies, which together protect essential proteins. Msn2 and 4 broadly reprogram transcription, triggering the response to oxidative stress, as well as biosynthesis of the protective sugar trehalose and glycolytic enzymes, while Hsf1 mainly induces the synthesis of molecular chaperones and reverses the transcriptional response upon prolonged mild heat stress (adaptation).

Integrated genomics and phenotype microarray analysis of Saccharomyces cerevisiae industrial strains for rice wine fermentation and recombinant protein production

The industrial potential of Saccharomyces cerevisiae has extended beyond its traditional use in fermentation to various applications, including recombinant protein production. Herein, comparative genomics was performed with three industrial S. cerevisiae strains and revealed a heterozygous diploid genome for the 98-5 and KSD-YC strains (exploited for rice wine fermentation) and a haploid genome for strain Y2805 (used for recombinant protein production). Phylogenomic analysis indicated that Y2805 was closely associated with the reference strain S288C, whereas KSD-YC and 98-5 were grouped with Asian and European wine strains, respectively. Particularly, a single nucleotide polymorphism (SNP) in FDC1, involved in the biosynthesis of 4-vinylguaiacol (4-VG, a phenolic compound with a clove-like aroma), was found in KSD-YC, consistent with its lack of 4-VG production. Phenotype microarray (PM) analysis showed that KSD-YC and 98-5 displayed broader substrate utilization than S288C and Y2805. The SNPs detected by genome comparison were mapped to the genes responsible for the observed phenotypic differences. In addition, detailed information on the structural organization of Y2805 selection markers was validated by Sanger sequencing. Integrated genomics and PM analysis elucidated the evolutionary history and genetic diversity of industrial S. cerevisiae strains, providing a platform to improve fermentation processes and genetic manipulation.

A module centered on the transcription factor Msn2 from arbuscular mycorrhizal fungus Rhizophagus irregularis regulates drought stress tolerance in the host plant

Arbuscular mycorrhizal (AM) fungi can form mutualistic endosymbiosis with > 70% of land plants for obtaining fatty acids and sugars, in return, AM fungi promote plant nutrients and water acquisition to enhance plant fitness. However, how AM fungi orchestrate its own signaling components in response to drought stress remains elusive. Here, we identify a transcription factor containing C2H2 zinc finger domains, RiMsn2 from Rhizophagus irregularis. To characterize the RiMsn2, we combined heterologous expression, subcellular localization in yeasts, and biochemical and molecular studies with reverse genetics approaches during the in planta phase. The results indicate that RiMsn2 is highly conserved across AM fungal species and induced during the early stages of symbiosis. It is significantly upregulated in mycorrhizal roots under severe drought conditions. The nucleus-localized RiMsn2 regulates osmotic homeostasis and trehalose contents of yeasts. Importantly, gene silencing analyses indicate that RiMsn2 is essential for arbuscule formation and enhances plant tolerance to drought stress. Results from yeasts and biochemical experiments suggest that the RiHog1-RiMsn2-STREs module controls the drought stress-responsive genes in AM fungal symbiont. In conclusion, our findings reveal that a module centered on the transcriptional activator RiMsn2 from AM fungus regulates drought stress tolerance in host plant.

The architecture of binding cooperativity between densely bound transcription factors

The binding of transcription factors (TFs) along genomes is restricted to a subset of sites containing their preferred motifs. TF-binding specificity is often attributed to the co-binding of interacting TFs; however, apart from specific examples, this model remains untested. Here, we define dependencies among budding yeast TFs that localize to overlapping promoters by profiling the genome-wide consequences of co-depleting multiple TFs. We describe unidirectional interactions, revealing Msn2 as a central factor allowing TF binding at its target promoters. By contrast, no case of mutual cooperation was observed. Particularly, Msn2 retained binding at its preferred promoters upon co-depletion of fourteen similarly bound TFs. Overall, the consequences of TF co-depletions were moderate, limited to a subset of promoters, and failed to explain the role of regions outside the DNA-binding domain in directing TF-binding preferences. Our results call for re-evaluating the role of cooperative interactions in directing TF-binding preferences.

Discovery and Validation of Context-Dependent Synthetic Mammalian Promoters

Cellular transcription enables cells to adapt to various stimuli and maintain homeostasis. Transcription factors bind to transcription response elements (TREs) in gene promoters, initiating transcription. Synthetic promoters, derived from natural TREs, can be engineered to control exogenous gene expression using endogenous transcription machinery. This technology has found extensive use in biological research for applications including reporter gene assays, biomarker development, and programming synthetic circuits in living cells. However, a reliable and precise method for selecting minimally-sized synthetic promoters with desired background, amplitude, and stimulation response profiles has been elusive. In this study, we introduce a massively parallel reporter assay library containing 6184 synthetic promoters, each less than 250 bp in length. This comprehensive library allows for rapid identification of promoters with optimal transcriptional output parameters across multiple cell lines and stimuli. We showcase this library's utility to identify promoters activated in unique cell types, and in response to metabolites, mitogens, cellular toxins, and agonism of both aminergic and non-aminergic GPCRs. We further show these promoters can be used in luciferase reporter assays, eliciting 50-100 fold dynamic ranges in response to stimuli. Our platform is effective, easily implemented, and provides a solution for selecting short-length promoters with precise performance for a multitude of applications.

Structural and functional properties of mSWI/SNF chromatin remodeling complexes revealed through single-cell perturbation screens

The mammalian SWI/SNF (mSWI/SNF or BAF) family of chromatin remodeling complexes play critical roles in regulating DNA accessibility and gene expression. The three final-form subcomplexes-cBAF, PBAF, and ncBAF-are distinct in biochemical componentry, chromatin targeting, and roles in disease; however, the contributions of their constituent subunits to gene expression remain incompletely defined. Here, we performed Perturb-seq-based CRISPR-Cas9 knockout screens targeting mSWI/SNF subunits individually and in select combinations, followed by single-cell RNA-seq and SHARE-seq. We uncovered complex-, module-, and subunit-specific contributions to distinct regulatory networks and defined paralog subunit relationships and shifted subcomplex functions upon perturbations. Synergistic, intra-complex genetic interactions between subunits reveal functional redundancy and modularity. Importantly, single-cell subunit perturbation signatures mapped across bulk primary human tumor expression profiles both mirror and predict cBAF loss-of-function status in cancer. Our findings highlight the utility of Perturb-seq to dissect disease-relevant gene regulatory impacts of heterogeneous, multi-component master regulatory complexes.

A transcription factor atlas of directed differentiation

Transcription factors (TFs) regulate gene programs, thereby controlling diverse cellular processes and cell states. To comprehensively understand TFs and the programs they control, we created a barcoded library of all annotated human TF splice isoforms (>3,500) and applied it to build a TF Atlas charting expression profiles of human embryonic stem cells (hESCs) overexpressing each TF at single-cell resolution. We mapped TF-induced expression profiles to reference cell types and validated candidate TFs for generation of diverse cell types, spanning all three germ layers and trophoblasts. Targeted screens with subsets of the library allowed us to create a tailored cellular disease model and integrate mRNA expression and chromatin accessibility data to identify downstream regulators. Finally, we characterized the effects of combinatorial TF overexpression by developing and validating a strategy for predicting combinations of TFs that produce target expression profiles matching reference cell types to accelerate cellular engineering efforts.

Rate thresholds in cell signaling have functional and phenotypic consequences in non-linear time-dependent environments

All cells employ signal transduction pathways to respond to physiologically relevant extracellular cytokines, stressors, nutrient levels, hormones, morphogens, and other stimuli that vary in concentration and rate in healthy and diseased states. A central unsolved fundamental question in cell signaling is whether and how cells sense and integrate information conveyed by changes in the rate of extracellular stimuli concentrations, in addition to the absolute difference in concentration. We propose that different environmental changes over time influence cell behavior in addition to different signaling molecules or different genetic backgrounds. However, most current biomedical research focuses on acute environmental changes and does not consider how cells respond to environments that change slowly over time. As an example of such environmental change, we review cell sensitivity to environmental rate changes, including the novel mechanism of rate threshold. A rate threshold is defined as a threshold in the rate of change in the environment in which a rate value below the threshold does not activate signaling and a rate value above the threshold leads to signal activation. We reviewed p38/Hog1 osmotic stress signaling in yeast, chemotaxis and stress response in bacteria, cyclic adenosine monophosphate signaling in Amoebae, growth factors signaling in mammalian cells, morphogen dynamics during development, temporal dynamics of glucose and insulin signaling, and spatio-temproral stressors in the kidney. These reviewed examples from the literature indicate that rate thresholds are widespread and an underappreciated fundamental property of cell signaling. Finally, by studying cells in non-linear environments, we outline future directions to understand cell physiology better in normal and pathophysiological conditions.

Yeast osmoregulation - glycerol still in pole position

In response to osmotic dehydration cells sense, signal, alter gene expression, and metabolically counterbalance osmotic differences. The main compatible solute/osmolyte that accumulates in yeast cells is glycerol, which is produced from the glycolytic intermediate dihydroxyacetone phosphate. This review covers recent advancements in understanding mechanisms involved in sensing, signaling, cell-cycle delays, transcriptional responses as well as post-translational modifications on key proteins in osmoregulation. The protein kinase Hog1 is a key-player in many of these events, however, there is also a growing body of evidence for important Hog1-independent mechanisms playing vital roles. Several missing links in our understanding of osmoregulation will be discussed and future avenues for research proposed. The review highlights that this rather simple experimental system—salt/sorbitol and yeast—has developed into an enormously potent model system unravelling important fundamental aspects in biology.

Ddp1 Cooperates with Ppx1 to Counter a Stress Response Initiated by Nonvacuolar Polyphosphate

In diverse cells from bacterial to mammalian species, inorganic phosphate is stored in long chains called polyphosphate (polyP). These nearly universal polymers, ranging from three to thousands of phosphate moieties in length, are associated with molecular functions, including energy homeostasis, protein folding, and cell signaling. In many cell types, polyphosphate is concentrated in subcellular compartments or organelles. In the budding yeast Saccharomyces cerevisiae, polyP synthesis by the membrane-bound vacuolar transporter chaperone (VTC) complex is coupled to its translocation into the lumen of the vacuole, a lysosome-like organelle, where it is stored at high concentrations. In contrast, the ectopic expression of the bacterial polyphosphate kinase (PPK) results in the toxic accumulation of polyP outside the vacuole. In this study, we used label-free mass spectrometry to investigate the mechanisms underlying this toxicity. We find that PPK expression results in the activation of a stress response mediated in part by the Hog1 and Yak1 kinases and the Msn2/Msn4 transcription factors as well as by changes in protein kinase A (PKA) activity. This response is countered by the combined action of the Ddp1 and Ppx1 polyphosphatases that function together to counter polyP accumulation and downstream toxicity. In contrast, the ectopic expression of previously proposed mammalian polyphosphatases did not impact PPK-mediated toxicity in this model, suggesting either that these enzymes do not function directly as polyphosphatases in vivo or that they require cofactors unique to higher eukaryotes. Our work provides insight into why polyP accumulation outside lysosome-like organelles is toxic. Furthermore, it serves as a resource for exploring how polyP may impact conserved biological processes at a molecular level.

Inferring the contribution of small RNAs to changes in gene expression in response to stress

A main strategy of bacteria adapting to environmental changes is the remodeling of their transcriptome. Changes in the transcript levels of specific genes are due to combined effects of various regulators, including small RNAs (sRNAs). sRNAs are post-transcriptional regulators of gene expression that mainly control translation, but also directly and indirectly affect the levels of their target transcripts. Yet, the relative contribution of an sRNA to the total change in the transcript level of a gene upon an environmental change has not been assessed. We present a design of differential gene expression analysis by RNA-seq that allows extracting the contribution of an sRNA to the total change in the transcript level of each gene in response to an environmental change by fitting a linear model to the data. We exemplify this for the sRNA RyhB in cells growing under iron limitation and show a variation among genes in the relative contribution of RyhB to the change in their transcript level upon iron limitation, from subtle to very substantial. Extracting the relative contribution of an sRNA to the total change in expression of genes is important for understanding the integration of regulation by sRNAs with other regulatory mechanisms in the cell.

Degree of Freedom of Gene Expression in Saccharomyces cerevisiae

The complexity of genome-wide gene expression has not yet been adequately addressed due to a lack of comprehensive statistical analyses. In the present study, we introduce degree of freedom (DOF) as a summary statistic for evaluating gene expression complexity. Because DOF can be interpreted by a state-space representation, application of the DOF is highly useful for understanding gene activities. We used over 11,000 gene expression data sets to reveal that the DOF of gene expression in Saccharomyces cerevisiae is not greater than 450. We further demonstrated that various degrees of freedom of gene expression can be interpreted by different sequence motifs within promoter regions and Gene Ontology (GO) terms. The well-known TATA box is the most significant one among the identified motifs, while the GO term "ribosome genesis" is an associated biological process. On the basis of transcriptional freedom, our findings suggest that the regulation of gene expression can be modeled using only a few state variables. IMPORTANCE Yeast works like a well-organized factory. Each of its components works in its own way, while affecting the activities of others. The order of all activities is largely governed by the regulation of gene expression. In recent decades, biologists have recognized many regulations for yeast genes. However, it is not known how closely the regulation links each gene together to make all components of the cell work as a whole. In other words, biologists are very interested in how many independent control factors are needed to operate an artificial "cell" that works the same as a real one. In this work, we suggested that only 450 control factors were sufficient to represent the regulation of all 5800 yeast genes.

Data-driven identification of inherent features of eukaryotic stress-responsive genes

Living organisms are continuously challenged by changes in their environment that can propagate to stresses at the cellular level, such as rapid changes in osmolarity or oxygen tension. To survive these sudden changes, cells have developed stress-responsive mechanisms that tune cellular processes. The response of Saccharomyces cerevisiae to osmostress includes a massive reprogramming of gene expression. Identifying the inherent features of stress-responsive genes is of significant interest for understanding the basic principles underlying the rewiring of gene expression upon stress. Here, we generated a comprehensive catalog of osmostress-responsive genes from 5 independent RNA-seq experiments. We explored 30 features of yeast genes and found that 25 (83%) were distinct in osmostress-responsive genes. We then identified 13 non-redundant minimal osmostress gene traits and used statistical modeling to rank the most stress-predictive features. Intriguingly, the most relevant features of osmostress-responsive genes are the number of transcription factors targeting them and gene conservation. Using data on HeLa samples, we showed that the same features that define yeast osmostress-responsive genes can predict osmostress-responsive genes in humans, but with changes in the rank-ordering of feature-importance. Our study provides a holistic understanding of the basic principles of the regulation of stress-responsive gene expression across eukaryotes.

Involvement of the High-Osmolarity Glycerol Pathway of Saccharomyces Cerevisiae in Protection against Copper Toxicity

Although essential for life, copper is also potentially toxic in concentrations that surpass physiological thresholds. The high-osmolarity glycerol pathway of yeast is the main regulator of adaptive responses and is known to play crucial roles in the responses to various stressors. The objective of this research is to determine whether the HOG pathway could be activated and to investigate the possible interplay of the HOG pathway and oxidative stress due to copper exposure. In this research, we demonstrate that copper could induce oxidative stress, including the elevated concentrations of reactive oxygen species (ROS) and malondialdehyde (MDA). Increased combination with GSH, increased intracellular SOD activity, and the up-regulation of relevant genes can help cells defend themselves against oxidative toxicity. The results show that copper treatment triggers marked and prolonged Hog1 phosphorylation. Significantly, oxidative stress generated by copper toxicity is essential for the activation of Hog1. Activated Hog1 is translocated to the nucleus to regulate the expressions of genes such as CTT1, GPD1, and HSP12, among others. Furthermore, copper exposure induced significant G1-phase cell cycle arrest, while Hog1 partially participated in the regulation of cell cycle progression. These novel findings reveal another role for Hog1 in the regulation of copper-induced cellular stress.

Gene by Environment Interactions reveal new regulatory aspects of signaling network plasticity

Phenotypes can change during exposure to different environments through the regulation of signaling pathways that operate in integrated networks. How signaling networks produce different phenotypes in different settings is not fully understood. Here, Gene by Environment Interactions (GEIs) were used to explore the regulatory network that controls filamentous/invasive growth in the yeast Saccharomyces cerevisiae. GEI analysis revealed that the regulation of invasive growth is decentralized and varies extensively across environments. Different regulatory pathways were critical or dispensable depending on the environment, microenvironment, or time point tested, and the pathway that made the strongest contribution changed depending on the environment. Some regulators even showed conditional role reversals. Ranking pathways' roles across environments revealed an under-appreciated pathway (OPI1) as the single strongest regulator among the major pathways tested (RAS, RIM101, and MAPK). One mechanism that may explain the high degree of regulatory plasticity observed was conditional pathway interactions, such as conditional redundancy and conditional cross-pathway regulation. Another mechanism was that different pathways conditionally and differentially regulated gene expression, such as target genes that control separate cell adhesion mechanisms (FLO11 and SFG1). An exception to decentralized regulation of invasive growth was that morphogenetic changes (cell elongation and budding pattern) were primarily regulated by one pathway (MAPK). GEI analysis also uncovered a round-cell invasion phenotype. Our work suggests that GEI analysis is a simple and powerful approach to define the regulatory basis of complex phenotypes and may be applicable to many systems.

The high osmolarity glycerol (HOG) pathway in fungi†

While fungi are widely occupying nature, many species are responsible for devastating mycosis in humans. Such niche diversity explains how quick fungal adaptation is necessary to endow the capacity of withstanding fluctuating environments and to cope with host-imposed conditions. Among all the molecular mechanisms evolved by fungi, the most studied one is the activation of the phosphorelay signalling pathways, of which the high osmolarity glycerol (HOG) pathway constitutes one of the key molecular apparatus underpinning fungal adaptation and virulence. In this review, we summarize the seminal knowledge of the HOG pathway with its more recent developments. We specifically described the HOG-mediated stress adaptation, with a particular focus on osmotic and oxidative stress, and point out some lags in our understanding of its involvement in the virulence of pathogenic species including, the medically important fungi Candida albicans, Cryptococcus neoformans, and Aspergillus fumigatus, compared to the model yeast Saccharomyces cerevisiae. Finally, we also highlighted some possible applications of the HOG pathway modifications to improve the fungal-based production of natural products in the industry.

Engineered yeast tolerance enables efficient production from toxified lignocellulosic feedstocks

Lignocellulosic biomass remains unharnessed for the production of renewable fuels and chemicals due to challenges in deconstruction and the toxicity its hydrolysates pose to fermentation microorganisms. Here, we show in Saccharomyces cerevisiae that engineered aldehyde reduction and elevated extracellular potassium and pH are sufficient to enable near-parity production between inhibitor-laden and inhibitor-free feedstocks. By specifically targeting the universal hydrolysate inhibitors, a single strain is enhanced to tolerate a broad diversity of highly toxified genuine feedstocks and consistently achieve industrial-scale titers (cellulosic ethanol of >100 grams per liter when toxified). Furthermore, a functionally orthogonal, lightweight design enables seamless transferability to existing metabolically engineered chassis strains: We endow full, multifeedstock tolerance on a xylose-consuming strain and one producing the biodegradable plastics precursor lactic acid. The demonstration of "drop-in" hydrolysate competence enables the potential of cost-effective, at-scale biomass utilization for cellulosic fuel and nonfuel products alike.

Loss of Smi1, a protein involved in cell wall synthesis, extends replicative life span by enhancing rDNA stability in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, replicative life span (RLS) is primarily affected by the stability of ribosomal DNA (rDNA). The stability of the highly repetitive rDNA array is maintained through transcriptional silencing by the NAD⁺-dependent histone deacetylase Sir2. Recently, the loss of Smi1, a protein of unknown molecular function that has been proposed to be involved in cell wall synthesis, has been demonstrated to extend RLS in S. cerevisiae, but the mechanism by which Smi1 regulates RLS has not been elucidated. In this study, we determined that the loss of Smi1 extends RLS in a Sir2-dependent manner. We observed that the smi1Δ mutation enhances transcriptional silencing at the rDNA locus and promotes rDNA stability. In the absence of Smi1, the stress-responsive transcription factor Msn2 translocates from the cytoplasm to the nucleus, and nuclear-accumulated Msn2 stimulates the expression of nicotinamidase Pnc1, which serves as an activator of Sir2. In addition, we observed that the MAP kinase Hog1 is activated in smi1Δ cells and that the activation of Hog1 induces the translocation of Msn2 into the nucleus. Taken together, our findings suggest that the loss of Smi1 leads to the nuclear accumulation of Msn2 and stimulates the expression of Pnc1, thereby enhancing Sir2-mediated rDNA stability and extending RLS in S. cerevisiae.

Oxidative Stress Responses and Nutrient Starvation in MCHM Treated Saccharomyces cerevisiae

In 2014, the coal cleaning chemical 4-methylcyclohexane methanol (MCHM) spilled into the water supply for 300,000 West Virginians. Initial toxicology tests showed relatively mild results, but the underlying effects on cellular biology were underexplored. Treated wildtype yeast cells grew poorly, but there was only a small decrease in cell viability. Cell cycle analysis revealed an absence of cells in S phase within thirty minutes of treatment. Cells accumulated in G1 over a six-hour time course, indicating arrest instead of death. A genetic screen of the haploid knockout collection revealed 329 high confidence genes required for optimal growth in MCHM. These genes encode three major cell processes: mitochondrial gene expression/translation, the vacuolar ATPase, and aromatic amino acid biosynthesis. The transcriptome showed an upregulation of pleiotropic drug response genes and amino acid biosynthetic genes and downregulation in ribosome biosynthesis. Analysis of these datasets pointed to environmental stress response activation upon treatment. Overlap in datasets included the aromatic amino acid genes ARO1, ARO3, and four of the five TRP genes. This implicated nutrient deprivation as the signal for stress response. Excess supplementation of nutrients and amino acids did not improve growth on MCHM, so the source of nutrient deprivation signal is still unclear. Reactive oxygen species and DNA damage were directly detected with MCHM treatment, but timepoints showed these accumulated slower than cells arrested. We propose that wildtype cells arrest from nutrient deprivation and survive, accumulating oxidative damage through the implementation of robust environmental stress responses.

Heat preadaptation improved the ability of Zygosaccharomyces rouxii to salt stress: a combined physiological and transcriptomic analysis

Zygosaccharomyces rouxii plays important roles in the brewing process of fermented foods such as soy sauce, where salt stress is a frequently encountered condition. In this study, effect of heat preadaptation on salt tolerance of Z. rouxii and the protective mechanisms underlying heat preadaptation were investigated based on physiological and transcriptomic analyses. Results showed that cells subjected to heat preadaptation (37 °C, 90 min) prior to salt stress aroused many physiological responses, including maintaining cell surface smooth and intracellular pH level, increasing Na⁺/K⁺-ATPase activity. Cells subjected to heat preadaptation increased the amounts of unsaturated fatty acids (palmitoleic C16:1, oleic C18:1, linoleic C18:2) and decreased the amounts of saturated fatty acids (palmitic C16:0, stearic C18:0) which caused the unsaturation degree (unsaturated/saturated = U/S ratio) increased by 2.4 times when compared with cells without preadaptation under salt stress. Besides, salt stress led to increase in contents of 5 amino acids (valine, proline, threonine, glycine, and tyrosine) and decrease of 2 amino acids (serine and lysine). When comparing the cells pre-exposed to heat preadaptation followed by challenged with salt stress and the cells without preadaptation under salt stress, the serine, threonine, and lysine contents increased significantly. RNA sequencing revealed that the metabolic level of glycolysis by Z. rouxii was weakened, while the metabolic levels of the pentose phosphate pathway and the riboflavin were enhanced in cells during heat preadaptation. Results presented in this study may contribute to understand the bases of adaptive responses in Z. rouxii and rationalize its exploitation in industrial processes.Key points• Heat preadaptation can improve high salinity tolerance of Z. rouxii.• Combined physiological and transcriptomic analyses of heat preadaptation mechanisms.• Provide theoretical support for the application of Z. rouxii.

Transcription factors and ABC transporters: from pleiotropic drug resistance to cellular signaling in yeast

Budding yeast Saccharomyces cerevisiae survives in microenvironments utilizing networks of regulators and ATP-binding cassette (ABC) transporters to circumvent toxins and a variety of drugs. Our understanding of transcriptional regulation of ABC transporters in yeast is mainly derived from the study of multidrug resistance protein networks. Over the past two decades, this research has not only expanded the role of transcriptional regulators in pleiotropic drug resistance (PDR) but evolved to include the role that regulators play in cellular signaling and environmental adaptation. Inspection of the gene networks of the transcriptional regulators and characterization of the ABC transporters has clarified that they also contribute to environmental adaptation by controlling plasma membrane composition, toxic-metal sequestration, and oxidative stress adaptation. Additionally, ABC transporters and their regulators appear to be involved in cellular signaling for adaptation of S. cerevisiae populations to nutrient availability. In this review, we summarize the current understanding of the S. cerevisiae transcriptional regulatory networks and highlight recent work in other notable fungal organisms, underlining the expansion of the study of these gene networks across the kingdom fungi.

Yeast Ppz1 protein phosphatase toxicity involves the alteration of multiple cellular targets

Control of the protein phosphorylation status is a major mechanism for regulation of cellular processes, and its alteration often lead to functional disorders. Ppz1, a protein phosphatase only found in fungi, is the most toxic protein when overexpressed in Saccharomyces cerevisiae. To investigate the molecular basis of this phenomenon, we carried out combined genome-wide transcriptomic and phosphoproteomic analyses. We have found that Ppz1 overexpression causes major changes in gene expression, affecting ~ 20% of the genome, together with oxidative stress and increase in total adenylate pools. Concurrently, we observe changes in the phosphorylation pattern of near 400 proteins (mainly dephosphorylated), including many proteins involved in mitotic cell cycle and bud emergence, rapid dephosphorylation of Snf1 and its downstream transcription factor Mig1, and phosphorylation of Hog1 and its downstream transcription factor Sko1. Deletion of HOG1 attenuates the growth defect of Ppz1-overexpressing cells, while that of SKO1 aggravates it. Our results demonstrate that Ppz1 overexpression has a widespread impact in the yeast cells and reveals new aspects of the regulation of the cell cycle.

Hog1-induced transcription of RTC3 and HSP12 is robust and occurs in cells lacking Msn2, Msn4, Hot1 and Sko1

The yeast MAP kinase Hog1 pathway activates transcription of several hundreds genes. Large-scale gene expression and DNA binding assays suggest that most Hog1-induced genes are regulated by the transcriptional activators Msn2/4, Hot1 and Sko1. These studies also revealed the target genes of each activator and the putative binding sites on their promoters. In a previous study we identified a group of genes, which we considered the bona fide targets of Hog1, because they were induced in response to expression of intrinsically active mutant of Hog1, in the absence of any stress. We previously analyzed the promoter of the most highly induced gene, STL1, and noticed that some promoter properties were different from those proposed by large-scale data. We therefore continue to study promoters individually and present here analyses of promoters of more Hog1's targets, RTC3, HSP12, DAK1 and ALD3. We report that RTC3 and HSP12 promoters are robust and are induced, to different degrees, even in cells lacking all four activators. DAK1 and ALD3 promoters are not robust and fully depend on a single activator, DAK1 on Sko1 and ALD3 on Msn2/4. Most of these observations could not be inferred from the large-scale data. Msn2/4 are involved in regulating all four promoters. It was assumed, therefore, that the promoters are spontaneously active in ras2Δ cells, in which Msn2/4 are known to be de-repressed. Intriguingly, the promoters were not active in BY4741ras2Δ cells, but were de-repressed, as expected, in ras2Δ cells of other genetic backgrounds. This study describes two phenomena. One, some Hog1's target promoters are most robust, backupped by many activators. Second, in contrast to most laboratory strains, the widely used BY4741 strain does not induce Msn2/4 activity when the Ras/cAMP cascade is downregulated.

Quantitative analysis of yeast MAPK signaling networks and crosstalk using a microfluidic device

Eukaryotic cells developed complex mitogen-activated protein kinase (MAPK) signaling networks to sense their intra- and extracellular environment and respond to various stress conditions. For example, S. cerevisiae uses five distinct MAP kinase pathways to orchestrate meiosis or respond to mating pheromones, osmolarity changes and cell wall stress. Although each MAPK module has been studied individually, the mechanisms underlying crosstalk between signaling pathways remain poorly understood, in part because suitable experimental systems to monitor cellular outputs when applying different signals are lacking. Here, we investigate the yeast MAPK signaling pathways and their crosstalk, taking advantage of a new microfluidic device coupled to quantitative microscopy. We designed specific micropads to trap yeast cells in a single focal plane, and modulate the magnitude of a given stress signal by microfluidic serial dilution while keeping other signaling inputs constant. This approach enabled us to quantify in single cells nuclear relocation of effectors responding to MAPK activation, like Yap1 for oxidative stress, and expression of stress-specific reporter expression, like pSTL1-qV and pFIG1-qV for high-osmolarity or mating pheromone signaling, respectively. Using this quantitative single-cell analysis, we confirmed bimodal behavior of gene expression in response to Hog1 activation, and quantified crosstalk between the pheromone- and cell wall integrity (CWI) signaling pathways. Importantly, we further observed that oxidative stress inhibits pheromone signaling. Mechanistically, this crosstalk is mediated by Pkc1-dependent phosphorylation of the scaffold protein Ste5 on serine 185, which prevents Ste5 recruitment to the plasma membrane.

An expanded cell wall damage signaling network is comprised of the transcription factors Rlm1 and Sko1 in Candida albicans

The human fungal pathogen Candida albicans is constantly exposed to environmental challenges impacting the cell wall. Signaling pathways coordinate stress adaptation and are essential for commensalism and virulence. The transcription factors Sko1, Cas5, and Rlm1 control the response to cell wall stress caused by the antifungal drug caspofungin. Here, we expand the Sko1 and Rlm1 transcriptional circuit and demonstrate that Rlm1 activates Sko1 cell wall stress signaling. Caspofungin-induced transcription of SKO1 and several Sko1-dependent cell wall integrity genes are attenuated in an rlm1Δ/Δ mutant strain when compared to the treated wild-type strain but not in a cas5Δ/Δ mutant strain. Genome-wide chromatin immunoprecipitation (ChIP-seq) results revealed numerous Sko1 and Rlm1 directly bound target genes in the presence of caspofungin that were undetected in previous gene expression studies. Notable targets include genes involved in cell wall integrity, osmolarity, and cellular aggregation, as well as several uncharacterized genes. Interestingly, we found that Rlm1 does not bind to the upstream intergenic region of SKO1 in the presence of caspofungin, indicating that Rlm1 indirectly controls caspofungin-induced SKO1 transcription. In addition, we discovered that caspofungin-induced SKO1 transcription occurs through self-activation. Based on our ChIP-seq data, we also discovered an Rlm1 consensus motif unique to C. albicans. For Sko1, we found a consensus motif similar to the known Sko1 motif for Saccharomyces cerevisiae. Growth assays showed that SKO1 overexpression suppressed caspofungin hypersensitivity in an rlm1Δ/Δ mutant strain. In addition, overexpression of the glycerol phosphatase, RHR2, suppressed caspofungin hypersensitivity specifically in a sko1Δ/Δ mutant strain. Our findings link the Sko1 and Rlm1 signaling pathways, identify new biological roles for Sko1 and Rlm1, and highlight the complex dynamics underlying cell wall signaling.

Candida albicans is the most common human fungal pathogen isolated in clinical settings. The echinocandin drug caspofungin is used to treat invasive candidiasis; however, the emergence of increasing echinocandin resistance underscores the need for new antifungal strategies. Elucidating the signaling mechanisms that govern caspofungin-induced tolerance has the potential to identify candidate proteins that could serve as novel therapeutic targets. Here, we expand the Rlm1 and Sko1 cell wall transcriptional network and find that Rlm1 indirectly regulates Sko1 signaling. Furthermore, we identify Sko1- and Rlm1-specific biological roles in caspofungin adaptation, such as osmoregulation and secretion. Lastly, we discover a protective role for glycerol in caspofungin tolerance. Overall, these findings provide mechanistic insight into the genetic and cellular bases underlying cell wall signaling in C. albicans.

Single-particle imaging of stress-promoters induction reveals the interplay between MAPK signaling, chromatin and transcription factors

Precise regulation of gene expression in response to environmental changes is crucial for cell survival, adaptation and proliferation. In eukaryotic cells, extracellular signal integration is often carried out by Mitogen-Activated Protein Kinases (MAPK). Despite a robust MAPK signaling activity, downstream gene expression can display a great variability between single cells. Using a live mRNA reporter, here we monitor the dynamics of transcription in Saccharomyces cerevisiae upon hyper-osmotic shock. We find that the transient activity of the MAPK Hog1 opens a temporal window where stress-response genes can be activated. We show that the first minutes of Hog1 activity are essential to control the activation of a promoter. Chromatin repression on a locus slows down this transition and contributes to the variability in gene expression, while binding of transcription factors increases the level of transcription. However, soon after Hog1 activity peaks, negative regulators promote chromatin closure of the locus and transcription progressively stops.

The coordinate actions of calcineurin and Hog1 mediate the stress response through multiple nodes of the cell cycle network

Upon exposure to environmental stressors, cells transiently arrest the cell cycle while they adapt and restore homeostasis. A challenge for all cells is to distinguish between stress signals and coordinate the appropriate adaptive response with cell cycle arrest. Here we investigate the role of the phosphatase calcineurin (CN) in the stress response and demonstrate that CN activates the Hog1/p38 pathway in both yeast and human cells. In yeast, the MAPK Hog1 is transiently activated in response to several well-studied osmostressors. We show that when a stressor simultaneously activates CN and Hog1, CN disrupts Hog1-stimulated negative feedback to prolong Hog1 activation and the period of cell cycle arrest. Regulation of Hog1 by CN also contributes to inactivation of multiple cell cycle-regulatory transcription factors (TFs) and the decreased expression of cell cycle-regulated genes. CN-dependent downregulation of G1/S genes is dependent upon Hog1 activation, whereas CN inactivates G2/M TFs through a combination of Hog1-dependent and -independent mechanisms. These findings demonstrate that CN and Hog1 act in a coordinated manner to inhibit multiple nodes of the cell cycle-regulatory network. Our results suggest that crosstalk between CN and stress-activated MAPKs helps cells tailor their adaptive responses to specific stressors.

In order to survive exposure to environmental stress, cells transiently arrest the cell division cycle while they adapt to the stress. Several kinases and phosphatases are known to control stress adaptation programs, but the extent to which these signaling pathways work together to tune the stress response is not well understood. This study investigates the role of the phosphatase calcineurin in the stress response and shows that calcineurin inhibits the cell cycle in part by stimulating the activity of the Hog1/p-38 stress-activated MAPK in both yeast and human cells. Crosstalk between stress response pathways may help cells mount specific responses to diverse stressors and to survive changes in their environment.

A genetic analysis reveals novel histone residues required for transcriptional reprogramming upon stress

Cells have the ability to sense, respond and adapt to environmental fluctuations. Stress causes a massive reorganization of the transcriptional program. Many examples of histone post-translational modifications (PTMs) have been associated with transcriptional activation or repression under steady-state growth conditions. Comparatively less is known about the role of histone PTMs in the cellular adaptive response to stress. Here, we performed high-throughput genetic screenings that provide a novel global map of the histone residues required for transcriptional reprogramming in response to heat and osmotic stress. Of note, we observed that the histone residues needed depend on the type of gene and/or stress, thereby suggesting a 'personalized', rather than general, subset of histone requirements for each chromatin context. In addition, we identified a number of new residues that unexpectedly serve to regulate transcription. As a proof of concept, we characterized the function of the histone residues H4-S47 and H4-T30 in response to osmotic and heat stress, respectively. Our results uncover novel roles for the kinases Cla4 and Ste20, yeast homologs of the mammalian PAK2 family, and the Ste11 MAPK as regulators of H4-S47 and H4-T30, respectively. This study provides new insights into the role of histone residues in transcriptional regulation under stress conditions.

Multi-kinase control of environmental stress responsive transcription

Cells respond to changes in environmental conditions by activating signal transduction pathways and gene expression programs. Here we present a dataset to explore the relationship between environmental stresses, kinases, and global gene expression in yeast. We subjected 28 drug-sensitive kinase mutants to 10 environmental conditions in the presence of inhibitor and performed mRNA deep sequencing. With these data, we reconstructed canonical stress pathways and identified examples of crosstalk among pathways. The data also implicated numerous kinases in novel environment-specific roles. However, rather than regulating dedicated sets of target genes, individual kinases tuned the magnitude of induction of the environmental stress response (ESR)–a gene expression signature shared across the set of perturbations–in environment-specific ways. This suggests that the ESR integrates inputs from multiple sensory kinases to modulate gene expression and growth control. As an example, we provide experimental evidence that the high osmolarity glycerol pathway is an upstream negative regulator of protein kinase A, a known inhibitor of the ESR. These results elaborate the central axis of cellular stress response signaling.

The InsP7 phosphatase Siw14 regulates inositol pyrophosphate levels to control localization of the general stress response transcription factor Msn2

The environmental stress response (ESR) is critical for cell survival. Yeast cells unable to synthesize inositol pyrophosphates (PP-InsPs) are unable to induce the ESR. We recently discovered a diphosphoinositol pentakisphosphate (PP-InsP₅) phosphatase in Saccharomyces cerevisiae encoded by SIW14 Yeast strains deleted for SIW14 have increased levels of PP-InsPs. We hypothesized that strains with high inositol pyrophosphate levels will have an increased stress response. We examined the response of the siw14Δ mutant to heat shock, nutrient limitation, osmotic stress, and oxidative treatment using cell growth assays and found increased resistance to each. Transcriptional responses to oxidative and osmotic stresses were assessed using microarray and reverse transcriptase quantitative PCR. The ESR was partially induced in the siw14Δ mutant strain, consistent with the increased stress resistance, and the mutant strain further induced the ESR in response to oxidative and osmotic stresses. The levels of PP-InsPs increased in WT cells under oxidative stress but not under hyperosmotic stress, and they were high and unchanging in the mutant. Phosphatase activity of Siw14 was inhibited by oxidation that was reversible. To determine how altered PP-InsP levels affect the ESR, we performed epistasis experiments with mutations in rpd3 and msn2/4 combined with siw14Δ. We show that mutations in msn2Δ and msn4Δ, but not rpd3, are epistatic to siw14Δ by assessing growth under oxidative stress conditions and expression of CTT1 Msn2-GFP nuclear localization was increased in the siw14Δ. These data support a model in which the modulation of PP-InsPs influence the ESR through general stress response transcription factors Msn2/4.

Overlapping responses between salt and oxidative stress in Debaryomyces hansenii

Debaryomyces hansenii is a halotolerant yeast of importance in basic and applied research. Previous reports hinted about possible links between saline and oxidative stress responses in this yeast. The aim of this work was to study that hypothesis at different molecular levels, investigating after oxidative and saline stress: (i) transcription of seven genes related to oxidative and/or saline responses, (ii) activity of two main anti-oxidative enzymes, (iii) existence of common metabolic intermediates, and (iv) generation of damages to biomolecules as lipids and proteins. Our results showed how expression of genes related to oxidative stress was induced by exposure to NaCl and KCl, and, vice versa, transcription of some genes related to osmotic/salt stress responses was regulated by H₂O₂. Moreover, and contrary to S. cerevisiae, in D. hansenii HOG1 and MSN2 genes were modulated by stress at their transcriptional level. At the enzymatic level, saline stress also induced antioxidative enzymatic defenses as catalase and glutathione reductase. Furthermore, we demonstrated that both stresses are connected by the generation of intracellular ROS, and that hydrogen peroxide can affect the accumulation of in-cell sodium. On the other hand, no significant alterations in lipid oxidation or total glutathione content were observed upon exposure to both stresses tested. The results described in this work could help to understand the responses to both stressors, and to improve the biotechnological potential of D. hansenni.

Nut1/Hos1 and Sas2/Rpd3 control the H3 acetylation of two different sets of osmotic stress-induced genes

Epigenetic information is able to interact with the cellular environment and could be especially useful for reprograming gene expression in response to a physiological perturbation. In fact the genes induced or repressed by osmotic stress undergo significant changes in terms of the levels of various histone modifications, especially in the acetylation levels of histone H3. Exposing yeast to high osmolarity results in the activation of stress-activated protein kinase Hog1, which plays a central role in gene expression control. We evaluated the connection between the presence of Hog1 and changes in histone H3 acetylation in stress-regulated genes. We found a parallel increase in the acetylation of lysines 9 and 14 of H3 in induced genes during stress, which was largely dependent on Hog1 at the genome-wide level. Conversely, we observed that acetylation decreased in repressed genes and was not dependent on Hog1. However, lack of Hog1 sometimes produced different, and even opposite, effects on the induction and acetylation of H3 of each gene. We also found that the acetylation state of lysine 9 of H3 was altered in the strains deficient in Nut1 HAT and Hos1 HDAC in the genes up-regulated during osmotic stress in an Msn2/Msn4-independent manner, while lysine 9 acetylation of H3 varied in the strains deficient in Sas2 HAT and Rpd3 HDAC for the Msn2/Msn4-dependent induced genes. The results presented here show new, unexpected participants in gene regulation processes in response to environmental perturbations.

Integrated TORC1 and PKA signaling control the temporal activation of glucose-induced gene expression in yeast

The growth rate of a yeast cell is controlled by the target of rapamycin kinase complex I (TORC1) and cAMP-dependent protein kinase (PKA) pathways. To determine how TORC1 and PKA cooperate to regulate cell growth, we performed temporal analysis of gene expression in yeast switched from a non-fermentable substrate, to glucose, in the presence and absence of TORC1 and PKA inhibitors. Quantitative analysis of these data reveals that PKA drives the expression of key cell growth genes during transitions into, and out of, the rapid growth state in glucose, while TORC1 is important for the steady-state expression of the same genes. This circuit design may enable yeast to set an exact growth rate based on the abundance of internal metabolites such as amino acids, via TORC1, but also adapt rapidly to changes in external nutrients, such as glucose, via PKA.

Genetic screen of the yeast environmental stress response dynamics uncovers distinct regulatory phases

Cells respond to environmental fluctuations by regulating multiple transcriptional programs. This response can be studied by measuring the effect of environmental changes on the transcriptome or the proteome of the cell at the end of the response. However, the dynamics of the response reflect the working of the regulatory mechanisms in action. Here, we utilized a fluorescent stress reporter gene to track the dynamics of protein production in yeast responding to environmental stress. The response is modulated by changes in both the duration and rate of transcription. We probed the underlying molecular pathways controlling these two dimensions using a library of ~1,600 single- and double-mutant strains. Dissection of the effects of these mutants and the interactions between them identified distinct modulators of response duration and response rate. Using a combination of mRNA-seq and live-cell microscopy, we uncover mechanisms by which Msn2/4, Mck1, Msn5, and the cAMP/PKA pathway modulate the response of a large module of stress-induced genes in two discrete regulatory phases. Our results and analysis show that transcriptional stress response is regulated by multiple mechanisms that overlap in time and cellular location.

Unsupervised GRN Ensemble

Inferring gene regulatory networks from expression data is a very challenging problem that has raised the interest of the scientific community. Different algorithms have been proposed to try to solve this issue, but it has been shown that different methods have some particular biases and strengths, and none of them is the best across all types of data and datasets. As a result, the idea of aggregating various network inferences through a consensus mechanism naturally arises. In this chapter, a common framework to standardize already proposed consensus methods is presented, and based on this framework different proposals are introduced and analyzed in two different scenarios: Homogeneous and Heterogeneous. The first scenario reflects situations where the networks to be aggregated are rather similar because they are obtained with inference algorithms working on the same data, whereas the second scenario deals with very diverse networks because various sources of data are used to generate the individual networks. A procedure for combining multiple network inference algorithms is analyzed in a systematic way. The results show that there is a very significant difference between these two scenarios, and that the best way to combine networks in the Heterogeneous scenario is not the most commonly used. We show in particular that aggregation in the Heterogeneous scenario can be very beneficial if the individual networks are combined with our new proposed method ScaleLSum.

PKA and HOG signaling contribute separable roles to anaerobic xylose fermentation in yeast engineered for biofuel production

Lignocellulosic biomass offers a sustainable source for biofuel production that does not compete with food-based cropping systems. Importantly, two critical bottlenecks prevent economic adoption: many industrially relevant microorganisms cannot ferment pentose sugars prevalent in lignocellulosic medium, leaving a significant amount of carbon unutilized. Furthermore, chemical biomass pretreatment required to release fermentable sugars generates a variety of toxins, which inhibit microbial growth and metabolism, specifically limiting pentose utilization in engineered strains. Here we dissected genetic determinants of anaerobic xylose fermentation and stress tolerance in chemically pretreated corn stover biomass, called hydrolysate. We previously revealed that loss-of-function mutations in the stress-responsive MAP kinase HOG1 and negative regulator of the RAS/Protein Kinase A (PKA) pathway, IRA2, enhances anaerobic xylose fermentation. However, these mutations likely reduce cells' ability to tolerate the toxins present in lignocellulosic hydrolysate, making the strain especially vulnerable to it. We tested the contributions of Hog1 and PKA signaling via IRA2 or PKA negative regulatory subunit BCY1 to metabolism, growth, and stress tolerance in corn stover hydrolysate and laboratory medium with mixed sugars. We found mutations causing upregulated PKA activity increase growth rate and glucose consumption in various media but do not have a specific impact on xylose fermentation. In contrast, mutation of HOG1 specifically increased xylose usage. We hypothesized improving stress tolerance would enhance the rate of xylose consumption in hydrolysate. Surprisingly, increasing stress tolerance did not augment xylose fermentation in lignocellulosic medium in this strain background, suggesting other mechanisms besides cellular stress limit this strain's ability for anaerobic xylose fermentation in hydrolysate.

A conserved role for transcription factor sumoylation in binding-site selection

Large numbers of eukaryotic transcription factors (TFs) are modified by SUMO post-translational modifications. Whereas the effect of TF sumoylation on the expression of target genes is largely context-dependent, it is not known whether the modification has a common function in regulating TFs throughout eukaryotic species. Here, I highlight four studies that used genome-wide chromatin-immunoprecipitation analysis (ChIP-seq) to examine whether sumoylation affects the selection of sites on the genome that are bound by human and yeast TFs. The studies found that impairing sumoylation led to deregulated binding-site selection for all four of the examined TFs. Predominantly, compared to wild-type forms, the sumoylation-deficient forms of the TFs bound to numerous additional non-specific sites, pointing to a common role for the modification in restricting TF binding to appropriate sites. Evidence from these studies suggests that TF sumoylation influences binding-site selection by modulating protein-protein interactions with other DNA-binding TFs, or by promoting conformational changes in the TFs that alter their DNA-binding specificity or affinity. I propose a model in which, prior to their sumoylation, TFs initially bind to chromatin with reduced specificity, which leads to spurious binding but also ensures that all functional sites become bound. Once the TFs are bound to DNA, sumoylation then acts to increase specificity and promotes release of the TFs from non-specific sites. The similar observations from these four genome-wide studies across divergent species suggest that binding-site selection is a general and conserved function for TF sumoylation.

The ability of transcription factors to differentially regulate gene expression is a crucial component of the mechanism underlying inversion, a frequently observed genetic interaction pattern

Genetic interactions, a phenomenon whereby combinations of mutations lead to unexpected effects, reflect how cellular processes are wired and play an important role in complex genetic diseases. Understanding the molecular basis of genetic interactions is crucial for deciphering pathway organization as well as understanding the relationship between genetic variation and disease. Several hypothetical molecular mechanisms have been linked to different genetic interaction types. However, differences in genetic interaction patterns and their underlying mechanisms have not yet been compared systematically between different functional gene classes. Here, differences in the occurrence and types of genetic interactions are compared for two classes, gene-specific transcription factors (GSTFs) and signaling genes (kinases and phosphatases). Genome-wide gene expression data for 63 single and double deletion mutants in baker's yeast reveals that the two most common genetic interaction patterns are buffering and inversion. Buffering is typically associated with redundancy and is well understood. In inversion, genes show opposite behavior in the double mutant compared to the corresponding single mutants. The underlying mechanism is poorly understood. Although both classes show buffering and inversion patterns, the prevalence of inversion is much stronger in GSTFs. To decipher potential mechanisms, a Petri Net modeling approach was employed, where genes are represented as nodes and relationships between genes as edges. This allowed over 9 million possible three and four node models to be exhaustively enumerated. The models show that a quantitative difference in interaction strength is a strict requirement for obtaining inversion. In addition, this difference is frequently accompanied with a second gene that shows buffering. Taken together, these results provide a mechanistic explanation for inversion. Furthermore, the ability of transcription factors to differentially regulate expression of their targets provides a likely explanation why inversion is more prevalent for GSTFs compared to kinases and phosphatases.

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.

Role of the highly conserved G68 residue in the yeast phosphorelay protein Ypd1: implications for interactions between histidine phosphotransfer (HPt) and response regulator proteins

Background

Many bacteria and certain eukaryotes utilize multi-step His-to-Asp phosphorelays for adaptive responses to their extracellular environments. Histidine phosphotransfer (HPt) proteins function as key components of these pathways. HPt proteins are genetically diverse, but share a common tertiary fold with conserved residues near the active site. A surface-exposed glycine at the H + 4 position relative to the phosphorylatable histidine is found in a significant number of annotated HPt protein sequences. Previous reports demonstrated that substitutions at this position result in diminished phosphotransfer activity between HPt proteins and their cognate signaling partners.

Results

We report the analysis of partner binding interactions and phosphotransfer activity of the prototypical HPt protein Ypd1 from Saccharomyces cerevisiae using a set of H + 4 (G68) substituted proteins. Substitutions at this position with large, hydrophobic, or charged amino acids nearly abolished phospho-acceptance from the receiver domain of its upstream signaling partner, Sln1 (Sln1-R1). An in vitro binding assay indicated that G68 substitutions caused only modest decreases in affinity between Ypd1 and Sln1-R1, and these differences did not appear to be large enough to account for the observed decrease in phosphotransfer activity. The crystal structure of one of these H + 4 mutants, Ypd1-G68Q, which exhibited a diminished ability to participate in phosphotransfer, shows a similar overall structure to that of wild-type. Molecular modelling suggests that the highly conserved active site residues within the receiver domain of Sln1 must undergo rearrangement to accommodate larger H + 4 substitutions in Ypd1.

Conclusions

Phosphotransfer reactions require precise arrangement of active site elements to align the donor-acceptor atoms and stabilize the transition state during the reaction. Any changes likely result in an inability to form a viable transition state during phosphotransfer. Our data suggest that the high degree of evolutionary conservation of residues with small side chains at the H + 4 position in HPt proteins is required for optimal activity and that the presence of larger residues at the H + 4 position would cause alterations in the positioning of active site residues in the partner response regulator.

Electronic supplementary material

The online version of this article (10.1186/s12858-019-0104-5) contains supplementary material, which is available to authorized users.

Molecular characterization of p38 MAPK from blunt snout bream (Megalobrama amblycephala) and its expression after ammonia stress, and lipopolysaccharide and bacterial challenge

p38 mitogen-activated protein kinase (MAPK) is an important protein which plays a key role in regulating the innate immunity, so exploring its molecular characterization is helpful in understanding the resistance against microbial infections in cultured fish. Here, a full-length cDNA of p38 MAPK was cloned from liver of blunt snout bream (Megalobrama amblycephala) which covered 2419 bp with an open reading frame of 1086 bp encoding 361 amino acids. p38 MAPK contained the characteristic structures of Thr-Gly-Tyr (TGY) motif and substrate binding site Ala-Thr-Arg-Trp (ATRW), which are conserved in MAPK family. To investigate p38 MAPK functions, two in vivo experiments were carried out to examine its expression following ammonia exposure and bacterial challenge. Also, an in vitro experiment was conducted to assess the role of p38 MAPK in inflammation of primary hepatocytes induced by lipopolysaccharide (LPS). The results showed the ubiquitous expression of p38 MAPK in all the tested tissues with varying levels. p38 MAPK mRNA expression was significantly up-regulated by ammonia stress and Aeromonas hydrophila challenge, and altered in a time-dependent manner. Moreover, the results indicated that the inflammatory response induced by LPS in hepatocytes is p38 MAPK dependent as knockdown of p38 MAPK using siRNA technology depressed the expression of IL-1β and IL-6. The findings in this study showed that p38 MAPK has anti-stress property, and plays key role in protection against bacterial infection and inflammation in blunt snout bream.

A constitutive active allele of the transcription factor Msn2 mimicking low PKA activity dictates metabolic remodeling in yeast

In yeast, protein kinase A (PKA) adjusts transcriptional profiles, metabolic rates, and cell growth in accord with carbon source availability. PKA affects gene expression mostly via the transcription factors Msn2 and Msn4, two key regulators of the environmental stress response. Here we analyze the role of the PKA-Msn2 signaling module using an Msn2 allele that harbors serine-to-alanine substitutions at six functionally important PKA motifs (Msn2A6) . Expression of Msn2A6 mimics low PKA activity, entails a transcription profile similar to that of respiring cells, and prevents formation of colonies on glucose-containing medium. Furthermore, Msn2A6 leads to high oxygen consumption and hence high respiratory activity. Substantially increased intracellular concentrations of several carbon metabolites, such as trehalose, point to a metabolic adjustment similar to diauxic shift. This partial metabolic switch is the likely cause for the slow-growth phenotype in the presence of glucose. Consistently, Msn2A6 expression does not interfere with growth on ethanol and tolerated is to a limited degree in deletion mutant strains with a gene expression signature corresponding to nonfermentative growth. We propose that the lethality observed in mutants with hampered PKA activity resides in metabolic reprogramming that is initiated by Msn2 hyperactivity.

Synthesizing Signaling Pathways from Temporal Phosphoproteomic Data

We present a method for automatically discovering signaling pathways from time-resolved phosphoproteomic data. The Temporal Pathway Synthesizer (TPS) algorithm uses constraint-solving techniques first developed in the context of formal verification to explore paths in an interaction network. It systematically eliminates all candidate structures for a signaling pathway where a protein is activated or inactivated before its upstream regulators. The algorithm can model more than one hundred thousand dynamic phosphosites and can discover pathway members that are not differentially phosphorylated. By analyzing temporal data, TPS defines signaling cascades without needing to experimentally perturb individual proteins. It recovers known pathways and proposes pathway connections when applied to the human epidermal growth factor and yeast osmotic stress responses. Independent kinase mutant studies validate predicted substrates in the TPS osmotic stress pathway.