# Determination of Total Saccharide Content inAuricularia auricula Based on Near-Infrared Spectroscopy

Authors: Gen Yang; Shanlei Li; Shixian Ji; Yunjie Wang; Jinmei Wang; Liqiang Ji; Changqin Li
Journal: Journal of Food Quality
 (2022)

---

## Abstract

In this paper, the content of total saccharide inAuricularia auricula from different regions was determined. Then, near-infrared (NIR) technology was used to collect the spectral information of the samples. The sample data were divided into calibration set and validation set. The best quantitative model of the total saccharide content of A. auricula was established by selecting the parameters such as spectral range, pretreatment method, and partial least square method (PLS) main factor number of the calibration set data. The validation set data were used to verify the reliability of this model. In this model, the original spectrum was used to preprocess by standard normal variate (SNV) + second derivative (SD) to eliminate the scattering effect caused by uneven particle distribution and the influence of noise on spectral data. The spectrum range was 4000–10000 cm−1, and the final choice of PLS main factor number was 11. Under this condition, the calibration set Rc2 of the model was 0.9092, the root mean square error of calibration (RMSEC) was 1.405, the root mean square error of prediction (RMSEP) was 1.507, and the residual predictive deviation (RPD) was 3.32. The validation samples were used to test the model, and the result showed that Rv2 = 0.9048 of the validation set. The result proved that the predicted value of the validation samples had a good linear relationship with the measured value. According to the T-test of the two sets of data in the validation set, the difference between the predicted value and the chemical value was not significant (P ≥ 0.05). The results were in line with the expected objectives. The established NIR quantitative model can be used to predict the total saccharide content of the black fungus sample to be tested.


---

## Body

## 1. Introduction

Auricularia auricula is the fruiting body of edible fungus used as medicine and food [1,2], which contains active ingredients such as cellulose, polysaccharide, and minerals such as calcium, iron, and phosphorus [3]. A. auricula polysaccharide (AAP) is one of the main effective components and has biological activities such as antitumor, lowering blood glucose, anticoagulant, antiaging, and enhancing immunity [4–7]. Therefore, A. auricula has great potential for medicinal development. At present, the methods of AAP extraction include water extraction, acid-base extraction, ultrasonic-assisted extraction, and enzyme extraction, but these methods usually take a lot of time and manpower and damage the original medicine to a certain extent [8]. Because of the complex composition of polysaccharides, it is difficult to determine them directly. In practice, the content of total saccharide is often used to replace the content of polysaccharide to evaluate the quality of food. So, it is necessary to find a rapid and accurate method to determine the total saccharide content of A. auricula.In recent years, near-infrared (NIR) spectroscopy has attracted more and more attention because of its fast detection speed, environmental protection, and no loss of samples. Especially since the 1990s, NIR had been widely used in industrial fields, such as applied physics, textile, architecture, history, pharmacy, and food [9–15]. In the near-infrared spectrum analysis, the commonly used modeling methods include multiple linear regression (MLR), principal component regression (PCR), partial least squares (PLS), and other linear correction methods, as well as artificial neural network (ANN) and support vector machine (SVM) and other nonlinear correction methods [16–18].In this paper, the PLS method was used as the modeling method to establish the near-infrared total saccharide prediction model ofA. auricula.


## 2. Materials and Methods

### 2.1. Materials

166 batches ofA. auricula were collected from different regions in China during 2018-2019, and the specific information was shown in Table 1. All samples were sealed and stored in a cool and dry place.Table 1
Sample information ofA. auricula.
NumberNameProvinceNumberNameProvince1A. auriculaYunnan33A. auriculaHeilongjiang2A. auriculaGuangdong34A. auriculaHeilongjiang3A. auriculaGuangdong35A. auriculaZhejiang4A. auriculaHeilongjiang36A. auriculaSichuan5A. auriculaYunnan37A. auriculaJilin6A. auriculaSichuan38A. auriculaJilin7A. auriculaJilin39A. auriculaHeilongjiang8A. auriculaHeilongjiang40A. auriculaZhejiang9A. auriculaHenan41A. auriculaHubei10A. auriculaHenan42A. auriculaShandong11A. auriculaFujian43A. auriculaZhejiang12A. auriculaHebei44A. auriculaHeilongjiang13A. auriculaYunnan45A. auriculaHeilongjiang14A. auriculaHebei46A. auriculaGuizhou15A. auriculaNingxia47A. auriculaSichuan16A. auriculaFujian48A. auriculaSichuan17A. auriculaShanxi49A. auriculaSichuan18A. auriculaJilin50A. auriculaSichuan19A. auriculaShandong51A. auriculaChongqing20A. auriculaAnhui52A. auriculaHeilongjiang21A. auriculaHenan53A. auriculaHeilongjiang22A. auriculaFujian54A. auriculaHeilongjiang23A. auriculaJilin55A. auriculaAnhui24A. auriculaHeilongjiang56A. auriculaHeilongjiang25A. auriculaHeilongjiang57A. auriculaJilin26A. auriculaHeilongjiang58A. auriculaAnhui27A. auriculaHeilongjiang59A. auriculaJilin28A. auriculaHeilongjiang60A. auriculaJilin29A. auriculaHeilongjiang61A. auriculaSichuan30A. auriculaHeilongjiang62–69A. auriculaShanghai31A. auriculaHeilongjiang70–166A. auriculaJilin32A. auriculaHeilongjiang


### 2.2. Instruments and Equipment

Spot light 400 mid/near-infrared spectrometer (PerkinElmer Inc, U.S.); rotary evaporator (EYELA Inc., Japan); AB135-s 100000 level electronic balance (Mettler Toledo Inc, Switzerland); UV-2600 ultraviolet visible spectrophotometer (Shimadzu Inc, Japan); DHG-9053A electric blast drying oven (Yi Heng, Shanghai, China); HH-S4 digital display constant temperature water bath (Xiang Tian, Changzhou, China); and SK-1 fast mixer (Hua Feng, Jintan, China).


### 2.3. Determination Method of Total Saccharide Content

#### 2.3.1. Preparation of Test Solution

The sample ofA. auricula was ground, passed through a 20 mesh sieve, and then 0.25 g of A. auricula powder was precisely weighed and placed in a 250 mL conical flask. 50 mL water and 15 mL concentrated hydrochloric acid were added, heated, and refluxed in a 100°C water bath for 3 h. After cooling to room temperature, the filter residue was washed with distilled water, the filtrate and washing solution were combined, and the volume of water was fixed to 250 mL to prepare the test solution.


#### 2.3.2. Drawing of Standard Curve

0.1, 0.2, 0.4, 0.6, 0.8, and 1.0 mL of 100 mg/L glucose standard solution were accurately taken to 10 mL test tube with plug, and then distilled water was added to make up to 10 mL. 1.0 mL of 5% phenol solution was added into the test solution, and then 5.0 mL of concentrated sulfuric acid was added. The reaction solution was mixed by using a rapid mixer. Then, the test tube was placed in a water bath at 30°C for 20 min, and an appropriate amount of reaction solution was taken to measure the absorbance at 490 nm. The standard curve was drawn, and the standard curve equation was calculated with the absorbance value as the ordinate and the glucose concentration (mg/mL) as the abscissa.


#### 2.3.3. Determination of Sugar Content in Samples

0.2 ml of the test solution was accurately piped into a 10 mL tube with plug, and distilled water was added to make up to 1.0 mL. The absorbance was measured according to the above steps. The total saccharide contentw was calculated by the standard curve, and the value was expressed in percentage (%). The calculation formula is as follows:(1)w=m1×v1×10−6m2×v2×1−ω×100,where v1 is the constant volume of sample, V2 is the volume of the sample solution taken during colorimetric determination, m1 is the sugar content in the sample solution which can be obtained from the standard curve, m2 is the sample quality, and ω is the water content of the sample.


#### 2.3.4. The Method of Near-Infrared Spectrum Collection

166 batches ofA. auricula powder samples (80 mesh) were spread separately on the integrating sphere of the NIR spectrometer, and the samples should be completely covered with the integrating sphere. The scanning spectrum range was 4000–10000 cm−1, the scanning times were 64, and the resolution was 8 cm−1. The original spectral data of sample were shown in Figure 1.Figure 1
NIRS spectrum of (A) auricula.


### 2.4. Method and Verification of Establishing Quantitative Model

The spectral data of samples were analyzed by Spectrum Quant software, and partial least square (PLS) method was selected as the modeling method. Four indexes, RMSEC (root mean square error of calibration), RMSEP (root mean square error of prediction),R2 (coefficient of determination), and RPD (residual predictive deviation), were used to evaluate the performance of the model. Generally, the larger the R2, the smaller the RMSEC and RMSEP, and the better the model prediction effect. And, the smaller the difference between RMSEC and RMSEP, the better the generalization ability of the prediction model. Generally, the ratio of RMSEC/RMSEP needs to be controlled between 0.8 and 1.2. RPD is the final evaluation index of the model, and its evaluation standard adopts the threshold segmentation method [19–21]. When RPD <2.0, it indicates that the model is very poor and cannot be applied; when RPD is 2.0–2.5, it indicates that the model can make a rough prediction and correlation evaluation of the samples; when RPD is 2.5–3.0, it indicates that the model is better and can be used for quantitative analysis of samples; when RPD> 3.0, it indicates that the model has excellent predictive ability and can be used for model analysis [22]. The calculation formula of RPD is as follows:(2)RPD=11−R2.


## 2.1. Materials

