Integration of lipidomics and metabolomics approaches for the discrimination of harvest time of green tea in spring season by using UPLC-Triple-TOF/MS coupled with chemometrics

The production season is one of the judgment standards of the green tea quality and spring tea is generally considered of higher quality. Moreover, early spring tea is usually more precious and sells for a higher price. Therefore, a multifaceted strategy that integrates lipidomics and metabolomics, based on UPLC-Triple-TOF/MS coupled with chemometrics, was developed to discriminate early spring green tea (ET) and late spring green tea (LT). Twenty-six lipids and forty-five metabolites were identified as characteristic components. As for characteristic lipids, most of glycerophospholipids and acylglycerolipids have higher contents in ET. By contrast, glycoglycerolipids, sphingolipids and hydroxypheophytin a were shown higher levels in LT samples. Most of the differential metabolites identified were more abundant in ET samples. LT samples have much higher catechin, procyanidin B2, and 3',8-dimethoxyapigenin 7-glucoside contents. Based on the integration of differential lipids and metabolites, the reconstructed orthogonal partial least squares discriminant analysis (OPLS-DA) model displayed 100% correct classification rates for harvest time discrimination of green tea samples. These results demonstrated that the integration of lipidomics and metabolomics approaches is a promising method for the discrimination of tea quality.


