- 1Key Laboratory for Research and Development of “Qin Medicine” of Shaanxi Administration of Traditional Chinese Medicine, Shaanxi University of Chinese Medicine, Xi’an, China
- 2College of Pharmacy and Shaanxi Qinling Application Development and Engineering Center of Chinese Herbal Medicine, Shaanxi University of Chinese Medicine, Xi’an, China
- 3State Key Laboratory of Research and Development of Characteristic Qin Medicine Resources (Cultivation), Shaanxi University of Chinese Medicine, Xianyang, China
Rheum officinale Baill., a medicinal herb rich in anthraquinones and tannins, exhibits region-specific variation in bioactive compound accumulation. Given the well-documented pharmacological properties of R. officinale Baill. including its laxative, antibacterial, anti-inflammatory and regulating intestinal function, this study integrated Illumina sequencing of rhizosphere microbiomes, HPLC quantification of 10 active components (e.g., rhein, physcion), and soil analysis across three cultivation regions Shaanxi (ZB), Hubei (HB), Chongqing (CQ). Results demonstrated significant regional disparities: ZB showed highest rhein and catechin levels, while CQ accumulated more physcion. Microbial diversity followed the order ZB>HB>CQ, with Proteobacteria and Ascomycota dominating bacterial and fungal communities, respectively. Soil pH, moisture (SWC), and Zn/Cu content strongly correlated with microbial structure. Notably, Rokubacteriales was significantly positively associated with anthraquinones accumulation. These findings suggest that soil properties modulate microbial communities, which in turn regulate secondary metabolite biosynthesis through nutrient cycling (e.g., nitrogen/phosphorus metabolism). This study elucidates the tripartite interaction of soil-microbe-metabolite networks in R. officinale Baill., providing insights for geoherbalism optimization. Future research will focus on optimizing R. officinale Baill. cultivation through soil condition management to enhance both quality and yield.
Introduction
Rheum officinale Baill. as a tall perennial herbaceous plant, its dried roots and rhizomes are known as the Chinese herbal medicine rhubarb (Chinese Pharmacopoeia Commission, 2025). The main chemical constituents include anthraquinones, anthrones, stilbenes, phenylbutanones, polysaccharides, tannins, and volatile oils (Cao et al., 2017). The laxative activity of R. officinale Baill. is primarily attributed to its anthraquinones derivatives (rhein, aloe-emodin) and anthrone glycosides (sennoside) (Duval et al., 2016; Takayama et al., 2012; Feng et al., 2013). The growing medicinal demand for R. officinale Baill. has outpaced wild supply, necessitating large-scale cultivation. While intensive farming system are now established in China’s major production regions, persistent monoculture challenges have emerged. Notably, soil-borne diseases (particularly root rot) have increased substantially due to continuous cropping issues, threatening industry sustainability (Lin et al., 2024; Xu et al., 2025). Research confirms these agricultural problems correlate strongly with soil microbiome alterations, where pathogenic bacterial/fungi and nutrient imbalances directly impact plant health (Pang et al., 2022; Bi et al., 2023). Beneficial secondary metabolites (such as medicinal plant active ingredients, microbial functional compounds, plant polyphenols, etc.) have exploded in market demand in recent years due to their high application value in medicine, food and health care, etc.; The contradiction between insufficient adaptation of production scale and demand expansion is continuing to give birth to deepening basic research in this field.
Rhizosphere microorganisms are often termed the plant’s ‘second genome’. In natural environment, plant roots harbor a diverse community of microorganisms within and surrounding the root system, collectively known as the rhizosphere microbiome, which includes bacteria and fungi (Luo et al., 2024). This microbiome accompanies plants throughout their life cycle, interacts closely with them, and plays a vital role in the growth, development, and secondary metabolite accumulation of medicinal plants (Yuan et al., 2022; Peng et al., 2020). Among these microorganisms are plant growth-promoting rhizobacteria (PGPR). PGPR constitute a core functional group that directly or indirectly sustain plant growth and soil health through mechanisms such as promoting nutrient uptake, secreting plant hormones, enhancing tolerance to abiotic stresses, and suppressing pathogenic microorganisms (Jabborova et al., 2024; Vaghela and Gohel, 2023). Certain rhizobacteria directly promote plant growth by secreting plant hormones or influencing the plant’s own hormone synthesis. Among them, phosphate-solubilizing bacteria (Bacillus, Pseudomonas, Enterobacter and Escherichia) represent the most dominant phosphorus-transforming groups in soil (Wang et al., 2023; Pang et al., 2024). Nitrogen-fixing bacteria such as Klebsiella and Rhodococcus can convert atmospheric nitrogen into ammonia, increasing soil nitrogen content and thereby promoting the growth of medicinal plants and the accumulation of their active ingredients (Shi et al., 2023). For instance, the latest research has found that Paraburkholderia, a probiotic flora in Coptis chinensis, can promote plant growth, activate immune responses, inhibit the main pathogens of coptis root rot, and significantly improve the overall health of plants (Cao et al., 2025), while Streptomyces TM32 from the rhizosphere of Curcuma longa exhibits potent antimicrobial effects against plant pathogens (Nakaew et al., 2019). Additionally, Burkholderia from the rhizosphere of Baphicacanthus cusia, as well as Pseudomonas and Pantoea strains from the rhizosphere of Salvia miltiorrhiza, can also directly participate in indigo biosynthesis and phenolic compound accumulation (Zeng et al., 2018; You et al., 2017).
Rhizosphere fungi are another important group within rhizosphere microorganisms, and their diversity plays a crucial role in stabilizing soil ecosystems, promoting plant growth, and directly participating in nutrient cycling between soil and plants, thereby making significant contributions to soil ecology (Adedayo and Babalola, 2013). Functionally, they are categorized into saprophytic fungi and symbiotic mycorrhizal fungi: beneficial fungi can enhance plant root activity, promote growth, and improve disease resistance, whereas harmful fungi may induce plant diseases or even lead to plant mortality (Finzi Razavi et al., 2016; Li Destri Nicosia et al., 2015). Among them, arbuscular mycorrhizal (AM) fungi form symbiotic associations with most terrestrial plants. By assisting hosts in absorbing water and mineral nutrients such as phosphorus (while acquiring organic nutrients from the hosts), they significantly influence plant secondary metabolism (Jiang et al., 2017; Welling et al., 2016; Kapoor et al., 2017). For instance, in Salvia miltiorrhiza, synthetic communities composed of fungi can significantly promote plant growth and improve the quality and yield of medicinal materials (Jia et al., 2024). Moreover, AM fungi are closely associated with the synthesis of phenolic compounds (Rosa-Mera et al., 2011). In addition, Aspergillus flavus can convert soluble arsenic into insoluble forms, reducing its toxicity to soil organisms to ensure plant growth (Mohd et al., 2019). In summary, the community diversity and functional specificity of rhizosphere microorganisms are the core research targets for elucidating the mechanisms underlying the growth regulation of medicinal plants and the accumulation of their active ingredients. Related research in this field has provided an important theoretical basis for the optimization of medicinal plant cultivation and the utilization of microbial resources.
This study investigated the relationship between rhizosphere microorganisms and R. officinale Baill. Using the second-generation NovaSeq platform for PE250 amplicon sequencing, we characterized bacterial and fungal communities in the 3-year-old plants. Our analysis focused on two key aspects: how these rhizosphere microorganisms influence soil physicochemical properties, and their impact on the accumulation of active components in medicinal plant R. officinale Baill.
Materials and methods
Sample collection
Field sample was conducted in November 2023 across three major Rheum officinale producing regions: Shaanxi, Hubei and Chongqing. At each location, we collected three-year-old plants along with their rhizosphere soil (Li et al., 2021). Geographic coordinates for all sampling sites (longitude and latitude information) were recorded using Global Positioning System (GPS) (Supplementary Table S1). The study employed a replicated sampling design with three sampling points per region, each comprising three biological R. officinale Baill. replicates, yielding a total of nine sample batches. Notably, Zhenba County in Shaanxi Province served as a key sampling site due to its provincial recognition as a genuine producing area. All samples were collected from established medicinal cultivation bases. Three-year-old plants showing no visible signs of damage or disease were carefully selected (Figure 1, Supplementary Figure S1).
For soil sampling, the soil adhering to the roots of R. officinal was collected by gentle shaking, and rhizosphere soil samples tightly bound to the fibrous roots were carefully collected using sterilized brushes (Supplementary Figure S2). The frozen soil samples were dispatched to Shanghai Paisno Biotechnology Co., Ltd. for the detection of microbial diversity. Meanwhile, the remaining soil samples were air-dried to remove impurities, passed through a 20-mesh sieve, and then subjected to the determination of soil physical and chemical properties.
Determination of 10 secondary metabolites in R. officinale Baill.
The high-performance liquid chromatography (HPLC) analysis was performed according to an established method for R. officinale Baill (Li et al., 2025), enabling simultaneous quantification of ten bioactive components: gcatechin, catechin, sennoside B, chrysophanol-8-O-glucoside, emodin-8-O-glucoside, aloe-emodin, rhein, emodin, chrysophanol and physcion. The fresh R. officinale Baill. herbs underwent the following steps. The workflow is summarized in the schematic below (Jiang et al., 2025), illustrating the key steps from sample preparation to chromatographic analysis (Supplementary Figure S3). First, the coarse outer skin was carefully scraped off. The peeled R. officinale Baill. was then sliced into 3.0 cm pieces. These slices were air-dried for seven consecutive days and subsequently dried in an electric blast drying oven at 50°C for seven days. Finally, the dried R. officinale Baill. slices were pulverized and sieved through a No. 4 sieve (with an inner diameter of sieve hole of 250 µm) to acquire samples appropriate for liquid phase detection. Precisely, 100 mg of each powdered sample was transferred to a 50 mL stoppered flask, mixed with 4.5 mL methanol, and weighed. After 30 minutes of ultrasonication (500W, 40kHz), samples were returned to room temperature, reweighted to compensate for solvent lose, then centrifuged at 11260 xg for 10 minutes. The supernatant was filtered through a 0.22 μm membrane. The reference substances were purchased from Chengdu Desite Biotechnology Co., Ltd. Individual reference standards were accurately weighed and dissolved in methanol to prepare10 mL stock solutions. A mixed working standard solution was prepared by combining 1.0 mL aliquots of each primary standard, followed by ten-fold dilution with methanol. All standard solutions were stored at 4 °C. Separation was achieved using a Maple Unitary C18 column (5 μm, 4.6 × 250 mm) with a binary mobile phase system: (A) methanol and (B) 0.2% phosphoric acid in water. The analysis was conducted at 30 °C with a flow rate of 1.0 mL/min, detection wavelength of 260 nm, and injection volume of 10 μL. Gradient elution was performed according to program Supplementary Table S2.
Determination of soil physicochemical properties
Soil pH was measured with a Sardolus PB-10 pH meter. Total Nitrogen (TN) content was determined through sulfuric acid-catalyzed digestion. Flame photometric was used to quantify total Potassium (TK) content. For phosphorus analysis, total Phosphorus (TP) and available Phosphorus (AP) were determined using the NaOH fusion method followed by molybdenum-antimony anti-colorimetric method. Soil organic carbon and organic matter were determined via the potassium dichromate-concentrated sulfuric acid external heating method. Ammonium nitrogen (NH4+-N) and Nitrate nitrogen (NO3-N) was analyzed utilizing the potassium chloride extraction method. Soil Moisture Content (SWC) was determined gravimetrically by oven drying. Micronutrients analysis such as iron (Fe), manganese (Mn), zinc (Zn) and copper (Cu) was performed using acid extraction followed by atomic absorption spectrophotometry (Bao, 2000). All soil physical and chemical properties analyses were conducted by Yangling Xinhua Ecological Technology Co., Ltd.
Soil DNA extraction, PCR amplification, Illumina sequencing, and data processing
Total soil DNA was extracted using the OMEGA Soil DNA Kit (D5635-02) (OMEGA Bio-tek, USA), with extraction quality verified by 0.8% agarose gel electrophoresis. For bacterial community analysis, the V3-V4 hypervariable region of 16S rRNA gene was amplified using universal primers 338F: (5’-ACTCCTACGGGAGGCAGCAG-3’) and 806R: (5’-GGACTACHVGGGTWTCTAAT-3’), with diluted DNA as template. Fungal community analysis targeted the ITS1 region using primers ITS1F(5’-GGAAGTAAAAGTCGTAACAAGG-3’) and ITS2R(5’-GCTGCGTTCTTCATCGATGC-3’) (Shanghai Paisano Biotechnology Co., Ltd., Shanghai). PCR products were purified using the Axygen Gel Recovery Kit (OMEGA, USA), followed by library preparation with the TruSeq Nano DNA LT Library Prep Kit (Illumina). Paired-end (2 × 250 bp) sequencing was performed on the Illumina MiSeq platform (Illumina, San Diego, CA) at Shanghai Personalbio Technology Co., Ltd.
Statistical analysis of data
The paired-end sequencing was conducted using the NovaSeq 6000 SP Reagent Kit (500 cycles) (Illumina Shanghai Personalbio Technology Co., Ltd., Shanghai). Raw sequence data were demultiplexed using the demux plugin, followed by primer removal with cutadapt plugin. Subsequent processing included quality filtering, denoising, read merged and chimera removal using DADA2 plugin to generate amplicon sequence variation (ASVs). The ASVs were using MAFFT, and phylogenetic trees with constructed with FastTree2. Bioinformatic analyses were primarily performed in QIIME2 and R (Bokulich et al., 2018). Alpha diversity indices (Chao1 index and Shannon index) were calculated from the ASV table and visualized as boxplots. Microbial taxonomic composition and relative abundance were analyzed using MEGAN and GraPhlAn (Asnicar et al., 2015; Huson et al., 2011). Beta diversity was assessed through principal coordinates analysis (PCoA) based on the Bray-curtis dissimilarity matrices, with community composition differences visualized in two-dimensional sorting plots. Statistical analyses were conducted by using SPSS26.0 for one-way ANOVA, Kruskal-Wallis and Duncan’s multiple comparison test (α=0.05). GraphPad Prism 8 software was used for plotting.
Results
The content of bioactive ingredients in R. officinale Baill. exhibits interregional variations
To investigate the phytochemical variability of medicinal R. officinale Baill. from different geographical origins, we quantitatively analyzed ten characteristic bioactive ingredients using HPLC (Supplementary Figure S4). Comparative analysis revealed significant compositional variations in six components including rhein (Figure 2A), chrysophanol-8-O-glucoside (Figure 2B), catechin (Figure 2C), emodin-8-O-glucoside (Figure 2D), gallic acid (Figure 2E) and physcion (Figure 2F), whereas the remaining four components exhibited relatively Table Sconcentrations without statistically significant difference (Supplementary Figure S5). Specifically, R. officinale Baill. from ZB exhibited significantly higher levels of rhein (3.50 × CQ, 1.85 × HB), chrysophanol-8-O-glucoside (2.89 × CQ, 1.54 × HB), catechin (3.30 × CQ, 1.48 × HB) and emodin-8-O-glucoside (1.36 × CQ, 1.61 × HB) but lower gallic acid content compared to other regions. Conversely, CQ showed the highest physcion content (1.63 × ZB, 1.56 × HB). No significant interregional differences were observed for sennoside B, emodin, aloe-emodin, or chrysophanol. These findings demonstrate distinct chemotypic profiles associated with geographical origins, which may influence the pharmacological properties of R. officinale Baill. products.

Figure 2. (A-F) is the numbering of the active ingredients in R. officinale Baill. The content of 6 effective components in R. officinale Baill. from 3 origins. The abscissa is the sample grouping, and the ordinate is the content of the components. Data are the mean of three replicates ± SE (standard error); different letters indicate significant differences at p < 0.05 according to analysis of variance (ANOVA).
Analysis of the composition and diversity of fungal and bacterial colonies in rhizosphere soil
To investigate the potential correlation between rhizosphere microbial communities and phytochemical variations in R. officinale Baill., we performed comprehensive analysis of bacterial and fungal populations in the rhizosphere soils from different cultivation regions. Microbial community analysis was performed on nine rhizosphere soil samples collected from three geographic regions using Illumina HiSeq high-throughput sequencing. The rarefaction curves of ITS (Supplementary Figure S6A) and 16S rRNA (Supplementary Figure S6B) sequencing demonstrated the continuous emergence of novel operational taxonomic units (OTUs) with increasing sequencing effort, indicating previously uncharacterized microbial diversity. Curve plateauing at sufficient sequencing depth, suggesting adequate coverage of extant microbial diversity with samples (Good’s coverage >99%). For fungi, the sequencing depth ranges from 0 to 60,000 reads, and the effective sequencing depth (saturation depth) is 10,000–60,000 reads. For bacteria, the sequencing depth ranges from 0 to 25,000 reads, with a saturation depth of 5,000–10,000 reads. The asymptotic behavior of the curves confirms that the sequencing depth achieved was sufficient to capture the majority of bacterial and fungal diversity present in the samples, while further sequencing would primarily reveal rare taxa.
Microbial community analysis was conducted at 97% sequence similarity threshold to determine operational OTUs in rhizosphere soil of R. officinale Baill. from three regions. Venn diagram analysis revealed distinct distribution patterns of shared and unique OTUs. For the fungal community (Figure 3A), there were a total of 1891 OTUs in the three regions, and only 52 OTUs were shared among all three regions, ZB region exhibited the highest diversity with 842 total OTUs, accounting for 44.53% of the total. While, HB region contained 643 total OTUs, accounting for 34.00%. And CQ region showed 616 total OTUs, accounting for 26.97%. For the bacterial community, a total of 13,360 OTUs were present in the three regions, of which 126 OTUs were common to all regions (Figure 3B). The samples from ZB demonstrated maximum diversity, with a total of 6364 OTUs, accounting for 47.63% of the total. Whereas, samples form HB contained 4625 OTUs, accounting for 34.62%, and samples from CQ showed only 2623 OTUs, accounting for 13.63%. The overall microbial diversity followed a consistent regional pattern: ZB>HB>CQ (p<0.05), indicating significantly higher microbial diversity in ZB rhizosphere soils compared to other regions.

Figure 3. Venn diagrams depicting fungi (A) and bacteria (B) in the rhizosphere soil of R. officinale Baill. across three regions. PCoA analysis of fungi (C) and bacteria (D) in rhizosphere soil of R. officinale Baill. among three regions. α diversity analysis of bacteria (E) and fungi (F) in the rhizosphere soil of R. officinale Baill. in three regions. The Venn diagram shows the number of OTUs common and unique to the three regions, with different colors representing different regions. Each point in the PCoA diagram represents a sample, and the dots of different colors indicate different groupings. The percentage in the coordinate axis bracket represents the proportion of the sample difference data (distance matrix) that the corresponding coordinate axis can interpret. Box plots visually compare diversity between groups. α diversity index reflects the species richness and uniformity in the sample. Higher index values indicate higher microbial community complexity.
Alpha diversity analysis was performed to assess microbial community characteristics across different regions. The Chao1 index, which estimates species richness, showed higher values corresponding to greater community richness. The Shannon index, incorporating both species richness and evenness, indicated higher community diversity with increasing values. Both indices followed the same trend: ZB (fungi: Chao1 = 377.92, Shannon=6.01; bacteria: Chao1 = 3274.80, Shannon=10.55) > HB (fungi: Chao1 = 289.33, Shannon=4.36; bacteria: Chao1 = 2247.59, Shannon=10.10) > CQ (fungi: Chao1 = 262.11, Shannon=5.72; bacteria: Chao1 = 1419.09, Shannon=8.73). According to bacterial communities, significant regional variations were detected in both Chao1(P<0.05) and Shannon indices (P<0.05) (Figure 3C). Chao1 index ranking: ZB>CQ>HB. Shannon index pattern: CQ>HB>ZB. Regarding fungal communities, no statistically significant differences were detected in either Chao1 (P >0.05) or Shannon (P >0.05) (Figure 3D) indices among regions. These results demonstrate distinct spatial patterns in microbial diversity, with fungal communities showing consistent richness-diversity.
Principal Component Analysis (PCA) was employed to assess beta diversity patterns among sampling regions. For fungal communities (Figure 3E), the first two principal coordinates (PCo1 and PCo2) explained 42.03% of total variation (PCo1:22%; PCo2 20.3%). Bacterial communities (Figure 3F) showed stronger separation, with PCo1 (36.6%) and PCo2 (23.3%) collectively accounting for 59.9% of observed variation. This indicates that the selected principal components are effective in capturing the major patterns of difference in the raw data. In the PCoA diagram, the soil samples of ZB, HB and CQ were distributed in different quadrants and were far away from each other, indicating that there was significant spatial heterogeneity in the composition of fungal and bacterial communities in rhizosphere soil in the three regions.
Composition of fungal and bacterial communities
We further analyzed bacterial and fungal communities in R. officinale Baill. rhizosphere soils in Baill. three regions. Microbial abundance data were normalized during analysis, and species were matched in the NCBI database. Based on the horizontal abundance data of soil fungi in ZB, HB and CQ, the predominant fungal phylum across all regions was Ascomycota, accounting for 73.08%-92.14% of total fungal sequences. Secondary phyla included Basidiomycota (2.27%-15.30%), Rozellomycota (0.22%-1.52%) and Mortierellomycota (0.33%-3.79%), while other phyla each represented was less than 0.1% (Figure 4A). This consistent dominance pattern suggests Ascomycota’s ecological adaptation to the studied rhizosphere environments. However, in accordance with the bacterial community structure, the regional variations in the composition bacterial phyla were more pronounced. The horizontal abundance distribution of soil bacteria phyla in the three regions was as follows: Proteobacteria (30.20%) and Acidobacteria (20.30%) were the dominant phyla in ZB samples, followed by Actinobacteriota (11.71%), Chloroflexi (4.41%), Verrucomicrobiota (3.26%), Bacteroidota (2.97%) and Gemmatimonadota (6.22%). Proteobacteria (37.39%) and Actinobacteriota (18.11%) were the dominant categories in HB samples, followed by Acidobacteria (16.11%), Chloroflexi (5.72%), Gemmatimonadota (5.98%), Verrucomicrobiota (2.23%) and Bacteroidota (5.23%). Proteobacteria (26.77%) and Acidobacteriota (21.27%) were the dominant phyla in the CQ sample, followed by Gemmatimonadota (9.31%), Verrucomicrobiota (3.01%), Chloroflexi (19.69%), Bacteroidota (1.08%) and Myxococcota (1.74%) (Figure 4A). The microbial communities were predominantly composed of known beneficial taxa, including plant-growth promoting Proteobacteria and organic matter-degrading Acidobacteria in bacteria, along with Ascomycota in fungi. The differential abundance patterns suggest region-specific environmental selection pressures shaping microbial community assembly.

Figure 4. Relative abundance of fungi at the phylum (A) and genus (B) levels in the cultivated soil of R. officinale Baill. from three origins. F, fungi; B, bacteria. The picture shows the top 10 taxa by relative abundance at the phylum and genus level. The colors distinguish between microbial taxa, while the x-axis indicates the sample name and the y-axis indicates the relative abundance (%).
The top 10 soil fungal genus in ZB, HB and CQ were Penicillium, Humicola, Fusarium, Longitudinalis, Staphylotrichum, Varicosporellopsis, Solicoccozyma, Oidiodendron, Saitozyma and Trichoderma (Figure 4B). Among them, Penicillium (5.59%) was the most abundant in ZB, Humicola (41.13%) was the most abundant in HB, and Penicillium (10.75%) was the most abundant in CQ. Notably, Fusarum showed significant regional variation (CQ:6.40% > HB:2.86 > ZB: 2.11%). Eight of the den dominant fungal genera are known beneficial microorganisms, with Trichoderma and Penicillium being particularly notable for their plant-growth promoting properties. The top 10 species of soil bacteria were JG30-KF-AS9, SC-I-84, Sphingomonas, Subgroup_2, Rokubacteriales, Vicinamibacteraceae, Gemmatimonas, RB41, KD4–96 and Candidatus_Solibacter. Rokubacteriales (4.84%) and SC-I-84 (4.67%) were the dominant genera in ZB, Sphingomonas (6.87%) and SC-I-84 (6.55%) were dominant in HB, and JG30-KF-AS9 (14.56%) and Subgroup_2 (7.29%) were dominant in CQ (Figure 4B). Ecologically significant beneficial genera including Rokubacteriales, Sphingomonas and Gemmatimonas showed higher relative abundance in ZB and HB (8.24%-8.89%) compared to CQ (2.72%). This distribution pattern potential functional differences in microbial-mediated soil processes across regions.
Correspondence analysis on soil microorganisms from three different R. officinale Baill. cultivation regions
The soil samples from the three regions were roughly divided into three quadrants, indicating that there were some differences in the composition of fungal and bacterial communities among different soil samples. Among the fungal communities, the ZB region was dominated by Penicillium and Humicola (Figure 5A). Penicillium and Humicola were significantly enriched in the HB region. In the CQ area, Solicoccozyma and Varicosporellopsis are dominant. Among the physicochemical properties, the arrow of Cu was the longest and pointed to the negative direction of the RDA1 axis, indicating that it had the greatest impact on the community. The direction was consistent with the distribution of CQ samples, Solicoccozyma and Varicosporellopsis, indicating that Cu significantly affected the enrichment of CQ samples and Solicoccozyma and Varicosporellopsis. The WS pointed to the positive direction of the RDA2 axis, which was close to the CQ sample, which may have a certain effect on the fungal community in the CQ area. Zn points to the positive direction of the RDA1 axis, which is correlated with the distribution trend of the ZB group. In the bacterial community, the Cu arrow was the longest and pointed to the negative direction of the RDA1 axis, which was the most influential factor (Figure 5B). The direction was highly consistent with the distribution of CQ samples and JG30-KF-AS9, indicating that Cu may significantly affect the enrichment of CQ samples and JG30-KF-AS9. The arrows of Zn and WS are shorter and have a weaker influence. Zn points to the positive direction of the RDA1 axis, which is correlated with the distribution trend of samples in the ZB group. Taken together, the influencing factors of fungal and bacterial communities were about the same.

Figure 5. Redundancy analysis (RDA) of soil fungi (A) and bacteria (B) communities of R. officinale Baill. cultivated in three different locations and their environmental factors. The percentage of coordinate axes in the graph represents the degree of difference in the raw data it can interpret, with dots representing different samples and different colors representing different groupings. Arrows indicate environmental factors, and the length represents the degree to which the corresponding environmental factors are related to microorganisms. The angle between the arrows indicates the correlation between the two, with the acute angle positively correlated and the obtuse angle negatively correlated.
Correlation analysis of soil microorganisms and bioactive components of R. officinale Baill.
Correlation analysis between soil factors and the active ingredient content in R. officinale Baill. was performed using SPSS software (version 27). It was performed between the top ten most abundant soil microorganisms (bacteria and fungi) and ten active ingredients in R. officinale Baill (Tohge, 2020). Firstly, the normality test of the data is carried out, and the Spearman correlation analysis method is selected to test the initial correlation coefficient and p-value according to the results. In order to control the risk of false positives caused by multiple comparisons, FDR correction was further used to process the p-value to more reliably identify the key soil factors affecting the formation of active ingredients in wolfberry. Fungal correlation showed that there were no fungi significantly associated with the ten active ingredients among the top 10 fungi in abundance (p>0.05) (Figure 6A). In the composition of the bacterial community, there was a significant correlation between two bacteria and the components (Figure 6B). Rokubacteriales was significantly positively correlated with rhein and emodin-8-O-glucoside (p<0.05), and JG30-KF-AS9 was significantly negatively correlated with rhein, catechin, and chrysophanol-8-O-glucoside (p<0.01), and positively correlated with physcion (p<0.05).

Figure 6. Correlation analysis between the content of active ingredients in R. officinale Baill. and the levels of fungi (A) and bacteria (B) in the cultivation soil. The abscissa in the figure represents the microbial species, and the ordinate is the ten active ingredients. The color of the square represents the level of the correlation coefficient.
Differences in physicochemical properties of rhizosphere soil of R. officinale Baill. in different production areas
To explore whether the differences in soil are related to the variations in the components of R. officinale Baill., we conducted test on various physical and chemical indicators of the soil.
Significant regional variations in soil physicochemical properties were observed among the three study areas (Table 1). Soil pH analysis revealed neutral conditions in ZB (pH 6.5-7.5), while HB and CQ exhibited acidic soils (pH 4.5-5.5). Trace element analysis demonstrated distinct distributions patterns: iron (Fe) concentrations were markedly elevated in HB (6.78 × ZB, 1.38 × CQ; p<0.01), whereas ZB showed the highest manganese (Mn) levels (5.91× HB, 1.20× CQ; p<0.01). CQ soils contained significantly greater Copper (Cu) (2.01 × ZB, 1.24 × HB; p<0.05), while Zinc (Zn) showed no significant interregional differences. Nutrient analysis revealed HB’s superior fertility, with significantly higher concentrations of total nitrogen (TN: 1.65 × ZB, 1.95 × CQ), organic matter (OM: 2.04 × ZB, 2.48 × CQ), available phosphorus (AP: 7.74 × ZB, 8.46 × CQ), and total phosphorus (TP: 2.57 × ZB, 1.92 × CQ; all p<0.01). CQ exhibited the highest total potassium (TK: 1.45 × ZB, 1.12 × HB; p<0.05) and ammonium nitrogen (NH4+-N: 6.48 × ZB, 1.53 × HB; p<0.01). ZB showed predominant nitrate nitrogen (NO3--N: 2.93 × HB, 8.01 × CQ; p<0.01), while water-soluble content (SWC) was significantly lower in HB compared to other region.
Correlation analysis of soil physical and chemical properties and bioactive components of R. officinale Baill.
The correlation between soil properties and the active ingredients of R. officinale Baill. was further discussed, and the p-value was FDR corrected (Figure 7; Supplementary Table S3). The results showed that Zn, Mn, TN, OM, AP, TP and SWC were not significantly correlated with the 10 active ingredients. Specifically, pH was significantly positively correlated with gallic acid, catechin, chrysophanol-8-O-glucoside, and rhein (P <0.05); Fe and gallic acid, chrysophanol-8-O-glucoside and rhein were significantly positively correlated. Cu, NH4+-N and NO3-N were significantly positively correlated with catechin, chrysophanol-8-O-glucoside, emodin-8-O-glucoside and rhein. TK and catechin, chrysophanol-8-O-glucoside, emodin-8-O-glucoside and rhein were significantly positively correlated. AK was significantly positively correlated with physcion.

Figure 7. Correlation analysis between the content of active ingredients in R. officinale Baill. and soil properties. The lower left line represents the Mantel test results of soil physical and chemical properties and active ingredients. The thickness of the line represents the overall correlation coefficient, and the color of the line represents the significance test result of the correlation coefficient. The triangle is a heat map of correlations between substances, with cell color blocks indicating the size of the correlation coefficient, red indicating positive correlation, and blue indicating negative correlation.
Discussion
A growing body of research has shown that R. officinale Baill. has a variety of pharmacological properties, including anti-inflammatory, anticancer, antibacterial, and antiviral, and can also be used as an expectorant. It has been recognized as effective in treating a variety of conditions, including stomach pain, nausea, and vomiting (Keshavarzi et al., 2021). The medicinal properties of R. officinale Baill. stem from the synergistic effects of several active ingredients (Supplementary Table S4). Gallic acid and catechin belong to tannins, sennoside B belongs to dianthranone, chrysophanol-8-O-glucoside and emodin-8-O-glucoside belong to conjugated anthraquinone, and aloe-emodin, rhein, emodin, chrysophanol and physcion belong to free anthraquinone. According to the results of HPLC analysis, there exist distinct regional patterns in bioactive compound accumulation. Specifically, ZB region exhibited the highest concentrations of tannins (gallic acid, catechin), dianthrones (sennoside B), and conjugated anthraquinones (chrysophanol-8-O-glucoside, emodin-8-O-glucoside). Conversely, free anthraquinone (aloe-emodin, rhein, emodin, chrysophanol, physcion) were minimally abundant in ZB but reached peak levels in CQ region. Studies have revealed that the pharmacological effect of tannins predominantly manifest as astringency and antidiarrheal properties (Xu et al., 2024). Regarding anthone, its pharmacological effects mainly include laxative and promoting of bile secretion. When combined with anthraquinone, it exhibits a potent laxative effect. Moreover, it also possesses antioxidant functions. Free anthraquinone is associated with various effect such as antibacterial, anti-inflammatory, antitumor and antioxidant effects (Xiang et al., 2020). The pronounced inter-regional variability in active ingredients profiles underscore the significant impact of geographical factors on secondary metabolite biosynthesis in R. officinale Baill (Komatsu et al., 2006). These suggest that origin-specific quality variations may influence clinical efficacy.
Soil physicochemical properties serve as critical determinants of medicinal plant yield and quality through their influence on plant growth, nutrient uptake, secondary metabolism and microbial interaction (Lv G. et al., 2024). Through the association of soil physicochemical properties with active components, it was discovered that pH exhibited a strong correlation with the content of ten active components. In contrast, TN, AK, and heavy metal ions in soil were only correlated with specific individual components. Thus, we hypothesized that soil pH might be the core factor driving the quality variation of R. officinale Baill., consistent with observations in other medicinal species. Study revealed that soil stands as one of the crucial determinants of Panax notoginseng production, and soil physicochemical properties and soil microbiome have high contributions to the biomass and saponins of Panax notoginseng (Xia et al., 2024). Specific pH conditions proved beneficial for the synthesis and accumulation of THSG and total anthraquinones in Polygonum multiflorum, thereby enhancing the quality of the medicinal material (Wu et al., 2023). Enhanced medicinal quality of Astragalus mongholicus correlates with improved soil quality, particularly pH, soil organic matter (SOM), and calcium content (STCa). These factors directly regulate the biosynthesis of key bioactive compounds including asragalosides (I, II, and IV) and calycosin (Sun et al., 2020). Similar studies in Lycium barbarum found that soil pH mediates the accumulation of fruit flavonoids and total sugars, demonstrating both direct and indirect regulatory effects (Liu et al., 2022). It is worth noting that correlation analysis only provides preliminary insights into the potential association between soil physicochemical factors, active ingredients and microorganisms, but cannot establish a causal mechanism between the two, so the results need to be interpreted with caution.
Rhizosphere microorganisms can directly or indirectly regulate the accumulation process of plant active ingredients by building an “interaction network” with plant roots (Lv J. et al., 2024). In terms of nutrient uptake enhancement, nitrogen-fixing microorganisms can convert N2 in the air into ammonia nitrogen available to plants, and efficiently supplement nitrogen for plants. In addition, these microorganisms can synthesize plant hormones such as indoleacetic acid (IAA) to promote root growth, and cytokinins to enhance leaf photosynthesis and fruit development. At the same time, microorganisms can also regulate the physical and chemical properties of rhizosphere soil and create a suitable microenvironment for the accumulation of active ingredients (Hacquard et al., 2017; Zhong et al., 2022). Analysis of soil samples from different R. officinale Baill. cultivation regions revealed significant correlations between microbial community characteristics and specific environmental, with copper (Cu), zinc (Zn) and water content (WC) emerging as primary determinant. Heavy metal contamination, originating from mineral mining, fossil fuel combustion, and industrial/agricultural activities, has been documented to alter soil physicochemical properties and disrupt microbial composition and community structure through cadmium (Cd), arsenic (As), and Cu accumulation (Perveen et al., 2017). Experimental evidence demonstrates a concentration-dependent relationship between Cu exposure and microbial diversity reduction. Biolog analysis indicate an inverse correlation between CuCl2 concentration (0, 10, 20, 100 and 300 μmol/L) and Shannon diversity indices, confirming Cu’s suppressive effect on microbial diversity (Woźniak et al., 2019; Kong et al., 2006). In addition, it has been found that Cu has a particularly significant effect on the composition of rhizosphere bacterial communities, which may be due to the fact that Cu can not only destroy the cellular structure and function of microorganisms, but also change the spatial conformation of proteins by binding to specific proteins, produce harmful clumps, interfere with the functions of key proteins, lead to cellular metabolic disorders and signal transduction obstruction, enter a toxic stress state, and eventually die (An et al., 2023). In this study, Cu exhibited the strongest pronounced effects observed in the CQ region. This spatial pattern suggests localized Cu accumulation may drive significant microbiome alteration in CQ soils. In contrast, Zn and water content demonstrated comparatively weaker impacts on microbial communities. These results are only derived from high-throughput data analysis. The specific mechanisms still need to be further verified through experiments.
Fusarium, a prevalent pathogenic genus among the dominant fungi taxa, exhibits significant phytotoxicity toward medicinal plants. For example, in Panax ginseng, Fusarium can cause root rot, leading to root decay, which affects the ginseng’s absorption of water and nutrients, and in severe cases, the entire plant can die (Jun et al., 2022). When Salvia miltiorrhiza is infected by Fusarium, the content of effective components such as tanshinone will decrease significantly, thus affecting the medicinal value of Danshen (Yang et al., 2020). In medicinal plants like Atractylodes macrocephala, root rot and wilt diseases often occur, leading to poor growth, reduced yield, and diminished quality (Huang et al., 2021). Our data revealed a negative correlation between Fusarium and catechin levels, suggesting potential suppression of catechin biosynthesis, possibly through interference with phenylpropanoid pathway enzymes. This provides us with a clue that when Fusarium infects the plants, it may reduce the medicinal quality of R. officinale Baill. by directly or indirectly affecting the synthesis of catechin. Subsequently, molecular techniques can be combined with plant pathology methods for verification. As a beneficial fungus, Trichoderma can solubilize poorly soluble nutrients such as phosphorus and potassium in the soil by secreting organic acids, enzymes, and other substances. It converts these nutrients into forms that can be absorbed by plants, thereby enhancing the medicinal plants’ ability to absorb and utilize these nutrients and promoting plant growth and development. For instance, during the cultivation of Astragalus membranaceus, it facilitates plant growth by dissolving soil nutrients, promoting root elongation, and generating hormone-like substances (He et al., 2022). Correlation analysis revealed a modest positive association between Trichoderma and the concentration of physcion and aloe-emodin. This pattern suggests that Trichoderma might contribute to enhanced synthesis or accumulation of physcion and aloe-emodin. Conversely, Trichoderma exhibited significant negative correlations with emodin-8-O-glucoside, rhein, and chrysophanol-8-O-glucoside, indicating potential inhibitory effects on these metabolites. Regarding Penicillium, evidence suggests this fungal genus promotes growth and secondary metabolism in medicinal plants. For instance, mono-inoculation of the endophytic fungus Penicillium steckii in Tripterygium wilfordii significantly increased triptolide and tripterine production (Song et al., 2020). Based on this functional analogy, Penicillium may similarly facilitate aloe-emodin in R. officinale Baill. All correlation analyses were interpreted with careful consideration of inherent observational limitations. While these analyses reveal significant associations between microbial communities and metabolic profiles, they do not establish direct causation. To validate these findings and elucidate mechanistic relationships, we propose employing complementary experimental approaches, including microbial isolation of key operational taxonomic units, controlled co-cultivation systems integrating isolated microbes with plant tissues, and Table Sisotope-assisted metabolic tracing to quantify microbial contributions. This multi-approach will bridge the gap between correlative observations and causative inferences.
Notably, among the top ten bacteria phyla in relative abundance, many remain understudied regarding their functional roles. The bacteria taxa most strongly associated with bioactive components in R. officinale Baill. include Rokubacteriales, JG30-KF-AS9, RB41, Subgroup_2 and. RB41 belongs to the phylum Acidobacteria, which contributes significantly to soil carbon cycling. Acidobacteria generally account for approximately 20% of soil bacterial communities (Jones et al., 2009) and carry out ecologically beneficial functions including plant polymers degradation, iron cycling, and single-carbon compound metabolism (Pankratov et al., 2012; Lu et al., 2010). RB41 and Rokubacteriales are significantly positively correlated with catechin, chrysophanol-8-O-glucoside, emodin-8-O-glucoside, and rhein, indicating that bacteria from these two genera may promote the accumulation of these four components. Subgroup_2 and JG30-KF-AS9 are extremely significantly positively correlated with physcion, suggesting that bacteria from these two genera may promote the synthesis and accumulation of physcion, while being negatively correlated to varying degrees with emodin-8-O-glucoside, rhein, catechin, and chrysophanol-8-O-glucoside, possibly inhibiting the synthesis and accumulation of these four components. Due to the limited sample size (N = 3), the statistical power of our analyses may be insufficient to detect subtle effects or interactions. Correlation analysis only provides preliminary insights into the potential association between soil physicochemical factors, active ingredients and microorganisms, but cannot establish a causal mechanism between the two, so the results need to be interpreted with caution.
While these correlation analysis provide preliminary insights into potential microbe-metabolite relationships in R. officinale Baill., their inherent limitations warrant cautious interpretation. Correlative data alone cannot establish causal mechanisms and may reflect indirect associations influenced by unmeasured environmental variables or microbial interactions. To validate the proposed roles of RB41, Rokubacteriales, Subgroup_2, and JG30-KF-AS9 in metabolite modulation, targeted experimental approaches are recommended. For instance, in vitro co-culture assays with isolated bacterial strains and R. officinale Baill. tissues/cell lines to directly quantify effects on specific metabolite accumulation.
Data availability statement
The datasets presented in this study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov, accession numbers PRJNA1288235 and PRJNA1320834.
Author contributions
YW: Methodology, Formal analysis, Investigation, Writing – original draft. FY: Writing – review & editing, Software, Conceptualization. YJL: Investigation, Data curation, Methodology, Writing – review & editing, Conceptualization. XH: Writing – review & editing. JG: Visualization, Writing – review & editing, Project administration, Investigation. MZ: Writing – review & editing, Methodology, Investigation. YY: Conceptualization, Writing – review & editing, Supervision. GZ: Project administration, Writing – review & editing, Funding acquisition. YML: Funding acquisition, Visualization, Writing – review & editing, Project administration.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This work was supported by the National Natural Science Foundation of China (82104334, 81973430); Youth Innovation Team Research Project of Shaanxi Provincial Education Department (24JP043); Shaanxi Provincial Department of Education Service Local Special Plan Project (24JC030); The Key Research and Development Program of Xianyang (L2023-ZDYF-SF-017).
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.
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 5and 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.1650792/full#supplementary-material
Supplementary Figure 1 | Schematic diagram of the sample preparation and HPLC analysis procedure.
Supplementary Figure 2 | The content of four effective components in R. officinale Baill. sourced from three locations shows no obviously variations: sennoside B (A), emodin (B), aloe-emodin (C), chrysophano (D). Data are the mean of three replicates ± SE (standard error); different letters indicate significant differences at p < 0.05 according to analysis of variance (ANOVA). ZB for Zhenba, Shaanxi. CQ for Chongqing. HB for Hubei.
Supplementary Figure 3 | Rhubarb plants from three medicinal rhubarb cultivation bases. A: ZB; B: HB; C: CQ.
Supplementary Figure 4 | Soil sample for fungal ITS and bacterial 16S sequencing from three regions. ZB for Zhenba, Shaanxi. CQ for Chongqing. HB for Hubei.
Supplementary Figure 5 | HPLC chromatograms of R. officinale Baill. from three production areas.
Supplementary Figure 6 | Dilution curves of fungal ITS (A) and bacterial 16S (B) sequencing in soil samples. The abscissa is the flattening depth, and the ordinate is the alpha diversity index calculated 10 times.
References
Adedayo, A. A. and Babalola, O. O. (2013). Fungi that promote plant growth in the rhizosphere boost crop growth. J. Fungi (Basel) 9, 239. doi: 10.3390/jof9020239
An, Q. R., Li, Y. Y., Zheng, N., Ma, J., Hou, S. N., Sun, S. Y., et al. (2023). Influence of cadmium andcopper mixtures to rhizosphere bacterial communities. SoilEcol. Lett. 5, 94–107. doi: 10.1007/s42832-021-0128-9
Asnicar, F., Weingart, G., Tickle, T. L., Huttenhower, C., and Segata, N. (2015). Compact graphical representation of phylogenetic data and metadata with GraPhlAn. PeerJ 3, e1029. doi: 10.7717/peerj.1029
Bi, Y. M., Zhang, X. M., Jiao, X. L., Li, J. F., Peng, N., Tian, G. L., et al. (2023). The relationship between shifts in the rhizosphere microbial community and root rot disease in a continuous cropping American ginseng system. Front. Microbiol. 14. doi: 10.3389/fmicb.2023.1097742
Bokulich, N. A., Kaehler, B. D., Ram, R. J., Matthew, D., Evan, B., Rob, K., et al. (2018). Optimizing taxonomic classification of marker-gene amplicon sequences with qiime 2’s q2-feature-classifier plugin. Microbiome 6, 90. doi: 10.1186/s40168-018-0470-z
Cao, Y. J., Pu, Z. J., Tang, Y. P., Shen, J., Chen, Y. Y., Kang, A., et al. (2017). Advances in bio-active constituents, pharmacology and clinical applications of rhubarb. Chin. Med. 12, 36. doi: 10.1186/s13020-017-0158-5
Cao, X., Yuan, Q., Hu, C., Zhang, H., Sun, X., Yan, B., et al. (2025). Wild wisdom meets cultivation: comparative rhizomicrobiome analysis unveils the key role of Paraburkholderia in growth promotion and disease suppression in Coptis chinensis. Microbiome 13, 150. doi: 10.1186/s40168-025-02136-4
Chinese Pharmacopoeia Commission (2025). Pharmacopoeia of the People’s Republic of China, vol. I Vol. 24 (Beijing: China Pharmaceutical Science and Technology Press).
Duval, J., Pecher, V., Poujol., M., and Lesellier, E. (2016). Research advances for the extraction, analysis and uses of anthraquinones: a review. Ind. Crops Prod. 94, 812–833. doi: 10.1016/j.indcrop.2016.09.056
Feng, T. S., Yuan, Z. Y., Yang, R. Q., Zhao, S., Lei, F., Xiao, X. Y., et al. (2013). Purgative components in rhubarbs: adrenergic receptor inhibitors linked with glucose carriers. Fitoterapia 91, 236–246. doi: 10.1016/j.fitote.2013.09.020
Finzi Razavi, B. S., Zarebanadkouki, M., Blagodatskaya, E., and Adrien, C. (2016). Rhizosphere shape of lentil and maize:spatial distribution of enzyme activities. Soil Biol. Biochem. 96, 229–237. doi: 10.1016/j.soilbio.2016.02.020
Hacquard, S., Spaepen, S., Garrido-Oter, R., and Schulze-Lefert, P. (2017). Interplay between innate immunity and the plant microbiota. Annu. Rev. Phytopathol. 55, 565–589. doi: 10.1146/annurev-phyto-080516-035623
He, C., Liu, C., Liu, H. F., Wang, W. Q., Hou, J. L., and Li, X. N. (2022). Dual inoculation of dark septate endophytes and Trichoderma viride drives plant performance and rhizosphere microbiome adaptations of Astragalus mongholicus to drought. Environ. Microbiol. 24, 324–340. doi: 10.1111/1462-2920.15878
Huang, X. G., Li, M. Y., Yan, X. N., Yang, J. S., Rao, M. C., and Yuan, X. F. (2021). The potential of Trichoderma brevicompactum for controlling root rot on Atractylodes macrocephala. Can. J. Plant Pathol. 43, 794–802. doi: 10.1080/07060661.2021.1933602
Huson, D. H., Mitra, S., Ruscheweyh, H. J., Weber, N., and Schuster, S. C. (2011). Integrative analysis of environmental sequences using MEGAN4. Genome Res. 21, 1552–1560. doi: 10.1101/gr.120618.111
Jabborova, D., Mamarasulov, B., Davranov, K., Enakiev, Y., Bisht, N., Singh, S., et al. (2024). Diversity and plant growth properties of rhizospheric bacteria associated with medicinal plants. Indian J. Microbiol. 64, 409–417. doi: 10.1007/s12088-024-01275-w
Jia, H. M., Zheng, C. W., Wu, Y. R., Wang, H., and Yan, Z. Y. (2024). Metabolomic approach reveals the mechanism of synthetic communities to promote high quality and high yield of medicinal plants—danshen (Salvia miltiorrhiza Bge.). Chem. Biol. Technol. Agric. 11, 120. doi: 10.1186/s40538-024-00651-4
Jiang, S., Li, H., Zhang, L., Mu, W., Zhang, Y., Chen, T., et al. (2025). Generic Diagramming Platform (GDP): a comprehensive database of high-quality biomedical graphics. Nucleic Acids Res. 52, 1670–1676. doi: 10.1093/nar/gkae973
Jiang, Y., Wang, W., Xie, Q., Liu, N., Liu, L., Wang, D., et al. (2017). Plants transfer lipids to sustain colonization by mutualistic mycorrhizal and parasitic fungi. Science 356, 1172–1175. doi: 10.1126/science.aam9970
Jones, R. T., Robeson, M. S., Lauber, C. L., Hamady, M., Knight, R., and Fierer, N. (2009). A comprehensive survey of soil acidobacterial diversity using pyrosequencing and clone library analyses. ISME J. 3, 442–453. doi: 10.1038/ismej.2008.127
Jun, W., Shi, F., Baohui, L. U., Lina, N. Y., Xue, W., Yan, J. Z., et al. (2022). Fusarium oxysporum f. sp. ginseng, a new forma specialis causing Fusarium root rot of Panax ginseng. Phytopathol. Mediterr. 61, 417–429. doi: 10.36253/phyto-13723
Kapoor, R., Anand, G., Gupta, P., and Mandal, S. (2017). Insight into the mechanisms of enhanced production of valuable terpenoids by arbuscular mycorrhiza. Phytochem. Rev. 16, 677–692. doi: 10.1007/s11101-016-9486-9
Keshavarzi, Z., Shakeri, F., Maghool, F., Jamialahmadi, T., Johnston, T. P., and Sahebkar, A. (2021). A review on the phytochemistry, pharmacology, and therapeutic effects of rheum ribes. Adv. Exp. Med. Biol. 1328 (56), 447–461. doi: 10.1007/978-3-030-73234-9_30
Komatsu, K., Nagayama, Y., Tanaka, K., Ling, Y., Cai, S. Q., Omote, T., et al. (2006). Comparative study of chemical constituents of rhubarb from different origins. Chem. Pharm. Bull. (Tokyo) 54, 1491–1499. doi: 10.1248/cpb.54.1491
Kong, W. D., Zhu, Y. G., Fu, B. J., Marschner, P., and He, J. Z. (2006). The veterinary antibiotic oxytetracycline and Cu influence functional diversity of the soil microbial community. Environ. pollut. 143, 129–137. doi: 10.1016/j.envpol.2005.11.003
Li, S. Z., Deng, Y., Du, X. F., Wu, Y., He, Q., Wang, Z., et al. (2021). Sampling cores and sequencing depths affected the measurement of microbial diversity in soil quadrats. Sci. Total Environ. 50, 144966. doi: 10.1016/j.scitotenv.2021.144966
Li, Y. M., Wang, X. R., Yang, G., Yan, T. T., Hu, X. C., Peng, L., et al. (2025). Profiling the accumulation of ten bioactive compounds in Rheum officinale across different growth years using integrated HPLC and transcriptomic approaches. Med. Plant Biol. 4, e022. doi: 10.48130/mpb-0025-0020
Li Destri Nicosia, M. G., Mosca, S., Mercurio, R., and Schena, L. (2015). Dieback of Pinus nigra Seedlings Caused by a Strain of Trichoderma viride. Plant Dis. 99, 44–49. doi: 10.1094/PDIS-04-14-0433-RE
Lin, Y., Cai, Y., Hu, F., Wang, J. Y., Jiang, H. L., Zhou, S. S., et al. (2024). Study on the correlation of Rhubarb root rot disease with soil nutrients and microbial community composition. J. Southwest Univ. Nat. Sci. Ed. 46, 70–83. doi: 10.13718/j.cnki.xdzk.2024.03.006
Liu, S. Y., Wang, Q. Q., Lei, Y. H., Wang, S. S., Chen, K. L., Li, Y., et al. (2022). Elucidating the interaction of rhizosphere bacteria and environmental factors in influencing active ingredient content of Lycium barbarum fruit in China. J. Appl. Microbiol. 132, 3783–3796. doi: 10.1111/jam.15502
Lu, S., Gischkat, S., Reiche, M., Akob, D. M., Hallberg, K. B., and Küsel, K. (2010). Ecophysiology of Fe-cycling bacteria in acidic sediments. Appl. Environ. Microbiol. 76, 8174–8183. doi: 10.1128/AEM.01931-10
Luo, C., He, Y., and Chen, Y. (2024). Rhizosphere microbiome regulation: Unlocking the potential for plant growth. Curr. Res. Microb. Sci. 8, 100322. doi: 10.1016/j.crmicr.2024.100322
Lv, G., Li, Z., Zhao, Z., Liu, H., Li, L., and Li, M. (2024). The factors affecting the development of medicinal plants from a value chain perspective. Planta 259, 108. doi: 10.1007/s00425-024-04380-8
Lv, J., Yang, S., Zhou, W., Liu, Z., Tan, J., and Wei, M. (2024). Microbial regulation of plant secondary metabolites: Impact, mechanisms and prospects. Microbiol. Res. 283, 127688. doi: 10.1016/j.micres.2024.127688
Mohd, S., Kushwaha, A. S., Shukla, J., Mandrah, K., Shankar, J., Arjaria, N., et al. (2019). Fungal mediated biotransformation reduces toxicity of arsenic to soil dwelling microorganism and plant. Ecotoxicol. Environ. Saf. 176, 108–118. doi: 10.1016/j.ecoenv.2019.03.053
Nakaew, N., Lumyong, S., Sloan, W. T., and Sungthong, R. (2019). Bioactivities and genome insights of a thermotolerant antibiotics-producing Streptomyces sp. TM32 reveal its potentials for novel drug discovery. Microbiologyopen 8, e842. doi: 10.1002/mbo3.842
Pang, F., Li, Q., Solanki, M. K., Wang, Z., Xing, Y. X., and Dong, D. F. (2024). Soil phosphorus transformation and plant uptake driven by phosphate-solubilizing microorganisms. Front. Microbiol. 15. doi: 10.3389/fmicb.2024.1383813
Pang, Z., Mao, X., Xia, Y., Xiao, J., Wang, X., Xu, P., et al. (2022). Multiomics reveals the effect of root rot on polygonati rhizome and identifies pathogens and biocontrol strain. Microbiol. Spectr. 10, e0238521. doi: 10.1128/spectrum.02385-21
Pankratov, T. A., Kirsanova, L. A., Kaparullina, E. N., Kevbrin, V. V., and Dedysh, S. N. (2012). Telmatobacter bradus gen. nov., sp. nov., a cellulolytic facultative anaerobe from subdivision 1 of the Acidobacteria, and emended description of Acidobacterium capsulatum Kishimoto et al., 1991. Int. J. Syst. Evol. Microbiol. 62, 430–437. doi: 10.1099/ijs.0.029629-0
Peng, Z., Guo, X. Z., Xu, Y., Liu, D. H., Wang, H. Y., Guo, L. P., et al. (2020). Advances in interaction between medicinal plants and rhizosphere microorganisms. Chin. J. Chin. Mater. Med. 45, 2023–2030. doi: 10.19540/j.cnki.cjcmm.20200302.116
Perveen, I., Raza, M. A., Sehar, S., Naz, I., Young, B., and Ahmed, S. (2017). Heavy metal contamination in water,soil,and milk of the industrial area adjacent to Swan River, Islamabad, Pakistan. HERA 23, 1564–1572. doi: 10.1080/10807039.2017.1321956
Rosa-Mera, C. J. D. L., Ferrera-Cerrato, R., Alejandro, A. A., Sánchez-Colín, M. J., and Muñoz-Muñiz, O. D. (2011). Arbuscular mycorrhizal fungi and potassium bicarbonate enhance the foliar content of the vinblastine alkaloid in Catharanthus roseus. Plant Soil. 349, 367–376. doi: 10.1007/s11104-011-0883-y
Shi, Z., Guo, X., Lei, Z., Wang, Y., Yang, Z., Niu, J., et al. (2023). Screening of high-efficiency nitrogen-fixing bacteria from the traditional Chinese medicine plant Astragalus mongolicus and its effect on plant growth promotion and bacterial communities in the rhizosphere. BMC Microbiol. 23, 292. doi: 10.1186/s12866-023-03026-1
Song, H., Liu, J., Song, P., Feng, L., Wu, C. Z., and Hong, W. (2020). Effects of inoculation with endophytic fungi and endophytic bacteria on growth and accumulation of secondary metabolites in tripterygium wilfordii. J. Trop. Subtrop. Bot. 28, 347–355. doi: 10.11926/jtsb.4176
Sun, H., Jin, Q., Wang, Q., Shao, C., Zhang, L. L., Guan, Y. M., et al. (2020). Effects of soil quality on effective ingredients of Astragalus mongholicus from the main cultivation regions in China. Ecol. Indic. 114, 106296. doi: 10.1016/j.ecolind.2020.106296
Takayama, K., Tsutsumi, H., Ishizu, T., and Okamura, N. (2012). The influence of rhein 8-O-β-D-glucopyranoside on the purgative action of sennoside a from rhubarb in mice. Biol. Pharm. Bull. 35, 2204–2208. doi: 10.1248/bpb.b12-00632
Tohge, T. (2020). From fruit omics to fruiting omics: systematic studies of tomato fruiting by metabolic networks. Mol. Plant 13, 1114–1116. doi: 10.1016/j.molp.2020.07.012
Vaghela, N. and Gohel, S. (2023). Medicinal plant-associated rhizobacteria enhance the production of pharmaceutically important bioactive compounds under abiotic stress conditions. J. Basic Microbiol. 63, 308–325. doi: 10.1002/jobm.202200361
Wang, C., Pan, G., Lu, X., and Qi, W. (2023). Phosphorus solubilizing microorganisms: potential promoters of agricultural and environmental engineering. Front. Bioeng. Biotechnol. 11. doi: 10.3389/fbioe.2023.1181078
Welling, M. T., Liu, L., Rose, T. J., Waters, D. L., and Benkendorff, K. (2016). Arbuscular mycorrhizal fungi: effects on plant terpenoid accumulation. Plant Biol. (Stuttg) 18, 552–562. doi: 10.1111/plb.12408
Woźniak, M., Gałązka, A., Tyśkiewicz, R., and Jaroszuk-Ściseł, J. (2019). Endophytic bacteria potentially promote plant growth by synthesizing different metabolites and their phenotypic/Physiological profiles in the biolog GEN III microPlate™ Test. Int. J. Mol. Sci. 20, 5283. doi: 10.3390/ijms20215283
Wu, Y., Leng, F., Liao, M., Yu, Y., Chen, Z., Wei, S., et al. (2023). Characterization of the physiological parameters, effective components, and transcriptional profiles of Polygonum multiflorum Thunb. Under pH stress. Plant Physiol. Biochem. 37, 108279. doi: 10.1016/j.plaphy.2023.108279
Xia, Q., Wang, B., Liu, Z., Wei, F. G., Yang, S. Z., Li, X. C., et al. (2024). Deciphering the key soil microbial taxa that contribute to saponin accumulation in a geo-authentic Sanqi ginseng production area: Evidence from four different varieties. Appl. Soil Ecol. 201, 105470. doi: 10.1016/j.apsoil.2024.105470
Xiang, H., Zuo, J., Guo, F., and Dong, D. (2020). What we already know about rhubarb: a comprehensive review. Chin. Med. 15, 88. doi: 10.1186/s13020-020-00370-6
Xu, R., Chen, W., Chen, S., Wang, X., Xu, J., Zhang, Y., et al. (2025). Unraveling the rhubarb (Rheum officinale baill.) root and rhizosphere microbial communities in response to pathogen exposure. Mol. Biotechnol. 32, 1–11. doi: 10.1007/s12033-025-01367-y
Xu, H., Wang, W., Li, X., Li, Y., Jiang, Y., Deng, C., et al. (2024). Botany, traditional use, phytochemistry, pharmacology and clinical applications of rhubarb (Rhei radix et rhizome): A systematic review. Am. J. Chin. Med. 52, 1925–1967. doi: 10.1142/S0192415X24500757
Yang, J., Wang, F., Wen, Y., Gao, S., Lu, C., Liu, Y., et al. (2020). First report of Fusarium proliferatum causing root rot disease in Salvia miltiorrhizae in China. Plant Dis. 105, 1210–1210. doi: 10.1094/PDIS-09-20-1908-PDN
You, H., Yang, S., Zhang, L., Hu, X., and Li, O. (2017). Promotion of phenolic compounds production in Salvia miltiorrhiza hairy roots by six strains of rhizosphere bacteria. Eng. Life Sci. 18, 160–168. doi: 10.1002/elsc.201700077
Yuan, Y., Zuo, J., Zhang, H., Meng, T. Z., and Sian, L. (2022). The Chinese medicinal plants rhizosphere: Metabolites, microorganisms, and interaction. Rhizosphere 22, 100540. doi: 10.1016/j.rhisph.2022.100540
Zeng, M., Zhong, Y., Cai, S., and Diao, Y. (2018). Deciphering the bacterial composition in the rhizosphere of Baphicacanthus cusia (NeeS) Bremek. Sci. Rep. 8, 15831. doi: 10.1038/s41598-018-34177-1
Keywords: 16S rRNA sequencing, fungal ITS high-throughput sequencing, Rheum officinale Baill., secondary metabolites, soil properties
Citation: Wang Y, Yan F, Liang Y, Hu X, Gao J, Zhang M, Yan Y, Zhang G and Li Y (2025) Comparative analysis of rhizosphere microbial communities and secondary metabolites in cultivated Rheum officinale from different regions of China. Front. Plant Sci. 16:1650792. doi: 10.3389/fpls.2025.1650792
Received: 20 June 2025; Accepted: 08 September 2025;
Published: 30 September 2025.
Edited by:
Zhen Wang, Yunnan Agricultural University, ChinaReviewed by:
Kankan Zhao, University of California, San Diego, United StatesPooja Sharma, Dr. Yashwant Singh Parmar University of Horticulture and Forestry, India
Copyright © 2025 Wang, Yan, Liang, Hu, Gao, Zhang, Yan, Zhang and Li. 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: Gang Zhang, amF5X2d1bWxpbmcyMDI1QDE2My5jb20=; Yimin Li, bGl5aW1pbnR0QDE2My5jb20=