Skip to main content

Engineering the xylose metabolism in Schizochytrium sp. to improve the utilization of lignocellulose



Schizochytrium sp. is a heterotrophic, oil-producing microorganism that can efficiently produce lipids. However, the industrial production of bulk chemicals using Schizochytrium sp. is still not economically viable due to high-cost culture medium. Replacing glucose with cheap and renewable lignocellulose is a highly promising approach to reduce production costs, but Schizochytrium sp. cannot efficiently metabolize xylose, a major pentose in lignocellulosic biomass.


In order to improve the utilization of lignocellulose by Schizochytrium sp., we cloned and functionally characterized the genes encoding enzymes involved in the xylose metabolism. The results showed that the endogenous xylose reductase and xylulose kinase genes possess corresponding functional activities. Additionally, attempts were made to construct a strain of Schizochytrium sp. that can effectively use xylose by using genetic engineering techniques to introduce exogenous xylitol dehydrogenase/xylose isomerase; however, the introduction of heterologous xylitol dehydrogenase did not produce a xylose-utilizing engineered strain, whereas the introduction of xylose isomerase did. The results showed that the engineered strain 308-XI with an exogenous xylose isomerase could consume 8.2 g/L xylose over 60 h of cultivation. Xylose consumption was further elevated to 11.1 g/L when heterologous xylose isomerase and xylulose kinase were overexpressed simultaneously. Furthermore, cultivation of 308-XI-XK(S) using lignocellulosic hydrolysates, which contained glucose and xylose, yielded a 22.4 g/L of dry cell weight and 5.3 g/L of total lipid titer, respectively, representing 42.7 and 30.4% increases compared to the wild type.


This study shows that engineering of Schizochytrium sp. to efficiently utilize xylose is conducive to improve its utilization of lignocellulose, which can reduce the costs of industrial lipid production.


Schizochytrium sp. is a heterotrophic marine fungus with an excellent lipid accumulation ability [1], reaching over 50% of the dry cell weight (DCW) under specific cultivation conditions [2]. The lipid composition of Schizochytrium sp. is simple, mainly including myristic acid (C14:0), palmitic acid (C16:0), docosapentaenoic acid (DPA) and docosahexaenoic acid (DHA) [3]. Myristic acid and palmitic acid can be used as biodiesel, and DHA is an essential fatty acid for human health. Moreover, in addition to fatty acids, Schizochytrium sp. is also considered as a good choice for squalene and carotenoids production [4, 5]. Therefore, reducing the fermentation cost of Schizochytrium sp. is great significance for their further development and application.

In recent years, industrial production based on Schizochytrium sp. has been expanding, especially for high-value products (e. g. DHA). For example, William et al. [6] screened a series of Schizochytrium sp. and Thraustochytrids strains in ocean, which were optimized for subsequent production by Omega Tech. Similarly, Schizochytrium sp. has been widely used in the production of DHA algae oil, with a production capacity of about 100,000 tonnes in 2004 [7]. However, in the industrial fermentation processes of Schizochytrium sp., the cost of raw materials still makes up a sizeable fraction of total costs, especially the organic carbon sources. According to reports, the cost of glucose makes up more than 73% of the overall cost when it is used as the primary carbon source [8]. As a result, numerous studies have tried to lower the cost of producing Schizochytrium sp. by substituting glucose with less expensive carbon sources. For example, when waste glycerol was used as the organic carbon source for the fermentation of Schizochytrium sp., 66.7 g/L of DCW and 17.3 g/L of DHA were obtained [9]. By optimizing the culture condition, Schizochytrium limacinum SR21 was able to produce 146 g/L of biomass and 82.3 g/L of lipids using acetic acid as the primary carbon source [10]. Yin et al. [11] used sugarcane molasses as a low-cost carbon source to replace glucose in fermentation medium for Schizochytrium sp., and the final titer of DHA reached 15.2 g/L after optimization of the fermentation conditions. These studies demonstrate that waste biomass resources can be used as an alternative carbon source by Schizochytrium sp. to produce lipids.

The annual global primary production of biomass is about 220 billion tonnes on dry weight basis per year [12]. It primarily comes from waste products from food industry, agriculture, municipal solid trash, and other sources. Utilizing these wastes to create some value-added products is advantageous for the environment as well as economically because dumping them in the ground or in landfills might result in major environmental issues [13]. With the synthesis of bioethanol from lignocellulose as a feedstock already becoming a well-established commercial technique, the production of high-value chemicals from lignocellulose has shown to be economically feasible [14, 15]. Lignocellulose is a promising renewable resource since it is plentiful and easily accessible [16,17,18,19]. However, there are two main limitations in the efficient utilization of lignocellulose by microorganisms. Firstly, growth inhibitors are generated during the pretreatment of lignocellulose, which can severely affect the fermentation performance. Secondly, the majority of microorganisms cannot efficiently utilize xylose for growth, despite the fact that it is the second most prevalent monosaccharide in lignocellulosic hydrolysates after glucose. To overcome the first obstacle, adaptive evolution strategies are frequently applied to increase the tolerance of microorganisms to lignocellulosic degradation products. Recently, Schizochytrium sp. mutants with great tolerance to hydrolysis by-products have been developed using adaptive laboratory evolution [20]. Moreover, detoxification techniques are often used to lower the inhibitor concentration in the hydrolysate [21]. To address the second obstacle, researchers have concentrated on constructing or enhancing xylose metabolic pathways to increase xylose utilization. To make xylose usable, it is frequently introduced into Corynebacterium glutamicum, Saccharomyces cerevisiae, Yarrowia lipolytica, and Propionibacterium freudenreichii via the heterologous xylose reductase–xylitol dehydrogenase pathway or xylose isomerase pathway [22,23,24,25]. However, there are still few reports on xylose utilization by Schizochytrium sp., necessitating further exploration.

The biomass of Thraustochytrium MAN37 was 3.1 g/L at 72 h and that of Aurantiochytrium sp. was 1.7 g/L at 120 h when cultured in a medium with xylose as the main carbon source, while thraustochytrid T18 could not grow in a medium containing xylose as the only carbon source [26,27,28]. These results indicate that different thraustochytrids species have differing capacities for utilizing xylose. Schizochytrium sp. HX-308 has a highly effective capacity for lipid synthesis. Its DCW and total lipid can reach 165.0 and 113.6 g/L, respectively, under specific circumstances [29]. Our previous research had shown that HX-308 has relatively high tolerance to furfural derivatives and organic acids found in lignocellulosic hydrolysates [30]. However, HX-308 is unable to use xylose to safeguard cell growth, which limits the utilization of lignocellulose. Here, we examined the functional activity of enzymes related to the xylose metabolism pathway of HX-308 in vitro. Furthermore, heterologous genes were introduced into HX-308 to obtain the engineered strain that can utilize xylose. Finally, the lignocellulose utilization ability of the resulting engineered strains was characterized in detail.


In vitro activity of enzymes encoded by native xylose metabolic pathway genes of HX-308

Fungi mainly utilize xylose via the redox pathway, while bacteria utilize the isomerase pathway [31, 32] (Fig. 1A). The xylulose produced by both pathways is converted into xylulose-5-phosphate by xylulokinase (XK) and channeled into the pentose phosphate pathway. In this study, multiple genes associated with the xylose metabolism pathway were identified from HX-308 by BLASTp alignment (from Hondaea fermentalgiana), including xylose isomerase (XI), xylose reductase (XR), as well as two XKs (XK3005 and XK7938) (Additional file 1). To further verify whether these four genes encode corresponding functional enzymes, the four encoded enzymes were cloned in Escherichia coli. However, only the XK3005 protein was successfully expressed as a soluble protein (Fig. 1B). Regrettably, the XI, XR and XK7938 proteins had no soluble expression even when the temperature was reduced from 28 to 16 °C [33]. It has been reported that fusing the enhanced green fluorescent protein (eGFP) to the target proteins could facilitate their expression [34,35,36]. Therefore, the eGFP was, respectively, fused to the C-terminus of the XI, XR and XK7938 proteins, and results showed that the fluorescence was successfully detected compared to the control strain (Fig. 1C). In this study, since the E. coli does not contain the XR–XDH pathway which avoids XI activity inhibition by xylitol, and does not need to add additional NADPH required in XR catalysis, then the related engineered E. coli extraction was used to test the enzyme function. As shown in Fig. 1D, when the resulting crude enzyme solution was incubated with xylose, the xylose content of the XI and XK3005 group did not decrease compared to the control, indicating that these two enzymes were inactive. However, the xylose consumption of the crude enzyme solution containing the XR and XK7938 proteins was, respectively, increased by 67.5 and 54.2%, compared to the control (Fig. 1D), confirming that these two proteins possess the predicted activity.

Fig. 1
figure 1

A Two xylose metabolism pathways in microorganisms. B Protein expression of four key genes involved in xylose metabolism of HX308. Lines 1–5: extracellular protein: Control, XK3005, XK7938, XI, XR; Lines 6–10: intracellular protein: Control, XK3005, XK7938, XI, XR. C Fluorescence diagram of strains containing three fusion proteins. D Xylose content of the crude enzyme solution system before and after incubation. XI xylose isomerase, XR xylose reductase, XDH xylitol dehydrogenase, XK xylulose kinase. Data represent the mean ± standard deviation (n = 3)

Engineering a complete xylose redox pathway in HX-308

The construction of a complete xylose metabolic pathway is essential for microbial utilization of xylose. Accordingly, we introduced the xylitol dehydrogenase (XDH) gene of Pichia pastoris (PsXDH), which is well worked in S. cerevisiae [37], to complement the xylose redox pathway in HX-308. We chose to construct the PsXDH expression cassette using the promoter of the housekeeping gene glyceraldehyde-3-phosphate dehydrogenase in order to effectively express PsXDH [38]. The expression cassette of the PsXDH gene, controlled by a strong promoter, was randomly integrated into the HX-308 genome by an Agrobacterium-mediated transformation method, resulting in the chromosomally engineered strain 308-XDH (Fig. 2A). However, there was no significant increase of biomass when 308-XDH was cultured in a medium with xylose as the sole carbon source. Moreover, when 308-XDH was cultured in a mixed sugar medium containing 20 g/L glucose and 30 g/L xylose, the xylose concentration in the medium still did not change significantly (Fig. 2B). Additionally, there was no significant difference in the amount of xylitol compared to the wild type. As a result, we speculate that 308-XDH's failure to utilize xylose may be due to the low XR activity of Schizochytrium sp. Notably, the transcript level of PsXDH was 66% lower at 60 h (after glucose depletion) than that of 24 h (before glucose depletion) (Fig. 2C). It is possible that PsXDH's ineffective expression is also to blame for 308-XDH's ineffective xylose use.

Fig. 2
figure 2

A PCR analysis of the 308-XDH. 1–5: XDH transformant, 6: wild type. B Time-course determination of residual glucose (g/L), residual xylose (g/L) and xylitol concentration (g/L) of HX-308 and 308-XDH cultivated in medium with 20 g/L glucose and 30 g/L xylose. C Quantification of the transcription levels of PsXDH in 24 h and 60 h cultivated in medium with 30 g/L glucose and 20 g/L xylose. 24 h represents the period of glucose consumption and 60 h represents the period of xylose consumption. Data represent the mean ± standard deviation (n = 3)

Overexpression of a heterologous XI gene to obtain a Schizochytrium sp. strain capable of xylose utilization

Although a heterologous XDH was successfully introduced, the engineered strain HX-308 was still unable to effectively utilize xylose. To construct a fully functional xylose isomerase pathway, we attempted to introduce a heterologous XI gene. The XI gene from the anaerobic fungus Piromyces sp. E2 was reported to have a good expression in S. cerevisiae, resulting in stronger enzyme activity than XI genes from most other sources [39, 40]. Therefore, the XI gene of Piromyces sp. E2 was introduced into HX-308 to obtain the engineered strain 308-XI (Fig. 3A). When 308-XI was cultured in a medium containing 20 g/l glucose and 30 g/l xylose, the xylose concentration decreased by 8.2 g/l after 60 h (Fig. 3B). And, the 11.3 g/L of DCW was obtained by 308-XI strain, which exceeded that of the wild type (8.81 g/L), indicating that xylose was used to convert into biomass. Furthermore, the mRNA levels of XI were 3.2-fold higher at 60 h (after glucose depletion) than at 24 h (before glucose depletion) (Fig. 3C). This indicates that overexpression of heterologous XI enabled xylose utilization by HX-308.

Fig. 3
figure 3

A PCR analysis of the 308-XI. 1–2: XI transformant, 3: wild type. B Time-course determination of residual glucose (g/L), residual xylose (g/L) and xylitol concentration (g/L) of HX-308 and 308-XI cultivated in medium with 20 g/L glucose and 30 g/L xylose. C Quantification of the transcription levels of XI in 24 h and 60 h cultivated in medium with 30 g/L glucose and 20 g/L xylose. 24 h represents the period of glucose consumption and 60 h represents the period of xylose consumption. Data represent the mean ± standard deviation (n = 3)

Overexpression of heterologous XK, a key limiting enzyme of the xylose metabolic pathway

Research has shown that XK may be a critical limiting enzyme in the xylose metabolic pathway. Overexpression of an XK gene in Y. lipolytica was a vital step in enabling the engineered strain to grow in a medium containing xylose as the sole carbon source [41]. Similarly, Thraustochytrid T18 was unable to grow in a medium with xylose as the sole carbon source until the XK gene from E. coli was introduced, combined with the overexpression of its endogenous XI [28]. In order to ensure that XK activity was strong enough to drive the ability of HX-308 to utilize xylose, in addition to native XK gene of HX-308, we selected another two kinds of XK genes from S. cerevisiae and E. coli that have been reported to possess high activity [28, 42]. Firstly, the catalytic activity of three kinds of XK proteins were tested in vitro, results showed that the activity of two heterologous was stronger than that of native XK (Additional file 1: Fig. S1). The crude enzyme solution containing XK proteins from S. cerevisiae and E. coli demonstrated an increase in xylose consumption of 24.1% and 28.8%, respectively, in comparison to the crude enzyme solution containing endogenous XK protein. Therefore, the two heterologous XK genes were selected for further genetic engineering. We then overexpressed XI from Piromyces sp. and XK from either S. cerevisiae or E. coli in HX-308, resulting in strains 308-XI-XK(S) and 308-XI-XK(E), respectively. When cultured in the medium containing 20 g/L glucose and 30 g/L xylose, the xylose consumption of 308-XI-XK(S) and 308-XI-XK(E) reached 11.1 and 9.4 g/L at 60 h, respectively (Fig. 4A), which are higher than that of the 308-XI strain, indicating that overexpression of XK further improve xylose utilization. Moreover, as shown in Fig. 4B, the final DCW of three engineered strains was all higher than that of the wild type, and the highest biomass of 308-XI-XK(S) was 15.7 g/L. It is worth mentioning that the conversion rate of xylose to biomass (g xyl/g dcw) by the 308-XI-XK(S) was 1.61, which was lower than that of the 308-XI-XK(E) (2.14 g xyl/g dcw), indicating that the 308-XI-XK(S) could produce more biomass under the condition of consuming the same amount of xylose compared another engineered strain. It has been reported that the conversion rate of xylose to biomass of the engineered S. cerevisiae and Y. lipolytica have reached a higher level [43,44,45] due to the abundance of gene editing tool, but it was better than that of Aurantiochytrium sp. [26]. Therefore, based on this result, we selected 308-XI-XK(S) for further investigation.

Fig. 4
figure 4

A Time-course determination of residual glucose (g/L), residual xylose (g/L) and xylitol concentration (g/L) of 308-XI-XK(S) and 308-XI-XK(E) cultivated in medium with 20 g/L glucose and 30 g/L xylose. B The DCW and the conversion rate of xylose to biomass (g xyl /g dcw) for HX-308, 308-XI, 308-XI-XK(S) and 308-XI-XK(E) at 60 h. Data represent the mean ± standard deviation (n = 3)

Carbon source utilization and lipid synthesis in 308-XI-XK(S) grown on corn straw hydrolysate

In order to further investigate the utilization of lignocellulose by 308-XI-XK(S), we conducted fermentation experiments in a medium with corn stover hydrolysate as the main carbon source. In the hydrolysate medium, we first set the initial sugar concentration to 30 g/L glucose and 10 g/L xylose. According to the results, xylose and glucose in the 308-XI-XK(S) medium were totally consumed at 60 h, while there was no significant xylose consumption in the wild type (Fig. 5A). At 60 h, the DCW and total lipids of 308-XI-XK(S) reached 22.4 and 5.3 g/L, representing 42.7 and 30.4% increases over the control strain, respectively (Fig. 5B). This demonstrates that overexpression of XI and XK can enhance the ability of Schizochytrium sp. to produce lipids from lignocellulose. A minimal amount of xylitol is still produced in 308-XI-XK(S) medium (Fig. 5A). Therefore, xylitol synthesis gene knockdown can be tried in an attempt to reduce xylose waste. Additionally, the fatty acid profile of 308-XI-XK(S) mainly included 10.4% C14:0, 16.2% C16:0, 17.7% DPA, and 50% DHA, similar to the wild type. Thus, the genetic modification technique used in this study increased lignocellulose utilization without changing the strains' fatty acid profiles.

Fig. 5
figure 5

Sugar consumption (A), DCW, TL and fatty acid composition (B) of HX-308 and 308-XI-XK(S) cultured in corn stover hydrolysate medium for 72 h. DCW dry cell weight, TL total lipid. Data represent the mean ± standard deviation (n = 3)


The cost of the culture medium has a significant impact on the economics of lipid production [46]. In order to lower the cost of fermentation substrates for Schizochytrium sp., it is crucial to use cheap carbon sources in place of glucose. To this end, sweet sorghum juice and coconut water have been investigated as fermentation substrates to lower the cost of lipid production [47, 48], but these strategies run the risk of competing with food. By contrast, non-food lignocellulose has been extensively exploited as a low-cost, easily accessible renewable material to lower the cost of fermentation for oleaginous microorganisms. Xylose is one of the main monosaccharides produced after pretreatment of lignocellulose, but most microorganisms cannot utilize xylose, limiting their effective use of lignocellulosic biomass.

The introduction of heterologous genes encoding efficient enzymes of the xylose utilization pathway is a common strategy to engineer microorganisms that cannot utilize xylose. S. cerevisiae lacks the enzyme necessary to produce xylulose, a crucial step in the xylose metabolic pathway, and is consequently unable to utilize xylose for growth. However, S. cerevisiae was successfully engineered to produce ethanol from xylose by introducing heterologous XR and XDH genes [32]. In fact, whole-genome sequencing and gene annotation of Schizochytrium sp. did not find a natural XDH gene. Although the heterologous XDH (PsXDH) was introduced, there was no significant change in xylitol. This phenomenon may be a result of the protein expression being restricted by the transcriptional level [49], that is, PsXDH cannot be efficiently expressed in HX-308. Another possibility is that the PsXDH protein successfully expressed in Schizochytrium sp. does not possess the corresponding functional activity. In contrast, the XI pathway makes more efficient use of xylose and reduces the production of xylitol, which is essential. It has been shown that the accumulation of large amounts of xylitol can inhibit normal host growth. When present at higher levels, it can lead to further cell lysis. The XI pathway avoids this phenomenon and reduces the accumulation of xylitol in HX-308 [28]. Additionally, the selection of an XDH with efficient activity in HX-308 is one of the optimal strategies to reduce xylitol and utilize xylose efficiently. Similar strategies have been successfully implemented in the introduction of XI. S. cerevisiae was engineered by introducing the XI gene from Clostridium thermosulfurogenes, and northern blot analysis revealed that the heterologous XI was efficiently transcribed. However, the recombinant strain was still unable to grow on xylose as the only carbon source [50]. In addition, Parachin et al. [51] enabled S. cerevisiae to grow on xylose by introducing a novel XI screened from a soil metagenome library. Therefore, in order to obtain Schizochytrium sp. that can utilize xylose, we can attempt to express more XDHs from various sources in Schizochytrium sp.

Schizochytrium sp. was able to use xylose when heterologous XI was overexpressed, and additional XK overexpression resulted in a further increase of the xylose utilization rate (Fig. 4). It was shown that higher XK activity worked in concert with XI to improve xylose utilization in Schizochytrium sp. The isomerase pathway has the advantage of having fewer metabolic steps than the redox pathway, and it is not necessary to consider cofactor balance when engineering strains based on the isomerase pathway. However, none of the strains were able to effectively utilize xylose as the sole carbon source in this study. It has been shown that the gene numbers of XI and XK is critical for the efficient use of xylose for growth. Several studies have shown that the gene numbers of XI and XK is critical for the efficient use of xylose for growth. For example, with increasing amounts of XI and XK, the ability of Thraustochytrid T18 to grow in a medium with xylose as the sole carbon source also increased [28]. Similarly, the same phenomenon exists in S. cerevisiae [52]. Therefore, this is likely to be the key to limiting the efficient use of xylose by 308-XI-XK(S) for normal cellular metabolism. However, 308-XI-XK(S) still utilizes xylose much less efficiently than glucose, and novel strategies may be needed to obtain Schizochytrium sp. strains with improved xylose utilization. Zhou et al. [52] combined genetic engineering and adaptive evolutionary strategies to modify S. cerevisiae, yielding a mutant strain with a high xylose utilization rate of 1.866 g/g/h.

Enabling microorganisms to use xylose is essential for the efficient production of bulk chemicals and biofuels from cheap and abundant lignocellulosic biomass. In this study, we improved the lignocellulose utilization of Schizochytrium sp. by enhancing its xylose utilization capacity. In addition, Schizochytrium sp. strains have recently been engineered for growth on other components present in lignocellulosic materials, such as cellobiose [53]. Future combinations of these approaches might produce more economically viable Schizochytrium sp., which would be beneficial not only for large-scale production of microbial lipids, but also for other fermentation products such as squalene, carotenoids, and astaxanthin.


In this study, we demonstrated that the endogenous XR and XK7938 genes of Schizochytrium sp. HX-308 have corresponding functional activities. In addition, after testing xylose metabolism-related genes from different sources, we engineered the strain 308-XI-XK(S), which is capable of efficiently producing lipids from xylose. When 308-XI-XK(S) was cultured in a medium with lignocellulose as the main carbon source, the strain’s xylose utilization rate reached 0.16 g/L/h, while its biomass and lipid production were significantly higher than those of the control. In conclusion, this study provides a reference for further improving the production of lipids by Schizochytrium sp. using lignocellulose as a cheap carbon source.

Materials and methods

Strains and plasmids

Schizochytrium sp. HX-308 (CCTCC M 209059) used in this study was stored in the China Center for Type Culture Collection [54]. Agrobacterium tumefaciens AGL-1 and E. coli BL21 (DE3) were kind gifts from Dr. Sheng Yang (CAS Center for Excellence in Molecular Plant Science). All plasmids used in this study are listed in Additional file 1: Table S1, and the DNA primers are listed in Additional file 1: Table S2. For efficient expression and execution of enzyme assays, XI, XR, XK3005 and XK7938 proteins were codon optimized for E. coli (Additional file 1: Table S4). All plasmids used to express recombinant proteins were constructed based on the pET-24a backbone using Gibson assembly. During plasmid construction, eGFP was fused to the C-terminus of the recombinant protein. The constructed plasmid was introduced into E. coli BL21 (DE3) by chemical transformation. All plasmids used to carry DNA elements were constructed based on pZPK-NeoR.

Culture conditions

Protein expression was conducted as described before [55]. The indicated E. coli BL21 (DE3) with recombinant plasmids were first activated in a 3 mL LB liquid medium at 37 °C and 220 rpm for 12 h. Then, 300 μL of the resulting seed culture was used to inoculate 250 mL shake flasks containing 30 mL of TB medium and incubated at 37 °C and 220 rpm. When the cell OD600 reached 2–4, IPTG was added to a final concentration of 0.3 mM to induce recombinant protein expression, after which the fermentation was continued at 28 °C. After 24 h, the cells were harvested by centrifugation at 5000 g for 10 min and 4 °C, resuspended in 100 mM Tris–HCl (pH 7.5), and lysed by ultrasonication. The crude cell lysate was centrifuged at 5000 g and 4 °C for 10 min to remove the cell debris and used for enzyme activity assays.

The culture medium and culture conditions for Schizochytrium sp. were the same as in our previous study [56]. The main components of Schizochytrium sp. medium consisted of glucose 50 g/L, yeast extract 10 g/L, sodium sulfate 12 g/L, magnesium sulfate 2 g/L, ammonium sulfate 4 g/L, potassium chloride 0.7 g/L, calcium chloride 0.2 g/L, sodium glutamate 10 g/L, and potassium dihydrogen phosphate 1 g/L. 1.0 mL of Schizochytrium sp. stored in glycerol was inoculated in a shaker containing 50 mL of medium and 250 mL of shaker with a baffle for activation. After 24 h of incubation, the activated cultures were inoculated into 500-mL shake flasks with baffles containing 100 mL of mixed sugar medium for incubation.

Identification of the xylose metabolism genes

For Schizochytrium sp. HX-308, the xylose metabolism genes were obtained by whole-genome sequencing and functional annotation. Further, through BLASTp in NCBI, the translated protein sequences show high identity with xylose metabolic pathway-related enzymes, including XI, XR and XK (NCBI accession number: txid2315210, from Hondaea fermentalgiana), whose functions have been identified. The accession numbers of the sequences reported in this article deposited in GenBank: OP428812 (XI), OP428813 (XR), OP428815 (XK3005) OP428814 (XK7938).

Transformation of HX-308

The transformation method of HX-308 was described previously [57]. Briefly, HX-308 was grown overnight, and then transferred to a fresh medium and incubated until the logarithmic phase. The cells were then collected by centrifuging at 1000 g for 5 min, and washed twice with sterile water. A. tumefaciens carrying the target plasmid was grown in YEB medium at 28 °C and 220 rpm for 12 h. The bacterial donor cells were collected by centrifugation at 4000 g for 5 min, resuspended in IM induction medium containing 50 μg/mL kanamycin, and pre-induced at 28 °C and 220 rpm for 8 h. The pre-induced cells were then collected by centrifugation at 4000 g for 5 min and washed twice with sterile water. HX-308 at a concentration of 1 × 107 and A. tumefaciens at a concentration of 1 × 108 were mixed and then spread on IM solid medium containing 50 μg/mL kanamycin. After incubation at 28 °C for 48 h, the co-cultures were washed with sterile water, spread on GPYS solid medium containing 300 g/mL cefotaxime sodium and 500 g/mL G418, and incubated for 3 to 5 days to screen for positive transformants. Finally, the presence of the transgene was confirmed by colony PCR.

Preparation of corn straw hydrolysate

Corn straw hydrolysate was prepared according to Yuan et al. [58]. In a 2 L flask, 500 g of a digestion mixture containing 10% w/w biomass and 2% w/w NaOH was heated at 121 °C for 30 min. After cooling, hydrochloric acid was used to bring the slurry's pH value up to 7.0, and the neutralized biomass was dried in a 60 °C oven until its moisture content was between 10 and 20%. Following pretreatment, cellulase and hemicellulase were used to hydrolyze the biomass at pH 4.8–5.0, 50 °C, and 220 rpm. The hydrolyzed slurry was centrifuged at 6000 g for 10 min after 72 h to obtain a liquid hydrolysate without solids.

Measurement of glucose, xylose, DCW, total lipid, fatty acids and inhibitor

The supernatant from the centrifuged fermentation broth was used to measure glucose concentrations as described previously [59]. The DCW was determined by gravimetric analysis after centrifuging a volume of fermentation broth at 6000 g for 5 min, removing the supernatant, washing the precipitate twice with water, and then drying in an oven at 60 °C to a constant weight. Lipid analysis was performed as previously reported [60]. Formic acid, acetic acid, and xylose concentrations were determined using a refractory index detector (RID-6A) and an Aminex HPX-87H column maintained at 65 °C. The mobile phase contained 5.0 mM H2SO4 and flowed at a rate of 0.6 mL/min. Furan and phenolic compounds were analyzed by HPLC using a UV/V is detector (SPD-20A) and a YMC-Pack ODS-A column kept at 35 °C. The mobile phase contained 50% v/v aqueous acetonitrile and had a flow rate of 1.0 mL/min at a detection wavelength of 220 nm. The primary components in the corn straw hydrolysate included the 90.1 g/L glucose, 29.3 g/L xylose, 1.3 g/L vanillin, 1.6 g/L acetic acid, 0.2 g/L furfural, and 0.4 g/L phenol. The conversion rate of xylose to biomass (g xyl /g dcw) for HX-308, 308-XI, 308-XI-XK(S) and 308-XI-XK(E) was calculated according to equation: Xylose consumption/(DCW20 g/L Glu+30 g/L Xyl- DCW20 g/L Glu).

Enzyme activity assays with crude cell lysates

The reaction system contained 1 g/L xylose and an appropriate amount of crude enzyme solution. After 12 h at 30 °C, the enzymatic reaction was stopped by boiling for 5 min, and the xylose consumption after was measured by HPLC.

Quantitative reverse-transcription PCR (qRT-PCR)

To observe changes in transcript levels of heterologous XDH and XI genes, total RNA samples were extracted from cells grown in xylose medium with and without glucose using a fungal RNA mini kit (Omega BioTek, USA). The cDNA synthesis and quantitative real-time RT-PCR (qRT-PCR) analysis were performed by Zoonbio Biotechnology Inc. (Nanjing, China). Gene expression levels were calculated using the 2−ΔΔCT method and normalized to the expression level of glucose-6-phosphate dehydrogenase. The primer sequences used in this study are shown in Additional file 1: Table S2.

One-sentence summary

The engineering of Schizochytrium sp. to efficiently utilize xylose is conducive to improving its utilization of lignocellulose, which can greatly reduce the costs of industrial lipid production.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files.



Dry cell weight


Docosapentaenoic acid


Docosahexaenoic acid


Xylose isomerase


Xylulose kinase


Xylitol dehydrogenase


Xylose reductase


Enhanced green fluorescent protein


  1. Jakobsen AN, Aasen IM, Josefsen KD, Strm AR. Accumulation of docosahexaenoic acid-rich lipid in thraustochytrid Aurantiochytrium sp strain T66: effects of N and P starvation and O2 limitation. Appl Microbiol Biotechnol. 2008;80:297–306.

    Article  CAS  PubMed  Google Scholar 

  2. Li J, Liu R, Chang G, Li X, Chang M, Liu Y, Jin Q, Wang X. A strategy for the highly efficient production of docosahexaenoic acid by Aurantiochytrium limacinum SR21 using glucose and glycerol as the mixed carbon sources. Biores Technol. 2015;177:51–7.

    Article  CAS  Google Scholar 

  3. Jiang Y, Fan K, Wong R, Chen F. Fatty acid composition and squalene content of the marine microalga Schizochytrium mangrovei. J Agric Food Chem. 2004;52:1196–200.

    Article  CAS  PubMed  Google Scholar 

  4. Zhang K, Chen L, Liu J, Gao F, He R, Chen W, Guo W, Chen S, Li D. Effects of butanol on high value product production in Schizochytrium limacinum B4D1. Enzyme Microb Technol. 2017;102:9–15.

    Article  CAS  PubMed  Google Scholar 

  5. Yang Q, Xie Z, Zheng X, Li K, Lu T, Lu Y, Chen C, Ling X. Genetic regulation and fermentation strategy for squalene production in Schizochytrium sp. Appl Microbiol Biotechnol. 2022;106:2415–31.

    Article  CAS  PubMed  Google Scholar 

  6. Barclay WR. Process for the heterotrophic production of microbial products with high concentrations of omega-3 highly unsaturated fatty acids. EP. 2009;EP1707061 A3.

  7. Brennan L, Owende P. Biofuels from microalgae-A review of technologies for production, processing, and extractions of biofuels and co-products. Renew Sustain Energy Rev. 2010;14:557–77.

    Article  CAS  Google Scholar 

  8. Du F, Wang YZ, Xu YS, Shi TQ, Huang H. Biotechnological production of lipid and terpenoid from thraustochytrids. Biotechnol Adv. 2021;48: 107725.

    Article  CAS  PubMed  Google Scholar 

  9. Kujawska N, Talbierz S, Dbowski M, Kazimierowicz J, Zieliński M. Optimizing docosahexaenoic acid (DHA) production by Schizochytrium sp. grown on waste glycerol. Energies. 2021;14:1685.

    Article  CAS  Google Scholar 

  10. Shafiq M, Zeb L, Cui G, Jawad M, Chi Z. High-density pH-auxostat fed-batch culture of Schizochytrium limacinum SR21 with acetic acid as a carbon source. Appl Biochem Biotechnol. 2020;192:1163–75.

    Article  CAS  PubMed  Google Scholar 

  11. Yin FW, Zhu SY, Guo DS, Ren LJ, Ji XJ, Huang H, Gao Z. Development of a strategy for the production of docosahexaenoic acid by Schizochytrium sp. from cane molasses and algae-residue. Bioresour Technol. 2019;271:118–24.

    Article  CAS  PubMed  Google Scholar 

  12. Kumari D, Singh R. Pretreatment of lignocellulosic wastes for biofuel production: a critical review. Renew Sustain Energy Rev. 2018;90:877–91.

    Article  CAS  Google Scholar 

  13. Taherzadeh MJ, Karimi K. Pretreatment of lignocellulosic wastes to improve ethanol and biogas production: a review. Int J Mol Sci. 2008;9:1621–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Isikgor FH, Becer CR. Lignocellulosic biomass: a sustainable platform for the production of bio-based chemicals and polymers. Polym Chem. 2015;6:4497–559.

    Article  CAS  Google Scholar 

  15. Cherubini F, Stromman AH. Chemicals from lignocellulosic biomass: opportunities, perspectives, and potential of biorefinery systems. Biofuels Bioproducts Biorefining-Biofpr. 2011;5:548–61.

    Article  CAS  Google Scholar 

  16. Chandra R, Takeuchi H, Hasegawa T. Methane production from lignocellulosic agricultural crop wastes: a review in context to second generation of biofuel production. Renew Sustain Energy Rev. 2012;16:1462–76.

    Article  CAS  Google Scholar 

  17. Chen X, Li Z, Zhang X, Hu F, Ryu DDY, Bao J. Screening of oleaginous yeast strains tolerant to lignocellulose degradation compounds. Appl Biochem Biotechnol. 2009;159:591.

    Article  CAS  PubMed  Google Scholar 

  18. Tsigie YA, Wang CY, Truong CT, Ju YH. Lipid production from Yarrowia lipolytica Po1g grown in sugarcane bagasse hydrolysate. Bioresour Technol. 2011;102:9216–22.

    Article  CAS  PubMed  Google Scholar 

  19. Dipesh K, Bhaskar S, John K. Utilization of lignocellulosic biomass by oleaginous yeast and bacteria for production of biodiesel and renewable diesel. Renew Sustain Energy Rev. 2017;73:654–71.

    Article  Google Scholar 

  20. Qi F, Zhang M, Chen Y, Jiang X, Lin J, Cao X, Huang J. A lignocellulosic hydrolysate-tolerant Aurantiochytrium sp. mutant strain for docosahexaenoic acid production. Bioresour Technol. 2017;227:221–6.

    Article  CAS  PubMed  Google Scholar 

  21. Xiong L, Huang C, Yang XY, Lin XQ, Chen XF, Wang C, Wang B, Zeng XA, Chen XD. Beneficial effect of corncob acid hydrolysate on the lipid production by oleaginous yeast Trichosporon dermatis. Prep Biochem. 2015;45:421–9.

    Article  CAS  Google Scholar 

  22. Prabhu AA, Ledesma-Amaro R, Lin CSK, Coulon F, Thakur VK, Kumar V. Bioproduction of succinic acid from xylose by engineered Yarrowia lipolytica without pH control. Biotechnol Biofuels. 2020;13:113.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Wei PL, Lin M, Wang ZQ, Fu HX, Yang HP, Jiang WY, Yang ST. Metabolic engineering of Propionibacterium freudenreichii subsp shermanii for xylose fermentation. Biores Technol. 2016;219:91–7.

    Article  CAS  Google Scholar 

  24. Kaneda J, Sasaki K, Gomi K, Shintani T. Heterologous expression of Aspergillus oryzae xylose reductase and xylitol dehydrogenase genes facilitated xylose utilization in the yeast Saccharomyces cerevisiae. Biosci Biotechnol Biochem. 2011;75:168–70.

    Article  CAS  PubMed  Google Scholar 

  25. Kang MK, Lee J, Um Y, Lee TS, Bott M, Park SJ, Woo HM. Synthetic biology platform of CoryneBrick vectors for gene expression in Corynebacterium glutamicum and its application to xylose utilization. Appl Microbiol Biotechnol. 2014;98:5991–6002.

    Article  CAS  PubMed  Google Scholar 

  26. Mariam I, Kareya MS, Nesamma AA, Jutur PP. Delineating metabolomic changes in native isolate Aurantiochytrium for production of docosahexaenoic acid in presence of varying carbon substrates. Algal Res Biomass Biofuels Bioprod 2021;55:102285.

  27. Miranda AF, Tran TLN, Abramov T, Jehalee F, Miglani M, Liu ZQ, Rochfort S, Gupta A, Cheirsilp B, Adhikari B, Purid M, Mouradova A. Marine protists and Rhodotorula yeast as bio-convertors of marine waste into nutrient-rich deposits for mangrove ecosystems. Protist. 2020;171: 125738.

    Article  CAS  PubMed  Google Scholar 

  28. Merkx-Jacques A, Rasmussen H, Muise DM, Benjamin JJR, Kottwitz H, Tanner K, Milway MT, Purdue LM, Scaife MA, Armenta RE, Woodhall DL. Engineering xylose metabolism in thraustochytrid T18. Biotechnol Biofuels. 2018;11:248.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Guo DH, Ji XJ, Ren LJ, Li GL, He H. Improving docosahexaenoic acid production by Schizochytrium sp. using a newly designed high-oxygen-supply bioreactor. AICHE J. 2017;63:4278–86.

    Article  CAS  Google Scholar 

  30. Wang LR, Zhang ZX, Wang YZ, Li ZJ, Huang PW, Sun XM. Assessing the potential of Schizochytrium sp. HX-308 for microbial lipids production from corn stover hydrolysate. Biotechnol J. 2022;17:e2100470.

    Article  PubMed  Google Scholar 

  31. Jeffries TW, Jin YS. Metabolic engineering for improved fermentation of xylose by yeasts. Appl Microbiol Biotechnol. 2004;63:495–509.

    Article  CAS  PubMed  Google Scholar 

  32. Danuza M, Viviane R, João A, Lidia M, Fernando T. Xylose fermentation by Saccharomyces cerevisiae: challenges and prospects. Int J Mol Sci. 2016;17(3):207.

    Article  Google Scholar 

  33. Groot N, Ventura S. Effect of temperature on protein quality in bacterial inclusion bodies. FEBS Lett. 2006;580:6471–6.

    Article  PubMed  Google Scholar 

  34. Waldo GS, Standish BM, Berendzen J, Terwilliger TC. Rapid protein-folding assay using green fluorescent protein. Nat Biotechnol. 1999;17:691–5.

    Article  CAS  PubMed  Google Scholar 

  35. Feilmeier BJ, Iseminger G, Schroeder D, Webber H, Phillips GJ. Green fluorescent protein functions as a reporter for protein localization in Escherichia coli. J Bacteriol. 2000;182:4068–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Hjelm A, Schlegel S, Baumgarten T, Klepsch M, Wickstrom D, Drew D, de Gier J-W. Optimizing E. coli-based membrane protein production using Lemo21(DE3) and GFP-fusions. Methods Mol Biol (Clifton, NJ). 2013;1033:381–400.

    Article  CAS  Google Scholar 

  37. Toivari MH, Aristidou A, Ruohonen L, Penttilä M. Conversion of xylose to ethanol by recombinant Saccharomyces cerevisiae: importance of xylulokinase (XKS1) and oxygen availability. Metab Eng. 2001;3:236–49.

    Article  CAS  PubMed  Google Scholar 

  38. Barber RD. GAPDH as a housekeeping gene: analysis of GAPDH mRNA expression in a panel of 72 human tissues. Physiol Genomics. 2005;21:389–95.

    Article  CAS  PubMed  Google Scholar 

  39. Marko K, Harhangi HR, Kristin SA, Winkler AA, Jetten MSM, De LWTAM, Den RJJJ, Op DCHJM, Van DJP, Pronk JT. High-level functional expression of a fungal xylose isomerase: the key to efficient ethanolic fermentation of xylose by Saccharomyces cerevisiae? Fems Yeast Research. 2004;4(1):69–78.

  40. Maris A, Winkler AA, Kuyper M, Laat W, Pronk JT. Development of efficient xylose fermentation in Saccharomyces cerevisiae: xylose isomerase as a key component. Adv Biochem Eng Biotechnol. 2007;108:179.

    PubMed  Google Scholar 

  41. Ledesma-Amaro R, Lazar Z, Rakicka M, Guo Z, Fouchard F, Coq CL, Nicaud JM. Metabolic engineering of Yarrowia lipolytica to produce chemicals and fuels from xylose. Metab Eng. 2016;38:115–24.

    Article  CAS  PubMed  Google Scholar 

  42. Inoue H, Hashimoto S, Matsushika A, Watanabe S, Sawayama S. Breeding of a xylose-fermenting hybrid strain by mating genetically engineered haploid strains derived from industrial Saccharomyces cerevisiae. J Ind Microbiol Biotechnol. 2014;41:1773–81.

    Article  CAS  PubMed  Google Scholar 

  43. Zhu YX, Zhang JT, Zhu L, Jia ZF, Li Q, Xiao W, Cao LM. Minimize the xylitol production in Saccharomyces cerevisiae by balancing the xylose redox metabolic pathway. Front Bioeng Biotechnol. 2021;9: 639595.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Kildegaard KR, Wang Z, Chen Y, Nielsen J, Borodina I. Production of 3-hydroxypropionic acid from glucose and xylose by metabolically engineered Saccharomyces cerevisiae. Metab Eng Commun. 2015;2:132–6.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Wu YF, Xu S, Gao X, Li M, Li DH, Lu WY. Enhanced protopanaxadiol production from xylose by engineered Yarrowia lipolytica. Microb Cell Fact. 2019;18:83.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Zbab C, Yzab C, Kai Z, Yfab C, Xcab C, Mlab C, Lyab C. High-value utilization of the waste hydrolysate of Dioscorea zingiberensis for docosahexaenoic acid production in Schizochytrium sp. Biores Technol. 2021;336: 125305.

    Article  Google Scholar 

  47. Liang Y, Sarkany N, Cui Y, Yesuf J, Trushenski J, Blackburn JW. Use of sweet sorghum juice for lipid production by Schizochytrium limacinum SR21. Biores Technol. 2010;101:3623–7.

    Article  CAS  Google Scholar 

  48. Unagul P, Assantachai C, Phadungruengluij S, Suphantharika M, Tanticharoen M, Verduyn C. Coconut water as a medium additive for the production of docosahexaenoic acid (C22:6 n3) by Schizochytrium mangrovei Sk-02. Biores Technol. 2007;98:281–7.

    Article  CAS  Google Scholar 

  49. Gonzalez-Pino MJ, Rangel-Aldao R, Slezynger TC. Expression of alpha- and beta-tubulin genes during growth of Trypanosoma cruzi epimastigotes. DNA Cell Biol. 1999;18:449.

    Article  CAS  PubMed  Google Scholar 

  50. Moes CJ, Pretorius IS, Zyl WHV. Cloning and expression of the Clostridium thermosulfurogenes D-xylose isomerase gene (xyLA) in Saccharomyces cerevisiae. Biotech Lett. 1996;18:269–74.

    Article  CAS  Google Scholar 

  51. Parachin NS, Gorwa-Grauslund MF. Isolation of xylose isomerases by sequence- and function-based screening from a soil metagenomic library. Biotechnol Biofuels. 2011;4:9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Hang Z, Cheng JS, Wang BL, Fink GR, Stephanopoulos G. Xylose isomerase overexpression along with engineering of the pentose phosphate pathway and evolutionary engineering enable rapid xylose utilization and ethanol production by Saccharomyces cerevisiae. Metab Eng. 2012;14:611–22.

    Article  Google Scholar 

  53. Juntila DJ. Genetic modification of the thraustochytrid Aurantiochytrium sp. 18W–13a for cellobiose utilization by secretory expression of β-glucosidase from Aspergillus aculeatus. Algal Res. 2019;40:101503.

    Article  Google Scholar 

  54. Ren LJ, Sun LN, Zhuang XY, Qu L, Ji XJ, Huang H. Regulation of docosahexaenoic acid production by Schizochytrium sp.: effect of nitrogen addition. Bioprocess Biosyst Eng. 2014;37:865–72.

    Article  CAS  PubMed  Google Scholar 

  55. Sun XM, Zhang ZX, Wang LR, Wang JG, Yang S. Downregulation of T7 RNA polymerase transcription enhances pET-based recombinant protein production via suppression of Escherichia coli BL21 (DE3) programmed cell death. Biotechnol Bioeng. 2021;118:153–63.

    Article  CAS  PubMed  Google Scholar 

  56. Sun XM, Ren LJ, Bi ZQ, Ji XJ, Zhao QY, He H. Adaptive evolution of microalgae Schizochytrium sp. under high salinity stress to alleviate oxidative damage and improve lipid biosynthesis. Bioresour Technol. 2018;267:438.

    Article  CAS  PubMed  Google Scholar 

  57. Huang PW, Xu YS, Sun XM, Shi TQ, Gu Y, Ye C, Huang H. Development of an efficient gene editing tool in Schizochytrium sp. and improving its lipid and terpenoid biosynthesis. Front Nutr. 2021;8:795651.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Yuan Y, Zhai R, Li Y, Chen X, Jin M. Correction to: Developing fast enzyme recycling strategy through elucidating enzyme adsorption kinetics on alkali and acid pretreated corn stover. Biotechnol Biofuels. 2018;11:327.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Ren LJ, Huang H, Xiao AH, Lian M, Jin LJ, Ji XJ. Enhanced docosahexaenoic acid production by reinforcing acetyl-CoA and NADPH supply in Schizochytrium sp. HX-308. Bioprocess Biosyst Eng. 2009;32:837.

    Article  CAS  PubMed  Google Scholar 

  60. Xu P, Li L, Zhang F, Stephanopoulos G, Koffas M. Improving fatty acids production by engineering dynamic pathway regulation and metabolic control. Proc Natl Acad Sci USA. 2014;111:11299.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


This work was supported by the National Natural Science Foundation of China (No. 22038007, No. 21908112), and the Nature Science Foundation of Jiangsu Province (No. BK20190706).

Author information

Authors and Affiliations



XMS and LRW planned the experiments. LRW, ZXZ, FTN, JL, PWH, and WM performed the experiments. LRW and FTN analyzed the data. LRW and ZXZ wrote the paper. XMS and QYZ revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xiao-Man Sun.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1: Fig S1

. Xylose content of the crude enzyme solution system before and after incubation. XK7938: xylulose kinase from Schizochytrium sp., XK(E): xylulose kinase from Escherichia coli, XK(S): xylulose kinase from Saccharomyces cerevisiae. Data represent the Kmean±standard deviation (n=3). Table S1. Plasmids used in this work. Table S2. Primers used in this work. All primers are synthesized by Synbio-tech Co., Ltd., China. Table S3. Native xylose metabolic pathway genes of Schizochytrium sp. HX-308. Table S4. Xylose metabolism pathway genes from HX-308 after codon optimization in E. coli.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, LR., Zhang, ZX., Nong, FT. et al. Engineering the xylose metabolism in Schizochytrium sp. to improve the utilization of lignocellulose. Biotechnol Biofuels 15, 114 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: