3.2. Statistical Analysis
Due to the extreme visual similarity between the
1H NMR spectra, previous exploratory analysis by PCA was performed to verify if the samples of edible oil-based nutraceuticals could be statistically differentiated. The PC1 x PC2 score chart indicated clear discrimination tendencies between the samples (
Figure 2A), with the positive sense of PC1, the statistical region responsible for discriminating the most significant number of samples. In contrast, the negative PC1 discriminated the soy, linseed, and chia samples. The accumulated variance explained by the first two components used in the construction of the graph was 86.9%.
Analyzing the PCA loadings (
Figure 2B), the principal chemical descriptors responsible for discriminating the samples presented in the score plot were identified. The positive scores of PC1, the statistical region where most of the sampling were discriminated, was strongly influenced by the signals of hydrogens with chemical shifts at
δ 0.88, assigned to the hydrogens of terminal methyl groups of oleic (ω-9) and SFAs. On the other hand, the trend observed along with the negative PC1 scores, the statistical region populated by soybean, linseed, and chia samples, was strongly influenced by the signals at
δ 5.38, 2.83, 2.08, and 0.98, assigned to vinylic, allylic, bis-allylic, and terminal methyl hydrogens of α-linolenic acid (ω-3), respectively.
Extrapolating the results from PCA, we can infer that the distinction between the oil-based nutraceuticals samples was strongly influenced by the different contents of MUFAs (ω-9) and PUFAs (ω-6 and ω-3) in addition to signals attributed to SFAs. Therefore, an adaptation between the equations described by Miyake et al. (1998) and the ERETIC method to calculate the relative percentages of the lipid descriptors (SFA, ω-9, ω-6, and ω-3) indicated as PCA discriminators for each nutraceutical sample evaluated in the present study [
31]. It is worth noting that the application of NMR in the quantification of fatty acids relies on other protocols, such as those described by Martinez-Yusta et al. (2014) and Santos et al. (2018) [
5,
32]. The equations used in the aforementioned calculations, as well as the integrations of the signal areas, are shown in
Figure 3.
Table 2 presents the results of the quantifications for each type of sample. Analyzing the percentages of ω-9 fatty acids and SFAs of the different matrices, it was found that nutraceuticals based on coconut, palm, andiroba, Brazil nut, and safflower oils, previously divided into positive PC1 scores (
Figure 2A) had low percentages of ω-3 fatty acids. On the other hand, these samples showed the highest levels of SFAs and ω-9 fatty acids. This information is consistent with the analysis of PCA loadings (
Figure 2B), which indicated as the main descriptor for the region, the signal in δ 0.88, previously assigned to terminal methyl groups of SFAs and ω-9 fatty acids.
The matrices with the highest percentages of ω-3 acids were the samples of chia and linseed. These samples were broken down into negative PC1 scores, as shown in
Figure 2A. Again, the descriptors (
δ 5.38, 2.83, 2.08, and 0.98) indicated in the loading plot (
Figure 2B), previously assigned to hydrogens of the
α-linolenic acid (ω-3) and others unsaturated chemical groups, support the observation. The large group in the center of the graph in
Figure 2A, composed of oil-based nutraceuticals sunflower, fish, soy, primrose, garlic, copaiba, and almond samples, showed intermediary levels of saturated and unsaturated fatty acids when compared with the other samples. No apparent correlation between the percentages of saturated and unsaturated fatty acids and the descriptors indicated in the loading graph was observed for the large group of samples centered on the PCA score graph (
Figure 2A). For comparative analytical purposes of the results from the NMR-PCA model presented, the profiles of oil-based nutraceuticals were also evaluated via the FTIR technique, which is recognized as a more accessible for the industry.
Figure 4 shows the score plot for PC1 x PC2, with accumulated variance explained by the first two components equal to 72.5%.
The data dispersion was similar to that obtained by the ¹H NMR data. Along the positive PC1 scores, safflower, palm, garlic, andiroba, and copaiba samples were discriminated, and the statistical region was strongly influenced by signals with wave numbers equal to 1417 cm-1 (C=CH), 1654 cm-1 (C =C), and 967 cm-1 (HC=CH). These signals are characteristic of absorptions (folding and stretching) of PUFAs. The coconut samples were discriminated in the negative scores of PC1, influenced mainly by signals referring to the stretching of C-O and C=O bonds, typical of ω-6 and ω-9 fatty acids, highlighting again the contents of these acids, previously determined by 1H NMR. As in the NMR analyses, chia and linseed samples were discriminated similarly. However, now it was along the negative PC2 scores, with the C=C link (ν 722 cm-1) as the main descriptor.
Regarding the similarity between the edible oils, the multivariate analysis of the FTIR data via HCA highlighted the remarkable similarity between the chia and linseed samples (
Figure 5), previously identified in the PCA treatment of the NMR data. The resulting dendrogram also indicated the existence of significant similarity between the samples of soybean, sunflower, and fish (pure and declared mixture). This similarity between the samples is the principal basis of the adulteration of fish oils, commonly practiced with these vegetable oils (soybean and sunflower) of lower market value. However, despite HCA presenting data in agreement with the PCA-NMR and PCA-FTIR models, the resulting dendrogram indicated a more remarkable similarity between copaiba and coconut oils, information not evidenced in the exploratory analyzes by PCA of NMR and FTIR data.
Although NMR allows the simultaneous obtaining of a qualitative (determination of constituents) and quantitative (relative lipid percentages) lipid screening, the FTIR technique showed concordant and rapid results, confirming analytical complementarity since similarities were observed in the discrimination of nutraceuticals by composition. The interpretation of the results can still be extrapolated in terms of clear distinctions between nutraceuticals regarding the origin of the oilseed matrices used, as well as peculiar results from the FTIR-HCA model that indicated the possibility of verifying the authenticity of the edible oils used as raw materials.