Monitoring transformation of two tropical lignocellulosics and their lignins after residence in Benin soils

Thermally assisted Hydrolysis and Methylation (THM), and 2D-heteronuclear single quantum coherence nuclear magnetic resonance (2D HSQC NMR) spectroscopy were used to monitor the transformation of ramial chipped wood (RCW) from Gmelina arborea and Sarcocephalus latifolius, together with their organosolv lignins, following soil incubation in Benin (West Africa). Mesh litterbags containing RCW were buried in soils (10 cm depth) and were retrieved after 0, 6, 12 and 18 months of field incubation. Chemical analysis showed that total carbohydrate content decreased, while total lignin content increased as RCW decomposition progressed. Ash and mineral content of RCW increased significantly after 18 months of decomposition in soil. Significant N-enrichment of the RCW was determined following 18 months incubation in soils, reaching 2.6 and 1.9 times the initial N-content for G. arborea and S. latifolius. Results of THM showed that the S + G sum, corresponding to lignins, increased with RCW residence time in the soils, in contrast to the response of compounds derived from carbohydrates, the sum of which decreased. Remarkably, lignin interunit linkages, most notably β-O-4′ aryl ethers, β-β′ resinol, β-5′ phenylcoumaran and p-PCA p-coumarate, survived after 18 months in the soil, despite their gradual decrease over the duration of the experiment.

