The Absence of the Transcription Factor Yrr1p, Identified from Comparative Genome Profiling, Increased Vanillin Tolerance Due to Enhancements of ABC Transporters Expressing, rRNA Processing and Ribosome Biogenesis in Saccharomyces cerevisiae

Enhancing the tolerance of Saccharomyces cerevisiae to inhibitors derived from lignocellulose is conducive to producing biofuel and chemicals using abundant lignocellulosic materials. Vanillin is a major type of phenolic inhibitor in lignocellulose hydrolysates for S. cerevisiae. In the present work, the factors beneficial to vanillin resistance in yeast were identified from the vanillin-resistant strain EMV-8, which was derived from strain NAN-27 by adaptive evolution. We found 450 SNPs and 44 genes with InDels in the vanillin-tolerant strain EMV-8 by comparing the genome sequences of EMV-8 and NAN-27. To investigate the effects of InDels, InDels were deleted in BY4741, respectively. We demonstrated that the deletion of YRR1 improved vanillin tolerance of strain. In the presence of 6 mM vanillin, deleting YRR1 increase the maximum specific growth rate and the vanillin consumption rate by 142 and 51%, respectively. The subsequent transcriptome analysis revealed that deleting YRR1 resulted in changed expression of over 200 genes in the presence of 5 mM vanillin. The most marked changes were the significant up-regulation of the dehydrogenase ADH7, several ATP-binding cassette (ABC) transporters, and dozens of genes involved in ribosome biogenesis and rRNA processing. Coincidently, the crude enzyme solution of BY4741(yrr1Δ) exhibited higher NADPH-dependent vanillin reduction activity than control. In addition, overexpressing the ABC transporter genes PDR5, YOR1, and SNQ2, as well as the RNA helicase gene DBP2, increased the vanillin tolerance of strain. Interestingly, unlike the marked changes we mentioned above, under vanillin-free conditions, there are only limited transcriptional differences between wildtype and yrr1Δ. This indicated that vanillin might act as an effector in Yrr1p-related regulatory processes. The new findings of the relationship between YRR1 and vanillin tolerance, as well as the contribution of rRNA processing and ribosome biogenesis to enhancing S. cerevisiae vanillin tolerance, provide novel targets for genetic engineering manipulation to improve microbes’ tolerance to lignocellulose hydrolysate.

Introduction

Saccharomyces cerevisiae is widely used in traditional ethanol production and shows promise as a cell factory for other chemical production because it is easily genetically manipulable, robust, and generally recognized as safe (Liu et al., 2009). In recent decades, much effort has been invested in producing ethanol and valuable chemicals with S. cerevisiae from lignocellulosic materials because these materials are the abundant renewable biomass resource on earth (Kim et al., 2015). The pretreatment process is critical for the effective use of lignocellulose because it is necessary to destroy the complicated structure of lignocellulose and to release monosaccharides. However, toxic chemicals, such as organic acids, furans, and phenolics, are unavoidably produced during pretreatment, and they inhibit the growth and metabolism of microorganisms in the subsequent fermentation process (Klinke et al., 2004; Liu, 2011). Among these toxic chemicals, the phenolics, which are produced from the segmental degradation of lignin, exhibit more toxicity than organic acids and furans to fermentation by S. cerevisiae (Klinke et al., 2004). Lignocellulose hydrolysate contains three kinds of phenolics, classified by the radicals they contain (para-hydroxyphenyl, guaiacyl, and syringyl). Vanillin is a typical guaiacyl phenolic (Klinke et al., 2004; Wang et al., 2016). It has been confirmed that in S. cerevisiae, vanillin triggers ROS accumulation, fractures the mitochondria and represses translation (Iwaki et al., 2013b; Nguyen et al., 2014a,b). To date, there are three strategies to improve the vanillin tolerance of S. cerevisiae. First, vanillin can be converted to the less-toxic vanillyl alcohol by certain reductases and dehydrogenases; therefore, identifying efficient reductases and dehydrogenases and overexpressing them has proven effective (Liu, 2011). Second, strains having higher ergosterol synthesis capacity show higher vanillin tolerance because higher levels of ergosterol enhance the fluidity and stability of the membrane (Endo et al., 2009). Finally, because vanillin causes cellular ROS accumulation (Nguyen et al., 2014a), high glutathione levels increase yeast robustness (Ask et al., 2013). Further studies are necessary to identify other effective strategies for improving the resistance of S. cerevisiae to vanillin.

In our previous study, a high vanillin-tolerant S. cerevisiae strain, EMV-8, was obtained via mutation induced by ethyl methanesulfonate (EMS) and then adaptive evolution in the hydrolysate. EMV-8 exhibited a high vanillin consumption rate and a high total antioxidant capacity (Shen et al., 2014). Furthermore, through transcriptome analysis of strain EMV-8 and subsequent gene overexpression tests, several dehydrogenases and reductases (ZWF1, ALD6, YNL134C, and YJR096W) were revealed to function in vanillin detoxification by converting vanillin to vanillyl alcohol (Wang et al., 2016). The functions of reductases and dehydrogenases in vanillin reduction were divided into two groups: ➀ have vanillin reduction activities (Adh6p, Adh7p, and proteins encoded by YNL134C and YJR096W); ➁ NADPH supply for vanillin reduction (Ald6p and Zwf1p). However, overexpressing these genes did not produce as much vanillin tolerance as observed in EMV-8. In this study, we sequenced the genome of EMV-8 and its parental strain, NAN-27. Several mutations were found. Among them, the absence of Yrr1p, a transcription factor, was demonstrated to improve tolerance to vanillin in strains with different genetic backgrounds. Then, the transcriptional differences between yrr1Δ and wildtype S. cerevisiae were compared to reveal the underlying mechanism.

