Mengyuan He
Qianqian Wang3
Yiming Wang
Junhua Zhang
94% of researchers rate our articles as excellent or good
Learn more about the work of our research integrity team to safeguard the quality of each article we publish.
Find out more
- 1School of Life Sciences, Ningxia University, Yinchuan, China
- 2School of Ecology and Environment, Ningxia University, Yinchuan, China
- 3State Key Laboratory of Soil and Sustainable Agriculture, Institute of Soil Science, Chinese Academy of Sciences, Nanjing, China
Long-term monoculture of Lycium barbarum significantly affects its productivity and soil health. Soil microbiota, which mediate the sustainable development of soil ecosystems, are influenced by the age of wolfberry plants. However, the comprehensive effects of long-term cultivation of L. barbarum on the soil microbial community are not yet fully understood. Here, we assessed the effects of stand age on the diversity, composition, assembly, and symbiotic networks of bacterial and fungal communities in the root zone soil of L. barbarum using high-throughput sequencing technology. The results showed that stand age significantly affected the α-diversity of bacterial and fungal communities, as evidenced by the tendency of their Shannon and Chao1 indices to increase and then decrease. At the same time, the structure of soil bacterial and fungal communities was significantly influenced by tree age. However, Proteobacteria (28.77%–32.81%) was always the most dominant bacterial phylum, and Ascomycetes (49.72%–55.82%) was always the most dominant fungal phylum. A number of genus-level biomarkers were also identified in soils associated with roots of trees of varying ages. Additionally, stochastic processes dominated the assembly of soil bacterial communities, whereas the balance between stochastic and deterministic processes in the assembly of fungal communities fluctuated with stand age. The complexity and stability of bacterial and fungal community networks were notably affected by tree age, particularly in networks from 10- and 15-year-old trees. The partial least squares path modeling (PLS-PM) analysis emphasized that stand age can indirectly regulate the diversity and network complexity of both bacterial and fungal communities by influencing soil physicochemical properties. Furthermore, the bacterial community, but not the fungal community, exhibited direct and strong regulation of network complexity. The study offers valuable data for improving the soil quality and fruit yield of L. barbarum under long-term continuous cropping, which has implications for the sustainable development of the L. barbarum industry.
1 Introduction
Wolfberry (Lycium barbarum) is a deciduous shrub belonging to the Solanaceae family, primarily found in China and other parts of Asia (Liu et al., 2022). It demonstrates remarkable tolerance to various abiotic stresses, including high salinity, cold, and drought (Zhao et al., 2015). Due to its medicinal and nutritional values, the cultivation of wolfberry has emerged as a key agro-industry contributing to regional economic growth (Li et al., 2023). However, continuous cropping for more than a decade tends to decrease the fruit yield and quality of L. barbarum, posing a significant challenge not only to sustainable wolfberry cultivation and economic development but also to the sustainability of the soil ecosystem (Zhang et al., 2018). This decline is primarily attributed to the deterioration of soil quality resulting from inadequate soil management practices, continuous pesticide application, and long-term monocropping systems (Zhu et al., 2022).
Microorganisms are an essential component of the soil ecosystem and play an active role in maintaining plant and soil health (Zhang et al., 2024c). Increases in the years of L. barbarum cultivation may have significant effects on the diversity and function of soil microbial communities (Zhang et al., 2023a; He et al., 2024b). For instance, long-term continuous cropping can alter soil physicochemical properties, potentially leading to a reduction in overall microbial community diversity (Song et al., 2024). Continuous cultivation may also result in significant alterations in the microbial community composition (Na et al., 2016; Zhang et al., 2018). Na et al. (2019) explored the relationship between soil microbial communities and the stand age of L. barbarum, discovering that long-term monocropping of L. barbarum fosters the enrichment of key phytopathogenic fungal genera, such as Gibberella, Alternaria, and Stemphylium. These fungi potentially pose threats to the health of the older L. barbarum plants. A recent study further revealed that changes in bacterial community composition, influenced by stand age, can enhance the enrichment of Fusarium spp., thereby jeopardizing the soil health of L. barbarum fields (Na et al., 2021). Despite the increasing recognition of soil microbial community structure in L. barbarum fields, the mechanisms underlying the generation and maintenance of dynamic changes in soil microbial communities due to tree age remain inadequately understood and discussed.
The effect of stand age on soil microbial communities is not only in terms of diversity and community structure, but also in terms of microbial community assembly and interactions (Wang et al., 2022; Zhang et al., 2023b). The process of microbial community assembly is widely recognized in the field of microbial ecology (Nemergut et al., 2013). As previously described, stochastic and deterministic processes, which emphasize “born free” and “environmental filtering”, respectively, may be facilitated or constrained by environmental changes across time and space, thereby influencing microbial community assembly processes (Fargione et al., 2003; Zhou and Ning, 2017). Currently, the assembly of microbial communities has been studied primarily in cash crops such as corn, soybean, wheat, and rice (Liu et al., 2020a; Cao et al., 2024; Liu et al., 2024). How continuous cultivation affects the assembly of soil microbial communities in the root zone of L. barbarum, i.e., the balance between stochastic and deterministic processes, is not known. Furthermore, microbial network analysis has been utilized as an important tool to study microbial interactions. Although it does not fully simulate microbial interactions in real environments, it can aid in understanding microbial community responses to external environmental changes (Wang et al., 2023). The complexity and stability of soil microbial networks are typically influenced by a variety of factors (Fierer, 2017). Regarding the effect of L. barbarum stand age on microbial networks, only in the study by Na et al. (2021) was it observed that the complexity of young stand networks was higher than that of old stand networks. Although this study spanned multiple tree ages, it only compared networks from younger and older stands and lacked research on smaller scales, specifically comparisons between multiple networks of different tree ages. It is critical to address this, as it is essential for exploring the dynamics of the soil microbiota in L. barbarum fields and understanding long-term changes in soil quality.
In this study, based on high-throughput sequencing technology, we depicted the distribution patterns, successional patterns, assembly mechanisms, and ecological networks of soil bacterial and fungal communities in L. barbarum fields with different stand ages, and we explored the relationship between tree age, soil physicochemical properties, microbial community diversity, and network complexity. It was hypothesized that i) the richest soil microbial communities are distributed in mid-aged stands, which exhibit the highest crop productivity compared with young and old stands; ii) the assembly processes of the soil microbial community and their ecological networks are influenced by increasing stand ages; and iii) the diversity of soil bacterial and fungal communities influences changes in their network complexity. The results of this study could provide guidance for the long-term management of L. barbarum fields.
2 Materials and methods
2.1 Study site and soil sample
The study site is located in Zhongning County (105°26′–106°07′E, 36°49′–37°50′N), Zhongwei City, Ningxia Hui Autonomous Region, China. Zhongning is renowned for having the best fruit quality and the largest production area of L. barbarum in Ningxia; the annual average altitude, annual precipitation, and mean temperature at the study site are 1,370 m, 244 mm, and 10.5°C, respectively (Wang et al., 2019). The major soil type is Anthrosols, and the soil texture is sandy loam.
L. barbarum field plots of four stand ages (3, 6, 10, and 15 years) were selected. Fertilizers, mostly Stanley compound fertilizer (N–P2O5–K2O), were applied twice a year at an average rate of 6,000–6,500 kg/hm2. In June 2022, six L. barbarum plants of similar growth in each plot were selected, and soil was collected from within its crown width (i.e., the root zone of the L. barbarum) after removing litter such as dead leaves from the surface soil. Each sample was separated into two parts. The first part of the sample was used to determine soil physicochemical properties, while the second part was stored at −80°C until it was used for the analysis of microbial communities.
2.2 Soil property analyses
The soil physicochemical parameters (Table 1) were obtained using the standardized method of Lu (2000). The details are as follows: soil pH and total salt (TS) content were measured in suspensions at a soil/water ratio of 1:5 (w/v) using a pH meter (PHS-3C; INESA Scientific Instrument Co., Nanjing, China) and a salinity meter (VZ8371BZ; AZ Instrument Corp., Taipei, China). Soil water content (SWC) was determined using the drying method, i.e., the ratio of the difference in weight of the soil sample before and after drying by 105℃ to the weight of the soil sample before drying. Additionally, ammonium–nitrogen (NH4+–N) and nitrate–nitrogen (NO3−–N) were quantified using the distillation method and ultraviolet spectrophotometry, respectively. The determination of available phosphorus (AP) and potassium (AK) contents was based on the Olsen method and ammonium acetate extraction–flame photometry, respectively. Alkali-hydrolyzable nitrogen (AN) and organic matter (OM) were analyzed using the alkaline diffusion method and the potassium dichromate oxidation–external heating method, respectively.
Table 1
Table 1. Soil physicochemical properties in wolfberry orchards of different ages.
2.3 Soil DNA extraction and sequencing
According to the manufacturer’s instructions, total genomic DNA was extracted from 0.5-g soil samples using a FastDNA Spin Kit for Soil (MP Biomedicals, Santa Ana, CA, USA). The V4 variable region of the bacterial 16S ribosomal RNA genes was amplified by PCR using the primer pair of 515F (5′-GTGCCAGCMGCCGCGGTAA-3′) and 806R (5′-GGACTACNNGGGTATCTAAT-3′) (Ding et al., 2015). For fungi, the internal transcribed spacer region 1 (ITS1) was amplified by PCR using the primer pair of ITS5-1737F (5′-GGAAGTAAAAGTCGTAACAAGG-3′) and ITS2-2043R (5′-GCTGCGTTCTTCATCGATGC-3′) (Lu et al., 2013). The 30 μL PCR reactions contained 2× Phusion Master Mix (New England Biolabs, Ipswich, MA, USA), 1 μM of each primer, 10 μL template DNA (5–10 ng), and ultraclean water to make up the volume. PCR conditions were as follows: one cycle at 98°C for 1 min, followed by 30 cycles at 98°C for 10 s, 50°C for 30 s, and 72°C for 30 s, and a final extension step at 72°C for 5 min. All PCR products were checked using 2% agarose gel electrophoresis, and the target bands of 200–300 bp were recovered using a DNA Purification Kit (Tiangen, Beijing, China).
The concentration of purified PCR products was measured using NanoDrop 1000 (NanoDrop Technologies, Wilmington, DE, USA). The samples were equally mixed according to their concentrations to construct libraries. Sequencing was accomplished on an Illumina NovaSeq 6000 platform (https://www.novogene.com/). Furthermore, the raw reads were spliced and filtered to obtain clean data. Then, the sequences were identified and distinguished based on the barcode sequences to obtain valid data. The raw sequences of the microorganisms from this study have been submitted to the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI), and the accession numbers for bacteria and fungi are PRJNA1215385 and PRJNA1215412, respectively.
2.4 Assembly of microbial communities
Neutral community models (NCMs), used to quantify the importance of stochastic processes in the assembly of bacterial and fungal communities, were carried out using the “hmisc”, “minpack.lm”, and “stats4” packages. Furthermore, a 50% Modified Stochasticity Ratio (MST) was used to assess the order of importance of deterministic and stochastic processes in microbial community assembly. Based on the Levins model, the habitat niche width of bacteria and fungi, obtained from the “spaa” package, were used to further explain the mechanisms by which microbial community assembly occurred (Fang et al., 2023).
2.5 Construction and analysis of microbial ecological networks
Based on the Gephi software (version 0.9.7), four ecological networks of bacteria and fungi were constructed across various stand ages of L. barbarum. In the construction process, Operational Taxonomic Units (OTUs) that satisfied both an average relative abundance of less than 0.1% and a representation in less than 1/3 of the total sample size were removed in order to control the complexity of the network and improve its accuracy. Then, using the “igraph” and “psych” packages, Spearman’s correlation of the screened OTUs across various stand ages of L. barbarum was calculated. Only strong (correlation coefficient |r| > 0.70) and significant relationships (significance p < 0.01) were retained in the networks. The stability of the microbial networks was assessed by calculating their robustness, which evaluates the reliability of the remaining nodes after randomly removing 50% of the nodes in the network. The above analyses refer to the description in Zhang et al. (2024a). In addition, the subgraph function in the “igraph” package was used to extract the subnetworks from each network, and then the average degree of each subnetwork was computed, which serves as a measure of the network’s complexity (Zhang et al., 2024b).
2.6 Data statistical analysis
R language (version 4.2.3) and SPSS (version 22.0) were used to process and statistically analyze the data. For the bacterial and fungal community, the “Vegan” package was utilized to analyze their α-diversity metrics (Shannon and Chao1 indices), conduct the β-diversity, and perform the Adonis test (He et al., 2024a). The β-diversity of microbial communities, which reflects dissimilarities in microbial community composition across different treatments, was visualized using principal coordinates analysis (PCoA). The β-diversity decomposition was further employed to evaluate the sources of differences in microbial community composition, i.e., turnover or nestedness, based on the “betapart” package (Wu et al., 2024b). The linear discriminant analysis effect size (LEfSe) method was implemented to identify differences in bacterial and fungal communities at the top 15 genera across various stand ages using the “microeco” package. The thresholds for the Linear Discriminant Analysis (LDA) score and the p-value were set at 2 and 0.1, respectively (Zhang et al., 2024c).
The partial least squares path modeling (PLS-PM) was used to identify the relationships between observed and latent variables, specifically to assess the effects of stand age, physicochemical properties, and microbial communities on network complexity using the “plspm” package. In this model, forward selection was conducted to obtain the following best revealed variables: stand age, physicochemical properties (including OM, AN, NH4+–N, NO3−–N, AP, and AK), bacterial or fungal community (represented by the Shannon index), and network complexity (measured by edges and average degree). Stand age was classified as an exogenous variable, physicochemical properties and microbial communities were classified as endogenous variables, and network stability was classified as a response variable. It is worth noting that variables with low loadings (loading < 0.7) were excluded, and the predictive performance of the model was evaluated using the goodness of fit (GOF) index.
3 Results
3.1 Diversity of soil bacterial and fungal communities
The effect of long-term continuous cultivation on the diversity and composition of soil microbial communities in the root zone of L. barbarum is evident. The α-diversity of bacterial (Figures 1A, B) and fungal (Figures 1E, F) communities, characterized by the Shannon and Chao1 indices, exhibited a tendency to increase and then decrease with increasing tree age. Specifically, soil bacterial and fungal communities in the root zone of 10-year-old L. barbarum possessed the highest diversity and abundance compared to those of other tree ages. The results of PCoAs based on the Bray–Curtis distance revealed that the β-diversity of soil bacterial (Adonis: R2 = 0.336, p < 0.001) and fungal communities (Adonis: R2 = 0.239, p < 0.001) was significantly differentiated between tree ages, confirming the significant effect of stand age on the composition of soil bacterial and fungal communities in the root zone of L. barbarum (Figures 1C, G). Meanwhile, β-diversity decomposition indices for soil bacterial and fungal communities in the root zone of L. barbarum of different ages were calculated to ascertain the primary reasons for the differences in microbial communities (Figures 1D, H). The results indicated that differences in soil bacterial and fungal communities among tree ages were primarily attributed to species replacement (bacteria, 90.30%; fungi, 64.20%).
Figure 1
Figure 1. The α-diversity index of bacterial (A, B) and fungal communities (E, F) in the root zone soils of Lycium barbarum of different ages. The PCoA analysis of bacterial (C) and fungal (G) communities. The β-diversity decomposition indices of bacterial (D) and fungal (H) communities based on Bray–Curtis distance. PCoA, principal coordinates analysis. *** indicate significant differences between groups at significance levels of 0.001. The below is same.
3.2 Composition of soil bacterial and fungal communities
The composition of soil bacterial and fungal communities in the root zone of L. barbarum of different ages exhibited dramatic variations at the phylum level (Figures 2A, B). Proteobacteria (28.77%–32.81%) and Ascomycota (49.72%–55.82%) were predominantly dominant among the bacterial and fungal phyla, respectively. As the stand age of L. barbarum increased, the relative abundance of the bacterial phylum Actinobacteria showed an increasing trend, whereas that of Gemmatimonadota exhibited a decreasing trend. Compared to that of bacteria, the abundance of fungal phyla varied more markedly in the root zone soils of L. barbarum of different ages; for example, 15-year-old L. barbarum root zone soils had lower abundances of Mortierellomycota, Basidiomycota, Blastocladiomycota, Aphelidiomycota, Rozellomycota, and Zoopagomycota.
Figure 2
Figure 2. Composition of soil bacterial (A) and fungal (B) phyla in the root zone of four tree-aged wolfberry orchards. LEfSe analysis of soil bacterial (C) and fungal (D) genera in the root zone of four tree-aged wolfberry orchards. LEfSe, linear discriminant analysis effect size.
LEfSe analyses, which assessed differences in the abundance of microbial taxa at the genus level in the soil of the root zone of L. barbarum of different ages, provided results that showed the most representative bacterial and fungal information in the soil of each stand age (Figures 2C, D) . For bacteria, the genera Sphingomonas, Tychonema_CCAP_1459-11B, Bryobacter, Ellin6067, Gemmatimonas, Flavisolibacter, Rhodobacter, and Skermanella exhibited superior relative abundance in the root zone soil of 3-year-old L. barbarum plants compared to the other tree age groups. In contrast, the genera Pseudomonas, Pontibacter, Arthrobacter, and Rheinheimera were more prevalent in 15-year-old L. barbarum than in 3-, 6-, and 10-year-old L. barbarum. Half of these genera belong to the phylum Proteobacteria, while the others belong to Firmicutes, Bacteroidota, Actinobacteria, Acidobacteriota, Cyanobacteria, and Gemmatimonadetes. Similarly, for fungi, there were differences in the dominant genera enriched in the root zone soil of L. barbarum of different stand ages, which belonged to the phyla Ascomycota, Basidiomycota, Blastocladiomycota, and Chytridiomycota. In particular, more genus-level biomarkers were found in 10-year-old trees, such as Bipolaris, Mycena, Entoloma, Lactifluus, Catenophlyctis, Cladorrhinum, Spizellomyces, and Stanjemonium.
3.3 Microbial community assembly processes and ecological niche width
In this study, a neutral model and a modified stochastic ratio null model were utilized to assess the assembly processes of soil bacterial (Figure 3A) and fungal (Figure 3B) communities in the root zone of L. barbarum of different ages. The results of the neutral model revealed the dominance of stochastic processes in the construction of bacterial communities across the plantation chronosequence, with stochastic processes explaining 79.7%–82.8% of the variance. Bacterial taxa exhibited higher Nm values at 10 (Nm = 39,657) and 15 years of age (Nm = 40,153) compared to those at 3 (Nm = 36,532) and 6 years (Nm = 35,083), suggesting that species dispersal within bacterial communities is higher at 10 and 15 years of age than at 3 and 6 years. The assembly process of the fungal community seemed to fluctuate between different tree ages, and the contribution of stochastic processes accounted for 47.70%–58.40% of the variance. The Nm value of the fungal community exhibited an overall decreasing tendency with increasing tree age, indicating that the species dispersal within the fungal community was limited by the stand age to a certain extent. MST analyses further confirmed the significance of stochastic processes in the assembly of bacterial communities (Figure 3C). Conversely, the assembly of fungal communities (Figure 3D) exhibited fluctuations due to stand age, which were driven by a combination of stochastic and deterministic processes. Additionally, the ecological widths of bacteria (Figure 3E) and fungi (Figure 3F) displayed an overall trend of increasing and then decreasing. Specifically, they exhibited the highest and lowest ecological niche widths at 10 and 15 years of age, respectively. It is noteworthy that bacteria demonstrated a broader range of ecological widths compared to fungi.
Figure 3
Figure 3. The weight of stochastic processes in bacterial (A) and fungal (B) community assembly across various stand ages of wolfberry orchards, predicted using NCMs. The blue solid and dashed lines represent the best fit to the neutral model and the 95% confidence intervals of the model predictions, respectively. R2 represents the overall fit to the neutral model. Nm indicates the assessment of within-community estimated migration. MST was used for assessing the ratios of stochastic and deterministic processes in the assembly of bacterial (C) and fungal (D) communities across various stand ages. Ecological niche widths of bacterial (E) and fungal (F) communities. NCMs, neutral community models. ** and *** indicate significant differences between groups at significance levels of 0.01 and 0.001, respectively. The below is same.
3.4 Molecular ecological networks of soil bacterial and fungal communities
We constructed ecological networks for bacterial and fungal communities and calculated their topological parameters to assess species interactions in soils from the root zone of four wolfberry plants of different ages (Figure 4; Table 2). The results revealed that the complexity of the bacterial and fungal networks exhibited an increasing and then decreasing trend as stand age increased. For bacteria, the 15-year-old network was clearly distinct from the other age networks, with the lowest complexity, having up to 27.34% fewer connections and a 24.89% lower average degree than the other age networks. For fungi, the network complexity followed the order of 15 years < 3 years < 6 years < 10 years. The 10-year-old fungal network was the most complex, primarily because it had the highest number of links (3,293) and average degree (16.548) compared to the 3-, 9-, and 15-old networks (edges, 836, 1544, and 442, respectively; average degree, 6.481, 10.722, and 4.533, respectively). Interestingly, the ecological networks of both bacteria and fungi were dominated by positive associations. Furthermore, we quantitatively assessed the robustness of the microbial ecological network by randomly removing 50% of the nodes, which served as an indicator of the network stability. The results confirmed that the stability of both the 15-year-old bacterial and fungal networks was the lowest (bacteria, 0.279; fungi, 0.276). Notably, the highest network modularity values were observed in the 15-year-old bacterial and fungal networks, but not in those of the other stand ages. In addition, further analysis revealed that the bacterial phyla Proteobacteria, Unidentified bacteria, and Acidobacteriota, as well as the fungal phyla Ascomycota and Unidentified fungi, were the key nodes for maintaining the respective bacterial and fungal networks.
Figure 4
Figure 4. Ecological networks of bacterial (A) and fungal (B) communities in the root zone soil of Lycium barbarum orchards at different ages (3, 6, 10, and 15 years). Nodes in the network are colored according to different bacterial and fungal phyla. Positive and negative correlations between different nodes are shown by the pink and blue lines, respectively.
Table 2
Table 2. Topological parameters of bacterial and fungal networks in the root zone soil of Lycium barbarum of different ages.
3.5 Relationships among stand age, diversity of the microbial community, and network complexity
To determine the relationship between the diversity of soil microbial communities and their network complexity, correlation analyses were performed (Figures 5A, D). The results revealed a significant positive correlation between the Shannon (bacteria, p < 0.01; fungi, p < 0.01) and Chao1 (bacteria, p < 0.001; fungi, p < 0.001) indices of microbial communities and the complexity of their respective networks. This confirmed the close relationship between the diversity of soil microbial communities in the root zone of L. barbarum of different ages and the complexity of their networks. Further, PLS-PM analyses were conducted to assess the effects of stand age, physicochemical properties, and microbial community diversity on the complexity of the microbial network (Figures 5B, E). The results indicated that stand age had a direct and strong significant effect (p < 0.05) on the physicochemical properties of the soil, but not on the diversity of bacterial and fungal communities. However, soil physicochemical properties exhibited significant and direct adverse effects on the diversity of both bacterial (path coefficient, −0.72, p < 0.01) and fungal (path coefficient, −0.78, p < 0.001) communities. Furthermore, soil physicochemical properties also had a significant negative regulatory effect on the complexity of both bacterial (path coefficient, −0.47, p < 0.05) and fungal (path coefficient, −0.65, p < 0.01) networks. It is noteworthy that a significant and positive relationship persisted between the diversity of bacterial communities (path coefficient, 0.51, p < 0.01) and the complexity of their networks when considering the effects of stand age, physicochemical properties, and microbial community on network complexity. However, this was not the case for fungal communities (path coefficient, 0.26, p > 0.05). In summary, stand age can indirectly regulate the diversity of bacterial and fungal communities, as well as the complexity of their networks, by influencing soil physicochemical properties. Additionally, the diversity of bacterial communities, but not fungal communities, exhibited a strong and positive regulatory effect on their network complexity. The results of the standardized effects analysis (Figures 5C, F) revealed that soil physicochemical properties had a significant impact on the complexity of both bacterial and fungal communities, emerging as the strongest predictors of bacterial and fungal network complexity. Following closely were the stand age and diversity of bacterial and fungal communities, which served as the second and third most influential factors, respectively, emphasizing their notable impact on network complexity.
Figure 5
Figure 5. Correlation analysis between microbial community diversity and ecological network complexity (A, bacteria; D, fungi). PLS analysis reveals direct and indirect effects of tree age, physicochemical properties, and microbial communities on network complexity (B, bacteria; E, fungi). Standardized effects of various factors on the complexity of bacterial (C) and fungal (F) networks. PLS, partial least squares. * and ** indicate significant differences between groups at significance levels of 0.05 and 0.01, respectively. The below is same.
4 Discussion
4.1 The diversity and composition of soil microbial communities in the root zone of L. barbarum are affected by stand age
Long-term monoculture of L. barbarum can result in changes in the soil microbial community. Microorganisms, as crucial components of soil ecosystems, contribute to the decomposition of organic matter and nutrient cycling and play a pivotal role in the creation of healthy soil habitats (Liu et al., 2020a). The results of this study confirmed the highest α-diversity among bacterial and fungal communities in middle-aged L. barbarum stands and the lowest in old-aged stands, respectively, thus supporting our first hypothesis. This may be attributed to the availability of soil nutrients following continuous monoculture over an extended period, where the levels of soil nutrients (e.g., C, N, and K) dictate soil microbial growth and activity (Zhang et al., 2023b). Although mid-aged stands with vigorous growth may consume a large amount of soil nutrients, farmers often apply higher levels of fertilizers at this stage to ensure high crop yields. At the same time, during periods of vigorous plant growth, the incidence of plant apomixis also increases, which, to some extent, raises the level of soil nutrients and enhances the soil microenvironment (Gong et al., 2024). The abundance of soil nutrients is also beneficial for the propagation of microbial communities. Numerous studies have established a positive correlation between soil nutrient levels and the diversity of bacteria and fungi (Semenov et al., 2018). For example, OM, AK, and AP were found to be the main factors influencing the distribution of soil bacterial communities in the rhizosphere of cucumber in a study by Huang et al. (2023).
Consistent with previous findings, our results indicated that Proteobacteria and Ascomycetes were the most abundant phyla among soil bacterial and fungal communities, respectively, in the root zone of L. barbarum. This may be attributed to their strong adaptability and survivability, which allow them to be widely distributed in various types of agricultural soils (Song et al., 2024). Other bacterial and fungal taxa also exhibited quantitative advantages at different tree ages, particularly fungi, which are known for their role in decomposing soil organic matter (such as plant and animal residues, including dead leaves) (Barberán et al., 2014). For example, Mortierellomycota, serving as an indicator of soil organic matter and nutrient levels, decomposes soil organic compounds like sugars and is significantly enriched in root zone soils of mid-aged (6-year-old and 10-year-old) L. barbarum (Fu et al., 2020). In addition, LEfSe analyses screened for markers of bacterial and fungal genera across different tree ages, suggesting that some specific microbial taxa play their key roles only under certain specific conditions.
4.2 Soil microbial community assembly and network complexity were affected by the stand age of L. barbarum
Previous studies have shown that long-term continuous monoculture significantly impacts microbial community assembly processes, leading to alterations in the diversity and composition of soil microbial communities (Liu et al., 2020b). Our results indicated that, although stochastic processes predominantly governed the assembly of bacterial communities, the influence of stochastic processes in the soil bacterial community of old-aged L. barbarum was significantly lower than that observed in other tree ages. This reflects the influence of increasing tree age on the structural changes in microbial communities under long-term continuous cultivation conditions (Zhou et al., 2024). Compared to bacterial communities, stochastic processes play a notably less important role in fungal communities, emphasizing the significant impact of changes in the soil environment under long-term cultivation on fungal communities. Generally, the theory of “body size plasticity” suggests that larger organisms are more susceptible to external environmental factors than smaller organisms; hence, fungi, being typically larger than bacteria, are more susceptible to external environmental factors than bacteria (Farjalla et al., 2012). In terms of ecological niche width, fungi generally occupy a narrower niche than bacteria, implying that fungi have access to fewer resources than bacteria. This contributes to the dominance of stochastic processes (rather than deterministic processes) in bacterial community assembly, as bacteria can tolerate and exploit a wider range of environmental conditions (He et al., 2024a).
Analysis of microbial ecological networks offers valuable insights into the complexity and stability of microbial community assemblages in response to aging. It is generally accepted that complex and stable networks contribute to the potential for soil ecosystem services, which are crucial for maintaining sustainable soil ecosystems (Huo et al., 2023). Alterations in bacterial and fungal ecological networks underscore the impact on associations between bacterial and fungal species, aligning with our second hypothesis. A previous study by Na et al. (2021) also found that networks in younger stands had more nodes, edges, and higher average degrees than those in older stands, confirming the presence of negative feedbacks between plants and soil beyond a certain number of planting years. Differences in microbial ecological networks imply that the interactions between microbial taxa present in a given environment are different and to some extent represent the degree of connectedness between microbial taxa. Interestingly, a high proportion of positive associations exist in all these networks, implying that cooperation rather than competition dominated the process of assembling soil bacterial and fungal communities in the root zone of L. barbarum. Generally, positive associations symbolize positive interactions and symbiosis, which can enhance resistance to the external environment (Coyte et al., 2015).
4.3 Relationships between stand age, physicochemical property, microbial community, and network complexity
Based on PLS-PM analyses, we found that tree age primarily modulates microbial community diversity and network complexity by influencing physicochemical properties, thereby emphasizing the strong influence of soil nutrient content on soil bacterial and fungal communities. Similar conclusions were reached by Yi et al. (2023), who found that soil physicochemical properties are significantly and strongly influenced by region and play a crucial role in the succession of microbial communities. Soil physicochemical factors have been recognized as key determinants of microbial community diversity and composition (Zhou and Zhou, 2023). Furthermore, stand age may also influence the complexity of bacterial and fungal networks indirectly through physicochemical properties, suggesting that soil nutrient content may facilitate interactions and symbioses within soil bacterial and fungal communities (Wu et al., 2024a).
Regarding symbiotic networks of bacteria and fungi, previous studies have primarily focused on the effects of abiotic factors, such as tree age and physicochemical properties, on their complexity and stability (Wang et al., 2023; Niu et al., 2024). For instance, He et al. (2024a) discovered that salinity gradients were dominant in shaping bacterial symbiotic networks, while Yang et al. (2023) found that elevation gradients had a similar influence on fungal symbiotic networks. Our results revealed that in addition to abiotic factors, biotic factors—particularly the diversity of the bacterial community—exert an extremely significant and strong influence on the complexity of its network. This implies that the presence or absence of specific bacterial community members not only affects their interactions but also plays a pivotal role in the formation of the entire network (Li et al., 2020; Yu et al., 2024). Although the diversity of fungal communities exhibited significant differences across tree ages, there was no notable effect on the complexity of their networks. This may be attributed to the overlapping ecological niches and functional redundancy within the fungal communities (Zeng et al., 2023). In conclusion, higher microbial diversity and network complexity generally indicate more stable ecological processes, which are crucial for the stability and sustainability of agroecosystems.
5 Conclusion
This study provides a comprehensive understanding of the effects of long-term monoculture of L. barbarum on soil microbial community diversity, composition, assembly, and network properties. The results indicated that tree age exerts a significant effect on bacterial and fungal communities. The diversity of these communities increased with tree age but decreased in older L. barbarum stands, emphasizing the selective pressure exerted by tree age on microbial populations. The composition of bacterial and fungal communities varied with tree age, with some phyla and genera showing dominance in soils of different ages. For instance, LEfSe analysis identified several bacterial and fungal genera as biomarkers that differentiate tree ages, emphasizing their role in the root zone soils of L. barbarum at various ages. In addition, the complexity and stability of bacterial and fungal networks exhibited a tendency to increase and then decrease with stand age, highlighting the disruption caused by old-aged L. barbarum cultivation to the linkages within the soil microbial community. While the assembly of bacterial communities was consistently dominated by stochastic processes as tree age increased, the assembly of fungal communities appeared to be more variable, underscoring the constraints imposed by specific environments of different tree ages on fungal community assembly. Tree age emerged as a pivotal factor influencing microbial community diversity and network complexity, with these effects primarily mediated by physicochemical properties. The findings of this study imply that nutrient management during L. barbarum cultivation should be enhanced, and old, low-productivity L. barbarum plants should be replaced in a timely manner. This is crucial for maintaining the health of soil ecosystems in L. barbarum gardens.
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: NCBI BioProject, accession numbers PRJNA1215385 (for bacteria) and PRJNA1215412 (for fungi).
Author contributions
MH: Formal Analysis, Software, Visualization, Writing – original draft. QW: Formal Analysis, Resources, Software, Validation, Writing – original draft. YW: Validation, Data curation, Methodology, Project administration, Supervision, Writing – review & editing. JZ: Data curation, Project administration, Supervision, Writing – review & editing, Funding acquisition, Investigation.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Key R&D Program of China (grant number 2021YFD1900602), the Natural Science Foundation of Ningxia (2023AAC03141), and the Leading Talents in the Technological Innovation Program of Ningxia (2022GKLRLX02).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Barberán, A., Ramirez, K. S., Leff, J. W., Bradford, M. A., Wall, D. H., Fierer, N., et al. (2014). Why are some microbes more ubiquitous than others? Predicting habitat breadth Soil bacteria. Ecol. Lett. 17, 794–802. doi: 10.1111/ele.12282
Cao, J., Yang, X., Guo, M. Y., Wang, C. (2024). Reclamation of abandoned cropland switches fungal community assembly from deterministic to stochastic processes. Sci. Total Environ. 951, 175494. doi: 10.1016/j.scitotenv.2024.175494
Coyte, K. Z., Schluter, J., Foster, K. R. (2015). The ecology of the microbiome: Networks, competition, and stability. Science. 350, 663–666. doi: 10.1126/science.aad2602
Ding, J. J., Zhang, Y. G., Deng, Y., Cong, J., Lu, H., Sun, X., et al. (2015). Integrated metagenomics and network analysis of soil microbial community of the forest timberline. Sci. Rep. 5, 7994. doi: 10.1038/srep07994
Fang, W. K., Fan, T. Y., Wang, S., Yu, X. K., Lu, A., Wang, X. M., et al. (2023). Seasonal changes driving shifts in microbial community assembly and species coexistence in an urban river. Sci. Total Environ. 905, 167027. doi: 10.1016/j.scitotenv.2023.167027
Fargione, J., Brown, C. S., Tilman, D. (2003). Community assembly and invasion: An experimental test of neutral versus niche processes. Proc. Natl. Acad. Sci. 100, 8916–8920. doi: 10.1073/pnas.1033107100
Farjalla, V. F., Srivastava, D. S., Marino, N. A. C., Azevedo, F. D., Dib, V., Lopes, P. M., et al. (2012). Ecological determinism increases with organism size. Ecology. 93, 1752–1759. doi: 10.1890/11-1144.1
Fierer, N. (2017). Embracing the unknown: disentangling the complexities of the soil microbiome. Nat. Rev. Microbiol. 15, 579–590. doi: 10.1038/nrmicro.2017.87
Fu, H. P., Zhou, P. Q., Wang, Y. J., Mo, Z. D., Li, Z. H., Ma, Q. P., et al. (2020). Effects of intercropping different green manures on fungal community characteristics in rhizosphere soil of tea plant. J. Tea Communication. 47, 406–415.
Gong, X. Q., Jarvie, S., Wen, J., Su, N., Yan, Y. Z., Liu., Q. F., et al. (2024). Compared with soil fungal diversity and microbial network complexity, soil bacterial diversity drives soil multifunctionality during the restoration process. J. Environ. Manage. 354, 120379. doi: 10.1016/j.jenvman.2024.120379
He, M. Y., Shen, C., Peng, S., Wang, Y. D., Sun, J. B., Zhang, J. H., et al. (2024a). The influence of soil salinization, induced by the backwater effect of the Yellow River, on microbial community dynamics and ecosystem functioning in arid regions. Environ. Res. 262, 119854. doi: 10.1016/j.envres.2024.119854
He, M. Y., Shen, C., Zhang, J. H., Wang, Y. D. (2024b). Effects of continuous cropping on the physiochemical properties, pesticide residues, and microbial community in the root soil of Lycium barbarum. Environ. Science. 45, 5578–5590. doi: 10.13227/j.hjkx.202311078
Huang, S. C., Yu, J. H., Hou, D., Yue, H. Z., Zhang, D. Q., Li, Y. J., et al. (2023). Response of soil microbial community diversity to continuous cucumber cropping in facilities along the Yellow River irrigation area. PloS One 18, e0289772. doi: 10.1371/journal.pone.0289772
Huo, X. Y., Ren, C. J., Wang, D. X., Wu, R. Q., Wang, Y. S., Li, Z. F., et al. (2023). Microbial community assembly and its influencing factors of secondary forests in Qinling Mountains. Soil Biol. Biochem. 184, 109075. doi: 10.1016/j.soilbio.2023.109075
Li, J., Li, C., Kou, Y., Yao, M. J., He, Z. L., Li, X. Z., et al. (2020). Distinct mechanisms shape soil bacterial and fungal co-occurrence networks in a mountain ecosystem. FEMS Microbiol. Ecol. 96, fiaa030. doi: 10.1093/femsec/fiaa030
Li, Y. K., Zou, N., Liang, X. J., Zhou, X., Guo, S. H., Wang, Y. J., et al. (2023). Effects of nitrogen input on soil bacterial community structure and soil nitrogen cycling in the rhizosphere soil of Lycium barbarum L. Front. Microbiol. 13. doi: 10.3389/fmicb.2022.1070817
Liu, W. J., Graham, E. B., Zhong, L. H., Zhang, J. W., Li, S. J., Lin, X. G., et al. (2020a). Long-term stochasticity combines with short-term variability in assembly processes to underlie rice paddy sustainability. Front. Microbiol. 11. doi: 10.3389/fmicb.2020.00873
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
Liu, X., Wang, Y. Z., Liu, Y. H., Chen, H., Hu, Y. L. (2020b). Response of bacterial and fungal soil communities to Chinese Fir (Cunninghamia lanceolate) long-term monoculture plantations. Front. Microbiol. 11. doi: 10.3389/fmicb.2020.00181
Liu, Y., Yang, C. D., Fu, X., Zhao, F. Z., Wang, J. (2024). Mulching drive changes in soil microbial community assembly processes and networks across aggregate fractions. Eur. J. Soil. Biol. 122, 103664. doi: 10.1016/j.ejsobi.2024.103664
Lu, L. H., Yin, S. X., Liu, X., Zhang, W. M., Gu, T. Y., Shen, Q. R., et al. (2013). Fungal networks in yield-invigorating and -debilitating soils induced by prolonged potato monoculture. Soil Biol. Biochem. 65, 186–194. doi: 10.1016/j.soilbio.2013.05.025
Na, X. F., Ma, C. X., Ma, S. L., Ma, X. R., Zhu, X. M., Xu, P. X., et al. (2019). Monocropping decouples plant–bacteria interaction and strengthens phytopathogenic fungi colonization in the rhizosphere of a perennial plant species. Plant Soil. 445, 549–564. doi: 10.1007/s11104-019-04311-7
Na, X. F., Ma, S. L., Ma, C. X., Liu, Z. Y., Xu, P. X., Zhu, H. B., et al. (2021). Lycium barbarum L. (goji berry) monocropping causes microbial diversity loss and induces Fusarium spp. enrichment at distinct soil layers. Appl. Soil Ecol. 168, 104107. doi: 10.1016/j.apsoil.2021.104107
Na, X. F., Zheng, G. Q., Peng, L., Lei, C. Y., Yang, H. Y., Ma, Y., et al. (2016). Microbial biodiversity in rhizosphere of Lycium bararum L. relative to cultivation history. Acta Pedologica Sinica. 53, 241–252. doi: 10.11766/trxb201503030643
Nemergut, D. R., Schmidt, S. K., Fukami, T., O'Neill, S. P., Bilinski, T. M., Stanish, L. F., et al. (2013). Patterns and processes of microbial community assembly. Microbiol. Mol. Biol. Rev. 77, 342–356. doi: 10.1128/mmbr.00051-12
Niu, Y. C., Kang, E. Z., Li, Y., Zhang, X. D., Yan, Z. Q., Li, M., et al. (2024). Non-flooding conditions caused by water table drawdown alter microbial network complexity and decrease multifunctionality in alpine wetland soils. Environ. Res. 254, 119152. doi: 10.1016/j.envres.2024.119152
Semenov, M. V., Chernov, T., Tkhakakhova, A. K., Zhelezova, A. D., Ivanova, E. A., Kolganova, T. V., et al. (2018). Distribution of prokaryotic communities throughout the Chernozem profiles under different land uses for over a century. Appl. Soil Ecol. 127, 8–18. doi: 10.1016/j.apsoil.2018.03.002
Song, R., Li, Y. S., Zhu, Z. W., Zhang, L., Wang, H., Li, H. (2024). Vineyard reclamation alters soil properties and microbial community in desertified land. Catena. 246, 108399. doi: 10.1016/j.catena.2024.108399
Wang, Y. J., Liang, X. J., Guo, S. J., Li, Y. K., Zhang, B., Yin, Y., et al. (2019). Evaluation of nutrients and related environmental factors for wolfberry (Lycium barbarum) fruits grown in the different areas of China. Biochem. Syst. Ecol. 86, 103916. doi: 10.1016/j.bse.2019.103916
Wang, K., Wang, X. X., Fei, H. Y., Wan, C. Y., Han, F. P. (2022). Changes in diversity, composition and assembly processes of soil microbial communities during Robinia pseudoacacia L. restoration on the Loess Plateau. China. J. Arid Land. 14, 561–575. doi: 10.1007/s40333-022-0064-2
Wang, X., Zhang, Q., Zhang, Z. J., Li, W. J., Liu, W. C., Xiao, N. J., et al. (2023). Decreased soil multifunctionality is associated with altered microbial network properties under precipitation reduction in a semiarid grassland. iMeta. 2, e106. doi: 10.1002/imt2.106
Wu, L. K., Ren, C. Y., Jiang, H., Zhang, W. Y., Chen, N., Zhao, X. N., et al. (2024a). Land abandonment transforms soil microbiome stability and functional profiles in apple orchards of the Chinese Losses Plateau. Sci. Total Environ. 906, 167556. doi: 10.1016/j.scitotenv.2023.167556
Wu, Y. H., Yang, Z. Y., Chen, S. R., Chen, L. Y., Sui, M. Z., Zhang, G. Q., et al. (2024b). Drivers of β-diversity and its component-dependence along the natural restoration in an extremely heterogeneous forest ecosystem. Catena. 243, 108224. doi: 10.1016/j.catena.2024.108224
Yang, Y., Qiu, K. Y., Xie, Y. Z., Li, X. C., Zhang, S., Liu, W. S., et al. (2023). Geographical, climatic, and soil factors control the altitudinal pattern of rhizosphere microbial diversity and its driving effect on root zone soil multifunctionality in mountain ecosystems. Sci. Total Environ. 904, 166932. doi: 10.1016/j.scitotenv.2023.166932
Yi, K., Li, Z. Q., Shang, D. S., Zhang, C. G., Li, M. L., Lin, D. Z., et al. (2023). Comparison of the response of microbial communities to region and rootstock disease differences in tobacco soils of southwestern China. Front. Microbiol. 14. doi: 10.3389/fmicb.2023.1333877
Yu, Z., Lee, C., Kerfahi, D., Li, N., Yamamoto, N., Yang, T., et al. (2024). Elevational dynamics in soil microbial co-occurrence: Disentangling biotic and abiotic influences on bacterial and fungal networks on Mt. Seorak. Soil Ecol. Lett. 6, 240246. doi: 10.1007/s42832-024-0246-2
Zeng, Q. C., Lebreton, A., Auer, L., Man, X. W., Jia, L. K., Wang, G. S., et al. (2023). Stable functional structure despite high taxonomic variability across fungal communities in soils of old-growth montane forests. Microbiome. 11, 217. doi: 10.1186/s40168-023-01650-7
Zhang, J. H., He, M. Y., Jia, K. L., Shen, C. (2023a). Effects of farmland abandonment on anthropogenic-alluvial soil microbiota and contaminant residues in Lycium barbarum fields. J. Appl. Microbiol. 134, lxad198. doi: 10.1093/jambio/lxad198
Zhang, C., Lei, S. L., Wu, H. Y., Liao, L. R., Wang, X. T., Zhang, L., et al. (2024a). Simplified microbial network reduced microbial structure stability and soil functionality in alpine grassland along a natural aridity gradient. Soil Biol. Biochem. 191, 109366. doi: 10.1016/j.soilbio.2024.109366
Zhang, J. H., Li, M., Jia, K. L., Zheng, G. Q., Long, X. E. (2018). Seasonal variation rather than stand age determines bacterial diversity in the rhizosphere of wolfberry (Lycium barbarum L.) associated with soil degradation. J. Soils Sed. 18, 1518–1529. doi: 10.1007/s11368-017-1854-6
Zhang, Y. K., Li, M., Zhang, J., Cui, J., Wang, X. Y., Song, M. M., et al. (2024c). Spatial soil heterogeneity rather than the invasion of Spartina alterniflora drives soil bacterial community assembly in an Eastern Chinese intertidal zone along an estuary coastline. Catena. 237, 107784. doi: 10.1016/j.catena.2023.107784
Zhang, C. Z., Xiang, X. J., Yang, T., Liu, X., Ma, Y. Y., Zhang, K. P., et al. (2024b). Nitrogen fertilization reduces plant diversity by changing the diversity and stability of arbuscular mycorrhizal fungal community in a temperate steppe. Sci. Total Environ. 918, 170775. doi: 10.1016/j.scitotenv.2024.170775
Zhang, X. Y., Xiong, S.-Y., Wu, X. K., Zeng, B. B., Mo, Y. M., Deng, Z. C., et al. (2023b). Dynamics of microbialcommunity structure, function and assembly mechanism with increasing stand age of Slash Pine (Pinus elliottii) Plantations in Houtian Sandy Area, South China. J. Microbiol. 61, 953–966. doi: 10.1007/s12275-023-00089-7
Zhao, J. H., Li, H. X., Xi, W. P., Wei, A., Niu, L. L., Cao, Y. L., et al. (2015). Changes in sugars and organic acids in wolfberry (Lycium barbarum L.) fruit during development and maturation. Food Chem. 173, 718–724. doi: 10.1016/j.foodchem.2014.10.082
Zhou, C. F., Gao, Q. Q., Tigabu, M., Wang, S. Z., Cao, S., Yu, Y. C. (2024). Continuous planting of Chinese fir monocultures significantly influences dissolved organic matter content and microbial assembly processes. Sci. Total Environ. 926, 171943. doi: 10.1016/j.scitotenv.2024.171943
Zhou, J., Ning, D. (2017). Stochastic community assembly: Does it matter in microbial ecology? Microbiol. Mol. Biol. Rev. 81, e00002–e00017. doi: 10.1128/mmbr.00002-17
Zhou, Y. J., Zhou, S. L. (2023). Role of microplastics in microbial community structure and functions in urban soils. J. Hazard. Mater. 459, 132141. doi: 10.1016/j.jhazmat.2023.132141
Keywords: Lycium barbarum, stand age, soil microorganism, ecological network, assembly mechanisms
Citation: He M, Wang Q, Wang Y and Zhang J (2025) Temporal dynamics of soil microbial symbioses in the root zone of wolfberry: deciphering the effects of biotic and abiotic factors on bacterial and fungal ecological networks. Front. Plant Sci. 16:1518439. doi: 10.3389/fpls.2025.1518439
Received: 31 October 2024; Accepted: 10 February 2025;
Published: 10 March 2025.
Edited by:
Amjad Iqbal, Abdul Wali Khan University Mardan, PakistanReviewed by:
Muhammad Irshad, Abdul Wali Khan University Mardan, PakistanXiaorong Wei, Northwest A&F University, China
Copyright © 2025 He, Wang, Wang and Zhang. 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: Junhua Zhang, emhhbmdqdW5odWE3MjhAMTYzLmNvbQ==