Soil nifH-harboring community assemblage varies across pecan cultivars

Introduction: This study focused on three pecan (Carya illinoinensis) cultivars (‘Pawnee’, ‘Mahan’, and ‘Jinhua’), systematically assessing variations in soil nitrogenase activity, characteristics traits of nifH-harboring microbial communities across these cultivars.

Methods: Using high-throughput sequencing technology, differences in the diversity, community composition, and network structure of nifH-harboring communities in the rhizosphere and bulk soils of pecan were examined across cultivars.

Results: Both cultivar type and soil compartment had significant effects on nitrogenase activity (p < 0.01). Among the three cultivars, ‘Mahan’ exhibited the highest soil nitrogenase activity in both its rhizosphere and bulk soils relative to the other two cultivars. Notably, rhizosphere soils across all cultivars displayed significantly stronger soil nitrogenase activity than their bulk soil counterparts. ‘Mahan’ harbored significantly higher microbial α-diversity (Sobs, Shannon, and Chao indices) than other cultivars (p < 0.05). β-diversity analysis revealed marked community divergence among cultivars, with the most pronounced differences observed in bulk soils. Specifically, the bulk soil of ‘Jinhua’ harbored a distinct microbial signature, marked by significant enrichment of Cyanobacteria and depletion of Alphaproteobacteria. Linear discriminant analysis effect size (LEfSe) further identified Rhizobiales and Burkholderiales as distinct biomarkers for the rhizosphere and bulk soils of ‘Jinhua’, respectively (LDA score > 4.0, p < 0.05). Microbial co-occurrence network analysis showed that the bulk soil of ‘Jinhua’ harbored the most complex microbial interaction network, characterized by the highest number of edges, and average connectivity. In co-occurrence networks, Azohydromonas, Bradyrhizobium, Azoarcus, Rhodomicrobium were found as the keystone taxa in maintaining network stability.

Discussion: This research elucidates the regulatory roles of pecan cultivars in shaping soil nitrogen fixation functions and microbial community assembly, providing valuable practical implications for precision microbiome management in pecan production.

1 Introduction

Nitrogen is a critical macronutrient essential for plant growth and development. However, its bioavailability in agricultural systems is frequently constrained by natural biogeochemical processes (Guo et al., 2023). Biological nitrogen fixation (BNF), mediated by specialized diazotrophic microorganisms, converts atmospheric dinitrogen (N₂) into biologically available ammonia. This process enhances soil fertility and reduces dependence on synthetic fertilizers in agricultural systems (Guo et al., 2023). The nifH gene, which encodes the iron protein subunit of nitrogenase, has been widely used as a reliable phylogenetic marker for investigating the diversity and distribution of diazotrophic communities in terrestrial ecosystems (Gaby John et al., 2018). Characterizing the diversity and abundance of soil nifH-harboring communities is essential for understanding the mechanisms of nitrogen cycling in ecosystems.

Pecan (Carya illinoinensis), a high-value nut tree cultivated worldwide in temperate regions, exhibits strong nitrogen-dependent growth responses that directly influence its yield potential (Wells, 2021; Pokhrel et al., 2025). While previous studies have extensively examined the agronomic effects of fertilization regimes on pecan productivity (Tong et al., 2024; Pokhrel et al., 2025), the functional role of indigenous diazotrophic communities associated with different cultivars remains poorly understood. It is well-documented that plant genotype can shape rhizosphere microbiome assembly through mechanisms such as root exudation (Rheault et al., 2020; Yang et al., 2022). For instance, studies in crops like sugarcane and cotton have demonstrated cultivar-specific differences in diazotrophic communities (Dong et al., 2018; Yang et al., 2022). Therefore, we hypothesize that distinct pecan cultivars may selectively enrich specific nifH-harboring microbial taxa, potentially influencing nitrogen fixation capacity.

Current knowledge of pecan-associated soil microbiomes has primarily focused on general bacterial and fungal communities (Ren et al., 2023; Ren et al., 2024), leaving a critical knowledge gap regarding diazotrophic assemblages. Specifically, the mechanisms underlying cultivar-specific recruitment of diazotrophic communities and their functional implications in pecan orchards have not been systematically investigated. Addressing this gap could provide insights into cultivar-specific microbial interactions that enhance nitrogen availability and reduce fertilizer dependency.

Specifically, we aimed to (1) characterize the diversity and composition of soil nifH-harboring microbial communities associated with different pecan cultivars and (2) quantitatively assess the soil nitrogen fixation potential of different pecan cultivars. We hypothesized that (i) distinct pecan cultivars harbor significantlybdifferent nifH-harboring microbial communities, and (ii) specific keystone microbial taxa drive cultivar-dependent variations in soil nitrogen fixation through deterministic community assembly processes. The findings are expected to provide a theoretical basis for developing microbiome-based strategies for sustainable pecan production, offering new insights for precision orchard management.

2 Materials and methods

2.1 Study site and soil sampling

This study was conducted in the Pecan Experimental Base of Nanjing Forestry University, Jurong city (J), Jiangsu Province, one of the main pecan cultivation regions in China with a subtropical humid monsoon climate and neutral-acid luvisol according to the soil taxonomy of the FAO. All experimental trees were consistent in key growth conditions: tree age (7 years), planting spacing (6 m × 8 m), management practices (uniform pruning, systematic weed control, and integrated pest control), and fertilization regime (annual application of organic compost). Rhizosphere (RS) and bulk soil (BS) samples were collected from three pecan cultivars, ‘Pawnee’ (P), ‘Mahan’ (M), and ‘Jinhua’ (J). Accordingly, six experimental groups were established: JPBS, JPRS, JMBS, JMRS, JJBS, and JJRS.

Six small plots (20 × 20 m²) were selected, with an inter-plot spacing of ≥ 50 m to avoid mutual interference. Each plot served as an independent biological replicate (n = 6). Within each plot, three uniformly distributed trees were chosen for soil mixture sampling (minimum spacing of 12 m between trees). From each tree, soil samples were collected from four directions (east, south, west, and north) and homogenized. Ultimately, each composite soil sample was composed of soil from 12 sampling points. Six composite soil samples were collected per group, with a total of 36 samples analyzed across all groups (total n = 36). Within a 0.5 m radius from the tree trunk, the surface humus was removed, and soil cores containing roots were collected using a soil auger. As pecan is a deep-rooted species, soil sampling was conducted at increased depths (0–60 cm) to ensure rhizosphere coverage. Following the identification of the pecan root system, rhizosphere soil was operationally defined as the soil firmly adhering to the root surface after gentle shaking to remove loose aggregates. This collected layer corresponds to soil within approximately 1–3 mm of the root surface. Concurrently, bulk soil was collected from areas located at least 10 cm away from any fine roots (0–60 cm depth). To avoid rhizosphere influence, the absence of roots was verified by visual inspection prior to sampling. All collected samples were immediately sealed in sterile polyethylene bags, clearly labeled, and stored in dry ice-cooled insulated containers to maintain biological integrity during transport to the laboratory. Upon arrival, each soil sample was sieved (2 mm mesh) and subdivided under aseptic conditions. One aliquot for DNA extraction and molecular analysis was stored at −80 °C. The other aliquot for soil physicochemical characterization was stored at 4 °C. prior to sampling, all tools were sterilized to eliminate potential microbial contamination.

2.2 Determination of soil nitrogenase activity

Soil biological nitrogen fixation (BNF) activity was assessed by measuring nitrogenase activity using the acetylene reduction assay (ARA), following a modified protocol based on the approach reported from published research (Qi et al., 2022). Briefly, 10 g of fresh soil (equivalent to 8.1 ± 0.06 g dry soil) was placed in a 20 mL penicillin bottle. All samples were maintained at their natural moisture content during the pre-incubation phase without additional water supplementation, but comparability within the experimental group was ensured by subjecting all samples to identical pre-treatment procedures. Samples were pre-incubated at 30 °C in the dark for 24 h with perforated breathable film to allow microbial acclimation. Following pre-incubation, the parafilm was replaced with a butyl rubber stopper and an aluminum-plastic cap to ensure airtight conditions. Using a syringe, 1 mL of gas was extracted from the bottle, and 1 mL of acetylene (C₂H₂) was injected (headspace volume: 10 mL). All bottles were gently shaken to ensure uniform gas distribution before incubation. Blank vials without soil substrate were included as negative controls in each batch to account for any ethylene background. Samples were then incubated at 30 °C in the dark under static conditions for 48 h. After incubation, 1 mL of headspace gas was sampled to measure the relative concentration of ethylene (C₂H₄).

Gas chromatography (GC) analysis was performed under the following optimized conditions: flame ionization detector (FID) temperature = 220 °C; capillary column (HP-PLOT Q, 30 m × 0.32 mm × 20 μm) temperature = 80 °C (isothermal); carrier gas = high-purity nitrogen (N₂, 99.999%) with a flow rate of 1.0 mL/min; fuel gas = high-purity hydrogen (H₂, 99.999%); auxiliary gas = compressed air. GC calibration was conducted using a series of C₂H₄ standard solutions (0.01–10 nmol/mL) to establish a standard curve (R² ≥ 0.999), with a detection limit of 0.01 nmol. Each sample was injected in triplicate to ensure analytical precision. Nitrogenase activity was expressed as the amount of C₂H₄ produced per gram of fresh soil per hour (nmol C₂H₄ g⁻¹ fresh soil h⁻¹).

2.3 High-throughput sequencing of nifH-harboring community

Total DNA was extracted from 0.5 g of soil using the FastDNA SPIN Kit for Soil (MP Biomedicals, CA, United States) following the manufacturer’s protocol. Extraction blanks (negative controls) were included in parallel to monitor potential contamination during the DNA extraction process. DNA integrity was assessed by 1% agarose gel electrophoresis, while purity and concentration were determined using a NanoDrop 2000 spectrophotometer (Thermo Scientific, Wilmington, United States). The nifH gene of nitrogen-fixing microorganisms was amplified using the primers: nifHF: 5’-AAAGGYGGWATCGGYAARTCCACCAC-3′, nifHR: 5’-TTGTTSGCSGCRTACATSGCCATCAT-3′. PCR conditions and amplification procedures followed previously described methods (Fan et al., 2019). Each sample was amplified in triplicate, and PCR products from the same sample were pooled. PCR negative controls (with no template) were also performed to ensure the absence of contamination during amplification. The target fragments were excised from a 2% agarose gel and purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, United States). Purified amplicons were sequenced on an Illumina MiSeq PE300 platform (Illumina, San Diego, United States) at Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China). The raw sequencing data have been deposited in the National Center for Biotechnology Information (NCBI), Sequence Read Archive (SRA) under accession number PRJNA907477.

Raw paired-end reads were first demultiplexed, quality-filtered using fastp (v0.20.0), and merged into contiguous sequences with FLASH (v1.2.7). Subsequently, operational taxonomic units (OTUs) were clustered at a 97% sequence similarity threshold using UPARSE (v7.1), with chimeric sequences identified and removed during this process. OTUs that aligned to chloroplast or mitochondrial sequences were further excluded to eliminate non-target contaminants. Following this, all samples were rarefied to the minimum number of sequences observed across the dataset to standardize sequencing depth. The cleaned sequence data were then aligned against the nifH functional gene database (Version 7.3)¹ using a 70% alignment threshold for taxonomic and functional annotation. This confidence threshold was selected as it is widely adopted in nifH gene studies to balance classification precision with the retention of sufficient sequences for ecological analysis, given the comparatively smaller and less comprehensive nature of the nifH gene reference database. Prior to downstream statistical analysis, the sequencing data were subjected to minimum sequence normalization to ensure consistency across samples.

2.4 Statistical analysis

After rarefaction, the coverage index of all samples reached above 95% (Supplementary Figure 1), indicating that the sequencing depth was sufficient to capture the majority of the diazotrophic microbial composition and could be used for subsequent diversity analyses. Because the data did not meet the assumptions for parametric tests, the Scheirer–Ray–Hare test was used to evaluate the effects of cultivar., niche, and their interaction on α-diversity. Post–hoc pairwise comparisons were performed using Kruskal-Wallis test, with Bonferroni correction applied to control the false-discovery rate in multiple testing. NMDS was carried out with 999 random starts using the vegan package (v2.6–6) to ensure the stability of the ordination. A two-factor PERMANOVA (adonis2 function, 999 permutations) was conducted to test the overall effects of cultivar., niche, and their interaction on community structure. Homogeneity of multivariate dispersions was tested with the betadisper function.

To evaluate differences in soil nitrogenase activity, analysis of two-way variance (ANOVA) and paired t-tests were applied, with nitrogenase activity further visualized using GraphPad Prism 8. For microbial co-occurrence analysis, Spearman’s correlation was calculated between OTUs, and all correlations were adjusted for multiple testing with the False Discovery Rate (FDR) correction to control false positives. Statistically robust correlations were used to construct co-occurrence networks, which were subsequently analyzed with the igraph package (v0.9.2). It should be noted that these networks reflect inferred associations based on correlation statistics; while they suggest potential ecological linkages, they do not constitute direct evidence of biological interaction or causality. Additional bioinformatics steps, including sequence quality control, OTU clustering, taxonomic annotation, compositional analysis, and visualization were performed on the Majorbio Cloud Platform.²

3 Results

3.1 Soil nitrogenase activity of different pecan cultivars

We found that both cultivar and soil compartment significantly affected nitrogenase activity (p < 0.01; Table 1). JMBS showed significantly higher nitrogenase activity than JJBS and JPBS (p < 0.05), whereas no significant difference was observed between JPBS and JJBS. JMRS exhibited significantly greater activity than JJRS and JPRS (p < 0.05), JPRS and JJRS were statistically indistinguishable. Notably, all three cultivars displayed significantly elevated nitrogenase activity in RS compared to BS (Figure 1). In summary, the ‘Mahan’ cultivar showed higher nitrogenase activity in both rhizosphere and bulk soils compared to other cultivars, with consistently greater activity observed in rhizosphere soil versus bulk soil.

Table 1

Table 1. Two-way ANOVA results of nitrogenase activity in rhizosphere soil and bulk soil of pecans of different varieties.

Figure 1

Figure 1. Nitrogenase activity in rhizosphere soil and bulk soil of pecans of different varieties. Boxplot elements: center line = median; box limits = Q3/Q1; whiskers = range. Significant differences (Tukey’s HSD after ANOVA, p < 0.05) are indicated by different lowercase letters.

3.2 Soil nifH-harboring community diversity of different pecan cultivars

High-throughput sequencing of the nifH-harboring community generated 849,512 high-quality sequences (344,954,731 bases; mean length = 406 bp) from rhizosphere and bulk soil samples across pecan cultivars (Supplementary Table 1). Following rarefaction to 8,413 sequences per sample, sequences clustered into 3,096 OTUs, representing 103 species, 64 genera, 47 families, 32 orders, 16 classes, and 9 phyla. The Scheirer-Ray-Hare test revealed significant differences in Sobs, Shannon, Simpson, and Chao indices among pecan varieties (p ≤ 0.05; Table 2), while no significant differences were observed between rhizosphere and bulk soils. Post-hoc comparisons demonstrated that: JM exhibited significantly higher Sobs and Chao indices than JJ, as well as a significantly higher Shannon index than JP (Table 3). Collectively, these α-diversity metrics suggest that pecan varieties exerted a stronger influence on the diversity of soil diazotrophic communities than ecological niche, with the ‘Mahan’ variety consistently demonstrating the highest diversity.

Table 2

Table 2. The Scheirer-Ray-Hare test results of α-diversity indexes of nifH-harboring community in rhizosphere soil and bulk soil of pecans of different varieties.

Table 3

Table 3. Multiple comparison results of the soil nifH-harboring community indexes in soil of pecans of different varieties.

Non-metric multidimensional scaling (NMDS) revealed highly significant differences in β-diversity among varieties (p < 0.01; Figure 2). Along NMDS2, JJBS samples were distinctly separated from JPBS and JMRS, while NMDS1 showed weaker separation. These results underscore those structural differences in nifH-harboring communities were predominantly associated with bulk soil rather than rhizosphere soil across pecan varieties.

Figure 2

Figure 2. NMDS analysis of the nifH-harboring community of in rhizosphere soil and bulk soil of pecans of different varieties. NMDS analysis based on Bray-Curtis distance. Different colors represent distinct sample groups. The closer two sample points are, the more similar their species composition. The x-axes and y-axes represent relative distances and have no actual units. Stress value: This metric evaluates the reliability of the NMDS ordination. Generally: Stress < 0.2: The NMDS result has reasonable explanatory power. Stress < 0.1: Indicates an excellent ordination with high confidence.

3.3 Composition and biomarker taxa of the soil nifH-harboring community of different pecan cultivars

As presented in Figure 3, Proteobacteria dominant the nifH-harboring community. Notably, JJBS exhibited a distinct microbial profile characterized by significantly lower abundance of Rhizobiales and higher representation of Cyanobacteria and Burkholderiales compared to other groups (Figure 3A). At the order level, Rhizobiales, Desulfuromonadales, Rhodospirillales, Burkholderiales, Rhodocyclales, Nostocales, Methylococcales, and Oscillatoriales were identified as the dominant bacterial orders, with JJBS demonstrating markedly different community structure from other groups (Figure 3B). Kruskal-Wallis H tests of predominant families revealed significant variations among cultivars. JMBS and JMRS showed significantly higher relative abundances of Myxococcaceae, Beijerinckiaceae, and unclassified Betaproteobacteria compared to other cultivars (p < 0.05). JJBS exhibited the highest relative abundance of Alcaligenaceae among all groups, while displaying the lowest relative abundances of unclassified Rhizobiales and Hyphomicrobiaceae (p < 0.05; Figure 4). LEfSe multilevel species discriminant analysis further identified Rhizobiales and Burkholderiales as biomarker taxa for JJRS and JJBS, respectively (LDA score >4.0, p < 0.05; Table 4).

Figure 3

Figure 3. The nifH-harboring community compositions of rhizosphere soil and bulk soil of pecans in phylum (A) and order (B) levels of pecans of different varieties. Taxonomic groups with relative abundances <0.01 were merged and classified as others.

Figure 4

Figure 4. Differences of dominant nifH-harboring taxa in rhizosphere soil and bulk soil of pecans of different varieties. Kruskal-Wallis test with FDR correction were used for assessing significance of variety differences. Only the top 10 diazotrophic families in relative abundance are shown (*p < 0.05).

Table 4

Table 4. Biomarkers (at order level) in rhizosphere soil and bulk soil of pecans of different varieties.

3.4 Network structure of the soil nifH-harboring community in different pecan cultivars

Network analysis of the top 500 most abundant OTUs revealed distinct topological patterns among pecan cultivars. Both rhizosphere soil and bulk soil of the ‘Jinhua’ cultivar exhibited significantly more complex networks compared to the other two cultivars (Figure 5; Table 5). Specifically, the JJBS network demonstrated the highest total connectivity and average node degree, along with the lowest modularity index, indicating exceptionally tight clustering of bacterial OTUs. These topological characteristics collectively contributed to the increased complexity of the network in JJBS. Comparative analysis showed that rhizosphere soil across all three cultivars consistently formed simpler networks than their bulk soil counterparts, as evidenced by significantly fewer edges, lower node degrees, and reduced overall network complexity (Figure 5; Table 5). Keystone species analysis based on betweenness centrality identified the top five most influential taxa in each group (abundance range: 7–310; Table 6). Taxonomic classification of these keystone taxa revealed distinct patterns across groups: In the JPRS group, the dominant taxa included unclassified Rhizobiales, Azohydromonas, and Bradyrhizobium. The JMBS group was characterized by Azoarcus and unclassified Betaproteobacteria, while JMRS showed a predominance of unclassified Rhizobiales. For the JJBS group, Rhodomicrobium and unclassified Alphaproteobacteria were prominent, and the JJRS group featured unclassified Rhizobiales and Alphaproteobacteria as key taxa (Table 6).

Figure 5

Figure 5. Correlation network analysis of the nifH-harboring community in rhizosphere soil and bulk soil of pecans of different varieties. Network of the top 500 OTUs by total abundance. Nodes (OTUs) are colored by phylum and sized by degree. Edges represent strong (Spearman’s ρ > 0.8) and significant (*p < 0.01) correlations (red: positive; green: negative). Sample size n = 36, covering 3 pecan varieties × 2 ecological niche × 6 replicates.

Table 5

Table 5. Topological properties of correlation network structure of the nifH-harboring community in rhizosphere soil and bulk soil of pecans of different varieties.

Table 6

Table 6. Keystone taxa in correlation network of the nifH-harboring community in rhizosphere soil and bulk soil of pecans of different varieties.

4 Discussion

4.1 Cultivar-specific regulation of nitrogen fixation potential and rhizosphere effects

The observed variation in nitrogenase activity among pecan cultivars suggested a potential role of root exudates in modulating the recruitment of free-living diazotroph. This interpretation aligned with established mechanisms in plant-microbe interactions. For instance, plant-secreted flavonoids were shown to enhance the metabolic activity of nitrogen-fixing microorganisms (Yu et al., 2021a; Yan et al., 2022). Similarly, studies in poplar systems demonstrated that host genetic regulation of flavonoid biosynthesis influenced rhizosphere microbiome assembly and promoted nitrogen assimilation through selective microbial enrichment (Wu et al., 2025). It was also recognized that the synthesis of secondary metabolites, including putative signaling molecules, varies according to plant genotype, developmental stage, and physiological status (Bag et al., 2022). Collectively, these findings supporte the hypothesis that cultivar-specific root chemistry mediates microbial nitrogen fixation efficiency by shaping the recruitment of functional microbes. However, direct evidence identifying flavonoids or other exudates as the primary signaling mediators awaits further experimental confirmation.

The energy-intensive nature of biological nitrogen fixation necessitates high-carbon environments (Smercina Darian et al., 2019). Root exudates, rich in low-molecular-weight carbon compounds (sugars, organic acids, and polysaccharides), provide ideal substrates to meet diazotrophic energetic demands (Smercina Darian et al., 2019). Crucially, recent studies suggested carbon substrate quality may be equally important as quantity in regulating fixation rates, with both root exudates and litter inputs potentially driving variation in nitrogen fixation capacity (Wang et al., 2023). In our study, the ‘Mahan’ cultivar demonstrated superior nitrogen fixation potential in both rhizosphere and bulk soils (Figure 1). Analysis of soil chemical properties revealed that the total carbon content in the rhizosphere soil of this cultivar was significantly higher than in other cultivars (Supplementary Table 2). This finding provides a direct mechanistic explanation for the enhanced nitrogen fixation capacity of the ‘Mahan’ cultivar from the perspective of carbon supply, suggesting that it may allocate more carbon resources to the soil through root exudates and litter inputs.

Notably, a consistent rhizosphere effect was observed across all pecan cultivars, as evidenced by significantly higher nitrogen fixation capacity in the rhizosphere than in bulk soils (Figure 1). This functional enhancement was further supported by the elevated total nitrogen (TN) content specifically in the rhizosphere of the ‘Mahan’ and ‘Jinhua’ cultivars (Supplementary Table 2), indicating not only increased nitrogen fixation activity but also actual nitrogen accumulation in the root-associated environment. The rhizosphere effect is likely driven by root-mediated processes such as localized nutrient inputs and the creation of unique physicochemical conditions that facilitate the establishment and function of nitrogen-fixing microorganisms. While previous studies have suggested that factors such as pH, phosphorus availability, micronutrients, and oxygen tension can influence nitrogen fixation (Smercina Darian et al., 2019), the present study provides direct evidence linking rhizosphere properties to both the process and outcome of biological nitrogen fixation in pecan cultivars.

4.2 Assembly of nifH-harboring communities: diversity, composition, and network complexity

4.2.1 Diversity and community differentiation of nifH-harboring communities among cultivars and niches

Our results revealed distinct cultivar-dependent patterns in nifH-harboring community assembly. The ‘Mahan’ cultivar demonstrated consistently higher α-diversity and nitrogen fixation potential in both rhizosphere and bulk soils (Figure 1; Table 3). We propose two putative mechanisms underlying this observed advantage:

1. Root exudate-mediated microbial recruitment.

Host plants can selectively enrich specific symbiotic diazotroph through root exudates. We speculate that ‘Mahan’ may secrete particular chemical signaling molecules such as flavonoids and benzoic acids, thereby promoting the diazotroph recruitment and activation. The mechanism has been previously documented in maize and sugarcane (Liu et al., 2020; Yu et al., 2021a). Notably, even in bulk soil outside the rhizosphere, ‘Mahan’ maintained higher nitrogen fixation capacity and nifH-harboring community diversity (Figure 1; Table 3), suggesting that its root exudates may possess stronger diffusion effects or be capable of inducing a microbial “memory effects” (Haichar et al., 2014; Sasse et al., 2018; Oyserman et al., 2022).

1. Keystone taxa-mediated niche construction contributing to functional stability.

Soils associated with ‘Mahan’ were markedly enriched in keystone taxa from Myxococcaceae and Beijerinckiaceae (Figure 4), which are known for their dual capabilities in nitrogen fixation and complex carbon metabolism (Borken et al., 2016; Wegner et al., 2019). We propose that this enrichment underpins efficient nitrogen fixation by securing a steady energy supply: the decomposition of complex carbon compounds by these microbes likely provides the necessary substrates and energy (e.g., carbon skeletons and ATP) to support the nitrogenase activity.

4.2.2 Dominant taxa and biomarkers in nifH-harboring communities

The predominance of Proteobacteria confirms their established role as primary nitrogen-fixing agents in terrestrial ecosystems (Ulrich et al., 2018; Meng et al., 2019; Zhou et al., 2021; Li et al., 2024). However, cultivar-specific partitioning emerged:

Notably, JJBS and JJRS recruited distinct nitrogen-fixing microbiomes, leading to divergent functional outcomes. JJBS exhibited a nifH-harboring community enriched with Cyanobacteria but depleted in Rhizobiales relative to JJRS. This shift toward cyanobacteria-mediated pathways may have compromised efficiency, as their photosynthetic oxygen inactivates nitrogenase while their growth intensifies competition for nutrients (Nawaz et al., 2025). Concurrently, JJRS, with its community favoring α-Proteobacteria (e.g., Rhizobiales) for symbiotic fixation, demonstrated superior nitrogenase activity. The key takeaway is that cultivar-specific microbial partitioning is a primary determinant of nitrogen fixation performance. This niche differentiation may be driven by cultivar-specific root exudates or microenvironmental conditions (Hartmann et al., 2009; Yu et al., 2021b; Ulbrich et al., 2022). In this study, the abundance of Alcaligenaceae in JJBS was significantly higher (Figure 4), which may reflect specific environmental selection pressure or host physiological state. Alcaligenaceae, typically reported as plant growth-promoting bacteria (PGPB; Abdel Latef et al., 2021), are known for heterotrophic nitrogen fixation and stress tolerance, preferring neutral to weakly alkaline environments, with abundance often increasing in high-nitrate conditions (Kalniņš et al., 2020).

In contrast, JMBS and JMRS specifically enriched for Myxococcaceae and Beijerinckiaceae (Figure 4). The enrichment of these two groups suggests that they may jointly participate in regulating soil nitrogen cycling-potentially through organic matter degradation (Myxococcaceae) or free-living nitrogen fixation driven by organic carbon (Beijerinckiaceae; Ulrich et al., 2018; Kalniņš et al., 2020; Wolińska et al., 2018; Carrascosa-Robles et al., 2024; Masuda et al., 2024).

LEfSe analysis further confirmed Rhizobiales and Burkholderiales are the most critical discriminatory biomarker taxa for JJRS and JJBS, respectively (Table 4). In particular, Burkholderiales have been recognized as signature nifH-harboring community members (Zhou et al., 2021), playing important roles in motility and nitrogen metabolism (Berkelmann et al., 2020). Many associative diazotroph within Burkholderiales prefer slightly acidic environments (Dludlu et al., 2018), which is closely related to N content (Zhou et al., 2021). Additionally, many members of the Burkholderiales group are renowned for aromatic compound biodegradation potential (Pérez-Pantoja et al., 2012).

Collectively, this evidence points to a clear ecological conclusion: different pecan cultivars adopt distinct nitrogen acquisition strategies by shaping specific rhizosphere microbial communities. JJRS relies on efficient symbiotic nitrogen fixation (Rhizobiales), whereas JJBS shifts toward associative or free-living nitrogen fixation pathways driven by Cyanobacteria and Burkholderiales, likely as an adaptive strategy to low-nitrogen soils. These cultivar-dependent differences in microbial recruitment hold important implications for guiding cultivar selection and formulating precision nutrient management practices. Whether these differences ultimately translate into variations in field-level nitrogen use efficiency warrants further investigation.

4.2.3 Microbial co-occurrence networks and keystone species

Notably, we observed that diazotrophic co-occurrence networks in rhizosphere soils had lower topological features (edge number, node degree) than those in bulk soils (Table 5; Figure 5). This simplification aligns with the established paradigm that strong host plant selection filters rhizosphere communities, reducing complexity and enriching core taxa. Our novel contribution is revealing how this pattern diverges at the cultivar level. Specifically, the exceptionally complex co-occurrence network in JJBS suggests cultivar-specific recruitment of a more interactive microbial consortium (Li et al., 2023), representing a distinct nitrogen acquisition strategy compared to JJRS.

Proteobacteria emerged as keystone species in the microbial networks, primarily due to their role in organic substrate mineralization, which supports the nitrogen fixation niche (Xiao et al., 2022). Their network-stabilizing role highlighted the importance of keystone species in regulating cross-kingdom interactions and maintaining nitrogen fixation under environmental fluctuations (Gao et al., 2021). In the JPRS network, Azohydromonas and Bradyrhizobium occupied central positions (Table 6), indicating a potential direct linke to host plant nitrogen acquisition efficiency. Bradyrhizobium is widely recognized for its symbiotic relationship with soybean (Sarao et al., 2025), yet it also functions as a dominant free-living diazotroph (Meng et al., 2019; Zhou et al., 2021; Li et al., 2022). Similarly, Azohydromonas was commonly detected prevalence in nifH-harboring community (Xie et al., 2021; Tang et al., 2023).

We identified Azoarcus and Rhodomicrobium as network hubs in JMBS and JJBS, respectively, (Table 6). Azoarcus is a well-characterized free-living diazotroph (Sarao et al., 2025), can form associative nitrogen-fixing relationships with plants to promote growth (Sakoda et al., 2019), and has been found enriched in rice endophytic communities (Wartiainen et al., 2008). The prominence of these non-symbiotic diazotrophs in ecological networks may reflect the ecological relevance of chemoheterotrophic nitrogen fixation pathways in bulk soils. Furthermore, the frequent detection of unclassified Rhizobiales as core species indicates certain limitations in current taxonomic databases for the functional interpretation of environmental microbes. For instance, this may lead to provisional functional annotations or obscure the specific ecological functions of certain uncultivated taxa. These observations further underscore the need to integrate multi-omics approaches in future research to more comprehensively elucidate the ecological roles and functional potential of these uncultivated microbial lineages.

4.3 Implications for pecan cultivation and future research directions

These findings provide novel insights for cultivar selection from the perspective of microbiome-mediated nitrogen acquisition optimization. Subsequent studies will integrate metabolomics with gnotobiotic cultivation experiments to establish causal relationships between root exudates and diazotroph recruitment. Our study demonstrated that cultivar selection modulated nitrogen fixation potential through microbial community restructuring. For pecan cultivation, promoting genotypes with strong rhizosphere effects (e.g., ‘Mahan’) may enhance soil nitrogen availability while reducing dependence on synthetic fertilizers.

However, the functional significance of taxonomic biomarkers requires further validation through metagenomics and stable isotope probing techniques to link community structure with in situ nitrogen fixation rates. Future investigations should also examine: (1) long-term interactions between cultivar-driven microbial communities and soil physicochemical properties, and (2) the feasibility of inoculating superior diazotroph to enhance nitrogen fixation capacity in low-performance cultivars.

In conclusion, this study revealed cultivar-specific strategies for recruiting and maintaining diazotrophic communities in pecan systems, with significant implications for precision agriculture and microbial resource management. The integration of plant genotype selection with soil microbiome engineering may offer a promising pathway for optimizing nutrient cycling and advancing sustainable pecan production.

5 Conclusion

This study demonstrates that pecan cultivars play a decisive role in shaping the structure of rhizosphere diazotrophic communities. ‘Mahan’ cultivar specifically enriched key nitrogen-fixing taxa, including Bradyrhizobium and Azoarcus, and reinforced the interaction stability of the microbial network. These findings provide a mechanistic, microbe-mediated explanation for the enhanced nitrogen-cycling potential associated with superior cultivars. Our results suggest that selecting cultivars with a high capacity for recruiting beneficial diazotrophs could serve as an ecologically informed strategy to reduce dependence on synthetic nitrogen fertilizers. Future work should prioritize elucidating the molecular drivers, such as cultivar-specific root exudate profiles, that govern microbial assembly, and conducting long-term field trials to validate whether the identified microbial biomarkers can reliably predict in-situ nitrogen fixation capacity. Such efforts will provide a robust scientific basis for microbiome-based management in pecan cultivation.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Author contributions

JL: Conceptualization, Data curation, Methodology, Visualization, Writing – original draft, Writing – review & editing. HW: Data curation, Writing – review & editing. YT: Data curation, Writing – review & editing. JB: Data curation, Writing – review & editing. PT: Funding acquisition, Supervision, Writing – review & editing. FP: Conceptualization, Funding acquisition, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was funded by the National Promotion Project for Forestry and Grassland Sci-Tech Achievements (2023133126), Science and Technology Innovation Project of Jiangxi Forestry Bureau (Innovation Special Project [2022] No. 39), National Key Research and Development Project of China (2021YFD1000403), Doctoral Scientific Research Start-up Fund of Jiujiang University (A00624).

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that Generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors 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.2025.1716240/full#supplementary-material

Footnotes

^http://fungene.cme.msu.edu/

^https://cloud.majorbio.com/

References

Abdel Latef, A. A., Omer, A. M., Badawy, A. A., Osman, M. S., and Ragaey, M. M. (2021). Strategy of salt tolerance and interactive impact of Azotobacter chroococcum and/or Alcaligenes faecalis inoculation on canola (Brassica napus L.) plants grown in saline soil. Plants 10:110. doi: 10.3390/plants10010110,

PubMed Abstract | Crossref Full Text | Google Scholar

Bag, S., Mondal, A., Majumder, A., Mondal, S. K., and Banik, A. (2022). Flavonoid mediated selective cross-talk between plants and beneficial soil microbiome. Phytochem. Rev. 21, 1739–1760. doi: 10.1007/s11101-022-09806-3,

PubMed Abstract | Crossref Full Text | Google Scholar

Berkelmann, D., Schneider, D., Meryandini, A., and Daniel, R. (2020). Unravelling the effects of tropical land use conversion on the soil microbiome. Environ. Microbiome 15:5. doi: 10.1186/s40793-020-0353-3,

PubMed Abstract | Crossref Full Text | Google Scholar

Borken, W., Horn, M. A., Geimer, S., Aguilar, N.a. B., and Knorr, K.-H. (2016). Associative nitrogen fixation in nodules of the conifer Lepidothamnus fonkii (Podocarpaceae) inhabiting ombrotrophic bogs in southern Patagonia. Sci. Rep. 6:39072. doi: 10.1038/srep39072,

PubMed Abstract | Crossref Full Text | Google Scholar

Carrascosa-Robles, Á., Pascual, J. A., Cuartero, J., De Santiago, A., Petropoulos, S. A., and Del Mar Alguacil, M. (2024). Optimizing purslane cultivation through legume intercropping and crop rotation: a study on yield and rhizosphere bacterial communities. Plant Soil. doi: 10.1007/s11104-024-07061-3,

PubMed Abstract | Crossref Full Text | Google Scholar

Dludlu, M. N., Chimphango, S. B. M., Stirton, C. H., and Muasya, A. M. (2018). Differential preference of Burkholderia and Mesorhizobium to pH and soil types in the Core cape subregion, South Africa. Gene 9:8. doi: 10.3390/genes9010002,

PubMed Abstract | Crossref Full Text | Google Scholar

Dong, M., Yang, Z., Cheng, G., Peng, L., Xu, Q., and Xu, J. (2018). Diversity of the bacterial microbiome in the roots of four Saccharum species: S. spontaneum, S. robustum, S. Barberi, and S. officinarum. Front. Microbiol. 9:267. doi: 10.3389/fmicb.2018.00267,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

Gaby John, C., Rishishwar, L., Valderrama-Aguirre Lina, C., Green Stefan, J., Valderrama-Aguirre, A., Jordan, I. K., et al. (2018). Diazotroph community characterization via a high-throughput nifH amplicon sequencing and analysis pipeline. Appl. Environ. Microbiol. 84, e01512–e01517. doi: 10.1128/AEM.01512-17,

PubMed Abstract | Crossref Full Text | Google Scholar

Gao, Y., Du, J., Bahar, M. M., Wang, H., Subashchandrabose, S., Duan, L., et al. (2021). Metagenomics analysis identifies nitrogen metabolic pathway in bioremediation of diesel contaminated soil. Chemosphere 271:129566. doi: 10.1016/j.chemosphere.2021.129566,

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, K. Y., Yang, J., Yu, N., Luo, L., and Wang, E. R. (2023). Biological nitrogen fixation in cereal crops progress, strategies, and perspectives. Plant Commun. 4:100499. doi: 10.1016/j.xplc.2022.100499,

PubMed Abstract | Crossref Full Text | Google Scholar

Haichar, F. E. Z., Santaella, C., Heulin, T., and Achouak, W. (2014). Root exudates mediated interactions belowground. Soil Biol. Biochem. 77, 69–80. doi: 10.1016/j.soilbio.2014.06.017

Crossref Full Text | Google Scholar

Hartmann, A., Schmid, M., Tuinen, D. V., and Berg, G. (2009). Plant-driven selection of microbes. Plant Soil 321, 235–257. doi: 10.1007/s11104-008-9814-y

Crossref Full Text | Google Scholar

Kalniņš, M., Bērziņš, A., Gudrā, D., Megnis, K., Fridmanis, D., Danilko, P., et al. (2020). Selective enrichment of heterotrophic nitrifiers Alcaligenaceae and Alcanivorax spp. from industrial wastewaters. AIMS Microbiol. 6, 32–42. doi: 10.3934/microbiol.2020002,

PubMed Abstract | Crossref Full Text | Google Scholar

Li, L. Z., Hu, Z. R., Tan, G., Fan, J. Q., Chen, Y. Q., Xiao, Y. S., et al. (2023). Enhancing plant growth in biofertilizer-amended soil through nitrogen-transforming microbial communities. Front. Plant Sci. 14:1259853. doi: 10.3389/fpls.2023.1259853,

PubMed Abstract | Crossref Full Text | Google Scholar

Li, K. X., Lin, H. M., Han, M., and Yang, L. M. (2024). Soil metagenomics reveals the effect of nitrogen on soil microbial communities and nitrogen-cycle functional genes in the rhizosphere of Panax ginseng. Front. Plant Sci. 15:1411073. doi: 10.3389/fpls.2024.1411073,

PubMed Abstract | Crossref Full Text | Google Scholar

Li, L. S., Xia, T. Z., and Yang, H. Q. (2022). Seasonal patterns of rhizosphere microorganisms suggest carbohydrate-degrading and nitrogen-fixing microbes contribute to the attribute of full-year shooting in woody bamboo. Front. Microbiol. 13:1033293. doi: 10.3389/fmicb.2022.1033293,

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, H., Brettell, L. E., Qiu, Z., and Singh, B. K. (2020). Microbiome-mediated stress resistance in plants. Trends Plant Sci. 25, 733–743. doi: 10.1016/j.tplants.2020.03.014,

PubMed Abstract | Crossref Full Text | Google Scholar

Masuda, Y., Mise, K., Xu, Z., Zhang, Z., Shiratori, Y., Senoo, K., et al. (2024). Global soil metagenomics reveals distribution and predominance of Deltaproteobacteria in nitrogen-fixing microbiome. Microbiome 12:95. doi: 10.1186/s40168-024-01812-1,

PubMed Abstract | Crossref Full Text | Google Scholar

Meng, H., Zhou, Z. C., Wu, R. N., Wang, Y. F., and Gu, J. D. (2019). Diazotrophic microbial community and abundance in acidic subtropical natural and re-vegetated forest soils revealed by high-throughput sequencing of nifH gene. Appl. Microbiol. Biotechnol. 103, 995–1005. doi: 10.1007/s00253-018-9466-7,

PubMed Abstract | Crossref Full Text | Google Scholar

Nawaz, T., Fahad, S., Gu, L., Xu, L., and Zhou, R. (2025). Harnessing nitrogen-fixing Cyanobacteria for sustainable agriculture: opportunities, challenges, and implications for food security. Nitrogen 6:16. doi: 10.3390/nitrogen6010016

Crossref Full Text | Google Scholar

Oyserman, B. O., Flores, S. S., Griffioen, T., Pan, X., Van Der Wijk, E., Pronk, L., et al. (2022). Disentangling the genetic basis of rhizosphere microbiome assembly in tomato. Nat. Commun. 13:3228. doi: 10.1038/s41467-022-30849-9,

PubMed Abstract | Crossref Full Text | Google Scholar

Pérez-Pantoja, D., Donoso, R., Agulló, L., Córdova, M., Seeger, M., Pieper, D. H., et al. (2012). Genomic analysis of the potential for aromatic compounds biodegradation in Burkholderiales. Environ. Microbiol. 14, 1091–1117. doi: 10.1111/j.1462-2920.2011.02613.x,

PubMed Abstract | Crossref Full Text | Google Scholar

Pokhrel, N., Ferguson, L., and Zhang, L. (2025). Importance of nitrogen and zinc fertility in pecan production. Horticulturae 11:82. doi: 10.3390/horticulturae11010082

Crossref Full Text | Google Scholar

Qi, J. H., Li, X. R., Jia, R. L., Yang, H. T., Shi, Y. F., Sun, J. Y., et al. (2022). Response of biocrust nitrogenase activity to the variation of rainfall regime in the Tengger Desert, northern China. Catena 212:106093. doi: 10.1016/j.catena.2022.106093

Crossref Full Text | Google Scholar

Ren, W., Zhang, L., Maness, N., Wang, X., Tang, M., and Xu, T. (2023). Changes in the diversity of pecan (Carya illinoinensis) rhizosphere microbial community with different nitrogen fertilization, a case study in Oklahoma pecan orchard. Sci. Hortic. 321:112365. doi: 10.1016/j.scienta.2023.112365

Crossref Full Text | Google Scholar

Ren, W., Zhang, L., Tondre, B., Wang, X. W., and Xu, T. Y. (2024). The rootstock genotype shapes the diversity of pecan (Carya illinoinensis) rhizosphere microbial community. Front. Microbiol. 15:1461685. doi: 10.3389/fmicb.2024.1461685,

PubMed Abstract | Crossref Full Text | Google Scholar

Rheault, K., Lachance, D., Morency, M. J., Thiffault, É., Guittonny, M., Isabel, N., et al. (2020). Plant genotype influences physicochemical properties of substrate as well as bacterial and fungal assemblages in the rhizosphere of balsam poplar. Front. Microbiol. 11:575625. doi: 10.3389/fmicb.2020.575625,

PubMed Abstract | Crossref Full Text | Google Scholar

Sakoda, M., Maika, M., Fumitaka, S., Nobuo, S., Yong, G., Yuji, M., et al. (2019). Azoarcus sp. strain KH32C affects rice plant growth and the root-associated soil bacterial community in low nitrogen input paddy fields. Soil Sci. Plant Nutr. 65, 451–459. doi: 10.1080/00380768.2019.1666648

Crossref Full Text | Google Scholar

Sarao, S. K., Boothe, V., Das, B. K., Gonzalez-Hernandez, J. L., and Brözel, V. S. (2025). Bradyrhizobium and the soybean rhizosphere: species level bacterial population dynamics in established soybean fields, rhizosphere and nodules. Plant Soil 508, 515–530. doi: 10.1007/s11104-024-06814-4

Crossref Full Text | Google Scholar

Sasse, J., Martinoia, E., and Northen, T. (2018). Feed your friends: do plant exudates shape the root microbiome? Trends Plant Sci. 23, 25–41. doi: 10.1016/j.tplants.2017.09.003,

PubMed Abstract | Crossref Full Text | Google Scholar

Smercina Darian, N., Evans Sarah, E., Friesen Maren, L., and Tiemann Lisa, K. (2019). To fix or not to fix: controls on free-living nitrogen fixation in the rhizosphere. Appl. Environ. Microbiol. 85, e02546–e02518. doi: 10.1128/AEM.02103-19,

PubMed Abstract | Crossref Full Text | Google Scholar

Tang, K., Liang, Y. A., Yuan, B., Meng, J. Y., and Feng, F. Y. (2023). Spatial distribution and core community of diazotrophs in biological soil crusts and subsoils in temperate semi-arid and arid deserts of China. Front. Microbiol. 14:1074855. doi: 10.3389/fmicb.2023.1074855,

PubMed Abstract | Crossref Full Text | Google Scholar

Tong, Y., Wang, Z., Gong, D., Huang, C., Ma, X., Ma, X., et al. (2024). Enhancing soil fertility and elevating pecan fruit quality through combined chemical and organic fertilization practices. Horticulturae 10:25. doi: 10.3390/horticulturae10010025,

PubMed Abstract | Crossref Full Text | Google Scholar

Ulbrich, T. C., Rivas-Ubach, A., Tiemann, L. K., Friesen, M. L., and Evans, S. E. (2022). Plant root exudates and rhizosphere bacterial communities shift with neighbor context. Soil Biol. Biochem. 172:108753. doi: 10.1016/j.soilbio.2022.108753

Crossref Full Text | Google Scholar

Ulrich, N., Kirchner, V., Drucker, R., Wright, J. R., Mclimans, C. J., Hazen, T. C., et al. (2018). Response of aquatic bacterial communities to hydraulic fracturing in northwestern Pennsylvania: a five-year study. Sci. Rep. 8:5683. doi: 10.1038/s41598-018-23679-7,

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, X., Pan, J., Liu, F., Chen, H., Jiao, Z., Liu, S., et al. (2023). Seasonal asynchrony in above- and below-ground phenology in a temperate forest: carbon allocation trade-off and plant-microbe interactions. Plant Soil 492, 573–586. doi: 10.1007/s11104-023-06200-6

Crossref Full Text | Google Scholar

Wartiainen, I., Eriksson, T., Zheng, W., and Rasmussen, U. (2008). Variation in the active diazotrophic community in rice paddy—nifH PCR-DGGE analysis of rhizosphere and bulk soil. Appl. Soil Ecol. 39, 65–75. doi: 10.1016/j.apsoil.2007.11.008

Crossref Full Text | Google Scholar

Wegner, C.-E., Gorniak, L., Riedel, S., Westermann, M., and Küsel, K. (2019). Lanthanide-dependent Methylotrophs of the family Beijerinckiaceae: physiological and genomic insights. Appl. Environ. Microbiol. 86, e01830–e01819. doi: 10.1128/AEM.01830-19,

PubMed Abstract | Crossref Full Text | Google Scholar

Wells, M. L. (2021). Pecan response to nitrogen fertilizer source. HortScience 56, 368–373. doi: 10.21273/HORTSCI15620-20

Crossref Full Text | Google Scholar

Wolińska, A., Kuźniar, A., Zielenkiewicz, U., Banach, A., and Błaszczyk, M. (2018). Indicators of arable soils fatigue – bacterial families and genera: a metagenomic approach. Ecol. Indic. 93, 490–500. doi: 10.1016/j.ecolind.2018.05.033

Crossref Full Text | Google Scholar

Wu, J. D., Liu, S. J., Zhang, H. Y., Chen, S. S., Si, J. N., Liu, L., et al. (2025). Flavones enrich rhizosphere Pseudomonas to enhance nitrogen utilization and secondary root growth in Populus. Nat. Commun. 16:1461. doi: 10.1038/s41467-025-56226-w,

PubMed Abstract | Crossref Full Text | Google Scholar

Xiao, D., He, X., Wang, G., Xu, X., Hu, Y., Chen, X., et al. (2022). Network analysis reveals bacterial and fungal keystone taxa involved in straw and soil organic matter mineralization. Appl. Soil Ecol. 173:104395. doi: 10.1016/j.apsoil.2022.104395

Crossref Full Text | Google Scholar

Xie, Z., Yu, Z., Li, Y., Wang, G., Tang, C., Mathesius, U., et al. (2021). Linking rhizospheric diazotrophs to the stimulation of soybean N2 fixation in a Mollisol amended with maize straw. Plant Soil 463, 279–289. doi: 10.1007/s11104-021-04904-1

Crossref Full Text | Google Scholar

Yan, D. W., Tajima, H., Cline, L. C., Fong, R. Y., Ottaviani, J., Shapiro, H. Y., et al. (2022). Genetic modification of flavone biosynthesis in rice enhances biofilm formation of soil diazotrophic bacteria and biological nitrogen fixation. Plant Biotechnol. J. 20, 2135–2148. doi: 10.1111/pbi.13894,

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, C. Z., Yue, H. C., Ma, Z., Feng, Z. L., Feng, H. J., Zhao, L. H., et al. (2022). Influence of plant genotype and soil on the cotton rhizosphere microbiome. Front. Microbiol. 13:1021064. doi: 10.3389/fmicb.2022.1021064,

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, R.-P., He, X., Baer, M., Beirinckx, S., Tian, T., Moya, Y.a. T., et al. (2021a). Plant flavones enrich rhizosphere Oxalobacteraceae to improve maize performance under nitrogen deprivation. Nat Plants 7, 481–499. doi: 10.1038/s41477-021-00897-y

Crossref Full Text | Google Scholar

Yu, R.-P., Lambers, H., Callaway, R. M., Wright, A. J., and Li, L. (2021b). Belowground facilitation and trait matching: two or three to tango? Trends Plant Sci. 26, 1227–1235. doi: 10.1016/j.tplants.2021.07.014,

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, L. T., Li, J. J., Pokhrel, G. R., Chen, J., Zhao, Y. L., Bai, Y., et al. (2021). nifH gene sequencing reveals the effects of successive monoculture on the soil diazotrophic microbial community in Casuarina equisetifolia plantations. Front. Plant Sci. 11:578812. doi: 10.3389/fpls.2020.578812,

PubMed Abstract | Crossref Full Text | Google Scholar