Abstract
Objective:
Continuous cropping obstacles and heavy chemical fertilizer use constrain Radix serratulae chinensis cultivation. This study evaluates the optimal concentration of plant growth-promoting rhizobacterium Bacillus amyloliquefaciens ART9 and investigates its underlying mechanisms via rhizosphere microbiome and host transcriptome analyses.
Methods:
Seedlings were treated with ART9 suspensions at OD600 of 0.4, 0.6, and 1.0. Growth parameters, soil nutrient contents, rhizosphere bacterial communities, and root transcriptomic profiles were assessed to compare treatment effects.
Results:
The OD600 0.6 treatment yielded the best performance, increasing plant height, stem diameter, and SPAD value by 48.56%, 14.09%, and 12.42%, respectively, over the control, and raising soil available nitrogen, phosphorus, and potassium by 177%, 179%, and 168%. Rhizosphere microbiota were significantly enriched with beneficial genera, notably Bacillus, Pseudomonas, and Sphingomonas. Transcriptomics revealed prominent enrichment in phenylpropanoid biosynthesis (map00940), plant hormone signal transduction (map04075), and circadian rhythm–plant (map04712) pathways. Upregulated genes included JAZ, ABF, PYL, PP2C, SPA1_2, and HY5, indicating coordinated regulation of hormone and light signaling.
Conclusion:
ART9 systematically promotes growth of Radix serratulae chinensis by reshaping the rhizosphere microbiome and modulating hormone/light signaling pathways, supporting its development as a sustainable biofertilizer for medicinal plant production.
1 Introduction
Radix serratulae chinensis(syn. Cimicifuga foetida) is a medicinal herb whose dried rhizome possesses exceptionally high medicinal value. In traditional Chinese medicine (TCM) clinical practice, Radix serratulae Chinensis is commonly used to treat wind-heat headache, toothache, oral ulcers, sore throat, rectal prolapse, and uterine prolapse (Zhang et al., 2024b). Modern pharmacological studies have further shown that Radix serratulae Chinensis is rich in two major classes of bioactive compounds—9, 19-cycloartane-type triterpenoid saponins and cinnamic acid derivatives—which exhibit anti-osteoporotic, antiviral, antitumor, antiallergic, and estrogen-like activities (Sun et al., 2017). Despite its medicinal value, traditional cultivation is often hindered by continuous cropping obstacles and soilborne diseases, which compromise yield and quality. Microbial fertilizers, which provide both biocontrol and growth promotion, offer a green and sustainable solution to these challenges, including excessive dependence on chemical pesticides and fertilizers, and inconsistent quality of medicinal materials.
Among microbial agents, plant growth-promoting rhizobacteria (PGPR), as a class of naturally occurring, safe, and highly effective beneficial microorganisms, have demonstrated tremendous potential for advancing the green and ecological cultivation of medicinal plants. PGPR can significantly promote the biosynthesis and accumulation of pharmaceutically active compounds through multiple mechanisms (Santoyo et al., 2016). For instance, via direct actions such as nitrogen fixation, phosphate solubilization, siderophore secretion, and phytohormone synthesis (e.g., indole-3-acetic acid), PGPR optimize root architecture and enhance water and nutrient uptake capacity, thereby providing sustained nutritional support to medicinal plants (; ). Previous studies have shown that PGPR inoculation can effectively increase biomass, root length, and dry matter accumulation in a wide range of medicinal plant species, laying the foundation for high-yield cultivation ().
Importantly, PGPR can markedly enhance the adaptability of medicinal plants to abiotic stresses. Under adverse conditions such as salinity, drought, or heavy metal contamination, PGPR can activate the plant antioxidant enzyme system (e.g., superoxide dismutase, catalase, and ascorbate peroxidase), promote the accumulation of osmotic regulators (e.g., proline and soluble sugars), and modulate ion homeostasis. These actions effectively maintain cell membrane stability and physiological function, ensuring normal growth and secondary metabolism in medicinal plants even under unfavorable environmental conditions (Patani et al., 2023; ). This property provides a feasible pathway for the ecological cultivation of authentic medicinal herbs on marginal lands. In terms of soil management, PGPR also play an active role: they can promote soil aggregate formation by secreting extracellular polysaccharides, thereby improving water and nutrient retention capacity; certain strains can also reduce the bioavailability of heavy metals through biosorption, metal chelation, and ACC deaminase activity, contributing to the green remediation of contaminated soils (; Oubohssaine et al., 2025). Furthermore, PGPR can induce systemic resistance in plants, suppress soilborne pathogens, and reduce reliance on chemical pesticides, thus safeguarding the safety and authenticity of Chinese medicinal materials at the source (; Wang et al., 2024).
Among PGPR, Bacillus amyloliquefaciens is particularly important because of its pronounced growth-promoting and stress-tolerance-enhancing effects in multiple crops. This bacterium secretes active substances such as indole-3-acetic acid (IAA) and spermidine, and induces upregulation of IPT gene expression in plant roots via secreted signaling molecules, thereby increasing endogenous cytokinin (CTK) biosynthesis (Upadhyay et al., 2022). These actions promote seed germination, root development, and chlorophyll synthesis, significantly enhancing biomass and physiological activity in crops such as tomato, rice, and cucumber (Shao et al., 2015). Recent research has progressed from phenotypic observations to molecular investigations. Transcriptomic analyses have revealed that B. amyloliquefaciens GKT04, isolated from banana roots, exhibits high inhibitory activity against Fusarium oxysporum f. sp. cubense race 4 (FOC4), demonstrating its considerable potential for the biocontrol of Fusarium wilt (Tian et al., 2021). B. amyloliquefaciens SN16–1 activates defense responses in tomatoes via systemic acquired resistance, which is dependent on the salicylic acid signaling pathway, thereby helping the host resist pathogenic fungi (Zhao et al., 2021). Similarly, in maize and cucumber systems, integrated multi-omics analyses have successfully linked differentially expressed genes with the synthesis of secondary metabolites, elucidating the molecular networks underlying environmental stress responses, environmental sensing, metabolic regulation and host–microbe interactions (Tian et al., 2021; Xu et al., 2020; Oleńska et al., 2020). Nevertheless, single-omics approaches are insufficient to comprehensively resolve the complex regulatory networks that govern PGPR–plant interactions. Given that the growth-promoting effects of B. amyloliquefaciens involve multidimensional responses encompassing hormonal regulation, metabolic reprogramming, and immune activation, there is an urgent need to integrate transcriptomic and microbiome datasets at multiple levels to systematically elucidate the underlying growth-promoting mechanisms of this bacterium. Multi-omics integrative analyses not only enable the precise identification of key functional genes and metabolic pathways but also clarify the intrinsic mechanisms governing the synthesis of plant secondary metabolites, stress resistance, and yield formation (Singh et al., 2022; Yang et al., 2022).
Unlike previous studies that merely reported growth promotion by Bacillus amyloliquefaciens, the present study offers three distinct novelties: (1) first determination of the optimal application concentration (OD600=0.6) of Bacillus amyloliquefaciens ART9 for promoting Radix serratulae chinensis growth, providing a key parameter for standardized field application; (2) first integration of rhizosphere microbiome and root transcriptome to systematically reveal the cascade regulatory network of ‘strain concentration–rhizosphere microecology–growth response’; (3) first discovery that ART9 promotes medicinal plant growth through coordinated regulation of phytohormone signaling and circadian rhythm pathways, a mechanism not yet reported in related species. This study provides a theoretical foundation and technical support for the development of green cultivation technologies for Radix serratulae chinensis, thereby promoting its standardized and industrialized development.
2 Materials and methods
2.1 Materials and experimental design
The plant growth-promoting strain Bacillus amyloliquefaciens ART9 used in this experiment was preserved and supplied by the Key Laboratory of Traditional Chinese Medicine Agriculture, Changsha Medical University. Strain ART9 was molecularly identified by 16S rRNA gene sequencing (GenBank accession: PV174444) and whole-genome draft analysis (BioProject: PRJNA1475832). Phylogenetic analysis showed 99.98% 16S rRNA sequence similarity to Bacillus velezensis reference strain (GCF_001461825.1), with an average nucleotide identity (ANI) of 98.5% to Bacillus amyloliquefaciens DSM7T. The strain was cultured to logarithmic phase and adjusted to OD600 of 0.4, 0.6, and 1.0, corresponding to CFU of (1.2 ± 0.3)×106, (2.5 ± 0.4)×106, and (4.8 ± 0.6)×106 CFU mL-1, respectively. The selected OD600 values were based on preliminary experiments: OD600<0.4 showed no significant growth-promoting effect, while OD600 >1.0 caused excessive cell density potentially leading to competitive inhibition. The strain was aerobically cultured in LB liquid medium at 28 °C with shaking at 180–220 rpm until reaching the logarithmic growth phase (OD600 = 0.4–1.0) for subsequent use. Uniform 30-day-old Radix serratulae chinensis seedlings were collected from the planting base in Guidong County, Chenzhou City, Hunan Province. After conventional field soil preparation, robust seedlings with consistent growth were transplanted. Root irrigation treatments with bacterial suspensions at OD600 values of 0.4, 0.6, and 1.0 were established, with 20 mL applied per seedling in three separate applications; sterile water was used as the control. Each treatment included three independent biological replicates, with 7 uniform seedlings planted within each replicate (21 seedlings in total per treatment). All plants were managed under conventional field conditions without the application of chemical fertilizers or pesticides. After treatment, agronomic traits such as plant height, stem diameter, and root morphology were determined according to previously reported methods to evaluate the growth-promoting effect of the strain.
2.2 Determination of agronomic traits and soil physicochemical properties
Plant height and stem diameter were directly measured using a digital Vernier caliper with a precision of 0.01 mm. Leaf SPAD values were measured at the same leaf position using a TYS-4N chlorophyll meter. Soil physicochemical properties were determined using standard laboratory protocols: pH (soil:water ratio 2.5:1, pH meter method); organic matter (potassium dichromate oxidation-external heating method); alkali-hydrolyzable nitrogen (alkaline hydrolysis diffusion method); available phosphorus (sodium bicarbonate extraction-molybdenum antimony anti-colorimetric method); available potassium (ammonium acetate extraction-flame photometry). Each sample was measured in triplicate, and results are presented as mean ± SD.
2.3 Rhizosphere soil microbiome sequencing and analysis
Rhizosphere soil samples were collected using the five-point sampling method, with three to six biological replicates per group. Total microbial DNA was extracted using a soil microbial total DNA extraction kit, and DNA integrity and purity were verified by 1% agarose gel electrophoresis and a NanoDrop spectrophotometer, respectively (Versmessen et al., 2024). The V3–V4 hypervariable region of the 16S rRNA gene was PCR-amplified using barcode-labeled specific primers (Kim et al., 2023). Qualified PCR products were used for library construction and subsequently subjected to paired-end sequencing on the Illumina NovaSeq 6000 platform. Raw sequencing reads were quality-controlled and merged using fastp (v0.23.4) and FLASH (v1.2.11), with low-quality sequences filtered out and chimeras removed. Operational taxonomic units (OTUs) were clustered at 97% sequence identity using USEARCH (v11). This 97% OTU threshold was adopted to align with the SILVA 138.2 database annotation standard and facilitate comparison with published PGPR rhizosphere studies. Taxonomic annotation was performed using the RDP Classifier (v2.11) against the SILVA 138.2/16S_bacteria database with a confidence threshold of 70%. After removing chloroplast and mitochondrial sequences, OTU abundance data were rarefied to 49, 248 sequences per sample. Alpha- and beta-diversity analyses were performed using Mothur (v1.30.2), differential species analysis was conducted in R (v3.3.1), and functional prediction was performed using PICRUSt2 (v2.2.0) with default parameters.
2.4 Root transcriptome sequencing and qRT-PCR validation of Radix serratulae Chinensis
On day 28 after bacterial suspension treatment, three uniformly growing Radix serratulae chinensis plants were sampled from the treatment (OD600 = 0.6) and control groups. Three biological replicates were set per treatment (each replicate consisted of roots pooled from 7 seedlings). Although the sample size is modest, based on previous transcriptomic studies in medicinal plants (Li et al., 2020), three biological replicates combined with stringent DEG thresholds (|log2FC|≥1, FDR<0.05) are sufficient to identify major regulatory pathways. Future studies should increase sample size to validate these findings. Root samples were immediately frozen in liquid nitrogen and stored at −80 °C for total RNA extraction, quality assessment, and Illumina NovaSeq 6000 PE150 sequencing, with ≥6 Gb clean data generated per sample (Modi et al., 2021). Raw reads were aligned to the Radix serratulae chinensis reference genome(https://ngdc.cncb.ac.cn/gwh/Assembly/486/show), but the alignment rate was only 0.31%–0.78%; we therefore used Trinity for de novo assembly. After quality control of raw sequencing data, Q20≥96.5%, Q30≥91.2%. A total of 28, 490 unigenes were assembled with N50 = 1, 256 bp. The average mapping rate of unigenes back to the transcriptome was 82.3%. BUSCO assessment showed 78.3% completeness. Annotation efficiency: NR 68.7%, KEGG 51.2%, GO 58.3%. Unigenes were annotated against the NR, SwissProt, KEGG, and GO databases. Differentially expressed genes (DEGs) were screened using DESeq2 with thresholds of |log2FC| ≥ 1 and FDR < 0.05 (Surachat et al., 2022). Venn diagram and principal component analysis (PCA) were subsequently performed, followed by KEGG enrichment analysis of DEGs related to growth, nutrient absorption, hormone signaling, secondary metabolism, and stress-response pathways, revealing that ART9 activated multiple growth-related pathways and providing candidate genes for further mechanistic investigation (Li et al., 2020). To verify the reliability of the transcriptome sequencing results, key DEGs associated with growth, hormone signaling, and metabolic pathways were selected for qRT-PCR validation. All primer pairs exhibited amplification efficiencies between 90% and 110% (Supplementary Table 1), with single-peak melting curves confirming amplification specificity and no primer-dimer formation. Reference genes GAPDH and 18S rRNA were evaluated using geNorm (stability M < 0.5) and used as internal controls (Supplementary Table 1). Reverse-transcribed cDNA was used as the template for qRT-PCR analysis. Each sample was run in three technical replicates. Relative gene expression levels were calculated using the 2-ΔΔCt method (Lv et al., 2021), and the consistency between qRT-PCR results and transcriptome data was analyzed, thereby providing data support and candidate genes for further elucidating the molecular mechanisms by which ART9 regulates the growth of Radix serratulae chinensis.
2.5 Integration of phenotypic and microbiome data
To systematically elucidate the coordinated changes among Radix serratulae chinensis phenotypic traits, soil physicochemical factors, and rhizosphere microbial community structure, a distance-based redundancy analysis (db-RDA) was employed for multidimensional data integration. Traditional redundancy analysis (RDA) imposes stringent requirements regarding data normality and linear relationships, whereas microbiome data are typically characterized by high dimensionality and sparsity. Db-RDA effectively overcomes these limitations by incorporating a distance matrix, thereby enabling a more accurate quantification of the extent to which environmental factors explain community variation.
2.6 Statistical analysis
All experiments were conducted using a completely randomized design. For agronomic phenotype measurements, the mean value of 7 seedlings within each biological replicate was taken as a single data point; final data are presented as mean ± standard deviation (SD) of three independent biological replicates (n = 3). All data were first subjected to Shapiro-Wilk normality test and Levene’s test for homogeneity of variances. Data meeting normality and homoscedasticity were analyzed by one-way ANOVA; otherwise, the Kruskal-Wallis non-parametric test was used. For microbiome differential analysis, the Benjamini-Hochberg method was applied to control the false discovery rate (FDR < 0.05). Following analysis of variance (ANOVA), Tukey’s HSD test was used for multiple comparisons, with statistical significance set at p ≤ 0.05. Data organization was performed using Microsoft Excel, and all statistical analyses and graphical processing were conducted using GraphPad Prism 10.1.2 software based on a one-way ANOVA.
3 Results
3.1 Concentration-dependent growth-promoting effects of ART9 on Radix serratulae Chinensis and soil properties
Bacillus amyloliquefaciens ART9 was used as the test strain and cultured in LB medium to OD600 values of 0.4, 0.6, and 1.0, representing three concentration gradients, to investigate its effects on the growth and development of Radix serratulae chinensis. The results showed that all ART9 treatments at different concentrations significantly promoted the growth of Radix serratulae chinensis, with the OD600 = 0.6 treatment group showing the optimal effect. Growth indices, including plant height, stem diameter, and SPAD values, as well as soil nitrogen (N), phosphorus (P), and potassium (K) contents (Table 1), were all significantly higher in the OD600 = 0.6 treatment group than in the other treatment groups and the control (CK), demonstrating a clear concentration-dependent peak response (Figure 1).
Table 1
| D600 | Plant height (cm) | Stem diameter (mm) | SPAD | N (mg kg-1) | P (mg kg-1) | K (mg kg-1) |
|---|---|---|---|---|---|---|
| Ck | 22.18 ± 3.10 | 5.11 ± 0.57 | 33.89 ± 3.18 | 7.90 ± 1.45 | 11.10 ± 2.21 | 28.20 ± 5.38 |
| 0.4 | 28.82 ± 3.34 | 5.78 ± 0.83 | 37.18 ± 2.14 | 14.80 ± 2.96 | 22.80 ± 6.66 | 51.20 ± 12.20 |
| 0.6 | 32.95 ± 3.96 | 5.83 ± 1.00 | 38.10 ± 2.85 | 21.90 ± 5.66 | 31.00 ± 7.95 | 75.50 ± 21.50 |
| 1.0 | 30.90 ± 4.68 | 5.22 ± 1.26 | 34.42 ± 4.00 | 14.60 ± 4.25 | 19.20 ± 4.12 | 46.90 ± 13.52 |
Effects of Bacillus amyloliquefaciens with different OD600 values on indexes of Radix serratulae Chinensis.
Data are presented as mean ± SD (n = 3 biological replicates, each with 7 seedlings). Different lowercase letters within the same column indicate statistically significant differences among treatments at p ≤ 0.05 (Tukey’s HSD test).
Figure 1
In terms of plant height, the OD600 = 0.6 ART9 treatment group attained the maximum value, representing a significant increase of 48.56% compared with that of CK. As the ART9 concentration increased from OD600 = 0.4 to 0.6, plant height showed a sustained and significant upward trend; however, when the concentration was further increased to OD600 = 1.0, the plant height declined. These results indicate that OD600 = 0.6 is the optimal concentration for promoting plant height growth in Radix serratulae chinensis. Regarding stem diameter, the OD600 = 0.6 treatment group ranked first among all treatment groups, showing a 14.09% increase compared with CK. This indicates that the ART9 bacterial suspension exerts a significant promoting effect on stem diameter in Radix serratulae chinensis, with the OD600 = 0.6 treatment producing the best results.
The trend in SPAD values, which reflect the relative chlorophyll content of the plant leaves, was consistent with that observed for plant height. The OD600 = 0.6 treatment group reached the highest SPAD value, representing an increase of 12.42% compared with that of CK. This result indicates that the OD600 = 0.6 ART9 bacterial suspension exerts the strongest promoting effect on chlorophyll synthesis in Radix serratulae chinensis, thereby effectively enhancing the photosynthetic potential of the plants.
3.2 Response of rhizosphere microbial community structure to ART9 concentration
Principal coordinate analysis (PCoA) at the OTU level revealed a clear separation trend between the microbial community structures of the control (CK) and treatment (S) groups (Figure 2). PC1 and PC2 explained 32.73% and 19.79% of the total variance, respectively, with a cumulative explanatory rate of over 50%. All CK group samples were concentrated in the negative region of the PC1 axis, whereas all S group samples were distributed in the positive region, with no overlap between the two groups. This indicates that ART9 treatment significantly altered the overall microbial community structure. Furthermore, the CK group sample points were tightly clustered, reflecting high within-group stability; although the S group sample points were slightly more dispersed, the between-group differences remained substantially greater than the within-group differences.
Figure 2
Venn diagram analysis further revealed differences in species composition between the two groups (Figure 2). A total of 3730 OTUs were detected in the CK group, of which 894 (19.80%) were unique to this group, whereas the S group yielded 3622 OTUs, including 786 (17.41%) group-specific OTUs. The two groups shared 2836 OTUs, accounting for 62.80% of the total OTU count, indicating that most of the core microbial taxa were common to both groups. However, the presence of unique OTUs suggests that ART9 treatment not only resulted in the loss or reduced abundance of certain species present in the CK group but also facilitated the emergence or significant enrichment of new species in the S group, reflecting a treatment-specific response of the microbial community at the species level.
Analysis of the rhizosphere microbial community composition (Figure 2C) showed that compared with the CK group, multiple microbial taxa were markedly enriched or stably detected in the ART9 treatment group (Figure 2). Specifically, Sphingomonas, Bradyrhizobium, Pseudomonas, Nitrospira, Chitinophagaceae, Gemmatimonadaceae, Cupriavidus, and the Terrabacteria group were enriched in the S group.
3.3 Transcriptional response of Radix serratulae Chinensis root gene expression to ART9 treatment
3.3.1 Gene expression profiling of Radix serratulae Chinensis
To analyze the changes in the gene expression profiles of Radix serratulae chinensis following inoculation with B. amyloliquefaciens ART9, Venn diagram analysis was first performed to examine the gene distribution between the two groups. The CK group contained considerably more unique genes (8186) than the R group (5669), with a ratio of 1.44:1. The two groups shared 24, 824 genes (Figure 3A), indicating that although gene expression exhibited extensive commonality, both groups possessed distinct and group-specific expression profiles. Principal component analysis (PCA) showed that PC1 and PC2 explained 30.04% and 25.18% of the expression variance, respectively. The CK and R group samples were clearly separated in the score plot (Figure 3B), confirming that ART9 treatment exerted a significant effect on the gene expression patterns of Radix serratulae chinensis, with between-group differences substantially exceeding within-group differences.
Figure 3
Volcano plot analysis was used to identify differentially expressed genes (DEGs). A total of 1753 significant DEGs were identified (|log2FC| ≥ 1, p < 0.05), including 778 upregulated genes (44.4%) and 975 downregulated genes (55.6%), with the number of downregulated genes being 1.25-fold higher than that of the upregulated genes (Figure 3C). To further characterize the most representative transcriptional changes, key candidate DEGs were summarized in Table 2. These results indicate that ART9 inoculation induced widespread transcriptional reprogramming in Radix serratulae chinensis, with downregulated genes accounting for a greater proportion than upregulated genes.
Table 2
| Metabolic pathway | Key upregulated genes | Key downregulated genes |
|---|---|---|
| Phenylpropanoid biosynthesis | None | PAL (K10775), CSE (K12356), CSE (K18368), COMT (K13066), E2.1.1.104 (K00588) |
| Plant hormone signal transduction | PP2C (K14497), ARR-B (K14491), PYL (K14496), ABF (K14432), JAZ (K13464), YUCCA (K11816) | AUX1 (K13946), ETR (K14509), EBF1_2 (K14515), ERF1 (K14516), CYCD3 (K14505), MYC2 (K13422), TGA (K14431), glnA (K01915), CALM (K02183) |
| Starch and sucrose metabolism | E2.7.1.4 (K00847), otsA (K00697), otsB (K01087) | pgm (K01835), AMY (K01176) |
| Circadian rhythm – plant | SPA1_2 (K16240), GI (K12124), HY5 (K16241), FKF1 (K12116) | CSNK2B (K03115) |
Key candidate genes identified in Radix serratulae chinensis roots.
3.3.2 Relative gene expression analysis
To further elucidate the growth-promoting mechanisms of Radix serratulae chinensis following inoculation with B. amyloliquefaciens ART9, KEGG pathway enrichment analysis was performed on the DEGs identified from the transcriptome data. Figure 3D presents the top 20 enriched pathways, revealing the metabolic changes induced by ART9 inoculation. Notably, the phenylpropanoid biosynthesis (map00940), plant hormone signal transduction (map04075), and circadian rhythm–plant (map04712) pathways were significantly enriched, with p-values of 0.001, 0.008, and 0.007, respectively. These pathways play critical roles in promoting root cell wall loosening, activating defense responses, and optimizing growth rhythmicity.
In addition, starch and sucrose metabolism (map00500), nitrogen metabolism (map00910), fatty acid biosynthesis (map00061), α-linolenic acid metabolism (map00592), and glutathione metabolism (map00480) pathways were significantly enriched. The coordinated activation of these pathways, all of which exhibited significant upregulation trends, collectively established a comprehensive growth-promoting network encompassing “energy supply–nitrogen utilization–membrane homeostasis–antioxidant defense, ” further highlighting the importance of organic compound synthesis in plant growth. These findings comprehensively elucidate the metabolic reprogramming mechanisms occurring in the roots of Radix serratulae chinensis following inoculation with B. amyloliquefaciens ART9 and identify the core pathway networks regulating plant growth and development.
3.3 qPCR validation
To further substantiate the growth-promoting mechanisms of Radix serratulae chinensis following inoculation with B. amyloliquefaciens ART9, qPCR was performed to validate six DEGs identified in the root transcriptome analysis. Compared with the control group, the inoculation treatment group exhibited significantly induced expression of all candidate genes, consistent with the trends observed in transcriptome sequencing data (Figure 4).
Figure 4
In the phytohormone signaling pathway, the jasmonate negative regulator JAZ (K13464) was upregulated approximately 2.2-fold, whereas the abscisic acid (ABA) signaling components ABF (K14432) and PYL (K14496) were upregulated approximately 9-fold and 3.5-fold, respectively, indicating the activation of both jasmonate and ABA signaling pathways. Protein phosphatase PP2C (K14497) exhibited the most pronounced upregulation, with an approximately 15-fold increase, suggesting its potential involvement in signal transduction through dephosphorylation. In addition, the core light signaling factors SPA1_2 (K16240) and HY5 (K16241) were significantly upregulated by approximately 5.3-fold, indicating the involvement of photomorphogenic signaling in the response to inoculation.
The qPCR results validated the reliability of the transcriptome data and confirmed JAZ, ABF, PYL, PP2C, SPA1_2, and HY5 as key genes responsive to inoculation. Among these, PP2C, ABF, SPA1_2, and HY5 were identified as the core regulatory factors. These findings provide an important experimental foundation for elucidating the molecular response mechanisms of Radix serratulae chinensis to ART9 inoculation and for the subsequent identification and application of key functional genes in the future.
3.4 db-RDA analysis of rhizosphere microbial communities at the OTU level in relation to environmental factor
To explore the association between rhizosphere microbial communities and phenotypic traits of Radix serratulae chinensis, distance-based redundancy analysis (db-RDA) was performed at the OTU level, and the results are shown in Figure 5. The ordination axes CAP1 and CAP2 explained 31.02% and 16.66% of the microbial community variation, respectively, with a cumulative explanatory rate of 47.68%, effectively capturing the major trends in community structure variation.
Figure 5
In terms of sample distribution, rhizosphere soil microbial communities exhibited clear intergroup separation between the PGPR treatment (s1) and control (ck1) groups. Samples from the s1 group were predominantly concentrated in the right region of the ordination space, whereas those from the ck1 group were distributed in the left region, indicating that PGPR treatment significantly altered the bacterial community structure in the rhizosphere soil of Radix serratulae chinensis, resulting in marked differences in community composition between the two groups.
Further analysis of the relationships between environmental factors and community structure revealed that the vectors representing plant height and soil nitrogen content (N) were highly consistent with the distribution of s1 group samples, with an extremely small angle between them, indicating a strong and positive correlation. In contrast, the vector representing chlorophyll content (SPAD) was oriented toward the distribution of the ck1 group samples. These results suggest that the divergence in rhizosphere microbial communities induced by PGPR treatment was primarily driven by changes in plant height and soil nitrogen content in Radix serratulae chinensis, and that these two factors exerted synergistic effects on community structure. The microbial communities in the control group were more closely associated with environmental factors related to chlorophyll synthesis.
The db-RDA results clearly demonstrated associations between rhizosphere microbial communities and key growth traits of Radix serratulae chinensis under PGPR treatment, providing correlative ecological evidence that may help elucidate the mechanisms by which PGPR promotes the growth of Radix serratulae chinensis through the regulation of rhizosphere microbial communities. However, these associations do not establish causation.
4 Discussion
PGPR are a class of beneficial bacteria that colonize plant roots and promote the growth of medicinal plants, including Radix serratulae Chinensis, through diverse direct and indirect mechanisms (). Screening highly effective PGPR strains and understanding their regulatory mechanisms are essential for advancing green cultivation and ensuring stable, high-quality yields of medicinal materials (Wahab et al., 2024; Zhang et al., 2026; ). Although numerous studies have demonstrated that PGPR can enhance crop growth by optimizing soil microbial communities and modulating root gene expression (Miljaković et al., 2020; Li et al., 2021), systematic investigation of the optimal concentration of Bacillus amyloliquefaciens ART9 for Radix serratulae Chinensis cultivation has been lacking. In the present field experiment, Radix serratulae chinensis seedlings treated with ART9 suspensions at OD600 = 0.6 exhibited significantly superior vigor compared with the other treatments (Figure 1). This study is the first to establish OD600 = 0.6 as the optimal application concentration of ART9 for promoting Radix serratulae Chinensis growth, thereby filling an important technical gap and providing a reliable basis for standardized inoculant application. These findings may directly guide farmers in applying bacterial irrigation during the seedling or transplanting stages to achieve maximal growth promotion with minimal input, thereby reducing cultivation costs and enhancing economic returns.
In addition to optimizing inoculant concentration, this study explored the molecular basis of growth promotion by analyzing genes involved in plant hormone signal transduction and circadian rhythm regulation. qRT-PCR validated the expression levels of six key genes: JAZ, ABF, PYL, PP2C, SPA1_2, and HY5 (Figure 4). Among these, JAZ, ABF, PYL, and PP2C are upregulated genes involved in the plant hormone signal transduction pathway, whereas SPA1_2 and HY5 are critical components of the circadian rhythm–plant pathway. The latter two form a dynamic and phase-complementary signaling module. HY5 acts as a positive regulator of light signaling and is activated at dawn to promote photosynthetic activity (Wang et al., 2021; Yadav et al., 2020; Jing et al., 2020). Conversely, SPA1 functions as a negative regulator, accumulating at dusk and suppressing HY5 activity through protein degradation and phosphorylation, thereby facilitating the transition to nocturnal metabolism. Their synergistic action enables plants to optimally allocate resources across the day–night cycle, balancing growth and stress responses, a mechanism particularly beneficial for Radix serratulae Chinensis growth. However, it should be noted that the transcriptomic data in this study reveal correlative changes in gene expression rather than direct causal mechanisms. Future studies should employ genetic manipulation (e.g., CRISPR knockout or overexpression of key genes such as HY5 and PP2C) or chemical intervention (e.g., ABA inhibitors) to validate the functional causality of these genes.
Within the plant hormone signal transduction pathway, the upregulation of JAZ, a negative regulator of jasmonic acid (JA) signaling, likely represents a feedback inhibition mechanism that prevents excessive activation of JA-mediated defence responses, thereby optimizing the growth–defence balance, consistent with previous studies (Xu et al., 2024; Li et al., 2025). In the abscisic acid (ABA) pathway, the marked upregulation of the transcription factor ABF and the receptor PYL indicates that ART9 treatment enhances the capacity of plants for ABA perception and response, thereby contributing to stress tolerance (). Notably, the protein phosphatase PP2C (K14497), a core negative regulator of ABA signaling, was induced approximately 15-fold (). PP2C upregulation (15-fold) appears contradictory to ABA signaling activation. A plausible explanation is that under non-stress conditions, ART9 treatment induces strong negative feedback expression of PP2C, thereby relieving ABA-mediated growth inhibition and promoting plant growth. This represents a typical ‘brake’ mechanism rather than direct activation of ABA signaling. Thus, the observed upregulation of ABF and PYL may reflect enhanced basal activity of ABA signaling, while high PP2C expression represents a balancing negative feedback. This dynamic reprogramming of hormone signaling likely enables plants to adapt to fluctuating environmental conditions (Pan et al., 2024). In the circadian rhythm pathway, the synchronous high expression of SPA1_2 and HY5 indicates that ART9 treatment optimizes light adaptability and coordinates growth metabolic rhythms, corroborating previous findings regarding the light signaling regulation of plant development (Mankotia et al., 2024; ). Collectively, these results confirm that ART9 specifically activates key genes involved in hormone signaling and circadian rhythm pathways, thereby providing a molecular basis for vigorous Radix serratulae Chinensis seedling growth.
Beyond molecular regulation, ART9 application profoundly affected the rhizosphere bacterial community. In the present study, the S group showed marked enrichment of Sphingomonas, Bradyrhizobium, Pseudomonas, Nitrospira, Chitinophagaceae, Gemmatimonadaceae, Cupriavidus, and the Terrabacteria group, all of which are dominant taxa commonly associated with plant growth-promoting functions (Figure 2C) (; ; Comeau et al., 2021; Zhang et al., 2024a). At the OTU level, db-RDA revealed a clear separation between the PGPR treatment (s1) and control (ck1) groups, indicating that PGPR inoculation significantly reshaped the rhizosphere bacterial community through rhizosphere colonization, niche competition, and modulation of microbial interaction network (Figure 5). Among the environmental factors analyzed, plant height and soil nitrogen content (N) emerged as the key drivers, with their vectors closely aligned with the distribution of s1 samples and a strong correlation between them. These findings suggest that PGPR-enriched bacterial communities directly enhance nitrogen uptake and utilization in Radix serratulae Chinensis, thereby promoting plant height growth. Enriched taxa such as Bradyrhizobium, Pseudomonas, and Cupriavidus possess nitrogen-fixing, phosphate-solubilizing, and growth-promoting substance-secreting capabilities, providing a microbiological basis for improved nitrogen transformation and assimilation. However, functional predictions in this study were based solely on PICRUSt2, which is computational inference. To validate the predicted microbial functions, future studies should perform metagenomic sequencing and biochemical assays for key functions (e.g., nitrogen fixation via acetylene reduction assay, phosphate solubilization via molybdenum antimony anti-colorimetric method). In contrast, the control group samples were more closely associated with SPAD values (chlorophyll content), suggesting that in the absence of PGPR treatment, the rhizosphere community was dominated by taxa associated with chlorophyll synthesis and photosynthetic efficiency. Although the ck1 group exhibited relatively high SPAD values, the s1 group achieved superior growth vigor through microbial community restructuring, indicating that PGPR promoted growth not by directly enhancing photosynthesis but by optimizing nutrient uptake and hormone regulation. However, functional predictions in this study were based solely on PICRUSt2, which is computational inference. To validate the predicted microbial functions, future studies should perform metagenomic sequencing and biochemical assays for key functions (e.g., nitrogen fixation via acetylene reduction assay, phosphate solubilization via molybdenum antimony anti-colorimetric method).
In this analysis, db-RDA analysis revealed associations between PGPR treatment promoted the enrichment of specific functional microbial groups, as well as their correlation with key phenotypic traits, such as plant height and nitrogen content. However, these associations do not directly imply causation. This positive community–phenotype association represents an important ecological mechanism underlying PGPR-mediated growth promotion (). Future studies should isolate the enriched dominant functional strains and conduct single-strain or synthetic community re-inoculation experiments in gnotobiotic systems to verify their individual contributions to plant height and nitrogen uptake, thereby establishing causality and identifying core functional strains and their molecular mechanisms. The cascade regulatory network of strain concentration, rhizosphere microecology, and growth response is proposed to explain the growth-promoting effect.
Collectively, this study makes three principal contributions to the literature. First, it identified OD600 = 0.6 as a candidate optimal concentration of ART9 for Radix serratulae Chinensis growth under the controlled plot conditions tested, providing a theoretical reference for future concentration optimization. Second, it reveals a dual regulatory mechanism at the molecular level, whereby ART9 coordinately modulates plant hormone signaling and circadian rhythm. Third, it integrates rhizosphere microbiome remodeling with phenotypic and transcriptomic data, suggesting that ART9 promotes growth by constructing a nitrogen utilization-centered functional microbial community rather than by directly increasing chlorophyll content (Timofeeva et al., 2022; Kalayu, 2019; Lidbury et al., 2016). Figure 6 presents a graphical summary of the overall mechanism. As shown in the graphical abstract, ART9 root irrigation at the optimal concentration reshapes the rhizosphere microbiome by enriching N-fixing and P-solubilizing genera (Sphingomonas, Bradyrhizobium, Pseudomonas, Nitrospira, Cupriavidus, Chitinophagaceae), which increases soil available nutrients. This microbiome remodeling is coupled with root transcriptional reprogramming that co-activates hormone signaling, including the ABA pathway (PYL, ABF, and PP2C as a feedback “brake”) and the JA pathway (JAZ), as well as light and circadian modules (HY5 and SPA1_2), ultimately promoting plant growth through a cascade regulatory network of “strain concentration- rhizosphere microecology-growth response”. The present study was carried out under controlled field plot experimental conditions to preliminarily explore the growth-promoting mechanism of strain ART9 on Radix serratulae chinensis. The current findings provide fundamental theoretical data and technical reference for subsequent large-scale field verification, strain agent development, and field application trials. Further multi-location, multi-year field validation is urgently required before agricultural popularization and commercial utilization.
Figure 6
Nevertheless, several limitations of this study should be acknowledged. The optimal application strategy requires further validation across different soil types, climatic conditions and growth stages. The specific molecular signals through which ART9 regulates host gene expression, either via secreted metabolites or alterations in the rhizosphere microenvironment, remain unclear. In addition, this study did not perform metabolomic analysis; therefore, statements regarding ‘metabolic reprogramming’ are primarily based on transcriptomic inference. In addition, this study focused on non-destructive growth parameters (plant height, stem diameter, and SPAD values) and did not collect destructive samples for biomass measurement (fresh weight and dry weight), which are important indicators for comprehensively evaluating growth-promoting effects. Future studies should include biomass measurements to provide a more complete assessment of ART9-mediated growth promotion. Furthermore, since Radix serratulae chinensis is a medicinal plant whose value depends primarily on the content of bioactive constituents (e.g., 9, 19-cycloartane-type triterpenoid saponins and cinnamic acid derivatives), this study focused only on growth promotion and did not quantify these active compounds. Moreover, the individual roles and interaction networks of enriched microbial taxa require experimental verification. To address these gaps, future studies should include multi-site and multi-year field trials, isolation and re-inoculation of dominant strains, identification of signaling molecules in sterile systems, untargeted metabolomics to directly detect changes in metabolites within the phenylpropanoid biosynthesis pathway (e.g., cinnamic acid, coumaric acid, monolignols), HPLC or LC-MS/MS quantification of key pharmacologically active components to assess whether ART9 treatment genuinely improves medicinal material quality and integrated transcriptomic–metabolomic analyses to comprehensively elucidate the global response network of Radix serratulae Chinensis roots.
5 Conclusions
This study systematically investigated the growth-promoting effects of Bacillus amyloliquefaciens ART9 on Radix serratulae chinensis and elucidated the underlying mechanisms through integrated microbiome and transcriptomic analyses. The optimal concentration of ART9 (OD600 = 0.6) was identified, which significantly enhanced plant height, stem diameter, chlorophyll content, and soil available N, P, and K. Multi-omics analyses revealed that ART9 reshaped the rhizosphere microbiome by enriching N-fixing and P-solubilizing genera (Sphingomonas, Bradyrhizobium, Pseudomonas, Nitrospira, Cupriavidus, and Chitinophagaceae), while root transcriptional reprogramming co-activated hormone signaling (ABA and JA pathways) and light/circadian modules (HY5 and SPA1_2). A cascade regulatory network of strain concentration-rhizosphere microecology -growth response is proposed to explain the growth-promoting effect. These findings provide a theoretical foundation and technical reference for the development of green cultivation technologies for Radix serratulae chinensis, though further multi-location, multi-year field validation is required before agricultural application.
Statements
Data availability statement
The raw sequencing data generated in this study have been deposited in the NCBI database under the following BioProject accession numbers: PRJNA1467762 (16S rRNA gene amplicon sequences), PRJNA1467520 (raw transcriptome data), and PRJNA1475832 (whole-genome draft of strain ART9).
Ethics statement
The plant material used in this study was commercially purchased. Experimental research on plants was conducted in compliance with relevant institutional, national, and international guidelines and legislation.
Author contributions
JW: Writing – original draft. JX: Writing – original draft. PF: Writing – original draft. JL: Data curation, Validation, Writing – original draft. BT: Validation, Writing – original draft, Data curation. YZ: Validation, Data curation, Writing – original draft. ZZ: Writing – review & editing. YF: Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This research was funded by the Hunan Natural Science Foundation (2026JJ90160), Science Research Project of Hunan Provincial Department of Education (24C0626), and the National Natural Science Foundation of China (NSFC 32272011).
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.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this 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.2026.1891133/full#supplementary-material
Supplementary Table 1qRT-PCR primers sequences used for validation of differentially expressed genes in Radix serratulae chinensis roots.
Supplementary Table 2KEGG pathway enrichment analysis of differentially expressed genes in Radix serratulae chinensis roots following ART9 treatment.
Supplementary Table 3Expression profiles and functional annotation of differentially expressed genes in Radix serratulae chinensis roots following ART9 treatment.
Supplementary Table 4OTU abundance table of rhizosphere microbial communities in ART9-treated and control groups.
References
1
AartiS.ParmarS.KumarA. (2024). Plant growth promoting rhizobacteria: their potential in sustainable agriculture. J. Sci. Innov. Nat. Earth4, 40–46. doi: 10.59436/jsianev4i2/270
2
AbubakarA. S.WuY.ChenF.ZhuA.ChenP.ChenK.et al. (2023). Comprehensive analysis of WUSCEL-related homeobox gene family in ramie (Boehmeria nivea) indicates its potential role in adventitious root development. Biology12, 1475. doi: 10.3390/biology12121475
3
ChenF.WangY.WangK.ChenJ.JinK.PengK.et al. (2023). Effects of Litsea cubeba essential oil on growth performance, blood antioxidation, immune function, apparent digestibility of nutrients, and fecal microflora of pigs. Front. Pharmacol.14, 1166022. doi: 10.3389/fphar.2023.1166022
4
ChenJ.WangY.WuY.ChuZ. (2024). Genome-wide identification and expression analysis of the PP2C gene family in Apocynum venetum and Apocynum hendersonii. BMC Plant Biol.24, 652. doi: 10.1186/s12870-024-05328-6
5
ChenS.WirthmuellerL.StauberJ.LoryN.HoltkotteX.LesonL.et al. (2016). The functional divergence between SPA1 and SPA2 in Arabidopsis photomorphogenesis maps primarily to the respective N-terminal kinase-like domain. BMC Plant Biol.16, 160. doi: 10.1186/s12870-016-0854-9
6
ChengM.ShiY.ChengY.HuH.LiuS.XuY.et al. (2024). Mulberry leaf polysaccharide improves cyclophosphamide-induced growth inhibition and intestinal damage in chicks by modulating intestinal flora, enhancing immune regulation and antioxidant capacity. Front. Microbiol.15, 1382639. doi: 10.3389/fmicb.2024.1382639
7
ComeauD.BalthazarC.NovinscakA.BouhamdaniN.JolyD. L.FilionM. (2021). Interactions between Bacillus spp. Pseudomonas spp. and Cannabis sativa promote plant growth. Front. Microbiol.12, 715758. doi: 10.3389/fmicb.2021.715758
8
DhimanI.NandniPoriaV.KumarS.YadavR.ShaikT.et al. (2026). Plant growth promoting rhizobacteria (PGPR) mediated amelioration of plant tolerance to abiotic stresses: drought, salinity, and heavy metals. Arch. Microbiol.208, 100. doi: 10.1007/s00203-025-04631-6
9
FidlerJ.GraskaJ.GietlerM.NykielM.PrabuckaB.Rybarczyk-PłońskaA.et al. (2022). PYR/PYL/RCAR receptors play a vital role in the abscisic-acid-dependent responses of plants to external or internal stimuli. Cells11, 1352. doi: 10.3390/cells11081352
10
Floc'hJ.-B.HamelC.LaterrièreM.TidemannB.St-ArnaudM.HijriM. (2021). Inter-kingdom networks of canola microbiome reveal Bradyrhizobium as keystone species and underline the importance of bulk soil in microbial studies to enhance canola production. Microb. Ecol.84, 1166–1181. doi: 10.1007/s00248-021-01905-6
11
GuptaR.KhanF.AlqahtaniF. M.HashemH. M.AhmadF. (2024). Plant growth–promoting rhizobacteria (PGPR) assisted bioremediation of heavy metal toxicity. Appl. Biochem. Biotechnol.196, 2928–2956. doi: 10.1007/s12010-023-04545-3
12
HeL.RenY.ZengW.WuX.ShenL.YuR.et al. (2021). Deciphering the endophytic and rhizospheric microbial communities of a metallophyte Commelina communis in different Cu-polluted soils. Microorganisms9, 1689. doi: 10.3390/microorganisms9081689
13
HeS.LiL.LvM.WangR.WangL.YuS.et al. (2024). PGPR: key to enhancing crop productivity and achieving sustainable agriculture. Curr. Microbiol.81, 377. doi: 10.1007/s00284-024-03893-5
14
JiaoX.TakishitaY.ZhouG.SmithD. L. (2021). Plant-associated rhizobacteria for biocontrol and plant growth enhancement. Front. Plant Sci.12, 634796. doi: 10.3389/fpls.2021.634796
15
JingY.GuoQ.LinR. (2020). The SNL-HDA19 histone deacetylase complex antagonizes HY5 activity to repress photomorphogenesis in Arabidopsis. New Phytol.229, 3221–3236. doi: 10.1111/nph.17114
16
KalayuG. (2019). Phosphate solubilizing microorganisms: promising approach as biofertilizers. Int. J. Agron.2019, 4917256. doi: 10.1155/2019/4917256
17
KimS.-A.LeeD. H.RyuB. H.HanN. S. (2023). Strain-specific barcode PCR and quantitative PCR assay for identification and enumeration of kimchi starter, Leuconostoc mesenteroides DRC1506. LWT Food Sci. Technol.174, 114460. doi: 10.1016/j.lwt.2023.114460
18
LiX.CooperN. G. F.O'TooleT. E.RouchkaE. C. (2020). Choice of library size normalization and statistical methods for differential gene expression analysis in balanced two-group comparisons for RNA-seq studies. BMC Genomics21, 75. doi: 10.1186/s12864-020-6502-7
19
LiY.ShaoJ.XieY.JiaL.FuY.XuZ.et al. (2021). Volatile compounds from beneficial rhizobacteria Bacillus spp. promote periodic lateral root development in Arabidopsis. Plant Cell Environ.44, 1663–1678. doi: 10.1111/pce.14021
20
LiY.WuF.LiC. (2025). Jasmonate signaling: integrating stress responses with developmental regulation in plants. J. Genet. Genomics52, 1490–1506. doi: 10.1016/j.jgg.2025.07.007
21
LidburyI. D. E. A.MurphyA. R. J.ScanlanD. J.BendingG. D.JonesA. M. E.MooreJ. D.et al. (2016). Comparative genomic, proteomic and exoproteomic analyses of three Pseudomonas strains reveals novel insights into the phosphorus scavenging capabilities of soil bacteria. Environ. Microbiol.18, 3535–3549. doi: 10.1111/1462-2920.13390
22
LvZ.HaoL.MaB.HeZ.LuoY.XinY.et al. (2021). Ciboria carunculoides suppresses mulberry immune responses through regulation of salicylic acid signaling. Front. Plant Sci.12, 658590. doi: 10.3389/fpls.2021.658590
23
MankotiaS.JakharP.SatbhaiS. B. (2024). HY5: a key regulator for light-mediated nutrient uptake and utilization by plants. New Phytol.241, 1929–1935. doi: 10.1111/nph.19516
24
MiljakovićD.MarinkovićJ.Balešević-TubićS. (2020). The significance of Bacillus spp. in disease suppression and growth promotion of field and vegetable crops. Microorganisms8, 1037. doi: 10.3390/microorganisms8071037
25
ModiA.VaiS.CaramelliD.LariM. (2021). The Illumina sequencing protocol and the NovaSeq 6000 system. Methods Mol. Biol.2242, 15–42. doi: 10.1007/978-1-0716-1099-2_2
26
OleńskaE.MałekW.WójcikM.SwiecickaI.ThijsS.VangronsveldJ. (2020). Beneficial features of plant growth-promoting rhizobacteria for improving plant growth and health in challenging conditions: a methodical review. Sci. Total Environ.743, 140682. doi: 10.1016/j.scitotenv.2020.140682
27
OubohssaineM.SbabouL.AuragJ. (2025). Enhancing ecosystem restoration and soil productivity through PGPR: a sustainable approach to bioremediation and biofertilization. Discov. Appl. Sci.7, 1379. doi: 10.1007/s42452-025-07510-3
28
PanG.YangX.HeJ.LiuZ.ChenF.ChenJ. (2024). Comprehensive analyses of the ARF gene family in cannabis reveals their potential roles in regulating cannabidiol biosynthesis and male flower development. Front. Plant Sci.15, 1394337. doi: 10.3389/fpls.2024.1394337
29
PataniA.PrajapatiD.AliD.KalasariyaH.YadavV. K.TankJ.et al. (2023). Evaluation of the growth-inducing efficacy of various Bacillus species on the salt-stressed tomato (Solanum lycopersicum). Front. Plant Sci.14, 1168155. doi: 10.3389/fpls.2023.1168155
30
SantoyoG.Moreno-HagelsiebG.del Carmen Orozco-MosquedaM.GlickB. R. (2016). Plant growth-promoting bacterial endophytes. Microbiol. Res.183, 92–99. doi: 10.1016/j.micres.2015.11.008
31
ShaoJ.LiS.ZhangN.CuiX.ZhouX.ZhangG.et al. (2015). Analysis and cloning of the synthetic pathway of the phytohormone indole-3-acetic acid in the plant-beneficial Bacillus amyloliquefaciens SQR9. Microb. Cell Fact.14, 130. doi: 10.1186/s12934-015-0323-4
32
SinghK. S.van der HooftJ. J. J.van WeesS. C. M.MedemaM. H. (2022). Integrative omics approaches for biosynthetic pathway discovery in plants. Nat. Prod. Rep.39, 1876–1896. doi: 10.1039/d2np00032f
33
SunQ.ZuoA.ZhangT. (2017). Advances in studies on chemical constituents in plants of Cimicifuga L. and their bioactivities and clinical application. Chin. Tradit. Herb. Drugs48, 3005–3016. doi: 10.7501/j.issn.0253-2670.2017.14.030
34
SurachatK.TaylorT. D.WattanamatiphotW.SukpisitS.JeenkeawpiamK. (2022). aTAP: automated transcriptome analysis platform for processing RNA-seq data by de novo assembly. Heliyon8, e10255. doi: 10.1016/j.heliyon.2022.e10255
35
TianD.SongX.LiC.ZhouW.QinL.WeiL.et al. (2021). Antifungal mechanism of Bacillus amyloliquefaciens strain GKT04 against Fusarium wilt revealed using genomic and transcriptomic analyses. Microbiol. Open10, e1192. doi: 10.1002/mbo3.1192
36
TimofeevaA.GalyamovaM.SedykhS. (2022). Prospects for using phosphate-solubilizing microorganisms as natural fertilizers in agriculture. Plants11, 2119. doi: 10.3390/plants11162119
37
UpadhyayS. K.SrivastavaA. K.RajputV. D.ChauhanP. K.BhojiyaA. A.JainD.et al. (2022). Root exudates: mechanistic insight of plant growth promoting rhizobacteria for sustainable crop production. Front. Microbiol.13, 916488. doi: 10.3389/fmicb.2022.916488
38
VersmessenN.Van SimaeyL.NegashA. A.VandekerckhoveM.HulpiauP.VaneechoutteM.et al. (2024). Comparison of DeNovix, NanoDrop and Qubit for DNA quantification and impurity detection of bacterial DNA extracts. PloS One19, e0305650. doi: 10.1371/journal.pone.0305650
39
WahabA.BibiH.BatoolF.MuhammadM.UllahS.ZamanW.et al. (2024). Plant growth-promoting rhizobacteria biochemical pathways and their environmental impact: a review of sustainable farming practices. Plant Growth Regul.104, 637–662. doi: 10.1007/s10725-024-01218-x
40
WangK.YinJ.ChenJ.MaJ.SiH.XiaD. (2024). Inhibition of inflammation by berberine: molecular mechanism and network pharmacology analysis. Phytomedicine128, 155258. doi: 10.1016/j.phymed.2024.155258
41
WangW.PaikI.KimJ.HouX.SungS.HuqE. (2021). Direct phosphorylation of HY5 by SPA1 kinase to regulate photomorphogenesis in Arabidopsis. The New phytologist. 230 (6), 2311–26. doi: 10.1111/nph.17332
42
XuL.SunX.WanX.LiH.YanF.HanR.et al. (2020). Identification of a Bacillus amyloliquefaciens H6 thioesterase involved in zearalenone detoxification by transcriptomic analysis. J. Agric. Food. Chem.68, 10071–10080. doi: 10.1021/acs.jafc.0c03954
43
XuX.HuJ.YuanZ. (2024). Stabilization or degradation? Post-translational modifications of JAZ proteins in plants. Mol. Plant17, 1002–1004. doi: 10.1016/j.molp.2024.06.010
44
YadavA.SinghD.LingwanM.YadukrishnanP.MasakapalliS. K.DattaS. (2020). Light signaling and UV-B-mediated plant growth regulation. J. Integr. Plant Biol.62, 1270–1292. doi: 10.1111/jipb.12932
45
YangL.YangY.HuangL.CuiX.LiuY. (2022). From single- to multi-omics: future research trends in medicinal plants. Brief. Bioinform.24, bbac485. doi: 10.1093/bib/bbac485
46
ZhangY.IshikawaM.KoshioS.YokoyamaS.DossouS.WangW.et al. (2024a). Effects of dietary supplementation with Bacillus subtilis natto on growth, digestive enzyme activity, immune response, and intestinal microorganisms of red sea bream, Pagrus major. Fishes9, 446. doi: 10.3390/fishes9110446
47
ZhangY.XieS.ZhangY.TangX.TangZ.LiuY.et al. (2026). A study on the metabolic characteristics of three rhizosphere bacteria promoting the growth of Cimicifuga dahurica Maxim. Ind. Crops Prod.240, 122622. doi: 10.1016/j.indcrop.2025.122622
48
ZhangQ.WeiW.JinX.LuJ.ChenS. J.OgajiO. D.et al. (2024b). Traditional uses, phytochemistry, pharmacology, quality control and clinical studies of Cimicifugae Rhizoma: a comprehensive review. Chin. Med.19, 1–18. doi: 10.1186/s13020-024-00937-7
49
ZhaoH.WangX.WangW. (2021). Bacillus amyloliquefaciens SN16-1-induced resistance system of the tomato against Rhizoctonia solani. Pathogens11, 35. doi: 10.3390/pathogens11010035
Summary
Keywords
Bacillus amyloliquefaciens ART9, plant hormone signal transduction, Radix Serratulae Chinensis, rhizosphere microbiome, transcriptomics
Citation
Wang J, Xu J, Fu P, Li J, Tang B, Zhang Y, Zeng Z and Fu Y (2026) Bacillus amyloliquefaciens ART9 modulates the growth of Radix serratulae Chinensis through rhizosphere microbiome reprogramming and root transcriptomic regulation. Front. Plant Sci. 17:1891133. doi: 10.3389/fpls.2026.1891133
Received
26 May 2026
Revised
22 June 2026
Accepted
24 June 2026
Published
10 July 2026
Volume
17 - 2026
Edited by
Yuri Shavrukov, Flinders University, Australia
Reviewed by
Chuansheng Mei, Institute for Advanced Learning and Research, United States
Mark Owusu Adjei, Nanjing Forestry University, China
Updates
Copyright
© 2026 Wang, Xu, Fu, Li, Tang, Zhang, Zeng and Fu.
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: Zhixiong Zeng, zengzx@126.com; Yafen Fu, fyf12152026@163.com
†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.