Materials and Methods

Medium and Culture Conditions

Yeast extract peptone dextrose (YPD) medium (10 g L^-1 yeast extract, 20 g L^-1 tryptone, 20 g L^-1 glucose) was used for activation, culture of host strains, and spot dilution growth (with 2% agar).

Synthetic complete medium (SD) or Sc-URA medium (1.7 g L^-1 yeast nitrogen base, Sangon, China, 5 g L^-1 ammonium sulfate, Sangon, China, CSM or CSM-URA, MP Biomedicals, Solon, OH, USA) was used for activation and batch fermentation of recombinant strains, supplying 20 g L^-1 glucose. Vanillin was used as an inhibitor and was added in the medium as indicated. All of the cultures were grown at 30°C unless otherwise indicated.

Strains and Plasmids

The S. cerevisiae industry strain NAN-27 (Zhang et al., 2010) and the vanillin-tolerant strain EMV-8, which was derived from NAN-27 (Shen et al., 2014), were the strains whose genomes were sequenced. Genetic manipulation including deletion and overexpression were mainly conducted in the lab strain BY4741 (MATa, his3Δ1 leu2Δ met5Δ ura3Δ, EUROSCARF, Germany). Deletion of YRR1 was also conducted in CEN.PK2-1C (MATa; ura3-52; trp1-289; leu2-3,112; his3Δ1; MAL2-8C; SUC2, EUROSCARF, Germany). The genes for overexpression were amplified from the mutant strain EMV-8. DBP2 was amplified from cDNA of BY4741. The gene deletion was performed using homologous recombination. The primers of this work are listed in Supplementary Table S1. The destruction cassette contained the KanMX expression cassette amplified from plasmid pUG6 (Escherichia coli plasmid with segment LoxP–KanMX4–LoxP) (Güldener et al., 1996) and the sequences that were homologous with the genes for knockout. The DNA fragments of the destruction cassettes were then transformed into the lab strain BY4741, and the mutants were screened in YPD medium with 800 mg L^-1 G418. pJFE3 (Shen et al., 2012) is a 2 μ plasmid with the TEF1 promoter, the PGK1 terminator, and URA3 as the selection marker. pJFE1 (developed in the laboratory) is a centromere plasmid with the TEF1 promoter, the PGK1 terminator, and the URA3 expression cassette as the selection marker.

Spot Dilution Growth Assay

A single colony was inoculated into 3 mL YPD medium and cultured for 12 h. The cultured cells were transferred into a fresh 10 mL of YPD medium with an initial OD₆₀₀ of 0.2. After overnight culture, the cells were harvested and washed three times with ddH₂O, followed by suspension in 1 mL of ddH₂O for 9 h to consume endogenous substances. The suspended cells were normalized to an OD₆₀₀ of 1.0. Then, 4 μl of 10-fold serial dilutions were spotted onto YPD or YPD with different concentrations of vanillin and incubated at 30°C for 1 days (for YPD) or 2–3 days (with vanillin in YPD).

Fermentation

A single colony was cultured into 3 mL of SD or SC-URA for 24 h. Then, the cultures were transferred into 20 mL of fresh medium with OD₆₀₀ at 0.2 for overnight culture. Then, the overnight-cultured cells were inoculated into 100-mL flasks with 40 mL of fermentation medium with an initial OD₆₀₀ of 0.2. Then, the fermentation was performed at 30°C and 200 rpm.

Calculation of Physiological Parameters

The biomass concentration was calculated with reference to the measured OD₆₀₀-dry weight correlation. For BY4741, one unit of OD₆₀₀ (detected by a BioPhotometer plus, Eppendorf, Germany) equals 0.18 g L^-1 of biomass. The maximum growth rates were estimated as the linear regression coefficients of the ln (OD₆₀₀) versus time during the exponential growth phase. Specific consumption rates of vanillin were calculated using the following equation:

$\mathit{r}\mathrm{=}\frac{{\mathit{A}}_{\mathit{n}}\mathit{-}{\mathit{A}}_{\mathit{m}}}{\frac{\mathrm{1}}{\mathrm{2}}{\mathrm{\Sigma }}_{\mathit{i}\mathit{=}\mathit{m}\mathrm{+}\mathrm{1}}^{\mathit{n}}\mathrm{(}{\mathit{B}}_{\mathit{i}}\mathrm{+}{\mathit{B}}_{\mathit{i}\mathit{-}\mathrm{1}}\mathrm{)}\mathrm{\times }\mathrm{(}{\mathit{t}}_{\mathit{i}}\mathit{-}{\mathit{t}}_{\mathit{i}\mathit{-}\mathrm{1}}\mathrm{)}}$

where r is the specific consumption rate during the phase from sampling point m to sampling point n; A, B, and t are the metabolite concentration, biomass concentration, and time, respectively, at sampling points n, i, and m, as previously described (Peng et al., 2012).

Analyses of Extracellular Fermentation Products

The concentrations of glucose and ethanol were measured using a Prominence LC-20A HPLC system (Shimadzu, Japan) furnished with an Aminex HPX-87H ion exchange column (Bio-Rad, USA) and refractive index detector RID-10A (Shimadzu, Japan), using 5 mmol L^-1 H₂SO₄ as the mobile phase with a flow rate of 0.6 mL/min at 45°C (Eliasson et al., 2000). A BioSil-C18 column (Bio-Rad, USA) with a mobile phase of 40% aqueous methanol at a flow rate of 0.6 mL/min at room temperature was used to test the concentrations of vanillin and vanillyl alcohol by HPLC. The peaks were detected via ultraviolet detection (SPD-M20A) at 210 nm (Ji et al., 2011).

