Comparison analysis of full-spectrum metabolomics revealed on the variation of potential metabolites of unscented, Chloranthus spicatus scented, and Osmanthus fragrans (Thunb.) Lour. scented Congou black teas

Introduction In recent years, scented black tea has attracted much attention due to its pleasant floral aroma and mellow flavor, but little research has been carried out on its flavor metabolic profile. Methods In this study, the flavor metabolic profiles of unscented, Chloranthus spicatus scented, and Osmanthus fragrans (Thunb.) Lour. scented Congou black teas were investigated using full-spectrum metabolomics analysis method, the first time that the flavor profiles of scented black tea were characterized in detail. Results and Discussion The results revealed that a total of 3,128 metabolites were detected in the three teas. Based on the criteria of variable importance in the project >1 and fold change ≥2 or ≤  0.5, 761 non-volatile metabolites and 509 volatile metabolites were filtered as differential metabolites. Many differential non-volatile metabolites belonged to flavonoids, phenolic acids, and terpenoids. Floral, fruity and herbaceous volatile metabolites were significantly up-regulated in Chloranthus spicatus scented Congou black tea while sweet and fruity volatile metabolites were significantly down-regulated in Osmanthus fragrans (Thunb.) Lour. scented Congou black tea. The results contribute to a better understanding of the scenting techniques on the flavor quality of scented black teas and provide some information on the flavor chemistry theory of scented black teas.


Introduction
Black tea, a fully fermented tea, is loved by consumers around the world for its health benefits (1)(2)(3)(4) and mellow taste.Black tea is widely produced in countries such as China, India, Sri Lanka, Kenya and Indonesia (5,6).It is popular and ranks first in the world in terms of production and consumption.The world black tea production accounted for around 65% of total tea production in 2017 (7).The Congou black tea (CK), a unique Chinese black tea variety, are with tightly striped shape, sweet and mellow taste, and floral and fruity aroma (8).The aroma of CK can be classified as sweet, floral or fruity, based on sensory characteristics (9).The general processes of producing CK are divided into two phases.The enzymatic phase includes the withering, rolling and fermentation processes, while the non-enzymatic phase include the drying process (10).The flavour quality is mainly formed in the enzymatic phase by the conversion of non-volatile metabolites, such as flavonoids, amino acids, fatty acids, and volatile metabolites derived from fatty acids, amino acids, terpenoids, carotenoids, and glycosidically bound volatiles (11).
In recent years, scented teas have been gaining attention for their pleasant floral aroma and mellow flavour.As a material with both ornamental and health benefits, such as protecting stomach and liver, dispersing phlegm, and antioxidant activity, Osmanthus fragrans (Thunb.)Lour. is ideal to be the raw material for the production of scented tea (12).Linalool, ionone and ocimene are aroma-contributing substances in it (13).Chloranthus spicatus is well-known for its elegance and is one of the "The Four Famous Chinese Camellias, " along with jasmine, Michelia alba DC. and Citrus aurantium.As one of the traditional Chinese scented teas, Chloranthus spicatus was also ideal for the production of scented tea.
The changes in various volatile metabolites during the scenting process are crucial to the formation of scented tea aroma.The scenting process decreased alcohols but increased ketones and esters in Osmanthus scented black tea (14).Nonyl aldehyde, γ-decalactone, 1-octene-3-ol, trans-2-hexenyl pentanoate, and phenethyl butyrate were key odorants contributing to its Osmanthus aroma.Currently, the changes of volatile metabolites during the scenting process of producing Chloranthus spicatus scented black tea is less studied.The changes of non-volatiles metabolites during the scenting process should have an effect on the formation of scented tea quality, but still largely remain unveiled.
Based on these facts, the present study characterized the flavor metabolites of CK, Osmanthus scented Congou black tea (OF), and Chloranthus spicatus scented Congou black tea (CH) by ultraperformance liquid chromatography tandem mass spectrometry (UPLC-MS/MS) and gas chromatography-mass spectrometry (GC-MS).The results aim to understand the effect of the scenting process on the metabolic profile of black tea and help broaden the ideas for the study of flavor quality of scented black tea.

