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
Xiaoxia Li1,2,3
Muhammad Awais1,2,3
Shuang Wang1,2,3,4Zhu Zhang1,2,3Shuning Zhao1,2,3
Yufeng Liu1,2,3Zhouping Sun1,2,3*
Hongdan Fu1,2,3*Tianlai Li1,2,3- 1College of Horticulture, Shenyang Agricultural University, Shenyang, China
- 2Key Laboratory of Protected Horticulture of Education Ministry and Liaoning Province, Shenyang, China
- 3National and Local Joint Engineering Research Center of Northern Horticultural Facilities Design and Application Technology, Shenyang, China
- 4College of Agriculture, Eastern Liaoning University, Dandong, Liaoning, China
Interactions between endophytes (endophytic bacteria and fungi) and plants are crucial in maintaining crop fitness in agricultural systems, particularly in relation to abundant and rare subcommunities involved in community construction. However, the influence of long-term fertilization on heterogeneous rhizosphere nitrogen and phosphorus environments and how these conditions affect the key subcommunities of root endophytes and their community assembly mechanisms remain unclear. We studied the 26th year of a field experiment conducted in a greenhouse with varying levels of nitrogen and phosphorus (CKP0, CKP1, CNP0, CNP1, ONP0, and ONP1) to assess the composition of tomato root endophytes and their impact on yield. We employed 16S rRNA and fungal ITS region amplicon sequencing to investigate the assembly mechanisms of abundant and rare endophytic subcommunities, network correlations, core subcommunity structures, and key species that enhance crop yield. The results indicated that organic manure and phosphorus fertilizers significantly increased the rhizosphere soil nitrogen content, phosphorus content, and phosphorus availability (labile P, moderately labile P, and non-labile P). These fertilizers also significantly affected the composition (based on Bray-Curtis distance) and community assembly processes (βNTI) of endophytic microbial subcommunities. The assembly of both bacterial and fungal subcommunities was primarily governed by dispersal limitation, with community structures being significantly regulated by the content of rhizosphere soil available nitrogen (AN) and moderately labile P (MLP). Rare bacterial and fungal subcommunities complemented the ecological niches of abundant subcommunities in the co-occurrence network, supporting community functions and enhancing network stability. Nitrogen-adapting abundant and rare bacterial subcommunities provided a stronger predictive correlation for tomato yield than phosphorus-adapting fungal subcommunities. Additionally, three core genera of rare endophytic bacteria such as Arthrobacter, Microbacterium, and Sphingobium were identified as potentially involved in improving crop yield improvement. These findings revealed the distinct assembly mechanisms of endophytic microbial subcommunities affected by fertilization, enhancing our understanding of better management practices and controlling endophytes to improve crop yield in intensive agricultural ecosystems.
Graphical abstract
Graphical Abstract.
1 Introduction
Endophytes are a group of living microorganisms that can successfully colonize various plant tissues without causing disease symptoms (Adeleke and Babalola, 2021; Santoyo et al., 2016). Plants develop complex coexistence and mutualistic relationships with various endophytes (Battu and Ulaganathan, 2020; Harman and Uphoff, 2019) due to the essential functions provided by these endophytes, such as nutrient acquisition (Bolan, 1991; Burragoni and Jeon, 2021; Xie et al., 2024), stress tolerance (Ali et al., 2014; Battu and Ulaganathan, 2020; Harman and Uphoff, 2019; Zhang et al., 2024), and disease suppression (Adeleke et al., 2021; Card et al., 2016; Omomowo and Babalola, 2019; Sena et al., 2024). The mutualistic benefits between plants and endophytes also depend on the genetic and environmental conditions of the host (Meng et al., 2022). In particular, plant root endophytes are inherited through host genes and successfully colonize after filtrating through rhizosphere microorganisms (Meng et al., 2022; Xie et al., 2024). Hence, plant root endophytes exhibit survival strategies similar to those of rhizosphere soil microorganisms, offering significant benefits for healthy plant growth and improved yield (Guo et al., 2024; Harman and Uphoff, 2019). Root endophyte symbiosis, which improves plant adaptability under stress conditions, is increasingly important. Notably, some root endophytes have developed strategies to adapt to the unique soil dynamic environments (Beltran-Garcia et al., 2021; Byregowda et al., 2022; Collavino et al., 2020; Edwards et al., 2015). For instance, root endophytic fungi play fundamental ecological roles in substituting phosphorus strategies in non-mycorrhizal crops or nitrogen availability (Mehta et al., 2019; Sindhu et al., 2024; Sun et al., 2022; Xie et al., 2024). These root endophytes served as a more accurate and effective microbial indicators for evaluating plant resistance to environmental disturbances (Rani et al., 2022; Sun et al., 2021). Although the understanding of microbial community composition in different soil environments has rapidly increased, it remains unclear whether the resistance of root endophytic bacteria and fungi to soil environments influences their beneficial effects on crop yield.
In the natural ecosystem, local microbial communities generally comprise a few dominant species (abundant taxa) and numerous low-abundance species (rare taxa) (Ma et al., 2024; Zheng et al., 2021). The assembly dynamics of abundant and rare taxa in different ecosystems are complex and inconsistent (Gschwend et al., 2021), especially in the context of soil nutrient imbalance in agricultural ecosystems, where they exhibit different distribution patterns and functional characteristics (Jiao et al., 2019a; Ma et al., 2022a; Wu et al., 2023). These differences are attributed to varying microbial survival strategies and substrate preferences (Long and Yao, 2020; Wan et al., 2021). Understanding community assembly process is vital for elucidating the coexistence status among soil microorganisms and preserving species diversity (Chen et al., 2020; Liu et al., 2024). The community assembly process has been categorized into deterministic process, including homogeneous and variable selection, and stochastic process, including homogenizing dispersal, dispersal limitation, and undominated (Chase and Myers, 2011; Stegen et al., 2015). Due to diverse ecological functions and collaborative capabilities of species, the assembly of abundant taxa communities often shows greater stochasticity (Xiong et al., 2020), which underpins the overall stability of the microbial community (Wan et al., 2021). Rare taxa, with high genetic diversity, can adapt to specific environments and play key roles in stabilizing microbial communities and ecosystem functions (Cui et al., 2023; Jia et al., 2022; Ma et al., 2024; Xue et al., 2020; Yu et al., 2024; Zhang et al., 2023). Some key taxa variations of the rare subcommunity have shown ecological potential for mediating microbial interactions, dominating community stability, and contributing to the complexity of symbiotic networks (Jousset et al., 2017; Liang et al., 2020; Xiong et al., 2020). Changes in the proportion and diversity of abundant and rare endophytic subcommunities inevitably affect the stability and ecological function of endophytic communities (Xi et al., 2024; Xiong et al., 2020; Zhang et al., 2024). Thus, understanding the distribution characteristics, assembly processes, and symbiosis patterns of abundant and rare subcommunities is crucial for examining how long-term fertilization affects endophytic microbial diversity and ecological stability.
The use of excessive inorganic nitrogen and phosphorus fertilizers in greenhouse cultivation, driven by yield and economic benefit, is a widespread practice (Wang et al., 2022; Zhang et al., 2017). This practice has led to excessive nutrient accumulation in the soil (Han et al., 2015; Tian et al., 2022) and exacerbated the negative effects on soil microbial community assembly (Bai et al., 2022; Banerjee et al., 2019; Liu et al., 2023b). Traditionally, manure has helped mitigate the harmful effects of chemical fertilizers on microbial communities, thereby contributing to the yield (Zhang et al., 2021). This positive and effective agronomic practice may be key to determining facility crop yield and maintaining agricultural sustainability. Increasing evidence suggests that environmental changes, such as variations in soil nitrogen and phosphorus levels, play a significant role in driving niche variation of root endophytic microbial species and determining the assembly process of root endophytic microbial communities (Bai et al., 2022; Long and Yao, 2020; Sindhu et al., 2024). Appropriate nitrogen levels are conducive to the dominant characteristics of endophytic fungi and endophytic bacterial communities, as well as higher diversity and richness (Liu et al., 2023b; Miranda-Carrazco et al., 2022; Sindhu et al., 2024). The composition and structural changes of rice root endophytic bacterial communities are affected by chemical phosphorus fertilizer (Long and Yao, 2020). Differences in the community structure of root endophytic bacteria and fungi in response to fertilization may be closely related to the host’s nutrient utilization. Therefore, it is essential to identify external environmental factors, particularly nitrogen and phosphorus availability, that influence the assembly of abundant and rare root endophytic subcommunities and their contribution to the yield. Understanding these factors can shed light on the co-evolution of different niche endophytic subcommunities with plants and help develop targeted strategies to manipulate key subcommunities and species to improve crop yields. In this study, we conducted a long-term experiment with nitrogen, phosphorus, and organic manure application in a solar greenhouse. The aims of this study were to (i) identify the relative contributions of the environmental factors in shaping root endophytes and assess how deterministic and stochastic processes influence abundant and rare subcommunities; and (ii) explore the species traits, ecological functions, and co-occurrence patterns of abundant and rare endophytic microbial taxa. We hypothesized the following: (i) The availability of nitrogen and phosphorus in the rhizosphere soil are potential key environmental variables driving stochastic process (dispersal limitation) in endophytic bacterial and fungal assembly, causing structural changes in subcommunities. (ii) Phosphorus-adapting endophytic fungi may occupy niches that complement nitrogen-adaption endophytic bacterial subcommunities, each contributing to main ecological functions. (iii) Rare bacterial subcommunities provide core species that maintain the stability of the endophytic microbial network and complement the functions of abundant bacterial subcommunities, ultimately enhancing crop yield potential.
2 Materials and methods
2.1 Field site and experimental design
The long-term targeted fertilization experiment was designed in 1996 as a solar greenhouse field at the Shenyang Agricultural University of Liaoning Province, China (41°480 N, 123°250 E). The experimental soil was classified as Glesol (FAO and ISRIC, 1988). The initial basic chemical properties were as follows: pH = 6.75, soil organic matter (SOM) = 24.30 g·kg−1, total P (TP) = 1.37 g·kg−1, available P (AP) = 70.80 mg·kg−1, labile P (LP) = 333.90 mg·kg−1, moderately labile P (MLP) = 364.73 mg·kg−1, non-labile P (NLP) = 1163.38 mg·kg−1, and total N (TN) = 1.16 g·kg−1. The fertilization experiment site has a cropping system of continuous tomato cultivation in spring and autumn every year. The fertilization experiment consisted of six treatments arranged in a completely randomized block design (1.5 * 1.0 * 0.8 m3 for each pool): (1) CKP0 (unfertilized), (2) CKP1 (chemical phosphorus fertilizer), (3) CNP0 (chemical nitrogen fertilizer), (4) CNP1 (chemical nitrogen plus phosphorus), (5) ONP0 (organic manure), and (6) ONP1 (organic manure plus phosphorus). Each treatment was replicated in three independent cement pools. The fertilizers application rate per year are shown in Supplementary Table S1. Chemical nitrogen fertilizer was applied twice during the tomato whole growing season. Chemical phosphorus fertilizers and organic manure (horse manure obtained 8.69% organic matter, 0.45% total N, and 0.30% total phosphorus) were evenly spread on the surface and incorporated into the plow layer by manual hoeing before planting. The fertilization experiment was conducted in the autumn season of the 26th year from August 2022 to December 2022. Each tomato (“Meisheng”) plant had a single branch after pruning, and there were four clusters of fruits on the branch and five fruits per cluster.
2.2 Root and rhizosphere soil sampling
Four rhizosphere soils were randomly collected from each pool and then merged into one rhizosphere soil sample (Bulgarelli et al., 2012). Four tomato root samples were randomly collected from each pool, washed, and merged into one root sample, which was stored in sterilized 1 × PBS (pH = 7.0). All samples were promptly transported to the laboratory. The rhizosphere soil samples were stored at −80°C, while root samples for processed further to prepare for the isolation of root endophytic fungi and bacteria. A total of 2.0 g mixed roots of each sample were surface sterilized in 75% ethanol (1 min), 3.25% sodium hypochlorite (3 min), and 75% ethanol (30 s), followed by three rinses in distilled water (Guo et al., 2000). Root samples were then stored in liquid nitrogen immediately. The key difference was that all solutions used were sterile and the last sterile washed water of each root sample was placed in a PDA medium for 1 week at 28°C to confirm the effectiveness of surface sterilization. A total of 36 experimental samples, including 18 rhizosphere soil samples and 18 root samples (6 treatments × 3 duplicates), were stored at −80°C for subsequent DNA extraction.
2.3 Soil chemical parameters determination
Soil pH was measured using 1:2.5 (w/v) soil-to-water suspension. Soil organic carbon (SOC) was determined using the H2SO4-K2Cr2O7 oxidation capacity method (Mebius, 1960). Available phosphorus (AP) was determined through 0.5 M NaHCO3 extraction, while total phosphorus (TP) was measured through H2SO4-HClO4 digestion. Soil phosphorus fractions were prepared and analyzed as reported by Hedley et al. (1982) and Tiessen and Moir (1993). These phosphorus fractions were classified into LP (resin-Pi, NaHCO3-Pi, and NaHCO3-Po), MLP (NaOH-Pi and NaOH-Po), and NLP (HCl-Pi, concentrated HCl-Pi, concentrated HCl-Po, and residual P) (Maranguit et al., 2017; Shi et al., 2023). Total nitrogen (TN) was quantified using an automatic Kjeldahl distillation–titration method (Cabrera and Beare, 1993). Available nitrogen (AN) was measured using the alkali-hydrolyzed diffusion method, while soil nitrate (NO3−-N) and ammonium (NH4+-N) were quantified using a flow-injection autoanalyzer (Skalar San++ CFA, Erkelenz, Germany) (Bornemisza et al., 1988).
2.4 DNA extraction and high-throughput sequencing
Total DNA was isolated from 0.5 g of pooled root samples (which were created by randomly collecting four surface-sterilized roots and mixing them into a single sample from each pool) using the CretMagTM Plant DNA Mini Kit: CretBiotech, Suzhou, China, following the manufacturer’s instructions. The DNA samples obtained from root samples were amplified using PCR with specific primer sets that included adaptors and barcodes. Fungal ITS1 region amplification was as follows: ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and ITS2 (5′- GCTGCGTTCTTCATCGATGC-3′). 16S rDNA gene V5-V7 region amplification was as follows: 799F (5′- AACMGGATTAGATACCCKG-3′) and 1193R (5′-ACGTCATCCCCACCTTCC-3′) (Liu et al., 2023b). Additionally, to ensure efficiency and accuracy, PCR amplification was performed using TaKaRa Premix Taq® Version 2.0 (TaKaRa Biotechnology Co., Dalian, China), and the reactions took place in the thermocycler PCR system BioRad S1000 (Bio-Rad Laboratory, CA) under the following thermal cycling program for fungi: 95°C/3 min (initial denaturation) followed by 34 cycles of 95°C/20 s, 56°C/20 s, 72°C/30 s, 72°C/5 min (final extension), and lastly held at 12°C. The following thermal cycling program for bacteria is as follows: 98°C/30 s (initial denaturation) followed by 32 cycles of 98°C/10 s, 53°C/20 s, 72°C/30 s, 72°C/2 min (final extension), and lastly hold at 12°C. The PCR mixtures for fungi contain 25 μL of 2 × ES Taq MasterMix(Dye), 2 μL of forward primer (10 μM), 2 μL of reverse primer (10 μM), 50 ng of template DNA, and ddH2O added to a final volume of up to 50 μL. The PCR mixtures for bacteria contain 25 μL of Q5, 5 μL of forward primer (10 μM), 5 μL of reverse primer (10 μM), 50 ng of template DNA, and ddH2O added to a final volume of up to 50 μL. PCR reactions were performed in triplicate. DNA integrity of isolated DNA was visually inspected using 1.5% agarose gel electrophoresis, and DNA concentration and volume were determined using the NanoDrop2000 (Thermo Fisher Scientific, United States) (For detailed information, refer Supplementary Table S2). During the PCR process, negative controls with water and positive control tests were set. After comparing the concentration of PCR products, the required volume of each sample was calculated based on the principle of equal mass, and the PCR products were mixed accordingly. Target strips were cut and purified using the E.Z.N.A.® Gel Extraction Kit (Omega, United States) and quantified using GeneTools Analysis Software (Version4.03.05.0, SynGene), following the manufacturer’s instructions. The cDNA library was constructed using the NEBNext® Ultra™ II DNA Library Prep Kit for Illumina® (New England Biolabs, United States), according to the manufacturer’s instruction. Purified amplicons were pooled in equimolar and paired-end sequenced by the Shanghai Biotree Biotech Co., Ltd. (Shanghai, China) on the Illumina HiSeq PE250 sequencing platform. The resulting paired sequence reads were merged, trimmed, filtered, and clustered into zero-radius operational taxonomic units using USEARCH (version 0.2.7) with a sequence similarity threshold of 0.97 (Edgar, 2016a; Rognes et al., 2016). The UNITE (ITS) database and the Ribosomal Database Project (RDP) (16S) database were used for taxonomical alignment at a threshold of 0.8 using the SINTAX algorithm (Edgar, 2016b). These databases were used to align multiple sequences, remove annotation for chloroplasts or mitochondria, analyze the phylogenetic relationship between different operational taxonomic units (OTUs), and investigate the differences between the dominant microbial community species across different groups. All of raw reads data are available in the NCBI Sequence Read Archive (SRA) database (Accession Number: PRJNA1182544).
2.5 Data analyses
A one-way ANOVA (p < 0.05) test was performed to analyze differences in the rhizosphere soil environment factors, yield, and the dominant microbial abundances among different fertilizations using IBM SPSS Statistics (version 22.0, IBM Corporation, Armonk, N.Y., United States) software. The analysis of microbial subcommunity characteristics were performed using R software (version 4.3.0), and some histograms were generated using Origin 2024. Abundant and rare OTUs were defined based on cutoffs above 0.1% and below 0.01% of total sequences, as described in previous studies (Jiao et al., 2017; Jiao and Lu, 2020; Zheng et al., 2021). The “abundant subcommunity” included OTUs with a mean relative abundance >0.1% in root endophytic bacterial and fungal communities across all samples, while the “rare subcommunity” included OTUs with a mean relative abundance <0.01%. The relative abundance of bacterial and fungal subcommunities at the phylum level was constructed using the “microtable” package in R (Liu et al., 2020). The Kruskal–Wallis test for α-diversity indexes (Shannon and Richness), non-metric multidimensional scaling (NMDS), and analysis of similarities (ANOSIM) were performed using the “vegan” package in R (Liu et al., 2015; Oksanen et al., 2012; Yu et al., 2024). The beta diversity was analyzed using the “betapart” package in R (Baselga and Orme, 2012). In this study, the Bray–Curtis-based Raup-Crick (RCbray) and the β-nearest taxon index (βNTI) were used to calculate the variations in phylogenetic and taxonomic diversity using the “comdistnt” function from the “picante” R package (Jiao and Lu, 2020). Quantitatively, the proportion of inferred bacterial and fungal subcommunity assembly ecological processes in root endophytes was analyzed using a phylogenetic bin-based null model analysis (with 999 randomizations), implemented via the “icamp.big” function in R (Ma et al., 2022a; Ning et al., 2020; Wu et al., 2023). Both homogeneous and variable selections were considered as forms of species sorting (deterministic process) (Jiao et al., 2019b). Based on the co-occurrence networks to explore the interrelationships of all OTUs in bacterial and fungal subcommunities, except for the rare and abundant subcommunities, the remaining OTUs were classified as “Others.” Robust correlations were identified using Spearman’s correlation coefficients (r) > 0.6 and false discovery rate-corrected p values <0.05 (p value corrected by FDR multiple test). These metrics were applied to construct random networks of root endophytes and four endophytic subcommunities using the “igraph” package in R (Hartman et al., 2018; Jiao et al., 2019b). The topological characteristics of the networks were calculated based on the following metrics: average clustering coefficient, average connectivity, average path distance, and modularity. Furthermore, the keystone species for endophytic subcommunity networks were identified based on the threshold values of Zi (measuring how well a node was connected to other nodes within its module) and Pi (measuring how well a node was connected to nodes between different modules), set at 2.5 and 0.62, respectively (Guimerà and Nunes Amaral, 2005; Lang et al., 2021; Liang et al., 2016; Wu et al., 2023). To identify the key root endophytic subcommunity and taxa predictors for yield, the “randomForest” package in R was used, and each predictor was assessed using the “rfPermute” package. A Mantel test was conducted to determine the correlation between community composition and rhizosphere soil environment factors. A structural equation model (SEM) was analyzed using IBM SPSS Amos 26 to quantify the direct and indirect effects of the rhizosphere environment factors (such as pH, AN, and MLP) and the α and β diversities of the bacterial subcommunities on the yield. All variables were standardized using Z transformation (mean = 0, standard deviation = 1) with the ‘scale’ function in the R package. The maximum likelihood method was employed to fit the model (Jin et al., 2020). The standardized total effects (STEs) of the rhizosphere soil nutrients and root endophytic bacterial subcommunities on the yield were calculated.
3 Results
3.1 Response of the rhizosphere soil physicochemical parameters to long-term fertilization
Long-term fertilization significantly increased the total nutrient content of the rhizosphere soil excluding NH4+-N (Table 1). The contents of SOC, TN, AN, and NO3−-N significantly increased in ONP0 and ONP1, by 175%–214%, 71%–84%, 362%–506%, and 956%–1,672%, respectively, compared with CKP0. Different fractions of phosphorus increased with the application of phosphorus and organic manure, with the highest increases in the contents of TP, AP, LP, and NLP observed in ONP1, by 284%, 426%, 330%, and 155%, respectively. The MLP content in ONP1, CKP1, and CNP1 significantly increased by 539%, 587%, and 756%, respectively, compared to CKP0. Soil pH significantly decreased across all treatments compared to CKP0, except in ONP0.
Table 1
Table 1. Effects of long-term fertilization on rhizosphere soil physicochemical parameters.
3.2 Response of the endophytic microbial subcommunity characteristics to long-term fertilization
In tomato roots, the total number of 3,122 bacterial OTUs and 487 fungal OTUs were observed (Supplementary Table S3). A high percentage of OTUs were recognized as rare taxa, constituting rare subcommunities of bacteria and fungi. These rare taxa represented an average of 84.30% and 55.44% of bacterial and fungal OTUs, respectively, but contributed only 6.04% and 0.91% to the average relative abundance in all samples. Conversely, abundant bacterial and fungal taxa constituted a small proportion of OTUs with mean values of 2.85% and 8.00% but accounted for 81.01% and 94.07% of the average relative abundance in all samples, respectively. The rare bacterial subcommunity exhibited higher α-diversity, in terms of both richness and Shannon indexes, compared to the abundant subcommunity under the same treatments. The rare fungal subcommunity showed only a higher Shannon index compared to abundant subcommunity under the same treatments (Figure 1). The Shannon index of abundant bacterial and fungal subcommunities significantly decreased in CNP0 and ONP0 compared to CKP0 (p < 0.05; Figures 1A,C). The richness of abundant and rare fungi in CNP1 and ONP0 also significantly decreased compared to CKP0 (p < 0.05; Figures 1G,H). There was no significant difference in the Shannon index of rare bacterial and fungal subcommunities or the richness of abundant and rare bacterial subcommunities across other treatments, compared to CKP0 (p < 0.05; Figures 1B,D–F).
Figure 1
Figure 1. Shannon and richness indexes of root endophytic abundant bacteria (A,E), rare bacteria (B,F), abundant fungi (C,G), and rare fungi (D,H) under different treatments. CKP0, Unfertilized; CKP1, Chemical phosphorus applied; CNP0, Chemical nitrogen applied; CNP1, Chemical nitrogen plus phosphorus; ONP0, Organic manure supply nitrogen; ONP1, Organic manure supply nitrogen plus phosphorus. *p < 0.05, as tested by multiple comparisons using the Kruskal–Wallis test, and only significant differences are displayed. The treatment denotations are the same as those in Figure 1.
Root endophyte OTUs were classified into 10 bacterial and five fungal phyla, excluding unclassified sequences. The abundant bacterial subcommunity was predominantly Proteobacteria (2.81%–54.50%), whereas the rare bacterial subcommunity was also dominated by Proteobacteria (22.15%–36.75%) (Figures 2A,B). The relative abundance of Proteobacteria in the abundant bacterial subcommunity was significantly decreased by nitrogen application (CNP0 and ONP0) compared to CKP0 (p < 0.05; Supplementary Table S4). The relative abundance of Actinobacteriota, Chloroflexi, and Gemmatimonadota in rare bacterial subcommunity was significantly increased in ONP1, CNP1, and CNP1, respectively, compared to CKP0 (p < 0.05; Supplementary Table S4). The relative abundance of Acidobacteriota was significantly increased in CKP1, CNP0, and CNP1 compared to CKP0 (p < 0.05; Supplementary Table S4). The abundant fungal subcommunity was dominated by Ascomycota (0.8%–13.36%) (Figure 2C), while the rare fungal subcommunity was dominated by Ascomycota (25.23%–51.43%) and Basidiomycota (8.8%–29.63%) (Figure 2D). There was no significant difference in the relative abundance of endophytic fungal phyla among different treatments (p < 0.05; Supplementary Table S4). The NMDS analysis of Bray–Curtis distances revealed a distinct separation of all OTUs composition between treatments (beta-diversity). This revealed a clear separation between treatments of the bacterial and fungal subcommunities (Figures 2E–H). The ANOSIM analysis further indicated that the composition of bacterial subcommunity significantly varied with N levels, whereas the composition of the fungal subcommunity significantly varied with P levels rather than among different treatments (p ≤ 0.05; Supplementary Table S5).
Figure 2
Figure 2. Subcommunity composition at the phylum level and non-metric multidimensional scaling (NMDS) plots of root endophytic abundant bacteria (A,E), rare bacteria (B,F), abundant fungi (C,G), and rare fungi (D,H) under different treatments. Unclassified bacterial and fungal OTUs were labeled as “Others.” The treatment denotations are the same as those in Figure 1.
3.3 Community assembly of root endophytic abundant and rare microbial subcommunities
β-nearest taxon index values were used to quantify the relative contributions of ecological processes in different fertilization treatments (Figures 3A–D). A significant discrepancy in βNTI values of rare fungal and abundant bacterial subcommunities was caused by different N and P levels, respectively (p ≤ 0.05). There was no significant discrepancy in βNTI of different microbial subcommunities under different treatments and were all between −2 and + 2. Furthermore, the null model analysis revealed that the ecological assembly of root endophytic abundant and rare subcommunities varied with different treatments. Phosphorus levels mainly governed the abundant bacterial subcommunity through dispersal limitation (23.42%–96.79%), which, and an undominated process (0%–76.58%), which gradually increased (p < 0.05; Figures 3A,E). Dispersal limitation (12.26%–50.22%) and undominated (39.82%–77.72%) processes mainly governed the rare bacterial subcommunity (Figure 3F). The abundant fungal subcommunity was entirely governed by stochastic processes, mainly shaped by the undominated process in CKP0, CNP0, and ONP0 and increase dispersal limitation process in P level (CKP1, CNP1, and ONP1) (Figure 3G). The rare fungal subcommunity was mainly governed by nitrogen levels through decreased dispersal limitation (13.73%–72.07%) and increased undominated process (24.82%–84.56%) (Figure 3H). In summary, the assembly of abundant bacterial, abundant fungal, and rare fungal subcommunities was governed by dispersal limitation and undominated processes, while stochastic processes dominated the rare bacterial subcommunity and gradually replaced deterministic processes completely under different treatments (Figure 3).
Figure 3
Figure 3. Root endophytic bacterial and fungal subcommunities assembly processes (based on βNTI and RCBray). βNTI of root endophytic abundant bacteria (A), rare bacteria (B), abundant fungi (C), and rare fungi (D) subcommunities. Proportions of deterministic and stochastic assembly processes in governing root endophytic abundant bacteria (E), rare bacteria (F), abundant fungi (G), and rare fungi (H) subcommunities under different treatments.
3.4 Keystone species and network topological properties of endophytic microbial subcommunities
Co-occurrence networks under different treatments revealed that the endophytic rare subcommunities, with a higher number of nodes, played a more important role than abundant subcommunities (Figures 4A–F). Rare bacterial nodes accounted for 53.88% to 58.80%, whereas abundant bacterial nodes accounted for only 3.92% to 5.57%. Rare fungal nodes accounted for 4.02% to 7.50%, whereas abundant fungal nodes accounted for 2.10% to 3.03%. A co-occurrence network was also constructed for root endophytic abundant and rare subcommunities that conformed to the power-law distribution (Supplementary Figure S2). We observed that the numbers of nodes, edges, and average degree in the root endophytic microbial network decreased with nitrogen application (CNP0, CNP1, ONP0, and ONP1), while the modularity increased and the graph density decreased with phosphorus application (CKP1, CNP1, and ONP1) (Supplementary Figure S1). The observed potential interactions at root endophytes were predominantly positive, accounting for 78.12%–87.38% of the total linkages (Figure 4). The proportion of negative correlations in the root endophytic network of ONP1 was the highest (21.88%) compared to CKP0 (Figure 4F).
Figure 4
Figure 4. Co-occurrence networks of bacterial and fungal coexistence communities under different fertilizations (A–F). The networks were established by calculating correlations among abundant, rare, and other OTUs.
Nitrogen application reduced the modularity and average path length of the abundant bacterial subcommunity network, while increasing the clustering coefficient and average path length and reducing the average degree in the rare bacterial subcommunity network. Phosphorus application increased the modularity and average path length in the network of the abundant fungal subcommunity, while decreasing the clustering coefficient. In the rare fungal subcommunity, phosphorus application increased the clustering coefficient and average degree, but decreased the modularity and average path length (Supplementary Table S6). The topological roles of the individual subcommunity network nodes were shown in a Zi-Pi plot (Figure 5). Peripherals were completely occupied by total nodes in abundant bacterial and abundant fungal subcommunities (Figures 5A,C). We found a total of 63 keystone species in the root endophytic rare bacterial subcommunity network; one was identified as a module hub, while the remaining species served as connectors (Figure 5B). These keystone species play a significant role (2.45%) in the network, and they mainly belong to Proteobacteria, Actinobacteriota, Bacteroidota, Firmicutes, Myxococcota, Chloroflex, Acidobacteriota, and Gemmatimonadota (Supplementary Table S7). These phyla were highly enriched under ONP1, followed by CKP1. Two keystone species belonging to Ascomycota played a 0.85% role in the rare fungal subcommunity network (Figure 5D; Supplementary Table S7).
Figure 5
Figure 5. Distributions of network roles by analyzing module features, within-module connectivities (Zi) and among-module connectivities (Pi), in the endophytic abundant bacterial (A), rare bacterial (B), abundant fungal (C), and rare fungal (D) co-occurrence networks of root, respectively.
3.5 Specific subcommunity composition and taxa as predictors of yield under long-term fertilization
Compared with the CKP0, fertilization treatments significantly increased the yield except CKP1 (p < 0.05; Figure 6A). The random forest regression model was used to identify the structure of the root endophytic subcommunity, including abundant bacteria, rare bacteria, abundant fungi, and rare fungi that could predict the yield across different fertilization treatments (Figure 6B). The structure of abundant bacteria and rare bacteria significantly contribute to yield (p < 0.05; Figure 6B). Moreover, the structure of the abundant bacterial subcommunity explained 22.03% of yield variation, whereas the rare bacterial subcommunity accounted for 3.09% (Figure 6B). Additionally, the structures of bacterial subcommunities were positively related to AN content (p < 0.05; Figure 6B). On the contrary, abundant fungi and rare fungi were positively associated with soil phosphorus especially MLP (p < 0.01; Figure 6B), and both of them have no significant positive contribution to yield (Figure 6A). These four subcommunities, except the abundant subcommunity, were negatively correlated to soil pH (p < 0.05; Figure 6B). The random forest regression model further predicted the top 10 genera of the rare subcommunity (54.93%) had a greater impact than on yield variation compared to those in the abundant subcommunity (19.32%), indicating that they are more important drivers to yield variation (Figure 6). No key role genus of abundant subcommunity was identified from network, though top 10 abundant genus could predict the yield variation (Figure 6C; Supplementary Table S7). Three keystone species Arthrobacter, Microbacterium, and Sphingobium have been found in the rare microbial subcommunity. These species contribute to increased yield in response to soil nutrient changes and also serve as connectors in the network (Figure 6D; Supplementary Table S7).
Figure 6
Figure 6. Changes in tomato yield (A). The random forest regression model shows the root endophytic subcommunity structures as drivers of yield. A correlation between the structure of microbial subcommunities and soil nutrients was identified using Spearman’s test (B). The random Forest regression model identifies the top 10 most important taxa of abundant (C) and rare (D) subcommunities at the genus level, along with their respective phylum, as key drivers of tomato yield. The bars illustrate the increase in mean squared error (%IncMSE) scores for the contribution of root endophytes at the genus level among six fertilization treatments, affecting the tomato yield. Significance levels are denoted with *p < 0.05, **p < 0.01, and ***p < 0.001.
Structural equation model was employed to quantify the correlations among key soil characters (AN, pH, and MLP) in response to amounts of N and P fertilizers, the most important root endophytic bacterial subcommunities, and yield (Figure 7). The soil AN, pH, and MLP explained much more variation of the β-diversity (r2 = 0.47) than α-diversity (r2 = 0.30) in abundant bacterial subcommunity, but much more variation of α-diversity (r2 = 0.54) than β-diversity (r2 = 0.38) in rare bacterial subcommunity. Soil AN had a positive direct and indirect effect by influencing the β-diversity of abundant and rare bacterial subcommunity on yield. It should be noted that the α-diversity of abundant and rare bacterial subcommunities showed no direct effect on yield. The contrasting results were observed for yield, being directly negatively affected by pH and MLP, while through the α-diversity of rare bacteria directly related to MLP (Figure 7A). The soil pH and MLP had a positive effect and AN had a negative effect on the α-diversity of abundant and rare subcommunities, but its effects were contrary on the β-diversity of abundant and rare bacterial subcommunities. The soil AN was found to have the largest standardized positive effect on yield, followed by the β-diversity of abundant and rare bacterial subcommunities, while the contribution of the abundant bacterial β-diversity was greater compared to that of the rare bacterial β-diversity (Figure 7B).
Figure 7
Figure 7. Structural equation model (SEM) illustrating how the root endophytic bacterial subcommunities influenced yield by responding to soil AN, pH, and MLP (A). The total standardized effects of these factors on yield, α-diversity (Shannon diversity), and β-diversity are shown (B). Red arrows, black arrows, and dashed arrows indicate the positive significant, negative significant, and nonsignificant relationships between different variables. The width of the arrows indicates the strength of the standardized path coefficient. Adjacent values near the arrows indicate path coefficients. r2 values indicate the proportion of variance explained by each variable. Significance levels are denoted with *p < 0.05, **p < 0.01, and ***p < 0.001.
4 Discussion
4.1 Long-term fertilization affected the construction processes in distribution characters of root endophytes subcommunities
The abundant endophytic bacterial and fungal subcommunities were prevalent in most fertilization treatments with an average relative abundance of 81.01% and 94.07%, respectively. Abundant species may occupy wider niches, competitively utilize more enriched resources, have high growth rates, and effectively adapt to the environment. All of these supported their persistence and their widespread distribution (Jiao and Lu, 2020; Liu et al., 2023a). In addition, the relative abundance of root endophytic rare bacterial and fungal subcommunities that contributed less to the total community abundance was unevenly distributed among different fertilization treatments. The different distribution patterns of root endophytic abundant and rare taxa may be attributed to their distinct life strategies. The rare subcommunities usually exhibited dormancy or low growth rates and competitiveness. Hence, root endophytic microbial subcommunities showed limitations but not extinction in nutrient-enriched environments caused by long-term fertilization (Jia et al., 2018; Jousset et al., 2017). Notably, the rare subcommunities had the highest abundance and diversity compared with abundant subcommunities, which supported the fact that the rare subcommunity was an important contributor to the microbial diversity (du et al., 2020; Li et al., 2021; Lynch and Neufeld, 2015). This trend may be attributed to the fact that a community has a limited capacity and excessive diversity was formed by species competition. Rare taxa can become abundant members of a community in response to favorable environmental conditions (Jia et al., 2022). Few rare species dominated the root samples in our study, but rare bacterial taxa were more susceptible to environmental heterogeneity and with ability to grow fast were affected by nutrient concentration increasing as conditionally rare taxa (Jia et al., 2022; Jiao et al., 2019a). This is possibly due to the selective pressure of resources in the rhizosphere soil environment driving the aggregation of root endophytic rare microbial communities (Liu et al., 2024). It is known that soil environment factors such as carbon, nitrogen, and phosphorus concentrations could across compartments affect root species abundance and interactions (Zhang et al., 2019). These findings suggest that the aggregation process of microbiome with different nutritional strategies can affect the diversity and abundance of endophytic microbial subcommunities (Li et al., 2021; Liu et al., 2023a). Similar evidence has been found in previous studies to explain this change (Dos Santos et al., 2022), where the alpha diversity of endophytic microorganisms increases with the application of fertilizers. However, the α-diversity of fungal subcommunities, especially the α-diversity of the rare fungal subcommunity, was mainly increased by the application of chemical nitrogen and phosphorus fertilizers in our study (p < 0.05). Similarly, we found that different N levels and P levels shaped complex rhizosphere soil nutrient heterogeneity that significantly affected the structure of root endophytic bacterial and fungal subcommunities (Ma et al., 2021; Ma et al., 2022b; Zhang et al., 2019). This may have mediated a richer diversity of the rare subcommunity and enhanced their involvement in nutrient cycling processes between environments (Liu et al., 2023a; Shu et al., 2021).
The microbial community assembly processes, which determine the presence and abundance of species, are generally divided into deterministic and stochastic processes (Feng et al., 2018; Jiao et al., 2019b). The null model suggested that the construction of root endophytic bacterial and fungal subcommunities are mainly dominated by stochastic processes especially dispersal limitation (Zhong et al., 2022). More importantly, the assembly process of abundant and rare subcommunities of endophytic bacteria and fungi exhibited different trade-off patterns, reflecting varying relative contributions in response to changes in nitrogen and phosphorus levels (Xing et al., 2023). The proportion of stochastic processes in microbial community assembly can also increase through the functional redundancy generated by biological interactions and spatial and environmental processes (Li and Gao, 2023). Under different fertilization treatments of nitrogen and phosphorus levels, rhizosphere environmental filtration emphasizes the stochastic process of dispersal limitation, which may be involved in determining the development of different sub community systems, which was a very novel finding. Previous studies have revealed that determinism seems to be related to low nutrient levels, and stochasticity increases with nutrient levels (Chase, 2010; Zhou et al., 2014). Therefore, the current nutritional strategy for microbial community assembly processes may be more important than the priority effects shaped by host selection (Li and Gao, 2023). In particular, as the subcommunity rarity increases, dispersal limitation was more conducive to driving the assembly and construction of root endophytic bacterial and fungal subcommunities in enough nitrogen and phosphorus environments, respectively. This validation showed that rare microbial subcommunities were not a random taxonomic cluster; instead, most of these rare subcommunities show restricted distribution or migratory diffusion (Cui et al., 2023; Liu et al., 2015; Long and Yao, 2020). Compared to rare taxa, more individuals in abundant taxa were likely to be involved in dispersal limitation (Liu et al., 2015), such as abundant bacterial subcommunities in our study. From this perspective, abundant taxa are likely to adopt growth yield or resource acquisition strategies that enable them to outcompete other microorganisms for broader niche space and nutrients (Cui et al., 2023; Delgado-Baquerizo et al., 2018; Dong et al., 2021; Jiao and Lu, 2020; Ma et al., 2022a). Consequently, we concluded that long-term fertilization indeed affected the specific niche differentiation of phosphorus-sensitive fungal and nitrogen-sensitive root endophytic bacterial, both abundant and rare subcommunities (Hamonts et al., 2018; Leach et al., 2017; Niu et al., 2017).
The rare fungal subcommunity exhibited a more intense and significant response to phosphorus than the abundant fungal subcommunity, possibly due to the long-life cycle, low community turnover rate, and differences in metabolic activity and diffusion ability of fungi (Jiao and Lu, 2020; Omomowo and Babalola, 2019). Higher nutrient levels increased the importance of stochastic processes (Chase, 2010; Chen et al., 2020), providing insights into the assembly patterns of root endophytic fungal subcommunities promoted by phosphorus. Nitrogen application maintained the control of the dispersal limitation on the enrichment of both abundant and rare bacterial subcommunities, significantly impacting their structure. Previous studies found that short-term high nitrogen levels have a greater effect on the interactions among endophytic fungal communities in maize than on endophytic bacteria (Bai et al., 2022). We found community assembly processes and ANOSIM analysis suggested that endophytic bacteria respond more strongly to nitrogen fertilizer than endophytic fungi. This difference may be partly because nutrient supplementation reduces the competitive advantage of fungi over bacteria (Michalska-Smith et al., 2022) and increases bacterial abundance (Wu et al., 2019), especially in rare bacterial subcommunity. On the other hand, bacteria rely more on fixed sources of nitrogen and carbon than fungi (Schmidt et al., 2014), and organic fertilizers also compensate for the relative lack of carbon when achieving nitrogen supply in this process. These indicated that these subcommunities were formed by more complex assembly processes, especially rare subcommunities, which may also be important metabolic strategies for microbial cells to cope with environmental pressures and be less susceptible to deterministic processes (Fujii et al., 2012; Lennon and Jones, 2011; Nemergut Diana et al., 2013). It is worth noting that stochastic process plays more importance in community assembly, leading to higher β-diversity of microorganisms and higher yield (Chase, 2010). The altered community aggregation characteristics of different subcommunities of root endophytic bacteria and fungi have improved the overall adaptability of microbial communities to their environment, potentially enhancing our understanding of ecosystem sustainability through various means. Rare subcommunities may serve as sensitive biological indicators for regulating their assembly mechanisms and reflecting environmental change risks, closely related to the composition of the entire endophytic microbial community. They may also be a key factor in enhancing the resistance of microbial networks through functional redundancy (Li et al., 2021; Lynch and Neufeld, 2015).
4.2 Co-occurrence network of root endophytes subcommunities and keystone species response to long-term fertilization
Network analysis has been identified as a powerful tool for exploring the role of microbial symbiotic interactions in coping with environmental stress (Michalska-Smith et al., 2022; Shu et al., 2021). Fertilizations have affected the network structural characteristics, providing insights into how symbiotic patterns of microbe form different endophytic microbial communities (Liu et al., 2023a) and helping to identify potential keystone taxa in network (Wu et al., 2022; Zhong et al., 2022). Our finding indicated that the frequency, intensity, and distribution patterns of OTU interactions within subcommunities of root endophytic bacteria and fungi vary in response to different fertilization treatments, particularly changed with nitrogen and phosphorus levels. These findings are consistent with previous studies, showing that the cross soil and plant compartments of nutrient concentrations significantly altered overall community dynamics by affecting niche differentiation and colonization of root endophytic bacteria and fungi (Michalska-Smith et al., 2022; Sun et al., 2022; Xu et al., 2022). This niche dispersal selectively recruits microbial subcommunities and formed stable coexistence through root interactions (Ernst et al., 2011). After phosphorus application (CKP1, CNP1, and ONP1), root endophytic microbial communities exhibited higher modularity and lower graph density topological parameters, with an increased proportion of negative cohesion, indicative of a more stable and complex symbiotic pattern. Positive correlations among microbial species have indicated cooperative relationships or niche overlap, while negative correlations have indicated competitive relationships or niche separation (Ernst et al., 2011; Ghoul and Mitri, 2016). Due to higher soil nutrient levels, the highest network complexity and microbial growth diversity have been observed in ONP1, suggesting high nutrient competition among taxa under ONP1 (Deng et al., 2016). This has been verified by the enhanced negative cohesion in the root endophytic symbiotic network under ONP1 (increased by 9.26% and 8.10% compared to CKP0 and ONP0, respectively). Such a complex symbiotic network increased community structure stability, enabling resistance to long-term chemical fertilizer-induced disturbances.
The topological characteristics of the subnetwork indicated that rare subcommunities exhibited more complex and stable symbiotic associations compared to abundant subcommunities (Li et al., 2021). We found that rare bacterial and fungal subcommunities have shown stronger clustering with nitrogen and phosphorus applications, respectively, reflecting the importance of ecological assembly processes and cooperative nutrient metabolism pathways among taxa in the endophytic microbiome (Cui et al., 2023; Fabiańska et al., 2019; Ma et al., 2022b; Sun et al., 2022). The tight associations in the symbiotic network may be the contribution of rare bacterial and rare fungal subcommunities, which were driven by dispersal limitation processes and achieved functional complementarity from different taxonomic compositions. Moreover, under nitrogen level treatments (CN, ON), although the number of nodes and interactions in the subnetwork have decreased, a higher proportion of negative cohesion helped to form more robust competitive relationships and more stable subnetwork, which can be explained by enrichment-dominated and depletion-supplemented patterns (Liu et al., 2024). The complex network structures indicated a more efficient potential for resource and information transmission within the community, allowing it to withstand environmental stress.
Studies have shown that chemical nitrogen fertilizer has significantly affected the structure of root endophytic bacterial communities (Miranda-Carrazco et al., 2022) and reduced the diversity and stability of soil microbial networks (Zhang et al., 2021). Our findings are similar to previous studies, but happened in rare bacterial subnetwork rather than abundant bacterial subnetwork, possibly because the subnetwork structure changes of abundant bacteria have not weakened the core subcommunities or species dominance. Compared to chemical nitrogen fertilizer, long-term organic manure application has likely supported taxonomic groups with similar functions, thereby increasing functional redundancy (Fan et al., 2019). This has a stronger positive effect on maintaining the complexity and stability of the root endophytic bacterial symbiotic network (Zhang et al., 2021), which further enhances the resistance toward the negative effect caused by chemical fertilizers. In general, the root endophytic abundant microbial subcommunities had high stability under different fertilization treatments, which indicated that it could play a dominant role in the symbiotic network. At the same time, the rare subcommunities and their key species may interact with abundant subcommunities under environmental disturbances, forming more stable symbiotic network and performing specific ecological functions (Jiao et al., 2017; Jiao et al., 2019a; Li et al., 2021).
Previous research has also provided evidence that the presence or loss of rare subcommunity and their key species under environmental disturbances can affect the stability of interactions within the entire community, leading to significant structural and functional changes (Cui et al., 2023; Ma et al., 2024). We found that the topological characters of nitrogen-adaption rare bacterial subcommunity tend to align with the structural changes of the entire endophytic microbiome. This further verify that rare endophytic bacterial subcommunity play a significantly more important role than endophytic fungal subcommunities (Liu et al., 2024). The rare bacterial subcommunity contained the most of key species occupying the nodes of module hubs and connectors in the network, and their abundance increases with the application of nitrogen and phosphorus fertilizers. These keystone species can disproportionately and significantly impact ecosystems, with their activity and abundance determining community integrity (Wu et al., 2023). At the phylum level, keystone species in rare bacterial subcommunities belong primarily to Proteobacteria, Actinobacteriota, Bacteroidota, Firmicutes, Myxococcota, Chloroflex, Acidobacteriota, Gemmatimonadota, and unclassified bacteria.
High nitrogen and phosphorus nutrient concentrations have been unfavorable for the growth of oligotrophic bacteria such as Acidobacteriota and Chloroflexi (Fierer et al., 2012). Long-term nitrogen and phosphorus fertilizer application has increased the relative abundance of copiotrophic bacteria, including Alphaproteobacteria, Gammaproteobacteria, Firmicutes (Ho et al., 2017), Bacteroidota, and Actinobacteriota (Dai et al., 2018; Fierer et al., 2012; Sun et al., 2023). Similarly, nitrogen and phosphorus fertilizers have significantly promoted the growth of Ascomycota, a hub in the root endophytic fungal subcommunity network, which is considered an r-strategy copiotrophic fungus (Guo et al., 2020). Ascomycota is often reported to harbor a high number of genes relevant to stress tolerance and nutrition uptake (e.g., nitrogen immobilization and phosphate transporter) (Egidi et al., 2019). Our results indicated that long-term nitrogen and phosphorus fertilizer applications enriched copiotrophic root endophytic nitrogen-adaption bacteria and phosphorus-adaption fungi, which may possess symbiotic ecological functions (Wu et al., 2023). This enrichment may benefit crop health and high yield in agricultural ecosystems. The most direct reason is that nitrogen and phosphorus fertilizers increased nutrient availability, which may differentially impact communication within bacterial subcommunities in the network, subsequently inducing the transformation of various key groups within the entire microbial community. Therefore, we concluded that long-term balanced fertilization with nitrogen and phosphorus, especially with organic manure and chemical phosphorus fertilizer (ONP1), may enrich more beneficial and important tomato root endophytes, thereby improving microbial network stability and crop yield.
4.3 The importance of endophytic bacterial subcommunity structure driven by AN and MLP to yield
Rhizosphere soil nutrient characteristics, particularly nitrogen, phosphorus, and their availability, are potential key factors influencing composition and aggregation of root endophytic microbial community (Liu et al., 2023b; Long and Yao, 2020; Meng et al., 2022; Sindhu et al., 2024). Studies indicated that soil changes are the initial step affecting root endophyte colonization. Although the impact of soil is relatively smaller than that of the host, nitrogen and phosphorus availability still influence endophytic microbial community aggregation and structure (Almario et al., 2017; Edwards et al., 2015). Previous research has shown that different nitrogen levels lead to differences in root endophytic bacteria (Beltran-Garcia et al., 2021). We discovered that different nitrogen supplements caused inconsistent changes in endophytic bacterial community diversity, but structures of endophytic bacteria had similar responses to same rhizosphere soil physicochemical parameters. Long-term different nitrogen fertilizer significantly increased AN content, which showed a significant positive correlation with the structure of both abundant and rare endophytic bacterial subcommunities. Additionally, the reduction in pH was strongly negatively correlated with the structure of rare endophytic bacterial subcommunities, suggesting that pH involved in determining specific endophytic bacterial communities (Lin et al., 2020; Zhang et al., 2024). Previous studies have indicated that endophytic fungal communities are closely related to phosphorus content (Wu et al., 2022). Chemical phosphorus fertilizer provided substantial amounts of labile and moderately labile phosphorus (Maranguit et al., 2017; Shi et al., 2023), and continuous changes in phosphorus availability promote the turnover of endophytic fungal subcommunity structures (Long and Yao, 2020). In non-mycorrhizal crops like tomatoes, it has been found that root endophytic fungi may serve as an alternative strategy for phosphorus utilization, supporting a significant response of endophytic fungal subcommunity to changes in rhizosphere soil MLP content.
Random forest model regression analysis revealed that the most important endophytic subcommunities for predicting crop yield were abundant endophytic bacterial subcommunities, followed by rare bacterial subcommunity, accounting for 22.03% and 3.09% of the variance in predictive factors, respectively. This finding differs from studies on maize-wheat/barley rotation systems (Xiong et al., 2020). In tomato roots affected by long-term fertilization, rare subcommunities structures and species that promote higher network stability may synergistically assist the positive effects of abundant subcommunities with environmental resistance on crop yield. Consequently, the associations between rhizosphere soil nitrogen and phosphorus availability and root endophytic bacterial and fungal subcommunities suggest a pattern of endophytic microbial structural succession, which optimizes microbial communities to influence the relative contributions of different subcommunities to crop yield. Notably, among the top 10 rare subcommunity predictors of yield, three key species were identified that Arthrobacter, Microbacterium, and Sphingobium. Interestingly, both Arthrobacter and Microbacterium belong to the Actinobacteriota, while Sphingobium belongs to Proteobacteria. These key rare taxa, as plant growth-promoting endophytes (PGPE), have potential to enhance plant resilience to environmental stress and improve crop health and yield, possessing specific functions such as nitrogen fixation and phosphorus solubilization (Machado et al., 2020; Walitang et al., 2017).
5 Conclusion
This study provided a systematic understanding of how root endophytic microbial subcommunities, adapted to different ecological niches under long-term fertilization, contribute to maintaining tomato yield. Long-term nitrogen and phosphorus fertilizers increased the proportion of negative cohesion in the root endophytic microbial co-occurrence network, forming ecological niche complementarity between rare bacterial subcommunities and fungal subcommunities. This effect was particularly pronounced in the ONP1 treatment, which led to a more complex, tighter, and stable network structure. The assembly of rare and abundant subcommunities of endophytic bacteria and fungi were dominated by stochastic processes. Nitrogen applications maintained the dispersal limitation ecological process, influencing the assembly process of endophytic rare and abundant bacterial subcommunities, while phosphorus applications controlled the assembly of rare and abundant subcommunities of endophytic fungi by increasing the proportion of the dispersal limitation process. The AN and MLP content in rhizosphere soil significantly affected the community structure changes of nitrogen-adaption root endophytic bacteria and phosphorus-adaption endophytic fungi. The nitrogen-adaption endophytic bacterial subcommunity, particularly the abundant bacterial subcommunity, had a higher explanatory power for tomato yield than the phosphorus-adaption fungal subcommunity. Arthrobacter, Microbacterium, and Sphingobium from the rare subcommunity can potentially serve as beneficial symbiotic endophytes to enhance crop yield. Overall, the application of organic manure combined with chemical phosphorus fertilizer favored the survival strategy of abundant subcommunities and the diversity of rare subcommunities. Ecological complementarity among different subcommunities contributes to forming a more complex and tighter symbiotic network, thereby driving functional redundancy in rare subcommunities potential for high crop yield. The findings have significant implications for sustainable fertilization practices to control root endophytic microbial communities for achieving high crop yield.
Data availability statement
The datasets from this study are available in online repositories, with repository names and accession numbers provided in the article.
Author contributions
XL: Conceptualization, Data curation, Formal analysis, Writing – original draft. MA: Writing – review & editing. SW: Formal analysis, Writing – review & editing. ZZ: Investigation, Formal analysis, Writing – original draft. SZ: Investigation, Visualization, Writing – review & editing. YL: Writing – review & editing. ZS: Funding acquisition, Supervision, Writing – review & editing. HF: Supervision, Writing – review & editing. TL: Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was financially supported by the National Key R&D Program of China (2022YFD1602110) and the China Agriculture Research System (CARS-23). The funders had no role in the execution of the study or the decision to publish the findings.
Acknowledgments
The authors thank Lijun Chen for his contributions in revising the manuscript.
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1487323/full#supplementary-material
References
Adeleke, B. S., and Babalola, O. O. (2021). The endosphere microbial communities, a great promise in agriculture. Int. Microbiol. 24, 1–17. doi: 10.1007/s10123-020-00140-2
Adeleke, B. S., Babalola, O. O., and Glick, B. R. (2021). Plant growth-promoting root-colonizing bacterial endophytes. Rhizosphere 20:100433. doi: 10.1016/j.rhisph.2021.100433
Ali, S., Charles, T. C., and Glick, B. R. (2014). Amelioration of high salinity stress damage by plant growth-promoting bacterial endophytes that contain ACC deaminase. Plant Physiol. Biochem. 80, 160–167. doi: 10.1016/j.plaphy.2014.04.003
Almario, J., Jeena, G., Wunder, J., Langen, G., Zuccaro, A., Coupland, G., et al. (2017). Root-associated fungal microbiota of nonmycorrhizal Arabis alpina and its contribution to plant phosphorus nutrition. Proc. Natl. Acad. Sci. USA 114, E9403–e9412. doi: 10.1073/pnas.1710455114
Bai, L., Zhang, X., Li, B., Sun, F., Zhao, X., Wang, Y., et al. (2022). Fungal communities are more sensitive to nitrogen fertilization than bacteria in different spatial structures of silage maize under short-term nitrogen fertilization. Appl. Soil Ecol. 170:104275. doi: 10.1016/j.apsoil.2021.104275
Banerjee, S., Walder, F., Büchi, L., Meyer, M., Held, A. Y., Gattinger, A., et al. (2019). Agricultural intensification reduces microbial network complexity and the abundance of keystone taxa in roots. ISME J. 13, 1722–1736. doi: 10.1038/s41396-019-0383-2
Baselga, A., and Orme, C. D. L. (2012). Betapart: an R package for the study of beta diversity. Methods Ecol. Evol. 3, 808–812. doi: 10.1111/j.2041-210X.2012.00224.x
Battu, L., and Ulaganathan, K. (2020). Whole genome sequencing and identification of host-interactive genes in the rice endophytic Leifsonia sp. ku-ls. Funct. Integr. Genomics 20, 237–243. doi: 10.1007/s10142-019-00713-z
Beltran-Garcia, M. J., Martínez-Rodríguez, A., Olmos-Arriaga, I., Valdes-Salas, B., Di Mascio, P., and White, J. F. (2021). Nitrogen fertilization and stress factors drive shifts in microbial diversity in soils and plants. Symbiosis 84, 379–390. doi: 10.1007/s13199-021-00787-z
Bolan, N. S. (1991). A critical review on the role of mycorrhizal fungi in the uptake of phosphorus by plants. Plant Soil 134, 189–207. doi: 10.1007/BF00012037
Bornemisza, E., Blakemore, LC, Searle, PL, and Daly, BK (1988). Methods for chemical analysis of soils. New Zealand soil bureau scientific report, 1987, 103 p.
Bulgarelli, D., Rott, M., Schlaeppi, K., Loren, V., van Themaat, E., Ahmadinejad, N., et al. (2012). Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. Nature 488, 91–95. doi: 10.1038/nature11336
Burragoni, S. G., and Jeon, J. (2021). Applications of endophytic microbes in agriculture, biotechnology, medicine, and beyond. Microbiol. Res. 245:126691. doi: 10.1016/j.micres.2020.126691
Byregowda, R., Prasad, S. R., Oelmüller, R., Nataraja, K. N., and Prasanna Kumar, M. K. (2022). Is endophytic colonization of host plants a method of alleviating drought stress? Conceptualizing the hidden world of endophytes. Int. J. Mol. Sci. 23:9194. doi: 10.3390/ijms23169194
Cabrera, M. L., and Beare, M. H. (1993). Alkaline persulfate oxidation for determining total nitrogen in microbial biomass extracts. Soil Sci. Soc. Am. J. 57, 1007–1012. doi: 10.2136/sssaj1993.03615995005700040021x
Card, S., Johnson, L., Teasdale, S., and Caradus, J. (2016). Deciphering endophyte behaviour: the link between endophyte biology and efficacious biological control agents. FEMS Microbiol. Ecol. 92, 1–19. doi: 10.1093/femsec/fiw114
Chase, J. M. (2010). Stochastic community assembly causes higher biodiversity in more productive environments. Science 328, 1388–1391. doi: 10.1126/science.1187820
Chase, J. M., and Myers, J. A. (2011). Disentangling the importance of ecological niches from stochastic processes across scales. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 366, 2351–2363. doi: 10.1098/rstb.2011.0063
Chen, G., Wang, W., Zhang, Y., Liu, Y., Gu, X., Shi, X., et al. (2020). Abundant and rare species may invoke different assembly processes in response to climate extremes: implications for biodiversity conservation. Ecol. Indic. 117:106716. doi: 10.1016/j.ecolind.2020.106716
Collavino, M. M., Cabrera, E. V. R., Bruno, C., and Aguilar, O. M. (2020). Effect of soil chemical fertilization on the diversity and composition of the tomato endophytic diazotrophic community at different stages of growth. Braz. J. Microbiol. 51, 1965–1975. doi: 10.1007/s42770-020-00373-3
Cui, W., Li, R., Fan, Z., Wu, L., Zhao, X., Wei, G., et al. (2023). Weak environmental adaptation of rare phylotypes sustaining soil multi-element cycles in response to decades-long fertilization. Sci. Total Environ. 871:162063. doi: 10.1016/j.scitotenv.2023.162063
Dai, Z., Su, W., Chen, H., Barberán, A., Zhao, H., Yu, M., et al. (2018). Long-term nitrogen fertilization decreases bacterial diversity and favors the growth of Actinobacteria and Proteobacteria in agro-ecosystems across the globe. Glob. Chang. Biol. 24, 3452–3461. doi: 10.1111/gcb.14163
Delgado-Baquerizo, M., Oliverio, A. M., Brewer, T. E., Benavent-González, A., Eldridge, D. J., Bardgett, R. D., et al. (2018). A global atlas of the dominant bacteria found in soil. Science 359, 320–325. doi: 10.1126/science.aap9516
Deng, Y., Zhang, P., Qin, Y., Tu, Q., Yang, Y., He, Z., et al. (2016). Network succession reveals the importance of competition in response to emulsified vegetable oil amendment for uranium bioremediation. Environ. Microbiol. 18, 205–218. doi: 10.1111/1462-2920.12981
Dong, Y., Wu, S., Deng, Y., Wang, S., Fan, H., Li, X., et al. (2021). Distinct functions and assembly mechanisms of soil abundant and rare bacterial taxa under increasing pyrene stresses. Front. Microbiol. 12:689762. doi: 10.3389/fmicb.2021.689762
Dos Santos, R. M., Cueva-Yesquén, L. G., Garboggini, F. F., Desoignies, N., and Rigobelo, E. C. (2022). Inoculum concentration and mineral fertilization: effects on the endophytic microbiome of soybean. Front. Microbiol. 13:900980. doi: 10.3389/fmicb.2022.900980
du, S., Dini-Andreote, F., Zhang, N., Liang, C., Yao, Z., Zhang, H., et al. (2020). Divergent co-occurrence patterns and assembly processes structure the abundant and rare bacterial communities in a salt marsh ecosystem. Appl. Environ. Microbiol. 86:e00322. doi: 10.1128/AEM.00322-20
Edgar, R C. (2016a). UNOISE2: improved error-correction for Illumina 16S and ITS amplicon sequencing. bioRxiv [Preprint]. doi: 10.1101/081257
Edgar, R C. (2016b). SINTAX: a simple non-Bayesian taxonomy classifier for 16S and ITS sequences. bioRxiv [Preprint]. doi: 10.1101/074161
Edwards, J., Johnson, C., Santos-Medellín, C., Lurie, E., Podishetty, N. K., Bhatnagar, S., et al. (2015). Structure, variation, and assembly of the root-associated microbiomes of rice. Proc. Natl. Acad. Sci. 112, E911–E920. doi: 10.1073/pnas.1414592112
Egidi, E., Delgado-Baquerizo, M., Plett, J. M., Wang, J., Eldridge, D. J., Bardgett, R. D., et al. (2019). A few Ascomycota taxa dominate soil fungal communities worldwide. Nat. Commun. 10:2369. doi: 10.1038/s41467-019-10373-z
Ernst, M., Neubert, K., Mendgen, K. W., and Wirsel, S. G. R. (2011). Niche differentiation of two sympatric species of Microdochium colonizing the roots of common reed. BMC Microbiol. 11:242. doi: 10.1186/1471-2180-11-242
Fabiańska, I., Sosa-Lopez, E., and Bucher, M. (2019). The role of nutrient balance in shaping plant root-fungal interactions: facts and speculation. Curr. Opin. Microbiol. 49, 90–96. doi: 10.1016/j.mib.2019.10.004
Fan, K., Delgado-Baquerizo, M., Guo, X., Wang, D., Wu, Y., Zhu, M., et al. (2019). Suppressed N fixation and diazotrophs after four decades of fertilization. Microbiome 7:143. doi: 10.1186/s40168-019-0757-8
FAO and ISRIC. (1988). Soil map of the world, revised legend. Food and agriculture organization of the united nations, Rome, Italy.
Feng, M., Adams, J. M., Fan, K., Shi, Y., Sun, R., Wang, D., et al. (2018). Long-term fertilization influences community assembly processes of soil diazotrophs. Soil Biol. Biochem. 126, 151–158. doi: 10.1016/j.soilbio.2018.08.021
Fierer, N., Lauber, C. L., Ramirez, K. S., Zaneveld, J., Bradford, M. A., and Knight, R. (2012). Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. ISME J. 6, 1007–1017. doi: 10.1038/ismej.2011.159
Fujii, M., Kojima, H., Iwata, T., Urabe, J., and Fukui, M. (2012). Dissolved organic carbon as major environmental factor affecting Bacterioplankton communities in Mountain Lakes of eastern Japan. Microb. Ecol. 63, 496–508. doi: 10.1007/s00248-011-9983-8
Ghoul, M., and Mitri, S. (2016). The ecology and evolution of microbial competition. Trends Microbiol. 24, 833–845. doi: 10.1016/j.tim.2016.06.011
Gschwend, F., Hartmann, M., Hug, A.-S., Enkerli, J., Gubler, A., Frey, B., et al. (2021). Long-term stability of soil bacterial and fungal community structures revealed in their abundant and rare fractions. Mol. Ecol. 30, 4305–4320. doi: 10.1111/mec.16036
Guimerà, R., and Nunes Amaral, L. A. (2005). Functional cartography of complex metabolic networks. Nature 433, 895–900. doi: 10.1038/nature03288
Guo, L. D., Hyde, K. D., and Liew, E. C. Y. (2000). Identification of endophytic fungi from Livistona chinensis based on morphology and rDNA sequences. New Phytol. 147, 617–630. doi: 10.1046/j.1469-8137.2000.00716.x
Guo, Z., Wan, S., Hua, K., Yin, Y., Chu, H., Wang, D., et al. (2020). Fertilization regime has a greater effect on soil microbial community structure than crop rotation and growth stage in an agroecosystem. Appl. Soil Ecol. 149:103510. doi: 10.1016/j.apsoil.2020.103510
Guo, B., Zhang, H., Liu, Y., Chen, J., and Li, J. (2024). Assembly of cereal crop fungal communities under water stress determined by host niche. Environ. Exp. Bot. 219:105650. doi: 10.1016/j.envexpbot.2024.105650
Hamonts, K., Trivedi, P., Garg, A., Janitz, C., Grinyer, J., Holford, P., et al. (2018). Field study reveals core plant microbiota and relative importance of their drivers. Environ. Microbiol. 20, 124–140. doi: 10.1111/1462-2920.14031
Han, J., Shi, J., Zeng, L., Xu, J., and Wu, L. (2015). Effects of nitrogen fertilization on the acidity and salinity of greenhouse soils. Environ. Sci. Pollut. Res. 22, 2976–2986. doi: 10.1007/s11356-014-3542-z
Harman, G. E., and Uphoff, N. (2019). Symbiotic root-endophytic soil microbes improve crop productivity and provide environmental benefits. Australas. Sci. 2019, 9106395–9106325. doi: 10.1155/2019/9106395
Hartman, K., van der Heijden, M. G. A., Wittwer, R. A., Banerjee, S., Walser, J.-C., and Schlaeppi, K. (2018). Cropping practices manipulate abundance patterns of root and soil microbiome members paving the way to smart farming. Microbiome 6:14. doi: 10.1186/s40168-017-0389-9
Hedley, M. J., Stewart, J. W. B., and Chauhan, B. S. (1982). Changes in inorganic and organic soil phosphorus fractions induced by cultivation practices and by laboratory incubations. Soil Sci. Soc. Am. J. 46, 970–976. doi: 10.2136/sssaj1982.03615995004600050017x
Ho, A., Di Lonardo, D. P., and Bodelier, P. L. E. (2017). Revisiting life strategy concepts in environmental microbial ecology. FEMS Microbiol. Ecol. 93:fix0006. doi: 10.1093/femsec/fix0006
Jia, X., Dini-Andreote, F., and Falcão Salles, J. (2018). Community assembly processes of the microbial rare biosphere. Trends Microbiol. 26, 738–747. doi: 10.1016/j.tim.2018.02.011
Jia, X., Dini-Andreote, F., and Salles, J. F. (2022). Unravelling the interplay of ecological processes structuring the bacterial rare biosphere. ISME Commun. 2:96. doi: 10.1038/s43705-022-00177-6
Jiao, S., Chen, W., and Wei, G. (2017). Biogeography and ecological diversity patterns of rare and abundant bacteria in oil-contaminated soils. Mol. Ecol. 26, 5305–5317. doi: 10.1111/mec.14218
Jiao, S., and Lu, Y. (2020). Abundant fungi adapt to broader environmental gradients than rare fungi in agricultural fields. Glob. Chang. Biol. 26, 4506–4520. doi: 10.1111/gcb.15130
Jiao, S., Wang, J., Wei, G., Chen, W., and Lu, Y. (2019a). Dominant role of abundant rather than rare bacterial taxa in maintaining agro-soil microbiomes under environmental disturbances. Chemosphere 235, 248–259. doi: 10.1016/j.chemosphere.2019.06.174
Jiao, S., Yang, Y., Xu, Y., Zhang, J., and Lu, Y. (2019b). Balance between community assembly processes mediates species coexistence in agricultural soil microbiomes across eastern China. ISME J. 14, 202–216. doi: 10.1038/s41396-019-0522-9
Jin, S., Vegelius, J., and Yang-Wallentin, F. (2020). A marginal maximum likelihood approach for extended quadratic structural equation modeling with ordinal data. Struct. Equ. Model. Multidiscip. J. 27, 864–873. doi: 10.1080/10705511.2020.1712552
Jousset, A., Bienhold, C., Chatzinotas, A., Gallien, L., Gobet, A., Kurm, V., et al. (2017). Where less may be more: how the rare biosphere pulls ecosystems strings. ISME J. 11, 853–862. doi: 10.1038/ismej.2016.174
Lang, M., Zou, W., Chen, X., Zou, C., Zhang, W., Deng, Y., et al. (2021). Soil microbial composition and phoD gene abundance are sensitive to phosphorus level in a Long-term wheat-maize crop system. Front. Microbiol. 11:605955. doi: 10.3389/fmicb.2020.605955
Leach, J. E., Triplett, L. R., Argueso, C. T., and Trivedi, P. (2017). Communication in the Phytobiome. Cell 169, 587–596. doi: 10.1016/j.cell.2017.04.025
Lennon, J. T., and Jones, S. E. (2011). Microbial seed banks: the ecological and evolutionary implications of dormancy. Nat. Rev. Microbiol. 9, 119–130. doi: 10.1038/nrmicro2504
Li, T., and Gao, J. (2023). Attribution of dispersal limitation can better explain the assembly patterns of plant microbiota. Front. Plant Sci. 14:1168760. doi: 10.3389/fpls.2023.1168760
Li, G.-l., Wu, M., Li, P.-f., Wei, S.-p., Liu, J., Jiang, C.-y., et al. (2021). Assembly and co-occurrence patterns of rare and abundant bacterial sub-communities in rice rhizosphere soil under short-term nitrogen deep placement. J. Integr. Agric. 20, 3299–3311. doi: 10.1016/S2095-3119(20)63462-1
Liang, Y., Xiao, X., Nuccio, E. E., Yuan, M., Zhang, N., Xue, K., et al. (2020). Differentiation strategies of soil rare and abundant microbial taxa in response to changing climatic regimes. Environ. Microbiol. 22, 1327–1340. doi: 10.1111/1462-2920.14945
Liang, Y., Zhao, H., Deng, Y., Zhou, J., Li, G., and Sun, B. (2016). Long-term oil contamination alters the molecular ecological networks of soil microbial functional genes. Front. Microbiol. 7:60. doi: 10.3389/fmicb.2016.00060
Lin, G.-Y., Lin, C.-Y., Chang, S.-J., and And Lin, W.-Y. (2020). The dynamics of endophytic bacterial community structure in Rice roots under different field management systems. Agronomy 10:1623. doi: 10.3390/agronomy10111623
Liu, Y., Chu, G., Stirling, E., Zhang, H., Chen, S., Xu, C., et al. (2023b). Nitrogen fertilization modulates rice seed endophytic microbiomes and grain quality. Sci. Total Environ. 857:159181. doi: 10.1016/j.scitotenv.2022.159181
Liu, C., Cui, Y., Li, X., and Yao, M. (2020). Microeco: an R package for data mining in microbial community ecology. FEMS Microbiol. Ecol. 97, 1–9. doi: 10.1093/femsec/fiaa255
Liu, S., Liu, R., Zhang, S., Shen, Q., Chen, J., Ma, H., et al. (2024). The contributions of sub-communities to the assembly process and ecological mechanisms of bacterial communities along the cotton soil–root continuum niche gradient. Microorganisms 12:869. doi: 10.3390/microorganisms12050869
Liu, L., Yang, J., Yu, Z., and Wilkinson, D. M. (2015). The biogeography of abundant and rare bacterioplankton in the lakes and reservoirs of China. ISME J. 9, 2068–2077. doi: 10.1038/ismej.2015.29
Liu, L., Zhu, L., Yan, R., Yang, Y., Adams, J. M., and Liu, J. (2023a). Abundant bacterial subcommunity is structured by a stochastic process in an agricultural system with P fertilizer inputs. Sci. Total Environ. 871:162178. doi: 10.1016/j.scitotenv.2023.162178
Long, X.-E., and Yao, H. (2020). Phosphorus input alters the assembly of Rice (Oryza sativa L.) root-associated communities. Microb. Ecol. 79, 357–366. doi: 10.1007/s00248-019-01407-6
Lynch, M. D. J., and Neufeld, J. D. (2015). Ecology and exploration of the rare biosphere. Nat. Rev. Microbiol. 13, 217–229. doi: 10.1038/nrmicro3400
Ma, J., Ma, K., Liu, J., and Chen, N. (2022a). Rhizosphere soil microbial community under ice in a high-latitude wetland: different community assembly processes shape patterns of rare and abundant microbes. Front. Microbiol. 13:783371. doi: 10.3389/fmicb.2022.783371
Ma, L., Niu, W., Li, G., Du, Y., Sun, J., Zhang, Q., et al. (2024). Crucial role of rare taxa in preserving bacterial community stability. Land Degrad. Dev. 35, 1397–1410. doi: 10.1002/ldr.4994
Ma, Y., Wang, D., Guo, X., Zhu, Y.-G., Delgado-Baquerizo, M., and Chu, H. (2022b). Root stoichiometry explains wheat endophytes and their link with crop production after four decades of fertilization. Sci. Total Environ. 846:157407. doi: 10.1016/j.scitotenv.2022.157407
Ma, Y., Weisenhorn, P., Guo, X., Wang, D., Yang, T., Shi, Y., et al. (2021). Effect of long-term fertilization on bacterial communities in wheat endosphere. Pedosphere 31, 538–548. doi: 10.1016/S1002-0160(19)60838-6
Machado, P. C., Andrade, P. H. M., Sousa, C. P. D., Souza, C. W. O. D., and And Lacava, P. T. (2020). In vitro characterization of endophytic bacteria associated with physic nut (Jatropha curcas L.) and their potential for plant-growth promotion and biocontrol / Caracterização in vitro de bactérias endofíticas associadas ao pinhão-manso (Jatropha curcas L.) e seu potencial de promoção de crescimento vegetal e biocontrole. Braz. J. Dev. 6, 88572–88589. doi: 10.34117/bjdv6n11-326
Maranguit, D., Guillaume, T., and Kuzyakov, Y. (2017). Land-use change affects phosphorus fractions in highly weathered tropical soils. Catena 149, 385–393. doi: 10.1016/j.catena.2016.10.010
Mebius, L. J. (1960). A rapid method for the determination of organic carbon in soil. Anal. Chim. Acta 22, 120–124. doi: 10.1016/S0003-2670(00)88254-9
Mehta, P., Sharma, R., Putatunda, C., and Walia, A. (2019). “Endophytic fungi: role in phosphate solubilization” in Advances in Endophytic Fungal Research: Present Status and Future Challenges. ed. B. P. Singh (Cham: Springer International Publishing), 183–209.
Meng, Y., Zhang, Q., Shi, G., Liu, Y., Du, G., and Feng, H. (2022). Can nitrogen supersede host identity in shaping the community composition of foliar endophytic fungi in an alpine meadow ecosystem? Front. Microbiol. 13:895533. doi: 10.3389/fmicb.2022.895533
Michalska-Smith, M., Song, Z., Spawn-Lee, S. A., Hansen, Z. A., Johnson, M., May, G., et al. (2022). Network structure of resource use and niche overlap within the endophytic microbiome. ISME J. 16, 435–446. doi: 10.1038/s41396-021-01080-z
Miranda-Carrazco, A., Navarro-Noya Yendi, E., Govaerts, B., Verhulst, N., and Dendooven, L. (2022). Nitrogen fertilizer application alters the root endophyte bacterial microbiome in maize plants, but not in the stem or rhizosphere soil. Microbiol. Spectr. 10, e0178522–e0101722. doi: 10.1128/spectrum.01785-22
Nemergut Diana, R., Schmidt Steven, K., Fukami, T., O'Neill Sean, P., Bilinski Teresa, M., Stanish Lee, F., et al. (2013). Patterns and processes of microbial community assembly. Microbiol. Mol. Biol. Rev. 77, 342–356. doi: 10.1128/mmbr.00051-12
Ning, D., Yuan, M., Wu, L., Zhang, Y., Guo, X., Zhou, X., et al. (2020). A quantitative framework reveals ecological drivers of grassland microbial community assembly in response to warming. Nat. Commun. 11:4717. doi: 10.1038/s41467-020-18560-z
Niu, B., Paulson, J. N., Zheng, X., and Kolter, R. (2017). Simplified and representative bacterial community of maize roots. Proc. Natl. Acad. Sci. USA 114, E2450–e2459. doi: 10.1073/pnas.1616148114
Oksanen, J., Blanchet, F., Kindt, R., Legendre, P., Minchin, P., O'Hara, R., et al. (2012). Vegan: community ecology package: ordination, diversity and dissimilarities, version 2.0-4. (Reprinted from: Reprint E dition).
Omomowo, O. I., and Babalola, O. O. (2019). Bacterial and fungal endophytes: tiny giants with immense beneficial potential for plant growth and sustainable agricultural productivity. Microorganisms 7:481. doi: 10.3390/microorganisms7110481
Rani, S., Kumar, P., Dahiya, P., Maheshwari, R., Dang, A. S., and Suneja, P. (2022). Endophytism: a multidimensional approach to plant–prokaryotic microbe interaction. Front. Microbiol. 13:861235. doi: 10.3389/fmicb.2022.861235
Rognes, T. F. T., Nichols, B., Quince, C., and Mahé, F. (2016). VSEARCH: a versatile open source tool for metagenomics. PeerJ 4:e2584. doi: 10.7717/peerj.2584
Santoyo, G., Moreno-Hagelsieb, G., del Carmen Orozco-Mosqueda, M., and Glick, B. R. (2016). Plant growth-promoting bacterial endophytes. Microbiol. Res. 183, 92–99. doi: 10.1016/j.micres.2015.11.008
Schmidt, S. K., Nemergut, D. R., Darcy, J. L., and Lynch, R. (2014). Do bacterial and fungal communities assemble differently during primary succession? Mol. Ecol. 23, 254–258. doi: 10.1111/mec.12589
Sena, L., Mica, E., Valè, G., Vaccino, P., and Pecchioni, N. (2024). Exploring the potential of endophyte-plant interactions for improving crop sustainable yields in a changing climate. Front. Plant Sci. 15:1349401. doi: 10.3389/fpls.2024.1349401
Shi, Q., Song, Q., Shan, X., Li, X., Wang, S., Fu, H., et al. (2023). Microorganisms regulate soil phosphorus fractions in response to low nocturnal temperature by altering the abundance and composition of the pqqC gene rather than that of the phoD gene. Biol. Fertil. Soils 59, 973–987. doi: 10.1007/s00374-023-01766-w
Shu, D., Guo, Y., Zhang, B., Zhang, C., Van Nostrand, J. D., Lin, Y., et al. (2021). Rare prokaryotic sub-communities dominate the complexity of ecological networks and soil multinutrient cycling during long-term secondary succession in China's loess plateau. Sci. Total Environ. 774:145737. doi: 10.1016/j.scitotenv.2021.145737
Sindhu, L., Song, Y., Li, L., Dong, Z., Yang, Q., Mangi, S. S., et al. (2024). Meta-omics revealed that nitrogen fertilization alters the endophytic fungi divergence in maize field ecosystem. Ecol. Indic. 160:111852. doi: 10.1016/j.ecolind.2024.111852
Stegen, J. C., Lin, X., Fredrickson, J. K., and Konopka, A. E. (2015). Estimating and mapping ecological processes influencing microbial community assembly. Front. Microbiol. 6:370. doi: 10.3389/fmicb.2015.00370
Sun, A., Jiao, X.-Y., Chen, Q., Wu, A.-L., Zheng, Y., Lin, Y.-X., et al. (2021). Microbial communities in crop phyllosphere and root endosphere are more resistant than soil microbiota to fertilization. Soil Biol. Biochem. 153:108113. doi: 10.1016/j.soilbio.2020.108113
Sun, K., Lu, F., Huang, P.-W., Tang, M.-J., Xu, F.-J., Zhang, W., et al. (2022). Root endophyte differentially regulates plant response to NO3− and NH4+ nutrition by modulating N fluxes at the plant–fungal interface. Plant Cell Environ. 45, 1813–1828. doi: 10.1111/pce.14304
Sun, C., Wang, R., Tang, G., Cai, S., Shi, H., Liu, F., et al. (2023). Integrated 16S and metabolomics revealed the mechanism of drought resistance and nitrogen uptake in rice at the heading stage under different nitrogen levels. Front. Plant Sci. 14:1120584. doi: 10.3389/fpls.2023.1120584
Tian, K., Xing, Z., Kalkhajeh, Y. K., Zhao, T., Hu, W., Huang, B., et al. (2022). Excessive phosphorus inputs dominate soil legacy phosphorus accumulation and its potential loss under intensive greenhouse vegetable production system. J. Environ. Manag. 303:114149. doi: 10.1016/j.jenvman.2021.114149
Tiessen, H., and Moir, J. (1993). “Characterization of available P by sequential extraction” in Soil Sampling and Methods of Analysis, 75–86. doi: 10.1201/9781420005271.ch25
Walitang, D. I., Kim, K., Madhaiyan, M., Kim, Y. K., Kang, Y., and Sa, T. (2017). Characterizing endophytic competence and plant growth promotion of bacterial endophytes inhabiting the seed endosphere of Rice. BMC Microbiol. 17:209. doi: 10.1186/s12866-017-1117-0
Wan, W., Gadd, G. M., Yang, Y., Yuan, W., Gu, J., Ye, L., et al. (2021). Environmental adaptation is stronger for abundant rather than rare microorganisms in wetland soils from the Qinghai-Tibet plateau. Mol. Ecol. 30, 2390–2403. doi: 10.1111/mec.15882
Wang, S., Zheng, S., Bian, T., Wu, T., Li, X., Fu, H., et al. (2022). Photosynthetic characteristics combined with metabolomics analysis revealed potential mechanisms of cucumber (Cucumis sativus) yield reduction induced by different phosphorus stresses. Sci. Hortic. 302:111156. doi: 10.1016/j.scienta.2022.111156
Wu, J., Liu, W., Zhang, W., Shao, Y., Duan, H., Chen, B., et al. (2019). Long-term nitrogen addition changes soil microbial community and litter decomposition rate in a subtropical forest. Appl. Soil Ecol. 142, 43–51. doi: 10.1016/j.apsoil.2019.05.014
Wu, D., Ma, Y., Yang, T., Gao, G., Wang, D., Guo, X., et al. (2022). Phosphorus and zinc are strongly associated with belowground fungal communities in wheat field under Long-term fertilization. Microbiol. Spectr. 10:e0011022. doi: 10.1128/spectrum.00110-22
Wu, C., Yan, B., Wei, F., Wang, H., Gao, L., Ma, H., et al. (2023). Long-term application of nitrogen and phosphorus fertilizers changes the process of community construction by affecting keystone species of crop rhizosphere microorganisms. Sci. Total Environ. 897:165239. doi: 10.1016/j.scitotenv.2023.165239
Xi, M., Wang, Y., Yang, J., Bi, X., Zhong, S., Duan, T., et al. (2024). Spatial distribution and community composition of endophytic fungi within Mussaenda pubescens stems. Fungal Biol. 128, 1815–1826. doi: 10.1016/j.funbio.2024.05.001
Xie, H., Chai, Y., Liu, Z., Hao, W., and Gai, J. (2024). Community assembly of endophytic bacteria and fungi differs in soil-root continuum of Carex cepillacea. Appl. Soil Ecol. 194:105206. doi: 10.1016/j.apsoil.2023.105206
Xing, W., Gai, X., Ju, F., and Chen, G. (2023). Microbial communities in tree root-compartment niches under cd and Zn pollution: structure, assembly process and co-occurrence relationship. Sci. Total Environ. 860:160273. doi: 10.1016/j.scitotenv.2022.160273
Xiong, C., He, J. Z., Singh, B. K., Zhu, Y. G., Wang, J. T., Li, P. P., et al. (2020). Rare taxa maintain the stability of crop mycobiomes and ecosystem functions. Environ. Microbiol. 23, 1907–1924. doi: 10.1111/1462-2920.15262
Xu, R., Li, K., Zhang, W., Chen, A., Hou, H., Wei, W., et al. (2022). Responses of endophytic bacterial communities in rice roots to phosphorus deficiency at the seedling stages. Eur. J. Soil Biol. 110:103405. doi: 10.1016/j.ejsobi.2022.103405
Xue, M., Guo, Z., Gu, X., Gao, H., Weng, S., Zhou, J., et al. (2020). Rare rather than abundant microbial communities drive the effects of long-term greenhouse cultivation on ecosystem functions in subtropical agricultural soils. Sci. Total Environ. 706:136004. doi: 10.1016/j.scitotenv.2019.136004
Yu, Z., Zhang, W., Duo, P., Zhang, Y., Tian, X., Sun, L., et al. (2024). Assembly processes of rare and abundant taxa are closely related to the bacterial and fungal functionality during desert vegetation reestablishment. Pedosphere. doi: 10.1016/j.pedsph.2024.03.004
Zhang, H., Huang, B., Dong, L., Hu, W., Akhtar, M. S., and Qu, M. (2017). Accumulation, sources and health risks of trace metals in elevated geochemical background soils used for greenhouse vegetable production in southwestern China. Ecotoxicol. Environ. Saf. 137, 233–239. doi: 10.1016/j.ecoenv.2016.12.010
Zhang, J., Na, M., Wang, Y., Ge, W., Zhou, J., and Zhou, S. (2024). Cadmium levels and soil pH drive structure and function differentiation of endophytic bacterial communities in Sedum plumbizincicola: a field study. Sci. Total Environ. 912:168828. doi: 10.1016/j.scitotenv.2023.168828
Zhang, M., Zhang, X., Zhang, L., Zeng, L., Liu, Y., Wang, X., et al. (2021). The stronger impact of inorganic nitrogen fertilization on soil bacterial community than organic fertilization in short-term condition. Geoderma 382:114752. doi: 10.1016/j.geoderma.2020.114752
Zhang, X., Zhao, W., Kou, Y., Fang, K., Liu, Y., He, H., et al. (2023). The contrasting responses of abundant and rare microbial community structures and co-occurrence networks to secondary forest succession in the subalpine region. Front. Microbiol. 14:1177239. doi: 10.3389/fmicb.2023.1177239
Zhang, X., Zhong, Z., Gai, X., Du, X., Bian, F., Yang, C., et al. (2019). Changes of root endophytic bacterial community along a Chronosequence of intensively managed lei bamboo (Phyllostachys praecox) forests in subtropical China. Microorganisms 7:616. doi: 10.3390/microorganisms7120616
Zheng, W., Zhao, Z., Lv, F., Wang, R., Wang, Z., Zhao, Z., et al. (2021). Assembly of abundant and rare bacterial and fungal sub-communities in different soil aggregate sizes in an apple orchard treated with cover crop and fertilizer. Soil Biol. Biochem. 156:108222. doi: 10.1016/j.soilbio.2021.108222
Zhong, Y., Sorensen, P. O., Zhu, G., Jia, X., Liu, J., Shangguan, Z., et al. (2022). Differential microbial assembly processes and co-occurrence networks in the soil-root continuum along an environmental gradient. iMeta 1:e18. doi: 10.1002/imt2.18
Keywords: endophyte, community assembly, abundant microbial taxa, rare microbial taxa, solar greenhouse
Citation: Li X, Awais M, Wang S, Zhang Z, Zhao S, Liu Y, Sun Z, Fu H and Li T (2025) The process of nitrogen-adaptation root endophytic bacterial rather than phosphorus-adaptation fungal subcommunities construction unveiled the tomato yield improvement under long-term fertilization. Front. Microbiol. 15:1487323. doi: 10.3389/fmicb.2024.1487323
Edited by:
Tereza Castellane, São Paulo State University, BrazilReviewed by:
Everlon Cid Rigobelo, São Paulo State University, BrazilAgnieszka Kuźniar, The John Paul II Catholic University of Lublin, Poland
Poonam Singh, LaserLeap Technologies (Portugal), Portugal
Copyright © 2025 Li, Awais, Wang, Zhang, Zhao, Liu, Sun, Fu and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhouping Sun, c3VuenBAc3lhdS5lZHUuY24=; Hongdan Fu, ZnVob25nZGFuQHN5YXUuZWR1LmNu