Genome Sequencing

Genomic DNA was extracted from EMV-8 and NAN-27 using a yeast genome DNA extraction kit (Sangon Biotech, China) after yeast cells were harvested during the mid-log phase. Three paired-end genomic DNA libraries (200, 500, and 6000-bp) were constructed from NAN-27, and a 500-bp genomic DNA library was constructed from EMV-8. The library quality was assessed using the Qubit@ 2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 system. Finally, the library was sequenced on an Illumina HiSeq platform, and 125-bp paired-end reads were generated (an Illumina Hiseq2500 system).

Transcriptome Analysis

The pre-cultured cells of BY4741 and BY4741(yrr1Δ) in SD medium were transferred into fresh SD or SD containing 5 mM vanillin with an initial OD₆₀₀ of 0.2. The cells were harvested during the log phase (OD₆₀₀≈1.6-2.0, vanillin content in the medium approximately 0.73∼0.78 g/L) and were quick frozen using liquid nitrogen. The UNIQ-10 Trizol RNA Purification Kit (Sangon Biotech, China) was used to extract total RNA. mRNA was isolated, fragmented, and then used as templates to synthesize cDNA. The short fragments were connected with adapters to obtain suitable fragments for PCR amplification. Ultimately, the libraries were sequenced by the Illumina HiSeq^TM 4000 (BGI Shenzhen, China). The differentially expressed genes were screened out according to the following criteria: Fold change ≥ 2 and FDR ≤ 0.001. All of the analyses were performed in biological duplicate. The transcriptome data were deposited in the NCBI Gene Expression Omnibus database (GEO accession number: GSE89854). Gene functions were annotated using the tool on Saccharomyces genome database (SGD)¹.

Enzyme Activity Assay

The strains were cultured in SD medium. The cells were collected when the OD₆₀₀ reached 4.0. Then, the cells were disrupted using a Fast Prep cell homogenizer (Thermo Savant, Germany) with 33 mM Na₃PO₄ buffer (pH = 7.0) containing 0.5-mm glass beads and 1 mM PMSF, which was used to inhibit any protease activity. The supernatant was collected by centrifugation at 4°C for the enzyme activity assay (Larroy et al., 2002).

The vanillin dehydrogenase activity was measured using the assay described by Larroy (Larroy et al., 2002). The 0.6-mL reaction mixtures used to measure enzyme activities toward vanillin contained 33 mM sodium phosphate buffer (pH 7.0), 0.5 mM NADPH, and 1 mM vanillin (using cuvettes with 0.2-cm path length). The molar absorption coefficient (ε₃₆₅) was 7.71 mM^-1 cm^-1 for vanillin plus NADPH. One unit (U) of enzyme activity is defined as the amount of enzyme that can reduce 1 μmol of NADPH plus vanillin in 1 min. The BCA protein assay reagent kit (Beyotime, China) was used to measure protein concentrations of crude enzyme solution. The specific enzyme activity (U mg^-1 protein) refers to the enzyme activity per milligram of protein.

Results

Identification of the Mutations Affecting Vanillin Tolerance by Genome Sequencing and Growth Assay Screening

To identify the key genetic targets related to vanillin resistance, comparative genome analysis was performed between the vanillin-tolerant strain EMV-8 and its parent strain, NAN-27. The genome sequencing results identified a total of 450 non-synonymous SNPs in CDS and 44 CDS with InDels in strain EMV-8 compared to its parent strain, NAN-27. In the present work, we focused on genes with InDels. Using PCR amplification and sequence analysis, only seven genes with InDels were confirmed in the evolved strain (Table 1).

TABLE 1

TABLE 1. Summary of InDels in the mutant strain EMV-8 compared with its parent strain, NAN-27.

To determine whether these InDel-containing genes influence vanillin tolerance in the mutant strain, these seven genes were separately knocked out in the lab strain BY4741. The growth profiles of the knock-out strains for these genes were tested using a spot dilution growth assay on media with or without vanillin (Figure 1).

FIGURE 1

FIGURE 1. Vanillin resistance profiles of BY4741 and its single-gene mutants. The deletion mutants were derived from BY4741 and are shown with their genotypes. Four microliters of each 10-fold dilution with an initial cell OD₆₀₀ = 1 was spotted on YPD (for 1 days) or YPD containing 4 or 5 mM vanillin (for 2–3 days).

Nonsense Mutation of YRR1 Increased the Vanillin Tolerance of S. cerevisiae, but Deleting its Paralog, YRM1, Did Not

Deleting YRR1, as by the addition of a stop codon at 409 bp that shortened the final protein product from 811 to 141 aa in EMV-8, notably enhanced the strain growth on the plates containing vanillin (Figure 1). Meanwhile, deleting other genes did not improve the vanillin tolerance of the strain. The maximum specific growth rates of BY4741(yrr1Δ) were, respectively, 23 and 142% faster than those of BY4741 with 5 and 6 mM vanillin in SD medium (Figure 2 and Table 3).

FIGURE 2

FIGURE 2. Vanillin-resistance profiles of yrr1Δ and its parent strain in liquid medium. The curves of growth (A) and vanillin consumption (B) when 5 or 6 mM vanillin was included in SD medium. Symbols: square, BY4741; triangle, BY4741(yrr1Δ); solid, SD with 5 mM vanillin; open, SD with 6 mM vanillin. (C) The growth curve in SD medium. Symbols: square, BY4741; circle, BY4741(yrr1Δ). The error bars indicate the standard deviation of independent triplicate experiments.

