Next Article in Journal
GLP-1 Receptor Agonist Improves Mitochondrial Energy Status and Attenuates Nephrotoxicity In Vivo and In Vitro
Next Article in Special Issue
Data-Driven Characterization of Metabolome Reprogramming during Early Development of Sorghum Seedlings
Previous Article in Journal
Benchmark Dataset for Training Machine Learning Models to Predict the Pathway Involvement of Metabolites
Previous Article in Special Issue
Metabolomic Reconfiguration in Primed Barley (Hordeum vulgare) Plants in Response to Pyrenophora teres f. teres Infection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolic Rewiring in Tea Plants in Response to Gray Blight Disease Unveiled by Multi-Omics Analysis

1
Tea Research Institute, Fujian Academy of Agricultural Sciences, Fuzhou 350013, China
2
Center for Plant Metabolomics, Haixia Institute of Science and Technology, Fujian Agriculture and Forestry University, Fuzhou 350002, China
3
State Key Laboratory of Ecological Pest Control for Fujian and Taiwan Crops, College of Plant Protection, Fujian Agriculture and Forestry University, Fuzhou 350002, China
4
Fuzhou Institute of Oceanography, Minjiang University, Fuzhou 350108, China
*
Authors to whom correspondence should be addressed.
Metabolites 2023, 13(11), 1122; https://doi.org/10.3390/metabo13111122
Submission received: 13 October 2023 / Revised: 25 October 2023 / Accepted: 31 October 2023 / Published: 1 November 2023
(This article belongs to the Special Issue Metabolomics and Plant Defence)

Abstract

:
Gray blight disease, which is caused by Pestalotiopsis-like species, poses significant challenges to global tea production. However, the comprehensive metabolic responses of tea plants during gray blight infection remain understudied. Here, we employed a multi-omics strategy to characterize the temporal transcriptomic and metabolomic changes in tea plants during infection by Pseudopestalotiopsis theae, the causal agent of gray blight. Untargeted metabolomic profiling with ultra-performance liquid chromatography–quadrupole time-of-flight mass spectrometry (UPLC-QTOFMS) revealed extensive metabolic rewiring over the course of infection, particularly within 24 h post-inoculation. A total of 64 differentially accumulated metabolites were identified, including elevated levels of antimicrobial compounds such as caffeine and (−)-epigallocatechin 3-gallate, as well as oxidative catechin polymers like theaflavins, theasinensins and theacitrins. Conversely, the synthesis of (+)-catechin, (−)-epicatechin, oligomeric proanthocyanidins and flavonol glycosides decreased. Integrated omics analyses uncovered up-regulation of phenylpropanoid, flavonoid, lignin biosynthesis and down-regulation of photosynthesis in response to the pathogen stress. This study provides novel insights into the defense strategies of tea plants against gray blight disease, offering potential targets for disease control and crop improvement.

1. Introduction

The tea plant (Camellia sinensis) is an important cash crop cultivated primarily for its leaves, which are used as raw material for tea production. China is the world’s largest tea producer and exporter, with an annual production of tea leaves exceeding 13 million tons in 2021, supporting the livelihoods of millions in the tea industry [1]. However, tea production faces persistent challenges from fungal diseases [2]. Pestalotiopsis-like species, known as the causal agents of gray blight disease in tea plants, are a highly destructive group of phytopathogens [3].
Pestalotiopsis-like species, which are classified into Pseudopestalotiopsis, Neopestalotiopsis and Pestalotiopsis, infect both tender shoots and mature leaves [4]. The initial symptoms of gray blight appear as small brown concentric spots on wounded leaves, which gradually expand into large necrotic lesions with black and brown colors. In severe cases, these lesions spread throughout the leaf, leading to defoliation [5]. Disease severity is exacerbated under warm and humid conditions, negatively impacting tea yield and quality [2]. Various Pestalotiopsis-like species have caused significant yield losses in major tea-producing countries like China, India and Japan [6,7,8]. Current disease management relies heavily on synthetic fungicides such as methyl benzimidazole carbamates and dithiocarbamates. However, these approaches come with high costs, environmental concerns and an increased risk of fungal resistance development [2,9,10]. Therefore, developing effective control measures requires a profound understanding of the defense mechanisms employed by the tea plant against the causal agent of gray blight disease.
Plants have evolved a sophisticated innate immune system to fend off pathogens through multilayered defense responses [11]. Physical structures such as thorns, trichomes, rigid cell walls and waxy cuticles form the first line of defense by preventing pathogen entry. In addition to these physical barriers, plants can recognize invading pathogens and initiate complex signaling cascades that trigger the targeted production of pathogenesis-related proteins, enzymes and specialized metabolites [12]. Among these metabolites are antimicrobial compounds, including constitutively expressed phytoanticipins and induced phytoalexins synthesized de novo in response to infection. These metabolites can directly inhibit microbial growth or elicit additional immune responses, serving to contain pathogens within the infected plant tissues [13]. Tea plants (Camellia sinensis), in particular, reply on the production of specialized metabolites (e.g., phenolics, alkaloids and terpenes) as a key defense strategy against biotic and abiotic stresses [14]. Phenolic compounds, such as (−)-epigallocatechin (EGC), (−)-epicatechin (EC) and their gallate esters, play a crucial role in leaf resistance. Among these phenolics, (−)-epigallocatechin 3-gallate (EGCG), the most abundant catechin in fresh tea leaves, is particularly noteworthy. EGCG has been shown to possess potent antibacterial, antifungal and antiviral properties [15,16]. Through a multifaceted mechanism of action, EGCG effectively combats fungi via disruption of cell structures, inhibition of key enzymes and synergistic potentiation of antifungal drugs, enabling robust activity against fungal invaders [16]. When infected by pathogens such as Colletotrichum fructicola or Pseudopestalotiopsis isolates, tea plants up-regulate EGCG or other catechin levels in correlation with enhanced resistance [17,18]. Nonetheless, the full scope of the specialized metabolic adaptation of tea plants during infection remains poorly understood.
In our previous study, we identified a novel pathogenic strain, Pseudopestalotiopsis theae CYF27, as the cause of gray blight disease in tea plants [19]. To further explore the molecular responses of tea plants to P. theae infection, we conducted an integrated multi-omics analysis. Through untargeted metabolomics, we characterized the temporal dynamics of the tea plant metabolome during P. theae infection. Additionally, through RNA sequencing, we uncovered transcriptional reprogramming events that underlie the immune response. This integrated omics approach provides new perspectives on the intricate host–pathogen interactions over time and may help inform the development of sustainable disease control strategies.

2. Materials and Methods

2.1. Chemical and Reagents

MS-grade acetonitrile and methanol were purchased from Thermo Fisher Scientific, Inc. (Pittsburgh, PA, USA), and formic acid was obtained from Honeywell Fluka (Seelze, Germany). Ultra-pure water was prepared with a Milli-Q purification system (Millipore, Bedford, MA, USA). Reference compounds were obtained from Sigma-Aldrich (St. Louis, MO, USA), ChemFaces (Wuhan, China), Yuanye Biotechnology Inc. (Shanghai, China) and BioBioPha Co., Ltd. (Kunming, China).

2.2. Plant Materials, Pathogen Inoculation and Sampling

One-year-old tea seedlings (C. sinensis cv. “Tieguanyin”) were obtained from a tea plantation in Anxi, Fujian, China (118°13′ E, 25°08′ N). The seedlings were cultivated hydroponically using Hoagland’s nutrient solution (pH 5.6) in a greenhouse at Fujian Agriculture and Forestry University (Fuzhou, China) under the following conditions: 25 ± 3 °C, 16 h photoperiod and 65 ± 5% relative humidity. To ensure the optimal growth of tea seedlings, the nutrient solution was aerated using an oxygen pump and replaced every seven days. After growing for over 30 days, seedlings that displayed uniform growth and were free from any signs of disease or insect infestation were selected as test materials.
The pathogenic strain P. theae CYF27, which has been previously reported [19], was used for inoculation in this study. CYF27 was first incubated at 25 °C on potato dextrose agar (PDA) medium for 20 days. PDA discs (6 mm) containing mycelia and conidia were excised to serve as the inocula. Prior to inoculation, tea leaves were rinsed with sterile water and air-dried to remove any residues from the surface. Each tea seedling received inoculation at three fully expanded leaves with two inoculation sites per leaf. Each site was gently scratched three times (~5 mm in length) with a sterile needle. A 5 mm disc of PDA medium bearing mycelia and conidia was placed onto the upper surface of the wounded tea leaf. The inoculated leaves were covered with cotton dampened in sterile water and wrapped in plastic. The seedlings were then transferred to an inoculation chamber (26 ± 3 °C, 16 h photoperiod, 90 ± 5% relative humidity) for 24 h of dark incubation. At 48 h post-inoculation, the coverings were removed to allow the seedlings to continue growing until harvest. The control samples were wounded tea leaves inoculated with blank PDA discs of the same size (hereafter referred as “CK”). For metabolome and transcriptome analyses, the entire fungal-inoculated leaves (hereafter referred as “PT”), along with the corresponding whole leaves from the CK samples, were harvested at 0, 1, 3 and 6 days post-inoculation (dpi) in three independent biological replicates. Each replicate consisted of three individual seedlings. Samples were snap-frozen in liquid nitrogen and subsequently stored at −80 °C for further analyses.

