Next Article in Journal
Alpha- and Beta-Cyclodextrin Inclusion Complexes with 5-Fluorouracil: Characterization and Cytotoxic Activity Evaluation
Next Article in Special Issue
The Traditional Chinese Medicine DangguiBuxue Tang Sensitizes Colorectal Cancer Cells to Chemoradiotherapy
Previous Article in Journal
A Simple Precursor for Highly Functionalized Fused Imidazo[4,5-b]pyridines and Imidazo[4,5-b]-1,8-naphthyridine
Previous Article in Special Issue
Organ- and Growing Stage-Specific Expression of Solanesol Biosynthesis Genes in Nicotiana tabacum Reveals Their Association with Solanesol Content
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pathway Analysis and Metabolites Identification by Metabolomics of Etiolation Substrate from Fresh-Cut Chinese Water Chestnut (Eleocharis tuberosa)

College of Food, Hainan University, Haikou 570228, China
*
Author to whom correspondence should be addressed.
Submission received: 9 October 2016 / Revised: 20 November 2016 / Accepted: 21 November 2016 / Published: 1 December 2016

Abstract

:
Fresh-cut Chinese water chestnuts (CWC) turn yellow after being peeled, reducing their shelf life and commercial value. Metabolomics, the systematic study of the full complement of small molecular metabolites, was useful for clarifying the mechanism of fresh-cut CWC etiolation and developing methods to inhibit yellowing. In this study, metabolic alterations associated with etiolation at different growth stages (0 day, 2 days, 3 days, 4 days, 5 days) from fresh-cut CWC were investigated using LC–MS and analyzed by pattern recognition methods (principal component analysis (PCA), partial least squares-discriminant analysis (PLS-DA), and orthogonal projection to latent structures-discriminant analysis (OPLS-DA)). The metabolic pathways of the etiolation molecules were elucidated. The main metabolic pathway appears to be the conversion of phenylalanine to p-coumaroyl-CoA, followed by conversion to naringenin chalcone, to naringenin, and naringenin then following different pathways. Firstly, it can transform into apigenin and its derivatives; secondly, it can produce eriodictyol and its derivatives; and thirdly it can produce dihydrokaempferol, quercetin, and myricetin. The eriodictyol can be further transformed to luteolin, cyanidin, dihydroquercetin, dihydrotricetin, and others. This is the first reported use of metabolomics to study the metabolic pathways of the etiolation of fresh-cut CWC.

1. Introduction

Metabolomics, a new field in systems biology, represents an emerging and powerful discipline concerned with the comprehensive analysis of small molecules (<1 kDa); it provides a powerful approach for the discovery of biomarkers in biological systems [1,2,3,4]. Metabolomics provides a quantitative analysis of all low-molecular-mass metabolites and lends insight into the relationship between metabolites and changes in physiology/pathophysiology status in biochemical networks and pathways by combining a range of different analytical technologies and calculation methods [5,6,7]. These techniques have recently demonstrated significant potential in many fields, such as toxicology [8], nutrition [9], environmental stress [10], global effects of genetic manipulation, disease diagnosis, drug mechanism and development, natural product discovery, and comparison of different growth stages [11,12]. Additionally, metabolomics has been used in food science to assess food quality, food safety, food microbiology, and food processing [13,14,15]. Furthermore, the global metabolite profile now enables much more detailed studies of noninvasive biomarkers that can detect different growth stages in plants.
Plant studies using metabolomics techniques aim for the simultaneous detection of all metabolites in plant tissues. Metabolomics detection technology includes nuclear magnetic resonance (NMR) spectroscopy, liquid chromatography-mass spectrometry (LC–MS) [16], and gas chromatography–mass spectrometry (GC–MS). Among these, GC–MS is mainly suited for compound classes appearing mainly in primary metabolic pathways, such as amino acids, fatty acids, carbohydrates, and organic acids [17]. LC–MS is more suitable for assessing the overall biochemical richness of plants. This technique analytically covers the large (semipolar) group of plant secondary metabolites including alkaloids, saponins, phenolic acids, phenylpropanoids, flavonoids, glucosinolates, polyamines, and their derivatives [18,19]. Therefore, LC–MS is preferred over the other techniques in the study of numerous structure databases [20]. At present, LC–MS has become one of the frequently used techniques in metabolomics studies because of its high sensitivity and reproducibility.
The Chinese water chestnut (CWC, Eleocharis tuberosa) is one of the most popular hydrophytic vegetables in China because of its unique taste [21]. However, it is usually peeled before being eaten, and fresh-cut CWC turns yellow after being peeled, reducing its shelf life and its commercial value. Recent studies have shown that CWC etiolation is linked to the presence of eriodictyol and naringenin [21]. In addition, most researches have shown that phenylalanine ammonia-lyase (PAL) was present in fresh-cut CWC tissue. With the extension of storage time, the activity of PAL increases and the etiolation of CWC also increases. Therefore, the metabolic pathways of eriodictyol and naringenin are believed to be associated with PAL activity [22]. PAL is a critical and rate-limiting enzyme in the phenylpropanoid metabolic pathway, a pathway that also contains eriodictyol and naringenin [23,24,25,26]. It is therefore probable that the etiolation process of fresh-cut CWC is related to the functioning of the phenylpropanoid metabolic pathway. However, there may be a variety of specific metabolic pathways related to yellowing, and few prior studies have investigated these pathways fully.
The synthesis of flavonoids begins with the shikimic acid metabolic pathway, which generates phenylalanine. This molecule then enters the phenylpropanoid metabolic pathway, and then, through the common enzyme chalcone synthase, enters the flavonoid synthesis pathway, which generates the various flavonoids through different branches [27,28]. Though the use of metabolomics to research flavonoids is relatively uncommon, both naringenin and eriodictyol are considered flavonoids, and the molecular weight of most flavonoids is below 1 kDa; therefore, the metabolites of the etiolation process are suitable for studying with metabolomics techniques. Literature about the process of involving metabolomics in the research of flavonoids is very common, especially in herbal aspects [29,30]. In terms of fruit and vegetables, the glycosides of luteolin and apigenin as well as quercetin conjugates were identified from Phoenix dactylifera (date palm) fruits using UPLC/PDA/ESI-TOF-MS methods [31,32]. Flavonoids in blueberry extracts were analyzed through HPLC-ESI-MS/MS, which identified that anthocyanins and flavonoid alcohols are the main components [32]. Therefore, there is precedence for identifying the characteristic flavonoids and metabolic pathways from plants.
In order to clarify the mechanism of fresh-cut CWC etiolation and yellowing, as well as to evaluate the safety of the etiolation metabolite, we sought to discover the metabolic pathways of the molecules responsible for etiolation. In this study, we used the LC–MS method combined with additional analytical technologies and calculation methods (principal component analysis (PCA), partial least squares-discriminant analysis (PLS-DA), orthogonal projection to latent structures-discriminant analysis (OPLS-DA)) to analyze the metabolic pathways of etiolation in fresh-cut CWC.

2. Results and Discussion

2.1. LC–MS Analysis of Metabolic Patterns

