Next Article in Journal
Residual Brewing Yeast as Substrate for Co-Production of Cell Biomass and Biofilm Using Candida maltosa SM4
Next Article in Special Issue
High Foam Phenotypic Diversity and Variability in Flocculant Gene Observed for Various Yeast Cell Surfaces Present as Industrial Contaminants
Previous Article in Journal
Techno-Economic Analysis of Integrating Soybean Biorefinery Products into Corn-Based Ethanol Fermentation Operations
Previous Article in Special Issue
Very High Gravity Bioethanol Revisited: Main Challenges and Advances
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Pichia kudriavzevii NBRC1279 and NBRC1664 to Simultaneous Saccharification and Fermentation for Bioethanol Production

1
Research Institute for Sustainable Chemistry, National Institute of Advanced Industrial Science and Technology (AIST), 3-11-32 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-0046, Japan
2
National Research Institute of Brewing, 3-7-1 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-0046, Japan
3
Department of Civil and Environmental Engineering, National Institute of Technology, Kure College, 2-2-11 Aga-minami, Kure, Hiroshima 737-8506, Japan
4
Graduate School of Integrated Sciences for Life, Hiroshima University, 1-3-1 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-8530, Japan
*
Author to whom correspondence should be addressed.
Submission received: 1 April 2021 / Revised: 18 May 2021 / Accepted: 21 May 2021 / Published: 27 May 2021
(This article belongs to the Special Issue Biofuel Fermentation from Renewable Biomass)

Abstract

:
Simultaneous saccharification and fermentation (SSF) is capable of performing enzymatic saccharification and fermentation for biofuel production in a single vessel. Thus, SSF has several advantages such as simplifying the manufacturing process, operating easily, and reducing energy input. Here, we describe the application of Pichia kudriavzevii NBRC1279 and NBRC1664 to SSF for bioethanol production. When each strain was incubated for 144 h at 35 °C with Japanese cedar particles, the highest ethanol concentrations were reached 21.9 ± 0.50 g/L and 23.8 ± 3.9 g/L, respectively. In addition, 21.6 ± 0.29 g/L and 21.3 ± 0.21 g/L of bioethanol were produced from Japanese eucalyptus particles when each strain was incubated for 144 h at 30 °C. Although previous methods require pretreatment of the source material, our method does not require pretreatment, which is an advantage for industrial use. To elucidate the different characteristics of the strains, we performed genome sequencing and genome comparison. Based on the results of the eggNOG categories and the resulting Venn diagram, the functional abilities of both strains were similar. However, strain NBRC1279 showed five retrotransposon protein genes in the draft genome sequence, which indicated that the stress tolerance of both strains is slightly different.

1. Introduction

Fossil fuels are the main energy sources in the world, and demand for these fuels is growing annually. In 2015, fossil fuel consumption was estimated at 78.4% of the global final energy consumption [1]. However, the use of fossil fuels emits greenhouse gasses, which are responsible for global warming and climate change. For these reasons, biofuels that would mitigate greenhouse gas emissions have been investigated as potential renewable energy sources to replace fossil fuels. Bioethanol is one of the most widely used biofuels globally. In 2015, bioethanol was used in 66 countries, and its production reached about 25.6 billion gallons [2].
Based on the source materials or the production method, bioethanol is generally classified as first, second, or third generation. Second-generation bioethanol is produced from lignocellulosic biomass, which is widely distributed on the earth and does not raise concerns about food sustainability [3]. The microbial process for producing second-generation bioethanol typically consists of three steps: pretreatment, enzymatic hydrolysis, and fermentation [4]. At the pretreatment step, lignocellulosic biomass is hydrolyzed into cellulose, hemicellulose and lignin by hydrothermal treatment or adding organic acids. Subsequently, at the enzymatic hydrolysis step, cellulose and hemicellulose are converted into hydrolysate containing mixed sugars. Finally, the hydrolysate is used as a carbon source in the fermentation step. However, many processes based on this method have not been put to practical use because they are not economically viable. To construct a cost-effective bioethanol production method, simultaneous saccharification and fermentation (SSF) has been studied since around the 2010s. In this method, enzymatic hydrolysis and fermentation are carried out simultaneously in the same container so that the manufacturing process can be simplified. Moreover, SSF is advantageous for its ease of operation and its reduced energy input. However, most of the enzymes used for the hydrolysis of lignocellulosic biomass show the optimum temperatures at around 50–70 °C. By contrast, the optimum temperature for fermentation (approximately 30 °C) using the budding yeast Saccharomyces cerevisiae, which is widely used for ethanol production, is lower than the optimal temperature for enzymatic hydrolysis. Therefore, a thermotolerant yeast is necessary to increase the productivity of bioethanol in SSF.
A thermotolerant yeast, Pichia kudriavzevii (previously known as Issatchenkia orientalis), has been isolated from various environments. P. kudriavzevii can catabolize a variety of different carbon sources and grows at around 45 °C as well as in high salt and sugar concentrations [5,6]. Based on those abilities, in previous studies, P. kudriavzevii species were applied to SSF for bioethanol production using soft biomass as the source material. For example, P. kudriavzevii SI is capable of producing 22.6 g/L ethanol from sugarcane bagasse [7]. Using cotton stalks, 19.8 g/L ethanol can also be produced by P. kudriavzevii HOP-1 [8]. In contrast, hard biomass has rigid cell walls and hence its enzymatic hydrolysis is difficult compared to soft biomass. Thus, there are relatively few reports of SSF using hard biomass. Moreover, no previous study has reported the application of NBRC strains of P. kudriavzevii to SSF for bioethanol production.
In the present study, we demonstrated that a combination of saccharification enzymes such as Optimash BG and Acremonium cellulase was effective for glucose production from cellulose and hemicellulose. Subsequently, we performed SSF using P. kudriavzevii NBRC1279 and NBRC1664 with particles from Japanese cedar or eucalyptus. To elucidate the characteristics of both strains, they were each subjected to genome sequencing and genome comparison.

