3.1. Oxidative Stability during Storage
The surface area of enzymes and substrates is increased by flour products, which further triggers widespread lipid oxidation and degrades the nutritional value and storage stability of quinoa flours. Therefore, we assessed the storage stability of quinoa flour by measuring FA contents, PVs, MDA values, and lipase activities during storage. Among them, FAs are usually used to measure the degree of lipid hydrolysis, PVs and MDA values are used to measure the degree of oxidation of quinoa flour, and lipase is the key enzyme in the process of lipid hydrolysis.
Lipase activity is an important factor contributing to the loss of quinoa quality. The changes in lipase activity of quinoa stored at various conditions are shown in
Figure 1A. It increased from 9.82 mg KOH/g at the beginning of storage to 14.83 (LL), 16.09 (NN), and 19.64 (HH) mg KOH/g at 120 days, respectively. It is clear that lipase activity increased with the prolongation of the storage time, indicating that quinoa lipids were continuously hydrolyzed during the storage process. At the same time, there were also significant differences between storage conditions, with lipase activity significantly higher (
p < 0.05) in the HH condition than on the 0th day from the 30th day onwards, whereas no significant change in lipase activity occurred in the LL and NN conditions at this time. Lipase activity also started to increase significantly at 60 days of storage under NN conditions, and it was not until 90 days of storage that lipase activity under LL conditions started to change significantly. These results suggested that LL has a delaying effect on the hydrolysis of fat during the storage of quinoa flour, suggesting that low temperatures inhibit enzyme activity, whereas an increase in lipase activity is induced with increasing storage temperatures, which in turn accelerates lipid hydrolysis and oxidation.
During storage, the fat in quinoa is hydrolyzed to fatty acids by the lipase enzymes, leading to an increase in the fatty acid. So, FA was a crucial storage indicator that also had an impact on lipid oxidation [
21]. As shown in
Figure 1B, the FA values during 120 d of storage under LL, NN, and HH conditions were determined. As the samples were exposed for longer times, a significant increase (
p < 0.05) in FA content was observed. It increased to more than three, five, and eight times the initial content at 120 days of storage, respectively. These results show that hydrolytic deterioration occurred in quinoa flour, which resulted in the production of free FAs. Compared with LL and NN conditions, the quinoa flour under the HH condition showed a significant increase in FAs (
p < 0.01). The fatty acid values of LL and NN conditions were 64 percent and 40 percent lower than those of HH conditions when stored for up to 120 days. It can be seen that the hydrolysis of fats in quinoa can be slowed down under a low-temperature storage environment.
PV is a measure of the primary oxidation of fats because hydroperoxides are the main compounds formed during fat oxidation [
22]. As shown in
Figure 1C, the PVs of the samples increased from 15.23 to 42.26, 71.95, and 100.12 mg/100 g after 120 d of treatment in LL, NN, and HH conditions, respectively. The peroxide values of HH conditions were 58% (
p < 0.01) and 28% (
p < 0.01) higher than those of LL and NN conditions, respectively. As expected, HH conditions had the highest level of PV compared to LL and NN conditions. Although the PV continued to increase after 60 d, there was a significant slowing down of the rate of increase compared to the previous period. This is mainly because the value of hydroperoxide is in dynamic equilibrium. As the oxidation reaction proceeds, the peroxide is further oxidized to secondary oxidation products, which affects the shelf-life of quinoa and its related products in turn.
Hydroperoxides are primary oxidation compounds, which will produce various secondary oxidation products (including aldehydes, alcohols, acids, and ketones) after decomposition, thus producing unpleasant tastes and smells. Malondialdehyde (MDA) is a stable and abundant secondary end product which is considered to be one of the most significant indicators of the state of lipid peroxidation [
23]. The accumulation of MDA indicates that lipids are oxidized and will cause damage to cells. Among them, the peroxidation of lipids significantly alters the physical properties of the lipid bilayer, which in turn affects its cytomembrane fluidity and permeability. And the accumulated MDA is thought to participate in the degradation of various biological functions by attaching to biomolecules such as proteins and nucleic acids [
24]. Therefore, MDA was also used to assess quinoa lipid oxidation during storage in this study. The MDA contents of all studied samples as an indication of fat oxidation in quinoa are presented in
Figure 1D. MDA values in LL and NN conditions were relatively lower during the first 30 days of storage, indicating that these flours did not initially display an oxidative rancidity. However, the MDA reached 80.66 μg/kg after 30 days of storage under HH conditions, which was significantly higher than that in LL and NN conditions (
p < 0.05). Meanwhile, the malondialdehyde content at the end of storage was only 71.32 μg/kg under LL conditions, which was significantly lower than that under NN and HH conditions (
p < 0.05), suggesting that high-temperature storage of quinoa flour increased the degree of membrane lipid peroxidation and the accumulation of MDA, which caused damage to cell membranes and cells. We also found that the changes in FAs and PVs were significantly greater than those in MDA throughout the storage process, which suggests that quinoa flour lipids were first hydrolyzed to produce fatty acids during storage, and then the fatty acids were further oxidized to produce primary and secondary oxidation products.
3.2. Global Lipidomics Analysis of Quinoa
Quinoa flour samples with three storage conditions at 0, 30, 60, 90, and 120 d were scanned three times by UPLC-Q-Exactive Orbitrap MS/MS in ESI+ and ESI− modes. QC samples were injected in every 12 samples and followed by two blank samples at the beginning and end of the analysis batch.
The representative base peak chromatogram of the sample on 0 day is shown in
Figure S1. There were 229 lipids found in Quinoa flour samples: 53 were found in the negative ion mode, and 176 were found in the positive ion mode. Quinoa lipids could be categorized into five categories by LIPID MAPS, including FA, GL, GP, SP, and ST. Furthermore, these five major classes could be further divided into 14 subcategories, including FA, TG, DG, PC, PE, PI, PS, PA, PG, LPC, LPE, Cer, Hex1Cer, and AcHexSiE (
Figure 2). The most abundant subcategory was Triglyceride (TG) (84 lipids), followed by diglyceride (DG) (26 lipids) and fatty acyls (FAs) (25 lipids). The above results illustrate the diversity and complexity of quinoa lipids. In addition to TG ([M+NH4]
+), DG ([M+NH4]
+), PI ([M+NH4]
+), and AcHexSiE ([M+NH4]
+), others were detected as [M+H]
+ or [M−H]
− adducts.
Table S2 provides a complete list of 229 lipids with their lipid names, theoretical
m/
z, detection
m/
z, RT information, and mass accuracy.
3.3. Quantification and Changes in Quinoa Lipids during Storage
In this study, an internal standard for isotopes was used to compare quantitative results. The sum of individual lipid contents was used to analyze and compare the content of similar or subclasses of lipids in quinoa flours, and heat maps (
Figure 3A,B) were used to visualize the dynamic changes in quinoa lipid contents at various storage periods.
The most abundant lipids in quinoa were GLs and GPs, followed by SPs, STs, and FAs (
Figure 3A). There were significant differences in various classes of lipids between unoxidized and oxidized quinoa samples. In the HH condition, the content of FAs was significantly increased from 1120.39 to 4432.24 nmol/g, and the contents of GLs and GPs were decreased from 320,040.07 and 8334.76 to 74,954.70 and 283.61 nmol/g, respectively. However, these differences were not significantly observed in LL and NN conditions. Lipid hydrolysis during quinoa storage is highly dependent on lipase and phospholipase [
25]. TG is hydrolyzed to FA and glycerol by lipase and then releases DG or MG, which continue to be hydrolyzed. And GPs are more susceptible to oxidation because of more unsaturated bonds [
26]. Therefore, FAs displayed a fluctuating upward trend, while GLs and GPs displayed a fluctuating downward trend during the storage period.
Next, we used a heat map to show the differences in lipid subclasses between fresh and stored quinoa flours at various storage stages (
Figure 3B). The lipid profiles of quinoa flours change continuously during storage. Slight changes in the lipid composition of quinoa during the initial storage period may be attributed to the important role of antioxidant components such as superoxide dismutase and carotenoids in inhibiting quinoa lipid oxidation [
27]. Later, with the increase in storage temperature and storage time, more changes took place in lipid compositions, and a significant change in lipid metabolite content was observed when stored in the HH condition for 120 days. This is due to the continuous involvement of lipids in quinoa flour in the corresponding metabolic reaction. Additionally, the oxidation process intensified with the degradation of antioxidant components in quinoa flour. The oxidation rate was also accelerated by lipid oxidation products with high oxidation capacities, such as free radicals and hydroperoxides [
28].
TG exudates from the damaged cell membranes during pre-treatment processes, such as the shelling and milling of quinoa, which is further hydrolyzed by contact with previously dormant but highly active lipases in the endosperm layer and germ tissue [
29]. We found 89 species of TG in quinoa in this study. The level of TG fluctuated up and down under LL and NN conditions, but finally decreased from 233,301.93 nmol/g to 203,824.95 nmol/g and 163,745.94 nmol/g, respectively. The decrease in TG content during storage was attributed to hydrolysis by lipase. And the level of TG increased, likely due to the action of glycerol-phosphate acyltransferase (the rate-limiting enzyme of TAG/GL synthesis) [
30]. However, the content of TGs in quinoa flour under HH conditions was directly reduced to 59,573.33 nmol/g. The above results indicate that there was higher lipase activity under HH conditions, which increases the hydrolysis of TG. The FAs, produced by hydrolysis, are further oxidized to produce substances with odors as well as reducing the nutritional value of quinoa.
DG is one of the major lipid subclasses in organisms and serves as the second messenger for a variety of cellular activities, as well as participating in the biosynthesis of phospholipids such as PC, PE, and PS [
20]. In this study, 26 species of DG were detected in quinoa. These changes followed a pattern that closely resembled that of TG. The results showed that the level of DG fluctuated up and down under LL and NN conditions, but finally decreased from 86,738.14 nmol/g to 54,125.28 nmol/g and 44,500.77 nmol/g, respectively. However, the content of DG in quinoa flour under HH conditions was directly reduced to 15,381.37 nmol/g. It is possible that DG accumulates due to the degradation of TG, whereas the decrease is due to hydrolysis by lipase.
Glycerophospholipids mainly include PAs, PCs, PEs, PGs, PSs, and PIs, which play a key role in maintaining cellular homeostasis and are the main structural components of cell membranes [
31]. In total, 5 LPCs, 5 LPEs, 3 PAs, 20 PCs, 16 PEs, 10 PGs, 10 PSs, and 12 PIs were identified in the study. The abundance of PCs, PEs, PSs, PIs, and PGs was significantly reduced under HH conditions, which may be associated with phospholipase A in quinoa. Phospholipase A, which consists of phospholipase A1 and phospholipase A2, hydrolyzes the respective acyl chains to generate FA and lysophospholipids (LPLs) [
32]. Similarly, these phenomena were obvious under HH conditions, but not obvious under LL and NN conditions. In conclusion, these findings show that quinoa’s lipid composition and content can be greatly impacted by storage conditions.
Ceramides and hexosylceramides are components of sphingolipids and glycosphingolipids, which play an important role in regulating cell differentiation, proliferation, apoptosis, aging, and other life activities [
33]. And these molecules are important structural compounds that make up the lipid bilayer of the cell membranes. In particular, ceramide is involved in mitochondrial-mediated apoptosis [
34]. The findings indicate a significant decrease in the content of Cer and Hex1Cer during storage under the HH condition, which indicates that storage conditions had a key influence on the degradation of Cer and Hex1Cer in quinoa. It is possible that the degradation of Cer may lead to mitochondrial dysfunction and apoptosis, which will affect the storage quality of quinoa.
We further analyzed the fatty acid composition of quinoa and identified 25 fatty acids in total. The unsaturated fatty acids of fresh quinoa flours were mainly composed of C18:1 oleic acid (17.68%), C18:2 linoleic acid (47.88%), and C18:3 linolenic acid (8.82%). Moreover, the proportion of saturated fatty acids was relatively low, and C16:0 palmitic acid (6.70%) was the main component (
Figure 3C). FAs increased significantly (
p < 0.05) in all quinoa samples during storage with the degradation of glycerolipids and glycerophospholipids (
Figure 3D). However, the percentage of linoleic acid C18:2 was slightly reduced after accelerated storage for 120 d, and the percentages of palmitic acid C16:0 increased. A decrease in PUFAs and an increase in saturated fats (SFAs) were observed after quinoa was oxidized, which was consistent with the observations of Sun et al. [
35]. An epidemiological study showed that intake of SFAs increased the risk of cardiovascular disease (CVD), but the intake of more PUFAs could slow down CVD [
36]. The increase in SFAs and decrease in PUFAs were not beneficial to human health.
3.4. Significantly Different Lipids between Fresh and Stored Quinoa
To further explore changes in lipid composition, we analyzed and visualized lipidomics data during quinoa storage by multivariate statistical analysis. We performed PCA analysis of the lipidomic data to understand the metabolic differences among groups and the variations in the samples within groups, and explored the main factors contributing to differences during storage.
Figure 4A shows the PCA score plot of quinoa lipid under different storage conditions. The distinction between 0 and HH conditions was clear, but the distances between LL and NN conditions were not clear. These findings indicated significant lipidome changes in HH conditions, and lipid oxidation was not significant under LL and NN conditions. However, the differences between samples of different oxidation stages were similar and could not be well distinguished. Therefore, the differences in quinoa lipid metabolism during different storage periods were additionally analyzed by PLS-DA, which resolved variations in the data more effectively than PCA. The PLS-DA score plot shows differences in lipid composition during quinoa oxidation. A clear separation was observed between unoxidized samples and other samples in HH conditions, but the differences among samples oxidized for 30 d, 60 d, 90 d, and 120 d in LL and NN conditions were small (
Figure 4B). The score plot of the PLS-DA mode shows a good interpretation rate (R
2X = 0.978 and R
2Y = 0.945) and better predictability (Q
2 = 0.864). Moreover, the robustness and predictability of the model were supported by the results of 200 permutation tests (
Figure 4C). The intercepts of Q
2 are expected to be below 0.00. According to the results, the PLS-DA model was reliable, predictive, and not overfit. The results confirmed that quinoa lipidomics profiles were significantly altered during storage in HH conditions.
To further elucidate the metabolic changes in quinoa during storage, we tried to explore the important differential lipid metabolites and associated pathways. Based on the PLS-DA model described above, we used the threshold
VIP > 1 and
p < 0.05 to screen for significantly different lipids during quinoa oxidation. The changes in these differential lipids during storage were visualized by the heatmap shown in
Figure S2. A total of 90 significantly different lipids, including 40 TGs, 21 DGs, 2 PGs, 3 PSs, 6 PEs, 9 PIs, 8 FAs, and 1 HexlCer, were identified in PLS-DA. As a result of oxidation, eight FAs were significantly up-regulated, including FA (18:1) and FA (18:2). All of the remaining significantly different lipids were down-regulated, with PS (18:2_18:2) having the highest fold-change value (3.74), followed by PS (18:0_19:1) (3.47), PS (18:0_19:0) (3.47), and TG (4:0_16:1_18:2) (3.39) (
Table S1). At the same time, the loading plot of PLS-DA (
Figure 4D) shows that the lipid species of the TG and DG classes were the most important variables contributing to the discrimination between quinoa flour, which can be used to distinguish fresh quinoa flour from oxidized quinoa flour. The lipid content of quinoa changes during storage, which may have an impact on its composition, nutritional value, and function. Thus, it is necessary to investigate the specific functions and regulatory mechanisms of specific lipids during oxidation.
3.5. Correlation Analysis between Lipase Activity and Lipids
We assessed the relationship between lipase activity and different subclasses of lipid content in quinoa during storage. The findings revealed a highly significant positive connection between lipase activity and FAs (R = 0.66) (
Figure 5A). However, it was negatively correlated with most other lipids, especially TGs (R = −0.82) and DGs (R = −0.79). These results also indicated that TGs and DGs were hydrolyzed by lipase to produce FAs during the storage of quinoa flour. At the same time, we found that there were significant correlations between the different lipids in quinoa. For example, highly significant correlations were found between DGs and PCs (R = 0.80), PEs (R = 0.84), and PSs (R = 0.82), with some studies suggesting that DGs were involved in their synthesis [
37]. Doblado-Maldonado et al. also found that the lipase activity was crucial to the deterioration of grain quality. It not only affects the lipid hydrolysis rancidity by lipase action [
38], but also affects subsequent fatty acid degradation and oxidation. Therefore, enzyme inactivation during the pre-treatment of quinoa would contribute to its oxidative stability. Ling et al. used radio-frequency (RF) heating treatment to inactivate lipase activity in wheat germ and found that it could improve the physicochemical properties of wheat germ and enhance its lipid stability [
39].
We conducted correlation heat map analysis on 90 lipids with substantial differences to discover the relationships between various lipids, and Pearson correlation coefficients were used to determine the correlation magnitude, which ranges from −1 to 1. Correlation size was assessed by the size of the circle, and negative or positive correlations were indicated by values closer to −1 or 1, respectively. We found that the majority of lipids were highly correlated with each other. Positive correlations were observed among lipids in the same subclasses and were stronger than the correlations between subclasses. As shown in
Figure 5B, significant positive correlations were found between TGs and DGs as well as between PEs, PSs, PGs, and PIs, whereas they were all negatively correlated with FA. These were because the metabolism of both glycerolipids and glycerophospholipids produces FAs. Therefore, there was a complex relationship between significantly different lipids. The different lipids also interact with each other, which in turn affects the oxidation process of quinoa.
3.6. Lipid Metabolism Pathways during Quinoa Storage
Metabolic pathway analysis is crucial for understanding the effect of storage time on the metabolic processes related to the lipid profiles of quinoa. In order to explore the lipid metabolic pathways associated with quinoa during storage, we mapped 90 species of significantly different lipids to the KEGG database to obtain the ID information. Moreover, MetaboAnalyst 5.0 was used to analyze metabolic pathways topologically. These significant lipids were mainly involved in 12 metabolic pathways, including the biosynthesis of fatty acids; the biosynthesis of unsaturated fatty acids; linoleic acid metabolism; glycerophospholipid metabolism; cutin, suberine, and wax biosynthesis; sphingolipid metabolism; GPI-anchor biosynthesis; glycerolipid metabolism; fatty acid elongation; inositol phosphate metabolism; the phosphatidylinositol signaling system; and fatty acid degradation. The results are shown as a bubble chart, and each bubble represents a metabolic pathway. Bubble size represents the enrichment value of the pathway and bubble color represent the influence value. As presented in
Figure 6A, these significantly different lipids were mainly enriched for linoleic acid metabolism, glycerophospholipid metabolism, and keratin, cork, and wax biosynthesis. PEs, PIs, PSs, PGs, and FAs were involved in the three lipid metabolism pathways above.
Based on KEGG pathway analysis,
Figure 6B shows the possible transformation of quinoa lipid molecules during storage. In the presence of phospholipase and lipase, hydrolyzed GP and GL could result in an increase in FAs, including FA (16:0), FA (18:1), and FA (18:2). During the hydrolysis of GL, TG, and DG, the major identified lipids is TG. At the same time, diacylglycerol acyltransferase (DGAT) could catalyze the covalent binding of DG with fatty acyl coenzyme A to produce TG. As mentioned in 3.3 above, the contents of TG and DG fluctuated up and down during the storage of quinoa. During the storage of quinoa, major lipids involved in the glycerophospholipid pathway include PE, PI, PS, and PG. DG could be involved in the synthesis of PE and PC as a lipid anchor in the presence of the corresponding phosphotransferases, while phospholipase C catalyzes the hydrolysis of PC and PE for the synthesis of DG [
40]. In addition, lysophospholipid (LPC and LPE) could be generated via the cleavage of glycerophospholipid (PC and PE) molecules by phospholipase A1 or phospholipase A2, whereas lysophospholipid and acyl-coenzyme A could form glycerophospholipid in the presence of lysophospholipid acyltransferase. Among them, phospholipases play a critical role in the metabolism of phospholipids and the signaling of cells [
41]. In the linoleic acid metabolism pathway, stored quinoa had a higher level of linoleic acid than fresh quinoa. GP and GL were degraded to produce linoleic acid (LA), which was subsequently oxidized to form primary oxidation products (ROOHs) under enzymatic and photosensitive oxidation involving lipoxygenase (LOX) and singlet oxygen (
1O
2), including 9-ROOH and 13-ROOH. In the auto-oxidation process, the methylene group of the double bond of unsaturated fatty acids removes an H and generates free radicals (R·). After entering the chain transfer stage, R· will generate peroxide radicals (ROO·) with the participation of oxygen due to instability, and ROO· will continue to react with LAs or other unsaturated fatty acids to generate new radicals (R·) and hydroperoxides (ROOHs) [
42]. However, these hydroperoxides are unstable and can be further oxidized to form secondary oxidation products such as aldehydes, alcohols, acids, and esters (R-CHOs, R-OHs, R-COOHs, and R-COO-R’s), which in turn produce an irritating odor that affects the edibleness of quinoa. During the storage of quinoa, lipid oxidation is a significant reaction that leads to spoilage and deterioration. Nevertheless, the identification of these oxidation byproducts was not achieved in the present investigation. Meanwhile, the upregulation of FAs during quinoa storage leads to disturbances in the biosynthetic metabolism of keratins, spermine, and waxes, which reduces the protective effect of these substances on quinoa. The above findings in this research further confirmed the ideas in the relevant literature.
In conclusion, significant changes in lipids occurred during the storage of quinoa, involving 12 metabolic pathways and various enzymes. KEGG pathway analysis of these significantly different lipids in fresh and oxidized quinoa is helpful to explore the metabolic process of quinoa lipids during storage. Therefore, it is necessary to further study the metabolic pathways and enzymes related to this process. For example, it is necessary to identify the oxidation and degradation of FA products in order to provide a deeper understanding of these lipid metabolic pathways and reveal the functions of associated enzymes.