To confirm that deletion of YRR1 was the decisive factor in vanillin detoxification, compensation of YRR1 was performed. YRR1 was inserted into the centromere plasmid pJFE1 under the control of the TEF1 promoter and the PGK1 terminator, producing the plasmid pJFE1-YRR1. Then, pJFE1-YRR1 was transformed into the BY4741(yrr1Δ), resulting in the YRR1-compensated strain BY4741(yrr1Δ, pJFE1-YRR1). Strains BY4741(yrr1Δ, pJFE1) and BY4741(pJFE1) were also constructed and used as controls. The growth test result showed that the lag phase was much longer in BY4741(yrr1Δ, pJFE1-YRR1) and BY4741(pJFE1), which expressed YRR1, than in BY4741 (yrr1Δ, pJFE1). The maximum growth rates of BY4741(pJFE1), BY4741(yrr1Δ, pJFE1-YRR1), and BY4741(yrr1Δ, pJFE1) were 0.0867 ± 0.0014, 0.0887 ± 0.0031, and 0.1214 ± 0.0012, respectively. The compensation of YRR1 prolonged the lag phase and the strain growth rate decreased by 27% in liquid Sc-ura with 6 mM vanillin (Figure 3).

FIGURE 3

FIGURE 3. Vanillin tolerance profiles of BY4741(pJFE1) (), BY4741(yrr1Δ, pJFE1-YRR1) (), and BY4741(yrr1Δ, pJFE1) () in liquid Sc-ura medium containing 6 mM vanillin. Data are shown as the average value ± standard deviation of the biological triplicates.

To investigate whether the deletion of YRR1 produced a general effect on vanillin tolerance, the strain CEN.PK2-1C (yrr1Δ), which has a genetic background different from that of BY4741, was constructed. This strain also exhibited obviously higher vanillin tolerance compared to the parent strain CEN.PK2-1C (Figure 4A). This result suggested that YRR1 deletion may be an effective method for enhancing the vanillin tolerance of S. cerevisiae.

FIGURE 4

FIGURE 4. Vanillin resistance profiles of (A) BY4741 and BY4741(yrm1Δ) and (B) CEN.PK2-1C and CEN.PK2-1C(yrr1Δ). Four microliters of each 10-fold dilution with an initial cell OD₆₀₀ = 1 was spotted on YPD (for 1 days) or YPD containing vanillin (for 2–3 days).

Yrr1p is a transcription factor containing a Zn₂-Cys₆ zinc-finger motif. YRR1 is a paralog of YRM1, which also contains a Zn₂-Cys₆ zinc-finger motif. To investigate whether Yrm1p has a function in vanillin tolerance similar to that of Yrr1p, the effect of YRM1 deletion was also studied. The deletion of YRM1 did not improve vanillin resistance (Figure 4B), demonstrating that this protein has a different function from Yrr1p in this context.

Transcriptome Analysis of BY4741 and BY4741(yrr1Δ) with and without Vanillin

To investigate the mechanism by which Yrr1p influences vanillin tolerance, the transcriptomic differences between BY4741 and BY4741(yrr1Δ) were analyzed in SD with 5 mM vanillin or without vanillin. The genes showing significant differential regulation under different conditions are listed in Table 2. Genes that are underlined are down-regulated; the others are up-regulated. The genes for which 2 ≥ fold change ≥ 1.5 are marked with “^∗”. Other genes listed in the table are those for which fold change ≥ 2.

TABLE 2

TABLE 2. The vanillin response genes and the genes significantly regulated by Yrr1p with and without vanillin stress.

TABLE 3

TABLE 3. Maximum specific growth rate and specific consumption rate for vanillin in yrr1Δ and BY4741.

When the cells were cultured in SD medium with no vanillin stress, the transcriptome differences between BY4741 and BY4741(yrr1Δ) were very small. Only eight genes showed obviously different expression (Figure 5), and these genes do not share any functions. Next, we analyzed how the strains responded to vanillin stress. BY4741 and BY4741(yrr1Δ) gave similar response, for each strain, over 200 genes showed differential expression when vanillin was present compared to vanillin free condition. When vanillin was present, dehydrogenases and ATP-binding cassette (ABC) transporters were up-regulated in both strains. This may be a general mechanism for resistance to vanillin stress in S. cerevisiae. Interestingly, genes involved in ribosome biogenesis and rRNA processing were significantly down-regulated in BY4741, while these genes and other related genes were significantly up-regulated in BY4741(yrr1Δ). This result suggested that the vanillin may decrease cell growth by repressing ribosome biogenesis and rRNA processing, thereby decreasing protein synthesis, and deleting YRR1 may restore these processes. The marked up-regulation of ribosome biogenesis and rRNA processing in BY4741(yrr1Δ) compared with BY4741 in the presence of 5 mM vanillin in the culture medium also supports this model. Based on these findings, more experiments were performed to further investigate how Yrr1p affects the vanillin tolerance of S. cerevisiae.

FIGURE 5

FIGURE 5. Crude enzyme activities using vanillin as the substrate. Error bars indicate the standard deviation of three independent replicates.

ADH7 Expression Was Induced by Vanillin, Especially in the Absence of YRR1

