Plant-Climate Interaction Effects: Changes in the Relative Distribution and Concentration of the Volatile Tea Leaf Metabolome in 2014–2016

Climatic conditions affect the chemical composition of edible crops, which can impact flavor, nutrition and overall consumer preferences. To understand these effects, we sampled tea (Camellia sinensis (L.) Kuntze) grown in different environmental conditions. Using a target/nontarget data analysis approach, we detected 564 metabolites from tea grown at two elevations in spring and summer over 3 years in two major tea-producing areas of China. Principal component analysis and partial least squares-discriminant analysis show seasonal, elevational, and yearly differences in tea from Yunnan and Fujian provinces. Independent of location, higher concentrations of compounds with aromas characteristic of farmers’ perceptions of high-quality tea were found in spring and high elevation teas. Yunnan teas were distinct from Fujian teas, but the effects of elevation and season were different for the two locations. Elevation was the largest source of metabolite variation in Yunnan yet had no effect in Fujian. In contrast seasonal differences were strong in both locations. Importantly, the year-to-year variation in chemistry at both locations emphasizes the importance of doing multi-year studies, and further highlights the challenge farmers face when trying to produce teas with specific flavor/health (metabolite) profiles.

Introduction

Greenhouse gasses accumulating in the atmosphere cause variability in local climates including prolonged heatwaves, droughts, heavy rains, flooding, cold, and frost, all of which damage crop production (IPCC, 2014; USGCRP, 2017; Fussel et al., 2018). The effect of extreme weather and climate variability on yield has been extensively studied for many crops (Kang et al., 2009; Lobell and Gourdji, 2012; Kurukulasuriya and Rosenthal, 2013). Most studies focus on one climate variable at a time; for example, temperature increases of up to 4°C are expected to decrease rice (Asia), wheat (India), and maize (US and Africa) yields by 20%–60% (NRC, 2011). Equally important to climate effects on yields is understanding how plant–climate interactions affect crop quality, including the secondary metabolites that contribute to the sensory and nutritional properties of plant materials. Changes in crop quality have been shown to influence consumer preference, acceptance, and, ultimately demand (Ahmed and Stepp, 2016). As attribution science (the science of attributing specific adverse weather conditions to climate change) advances (Otto, 2016; Imada et al., 2018; Knutson et al., 2018), understanding how plant–climate interactions influences human and natural systems is critical in developing long-term, sustainable agro-ecosystems.

Tea (Camellia sinensis (L.) Kuntze) is a long-lived crop grown in subtropical regions that are vulnerable to climate variability and is predicted to be severely impacted by climate change (IPCC, 2014; Han et al., 2018). The arrival of the East Asian Monsoon rains, an extreme weather event, is occurring earlier, lasting longer and shortening the spring harvest for high quality tea (Ahmed et al., 2014; Kowalsick et al., 2014; Boehm et al., 2016). Since tea plants are mostly grown in mountainous areas, elevational differences can also affect tea quality: higher quality generally occurs at higher elevation (Han et al., 2017; Kfoury et al., 2018b). Therefore, it is of great interest to understand the effect of season and elevation over multiple years in order to develop strategies for sustaining high quality tea in the era of climate change (FAO, 2015; Han et al., 2017; Han et al., 2018). Our previous work on tea in Yunnan Province, China reveals striking changes in the distribution and concentration of metabolites in response to differences in precipitation and temperature within a 1-year period (Ahmed et al., 2014; Kowalsick et al., 2014; Robbat Jr et al., 2017; Kfoury et al., 2018a; Kfoury et al., 2018b).

Tea has a complex secondary metabolite profile comprised of hundreds of volatile and non-volatile compounds that contribute to quality by impacting flavor, nutritional, and health attributes (Ahmed and Stepp, 2012). Of the ∼450 volatile metabolites detected in any given tea sample, our previous research found that two-thirds increased or decreased in concentration; most by more than 50%, more than 100 by 100%, and some by 1,000% due to differences in spring and summer rainfall (∼0 vs. 400 mm) and temperature (22°C vs. 28°C) as well as elevational differences in temperature (5°C) (Kowalsick et al., 2014; Kfoury et al., 2018b). We found a greater number of metabolites with higher concentrations in the spring (cooler temperature/no rain) and high elevation (cooler temperature) teas that exhibited sweet, floral, honey-like aromas compared to hay, grassy, earthy notes in the summer (warmer temperature/monsoon rain) and low elevation (warmer temperature) teas. These results are consistent with farmer sensory perceptions of high- and low-quality tea (Ahmed and Stepp, 2012; Ahmed et al., 2014; Han et al., 2017). High elevation teas also contained more and higher concentrations of volatile metabolites that possess analgesic, antianxiety, antibacterial, anticancer, antidepressant, antifungal, anti-inflammatory, antioxidant, antistress, and cardioprotective properties (Kfoury et al., 2018b). This finding is extremely important since time and location of harvest can affect the medicinal properties of tea extracts used in clinical trials (Tsao et al., 2009; Wang et al., 2014; Yan et al., 2017; Balsan et al., 2019).

Less dramatic, although equally important, are the differences in catechin (polyphenolic antioxidants) and methylxanthine (stimulants) concentrations. Although these metabolites were higher in spring tea, the monsoon rains induced higher total phenolic concentrations and antioxidant potential in summer tea (Ahmed et al., 2014). Also, the low elevation tea contained higher concentrations of caffeine, epicatechin gallate, gallocatechin, and catechin, while other catechins and methylxanthines were indifferent to elevation (Kfoury et al., 2018b). Importantly, many tea famers experience up to a 50% decrease in revenue for summer and low elevation tea most likely due to an increase in unpleasant aromatics and compounds associated with bitterness and astringency (Ahmed et al., 2014; Han et al., 2017).

Clearly, both elevation and season have important effects on the chemistry and quality of tea. Whether these differences hold true across years is understudied. Here, we performed a multi-year study to explore the effects of these drivers on tea quality. Specifically, we study the variation of tea quality based on the differences in metabolite distribution and concentration in tea due to seasonal and elevation differences over a 3-year period in Yunnan and Fujian Provinces, China. To understand how changes in climate will affect plant volatiles, it is necessary to obtain the total, detectable metabolome. Here, we employed a targeted/untargeted data analysis approach to analyze gas chromatography/mass spectrometry (GC/MS) data of tea extracted by stir bar sorptive extraction.

