Next Article in Journal
Poultry-Based Amendments and Cover Crop Residues Enhance Nutrient Cycling and Soil Health in Greenhouse Conditions
Previous Article in Journal
Selective Retention of Cross-Fertilised Fruitlets during Premature Fruit Drop of Hass Avocado
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determination of Saffron Flower Metabolites by Near-Infrared Spectroscopy for Quality Control

by
Jorge F. Escobar-Talavera
,
María Esther Martínez-Navarro
,
Gonzalo L. Alonso
* and
Rosario Sánchez-Gómez
*
Cátedra de Química Agrícola, Escuela Técnica Superior de Ingeniería Agronómica y de Montes y Biotecnología (ETSIAMB), Universidad de Castilla-La Mancha, Avda. de España s/n, 02071 Albacete, Spain
*
Authors to whom correspondence should be addressed.
Horticulturae 2024, 10(6), 593; https://doi.org/10.3390/horticulturae10060593
Submission received: 30 April 2024 / Revised: 30 May 2024 / Accepted: 3 June 2024 / Published: 6 June 2024
(This article belongs to the Section Postharvest Biology, Quality, Safety, and Technology)

Abstract

:
Saffron, obtained by dehydrating the stigmas of the Crocus sativus flower, is a spice of great importance. In saffron, the flower emerges before leaf formation, thanks to the nutritional reserves of the corm. Early knowledge of metabolite levels such as crocins, picrocrocin, safranal, anthocyanins, or kaempferols in flowers serves as a guide to evaluate the quality of the corm (coloring power, flavor, aroma, or antioxidant capacity, among others). In this study, near-infrared spectroscopy (NIR) was calibrated and validated to determine the main saffron metabolites, both in stigmas and in floral residue. To achieve this, saffron flowers from different locations of the Denomination of Origin (D.O.) “Azafrán de La Mancha” (Castilla-La Mancha, Spain) were analyzed using NIR spectroscopy. Prior to this, samples were analyzed by RP-HPLC-DAD, where the concentration of all cited metabolites was determined. The development of a predictive model through NIR calibration and validation was successful, achieving high R2 values, especially in the case of the sum of crocins and kaempferol-3-O-β-sophoroside. Using these predictive models, it is possible to determine the quality of saffron corm by analyzing the flower.

1. Introduction

The plant Crocus sativus L. has been cultivated since ancient times to obtain the dehydrated stigmas of its flowers, which constitute the spice saffron, known for its culinary use and various bioactive properties [1]. As with almost all spices, it cannot be determined whether it was first used as a medicinal plant or as a culinary condiment. Saffron is valued in cooking for its high content of crocins, substances responsible for its color; picrocrocin, the compound responsible for its bitter taste; and safranal, the main substance responsible for its exquisite aroma [2,3]. When ingested, crocetin, the hydrolysis product of crocins, is detected in the blood [4]. It has been shown that this molecule, in its trans form, acts as a blood oxygenator and is capable of crossing the blood–brain barrier [5]. Together with safranal, crocetin is attributed with the bioactive properties of saffron [6]
The floral remnants are a by-product consisting of tepals and stamens [7] which mainly contain anthocyanins and flavonols. These compounds are applied in the cosmetics and parapharmacy industries due to their antioxidant and bioactive properties [8].
In autumn, the flowers emerge from the corm (thickened stem) (Figure 1) before the leaves, a condition known as hysteranthy. The corm is planted at the beginning of summer, and during the autumn, winter, and spring months, daughter corms develop and grow, which then bloom the following autumn, while the mother corm dies. In other words, the plant remains alive for two years, with the first year dedicated to growth and development and the second year to offspring production. Normally, in plots where Crocus sativus L. is cultivated, the corms are kept buried until the plantation density is high, and soil nutrient competition prevents the formed corms from growing. This timeframe depends on the fertility characteristics of the production areas. Thus, in Spain, the same corms remain buried for 4 to 5 years. Crocus sativus L. is a sterile plant and does not reproduce by seeds but rather through budding [9]. Corms obtained from previous-year plantations act as seeds for new plantations, limiting the rapid spread of cultivation due to their availability. Commercially, corm suppliers certify their good sanitary conditions, but it is not possible guarantee their quality.
Currently, there is no analytical methodology to assure corm quality, since it can only be measured through the content of metabolites in the product, in both stigmas and floral remnants. Therefore, to indirectly control the quality of a corm, it is necessary to determine the concentration of the main metabolites such as crocins, picrocrocin [10], flavonols, and anthocyanins [11] in the saffron flowers directly emerging from the corm. This is not the case for safranal, since it is a compound that forms during the dehydration process and therefore does not directly derive from the corm.
The non-targeted near-infrared spectroscopy (NIR) technique can swiftly provide the metabolic fingerprint of agricultural products. Combining this technique with chemometrics can facilitate the data analysis, discrimination, and classification of agricultural products, as well as the prediction of concentrations and parameters related to their quality [12,13,14]. Supervised chemometric methods are widely known, including Linear Discriminant Analysis (LDA), Partial Least Squares Discriminant Analysis (PLS-DA), partial least squares regression (PLS-R), Soft Independent Modeling of Class Analogy (SIMCA), Support Vector Machine (SVM), Random Forests (RFs), Artificial Neural Networks (ANNs), K-Nearest Neighbors (K-NN), Linear Learning Machine (LLM), and Bayes Linear Discriminant Analysis [15,16], among others. When assessing the quality of an agricultural product, calibration is necessary, which can be performed through univariate, multiple regression, principal component, least squares, cross-validation, and independent test set calibrations, as well as through minimum angle regression and elastic net. The goodness of the models can be evaluated using certain parameters such as the coefficient of regression (R2), the Prediction Residual Sum of Squares (PRESS), the residual predictive deviation (RPD), the Standard Error of Prediction (SEP), or the standard error of cross-validation (SECV) [13], principally.
Nowadays, the analysis of the main metabolites of the Crocus sativus L. flowers, which encompasses stigmas or floral residue, involves analysis by chromato-graphic techniques such as Reversed-Phase High-Performance Liquid Chromatography with Diode Array Detection (RP-HPLC-DAD) [11,17]. As is known, these types of techniques are (a) time-consuming; (b) labor-intensive methods; and (c) require sample treatment. Therefore, scaling them up to a large number of samples is not an easy task. In this line, replacing their use by employing the NIR spectroscopy technique for the estimation of any of the main metabolites presents a better option for routine analysis in time due to the following benefits: (a) a fast technique; (b) no sample preparation required, since solids and liquids can be used in pure forms; (c) a low cost per sample, as no chemicals or solvents are needed; (d) an environmentally friendly technique, due to no waste being generated. Taking all of the above into account, the aim of this study is to calibrate NIR equipment to evaluate the quality of Crocus sativus L. corms by calibrating the content of the main metabolites in their flowers.

