Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Plant Sci., 14 January 2026

Sec. Plant Physiology

Volume 16 - 2025 | https://doi.org/10.3389/fpls.2025.1750493

This article is part of the Research TopicInduced Resistance for a Sustainable Future of AgricultureView all 6 articles

Multi-omics analysis of molecular mechanisms driving the grafting- enhanced resistance of tea plants to Colletotrichum camelliae

Yue-Xin Li,&#x;Yue-Xin Li1,2†Kai-Qin Lin&#x;Kai-Qin Lin1†An-Ran WangAn-Ran Wang1Jie WeiJie Wei1Dong-Hai Yan*Dong-Hai Yan1*De-Gang Zhao,*De-Gang Zhao3,4*
  • 1Tea Research Institute, Guizhou Academy of Agricultural Sciences, Guiyang, China
  • 2Guizhou Provincial Government, Guizhou Key Laboratory of Molecular Breeding for Characteristic Horticultural Crops, Guiyang, China
  • 3The Key Laboratory of Plant Resources Conservation and Germplasm Innovation in Mountainous Region (Ministry of Education), College of Life Sciences/College of Tea Sciences, Guizhou University, Guiyang, China
  • 4Plant Conservation & Breeding Technology Center, Biotechnology Institute of Guizhou Province, Guizhou Academy of Agricultural Sciences, Guiyang, China

Background: The fungal pathogen Colletotrichum camelliae causes a devastating disease that severely limits tea plant (Camellia sinensis) yield and quality. Grafting onto resistant rootstocks offers a sustainable control strategy, yet resistant rootstocks confer scion protection remains obscure.

Methods: Susceptible cultivar ‘Qianmei 818’ was grafted onto resistant ‘Qianmei 419’. Profiled systemic defenses using integrated RNA-seq, sRNA-seq, and metabolomics, complemented by phytohormone and defense-enzyme assays and qRT-PCR validation.

Results: Hetero-grafting conferred near-complete resistance, reducing lesion diameters by 98.71% compared to ungrafted controls, with elevated PAL activity and accumulation of IAA, GA3, and MeSA. Multi-omics analyses identified 1205 differentially expressed genes, 157 differentially expressed miRNAs, and 791 differential metabolites. Pathway integration indicated extensive reprogramming of phenylpropanoid biosynthesis, sulfur metabolism, and plant hormone signaling. Notably, specific miRNA-mRNA regulatory modules, such as downregulation of csi-miR395b-3p and a novel miR397 paralleled up-regulation of their targets in sulfur assimilation (CsAPS1) and lignin biosynthesis (CsCCoAOMT, CsCCR2), respectively, linking miRNA control to reinforcement of structural and biochemical defenses.

Conclusions: Resistant rootstocks activate scion-wide defense networks through miRNA-mediated transcriptome-metabolome remodeling, achieving robust resistance while maintaining tea quality. The elucidated modules provide actionable targets and genetic resources for breeding and grafting strategies toward sustainable disease management.

1 Introduction

Camellia sinensis, commonly known as tea plant, is one of the most widely cultivated and economically important crops globally. A prevalent fungal disease caused by Colletotrichum camelliae (C. camelliae) affects tea plants across major tea-producing regions in China, leading to substantial yield losses and compromised tea quality, posing a significant threat to the sustainable development of the tea industry (Zhang et al., 2024). Although chemical treatments provide short-term disease suppression, the resulting pesticide residues raise significant concerns regarding tea safety and environmental health (Sun et al., 2017). Thus, developing environmentally sustainable disease management strategies and enhancing the inherent resistance of tea plants are imperative for the long-term viability of the tea sector (Chen et al., 2018; Zhao et al., 2023).

Grafting is an established horticultural practice widely employed to improve plant resistance and performance across diverse crops (Goldschmidt, 2014; Kyriacou et al., 2017). In various species, grafting onto resistant rootstocks confers systemic disease resistance to susceptible scions through multiple mechanisms, including the long-distance transport of mobile signals such as hormones (e.g., salicylic acid and jasmonic acid), proteins, mRNAs, and small RNAs (including miRNAs and siRNAs) via the phloem (Guan et al., 2012; Jensen et al., 2012; Jeynes-Cupper and Catoni, 2023). These systemic signals are capable of reprogramming gene expression in the scion, promoting the accumulation of defense-related secondary metabolites, activating systemic acquired resistance, and inducing epigenetic modifications that collectively strengthen the plant’s innate immunity against biotic stress (Spanò et al., 2020). In tea cultivation, mounting evidence demonstrates that judicious selection of rootstocks can effectively enhance scion resistance, yield, and quality (Zhang, 1979; Wang, 2007; Karunakaran and Ilango, 2019). Notably, rootstock grafting can systemically reprogram the metabolic landscape of tea scions, leading to the enhanced accumulation of defense-related secondary metabolites, which are critical components of plant innate immunity against biotic stressors (Deng et al., 2017; Liu et al., 2025). This graft-mediated metabolic reprogramming underscores the potential of utilizing resistant rootstocks as a viable strategy to fortify scions against fungal pathogens such as C. camelliae.

Despite these advances, the precise long-distance signal transduction networks—especially the identity, mobility, and regulatory roles of rootstock-derived small RNAs and associated RNA-silencing machinery—in graft-induced resistance to C. camelliae remain largely unknown. In preliminary work, leveraging novel cultivars developed by the Guizhou Provincial Tea Research Institute, we established an effective grafting system pairing the C. camelliae susceptible cultivar ‘Qianmei 818’ as scion with the resistant cultivar ‘Qianmei 419’ as rootstock. Field trials revealed that grafting resulted in a pronounced reduction in C. camelliae incidence in ‘Qianmei 818’ scions while preserving their characteristic high epigallocatechin gallate (EGCG) content. This grafting system thus provides an ideal platform for dissecting the molecular signaling events underlying grafting-induced disease resistance in tea. To further elucidate these intricate molecular mechanisms, the present study implements a comprehensive multi-omics approach, integrating transcriptome (RNA-Seq), metabolome, and small RNA (sRNA-Seq) analyses. This integrated analysis aims to explore the intricate biological interactions and signal transduction pathways occurring between the scion and rootstock, thereby providing critical theoretical insights and practical guidance for enhancing resistance to C. camelliae in tea plants.

2 Results

2.1 Grafting strongly reduces lesion development caused by C. camelliae

The susceptibility of tea plants to C. camelliae was evaluated across different grafting combinations, revealing markedly distinct outcomes. The ungrafted susceptible scion cultivar ‘818’ developed the largest lesion diameter, averaging 11.67 ± 0.74 mm. In contrast, the resistant rootstock ‘419’ demonstrated strong resistance, with lesion diameters limited to 1.13 ± 0.07 mm. Self-grafted ‘818/818’ plants exhibited a significant 50.9% reduction in lesion size compared to ungrafted ‘818’. Remarkably, hetero-grafted ‘818/419’ plants showed an even greater decrease, with lesion diameters reduced by 98.7% relative to ungrafted ‘818’ and by 80.3% compared to ‘818/818’ (all differences statistically significant; p < 0.05) (Figures 1A, B). These findings robustly confirm that grafting substantially enhances resistance to C. camelliae in the susceptible cultivar, with hetero-grafting onto resistant ‘419’ rootstock conferring superior protective effects. This underscores the potential of resistant rootstocks to modulate scion disease resistance through grafting.

Figure 1
Panel A shows four types of plant grafts with varying stem and leaf colors, accompanied by root illustrations. Below, four leaves display varying lesion sizes labeled 818, 419, 818/818, and 818/419. Panel B depicts a bar graph comparing lesion diameters in millimeters for these labels, with 818 having the largest and 818/419 the smallest.