In this study, five time points from the yellowing process of CWC were studied using an LC–MS metabolomics approach. For the analysis, we compared 0 day (see Methods) with 2 days (0–2), 2 days with 3 days (2–3), 3 days with 4 days (3–4), and 4 days with 5 days (4–5) for screening and identification of the metabolites. According to peak height data, we detected 1046 positive-ion (ESI+) peaks and 1962 negative-ion (ESI−) peaks by Mass Profiler using the same acquiring method of 0.5–15 min of retention time, and then the normalized data were fed to SIMCA-P. The parameters (R2X, R2Y, and Q2Y) obtained from SIMCA-P indicated that the model was able to distinguish between each of the yellowing time points in both ion modes.
To further evaluate the capability of the LC–MS-based metabolomics approach to distinguish between all of the time points, multivariate data analyses (PCA, PLS-DA, OPLS-DA) were performed. In the PCA analysis, we identified five principal components in positive mode (R2X = 0.464, Q2 = 0.124) (Figure 1). For this nonsupervisory PCA model, the R2X parameter-which describes the rate of the model-can be used to determine the model’s quality. Generally, a model is reliable when R2X > 0.4 [33]. Thus, the current PCA model was suited to explain the metabolic differences between each group. Every point in Figure 2 represents a sample with the fresh-cut CWC; from the PCA analysis it can be visually seen that there is a significant difference between each pair of time points, except for between 3 days and 4 days. Furthermore, the between-group differences were significantly greater than the intergroup differences, indicating that substantial biochemical perturbation occurred in all groups. PCA as an unsupervised analysis method, responsive to the primitive state of the data, can be used to observe the relationship between each group, but it cannot account for intergroup errors or eliminate the random error unrelated to this research. In order to obtain information on which metabolites were the cause of the differences, we conducted PLS-DA and OPLS-DA analyses for each pair of time points.
The supervised PLS-DA analysis revealed greater differences in the metabolic signature between each pair of time points, focusing on intergroup differences and minimizing between-group differences to better grasp the overall characteristics and variation of the data. For each pair of time points studied, this analysis was performed in both the positive mode and negative mode. Results for each pair are as follows (Figure 2): in positive mode, groups 0–2 obtained one principal component (R2X = 0.253, R2Y = 0.97, Q2 = 0.845); groups 2–3, one principal component was obtained (R2X = 0.328, R2Y = 0.923, Q2 = 0.68); groups 3–4 obtained one principal component (R2X = 0.179, R2Y = 0.897, Q2 = 0.588); groups 4–5, one principal component was obtained (R2X = 0.198, R2Y = 0.867, Q2 = 0.514); in negative mode, groups 0–2 obtained one principal component (R2X = 0.278, R2Y = 0.945, Q2 = 0.829); 2–3, one principal component was obtained (R2X = 0.298, R2Y = 0.913, Q2 = 0.787); groups 3–4 obtained two principal components (R2X = 0.321, R2Y = 0.991, Q2 = 0.806); groups 4–5, one principal component was obtained (R2X = 0.259, R2Y = 0.879, Q2 = 0.712). The model parameters (R2X and Q2) indicate that the current PLS-DA models were reliable and suitable for interpretation between each pair of time points of metabolic differences and finding the differences. Meanwhile, the left of the arrange experiments were randomly generated depending on whether the value of R2 and Q2 was less than the original value from the right. This indicates that the predictive power of the original model is greater than any randomized y variables, thus, the model is effective and it can continue to be applied to the different metabolites.
Typically, OPLS-DA analysis is conducted to search for the different components that lead to the maximum differences between two samples. In order to further obtain the metabolic differences, this study applied an OPLS-DA approach; separating samples into two blocks. As can be seen from Figure 3, compared with PCA analysis, each group received the highest degree of separation, which is useful to determine the metabolic differences.

2.2. Identification of Metabolites

All metabolites from each total ion chromatogram were extracted and aligned using Mass Profiler software. All information, including retention time, peak intensity, and exact mass, was supplied by the LC-Q/TOF-MS platform. The metabolites (VIP (variable importance in the projection) > 1) were analyzed by t-test, and metabolites with a p-value < 0.05 were identified as metabolic differences.

2.2.1. Identification of Metabolites in 0–2 Days

Comparing the samples of 0 day with 2 days, the main metabolite differences are shown in Table 1 and Table 2. In the positive mode, 5′-deoxy-5′-(methylthio) adenosine, ferulic acid, and linolenic acid showed significant increases, while N-palmitoyl threonine and stearidonic acid decreased. Among these, the linolenic acid at 2 days compared with 0 day increased by about 19 times. In the negative mode, salicylaldehyde, salicylic acid, naringenin, eriodictyol, pentahydroxyflavanone, apigenin-glucuronide, and luteolin 5-methyl ether 7-glucoside were increased significantly. The naringenin has increased by about 4820 times and the eriodictyol was increased by about 700 times, while the PE (phosphatidyl ethanolamine) (18:3), PE (18:1), and quercetin 3-O-(6″-malonyl-glucoside) 7-O-glucoside at 2 days compared with 0 day were decreased by 600 times. Thus, the metabolites between 0 day and 2 days have significant differences, and these metabolites were related to the pathways of fresh-cut CWC yellowing.

2.2.2. Identification of Metabolites in 2–3 Days

Comparing the samples of 2 days with 3 days, the main metabolite differences are shown in Table 3 and Table 4. In the positive mode, cyanidin 3-O-(6-O-malonyl-β-d-glucoside) and dihydroxy-6,7,4′,5′-tetramethoxyflavone increased significantly, while aconitic acid and MG (18:2) decreased. In the negative mode, phenylacetic acid, vanillic acid, quinic acid, quercetin and its derivatives, dihydroxy-dimethoxy-methylenedioxy, the derivatives of naringenin, malvidin, and luteolin were increased significantly, all of them showing an average increase of about 280 times. In contrast, ascorbic acid and arginine at 3 days, compared with 2 days, were decreased by 4 times. The major differences were identified to further explore the pathway of CWC yellowing.

2.2.3. Identification of Metabolites in 3–4 Days

Comparing the samples of 3 days with 4 days, the main metabolite differences are shown in Table 5 and Table 6. Overall, the metabolic differences between 3 days and 4 days were few. In the positive mode, dihydroxy-7,2′-dimethoxyflavanone at 4 days compared with 3 days was increased by about 3.5 times. On the other hand, trihydroxy-3-methoxyflavone decreased by about 12 times. From the negative mode, kolaflavanone, eriodictyol, the derivative of naringenin, cyanidin, and myricetin were increased. In contrast, linolenic acid and phenylacetic acid at 4 days, compared with 3 days were decreased.

2.2.4. Identification of Metabolites in 4–5 Days

Comparing the samples of 4 days with 5 days, the main metabolite differences are shown in Table 7 and Table 8. In the positive mode, dehydrophytosphingosine at 5 days, compared with 4 days, increased by about 5400 times. From the negative mode, naringenin-4’-O-β-d-glucuronide changes significantly. At 5 days, compared with the 4 days, it fell by about 110 times; other metabolic changes were small.

2.3. Biological Pathway and Function Analysis

Using pattern recognition analysis on the metabolites, we were able to achieve a clear separation between each pair of time points. To identify the pathways involved in fresh-cut CWC etiolation, previously identified metabolites were analyzed using MetPA (Metabolomic Pathway Analysis) software. MetPA is a web-based tool that derives much of its power from the KEGG metabolic pathways database. In metabolic networks, changes in more “important” positions result in a more substantial effects on the pathway compared to changes occurring in marginal or relatively isolated positions. In our analysis, we identified several pathways that appear to be involved in etiolation, including the glycolysis pathway, fatty acid metabolic pathway, citrate cycle pathway, and flavonoid biosynthesis pathway. The green parts in Figure 4 represent the differences we have detected. In the glycolysis pathway, the α-d-glucose is converted to α-d-glucose-6p and then further transformed into glyceraldehyde-3P, while the α-d-glucose-6p is converted to histidine; glyceraldehyde-3P and glycerate-1,3P2 are convert to each other, and at the same time converted to 2-C-methyl-d-4-erythritol-4P. Glycerate-3P and phosphoenol-pyruvate are convert to each other and generate phosphoserine. Pyruvate is created from phosphoenol-pyruvate, and converts to octanoate (8:0), triggered the citrate cycle pathway. Furthermore, pyruvic, leucine, and l-alanine interactions in the flavonoid biosynthesis pathway have a certain effect.
The major substance that causes the fresh-cut CWC yellowing are naringenin and eriodictyol, both found in the flavonoid biosynthesis pathway. Therefore, the flavonoid biosynthesis pathway is the most important pathway for CWC yellowing. As we can see from Figure 5, phenylalanine is converted to p-coumaroyl-CoA, followed by conversion to naringenin chalcone, to naringenin, and then naringenin enters different pathways. Firstly, it can transform into apigenin and its derivatives; secondly, it can produce eriodictyol and its derivatives; and thirdly it can produce dihydrokaempferol, quercetin, and myricetin. From here, eriodictyol can be further transformed to luteolin, cyanidin, dihydroquercetin, dihydrotricetin, and others. According to the metabolism pathway, quercetin, myricetin, apigenin, luteolin, and other substances were also yellowing substances. Therefore, the yellowing of fresh-cut CWC is caused by several yellow substances, and our study indicated that the flavonoid biosynthetic pathway is the main mechanism for fresh-cut CWC yellowing. Detailed construction of the yellowing metabolism pathway is shown in Figure 5.

3. Materials and Methods

3.1. Chemicals, Reagents and Materials

Acetonitrile (LC/MS grade) and methanol (HPLC grade) were purchased from Merck (Damstadt, Gemany). Formic acid was purchased from CNW (Düsseldorf, Gemany). Leucine-enkephalin was obtained from Sigma Company (San Francisco, CA, USA). Distilled water was purchased from Watsons (Hong Kong, China). NaClO were purchased from “Xiya” Company (Chengdu, China).
Fresh CWCs were purchased in Haikou, China and stored in the laboratory at 17 °C. Fruits without physical damage or diseases and with uniform size were selected for analysis. CWCs were treated as reported by You [34] and Peng [35]. Fresh CWCs (100 kg) were washed, peeled using a sharp, stainless steel knife, and chopped into small, thick slices. The slices were then disinfected with 0.1% NaClO for 15 min before airing. The CWCs were packed with plastic packaging machines and preserved at 17 °C for 5 days.
To assess the primary metabolites as the CWCs turn yellow, CWCs were randomly selected from five different yellowing stages (Figure 5): one hour post-cutting (0 day), two days post-cutting (2 days), three days post-cutting (3 days), four days post-cutting (4 days), and five days post-cutting (5 days). Slices taken from the yellowed surface were then used as samples. Since the sample between 0 day and 1 d showed little difference, the sample of 1 d was removed from further analyses. Six samples from each of the five time stages (0.5 g each) were transferred to a 1.5 mL centrifuge tube, mixed with 1 mL solvent (MeOH/H2O, 7:3 v/v) by smashing for 1 min, extracted by ultrasonication for 10 min, and then separated by centrifugation at 12,000 rpm at 4 °C for 15 min. Finally, 100 μL aliquots of the supernatant were transferred to autosampler vials.

3.2. Metabolic Profiling and Data Analysis

Chromatographic analysis was performed on an Agilent C18 column (Agilent, 100 mm × 2.1 mm, 1.8 μm) using an Agilent LC system (Agilent, 1290 Infinity LC, 6530 UHD, Palo Alto, CA, USA). The column temperature was maintained at 40 °C and eluted at a flow rate of 0.4 mL/min. The mobile phase consisted of a linear gradient system of (A) 0.1% formic acid in water and (B) 0.1% formic acid in acetonitrile. The gradient elution process: 0–3 min, 20%–50% B; 3–13 min, 50%–95% B; 13–15 min, 95% B. The injection volume was 4 μL and the autosampler temperature was 4 °C.
The eluent was introduced into the mass spectrometer (Agilent 6530 Accurate-Mass Q-TOF/MS). The MS system was operated using positive-ion (ESI+) and negative-ion (ESI−) modes. Conditions of positive-ion analysis were as follows: nitrogen was used as the nebulizing gas; the cone gas flow was set at 50 L/h, and the desolvation gas flow was set at 600 L/h. The capillary voltage was 4.0 kV, the sampling cone voltage was 35.0 kV, the source temperature was 100 °C, the desolvation temperature was 350 °C, and the extraction cone voltage was 4.0 V. In negative-ion mode, the cone gas flow was set at 50 L/h and the desolvation gas flow was set at 700 L/h. The capillary voltage was 3.5 kV, the sampling cone voltage was 50 kV, the source temperature was 100 °C, the desolvation temperature was 300 °C, and the extraction cone voltage was 4.0 V. The scan time was set to 0.03 s with a 0.02 s interscan time. Data were collected from 50 to 1000 m/z. For accurate and repeatable mass acquisition, a lock-mass of leucine-enkephalin was used via a lock spray interface monitor for positive-ion mode ([M + H]+ = 556.2771 Da) and negative-ion mode ([M − H] = 554.2615 Da) to ensure accuracy during the MS analysis.
All data were processed using Mass Profiler in the MassHunter software platform (Agilent) for LC/MS data preprocessing, followed by manual postediting in Microsoft Excel. We removed impurity peaks due to column bleeding and sample preparation. The intensity of each ion was normalized with respect to the total ion count to generate a data matrix that consisted of the rt_mass (retention time_mass), the quantity of observe, and the normalized peak. After that, the multivariate data matrix was introduced into SIMCA-P 13.0 software (Umetrics, Umeå, Sweden) for PCA. PCA data were visualized by plotting the PC scores, with each point representing one rt_mass pair. Thus, the loading plot gives an indication of the metabolites that most strongly influence the patterns in the score plot. In order to further discriminate the metabolites and biomarkers, the data were analyzed using PLS-DA and OPLS-DA methods. Variables that significantly discriminated between groups were considered potential biomarkers and were assigned with p-values using a student’s test (t-test); p-values less than 0.05 were considered statistically significant [36]. The identities of specific metabolites were confirmed by comparison of their mass spectra and mass-to-charge ratio (m/z) with free online databases, including METLIN (http://metlin.scripps.edu/).
The construction, interaction, and pathway analyses of potential biomarkers were performed with MetaboAnalyst 2.0, with database sources, including METLIN and KEGG (http://www.genome.jp/kegg/), to help identify pathways that were most significantly altered.

4. Discussion

The yellowing of fresh-cut CWC is one of the leading causes of reduced shelf life and commercial value, and it is known to be enzymatically different from the browning of other vegetables and fruits [37,38,39]. Despite attempts to reduce CWC yellowing [37,40], we still lack methods to fundamentally inhibit etiolation, and the mechanism of yellowing remains unclear. The advent of the metabolomics field represents a paradigm shift in metabolic research, away from methods that focus on a limited number of single pathways and towards approaches that capture the entire complexity of metabolic networks [41]. Emerging techniques in metabolomics provide a powerful platform for discovering novel biomarkers and biochemical pathways [42]. This may allow us to develop techniques that inhibit CWC yellowing, potentially leading to future advances in inhibition of these pathways.
In this study, samples from different time points of CWC yellowing were analyzed by LC-Q/TOF-MS, coupled with multivariate data analyses; the major metabolic differences were obtained. All samples were successfully modeled, and we identified significant differences in PCA scores between each pair of time points. The pattern was verified as authentic by the PLS-DA analysis. Then, the difference between each pair of time points was analyzed by OPLS-DA combined with t-test, and the metabolites were obtained. It was found that metabolic differences include organic acids (such as salicylic acid, citric acid, ferulic acid, leucine, l-aspartic acid, lysine etc.) and flavonoids (like eriodictyol, naringenin, quercetin, myricetin, apigenin, luteolin etc.), while sucrose mainly existed in negative mode. Among those mentioned, two of the flavonoids were yellowing substances. The content changes of eriodictyol and naringenin is outstanding, which showed an upward trend that was proportional to the degree of fresh-cut CWC yellowing during the storage time. The change of myricetin showed a moderate upward trend during storage. In comparison with eriodictyol and naringenin, the change in apigenin and luteolin showed less of an upward trend and did not rise after that time, and the content of quercetin showed fluctuations in growth, which is currently under further investigation. It was demonstrated that the yellowing phenomenon of fresh-cut CWC was caused by a series of yellowing substances, mainly eriodictyol and naringenin, which is consistent with the study of Pan [21]; the content of flavonoids that cause yellowing was proportional to the degree of fresh-cut CWC yellowing. Nowadays, flavonoids are known to be key antioxidant substances that have antiaging, disease resistance, antimicrobial, anti-inflammatory, and antitumor activities, and are becoming more and more abundant in modern life. This study aims to lay the foundation for the further study of fresh-cut CWC yellowing by investigating the substances causing the phenomenon [43]. It was proved that yellowing substances of fresh-cut CWC were connected with the glycolysis pathway, fatty acid metabolism pathway, citrate cycle pathway, and flavonoid biosynthetic pathway by identifying metabolic differences combined with METLIN and KEGG data. Among the metabolites identified, organic acids began to accumulate in the initial stage of the fresh-cut CWC yellowing and were consumed as respiratory matrixes of the glycolysis pathway and citrate cycle pathway in their last stages. Organic acids also participate in respiration and metabolic processes of phenols, amino acids, esters, and flavonoids [44], explaining that organic acids increased gradually at the early time points and descend gradually at the end.
The flavonoids causing the yellowing of fresh-cut CWC mainly derived from the flavonoid biosynthetic pathway. It was proved that phenylalanine was converted into p-coumaric-CoA in the first step, which consistent with the study of Burbulis [45]. Naringenin, which is a significant checkpoint in the flavonoid biosynthetic pathway, is converted into the second key substance, eriodictyol, which was converted into dihydrokaempferol, then quercetin, and finally, myricetin. In the metabolic pathway, luteolin, apigenin and its derivatives, quercetin, and myricetin were also yellowing substances, in addition to naringenin and eriodictyol, which revealed that the yellowing is the result of the joint action of flavonoids, mainly naringenin and eriodictyol. This suggests that the etiolation molecules are a mixture of metabolites, but it is clear that the main substances are naringenin and eriodictyol. According to other research reports, in the tissue of grape and corn, 4-coumarin-CoA is generated by phenylalanine in the phenylpropanoid metabolic pathway, which is converted into naringenin chalcone, then naringenin, eriodictyol, and pentahydroxyflavanone were generated in the end [46]; these reports are consistent with this study. The pathway of fresh-cut CWC yellowing is complex, and each pathway contains multiple metabolites; there are rare reports about the relationship between the intermediate products and their interaction mechanism that leads to yellow that it is a complex pathway. In the further work, an in vitro model will be set up to simulate the process of flesh-cut CWC yellowing to confirm the reliability of the metabolic pathways.

5. Conclusions

To the best of our knowledge, this is the first report of a metabolomics study of the etiolation of fresh-cut CWC. A rapid, reliable, and sensitive LC-Q/TOF-MS method was validated for qualitative analysis of the different stages of CWC yellowing. Metabolites from each time point could be distinguished from each other by this method, which demonstrates that a significant metabolic perturbation is occurring. These 3008 peaks (1046 peaks in positive mode and 1962 in negative mode) could be useful as biomarkers for the different stages of yellowing. Interestingly, the flavanoid metabolic pathway emerged as a key pathway in the etiolation of fresh-cut CWC. This study also identified eriodictyol, naringenin, and sucrose as important features in this process, verifying the results of our previous study [21]. This investigation is a good starting point for providing a better theoretical understanding of how to efficiently inhibit the yellowing of CWC. However, further studies are needed for a deeper understanding of the metabolic pathways and enzymes that are associated with CWC yellowing.

Acknowledgments

This project was supported by the National Natural Science Foundation of China (Grant No. 31360414). We sincerely thank Shanghai Sensichip Infotech Co. Ltd. for their technical support.

Author Contributions

Author Yi-Xiao Li and Yong-Gui Pan designed the research and discussed the data; Yi-Xiao Li, Feng-Ping He, Shang-Bin Li and Meng-Qi Yuan performed the experiments; Yong-Gui Pan provided the sample; Yi-Xiao Li and Yong-Gui Pan analyzed the data; Yi-Xiao Li wrote the paper; Yi-Xiao Li and Yong-Gui Pan modified the language of the paper. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Armitage, E.G.; Barbas, C. Metabolomics in cancer biomarker discovery: Current trends and future perspectives. J. Pharm. Biomed. Anal. 2014, 87, 1–11. [Google Scholar] [CrossRef] [PubMed]
  2. Li, P.L.; Sun, H.G.; Hua, Y.L.; Ji, P.; Zhang, L.; Li, J.X.; Wei, Y. Metabolomics study of hematopoietic function of Angelica sinensis on blood deficiency mice model. J. Ethnopharmacol. 2015, 166, 261–269. [Google Scholar] [CrossRef] [PubMed]
  3. Zhao, Y.Y.; Zhang, L.; Long, F.Y.; Cheng, X.L.; Bai, X.; Wei, F.; Lin, R.C. UPLC-Q-TOF/HSMS/MSE-based metabonomics for adenine-induced changes in metabolic profiles of rat faeces and intervention effects of ergosta-4,6,8(14),22-tetraen-3-one. Chem. Biol. Interact. 2013, 201, 31–38. [Google Scholar] [CrossRef] [PubMed]
  4. Gao, X.; Guo, M.; Zhao, B.; Peng, L.; Su, J.; Bai, X.; Li, J.; Qiao, Y. A urinary metabonomics study on biochemical changes in yeast-induced pyrexia rats: A new approach to elucidating the biochemical basis of the febrile response. Chem. Biol. Interact. 2013, 204, 39–48. [Google Scholar] [CrossRef] [PubMed]
  5. Zhang, A.; Sun, H.; Han, Y.; Yuan, Y.; Wang, P.; Song, G.; Yuan, X.; Zhang, M.; Xie, N.; Wang, X. Exploratory urinary metabolic biomarkers and pathways using UPLC-Q-TOF-HDMS coupled with pattern recognition approach. Analyst 2012, 137, 4200–4208. [Google Scholar] [CrossRef] [PubMed]
  6. Holmes, E.; Loo, R.L.; Stamler, J.; Bictash, M.; Yap, I.K.; Chan, Q.; Ebbels, T.; De Iorio, M.; Brown, I.J.; Veselkov, K.A. Human metabolic phenotype diversity and its association with diet and blood pressure. Nature 2008, 453, 396–400. [Google Scholar] [CrossRef] [PubMed]
  7. Wang, F.; Liu, X.; Liu, C.; Liu, Z.; Sun, L. Effects of antibiotic antitumor drugs on nucleotide levels in cultured tumor cells: an exploratory method to distinguish the mechanisms of antitumor drug action based on targeted metabolomics. Acta Pharm. Sin. B 2015, 5, 223–230. [Google Scholar] [CrossRef] [PubMed]
  8. Mc Auley, M.T.; Choi, H.; Mooney, K.; Paul, E.; Miller, V.M. Systems biology and synthetic biology: A new epoch for toxicology research. Adv. Toxicol. 2015, 2015, 575403. [Google Scholar] [CrossRef]
  9. Freemark, M. Metabolomics in nutrition research: Biomarkers predicting mortality in children with severe acute malnutrition. Food Nutr. Bull. 2015, 36 (Suppl. S1), 88S–92S. [Google Scholar] [CrossRef]
  10. Masclaux-Daubresse, C.; Clément, G.; Anne, P.; Routaboul, J.-M.; Guiboileau, A.; Soulay, F.; Shirasu, K.; Yoshimoto, K. Stitching together the multiple dimensions of autophagy using metabolomics and transcriptomics reveals impacts on metabolism, development, and plant responses to the environment in Arabidopsis. Plant Cell Online 2014, 26, 1857–1877. [Google Scholar] [CrossRef] [PubMed]
  11. Wang, X.; Sun, H.; Zhang, A.; Sun, W.; Wang, P.; Wang, Z. Potential role of metabolomics apporoaches in the area of traditional Chinese medicine: As pillars of the bridge between Chinese and Western medicine. J. Pharm. Biomed. Anal. 2011, 55, 859–868. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, A.; Sun, H.; Yan, G.; Wang, P.; Wang, X. Mass spectrometry-based metabolomics: Applications to biomarker and metabolic pathway research. Biomed. Chromatogr. 2016, 30, 7–12. [Google Scholar] [CrossRef] [PubMed]
  13. Chin, E.; Slupsky, C.; Weimer, B.; Slupsky, C. Applications of metabolomics in food science: Food composition and quality, sensory and nutritional attributes. In Metabolomics in Food and Nutrition; Woodhead Publ. Ltd.: Cambridge, UK, 2013; pp. 217–230. [Google Scholar]
  14. McGhie, T.K.; Rowan, D.D. Metabolomics for measuring phytochemicals, and assessing human and animal responses to phytochemicals, in food science. Mol. Nutr. Food Res. 2012, 56, 147–158. [Google Scholar] [CrossRef] [PubMed]
  15. Cevallos-Cevallos, J.M.; Reyes-De-Corcuera, J. Metabolomics in food science. Adv. Food Nutr. Res. 2011, 67, 1–24. [Google Scholar]
  16. Cai, Y.; Gao, Y.; Tan, G.; Wu, S.; Dong, X.; Lou, Z.; Zhu, Z.; Chai, Y. Myocardial lipidomics profiling delineate the toxicity of traditional Chinese medicine Aconiti Lateralis radix praeparata. J. Ethnopharmacol. 2013, 147, 349–356. [Google Scholar] [CrossRef] [PubMed]
  17. Gullberg, J.; Jonsson, P.; Nordström, A.; Sjöström, M.; Moritz, T. Design of experiments: an efficient strategy to identify factors influencing extraction and derivatization of Arabidopsis thaliana samples in metabolomic studies with gas chromatography/mass spectrometry. Anal. Biochem. 2004, 331, 283–295. [Google Scholar] [CrossRef] [PubMed]
  18. De Vos, R.C.; Moco, S.; Lommen, A.; Keurentjes, J.J.; Bino, R.J.; Hall, R.D. Untargeted large-scale plant metabolomics using liquid chromatography coupled to mass spectrometry. Nat. Protoc. 2007, 2, 778–791. [Google Scholar] [CrossRef] [PubMed]
  19. T’Kindt, R.; Morreel, K.; Deforce, D.; Boerjan, W.; van Bocxlaer, J. Joint GC-MS and LC-MS platforms for comprehensive plant metabolomics: Repeatability and sample pre-treatment. J. Chromatogr. B 2009, 877, 3572–3580. [Google Scholar] [CrossRef] [PubMed]
  20. Yin, P.; Xu, G. Current state-of-the-art of nontargeted metabolomics based on liquid chromatography-mass spectrometry with special emphasis in clinical applications. J. Chromatogr. A 2014, 1374, 1–13. [Google Scholar] [CrossRef] [PubMed]
  21. Pan, Y.-G.; Li, Y.-X.; Yuan, M.-Q. Isolation, purification and identification of etiolation substrate from fresh-cut Chinese water-chestnut (Eleocharis tuberosa). Food Chem. 2015, 186, 119–122. [Google Scholar] [CrossRef] [PubMed]
  22. Pan, Y.; Chen, W. Study on the enzyme activity correlated to browning in the fresh-cut Chinese water-chestnut. Anhui Agric. Sci. 2007, 35, 9856–9857. [Google Scholar]
  23. Reuben, S.; Rai, A.; Pillai, B.V.; Rodrigues, A.; Swarup, S. A bacterial quercetin oxidoreductase QuoA-mediated perturbation in the phenylpropanoid metabolic network increases lignification with a concomitant decrease in phenolamides in Arabidopsis. J. Exp. Bot. 2013, 64, 5183–5194. [Google Scholar] [CrossRef] [PubMed]
  24. Yin, R.; Messner, B.; Faus-Kessler, T.; Hoffmann, T.; Schwab, W.; Hajirezaei, M.-R.; von Saint Paul, V.; Heller, W.; Schäffner, A.R. Feedback inhibition of the general phenylpropanoid and flavonol biosynthetic pathways upon a compromised flavonol-3-O-glycosylation. J. Exp. Bot. 2012, 63, 2465–2478. [Google Scholar] [CrossRef] [PubMed]
  25. Lin, Y.; Jain, R.; Yan, Y. Microbial production of antioxidant food ingredients via metabolic engineering. Curr. Opin. Biotechnol. 2014, 26, 71–78. [Google Scholar] [CrossRef] [PubMed]
  26. Yuan, Y.; Song, L.; Li, M.; Liu, G.; Chu, Y.; Ma, L.; Zhou, Y.; Wang, X.; Gao, W.; Qin, S. Genetic variation and metabolic pathway intricacy govern the active compound content and quality of the Chinese medicinal plant Lonicera japonica thunb. BMC Genom. 2012, 13, 195. [Google Scholar] [CrossRef] [PubMed]
  27. Samanta, A.; Das, G.; Das, S.K. Roles of flavonoids in plants. Int. J. Pharm. Sci. Technol. 2011, 6, 12–35. [Google Scholar]
  28. Ibrahim Alzand, K.; Abdalkarim Mohamed, M. Flavonoids: Chemistry, Biochemistry and Antioxidant activity. J. Pharm. Res. 2012, 5, 4013. [Google Scholar]
  29. Song, H.-H.; Ryu, H.W.; Lee, K.J.; Jeong, I.Y.; Kim, D.S.; Oh, S.-R. Metabolomics investigation of flavonoid synthesis in soybean leaves depending on the growth stage. Metabolomics 2014, 10, 833–841. [Google Scholar] [CrossRef]
  30. Carvalho, E.; Franceschi, P.; Feller, A.; Palmieri, L.; Wehrens, R.; Martens, S. A targeted metabolomics approach to understand differences in flavonoid biosynthesis in red and yellow raspberries. Plant Physiol. Biochem. 2013, 72, 79–86. [Google Scholar] [CrossRef] [PubMed]
  31. Farag, M.A.; Mohsen, M.; Heinke, R.; Wessjohann, L.A. Metabolomic fingerprints of 21 date palm fruit varieties from Egypt using UPLC/PDA/ESI-qTOF-MS and GC–MS analyzed by chemometrics. Food Res. Int. 2014, 64, 218–226. [Google Scholar] [CrossRef]
  32. Yuan, S.; Yao, S.; Geng, Y.; Cai, P.; Zhang, W. Identification of Anthocyanins and Flavonols in Extract of Blueberry by Using HPLC-ESI-MS/MS. Acta Chem. Scand. 2009, 67, 318–322. [Google Scholar]
  33. Li, J.; Yang, L.; Li, Y.; Tian, Y.; Li, S.; Jiang, S.; Wang, Y.; Li, X. Metabolomics study on model rats of chronic obstructive pulmonary disease treated with Bu-Fei Jian-Pi. Mol. Med. Rep. 2015, 11, 1324–1333. [Google Scholar] [CrossRef] [PubMed]
  34. You, Y.; Jiang, Y.; Sun, J.; Liu, H.; Song, L.; Duan, X. Effects of short-term anoxia treatment on browning of fresh-cut Chinese water chestnut in relation to antioxidant activity. Food Chem. 2012, 132, 1191–1196. [Google Scholar] [CrossRef]
  35. Pen, L.; Jiang, Y. Effects of chitosan coating on shelf life and quality of fresh-cut Chinese water chestnut. LWT Food Sci. Technol. 2003, 36, 359–364. [Google Scholar] [CrossRef]
  36. Zaykin, D.V.; Zhivotovsky, L.A.; Westfall, P.H.; Weir, B.S. Truncated product method for combining p-values. Genet. Epidemiol. 2002, 22, 170–185. [Google Scholar] [CrossRef] [PubMed]
  37. Zvitov-Ya’ari, R.; Nussinovitch, A. Browning prevention in rehydrated freeze-dried non-blanched potato slices by electrical treatment. LWT Food Sci. Technol. 2014, 56, 194–199. [Google Scholar] [CrossRef]
  38. Gomes, M.H.; Vieira, T.; Fundo, J.F.; Almeida, D.P. Polyphenoloxidase activity and browning in fresh-cut ‘Rocha’pear as affected by pH, phenolic substrates, and antibrowning additives. Postharvest Biol. Technol. 2014, 91, 32–38. [Google Scholar] [CrossRef]
  39. Sagar, V.; Kumar, P.S. Recent advances in drying and dehydration of fruits and vegetables: A review. J. Food Sci. Technol. 2010, 47, 15–26. [Google Scholar] [CrossRef] [PubMed]
  40. Zhou, D.; Li, L.; Wu, Y.; Fan, J.; Ouyang, J. Salicylic acid inhibits enzymatic browning of fresh-cut Chinese chestnut (Castanea mollissima) by competitively inhibiting polyphenol oxidase. Food Chem. 2015, 171, 19–25. [Google Scholar] [CrossRef] [PubMed]
  41. Pah, A.; Guimera, R.; Mustoe, A.; Amaral, L. Use of a global metabolic network to curate organismal metabolic networks. Sci. Rep. 2013, 3, 1695. [Google Scholar] [CrossRef] [PubMed]
  42. Bai, Y.; Zhang, H.; Sun, X.; Sun, C.; Ren, L. Biomarker identification and pathway analysis by serum metabolomics of childhood acute lymphoblastic leukemia. Clin. Chim. Acta 2014, 436, 207–216. [Google Scholar] [CrossRef] [PubMed]
  43. Pereira, D.M.; Valentão, P.; Pereira, J.A.; Andrade, P.B. Phenolics: From chemistry to biology. Molecules 2009, 14, 2202–2211. [Google Scholar] [CrossRef]
  44. Zhang, J. Metabolic Profiling of the Major Fruit Quality Components during the Development of Fragaria ananassa “Fengxiang” Fruits; Anhui Agriculture University: Hefei, China, 2009. [Google Scholar]
  45. Burbulis, I.E.; Winkel-Shirley, B. Interactions among enzymes of the Arabidopsis flavonoid biosynthetic pathway. Proc. Natl. Acad. Sci. USA 1999, 96, 12929–12934. [Google Scholar] [CrossRef] [PubMed]
  46. Ferreyra, M.F.; Rius, S.P.; Casati, P. Flavonoids: Biosynthesis, biological functions, and biotechnological applications. Front. Plant Sci. 2012, 3, 1–15. [Google Scholar]
  • Sample Availability: Not Available.
Figure 1. Statistical principal component analysis (PCA) analysis of the 5 time points.
Figure 1. Statistical principal component analysis (PCA) analysis of the 5 time points.
Molecules 21 01648 g001
Figure 2. Effective verification of fresh-cut Chinese water chestnut during different storage times. (A) 0–2 comparison, ESI+; (B) 2–3 comparison, ESI+; (C) 3–4 comparison, ESI+; (D) 4–5 comparison, ESI+; (E) 0–2 comparison ESI−; (F) 2–3 comparison, ESI−; (G) 3–4 comparison, ESI−; (H) 4–5 comparison, ESI−.
Figure 2. Effective verification of fresh-cut Chinese water chestnut during different storage times. (A) 0–2 comparison, ESI+; (B) 2–3 comparison, ESI+; (C) 3–4 comparison, ESI+; (D) 4–5 comparison, ESI+; (E) 0–2 comparison ESI−; (F) 2–3 comparison, ESI−; (G) 3–4 comparison, ESI−; (H) 4–5 comparison, ESI−.
Molecules 21 01648 g002
Figure 3. Orthogonal projection to latent structures-discriminant analysis (OPLS-DA) results of all samples. (A) 0–2 comparison, scanned by ESI+; (B) 2–3 comparison, ESI+; (C) 3–4 comparison, ESI+; (D) 4–5 comparison, ESI+. “*” means the sample separation significantly.
Figure 3. Orthogonal projection to latent structures-discriminant analysis (OPLS-DA) results of all samples. (A) 0–2 comparison, scanned by ESI+; (B) 2–3 comparison, ESI+; (C) 3–4 comparison, ESI+; (D) 4–5 comparison, ESI+. “*” means the sample separation significantly.
Molecules 21 01648 g003
Figure 4. Metabolomic Pathway Analysis (MetPA) construction of the metabolic pathways involved in the etiolation of fresh-cut CWC.
Figure 4. Metabolomic Pathway Analysis (MetPA) construction of the metabolic pathways involved in the etiolation of fresh-cut CWC.
Molecules 21 01648 g004
Figure 5. The yellowing process from the fresh-cut Chinese water chestnut. Each day has 5 replicate samples.
Figure 5. The yellowing process from the fresh-cut Chinese water chestnut. Each day has 5 replicate samples.
Molecules 21 01648 g005
Table 1. The main identification obtained between 0 day and 2 days in ESI+ (p < 0.05).
Table 1. The main identification obtained between 0 day and 2 days in ESI+ (p < 0.05).
VIPRT (min)MassNameT-TestFold (2/0) *
1.24813.14278.2252Linolenic Acid0.001+6.567
1.4084.91194.0582ferulic acid0.000+4.237
1.7794.06297.095′-Deoxy-5′-(methylthio)adenosine0.001+4.040
1.4796.60548.1175Chrysoeriol 7-O-(6″-malonyl-glucoside)0.000+2.683
1.20410.19357.2886N-palmitoyl threonine0.009−2.457
1.8438.06276.2095Stearidonic Acid0.006−2.086
1.1374.33137.0477Anthranilic acid0.000−1.716
1.1814.39264.137Abscisic Acid0.038−1.408
1.6680.82115.064Proline0.007+0.275
* 0 day and 2 days is the logarithm of the ratio of the value of value (base 2), “+” means 2 days increase relative to 0 day, “−” means decrease. VIP: Variable Importance in the Projection; RT: retention time.
Table 2. The main identification obtained between 0 day and 2 days in ESI− (p < 0.05).
Table 2. The main identification obtained between 0 day and 2 days in ESI− (p < 0.05).
VIPRT (min)MassNameT-TestFold (2/0) *
1.6127.57272.0694Naringenin0.001+28.844
1.7785.82288.0642Eriodictyol0.000+26.056
1.7736.11462.1179Luteolin 5-methyl ether 7-glucoside0.000+24.638
1.32511.24479.3027PE (18:1)0.017−24.315
1.75913.13278.2251Linolenic Acid0.00024.309
1.1533.47712.1513Quercetin 3-O-(6″-malonyl-glucoside) 7-O-glucoside0.046−24.296
1.6244.87304.0593Pentahydroxyflavanone0.001+24.082
1.4594.97450.1182Eriodictyol 5-O-glucoside0.006+23.416
1.39810.03475.271PE (18:3)0.010−23.115
1.4826.57446.0862Apigenin-glucuronide0.005+22.622
1.45910.68477.2872PE (18:2)0.006−5.758
1.4781.31180.064Glucose0.005−3.408
1.6985.34448.1027Naringenin-4′-O-β-d-Glucuronide0.000+3.061
1.3395.21670.1768Limocitrin 3,7-diglucoside0.015+2.829
1.4804.48176.06872-Isopropylmaleate0.005+2.556
1.2451.11668.1634Euphorbianin0.027+2.481
1.6344.91180.0438Caffeic Acid0.001+2.355
1.2355.79490.1128Cyanidin 3-(6-acetylgalactoside)0.029+2.216
1.7164.58192.0629Quinic acid0.000+1.921
1.5155.58207.0886N-Acetyl-l-phenylalanine0.003−1.786
1.5874.95354.0962Chlorogenic Acid0.0021.717
1.6088.12492.1279Malvidin 3-O-glucoside0.001−0.310
1.3695.69168.042Vanillic acid0.012−0.247
1.3696.30272.0687Naringeninchalcone0.012−0.247
1.3699.73356.1271Laurifolin0.012−0.247
1.3695.78534.1009Cyanidin 3-O-(6-O-malonyl-β-d-glucoside)0.012−0.247
1.3695.23535.1341Malvidin 3-(6-acetylglucoside)0.012−0.247
1.3696.31558.1173Dihydromorelloflavone0.012−0.247
1.3696.97754.1943Luteolin 7-glucuronide-3′,4′-dirhamnoside0.012−0.247
1.3698.24784.1795Apigenin 7-Glucosyl-(1->2)-glucuronide-4′-glucuronide0.012−0.247
* 0 day and 2 days is the logarithm of the ratio of the value of value (base 2), “+” means 2 days increase relative to 0 day, “−” means decrease. VIP: Variable Importance in the Projection; RT: retention time.
Table 3. The main identification obtained between 2 days and 3 days in ESI+ (p < 0.05).
Table 3. The main identification obtained between 2 days and 3 days in ESI+ (p < 0.05).
VIPRT (min)MassNameT-TestFold (3/2) *
1.48210.60354.2777MG (18:2)0.023−3.108
1.3336.11534.1018Cyanidin 3-O-(6-O-malonyl-β-d-glucoside)0.010+2.596
1.0811.04174.0166Aconitic acid0.007−2.281
1.0375.95374.0975Dihydroxy-6,7,4′,5′-tetramethoxyflavone0.022+2.078
1.6566.34316.0954Dihydroxy-dimethoxyflavanone0.018+0.879
1.4490.82115.064Proline0.019−0.282
* 2 days and 3 days is the logarithm of the ratio of the value of value (base 2), “+” means 3 days increases relative to 2 days, “−” means decrease. VIP: Variable Importance in the Projection; RT: retention time.
Table 4. The main identification obtained between 2 days and 3 days in ESI− (p < 0.05).
Table 4. The main identification obtained between 2 days and 3 days in ESI− (p < 0.05).
VIPRT (min)MassNameT-TestFold (3/2) *
1.2894.578192.0634Quinic acid0.019+25.997
1.5583.970191.0619N-Acetyl-dl-methionine0.002+25.741
1.8106.100448.1022Naringenin-4′-O-β-d-Glucuronide0.000+24.727
1.8305.969332.0543Quercetagetin-methyl ether0.000+24.325
1.2959.730356.1271Laurifolin0.018+23.249
1.7395.783534.1009Luteolin 5-(6″-malonylglucoside)0.000+22.968
1.5885.234535.1341Malvidin 3-(6-acetylglucoside)0.001+22.258
1.2645.144358.068Dihydroxy-dimethoxy-methylenedioxy0.022+21.350
1.2007.953610.0885Gallocatechin 7,4′-di-O-gallate0.032+21.183
1.36110.026475.271PE (18:3)0.011+21.183
1.7795.245318.0388Myricetin0.000+6.475
1.25411.404532.1541Flavonol 3-O-d-xylosylglucoside0.023+5.586
1.4225.464534.0963Cyanidin 3-O-(6-O-malonyl-β-d-glucoside)0.007+4.072
1.7375.877218.1155Pantothenic Acid0.000+3.618
1.7577.131302.0438Quercetin0.000+3.545
1.3238.530286.0838Naringenin 5-methyl ether0.015+3.037
1.6064.247304.0591Pentahydroxyflavanone0.001+2.656
1.3426.523476.1319Naringenin 7-O-beta-d-glucoside 6″-acetate0.013+2.508
1.3995.708450.1176Eriodictyol-O-glucoside0.008+2.388
1.5771.056176.0325Ascorbic acid0.001−2.228
1.39513.135278.2251Linolenic Acid0.009+1.268
1.1235.787490.1128Cyanidin 3-(6-acetylgalactoside)0.049+1.194
1.3384.912180.0438Caffeic Acid0.013−1.091
1.3120.721422.0835Sucrose-6-phosphate0.016−0.608
* 2 days and 3 days is the logarithm of the ratio of the value of value (base 2), “+” means 3 days increases relative to 2 days, “−” means decrease. VIP: Variable Importance in the Projection; RT: retention time.
Table 5. The main identification obtained between 3 days and 4 days in ESI+ (p < 0.05).
Table 5. The main identification obtained between 3 days and 4 days in ESI+ (p < 0.05).
VIPRT (min)MassNameT-TestFold (4/3) *
1.3837.64300.064Trihydroxy-3-methoxyflavone0.010−3.618
1.0396.34316.0954Dihydroxy-7,2′-dimethoxyflavanone0.000+1.796
1.15510.43354.2779MG (18:2)0.037−1.697
1.89911.12356.2934MG (18:1)0.036−1.169
1.3334.02117.0579Indole0.010+0.508
2.04110.39198.1623Lauroleic acid0.029+0.304
* 3 days and 4 days is the logarithm of the ratio of the value of value (base 2), “+” means 4 days increases relative to 3 days, “−” means decrease. VIP: Variable Importance in the Projection; RT: retention time.
Table 6. The main identification obtained between 3 days and 4 days in ESI− (p < 0.05).
Table 6. The main identification obtained between 3 days and 4 days in ESI− (p < 0.05).
VIPRT (min)MassNameT-TestFold (4/3) *
2.1208.67588.1278Kolaflavanone0.000+5.933
2.1049.73356.1271Laurifolin0.000+3.667
1.5161.09433.1176Cyanidin 3-rhamnoside0.034+3.486
1.5791.32216.0402-C-Methyl-d-erythritol 4-phosphate0.026−3.477
2.2154.87138.0319Salicylic acid0.000+2.569
1.8715.50302.0797Ferreirin0.004+2.537
1.5180.68250.153Xanthoxin0.034−2.407
2.1737.52316.0956Eriodictyol 7,3′-dimethyl ether0.000+2.239
2.2415.69168.042Vanillic acid0.000+2.192
1.6647.95610.0885Epigallocatechin 5,7-di-O-gallate0.017+2.162
2.0285.14358.068Dihydroxy-dimethoxy-methylenedioxy0.001+1.843
2.0963.97191.0619N-Acetyl-dl-methionine0.001+1.830
1.5544.97450.1182Eriodictyol 6-O-glucoside0.029+1.777
1.8785.48448.1021Naringenin 5-O-glucuronide0.004+1.723
1.7657.94316.0593Quercetin 3-methyl ether0.009+1.717
1.5655.82288.0642Eriodictyol0.027+1.613
1.6838.53286.0838Naringenin 5-methyl ether0.015+1.450
1.6216.60136.0522Phenylacetic acid0.021−1.387
2.0325.25318.0388Myricetin0.001+1.256
1.7844.58192.0634Quinic acid0.008+1.190
1.64013.13278.2251Linolenic Acid0.019−1.021
* 3 days and 4 days is the logarithm of the ratio of the value of value (base 2), “+” means 4 days increases relative to 3 days, “−” means decrease. VIP: Variable Importance in the Projection; RT: retention time.
Table 7. The main identification obtained between 4 days and 5 days in ESI+ (p < 0.05).
Table 7. The main identification obtained between 4 days and 5 days in ESI+ (p < 0.05).
VIPRT (min)MassNameT-TestFold (5/4) *
1.22211.44315.278Dehydrophytosphingosine0.023+22.358
1.54713.14278.2252Linolenic Acid0.000+2.291
1.10610.43354.2779MG(18:2)0.017+1.874
1.3376.98708.189Patuletin 3,7-(3-acetylrhamnoside)0.004+1.359
1.1107.64300.064Trihydroxy-3-methoxyflavone0.032+1.167
1.3170.66155.0697Histidine0.003-0.467
* 4 days and 5 days is the logarithm of the ratio of the value of value (base 2), “+” means 5 days increases relative to 4 days, “−” means decrease. VIP: Variable Importance in the Projection; RT: retention time.
Table 8. The main identification obtained between 4 days and 5 days in ESI− (p < 0.05).
Table 8. The main identification obtained between 4 days and 5 days in ESI− (p < 0.05).
VIPRT (min)MassNameT-TestFold (5/4) *
1.4426.18448.1026Naringenin-4′-O-β-d-Glucuronide0.015−23.434
1.2591.06176.0325Ascorbic acid0.042−3.930
1.8284.39130.063Ketoleucine0.000+2.189
1.84013.13278.2251Linolenic Acid0.000+2.116
1.6407.07264.1366Abscisic Acid0.003+2.097
1.9545.23535.1341Malvidin 3-(6-acetylglucoside)0.000+2.016
1.5035.21670.1768Limocitrin 3,7-diglucoside0.010+1.913
1.5154.78122.0368Salicylaldehyde0.009−1.547
1.43613.73280.2405Linoleic acid0.0161.395
1.6941.34131.0952Leucine0.002−1.374
1.5807.90302.0799Ferreirin0.005+1.275
1.2383.96354.0963Chlorogenic Acid0.046−1.241
1.2286.74270.0534Trihydroxyflavone0.048−1.213
1.8435.69168.042Vanillic acid0.000+1.072
1.8070.69155.0699Histidine0.000−0.979
1.9512.22148.0526Cinnamic acid0.000−0.827
1.2809.73356.1271Laurifolin0.038+0.821
1.4764.58376.0789Tetrahydroxy-3,6,2′-trimethoxyflavone0.012−0.815
1.2598.67588.1278Kolaflavanone0.042+0.763
1.3525.89194.0581ferulic acid0.026−0.639
1.5434.01117.0581Indole0.007−0.609
1.4930.74147.0545Glutamate0.011−0.503
1.6104.58422.0842Sucrose-6-phosphate0.004−0.513
* 4 days and 5 days is the logarithm of the ratio of the value of value (base 2), “+” means 5 days increases relative to 4 days, “−” means decrease. VIP: Variable Importance in the Projection; RT: retention time.

Share and Cite

MDPI and ACS Style

Li, Y.-X.; Pan, Y.-G.; He, F.-P.; Yuan, M.-Q.; Li, S.-B. Pathway Analysis and Metabolites Identification by Metabolomics of Etiolation Substrate from Fresh-Cut Chinese Water Chestnut (Eleocharis tuberosa). Molecules 2016, 21, 1648. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules21121648

AMA Style

Li Y-X, Pan Y-G, He F-P, Yuan M-Q, Li S-B. Pathway Analysis and Metabolites Identification by Metabolomics of Etiolation Substrate from Fresh-Cut Chinese Water Chestnut (Eleocharis tuberosa). Molecules. 2016; 21(12):1648. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules21121648

Chicago/Turabian Style

Li, Yi-Xiao, Yong-Gui Pan, Feng-Ping He, Meng-Qi Yuan, and Shang-Bin Li. 2016. "Pathway Analysis and Metabolites Identification by Metabolomics of Etiolation Substrate from Fresh-Cut Chinese Water Chestnut (Eleocharis tuberosa)" Molecules 21, no. 12: 1648. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules21121648

Article Metrics

Back to TopTop