2. Materials and Methods

2.1. Chemicals and Reagents

The solvents and reagents employed were all HPLC purity or analytical grade. Acetonitrile (ACN) (CAS: 75-05-8), trifluoroacetic acid (TFA) (CAS: 7605-1), and hydrochloric acid (HCl) (CAS: 7647-01-0) were supplied from Panreac (Barcelona, Spain). Ultrahigh-purity water was produced using a Milli-Q system (Millipore, Bedford, MA, USA).
For standards, crocetin esters, trans-crocetin di-(β-D-gentiobiosyl) ester (trans-4-GG) (CAS: 55750-85-1), and trans-crocetin (β-D-glucosyl)-(β-D-gentiobiosyl) ester (trans-3-Gg) (CAS: 55750-84-0) with purity levels of ≥99% were obtained from Phytolab GmbH & Co. KG (Vestenbergsgreuth, Germany); picrocrocin (synthesized in the laboratory according to Sánchez et al. [3]) and kaempferol-3-O-β-sophoroside (CAS: 19895-95-5) were obtained from BLD Pharmatech (Shanghai, China).
For standards’ preparation, a stock solution was prepared, from which, dilutions were made to obtain the different calibration points. For the dilutions, water/hydrochloric acid (100:1, v/v) was used in the same way as for the extractions. Storage took place at −20 °C.

2.2. Plant Material

A total of 173 samples of Crocus sativus flowers from different locations within the D.O. “Azafrán de La Mancha” (Castilla-La Mancha, Spain) were collected during the 2022 and 2023 harvest seasons. The fresh Crocus sativus L. flowers were subjected to a freeze-drying process in the LyoAlfa 6-50 freeze-dryer (Telstar, Terrasa, Spain) for 5 days to ensure complete lyophilization (constant weight). The freeze-dryer conditions were −50 ± 2 °C and 10−3 mbar. Subsequently, dried Crocus sativus L. flowers were ground into a homogeneous powder using a mortar and stored at room temperature (18 ± 3 °C) in a chamber with silica gel until analysis.
Moisture content analysis was conducted on the samples using a moisture balance equipped with a halogen lamp, specifically the XM-120 T model (Cobos, Barcelona, Spain), operating at a temperature of 105 °C. When moisture loss was less than 0.1% in 180 s, it was considered that the samples had reached constant mass.

2.3. Preparation of Extracts of Crocus sativus L. Flowers

Extracts of Crocus sativus L. flowers were prepared according to Moratalla-Lopez et al.’s method [11]. Briefly, 200 mg of lyophilized and grounded powder from each sample was mixed with 25 mL of water/HCl (100:1, v/v) and stirred in the dark at 500 rpm for 1 h. Afterward, extracts underwent centrifugation at 3500 rpm for 5 min (Selecta, Barcelona, Spain). The supernatant was then filtered through a hydrophilic polytetrafluoroethylene (PTFE) membrane with pores of 0.45 μm (Millipore, Bedford, MA, United States). The resulting filtered extract was carefully transferred to a vial for Reversed-Phase High-Performance Liquid Chromatography with Diode Array Detection (RP-HPLC-DAD) analysis.

2.4. Chromatographic Conditions

The content of Crocus sativus L. flowers’ metabolites was determined according to a method developed by Moratalla-López et al. [11]. The RP-HPLC-DAD analyses were performed using an Agilent 1200 HPLC chromatograph (Palo Alto, CA, United States). The chromatographic column employed was a Develosil ODS-HG-5 with dimensions of 250 × 4.6 mm, 5 µm, (Teknokroma, Sant Cugat Del Vallès, Barcelona, Spain), maintained at a temperature of 40 °C. The mobile phase consisted of ultrahigh-purity water/trifluoroacetic acid (TFA) (99.5:0.5, v/v) (A) and acetonitrile (ACN) (B). The elution gradient was as indicated in Table 1. The flow rate was set at 1.0 mL/min, and the injection volume was 30 μL.
For all the compounds, identification was performed using the DAD by comparison with the corresponding UV–visible spectra and the retention time of their pure standards. Thus, the DAD detector was set at 440 nm to identify crocetin esters (trans-crocetin (β-D-triglucosyl)-(β-D-gentiobiosyl) ester, trans-5-tG; trans-crocetin (β-D-neapolitanosyl)-(β-D-glucosyl), trans-5-ng; trans-crocetin di-(β-D-gentiobiosyl) ester, trans-4-GG; trans-crocetin (β-D-neapolitanosyl)-(β-D-glucosyl) ester, trans-4-ng; trans-crocetin (β-D-gentiobiosyl), trans-2-G; trans-crocetin (β-D-glucosyl)-(β-D-gentiobiosyl) ester, trans-3-Gg; cis-crocetin di-(β-D-gentiobiosyl) ester, cis-4-GG; cis-crocetin (β-D-glucosyl)-(β-D-gentiobiosyl) ester, cis-3-Gg; cis-crocetin di-(β-D-gentiobiosyl) ester, cis-4-GG; cis-crocetin di-(β-D-glucosyl) ester, cis-2-gg; and, cis-crocetin (β-D-gentiobiosyl) ester, cis-2-G), at 250 nm for the detection of picrocrocin, at 266 nm for the detection of flavonols (kaempferol-3-O-β-sophoroside-7-O-β-glucoside, K-3-O-s-7-O-g; kaempferol-3-O-β-sophoroside, K-3-O-s; kaempferol aglycone, K), and at 520 nm for the identification of anthocyanins (delphinidin-3,5-di-O-β-glucoside, D-3,5-di-O-g; petunidin-3,5-di-O-β-glucoside, P-3,5-di-O-g; delphinidin-3-O-β-glucoside, D-3-O-g; malvidin-3,5-diO-β-glucoside, M-3,5-di-O-g; and petunidin-3-O-β-glucoside, P-3-O-g). Quantification was based on the calibration curves of the respective commercial standards at five different concentrations achieved by a UV–visible signal (R2 = 0.9900–0.9999). All analyses were performed in duplicate, with two measurements taken for each replicate [2,18].