Figure 1. Effect of grafting on the susceptibility of tea plants to C. camelliae. (A) Symptoms 3 day after inoculation of leaves from four different grafting combinations (Susceptible cultivar ‘818’ alone; resistant cultivar ‘419’ alone; ‘818’ grafted on ‘818’; ‘818’ grafted on ‘419’). The 5-mm-diameter mycelial plugs used for inoculation were removed prior to the examination of the leaves. (B) Average lesion diameters 3 days after leaf inoculation. Different letters above the bars indicate significantly different lesion diameters according to Duncan’s multiple comparison test (p < 0.05). The error bars represent the standard error of the mean (n=30).

2.2 Grafting reshapes the production of defense enzymes, phytohormones, and tea quality-related metabolites

To elucidate the physiological and biochemical bases underpinning grafting-enhanced disease resistance, we systematically analyzed key components in the scion leaves of ‘818’, ‘419’, and grafted plants (‘818/818’ and ‘818/419’).

2.2.1 Defense enzyme activities

Phenylalanine ammonia-lyase (PAL), peroxidase (POD) and superoxide dismutase (SOD) are widely recognized as key biochemical indicators for evaluating disease resistance mechanisms in grafted plants (Singh et al., 2022). In our study, similar patterns were observed for the leaf content in the three enzymes, with statistically significant (p < 0.05) and markedly different values for heterografted ‘818/419’and for the three other plant groups (Figure 2A). Compared to ungrafted ‘818’, ungrafted ‘419’ and self-grafted ‘818/818’, the values for ‘818/419’ were 18.23-24.29% higher for PAL activity and they were 67.60-63.51% and 51.02-46.64% lower for POD and SOD activities, respectively. These contrasted enzyme activity profiles indicate that hetero-grafting with the resistant rootstock ‘419’ significantly elevates PAL activity but suppresses POD and SOD activity, potentially enhancing phenylpropanoid metabolism and contributing to improved defense responses in tea plants.

Figure 2
Nine bar graphs display various biochemical measurements in different grafting combinations of tea plants. Section (A) shows PAL, POD, and SOD enzyme activity; section (B) presents IAA, GA3, JA, SA, and MeSA levels; section (C) features EGCG, EC, ECG, EGC, L-Theanine, and Caffeine concentrations. Each graph uses different colored bars to indicate combinations, with labeled mean differences using letters above bars.

Figure 2. Physiological and biochemical effects of grafting on tea plants. (A) Activity of defense-related enzymes phenylalanine ammonia-lyase (PAL), peroxidase (POD) and superoxide dismutase (SOD) in tea leaves (in units per gram of leaf tissue). (B) Leaf tissue content in endogenous hormones Indole-3-acetic acid (IAA), gibberellin A3 (GA3), salicylic acid (SA), methyl salicylate (MeSA) and jasmonic acid (JA) (in µg per gram of leaf tissue). (C) Leaf tissue content in quality-related compounds epigallocatechin gallate (EGCG), epicatechin (EC), epigallocatechin (EGC), epicatechin gallate (ECG), L-theanine and caffeine (in mg per gram of leaf tissue). For a given compound, different letters above the bars indicate significant differences between grafting treatments according to Duncan’s multiple comparison test (p < 0.05). The error bars represent the standard error of the mean (n=3).

2.2.2 Phytohormone profiles

Plant hormones play crucial roles in regulating graft healing and defense responses. We focused on five phytohormones, Indole-3-acetic acid (IAA), gibberellin A3 (GA3), salicylic acid (SA), methyl salicylate (MeSA) and jasmonic acid (JA), because of their established roles in promoting growth, wound healing, and stress defense in grafted plants (Sharma and Zheng, 2019; Zhang et al., 2022; Wang et al., 2020a). Our analysis revealed that hetero-grafting with the resistant rootstock ‘419’ profoundly altered the hormonal landscape in the scion (Figure 2B). The hetero-graft ‘818/419’ exhibited significantly elevated levels of IAA and GA3, suggesting enhanced growth and wound-healing capacity. Regarding defense-related hormones, a notable shift was observed in the salicylate pathway: while basal SA levels were lower in ‘818/419’, its volatile derivative MeSA accumulated to the highest level. Concurrently, JA content was also reduced in the hetero-graft compared to the susceptible ‘818’ and self-grafted ‘818/818’. These results indicate that the resistant rootstock ‘419’ reprograms scion hormone homeostasis, promoting growth regulators (IAA, GA3) and channeling SA into a mobile defense signal (MeSA), while attenuating basal JA and SA pathways.

2.2.3 Quality-related metabolites

The cultivar ‘818’ is appreciated by tea consumers for its high epigallocatechin gallate (EGCG) content. We measured EGCG and other key quality-related compounds such as epicatechin (EC), epigallocatechin (EGC), epicatechin gallate (ECG), L-theanine and caffeine in ‘818’, ‘419’, and grafted plants. Cultivar ‘818’ exhibited significantly higher EGCG, EC, EGC, and L-theanine compared with ‘419’ (p < 0.05), whereas ‘419’ contained elevated ECG and caffeine. Self-grafting further increased EGCG, EC, and L-theanine but reduced ECG and caffeine in ‘818/818’. In hetero-grafted ‘818/419’ plants, the EGCG content was similar to levels of ungrafted ‘818’, while those of EC and EGC were significantly higher than in all other groups (p < 0.05). L-theanine also increased in ‘818/419’ compared to the ungrafted groups (p < 0.05). Meanwhile, ECG and caffeine exceeded both ‘818’ and ‘818/818’ levels (p < 0.05) (Figure 2C). These results indicate that hetero-grafting integrates key quality traits from both rootstock and scion, synergistically enhancing levels of EC, EGC, and L-theanine while preserving EGCG and enriching ECG and caffeine, thereby improving overall tea quality.

2.3 Grafting induces transcriptomic remodeling in the scion

To explore molecular mechanisms of grafting-enhanced disease resistance, transcriptome sequencing was performed on scions from ungrafted ‘818’, self-grafted ‘818/818’, and hetero-grafted ‘818/419’. High-quality data totaled 56.77 Gb, with individual samples yielding ≥5.96 Gb clean reads, Q30>96.6%, and genome mapping rates >85%. Differential expression analysis identified 229 differentially expressed genes (DEGs; 75 upregulated, 154 downregulated) in ‘818/818’ vs ‘818’; 565 DEGs (155 upregulated, 410 downregulated) in ‘818/419’ vs ‘818/818’; and 1205 DEGs (299 upregulated, 906 downregulated) in ‘818/419’ vs ‘818’ (Supplementary Figure S1A), indicating stronger regulatory effects in plants grafted on resistant rootstocks.

Gene ontology (GO) enrichment revealed that grafting primarily modulated genes involved in metabolic and cellular processes, cellular components, catalytic activity, and binding, whereas resistant rootstocks broadened activation to include metabolic processes, signal transduction, and transcriptional regulation (Figures 3A, B). Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed DEGs in ‘818 vs 818/818’ enriched pathways such as starch/sucrose metabolism, galactose metabolism, unsaturated fatty acid biosynthesis, phenylpropanoid biosynthesis, and peroxisome function (Figures 3C, D). In contrast, ‘818/818 vs 818/419’ DEGs were enriched in glutathione, nitrogen, sulfur metabolism, linoleic acid metabolism, pentose phosphate pathway, carbon metabolism, photosynthesis, and stress-response pathways including chloroplast light-harvesting complex assembly and cysteine/methionine metabolism. Thus, grafting mainly induced metabolic remodeling and cellular adaptation facilitating grafting stress tolerance and transport, while resistant rootstocks potentiated defense-related metabolic and signaling pathways in the scion.

Figure 3
Four-panel data visualization comparing gene ontology (GO) and KEGG pathways between different samples. Panel A shows GO biological processes, cellular components, and molecular functions with upregulated and downregulated genes. Panel B displays similar GO categories for a different comparison. Panel C illustrates KEGG pathways for various metabolic processes, using color gradients to indicate q-values. Panel D provides another KEGG comparison with pathways such as glutathione metabolism and photosynthesis, with a similar q-value gradient. Each panel uses distinct colors to differentiate between upregulated (blue/orange) and downregulated (green) genes.

