Irsa Ejaz1
Siyang He1†
Wei Li1†Naiyue Hu1†
Chaochen Tang1,2
Songbo Li1Meng Li1,3,4Boubacar Diallo1,5,6
Guanghui Xie1,6
Kang Yu7*- 1College of Agronomy and Biotechnology, China Agricultural University, Beijing, China
- 2Crops Research Institute, Guangdong Academy of Agricultural Sciences, Guangzhou, China
- 3College of Bioscience and Biotechnology, Hunan Agricultural University, Changsha, China
- 4Hunan Branch, National Energy R&D Center for Non-food Biomass, Hunan Agricultural University, Changsha, China
- 5Department of Agriculture, High Institute Agronomic and Veterinary, Faranah, Guinea
- 6National Energy R&D Center for Non-food Biomass, China Agricultural University, Beijing, China
- 7School of Life Sciences, Technical University of Munich, Freising, Germany
Near-infrared spectroscopy (NIR) is a non-destructive, fast, and low-cost method to measure the grain quality of different cereals. However, the feasibility for determining the critical biochemicals, related to the classifications for food, feed, and fuel products are not adequately investigated. Fourier-transform (FT) NIR was applied in this study to determine the eight biochemicals in four types of sorghum samples: hulled grain flours, hull-less grain flours, whole grains, and grain flours. A total of 20 hybrids of sorghum grains were selected from the two locations in China. Followed by FT-NIR spectral and wet-chemically measured biochemical data, partial least squares regression (PLSR) was used to construct the prediction models. The results showed that sorghum grain morphology and sample format affected the prediction of biochemicals. Using NIR data of grain flours generally improved the prediction compared with the use of NIR data of whole grains. In addition, using the spectra of whole grains enabled comparable predictions, which are recommended when a non-destructive and rapid analysis is required. Compared with the hulled grain flours, hull-less grain flours allowed for improved predictions for tannin, cellulose, and hemicellulose using NIR data. This study aimed to provide a reference for the evaluation of sorghum grain biochemicals for food, feed, and fuel without destruction and complex chemical analysis.
Introduction
Sorghum (Sorghum bicolor L.) is the fifth most commanding cereal crop in the world and its grain production has reached up to 57.50 million tons in 2020 (de Morais Cardoso et al., 2017; Food and Agriculture Organization of the United Nations, 2020; Stamenković et al., 2020). Sorghum grain is known for its nutritional quality and there is a worldwide growing market of sorghum grain for use as human food and consumed in the preparation of many foods (Ratnavathi and Patil, 2014; Bader Ul Ain et al., 2019; Sihono et al., 2019; Palacios et al., 2021). It is an excellent feed for animals and its feeding value is generally considered more than 95% of the feeding value of yellow dent maize (Waniska et al., 2016). In addition, it is a unique energy crop that can be used through numerous production routes—for instance, converting starch into ethanol, lignocelluloses into bio-oil, biochar, and biohydrogen, and fat into biodiesel (Mirfakhar et al., 2020; Stamenković et al., 2020). Sorghum cultivars can be classified into four classes: sorghum, white sorghum, tannin sorghum, and mixed sorghum according to the US Grain Standards (Waniska et al., 2016). The grains from different cultivars differ in the presence or absence of thick outer-covering/pericarp called hulled or hull-less grains (Waniska and Rooney, 2000; Guindo et al., 2016). Sorghum for food must have lower fat but higher starch and protein for good digestibility, whereas for feed, it must have higher protein content but lower concentrations of tannin and insoluble dietary fibers especially hemicellulose (Waniska et al., 2016). The insoluble dietary fibers are desirable by many biofuel industries due to its property of making a low viscous solution (Fuller et al., 2016; Qiu et al., 2017; Stamenković et al.,2020). Therefore, a precise classification of sorghum grains requires the determination of various chemical components.
Sorghum grain hulls comprise of condensed tannins that have cross-linkages with starch, protein, and lignocellulosic components (Evers and Millar, 2002; Siwela et al., 2007; Barros et al., 2012; Gyori, 2017). Cellulose, hemicellulose, and lignin are important soluble and insoluble dietary fibers that have complex structural polymerization and are difficult to be determined directly. Unlike other cereal grains, cellulose, hemicellulose, and lignin in sorghum grains are washed out with other sugars during their extraction process (Jin et al., 2017; Li et al., 2017; Krasznai et al., 2018). To cope with this problem, acid hydrolysis has to be done several times, which, however, is not desirable due to specialized laboratory facilities (Ricardo Soccol et al., 2011; Heredia-Olea et al., 2012). In addition, these wet-chemical methods are strenuous, time consuming, and expensive for large samples, making it infeasible for real-time sorting of grains (Orman and Schumann, 1991; Guindo et al., 2016; Yang et al., 2016; Caporaso et al., 2018; Srivastava et al., 2018). Therefore, inexpensive and immediate quantification methods are needed to sort out sorghum grains for food, feed, and fuel.
Near-infrared (NIR) spectroscopy is a non-destructive, rapid, and low-cost method to estimate the biochemicals of cereal grains and differentiate them based on their chemical composition (Font et al., 2006; Hell et al., 2016; Caporaso et al., 2018). It is a feasible alternative to the time-inefficient and resource-intensive conventional methods of analysis, such as the Kjeldahl or high-performance liquid chromatography (HPLC) (Beć et al., 2021). A recent trend is to miniaturize NIR instruments that reduce the weight and cost (Wiedemair et al., 2019). However, Benchtop Fourier-Transform (FT) NIRS operates over the entire wavelength region with a high spectral resolution, and a good signal-to-noise ratio for yielding fast and accurate measurements (Beć et al., 2021). The NIR-based prediction models have been developed to estimate the chemical composition in grains of various crops (Pohl and Senn, 2011; Ferreira et al., 2015; Li et al., 2015; Kamboj et al., 2017). Previously reported NIR model calibrations using the whole grain spectra produce different results from using flour spectra and had moderate accuracy (De Alencar Figueiredo et al., 2006). Improved predictions are achievable by using the spectra of whole grains, when compared with the use of spectra collected from the flour samples (Beloshapka et al., 2016; Hu et al., 2021). It is also reported that the influence of grain hulls on the prediction of grain chemical content varies for different cultivars, and that milled grain flours might improve the prediction for grain chemical content (Wiedemair et al., 2019). However, due to the variation of hulls presence or absence in sorghum grains, the influence of hulls on the spectroscopic determination of cellulose, hemicellulose, lignin, and other components in sorghum grains is still unclear. The previous studies on how to effectively estimate chemical composition in different sorghum grain types are insufficient.
Thus, the objectives of this research were: (1) to evaluate the feasibility of using FT-NIR spectroscopy to determine a variety of chemical components, especially cellulose, hemicellulose, and lignin in sorghum grains; and (2) to evaluate the influence of four sample types (whole grains, flours, hulled grain flours, and hull-less grain flours) on the prediction of chemical components to improve the grain sorting efficiency for human food, animal feed, and biofuel.
Materials and Methods
Sampling and Measurements
A total of 20 hybrids of sorghum were grown at the experimental stations of China Agricultural University in Zhuozhou (Hebei province) and Jiexiu (Shanxi province), China, in 2017. In total 98 samples were collected from both the experimental stations (Supplementary Table 1). All the samples were first oven-dried at 45°C for 10 days to constant weight and then scanned as whole grains to acquire absorbance spectra using Thermo Antaris II FT-NIR spectrometer (Thermo Scientific Inc., Madison, WI, USA). Subsequently, the sorghum grain samples were milled to flours using a grinder, sifted through a 40–80 mesh screen, and scanned again as hulled or hull-less grain flours. The flours were stored at 4°C for chemical analysis. Starch, protein, fat, tannin, cellulose, hemicellulose, lignin, and ash contents were determined using these flour samples (Table 1).
Table 1. Descriptive statistics [minimum (min), maximum (max), mean, SD, and coefficient of variation (CV)] of laboratory analytical data comprising starch, protein, fat, tannin, cellulose, hemicellulose, lignin, and ash contents of distinct types of sorghum grains (whole grains, whole grain flours, hulled grain flours, and hull-less grain flours), used as reference data for calibration (C) and validation (V) subsets to construct PLSR models.
Determination of Chemical Components and Sample Classification
Starch content in sorghum grains was determined by the Anthrone colorimetric method and assayed using a UV-visible (UV-VIS) spectrometer (TU-1901, Beijing Purkinje Instruments Co., Ltd., Beijing, China) (Li et al., 2014). Crude protein was determined by Kjeldahl method through the measurement of total nitrogen. The content of crude fat was determined by ether extraction (Padmore, 1990). Tannin was extracted with dimethyl amide solution and tannic acid was used as a standard to determine the content of sorghum grain tannin (ISO, 1988). Cellulose, hemicellulose, and lignin were determined according to the modified National Renewable Energy Laboratory standards (NREL/TP-510-42618, Revised August 2012). Acid hydrolysis was done using the sulfuric acid (H2SO4). Acid insoluble lignin in the acid hydrolysis solution was determined to calculate the content of cellulose, hemicellulose, and lignin according to the following Equations (1–3), respectively.
Additionally, we used perchloric acid as sorghum grains had a greater amount of starch and lignocellulosic components that were difficult to determine without using perchloric acid for acidification. After acidification, NREL standard methods were followed to determine cellulose, hemicellulose, and lignin content in sorghum grains (Sluiter et al., 2012). The ash content was determined by using a muffle furnace (VULCAN 3-550, Dentsply International Inc., York, PA, USA).
The samples were classified into three categories as food, feed, and fuel based on the available chemical content in the grains. The samples (n = 98) with equal or more than 15 g kg−1 (1.5%) tannin and 50 g kg−1 (5%) hemicellulose content were categorized as fuel samples (Waniska et al., 2016) (Table 1). The food and feed samples contained tannin and hemicellulose content < 15 and 50 g kg−1, respectively. The samples having starch < 650 g kg−1 (65%) were classified as feed and ≥ 650 g kg−1 (65%) were classified as food samples (Pontieri et al., 2010; Waniska et al., 2016; Beloshapka et al., 2016). Following the chemical properties, 14 samples were classified for food, 21 samples for feed, 38 samples for fuel, and 25 samples classified for both feed and fuel. The sample classification is shown in Supplementary Table 1.
FT-NIR Spectroscopy and Data Processing
The FT-NIR spectra were scanned using the Thermo Antaris II FT-NIR spectrometer (Thermo Scientific Inc., Madison, WI, USA) prepared with an absorbance accessory. The grain and flour samples were put in a special Petri dish (black sides and transparent glass bottom) and fitted in the FT-NIR machine. Technical Quality (TQ) Analyst-pro software (Thermo Scientific Inc., Madison, WI, USA) was used to run the machine. Each spectrum had 64 scans at 4 cm−1 resolutions with the wavenumber range among 4,000–10,000 cm−1, including 1,557 spectral variables.
All obtained absorption spectra were analyzed using ChemDataSolution Version 3.1.0 (Dalian Chem. Data Solution Technology Co. Ltd., Dalian, China) (Figure 1). The raw spectra with the whole spectral range were preprocessed to remove random noise and small peaks by the combination of pretreatments: Multiplicative scatter correction (MSC) (Geladi and Macdougall, 1985) and standard normal variate (SNV) transformation (Brown et al., 2000) were used to remove the variations caused by instrument settings, sample, and environmental conditions. Norris derivative by means of Savitzky Golay (SG) method with 52 smoothing points (Brown and Wentzell, 2000) were used to resolve the spectra peak overlap and evacuate linear baseline drift with the selection of whole wavelength region (4,000–10,000 cm−1). The purpose of the pretreatments was to remove multiplicative and additive effects due to the machine settings and variations caused by the sample and environmental conditions. Based on different pretreatment combinations and performance, the models were cross-validated and the partial least square (PLS) factors were identified in the models for cross-validation (Geladi and Macdougall, 1985) (Table 3). To select the optimal number of PLS factors, the “leave-one-out” cross-validation method was used for developing PLS models to avoid over-fitting (Vehtari et al., 2017). Following the preprocessing treatments, the sample choice for the calibration and validation was carried out by the Kennard–Stone (KS) algorithm and Sample-set partitioning based on joint X–Y distances (SPXY) with the ratio of 3:1 (calibration:validation), respectively (Kennard and Stone, 1969) (Table 2). While KS algorithm only concerns the similarity between independent variables of the two subsets, SPXY combines independent and dependent variables (Galvão et al., 2005). All the spectral data were mean-centered before calibration and the competitive adaptive reweighted sampling (CARS) method was used to select the wavelength (Li et al., 2009; Wang et al., 2015). The regression coefficients of PLS models were used as an index for evaluating the importance of each wavelength (Li et al., 2009). The PLS regression was used to construct the calibration models for the biochemical evaluation of grains (Figure 1).
Figure 1. The morphological classification of sorghum grains by means of hulled and hull-less grains, also used as whole grains and flours by means of Fourier Transform near-infrared spectroscopy (FT-NIR) spectral data collection followed by reference analysis. The spectra were processed by combination of pre-treatments, cross-validated, calculation of principal components (PCs) by principal component analysis (PCA), divided into calibration and validation subsets and ultimately, partial least square regression (PLSR) models were built using FT-NIR spectral and reference data.
Table 2. The PLSR followed by mathematical treatments of spectra to remove random noise and to construct FT-NIR models for biochemical components (starch, protein, fat, tannin, cellulose, hemicellulose, lignin, and ash) of whole grains, whole grain flours, hulled, and hull-less grain flours of sorghum.
First, all the values of each biochemical were built-in the model to gain an understanding of the initial standard error and their correlation coefficients. Outliers were detected from the first calibrated model by calculating the CI (95%) of the first two principal components (PCs) from principal component analysis (PCA) plots in Moli Software (Dalian Chem. Data Solution Technology Co. Ltd., Dalian, China) and the samples out of this interval were removed individually from the sample set of each relevant model. The models were recalibrated by following the above spectra preprocessing methods to lower the number of PCs.
Performance Index of Partial Least Squares Regression Models
The performance of multivariate calibrations was evaluated according to the determination coefficient (R2), root mean square error of calibration (RMSEc), root mean square error of validation (RMSEv), and the ratio of prediction to deviation (RPD) for the calibration and validation subsets. These parameters were used to determine the potential of the partial least squares regression (PLSR) models for the determination of biochemical components in sorghum grains (Wu et al., 2015; Yang et al., 2016) (Table 3). In model development and applications, R2 > 0.95 and RPD > 4 indicates that the model could predict efficiently, 0.95 > R2 > 0.9 and 4 > RPD > 3 regarded as successful calibration, 0.9 > R2 > 0.8 and 3 > RPD > 2.25 considered satisfactory, and 0.8 > R2 > 0.7 and 2.25 > RPD > 1.75 is regarded good for preliminary screenings (Malley et al., 2004; Arana et al., 2005) (Table 3).
Table 3. Statistics of PLSR models constructed using FT-NIR spectroscopy to determine starch, protein, fat, tannin, cellulose, hemicellulose, lignin, and ash in distinct types of sorghum grains (whole grains, whole grain flours, hulled, and hull-less grain flours) used for calibration, validation, and cross-validation subsets.
Results
Chemical Properties and Grain Grading
The minimum contents of tannin and high starch in food grains were 11.66 and 695.57 g kg−1, respectively found in 622A × J7645Z hybrid (Table 1; Supplementary Table 1). The minimum hemicellulose content in the food samples was 29.92 g kg−1 in 622A × J7645Z (Table 1; Supplementary Table 1). The Awanlek hybrid contains minimum tannin and starch content of 6.44 and 597.11 g kg−1, respectively, useful for feed (Table 1; Supplementary Table 1). Fuel samples contain maximum cellulose and hemicellulose content (125.86 and 80.35 g kg−1, respectively) (X098 and NW3) (Table 1; Supplementary Table 1). J7645Z hybrid contains maximum starch content that is useful for fuel samples (730.68 g kg−1) (Table 1; Supplementary Table 1). The maximum content of protein, fat, lignin, and ash in fuel samples were 169.44, 44.91, 75.99, and 30.59 g kg−1, respectively. The samples used for both feed and fuel, contain low tannin and starch (useful for feed) but high cellulose, hemicellulose, and lignin content (useful for fuel). The minimum starch and tannin content in combined feed/fuel samples were 498.65 and 13.85 g kg−1, respectively (AMP450 × NW3) (Table 1; Supplementary Table 1). Moreover, the average chemical properties for different morphological types of sorghum grains revealed that hull-less grain flours had high levels of starch content (647.4 g kg−1) compared with the whole grains/flours (585.6 g kg−1) and hulled grain flours (548.1 g kg−1) (Table 1; Supplementary Table 1). By contrast, the content of protein and fat were not significantly different for all types of datasets. Tannin content was found higher in hulled flours (15.0 g kg−1) and lower in hull-less flours (11.1 g kg−1) compared to whole grains/flours samples (13.6 g kg−1) (Table 1). The average cellulose content was the highest in hull-less flours (124.3 g kg−1) and the least in hulled flours (116.6 g kg−1) while on the contrary, average hemicellulose content was the highest in hulled flours (66.4 g kg−1) and the least in hull-less flours (35.3 g kg−1) compared with whole grains/flours. Similarly, the average lignin and ash content were also high in hulled flours (68.2 and 26.3 g kg−1, respectively) (Table 1).
The coefficient of variation (CV) for the biochemical components in different sorghum morphological types (grains, flours, hulled flours, and hull-less flours) used as reference data for calibration and validation subsets are shown in Table 1. Biochemical data used in the whole grains and flours model showed that hemicellulose and lignin have the highest CVs for both the calibrations (30.52 and 28.48%, respectively) and validation subsets (21.42 and 19.12%, respectively). Tannin has CV for both calibration and validation (30.27 and 9.08%, respectively), followed by cellulose (26.06 and 12.66%, respectively), fat (24.32 and 15.07%, respectively), ash (18.94 and 14.0%, respectively), starch (16.67 and 12.95%, respectively), and the smallest by protein (10.89 and 9.66%, respectively). The biochemical data used in the hulled and hull-less grain flour models showed similar ranges of CV (Table 1).
Spectral Pretreatment and Models Performance
The mean, minimum, and maximum absorbance of sorghum grain samples are shown in Figures 2, 3. The minimum absorption was found between 10,000 and 9,000 cm−1 wavelength regions. The maximum absorption peaks were found in the 9,000–4,000 cm−1 wavelength regions and showed a common phenomenon in whole grains, whole grain flours, hulled, and hull-less flours. The absorbance of whole-grain flours was the least among all types of samples. Whole grain flours absorbed at around 10,000 cm−1, with values ranging from 0.25 to 0.4 (Figure 2) that was adjacent to the absorption of hull-less and hulled grain flour which were absorbed between 0.4 and 0.43, respectively (Figure 3).
Figure 2. Summary of spectral absorbance and characteristic wavelengths distribution measurements for sorghum whole grains vs. flours sample set (n = 98). The range of coefficient variation (CV, blue line), mean absorbance (black line), minimum/maximum (green dotted line) spectral absorbance, total range of absorbance measurements (red line) for each wavelength (cm−1), and SD+/– SD between sample sets were calculated.
Figure 3. Summary of spectral absorbance and characteristic wavelengths distribution measurements for sorghum hulled (n = 61) vs. hull-less grain flours sample set (n = 37). The CV, coefficient of variation (blue line) (SD/mean), mean absorbance (black line), minimum/maximum (green dotted line) spectral absorbance, total range of absorbance measurements (red line) for each wavelength (cm−1), and +/– SD between sample sets were calculated.
The development of PLS regression models for each biochemical component of sorghum grains were performed to minimize the effect of the change in the baseline drift. Most of the models performed best in a combination of second derivative SG52, CARS, and SPXY methods, except for lignin in whole grains and starch, protein, tannin, and ash in hull-less grain flours that performed best with first derivative and KS sample selection method. All the models were constructed by critically evaluating different spectral pretreatments that gave the higher R2c (≥0.96), R2v (≥0.72), R2cv (≥0.84), RPDv (≥1.9), and RPDcv (≥2.5) (Table 3). A T-test was applied to determine the prediction accuracy of the models for different sample formats (Supplementary Figures 1, 2).
Comparison of the Models of Whole Grains and Flours
Table 3 shows model statistics for all biochemical components in the models of whole grains and flours. In the whole grain models, the R2c for all the chemical components were ≥0.97, R2v were ≥0.72 and R2cv were ≥0.90, while RPDc were ≥1.94 and RPDcv were ≥3.18. Based on the highest R2c (0.99) and RPDc (6.18), hemicellulose model was considered best among the other chemical components in whole grain models (Table 3). In flours models, the R2c for all the chemical components were ≥0.96, R2v were ≥0.90, and R2cv were ≥0.84, while RPDc were ≥3.05 and RPDcv were ≥2.54. Hemicellulose and lignin were considered best among the other chemical components in flours models based on the highest R2c (0.99 for both) and RPDc (6.89 and 6.92, respectively) (Table 3). Overall, the model validation for flours was better than the whole grains. Figures 4, 5 show the scatter plots for predicted and measured values of different biochemical components in the whole grains and flours, respectively.
Figure 4. Comparison of scatter plots of measured vs. predicted values (g kg−1) of the whole grains for starch, protein, fat, tannin, cellulose, hemicellulose, lignin, and ash of sorghum grains for the external validation subsets based on PLSR models. The black and gray lines in plot showed mean, maximum, and minimum values while dotted lines showed average chemical composition between calibration and validation sample set.
Figure 5. Comparison of scatter plots of measured vs. predicted values (g kg−1) of the whole grain flours for starch, protein, fat, tannin, cellulose, hemicellulose, lignin, and ash of sorghum grains for the external validation subsets based on PLSR models. The black and gray lines in plot showed mean, maximum, and minimum values while dotted lines showed average chemical composition between calibration and validation sample set.
Comparison of the Models of Hulled and Hull-Less Flours
The statistics for all the biochemical components in the models of hulled and hull-less grain flours are shown in Table 3. In hulled grain flours models, the R2c for all the chemical components were ≥0.98, R2v were ≥0.90, and R2cv were ≥0.92, while RPDc were ≥2.74 and RPDcv were ≥3.48. Based on the highest R2c (0.99) and RPDc (7.17), protein model was considered the best among other chemical components in hulled grain flours models (Table 3). In hull-less grain flours models, the R2c for all the chemical components were ≥0.97, R2v were ≥0.93 and R2cv were ≥0.88, while RPDc were ≥2.88 and RPDcv were ≥ 3.01. Fat, tannin, and cellulose were considered the best among other chemical components in flours models based on the highest R2c (0.99, 0.97, and 0.99, respectively) and RPDc (6.28, 8.46, and 6.1, respectively) (Table 3). Figures 6, 7 show the scatter plots for predicted and measured values of different biochemical components in the hulled and hull-less and flours, respectively. The hulled grain flours allowed for the best prediction for protein, whereas the hull-less grain flours allowed for the best prediction for tannin, cellulose, and hemicellulose (Figures 6, 7).
Figure 6. Comparison of scatter plots of measured vs. predicted values (g kg−1) of the hulled grain flours for starch, protein, fat, tannin, cellulose, hemicellulose, lignin, and ash of sorghum grains for the external validation subsets based on PLSR models. The black and gray lines in plot showed mean, maximum, and minimum values while dotted lines showed average chemical composition between calibration and validation sample set.
Figure 7. Comparison of scatter plots of measured vs. predicted values (g kg−1) of the hull-less grain flours for starch, protein, fat, tannin, cellulose, hemicellulose, lignin, and ash of sorghum grains for the external validation subsets based on PLSR models. The black and gray lines in plot showed mean, maximum, and minimum values while dotted lines showed average chemical composition between calibration and validation sample set.
Discussion
Chemical Properties and Food, Feed, and Fuel Grading
In this study, sorghum grains were graded based on their respective chemical content for food, feed, and fuel. Out of the 20 hybrids of sorghum grains used in this study, hull-less hybrids which include 622A × Awanlek, 622A × J7645Z, and 624A × J7645Z contained high starch and cellulose while low tannin and hemicellulose content and were considered suitable for food. The protein content was high in the food samples compared with the feed samples. Hull-less hybrids (Awanlek, 624A × Awanlek, AMP450 × Awanlek, and AMP450 × J7645Z) having low tannin, starch, cellulose, and hemicellulose content were used for feed. Similar to previous studies, the content of fat, lignin, and ash in both food and feed samples (hull-less) was nearly in the same range (Pontieri et al., 2010; Waniska et al., 2016; Beloshapka et al., 2016). Due to the non-digestible nature of hemicellulose and the bitter taste of tannin, these components are not desirable for food and feed, which leads to conclude that hull-less samples can be used as food (Table 1). The overall content of lignocelluloses in feed was higher than the food but significantly lower than the fuel grain samples. Based on required chemical components, whole grains/flours can be used as animal feed.
Moreover, owing to the high content of tannin, lignocelluloses, and ash, hulled grain hybrids (NW2, NW3, X098, 624A × NW2, AMP450 × NW2, 624A × X098, and 622A × NW3) were used as fuel. However, 622A × NW2, 622A × X098, AMP450 × X098, 624A × NW3, and AMP450 × NW3 hybrids were used for both feed and fuel due to desirable chemical content for feed and fuel (low tannin and high hemicellulose). The insoluble dietary fibers especially hemicelluloses are not seen as desirable as food (Fuller et al., 2016; Qiu et al., 2017; Stamenković et al., 2020). According to the previous studies, the hulled grains with an increasing tannin, cellulose, hemicellulose, lignin, and ash were used as an ideal additive for biofuel production (Table 1) (Eggum et al., 1983; Corredor et al., 2005; Stamenković et al., 2020). In this study, insoluble dietary fibers were found higher in the hulled grains with increasing content of starch, tannin, and ash and proposed as a supreme stabilizer for biofuel. Similar to the previous studies, hulled grains contained less starch while high fibers content (hemicellulose and lignin), tannin and ash compared with hull-less grains (Eggum et al., 1983).
A large CV in the sample datasets for starch, protein, fat, tannin, cellulose, hemicellulose, lignin, and ash was observed in this study. The previous studies found that the CV of starch, protein, fat, tannin, cellulose, hemicellulose, lignin, and ash in sorghum grains for bulk number of samples were 5.32, 13.2, 16.3, 9.09, 32, 30, 19, and 9.41%, respectively (De Alencar Figueiredo et al., 2006; Hill et al., 2012; Li et al., 2014; Wang et al., 2020). Compared with the previous studies, the higher CVs in this study indicate a wide variation among the sorghum grain samples and a large potential of model robustness when more samples from more hybrids and from more regions to be included.
PLSR Model Performance
For each chemical component and morphological type (whole grains, whole grain flours, hulled grain flours, and hull-less grain flours) of sorghum grain samples, spectral preprocessing yielded different results in the PLSR models. The predicted and reference values of whole grains, whole grain flours, hulled grain flours, and hull-less grain flours were consistent (Figures 4–7, respectively). Several physical phenomena contribute to the additive and multiplicative scatter effects, caused by the differences in the grain surface structure, i.e., hulled or hull-less grains. The spectral data were preprocessed to eliminate the influence of light scattering from the non-homogeneity of the sample particle distribution (Pohl and Senn, 2011). MSC can be considered as a suitable method when working with the samples constituted by particle size because it varies according to the grain hardness and uniformity to random noise (Sampaio et al., 2018). MSC has a tendency to produce outliers, while SNV induce curved structures derived from the treated spectra (Fearn et al., 2009). Both MSC and SNV methods assisted in getting models performance with different combinations of treatments for different chemicals.
The variations within and across the calibration and validation processes were maximized and repeated several times using the SPXY and KS sample selection methods (Vehtari et al., 2017). To demonstrate the effect of SPXY and KS optimization on the partitioning of calibration and validation subsets, PCA was applied and results indicate that the SPXY and KS methods could efficiently optimize NIRS results as well as the analytical properties of both the calibration and validation subsets (Table 3). These distributions and optimizations led to the good performance of FT-NIR models.
Lignin in whole grain flours showed reasonably accurate results using the first derivative and KS sample selection method while all other chemical components in whole grains showed precise results with a combination of the second derivative and SPXY sample selection method. Likewise, starch, protein, tannin, and ash in hull-less flours showed good results using the first derivative while all other morphological types of grains showed accurate results using the second derivative. Starch, fat, and tannin in whole grain flours, fat, and cellulose in hulled grain flours, and cellulose in hull-less grain flours showed reasonably precise results using KS sample selection method, compared with other chemicals in all the morphological types of sorghum grains (Table 2). Similar to the previous studies, spectral pretreatment especially derivative and spectra smoothing improved the models' performance and quality (Wolfrum et al., 2013).
The values of R2 and RPD for all the models were high enough to use them for future predictions (Malley et al., 2004). In contrary to the previous studies, the prediction accuracy (paired t-test between the measured and predicted values) of both the models of whole grains and flours was not different (Flinn et al., 1998; Garnsworthy et al., 2000; De Alencar Figueiredo et al., 2006) (Table 3; Figures 4, 5). However, it is suggested that the whole grains can be used instead of grinding because the goodness of fit of both grains and flours models showed. Moreover, the stiff texture of the hulls has a significant effect on the performance of calibration of biochemical components, indicating that the hulls may distort spectral information (De Alencar Figueiredo et al., 2006; Guindo et al., 2016) (Table 3; Figures 6, 7). To determine the prediction accuracy of the models, t-test was applied between the morphological types of sorghum grains, i.e., grains vs. flours, hulls vs. hull-less. It showed that the flour models were significantly better for prediction of protein, tannin, and ash while the hull-less models were significantly better for prediction of fat, hemicellulose, and lignin (Supplementary Figures 1, 2). However, the performance index of all the models demonstrated that these models performed well in the sample set and could be used for the estimation of future screenings and predictions, efficiently (Table 3).
The previous research showed the potential for analysis for the major chemical components and forms of grains according to their needs and interest. Nevertheless, models determining the content of the chemical components especially cellulose, hemicellulose, and lignin in sorghum grain flours are rarely investigated, nor systematically compared for different grain types. FT-NIR spectroscopy for grading grains based on grain quality parameters is preferable because it is non-destructive and non-hazardous. RPD and RPDcv showed that all the models could successfully predict the biochemicals, suggesting that the FT-NIR spectroscopy is a promising tool for predictions of sorghum grain biochemicals used for food, feed, and fuel (Table 3). The models presented in this paper can be applied to the industry level for producers to select grains with outstanding potential for food, feed, and fuel productions. This is a high throughput method with high accuracy as compared with the wet chemical determination methods to determine the desired chemical components for human food, animal feed, and biofuel. To improve model performance, it is necessary to extend the sample size over years so that prediction could be improved because the environmental factors induce high influence on the properties and chemical components of sorghum grains.
Conclusion
Biochemical components of sorghum were successfully predicted for enhancing grain sorting efficiency for food, feed, and fuel using FT-NIR spectroscopy. The PLS regression models allowed to screen and predict chemical components especially lignocellulose (cellulose, hemicellulose, and lignin) predictions that were rarely reported. Cellulose, and hemicellulose and lignin concentrations could be predicted using NIR data measured on whole grains, with high accuracy (R2 > 0.88); thus, whole grains are recommended instead of flours to save grinding time and labor. NIR-based PLS models for varied compositions of hulled and hull-less grain flours compared to the wet chemical procedures are desirable. The established PLSR models could enable food, feed, and fuel producers to efficiently evaluate a large number of samples by predicting the required biochemical components in sorghum grains.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
IE: conceptualization, methodology, and writing—original draft. SH: formal analysis and investigation. WL: data analysis, visualization, and results interpretation. NH: results interpretation and writing. CT, SL, and BD: resources. ML: software. GX: data curation and funding acquisition. KY: supervision, data analysis, results interpretation, writing, and funding acquisition. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by the Department of Energy Conservation and Technology Equipment of China's National Energy Administration [Science and Technology Department, No. (2012) 32] and funded by China Datang New Energy Co. LTD, Henan Tianguan Group Co., Ltd., and the CAU Basic Scientific Research Fund (2020RC037).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.720022/full#supplementary-material
References
Arana, I., Jarén, C., and Arazuri, S. (2005). Maturity, variety and origin determination in white grapes (Vitis vinifera L.) using near infrared reflectance technology. J. Near Infrared Spectrosc. 13, 349–357. doi: 10.1255/jnirs.566
Bader Ul Ain, H., Saeed, F., Khan, M. A., Niaz, B., Khan, S. G., Anjum, F. M., et al. (2019). Comparative study of chemical treatments in combination with extrusion for the partial conversion of wheat and sorghum insoluble fiber into soluble. Food Sci. Nutr. 7, 2059–2067. doi: 10.1002/fsn3.1041
Barros, F., Awika, J. M., and Rooney, L. W. (2012). Interaction of tannins and other sorghum phenolic compounds with starch and e ff ects on in vitro starch digestibility. J. Agric. Food Chem. 60, 11609–11617. doi: 10.1021/jf3034539
Beć, K. B., Grabska, J., and Huck, C. W. (2021). Principles and applications of miniaturized near-infrared (NIR) spectrometers. Chem. A Eur. J. 27, 1514–1532. doi: 10.1002/chem.202002838
Beloshapka, A., Buff, P., Fahey, G., and Swanson, K. (2016). Compositional analysis of whole grains, processed grains, grain co-products, and other carbohydrate sources with applicability to pet animal nutrition. Foods 5:23. doi: 10.3390/foods5020023
Brown, C. D., Vega-Montoto, L., and Wentzell, P. D. (2000). Derivative preprocessing and optimal corrections for baseline drift in multivariate calibration. Appl. Spectrosc. 54, 1055–1068. doi: 10.1366/0003702001950571
Caporaso, N., Whitworth, M. B., and Fisk, I. D. (2018). Near-infrared spectroscopy and hyperspectral imaging for non-destructive quality assessment of cereal grains. Appl. Spectrosc. Rev. 53, 667–687. doi: 10.1080/05704928.2018.1425214
Corredor, D. Y., Bean, S., Schober, T., and Wang, D. (2005). Effect of decortication of sorghum on ethanol production and composition of DDGS. 2005 ASAE Annu. Int. Meet. doi: 10.13031/2013.19582
De Alencar Figueiredo, L. F., Davrieux, F., Fliedel, G., Rami, J. F., Chantereau, J., Deu, M., et al. (2006). Development of NIRS equations for food grain quality traits through exploitation of a core collection of cultivated sorghum. J. Agric. Food Chem. 54, 8501–8509. doi: 10.1021/jf061054g
de Morais Cardoso, L., Pinheiro, S. S., Martino, H. S. D., and Pinheiro-Sant’Ana, H. M. (2017). Sorghum (Sorghum bicolor L.): Nutrients, bioactive compounds, and potential impact on human health. Crit. Rev. Food Sci. Nutr. 57, 372–390. doi: 10.1080/10408398.2014.887057
Eggum, B. O., Monowar, L., Bach Knudsen, K. E., Munck, L., and Axtell, J. (1983). Nutritional quality of sorghum and sorghum foods from Sudan. J. Cereal Sci. 1, 127–137. doi: 10.1016/S0733-5210(83)80030-7
Evers, T., and Millar, S. (2002). Cereal grain structure and development: some implications for quality. J. Cereal Sci. 36, 261–284. doi: 10.1006/jcrs.2002.0435
Fearn, T., Riccioli, C., Garrido-Varo, A., and Guerrero-Ginel, J. E. (2009). On the geometry of SNV and MSC. Chemom. Intell. Lab. Syst. 96, 22–26. doi: 10.1016/j.chemolab.2008.11.006
Ferreira, D. S., Pallone, J. A. L., and Poppi, R. J. (2015). Direct analysis of the main chemical constituents in Chenopodium quinoa grain using Fourier transform near-infrared spectroscopy. Food Control 48, 91–95. doi: 10.1016/j.foodcont.2014.04.016
Flinn, P. C., Black, R. G., Iyer, L., Brouwer, J. B., and Meares, C. (1998). Estimating the food processing characteristics of pulses by near infrared spectroscopy, using ground or whole samples. J. Near Infrared Spectrosc. 6, 213–220. doi: 10.1255/jnirs.139
Font, R., del Río-Celestino, M., and de Haro-Bailón, A. (2006). The use of near-infrared spectroscopy (NIRS) in the study of seed quality components in plant breeding programs. Ind. Crops Prod. 24, 307–313. doi: 10.1016/j.indcrop.2006.06.012
Food and Agriculture Organization of the United Nations (2020). FAOSTAT Statistical Database. Rome: FAO.
Fuller, S., Beck, E., Salman, H., and Tapsell, L. (2016). New horizons for the study of dietary fiber and health: a review. Plant Foods Hum. Nutr. 71, 1–12. doi: 10.1007/s11130-016-0529-6
Galvão, R. K. H., Araujo, M. C. U., Jos,é, G. E., Pontes, M. J. C., Silva, E. C., and Saldanha, T. C. B. (2005). A method for calibration and validation subset partitioning. Talanta 67, 736–740. doi: 10.1016/j.talanta.2005.03.025
Garnsworthy, P. C., Wiseman, J., and Fegeros, K. (2000). Prediction of chemical, nutritive and agronomic characteristics of wheat by near infrared spectroscopy. J. Agric. Sci. 135, 409–417. doi: 10.1017/S0021859699008382
Geladi, P., and Macdougall, D. (1985). Linearization and scatter-correction for near-infrared reflectance spectra of meat. Appl. Spectr. 39, 491–500. doi: 10.1366/0003702854248656
Guindo, D., Davrieux, F., Teme, N., Vaksmann, M., Bastianelli, D., Verdeil, J., et al. (2016). Pericarp thickness of sorghum whole grain is accurately predicted by NIRS and can affect the prediction of other grain quality parameters. J. Cereal Sci. 69, 218–227. doi: 10.1016/j.jcs.2016.03.008
Gyori, Z. (2017). “Chapter 11 - Corn: grain-quality characteristics and management of quality requirements,” in Woodhead Publishing Series in Food Science, Technology and Nutrition, eds C. Wrigley, I. Batey, and D. B. T.-C. G. (Second E. Miskelly (Woodhead Publishing), 257–290. doi: 10.1016/B978-0-08-100719-8.00011-5
Hell, J., Prückler, M., Danner, L., Henniges, U., Apprich, S., Rosenau, T., et al. (2016). A comparison between near-infrared (NIR) and mid-infrared (ATR-FTIR) spectroscopy for the multivariate determination of compositional properties in wheat bran samples. Food Control 60, 365–369. doi: 10.1016/j.foodcont.2015.08.003
Heredia-Olea, E., Pérez-Carrillo, E., and Serna-Saldívar, S. O. (2012). Effects of different acid hydrolyses on the conversion of sweet sorghum bagasse into C5 and C6 sugars and yeast inhibitors using response surface methodology. Bioresour. Technol. 119, 216–223. doi: 10.1016/j.biortech.2012.05.122
Hill, H., Slade Lee, L., and Henry, R. J. (2012). Variation in sorghum starch synthesis genes associated with differences in starch phenotype. Food Chem. 131, 175–183. doi: 10.1016/j.foodchem.2011.08.057
Hu, N., Li, W., Du, C., Zhang, Z., Gao, Y., Sun, Z., et al. (2021). Predicting micronutrients of wheat using hyperspectral imaging. Food Chem. 343:128473. doi: 10.1016/j.foodchem.2020.128473
ISO (1988). Sorghum — Determination of tannin content. Technical Committee ISO/TC 34, Agricultural food products. Available online at: iso.org/obp/ui/#iso:std:17487:en
Jin, X., Chen, X., Shi, C., Li, M., Guan, Y., Yeon, C., et al. (2017). Bioresource technology determination of hemicellulose, cellulose and lignin content using visible and near infrared spectroscopy in Miscanthus sinensis. Bioresour. Technol. 241, 603–609. doi: 10.1016/j.biortech.2017.05.047
Kamboj, U., Guha, P., and Mishra, S. (2017). Characterization of chickpea flour by near infrared spectroscopy and chemometrics characterization of chickpea flour by near infrared spectroscopy and chemometrics. Anal. Lett. 50, 1754–1766. doi: 10.1080/00032719.2016.1247163
Kennard, R. W., and Stone, L. A. (1969). Computer aided design of experiments. Technometrics 11, 137–148. doi: 10.1080/00401706.1969.10490666
Krasznai, D. J., Champagne Hartley, R., Roy, H. M., Champagne, P., and Cunningham, M. F. (2018). Compositional analysis of lignocellulosic biomass: conventional methodologies and future outlook. Crit. Rev. Biotechnol. 38, 199–217. doi: 10.1080/07388551.2017.1331336
Li, H., Liang, Y., Xu, Q., and Cao, D. (2009). Key wavelengths screening using competitive adaptive reweighted sampling method for multivariate calibration. Anal. Chim. Acta 648, 77–84. doi: 10.1016/j.aca.2009.06.046
Li, J., Danao, M. C., Chen, S., Li, S., and Brown, P. J. (2015). Prediction of starch content and ethanol yields of sorghum grain using near infrared spectroscopy. J. Near Infrared Spectr. 92, 85–92. doi: 10.1255/jnirs.1146
Li, M., Feng, S., Wu, L., Li, Y., Fan, C., Zhang, R., et al. (2014). Sugar-rich sweet sorghum is distinctively affected by wall polymer features for biomass digestibility and ethanol fermentation in bagasse. Bioresour. Technol. 167, 14–23. doi: 10.1016/j.biortech.2014.04.086
Li, M., Wang, J., Du, F., Diallo, B., and Xie, G. H. (2017). High-throughput analysis of chemical components and theoretical ethanol yield of dedicated bioenergy sorghum using dual-optimized partial least squares calibration models. Biotechnol. Biofuels 10, 1–16. doi: 10.1186/s13068-017-0892-z
Malley, D. F., Martin, P. D., and Ben-Dor, E. (2004). Near infrared spectroscopy in agriculture. Appl. Anal. soils 44, 729–784.
Mirfakhar, M., Asadollahi, M. A., Amiri, H., and Karimi, K. (2020). Co-fermentation of hemicellulosic hydrolysates and starch from sweet sorghum by Clostridium acetobutylicum: a synergistic effect for butanol production. Ind. Crops Prod. 151:112459. doi: 10.1016/j.indcrop.2020.112459
Orman, B. A., and Schumann, R. A. (1991). Comparison of near-infrared spectroscopy calibration methods for the prediction of protein, oil, and starch in maize grain. J. Agric. Food Chem. 39, 883–886. doi: 10.1021/jf00005a015
Padmore, J. M. (1990). “Fat (crude) or ether extract in animal feed. Method 920.39,” in Off. methods Anal. Assoc. Off. Anal. Chem. 15th ed. Arlington, VA.
Palacios, C. E., Nagai, A., Torres, P., Rodrigues, J. A., and Salatino, A. (2021). Contents of tannins of cultivars of sorghum cultivated in Brazil, as determined by four quantification methods. Food Chem. 337:127970. doi: 10.1016/j.foodchem.2020.127970
Pohl, F., and Senn, T. (2011). A Rapid and sensitive method for the evaluation of cereal grains in bioethanol production using near infrared reflectance spectroscopy. Bioresour. Technol. 102, 2834–2841. doi: 10.1016/j.biortech.2010.10.058
Pontieri, P., Maro, A., Di T.amburino, R., Stefano, M., De Tilley, M., Bean, S. R., et al. (2010). Chemical composition of selected food-grade sorghum varieties grown under typical Mediterranean Conditions 55, 139–143. Available online at: https://www.ars.usda.gov/ARSUserFiles/30200510/2010%20-%20Chemical%20Composition%20of%20Selected.pdf
Qiu, S., Yadav, M. P., and Yin, L. (2017). Characterization and functionalities study of hemicellulose and cellulose components isolated from sorghum bran, bagasse and biomass. Food Chem. 230, 225–233. doi: 10.1016/j.foodchem.2017.03.028
Ratnavathi, C. V., and Patil, J. V. (2014). Sorghum utilization as food. J. Nutr. Food Sci. 04, 1–8. doi: 10.4172/2155-9600.1000247
Ricardo Soccol, C., Faraco, V., Karp, S., Vandenberghe, L. P. S., Thomaz-Soccol, V., Woiciechowski, A., et al. (2011). “Chapter 5 - Lignocellulosic bioethanol: current status and future perspectives,” in Biofuels, eds. A. Pandey, C. Larroche, S. C. Ricke, C.-G. Dussap, and E. Gnansounou (Amsterdam: Academic Press) 101–122.
Sampaio, P. S., Soares, A., Castanho, A., Almeida, A. S., Oliveira, J., and Brites, C. (2018). Optimization of rice amylose determination by NIR-spectroscopy using PLS chemometrics algorithms. Food Chem. 242, 196–204. doi: 10.1016/j.foodchem.2017.09.058
Sihono, H. S., Wijaya, M. I., Marina, Y. M., Puspitasari, W., and Wahyono, T. (2019). Exploring the potential of sorghum for food, animal feed and bioenergy through induced mutation breeding. Int. Conf. Nucl. Capacit. Build. Educ. Res. Appl. 1, 22–28.
Siwela, M., Taylor, J. R. N., de Milliano, W. A. J., and Duodu, K. G. (2007). Occurrence and location of tannins in finger millet grain and antioxidant activity of different grain types. Cereal Chem. 84, 169–174. doi: 10.1094/CCHEM-84-2-0169
Sluiter, A., Hames, B., Ruiz, R. O., Scarlata, C., Sluiter, J., Templeton, D., et al. (2012). Determination of structural carbohydrates and lignin in biomass. Biomass Anal. Technol. Team Lab. Anal. Proced. 1–14.
Srivastava, S., Mishra, G., and Mishra, H. N. (2018). FTNIR-A robust diagnostic tool for the rapid detection of Rhyzopertha dominica and Sitophilus oryzae infestation and quality changes in stored rice grains. Food Bioprocess Technol. 11, 785–796. doi: 10.1007/s11947-017-2048-3
Stamenković, O. S., Siliveru, K., Veljkovi,ć, V. B., Banković-Ili,ć, I. B., Tasić, M. B., Ciampitti, I. A., et al. (2020). Production of biofuels from sorghum. Renew. Sustain. Energy Rev. 124:109769. doi: 10.1016/j.rser.2020.109769
Vehtari, A., Gelman, A., and Gabry, J. (2017). Practical Bayesian model evaluation using leave-one-out cross-validation and WAIC. Stat. Comput. 27, 1413–1432. doi: 10.1007/s11222-016-9696-4
Wang, X., Han, X., Li, L., and Zheng, X. (2020). Optimization for quantification of sorghum tannins by ferric ammonium citrate assay. Grain Oil Sci. Technol. 3, 146–153. doi: 10.1016/j.gaost.2020.07.001
Wang, Z. X., He, Q. P., and Wang, J. (2015). Comparison of variable selection methods for PLS-based soft sensor modeling. J. Process Control 26, 56–72. doi: 10.1016/j.jprocont.2015.01.003
Waniska, R. D., and Rooney, L. W. (2000). “Structure and chemistry of the sorghum caryopsis,” in Sorghum: Production, Agronomy, Chemistry and Utilization, eds. W. Smith and R. A. Frederiksen (New York, NY: Wiley & Sons) 649–688.
Waniska, R. D., Rooney, L. W., and McDonough, C. M. (2016). Sorghum: utilization. Encycl. Food Grains Second Ed. 3–4, 116–123. doi: 10.1016/B978-0-12-394437-5.00128-5
Wiedemair, V., Mair, D., Held, C., and Huck, C. W. (2019). Investigations into the use of handheld near-infrared spectrometer and novel semi-automated data analysis for the determination of protein content in different cultivars of Panicum miliaceum L. Talanta 205:120115. doi: 10.1016/j.talanta.2019.120115
Wolfrum, E., Payne, C., Stefaniak, T., and Rooney, W. (2013). Multivariate Calibration Models for Sorghum Composition using Near-Infrared Spectroscopy. Technical Report NREL/TP-510056838. Golden: National Renewable Energy Laboratory (NREL). Available online at: nrel.gov/docs/fy13osti/56838.pdf
Wu, L., Li, M., Huang, J., Zhang, H., Zou, W., Hu, S., et al. (2015). A near infrared spectroscopic assay for stalk soluble sugars, bagasse enzymatic saccharification and wall polymers in sweet sorghum. Bioresour. Technol. 177, 118–124. doi: 10.1016/j.biortech.2014.11.073
Keywords: FT-NIR, sorghum grains, biochemical composition, food, feed, fuel, PLSR
Citation: Ejaz I, He S, Li W, Hu N, Tang C, Li S, Li M, Diallo B, Xie G and Yu K (2021) Sorghum Grains Grading for Food, Feed, and Fuel Using NIR Spectroscopy. Front. Plant Sci. 12:720022. doi: 10.3389/fpls.2021.720022
Received: 03 June 2021; Accepted: 19 August 2021;
Published: 16 September 2021.
Edited by:
Panmanas Sirisomboon, King Mongkut's Institute of Technology Ladkrabang, ThailandReviewed by:
Christian Wolfgang Huck, University of Innsbruck, AustriaKhwantri Saengprachatanatug, Khon Kaen University, Thailand
Copyright © 2021 Ejaz, He, Li, Hu, Tang, Li, Li, Diallo, Xie and Yu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kang Yu, a2FuZy55dSYjeDAwMDQwO3R1bS5kZQ==
†These authors have contributed equally to this work