Skip to main content

ORIGINAL RESEARCH article

Front. Nutr., 14 April 2022
Sec. Food Chemistry
This article is part of the Research Topic The Effects of Food Processing on Food Components and Their Health Functions, Volume II View all 10 articles

Assessment of Variations in Round Green Tea Volatile Metabolites During Manufacturing and Effect of Second-Drying Temperature via Nontargeted Metabolomic Analysis

\nHuajie Wang,&#x;Huajie Wang1,2Yaya Yu&#x;Yaya Yu1Wen Ouyang,Wen Ouyang1,3Yongwen JiangYongwen Jiang1Jinjin WangJinjin Wang1Jinjie Hua
Jinjie Hua1*Haibo Yuan
Haibo Yuan1*
  • 1Key Laboratory of Tea Biology and Resources Utilization, Ministry of Agriculture, Tea Research Institute, Chinese Academy of Agricultural Sciences, Hangzhou, China
  • 2Department of Tea Science, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou, China
  • 3State Key Laboratory of Tea Plant Biology and Utilization, College of Tea and Food Science and Technology, Anhui Agricultural University, Hefei, China

Round green tea (RGT) is widely distributed and presents a high yield in China. The quality of RGT can be determined by its aroma; however, the transformation and formation of volatile metabolites during RGT processing remain unclear. In this study, 173 volatile compounds (nine categories) were identified totally from RGT via gas chromatography-mass spectrometry with infrared-assisted headspace-solid phase microextraction. These substances exhibited different changing trends during various procedures, with the most intense transformation occurring during fixation, followed by pan-frying and second drying; moreover, 51 substances were screened, mainly containing fatty acid-derived volatiles (i.e., (E)-2-hexen-1-ol, Hexanal, pentanal, hexanal) and glycoside-derived volatiles (i.e., linalool, geraniol, benzyl alcohol, benzaldehyde), and their evolution during processing was clarified. Furthermore, the effect of the second-drying temperature on volatile compound metabolism was clarified, and 90°C was the best temperature for RGT aroma. This research lays a foundation for in-depth quality control and the aroma formation mechanism of RGT.

Introduction

Green tea has gained immense interest among consumers worldwide due to its green leaves and soup color, high fragrance, brisk and mellow taste, and various physiological activities such as anticancer, antioxidant, antiradiation, and antibacterial properties (14). Green tea can be classified into flat-shape, curled-shape, round (bead)-shape, needle (bud)-shape, and strip-shape, based on its appearance type (5). Round green tea (RGT), such as Yongxi Huoqing green tea in Anhui province, Fenghua curled tea and Pingshui Rizhu tea in Zhejiang Province, and Emerald tea in Guizhou Province, are widely consumed in China owing to their unique quality characteristics including tight, heavy, and sturdy round shape, green glossy color, high fragrance, and brisk and mellow taste (6, 7). Moreover, RGT accounts for a large proportion of green tea exported from China and is consumed by people worldwide.

Aroma is a key factor in evaluating the quality of green tea and affects its market value. The volatile compounds in green tea mainly include alcohols, aldehydes, ketones, esters, aromatic hydrocarbons, and terpenes (8). In tea, the volatile compounds arise from two sources—fresh tea leaves and substances formed during processing. These substances can be classified into four categories from different precursors: amino acid-derived volatiles, wherein amino acids undergo a Strecker reaction under the action of heat, and Maillard reaction with sugars to produce aldehydes, ketones, furans, pyrazines, and high-boiling volatile compounds (9, 10); glycoside-derived volatiles and glycoside-derived Non-volatiles, which release bound volatile compounds such as linalool and geraniol under the action of glycosidase or heat (11); lipid-derived volatiles, wherein lipids are first degraded to produce linolenic acid and other fatty acids, followed by production of lower-pointing alcohols and aldehydes (12); and carotenoid-derived volatiles, wherein carotenoids are precursors of volatile compounds such as ionone, 2,6,10,10-tetramethyl-1-oxaspiro[4.5]dec-6-en-8-one, and β-tanshinone (13).

The combination of processing technologies directly affects the formation pathways of volatile compounds and the type and quality of the green tea aroma (14). Moreover, tea leaves with distinct shapes made using different processing techniques have different cell breakage rates and material leaching rates, resulting in various aroma types, for example, flat and striped teas often exhibit floral and chestnut aromas, and needle-shaped teas often reveal delicate aromas (15, 16). Previous studies have reported the aroma of flat striped, needle-shaped, curly-shaped green teas (7, 17, 18); however, little research has been conducted on the aroma of RGT. The processing technology of RGT includes spreading, fixation, rolling, pan-frying, roll-roasting, and second-hot air drying. Among these, spreading promotes the loss of low-boiling compounds, such as (z)-3-hexen-1-ol and (E)-2-hexenal (19). Fixation may increase the formation of high boiling point substances, such as 3-methyl-butanal, linalool, and trans-β-ionone, under high temperature (8). Pan frying, the key process for shape formation and internal quality of RGT, is generally performed twice or thrice, with an initial increase in temperature and time, followed by a gradual decrease. Repeated pan frying helps to form round, tight, and heavy shape, and can trigger the biochemical reaction in the leaves under prolonged heat, which is conducive to the transformation and retention of quality components and lays the foundation for the formation of excellent quality. Roll-roasting and second-hot air drying promotes further emission of grass-flavored substances, formation of high-boiling substances, and the Maillard reaction of amino acids and soluble sugars to form heterocyclic substances such as furans and pyrazines at high temperatures (20).

The current research on RGT only focuses on the effect of processing technology on the sensory evaluation results and the total amount of physical and chemical components (21, 22). Non-targeted metabolomics is a novel analytical tool that allows the simultaneous determination of numerous substance components and comprehensive analysis of metabolites changes. It has been widely used in the analysis of Non-volatile and volatile metabolites of tea in recent years (2325). However, few studies have reported the formation of aroma quality of RGT and volatile metabolite variations during manufacturing; moreover, the effect of RGT processing technology on the metabolism of volatile compounds remains to be elucidated. Therefore, this study uses the “Jiu Keng” fresh tea leaves as the raw material to process RGT and to obtain samples of the whole process (spreading, fixation, rolling, pan-frying, roll-roasting, and second drying). In addition, three second-drying temperatures (70°C, 90°C, and 110°C) were set to make RGT. Volatile compounds of these tea samples were detected using infrared-assisted-solid-phase microextraction (IRAE-HS-HPME) and gas chromatography-mass spectrometry (GC-MS) platforms. Multiple statistical analysis methods were used to elucidate the dynamic evolution of volatile compounds and the key metabolic pathways of volatiles biosynthesis during the processing of RGT, and the influence of different second-drying temperatures on volatile metabolites conversion. This study can provide a theoretical reference and objective basis for the mechanism underlying the formation of the aroma of high-quality RGT.

Materials and Methods

Experimental Materials

Fresh tea leaves (FTL) of JiuKeng (Camellia sinensis L.) were collected from the Zhejiang Kaihua tea production area (29°10′ N, 118°10′ E; at an elevation of 700 m a.s.l.) in April 2020. The moisture content of the FTL was ~78%. Each shoot comprised two leaves and one bud.

All the chromatography-grade volatile standards were purchased from Shanghai Yuanye Biological Technology Company (Shanghai, China). The purity of these chemical standards was ≥98%.

Green Tea Processing

(1) FTL (25 kg) were spread 2–3 cm thickness on each wooden net plates (1 m × 1 m) for 12 h at the temperature 19–22°C with a relative air humidity of 60% environment until the moisture reached 70.0%.

(2) Fixation was performed using a roller fixation machine (80-cm type, electromagnetic roller fixation machine, Ningbo Yaojiangyuan Machinery Co., Ltd., Ningbo, China), the temperature of rolling was 300°C, the speed of rolling was 22 rpm, the fixation time was 120 s, and 120 kg of tea leaves were treated per hour. The fixated tea leaves (FTS) spreading and resurgence was conducted at 0.5 h on the wooden net plates.

(3) The cooled fixated tea leaves (FLS) were subjected to rolling for 30 min in a rolling machine (6CR-50 type, Zhejiang Shangyang Machinery Co., Ltd., Quzhou, China) using the following sequences: no-pressure rolling (5 min), light-pressure rolling (10 min), intermediate-pressure rolling (10 min), and no-pressure rolling (5 min).