Materials and Methods

Sample Collection

Tea leaves were collected from four communities located at high and low elevations in two major tea-producing provinces of China during the spring and summer harvest seasons from 2014–2016. This includes one high (1651 m) and one low (1162 m) elevation tea-producing community in Menghai County in Yunnan Province, which grow the large-leaf variety of tea (Camellia sinensis var. assamica). The spring sampling occurred in March, with summer sampling in June. The other tea-producing communities are in Anxi County in Fujian Province and included one high (690 m) and one low (112 m) elevation site, which grow the small-leaf variety (Camellia sinensis var. sinensis). The spring sampling occurred in May, with the summer sampling in July. Samples consisted of the terminal bud plus two adjacent leaves from five plants in three plots within three tea farms in each community, collected over three consecutive days. Since results of a previous study showed no significant plot-to-plot differences on the same farm (Ahmed et al., 2014), samples from each plot were pooled and homogenized to produce daily samples (n = 3) and treated as independent replicates within each sampling period at each site. A microwave oven was used in the field to stop enzymatic oxidation (Ahmed et al., 2014; Kowalsick et al., 2014). The dried leaves were shipped to Tufts University, where they were stored at −20°C until analyzed.

Sample Preparation

Aqueous infusions were prepared by brewing 3 g of tea in 30 ml of deionized water at 90°C and cooling to room temperature for 30 min in a closed container. 10ml aliquots were syringe filtered (0.45-µm polytetrafluoroethylene, Fisher Scientific, Pittsburgh, PA) into 10ml Teflon-sealed vials. Organics sorbed into a 0.5-mm thick × 10-mm long polydimethylsiloxane stir bar (Gerstel, Mülheiman der Ruhr, Germany) stirred at 1,200 rpm for 1 h. Stir bars were removed from the vials, rinsed with deionized water, dried with a lint-free wipe, and placed into glass desorption tubes for analysis.

GC/MS Analysis

Analyses were performed on an Agilent (Santa Clara, CA) 6890/5975 GC/MS equipped with a MultiPurpose Sampler (Gerstel). The thermal desorption unit (TDU, Gerstel) provided splitless transfer of the sample from the stir bar into a programmable temperature vaporization inlet (CIS, Gerstel). The TDU was heated from 40°C (0.70 min) to 275°C (3 min) at 600°C/min under 50 ml/min of helium. After 0.1 min, the CIS, operating in solvent vent mode, was heated from −100°C to 275°C (5 min) at 12°C/s. The GC column (30 m × 250 µm × 0.25 µm RXI-5MS, Restek, Bellefonte, PA) was heated from 40°C (1 min) to 280°C at 5°C/min with 1.2 ml/min of constant helium flow. The MS operated in full scan mode between 40 and 350 m/z, with an electron impact ionization energy of 70 eV. The ion source and quadrupole temperatures were 230°C and 150°C, respectively. A standard mixture of C₇-C₃₀n-alkanes (Sigma-Aldrich, St. Louis, MO) was used to calculate the retention index of each compound. Concentration differences were calculated as relative peak area differences for each compound using naphthalene-d⁸ (Restek) as an internal standard. A total of 300 reference standards were purchased from Sigma-Aldrich, Fisher Scientific, Alfa Aesar (Ward Hill, MA), TCI (Tokyo, Japan), Acros Organics (Pittsburgh, PA) and MP Biomedicals (Santa Ana, CA) to confirm compound identity.

Data Analysis

Ion Analytics (Gerstel) data analysis software was used to analyze the samples based on a target/nontarget approach as previously described in (Robbat Jr et al., 2017; Kfoury et al., 2018a).

Statistical Analysis

Principal component analysis (PCA) and partial least squares-discriminant analysis (PLS-DA) were performed on autoscaled (mean-centered and unit-variance scaled) data using MetaboAnalyst 4.0 (Chong et al., 2018). Permutational multivariate analysis of variance (PERMANOVA) was conducted using 999 permutations using the vegan package in R (R, 2014; Oksanen et al., 2018). PCA was used to visualize group differences with confirmation made by PERMANOVA. Because current implementations of PLS-DA cannot account for factorial experimental designs, PLS-DA was used only to identify important metabolites for differentiating levels of one variable (e.g., elevation) across all levels of other variables (e.g., across all years and seasons). The quality of the PLS-DA model is described by R² and Q² values (Eriksson et al., 2006). R² measures the degree of fit of data to the model. A seven-fold cross validation was used to produce Q², which measures the predictive power of the model. A feature of PLS-DA is the ability to summarize predictor variable importance across all predictive components with a variable importance in projection (VIP) score (Chong and Jun, 2005). Metabolites with a VIP > 1.0 and statistically different among levels of a response variable (Mann-Whitney test, p < 0.05) were considered important discriminators. Metabolite diversity was calculated as a Simpson’s diversity index (1-D) and the effects of season, elevation year, and their interactions on metabolite diversity were assessed using an ANOVA (p < 0.05).

Climate Data

Daily maximum and minimum temperatures were obtained from the Climate Prediction Center’s Global Daily Temperature dataset with each grid cell spanning 0.5° latitude × 0.5° longitude (CPC, 2018). Daily precipitation totals were obtained from the Integrated Multi-satellitE Retrievals (IMERG) for Global Precipitation Measurement (GPM) dataset (Huffman et al., 2014), with a spatial resolution of 0.1° latitude × 0.1° longitude. Values are derived from passive microwave sensors from the GPM constellation and ground-based rain gauge observations. Since the IMERG data were not available prior to March 12, 2014, some spring 2014 days are not included in the analysis. The coordinates for each site and corresponding data were imported into a Geographic Information System. Temperature and precipitation values at each site were extracted from corresponding grid cells.

Results

Yunnan teas, including the ones used in this study, were used to create a Yunnan-specific database of ∼600 compounds. The database served as the initial list of target compounds when we analyzed tea from Fujian based on our targeted/untargeted workflow. For example, spectral deconvolution of the total ion current chromatogram (TIC, Figure 1A) of high elevation, spring tea yielded 444 target compounds (Figure 1B), where each colored peak indicates a single compound. Subtraction of the mass spectra for these target compounds revealed 32 compounds (nontargets) specific to this sample (Figure 1C). By adding these compounds to the database, they become target compounds for the next sample. Subtraction of their mass spectra in subsequent samples reveals new metabolites that are specific to that sample. Based on this approach, we detected 518 volatile metabolites, 58 of them unique to plants grown in Fujian from 2014–2016. Similarly, 506 metabolites were detected, 46 unique, in plants grown in Yunnan over the same time. It should be noted that none of the compounds detected, including unknowns, are known pesticides in the NIST and Wiley databases.

FIGURE 1

Figure 1 Total ion current (TIC) chromatogram of high elevation spring tea from Fujian (A) and reconstructed ion current (RIC) chromatograms of target (B) and non-target (C) compounds. Each colored peak in the RIC corresponds to a specific compound in the sample.

Although Yunnan and Fujian plants produced 460 common compounds (∼90%), differences in their concentration as well as unique metabolites yielded a PCA score plot, Figure 2, that separated the samples by location (and/or variety), which we expected due to differences in terroir and farmer practices. Interestingly, the monoterpenes are negatively correlated with PC1 and, therefore, more concentrated in Yunnan tea. Whereas, the sesquiterpenes correlate positively with PC1 and are higher in concentration in Fujian tea. Because of the inherent differences in tea leaf chemistry, we analyzed the data from each province separately to understand the effect of climate on plant metabolites.

FIGURE 2

Figure 2 Principal component analysis (PCA) score plot of Yunnan and Fujian tea.

Plant-Climate Interactions: Yunnan Tea

Table 1 lists the 10-day cumulative rainfall and average temperature prior to each harvest from 2014–2016. We selected this time period based on previous studies, where differences in metabolite chemistry were observed 5 days after the onset of the East Asian Monsoon rains (Ahmed et al., 2014; Kowalsick et al., 2014). Spring rainfall was ≤ 0.1 mm and summer 58 ± 14 mm, with the seasonal difference in temperature ∼3.5 ± 1°C. Although the two elevational sites, from which the samples were harvested, fall within the same latitude/longitude grid cell for temperature, we estimated a 3°C difference between the high (cooler) and low (warmer) elevations using the adiabatic lapse rate (NOAA, 2017).

TABLE 1

Table 1 10-day cumulative rainfall and average temperature prior to each harvest (Huffman et al., 2014; CPC, 2018).

Figure 3 shows the PCA score plots of PC1 vs. PC2 (A) and PC1 vs. PC3 (B), which explain 43% of the variation in metabolite chemistry. PCA separated the metabolite profiles by elevation (circles vs. triangles) on PC1 (except for 2014), season (open vs. closed shapes) on PC2, and year, 2014 (red) from 2015/2016 (black/blue), on PC3. The separations were confirmed by three-way PERMANOVA showing a significant main effect of elevation (F_(1,24) = 34.57, p = 0.001), season (F_(1,24) = 11.23, p = 0.001), and year (F_(2,24) = 13.05, p = 0.001). PERMANOVA also revealed significant interactive effects between year and elevation (F_(2,24) = 13.29, p = 0.001) and year and season (F_(2,24) = 4.01, p = 0.010), but not season and elevation (F_(2,24) = 0.78, p = 0.425). These interactive effects indicate that there is a different elevational and seasonal effect in at least one of the 3 years. As seen in the score plot (Figure 3A), 2014 samples do not separate by elevation along PC1, unlike 2015 and 2016 samples. In addition, 2015 samples separate differently by season than 2014/2016 samples along PC2.

FIGURE 3

Figure 3 Principal component analysis (PCA) score plots of Yunnan tea. (A) PC1 vs. PC2. (B) PC1 vs. PC3.

Metabolite diversity was described by the Simpson’s diversity index (Simpson, 1949), which accounts for the number of compounds and their relative abundances in the samples. An index value of 0 means complete dominance of one compound and 1, complete evenness of all compounds. As the number and evenness of compounds increases, namely more compounds with similar relative peak areas, so does the diversity index. Metabolite diversity index ranged from 0.81 to 0.95 and was significantly affected by elevation (F_(1,24) = 10.85, p = 0.003), season (F_(1,24) = 30.42, p < 0.001), and year (F_(2,24) = 19.60, p < 0.001). There were significant interactions between season and year (F_(2,24) = 7.42, p = 0.003) and between elevation and year (F_(2,24) = 4.94, p = 0.016), but not between season and elevation (F_(2,24) = 0.27, p = 0.609). These results match the PERMANOVA for the metabolite relative peak area data. In general, metabolite diversity was greatest in low elevation sites and summer, meaning these plants produced a greater number of compounds having similar metabolite concentrations in response to the conditions studied.

We used PLS-DA to identify metabolites that largely vary between elevations, seasons, and years, independently. The models for elevation (R² = 0.867, Q² = 0.659), season (R² = 0.884, Q² = 0.736), and year (R² = 0.854, Q² = 0.782) had high explanatory and predictive power. We identified 138 metabolites that discriminated high from low elevation (Supplementary Table S1), 129 metabolites that differentiated spring from summer (Supplementary Table S2), and 155 metabolites that separated years (Supplementary Table S3), where the main separation effect is between 2014 and 2015/2016. Some metabolites are affected by more than one climate variable resulting in 353 discriminator compounds.

Plant-Climate Interactions: Fujian Tea

Table 1 also lists the climate data for Fujian Province. In contrast to Yunnan, a larger temperature increase occurs from spring to summer (∼6.1 ± 1°C), but rainfall patterns were erratic from year to year. The high and low elevation sites are ∼32 km apart from each other and can therefore be distinguished by latitude/longitude climate grids.

Figure 4 shows the score plots of PC1 vs. PC2 (A) and PC1 vs. PC3 (B). Like Yunnan, the first three axes explain 43% of the variation in the data. The score plot revealed metabolite profiles separated by year (red, black, and blue) on PC2, but no clear seasonal or elevational separation was found. The three-way PERMANOVA confirmed the yearly separation (F_(2,24) = 7.61, p = 0.001) and lack of elevational separation (F_(1,24) = 1.96, p = 0.134), but revealed a significant seasonal separation (F_(1,24) = 3.79, p = 0.024). Upon further inspection, a seasonal separation (open vs. closed shapes) along PC4 was found (Supplementary Figure S1). There was also a significant interaction between year and season (F = 4.345, df = 2, p = 0.004) but not year and elevation (F = 1.686, df = 2, p = 0.144) or season and elevation (F = 1.470, df = 1, p = 0.205).

FIGURE 4

Figure 4 Principal component analysis (PCA) score plots of Fujian tea. (A) PC1 vs. PC2. (B) PC1 vs. PC3.

Metabolite diversity in Fujian varied significantly by year (F_(2,24) = 5.10, p = 0.014) and season (F_(1,24) = 4.32, p = 0.049), but not elevation (F_(1,24) = 0.003, p = 0.953). None of the two-way or three-way interactions were significant (all p > 0.1). Metabolite diversity was highest in 2014, which had the largest temperature difference between spring and summer (∼8°C), and higher in summer compared to spring.

PLS-DA was used to identify metabolites that largely vary between seasons and years, independently. The models for season (R² = 0.913, Q² = 0.711) and year (R² = 0.901, Q² = 0.817) had high explanatory and predictive power. We identified 101 metabolites that differentiated spring from summer (Supplementary Table S4), and 133 metabolites that separated years (Supplementary Table S5). Of the 518 metabolites detected in Fujian teas, 226 discriminate the seasonal and/or yearly differences.

Discussion

Elevational Variation

While the elevation difference (∼500 m) is the same in the two provinces, it was the largest source of metabolite variation in Yunnan yet had no effect in Fujian. The Yunnan sites are ∼1,000 m higher in elevation, meaning the Fujian sites are ∼6°C warmer than even the low elevation site in Yunnan (NOAA, 2017). The much warmer temperatures at the Fujian sites could account for the lack of elevational variation seen. In addition, while we classified the site located at 690 m in Fujian as high elevation based on being 500 m higher than the other Fujian site and the elevations found in the area, an elevation of 690 m is not generally regarded as high elevation within the tea industry for high-quality tea.

In Yunnan, four compounds (43, 165, 168, 193—all unknowns) were unique to high elevation samples and four compounds (75, verbenone, 4-methyldecane, 4-ethylbenzaldehyde) were unique to low elevation samples. These compounds could be important biomarkers for distinguishing high and low elevation teas and if the unknowns prove to be of sensory or health importance, identification is necessary. In addition, high elevation teas contain greater concentrations of compounds that have been described as sweet, floral, and fruity (Perflavory, 2015), such as Z-jasmone, maltol, methyl 4-methylbenzoate, and bergamal. In contrast, the low elevation teas contain higher concentrations of green, earthy, waxy, and camphoraceous compounds such as 2-phenyl-2-propanol, nonanal, dodecanal, and verbenone (Perflavory, 2015). These results correspond with farmer perceptions that high elevation teas are generally of higher quality due to the sweet, floral, honey-like, and fruity aromatics (Owuor et al., 1990; Ahmed and Stepp, 2012; Han et al., 2017).

In contrast to previous work on Yunnan tea where only discriminating compounds with higher concentrations in high elevation tea had reported health benefits (Kfoury et al., 2018b), in this study compounds with positive correlations to both high and low elevation teas have reported health-beneficial properties (Supplementary Table S1). Nutraceutical compounds that distinguish high from low elevation tea include (E)-caryophyllene (anticancer, antidepressant, anti-inflammatory), isoeugenol (antibacterial, antioxidant), epi-α-cadinol (antibacterial, anti-inflammatory), (E)-nerolidol (antianxiety, antimalarial, anticancer), and α-pinene (antiviral, analgesic) (Astani et al., 2010; Tung et al., 2010; Russo, 2011; Guimarães et al., 2013; Bahi et al., 2014; Fidyt et al., 2016; Guerrini et al., 2016; Marteau et al., 2016; Shan et al., 2016; Resende et al., 2017). Health beneficial compounds that differentiate low from how elevation tea exhibit antibacterial (verbenone, decanal, undecanal, dodecanal), antifungal (nonanal), and antiseptic (2-phenoxyethanol) properties (Buhrer et al., 2002; Battinelli et al., 2006; Fujita et al., 2015; Iscan, 2017). Although the volatile metabolites represent a small fraction of the total tea mass, findings indicate that volatile tea extracts have proven health benefits (Takahashi et al., 2014; Li et al., 2017). However, further studies are needed to determine if the compounds found here are present in adequate concentrations to provide these purported health benefits.

Seasonal Variation

Although rainfall and temperature generally increase from spring to summer in both Yunnan and Fujian, Yunnan experiences a striking increase in rain and smaller changes in temperature compared to Fujian, which has a larger temperature increase and erratic rainfall from year to year. The leaves from both provinces had significant chemical variations due to seasonal changes. Interestingly, nine compounds (N-ethylsuccinimide, isomenthone, (E)-β-ocimene, menthone, 160, flouranthene, 1-ethyl-1H-pyrrole, 76, 1-ethyl-1H-pyrrole-2-carboxaldehyde) were discriminators of spring tea and five compounds (36, 167, 181, 224, 2-phenoxyethanol) discriminate summer tea in both Yunnan and Fujian (Supplementary Tables S2 and S4). In addition, 1-ethyl-1H-pyrrole was unique to spring tea in both locations and could be a key spring tea biomarker, independent of location, and variety.

Regardless of the inherent differences between the provinces, more compounds that have been described as sweet, floral, honey-like, and fruity are discriminators of spring tea (Perflavory, 2015). These include (E)-β-ocimene, 2,2,6-trimethylcyclohexanone, (Z)-jasmone, phenylethyl acetate, and 5-methylfurfural. On the other hand, 1-octen-3-ol, (2E,4E)-heptadienal, 2-pentylfuran, α-copaene, and cis-calamenene are higher in concentration in summer tea and are characterized as green, herbal, earthy, and woody (Supplementary Tables S2 and S4). In addition, compounds significant to both seasons have reported health benefits (Supplementary Tables S2 and S4). Spring-associated compounds include menthone (antibacterial, anti-inflammatory), eucalyptol (antibacterial, cardioprotective), indole (antibacterial, antifungal), and coumarin (anti-inflammatory, anticancer) (Egan et al., 1990; Muroi and Kubo, 1993; Nidiry and Babu, 2005; Santos et al., 2011; Shan et al., 2016; Iscan, 2017). Summer-associated compounds include (E)-β-ionone (anticancer, antibacterial), borneol (anti-inflammatory, analgesic), methyl anthranilate (antifungal), and quinoline (antimalarial, anticonvulsant, anticancer) (Muroi and Kubo, 1993; Nidiry and Babu, 2005; Oz et al., 2015; Anjali and Singh, 2016; Furtado Kelly et al., 2016).

Interestingly, in Yunnan, many of the discriminating metabolites are in both spring and high elevation or in the summer and low elevation teas. For example, spring and high elevation teas contain significantly higher amounts of 2-hydroxy-5-methylacetophenone, isoeugenol, 4-methylbenzaldehyde, and norfuraneol that are described as sweet, fruity, floral, honey-like compounds (Perflavory, 2015). In contrast, summer and low elevation teas have significantly higher amounts of 2,6-dimethyl-3,7-octadiene-2,6-diol, 2-phenoxyethanol, octanal, and nonanal that are described as herbal, green, fatty, and metallic (Perflavory, 2015). These findings agree with farmers’ perceptions that high elevation and spring teas are higher in aromatic quality possessing sweet, fruity, floral, honey-like characteristics (Ahmed and Stepp, 2012; Ahmed et al., 2014; Han et al., 2017). Most likely, these metabolites are induced by temperature since spring/high elevations are cooler than summer/low elevations.

Yearly Variation and Interactive Effects

In Yunnan, both PCA and PLS-DA show the yearly separation is between 2014 and 2015/2016. Based on the climate data, 2014 experienced the coolest spring and warmest summer temperatures compared to 2015 and 2016. Interestingly, many of the higher concentration compounds in the 2014 tea were also higher in concentration in either spring or summer tea. For example, compounds such as (E)-β-ocimene, (Z)-herboxide, (3Z)-hexenyl acetate, and safranal were higher in concentration in 2014 and spring teas and are described as sweet, herbal, fruity aromas (Perflavory, 2015). Compounds such as (4Z)-heptenal, pyridine, and 2-methyldecane were higher in concentration in 2014 and summer teas and are characterized as fatty, green, and fishy (Perflavory, 2015). These compounds are likely being induced by the cooler spring or warmer summer temperatures. In contrast, compounds higher in concentration in 2015/2016 are not affected by seasonal variations such as 4-ethyl-2-methoxyphenol (smoky, phenolic), catechol (no aroma), 4-methyl-3-penten-2-one (sweet, earthy), and γ-heptalactone (sweet, nutty) (Perflavory, 2015), which are likely not influenced by temperature changes. Another indication that 2014 differed from 2015/2016 is that 2014 samples do not separate by elevation, which is confirmed by the interactive effect between year and elevation. Compounds such as hexanoic acid, 2-nonanone, decanal, (E)-herboxide, 4-methylbenzaldehyde, and 2-methoxy-4-vinylphenol exhibit no change in concentration between elevations in 2014, but are significantly higher in concentration at one elevation or the other in 2015/2016 tea.

Also, 2014 and 2016 samples separate similarly by season whereas 2015 samples do not follow the same pattern along PC2, which is confirmed by the interactive effect between year and season. While several compounds exhibit the opposite change in concentration from spring to summer, others show a more or less enhanced concentration change in 2015 compared to 2014/2016. For example, (E)-β-ocimene, methyl salicylate, and 7-methoxycoumarin exhibit opposite concentration changes from spring to summer in 2015 compared to 2014/2016. In addition, hexanol, isomenthone, and camphor have a greater concentration difference between spring and summer whereas 2E-hexenal, benzenacetonitrile and 4-methyl-3-penten-2-one have a smaller difference in 2015. The climate data do not provide any indication as to what is causing these interactive effects. However, (E)-β-ocimene, methyl salicylate, 2E-hexenal, and 4-methyl-3-penten-2-one are induced by herbivory in tea, so these year-to-year differences could reflect differences in the timing of pest outbreaks (Cai et al., 2014; Kfoury et al., 2017; Scott et al., 2019).

Year-to-year variations account for the greatest source of variation in Fujian and is likely due to the inconsistent rainfall patterns seen from year-to-year (Table 1). The significant interactive effect between year and season confirms this finding, not found between year and elevation or season and elevation. In 2016, the plants experienced an extremely dry summer, which is the reverse of the previous two years. As a result, many metabolite concentrations increased/decreased in the opposite direction from 2014 and 2015. For example, safranal (sweet, herbal), norfuraneol (sweet, caramel), cyclohexanone (minty), o-xylene (geranium), isoeugenol (floral, clove), and geranyl acetone (floral, fruity) are higher in concentration in spring tea in 2014 and 2015, but higher in summer tea in 2016. Similarly, 2-ethylhexanol (green, oily), camphor (camphor, medicinal), methyl hexanoate (fruity, fatty), biphenyl (floral, green), and 5-ethyl-2(5H)-furanone (no aroma) are higher in concentration in 2014/2015 summer, compared to the 2016 spring tea. It is likely that the changes in concentration of these and many other metabolites that behave similarly are the result of changes in rainfall.

Conclusion

In this work, we demonstrated that our target/nontarget approach provides efficient and comprehensive analysis of a complex, natural product, namely, tea. Findings show the differences in season, elevation and year cause significant and interactive alterations in tea chemistry. Independent of location, cooler temperatures in the spring and at high elevation concomitant with lower rainfall, results in higher concentrations of compounds with aromas characteristic consistent with farmers’ perceptions of high-quality tea. Although many of the metabolites identified have reported health-beneficial properties studies are needed to assess efficacy and identify unknowns that are discriminator compounds that might also contribute to quality. Given the interactive effects found between years and season/elevation, future studies should be cautious in drawing conclusions based on only 1 year of sampling and more multi-year studies should be conducted to assist farmers with the challenge of producing a consistently high quality product. More natural-human systems investigations are needed to assess plant response under ever changing environmental conditions.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Author Contributions

AR, CO, SA, SC, TG, and JS acquired funding and conceptualized the project. DX, CL, SA, and WH supported with in-field sample collections and logistics. CM acquired and analyzed the climate data. NK acquired and analyzed the metabolite data. NK and ES performed statistical analysis. NK, AR, ES, and CO wrote the manuscript draft. All authors contributed to manuscript revision and approved the submitted version.

Funding

This work was supported by the National Science Foundation, Grant BCS-1313775.

Conflict of Interest

AR is the founder of Ion Analytics.

The remaining 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.

Acknowledgments

The authors appreciate the support of Gerstel GmbH, Gerstel USA, Agilent Technologies, and Ion Analytics for providing the instruments and software to analyze samples.

Supplementary Material

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

References

Ahmed, S., Stepp, J. R. (2012). “Green Tea: The Plants, Processing, Manufacturing and Production,” in Tea in Health and Disease Prevention. Ed. Preedy, V. (London, UK: Academic Press).

Google Scholar

Ahmed, S., Stepp, J. R. (2016). Beyond yields: Climate change effects on specialty crop quality and agroecological management. Elem. Sci. Anth. 4, 000092. doi: 10.12952/journal.elementa.000092

CrossRef Full Text | Google Scholar

Ahmed, S., Stepp, J. R., Orians, C., Griffin, T., Matyas, C., Robbat, A., Jr., et al. (2014). Effects of extreme climate events on tea (Camellia sinensis) functional quality validate indigenous farmer knowledge and sensory preferences in tropical China. PloS One 9 (10), e109126. doi: 10.1371/journal.pone.0109126

PubMed Abstract | CrossRef Full Text | Google Scholar

Anjali, P. D., Singh, D. (2016). Quinoline: A diverse therapeutic agent. Int. J. Pharm. Sci. Res. 7, 1–13. doi: 10.13040/IJPSR.0975-8232.7(1).1-13

CrossRef Full Text | Google Scholar

Astani, A., Reichling, J., Schnitzler, P. (2010). Comparative study on the antiviral activity of selected monoterpenes derived from essential oils. Phytother. Res. 24 (5), 673–679. doi: 10.1002/ptr.2955

PubMed Abstract | CrossRef Full Text | Google Scholar

Bahi, A., Al Mansouri, S., Al Memari, E., Al Ameri, M., Nurulain, S. M., Ojha, S. (2014). beta-Caryophyllene, a CB2 receptor agonist produces multiple behavioral changes relevant to anxiety and depression in mice. Physiol. Behav. 135, 119–124. doi: 10.1016/j.physbeh.2014.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Balsan, G., Pellanda, L. C., Sausen, G., Galarraga, T., Zaffari, D., Pontin, B., et al. (2019). Effect of yerba mate and green tea on paraoxonase and leptin levels in patients affected by overweight or obesity and dyslipidemia: a randomized clinical trial. Nutr. J. 18 (1), 5. doi: 10.1186/s12937-018-0426-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Battinelli, L., Daniele, C., Cristiani, M., Bisignano, G., Saija, A., Mazzanti, G. (2006). In vitro antifungal and anti-elastase activity of some aliphatic aldehydes from Olea europaea L. fruit. Phytomedicine 13, 558–563. doi: 10.1016/j.phymed.2005.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Boehm, R., Cash, S., Anderson, B., Ahmed, S., Griffin, T., Robbat, A., et al. (2016). Association between Empirically Estimated Monsoon Dynamics and Other Weather Factors and Historical Tea Yields in China: Results from a Yield Response Model. Climate 4 (2), 20. doi: 10.3390/cli4020020

CrossRef Full Text | Google Scholar

Buhrer, C., Bahr, S., Siebert, J., Wettstein, R., Geffers, C., Obladen, M. (2002). Use of 2% 2-phenoxyethanol and 0.1% octenidine as antiseptic in premature newborn infants of 23-26 weeks gestation. J. Hosp. Infect. 51, 305–307. doi: 10.1053/jhin.2002.1249

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, X.-M., Sun, X.-L., Dong, W.-X., Wang, G.-C., Chen, Z.-M. (2014). Herbivore species, infestation time, and herbivore density affect induced volatiles in tea plants. Chemoecology 24 (1), 1–14. doi: 10.1007/s00049-013-0141-2

CrossRef Full Text | Google Scholar

Chong, I.-G., Jun, C.-H. (2005). Performance of some variable selection methods when multicollinearity is present. Chemometr. Intell. Lab. Syst. 78 (1-2), 102–112. doi: 10.1016/j.chemolab.2004.12.011

CrossRef Full Text | Google Scholar

Chong, J., Soufan, O., Li, C., Caraus, I., Li, S., Bourque, G., et al. (2018). MetaboAnalyst 4.0: towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. gky310 46, W486-W494. doi: 10.1093/nar/gky310

PubMed Abstract | CrossRef Full Text | Google Scholar

CPC (2018). CPC Global Temperature data provided by the NOAA/OAR/ESRl PSD, Available: https://www.esrl.noaa.gov/psd/data/gridded/data.cpc.globaltemp.html.

Google Scholar

Egan, D., O’Kennedy, R., Moran, E., Cox, D., Prosser, E., Thornes, R. D. (1990). The Pharmacology, Metabolism, Analysis, and Applications of Coumarin and Coumarin-Related Compounds. Drug Metab. Rev. 22, 503–529. doi: 10.3109/03602539008991449

PubMed Abstract | CrossRef Full Text | Google Scholar

Eriksson, L., Johansson, E., Kettaneh-Wold, N., Trygg, J. (2006). Multi- and Megavariate Data Analysis: Part II: Advanced Applications and Method Extensions. Malmo, Sweden: Umetric Inc.

Google Scholar

FAO (2015). Socio-economic implications of climate change for tea producing countries. Rome: Food and Agriculture Organization of the United Nations.

Google Scholar

Fidyt, K., Fiedorowicz, A., Strzadala, L., Szumny, A. (2016). beta-caryophyllene and beta-caryophyllene oxide-natural compounds of anticancer and analgesic properties. Cancer Med. 5 (10), 3007–3017. doi: 10.1002/cam4.816

PubMed Abstract | CrossRef Full Text | Google Scholar

Fujita, K.-i., Chavasiri, W., Kubo, I. (2015). Anti-Salmonella Activity of Volatile Compounds of Vietnam Coriander. Phytother. Res. 29 (7), 1081–1087. doi: 10.1002/ptr.5351

PubMed Abstract | CrossRef Full Text | Google Scholar

Furtado Kelly, S., de Oliveira Andrade, F., Campos, A., Rosim Mariana, P., Vargas-Mendez, E., Henriques, A. (2016). β-ionone modulates the expression of miRNAs and genes involved in the metastatic phenotype of microdissected persistent preneoplastic lesions in rats submitted to hepatocarcinogenesis. Mol. Carcinog. 56 (1), 184–196. doi: 10.1002/mc.22483

PubMed Abstract | CrossRef Full Text | Google Scholar

Fussel, H.-M., Lourenco, T. C., Hilden, M., Leitner, M., Marx, A., Prutsch, A., et al. (2018). National climate change vulnerability and risk assessments in Europe. Copenhagen, Denmark:European Environment Agency.

Google Scholar

Guerrini, A., Sacchetti, G., Grandini, A., Spagnoletti, A., Asanza, M., Scalvenzi, L. (2016). Cytotoxic effect and TLC bioautography-guided approach to detect health properties of amazonian hedyosmum sprucei essential oil. Evid.-based Complement. Altern. Med. 8, 1–8. doi: 10.1155/2016/1638342

CrossRef Full Text | Google Scholar

Guimarães, A. G., Quintans, J. S. S., Quintans-Júnior, L. J. (2013). Monoterpenes with analgesic activity—a systematic review. Phytother. Res. 27 (1), 1–15. doi: 10.1002/ptr.4686

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, W.-Y., Huang, J.-G., Li, X., Li, Z.-X., Ahammed, G. J., Yan, P., et al. (2017). Altitudinal effects on the quality of green tea in east China: a climate change perspective. Eur. Food Res. Technol. 243, 323–330. doi: 10.1007/s00217-016-2746-5

CrossRef Full Text | Google Scholar

Han, W.-Y., Li, X., Ahammed, G. J. (2018). Stress Physiology of Tea in the Face of Climate Change. Singapore: Springer.

Google Scholar

Huffman, G., Bolvin, D., Braithwaite, D., Hsu, K., Joyce, R., Xie, P. (2014). Integrated Multi-satellitE Retrievals for GPM (IMERG), NASA/GSFC. Available: https://disc.gsfc.nasa.gov/.

Google Scholar

Imada, Y., Shiogama, H., Takahashi, C., Watanabe, M., Mori, M., Kamae, Y., et al. (2018). Climate Change Increased the Likelihood of the 2016 Heat Extremes in Asia. Bull. Am. Meteorol. Soc 99 (1), S97–S101. doi: 10.1175/bams-d-17-0109.1

CrossRef Full Text | Google Scholar

IPCC (2014). Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part A: Global and Sectoral Aaspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, United Kingdom and New York, NY:Cambridge University Press.

Google Scholar

Iscan, G. (2017). Antibacterial and Anticandidal Activities of Common Essential Oil Constituents. Rec. Nat. Prod. 11 (4), 374–388. doi: 10.1016/S2221-1691(12)60220-2

CrossRef Full Text | Google Scholar

Kang, Y., Khan, S., Ma, X. (2009). Climate change impacts on crop yield, crop water productivity and food security - A review. Prog. Nat. Sci. 19, 1665–1674. doi: 10.1016/j.pnsc.2009.08.001

CrossRef Full Text | Google Scholar

Kfoury, N., Scott, E., Orians, C., Robbat, A. (2017). Direct contact sorptive extraction: a robust method for sampling plant volatiles in the field. J. Agric. Food Chem. 65 (38), 8501–8509. doi: 10.1021/acs.jafc.7b02847

PubMed Abstract | CrossRef Full Text | Google Scholar

Kfoury, N., Baydakov, E., Gankin, Y., Robbat, A. (2018a). Differentiation of key biomarkers in tea infusions using a target/nontarget gas chromatography/mass spectrometry workflow. Food Res. Int. 113, 414–423. doi: 10.1016/j.foodres.2018.07.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Kfoury, N., Morimoto, J., Kern, A., Scott, E. R., Orians, C. M., Ahmed, S., et al. (2018b). Striking changes in tea metabolites due to elevational effects. Food Chem. 264, 334–341. doi: 10.1016/j.foodchem.2018.05.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Knutson, T. R., Kam, J., Zeng, F., Wittenberg, A. T. (2018). CMIP5 Model-based assessment of anthropogenic influence on record global warmth during 2016 . Bull. Am. Meteorol. Soc 99 (1), S11–S15. doi: 10.1175/BAMS-D-17-0104.1

CrossRef Full Text | Google Scholar

Kowalsick, A., Kfoury, N., Robbat, A., Jr., Ahmed, S., Orians, C., Griffin, T., et al. (2014). Metabolite profiling of Camellia sinensis by automated sequential, multidimensional gas chromatography/mass spectrometry reveals strong monsoon effects on tea constituents. J. Chromatgr. A 1370, 230–239. doi: 10.1016/j.chroma.2014.10.058

CrossRef Full Text | Google Scholar

Kurukulasuriya, P., Rosenthal, S. (2013). “Climate change and agriculture: a review of impacts and adaptations,” in Climate Change series (Washington, DC:World Bank).

Google Scholar

Li, Y.-S., Kawasaki, Y., Tomita, I., Kawai, K. (2017). Antioxidant properties of green tea aroma in mice. J. Clin. Biochem. Nutr. 61 (1), 14–17. doi: 10.3164/jcbn.16-80

PubMed Abstract | CrossRef Full Text | Google Scholar

Lobell, D. B., Gourdji, S. M. (2012). The Influence of Climate Change on Global Crop Productivity. Plant Physiol. 160 (4), 1686–1697. doi: 10.1104/pp.112.208298

PubMed Abstract | CrossRef Full Text | Google Scholar

Marteau, C., Guitard, R., Penverne, C., Favier, D., Nardello-Rataj, V., Aubry, J.-M. (2016).Boosting effect of ortho-propenyl substituent on the antioxidant activity of natural phenols. Food Chem. 196, 418–427. doi: 10.1016/j.foodchem.2015.09.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Muroi, H., Kubo, I. (1993). Combination effects of antibacterial compounds in green tea flavor against streptococcus-mutans. J. Agric. Food Chem. 41 (7), 1102–1105. doi: 10.1021/jf00031a017

CrossRef Full Text | Google Scholar

Nidiry, E. S., Babu, C. S. B. (2005). Antifungal activity of tuberose absolute and some of its constituents. Phytother. Res. 19 (5), 447–449. doi: 10.1002/ptr.1630