2.5. Near-Infrared Spectroscopy Measurements

The powdered saffron flower samples were deposited to completely cover the surface of a 1 cm2 quartz Petri dish (Perkin-Elmer, Norwalk, CT, USA) and analyzed by Perkin Elmer Spectrum One FT-NIR equipment (Norwalk, CT, USA) coupled with a near-infrared reflectance accessory (NIRA). Data acquisition was carried out over a wavelength range of 10,000–4000 cm−1, with a fixed resolution of 16 cm−1. All samples were scanned in duplicate.

2.6. Data Analysis

The model calibration analysis was performed with the average of two replicated spectra for each Crocus sativus L. flower sample. In terms of data processing, a multivariate analysis approach was employed to conduct both quantitative and qualitative assessments. The tracking procedure in the partial least squares regression (PLS-R) modeling generated a systematic profile that incorporated various spectral pretreatment methods to mitigate diverse adverse effects stemming from the physical properties of the sample, technical errors during measurements, or simply instrumental noise. Spectrum-One software version 1.00 was utilized to perform data pretreatment and establish PLS-R models for the studied content. The pretreatment combinations of the spectrum varied depending on the calibrated metabolites, including standard normal variate (SNV), which was applied to all spectra used in the calibration models, which first normalizes a spectrum by calculating the average intensity value and then subtracting this value from the spectrum; smoothing (5, 13, and 19 points); and Savitzky–Golay first and second derivative (5, 9, and 13 points) pre-processing for all measurements to obtain the most robust model.
The model’s performance was examined using the coefficient of determination for calibration, R2c (1), which indicates the proportion of variance in the dependent variable that the independent one can explain, and the residual predictive deviation, RPD (2); by providing a metric of model validity, higher values correspond to better predictive capacity in the model, as previously described by Luo et al. [19]
R 2 c = i ( y i f i ) 2 / i ( y i y ¯ ) 2
R P D = [ S t a n d a r d   d e v i a t i o n   o f   m e a s u r e d   e x t r a c t s / R M S E ]
where RMSE is the root mean square error.
Data pre-processing methods and the selection of wavenumber ranges resulted in high predictability and the precise estimation of Crocus sativus L. flower metabolites.
For the calibration process, cross-validation was employed to determine the most suitable number of latent variables and prevent excessive fitting of the equation. Alongside r2, the ratio of standard deviation to standard error of cross-validation, RPDc (3), was used to evaluate the overall accuracy of each equation’s fit according to Suchat et al. [20]
R P D c = S D c a l S E C V
The RPDp (4) index was also calculated. This index represents the relationship between the standard deviation of the reference data and the standard error of cross-validation or test set validation, but specifically using the predicted reference values instead of the actual values. This index provides a measure of the predictive capability of the model relative to the variability in the reference data, which can help determine the usefulness of the model for specific applications. A higher RPDp value indicates the better predictive capability of the model according to Suchat et al. [20].
R P D p = S D c a l S E P
where SDcal refers to the standard deviation obtained from the mean of each parameter of the samples and SEP refers to the Standard Error of Prediction in calibration.
The RMSEC (Root Mean Square Error of Calibration) (5) was also calculated, which is a measure used to evaluate the quality of a calibration model. The lower the RMSEC value, the better the predictive capability of the model. This value represents the standard deviation of the residuals between the predicted values and the actual values in the calibration dataset.
R M S E C = i = 0 n f i y i 2 / n 1 / 2

3. Results and Discussion

3.1. Crocus sativus Flower Metabolite Content (RP-HPLC-DAD) for Set of Samples

The entire set of Crocus sativus L. flower samples (173) used in the present study, from all the localities included under the D.O. “Azafrán de La Mancha”, were of two ages (2022–2023) and submitted to different agronomic practices to obtain a more representative trial. The descriptive statistic data for the main metabolite content are shown in Table 2.
In Table 2, the concentration ranges, mean, standard deviation, and coefficient of variation (in %, calculated by dividing the standard deviation by the mean of each compound, per 100) of each compound for the groups analyzed can be observed: crocins, picrocrocin, flavonols, and anthocyanins. The results obtained in the RP-HPLC-DAD analysis show that the upper concentration ranges of the main compounds increased compared to the values found in the literature using the same method. For example, for the major crocin, trans-4-GG, the highest concentration value found in the analyses was 33.89 g/kg, which was higher than the levels found in the study of Moratalla et al. [11], where the maximum concentration value found for this crocin was 18.98 g/kg. The same was observed in the case of one of the other major crocins, trans-3-Gg, where the maximum concentration value found was 25.06 g/kg, while the value found previously by Moratalla et al. [11] was 13.51 g/kg.
Since the sample analyzed in this work was fresh flowers and not saffron spice, to compare the metabolite content, it was necessary to perform a conversion, as certain metabolites such as crocins are more concentrated in the spice. To carry out the conversion, the value obtained in the flowers was multiplied by 9.2 [8]. Thus, for example, the value of 33.89 g/kg in trans-4-GG crocin would be equivalent to 315.17 g/kg of the same crocin in the spice. The most relevant kaempferol was k-3-O-β-s, since its concentration was much higher than the rest of the kaempferols. Due to its high concentration, it is the only kaempferol among those analyzed for which proper calibration and validation were achieved. The picrocrocin average was higher than that determined by Kabiri et al. [21], as multiplying the values found in the present study by 9.2 to determine the corresponding content in spice revealed a value of 277 g/kg, higher than the one of 120 g/kg determined by Kabiri et al. [21]. As for the content of k-3-O-β-s, a higher concentration (35.8 g/kg, reaching values of 50.2 g/kg) was determined compared to that determined by Carmona et al. [22] (3.13 g/kg), mainly because the sample in that study was saffron spice and not the whole flower. Since it is in the floral residue where the highest concentration of kaempferol is found, the result is a concentration of k-3-O-β-s more than ten times higher than in the spice. If the maximum value obtained for k-3-O-β-s (50.2 g/kg) is divided by 0.78, the equivalent concentration that the tepals of flowers would have can be determined [23], rendering us able to compare such values with those obtained by Vago et al. [24]. In this way, the higher concentration of k-3-O-s would be 64.3 mg/g in tepals, which would be a value consistent with those determined by Vago et al. [24]. In relation to the coefficients of variation (CV) obtained (Table 2), in most cases, values greater than 30% were observed, which could demonstrate the great variability found in the different flowers.

