Next Article in Journal
Free and Bound Phenolic Profiles and Antioxidant Activities in Melon (Cucumis melo L.) Pulp: Comparative Study on Six Widely Consumed Varieties Planted in Hainan Province
Previous Article in Journal
Quality Characteristics of Raspberry Fruits from Dormancy Plants and Their Feasibility as Food Ingredients
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Foods
Volume 12
Issue 24
10.3390/foods12244445
foods-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolomics Study of the Effect of Transcription Factor NOR-like1 on Flavonoids in Tomato at Different Stages of Maturity Using UPLC-MS/MS

by
Di Guan
1,
Ying Zhao
1,
Xiaodan Zhao
1,* and
Daqi Fu
2
1
School of Food and Health, Beijing Technology and Business University, Beijing 100048, China
2
College of Food Science and Nutritional Engineering, China Agricultural University, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Foods 2023, 12(24), 4445; https://doi.org/10.3390/foods12244445
Submission received: 12 November 2023 / Revised: 4 December 2023 / Accepted: 8 December 2023 / Published: 11 December 2023
(This article belongs to the Topic Edible Plants as Sources of Polyphenols: Cultivation Techniques, Extraction, and Nutraceutical Applications)
Browse Figures
Versions Notes

Abstract

:
Tomato fruits are rich in flavonoids. This study explores the effect of transcription factor SlNOR-like1 on the accumulation of flavonoids in tomato fruits at different ripening stages. We used ultra-pressure liquid chromatography–tandem mass spectrometry (UPLC-MS/MS) to analyze wild-type (WT) and NOR-like1 CRISPR/Cas9-edited (NOR-like1) tomato fruits. A total of 50 flavonoid metabolites were accurately identified and determined in tomatoes. The flavonoid metabolic differences were observed among the different tomato sample groups using PCA and OPLS-DA analysis. There were 16 differential flavonoids (13 upregulated and 3 downregulated) identified between WT-GR (WT tomato at the green-ripening stage) and NOR-like1-GR (NOR-like1 tomato at the green-ripening stage), 9 differential flavonoids (six upregulated and three downregulated) identified between WT-BR3 (WT tomato at the color-breaking stage) and NOR-like1-BR3 (NOR-like1 tomato at the color-breaking stage), and 12 differential flavonoids (11 upregulated and 1 downregulated) identified between WT-BR9 (WT tomato at the red-ripening stage) and NOR-like1-BR9 (NOR-like1 tomato at the red-ripening stage). Rutin, nicotiflorin, naringenin chalcone, eriodictyol, and naringenin-7-glucoside were the five flavonoids with the highest content in the ripening stages (BR3 and BR9) in both WT and NOR-like1 tomato fruits. The overall flavonoid contents in WT tomato fruits changed little from GR to BR3 and decreased from BR3 to BR9; meanwhile, in the NOR-like1 tomato fruits, the total amounts of the flavonoids exhibited an increasing trend during all three ripening stages. The accumulation pattern of flavonoid metabolites in NOR-like1 tomato fruits differed from that in WT tomato fruits, especially in the later ripening process of BR9. The transcription factor SlNOR-like1 has an impact on the accumulation of flavonoids in tomato fruits. The results provide a preliminary basis for subsequent research into its regulatory mechanism and will be helpful for attaining future improvements in the nutritional quality and postharvest treatment of tomato fruits.

1. Introduction

Tomatoes contain many flavonoids. As typical and representative polyphenolic compounds, flavonoids have very important functional effects. Flavonoids can effectively improve the stress resistance of plants and resist a variety of biological and abiotic stresses [1,2] occurring in the process of plant growth and development; they also have antioxidant and other effects in animal experiments and human disease studies [3]. Tomato is an important vegetable crop throughout the world, and the study of its methods of genome repair is of great significance for the stable quality of tomato [4]. The ripening of tomatoes leads to significant physiological and biochemical changes, and primary metabolites play crucial roles during fruit development and ripening, improving tomato taste and flavor [5].
Flavonoids are secondary metabolites that are modified by a wide range of chemical reactions and are ubiquitous in the biosynthetic pathways in plants [6]. The synthesis of flavonoid compounds, including chalcones, flavonoids, flavonols, anthocyanins, and proanthocyanidins, is a complex process, and their metabolic pathways have been extensively studied using biochemical and molecular biology techniques (Figure 1) [7,8,9].
In recent years, many scholars have conducted extensive research on flavonoid substances in fruits, including their composition, content, detection methods, and related genes [10,11,12,13]. At present, the influence of transcription factors on flavonoids has mainly been focused on the MYB, bHLH, and WD40 transcription factors; research into NAC transcription factors is only in the initial stage, and the function and regulation of NAC genes is relatively poorly studied in tomatoes [14]. Previous studies have confirmed that NAC has different expression patterns in different species and have isolated and identified NAC transcription factors in plants to study their effects on plant growth and development, biotic and abiotic stress, and plant metabolites.
The study of transcription factor families is an important area of post-genomic research. The biosynthetic pathways of plant flavonoids comprise actions of structural genes and regulatory genes, and the enzymes in the synthetic pathway are catalyzed by the proteins encoding synthetic structural genes, including FLS, F3H, CHS, CHI, C4H, etc. Moreover, they are regulated by transcription factors [15].
Huang et al. (2022) identified 183 NAC transcription factors from the radish OZ-16 genome [16]. The expression patterns of 10 RsNAC genes were verified by means of qPCR, which provided comprehensive information on the radish NAC gene family and identified candidate RsNACs associated with taproot discoloration. Sun et al. (2019) discovered an NAC transcription factor, MdNAC52, which accelerates the expression of this gene when the color of apples deepens, and the overexpression of the MdNAC52 gene can stimulate the accumulation of anthocyanins [17]. Duan et al. (2020) identified a total of 111 putative NAC transcription factors based on celery transcriptomic and genomic databases. Under field conditions, the overexpression of SNAC1 significantly improved the drought tolerance of transgenic rice without affecting yield [18,19]. Rice overexpressing OsNAC10 resulted in higher yields under both normal and field drought conditions [20]. Guo et al. (2022) identified IbNAC genes associated with abiotic stress through RNA-seq detection and qRT-PCR analysis, and bioinformatics data supported the prospect of genetically engineering candidate regulators, such as IbNAC006, IbNAC029, IbNAC138, and IbNAC143, to improve the stress tolerance of sweet potatoes [21]. Soybean GmNAC06 transcription factor and capsicum CaNAC46 transcription factor were found to be positive regulators of salt stress tolerance [22,23]. Similarly, SINAC4 and SINAC1I genes are involved in the regulation of salt tolerance and early tolerance in tomatoes [24,25].
NOR-like1 is a novel positive regulator of tomato fruit ripening and plays an important role in transcriptional regulatory networks. CRISPR/Cas9 gene editing technology was used to obtain NOR-like1 gene knockout mutants, and it was found that the fruit ripening of NOR-like1 mutant fruits was delayed by 14 days compared with wild-type tomatoes; moreover, the ripening process after color breaking was also significantly inhibited [26]. Compared with the wild-type, fruit softening was significantly inhibited in CR-NOR (CRISPR/Cas9-edited NAC-NOR) but was slightly affected in the NOR mutant and accelerated in overexpressed NAC-NOR (oe-NOR) fruit [27].
In our previous study, we preliminarily explored the effects of the NOR-like1 transcription factor on the profile of metabolites at different stages of tomato maturity using transcriptomics and untargeted metabolomics and found that NOR-like1 significantly affected the changes in flavonoids at different developmental stages [28]. Therefore, the specific influence of the transcription factor SlNORlike1 on tomato flavonoid compounds needs to be further studied. In this study, fruits of wild-type and CRISPR/Cas9-edited lines were collected at different ripening stages: green ripening (GR), color breaking (BR3), and red ripening (BR9). We quantitatively analyzed the flavonoids in the samples. Regression curves of flavonoid standards were established. The flavonoids identified in the sample were fitted to the standard curve, and a total of 50 flavonoids were detected. Then, we analyzed these 50 flavonoids using multivariate analysis. From a paired comparison, we screened out the differential flavonoids between WT and NOR-like1 tomato fruits at each ripening stage and explored the trends in the accumulated flavonoids in WT and NOR-like1 during ripening to understand the effect of NOR-like1 on flavonoids during tomato ripening.

2. Materials and Methods

2.1. Plant Materials and Sample Preparation