Adh7p, a NADPH-dependent medium-chain alcohol dehydrogenase, has been reported to have reductase activity and to convert vanillin to vanillyl alcohol (Nguyen et al., 2015). In both BY4741 and BY4741(yrr1Δ), the expression of ADH7 was up-regulated in the presence of vanillin, which was consistent with the report that ADH7 was inducible under severe vanillin stress (Nguyen et al., 2015). However, in BY4741(yrr1Δ), the increase in ADH7 was eight times that in BY4741. The deletion of YRR1 significantly increased the specific consumption rate of vanillin by 39 and 51% in SD with 5 and 6 mM vanillin, respectively (Figure 2). In BY4741(yrr1Δ), 6 mM vanillin was completely reduced within 26 h, while at that time point, only ∼10% of the vanillin was reduced by BY4741. Furthermore, in the enzyme activity assay, the cell-free extract of the BY4741(yrr1Δ) exhibited 95% higher NADPH-dependent reductase activity toward vanillin than did BY4741, while the cell-free extract of the BY4741(yrr1Δ) exhibited no NADH-dependent reductase activity compared to the control (Figure 5). Together, these results suggested that the up-regulation of ADH7 might be one reason why deleting YRR1 enhanced the strain tolerance to vanillin. To verify whether this is the only reason, we constructed the BY4741(adh7Δyrr1Δ). Interestingly, the adh7Δyrr1Δ strain still grew better than BY4741 in the presence of vanillin (Figure 6). Therefore, there are other mechanisms related to vanillin tolerance in the BY4741(yrr1Δ).

FIGURE 6

FIGURE 6. Vanillin resistance profiles of BY4741, yrr1Δ, and yrr1Δadh7Δ. Each test was repeated three times.

Deleting YRR1 Increased the Vanillin Tolerance of the Strain via the Increased Expression of Multidrug Transporters

The multidrug resistance proteins Pdr5p, Yor1p, and Snq2p occur in the plasma membrane and belong to the ABC transporter superfamily, members of which export a diverse range of intracellular compounds using the energy supplied by ATP hydrolysis (Ananthaswamy et al., 2010; Guo et al., 2012). The transcriptional levels of PDR5, YOR1, and SNQ2 in both BY4741 and BY4741(yrr1Δ) were significantly up-regulated under vanillin stress, and the degree of up-regulation of PDR5 in BY4741(yrr1Δ) was notably greater than that in BY4741. Moreover, in our previous transcriptome analysis comparing EMV-8 with NAN-27, several multidrug transporters were expressed at markedly higher levels in EMV-8 than in NAN-27, including Pdr10p (6.2 times) and Yor1p (3.6 times). Next, the abovementioned genes related to multidrug resistance were overexpressed in the lab strain BY4741. In Sc-ura medium, the strain overexpressing PDR5 and PDR10 grew much more slowly than the control BY4741(pJFE3) (Figure 7C). Nonetheless, when vanillin was present, the strain overexpressing PDR5, YOR1, and SNQ2 grew markedly faster than control (Figures 7A,B) and had a shorter lag phase which was defined as the time before OD₆₀₀ reached 1.0 (Table 4).

TABLE 4

TABLE 4. Maximum specific growth rate and specific consumption rate for vanillin of recombinant strains.

FIGURE 7

FIGURE 7. The characteristics of strains overexpressing PDR5 (); PDR10 (); SNQ2 (); and YOR1 (), as well as the control BY4741(pJFE3) () in Sc-ura with 6 mM vanillin (A,B) and in Sc-ura without vanillin (C). Data are shown as the mean ± standard deviation of independent triplicate experiments.

Deleting YRR1 Increased the Vanillin Tolerance of the Strain via the Enhancement of Ribosome Biogenesis and rRNA Processing

Ribosomes are responsible for protein synthesis in the cell (Pelava et al., 2016). Reduced or abnormal ribosome production occurs in response to cellular stress and leads to decreased growth and metabolic capacity of the strain (Nomura, 1999). Our transcriptional analysis result indicated that ribosome biogenesis and rRNA processing were repressed by vanillin in wildtype S. cerevisiae. Deleting YRR1 increased the expression of genes related to ribosome biogenesis and rRNA processing when vanillin was present. Nine RNA helicases, including DBP2, DBP3, DBP9, DBP10, PRP43, MTR4, HAS1, DHH1, and DRS1. It was reported that the proteins encoding by these genes are participated in ribosome biogenesis, rRNA processing, and translation initiation (Rodriguez-Galan et al., 2013), were identified from dozens of markedly up-regulated genes in the comparison between BY4741 and BY4741(yrr1Δ). To further confirm whether enhancing ribosome biogenesis and rRNA processing can increase the vanillin tolerance of a strain, the gene DBP2, which encodes an ATP-dependent RNA helicase that functions in mRNP assembly (Ma et al., 2013) and is the gene most significantly down-regulated in BY4741 under vanillin stress but up-regulated in yrr1Δ, was overexpressed in BY4741. Although the maximum specific growth rate of BY4741(DBP2) was decreased by 15% compared to the control under vanillin-free conditions, in the presence of vanillin, the growth and metabolic characteristics of BY4741(DBP2) were better than those of the control. In SD medium containing 6 mM vanillin, the maximum specific growth rate of strain overexpressing DBP2 was 7% higher than that of the control. The OD₆₀₀ of the DBP2-overexpressing strain reached 1.0 and 2.0 after only 24 and 30 h, respectively, while reach the same biomass in the control required 30 and 40 h, respectively. The specific vanillin consumption rate of the strain overexpressing DBP2 was 0.062 ± 0.002 g g^-1 h^-1, which was 59% higher than that of the control (0.039 ± 0.001 g g^-1 h^-1, Figure 8).

FIGURE 8

FIGURE 8. The characteristics of a strain overexpressing DBP2 () and the control strain BY4741(pJFE1) () in Sc-ura containing 6 mM vanillin (A,B) and in Sc-ura without vanillin (C). Data are shown as the mean ± standard deviation of the biological triplicates.

Discussion