Figure 3. Functional enrichment of DEGs from ‘818 vs 818/818’ and ‘818/818 vs 818/419’. (A, B) Gene Ontology (GO) enrichment; (C, D) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment.

2.4 Grafting modulates miRNA expression and target pathways

Small RNA sequencing of scions from ‘818’, ‘818/818’, and ‘818/419’ yielded 109.68 million clean reads, with Q30 > 85%, and sequence lengths predominantly 21–24 nt, characteristic of plant miRNAs. A total of 615 miRNAs were identified (88 known, 527 novel). Differential expression analysis revealed 126 differentially expressed miRNA (DEMs;48 upregulated, 78 downregulated) in ‘818 vs 818/818’ and 31 DEMs (15 upregulated, 16 downregulated) in ‘818/818 vs 818/419’ (Supplementary Figure S1B). KEGG enrichment of DEM targets indicated involvement in plant hormone signaling, pathogen interaction, homologous recombination, and mRNA surveillance in ‘818 vs 818/818’ (Figure 4A); and additionally sulfur metabolism, steroid biosynthesis, diarylheptanoid/gingerol biosynthesis, endocytosis, and photosynthesis regulation in ‘818/818 vs 818/419’ (Figure 4B). These data suggest that grafting mediates physiological adaptation and defense enhancement via miRNA-regulated networks governing defense, tissue repair, and metabolism.

Figure 4
Two bar charts labeled A and B display KEGG pathway analyses. Chart A compares pathways with gene numbers ranging from zero to twenty-five, highlighting notable ones like “plant-pathogen interaction” and “endocytosis.” Chart B displays gene numbers up to ninety, emphasizing “plant hormone signal transduction” and “plant-pathogen interaction.” Both charts use color gradients representing q-values from red (0.2) to blue (0.6 or 0.75).

Figure 4. KEGG enrichment of miRNA target genes in different tea grafting combinations. (A) Comparison between ungrafted ‘818’ and self-grafted ‘818/818’; (B) Comparison between self-grafted ‘818/818’ and hetero-grafted ‘818/419’. Color represents adjusted q-value.

2.5 Metabolomic reprogramming associated with enhanced resistance

Widely targeted metabolomics analysis detected 1870 metabolites in scion leaves. Differential analysis identified 430 differential metabolites (261 upregulated, 169 downregulated) in ‘818 vs 818/818’, and 361 (171 upregulated, 190 downregulated) in ‘818/818 vs 818/419’ (Supplementary Figure S1C). KEGG enrichment showed upregulated metabolites in ‘818 vs 818/818’ enriched in monoterpenoid biosynthesis, phenylalanine metabolism, and ubiquinone/terpenoid-quinone biosynthesis; downregulated ones were enriched in carboxylic acid metabolism, the glyoxylate cycle, glycerophospholipid metabolism, and beta-alanine metabolism (Figure 5A). In the ‘818/818 vs 818/419’ comparison, upregulated metabolites were enriched in glycerophospholipid metabolism, caffeine metabolism, and phosphatidylinositol signaling; downregulated metabolites were involved in acetyl-CoA and lipid pathways (Figure 5B). These results indicate that grafting modulated primary energy metabolism and enhanced accumulation of secondary metabolites related to quality and disease resistance. Resistant rootstock ‘419’ further amplified disease-relevant metabolic pathways, promoting biosynthesis and accumulation of defense compounds that improve resistance and quality.

Figure 5
(A) Dot plot showing KEGG pathways with differential abundance scores for 818 versus 8181818. (B) Similar dot plot for 8181818 versus 8181419. (C) Venn diagram comparing metabolites and genes for 818 versus 8181818 shows overlapping and unique elements. (D) Bar chart of KEGG pathways, highlighting gene and metabolite counts for 818 versus 8181818. (E) Venn diagram for 8181818 versus 8181419 comparing genes and metabolites. (F) Bar chart of top ten KEGG pathways for 8181818 versus 8181419, indicating gene and metabolite numbers.

Figure 5. Integrated KEGG pathway enrichment analysis of differential metabolites and co-expressed genes in different tea grafting combinations. (A, B) Visualization of KEGG pathway enrichment based on differential metabolites in ‘818 vs 818/818’and ‘818/818 vs 818/419’. For each pathway, the horizontal position of the dot represents the Differential Abundance (DA) Score (ranging from -1 to 1), indicating the overall up- or down-regulation trend of all annotated metabolites within that pathway. The distance from the center line reflects the absolute value of the DA Score. Dot size corresponds to the number of differential metabolites mapped to the pathway. Dot and line color reflects the statistical significance (p-value), with a gradient from blue (less significant) to red (more significant); (C, E) Venn diagrams illustrating overlapping KEGG-enriched pathways between ‘818 vs 818/818’ and ‘818/818 vs 818/419’ for differential metabolites and co-expressed differential genes and metabolites; (D, F) Top 10 KEGG-enriched pathways for co-expressed differential genes and metabolites in ‘818 vs 818/818’and ‘818/818 vs 818/419’. Bar height indicates pathway activity level within the measured samples.

2.6 Integrated multi-omics reveals coordinated defense networks

Integrative analysis of transcriptomic and metabolomic data revealed that co-expressed DEGs and metabolites in the comparison ‘818 vs 818/818’ were enriched in 34 KEGG pathways, predominantly starch and sucrose metabolism, phenylalanine metabolism, carbon metabolism, and isoflavonoid biosynthesis (Figures 5C, D). In ‘818/818 vs 818/419’, 49 enriched pathways were identified, including plant hormone signal transduction, carbon metabolism, glutathione metabolism, amino acid biosynthesis, and photosynthesis regulation (Figures 5E, F). These results suggest that grafting modulates basal metabolism and secondary metabolite synthesis, while resistant rootstocks activate comprehensive defense-related gene expression and metabolite accumulation, enhancing the plant’s immune response. Complementary miRNA sequencing data corroborated the observed enrichment of target genes in hormone signaling and glutathione metabolism pathways, underscoring the miRNA-mediated post-transcriptional regulation coordinating these responses.

2.7 qRT-PCR validation of miRNA-mRNA pairs

To verify the sequencing data, we selected seven differentially expressed miRNAs and ten associated target genes for qRT-PCR analysis. The expression patterns were largely consistent with the high-throughput sequencing results, revealing predominant inverse correlations between the miRNAs and their putative targets (Figure 6; Supplementary Tables S2, S3). All seven validated miRNAs showed significant downregulation in the resistant hetero-graft ‘818/419’ compared to the susceptible ungrafted ‘818’ (p < 0.05). Among them, csi-miR395b-3p, novel-miR397, nta-miR160a, mdm-miR171h, novel-miR170, and vvi-miR535c were also downregulated in self-grafted ‘818/818’ and the resistant rootstock ‘419’, whereas vvi-miR159c was specifically suppressed in ‘818/419’.

Figure 6
Bar charts illustrating the relative expression levels of various miRNAs and genes in tea plant grafting combinations. Section A shows miRNA expression (csi-miR395b-3p, mdm-miR171h, novel-miR397, etc.) across different combinations, with a flow diagram linking each miRNA to specific genes (like CsAPS1, CsSCL6). Section B displays charts for gene expression (CsAPS1, CsCCoAOMT, CsCCR2, etc.). Different grafting combinations are represented, with varied expression levels indicated by bars and statistically significant differences marked by letters above each bar.

Figure 6. Experimental validation of key omics findings by qRT-PCR. (A) Relative expression levels of selected differentially expressed miRNAs; (B) Relative expression levels of selected differentially expressed genes. Expression levels were validated in the relevant tea grafting combinations. Different letters above the bars indicate significantly different lesion diameters according to Duncan’s multiple comparison test (p < 0.05).

Conversely, the expression of most target genes was significantly upregulated. For example, downregulation of csi-miR395b-3p was associated with strong induction of its target involved in sulfur assimilation, CsAPS1, in both ‘419’ and ‘818/419’. Suppression of novel-miR397 correlated with marked upregulation of three lignin biosynthesis-related genes—CsCCoAOMT, CsCCR2, and CsDLO2—in the resistant and hetero-grafted plants. Downregulation of nta-miR160a coincided with elevated levels of CsLAX2 (an auxin transporter) and CsCYP78A9, indicating enhanced auxin signaling. Similarly, downregulation of vvi-miR159c and mdm-miR171h aligned with upregulation of their respective putative targets: the transcription factor CsMYB30 and the signaling gene CsSCL6.

Notably, the self-grafted ‘818/818’ typically displayed intermediate expression changes, whereas hetero-grafting with ‘419’ elicited the most substantial transcriptional reprogramming. These validated miRNA-mRNA pairs constitute a core post-transcriptional regulatory network that underlies the grafting-enhanced defense response.

3 Discussion

Grafting technology, as an effective strategy to enhance plant resistance, has shown substantial application value across various crops and economic plants (Goldschmidt, 2014; Kyriacou et al., 2017). In this study, we elucidated how the resistant rootstock ‘Qianmei 419’ (‘419’) improves resistance to C. camelliae in the susceptible tea scion ‘Qianmei 818’ (‘818’) by integrating transcriptomic, metabolomic, and small RNA omics analyses. Our study showed that disease severity (assessed by the average lesion diameter of inoculated leaves) was 98.71% lower in hetero-grafted ‘818/419’ plants than is ungrafted ‘818’. However, the complex molecular mechanisms underlying this grafting-conferred resistance required further clarification. By providing combined physiological, biochemical, and multi-omics data, this work reveals a coordinated network involving hormone signaling, miRNA regulation, and metabolic reprogramming during grafting, thereby providing a theoretical framework and practical insights for breeding disease-resistant tea cultivars.

3.1 Hormonal and enzymatic bases of graft enhanced resistance to C. camelliae

Grafting markedly altered the scion’s hormonal profile, especially in the heterograft ‘818/419’, where IAA and GA3 contents increased significantly. These growth-promoting hormones likely support wound healing and vigor, thereby contributing to enhanced resistance (Sharma and Zheng, 2019; Zhang et al., 2022). In parallel, MeSA, a mobile derivative of SA associated with defense priming (Lan et al., 2025), also accumulated strongly in both self- and hetero-grafted plants, suggesting that grafting reprograms SA metabolism toward enhanced defensive readiness. At the enzymatic level, heterografting with the resistant rootstock elevated PAL activity, a key enzyme for lignin and flavonoid biosynthesis (Dixon et al., 2002), while reducing SOD and POD activities, indicating a shift from broad antioxidant responses toward specialized defense pathways. Taken together, these findings suggest that grafting, particularly onto resistant rootstocks, enhances tea plant resistance to C. camelliae by promoting IAA- and GA3-mediated growth, strengthening SA-derived MeSA defense signaling, and activating PAL-driven secondary metabolism.

3.2 Resistant rootstocks mediate C. camelliae resistance through miRNA-regulated gene expression

The miRNA-mRNA expression patterns described above provide a mechanistic basis for understanding how resistant rootstocks influence scion defense responses. Building on these results, small RNA sequencing and qPCR validation further revealed that grafting, particularly hetero-grafting with resistant rootstock ‘419’, induces pronounced shifts in scion miRNA expression profiles, accompanied by reciprocal changes in their target or putative target genes. For example, csi-miR395b-3p was significantly downregulated in ‘419’ and ‘818/419’, while its target CsAPS1 and related sulfur metabolism gene CsLCD were upregulated, a result consistent with the known role of miR395’s in sulfur homeostasis (Kawashima et al., 2009; Tao et al., 2012). Enhanced sulfur metabolism is known to improve antioxidant capacity via glutathione synthesis and production of sulfur-containing defense compounds (Bednarek et al., 2009; Hacham et al., 2025).

Our study resulted in the identification of a novel miRNA, that we propose to name ‘novel-miR397’ (Supplementary Table S4). The downregulation of this novel-miR397 correlated with increased expression of three genes, CsCCoAOMT, CsCCR2, and CsDLO2, that are known to be linked with reinforcing lignin biosynthesis and structural barriers (Huang et al., 2021; Xue et al., 2019; Wang et al., 2019). Reduced nta-miR160a was linked to elevated CsLAX2 and CsCYP78A9, indicating enhanced auxin transport and organ development (Hao et al., 2022; Moreno-Piovano et al., 2017; Sotelo-Silveira et al., 2013). Downregulation of vvi-miR159c corresponded with higher expression of CsMYB30, a gene known to code for a stress-responsive transcription factor (Millar et al., 2019; Fichman et al., 2020), while reduced mdm-miR171h matched increased expression of CsSCL6, a gene involved in hormone signaling (Wang et al., 2020b; Yang et al., 2010). These coordinated miRNA–mRNA shifts, with intermediate changes in self-grafted ‘818/818’, support a systemic, miRNA-driven regulatory mechanism in graft-mediated resistance, a finding consistent with the role proposed by Melnyk et al. (2011b) for long-distance miRNA transport.

3.3 Graft-induced transcriptional and metabolic reprogramming enhances resistance while preserving scion quality traits

Integrated metabolomic and transcriptomic analyses demonstrated that grafting, especially hetero-grafting with resistant rootstock ‘419’, reprograms scion metabolism and gene expression. In ‘818/419’, upregulated metabolites were enriched in glycerophospholipid metabolism, phosphatidylinositol signaling, and caffeine biosynthesis, processes linked in previous studies to membrane stability, immune signal transduction, and pathogen inhibition (Guo et al., 2023; Abd-El-Haliem and Joosten, 2017; Zhang et al., 2010). Integrated metabolomic and transcriptomic analyses revealed that in ‘818/419’ scion leaves, differential metabolites and genes in pathways such as plant hormone signal transduction, glutathione metabolism, and phenylpropanoid biosynthesis exhibited synergistic expression, demonstrating a tight transcriptional-metabolic coupling. This further demonstrates that the increased expression of key genes such as CsAPS1, CsCCoAOMT, CsCCR2, CsDLO2, CsLAX2 and CsLCD promotes plant hormone and lignin synthesis, as well as signal transduction. These changes parallel the downregulation of their regulatory miRNAs, reinforcing the link between post-transcriptional control and metabolic adaptation.

Notably, the defense-related metabolic reprogramming induced by the resistant rootstock did not compromise the key quality-determining metabolites of the scion. ‘818/419’ was shown in our study to hold important traits associated with tea quality. It retained the characteristic high EGCG content of scion ‘818’, while showing increased EC, EGC, and L-theanine. Furthermore, it preserved two desirable traits of rootstock ‘419’: high ECG and caffeine contents. These findings demonstrate that grafting onto a resistant rootstock triggers extensive miRNA-directed regulatory cascades, broad transcriptional reprogramming, and differential enrichment of defense-associated pathways, thereby conferring robust phylloxera resistance while maintaining the scion’ s canonical levels of key quality-associated metabolites and permitting the accumulation of additional compounds contributed by the rootstock.

3.4 An integrated model of grafting-enhanced resistance to C. camelliae

This study integrates multi-omics and phenotypic data to propose a hierarchical, systemic resistance network orchestrated by grafting-responsive miRNAs. The enhanced resistance in ‘818/419’ scions is a coherent cascade initiated by rootstock-derived signals (Figure 7). Specifically, the differential expression of key miRNAs (e.g., csi-miR395b-3p, novel-miR397) in the scion fine-tunes the transcriptome, leading to the upregulation of target genes governing sulfur assimilation (CsAPS1), lignin biosynthesis (CsCCoAOMT, CsCCR2), and auxin transport (CsLAX2). This transcriptional reprogramming subsequently drives a metabolic shift toward defense-ready states, characterized by the accumulation of MeSA, IAA for signaling, enhanced lignin deposition for structural reinforcement, and elevated glutathione for redox homeostasis. Consequently, the tea plant’s strategy pivots from a broad antioxidant response (reflected in altered SOD/POD activities) to a targeted deployment of specialized physical and biochemical barriers, which directly translates into the dramatic 98.71% reduction in lesion diameter.

Figure 7
Diagram illustrating a grafting system. On the left, a plant labeled “Qianmei 818” (susceptible) is grafted onto “Qianmei 419” (resistant). The center panel details molecular mechanisms, showing miRNA downregulation and target gene upregulation like CsAPS1 and CsMYB30. Systemic signaling leads to hormonal changes, affecting pathways and tea quality. The right panel shows phenotypic outcomes, with “Colletotrichum camelliae” causing a large lesion on ungrafted leaves, while grafted leaves show enhanced resistance with a small lesion.

Figure 7. Multi-omics model of grafting enhanced resistance to C. camelliae in tea plants. Grafting the susceptible scion ‘Qianmei 818’ onto the resistant rootstock ‘Qianmei 419’ triggers systemic signaling that reprograms the scion. (1) Key miRNAs shift (e.g., downregulation of csi-miR395b-3p and novel-miR397); (2) Putative target genes are upregulated, activating defense pathways including sulfur metabolism (CsAPS1, CsLCD), lignin biosynthesis (CsCCoAOMT, CsCCR2, CsDLO2), and hormone signaling/transport (CsLAX2, CsCYP78A9); (3) Phytohormone homeostasis shifts, with increased IAA, GA3, and MeSA, and decreased JA and SA; (4) Defense enzyme activities are modulated (higher PAL; lower POD and SOD); (5) The metabolome is remodeled, enriching glycerophospholipid metabolism and phenylpropanoid biosynthesis while preserving or enhancing quality-related metabolites (EGCG, caffeine, L-theanine). Together, these changes strengthen structural barriers (lignin), augment chemical defenses, and prime immunity, resulting in markedly reduced lesion development following C. camelliae infection.

Grafting-induced resistance conferred by resistant rootstocks is well-documented in various crops, including grapevine (Cookson et al., 2013), apple (Jensen et al., 2012), watermelon (Zhang et al., 2025), and citrus (Wang et al., 2025). The enhanced resistance in these grafted systems is primarily attributed to long-distance signaling mechanisms. In grafting systems, mobile signals such as miRNAs, hormones, and metabolites are likely transported from rootstock to scion through the phloem via source-sink dynamics, thereby reprogramming scion responses. In this study, we postulate that similar long-distance signaling enhances scion defense against C. camelliae in grafted tea plants. Numerous studies have established that 21–24 nt small RNAs, including many of the miRNAs identified here, can function as mobile silencing signals that traverse graft junctions in a selective, Argonaute-dependent manner (Melnyk et al., 2011a; Li et al., 2021). Once in the recipient scion tissue, these mobile miRNAs-transported as single-or double-stranded forms bound to Argonaute proteins-trigger target mRNA cleavage, translational repression, or epigenetic modifications such as DNA methylation. This process is frequently amplified by scion-localized RNA-dependent RNA polymerases and mirrors endogenous long-distance signaling events, such as the well-characterized shoot-to-root movement of miR399 during phosphate homeostasis (Lin et al., 2008; Pant et al., 2008). The observed changes in miRNA expression levels in the scions of grafted tea plants may therefore be associated with the mobility of miRNAs between rootstock and scion. Additionally, hormones such as MeSA, which accumulated in our hetero-grafts, represent another class of mobile signals, as volatile or phloem-transported SA derivatives, they prime systemic acquired resistance across distances (Park et al., 2007). Similarly, defense-related metabolites (e.g., from phenylpropanoid and glutathione pathways) may be conveyed with photoassimilates via the phloem, thereby enhancing scion defense compounds. Currently, our data provide indirect evidence for the long-distance transport of these substances in grafted tea plants. Future studies employing tracer assays or direct phloem profiling are needed to validate these transport mechanisms. Notably, while our multi-omics approach provides robust correlative support for miRNA-mRNA networks, functional validations such as miRNA or gene overexpression and silencing are crucial to establish causality and delineate the mechanistic details of grafting-conferred immunity.

4 Conclusion

Through an integrative multi-omics approach, this study delineates the molecular mechanisms underlying grafting-enhanced C. camelliae resistance in tea plants. The resistant rootstock ‘419’ remodels hormonal homeostasis, modulates key miRNAs and their targets, and reprograms metabolic pathways, thereby significantly improving disease resistance and simultaneously enhancing tea quality in the susceptible scion ‘818’. These insights provide a theoretical foundation for disease-resistant breeding in tea and support the advancement of sustainable, environmentally friendly disease control strategies. Future work will focus on elucidating the mechanisms governing long-distance miRNA transport between rootstock and scion and evaluating grafted tea plant resistance under complex field conditions to optimize practical applications of grafting technology in tea cultivation.

5 Materials and methods

5.1 Plant materials and grafting scheme

Two-year-old tea plants (Camellia sinensis) cultivated by the Guizhou Provincial Tea Research Institute were used in this study. The C. camelliae susceptible cultivar ‘Qianmei 818’ (hereafter ‘818’) served as the scion, while the C. camelliae resistant cultivar ‘Qianmei 419’ (hereafter ‘419’) was employed as the rootstock for grafting. Grafting was conducted in spring using the cleft grafting method described by Liu et al. (2017). Two graft combinations were established: self-grafted (‘818/818’; i.e., ‘818’ scion onto ‘818’ rootstock) and hetero-grafted (‘818/419’; i.e., ‘818’ scion onto ‘419’ rootstock). Following grafting, seedlings underwent standard field management and were cultivated for one year to ensure stable growth prior to downstream experiments. This one-year period was chosen to ensure complete graft union formation, stable systemic signal transmission between rootstock and scion, and consistent physiological maturity prior to inoculation experiments, as preliminary field trials indicated that shorter periods resulted in incomplete healing and variable resistance phenotypes.

5.2 Detached leaf inoculation assay for C. camelliae

The resistance of grafted tea plants to C. camelliae was assessed using the detached leaf inoculation assay of Yoshida and Takeda (2006). The first fully expanded leaves (third leaf from the apex) were excised from healthy plants belonging to four groups: ungrafted ‘818’, ungrafted ‘419’, self-grafted ‘818/818’, and hetero-grafted ‘818/419’. The leaves were surface-sterilized by wiping with 75% ethanol and rinsed with sterile distilled water. They were blotted dry, then placed in 90-mm-diameter Petri dishes on water-soaked sterile filter paper, with their adaxial side up. The leaves were then immediately inoculated with C. camelliae stain GC-2. This strain was originally isolated from anthracnose-infected tea leaves in Guizhou Province, morphologically and phylogenetically identified as C. camelliae based on multi-gene (ITS, GAPDH, ACT, CHS) analysis combined with pathogenicity tests (as detailed in Chen, 2023, Master’s thesis, Guizhou University), and preserved in the Plant Protection Laboratory, Guizhou institute of tea science. Long-term stock cultures are maintained in 20% glycerol at -80°C, with working cultures stored on potato dextrose agar (PDA) slants at 4°C. Mycelial plugs, 5 mm in diameter, were excised from the growing margin of a 7-day old colony incubated on potato dextrose agar (PDA) at 25°C in the dark. On each leaf, a mycelial plug was deposited on the left side of the midrib, with the mycelium in direct contact with the plant tissue. A plug of non-inoculated PDA was placed on the right side of the midrib as a control. The Petri dishes were then covered with moist cotton, and incubated at 28°C in darkness. The moisture of the cotton was monitored during incubation and water was added as needed to maintain saturated humidity around the leaves. Disease development was assessed by measuring the diameter of the lesions 3 days after inoculation. To facilitate the observation of lesions with low development (diameters smaller than 5 mm) the mycelial plugs were removed from the leaves prior to their examination. Each treatment included 10 replicate leaves, and the whole experiment was carried out three times independently.

5.3 Physiological parameter measurement

Healthy shoots comprising one bud and two leaves were excised from plants with consistent growth status in all four treatment groups (‘818’, ‘419’, ‘818/818’, and ‘818/419’). The shoots were immediately flash-frozen in liquid nitrogen and stored at -80°C until analyses were carried out. The activity of defense-related enzymes peroxidase (POD), superoxide dismutase (SOD) and phenylalanine ammonia-lyase (PAL) was quantified using commercial colorimetric assay kits (Suzhou Kaiming, China) following the manufacturer’s instructions.

The levels of plant hormones including indole-3-acetic acid (IAA), gibberellic acid (GA3), salicylic acid (SA), and jasmonic acid (JA), alongside major catechins (epigallocatechin gallate [EGCG], epicatechin [EC], epigallocatechin [EGC], and epicatechin gallate [ECG]), L-theanine and caffeine, were determined via high-performance liquid chromatography (HPLC). Metabolites were extracted using methanol, and quantification was performed based on calibration curves generated from purified standards, with relative standard deviations maintained below 15%. Concentrations (µg/g) were calculated as follows:

Concentration = (Instrument Reading x Final Extract Volume [mL] x Dilution Factor)/Sample Mass [g].

For each of the four plant treatment groups, we analyzed three biological replicate samples, each sample being composed of shoots from five plants.

5.4 Multi-omics data analysis

5.4.1 Sample collection, RNA, and metabolite extraction

Samples of shoots comprising one bud and the two first leaves were collected from plants of each of three plant groups (‘818’ ungrafted, ‘818/818’ self-grafted, and ‘818/419’ hetero-grafted), immediately flash-frozen in liquid nitrogen, and stored at -80°C. Total RNA extraction was performed using a Plant Total RNA Extraction Kit (Tiangen Biotech Co., Ltd., Beijing, China), following the manufacturer’s protocol. RNA quality and concentration were evaluated using a NanoDrop 2000 spectrophotometer and Agilent 2100 Bioanalyzer. Only high-quality RNA samples were subjected to subsequent mRNA and micro RNA (miRNA) sequencing. Metabolites were extracted using a methanol: acetonitrile: water solvent mixture (1:2:1, v/v/v). Ultimately, each group was represented by three biological replicates, each containing pooled material from five plants.

5.4.2 RNA-Seq and sRNA-Seq analysis

Separate mRNA and sRNA libraries were constructed and sequenced on the Illumina platform. Clean mRNA reads were aligned to the tea plant reference genome Shuchazao V2 (Xia et al., 2020) with HISAT2 v2.0.4 (-{{-}}-dta -p 6 -{{-}}-max-intronlen 5000000), and transcripts were assembled using StringTie v2.2.1. Differentially expressed genes (DEGs) were identified by DESeq2 (|log2FoldChange| > 1.5, p < 0.05, Benjamini-Hochberg FDR). Functional annotation was performed against Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases.

For small RNA data, reads were mapped with Bowtie (-v 0 -S -f) for known miRNAs and miRDeep2 (-g 50000 -l 250 -d -m 10 -v -P -n d) for novel miRNAs. Differentially expressed miRNAs (DEMs) were identified using DESeq2 (|log2FoldChange| > 0.58, p < 0.05, FDR). Target genes of DEMs were predicted using TargetFinder v1.6, and annotated via BLAST searches against public databases, followed by KEGG pathway enrichment analyses.

5.4.3 Widely targeted metabolomics analysis

Metabolite profiling was conducted using a widely targeted metabolomics approach on an Acquity I-Class PLUS ultra-performance liquid chromatography (UPLC) system coupled with an AB Sciex QTRAP 6500+ mass spectrometer. Metabolite identification and quantification were performed against the Biomarker Technologies proprietary database GB-PLANT supplemented by public databases (KEGG Compound, HMDB, and Lipidmaps). Raw data were processed using Analyst v1.7.2 for peak integration, calibration, and normalization. Differentially accumulated metabolites (DAMs) were identified via OPLS-DA with thresholds of Variable Importance in Projection (VIP) scores ≥ 1,fold change>1 and p < 0.05. Metabolites were subsequently annotated and mapped to metabolic pathways using the KEGG database. Integrated multi-omics analyses combining widely targeted metabolomics with transcriptomic and small RNA datasets were performed to elucidate the molecular regulatory networks underpinning grafting-induced anthracnose resistance.

5.5 qRT-PCR validation

To validate RNA-Seq and sRNA-Seq findings, representative differentially expressed miRNAs and their target mRNAs were selected for quantitative real-time PCR (qRT-PCR). miRNA cDNA synthesis employed stem-loop primers, whereas mRNA was reverse transcribed via a universal reverse transcription kit (TAKARA Bio Inc., Japan). qRT-PCR was performed using SYBR Green on an ABI 7500 Real-Time PCR System. Pc-222-3p and CsGAPDH served as endogenous controls for miRNA and mRNA quantification, respectively. Relative expression was calculated using the 2–ΔΔCt method. Each reaction included three biological replicates and three technical replicates. Primer sequences are listed in Supplementary Table S1.

5.6 Statistical analysis and data integration

Data were analyzed by one-way analysis of variance (ANOVA) using SPSS 22.0. Significant differences between groups were determined by Duncan’s multiple range test at a significance threshold of p < 0.05.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions

Y-XL: Writing – original draft, Writing – review & editing, Data curation, Formal analysis, Funding acquisition, Project administration. K-QL: Writing – original draft, Funding acquisition, Project administration. A-RW: Data curation, Formal Analysis, Investigation, Writing – review & editing. JW: Data curation, Formal analysis, Investigation, Writing – review & editing. D-HY: Investigation, Project administration, Supervision, Writing – review & editing. D-GZ: Funding acquisition, Project administration, Supervision, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the Guizhou Province High-level Innovative Talent Training Program Project ([2016]4003), the Guizhou Academy of Agricultural Sciences Talent Special Project (Grant No. 2023-02), the Guizhou Province Science and Technology Plan Project (Qian Ke He Ji Chu QN[2025]232), the Qian Nong Ke Doctoral Fund ([2025]06), as well as the GZMARS-Tea Project (GZYCYJSTX-02).

Acknowledgments

Sequencing in this study was performed by Biomarker Technologies Co., Ltd. (Beijing, China).

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer WJ declared a past co-authorship with the author(s) D-GZ to the handling editor.

Generative AI statement

The author(s) declared that generative AI was used in the creation of this manuscript. Generative AI was used in the preparation of this manuscript solely for the purpose of language polishing and refinement. The AI tool (e.g., ChatGPT) was employed to assist with grammar checking, sentence structure improvement, and ensuring fluency in English. All scientific content, data analysis, interpretation of results, and conclusions are solely the responsibility of the authors. The authors have thoroughly reviewed and edited all AI-generated suggestions and take full responsibility for the final content of the manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

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

Supplementary material

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

References

Abd-El-Haliem, A. M. and Joosten, M. H. A. J. (2017). Plant phosphatidylinositol-specific phospholipase C at the center of plant innate immunity. J. Integr. Plant Biol. 59, 164–179. doi: 10.1111/jipb.12520

PubMed Abstract | Crossref Full Text | Google Scholar

Bednarek, P., Piślewska-Bednarek, M., Svatoš, A., Schneider, B., Doubský, J., Mansurova, M., et al. (2009). A glucosinolate metabolism pathway in living plant cells mediates broad-spectrum antifungal defense. Science 323, 101–106. doi: 10.1126/science.1163732

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, W. (2023). Activity of three plant-derived volatile components against Colletotrichum camelliae in tea plants (Guiyang: Guizhou University).

Google Scholar

Chen, S., Ren, Y. F., Li, D. X., Wang, Y., Jin, L. H., and Tan, X. F. (2018). Research and application pathways of sustainable control techniques for tea diseases in Guizhou. China Plant Prot. Guide. 38, 85–90.

Google Scholar

Cookson, S. J. and Nathalie, O. (2013). Grafting with rootstocks induces extensive transcriptional re-programming in the shoot apical meristem of grapevine. BMC Plant Biol. 13, 147–161. doi: 10.1186/1471-2229-13-147

PubMed Abstract | Crossref Full Text | Google Scholar

Deng, W. W., Fan, Y. B., Gu, C. C., Li, D. X., and Wan, X. C. (2017). Changes in morphological characters and secondary metabolite contents in leaves of grafting seedlings with Camellia sinensis as scions and C. oleifera as Stocks. J. Trop. Subtropical Bot. 25, 35–42. doi: 10.11926/jtsb.3641

Crossref Full Text | Google Scholar

Dixon, R. A., Achnine, L., Kota, P., Liu, C. J., Reddy, M. S. S., and Wang, L. (2002). The phenylpropanoid pathway and plant defense—a genomics perspective. Mol. Plant Pathol. 3, 371–390. doi: 10.1046/j.1364-3703.2002.00131.x

PubMed Abstract | Crossref Full Text | Google Scholar

Fichman, Y., Sara, I. Z., Soham, S., David, J. B., Ronald, J. M., and Rajeev, K. A. (2020). MYB30 orchestrates systemic reactive oxygen signaling and plant acclimation. Plant Physiol. 184, 666–675. doi: 10.1104/pp.20.00859

PubMed Abstract | Crossref Full Text | Google Scholar

Goldschmidt, E. E. (2014). Plant grafting: new mechanisms, evolutionary implications. Front. Plant Sci. 5. doi: 10.3389/fpls.2014.00727

PubMed Abstract | Crossref Full Text | Google Scholar

Guan, W., Zhao, X., Hassell., R., and Thies., J. (2012). Defense mechanisms involved in disease resistance of grafted vegetables. HortScience 47, 164–170. doi: 10.21273/HORTSCI.47.2.164

Crossref Full Text | Google Scholar

Guo, A., Yang, Y., Wu, J., Qin, N., Hou, F., and Gao, Y. (2023). Lipidomic and transcriptomic profiles of glycerophospholipid metabolism during Hemerocallis citrina Baroni flowering. BMC Plant Biol. 23, 1–15. doi: 10.1186/s12870-022-04020-x

PubMed Abstract | Crossref Full Text | Google Scholar

Hacham, Y., Kaplan, A., Cohen, E., Gal, M., and Amir, R. (2025). Sulfur metabolism under stress: Oxidized glutathione inhibits methionine biosynthesis by destabilizing the enzyme cystathionine γ-synthase. J. Integr. Plant Biol. 67, 87–100. doi: 10.1111/jipb.13799

PubMed Abstract | Crossref Full Text | Google Scholar

Hao, K., Wang, Y., Zhu, Z., Wu, Y., Chen, R., and Zhang, L. (2022). miR160: An indispensable regulator in plant. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.833322

PubMed Abstract | Crossref Full Text | Google Scholar

Huang, S., Zhou, J., Gao, L., and Tang, Y. (2021). Plant miR397 and its functions. Funct. Plant Biol. 48, 361–370. doi: 10.1071/fp20342

PubMed Abstract | Crossref Full Text | Google Scholar

Jensen, P., Halbrendt, N., Fazio, G., Makalowska, I., Altman, N., and Praul, C. (2012). Rootstock-regulated gene expression patterns associated with fire blight resistance in apple. BMC Genomics 13, 9–26. doi: 10.1186/1471-2164-13-9

PubMed Abstract | Crossref Full Text | Google Scholar

Jeynes-Cupper, K. and Catoni, M. (2023). Long-distance signalling and epigenetic changes in crop grafting. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1121704

PubMed Abstract | Crossref Full Text | Google Scholar

Karunakaran, R. and Ilango, R. V. J. (2019). Grafting influence on productivity and drought tolerance of tea clones. J. Agric. Sci. 157, 1–9. doi: 10.1017/s0021859619000480

Crossref Full Text | Google Scholar

Kawashima, C. G., Yoshimoto, N., Maruyama-Nakashita, A., Tsuchiya, Y. N., Saito, K., and Takahashi, H. (2009). Sulphur starvation induces the expression of microRNA-395 and one of its target genes but in different cell types. Plant J. 57, 313–321. doi: 10.1111/j.1365-313x.2008.03690.x

PubMed Abstract | Crossref Full Text | Google Scholar

Kyriacou, M. C., Leskovar, D. I., and Colla, G. (2017). Vegetable grafting: A toolbox for securing yield stability under multiple stress conditions. Front. Plant Sci. 8, 2255. doi: 10.3389/fpls.2017.02255

PubMed Abstract | Crossref Full Text | Google Scholar

Lan, L., Zhang, L., Cao, L., and Wang, S. (2025). WRKY1-mediated interconversion of MeSA and SA in neighbouring apple plants enhances defense against Powdery Mildew. Plant Cell Environ. 48, 3105–3120. doi: 10.1111/pce.15323

PubMed Abstract | Crossref Full Text | Google Scholar

Li, S., Wang, X., Xu, W., Liu, T., Cai, C., and Chen, L. (2021). Unidirectional movement of small RNAs from shoots to roots in interspecific heterografts. Nat. Plants 7, 23–33. doi: 10.1038/s41477-020-00829-2

PubMed Abstract | Crossref Full Text | Google Scholar

Lin, S., Chiang, S., Lin, W., Chen, J., Tseng, C., and Wu, P. (2008). Regulatory network of microRNA399 and PHO2 by systemic signaling. Plant Physiol. 147, 732–746. doi: 10.1104/pp.108.116269

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, Y., Liu, J., Tian, Y., Ye, S., Pang, D., and Chen, L. (2025). Multi-Omics analysis revealed the accumulation of flavonoids and shift of fungal community structure caused by rea grafting (Camellia sinensis L.). Plants 14, 1176. doi: 10.3390/plants14081176

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, S., Wei, J., Chen, J., Li, Y., and Lin, K. (2017). Research on fast propagation of tea-specific germplasm resources by grafting. China Agric. Bull. 33, 78–82.

Google Scholar

Melnyk, C. W., Molnar, A., Bassett, A., and Baulcombe, D. C. (2011b). Mobile 24 nt small RNAs direct transcriptional gene silencing in the root meristems of Arabidopsis thaliana. Curr. Biol. 21, 1678–1683. doi: 10.1016/j.cub.2011.08.065

PubMed Abstract | Crossref Full Text | Google Scholar

Melnyk, C. W., Molnar, A., and Baulcombe, D. C. (2011a). Intercellular and systemic movement of RNA silencing signals. EMBO J. 30, 3553–3563. doi: 10.1038/emboj.2011.274

PubMed Abstract | Crossref Full Text | Google Scholar

Millar, A. A., Lohe, A., and Wong, G. (2019). Biology and function of miR159 in plants. Plants 8, 255–272. doi: 10.3390/plants8080255

PubMed Abstract | Crossref Full Text | Google Scholar

Moreno-Piovano, G. S., Moreno, J. E., Cabello, J. V., Arce, A. L., Otegui, M. E., and Chan, R. L. (2017). A role for LAX2 in regulating xylem development and lateral-vein symmetry in the leaf. Ann. Bot. 120, 577–590. doi: 10.1093/aob/mcx091

PubMed Abstract | Crossref Full Text | Google Scholar

Pant, B. D., Buhtz, A., Kehr, J., and Scheible, W. R. (2008). MicroRNA399 is a long-distance signal for the regulation of plant phosphate homeostasis. Plant J. 53, 731–738. doi: 10.1111/j.1365-313X.2007.03363.x

PubMed Abstract | Crossref Full Text | Google Scholar

Park, S. W., Kaimoyo, E., Kumar, D., Mosher, S., and Klessig, D. F. (2007). Methyl salicylate is a critical mobile signal for plant systemic acquired resistance. Science 318, 113–116. doi: 10.1126/science.1147113

PubMed Abstract | Crossref Full Text | Google Scholar

Sharma, A. and Zheng, B. (2019). Molecular responses during plant grafting and its regulation by auxins, cytokinins, and gibberellins. Biomolecules 9, 397–414. doi: 10.3390/biom9090397

PubMed Abstract | Crossref Full Text | Google Scholar

Singh, V. K., Upadhyay, A. K., Kumar, A., and Singh, T. P. (2022). How do plants defend themselves against pathogens-Biochemical mechanisms and genetic interventions. Physiol Mol Biol Plants 28, 893–911. doi: 10.1007/s12298-022-01146-y

PubMed Abstract | Crossref Full Text | Google Scholar

Sotelo-Silveira, M., Cucinotta, M., Chauvin, A. L., Montes, R. A. C., Colombo, L., and Marsch-Martínez, N. (2013). Cytochrome P450 CYP78A9 is involved in Arabidopsis reproductive development. Plant Physiol. 162, 779–799. doi: 10.1104/pp.113.218214

PubMed Abstract | Crossref Full Text | Google Scholar

Spanò, R., Ferrara, M., Montemurro, C., Mulè, G., Gallitelli, D., and Mascia, T. (2020). Grafting alters tomato transcriptome and enhances tolerance to an airborne virus infection. Sci. Rep. 10, 2538–2551. doi: 10.1038/s41598-020-59421-5

PubMed Abstract | Crossref Full Text | Google Scholar

Sun, Q. Y., Zhang, J. X., Xia, X. J., and Wang, Y. (2017). Current status and prospects of biological control of tea anthracnose. Tea Industry Commun. 39, 39–72. doi: 10.16015/j.cnki.jteabusiness.2017.0023

Crossref Full Text | Google Scholar

Tao, S., Li, J., Gu, X., Wang, Y., Xia, Q., and Qin, B. (2012). Quantitative analysis of ATP sulfurylase and selenocysteine methyltransferase gene expression in different organs of tea plant (Camellia sinensis). Am. J. Plant Sci. 3, 44–51. doi: 10.4236/ajps.2012.31004

Crossref Full Text | Google Scholar

Wang, W. J. (2007). Preliminary report on the relationship between different rootstocks grafted with Tieguanyin and tea quality. Tea Sci. Technol. 4, 6–8. doi: 10.3969/j.issn.1007-4872.2007.04.003

Crossref Full Text | Google Scholar

Wang, Y., Feng, C., Zhai, Z., Peng, X., and Li, T. (2020b). The apple microR171i-SCARECROW-LIKE PROTEINS 26.1 module enhances drought stress tolerance by integrating ascorbic acid metabolism. Plant Physiol. 184, 1555–1571. doi: 10.1104/pp.20.00476

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, J., Li, D., Chen, N., Chen, J., Mu, C., and Yin, K. (2020a). Plant grafting relieves asymmetry of jasmonic acid response induced by wounding between scion and rootstock in tomato hypocotyl. PloS One 15, e0241317–e0241334. doi: 10.1371/journal.pone.0241317

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, P., Liu, F., Sun, Y., Liu, X., and Jin, L. (2025). Physiological and molecular insights into citrus rootstock- scion interactions: compatibility, signaling, and impact on growth, fruit quality and stress responses. Horticulturae 11, 1110–1128. doi: 10.3390/horticulturae11091110

Crossref Full Text | Google Scholar

Wang, Y. X., Teng, R. M., Wang, W. L., Wang, Y., Shen, W., and Zhuang, J. (2019). Identification of genes revealed differential expression profiles and lignin accumulation during leaf and stem development in tea plant (Camellia sinensis (L.) O. Kuntze). Protoplasma 256, 359–370. doi: 10.1007/s00709-018-1299-9

PubMed Abstract | Crossref Full Text | Google Scholar

Xia, E. H., Tong, W., Hou, Y., An, Y., Chen, L., and Wu, Q. (2020). The reference genome of tea plant and resequencing of 81 diverse accessions provide insights into genome evolution and adaptation of tea plants. Mol. Plant 13, 1013–1026. doi: 10.1016/j.molp.2020.04.010

PubMed Abstract | Crossref Full Text | Google Scholar

Xue, C., Yao, J. L., Qin, M. F., Zhang, M. Y., Allan, A. C., and Wang, D. F. (2019). PbrmiR397a regulates lignification during stone cell development in pear fruit. Plant Biotechnol. J. 17, 103–117. doi: 10.1111/pbi.12950

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, H. Q. (2010). MicroRNA171c-targeted SCL6-II, SCL6-III, and SCL6-IV genes regulate shoot branching in Arabidopsis. Mol Plant. 3, 794–806. doi: 10.1093/mp/ssq042

PubMed Abstract | Crossref Full Text | Google Scholar

Yoshida, K. and Takeda, Y. (2006). Evaluation of anthracnose resistance among tea genetic resources by wound-inoculation assay. Jpn. Agric. Res. Q. 40, 379–386. doi: 10.6090/jarq.40.379

Crossref Full Text | Google Scholar

Zhang, J., Liu, H., Yao, J., Ma, C., Yang, W., and Lei, Z. (2024). Plant-derived citronellol can significantly disrupt cell wall integrity maintenance of Colletotrichum camelliae. Pestic. Biochem. Physiol. 219, 105893. doi: 10.1016/j.pestbp.2024.106087

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, H. Y., Qi, L., and Zhang, Z. Z. (2010). Inhibitory effects of caffeine on major foliar pathogenic fungi of tea plants. J. Nanjing Agric. Univ. 33, 63–67. doi: 10.7666/d.y1597354

Crossref Full Text | Google Scholar

Zhang, Y. L. (1979). Discussion on grafting in tea hybrid breeding. Tea Commun. 6, 34–36.

Google Scholar

Zhang, Y., Fu, X., Feng, Y., Zhang, X., Bi, H., and Ai, X. (2022). Abscisic acid mediates salicylic acid induced chilling tolerance of grafted cucumber by activating H2O2 biosynthesis and accumulation. Int. J. Mol. Sci. 23, 16057. doi: 10.3390/ijms232416057

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Z., Sun, J., Zhai, W., Zhu, C., Liu, Z., and Zhang, D. (2025). Integrated transcriptomic and metabolomic analysis reveals flavonoid-mediated cold adaptation in rootstock ‘Totosga’. Scientia Hortic. 351, 114414–114427. doi: 10.1016/j.scienta.2025.114414

Crossref Full Text | Google Scholar

Zhao, Y., Ma, Y., Sheth, S., and Wang, Y. (2023). Current status and prospects of green control of major foliar fungal diseases in tea plants in China. Plant Prot. 49, 133–144 + 166. doi: 10.16688/j.zwbh.2023479

Crossref Full Text | Google Scholar

Keywords: Camellia sinensis, Colletotrichum camelliae, disease resistance, grafting, multi-omics

Citation: Li Y-X, Lin K-Q, Wang A-R, Wei J, Yan D-H and Zhao D-G (2026) Multi-omics analysis of molecular mechanisms driving the grafting- enhanced resistance of tea plants to Colletotrichum camelliae. Front. Plant Sci. 16:1750493. doi: 10.3389/fpls.2025.1750493

Received: 20 November 2025; Accepted: 22 December 2025; Revised: 17 December 2025;
Published: 14 January 2026.

Edited by:

Choong-Min Ryu, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Republic of Korea

Reviewed by:

Wenbin Jiang, Southwest University, China
Rui Han, Nanjing Agricultural University, China

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

*Correspondence: Dong-Hai Yan, ZG9uZ2hhaUAxMjYuY29t; De-Gang Zhao, ZGd6aGFvQGd6dS5lZHUuY24=

These authors have contributed equally to this work and share first authorship

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