Wild-type tomato Ailsa Craig (AC) and a NOR-like1 tomato transgenic line employing CRISPR/Cas9 gene-editing techniques were both grown in a greenhouse at China Agricultural University. Fruits of wild-type and NOR-like1 fruits were collected at different ripening stages (GR, BR3, and BR9). The sampling method was as described in the previous study [28]. Thus, samples were divided into six groups: WT-GR, WT-BR3, WT-BR9, NOR-like1-GR, NOR-like1-BR3, and NOR-like1-BR9. The samples were ground into powder after freeze-drying. Then, 20 mg of powder was weighed and extracted with 0.5 mL of 70% methanol. [2H6]-Daidzein, [2H3]-rutin, and [13C3]-(+/−)-gallocatechin were used as internal standards, and 10 μL of internal standards (4000 nmol/L) was added into the extract for quantitation. After ultrasonic extraction for 30 min, the extract was centrifuged at 12,000× g at 4 °C for 5 min and the supernatant was filtered through a 0.22 μm membrane for further liquid chromatography–tandem mass spectrometry (LC-MS/MS) analysis.

2.2. Ultra-Performance Liquid Chromatography (UPLC) Conditions

The sample extracts were analyzed using a UPLC-ESI-MS/MS system, ExionLC™ AD (SCIEX, Framingham, MA, USA) coupled with a QTRAP® 6500+ (SCIEX, Framingham, MA, USA) mass spectrometer, which was based on the description by Zhang et al. [29]. Using [2H6] Daidzein, [2H3] rutin, [13C3]-(+/−)-gallocatechin as internal standard substances, the mixed concentration of these internal standard substances is 4000 nmol/L, with an added scalar of 10 μL and an injection volume of 2 μL. All samples were analyzed using Waters ACQUITY UPLC HSS T3 C18 (100 mm × 2.1 mm i.d., 1.8 µm) at 40 °C. The solvent system was water with 0.05% formic acid as the A phase and acetonitrile with 0.05% formic acid as the B phase. At the beginning of the reaction, gradient elution was set as A/B at 90:10 (V/V); then, A/B was adjusted to 80:20 (V/V) after 1 min, 30:70 (V/V) after 9 min, 5:95 (V/V) after 12.5 min, and 90:10 (V/V) after 13.6 min until the end of the reaction. The gradient elution program protocol was performed as described by Wang et al. [30]. The eluent was then connected to an ESI-triple quadrupole-linear ion trap (QTRAP)-MS.

2.3. Electrospray Ionization (ESI)–Tandem Mass Spectrometry (MS/MS) (ESI-MS/MS) Conditions

Linear ion trap (LIT) and triple quadrupole (QQQ) scans were conducted on a triple quadrupole-linear ion trap mass spectrometer (QTRAP) equipped with an ESI operating in positive and negative ion mode. The ESI source operation parameters were as follows: an ion source, ESI+/−; the temperature of the electric spray ion source (ESI) was 550 °C; the mass spectrum voltage was 5500 V in the positive ion mode and −4500 V in the negative ion mode; and the curtain gas pressure (CUR) was 35 psi. The Metware Database (MWDB), self-built based on standards, was applied to qualitatively analyze the data from mass spectrometry detection. Quantitative analysis was completed using the Multiple Reaction Monitoring (MRM) mode of triple-quadrupole mass spectrometry. The mass spectrometry data were processed using Analyst 1.6.3 software (SCIEX, Toronto, ON, Canada). MultiQuant 3.0.3 software (SCIEX, Toronto, ON, Canada) was used to quantify all metabolites.

2.4. Multivariate Statistical Analysis of Flavonoid Metabolites

Unsupervised principal component analysis (PCA) was performed using the built-in statistical prcomp function of the R software of V3.5.1 version (www.r-project.org/; accessed on 10 February 2022); the data were unit variance scaled before unsupervised PCA was conducted.
The metabolite content data in the cluster analysis were standardized using unit variance scaling (UV), and all samples were subjected to cluster heat map analysis. The R program script was used to draw the cluster heat map.
OPLS-DA was used to screen differential variables by removing uncorrelated differences. Firstly, log2-centralization was performed on the data; then, the MetaboAnalystR package (version V1.0.1) in the R language (version V3.5.1) was used for OPLSR (the R package script used for volcanic maps is ggplot2 v3.36), the Anal function was used for OPLS-DA analysis, and, finally, ggplot2 was used to draw the map to obtain OPLS-DA analysis results. Based on the results of OPLS-DA, the variable importance in projection (VIP) and fold change in univariate analysis were combined to perform further screening for differential metabolites with threshold values of fold changes of ≥2 or ≤0.5 and VIP ≥ 1. The annotation results of the significantly different metabolites are classified by pathway enrichment in the KEGG database (https://www.kegg.jp/), illustrating all possible metabolic pathways.

3. Results and Discussion

In this study, we used metabolomics based on UPLC-MS/MS to determine the flavonoid profiles of tomato fruits. Methods for the accurate determination of flavonoid standards and their regression curves were established. A total of 50 flavonoid metabolites were identified in the tomato samples. Total ion current diagrams detected by mass spectrometry of QC samples are shown in Figure 2. Detailed information on the fifty flavonoid metabolites identified in the tomato fruit samples are in Table S1. MS/MS spectra of the fifty flavonoids identified in the tomato samples are in Supplementary File S1. UHPLC-MS/MS parameters are in Supplementary File S2. Calibration curves and quantitative details of the fifty flavonoids are shown in Supplementary File S3.

3.1. Comprehensive Analysis of Flavonoid Compounds

The hierarchical cluster analysis (HCA) results of samples and metabolites are presented as heatmaps with dendrograms. In HCA, the normalized signal intensities of metabolites (unit variance scaling) are visualized as a color spectrum. We performed a qualitative and quantitative determination of the metabolites of the flavonoids in tomato fruits; 50 flavonoid metabolites in the fruit samples were identified and quantitatively determined, including 6 chalcones, 8 flavonols, 4 flavanones, 5 flavanonols, 1 flavone glycosides, 8 flavones, 10 flavonols, 6 isoflavones, 1 xanthone, and 1 other flavonoid (Figure 3).

3.2. Differential Flavonoid Metabolites Based on Principal Component Analysis (PCA)

In the present study, the three principal components—PC1, PC2, and PC3—accounted for 82.69%, 6.2%, and 3.59% of the variability in the samples, respectively (Figure 4). The six different sample groups were clearly separated. The PCA results demonstrate a significant trend for the separation of the flavonoid metabolomes between the groups and indicate the difference in metabolites between each group of samples.

3.3. Screening and Analysis of Differential Flavonoid Metabolites

To analyze the trends in the changes in flavonoid metabolites during the fruit ripening in WT and NOR-like1 tomato fruit, differential flavonoid metabolites were screened and analyzed through a paired comparison. All flavonoids detected in tomato fruits were screened to identify differential metabolites with fold a change ≥2 or ≤0.5 and VIP ≥ 1, which are represented in volcano maps. The results of the paired comparison are shown in Table 1. The contents of 50 flavonoids at different stages are detailed in Table 2.

3.3.1. Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA) of Different Tomato Sample Groups

Pairwise comparisons were carried out using Orthogonal Partial Least Squares–Discriminant Analysis (OPLS-DA) models to clarify the metabolic differences observed among the different tomato sample groups. The parameters of the OPLS-DA evaluation model include R2X, R2Y, and Q2, where R2X and R2Y represent the explanatory power of the constructed model to the X and Y matrices. Q2 represents the predictive ability of the model. Q2 > 0.5 is generally considered to indicate an acceptable model and, theoretically, the closer the values of R2 and Q2 to 1, the better the model [31]. The OPLS-DA results of the pairwise comparison between groups are shown in Figure 5. In the results, the R2 values of all the comparison groups were above 0.5, which indicated a good fit for the model. Different tomato groups were distributed in different zones, implying that there were significant differences in their metabolites.

3.3.2. Changes in Flavonoid Metabolites in WT Tomato Fruits at Different Ripening Stages