Yrr1p, a Zn₂-Cys₆ zinc-finger transcription factor, has been reported to be an activator of genes related to multidrug resistance (Zhang et al., 2001; Le Crom et al., 2002; Teixeira et al., 2008; Kodo et al., 2013). It has been reported that Yrr1p contributes to tolerance of toxicants such as 4-nitroquinoline-N-oxide (4-NQO) (Cui et al., 1998). In the present work, we found that the absence of Yrr1p provided yeast strains with vanillin resistance characteristics.

The transcription result showed that only eight genes (Table 5) were significantly differentially regulated after YRR1 deletion when the cells were cultured under vanillin-free conditions. Furthermore, changing the expression levels of these genes had no effect on the vanillin tolerance of strain (data not shown). However, the situation was quite different when strains were cultured under vanillin stress. When vanillin was present, the absence of Yrr1p led to a marked effect on the transcriptome, with at least three related functions: vanillin detoxification, excretion of vanillin or other toxic molecules generated because of vanillin stress, and the promotion of the biological processes of rRNA processing and ribosome biogenesis. In these up-regulated genes, the dehydrogenase gene ADH7 was the most up-regulated in comparison of BY4741 and BY4741(yrr1Δ) when vanillin existed. Adh7p and its homologous gene Adh6p were considered to reduce vanillin to vanillyl alcohol (Liu, 2011; Wang et al., 2016). However, the deletion of ADH7 in BY4741(yrr1Δ) didn’t affect its tolerance to vanillin. Hence, ADH7 was not the main reason why deleting YRR1 enhanced the strain tolerance to vanillin. Genes for dehydrogenases and ABC transporters were up-regulated under vanillin stress. When Yrr1p was absent, their expression levels were even higher. Genes involved in ribosome biogenesis and rRNA processing, which are important for protein synthesis, were down-regulated under vanillin stress but up-regulated when Yrr1p was absent. These phenomena suggested that vanillin might be an effector molecule in this context. Vanillin may induce the expression of some genes related to rRNA processing and ribosome biogenesis. However, these genes are directly or indirectly repressed by Yrr1p, only when Yrr1p is absent, their expression are activated by vanillin. The regulatory mechanism of Yrr1p under vanillin stress remain to be discovered.

TABLE 5

TABLE 5. The genes significantly regulated after deletion of YRR1 without vanillin stress.

ATP-binding cassette transporters have been considered to play a key role in increasing strains’ tolerance to chemicals by exporting inhibitors (Liu, 2011). In the present study, overexpression of PDR5, YOR1, and SNQ2 was demonstrated to improve vanillin tolerance. However, overexpressing these large proteins also burdened the cells. The strain overexpressing these genes grow more slowly than wildtype under vanillin-free conditions or after vanillin was depleted. Therefore, their expression should be controlled at a suitable level.

When cells are stimulated by harmful conditions, such as starvation, heat shock, oxidative stress, or acidic stress, translation is depressed (Nomura, 1999; Shenton et al., 2006; Hoyle et al., 2007; Iwaki and Izawa, 2012; Iwaki et al., 2013a; Dever et al., 2016). It has been reported that under the vanillin stress, messenger ribonucleoprotein (mRNP) granule formation (P-bodies and stress granules) may form because the initiation of translation was inhibited and the polysomes were disassembled (Iwaki et al., 2013b). Our transcription results also confirmed that ribosome biogenesis and rRNA processing, which are considered the source of all protein synthesis (Pelava et al., 2016), were significantly inhibited by vanillin. Furthermore, up-regulating ribosome biogenesis and rRNA processing, whether by deleting YRR1 or overexpressing DBP2, improves cell growth under vanillin stress. Only overexpressing single or several genes did not elevate the protein synthesis efficiently, whereas deleting YRR1 holistically enhanced the ribosome biogenesis and rRNA processing.

Conclusion

We compared the genome sequence of a vanillin-tolerant strain and its parent. Then, we found that the loss of Yrr1p is an efficient strategy for enhancing the vanillin tolerance of S. cerevisiae, likely by improving the strain’s capacity for vanillin detoxification, vanillin excretion, rRNA processing and ribosome biogenesis. Furthermore, the up-regulation of rRNA processing and ribosome biogenesis occurred only in the presence of vanillin, suggesting that vanillin might be an effector molecule in this regulatory process. To our knowledge, this is the first report of a relationship between YRR1 and vanillin tolerance and the first report of enhancing S. cerevisiae vanillin tolerance via up-regulating rRNA processing and ribosome biogenesis. These new findings will provide novel targets for constructing strains tolerant of lignocellulosic materials and valuable insights into the novel functions of YRR1 in S. cerevisiae.

Author Contributions

XW, YS, and XB conceived and designed the study; XW and ZL participated in the design of experiments and data collection; XW performed genome sequencing and RNA-seq data analysis; XW and YS analyzed data and drafted the manuscript. YS, XB, and JH supervised and coordinated the overall study. All authors read and approved the final manuscript.

Funding

This work was supported by grants from the National High Technology Research and Development Program of China (2014AA021903), the National Key Technology R&D Program of China (2014BAD02B07), the Project of the National Energy Administration of China (NY20130402), the State Key Laboratory of Motor Vehicle Biofuel Technology (KFKT2013002, 2013004), and the Shandong Key Laboratory of Straw Biorefinement Technologies.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Abbreviations

NADH, Reduced form of nicotinamide adenine dinucleotide; NADPH, Reduced form of nicotinamide adenine dinucleotide phosphate; ROS, reactive oxygen species.

Footnotes

1. ^http://sgd-beta.stanford.edu/

References

Ananthaswamy, N., Rutledge, R., Sauna, Z. E., Ambudkar, S. V., Dine, E., Nelson, E., et al. (2010). The signaling interface of the yeast multidrug transporter Pdr5 adopts a cis conformation, and there are functional overlap and equivalence of the deviant and canonical Q-loop residues. Biochemistry 49, 4440–4449. doi: 10.1021/bi100394j

