CRISPRi screens reveal genes modulating yeast growth in lignocellulose hydrolysate
Biotechnology for Biofuels volume 14, Article number: 41 (2021)
Baker’s yeast is a widely used eukaryotic cell factory, producing a diverse range of compounds including biofuels and fine chemicals. The use of lignocellulose as feedstock offers the opportunity to run these processes in an environmentally sustainable way. However, the required hydrolysis pretreatment of lignocellulosic material releases toxic compounds that hamper yeast growth and consequently productivity.
Here, we employ CRISPR interference in S. cerevisiae to identify genes modulating fermentative growth in plant hydrolysate and in presence of lignocellulosic toxins. We find that at least one-third of hydrolysate-associated gene functions are explained by effects of known toxic compounds, such as the decreased growth of YAP1 or HAA1, or increased growth of DOT6 knock-down strains in hydrolysate.
Our study confirms previously known genetic elements and uncovers new targets towards designing more robust yeast strains for the utilization of lignocellulose hydrolysate as sustainable feedstock, and, more broadly, paves the way for applying CRISPRi screens to improve industrial fermentation processes.
The baker’s yeast S. cerevisiae is the most frequently used eukaryotic cell factory [9, 35, 53, 59]. The competitive production of biofuels and many other value-added compounds with yeast requires the use of cheap substrates that do not compete with food, feed and arable land. Lignocellulosic materials represent an economic and environmentally sustainable alternative feedstock. Spruce softwood is a promising lignocellulose source in the northern hemisphere [85, 87] and an abundant side product of the lumber, pulp and paper industry . Lignocellulose has a complex structure, largely consisting of cellulose, hemicellulose and lignin. The extraction of fermentable mono-saccharide sugars (glucose) and hemicellulose polymers (pentose and hexose sugars) from cellulose biomass requires hydrolysis pretreatment. During this processing step, toxic compounds are released to the soluble hydrolysate which represent a major challenge in using this feedstock in an industrial setting [2, 5, 64, 66, 82]. These compounds can be classified in three main groups: aliphatic acids, furan aldehydes and phenolic/aromatic derivatives . Aliphatic acids lower the intracellular pH and interfere with nutrient uptake . Furans inhibit dehydrogenases and raise levels of reactive oxygen species [1, 7]. Phenolics perturb plasma membrane composition and potential, resulting in a disruption of cell signalling and sorting processes . Understanding the impact of these substances on yeast growth and the cellular mechanisms of tolerance will enhance the use of lignocellulose-based biotechnological processes.
Several efforts have been made to characterize transcriptome changes  and to improve substrate utilization and tolerance to lignocellulosic inhibitors in industrial yeast strains [13, 73, 86], including deletion collection screens for synthetic and straw lignocellulosic hydrolysates which linked tolerance mechanisms to ATPase activity and pH, pentose phosphate metabolism, lipid metabolism and the biosynthesis of amino acids [65, 77]. The traditional generation of deletion collections for industrial strains is laborious, since it requires the change of genomic sequence in multiple alleles, and cannot assess the effects of transcript down-regulation. The emerging CRISPR-based knock-out, interference and activation systems thus offer genetic screens with broader phenotypic scope.
Here, we established CRISPR interference (CRISPRi) screens to identify genetic functions that tune yeast growth in spruce hydrolysate. CRISPRi is emerging as a powerful tool to study genotype–phenotype relations via precise transcriptional repression [25, 26, 55, 56, 78]. We employed an in-house developed single plasmid CRISPRi system  which is inducible by anhydrotetracycline (ATc)  to repress S. cerevisiae transcription factors (TFs, n = 161)  and protein kinases (PKs, n = 129) , key players in the regulation of cellular mechanisms and pathways. We detect genes capable of modulating yeast growth in hydrolysate, as well as in presence of toxic lignocellulose compounds to understand toxicity-dependent effects. We confirm previously known links to growth in hydrolysate and discover novel genetic associations that can be directly applied to advance sustainable bioprocesses (Fig. 1a).
Characterization of yeast fermentation in lignocellulosic hydrolysate and in presence of growth inhibitor compounds
In this study, we make use of CRISPR interference to profile genetic functions affecting growth of the diploid BY4743 yeast in spruce hydrolysate-containing media and in presence of growth-inhibitory compounds typically present in spruce and other hydrolysates (Fig. 1a). The laboratory BY4743 strain is well characterized, genetically amenable and comes with lots of available high-quality datasets on genome sequencing, annotation and function, in contrast to most process-specialized and frequently polyploid industrial yeast strains . We characterized growth of BY4743 across different media in order to identify conditions suitable for genetic screens, and further quantified glucose consumption and metabolite secretion profiles under the selected screening conditions.
To assess the general feasibility of employing BY4743 for fermentation, we were interested in the strains’ ability to grow the presence of high amounts of ethanol (EtOH). As expected, high EtOH concentrations of 5% or 10% decreased maximum growth rate by approximately 30% or 50%, respectively (Additional file 1: Fig. S1a). We then went on to characterize tolerance to lignocellulosic hydrolysate and to a cocktail of eight selected growth inhibitors (inhibitor cocktail, referred to as IC) which are commonly found in lignocellulosic hydrolysates, based on previous studies (Table 1 and Methods). Comparing growth in different dilutions of an industrially derived and widely used spruce hydrolysate [15, 36, 46, 47, 81, 86], we found that the BY4743 strain grew well in media supplemented with up to 14% hydrolysate without considerable changes in growth rate, while no growth was observed in hydrolysate concentrations of 20% and more (Additional file 1: Fig. S1b). Culturing yeast in 50-ml flasks, we measured decreased growth in the 10% hydrolysate dilution (Fig. 1b). We further measured growth in SCM supplemented with different concentrations of inhibitor cocktail. Growth profiles in 15–35% of inhibitor cocktail were similar to those in SCM in well-plate format (Additional file 1: Fig. S1c). For 50-ml flask cultures, we observed decreased growth in 45% IC with comparable impact to 10% hydrolysate (Fig. 1b) and therefore selected these conditions for further experiments, including the CRISPRi-based competition assays (Fig. 1a).
We next quantified the concentration of extracellular metabolites during yeast cultivation in synthetic complete media (SCM), 10% hydrolysate and 45% inhibitor compound cocktail (IC). Remarkably, ethanol yields were increased by ~ 65% in hydrolysate-supplemented medium compared to SCM and IC conditions (Fig. 1c). The increased ethanol production can partially be explained by the ~ 25% higher initial glucose levels (Additional file 2: Table S1) and the breakdown of large sugar polymers during hydrolysis pretreatment that presumably increased the concentration of fermentable carbohydrates , such as galactose and mannose , see Methods). We further quantified acetic acid, a common growth inhibitor, and found its concentration to decrease by approximately 50% during fermentation, while glucose was fully consumed (Additional file 3: Fig. S2). Interestingly, this suggests that acetic acid, and potentially also other growth-inhibitory compounds, can be metabolized under the given conditions. Taken together, the quantified dynamics in growth and metabolite abundances during fermentations in hydrolysate and in the presence of an lignocellulosic inhibitor mixture which showed that growth-inhibitory substances can be metabolized and enabled us to derive optimal conditions for functional genomics assays.
CRISPRi effects are reproducible and capture positive controls
To identify genes that modulate yeast growth in the presence of lignocellulose hydrolysate and inhibitor cocktail, we employed a plasmid-based CRISPRi system that allows for continuous expression of a nuclease-deactivated (d)Cas9 protein fused to the potent Mxi1 repressor domain , as well as the inducible expression of a gRNA from an anhydrotetracycline (ATc)-controlled promoter. For gene dosage screens, we used a library of 1573 gRNAs to repress 161 transcription factors and 129 protein kinases which are highly involved in the regulation of growth adaptation and cellular signalling. While the centromeric plasmids enable repression of only a single gene in each cell, we designed libraries with up to six gRNA locations per gene target .
Repression effects on cellular fitness were assayed in different media (SCM, SCM 10% hydrolysate and SCM 45% IC), performing selections of ~ 25 generations for all screens to enable direct comparison of effects on doubling time (Fig. 2a). Each condition was assayed in triplicate in CRISPRi-inducing (+ ATc) as well as in reference (-ATc) conditions. After selection in different media, plasmids were extracted and gRNA barcodes quantified by sequencing to compare strain abundances between + ATc and -ATc populations. Multidimensional scaling (MDS) of read count samples indicate high similarity between replicates and allow to estimate sources of variability of CRISPRi effect size (Fig. 2b). Samples of the inhibitor cocktail (SCM + 45% IC) were positioned closely to hydrolysate samples (SCM + 10% Hydrolysate), suggesting higher similarity to this condition than compared with SCM (Fig. 2b). For a detailed comparison of fitness effects, we provide correlations of read counts (Additional file 4: Fig. S3), guide RNA (Additional file 5: Fig. S4) and gene log2-fold changes (log2FCs) across all screen conditions (Additional file 6: Fig. S5). Read counts of the three biological replicates profiled for each condition are highly correlated, giving confidence in the derived effects (Additional file 4: Fig. S3). Notably, we further found gRNA fold changes highly correlated with those of a previous study where the TF and PK libraries have been profiled individually in SCM with different transformation batches and without oxygen limitation , indicating high reproducibility of CRISPRi effects (Pearson R2 = 0.67, p-value < 2.2e−16, Fig. 2c). Accordingly, 67% of variation between CRISPRi effects observed in one of the screens were explained by the other. Interestingly, we measured increased fitness for the repression of three genes (HAP1, RIM11 and RME1) that were not significant in the previous study , presumably due to the oxygen-limited conditions applied here to focus on fermentation. This is supported by the haeme-activated protein 1 (Hap1) TF which represses Rox1 in non-stress conditions . In response to hypoxia, Hap1 is inhibited to de-repress Rox1 and induce hypoxic stress gene expression [39, 49, 88] which would be enhanced in CRISPRi strains. The Rme1 TF and the Rim11 PK (one of four glutathione synthetase kinase 3 homologs) regulate meiosis, and their knock-down effects may hint to growth-antagonizing roles in oxygen-limited environments. As anticipated from a competitive fitness assay, genes essential for viability (based on the Saccharomyces Genome Database SGD) , showed higher depletion compared to non-essential genes and represent positive controls that validate the experimental setup (Fig. 2d).
Having validated gene dosage effects in SCM and their reliability across studies, we screened CRISPRi populations in 10% spruce hydrolysate (Fig. 2a). This revealed fitness effects (gene fold change with false discovery rate (FDR) < 0.05 and at least two gRNAs with absolute log2FC ≥ 1 and FDR < 0.05) for the repression of 24 genes (Fig. 3a and Additional file 7: Table S2). Ten of these also showed growth effects in SCM (brown dots in Fig. 3a). Four more genes also caused general fitness defects although missing the strict significance requirements in the SCM screen. This left ten genes with hydrolysate-associated roles and minor or no impact on growth (labelled in Fig. 3a). Repression of seven genes decreased (YAP1, HAA1, HOG1, PBS2, UGA3, CDC15 and UME6) and of three genes increased growth (DOT6, SKO1 and BUB1) specifically in hydrolysate. Most of these have well-known roles in adaptation to diverse kinds of stress.
Repression of YAP1 results in decreased growth in SCM + 10% hydrolysate medium (Fig. 3a). The Yap1 TF induces gene expression in response to oxidative stress [48, 60] and is known to be activated by toxic compounds including furans  and phenolic molecules . In line with CRISPRi effects, the yap1Δ mutant has decreased fitness in hydrolysate generated from Miscanthus plants , and a recent study also demonstrated that YAP1 overexpression increases growth in spruce hydrolysate , thus providing a direct validation of screen effects and illustrating how gene functions can be applied for strain design in biotechnology.
Likewise, we measured decreased fitness in hydrolysate for HAA1 in CRISPRi screening (Fig. 3a). The Haa1 TF is required for adaptation to mildly acidic environments , so that effects are most likely explained by the acidic pH of spruce hydrolysate (pH 4.5 of the 10% hydrolysate medium, compared to pH 5.5 of SCM). Akin to YAP1, the overexpression of HAA1 in an industrial yeast strain increased growth rate in hardwood hydrolysate and additionally improved ethanol production .
We further found reduced hydrolysate fitness upon repression of UME6 (Fig. 3a), which encodes a regulator of meiotic and translation-related genes . Accordingly, the ume6Δ mutant is sensitive to oxidative , hyperosmotic  high temperature stress  and to more than 20 diverse chemicals some of which may have properties similar to lignocellulosic compounds (based on SGD ). The ume6Δ strain has also been measured with decreased fitness in Miscanthus plant-derived hydrolysate .
We also measured decreased growth in hydrolysate upon repression of HOG1 (Fig. 3a). The mitogen-activated protein kinase (MAPKs) is a central component of the high osmolarity glycerol (HOG) pathway which mediates adaptation to hyperosmotic environments. We confirmed effects on HOG1 by dilution spot plating of an individual CRISPRi strain expressing a HOG1-targeting gRNA (AGGATCTTCGAAGGGAAGGA) with strong effect size in the screen (log2FCs of -3.35 in SCM + 10% Hydrolysate, compared to -0.62 in SCM) (Fig. 3b). In line with these results, the deletion of HOG1 in different strain backgrounds has been reported with reduced fitness in hydrolysate derived from corn stover . In addition, we identified repression of the HOG1-activating kinase PBS2 to result in decreased fitness in spruce hydrolysate, and mild growth defects are also known from the pbs2Δ mutant in Miscanthus plant hydrolysate . HOG pathway effects could be due to the high osmolarity of spruce hydrolysate as a result of dissolved compatible solutes which are released from lignocellulosic material during hydrolysis, including salts and carbohydrates.
In contrast, repression of the HOG-downstream TF SKO1 increased cell growth in the presence of hydrolysate. Sko1 is constitutively nuclear and bound to promoters for repression in non-stress condition [62, 68, 72]. Upon hyperosmotic stress which likely prevails in hydrolysate, Hog1 is activated and transitions into the nucleus where it phosphorylates Sko1 which is then exported to the cytosol, activating transcription via de-repression [62, 68, 72]. Interestingly, genes induced in lignocellulose hydrolysate are enriched for Sko1 target genes  which supports the relevance of Sko1-controlled transcripts in hydrolysate.
We notably observed increased growth in hydrolysate for repression of DOT6 (Fig. 3a) which encodes a transcriptional repressor that responds to osmotic and oxidative stress [52, 62]. The stimulated fitness in hydrolysate upon DOT6 repression can, akin to SKO1, be explained through de-repression. Proteins encoded by genes that affect growth in hydrolysate form physical interaction networks with bundled interactions at components of the osmotic (Hog1, Pbs2, Sko1) and oxidative stress response (Yap1, Dot6) and connect to growth regulators (Additional file 8: Fig. S6).
The mechanisms underlying further hydrolysate-specific effects of BUB1, CDC15 and UGA3 seem not directly linked to their known functions and offer leads for future investigation (Fig. 3a). CDC15 and BUB1 encode protein kinases involved in cell cycle control, e.g. with Bub1 hindering cell cycle progression if the spindle apparatus is damaged [29, 30]. UGA3 encodes a transcriptional activator of the gamma-aminobutyrate (GABA) pathway. Notably, the uga3 mutant was found with increased cellular fitness in presence of ethanol .
To better understand the cellular processes remodelled in hydrolysate, we determined the target genes of TFs that modulated fitness in hydrolysate using chromatin immunoprecipitation (ChIP) data  which were enriched for functions in translation, RNA binding, cyclic compound binding, the regulation of carbon metabolism (glycolysis and gluconeogenesis) and nuclear export of non-coding RNAs (Additional file 9: Fig. S7). We further searched for the phosphorylation targets of PKs modulating hydrolysate fitness, and found that these were enriched for functions in general kinase activity, Mitogen-activated protein kinase (MAPK) and target of rapamycin (TOR) signalling, mitophagy, stress granule components, chronic cell aging, as well as the cellular responses to osmotic stress, organic substances and acidic chemicals (Additional file 10: Fig. S8). Toxic chemicals, low pH and hyperosmotic conditions thus represent major stressors that yeast cells genetically brace against during growth in spruce hydrolysate.
Contributions of lignocellulosic growth inhibitors
We next screened for growth effects in the presence of eight lignocellulosic inhibitor compounds to determine genetic effects linked to their toxicity (IC composition in Table 1). After 25 generations selection in SCM + 45% Inhibitor Cocktail, we identified seven genes with inhibitor-specific functions compared to SCM, five of which overlap with hydrolysate-specific effects (Fig. 3c, Additional file 7: Table S2). Fitness effects observed for repression of HAA1, HOG1, PBS2, SKO1 and UGA3 in hydrolysate were not measured in presence of inhibitor compounds and are therefore likely independent of the toxicity caused by the used substances and concentrations (Fig. 3g). In contrast, the strong hydrolysate-specific growth effects of BUB1, DOT6, UME6 and YAP1 CRISPRi strains were reproduced with the inhibitor cocktail and were thus caused by one or multiple compounds in the cocktail (Fig. 3g). In addition, we find that the repression of STB5 reduced growth in inhibitor-supplemented and hydrolysate media although the gene barely missed the strict significance requirements in the hydrolysate screen (Fig. 3g). STB5 encodes a zinc cluster activator of pleiotropic multidrug resistance genes . In accord, stb5Δ strains are reported with decreased fitness in Miscanthus plant hydrolysate , while its overexpression was shown to increase growth in spruce hydrolysate . Three genes with inhibitor-specific effects are transcribed from bidirectional promoters which are hard to dissect by CRISPRi (CKA2|SLD7, STE7|DHH1, FKH2|YNL067W-A) since effects can be caused by the perturbation of either of the two genes, or their combined impact. Two of these loci (CKA2|SLD7, STE7|DHH1) were also significant in the hydrolysate screen. Clustering of gene fold changes across screens allows for the identification of specific and shared gene functions between conditions (Additional file 11: Fig. S9). Finally, we confirmed five strong loss of function fitness effects of our CRISPRi screens by using gene deletion strains of the prototrophic haploid library (BY4741 background) . We found hydrolysate growth decreased in hog1Δ, pbs2Δ and stb5Δ strains, and increased in dot6Δ and sko1Δ strains, recapitulating CRISPRi screen results (Additional file 12: Fig. S10).
Here, we employed CRISPR interference screens to identify regulatory genes capable of adjusting the growth of baker’s yeast in spruce hydrolysate and in the presence of lignocellulosic toxins. This allowed us to explain contributions of toxicity in hydrolysate growth conditions, and how genetic screens can be utilized towards overcoming current challenges in hydrolysate fermentation.
CRISPRi perturbations are powerful to probe genotype–phenotype relations in high throughput, with low cost and labour, and with high reproducibility. Our single plasmid CRISPRi system has been deeply characterized [33, 34, 78, 79], is freely available on addgene (#73796) and supported by an in-house developed gRNA design platform (http://lp2.github.io/yeast-crispri/)  and a customizable computational analysis pipeline (Methods) . The plasmid system can be transformed into any strain background, including polyploid and industrial strains, to devise strategies for improving process performance.
The presented CRISPRi screens on yeast growth in hydrolysate complement genetic screens with deletion mutants [65, 77]. We found that hydrolysate-specific functions are frequently connected to stress adaptation in the responses to oxidative (Yap1, Stb5), osmotic (Hog1, Pbs2, Sko1), acidic (Haa1) and general stress (Dot6, Ume6). These mirror the environmental conditions prevailing in the hydrolysate which are presumably perceived as stressful by yeast cells. Overcoming them can partially be tackled by the pretreatment of hydrolysate by, e.g. de-salting, pH adjustment or by addition of reducing agents to quench reactive oxygen species. Alternatively, or in addition, yeast strains can be genetically engineered to enhance growth, stress tolerance and productivity. Here we demonstrate that CRISPRi screens provide the opportunity to identify gene targets for strain optimization and, notably, multiple of our hits were already known to affect yeast growth in hydrolysate and have been successfully applied, including Yap1, Stb5 and Haa1. Overexpression of these three genes has been shown to increase tolerance to hydrolysate [13, 86] which not only validates their genetic functions, but also demonstrates their potential to improve bioprocesses.
To our knowledge, the presented study presents the first functional genetic screens of yeast in spruce hydrolysate. These screens suggest novel genes, expression of which can be modified to optimize hydrolysate fermentation, such as by overexpression of CDC15 and UGA3, or by decreasing expression levels of BUB1, DOT6 and SKO1. In addition, our finding that repression of the two key HOG signalling components (PBS2 and HOG1) result in growth defects in hydrolysate suggests an important regulatory role of this pathway during hydrolysate fermentation. We employed five haploid BY4741 deletion collection strains  to validate a handful of CRISPRi screen effects, demonstrating that hydrolysate fitness is enhanced in dot6Δ and sko1Δ mutants, and reduced in hog1Δ, pbs2Δ and stb5Δ backgrounds (Additional file 12: Fig. S10). These preliminary results warrant further investigation into the role of these genes in hydrolysate tolerance.
While our aim was to characterize genetics underlying fermentative growth, employing the screen under other selection conditions, such as low pH , high ethanol , and high temperatures [6, 17, 31, 34, 83], could be used to further improve specific processes.
Since growth-inhibiting agents are described as major challenge for utilizing plant hydrolysates for fermentation, we screened for fitness effects in presence of eight such compounds covering aliphatic acids, furan aldehydes and phenolic/aromatic molecules . While our cocktail was defined to largely cover these prevalent impacts, there are presumably additional lignocellulosic substances with yet other mechanisms of toxicity. Our finding that at least one-third of hydrolysate-specific gene functions were explained by toxicity effects therefore represents a lower bound and likely an underestimation. Due to the high prevalence of toxicity effects, relieving them through genetic regulation is an attractive avenue to facilitate the fermentation of lignocellulosic material.
Taken together, we show how CRISPRi screens can be used to identify genetic elements underlying complex environmental conditions encountered by cells in industrial bioprocesses. Our study pinpoints genetic functions that can be engineered to facilitate utilization of lignocellulose as feedstock for yeast fermentation, and thereby hopefully motivates and contributes to the establishment of environmentally sustainable procedures in biotechnology.
Materials and methods
Yeast strains, bacterial strains and plasmids
Two gRNA libraries targeting sets of 129 protein kinases  with 688 gRNAs and 161 transcription factors  with 885 gRNAs were used as described in Jann et al. . These DNA oligonucleotide libraries were integrated in the dCas9-Mxi1 plasmid system at the NotI restriction site by Gibson Assembly  as described in Smith et al. , transformed into the diploid BY4743 strain background as outlined by Gietz and Schiestl  and pooled together, yielding a final library with 1573 gRNAs. Growth phenotypes of five mutants (dot6Δ, hog1Δ, pbs2Δ, sko1Δ, stb5Δ) from the haploid BY4741 deletion collection  and the BY4741 wild type were measured. (Additional file 13: AF1 lists all chemical compounds, oligonucleotides, plasmids, bacterial and yeast strains, as well as all gRNAs used in this study).
Filter-sterilized synthetic complete uracil-dropout medium (SCM-Ura, for simplicity referred to as “SCM” in the manuscript) containing 20 g/L glucose, 6.7 g/L yeast nitrogen base without amino acids and with ammonium sulphate and 2 g/L amino acids uracil-dropout mix, pH 5.5, was used for fermentation. SCM-Ura was supplemented with hydrolysate or inhibitor cocktail mixture as indicated. Spruce hydrolysate was kindly provided by SEKAB (Örnsköldsvik, Sweden).
The sawdust raw material used to generate the hydrolysate contains 13 g/kg dry solid/dry weight (DS) arabinose, 20 g/kg DS galactose, 408 g/kg DS glucose, 50 g/kg DS xylose, 108 g/kg DS mannose, 0.26% ash (525 °C), 2.3% cyclohexane/acetone soluble matter and 27% DS lignin (Klason method) as reported by SEKAB. The sawdust with a dry matter content of 50% was pretreated with sulphur dioxide, enzymatically hydrolysed and filtered to remove solids.
The resulting hydrolysate is of defined composition and has been chemically characterized with 83 g/L glucose and higher order carbohydrates (26 g/L mannose, 9 g/L xylose, and less than 4 g/L galactose and arabinose), as well as toxic compounds that include, among others, 4.7 g/L acetic acid, 3.4 g/L 5-hydroxymethyl-furfural, 1.2 g/L furfural and less than 1 g/L of phenolic derivatives, levulinic and formic acid . In CRISPRi screening, SCM media was supplemented with 10% of the described hydrolysate of pH 4.5.
We defined an inhibitor cocktail (Table 1) based on previous characterizations of spruce and other plant hydrolysates, focussing particularly on medium produced by SEKAB [15, 36, 46, 47, 81, 86] or generated under similar technical and chemical conditions [51, 54, 67]. The inhibitor mix was generated to impose growth defects of all three major inhibitor compound types (aliphatic acids, furan aldehydes and phenolic/aromatic molecules) and their concentrations are not proportional to the ones of the used hydrolysate from SEKAB. Inhibitor compounds were chosen based on described effect size and occurrence in hydrolysates based on literature. Inhibitor compounds were purchased from Sigma-Aldrich, and mixed to a 5× concentrated cocktail before addition to SCM (Table 1). Coniferyl aldehyde, vanillin and p-hydroxybenzaldehyde were diluted in 3 ml 0.1 M NaOH prior to addition. Inhibitor-supplemented SCM-Ura was adjusted to pH 5 with 0.1 M NaOH.
Quantitative determination of acetate, ethanol and glucose was performed by high-performance liquid chromatography (HPLC). The HPLC system was equipped with a refractive index detector (Alliance HPLC with 2414 RID, Waters, Eschborn, Germany) and a Rezex ROA-Organic Acid H + (8%) column (Phenomenex, Aschaffenburg, Germany) which was held at 65 °C during all measurements. 0.5 mM sulfuric acid was used as mobile phase carrying the samples at a flow rate of 0.5 ml/min. Samples were stored at 6 °C prior to each run which took 15 min. Hydrolysate-containing samples were diluted four times. Compounds were quantified by comparing the metabolite peak in the sample with a mixture of standards with known concentrations of each metabolite.
Dilution spot plating
CRISPRi cultures were grown in SCM and SCM + 10% hydrolysate, either with or without addition of 250 ng/mL ATc for 24 h at 30 °C. A dilution series was prepared in medium without ATc, and 10 µL of the dilutions were plated on SCM or SCM + 10% hydrolysate agar plates. Photographs were taken after 2 days incubation at 30 °C to report on colony size.
For HPLC experiments, yeast cultures were grown in 50 mL in shake flasks at 30 °C with 175 rpm. Flasks were sealed with rubber plugs to mimic oxygen-limiting conditions, stimulating fermentation rather than respiration. Plugs were pierced with sterile needles (27, BD Microlance, Becton Dickinson) and cotton plugs to enable CO2 exchange. Plate reader assays were performed under oxygen-limited conditions in 96-well plates (Nunclon Delta Surface, Thermo Scientific) with Synergy HTX Multi-Detection Microplate Readers (BioTek) at 30 °C. Plates were inoculated with an optical density at wavelength 600 nm (OD600) of 0.005. OD600 was measured in 15 min intervals and shaken at 800 rpm for 10 s before measurements. For screens, yeast was cultured in sealed falcon tubes as described below.
Competitive growth assays
For screens, defrosted hydrolysate and inhibitor cocktail were diluted in SCM-Ura to 10% and 45%, respectively. Screens were performed in triplicate samples, each with addition of 250 ng/mL ATc to induce gRNA expression and without ATc as reference. Yeast cultures were profiled in 15 mL falcon tubes containing 11 mL medium that was inoculated with overnight cultures at OD 0.005. Falcon tubes were sealed to mimic oxygen-limited conditions, and a needle (20 G, BD Microlance) stuffed with cotton was pierced through the lid. Cells were then grown in a shaking incubator at 30 °C and 180 rpm. The tubes were opened only to assess growth stage. Cells were transferred to fresh medium in late mid-exponential phase (before reaching OD600 = 1). During the screen, two such transfers were performed to maintain exponential growth. Cell samples were centrifuged, and pellets were washed and used to extract plasmid DNA.
Plasmid DNA was purified using the Miniprep kit (QIAprep Spin, Qiagen, Hilden, Germany) with modified protocol. Cell pellets were resuspended in P1 solution (accordingly to kit manufacturer) and incubated with 9U lyticase (Sigma-Aldrich) at 37 °C for 30 min, followed by harsh vortexing for 2 min. The remaining steps were performed following the provided kit protocol. Illumina sequencing adapters and inline barcodes were introduced to DNA barcodes by PCR, creating a specific double index of samples. All PCR products were analysed by gel electrophoresis and then purified (QIAquick, Quiagen). Samples were pooled to similar amounts. The resulting sequencing library was concentrated by performing another PCR purification, and a 1% agarose gel electrophoresis was performed for size-selection and purification. DNA quality was checked with a DNA-Bioanalyzer, and a diluted sample was sequenced on an Illumina Next-Seq500 machine in paired-end mode with 75 bases read length.
Data processing and analysis
Raw Illumina sequencing reads were demultiplexed using Jemultiplexer. Trimmed reads were aligned to a gRNA sequence reference genome with Burrows-Wheeler aligner to compute read counts. Read counts were processed in R code with the edgeR package, using a generalized linear model to compute the log2 guide and gene fold change (log2FC) between + ATc and -ATc populations that both went through selection. Significant genes were required to have a gene log2FC with FDR < 0.05 and at least two supporting gRNAs with FDR < 0.05 and absolute log2FC ≥ 1.
Gene ontology-enrichment was performed using the gProfiler2 R package . TF target genes were determined using Chromatin Immuno-Precipitation on chip (ChIP-chip) data from Gonçalves et al. . Genes bound by two or more TFs from the set of significant modulators identified in CRISPRi screens were used for enrichment analysis with the S. cerevisiae genome as statistic background. Phosphorylation targets of protein kinases were determined with data from phosphogrid 2.0 . The PK targets of significant modulators of CRISPRi screens were used for enrichment analysis, using all protein kinase targets measured in the phosphoproteomics data set as statistic background.
Plots were generated in R (V. 3.4.1)  with the ggplot2 (V. 3.1.0) , ggally (V. 1.3)  and pheatmap (V. 1.0.10) packages . In boxplots, the middle line denotes the median, and lower and upper hinges denote the first and third quartiles, respectively. Figures were designed in Adobe Illustrator 2019.
Availability of data and materials
Demultiplexed sequencing data, read counts, gRNA fold changes and gene fold changes were deposited at Gene Expression Omnibus and are accessible under GSE155590 with the token chkxawwunlgrjeb. Computed gRNA barcode and gene level fold changes are provided in the Additional file 14: AF2 and Additional file 15: AF3, respectively.
Clustered Regularly Interspaced Short Palindromic Repeats interference
(Dead) Clustered Regularly Interspaced Short Palindromic Repeats-associated protein 9
High-performance liquid chromatography
Inhibitor compound cocktail
Synthetic complete medium
Allen SA, Clark W, McCaffery JM, Cai Z, Lanctot A, Slininger PJ, Lewis Liu Z, Gorsich SW. Furfural induces reactive oxygen species accumulation and cellular damage in Saccharomyces cerevisiae. Biotechnol Biofuels. 2010;3(1):2.
Almeida JR, Modig T, Petersson A, Hähn-Hägerdal B, Lidén G, Gorwa-Grauslund MF. Increased tolerance and conversion of inhibitors in lignocellulosic hydrolysates by Saccharomyces cerevisiae. J Chem Technol Biotechnol Int Res Process Environ Clean Technol. 2007;82(4):340–9.
Alriksson B, Cavka A, Jönsson LJ. Improving the fermentability of enzymatic hydrolysates of lignocellulose through chemical in-situ detoxification with reducing agents. Biores Technol. 2011;102(2):1254–63.
Akache B, MacPherson S, Sylvain MA, Turcotte B. Complex interplay among regulators of drug resistance genes in Saccharomyces cerevisiae. J Biol Chem. 2004;279(27):27855–60.
Azhar SHM, Abdulla R, Jambo SA, Marbawi H, Gansau JA, Faik AAM, Rodrigues KF. Yeasts in sustainable bioethanol production: a review. Biochem Biophys Rep. 2017;10:52–61.
Banat IM, Nigam P, Singh D, Marchant R, McHale AP. Ethanol production at elevated temperatures and alcohol concentrations: Part I-Yeasts in general. World J Microbiol Biotechnol. 1998;14(6):809–21.
Banerjee N, Bhatnagar R, Viswanathan L. Inhibition of glycolysis by furfural in Saccharomyces cerevisiae. Eur J Appl Microbiol Biotechnol. 1981;11(4):226–31.
Bastian M, Heymann S, Jacomy M. Gephi: an open source software for exploring and manipulating networks. In: Third international AAAI conference on weblogs and social media; 2009.
Berry DR, Russell I, Stewart GC. Yeast biotechnology. Berlin: Springer Science & Business Media; 2012.
Breitkreutz A, Choi H, Sharom JR, Boucher L, Neduva V, Larsen B, Lin ZY, Breitkreutz BJ, Stark C, Liu G, Ahn J. A global protein kinase and phosphatase interaction network in yeast. Science. 2010;328(5981):1043–6.
Brown JA, Sherlock G, Myers CL, Burrows NM, Deng C, Wu HI, McCann KE, Troyanskaya OG, Brown JM. Global analysis of gene function in yeast by quantitative phenotypic profiling. Mol Syst Biol. 2006;2(1):2006–0001.
Cherry JM, Adler C, Ball C, Chervitz SA, Dwight SS, Hester ET, Jia Y, Juvik G, Roe TY, Schroeder M, Weng S, Botstein D. SGD: Saccharomyces genome database. Nucleic Acids Res. 1998;26(1):73–9.
Cunha JT, Costa CE, Ferraz L, Romaní A, Johansson B, Sá-Correia I, Domingues L. HAA1 and PRS3 overexpression boosts yeast tolerance towards acetic acid improving xylose or glucose consumption: unravelling the underlying mechanisms. Appl Microbiol Biotechnol. 2018;102(10):4589–600.
De Boer CG, Hughes TR. YeTFaSCo: a database of evaluated yeast transcription factor sequence specificities. Nucleic Acids Res. 2012;40(D1):D169–79.
Demeke MM, Dietz H, Li Y, Foulquié-Moreno MR, Mutturi S, Deprez S, Abt TD, Bonini BM, Liden G, Dumortier F, Verplaetse A, Boles E, Verplaetse A. Development of a D-xylose fermenting and inhibitor tolerant industrial Saccharomyces cerevisiae strain with high performance in lignocellulose hydrolysates using metabolic and evolutionary engineering. Biotechnol Biofuels. 2013;6(1):89.
Dudley AM, Janse DM, Tanay A, Shamir R, Church GM. A global view of pleiotropy and phenotypically derived gene function in yeast. Mol Syst Biol. 2005;1(1):2005–0001.
Edgardo A, Carolina P, Manuel R, Juanita F, Baeza J. Selection of thermotolerant yeast strains Saccharomyces cerevisiae for bioethanol production. Enzyme Microbial Technol. 2008;43(2):120–3.
Farzadfard F, Perli SD, Lu TK. Tunable and multifunctional eukaryotic transcription factors based on CRISPR/Cas. ACS Synth Biol. 2013;2(10):604–13.
Fernandes AR, Mira NP, Vargas RC, Canelhas I, Sá-Correia I. Saccharomyces cerevisiae adaptation to weak acids involves the transcription factor Haa1p and Haa1p-regulated genes. Biochem Biophys Res Commun. 2005;337(1):95–103.
Francois JM, Alkim C, Morin N. Engineering microbial pathways for production of bio-based chemicals from lignocellulosic sugars: current status and perspectives. Biotechnol Biofuels. 2020;13(1):1–23.
Fujita K, Matsuyama A, Kobayashi Y, Iwahashi H. The genome-wide screening of yeast deletion mutants to identify the genes required for tolerance to ethanol and other alcohols. FEMS Yeast Res. 2006;6(5):744–50.
Giaever G, Chu AM, Ni L, Connelly C, Riles L, Véronneau S, Dow S, Lucau-Danila A, Anderson K, André B, Arkin AP, Astromoff A, El Bakkoury M, Bangham R, Benito R, Brachat S, Campanaro S, Curtiss M, Davis K, Deutschbauer A, Entian K, Flaherty P, Foury F, Garfinkel DJ, Gerstein M, Gotte D, Güldener U, Hegemann JH, Hempel S, Herman Z, Jaramillo DF, Kelly DE, Kelly SL, Kötter P, LaBonte D, Lamb DC, Lan N, Liang H, Liao H, Liu L, Luo C, Lussier M, Mao R, Menard P, Ooi SL, Revuelta JL, Roberts CJ, Rose M, Ross-Macdonald P, Scherens B, Schimmack G, Shafer B, Shoemaker DD, Sookhai-Mahadeo S, Storms RK, Strathern JN, Valle G, Voet M, Volckaert G, Wang C, Ward TR, Wilhelmy J, Winzeler EA, Yang Y, Yen G, Youngman E, Yu K, Bussey H, Boeke JD, Snyder M, Philippsen P, Davis RW, Johnston M. Functional profiling of the Saccharomyces cerevisiae genome. Nature. 2002;418(6896):387–91.
Gibson DG, Young L, Chuang RY, Venter JC, Hutchison CA, Smith HO. Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat Methods. 2009;6(5):343–5.
Gietz RD, Schiestl RH. High-efficiency yeast transformation using the LiAc/SS carrier DNA/PEG method. Nat Protoc. 2007;2(1):31–4.
Gilbert LA, Larson MH, Morsut L, Liu Z, Brar GA, Torres SE, Stern-Ginossar N, Brandmann O, Whitehead EH, Doudna JA, Lim WA, Weissman JS, Qi LS. CRISPR-mediated modular RNA-guided regulation of transcription in eukaryotes. Cell. 2013;154(2):442–51.
Gilbert LA, Horlbeck MA, Adamson B, Villalta JE, Chen Y, Whitehead EH, Guimaraes C, Panning B, Ploegh HL, Bassik MC, Qi LS, Kampmann M, Weissman JS. Genome-scale CRISPR-mediated control of gene repression and activation. Cell. 2014;159(3):647–61.
Ginestet C. ggplot2: elegant graphics for data analysis. J R Stat Soc Ser A Stat Soc. 2011;174:245–6.
Gonçalves E, Nakic ZR, Zampieri M, Wagih O, Ochoa D, Sauer U, Beltrao P, Saez-Rodriguez J. Systematic analysis of transcriptional and post-transcriptional regulation of metabolism in yeast. PLoS Comput Biol. 2017;13(1):e1005297.
Goto GH, Mishra A, Abdulle R, Slaughter CA, Kitagawa K. Bub1-mediated adaptation of the spindle checkpoint. PLoS Genet. 2011;7(1):e1001282.
Hartwell LH, Culotti J, Reid B. Genetic control of the cell-division cycle in yeast, I. Detection of mutants. Proc Natl Acad Sci. 1970;66(2):352–9.
Hasunuma T, Kondo A. Consolidated bioprocessing and simultaneous saccharification and fermentation of lignocellulose to ethanol with thermotolerant yeast strains. Process Biochem. 2012;47(9):1287–94.
Jacomy M, Venturini T, Heymann S, Bastian M. ForceAtlas2, a continuous graph layout algorithm for handy network visualization designed for the Gephi software. PLoS ONE. 2014;9(6):e98679.
Jaffe M, Dziulko A, Smith JD, Onge RPS, Levy SF, Sherlock G. Improved discovery of genetic interactions using CRISPRiSeq across multiple environments. Genome Res. 2019;29(4):668–81.
Jann C, Johansson A, Smith JD, Parts L, Steinmetz LM. Gene dosage screens in yeast reveal core signalling pathways controlling heat adaptation. bioRxiv. 2020. https://doi.org/10.1101/2020.08.26.267674.
Jansen ML, Bracher JM, Papapetridis I, Verhoeven MD, de Bruijn H, de Waal PP, van Maris A, Klaassen P, Pronk JT. Saccharomyces cerevisiae strains for second-generation ethanol production: from academic exploration to industrial implementation. FEMS Yeast Res. 2017;17(5):fox044.
Johansson E, Xiros C, Larsson C. Fermentation performance and physiology of two strains of Saccharomyces cerevisiae during growth in high gravity spruce hydrolysate and spent sulphite liquor. BMC Biotechnol. 2014;14(1):47.
Jönsson LJ, Alriksson B, Nilvebrant N-O. Bioconversion of lignocellulose: inhibitors and detoxification. Biotechnology for Biofuels. 2013;6(1):16.
Jöonsson LJ, Alriksson B, Nilvebrant N-O. Bioconversion of lignocellulose: inhibitors and detoxification. Biotechnol Biofuels. 2013;6(1):16.
Keng T. HAP1 and ROX1 form a regulatory pathway in the repression of HEM13 transcription in Saccharomyces cerevisiae. Mol Cell Biol. 1992;12(6):2616–23.
Keweloh H, Weyrauch G, Rehm H-J. Phenol-induced membrane changes in free and immobilized Escherichia coli. Appl Microbiol Biotechnol. 1990;33(1):66–71.
Kim D, Hahn JS. Roles of the Yap1 transcription factor and antioxidants in Saccharomyces cerevisiae’s tolerance to furfural and 5-hydroxymethylfurfural, which function as thiol-reactive electrophiles generating oxidative stress. Appl Environ Microbiol. 2013;79(16):5069–77.
Koivistoinen OM, Kuivanen J, Barth D, Turkia H, Pitkänen JP, Penttilä M, Richard P. Glycolic acid production in the engineered yeasts Saccharomyces cerevisiae and Kluyveromyces lactis. Microb Cell Fact. 2013;12(1):82.
Kolberg L, Raudvere U, Kuzmin I, Vilo J, Peterson H. gprofiler2 - an R package for gene list functional enrichment analysis and namespace conversion toolset g:Profiler. F1000Research. 2020;9:709.
Kolde, R. Package `pheatmap’. Bioconductor. 2012.
Kolde R. Pheatmap: Pretty Heatmaps R. Package. Bioconductor; CRAN-project.org/package=pheatmap; 2012.
Koppram R, Albers E, Olsson L. Evolutionary engineering strategies to enhance tolerance of xylose utilizing recombinant yeast to inhibitors derived from spruce biomass. Biotechnol Biofuels. 2012;5(1):32.
Koppram R, Olsson L. Combined substrate, enzyme and yeast feed in simultaneous saccharification and fermentation allow bioethanol production from pretreated spruce biomass at high solids loadings. Biotechnol Biofuels. 2014;7(1):54.
Kuge S, Jones N, Nomoto A. Regulation of YAP-1 nuclear localization in response to oxidative stress. EMBO J. 1997;16(7):1710–20.
Kwast KE, Burke PV, Poyton RO. Oxygen sensing and the transcriptional regulation of oxygen-responsive genes in yeast. J Exp Biol. 1998;201(8):1177–95.
Lardenois A, Becker E, Walther T, Law MJ, Xie B, Demougin P, Strich R, Primig M. Global alterations of the transcriptional landscape during yeast growth and development in the absence of Ume6-dependent chromatin modification. Mol Genet Genomics. 2015;290(5):2031–46.
Larsson S, Reimann A, Nilvebrant NO, Jönsson LJ. Comparison of different methods for the detoxification of lignocellulose hydrolyzates of spruce. Appl Biochem Biotechnol. 1999;77(1–3):91–103.
Lippman SI, Broach JR. Protein kinase A and TORC1 activate genes for ribosomal biogenesis by inactivating repressors encoded by Dot6 and its homolog Tod6. Proc Natl Acad Sci. 2009;106(47):19928–33.
Lopes DD, Rosa CA, Hector RE, Dien BS, Mertens JA, Ayub MAZ. Influence of genetic background of engineered xylose-fermenting industrial Saccharomyces cerevisiae strains for ethanol production from lignocellulosic hydrolysates. J Ind Microbiol Biotechnol. 2017;44(11):1575–88.
Martín C, Wu G, Wang Z, Stagge S, Jönsson LJ. Formation of microbial inhibitors in steam-explosion pretreatment of softwood impregnated with sulfuric acid and sulfur dioxide. Biores Technol. 2018;262:242–50.
McGlincy NJ, Meacham ZA, Swain K, Muller R, Baum R, Ingolia N. A genome-scale CRISPR interference guide library enables comprehensive phenotypic profiling in yeast. BiorXiv. 2020. https://doi.org/10.1101/2020.03.11.988105.
Momen-Roknabadi A, Oikonomou P, Tavazoie S. An inducible CRISPR-interference library for genetic interrogation of Saccharomyces cerevisiae biology. bioRxiv. 2020;3:1–2.
Mülleder M, Capuano F, Pir P, Christen S, Sauer U, Oliver SG, Ralser M. A prototrophic deletion mutant collection for yeast metabolomics and systems biology. Nat Biotechnol. 2012;30(12):1176–8.
Nguyen TTM, Iwaki A, Ohya Y, Izawa S. Vanillin causes the activation of Yap1 and mitochondrial fragmentation in Saccharomyces cerevisiae. J Biosci Bioeng. 2014;117(1):33–8.
Ostergaard S, Olsson L, Nielsen J. Metabolic engineering of Saccharomyces cerevisiae. Microbiol Mol Biol Rev. 2000;64(1):34–50.
Okazaki S, Tachibana T, Naganuma A, Mano N, Kuge S. Multistep disulfide bond formation in Yap1 is required for sensing and transduction of H2O2 stress signal. Mol Cell. 2007;27(4):675–88.
Pampulha M, Loureiro-Dias M. Combined effect of acetic acid, pH and ethanol on intracellular pH of fermenting yeast. Appl Microbiol Biotechnol. 1989;31(5–6):547–50.
Pascual-Ahuir A, Posas F, Serrano R, Proft M. Multiple levels of control regulate the yeast cAMP-response element-binding protein repressor Sko1p in response to stress. J Biol Chem. 2001;276(40):37373–8.
Patel A, Arora N, Sartaj K, Pruthi V, Pruthi PA. Sustainable biodiesel production from oleaginous yeasts utilizing hydrolysates of various non-edible lignocellulosic biomasses. Renew Sustain Energy Rev. 2016;62:836–55.
Parawira W, Tekere M. Biotechnological strategies to overcome inhibitors in lignocellulose hydrolysates for ethanol production. Crit Rev Biotechnol. 2011;31(1):20–31.
Pereira FB, Romaní A, Ruiz HA, Teixeira JA, Domingues L. Industrial robust yeast isolates with great potential for fermentation of lignocellulosic biomass. Biores Technol. 2014;161:192–9.
Pérez J, Munoz-Dorado J, De la Rubia TDLR, Martinez J. Biodegradation and biological treatments of cellulose, hemicellulose and lignin: an overview. Int Microbiol. 2002;5(2):53–63.
Petersson A, Lidén G. Fed-batch cultivation of Saccharomyces cerevisiae on lignocellulosic hydrolyzate. Biotech Lett. 2007;29(2):219–25.
Proft M, Serrano R. Repressors and upstream repressing sequences of the stress-regulated ENA1 gene in Saccharomyces cerevisiae: bZIP protein Sko1p confers HOG-dependent osmotic regulation. Mol Cell Biol. 1999;19(1):537–46.
Qian W, Ma D, Xiao C, Wang Z, Zhang J. The genomic landscape and evolutionary resolution of antagonistic pleiotropy in yeast. Cell Rep. 2012;2(5):1399–410.
R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing; 2018.
Ragauskas AJ, Nagy M, Kim DH, Eckert CA, Hallett JP, Liotta CL. From wood to fuels: integrating biofuels and pulp production. Ind Biotechnol. 2006;2(1):55–65.
Rep M, Proft M, Remize F, Tamás M, Serrano R, Thevelein JM, Hohmann S. The Saccharomyces cerevisiae Sko1p transcription factor mediates HOG pathway-dependent osmotic regulation of a set of genes encoding enzymes implicated in protection from oxidative damage. Mol Microbiol. 2001;40(5):1067–83.
Romaní A, Pereira F, Johansson B, Domingues L. Metabolic engineering of Saccharomyces cerevisiae ethanol strains PE-2 and CAT-1 for efficient lignocellulosic fermentation. Biores Technol. 2015;179:150–8.
Sadowski I, Breitkreutz BJ, Stark C, Su TC, Dahabieh M, Raithatha S, Bernhard W, Oughtred R, Dolinski K, Barreto K, Tyers M. The PhosphoGRID Saccharomyces cerevisiae protein phosphorylation site database: version 2.0 update. Database. 2013. https://doi.org/10.1038/s42003-020-01452-9.
Sardi M, Rovinskiy N, Zhang Y, Gasch AP. Leveraging genetic-background effects in Saccharomyces cerevisiae to improve lignocellulosic hydrolysate tolerance. Appl Environ Microbiol. 2016;82(19):5838–49.
Schloerke B, Crowley J, Cook D, Hofmann H, Hadley W, Briatte F, Marbach M, Thoen E, Elberg A, Larmarange J. GGally: Extension to ‘ggplot2’. R Package version 1.3.2 https://CRAN.R-project.org/package=GGally. 2017.
Skerker JM, Leon D, Price MN, Mar JS, Tarjan DR, Wetmore KM, Deutschbauer AM, Baumohl JK, Bauer S, Ibánez AB, Mitchell VD. Dissecting a complex chemical stress: chemogenomic profiling of plant hydrolysates. Mol Syst Biol. 2013;9(1):674.
Smith JD, Suresh S, Schlecht U, Wu M, Wagih O, Peltz G, Davis RW, Steinmetz LM, Parts L, Onge RPS. Quantitative CRISPR interference screens in yeast identify chemical-genetic interactions and new rules for guide RNA design. Genome Biol. 2016;17(1):45.
Smith JD, Schlecht U, Xu W, Suresh S, Horecka J, Proctor MJ, Aijar R, Bennet RAO, Chu A, Li YF, Roy K, Davis RW, Steinmetz LM, Hyman RW, Levy SF, Onge RPS. A method for high-throughput production of sequence-verified DNA libraries and strain collections. Mol Syst Biol. 2017;13(2):913.
Snel B, Lehmann G, Bork P, Huynen MA. STRING: a web-server to retrieve and display the repeatedly occurring neighbourhood of a gene. Nucleic Acids Res. 2000;28(18):3442–4.
Soudham VP, Brandberg T, Mikkola JP, Larsson C. Detoxification of acid pretreated spruce hydrolysates with ferrous sulfate and hydrogen peroxide improves enzymatic hydrolysis and fermentation. Biores Technol. 2014;166:559–65.
Sun Y, Cheng J. Hydrolysis of lignocellulosic materials for ethanol production: a review. Biores Technol. 2002;83(1):1–11.
Van Vleet JH, Jeffries TW. Yeast metabolic engineering for hemicellulosic ethanol production. Curr Opin Biotechnol. 2009;20(3):300–6.
Wagner ER, Myers KS, Riley NM, Coon JJ, Gasch AP. PKA and HOG signaling contribute separable roles to anaerobic xylose fermentation in yeast engineered for biofuel production. PLoS ONE. 2019;14(5):e0212389.
Willför S, Holmbom B. Isolation and characterisation of water soluble polysaccharides from Norway spruce and Scots pine. Wood Sci Technol. 2004;38(3):173–9.
Wu G, Xu Z, Jönsson LJ. Profiling of Saccharomyces cerevisiae transcription factors for engineering the resistance of yeast to lignocellulose-derived inhibitors in biomass conversion. Microb Cell Fact. 2017;16(1):199.
Zhu JY, Pan XJ, Wang GS, Gleisner R. Sulfite pretreatment (SPORL) for robust enzymatic saccharification of spruce and red pine. Biores Technol. 2009;100(8):2411–8.
Zitomer RS, Lowry CV. Regulation of gene expression by oxygen in Saccharomyces cerevisiae. Microbiol Mol Biol Rev. 1992;56(1):1–11.
The authors wish to thank Antonius J.A. van Maris (Royal Institute of Technology, Stockholm, Sweden) for his supervisory support of the Master’s Thesis of F.G. Adnan Cavka from SEKAB (Örnsköldsvik, Sweden) is gratefully acknowledged for the kind provision of pretreated lignocellulose hydrolysate and technical advice. The authors very much appreciated the useful discussions with Maurizio Bettiga (Chalmers University of Technology, Gothenburg, Sweden) regarding the experimental setup and Justin Smith (Stanford University) for CRISPRi construct setup.
Open Access funding enabled and organized by Projekt DEAL. This study was supported by an Advanced Investigator grant from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (AdG-742804 to L.M.S.), by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation; project STE 1422/4-1 to L.M.S.), and a Joachim Herz Add-On Fellowship (to C.J.).
Ethics approval and consent to participate
Consent for publication
All authors approved the manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Yeast growth across media in 96-well plate format. The optical density at 600 nm (OD600, on y-axis) was quantified over time (x-axis) during growth of BY4743 in synthetic complete media with 2% glucose (SCM), as well as in SCM that was supplemented with different concentrations of (a) ethanol, (b) spruce hydrolysate or (c) inhibitor compound cocktail. The respective supplement concentrations are indicated individually. For the IC mixture, the 1x IC stock was diluted to the indicated percentages. The curves denote the average of n = 4 wells measured in 96-well format, normalized by subtraction of media background.
HPLC measurements. Metabolites concentrations (g L−1), measured by HPLC, of yeast cultures grown in SCM, SCM+10% hydrolysate and SCM+45% inhibitor cocktail. Glucose, ethanol and acetic acid concentrations were measured in the growth medium (Initial) and at the end of fermentations (Final). Three biological triplicates were performed for the three tested conditions, error represents the standard deviation between replicates.
Acetic acid metabolization. Changes of acetic acid (in g L−1) concentration during fermentation in different growth conditions (indicated in figure legend) at cultivation start and end points, measured by HPLC. Initial data corresponds to media used to inoculate, while each point in “Final” correspond to acetic acid concentration of three biological replicates.
Read count correlation. Spearman correlations of read count samples across screens.
Guide RNA fold changes across conditions. Scatter plots with dots denoting gRNAs, density distributions and Pearson correlations of gRNA log2 fold changes across screen conditions.
Gene fold changes across conditions. Scatter plots with dots denoting genes, density distributions and Pearson correlations of gene log2 fold changes across screen conditions. Line denotes smoothed linear fits.
CRISPRi effects across media. Significant genes (gene level fold change with FDR < 0.05 and at least two gRNAs with absolute log2FC ≥ 1 and FDR < 0.05) are shown across screens with their mean log2FC, their maximum gRNA log2FC and ID to specify if a TF or PK is targeted. For target genes transcribed from bidirectional promoters, both genes are reported (separated with a vertical dash). Tables are ordered by gene log2FC. Rows of essential genes as defined by in-viable knock-out mutants  are in green color. One and two asterisks (*, **) behind a gene name indicate that repression caused hydrolysate-specific or inhibitor-specific effects, respectively (not measured in SCM).
Protein–protein interaction network between modulators of hydrolysate growth. Experimental protein–protein interactions of genes modulating cellular fitness in hydrolysate, obtained from STRING . Dots denote genes, coloured by gradients from light to dark by increased strength in either positively (green) or negatively (red) modulating hydrolysate fitness, obtained from screen log2-fold changes. Dot and gene label size denote the multiple-testing adjusted FDRs obtained in the screen. Line thickness indicates confidence of the physical interaction obtained from the STRING database. Network visualization was performed with Gephi , using the Force Atlas 2 algorithm for clustering with standard parameters .
Hydrolysate-specific TF target gene functions. GO enrichment of TF target genes determined from ChIP-chip (Gonçalves et al., 2017) of TFs which modulate hydrolysate growth, generated using the gProfiler2 R package (Reimand et al. 2019).
Hydrolysate-specific PK interactor functions. GO enrichment of PK phosphorylation targets determined from Phospho-proteomics data  of PKs which modulate hydrolysate growth, generated using the gProfiler2 R package (Reimand et al. 2019).
CRISPRi effects across screens. Log2 gene fold changes compared between SC medium, SCM + 10% Hydrolysate and SCM + 45% Inhibitor Cocktail. The heatmap was generated with the pheatmap R package (Kolde 2019).
Growth profile in SCM and in SCM+10% Hydrolysate of prototrophic gene deletion strains. The optical density at 600 nm (OD600, on y-axis) was quantified over time (hour, x-axis) during growth of prototrophic BY4741 WT (grey) and the prototrophic BY4741 deletion strains (orange) in SCM and in SCM supplemented with 10% spruce hydrolysate. The curves denote the average of n = 3 wells measured in 96-well format, normalized by subtraction of media background.
Additional file 13: AF1. List of all chemical compounds, oligonucleotides, plasmids, bacterial and yeast strains, as well as all gRNAs sequences used in this study.
About this article
Cite this article
Gutmann, F., Jann, C., Pereira, F. et al. CRISPRi screens reveal genes modulating yeast growth in lignocellulose hydrolysate. Biotechnol Biofuels 14, 41 (2021). https://doi.org/10.1186/s13068-021-01880-7