The expression of metabolites varied at different growth stages of wild tomato fruits. Twelve differential metabolites (four upregulated and eight downregulated) were identified between WT-BR3 and WT-BR9 (Figure 6a). Nineteen differential metabolites (eighteen upregulated and one downregulated) were identified between WT-GR and WT-BR3 (Figure 6b).
The Venn diagram shows that there were seven common differential metabolites in the two combinations: 6,2′-dihydroxyflavone, dihydrokaempferol, phloretin, naringenin chalcone, trilobatin, isoliquiritigenin, and phlorizin (Figure 6c).
The individual metabolite contents are shown in Table 2. The total contents of all the flavonoids in WT-GR, WT-BR3, and WT-BR9 were 281.02 mg/kg DW, 276.25 mg/kg DW, and 223.37 mg/kg DW, respectively. Rutin, nicotiflorin, naringenin chalcone, eriodictyol, and naringenin-7-glucoside were the five flavonoids with the highest content in the ripening stages (BR3 and BR9) of the WT tomato fruits.
Rutin and nicotiflorin were the two richest components in WT-GR. Rutin changed insignificantly, while nicotiflorin decreased from GR to BR3. The content of naringenin chalcone increased to 28.46 mg/kg DW in BR3, indicating its biosynthesis was very active in this period.
The contents of some other flavonoids also changed during the ripening stages. Phlorizin and trilobatin increased significantly throughout the entire ripening process. The contents of nine flavonoids (luteolin, liquiritin, eriodictyol, narirutin, engeletin, genistein, glycitin, quercetin, and taxifolin) increased significantly from the GR stage to the BR3 stage but changed insignificantly from the BR3 stage to the BR9 stage. The contents of 6,2′-dihydroxyflavone, dihydrokaempferol, phloretin, naringenin chalcone, and isoliquiritigenin increased significantly from the GR stage to the BR3 stage but declined significantly from the BR3 stage to the BR9 stage. These results showed that the contents of a variety of flavonoids increased significantly in the GR stage to the BR3 stage, which indicated that flavonoid metabolite biosynthesis is very active in this period.
In contrast, in the later ripening process BR9, the differential flavonoid metabolites in WT-BR9 tomato fruits were mostly downregulated compared with those in WT-BR3. Astragalin declined significantly during the entire ripening process of WT tomato fruits. Narcissin, baimaside, and sieboldin changed little from GR to BR3 but decreased from BR3 to BR9. In WT-BR3, rutin, nicotiflorin, and naringenin chalcone were the three predominant components, while their contents all decreased from BR3 to BR9. Taken together, the decreased contents of the three main flavonoids and other downregulated flavonoid metabolites lead to a lower total flavonoid content in WT-BR9 than in WT-BR3.

3.3.3. Change in Flavonoid Metabolites in NOR-like1 Tomato Fruits at Different Ripening Stages

To investigate the changes in flavonoid metabolites among different ripening stages of NOR-like1 tomato fruits, differentially accumulated flavonoid metabolites were analyzed among NOR-like1-GR vs. NOR-like1-BR3 and NOR-like1-BR3 vs. NOR-like1-BR9. Eleven differential metabolites (eight upregulated and three downregulated) were identified between NOR-like1-BR3 and NOR-like1-BR9 (Figure 7a). Twenty-four differential metabolites (seventeen upregulated and seven downregulated) were identified between NOR-like1-GR and NOR-like1-BR3 (Figure 7b).
The Venn diagram shows that two combinations have five common differential metabolites (genistein, naringenin-7-glucoside, phloretin, trilobatin, and phlorizin) (Figure 7c).
The total contents of all the flavonoids in NOR-like1-GR, NOR-like1-BR3, and NOR-like1-BR9 were 176.87 mg/kg DW, 214.39 mg/kg DW, and 245.93 mg/kg DW, respectively. Rutin, nicotiflorin, naringenin chalcone, eriodictyol, and naringenin-7-glucoside were also the five flavonoids with the highest content in the ripening stages (BR3 and BR9) in NOR-like1 tomato fruits.
Rutin and nicotiflorin were the two dominant flavonoids among the NOR-like1-GR fruits. From GR to BR3, the contents of rutin and naringenin chalcone increased to 125.43 mg/kg DW and 29.6 mg/kg DW in BR3, respectively, while nicotiflorin decreased from 68.36 mg/kg DW to 46.47 mg/kg DW. Moreover, there were seventeen upregulated flavonoid metabolites between NOR-like1-GR and NOR-like1-BR3 (Figure 7b). The results showed that flavonoid metabolites accumulated in large amounts from the GR stage to the BR3 stage in NOR-like1 tomato fruits. During this period, flavonoid metabolite biosynthesis is very active. The increased contents of the two main flavonoids (rutin and naringenin chalcone) and other upregulated flavonoid metabolites lead to a higher total of flavonoid content in NOR-like1-BR3 than that in NOR-like1-GR.
From BR3 to BR9, the contents of rutin, nicotiflorin, eriodictyol, and naringenin-7-glucoside increased to 145.01 mg/kg DW, 60.35 mg/kg DW, 11.84 mg/kg DW, and 3.72 mg/kg DW, respectively. The content of these four components comprised more than 90% of the total content of all the flavonoids in NOR-like1-BR9. Naringenin chalcone decreased slightly from 29.6 to 21.6 mg/kg DW. The contents of the other flavonoids also changed from BR3 to BR9. Among them, eight were upregulated and three were downregulated. The increased contents of the four dominant flavonoids and the eight upregulated flavonoids led to an increase in the total flavonoid content in NOR-like1-BR9 than that in NOR-like1-BR3.

3.3.4. Differential Flavonoid Metabolites between WT and NOR-like1 Tomato Fruits

In order to explore the differences in the flavonoid metabolites at different growth stages of the two varieties of tomato fruits, differential metabolites were screened using a pairwise comparison. Figure 8a–c presents Volcano plots showing the levels of differential metabolites between WT and NOR-like1 at different maturity stages. Detailed information about the contents of flavonoid metabolites detected in each tomato fruit sample is shown in Table 2.
During the BR3 stage, nine differential flavonoids (six upregulated and three downregulated) were identified between WT-BR3 and NOR-like1-BR3 (Figure 8a). During the BR9 stage, 12 differential metabolites (11 upregulated and 1 downregulated) were identified between WT-BR9 and NOR-like1-BR9 (Figure 8b). During the GR stage, 16 differential metabolites (13 upregulated and 3 downregulated) were identified between WT-GR and NOR-like1-GR (Figure 8c).
The Venn diagram shows that there were two common differential metabolites (glycitin and apigenin-7-glucuronide) in the three combinations (Figure 8d), three common differential metabolites (glycitin, apigenin-7-glucuronide, and cynaroside) in the combination of WT-GR_vs_NOR-like1-GR and WT-BR3_vs_NOR-like1-BR3, four common differential metabolites (glycitin, naringenin chalcone, apigenin-7-glucuronide, and calycosin-7-O-β-D-glucoside) in the combination of WT-BR9_vs_NOR-like1-BR9 and WT-BR3_vs_NOR-like1-BR3, and four common differential metabolites (glycitin, apigenin-7-glucuronide, narirutin, and (-)-catechin) in the combination of WT-GR_vs_NOR-like1-GR and WT-BR9_vs_NOR-like1-BR9. Apigenin-7-glucuronide was detected only in NOR-like1 tomato fruits in all three maturation stages and not in the WT. The contents of glycitin in NOR-like1 tomato fruits were significantly higher than those in the WT in all three maturation stages.
Compared to the WT tomato fruits, most of the flavonoids were significantly changed in NOR-like1 tomato fruits during the ripening stages, indicating that SlNOR-like1 does regulate the biosynthesis of flavonoids. In our previous study, by using untargeted metabolomics analysis to study the metabolic profile of NOR-like1 gene-edited tomato fruits with WT tomato fruits for comparison, it was found that flavonoid metabolites changed significantly during the whole ripening stages [28], which corresponded to the results here.
The total amount of all the flavonoids in NOR-like1 tomato fruits was 176.87 mg/kg DW in the GR stage and increased to 214.39 mg/kg DW in BR3 and 245.93 mg/kg DW in BR9. The total amount of all the flavonoids in WT tomato fruits was 281.02 mg/kg DW in the GR stage, which slightly decreased to 276.27 mg/kg DW in BR3 and then to 223.37 mg/kg DW in BR9. The total amount of all the flavonoids in NOR-like1 tomato fruits was lower than in WT tomato fruits during the GR stage and BR3 stage. However, during the BR9 stage, the total amount of all the flavonoids in NOR-like1 tomato fruits was higher than in WT tomato fruits. The dynamic changes of the flavonoids differed among the two different tomato strains.
To explore the reason causing such a difference, we compared the contents of the predominant flavonoid compounds and their changes in two different tomato strains and found that rutin, nicotiflorin, naringenin chalcone, eriodictyol, and naringenin-7-glucoside were the five flavonoids with the highest content in the ripening stages (BR3 and BR9) in both WT and NOR-like1 tomato fruits. However, the accumulation patterns of some flavonoids differed between the WT and NOR-like1 tomato fruits. Rutin was the most abundant and significantly variable. Rutin has been extensively detected in tomatoes in multiple literature studies [32,33,34,35,36,37,38,39]. In our study, the content of rutin in WT tomato fruits was highest during the GR stage (175.26 mg/kg DW), and its content in NOR-like1-GR tomato fruits was 101.18 mg/kg DW. However, during the subsequent growth process, the content of rutin gradually increased in NOR-like1 tomato fruits but decreased in the WT. In BR9, the content of rutin in NOR-like1 tomato fruits was 145.01 mg/kg DW and, in WT tomato fruits, was 147.79 mg/kg DW, which was very close. In our research, rutin was the most abundant flavonoid in both WT tomato fruits and NOR-like1 tomato fruits; however, the accumulation patterns in the two varieties were different.
Naringenin chalcone and rutin were detected in large quantities in different varieties of tomatoes [40,41]. Similarly, in our study, we found that the content of naringenin chalcone and rutin in WT tomato fruits was different from that in NOR-like1 tomato fruits but were both very high. The content of naringenin chalcone in NOR-like1 tomato fruit was highest during the BR3 stage (29.6 mg/kg DW); the content of naringenin chalcone in WT-BR3 tomato fruits was 28.46 mg/kg DW. When in the BR9 growth stage, the content of naringenin chalcone in NOR-like1 tomato fruits was 21.6 mg/kg DW, which was much higher than that in WT tomato fruit (10.67 mg/kg DW). The content of naringenin chalcone decreased from the BR3 stage to the BR9 stage. In our study, six types of chalcones showed significant changes, including isoliquirigenin, philorizin, trilobatin, philoretin, naringenin chalcone, and sieboldin. Among them, naringenin chalcone showed the most significant changes and had the highest content. Iijima et al. (2008) demonstrated that naringenin chalcone was the major chalcone in tomato fruit, revealing the highest accumulation levels at the breaker stage, which then decreased during ripening [42]. Similarly, Raffo et al. (2002) also mentioned that the content of naringenin gradually decreased with maturity [43]. Naringenin chalcone and rutin are considered to be the most abundant substances based on a description by Muir et al. (2001) [44].
The content of nicotiflorin in WT tomato fruits was highest during the GR stage (102.77 mg/kg DW), and the total content of nicotiflorin in NOR-like1-GR tomato fruits was 68.37 mg/kg DW. When in the BR3 growth stage, the content of nicotiflorin in NOR-like1 tomato fruit was 46.47 mg/kg DW, and the content of nicotiflorin in WT-BR3 tomato fruit was 65.96 mg/kg DW. When in the BR9 growth stage, the content of nicotiflorin NOR-like1 tomato fruits was 60.35 mg/kg DW; meanwhile, the content of nicotiflorin WT tomato fruits decreased to 53.63 mg/kg DW, lower than the NOR-like1 tomato fruits. Nicotiflorin was also extensively detected in tomatoes in a study by Le Gall et al. (2003) [45].
The content of eriodictyol increased during the GR stage to the BR3 stage, with 4.89 mg/kg DW in WT-BR3 and 7.18 mg/kg DW in NOR-like1-BR3, decreased slightly to 4.66 mg/kg DW in WT-BR9, but increased to 11.84 mg/kg DW in NOR-like1-BR9. For tomatoes in the same growth period, the content of eriodictol in NOR-like1 tomato fruits was significantly higher than that in the WT. Mintz Oron et al. (2008) found that tomatoes contain a large amount of chalconingenin, naringenin, eriodictyol, and nicotiflorin [46].
From the above data comparison, the accumulation pattern of flavonoid metabolites in NOR-like1 tomato fruits differed from that in WT tomato fruits, especially in the later ripening process BR9. The contents of the four dominant flavonoids (rutin, nicotiflorin, eriodictyol, and naringenin-7-glucoside) increased, and there were eight upregulated flavonoids in NOR-like1-BR9 than that in NOR-like1-BR3. The differential flavonoid metabolites in WT-BR9 tomato fruits were mostly downregulated compared with those in WT-BR3. As the three predominant components in WT-BR3, the contents of rutin, nicotiflorin, and naringenin chalcone decreased from BR3 to BR9. The transcriptomics results in our previous studies [28] showed that the expression levels of key genes involved in the biosynthesis of flavonoid metabolites were upregulated significantly during BR9 in NOR-like1 tomatoes compared to WT tomatoes, which may provide some support to show that NOR-like1 influences the metabolomics of flavonoid metabolites in tomato fruits at the level of gene expression.