PubMed Abstract | CrossRef Full Text | Google Scholar

Ask, M., Mapelli, V., Höck, H., Olsson, L., and Bettiga, M. (2013). Engineering glutathione biosynthesis of Saccharomyces cerevisiae increases robustness to inhibitors in pretreated lignocellulosic materials. Microb. Cell Fact. 12:87. doi: 10.1186/1475-2859-12-87

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, Z., Shiraki, T., Hirata, D., and Miyakawa, T. (1998). Yeast gene YRR1, which is required for resistance to 4-nitroquinoline N-oxide, mediates transcriptional activation of the multidrug resistance transporter gene SNQ2. Mol. Microbiol. 29, 1307–1315. doi: 10.1046/j.1365-2958.1998.01027.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dever, T. E., Kinzy, T. G., and Pavitt, G. D. (2016). Mechanism and regulation of protein synthesis in Saccharomyces cerevisiae. Genetics 203, 65–107. doi: 10.1534/genetics.115.186221

PubMed Abstract | CrossRef Full Text | Google Scholar

Eliasson, A., Christensson, C., Wahlbom, C. F., and Hahn-Hägerdal, B. (2000). Anaerobic xylose fermentation by recombinant Saccharomyces cerevisiae carrying XYL1, XYL2, and XKS1 in mineral medium chemostat cultures. Appl. Environ. Microbiol. 66, 3381–3386. doi: 10.1128/AEM.66.8.3381-3386.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Endo, A., Nakamura, T., and Shima, J. (2009). Involvement of ergosterol in tolerance to vanillin, a potential inhibitor of bioethanol fermentation, in Saccharomyces cerevisiae. FEMS Microbiol. Lett. 299, 95–99. doi: 10.1111/j.1574-6968.2009.01733.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Güldener, U., Heck, S., Fielder, T., Beinhauer, J., and Hegemann, J. H. (1996). A new efficient gene disruption cassette for repeated use in budding yeast. Nucleic Acids Res. 24, 2519–2524. doi: 10.1093/nar/24.13.2519

CrossRef Full Text | Google Scholar

Guo, X., Li, J., Wang, T., Liu, Z., Chen, X., Li, Y., et al. (2012). A mutation in intracellular loop 4 affects the drug-efflux activity of the yeast multidrug resistance ABC transporter Pdr5p. PLoS ONE 7:e29520. doi: 10.1371/journal.pone.0029520

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoyle, N. P., Castelli, L. M., Campbell, S. G., Holmes, L. E. A., and Ashe, M. P. (2007). Stress-dependent relocalization of translationally primed mRNPs to cytoplasmic granules that are kinetically and spatially distinct from P-bodies. J. Cell Biol. 179, 65–74. doi: 10.1083/jcb.200707010

PubMed Abstract | CrossRef Full Text | Google Scholar

Iwaki, A., and Izawa, S. (2012). Acidic stress induces the formation of P-bodies, but not stress granules, with mild attenuation of bulk translation in Saccharomyces cerevisiae. Biochem. J. 446, 225–233. doi: 10.1042/Bj20120583

PubMed Abstract | CrossRef Full Text | Google Scholar

Iwaki, A., Kawai, T., Yamamoto, Y., and Izawa, S. (2013a). Biomass conversion inhibitors furfural and 5-hydroxymethylfurfural induce formation of messenger RNP granules and attenuate translation activity in Saccharomyces cerevisiae. Appl. Environ. Microbiol. 79, 1661–1667. doi: 10.1128/Aem.02797-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Iwaki, A., Ohnuki, S., Suga, Y., Izawa, S., and Ohya, Y. (2013b). Vanillin inhibits translation and induces messenger ribonucleoprotein (mRNP) granule formation in Saccharomyces cerevisiae: application and validation of high-content, image-based profiling. PLoS ONE 8:e61748. doi: 10.1371/journal.pone.0061748

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, L., Shen, Y., Xu, L. L., Peng, B. Y., Xiao, Y. Z., and Bao, X. M. (2011). Enhanced resistance of Saccharomyces cerevisiae to vanillin by expression of lacA from Trametes sp AH28-2. Bioresour. Technol. 102, 8105–8109. doi: 10.1016/j.biortech.2011.06.057

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S. K., Jin, Y. S., Choi, I. G., Park, Y. C., and Seo, J. H. (2015). Enhanced tolerance of Saccharomyces cerevisiae to multiple lignocellulose-derived inhibitors through modulation of spermidine contents. Metab. Eng. 29, 46–55. doi: 10.1016/j.ymben.2015.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Klinke, H. B., Thomsen, A. B., and Ahring, B. K. (2004). Inhibition of ethanol-producing yeast and bacteria by degradation products produced during pre-treatment of biomass. Appl. Microbiol. Biotechnol. 66, 10–26. doi: 10.1007/s00253-004-1642-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kodo, N., Matsuda, T., Doi, S., and Munakata, H. (2013). Salicylic acid resistance is conferred by a novel YRR1 mutation in Saccharomyces cerevisiae. Biochem. Biophys. Res. Commun. 434, 42–47. doi: 10.1016/j.bbrc.2013.03.069

PubMed Abstract | CrossRef Full Text | Google Scholar

Larroy, C., Fernández, M. R., González, E., Parés, X., and Biosca, J. A. (2002). Characterization of the Saccharomyces cerevisiae YMR318C (ADH6) gene product as a broad specificity NADPH-dependent alcohol dehydrogenase: relevance in aldehyde reduction. Biochem. J. 361, 163–172. doi: 10.1042/0264-6021:3610163

