Widely Targeted UHPLC-MS/MS Metabolomic Analysis on the Chemical Variation in Blueberry-Filled Pastries During Processing

Introduction

Blueberries contain various phytochemicals with potential health effects, including anthocyanins, flavonols, vitamins, lutein, sugars, and organic acid (1–3). The consumption of polyphenol-rich blueberries and their polyphenol extracts is associated with the protection and management of non-communicable diseases such as cancers, cardiovascular diseases, diabetes, osteoporosis, and neurodegenerative diseases (1). However, the perishable nature and seasonal availability of berries limit the all-year-round supply and consumption. Food processing allows the all-year-round provision and provides various types of blueberry products to compensate for the different requirements and appreciations by the customers. One large group of products are bakery food products filled with blueberries or blueberry containing materials, such as pies, pastries, muffins, cakes, and cookies. Such food products are becoming increasingly popular and are consumed daily in many areas of the world. Processing can prolong the shelf life of blueberry products and increase the diversity of products for customers' choice; however, nutritional compounds, such as polyphenols, may be destroyed to some extent during processing, especially during thermal processing, which in turn affects the health-promoting properties of the products (4–6).

However, knowledge of the effect of thermal processing on the phytochemical profile of the products is limited to a few compounds due to the conventional analytical methods focusing only on a certain group of chemicals (3, 7, 8), although the biological importance is determined by the whole phytochemical profile. Thus, a technique should be developed to fully elucidate the components in a food product for evaluating the nutritional changes during thermal processing. Metabolomic analysis is used to provide a comprehensive chemical profile of as many small molecules as possible in a system (plant, cell, or food material), which allows for a more thorough and encompassing analysis of molecular composition than traditional methods. The application of metabolomic study on the tracking of dynamic compositional changes during food processing could provide a holistic view on the chemical changes and correlated reactions under various treatments of food. Thus, in recent years, metabolomic analysis has been applied in food products, such as in the investigation of nectar beverage of black raspberries (9), roasted coffee beans (10), wines (11), and extra-virgin olive oil (12).

In our previous study (8), the variation in individual flavonoids, including flavonol glycosides and anthocyanins, was investigated in two thermal processing procedures of blueberry-filled pastries. The results showed that stir-frying largely decreased the contents of flavonoids, especially anthocyanins, in blueberry filling, and baking caused less variation in flavonoid profiles. However, the changes in other phenolic compounds, which may be related to the degradation of flavonoids, and the components aside from polyphenols in the fillings remain largely undiscovered. The chemical changes and the underlying mechanisms during food production should be further studied. Food matrix is such a complex mixture that phenolic compounds, with their chemical reactivity, may undergo reactions with various compounds other than their own degradation. Therefore, a holistic view on the variation of phytochemicals in blueberries during processing is important to provide insights into the changes in bioactive compounds, and encourage and give hints for scientists to further evaluate the underlying chemical reaction mechanisms of key varying compounds. Such information on the effect of thermal processing on the phytochemical profiles in real food systems, such as blueberry pastry as in this case, remains largely limited.

The current study aims to provide useful information on the chemical changes in blueberry-containing products during thermal processing by using a widely targeted metabolomic approach. The results will present clues on the degradation of thermal-sensitive components and the interaction between bioactive compounds of the products during processing. The work also provides further insights into the utilization of widely targeted metabolomic techniques in the screening of key chemical changes and the revelation of underlying chemical mechanisms.

Materials and Methods

Materials and Chemicals

Fresh blueberries (Vaccinium spp.) were imported from Chile, and all the other ingredients of blueberry-filled pastries were purchased from the local market. Methanol, acetonitrile, and acetic acid were of HPLC grade and purchased from Merck (Darmstadt, Hesse, Germany). Double deionized ultrapure water was prepared using an ultrapure-water purification system (Millipore, Bedford, MA, USA). Chemical standards were purchased from Sigma–Aldrich (St. Louis, MO, USA) and BioBioPha (Kunming, Yunnan, China).

Sample Preparation for Metabolomics Analysis

Samples from the same batch of blueberry pastry filling samples prepared according to our previous study (8) were subjected to metabolomic analysis for the understanding of phytochemical changes under different processing procedures, including the raw filling, stir-fried, and baked samples (170°C, 23 min). All samples were lyophilized and crushed using a mixer mill (MM 400, Retsch, Haan, Germany) with a zirconia bead for 90 s at 30 Hz. Exactly 100 mg sample was then weighed accurately, extracted with 1 mL of 70% aqueous methanol overnight at 4°C by vigorously vortexing for three times. After centrifugation at 10,000 × g for 10 min, the extracts were filtered (0.22 μm) into a vial for UHPLC-MS/MS analysis. The quality control (QC) samples were prepared in triplicate by evenly mixing all the sample extracts tested and injected every six sample injections to monitor the measure repeatability.

UHPLC-MS/MS Analysis

Exactly 2 μL of sample were injected and analyzed using a Shim-pack UFLC SHIMADZU CBM30A system (Shimadzu Corporation, Kyoto, Japan) interfaced to an Applied Biosystems 6500 QTRAP mass spectrometer (AB Sciex, Foster, CA, USA). The sample was separated on a Waters ACQUITY UPLC HSS T3 C18 column (100 mm × 2.1 mm, particle size 1.8 μm) (Waters Corporation, Milford, MA, USA) by using water (0.04% acetic acid) as solvent A and acetonitrile (0.04% acetic acid) as solvent B. The eluting gradient programme was as follows: 0.0–11.0 min, 5%−95% B; 11.0–12.0 min, 95% B; 12.0–12.1 min, 95%−5% B; and 12.1–15.0 min, 5% B. The flow rate was 0.4 mL/min. The temperature of column oven was kept at 40°C during the analysis. The effluent was connected to an ESI-triple quadrupole-linear ion trap (QQQ-LIT) mass spectrometer equipped with an ESI Turbo Ion-Spray interface operating in both positive and negative ion modes and controlled by Analyst 1.6.3 software. The ESI source operation parameters were as follows: source temperature 500°C; ion spray voltage (IS), 5500 V; and ion source gas I (GSI)/gas II (GSII)/curtain gas (CUR), 55/60/25 psi. Instrument tuning and mass calibration were performed with 10 and 100 μmol/L polypropylene glycol solutions in QQQ and LIT modes, respectively. QQQ scans were acquired as multiple reaction monitoring (MRM) experiments with collision gas set to 5 psi. Declustering potential (DP) and collision energy (CE) for individual MRM transitions were obtained with further DP and CE optimisation. A specific set of MRM transitions was monitored for each period according to the chemicals eluted within this period.

Qualification and Quantification of Chemicals

For qualification of the chemicals, the primary and secondary mass spectrometry data were analyzed using MVDB V2.0 Database (Metware Biotechnology Co., Ltd., Wuhan, China) and the public database of metabolite information. The mass error comparing with the database was limited within ±0.4 Da, and the error for retention time was allowed within ±0.1 min. The interference from isotope signals, duplicate signals of K⁺, Na⁺, and ${\text{NH}}_4^{+}$ ions, and duplicate signals of fragment ions derived from other larger molecules were excluded. Structural analysis of chemicals was conducted in reference to the existing mass spectrometry databases, such as MassBank (http://www.massbank.jp), KNAPSAcK (http://kanaya.naist.jp/KNApSAcK), HMDB (http://www.hmdb.ca), and METLIN (http://metlin.scripps.edu/index.php). The chemicals were quantified in the MRM mode using QQQ mass spectrometry. After chemical data from different samples were obtained, the peak areas of the mass spectra of all substances were integrated, and the mass spectrum peaks of the same chemical in different samples were subjected to integration correction.

Statistical Analysis

All samples were analyzed in triplicate. To investigate the effect of different thermal processing procedures on the composition of chemicals, hierarchical cluster analysis (HCA), principal component analysis (PCA), partial least squares-discriminant analysis (PLS-DA), and univariate analysis (UVA) on different samples were conducted with the MetaboAnalyst 4.0 platform. The thresholds for screening differential chemicals were set as follows: fold change ≥ 2 or ≤ 0.5 and VIP ≥ 1.

Results and Discussion

Our previous study (8) showed that the flavonoid content varied considerably during the thermal processing of blueberry-filled pastries, and the procedures of stir-frying and baking contributed to the different levels of the loss of anthocyanins and flavonol glycosides. However, the degradation products of these flavonoids and the variation of the other phytochemicals were left undetected. Pastry filling is a complex food mixture that contains many components other than flavonoids, and various reactions may occur in addition to the simple degradation of flavonoids. A holistic view on the phytochemical changes during the production of blueberry-filled pastries will help us evaluate the nutritional changes during thermal processing. The metabolomic approach has been largely applied to investigate the metabolite changes during plant development, such as fruit ripening and seed germination (13–15).

This work applied metabolomic analysis to investigate the phytochemical changes of blueberry-filled pastry products during different thermal processing stages. All the samples were lyophilized and analyzed on the basis of dry weight mass to eliminate the water deviation between the samples in different thermal processing procedures.

Full List of Chemical Components in Raw and Thermal-Treated Samples

A total ion current chromatogram of QC sample (an even mixture of all the samples involved in this study) and a multi-peak detection plot of chemicals in the MRM mode of the same sample are illustrated in Supplementary Figure 1. The chemicals in the samples were qualified and quantified by mass spectrometry by using the MVDB V2.0 database. The appearance and abundance of different substances in various samples were evaluated by MRM analysis. As shown in Supplementary Figures 1C,D, each peak displaying a distinct color represents one compound in the sample analyzed.

A total of 630 compounds were identified in the current study in the blueberry fillings of pastries collected at different processed stages, including 148 flavonoids, 21 carbohydrates, 69 organic acids, 24 quinates and their derivatives, 19 benzoic acid derivatives, 17 coumarins, 29 hydroxycinnamoyl derivatives, 5 proanthocyanidins, 2 flavonolignans, 28 amino acids and 56 of their derivatives, 18 fatty acids, 16 glycerolipids, 35 glycerophospholipids, 4 nicotinic acid derivatives, 57 nucleotides and their derivatives, 16 phenolamides, 7 indole derivatives, 3 pyridine derivatives, 4 cholines, 7 alcohols and polyols, 7 alkaloids, 7 tryptamine derivatives, 3 terpenoids, and 28 other chemicals. Amongst them, the largest group was flavonoid, which consisted of 13 anthocyanins, 10 catechin derivatives, 17 flavanones, 37 flavones, 26 flavone C-glycosides, 38 flavonols, and 7 isoflavones. The name and molecular formula of compounds, their precursor and product ions for qualification and quantification with MRM, and the peak integration values of the compounds are all listed in Supplementary Table 1. To the best of our knowledge, this study was the first to utilize the targeted UHPLC-MS/MS-based metabolomic approach in analyzing the chemical variations in fruit fillings of bakery products under different thermal treatments. Such a large amount of the compositional information provides a comprehensive understanding of the chemical changes in berries commonly enriched with bioactive polyphenols when subjected to thermal processing.

Difference in the Composition Amongst Raw, Stir-Fried, and Baked Blueberry Fillings

PCA and HCA were applied to analyze the data obtained for a full view of the differences between the blueberry filling samples obtained at different processing stages. The PCA score scatter plot (Figure 1) showed that samples obtained at different processing stages could be well-distinguished. The QC samples, as a mixture of all the three samples investigated, were located in the middle of the plot. The PC1 explained 56.11% of the variances of chemical data of all the samples investigated and predominantly distinguished the raw filling samples from the thermally treated samples, such as stir-fried and baked samples. PC2 explained only 13.51% of the variances of chemical data and distinguished mainly the stir-fried filling samples from the baked filling samples. This finding indicated that stir-frying and baking dramatically changed the composition of blueberry fillings. This finding was further confirmed by the heatmap of the hierarchical clustering analysis of differential chemicals amongst samples of raw, fried, and baked fillings (Figure 2). A total of 288 differential chemicals were screened by UVA. The HCA of the 288 differential chemicals in the three samples showed clear grouping patterns. The color sequence from red to green indicates a decrease in chemical contents. The abscissa of the thermogram indicates the different sample groups, and the ordinate indicates the differential chemicals identified between groups. According to the Venn diagram shown in Figure 3, 197 and 238 chemicals significantly differed (fold change ≤ 0.5 or ≥ 2, and VIP ≥ 1) between raw and fried fillings and between raw and baked fillings, respectively. However, the number of differential chemicals in fried filling vs. baked filling was only 75, which confirmed our previous finding that stir-frying was the major step contributing to the most variation in chemicals during blueberry pastry preparation (8). Consistently, the heatmap of HCA showed larger differences in abundance between raw and fried fillings than those between baked and fried fillings (Figure 2).

FIGURE 1

Figure 1. PCA scores scatter plot showing differences in composition between blueberry filling samples obtained at different thermal processing stages. mix, quality control sample; raw, fried, and baked refer to raw, fried, and baked fillings, respectively.

FIGURE 2

Figure 2. Heatmap of the hierarchical clustering analysis of differential chemicals amongst samples of raw, fried, and baked fillings. raw, fried, and baked refer to raw, fried, and baked fillings, respectively.

FIGURE 3

Figure 3. Venn diagram of the differential compounds of blueberry fillings at different processing stages. raw, fried, and baked refer to raw, fried, and baked fillings, respectively.

Amongst 288 differential compounds, 117 belonged to the group of phenolic compounds (Table 1), including 18 flavones (out of 37 flavones detected), 13 anthocyanins (13), 21 flavonols (38), 12 flavone C-glycosides (26), 8 flavanones (17), 14 benzoic acid derivatives (19), 14 hydroxycinnamoyl derivatives (29), 4 catechin derivatives (10), 5 coumarins (17), 6 quinates and their derivatives (24), 1 isoflavone (7), and 1 proanthocyanidin (5). The other compounds (Table 2) included 34 nucleotides and their derivatives (57), 29 glycerophospholipids (35), 21 organic acids (69), 5 amino acids (28), 22 amino acid derivatives (56), 13 glycerolipids (16), 8 carbohydrates (21), 8 phenolamides (16), 6 fatty acids (18), 3 tryptamine derivatives (7), 3 indole derivatives (7), 2 alkaloids (7), 2 terpenoids (3), 2 nicotinic acid derivatives (4), 2 alcohols and polyols (7), 1 pyridine derivative (3), and 1 choline (4). However, the flavonolignans detected in the fillings did not vary during the processing of pastries.

TABLE 1

Table 1. List of differential chemicals belonging to phenolic compounds and their variations under different thermal processing procedures.

TABLE 2

Table 2. List of differential chemicals other than phenolic compounds and their variations under different thermal processing procedures.

Variation in the Phenolic Composition of Blueberry Fillings During Thermal Processing

All the anthocyanins showed decreasing trends during the preparation of blueberry-filled pastries, and the most variation occurred at the stir-frying stage. This result was in accordance with our previous finding (8) in which anthocyanins were the most sensitive toward stir-frying and less sensitive to baking, which should be most likely attributed to the higher temperature and longer period of heating occurred during stir-firing than baking. Moreover, the low concentration of oxygen in the sealed fillings might decrease the reduction in flavonoids under baking conditions (16), while abundant air perfused into the fillings and accelerated the degradation of flavonoids during stir-frying.

Except for cyanidin 3-O-malonylhexoside, the contents of all the anthocyanins decreased by 52–100% during stir-frying in the current study. However, during baking, cyanidin 3-O-malonylhexoside along with four other anthocyanins, namely cyanidin 3-O-galactoside, malvidin 3-O-galactoside, delphinidin 3-O-galactoside, and petunidin 3-O-galactoside, decreased by 52–67%. Cyanidin 3,5-O-diglucoside was neither detected in the fried filling nor in the baked filling. Therefore, it was totally degraded after thermal treatment. Cyanidin 3,5-O-diglucoside was identified by Hou et al. (17) in the anthocyanin-rich extract from black rice and was found to be comparable to the other anthocyanins detected in terms of thermal sensitivity because of the low values of activation energies (E_a). Cevallos-Casals and Cisneros-Zevallos (18) suggested that sugars and their degradation products might accelerate the degradation of anthocyanins because of the association of the degradation rate of anthocyanin with the degradation rate of sugars to furfural-type compounds derived from the Maillard reaction. 2-Furoic acid, which is an oxidation product of furfural that originated from the decomposition of sugars, and ascorbic acid (19, 20) was found to be produced in the fillings after thermal processing (Table 1). Moreover, Louarme and Billaud (21) demonstrated that the production of 2-furoic acid depended on the oxidative degradation reactions of ascorbic acid rather than the thermal degradations, which might explain the remarkable occurrence of the compound during stir-frying (high oxygen exposure and long heating period) but mild changes under baking (low oxygen exposure in the sealed fillings and short period of heating).

Some flavonol aglycones detected in the current study increased considerably after thermal treatment by stir-frying, accompanying the decreases in the contents of their glycosides (Table 1). For example, the abundance of quercetin increased by 2.5 times, whereas some of its glycosides, such as quercetin-3,4′-O-diglucoside, quercetin 7-O-malonylhexosyl-hexoside, and quercetin 3-arabinoside, decreased remarkably after stir-frying. Kaempferol glycosides, such as kaempferol 3-O-galactoside, kaempferol 3-O-rutinoside, kaempferol 3-O-robinobioside, kaempferol 3,7-dirhamnoside, and kaempferol-3-O-robinoside-7-O-rhamnoside, decreased after thermal processing, whereas the aglycone kaempferol increased. The decreasing trends of the flavonol glycosides could be attributed to the deglycosylation that occurred during thermal processing as discussed in our previous study (8).

Changes in other phenolic compounds in blueberry filling during pastry preparation can be also determined by the widely targeted metabolomic analysis with UHPLC-MS/MS. For example, protocatechuic acid in the blueberry filling increased by four times after stir-frying. Protocatechuic acid was reported to be thermally degraded from cyanidin by the liberation of the catecholic B-ring or from flavonoids such as quercetin (9). It shared the same ortho-dihydroxyphenyl chemical moiety with cyanidin, which was reported to be critical for the biological properties of anthocyanins (22, 23). Thus, the degradation of cyanidin to protocatechuic acid might not result in the reduction of bioactivity. Protocatechuic acid was demonstrated to possess chemopreventive activity against several different types of cancers in animal and cell studies (24, 25). Coumarins are another group of phenolic compounds that possess various biological and therapeutic properties such as anti-oxidant, anti-microbial, anti-viral, anti-diabetic, anti-coagulant, estrogenic, vasodilator, anti-convulsant, anti-inflammatory, anti-hypertensive, and anti-cancer activities (26). In blueberry fillings, the amount of 6-methoxy-7,8-dihydroxycoumarin increased remarkably by 13 times after stir-frying, and by 34 times by the end of thermal processing. The degradation of anthocyanins might contribute to the increase (27). Moreover, 4-hydroxybenzoic acid was also reported to be the thermal degradation product of anthocyanin pelargonidin-3-glucoside (28). Although its content did not increase significantly in the filling after stir-frying, the abundance of pelargonidin-3-glucoside decreased remarkably. However, the final abundance of the compound increased significantly in the baked product compared with the raw filling. Other benzoic acid derivatives, such as gallic acid, vanillin, syringaldehyde, methyl gallate, 2,5-dihydroxybenzoic acid, 4-hydroxybenzaldehyde, anthranilic acid, and 2,3-dihydroxybenzoic acid, showed elevations in the filling after stir-frying. The increases in the amounts of vanillin, syringaldehyde, phenol, and phenol derivatives during heat treatment might be attributed to the thermal degradation of lignin (29). The thermal decomposition of the lignin polymer commonly started with the cleavage of the α-ether and β-ether bonds. As a result, a mixture of phenol-, guaiacyl-, and syringyl-type derivatives with their substituents in the aromatic ring was released (30); the attack of oxygen might lead to the formation of aromatic aldehydes and ketones (31). Amongst the hydroxycinnamoyl derivatives detected, 12 varied under the stir-frying treatment. Ferulic acid, cinnamic acid, p-coumaric acid, and sinapic acid showed increasing trends similar to that of benzoic acid derivatives. Moreover, the catechin derivatives, protocatechuic acid, protocatechuic aldehyde, and epigallocatechin showed increasing trends during stir-frying. The increases in these phenolic acids might be explained by the thermal degradation of flavonoids and the release of bound phenolic compounds from the insoluble polymers in the fruits under heat treatments (32, 33).

Amongst the flavanones, the abundance of naringenin and its chalcone increased after stir-frying, whereas that of naringenin O-malonylhexoside decreased. Lou et al. (34) observed increases in the amounts of naringenin, tangeretin, and gallic acid in the immature calamondin after heat treatment at 150°C for 1.5 h and attributed these increases to the liberation of the compound from the immature calamondin peel during heating. Heat treatment could cause the degradation of cell wall structure to release the bound phenolic acids. Accordingly, the contents of tangeretin and gallic acid in the blueberry filling also increased after stir-frying (Table 1). However, none of the flavanones displayed notable changes in fillings during baking. Interestingly, sakuranetin, which is a methylation product of naringenin and has not been detected in the raw fillings, was produced in fried and baked fillings. Although the biosynthesis of sakuranetin from naringenin in plants (35) and the reversable O-demethylation of sakuranetin to naringenin in living tissues under the presence of cytochrome P450 monooxygenases have been well-studied (36), the formation of these components and their interactions and correlations are difficult to speculate without further investigations.

Flavone C-glycosides displayed no variation during baking; during stir-frying, seven flavone C-glycosides showed significant decreases in abundance (Table 1). By contrast, two flavone C-glycosides, namely, 8-C-hexosyl chrysoeriol O-hexoside and hesperetin C-hexoside O-hexoside, increased significantly in the fillings during stir-frying. Moreover, 8-C-hexosyl chrysoeriol O-hexoside was only observed after stir-frying. Therefore, it is a novel degradation product that mostly originated from flavonoids. Amongst the flavones, tricin 5-O-hexosyl-O-hexoside and tricin O-hexosyl-O-syringin alcohol were novel degradation products that appeared after stir-frying. The other flavones either decreased or remained constant during stir-frying.

Variation in Other Chemicals in Blueberry Fillings During Thermal Processing

Other than flavonoids, the other chemicals in the blueberry fillings were also fully investigated by the widely targeted metabolomic analysis with UHPLC-MS/MS to obtain a comprehensive view on the variation in whole chemical compounds in blueberry filling during pastry preparation and gain further insight into the underlying mechanisms.

During thermal processing, the abundance of sucrose decreased by 32% at the stage of stir-frying, while the amount of glucose remained unchanged after stir-frying. The contents of the oxidized products of glucose, gluconic acid and glucono-1,5-lactone, increased by 52 and 123%, respectively (Table 2). However, only the variation in glucono-1,5-lactone was statistically significant. The increases in the amount of gluconic acid and glucono-1,5-lactone might be attributed to the oxidation of glucose induced by H₂O₂ in presence of ferrous ion (37), in which H₂O₂ can be produced from ascorbic acid catalyzed by Cu(II) (38).

Amongst all the 28 amino acids detected, only L-methionine, L-glutamine, and L-(+)-ornithine decreased significantly in the blueberry filling after stir-frying. Baking only affected the amount of L-glutamine. In comparison with the raw filling, baked filling showed higher levels of L-histidine, L-(+)-arginine, L-methionine, and L-glutamine. Pyrrole-2-carboxylic acid, which is a product of the Maillard browning reaction between xylose and amino acid/peptides, increased sharply by 667% after stir-frying (39, 40).

Fang et al. (41) found that triphosphates and diphosphates of nucleotides could be degraded into monophosphate analogs or even nucleosides during heating. This finding might explain the considerable increases in the amounts of nucleosides and their monophosphates in blueberry fillings after heating. The abundance of adenosine, adenosine 3′-monophosphate, inosine 5′-monophosphate, adenosine 5′-monophosphate, guanosine monophosphate, cytidine 5′-monophosphate, and guanosine 3′,5′-cyclic monophosphate increased by 108, 354, 307, 260, 317, 123, and 2206%, respectively, after stir-frying. 2′-Deoxycytidine-5′-monophosphate appeared as a new product in the filling after stir-frying. Surprisingly, the degradation products of methylated nucleosides extensively increased; for instance, 1-methyladenosine increased by 855%, and 5-methyluridine and 1-methylxanthine were produced only after the thermal treatment of stir-frying. The results demonstrated that nucleosides and even free nucleotides would be degraded during heat treatment (41). Moreover, the thermal processing, especially at acidic conditions, might even elicit DNA degradation. Bitskinashvili et al. (42) reported that the combined thermal-acid treatment at 100°C and pH 2–4 could cause considerable degradations of maize and wheat DNA.

Amongst the seven alkaloids, only two showed significant variation after thermal processing, and both displayed increasing trends. However, three tryptamine derivatives were detected to be significantly different and showed reducing trends under the thermal treatment of stir-frying. The other components, such as organic acids, lipids, and amino acid derivatives, showed individually different variation behaviors under thermal processing (Table 2). Amongst the organic acids, only methylglutaric and D-pantothenic acids decreased in terms of abundance upon thermal treatments. Citraconic acid and 2-furoic acid occurred as new products after the thermal treatment of stir-frying. Glutaric acid, which is a degradation product that originated from sugar and lignin under heat treatment, increased by 1.3 times in the filling after stir-frying. D-Erythronolactone, which is a lactone formed by the reaction of D-xylose and calcium hydroxide (43), increased by 186 and 106%, respectively, after stir-frying and baking. Butter was incorporated into the pastry fillings and exposed to high temperature during stir-frying. Heating of butter increased its lactone content (43). As a result, the tridecalactone, which is a component in butter, was observed after stir-frying.

Conclusions

In this study, we applied a widely targeted metabolomic approach to investigate the global chemical changes during different thermal processing procedures, namely, stir-frying and baking, which were involved in the preparation of traditional blueberry filled pastry. A total of 630 chemicals were detected in blueberry fillings, and 288 were screened to be differential chemicals as samples underwent different thermal treatments. In comparison with baking, stir-frying contributed the most to the deviation of the abundance of components in blueberry fillings. Anthocyanin was the most sensitive toward thermal treatments. The other classes of phytochemicals, such as glycerophospholipids, nucleotides and their derivatives, benzoic acid derivatives, flavonols, flavones, flavanones, flavone C-glycosides, hydroxycinnamoyl derivatives, and phenolamides, were also comparably sensitive toward thermal processing. The variation in the chemicals and the novel compounds that appeared after thermal processing detected by this widely targeted metabolomic analysis could provide important insights into the degradation mechanism of phytochemicals and interactions between food ingredients. However, metabolomic analysis for the total understanding of the variation in chemical features and underlying mechanisms during food preparation and corresponding reactions remains a major challenge. The currently developed metabolomic databases are mostly specific to food products, such as FooDB, which is a subset of the Human Metabolite Database (44). This database contains information about thousands of phytochemicals but lacks information on the corresponding chemical degradation products or interaction products. Additional efforts should be applied on the development of metabolomic approaches and their application in investigating food processing.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author Contributions

JZ, WH, SO, and PL contributed conception and design of the study. JZ, ZW, NY, and KZ performed the experiments. JZ, ZW, and PL organized the database and performed the statistical analysis. JZ and SO wrote the first draft of the manuscript. PL and ZW revised the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China [Nos. 31701607 and 31972180]; the Pearl River Talent Program (No. 20170096) and the Science and Technology Planning Project (Grant No. 2018B050502008) by Guangdong Science and Technology Department, Guangdong, China; and the Open Fund of Key Laboratory of Biotechnology and Bioresources Utilization (Dalian Minzu University), Ministry of Education (No. KF2018003), China.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2020.569172/full#supplementary-material

References

1. Grace MH, Xiong J, Esposito D, Ehlenfeldt M, Lila MA. Simultaneous LC-MS quantification of anthocyanins and non-anthocyanin phenolics from blueberries with widely divergent profiles and biological activities. Food Chem. (2019) 227:336–46. doi: 10.1016/j.foodchem.2018.10.101

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Nadulski R, Maslowski A, Mazurek A, Sobczak P, Szmigielski M, Zukiewicz-Sobczak W, et al. Vitamin C and lutein content of northern highbush blueberry (Vaccinium corymbosum L.) juice processed using freezing and thawing. J Food Meas Charact. (2019) 13:2521–8. doi: 10.1007/s11694-019-00172-x

CrossRef Full Text | Google Scholar

3. Nemzer B, Vargas L, Xia XY, Sintara M, Feng H. Phytochemical and physical properties of blueberries, tart cherries, strawberries, and cranberries as affected by different drying methods. Food Chem. (2018) 262:242–50. doi: 10.1016/j.foodchem.2018.04.047

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Cesa S, Carradori S, Bellagamba G, Locatelli M, Casadei MA, Masci A, et al. Evaluation of processing effects on anthocyanin content and colour modifications of blueberry (Vaccinium spp.) extracts: comparison between HPLC-DAD and CIELAB analyses. Food Chem. (2017) 232:114–23. doi: 10.1016/j.foodchem.2017.03.153

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Klopsch R, Baldermann S, Voss A, Rohn S, Schreiner M, Neugart S. Bread enriched with legume microgreens and leaves-ontogenetic and baking-driven changes in the profile of secondary plant metabolites. Front Chem. (2018) 6:322. doi: 10.3389/fchem.2018.00322

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Zhao Y, Wu XL, Yu LL, Chen P. Retention of polyphenols in blueberries (Vaccinium corymbosum) after different cooking methods, using UHPLC–DAD–MS based metabolomics. J Food Compos Anal. (2017) 56:55–66. doi: 10.1016/j.jfca.2016.12.003

CrossRef Full Text | Google Scholar

7. Hager A, Howard LR, Prior RL, Brownmiller C. Processing and storage effects on monomeric anthocyanins, percent polymeric color, and antioxidant capacity of processed black raspberry products. J Food Sci. (2008) 73:H134–H140. doi: 10.1111/j.1750-3841.2008.00855.x

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Yang N, Qiu RX, Yang S, Zhou KN, Wang CT, Ou SY, et al. Influences of stir-frying and baking on flavonoid profile, antioxidant property, and hydroxymethylfurfural formation during preparation of blueberry-filled pastries. Food Chem. (2019) 287:167–75. doi: 10.1016/j.foodchem.2019.02.053

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Teegarden MD, Schwartz SJ, Cooperstone JL. Profiling the impact of thermal processing on black raspberry phytochemicals using untargeted metabolomics. Food Chem. (2019) 274:782–8. doi: 10.1016/j.foodchem.2018.09.053

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Pérez-Míguez R, Sánchez-López E, Plaza M, Castro-Puyana M, Marina ML. A non-targeted metabolomic approach based on reversed-phase liquid chromatography–mass spectrometry to evaluate coffee roasting process. Anal Bioanal Chem. (2018) 410:7859–70. doi: 10.1007/s00216-018-1405-z

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Nikolantonaki M, Julien P, Coelho C, Roullier-Gall C, Ballester J, Schmitt-Kopplin P, et al. Impact of glutathione on wines oxidative stability: a combined sensory and metabolomic study. Front Chem. (2018) 6:182. doi: 10.3389/fchem.2018.00182

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Blasi F, Rocchetti G, Montesano D, Lucini L, Chiodelli G, Ghisoni S, et al. Changes in extra-virgin olive oil added with Lycium barbarum L. carotenoids during frying: chemical analyses and metabolomic approach. Food Res Int. (2018) 105:507–16. doi: 10.1016/j.foodres.2017.11.061

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Gu EJ, Kim DW, Jang GJ, Song SH, Lee JI, Lee SB, et al. Mass-based metabolomic analysis of soybean sprouts during germination. Food Chem. (2017) 217:311–9. doi: 10.1016/j.foodchem.2016.08.113

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Santos MCB, Lima LRD, Nascimento FR, do Nascimento TP, Cameron LC, Ferreira MSL. Metabolomic approach for characterization of phenolic compounds in different wheat genotypes during grain development. Food Res Int. (2019) 124:118–28. doi: 10.1016/j.foodres.2018.08.034

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Yuan YF, Zhao YP, Yang JL, Jiang YM, Lu F, Jia YX, et al. Metabolomic analyses of banana during postharvest senescence by ¹H-high resolution-NMR. Food Chem. (2017) 218:406–12. doi: 10.1016/j.foodchem.2016.09.080

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Buchner N, Krumbein A, Rohn S, Kroh LW. Effect of thermal processing on the flavonols rutin and quercetin. Rapid Commun Mass Sp. (2006) 20:3229–35. doi: 10.1002/rcm.2720

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Hou ZH, Qin PY, Zhang Y, Cui SH, Ren GX. Identification of anthocyanins isolated from black rice (Oryza sativa L.) and their degradation kinetics. Food Res Int. (2013) 50:691–7. doi: 10.1016/j.foodres.2011.07.037

CrossRef Full Text | Google Scholar

18. Cevallos-Casals BA, Cisneros-Zevallos L. Stability of anthocyanin-based aqueous extracts of Andean purple corn and red-fleshed sweet potato compared to synthetic and natural colorants. Food Chem. (2004) 86:69–77. doi: 10.1016/j.foodchem.2003.08.011

CrossRef Full Text | Google Scholar

19. Ortu E, Caboni P. Levels of 5-hydroxymethylfurfural, furfural, 2-furoic acid in sapa syrup, Marsala wine and bakery products. Int J Food Prop. (2017) 20:S2543–51. doi: 10.1080/10942912.2017.1373668

CrossRef Full Text | Google Scholar

20. Zhang YY, Song Y, Hu XS, Liao XJ, Ni YY, Li QH. Effects of sugars in batter formula and baking conditions on 5-hydroxymethylfurfural and furfural formation in sponge cake models. Food Res Int. (2012) 49:439–45. doi: 10.1016/j.foodres.2012.07.012

CrossRef Full Text | Google Scholar

21. Louarme L, Billaud C. Evaluation of ascorbic acid and sugar degradation products during fruit dessert processing under conventional or ohmic heating treatment. Lwt-Food Sci Technol. (2012) 49:184–7. doi: 10.1016/j.lwt.2011.12.035

CrossRef Full Text | Google Scholar

22. Hou DX, Kai K, Li JJ, Lin SG, Terahara N, Wakamatsu M, et al. Anthocyanidins inhibit activator protein 1 activity and cell transformation: structure-activity relationship and molecular mechanisms. Carcinogenesis. (2003) 25:29–36. doi: 10.1093/carcin/bgg184

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Hou DX, Yanagita T, Uto T, Masuzaki S, Fujii M. Anthocyanidins inhibit cyclooxygenase-2 expression in LPS-evoked macrophages: structure-activity relationship and molecular mechanisms involved. Biochem Pharmacol. (2005) 70:417–25. doi: 10.1016/j.bcp.2005.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Tanaka T, Tanaka T, Tanaka M. Potential cancer chemopreventive activity of protocatechuic acid. Int J Clin Exp Med. (2011) 3:27–33. doi: 10.1016/j.jecm.2010.12.005

CrossRef Full Text | Google Scholar

25. Teegarden MD, Knobloch TJ, Weghorst CM, Cooperstone JL, Peterson DG. Storage conditions modulate the metabolomic profile of a black raspberry nectar with minimal impact on bioactivity. Food Funct. (2018) 9:4593–601. doi: 10.1039/C8FO00639C

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Hussain MI, Syed QA, Khattak MNK, Hafez B, Reigosa MJ, El-Keblawy A. Natural product coumarins: biological and pharmacological perspectives. Biologia. (2019) 74:863–88. doi: 10.2478/s11756-019-00242-x

CrossRef Full Text | Google Scholar

27. Zhao MY, Luo YH, Li Y, Liu X, Wu JH, Liao XJ, et al. The identification of degradation products and degradation pathway of malvidin-3-glucoside and malvidin-3,5-diglucoside under microwave treatment. Food Chem. (2013) 141:3260–7. doi: 10.1016/j.foodchem.2013.05.147

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Sadilova E, Stintzing FC, Carle R. Thermal degradation of acylated and non-acylated anthocyanins. J Food Sci. (2006) 71:C504–12. doi: 10.1111/j.1750-3841.2006.00148.x

CrossRef Full Text | Google Scholar

29. Mancilla-Margalli NA, López MG. Generation of maillard compounds from inulin during the thermal processing of Agave tequilana Weber var. Azul J Agric Food Chem. (2002) 50:806–12. doi: 10.1021/jf0110295

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Lourenço A, Gominho J, Pereira H. Chemical characterization of lignocellulosic materials by analytical pyrolysis. In: P Kusch, editor. Analytical Pyrolysis. London, UK: IntechOpen Limited. (2018). p. 1–22. doi: 10.5772/intechopen.80556

CrossRef Full Text | Google Scholar

31. Wong Z, Chen KF, Li J. Formation of vanillin and syringaldehyde in an oxygen delignification process. Bioresources. (2010) 5:1509–16.

Google Scholar

32. Choi MY, Chai C, Park JH, Lim J, Lee J, Kwon SW. Effects of storage period and heat treatment on phenolic compound composition in dried Citrus peels (Chenpi) and discrimination of Chenpi with different storage periods through targeted metabolomic study using HPLC-DAD analysis. J Pharmaceut Biomed. (2011) 54:638–45. doi: 10.1016/j.jpba.2010.09.036

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Jeong SM, Kim SY, Kim DR, Jo SC, Nam KC, et al. Effect of heat treatment on the antioxidant activity of extracts from citrus peels. J Agric Food Chem. (2004) 52:3389–93. doi: 10.1021/jf049899k

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Lou SN, Lin YS, Hsu YS, Chiu EM, Ho CT. Soluble and insoluble phenolic compounds and antioxidant activity of immature calamondin affected by solvents and heat treatment. Food Chem. (2014) 161:246–53. doi: 10.1016/j.foodchem.2014.04.009

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Rakwal R, Agrawal GK, Yonekura M, Kodama O. Naringenin 7-O-methyltransferase involved in the biosynthesis of the flavanone phytoalexin sakuranetin from rice (Oryza sativa L.). Plant Sci. (2000) 155:213–21. doi: 10.1016/S0168-9452(00)00223-5

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Katsumata S, Hamana K, Horie K, Toshima H, Hasegawa M. Identification of sternbin and naringenin as detoxified metabolites from the rice flavanone phytoalexin sakuranetin by Pyricularia oryzae. Chem Biodivers. (2017) 14:e1600240. doi: 10.1002/cbdv.201600240

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Rinsant D, Chatel G, Jérôme F. Efficient and selective oxidation of D-glucose into gluconic acid under low-frequency ultrasonic irradiation. Chemcatchem. (2014) 6:3355–9. doi: 10.1002/cctc.201402604

CrossRef Full Text | Google Scholar

38. Zhou P, Zhang J, Zhang Y, Liu Y, Liang J, Liu B, et al. Generation of hydrogen peroxide and hydroxyl radical resulting from oxygen-dependent oxidation of L-ascorbic acid via copper redox-catalyzed reactions. RSC Adv. (2016) 6:38541–7. doi: 10.1039/C6RA02843H

CrossRef Full Text | Google Scholar

39. Yen GC, Lee TC. Mutagen formation in the reaction of maillard browning products, 2-acetylpyrrole and its analogues, with nitrite. Food Chem Toxicol. (1986) 24:1303–8. doi: 10.1016/0278-6915(86)90062-1

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Zou TT, Kang L, Yang C, Song HL, Liu Y. Flavour precursor peptide from an enzymatic beef hydrolysate Maillard reaction-II: mechanism of the synthesis of flavour compounds from a sulphur-containing peptide through a Maillard reaction. Lwt-Food Sci Technol. (2019) 110:8–18. doi: 10.1016/j.lwt.2019.04.022

CrossRef Full Text | Google Scholar

41. Fang ML, Ivanisevic J, Benton HP, Johnson CH, Patti GJ, Hoang LT, et al. Thermal degradation of small molecules: a global metabolomic investigation. Anal Chem. (2015) 87:10935–41. doi: 10.1021/acs.analchem.5b03003

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Bitskinashvili K, Gabriadze I, Kutateladze T, Vishnepolsky B, Mikeladze D, Datukishvili N. Effects of thermal-acid treatment on degradation and amplification of wheat and maize DNA. J Food Nutr Res. (2018) 57:242–51.

Google Scholar

43. Maga JA. Lactones in foods. Crit Rev Food Sci. (1976) 8:1–56. doi: 10.1080/10408397609527216

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Wishart DS, Jewison T, Guo AC, Wilson M, Knox C, Liu YF, et al. HMDB 3.0-the human metabolome database in 2013. Nucl Acids Res. (2013) 41:801–7. doi: 10.1093/nar/gks1065

CrossRef Full Text | Google Scholar