Chemicals
Methanol, acetonitrile and hexane (all chromatography grade) for the tests were obtained from Merck (Darmstadt, Germany).Formic acid was purchased from Aladdin (Shanghai, China).Sodium chloride was produced from Sinopharm (Shanghai, China).Authentic standards used in the full-spectrum metabolomics analysis were obtained from MetWare (Wuhan, China).

Preparation of tea samples
The CK used in the experiment was purchased from Jingdezhen Jin Gui Yuan Agricultural Development Co (Jiangxi, China).The black tea was made according to the traditional Congou black teas process, which consisted mainly of indoor withering, rolling, fermentation and drying.The CH and OF were scented with CK as the raw material, and the specific scenting process was as follows (15): 1. Flower spreading: Pick Chloranthus spicatus and Osmanthus buds, ensure their integrity during the picking process, remove impurities such as flower tips, stalks and leaves and then stall them indoors.Indoor environment: ventilation, humidity of 65-75%, temperature of 20-24°C, thickness of 15 cm, spreading time of 30 min.2. Mixing and resting of tea flowers: 1/3 of the CK was used to produce CH and the other 1/3 was used to produce OF.For CH, the flowers are mixed with CK in a ratio of 3:100 (flowers: dried tea, w:w) for 100 d.For OF, the flowers are mixed with CK in a ratio of 3:10 (flowers: dried tea, w: w) for 2 days, with 1 renewal, the process lasts about 4 days in total.3. Drying: The scented teas were baked at 85°C (6CTH-30, Zhejiang Chunjiang Tea Machinery Co., Ltd.) for 412 min, cooled for 40 min, baked at 75°C for 90 min and cooled for 120 min before removing the petals and debris, resulting in scented Congou black teas.Three replicates were set for each treatment.The samples from three treatments were stored in a refrigerator at −80°C until assayed.The three samples were mixed and sampled three times for subsequent testing to ensure homogeneity and traceability.

UPLC-MS/MS based analysis of non-volatile metabolites 2.3.1. Sample pretreatment
The three samples were vacuum freeze-dried in a lyophilizer (Scientz-100F) before being ground (30 Hz, 1.5 min) into powder by a grinder (MM 400, Retsch).Afterwards, 50 mg of sample powder was accurately weighed and 1.2 mL of 70% methanolic aqueous internal standard extraction solution pre-chilled at −20°C was added.The sample was vortexed once every 30 min for 30 s, for a total of 6 times.Then centrifuged (12,000 rpm, 3 min), the supernatant was aspirated and filtered through 0.22 μm microporous membrane for determination.

UPLC-MS/MS conditions
An UPLC system equipped with a tandem mass spectrometer was used to analyze the non-volatile metabolites in the samples.(1) Liquid phase conditions: An Agilent SB-C18 1.8 μm, 2.1 mm * 100 mm column was equipped with a flow rate of 0.35 mL/min, a column temperature of 40°C, and an injection volume of 2 μL.The mobile phase A was 0.1% formic acid/ultrapure water and mobile phase B was 0.1% formic acid/ acetonitrile.The gradient elution method was as follows: mobile phase B was 5% at 0 min; increased linearly to 95% at 9 min, maintained for 1 min; decreased to 5% at 10-11.10 min and equilibrated at that ratio until 14 min.(2) Mass spectrometry conditions: electrospray ionization temperature (ESI) of 500°C, ion spray (IS) voltage of 5,500 V (positive ion mode) and − 4,500 V (negative ion mode); ion source gases I (GSI), gases II (GS II), and curtain gas (CUR) of 50, 60 and 25 psi, respectively; collision-induced ionization parameters of high.The multiple reaction monitoring (MRM) mode was used for triple quadrupole (QQQ) scanning, and the collision gas (nitrogen) was medium.The declustering potential (DP) and collision energy (CE) of individual MRM ion pairs were evaluated by further DP and CE optimization.A specific set of MRM ion pairs was monitored in each period depending on the metabolites eluted during the period (16).