3.2. NIR Interpretation for Crocus sativus L. Flower Metabolites

The analysis of the main metabolites of the Crocus sativus L. flowers, which reveals the quality of the corms, involves the analysis of either the stigma or the floral residue, having to resort in both cases to chromatographic techniques such as RP-HPLC-DAD [11,17]. This technique, among others, is a time-consuming, labor-intensive method and requires sample treatment and is not easily scaled up to hundreds of samples. Replacing its use by employing the near-infrared (NIR) spectroscopy technique for the estimation of any of the main metabolites presents a better option for routine analysis in time. Samples were fully scanned in the range of 10,000–4000 cm−1 (750–2500 nm) to fine-tune the scan to generate the best PLS-R model.
Figure 2 shows the total spectra obtained by NIR spectroscopy, represented by different colors, and, as can be seen, all the samples show a similar spectrum, with higher or lower absorbance, depending on the regions. In NIR spectroscopy, there are different measurement ranges associated with different types of chemical bonds. Among the distinctive features observed, there are three prominent bands at 8331, 6853, and 5164 cm−1, reflecting the fundamental and overtone vibrations of the O-H group present in water molecules. Another significant band at 5798 cm−1 corresponds to the fundamental vibration of the C-H group, typically indicative of the presence of lipids, oils, and aromatic compounds [25]. It is worth noting the signal centered at 4670 cm−1, frequently associated with C-H groups found in aromatic compounds, especially those with phenolic characteristics [26].

3.3. NIR Calibration and Validation

As seen in Table 2, many compounds were quantified in the flower samples. However, it was not possible to obtain positive calibration–validation results. Thus, the best results were obtained for the following combinations/compounds: (a) Σcis-crocins, which includes cis-4-GG, cis-3-Gg, cis-2-gg, and cis-2-G; (b) Σtrans-crocins, which encompasses trans-5-tG, trans-5-ng, trans-4-GG, trans-4-ng, trans-3-Gg, and trans-2-g; (c) Σ crocins, as the sum of the previous cis and trans; (d) picrocrocin; and (e) kaempferol-3-O-β-sophoroside, which, due to its high concentration among the flavonols analyzed, proper calibration and validation had been achieved. Related to the anthocyanin compound, it was not possible to achieve good calibration for any of those analyzed, nor for all of them.
Table 3 shows the statistics of the calibration set of the metabolites above in all samples analyzed. For the calibration model’s development, the entire infrared spectral region (10,000–4000 cm−1) was considered for spectral acquisition after eliminating the redundant spectra. This was mainly because the target molecules in this study contain several types of bonds, which makes it difficult to select a narrower range than the one chosen for calibration–validation. As shown in Table 2, a wide-ranging concentration value was found in the Crocus sativus L. flowers for each metabolite, indicating good scattering for such model development. In the initial attempt, it was sought to calibrate the different compounds by focusing on the ranges corresponding to the mentioned peaks. However, the calibrations yielded highly deficient results, with excessively low R2 values. This outcome can likely be attributed, together with the above, to the intrinsic nature of the samples, since saffron flower is a heterogeneous sample, containing a myriad of compounds. These compounds may vary in composition and concentration, leading to potential interference in the detection of the target compounds. The presence of such diverse compounds within the saffron flower could introduce complexities and challenges in accurately quantifying the different metabolites. Factors such as the varying chemical properties, solubilities, and reactivities of these compounds may contribute to the observed difficulties in achieving robust calibration. Therefore, further refinement and optimization of the calibration methodology may be necessary to account for the complexities inherent in the saffron flower sample and to improve the accuracy and reliability of calibration.
The more accurate calibration model, for each metabolite or sum of them, obtained with NIR spectral data against their RP-HPLC-DAD determination is summarized in Table 4. Thus, the final model selected and presented was based on the analyses of some error parameters: higher RPD and lower standard errors given by RMSECV (root mean square error of cross-validation) and RMSEP (root mean square error of prediction) (Table 4).
PLS-R was used to perform the calibration model with the more appropriate pretreatments as above to increase the performance of the predictive models in the selected spectral range, which are shown in Table 3, where different spectral ranges were identified for each model depending on the specific metabolite analyzed. Notably, for the major kaempferol in Crocus sativus L. flowers (kaempferol-3-O-β-sophoroside), the entire range of measurements (10,000–4000 cm−1) was considered. However, for crocins and picrocrocin, a narrower spectral range of 9000.9–4000 cm−1 was deemed suitable for calibration and validation.
In the calibration process for picrocrocin, a first derivative with 9 points was applied, followed by smoothing with 5 points. This method resulted in an R2cv value of 82.9%. Similarly, for the sum of crocins, which includes both cis and trans-crocins, the same spectral range was used. Employing a first derivative with 13 points and a smoothing of 13 points, an R2cv of 98.83% was achieved. However, it is worth noting that individual calibrations were conducted for cis and trans-crocins, since it is known that their bioactivity is different [27]; so, for this reason, it was considered necessary to perform a separate calibration of them. In this case, the use of higher derivatives and additional smoothing points to mitigate noise during derivation was requited. Thus, trans-crocins were calibrated within the same spectral range as crocins but with a higher derivative (second with 13 points) and more smoothing (19 points), yielding an R2cv of 84.05%. On the other hand, calibration for cis-crocins was also within the same spectral range but with different smoothing (13 points), resulting in an R2cv of 90.81%. If the value obtained for the calibration of the sum of crocins is compared with other calibration values used in the food industry, it can be observed that the R2 value obtained in the present work is equal to or even higher than some of the R2 values of those calibrations. For example, Kawasaki et al. (2008) [28] obtained a calibration for the fat content of raw milk with a validation R2 of 0.95. Another example is the calibration of the sugar content in pears conducted by Xu, Qi, Sun, Fu, and Ying (2012) [29], which obtained a validation R2 value of 0.87. Finally, for the calibration of kaempferol-3-O-β-sophoroside, a first derivative with 5 points was applied, followed by smoothing with 13 points, resulting in an R2cv of 80.51%.
In previous studies related to NIR calibration in saffron spice, good results to determine the content of the main crocins and picrocrocin was achieved [30]. Since saffron flower is a sample containing many more interferences than the spice and has different analyte concentrations due to being a more complex matrix, the calibrations for the spice were not suitable for flower calibration. In the present research, it was started with the range of 10,000–4000 cm−1 for sample measurement using NIR, achieving R2cv values between 0.82 and 0.98.
R2cv values between 0.6 and 0.8 provide an indication of an approximate predictive model, while values between 0.8 and 0.9 suggest that the model is suitable. However, when the R2c surpasses the threshold of 0.9, it can be said that the model is excellent in terms of its predictive capability [31]. The R2cv value obtained for total crocins’ validation was 0.98 (Table 4), suggesting that the method is excellent in terms of predictive capacity. The validations of the sum of cis-crocins and the sum of trans-crocins resulted in R2cv values of 0.90 and 0.84, respectively, which suggest a good predictive model, being excellent for cis-crocins. Related to picrocrocin, the R2cv value obtained was 0.82, indicating that the predictive capacity of this method is good, and, for kaempferol-3-O-β-sophoroside, an R2cv of 0.80 was obtained, also considered indicative of a correct predictive model.
The Ratio of Performance to Deviation (RPD) index, which reports the ratio between the standard deviation of the reference data of the validation set and the standard error of cross-validation prediction or the test set validation, has been employed as an essential measure to assess the accuracy of the model. This parameter is crucial for determining the applicability of the model. Different threshold values for the accuracy of the model given by RPD have been found in the bibliograph. According to Prieto et al. [32], RPD values exceeding 3 are deemed suitable for screening purposes, while those surpassing 5 are considered appropriate for quality control. Values surpassing 8 are deemed exceptional for all analytical tasks. In the validation of the sum of crocins, RPD values exceeding 8 were attained. Similarly, in the validations of picrocrocin and trans-crocins, both RPD values exceeded 5, with picrocrocin registering a higher RPD value of 7.27. These results imply that both validations are suitable for a range of applications, including analytical purposes and quality control [32]. However, RPD values obtained for the cis-crocins and kaempferol-3-O-β-sophoroside fell below 3, indicating that their predictive capabilities are not as robust as the previous ones. Nonetheless, given their high R2cv values, such models may still serve a useful purpose in preliminary screening processes.
As observed in Figure 3, differences are observed in the peaks at 7000 cm−1 and 5300 cm−1 when overlaying the spectra of three selected samples with different concentrations of total crocins. As it can be observed in the two peaks indicated, there was a difference in the area under the curve which could have been due to the differences in the total crocin amounts present in the three samples.
Figure 4 shows scatter plots of the values of different compounds from saffron flower obtained by RP-HPLC-DAD versus predicted by NIR. Regression between laboratory values and those predicted by NIR resulted in good correlations: 0.98 for Σ crocins, 0.96 for Σ trans-crocins, 0.85 for Σ cis-crocins, 0.95 for picrocrocin, and 0.81 for kaempferol-3-O-β-sophoroside. The differences between results from the different metabolites could not have been due to their concentration in the flower samples, since, although for the crocins the predictive values obtained were good and this can be related to their higher content in the flowers, in the case of the kaempferol-3-O-β-sophoroside, such favorable predictive values were not observed.

