Physical–chemical properties of cell wall interface significantly correlated to the complex recalcitrance of corn straw

Background Tissue heterogeneity significantly influences the overall saccharification efficiency of plant biomass. However, the mechanisms of specific organ or tissue recalcitrance to enzymatic deconstruction are generally complicated and unclear. A multidimensional analysis of the anatomical fraction from 12 corn cultivars was conducted to understand the essence of recalcitrance. Results The results showed that leaf, leaf sheath, stem pith and stem rind of corn straw exhibited remarkable heterogeneity in chemical composition, physical structure and cell type, which resulted in the different saccharification ratio of cellulose. The high saccharification ratio ranging from 21.47 to 38.96% was in stem pith, whereas the low saccharification ratio ranging from 17.1 to 27.43% was in leaf sheath. High values of lignin, hemicelluloses, degree of polymerization and crystallinity index were critical for the increased recalcitrance, while high value of neutral detergent soluble and pore size generated weak recalcitrance. Interestingly, pore traits of cell wall, especial for microcosmic interface structure, seemed to be a crucial factor that correlated to cellulase adsorption and further affected saccharification. Conclusions Highly heterogeneity in cell wall traits influenced the overall saccharification efficiency of biomass. Furthermore, the holistic outlook of cell wall interface was indispensable to understand the recalcitrance and promote the biomass conversion. Graphic abstract Supplementary Information The online version contains supplementary material available at 10.1186/s13068-021-02047-0.


Background
Governments and energy corporations around the world are actively boosting the efficient utilization of abundant, renewable and bio-degradable feedstocks, due to the exhaustion of fossil resources and the severe environmental pollution [1]. The output of lignocellulosic biomass can reach 181.5 billion tonnes annually [2]. Lignocellulosic biomass refers to the biological feedstocks, mainly including agricultural products, forestry wastes, and industrial residues. These feedstocks can be potentially converted to bioenergy, biochemicals and other bioproducts based on the biorefinery [3]. An integrated biorefinery is a sustainable facility involving in various physical, thermochemical and biochemical approaches, etc. [4]. It is largely acknowledged that saccharification of cellulose with economic competitiveness is deemed as one of the core requirements in biomass-based biorefinery.
Enzyme-catalyzed saccharification of cellulose is the most optimal approach, yet, the efficiency in this process is commonly hampered by various factors. Series of physical and chemical factors that hinder the enzymatic saccharification of cellulose are deemed as "recalcitrance", which is referred as (but are not limited to) composition (consisting of cellulose, lignin, hemicelluloses and pectin), polymer properties (such as lignin monomer components, cellulose crystalline index and degree of polymerization), and interactive network among these polymers [5]. Assessing and destroying the recalcitrance are prerequisite to circumvent the technical and economic obstacle in recent circular bioeconomy.
Plant biomass has evolved various functional and morphological fractions involving in complicated constituents and structures, which is named as ''heterogeneity". In terms of hemicelluloses and lignin, the amount, distribution and their decoration pattern add complexity to the origins and tissues of plant biomass [6]. Crowe also indicates lignin content commonly exhibits a negative correlation with saccharification ratios, yet, is not the sole contributor to recalcitrance of switchgrass [7]. Thus, a consensus can be reached that the high heterogeneity of plant biomass leads to the significant complexity of recalcitrance, which further influences the overall conversion efficiency of plant biomass. A further complication is that the complexity commonly causes empirically implemented enzymatic saccharification process and drastic pretreatment for recalcitrance more than necessary [8]. Hence, this immense heterogeneity is considered as the chief reason why the precise mechanisms of recalcitrance are still ambiguous and also closely influence the optimal choice of pretreatment and enzymatic hydrolysis strategy. Additionally, it is seemingly difficult to systematically compare these qualities to evaluate the intrinsic recalcitrance among different feedstocks due to the heterogeneity, especially for different studies.
Dried plant biomass is primarily composed of cell wall. The essence of enzymatic saccharification is the hydrolysis of cellulose in cell wall [9]. The enzymatic saccharification correlates to the properties of cell wall including the relative abundances of composition, fine proportion, their interconnecting among different matrices and multiple structures [10]. The basic types of cells among various plants are overall similar. Varied behaviors of different cell types in enzymatic hydrolysis reflect the tissue-specific recalcitrance that influenced by the nature of cell wall [11,12]. Thus, comprehending the resistant factors with viewpoint of cell wall might provide a deep insight into the recalcitrance of plant biomass. Allowing the enzyme access to the cell wall is the primary step in enzymatic hydrolysis of cellulose, which significantly influences the latter saccharification process [13]. However, the susceptibility of different cell wall, especially the accessibility of bio-interface to enzyme is still unclear.
Corn is the main agricultural crop on the earth. Corn straw (CS) is the main residue of corn. Based on the desirable agronomic and biochemical properties in terms of high biomass yield and cellulose content, and C4 carbon fixation, CS has attracted huge attention among researchers for biomass biorefinery. From visual angle, it is easily found that leaf, leaf sheath, stem rind and stem pith have different morphologies. Thus, CS is especially suitable model to clarify tissue-specific recalcitrance for its greater diversity of cells than other plant biomass. Although enzymatic saccharification of cellulose from CS has been studied by many researchers, little investigation concerning the relevance of enzymatic saccharification to heterogeneity in composition and architecture of specific organs or tissues has been reported. In an effort to clarify the fundamental knowledge about origins of recalcitrance, enzymatic saccharification and physicalchemical properties of cell wall in corn cultivars with different macroscopic phenotypes were determined. And the relationship between multidimensional properties of cell wall and the enzymatic digestibility of cellulose were discussed. This paper focused on understanding the recalcitrance of cell wall interface, which will also help in selecting of corn varieties with improved biorefining capabilities and setting downstream pretreatment.

Saccharification performance
The 144 saccharification ratio values of the four anatomical fractions from 12 corn cultivars were used to generate a heat map contained in a matrix with legend color bar. Dendrograms were selected to describe the similarity between clusters. Hierarchical clustering analysis confirmed the initial idea that some cultivars necessarily had a comparable performance in saccharification (Fig. 1). Additionally, there were no obvious correlations between the saccharification ratio and the macroscopic phenotypes of corn straw such as height, dry weight, etc. The general varying patterns for saccharification ratio among the same fraction within the different cultivar was similar, suggesting that biomass saccharification obviously correlated to the specific organs or tissues (Fig. 2). The highest saccharification ratio was obtained from stem pith, and followed by leaf, stem rind and leaf sheath in sequence, despite the difference between the stem pith and leaf was not significant. Specially, the saccharification ratio of stem pith ranged from 21.47 to 38.96%, whereas it ranged from 17.1 to 27.43% among leaf sheath.
Enzymatic digestibility expressed as saccharification ratio can be used to reflect the recalcitrance. Higher enzymatic digestibility indicates weaker recalcitrance [14]. The outermost leaf sheath was the most recalcitrant fraction, and followed by the stem rind. The leaf was the least recalcitrant fraction except for the stem pith. Hence, recalcitrance exhibited obvious heterogeneity among organs or tissues. Li also reports that different fractions of sorghum with structural or chemical variations generate special recalcitrance [15]. This lack of agreement about recalcitrance of leaf, leaf sheath, stem rind and stem pith indicated the complication of cell wall assembly, which might further translate dictate recalcitrance of cell wall to sugar production. Saccharification of cellulose is an enzymatic hydrolysis of cell wall in plant, exhibiting as glucose release from cellulose. Thus, the relationships between the enzymatic digestibility and physiochemical characteristics of cell wall were analyzed.

Associations between cell wall properties and their impact on enzymatic digestibility
Representative DK and YH were chosen to the further study, because they were expected to exhibit extreme and contrasting phenotypes. DK had the highest biomass and height among the 12 cultivars, yet, YH represented the opposite phenotypes (the lowest biomass and height). Physical-chemical properties were determined for each leaf, leaf sheath, stem rind and stem pith, and used to evaluate the relationship between cell wall and recalcitrance. The values of enzymatic digestibility were different for the stem pith and stem rind across DK and YH, yet similar for the leaf sheath (Additional file 3: Table S2). Overall, enzymatic digestibility of leaf, stem rind and stem pith were significantly higher in DK than in YH (P < 0.01). Detailed information concerning the different physical-chemical properties and significance analysis is provided in Tables 1 and 2. The general varying patterns for physical-chemical properties within the same anatomical fraction were similar. The crystallinity index, degree of polymerization, and ash values were higher in stem rind and leaf sheath than in stem pith and leaf. Leaf sheath showed the least content of NDS, whereas leaf contained the highest content of NDS, due to the large amounts of soluble sugars and protein. Similar trend is also observed in sorghum stem [16]. Higher crystallite dimension and DB/DO ratios were in stem pith or leaf. Moreover, higher cellulose, hemicelluloses and lignin values were in stem rind or leaf sheath. The collective data confirmed that characteristics of cell wall from different organs or tissues were related to the varied recalcitrance.
A correlation matrix was performed to evaluate how the recalcitrance was associated with cell wall properties (Fig. 3). The results showed that the contents of cellulose and NDS positively impacted the enzymatic digestibility. The positive effect of cellulose is not surprising for it probably reflected a greater availability of substrate for enzymatic hydrolysis. Although NDS might interact with hemicellulose and lignin to influence the enzymatic digestibility of cellulose, higher content of NDS commonly indicated the lower relative content of other relative recalcitrance [17]. The influences of ash and crystallite dimension among the subset of samples on enzymatic digestibility were not statistically significant. The negative impact of lignin, degree of polymerization and crystallinity index values on enzymatic hydrolysis were in accordance with the recalcitrance model reported in the previous literature [18,19]. The negative impact of hemicelluloses was not statistically significant. Because, the physical barriers of hemicelluloses to block cellulose is highly correlated to the ultrastructure of cellulose especially for the accessibility of cellulose surface [20]. Additionally, the leaf sheath had intermediate lignin content and crystallinity index compared to the other tissues yet the lowest DO/DB ratio and digestibility, suggesting the DO/DB ratio had a positive effect on enzymatic digestibility. The DO/DB ratio determined by Simon's staining can be used to describe the relative interior and exterior accessibility of pore [21]. Enzyme accesses to the cellulose core needs to be in intimate and prolonged contact [22]. Pore properties involving in volume, specific surface, tortuosity, size and fractal dimension might all affect accessibility of cellulose. Thus, pore characteristics and size distribution of various samples were determined.

Desirable pore traits of cell wall for enzymatic digestibility
Generally, the leaf had the highest average pore diameter, whereas the stem pith had the highest total pore area and porosity (Table 3). These results were significantly positively correlated with enzymatic digestibility, indicating that higher pore size and more pore number could facilitate enzymatic hydrolysis. Differential intrusion volume vs. diameter was also plotted (Fig. 4). The pore size distribution trends of DK and YH were similar. The different pore especially for the class corresponding to the pore size ranging from 5 to 15 nm obviously correlated to the enzymatic digestibility. Plant biomass can be deemed as special porous material, consisting of highly ordered porous structure at various levels, such as cell lumen, inter-cellular spaces, pit in cell wall, and space among chemical polymers. Arantes suggests that pore size distribution is mainly responsible for the efficiency of enzymatic hydrolysis [23]. Sun suggests that the specific structure of pore not the pore volume impacts the efficiency of enzymatic hydrolysis [24]. A large fraction of the pore volume might be ascribed to pores with small size. Pores with size bigger than 5.1 nm are sufficient for accessibility of cellulase, but small pores might resist mass transfer, thus hindering enzymatic hydrolysis [25]. In sum, catalytic cleavage of the cellulose was not the actual rate-limiting step, but rather the limited accessibility of the cellulose chains to the enzyme within the substrate. Pore size instead of volume could accurately represent the relationship between the porous structure and accessibility of the cellulose.

Morphology variation and cellulase adsorption
The labeled cellulase with the specific function of recognizing and adsorbing on cellulose was selected to observe the cellulase adsorption of different cell wall.
The imaging results displayed that all the fractions of CS gave the complete structure after sectioning, except for the leaf sheath (Fig. 5). The transverse view exhibited the typical tissue structure of monocotyledons consisting of parenchyma cells and scattered vascular bundles. Parenchyma cells accounted for the half of stem rind and stem pith roughly. Additionally, leaf constituted of the mesophyll, vascular bundles and epidermal cells. Chemical mapping of chromophores on both the surface and inner fractions of cell wall was revealed by CLSM. The red fluorescence emission of DyLight-labeled enzyme was not affected by blue lignin autofluorescence. The variations of red fluorescence from the same fraction of DK or YH were similar, yet significant different within different fractions of the same cultivar. Precisely, the labeled cellulase interacted with the surface at or within cell wall of stem pith strongly, yet very weakly in the leaf sheath. The labeled cellulase was bound to the mesophyll and vascular bundles distinctly but not to the adjacent epidermal cells in leaf. Epidermal cells with special structure encrusted with waxes and cutin mainly envelop leaves, stem rind and leaf sheath. The cuticle is overlaid by wax deposits that play a role in resistance to microbes and insects, and water retention in the nature [26]. Thus, dermal cells provided rigidity, which was significantly resistant to enzyme adsorption. From the view of plant anatomy, vascular bundles surrounding sclerenchyma tissue often have thickened and lignified secondary walls, whereas parenchyma cells only have primary walls and commonly do not have lignified secondary walls. Interestingly, more cellulase was adsorbed on the highly lignified vascular bundles than the less lignified parenchyma in the stem pith. However, Ding reports that general digestibility is significantly negatively correlated with lignin content of cell wall [27]. This might be caused by the non-productive binding of lignin in vascular bundles. Donaldson shows a moderate correlation between the enzyme adsorption and the cell wall histochemistry [28]. Our previous results also demonstrate that the lignin distribution, and interactions between lignin and cellulose could reduce the digestibility of cellulose [29]. These results revealed that the complicated physical-chemical properties of cell wall interface and its inter-relationship significantly affected the rankings of cellulase adsorption and also correlated to recalcitrance. The recalcitrance of sugar cane becomes weak gradually from the outer to the inner fractions [30]. Recalcitrance of cell wall in Miscanthus correlates to the tissue-specific distribution of lignin and development stage [31]. The heterogeneity of plant biomass highly depended on their origin, organ type, tissue trait, cell wall architecture, monomer composition, chemical bonds, its stage of development, etc., which further caused heterogeneity of the recalcitrance. Multiple architectures of cell wall, especially for the interface made the glycosyl hydrolase more difficult attack the plant cell wall. To understand the recalcitrance at cell wall interface might be the vital step for the selections of pretreatment conditions and enzymatic hydrolysis strategies.

Conclusions
The compiled data revealed that the heterogeneity in structure and physicochemical properties of cell wall were distinct among organs or tissues of CS, and significantly correlated to the recalcitrance. Critical factors of cell wall had obvious impacts on cellulase adsorption that furthermore determined enzymatic digestion. Among a set of variables, pore size was shown as the predominant factor that affected the access of cellulase to polysaccharides within the cell wall. Thus, a holistic view of microscopic interface of cell wall was proposed, considering that concrete variable have interactive impacts on cellulase adsorption depending on overall interface structure assembly. Furthermore, a combination of many variables might contribute to saccharification efficiency, indicating that biorefineries should be the optimal choice to better exploit heterogeneity and optimize biomass flows.

Materials
Twelve corn (Zea mays L.) cultivars with different macroscopic phenotypes including dry weight, height and diameter were grown in experimental field of Helin, Inner Mongolia, China (Additional file 2: Table S1). The total aboveground biomass was harvested for this experiment when it was mature. The fresh corn straw without internode was selectively separated into four anatomical fraction including leaf, leaf sheath, stem rind and stem pith by hand, based on their macroscopic appearance and function (Additional file 1: Figure S1). After that the sample was washed by tap water and air dried. The tested samples were random mixed for sampling uniformity and further divided into two fractions for different analysis requirement. One fraction was ground to powder with size ranging from 0.84 to 0.42 mm through high-speed rotary cutting mill (BJ-300A, Baijie, China). The other fraction was sectioned and used for microscopy. Based on the results of saccharification ratio, the DK and YH were selected for the further physical-chemical analysis of cell wall. Cellulase (Trichoderma viride G) was bought from Shanghai Yuanye Co., Ltd. The activity of cellulase presenting a filter paper unit of 110.2 U/g was obtained following the classic procedure described by Ghose [32].

Assessment of enzymatic digestibility
Enzymatic digestibility was determined via cellulasemediated hydrolysis and quantitation of reducing sugars generated from cellulose and hemicelluloses. The enzymatic hydrolysis was performed in 50 mM sodium citrate (pH 4.8) with solid consistency of 5% (w/w), enzyme loading of 20 FPU/g and at 50 ○ C for 48 h. The reducing sugars were assayed according to the 3,5-dinitrosalicylic acid method described by Miller [33,34]: where m g represented amount (g) of the generated reducing sugars during enzymatic hydrolysis, and m H-Cel and m Cel represented weight of hemicelluloses and cellulose in the tested sample (g), respectively.

Chemical composition analysis
Chemical characterization analysis was performed according to a protocol recommended by Van Soest [35]. Various dry samples were treated by neutral buffered detergent solution, hydrochloric acid (2 M) and sulfuric acid (7.34 M) to determine neutral detergent fiber (NDF), acid detergent fiber (ADF), and acid detergent lignin (ADL) in sequence. The ash content was assayed by heating the ADL at 550 °C for 180 min. The contents of neutral detergent soluble (NDS), hemicelluloses, cellulose, and lignin were calculated by subtracting the corresponding values from the initial dry weight, NDF, ADF, ADL and ash fractions in sequence.

Crystallinity analysis
Powder X-ray diffractometer (X'Pert PRO, PANalytical, Netherlands) was used to measure the Segal crystallinity index (CrI) and crystallite dimension (CD). Scans were performed in triplicate at 1°/min from 5° to 40° 2θ with a step size of 0.02°. The CrI was calculated from the XRD patterns based on the empirical peak-height method (Eq. 2) [36]. The dimensions of the crystallites in various samples were determined through Scherrer method (Eq. 3): I total was the diffraction intensity of major peak at 2θ between 22° and 23° for cellulose Iβ. I amor was the diffraction intensity of the minima between the major peak and secondary peaks.
where K was a constant value of 0.94; λ was the X-ray wavelength of 0.1542 nm; β was peak width of diffraction band at the half maximum height; τ represented the crystallite dimension; and θ was the Bragg angle fitted from the major peak by Inc Jade 6.5 software.

Degree of polymerization determination
Dry powders were extracted by 4 M KOH solution containing sodium borohydride (1.0 mg/mL) at 25 °C for 1 h, with solid-to-liquid ratio of 1:10. After being extracted for repeated twice, the alkaline slurry was further washed by distilled water and extracted with 10 mL 8% NaClO 2 solution (containing 1.5% glacial acetic acid) at 25 °C for 1 h to obtain crude cellulose. The degree of polymerization of the crude cellulose was determined by an Ubbelohde viscosity meter with the modified viscosity method at ambient temperature (around 25 °C) [37].

Simons' stain
The overall porosity was measured by Simons' stain [38]. The DB and DO dyes using the same amount of 10 mg/ml were mixed with the ratio of 1:1. Precisely 0.1 g of dry powders and increasing volumes of the mixed dye solutions (0.25, 0.5, 0.75, 1.00, 1.50 ml) were added into the centrifuge tubes. Phosphate buffered saline solution (pH 6, 0.2 M PO 4 , 1.40 M NaCl) was added into each tube to make the final volume up to 10 mL. After incubating at 70 °C for 6 h with gentle shaking, the tubes were centrifuged at 6200×g for 10 min. The absorbance of the supernatant solution in the tube was measured at the maximum wave length of (2) CrI = (I total − I amor )/I total .

Determination of porous parameters of corn straw
The porous parameters of various samples were evaluated by mercury porosimeter (AutoPore IV 9500, Micromeritics Instrument Corporation, USA), with a contact angle of 130° and pressure ranging 1 × 10 -1 to 6.1 × 10 4 psia. Intrusion pressure was directly converted to the relevant pore size through the Washburn equation [39].

Sectioning and cellulase labeling
To visualize the enzyme adsorption, pieces (around 5 mm) were cut from the middle portion of each fresh sample and embedded by tissue freezing medium (Leica, 4 fl) under -20 °C. The transverse slices (8 μm in thickness) of different samples were prepared using Leica cryostat (Leica, RM2015) for confocal laser scanning microscopy (CLSM, LSM-710, ZEISS, Germany). The cellulase labeling procedure was carried out as the manufacturer's instructions adopted by Donaldson with a slight modification [28]. Concretely, the freeze-dried cellulase powder was dissolved in borate-phosphate buffered saline (0.05 M, pH 7.5) and mixed with DyLight 633 ™ in the tube covered with aluminum foil at ambient temperature (around 25 °C) for 60 min, with gentle shaking [40]. And then the mixture was passed through the supplied resin and centrifuged multiple times (6200×g, 15 min) to separate the unlabeled fluorophore. The purified labeled cellulase was kept at 4 °C avoiding light. Cellulase (the given average molecular weight of 65 kDa), with fluorophore-to-protein molar ratio of 7.48:1, was used to label the section of samples for 60 min. The distribution of the label cellulase on cell wall was described by CLSM, using the identical instrument settings: 10 × /1.40 NA Plan-Apochromatic objective lens, pinhole size of 1 AU, 405 nm and 633 nm sequential excitation, and 410-480 nm (blue lignin autofluorescence), and 650-750 (red DyLight 633 ™ labeled cellulase fluorescence) emission. Lignin fluorescence and labeled enzyme were visualized by sequential excitation, without any significant bleedthrough of signals. All images were shown under the uniform condition with optimal intensity.

Statistical analysis
Heat map and dendrogram were plotted using the original "heatmap" function shipped with R installation, based on the hierarchical clustering results of saccharification ratio. Biological triplicates were selected for each corn cultivar. Chemical and physical analysis was carried out in technical triplicates to ensure reproducibility. The data were expressed as the mean of triplicates ± standard