3.4. Enrichment Analysis of Differential Flavonoid Metabolites

The results of further KEGG enrichment analysis are shown in Figure 9a–f. It is revealed that the significant enrichment of the differential flavonoid metabolites was distributed in the following pathways: flavonoid biosynthesis, flavone and flavonol biosynthesis, the biosynthesis of secondary metabolites, metabolic pathways, and isoflavonoid biosynthesis. The most enriched pathway was the ‘biosynthesis of secondary metabolites’. This was followed by ‘flavonoid biosynthesis’.
In NOR-like1-BR3_vs_NOR-like1-BR9 (Figure 9a), phlorizin and phloretin were enriched in the flavonoid biosynthesis pathway; genistein was enriched in the isoflavonoid biosynthesis. In NOR-like1-GR_vs_NOR-like1-BR3 (Figure 9b), quercitrin, quercetin, and luteolin were enriched in flavone and flavonol biosynthesis; isoliquiritigenin, phlorizin, phloretin, naringenin chalcone, (-)-epicatechin, (-)-epigallocatechin, eriodictyol, dihydrokaempferol, taxifolin, luteolin, and quercetin were related to ‘flavonoid biosynthesis’; genistein, glycitin, and 2′-hydroxygenistein were enriched in isoflavonoid biosynthesis. In WT-BR3_vs_WT-BR9 (Figure 9c), baimaside was enriched in flavone and flavonol biosynthesis. Isoliquiritigenin, phlorizin, phloretin, naringenin chalcone, and dihydrokaempferol were enriched in flavonoid biosynthesis. In WT-GR_vs_WT-BR3 (Figure 9d), luteolin, quercetin, and astragalin were enriched in flavone and flavonol biosynthesis. Isoliquiritigenin, phlorizin, phloretin, naringenin chalcone, eriodictyol, dihydrokaempferol, taxifolin, luteolin, and quercetin were enriched in flavonoid biosynthesis. Genistein, glycitin, and 2′-hydroxygenistein were related to ‘isoflavonoid biosynthesis’. In WT-GR_vs_NOR-like1-GR (Figure 9e), astragalin and quercitrin were enriched in flavone and flavonol biosynthesis. Phlorizin, phloretin, naringenin chalcone, (-)-epicatechin, and (-)-epigallocatechin were enriched in flavonoid biosynthesis. The regulation of ononin and glycitin was related to isoflavonoid biosynthesis. In WT-BR9_vs_NOR-like1-BR9 (Figure 9f), luteolin was enriched in flavone and flavonol biosynthesis. Isoliquiritigenin, naringenin chalcone, eriodictyol, and luteolin were enriched in flavonoid biosynthesis. Glycitin was related to ‘isoflavonoid biosynthesis’. In all the above comparison groups, flavonoid biosynthesis was significantly enriched, which may indicate that the differential flavonoids in these enrichment pathways may lead to changes in flavonoid accumulation between the two types of tomato fruits during the whole ripening process.

4. Conclusions

Flavonoid quantitative metabolomics was used to accurately quantify and analyze the changes in the flavonoid metabolites of wild-type tomato fruits (WT) and gene-edited tomato fruits (NOR-like1) during the tomato ripening process. A total of 50 flavonoid metabolites were identified from the WT and NOR-like1 tomato fruits using UPLC-MS/MS. Rutin, nicotiflorin, naringenin chalcone, eriodictyol, and naringenin-7-glucoside were the five flavonoids with the highest content in the ripening stages (BR3 and BR9) in both WT and NOR-like1 tomato fruits. The overall flavonoid contents in WT tomato fruits changed little from GR to BR3 and decreased from BR3 to BR9; meanwhile, in NOR-like1 tomato fruits, the total amounts of the flavonoids exhibited an increasing trend during all three ripening stages. The accumulation pattern of flavonoid metabolites in NOR-like1 tomato fruits differed from that in WT tomato fruits, especially in the later ripening process BR9. These findings highlight significant differences in the profile of flavonoid metabolites between WT and NOR-like1 tomato samples and showed that NOR-like1 regulated the biosynthesis of flavonoids and changed the accumulation pattern of flavonoids. Overall, our research provides new insights into the dynamic changes in flavonoid profiles caused by transcription factor NOR-like1 during ripening. Our future studies may focus on the flavonoids most significantly affected, extensively investigate the interaction mechanism between transcription factor NOR-like1 and the related genes, and elucidate the regulatory mechanism of NOR-like1 on flavonoids.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/foods12244445/s1, Table S1. Detailed information on the fifty flavonoid metabolites identified in the tomato fruit samples. Supplementary File S1. MS/MS spectra of the fifty flavonoids identified in the tomato samples. Supplementary File S2. UHPLC-MS/MS parameters of the fifty flavonoids. Supplementary File S3. Calibration curves and quantitative details of the fifty flavonoids.

Author Contributions

Conceptualization, supervision, funding acquisition: X.Z.; formal analysis, writing—original draft preparation, visualization: D.G.; resources: Y.Z.; investigation resources: D.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (National Natural Science Foundation of China. No. 31972129).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data supporting the results of this study are included in the present article.

Acknowledgments

We appreciate Wuhan Metware Biotechnology Co. for providing technical support and China Agricultural University for providing the tomato material.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gao, J.; Song, Y.; Mao, J. The determination of flavonoids in some fruit and vegetables. China Food Addit. 2007, 18, 211–213. [Google Scholar] [CrossRef]
  2. Treutter, D. Significance of flavonoids in plant resistance: A review. Environ. Chem. Lett. 2006, 4, 147–157. [Google Scholar] [CrossRef]
  3. Bui, T.T.; Nguyen, T.H. Natural product for the treatment of Alzheimer’s disease. J. Basic Clin. Physiol. Pharmacol. 2017, 28, 413–423. [Google Scholar] [CrossRef] [PubMed]
  4. Tan, T.; Li, Y.; Tang, B.; Chen, Y.; Chen, X.; Xie, Q.; Hu, Z.; Chen, G. Knockout of SlALKBH2 weakens the DNA damage repair ability of tomato. Plant Sci. 2022, 319, 111266. [Google Scholar] [CrossRef] [PubMed]
  5. Tao, X.; Wu, Q.; Huang, S.; Zhu, B.; Chen, F.; Liu, B.; Cai, L.; Mao, L.; Luo, Z.; Li, L.; et al. Exogenous abscisic acid regulates primary metabolism in postharvest cherry tomato fruit during ripening. Sci. Hortic. 2022, 299, 111008. [Google Scholar] [CrossRef]
  6. Monribot-Villanueva, J.L.; Altúzar-Molina, A.; Aluja, M.; Zamora-Briseño, J.A.; Elizalde-Contreras, J.M.; Bautista-Valle, M.V.; Arellano de Los Santos, J.; Sánchez-Martínez, D.E.; Rivera-Reséndiz, F.J.; Vázquez-Rosas-Landa, M.; et al. Integrating proteomics and metabolomics approaches to elucidate the ripening process in white Psidium guajava. Food Chem. 2022, 367, 130656. [Google Scholar] [CrossRef] [PubMed]
  7. Yonekura-Sakakibara, K.; Higashi, Y.; Nakabayashi, R. The Origin and Evolution of Plant Flavonoid Metabolism. Front. Plant Sci. 2019, 10, 943. [Google Scholar] [CrossRef] [PubMed]
  8. Xiao, J.; Muzashvili, T.S.; Georgiev, M.I. Advances in the biotechnological glycosylation of valuable flavonoids. Biotechnol. Adv. 2014, 32, 1145–1156. [Google Scholar] [CrossRef]
  9. Liu, Z.-Q. What about the progress in the synthesis of flavonoid from 2020? Eur. J. Med. Chem. 2022, 243, 114671. [Google Scholar] [CrossRef]
  10. Zhang, X.; Qi, X.; Chen, Y.; Yao, H.; Zheng, Y. Research Advances of Flavonoids in Strawberry. North. Hortic. 2020, 44, 128–133. [Google Scholar] [CrossRef]
  11. Tejada, S.; Pinya, S.; Martorell, M.; Capó, X.; Tur, J.A.; Pons, A.; Sureda, A. Potential Anti-inflammatory Effects of Hesperidin from the Genus Citrus. Curr. Med. Chem. 2018, 25, 4929–4945. [Google Scholar] [CrossRef]
  12. Shen, N.; Wang, T.; Gan, Q.; Liu, S.; Wang, L.; Jin, B. Plant flavonoids: Classification, distribution, biosynthesis, and antioxidant activity. Food Chem. 2022, 383, 132531. [Google Scholar] [CrossRef] [PubMed]
  13. Staszowska-Karkut, M.; Materska, M. Phenolic Composition, Mineral Content, and Beneficial Bioactivities of Leaf Extracts from Black Currant (Ribes nigrum L.), Raspberry (Rubus idaeus), and Aronia (Aronia melanocarpa). Nutrients 2020, 12, 463. [Google Scholar] [CrossRef] [PubMed]
  14. Tohge, T.; de Souza, L.P.; Fernie, A.R. Current understanding of the pathways of flavonoid biosynthesis in model and crop plants. J. Exp. Bot. 2017, 68, 4013–4028. [Google Scholar] [CrossRef] [PubMed]
  15. Song, X.; Wei, J.; Di, S.; Pang, Y. Recent Advances in the Regulation Mechanism of Transcription Factors and Metabolic Engineering of Anthocyanins. Chin. Bull. Bot. 2019, 54, 133–156. [Google Scholar] [CrossRef]
  16. Huang, Y.; Cui, L.; Chen, W.; Liu, Z.; Yuan, W.; Zhu, F.; Jiao, Z.; Zhang, Z.; Deng, X.; Wang, L.; et al. Comprehensive analysis of NAC transcription factors and their expressions during taproot coloration in radish (Raphanus sativus L.). Sci. Hortic. 2022, 299, 111047. [Google Scholar] [CrossRef]
  17. Sun, Q.; Jiang, S.; Zhang, T.; Xu, H.; Fang, H.; Zhang, J.; Su, M.; Wang, Y.; Zhang, Z.; Wang, N.; et al. Apple NAC transcription factor MdNAC52 regulates biosynthesis of anthocyanin and proanthocyanidin through MdMYB9 and MdMYB11. Plant Sci. 2019, 289, 110286. [Google Scholar] [CrossRef] [PubMed]
  18. Duan, A.-Q.; Yang, X.-L.; Feng, K.; Liu, J.-X.; Xu, Z.-S.; Xiong, A.-S. Genome-wide analysis of NAC transcription factors and their response to abiotic stress in celery (Apium graveolens L.). Comput. Biol. Chem. 2020, 84, 107186. [Google Scholar] [CrossRef]
  19. Hu, H.; Dai, M.; Yao, J.; Xiao, B.; Li, X.; Zhang, Q.; Xiong, L. Overexpressing a NAM, ATAF, and CUC (NAC) transcription factor enhances drought resistance and salt tolerance in rice. Proc. Natl. Acad. Sci. USA 2006, 103, 12987–12992. [Google Scholar] [CrossRef]
  20. Jeong, J.S.; Kim, Y.S.; Baek, K.H.; Jung, H.; Ha, S.-H.; Do Choi, Y.; Kim, M.; Reuzeau, C.; Kim, J.-K. Root-Specific Expression of OsNAC10 Improves Drought Tolerance and Grain Yield in Rice under Field Drought Conditions. Plant Physiol. 2010, 153, 185–197. [Google Scholar] [CrossRef]
  21. Guo, F.; Liu, S.; Zhang, C.; Dong, T.; Meng, X.; Zhu, M. Genome-wide systematic survey and analysis of NAC transcription factor family and their response to abiotic stress in sweetpotato. Sci. Hortic. 2022, 299, 111048. [Google Scholar] [CrossRef]
  22. Li, M.; Chen, R.; Jiang, Q.; Sun, X.; Zhang, H.; Hu, Z. GmNAC06, a NAC domain transcription factor enhances salt stress tolerance in soybean. Plant Mol. Biol. 2021, 105, 333–345. [Google Scholar] [CrossRef] [PubMed]
  23. Ma, J.; Wang, L.-Y.; Dai, J.-X.; Wang, Y.; Lin, D. The NAC-type transcription factor CaNAC46 regulates the salt and drought tolerance of transgenic Arabidopsis thaliana. BMC Plant Biol. 2021, 21, 11. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, L.; Hu, Z.; Zhu, M.; Zhu, Z.; Hu, J.; Qanmber, G.; Chen, G. The abiotic stress-responsive NAC transcription factor SlNAC11 is involved in drought and salt response in tomato (Solanum lycopersicum L.). Plant Cell Tissue Organ Cult. (PCTOC) 2017, 129, 161–174. [Google Scholar] [CrossRef]
  25. Zhu, M.; Chen, G.; Zhang, J.; Zhang, Y.; Xie, Q.; Zhao, Z.; Pan, Y.; Hu, Z. The abiotic stress-responsive NAC-type transcription factor SlNAC4 regulates salt and drought tolerance and stress-related genes in tomato (Solanum lycopersicum). Plant Cell Rep. 2014, 33, 1851–1863. [Google Scholar] [CrossRef] [PubMed]
  26. Gao, Y.; Wei, W.; Zhao, X.; Tan, X.; Fan, Z.; Zhang, Y.; Jing, Y.; Meng, L.; Zhu, B.; Zhu, H.; et al. A NAC transcription factor, NOR-like1, is a new positive regulator of tomato fruit ripening. Hortic. Res. 2018, 5, 75. [Google Scholar] [CrossRef] [PubMed]
  27. Gao, Y.; Wei, W.; Fan, Z.; Zhao, X.; Zhang, Y.; Jing, Y.; Zhu, B.; Zhu, H.; Shan, W.; Chen, J.; et al. Re-evaluation of the nor mutation and the role of the NAC-NOR transcription factor in tomato fruit ripening. J. Exp. Bot. 2020, 71, 3560–3574. [Google Scholar] [CrossRef] [PubMed]
  28. Yang, X.; Zhao, X.; Fu, D.; Zhao, Y. Integrated Analysis of Widely Targeted Metabolomics and Transcriptomics Reveals the Effects of Transcription Factor NOR-like1 on Alkaloids, Phenolic Acids, and Flavonoids in Tomato at Different Ripening Stages. Metabolites 2022, 12, 1296. [Google Scholar] [CrossRef]
  29. Zhang, H.; Du, X.; Yu, J.; Jin, H.; Liu, N. Comparative Metabolomics study of flavonoids in the pericarp of different coloured bitter gourds (Momordica charantia L.). Physiol. Mol. Biol. Plants 2022, 28, 1347–1357. [Google Scholar] [CrossRef]
  30. Wang, Y.; Cheng, J.; Jiang, W.; Chen, S. Metabolomics study of flavonoids in Coreopsis tinctoria of different origins by UPLC-MS/MS. PeerJ 2022, 10, e14580. [Google Scholar] [CrossRef]
  31. Li, J.; Hossain, M.S.; Ma, H.; Yang, Q.; Gong, X.; Yang, P.; Feng, B. Comparative metabolomics reveals differences in flavonoid metabolites among different coloured buckwheat flowers. J. Food Compos. Anal. 2020, 85, 103335. [Google Scholar] [CrossRef]
  32. Hunt, G.M.; Baker, E.A. Phenolic constituents of tomato fruit cuticles. Phytochemistry 1980, 19, 1415–1419. [Google Scholar] [CrossRef]
  33. Iijima, Y.; Nakamura, Y.; Ogata, Y.; Tanaka, K.; Sakurai, N.; Suda, K.; Suzuki, T.; Suzuki, H.; Okazaki, K.; Kitayama, M.; et al. Metabolite annotations based on the integration of mass spectral information. Plant J. 2008, 54, 949–962. [Google Scholar] [CrossRef] [PubMed]
  34. Krishnan, M.; Babu, S.; Thomas, S.A.; Surulivel, J.S.; Ayyanar, K. Molecular docking analysis of VEGF with compounds from tomato. Bioinformation 2022, 18, 478–481. [Google Scholar] [CrossRef] [PubMed]
  35. Spencer, J.P.; Kuhnle, G.G.; Hajirezaei, M.; Mock, H.P.; Sonnewald, U.; Rice-Evans, C. The genotypic variation of the antioxidant potential of different tomato varieties. Free Radic. Res. 2005, 39, 1005–1016. [Google Scholar] [CrossRef] [PubMed]
  36. Slimestad, R.; Fossen, T.; Verheul, M.J. The flavonoids of tomatoes. J. Agric. Food Chem. 2008, 56, 2436–2441. [Google Scholar] [CrossRef] [PubMed]
  37. Vallverdú-Queralt, A.; Jáuregui, O.; Di Lecce, G.; Andrés-Lacueva, C.; Lamuela-Raventós, R.M. Screening of the polyphenol content of tomato-based products through accurate-mass spectrometry (HPLC-ESI-QTOF). Food Chem. 2011, 129, 877–883. [Google Scholar] [CrossRef] [PubMed]
  38. Le Gall, G.; DuPont, M.S.; Mellon, F.A.; Davis, A.L.; Collins, G.J.; Verhoeyen, M.E.; Colquhoun, I.J. Characterization and content of flavonoid glycosides in genetically modified tomato (Lycopersicon esculentum) fruits. J. Agric. Food Chem. 2003, 51, 2438–2446. [Google Scholar] [CrossRef]
  39. Larbat, R.; Le Bot, J.; Bourgaud, F.; Robin, C.; Adamowicz, S. Organ-specific responses of tomato growth and phenolic metabolism to nitrate limitation. Plant Biol. 2012, 14, 760–769. [Google Scholar] [CrossRef]
  40. Morelli, C.F.; Cutignano, A.; Speranza, G.; Abbamondi, G.R.; Rabuffetti, M.; Iodice, C.; De Prisco, R.; Tommonaro, G. Taste Compounds and Polyphenolic Profile of Tomato Varieties Cultivated with Beneficial Microorganisms: A Chemical Investigation on Nutritional Properties and Sensory Qualities. Biomolecules 2023, 13, 117. [Google Scholar] [CrossRef]
  41. Wang, H.; Zhang, Z.; Song, J.; Tian, M.; Li, R.; Cui, X. Phenolic compound identification in tomato fruit by UPLC-QTOF-MS. LWT 2023, 182, 114791. [Google Scholar] [CrossRef]
  42. Iijima, Y.; Suda, K.; Suzuki, T.; Aoki, K.; Shibata, D. Metabolite Profiling of Chalcones and Flavanones in Tomato Fruit. J. Jpn. Soc. Hortic. Sci. 2008, 77, 94–102. [Google Scholar] [CrossRef]
  43. Raffo, A.; Leonardi, C.; Fogliano, V.; Ambrosino, P.; Salucci, M.; Gennaro, L.; Bugianesi, R.; Giuffrida, F.; Quaglia, G. Nutritional value of cherry tomatoes (Lycopersicon esculentum Cv. Naomi F1) harvested at different ripening stages. J. Agric. Food Chem. 2002, 50, 6550–6556. [Google Scholar] [CrossRef] [PubMed]
  44. Muir, S.R.; Collins, G.J.; Robinson, S.; Hughes, S.; Bovy, A.; Ric De Vos, C.H.; van Tunen, A.J.; Verhoeyen, M.E. Overexpression of petunia chalcone isomerase in tomato results in fruit containing increased levels of flavonols. Nat. Biotechnol. 2001, 19, 470–474. [Google Scholar] [CrossRef] [PubMed]
  45. Le Gall, G.; Colquhoun, I.J.; Davis, A.L.; Collins, G.J.; Verhoeyen, M.E. Metabolite profiling of tomato (Lycopersicon esculentum) using 1H NMR spectroscopy as a tool to detect potential unintended effects following a genetic modification. J. Agric. Food Chem. 2003, 51, 2447–2456. [Google Scholar] [CrossRef]
  46. Mintz-Oron, S.; Mandel, T.; Rogachev, I.; Feldberg, L.; Lotan, O.; Yativ, M.; Wang, Z.; Jetter, R.; Venger, I.; Adato, A.; et al. Gene expression and metabolism in tomato fruit surface tissues. Plant Physiol. 2008, 147, 823–851. [Google Scholar] [CrossRef]
