Skip to main content

Exploring the stress response mechanisms to 2-phenylethanol conferred by Pdr1p mutation in Saccharomyces cerevisiae

Abstract

Background

The 2-phenylethanol (2-PE) tolerance phenotype is crucial to the production of 2-PE, and Pdr1p mutation can significantly increase the tolerance of 2-PE in Saccharomyces cerevisiae. However, its underlying molecular mechanisms are still unclear, hindering the rational design of superior 2-PE tolerance performance.

Results

Here, the physiology and biochemistry of the PDR1_862 and 5D strains were analyzed. At 3.5 g/L 2-PE, the ethanol concentration of PDR1_862 decreased by 21%, and the 2-PE production of PDR1_862 increased by 16% than those of 5D strain. Transcriptome analysis showed that at 2-PE stress, Pdr1p mutation increased the expression of genes involved in the Ehrlich pathway. In addition, Pdr1p mutation attenuated sulfur metabolism and enhanced the one-carbon pool by folate to resist 2-PE stress. These metabolic pathways were closely associated with amino acids metabolism. Furthermore, at 3.5 g/L 2-PE, the free amino acids content of PDR1_862 decreased by 31% than that of 5D strain, among the free amino acids, cysteine was key amino acid for the enhancement of 2-PE stress tolerance conferred by Pdr1p mutation.

Conclusions

The above results indicated that Pdr1p mutation enhanced the Ehrlich pathway to improve 2-PE production of S. cerevisiae, and Pdr1p mutation altered the intracellular amino acids contents, in which cysteine might be a biomarker in response to Pdr1p mutation under 2-PE stress. The findings help to elucidate the molecular mechanisms for 2-PE stress tolerance by Pdr1p mutation in S. cerevisiae, identify key metabolic pathway responsible for 2-PE stress tolerance.

Background

2-Phenylethanol (2-PE) is an aromatic alcohol with a rose-like fragrance that is widely used in various fields such as food, medicine, and fragrance industries [1]. With increasing concern about climate change, it is of great significance to develop green, sustainable, and efficient technologies to produce 2-PE [2]. At present, Saccharomyces cerevisiae is one of the most promising microorganisms for the biosynthesis of 2-PE [3]. However, 2-PE has a strong inhibitory effect on S. cerevisiae as its structure [4]. It results in the maximum concentration of 2-PE below 6.3 g/L, which is far below the level required for industrial fermentation.

Pdr1p is a central multidrug-resistant transcription factor in Candida glabrata [5]. C. glabrata acquires azole resistance through gain-of-function (GOF) mutations in Pdr1p [6, 7]. GOF mutation usually is single amino acid substitution mutations in the gene encoding Pdr1p. These mutations yield a GOF phenotype and lead to the elevated transcription of downstream target genes [8]. GOF alleles of Pdr1p have been extensively studied in azole resistance of C. glabrata [9]. In addition, the mutant allele, Pdr1p (R821H), conferred resistance to high iron, copper and manganese medium in S. cerevisiae [10]. And Pdr1pF815S also showed increased resistance to organic solvent in S. cerevisiae [11]. Currently, it has been demonstrated that Pdr1p mutation (C862R) can substantially increase the 2-PE stress tolerance of S. cerevisiae [12]. But little is known of how Pdr1p is normally regulated [13].

Therefore, the present study investigated the growth, physiological, and biochemical parameters of Pdr1p mutation in S. cerevisiae under 2-PE stress. Furthermore, transcriptomic analysis was used to reveal the global gene response in response to Pdr1p mutation at 2-PE stress. The findings shed light on the molecular mechanisms for Pdr1p mutation in the 2-PE stress response, which would provide the target for constructing stress-resistant yeast cell factories.

Results

Changes in growth and metabolism conferred by Pdr1p mutation at 2-PE stress

To investigate the effect of the Pdr1p mutation on the growth and fermentation of S. cerevisiae under 2-PE stress, in this study, PDR1_862 strain and 5D strain were cultured under 3.5 g/L and 0 g/L 2-PE, respectively. As shown in Fig. 1A, no differences were observed in the growth of PDR1_862 and 5D strains at 0 g/L 2-PE. Under 3.5 g/L 2-PE stress, the OD600 of PDR1_862 and 5D strains was 21.25 and 13.8, respectively, which were 27% and 50% lower than the OD600 in the absence of 2-PE stress, respectively (Fig. 1B). Notably, 3.5 g/L 2-PE substantially inhibited the growth of the yeast strains. However, the maximum OD600 of strain PDR1_862 under 3.5 g/L 2-PE stress was 54% higher than that of 5D, and its growth rate was 11% higher than that of 5D (Table 1). Similarly, the maximum biomass dry weight of strain PDR1_862 under 3.5 g/L 2-PE stress (7.16 mg/mL) increased by 45.05% than that of the 5D strain. These results indicated that the Pdr1p mutation markedly enhances the growth of S. cerevisiae under 2-PE stress. However, no difference was observed in the glucose consumption between PDR1_862 and 5D strains under 2-PE stress. Despite this, some differences were noted in their ability to synthesize ethanol: under 3.5 g/L 2-PE stress, the maximum ethanol concentration and ethanol synthesis rate of the PDR1_862 mutant strain were 21% and 22% lower than those of the 5D strain, respectively. Therefore, the growth of the PDR1_862 strain increased while the rate of ethanol synthesis decreased under 2-PE stress.

Fig. 1
figure 1

The OD600 and biomass dry weight of PDR1_862 and 5D strains at 0 g/L (A) and 3.5 g/L 2-PE (B), and 2-PE concentration (C) of PDR1_862 and 5D strains. * Represents p < 0.05, as determined by Student’s t-test

Table 1 Effect of the mutant of Pdr1p on the fermentation parameters

Effect of Pdr1p mutation on 2-PE content of S. cerevisiae

Pdr1p mutation enhanced the adaptability of S. cerevisiae to 2-PE stress, which led us to question if this enhancement would help to improve the ability of S. cerevisiae to produce 2-PE from l-Phe through Ehrlich pathway. Therefore, this paper compared the differences in the ability of PDR1_862 mutant and 5D strains to synthesize 2-PE. As shown in Fig. 1C, there was no difference in 2-PE concentration between PDR1_862 and 5D during early fermentation (0–6 h). However, at 12 h, the 2-PE concentration of PDR1_862 increased by 10% than that of 5D strain when fermented for 24 h, the 2-PE concentration of PDR1_862 significantly increased by 17% than that of 5D strain. After 36 h, the 2-PE concentration produced by strain PDR1_862 increased to 2.22 g/L, which was 16% higher than that of the 5D strain (1.92 g/L). These results demonstrated that the increased resistance to 2-PE stress in strain PDR1_862 was conducive to the improvement of 2-PE concentration.

In summary, the Pdr1p mutation reduced ethanol synthesis rate of yeast cells under 2-PE stress but enhanced the growth and the ability of the yeast cells to synthesize 2-PE.

Global gene response at the transcriptional level in response to Pdr1p mutation at 2-PE stress

To further investigate the underlying molecular mechanisms of Pdr1p mutation in the 2-PE stress response, RNA-Seq was conducted to compare global gene expression in S. cerevisiae PDR1_862 and 5D strains at 0 g/L and 3.5 g/L 2-PE. As shown in Fig. 2A, the samples exhibited high intragroup repeatability and large intergroup differences. A total 1524 genes were differentially expressed in the four groups, and 1420 DEGs were identified in the 5D vs 5DPE (749 upregulated and 671 downregulated); 711 DEGs were identified in the PDR vs PDRPE (409 upregulated and 302 downregulated); 120 DEGs were found in the 5DPE vs PDRPE (35 upregulated and 85 downregulated); 7 genes were upregulated between the 5D and PDR group (Fig. 2B). Among the DEGs, 4 genes were identified in all four groups, and eight DEGs were present only in 5DPE vs PDRPE (Fig. 2C), with three genes (SNZ1, GCV1, and TNA1) significantly upregulated and five genes (SNQ2, CPA1, ROX1, MET32, and a gene with unknown function) significantly downregulated. These genes might be the targets for that Pdr1p mutation conferred to 2-PE stress. Among the eight genes, some genes, such as GCV1, MET32, CPA1, have been previously shown to be associated with the amino acid metabolism [14,15,16].

Fig. 2
figure 2

Global analysis of transcriptomes and DEGs. A PCA analysis. B Differential statistics between groups, red represents upregulation, green represents downregulation. C Venn analysis of differential genes in different groups

To further analyze pathways enriched by the DEGs, the total DEGs were subjected to KEGG pathway enrichment analysis. KEGG analysis revealed that the DEGs in the 5D vs 5DPE group were enriched in ribosomes, cyanoamino acid metabolism, purine metabolism, thiamine metabolism, fructose and mannose metabolism, nitrogen metabolism, amino sugar and nucleotide sugar metabolism, aminoacyl-tRNA biosynthesis, propionate metabolism, pyruvate metabolism, glycolysis/gluconeogenesis, unsaturated fatty acid biosynthesis, fatty acid elongation. However, Pdr1p mutation significantly altered the metabolic pathways in S. cerevisiae under 2-PE stress. The significant metabolic pathways enriched in 5DPE vs PDRPE were sulfur metabolism, folate one-carbon pool, yeast meiosis, glyoxylate and dicarboxylic acid metabolism, ABC transporter protein, glycine/serine/threonine metabolism (Fig. 3B). Among the significant enrichment pathways, sulfur metabolism, glyoxylate and dicarboxylic acid metabolism, and glycine/serine/threonine metabolism were closely related to amino acid metabolism (Fig. 3C). Therefore, the amino acid metabolism might subject to perturbation by Pdr1p mutation, and the intracellular amino acid concentrations would be changed. Of these metabolic pathways, sulfur metabolism displayed the most significant enrichment, which was inextricably bound to cysteine metabolism and glycine, serine, and threonine metabolism. Amino acid metabolism involved in alleviating stress [17, 18], such as acid stress, oxidative stress.

Fig. 3
figure 3

The KEGG enrichment analysis of the DEGs in 5DPE_vs_PDRPE (A) and 5DPE_vs_PDRPE (B), and the network of the enrichment pathways (C)

Suppression of sulfur metabolism conferred by Pdr1p mutation

Compared to those in the absence of 2-PE, the transcription levels of the ATP sulfurylase (MET3, log2FC =  + 2.48), APS kinase (MET14, log2FC =  + 3.40), PAPS reductase (MET16, log2FC =  + 2.94), and sulfite reductase (MET10, log2FC =  + 1.67) genes involved in sulfur metabolism in 5D strain were significantly upregulated (Fig. 4); for strain PDR1_862, there was no difference of the above genes between 0 and 3.5 g/L 2-PE. However, compared to those of 5D strain at 3.5 g/L 2-PE, the transcription levels of the ATP sulfurylase (MET3, log2FC = − 2.04), APS kinase (MET14, log2FC = − 2.41), PAPS reductase (MET16, log2FC = − 1.25), and sulfite reductase (MET10, log2FC = − 1.33) genes were significantly downregulated in PDR1_862 strain. In conclusion, 2-PE stress enhanced transcription of the genes related to sulfur metabolism pathway in S. cerevisiae, and this transcriptional enhancement was attenuated by the Pdr1p mutation.

Fig. 4
figure 4

Changes of differentially expressed genes in sulfur metabolism and one-carbon pool by folate

In addition, hydrogen sulfide (produced from sulfur metabolism) and O-acetylhomoserine were converted into cysteine by enzymes such as O-acetylserine sulfhydrylase (MET17) [19, 20]. 2-PE stress significantly upregulated the transcription of MET17 (log2FC = + 1.67) in strain 5D. However, compared with that of 5D strain, the transcription level of MET17 was significantly downregulated in the PDR1_862 strain under 2-PE stress. Therefore, under 2-PE stress, Pdr1p mutation might reduce the synthesis of l-cysteine, which is an important sulfur-containing amino acid in yeast cells.

Enhancement of one-carbon pool by folate conferred by Pdr1p mutation

Folate concentrations of yeast cells are limiting, which can mediate one-carbon metabolism such as serine glycine interconversion and de novo purine synthesis. Compared to those in the absence of 2-PE, the transcription levels of formyltransferase (ADE17, log2FC = − 162), glycine decarboxylase (GCV1, log2FC = − 0.16), glycine hydroxymethyltransferase (SHM2, log2FC = − 1.64), and methenyltetrahydrofolate cyclohydrolase (MIS1, log2FC = − 0.21), which are all involved in the folate one-carbon pool, were downregulated in 5D strain under 2-PE stress. However, the transcription levels of formyltransferase (ADE17, log2FC = + 1.03), glycine decarboxylase (GCV1, log2FC = + 1.24), and glycine hydroxymethyltransferase (SHM2, log2FC = + 1.30) were significantly upregulated in the PDR1_862 strain compared with those in the 5D strain under 2-PE stress. 2-PE inhibited the expression of genes in the folate one-carbon pool pathway in S. cerevisiae, the Pdr1p mutation alleviated the transcriptional suppression of these genes. Folate one-carbon pool pathway was associated with the interconversion between glycine and serine, which also related to other amino acid metabolism (Fig. 4).

Enhancement of the Ehrlich pathway conferred by Pdr1p mutation

Notably, 2-PE stress significantly reduced the transcription of aromatic amino acid aminotransferase (ARO9, log2FC = − 2.05) and aromatic amino acid decarboxylase (ARO10, log2FC = − 2.21) of the Ehrlich pathway in S. cerevisiae, whereas no significant change was found in the transcription of the above two genes in the PDR1_862 strain (Fig. 5). However, the transcription levels of ARO9 and ARO10 were significantly upregulated 2.39- and 2.30-fold, respectively, in 5DPE vs PDRPE. Thus, 2-PE stress repressed the transcription of the genes encoded in the Ehrlich pathway in S. cerevisiae, but the Pdr1p mutation increased the transcription level of the genes in this pathway. Pdr1p mutation enhance 2-PE synthesis through the Ehrlich pathway in S. cerevisiae.

Fig. 5
figure 5

Changes of differentially expressed genes in Ehrlich pathway

Validation of selected genes by qPCR

To validate RNA-seq data, 12 genes were selected to analyze the expression level by qPCR. As shown in Fig. 6, 2-PE stress increased the expression levels of MET3, MET10, MET14, MET16, and MET17 in S. cerevisiae, and the transcription levels of these genes was still upregulated by 2-PE stress after Pdr1p mutation. However, their expression was significantly downregulated compared with that of the 5D strain. 2-PE stress significantly downregulated the expressions of the ADE17, SHM2, GCV1, and GCV2 genes. Although 2-PE stress still downregulated the transcription of these genes after Pdr1p mutation, they were significantly upregulated compared with those of the 5D strain. The expression pattern of the abovementioned genes was consistent with that of the RNA-seq data.

Fig. 6
figure 6

Transcription levels of wild strain and PDR1_862 strain at 0 g/L 2-PE and 3.5 g/L 2-PE. * Represents p < 0.05, ** represents p < 0.01, *** represents p < 0.001, ns represents no significance

Effect of Pdr1p mutation on 2-PE content of intracellular free amino acids

The above results indicated that the amino acid metabolism was subjected to disturbance. Next, intracellular free amino acids contents of the PDR1_862 and 5D strains at 0 g/L and 3.5 g/L 2-PE were analyzed. Table 2 shows the results of intracellular free amino acids. The total intracellular amino acid content of the 5D strain at 3.5 g/L 2-PE decreased by 13% than that at 0 g/L 2-PE, whereas the contents of serine, glutamate, alanine, and cysteine increased by 14%, 59%, 20%, and 37% than those in the absence of 2-PE stress, respectively. Similarly, the total intracellular amino acid content of the PDR1_862 strain at 3.5 g/L 2-PE decreased by 14% than that at 0 g/L 2-PE; whereas, the levels of aspartic acid, serine, glutamate, alanine, lysine, and arginine were 1.3, 2.03, 3.00, 2.46, 2.84, and 2.26 times higher, respectively, than those in the absence of 2-PE stress. However, under 2-PE stress, compared with that of the 5D strain, the intracellular contents of all free amino acids in the PDR1_862 strain were reduced (except for that of lysine), and the total intracellular amino acid content of the PDR1_862 strain was 31% lower than that of the 5D strain. These results demonstrated that 2-PE stress impaired intracellular amino acids homeostasis in S. cerevisiae, and Pdr1p mutation further reduced the intracellular amino acids concentration.

Table 2 Free amino acids content of PDR1_862 and 5D strains at 0 g/L and 3.5 g/L 2-PE

Subsequently, to further analyze the differences in the free amino acid contents among the four groups, principal component analysis was performed. As shown in Fig. 7A, the first two principal components (PC1, PC2) accounting for 84% of all samples revealed most information of the 17 FAAs. Among them, PC1 and PC2 accounted for 47.4% and 36.4% of the total variance, respectively. The scatter points showed that the samples were classified into four different groups, indicating that Pdr1p mutation had a significant effect on the contents of 17 free amino acids.

Fig. 7
figure 7

Comparison of free amino acids between PDR1_862 and 5D strain at 0 g/L and 3.5 g/L 2-PE stress. A and B PCA score plot and the score plot by PLS-DA analysis of free amino acids. C The VIP values of the amino acids

In order to screen potential marker amino acids among the 17 free amino acids, and partial least squares discriminant analysis (PLS-DA) was used. As shown in Fig. 7B, the scores of PLS-DA were 83%, and the samples were well separated. The marker amino acids were identified according with variable importance in the projection (VIP), and VIP > 1.0 is usually used as the selection criterion. The VIP values of cysteine, serine, glutamate, and histidine were 1.31, 1.26, 1.23, and 1.06, respectively (Fig. 7C), suggesting that cysteine was one of the most important amino acids in response to 2-PE stress, followed by serine, glutamate, and histidine. l-Cysteine metabolism plays a vital role in the metabolism of other intracellular amino acids and is closely associated with glycine, serine, and threonine metabolism; pyruvate metabolism; arginine and proline metabolism; valine, leucine, and isoleucine biosynthesis; and alanine, aspartic acid, and glutamate metabolism. It has been demonstrated that cysteine plays a significant role in the resistance against oxidative stress [21].

2-PE molecular docking

Although the mechanism by which Pdr1p mutation enhanced 2-PE stress tolerance is well understood, the molecular mechanisms of how mutant Pdr1p responds to 2-PE remain unknown. It has been demonstrated that Pdr1p directly binds to various xenobiotics and drugs [5]. We hypothesized that Pdr1p might also directly bind to 2-PE, and the Pdr1p mutation would change its molecular conformational relationship with 2-PE. In this study, molecular docking was performed to construct the complex model of Pdr1p (Fig. 8A). A certain number of hydrogen bonds were formed between the Pdr1p and the hydroxyl group of 2-PE. 2-PE bound to the internal cavity of Pdr1p with a low free energy, resulting in the formation of a relatively stable complex between 2-PE and Pdr1p. The binding energy between Pdr1p or Pdr1p C862R and 2-PE was − 5.4 kcal/mol. Compared to Pdr1p, the interaction between the mutant of C862R and 2-PE remained unchanged (Fig. 8B). The results showed that Pdr1p and 2-PE bind to form a complex in response to 2-PE stress, but Pdr1p mutation has no effect on the molecular conformational relationship of the complex. Pdr1p might have multiple interfaces to control activity of the Pdr1p mediating stress resistance [22].

Fig. 8
figure 8

Molecular docking results between 2-PE and Pdr1p (A)/Pdr1p C862R (B)

Discussion

In this sturdy, at 2-PE stress, the Pdr1p mutation increased the transcription level of the genes in the Ehrlich pathway, 2-PE content of PDR1_862 strain increased by 16% than that of 5D strain. However, the ethanol concentration of PDR1_862 strain decreased than that of the 5D strain under 2-PE stress. 2-PE was biosynthesized through the Ehrlich pathway, and NADH was needed for final reduction step in the pathway [23], while ethanol was synthesized at the expense of NADH. Therefore, under 2-PE stress, Pdr1p mutation reduced ethanol synthesis, which might increase the NADH availability to synthesize 2-PE. It has been demonstrated that the efficiency of 2-PE biosynthesis is often affected by an insufficient supply of NADH [24]. Wang et al. developed a NADH self-sufficient system for 2-PE synthesis, which resulted in 3.8-fold increase of 2-PE biocatalyst efficiency [25].

In addition, four molecules of NAD(P)H are consumed when one molecule of sulfate is reduced to produce hydrogen sulfide during sulfur metabolism [26]. 2-PE stress enhanced the transcription levels of the genes related to sulfur metabolism in 5D strain, which might increase the flux of sulfur metabolism, and in turn, increased the consumption of NAD(P)H. However, under 2-PE stress, the Pdr1p mutation attenuated their transcription and decreased the consumption of NAD(P)H. This suggested that Pdr1p mutation actually saves cellular NAD(P)H consumption through weakening sulfur metabolism, which might increase NAD(P)H availability for 2-PE biosynthesis.

The regulation of 2-PE stress tolerance by Pdr1p mutation in S. cerevisiae is global and complex. Pdr1p mutation could substantially affect multiple intracellular metabolic pathways, such as l-Phe metabolism, folate one-carbon pool metabolism, and sulfur metabolism. These metabolic pathways are closely related to intracellular amino acid metabolism [14, 27]. In this study, Pdr1p mutation altered the intracellular amino acids contents. Among the amino acids, l-cysteine was the key amino acid in response to Pdr1p mutation under 2-PE stress. Amino acid metabolism involved in alleviating stress [17, 18], such as acid stress, oxidative stress. It has been demonstrated that 2-PE stress induces intracellular ROS accumulation in S. cerevisiae, Pdr1p mutation significantly reduces intracellular ROS. l-Cysteine plays an important role in the defence against oxidative stress. Therefore, l-cysteine might be the biomarker of Pdr1p mutation in response to 2-PE stress. Exogenous amino acids might enhance the ability of S. cerevisiae to survive at 2-PE stress.

S. cerevisiae, as one of the most promising microorganisms for 2-PE biosynthesis, is a popular strain for the research on the 2-PE biosynthesis. 2-PE has a strong inhibitory effect on the yeast cells, which makes it difficult to realize scale-up production of 2-PE, while currently, the limited understanding of the 2-PE stress tolerance mechanism in yeast cells hinders the rational design and industrial application of yeast cells with high 2-PE productivity [28]. The Pdr1p mutation enhanced growth phenotype under 2-PE stress, simultaneously improved the production of 2-PE, which is important for achieving industrial production of 2-PE. This would have significant implications for the targeted design, integration of stress-resistant genetic circuits, and reprogramming of S. cerevisiae cell factories for the biosynthesis of 2-PE.

Conclusion

The Pdr1p mutation enhanced growth phenotype under 2-PE stress, simultaneously improved the production of 2-PE. Based on the results of transcriptome analysis, the Pdr1p mutation attenuated the sulfur metabolism, enhanced folate one-carbon pool and Ehrlich pathways. These pathways are associated with the amino acid metabolism. Pdr1p mutation was shown to be indeed affected the intracellular amino acids content. Pdr1p mutation reduced the intracellular free amino acid concentration under 2-PE stress, among the amino acids, cysteine was the key amino acids in response to Pdr1p mutation under 2-PE stress.

Materials and methods

Strains and medium

The S. cerevisiae PDR1_862 and S. cerevisiae CEN.PK113-5D (5D strain) were used in this study [12]. YEPD medium (20 g/L peptone, 10 g/L yeast extract, 20 g/L glucose) was used to cultivate S. cerevisiae, and 2-PE were added at concentration of 3.5 g/L or 0 g/L if needed. The 2-PE synthetic medium consisted of glucose (30 g/L), l-phenylalanine (4.5 g/L), MgSO4–7H2O (0.5 g/L), and YNB (without (NH4)2SO4) (1.7 g/L).

RNA extraction

The PDR1_862 mutant and 5D strain were cultured at 30 °C, 200 rpm for 12 h, the cultures were inoculated into YEPD containing 0 g/L or 3.5 g/L 2-PE with an initial OD600 of 0.5. Three biological replicates were set up. The cultures were collected, and total RNA was extracted. Then RNA quality was determined by 5300 Bioanalyser (Agilent) and quantified using the ND-2000 (NanoDrop Technologies). Only high-quality RNA sample (OD260/280 = 1.8–2.2, OD260/230 ≥ 2.0, RIN ≥ 6.5, 28S:18S ≥ 1.0, > 1 μg) was used to construct sequencing library. RNA purification, reverse transcription, library construction and sequencing were performed at Shanghai Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China).

Different expression analysis

To identify differential expression genes (DEGs) between two different samples, differential expression analysis was performed using the DESeq2 or DEGseq. DEGs with |log2FC|1 and FDR ≤ 0.05 (DESeq2) or FDR 0.001(DEGseq) were considered to be significantly different expressed genes. In addition, functional-enrichment analysis including GO and KEGG were performed to identify which DEGs were significantly enriched in GO terms and metabolic pathways at Bonferroni corrected p-value ≤ 0.05 compared with the whole-transcriptome background [29].

Detection of gene transcription levels

The total RNA of PDR1_862 mutant and 5D strain under 0 g/L or 3.5 g/L 2-PE were used to verify the gene transcription levels. The genomic DNA was removed from the total RNA via incubation at 42 °C for 2 min using the 4 × gDNA wiper Mix. The RNA was then reverse-transcribed into cDNA using 5 × HiScript II QRT SuperMix. This cDNA was used as a template for PCR amplification in a CFX96TM Real-Time PCR Detection System (Bio-Rad, USA) using ChamQ universal SYBR qPCR Master Mix, with ENO1 as the internal reference [30].

Intracellular amino acids determination

The yeast cells were added to 8 mL of 40% ethanol/water solution + glass beads (cells:glass beads = 1:5) and then vortexed for 15 min. The mixture was centrifuged (8000×g, 5 min). A certain amount of sample was obtained via pipetting, and 4 mL of pure hydrochloric acid was added in a 1:1 volume ratio. After nitrogen-blowing and hydrolysis, the obtained sample was passed through a 0.22-μm filter. The contents of free amino acids were then detected using an amino acid analyzer (L8900) with a column size of 200 mm × 46 mm and a cationic exchange resin column for ion-exchange chromatography. The flow rate was 0.2 mL/min, the column was 55 °C, and the reaction chamber was 138 °C.

Structural prediction and analysis of Pdr1p

Discovery Studio 2019 was used to construct the Pdr1p model. Then CB-DOCK2 will conduct cavity detection based on artificial neural networks, and utilizes Autodock Vina for docking [31]. The analysis of the interaction forces between 2-PE and Pdr1p will be completed by the PLIP online server.

Analytical methods

The cell density (OD600) was measured by spectrophotometer (752N, China). The concentration of 2-PE was measured by HPLC (Thermo Fisher Scientific, USA) according to Huili Xia et al. [12]. Ethanol content was detected by Biosensor Analyzer (Sieman Technology Co., Ltd, Shenzhen China).

Statistical and reproducibility

The data analysis and graphing were performed by GraphPad Prism 8 (GraphPad Software, USA). All experiments were conducted in three biological replicates.

Availability of data and materials

All data generated or analyzed during this study are included in this manuscript.

References

  1. Martínez-Avila O, Sánchez A, Font X, Barrena R. Bioprocesses for 2-phenylethanol and 2-phenylethyl acetate production: current state and perspectives. Appl Microbiol Biotechnol. 2018;102:9991–10004.

    Article  PubMed  Google Scholar 

  2. Qian X, Yan W, Zhang W, Dong W, Ma J, Ochsenreither K, Jiang M, Xin F. Current status and perspectives of 2-phenylethanol production through biological processes. Crit Rev Biotechnol. 2019;39:235–48.

    Article  CAS  PubMed  Google Scholar 

  3. Wang Y, Zhang H, Lu X, Zong H, Zhuge B. Advances in 2-phenylethanol production from engineered microorganisms. Biotechnol Adv. 2019;37:403–9.

    Article  CAS  PubMed  Google Scholar 

  4. Hua D, Xu P. Recent advances in biotechnological production of 2-phenylethanol. Biotechnol Adv. 2011;29:654–60.

    Article  CAS  PubMed  Google Scholar 

  5. Thakur JK, Arthanari H, Yang F, Pan SJ, Fan X, Breger J, Frueh DP, Gulshan K, Li DK, Mylonakis E, et al. A nuclear receptor-like pathway regulating multidrug resistance in fungi. Nature. 2008;452:604–9.

    Article  CAS  PubMed  Google Scholar 

  6. Simonicova L, Moye-Rowley WS. Characterizing Candida glabrata Pdr1, a hyperactive transcription factor involved in azole resistance. Methods Mol Biol. 2023;2658:169–79.

    Article  CAS  PubMed  Google Scholar 

  7. Khakhina S, Simonicova L, Moye-Rowley WS. Positive autoregulation and repression of transactivation are key regulatory features of the Candida glabrata Pdr1 transcription factor. Mol Microbiol. 2018;107:747–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Mitchell AP, Simonicova L, Moye-Rowley WS. Functional information from clinically-derived drug resistant forms of the Candida glabrata Pdr1 transcription factor. PLoS Genet. 2020;16: e1009005.

    Article  Google Scholar 

  9. Gale AN, Pavesic MW, Nickels TJ, Xu Z, Cormack BP, Cunningham KW. Redefining pleiotropic drug resistance in a pathogenic yeast: Pdr1 functions as a sensor of cellular stresses in Candida glabrata. mSphere. 2023;8: e0025423.

    Article  PubMed  Google Scholar 

  10. Matsui K, Teranishi S, Kamon S, Kuroda K, Ueda M. Discovery of a modified transcription factor endowing yeasts with organic-solvent tolerance and reconstruction of an organic-solvent-tolerant Saccharomyces cerevisiae strain. Appl Environ Microbiol. 2008;74:4222–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Carvajal E, van den Hazel HB, Cybularz-Kolaczkowska A, Balzi E, Goffeau A. Molecular and phenotypic characterization of yeast PDR1 mutants that show hyperactive transcription of various ABC multidrug transporter genes. Mol Gen Genet. 1997;256:406–15.

    Article  CAS  PubMed  Google Scholar 

  12. Xia H, Kang Y, Ma Z, Hu C, Yang Q, Zhang X, Yang S, Dai J, Chen X. Evolutionary and reverse engineering in Saccharomyces cerevisiae reveals a Pdr1p mutation-dependent mechanism for 2-phenylethanol tolerance. Microb Cell Fact. 2022;21:269.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Vu BG, Thomas GH, Moye-Rowley WS, Goldman GH. Evidence that Ergosterol biosynthesis modulates activity of the Pdr1 transcription factor in Candida glabrata. MBio. 2019;10:10–1128.

    Article  Google Scholar 

  14. Lan X, Field MS, Stover PJ. Cell cycle regulation of folate-mediated one-carbon metabolism. Wiley Interdiscip Rev Syst Biol Med. 2018;10: e1426.

    Article  PubMed  Google Scholar 

  15. Shrivastava M, Feng J, Coles M, Clark B, Islam A, Dumeaux V, Whiteway M. Modulation of the complex regulatory network for methionine biosynthesis in fungi. Genetics. 2021;217:iyaa049.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Tan Q, Zhao X, He H, Zhang J, Yi T. Carbamoyl phosphate synthetase subunit Cpa1 interacting with Dut1, controls development, arginine biosynthesis, and pathogenicity of Colletotrichum gloeosporioides. Fungal Biol. 2021;125:184–90.

    Article  CAS  PubMed  Google Scholar 

  17. Xu J, Zhao N, Meng X, Zhang T, Li J, Dong H, Wei X, Fan M. Contribution of amino acids to Alicyclobacillus acidoterrestris DSM 3922T resistance towards acid stress. Food Microbiol. 2023;113: 104273.

    Article  CAS  PubMed  Google Scholar 

  18. Suzuki T, Sugiyama M, Wakazono K, Kaneko Y, Harashima S. Lactic-acid stress causes vacuolar fragmentation and impairs intracellular amino-acid homeostasis in Saccharomyces cerevisiae. J Biosci Bioeng. 2012;113:421–30.

    Article  CAS  PubMed  Google Scholar 

  19. Wang H, Li Y, Che Y, Yang D, Wang Q, Yang H, Boutet J, Huet R, Yin S. Production of l-methionine from 3-methylthiopropionaldehyde and O-acetylhomoserine by catalysis of the yeast O-acetylhomoserine sulfhydrylase. J Agric Food Chem. 2021;69:7932–7.

    Article  CAS  PubMed  Google Scholar 

  20. Chan KX, Phua SY, Van Breusegem F. Secondary sulfur metabolism in cellular signalling and oxidative stress responses. J Exp Bot. 2019;70:4237–50.

    Article  CAS  PubMed  Google Scholar 

  21. Tikhomirova A, Rahman MM, Kidd SP, Ferrero RL, Roujeinikova A. Cysteine and resistance to oxidative stress: implications for virulence and antibiotic resistance. Trends Microbiol. 2024;32(1):93–104.

    Article  CAS  PubMed  Google Scholar 

  22. Moye-Rowley WS. Multiple interfaces control activity of the Candida glabrata Pdr1 transcription factor mediating azole drug resistance. Curr Genet. 2019;65:103–8.

    Article  CAS  PubMed  Google Scholar 

  23. Xia H, Shangguan L, Chen S, Yang Q, Zhang X, Yao L, Yang S, Dai J, Chen X. Rapamycin enhanced the production of 2-phenylethanol during whole-cell bioconversion by yeast. Appl Microbiol Biotechnol. 2022;106:6471–81.

    Article  CAS  PubMed  Google Scholar 

  24. Hwang JY, Park J, Seo JH, Cha M, Cho BK, Kim J, Kim BG. Simultaneous synthesis of 2-phenylethanol and L-homophenylalanine using aromatic transaminase with yeast Ehrlich pathway. Biotechnol Bioeng. 2009;102:1323–9.

    Article  CAS  PubMed  Google Scholar 

  25. Wang P, Yang X, Lin B, Huang J, Tao Y. Cofactor self-sufficient whole-cell biocatalysts for the production of 2-phenylethanol. Metab Eng. 2017;44:143–9.

    Article  CAS  PubMed  Google Scholar 

  26. Mendoza-Cózatl D, Loza-Tavera H, Hernández-Navarro A, Moreno-Sánchez R. Sulfur assimilation and glutathione metabolism under cadmium stress in yeast, protists and plants. FEMS Microbiol Rev. 2005;29:653–71.

    Article  PubMed  Google Scholar 

  27. Kawano Y, Suzuki K, Ohtsu I. Current understanding of sulfur assimilation metabolism to biosynthesize L-cysteine and recent progress of its fermentative overproduction in microorganisms. Appl Microbiol Biotechnol. 2018;102:8203–11.

    Article  CAS  PubMed  Google Scholar 

  28. Wang Y, Zhang Z, Lu X, Zong H, Zhuge B. Transcription factor Hap5 induces gsh2 expression to enhance 2-phenylethanol tolerance and production in an industrial yeast Candida glycerinogenes. Appl Microbiol Biotechnol. 2020;104(9):4093–107.

    Article  CAS  PubMed  Google Scholar 

  29. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Dai J, Li K, Song N, Yao W, Xia H, Yang Q, Zhang X, Li X, Wang Z, Yao L, et al. Zygosaccharomyces rouxii, an aromatic yeast isolated from chili sauce, is able to biosynthesize 2-phenylethanol via the shikimate or Ehrlich pathways. Front Microbiol. 2020;11: 597454.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Liu Y, Yang X, Gan J, Chen S, Xiao ZX, Cao Y. CB-Dock2: improved protein-ligand blind docking by integrating cavity detection, docking and homologous template fitting. Nucleic Acids Res. 2022;50:W159-w164.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundations of China (Grant Nos. 31871789 and 41876114), the key project of Hubei Provincial Department of Education (T2022011), the Science Foundation of Donghai Laboratory (No. DH-2022KF0218), the Natural Science Foundation of Zhejiang Province (No. LY18D060007), the Natural Science Foundation of Hubei Province (No. 2024AFB803), the Natural Science Foundation of Henan Province (No: 242300421581), and the Key Scientific and Technological Special of Henan Province (No: 191110110600).

Author information

Authors and Affiliations

Authors

Contributions

HL X, J D and X C designed the whole and wrote the final manuscript. DQ L and N S carried out all experiments, data collection, and manuscript editing. Ll SG and R Z participated in data analysis.

Corresponding author

Correspondence to Jun Dai.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xia, H., Song, N., Liu, D. et al. Exploring the stress response mechanisms to 2-phenylethanol conferred by Pdr1p mutation in Saccharomyces cerevisiae. Biotechnol Biofuels 17, 109 (2024). https://doi.org/10.1186/s13068-024-02559-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-024-02559-5

Keywords