166 batches ofA. auricula were collected from different regions in China during 2018-2019, and the specific information was shown in Table 1. All samples were sealed and stored in a cool and dry place.Table 1
Sample information ofA. auricula.
NumberNameProvinceNumberNameProvince1A. auriculaYunnan33A. auriculaHeilongjiang2A. auriculaGuangdong34A. auriculaHeilongjiang3A. auriculaGuangdong35A. auriculaZhejiang4A. auriculaHeilongjiang36A. auriculaSichuan5A. auriculaYunnan37A. auriculaJilin6A. auriculaSichuan38A. auriculaJilin7A. auriculaJilin39A. auriculaHeilongjiang8A. auriculaHeilongjiang40A. auriculaZhejiang9A. auriculaHenan41A. auriculaHubei10A. auriculaHenan42A. auriculaShandong11A. auriculaFujian43A. auriculaZhejiang12A. auriculaHebei44A. auriculaHeilongjiang13A. auriculaYunnan45A. auriculaHeilongjiang14A. auriculaHebei46A. auriculaGuizhou15A. auriculaNingxia47A. auriculaSichuan16A. auriculaFujian48A. auriculaSichuan17A. auriculaShanxi49A. auriculaSichuan18A. auriculaJilin50A. auriculaSichuan19A. auriculaShandong51A. auriculaChongqing20A. auriculaAnhui52A. auriculaHeilongjiang21A. auriculaHenan53A. auriculaHeilongjiang22A. auriculaFujian54A. auriculaHeilongjiang23A. auriculaJilin55A. auriculaAnhui24A. auriculaHeilongjiang56A. auriculaHeilongjiang25A. auriculaHeilongjiang57A. auriculaJilin26A. auriculaHeilongjiang58A. auriculaAnhui27A. auriculaHeilongjiang59A. auriculaJilin28A. auriculaHeilongjiang60A. auriculaJilin29A. auriculaHeilongjiang61A. auriculaSichuan30A. auriculaHeilongjiang62–69A. auriculaShanghai31A. auriculaHeilongjiang70–166A. auriculaJilin32A. auriculaHeilongjiang


## 2.2. Instruments and Equipment

Spot light 400 mid/near-infrared spectrometer (PerkinElmer Inc, U.S.); rotary evaporator (EYELA Inc., Japan); AB135-s 100000 level electronic balance (Mettler Toledo Inc, Switzerland); UV-2600 ultraviolet visible spectrophotometer (Shimadzu Inc, Japan); DHG-9053A electric blast drying oven (Yi Heng, Shanghai, China); HH-S4 digital display constant temperature water bath (Xiang Tian, Changzhou, China); and SK-1 fast mixer (Hua Feng, Jintan, China).


## 2.3. Determination Method of Total Saccharide Content

### 2.3.1. Preparation of Test Solution

The sample ofA. auricula was ground, passed through a 20 mesh sieve, and then 0.25 g of A. auricula powder was precisely weighed and placed in a 250 mL conical flask. 50 mL water and 15 mL concentrated hydrochloric acid were added, heated, and refluxed in a 100°C water bath for 3 h. After cooling to room temperature, the filter residue was washed with distilled water, the filtrate and washing solution were combined, and the volume of water was fixed to 250 mL to prepare the test solution.


### 2.3.2. Drawing of Standard Curve

0.1, 0.2, 0.4, 0.6, 0.8, and 1.0 mL of 100 mg/L glucose standard solution were accurately taken to 10 mL test tube with plug, and then distilled water was added to make up to 10 mL. 1.0 mL of 5% phenol solution was added into the test solution, and then 5.0 mL of concentrated sulfuric acid was added. The reaction solution was mixed by using a rapid mixer. Then, the test tube was placed in a water bath at 30°C for 20 min, and an appropriate amount of reaction solution was taken to measure the absorbance at 490 nm. The standard curve was drawn, and the standard curve equation was calculated with the absorbance value as the ordinate and the glucose concentration (mg/mL) as the abscissa.


### 2.3.3. Determination of Sugar Content in Samples

0.2 ml of the test solution was accurately piped into a 10 mL tube with plug, and distilled water was added to make up to 1.0 mL. The absorbance was measured according to the above steps. The total saccharide contentw was calculated by the standard curve, and the value was expressed in percentage (%). The calculation formula is as follows:(1)w=m1×v1×10−6m2×v2×1−ω×100,where v1 is the constant volume of sample, V2 is the volume of the sample solution taken during colorimetric determination, m1 is the sugar content in the sample solution which can be obtained from the standard curve, m2 is the sample quality, and ω is the water content of the sample.


### 2.3.4. The Method of Near-Infrared Spectrum Collection

166 batches ofA. auricula powder samples (80 mesh) were spread separately on the integrating sphere of the NIR spectrometer, and the samples should be completely covered with the integrating sphere. The scanning spectrum range was 4000–10000 cm−1, the scanning times were 64, and the resolution was 8 cm−1. The original spectral data of sample were shown in Figure 1.Figure 1
NIRS spectrum of (A) auricula.


## 2.3.1. Preparation of Test Solution

The sample ofA. auricula was ground, passed through a 20 mesh sieve, and then 0.25 g of A. auricula powder was precisely weighed and placed in a 250 mL conical flask. 50 mL water and 15 mL concentrated hydrochloric acid were added, heated, and refluxed in a 100°C water bath for 3 h. After cooling to room temperature, the filter residue was washed with distilled water, the filtrate and washing solution were combined, and the volume of water was fixed to 250 mL to prepare the test solution.


## 2.3.2. Drawing of Standard Curve

0.1, 0.2, 0.4, 0.6, 0.8, and 1.0 mL of 100 mg/L glucose standard solution were accurately taken to 10 mL test tube with plug, and then distilled water was added to make up to 10 mL. 1.0 mL of 5% phenol solution was added into the test solution, and then 5.0 mL of concentrated sulfuric acid was added. The reaction solution was mixed by using a rapid mixer. Then, the test tube was placed in a water bath at 30°C for 20 min, and an appropriate amount of reaction solution was taken to measure the absorbance at 490 nm. The standard curve was drawn, and the standard curve equation was calculated with the absorbance value as the ordinate and the glucose concentration (mg/mL) as the abscissa.


## 2.3.3. Determination of Sugar Content in Samples

0.2 ml of the test solution was accurately piped into a 10 mL tube with plug, and distilled water was added to make up to 1.0 mL. The absorbance was measured according to the above steps. The total saccharide contentw was calculated by the standard curve, and the value was expressed in percentage (%). The calculation formula is as follows:(1)w=m1×v1×10−6m2×v2×1−ω×100,where v1 is the constant volume of sample, V2 is the volume of the sample solution taken during colorimetric determination, m1 is the sugar content in the sample solution which can be obtained from the standard curve, m2 is the sample quality, and ω is the water content of the sample.


## 2.3.4. The Method of Near-Infrared Spectrum Collection

166 batches ofA. auricula powder samples (80 mesh) were spread separately on the integrating sphere of the NIR spectrometer, and the samples should be completely covered with the integrating sphere. The scanning spectrum range was 4000–10000 cm−1, the scanning times were 64, and the resolution was 8 cm−1. The original spectral data of sample were shown in Figure 1.Figure 1
NIRS spectrum of (A) auricula.


## 2.4. Method and Verification of Establishing Quantitative Model

The spectral data of samples were analyzed by Spectrum Quant software, and partial least square (PLS) method was selected as the modeling method. Four indexes, RMSEC (root mean square error of calibration), RMSEP (root mean square error of prediction),R2 (coefficient of determination), and RPD (residual predictive deviation), were used to evaluate the performance of the model. Generally, the larger the R2, the smaller the RMSEC and RMSEP, and the better the model prediction effect. And, the smaller the difference between RMSEC and RMSEP, the better the generalization ability of the prediction model. Generally, the ratio of RMSEC/RMSEP needs to be controlled between 0.8 and 1.2. RPD is the final evaluation index of the model, and its evaluation standard adopts the threshold segmentation method [19–21]. When RPD <2.0, it indicates that the model is very poor and cannot be applied; when RPD is 2.0–2.5, it indicates that the model can make a rough prediction and correlation evaluation of the samples; when RPD is 2.5–3.0, it indicates that the model is better and can be used for quantitative analysis of samples; when RPD> 3.0, it indicates that the model has excellent predictive ability and can be used for model analysis [22]. The calculation formula of RPD is as follows:(2)RPD=11−R2.


## 3. Results

### 3.1. Determination Results of Total Saccharide Content of A. auricula

According to formula (1), the measurement results of the total saccharide content in 166 batches were shown in Table 2.Table 2
The results of total saccharide content ofA. auricula.
NumberTotal saccharide content (%) (x ± s, n = 3)114.54 ± 0.37215.28 ± 0.12312.03 ± 1.45412.60 ± 1.20518.37 ± 3.31613.60 ± 2.86710.79 ± 1.25811.07 ± 0.97915.42 ± 3.021018.60 ± 1.161113.14 ± 2.07129.26 ± 0.621311.82 ± 1.361411.37 ± 0.661511.92 ± 1.281612.23 ± 3.561719.78 ± 0.981812.47 ± 2.201918.22 ± 2.032019.36 ± 0.802114.43 ± 0.262217.86 ± 2.132317.82 ± 3.902413.25 ± 1.542514.04 ± 0.27264.52 ± 0.57277.91 ± 3.66286.08 ± 1.00294.01 ± 0.08305.94 ± 1.99319.56 ± 0.44326.31 ± 3.79339.62 ± 1.313411.49 ± 2.603512.33 ± 1.343615.61 ± 0.973712.86 ± 3.313811.91 ± 1.393912.77 ± 2.11406.99 ± 0.014114.91 ± 1.00423.83 ± 1.57439.02 ± 0.47445.41 ± 0.06455.54 ± 0.25464.98 ± 1.04473.33 ± 2.284810.07 ± 1.89494.62 ± 0.135010.67 ± 0.645111.57 ± 2.88525.75 ± 1.31534.63 ± 0.145410.01 ± 2.13552.85 ± 1.50563.59 ± 1.36572.99 ± 1.135811.85 ± 0.41595.25 ± 1.536011.60 ± 1.126120.93 ± 0.416219.19 ± 3.296319.47 ± 3.696414.07 ± 2.52658.33 ± 0.66663.48 ± 1.02677.09 ± 0.76684.78 ± 1.80696.27 ± 0.267014.66 ± 2.937125.24 ± 0.32728.51 ± 1.53735.16 ± 0.197414.52 ± 2.937512.22 ± 0.18765.11 ± 1.917713.56 ± 0.77810.54 ± 0.27797.35 ± 0.06809.38 ± 1.2816.51 ± 0.05828.90 ± 0.468321.48 ± 0.68414.83 ± 1.1854.05 ± 1.488614.43 ± 0.58715.68 ± 0.75889.98 ± 2.36896.65 ± 1.93905.85 ± 1.47915.16 ± 0.049212.29 ± 0.96937.59 ± 0.97948.41 ± 0.22956.69 ± 0.369614.54 ± 1.07979.51 ± 0.73981.30 ± 0.05995.55 ± 0.341007.88 ± 0.521015.73 ± 1.310217.33 ± 1.4810311.61 ± 0.081044.09 ± 0.9510510.52 ± 0.1910612.61 ± 3.131078.88 ± 0.3710818.72 ± 0.0410911.02 ± 1.641109.66 ± 1.7311112.62 ± 0.581122.61 ± 0.6411310.24 ± 0.311143.27 ± 0.231157.62 ± 0.1511612.23 ± 1.701172.67 ± 0.171188.75 ± 3.0711914.86 ± 0.481206.10 ± 0.2912111.01 ± 1.7312213.83 ± 3.3112311.26 ± 2.651245.97 ± 2.5312511.53 ± 2.2012617.25 ± 3.111276.89 ± 0.151284.11 ± 1.281298.19 ± 2.311309.59 ± 2.0213115.18 ± 1.321323.81 ± 2.5113314.84 ± 0.8013414.39 ± 0.4613514.51 ± 1.451363.66 ± 0.1713712.90 ± 0.791382.56 ± 0.151396.38 ± 0.271404.86 ± 0.1114110.27 ± 0.341429.79 ± 0.6714315.98 ± 1.8614414.39 ± 0.811452.11 ± 1.5314614.10 ± 0.214710.61 ± 0.321486.68 ± 1.6814912.00 ± 1.5615016.17 ± 1.291516.48 ± 0.1615214.46 ± 0.6415311.69 ± 1.711549.96 ± 3.0315512.10 ± 0.4315613.69 ± 0.7415710.38 ± 3.2315813.77 ± 0.7815912.01 ± 0.4116011.95 ± 0.5016113.14 ± 0.1816213.85 ± 0.701638.44 ± 2.141644.61 ± 0.361657.04 ± 0.6416615.71 ± 0.09


### 3.2. Removal of Model Abnormal Samples

In this study, software was used to identify the abnormal samples automatically. The PCA score map was used to cluster analyse the spectral data ofA. auricula samples, and the differences among the samples were counted. Diagnosis of abnormal data by test of studentized residual and leverage values method [23]. Finally, 47 batches of abnormal samples were removed from the model. The modeling process was shown in Figures 2 and 3.Figure 2
Test of studentized residual.Figure 3
Leverage values of samples.


### 3.3. Division of Calibration Set and Verification Set

The calibration set and verification set were selected from all samples randomly. The following two conditions must be met during the selection process: ① ensure that the total saccharide content of the verification set samples covers the content gradient of all samples. ② Ensure that the content range of the selected verification set samples does not exceed the calibration set.In Table3, the samples of the verification set were completely distributed within the range of the calibration set, indicating that the division of the sample set meets the modeling requirements.Table 3
Range of total saccharide content in calibration samples and validation samples.
Sample setNumber of samplesMaximum (%)Minimum (%)Correction set10419.18962.5641Validation set1517.33003.5860


### 3.4. Selection of Spectral Range for Modeling

The calibration model is established by choosing appropriate spectral range which has strong prediction ability, and good robustness can often be obtained [24]. In Figure 1, the sample had obvious characteristic peak absorption at 5100 cm−1 and 7000 cm−1. Meanwhile, R2, RMSEC, RMSEP, and RPD were used as model evaluation indicators to establish the model in four spectral intervals. The appropriate spectral range was selected by comparing the modeling results of different spectral intervals. The results are shown in Table 4.Table 4
Experimental results of spectral range screening.
Spectral rangeR2RMSECRMSEPRPD4200–5380 cm−10.53843.1683.4881.475390–6060 cm−10.47323.3843.7141.386070–7500 cm−10.82551.9482.0812.394000–10000 cm−10.90921.4051.5073.32In Table4, when the spectral range was 4000–10000 cm−1, the values of RMSEC and RMSEP were the smallest, and the model effect was the best.


### 3.5. Selection of Spectral Pretreatment Methods

It is necessary to preprocess the original spectrogram. The robustness and accuracy of the model can be improved by preprocessing the spectrogram. Common preprocessing methods include S-G convolution smoothing, FD (first derivative), SD (second derivative), SNV (standard normal variable transformation), and MSC (multivariate scattering correction). In this study, R2, RMSEC, and RMSEP were used as model evaluation indicators to compare the modeling results of different pretreatment methods and selected the appropriate spectral pretreatment method.In Table5, the best pretreatment method of total saccharide model of A. auricula was “SNV + SD.” The spectrum of A. auricula after pretreatment was shown in Figure 4.Table 5
Experimental results of spectral pretreatment.
Pretreatment methodR2RMSECRMSEPRPDMSC + FD0.74102.3732.6021.96MSC + SD0.90831.4121.5143.30SNV + FD0.77882.1932.4152.12SNV + SD0.90921.4051.5073.32Figure 4
Spectrum ofA. auricula after pretreatment.


### 3.6. Determination of the Main Factor of PLS

When the method of partial least squares was used to establish the model, the appropriate number of PLS main factors was related to the quality of the model directly. If the number of PLS main factors was too small, it would cause the loss of useful information in the original spectrum, which would often reduce the predictive ability of the model. If the number of PLS main factors was too much, some information that was not related to the sample would be included. This situation would lead to a larger prediction error of the model [25]. According to the comprehensive evaluation index R2 and RMSEP, the number of main factors used in PLS modeling was 11.


## 3.1. Determination Results of Total Saccharide Content of A. auricula

According to formula (1), the measurement results of the total saccharide content in 166 batches were shown in Table 2.Table 2
The results of total saccharide content ofA. auricula.
NumberTotal saccharide content (%) (x ± s, n = 3)114.54 ± 0.37215.28 ± 0.12312.03 ± 1.45412.60 ± 1.20518.37 ± 3.31613.60 ± 2.86710.79 ± 1.25811.07 ± 0.97915.42 ± 3.021018.60 ± 1.161113.14 ± 2.07129.26 ± 0.621311.82 ± 1.361411.37 ± 0.661511.92 ± 1.281612.23 ± 3.561719.78 ± 0.981812.47 ± 2.201918.22 ± 2.032019.36 ± 0.802114.43 ± 0.262217.86 ± 2.132317.82 ± 3.902413.25 ± 1.542514.04 ± 0.27264.52 ± 0.57277.91 ± 3.66286.08 ± 1.00294.01 ± 0.08305.94 ± 1.99319.56 ± 0.44326.31 ± 3.79339.62 ± 1.313411.49 ± 2.603512.33 ± 1.343615.61 ± 0.973712.86 ± 3.313811.91 ± 1.393912.77 ± 2.11406.99 ± 0.014114.91 ± 1.00423.83 ± 1.57439.02 ± 0.47445.41 ± 0.06455.54 ± 0.25464.98 ± 1.04473.33 ± 2.284810.07 ± 1.89494.62 ± 0.135010.67 ± 0.645111.57 ± 2.88525.75 ± 1.31534.63 ± 0.145410.01 ± 2.13552.85 ± 1.50563.59 ± 1.36572.99 ± 1.135811.85 ± 0.41595.25 ± 1.536011.60 ± 1.126120.93 ± 0.416219.19 ± 3.296319.47 ± 3.696414.07 ± 2.52658.33 ± 0.66663.48 ± 1.02677.09 ± 0.76684.78 ± 1.80696.27 ± 0.267014.66 ± 2.937125.24 ± 0.32728.51 ± 1.53735.16 ± 0.197414.52 ± 2.937512.22 ± 0.18765.11 ± 1.917713.56 ± 0.77810.54 ± 0.27797.35 ± 0.06809.38 ± 1.2816.51 ± 0.05828.90 ± 0.468321.48 ± 0.68414.83 ± 1.1854.05 ± 1.488614.43 ± 0.58715.68 ± 0.75889.98 ± 2.36896.65 ± 1.93905.85 ± 1.47915.16 ± 0.049212.29 ± 0.96937.59 ± 0.97948.41 ± 0.22956.69 ± 0.369614.54 ± 1.07979.51 ± 0.73981.30 ± 0.05995.55 ± 0.341007.88 ± 0.521015.73 ± 1.310217.33 ± 1.4810311.61 ± 0.081044.09 ± 0.9510510.52 ± 0.1910612.61 ± 3.131078.88 ± 0.3710818.72 ± 0.0410911.02 ± 1.641109.66 ± 1.7311112.62 ± 0.581122.61 ± 0.6411310.24 ± 0.311143.27 ± 0.231157.62 ± 0.1511612.23 ± 1.701172.67 ± 0.171188.75 ± 3.0711914.86 ± 0.481206.10 ± 0.2912111.01 ± 1.7312213.83 ± 3.3112311.26 ± 2.651245.97 ± 2.5312511.53 ± 2.2012617.25 ± 3.111276.89 ± 0.151284.11 ± 1.281298.19 ± 2.311309.59 ± 2.0213115.18 ± 1.321323.81 ± 2.5113314.84 ± 0.8013414.39 ± 0.4613514.51 ± 1.451363.66 ± 0.1713712.90 ± 0.791382.56 ± 0.151396.38 ± 0.271404.86 ± 0.1114110.27 ± 0.341429.79 ± 0.6714315.98 ± 1.8614414.39 ± 0.811452.11 ± 1.5314614.10 ± 0.214710.61 ± 0.321486.68 ± 1.6814912.00 ± 1.5615016.17 ± 1.291516.48 ± 0.1615214.46 ± 0.6415311.69 ± 1.711549.96 ± 3.0315512.10 ± 0.4315613.69 ± 0.7415710.38 ± 3.2315813.77 ± 0.7815912.01 ± 0.4116011.95 ± 0.5016113.14 ± 0.1816213.85 ± 0.701638.44 ± 2.141644.61 ± 0.361657.04 ± 0.6416615.71 ± 0.09


## 3.2. Removal of Model Abnormal Samples

In this study, software was used to identify the abnormal samples automatically. The PCA score map was used to cluster analyse the spectral data ofA. auricula samples, and the differences among the samples were counted. Diagnosis of abnormal data by test of studentized residual and leverage values method [23]. Finally, 47 batches of abnormal samples were removed from the model. The modeling process was shown in Figures 2 and 3.Figure 2
Test of studentized residual.Figure 3
Leverage values of samples.


## 3.3. Division of Calibration Set and Verification Set

The calibration set and verification set were selected from all samples randomly. The following two conditions must be met during the selection process: ① ensure that the total saccharide content of the verification set samples covers the content gradient of all samples. ② Ensure that the content range of the selected verification set samples does not exceed the calibration set.In Table3, the samples of the verification set were completely distributed within the range of the calibration set, indicating that the division of the sample set meets the modeling requirements.Table 3
Range of total saccharide content in calibration samples and validation samples.
Sample setNumber of samplesMaximum (%)Minimum (%)Correction set10419.18962.5641Validation set1517.33003.5860


## 3.4. Selection of Spectral Range for Modeling

The calibration model is established by choosing appropriate spectral range which has strong prediction ability, and good robustness can often be obtained [24]. In Figure 1, the sample had obvious characteristic peak absorption at 5100 cm−1 and 7000 cm−1. Meanwhile, R2, RMSEC, RMSEP, and RPD were used as model evaluation indicators to establish the model in four spectral intervals. The appropriate spectral range was selected by comparing the modeling results of different spectral intervals. The results are shown in Table 4.Table 4
Experimental results of spectral range screening.
Spectral rangeR2RMSECRMSEPRPD4200–5380 cm−10.53843.1683.4881.475390–6060 cm−10.47323.3843.7141.386070–7500 cm−10.82551.9482.0812.394000–10000 cm−10.90921.4051.5073.32In Table4, when the spectral range was 4000–10000 cm−1, the values of RMSEC and RMSEP were the smallest, and the model effect was the best.


## 3.5. Selection of Spectral Pretreatment Methods

It is necessary to preprocess the original spectrogram. The robustness and accuracy of the model can be improved by preprocessing the spectrogram. Common preprocessing methods include S-G convolution smoothing, FD (first derivative), SD (second derivative), SNV (standard normal variable transformation), and MSC (multivariate scattering correction). In this study, R2, RMSEC, and RMSEP were used as model evaluation indicators to compare the modeling results of different pretreatment methods and selected the appropriate spectral pretreatment method.In Table5, the best pretreatment method of total saccharide model of A. auricula was “SNV + SD.” The spectrum of A. auricula after pretreatment was shown in Figure 4.Table 5
Experimental results of spectral pretreatment.
Pretreatment methodR2RMSECRMSEPRPDMSC + FD0.74102.3732.6021.96MSC + SD0.90831.4121.5143.30SNV + FD0.77882.1932.4152.12SNV + SD0.90921.4051.5073.32Figure 4
Spectrum ofA. auricula after pretreatment.


## 3.6. Determination of the Main Factor of PLS

When the method of partial least squares was used to establish the model, the appropriate number of PLS main factors was related to the quality of the model directly. If the number of PLS main factors was too small, it would cause the loss of useful information in the original spectrum, which would often reduce the predictive ability of the model. If the number of PLS main factors was too much, some information that was not related to the sample would be included. This situation would lead to a larger prediction error of the model [25]. According to the comprehensive evaluation index R2 and RMSEP, the number of main factors used in PLS modeling was 11.


## 4. Verification and Evaluation of the Model

### 4.1. External Verification

The near-infrared quantitative model was used to predict the total saccharide content of the validation set samples, and the predicted value of the quantitative model and the measured chemical value were shown in Table6 and Figure 5. In Figure 5, R2 = 0.9048 was higher than 0.9 of the verification set, which indicated that there was a good linear relationship between the predicted value and the measured value of the model. Then, the statistical software of SPSS was used to test the predicted results with paired T-test. The set confidence was 95%. The results showed that there was no significant difference between the predicted value and the measured chemical value of the quantitative model (P = 0.182 ≥ 0.05).Table 6
Prediction results of total saccharide content of validation samples.
NumberChemical value (%)Predicted value (%)3512.3311.86563.594.286011.6012.35677.096.79728.516.617713.5613.29948.415.8810217.3318.271309.598.201429.7911.3415214.4613.2315710.3810.971644.614.9516615.7113.458715.6814.27Figure 5
The predicted vs. chemical value of validation samples.


### 4.2. Evaluation of the Model

The parameters of the quantitative model were shown in Table7, and the linear scatter diagram of the chemical value and the predicted value of the model was shown in Figure 6. The best near-infrared model of total saccharide content of A. auricula was obtained by the screening of modeling parameters.Table 7
Parameters of the optimal NIR models for total saccharide.
ProjectsNIR model parametersPretreatment methodSNV + SDSpectral range (cm−1)4000–10000PLS main factor number11R20.9092RMSEC1.405RMSEP1.507RPD3.32Figure 6
The predicted vs. measured value of the optimal model for total saccharide.


## 4.1. External Verification

The near-infrared quantitative model was used to predict the total saccharide content of the validation set samples, and the predicted value of the quantitative model and the measured chemical value were shown in Table6 and Figure 5. In Figure 5, R2 = 0.9048 was higher than 0.9 of the verification set, which indicated that there was a good linear relationship between the predicted value and the measured value of the model. Then, the statistical software of SPSS was used to test the predicted results with paired T-test. The set confidence was 95%. The results showed that there was no significant difference between the predicted value and the measured chemical value of the quantitative model (P = 0.182 ≥ 0.05).Table 6
Prediction results of total saccharide content of validation samples.
NumberChemical value (%)Predicted value (%)3512.3311.86563.594.286011.6012.35677.096.79728.516.617713.5613.29948.415.8810217.3318.271309.598.201429.7911.3415214.4613.2315710.3810.971644.614.9516615.7113.458715.6814.27Figure 5
The predicted vs. chemical value of validation samples.


## 4.2. Evaluation of the Model

The parameters of the quantitative model were shown in Table7, and the linear scatter diagram of the chemical value and the predicted value of the model was shown in Figure 6. The best near-infrared model of total saccharide content of A. auricula was obtained by the screening of modeling parameters.Table 7
Parameters of the optimal NIR models for total saccharide.
ProjectsNIR model parametersPretreatment methodSNV + SDSpectral range (cm−1)4000–10000PLS main factor number11R20.9092RMSEC1.405RMSEP1.507RPD3.32Figure 6
The predicted vs. measured value of the optimal model for total saccharide.


## 5. Discussion

The quantitative analysis of edible mushroom mainly focused on the determination of water, protein, polysaccharides, triterpene, and other components. Although the chemical standard analysis method can detect the total saccharide content, its operation steps are tedious and the efficiency is low; there are some limitations in the practical application process. In the research on the detection of effective components inCordyceps militaris, Wang Di [26] proved that it was feasible to use near-infrared technology combined with PLS algorithm to achieve rapid prediction of effective components in edible mushroom. Li Junshan [27] directly measured the water content in Poria cocos by near-infrared technology combined with PLS method and established a prediction model. The optimal number of the latent variables for the model was 8. The coefficient of determination (R2) for the model was 0.998. The RMSECV was 0.453, and RMSEP was 0.366. The results of the model showed that the method of determining the water content of Poria cocos by near-infrared technology was faster and simpler than other methods. Wu Lun [28] determined the saponins content of the Honey-Fried Processing of Rhizoma Cimicifugae by near-infrared technology. Through the study of the band range, the spectrum used for modeling was finally determined, and the intervals are 5200–6700 cm−1 and 7700–8800 cm−1. This study showed that the selection of spectral range played an important role in the prediction effect of the model. Li Xiaowen [29] determined the total saccharide content of Coprinus comatus by near-infrared technology. 80 samples of C. comatus with different maturity in the same growth cycle were collected, and the PLS model was established after pretreatment such as normalization, second derivative, and smoothing. The results showed that the correlation coefficients of calibration set and prediction set were 0.9989 and 0.9823, respectively. The total saccharide content of C. comatus was detected quickly and accurately.According to the documentation, the combination of near-infrared spectroscopy and chemometrics to establish a prediction model can quickly and accurately predict the content of the tested substance, and it is suitable for the rapid detection of nutrients in edible mushroom.


## 6. Conclusion

In this study, the near-infrared model of the total saccharide ofA. auricula established had 104 calibration sets and 15 validation sets, the best pretreatment method was SNV +SD, the best spectral range was 4000–10000 cm−1, and the best PLS main factor number was 11. R2, RMSEC, RMSEP, and RPD of the quantitative model were 0.9092, 1.405, 1.507, and 3.32, respectively. The predicted value of the model had a good linear relationship with the measured value. According to the T-test of the validation set, the difference between the predicted value and the measured chemical value was not significant (P ≥ 0.05).The results showed that the quantitative model had good fitting ability, and it is feasible to determine the total saccharide content ofA. auricula by near-infrared spectroscopy. This method is simple and accurate, which can effectively determine the total saccharide content of A. auricula. It can provide technical reference for further expanding the application of near-infrared spectroscopy in the determination of total saccharide content, and it is expected to be used in the detection and analysis of the quality of traditional Chinese medicine.



---

*Source: 8858235-2022-02-24.xml*