Foods 12 04445 g001
Figure 1. Representative flavonoid derivatives and relationships among their structures.
Figure 1. Representative flavonoid derivatives and relationships among their structures.
Foods 12 04445 g001
Foods 12 04445 g002
Figure 2. Mass spectrometry quality control (QC) total ion current (TIC) superposition diagram: (a) negative-ion TIC diagram, (b) positive-ion mode TIC diagram.
Figure 2. Mass spectrometry quality control (QC) total ion current (TIC) superposition diagram: (a) negative-ion TIC diagram, (b) positive-ion mode TIC diagram.
Foods 12 04445 g002
Foods 12 04445 g003
Figure 3. Heatmap visualization of all 50 detected flavonoid metabolites. Each sample has its own column, while each flavonoid is represented by one row. The scale represents the content level obtained after standardization (the redder the color, the higher the content).
Figure 3. Heatmap visualization of all 50 detected flavonoid metabolites. Each sample has its own column, while each flavonoid is represented by one row. The scale represents the content level obtained after standardization (the redder the color, the higher the content).
Foods 12 04445 g003
Foods 12 04445 g004
Figure 4. PCA analysis of the flavonoid metabolites in the tomato fruits of the six sample groups (3D image). Points with different colors represent different sample groups.
Figure 4. PCA analysis of the flavonoid metabolites in the tomato fruits of the six sample groups (3D image). Points with different colors represent different sample groups.
Foods 12 04445 g004
Foods 12 04445 g005
Figure 5. OPLS-DA model plots and loading plots of pairwise comparisons between the different tomato sample groups: (a) NOR-like1-BR3 vs. NOR-like1-BR9 (R2X = 0.786, R2Y = 0.993, Q2 = 0.63); (b) NOR-like1-GR vs. NOR-like1-BR3 (R2X = 0.626, R2Y = 0.999, Q2 = 0.822); (c) WT-BR3 vs. WT-BR9 (R2X = 0.594, R2Y = 0.989, Q2 = 0.771); (d) WT-GR vs. WT-BR3 (R2X = 0.75, R2Y = 0.999, Q2 = 0.943); (e) WT-BR3 vs. NOR-like1-BR3 (R2X = 0.755, R2Y = 0.955, Q2 = 0.632); (f) WT-BR9 vs. NOR-like1-BR9 (R2X = 0.585, R2Y = 0.998, Q2 = 0.752); (g) WT-GR vs. NOR-like1-GR (R2X = 0.726, R2Y = 0.999, Q2 = 0.891).
Figure 5. OPLS-DA model plots and loading plots of pairwise comparisons between the different tomato sample groups: (a) NOR-like1-BR3 vs. NOR-like1-BR9 (R2X = 0.786, R2Y = 0.993, Q2 = 0.63); (b) NOR-like1-GR vs. NOR-like1-BR3 (R2X = 0.626, R2Y = 0.999, Q2 = 0.822); (c) WT-BR3 vs. WT-BR9 (R2X = 0.594, R2Y = 0.989, Q2 = 0.771); (d) WT-GR vs. WT-BR3 (R2X = 0.75, R2Y = 0.999, Q2 = 0.943); (e) WT-BR3 vs. NOR-like1-BR3 (R2X = 0.755, R2Y = 0.955, Q2 = 0.632); (f) WT-BR9 vs. NOR-like1-BR9 (R2X = 0.585, R2Y = 0.998, Q2 = 0.752); (g) WT-GR vs. NOR-like1-GR (R2X = 0.726, R2Y = 0.999, Q2 = 0.891).
Foods 12 04445 g005
Foods 12 04445 g006
Figure 6. The differential metabolite analysis results in different ripening stages in WT tomato fruit: (a,b) The volcano plot shows the differential metabolite expression levels between WT-GR vs. WT-BR3 and WT-BR3 vs. WT-BR9, respectively. (c)The Venn diagram shows the differential metabolites from different ripening stages in WT tomato fruits.
Figure 6. The differential metabolite analysis results in different ripening stages in WT tomato fruit: (a,b) The volcano plot shows the differential metabolite expression levels between WT-GR vs. WT-BR3 and WT-BR3 vs. WT-BR9, respectively. (c)The Venn diagram shows the differential metabolites from different ripening stages in WT tomato fruits.
Foods 12 04445 g006
Foods 12 04445 g007
Figure 7. The differential metabolite analysis results at different ripening stages in NOR-like1 tomato fruits: (a,b) Volcano plot showing the differential metabolite expression levels among NOR-like1-GR vs. NOR-like1-BR3 and NOR-like1-BR3 vs. NOR-like1-BR9, respectively. (c) Venn diagram showing the differential metabolites from different ripening stages in NOR-like1 tomato fruits.
Figure 7. The differential metabolite analysis results at different ripening stages in NOR-like1 tomato fruits: (a,b) Volcano plot showing the differential metabolite expression levels among NOR-like1-GR vs. NOR-like1-BR3 and NOR-like1-BR3 vs. NOR-like1-BR9, respectively. (c) Venn diagram showing the differential metabolites from different ripening stages in NOR-like1 tomato fruits.
Foods 12 04445 g007
Foods 12 04445 g008
Figure 8. Analysis of differential metabolites between WT and NOR-like1 tomato fruits at the same ripening stage: (a–c) The volcano plot shows the expression levels of the differential metabolites among WT-BR3 vs. NOR-like1-BR3, WT-BR9 vs. NOR-like1-BR9, and WT-GR vs. NOR-like1-GR, respectively. (d) The Venn diagram shows the differential metabolites between WT and NOR-like1 tomato fruits at the same ripening stage.
Figure 8. Analysis of differential metabolites between WT and NOR-like1 tomato fruits at the same ripening stage: (a–c) The volcano plot shows the expression levels of the differential metabolites among WT-BR3 vs. NOR-like1-BR3, WT-BR9 vs. NOR-like1-BR9, and WT-GR vs. NOR-like1-GR, respectively. (d) The Venn diagram shows the differential metabolites between WT and NOR-like1 tomato fruits at the same ripening stage.
Foods 12 04445 g008
Foods 12 04445 g009aFoods 12 04445 g009b
Figure 9. KEGG enrichment analysis of the differential metabolites in the comparisons of (a) NOR-like1-BR3 vs. NOR-like1-BR9, (b) NOR-like1-GR vs. NOR-like1-BR3, (c) WT-BR3 vs. WT-BR9, (d) WT-GR vs. WT-BR3, (e) WT-GR vs. NOR-like1-GR, and (f) WT-BR9 vs. NOR-like1-BR9.
Figure 9. KEGG enrichment analysis of the differential metabolites in the comparisons of (a) NOR-like1-BR3 vs. NOR-like1-BR9, (b) NOR-like1-GR vs. NOR-like1-BR3, (c) WT-BR3 vs. WT-BR9, (d) WT-GR vs. WT-BR3, (e) WT-GR vs. NOR-like1-GR, and (f) WT-BR9 vs. NOR-like1-BR9.
Foods 12 04445 g009aFoods 12 04445 g009b
Table 1. Results of pairwise comparison of fifty flavonoid metabolites identified in the six tomato sample groups.
Table 1. Results of pairwise comparison of fifty flavonoid metabolites identified in the six tomato sample groups.
IndexCompoundsWT-GR vs. NOR-like1-GRWT-BR3 vs. NOR-like1-BR3WT-BR9 vs. NOR-like1-BR9NOR-like1-GR vs.
NOR-like1-BR3		NOR-like1-BR3 vs.
NOR-like1-BR9		WT-GR vs. WT-BR3WT-BR3 vs. WT-BR9
VIPLog2FCVIPLog2FCVIPLog2FCVIPLog2FCVIPLog2FCVIPLog2FCVIPLog2FC
1Astilbin--------------
2Miquelianin--------------
9Genistein------1.48inf1.411.381.38inf--
106Typhaneoside--1.38−1.27----------
11Hydroxysafflor yellow A----1.241.111.191.35------
115Quercitrin1.082.77----1.24−4.87------
117Engeletin--1.35−1.18------1.39inf--
118Narcissin------------1.04−1.71
119Astragalin1.25−1.22--------1.19−1.19--
124Silychristin--------------
126Eriodictyol----1.371.341.356.46--1.388.24--
13Isoorientin----1.322.84--1.83inf----
135Theaflavin 3,3′-digallate--------------
1366,2′-Dihydroxyflavone----1.221.051.05inf--1.39inf1.27−1.38
138Baimaside------------1.15−6.44
139Dihydrokaempferol------1.476.66--1.387.251.57−1.78
145Procyanidin B21.14.5----1.02−3.51------
147(-)-Epicatechin1.293.11----1.32−3.74------
150(-)-Catechin gallate1.356.76----1.2−6.36------
152Ononin1.13−inf------------
153Glycitin1.63inf1.771.941.721.781.121.62--1.39inf--
157Isomangiferin--1.631.49----------
160Apigenin 7-glucoside--------1.83inf--1.2inf
169Quercetin 3-O-(6″-galloyl)-β-D-galactopyranoside--------------
176Naringenin−7-glucoside------1.394.91.461.391.325.07--
177Phloretin1.08inf----1.347.641.07−1.111.39inf1.48−2.26
1832′-Hydroxygenistein------1.05inf--1.38inf--
184Afzelechin--------------
191Naringenin chalcone1.154.01--1.311.021.387.85--1.3811.81.41−1.42
195(-)-Gallocatechin1.152.18----1.13−2.11------
197Isorhamnetin 3-O-glucoside--------------
201Trilobatin------1.48inf1.261.41.39inf1.671.68
202Sieboldin--1.361.04----1.823.19--1.494.91
23Quercetin------1.05inf--1.39inf--
25Apigenin-7-glucuronide1.15inf1.98inf1.24inf--1.27−1.34----
28Scutellarin--------------
30Isoliquiritigenin----1.652.21.05inf--1.38inf1.61_2.77
42Nicotiflorin--------------
48Narirutin--1.89−5.141.27−1.89--1.724.391.3917.52--
50Luteolin----1.021.361.05inf--1.39inf--
55(-)-Catechin--1.811.431.571.33--------
56Taxifolin------1.48inf--1.38inf--
57Rutin--------------
58Isorhamnetin---−----------
65Phlorizin1.35−1.60----1.353.181.611.851.282.051.61.05
69Liquiritin------1.361.37--1.211.5--
80(-)-Gallocatechin gallate1.44.27----1.13−3.33------
84(-)-Epigallocatechin1.423.19----1.29−3.36------
90Cynaroside1.161.731.38inf----1.27−inf----
96Calycosin-7-O-β-D-glucoside1.64inf--1.74inf------1.23−inf
* Log2 FC ≥ 1 or ≤−1 and VIP ≥ 1 were considered to denote significant difference. “Inf” means infinity; “-“ in the table indicates insignificant difference.
Table 2. The individual contents (mg/kg DW) of 50 flavonoids at different stages in the fruit of WT and NOR-like1 tomatoes.
Table 2. The individual contents (mg/kg DW) of 50 flavonoids at different stages in the fruit of WT and NOR-like1 tomatoes.
IndexCompoundsWT-(mg/kg)WT-BR3
(mg/kg)		WT-BR9
(mg/kg)		NOR-like1-GR
(mg/kg)		NOR-like1-BR3
(mg/kg)		NOR-like1-BR9
(mg/kg)
1Astilbin0.021 ± 0.010.018 ± 00.013 ± 00.03 ± 0.010.024 ± 00.022 ± 0
2Miquelianin0.011 ± 00.018 ± 00.021 ± 0.010.013 ± 00.014 ± 00.06 ± 0.09
9Genistein-0.013 ± 00.018 ± 0-0.011 ± 0.010.03 ± 0.01
106Typhaneoside0.012 ± 0.010.09 ± 0.040.03 ± 0.030.018 ± 0.030.04 ± 0.070.04 ± 0.04
11Hydroxysafflor yellow A0.08 ± 0.060.17 ± 0.090.10 ± 0.030.08 ± 0.040.21 ± 0.020.22 ± 0.11
115Quercitrin0.015 ± 0.020.004 ± 00.002 ± 00.1 ± 0.090.004 ± 0.010.006 ± 0.01
117Engeletin-0.016 ± 0.010.008 ± 0.01-0.007 ± 0.010.008 ± 0.01
118Narcissin1.36 ± 0.311.76 ± 1.870.54 ± 0.120.8 ± 0.161.14 ± 0.280.77 ± 0.08
119Astragalin0.27 ± 0.060.12 ± 0.040.09 ± 0.070.12 ± 0.060.1 ± 0.020.23 ± 0.24
124Silychristin-0.013 ± 0.010.03 ± 0.010.003 ± 0.010.003 ± 0.010.017 ± 0.02
126Eriodictyol0.016 ± 04.89 ± 2.314.67 ± 2.540.08 ± 0.067.19 ± 5.1611.84 ± 3.07
13Isoorientin--0.005 ± 0.010.008 ± 0.01-0.04 ± 0.02
135Theaflavin 3,3′-digallate-0.05 ± 0.09----
1366,2′-Dihydroxyflavone-0.019 ± 00.007 ± 0-0.014 ± 0.010.015 ± 0
138Baimaside0.59 ± 10.5 ± 0.530.006 ± 00.12 ± 0.080.07 ± 0.10.023 ± 0.03
139Dihydrokaempferol0.004 ± 00.63 ± 0.10.19 ± 0.060.005 ± 00.51 ± 0.260.32 ± 0.15
145Procyanidin B20.016 ± 0.020.023 ± 0.020.04 ± 0.030.37 ± 0.170.03 ± 0.030.028 ± 0.05
147(-)-Epicatechin0.05 ± 0.030.04 ± 0.020.05 ± 0.020.45 ± 0.440.03 ± 0.010.04 ± 0.03
150(-)-Catechin gallate0.015 ± 0.010.01 ± 0.010.021 ± 0.031.57 ± 2.430.019 ± 0.020.017 ± 0.01
152Ononin0.04 ± 0.05-----
153Glycitin-0.05 ± 0.010.04 ± 00.06 ± 0.050.19 ± 0.10.13 ± 0.03
157Isomangiferin0.005 ± 00.002 ± 00.003 ± 00.005 ± 00.005 ± 00.006 ± 0
160Apigenin 7-glucoside--0.011 ± 0.01--0.013 ± 0.01
169Quercetin 3-O-(6″-galloyl)-β-D-galactopyranoside0.003 ± 0-----
176Naringenin-7-glucoside0.08 ± 0.082.58 ± 0.754.05 ± 1.640.05 ± 0.031.42 ± 0.913.72 ± 0.31
177Phloretin-1.80 ± 0.440.37 ± 0.170.005 ± 0.011.07 ± 0.680.50 ± 0.27
1832′-Hydroxygenistein-0.006 ± 00.003 ± 0-0.005 ± 00.005 ± 0
184Afzelechin0.04 ± 0.040.013 ± 0.020.019 ± 0.020.02 ± 0.04-0.009 ± 0.02
191Naringenin chalcone0.008 ± 028.46 ± 5.5710.67 ± 4.940.13 ± 0.1529.60 ± 22.1921.62 ± 5.88
195(-)-Gallocatechin0.04 ± 0.040.05 ± 0.010.08 ± 0.030.20 ± 0.190.05 ± 0.010.04 ± 0.01
197Isorhamnetin 3-O-glucoside0.004 ± 00.004 ± 0.010.002 ± 00.002 ± 00.001 ± 0-
201Trilobatin-0.008 ± 00.03 ± 0-0.007 ± 00.018 ± 0.01
202Sieboldin0.008 ± 0.010.007 ± 0.010.22 ± 0.050.017 ± 00.015 ± 0.010.14 ± 0.04
23Quercetin-0.13 ± 0.040.09 ± 0.05-0.09 ± 0.080.09 ± 0.03
25Apigenin-7-glucuronide---0.002 ± 00.007 ± 00.003 ± 0
28Scutellarin---0.019 ± 0.03-0.005 ± 0.01
30Isoliquiritigenin-0.005 ± 00.001 ± 0-0.006 ± 0.010.003 ± 0
42Nicotiflorin102.77 ± 14.5665.96 ± 25.7653.63 ± 13.2868.37 ± 45.1746.47 ± 13.7560.35 ± 23.33
48Narirutin-0.04 ± 0.010.09 ± 0.05-0.001 ± 00.025 ± 0.02
50Luteolin-0.009 ± 00.013 ± 0.01-0.015 ± 0.010.03 ± 0.01
55(-)-Catechin0.07 ± 0.030.06 ± 0.010.05 ± 00.34 ± 0.430.16 ± 0.060.12 ± 0.04
56Taxifolin-0.04 ± 0.010.06 ± 0.02-0.05 ± 0.040.09 ± 0.01
57Rutin175.26 ± 44.54168.46 ± 14.85147.79 ± 18.69101.18 ± 47.22125.43 ± 29.2145.01 ± 15.91
58Isorhamnetin-0.004 ± 0.01----
65Phlorizin0.012 ± 0.010.05 ± 0.010.11 ± 0.010.004 ± 00.04 ± 0.020.13 ± 0.02
69Liquiritin0.004 ± 00.011 ± 00.013 ± 0.010.004 ± 00.01 ± 00.013 ± 0
80(-)-Gallocatechin gallate0.09 ± 0.070.06 ± 0.040.11 ± 0.121.66 ± 1.570.17 ± 0.190.06 ± 0.02
84(-)-Epigallocatechin0.09 ± 0.060.06 ± 0.030.10 ± 0.110.85 ± 0.520.08 ± 0.060.04 ± 0.02
90Cynaroside0.05 ± 0.09--0.17 ± 0.030.06 ± 0.05-
96Calycosin-7-O-β-D-glucoside-0.009 ± 0.01-0.025 ± 0.010.025 ± 0.020.03 ± 0
total amount281.02 ± 28.02276.25 ± 22.95223.37 ± 33.78176.87 ± 66.48214.39 ± 39.51245.93 ± 35.24
* The result is indicated as the mean ± SD (n = 3); “-” in the table indicates not detected.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guan, D.; Zhao, Y.; Zhao, X.; Fu, D. Metabolomics Study of the Effect of Transcription Factor NOR-like1 on Flavonoids in Tomato at Different Stages of Maturity Using UPLC-MS/MS. Foods 2023, 12, 4445. https://doi.org/10.3390/foods12244445

AMA Style

Guan D, Zhao Y, Zhao X, Fu D. Metabolomics Study of the Effect of Transcription Factor NOR-like1 on Flavonoids in Tomato at Different Stages of Maturity Using UPLC-MS/MS. Foods. 2023; 12(24):4445. https://doi.org/10.3390/foods12244445

Chicago/Turabian Style

Guan, Di, Ying Zhao, Xiaodan Zhao, and Daqi Fu. 2023. "Metabolomics Study of the Effect of Transcription Factor NOR-like1 on Flavonoids in Tomato at Different Stages of Maturity Using UPLC-MS/MS" Foods 12, no. 24: 4445. https://doi.org/10.3390/foods12244445

APA Style

Guan, D., Zhao, Y., Zhao, X., & Fu, D. (2023). Metabolomics Study of the Effect of Transcription Factor NOR-like1 on Flavonoids in Tomato at Different Stages of Maturity Using UPLC-MS/MS. Foods, 12(24), 4445. https://doi.org/10.3390/foods12244445

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop