The process of nitrogen-adaptation root endophytic bacterial rather than phosphorus-adaptation fungal subcommunities construction unveiled the tomato yield improvement under long-term fertilization.

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 , , and 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.