(4) The rolled leaves (RTL) were pan-fired thrice in a double-pan roasting machine (6CCGQ-50 type, Zhejiang Shangyang Machinery Co., Ltd., Quzhou, China). First, RTL was treated for 20 min under a pan temperature of 210°C with a fried speed of 38 r/min until the leaves achieved a moisture content of 27.5%, as the pan-fried 1 leaves (PF1); thereafter, PF1 was treated for 40 min at a pan temperature of 180°C with a fried speed of 32 rpm until the leaves reached a moisture content of 15.0%, as the pan-fried 2 leaves (PF2); and PF2 was treated for 25 min at a pan temperature of 200°C with a fried speed of 32 rpm until the leaves had a moisture content of 10.0% (pan-fried 3 leaves (PF3).

(5) The pan-fried three leaves (PF3) were depilated for 15 min in a depilating machine (100-type, Zhejiang Shangyang Machinery Co., Ltd., Quzhou, China).

(6) The depilated leaves (DTL) were dried at 100°C for 15 min in a six-edge roasting machine (6CCP-60 type, Zhejiang Shangyang Machinery Co., Ltd., Quzhou, China) until 7.50% of leaf water content was achieved; thereafter, the leaves were spread out for a further 30 min, as the rotary pot first-dried leaves (RPL).

(7) Eventually, the samples were independently dried at 70°C, 90°C, and 110°C for 30 min until complete dryness was achieved in a box hot air-drying machine (Zhejiang Shangyang Machinery Co., Ltd., Quzhou, China); these dried at 70°C, 90°C, 110°C in box hot-air as the second-dried leaves (BHL70, BHL90, BHL110).

The overall production process of RGT is illustrated in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1. Flow chart of round green tea processing.

Analysis of the Volatile Compounds Present in the RGT Samples

Sample Pretreatment

The green tea samples (0.5 g) were accurately weighed (0.001 g), and 20 mg/L ethyl caprate (10 μL) were added to a 20 mL headspace vial along with boiled distilled water (1 mL). The headspace vial cap was immediately tightened, and the protective cap was impaled using a DVB/CAR/PDMS fiber (2cm) head with a manual handle (SPME, Supelco, Commonwealth of Pennsylvania, USA). Thereafter, the vial was heated over a 100 W infrared (IR) device (Qiyi Lighting Company, Zhejiang, China) for 15 min, and the fiber head was inserted into the GC-MS inlet and desorbed at 250°C for 5 min. Each sample was subjected thrice to IR-assisted headspace-solid-phase microextraction (8).

GC-MS Analysis

The gas chromatograph (7890B−7000C, Agilent Technologies, Palo Alto, CA, USA) equipped with a HP-5 ms ultra-inert capillary column (30 m × 0.25 mm × 0.25 μm) was used in the splitless injection mode with a high-purity helium (99.999%) flow at 1.0 mL/min for GC-MS analysis. The GC injector temperature was maintained at 250°C, and the temperature program of following column was employed: initial temperature, 50°C (held for 5 min), increased to 150°C at a rate of 4°C/min (held for 2 min), and increased to 270°C at a rate of 10°C/min (held for 6 min) (26, 27).

Qualitative and Quantitative Analyses

For qualitative analysis, the volatile metabolites with a similarity of >80% to the NIST11 standard library (22) were screened using the Agilent Mass Hunter unknown analysis program. Kovàts retention indices for each compound were obtained by calculating the linear formula of n-alkanes (C7-C40) and by comparing with theoretical retention indices referred to the literature value with the same or equivalent chromatographic column on the NIST Chemistry WebBook (http://webbook.nist.gov/chemistry/), based on the difference of RI within 30 (8). Moreover, some key aroma substances, such as 1-hexanol, (E)-4-heptenal, octanal, 1-decanol, linalool, naphthalene, decanal, geraniol, indole, α-ionone, cedrol, trans-β-ionone, hexadecenoic acid ethyl ester, and linoleic acid ethyl ester were identified further with authentic standards. For quantitative analysis, the mass concentrations of the volatile metabolites were calculated regarding the internal standard method using equation:

Ci = (Cis × Ai)/Ais    (1)

Ci represents the concentration of any volatile metabolites (μg/L), Cis is the concentration of the internal standard (μg/L) (ethyl caprate, 20 mg/L), Ai is the chromatographic peak area of any metabolites, and Ais represents the chromatographic peak area of the internal standard (27).

OAV

The OAV is the ratio of the concentration of volatile compounds in a sample to the odor threshold (28), calculated according to equation:

OAV = c/OT    (2)

c is the concentration of the metabolites in the green tea sample (μg/L), and OT represents the odor threshold of the corresponding components in water (μg/L).

Statistical Analysis

The tests were repeated in triplicate, and the results of each test were expressed as the average of three replicates. The analyze the principal components affecting the volatile metabolites according to the different green tea processing using SIMCA P13 software (Umetrics, Umea, Sweden). MEV 9.0 software (https://mev-tm4-org.caas.cn) was used to generate heatmaps of key differential metabolites. SAS software (version 9.4; SAS Institute Inc., Cary, USA) was used to analyze the significant differences in the concentration of volatile compounds according to different treatments.

Results

Analysis of Volatile Compounds During the Processing of Green Tea

Overview of the Profile of Volatile Compounds in RGT

In total, 173 volatile compounds were identified (Supplementary Table S1) in this study, including 32 alcohols, 14 aldehydes, 17 ketones, 31 esters, 16 terpenes, 15 aromatic hydrocarbons, 28 alkanes, 3 furans, and 17 other types.

Moreover, a Non-supervised metrology PCA tool was used to conduct a comprehensive analysis of these 173 volatile compounds, as depicted in the score plot (Figure 2A). All tea samples revealed marked separation and clustering. FTL and STL are on the left side of the score plot, whereas the samples of FTS and subsequent samples are on the right. The trajectory plot (Figure 2B) clearly depicts stepwise alterations and distinct differences in volatile compounds from fresh leaves to finished leaves. Furthermore, a great change was observed in volatile compounds after fixation, and a further shift occurred in the pan-fried and box hot air drying; however, minor variations were observed during rolling, depilation, and rotary pot first-drying.

FIGURE 2
www.frontiersin.org

Figure 2. (A) Principal Component Analysis (PCA) score scatter plot composed of all volatile compounds in the round green tea process. (B) The dynamic metabolic trajectory plot of volatile compounds, drawn by scatting of the Principal Component Analysis (PCA) scores t[1] vs. t[2]. (C) Variation trend chart of different types of volatile compounds content in different processes of round green tea. (D) Variation of the proportion of different types of volatile compounds in different process of round green tea.

Dynamic Changes of Volatile Categories During RGT Processing

Nine categories of volatile compounds underwent remarkable changes during RGT processing. Changes in the content of different categories of volatile compounds are illustrated in Figure 2C; the contents of alcohols and aldehydes slightly increased in STL and markedly decreased in FTS; that of aldehydes significantly increased in RTL and PF2 and then decreased (p < 0.05); and the ketones markedly decreased in STL, FTS, and BHL, and increased in FLS and RPL. The ester content decreased throughout the process, and the sharpest decline was observed after fixation; terpenes increased in STL, significantly decreased in FTS, and slightly fluctuated in FLS (p < 0.05); and aromatic hydrocarbons revealed an alternate decreasing and increasing trend in FTS. The content of furans initially increased before the third pan-fried process and then decreased in DTL and BHL. The contents of alkanes and others (such as benzoic acid, hydrazide, benzyl nitrile, and 2-hexyl-5-pentylpyrrolidine) showed downward trends, particularly in fixation. The total volatile compound content revealed following trend: FTL > STL > RTL > FLS > FTS > PF2 > PF3 > PF1 > RPL > DTL > BHL (Supplementary Table S1). From the proportion of categories (Figure 2D), alcohols accounted for the highest proportion in the whole process (30.9–37.5%); moreover, the proportion of aldehydes increased significantly in RTL, PF2, and BHL (p < 0.05), whereas that of esters decreased in the process (p < 0.05). The proportion of ketones increased significantly in the FTS and decreased in the BHL (p < 0.05). The terpenes increased in the STL and then decreased. It is speculated that the activity of glycoside hydrolase increased during spreading, thereby promoting the release of terpene volatile compounds. Furans and aromatic hydrocarbons occupied a small proportion in the whole process (0.9–4.8%).

Differential Compound Screening and Evolution Analysis During RGT Processing

Analysis of the volatile compounds indicate that volatile compounds undergo drastic changes during processing. To clearly distinguish the important compounds in the process and to determine the law of their evolution, 51 substances were screened via variable projection importance (VIP) in the PLS-DA analysis (VIP > 1) and one-way analysis of variance (p < 0.05), as depicted in the loading diagram of PLS-DA (Figure 3A). Moreover, 51 substances mainly include alcohols, aldehydes, esters, ketones, aromatics hydrocarbons, terpenes, and furans, such as (4E)-4-hexen-1-ol, acetate, hexanal, heptanal, octanal, benzene acetaldehyde, 3-methyl-butanal, 2-methyl-butanal, benzene acetaldehyde, 2,5-dimethyl-benzaldehyde, 2-heptanone, α-cadinol, 3,7-dimethyl-1,5,7-octatrien-3-ol, 1-penten-3-ol, 1-pentanol, (E)-2-hexen-1-ol, acetate, styrene, (Z)-hexanoic acid, 3-hexenyl ester, α-ionone, and isophorone (Figure 3B).

FIGURE 3
www.frontiersin.org

Figure 3. (A) Loading diagram of Partial Least Squares Discrimination Analysis (PLS-DA) model. (B) Heat map of 51 difference volatile compounds.

Evolution of Key Different Compounds and Their Related Compounds During Processing

From the 51 compounds obtained in section Differential Compound Screening and Evolution Analysis During RGT Processing, the metabolic transformation rules of 34 key volatile compounds with definite aroma characteristics (9, 13) were revealed. Combined with the transformation and metabolism-related substances, this study clarified the evolutionary law of total 48 volatile compounds in the process of RGT, which contained 21 fatty acid-derived volatile substances, 18 glycoside-derived volatile substances, 7 amino acid-derived volatile substances, and 3 carotenoid-derived volatiles (Supplementary Figure S1).

Fatty Acid-Derived Volatiles

Lipids are important precursors for the formation of volatile compounds in green tea. The degradation products of unsaturated fatty acids or saturated fatty acids, such as α-linolenic acid and linoleic acid, can form C6-aldehydes under the action of heat or lipoxidase (LOX), and then further undergo redox reaction, isomerization, and may form esters, as well as chemical reactions to generate corresponding ketones, and alcohols (8). Figure 4A illustrates such compound changes during processing of RGT: under the action of lipoxidase, the hexanal, pentanal and heptanal derived from linoleic acid showed an upward trend, whereas (E)-2-hexen-1-ol, acetate decreased sharply. Hexanal, 1-hexanol, 1-pentanol, 1-heptanol and (E)-3-hexen-1-ol derived from linoleic acid and α-linolenic acid decreased significantly after fixating, whereas (E)-2-hexen-1-ol, acetate showed increased trend. During the spreading process, under the action of hydrolases, oxidases, and isomerases, hydrolysis of esters is increased, and the degradation of lipids results in the formation of corresponding alcohols and aldehydes, under the heat of fixation, the oxidation and esterification were promoted, these results are consistent with the existing findings (19, 29). High boiling point compounds, such as pentanal, heptanal, hexanal, and 2-heptanone, 6-methyl-5-hepten-2-one, and 1-penten-3-ol (Supplementary Figure S1) increased significantly during pan-fried process (p < 0.05), likely because the process of rolling promoted the fragmentation of tea tissue cells, resulting in a large amount of unsaturated fatty acids overflowing. Thus, these unsaturated fatty acids underwent strong thermal oxidation reactions under long-term high temperature treatment of pan-fried to form alcohols, aldehydes with high boiling points. Furthermore, heat fixation and pan frying increased the volatilization of low-boiling alcohols and induce the formation of corresponding aldehydes, ketones esters by oxidation, isomerization, esterification reactions.

FIGURE 4
www.frontiersin.org

Figure 4. Evolution of compounds from different sources: (A) Fatty derived volatiles, (B) Glycoside-derived volatiles, (C) Amino acid-derived volatiles, and (D) Carotenoid-derived volatiles.

Glycoside-Derived Volatiles

Glycoside-derived volatiles generally exhibit floral, sweet, and fruity aromas and form an important component of green tea aroma (11). During the process, monoterpene alcohols (linalool, and geraniol) or aromatic alcohols (benzyl alcohol and phenylethanol) can be liberated by glycosidic enzymes or heat reaction. In this study, the levels monoterpene alcohols, such as linalool, geraniol markedly decreased after fixation (Figure 4B), whereas the dehydrogenation compound (3,7-dimethyl-1,5,7-octatrien-3-ol), isomerisation compound (nerol) considered important components of the floral and fruity aroma of green tea (28), remarkably increased. The levels of methyl salicylate, alfa-copaene, cubenene, α-cadinol, etc., (Supplementary Figure S1) markedly decreased in FTS, whereas those of citral, as the oxide of nerol, which is the isomer of geraniol and linalool, were significantly increased in PF2 and BHL, particularly in PF2 (p < 0.05). This indicates that the isomerization of linalool, geraniol and the oxidation reaction of nerol may generate citral under long-term high-temperature treatment of PF2 (11). The level of aromatic compounds, such as benzyl alcohol, benzaldehyde, styrene, phenylethyl alcohol markedly decreased after fixation, while corresponding oxides (benzene acetaldehyde) was increased. Benzoic acid, ethyl ester was increased significantly after rolling process, which may relate to the increase of cell fragmentation rate under the action of external force and high leaching rate. The high temperature treatment of RPL promoted the oxidation and methylation reaction to form acetophenone, 2,5-dimethyl-benzaldehyde. Many glycoside-derived volatiles disappeared during the fixation process, which is consistent with the results of previous studies (11). It is speculated that in a high-temperature environment, some volatiles are lost and recombine with glycosides to form nonvolatile compounds; furthermore, they undergo oxidation, isomerization, and other reactions (8).

Amino Acid-Derived Volatiles

Amino acid-derived volatiles are generally derived from the Strecker reaction of amino acids and the Maillard reaction of amino acids with sugars (9). Acid-derived volatiles were identified in this study (Figure 4C). Among these, 2-methyl-butanal and 3-methyl-butanal increased significantly after spreading (p < 0.05), which is related to the reaction of leucine and isoleucine with glutamine under the action of transaminase and decarboxylase. After fixation, the levels of 3-methyl-butanal, 2-methyl-butanal, and 2-methyl-butanol were significantly reduced (p < 0.05), whereas those of 3-methyl-butanol significantly increased (p < 0.05). Moreover, at high temperatures, oxidation reaction was triggered, and strong reduction reaction and isomerization reaction also occurred (10); however, their contents in PF3 and BHL have rebounded, which is related to the high temperature and long-term treatment at this stage to promote the Strecker reaction of amino acids and oxidation and isomerization. 2-methoxy-furan, cis-2-(2-pentenyl) furan, and 2-propyl-furan are the Maillard reaction products; the levels of 2-methoxy-furan increased in FTS, PF2, BHL, particularly in BHL, which is consistent with the results of the existing research (9), and 2-propyl-furan exhibited an upward trend in most processes but decreased in DTL and BHL. Moreover, cis-2-(2-pentenyl) furan increased in STL, FTS, and RPL, whereas it decreased in BHL significantly (p < 0.05), which contradicts the existing research that high temperature and long-term drying is conducive to the Maillard reaction (30), and may be related to the further transformation of these two substances at high temperatures.

Carotenoid-Derived Volatiles

Carotenoids can generate ionone, 2,6,10,10-tetramethyl-1-oxaspiro [4.5] dec-6-en-8-one, and other floral and fruity substances under enzymatic oxidation and heat, and are important precursors of the formation of green tea aroma (9). The evolution of β-ionone, α-ionone and isophorone is depicted in Figure 4D; β-ionone, α-ionone increased during FLS and RPL, and the decrease was most obvious at BHL (29). Isophorone is another carotenoid-derived volatile (31), which exhibited decreased trend before PF2, but increased mostly in PF2 and further decreased. The rebound in PF1 and PF2 is related to the increase in cell disruption and the overflow of its precursor substances, which further promotes its synthesis under the action of high temperature (32).

Effect of Box Hot Air Second-Drying Temperature on Volatile Compounds

In section Analysis of Volatile Compounds During the Processing of Green Tea discussion, it was found that the hot air drying process is crucial in the formation of aroma during RGT processing. Therefore, this study will further discuss the influence of different second-drying temperatures on the volatile compounds of RGT.

Volatile Profile of Green Tea Under Different Box Hot Air Second-Drying Temperature

In total, 247 volatile compounds were obtained in different second-drying temperature samples (Supplementary Table S2). The comparison of different temperature samples revealed that alcohols and furans were significantly higher in BHL70 (p < 0.05), whereas alkanes, terpenes, aromatic hydrocarbons, and pyrazines were higher in BHL110 than in other samples, and the contents of various compounds (such as alkanes, terpenes, aromatic hydrocarbons, alcohols, and furans) were moderate in BHL90 (Figure 5A). Substances with green aroma, such as (E)-3-hexen-1-ol, 1-hexanol, 1-decanol, hexanal, (E)-2-hexenal, and (E)-4-heptenal were higher in BHL70; naphthalene, 2-methoxy-furan, (E)-2-pentenal, 3-octen-2-one, and 2,6,10,10-tetramethyl-1-oxaspiro [4.5] dec-6-ene were higher in BHL90; whereas 3-ethyl-1H-pyrrole, 2-ethyl-6-methyl-pyrazine, 3-ethyl-2,5-dimethyl-pyrazine, and roast aroma were higher in BHL110 than BHL90 and BHL70 (Supplementary Table S2).

FIGURE 5
www.frontiersin.org

Figure 5. (A) Histogram of different types of volatile compounds of round green tea by different second-drying temperature. (B) Partial Least Squares Discrimination Analysis (PLS-DA) score. (C) Loading diagram for volatile compounds of round green tea by different second-drying temperature. (D) Heat map of 117 difference volatile compounds. Different letters indicate the significance of the difference of the same type substance at different temperatures through LSD test (p < 0.05).

PLS-DA

To further clarify the effect of the second-drying temperature on volatile compounds and obtain key differential compounds, PLS-DA analysis (33) was conducted in this study based on 247 volatile compounds. The score map is illustrated in Figure 5B, and three groups of samples were distinguished. While detecting t [1], BHL70 and BHL90 were on the right side of the figure, and BHL110 was on the left side of the figure. The model parameters were R2X = 0.941, R2Y = 0.998, and Q2 = 0.988, indicating that the model is stable and reliable and has strong predictive ability.

To obtain the key volatile compounds that distinguish samples at different second-drying temperatures, variable importance of projection (VIP) and difference significance analysis were used to screen differential substances. In total, 117 substances (load diagram, Figure 5C), such as benzaldehyde, octanal, D-Limonene, naphthalene, hexanal, 3-ethyl-1H-Pyrrole, (E)-3-Hexen-1-ol, 2,5-dimethyl-pyrazine, and trans-β-ionone, were obtained. Supplementary Figure S2 illustrates the changes in 117 substances under different second-drying temperatures. The heatmap (Figure 5D) showed the changes of 117 differential compounds under different second-drying treatments clearly. Octanal, (E, E)-2,4-heptadienal, 1-Octanol, nonanal, and (E)-4-heptenal, with a grassy aroma, were higher in BHL70; 2-methoxy-furan, (E)-2-pentenal, 3-octen-2-one, and 2,6,10,10-tetramethyl-1-oxaspiro [4.5] dec-6-ene were higher in BHL90; and 3-ethyl-1H-pyrrole, 2-ethyl-6-methyl-pyrazine, and 3-ethyl-2,5-dimethyl-pyrazine were higher in BHL110. This result follows the previous studies reporting that significantly lower temperatures lead to greater retention of grass-flavored substances and insufficient formation of high-boiling aromatic substances (p < 0.05), whereas higher temperature can significantly promote the Maillard reaction of amino acids and sugars to produce too many pyrazines, which imparts a baking flavor to the tea (34).

OAV

OAV analysis is frequently applied to evaluate the contribution of volatile compounds in aroma (27, 35). Odor thresholds and characteristics of the 41 patients (Supplementary Table S3) consuming the 117 substances were obtained based on previous studies (8, 3639). Of these, 21 substances (Table 1), with important contributions to aroma, were screened out according to the principle of OAV > 1. In particular, the OAV of trans-β-ionone is >9,000, and 1-octanol is >450, both of which are important components of aroma composition (8, 38). From the perspective of different second-drying temperature treatments, the OAVs of 1-hexanol, (E)-4-heptenal, octanal, 1-octanol, nerolidol 2, benzaldehyde with green, and flower aroma were higher in BHL70 than in the other treatments, whereas naphthalene, 2,6,10,10-tetramethyl-1-oxaspiro [4.5] dec-6-ene and decanal significantly contributed to BHL90 (p < 0.05). These compounds, such as 3,7-dimethyl-1,5,7-octatrien-3-ol, 2,5-dimethyl-pyrazine, ethyl-pyrazine, 2-ethyl-6-methyl-pyrazine, and 3-ethyl-2,5-dimethyl-pyrazine with roasted aroma are only present in BHL110; however, due to the high threshold of these substances, their OAV <1, but the relatively low OAV values of other substances highlight the baking flavor.

TABLE 1
www.frontiersin.org

Table 1. Aroma characteristics and OAV values of 21 different volatile compounds (OAV > 1) in different second-drying temperature samples.

Discussion

Processing technology and second-drying temperature are key factors in the formation of RGT aroma (8, 39). In this study, the concentration of all the identified volatile compounds have undergone substantial changes during tea processing, particularly during the fixation process, followed by pan-frying and second-drying, which is consistent with previous reports (7, 8). Under the high-temperature fixation treatment, many low-boiling point alcohols, aldehydes, and esters are lost and transformed, and violent degradation, isomerization, redox, and other reactions occur, consequently forming high-boiling point compounds. During the pan-fried process, the tea leaves and pot wall are repeatedly rubbed. Therefore, the cell fragmentation rate increases, and the volatile compounds and their precursors are further released from the cells. Moreover, high temperature increases thermochemical reactions and promotes the formation of furans, aldehydes, and aromatic hydrocarbons with high-boiling point. During the box hot air second-drying, few low-boiling volatile compounds were further lost. This could be related to the large consumption of precursor compounds in fixation and pan-frying process. The levels of fatty acid-derived volatiles and glycoside-derived volatiles have undergone considerable changes during RGT processing, especially during the fixation process. Under high temperature condition during the fixation process, the temperature of tea leaves changes from 20 ~ 25°C to 70 ~ 80°C, the moisture content of tea leaves changes from 70 to 30% ~ 40%, the state of enzymes changes from active to inactive, and a series of enzymatic chemical reactions and thermochemical reactions take place in the tea leaves in 2 ~ 5 min. High temperature can promote the degradation of lipids to generate free fatty acids, which generate alcohols and aldehydes volatiles under the effect of heat-induced oxidation. The increase in glycoside-derived volatiles in fixation process may be due to the simultaneous enzymatic hydrolysis and thermochemical degradation of glycosides. This result is consistent with the previous reports that lipids and glycosides undergo drastic changes during fixation (11, 29, 40).

The second-drying temperature also had a great effect on the levels of volatile compounds. When the temperature was too low (70°C), the dry-heat reaction was insufficient and several green grass-flavor alcohols and aldehydes were retained (11). Therefore, the tea sample exhibited a faint scent. However, when the temperature was too high (110°C), the dry-heat reactions (Maillard reaction, Caramelization reaction, etc.) were too intense and several baking-flavor volatiles were generated (34). At the appropriate temperature (90°C), the content of each compound was moderate and the proportion of various types of volatile compounds was coordinated. In brief, the temperature at which biochemical reactions occur at a moderate rate to avoid overfired flavor is the best temperature to produce high-quality RGT.

Conclusions

In this study, the evolution of volatile compounds during the processing of RGT was investigated using IR-HS-SPME technology combined with GC-MS, and the effects of different second-drying temperatures on the formation of volatile compounds were clarified. Fixation is considered the most drastic process of volatile metabolite conversion, followed by pan-fried and box hot air second-drying. Moreover, 51 key volatile substances were screened out, and fatty acid-derived volatiles and glycoside-derived volatiles changed most dramatically during processing. Substances with excellent aroma flavor, such as nonanal, decanal, tau-cadinol, naphthalene, and 2,6,10,10-tetramethyl-1-oxaspiro [4.5] dec-6-ene were significantly higher under BHL90 treatment (p < 0.05) and 90°C was the most suitable temperature for the second-drying processing of RGT. This study provides a theoretical basis and technical guidance for the processing of high-quality RGT and lays a foundation for in-depth exploration of the aroma formation mechanism during green tea processing. As only the precursors of part volatile metabolites are known, the number of classified volatile compounds discussed in this study is limited. In future research, we need to assess the volatile precursor substances, analyze the relationship between Non-volatile compounds and volatile compounds, and explore the mechanism of RGT aroma formation.

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

HW: conceptualization, investigation, formal analysis, software, writing–original draft, and editing. YY: investigation, resources, and writing–review and editing. WOY and JW: investigation. YJ: data curation software. JH: project administration, writing–review, and editing. HY: funding acquisition and supervision. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Science and Technology Innovation Project of the Chinese Academy of Agricultural Sciences (CAAS-ASTIP-TRICAAS), the open fund of State Key Laboratory of Tea Plant Biology and Utilization (SKLTOF20210110), and the National Science Foundation of China (31972465).

Conflict of Interest

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

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

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

Abbreviations

FTL, Fresh tea leaves; STL, Spread tea leaves; FTS, Fixated tea leaves; FLS, Spread Fixated tea leaves; RTS, Rolled tea leaves; PF1, Pan-fried 1 leaves; PF2, Pan-fried 2 leaves; PF3, Pan-fried 3 leaves; DTL, Depilated tea leaves; RPL, Rotary pot first-dried leaves; BHL, Box hot air second-dried leaves; BHL70, The temperature of box hot air second-dried is 70°C; BHL90, The temperature of box hot air second-dried is 90°C; BHL110, The temperature of box hot air second-dried is 110°C; RGT, Round green tea; IRAE-HS-SPME, Infrared-assisted headspace-solid phase microextraction; GC-MS, Gas chromatography-mass spectrometry; MEV, Multiple experiment viewer; PLS-DA, Partial least-squares discriminant analysis; VIP, Variable importance projection.

References

1. Aoshima H, Ayabe S. Prevention of the deterioration of polyphenol-rich beverages. Food Chem. (2007) 100:350–5. doi: 10.1016/j.foodchem.2005.09.052

CrossRef Full Text | Google Scholar

2. Yan HM, Zhong YZ, Duan YH, Chen QH, Li FN. Antioxidant mechanism of tea polyphenols and its impact on health benefits. Anim Nutr. (2020) 6:115–23. doi: 10.1016/j.aninu.2020.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Bag A, Mondal A, Majumder A, Banik A. Tea and its phytochemicals: hidden health benefits & modulation of signaling cascade by phytochemicals. Food Chem. (2021) 371:131098. doi: 10.1016/j.foodchem.2021.131098

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Shang A, Li JH, Zhou DD, Gan RY, Li HB. Molecular mechanisms underlying health benefits of tea compounds. Free Radic Biol Med. (2021) 172:181–200. doi: 10.1016/j.freeradbiomed.2021.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Bai XB, Li B, Deng HX. Research progress of green tea shaping equipment in China. Chin Tea Process. (2021) 2021:5–11. doi: 10.15905/j.cnki.33-1157/ts.2021.02.001

CrossRef Full Text | Google Scholar

6. Yu FL, Guo L, Zhu Y, Dai WD, Zhang Y, Du QZ, et al. Variation of aroma components of particle type green tea during the manufacturing process. Food Sci Technol Res. (2015) 40:45–50+56. doi: 10.13684/j.cnki.spkj.2015.10.010

CrossRef Full Text | Google Scholar

7. Hu X, Wei YM, Fang SM, Wang YJ, Xu SS, et al. Study on continuous shaping technology of granular green tea. Chin Tea Process. (2020) 2020:27–34. doi: 10.15905/j.cnki.33-1157/ts.2020.03.006

CrossRef Full Text | Google Scholar

8. Wang HJ, Hua JJ, Jiang YW, Yang YQ, Wang JJ, Yuan HB. Influence of fixation methods on the chestnut-like aroma of green tea and dynamics of key aroma substances. Food Res Int. (2020) 136:109479. doi: 10.1016/j.foodres.2020.109479

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Ho CT, Zheng X, Li SM. Tea aroma formation. Food Sci Hum Well. (2015) 4:9–27. doi: 10.1016/j.fshw.2015.04.001

CrossRef Full Text | Google Scholar

10. Wan XC. Biochemistry of tea. Third Edition. China Agricultural Publishing House, Beijing (2003). p. 160–80.

Google Scholar

11. Li PL, Zhu Y, Lu ML, Yang C, Xie DC, Tan JF, et al. Variation patterns in the content of glycosides during green tea manufacturing by a modification-specific metabolomics approach: Enzymatic reaction promoting an increase in the glycosidically bound volatiles at the pan firing stage. Food Chem. (2019) 279:80–7. doi: 10.1016/j.foodchem.2018.11.148

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Li J, Hua JJ, Yuan HB, Deng YL, Zhou QH, Yang YQ, et al. Investigation on green tea lipids and their metabolic variations during manufacturing by nontargeted lipidomics. Food Chem. (2021) 339:128114. doi: 10.1016/j.foodchem.2020.128114

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Chen QC, Zhu Y, Liu YF, Liu Y, Dong CW, Lin Z, et al. Black tea aroma formation during the fermentation period. Food Chem. (2022) 374:131640. doi: 10.1016/j.foodchem.2021.131640

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Dong C, He C, Zhu JY, Qu FF, Ran MQ, Ai ZY, et al. Effect of shaping methods on the quality of E-Cha No.10 green tea products. J Yunnan Agric Univ. (2020) 35:139–48. doi: 10.12101/j.issn.1004-390X(n).201805002

CrossRef Full Text | Google Scholar

15. Zhou TS, Mi XL, Wang YC, Yu YB, Li SL, Qi YG. Effect of different processing techniques on the quality of “Shaancha 1” green tea. Food science. (2017) 38:148–54. doi: 10.7506/spkx1002-6630-201703025

CrossRef Full Text | Google Scholar

16. Chen L, Ji KW. Research progress on the relationship between tea shape and quality (continued). TEA IN FUJIAN. (2000) 2000:2–4. doi: 10.3969/j.issn.1005-2291.2000.02.001

CrossRef Full Text | Google Scholar

17. Fan QY, Lan XP, Li WC, Wang JP, Tang XL. Effects of different processing techniques on the quality of green tea by Jingbai 2#. Chin Tea Process. (2017) 2017:54–6+71. doi: 10.15905/j.cnki.33-1157/ts.2017.z2.011

CrossRef Full Text | Google Scholar

18. Flaig M, Qi S, Wei GD, Yang XG, Schieberle P. Characterization of the key odorants in a high-grade chinese green tea beverage (Camellia sinensis; Jingshan cha) by means of the sensomics approach and elucidation of odorant changes in tea leaves caused by the tea manufacturing process. J Agric Food Chem. (2020) 68:5168–79. doi: 10.1021/acs.jafc.0c01300

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Yu XL. Effects of different withering methods on components metabolism related to color, aroma and taste quality in green tea. Hua Agric Univ. (2020) 133–65. doi: 10.27158/d.cnki.ghznu.2020.000041

CrossRef Full Text | Google Scholar

20. Zhu YM, Dong JJ, Jin J, Liu JH, Zheng XQ, Lu JL, et al. Roasting process shaping the chemical profile of roasted green tea and the association with aroma features. Food Chem. (2021) 353:129428. doi: 10.1016/j.foodchem.2021.129428

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Guo YL, Lai LL. Process optimization of rolled and curled green tea made of big leaves species. Chin J Tropic Crops. (2011) 32:2152–6. doi: 10.3969/j.issn.1000-2561.2011.11.030

CrossRef Full Text | Google Scholar

22. Guo L, Yu FL, Zhao F, Zhang Y, Zhu Y, Dai WD, et al. Correlation of physical properties of granular green tea with its appearance quality. Food Sci. (2020) 41:25–30. doi: 10.7506/spkx1002-6630-20190930-367

CrossRef Full Text | Google Scholar

23. Xue JJ, Liu PP, Guo GY, Wang WW, Zhang JY, Wang W, et al. Profiling of dynamic changes in non-volatile metabolites of shaken black tea during the manufacturing process using targeted and non-targeted metabolomics analysis. LWT. (2022) 156:113010. doi: 10.1016/j.lwt.2021.113010

CrossRef Full Text | Google Scholar

24. Chen S, Liu HH, Zhao XM, Li XL, Shan WN, Wang XX, et al. Non-targeted metabolomics analysis reveals dynamic changes of volatile and non-volatile metabolites during oolong tea manufacture. Food Res Intern. (2020) 128:108778. doi: 10.1016/j.foodres.2019.108778

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Zhou JT, Yu XL, He C, Qiu AD, Li YC, Shu QN, et al. Withering degree affects flavor and biological activity of black tea: a non-targeted metabolomics approach. LWT. (2020) 130:109535. doi: 10.1016/j.lwt.2020.109535

CrossRef Full Text | Google Scholar

26. Yang YQ, Zhang MM, Hua JJ, Deng YL, Jiang YW, Li J, et al. Quantitation of pyrazines in roasted green tea by infrared-assisted extraction coupled to headspace solid-phase microextraction in combination with GC-qqq-MS/MS. Food Res Intern. (2020) 134:109167. doi: 10.1016/j.foodres.2020.109167

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Zhang MM, Yang YQ, Yuan HB, Hua JJ, Deng YL, Jiang YW, et al. Contribution of addition theanine/sucrose on the formation of chestnut-like aroma of green tea. LWT. (2020) 129:109512. doi: 10.1016/j.lwt.2020.109512

CrossRef Full Text | Google Scholar

28. Zhu JC, Niu YW, Xiao ZB. Characterization of the key aroma compounds in Laoshan green teas by application of odour activity value (OAV), gas chromatography-mass spectrometry-olfactometry (GC-MS-O) and comprehensive two-dimensional gas chromatography mass spectrometry (GC × GC-qMS). Food Chem. (2021) 339:128136. doi: 10.1016/j.foodchem.2020.128136

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Wang HJ, Hua JJ, Yu QY, Li J, Wang JJ, Deng YL, et al. Widely targeted metabolomic analysis reveals dynamic changes in non-volatile and volatile metabolites during green tea processing. Food Chem. (2021) 363:130131. doi: 10.1016/j.foodchem.2021.130131

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Yang ZY, Baldermann S, Watanabe N. Recent studies of the volatile compounds in tea. Food Res Intern. (2013) 53:585–99. doi: 10.1016/j.foodres.2013.02.011

CrossRef Full Text | Google Scholar

31. Ravichandran R. Carotenoid composition, distribution and degradation to flavour volatiles during black tea manufacture and the effect of carotenoid supplementation on tea quality and aroma. Food Chem. (2002) 78:23–8. doi: 10.1016/S0308-8146(01)00303-X

CrossRef Full Text | Google Scholar

32. Guo L, Zhao F, Chen JZ, Wang SQ, Yu FL, Zhang Y, et al. Analysis of aroma characteristic components of granular green tea. Tea Comm. (2021) 48:287–92. doi: 10.3969/j.issn.1009-525X.2021.02.014

CrossRef Full Text | Google Scholar

33. Zhong JY, Chen N, Huang SC, Fan XW, Zhang Y, Ren DB, et al. Chemical profiling and discrimination of green tea and Pu-erh raw tea based on UPLC–Q–Orbitrap–MS/MS and chemometrics. Food Chem. (2020) 326:126760. doi: 10.1016/j.foodchem.2020.126760

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Yang YQ, Chen JY, Jiang YW, Qian MC, Deng YL, Xie JL, et al. Aroma dynamic characteristics during the drying process of green tea by gas phase electronic nose and gas chromatography-ion mobility spectrometry. LWT. (2022) 154:112691. doi: 10.1016/j.lwt.2021.112691

CrossRef Full Text | Google Scholar

35. Guo XY, Ho CT, Schwab W, Wan XC. Aroma profiles of green tea made with fresh tea leaves plucked in summer. Food Chem. (2021) 363:130328. doi: 10.1016/j.foodchem.2021.130328

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Wang MQ, Ma WJ, Shi J, Zhu Y, Lin Z, Lv HP. Characterization of the key aroma compounds in Longjing tea using stir bar sorptive extraction (SBSE) combined with gas chromatography-mass spectrometry (GC–MS), gas chromatography-olfactometry (GC–O), odor activity value (OAV), and aroma recombination. Food Res Intern. (2020) 130:108908. doi: 10.1016/j.foodres.2019.108908

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Guo XY, Ho CT, Wan XC, Zhu H, Liu QO, Wen Z. Changes of volatile compounds and odor profiles in Wuyi rock tea during processing. Food Chem. (2021) 341:128230. doi: 10.1016/j.foodchem.2020.128230

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Zhu Y, Lv HP, Shao CY, Kang SY, Zhang Y, Guo L, et al. Identification of key odorants responsible for chestnut-like aroma quality of green teas. Food Res Intern. (2018) 108:74–82. doi: 10.1016/j.foodres.2018.03.026

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Liao XL, Yan JN, Wang B, Meng Q, Zhang LY, Tong HR. Identification of key odorants responsible for cooked corn-like aroma of green teas made by tea cultivar ‘Zhonghuang 1′. Food Res Intern. (2020) 136:109355. doi: 10.1016/j.foodres.2020.109355

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Cui JL, Katsuno T, Totsuka K, Ohnishi T, Takemoto H, Mase N, et al. Characteristic fluctuations in glycosidically bound volatiles during tea processing and identification of their unstable derivatives. J Agric Food Chem. (2016) 64:1151–7. doi: 10.1021/acs.jafc.5b05072

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: round green tea, volatile metabolites, second-drying temperature, fatty acid-derived volatiles, glycoside-derived volatiles, nontargeted metabolomic analysis

Citation: Wang H, Yu Y, Ouyang W, Jiang Y, Wang J, Hua J and Yuan H (2022) Assessment of Variations in Round Green Tea Volatile Metabolites During Manufacturing and Effect of Second-Drying Temperature via Nontargeted Metabolomic Analysis. Front. Nutr. 9:877132. doi: 10.3389/fnut.2022.877132

Received: 16 February 2022; Accepted: 21 March 2022;
Published: 14 April 2022.

Edited by:

Jinkai Zheng, Institute of Food Science and Technology (CAAS), China

Reviewed by:

Tao Feng, Shanghai Institute of Technology, China
Ya-Fang Shang, Hefei University of Technology, China

Copyright © 2022 Wang, Yu, Ouyang, Jiang, Wang, Hua and Yuan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jinjie Hua, huajinjie@tricaas.com; Haibo Yuan, 192168092@tricaas.com

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.