ents of studied samples varied significantly as shown by polynomial contrasts (linear and quadratic) with the progress of RCW decomposition in soils for G. arborea and S. latifolius. Total extractive contents decreased by more than 60% (compared to the initial content) for the two species after 18 months of soil incubation. This result is consistent with those described in other published reports, which indicated that extractive contents progressively decreased during wood decomposition 34,35 . Extractives, given that they are "free molecules, " most notably phenolic products, are the first to be released during decomposition. In both species, total carbohydrate Remaining mass, elemental analysis and degradability indices of RCW of G. arborea and S. latifolius during soil incubation. Figure 2 summarizes the data on remaining mass ( Fig. 2A), elemental composition and degradability indices (Fig. 2B) and mineral content (Fig. 2C) of the study RCW. RCW mass losses www.nature.com/scientificreports/ decelerated over time for S. latifolius and G. arborea and had linear trend after 6, 12 and 18 months in the soils ( Fig. 2A). The total mass loss of RCW throughout the decomposition period was higher for G. arborea (80%) than for S. latifolius (70%). The percentage of maximum mass loss was recorded between months 12 and 18 for S. latifolius RCW and between months 0 and 6 for G. arborea RCW. This difference in mass loss between the two species could be associated with the high amount of labile C in the carbohydrate and extractive components of G. arborea RCWand high quantities of recalcitrant C that was associated with the lignin of S. latifolius RCW, as reported in our previous work 2 . The C content for both species showed a relative decrease in time, while the O contents generally exhibited a slight increase, which was probably related to oxidative microbial transformation. More interestingly, the N content of RCW increased substantially, reaching 2.6 and 1.9 times the initial content for G. arborea and S. latifolius, respectively. Progressive N-enrichment of RCW indicates that its application could improve the release and availability of N in soils, which is fundamental to maintaining soil fertility and important for plant nutrition. These results are in agreement with published reports, which also confirmed an increase of N content in decomposing biomass in soils 29,45,46 . This could be explained, at least in part, by an immobilization stage involving microorganisms, the enzymes that they secrete, and the protein cores of the former which contain nitrogen 45 . Pei and coworkers demonstrated that N released during decomposition can be immobilized by microorganisms, stabilized in the soil organo-mineral particles or absorbed to sustain crop productivity 38,46,47 . C/N ratios are often used as degradability indices to explain the decomposition patterns of lignocellulosic materials under study. The initial RCW had a C/N ratio between 82.76:1 and 94.92:1 and, theoretically, would result in strong nitrogen immobilization. The C/N ratio of RCW progressively decreased throughout the 18 months of the study, achieving a value of almost one-third (around 33.32 ± 0.2) and half (40.65 ± 0.23) of the initial ratios that have determined respectively for G. arborea and S. latifolius RCW. There is a positive correlation (0.94, P < 0.001) between C/N ratio and the remaining mass; the C/N ratio of RCW at 18 months is slightly higher than 25:1, indicating N-immobilization, but to a lesser degree than in the initial stages of decay. For both species, there was a negative correlation between C/N ratio and N content (− 0.95, P < 0.001), indicating that C/N ratio was a key factor affecting N enrichment as reported by Pei et al. 46 .
Mineral contents of RCW ( Fig. 2C) varied in various linear and quadratic trends with the RCW decomposition progress in soils for both species. Phosphorus and potassium contents significantly decreased, while calcium and manganese contents increased in RCW after 18 months of decomposition in soils for both G. arborea and S. latifolius RCW. The increase in manganese content could be related to the activity of white rot fungi since their enzymatic system comprises normally manganese-dependent lignin peroxidase 41,48 . Magnesium content decreased with RCW decay in G. arborea, while the opposite result was observed for S. latifolius. This relative variation in mineral contents reflects different dynamics of their release during RCW mineralization in soils, which is intended to improve nutrient status of the soils.
Analytical pyrolysis Py-GC/MS of the RCW after residence in soils. Thermally assisted Hydrolysis and Methylation was conducted to follow changes in the composition of RCW from G. arborea and S. latifolius as these materials decomposed in the soil. Pyrograms of RCW from G. arborea and S. latifolius are shown in Figure S1 (supporting data). The compounds that are reported in Table 1 revealed that the syringyl guaiacyl (S + G) sum corresponding to lignins increased, while the sum of the compounds corresponding to carbohydrate content decreased. Several papers revealed that pyrolysis technique is less effective in detecting carbohydrate units in respect to the lignin component 49,50 . This is probably the result of pyrolytic rearrangement of polysaccharides in complex matrices combined with the derivatizing agent (TMAH) which prevent chemical access and are believed to negatively interfere in the diagnostic release of carbohydrates derivatives and polysaccharides components. However, these trends of decreasing carbohydrate content are also indicated by CH/L ratio (carbohydrate/lignin ratio), confirming the results of global chemical analyses that are presented in Fig. 1. The main guaiacyl (G) moieties that were released by RCW samples from G. arborea have been identified as peaks 10R, 16R and 22R (Table 1), while the main syringyl (S) moieties of the same samples were identified as peaks 25R and 26R. Results showed that G. arborea samples contain GS-type lignin derivatives, dominated by S units, as confirmed by their lignin analyses presented below and reported in our previous study 2 . The S/G ratio of G. arborea samples increased slightly with RCW decomposition in soils, except for RCW18. This result could be explained by the fact that lignins from G. arborea are S-dominant. However, the moieties consisted of S units that are known to cleave more rapidly. Due to the easier cleavage of the ether bond by which the S-units are preferentially linked, the S unit content remained high. The very high initial S/G in G. arborea samples could perhaps explain the slight increase of the S/G ratio in the RCW samples even after 18 months incubation in soils. As for the RCW samples from S. latifolius, the main G moieties were identified (Table 1) in addition to those that were identified in G. arborea samples. The major compounds based on S units are the same as those identified in G. arborea samples, but with a lower relative abundance. The relative abundance of lignin derivatives units (S, G) ( Table 2) showed that RCW from S. latifolius has GS-type lignin that is dominated by G units. The trend is linear (p < 0.0001) and, there is an indication of decrease in the S/G ratio, indicating greater recalcitrance of G-lignin. The lower initial value of the S/G ratio in S. latifolius compared to G. arborea is consistent with the results of other analytical methods (higher Klason lignin + acid solubles, higher C content) as has been reported in previous studies 28,30,31,51-53 . Yield, purity, and elemental analysis of organosolv lignins extracted from RCW as degradation proceeded in soils. It has been demonstrated previously 2 that the organosolv process is well-suited to get access to high purity lignins from the RCW of G. arborea and S. latifolius, thereby increasing their suitability for studies of chemical structures. www.nature.com/scientificreports/ RCW was greater than 15% and did not vary substantially with the progress of RCW decomposition within the experimental timeframe. However, the yield of organosolv lignin from S. latifolius RCW samples decreased with RCW incubation time, as reported in Table 3. For both species, lignin recovery decreased with RCW decomposition and this decrease is much greater for lignin of S. latifolius than G. arborea. Thermally assisted hydrolysis and methylation (THM) determined that more condensed structures accumulated in RCW lignins with increasing duration of incubation (Table 4), which could explain the difficulty of extracting greater proportion of lignins by the organosolv process. Nevertheless, the selected organosolv process seems appropriate, given that it yielded, from both species, organosolv lignins with purities higher than 95% (sum of Klason plus acid soluble lignin of the isolated lignins) and, therefore, suitable for structural studies. The elemental composition of the isolated organosolv lignins is presented in Table 3, revealing significant variation in C, H, O, and N contents. The presence of N in the lignins under study could be due to the potential formation of protein-lignin complexes during the organosolv pulping process 54 . Yet, this could also be explained by nitrogen presence in the native lignins, as nitrogen content has been determined in organosolv lignins extracted following all litter bag retrievals (0, 6, 12, 18 months) ( Table 3). These findings indicate that lignin-protein chemical bonds are strong and, therefore, difficult to remove through pulping 2,54 . It should be noted that the N content of the lignins significantly increased as RCW decomposition progressed in soils. This was demonstrated for all RCW samples (Fig. 2), indicating the potential reactions between lignins and proteins from enzymes or microorganisms as the likely sites of nitrogen accumulation upon RCW degradation in soils.

Py-GC/MS of RCW G. arborea and S. latifolius lignins.
Organosolv lignins from G. arborea and S. latifolius were also analyzed by the Thermally assisted Hydrolysis and Methylation (THM), as was the case for the RCWs (Table 1). Figure 3 shows the chromatograms of pyrolytic products that were detected by THM for the  Table S1 (Supporting data). Based upon the figure and Table 4, the distribution of the pyrolytic products from the lignins showed linear and quadratic variation in their chemical structure as the RCW degradation progressed in soils. The main S-unit moieties that were present in the lignins were identified as peaks 15L, 16L, 18L, 22L, 25L, 26L, 30L, 37L (L refers to lignin). While the G-unit moieties were identified as 7L, 10L, 11L, 17L, 19L,

Lignin yield (%) Lignin recovery (%) TL (%) KL (%) ASL (%) Glucose (%) Ash (%) C (%) N (%) H (%) O (%)
RCW lignin (G. arborea)    Table 4. S-unit abundance increased, while those of G-units decreased slightly for lignin samples as RCW decomposed in soils. p-hydroxyphenyl (H)-unit abundance is very low (0.33-0.94%), indicating that lignins are G-S type as determined previously. This distribution of lignin derivatives is similar to that observed following THM of lignin from several hardwoods 28,31 . Carbohydrate pyrolysis products such as anhydrosugars (1,6-anhydro-β-d-glucopyranose), pyrans and furans 28,52 were not detected (Fig. 3), confirming the purity of the organosolv lignin. This reults also confirms the results regarding the high purity of the lignins under study (low carbohydrate content) ( Table 3). The validity of choosing the organosolv process was affirmed, given that it led to high purity lignin preparations 33,55 . It shoud be noted that the recovery of lignin by organosolv process steadily decreased for lignins, which had been isolated from the RCW after 6, 12 and 18 months (Table 3). Nevertheless, the isolated lignins were determined to be of high purity, a compromise that must be made regardless of the method of lignin extraction that is chosen. The S/G ratio of the studied lignins steadily increased for the lignins isolated from RCW after long incubations in the soils, which could reflect the fact that a more easily soluble part of partially decomposed lignin has been recovered preferentially by the organosolv method that was applied with the same parameters throughout this study. This result could also indicate progressive depolymerization of lignin during decomposition, since it is well known that with increasing decomposition of woody biomass, the S/G ratio is usually increasing 15,21 .

Polymer properties of lignins isolated from the RCW samples. Values of the number average (Mn),
weight-average (Mw), and polydiversity index (PDI) are presented in Table 5. The molecular weight distribution curves are shown in Figure S2 (Supporting data). Mn values varied significantly with time. No significant difference for Mw and PDI was apparent over time or between species, but lignin that was isolated from RCW after 18 months of incubation in soil appears much less stable than that of the original material given that lower polydispersity is indicative of good physicochemical stability of lignin 56 . Mw and PDI both vary depending upon the lignin isolation process, molecular weight distribution methods, and the type of plant material and its degradation state 2,57 . In considering the Mw data of all samples as a function of their incubation time in soil, an increase (although not significant) was noted in the average molecular weight (Mw) of all lignin samples, from 1703 to 2028 g/mol for G. arborea and from 1692 to 2158 g/mol for S. latifolius. This response clearly indicates that lignin polymerization reactions with enzymes of microorganisms are taking place, as observed in the RCW analysis.

