Carotenoids Extraction from Orange Peels Using a Thymol-Based Hydrophobic Eutectic Solvent

: The food industry produces substantial amounts of waste, which can cause a lot of environmental issues. However, such waste is also a valuable source of bioactive substances that can potentially be used either by the food industry or other types of industries, in the production of medicines, nutraceuticals, cosmetics, etc. The present study proposes a novel approach to extract such bioactive compounds from orange peel waste using hydrophobic eutectic solvents synthesized with thymol and fatty acids (hexanoic and octanoic acid). A response surface methodology was employed to optimize the extraction conditions and achieve maximum recovery of carotenoids. The optimal hydrophobic eutectic solvent consisted of thymol and hexanoic acid at a molar ratio of 2:1, and the optimum extraction was achieved using a solvent-to-solid ratio of 12:1 and a temperature of 20 ◦ C for 78 min; this resulted in a recovery of 259.45 µ g of total carotenoids per g of dry matter, which is a signiﬁcantly higher recovery compared to common organic solvents. Based on the above, it is demonstrated that hydrophobic eutectic solvents is a promising solvent that can be used to extract bioactive compounds from orange peel waste.


Introduction
The food industry generates a substantial amount of solid and liquid waste as a result of the food production process.With regard to the orange-related industries, the by-products account for about 50% of the fresh fruit weight [1].Out of these waste, 60-65% consists of peels [2], while the rest consists of seeds, pulp, membranes, etc.The disposal of orange peels is a potential environmental pollution issue [3].However, if valorized properly, they can serve as a good source of valuable constituents such as pectin, carotenoids, and flavonoids.These substances can be used as food additives (e.g., colorants, flavorings, and antioxidants), thereby promoting acceptance of food products or improving human health [4], as well as can be utilized in pharmaceuticals, cosmetics, and other industries.To this end, various methods are being applied in order to extract beneficial compounds, such as carotenoids, from orange peels [5].However, most of the extraction processes use organic solvents such as acetone, hexane, ethanol, etc.According to published research, acetone can provide an extraction yield of 7.88 ± 0.59 µg/g dry weight [6], hexane can provide an yield of 113.5 ± 0.4 mg/100 g fresh weight, and ethanol can provide an yield of 110.3 ± 0.3 mg/100 g fresh weight [7].These solvents pose risks to human health due to their flammability and toxicity, and are harmful to the environment [8].To address the issue of organic solvent usage, as well as to promote more "green" practices, in 2010, Anastas and Eghbali [9] proposed 12 principles for greener procedures in chemistry and chemical engineering, highlighting the importance of employing so-called green solvents.
One relatively new and promising category of green solvents is eutectic solvents (ESs), also known as deep eutectic solvents (DESs).A eutectic mixture is a mixture of two or more compounds, with relatively high melting points, that exhibits a significantly lower melting point than the individual components.The decrease in the melting point of DESs mainly stems from the formation of a hydrogen bond network between the hydrogen bond donor (HBD) and the hydrogen bond acceptor (HBA).As such, even though the HBDs and HBAs used are typically solid at room temperature, the resulting mixture becomes liquid and can be utilized as a solvent.Although this property of eutectic mixtures is not new, it was not until 2002 that Abbott et al. [10] showed that eutectic mixtures of choline chloride and urea can be used as green solvents with unique properties.DESs can be easily synthesized via simple methods [11].DESs exhibit several advantages over conventional organic solvents, including low vapor pressure, easy and cost-effective preparation, and non-toxicity [12].More importantly, DESs are considered to be tunable solvents as their physicochemical properties can be tailored by properly selecting the HBA and HBD, by adjusting their molar ratio, or even by the introduction of a third compound in the eutectic mixture, rendering it more suitable for the extraction of specific compounds [13].Finally, since the synthesis of DESs is straightforward, with the synthesis yield being 100% and no requirement of purification steps, DESs are excellent candidates for large-scale industrial applications [14].Although DESs have been used for the extraction of various compounds, attempts were made to synthesize new DESs based on natural components, already present in food matrices, so as to negate the need to separate the target compounds from the DES.As such, the natural deep eutectic solvents (NADESs) were developed.NADESs exhibit excellent extraction properties, and in some cases, they exhibit even more favorable characteristics, such as the deceleration of anthocyanin degradation [15].
However, DES and NADES have some disadvantages that limit their extraction potential.One major disadvantage is their high viscosity, which limits their extraction efficiency and can lead to challenges in industrial-scale applications due to heat and mass transfer limitations during extraction or dissolution [16].Although the viscosity of NADES can be significantly reduced by adding water, the increase in water content in the NADES weakens the hydrogen bonding interactions, and might even completely disrupt them if the water content exceeds 50% [17].Another limitation of (NA)DES is their hydrophilicity resulting in inefficient applicability for non-polar substances.In 2015, Van Osch et al. [18] created a DES based on a quaternary ammonium salt and decanoic acid, and used it for the extraction of volatile fatty acids from aqueous solutions.Since this solvent was hydrophobic, a new category of DES, called hydrophobic DES (HDES) was created, paving the way for a series of applications such as the extraction of pesticides from aqueous environments [19], enzyme activation [20], dye-sensitized solar cells [21], etc.
The aim of this study was to conduct the synthesis of HDES based on thymol and fatty acids (hexanoic and octanoic) for extracting carotenoids from orange peel waste.In recent years, many studies have explored the use of NADES for the extraction of bioactive compounds, such as phenolic compounds [22][23][24][25]; however, fewer studies have focused on the extraction of carotenoids [26,27].Since carotenoids are lipophilic substances, a hydrophobic DES based on thymol could prove to be a good extraction solvent, and to the best of our knowledge, a thymol-based HDES has never been used for the extraction of carotenoids from orange peels.To this end, various syntheses were evaluated, so as to select the optimum HDES.Next, the extraction parameters were optimized, in order to further enhance the extraction process, using response surface methodology (RSM) [28,29].As carotenoids are soluble in acetone, ethyl acetate, and hexane [30], these solvents were used for comparison.

Orange Sample Preparation
Fresh oranges of the Merlin (Washington navel) cultivar (Citrus sinensis L.), collected at the maturity stage, were purchased from selected orchards in the Argos area (Peloponnese, Greece) in December 2022.The physical and chemical parameters of the oranges were evaluated to ensure that they were at the maturity stage.The color of the fruit limb/twig was green, the fruit rind (epicarp) color was orange, and the weight of the fruits ranged from 150 to 200 g.Furthermore, the total soluble solids (TSSs) ( • Brix) were between 10 and 12, the titratable acidity (TA) (expressed as % citric acid) was >0.8, and the TSS/TA ratio was >12.
The oranges were washed with tap water and dried with a paper towel.The peels were removed manually, cut into smaller pieces (~2 × 2 cm), and placed in a Biobase BK-FD10P freeze-dryer (Jinan, China) for 24 h in order to remove water.The freeze-dried peels were then pulverized and placed in sieves, so as to separate them according to size.The powdered orange peels with an average particle diameter of 470 µm were used for the preparation of the extracts (the amount of the powder from the smaller average particle diameter was substantially smaller compared to the used particles, and therefore was not used as we opted for the extraction of compounds from the major fraction of the powder).

HDES Synthesis
A total of 9 HDES were synthesized.Thymol was used as the HBA, whereas hexanoic acid and octanoic acid were used as the HBDs.Since the fatty acids have hydroxyl groups (−OH), they can either behave as HBAs or HBDs, so the combination of the two fatty acids was also tested.In the case of an HDES consisting of two fatty acids, the one with the longest chain behaves as the HBD [31].For each combination of HBA and HBD, molar ratios of 2:1, 1:1, and 1:2 were chosen, as shown in Table 1.To synthesize the HDES, appropriate amounts of HBA and HBD were mixed in 25 mL glass bottles and heated at 70 • C under stirring at 350 rpm until a clear homogeneous liquid was formed.The HDES were then allowed to cool to room temperature, and inspected for the formation of crystals after 24 h.
The experimental methodology for density measurement involved precise sample preparation, temperature calibration, and density determination.All HDES samples were carefully transferred into pre-weighed density bottles, and the masses were recorded.Reproducibility was ensured through multiple measurements.

Extraction Process
Initially, all HDES were assessed to examine their efficiency for the extraction.To this end, the HDES and the orange peels were mixed at a solvent-to-solid ratio of 10:1, and the mixture was stirred at 500 rpm for 60 min at room temperature.For means of comparison, extracts were also prepared using three conventional organic solvents (i.e., hexane, ethyl acetate, and acetone).
Using the optimum HDES, the overall extraction process was further optimized.For the optimization of the extraction parameters, an RSM was employed.The parameters to be tested were solvent-to-solid ratio, temperature, and time.An experiment using a Box-Behnken design with 15 design points, including 3 central points, served as the foundation for the optimization, and the response was the extraction yield in total carotenoids (Y TCn ).The process variables were established in three levels (Table 2).The analysis of variance (ANOVA) and the summary-of-fit tests were used to evaluate the overall model significance (R 2 , p) at a minimum level of 95%.

Total Carotenoid Content (TCC)
A method developed by Biswas et al. [32] was implemented for the determination of the total carotenoid content of the extracts.Briefly, 200 µL of the sample was mixed with 800 µL of ethanol, and the solution was vigorously shaken for 30 s.The absorbance was read at 428 nm using a Shimadzu spectrophotometer (UV-1700, Shimadzu Europa GmbH, Duisburg, Germany).A calibration curve, created with β-carotene (1-10 µg/mL) as a standard compound, was used for the determination of the TCC.

Antiradical Activity (DPPH Assay)
The antiradical activity of the extracts was determined via the 2,2-diphenyl-1-picrylhydrazyl (DPPH) assay using a previously described method [33].An aliquot of 200 µL of the sample was mixed with 800 µL of ethanol, and the solution was vigorously shaken for 30 s.Briefly, 25 µL of the solution was mixed with 975 µL of DPPH solution (100 µM).After mixing, the absorbance of the solution was read at 515 nm (A 515(i) ).The solutions were incubated for 30 min in the absence of light, and the absorbance was read at 515 nm (A 515(f) ).The capacity to scavenge the DPPH radical was expressed as The antiradical activity (A AR ) was expressed as µmol ascorbic acid equivalents (AAE), using an ascorbic acid calibration curve (C AA , 50-1.000 µM) and the following Equation (2): where V is the volume of the extraction medium (in L) and w is the weight of the dry matter (in g).

Ferric Reducing Antioxidant Power (FRAP) Assay
The ferric-reducing power was determined using a previously described method [34].An aliquot of 200 µL of the sample was mixed with 800 µL of ethanol, and the solution was vigorously shaken for 30 s.Briefly, 50 µL of the solution was mixed with 50 µL of FeCl 3 solution (4 mM in 0.05 M HCl), and the solutions were incubated at 37 • C for 30 min.After incubation, 900 µL of 2,4,6-tris(2-pyridyl)-s-triazine (TPTZ) solution (1 mM in 0.05 HCl) was added, and 5 min later, the absorbance was read at 620 nm.A calibration curve, with ascorbic acid (C AA , 50-500 µM) as a standard compound, was used, and the P R was determined as µmol ascorbic acid equivalents (AAE) per g of dm (dry matter), using the following Equation (3): where V is the volume of the extraction medium (in L) and w is the weight of the dry sample (in g).

Color Analysis
The color analysis of the extracts was conducted through measuring the absorbance in three different wavelengths (420, 520, and 620 nm), as well as using a Lovibond CAM-System 500.

Colorimeter Method
A Lovibond CAM-System 500 Imaging Colorimeter was used for the color (CIE L*, a*, b*) analysis of the extracts.The parameter L* indicates lightness, a* redness, and b* yellowness [35].The values of these parameters were used to calculate Chroma (C * Ab ) and hue-angle (h o Ab ), according to Equations ( 4) and ( 5), respectively.

Absorbance Method
An aliquot of 200 µL of the sample was mixed with 800 µL of ethanol, and the solution was vigorously shaken for 30 s.The absorbance of the solution was read at 420, 520, and 620 nm.The sum of the absorbance at the above wavelengths provides the color intensity (CI), and the ratio of the absorbance at 420 nm over the absorbance at 520 nm provides the hue (Equations ( 6) and ( 7)) [36].Also, the color composition was calculated, i.e., the contribution of the three components (yellow, red and blue, expressed as a percentage), using the following Equations ( 8)-( 10): Blue (%) = A 620 CI × 100 (10)

Statistical Analysis
The experimental design, statistical analysis related to the response surface methodology, and distribution analysis were all created using the JMP ® Pro 16 (SAS, Cary, NC, USA) software.The results were expressed as mean values, and the standard deviation was also calculated.

Results and Discussion
In this study, we investigated the use of eutectic solvents derived from thymol and fatty acids (hexanoic and octanoic acid) at various molar ratios (2:1, 1:1, and 1:2).While the term "DES" typically refers to a specific eutectic composition, we employ it here to encompass a range of compositions within the vicinity of the eutectic point.These ratios were examined, since they are common to form eutectic mixtures for both HBD and HBA, in order to evaluate the potential of the eutectic mixtures as solvents for carotenoid extraction.The relatively high density and viscosity values of DES are limiting factors in their applicability as extraction solvents.Since hydrophobic DES generally have lower densities and viscosities compared to hydrophilic ones [37][38][39], we proposed the use of a HDES synthesized from thymol and fatty acids to be used as a solvent for the extraction of carotenoids from orange peels.

Choice of Solvent
A total of nine solvents based on thymol and fatty acids (hexanoic and octanoic) were synthesized at three molar ratios.These HDES, along with three conventional organic solvents were used to extract carotenoids from orange peel waste.The results are shown in Table 3.The absorbance spectra (not depicted) showed a peak at 428 nm, indicating the presence of violaxanthin in the extracts, which is the major carotenoid in orange peels [40].Statistically significant differences (p < 0.05) are denoted with different superscript letters (e.g., a-c ).
Although carotenoids are mainly non-polar substances, the epoxy functional group in violaxanthin increases its polarity [41], so more polar solvents would achieve a higher extraction yield.According to Martins et al. [42], the aromatic ring in the structure of thymol results in high values of the solvatochromic parameter π*, which is related to the polarizability/dipolarity of its mixtures; so in our case, the thymol-based HDES should present high extraction yields for the polar carotenoids in the orange peels.
As can be seen in Table 3, the highest extraction yield was achieved with Thy/Hex 2:1, followed by Thy/Oct 2:1.The majority of the other HDES had comparable extraction yields, except for Thy/Oct 1:1, which proved to be the least efficient.When octanoic acid was used as the HBD, the extraction yield was lower than when hexanoic acid was the HBD.This can be attributed to the fact that when the number of carbon atoms in the fatty acid increases, the polarity of HDES decreases [43].
As far as the organic solvents are concerned, hexane and acetone exhibited the lowest extraction yields 141.58 ± 3.49 and 128.15 ± 6.9 µg CtE/g dm, respectively, while ethyl acetate yielded much better results (171.2 ± 8.16 µg CtE/g dm), similar to most of the HDES (max 184.94 ± 4.69 µg CtE/g dm for Thy/Hex 2:1).This is because ethyl acetate presents higher solubility for polar carotenoids compared to hexane, as predicted via the Conductorlike Screening Model for Real Solvents (COSMO-RS) analysis, using COSMOthermX program (version C30 release 13.01) [44].HDES synthesized exclusively of fatty acids had very good extraction yields, and is reported to be the most efficient one in a number of recent publications, i.e., a mixture of capric acid and lauric acid at a molar ratio of 1:2 for the extraction of carotenoids from tomato achieving an extraction of 7.90 ± 0.15 mg lycopene/100 g fresh weight [45], and a mixture of caprylic acid and capric acid at a molar ratio of 3:1 for the extraction of 151.41 mg β-carotene/mL from pumpkin [27].Fan et al. [46] found that an HDES composed of thymol and fenchyl alcohol proved to be the most efficient for the recovery of lutein from microalgae, achieving a yield of 6.26 ± 0.40 mg/g.Based on the above, the synthesized HDES can be used as an extraction solvent, achieving comparable or enhanced extraction yields as compared to organic solvents.Thus, Thy/Hex 2:1 was selected as the optimal HDES.

Extraction Optimization
The Box-Behnken Design (BBD) was chosen because it is specifically constructed to fit a second-order model, which is the main focus of most RSM investigations.Additionally, a BBD frequently requires fewer experimental runs.Thus, a BBD was applied to test the effect of solvent-to-solid ratio (X 1 ), extraction time (X 2 ), and temperature (X 3 ) on total carotenoid content (TCC).Fifteen experiments were conducted, as seen in Table 4, and the results were analyzed using ANOVA (Figure 1) to evaluate the statistical significance of the model.Model fitting was assessed with the square coefficient of correlation (R 2 ), which was over 0.96, and suggested a satisfactory agreement between the measured and predicted values.Figure 2 shows the generated response surface plots.The optimum extraction conditions, as calculated using the statistical analysis, were a solvent-to-solid ratio of 12:1, an extraction time of 78 min, and a temperature of 20 • C.   The final extractions were performed at the optimum conditions using Thy/Hex 2:1, whereas hexane, ethyl acetate, and acetone were also examined for comparison.The final extractions were performed at the optimum conditions using Thy/Hex 2:1, whereas hexane, ethyl acetate, and acetone were also examined for comparison.The interaction between solventto-solid ratio (X 1 ) and extraction time (X 2 ), (B) the interaction between solvent-to-solid ratio (X 1 ) and temperature (X 3 ), and (C) the interaction between extraction time (X 2 ) and temperature (X 3 ).
The final extractions were performed at the optimum conditions using Thy/Hex 2:1, whereas hexane, ethyl acetate, and acetone were also examined for comparison.

TCC of the Extracts
The determination of the total carotenoids (Table 5) showed that HDES achieved a better extraction yield than all of the conventional organic solvents, with ethyl acetate being the most efficient of the three solvents, followed by hexane and acetone, which was the least efficient.It is significant that the quantity of the TCC extracted via HDES was 59% higher (statistically significant at p < 0.05) compared to hexane, which is the most commonly employed organic solvent for the extraction of carotenoids; it was also 37% and 76% higher compared to ethyl acetate and acetone, respectively.These data are supported by a number of recent publications which report that, as compared to hexane, HDES achieved higher extraction yields of carotenoids from orange peels [26], tomato pomace [47], and kale [48].Additionally, the TCC of 259.45 ± 3.46 µg/g dm that was obtained under the optimum extraction conditions was very close to the predicted value of 255.93 ± 21.61 µg/g dm, which suggested a good prediction ability of the model.Even though the extraction yield of 259.45 ± 3.46 µg/g dm may not seem substantial, considering the market price of USD 300-3000 per kg of β-carotene and the market volume growth [49], along with the cost of other methods of carotenoid production, e.g., the cost of microbial growth medium for the microbial production of pigments [50], the use of HDES for carotenoid extraction could prove to be more efficient and sustainable.Within each column, statistically significant differences (p < 0.05) are denoted with different superscript letters (e.g., a-d ).

Antioxidant Properties of the Extracts
As seen from the results in Table 5, the extract obtained with HDES exhibited better antioxidant properties (statistically significant at p < 0.05), quantified using the FRAP assay and the DPPH assay, compared to the three other solvents.
More specifically, the FRAP assay results showed that the extracts obtained with ethyl acetate and acetone exhibited similar reducing power, whereas the hexane extract had almost double reducing power.The extract obtained with Thy/Hex 2:1 had a much higher reducing power compared to ethyl acetate and acetone (152% and 174%, respectively), and was 38% higher compared to hexane.
The comparison of the DPPH assay results further validated that the extract obtained using hexane had better antioxidant activity than the ones obtained with ethyl acetate and acetone, but again, the Thy/Hex 2:1 extract was found to achieve an increased scavenging activity, from 198% up to 335% (statistically significant for p < 0.05), compared to the organic solvents.
These results are in accordance with recent publications.Alsaud et al. [51] measured the antioxidant activity of Manuka leaves extract, and found that the extract obtained using HDES showed higher antioxidant activity than the one obtained with hexane.In a study of orange peel, Viñas-Ospino et al. [7] showed that extracts obtained using hydrophobic DES had higher antioxidant activity compared to conventional solvents and hydrophilic DES.According to some researchers, in addition to being better at extracting substances with antioxidant activities, (NA)DES can also enhance the antioxidant activity of the extracts and improve their storage stability [52,53].In our case, the higher antioxidant activity of the Thy/Hex 2:1 extract could be attributed to the solvent itself, since the antioxidant properties of thymol have been well studied in the past [54][55][56][57].

Color Analysis of the Extracts
The values of the color parameters obtained with the Lovibond colorimeter are shown in Table 6.The extracts obtained with the three organic solvents had very similar color saturation and hue values (expressed by the chroma and hue-angle, respectively).Their hue-angles were very close to 90 degrees (yellow), whereas the extract obtained with Thy/Hex 2:1 was the farthest (70.8 degrees), which means less yellow and redder (90 degrees hue angle is yellow and zero degrees is red), as can also be seen in Figure 3. Within each column, statistically significant differences (p < 0.05) are denoted with different superscript letters (e.g., a-c ).
The values of the color parameters obtained with the Lovibond colorimeter are shown in Table 6.The extracts obtained with the three organic solvents had very similar color saturation and hue values (expressed by the chroma and hue-angle, respectively).Their hue-angles were very close to 90 degrees (yellow), whereas the extract obtained with Thy/Hex 2:1 was the farthest (70.8 degrees), which means less yellow and redder (90 degrees hue angle is yellow and zero degrees is red), as can also be seen in Figure 3.  Within each column, statistically significant differences (p < 0.05) are denoted with different superscript letters (e.g., a-c ).
The hue-angle (h°) data support our earlier findings about the carotenoid contents in the extracts, as according to Kishimoto et al. [58], the redness of the extract increases as the carotenoid content increases, which is the case here.As an alternative method to analyze the color parameters, the absorbance of the extracts (after proper dilution) was measured at 420, 520, and 620 nm.The calculated values for color intensity (CI), hue (h), and the percentage of each color are given in Table 7.The hue-angle (h • ) data support our earlier findings about the carotenoid contents in the extracts, as according to Kishimoto et al. [58], the redness of the extract increases as the carotenoid content increases, which is the case here.
As an alternative method to analyze the color parameters, the absorbance of the extracts (after proper dilution) was measured at 420, 520, and 620 nm.The calculated values for color intensity (CI), hue (h), and the percentage of each color are given in Table 7.It can be seen that the extracts obtained with the organic solvents had similar color attributes, whereas the one obtained with HDES differed.The Thy/Hex 2:1 extract had a lower percentage of yellow and a higher percentage of red, which again supports our earlier results about the higher concentration of carotenoids.This finding is also corroborated by the correlation analysis, as described in the next section.The lower percentage of yellow may arise from the decreased co-extraction of other pigments such as anthocyanins or chlorophyll derivatives, that may influence the observed colors.

Principal Component Analysis (PCA) and Multivariate Correlation Analysis (MCA)
In order to minimize the dimensionality of the multivariate data and obtain an improved understanding of the results, a principal component analysis (PCA) was also performed.Figure 4 shows that the two main components could explain 92.9% of the variation (eigenvalues > 1), which was considered to be a statistically significant parameter (p < 0.0001).PC1 demonstrated a positive correlation with total carotenoids, DPPH, FRAP, color intensity (CI), and redness (a*), and a negative correlation with other color indices (L*, b*, C*, h • , and hue).PC1 also explained 82.2% of the variability.With the exception of several color characteristics (a*, b*, C*, and hue), PC2 can account for 10.7% of the variance, and shows a positive correlation with all of the variables examined.According to the PCA plots in Figure 4, the loading direction of DPPH, FRAP, and the total carotenoids is the same with the Thy/Hex 2:1 sample parameter; however, the loading direction of C* is different.Because of this, it can be assumed from Figure 5 that the total carotenoids are positively correlated (>0.9) with the antioxidant parameters, but negatively correlated (−0.9) with C*.The stronger the correlation between the measured variables, the closer it is to 1. Additionally, the best correlation was found between CI and total carotenoids (0.99), which was found to be a statistically significant factor (p < 0.0001).

Conclusions
In this study, a series of hydrophobic deep eutectic solvents was synthesized to test their efficiency in the recovery of carotenoids from orange peels.The initial extractions indicated that HDES composed of thymol and hexanoic acid at a molar ratio of 2:1 achieved the highest extraction yield, and it was used for further optimization of the extraction parameters.Under optimum conditions, the extract obtained with Thy/Hex 2:1 had a higher carotenoid concentration compared to hexane, acetone, and ethyl acetate; also, the extracted compounds exhibited enhanced antioxidant properties, as evidenced by the DPPH and FRAP assay results.Based on the above results and the fact that the process of HDES preparation is economical, easy, and 100% atom efficient, we can conclude that HDES prepared from thymol and hexanoic acid provides a viable and sustainable alternative to conventional organic solvents for the extraction of carotenoids from orange peels.

Figure 1 .
Figure 1.Plot of predicted vs. actual values of the response (YTCn) (plot A), and desirability function (plot B), describing the effect of three independent variables considered (solvent-to-solid ratio X1, extraction time X2, temperature X3) on the YTCn upon simultaneous variation.Asterisks and colored values denote statistically significant values, while inset tables include statistics relevant to the evaluation of the resulting model.

Figure 1 .Figure 1 .
Figure 1.Plot of predicted vs. actual values of the response (Y TCn ) (plot A), and desirability function (plot B), describing the effect of three independent variables considered (solvent-to-solid ratio X 1 , extraction time X 2 , temperature X 3 ) on the Y TCn upon simultaneous variation.Asterisks and colored values denote statistically significant values, while inset tables include statistics relevant to the evaluation of the resulting model.

Figure 2 .
Figure 2. Response surface plots for the total carotenoid yields.(A)The interaction between solventto-solid ratio (X 1 ) and extraction time (X 2 ), (B) the interaction between solvent-to-solid ratio (X 1 ) and temperature (X 3 ), and (C) the interaction between extraction time (X 2 ) and temperature (X 3 ).

Figure 4 .
Figure 4. Principal component analysis (PCA) for the measured variables.The axis scores fo and PC2 are displayed.The inset table includes the eigenvalues.Asterisks and colored value note statistically significant values.

Figure 4 .
Figure 4. Principal component analysis (PCA) for the measured variables.The axis scores for PC1 and PC2 are displayed.The inset table includes the eigenvalues.Asterisks and colored values denote statistically significant values.

Figure 4 .
Figure 4. Principal component analysis (PCA) for the measured variables.The axis scores for and PC2 are displayed.The inset table includes the eigenvalues.Asterisks and colored values note statistically significant values.

Figure 5 .
Figure 5. Multivariate correlation analysis of measured variables.Figure 5. Multivariate correlation analysis of measured variables.

Figure 5 .
Figure 5. Multivariate correlation analysis of measured variables.Figure 5. Multivariate correlation analysis of measured variables.

Table 2 .
Levels of independent variables in actual and coded values for the optimization of the extraction of total carotenoids.

Table 3 .
Extraction yield of total carotenoids (Y TCn ) obtained via the prepared HDES, and conventional solvents expressed as µg of total carotenoids per g of dry matter.

Table 4 .
Experimental design, measured and predicted responses of the dependent variable expressed as µg of total carotenoids per g of dry matter.

Table 5 .
Extraction yield of total carotenoids and antioxidant properties of the extracts obtained with different solvents at optimum extraction conditions.

Table 6 .
Color parameters of the extracts.Measurements were obtained with a colorimeter.

Table 6 .
Color parameters of the extracts.Measurements were obtained with a colorimeter.

Table 7 .
Color parameters of the extracts.Measurements obtained with UV-Vis spectrophotometer.

Table 7 .
Color parameters of the extracts.Measurements obtained with UV-Vis spectrophotometer.