Exploration of High-Nutritional-Quality Vegetable Oil Blend with Enhanced Oxidative Stability as a Frying Medium Substitute for Palm Oil

: Blending is a commonly utilized technique for enhancing the oxidative stability, nutritional quality, and physicochemical properties of vegetable oils. This study explored the potential of a vegetable oil blend consisting of common seed oils (sunflower, soybean, rapeseed, cottonseed


Introduction
Oils and fats are an important source of energy to humans, since they supply calories and provide the body with fat-soluble vitamins.Additionally, they serve as precursors to various hormones and supply essential fatty acids [1].The fatty acid compositions of some commonly used vegetable oils vary significantly [2].For instance, coconut and palm oils have a high content of saturated fats (i.e., ~90 and ~44%, respectively), whereas seed oils (i.e., corn, sunflower, or soybean oil) have mostly polyunsaturated fatty acids.There is a strong relation between the economic, organoleptic, and functional qualities of foods and food products and the types of vegetable oils or fats used [3].Olive oil and other vegetable oils, such as mid-oleic sunflower oil and rapeseed oil containing cis-monounsaturated and polyunsaturated fatty acids, are known to provide health benefits to consumers, such as protection from cardiovascular diseases [4].In fact, the monounsaturated-topolyunsaturated fatty acid ratios of conventional vegetable oils such as soybean (1.5:3.8),rapeseed (9.7:3.9),sunflower (2.5:5.9),cottonseed (2.0:5.7), and corn (2.0:3.0)oils were previously documented [5].In terms of thermal stability, it is widely known that the higher the amount of double bonds in a fatty acid, the greater the susceptibility to oxidation [6].On the contrary, trans fatty acids shows exceptional thermal stability, but a high consumption rate could dysregulate biological pathways that affect cardiometabolic health [7].
The oil obtained from the palm tree (Elaeis guineensis) has gained popularity as a major edible oil in recent decades in the food industry.Palm oil has been used extensively in Western Europe as a substitute for partially hydrogenated oils in food applications given that it is solid at ambient temperature and its fractions provide a broad spectrum of melting profiles [8].The high content of saturated fatty acids is a concern to human health when compared to other conventional vegetable oils.It should be noted that the yield of palm oil per year is the highest of any edible oil.However, palm oil is also a topic of debate, since it causes negative environmental issues such as the destruction of vast rainforests.A partial or total reduction of this oil may provide a long-term solution to deforestation while simultaneously having positive economic and environmental effects [9].To that end, legislation in an increasing number of nations is imposing restrictions, which has led to a growing interest in finding alternatives to palm oil [10].
The blending of vegetable oils possessing distinct properties represents a straightforward approach to developing novel products that possess the intended sensory and oxidative characteristics [11].A single vegetable oil may not meet consumer standards for oxidative stability, physical and chemical properties, or nutritional quality.The blending of oils alters their triacylglycerol profiles, which subsequently impacts various physical properties of the oils, including their density, viscosity, cloud point, solid fat content, and sensory quality [12].Furthermore, oils having different amounts of unsaponifiable material, antioxidant compounds, and fatty acid compositions would fry differently and reduce hazardous compounds [13].Bioactive compounds such as tocopherols, polyphenols, chlorophylls, and carotenoids could also contribute to the oxidative stability of oils, since they scavenge free radicals [6].It was previously mentioned that the formation of hazardous polar compounds, such as triglyceride oligomers or diglycerides, during deep frying was more significantly influenced by the type of oil rather than by the type of food.Trans fatty acid synthesis is directly related to the oxidative stability of oils; thus, elements that speed up lipid oxidation will also raise the trans fatty acid content [14].These reactions may significantly impair the quality of oils as the frying process progresses [15].Therefore, by identifying the properties of oils and utilizing suitable oils for blending, it is possible to enhance the nutritional and sensory attributes of the product [16].Given that the hydrogenation of fatty acids produces dangerous trans fatty acids and lowers the quality of the oil, specific oil blending could improve the oil stability under frying and replace these types of oils in the food industry [17].
Palm oil has been used for the frying process due to its high stability at high temperatures [18].Several seed oils, produced from plants that can be cultivated in various climates worldwide, with high concentrations of polyunsaturated fatty acids, such as corn oil or sunflower oil, are also utilized for frying in the food sector [19,20].To the best of the authors' knowledge, studies examining the blending of vegetable oils using multiple different seed oils are limited.To that end, this study aimed to create a vegetable oil blend with five widely commercially available seed oils used in the food sector (i.e., sunflower, soybean, corn, rapeseed, and cottonseed oils) as a palm oil substitute through partial least squares.The nutritional quality and oxidative stability of the blend were compared to palm oil.The results of this study could indicate a suitable oil blend of high nutritional quality that could be widely used in the food industry.

Vegetable Oil Sample Procurement
Fresh, refined vegetable oils, including soybean oil, sunflower oil, rapeseed oil, cottonseed oil, corn oil, and palm oil, were all purchased from a local market (Karditsa, Greece).All samples were stored in the refrigerator at a temperature < 4 • C. To achieve thermal equilibrium between them and room temperature, they were kept in the dark (20 • C) prior to analysis and were immediately analyzed after opening.Each prepared blend was homogenized by stirring for at least 20 min with a Heidolph magnetic stirrer (Heidolph Instruments GmbH & Co. KG, Schwabach, Germany).

Blended Vegetable Oil Mixture Preparations and Thermal Treatment
A diverse range of blends was produced by combining pure edible vegetable oil samples at different concentrations.The variables X 1 -X 5 represent each vegetable oil, excluding palm oil, and their respective ratio is indicated in the Coded Variable Levels.Further details are available in Table 1 below.A quantity of 10 g of each oil blend and palm oil was placed in Duran bottles and was incubated at 60 • C for 28 days in a water bath in the dark, as per our previous study [21] and other similar approaches [22][23][24].To perform the oil stability analyses, sampling was conducted before incubation and 14 and 28 days after incubation.To shed more light on the quality of the blended oils, pure substrate oils were also included in the study.

Main Effects Screening Design
Factorial experiments benefit substantially from orthogonal arrays with mixed levels when studying factors with different numbers of levels [25].However, many mixed-level orthogonal arrays require a large number of experiments, which presents a major concern.By slightly easing the requirement of orthogonality, nearly orthogonal arrays reduce the number of experiments needed, making them a good fit for such situations [26].To that end, the Main Effects Screening design within a Response Surface Methodology (RSM) was utilized to effectively explore the design space with fewer experiments in order to investigate the optimal ratios of the oil samples and to assess their antioxidant activities and thermal stabilities.To achieve this, the ratio of the oils was adjusted.The optimization method was accomplished using a Screening design of 25 design points with different compositions.The design also included five pure vegetable oils along with pure palm oil.Summary-of-fit tests, the analysis of variance (ANOVA), the overall model significance (R 2 , p-value), and the significance of the model (equations) coefficients were evaluated at a minimum level of 95%, with Equation (1) describing the polynomial plot: where the independent variables are represented as X i and X j , and the predicted response variable is defined as Y k .The terms β 0 , β i , β ii , and β ij describe the model's linear, quadratic, and interaction terms, represented by the intercept and regression coefficients, respectively.
2.5.Oil Quality Tests 2.5.1.Accelerated Oxidation through Rancimat Apparatus The oxidation stability assessment was conducted using a Rancimat 743 apparatus from Metrohm LTD.(Herisau, Switzerland) and was similar to a previous approach [27].The device was connected to 873 Biodiesel Rancimat software (version 1.00).The mass of the oil sample used was 3 ± 0.01 g, with the device being set to a temperature of 110 • C and an air flow rate of 15 L/h.The assessment was conducted on days 0, 14, and 28 of incubation.The beginning of oxidative degradation, or rancidity, was identified as the point of abrupt inflation in conductivity from the plotted graph.

Determination of Antiradical Activity (A AR )
The antiradical activity was evaluated using an established methodology [28].Briefly, 1.0 g of oil was dissolved in 10 mL of ethyl acetate.Then, 50 µL of each oil sample was mixed with 1950 µL DPPH • solution (100 µM in ethyl acetate) and vigorously shaken for 10 s.The absorbance was immediately recorded at 515 nm after mixing (A 515(i) ) and after exactly 30 min (A 515(f) ).The % inhibition of DPPH • was calculated according to the following Equation (2): The antiradical activity was calculated as µmol of Trolox equivalent antioxidant capacity (TEAC) per kilogram of oil (µmol TEAC/kg oil) through a calibration curve (50-500 µM Trolox).
The percentage of oxidative stability of oils (%OS) was calculated in both Rancimat induction period (IP) and DPPH values on days 0, 14, and 28 of incubation using the following equation: Oxidative Stability (%)= Conjugated dienes (CDs) and conjugated trienes (CTs) were measured following our previously described method [29].Cyclohexane was added to a 5 mL volumetric flask along with 0.01 g of oil sample.The concentrations of CDs were detected at 232 nm, while the concentrations of CTs were detected at 270 nm.
The absorbances at 232 nm and 270 nm are denoted as A 232 and A 270 , respectively; ε is the molar absorptivity of linoleic acid hydroperoxide, which is 2.525 × 10 4 M −1 cm −1 , and l is the path length of the cuvette (1 cm).To express the concentrations of CDs (C CD ) and CTs (C CT ) in mmol per kg of oil, the factor 5 × 10 3 accounts for the volume of the solvent (5 mL) used to dissolve the oil sample, and w is the weight of the oil sample in grams.
The percentage of oxidative power (%OP) of oils was calculated in both CD and CT values on days 0, 14, and 28 of incubation using the following equation: Oxidative Power (%) = CD,CT value 28 d − CD, CT value 0 d CD,CT value 0 d × 100 (6)

ATR-FTIR Spectroscopy Acquisition
The FTIR spectra were obtained using a Shimadzu Fourier-Transform Infrared Pres-tige21 spectrophotometer (Shimadzu Corporation, Kyoto, Japan) equipped with an Attenuated Total Reflectance (ATR) accessory with a trough plate composed of ZnSe crystal.The instrument was equipped with a deuterated triglycine sulphate doped with an L-alanine (DLATGS) detector, a KBr beam splitter, high-energy throughput optical elements, and a bright ceramic light source.A total of 32 scans were taken from 0.8 mL of each oil sample, ranging from 4000 to 400 cm −1 at a resolution of 4 cm −1 .The ATR crystal was thoroughly cleaned before and after the measurement with a tissue soaked in acetone to ensure that no oil residue was present.For each sample, all measurements were conducted in triplicate.The FTIR spectra of all samples were processed through Shimadzu IRsolution (version 1.60) software.The spectra of the samples were recorded for an analysis of the baseline and signal-to-noise ratios.

Statistical Analysis
The samples were measured in triplicate for each assay, and the results are shown as means of triplicate determinations.A one-way analysis of variance (ANOVA) test was used to determine the statistical significance of the differences between mean values; a multiple comparisons test (Tukey HSD) was used to determine the statistical significance at p < 0.05.Using JMP ® Pro 16 software (SAS, Cary, NC, USA), the experimental design for the Response Surface Methodology and all related statistics was completed.Pareto plot analysis, partial least squares (PLS) analysis, multiple factor analysis (MFA), and multivariate correlation analysis (MCA) were all performed using the same software.

Results and Discussion
Oil blending is a traditionally employed technique to modify fatty acid profiles while improving the nutritional quality in oils.Blended oils with better properties and increased oxidative stability are a common product of the food industry.De Marco et al. [30] studied the possibility of preparing a binary mixture of sunflower and palm oils in various ratios, creating a mixture richer in tocopherols and UFAs, which would retain its liquid form under ambient conditions and be oxidative and more stable in deep frying conditions.However, they used only one oil to partially substitute palm oil.A similar approach was conducted by Siddique et al. [31], who blended rapeseed, soybean, and sunflower oils with palm oil in order to achieve a stable blend with a low melting point.Several studies also include the use of unconventional oils (i.e., black cumin oil, coriander oil, or pomegranate seed oil) to improve the nutritional quality of vegetable oils [19,32,33].In our study, we opted to use several cheap and conventional vegetable oils with a wide range of MUFAs (~19-62%) and PUFAs (~30-56%) in order to enhance the nutritional quality of the blend and SFAs (~8-28%) to maintain thermal stability.We aimed to examine the potential oxidation-resistant capacity of the optimum blended vegetable oil sample in comparison to palm oil, which has considerably higher long-term oxidative stability [34].The specific oils were chosen as alternatives to palm oil, since they are the most prevalent oils produced in the European Union [22].

Alterations in Antiradical Activity and Oxidative Stability
The Trolox-equivalent antiradical activity (TEAC) and oxidative stability were assessed before and after incubation to determine whether alterations in the oxidative condition of each sample affected their radical-inhibiting capacity.To that end, sampling was conducted initially and after 14 and 28 days following thermal incubation at 60 • C. The results of the Rancimat procedure and DPPH • inhibition activity are shown in Table 2.It was observed that the induction period ranged from 1.12-10.8h in unexposed oil samples.Design point 7, which consisted of sunflower, cottonseed, rapeseed, and corn oils in high proportions, had the lowest value at the induction period but showed moderately high (67.9%)oxidative stability, whereas palm oil had the highest value, indicating high oxidation resistance.It was revealed that a high content of both sunflower and rapeseed oils was not preferable in the blended oil (vide infra).However, design point 17, mainly consisting of rapeseed oil, showed the highest oxidative stability of 87.1%.Despite the low pre-incubation resistance of 2.69 h at 60 • C, it appeared to resist oxidation, demanding 2.34 h for full oxidation after 28 days.Regarding DPPH • inhibition activity, it ranged from 66.51 to 179.07 µmol TEAC/kg of oil.The lowest value was measured in palm oil, whereas design point 13 had the highest value.By the end of the incubation period, the radical inhibition ability decreased in all samples, as expected.However, the decrease ratio was different for each sample, indicating that bioactive compounds naturally present in vegetable oils could explain their inherent radical-scavenging potential [35].

Generation of Oxidation By-Products Assessment
To simulate long-term stability, the conjugated dienes and trienes were measured until day 28 of incubation.Monitoring the CD values helps assess the extent of lipid oxidation and rancidity.CTs are less stable than CDs and could decompose after a certain period when compared to CDs.The results from the CD and CT values along with the % oxidative power (%OP) are shown in detail in Table 3.It can be observed that CDs ranged from 6.77 to 22.38 mmol/kg oil, whereas CTs ranged from 1.28 to 8.48 mmol/kg oil.It should be highlighted that lower CD values indicate better oxidative stability.In addition, a low value of %OP possibly suggested preferable properties, but negative values were not included in the RSM, as a decrease from initial CT values signifies the presence of oxidation products and the degradation of unsaturated fatty acids [36].Interestingly, it was observed that sunflower oil (i.e., design point 27) had the least resistance to oxidation, while design point 3 showed the highest protective effect from oxidation.This was probably caused by the chemical composition of the blended oil when compared to sole sunflower oil.Blending vegetable oils resulted in a protective effect, achieving lower CD and CT values also from palm oil.Similar results were obtained by Ramroudi et al. [12], where sunflower, sesame, and corn oils and their blends were exposed to deep frying for 3 days.It was observed that corn oil was the most stable oil, with the blends also achieving promising results.

Impact of Variables on Oxidative Stability Assays through Pareto Plot Analysis
Pareto plots are a valuable tool in statistical analysis, particularly in the field of quality control and improvement [37].They help in identifying the most significant factors in a dataset, which in this context, are the transformed estimates for the %OS and %OP in various assays.In Figure 1, the reference gold line (or gold reference rectangle) mentioned is crucial, as it denotes the level of significance, commonly set at p < 0.05, which helps in distinguishing between statistically significant effects and those that are not.This differentiation is essential for making informed decisions based on the data.For instance, in the Rancimat and DPPH assays, which measure the %OS, and the CD and CT assays, which assess the %OP, the Pareto plot can visually communicate which factors have the most substantial impact on the outcomes and should therefore be prioritized for optimization or further investigation.The use of such plots simplifies complex data and aids in focusing efforts where they can be most effective, ensuring resources are not wasted on insignificant variables [38].Regarding the variables examined through the Pareto plot, it was observed that the X 2 and X 3 variables had a positive impact on oxidative stability under accelerated conditions using the Rancimat apparatus.As for DPPH • inhibition activity, it was concluded that corn oil highly contributed to the antioxidant activity of blended vegetable oils.This trend was previously noticed, as pure corn oil had the highest antioxidant activity of all pure substrate oils.In addition, the specific oil kept the highest inhibition activity after 28 days of incubation.On the other hand, it was observed that the variables X 2 and X 4 (i.e., sunflower and cottonseed oils) had a negative impact on the CD values, indicating that their blend produced high amounts of primary oxidation products.This may also explain the negative effect on the formation of conjugated trienes.

Multiple Factor Analysis (MFA)
A method that expands the Principal Components Analysis to datasets with several variables assessed on the same items is Multiple Factor Analysis (MFA).When the variables are transformed into orthogonal factors, it becomes possible to compare the variables.These elements show how things are similar and different according to the variable's evaluations.To investigate the connections between the measured variables, we employed

Multiple Factor Analysis (MFA)
A method that expands the Principal Components Analysis to datasets with several variables assessed on the same items is Multiple Factor Analysis (MFA).When the variables are transformed into orthogonal factors, it becomes possible to compare the variables.These elements show how things are similar and different according to the variable's evaluations.To investigate the connections between the measured variables, we employed MFA, the outcomes of which are shown in Figure 2. The factor scores for each measurement variable on the first two dimensions are also shown.These dimensions account for 30 and 20.6% of the total variance, respectively.The analyzed variables (X 1 -X 5 ) are also shown in the figure, with their levels indicated.As a result of their connection in the factor space, the plot shows blocks of elements that are familiar (or similar).It was observed that the position of pure oil samples varied significantly.Samples X 1 and X 4 were positioned nearby, indicating a similar trend between the two oils in the investigated assays.For instance, a negative correlation of both oil samples was observed with %OS from the accelerated oxidation process using Rancimat at 0 and 14 d.On the other hand, a positive correlation between sample X 5 (i.e., corn oil) and DPPH • was observed, an outcome that could be attributed to the chemical composition of the specific oil sample.A measure of the contribution of each set of variables to each dimension, block parti inertia, is shown in Figure 3.The block partial inertia is multiplied by the eigenvalue an divided by 100 to calculate it.This facilitates the evaluation of the extent to which ea block contributes to the overall structure of the data.Dimension 1 accounts for 30% of t total variance within the dataset, representing the most significant pattern or trend.Co versely, Dimension 2 explains 20.6% of the variance, capturing the second most significa pattern, which is orthogonal to the first dimension.The inset tables display the variab loadings, which represent the correlations between the variables and the components, A measure of the contribution of each set of variables to each dimension, block partial inertia, is shown in Figure 3.The block partial inertia is multiplied by the eigenvalue and divided by 100 to calculate it.This facilitates the evaluation of the extent to which each block contributes to the overall structure of the data.Dimension 1 accounts for 30% of the total variance within the dataset, representing the most significant pattern or trend.Conversely, Dimension 2 explains 20.6% of the variance, capturing the second most significant pattern, which is orthogonal to the first dimension.The inset tables display the variable loadings, which represent the correlations between the variables and the components, as well as the contribution of each original parameter to the principal component.The dots of varying sizes represent the significance of the data points.Larger dots signify higher values or greater importance, while smaller dots denote lower values or lesser importance.The relative relevance of the variables in each set is shown by the colors in both Figures 2 and 3, which correspond to the block parameters.

Multivariate Correlation Analysis (MCA)
To further understand the relationship between the variables, an MCA was also ducted.The ability to measure the degree of positive or negative correlation betwee variables is the primary benefit of this analysis compared to the previous.The follo caption indicates that the color map used in this context uses a color scale to repr correlation values ranging from -1 to 1. Figure 4 depicts the findings from this ana An expected trend included the positive correlation with CD and CT values.A po correlation (>0.6) in CD and CT values at 0 days was observed, whereas a strong po correlation (>0.8) was confirmed at 14 days.Interestingly, no strong negative correl between each variable was observed.Regarding thermal-oxidative stability and bioa compound correlation, it should be noted that Rancimat and DPPH variables did not a strong positive correlation.This probably means that the observed antiradical ac was due to bioactive compounds rather than the fatty acid composition.However oxidation resistance observed in Rancimat may not be attributed to the presence of oxidant compounds but rather to the structural characteristics of fatty acids.This i cause the chemical bonds within trans fatty acids show greater oxidative stability pared to the corresponding cis bonds [39].

Multivariate Correlation Analysis (MCA)
To further understand the relationship between the variables, an MCA was also conducted.The ability to measure the degree of positive or negative correlation between the variables is the primary benefit of this analysis compared to the previous.The following caption indicates that the color map used in this context uses a color scale to represent correlation values ranging from −1 to 1. Figure 4 depicts the findings from this analysis.An expected trend included the positive correlation with CD and CT values.A positive correlation (>0.6) in CD and CT values at 0 days was observed, whereas a strong positive correlation (>0.8) was confirmed at 14 days.Interestingly, no strong negative correlation between each variable was observed.Regarding thermal-oxidative stability and bioactive compound correlation, it should be noted that Rancimat and DPPH variables did not have a strong positive correlation.This probably means that the observed antiradical activity was due to bioactive compounds rather than the fatty acid composition.However, the oxidation resistance observed in Rancimat may not be attributed to the presence of antioxidant compounds but rather to the structural characteristics of fatty acids.This is because the chemical bonds within trans fatty acids show greater oxidative stability compared to the corresponding cis bonds [39].

Optimal Blended Vegetable Oil Composition
An inexpensive technique is the use of conventional oils in mixtures [40].Blending vegetable oils has become a widely accepted practice world widely.For instance, combinations of sunflower oil and canola oil or palm oil [30,41], hydrogenated soybean oil and corn oil, or high-oleic sunflower oil and soybean oil are some examples [42,43].To meet industrial and consumer demands, it has also been acceptable to mix conventional edible oil (i.e., sunflower or corn oil) with uncommon oils like rice bran oil.[40].The nutritional quality and stability of oils can be enhanced by blending different vegetable oils, which changes the fatty acid composition and increases the quantities of natural bioactive compounds in the blends [44].
The impact of the examined oil combinations (variables X1, X2, X3, X4, and X5) was explored using a PLS analysis.Visualizing the impact of each vegetable oil on both oxidative stability and inhibitory activity testing, Figure 5A shows the application of the PLS technique to generate a correlation loading plot.The impact of each oil on thermal stability and antiradical activity was revealed in Figure 5B.An increased variable significance for the projection (VIP) factor, particularly beyond 0.8, suggests that this variable makes a larger contribution.According to the results, variables X1, X2, and X3 were observed to have a negative impact on both the inhibition activity and oxidative stability in ascending concentrations.To this end, it was considered preferable for the three specific oils to be present in a low ratio.On the other hand, it was observed that variables X4 and X5 (i.e., cottonseed and corn oils) highly affected the oxidative stability through accelerated oxi-

Optimal Blended Vegetable Oil Composition
An inexpensive technique is the use of conventional oils in mixtures [40].Blending vegetable oils has become a widely accepted practice world widely.For instance, combinations of sunflower oil and canola oil or palm oil [30,41], hydrogenated soybean oil and corn oil, or high-oleic sunflower oil and soybean oil are some examples [42,43].To meet industrial and consumer demands, it has also been acceptable to mix conventional edible oil (i.e., sunflower or corn oil) with uncommon oils like rice bran oil.[40].The nutritional quality and stability of oils can be enhanced by blending different vegetable oils, which changes the fatty acid composition and increases the quantities of natural bioactive compounds in the blends [44].
The impact of the examined oil combinations (variables X 1 , X 2 , X 3 , X 4 , and X 5 ) was explored using a PLS analysis.Visualizing the impact of each vegetable oil on both oxidative stability and inhibitory activity testing, Figure 5A shows the application of the PLS technique to generate a correlation loading plot.The impact of each oil on thermal stability and antiradical activity was revealed in Figure 5B.An increased variable significance for the projection (VIP) factor, particularly beyond 0.8, suggests that this variable makes a larger contribution.According to the results, variables X 1 , X 2 , and X 3 were observed to have a negative impact on both the inhibition activity and oxidative stability in ascending concentrations.To this end, it was considered preferable for the three specific oils to be present in a low ratio.On the other hand, it was observed that variables X 4 and X 5 (i.e., cottonseed and corn oils) highly affected the oxidative stability through accelerated oxidation and DPPH • inhibition activity.As such, an increased content of variables X 4 and X 5 was found preferable to enhance the mentioned activities.It was previously shown that the X 5 variable had the highest %OS in the DPPH • assay.A possible explanation for this trend could lie mostly in the fatty acid composition, as previously revealed using MCA.Regarding the exclusion of rapeseed oil, despite the higher monounsaturated fatty acid content than polyunsaturated fatty acids [5], this ratio may not fully compensate for the high content of carotenoids, which are also vulnerable to oxidation [45].As such, corn and cottonseed oils were shown to have the ideal composition of fatty acids and bioactive compounds to create a nutritious and thermal stable oil blend.
Lipidology 2024, 1, FOR PEER REVIEW 13 the X5 variable had the highest %OS in the DPPH • assay.A possible explanation for this trend could lie mostly in the fatty acid composition, as previously revealed using MCA.
Regarding the exclusion of rapeseed oil, despite the higher monounsaturated fatty acid content than polyunsaturated fatty acids [5], this ratio may not fully compensate for the high content of carotenoids, which are also vulnerable to oxidation [45].As such, corn and cottonseed oils were shown to have the ideal composition of fatty acids and bioactive compounds to create a nutritious and thermal stable oil blend.Consequently, the optimum sample underwent the same assays as palm oil and was compared, the ratio and results of which are shown in Table 4. Figure 6 illustrates the relationship between the induction period and temperature through a polynomial regression fit.This graphical representation is crucial for understanding the kinetics of the reaction under study.The curve obtained from the polynomial regression provides insights into the rate at which the induction period changes with the temperature, offering valuable data for predictive modeling and analysis.Low-temperature (90-110 °C) incubation declares palm oil as more resistant to oxidation.Since saturated fatty acids are more resistant to oxidation than polyunsaturated fatty acids [46,47], this outcome is anticipated.On the contrary, they show a similar pattern at temperatures > 120 °C, providing further evidence that their stability is comparable in this kind of exposure.Furthermore, the inhibition activity was significantly higher (p < 0.05) for the optimal compared to the palm oil by ~40%, rendering it healthier when consumed raw.This outcome could be attributed to the prevalent polyunsaturated fatty acids found in the optimal oil sample.However, the presence of these fatty acids led to higher CD values by 17% before the accelerated oxidation process, since they are more easily oxidized than saturated fatty acids.The oil blending effect led to an induction period until full rancidity at 5.38 h, which was not enough to match the value of palm oil (10.31 h).Consequently, the optimum sample underwent the same assays as palm oil and was compared, the ratio and results of which are shown in Table 4. Figure 6 illustrates the relationship between the induction period and temperature through a polynomial regression fit.This graphical representation is crucial for understanding the kinetics of the reaction under study.The curve obtained from the polynomial regression provides insights into the rate at which the induction period changes with the temperature, offering valuable data for predictive modeling and analysis.Low-temperature (90-110 • C) incubation declares palm oil as more resistant to oxidation.Since saturated fatty acids are more resistant to oxidation than polyunsaturated fatty acids [46,47], this outcome is anticipated.On the contrary, they show a similar pattern at temperatures > 120 • C, providing further evidence that their stability is comparable in this kind of exposure.Furthermore, the inhibition activity was significantly higher (p < 0.05) for the optimal compared to the palm oil by ~40%, rendering it healthier when consumed raw.This outcome could be attributed to the prevalent polyunsaturated fatty acids found in the optimal oil sample.However, the presence of these fatty acids led to higher CD values by 17% before the accelerated oxidation process, since they are more easily oxidized than saturated fatty acids.The oil blending effect led to an induction period until full rancidity at 5.38 h, which was not enough to match the value of palm oil (10.31 h).Table 4. Comparison of palm oil and the optimal sample (X 1 : 2, X 2 : 1, X 3 : 0, X 4 : 4, and X 5 : 4) in the various analyses.

Parameters Palm Oil Optimal
Rancimat (h) Within each row, statistically significant differences (p < 0.05) are denoted with lowercase letters (e.g., a,b).Table 5 presents the statistical analysis results conducted to examine the relationship between the induction period and temperature.A polynomial regression fit was employed, providing a detailed view of how temperature variations can influence the induction phase.The data indicate a significant correlation between these variables, suggesting that the temperature plays a crucial role in the stability of these oils during the induction period, especially that above 120 °C. Figure 6 indicates a fit for both samples, with the optimal one also having a lower value in the root mean square error (RMSE).It was also observed that the palm oil curve had a steeper slope than the blended oil.

ATR-FTIR Spectra Analysis
ATR-FTIR spectroscopy is well known as a quick, easy, and inexpensive way to measure qualitative and quantitative data from fatty acids [48].A further comparison in the chemical compositions of the two oil samples was conducted through ATR-FTIR Table 5 presents the statistical analysis results conducted to examine the relationship between the induction period and temperature.A polynomial regression fit was employed, providing a detailed view of how temperature variations can influence the induction phase.The data indicate a significant correlation between these variables, suggesting that the temperature plays a crucial role in the stability of these oils during the induction period, especially that above 120 • C. Figure 6 indicates a fit for both samples, with the optimal one also having a lower value in the root mean square error (RMSE).It was also observed that the palm oil curve had a steeper slope than the blended oil.

ATR-FTIR Spectra Analysis
ATR-FTIR spectroscopy is well known as a quick, easy, and inexpensive way to measure qualitative and quantitative data from fatty acids [48].A further comparison in the chemical compositions of the two oil samples was conducted through ATR-FTIR analysis, the stacked spectra of which are illustrated in Figure 7.The major wavenumbers in the ATR-FTIR spectra of the two oils typically correspond to the functional groups that are most abundant or reactive in the samples (Table 6), providing valuable insights on sample composition and aiding in more accurate analysis and quality control in the food industry.
analysis, the stacked spectra of which are illustrated in Figure 7.The major wavenumbers in the ATR-FTIR spectra of the two oils typically correspond to the functional groups that are most abundant or reactive in the samples (Table 6), providing valuable insights on sample composition and aiding in more accurate analysis and quality control in the food industry.The spectra show the characteristic bands of the fatty acids in the two oil samples.Firstly, it can be seen that from ~2800-2900 cm -1 are the characteristic symmetric and nonsymmetric C-H bond vibrations of the aliphatic part of the fatty acids.It was clearly shown that the two bands at 2922 and 2853 cm -1 were more intense in palm oil, as they refer to the bonds between C and H in saturated fatty acids.This finding was expected, as palm oil has a greater amount of saturated fatty acids that significantly contribute to its thermal stability.It should be noted that the difference in the specific bands between the two oils was not significant, meaning that the optimum blend also possesses important oxidative stability.On the other hand, however, at 3006 cm -1 the stretching symmetric  The spectra show the characteristic bands of the fatty acids in the two oil samples.Firstly, it can be seen that from ~2800-2900 cm −1 are the characteristic symmetric and non-symmetric C-H bond vibrations of the aliphatic part of the fatty acids.It was clearly shown that the two bands at 2922 and 2853 cm −1 were more intense in palm oil, as they refer to the bonds between C and H in saturated fatty acids.This finding was expected, as palm oil has a greater amount of saturated fatty acids that significantly contribute to its thermal stability.It should be noted that the difference in the specific bands between the two oils was not significant, meaning that the optimum blend also possesses important oxidative stability.On the other hand, however, at 3006 cm −1 the stretching symmetric vibration of the cis double bonds is distinguishable, whereas in the optimal sample, the absorbance is higher, demonstrating the existence of more cis double bonds due to the presence of more unsaturated fatty acids.This result is supported by the absorptions at 912 and 721 cm −1 , where rocking and bending vibrations of -HC=CH-are highlighted, with the absorptions being higher in the optimal sample by up to 40% compared to palm oil.The bands referring to the double bonds between C (i.e., 912, and 3006 cm −1 ) revealed that the optimum sample had the highest amount of unsaturated fatty acids, which have health-promoting benefits.Similar results were found in a study by Poiana et al. [49], who investigated olive oil adulterations with other edible oils.Finally, a band in 1745 cm −1 was observed at similar absorbances in both oil samples, which indicates the presence of an ester carbonyl functional group.

Limitations and Future Perspectives
The objective of this study was to conduct an initial investigation to identify the most suitable combination of seed oils with a balanced ratio of fatty acids and bioactive compounds, thereby classifying it as a nutritious oil.Furthermore, this blended oil would demonstrate comparable oxidative stability to palm oil, potentially leading to a gradual decrease in the reliance on palm oil and its detrimental impact on the natural environment.Nevertheless, it is important to emphasize that future research necessitates the evaluation of additional bioactive compounds, such as tocopherols, carotenoids, and sterols, to elucidate the impacts of these compounds.Finally, for a better understanding of the oxidative stability of oils at high temperatures, a future study could also compare the optimal blend with palm oil at deep-frying temperatures.

Conclusions
To summarize, the global reliance on deep-fried foods necessitates the development of oils of higher nutritional quality.Palm oil is inadequate to meet this demand, as it has plenty of saturated fatty acids.In fact, it may contribute to severe health complications by producing hazardous oxidation by-products.On the other hand, the wide range of monounsaturated and polyunsaturated fatty acids could provide a healthier alternative to several sectors of the food industry.Furthermore, considering that oilseeds could be cultivated in different climates worldwide and produced on a large scale, it could lead to the substitution of palm oil, providing a sustainable way to reduce the detrimental environmental effects of deforestation.A novel blended oil consisting of common vegetable oils such as soybean, sunflower, cottonseed, and corn oils in a 2:1:4:4 ratio, respectively, was prepared and compared to palm oil through several assays (i.e., accelerated oxidation conditions through Rancimat, antiradical inhibition activity, and the measurement of oxidation by-products).The high thermal stability of palm oil due to saturated fatty acids was unmatched by the blended oil; however, its high nutritional quality was revealed through its high antioxidant capacity (~40% higher).This finding was also supported by ATR-FTIR spectra, which revealed the high content of trans fatty acids in palm oil and cis fatty acids in blended oil, rendering it more suitable for consumption.This work could provide the basis for the large-scale production of antioxidant-rich nutritional oil blends using common vegetable oils with substantial oxidative stability as a substitute for palm oil in the food industry.

Lipidology 2024, 1 ,Figure 1 .
Figure 1.Pareto plots of transformed estimates for oxidative stability (%OS) and oxidative power (%OP) in various assays.Plot (A) is Rancimat, and plot (B) is DPPH assays, which measure %OS; plot (C) is the CD assay, and plot (D) is the CT assay, which assess %OP.A gold reference rectangle is drawn on the plot to indicate the significance level (p < 0.05).

Figure 1 .
Figure 1.Pareto plots of transformed estimates for oxidative stability (%OS) and oxidative power (%OP) in various assays.Plot (A) is Rancimat, and plot (B) is DPPH assays, which measure %OS; plot (C) is the CD assay, and plot (D) is the CT assay, which assess %OP.A gold reference rectangle is drawn on the plot to indicate the significance level (p < 0.05).

idology 2024, 1 ,
FOR PEER REVIEW correlation between sample X5 (i.e., corn oil) and DPPH • was observed, an outcome th could be attributed to the chemical composition of the specific oil sample.

Figure 2 .
Figure 2. Multiple factor analysis for the measured variables; the plot displays the factor scores each variable.Inset tables include variable loadings and colors representing the variables' block p rameters.

Figure 2 .
Figure 2. Multiple factor analysis for the measured variables; the plot displays the factor scores of each variable.Inset tables include variable loadings and colors representing the variables' block parameters.

gy 2024, 1 ,Figure 3 .
Figure 3. Multiple factor analysis for the measured variables; the plot displays the factor sco each variable and shows the block partial inertias.The inset table includes block partial inertia colors representing the variables' block parameters.

Figure 3 .
Figure 3. Multiple factor analysis for the measured variables; the plot displays the factor scores of each variable and shows the block partial inertias.The inset table includes block partial inertias and colors representing the variables' block parameters.

Figure 4 .
Figure 4. Multivariate correlation analysis of measured variable parameters.

Figure 4 .
Figure 4. Multivariate correlation analysis of measured variable parameters.

Figure 5 .
Figure 5.The PLS prediction profiler of each variable and the desirability function for optimizing blended vegetable oils are shown in plot (A), while the Variable Importance Plot (VIP) option graph with the VIP values for each predictor variable is shown in the plot (B) table.A blue dashed line in the VIT at 0.8 indicates the significance level of each variable.

Figure 5 .
Figure 5.The PLS prediction profiler of each variable and the desirability function for optimizing blended vegetable oils are shown in plot (A), while the Variable Importance Plot (VIP) option graph with the VIP values for each predictor variable is shown in the plot (B) table.A blue dashed line in the VIT at 0.8 indicates the significance level of each variable.

Figure 6 .
Figure 6.The relationship between the induction period and temperature through a polynomial regression fit.

Figure 6 .
Figure 6.The relationship between the induction period and temperature through a polynomial regression fit.

Figure 7 .
Figure 7. Stacked ATR-FTIR spectra of palm oil (orange line) and the optimal sample (green line).Blue values indicate the major wavenumbers in specific functional groups.

Figure 7 .
Figure 7. Stacked ATR-FTIR spectra of palm oil (orange line) and the optimal sample (green line).Blue values indicate the major wavenumbers in specific functional groups.

Table 1 .
The actual and coded levels of the independent variables were used to optimize the best blend of vegetable oils using the Main Effects Screening design.
(CDs)and triene (CTs) after storage; these values were excluded from RSM.

Table 5 .
Statistical analysis between the induction period and temperature for both palm oil and optimal sample.

Table 5 .
Statistical analysis between the induction period and temperature for both palm oil and optimal sample.

Table 6 .
Major wavenumbers in the ATR-FTIR spectra of palm oil and the optimal sample.
b 0.227 ± 0.017 a C-H stretching symmetric vibration of the cis double bonds, =CH 2922 1.421 ± 0.053 a 1.276 ± 0.043 b Asymmetric stretching vibration of C-H of aliphatic CH2 group 2853 1.076 ± 0.052 a 0.931 ± 0.020 b Symmetric stretching vibration of C-H of aliphatic

Table 6 .
Major wavenumbers in the ATR-FTIR spectra of palm oil and the optimal sample.Within each row, statistically significant differences (p < 0.05) are denoted with lowercase letters (e.g., a, b).