PubMed Abstract | CrossRef Full Text | Google Scholar

NOAA (2017). Environmental Temperature Lapse Rate. NOAA Storm Prediction Center: National Weather Service. Available: http://www.spc.noaa.gov/exper/soundings/help/lapse.html.

Google Scholar

NRC (2011). Climate Stabilization Targets: Emissions, Concentrations, and Impacts over Decades to Millennia. Washington, DC: The National Academies Press.

Google Scholar

Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D. (2018). Vegan: Community Ecology Package. R. https://cran.r-project.org, https://github.com/vegandevs/vegan.

Google Scholar

Otto, F. E. L. (2016). The art of attribution. Nat. Clim. Change 6, 342. doi: 10.1038/nclimate2971

CrossRef Full Text | Google Scholar

Owuor, P. O., Obaga, S. O., Othieno, C. O. (1990). The Effects of Altitude on the Chemical Composition of Black Tea. J. Sci. Food Agric. 50 (1), 9–17. doi: 10.1002/jsfa.2740500103

CrossRef Full Text | Google Scholar

Oz, M., Lozon, Y., Sultan, A., Yang, K.-H. S., Galadari, S. (2015). Effects of monoterpenes on ion channels of excitable cells. Pharmacol. Ther. 152, 83–97. doi: 10.1016/j.pharmthera.2015.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Perflavory (2015). The Good Scents Company. Oak Creek, WI.

Google Scholar

R. (2014). R: A language and environment for statistical computing. Vienna, Austria: Foundation for Statistical Computing.

Google Scholar

Resende, D. B., Martins, H.H.D.A., Souza, T.B.D., Carvalho, D. T., Piccoli, R. H., Schwan, R. F., et al. (2017). Synthesis and in vitro evaluation of peracetyl and deacetyl glycosides of eugenol, isoeugenol and dihydroeugenol acting against food-contaminating bacteria. Food Chem. 237, 1025–1029. doi: 10.1016/j.foodchem.2017.06.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Robbat, A., Kfoury, N., Baydakov, E., Gankin, Y. (2017). Optimizing targeted/untargeted metabolomics by automating gas chrromatography/mass spectrometry workflows. J. Chromatgr. A 1505, 96–105. doi: 10.1016/j.chroma.2017.05.017

CrossRef Full Text | Google Scholar

Russo, E. B. (2011). Taming THC: potential cannabis synergy and phytocannabinoid-terpenoid entourage effects. Br. J. Pharmacol. 163 (7), 1344–1364. doi: 10.1111/j.1476-5381.2011.01238

PubMed Abstract | CrossRef Full Text | Google Scholar

Santos, M. R. V., Moreira, F. V., Fraga, B. P., Souza, D.P.D., Bonjardim, L. R., Quintans-Junior, L. J. (2011). Cardiovascular effects of monoterpenes: a review. Rev. Bras. Farmacogn. 21, 764–771. doi: 10.1590/S0102-695X2011005000119

CrossRef Full Text | Google Scholar

Scott, E. R., Li, X., Kfoury, N., Morimoto, J., Han, W.-Y., Ahmed, S., et al. (2019). Interactive effects of drought severity and simulated herbivory on tea (Camellia sinensis) volatile and non-volatile metabolites. Environ. Exper. Bot. 157, 283–292. doi: 10.1016/j.envexpbot.2018.10.025

CrossRef Full Text | Google Scholar

Shan, M. Q., Qian, Y., Yu, S., Guo, S. C., Zhang, L., Ding, A. W., et al. (2016). Anti-inflammatory effect of volatile oil from Schizonepeta tenuifolia on carrageenin-induced pleurisy in rats and its application to study of appropriate harvesting time coupled with multi-attribute comprehensive index method. J. Ethnopharmacol. 194, 580–586. doi: 10.1016/j.jep.2016.10.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Simpson, E. H. (1949).Measurement of Diversity. Nature 163 (4148), 688–688. doi: 10.1038/163688a0

CrossRef Full Text | Google Scholar

Takahashi, A., Watanabe, T., Fujita, T., Hasegawa, T., Saito, M., Suganuma, M. (2014). Green tea aroma fraction reduces β-amyloid peptide-induced toxicity in Caenhorhabditis elegans transfected with human β-amyloid minigene. Biosci. Biotechnol. Biochem. 78 (7), 1206–1211. doi: 10.1080/09168451.2014.921553

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsao, A. S., Liu, D., Martin, J., Tang, X., Lee, J. J., El-Naggar, A. K., et al. (2009). Phase II randomized, placebo-controlled trial of green tea extract in patients with high-risk oral premalignant lesions. Cancer Prev. Res. 2 (11), 931–941. doi: 10.1158/1940-6207.CAPR-09-012

CrossRef Full Text | Google Scholar

Tung, Y. T., Yen, P. L., Lin, C. Y., Chang, S. T. (2010). Anti-inflammatory activities of essential oils and their constituents from different provenances of indigenous cinnamon (Cinnamomum osmophloeum) leaves. Pharm. Biol. 48 (10), 1130–1136. doi: 10.3109/13880200903527728

PubMed Abstract | CrossRef Full Text | Google Scholar

USGCRP (2017). Climate Science Special Report: Fourth National Climate Assessment, Volume 1. (eds.) Wuebbles, DJ, Fahey, DW, Hibbard, KA, Dokken, DJ, Stewart, BC, Maycock, TK Washington, DC:U.S. Global Change Research Program.

Google Scholar

Wang, X., Tian, J., Jiang, J., Li, L., Ying, X., Tian, H., et al. (2014). Effects of green tea or green tea extract on insulin sensitivity and glycaemic control in populations at risk of type 2 diabetes mellitus: a systematic review and meta-analysis of randomised controlled trials. J. Hum. Nutr. Diet. 27 (5), 501–512. doi: 10.1111/jhn.12181

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, Y., Sui, X., Yao, B., Lavie, C. J., Blair, S. N. (2017). Is There a Dose–Response Relationship between Tea Consumption and All-Cause, CVD, and Cancer Mortality? J. Am. Coll. Nutr. 36 (4), 281–286. doi: 10.1080/07315724.2016.1261054

PubMed Abstract | CrossRef Full Text | Google Scholar