2.3. Metabolite Extraction and UPLC-QTOFMS Analysis

To prepare samples for metabolite analysis, the freeze-dried samples were ground to fine powders with pre-chilled mortars and pestles. Approximately 30 mg of leaf powder was weighted and extracted with 1.0 mL of 70% aqueous methanol. The resulting mixture was subjected to sonication for 20 min, followed by centrifugation at 12,000× g for 10 min. The supernatant was then filtered using a 0.22 μm PVDF filter (Millipore).
The metabolomics data were acquired on a Waters Acquity UPLC system (Milford, MA, USA) coupled with a Waters SYNAPT G2-Si HDMS QTOF mass spectrometer (Manchester, UK). The instrument was operated in electrospray ionization (ESI) mode and controlled with MassLynx 4.2 software (Waters). UPLC separation was carried out on a Waters Acquity UPLC HSS T3 column (1.8 μm, 2.1 × 100 mm). The mobile phase consisted of water containing 0.1% formic acid (Solvent A) and acetonitrile containing 0.1% formic acid (Solvent B). Gradient elution was performed as per the following program: 0–2 min (1–7% B), 2–13 min (7–40% B), 13–14 min (40–99% B), 14–18 min (held at 99% B) and 5 min of re-equilibration time before the next injection. The CK and PT samples were injected randomly. The ESI parameters were as follows: source temperature, 100 °C; desolvation temperature, 350 °C; cone gas flow, 50 L/h; desolvation gas flow, 800 L/h; capillary voltage, 1.28 kV; sampling cone voltage, 40 V. Data were acquired in both negative and positive ionization modes, operating in a full-scan mode over a mass range of 50–1200 m/z. The MSe data were collected in the continuum mode with a collision energy ramp ranging from 10 to 50 eV. Online calibration of MS data was achieved by continually infusing leucine enkephalin (Waters, Milford, MA, USA) at a rate of 2 ng/min.

2.4. Metabolite Annotation and Multivariate Analysis

Progenesis QI software (v2.4, Nonlinear Dynamics) was used to process the raw data acquired in the positive and negative ionization modes separately. We detailed the process of data processing and metabolite identification in a prior study [20]. Briefly, raw data were imported into Progenesis QI for peak alignment, picking and normalization (with normalization to all compounds) using the default settings. To achieve the fusion of precursor ions, different adduct ion forms were grouped. These included [M−H], [2M−H] and [M+FA−H] in ESI and [M+H]+, [2M+H]+ and [M+Na]+ in ESI+. The generated matrix involving the information of retention time, m/z and normalized peak abundance from each mode was exported and combined for subsequent multivariate statistical analyses in Simca-P 14.1 software (Umetrics, Umea, Sweden). After Pareto scaling, principal component analysis (PCA) and orthogonal partial least-squares discriminant analysis (OPLS-DA) were applied to identify metabolites with significant differences between the CK and PT groups at each time point, based on a fold change (FC) threshold > 2, variable importance in projection (VIP) value > 1 and p-value < 0.05. The differential metabolites were annotated by referencing an in-house MS database as well as public databases such as Metlin, MassBank and HMDB [20,21,22,23].

2.5. Transcriptome Analysis

