Skip to main content

Assessing the biodiversity of rhizosphere and endophytic fungi in Knoxia valerianoides under continuous cropping conditions

Abstract

Background

Rhizosphere and endophytic fungi play important roles in plant health and crop productivity. However, their community dynamics during the continuous cropping of Knoxia valerianoides have rarely been reported. K. valerianoides is a perennial herb of the family Rubiaceae and has been used in herbal medicines for ages. Here, we used high-throughput sequencing technology Illumina MiSeq to study the structural and functional dynamics of the rhizosphere and endophytic fungi of K. valerianoides.

Results

The findings indicate that continuous planting has led to an increase in the richness and diversity of rhizosphere fungi, while concomitantly resulting in a decrease in the richness and diversity of root fungi. The diversity of endophytic fungal communities in roots was lower than that of the rhizosphere fungi. Ascomycota and Basidiomycota were the dominant phyla detected during the continuous cropping of K. valerianoides. In addition, we found that root rot directly affected the structure and diversity of fungal communities in the rhizosphere and the roots of K. valerianoides. Consequently, both the rhizosphere and endophyte fungal communities of root rot-infected plants showed higher richness than the healthy plants. The relative abundance of Fusarium in two and three years old root rot-infected plants was significantly higher than the control, indicating that continuous planting negatively affected the health of K. valerianoides plants. Decision Curve Analysis showed that soil pH, organic matter (OM), available K, total K, soil sucrase (S_SC), soil catalase (S_CAT), and soil cellulase (S_CL) were significantly related (p < 0.05) to the fungal community dynamics.

Conclusions

The diversity of fungal species in the rhizosphere and root of K. valerianoides was reported for the first time. The fungal diversity of rhizosphere soil was higher than that of root endophytic fungi. The fungal diversity of root rot plants was higher than that of healthy plants. Soil pH, OM, available K, total K, S_CAT, S_SC, and S_CL were significantly related to the fungal diversity. The occurrence of root rot had an effect on the community structure and diversity of rhizosphere and root endophytic fungi.

Peer Review reports

Introduction

The rhizosphere is defined as a narrow area of soil around the roots that has the most significant effect on plant nutrition and growth [1, 2]. Fungal communities colonize both the plants and the rhizosphere soil affecting plant’s fitness in various ways, such as influencing nutrient uptake, inducing soil-borne diseases, and affecting the activity of plant pathogens [3]. Interactions between plant roots and rhizosphere microbiomes collectively form the plant-root microorganism complex [4]. Some agronomic practices can affect plant-microbial interactions, including continuous cropping [5,6,7], crop rotation [8], intercropping [9], irrigation and fertilization [10, 11], as well as the application of pesticide-fertilizer combinations [12]. Continuous cropping is widely practiced in Chinese agriculture, which involves the cultivation of the same crop in the same soil year after year [13]. Continuous cropping for many years causes loss in crop yields, deterioration in quality, increased susceptibility to pests and diseases, and reduction in soil fertility [3]. The harmful effects of continuous cropping have been reported on various rhizome medicinal plants, including Salvia miltiorrhiza [14], Panax notoginseng [15], and Coptis chinensis [16].

Knoxia valerianoides, also known as ‘Zi Daji’ in China, is a perennial herb belonging to the family Rubiaceae, and its root is a traditional Chinese medicine (TCM) owing to its medicinal properties [17]. The roots of K. valerianoides are the main ingredient in the Chinese patent medicine Zijinding [18]. The major organic compounds in the extracts of K. valerianoides are anthraquinones and triterpenoids [19, 20]. Studies have demonstrated that continuous cropping leads to the decrease in the yield and quality of S. miltiorrhiza [14], which may be attributed to the composition change of the rhizosphere soil microbial community [21]. Changes in soil microflora are considered the primary cause of yield decline in rhizome medicinal plants under continuous cropping conditions [22]. Root rot is a common disease affecting the cultivation of various root and rhizome medicinal plants, and has been reported on numerous occasions [23, 24]. In recent years, with the expansion of Chinese herbal medicine cultivation, root rot has been reported to increase with each passing year in continuous planting conditions. P. notoginseng is particularly susceptible to root rot disease, with an annual incidence ranging from 5 to 20%, causing severe losses sometimes up to 70% [25]. Similarly, in Canada 20-30% of ginseng crops are lost annually to root rot [26]. Studies have shown that most pathogens of root rot are fungi, including Fusarium oxysporum, F. solani, and Rhizoctonia solani [27,28,29], with Fusarium as the predominant pathogenic factor.

Rhizosphere fungi are involved in important ecological processes and thus are key players in the promotion of plant growth and maintenance of plant health [30, 31]. Some rhizosphere fungi are associated with pathogen growth inhibition and thus act as biocontrol agents; however, imbalances in rhizosphere fungal community structure have been associated with the development of root rot [32, 33]. Furthermore, endophytic fungi, characterized by their non-pathogenic nature, have shown promise as biological control agents [34]. The presence and abundance of endophytic fungi significantly influence the overall health of host plants, primarily through their positive impact on the plant’s immune system [35, 36]. Recent studies have indicated that continuous cropping practices can disrupt the balance of plant root microbial complex by the accumulation of pathogenic fungi in the rhizosphere [3, 37].

There are several reports exploring the plant root microbial complex of rhizome medicinal plants; however, there is no such report on K. valerianoides. In this study, we aimed to fill this knowledge gap. The primary objective was to investigate the diversity and composition of fungal communities within the plant root microbial complex of K. valerianoides using high-throughput sequencing technology. Second objective was to characterize rhizosphere and endophytic fungal communities of healthy and infected K. valerianoides under continuous cropping. Additionally, we aimed to analyze the soil physico-chemical properties and enzymatic activity, and determine their relationship with rhizosphere and endophytic fungal community structures. The overarching goal of this study was to gain insights into the pathogenesis of K. valerianoides root rot.

Materials & methods

Study site and experimental design

The field experiment was performed in Xiangyun County, Dali Bai Autonomous Prefecture of Yunnan Province (25°25′N, 100°40′E), a prominent region for the cultivation of K. valerianoides in China. In August 2022, root rot was observed on K. valerianoides root with about 15% incidence. The rhizosphere soil of 1, 2, and 3-years old K. valerianoides plants were collected by using a diagonal method [5]. A total of 6 plots were selected for sampling. Within each sample plot, rhizosphere soil samples of 5 healthy or root rot-infected K. valerianoides plants were combined to make a single sample, and the unplanted soil of the same plot served as control. The average depth of the main root and fibrous root of K. valerianoides was determined to be 0–30 cm. Consequently, sampling was conducted at a depth of 0–30 cm. During sampling, the plants were carefully excavated, and the loose soil that naturally fell off was discarded. The soil left on the root surface was considered the rhizosphere soil [5]. After collection, each sample was sieved using a 4 mm soil sieve, and the remaining root sample was quickly placed in a dry ice box storage. The sampling sites were all abandoned recreational land before the planting of K. valerianoides; after planting, they were managed using a mixed model of organic fertiliser and chemical fertiliser (continuous planting period of 3–4 years).

The root surface disinfection steps are as follows: The root surface was rinsed with flowing tap water, followed by three rinses with sterile distilled water. The root of K. valerianoides was then divided into small sections with a sterile scalpel for surface disinfection. Each tissue sample was immersed in 70% ethanol for 4 min, washed with fresh 2% NaClO solution for 1–2 min, subsequently immersed in 70% ethanol for 1 min, and finally rinsed three times with sterile distilled water [38].

Five replicated soil or root samples were combined to form a single sample. These samples included healthy rhizosphere soil from 1-year-old (S_H1), 2-year-old (S_H2), and 3-year-old (S_H3) K. valerianoides plants. Similarly, root rot affected rhizosphere soil from 2-year-old (S_D2), and 3-year-old (S_D3) K. valerianoides plants. A sample from unplanted soil of the same plot (S_CK) was used as a control. Additionally, endophyte samples of healthy roots from 1-year-old (R_H1), 2-year-old (R_H2), and 3-year-old (R_H3) K. valerianoides plants were collected, along with endophytic samples of root rot affected 2-year-old (R_D2), and 3-year-old (R_D3) K. valerianoides plants; while samples from 1-year-old plants (R_H1) were used as control. All the samples were stored at -80℃ until further analyses.

Soil physico-chemical properties and enzyme activity

The physico-chemical properties and enzyme activity of all soil samples were analyzed following established protocols. Soil pH was measured using a glass electrode pH meter (PHS-3E, China) in a 2.5:1 water: soil (v/w) suspension. Soil organic matter (OM) content was determined by K2Cr2O7 redox titration [39], soil hydrolytic nitrogen (HN) was quantified by alkaline hydrolysis, while soil available phosphorus (AP) was extracted with sodium bicarbonate and determined by molybdenum blue method [40]. The measurement of soil available potassium (AK) followed standard procedures outlined by MC Clean [41]. Established protocols were followed to measure total nitrogen (TN) [42], total phosphorus (TP), and total potassium (TK) [43] in the soil samples. To determine the enzymatic activity of soil samples; namely soil catalase (S_CAT), soil urease (S_UE), soil sucrase (S_SC), and soil cellulase (S_CL), specific kits (G0303F, G0301F, G0302F, G0308F, Grace) were employed. S_CAT, S_UE, and S_SC are soil extracellular enzymes whose activities are important bioactivity indicators for soil quality tests. S_CL play an important role in the degradation and decomposition of organic matter in soil. The assay procedures were performed following the kit instructions, and the measurements were taken using an ultraviolet spectrophotometer (Shanghai Jinghua 754).

DNA extraction, PCR amplification, and sequencing

Total microbial genomic DNA was extracted from the rhizosphere soil, blank control soil samples, and plant root samples using the PowerSoil Pro DNA Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. The quality and concentration of DNA were assessed by 1.0% agarose gel electrophoresis and a NanoDrop® ND-2000 spectrophotometer (Thermo Scientific Inc., USA) and kept at -80 ℃ until further analyses. For the amplification of the hypervariable region ITS of the fungal ITS gene, primer pairsITS1-F (5’- CTTGGTCATTTAGAGGAAGTAA-3’) and ITS2-R (5’-GCTGCGTTCTTCATCGATGC-3’) were utilized with an ABI GeneAmp® 9700 PCR thermocycler (ABI, CA, USA) [44]. The PCR reaction mixture consisted of 5×Fast Pfu buffer (4 µL), 2.5 mM dNTPs (2 µL), 5 µM each primer (0.8 µL), Fast Pfu polymerase (0.4 µL), template DNA (10 ng), and ddH2O was added to make a final volume of 20 µL. The PCR amplification procedure involved initial denaturation at 95 ℃ for 3 min, followed by 27 cycles of denaturing at 95 ℃ for 30 s, annealing at 55 ℃ for 30 s, and extension at 72 ℃ for 45 s, with a final single extension at 72 ℃ for 10 min, and subsequent cooling at 4 ℃. All samples were amplified in triplicate, and the resulting amplicons were extracted from 2% agarose gel and purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA), following the manufacturer’s instructions. The quantification of the extracted amplicon was performed using a Quantus™ fluorescence meter (Promega, USA). The sequencing library was prepared by Illumina according to their established protocols. Subsequently, the prepared library underwent qualification to assess its quality and integrity. For sequencing, the Illumina MiSeq PE300 platform (Illumina, San Diego, USA) was employed, and the sequencing service was entrusted to Shanghai Majorbio Bio-Pharm Technology Co., Ltd., following their standard procedures.

Data processing and statistical analyses

The sequencing raw reads obtained from the Illumina MisSeq PE300 platform were subjected to quality control using fastp (version 0.19.6) to remove low-quality reads. The remaining high-quality reads were further processed for merging using FLASH (version 1.2.11) to obtain longer sequences [45, 46]. The optimized sequences were then clustered into operational taxonomic units (OTUs) using UPARSE 7.1 with a 97% sequence similarity threshold [47]. To account for variations in sequencing depth, the number of sequences in each sample was normalized to 20,000 reads, ensuring comparable coverage among samples while maintaining an average sequence coverage of 99.9%. The taxonomic classification of each representative OTU sequence was performed using the RDP Classifier (version 2.11, http://sourceforge.net/projects/rdp-classifier/) against the ITS fungi database (release 8.0, http://unite.ut.ee/index.php) using a confidence threshold of 0.7. Alpha diversity metrics including Sobs, Chao, Shannon, Simpson, and Coverage were calculated by mothur software (version 1.30.2, http://www.mothur.org/wiki/Calculators), and the difference in Alpha diversity between groups was analyzed by Student’s t-test [48]. Statistical analyses and data visualization were performed in R language (version 3.3.1).

Non-metric multidimensional scaling (NMDS) analysis based on the abund-jaccard distance algorithm was employed to explore similarities in microbial community structure among samples, and a heatmap diagram was generated to visualize the differences in species composition. Additionally, Decision Curve Analysis (DCA) was performed to investigate the impact of soil environmental factors on the fungal community structure in the rhizosphere roots. Data correlation analysis was carried out using SPSS 19.0.

Results

Soil physico-chemical properties and enzyme activity

The analysis of soil physico-chemical properties and enzyme activity revealed notable effects of continuous cropping on the soil. Specifically, the pH of the continuous cropping soil showed a significant increase (p < 0.05) compared to the unplanted soil (S_CK). Furthermore, the continuous cropping soil exhibited significantly higher levels of OM, hydrolyzable nitrogen (N), available phosphorus (P), total nitrogen (N), and total phosphorus (P) in comparison to the unplanted soil. Additionally, the enzyme activities of soil sucrase (S_SC) and soil urease (S_UE) were significantly higher in the continuous cropping soil than in the unplanted soil. S_SC increased with the increase of planting years. Available potassium (K) and S_UE showed that healthy plants were higher than infected plants, while S_CAT, pH and S_CL showed the opposite (Table 1).

Table 1 Soil physico-chemical properties and enzyme activity of rhizospheric soils in Xiangyun county, China (n = 6). Different letters in the same column indicate significant differences at p < 0.05

Composition and α-Diversity of rhizosphere and endophytic fungi

Composition and α-diversity analyses were performed to investigate the fungal communities in the rhizosphere and endophytic compartments. High-throughput sequencing was performed on collected samples, including six rhizosphere soil fungi (S_H1, S_H2, S_H3, S_D2, S_D3, S_CK) and five root endophytic fungi samples (R_H1, R_H2, R_H3, R_D2, R_D3) with six biological replicates per sample. A total of 3,211,327 ITS sequences were obtained after quality control, resulting in 769,901,943 bases with an average length of 239 bp. Subsequent normalization analysis yielded 33,872 sequences per sample, generating 5,312 OTUs. Of these, 1,953 OTUs were found in the root samples of K. valerianoides, 522 OTUs were shared among the rhizosphere soil samples of K. valerianoides, and 12 OTUs were shared among the root samples of K. valerianoides (Fig. 1A and B).

Fig. 1
figure 1

A and B Venn diagrams at the OTU level in root (A) and soil (B) groups. Rarefaction curves(C) of fungal communities based on observed operational taxonomic units (OTUs) for 11 plant and soil samples under a Knoxia valerianoides continuous cropping system. Repeat 6 times for each sample

The dilution curve was plotted using randomly selected sequencing data to assess Alpha diversity, with sequencing data as the abscissa and Alpha diversity index value as the ordinate (Fig. 1C). The curve of each sample displayed a plateau, indicating that the obtained reads adequately represent the fungal community’s diversity, and the coverage of all samples exceeded 99.31%.

Alpha diversity was estimated by calculating richness indices such as Sobs, Chao, Ace, Simpson, and others (Table 2). Overall, rhizosphere soil fungi exhibited higher richness and diversity indices compared to root endophytic fungi. Within the same planting period, the diversity indexes (Sobs, Chao, Simpson, and Ace) of endophytic fungi in root rot-infected roots were higher compared to the healthy plants, with R_D3 showing higher diversity than R_D2. Healthy plant R_H1 showed higher richness and diversity, while the Ace index decreased gradually with the increasing continuous cropping period. Similar patterns were observed for rhizosphere fungal diversity indices in relation to the root rot disease. The diversity indexes Ace and Chao of the healthy rhizosphere fungi increased with the age of the plant and were significantly higher (p < 0.05) than S_H1. The rhizosphere fungal diversity indices of root rot-affected plants, including Sobs, Ace, and Chao were significantly higher than those of S_CK. The field survey found that the incidence of 1–3 year old K. valerianoides root rot gradually increased with the increase of planting years, which was basically consistent with the change of fungi diversity with years.

Table 2 The diversity indexes of endophytic fungi and rhizosphere fungi were tested by LSD in One-way ANOVA at Genus level (n = 6). Different letters in the same column indicate significant differences at p < 0.05 among the samples from same source

Structure of the fungal communities

The structure of the fungal communities was investigated, revealing a total of 883 genera, 365 families, 152 orders, 61 classes, and 16 phyla in the rhizosphere soil and root samples of K. valerianoides. As illustrated in Fig. 2, the dominant phyla were Ascomycota, Basidiomycota, Olpidiomycota, and Mortierellomycota. Notably, the abundance of Ascomycota in R_H3 and R_D3 was the highest, accounting for 90.86% and 98.88%, respectively.

Fig. 2
figure 2

Relative species abundance at the phylum level in fungi communities, visualized horizontally. Based on the species annotation results, we choose to examine the distribution of relative species abundance among the top-ranked phyla; the relative abundance distribution of the top ten species is shown

Hierarchical cluster analysis was performed on the top 35 genera based on the fungal abundance in 11 plant and soil samples, and the results are presented in Fig. 3A. The analysis revealed that all rhizosphere soil samples were clustered together into one category, while the root samples were divided into three distinct categories. The relative abundance of certain fungal genera in soil and root samples differed significantly across different continuous cropping years. Solicoccozyma (8%) was the dominant genus in S_CK, but its relative abundance decreased with increasing continuous cropping time. Mortierella and Didymellaceae were dominant genera in S_H3 (10%), S_D2 (8%), and S_D3 (9%, 12%), with S_D3 exhibiting significantly higher abundance than S_D2. Chaetomiaceae was the dominant genus in S_H1 (33%), Olpidium in S_H2 (12%) and R_H2 (73%), Phialophora in R_H3 (64%) and R_D3 (41%), Exophiala in R_D2 (28%), and Issatchenkia in R_H1(21%).

Fig. 3
figure 3

The top 35 most abundant fungal communities across 11 plant and soil samples (A) and non-metric multidimensional scaling (NMDS) based on the operational taxonomic unit (OTUs) composition (B). Repeat 6 times for each sample

To visually represent the differences between different treatments, NMDS analysis based on abund-jaccard dissimilarity distance was performed on the sample data. The results demonstrated significant differences among different samples (R = 0.7925, p = 0.001), and the stress value was 0.152, indicating that continuous cropping exerted a significant impact on the composition of the fungal community (Fig. 3B).

Correlation between soil properties, enzyme activity, and the fungal communities

The correlation between soil properties, enzyme activity, and fungal community composition was examined using DCA. The results revealed that the response of fungal community composition to soil properties followed a unimodal model (Axis_lengths = 8.05). Consequently, Canonical Correspondence Analysis (CCA) was employed to investigate the relationship between the fungal community and soil environmental factors. It was observed that the soil environmental factors such as pH, total P, total K, total N, and S_CL exhibited correlations among themselves, as did OM, hydrolytic N, available K, S_UE, and S_CL. Furthermore, pH (p = 0.001), OM (p = 0.001), available K (p = 0.001), total K (p = 0.001), total P (p = 0.018), S_CAT (p = 0.001), S_SC (p = 0.001) and S_CL (p = 0.001) demonstrated significant correlations with fungal community structure (Fig. 4A). The Spearman correlation heatmap indicated that the relative abundance of the genus Olpidium exhibited significant negative correlations with S_UE, S_CL, available K, OM and hydrolytic N while displaying positive correlations with total K, pH, and S_CAT. Moreover, the relative abundance of the genus Dactylonectria showed a significant positive correlation with total K, pH, and S_CAT. Genus Phialophora exhibited significant positive correlations with OM, hydrolytic N, and S_SC (Fig. 4B).

Fig. 4
figure 4

Canonical correspondence analysis (CCA) between the soil properties, enzyme activity and fungal community (A). Spearman correlation heatmap between the soil properties, enzyme activity and fungal community (B)

Discussion

The sustainable production of K. valerianoides necessitates a profound understanding of the fungal community’s structure and diversity. In this study, we employed high-throughput sequencing technology to assess the impact of continuous planting on the rhizosphere soil and root endophytic fungal communities of K. valerianoides. Our findings revealed that continuous cropping exerted a substantial influence on the diversity, structure, and composition of rhizosphere soil and endophytic fungal communities. Following continuous cropping, the Alpha diversity indices of Sobs, Ace, Chao, and Simpson exhibited an increase in the rhizosphere soil, implying that continuous cropping of K. valerianoides enhanced the diversity and richness of fungi in rhizosphere soil. This observation is consistent with previous studies involving P. notoginseng where, the diversity of fungi in the rhizosphere soil increased with the increase of continuous cropping period [37]. Conversely, the diversity index Ace and Chao of endophytic fungi in healthy plant roots exhibited a gradual decrease with increasing period of continuous cropping, along with a decline in the Shannon index. Root rot has traditionally been considered the primary concern in continuous cropping practices for cultivating rhizome medicinal plants [49]. In comparison to the infected plants, we observed relatively lower fungal diversity in healthy plants and rhizosphere soil samples. Consequently, the reduction in the diversity of rhizosphere and endophytic fungi, especially the decrease in beneficial taxa [50], may contribute to the occurrence of diseases during the process of continuous cropping.

Continuous cropping of K. valerianoides significantly influenced the community composition of rhizosphere soil fungi. Ascomycota and Basidiomycota, known as key decomposers among soil fungal communities [49], were identified as the dominant phyla in all rhizosphere soil samples, confirming their ecological importance. This observation is consistent with the previous findings [30]. At the genus level, Fusarium and Gibberella exhibited an increase in abundance with prolonged continuous cropping periods. Fusarium species are notorious plant pathogens [51], including important pathogenic species of F. oxysporum and F. solani [27]. Prior studies have indicated that the disturbance of rhizosphere microbial communities, particularly the proliferation of potential pathogens, could contribute to the transition of the rhizosphere microenvironment from ‘health’ to ‘disease’ [52]. Thus, the observed changes in rhizosphere soil fungal communities, including the increased abundance of Fusarium and Gibberella, suggest a potential association between continuous cropping and the development of pathogenic conditions.

In recent years, the significance of endophytic fungi in plant tissues for plant health has garnered considerable attention. During continuous cropping, the diversity of endophytic fungi in roots was found to be lower than that in rhizosphere soil. Furthermore, a decline in endophytic fungal diversity indices Chao, Ace, and Sobs was observed in infected root samples as the duration of planting increased, indicating a reduction in endophytic fungal abundance. Ascomycota and Basidiomycota were identified as the dominant phyla in both healthy and infected K. valerianoides plants, mirroring the fungal distribution observed in rhizosphere soil. Ascomycetes was the most prevalent taxon, as it plays an important ecological role as a decomposer [53]. Previous studies have demonstrated that endophytic fungi primarily originate from soil fungi and gain entry into plants through roots, stems, and leaves [37].

Given that K. valerianoides roots are used in traditional medicine, it is crucial to note that certain coexisting endophytic fungi can significantly impact the formation of plant metabolites, consequently influencing the quality of medicinal plant raw materials [54]. Similarly, in case of marijuana (Cannabis sativa L.), the plants with pathogen-infected roots produced significantly reduced crop yield and quality [55]. At the genus level, there was a notable increase in the abundance of Fusarium, Pyrenochaeta, and Exophiala in root rot-infected samples. The proliferation of pathogens associated with root rot, such as F. oxysporum, F. solani, and Phytophthora cactorum, might have contributed to the reduction in endophytic fungi. This may be related to the fact that K. valerianoides planting improves soil organic matter content. These findings further substantiate the long-standing recognition of the genus Fusarium as an important plant pathogen [51]. Collectively, these outcomes serve as a foundation for future investigations concerning microbial community changes in K. valerianoides and their potential role as an indicator of plant health during continuous cropping.

The comprehensive changes in soil properties are intimately associated with continuous cropping challenges and crop productivity [22, 56]. Prolonged monoculture practices have had a detrimental impact on sustainable productivity, leading to a decrease in soil pH and OM content [57]. Long-term fertilization practices (32 years) can influence soil enzyme activity as well as physical and chemical properties [58]. Previous studies have indicated a correlation between the occurrence of ginseng root rot OM, available potassium, alkali-hydrolyzed nitrogen, pH, sucrase, and catalase activity in rhizosphere soil [59].

Microbial diversity encompasses the species composition and interspecific variations within microbial communities, reflecting the richness of soil microbial life, and serving as an indicator of soil quality and restoration potential [60]. Soil enzymes are an important biological indicator of the prevalence of material and energy metabolism in soil and soil quality. S_CAT, S_UE, S_SC and S_CL play important roles in the cycling and transformation of soil nutrients such as nitrogen and carbon, as well as in plant protection [61]. The S_UE activity of infected plants was lower than that of healthy plants, which was consistent with the results of previous studies [62]. In our study, soil pH, OM, available K, total K, S_CAT, S_SC, and S_CL exhibited the strongest correlations with fungal community structure. Soil pH has a profound influence on the chemical and physical properties of soil and serves as a primary factor shaping the microbial community composition in rhizosphere soil [63]. Notably, the rhizosphere soil of root rot-infected K. valerianoides exhibited higher pH levels compared to healthy soil. Furthermore, HN played a crucial role in fungal community composition under continuous cropping. The nitrogen fertilization system significantly impacts the composition of the soil fungal community, and the application of organic fertilizer can mitigate the abundance of soil fungal pathogens such as Fusarium [64]. Moreover, low irrigation practices can enhance soil enzyme activity [65], while high saline-alkali stress leads to a decrease in soil enzyme activity [66]. In our study, a significant correlation was observed between soil pH and S_CL, which may also contribute to the occurrence of root rot of K. valerianoides. It is speculated that pH and S_CL could be important predictors of fungal community composition in continuous cropping soil.

Conclusions

In conclusion, continuous cropping of K. valerianoides significantly altered the structure and diversity of both the rhizosphere soil and root endophytic fungal communities. Continuous cropping led to an increase in the abundance of pathogenic fungi, consequently, promoting soil-borne diseases in K. valerianoides. The rhizosphere soil fungal community of root rot-infected plants exhibited higher richness compared to the control (S_CK). Severe root rot infection displayed a negative effect on the diversity of rhizosphere and root fungal communities. Moreover, continuous cropping resulted in a significant increase in soil pH, OM, hydrolyzable nitrogen, available phosphorus, total N, total P, and soil urease activity. Additionally, soil sucrose activity showed an increasing trend with the duration of the continuous cropping period. It is a compelling evidence to suggest important relationships of soil physico-chemical properties and enzyme activity with the structural dynamics of fungal communities. Further investigations on the rhizosphere and endophytic fungi of both healthy and infected plants are needed for comprehensive understanding of pathogenesis of root rot disease in K. valerianoides.

Data availability

Original data have been deposited into the NCBI SRA database with the accession number PRJNA981433.

References

  1. Mendes R, Garbeva P, Raaijmakers JM. The rhizosphere microbiome: significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms. FEMS Microbiol Rev. 2013;37.

  2. Lakshmanan V, Selvaraj G, Bais HP. Functional soil microbiome: belowground solutions to an aboveground problem. Plant Physiol. 2014;166:689–700.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Zhu B, Wu J, Ji Q, Wu W, Dong S, Yu J, et al. Diversity of rhizosphere and endophytic fungi in Atractylodes macrocephala during continuous cropping. PeerJ. 2020;8:e8905.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Philippot L, Raaijmakers JM, Lemanceau P, van der Putten WH. Going back to the roots: the microbial ecology of the rhizosphere. Nat Rev Microbiol. 2013;11:789–99.

    Article  CAS  PubMed  Google Scholar 

  5. Tan Y, Cui Y, Li H, Kuang A, Li X, Wei Y, et al. Diversity and composition of rhizospheric soil and root endogenous bacteria in Panax notoginseng during continuous cropping practices. J Basic Microbiol. 2017;57:337–44.

    Article  CAS  PubMed  Google Scholar 

  6. Xu J, Ou X-H, Jiang W-K, Yuan Q-S, Wang Y-H, Yang J, et al. [Effect of Gastrodiae elata-phallus impudicus sequential planting pattern on soil microbial community structure]. China J Chin Materia Med. 2020;45:463–71.

    Google Scholar 

  7. Wang T, Yang K, Ma Q, Jiang X, Zhou Y, Kong D, et al. Rhizosphere microbial community diversity and function analysis of cut chrysanthemum during continuous monocropping. Front Microbiol. 2022;13:801546.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Lin G, Wan D, Yang X, Zhao K, Zhu Y, Yan Z. [Characteristics of soil microbial variation during crop rotation period at cultivation area of Salvia Miltiorrhiza in Zhongjiang of Sichuan province]. China J Chin Materia Med. 2009;34:3184–7.

    Google Scholar 

  9. Noushahi HA, Zhu Z, Khan AH, Ahmed U, Haroon M, Asad M, et al. Rhizosphere microbial diversity in rhizosphere of Pinellia ternata intercropped with maize. 3 Biotech. 2021;11:469.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Liu Z, Liu P, Jia X, Cheng Y, Dong S, Zhao B, et al. [Effects of irrigation and fertilization on soil microbial properties in summer maize field]. Chin J Appl Ecol. 2015;26:113–2.

    Google Scholar 

  11. Xie Y, Wang Z, Cheng X, Qiu R, Hamoud YA, Hong C, et al. Dissecting the combined effects of cultivar, fertilization, and irrigation on rhizosphere bacterial communities and nitrogen productivity in rice. Sci Total Environ. 2022;835:155534.

    Article  CAS  PubMed  Google Scholar 

  12. Huang W, Lu Y, Chen L, Sun D, An Y. Impact of pesticide/fertilizer mixtures on the rhizosphere microbial community of field-grown sugarcane. 3 Biotech. 2021;11:210.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Shipton PJ. Monoculture and soilborne plant pathogens. Annu Rev Phytopathol. 1977;15:387–407.

    Article  Google Scholar 

  14. Liu W, Zhang L, Zhang Y-Y, Zhou J, Sun Q-L, Wang X, et al. [Influence of continuous cropping years on yield and active compounds in Salvia miltiorrhiza f. Alba]. China J Chin Materia Med. 2013;38:4252–6.

    Google Scholar 

  15. Bao L, Liu Y, Ding Y, Shang J, Wei Y, Tan Y, et al. Interactions between phenolic acids and microorganisms in rhizospheric soil from continuous cropping of Panax notoginseng. Front Microbiol. 2022;13:791603.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Alami MM, Xue J, Ma Y, Zhu D, Abbas A, Gong Z, et al. Structure, function, diversity, and composition of fungal communities in rhizospheric soil of Coptis chinensis franch under a successive cropping system. Plants (Basel). 2020;9:E244.

    Google Scholar 

  17. Chen X-J, Pu X-Y, Pu X-M, Li X, Liu Z-B, Mei M-J, et al. Extracts of Knoxia Roxburghii (Spreng.) M. A. Rau induce apoptosis in human MCF-7 breast cancer cells via mitochondrial pathways. Molecules. 2022;27:6435.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Pu X-Y, Chen X-J, Wang X-G, Li G-D, Zhang X-M, Qiu B. Estimation of genome size of Knoxia roxburghii by flow cytometry and genome survey. Mol Plant Breed. 2022; 1–13.

  19. Zhao F, Wang S-J, Lin S, Zhu C-G, Yuan S-P, Ding X-Y, et al. Anthraquinones from the roots of Knoxia valerianoides. J Asian Nat Prod Res. 2011;13:1023–9.

    Article  CAS  PubMed  Google Scholar 

  20. Zhao F, Wang S, Lin S, Zhu C, Yue Z, Yu Y, et al. Natural and unnatural anthraquinones isolated from the ethanol extract of the roots of Knoxia valerianoides. Acta Pharm Sinica B. 2012;2:260–6.

    Article  Google Scholar 

  21. Jie W, Yang D, Yao Y, Guo N. Effects of rhizophagus intraradices on soybean yield and the composition of microbial communities in the rhizosphere soil of continuous cropping soybean. Sci Rep. 2022;12:17390.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Li M, Chen Z, Qian J, Wei F, Zhang G, Wang Y, et al. Composition and function of rhizosphere microbiome of Panax notoginseng with discrepant yields. Chin Med. 2020;15:85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Latha P, Anand T, Prakasam V, Jonathan EI, Paramathma M, Samiyappan R, Combining. Pseudomonas, Bacillus and Trichoderma strains with organic amendments and micronutrient to enhance suppression of collar and root rot disease in physic nut. Applied Soil Ecology. 2011;49 none:215–23.

  24. Gordo SMC, Pinheiro DG, Moreira ECO, Rodrigues SM, Poltronieri MC, de Lemos OF, et al. High-throughput sequencing of black pepper root transcriptome. BMC Plant Biol. 2012;12:168.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Wu HZ, Zeng Y, Guo L, Huang L, Chen B. Molecular characterization of microbial communities in the rhizosphere soils and roots of diseased and healthy Panax notoginseng. Antonie Van Leeuwenhoek. 2015;108:1059–74.

    Article  PubMed  Google Scholar 

  26. Walsh JP, DesRochers N, Renaud JB, Seifert KA, Yeung KK-C, Sumarah MW. Identification of N,N’,N″-triacetylfusarinine C as a key metabolite for root rot disease virulence in American ginseng. J Ginseng Res. 2021;45:156–62.

    Article  PubMed  Google Scholar 

  27. Latiffah Z, Hayati M, Baharuddin S, Maziah Z. Identification and pathogenicity of Fusarium species associated with root rot and stem rot of Dendrobium. Asian J Plant Pathol. 2010;3:14–21.

    Article  Google Scholar 

  28. Scherm B, Balmas V, Spanu F, Pani G, Delogu G, Pasquali M, et al. Fusarium Culmorum: causal agent of foot and root rot and head blight on wheat: the wheat pathogen Fusarium Culmorum. Mol Plant Pathol. 2013;14:323–41.

    Article  CAS  PubMed  Google Scholar 

  29. Chittem K, Mathew FM, Gregoire M, Lamppa RS, Goswami RS. Identification and characterization of Fusarium spp. associated with root rots of field pea in North Dakota. Eur J Plant Pathol. 2015;143:641–9.

    Article  Google Scholar 

  30. Miao C-P, Mi Q-L, Qiao X-G, Zheng Y-K, Chen Y-W, Xu L-H, et al. Rhizospheric fungi of Panax notoginseng: diversity and antagonism to host phytopathogens. J Ginseng Res. 2016;40:127–34.

    Article  PubMed  Google Scholar 

  31. Dong L-L, Niu W-H, Wang R, Xu J, Zhang L-J, Zhang J, et al. [Changes of diversity and composition of fungal communities in rhizosphere of Panax ginseng]. China J Chin Materia Med. 2017;42:443–9.

    Google Scholar 

  32. Yu M, Jiang J-L, Ren X-M, Li L, Jiao C-J, Yang L-J, et al. [Research on relationship between occurrence of root rot and changes of fungal communities in rhizosphere of Panax quinquefolius]. China J Chin Materia Med. 2018;43:2038–47.

    Google Scholar 

  33. Silva-Valderrama I, Toapanta D, Miccono M, de LA, Lolas M, Díaz GA, Cantu D, et al. Biocontrol potential of grapevine endophytic and rhizospheric fungi against trunk pathogens. Front Microbiol. 2020;11:614620.

    Article  PubMed  Google Scholar 

  34. Collinge DB, Jensen B, Jørgensen HJ. Fungal endophytes in plants and their relationship to plant disease. Curr Opin Microbiol. 2022;69:102177.

    Article  CAS  PubMed  Google Scholar 

  35. Zhu B, Wu L, Wan H, Yang K, Si J, Qin L. Fungal elicitors stimulate biomass and active ingredients accumulation in Dendrobium catenatum plantlets. Biologia. 2018;73:917–26.

    Article  CAS  Google Scholar 

  36. Yan L, Zhu J, Zhao X, Shi J, Jiang C, Shao D. Beneficial effects of endophytic fungi colonization on plants. Appl Microbiol Biotechnol. 2019;103:3327–40.

    CAS  PubMed  Google Scholar 

  37. Tan Y, Cui Y, Li H, Kuang A, Li X, Wei Y, et al. Rhizospheric soil and root endogenous fungal diversity and composition in response to continuous Panax notoginseng cropping practices. Microbiol Res. 2017;194:10–9.

    Article  PubMed  Google Scholar 

  38. Zhang C, Ma X, Zhu R, Liu Z, Gu M, Zhang J, et al. Analysis of the endophytic bacteria community structure and function of Panax notoginseng based on high-throughput sequencing. Curr Microbiol. 2020;77:2745–50.

    Article  CAS  PubMed  Google Scholar 

  39. Schnitzer M. Total carbon, organic matter, and carbon. 1982.

  40. Dane JH, Hopmans JH. Methods of soil analysis. 2002.

  41. Mclean EO, Watson ME. Soil measurements of plant-available potassium. Potassium in Agriculture; 1985.

  42. Singh BK, Bardgett RD, Smith P, Reay DS. Microorganisms and climate change: terrestrial feedbacks and mitigation options. Nat Rev Microbiol. 2010;8:779–90.

    Article  CAS  PubMed  Google Scholar 

  43. Lopes D, Rego F, Nunes C et al. Growth, soil properties and foliage chemical analysis comparison between pure and mixed stands of Castanea sativa Mill. and Pseudotsuga menziesii (Mirb.) Franco, in Northern Portugal. (Special Issue: Mixed and pure forests in a changing world.). Forest Systems. 2011;20:496–507.

  44. Adams RI, Miletto M, Taylor JW, Bruns TD. Dispersal in microbes: fungi in indoor air are dominated by outdoor air and show dispersal limitation at short distances. ISME J. 2013;7:1262–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27:2957–63.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Chen S, Zhou Y, Chen Y, Gu J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34:884–90.

    Article  Google Scholar 

  47. Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10:996–8.

    Article  CAS  PubMed  Google Scholar 

  48. Schloss, Sl W, Eb TR Jr. H, Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75.

  49. Dai L, Singh SK, Gong H, Tang Y, Peng Z, Zhang J, et al. Rhizospheric microbial consortium of Lilium lancifolium Thunb. Causes lily root rot under continuous cropping system. Front Microbiol. 2022;13:981615.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Ji L, Nasir F, Tian L, Chang J, Sun Y, Zhang J, et al. Outbreaks of root rot disease in different aged American ginseng plants are associated with field microbial dynamics. Front Microbiol. 2021;12:676880.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Moss MO, Smith JE. The applied mycology of Fusarium. Symposium of the british mycological society held at queen mary college, London, september 1982. Applied Mycology of Symposium of the British Mycological Society. 1984.

  52. Wang R, Wang Y, Yang Q, Kang C, Li M. Unraveling the characteristics of the microbial community and potential pathogens in the rhizosphere soil of Rehmannia glutinosa with root rot disease. Appl Soil Ecol. 2018;130:S0929139318301732.

    Article  Google Scholar 

  53. Miadlikowska J, Kauff F, Hofstetter V, Fraker E, Grube M, Hafellner J, et al. New insights into classification and evolution of the Lecanoromycetes (Pezizomycotina, Ascomycota) from phylogenetic analyses of three ribosomal RNA- and two protein-coding genes. Mycologia. 2006;98:1088–103.

    Article  CAS  PubMed  Google Scholar 

  54. Jia M, Chen L, Xin H-L, Zheng C-J, Rahman K, Han T, et al. A friendly relationship between endophytic fungi and medicinal plants: a systematic review. Front Microbiol. 2016;7:906.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Punja ZK, Collyer D, Scott C, Lung S, Holmes J, Sutton D. Pathogens and molds affecting production and quality of Cannabis sativa L. Front Plant Sci. 2019;10:1120.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Koskei K, Munyasya AN, Wang Y-B, Zhao Z-Y, Zhou R, Indoshi SN, et al. Effects of increased plastic film residues on soil properties and crop productivity in agro-ecosystem. J Hazard Mater. 2021;414:125521.

    Article  CAS  PubMed  Google Scholar 

  57. Zhao Q, Xiong W, Xing Y, Sun Y, Lin X, Dong Y. Long-term coffee monoculture alters soil chemical properties and microbial communities. Sci Rep. 2018;8:6116.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Ma N, Li T, Wu C, Zhang E. [Effects of long-term fertilization on soil enzyme activities and soil physicochemical properties of facility vegetable field]. Chin J Appl Ecol. 2010;21:1766–71.

    CAS  Google Scholar 

  59. Ma Y, Xu H-Y, Song M-H, Feng Z-W, Shao Z-W, Lu B-H et al. Relationship analysis between ginseng root rot diseases and main chemical properties, enzyme activities of rhizosphere soil. J Jilin Agricultural Univ. 2021;1–10.

  60. Zhou X, Yu G, Wu F. Responses of soil microbial communities in the rhizosphere of cucumber (Cucumis sativus L.) to exogenously applied p-hydroxybenzoic acid. J Chem Ecol. 2012;38:975–83.

    Article  PubMed  Google Scholar 

  61. Guan H-L, Yang J-Z, Chen Y-J, Cui X-M, Wany, Zhang Y-F. Change of rhizospheric microbe colony in cultivated soil and its correlation to root rot disease in Panax notoginseng. Soils. 2010;42:378–84.

    Google Scholar 

  62. Wang R, Zhang H, Sun L, Qi G, Chen S, Zhao X. Microbial community composition is related to soil biological and chemical properties and bacterial wilt outbreak. Sci Rep. 2017;7:343.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Lauber C, Hamady M, Knight R, Fierer N. Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl Environ Microbiol. 2009;75.

  64. Wen Y-C, Li H-Y, Lin Z-A, Zhao B-Q, Sun Z-B, Yuan L, et al. Long-term fertilization alters soil properties and fungal community composition in fluvo-aquic soil of the North China Plain. Sci Rep. 2020;10:7198.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Muhammad I, Yang L, Ahmad S, Zeeshan M, Farooq S, Ali I, et al. Irrigation and nitrogen fertilization alter soil bacterial communities, soil enzyme activities, and nutrient availability in maize crop. Front Microbiol. 2022;13:833758.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Yang D, Tang L, Cui Y, Chen J, Liu L, Guo C. Saline-alkali stress reduces soil bacterial community diversity and soil enzyme activities. Ecotoxicology. 2022;31:1356–68.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank MogoEdit (https://www.mogoedit.com) for its English editing during the preparation of this manuscript.

Funding

This research was financially supported by the major science and technology special program of Yunnan Province (202102AA310037; 202102AA310045-04); Yunnan Provincial Science and Technology Department Applied Basic Research Joint Special Funds of Chinese Medicine (202001AZ070001-015). Author BQ has received research support. The funders had no infuence on the contents of this manuscript.

Author information

Authors and Affiliations

Authors

Contributions

All authors contributed to the study conception and design. CL and LZ contributed equally to this work. BQ designed and revised the manuscript. HL performed the experiments and collected data. CL and LZ wrote the manuscript and analyzed data. XH and JD revised the manuscript.

Corresponding authors

Correspondence to Xiahong He, Jiahong Dong or Bin Qiu.

Ethics declarations

Ethics approval and consent to participate

All methods including the plant and soil material used in this study were carried out in accordance with relevant guidelines and regulations. Permission was obtained from the Knoxia growers for collecting K. valerianoides samples from the participating planting areas in this study with provisional annuity. Knoxia growers approved the manuscript to be published in peer reviewed journals.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, C., Zhang, L., Li, H. et al. Assessing the biodiversity of rhizosphere and endophytic fungi in Knoxia valerianoides under continuous cropping conditions. BMC Microbiol 24, 195 (2024). https://doi.org/10.1186/s12866-024-03357-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12866-024-03357-7

Keywords