Identification of non-volatile metabolites
Based on the MWDB (metware database), the substance characterization was performed based on the secondary spectrum information of the substance, and the analysis removed the isotopic signals, the duplicate signals containing K + , Na + , NH 4 + , and the duplicate signals of fragment ions that were themselves other larger molecular weight substances.After obtaining the metabolite spectral analysis data of different samples, the peaks of all substances were integrated by peak area using MultiQuant software, and the mass spectra of the same metabolite in different samples among them were corrected for the integration of the peaks.

GC-MS based analysis of volatile metabolites 2.4.1. Sample pretreatment
The sample was removed from the refrigerator at −80°C and firstly ground by liquid nitrogen, followed by vortexing to make a homogeneous mixture.The 500 mg of the sample was accurately weighed in a headspace vial, subsequently 2 mL of saturated sodium chloride solution were added, the vial was sealed, and then the sample was extracted by fully automated headspace solid phase microextraction (HS-SPME) for GC-MS analysis.

GC-MS conditions
The headspace vial containing the tea samples was shaken and equilibrated at 60°C for 5 min, followed by insertion of a 120 μm divinylbenzene/carboxen/polydimethylsiloxane (DVB/CWR/PDMS) extraction head into the vial for 15 min (stirring while extracting was used in the extraction process with a stirring rate of 1,000 rpm, thus stability and homogeneity of the extraction process and reducing human error), followed by desorption at the inlet end at 250°C for 5 min for subsequent GC-MS separation and identification.It should be noted that the extraction head should be heated at the inlet end at 250°C for 5 min before sampling.(1) Chromatographic conditions: DB-5MS capillary column (30 m × 0.25 mm × 0.25 μm, Agilent J&W Scientific, Folsom, CA, United States) supplemented with high purity helium (purity not less than 99.999%) as carrier gas.The program temperature was set as follows: 40°C for 3.5 min, then 10°C/min to 100°C, then 7°C/min to 180°C, and finally 25°C/min to 280°C for 5 min.(2) Mass spectrometry conditions: electron bombardment ion source (EI); ion source, quadrupole and mass spectrometry interface temperatures of 230°C, 150°C and 280°C, respectively; electron energy of 70 eV, scanning mode of selected ion detection mode (SIM), qualitative and quantitative ion precision scanning (GB 23200.8-2016).

Identification of volatile metabolites
Based on multi-species, literatures, authentic standards, and retention index (RI), the database was established independently, containing the identified retention time (RT) as well as qualitative and quantitative ions for the selected ion detection mode for precise scanning, with one quantitative ion and two to three qualitative ions selected for each compound, respectively.All ions to be detected in each group were detected separately in the order of of their peaks, in separate time periods.If the retention time of the detected peaks was consistent with the authentic standard and the selected ions all appeared in the mass spectra of the samples after background deduction, the substance would be ascertained.The odor descriptions for the qualitative aroma substances were taken from websites such as http://www.thegoodscentscompany.com, http://perflavory.com/,http://www.odour.org.uk, http://foodflavorlab.cn.The quantitative ions were selected by MassHunter software for the integration and calibration of the peaks in order to enhance the accuracy of quantification.

Statistical analysis
All samples were set up with three biological replicates and statistical analysis was performed using IBM SPSS statistics (version 25.0, SPSS Inc., Chicago, IL).The principal component analysis (PCA), partial least squares-discriminant analysis (PLS-DA), clustering heat map, K-Means plot, and key metabolite kyoto encyclopedia of genes and genomes (KEGG) enrichment pathway maps were done by R software. 1 The identified metabolites were first annotated through the KEGG compound database, 2 and subsequently the annotated completed metabolites were mapped to the KEGG pathway database. 3 The pathways mapped to metabolites with significant regulation were subsequently sent to metabolite set enrichment analysis (MSEA), whose significance was determined by p-values of hypergeometric tests.

Overall profile of metabolites in the three black teas
Metabolomics allows the detection and screening of significant differential metabolites from different biological samples, which can be used as a basis to further elucidate the metabolic processes and mechanisms of change in biological samples.In this study, a fullspectrum metabolomics analysis based on dual platforms of UPLC-MS/MS and GC-MS was carried out.A total of 3,128 metabolites, including 2,396 non-volatile metabolites and 732 volatile metabolites, were detected in three black tea samples.
The heatmap distribution in Figure 1B indicated that the metabolite contents in the three black tea samples were significantly different.To further validate the above results, PCA, as an unsupervised pattern recognition method for multidimensional data statistical analysis, was introduced to analyze the differences of overall metabolites (Figure 1C), non-volatile metabolites (Figure 1D), and volatile metabolites (Figure 1E) among the three black tea samples.The results revealed that the three samples achieved good separation in each PCA model, and the sum of the first and second principal components of the above three data sets explained more than 50% of the whole data set, reaching 63.81, 50.41 and 79.91%, respectively, indicating that the three black tea metabolites reached significant

Important differential metabolites among the three black teas
1,270 differential metabolites were initially identified based on the screening criteria of variable importance in the project (VIP)>1 and fold change (FC) ≥2 or ≤ 0.5, among which 761 non-volatile metabolites (Supplementary Table S1) and 509 volatile metabolites (Supplementary Table S2) were identified.
The Venn plot indicated the relationships of differential metabolites between different black tea samples.The results showed that there were 59 differential metabolites between the three comparison groups (CK vs. OF, CK vs. CH, OF vs. CH) and 82, 50, 96 differential metabolites were unique to each group (Figure 2A).To further verify whether the differential metabolites among the three black tea samples reached significant differences, an unsupervised model statistical method, i.e., PCA (Figure 2B), and a supervised model statistical method, i.e., PLS-DA (Figure 2C), were used simultaneously to screen for differential metabolites.It was found that both models could achieve good separation.In 200 permutations of validation of the PLS-DA model (Figure 2D), R2X = 0.733, R2Y = 1 (p < 0.005), Q2 = 0.979 (p < 0.005), indicating a good model fit.
The Volcano plot was mainly used to show the relative abundances of metabolites between two samples.Each dot in the plot indicated a differential metabolite.Green dots indicated down-regulated substances, red dots indicated up-regulated substances, and gray indicated substances with insignificant differences.The results showed that compared to CK, 197 substances were up-regulated and 85 substances were down-regulated in CH (Figure 2E), while 92 substances were up-regulated and 193 substances were downregulated in OF (Figure 2F), which might be related to the different flowers used in the scenting process.
The K-means plot focused on classifying the trends of relative abundances of important differential metabolites under different comparison groups to facilitate subsequent targeted analysis.According to Figure 2G, three different trends in the levels of these important differential metabolites were detected in CK, CH and OF.Among them, 342 metabolites were lower in OF than in CK and 371 metabolites were higher in CH than in CK.It was obvious that more substances were up-regulated in CH and down-regulated in OF after CK was scented, suggesting that the type of flower had a great influence on the metabolites of CK.
Subsequently, the top 20 KEGG metabolic pathways in the three samples were screened based on differential abundance (DA) scores and p-values sizes (Figures 2H,I).The sesquiterpenoid and triterpenoid biosynthesis, limonene and pinene degradation, and tyrosine metabolism reached significant levels in the metabolite expression of CH (compared to CK), with the first two pathways being significantly up-regulated and the last one remaining basically stable.Sesquiterpenoids and triterpenoids were generally biosynthesized via the cytoplasmic mevalonate (MVA) pathway in plants (18)(19)(20).Triterpenoids usually had a variety of stereoisomers and positional isomers, but the main compounds were the madecassic acid isomer and asiatic acid isomer (21).In OF, the expression of metabolic pathways such as terpenoid backbone biosynthesis, flavone and flavonol biosynthesis, and tyrosine metabolism were significantly down-regulated overall.Terpenoids were derived from the same C5 isoprene backbone produced from isopentenyl pyrophosphates (IPP) and dimethyl allyl pyrophosphates (DMAPP) and were used in the synthesis of terpene backbones by rearrangement and cyclisation reactions (22), where linalool, nerol and alpha-terpineol could be produced by further reactions in the biosynthesis of terpene backbones (23).The results revealed differences in flavor quality between different scented black teas (CH and OF) from a metabolomics perspective.However, the most critical contributors to this difference were still unknown, both in terms of non-volatile and volatile substances.

Analysis of key differential non-volatile metabolites among the three black teas 3.3.1. Analysis of key differential categories
Based on the percentage of metabolite categories, all non-volatile metabolites were further screened and a total of 173 metabolites were selected (Table 1).In terms of categories, the metabolites in the three black teas mainly consisted of amino acids and derivatives, phenolic acids, flavonoids, lipids, accounting for 7.69, 20.12, 52.66 and 19.53%, respectively.In terms of peak area, flavonoids accounted for a relatively large proportion of all key metabolites, reaching 96.30% in CK, and compared with CK, the percentage of peak area of flavonoids increased by 0.39 and 6.72% in CH and OF, respectively.The variation of flavonoids was determined by the thermal stability during the drying stage (17).Compared with CK, the peak area of lipids in CH was reduced by 0.43%.The lipids were important precursors of volatile metabolites in tea leaves (24).It was reported that lipids were degraded by enzymatic hydrolysis or oxidation of glycolipids and phospholipids to produce flavor volatiles during black tea processing (25), which was in accordance with our experimental results.The total peak area of amino acids and derivatives in OF was reduced by 5.18% compared to CK (Supplementary Table S1), which indicated that the Maillard reaction, deamidation polymerization and caramelization might occur during the scenting process (26).
Overall, there was a downward trend in the non-volatile metabolites of black tea after scenting.Compared to CK, phenolic acids were significantly up-regulated in CH (Log2FC = 31.66)whilst flavonoids were significantly down-regulated (Log2FC = −39.00).Amino acids and derivatives were significantly up-regulated in OF (Log2FC = 17.60) whilst flavonoids were significantly down-regulated (Log2FC = −30.43).

Analysis of key differential compounds
Based on the values of magnitude of the difference multiplier (Log2FC) in each comparison group (CK vs. CH, CK vs. OF), the substances ranked the top 20 of the key non-volatile differential metabolites were prioritized for further discussion.As illustrated in Figures 3A,B, the Log2FC of the selected metabolites ranged from −22.58 to 18.83 in CK vs. CH and − 14.62 to 18.03 in CK vs. OF.Based on |Log2FC| > 15, 7 and 4 substances were screened as key differential compounds in CH (Figure 3D) and OF (Figure 3C), respectively.
The results revealed that the 11 substances were divided into two main groups, which were flavonoids and phenolic acids.Combined with the relevant information in Table 1, it was concluded that compared to CK, forsythoside H, alpha-asarone,        ) Vahl, was a type of forsythoside and might possess anti-inflammatory activity.Alpha-asarone was a polyphenol with anti-bacterial and anti-inflammatory activities.

Analysis of key differential volatile metabolites among the three black teas 3.4.1. Analysis of key differential categories
Based on the percentage of metabolite categories, all volatile metabolites were further screened and 110 metabolites with aroma contribution were selected, including 11 alcohols, 9 ketones, 29 esters, 31 terpenoids, 14 heterocyclic compounds, 4 aldehydes, 5 aromatics, 3 acids, 3 nitrogen compounds, and 1 phenol (Table 2).It was consistent with a previous study which revealed that terpenoids were the most important aroma components in CK (28).The sum of relative abundances of esters, terpenoids, heterocyclic compounds, alcohols, ketones and aromatics was 97.08% in CK.Compared to CK, the relative abundances of heterocyclic compounds and alcohols in CH decreased by 19.62 and 11.49%, respectively, while aromatics increased by 21.23%.In OF, the relative abundances of heterocyclic compounds and esters increased by 17.21 and 13.59%, respectively, while the relative abundances of terpenoids and ketones showed a decreasing trend, reaching 14.87 and 10.31%, respectively (Supplementary Table S2).A previous study found a significant increase in both ketones and esters during the scenting of Osmanthus black tea (14), which was not consistent with the results of this study and may be influenced by the variety of tea leaves, tea processing parameters and scenting process.

Analysis of key differential compounds
Based on the magnitude of Log2FC in each comparison group (CK vs. CH, CK vs. OF), the top 20 substances in the key volatile differential metabolites were selected preferentially and their FCs were shown in Figures 4A,B.It was found that the selected metabolites were significantly up-regulated in CH compared to CK, while those in OF were significantly down-regulated.Based on |Log2FC| > 6, six and five substances were screened as key differential compounds in CH (Figure 4C) and OF (Figure 4D), respectively.
Among them, the terpenoids and ester accounted for a greater number, with six and three varieties, respectively.Acids and heterocyclic compounds each accounted for one.Combined with the relevant information in Table 2, the relative abundances of cis-beta-farnesene, bornyl butyrate, gamma-decalactone, d-nerolidol, gamma-selinene, and trans-alpha-bergamotene were significantly increased in CH (compared to CK) and FC values reached 51533.83,4183.31,400.92, 116.77, 64.68, 64.57, respectively.The cis-beta-farnesene and trans-alpha-bergamotene were significantly associated with floral aroma (28, 29).Transalpha-bergamotene was a citrus odorant which was considered to be the most important sesquiterpene in the characteristic odour of Polygonum minus (30).Bornyl butyrate, which presented an herbaceous aroma (31), usually presented in the Aloysia triphylla and Aloysia citrodora. 4Gamma-decalactone, which delivered a pleasant peachy aroma, was regarded as one of the characteristic components contributing to the Osmanthus aroma of Osmanthus black tea (14).Nerolidol was considered to be one of the skeleton aroma components of CK (32).The decrease of nerolidol after scenting might be related to the drying process, as it was previously reported that there was a significant positive correlation between the content of nerolidol and the degree of roasting (33).Gamma-selinene, as a sesquiterpene with floral aroma, was detected in the essential oils from Chloranthus spicatus flowers (34).
Nevertheless, in OF, sabinene, neryl isobutyrate, 2-methyl-2pentenoic acid, beta-ionone epoxide, and coumarin were significantly down-regulated (FC ≤ 0.01).Sabinene, a representative of the bicyclic monoterpenes, conveyed a typical citrus note (35).Beta-ionone epoxide (sweet berry aroma), an apo carotenoid monoterpenoid, was found decreased in summer green tea (Camellia sinensis L.) after yellowing (36) and steamed green tea after baking (37), which was consistent with the results of this study.Neryl isobutyrate rendered a sweet rose scent and 2-methyl-2-pentenoic acid offered a fruity note (31).Chen (38) and Su et al. (32) both revealed that coumarin was one of the key substances in the formation of the typical aroma of special grade Keemun black tea.The coumarin was described as a sweet grass and cherry blossom-like aroma, which was formed by intramolecular esterification of hydroxycinnamic acid following hydrolysis of 2-coumaric acid primeveroside (18).The reduction of coumarin in OF might be related to the drying process (39).
In conclusion, terpenes, including monoterpenes and sesquiterpenes, were the main differential category.The increase of cis-beta-farnesene and bornyl butyrate might contribute to an increased in floral, fruity and herbaceous aromas in CH.The decrease of sabinene and neryl isobutyrate might be related to the insignificant sweet fruity aroma in OF.The result was generally consistent with the results of KEGG metabolic pathways analysis (Figures 2H,I).As such, potential differential metabolites between the CK vs. CH and CK vs. OF comparator groups were characterized in scented teas with CK as the tea raw material.

Conclusion
The full-spectrum metabolomics (UPLC-MS/MS and GC-MS) adopted in this study permits a more comprehensive and objective characterization of the differences in non-volatile and volatile flavor metabolites among unscented and scented Congou  black teas.The results indicated that the scenting process and the types of flowers had a significant effect on the flavor quality of Congou black tea.Flavonoids, phenolic acids and terpenoids were screened as key differential categories by PCA and PLS-DA, with most significant changes in terpenoids in different scented teas.Floral, fruity and herbaceous volatiles were significantly up-regulated in CH, whilst the sweet fruity volatiles were significantly down-regulated.This study is the first to characterize the flavor profiles of CH and OF in some detail, the results enrich the research on flavor metabolisms of scented black tea and has certain implications for the processing and quality control of scented black tea.Future research needs to broaden the range of other types of flowers and to extend this study to the other five major tea categories, so as to find the most suitable flowers for each tea category to better guide the production of scented teas.

FIGURE 1
FIGURE 1Analysis of all metabolites in three types of black tea (CK, CH and OF).(A) Ring diagram of total metabolite class composition; (B) Heat map of total metabolite content clustering (processed by Z-score); (C) PCA analysis of total metabolites; (D) PCA analysis of non-volatile metabolites; (E) PCA analysis of volatile metabolites.CK, Congou black tea; OF, Osmanthus black tea; CH, Chloranthus spicatus black tea.Same below.

FIGURE 2
FIGURE 2 Analysis of important differential metabolites in three types of black tea (CK, CH and OF).(A) Venn plot; (B) PCA plot; (C) PLS-DA plot; (D) PLS-DA 200-permutation plot; (E,F) Volcano plot (E.CK vs. CH; F. CK vs. OF); (G) K-means plot; (H,I) KEGG differential abundance score plot (H.CK vs. CH; I. CK vs. OF; where the color and size of the dots indicated the size of the p-values and number of metabolites, respectively; DA score indicated the overall change of all metabolites in the metabolic pathway, 1: up-regulation, −1: down-regulation, and the length of the line segment indicated the DA score).

a
Molecular weight of the parent ion of the substance after ionization by the electrospray ion source.b Characteristic fragment ions.c Ionization mode (M + H is positively charged, M-H is negatively charged).d An * indicates the possible presence of isomers.eMetabolite identification criteriawhere 1 indicates a score of 0.7 or higher for the sample substance secondary mass spectrum (all fragment daughter ions of the substance), RT and database substance match; 2indicates a score of 0.5-0.7 for the sample substance secondary mass spectrum (all fragment daughter ions of the substance), RT and database substance match; 3: sample substance Q1, Q3, RT, DP, CE and database substance match.f CH_FC and OF_FC indicate the fold change of CK vs FC and CK vs OF metabolites, respectively.g Fold change were logarithmic with a base of 2. h Metabolite up-and down-regulation type.

TABLE 1
Key differential non-volatile metabolites changes among the three black teas.

TABLE 2
Key differential volatile metabolites changes among the three black teas.
a Retention index of metabolites on a non-polar column.b Quantitative ions.c Qualitative ions.d e CH_FC and f Fold change were logarithmic with a base of 2. g Metabolite up-and down-regulation type.