Yuguang Xiao
Xiaoshu Zhang2
Yanjie Li
Shu Diao- 1Research Institute of Subtropical Forestry, Chinese Academy of Forestry, Hangzhou, China
- 2School of Civil Engineering and Architecture, Xinxiang University, Xinxiang, China
- 3Research Institute of Landscape Plants, Guizhou Academy of Forestry, Guiyang, China
Currently the determination of cyanidin 3-rutinoside content in plant petals usually requires chemical assays or high performance liquid chromatography (HPLC), which are time-consuming and laborious. In this study, we aimed to develop a low-cost, high-throughput method to predict cyanidin 3-rutinoside content, and developed a cyanidin 3-rutinoside prediction model using near-infrared (NIR) spectroscopy combined with partial least squares regression (PLSR). We collected spectral data from Michelia crassipes (Magnoliaceae) tepals and used five different preprocessing methods and four variable selection algorithms to calibrate the PLSR model to determine the best prediction model. The results showed that (1) the PLSR model built by combining the blockScale (BS) preprocessing method and the Significance multivariate correlation (sMC) algorithm performed the best; (2) The model has a reliable prediction ability, with a coefficient of determination (R2) of 0.72, a root mean square error (RMSE) of 1.04%, and a residual prediction deviation (RPD) of 2.06. The model can be effectively used to predict the cyanidin 3-rutinoside content of the perianth slices of M. crassipes, providing an efficient method for the rapid determination of cyanidin 3-rutinoside content.
1 Introduction
Michelia crassipes Y.W. Law is an evergreen shrub or small tree, the only purple-flowered species in the genus Michelia, sporadically distributed in Guangdong, Hunan, Guangxi, Jiangxi, Guizhou and other provinces of China, and grows in dense forests on mountain slopes and in ravines at an altitude of 300-1000 m (Committee FoCE, 1996; Liu et al., 2002; Yang et al., 2003). The flower color of genus Michelia is mostly white or yellowish, while the tepals of M. crassipes are purplish-red or deep purple, so it is often used as an important parent for the improvement of the flower color of genus Michelia and is an excellent resource for flower viewing and flower color breeding (Liao, 2007; Shao et al., 2015b; Shao et al., 2015a; Shao et al., 2016; Chai et al., 2018).
Anthocyanins are a class of flavonoid that are widely found in plants in nature. These anthocyanins are multi-functional and can play an important role in protecting against ultraviolet radiation, coping with drought and fighting pathogens (Tohge and Fernie, 2017). As a water-soluble natural pigment, anthocyanins appear blue in alkaline cellular fluids and red under acidic conditions. Therefore, many plant organs such as flowers, leaves, and fruits appear purple, red, or blue, with a positive correlation between the shade of color and anthocyanin content (Tanaka et al., 2008; Li et al., 2014). Cyanidin 3-rutinoside (Cy3R), the main component of anthocyanins in the tepals of M. crassipes, plays an important role in the formation of purple color in the tepals of M. crassipes (Liu et al., 2020b). Previous studies have found that M. crassipes exhibits significant genetic diversity, with tepals of different individuals differing in color, all showing a purple hue (He et al., 2018; Xiao et al., 2023). The correlation between flower color phenotype and Cy3R content is expected to provide important basic information for revealing the mechanism of flower color formation in plants and related genetic analysis.
There are many traditional methods used to detect anthocyanins content in plant tissue, such as microwave method, pH differential method and high performance liquid chromatography (Lee et al., 2005; Chen et al., 2007; Rong et al., 2016). The results of these traditional methods are accurate, but they are time-consuming and cumbersome as they require a lot of labor and material resources during the experimental process (Dzhanfezova et al., 2020). In recent years, High-performance liquid chromatography (HPLC) has begun to be gradually used for the determination of anthocyanins content (Kim and Lee, 2020; Thuy et al., 2021), which is fast and simple to operate, but requires expensive instrumentation and cannot be quickly detected in the field (Liu et al., 2022). In addition, all of these methods require sample destruction, which makes it difficult to achieve non-destructive detection and has a certain impact on the environment (Firmani et al., 2019). Therefore, it is of great significance to develop simpler, rapid, and non-destructive methods for the determination of anthocyanins content.
Near-Infrared (NIR) spectroscopy is a fast, easy-to-use and non-destructive detection technique (Wetzel, 1998; Zhang et al., 2023) which utilizes the spectral information in the near-infrared wavelength band (800 - 2500 nm) to obtain chemical and structural information about a specimen (Rinnan and Rinnan, 2007). The origin of this technique dates back to the late 1850s (Butler, 1983). With continuous development and maturation, NIR spectroscopy is now widely used in the fields of food, medicine, agriculture and industry (Biancolillo et al., 2019; Abu-Khalaf and Hmidat, 2020; Prananto et al., 2020; Rossi and Lozano, 2020; Li et al., 2023; Trenfield et al., 2023). In recent years, researchers have begun to apply NIR spectroscopy to forestry. For example, Y Zhang, Q Luan, J Jiang and Y Li (Zhang et al., 2021) utilized near-infrared (NIR) spectroscopy combined with partial least squares regression (PLSR) to predict the malondialdehyde (MDA) content of slash pine needles in a real-time and rapid manner to understand plant stress. In addition, Zhang et al. (2023) utilized near-infrared (NIR) spectroscopy to non-destructively detect the sugar content of peach under various conditions.
NIR spectroscopic data can be obtained from NIR instruments. These data contain a lot of information about the physical and chemical properties of the molecules (Czarnecki et al., 2021). These data provide a valuable resource for analysis, but they are also accompanied by noise interference (Liu et al., 2020a). To effectively eliminate noise, preprocessing spectral data becomes a critical step in constructing chemometric models (Katsumoto et al., 2001). In addition, choosing appropriate variables (bands) can significantly improve the model performance (Ma et al., 2018). However, no studies have been reported on the prediction of anthocyanin content of M. crassipes tepals.
Therefore, the aim of this study was to (1) establish a model for predicting the content of cyanidin 3-rutinoside in M. crassipes tepals with the help of near-infrared spectroscopy combined with chemometrics; and (2) compare the model performance of different combinations of spectral preprocessing and variable selection methods. The established model for predicting the content of cyanidin 3-rutinoside can not only realize the rapid acquisition of the flower color phenotype of M. crassipes, but also provide a reference for the rapid and non-destructive detection of the content of cyanidin 3-rutinoside in other plant species.
2 Materials and methods
2.1 Plant materials
The plant materials used in this experiment were obtained from the germplasm resource nursery of the Chinese Academy of Forestry Research Institute of Subtropical Forestry (30° 3’ N, 119° 57’ E) and Guizhou Academy of Forestry (26° 30’ N, 106° 44’ E). Based on the results of the previous flower color survey of M. crassipes resources in the two locations, M. crassipes individuals with large differences in flower color were randomly selected. Samples were collected in the morning of April-May 2023 when the weather was clear. M. crassipes flowers at the bud stage (flower buds enlarged, bracts dehiscent, showing purple tepals) and at blooming stage (both rounds of tepals unfolded, with a large amount of pollen dispersed, but not browning and withering) were plucked together with their pedicels, and then wrapped around the pedicels at the point of fracture with wet paper towels, and carefully put into air-filled self-sealing bags, to prevent the petals from falling off by squeezing (Fu and Dai, 2016; Yuan et al., 2023). A total of 66 samples were brought back to the laboratory for NIR spectroscopy. The collected samples were stored in a refrigerator at -80°C for the subsequent determination of cyanidin 3-rutinoside content.
2.2 Determination of monomeric anthocyanin content
Spectrophotometric method, is considered as a valid alternative to HPLC method due to its simplicity, rapidity and economy (Lee et al., 2008). This method is similar to HPLC method in terms of accuracy of results (Lao and Giusti, 2016), therefore, spectrophotometric method was used in this study for the determination of anthocyanin content. Pre-prepared 1% hydrochloric acid-methanol solution for anthocyanin extraction was obtained as follows: 3 ml of 36% concentrated hydrochloric acid was aspirated with a pipette gun and fixed to 100 ml with methanol (Lin et al., 2011). Accurately weighed 0.25 g of the sample was cut into 10 ml centrifuge tubes, replenished with 1% hydrochloric acid-methanol solution to 8 ml, and extracted at a low temperature and protected from light at 4 °C for 48 h, during which time the centrifuge tubes were shaken 2-3 times. A 96-well plate was prepared with 1% hydrochloric acid-methanol solution as a blank control, and 200 μl of anthocyanin extract was taken, and the absorbance value was read at 530 nm with a microplate reader (SpectraMax iD5, Molecular Devices, USA), and three replicates were set for each sample. The standard curve was plotted by gradient dilution with cyanidin 3-rutinoside standard (≥95%) (Shanghai Yuanye Biotechnology Co., Ltd.). The content of cyanidin 3-rutinoside was calculated using the following formula:
Where: C = content of cyanidin 3-rutinoside (mg ml-1) in the measuring tube obtained from the standard curve; VT = total volume of anthocyanin extract (ml) = 8; V1 = volume of anthocyanin crude extract used in the addition of the sample (ml); W = fresh weight of the sample (g).
2.3 NIR spectrum measurements
Spectral raw data were determined using a portable near-infrared spectral analyzer (LF-2500, Spectral evolution, USA). The spectral range was 1000-2500 nm with a resolution of 6 nm. The outer petals of the collected petals were placed on the background board, and the handheld fiber-optic contact probe was used to directly scan the petals at different flower colors. In order to minimize noise contamination and to ensure accuracy, the probe was closely attached to the petal surface during the measurement, while standard whiteboard correction was performed in time. A total of 129 spectral data were measured. From the 129 spectral data, 103 data were randomly selected as the calibration set and 26 data as the validation set.
2.4 Spectral analysis methods
Spectra typically have a relatively low signal-to-noise ratio in this region of 2400-2500 nm, and this spectral region was removed in order to eliminate the effect of noise (Xu et al., 2018; Guo et al., 2021). Preprocessing of spectral data is necessary to further minimize the effects of instruments, probe offsets, and surroundings on spectral data and to maximize the spectral differences (Osborne et al., 1993; Qiu et al., 2022). In this study, six preprocessing methods were applied, namely Standard normal variate (SNV), Block scale (BS), Detrended variable (DET), and Block normalization (BN), Removal of polynomial trends and standard normal transformation (DET-SNV), Block scale and standard normal transformation (BS-SNV). Four variable selection methods are also applied: bounded variable elimination (bve) (Eén and Biere, 2005; Soos et al., 2020), genetic algorithm (ga) (Molajou et al., 2021), regularized elimination procedure, and rep) (Mehmood et al., 2011), Significance multivariate correlation (sMC) (Tran et al., 2014).
As a classical linear multivariate analysis algorithm, PLSR has been widely used in the field of spectral data modeling (Cheng and Sun, 2017). When the number of independent variables is large and multicollinearity exists among these independent variables, the use of traditional multiple regression methods may lead to a decrease in the predictive performance of the model (Ma et al., 2023; Yang et al., 2023). Also, in the face of a limited number of samples, traditional methods may increase the risk of overfitting. However, PLSR methods can address these challenges more effectively and provide a better way to solve the above problems. Therefore, in this study, we completed the construction of a prediction model for the content of cyanidin 3-rutinoside based on PLSR in combination with the above preprocessing methods. The number of latent variables (LVs) was optimized by Leave-one-out cross-validation (LOOCV). Meanwhile, we used the coefficient of determination (R2), the root mean square error (RMSE), residual prediction deviation (RPD) and number of LVs as metrics to evaluate the model performance (Jin et al., 2020; Hssaini et al., 2022). Among these metrics, the closer the R2 value is to 1, the better and more stable the model fit is. Whereas, the closer the RMSE value is to 0, the higher the RPD value is, the superior predictive performance of the model is indicated, and the number of LVs is less than 10 as much as possible to avoid overfitting the model (Guo et al., 2021; Hssaini et al., 2022). Identification of the spectral regions that have a significant impact on the model was performed by building the PLSR model in eight independent sessions. In each modeling, the dataset was randomly assigned and divided into a calibration set and a validation set in an 8:2 ratio.
2.5 Software tools
All data were completed analyzed on R software (v4.3.1). The R packages “pls” and “enpls” were used to construct the PLSR model; The “prospectr” package was used to manipulate NIR spectral data (Wehrens and Mevik, 2007; Stevens and Ramirez-Lopez, 2014; Xiao et al., 2019). All plotting was performed using the “ggplot2” package (Wickham, 2011).
3 Results
3.1 Features of spectra
Selected raw spectra of eight representative M. crassipes tepals are shown in Figure 1A. The spectra after SNV, BS, BN, DET, BS-SNV and DET-SNV pretreatment are shown in Figures 1B–G, respectively. By observing the raw spectra, it was found that the samples exhibited significant absorption characteristic peaks near the bands of about 1400 nm and 2100 nm, and this observation was similar to the spectra after applying SNV, BN, DET, and DET-SNV preprocessing. However, the spectra after applying the BS and BS-SNV treatments show a greater number of peaks with sharper morphology, exhibiting more pronounced volatility. Additional absorption peaks were observed even in the originally relatively smooth spectral region.
Figure 1 Spectra of M. crassipes tepals; (A) raw spectra; (B) SNV; (C) BS; (D) BN; (E) DET; (F) BS-SNV; (G) DET-SNV.
3.2 Statistical values for cyanidin 3-rutinoside
The quantitative analysis conducted in this study on the concentration of cyanidin 3-rutinoside within the tepals of M. crassipes is graphically represented in Figure 2, where the minimum value was 1.89, the maximum value was 10.83, and the mean value was 5.25 with a standard deviation of 2.11. The determined values of Cy3R content showed a wide range of variation, a result that facilitates the calibration of the model.
Figure 2 The boxplot of cyanidin 3-rutinoside content for M. crassipes samples. The boxes represent the interquartile range, the lines inside the boxes represent the medians, and the whiskers denote the lowest and highest values within 1.5 times the interquartile range. Each point indicates a value of cyanidin 3-rutinoside content.
3.3 Model performance
The effects of six different spectral data preprocessing methods with four variable selection strategies in PLSR models are summarized in Table 1, including performance metrics for both the calibration and validation sets. Among all models, the calibration set has an average R2 value of 0.68 and an average RMSE value of 1.18%, with the highest values of 0.68 (R2) and 1.20% (RMSE), and the lowest values of 0.67 (R2) and 1.16% (RMSE); while the validation set has an average R2 value of 0.73 and an average RMSE value of 1.03%, with the highest values of 0.75 (R2) and 1.09% (RMSE), and the lowest values were 0.69 (R2) and 1.01% (RMSE). In addition, the mean value of RPD values for all models was 1.65 with the highest value of 2.06 and the lowest value of 1.34; the number of LVs ranged between 3 and 15, with 13 models having a number of LVs greater than 10, which may be an overfitting phenomenon.
Table 1 Comparison of R2, RMSE and RPD values of calibration and validation sets of PLSR prediction models based on different spectral preprocessing and variable selection methods.
The performance of the models with SNV, DET and DET-SNV preprocessing methods was improved compared to the models without data preprocessing. Without the variable selection method, the model built by the DET-SNV preprocessing method had the highest performance with a calibration set R2 and RMSE of 0.68 and 1.18%, respectively, and an RPD value of 1.68. This was followed by the SNV, DET, and BN preprocessing methods. The BS and BS-SNV preprocessing methods had the worst model performance, with a calibration set R2 and RMSE were 0.67 and 1.19%, respectively.
When combining the four variable selection methods with all the preprocessing methods, the model performance was essentially similar. However, when combining the sMC variable selection methods with the BS preprocessing methods, the PLSR model performed best, with R2 and RMSE of 0.68 and 1.18% for the calibration set, and 0.72 and 1.04% for the validation set, with an RPD value of 2.06, and a number of LVs of 9.
3.4 Establishment of a predictive model for cyanidin 3-rutinoside content
Based on the results in Table 1, we used the BS preprocessing method and the sMC variable selection algorithm to construct a PLSR model for the prediction of Cy3R content. The constructed Cy3R prediction model was used to estimate the Cy3R content in the validation set, and the estimated values were compared with the actual chemical assay results. As shown in Figure 3A, we can observe that the relationship between the estimated values and the actual measured values is closer to a linear regression line, which means that the predicted values in the validation set are closer to the actual values and perform better relative to Figure 3B. Therefore, compared to the original full-spectrum model, the Cy3R prediction model utilizes only 9% of the spectral bands to achieve superior prediction results. Figure 4 illustrates the distribution of residuals for the two models. most of the residual values for the Cy3R prediction model fall in the range of -1 to 1, and only a few residual values are distributed between -2 and 2. Compared to that, most of the residuals of the original full-spectrum model are distributed in the range of -2 to 2. This indicates that the prediction performance of the Cy3R prediction model is more stable and accurate. Figure 5 shows the eight randomly selected key variables for the Cy3R prediction model when using the sMC variable selection method. Among them, the variables in the bands at 1094.2, 1113, 1383.5, 1874.7, and 2385.7 nm have extremely important effects on the construction of the prediction model. These bands play a key role in the modeling process and help to improve the accuracy and reliability of the predictions.
Figure 3 Scatterplot of predicted Cy3R content of M. crassipes tepals based on (A) block-scale-significance multivariate correlation (BS-sMC) algorithm combined with partial least squares regression (PLSR) modeling and (B) original full-length spectral PLSR modeling. The black dashed line indicates the predicted Cy3R values vs. measured values; the gray solid line is the linear regression line of the model; the error bars for each scatter indicate the prediction error obtained by eight random calibrations of the model.
Figure 4 Residual plots of predicted M. crassipes Cy3R content based on (A) the BS-sMC algorithm PLSR model and (B) the original full-length spectral PLSR model. The error bars of the predicted values represent the SDs derived from the eight simulation models.
Figure 5 Spectral effects of PLS model with 8 random runs.
4 Discussion
M. crassipes, as an excellent ornamental plant, usually needs to obtain a large amount of trait information during the selection and breeding process. Flower color is an important trait in ornamental plants, which is mainly affected by anthocyanin content (Zhang et al., 2022). Determination of the correlation between flower color phenotype and pigment composition can also provide an important basis for the study of flower color formation mechanism (Fu and Dai, 2016). Although the traditional determination of anthocyanin composition and content has accurate and reliable results, it is time-consuming and destructive to the plant, and it is not possible to monitor the long-term dynamics of a physiological index. Therefore, the aim of this study was to establish a PLSR model using NIR spectroscopy to estimate and predict the Cy3R content of M. crassipes tepals, which provides a reference for high-throughput analysis of plant phenotypes. In selective breeding, it is beneficial to obtain the required phenotypic trait information quickly and accelerate the breeding process. One of the most basic and widely used modeling methods for predicting plant physiological content in near-infrared spectroscopy is the partial least squares method. For example, Reuben et al. concluded that the PLSR model could accurately predict the total anthocyanin content of the peel (Buenafe et al., 2022). Olaoluwa et al. accurately predicted avocado ripeness parameters using NIR spectroscopy combined with the PLSR model, and their predictive model for both dry matter and moisture content achieved an R2 of 0.92, with RPD values of 2.19 and 2.06, respectively (Olarewaju et al., 2016).
The findings of this study indicate that the utilization of various spectral preprocessing techniques does not uniformly enhance the performance of the models developed. In fact, certain preprocessing methods may result in a diminution of predictive accuracy, aligning with the outcomes reported by Vašát et al. (2017). In this study, we investigated the relationship between NIR spectra and Cy3R content. We compared the performance of prediction models constructed by six different spectral preprocessing methods and four variable selection algorithms in combination with PLSR, and finally confirmed the combination of the BS preprocessing method and the sMC variable selection method as the best prediction model. The R2 and RMSE of this model were 0.72 and 1.04%, respectively. These values were lower than the results of Liu et al.’s study (R2 = 0.90, RMSE = 0.30%, RPD = 3.19) for the anthocyanin content of Prunus cerasifera leaves (Liu et al., 2019). This difference may stem from the different locations where the spectral data were collected. The leaves of Prunus cerasifera are relatively large and more easily spreadable, making spectral data collection relatively easy. However, in contrast, M. crassipes tepals have a smaller surface area and are irregularly shaped, making them less likely to spread. Therefore, when collecting spectral data from tepals, the fiber-optic probe may not be able to fit completely on their surfaces, which introduces potentially interfering information and reduces the accuracy of the Cy3R content prediction model. In addition, tepals have high moisture content, which may also further reduce the accuracy of the model (Agelet and Hurburgh, 2014; Manzoor et al., 2022).
Models with high R2 and low RMSE usually indicate that the difference between the model’s predicted values and the actual measured values is small. However, previous studies have shown that the RPD value is an important indicator for confirming whether a model is reliable or not (Saeys et al., 2005; Davey et al., 2009; Magwaza et al., 2012). It is generally accepted that an RPD value of less than 1.5 implies that the model is unreliable, a model with an RPD value between 1.5 and 2.0 is suitable for rough estimation only, a model with an RPD value between 2.0 and 2.5 is suitable for quantitative prediction, a model with an RPD value between 2.5 and 3.0 is considered good, and a model with an RPD value of more than 3.0 is highly satisfactory (Malley et al., 2000; Saeys et al., 2005; Zimmermann et al., 2007; Magwaza et al., 2012; Olarewaju et al., 2016). In this study, the predictive model built by the BS-sMC combination had an RPD value of 2.06 even though the difference in the R2 and RMSE values of the models built by the combination of other different preprocessing and variable selection methods was not considered significant. This means that the model is suitable for quantitative prediction and can be reliably used for prediction of Cy3R content in tepals. This finding proves its potential value in practical applications.
The collection of spectral data is always unavoidably contaminated by environmental noise, so it is important to select effective spectral information (Guo et al., 2020). Appropriate preprocessing of spectral data and variable selection can effectively improve the accuracy of the model and make the modeling task easier (Mishra et al., 2020). The central goal of the BS preprocessing approach is to equalize the effects between different blocks, which may have different scales and number of variables, through block scaling and block variance scaling. This helps to avoid any one block having a dominant influence on the modeling results (Mishra et al., 2021). Analyzing the spectrograms, it is observed that the spectra preprocessed using the BS method exhibit a heightened number of absorption peaks in comparison to spectra treated with alternative preprocessing techniques. This observation might suggest that the BS preprocessing aids in uncovering subtle spectral variances, previously obscured by noise, thereby augmenting the detectability of potential characteristic bands within the spectral data (Vašát et al., 2017). These additional characteristic bands are potentially valuable because they can provide additional quantitative information to the PLSR model. The enrichment of the data has the potential to enhance the stability of the model, as reflected in the significant improvement in the model RPD values. In addition, we found that the sMC algorithm is very effective in variable selection and helps to build a reliable predictive model. This algorithm has been successfully used in other studies to predict different chemical compositions, such as chlorophyll content of Sassafras tzumu leaves and malondialdehyde content of slash pine needles (Li et al., 2019; Zhang et al., 2021). sMC algorithm also revealed several important spectral features related to Cy3R in this study, including wavelengths of 1094.2, 1113.0, 1383.5, 1874.7, and 2385.7 nm. As reported by Kokaly et al. phenolics will exhibit spectral features in the range of 1000-1500 nm, with the larger phenolic compounds exhibiting spectral features near 1470 nm, which is caused by the presence of O-H bonds in their molecular structure (Kokaly and Skidmore, 2015). In addition, we observe that the residual values of the model are more tightly distributed within the horizontal bands. This suggests that our predictive model is more suitable for practical applications, as the narrower distribution bands imply better fitting accuracy and higher prediction accuracy. These results further validate the reliability and practicality of our established model.
The model constructed in this study utilizing near-infrared spectroscopy demonstrated promising predictive capabilities; however, there remains scope for further optimization of its performance. Importantly, the dataset acquired reflects merely a single temporal snapshot within a specific year, and the influence of environmental variables (e.g., light and temperature) on the phytochemical composition may introduce additional uncertainty into the predictive model. To enhance the model’s accuracy and reliability, future endeavors will encompass a repeatability assessment and a planned substantial increase in the sample size. These steps will facilitate more comprehensive inversion studies and the subsequent validation of the model’s predictions against laboratory analytical results.
5 Conclusions
In this study, a model for predicting the content of cyanidin 3-rutinoside in M. crassipes tepals was successfully constructed using NIR spectroscopy and PLSR. This model provides a non-destructive method for the rapid determination of cyanidin 3-rutinoside content in M. crassipes tepals. It is worth mentioning that the reliability of the model can be enhanced by using spectral preprocessing and variable selection methods. We clearly demonstrated that the PLSR model based on the combination of the BS preprocessing method and the sMC variable selection method exhibited the best performance. This study not only furnishes essential data for elucidating the biochemical mechanisms underlying flower color formation but also pioneers new pathways for the high-throughput quantitative analysis of flower color phenotypic traits. Moreover, the development of an efficacious predictive model for chemical composition markedly contributes an invaluable reference for the detection and analysis of cyanidin-3-rutinoside content across a broad spectrum of plant research domains, particularly in other plant species.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
YX: Data curation, Formal analysis, Investigation, Writing – original draft. XZ: Investigation, Writing – original draft. JL: Resources, Supervision, Writing – review & editing. HL: Resources, Writing – review & editing. JJ: Supervision, Writing – review & editing. YL: Conceptualization, Data curation, Writing – review & editing. SD: Conceptualization, Funding acquisition, Project administration, Resources, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The research was supported by Zhejiang Science and Technology Major Program on Agricultural New Variety Breeding” (2021C02071-3).
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.
References
Abu-Khalaf, N., Hmidat, M. (2020). Visible/Near Infrared (VIS/NIR) spectroscopy as an optical sensor for evaluating olive oil quality. Comput. Electron. Agric. 173, 105445. doi: 10.1016/j.compag.2020.105445
Agelet, L. E., Hurburgh, C. R., Jr. (2014). Limitations and current applications of Near Infrared Spectroscopy for single seed analysis. Talanta 121, 288–299. doi: 10.1016/j.talanta.2013.12.038
Biancolillo, A., Firmani, P., Bucci, R., Magrì, A., Marini, F. (2019). Determination of insect infestation on stored rice by near infrared (NIR) spectroscopy. Microchemical J. 145, 252–258. doi: 10.1016/j.microc.2018.10.049
Buenafe, R. J., Tiozon, J., Boyd, L. A., Sartagoda, K. J., Sreenivasulu, N. (2022). Mathematical modeling to predict rice's phenolic and mineral content through multispectral imaging. Food Chem. Adv. 1, 100141. doi: 10.1016/j.focha.2022.100141
Chai, Y., Hu, X., Zhang, D., Liu, X., Liu, C., Jin, X. (2018). Studies on Compatibility of lnterspecific Hybridization Between Michelia crassipes and M. Figo, M. maudiae, M. platypetala. Acta Hortic. Sin. 45, 1970–1978. doi: 10.16420/j.issn.0513-353x.2017-0779
Chen, F., Sun, Y., Zhao, G., Liao, X., Hu, X., Wu, J., et al. (2007). Optimization of ultrasound-assisted extraction of anthocyanins in red raspberries and identification of anthocyanins in extract using high-performance liquid chromatography–mass spectrometry. Ultrasonics Sonochemistry 14, 767–778. doi: 10.1016/j.ultsonch.2006.12.011
Cheng, J.-H., Sun, D.-W. (2017). Partial least squares regression (PLSR) applied to NIR and HSI spectral data modeling to predict chemical properties of fish muscle. Food Eng. Rev. 9, 36–49. doi: 10.1007/s12393-016-9147-1
Czarnecki, M. A., Beć, K. B., Grabska, J., Hofer, T. S., Ozaki, Y. (2021). “Overview of application of NIR spectroscopy to physical chemistry,” in Near-Infrared Spectroscopy: Theory, Spectral Analysis, Instrumentation, and Applications (Singapore: Springer), 297–330.
Davey, M. W., Saeys, W., Hof, E., Ramon, H., Swennen, R. L., Keulemans, J. (2009). Application of visible and near-infrared reflectance spectroscopy (Vis/NIRS) to determine carotenoid contents in banana (Musa spp.) fruit pulp. J. Agric. Food Chem. 57, 1742–1751. doi: 10.1021/jf803137d
Dzhanfezova, T., Barba-Espín, G., Müller, R., Joernsgaard, B., Hegelund, J. N., Madsen, B., et al. (2020). Anthocyanin profile, antioxidant activity and total phenolic content of a strawberry (Fragaria× ananassa Duch) genetic resource collection. Food Bioscience 36, 100620. doi: 10.1016/j.fbio.2020.100620
Eén, N., Biere, A. (2005). Effective preprocessing in SAT through variable and clause elimination. In Bacchus, F., Walsh, T. (eds.) Theory and Applications of Satisfiability Testing. (Berlin: Springer). doi: 10.1007/11499107_5
Firmani, P., De Luca, S., Bucci, R., Marini, F., Biancolillo, A. (2019). Near infrared (NIR) spectroscopy-based classification for the authentication of Darjeeling black tea. Food Control 100, 292–299. doi: 10.1016/j.foodcont.2019.02.006
Fu, J., Dai, S. (2016). Analysis of color phenotypic and pigment contents of chrysanthemum based on hyperspectral imaginc. J. Beijing Forestry Univ. 38, 88–98. doi: 10.13332/j.1000-1522.20150483
Guo, P., Li, T., Gao, H., Chen, X., Cui, Y., Huang, Y. (2021). Evaluating calibration and spectral variable selection methods for predicting three soil nutrients using Vis-NIR spectroscopy. Remote Sens. 13, 4000. doi: 10.3390/rs13194000
Guo, Y., Liu, C., Ye, R., Duan, Q. (2020). Advances on water quality detection by uv-vis spectroscopy. Appl. Sci. 10, 6874. doi: 10.3390/app10196874
He, X., Jiao, Z., Zheng, J., Dou, Q., Huang, L., Wang, B., et al. (2018). Construction of fingerprint of Michelia germplasm by fluorescent SSR markers. Mol. Plant Breed. 16, 4705–4714. doi: 10.13271/j.mpb.016.004705
Hssaini, L., Razouk, R., Bouslihim, Y. (2022). Rapid prediction of fig phenolic acids and flavonoids using mid-infrared spectroscopy combined with partial least square regression. Front. Plant Sci. 13, 13. doi: 10.3389/fpls.2022.782159
Jin, X., Li, S., Zhang, W., Zhu, J., Sun, J. (2020). Prediction of soil-available potassium content with visible near-infrared ray spectroscopy of different pretreatment transformations by the boosting algorithms. Appl. Sci. 10, 1520. doi: 10.3390/app10041520
Katsumoto, Y., Berry, J., Ozaki, Y. (2001). Modern pretreatment methods in NIR spectroscopy. Near Infrared Anal. 2, 29–36.
Kim, I., Lee, J. (2020). Variations in anthocyanin profiles and antioxidant activity of 12 genotypes of mulberry (Morus spp.) fruits and their changes during processing. Antioxidants 9, 242. doi: 10.3390/antiox9030242
Kokaly, R. F., Skidmore, A. K. (2015). Plant phenolics and absorption features in vegetation reflectance spectra near 1.66μm. Int. J. Appl. Earth Observation Geoinformation 43, 55–83. doi: 10.1016/j.jag.2015.01.010
Lao, F., Giusti, M. M. (2016). Quantification of purple corn (Zea mays L.) anthocyanins using spectrophotometric and HPLC approaches: Method comparison and correlation. Food Analytical Methods 9, 1367–1380. doi: 10.1007/s12161-015-0318-0
Lee, J., Durst, R. W., Wrolstad, R. E., Eisele, T., Giusti, M. M., Hach, J., et al. (2005). Determination of total monomeric anthocyanin pigment content of fruit juices, beverages, natural colorants, and wines by the pH differential method: collaborative study. J. AOAC Int. 88, 1269–1278. doi: 10.1093/jaoac/88.5.1269
Lee, J., Rennaker, C., Wrolstad, R. E. (2008). Correlation of two anthocyanin quantification methods: HPLC and spectrophotometric methods. Food Chem. 110, 782–786. doi: 10.1016/j.foodchem.2008.03.010
Li, W., Wang, B., Wang, M., Chen, M., Yin, J. M., Kaleri, G. M., et al. (2014). Cloning and characterization of a potato StAN11 gene involved in anthocyanin biosynthesis regulation. J. Integr. Plant Biol. 56, 364–372. doi: 10.1111/jipb.12136
Li, Y., Sun, Y., Jiang, J., Liu, J. (2019). Spectroscopic determination of leaf chlorophyll content and color for genetic selection on. Sassafras tzumu. Plant Methods 15, 1–11. doi: 10.1186/s13007-019-0458-0
Li, W., Du, Y., Wang, Y., Liu, Z., Zhen, J., Du, W., et al. (2023). Research on 0n-line eficient near-infrared spectral ecognition and autoatic sorting technology of waste textiles based on convolutional neural network. Spectrosc. Spectral Anal. 43, 2139–2145.
Liao, J. (2007). Selective breeding of new cultivars for the Michelia crassipes law. J. Northwest Forestry Univ. 02, 76–78.
Lin, W., Wang, D., Wang, H., Ji, G., Liu, W. (2011). A study on correlation and regression analysis of anthocyanin contents and color lndices in Chinese cabbage. J. Qingdao Agric. Univ. 28, 201–204.
Liu, D., Ouyang, S., Zeng, G. (2002). Rare ornamental flowering trees - M. crassipes. Plants 2, 18–19.
Liu, X., Liu, C., Shi, Z., Chang, Q. (2019). Comparison of prediction power of three multivariate calibrations for estimation of leaf anthocyanin content with visible spectroscopy in. Prunus cerasifera. PeerJ 7, e7997. doi: 10.7717/peerj.7997
Liu, C., Yang, S. X., Li, X., Xu, L., Deng, L. (2020a). Noise level penalizing robust Gaussian process regression for NIR spectroscopy quantitative analysis. Chemometrics Intelligent Lab. Syst. 201, 104014. doi: 10.1016/j.chemolab.2020.104014
Liu, C., Yu, Q., Li, Z., Jin, X., Xing, W. (2020b). Metabolic and transcriptomic analysis related to flavonoid biosynthesis during the color formation of Michelia crassipes tepal. Plant Physiol. Biochem. 155, 938–951. doi: 10.1016/j.plaphy.2020.06.050
Liu, X., Yu, J., Liu, C., Deng, X. (2022). Estimation of leaf anthocyanin content in Prunus cerasifera based on color indices and BP neural network. J. Northwest Forestry Univ. 37, 145–152.
Ma, Y., He, H., Wu, J., Wang, C., Chao, K., Huang, Q. (2018). Assessment of polysaccharides from mycelia of genus Ganoderma by mid-infrared and near-infrared spectroscopy. Sci. Rep. 8, 10. doi: 10.1038/s41598-017-18422-7
Ma, W., Xin, Z., Han, C., Sang, Z., Shang, K., Li, Y. (2023). Modeling of soil nitrogen content in maize field of Shanxi province by visible near-lnfrared spectroscopy. J. Shanxi Agric. Sci. 51, 750–755.
Magwaza, L. S., Opara, U. L., Terry, L. A., Landahl, S., Cronje, P. J., Nieuwoudt, H., et al. (2012). Prediction of ‘Nules Clementine’ mandarin susceptibility to rind breakdown disorder using Vis/NIR spectroscopy. Postharvest Biol. Technol. 74, 1–10. doi: 10.1016/j.postharvbio.2012.06.007
Malley, D., Lockhart, L., Wilkinson, P., Hauser, B. (2000). Determination of carbon, carbonate, nitrogen, and phosphorus in freshwater sediments by near-infrared reflectance spectroscopy: Rapid analysis and a check on conventional analytical methods. J. Paleolimnology 24, 415–425. doi: 10.1023/A:1008151421747
Manzoor, M. F., Hussain, A., Naumovski, N., Ranjha, M. M. A. N., Ahmad, N., Karrar, E., et al. (2022). A narrative review of recent advances in rapid assessment of anthocyanins in agricultural and food products. Front. Nutr. 9, 901342. doi: 10.3389/fnut.2022.901342
Mehmood, T., Martens, H., Sæbø, S., Warringer, J., Snipen, L. (2011). A Partial Least Squares based algorithm for parsimonious variable selection. Algorithms Mol. Biol. 6, 1–12. doi: 10.1186/1748-7188-6-27
Mishra, P., Roger, J.-M., Jouan-Rimbaud-Bouveresse, D., Biancolillo, A., Marini, F., Nordon, A., et al. (2021). Recent trends in multi-block data analysis in chemometrics for multi-source data integration. TrAC Trends Analytical Chem. 137, 116206. doi: 10.1016/j.trac.2021.116206
Mishra, P., Woltering, E., El Harchioui, N. (2020). Improved prediction of ‘Kent’mango firmness during ripening by near-infrared spectroscopy supported by interval partial least square regression. Infrared Phys. Technol. 110, 103459. doi: 10.1016/j.infrared.2020.103459
Molajou, A., Nourani, V., Afshar, A., Khosravi, M., Brysiewicz, A. (2021). Optimal design and feature selection by genetic algorithm for emotional artificial neural network (EANN) in rainfall-runoff modeling. Water Resour. Manage. 35, 2369–2384. doi: 10.1007/s11269-021-02818-2
Olarewaju, O. O., Bertling, I., Magwaza, L. S. (2016). Non-destructive evaluation of avocado fruit maturity using near infrared spectroscopy and PLS regression models. Scientia Hortic. 199, 229–236. doi: 10.1016/j.scienta.2015.12.047
Osborne, B. G., Fearn, T., Hindle, P. H. (1993). Practical NIR spectroscopy with applications in food and beverage analysis. (Harlow: Longman Scientific and Technical) 227.
Prananto, J. A., Minasny, B., Weaver, T. (2020). Near infrared (NIR) spectroscopy as a rapid and cost-effective method for nutrient analysis of plant leaf tissues. Adv. Agron. 164, 1–49. doi: 10.1016/bs.agron.2020.06.001
Qiu, L., Zhang, M., Mujumdar, A. S., Chang, L. (2022). Convenient use of near-infrared spectroscopy to indirectly predict the antioxidant activitiy of edible rose (Rose chinensis Jacq “Crimsin Glory” HT) petals during infrared drying. Food Chem. 369, 130951. doi: 10.1016/j.foodchem.2021.130951
Rinnan, R., Rinnan, Å. (2007). Application of near infrared reflectance (NIR) and fluorescence spectroscopy to analysis of microbiological and chemical properties of arctic soil. Soil Biol. Biochem. 39, 1664–1673. doi: 10.1016/j.soilbio.2007.01.022
Rong, C., Zhang, X., Li, T., Hu, J., Zhang, J., Li, Z. (2016). Optimization of microwave-assisted extraction of proanthocyanidins from Rosa davurica pall. Seeds by response surface methodology. Food Sci. 37, 41–46.
Rossi, G. B., Lozano, V. A. (2020). Simultaneous determination of quality parameters in yerba mate (Ilex paraguariensis) samples by application of near-infrared (NIR) spectroscopy and partial least squares (PLS). Lwt 126, 109290. doi: 10.1016/j.lwt.2020.109290
Saeys, W., Mouazen, A. M., Ramon, H. (2005). Potential for onsite and online analysis of pig manure using visible and near infrared reflectance spectroscopy. Biosyst. Eng. 91, 393–402. doi: 10.1016/j.biosystemseng.2005.05.001
Shao, W., Jiang, J., Dong, R. (2015a). A new Michelia cultivar 'Mengzi'. Acta Hortic. Sin. 42, 1863–1864. doi: 10.19433/j.cnki.1006-9119.2015.22.010
Shao, W., Jiang, J., Dong, R. (2016). A new Michelia cultivar 'Mengxing'. Acta Hortic. Sin. 43, 1219–1220. doi: 10.16420/j.issn.0513-353x.2015-0888
Shao, W., Jiang, J., Dong, R., Luan, Q. (2015b). A new variety, Michelia 'Mengyuan'. SCIENTIA Silvae SINICAE 51, 155.
Soos, M., Gocht, S., Meel, K. S. (2020). Tinted, detached, and lazy CNF-XOR solving and its applications to counting and sampling. In: Lahiri, S., Wang, C. (eds) Computer Aided Verification. (Cham: Springer).
Tanaka, Y., Sasaki, N., Ohmiya, A. (2008). Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 54, 733–749. doi: 10.1111/j.1365-313X.2008.03447.x
Thuy, N. M., Minh, V. Q., Ben, T. C., Thi Nguyen, M. T., Ha, H. T. N., Tai, N. V. (2021). Identification of anthocyanin compounds in butterfly pea flowers (Clitoria ternatea L.) by ultra performance liquid chromatography/ultraviolet coupled to mass spectrometry. Molecules 26, 4539. doi: 10.3390/molecules26154539
Tohge, T., Fernie, A. R. (2017). Leveraging natural variance towards enhanced understanding of phytochemical sunscreens. Trends Plant Sci. 22, 308–315. doi: 10.1016/j.tplants.2017.01.003
Tran, T. N., Afanador, N. L., Buydens, L. M., Blanchet, L. (2014). Interpretation of variable importance in partial least squares with significance multivariate correlation (sMC). Chemometrics Intelligent Lab. Syst. 138, 153–160. doi: 10.1016/j.chemolab.2014.08.005
Trenfield, S. J., Xu, X., Goyanes, A., Rowland, M., Wilsdon, D., Gaisford, S., et al. (2023). Releasing fast and slow: Non-destructive prediction of density and drug release from SLS 3D printed tablets using NIR spectroscopy. Int. J. Pharmaceutics: X 5, 100148. doi: 10.1016/j.ijpx.2022.100148
Vašát, R., Kodešová, R., Klement, A., Borůvka, L. (2017). Simple but efficient signal pre-processing in soil organic carbon spectroscopic estimation. Geoderma 298, 46–53. doi: 10.1016/j.geoderma.2017.03.012
Wehrens, R., Mevik, B.-H. (2007). The pls package: principal component and partial least squares regression in R. Journal of Statistical Software 18 (2), 1–23. doi: 10.18637/jss.v018.i02
Wetzel, D. L. (1998). “Analytical near infrared spectroscopy,” in Developments in Food Science, vol. 39. (Elsevier), 141–194. doi: 10.1016/S0167-4501(98)80009-5
Wickham, H. (2016). ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York. doi: 10.1007/978-3-319-24277-4_9
Xiao, N., Cao, D., Li, M., Xu, Q. (2019) enpls: Ensemble Partial Least Squares Regression. Available online at: cloudr-project.org.
Xiao, Y., Jiang, X., Lu, C., Liu, J., Diao, S., Jiang, J. (2023). Genetic diversity and population structure analysis in the Chinese endemic species Michelia crassipes based on SSR markers. Forests 14, 508. doi: 10.3390/f14030508
Xu, D., Ma, W., Chen, S., Jiang, Q., He, K., Shi, Z. (2018). Assessment of important soil properties related to Chinese Soil Taxonomy based on vis–NIR reflectance spectroscopy. Comput. Electron. Agric. 144, 1–8. doi: 10.1016/j.compag.2017.11.029
Yang, C., Chen, J., Fang, X. (2003). Yang l: Excellent garden tree species Michelia crassipes. Guizhou Forestry Sci. Technol. 6, 16–18.
Yang, F., Li, R., Feng, H., Li, T., Wang, G. (2023). Comparison of Hyperspectral Remote Sensing nversion Methods for Plant Nitrogen Content in different growth stages. J. Northeast Agric. Sci. 48, 118–124. doi: 10.16423/j.cnki.1003-8701.2023.03.025
Yuan, J., Jin, X., Zhang, Z., Xiao, Y., Yu, Q., Hu, X. (2023). Volatility Components of Michelia crassipes Tepals at Different Flowering Stages. Acta Hortic. Sin. 50, 1095–1109. doi: 10.16420/j.issn.0513-353x.2022-0176
Zhang, Y., Luan, Q., Jiang, J., Li, Y. (2021). Prediction and utilization of malondialdehyde in exotic pine under drought stress using near-infrared spectroscopy. Front. Plant Sci. 12, 735275. doi: 10.3389/fpls.2021.735275
Zhang, H., Wang, H., Li, Y., Gao, J., Yuan, X., Wang, L., et al. (2022). The Chemical Composition and Transcriptome Analysis Reveal the Mechanism of Color Formation in Rosa Hybrida cv 'Double delight'. Chin. Bull. Bot. 57, 649–660.
Zhang, X., Zhu, Y., Zhao, Y., Chen, M., Sun, Q., Xie, B., et al. (2023). Optimization of nondestructive testing method for soluble solid content of peach based on visible/near intrared spectroscopy. Acta Agriculturae Zhejiangensis 35, 1617–1625. doi: 10.3969/j.issn.1004-1524.20220862
Keywords: model calibration, NIR spectroscopy, regression algorithm, cyanidin 3-rutinoside, Michelia crassipes
Citation: Xiao Y, Zhang X, Liu J, Li H, Jiang J, Li Y and Diao S (2024) Prediction of cyanidin 3-rutinoside content in Michelia crassipes based on near-infrared spectroscopic techniques. Front. Plant Sci. 15:1346192. doi: 10.3389/fpls.2024.1346192
Received: 01 December 2023; Accepted: 17 April 2024;
Published: 03 May 2024.
Edited by:
Tanima Bhattacharya, Lincoln University College, MalaysiaReviewed by:
Mun Fei Yam, University of Science Malaysia (USM), MalaysiaLahcen Hssaini, Institut National de la Recherche Agronomique, Morocco
Copyright © 2024 Xiao, Zhang, Liu, Li, Jiang, Li and Diao. 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: Shu Diao, diaoshu0802@163.com