- 1College of Horticulture, Fujian Agriculture and Forestry University, Fuzhou, China
- 2College of Tea and Food Science, Wuyi University, Wuyishan, China
- 3College of Life Science, Longyan University, Longyan, China
Cuttings are an important way of propagating tea trees (Camellia sinensis). In this study, Dahongpao mother tree (MD) and cutting Dahongpao (PD) were used as research objects and their rhizosphere soil were collected and performed metabolomics analysis. At the same time, soil nutrient content, microbial physiological indexes, and microbial carbon source utilization were determined, which in turn obtained the effect of cuttings on metabolites, microorganisms, and nutrient cycling in rhizosphere soil of tea trees. The results showed that available nitrogen, available phosphorus and available potassium in the rhizosphere soil of MD were significantly higher (p < 0.05) than in PD. Secondly, microbial biomass carbon, microbial biomass nitrogen, microbial respiration, bacterial number, fungal number, and actinomycete number were also significantly higher in rhizosphere soil of MD than in PD. There were six groups of rhizosphere soil characteristic metabolites that differentiated MD from PD, of which the content of acid, amine, phenol, heterocyclic compound, alcohol and lipid was significantly higher in MD compared to PD, while carbohydrate content was significantly less in MD. There were five groups of rhizosphere soil microorganisms that differentiated MD from PD, in which microorganisms with carboxylic acid, amines, fatty acid and phenolic acid as carbon sources were significantly larger in MD than in PD, whereas microorganisms with carbohydrates as carbon sources were significantly smaller in MD than in PD. It can be seen that the number and content of rhizosphere soil characteristic metabolites were higher in MD than in PD. This enhanced the number of microorganisms with different carbon source utilization rates, increased microbial diversity and abundance, promoted nutrient transformation, increased the content of available nutrients, which in turn facilitated the growth of tea trees. This study provides an important reference for the use of metabolites to regulate soil microbial colonization, improve soil nutrient transformation, and maintain healthy growth of tea trees.
1 Introduction
There are numerous methods of artificial asexual propagation of plants, and cuttings are the predominant method of seedling production due to their rapid and efficient characteristics (Megersa, 2017). Cuttings are propagated asexually from mature and good parent plants, whose genomes are theoretically identical to those of the parents and can efficiently inherit the basic characteristics of the parents (Nybom and Lācis, 2021). However, plants from cuttings are susceptible to disturbances from external environmental factors during the growing process, leading to changes in the intensity of gene expression within their different tissues, which in turn affects their metabolic functions (Gil et al., 2020). Second, from the point of view of growth and resistance, even if parents and cuttings are planted in the same environment, there are significant differences in both age and growth, which in turn lead to significant differences in growth rate and resistance strength (Hewitt, 2020; Xiao et al., 2023). It can be seen that although plant cuttings are asexually propagated, there are still some differences between them and their parents, which in turn may affect their growth.
Soil is the medium in which plants are grown, and the rhizosphere is the area in which plants are in close contact with the soil. Changes in the community structure and function of rhizosphere soil microorganisms directly affect the rhizosphere micro-ecosystem, which in turn affects plant growth (Zhao et al., 2021). Singh (2021) found that rhizosphere soil microbial diversity remained significantly different when parent plants and their asexually propagated cuttings were planted in the same plot and managed in the same way. Wang et al. (2024a) subjected Gastrodia to multi-generation asexual propagation and analyzed the effect of multi-generation asexually propagated cuttings on soil microorganisms and found that the microbial diversity of the rhizosphere soil of Gastrodia declined significantly with the increase in the number of generations of propagation, and the microbial community structure underwent a significant change, which was manifested by the gradual decrease of probiotic bacteria and the continuous increase of pathogenic bacteria. Lin et al. (2024) investigated microbial changes in the rhizosphere soil after asexually propagated sugarcane planting and found significant changes in rhizosphere soil microbial functions, particularly a significant decrease in microbes related to nutrient transformation, which in turn led to a decrease in soil enzyme activities and a reduction in soil nutrient biotransformation. Adomako et al. (2022) explored the rhizosphere soil nutrient transformation capacity of asexually propagated Solidago canadensis and found that the nitrogen-phosphorus ratio of the rhizosphere soil of the asexually propagated cuttings was significantly altered, soil nutrient ratios were imbalanced, and the productivity of Solidago canadensis was significantly reduced, as compared to the parent. It can be seen that asexually propagated cuttings, although inheriting the characteristics of the parent, have undergone changes in their rhizosphere soil microbial diversity and function during planting, which in turn may have affected the biotransformation of soil nutrients and plant growth.
Tea tree is an important economic plant, and Dahongpao mother tree (Camellia sinensis) is an icon of the tea industry in Wuyishan City, Fujian, China, and has been listed by the local government as a key protected object (Chen, 2015). In the mid-1980s, a batch of cuttings were successfully produced by asexual propagation of cuttings for the first time using Dahongpao mother tree as the mother tree (Hou and Wu, 2019). From 2006 to now, in order to effectively protect Dahongpao mother tree, the local government has banned its harvesting for tea production (Ng et al., 2018). As a result, most Dahongpao tea sold on the market today is processed from tea leaves from cuttings of Dahongpao mother tree that were raised into tea trees. Ye et al. (2024) analyzed the quality of Dahongpao mother tree and cutting Dahongpao and found that compared to the parent plant, cuttings had significantly lower aromas such as floral, fruity, green and woody, and significantly lower taste characteristics such as fresh and brisk taste and mellowness. Jia et al. (2024b) analyzed the rhizosphere soil microbial diversity of Dahongpao mother tree and their cuttings and found that the rhizosphere soil microbial diversity of cuttings was significantly reduced, and function was reduced, which in turn affected nutrient uptake by tea trees, and tea leaf quality was reduced. It has been reported that the accumulation of root secretions in rhizosphere soil during plant cultivation is highly susceptible to altering the structure and function of soil microbial communities, which in turn affects soil nutrient transformations and influences plant growth and quality (Shen and Lin, 2021; Tan et al., 2022). The type and content of rhizosphere soil metabolites significantly influence the number and function of soil microorganisms that are addicted to different types of metabolites (Withers et al., 2020). Soil metabolites regulate changes in microbial function, thus affecting soil nutrient transformation and plant growth (Brown et al., 2021; Sun et al., 2022; Wu et al., 2022). It is hypothesized that the changes in growth and quality of cutting Dahongpao may be related to changes in the type or quantity of their rhizosphere soil metabolites. It is of great significance to deeply reveal the metabolite differences between Dahongpao mother trees and cutting Dahongpao, and to search for characteristic metabolites and their effects on soil microbes and nutrient cycling, in order to utilize the metabolites to regulate the reproduction of soil microbes, to regulate the transformation of soil nutrients, and to maintain the healthy growth of tea trees.
Accordingly, in this study, Dahongpao mother tree (MD) and the first asexually propagated cuttings of Dahongpao (PD) in the 1980s were used as objects, and their rhizosphere soils were collected to determine available nutrient content, microbial number, and physiological indexes. Carbon source utilization by rhizosphere soil microorganisms of MD and PD was determined using BIOLOG ECO microplate method. At the same time, metabolomics technology was used to determine rhizosphere soil metabolite compositions and their contents of MD and PD, and to screen for characteristic metabolites distinguishing MD and PD. On this basis, the interactions between soil characteristic metabolites, microorganisms and soil nutrients were further analyzed, with a view to laying a theoretical foundation for metabolites to regulate soil microbial reproduction, improve soil nutrient transformation, and maintain the healthy growth of tea trees.
2 Materials and methods
2.1 Collection of experimental samples
The sampling site of this study was located in Jiulongke Scenic Area (117°57′19.098″ E, 27°40′17.8212″ N), Wuyishan, Fujian Province, China. In May 2023, rhizosphere soils were collected from MD and PD with three independent replicates of each sample. Among them, MD is about 390 years old, while PD is about 40 years old, and both tea trees are planted in the same plot at a distance of about 20 m. The soil is red gravelly rock. The average annual temperature of the planting site is about 13°C, the precipitation is above 2000 mm and the relative humidity is up to 85%. The sampling method of tea tree rhizosphere soil was as follows: 2 plants of each tea tree were randomly selected, and the surface soil was shoveled layer by layer until the root system appeared (about 20 cm), and the soil attached to the root system was collected, and thoroughly mixed, which was rhizosphere soil. The collected rhizosphere soil was immediately to the ice box, part of it was air-dried and used for the determination of available nitrogen, phosphorus and potassium content, and part of it was used for the determination of microbial population, microbial physiological indexes, microbial carbon source utilization and soil metabolites, with three independent replicates for each sample.
2.2 Determination of soil available nitrogen, phosphorus and potassium content
The available nitrogen, phosphorus and potassium content of tea tree rhizosphere soil was determined by Wang et al. (2024b) with three independent replicates for each sample. Briefly, soil available nitrogen was determined by leaching with NaOH solution, and leachate was titrated with hydrochloric acid and then converted to content. Soil available phosphorus was extracted by NaHCO3, the extract was added molybdenum-antimony resistance for color development, and then determined by colorimetric method, finally converted to content. Soil available potassium was extracted using ammonium acetate, and the extract was measured directly by flame photometer and then converted to content.
2.3 Determination of physiological indexes of soil microorganisms
Microbial biomass carbon, biomass nitrogen and respiration of tea tree rhizosphere soil were determined with reference to Schnecker et al. (2023) with three independent replicates of each sample. Briefly, microbial biomass carbon and nitrogen were determined by chloroform fumigation extraction, i.e., soil samples were fumigated in chloroform for 24 h, extracted with 1 M KCl, and then determined by a TOC/TN analyzer (TOC-l CPH/CPN, Shimadzu, Kyoto, Japan). Where microbial biomass carbon is calculated as (fumigated organic carbon - unfumigated organic carbon)/0.38 and microbial biomass nitrogen is calculated as (fumigated total nitrogen - unfumigated total nitrogen)/0.54. The intensity of soil microbial respiration was measured by the alkali absorption method in mg CO2/kg·h, i.e., it was assessed based on the amount of CO2 released per kilogram of soil per hour.
2.4 Determination of soil bacteria, fungi and actinomycetes number
Quantitative analysis of bacteria, fungi and actinomycetes in tea tree rhizosphere soil was performed by qRT-PCR with three independent replicates per sample, as described in Ye et al. (2023). Briefly, 0.5 g of fresh soil was taken and soil DNA was extracted using the Bio-Fast Soil Genomic DNA (BioFlux, Hangzhou, China) extraction kit, and DNA was purified using the gel recovery kit of TianGen Biotech Co., Ltd.The primers used for bacterial quantification were F27 (5′-AGAGTTTGATCMTGCCTCAG-3′) and R1492 (5′-TACHHYTACCTTGTTACGACTT-3′). The bacterial PCR program was set to 95°C pre-denaturation for 4 min, 94°C for 1 min, 55°C for 1 min, 72°C for 1 min, and 35 cycles. The primers used for fungal quantification were 5.8S (5′-CGCTGCGTTCTTCATCG-3′) and ITSIF (5′-CTTGGTCATTTAGAGGAAGTAA-3′). The fungal PCR program was set to 95°C pre-denaturation for 15 s, 94°C for 35 s, 53°C for 30 s, 72°C for 30 s with 35 cycles. The primers used for actinomycetes quantification were act920f (5′-TACGGCCGCAAGGCTA-3′) and act1200r (5′-TCRTCCCCACCTTCCTCCG-3′). The actinomycetes PCR program was set to pre-denaturation at 95°C for 5 min, 95°C for 15 s, 65°C for 30 s, 72°C for 15 s, and 30 cycles.
2.5 Extraction and determination of soil metabolites
Extraction and derivatization of rhizosphere soil metabolites of tea trees were performed (Fu et al., 2022). Fresh soil samples were vacuum freeze-dried and ground to powder with three independent replicates for each sample. 0.5 g of the sample was weighed, and 1 mL of extraction solution (methanol:isopropanol:water in the ratio of 3:3:2, v/v) was added, and after shaking for 3 min at room temperature, the sample was placed in an ice bath and ultrasonicated for 20 min, and centrifuged for 3 min at 12,000 r/min at 4°C, and the supernatant was collected. The supernatant was added with 0.02 mL of internal standard (10 μg/mL), mixed well, blown dry under nitrogen, added with 0.1 mL of methoxamine pyridine solution (0.015 g/mL), oximilized at 37°C for 2 h, then added with 0.1 mL of BSTFA (containing 1% TMCS), and the reaction was carried out at 37°C for 30 min to obtain the derivatization solution. The derivatization solution was diluted to 1 mL and passed through a 0.22 μm organic filter membrane for GC-MS testing.
The GC-MS equipment used for the determination of soil metabolites was an Agilent 8890 + 5977B gas chromatography-mass spectrometry instrument (Agilent, Palo Alto, California, USA), and the chromatographic column used was a DB-5MS (30 m × 0.25 mm × 0.25 μm, J&W Scientific, USA). The parameters of GC-MS were set as follows: carrier gas was helium, injection volume was 1 μL, front inlet mode was 5:1, flow rate was 1.2 mL/min; Oven temperature ramp was held at 40°C for 1 min, raised to 100°C at a rate of 20°C/min, raised to 300°C at a rate of 15°C/min, and held at 300°C for 5 min; Transfer line temperature was 280°C, ion source temperature was 230°C, quad temperature was 150°C, and electron energy was 70 eV. Qualitative and quantitative methods for soil metabolites were performed through selecting 2 ~ 3 qualitative ions and 1 quantitative ion for each compound during the determination. and then compared with the NIST20 mass spectrometry database. A compound is qualitative when the selected ions, net of background, all appear in the mass spectrum and the retention time is consistent with the reference value. On this basis, the compound can be quantified by integrating and correcting ions according to the chosen quantification ion (Yuan et al., 2022).
2.6 Determination of microbial carbon source utilization
The BIOLOG ECO microplate has 96 wells containing 31 carbon sources, 3 wells per carbon source, i.e. 3 replicates, and 3 blank control wells. The 31 carbon sources can be categorized into six groups, namely carbohydrate, carboxylic acid, phenolic acid, fatty acid, amines and amino acid, respectively (Li et al., 2024). In this study, the utilization of different carbon sources by rhizosphere soil microorganisms of MD and PD was determined using BIOLOG ECO microplate method with reference to Wang et al. (2023). Briefly, 10 g of fresh rhizosphere soil was taken in a conical flask, 90 mL of sterile saline was added, sealed, shaken, and then placed in a shaker at 120 r/min for 10 min. 5 mL of supernatant was taken and diluted 10-fold with sterile water and left to stand, then 5 mL of supernatant was taken again and diluted 10-fold with sterile water to obtain a dilution of 1:1000 for BIOLOG ECO microplate assay. For the BIOLOG ECO microplate assay, 150μL of dilution was added to each well, and a blank well with an equal volume of sterile water was added as a control. The BIOLOG ECO microplate was incubated in a constant temperature incubator at 28°C, protected from light for 7 days, and then the absorbance was measured at 590 nm. The utilization of different carbon sources by rhizosphere soil microorganisms was expressed as Average well color development (AWCD) per pore. AWCD = [∑(C-R)]/31, where C is the absorbance measured after 7 days of incubation in each well and R is the absorbance of the control well.
2.7 Statistical analysis
The data obtained in this study were first performed using Excel 2020 for a preliminary statistical analysis of the raw data. Data variances were analyzed using variance analysis (ANOVA) and paired Student ‘s t-tests. The data were plotted using Rstudio software (version 4.2.3), with box and violin plots produced using the R package gghalves 0.1.4, principal component plots produced using the R package ggbiplot 0.55, and bubble heat maps produced using the R package ggplot2 3.5.1. The R package used for orthogonal partial least squares discrimination analysis (OPLS-DA) model construction for MD and PD is ropls and mixOmics, the R package used for bubble feature maps is ggplot2 3.4.0, the R package used for TOPSIS is dplyr 1.1.4, and vegan 2.6.4 for redundancy analysis. The R package used for constructing partial least squares structural squation modeling (PLS-SEM) equations for the different indexes was plspm 0.4.9, and the R package used for correlation-interaction network maps was linkET 0.0.7.1.
3 Results and discussion
3.1 Analysis of soil available nutrient content, microbial number and physiological indexes
Plant growth requires nutrient uptake from the soil. Therefore, high or low soil nutrient content, especially available nutrient content, directly or indirectly affects nutrient uptake and accumulation by the plant root system, which influences plant growth (Soares et al., 2019). In this study, available nutrient content of rhizosphere soil of Dahongpao mother tree (MD) and cutting Dahongpao (PD) and found (Figure 1A) that available nitrogen, available phosphorus, and available potassium contents of rhizosphere soil of MD were significantly higher than those of PD (p < 0.05). Soil available nutrient content is closely related to the biotransformation capacity of soil to nutrients, and the higher the transformation capacity, the more favorable the release of nutrients, which in turn increases the available nutrient content (Macik et al., 2020). The transformation of soil nutrients requires the participation of microorganisms, which can change the soil environment, alter the biotransformation capacity of nutrients, and affect plant growth (Zhong et al., 2020; Bai et al., 2021). In this study, further analysis of rhizosphere soil microbial physiological indexes of MD and PD revealed (Figure 1B) that microbial biomass carbon content, microbial biomass nitrogen content, and microbial respiration of MD rhizosphere soil were 178.66 mg/kg, 69.13 mg/kg, and 18.16 mg CO2/kg·h, respectively, while those of PD were 142.10 mg/kg, 46.76 mg/kg, and 13.49 mg CO2/kg·h, respectively, and MD was significantly larger than PD (p < 0.05). Secondly, the analysis of bacterial, fungal and actinomycete numbers showed (Figure 1C) that those of MD were all significantly higher (p < 0.05) than PD, as evidenced by the fact that the bacterial, fungal and actinomycete numbers in the rhizosphere soil of MD were 16.29×109 cell/g·soil, 5.12×109 cell/g·soil and 6.40×109 cell/g·soil, respectively, while those of PD were 10.32×109 cell/g·soil, 3.14×109 cell/g·soil and 3.21×109 cell/g·soil, respectively. The increase of soil microbial biomass carbon and nitrogen facilitates the mineralization of soil nutrients and improves the biotransformation of soil nutrients (Ma et al., 2020), while microbial respiration characterizes the number of microorganisms in the soil, and the stronger microbial respiration is, the higher the number of microorganisms is, and the more diverse the microorganisms are (Chen et al., 2021). It can be seen that compared with PD, the number of soil microorganisms and their respiration intensity in the rhizosphere of MD were significantly enhanced, which increased the microbial biomass carbon and nitrogen of the soil, enhanced the mineralization capacity of soil nutrients, and promoted soil nutrient biotransformation, which was more conducive to the improvement of soil available nutrient content, and then promoted the growth of tea trees.

Figure 1. Analysis of available nutrient content and microbial number of rhizosphere soil of tea trees. MD, Dahongpao mother tree; PD, Cutting Dahongpao; (A) Analysis of available nutrient content; (B) Analysis of microbial physiological indexes; (C) Analysis of microbial number.
3.2 Soil metabolite analysis
Soil metabolomics can effectively analyze low molecular weight compounds in plant rhizosphere soils and thus assess soil texture (Withers et al., 2020). And, soil metabolomics characterizes the metabolic state of soil biomes and can be used to assess soil function (Wu et al., 2022). Therefore, in this study, soil metabolomics technique was used to analyze the rhizosphere soil metabolites of MD and PD, and the results showed (Figure 2A) that the rhizosphere soil metabolite content of MD was significantly higher (p < 0.05) than that of PD. PCA analysis showed (Figure 2B) that the soil metabolites of MD and PD were significantly different, and the two principal components could effectively differentiate MD from PD with a total contribution of 91.76%. Further analysis of rhizosphere soil metabolites in MD and PD revealed (Figure 2C) that a total of 187 metabolites were detected, which can be categorized into 18 groups (including others), of which, 16 groups of metabolites were significantly higher in MD than in PD, namelyterpenes, phenol, nitrogen compounds, lipid, ketone, hydrocarbons, heterocyclic compound, ester, aromatics, amino acid, amine, aldehyde, alcohol, acid and others, whereas 2 groups of metabolite were significantly less in MD than in PD, namely organic acid and carbohydrate. It can be seen that there was a significant difference in rhizosphere soil metabolites between MD and PD, with MD having significantly higher soil metabolite content than PD. Soil metabolite content and types play important roles in the regulation of soil microbial community structure and function (Sikder and Vestergard, 2020). The abundance and diversity of soil metabolites contribute to the diversity of soil microbial communities, which in turn enriches soil microbial functions, enhances soil nutrient biotransformation, and promotes plant growth (Ni et al., 2021; Ortiz and Sansinenea, 2022). It can be seen that the total amount of metabolites and the content of different metabolic types were higher in rhizosphere soil of MD compared with PD, which is more conducive to promoting the proliferation of different types of microorganisms, and more conducive to increasing the richness and diversity of soil microbial community, which in turn promotes the transformation of soil nutrients and influences the growth of tea trees.

Figure 2. Analysis of metabolite content of rhizosphere soil of tea tree. MD, Dahongpao mother tree; PD, Cutting Dahongpao; (A) Total amount analysis of MD and PD rhizosphere soil metabolites; (B) PCA analysis of MD and PD rhizosphere soil metabolites; (C) Content analysis of MD and PD rhizosphere soil metabolites after classification.
3.3 Screening and content analysis of soil characteristic metabolites
The OPLS-DA model has an important role in screening key metabolites among different samples, which can be used to screen key differential metabolites by the metabolite’s variable importance projection value (VIP) in distinguishing between different samples (Jia et al., 2024a). However, the model needs to be evaluated after construction, and only the model with significant levels of fit and predictability can be used for screening and analysis (Rivera-Pérez et al., 2022; Li et al., 2023b). Accordingly, based on the above analysis, this study constructed an OPLS-DA model of MD and PD with detected soil metabolites and their contents, and screened for key soil metabolites that differentiated MD from PD. The result showed (Figure 3A) that the model constructed by MD and PD was tested with a goodness-of-fit R2Y value of 0.999 and a predictability Q2 value of 0.988, both at significant levels (p < 0.05). The OPLS-DA score plot showed (Figure 3A) that MD and PD were effectively differentiated, with a difference of 85.30% between groups and only 4.91% within groups. The VIP values of different metabolites in distinguishing MD from PD were obtained by S-Plot plots, and a total of 138 key metabolites with VIP greater than 1 were obtained (Figure 3A). The 138 metabolites were further analyzed using bubble feature maps, and 36 characteristic metabolites were screened (Figure 3B). Characteristic metabolite categorization analysis showed (Figure 3C) that the 36 characteristic metabolites could be classified into 12 groups (including others), of which 11 groups of metabolites were significantly greater in MD than PD, namely terpenes, phenol, lipid, heterocyclic compound, ester, aromatics, amino acid, amine, alcohol, acid and others, while only carbohydrate content in MD was significantly lower than PD. TOPSIS was used to analyze the contribution of each of the 11 groups of metabolites in differentiating MD from PD, and the result showed (Figure 3D) that only 7 groups of metabolites contributed more than 10% in differentiating MD from PD, namely acid (96.62%), carbohydrate (53.98%), amine (24.35%), phenol (22.36%), heterocyclic compound (20.72%), alcohol (17.02%) and lipid (16.41%). Tea tree has been reported to be an acidophilic plant with low soil pH, which is suitable for the colonization of acidophilic microorganisms (Yang et al., 2024). Both phenol and alcohol can be converted to acids upon oxidation (De Nobili et al., 2020). Moderate amounts of Amine stimulate nitrogen conversion by soil microorganisms and further convert amine into plant-available and available nitrogen (Paśmionka et al., 2021). Furthermore, when the content of heterocyclic compound in soil is high, it can promote the propagation of bacteria and fungi in soil, increase the efficiency of carbon utilization by microorganisms, enhance the conversion capacity of soil carbon, and promote the growth of plants (Li et al., 2023a). Lipid serves as an important raw material for microbial colonization, which reduces the depletion of available phosphorus in the soil and promotes phosphorus uptake by plants (Warren, 2020). And carbohydrate is an important raw material for microbial reproduction and a major carbon source (Xie et al., 2019). It can be seen that there were significant differences between MD and PD in rhizosphere soil metabolite content, especially characteristic metabolites. The carbohydrate content was significantly higher in the rhizosphere soil of PD compared to MD, which could lead to colonization of rhizosphere soil microorganisms that primarily used carbohydrate as a carbon source. The content of acid, amine, phenol, heterocyclic compound, alcohol and lipid in the rhizosphere soil of MD were significantly higher than that of PD, which could promote the propagation of microorganisms that used different types of carbon sources as raw materials, and was more conducive to increasing the diversity and abundance of the soil microbial community, which was in turn conducive to promoting the transformation of different types of nutrients in the soil.

Figure 3. Screening of characteristic metabolites of tea tree rhizosphere soil. MD, Dahongpao mother tree; PD, Cutting Dahongpao; (A) Construction of OPLS-DA model for MD and PD and screening for key differential metabolites; (B) Screening for characteristic metabolites distinguishing MD from PD by bubble feature map; (C) Content analysis after categorization of characteristic metabolites; (D) TOPSIS analysis of the contribution of different groups of metabolites in distinguishing MD from PD.
3.4 Microbial carbon source utilization analysis
During growth, plants can maintain their own growth by releasing root secretions that regulate the type and amount of rhizosphere soil metabolites in order to adapt to changes in the environment, which in turn affects microbial colonization and alters the structure of the soil microbial community and soil nutrient content (Barra and Terenzi, 2021; Bi et al., 2022). In the present study, significant changes in metabolites were found in rhizosphere soil of MD and PD, which were hypothesized to possibly affect the colonization of the corresponding microorganisms. Accordingly, in this study, the BIOLOG microplate method was used to analyze the number of rhizosphere soil microorganisms utilizing different types of carbon sources in MD and PD, and the results showed (Figure 4A) that the overall utilization of carbon sources by rhizosphere soil microorganisms was significantly higher in MD than in PD, and for the utilization of different types of carbon sources, the number of microorganisms using carboxylic acid, phenolic acid, fatty acid, amino acid and amines as carbon sources was significantly higher in MD compared to PD, whereas the number of microorganisms using carbohydrate as carbon source was significantly lower in MD. PCA analysis with different carbon source utilization rates of microorganisms found (Figure 4B) that principal component 1 and principal component 2 could effectively differentiate MD from PD, with 97.70% contribution from principal component 1 and 1.11% from principal component 2, giving an overall contribution of 98.81%. And it was found that microorganisms with carboxylic acid, phenolic acid, fatty acid, amino acid and amines as carbon sources were significantly associated with MD, whereas microorganisms with carbohydrates as carbon sources were significantly associated with PD. TOPSIS was further used to analyze the contribution of different microorganisms in distinguishing MD from PD, and the results showed (Figure 4C) that only five types of microorganisms that contributed more than 1% to distinguishing MD from PD were microorganisms that used carbohydrate, carboxylic acid, amines, fatty acid and phenolic acid as carbon source, respectively. Phenol can be oxidized to phenolic acid (Rashmi and Negi, 2020), alcohol can be oxidized to carboxylic acid (Pradhan et al., 2020), and lipid can be converted to fatty acid (Dyall et al., 2022). And heterocyclic compound is conducive to promoting soil microbial colonization and enhancing soil nutrient biotransformation (Li et al., 2023a). It can be seen that the content of acid, amine, phenol, heterocyclic compound, alcohol and lipid in soil metabolites was significantly higher in MD than in PD, which in turn was more favorable to promote microbial colonization with carboxylic acid, amines, fatty acid and phenolic acid as carbon sources, whereas PD was more favorable to promote microbial colonization with carbohydrates as carbon sources. The number of microorganisms utilizing different carbon sources in rhizosphere soil of MD was significantly higher than that of PD, and the microbial diversity and abundance of soil was higher, which in turn was more conducive to promoting soil nutrient transformation and tea tree growth.

Figure 4. Analysis of microbial carbon source utilization in rhizosphere soil of tea trees. MD, Dahongpao mother tree; PD, Cutting Dahongpao; (A) Quantitative analysis of microorganisms utilizing different carbon sources; (B) PCA analysis of rhizosphere soils microorganisms utilizing different carbon sources in MD and PD; (C) TOPSIS analysis of the contribution rate of microorganisms utilizing different carbon sources in distinguishing MD from PD.
3.5 Interaction analysis
In this study, characteristic metabolites in rhizosphere soil of tea trees, microorganisms utilizing different carbon sources were further analyzed for interaction effects with soil available nutrients and microbial physiological indexes. Redundancy analysis showed (Figure 5A) that among soil characteristic metabolites, acid, amine, phenol, heterocyclic compound, alcohol and lipid were significantly associated with MD, and soil microorganisms with carboxylic acid, phenolic acid, fatty acid, amino acid and amines as carbon sources were significantly associated with MD. In contrast, both characteristic metabolites and microorganisms utilizing different carbon sources were correlated with soil available nutrient content and microbial physiological indexes. Correlation interaction network analysis showed (Figure 5B) that soil characteristic metabolites, microorganisms utilizing different carbon sources, soil available nutrients, and microbial physiological indexes associated with MD were significantly and positively correlated. The PLS-SEM equations of different indexes were further constructed, and the results showed (Figure 5C) that rhizosphere soil characteristic metabolites positively regulated microorganisms using different carbon sources in the soil (1.00**), positively regulated soil microbial numbers and physiological indexes (0.98**), and thus positively regulated soil available nutrient content (0.97**). It can be seen that the higher number and content of characteristic metabolites in rhizosphere soil of MD compared to PD was more conducive to increasing the microorganisms using different carbon sources in the soil, which in turn increased the diversity and abundance of soil microorganisms, promoted soil nutrient transformation, improved soil available nutrient content, and promoted tea tree growth.

Figure 5. Interaction analysis of soil available nutrients, microorganisms and metabolites. MD, Dahongpao mother tree; PD, Cutting Dahongpao; (A) redundancy analysis; (B) correlation-interaction network analysis; (C) PLS-SEM equation analysis.
4 Conclusion
In this study, Dahongpao mother tree (MD) and cutting Dahongpao (PD) were analyzed for rhizosphere soil available nutrient content, microbial physiological indexes, microbial carbon source utilization and soil metabolites. It was found (Figure 6) that MD had significantly higher levels of acid, amine, phenol, heterocyclic compound, alcohol and lipid in rhizosphere soil metabolites than PD, which in turn was more conducive to promoting rhizosphere microbial colonization using carboxylic acid, amines, fatty acid and phenolic acid as carbon sources, whereas PD was more conducive to promoting microbial colonization using carbohydrates as carbon sources. Interaction analysis showed that rhizosphere soil characteristic metabolites positively regulated soil microorganisms using different carbon sources, positively regulated soil microbial numbers and physiological indexes, and thus positively regulated soil available nutrient content. In conclusion, the higher content of characteristic metabolites in the rhizosphere soil of MD compared with PD was more conducive to stimulating the number of microorganisms utilizing different carbon sources, which in turn increased the diversity and abundance of soil microorganisms, and was more conducive to facilitating the transformation of soil nutrients, increasing the content of soil available nutrients, and promoting the growth of tea trees. This study analyzes the effects of cuttings on the rhizosphere soil microecosystem of tea trees from the perspectives of soil metabolites and microorganisms, which is of great guiding significance for the use of metabolites to regulate the propagation of soil microorganisms, improve soil nutrient transformation, and maintain the healthy growth of tea trees.

Figure 6. Mechanism analysis of tea rhizosphere soil metabolites regulating soil microbial community affecting soil nutrient cycling.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Author contributions
WC: Conceptualization, Data curation, Formal analysis, Writing – original draft. SZ: Conceptualization, Data curation, Formal analysis, Writing – original draft. YHW: Formal analysis, Investigation, Visualization, Writing – original draft. LH: Formal analysis, Investigation, Visualization, Writing – original draft. MQ: Formal analysis, Investigation, Visualization, Writing – original draft. YLW: Formal analysis, Investigation, Visualization, Writing – original draft. YL: Formal analysis, Investigation, Visualization, Writing – original draft. QZ: Data curation, Methodology, Writing – review & editing. TW: Data curation, Methodology, Writing – review & editing. XJ: Data curation, Methodology, Writing – review & editing. HW: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing. JY: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by Natural Science Foundation of Fujian Province (2024N0009, 2024J01861, 2024J01866, 2021J011137); Construction of first-class undergraduate specialty (tea science) in Fujian Province (SJZY2019004); Nanping City Science and Technology Plan Project (N2024Z010, NP2021KTS06, NP2021KTS05, N2021Z012); Innovation and Entrepreneurship Training Program for College Students (202411312002, S202210397033, 202310397003, 202310397013); Faculty and students co-creation team of Wuyi University (2021-SSTD-01, 2021-SSTD-05).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
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.1508622/full#supplementary-material
Supplementary Table 1 | Metabolites in rhizosphere soil of MD and PD.
References
Adomako, M. O., Xue, W., Du, D. L., and Yu, F. H. (2022). Soil microbe-mediated N:P stoichiometric effects on Solidago canadensis performance depend on nutrient levels. Microb Ecol. 83, 960–970. doi: 10.1007/s00248-021-01814-8
Bai, X., Dippold, M. A., An, S., Wang, B., Zhang, H., and Loeppmann, S. (2021). Extracellular enzyme activity and stoichiometry: The effect of soil microbial element limitation during leaf litter decomposition. Ecol. Indic. 121, 107200. doi: 10.1016/j.ecolind.2020.107200
Barra, De NobiliA. and Terenzi, V. (2021). Rhizosphere microbial communities and heavy metals. Microorganisms. 9, 1462. doi: 10.3390/microorganisms9071462
Bi, B., Yuan, Y., Zhang, H., Wu, Z., Wang, Y., and Han, F. (2022). Rhizosphere soil metabolites mediated microbial community changes of Pinus sylvestris var. mongolica across stand ages in the Mu Us Desert. Appl. Soil Ecol. 169, 104222. doi: 10.1016/j.apsoil.2021.104222
Brown, R. W., Chadwick, D. R., Zang, H., and Jones, D. L. (2021). Use of metabolomics to quantify changes in soil microbial function in response to fertiliser nitrogen supply and extreme drought. Soil Biol. Biochem. 160, 108351. doi: 10.1016/j.soilbio.2021.108351
Chen, D. H. (2015). Study on Wuyi rock tea (Dahongpao) (Beijing, China: China Science and Technology Press).
Chen, L. F., He, Z. B., Wu, X. R., Du, J., Zhu, X., Lin, P. F., et al. (2021). Linkages between soil respiration and microbial communities following afforestation of alpine grasslands in the northeastern Tibetan Plateau. Appl. Soil Ecol. 161, 103882. doi: 10.1016/j.apsoil.2021.103882
De Nobili, M., Bravo, C., and Chen, Y. (2020). The spontaneous secondary synthesis of soil organic matter components: A critical examination of the soil continuum model theory. Appl. Soil Ecol. 154, 103655. doi: 10.1016/j.apsoil.2020.103655
Dyall, S. C., Balas, L., Bazan, N. G., Brenna, J. T., Chiang, N., da Costa Souza, F., et al. (2022). Polyunsaturated fatty acids and fatty acid-derived lipid mediators: Recent advances in the understanding of their biosynthesis, structures, and functions. Prog. Lipid Res. 86, 101165. doi: 10.1016/j.plipres.2022.101165
Fu, Q., Lai, J. L., Ji, X. H., Luo, Z. X., Wu, G., and Luo, X. G. (2022). Alterations of the rhizosphere soil microbial community composition and metabolite profiles of Zea mays by polyethylene-particles of different molecular weights. J. Hazard. Mater. 423, 127062. doi: 10.1016/j.jhazmat.2021.127062
Gil, C. S., Jung, H. Y., Lee, C., and Eom, S. H. (2020). Blue light and NAA treatment significantly improve rooting on single leaf-bud cutting of Chrysanthemum via upregulated rooting-related genes. Sci. Hortic. 274, 109650. doi: 10.1016/j.scienta.2020.109650
Hewitt, A. (2020). Genetic and environmental factors in the trade-off between sexual and asexual reproduction of a rare clonal angiosperm. Austral. Ecol. 45, 187–194. doi: 10.1111/aec.12846
Hou, D. W. and Wu, Y. L. (2019). Research on the brand construction of Dahongpao tea industry. Issues For. Econ. 39, 413–420. doi: 10.16832/j.cnki.1005-9709.2019.04.011
Jia, M., Cheng, P., Wang, Y., Pang, X., Li, M., Hong, L., et al. (2024a). Effects of aviation mutagenesis on soil chemical indexes, enzyme activities, and metabolites of Dahongpao (Camellia sinensis) tea trees. Plants. 13, 1291. doi: 10.3390/plants13101291
Jia, X., Li, M., Zhang, Q., Jia, M., Hong, L., Zhang, S., et al. (2024b). Analysis of rhizosphere soil microbial diversity and its functions between Dahongpao mother tree and cutting Dahongpao. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1444436
Li, J., Chen, L., Zhang, J., Zhang, C., Ma, D., Zhou, G., et al. (2023a). Organic amendments with high proportion of heterocyclic compounds promote soil microbiome shift and microbial use efficiency of straw-C. Front. Microbiol. 14. doi: 10.3389/fmicb.2023.1087709
Li, W., Zhuang, J., Xu, P., Xu, X., and BiBi, Z. (2024). Soil microbial carbon source utilization under different land use patterns along the Yellow River, China. Appl. Ecol. Env. Res. 22, 1339–1353. doi: 10.15666/aeer/2202_13391353
Lin, X., Yang, D., Zhu, Y., Qin, Y., Liang, T., Yang, S., et al. (2024). Changes in root metabolites and soil microbial community structures in rhizospheres of sugarcanes under different propagation methods. Microb. Biotechnol. 17, e14372. doi: 10.1111/1751-7915.14372
Li, Z., Sun, X., Xu, T., Dai, W., Yan, Q., Li, P., et al. (2023b). Insight into the dynamic variation and retention of major aroma volatile compounds during the milling of Suxiang japonica rice. Food Chem. 405, 134468. doi: 10.1016/j.foodchem.2022.134468
Ma, Q., Wen, Y., Wang, D., Sun, X., Hill, P. W., Macdonald, A., et al. (2020). Farmyard manure applications stimulate soil carbon and nitrogen cycling by boosting microbial biomass rather than changing its community composition. Soil Biol. Biochem. 144, 107760. doi: 10.1016/j.soilbio.2020.107760
Macik, M., Gryta, A., and Frąc, M. (2020). Biofertilizers in agriculture: An overview on concepts, strategies and effects on soil microorganisms. Adv. Agron. 162, 31–87. doi: 10.1016/bs.agron.2020.02.001
Megersa, H. G. (2017). Propagation methods of selected horticultural crops by specialized organs: review. J. Hortic. 4, 100198. doi: 10.4172/2376-0354.1000198
Ng, K. W., Cao, Z. J., Chen, H. B., Zhao, Z. Z., Zhu, L., and Yi, T. (2018). Oolong tea: A critical review of processing methods, chemical composition, health effects, and risk. Crit. Rev. Food Sci. 58, 2957–2980. doi: 10.1080/10408398.2017.1347556
Ni, H., Jing, X., Xiao, X., Zhang, N., Wang, X., Sui, Y., et al. (2021). Microbial metabolism and necromass mediated fertilization effect on soil organic carbon after long-term community incubation in different climates. ISME J. 15, 2561–2573. doi: 10.1038/s41396-021-00950-w
Nybom, H. and Lācis, G. (2021). Recent large-scale genotyping and phenotyping of plant genetic resources of vegetatively propagated crops. Plants. 10, 415. doi: 10.3390/plants10020415
Ortiz, A. and Sansinenea, E. (2022). The role of beneficial microorganisms in soil quality and plant health. Sustainability 14, 5358. doi: 10.3390/su14095358
Paśmionka, I. B., Bulski, K., and Boligłowa, E. (2021). The Participation of microbiota in the transformation of nitrogen compounds in the soil—A review. Agronomy. 11, 977. doi: 10.3390/agronomy11050977
Pradhan, D. R., Pattanaik, S., Kishore, J., and Gunanathan, C. (2020). Cobalt-catalyzed acceptorless dehydrogenation of alcohols to carboxylate salts and hydrogen. Org. Lett. 22, 1852–1857. doi: 10.1021/acs.orglett.0c00193
Rashmi, H. B. and Negi, P. S. (2020). Phenolic acids from vegetables: A review on processing stability and health benefits. Food Res. Int. 136, 109298. doi: 10.1016/j.foodres.2020.109298
Rivera-Pérez, A., Romero-González, R., and Frenich, A. G. (2022). Fingerprinting based on gas chromatography-orbitrap high-resolution mass spectrometry and chemometrics to reveal geographical origin, processing, and volatile markers for thyme authentication. Food Chem. 393, 133377. doi: 10.1016/j.foodchem.2022.133377
Schnecker, J., Baldaszti, L., Gündler, P., Pleitner, M., Sandén, T., Simon, E., et al. (2023). Seasonal dynamics of soil microbial growth, respiration, biomass, and carbon use efficiency in temperate soils. Geoderma. 440, 116693. doi: 10.1016/j.geoderma.2023.116693
Shen, F. T. and Lin, S. H. (2021). Priming effects of cover cropping on bacterial community in a tea plantation. Sustainability. 13, 4345. doi: 10.3390/su13084345
Sikder, M. M. and Vestergard, M. (2020). Impacts of root metabolites on soil nematodes. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.01792
Singh, P. (2021). Propagation of 14 native prairie forbs by sexual and asexual methods. Native Plants J. 22, 345–354. doi: 10.3368/npj.22.3.345
Soares, J. C., Santos, C. S., Carvalho, S. M., Pintado, M. M., and Vasconcelos, M. W. (2019). Preserving the nutritional quality of crop plants under a changing climate: importance and strategies. Plant Soil. 443, 1–26. doi: 10.1007/s11104-019-04229-0
Sun, Y., Duan, C., Cao, N., Ding, C., Huang, Y., and Wang, J. (2022). Biodegradable and conventional microplastics exhibit distinct microbiome, functionality, and metabolome changes in soil. J. Hazard. Mater. 424, 127282. doi: 10.1016/j.jhazmat.2021.127282
Tan, X., Xie, H., Yu, J., Wang, Y., Xu, J., Xu, P., et al. (2022). Host genetic determinants drive compartment-specific assembly of tea plant microbiomes. Plant Biotechnol. J. 20, 2174–2186. doi: 10.1111/pbi.13897
Wang, X., Gao, Y., Zang, P., Zhang, G., Yang, X., and Liu, Q. (2024a). Decrease in beneficial bacteria and increase in harmful bacteria in Gastrodia seedlings and their surrounding soil are mainly responsible for degradation of Gastrodia asexual propagation. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1334958
Wang, Y. H., Hong, L., Li, J. J., Zhang, Q. X., Wang, A. Q., Lin, S. X., et al. (2024b). Analysis of growth inhibition of continuously planted Casuarina equisetifolia in relation to characteristic soil microbial functions and nutrient cycling. Appl. Soil Ecol. 202, 105607. doi: 10.1016/j.apsoil.2024.105607
Wang, Y., Li, J., Li, M., Jia, X., Cai, Y., Hu, M., et al. (2023). Effect of continuous planting on Casuarina equisetifolia rhizosphere soil physicochemical indexes, microbial functional diversity and metabolites. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1288444
Warren, C. R. (2020). Soil microbial populations substitute phospholipids with betaine lipids in response to low P availability. Soil Biol. Biochem. 140, 107655. doi: 10.1016/j.soilbio.2019.107655
Withers, E., Hill, P. W., Chadwick, D. R., and Jones, D. L. (2020). Use of untargeted metabolomics for assessing soil quality and microbial function. Soil Biol. Biochem. 143, 107758. doi: 10.1016/j.soilbio.2020.107758
Wu, C., Ma, Y., Wang, D., Shan, Y., Song, X., Hu, H., et al. (2022). Integrated microbiology and metabolomics analysis reveal plastic mulch film residue affects soil microorganisms and their metabolic functions. J. Hazard. Mater. 423, 127258. doi: 10.1016/j.jhazmat.2021.127258
Xiao, W., Zhang, Z., Wang, H., Han, G., Yan, Z. Y., and He, D. (2023). Recombination of endophytic bacteria in asexual plant Ligusticum chuanxiong Hort. caused by transplanting. PeerJ. 11, e15579. doi: 10.7717/peerj.15579
Xie, X. G., Zhang, F. M., Wang, X. X., Li, X. G., and Dai, C. C. (2019). Phomopsis liquidambari colonization promotes continuous cropping peanut growth by improving the rhizosphere microenvironment, nutrient uptake and disease incidence. J. Sci. Food Agr. 99, 1898–1907. doi: 10.1002/jsfa.9385
Yang, Z., Yao, B., Li, R., Yang, W., Dong, D., Ye, Z., et al. (2024). Systematic review assessing the effects of amendments on acidic soils pH in tea plantations. PeerJ. 12, e17653. doi: 10.7717/peerj.17653
Ye, J., Wang, Y., Lin, S., Wang, Y., Chen, P., Hong, L., et al. (2023). Metabolomics analysis of the effect of acidification on rhizosphere soil microecosystem of tea tree. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1137465
Ye, J., Zhang, Q., Li, M., Wang, Y., Jia, M., Hong, L., et al. (2024). Tea quality of the mysterious “Dahongpao Mother Tree”(Camellia sinensis). Foods. 13, 1548. doi: 10.3390/foods13101548
Yuan, H., Cao, G., Hou, X., Huang, M., Du, P., Tan, T., et al. (2022). Development of a widely targeted volatilomics method for profiling volatilomes in plants. Mol. Plant 15, 189–202. doi: 10.1016/j.molp.2021.09.003
Zhao, M., Zhao, J., Yuan, J., Hale, L., Wen, T., Huang, Q., et al. (2021). Root exudates drive soil-microbe-nutrient feedbacks in response to plant growth. Plant Cell Environ. 44, 613–628. doi: 10.1111/pce.13928
Keywords: tea tree, cuttings, soil metabolites, microorganisms, nutrient cycling
Citation: Cheng W, Zhang S, Wang Y, Hong L, Qiu M, Wang Y, Luo Y, Zhang Q, Wang T, Jia X, Wang H and Ye J (2025) Dahongpao mother tree affects soil microbial community and nutrient cycling by increasing rhizosphere soil characteristic metabolite content. Front. Plant Sci. 16:1508622. doi: 10.3389/fpls.2025.1508622
Received: 09 October 2024; Accepted: 30 April 2025;
Published: 26 May 2025.
Edited by:
Shixiao Yu, Sun Yat-sen University, ChinaReviewed by:
Islam A. Abd El-Daim, Aberystwyth University, United KingdomSaloni Sharma, University of Granada, Spain
Copyright © 2025 Cheng, Zhang, Wang, Hong, Qiu, Wang, Luo, Zhang, Wang, Jia, Wang and Ye. 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: Haibin Wang, dzEzNTk5MDg0ODQ1QHNpbmEuY29t; Jianghua Ye, amh5ZTE5OThAMTI2LmNvbQ==
†These authors have contributed equally to this work
‡ORCID: Xiaoli Jia, orcid.org.0000-0001-5076-0604