2D-HSQC NMR analysis of the organosolv lignins isolated from the RCW of G. arborea and S. latifolius.
To obtain additional information regarding the progress of changes in the studied lignin structures after incubation in the soil, 2D-HSQC NMR analyses were performed based upon methods that have been described in previous studies 20, [57][58][59] . Figures 4 and 5 present the aliphatic oxygenated region (δ C /δ H 50-90/2.5-5.8) and the aromatic/unsaturated regions (δ C /δ H 90-155/6.0-8.0) of the 2D-HSQC NMR spectra of the organo-  (Table 6). It is remarkable that the β-O-4′ moiety (Fig. 6), which could be regarded as the hallmark of lignin identity, is preserved. HSQC spectra revealed the presence of other characteristic signals which correspond to C α -H α (71. In the aromatic unsaturated region of the HSQC spectra (Figs. 4, 5, bottom panel), the main correlation peaks that were found correspond to the aromatic rings of different lignin units (S, G), and to p-hydroxycinnamates (p-coumarates, pCA). A strong characteristic signal of S and G units was identified by their correlation peaks at δ C /δ H 104.2/6.61, 110.7/6.98, 115.5/6.67, 119.9/6.91 ppm corresponding to S 2,6 , G 2 , G 5 and G 6 , respectively. Also, the signal corresponding to Cα-oxidized S-units (S′ 2,6 ) at δ C /δ H 107.4/7.28 ppm was only observed in lignins from G. arborea samples. All study lignins were determined to contain a signal of p-coumarate (pCA), which was confirmed by important correlation (δ C /δ H ) at 115.7/6.69 ppm (pCA 3,5 ), while signals corresponding to p-coumarate in C 7 -H 7 (pCA 7 ) (137.3/7.87 ppm) and pCA2,6 (124.7/7.44 ppm) were only found in the spectral data for lignins from S. latifolius.  www.nature.com/scientificreports/ The relative content of the main lignin substructures, the molar content of S and G units, as well as the molar S/G ratio of the study lignins was evaluated from the contour integration volume in the 2D HSQC spectra. The results are presented in Table 6 20,57-59 . The semiquantitative results of the lignins under study demonstrate a predominance of β-O-4′ aryl ether linkages, which decreased from 63.34 to 37.18% for G. arborea and from 65.67 to 34.92% for S. latifolius as the duration of soil incubation increased. At the same time, the relative content of β-β′ linkages decreased from 9.83% to 2.51% in G. arborea and from 20.1% to 7.49% in S. latifolius, respectively. Furthermore, β-5′ substructure was only determined for S. latifolius samples and it decreased from 7.77 to 2.98% with increasing length of incubation.
The results indicate a progressive degradation of lignins with the increase in residence time of RCWs in soils, which seems to be more pronounced for lignin of G. arborea than that of S. latifolius (Table 6). Previous research has indicated that the cleavage of β-O-4′ linkages was the main mechanism of lignin depolymerization during wood decomposition 41,59 . This supports our result suggesting that RCW decomposition in soils is linked to the decline in β-O-4′ aryl ether linkages; nevertheless, they remain important moieties in the study lignins. Moreover, the S/G ratio of lignin is considered to be a major indicator of lignin depolymerization by white-rot fungi in soils 15,21,41 . The S/G ratios that were calculated increased from 1.68:1 to 2.07:1 and from 0.83:1 to 1.01:1 for RCW lignins that were isolated respectively from G. arborea and S. latifolius. Increased S/G ratios indicated that effective degradation of lignin was occurring during RCW decomposition, which agrees with thermochemolysis using tetramethylammonium hydroxide (TMAH) results. It is also important to note the incorporation of p-hydroxycinnamates in lignin of Sarcocephalus latifolius RCW, which agrees in part with results of the thermochemolysis using TMAH that are presented in Table 4.

Conclusion
The study presented here provides insights into the transformation of the RCWs of G. arborea and S. latifolius and the fate of their lignins following incubation in the soil at two different sites in Benin. As RCW decomposition advanced, total carbohydrate and extractive contents decreased, while lignin content of RCW increased. The significant increase of ash contents of the RCW after 18 months of decomposition (up to 4 times) could suggest a potential contribution of microbial activity to mineralization, leading to release of minerals favourable for soil fertility and crop nutrition. The important enrichment in nitrogen of the RCW following 18 months incubation in soils, reaching 2.6 and 1.9 times the initial N-content for G. arborea and S. latifolius, respectively, could also be linked to microbial activity. This N-enrichment of the RCW indicates that their application could improve the release and availability of N in soils, which is fundamental to maintaining soil fertility and important for plant nutrition. Analytical pyrolysis showed an increase in the S + G products sum, which clearly indicated to the relative lignin enrichment of the RCW from G. arborea and S. latifolius after residence in soils. Furthermore, the S/G ratio of RCW from G. arborea increased slightly, while the opposite trend was observed for S. latifolius after long-term RCW decomposition. Organosolv lignins that were isolated from RCWs are G-S type and their purity was determined to be greater than 95%. An increase in S/G ratio was determined for isolated lignins by Py-GC-MS analysis, which was more prominent for G. arborea than for the S. latifolius lignin. This response indicated a pattern of slower decomposition in S. latifolius samples in soil compared to G. arborea samples, that could be favorable to soil organic carbon stabilization. Remarkably, the β-O-4′ substructures were the main substructures of all studied lignins; which remained to be important throughout the process of RCW decay. This is a significant finding, as it could indicate to the survival of the β-O-4′ substructures in lignins, after longer residence times in soils. These moieties were accompanied by β-β′, and pCA-containing moieties in both lignin samples, while β-5′ moieties were only detected in S. latifolius lignin, thereby confirming the great recalcitrance of these structures. The combination of analytical pyrolysis with spectroscopic and chromatographic methods that were applied in this study provided important insights regarding the progress of lignocellulosic decomposition and the fate of its lignins following incubation in tropical soils.

Materials and methods
RCW sampling and preparation. Ramial chipped wood (RCW) of S. latifolius and G. arborea was collected in May 2018, respectively from a 10-year-old plantation that was located near the Kika Agricultural Technical High School (09°17′09″N, 2°45′06″E; Northeastern Benin) and a 5-year-old plantation at the Mèdji Agricultural High School in Sékou (06°41′64″N, 02°34′82″E; Atlantic coastal Benin). Permission to the collection of www.nature.com/scientificreports/ young branches from both S. latifolius and G. arborea complies with institutional, national, and international guidelines and legislation (Certificate issued by the environmental NGO APRETECTRA, Cotonou, Benin). The selection of these wood species was based on their availability and proximity to the experimental sites applying the RCWs. RCWs were applied in an experimental design that is part of a larger agronomic trial at two sites (Sekou and Kika) with contrasting climate conditions in Benin. On one hand, Sekou site belongs to the Guinean region which is humid area with a bimodal rainfall patterm (average of 1200-1500 mm/year) 61 62 ) and Alfisol (US soil classification) for which we have determined a pH of 6.6, C and N content of 10.72 g/kg and 0.97 g/kg respectively, bulk density of 1.24 g/cm 3 , and a silty-sandy texture at the first 30 cm depth. One should note that the site choice and consequently of their soils has been dictated by the larger framework of the project which involved the participation of two agricultural schools and their students in the project's follow-up. For each species, 100 g of RCW (with 40.2% initial humidity) were placed in 20 × 20 cm nylon mesh litter bags (with 1 mm-mesh size, therefore porous and permitting contact with soil) 14,63 , which were incorporated into the soils at a depth between 1 and 10 cm, then the top of litter bags was levelled with soil surface. Sixty-four bags of each RCW species (with initial humidity determined at 40.2%) were initially buried in the soils in the experimental design. Buried RCWs are therefore in contact with the soil, which favored the decomposition process. We retrieved the nylon litter bags from the soils at six-month intervals for the total 18 month duration of the experiment. RCW samples were collected at four time periods: RCW0 (initial RCW sample not incorporated in soils); RCW6, RCW12, and RCW18 are RCW samples removed respectively after 6, 12 and, 18 months of residence in soils. At each sampling, the bags with RCW were carefully washed and rinsed in distilled water to separate them from the mineral soil particles as thoroughly as possible before oven-drying at 40 °C. The oven-dried RCW samples were weighed, then ground and sieved through a shaking screen. The sawdust fraction between 40 and 60 mesh was collected for chemical analyses.
RCW chemical composition determination. Total extractive content of RCW samples was determined by Soxhlet extraction for 6 h with anhydrous toluene-ethanol (0.427:1.000, v/v), followed by hot water extraction. Klason and acid-soluble lignin contents of G. arborea and S. latifolius RCW were determined using standard methods 64 . Carbohydrates were determined using the protocol of National Renewable Energy Laboratory (NREL, Golden, CO, USA) 65 . The analysis was conducted by high performance liquid chromatography using a refractive index detector (HPLC-RID), and performed on an Agilent Technologies 1200 series system equipped with a Rezex RHM-Monosaccharide H + (8%) column (300 × 7.8 mm). Ash contents were determined according to the ASTM standard method 66 . Mineral contents (Ca, Mg, K, P, Mn) of study samples were determined using an optical emission spectrometer with excitation by inductively coupled plasma (ICP) on an Optimum Machine DV 4300 (Perkin Elmer, Norwalk, CT, USA) with a Scott-type of nebulizer. Elemental analysis of RCW samples were performed using a Perkin Elmer 2400 Series II (CHNO analyzer). The contents (%) of C, H and N in RCW sawdust are based on ash-free dry mass, while O content was computed by difference, i.e., O = 100 -(C + H + N). Each result corresponds to the mean of three replicates.
Organosolv lignin extraction. The organosolv lignins were isolated from the RCW samples after 0, 6, 12 and 18 month soil incubations after pre-extraction of the biomass with ethanol-water mixture, followed by the pulping according to the procedure that was developed in our laboratory and demonstrated to provide high purity lignin (containing low residual sugars and minerals) 2 . Briefly, 100 g of extractive-free RCW biomass were treated with 1 L of ethanol-water mixture (1:1, v/v), and catalyzed with 6 mmol ferric chloride (FeCl 3 ·6H 2 O) in a Parr reactor series 4842 (2 L) at 180 °C for 90 min. The residual liquor was obtained by vacuum filtration to separate it from the solid organosolv pulp, evaporated to remove ethanol, and then precipitated by acidification (2 M HCl, pH = 1.5). Finally, the organosolv lignin was recovered after filtration and oven-dried at 40 °C overnight. Lignin yield and recovery were calculated using Eqs. (1) and (2), respectively.
Py-GC/MS analyses. RCW  www.nature.com/scientificreports/ methylation and thus overall volatility of pyrolysis products but also improved the extent of the pyrolysis process which yielded higher quantities of pyrolysis products available for detection. A 0.4 mg portion of the samples was impregnated with a 10 µL of TMAH in methanol at 25% prior to loading them into ferromagnetic tubes for drying for 30 s at 100 °C. The sample (RCW, lignin) was pyrolyzed according to the following program: both the GC injector and the temperature of the pyrolyzer transfer line were set at 250 °C and maintained for 10 s. The temperature of the transfer line was then increased to 550 °C at a rate of 20 °C/s and held for 10 s at the temperature. The carrier gas (helium) was operated at the flow rate of 1.0 mL/min. The mass spectrometer was set in electron impact mode (EI, 70 eV, m/z = 35-400) at 1 s per scan. Each chromatogram peak was identified with the mass spectral library of the instrument (NIST) and completed using data that were available from published work 28,30,31,[51][52][53] . Three replicates were examined for each sample. Relative abundance of each compound was calculated by dividing its GC peak area by the summed areas of all selected pyrolysis products.
2D NMR analyses of lignins. About 80 mg of organosolv lignin sample was dissolved in 0.7 mL of dimethylsulfoxide-d 6 (DMSO-d 6 ); its nuclear magnetic resonance spectra were recorded on a Varian NMR spectrometer at 500 MHz at 25 °C in DMSO-d 6 according to published methods 2,19,37,60 . The program used was the Bruker standard pulse program "hsqcetgpsi". The following parameters were using per sample for a total data acquisition of 5 h: 2048 data points acquired from 10 to 0 ppm in F2 followed by a 1 s recycle delay, 160 to 0 ppm in F1 (13C) with 256 increments of 64 scans. Semi-quantitative analysis was performed using volume integration of contours in the 2D HSQC spectra to calculate relative numbers of inter-linkages and the syrinyl/guaiacyl (S/G) ratio of the study lignins.

Molar mass distribution.
The molar mass distribution of the study organosolv lignins was determined by gel permeation chromatography (GPC) on an Agilent 1200 series system under the following conditions: column, PL gel 5 μm Mixed-D 300 × 7.5 nm; mobile phase, tetrahydrofuran (THF); polystyrene, standard (580-28,770 Da, Agilent). Prior to injection, 40 mg of sample were dissolved in 4 mL of tetrahydrofuran and filtered through 0.45 μm porosity filter. The chromatographic analysis was performed in triplicate per sample. From the molar mass distribution, the number average molecular weight (Mn), the weight average molecular weight (Mw) and the polydispersity index (PDI = Mw/Mn) were determined.
Statistical data analysis. Differences in chemical properties between RCW samples of study species that were recovered from soils, together with the polymer properties of isolated organosolv lignins, were performed using one-way analysis of variance. Orthogonal polynomial contrasts were used to determine regression polynomial significance (linear, quadratic) at 5%. All analyses were performed using the Statistical Analysis System (SAS) version 9.4 67 .