3.1 Single factor analysis
To study the effect of the extraction temperature of 70, 75, 80, 85 and 90°C on the extraction yield of Sap when the extraction power (120 W), time (15 min) and liquid to material ratio (30:1). The result shows that the yield of Sap increased when temperature increased from 70 to 80 ℃, but then decreased (Fig. S2A). This may be because that the solubility of solute in solvent increases with the increase of temperature [35]. In addition, the viscosity of solvent decreased, which improved the diffusivity of solvent in the tissue [36]. With further increase in extraction temperature, the cavitation effect is gradually weakened, which led to the yield decreased. Therefore, the extraction temperature was selected 80 ℃.
Extraction was carried out at different extraction time (5, 10, 15, 20 and 25 min) when the extraction temperature, power, and material to liquid ratio was 80℃, 120 W, 1:30, respectively. The yield of Sap increased when extraction time increased from 5 to 15 min and then the yield decreased (Fig. S2B). These results can be explained that the increasing ultrasonic time led to solvent permeated into the dried raw materials, which promote the release and diffusion of Sap [37]. However, the ultrasound treatment for the prolonged may result in the destruction of Sap structure and decrease the extraction rate [38]. Therefore, the extraction time was selected 15 min as the central point of the RSM experiments.
To study the effect of the different liquid to material ratio (10, 20, 30, 40 and 50) on the extraction yield of Sap when temperature, time and power was set at 80 ℃, 15 min, 120 W, respectively. The yield of Sap increased when LMR increased from 10 to 40 mL/g and then decreased (Fig. S2C). The extraction rate of Sap increase, on the one hand, can be attributed to the increase of cavitation effect with the decrease of solvent concentration, resulting in the diffusion of polysaccharides show more quickly [39]. On the other hand, the concentration difference between the intracellular and extracellular of plant cells can also promote the diffusion of polysaccharides [2]. However, the excess solvents may cause the decrease in distribution of energy of the ultrasonic wave and has a negative effect on the yield of polysaccharides [2, 38]. Thus, the LMR was set as 40 mL/g.
The effects of extraction power on Sap yield were investigated when the experiments conditions were fixed at 80 ℃, 15 min, 40:1 and the extraction power was set at 60, 90, 120, 150 and 180 W. The result shows that the yield of Sap increased with the increase in ultrasonic power from 60 to 150 W and then decreased (Fig. S2D). This may be due to the bubble size increases with the increase in the ultrasonic power, resulting in the implosion also intensifies [40]. The cells are broken and can effectively improve the diffusivity and extraction yield. In addition to that the high ultrasound intensity may damage the structure of polysaccharides reducing the yield [41]. Therefore, the extraction power is selected 150 W.
3.2. Optimization of Sap extracting conditions
Based on the results of single factor experiments. The main independent variables (A: ultrasonic power, B: ultrasonic time, and C: ultrasonic temperature) on the yield of Sap was displayed in Table 1. The yields of Sap ranged from 13.99–21.09% in the 17 experiments, such relatively great difference in extraction yield could be because the surface damage degree of the sample is different with the change of ultrasonic condition, and the amount of polysaccharide from the particles to the solution is also changed [2]. Thus, the optimized ultrasonic extraction parameter are required to improve the yield of Sap. In this study, multiple regression analysis and a second-order polynomial equation, generated by Design-Expert 12.0 software, were used to explain the relationship among three factors and responses. The final simplified equation for Sap extraction yield Y1 in terms of coded parameters is:
$${\text{Y}}_{\text{1}}\text{=20.06+1.59 }\text{A}\text{+0.65 }\text{B}\text{-1.38 }\text{C}\text{+0.12 }\text{AB}\text{-1.18 }\text{AC}\text{+0.35 }\text{BC}\text{-1.29 }{\text{A}}^{\text{2}}\text{-1.99 }{\text{B}}^{\text{2}}\text{-1.91 }{\text{C}}^{\text{2}}$$
Where Y1, A, B, C is extraction yield, ultrasonic power (W), ultrasonic time (min), and ultrasonic temperature (℃), respectively.
The fitting degree of model plays an important role in statistics. The Variance (ANOVA) and correlation coefficients (For coefficients R2 and adjusted R2) were the main methods of expression. Y1 is highly significant with a high F value (62.77) and low p-values (p < 0.0001), which proved that the regression model is significant (Table S2),. Moreover, the value of correlation coefficient R2 (0.9878) and adjusted R2Adj (0.972) are almost 1, which can explain the model is a high degree of correlation. The affecting of parameters on extraction process: A > C > B (p > 0.05), which might be due to the large polysaccharide particles are decomposed into smaller particles led to improving the yield under the diffusion of strong splitting force [2]. In terms of the interactions, the synergistic (positive) effect between the ultrasonic power and temperature (AC) is the most prominent (p < 0.01), which can improve the polysaccharide yield. This may be because both factors (AC) possess the capacity of destructive. As the power and temperature increases, yielded higher overall driving force, which led to the effect of extraction time is more significant. This means that the mild extraction conditions (low level) may yield a positive effect. When the extraction conditions are harsher, the effect of AC as the dominant factor. On the other hand, the coefficient variation (C.V.: 2.24%) and Adeq. precision (23.79) demonstrated experimental values shows high reliability and precision. The linear coefficient (AB, BC: >0.05) were non-significant, and the term coefficients (A2, B2, C2: <0.05) was found to be at a significant.
3.3. Response surface analysis
The 3D response surface and two-dimension contour plots are mainly applied to express the interactions between the dependent /independent variables when other variables are fixed at zero level. The trend of response values kept increasing with increase values of variables, and then decreased, that illustrates the response surface is steady (Fig. S3). The effect of interactions extraction time and power had no significant effect on extraction rates when the temperature was fixed. The yield of Sap gradually increased when time increased from 10 to 15.68 min and power increased from 120 to 177.11 W (Fig. S3A and S3B). The further increase of extraction time and power had the negative effect on Sap yield. The extraction temperature and power exhibited quadratic effects on the Sap yields (Fig. S3C and S3D). The extraction rate increased significantly with the increased of extraction temperature over the range of 75 ℃-76.86 ℃ and power increased over the range of 120-177.11 W. The extraction yield didn’t increase when extraction time and power reached the maximum yield, which may be due to the polysaccharides degraded under the synergistic effect of higher temperature and higher ultrasonic power. The interactions of extraction temperature and time on the extraction yield of Sap. The extraction rate increased when the temperature and time increased in the ranges of 75 ℃ to 76.86 ℃ and 10 to 15.68 min, respectively (Fig. S3E and S3F). Then, the yield declined gradually when the extraction temperature and time exceeded the ranges. This could be because polysaccharides were decomposed under higher temperature and longer extraction time. Based on the prediction of mathematical model, the optimal conditions for obtaining the maximum yield of Sap (21.24%) were as follows: 177.11 W, 15.68 min, 40 mL/g, 76.86 ℃. To simplify the actual operation parameters, extraction power, temperature, time, and liquid to material ratio set at 180 W, 77 ℃, 16 min; 40 mL/g.
3.4. Physicochemical properties of Sap
After analysis, the total carbohydrate content of the Sap extracted from the S sagittifolia L. was 79.26 ± 0.39% as determined by the phenol‑sulfuric acid colorimetric method (Table 2). Moreover, the polyphenolic composition has been shown to affect the antioxidant activity of active substance, we also investigated whether the polyphenolic composition of Sap affects their antioxidant capacity. But no phenolic content was detected in Sap. Protein content was determined to be 2.67 ± 0.06%, these findings are consistent with the UV weak absorption peak was recorded at 280 nm for protein (Fig .1A). Additionally, the uronic acid was also detected, the content of uronic acid in Sap 1.33 ± 0.24%. Mw play an important role in biological activities of polysaccharides. The Mw and the molecular weight distribution coefficients (Mw/Mn) of Sap was 120.5 kDa and 1.591 (Fig. 1B). Moreover, it was shown that when the polydispersity index of the samples are smaller than 2 could imply that samples is less likely to form large sized agglomerate in aqueous solution [18]. The monosaccharide composition of Sap was determined to be mainly composed of rhamnose, arabinose, mannose, glucose, galactose and the molar ratio was 15.7: 12.9: 1.6: 36.4: 33.4 in Sap (Fig. 1C, Table 2). The Sap has the characteristic absorption peaks of carbohydrates as determined by FT-IR (Fig. 1D). In detail, the signal at 3422 cm− 1 (O-H stretching vibration) and 2933 cm− 1 (C-H stretching vibration) represents the characteristic groups of polysaccharides [24]. In particular, the absorption at 1735 cm− 1 and 1654 cm− 1 is related to the uronic acid (C = O stretching vibrations) and proteins, respectively. The band at 1000–1200 cm− 1 are attributed to the bond stretching vibrations of C-O-C and O-C-O bonds reflects the presence of pyranose [42]. The signal at 900 cm− 1 and 836 cm− 1 represents the feature of β-glycosidic linkages and α-type glycosidic linkages, respectively [43].
Table 2
Total carbohydrate (%)
|
Total phenolic
|
Protein
(%)
|
Uronic acid (%)
|
Mw
(kDa)
|
Mn
|
Mw/Mn
|
Monosaccharide compositions (mol%)
|
Rha
|
Ara
|
Man
|
Glc
|
Gal
|
79.26 ± 0.39
|
N. D
|
2.67 ± 0.06
|
1.33 ± 0.24
|
120.5
(± 0.804%)
|
75.76
(± 0.317%)
|
1.591
(± 0.865%)
|
15.7
|
12.9
|
1.6
|
36.4
|
33.4
|
Note: Each value represents the mean ± standard deviation (n = 3). N.D.: Not detected or below the limit of quantification. (Rha:Rhamnose, Ara:Arabinose, Man:Mannose, Glc:Glucose, Gal:Galactose) |
The morphology and chain conformation are a critical indicator of polysaccharides involved in a variety of biological activities. We made use of AFM to visualize Sap has a compact globular architecture and a vast number of uneven lumps structure (Fig. 2A). The characteristic structures of Sap were examined by SEM image of 200 and 1000-fold magnification and its displayed as non-uniform and broken microsphere fragment’s structure (Fig. 2B). The action of ultrasonic cavitation may be destructive to the microstructure of samples to expose more active sites, which may enhance biological activity.
3.5. Analysis of spectral preprocessing
The DPPH, ABTS and hydroxyl radical scavenging rate increased with the increase of Sap concentration (Fig. S3). The free radical scavenging rate trends of Sap at varied concentrations was in agreement with the reported results [18]. It can be observed that the original spectra of Sap changed when different preprocessing were employed (Fig. 3, Table 3). The correlation coefficients of the models of DPPH antioxidant activities were significantly improved by using MSC preprocessing. The reason for this phenomenon is because the scattered light of different particle sizes was removed by linear fitting each individual and reference spectrum. Therefore, the optimum PLS model was achieved when 8 PLS factors and MSC spectra preprocessing were included in this work for DPPH (Rp = 0.9568 and RMSEP = 2.319%). The prediction by PLS model based on original spectra was superior to that of spectral preprocessing with higher Rp of 0.9651, 0.9602 and lower RMSEP of 6.160% and 4.265% for ABTS and hydroxyl radical, respectively.
Table 3
The results of PLS model based on different spectral preprocessing methods.
Parameter
|
spectral type
|
|
Calibration set
(n = 90)
|
Cross-validation
|
Prediction set
(n = 30)
|
LVs
|
Rc
|
RMSEC
|
Rcv
|
RMSECV
|
Rp
|
RMSEP
|
|
Original
|
7
|
0.9686
|
2.01
|
0.9392
|
2.902
|
0.9410
|
2.710
|
DPPH
|
SNV
|
8
|
0.9850
|
1.459
|
0.9085
|
3.573
|
0.9525
|
2.430
|
MSC
|
8
|
0.9853
|
1.442
|
0.9049
|
3.643
|
0.9568
|
2.319
|
Baseline
|
8
|
0.9697
|
2.062
|
0.9413
|
2.855
|
0.9386
|
2.756
|
Derivative I
|
6
|
0.9775
|
1.779
|
0.9116
|
3.475
|
0.9060
|
3.535
|
Derivative II
|
4
|
0.9317
|
3.068
|
0.8232
|
4.796
|
0.8678
|
4.067
|
Original
|
6
|
0.9689
|
5.497
|
0.9457
|
7.231
|
0.9651
|
6.160
|
SNV
|
6
|
0.9696
|
5.435
|
0.9349
|
7.901
|
0.9440
|
7.652
|
ABTS
|
MSC
|
5
|
0.9579
|
6.381
|
0.9217
|
8.649
|
0.9319
|
8.482
|
Baseline
|
6
|
0.9635
|
5.946
|
0.9384
|
7.700
|
0.9604
|
6.413
|
Derivative I
|
5
|
0.9799
|
4.431
|
0.9501
|
6.946
|
0.9509
|
7.225
|
Derivative II
|
4
|
0.9519
|
6.81
|
0.8744
|
10.79
|
0.8832
|
11.29
|
Original
|
6
|
0.9679
|
3.891
|
0.9486
|
4.904
|
0.9602
|
4.265
|
SNV
|
6
|
0.9680
|
3.886
|
0.9336
|
5.559
|
0.9465
|
4.966
|
MSC
|
6
|
0.9689
|
3.831
|
0.9252
|
5.908
|
0.9430
|
5.177
|
Hydroxyl radical
|
Baseline
|
6
|
0.9647
|
4.079
|
0.9418
|
5.217
|
0.9550
|
4.569
|
Derivative I
|
5
|
0.9696
|
3.789
|
0.9292
|
5.722
|
0.9253
|
6.014
|
Derivative II
|
4
|
0.9394
|
5.307
|
0.8392
|
8.431
|
0.8570
|
8.432
|
3.6. NIR calibration results
3.6.1 Analysis of DPPH
As shown in Table 5, the model for DPPH scavenging activity showed an Rc of 0.9853, RMSECV of 1.442%, Rp of 0.9568 and RMSEP of 2.319% in the MSC-PLS model. The iPLS based model for DPPH scavenging activity, with 20 intervals for the full spectrum. Then the model of each subinterval was established. The predictive MSC-iPLS model of DPPH, the model built on the 19 subintervals of the whole spectral interval was the best, which corresponds to the spectral interval of 1358.34-1439.97 nm, with the highest value of Rp =0.81 and RMSEP = 4.82% in the prediction.
The results showed that Si-PLS model of DPPH scavenging activity showed that the optimal number of LVs selected from the 4 intervals from 15 (Table
4). The optimal PCs was 9, and combinations of intervals selected were [
1,
5,
10,
13], which corresponding to 868.6-977.4, 1303.8-1412.6, 1847.8-2065.4 and 2174.2–2283.0 nm in the spectral regions (Fig.
4A), the Si-PLS model yielded
Rc= 0.9459,
RMSEC = 2.76%,
Rp= 0.9455 and
RMSEP = 2.70%. The spectra regions are features of various unique constituents, mainly related to the combination of some functional group and overtones and vibrational modes such as -CH, -NH, -SH and -OH groups [
44,
45]. We observed that the spectra regions 868.6-977.4 nm commonly associated to third C-H overtone of polysaccharides. At 1303.8-1412.6 nm, speculate that the variation of is related to the O-H first overtone. At 1847.8-2065.4 nm is due to the O-H stretching vibration and at 2174.2–2283.0 nm are derived from the C-H stretching vibration and CH
2 deformation. Compared with Si-PLS model of full spectra, MSC-PLS model showed good prediction effect. Figure
4D is the scatter plots that show DPPH estimated by MSC-PLS model in the prediction sets.
Table 4
Spectral subinterval of Si-PLS calibration model.
Parameters
|
Number of subintervals
|
PCs
|
Selected subintervals
|
RC
|
RMSEC
|
RP
|
RMSEP
|
DPPH
|
15
|
9
|
[ 1, 5, 10, 13]
|
0.9459
|
2.76
|
0.9455
|
2.70
|
|
16
|
11
|
[ 5, 10, 11, 16]
|
0.9321
|
3.08
|
0.8837
|
4.06
|
17
|
7
|
[ 1, 8, 11, 16]
|
0.9378
|
2.94
|
0.9361
|
2.95
|
18
|
9
|
[ 6, 12, 15, 16]
|
0.931
|
3.10
|
0.8855
|
3.81
|
19
|
10
|
[ 2, 5, 6, 12]
|
0.9347
|
3.01
|
0.8842
|
3.8
|
20
|
9
|
[ 7, 13, 16, 17]
|
0.9375
|
2.95
|
0.9132
|
3.31
|
21
|
9
|
[ 7, 13, 14, 18]
|
0.9246
|
3.24
|
0.9197
|
3.22
|
22
|
10
|
[ 6, 7, 14, 18]
|
0.945
|
2.78
|
0.8928
|
3.84
|
23
|
12
|
[ 7, 15, 16, 23]
|
0.9462
|
2.75
|
0.8389
|
5.54
|
24
|
9
|
[ 8, 15, 16, 20]
|
0.9346
|
3.01
|
0.9306
|
3.17
|
25
|
9
|
[ 8, 16, 17, 21]
|
0.9388
|
3.02
|
0.9259
|
3.21
|
26
|
7
|
[ 1, 13, 16, 24]
|
0.9404
|
2.72
|
0.9360
|
2.97
|
27
|
10
|
[ 8, 17, 18, 22]
|
0.9357
|
3.01
|
0.9249
|
3.22
|
28
|
9
|
[ 9, 18, 19, 23]
|
0.9354
|
3.00
|
0.9163
|
3.38
|
29
|
6
|
[ 1, 10, 14, 17]
|
0.9248
|
3.21
|
0.9251
|
3.03
|
30
|
8
|
[ 9, 18, 19, 25]
|
0.9345
|
3.01
|
0.9226
|
3.07
|
ABTS
|
15
|
12
|
[ 2, 6, 10, 14]
|
0.9771
|
4.73
|
0.9214
|
9.06
|
|
16
|
10
|
[ 1, 6, 11, 16]
|
0.9717
|
5.25
|
0.9617
|
7.57
|
17
|
10
|
[ 2, 7, 11, 17]
|
0.9716
|
5.27
|
0.9542
|
7.46
|
18
|
9
|
[ 2, 7, 12, 17]
|
0.9753
|
4.92
|
0.9351
|
8.12
|
19
|
8
|
[ 2, 7, 9, 12]
|
0.9727
|
5.16
|
0.9607
|
6.4
|
20
|
12
|
[ 2, 8, 17, 20]
|
0.9702
|
5.42
|
0.9514
|
7.02
|
21
|
9
|
[ 2, 8, 9, 21]
|
0.971
|
5.32
|
0.9534
|
7.3
|
22
|
9
|
[ 2,8, 9, 22]
|
0.9704
|
5.37
|
0.9507
|
7.45
|
23
|
6
|
[ 2, 9, 12, 15]
|
0.9682
|
5.56
|
0.9457
|
7.58
|
24
|
7
|
[ 2, 7, 9, 15]
|
0.9704
|
5.37
|
0.9408
|
8.09
|
25
|
9
|
[ 1, 12, 17, 24]
|
0.9693
|
5.47
|
0.9667
|
6.75
|
26
|
11
|
[ 3, 10. 17. 24]
|
0.9664
|
5.72
|
0.9016
|
10.6
|
27
|
10
|
[ 3, 10, 12, 17]
|
0.9725
|
5.19
|
0.9352
|
8.31
|
28
|
9
|
[ 2, 18, 19, 27]
|
0.967
|
5.66
|
0.9416
|
9.24
|
29
|
10
|
[ 3, 11, 19, 27]
|
0.9669
|
5.68
|
0.9291
|
8.89
|
30
|
9
|
[ 3, 11, 13, 19]
|
0.9706
|
5.41
|
0.9316
|
8.6
|
Hydroxyl radical
|
15
|
10
|
[ 1, 7, 10, 15]
|
0.9644
|
4.11
|
0.9542
|
4.67
|
|
16
|
9
|
[ 1, 8, 11, 15]
|
0.9507
|
4.83
|
0.9614
|
4.29
|
17
|
9
|
[ 1, 8, 11, 17]
|
0.9693
|
3.81
|
0.9617
|
4.23
|
18
|
8
|
[ 2, 5, 7, 9]
|
0.969
|
3.83
|
0.958
|
4.37
|
19
|
9
|
[ 1, 9, 12, 19]
|
0.9696
|
3.79
|
0.9437
|
5.15
|
20
|
8
|
[ 2, 8, 9, 17]
|
0.9691
|
3.82
|
0.9532
|
4.7
|
21
|
8
|
[ 2, 6, 8, 10]
|
0.9702
|
3.75
|
0.9526
|
4.65
|
22
|
11
|
[ 2, 6, 13, 19]
|
0.9681
|
3.89
|
0.9043
|
7.21
|
23
|
9
|
[ 2, 9, 10, 16]
|
0.9688
|
3.84
|
0.9652
|
4.22
|
24
|
7
|
[ 2, 7, 9, 12]
|
0.9709
|
3.71
|
0.955
|
4.55
|
25
|
10
|
[ 2, 6, 11, 17]
|
0.9669
|
3.98
|
0.9537
|
4.93
|
26
|
8
|
[ 2, 6, 17, 18]
|
0.9655
|
4.04
|
0.9426
|
6.13
|
27
|
8
|
[ 2, 10, 11, 18]
|
0.9705
|
3.74
|
0.9568
|
4.84
|
28
|
9
|
[ 2, 3, 12, 19]
|
0.9714
|
0.368
|
0.9651
|
4.91
|
29
|
7
|
[ 2, 8, 11, 14]
|
0.9646
|
4.09
|
0.9582
|
4.46
|
30
|
8
|
[ 2, 3, 13, 20]
|
0.9668
|
3.96
|
0.9532
|
5.01
|
Table 5
Comparison of different algorithms.
Parameter
|
Spectral type
|
LVs
|
Rc
|
RMSEC
|
Rp
|
RMSEP
|
DPPH
|
MSC-PLS
|
8
|
0.9853
|
1.442
|
0.9568
|
2.319
|
MSC-iPLS
|
5
|
0.8376
|
4.61
|
0.8100
|
4.82
|
MSC-Si-PLS
|
9
|
0.9459
|
2.76
|
0.9455
|
2.70
|
ABTS
|
PLS
|
6
|
0.9689
|
5.497
|
0.9651
|
6.160
|
iPLS
|
6
|
0.8836
|
10.40
|
0.8747
|
12.00
|
Si-PLS
|
9
|
0.9693
|
5.47
|
0.9667
|
6.75
|
Hydroxyl radical
|
PLS
|
6
|
0.9679
|
3.891
|
0.9602
|
4.265
|
iPLS
|
5
|
0.8737
|
7.53
|
0.8712
|
7.59
|
Si-PLS
|
9
|
0.9688
|
3.84
|
0.9652
|
4.22
|
3.6.2 Analysis of ABTS
As for ABTS scavenging activity, PLS model yielded Rc= 0.9689, RMSECV = 5.497% and Rp= 0.9651 and RMSEP = 6.16%. In iPLS model, the optimal model was built on the 19 subintervals of the full spectral interval, which corresponds to the spectral interval of 1848.08-1929.70 nm, with Rc= 0.8836, RMSECV = 10.40% and Rp= 0.8747, RMSEP = 12% (Table 5). As shown in Table 4, the optimal spectral intervals selected by Si-PLS model for the prediction of ABTS. The optimal parameter of Si-PLS model was a combination of 4 from 25 intervals, the optimal LVs was 9, and combination of intervals selected were [1, 12, 17, 24], corresponding to 868.6-933.9, 1586.9-1652.2, 1913.4-1978.7 and 2370.5-2435.8 nm (Fig. 4B). The spectra regions at 868.6-933.9 nm are connected with the third C-H overtone. 1586.9-1652.2 nm is due to the first C-H overtone. The spectra regions around 1913.4-1978.7 nm were assigned to O-H stretching vibration, and that at 2370.5-2435.8 nm were assigned to C-H stretching vibration and CH2 deformation. Figure 4E was the scatter plot showing ABTS estimated by Si-PLS model in the prediction sets. The highest Rp= 0.9667 and lowest RMSEP = 6.75%.
3.6.3 Analysis of Hydroxyl radical
As for the hydroxyl radical inhibition activity, the PLS model yielded Rc= 0.9679, RMSECV = 3.891% and Rp= 0.9602, RMSEP = 4.265%. In the iPLS model, the best model was built on 19 subintervals of the whole spectral interval, which corresponds to the spectral interval of 1358.34-1439.97 nm, with Rc= 0.8737, RMSECV = 7.53% and Rp= 0. 8712, RMSEP = 7.59% (Table 5). The optimal parameter of Si-PLS model was a combination of 4 intervals from 23 interval, the optimal LVs was 9, and the combinations of intervals selected were [2, 9, 10, 16], corresponding to 939.5-1010.4 nm, 1435.8-1506.7 nm, 1506.7-1577.6 nm and 1932.1–2003.0 nm (Fig. 4C), respectively. The spectra regions 939.5-1010.4 nm is mainly attributed to the second O-H overtone. 1435.8-1506.7 nm and 1506.7-1577.6 nm is due to the CH3, CH2, CH of polysaccharides. The spectra regions around 1932.1–2003.0 nm contributed from O-H stretching vibration. The scatter plot showing hydroxyl radical estimated by Si-PLS model in the prediction sets (Fig. 4F). Here, the Rp= 0.9652 and RMSEP = 4.22%. Compared with PLS and iPLS method, the Si-PLS models of values of Rp and Rc are higher, and RMSECV and RMSEP are lower.