. Introduction
Tea is considered as the most widely consumed beverage except water in the world. Green tea is the most prevalent, particularly in East Asian countries, due to its exclusive satisfactory sensory qualities and health benefits (Ratnani and Malik, 2022). The price of green tea is related to its quality, good quality usually corresponds to a high price (Li et al., 2021b). Some profitseeking merchants label low-quality tea as high-quality tea, which is harmful to consumer protection and the development of the tea market. Therefore, tea quality evaluation methods are urgently needed.
The quality of tea is influenced by many factors, including variety, climate, soil, harvest season, manufacturing process, and storage (Yin et al., 2022;Zhao et al., 2022). In many cases, the production season is the most important judgment standard of the green tea quality, which has become a decisive factor in ascertaining its price. It was reported that spring tea usually contains higher levels of free amino acids and moderate levels of catechins, which contribute to .
umami and mellow flavor of tea (Xu et al., 2012;Liu et al., 2015). On the contrary, higher levels of procyanidin B2 and gallotannins enhance the bitterness and astringency of teas (Jiang et al., 2022). Therefore, it's meaningful to understand the variation of chemical constituents of green tea related to season harvest period for the development of tea quality evaluation methods. Several studies have investigated the differences in chemical constituents of tea harvested in different seasons. These studies suggest that chemical components such as amino acids and catechins exhibited sharp seasonal fluctuations (Dai et al., 2015;Ma et al., 2022). The contents of amino acids decreased whereas carbohydrates, flavonoids and their glycosides, tricarboxylic acid cycle and photorespiration pathways were strongly reinforced in the late spring season (Liu et al., 2016;Ma et al., 2022). In our previous study, we also found that there are higher contents of gallocatechin and procyanidin B1 in summer tea (Peng et al., 2021). However, most of these studies were focused on the measurements of hydrosoluble metabolites, and there are few studies on hydrophobic compounds. In fact, some hydrophobic compounds also have an important influence on the flavor quality of tea. For example, lipid-soluble chlorophylls are major pigments in fresh tea leaves that constitute an essential basis of the unique color of green tea (Yu et al., 2019). Furthermore, the oxidation and degradation products of unsaturated fatty acids produce a number of fragrant volatiles, which are responsible for the fresh odor of tea as well (Ho et al., 2015). Herein, it is meaningful to study hydrosoluble and hydrophobic compounds simultaneously so as to gain a comprehensive understanding of the differences in metabolites of green tea with different qualities.
Lipidomics and metabolomics are powerful tools that enable comprehensive investigation of the lipids and metabolites possible in tea samples. Li et al. (2017Li et al. ( , 2021a reveals dynamic changes in tea lipids and their metabolic pathways during the manufacturing process of black tea and green tea using LC-MS combined with chemometric methods. However, the ability of single lipidomics or metabolomics to discover crucial compounds responsible for the quality evaluation and authenticity identification is limited. Metabolomics mainly focuses on polar metabolites. Lipidomics strives to comprehensively identify and quantify all kinds of lipid molecular species (Wenk, 2010). Lipids need to be treated and analyzed separately (i.e., requiring different solvent systems) from small-molecule metabolites due to their hydrophobic nature. Therefore, the integration of lipidomics and metabolomics can greatly expand the scope of composition analysis (Liu and Xu, 2018;Chen et al., 2021;Jia et al., 2021), which was a powerful strategy covering both polar and non-polar compounds. The current study aims to roundly expound the differentiation of chemical compositions between the early spring green tea (ET) and late spring green tea (LT) in combination with UPLC-Triple-TOF/MS-based lipidomics and metabolomics. Through the integration of lipidomics and metabolomics, more differential compounds between ET and LT may be screened based on multivariate statistical analysis methods, so as to be used for the picking time discrimination of spring tea.
. Materials and methods . . Reagents and materials HPLC-grade methanol (MeOH) was supplied by J. T. Baker (Danville, PA), and HPLC-grade methyl tert-butyl ether (MTBE) was purchased from Aladdin Biochemical Technology Co., Ltd (Shanghai, China). MS-grade of Acetonitrile (ACN) and HPLCgrade isopropanol (IPA) was obtained from Thermo Fisher Scientific Co., Ltd (Shanghai, China). MS-grade of Formic acid (FA) and ammonium acetate were available from ACS chemistry Company (Shanghai, China). Distilled water was obtained from Watsons (Hong Kong, China). Procyanidin B2, Myricitrin, Diosmetin-7-O-β-Dglucopyranoside and 4-OH-coumarin were purchased from Shanghai Yuanye Bio-Technology Co., Ltd (Shanghai, China). Theanine was obtained from Shanghai Anpu Experimental Technology Co., Ltd (Shanghai, China). The purity of all standards was >98%. Accurately weigh 1 mg of each standard and dissolve to 10 mL of methanol. A series of standard stock solutions at a concentration of 0.1 mg/mL were obtained. Then, a mixed standard solution was prepared by mixing 100 µL of each stock solution. All the solutions were stored in a brown volumetric flask at 4 • C before analysis.

. . Green tea samples and preparation
A total of 20 Laoshan green tea samples were collected from Qingdao Laoshan Wanggezhuang Tea Factory (Shandong Province, China). Laoshan green tea is the northernmost green tea in China, which is usually picked from mid-to-late April to late September. Among them, spring tea is picked from April to June. In this work, Grain Rain, the Chinese solar term, is the time point for dividing early spring tea and late spring tea. Ten early spring green tea (ET) were picked on April 16 th , and the late spring green teas (LT) were picked on May 24 th (LT) in 2020. These tea samples are taking one bud with two leaves and belong to the same cultivars. The fresh leaves were processed into tea products through spreading out, deenzyming, rolling and drying in accordance with flat tea process technology. Samples were ground into powder using a grinding miller and passed through a 60-mesh sieve and then stored at −18 • C. Each sample was in triplicate. The lipids and metabolites of tea were extracted on the basis of the previously reported procedure with minor modifications (Li et al., 2017;Peng et al., 2021). Briefly, 20 mg tea powder was accurately weighed in a 1 mL centrifuge tube. 300 µL of MeOH and 1 mL of MTBE were added to the tube. After continuous vortex mixing for 1 h, 300 µL distilled water was added to form a twophase system. After standing for 1 min at room temperature, the twophase system was separated by centrifugation for 10 min at 12,000 rmp. The supernatants were lipid metabolites. Then, 200 µL of the upper phase was transferred to a clean centrifuge tube, dried naturally at room temperature and stored at −80 • C until analyzed. Each dried sample was dissolved in 200 µL of ACN/IPA/H 2 O (65:30:5, v/v/v) and filtered through 0.22 µm PTME hydrophobic membrane. Finally, they were transferred into sample vials for analysis. Quality control (QC) samples were prepared by mixing 50 µL of lipid extracts of tea samples to evaluate the lipidomics analysis process.
For the extraction of characteristic metabolites of tea samples, 20 mg of tea powder was accurately weighed into a centrifuge tube. Then, 1 mL of 70% methanol was added. The mixtures were extracted with ultrasonic for 10 min at 60 • C and centrifuged at 12,000 rpm for 10 min. After cooling to room temperature, the supernatants were diluted 5-fold with distilled water and filtered through a 0.22 µm PTME hydrophilic membrane. All stock solutions were stored at 4 • C prior to UPLC-Triple-TOF/MS . /fsufs. .
analysis. QC samples were prepared in the same way as in lipidomics analysis. The QC samples were inserted at every ten tea samples during analysis, and one QC was inserted before and after injection, respectively.

. . UPLC-triple-TOF/MS analysis
The analysis of tea samples was carried out on an ultraperformance liquid chromatography (UPLC)-Triple-time-of-flight mass spectrometer system (Triple TOF TM 5600 + , AB SCIEX, USA) in both positive and negative modes. The LC separation was performed on an ExionLC TM system consisting of a binary high pressure mixing gradient pump with degasser, a thermostated auto-sampler, a controller, a photodiode array (PDA) detector, a reservoir tray, and a column oven. The Triple TOF TM 5600 + mass spectrometer was equipped with a DuoSpray TM source.
For lipidomics analysis, a Kinetex XB-C18 column (2.6 µm, 100 × 2.1 mm, Phenomenex, USA) was used for chromatographic separations with a column temperature of 50 • C. Gradient elution was carried out using mobile phase A ACN/H 2 O (60:40, v/v) and B IPA/ACN (90:10, v/v), both containing 5 mM ammonium acetate, at a flow rate of 0.26 mL/min. The optimal gradient elution program was set as follows: 0-2 min, 32% B; 2-4 min, 32-60% B; 4-13 min, 60-97% B; 13-17 min. 97% B; 17-17.1 min, 97-32% B; 17.1-20 min, 32% B, a total elapsed time was 20 min. The samples were kept at 4 • C in the auto-sampler and 5 µL was injected for analysis. Full scan (TOF-MS) and data-independent acquisition (DIA) modes were used simultaneously to acquire MS and MS/MS spectra. The TOF-MS and MS/MS analysis were acquired in positive mode with the scan ranges of 200-1,200 and 100-1,200 Da for the precursor ions and product ions, respectively. MS acquisition parameters were set as follows: ionization source temperature, 500 • C; ionspray voltage floating, 5.5 kV; gas 1 pressure, 55 psi, and gas 2 pressure, 55 psi; curtain gas pressure, 35 psi; declustering potential, 70 V; collision energy was of 45 V with the spread of ±15 V. An automatic calibration was performed every 10 samples by the calibration delivery system (CDS, APCI calibration solution) to ensure mass accuracy during the measurement.
For metabolomics analysis, chromatographic separations were performed on a Kinetex F5 column (2.6 µm, 100 × 2.1 mm, Phenomenex, USA) with a column temperature of 30 • C. Gradient elution was carried out using mobile phase (A) water with 0.1% (v/v) formic acid and (B) ACN at a flow rate of 0.3 mL/min. The gradient elution program was operated as follows: 0-0.5 min, 5% B; 0.5-3 min, 5-40% B; 3-9 min, 40-90% B; 9-12 min, 95% B; 12-12.1 min, 95-5% B; 12.1-15 min, 5% B. The injection volume of the sample was set at 5 µL. The temperature of the auto-sampler was kept at 10 • C. Full scan (TOF-MS) and data-independent acquisition (DIA) modes were used simultaneously to acquire MS and MS/MS spectra. The TOF-MS and MS/MS analysis were acquired in both positive and negative ionization modes with the scan ranges of 100-1,000 and 50-1,000 Da for the precursor ions and product ions, respectively. The MS parameters in positive and negative ionization modes were set as follows: ionization source temperature, 550 • C; capillary voltage, +5.5 kV (ESI+) and −4.5 kV (ESI-); gas 1 pressure, 55 psi, and gas 2 pressure, 55 psi; curtain gas pressure, 25 psi; declustering potential, 70 V; collision energy was of 35 V with the spread of ±15 V. As in lipid analysis, automatic calibration was also performed every 10 samples to ensure mass accuracy during the measurement.

. . Data processing and chemometric analysis
The raw LC-MS data were analyzed by the MarkerView TM 1.3 software (AB Sciex LLC., Framingham, MA, USA). The process parameters were shown in Support Information. The retention time range for peak finding in lipidomics analysis and metabolomics analysis were 0.5 and 17, 0.5 and 13 min, respectively. The subtraction multiple factor was set to 1.5 and the noise threshold was set to 100. The retention time tolerance and mass tolerance were set to be 0.5 min and 10 ppm for peak alignment, respectively. Two conditions were used to filter the extraction peaks: one is to remove peaks in <30 samples and the other is to set the maximum number of peaks to be 5,000 and 4,000 for lipidomics and metabolomics analysis, respectively. The data was processed using total area sum normalization before exporting.
A table of peak values was exported for principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA) by using SIMCA-P software (MKS Umetrics AB, Umeå, Sweden). The quality of models was defined by the parameters of R 2 X cum , R 2 Y cum , and Q 2 cum , in which R 2 X cum and R 2 Y cum represent the fraction of the variance of the X matrix and Y matrix, respectively, while Q 2 cum indicates the predictive accuracy of the model (Szeto et al., 2010). The model is considered reliable when three indicators (R 2 X cum , R 2 Y cum , and Q 2 cum ) higher than 0.5 and close to 1 indicate an excellent model. A permutation test with 200 permutations was conducted to assess the validity of the model (Mahadevan et al., 2008). HCA heatmap analysis was performed with R program (version 3.6.2) using Manhattan distance measurement method for differential lipids and Euclidean distance measurement method for differential metabolites, respectively.
/fsufs. . . Results and discussion . . Lipidomics and metabolomics analysis based on UPLC-triple-TOF/MS As shown in Figure 1, the TICs of lipids and metabolites of green tea picked in April and May were very similar. Further statistical analysis is necessary. After data preprocessing, a total of 1,329 ions were extracted in lipidomics analysis. In metabolomics analysis, a total of 2,737 and 2,658 ions were acquired in positive and negative modes, respectively. With the combination of lipidomics and metabolomics, plenty of chemical information was obtained for further statistical analysis.

. . Lipidomics character of early and late spring Laoshan tea
As displayed in the PCA score plot (Figure 2A), the QC samples were gathered tightly, indicating the good stability and precision of the lipid profile acquisition method. ET and LT samples can be separated in view of the lipid profiles, demonstrating the lipids differences between ET and LT. Subsequently, a supervised method OPLS-DA was applied to study the lipidomics data. As shown in Figure 2B, excellent classification result was achieved for ET and LT samples. The R 2 X cum , R 2 Y cum , and Q 2 cum of the OPLS-DA model based on lipidomics data were 0.595, 0.977, and 0.912 (Table 1), respectively, indicating the wonderful predictive ability of the model. The results of 200 times permutation test ( Figure 3A) also exhibited that the model is reliable and not over-fitted (Q 2 < 0.05, R 2 values < original values). Based on the OPLS-DA model and univariate analysis (p-value and fold change), lipids with VIP > 1, p-values < 0.05 and fold change >2 or <0.5 were regarded as differential lipids. As a result, 26 lipid ions were screened out.
The 26 lipid ions were identified as follows: these lipids were preliminarily discriminated by MS-DIAL and MS-FINDER software based on MS and MS 2 information. The potential name of the compound was predicted by MS-DIAL software, and MS-FINDER provides the formula of the compound. Then, the chemical structure of the compound was searched in databases such as LipidMaps, PubChem, MassBank, HMDB, MassBank of North America and ChemSpider and downloaded as a mol format file. The file was used to generate theoretical simulated MS/MS fragments by using the PeakView software. The theoretical fragments are then matched with the fragment pane of the corresponding experimental MS/MS fragments. Match scores >70 are acceptable. At last, 25 lipid molecules were identified and divided into 5 lipid species consisting of 1 chlorophyll derivative, 7 glycoglycerolipids, 1 sphingolipid, 12 glycerophospholipids, 4 acylglycerolipids. Table 2 summarized the detailed information of these identified chemical components.

. . Metabolomics character of early and late spring Laoshan tea
Good aggregation of QC samples in the PCA score plot ( Figures 4A, B) indicated that the metabolomics acquisition method was feasible. Although a certain degree of separation of ET and LT samples can be reflected, the overlap between partial samples makes the PCA models difficult to use for satisfactory classification. The supervised method was further used for classification and discrimination. The OPLS-DA score plots showed a clear separation of the two groups of tea samples in both positive and negative ionization modes ( Figures 4C, D). The Q 2 values (Table 1) of the OPLS-DA model based on the two modes were close to 1, reflecting excellent predictive ability. 200 times permutation tests (Figures 3B, C) with Q 2 intercept of −0.407 and −0.703 also presented that the OPLS-DA models have good reliability without overfitting. Based on the OPLS-DA models, secondary metabolites with VIP > 2, pvalues < 0.05 and fold change >2 or <0.5 were selected as differential .
/fsufs. .  compounds. As a result, 144 positive ions and 99 negative ions were screened out.

. . Integrated lipidomics and metabolomics for green tea samples analysis
In order to determine whether the identified characteristic compounds are representative to distinguish harvest periods of green teas, three reconstructed OPLS-DA models were generated based on the 26 differential lipids ( Figure 5A), 19 differential markers in positive ion mode ( Figure 5B) and 26 differential markers in negative ion mode ( Figure 5C), respectively. The values of R 2 X cum , R 2 Y cum , and Q 2 cum in the reconstructed OPLS-DA models were listed in Table 1. Compared to the original models, the reconstructed models have higher R 2 X cum values, which mean the removal of redundant information. The Q 2 cum values of the reconstructed OPLS-DA models were close to that of the original models, demonstrating the usefulness of these     screened differential components. Compared with the differential lipids or metabolites alone based model, the integrated model has better differentiation for the two groups of samples ( Figure 5D) with R 2 X cum 0.590, R 2 Y cum 0.961, and Q 2 cum 0.893. Permutation tests further indicated that these OPLS-DA models were reliable ( Figure 6). External validation was performed using 20% of the total samples as a testing set. The correct classification rates of the three . /fsufs. . reconstructed models for the two groups of tea samples were 100% (Table 1). The corresponding Y pred values were listed in Table 4. As shown in Table 4, the classification of M2 is ambiguous in the reconstructed OPLS-DA model based on differential lipids, while larger Y Pred (May) values (Table 4) were obtained by the integrated model, indicating the improvement of the robustness of the model. The above results imply that the selected differential lipids and metabolites are very reliable. What's more, the integration of differential compounds is effective for the discrimination of green tea harvested at different times in spring.

. . Changes of di erential lipids in green tea with harvest time
The variation of differential lipids with harvest periods was visualized by using a heatmap ( Figure 7A). Each row represents a differential lipid, and each column represents a tea sample. The color grading from green to red indicates the relative intensity of differential lipids changes from low to high levels. The tea samples were divided into two groups correspondingly except for two samples ( Figure 7A-a). In ET samples, the concentrations of most glycerophospholipids and several acylglycerolipids were higher, while glycoglycerolipids and several acylglycerolipids, sphingolipids, chlorophylls and derivatives are richer in LT. These characteristic The concentration of Hydroxypheophytin (HydPhe) a exhibited a much higher level in LT ( Figure 7A-b), which might be due to that chlorophyll is unstable in high temperature conditions. A relatively low temperature environment was more conducive to chlorophyll preserving and greenness retaining (Yu et al., 2019), while the high temperatures of picking tea in late spring periods may result in a degradation reaction of chlorophyll and further lead to the formation of more HydPhe a (Wei et al., 2011). Of note, chlorophylls generally exhibited bright green color, while the HydPhe displayed an undesirable dark and olive-brown color, resulting in LT color quality inferior to ET in view of the sensory quality.
SQDG is one of the predominant membrane lipids in thylakoids that existed on the photosynthetic membrane of higher plants .

FIGURE
HCA heatmap of di erential lipids (A) and di erential metabolites (B) between ET and LT samples. (Chang et al., 2014), which is the only lipid on the photosynthetic membrane that is negatively charged under physiological conditions (Benning, 1998). In this study, 4 species of SQDGs exhibited a higher level in LT samples ( Figure 7A-c). The increase of SQDGs is probably associated with the increase of photosynthesis in May because of more sunshine and higher temperatures, which is in accordance with . /fsufs. . the previous study that reported high temperature action caused an SQDG accumulation in wheat leaves (Taran et al., 2000). Glycerophospholipids, as major components of the plasma membrane, exhibited a higher content in ET. A previous study showed that most of PCs and PEs performed negatively correlated with the temperature and sunshine exposure . The lower contents of glycerophospholipids (PCs, PEs, PIs) in LT could be attributed to lipids degradation caused by the high-temperature action and long-sunshine exposure. In addition, other extraplastidial membrane lipids also performed significant differences between the two groups of tea samples. Ceramide [Cer 20:1/24:2(3O)], a major plant sphingolipid, showed a lower level in ET samples in Figure 7A-d. Acylglycerolipids included TGs and DGs were primarily efficient storage form of energy (Hu et al., 2008). Among them, DG (16:1/18:2) and DG (18:3/20:5) in group II exhibited a higher content in ET compared to that in LT, TG (15:3/16:5/20:6) and DG (22:1/18:2) showed opposite variation trend ( Figure 7A-e).

. . Changes of di erential metabolites in green tea with harvest time
The tea samples were divided into two categories clearly in the heatmap ( Figure 7B-a) based on differential metabolites. As for the identified 45 metabolites, 42 of them exhibited higher contents in ET; only 3 metabolites have higher levels in LT. Free amino acids are closely related to the aroma and taste of tea and existed positively correlated with the quality of tea (Lee et al., 2013). Theanine was confirmed as the most abundant free amino acid in tea (Deng and Ashihara, 2015). In Figure 7B-b, higher levels of theanine and uridylic acid were contained in ET. This result corresponded to the report of Liu et al. (2016), in which it reported that within the relatively short duration of the spring tea season, the levels of theanine showed obviously decreased during late spring harvest periods.
Flavonoids and flavonol/flavone glycosides were widespread secondary metabolites in plants with many significant biological functions. These compounds mainly existed in tea based on glycoside forms of kaempferol, myricetin and quercetin accounting for 2-3% (Lee et al., 2008). It was reported that flavonoids and their glucosides confer extremely low taste of mouth-dry and astringency (Scharbert et al., 2004;Scharbert and Hofmann, 2005). In the current study, the contents of flavonoids and flavonol/flavone glycosides were notably changed with respect to the harvest period of tea in the spring season. The contents of several flavone glycosides, such as myricetin 3-rutinoside-7-glucoside, myricetin 7-(6 ′′ -galloylglucoside), kaempferol-3-O-rutinoside, kaempferol 3-(2 ′′ , 6 ′′ -di-(E)-p-coumaryglucoside) in ET samples were higher than those in LT samples ( Figure 7B-c). In addition, myricetin, myricitrin, rutin, and other flavonoids in early spring tea samples were higher than those in late spring tea samples. However, 3' ,8-dimethoxyapigenin 7-glucoside and procyanidin B2 present significant increase in LT, which were believed to be the contributor of the astringent taste of tea (Dai et al., 2015). A possible reason is that higher temperatures and sun exposure in the late spring season enhanced the activity of some enzymes and the expression of structural genes encoding of flavonoids, and therefore result in the contents increase of some flavonoids (Sharma et al., 2010;Cheynier et al., 2013). In a word, flavonoids and flavonol/flavone glycosides exhibited diversity and concentration variations in green tea with respect to different harvesting times. The results were in consistent with that reported in the reference (Zeng et al., 2020).
Catechins are the main components in green tea that account for 60-80% of all polyphenols and contribute to the bitter and astringency taste of tea. In this study, catechin exhibited a higher level in LT ( Figure 7B-d), this result was consistent with the findings among the teas from different periods of spring (Liu et al., 2016). Phenolic acids and derivatives are a kind of aromatic compounds containing carboxyl and hydroxyl groups that play an important role in tea taste. Yang et al. (2018) found that the content of trigalloyl-glucose has a strong positive correlation with the grade of white tea. In the present study, 1,4,6-trigalloyl-β-D-glucopyranose was showed a higher level in ET compared to that in LT, shown in Figure 7B-e, which was in accordance with the previous reports (Tan et al., 2017). It was reported that the white tea harvest in early spring contained more levels of trigalloyl-glucose than that harvest in the late spring season.

. Conclusions
In this study, UPLC-Triple-TOF/MS-based non-targeted lipidomics and metabolomics methods were used to explore the difference of lipids and metabolites between ET and LT. With the aid of chemometric tools, characteristic compounds between the two groups of tea samples were revealed, and the results demonstrated that the differential components can distinguish ET and LT. A total of 70 characteristic components including 25 lipids and 45 secondary metabolites were identified. The correct classification rates of the reconstructed model based on the combination of differential lipids and metabolites for the two groups of tea samples was 100%. These results indicated that lipidomics and metabolomics coupled with chemometrics is a useful strategy to reveal the differences between ET and LT, and has the potential to be used for the quality evaluation of tea.

Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions H-NC: conceptualization, investigation, software, visualization, and writing-original draft. H-WG: funding acquisition, resources, supervision, software, validation, and review and editing. Z-QL: data curation and formal analysis. WS: supervision and writing-review and editing. BD, ZL, YC, and WL: writing-review and editing. X-LY: funding acquisition, methodology, project administration, and writing-review and editing. HF: conceptualization, supervision, funding acquisition, and validation. All authors contributed to the article and approved the submitted version.