PubMed Abstract | CrossRef Full Text | Google Scholar

Le Crom, S., Devaux, F., Marc, P., Zhang, X., Moye-Rowley, W. S., and Jacq, C. (2002). New insights into the pleiotropic drug resistance network from genome-wide characterization of the YRR1 transcription factor regulation system. Mol. Cell. Biol. 22, 2642–2649. doi: 10.1128/MCB.22.8.2642-2649.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z. L. (2011). Molecular mechanisms of yeast tolerance and in situ detoxification of lignocellulose hydrolysates. Appl. Microbiol. Biotechnol. 90, 809–825. doi: 10.1007/s00253-011-3167-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z. L., Ma, M. G., and Song, M. Z. (2009). Evolutionarily engineered ethanologenic yeast detoxifies lignocellulosic biomass conversion inhibitors by reprogrammed pathways. Mol. Genet. Genomics 282, 233–244. doi: 10.1007/s00438-009-0461-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, W. K., Cloutier, S. C., and Tran, E. J. (2013). The DEAD-box protein Dbp2 functions with the RNA-binding protein Yra1 to promote mRNP assembly. J. Mol. Biol. 425, 3824–3838. doi: 10.1016/j.jmb.2013.05.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, T. T. M., Iwaki, A., and Izawa, S. (2015). The ADH7 Promoter of Saccharomyces cerevisiae is vanillin-inducible and enables mRNA translation under severe vanillin stress. Front. Microbiol. 6:1390. doi: 10.3389/fmicb.2015.01390

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, T. T. M., Iwaki, A., Ohya, Y., and Izawa, S. (2014a). Vanillin causes the activation of Yap1 and mitochondrial fragmentation in Saccharomyces cerevisiae. J. Biosci. Bioeng. 117, 33–38. doi: 10.1016/j.jbiosc.2013.06.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, T. T. M., Kitajima, S., and Izawa, S. (2014b). Importance of glucose-6-phosphate dehydrogenase (G6PDH) for vanillin tolerance in Saccharomyces cerevisiae. J. Biosci. Bioeng. 118, 263–269. doi: 10.1016/j.jbiosc.2014.02.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Nomura, M. (1999). Regulation of ribosome biosynthesis in Escherichia coli and Saccharomyces cerevisiae: diversity and common principles. J. Bacteriol. 181, 6857–6864.

PubMed Abstract | Google Scholar

Pelava, A., Schneider, C., and Watkins, N. J. (2016). The importance of ribosome production, and the 5S RNP–MDM2 pathway, in health and disease. Biochem. Soc. Trans. 44, 1086–1090. doi: 10.1042/BST20160106

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, B., Shen, Y., Li, X., Chen, X., Hou, J., and Bao, X. (2012). Improvement of xylose fermentation in respiratory-deficient xylose-fermenting Saccharomyces cerevisiae. Metab. Eng. 14, 9–18. doi: 10.1016/j.ymben.2011.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodriguez-Galan, O., Garcia-Gomez, J. J., and de la Cruz, J. (2013). Yeast and human RNA helicases involved in ribosome biogenesis: current status and perspectives. Biochim. Biophys. Acta 1829, 775–790. doi: 10.1016/j.bbagrm.2013.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Y., Chen, X., Peng, B. Y., Chen, L. Y., Hou, J., and Bao, X. M. (2012). An efficient xylose-fermenting recombinant Saccharomyces cerevisiae strain obtained through adaptive evolution and its global transcription profile. Appl. Microbiol. Biotechnol. 96, 1079–1091. doi: 10.1007/s00253-012-4418-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Y., Li, H., Wang, X., Zhang, X., Hou, J., Wang, L., et al. (2014). High vanillin tolerance of an evolved Saccharomyces cerevisiae strain owing to its enhanced vanillin reduction and antioxidative capacity. J. Ind. Microbiol. Biotechnol. 41, 1637–1645. doi: 10.1007/s10295-014-1515-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Shenton, D., Smirnova, J. B., Selley, J. N., Carroll, K., Hubbard, S. J., Pavitt, G. D., et al. (2006). Global translational responses to oxidative stress impact upon multiple levels of protein synthesis. J. Biol. Chem. 281, 29011–29021. doi: 10.1074/jbc.M601545200

PubMed Abstract | CrossRef Full Text | Google Scholar

Teixeira, M. C., Dias, P. J., Simoes, T., and Sa-Correia, I. (2008). Yeast adaptation to mancozeb involves the up-regulation of FLR1 under the coordinate control of Yap1, Rpn4, Pdr3, and Yrr1. Biochem. Biophys. Res. Commun. 367, 249–255. doi: 10.1016/j.bbrc.2007.12.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Liang, Z., Hou, J., Bao, X., and Shen, Y. (2016). Identification and functional evaluation of the reductases and dehydrogenases from Saccharomyces cerevisiae involved in vanillin resistance. BMC Biotechnol. 16:31. doi: 10.1186/s12896-016-0264-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Cui, Z., Miyakawa, T., and Moye-Rowley, W. S. (2001). Cross-talk between transcriptional regulators of multidrug resistance in Saccharomyces cerevisiae. J. Biol. Chem. 276, 8812–8819. doi: 10.1074/jbc.M010686200

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Shen, Y., Shi, W., and Bao, X. (2010). Ethanolic cofermentation with glucose and xylose by the recombinant industrial strain Saccharomyces cerevisiae NAN-127 and the effect of furfural on xylitol production. Bioresour. Technol. 101, 7104–7110. doi: 10.1016/j.biortech.2010.03.129

PubMed Abstract | CrossRef Full Text | Google Scholar