4. Conclusions

Good predictive models have been developed in the NIR calibrations of the principal saffron flower metabolites. Excellent predictive models have been obtained for the calibrations of Σ crocins and picrocrocin, with an RPD close to 8, which qualifies the model for quality control purposes. Although the calibrations of Σcis-crocins (cis-4-GG + cis-3-Gg + cis-2-gg + cis-2-G), Σtrans-crocins (trans-5-tG + trans-5-ng + trans-4-GG+trans-4-ng + trans-3-Gg + trans-2-g), and kaempferol-3-O-β-sophoroside did not achieve as high R2cv or RPD values, the models are still suitable for screening purposes.
Measuring parameters such as the sum of total crocins or picrocrocin in saffron flower allows for an indirect measurement of corm quality since the metabolites present in the flower originate from the nutrients stored in the corm the previous year.

Author Contributions

Investigation, J.F.E.-T. and M.E.M.-N.; Visualization, M.E.M.-N., G.L.A., and R.S.-G.; Writing—Original Draft Preparation, J.F.E.-T., G.L.A., and R.S.-G.; Writing—Review and Editing, G.L.A. and R.S.-G.; Methodology, M.E.M.-N., G.L.A., and R.S.-G.; Supervision and Project Administration, G.L.A. and R.S.-G. All authors have read and agreed to the published version of the manuscript.

Funding

The authors wish to thank the Ministry of Education, Universities and Research of the Community Board of Castilla-La Mancha, and the European Regional Development Fund (FEDER) for financing this work through the AZUVOL II project (ref.: SBPLY/21/180501/000014) and the Universidad de Castilla-La Mancha for the financed Project 2023-GRIN-34180.

Data Availability Statement

The data presented in this study are available on request from the corresponding author due to privacy.

Acknowledgments

The authors wish to thank D.O. “Azafrán de La Mancha” (Castilla-La Mancha, Spain) and the farmers involved for their support in this research in the supply of the samples.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. José Bagur, M.; Alonso Salinas, G.L.; Jiménez-Monreal, A.M.; Chaouqi, S.; Llorens, S.; Martínez-Tomé, M.; Alonso, G.L. Saffron: An Old Medicinal Plant and a Potential Novel Functional Food. Molecules 2017, 23, 30. [Google Scholar] [CrossRef] [PubMed]
  2. Valle García-Rodríguez, M.; Serrano-Díaz, J.; Tarantilis, P.A.; López-Córcoles, H.; Carmona, M.; Alonso, G.L. Determination of Saffron Quality by High-Performance Liquid Chromatography. J. Agric. Food Chem. 2014, 62, 8068–8074. [Google Scholar] [CrossRef] [PubMed]
  3. Sánchez, A.M.; Carmona, M.; del Campo, C.P.; Alonso, G.L. Solid-Phase Extraction for Picrocrocin Determination in the Quality Control of Saffron Spice (Crocus sativus L.). Food Chem. 2009, 116, 792–798. [Google Scholar] [CrossRef]
  4. Chryssanthi, D.G.; Lamari, F.N.; Georgakopoulos, C.D.; Cordopatis, P. A New Validated SPE-HPLC Method for Monitoring Crocetin in Human Plasma-Application after Saffron Tea Consumption. J. Pharm. Biomed. Anal. 2011, 55, 563–568. [Google Scholar] [CrossRef] [PubMed]
  5. Okonkwo, D. Trans-Sodium Crocetinate Increases Oxygen Delivery to Brain Parenchyma in Rats on Oxygen Supplementation. Neurosci. Lett. 2003, 352, 97–100. [Google Scholar] [CrossRef] [PubMed]
  6. Vignolini, P.; Heimler, D.; Pinelli, P.; Ieri, F.; Sciullo, A.; Romani, A. Characterization of By-Products of Saffron (Crocus sativus L.) Production. Nat. Prod. Commun. 2008, 3, 1959–1962. [Google Scholar] [CrossRef]
  7. Serrano-Díaz, J.; Sánchez, A.M.; Maggi, L.; Martínez-Tomé, M.; García-Diz, L.; Murcia, M.A.; Alonso, G.L. Increasing the Applications of Crocus sativus Flowers as Natural Antioxidants. J. Food Sci. 2012, 77, C1162–C1168. [Google Scholar] [CrossRef] [PubMed]
  8. Wu, Y.; Gong, Y.; Sun, J.; Zhang, Y.; Luo, Z.; Nishanbaev, S.Z.; Usmanov, D.; Song, X.; Zou, L.; Benito, M.J. Bioactive Components and Biological Activities of Crocus sativus L. Byproducts: A Comprehensive Review. J. Agric. Food Chem. 2023, 71, 19189–19206. [Google Scholar] [CrossRef] [PubMed]
  9. Gismondi, A.; Fanali, F.; Martínez Labarga, J.M.; Caiola, M.G.; Canini, A. Crocus sativus L. Genomics and Different DNA Barcode Applications. Plant Syst. Evol. 2013, 299, 1859–1863. [Google Scholar] [CrossRef]
  10. Eghbali, S.; Farhadi, F.; Askari, V.R. An Overview of Analytical Methods Employed for Quality Assessment of Crocus sativus (Saffron). Food Chem. X 2023, 20, 100992. [Google Scholar] [CrossRef]
  11. Moratalla-López, N.; Sánchez, A.M.; Lorenzo, C.; López-Córcoles, H.; Alonso, G.L. Quality Determination of Crocus sativus L. Flower by High-Performance Liquid Chromatography. J. Food Compos. Anal. 2020, 93, 103613. [Google Scholar] [CrossRef]
  12. Bwambok, D.K.; Siraj, N.; Macchi, S.; Larm, N.E.; Baker, G.A.; Pérez, R.L.; Ayala, C.E.; Walgama, C.; Pollard, D.; Rodriguez, J.D.; et al. Qcm Sensor Arrays, Electroanalytical Techniques and Nir Spectroscopy Coupled to Multivariate Analysis for Quality Assessment of Food Products, Raw Materials, Ingredients and Foodborne Pathogen Detection: Challenges and Breakthroughs†. Sensors 2020, 20, 6982. [Google Scholar] [CrossRef]
  13. Cebi, N.; Bekiroglu, H.; Erarslan, A. Nondestructive Metabolomic Fingerprinting: FTIR, NIR and Raman Spectroscopy in Food Screening. Molecules 2023, 28, 7933. [Google Scholar] [CrossRef]
  14. Tian, W.; Li, Y.; Guzman, C.; Ibba, M.I.; Tilley, M.; Wang, D.; He, Z. Quantification of Food Bioactives by NIR Spectroscopy: Current Insights, Long-Lasting Challenges, and Future Trends. J. Food Compos. Anal. 2023, 124, 105708. [Google Scholar] [CrossRef]
  15. González-Domínguez, R.; Sayago, A.; Fernández-Recamales, Á. An Overview on the Application of Chemometrics Tools in Food Authenticity and Traceability. Foods 2022, 11, 3940. [Google Scholar] [CrossRef]
  16. Wang, H.P.; Chen, P.; Dai, J.W.; Liu, D.; Li, J.Y.; Xu, Y.P.; Chu, X.L. Recent Advances of Chemometric Calibration Methods in Modern Spectroscopy: Algorithms, Strategy, and Related Issues. TrAC-Trends Anal. Chem. 2022, 153, 116648. [Google Scholar] [CrossRef]
  17. Stelluti, S.; Caser, M.; Demasi, S.; Scariot, V. Sustainable Processing of Floral Bio-Residues of Saffron (Crocus sativus L.) for Valuable Biorefinery Products. Plants 2021, 10, 523. [Google Scholar] [CrossRef]
  18. Predieri, S.; Magli, M.; Gatti, E.; Camilli, F.; Vignolini, P.; Romani, A. Chemical Composition and Sensory Evaluation of Saffron. Foods 2021, 10, 2604. [Google Scholar] [CrossRef]
  19. Luo, Z.; Thorp, K.R.; Abdel-Haleem, H. A High-Throughput Quantification of Resin and Rubber Contents in Parthenium Argentatum Using near-Infrared (NIR) Spectroscopy. Plant Methods 2019, 15, 1–14. [Google Scholar] [CrossRef]
  20. Suchat, S.; Pioch, D.; Palu, S.; Tardan, E.; van Loo, E.N.; Davrieux, F. Fast Determination of the Resin and Rubber Content in Parthenium Argentatum Biomass Using near Infrared Spectroscopy. Ind. Crops Prod. 2013, 45, 44–51. [Google Scholar] [CrossRef]
  21. Kabiri, M.; Rezadoost, H.; Ghassempour, A. A Comparative Quality Study of Saffron Constituents through HPLC and HPTLC Methods Followed by Isolation of Crocins and Picrocrocin. LWT 2017, 84, 1–9. [Google Scholar] [CrossRef]
  22. Carmona, M.; Sánchez, A.M.; Ferreres, F.; Zalacain, A.; Tomás-Barberán, F.; Alonso, G.L. Identification of the Flavonoid Fraction in Saffron Spice by LC/DAD/MS/MS: Comparative Study of Samples from Different Geographical Origins. Food Chem. 2007, 100, 445–450. [Google Scholar] [CrossRef]
  23. Serrano-Díaz, J.; Sánchez, A.M.; Martínez-Tomé, M.; Winterhalter, P.; Alonso, G.L. A Contribution to Nutritional Studies on Crocus sativus Flowers and Their Value as Food. J. Food Compos. Anal. 2013, 31, 101–108. [Google Scholar] [CrossRef]
  24. Vago, R.; Trevisani, F.; Vignolini, P.; Vita, C.; Fiorio, F.; Campo, M.; Ieri, F.; Di Marco, F.; Salonia, A.; Romani, A.; et al. Evaluation of Anti-Cancer Potential of Saffron Extracts against Kidney and Bladder Cancer Cells. Food Biosci. 2024, 57, 103501. [Google Scholar] [CrossRef]
  25. Frizon, C.N.T.; Oliveira, G.A.; Perussello, C.A.; Peralta-Zamora, P.G.; Camlofski, A.M.O.; Rossa, Ü.B.; Hoffmann-Ribani, R. Determination of Total Phenolic Compounds in Yerba Mate (Ilex Paraguariensis) Combining near Infrared Spectroscopy (NIR) and Multivariate Analysis. LWT 2015, 60, 795–801. [Google Scholar] [CrossRef]
  26. Li, W.; Qu Haibin, H. Rapid Quantification of Phenolic Acids in Radix Salvia Miltrorrhiza Extract Solutions by FT-NIR Spectroscopy in Transflective Mode. J. Pharm. Biomed. Anal. 2010, 52, 425–431. [Google Scholar] [CrossRef]
  27. Moratalla-López, N.; Bagur, M.J.; Lorenzo, C.; Martínez-Navarro, M.E.; Rosario Salinas, M.; Alonso, G.L. Bioactivity and Bioavailability of the Major Metabolites of Crocus sativus L. Flower. Molecules 2019, 24, 2827. [Google Scholar] [CrossRef]
  28. Kawasaki, M.; Kawamura, S.; Tsukahara, M.; Morita, S.; Komiya, M.; Natsuga, M. Near-Infrared Spectroscopic Sensing System for on-Line Milk Quality Assessment in a Milking Robot. Comput. Electron. Agric. 2008, 63, 22–27. [Google Scholar] [CrossRef]
  29. Xu, H.; Qi, B.; Sun, T.; Fu, X.; Ying, Y. Variable Selection in Visible and Near-Infrared Spectra: Application to on-Line Determination of Sugar Content in Pears. J. Food Eng. 2012, 109, 142–147. [Google Scholar] [CrossRef]
  30. Zalacain, A.; Ordoudi, S.A.; Díaz-Plaza, E.M.; Carmona, M.; Blázquez, I.; Tsimidou, M.Z.; Alonso, G.L. Near-Infrared Spectroscopy in Saffron Quality Control: Determination of Chemical Composition and Geographical Origin. J. Agric. Food Chem. 2005, 53, 9337–9341. [Google Scholar] [CrossRef]
  31. Zornoza, R.; Guerrero, C.; Mataix-Solera, J.; Scow, K.M.; Arcenegui, V.; Mataix-Beneyto, J. Near Infrared Spectroscopy for Determination of Various Physical, Chemical and Biochemical Properties in Mediterranean Soils. Soil. Biol. Biochem. 2008, 40, 1923–1930. [Google Scholar] [CrossRef] [PubMed]
  32. Prieto, N.; Pawluczyk, O.; Dugan, M.E.R.; Aalhus, J.L. A Review of the Principles and Applications of Near-Infrared Spectroscopy to Characterize Meat, Fat, and Meat Products. Appl. Spectrosc. 2017, 71, 1403–1426. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Flowers emerging from a cluster of corms derived from a single corm after three years in the same cultivation soil.
Figure 1. Flowers emerging from a cluster of corms derived from a single corm after three years in the same cultivation soil.
Horticulturae 10 00593 g001
Figure 2. Near-infrared spectrum of Crocus sativus flower samples.
Figure 2. Near-infrared spectrum of Crocus sativus flower samples.
Horticulturae 10 00593 g002
Figure 3. First derivative of 3 spectra with different crocin contents (high, medium, and low). The rounded areas correspond to the peaks with the highest signal.
Figure 3. First derivative of 3 spectra with different crocin contents (high, medium, and low). The rounded areas correspond to the peaks with the highest signal.
Horticulturae 10 00593 g003
Figure 4. Scatter plot obtained with RP-HPLC-DAD versus values predicted by the NIR in saffron flower samples using PLS-R: (a) Σ cis-crocins; (b) Σ trans-crocins; (c) Σ crocins; (d) picrocrocin; (e) kaempferol-3-O-β-sophoroside.
Figure 4. Scatter plot obtained with RP-HPLC-DAD versus values predicted by the NIR in saffron flower samples using PLS-R: (a) Σ cis-crocins; (b) Σ trans-crocins; (c) Σ crocins; (d) picrocrocin; (e) kaempferol-3-O-β-sophoroside.
Horticulturae 10 00593 g004aHorticulturae 10 00593 g004b
Table 1. Chromatographic conditions of RP-HPLC-DAD.
Table 1. Chromatographic conditions of RP-HPLC-DAD.
MinutesAcetonitrile (%)Water/Trifluoroacetic Acid (%)
00100
303070
40–458020
50–552080
55–590100
Table 2. Descriptive statistics for Crocus sativus L. flower metabolites (RP-HPLC-DAD data).
Table 2. Descriptive statistics for Crocus sativus L. flower metabolites (RP-HPLC-DAD data).
CompoundRange (g kg−1 Flower Dry Weight)MeanSDCV (%)
cis-4-GG0.56–2.301.310.5844
cis-2-G0.03–0.520.260.0915
cis-3-Gg0.58–2.282.430.7932
cis-2-gg0.03–0.070.030.0133
trans-5-tG0.15–0.510.280.1316
trans-5-ng0.22–0.350.240.1250
trans-4-GG8.32–33.8916.848.5450
trans-4-ng0.17–2.300.200.0840
trans-3-Gg3.94–25.069.564.4746
trans-2-G 0.03–0.380.520.3159
Picrocrocin2.03–49.530.146.4421
K-3-O-β-s0.42–50.235.287.6521
K-3-O-s-7-O-g3.42–7.405.581.1520
K0.01–0.190.090.0666
D-3,5-di-O-g1.09–9.16.211.0416
P-3,5-di-O-g0.27–0.450.30.1136
D-3-O-g0.76–2.231.080.2725
M-3,5-di-O-g0.07–4.642.671.2446
P-3-O-g0.24–0.260.210.0942
cis-4-GG: cis-crocetin di-(β-D-gentiobiosyl) ester; cis-2-G: cis-crocetin (β-D-gentiobiosyl) ester); cis-3-Gg: cis-crocetin (β-D-glucosyl)-(β-D-gentiobiosyl) ester; cis-2-gg: cis-crocetin di-(β-D-glucosyl) ester; trans-5-tG: trans-crocetin (β-D-triglucosyl)-(β-D-gentiobiosyl) ester; trans-5-ng: trans-crocetin (β-D-neapolitanosyl)-(β-D-; trans-4-GG: trans-crocetin di-(β-D-gentiobiosyl) ester; trans-4-ng: trans-crocetin (β-D-neapolitanosyl)-(β-D-glucosyl) ester; trans-3-Gg: trans-crocetin (β-D-glucosyl)-(β-D-gentiobiosyl) ester; trans-2-G; trans-crocetin (β-D-glucosyl)-(β-D-gentiobiosyl) ester; K-3-β-O-s: kaempferol-3-O-β-sophoroside; K-3-O-s-7-O-g: kaempferol-3-O-β-sophoroside-7-O-β-glucoside; K: kaempferol aglycone; D-3,5-di-O-g: delphinidin-3,5-di-O-β-glucoside; P-3,5-di-O-g: petunidin-3,5-di-O-β-glucoside; D-3-O-g: delphinidin-3-O-β-glucoside; M-3,5-di-O-g: malvidin-3,5-diO-β-glucoside; P-3-O-g: petunidin-3-O-β-glucoside.
Table 3. NIR statistical parameters for the calibration data of the measured Crocus sativus flower metabolites.
Table 3. NIR statistical parameters for the calibration data of the measured Crocus sativus flower metabolites.
NSpectral Range
(cm−1)
Pre-Process *R2c
Σ  cis-Crocins1739000.9–40002.13.131.00
Σ  trans-Crocins1739000.9–40002.13.190.97
Σ Crocins1739000.9–40001.13.130.99
Picrocrocin779000.9–40001.9.50.95
K-3-O-β-s7710,000–40001.5.130.99
* The three-digit code (a.b.c.), where ‘a’ refers to the number of the derivative (1 and 2), ‘b’ is the interval over which the derivative is calculated (5, 9, and 13), and ‘c’ corresponds to the number of data points in a running mean or smoothing (5, 13, and 19). R2c: coefficient of determination in calibration; RMSEC: standard error of calibration. Σ cis-Crocins: cis-4-GG + cis-3-Gg + cis-2-gg + cis-2-G; Σ trans-Crocins: trans-5-tG + trans-5-ng + trans-4-GG + trans-4-ng + trans-3-Gg + trans-2-g. K-3-β-O-s: kaempferol-3-O-β-sophoroside.
Table 4. NIR statistical parameters for the cross-validation data of the measured Crocus sativus L. flower metabolites.
Table 4. NIR statistical parameters for the cross-validation data of the measured Crocus sativus L. flower metabolites.
CompoundR2cvRMSECVRMSECRPDRPDcRPDp
Σ  cis-Crocins0.900.480.112.472.311.25
Σ  trans-Crocins0.842.214.016.512.451.16
Σ Crocins0.981.240.788.838.835.45
Picrocrocin0.823.15.27.29.237.27
K-3-O-β-s0.803.460.0012.212.181.22
R2CV: coefficient of determination in cross validation RMSECV: root mean squared error of cross validation; RMSEC: root mean squared error of calibration; RPD: calculated according to Luo et al. [19]; RPDc: residual predictive deviation (SD/SECV); RPDp: SDcal/SEP. Σ cis-Crocins: cis-4-GG + cis-3-Gg + cis-2-gg + cis-2-G; Σ trans-Crocins: trans-5-tG + trans-5-ng + trans-4-GG + trans-4-ng + trans-3-Gg + trans-2-g. K-3-β-O-s: kaempferol-3-O-β-sophoroside.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Escobar-Talavera, J.F.; Martínez-Navarro, M.E.; Alonso, G.L.; Sánchez-Gómez, R. Determination of Saffron Flower Metabolites by Near-Infrared Spectroscopy for Quality Control. Horticulturae 2024, 10, 593. https://doi.org/10.3390/horticulturae10060593

AMA Style

Escobar-Talavera JF, Martínez-Navarro ME, Alonso GL, Sánchez-Gómez R. Determination of Saffron Flower Metabolites by Near-Infrared Spectroscopy for Quality Control. Horticulturae. 2024; 10(6):593. https://doi.org/10.3390/horticulturae10060593

Chicago/Turabian Style

Escobar-Talavera, Jorge F., María Esther Martínez-Navarro, Gonzalo L. Alonso, and Rosario Sánchez-Gómez. 2024. "Determination of Saffron Flower Metabolites by Near-Infrared Spectroscopy for Quality Control" Horticulturae 10, no. 6: 593. https://doi.org/10.3390/horticulturae10060593

APA Style

Escobar-Talavera, J. F., Martínez-Navarro, M. E., Alonso, G. L., & Sánchez-Gómez, R. (2024). Determination of Saffron Flower Metabolites by Near-Infrared Spectroscopy for Quality Control. Horticulturae, 10(6), 593. https://doi.org/10.3390/horticulturae10060593

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop