Infrared guided smart food formulation: an innovative spectral reconstruction strategy to develop anticipated and constant apple puree products

Apple puree accounts for the second largest market of fruit puree, with a global market value of about 2,000 million USD annually^[1]. It is used as the basic ingredient of jams, preserves or compotes, which are popular among people of all ages, especially for babies and elders^[2]. According to previous work, a large diversity of apple cultivars^[3] and processing conditions (cooking parameters, grinding intensity and refining levels, etc.)^[4−6] can introduce strong chemical, textural and rheological variations on processed purees. However, the current puree processing systems are not always adapted and optimized to the raw apples but meet the required constant quality standards of final products. Further, the increasing demands of various anticipated products for current consumers have put much stress on the development of personalized puree products from industrial manufacturers. Therefore, it would be highly beneficial to develop innovative puree production strategies, that can consider a large variability of raw materials and provide new solutions to reach the anticipated and constant taste and texture of purees.

Puree formulation is one of the most economical and efficient strategies for manufacturers to adjust the texture and taste of final puree products depending on a mixture of different proportions of single apple varieties^[7]. However, fruit manufacturers do not often have the access to good choices to determine how to formulate puree products with large variability of raw materials. The challenge is therefore to develop innovative strategies to provide specific guidance for the formulation of final puree products based on the information of single cultivar puree, to reach their anticipated and constant taste and texture.

Visible and near-infrared (Vis-NIR) spectroscopy has been applied as a simple and rapid technique to give considerable predictions of chemical, physical and textural properties of raw apples^[8−10] and processed puree^[5,11−13]. In our previous work, a spectral reconstruction strategy based on the concentration profile of multivariate curve resolution-alternative least squares (MCR-ALS) was first developed^[13]. Based on that, the mid-infrared (MIR) spectra of single-cultivar apple purees can be used to reconstruct the spectra of differently formulated purees, then the multivariate regression models using the reconstructed spectra of formulated purees can successfully predict their quality characteristics (soluble sugars, titratable acidity, pH and viscosity, etc.)^[13]. However, this strategy was not available to reconstruct the Vis-NIR spectra of formulated purees from their corresponding spectra of single-cultivar purees, because of the unacceptable concertation profiles of MCR-ALS. One possible reason could be the large sets of spectral variables that were strongly collinear and noisy, which may affect the success of spectral reconstruction.

To partly compensate for these effects, spectral variable selection is an important approach, as it tends to parsimonious data representation and can result in multivariate models with greater predictive ability^[14]. Particularly, several spectral variable selection methods, such as competitive adaptive reweighted sampling (CARS)^[15], successive projections algorithm (SPA)^[16], and uninformative variable elimination (UVE)^[17], can significantly improve the Vis-NIR prediction accuracy. The challenging work here was to reconstruct the spectra of final formulated purees according to the most relative Vis-NIR spectral variables of single cultivar apple purees by MCR-ALS. Then, the predictive models of formulated puree quality traits (physical and chemical) using the reconstructed spectra dataset could predict the properties of formulated puree products based on the relative Vis-NIR spectra of composed single-cultivar purees. If so, this new strategy opens a new possibility to provide practical and suitable strategies for the multicriteria optimization of puree formulation with anticipated and constant quality.

This study intended to develop a smart food formulation model based on the rapid and high-throughput spectral information of their individual composed components, which could provide various formulation guidance and monitor their quality parameters of final products. To reach this objective, Vis-NIR technique coupled with spectral variable selection methods (CARS, SPA, UVE) were applied both on the different formulated purees and their corresponding single-cultivar purees to highlight their featured spectral variables; then the selected spectra variables of single-cultivar purees to reconstruct spectra of formulated purees; and finally investigated the possibility to develop regression models to evaluate the quality parameters of final formulated purees.

The experiment was conducted on four apple varieties: 'Golden Delicious'(GD), 'Granny Smith'(GS), 'Braeburn'(BR), and 'Royal Gala'(GA). They were harvested at commercial maturity from La Pugère experimental orchard (Mallemort, Bouches du Rhône, France) in 2019, and stored at 4 °C and around 90% relative humidity for up to 2 months to ensure starch regression. A multi-functional processing system (Roboqbo, Qb8-3, Bentivoglio, Italy) was used to process apple purees following a Hot Break recipe: cooked at 95 °C for 5 min at a 1,500 rpm grinding speed, then cooled down to 65 °C while maintaining the grinding speed. During three successive weeks, around 2 kg of each apple cultivar was processed into single-cultivar purees, then conditioned in two hermetically sealed cans: one was cooled in a cold room (4 °C) before formulation, while the other was stored at –20 °C for biochemical measurement of individual sugars (fructose, sucrose, and glucose) and malic acid.

Puree formulation

After puree processing, the four single cultivar purees were formulated by two of each of them, into six different experimental groups named 'A' (GS × GS), 'B' (GD × BR), 'C' (GD × GA), 'D' (GS × BR), 'E' (GS × GA) and 'F' (BR × GA), respectively (Fig. 1). Each experimental group (A−F) included nine samples with different formulated proportions of weight, which were divided into two subsets: the first included six proportions (10%:90%, 25%:75%, 50%:50%, 75%:25%, 90%:10%, 95%:5%) for the modeling set, while the second included three proportions (80%:20%, 33%:67%, 14%:86%) for the external prediction set. Finally, all the single and formulated purees were prepared for both the Vis-NIR spectral measurements and quality characterizations.

Figure 1. 

Experimental scheme of puree reformation, quality characterizations, and spectral acquisition.

Determination of puree quality traits

Physical characterizations

The puree color was determined three times using a CR-400 chromameter (Minolta, Osaka, Japan) and expressed in the CIE 1976 L*a*b* color space (illuminant D65, 0° view angle, illumination area diameter 8 mm). Puree rheological measurements were carried out using a Physica MCR-301 controlled stress rheometer (Anton Paar, Graz, Austria) and a 6-vane geometry (FL100/6W) with a gap of 3.46 mm, at 22.5 °C. The flow curves were performed after a pre-shearing period of 1 min at a shear rate of 50 s⁻¹, followed by 5 min at rest. The viscosity was then measured at a controlled shear rate range of [10; 250] s⁻¹ on a logarithmic ramp. The values of viscosity at 50 and 100 s⁻¹ (η₅₀ and η₁₀₀ respectively) were kept as final indicators of the puree viscosity linked to sensory characteristics during consumption^[18].

Biochemical characterizations

The dry matter content (DMC) was estimated from the weight of freeze-dried samples upon reaching a constant weight (freeze-drier, 5 d). Titratable acidity (TA) was determined by titration up to pH 8.1 with 0.1 mol/L NaOH and expressed in mmol H⁺ kg⁻¹ of fresh weight (FW) using an autotitrator (Methrom, Herisau, Switzerland). Individual sugars and malic acid were quantified using colorimetric enzymatic kits (R-biopharm, Darmstadt, Germany), respectively. The contents of glucose, fructose, sucrose, and malic acid were expressed in g kg⁻¹ FW. The total sugar content (TSC) of each puree sample was presented as the sum of characterized glucose, sucrose and fructose. The individual sugars (fructose, glucose, sucrose) and malic acid contents of formulated puree samples were calculated based on the measured values of processed single cultivar purees.

Spectrum acquisition

A multi-purpose analyzer spectrometer (Bruker Optics®, Wissembourg, France) with OPUS software Version 5.0 (Bruker Optics®) was used to acquire the Vis-NIR spectral data of purees at 23 °C, which can provide diffuse reflectance measurements with a spectral resolution of 8 cm^-1 from 500 to 2,500 nm. Totally 32 scans were recorded and averaged for each spectrum of purees. Purees were transferred into 10 mL glass vials (5 cm height × 18 mm diameter) which were placed on the automated sample wheel of the spectrophotometer. A reference background measurement was automatically activated before each data set acquisition using an internal Spectral on reference. Each puree sample was measured three times on different aliquots. Totally, 36 spectra of single-cultivar purees (4 cultivars × 3 testing weeks × 3 replicates) and 486 spectra of their formulated purees (6 experimental groups × 9 formulated proportions × 3 testing weeks × 3 replicates) (Fig. 1).

Statistical analyses and chemometrics

Statistical analyses of reference data

After checking the normal distribution with a Shapiro-Wilk test (α = 0.05), the reference data of processed purees were presented as mean values and standard deviation values (SD) in Tables 1 & 2. Analysis of variance (ANOVA) and Pearson correlation analysis were carried out to determine the significant differences and internal correlations of puree quality traits of the different single apple cultivars and formulated puree groups using XLSTAT (version 2018.5.52037, Addinsoft SARL, Paris, France) data analysis toolbox, which was described in our previous work^[13]. The physical (a* and b* values), rheological (viscosity η₅₀), and biochemical (SSC, TA, and fructose) parameters of formulated purees were displayed in the boxplots using R software (version 2.6.2)^[19] (Fig. 2).

Figure 2. 

Boxplot of colors (a* and b*), rheological parameters (η₅₀), total sugars (TSC), titratable acidity (TA) and fructose of different formulated puree groups.

Spectral analyses

Spectral discrimination and multivariate analyses were performed with MATLAB 7.5 (Mathworks Inc. Natick, MA, USA) software using the SAISIR package^[20]. Principal component analysis (PCA) was carried out on the single-cultivar puree spectra to evaluate their variability and point out the contributed wavelengths. ANOVA was performed on the Vis-NIR spectra of all formulated purees to analyze their variations during puree formulation.

Spectra variable selection

Three spectral variable selection methods, including competitive adaptive reweighted sampling (CARS), successive projections algorithm (SPA) and uninformative variable elimination (UVE) have been applied respectively on the formulated puree spectra matrix D (n × λ), which was made up with the number of samples (n) and the intensity at each wavelength (λ with 2722 spectral variables from 400 to 2,500 nm) (Fig. 3), to extract the most informative wavelengths for prediction models of each puree quality traits. After comparison, the specific spectral variables related to all puree quality traits were extracted by CARS and composed as the matrix D' (n × λ′) of all formulated purees, consisting of the same number of samples (n) and the intensity at selected variables (λ′ with 671 spectral variables from 400 to 2,500 nm).

Figure 3. 

Process of Vis-NIR spectral data by multivariate resolution alternative least square (MCR-ALS) and spectral reconstruction of formulated purees.

Spectra reconstruction

Multivariate curve resolution-alternative least square (MCR-ALS) is an effective multivariate self-modeling curve resolution method, which can simultaneously elucidate the pure spectra of different species present in processed products and their concentration profiles^[21]. As displayed in Fig. 3, the Vis-NIR spectra of formulated purees were reconstructed from their composed single-cultivar purees based on our previously developed method^[13]. The S^T matrix (s × λ′) is the spectroscopic matrix describing the selected spectral variables (λ′) of all single-cultivar purees (s). The matrix D′ can be mathematically decomposed into the individual contributions related to the spectral information of matrix S^T according to Eqn (1) and is interactively transformed using an alternative least square (ALS) procedure as Eqn (2).

Matrix C (n× q) is the concentration matrix describing the contribution of every single-cultivar puree (q) in reconstructed purees (n). E is the error matrix that provides the data variation not explained by their contributions. The matrix (S^T)⁺ is the pseudo-inverse matrix of S^T. A general constraint used in the curve resolution method is the non-negativity on the concentration profiles.

Once the concentration profiles (matrix C) for each single-cultivar spectrum of Golden Delicious (C_GD), Granny Smith (C_GS), Braeburn (C_BR), and Royal Gala (C_GA), were obtained, they were used to reconstruct a new spectroscopic matrix R (n × k) for monitoring all formulated purees. Each row R_i (i = 1,…n) was made up of a reconstructed spectrum. And each column R_j (j = 1,…λ′) gave the reconstructed spectral intensity at each selected Vis-NIR wavelength based on the corresponding pure puree spectra of Golden Delicious (λ′_GD), Granny Smith (λ′_GS), Braeburn (λ′_BR) and Royal Gala (λ′_GA), following Eqn (3).

Spectral multivariate regression

Spectral pre-processing and multivariate regression were performed with MATLAB 7.5 (Mathworks Inc. Natick, MA, USA) software with the 'PLS' toolbox and displayed in Fig. 4. For all spectral datasets, standard normal variate (SNV) and derivative transform calculation (Savitzky–Golay method, window size = 11, 21, 31, 41) of the first or second order were compared before multivariate regression. SNV pre-processing applied on the Vis-NIR spectral data showed the best performances to predict puree quality and was then systematically used.

Figure 4. 

Overview of the applied methodology of Vis-NIR spectra pre-processing, spectral variable selection, spectral reconstruction and multivariate regression.

The partial least square (PLS) regression models were developed to predict the quality characteristics of formulated purees based on: i) their full spectral variables (FULL) and the selected spectral matrices from CARS, UVE, and SPA, respectively (Table 1); and ii) the reconstructed Vis-NIR spectral matrix (Table 2). All aforementioned spectral matrices correspond to the same reference dataset. The 324 spectra of the formulated purees (6 groups × 6 proportions × 3 weeks × 3 replicates) were used for modeling calibration. Then, the calibrated models were further validated with the external prediction set of 162 puree spectra (6 groups × 3 proportions × 3 weeks × 3 replicates) (see Fig. 1). The optimal numbers of latent variables (LVs) for PLSR models were selected by the Venetian blinds cross-validation method. The prediction ability of developed models was described by the determination coefficients of calibration ($\rm R_c^2$) and external prediction ($\rm R_p^2$), root mean square error of calibration (RMSEC) and prediction (RMSEP), residual predictive deviation (RPD) value as described by Nicolai et al.^[22]. And the most contributed spectral variables (VIP) for each puree characteristic during the modeling training were obtained and analyzed.

After processing, four different single cultivar purees presented significant (p < 0.05) differences in physical, chemical and rheological properties, which were in line with our previous results^[23] (Fig. 2). Particularly, the BR purees had the highest redness (a* values = 1.76 ± 0.27) and lowest yellowness (b* values = 10.18 ± 0.45) among these four puree groups. Four apple cultivars presented a significant difference (p < 0.05) in puree viscosity (η₅₀ and η₁₀₀). Particularly, the GD, GS and BR presented a significantly higher (p < 0.0001) viscosity than GA purees, which were addressed as their bigger particle sizes and promoted cell adhesion with more branched pectin^[24,25]. The TSC and DMC of GD purees were significant (p < 0.0001) higher than the other three groups, while the GS purees gave the highest acidic values (TA = 89.1 ± 1.3 meq/kg FW, and malic acid = 8.1 ± 0.5 g/kg FW). Besides, GD and GA purees presented the significant largest differences (p < 0.0001) in all individual sugar content (fructose, glucose and sucrose).

After puree formulation, three groups with the admixture of BR purees (B, D, F) can introduce a more intensive variation of color parameters (a* and b* values) than other groups (A, C, E) (Fig. 2a & b). Besides, the formulated purees prepared with GS purees (A, D, E) presented a relatively large variability of viscosity at a shear rate of 50 s⁻¹ (η₅₀) (Fig. 2c) than other groups (B, C, F), which is commonly used to describe the in-mouth texture perception of fluid foods^[26]. The variability of TSC in different puree groups was ranked as C > E > F > A > B > D, and the addition of GA purees (C, E, F) resulted in a relatively large range of TSC (Fig. 2). However, the variability of fructose in formulated purees, which is the dominant individual sugar in purees^[27], was different from the results of TSC and DMC. Three puree groups formulated with GS (C, E, F) gave a relatively larger variation of fructose than other groups (A, B, D) (Fig. 2f). The formulation of GS and BR purees (E) resulted in the largest variations of TA values among the six groups, but the limited changes of TA in GD and BR formulated purees (B) (Fig. 2e).

Consequently, the different formulation strategies based on these four different single cultivar purees can provide a large variability of appearance, and rheological and chemical properties.

Spectra analysis of formulated purees and their composed single cultivar purees

PCA was performed on the Vis-NIR spectra of four single cultivar purees to highlight their most variable wavelength regions (Fig. 5). The first principal component (PC1) discriminated against 'GD' and 'BR' purees and accounted for 54.2% of the total variation (Fig. 5a). The specific wavelengths at 524–528 nm, 672 nm in the Vis range and around 800–1,250 nm in the NIR range were the main contributors to the PC1 (Fig. 5b). The two bands at 528 and 672 nm are explained by the anthocyanin and chlorophyll contents of apples, respectively^[28,29]. The NIR region at 800–1,250 nm is known as the absorption of apple carbohydrates and water content^[8,10] and is already used for apple cultivar discrimination^[30]. The 'GS' and 'GA' purees could be separated by the second principal component (PC2), with an explanation rate of 33.4%. The PC2 score was highly correlated to the visible spectra wavelengths at 612 and 672 nm, indicating the large differences in greenness between 'GS' and 'GA' purees^[31].

Figure 5. 

Principal component analysis (PCA) of the SNV pre-treated Vis-NIR spectra of four single cultivar purees. (a) Discrimination map of four apple cultivars. (b) Beta-coefficients of first principal component (PC1). (c) Beta-coefficients of second principal component (PC2).

ANOVA was performed on the Vis-NIR spectra of all formulated apple purees respectively (Fig. 6a & b), to point out the wavelength that varied at different formulation strategies. According to the F-values, the variability of formulation was much higher at the visible spectral region (F-values from 7.1 to 147.6) (Fig. 6a) than at the NIR region (F-values from 2.4 to 15.3) (Fig. 6b). This result demonstrated that the color varied more intensively than the chemical parameters among different formulated apple purees, from the spectroscopic point of view. Particularly, the most variable visible spectral wavelengths were located at 528, 614, and 672 nm, which is similar to the PCA results of four single-cultivar purees (Fig. 5). At the NIR region, the most informative wavelengths were located at 800−1,400 nm, in particular with 916, 1,070, and 1,270 nm, respectively. The NIR region at around 916 nm was related to the C-H and O-H bands of sucrose^[32], and the typical peaks at 1,070 and 1,075 contributed to the soluble sugars of apples^[33] and processed juices^[34].

Figure 6. 

ANOVA results of the SNV pre-treated (a) visible (500–780 nm) spectra and (b) NIR (780–2,500 nm) spectra of all formulated apple purees.

Accordingly, the specific Vis-NIR spectral wavelengths at 524−528 nm, 614, 672, 916, 1,070 and 1,270 nm were shown to be potentially linked to the different chemical and physical variations among different single-cultivar purees and their formulated products.

Predicting characteristics of formulated purees based on Vis-NIR spectra

In this part, PLS regression coupled with all the spectral variables (FULL) or the selected spectral variables based on SPA, CARS and UVE were applied to compare their ability to predict color, rheological and biochemical characteristics of formulated purees (Table 1). As expected, the decreases of determination coefficients between the calibration set ($\rm R_c^2$) and the external validation set ($\rm R_p^2$) for all developed PLS models. According to the RPD values over 2.5^[22], acceptable to good predictions were obtained for a*, viscosity (η₅₀), TA, pH, glucose and malic acids.

For all the color parameters (L*, a* and b*), PLS models coupled with three different spectral variable selection methods gave better predictions than those with full spectral variables. The excellent predictions of a* values were obtained for all PLS models, with an RPD decreasing order based on CARS (RPD = 5.56), UVE (RPD = 5.38) and SPA (RPD = 5.24) methods, respectively. And the selected spectral variables for a* values of purees were mainly dominated in the visible spectral region at around 672 and 614 nm, which were described previously in this study.

Apparent puree viscosity at a share rate value of 50 s⁻¹ (η₅₀) can be considerably predicted by UVE-PLS models, with an $\rm R_p^2$ of 0.87 and an RPD value of 2.73. Besides, the UVE-PLS predictions of puree viscosity at a single shear rate value of 100 s⁻¹ (η₁₀₀) ($\rm R_p^2$ = 0.85 and RPD = 2.49) observed here was much better than its perdition from the full spectral variables ($\rm R_p^2$ = 0.79 and RPD = 2.17). The specific spectral variables of puree viscosity were identified at 672−677 nm, 1,043−1,079 nm and 1,379−1,381 nm, mainly presenting the chlorophyll^[29], the C-H and O-H vibration of soluble (mainly sugars and acids)^[35] and water contents^[36]. These results evidenced the spectral variable selection on Vis-NIR spectra gives obvious improvements to estimate the puree viscosity at different formulation strategies. Different from the poor predictions of puree viscosity using the full spectral variables of NIR^[27] and Vis-NIR^[13], the possible reasons for these improvements of prediction results could be: i) the exclusion of non-relevant or unimportant spectral variables which disturbed the calibration modeling, and ii) the good collinear relationships (r² = 0.72) between puree color a* values and viscosity in different formulated purees, which also has been observed between raw and processed apples from our previous work^[23].

For biochemical characteristics, all the developed PLS models did not give satisfactory predictions for the DMC ($\rm R_P^2$ < 0.66, RPD < 1.66) of the formulated purees. These poor results were relatively lower than our previous predictions by Vis-NIR of SSC ($\rm R_P^2$ = 0.80, RPD = 2.30) and DMC ($\rm R_P^2$ = 0.73, RPD = 1.9) of individually processed purees^[23]. The main reason was probably related to the lower variations of DMC among these formulated purees (SD = 0.006 g/g) in comparison with the previous study taking into account the variability of individual apples (SD = 0.008 g/g). TSC, presenting all the individual sugars (sucrose, fructose and glucose), was satisfactorily predicted by CARS-PLS with 101 selected spectral variables ($\rm R_P^2$ = 0.92, RPD = 3.66). However, the Vis-NIR prediction was only acceptable for glucose (RPD > 2.42) but not for fructose (RPD < 1.61) and sucrose (RPD < 1.60).

Concerning the different characteristics of puree acidity including pH, TA and malic acid, both CARS-PLS and UVE-PLS models provided considerable predictions with $\rm R_p^2$ > 0.87 and RPD > 2.69. It can be noticed that better puree acidity predictions can be obtained using these two methods than the full spectral variables, which was also highlighted on apples^[15]. Particularly, the most contributed spectral variables were located at 664−668 nm and 1,065−1,165 nm, which presented the vibration of O-H from carboxylic acids in fruits^[35].

Consequently, PLS models coupled with three spectral variable selection methods offered better predictions of puree physical, rheological and biochemical characteristics than using full Vis-NIR spectral variables alone. In comparison with three spectral variable selection methods, CARS always extracted the intermediate numbers of spectral variables (higher than SPA and much lower than UVE), but reached similar PLS prediction results as the UVE method with thousands of featured variables. In summary, PLS models coupled with selected VIS-NIR spectral variables have promising to well estimate the* color parameter, viscosity (η₅₀), TSC, TA, pH, glucose and malic acids in formulated apple purees.

Predicting characteristics of formulated purees based on the reconstructed VIS-NIR spectra

Based on the CARS-PLS prediction results, a total 671 Vis-NIR spectral variables were extracted from the 2,722 spectral variables of full wavelengths. According to our method described, MCR-ALS was applied to the selected Vis-NIR spectral variables of all formulated purees and of the four single-cultivar purees, to compute the concentration profiles of relevant single-cultivar composition. In total, 486 spectra of formulated purees were reconstructed (including 324 spectra of the calibration set and 162 spectra of the validation set, respectively) based on their corresponding 36 spectra of single-cultivar purees.

PLS models were developed using these reconstructed Vis-NIR spectra and the reference data of all formulated purees. A good prediction of a* values of puree can be obtained with $\rm R_p^2$ and RPD of 0.92 and 3.30, but not for L* ($\rm R_p^2$ = 0.62) and b* values ($\rm R_p^2$ = 0.42). However, the PLS models based on the reconstructed Vis-NIR spectra were not precise enough to evaluate puree rheological parameters (η50 and η100) ($\rm R_p^2$ < 0.82, RPD < 2.22), according to their PRD values less than 2.5^[22]. This result indicated the difficulties in predicting the viscosity of formulated puree from their Vis-NIR spectra of single-cultivar purees. Besides, the unsatisfactory predictions of DMC of formulated purees ($\rm R_p^2$ = 0.46, RPD = 1.38), were in line with their PLS results based on the real spectra (in Table 1). Although the PLS results of TSC based on the reconstructed Vis-NIR spectra (in Table 2) were less accurate than the real spectral prediction (in Table 1), the PLS models had an acceptable ability to estimate it for all formulated purees (RMSEP = 4.8 g/kg FW, RPD = 2.64). Considering the global acidity parameters, although the slightly lower RMSEP and RPD values than the results obtained directly on the real Vis-NIR spectra (RMSEP = 0.5 meq/kg, RPD = 3.37 for TA, and RMSEP = 0.09, RPD = 2.73 for pH), the PLS models had an acceptable ability to estimate them for all formulated purees (RMESP = 0.6 meq/kg, RPD = 2.55 for TA, and RMSEP = 0.10, RPD = 2.47 for pH). Particularly, similar VIP wavelengths were observed in reconstructed spectra and directly on formulated purees described above, mainly the 664−668 nm and 1,065−1,165 nm, which were described earlier. For individual sugars and acids, PLS models provided the acceptable predictions of malic acid ($\rm R_p^2$ = 0.86, RPD < 2.67), but not for all individual sugars ($\rm R_p^2$ < 0.82, RPD < 2.25).

The concentration profile of MCR-ALS opened a potential way to directly estimate the a* value, TA, pH, malic acid, and glucose for formulated purees based only on the CARS selected spectral variables of single-cultivar purees. Overall, the PLS results of these aforementioned quality parameters based on the reconstructed spectra of formulated purees presented a relatively lower prediction accuracy than directly on real puree spectra, because of the non-negativity of the concentration profiles which could constrain the spectral reconstruction^[37]. Compared to our previous prediction models obtained on the mid-infrared spectra of purees^[3], these results provided further evidence of our new chemometric strategy that using the concentration profile of MCR-ALS to reconstruct the Vis-NIR spectra of formulated purees from their corresponding spectra of single-cultivar purees. Besides, what stands out in this work is the necessity of spectral variable selection on the puree spectra to obtain considerable concentration profiles for the spectral reconstruction of formulated purees.

It's the first time to reconstruct the spectra of formulated products using the concentration profile of MCR-ALS from selected Vis-NIR spectral variables of composed raw materials. Innovatively, this strategy opens the possibility to guide the production of constant and anticipated purees by simply scanning the single-cultivar apple purees in the apple industry. For instance, after acquiring Vis-NIR spectra of the four single-cultivar purees, our developed PLS models could: i) provide several strategies to formulate purees with defined tastes (e.g. 127.3 ± 5.7 g/kg FW of TSC and 7.8 ± 0.2 meq/kg FW of TA, which might be reached with the formulate solutions as 80% GS–20% GA, 33.3% GD–66.6% GS and 80% GS–20% BR purees), depending the might be used in industry; or ii) simulate and optimize puree formulation for anticipated products development depending on the market, such as 20% GD−80% GS purees (low sweetness, high acidity), 33.3% GD–66.6% GA purees (high sweetness, low acidity) and GS 14%–GA 86% (low sweetness, low acidity).

This study firstly demonstrated a better ability of Vis-NIR spectroscopy coupled with advanced chemometric methods (CARS, UVE and SPA variable selection and PLS regression) to estimate the physical (a* value), rheological (η₅₀ and η₁₀₀) and chemical compositions of apple purees than the use of full wavelengths.

Further, an innovative spectral reconstruction strategy based on the MCR-ALS and spectral variable selection was developed to provide practical and suitable strategies for the multicriteria optimization of puree formulation with anticipated and constant quality (a* value, TSC, TA, glucose, malic acid) from their composed single-cultivar purees. As far as we know, this was the first report providing the potential formulation strategy to develop anticipated and constant final fruit products using the Vis-NIR spectral information of the initial purees based on the spectral reconstruction approach. Further this new chemometric strategy has the potential to provide production guidance for other food formulations, such as multifruit juices, blend oils, and even admixed flavoring agent etc, based on the Vis-NIR spectra acquired directly on their composed raw materials.