High-temperature ethanol fermentation from pineapple waste hydrolysate and gene expression analysis of thermotolerant yeast Saccharomyces cerevisiae

High-temperature ethanol fermentation by thermotolerant yeast is considered a promising technology for ethanol production, especially in tropical and subtropical regions. In this study, optimization conditions for high-temperature ethanol fermentation of pineapple waste hydrolysate (PWH) using a newly isolated thermotolerant yeast, Saccharomyces cerevisiae HG1.1, and the expression of genes during ethanol fermentation at 40 °C were carried out. Three independent variables, including cell concentration, pH, and yeast extract, positively affected ethanol production from PWH at 40 °C. The optimum levels of these significant factors evaluated using response surface methodology (RSM) based on central composite design (CCD) were a cell concentration of 8.0 × 107 cells/mL, a pH of 5.5, and a yeast extract concentration of 4.95 g/L, yielding a maximum ethanol concentration of 36.85 g/L and productivity of 3.07 g/L. Gene expression analysis during high-temperature ethanol fermentation using RT–qPCR revealed that the acquisition of thermotolerance ability and ethanol fermentation efficiency of S. cerevisiae HG1.1 are associated with genes responsible for growth and ethanol stress, oxidative stress, acetic acid stress, DNA repair, the pyruvate-to-tricarboxylic acid (TCA) pathway, and the pyruvate-to-ethanol pathway.

Increasing energy demands encourage scientists to find low-cost, clean, renewable, and sustainable alternative energy sources [1][2][3] . A comparative study of literature on various alternative fuels, such as ethanol, vegetable oils, microbial oils, biomass, glycerol, biodiesel, and hydrogen, has been reported 4,5 . Commercial ethanol for biofuel is produced from feedstocks such as sugarcane, corn, and cassava. These raw materials, which are also food for human needs and animal feed, are competitively priced 6 . Agricultural wastes, particularly lignocellulosic materials, have been considered promising for second-generation bioethanol production. Pineapple peel, core, stem, and leaves are byproducts of pineapple processing (approximately 50% (w/w) of the pineapple weight) 7 . These byproducts are highly biodegradable and rich in proteins and carbohydrates, which are promising raw and abundant materials for ethanol production 8,9 . Thailand and Vietnam are the top countries in pineapple production, producing 2.21 and 0.59 million metric tons, accounting for 8.91% and 2.38% of the world's production, respectively 10 .
In summer, the temperature in Thailand and Vietnam dramatically increases, which will increase with global warming. Furthermore, the temperature inside a bioreactor may rise from 30 °C to approximately 40 °C during ethanol fermentation 11 . High temperatures inhibit cell growth and the metabolic activity of yeast cells, resulting in a reduction in ethanol yield and volumetric ethanol productivity 12,13 . Therefore, the use of thermotolerant microorganisms is a promising approach to solving the problem of ethanol production at high temperatures.  27,32 were determined for their effect on ethanol production by S. cerevisiae HG1.1. The ethanol fermentation was conducted in triplicate using a 250-mL Erlenmeyer flask containing 100 mL of PWH (pH 5.0) supplemented with various nitrogen sources at different concentrations and an initial yeast cell concentration of 5.0 × 10 6 cells/mL. All flasks were incubated at 40 °C on a rotary shaker at 100 rpm. Samples were withdrawn every 12 h and subjected to ethanol and total sugar analyses.
Optimization of ethanol production at high temperature. Based PO 4 ), and yeast extract were chosen 12,23,27,28,32,33 . The significant independent factors positively affecting ethanol production from PWH by S. cerevisiae HG1.1 were screened and selected using PBD. The codes and actual values of the independent factors are presented in Table 1. The batch ethanol fermentation experiments were performed in triplicate using a 250-mL Erlenmeyer flask containing 100 mL of PWH (pH 5.0). The ethanol concentration was set as the response variable in this study. The significant independent variables selected based on PBD were subjected to an optimization experiment using the RSM based on the CCD. The confirmatory experiment was carried out using the optimized conditions from the response surface analysis.
RT-qPCR analysis of gene expression in S. cerevisiae HG1.1 under high-temperature ethanol fermentation. The yeast inoculum was transferred into a 250-mL Erlenmeyer flask containing 100 mL PWH (pH 5.5) supplemented with 4.95 g/L yeast extract with an initial cell concentration of 8.0 × 10 7 cells/ mL. All flasks were incubated on a rotary shaker at 100 rpm under four different fermentation conditions: (1) unstressed condition (flasks were incubated at 30 °C for 9 h); (2) heat shock condition (flasks were incubated at 30 °C for 9 h, then shifted to 40 °C for 30 min); (3) short-term heat stress (flasks were incubated at 30 °C for 9 h, then shifted to 40 °C for 3 h); and (4) long-term heat stress (flasks were incubated at 40 °C for 9 h). Yeast cells were harvested at specific time points (i.e., 9 h for unstressed, 9 h and 30 min for heat shock, 12 h for short-term heat stress, and 9 h for long-term heat stress) by centrifugation at 5,000 rpm and 4 °C for 5 min, and then subjected to total RNA isolation using an RNA extraction kit (GF-1 Total RNA extraction kit, Vivantis, USA) with some modifications as described by Techaparin et al. 34 . The RNA concentration in each sample was measured and adjusted using a BioDrop μLITE (BioDrop Ltd, UK). RT-qPCR was performed in triplicate on a 7500 Fast Real-Time PCR System using the qPCRBIO SyGreen One-Step Detect Lo-ROX (PCR Biosystems, London, UK). The reactions were conducted with a total volume of 20 μL containing 1 μL RNA sample (100 ng RNA), 0.8 μL of each specific forward and reverse primer, 1 μL 20 × RTase, 10 μL 2 × qPCRBIO Sygreen One-Step mix, and 6.4 μL RNase-free water. The thermal cycling conditions were as follows: 45 °C for 30 min; 95 °C for 2 min; 40 cycles of 95 °C for 15 s, and 60 °C for 1 min. A list of primer pairs used for RT-qPCR is shown in Table 2. The RNasefree water was used instead of the RNA template for the negative control. The actin gene (ACT1) was used as an internal control. The relative gene expression was calculated using the 2 −ΔΔCT method in which the target gene amount was adjusted to the reference gene (ACT1 gene).
Analytical methods and data analysis. Viable cell concentration was determined by a haemacytometer using the methylene blue staining technique 35 . The total sugars were analyzed by the phenol sulfuric acid method 36 using a spectrophotometer (UV-1601, Shimadzu). The ethanol concentration was determined by gas chromatography (GC-14B, Shimadzu) using a packed column of polyethylene glycol (PEG-20 M) with a flame ionization detector 37 . The following equations were used to calculate the fermentation parameters: ethanol yield (Y p/s , g/g) = P E /[S 0 − S t ]; volumetric ethanol productivity (Q p , g/L.h) = P E /t; yield efficiency (E y , %) = [Y p/s /0.511] × 100; sugar consumption (S c , %) = [S t /S 0 ] × 100, where P E is ethanol concentration (g/L), S 0 is initial sugar concentration (g), S t is sugar concentration (g) at time t, and t is fermentation time (h). The data are expressed as the mean ± standard deviation (SD). Analysis of variance was used to evaluate the differences among the treatments using Duncan's multiple range tests (DMRT). The statistical analysis was carried out using Statgraphics Centurion XV (Statpoint Technologies Inc., USA).  38 reported that raw pineapple peel comprised 9.43% hemicellulose, 20.44% cellulose, and 41.21% lignin. Choonut et al. 8 showed that 51.13% hemicellulose, 37.68% cellulose, and 10.24% lignin were detected in pineapple peel after hot water pretreatment at 100 °C for 240 min. The cellulose content of pineapple waste used in this study was lower than that of other agricultural wastes, such as rice straw (32-47% cellulose) 39 and corn stover (38-40% cellulose) 40 . However, the cellulose and hemicellulose contents were greater than those of yam peel (5.7% cellulose, 5.1% hemicellulose) and cassava peel (12.7% cellulose, 5.5% hemicellulose) 41 . Glucose and fructose were the principal sugars found in PWH, accounting for 41.11 and 40.87 g/L, respectively, while sucrose and maltose were not detected. Xylose and arabinose were also detected at 5.34 and 4.42 g/L, respectively. Rattanapoltee and Kaewkannetra 31 reported that only 18.41 g/L glucose and 24.55 g/L fructose were present in pineapple peel hydrolysate. The total sugar concentration was only 55.91 g/L. Niwaswong et al. 38 reported 82.10 g/L reducing sugars by dilute acid hydrolysis of pineapple peel waste. Formic acid, acetic acid, and furfural are considered inhibitors that are derived from acid hydrolysis. The acetic acid concentration of PWH was 8.39 g/L, whereas the concentrations of formic acid and furfural were 0.96 g/L and 0.36 mg/L, respectively, which were lower than those reported by Rattanapoltee and Kaewkannetra 31 .
The total sugars of PWH were 103.03 g/L, which was as high as the total sugars found in orange peel hydrolysate (101 g/L) 33 and banana peel hydrolysate (155 g/L) 42 . Furthermore, the PWH also contains some minerals, such as nitrogen (686 mg/L), phosphorus (274 mg/L), magnesium (126 mg/L), manganese (34 mg/L), and zinc (5 mg/L), which are essential for yeast growth and metabolic activity. Based on the sugar and mineral Effect of inorganic nitrogen sources on ethanol production. The supplementation of inorganic nitrogen sources into the PWH did not significantly increase the final ethanol concentrations compared to the control treatment without inorganic nitrogen supplementation. The highest ethanol concentrations of 34.60, 34.56, and 33.55 g/L were achieved from the medium supplemented with NH 4 NO 3 , (NH 4 ) 2 SO 4 , and (NH 4 ) 2 HPO 4 , respectively, which were not significantly different from the control treatment (33.54 g/L). Furthermore, supplementation with CO(NH 2 ) 2 resulted in a lower ethanol concentration than the control (Table 3). Due to the low sugar content of PWH, S. cerevisiae HG1.1 quickly converted all sugars to reach the maximum ethanol concentration with a low consumption level of nitrogen sources. On the other hand, PWH may contain sufficient nitrogen sources for yeast growth and metabolic activity. Generally, nitrogen is essential when the fermentation process is carried out at a high initial sugar concentration (for example, greater than 200 g/L) 43 . The results in the present study coincide with those reported by Charoensopharat et al. 32 and Arora et al. 44 . The C/N ratio of the fermentation medium also played a crucial role in ethanol production. One study reported the optimum C/N ratio of 7.9 for ethanol production from sago starch using recombinant S. cerevisiae YKU 131 45 , while the other reported the value of 35.2 for ethanol production from tapioca starch using coculture of Aspergillus niger and S. cerevisiae 46 . In this study, the C/N ratio of the fermentation medium was not determined; thus, further study is needed to clarify this hypothesis. Based on the results in this study, inorganic nitrogen was not selected as the independent variable in the experiment on screening factors that affected ethanol production using a Plackett-Burman design (PBD).  (Table 4). Three factors, including initial cell concentration (A), pH (B), and yeast extract (G), were the most significant variables in ethanol production from PWH by the thermotolerant yeast S. cerevisiae HG1.1, with p values < 0.05 ( Table 5). Analysis of the adequate levels of these three crucial factors showed that the most influential factor was pH (p value was 0.0002), and yeast extract was the most negligible influential factor (p value was 0.0485). The selected model was significant (p value < 0.005), with a high confidence level based on the values of R-squared (0.9837) and adjusted R-squared (0.9551). Based on the t value limit on the Pareto chart ( Fig. 1), three variables, including initial cell concentration (A), pH (B), and yeast extract (G), were Table 3. Kinetic parameters of ethanol production from PWH supplemented with various inorganic nitrogen sources and concentrations at 40 °C for 18 h. P E : ethanol concentration (g/L); Y p/s : ethanol yield (g/g); E y : yield efficiency (%); and S c : sugar consumption (%). Mean values with different letters in the same column are significantly different (p value < 0.05) based on DMRT analysis.

Concentration (g/L)
Parameter (Mean ± SD) www.nature.com/scientificreports/ considered significant variables. All three variables positively affected ethanol production from PWH using S. cerevisiae HG1.1. In PWH with low total sugars (ca. 103 g/L), both inorganic nitrogen sources and other salts were unnecessary. However, the initial cell concentration strongly affects the ethanol production rate. Higher initial cell concentrations can promote the fermentation rate and ethanol production efficiency. Techaparin et al. 28 reported that when the initial cell concentration increased from 1.0 × 10 7 to 3.0 × 10 8 cells/mL, the ethanol concentrations from sweet sorghum juice were raised from 64.79 to 84.32 g/L using S. cerevisiae KKU-VN8. Greater than ten times higher ethanol productivity from hydrolyzed sugarcane bagasse was achieved when the inoculum size of S. cerevisiae ITV-01 was increased from 0.2 to 10 g/L 47 .

P E (g/L) Y p/s (g/g) E y (%) S c (%)
Yeast growth and fermentation activity are affected directly by the pH of the fermentation medium. The enzymes involved in the yeast growth and ethanol production pathway may be inactivated at a low pH level 48 . Although S. cerevisiae can grow well at pH values between 4.0 and 6.0, the optimum pH for ethanol production is approximately 5.0-5.5. Singh and Bishnoi 49 demonstrated that pH 5.5 was the optimum pH value from statistical optimization of ethanol production from pretreated wheat straw hydrolysate using S. cerevisiae MTCC 174. Izmirlioglu and Demirci 50 also found that pH 5.5 was the optimum value for ethanol production from potato mash waste using S. cerevisiae ATCC 24,859, which yielded a 30.99 g/L ethanol concentration. This pH value was also the optimum condition for ethanol production using K. marxianus NIRE-K3 at 45 °C, providing a 93.2% yield efficiency and 0.48 g/g ethanol yield 44 .
Yeast extract is widely used as the primary organic nitrogen source in several ethanol fermentation processes. It has been recognized as having a highly positive effect on ethanol production 44,51 . In the present study, 4.95 g/L yeast extract was the optimum concentration for ethanol production from PWH by S. cerevisiae HG1.1. Different optimum concentrations of yeast extract for ethanol production have also been reported. For instance, Schnierda et al. 51 demonstrated that 9.43 and 9.24 g/L ethanol were attained from molasses-based medium (20 g/L sugar) supplemented with 0.5 g/L total yeast assimilable nitrogen by S. cerevisiae EC1118 and I. orientalis Y1161,  www.nature.com/scientificreports/ respectively. Yeast extract at 3.0 g/L was determined to be the optimum concentration for ethanol production by S. cerevisiae NP01 using a fermentation medium containing 280 g/L sucrose 52 .
In comparison, yeast extract at 9.0 g/L promoted ethanol production from sweet sorghum juice containing 270 g/L total sugars using S. cerevisiae NP01 53 . Yeast extract is essential for efficient ethanol fermentation, especially under very high gravity fermentation conditions and high-temperature fermentation processes, but the most challenging factor is the high cost. Therefore, many scientists have tried to replace other low-cost organic nitrogen sources, such as dried spent yeast, corn-steep liquor, poultry meal, and feather meal. They have also demonstrated their potential application in the production of ethanol and other biochemicals 32,52 . Optimization conditions for ethanol production by S. cerevisiae HG1.1 using CCD. The experimental design codes and actual values of the significant independent factors, including initial cell concentration (5.0 × 10 6 to 1.0 × 10 8 cells/mL), pH (4.0 to 6.5), and yeast extract (3.0 to 12.0 g/L), are shown in Table 6. The observed ethanol concentrations from the CCD with 20 experimental runs were 19.10-33.54 g/L, and the predicted ethanol concentrations were 19.49-33.78 g/L ( Table 7). The ethanol productivities were 1.59-2.80 g/L.h. The quadratic polynomial regression model and a second-order polynomial equation to predict the final ethanol concentration (PE) as a function of the fermentation variables were established, and the prediction equation was as follows: The results revealed that the model was statistically significant (p value < 0.0001) ( Table 8). The model was reliable because the p value of lack of fit was not statistically significant (p value > 0.005), the R-squared was 0.9808, and the adjusted R-squared was 0.9635, which was close to the R-squared. The standard deviation and coefficient values were only 0.71 and 2.40%, respectively. ANOVA also demonstrated that all these factors strongly   Fig. 2. The most fixed model was achieved when the yeast extract value was fixed at 7.50 g/L, and the cell concentration and pH levels were varied. The ethanol concentration was strongly affected by both cell concentration and pH. The maximum ethanol concentration of 33.54 g/L was achieved after 15 h of fermentation at the center pH value (pH 5.5) and cell concentration of 1.0 × 10 8 cells/mL. The maximum ethanol productivity (2.80 g/L.h) was also attained. Based on the three-factor quadratic polynomial equation, the maximum predicted ethanol concentration was 33.67 g/L under the optimum conditions: cell concentration of 8.0 × 10 7 cells/mL, pH of 5.4, and yeast extract concentration of 4.9 g/L.    (Table 9). Figure 3 shows the time profile of ethanol production from PWH at 40 °C using S. cerevisiae HG1.1. The ethanol concentration quickly reached the maximal value (36.85 g/L) after 12 h of fermentation, corresponding to the dramatic decrease in total sugars (from 102.98 to 25.93 g/L). The ethanol content was slightly decreased after it reached the maximum concentration due to the oxidation of ethanol by yeast when the sugar in the fermentation medium was depleted. The remaining sugars, mostly C-5 sugars, such as xylose and arabinose, were almost unchanged since S. cerevisiae could not consume this type of sugar. The remaining total sugars in the fermented medium were 21.79 g/L. Although PWH contained some fermentation inhibitors, such as acetic acid (8.23 g/L), formic acid (0.96 g/L), and furfural (0.68 mg/L), the growth and fermentation activity of S. cerevisiae HG1.1 were not affected. The ethanol concentration, productivity, and yield efficiency achieved from PWH by S. cerevisiae HG1.1 at 40 °C were relatively high compared to several previous studies summarized in Table 10. This finding suggests that pineapple waste is a promising agricultural waste for second-generation bioethanol production.

Run Cell conc. (cells/mL) pH Yeast extract (g/L)
Ethanol (P E , g/L) www.nature.com/scientificreports/ sion studies were carried out using the synthetic medium. Only a few studies have shown the gene expression pattern using lignocellulosic materials as a feedstock. Thus, this study evaluated the expression of some groups of genes related to growth and ethanol production pathway, ethanol, oxidative and acetic acid stress, and DNA repair. The expression levels of sixteen genes responsible for growth and ethanol stress (ATP6, OLE1, ERG8), oxidative stress (GLR1, SOD1), DNA repair (RAD14, MRE11, POL4), the pyruvate-to-TCA pathway (PDA1, CIT1, LYS21), the pyruvate-to-ethanol pathway (PDC1, ADH1, ADH2), and acetic acid stress (ACS1, ALD2) in S. cerevisiae HG 1.1 were successfully evaluated using RT-qPCR. As shown in Table 11, ERG8, RAD14, ADH2, and ALD2 were up-regulated under a control growth condition (30 °C) and down-regulated under heat stress. ATP6, GLR1, SOD1, CIT1, LYS21, ADH1, and ACS1 genes were highly expressed under heat shock at 40 °C for 30 min Table 10. Comparative study on ethanol production from waste hydrolysate using S. cerevisiae and other yeast strains. T: temperature; S 0 : initial sugar concentration; P E : ethanol concentration; Q p : ethanol productivity; E y : yield efficiency; *: calculated by the authors (not given directly in the reference); -: not determined in the referenced study.

Strain
Substrate T (°C) S 0 (g/L) P E (g/L) Q p (g/L.h) E y (%) References  www.nature.com/scientificreports/ and markedly decreased under short-and long-term stress conditions. In contrast, the expression levels of OLE1, MRE11, POL4, PDA1, and PDC1 were highly increased when yeast cells were shifted from 30 °C to 40 °C for 3 h (short-term stress) and then dramatically decreased when cells were shifted from 30 °C to 40 °C for 9 h (long-term stress). Under heat stress, the expression levels of the PDC1 gene were more remarkable than that under a control condition. Relatively low expression levels of the RAD14 gene responsible for DNA repair and the ALD2 gene responsible for acetic acid stress were detected in yeast cells under long-term stress. It should be noted from the present study that genes involved in the same stress condition exhibited different expression patterns, suggesting their unique expression profile during high-temperature ethanol fermentation. ATP6, OLE1, and ERG8 are essential genes responsible for yeast growth and ethanol stress 63 . In response to heat and ethanol stresses, more ATP is needed for several biosynthesis processes that produce critical components to protect microbial cells, such as trehalose, glycogen, unsaturated fatty acids and heat shock proteins 63 . As shown in the present study, the expression of ATP6 was triggered by a heat shock condition, and its expression was slightly decreased under short-and long-term heat stresses. ATP6 mitochondrially encoded subunit a of the F0 sector of mitochondrial F1F0 ATP synthase. It integrates into the F0F1-ATPase complex and completes the process to yield a functional ATPase 64 . Thus, a high level of ATP6 expression in S. cerevisiae HG1.1 might be correlated with ATP production under heat shock conditions. In S. cerevisiae sun049T and K. marxianus DMKU 3-1042, the ATP6 gene is also highly up-regulated during high-temperature ethanol production at 38°C 65 and 45°C 60 , respectively.
OLE1, encoding a fatty acid desaturase, synthesizes monounsaturated fatty acids, such as palmitoleic acid and oleic acid, from saturated fatty acids, such as palmitic acid and stearic acid 66 . These unsaturated fatty acids and ergosterol maintain membrane fluidity as an adaptive response to the physiochemical interaction of both temperature and ethanol stresses 67 . It has been reported in S. cerevisiae that the overexpression of the OLE1 gene enhanced acetic acid tolerance and other stresses, such as ethanol, H 2 O 2 , NaCl, benzoic acid, diamide, and menadione 68 . In this study, the OLE1 gene of S. cerevisiae HG1.1 was highly expressed under heat shock and short-term heat stress, which was different from that of Qiu and Jiang 69 , who demonstrated that the OLE1 gene of S. cerevisiae M1 was approximately 2.2-3.0-fold overexpressed at 30 °C under very high gravity ethanol production. The overexpression of the OLE1 gene at 30 °C was also observed in S. cerevisiae K-9 under shaking and static sake fermentation 70 . Based on the present study, OLE1 is a heat-shock responsive gene in S. cerevisiae HG1.1. The ERG8 gene encodes phosphomevalonate kinase, which converts phosphomevalonate to diphosphomevalonate using ATP in ergosterol biosynthesis 71 . In the present study, the ERG8 gene was down-regulated under heat stress conditions, which coincided with Rossignol et al. 72 , who pointed out that most of the genes encoding proteins involved in ergosterol biosynthesis in S. cerevisiae EC1118, including ERG8, were down-regulated at high-temperature fermentation. In S. cerevisiae YZ1 and YF3, low levels of ERG8 expression were also observed under high-temperature conditions (42 °C), resulting in a reduction in ergosterol accumulation 73 .
DNA damage, including base disruption, base loss, and strand breaks, is not only induced by exposure to environmental agents, such as heat, UV rays, ROS, and oxidizing agents but also spontaneously generated during cellular metabolism 18 . RAD14, MRE11 and POL4 are common genes that encode proteins or enzymes involved in the DNA repair of yeast 74,75 . The expression of the RAD14 gene in S. cerevisiae HG1.1 was decreased under heat stress, similar to that reported by Boiteux and Jinks-Robertson 76 . The RAD14 gene is recognized as a DNA damage binding factor for nucleotide excision repair in the UV-damaged DNA of S. cerevisiae. This gene is not induced by heat stress. The other two genes, i.e., MRE11 and POL4, were up-regulated under heat shock and short-term heat stress conditions, and their expression was slightly reduced under long-term heat stress. The present results were similar to those of the MRE11 and POL4 genes in K. marxianus DMKU 3-1042, in which both genes were up-regulated under high-temperature stress 60 . It was proposed from the present study that MRE11 and POL4 are involved in the DNA repair of S. cerevisiae HG1.1 under heat stress conditions. Oxidative stress, by increasing the accumulation of reactive oxygen species (ROS), such as superoxide anions, hydrogen peroxide, and hydroxyl radicals, has been shown to cause denaturation of macromolecules, such as DNA, RNA, proteins, and lipids, in yeast cells. Several genes, including SOD1 and GLR1, are responsible for oxidative stress in yeasts. The SOD1 gene encodes superoxide dismutase, while the GLR1 gene encodes glutathione reductase. Both proteins have been shown to scavenge the superoxide anion radical and hydrogen peroxide, which can then be converted to water by the action of catalases or peroxidases 19 . The up-regulation of SOD1 and GLR1 genes under heat stress might correlate with high ROS accumulation. The overproduction of superoxide dismutase and glutathione reductase in S. cerevisiae HG1.1 might be needed to convert oxidative substrates such as superoxide anion radicals to hydrogen peroxide and finally to H 2 O. Overexpression of SOD1 and GLR1 has also been reported in K. marxianus DMKU 3-1042 and S. cerevisiae YZ1 and YF3 under heat stress 60,77 . In S. cerevisiae M1, up-regulation of the SOD1 gene and down-regulation of the GLR1 gene under high ethanol and high osmotic pressure have been reported 69 . Based on this information, the SOD1 gene, but not GLR1, can be activated by heat, ethanol, and osmotic stresses, depending on the yeast species.
PDA1 and CIT1 are involved in the pyruvate-to-TCA pathway. These genes in S. cerevisiae HG1.1 were upregulated under heat stress conditions, particularly under heat shock and short-term heat stress. Under long-term heat stress, the expression levels of both genes were slightly reduced. The expression of the PDA1 gene of S. cerevisiae HG1.1 was similar to that of S. cerevisiae Y-50316 63 . However, they differed somewhat from those reported in S. cerevisiae M1, where the PDA1 gene was down-regulated while CIT1 was up-regulated under heat stress 69 . In S. cerevisiae Y-50316, the expression of the PDA1 gene is activated not only by heat but also by ethanol stress 63 . A high expression level of the CIT1 gene has also been reported in S. cerevisiae when cells are exposed to a high temperature of 35 °C for 10 min. The increasing expression level of CIT1 increased the conversion of acetyl-CoA into the TCA pathway, leading to the accumulation of metabolic intermediates involved in the stress response 16 .
LYS21 encodes homocitrate synthase, which functions to synthesize homocitrate from acetyl-CoA and oxoglutarate. Homocitrate is a precursor for the biosynthesis of L-lysine, which plays an essential protective role www.nature.com/scientificreports/ in response to oxidative stress induced by hydrogen peroxide in S. cerevisiae 78 . Furthermore, the homocitrate synthase enzyme is also associated with the mechanism of DNA repair in the nucleus 79 . In the present study, the expression of LYS21 was enhanced under heat shock and short-term heat stress, and its expression slightly decreased after exposure to long-term heat stress. In K. marxianus DMKU 3-1042, the LYS21 gene is also upregulated under heat stress at 45°C 60 . Therefore, it was proposed from this finding that the S. cerevisiae HG1. 1 The LYS21 gene may be involved in DNA repair under heat stress. In yeast cells, PDC1, ADH1, and ADH2 are involved in a pyruvate-to-ethanol pathway. These genes are highly expressed in the stationary growth phase of K. marxianus DBKKU Y-102 under heat stress at 45°C 32 . Down-regulation of the ADH2 gene has been reported in S. cerevisiae KKU-VN8 under heat stress at 40°C 34 . In S. cerevisiae Y-50316, the expression of ADH1 and ADH2 is also induced by ethanol stress 63 . In this study, PDC1 and ADH1, but not ADH2, were up-regulated under heat stress at 40 °C, suggesting that heat stress could trigger the expression of PDC1 and ADH1 genes while suppressing the expression of ADH2 in S. cerevisiae HG1.1 during high-temperature ethanol fermentation. A high ethanol concentration produced by S. cerevisiae HG1.1 at 40 °C might also be correlated with the overexpression of PDC1 and ADH1 genes.
Several genes, including ACS1 (encoded acetate-CoA ligase) and ALD2 (encoded aldehyde dehydrogenase), are responsible for acetic acid stress. The expression of these genes in S. cerevisiae HG1.1 under heat stress was investigated in this study. The results revealed that ACS1 was up-regulated under heat shock and short-term heat stress, whereas ALD2 was down-regulated under all stress conditions. The up-regulation of ACS1 in S. cerevisiae HG1.1 under heat stress may lead to a high formation of acetic acid from acetyl-CoA but not from acetaldehyde because aldehyde dehydrogenase also utilizes NAD(P) + . The conversion of acetyl-CoA to acetic acid might generate more ATP, which can be used as an energy source for the biosynthesis of essential components or enzymes critical for yeast adaptation under heat stress. In K. marxianus DMKU 3-1042, ACS1 and ALD2 are up-regulated under heat stress at 45°C 60 , while they are highly expressed in S. cerevisiae M1 under normal growth condition (30 °C) 69 , suggesting that their expression profiles depend on the yeast species.

Conclusion
The maximum ethanol concentration of 36.85 g/L, the productivity of 3.07 g/L.h, and yield efficiency of 93.61% were achieved after fermentation of PWH using S. cerevisiae HG1.1 at 40 °C under the optimum yeast inoculum concentration of 8.0 × 10 7 cells/mL, pH of 5.5, and yeast extract concentration of 4.95 g/L. The expression of genes during high-temperature ethanol fermentation using RT-qPCR revealed that most of the genes, except ERG8, RAD14, ADH2, and ALD2, were up-regulated under heat stress conditions, particularly under heat shock and short-term heat stress. Interestingly, up-regulation of the SOD1, PDA1, CIT1, PDC1, and ADH1 genes was observed under all stresses compared to a control treatment (unstress). Although the gene expression profiles were distinctive depending on the nature and characteristics of the yeasts, the thermotolerance acquisition and fermentation efficiency of S. cerevisiae HG1.1 during high-temperature ethanol fermentation correlated with genes responsible for growth and ethanol stress, oxidative stress, acetic acid stress, DNA repair, pyruvate-to-TCA, and pyruvate-to-ethanol pathway. These results provide useful information for further advanced research to explore the regulation of the genes to benefit this potential thermotolerant yeast for producing ethanol or other valuable bio-products under high-temperature fermentation conditions. Submission declaration and verification. Submission of an article implies that the work described has not been published previously in any form.

Data availability
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.