To explore the temporal changes in transcriptional response during P. theae infection, we performed RNA-seq analysis on pathogen-inoculated leaf samples harvested at 0, 1, 3 and 6 dpi (designated as PT0, PT1, PT3 and PT6, respectively). Briefly, total RNA was extracted using the RNAprep Pure Plant Kit (Tiangen, Beijing, China) according to the manufacturer’s instructions. Sequencing libraries were prepared using the NEBNext Ultra RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The resulting library was then paired-end-sequenced (150 bp reads) on an Illumina HiSeq2500 platform (Novogene Biotech, Beijing, China) following standard protocols. The raw sequence data were deposited in the Genome Sequence Archive (GSA) database under the accession no. CRA011958.
After quality trimming, the clean reads were mapped to the reference transcriptome (Camellia sinensis cv. “Suchazao”) from the TPIA database (http://tpia.teaplants.cn, accessed on 16 June 2022) using the Salmon program (v0.9.0) [24] after the exclusion of reads that contained P. theae sequences. Gene expression was quantified as transcripts per million (TPM). DESeq2 was employed to identify differentially expressed genes (DEGs) with a false discovery rate (FDR) < 0.05 and |log2fold change| > 1. Gene ontology (GO) enrichment analysis were performed using the clusterProfiler package [25]. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed using the KOBAS software [26].

2.6. Quantitative Real-Time PCR (qRT-PCR)

To validate the expression levels of DEGs, qRT-PCR was performed using SYBR Premix Ex TaqTM II (Takara, Dalian, China) on a CFX96TM real-time PCR system (Bio-Rad, Hercules, CA, USA) as previously described [27]. Primers for 12 randomly chosen DEGs were designed using the Primer Premier software (v5.00) (Table S1). The relative expression of genes was determined using the 2−∆∆Ct method [28] with the CsGAPDH gene (KA295375.1) as an endogenous control. Each biological sample was analyzed in triplicate.

2.7. Statistical Analysis

Data were analyzed using GraphPad Prism 8, TBtools (v.1.131) [29] and Microsoft Excel 2010. Results were reported as the mean ± standard deviation (SD) obtained from at least three replicates. The statistical discrepancy between treatments was evaluated using the nonparametric Mann–Whitney U test or one-way ANOVA with a post hoc Fisher’s least significant difference test as appropriate. In all cases, differences were considered significant at p < 0.05.

3. Results and Discussion

3.1. Morphological Analysis of Gray Blight Symptom Development in Tea Leaves

The onset and progression of disease symptoms on tea leaves inoculated with P. theae were monitored over the infection process (Figure 1A). At 0 and 1 dpi, no symptoms were observed on the wounded leaves. At 3 dpi, the symptoms began to manifest, appearing as small brown lesions specifically at the sites where the leaves had been pierced. These lesions gradually enlarged as the infection developed. By 6 dpi, the individual lesions merged together, forming prominent necrotic brown spots across the injured areas. Quantitative analysis showed that the lesion size increased by 4.5-fold from 3 to 6 dpi (Figure 1B). Interestingly, a similar study conducted on the “Suchazao” tea cultivar using a different Pseudopestalotiopsis isolate reported a similar lesion onset at 4 dpi, but the subsequent progression was much slower. In fact, it took 13 dpi for the lesions in “Suchazao” to coalesce into large brown spots resembling those observed at 6 dpi in our study [18]. This disparity implies that the rate at which symptoms spread may vary depending on the specific combination of the tea cultivar and the pathogenic Pseudopestalotiopsis isolate under investigation.

3.2. Metabolic Analysis of Tea Plant Response to P. theae Infection

The changes in catechin contents in tea leaves afflicted by gray blight disease have been explored previously [18]. However, the impact of P. theae infection on other defense-related metabolites (e.g., phytoanticipins and phytoalexins) in tea plants remains largely unknown. To gain comprehensive insights into how tea metabolism is altered during pathogenesis, we conducted non-targeted metabolomics profiling of tea leaves subjected to P. theae inoculation or mock inoculation at 0, 1, 3 and 6 dpi. The representative total ion chromatograms of tea leaf metabolic fingerprints at 1 dpi in the negative mode are presented in Figure S1. After excluding fungal metabolites, we detected a total of 5362 metabolic features in the negative ionization mode and 3048 in the positive mode across all samples. Signals occurring at low abundance (maximum abundance < 1000) were removed, retaining 677 and 295 features in the negative and positive modes, respectively. Chemometrics analyses were then performed on the merged data from both modes. Our analysis unveiled dynamic metabolic shifts following fungal inoculation, as indicated by the log2-fold differences in the relative metabolite abundance between the infected groups and the control groups (Figure 2A). As expected, the metabolite profiles of tea samples at 0 dpi were similar between the CK and PT groups. However, at 1 dpi, the infected group displayed the most intense metabolic response among all time points examined, with sharp increases in both up-regulated and down-regulated metabolic features. Studies have shown that the early plant–pathogen interactions typically take place within 24 h post-inoculation, during which pathogens can fully invade and colonize the host tissues [30]. Hence, the direct contact between the tea plant and P. theae likely triggered the activation of the plant’s defense system at 1 dpi. When comparing the pathogen-inoculated group to the mock-inoculated group, there was an increased abundance of 222, 212 and 184 metabolic features and a decreased abundance of 231, 172 and 136 features at 1, 3 and 6 dpi, respectively (Figure 2A). From 1 to 6 dpi, the number of up-regulated features in the PT group compared to CK only slightly declined, whereas the number of down-regulated features decreased by 41%. This sustained up-regulation and lessening down-regulation of metabolic features in response to fungal infection suggest complex and dynamic interactions between the tea plant and the pathogen.
PCA analysis showed that, except for samples harvested at 0 dpi, the pathogen-inoculated samples were clearly distinguished from the mock-inoculated samples, with the first two principal components (PCs) explaining 73.8% of the total variance (Figure 2B). This again highlights the significant impact of P. theae infection on metabolite rewiring in tea plants. Within the PT group, one biological replicate harvested at 6 dpi deviated from the other two replicates, which made clustering by inoculation duration less distinct in the PCA plot.
To further elucidate the differences between the fungal-inoculated and mock-inoculated samples, we employed OPLS-DA analysis, a robust supervised multivariate method known for its ability to identify statistically significant variables [31]. The OPLS-DA score plots visually depict the separation between the sample groups based on their metabolite profiles at each time point (Figure S2). This analysis enabled us to observe distinct cluster patterns and confirm the significant metabolic alterations associated with the fungal infection. A total of 64 metabolites displayed statistically significant differences (p < 0.05, VIP > 1 and |FC| > 2), as shown in Table 1. These differentially accumulated metabolites (DAMs) were classified into 10 categories: polymerized catechin derivatives (23), flavonol glycosides (13), flavanols (8), flavone and glycosides (2), flavanone glycosides (1), xanthine alkaloids (1), amino acids (1), hydrolysable tannins (1), phenolic acids and derivatives (1) and unknown compounds (13). Flavonoid compounds made up the vast majority (72%) of the identified DAMs, pointing to their key roles in controlling gray blight infection.
A Venn diagram was constructed to visualize the overlap and differences in annotated DAMs between the treatment and control groups (Figure 2C). The number of DAMs showed a notable decrease as the infection progressed, with 45 identified at 1 dpi, 41 at 3 dpi and only 10 at 6 dpi. This indicates that the tea plant’s metabolic response was most pronounced immediately after pathogen invasion but tapered over time. This decline could signify a transition from acute defense to a more balanced state. Alternatively, the pathogen’s manipulation or suppression of certain metabolic pathways may contribute to the decrease. The metabolic alterations coincided with the phenotypic observation of progressively enlarged necrotic spots as the inoculation duration increased (Figure 1). Among the DAMs shared across all time points (Figure 2C and Table 1), four metabolites—caffeine, theaflavin 3,3′-digallate, EGCG and theacitrin C—showed a persistent induction in response to the P. theae challenge compared to CK. As a naturally occurring alkaloid prevalent in plants, caffeine functions as a broad-spectrum deterrent against herbivores and pathogens through its antimicrobial properties [32]. Notably, previous studies have demonstrated the ability of caffeine to inhibit the growth of the tea pathogens P. theae and C. fructicola even at low concentrations [17,33]. This underscores the probable essential role of caffeine in the innate defense response of tea plants to fungal infection. Similarly, reports from various plant–fungus interaction systems have documented a sustained increase in different flavonoid compounds as an integral component of innate immunity [34]. The consistent accumulation of these specialized metabolites in infected tea leaves throughout the infection process may suggest a strategic modulation of metabolism that imparts long-lasting resistance against pathogen challenges. Further evaluation is warranted to explore the potential of these metabolites as candidate markers of resistance and to elucidate how tea plants coordinate chemical defense.

3.3. Metabolic Reprogramming of Flavonoid Biosynthesis in Tea Plants against P. theae Infection

The heatmap visualization of the annotated DAMs reveals fascinating patterns of flavonoid metabolism in tea plants upon pathogen attack (Figure 3A–C). During early infection (1 and 3 dpi), the levels of EGCG, (−)-epigallocatechin (EGC), (−)-epicatechin gallate (ECG) and (−)-epigallocatechin 3-O-(3-O-methyl)gallate (EGCG3″Me) notably increased, while the EC and C levels decreased. Tea polyphenols, with EGCG, EGC, ECG and EC being the dominant compounds in tea leaves, demonstrate effective antifungal properties by inhibiting the growth and spore germination of various phytopathogens [16]. Additionally, studies conducted in Kenya and India have indicated an inverse correlation between the content of tea polyphenols and susceptibility to the pathogen P. theae [35]. These studies suggest the potential of tea polyphenols as antifungal agents for disease control, although the contribution of each polyphenol component to the antifungal activity is yet to be established. Of note, our results differ from those reported by Wang et al., who observed increased EC and GC contents in gray blight-inflicted tea leaves at 4 dpi, with minimal changes in EGCG and ECG [18]. The reasons for the disparities remain unclear but could relate to differences in P. theae isolates and tea cultivars used in the two studies. Therefore, further studies are needed to comprehend the influence of these factors on catechin biosynthesis during fungal infection and to clarify the antifungal effects of individual catechins against P. theae.
In addition to alterations in catechin metabolism, the concentrations of various proanthocyanidin oligomers, including prodelphinidin B, procyanidin dimers, EC-GC dimer, (E)GC-(E)CG dimer and (E)GC-(E)C-(E)C trimer, decreased in tea leaves. This decline could be attributed to the depletion of their precursors, EC and C. There was also a marked reduction in the biosynthesis of flavonol glycosides, such as several derivatives of quercetin, kaempferol and myricetin, which are typically abundant in tea leaves. On the contrary, the concentrations of several oxidative polymerized products of catechins, namely theaflavins (theaflavin, theaflavin 3-gallate, theaflavin 3′-gallate and theaflavin 3, 3′-digallate), theasinensins (theasinensin A, D and F) and theacitrins (theacitrin A and C), increased alongside the elevated levels of their precursors, EGCG and EGC (Figure 3). There are two potential explanations for this accumulation. Firstly, these oxidative polymers might possess additional antifungal or antioxidant properties that directly contribute to the stress response against fungal infection. Alternatively, fungal infection could trigger oxidative stress in tea leaves through the enhanced activities of enzymes like peroxidase and polyphenol oxidase. Such oxidative stress may result in increased catechin polymerization. Although theaflavins, theasinensins and theacitrins are typically found in higher quantities during the fermentation of oolong and black teas, it is worth noting that a study detected increased levels of theaflavins in fresh tea leaves infested by tea green leafhoppers [36,37,38,39]. Nonetheless, the induction of oxidative catechin polymers as a response to biotic stresses has received little exploration. Understanding whether these polymers function directly as defensive compounds or are formed indirectly due to oxidative stress would enhance our comprehension of the metabolic defenses employed by tea plants.
To analyze the variation trends of catechin monomers and polymers during the course of infection, we plotted the relative contents of these compounds in fungal-infected tea samples (Figure 4). EGCG consistently showed a significant increase at 1 dpi (p < 0.05) and maintained higher levels until 6 dpi. The level of EC was significantly lower at 1 dpi compared to 0 dpi (p < 0.05). C levels also tended downward, albeit the change was not statistically significant. As for oxidized catechin polymers, the aforementioned theaflavins, theasinensins and theacitrins all showed a marked surge at 1 dpi and primarily remained elevated at subsequent time points despite some fluctuations.
Collectively, the results show that tea plants strategically reconfigure their flavonoid profiles in response to infection by prioritizing the production of potent antimicrobial flavanols (e.g., EGCG) while reallocating resources away from potentially less bioactive compounds (e.g., EC, C and flavonol glycosides). This orchestrated adjustment likely strengthens the defense against fungal pathogens. Our work offers unique perspectives into the intricate flavonoid-mediated responses of tea plants to control infection.

3.4. Transcriptomic Analysis of Tea Plant Response to P. theae Infection

To further understand the global transcriptome responses modulated by P. theae infection, RNA sequencing was performed on P. theae-inoculated tea leaves at 0, 1, 3 and 6 dpi. Detailed information on RNA sequencing is provided in Table S2. Overall, 330,641,626 clean reads were generated from 12 samples, and these reads were mapped to the reference genome of “Suchazao”, with alignment rates ranging between 86.00% and 89.61%. PCA analysis of the transcriptome data demonstrated clear clustering of biological replicates within each time point, providing evidence of P. theae-induced changes in gene expression (Figure 5A). The first two PCs collectively accounted for 64.6% of the total variance. Samples at 1 dpi were positioned at the far left of the score plot, whereas samples at 0 dpi were located at the far right along the first PC. This is concordant with the metabolomics data showing that the most drastic variation was induced at 1 dpi. A total of 29,037 transcripts were detected in all samples (Table S3). After filtering by FDR < 0.05 and absolute log2 (fold change) > 1, we identified 5967 DEGs (4547 up and 1420 down), 1485 DEGs (1269 up and 216 down) and 2408 DEGs (1708 up and 700 down) in PT_0 dpi vs. PT_1 dpi, PT_0 dpi vs. PT_3 dpi and PT_0 dpi vs. PT_6 dpi comparisons, respectively, suggesting that P. theae infection caused substantial transcriptional reprogramming in the tea plants (Figure 5B). In all comparisons, a higher number of genes were up-regulated compared to those that were down-regulated, which is consistent with the changes at the metabolomic level. In addition, many DEGs were unique to each time point, with 2845, 419 and 322 genes uniquely up-regulated and 1001, 88 and 278 genes uniquely down-regulated at 1, 3 and 6 dpi, respectively (Figure 5C). The reliability of the RNA-seq data was validated through qRT-PCR analysis, which confirmed the expression patterns of 12 randomly selected DEGs in the phenylpropanoid, flavonoid, shikimate and lignin pathways (Figure S3).
GO classification of DEGs was performed to categorize the unigenes into biological processes, molecular functions and cellular components (Figure S4). Within the biological process category, the DEGs were found to be abundant in GO terms such as “cellular process”, “single organism process”, “metabolic process” and “response to stimulus”. The molecular function category mainly comprised DEGs associated with “binding” and “catalytic activity”. In the cellular component category, highly represented groups included “cell”, “cell part”, “organelle” and “membrane”. We conducted KEGG pathway enrichment analysis on DEGs to investigate the primary metabolic pathways involved in tea plant response. The top 10 enriched up-regulated and down-regulated pathways at 1, 3 and 6 dpi are listed in Figure 6. Key enriched pathways included phenylpropanoid biosynthesis, flavonoid biosynthesis, photosynthesis, amino acid metabolism, starch and sucrose metabolism, plant–pathogen interaction and plant hormone signal transduction. Three pathways—“phenylpropanoid biosynthesis”, “flavonoid biosynthesis” and “photosynthesis-antenna proteins”—were commonly enriched across all time points. Previous studies have demonstrated that pathogen-infected plant leaves often experience a decrease in photosynthetic rates due to disturbances in the photosynthetic machinery [40]. In line with this, our studies also found that P. theae infection, particularly at 1 dpi, disrupted pathways related to photosynthesis, such as “photosynthesis-antenna proteins”, “carbon fixation in photosynthetic organisms” and “porphyrin and chlorophyll metabolism” (Table S4). The disruption of photosynthesis may have broader implications for plant defense, as the chloroplast plays a pivotal role in fortifying plant immunity through the synthesis of a diverse array of molecules, hormones and proteins that possess antimicrobial properties [41].

3.5. Transcriptomic Rewiring of Phenylpropanoid and Flavonoid Metabolism in Tea Plants in Response to P. theae Infection

KEGG pathways for phenylpropanoid and flavonoid biosynthesis were significantly enriched following infection (Figure 6). This correlates well with the elevated flavonoid levels from the metabolomics data, indicating the importance of these pathways in tea plant response to gray blight disease. Similar enrichment of phenylpropanoid and flavonoid biosynthesis was also observed in a prior study involving the interaction between the tea cultivar “Suchazao” and a pathogenic Pseudopestalotiopsis isolate [18]. Therefore, we closely examined the transcriptional profiles of genes involved in phenylpropanoid and flavonoid metabolism in tea leaves.
The role of phenylpropanoids as inducible antimicrobial compounds and signaling molecules in plant–pathogen interactions is well established [42]. Similarly, we observed pronounced increases in the transcription levels of genes encoding phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H) and 4-coumarate-CoA ligase (4CL) within the general phenylpropanoid pathway (Figure 7A). PAL, which catalyzes the initial committed step critical for regulating phenylpropanoid metabolism, showed up-regulation of all six transcripts at 1 dpi, with sustained high expression at late infection stages. C4H and 4CL function sequentially downstream to produce p-coumaroyl CoA, an essential precursor for the biosynthesis of various phenolic compounds, such as flavonoids and lignin monomers [43]. They showed similar expression patterns to PAL (Figure 7A,B).
With only some exceptions, genes in the flavonoid pathway were predominantly up-regulated in response to the pathogen at 1 dpi (Figure 7A). Metabolite profiling demonstrated differential accumulation of flavanols and flavonol glycosides. However, differential expression of genes responsible for directing metabolic flux towards specific branches of flavanols or flavonols was not evident. This disconnect between metabolite and transcriptional changes hints at the possibility of additional post-transcriptional controls in regulating flavonoid biosynthesis in tea plants challenged by P. theae.
Also being generated via the phenylpropanoid pathway, lignin comprises a highly branched polymer of monolignols that functions as a fundamental building block of plant cell walls. During the early defense response to pathogen infection (1 and 3 dpi), most lignin biosynthetic genes in tea leaves were up-regulated, though there were some instances of down-regulation for certain transcripts (Figure 7A). This is consistent with the well-documented role of lignin in enhancing plant tolerance to various biotic and abiotic stresses across plant species [44]. Aside from serving as a barricade for pests and pathogens, mounting evidence shows that cell wall lignin dynamically adjusts its concentration and composition in response to environmental stimuli, providing protection in a wide range of contexts [44,45,46]. The resilience of lignin against microbial degradation suggests that increased lignification may fortify the infected cells and restrict pathogen spread [45]. It should be noted that variations in lignin accumulation under fungal infection were not investigated in our study, and future qualitative and quantitative analyses of lignin dynamics will provide additional insights. Intriguingly, peroxidase and laccase transcripts universally showed up-regulation at 1 dpi (Figure 7A). These enzymes are actively involved in the polymerization of monolignols into lignin polymers [47]. Isolation and biochemical analysis of individual isoenzymes in the future experiments will disentangle their roles in lignin formation and pathogen protection.

4. Conclusions

Through integrated “omics” analyses, this study uncovered the dynamic global transcriptomic and metabolomic changes in tea plants challenged with P. theae. By employing untargeted metabolomics coupled with multivariate testing, we identified 64 differentially accumulated metabolites between infected and control groups across multiple time points, highlighting the significant impact of this pathogen on tea plant metabolism. Rapid and substantial metabolic rewiring was observed at 1 dpi, being characterized by enhanced synthesis of caffeine, EGCG and oxidative catechin polymers, as well as a reduction in EC, proanthocyanidins and flavonol glycosides. Combined analysis of metabolomic and transcriptomic data implicates the involvement of phenylpropanoid, flavonoid and lignin pathways in tea plant resistance. Notably, we reveal for the first time the differential modulation of different branches of the flavonoid biosynthetic pathway in the tea plant in response to fungal challenge. These findings shed light on the complex and dynamic interactions between the tea plant and the causal agent of gray blight. Furthermore, our study enables the identification of key defensive genes, metabolites and pathways with potential for disease resistance breeding. Given that resistance can vary among tea cultivars, future comparative studies exploring varieties with varying degrees of susceptibility will provide deeper insights into the resistance mechanisms and aid in the development of more effective strategies for disease management and crop improvement.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/metabo13111122/s1, Figure S1: representative total ion chromatograms of mock-inoculated and fungal-inoculated tea leaves at 1 dpi in the negative ionization mode; Figure S2: OPLS-DA analysis of metabolomics datasets; Figure S3: qRT-PCR validation of RNA-seq DEGs; Figure S4: gene ontology (GO) enrichment analysis of differentially expressed genes in gray blight-infected tea plants at 1, 3 and 6 dpi; Table S1: list of primers for qRT-PCR; Table S2: the quality of RNA-seq data; Table S3: gene detected in all samples; Table S4: KEGG enrichment analysis.

Author Contributions

Conceptualization, X.Y. and Z.W.; methodology, Z.D., X.W. and C.Z.; formal analysis, S.Z.; investigation, S.Z. and X.Y.; data curation, S.Z.; writing—original draft preparation, S.Z. and X.Y.; writing—review and editing, X.Y. and Z.W.; supervision, X.Y. and Z.W.; funding acquisition, X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fujian Agriculture and Forestry University (FAFU) Construction Project for Technological Innovation and Service System of Tea Industry Chain (K1520005A02) and Fundamental Research Project of Fujian Provincial Research Institute for Public Welfare, China (Grant No. 2023R1027002).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We would like to thank Hongli Hu of Fujian Agriculture and Forestry University for providing the strain used in our study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Food and Agriculture Organization Statistics. FAOSTAT (Online). 2021. Available online: http://www.fao.org/faostat/en/#data/QC (accessed on 16 June 2022).
  2. Pandey, A.K.; Sinniah, G.D.; Babu, A.; Tanti, A. How the global tea industry copes with fungal diseases-challenges and opportunities. Plant Dis. 2021, 105, 1868–1879. [Google Scholar] [CrossRef]
  3. Wang, S.; Mi, X.; Wu, Z.; Zhang, L.; Wei, C. Characterization and pathogenicity of Pestalotiopsis-like species associated wth gray blight disease on Camellia sinensis in Anhui Province, China. Plant Dis. 2019, 103, 2786–2797. [Google Scholar] [CrossRef] [PubMed]
  4. Maharachchikumbura, S.S.N.; Guo, L.-D.; Cai, L.; Chukeatirote, E.; Wu, W.P.; Sun, X.; Crous, P.W.; Bhat, D.J.; McKenzie, E.H.C.; Bahkali, A.H.; et al. A multi-locus backbone tree for Pestalotiopsis, with a polyphasic characterization of 14 new species. Fungal Divers. 2012, 56, 95–129. [Google Scholar] [CrossRef]
  5. Chen, Y.; Zeng, L.; Shu, N.; Jiang, M.; Wang, H.; Huang, Y.; Tong, H. Pestalotiopsis-like species causing gray blight disease on Camellia sinensis in China. Plant Dis. 2018, 102, 98–106. [Google Scholar] [CrossRef] [PubMed]
  6. Takeda, Y. Genetic analysis of tea gray blight resistance in tea plants. Jpn. Agric. Res. Q. 2002, 36, 143–150. [Google Scholar] [CrossRef]
  7. Joshi, S.; Sanjay, R.; Baby, U.; Mandal, A.K. Molecular characterization of Pestalotiopsis spp. associated with tea (Camellia sinensis) in southern India using RAPD and ISSR markers. Indian J. Biotechnol. 2009, 8, 377–383. [Google Scholar]
  8. Liu, F.; Hou, L.; Raza, M.; Cai, L. Pestalotiopsis and allied genera from Camellia, with description of 11 new species from China. Sci. Rep. 2017, 7, 866. [Google Scholar] [CrossRef]
  9. Pandey, A.K.; Hubbali, M.; Vandana; Dutta, P.; Babu, A. Characterization and identification of fungicide insensitive Pestalotiopsis-like species pathogenic to tea crop in India. World J. Microbiol. Biotechnol. 2022, 39, 34. [Google Scholar] [CrossRef]
  10. Yamada, K.; Sonoda, R.; Ishikawa, K. Population genetic structure of QoI-resistant Pestalotiopsis longiseta isolates causing tea gray blight. Plant Dis. 2016, 100, 1686–1691. [Google Scholar] [CrossRef] [PubMed]
  11. Nishad, R.; Ahmed, T.; Rahman, V.J.; Kareem, A. Modulation of plant defense system in response to microbial interactions. Front. Microbiol. 2020, 11, 1298. [Google Scholar] [CrossRef]
  12. Lindsay, W.P.; Lamb, C.J.; Dixon, R.A. Microbial recognition and activation of plant defense systems. Trends Microbiol. 1993, 1, 181–187. [Google Scholar] [CrossRef]
  13. González-Lamothe, R.; Mitchell, G.; Gattuso, M.; Diarra, M.S.; Malouin, F.; Bouarab, K. Plant antimicrobial agents and their effects on plant and human pathogens. Int. J. Mol. Sci. 2009, 10, 3400–3419. [Google Scholar] [CrossRef]
  14. Zeng, L.; Zhou, X.; Liao, Y.; Yang, Z. Roles of specialized metabolites in biological function and environmental adaptability of tea plant (Camellia sinensis) as a metabolite studying model. J. Adv. Res. 2021, 34, 159–171. [Google Scholar] [CrossRef] [PubMed]
  15. Friedman, M. Overview of antibacterial, antitoxin, antiviral, and antifungal activities of tea flavonoids and teas. Mol. Nutr. Food Res. 2007, 51, 116–134. [Google Scholar] [CrossRef]
  16. Yang, Y.; Zhang, T. Antimicrobial activities of tea polyphenol on phytopathogens: A review. Molecules 2019, 24, 816. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, L.; Wang, Y.; Cao, H.; Hao, X.; Zeng, J.; Yang, Y.; Wang, X. Transcriptome Analysis of an Anthracnose-Resistant Tea Plant Cultivar Reveals Genes Associated with Resistance to Colletotrichum camelliae. PLoS ONE 2016, 11, e0148535. [Google Scholar] [CrossRef]
  18. Wang, S.; Liu, L.; Mi, X.; Zhao, S.; An, Y.; Xia, X.; Guo, R.; Wei, C. Multi-omics analysis to visualize the dynamic roles of defense genes in the response of tea plants to gray blight. Plant J. 2021, 106, 862–875. [Google Scholar] [CrossRef]
  19. Zheng, S.; Chen, R.; Wang, Z.; Liu, J.; Cai, Y.; Peng, M.; Zhang, T.; Li, Y.; Wang, B.; Bao, J.; et al. High-quality genome assembly of Pseudopestalotiopsis theae, the pathogenic fungus causing tea gray blight. Plant Dis. 2021, 105, 3723–3726. [Google Scholar] [CrossRef]
  20. Yu, X.; Xiao, J.; Chen, S.; Yu, Y.; Ma, J.; Lin, Y.; Li, R.; Lin, J.; Fu, Z.; Zhou, Q.; et al. Metabolite signatures of diverse Camellia sinensis tea populations. Nat. Commun. 2020, 11, 5586. [Google Scholar] [CrossRef] [PubMed]
  21. Montenegro-Burke, J.R.; Guijas, C.; Siuzdak, G. METLIN: A tandem mass spectral library of standards. Methods Mol. Biol. 2020, 2104, 149–163. [Google Scholar] [CrossRef]
  22. Horai, H.; Arita, M.; Kanaya, S.; Nihei, Y.; Ikeda, T.; Suwa, K.; Ojima, Y.; Tanaka, K.; Tanaka, S.; Aoshima, K.; et al. MassBank: A public repository for sharing mass spectral data for life sciences. J. Mass Spectrom. 2010, 45, 703–714. [Google Scholar] [CrossRef] [PubMed]
  23. Wishart, D.S.; Guo, A.; Oler, E.; Wang, F.; Anjum, A.; Peters, H.; Dizon, R.; Sayeeda, Z.; Tian, S.; Lee, B.L.; et al. HMDB 5.0: The Human Metabolome Database for 2022. Nucleic Acids Res. 2022, 50, D622–D631. [Google Scholar] [CrossRef] [PubMed]
  24. Patro, R.; Duggal, G.; Love, M.I.; Irizarry, R.A.; Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 2017, 14, 417–419. [Google Scholar] [CrossRef] [PubMed]
  25. Yu, G.; Wang, L.G.; Han, Y.; He, Q.Y. clusterProfiler: An R package for comparing biological themes among gene clusters. Omics A J. Integr. Biol. 2012, 16, 284–287. [Google Scholar] [CrossRef] [PubMed]
  26. Mao, X.; Cai, T.; Olyarchuk, J.G.; Wei, L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 2005, 21, 3787–3793. [Google Scholar] [CrossRef]
  27. Zhao, X.; Chen, S.; Wang, S.; Shan, W.; Wang, X.; Lin, Y.; Su, F.; Yang, Z.; Yu, X. Defensive responses of tea plants (Camellia sinensis) against tea green leafhopper attack: A multi-omics study. Front. Plant Sci. 2020, 10, 1705. [Google Scholar] [CrossRef]
  28. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  29. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  30. Shen, Y.; Liu, N.; Li, C.; Wang, X.; Xu, X.; Chen, W.; Xing, G.; Zheng, W. The early response during the interaction of fungal phytopathogen and host plant. Open Biol. 2017, 7, 170057. [Google Scholar] [CrossRef]
  31. Wiklund, S.; Johansson, E.; Sjöström, L.; Mellerowicz, E.; Edlund, U.; Shockcor, J.; Gottfries, J.; Moritz, T.; Trygg, J. Visualization of GC/TOF-MS-based metabolomics data for identification of biochemically interesting compounds using OPLS class models. Anal. Chem. 2008, 80, 115–122. [Google Scholar] [CrossRef]
  32. Kim, Y.-S.; Choi, Y.E.; Sano, H. Plant vaccination: Stimulation of defense system by caffeine production in planta. Plant Signal. Behav. 2010, 5, 489–493. [Google Scholar] [CrossRef] [PubMed]
  33. Zhang, Z. Antifungal activity of caffeine against fungal pathogens of tea plant. J. Nanjing Agric. Univ. 2010, 2, 63–67. [Google Scholar]
  34. Ramaroson, M.-L.; Koutouan, C.; Helesbeux, J.-J.; Clerc, V.; Latifa, H.; Geoffriau, E.; Briard, M. Role of phenylpropanoids and flavonoids in plant resistance to pests and diseases. Molecules 2022, 27, 8371. [Google Scholar] [CrossRef]
  35. Chen, Z.M.; Sun, X.L.; Dong, W.X. Genetics and chemistry of the resistance of tea plant to pests. In Global Tea Breeding; Chen, L., Apostolides, Z., Chen, Z.-M., Eds.; Zhejiang University Press: Zhejiang, China, 2012; pp. 343–360. [Google Scholar]
  36. Balentine, D.A.; Wiseman, S.A.; Bouwens, L.C. The chemistry of tea flavonoids. Crit. Rev. Food Sci. Nutr. 1997, 37, 693–704. [Google Scholar] [CrossRef]
  37. Davis, A.L.; Lewis, J.R.; Cai, Y.; Powell, C.; Davis, A.P.; Wilkins, J.P.G.; Pudney, P.; Clifford, M.N. A polyphenolic pigment from black tea. Phytochemistry 1997, 46, 1397–1402. [Google Scholar] [CrossRef]
  38. Matsuo, Y.; Li, Y.; Watarumi, S.; Tanaka, T.; Kouno, I. Production and degradation mechanism of theacitrin C, a black tea pigment derived from epigallocatechin-3-O-gallate via a bicyclo[3.2.1]octane-type intermediate. Tetrahedron 2011, 67, 2051–2059. [Google Scholar] [CrossRef]
  39. Liao, Y.; Yu, Z.; Liu, X.; Zeng, L.; Cheng, S.; Li, J.; Tang, J.; Yang, Z. Effect of major tea insect attack on formation of quality-related nonvolatile specialized metabolites in tea (Camellia sinensis) leaves. J. Agric. Food. Chem. 2019, 67, 6716–6724. [Google Scholar] [CrossRef]
  40. Duan, S.; Jin, J.; Gao, Y.; Jin, C.; Mu, J.; Zhen, W.; Sun, Q.; Xie, C.; Ma, J. Integrated transcriptome and metabolite profiling highlights the role of benzoxazinoids in wheat resistance against Fusarium crown rot. Crop J. 2022, 10, 407–417. [Google Scholar] [CrossRef]
  41. Lu, Y.; Yao, J. Chloroplasts at the crossroad of photosynthesis, pathogen infection and plant defense. Int. J. Mol. Sci. 2018, 19, 3900. [Google Scholar] [CrossRef]
  42. Naoumkina, M.A.; Zhao, Q.; Gallego-Giraldo, L.; Dai, X.; Zhao, P.X.; Dixon, R.A. Genome-wide analysis of phenylpropanoid defence pathways. Mol. Plant Pathol. 2010, 11, 829–846. [Google Scholar] [CrossRef]
  43. Volpi e Silva, N.; Mazzafera, P.; Cesarino, I. Should I stay or should I go: Are chlorogenic acids mobilized towards lignin biosynthesis? Phytochemistry 2019, 166, 112063. [Google Scholar] [CrossRef] [PubMed]
  44. Yadav, S.; Chattopadhyay, D. Lignin: The building block of defense responses to stress in plants. J. Plant Growth Regul. 2023, 42, 6652–6666. [Google Scholar] [CrossRef]
  45. Moura, J.C.; Bonine, C.A.; de Oliveira Fernandes Viana, J.; Dornelas, M.C.; Mazzafera, P. Abiotic and biotic stresses and changes in the lignin content and composition in plants. J. Integr. Plant Biol. 2010, 52, 360–376. [Google Scholar] [CrossRef] [PubMed]
  46. Ithal, N.; Recknor, J.; Nettleton, D.; Maier, T.; Baum, T.J.; Mitchum, M.G. Developmental transcript profiling of cyst nematode feeding cells in soybean roots. Mol. Plant-Microbe Interact. 2007, 20, 510–525. [Google Scholar] [CrossRef]
  47. Boerjan, W.; Ralph, J.; Baucher, M. Lignin Biosynthesis. Annu. Rev. Plant Biol. 2003, 54, 519–546. [Google Scholar] [CrossRef]
Figure 1. Development of gray blight disease on tea leaves over time. Intact leaves were inoculated with mycelia and conidia mixture of P. theae. (A) Representative symptoms on tea leaves using intact plant inoculation. Scale bar = 1 cm. (B) Statistics of disease development on tea leaves assessed at 0, 1, 3 and 6 dpi following fungal inoculation. The combined leaf disease area (cm2) from two inoculations per leaf was calculated as a percentage of the total leaf area using ImageJ software.
Figure 1. Development of gray blight disease on tea leaves over time. Intact leaves were inoculated with mycelia and conidia mixture of P. theae. (A) Representative symptoms on tea leaves using intact plant inoculation. Scale bar = 1 cm. (B) Statistics of disease development on tea leaves assessed at 0, 1, 3 and 6 dpi following fungal inoculation. The combined leaf disease area (cm2) from two inoculations per leaf was calculated as a percentage of the total leaf area using ImageJ software.
Metabolites 13 01122 g001
Figure 2. Time course analysis of tea metabolome reconfiguration in response to gray blight infection. (A) Overview of altered metabolic features. The violin plot shows the fold change values for the fungal-inoculated group relative to the mock-inoculated group at each time point. (B) Principal component analysis of tea leaves inoculated with P. theae or the mock control. (C) Venn diagram of differential metabolites between fungal-inoculated and mock-inoculated groups at each time point. PT, tea leaves inoculated with P. theae. CK, tea leaves inoculated with the mock control.
Figure 2. Time course analysis of tea metabolome reconfiguration in response to gray blight infection. (A) Overview of altered metabolic features. The violin plot shows the fold change values for the fungal-inoculated group relative to the mock-inoculated group at each time point. (B) Principal component analysis of tea leaves inoculated with P. theae or the mock control. (C) Venn diagram of differential metabolites between fungal-inoculated and mock-inoculated groups at each time point. PT, tea leaves inoculated with P. theae. CK, tea leaves inoculated with the mock control.
Metabolites 13 01122 g002
Figure 3. Heatmaps of annotated metabolites showing significant differences between fungal-inoculated and mock-inoculated tea leaves at 1 dpi (A), 3 dpi (B) and 6 dpi (C). Red and blue colors signify increased and decreased levels of metabolites. PT, tea leaves inoculated with P. theae. CK, tea leaves inoculated with the mock control.
Figure 3. Heatmaps of annotated metabolites showing significant differences between fungal-inoculated and mock-inoculated tea leaves at 1 dpi (A), 3 dpi (B) and 6 dpi (C). Red and blue colors signify increased and decreased levels of metabolites. PT, tea leaves inoculated with P. theae. CK, tea leaves inoculated with the mock control.
Metabolites 13 01122 g003
Figure 4. Quantitative analysis of catechins and polymerized catechin derivatives extracted from metabolomics data. The y-axis denotes peak area values obtained through normalization to the “normalize to all compounds function” in Progenesis QI. The relative abundance of each metabolite is presented as mean ± SD (n = 3). Different letters denote significant differences in metabolite contents at different time points, as determined by Fisher’s least significant difference test at p < 0.05.
Figure 4. Quantitative analysis of catechins and polymerized catechin derivatives extracted from metabolomics data. The y-axis denotes peak area values obtained through normalization to the “normalize to all compounds function” in Progenesis QI. The relative abundance of each metabolite is presented as mean ± SD (n = 3). Different letters denote significant differences in metabolite contents at different time points, as determined by Fisher’s least significant difference test at p < 0.05.
Metabolites 13 01122 g004
Figure 5. Time course analysis of tea transcriptome reprogramming in response to gray blight infection. (A) Principal component analysis of fungal-inoculated tea leaves at 0, 1, 3 and 6 dpi. (B) Number of differentially expressed genes (DEGs) at each time point. Teal and orange columns indicate up- and down-regulated genes, respectively. (C) Venn diagrams illustrating the number of DEGs up- or down-regulated due to fungal treatment over the time course.
Figure 5. Time course analysis of tea transcriptome reprogramming in response to gray blight infection. (A) Principal component analysis of fungal-inoculated tea leaves at 0, 1, 3 and 6 dpi. (B) Number of differentially expressed genes (DEGs) at each time point. Teal and orange columns indicate up- and down-regulated genes, respectively. (C) Venn diagrams illustrating the number of DEGs up- or down-regulated due to fungal treatment over the time course.
Metabolites 13 01122 g005
Figure 6. KEGG pathway enrichment analysis of DEGs following fungal treatment. The up-regulated pathways are in yellow, while the down-regulated pathways are in blue.
Figure 6. KEGG pathway enrichment analysis of DEGs following fungal treatment. The up-regulated pathways are in yellow, while the down-regulated pathways are in blue.
Metabolites 13 01122 g006
Figure 7. Expression profiles of genes involved in phenylpropanoid, lignin and flavonoid pathways in tea leaves after fungal treatment. (A) Heatmap comparison of the TPM values of related genes in three pathways based on RNA-seq data. (B) Simplified diagram of phenylpropanoid, lignin and flavonoid pathways. PAL, phenylalanine ammonia lyase; C4H, cinnamate 4-hydroxylase; 4CL, 4-coumarate-CoA ligase; CHS, chalcone synthase; CHI, chalcone isomerase; F3H, flavanone 3-hydroxylase; F3′H, flavanone 3′-hydroxylase; F3′5′H, flavanone 3′,5′-hydroxylase; FLS, flavonol synthase; DFR, dihydroflavonol 4-reductase; LAR, leucoanthocyanidin reductase; ANS, anthocyanidin synthase; ANR, anthocyanidin reductase; SCPL1A, type 1A serine carboxypeptidase-like acyltransferases; CCR, cinnamoyl-CoA reductase; HCT, hydroxycinnamoyl-CoA:shikimate/quinate hydroxycinnamoyl transferase; C3′H, p-coumaroyl quinate/shikimate 3′-hydroxylase; CCoAOMT, caffeoyl-CoA O-methyltransferase; CAD, cinnamyl alcohol dehydrogenase; COMT, caffeic acid O-methyltransferase; F5H, ferulate 5-hydroxylase; PER, peroxidase; LAC, laccase.
Figure 7. Expression profiles of genes involved in phenylpropanoid, lignin and flavonoid pathways in tea leaves after fungal treatment. (A) Heatmap comparison of the TPM values of related genes in three pathways based on RNA-seq data. (B) Simplified diagram of phenylpropanoid, lignin and flavonoid pathways. PAL, phenylalanine ammonia lyase; C4H, cinnamate 4-hydroxylase; 4CL, 4-coumarate-CoA ligase; CHS, chalcone synthase; CHI, chalcone isomerase; F3H, flavanone 3-hydroxylase; F3′H, flavanone 3′-hydroxylase; F3′5′H, flavanone 3′,5′-hydroxylase; FLS, flavonol synthase; DFR, dihydroflavonol 4-reductase; LAR, leucoanthocyanidin reductase; ANS, anthocyanidin synthase; ANR, anthocyanidin reductase; SCPL1A, type 1A serine carboxypeptidase-like acyltransferases; CCR, cinnamoyl-CoA reductase; HCT, hydroxycinnamoyl-CoA:shikimate/quinate hydroxycinnamoyl transferase; C3′H, p-coumaroyl quinate/shikimate 3′-hydroxylase; CCoAOMT, caffeoyl-CoA O-methyltransferase; CAD, cinnamyl alcohol dehydrogenase; COMT, caffeic acid O-methyltransferase; F5H, ferulate 5-hydroxylase; PER, peroxidase; LAC, laccase.
Metabolites 13 01122 g007
Table 1. Identification or tentative identification of significantly differential metabolites in response to P. theae infection.
Table 1. Identification or tentative identification of significantly differential metabolites in response to P. theae infection.
No.Metabolite AssignmentRT (min)Metabolite ClassAdductsFormulaTheoretical m/zMeasured m/zPPM ErrorMS/MS FragmentsTime Point
1theanine *1.40amino acids[M−H]C7H14N2O3173.0926173.09312.89155.0826, 128.09453 dpi
2theogallin *2.92phenolic acids and derivatives[M−H]C14H16O10343.0665343.06722.04191.0568, 169.0143, 125.02443 dpi
3prodelphinidin B4.14polymerized catechin derivatives[M−H]C30H26O14609.1244609.12500.99441.0826, 423.0716, 305.0665, 125.02421, 3 dpi
4EGC-(2->7,4->8)-EGCG4.16polymerized catechin derivatives[M−H]C37H28O18759.1197759.1191−0.79605.09291, 3 dpi
5catechin-4-ol 3-O-hexoside4.53flavanols[M−H]C21H24O12467.1190467.1188−0.43305.0659, 125.02463 dpi
6prodelphinidin trimer GC-C-C isomer 14.74polymerized catechin derivatives[M−H]C45H38O20897.1878897.18790.11771.15693 dpi
7EC-GC dimer4.81polymerized catechin derivatives[M+H]+C30H26O13595.1452595.1440−2.02279.09223 dpi
8theacitrin A4.86polymerized catechin derivatives[M−H]C37H28O18759.1197759.1194−0.40741.1062, 571.08691 dpi
9prodelphinidin trimer GC-C-C isomer 24.86polymerized catechin derivatives[M−H]C45H38O20897.1878897.1875−0.33771.1583 dpi
10epigallocatechin *4.91flavanols[M−H]C15H14O7305.0661305.06661.64179.0352, 125.02511, 6 dpi
11unknown4.93unknown[M−H]C45H40O20899.2035899.20350.00ND3 dpi
12(E)GC+(E)C+(E)C isomer 15.14polymerized catechin derivatives[M−H]C45H38O19881.1929881.1926−0.34423.0722, 305.0671, 287.0561, 125.02483 dpi
13catechin *5.33flavanols[M−H]C15H14O6289.0712289.07182.08179.0347, 137.0241, 123.04511 dpi
14strictinin5.34hydrolysable tannins[M−2H]2−C27H22O18633.0728633.07300.32316.0328, 300.9989, 275.01941 dpi
15(E)GC+(E)C+(E)C isomer 25.56polymerized catechin derivatives[M−H]C45H38O19881.1929881.1927−0.23305.0649, 287.05503 dpi
16caffeine *5.56alkaloids[M+H]+C8H10N4O2195.0882195.08935.64138.0671, 110.07200, 1, 3, 6 dpi
17procyanidin dimer5.68polymerized catechin derivatives[M−H]C30H26O12577.1346577.13500.69451.1026, 425.0873, 407.0768, 289.07153 dpi
18theasinensin A5.71polymerized catechin derivatives[M−2H]2−C44H34O22913.1463913.1446−1.86743.1231, 591.1140, 169.0142, 125.02431 dpi
19unknown5.72unknown[M−2H]2−C26H43O36930.1456930.1450−0.65ND1, 3 dpi
20(E)GC-(E)CG dimer5.86polymerized catechin derivatives[M−H]C37H30O17745.1405745.1399−0.81593.1121, 423.0715, 169.0145, 125.02491, 3 dpi
21theasinensin D5.90polymerized catechin derivatives[M−H]C44H34O22913.1463913.1435−3.07423.0709, 285.0396, 169.0140, 125.02413, 6 dpi
22procyanidin trimer6.00polymerized catechin derivatives[M−H]C45H38O18865.1980865.19800.00739.1656, 713.1496, 695.1396, 413.0866, 245.04723 dpi
23carthamidin diglucoside6.03flavonol glycosides[M−H]C27H32O16611.1612611.16282.62491.1191, 449.1241, 397.07711, 3 dpi
24isovitexin glucoside6.03flavone glycosides[M+H]+C27H30O15595.1663595.1656−1.18433.1122, 313.07111 dpi
25(E)GC+(E)C+(E)C isomer 36.10polymerized catechin derivatives[M−H]C45H38O19881.1929881.1922−0.79423.0712, 305.0658, 287.0557, 125.02403 dpi
26unknown6.10unknown[M−H]C38H32O19791.1460791.1453−0.88ND1 dpi
27diglucopyranosyl trihydroxyflavanone6.12flavanone glycosides[M−H]C27H32O15595.1663595.16670.67475.1249, 433.1353, 313.09341, 3 dpi
28epicatechin *6.23flavanols[M−H]C15H14O6289.0712289.07223.46245.0820, 203.0710, 123.04511 dpi
29unknown6.30unknown[M−H]C19H30O8385.1862385.1854−2.08153.09161 dpi
30epigallocatechin gallate *6.35flavanols[M−H]C22H18O11457.0771457.07883.72305.0669, 169.0152, 125.02481, 3, 6 dpi
31theacitrin C6.36polymerized catechin derivatives[M−H]C44H32O22911.1307911.13171.10453.0459, 231.02781, 3, 6 dpi
32unknown6.37unknown[M−H]C13H18O4227.0344227.03471.32ND1, 3 dpi
33unknown6.38unknown[M−H]C19H32O8387.2019387.2006−3.36ND1, 3 dpi
34unknown6.38unknown[M−H]C38H36O19795.1773795.17730.00ND1 dpi
35unknown6.50unknown[M−H]C38H32O19791.1460791.1457−0.38ND1 dpi
36(E)CG+(E)C isomer 16.70polymerized catechin derivatives[M−H]C37H30O16729.1456729.14570.14577.1154, 407.0759, 289.0709, 125.02423 dpi
37unknown6.72unknown[M−H]C38H32O19791.1460791.1459−0.13ND1 dpi
38theasinensin F6.75polymerized catechin derivatives[M−2H]2−C44H34O21897.1514897.1505−1.00575.11863 dpi
39(E)CG+(E)C isomer 26.80polymerized catechin derivatives[M−H]C37H30O16729.1456729.1446−1.37577.1193, 407.0758, 289.0710, 125.02491, 3 dpi
40myricetin 3-neohesperidoside6.96flavonol glycosides[M−H]C27H30O17625.1404625.14050.16316.02241, 3 dpi
41myricetin 3’-glucoside7.15flavonol glycosides[M−H]C21H20O13479.0826479.08290.63316.02251 dpi
42unknown7.35unknown[M−H]C39H50O25917.2563917.25690.65ND1 dpi
43quercetin 3-O-glucosyl rutinoside7.39flavonol glycosides[M−H]C33H40O21771.1984771.19880.52609.1451, 463.0621, 301.0350, 300.02781 dpi
44epigallocatechin 3-(3-O-methylgallate) *7.43flavanols[M−H]C23H20O11471.0927471.09321.06287.0556, 269.04531, 6 dpi
45digalloylprocyanidin dimer7.53polymerized catechin derivatives[M−H]C44H34O20881.1565881.15680.34729.1407, 711.1342, 169.01413 dpi
46rutin * 7.73flavonol glycosides[M−H]C27H30O16609.1456609.14570.16533.1294, 300.0277, 271.0251, 255.0299, 243.02981, 3, 6 dpi
47kaempferol 3-O-galactosyl rutinoside7.75flavonol glycosides[M−H]C33H40O20755.2035755.2030−0.66533.1287, 285.04031 dpi
48epicatechin 3-O-gallate *7.85flavanols[M−H]C22H18O10441.0822441.08332.49331.0458, 289.0720, 169.0146, 125.02473 dpi
49tricetin7.92flavones[M+H]+C15H10O7303.0505303.0504−0.33285.04211, 6 dpi
50quercetin 3-O-galactoside7.93flavonol glycosides[M−H]C21H20O12463.0877463.08810.86300.02761 dpi
51kaempferol 3-O-glucosyl rutinoside8.02flavonol glycosides[M−H]C33H40O20755.2035755.20360.13285.04041 dpi
52quercetin 3-O-glucoside *8.04flavonol glycosides[M−H]C21H20O12463.0877463.0872−1.08301.0344, 300.02721, 3 dpi
53capilliposide I isomer8.17flavonol glycosides[M−H]C48H56O271063.29311063.2919−1.13917.2332, 755.1826, 609.1454, 377.0873, 301.0343, 300.02703 dpi
54unknown8.39unknown[M−H]C24H42O11505.2649505.2645−0.79551.27061, 3 dpi
55kaempferol 3-O-rutinoside *8.46flavonol glycosides[M−H]C27H30O15593.1506593.1504−0.341187.3071, 285.03981, 3 dpi
56unknown8.58unknown[M−H]C20H36O9419.2281419.2269−2.86355.10541 dpi
57epiafzelechin 3-gallate8.94flavanols[M−H]C22H18O9425.0873425.08770.94273.0763, 255.06526 dpi
58sarmenoside II isomer 110.31flavonol glycosides[M−H]C42H46O22901.2402901.24040.22755.1822, 609.1461, 301.0351, 300.02741, 3 dpi
59theaflavin *10.59polymerized catechin derivatives[M−H]C29H24O12563.1190563.11930.53425.0883, 269.0456, 137.02461 dpi
60sarmenoside II isomer 210.91flavonol glycosides[M−H]C42H46O22901.2402901.24020.00755.1823, 609.1449, 301.0346, 300.02701, 3 dpi
61theaflavin 3-gallate *10.97polymerized catechin derivatives[M−H]C36H28O16715.1299715.13020.42563.1178, 423.2219, 125.02441, 3 dpi
62theaflavin 3,3’-digallate *11.17polymerized catechin derivatives[M−H]C43H32O20867.1409867.14130.46697.1193, 389.0657, 178.84231, 3, 6 dpi
63theaflavin-3’-gallate *11.21polymerized catechin derivatives[M−H]C36H28O16715.1299715.1287−1.68563.1173, 125.02381, 3 dpi
64unknown13.78unknown[M−H]C33H54O15689.3384689.3376−1.20671.3262, 653.31620 dpi
* Confirmed by comparison with authentic standards. ND, not determined.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zheng, S.; Du, Z.; Wang, X.; Zheng, C.; Wang, Z.; Yu, X. Metabolic Rewiring in Tea Plants in Response to Gray Blight Disease Unveiled by Multi-Omics Analysis. Metabolites 2023, 13, 1122. https://doi.org/10.3390/metabo13111122

AMA Style

Zheng S, Du Z, Wang X, Zheng C, Wang Z, Yu X. Metabolic Rewiring in Tea Plants in Response to Gray Blight Disease Unveiled by Multi-Omics Analysis. Metabolites. 2023; 13(11):1122. https://doi.org/10.3390/metabo13111122

Chicago/Turabian Style

Zheng, Shiqin, Zhenghua Du, Xiaxia Wang, Chao Zheng, Zonghua Wang, and Xiaomin Yu. 2023. "Metabolic Rewiring in Tea Plants in Response to Gray Blight Disease Unveiled by Multi-Omics Analysis" Metabolites 13, no. 11: 1122. https://doi.org/10.3390/metabo13111122

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

Article Metrics

Back to TopTop