2. Materials and Methods

2.1. Effect of Temperature on Enzymatic Hydrolysis of Cellulose and Hemicellulose

Using filter paper (3MM Chr; Whatman, Kent, UK) as a substrate, the optimal temperatures were evaluated by determining the glucose concentration at temperatures ranging from 30 °C to 65 °C. The reaction mixture (1 mL) contained 50 mM acetate buffer (pH 5.0), 10 FPU/L Acremonium cellulase (Meiji Seika Pharma, Tokyo, Japan), 20 µL/L Optimash BG (Genencor, Palo Alto, CA, USA), and 3MM Chr (1 cm × 1 cm). After the reaction mixture was incubated for 60 min at each temperature, the glucose concentration was determined by high performance liquid chromatography (HPLC). The determination of glucose is described in the following section. All data were obtained from at least triplicate experiments with the same enzymes.

2.2. Simultaneous Saccharification and Fermentation using P. kudriavzevii NBRC1279 and NBRC1664

The wood chips from Japanese cedar or eucalyptus were ground to <0.2 µm size particles by a cutter mill (MKCM-3; Masuko Sangyo, Saitama, Japan), and were dried in vacuo at 40 °C. The eucalyptus chips consisted of materials from six species (main component: Eucalyptus globulus). Japanese cedar woody particles were composed of 38.2% glucan, 6.0% xylan, 8.4% mannan, and 34.6% lignin. Japanese eucalyptus woody particles were composed of 40.0% glucan, 10.4% xylan, and 28.8% lignin. These lignocellulosic biomasses were used as substrates for SSF.
Using the dried particles of Japanese cedar or eucalyptus as source material, the SSF experiments were performed. P. kudriavzevii NBRC1279 and NBRC1664 were separately pre-cultured aerobically for 16 h at 30 °C in YPD medium, which included 10 g/L yeast extract, 20 g/L peptone, and 20 g/L glucose. Subsequently, the pre-cultures were separately collected by centrifugation, and washed twice with sterilized water. Finally, the washed cultures were diluted to OD600 = 2.0 with 30 mL of SSF medium. SSF experiments were performed in a 50 mL screw-capped vial with a magnetic stirrer under the following conditions: 3 g/L pretreated woody particles, 5 g/L yeast extract, 5 g/L peptone, 2 g/L NH4Cl, 1 g/L KH2PO4 and 0.3 g/L MgSO4·7H2O in 20 mM citrate buffer (pH 5.0) with 1000 FPU/L Acremonium cellulase and 2 mL/L Optimash BG at 120 rpm for 96 h at each temperature.
The effect of culture temperature on bioethanol production in SSF was examined at 30–45 °C. All data were obtained from at least triplicate experiments.

2.3. Quantification of Ethanol and Sugars

After each culture was clarified by centrifugation and filtration, the resultant supernatant was subjected to HPLC with an Aminex HPX-87H cationic exchange column (Bio-Rad Labs, Richmond, CA, USA). Concentrations of ethanol, glucose, and xylose were determined using a refractive index indicator. The chromatographic conditions were as follows: mobile phase, 1.5 mM H2SO4; flow rate, 0.6 mL/min; and the column oven temperature, 80 °C.

2.4. Genomic DNA Preparation

P. kudriavzevii NBRC1279 and NBRC1664 were separately cultured aerobically for 16 h at 30 °C in YPD medium. Cells were then separately harvested by centrifugation and washed twice with sterile water. Genomic DNA was extracted using the MasterPure Yeast DNA Purification Kit (Epicenter, Madison, WI, USA) and purified using AMPure XP (Beckman Coulter, Brea, CA, USA), according to the manufacturers’ instructions. The concentration and purity of genomic DNA were measured by a Quant-iT dsDNA BR assay kit (Invitrogen, Waltham, MA, USA) and a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA), respectively.

2.5. Genome Sequencing

To construct genome sequence libraries, genomic DNAs from P. kudriavzevii NBRC1279 (5 µg) and NBRC1664 (5 µg) were separately fragmented into approximately 20-kb pieces using a g-TUBE (Covaris, Brighton, UK). The resultant fragments were ligated to the sequencing adapters using a SMRTbell Template Prep Kit 1.0 (Pacific Biosciences, Menlo Park, CA, USA), yielding SMRTbell libraries. These SMRTbell libraries were selected using a BluePippin DNA Size-Selection System (15-kb size cutoff) (Sage Science, Beverly, MA, USA). The concentration and library size of SMRTbell libraries were measured using an Agilent 2200 TapeStation (Agilent Technologies, Santa Clara, CA, USA). The SMRTbell libraries were then bound to polymerases and sequencing primers using the DNA/Polymerase Binding Kit P6 v2 (Pacific Biosciences) to yield genome sequencing templates. Calculation of the concentration of the polymerase-template complex for binding and annealing reaction was performed by Binding Calculator v2.3.1.1 (Pacific Bio-sciences). Polymerase-template complexes were bound to magnetic beads using a MagBead Kit (Pacific Biosciences) and loaded on a total of 16 SMRT Cell 8 Pac v3 (Pacific Biosciences). Genome sequencing was then performed on a PacBio RS II (Pacific Biosciences).

2.6. Genome Assembly, Annotation, and Comparison

