Skip to main content

MdFRK2-mediated sugar metabolism accelerates cellulose accumulation in apple and poplar

Abstract

Background

Cellulose is not only a common component in vascular plants, but also has great economic benefits for paper, wood, and industrial products. In addition, its biosynthesis is highly regulated by carbohydrate metabolism and allocation in plant. MdFRK2, which encodes a key fructokinase (FRK) in apple, showed especially high affinity to fructose and regulated carbohydrate metabolism.

Results

It was observed that overexpression of MdFRK2 in apple decreased sucrose (Suc) and fructose (Fru) with augmented FRK activity in stems, and caused the alterations of many phenotypic traits that include increased cellulose content and an increase in thickness of the phloem region. To further investigate the involved mechanisms, we generated FRK2-OE poplar lines OE#1, OE#4 and OE#9 and discovered (1) that overexpression of MdFRK2 resulted in the huge increased cellulose level by shifting the fructose 6-phosphate or glucose 6-phsophate towards UDPG formation, (2) a direct metabolic pathway for the biosynthesis of cellulose is that increased cleavage of Suc into UDP-glucose (UDPG) for cellulose synthesis via the increased sucrose synthase (SUSY) activity and transcript levels of PtrSUSY1, (3) that the increased FRK activity increases the sink strength overall so there is more carbohydrate available to fuel increased cambial activity and that resulted in more secondary phloem. These results demonstrated that MdFRK2 overexpression would significantly changes the photosynthetic carbon flux from sucrose and hexose to UDPG for increased cellulose synthesis.

Conclusions

The present data indicated that MdFRK2 overexpression in apple and poplar changes the photosynthetic carbon flux from sucrose and hexose to UDPG for stem cellulose synthesis. A strategy is proposed to increase cellulose production by regulating sugar metabolism as a whole.

Background

Cell walls are one of the most important features of plant cells. About 70% of photosynthetic products are converted to polymers that accumulate in cell walls, and cell wall biomass is the most abundant renewable resource on earth [1]. Since the non-renewable nature of fossil energy has become a major challenge to sustainable human development, exploring how to convert cell wall biomass into usable energy is of great importance for human society.

The cell wall structure is complex and primarily consists of various polysaccharides, such as cellulose, hemicellulose, and pectin [2]. Cellulose, the most abundant renewable polymer in the biosphere [3], is already an important source of raw materials for textiles and paper products and has great potential for use in renewable biofuels or food production. Cellulose is a major component of the cell wall, where it exists in the form of microfibrils. Cellulose is composed of β-1,4-glucan chains linked by intermolecular hydrogen bonds and is made by cellulose synthase A (CESA) complexes (CSCs) at the plasma membrane that are believed to form sixfold-symmetrical rosettes [4]. CesA encodes a glycosyltransferase that belongs to the GT2 family and has been reported to have key roles in cellulose biosynthesis [5]. Indirect evidence suggested that CESA-4 and CESA-8 are associated with secondary wall cellulose synthesis in Arabidopsis [6]. Overexpression of a CesA mutant gene, CesA7fra5, in Arabidopsis led to a strong cellulose deficiency [7]. In addition, regulation of cellulose synthesis has a direct impact on plant growth and development. It is remarkable to note that reduced cellulose synthesis affects the softening of apples [8].

It is clear from studies of cellulose regulation that carbohydrate supply are essential for normal cellulose synthesis [9, 10]. Phloem is a transport system connecting source and sink, and it is required for the distribution of photosynthetic products (sugars) from source to sink organs [11]. In higher plants, photosynthetically active leaves export a large amount of carbon assimilates, primarily in the form of sucrose (Suc), via the phloem in support of sink development. Within the cytosol of sink cells, sucrose synthase (SUSY, EC 2.4.1.13) cleaves Suc to UDP-glucose (UDPG), which is the sole precursor for cellulose [12, 13]. SUSY activity has an important role in sink tissue metabolism and is a major area of interest within the field of cellulose biosynthesis [14, 15]. Another pathway of UDPG production from Suc is via neutral invertases (NINV, EC 3.2.1.26), hexose kinase enzymes, hexose phosphate interconversion enzymes and UDP-glucose pyrophosphorylase (UGPase, EC 2.7.7.9). Rende et al. [16] suggested that suppression of NINV resulted in a 38–55% reduction in NINV activity, a lower UDPG level, and a 9–13% decrease in cellulose content in poplar wood tissue.

Transport of Suc to the sink occurs via the phloem, and then in sink tissues NINV converts Suc into Glc and Fru, or SUSY converts Suc into UDPG and Fru [13, 17]. The resulting Glc and Fru can be phosphorylated by hexokinase (HxK, EC 2.7.1.1) and fructokinase (FRK, EC 2.7.1.4) to generate hexose phosphates, which are then converted to UDPG for cellulose biosynthesis [18]. FRK is the main fructose-phosphorylating enzymes in planta [19], and that FRK2 orthologs play a vital role in controlling Fru utilization and carbon partition in hybrid aspen (Populus tremula × tremuloides) [20], apple (Malus × domestica.) [17], peach (Prunus persica) [21], and tomato (Solanum lycopersicon cv. MP1) [22]. Roach et al. showed that RNA interference of FRK2 (RNAi-FRK2) in hybrid aspen decreased cell wall fiber thickness and the proportion of cellulose in the cell wall, implying that FRK2 facilitates carbon allocation to cellulose in the wood [20]. Apple and other Rosaceae fruit trees transport more of their photosynthates as fructose, with more than 80% of the total carbon flux transporting as Fru into sink cells [23]. Our previous research demonstrated, compared with other plants, that MdFRK2 is highly expressed in apple sink organs and that the MdFRK2 protein not only has a high-affinity for Fru but also a high enzymatic activity [17, 24], and overexpression MdFRK2 in apple had a significant change in carbohydrate metabolism in the leaves [17]. Therefore, we inferred that augmented MdFRK2 in plant sink cells might play a strong role in regulating sucrose metabolism and push more carbon forward for cellulose biosynthesis.

In this work, we observed that the transgenic apple with overexpressing MdFRK2, as our previous report [17], had significantly increased cellulose content and thickened the secondary phloem. To further investigate the involved mechanisms, we generated MdFRK2-OE poplar lines. The increased Fru phosphorylation in MdFRK2-transgenic poplar lines could promote the Suc cleavage via SUSY, and then generated more UDPG, which is the sole precursor for cellulose biosynthesis. MdFRK2-OE in poplar resulted in an increase in thickness of the secondary phloem region, but had no effects on stem height or diameter. We proposed a model that serves as a roadmap for future work to better understand the molecular network responsible for cellulose biosynthesis. These results demonstrated that MdFRK2 overexpression in apple and poplar changes the photosynthetic carbon flux from Suc and hexose to UDPG for cellulose synthesis and will contribute to research on cellulose biosynthesis.

Results

Overexpression of MdFRK2 increased cellulose level in transgenic apple

We generated MdFRK2-OE transgenic apple lines (OE-4 and OE-9) using the cauliflower mosaic virus 35S (CaMV 35S) promoter (Fig. 1a). Compared to wild type (WT), MdFRK2 expression in the stems increased by 12.55- and 13.79-fold in lines OE-4 and OE-9, respectively (Fig. 1b). The activity of FRK was markedly increased in the OE lines (Fig. 1c). Cellulose contents of MdFRK2-OE transgenic apple stems were increased 0.12- and 0.15-fold, respectively, over WT (Fig. 1d). Light microscopy revealed that the MdFRK2-OE lines had significantly increased in thickness of the cambium and secondary phloem region compared to WT (Fig. 1e). These results showed that MdFRK2 is a major player in cellulose synthesis in apple plants.

Fig. 1
figure 1

Molecular and microscopic characterization in stems of wild type (WT) and transgenic apple plants (OE-4 and OE-9). a Phenotypes of WT and transgenic lines (OE-4 and OE-9). b, c Quantitative RT-PCR of MdFRK2 expression and FRK (fructokinase) activity in stems from WT, OE-4 and OE-9. In (b), the MdActin gene was used as internal control. d Content of cellulose in stems from WT and OE-4 and OE-9. Bars represent the mean value ± SE (n ≥ 3). The asterisk indicates significant differences at P ≤ 0.05. e Microstructure of petiole sections in WT and transgenic line OE-4. 1: pith; 2: parenchyma; 3: bast fiber; 4: phloem; 5: xylem. Scale bars = 200 μm

Alteration of soluble carbohydrate concentration in MdFRK2 transgenic apple

To decipher how changes in the level of soluble carbohydrates affects cellulose synthesis with MdFRK2 overexpression, the Fru concentrations in OE-4 and OE-9 transgenic lines were measured. Fru levels decreased by 15.00% in OE-4 and 18.28% in OE-9 (Fig. 2). It is worth mentioning here that the Suc concentration in the OE-4 transgenic line decreased significantly, by 7.50%, while Glc concentrations significantly decreased by 18.68% and 15.26%, respectively (Fig. 2). These results indicated that the links between MdFRK2 and cellulose synthesis might be related to sugar metabolism in the sink.

Fig. 2
figure 2

Sugar concentrations in the stems of the transgenic apple plants overexpressing MdFRK2. Sucrose, fructose, glucose and galactose concentrations were measured in wild-type WT and transgenic OE-4 and OE-9 lines. Bars represent the mean value ± SE (n ≥ 3). Asterisks indicate significant differences at P ≤ 0.05

Heterologous expression of MdFRK2 in poplar plants

To further understand the cellulose synthesis-related functions of MdFRK2 in stems, we heterologously expressed it in poplar (Populus clone 717) under the expression of CaMV promoter. Poplar is rich in cellulose and is an important source material for everyday products such as cloth, paper and biofuels [25]. Analysis of the DNA and mRNA levels revealed that three heterologous lines (OE#1, OE#4 and OE#9) (Fig. 3a) exhibited increases in MdFRK2 transcript levels in stems relative to the levels in WT controls (Fig. 3b, c). All three lines also exhibited significantly increased enzyme activity. Interestingly, the FRK activity of mature leaves did not differ between the transgenic and wild-type plants (Additional file 1: Fig. S1). Lines OE#1, OE#4 and OE#9 displayed 0.70, 0.57 and 0.49-fold increases, respectively, in FRK activity in the stems, relative to the levels of the untransformed WT controls (Fig. 3d). Regardless of these observed differences, there were no significant changes in stem height and diameter in transgenic lines (Additional file 3: Table S1). These data showed that apple FRK2 functioned very specifically in sugar metabolism in poplar.

Fig. 3
figure 3

Characterization of wild type and transgenic poplar overexpressing apple FRK2. a Phenotype of wild type (WT, Populus clone 717) and transgenic lines (OE#1, OE#4 and OE#9). b DNA levels in the stems of WT and transgenic lines OE#1, OE#4 and OE#9. M, marker; P, plasmid. c, d Quantitative RT-PCR of MdFRK2 expression and FRK (fructokinase) activity in stems from WT and transgenic lines OE#1, OE#4 and OE#9. In (c), the PtrActin gene was used as internal control. Bars represent the mean value ± SE (n ≥ 3). An asterisk indicates a significant difference at P ≤ 0.05

Heterologous expression of MdFRK2 altered soluble sugar concentration of poplar plants

To examine the effect of overexpression of MdFRK2 on carbohydrate metabolism in poplar, we determined sugar concentrations in transgenic lines using gas chromatography mass spectrometry GC/MS (Fig. 4). The high MdFRK2 transcript levels in stems resulted in lower concentrations of Fru in the three transgenic lines, with levels reduced to 13.87% of the control level for OE#1, 15.76% for OE#4 and 15.80% for OE#9. The concentrations of Suc in these three transgenic lines were decreased by 17.65%, 15.16% and 15.84%, respectively. In the transgenic lines, the concentration of Glc was also decreased. These data demonstrated that MdFRK2 modulate sucrose and hexose metabolism in a heterologous species.

Fig. 4
figure 4

Sugar concentrations in the stems of transgenic poplar plants stems overexpressing MdFRK2. Sucrose, fructose, glucose and galactose concentrations were measured in wild type and transgenic OE#1, OE#4 and OE#9. Bars represent the mean value ± SE (n ≥ 3). An asterisk indicates significant differences at P ≤ 0.05

Carbohydrate metabolism pathway in stems of transgenic poplar

To determine why the sugar concentrations changed in the transgenic poplar lines expressing MdFRK2, enzyme activity and expression of genes related to sugar metabolism were assessed (Fig. 5a). The activity of cell wall invertase (CWINV) was not statistically different between the wild type and transgenic lines. However, the activities of neutral invertase (NINV) and sucrose synthase (SUSY), both of which are related to Suc breakdown, were significantly increased (Additional file 2: Figure S2). A similar pattern was observed for HxK and FRK activities. These results showed that the breakdown of Suc and Fru in sink cells is through the activities of SUSY and FRK, respectively. The changes in Suc and Fru were similar to the pattern of enzyme activities involved in Suc and Fru metabolism.

Fig. 5
figure 5

Enzyme activities and the expression of genes related to carbohydrate metabolism in the transgenic poplars (OE#1, OE#4 and OE#9) stems overexpressing MdFRK2. Activities of enzymes (a) and relative mRNA expression of genes (b) related to sugar metabolism in stems of wild type (WT, Populus clone 717) and transgenic lines (OE#1, OE#4 and OE#9). CWINV, cell wall invertase; NINV, neutral invertase; SUSY, sucrose synthase; HxK, hexokinase. Quantitative RT-PCR was performed with gene-specific primers using PtrActin as an internal control. Bars represent the mean value ± SE (n ≥ 3). An asterisk indicates significant differences at P ≤ 0.05

In addition, the transcript abundance of genes encoding these enzymes were investigated (Fig. 5b). Gene involved in Suc degradation (PtrSUSY1) was significantly upregulated in the transgenic poplar. In accord with the reduced Glc concentrations (Fig. 4), Transcripts for PtrHxK1, were increased. However, the expression level of PtrCWINV2 was unchanged relative to control levels in transgenic lines. Taken together, these findings further suggested that the decreased Suc and Fru concentrations in MdFRK2-transgenic Populus were due to increased cleavage of Suc and Fru phosphorylation into hexose phosphates via the increased SUSY and FRK activities and transcript levels of PtrSUSY1 and MdFRK2, respectively.

Heterologous expression of MdFRK2 accelerated UDPG accumulation of poplar plants

SUSY directly produces UDPG, which is the substrate for cellulose synthesis in sink organs. To confirm that alteration of UDPG in the transgenic lines was caused by increased SUSY or FRK activity, the UDPG concentration was detected (Fig. 6). In poplar overexpressing MdFRK2, the concentrations of UDPG were significantly increased, by 0.95-, 0.70- and 0.69-fold, respectively. Furthermore, the levels of fructose 6-phosphate (F6P), glucose 6-phsophate (G6P) and glucose 1-phsophate (G1P) were increased greatly in all three OE lines relative to the levels in WT. These results suggested that MdFRK2 overexpression could accelerate UDPG accumulation and that the capacity for cellulose synthesis via UDPG would be increased.

Fig. 6
figure 6

Hexose phosphate concentrations in stems of transgenic poplars (OE#1, OE#4 and OE#9) overexpressing MdFRK2. F6P, fructose 6-phosphate; G6P, glucose 6-phsophate; G1P, glucose 1-phsophate; UDPG, UDP-glucose. Bars represent the mean value ± SE (n ≥ 3). An asterisk indicates significant differences at P ≤ 0.05

Heterologous expression of MdFRK2 accelerated cellulose accumulation in Populus plants

After seeing increased UDPG levels in stems of these three transgenic lines with up-regulated MdFRK2 expression, the cellulose content of all transgenic lines were measured (Fig. 7a). The contents of cellulose increased 0.49-, 0.21- and 0.24-fold in the three transgenic lines as compared with WT, respectively. Accordingly, the three transgenic poplar lines contained increased hemicellulose contents, by 6.90% in OE#1, 15.64% in OE#4 and 14.84% in OE#9 (Fig. 7b). The pectin contents in the transgenic poplar stems decreased 0.22-, 0.19- and 0.15-fold compared to WT (Fig. 7c), whereas the lignin content showed no significant changes (Fig. 7d).

Fig. 7
figure 7

Contents of cellulose, hemicellulose, pectin, lignin and the expression of genes that are involved in the biosynthesis of cellulose (PtrCESA4, and PtrCESA8), hemicellulose (PtrFRA8, and PtrIRX9), and lignin (PtrPAL1, and PtrC4H1) in stems of the transgenic poplar lines (OE#1, OE#4 and OE#9) overexpressing MdFRK2. Contents of cellulose, hemicellulose, pectin, lignin (ad) and relative mRNA expression of genes (e) that are involved in the biosynthesis of cellulose (PtrCESA4, and PtrCESA8), hemicellulose (PtrFRA8, and PtrIRX9), and lignin (PtrPAL1, and PtrC4H1) in stems of wild type (WT, Populus clone 717) and transgenic lines (OE#1, OE#4 and OE#9). Quantitative RT-PCR was performed with gene-specific primers using PtrActin as an internal control. Bars represent the mean value ± SE (n ≥ 3). An asterisk indicates significant differences at P ≤ 0.05

To further identify key genes contributing to changes in cellulose levels in the transgenic lines, the expression levels of genes related to these traits were measured using qRT-PCR. The results suggested that the expression levels of genes involved in the biosynthesis of cellulose (PtrCesA4 and PtrCesA8), were highly up-regulated in OE#1, OE#4 and OE#9 transgenic lines compared to in wild type (Fig. 7e). However, the expression levels of hemicellulose (PtrFRA8 and PtrIRX9) and lignin biosynthetic genes (PtrPAL1 and PtrC4H1) did not obviously differ. These results demonstrated that cellulose is a major sink for FRK2-metabolized carbon.

Heterologous expression of MdFRK2 changed secondary phloem in Populus plants

As reported previously [26], the cambium, is a meristem between xylem and phloem, produces secondary phloem. To test whether the increased cellulose level altered the secondary phloem of transgenic lines, safranin and fast green were used to stain cellulose and lignin, respectively. Examination by light microscopy revealed that OE#1, OE#4 and OE#9 transgenic lines had significantly increased in thickness of the   secondary phloem region (Fig. 8a–d), which was increased by 8.33%, 10.96% and 12.94% compared to wild type, respectively (Fig. 8e). These results indicated that increased FRK activity increases the sink strength overall so there is more carbohydrate available to fuel increased cambial activity and that resulted in more secondary phloem.

Fig. 8
figure 8

Secondary phloem of the transgenic poplar lines (OE#1, OE#4 and OE#9) overexpressing MdFRK2. a–d Cellulose (green color) in stem sections stained with safranin and fast green. Scale bars = 200 μm. a Wild type (WT, Populus clone 717); b: transgenic line OE#1; c: transgenic line OE#4; d: transgenic line OE#9. ph: phloem (black arrow); xy: xylem (black arrow); C: cambium (black arrow); red lines: the cambium and secondary phloem region. An increase in thickness e and area f of the secondary phloem region in transgenic lines (OE#1, OE#4 and OE#9) in comparison with WT. Bars represent the mean value ± SE (n ≥ 3). The asterisk indicates a significant difference at P ≤ 0.05

Discussion

MdFRK2 is essential for cellulose biosynthesis in apple plants

Photosynthetically active leaves produce many forms of carbohydrates, such as Suc, which serve as energy and carbon sources for cellulose synthesis in the sink of most plant species. There have been many reports suggested that accelerated Suc decomposition leads to increased cellulose content in SUSY overexpression plants [27,28,29]. More recently, FRK2 is generally considered a main regulator that functions in carbohydrate metabolism in the leaves [17]. However, it remains unknown how MdFRK2 plays a role in regulating sugar metabolism for cellulose biosynthesis.

In this report, two MdFRK2-OE transgenic apple plants, OE-4 and OE-9, were observed (Fig. 1a). The significant increases in MdFRK2 transcript levels and FRK enzyme activity in stems suggest that MdFRK2 may have a special function in stem development. Similarly, the influence of the FRK2 gene on cell-wall biosynthesis was reported in aspen, suggesting that the reduction in FRK2 activity in FRK2-RNAi hybrid aspen primarily affected cellulose [20]. Thus, we hypothesized that apple FRK2 plays direct roles in cellulose synthesis. As we had expected, elevated MdFRK2 transcript levels significantly increased the cellulose content and resulted in an increase in thickness of the cambium and secondary phloem region in the OE lines. Our data demonstrated that MdFRK2 aids cellulose synthesis in apple plants. In addition, the previous studies have reported that hybrid aspen FRK2 could effectively control carbon flux into cell walls for cellulose biosynthesis [20]. In the present study, overexpression of MdFRK2 resulted in significantly lower concentrations of soluble sugar in the MdFRK2-OE apple lines (Fig. 2), which is similar to the results from previous studies of apple and aspen as discussed. In light of these results, we speculated that the links between MdFRK2 and cellulose synthesis might be related to sugar metabolism in the sink.

Cellulose is a major sink for FRK2-metabolized carbon

Poplar is more richness in cellulose than apple, and poplar is an important source for the raw materials used in everyday products such as cloth, paper and biofuels [30]. To study the functions of MdFRK2 from apple, we heterologously expressed MdFRK2 in Populus and obtained three lines overexpressing the fructokinase gene (Fig. 3a, b). As expected, the mRNA levels and enzyme activity in stems increased significantly (Fig. 3c, d, Additional file 1: Fig. S1) due to MdFRK2 overexpression. This result was expected, as it was previously reported that overexpression of MdFRK2 in apple significantly increased MdFRK2 transcription and FRK enzyme activity [17]. Surprisingly, this increased FRK2 activity resulted in altered cell wall cellulose content, with increases of 21% to 24% over control levels (Fig. 7a), without affecting plant growth (Additional file 3: Table S1). Taken together, these results are consistent with previous reports that the function of apple FRK2 is very specific for cellulose biosynthesis in poplar.

It is particularly notable that overexpression of MdFRK2 in transgenic poplars caused significant decreases in soluble sugar and increases in hexose-phosphates levels, especially UDPG (Figs. 4, 6). Moreover, the key enzymes and the expression levels of main genes related to soluble sugar metabolism also showed a significant increase (Fig. 5). Thus, we suggest that cellulose synthesis might be related to modulation of Suc and hexose metabolism via MdFRK2. Here we present lines of direct evidence to explain which mechanisms might be involved in this. First, cellulose is generated from the precursor UDP-glucose (UDPG), which can be derived from the cleavage of Suc by SUSY to directly yield UDPG and Fru [31]. In this context, our data are in line with previous reports that showed that the increased cellulose in MdFRK2-transgenic poplars was due to increased cleavage of Suc into UDPG via the increased SUSY activity and transcript levels of PtrSUSY1 (Figs. 4, 5 and 6, Additional file 2: Fig. S2). Similar results have been found in cotton [15], hybrid aspen [20], and wheat (Triticum aestivum L.) [32], with overexpression of SUSY being associated with increased cellulose synthesis. These data showed a direct metabolic pathway for the biosynthesis of UDPG to cellulose.

Second, the resulting Fru and Glc can be readily phosphorylated by the increased activities of FRK and HxK enzymes to generate hexose-phosphates that are then converted to UDPG for cellulose synthesis. Here, MdFRK2-overexpressing transgenic poplar showed that transcripts of other FRKs were not influenced (Fig. 5b), suggesting that the increased FRK activity is, therefore, due to an increase in MdFRK2 activity. This increase correlated with decreased Fru and increased F6P, UDPG and cellulose levels (Figs. 4, 6, 7a). It is, therefore, plausible that high FRK activity and the resulting decrease in Fru levels could lead to an increase in cellulose by shifting the F6P towards UDPG formation. A similar conclusion comes from a study of FRK2-RNAi in hybrid aspen [20]. Furthermore, it is important to note that G6P, a hexose-phosphate, is an important intermediate in cellulose synthesis [33]. In MdFRK2-transgenic poplar plants, the resulting Glc in carbon metabolism will be phosphorylated by the increased HxK enzyme activity (correlated to increased transcript levels of PtrHxK1) to produce G6P (Fig. 4, 5 and 6). G6P has three possible fates, namely, (1) entering glycolysis/TCA cycle, (2) being used for starch synthesis, (3) being converted to UDPG for cellulose synthesis. These findings provide novel insights into the relationship between MdFRK2 and cellulose synthesis [13, 34].

Besides cellulose, hemicelluloses and pectins influence both crystallinity and network connectivity of cellulose microfibrils in cell walls [5]. In this study, our data showed a significant increase in hemicellulose content in transgenic poplar and a reduction in the pectin concentration when compared with WT (Fig. 7b, c). The data support the idea that increased SUSY activity in transgenic tobacco could distribute more carbon to hemicellulose synthesis, but less to pectins, as previously described [31]. A possible explanation for this might be that the increased MdFRK2 activity provided more carbohydrate substrates to hemicellulose, which is in line with the described mechanism that links MdFRK2 with regulated cellulose or hemicellulose production and might be related to its ability to regulate carbon distribution. It is interesting to note that the expression of genes (PtrFRA8 and PtrIRX9) related to hemicellulose [35] were unchanged in relative to control levels in transgenic poplars (Fig. 7e).

Potential basis for cellulose biosynthesis in apple and poplar

A multiple transmembrane spanning protein, CESA, directly influences cellulose content. UDPG is the substrate of the CESAs [36]. Indirect evidence suggests that CESA4 and CESA8 seem to be absolutely necessary for secondary wall cellulose synthesis [6]. In this study, the expression of genes involved in cellulose biosynthesis (PtrCesA4 and PtrCesA8) were significantly increased in MdFRK2-overexpressing transgenic poplars (Fig. 7e). Based on these data, we believe that there is a positive correlation between cellulose synthase and CESA gene expression. The importance of plant FRK2 for the development of vascular tissues has been documented in sinks in hybrid aspen [20] and tomato [37]. In our report, overexpression of MdFRK2 resulted in augmented cellulose levels in transgenic poplar and an increase in thickness of the secondary phloem region, which increased 10.96% to 12.94% compared to wild type (Figs. 7a, 8ae). Together with the increased cellulose levels, these results supported the view that with increased FRK activity increases the sink strength overall so there is more carbohydrate available to fuel increased cambial activity and that resulted in more secondary phloem, suggesting that MdFRK2 plays a general role in carbon partitioning to UDPG for cellulose synthesis.

Conclusion

The results reported here suggest that (1) overexpression of MdFRK2 could increase carbon flux to UDPG for cellulose biosynthesis, which can then be used for phloem development in apple and poplar plants; (2) a direct metabolic pathway for the biosynthesis of UDPG is through increased cleavage of Suc into UDPG via the increased SUSY activity and transcript levels of SUSY1; (3) an alternative pathway of UDPG production from Suc via NINV; and (4) overexpression of MdFRK2 results in increased cellulose content by shifting F6P or G6P towards UDPG formation (Fig. 9). These results provide insights into the role of MdFRK2 in controlling cellulose or hemicellulose production that might be related to its ability to regulate sugar metabolism.

Fig. 9
figure 9

Proposed model of MdFRK2-regulated cellulose synthesis in apple and poplar. In stems, a direct metabolic pathway for the biosynthesis of cellulose is through the increased cleavage of Suc into UDPG via increased SUSY activity, which is an alternative pathway for UDPG production from Suc via NINV. One observation of this study is that overexpression of MdFRK2 resulted in increased cellulose content by shifting F6P or G6P towards UDPG formation. With increased FRK activity increases the sink strength overall so there is more carbohydrate available to fuel increased cambial activity and that resulted in more secondary phloem. Blue boxes indicate a reduction in soluble sugar concentrations. Yellow boxes indicate increases in soluble sugar or hexose phosphate concentrations, enzyme activity or relative mRNA expression of genes. White boxes, no change in relative mRNA expression of genes that are involved in the biosynthesis of lignin (PAL1, and C4H1). Purple boxes indicate the relative mRNA expression of MdFRK2

Materials and methods

Plant materials

Tissue-cultured WT and and MdFRK2-transformed ‘GL-3’ apple plantlets were grown on Murashige and Skoog (MS) medium supplemented with 0.2 mg L−1 IAA, 0.3 mg L−1 6-BA and 25 mg · L − 1 kanamycin for 4 weeks. They were then rooted in MS medium (MS + 0.5 mg IBA and 0.5 mg IAA) for 2 months.

Untransformed poplar WT (Populus clone 717) and MdFRK2-transformed ‘Populus clone 717’ poplar plantlets were grown on MS medium containing 0.25 g L−1 MES, 0.1 g L−1 Inositol and 0.3 g L−1 l-glutamine for 2 months.

After rooting, all genotypes were transferred to a culture room at Northwest A&F University, Yangling, Shaanxi, China, under conditions previously reported by Wang et al. [38]. After these plants had grown for 2 months, stems were collected and immediately frozen in liquid nitrogen and stored at −80 ºC.

Cloning of MdFRK2

The MdFRK2 sequence (MD04G1042400) was retrieved from the Malus Genome Database (http://www.rosaceae.org). Specific primers were designed for gene cloning (Additional file 2: Table S2). Total RNA was extracted from young fruits of ‘Gala’ apple, and cDNA was synthesized using PrimeScript™ II Reverse Transcriptase (Takara, Dalian, China).

Vector constructs and Populus transformation

The coding region of MdFRK2 was cloned into the gateway vector PGWB402 with the CaMV 35S promoter. The recombinant plasmid was transformed into Agrobacterium tumefaciens strain EHA105 for transformation. The transformation of Populus was done according to the procedure of Dai et al. [39]. Afterward, Popular plants were grown at 24 °C under a 15-h photoperiod supplemented with fluorescent light at 60 μmol m−2 s−1.

Transgenic plants were screen by kanamycin resistance and PCR analysis of extracted DNA. From the four overexpression transgenic lines obtained, we selected three for further analysis. Untransformed WT served as the control plants.

Sugar concentration measurement

As recently described [40], soluble sugars and hexose phosphates were extracted and then derivatized with methoxyamine hydrochloride and Nmethyl-N-trimethylsilyl-trifluoroacetamide. After derivatization, the metabolites were analyzed using a Shimadzu GC/MS-2010SE (Shimadzu Corporation, Tokyo, Japan).

Enzyme assay and expression analysis

CWINV, SUSY, NINV, FRK, and HxK in steam samples were extracted as described [13]. Soluble proteins were measured using Coomassie blue, and enzyme activities were expressed on a protein basis.

Total RNA from frozen tissues was extracted with RNAprep plant kit (CWBIO, Beijing, China), and cDNA was synthesized using PrimeScript™ II Reverse Transcriptase (Takara, Dalian, China).

Gene-specific primers were designed in NCBI and used for qRT-PCR. The analysis of PCR products was done according to procedure [37]. Data were analyzed using the ΔΔCT method. Primers used in this study are listed in Additional file 4: Table S2.

Histological analysis

Stems and petiole were cut into 0.5-cm segments and fixed in FAA stationary liquid at 4 °C for 3 days, drawn under a vacuum for 1 h, dehydrated in a graded ethanol series (2-h each), cleared in dimethylbenzene: ethanol (50: 50, v/v), and then cleared in 100% dimethylbenzene twice, for 2 h each time. Tissues were embedded in dimethylbenzene: wax (50: 50, v/v) overnight at 60 °C for 6 h, and then embedded in paraffin. The embedded tissues were sectioned to 4-μm thickness for staining using a Leica RM 2235 microtome (Leica) and adhered to microscope slides (Thermo Fisher) at 42 °C for 15 min. Afterward, some sections were stained with safranin and fast green for lignified cell wall and cell wall of phloem observation, which, respectively, take on red and green color under a light microscope. Finally, imaging was performed with an Olympus BX51 light microscope. The vascular bundle areas of three different tissue samples were quantified in a randomly selected area using Image J software (http://rsbweb.nih.gov/ij/).

Cellulose, hemicellulose, pectin and lignin content analyses

Cellulose was extracted according to the procedure used by Sun et al. [41] with some changes. Samples were analyzed for hemicellulose and pectin as previously reported [42]. Briefly, the stems were ground in liquid nitrogen, followed by one extraction with 70% ethanol and three extractions with chloroform: methanol (1:1, v/v). Aliquots (0.50 g and 0.10 g) of alcohol-insoluble residue were prepared for determination of cell wall components. First, the 0.50-g aliquot was homogenized in 90% DMSO, after which it was resuspended in 50 mM CDTA, 50 mM Na2CO3 and 24% KOH and centrifuged at 8000 r/min for 5 min; the residue was washed three times with 100% acetone and then freeze-dried; the solids were considered as cellulose. The 0.10-g aliquot was treated with 50 mL of 72% sulfuric acid at room temperature for 1 h. Sulfuric acid was then diluted to 4%, and the mixture was heated at 121 °C for 1 h. Released monosaccharides were analyzed using an ICS-3000 HPLC system (Thermo Fisher Scientific). Hemicellulose was calculated using the amount of Glc released by 4% sulfuric acid. Pectin concentration was calculated based on the amount of galacturonic acid (GalA). The lignin concentration was measured according to our previous method [43], in which extracted ground stem tissue (0.1 g) was treated with 3 ml of 72% H2SO4 and stirred every 10 min for 2 h. The sample was then diluted with 112 ml of deionized water and autoclaved at 121 °C for 1 h. The acid-soluble lignin component was determined at 205 nm by spectrophotometry.

Statistical analysis

SPSS Statistics 21 (SPSS, Inc., Chicago, IL, United States) was used to analyze all data in this study. Data was graphed with Sigma Plot 12.0 software. Data were analyzed using an independent t test, with a significance level of P ≤ 0.05. Values are presented as the mean ± standard error (SE) in at least biological triplicate for each measurement.

Availability of data and materials

The data sets supporting the conclusions of this article are included within the article and its additional files.

Abbreviations

Suc:

Sucrose

Glc:

Glucose

Fru:

Fructose

OE:

Overexpression

UDPG:

UDP-glucose

SUSY:

Sucrose synthase

F6P:

Fructose 6-phosphate

G6P:

Glucose 6-phsophate

CESA:

Cellulose synthase A

NINV:

Neutral invertases

HxK:

Hexokinase

FRK:

Fructokinase

CaMV 35S:

Cauliflower mosaic virus 35S

WT:

Wild type

GC/MS:

Gas chromatography mass spectrometry

CWINV:

Cell wall invertase

G1P:

Glucose 1-phsophate

ph:

Phloem

xy:

Xylem

C:

Cambium

References

  1. Kim SJ, Chandrasekar B, Rea AC, Danhof L, Zemelis-Durfee S, Thrower N, Shepard ZS, Pauly M, Brandizzi F, Keegstra K. The synthesis of xyloglucan, an abundant plant cell wall polysaccharide, requires CSLC function. Proc Natl Acad Sci. 2020;117(33):20316–24.

    Article  CAS  PubMed  Google Scholar 

  2. Silva-Sanzana C, Estevez JM, Blanco-Herrera F. Influence of cell wall polymers and their modifying enzymes during plant–aphid interactions. J Exp Bot. 2020;71(13):3854–64.

    Article  CAS  PubMed  Google Scholar 

  3. Hickey RJ, Pelling AE. Cellulose biomaterials for tissue engineering. Front Bioeng Biotechnol. 2019;7:45.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Polko JK, Kieber JJ. The regulation of cellulose biosynthesis in plants. Plant Cell. 2019;31(2):282–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. McFarlane HE, Döring A, Persson S. The cell biology of cellulose synthesis. Annu Rev Plant Biol. 2014;65:69–94.

    Article  CAS  PubMed  Google Scholar 

  6. Persson S, Paredez A, Carroll A, Palsdottir H, Doblin M, Poindexter P, Khitrov N, Auer M, Somerville CR. Genetic evidence for three unique components in primary cell-wall cellulose synthase complexes in Arabidopsis. Proc Natl Acad Sci. 2007;104(39):15566–71.

    Article  CAS  PubMed  Google Scholar 

  7. Endler A, Persson S. Cellulose synthases and synthesis in Arabidopsis. Mol Plant. 2011;4(2):199–211.

    Article  CAS  PubMed  Google Scholar 

  8. Batisse C, Buret M, Coulomb PJ. Biochemical differences in cell wall of cherry fruit between soft and crisp fruit. J Agric Food Chem. 1996;44(2):453–7.

    Article  CAS  Google Scholar 

  9. Ivakov A, Flis A, Apelt F, Fünfgeld M, Scherer U, Stitt M, Kragler F, Vissenberg K, Persson S, Suslov D. Cellulose synthesis and cell expansion are regulated by different mechanisms in growing Arabidopsis hypocotyls. Plant Cell. 2017;29(6):1305–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Lampugnani ER, Flores-Sandoval E, Tan QW, Mutwil M, Bowman JL, Persson S. Cellulose synthesis–central components and their evolutionary relationships. Trends Plant Sci. 2019;24(5):402–12.

    Article  CAS  PubMed  Google Scholar 

  11. Patrick JW. Phloem unloading: sieve element unloading and post-sieve element transport. Annu Rev Plant Biol. 1997;48(1):191–222.

    Article  CAS  Google Scholar 

  12. Ruan YL. Signaling role of sucrose metabolism in development. Mol Plant. 2012;5:763–5.

    Article  CAS  PubMed  Google Scholar 

  13. Li M, Feng F, Cheng L. Expression patterns of genes involved in sugar metabolism and accumulation during apple fruit development. PLoS ONE. 2012;7:e33055.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Amor Y, Haigler CH, Johnson S, Wainscott M, Delmer DP. A membrane-associated form of sucrose synthase and its potential role in synthesis of cellulose and callose in plants. Proc Natl Acad Sci. 1995;92(20):9353–7.

    Article  CAS  PubMed  Google Scholar 

  15. Ruan YL, Llewellyn DJ, Furbank RT. Suppression of sucrose synthase gene expression represses cotton fiber cell initiation, elongation, and seed development. Plant Cell. 2003;15(4):952–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Rende U, Wang W, Gandla ML, Jönsson LJ, Niittylä T. Cytosolic invertase contributes to the supply of substrate for cellulose biosynthesis in developing wood. New Phytol. 2017;214(2):796–807.

    Article  CAS  PubMed  Google Scholar 

  17. Yang J, Zhu L, Cui W, Zhang C, Li D, Ma B, Cheng L, Ruan Y, Ma F, Li M. Increased activity of MdFRK2, a high-affinity fructokinase, leads to upregulation of sorbitol metabolism and downregulation of sucrose metabolism in apple leaves. Horticulture Res. 2018;5(1):1–12.

    Article  CAS  Google Scholar 

  18. Verbančič J, Lunn JE, Stitt M, Persson S. Carbon supply and the regulation of cell wall synthesis. Mol Plant. 2018;11(1):75–94.

    Article  PubMed  CAS  Google Scholar 

  19. Riggs JW, Cavales PC, Chapiro SM, Callis J. Identification and biochemical characterization of the fructokinase gene family in Arabidopsis thaliana. BMC Plant Biol. 2017;17(1):1–18.

    Article  CAS  Google Scholar 

  20. Roach M, Gerber L, Sandquist D, Gorzsás A, Hedenström M, Kumar M, Steinhauser MC, Feil R, Daniel G, Stitt M, Sundberg B, Niittylä T. Fructokinase is required for carbon partitioning to cellulose in aspen wood. Plant J. 2012;70(6):967–77.

    Article  CAS  PubMed  Google Scholar 

  21. Desnoues E, Génard M, Quilot-Turion B, Baldazzi V. A kinetic model of sugar metabolism in peach fruit reveals a functional hypothesis of a markedly low fructose-to-glucose ratio phenotype. Plant J. 2018;94(4):685–98.

    Article  CAS  PubMed  Google Scholar 

  22. Damari-Weissler H, Rachamilevitch S, Aloni R, German MA, Cohen S, Zwieniecki MA, Michele HN, Granot D. LeFRK2 is required for phloem and xylem differentiation and the transport of both sugar and water. Planta. 2009;230(4):795–805.

    Article  CAS  PubMed  Google Scholar 

  23. Li M, Li P, Ma F, Dandekar AM, Cheng L. Sugar metabolism and accumulation in the fruit of transgenic apple trees with decreased sorbitol synthesis. Horticult Res. 2018;5(1):1–11.

    CAS  Google Scholar 

  24. Yang J, Zhan R, Jin Y, Song J, Li D, An G, Li M. Functional analysis of the promoter of the MdFRK2 gene encoding a high-affinity fructokinase in apple (Malus × domestica). Scientia Horticult. 2020;265:109088.

    Article  CAS  Google Scholar 

  25. Zhang M, Xu X, Zheng Y, Zhang Y, Deng X, Luo S, Wu Q, Xu J, Zhang S. Expression of a plastid-localized sugar transporter in the suspensor is critical to embryogenesis. Plant Physiol. 2021;185(3):1021–38.

    Article  CAS  PubMed  Google Scholar 

  26. Zhong Y, Mellerowicz EJ, Lloyd AD, Leinhos V, Riding RT, Little CH. Seasonal variation in the nuclear genome size of ray cells in the vascular cambium of Fraxinus americana. Physiol Plant. 1995;93(2):305–11.

    Article  CAS  Google Scholar 

  27. Delmer DP, Haigler CH. The regulation of metabolic flux to cellulose, a major sink for carbon in plants. Metab Eng. 2002;4(1):22–8.

    Article  CAS  PubMed  Google Scholar 

  28. Coleman HD, Yan J, Mansfield SD. Sucrose synthase affects carbon partitioning to increase cellulose production and altered cell wall ultrastructure. Proc Natl Acad Sci. 2009;106(31):13118–23.

    Article  CAS  PubMed  Google Scholar 

  29. Ahmed M, Akhtar S, Fanglu M, Hasan MM, Shahid AA, Yanang X, Sarwar MB, Rao AQ, Husnain T, Wang X. Sucrose Synthase (SuSy) gene expression: an indicator for cotton fiber initiation and early development. Russian J Plant Physiol. 2019;66(1):41–9.

    Article  CAS  Google Scholar 

  30. Zhang B, Gao Y, Zhang L, Zhou Y. The plant cell wall: biosynthesis, construction, and functions. J Integr Plant Biol. 2021;63(1):251–72.

    Article  PubMed  Google Scholar 

  31. Li M, Wang S, Liu Y, Zhang Y, Ren M, Liu L, Lu T, Wei H, Wei Z. Overexpression of PsnSuSy1, 2 genes enhance secondary cell wall thickening, vegetative growth, and mechanical strength in transgenic tobacco. Plant Mol Biol. 2019;100(3):215–30.

    Article  CAS  PubMed  Google Scholar 

  32. Kaur S, Zhang X, Mohan A, Dong H, Vikram P, Singh S, Zhang Z, Gill KS, Dhugga KS, Singh J. Genome-wide association study reveals novel genes associated with culm cellulose content in bread wheat (Triticum aestivum, L.). Front Plant Sci. 2017;8:1913.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Coleman HD, Ellis DD, Gilbert M, Mansfield SD. Up-regulation of sucrose synthase and UDP-glucose pyrophosphorylase impacts plant growth and metabolism. Plant Biotechnol J. 2006;4(1):87–101.

    Article  CAS  PubMed  Google Scholar 

  34. Griffiths CA, Paul MJ, Foyer CH. Metabolite transport and associated sugar signaling systems underpinning source/sink interactions. BBA-Bioenerg. 2016;1857(10):1715–25.

    Article  CAS  Google Scholar 

  35. Wu AM, Hörnblad E, Voxeur A, Gerber L, Rihouey C, Lerouge P, Marchant A. Analysis of the Arabidopsis IRX9/IRX9-L and IRX14/IRX14-L pairs of glycosyltransferase genes reveal critical contributions to biosynthesis of the hemicellulose glucuronoxylan. Plant Physiol. 2010;153(2):542–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Li Y, Xu T, Tschaplinski TJ, Engle NL, Yang Y, Graham DE, He ZL, Zhou J. Improvement of cellulose catabolism in Clostridium cellulolyticum by sporulation abolishment and carbon alleviation. Biotechnol Biofuels. 2014;7(1):1–13.

    Article  CAS  Google Scholar 

  37. Stein O, Secchi F, German MA, Damari-Weissler H, Aloni R, Holbrook NM, Zwieniecky MA, Granot D. The tomato cytosolic fructokinase FRK1 is important for phloem fiber development. Biol Plant. 2018;62(2):353–61.

    Article  CAS  Google Scholar 

  38. Wang Z, Wei X, Yang J, Li H, Ma B, Zhang K, Zhang Y, Cheng L, Li M. Heterologous expression of the apple hexose transporter MdHT22 altered sugar concentration with increasing cell wall invertase activity in tomato fruit. Plant Biotechnol J. 2020;18:540–52.

    Article  CAS  PubMed  Google Scholar 

  39. Dai H, Li W, Han G, Yang Y, Ma Y, Li H, Zhang Z. Development of a seedling clone with high regeneration capacity and susceptibility to Agrobacterium in apple. Sci Hortic. 2013;164:202–8.

    Article  CAS  Google Scholar 

  40. Zhu L, Li B, Wu L, Li H, Wang Z, Wei X, Ma B, Zhang Y, Ma F, Ruan Y, Li M. MdERDL6-mediated glucose efflux to the cytosol promotes sugar accumulation in the vacuole through up-regulating TSTs in apple and tomato. Proc Natl Acad Sci. 2021;118:1.

    Article  Google Scholar 

  41. Sun Q, Foston M, Meng X, Sawada D, Pingali SV, O’Neill HM, Li HJ, Wyman CE, Langan P, Ragauskas AJ, Kumar R. Effect of lignin content on changes occurring in poplar cellulose ultrastructure during dilute acid pretreatment. Biotechnol Biofuels. 2014;7(1):1–14.

    Article  CAS  Google Scholar 

  42. Chen S, Jia H, Zhao H, Liu D, Liu Y, Liu B, Somerville CR. Anisotropic cell expansion is affected through the bidirectional mobility of cellulose synthase complexes and phosphorylation at two critical residues on CESA3. Plant Physiol. 2016;171(1):242–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Maloney VJ, Park JY, Unda F, Mansfield SD. Sucrose phosphate synthase and sucrose phosphate phosphatase interact in planta and promote plant growth and biomass accumulation. J Exp Bot. 2015;66(14):4383–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

This work was supported by the Program for the National Natural Science Foundation of China (No. 31872043) and the Training Program Foundation for the Young Talents of Northwest A&F University. We are grateful to Dr. Jing Zhang and Miss JIng Zhao (Horticulture Science Research Center, Northwest A&F University, Yangling, China) for providing professional technical assistance with LC–MS/MS analysis.

Funding

This work was supported by the Program for the National Natural Science Foundation of China (No. 31872043) and the Training Program Foundation for the Young Talents of Northwest A&F University.

Author information

Authors and Affiliations

Authors

Contributions

SJ and LZ wrote the article, planned and performed the experiments. CZ and BM were involved in planning the experiments and writing the article. FM and ML participated in the discussion of all the results. ML also critically reviewed the text. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Mingjun Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

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.

FRK enzyme activity in the transgenic poplars (OE#1, OE#4 and OE#9) leaves overexpressing MdFRK2.

Additional file 2: Fig. S2.

SUSY activity in stems of transgenic poplars (OE#1, OE#4 and OE#9) overexpressing MdFRK2.

Additional file 3: Table S1.

Phenotypic characteristics of the transgenic poplars (OE#1, OE#4 and OE#9) overexpressing MdFRK2.

Additional file 4: Table S2.

Primers used in this study.

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Su, J., Zhang, C., Zhu, L. et al. MdFRK2-mediated sugar metabolism accelerates cellulose accumulation in apple and poplar. Biotechnol Biofuels 14, 137 (2021). https://doi.org/10.1186/s13068-021-01989-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-021-01989-9

Keywords