Raw data from the 16 SMRT Cells containing genomic DNAs of P. kudriavzevii NBRC1279 and NBRC1664 were assembled de novo using SMRT Analysis v2.3.0 (Pacific Biosciences; [9]) to filter subreads and CCS reads. Prediction of the open reading frame and annotation were performed with MAKER2 v2.31.3 [10], AUGUSTUS v3.0 [11], and NCBI BLAST v2.2.29+ [12]. The tRNA gene was detected using tRNAScan-SE v1.23 [13]. The Venn diagram based on their predicted coding sequences of P. kudriavzevii NBRC1279 and NBRC1664 was constructed using OrthoVenn2 [14]. The protein coding sequences of P. kudriavzevii NBRC1279 and NBRC1664 were compared against the eggNOG database.

3. Results and Discussion

3.1. Use of Optimash BG and Acremonium Cellulase for the Enzymatic Hydrolysis of Cellulose and Hemicellulose

Optimash BG, a mixture of β-glucanase and xylanase, can hydrolyze cellobiose, cellulose and hemicellulose, respectively. Acremonium cellulase has β-glucosidase activity and catalyzes the hydrolysis of cellobiose to produce bimolecular glucose. Thus, we considered that the simultaneous use of Optimash BG and Acremonium cellulase is an effective method for glucose production from lignocellulosic feedstocks.
To enhance sugar production by enzymatic hydrolysis using Optimash BG and Acremonium cellulase, we evaluated the effect of temperature at the optimal pH (5.0) for the enzymatic activity. Filter paper 3MM Chr is made of cellulose and hemicellulose, so it was used as a substrate. When the filter paper was hydrolyzed by the enzymes, glucose was obtained as the main sugar (Figure 1). When the effect of temperature on sugar production was assessed, the glucose concentration was increased with the increase in temperature, and the concentration reached the maximum (517 ± 3.00 mg/L) at 60 °C (Figure 1) with the productivity of 8.62 mg/(L·min) (Table 1).
Although simultaneous use of Optimash BG and Acremonium cellulase showed the maximum productivity at 60 °C, the rate when incubated at 65 °C was decreased to 57% (Table 1). We considered that this decrease resulted from the inactivation of Acremonium cellulase. Several cellulases from the genus Acremonium are inactivated at around 50 °C and lose their activities completely at around 70 °C [15,16].

3.2. Application of P. kudriavzevii NBRC1279 and NBRC1664 to SSF for Bioethanol Production

Japanese cedar and eucalyptus show more than 80% of holocellulose yields [17], which releases glucose and other sugars by enzymatic hydrolysis. Thus, Japanese cedar and eucalyptus were used as source materials for SSF.
When the filter paper was hydrolyzed by Optimash BG and Acremonium cellulase, glucose was obtained as the main sugar (Figure 1), which indicated that the combination of both enzymes was effective for glucose production from cellulose and hemicellulose. Thus, we next applied the combination of Optimash BG, Acremonium cellulase, and P. kudriavzevii NBRC strains for SSF to produce bioethanol.
When we started this study, the available P. kudriavzevii strains were only limited to NBRC1279 and NBRC1664; therefore, both strains were used. Moreover, both strains were capable of producing ethanol in YPD medium with glucose as the sole carbon source (Figure S1, see supplementary data); therefore, P. kudriavzevii NBRC1279 and NBRC1664 were applied to SSF. Using Japanese cedar particles as the source material, P. kudriavzevii NBRC1279 and NBRC1664 were incubated for 144 h at 35 °C, and the highest ethanol concentrations were reached 21.9 ± 0.50 g/L and 23.8 ± 3.9 g/L, respectively (Figure 2A,B). By contrast, the glucose concentrations when incubated at 35 °C were below the detection limit (Figure 2C,D). Galactose, mannose, arabinose and xylose were also generated, and their amounts were increased as time advanced (Tables S1 and S2), which was caused by carbon catabolite repression [18]. Carbon catabolite repression, a molecular mechanism, suppresses the expression of the enzymes involved in the use of alternate carbon sources in the presence of glucose. Those results indicated that the glucose produced by the enzymatic hydrolysis was consumed as the main carbon source by both strains, and then bioethanol was produced. When both strains were incubated for 144 h at 45 °C, the glucose concentrations were more than 30.0 g/L, but the ethanol concentrations were less than 2.0 g/L. These results indicated that the enzymes had sufficient activities, but that the ethanol productivities of both strains were inert.
When P. kudriavzevii NBRC1279 and NBRC1664 were incubated with Japanese eucalyptus particles, the optimum temperatures were 30 °C, and the highest ethanol concentrations were 21.6 ± 0.29 g/L and 21.3 ± 0.21 g/L, respectively (Figure 3A,B). The glucose concentrations were also below the detection limit (Figure 3C,D). As in the case of the SSF using Japanese cedar particles, galactose, mannose, arabinose and xylose were unconsumed (Tables S3 and S4). The profiles of glucose and ethanol concentrations at 45 °C were similar to those using Japanese cedar particles (Figure 2C,D).
When P. kudriavzevii NBRC1279 and NBRC1664 were applied to SSF for bioethanol production, both strains showed rapid ethanol productivities in the first 24 h, and then their rates decreased slowly (Figure 2 and Figure 3). When P. kudriavzevii NBRC1664 was incubated with Japanese cedar particles at 37–42 °C, more than 17.5 g/L ethanol was produced in 24 h (Figure 2B). However, the ethanol productivity of P. kudriavzevii NBRC 1279 was slightly lower than that of strain NBRC1664, and 11.6 g/L ethanol was produced in 24 h at 42 °C (Figure 2A). When Japanese eucalyptus particles were used as the source material, a similar trend in production profiles was observed. After 24 h of incubation at 42 °C, P. kudriavzevii NBRC1664 produced 12.7 g/L ethanol, whereas strain NBRC1279 produced 7.0 g/L ethanol (Figure 3A,B). These results demonstrated that P. kudriavzevii NBRC1279 and NBRC1664 have similar ethanol productivities at 30–37 °C, but strain NBRC1664 has superior productivity at 40–42 °C. Thus, when comparing the two strains, we considered that strain NBRC1664 has industrial potential for the application of SSF. When the hydrolysate is prepared from Japanese eucalyptus particles, acetate, furfural and 5-hydroxymethylfurfural are also generated as byproducts [19]. The organic acid and aldehydes inhibit enzymatic activity and microbial growth. In other words, as the temperature increases, the enzyme activity is increased, but the amounts of by-products are also increased. Thus, the ethanol productivity using Japanese eucalyptus particles was lower than that using Japanese cedar particles, and the error bars were larger (Figure 2 and Figure 3). In fact, a similar trend has been observed in the ethanol production using S. cerevisiae and Kluyveromyces marxianus strains [20].
Several other P. kudriavzevii strains have been previously used in SSF for bioethanol production on a laboratory scale (Table 2). At least when comparing the results shown in Table 2, the highest concentrations and productivities in this study were lower than those in previous studies, although it is difficult to accurately compare against results reported previously because the source materials and conditions in SSF are different. However, our method did not require the pretreatment of the source material, such as acid-impregnated steam explosion (AISE) [7], and special expertise for bioethanol production, which could simplify the production process.
Using a 70 L scale fermenter, we previously performed the production of bioethanol from hydrolysate prepared from Japanese eucalyptus particles [21]. In this case, the ethanol concentration reached 53.5 g/L, which was equivalent to that observed in laboratory-scale experiments. However, this method was required for the three-step processes of hydrothermal pretreatment (for 2 days), enzymatic hydrolysis (for 3 days), and fermentation (for 3 days). This processing time leads to increased cost. As in the previous method, assuming that the ethanol productivity of SSF is not decreased by the scale-up of the fermenter, our method significantly reduces the processing time. Thus, our method would be much more cost-effective and practical for bioethanol production from lignocellulosic biomass, which may be advantageous for industrial use.
After the enzymatic hydrolysis of filter paper, the maximum concentration of glucose was yielded at 60 °C (Figure 1), and the productivity was more than 80% at 50–60 °C (Table 1). However, the ethanol productivities of P. kudriavzevii NBRC1279 and NBRC1664 were significantly reduced when incubated at 45 °C. Thus, it is necessary to resolve the difference between the optimal temperature for enzymatic hydrolysis and the fermentation temperature for bioethanol production for the practical application of SSF. To reduce the difference, we are planning to overexpress CsHSP [22] or HSP20 [23] genes, which enhance heat resistance in Pichia cells, in P. kudriavzevii NBRC1279 and NBRC1664.

3.3. Difference in P. kudriavzevii NBRC1279 and NBRC1664 Based on Draft Genome Sequences

As mentioned above, strain NBRC1664 showed higher ethanol productivities at 40–42 °C than strain NBRC1279 (Figure 2 and Figure 3). To elucidate the differences in ethanol productivity and ethanol tolerance, P. kudriavzevii NBRC1279 and NBRC1664 were subjected to genome sequencing, and then genome comparison was carried out.
General genome features of P. kudriavzevii NBRC1279 and NBRC1664 are listed in Table 3. Raw data from P. kudriavzevii NBRC1279 included 192,636 reads with 201-coverage. The genome sequence was 12,851,201 bases, and the GC content was 37.6%. The assembly generated 79 contigs with a maximum length of 2,854,910 bases and an N50 contig size of 1,137,449 bases. After the prediction of the open reading frames and proteins, 4897 putative open reading frames were identified in the sequence of strain NBRC1279. We also identified 460 putative tRNA genes. Raw data from P. kudriavzevii NBRC1664 yielded 111,234 reads with 102-coverage, and the genome sequence was identified as 12,362,690 bases with a GC content of 37.8%. The genome sequence of P. kudriavzevii NBRC1664 was assembled into 76 contigs with a maximum length of 2,539,505 bases and an N50 contig size of 1,482,614 bp. Furthermore, 4815 putative open reading frames were identified.
To clarify the genetic difference between P. kudriavzevii NBRC1279 and NBRC1664, we applied the predicted coding sequences from both strains to eggNOG categories and compared them (Table 4). Among these categories, P. kudriavzevii NBRC1279 exhibited greatest proportions to function unknown (17.18%); general function prediction only (11.86%); translation, ribosomal structure, and biogenesis (8.54%); posttranslational modification, protein turnover, and chaperones (8.17%); and intracellular trafficking, secretion, and vesicular transport (7.57%). P. kudriavzevii NBRC1664 exhibited similar proportions comparable to those of P. kudriavzevii NBRC1279 (Table 4).
Subsequently, we constructed the Venn diagram based on the predicted coding sequences in the genomes (Figure 4). The Venn diagram showed 4113 orthologous genes and indicated high similarity between both P. kudriavzevii strains. In the Venn diagram, 14 and 7 paralogous genes were present in draft genome sequences of P. kudriavzevii NBRC1279 and NBRC1664, respectively (Table 5). The comparison of paralogous genes revealed that five retrotransposon protein genes are contained only within the draft genome of P. kudriavzevii NBRC1279. A retrotransposon protein gene can insert into another locus within the genome via an RNA intermediate. For example, the integration of a transposable element in Schizosaccharomyces pombe enhances the expression of stress response genes, which leads to improved growth under heat [24], oxidative [24], and heavy metal [25] stresses. Moreover, the upregulation of retrotransposon-related genes enhances ethanol tolerance in S. cerevisiae CECT10094 cells [26]. Thus, P. kudriavzevii NBRC1279 may have a different capacity of stress tolerances compared with strain NBRC1664. Further work such as transcriptome analysis is required to determine whether an association exists between the retrotransposon proteins and the stress tolerance in strain NBRC1279.
Based on the results of the eggNOG categories and the Venn diagram of predicted coding sequences, the enzymes required for ethanol production were similar (Table 4). Thus, we considered that the difference in ethanol productivities was due to the difference in the expression levels of enzymes and enzyme activities related to ethanol production. The draft genome sequences of P. kudriavzevii NBRC1279 and NBRC1664 are also useful to identify genes that may confer stress tolerance to microorganisms. In fact, we have previously succeeded in conferring tolerance to acid and salt in S. cerevisiae by overexpressing a GPI-anchored protein gene (IoGAS1) from P. kudriavzevii NBRC1279 [27,28].

3.4. Genome Properties for Ethanol Tolerance

To maintain ethanol productivity during SSF, the heat and ethanol tolerances of the host are key factors because high temperature [29] and high concentrations of ethanol [30,31] inhibit the growth of yeast. The thermotolerance mechanism of yeast has been demonstrated by gene expression profiles [29] and proteomics analysis [32,33]. Previous study has also shown that in the bioethanol production from steam-pretreated softwood based on SSF using the S. cerevisiae strain, increased ethanol concentration leads to growth inhibition [31]. Therefore, we considered that determining the ethanol tolerance-related genes of P. kudriavzevii is important because genetic engineering may be needed to avoid the toxicity of the high concentrations of ethanol produced in the cells [34]. Thus, we tried to confirm the ethanol tolerance capacities of strains NBRC1279 and NBRC1664 based on the draft genome sequences.
Transcriptome profiling of P. kudriavzevii CBS 12547 under ethanol stress has suggested that ethanol tolerance is achieved by various factors [35]. The upregulation of multiple genes involved in ergosterol biosynthesis (ERG2, ERG3, ERG27) and trehalose metabolism (TPS1) enhances the accumulation of ergosterol and trehalose, respectively, which is effective for membrane protection of P. kudriavzevii. Genes associated with heat stress response (LRE1, WSC1, SGT2) and membrane biogenesis (RRT12, GAS4, FLO1, IFF6) have also been upregulated, and these maintain cell wall integrity under ethanol stress. The upregulation of genes associated with molecular chaperones (HSP42, HSP78, HSP104) and the ubiquitin–proteasome system (UBP16, BUL2, TOM1, HUL4, BRE1, CUE2) avoids the aggregation of unfolded or misfolded proteins. In the genome sequence of P. kudriavzevii NBRC1279, upregulated genes and homologous genes were conserved. Although LRE1 was not conserved, other upregulated genes were also conserved in the genome sequence of P. kudriavzevii NBRC1664. Thus, we considered that P. kudriavzevii NBRC1279 and NBRC1664 have ethanol tolerances similar to that of P. kudriavzevii CBS 12547.

4. Conclusions

In this study, we confirmed glucose production when using both Optimash BG and Acremonium cellulase. Subsequently, we demonstrated SSF for bioethanol production using P. kudriavzevii NBRC1279 and NBRC1664 as well as Optimash BG and Acremonium cellulase. When the highest concentrations and productivities in this study were compared with those of previous studies, those of this study were lower. However, our method did not require pretreatment and special expertise, and is thus much more cost-effective and advantageous for industrial use. To elucidate the differences in the ethanol productivity and ethanol tolerance of both strains, genome sequencing and genome comparison were carried out. The draft genome sequences of P. kudriavzevii NBRC1279 and NBRC1664 revealed that both strains have similar ethanol tolerances compared with a related strain. However, further analyses, such as transcriptome analysis, are required to reveal the reasons for the different ethanol productivities between both strains at the genetic level. Based on the Venn diagram, strain NBRC1279 appeared to have a different capacity of stress tolerances from strain NBRC1664.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/fermentation7020083/s1, Figure S1: Time-dependent changes in ethanol (white circle) and glucose (black square) produced by P. kudriavzevii NBRC1279 (A) and NBRC1664 (B). Tables S1 and S2: The other sugar concentration during bioethanol production with P. kudriavzevii NBRC1279 and NBRC1664 using Japanese cedar particles. Tables S3 and S4: The other sugar concentration during bioethanol production with P. kudriavzevii NBRC1279 and NBRC1664 using Japanese eucalyptus particles.

Author Contributions

Conceptualization, H.A., A.M.; methodology, H.A., T.G., Y.I., T.S.; validation, H.A., Y.I.; formal analysis, H.A., T.G.; investigation, H.A., Y.I., Z.-i.K.; resources, A.M.; data curation, H.A., Y.I.; writing—original draft preparation, H.A.; writing—review and editing, H.A., T.G., Y.I., T.S., Z.-i.K., A.M.; visualization, H.A.; supervision, H.A., A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Draft genome sequences of P. kudriavzevii NBRC1279 and NBRC1664 were deposited in DDBJ/EMBL/GenBank databases under accession numbers BHFO01000001–BHFO01000079 and BHFP01000001–BHFP01000076, respectively.

Acknowledgments

We are grateful to all members of the Bio-conversion Research Group at our Institute (Research Institute for Sustainable Chemistry, National Institute of Advanced Industrial Science and Technology (AIST)) for their technical assistance and valuable discussion.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. REN21. Renewables 2017 Global Status Report; Renewable Energy Policy Network for the 21st Century: Paris, France, 2017; pp. 28–39. [Google Scholar]
  2. Rosales-Calderon, O.; Arantes, V. A review on commercial-scale high-value products that can be produced alongside cellulosic ethanol. Biotechnol. Biofuels 2019, 12, 240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Branco, R.H.R.; Serafim, L.S.; Xavier, A.M.R.B. Second generation bioethanol production: On the use of pulp and paper industry wastes as feedstock. Fermentation 2019, 5, 4. [Google Scholar] [CrossRef] [Green Version]
  4. Robak, K.; Balcerek, M. Review of second generation bioethanol production from residual biomass. Food Technol. Biotechnol. 2018, 56, 174–187. [Google Scholar] [CrossRef] [PubMed]
  5. Boidin, J.; Pignal, M.C.; Besson, M. Le genre Pichia sensu lato (Deuxième Contribution). Bull. Soc. Mycol. France 1965, 81, 566–606. [Google Scholar]
  6. Kurtzman, C.P.; Robnett, C.J.; Basehoar-Powers, E. Phylogenetic relationships among species of Pichia, Issatchenkia and Williopsis determined from multigene sequence analysis, and the proposal of Barnettozyma gen. nov., Lindnera gen. nov. and Wickerhamomyces gen. nov. FEMS Yeast Res. 2008, 8, 939–954. [Google Scholar] [CrossRef] [Green Version]
  7. Yuan, S.F.; Guo, G.L.; Hwang, W.S. Ethanol production from dilute-acid steam exploded lignocellulosic feedstocks using an isolated multistress-tolerant Pichia kudriavzevii strain. Microb. Biotechnol. 2017, 10, 1581–1590. [Google Scholar] [CrossRef]
  8. Oberoi, H.S.; Babbar, N.; Sandhu, S.K.; Dhaliwal, S.S.; Kaur, U.; Chadha, B.S.; Bhargav, V.K. Ethanol production from alkali-treated rice straw via simultaneous saccharification and fermentation using newly isolated thermotolerant Pichia kudriavzevii HOP-1. J. Ind. Microbiol. Biotechnol. 2012, 39, 557–566. [Google Scholar] [CrossRef]
  9. Chin, C.S.; Alexander, D.H.; Marks, P.; Klammer, A.A.; Drake, J.; Heiner, C.; Clum, A.; Copeland, A.; Huddleston, J.; Eichler, E.E.; et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat. Methods 2013, 10, 563–569. [Google Scholar] [CrossRef]
  10. Holt, C.; Yandell, M. MAKER2: An annotation pipeline and genome-database management tool for second-generation genome projects. BMC Bioinform. 2011, 12, 491. [Google Scholar] [CrossRef] [Green Version]
  11. Stanke, M.; Diekhans, M.; Baertsch, R.; Haussler, D. Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics 2008, 24, 637–644. [Google Scholar] [CrossRef] [Green Version]
  12. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  13. Lowe, T.M.; Eddy, S.R. tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25, 955–964. [Google Scholar] [CrossRef]
  14. Xu, L.; Dong, Z.; Fang, L.; Luo, Y.; Wei, Z.; Guo, H.; Zhang, G.; Gu, Y.Q.; Coleman-Derr, D.; Xia, Q.; et al. OrthoVenn2: A web server for whole-genome comparison and annotation of orthologous clusters across multiple species. Nucleic Acids Res. 2019, 47, W52–W58. [Google Scholar] [CrossRef] [Green Version]
  15. Pitson, S.M.; Seviour, R.J.; McDougall, B.M.; Woodward, J.R.; Stone, B.A. Purification and characterization of three extracellular (1→3)-β-D-glucan glucohydrolases from the filamentous fungus Acremonium persicinum. Biochem. J. 1995, 308, 733–741, Jayus. [Google Scholar] [CrossRef]
  16. McDougall, B.M.; Seviour, R.J. Purification and characterization of the (1→3)-β-glucanases from Acremonium sp. IMI 383068. FEMS Microbiol. Lett. 2004, 230, 259–264. [Google Scholar]
  17. Ono, Y.; Nakamura, Y.; Zhou, Y.; Horikawa, Y.; Isogai, A. Linear and branched structures present in high-molar-mass fractions in holocelluloses prepared from chara, haircap moss, adiantum, ginkgo, Japanese cedar, and eucalyptusLinear and branched structures present in high-molar-mass fractions in holocelluloses prepared from chara, haircap moss, adiantum, ginkgo, Japanese cedar, and eucalyptus. Cellulose 2021, in press. [Google Scholar]
  18. Gancedo, J.M. Yeast carbon catabolite repression. Microbiol. Mol. Biol. Rev. 1998, 62, 334–361. [Google Scholar] [CrossRef] [Green Version]
  19. Fujimoto, S.; Inoue, S.; Yoshida, M. High solid concentrations during the hydrothermal pretreatment of eucalyptus accelerate hemicellulose decomposition and subsequent enzymatic glucose production. Bioresour. Technol. Rep. 2018, 4, 16–20. [Google Scholar] [CrossRef]
  20. Goshima, T.; Tsuji, M.; Inoue, H.; Yano, S.; Hoshino, T.; Matsushika, A. Bioethanol production from Lignocellulosic biomass by a novel Kluyveromyces marxianus strain. Biosci. Biotechnol. Biochem. 2013, 77, 1505–1510. [Google Scholar] [CrossRef] [Green Version]
  21. Fujii, T.; Murakami, K.; Endo, T.; Fujimoto, S.; Minowa, T.; Matsushika, A.; Yano, S.; Sawayama, S. Bench-scale bioethanol production from eucalyptus by high solid saccharification and glucose/xylose fermentation method. Bioprocess Biosyst. Eng. 2014, 37, 749–754. [Google Scholar] [CrossRef] [Green Version]
  22. Wang, M.; Zou, Z.; Li, Q.; Xin, H.; Zhu, X.; Chen, X.; Li, X. Heterologous expression of three Camellia sinensis small heat shock protein genes confers temperature stress tolerance in yeast and Arabidopsis thaliana. Plant Cell Rep. 2017, 36, 1125–1135. [Google Scholar] [CrossRef]
  23. Guo, L.M.; Li, J.; He, J.; Liu, H.; Zhang, H.M. A class I cytosolic HSP20 of rice enhances heat and salt tolerance in different organisms. Sci. Rep. 2020, 10, 1383. [Google Scholar] [CrossRef]
  24. Feng, G.; Leem, Y.E.; Levin, H.L. Transposon integration enhances expression of stress response genes. Nucleic Acids Res. 2013, 41, 775–789. [Google Scholar] [CrossRef] [Green Version]
  25. Esnault, C.; Lee, M.; Ham, C.; Levin, H.L. Transposable element insertions in fission yeast drive adaptation to environmental stress. Genome Res. 2019, 29, 85–95. [Google Scholar] [CrossRef] [Green Version]
  26. Navarro-Tapia, E.; Nana, R.K.; Querol, A.; Pérez-Torrado, R. Ethanol cellular defense induce unfolded protein response in yeast. Front. Microbiol. 2016, 7, 189. [Google Scholar] [CrossRef]
  27. Matsushika, A.; Negi, K.; Suzuki, T.; Goshima, T.; Hoshino, T. Identification and characterization of a novel Issatchenkia orientalis GPI-anchored protein, IoGas1, is required for resistance to low pH and salt stress. PLoS ONE 2016, 11, e0161888. [Google Scholar] [CrossRef]
  28. Wada, K.; Fujii, T.; Akita, H.; Matsushika, A. IoGAS1, a GPI-anchored protein derived from Issatchenkia orientalis, confers tolerance of Saccharomyces cerevisiae to multiple acids. Appl. Biochem. Biotechnol. 2019, 190, 1349–1359. [Google Scholar] [CrossRef]
  29. Auesukaree, C.; Koedrith, P.; Saenpayavai, P.; Asvarak, T.; Benjaphokee, S.; Sugiyama, M.; Kaneko, Y.; Harashima, S.; Boonchird, C. Characterization and gene expression profiles of thermotolerant Saccharomyces cerevisiae isolates from Thai fruits. J. Biosci. Bioeng. 2012, 114, 144–149. [Google Scholar] [CrossRef]
  30. Zhang, Q.; Wu, D.; Lin, Y.; Wang, X.; Kong, H.; Tanaka, S. Substrate and product inhibition on yeast performance in ethanol fermentation. Energy Fuels 2015, 29, 1019–1027. [Google Scholar] [CrossRef]
  31. Alkasrawi, M.; Galbe, M.; Zacchi, G. Recirculation of process streams in fuel ethanol production from softwood based on simultaneous saccharification and fermentation. Appl. Biochem. Biotechnol. 2002, 98, 849–861. [Google Scholar] [CrossRef]
  32. Shui, W.; Xiong, Y.; Xiao, W.; Qi, X.; Zhang, Y.; Lin, Y.; Guo, Y.; Zhang, Z.; Wang, Q.; Ma, Y. Understanding the mechanism of thermotolerance distinct from heat shock response through proteomic analysis of industrial strains of Saccharomyces cerevisiae. Mol. Cell Proteom. 2015, 14, 885–897. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Xiao, W.; Duan, X.; Lin, Y.; Cao, Q.; Li, S.; Guo, Y.; Gan, Y.; Qi, X.; Zhou, Y.; Guo, L.; et al. Distinct proteome remodeling of industrial Saccharomyces cerevisiae in response to prolonged thermal stress or transient heat shock. J. Proteome Res. 2018, 17, 1812–1825. [Google Scholar] [CrossRef] [PubMed]
  34. Jhariya, U.; Dafale, N.A.; Srivastava, S.; Bhende, R.S.; Kapley, A.; Purohit, H.J. Understanding ethanol tolerance mechanism in Saccharomyces cerevisiae to enhance the bioethanol production: Current and future prospects. Bioenerg. Res. 2021, 14, 670–688. [Google Scholar] [CrossRef]
  35. Miao, Y.; Xiong, G.; Li, R.; Wu, Z.; Zhang, X.; Weng, P. Transcriptome profiling of Issatchenkia orientalis under ethanol stress. AMB Express 2018, 8, 39. [Google Scholar] [CrossRef]
Figure 1. Effect of temperature on enzymatic hydrolysis. The concentrations of glucose and xylose are indicated as white and gray bars, respectively. Error bars indicate SE (n = 3).
Figure 1. Effect of temperature on enzymatic hydrolysis. The concentrations of glucose and xylose are indicated as white and gray bars, respectively. Error bars indicate SE (n = 3).
Fermentation 07 00083 g001
Figure 2. Comparison of the time-courses for P. kudriavzevii NBRC1279 and NBRC1664 during bioethanol production using Japanese cedar particles by SSF. (A) Ethanol production with strain NBRC1279; (B) Ethanol production with strain NBRC1664; (C) Glucose production with strain NBRC1279; (D) Glucose production with strain NBRC1664. Error bars indicate SE (n = 3).
Figure 2. Comparison of the time-courses for P. kudriavzevii NBRC1279 and NBRC1664 during bioethanol production using Japanese cedar particles by SSF. (A) Ethanol production with strain NBRC1279; (B) Ethanol production with strain NBRC1664; (C) Glucose production with strain NBRC1279; (D) Glucose production with strain NBRC1664. Error bars indicate SE (n = 3).
Fermentation 07 00083 g002
Figure 3. Comparison of the time-courses for P. kudriavzevii NBRC1279 and NBRC1664 during bioethanol production using Japanese eucalyptus particles by SSF. (A) Ethanol production with strain NBRC1279; (B) Ethanol production with strain NBRC1664; (C) Glucose production with strain NBRC1279; (D) Glucose production with strain NBRC1664. Error bars indicate SE (n = 3).
Figure 3. Comparison of the time-courses for P. kudriavzevii NBRC1279 and NBRC1664 during bioethanol production using Japanese eucalyptus particles by SSF. (A) Ethanol production with strain NBRC1279; (B) Ethanol production with strain NBRC1664; (C) Glucose production with strain NBRC1279; (D) Glucose production with strain NBRC1664. Error bars indicate SE (n = 3).
Fermentation 07 00083 g003
Figure 4. The Venn diagram of predicted coding sequences of P. kudriavzevii NBRC1279 and NBRC1664.
Figure 4. The Venn diagram of predicted coding sequences of P. kudriavzevii NBRC1279 and NBRC1664.
Fermentation 07 00083 g004
Table 1. Glucose productivity.
Table 1. Glucose productivity.
Temperature (°C)Glucose Productivity [mg/(L·min)]
251.24 ± 0.08
302.01 ± 0.04
352.33 ± 0.08
403.69 ± 0.09
454.69 ± 0.12
506.86 ± 0.45
558.10 ± 0.59
608.62 ± 0.05
654.89 ± 0.31
Table 2. Comparison of bioethanol production.
Table 2. Comparison of bioethanol production.
StrainSource MaterialPretreatment MethodCulture Temperature(°C)Concentration(g/L)Productivity
[(g/(L·h)]
Reference
P. kudriavzevii HOP-1Rice strawAlkali treatment4024.31.10[8]
P. kudriavzevii SIRice strawAISE4233.41.07[7]
P. kudriavzevii NBRC1279Japanese cedar-3521.90.76 *This study
P. kudriavzevii NBRC1279Japanese eucalyptus-3021.60.72 *This study
P. kudriavzevii NBRC1664Japanese cedar-3523.80.74 *This study
P. kudriavzevii NBRC1664Japanese eucalyptus-3021.30.70 *This study
* The value was calculated from the concentration after 24 h of incubation.
Table 3. Genome features of P. kudriavzevii NBRC1279 and NBRC1664.
Table 3. Genome features of P. kudriavzevii NBRC1279 and NBRC1664.
PropertiesNBRC1279NBRC1664
Genome length (bp)12,851,20112,362,690
GC content (%)37.637.8
Contig numbers7976
Coding sequence numbers48974815
tRNA460341
Table 4. eggNOG categories of protein coding sequences in P. kudriavzevii NBRC1279 and NBRC1664.
Table 4. eggNOG categories of protein coding sequences in P. kudriavzevii NBRC1279 and NBRC1664.
NBRC1279NBRC1664
ClassDescriptionCountProportion (%)CountProportion (%)
Information storage and processing
J        Translation, ribosomal structure, and biogenesis3458.543498.85
A        RNA processing and modification541.34501.27
K        Transcription2375.872355.96
L        Replication, recombination, and repair2415.972085.27
B        Chromatin structure and dynamics390.97401.01
Cellular processes and signaling
D        Cell cycle control, cell division, chromosome partitioning360.89380.96
Y        Nuclear structure00.0000.00
V        Defense mechanisms330.82370.94
T        Signal transduction mechanisms1714.231764.46
M        Cell wall/membrane/envelope biogenesis451.11471.19
N        Cell motility70.1750.13
Z        Cytoskeleton551.36551.39
W        Extracellular structures00.0000.00
U        Intracellular trafficking, secretion, and vesicular transport3067.572977.53
O        Posttranslational modification, protein turnover, chaperones3308.173378.54
Metabolism
C        Energy production and conversion1954.831724.36
G        Carbohydrate transport and metabolism1583.911523.85
E        Amino acid transport and metabolism2055.072055.20
F        Nucleotide transport and metabolism711.76681.72
H        Coenzyme transport and metabolism932.30862.18
I        Lipid transport and metabolism1042.571002.54
P        Inorganic ion transport and metabolism992.45982.48
Q        Secondary metabolites biosynthesis, transport, and catabolism431.06471.19
Poorly characterized
R        General function prediction only47911.8646911.89
S        Function unknown69417.1867317.06
Table 5. Paralogous genes of P. kudriavzevii NBRC1279 and NBRC1664.
Table 5. Paralogous genes of P. kudriavzevii NBRC1279 and NBRC1664.
NBRC1279NBRC1664
Aldo/keto reductase family proteinCytochrome c oxidase subunit 2
Leptin receptor gene-related proteinHmra1p
Molecular chaperones HSP70/HSC70Metalloendoprotease HMP1
Monocarboxylate transporterNADH dehydrogenase subunits 2, 5, and related proteins
Peptide methionine sulfoxide reductasePhenylalanyl-tRNA synthetase, β-subunit
Predicted transporterPhospholipase
Retrotransposon protein1Uncharacterized conserved protein
Retrotransposon protein2
Retrotransposon protein3
Retrotransposon protein4
Retrotransposon protein5
Solute carrier family 25
Ubiquitin-conjugating enzyme E2
Uncharacterized protein
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Akita, H.; Goshima, T.; Suzuki, T.; Itoiri, Y.; Kimura, Z.-i.; Matsushika, A. Application of Pichia kudriavzevii NBRC1279 and NBRC1664 to Simultaneous Saccharification and Fermentation for Bioethanol Production. Fermentation 2021, 7, 83. https://0-doi-org.brum.beds.ac.uk/10.3390/fermentation7020083

AMA Style

Akita H, Goshima T, Suzuki T, Itoiri Y, Kimura Z-i, Matsushika A. Application of Pichia kudriavzevii NBRC1279 and NBRC1664 to Simultaneous Saccharification and Fermentation for Bioethanol Production. Fermentation. 2021; 7(2):83. https://0-doi-org.brum.beds.ac.uk/10.3390/fermentation7020083

Chicago/Turabian Style

Akita, Hironaga, Tetsuya Goshima, Toshihiro Suzuki, Yuya Itoiri, Zen-ichiro Kimura, and Akinori Matsushika. 2021. "Application of Pichia kudriavzevii NBRC1279 and NBRC1664 to Simultaneous Saccharification and Fermentation for Bioethanol Production" Fermentation 7, no. 2: 83. https://0-doi-org.brum.beds.ac.uk/10.3390/fermentation7020083

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop