Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 18 May 2023
Sec. Terrestrial Microbiology
This article is part of the Research Topic Soil Microbiome Community and Functional Succession Mechanism Driven by Different Factors in Agricultural Ecology View all 43 articles

Long-term fertilization coupled with rhizobium inoculation promotes soybean yield and alters soil bacterial community composition

Wanling WeiWanling Wei1Dawei GuanDawei Guan1Mingchao Ma,Mingchao Ma1,2Xin Jiang,Xin Jiang1,2Fenliang FanFenliang Fan1Fangang MengFangang Meng3Li Li,Li Li1,2Baisuo ZhaoBaisuo Zhao2Yubin ZhaoYubin Zhao1Fengming Cao,Fengming Cao1,2Huijun ChenHuijun Chen1Jun Li,
Jun Li1,2*
  • 1Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, China
  • 2Laboratory of Quality & Safety Risk Assessment for Microbial Products (Beijing), Ministry of Agriculture, Beijing, China
  • 3Soybean Research Institute, Jilin Academy of Agricultural Sciences, Jilin, China

Microbial diversity is an important indicator of soil fertility and plays an indispensable role in farmland ecosystem sustainability. The short-term effects of fertilization and rhizobium inoculation on soil microbial diversity and community structure have been explored extensively; however, few studies have evaluated their long-term effects. Here, we applied quantitative polymerase chain reaction (qPCR) and amplicon sequencing to characterize the effect of 10-year fertilizer and rhizobium inoculation on bacterial communities in soybean bulk and rhizosphere soils at the flowering–podding and maturity stages. Four treatments were examined: non-fertilization control (CK), phosphorus and potassium fertilization (PK), nitrogen and PK fertilization (PK + N), and PK fertilization and Bradyrhizobium japonicum 5821 (PK + R). Long-term co-application of rhizobium and PK promoted soybean nodule dry weight by 33.94% compared with PK + N, and increased soybean yield by average of 32.25%, 5.90%, and 5.00% compared with CK, PK, and PK + N, respectively. The pH of PK + R was significantly higher than that of PK and PK + N at the flowering–podding stage. The bacterial abundance at the flowering–podding stage was positively correlated with soybean yield, but not at the maturity stage. The significant different class Gemmatimonadetes, and the genera Gemmatimonas, and Ellin6067 in soil at the flowering–podding stage were negatively correlated with soybean yield. However, the bacterial community at class and genus levels at maturity had no significant effect on soybean yield. The key bacterial communities that determine soybean yield were concentrated in the flowering–podding stage, not at maturity stage. Rhizosphere effect, growth period, and treatment synergies resulted in significant differences in soil bacterial community composition. Soil organic matter (OM), total nitrogen (TN), pH, and available phosphorus (AP) were the main variables affecting bacterial community structure. Overall, long-term co-application of rhizobium and fertilizer not only increased soybean yield, but also altered soil bacterial community structure through niche reconstruction and microbial interaction. Rhizobium inoculation plays key role in reducing nitrogen fertilizer application and promoting sustainable agriculture practices.

1. Introduction

Chemical fertilizer application is a common strategy for enhancing soil fertility and crop yield in agricultural production. Effective fertilization can promote soil ecosystem health. However, irrational fertilization can lead to a series of environmental problems, including soil biodiversity and productivity reduction (Guo et al., 2010; Yang et al., 2018), environmental pollution (Probert et al., 1998; Drijber et al., 2000), and soil acidification (Ahmad et al., 2013; Xiao et al., 2021; Xie et al., 2021). At present, measures such as fallowing, organic substitution (Dai et al., 2019; Ji et al., 2020; Zhang et al., 2021), green manure planting (Zhang et al., 2017; Ma et al., 2021), and biofertilizer promotion (Wu et al., 2005; Castro et al., 2020; Zhou et al., 2021) have been implemented in agriculture to reduce the detrimental impacts of fertilization. Biofertilizers have become a research hotspot because of their environmental friendliness and effectiveness. They can activate fixed nutrients in the soil, stimulate growth and absorption by the root system, and reduce fertilization rates with the corresponding nutrient elements (Abd-Alla et al., 2014; Koskey et al., 2017). With such soil adjustments, beneficial bacteria can be increased, harmful bacteria eliminated, and crop diseases reduced (Igiehon and Babalola, 2017; Schütz et al., 2018). Soybean rhizobia are biofertilizers and facilitate nitrogen fixation and emission reduction (Soumare et al., 2020).

The Gram-negative soil bacterium rhizobium is widely distributed and can stimulate legumes to produce root nodules and symbiotic nitrogen fixation (Ditta et al., 1980; Quandt and Hynes, 1993; Han et al., 2020). Rhizobium has positive effects on crop growth and nutrient absorption, including promoting the absorption of phosphorus and iron, and enhancing the production of plant hormones that promote crop growth. They can additionally enlarge the abundance of beneficial microorganisms, minimize the prevalence of pests and diseases, and expand crop yield (Zhuang et al., 2007; Ahemad and Kibret, 2014). Trabelsi and Mhamdi (2013) reported that inoculation with two local rhizobia strains increased potato yield by 32% and decreased potato wireworm infestation by 56%. The plant hormone gibberellin produced by rhizobia enlarges the size of nodules in host legumes (Nett et al., 2022). Furthermore, inoculation with rhizobia can promote yields in crops, such as soybean (Ronner et al., 2016), milk vetch (Liu et al., 2022), faba bean (Pereira et al., 2019), maize (Marks et al., 2015), wheat (Yanni et al., 2016), and rice (Mehboob et al., 2009). In addition, the effect of rhizobia inoculum alone is reportedly weaker than that of rhizobia combined with other strains or chemical fertilizer. The nodulation rate and yield of soybean could be improved by co-inoculation of Bradyrhizobium japonicum and Azospirillum brasilense (Hungria et al., 2013; Barbosa et al., 2021; Moretti et al., 2021). Combination of Rhizobium sp. and chemical fertilizer (N15P15K15) achieved the best peanut nutrition and production parameters (Adjile et al., 2020). Rhizobium has been demonstrated to be an environmentally friendly substitute to nitrogen fertilizer based on the positive effects above (Gitonga et al., 2021; Halwani et al., 2021; Pires et al., 2021; Kumawat et al., 2022). Although rhizobium can form symbiotic relationships with legumes, and provide a considerable number of nitrogen to plant, continuous inoculation ensures the long-term availability of nitrogen in the soil. Therefore, it is an ecologically manageable choice for enhancing agricultural soil environment (Korir et al., 2017; Thilakarathna and Raizada, 2017; Zilli et al., 2021).

Several long-term studies have reported that fertilization alters soil microbial community composition. Zhou et al. (2015) observed that 34 years of fertilization decreased bacterial diversity. Long-term fertilization enhanced bacterial abundance and modified bacterial composition in rhizosphere soil (Wang et al., 2018). Inoculation with rhizobia can alter soil microbial community structure; however, most of the current studies are based on short-term time scales. One-year inoculation of rhizobia and application of appropriate nitrogen fertilizer in the field increased the bacterial richness in the rhizosphere soil (Trabelsi et al., 2011). Furthermore, rhizobium inoculation may contribute to the rotational benefits of legumes in potato cropping systems not only by providing fixed nitrogen, but also by increasing microbial diversity and structure, potentially stimulating plant growth promoting rhizobacteria and enhancing disease control (Trabelsi et al., 2012). The effects of Achillea millefolium EO (Essential oils) and three different rhizobia on soybean were studied in a greenhouse experiment. The results showed that compared with the control, the bacterial colony forming units decreased after EO application and increased after inoculation with rhizobia (Turan et al., 2019).

Soybeans [Glycine max (L.) Merr.] are native to China and their seeds are rich in protein (64%) and oil (30%; Niwińska et al., 2020; Kumawat et al., 2022). Inoculating soybean crop soil with appropriate rhizobia can supplement a high number of effective rhizobia. This can further enhance the soil fertility and increase soybean yields (Hungria et al., 2013; Igiehon et al., 2021). Microbe is a vital indicator of soil fertility and performs an integral function in farmland ecosystem sustainability. To date, the short-term effects of fertilization and rhizobium inoculation on soil bacterial diversity and community composition have been explored extensively; however, few studies have evaluated their long-term effects. In the present study, a soybean field in Northeast China with a history of 10 years of chemical fertilization and inoculation with rhizobium was examined to determine the effects of the unique fertilization strategies on bacterial community abundance and composition at the flowering–podding and maturity stages, combined with amplicon sequencing and quantitative PCR (qPCR). We hypothesized that long-term fertilization and co-inoculation with rhizobium would alter the soil bacterial community structure. Inoculation with rhizobia would improve soybean yield by increasing nodule dry weight. Our results would provide novel insights into the effects of long-term fertilization and rhizobium inoculation on soybean yield and soil bacterial community structure, and provide a theoretical basis for microbial fertilizer development and utilization.

2. Materials and methods

2.1. Experimental site and soil sampling

Since 2011, a long-term fertilization experiment was carried out in the modern agricultural industrial technology demonstration base, located at the Jilin Academy of Agricultural Sciences, Gongzhuling County, Jilin Province, China (43° 52′ 88′′ N, 124° 80′ 55′′ E, 42 m elevation) with typical chernozem. The cropping system was an annual rain-fed rotation system and the main crop was soybean, which was planted continuously throughout the year. Bradyrhizobium japonicum 5821 was isolated from the soybean root nodule of Kenjiandou 28, Jiusan Farm, Nenjiang. 5 ml of rhizobia solution (concentration was 5 × 109 CFU ml−1) was mixed with 1 kg soybean seeds. Sow 45 kg of soybean seeds per hectare. The soybean variety used was Jiyu 86.

Soil samples were collected at the flowering–podding (16 July), and maturity (29 September) stages in 2021. In this study, four treatments with three replicates were examined, and each replicate contained five random soybean plants. Namely (1) CK, non-fertilization control; (2) PK, phosphorus (75 kg P2O5 ha−1), and potassium fertilization (75 kg K2O ha−1); (3) PK + N, PK chemical fertilizers plus urea (60 kg N ha−1); and (4) PK + R, PK chemical fertilizers plus B. japonicum 5821. Topsoil (0–20 cm) and borrowed approximately 30 cm from the plants was collected as bulk soil. The rhizosphere soil of five plants was randomly collected, the excess soil was shaken off, and the soil close to the roots was gently collected with a brush and mixed into a composite sample. The soil sample was thoroughly homogenized by removing the weeds and gravel with a 2-mm sieve. Part of the soil samples were naturally air-dried and stored at 4 degrees, respectively, for physicochemical detection, and the rest were stored at −80 degrees for molecular experiments. Soybean nodules were collected from soybean roots at the flowering–podding stage, and the nodules were placed in an 80°C oven until reaching constant weight.

2.2. Soil physicochemical analysis

The ratio of soil to distilled water was 1:2.5 (weight/volume) was used to determine soil pH. Organic matter (OM) was measured according to loss on ignition of dried weight in a muffle furnace at 550°C for 6 h. CNS-2000 analyzer (LECO, St. Joseph, MI, USA) was used to estimate the TN content by burning of air-dried soil which was passed through a 0.15 mm sieve. Soil available N (AN) was measured by the diffusion plate alkaline hydrolysis method, and H2SO4 titration was used to determine its content (Zhou et al., 2019). Available phosphorus (AP) was extracted by 0.5 M NaHCO3 and the molybdenum blue colorimetric method was used for analysis (Olsen et al., 1954). Available potassium (AK) was extracted with 1 M ammonium acetate and determined by flame photometer (FP640, INASA, China).

2.3. Extraction of soil DNA and 16S rRNA gene quantitation

Soil DNA (1 g of fresh sample) was extracted using the DNeasy® PowerSoil® Kit (Qiagen, Hilden, Germany). The DNA concentration and purity were evaluated using a NanoDrop ND-1000 UV–Vis Spectrophotometer (Thermo Fisher Scientific, Rockwood, TN, USA) and 1% (w/v) agarose gel electrophoresis. The copy number of 16S rRNA gene (V4 fragment) was determined by qPCR using an ABI 7500 thermal cycler (Applied Biosystems, Waltham, MA, USA) with the primers 515FmodF (5′-GTGYCAGCMGCCGCGGTAA-3′) and 806RmodR (5′-GGACTACNVGGGTWTCTAAT-3′; Wang et al., 2018). The construction of reaction system and extraction of plasmids refer to Zhou et al. (2015). Three replicates of qPCR were performed in each group. The specificity of the amplified 16S rRNA gene was evaluated using melt curve with fluorescence measurement at temperatures ranging from 60 to 95°C. The parameter Ct (threshold period) received by ABI 7500 (version 1.0.6) was used to determine the copies of 16S rRNA gene (Zhou et al., 2019).

2.4. Amplicon sequencing of 16S rRNA gene

The purified DNA was amplified using primers 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′; Ding et al., 2020). They were sequenced on the MiSeq PE300 platform (Illumina, San Diego, CA, USA) at the Sanger Biotech Co., Ltd., Shanghai, China. In the reads of the original 16S rRNA gene, shorter sequences and those with ambiguous bases were discarded. QIIME (v. 1.9.1) was used to identify and remove chimeric and noisy sequences. Use UPARSE (v. 11) to cluster operational taxonomic units (OTUs) at 97% similarity cut-off points. A unique classification of OTUs was confirmed based on a comparison with the SILVA database (Release 138, http://www.arb-silva.de). Total potential OTUs and bacterial diversity were estimated using Mothur software (1.30.2). The original sequences were uploaded to the NCBI Sequence Read Archive under study PRJNA859097 (Chen et al., 2021; Zhang et al., 2022).

2.5. Statistical analysis

SPSS 24 (SPSS, Chicago, USA) was used to conduct one-way analysis of variance to analyze the differences in basic properties, bacterial abundance and diversity, and Duncan’s test at p < 0.05 was used to compare the significance between the soil treatments. The difference of physicochemical properties between bulk and rhizosphere soils was analyzed by T-test. The effects of fertilization, rhizosphere effect, and growth stage interaction on bacterial abundance and community diversity were analyzed by multiway analysis of variance (ANOVO) (Chen et al., 2021). Boxplot and regression analysis were carried out with the “ggplot2” package; principal coordinate analysis (PCoA) and Mantel test were proceeded by “vegan” package in R software (v 3.6.1). Redundancy analysis (RDA) was performed using CANOCO (version 5.0) to visualize the effects of soil physicochemical factors on bacterial OTU composition. The relationships among the soil properties, nodule dry weight, soybean yield, and bacterial community composition were detected by calculating Spearman correlation coefficients.

AMOS software (IBM® SPSS® Amos 26.0.0) was used to conduct structural equation modeling (SEM) to account for the direct and indirect relationships among soil main physicochemical factors, nodule dry weight, soybean yield, and bacterial community structure. The first principal component (PC1) of PCoA was used as the index of bacterial community composition, and the Shannon index was used as the index of bacterial diversity. The best fitting model was obtained based on the maximum likelihood of fit, namely, p-values, the goodness of Chi-square test (χ2) and fit index (GFI), and the approximate root mean square error (RMSEA; Lefcheck, 2015; Kwok et al., 2018).

3. Results

3.1. Variations in environmental factors among treatments

Long-term fertilization and inoculation with rhizobium changed the soil physicochemical properties, and soybean nodulation rate and yield (Table 1). Long-term nitrogen fertilization (PK + N) reduced the root nodule dry weight by 69.02% compared with PK, by 33.94% when compared with PK + R during the flowering–podding stage. The soybean yield was the highest in PK + R (3,025.33 kg ha−1). Inoculation with B. japonicum 5821 increased soybean yield by 737.73 kg ha−1 (32.25%) on average when compared with CK, by 168.66 kg ha−1 (5.90%) on average when compared with PK, and by 144.00 kg ha−1 (5.00%) on average when compared with PK + N. Long-term fertilization (PK, PK + N, PK + R) substantially decreased the soil pH, particularly in PK at flowering-podding stage and in PK+ N at maturity stage. The pH of the PK + R was notably greater than that of the PK and PK + N. Long-term application of PK+ N and PK + R significantly increased AP and AK content in soil at both stages. However, the effects of different treatments on TN and AN were inconsistent. Their contents fluctuated within a small range among different treatments, which might be related to the strong nitrogen fixation potential of the soybean roots. Additionally, at the flowering–podding stage, the pH of rhizosphere soil in all treatments was significantly higher than that of bulk soil. There was no significant difference between the bulk and rhizosphere contents of OM and AN in all treatments. Except PK + N, AK in other treatments was significantly greater in the rhizosphere than in the bulk soil. At the maturity stage, OM and AK in rhizosphere soil were significantly higher than that in bulk soil under fertilization conditions (p < 0.05). TN in bulk soil was significantly lower than that in rhizosphere soil except nitrogen fertilizer treatment.

TABLE 1
www.frontiersin.org

Table 1. Properties of the bulk soil and rhizosphere samples from soybean cultivated at the flowering–podding stage and maturity stage under different fertilization levels.

3.2. Variations in abundance and richness of bacteria

Results of the three-way ANOVA confirmed that the rhizosphere effect (p < 0.001) and treatment (p < 0.01) had remarkable influence on 16S rRNA gene abundance, whereas growth stage had no remarkable influence (Figure 1). The interactions among growth stage, rhizosphere effect, and treatment were noteworthy (p < 0.001). The abundances of soil bacteria were 14.91× 108 to 46.07 × 108 copies g−1, which in CK were obviously less than other treatments at the flowering–podding stage (Figure 1A), whereas those in PK + R were apparently greater than other treatments in bulk soil and obviously less than other treatments in rhizosphere soil at the maturity stage (Figure 1B).

FIGURE 1
www.frontiersin.org

Figure 1. Bacterial abundance (A,B) and richness (C,D) in soybean field trials at the flowering–podding stage (A,C) and maturity stage (B,D). CK: non-inoculated control in soil; PK, superphosphorus and potassium chloride; PK + N, PK chemical fertilizers plus urea; PK + R, PK chemical fertilizers plus Bradyrhizobium japonicum 5821. Different letters above bars indicate significant differences (one-way ANOVA, p < 0.05, Duncan’s multiple-range test) among different treatments at each growth stage. The overall effects of growth stage (G), rhizosphere effect (R), and treatment (T) on bacterial abundance and Shannon index were evaluated by three-way ANOVA, with the results shown at the top of the figure. *0.01 < p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001.

A complete of 7,132,055 high-quality reads had been obtained from 48 soil samples. The coverage values ranged from 96.05% to 96.52%. Three-way ANOVA confirmed that the growth stage (p < 0.001) and rhizosphere effect (p < 0.05) had obvious significance on the Shannon index, whereas the treatment had no remarkable influence. The interactions among growth stage, rhizosphere effect, and treatment were noteworthy (p < 0.05; Figure 1). At the flowering–podding stage, the Shannon index in the PK treatment was significantly higher than that in the PK + N treatment in the bulk soil, and that in CK and PK + N were greater than other treatments in the rhizosphere soil (Figure 1C). At maturity, there was no remarkable difference in Shannon index among any treatments in the bulk soil. In the rhizosphere soil, Shannon index in PK + R was 2.91%, 5.05%, and 5.05% higher than those in CK, PK, and PK + N, respectively (Figure 1D).

Linear regression analysis showed that bacterial abundance in bulk and rhizosphere soil at the flowering–podding stage was positively correlated with soybean yield (Figures 2A,C), whereas they had no significant correlation with soybean yield at the maturity stage (Figures 2E,G). The bacterial richness in bulk and rhizosphere soil at both stages was not correlated with soybean yield (Figures 2B,D,F,H). There was no significant correlation between soybean yield and root nodule dry weight (Supplementary Figure S1).

FIGURE 2
www.frontiersin.org

Figure 2. Linear regression relationships between bacterial abundance (A,C,E,G) and richness (B,D,F,H), and soybean yield, in bulk (A,B,E,F) and rhizosphere soil (C,D,G,H), at the flowering–podding stage (A–D) and at the maturity stage (E–H).

3.3. Variations in community composition of bacteria

PCoA analysis of the bacterial community structure showed that there was a significant separation among the clusters at both growth stages (p = 0.001; Figure 3). At the flowering–podding stage, 32.01% could be explained by PC1 and 29.73% by PC2 (Figure 3A). Meanwhile, at the maturity stage, 43.23% could be explained by PC1 and 7.92% by PC2 (Figure 3B). The different treatments presented remarkable separation at the flowering–podding and maturity stages. Rhizosphere effects also led to prominent separation of the bacterial community structure (Figures 3A,B). Three-way ANOVA showed that the growth stage (p < 0.001) and rhizosphere effect (p < 0.05) had significant influence on the bacterial community composition, whereas the treatment had no remarkable influence. The effects of different treatments on bacterial beta diversity were significantly different at the same stage and space scale (p = 0.001, Supplementary Table S1). Specifically, at the flowering-podding stage, the R2 value was 0.7963 in bulk soil and 0.9506 in rhizosphere soil. At the maturity stage, the R2 value was 0.7037 in bulk soil and 0.7562 in rhizosphere soil.

FIGURE 3
www.frontiersin.org

Figure 3. Principal coordinate analysis (PCoA) ordinations of bacterial community composition from the bulk soil and the rhizosphere of soybean under different fertilization levels at the flowering–podding stage (A) and the maturity stage (B). Differences in bacterial beta diversity among different fertilization treatments were determined through PERMANOVA based on the Bray–Curtis distance matrix. CK: non-inoculated control in soil; PK, superphosphorus and potassium chloride; PK + N, PK chemical fertilizers plus urea; PK + R, PK chemical fertilizers plus Bradyrhizobium japonicum 5821. The effects of growth stage (G), rhizosphere effect (R), and treatment (T) on bacterial community composition were evaluated by three-way ANOVA, with the results shown at the top of the figure. *0.01 < p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001.

The predominant bacterial classes in all samples were Alphaproteobacteria, Actinobacteria, Thermoleophilia, Gammaproteobacteria, Acidobacteria, Vicinamibacteria, and Gemmatimonadetes, accounting for 65.44–74.24% of the total sequences at the flowering–podding stage (Figure 4A) and 69.11–74.42% at the maturity stage (Figure 4B). Subsequently, we analyzed the changes in dominant bacteria communities at class level caused by fertilization and rhizobium (Supplementary Table S2). Notably, the application of the same fertilizer or rhizobium in bulk or rhizosphere soil during different growth periods resulted in different dominant bacteria in the community. For instance, when applicated with PK + R, in bulk soil at the flowering–podding stage, the relative abundances of the classes Thermoleophilia, Chloroflexia, and Bacilli were significantly (p < 0.05) increased, while the relative abundance of the classes Alphaproteobacteria, Gammaproteobacteria, Gemmatimonadetes, Bacteroidia, Saccharimonadia, and Holophagae were significantly decreased. In rhizosphere soil at the flowering–podding stage, the relative abundances of the classes Actinobacteria, Thermoleophilia, Chloroflexia, Ktedonobacteria, and Bacilli were significantly (p < 0.05) increased, while the relative abundance of the classes Gammaproteobacteria, Gemmatimonadetes, Bacteroidia, Polyangia, Saccharimonadia, and Holophagae were significantly decreased. In bulk soil at the maturity stage, the relative abundance of the class Bacteroidia and Gemmatimonadetes was significantly decreased. In rhizosphere soil at the maturity stage, the relative abundances of the classes Alphaproteobacteria, Gammaproteobacteria, and Bacteroidia were significantly (p < 0.05) increased, while the relative abundance of the classes Thermoleophilia, and Gemmatimonadetes were significantly decreased. The relative abundance of the class Gemmatimonadetes was significantly decreased in all treatments during the both growth stages.

FIGURE 4
www.frontiersin.org

Figure 4. Relative abundance of dominant bacteria at the class level (relative abundance >1%) (A,B) for each treatment at the flowering–podding stage (A) and at the maturity stage (B). CK: non-inoculated control in soil; PK, superphosphorus and potassium chloride; PK + N, PK chemical fertilizers plus urea; PK + R, PK chemical fertilizers plus Bradyrhizobium japonicum 5821.

The dominant genus was norank_f_norank_o_Gaiellales (2.51–7.71%) in bulk and rhizosphere soil at both stages, and their relative abundances varied significantly across different treatments (Figure 5). When applicated with PK + R, in bulk soil at the flowering–podding stage, the relative abundances of the genera norank_f_norank_o_Gaiellales, Gaiella, Nocardioides, and Blastococcus were significantly (p < 0.05) increased, while the relative abundances of the genera Gemmatimonas, norank_f_SC-I-84, Sphingomonas, and Ellin6067 were significantly decreased. In rhizosphere soil at the flowering–podding stage, the relative abundances of the genera Blastococcus, Nocardioides, and norank_f_norank_o_Gaiellales were significantly (p < 0.05) increased, while the relative abundance of the genera norank_f_SC-I-84, Ellin6067, norank_f_norank_o_Saccharimonadales, and Gemmatimonas were significantly decreased. In bulk soil at the maturity stage, the relative abundances of the genera norank_f_norank_o_Vicinamibacterales, Arthrobacter, and norank_f_Vicinamibacteraceae were significantly (p < 0.05) increased. Overall, co-application of fertilization and rhizobium significantly increased the relative abundance of the genera norank_f_norank_o_Gaiellales, Nocardioides, and Blastococcus, and decreased the relative abundance of Gemmatimonas, norank_f_SC-I-84, and Ellin6067 in bulk and rhizosphere soil at the flowering–podding stage. In addition, fertilization resulted in a significant reduction of norank_f_norank_o_Elsterales and norank_f_Xanthobacteraceae in bulk soil at flowering–podding stage. Application of nitrogen fertilizer resulted in significant decrease in Bacillus in bulk soil, and Mycobacterium and norank_f_norank_o_C0119 in rhizosphere soil at maturity stage (Figure 5; Supplementary Table S3). The relative abundances of Bradyrhizobium in bulk and rhizosphere soil at the flowering–podding and maturity stages are shown in Supplementary Figure S2. Specifically, the relative abundance of Bradyrhizobium was lower in PK + R compared to PK + N in bulk soil at both the flowering-podding and maturity stages. However, in rhizosphere soil at the flowering-podding stage, the relative abundance of Bradyrhizobium was higher in PK + R compared to PK + N, while the opposite trend was observed at the maturity stage.

FIGURE 5
www.frontiersin.org

Figure 5. Relative abundances of top 20 dominant genera in different treatments at flowering–podding (A,B) and maturity (C,D) stages in bulk (A,C) and rhizosphere (B,D) soil. CK: non-inoculated control in soil; PK, superphosphorus and potassium chloride; PK + N, PK chemical fertilizers plus urea; PK + R, PK chemical fertilizers plus Bradyrhizobium japonicum 5821.

Spearman correlation analysis revealed the correlation between soybean yield and major bacterial communities at class and genus levels (Figure 6). In bulk soil, at the flowering–podding stage, the class KD4–96 (Figure 6A) and the genus norank_f_norank_o_norank_c_KD4-96 (Figure 6E) were positively correlated with soybean yield. Conversely, the classes Gammaproteobacteria, Gemmatimonadetes, Bacteroidia, and Nitrospiria (Figure 6A), and the genera Gemmatimonas, Candidatus_Solibacter, Bryobacter, and Ellin6067 were negatively correlated with soybean yield (Figure 6E). The bacterial communities in rhizosphere soil at the flowering–podding stage showed that the relative abundances of the classes Chloroflexia, Bacilli, and Ktedonobacteria (Figure 6B), and the genera Bradyrhizobium, Arthrobacter, and Nocardioides (Figure 6F) were positively correlated with soybean yield. The classes Gemmatimonadetes, Acidobacteria, Polyangia, and Holophagae (Figure 6B), and the genera Gemmatimonas, norank_f_norank_o_Acidobacteriales, norank_f_SC-I-84, Ellin6067, and Candidatus_Solibacter (Figure 6F) were negatively correlated with soybean yield. In summary, the significant different class Gemmatimonadetes, and the genera Gemmatimonas and Ellin6067 in soil at the flowering–podding stage were negatively correlated with soybean yield. However, the soil bacterial community at both class and genus levels at maturity was not correlated with soybean yield (Figures 6C,D,G,H).

FIGURE 6
www.frontiersin.org

Figure 6. Spearman correlation coefficients between soybean yield and dominant bacteria at class (relative abundance >1%) (A–D) and genus level (20 most abundant) (E–H) in bulk soil (A,C,E,G) and rhizosphere soil (B,D,F,H) at the flowering–podding stage (A,B,E,F) and maturity stage (C,D,G,H). Bold font indicates classes and genera with significant differences among different treatments. *0.01 < p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001.

3.4. Factors driving bacteria variation in black soil

The RDA analysis showed the effects of soil properties on bacterial community composition at different stages (Figure 7). Overall, 29.41% of variation in the bacterial composition in these treatments at the flowering–podding stage was explained by RDA1 and 17.82% by RDA2 (Figure 7A). At the maturity stage, 68.14% of variation was explained by RDA1 and 0.51% by RDA2 (Figure 7B). Treatments at both stages in the bulk soil were separated along RDA2, whereas the treatments in rhizosphere soil were isolated from the bulk soil along RDA1. AK and TN were significantly related to bacterial communities at flowering-podding stage, and AK, OM, TN, and pH at maturity stage. This indicates the presence of a specific bacterial community composition between the bulk and rhizosphere soil treatments. Mantel test was proceeded to further investigate the influence of environmental factors on bacterial composition. There was a significant positive correlation between pH and bacterial community in the rhizosphere soil at the maturity stage (p < 0.05, Supplementary Table S4).

FIGURE 7
www.frontiersin.org

Figure 7. Redundancy analysis (RDA) profile constructed from the OTU composition of bacteria and soil properties in the bulk and rhizosphere soil samples under different fertilization levels at the flowering–podding stage (A) and the maturity stage (B). The position and length of the arrows indicate the direction and strength of the influence of soil variables on bacterial communities, respectively. The significant variables are shown by red arrows (p < 0.01). CK, non-inoculated control in soil; PK, superphosphorus and potassium chloride; PK + N, PK chemical fertilizers plus urea; PK + R, PK chemical fertilizers plus Bradyrhizobium japonicum 5821.

3.5. Integrated responses of soil properties and soybean yield on bacterial structure

Spearman’s correlation analyses between soil properties, nodule dry weight, soybean yield, and bacterial structure are shown in Supplementary Table S5. Nodule dry weight was negatively related to pH and positively related to AP and AK contents. Soybean yield was positively related to AP. Bacterial abundance was positively related to AN and AK, whereas bacterial diversity was negatively related to OM, AN, and AK. Bacterial composition was positively related to OM, AN, and AK, and negatively related to TN and bacterial abundance.

We assessed the effects of bacterial structure on main physicochemical factors and the soybean yield by using an SEM model based on Supplementary Table S5 (Figure 8A). This model fits our causal hypothesis (χ2 = 24, df = 18, p = 0.16, GFI = 0.92, RMSEA = 0.08). The effect of AP (path coefficient = 0.58) on soybean yield was significantly positive. The path coefficient of pH on bacterial composition was 0.66, which was higher than that of AP (0.42) and OM (0.21). The negative path coefficient of TN on bacterial composition was −0.42, which was higher than the path coefficient on bacterial diversity (−0.39). Bacterial composition positively affected nodule dry weight, and the path coefficient was 0.48, whereas OM and pH negatively affected nodule dry weight, the path coefficient was −0.51 and − 0.46. The final model explained the weightiness of different components, with soybean yield accounting for 34%, nodular dry weight for 33%, bacterial abundance for 6%, bacterial composition for 75%, and bacterial diversity for 29%. Details of the standardized direct and indirect effects for the SEM models are shown in Figure 8B.

FIGURE 8
www.frontiersin.org

Figure 8. Structural equation model (SEM) showing the hypothesized causal relationships among soil properties (OM, TN, pH, AP), nodule dry weight, soybean yield and bacterial abundance, diversity, and composition (A). This model resulted in a good fit to the data, with a model χ2 = 24, df = 18, p = 0.16, GFI = 0.92, RMSEA = 0.08. Red arrows indicate significant positive correlations, while blue indicates significant negative relationships (p < 0.05). R2 values represent the proportion of the variance explained for each endogenous variable. The direct and indirect effects of factors on bacterial composition were determined using SEM (B).

4. Discussion

4.1. Long-term co-application of Bradyrhizobium japonicum 5821 and fertilizer promotes soybean yield and alters soil properties

Rhizobium inoculation has been suggested to promote soybean nodules and yields (Albareda et al., 2009; Hungria et al., 2013; Zhong et al., 2019; Ulzen et al., 2020). In our study, long-term inoculation with B. japonicum 5821 and application of PK fertilizer promoted soybean nodule dry weight by 33.94% when compared with PK + N, and increased soybean yield when compared with CK, PK, and PK + N (Table 1). This was similar to the study, soybean yield increased by 10.1% (180 kg hm−1) after 8 years of prolonged field inoculation with Bradyrhizobium and NPK fertilizer application in central India (Rawat et al., 2013). However, most previous studies have been carried out on other crops based on short-term field trials, not at long-term scale. For example, 2 years of field experiments showed that inoculating rhizobia alone was not enough to promote cowpea yield; only rhizobia inoculation combined with phosphate and potassium fertilizer could promote the yield (Chiamaka, 2015). It can be seen that long-term inoculation could increase soybean yield, but short-term inoculation did not necessarily increase yield, which was related to the cumulative effects of continuous inoculation, suggesting that continuous inoculation was necessary for production. In addition, Bradyrhizobium spp. inoculation and N fertilizer application for two consecutive years increased soybean yield by 130 kg ha−1, and promoted soybean dry weight and nitrogen content, but N fertilizer application without rhizobia inoculation only increased plant dry weight, the effect on soybean yield and nitrogen content was insignificant (Ruiz Diaz et al., 2009). Consequently, rhizobium combined with appropriate fertilizer can improve the yield of leguminous crops by improving soil fertility and crop root environment, and providing sufficient nitrogen for crop growth (Quandt and Hynes, 1993; Nett et al., 2022). And B. japonicum 5821 inoculation instead of nitrogen fertilizer increased soybean yield.

Long-term nitrogen fertilization (PK + N) reduced root nodule dry weight by 69% when compared with PK, and by 33.94% when compared with PK + R, because soil biological nitrogen fixation can be inhibited by excessive nitrogen fertilizer (Liu et al., 2019; Smercina et al., 2019; Zheng et al., 2019). Excessive nitrogen application can also negatively affect soybean growth (Zhou et al., 2006). Nitrogen application weakened the symbiotic nitrogen fixation ability of soybean with different genotypes (Reinprecht et al., 2020). Nitrogen fertilizer significantly reduced rhizobium abundance treated by ambient CO2 (Liu et al., 2021). In addition, PK + R had higher soybean yield and nodulation compared with PK + N. The synergistic effects of nodulation and nitrogen fixation on yield indicated that rhizobium inoculation reduced nitrogen availability, this is in line with preceding research (Turan et al., 2019; Sanyal et al., 2020). Therefore, rhizobia inoculation is a greener and more effective agricultural management measure than nitrogen fertilizer application.

The difference of pH between bulk and rhizosphere soil is insignificant in each treatment at the flowering–podding stage, which confirmed the cumulative effect and precise effect of long-term targeted inoculation of rhizobia. Long-term fertilization resulted in soil acidification (Barak et al., 1997; Schroder et al., 2011; Yang et al., 2018). Our study also showed that long-term fertilization resulted in a significant decrease in soil pH (Table 1). The pH of PK + R was obviously greater than that of PK and PK + N at the flowering–podding stage, suggesting that inoculation with B. japonicum 5821 can prevent soil acidification. This observation is consistent with preceding research (Watkin et al., 2000; Makoi et al., 2013; Alemayehu and Dechassa, 2022). Because rhizobia was more suitable for survival in near neutral pH environment than low pH environment. Inoculating rhizobium significantly increased soil pH, further increased the availability of Ca, Na, Fe, Cu, Zn, and Mn nutrients in rhizosphere soil (Bambara and Ndakidemi, 2010), also considerably extended the absorption of nutrient elements such as P, K, Ca, and Mg in plants (Makoi et al., 2013). It is suggested that inoculation with rhizobium inhibits soil acidification by increasing nutrient availability in soil and extending absorption of nutrient in plants. The total and available nutrients of rhizosphere soil were higher than that of bulk soil at the maturity stage, which is supported by recent studies (Chen et al., 2018, 2019; Li et al., 2019). This may be related to the carbohydrate and amino acid substances in the rhizosphere secretions promoted the contents of various nutrients in the rhizosphere soil, resulted in a significantly different rhizosphere microenvironment from that in bulk soil (Ai et al., 2013). Furthermore, the pH of PK + R was observably less than that of CK at the flowering–podding stage, indicating certain limitations on the effect of rhizobia inoculation on soil pH. Thus, inoculation with rhizobia can prevent soil acidification to some extent.

4.2. Effects of long-term co-application of Bradyrhizobium japonicum 5821 and fertilizer on bacterial community composition

Rhizobia inoculation not only affected the growth of aboveground crops, but also soil microbiota. In the present study, a combination of amplicon sequencing and qPCR analyses revealed that fertilization and inoculation with rhizobia had far-reaching influence on bacterial abundance and richness (Figure 1). We observed that at the maturity stage, bacterial abundance (in bulk soil) and richness (in rhizosphere soil) of PK + R were significantly higher than those of other treatments, indicating that co-application of B. japonicum 5821 and PK chemical fertilizers increased bacterial abundance and richness. This is consistent with previous reports on single inoculation of rhizobia or single fertilization. Fall et al. (2016) observed that rhizobium inoculated with native trees of Senegalia senegal (L.) Britton in a gum Arabic production area improved soil microbial biomass and functional diversity. Phaseolus beans inoculated with two native rhizobium significantly increased bacterial richness (Trabelsi et al., 2011). Co-inoculation of alfalfa with rhizobia slightly promoted microbial diversity in rhizosphere soil (Ju et al., 2019). Wang et al. (2018) characterized rhizosphere and bulk soil bacterial communities in a 36-year fertilizer experiment and found that application of N fertilizer decreased bacterial abundance, whereas MNPK (horse plus NPK) fertilizer enhanced bacterial abundance in the maize soil. Long-term application of high phosphorous fertilizer reduced bacterial diversity in wheat rhizosphere soil (Liu et al., 2020). Liu et al. (2021) reported that inoculation with Bradyrhizobium diazoefficiens USDA 110 reduced the diversity of soybean rhizosphere microbes. These inconsistent results may be on account of the interaction between soil and plants, which leads to complex soil environments and diverse microorganisms (Wu et al., 2011; Francioli et al., 2016; Hu et al., 2018; Semenov et al., 2020). Furthermore, the bacterial abundance in PK + R was observably lower than in other treatments in rhizosphere soil at maturity (Figure 1B). First, this may be because the inoculation of rhizobia inhibits the growth of other bacteria, leading to the enrichment of rhizobia in the soybean rhizosphere soil. Secondly, rhizobium can symbiosis nitrogen fixation with soybean, and soybean growth slows down after entering the maturity stage compared with the earlier stages, resulting in a decrease in bacterial abundance (Sohn et al., 2021). Third, crop roots mainly stimulate specific rhizosphere groups, resulting in communities that become increasingly different from the bulk soil, often with lower diversity (Shi et al., 2015; Fan et al., 2017; Nuccio et al., 2020).

Co-application of B. japonicum 5821 and PK fertilizer led to significant diversification of soil bacterial community composition (Supplementary Tables S2, S3). The co-application decreased the relative abundance of the class Gemmatimonadetes in all treatments during both growth stages (Figure 4); and significantly increased the relative abundance of the genera norank_f_norank_o_Gaiellales, Nocardioides, and Blastococcus, and decreased the relative abundance of Gemmatimonas, norank_f_SC-I-84, and Ellin6067 in bulk and rhizosphere soil at the flowering–podding stage (Figure 5). Such variation in bacterial composition caused by long-term inoculation was similar to those in the short-term scale. Sun et al. (2009) carried out 1-year field inoculation experiments, found that inoculation with Sinorhizobium meliloti CCBAU01199 increased Alphaproteobacteria and Betaproteobacteria relative abundance, decreased Gammaproteobacteria, Deltaproteobacteria, Firmicutes, and Actinobacteria abundance in alfalfa rhizosphere soil. During in situ restoration of a vanadium titanomagnetite tailings dam using Pongamia pinnata for 2 years, the abundance of groups under the phylum Proteobacteria increased in rhizosphere flora, and OTUs associated with rhizobia were preferably enriched (Yu et al., 2019). Furthermore, fertilization has been shown to alter soil microbial community composition (Geisseler and Scow, 2014; Ma et al., 2018; Zhou et al., 2019; Zhu et al., 2021). Nitrogen directly affect the bacterial community composition and soil factors, NPK directly affect the fungi community composition (Cassman et al., 2016). We can imply that long-term co-application of B. japonicum 5821 and PK fertilizer induced varied changes in the bacterial community structures (Zhu et al., 2018; Chee-Sanford et al., 2019; Kalam et al., 2022).

From our results, soil bacterial abundance at the flowering–podding stage were positively correlated with soybean yield, but not at the maturity stage (Figure 2). This is related to the different ecological functions of bacterial community in plant development stage. Root exudates are the communication link between plants and soil bacterial communities. Plants at different development stages release root exudates to change the assembly of plant microbiome (Ajilogba et al., 2022). And soil microbial activity at the flowering–podding stage of soybean was more vigorous than that at the maturity stage (Xu et al., 2009; Sohn et al., 2021). Additionally, Spearman correlation results revealed that at the flowering–podding stage, the class KD4–96 and the genus norank_f_norank_o_norank_c_KD4-96 in bulk soil, the classes Chloroflexia, Bacilli, and Ktedonobacteria, and the genera Bradyrhizobium, Arthrobacter, and Nocardioides in rhizosphere soil, were positively correlated to soybean yield. The significant different class Gemmatimonadetes, the genera Gemmatimonas and Ellin6067 in bulk and rhizosphere soil were negatively correlated with soybean yield (Figures 6A,B,E,F). The results showed that these bacteria had significant effect on soybean yield. This was similar to the report by Niraula et al. (2022) who observed that the rhizosphere soil microorganisms during soybean R1 (beginning of flowering)—R2 (blooming) stages, the class Anaerolineae, family Micromonosporaceae, and genera Plantomyces, Nitrospira, and Rhizobium have important effects on soybean yield. But the soil bacterial community at both class and genus levels at maturity was not significantly correlated with soybean yield (Figures 6C,D,G,H), this further indicated that the key bacterial communities determine soybean yield were concentrated in the early stages of soybean growth.

Mantel test and Spearman correlation results showed that soil OM, TN, pH, and AP were the dominant variables affecting bacterial community composition (Supplementary Tables S4, S5). The SEM showed significant effects of the four variables and bacterial diversity on bacterial composition. Nodule dry weight was negatively affected by OM and pH, and soybean yield was positively affected by AP (Figure 8). Preceding reports agree with our results (Yan et al., 2014, 2019; Liu et al., 2021).

5. Conclusion

In the present study, we evaluated the effects of four treatments (CK, PK, PK + N, and PK + R) on the bacterial composition of soybean grown in the black soil of Northeast China at the flowering–podding and maturity stages. Long-term inoculation with B. japonicum 5821 and application of PK fertilizer increased soybean nodule dry weight and soybean yield and altered soil properties. Co-application of B. japonicum 5821 and PK increased bacterial abundance in soybean bulk soil, and reduced bacterial abundance in rhizosphere soil at the maturity stage. The classes Alphaproteobacteria, Actinobacteria, Thermoleophilia, Gammaproteobacteria, Acidobacteria, Vicinamibacteria, and Gemmatimonadetes were the dominant bacteria across all the soil samples. Co-inoculation with B. japonicum 5821 and PK fertilizer strongly altered the bacterial community composition. The key bacterial communities that determine soybean yield were concentrated in the flowering–podding stage, not at maturity stage. Soil OM, TN, pH, and AP were the dominant variables affecting bacterial composition. The results demonstrate that long-term inoculation of rhizobia has the potential to promote soybean productivity and nitrogen fixation ability, and to improve soil fertility.

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 at: NCBI—PRJNA859097.

Author contributions

WW, MM, XJ, and JL designed the study. FM provided resources. WW, DG, LL, YZ, FC, and HC performed the experiments. WW analyzed all the data, prepared the figures and tables, and wrote the first draft of the manuscript. FF, BZ, and JL edited the manuscript and checked the language. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Major Science and Technology Project of Yunnan Province (202202AE090025), the National Key Research and Development Program of China (2021YFD1700200), the Earmarked Fund for Modern Agro-industry Technology Research System (CARS-04), and the National Natural Science Foundation of China (32201320).

Acknowledgments

We sincerely thank FM of the Soybean Research Institute, Jilin Academy of Agricultural Sciences for providing us with valuable soil samples in long-term localized fertilization and inoculation with rhizobia and DG for helping in collecting soil samples.

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.2023.1161983/full#supplementary-material

References

Abd-Alla, M., El-Enany, A., Nafady, N., Khalaf, D., and Morsy, F. (2014). Synergistic interaction of Rhizobium leguminosarum bv. Viciae and arbuscular mycorrhizal fungi as a plant growth promoting biofertilizers for faba bean (Vicia faba L.) in alkaline soil. Microbiol. Res. 169, 49–58. doi: 10.1016/j.micres.2013.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Adjile, A., Didagbe, O., Ayelo, C., Boko, F., Houngnandan, P., and Mongbo, R. (2020). Inoculation of Rhizobiums sp strains to improve soil fertility: A peanut trial in Covè and Ouessè (Benin). Agric. Sci. 11, 440–455. doi: 10.4236/as.2020.114027

CrossRef Full Text | Google Scholar

Ahemad, M., and Kibret, M. (2014). Mechanisms and applications of plant growth promoting rhizobacteria: current perspective. J. King. Saud. Univ. Sci. 26, 1–20. doi: 10.1016/j.jksus.2013.05.001

CrossRef Full Text | Google Scholar

Ahmad, W., Singh, B., Dijkstra, F., and Dalal, R. (2013). Inorganic and organic carbon dynamics in a limed acid soil are mediated by plants. Soil. Boil. Biochem. 57, 549–555. doi: 10.1016/j.soilbio.2012.10.013

CrossRef Full Text | Google Scholar

Ai, C., Liang, G., Sun, J., Wang, X., He, P., and Zhou, W. (2013). Different roles of rhizosphere effect and long-term fertilization in the activity and community structure of ammonia oxidizers in a calcareous fluvo-aquic soil. Soil. Boil. Biochem. 57, 30–42. doi: 10.1016/j.soilbio.2012.08.003

CrossRef Full Text | Google Scholar

Ajilogba, C. F., Olanrewaju, O. S., and Babalola, O. O. (2022). Plant growth stage drives the temporal and spatial dynamics of the bacterial microbiome in the rhizosphere of Vigna subterranea. Front. Microbiol. 13:825377. doi: 10.3389/fmicb.2022.825377

PubMed Abstract | CrossRef Full Text | Google Scholar

Albareda, M., Rodríguez-Navarro, D., and Temprano, F. (2009). Soybean inoculation: dose, N fertilizer supplementation and rhizobia persistence in soil. Field Crop Res. 113, 352–356. doi: 10.1016/j.fcr.2009.05.013

CrossRef Full Text | Google Scholar

Alemayehu, D., and Dechassa, N. (2022). Inoculating faba bean seed with rhizobium bacteria increases the yield of the crop and saves farmers from the cost of applying phosphorus fertilizer. Int. J. Plant. Prod. 16, 261–273. doi: 10.1007/s42106-022-00187-3

CrossRef Full Text | Google Scholar

Bambara, S., and Ndakidemi, P. A. (2010). Changes in selected soil chemical properties in the rhizosphere of Phaseolus vulgaris L. supplied with Rhizobium inoculants, molybdenum and lime. Sci. Res. Essays 5, 679–684. doi: 10.1073/pnas.1000051107

CrossRef Full Text | Google Scholar

Barak, P., Jobe, B. O., Krueger, A. R., Peterson, L. A., and Laird, D. A. (1997). Effects of long-term soil acidification due to nitrogen fertilizer inputs in Wisconsin. Plant Soil 197, 61–69. doi: 10.1023/A:1004297607070

CrossRef Full Text | Google Scholar

Barbosa, J., Hungria, M., Sena, J., Poggere, G., dos Reis, A., and Corrêa, R. (2021). Meta-analysis reveals benefits of co-inoculation of soybean with Azospirillum brasilense and Bradyrhizobium spp. Brazil. Appl. Soil. Ecol. 163:103913. doi: 10.1016/j.apsoil.2021.103913

CrossRef Full Text | Google Scholar

Cassman, N., Leite, M., Pan, Y., de Hollander, M., van Veen, J., and Kuramae, E. (2016). Plant and soil fungal but not soil bacterial communities are linked in long-term fertilized grassland. Sci. Rep. 6:23680. doi: 10.1038/srep23680

PubMed Abstract | CrossRef Full Text | Google Scholar

Castro, J., Calijuri, M., Ferreira, J., Assemany, P., and Ribeiro, V. (2020). Microalgae based biofertilizer: A life cycle approach. Sci. Total Environ. 724:138138. doi: 10.1016/j.scitotenv.2020.138138

PubMed Abstract | CrossRef Full Text | Google Scholar

Chee-Sanford, J., Tian, D., and Sanford, R. (2019). Consumption of N2O and other N-cycle intermediates by Gemmatimonas aurantiaca strain T–27. Microbiology 165, 1345–1354. doi: 10.1099/mic.0.000847

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Li, K., Shang, J., Wu, Y., Chen, T., Wanyan, Y., et al. (2021). Plant growth-promoting bacteria improve maize growth through reshaping the rhizobacterial community in low-nitrogen and low-phosphorus soil. Biol. Fertil. Soils 57, 1075–1088. doi: 10.1007/s00374-021-01598-6

CrossRef Full Text | Google Scholar

Chen, S., Zhu, Y., Shao, T., Long, X., Gao, X., and Zhou, Z. (2019). Relationship between rhizosphere soil properties and disease severity in highbush blueberry (Vaccinium corymbosum). Appl. Soil Ecol. 137, 187–194. doi: 10.1016/j.apsoil.2019.02.015

CrossRef Full Text | Google Scholar

Chen, X., Ding, Z., Tang, M., and Zhu, B. (2018). Greater variations of rhizosphere effects within mycorrhizal group than between mycorrhizal group in a temperate forest. Soil. Boil. Biochem. 126, 237–246. doi: 10.1016/j.soilbio.2018.08.026

CrossRef Full Text | Google Scholar

Chiamaka, E. (2015) Growth and yield response of cowpea (Vigna unguiculata [L] Walp) to NPK fertilizer and rhizobia inoculation in the Guinea and Sudan savanna zones of Ghana. dissertation/doctoral’s thesis. Kumasi Ghana: Kwame Nkrumah University of Science and Technology.

Google Scholar

Dai, X., Zhou, W., Liu, G., Liang, G., He, P., and Liu, Z. (2019). Soil C/N and pH together as a comprehensive indicator for evaluating the effects of organic substitution management in subtropical paddy fields after application of high-quality amendments. Geoderma 337, 1116–1125. doi: 10.1016/j.geoderma.2018.11.023

CrossRef Full Text | Google Scholar

Ding, C., Adrian, L., Peng, Y., and He, J. (2020). 16S rRNA gene-based primer pair showed high specificity and quantification accuracy in detecting freshwater Brocadiales anammox bacteria. FEMS Microbiol. Ecol. 96:3. doi: 10.1093/femsec/fiaa013

PubMed Abstract | CrossRef Full Text | Google Scholar

Ditta, G., Stanfield, S., Corbin, D., and Helinski, D. (1980). Broad host range DNA cloning system for gram-negative bacteria: construction of a gene bank of Rhizobium meliloti. Proc. Natl. Acad. Sci. U. S. A. 77, 7347–7351. doi: 10.1073/pnas.77.12.7347

PubMed Abstract | CrossRef Full Text | Google Scholar

Drijber, R., Doran, J., Parkhurst, A., and Lyon, D. (2000). Changes in soil microbial community structure with tillage under long-term wheat-fallow management. Soil. Boil. Biochem. 32, 1419–1430. doi: 10.1016/s0038-0717(00)00060-2

CrossRef Full Text | Google Scholar

Fall, D., Bakhoum, N., Nourou Sall, S., Zoubeirou, A., Sylla, S., and Diouf, D. (2016). Rhizobial inoculation increases soil microbial functioning and gum arabic production of 13–year-old Senegalia senegal (L.) Britton, trees in the north part of Senegal. Front. Plant Sci. 7:1355. doi: 10.3389/fpls.2016.01355

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, K., Cardona, C., Li, Y., Shi, Y., Xiang, X., Shen, C., et al. (2017). Rhizosphere-associated bacterial network structure and spatial distribution differ significantly from bulk soil in wheat crop fields. Soil. Boil. Biochem. 113, 275–284. doi: 10.1016/j.soilbio.2017.06.020

CrossRef Full Text | Google Scholar

Francioli, D., Schulz, E., Lentendu, G., Wubet, T., Buscot, F., and Reitz, T. (2016). Mineral vs. organic amendments: microbial community structure, activity and abundance of agriculturally relevant microbes are driven by long-term fertilization strategies. Front. Microbiol. 7:1446. doi: 10.3389/fmicb.2016.01446

PubMed Abstract | CrossRef Full Text | Google Scholar

Geisseler, D., and Scow, K. (2014). Long-term effects of mineral fertilizers on soil microorganisms – A review. Soil. Boil. Biochem. 75, 54–63. doi: 10.1016/j.soilbio.2014.03.023

CrossRef Full Text | Google Scholar

Gitonga, N., Njeru, E., Cheruiyot, R., and Maingi, J. (2021). Bradyrhizobium inoculation has a greater effect on soybean growth, production and yield quality in organic than conventional farming systems. Cogent. Food. Agr. 7:1. doi: 10.1080/23311932.2021.1935529

CrossRef Full Text | Google Scholar

Guo, J., Liu, X., Zhang, Y., Shen, J., Han, W., Zhang, W., et al. (2010). Significant acidification in major Chinese croplands. Science 327, 1008–1010. doi: 10.1126/science.1182570

PubMed Abstract | CrossRef Full Text | Google Scholar

Halwani, M., Reckling, M., Egamberdieva, D., Omari, R., Bellingrath-Kimura, S., Bachinger, J., et al. (2021). Soybean nodulation response to cropping interval and inoculation in European cropping systems. Front. Plant Sci. 12:638452. doi: 10.3389/fpls.2021.638452

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, Q., Ma, Q., Chen, Y., Tian, B., Xu, L., Bai, Y., et al. (2020). Variation in rhizosphere microbial communities and its association with the symbiotic efficiency of rhizobia in soybean. ISME J. 14, 1915–1928. doi: 10.1038/s41396-020-0648-9

CrossRef Full Text | Google Scholar

Hu, Y., Xia, Y., Sun, Q., Liu, K., Chen, X., Ge, T., et al. (2018). Effects of long-term fertilization on phoD-harboring bacterial community in karst soils. Sci. Total Environ. 628-629, 53–63. doi: 10.1016/j.scitotenv.2018.01.314

PubMed Abstract | CrossRef Full Text | Google Scholar

Hungria, M., Nogueira, M., and Araujo, R. (2013). Co-inoculation of soybeans and common beans with rhizobia and azospirilla: strategies to improve sustainability. Biol. Fertil. Soils 49, 791–801. doi: 10.1007/s00374-012-0771-5

CrossRef Full Text | Google Scholar

Igiehon, N., and Babalola, O. (2017). Biofertilizers and sustainable agriculture: exploring arbuscular mycorrhizal fungi. Appl. Microbiol. Biotechnol. 101, 4871–4881. doi: 10.1007/s00253-017-8344-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Igiehon, N., Babalola, O., Cheseto, X., and Torto, B. (2021). Effects of rhizobia and arbuscular mycorrhizal fungi on yield, size distribution and fatty acid of soybean seeds grown under drought stress. Microbiol. Res. 242:126640. doi: 10.1016/j.micres.2020.126640

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, L., Ni, K., Wu, Z., Zhang, J., Yi, X., Yang, X., et al. (2020). Effect of organic substitution rates on soil quality and fungal community composition in a tea plantation with long-term fertilization. Biol. Fertil. Soils 56, 633–646. doi: 10.1007/s00374-020-01439-y

CrossRef Full Text | Google Scholar

Ju, W., Liu, L., Fang, L., Cui, Y., Duan, C., and Wu, H. (2019). Impact of co-inoculation with plant-growth-promoting rhizobacteria and rhizobium on the biochemical responses of alfalfa-soil system in copper contaminated soil. Ecotox. Environ. Safe. 167, 218–226. doi: 10.1016/j.ecoenv.2018.10.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalam, S., Basu, A., and Podile, A. (2022). Difficult-to-culture bacteria in the rhizosphere: the underexplored signature microbial groups. Pedosphere 32, 75–89. doi: 10.1016/s1002-0160(21)60062-0

CrossRef Full Text | Google Scholar

Korir, H., Mungai, N., Thuita, M., Hamba, Y., and Masso, C. (2017). Co-inoculation effect of rhizobia and plant growth promoting rhizobacteria on common bean growth in a low phosphorus soil. Front. Plant Sci. 8:141. doi: 10.3389/fpls.2017.00141

PubMed Abstract | CrossRef Full Text | Google Scholar

Koskey, G., Mburu, S., Njeru, E., Kimiti, J., Ombori, O., and Maingi, J. (2017). Potential of native rhizobia in enhancing nitrogen fixation and yields of climbing beans (Phaseolus vulgaris L.) in contrasting environments of eastern Kenya. Front. Plant Sci. 8:443. doi: 10.3389/fpls.2017.00443

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumawat, K., Singh, I., Nagpal, S., Sharma, P., Gupta, R., and Sirari, A. (2022). Co-inoculation of indigenous pseudomonas oryzihabitans and Bradyrhizobium sp. modulates the growth, symbiotic efficacy, nutrient acquisition, and grain yield of soybean. Pedosphere 32, 438–451. doi: 10.1016/s1002-0160(21)60085-1

CrossRef Full Text | Google Scholar

Kwok, O., Cheung, M., Jak, S., Ryu, E., and Wu, J. (2018). Editorial: recent advancements in structural equation modeling (SEM): from both methodological and application perspectives. Front. Psychol. 9:1936. doi: 10.3389/fpsyg.2018.01936

PubMed Abstract | CrossRef Full Text | Google Scholar

Lefcheck, J. (2015). PIECEWISESEM: piecewise structural equation modelling in R for ecology, evolution, and systematics. Methods Ecol. Evol. 7, 573–579. doi: 10.1111/2041-210x.12512

CrossRef Full Text | Google Scholar

Li, J., Zhou, L., and Lin, W. (2019). Calla lily intercropping in rubber tree plantations changes the nutrient content, microbial abundance, and enzyme activity of both rhizosphere and non-rhizosphere soil and calla lily growth. Ind. Crop. Prod. 132, 344–351. doi: 10.1016/j.indcrop.2019.02.045

CrossRef Full Text | Google Scholar

Liu, J., Han, J., Zhu, C., Cao, W., Luo, Y., Zhang, M., et al. (2021). Elevated atmospheric CO2 and nitrogen fertilization affect the abundance and community structure of rice root-associated nitrogen-fixing bacteria. Front. Microbiol. 12:628108. doi: 10.3389/fmicb.2021.628108

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Ma, Q., Hui, X., Ran, J., Ma, Q., Wang, X., et al. (2020). Long-term high-P fertilizer input decreased the total bacterial diversity but not phoD-harboring bacteria in wheat rhizosphere soil with available-P deficiency. Soil. Boil. Biochem. 149:107918. doi: 10.1016/j.soilbio.2020.107918

CrossRef Full Text | Google Scholar

Liu, X., Li, Q., Li, Y., Guan, G., and Chen, S. (2019). Paenibacillus strains with nitrogen fixation and multiple beneficial properties for promoting plant growth. PeerJ 7:e7445. doi: 10.7717/peerj.7445

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Guo, Z., and Shi, H. (2022). Rhizobium symbiosis leads to increased drought tolerance in chinese milk vetch (Astragalus sinicus L.). Agronomy 12:725. doi: 10.3390/agronomy12030725

CrossRef Full Text | Google Scholar

Liu, Y., Ma, B., Chen, W., Schlaeppi, K., Erb, M., Stirling, E., et al. (2021). Rhizobium symbiotic capacity shapes root-associated microbiomes in soybean. Front. Microbiol. 12:709012. doi: 10.3389/fmicb.2021.709012

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, B., Lv, X., Cai, Y., Chang, S., and Dyck, M. (2018). Liming does not counteract the influence of long-term fertilization on soil bacterial community structure and its co-occurrence pattern. Soil. Boil. Biochem. 123, 45–53. doi: 10.1016/j.soilbio.2018.05.003

CrossRef Full Text | Google Scholar

Ma, D., Yin, L., Ju, W., Li, X., Liu, X., Deng, X., et al. (2021). Meta-analysis of green manure effects on soil properties and crop yield in northern China. Field Crop Res. 266:108146. doi: 10.1016/j.fcr.2021.108146

CrossRef Full Text | Google Scholar

Makoi, J. H., Bambara, S., and Ndakidemi, P. A. (2013). Rhizobium inoculation and the supply of molybdenum and lime affect the uptake of macroelements in common bean (‘P. vulgaris L.’) plants. Aust. J. Crop. Sci. 7, 784–793. doi: 10.3316/informit.365376489416920

CrossRef Full Text | Google Scholar

Marks, B., Megías, M., Ollero, F., Nogueira, M., Araujo, R., and Hungria, M. (2015). Maize growth promotion by inoculation with Azospirillum brasilense and metabolites of Rhizobium tropici enriched on lipo-chitooligosaccharides (LCOs). AMB Express 5:71. doi: 10.1186/s13568-015-0154-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Mehboob, I., Naveed, M., and Zahir, Z. (2009). Rhizobial association with non-legumes: mechanisms and applications. Crit. Rev. Plant Sci. 28, 432–456. doi: 10.1080/07352680903187753

CrossRef Full Text | Google Scholar

Moretti, L., Crusciol, C., Bossolani, J., Calonego, J., Moreira, A., Garcia, A., et al. (2021). Beneficial microbial species and metabolites alleviate soybean oxidative damage and increase grain yield during short dry spells. Eur. J. Agron. 127:126293. doi: 10.1016/j.eja.2021.126293

CrossRef Full Text | Google Scholar

Nett, R., Bender, K., and Peters, R. (2022). Production of the plant hormone gibberellin by rhizobia increases host legume nodule size. ISME J. 16, 1809–1817. doi: 10.1038/s41396-022-01236-5

CrossRef Full Text | Google Scholar

Niraula, S., Rose, M., and Chang, W. (2022). Microbial co-occurrence network in the rhizosphere microbiome: its association with physicochemical properties and soybean yield at a regional scale. J. Microbiol. 60, 986–997. doi: 10.1007/s12275-022-2363-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Niwińska, B., Witaszek, K., Niedbała, G., and Pilarski, K. (2020). Seeds of n-GM soybean varieties cultivated in Poland and their processing products as high-protein feeds in cattle nutrition. Agriculture 10:174. doi: 10.3390/agriculture10050174

CrossRef Full Text | Google Scholar

Nuccio, E., Starr, E., Karaoz, U., Brodie, E., Zhou, J., Tringe, S., et al. (2020). Niche differentiation is spatially and temporally regulated in the rhizosphere. ISME J. 14, 999–1014. doi: 10.1038/s41396-019-0582-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Olsen, S., Watanabe, F., Cosper, H., Larson, W., and Nelson, L. (1954). Residual phosphorus availability in long-time rotations on calcareous soils. Soil Sci. 78, 141–152. doi: 10.1097/00010694-195408000-00008

CrossRef Full Text | Google Scholar

Pereira, S., Mucha, Â., Gonçalves, B., Bacelar, E., Látr, A., Ferreira, H., et al. (2019). Improvement of some growth and yield parameters of faba bean (Vicia faba) by inoculation with Rhizobium laguerreae and arbuscular mycorrhizal fungi. Crop Pasture Sci. 70, 595–605. doi: 10.1071/cp19016

CrossRef Full Text | Google Scholar

Pires, G., Eloá de Lima, M., Zanchi, C., Moretti de Freitas, C., Andrade de Souza, J., Andrea de Camargo, T., et al. (2021). Arbuscular mycorrhizal fungi in the rhizosphere of soybean in integrated crop livestock systems with intercropping in the pasture phase. Rhizosphere 17:100270. doi: 10.1016/j.rhisph.2020.100270

CrossRef Full Text | Google Scholar

Probert, M., Dimes, J., Keating, B., Dalal, R., and Strong, W. (1998). APSIM's water and nitrogen modules and simulation of the dynamics of water and nitrogen in fallow systems. Agric. Syst. 56, 1–28. doi: 10.1016/s0308-521x(97)00028-0

CrossRef Full Text | Google Scholar

Quandt, J., and Hynes, M. (1993). Versatile suicide vectors which allow direct selection for gene replacement in gram-negative bacteria. Gene 127, 15–21. doi: 10.1016/0378-1119(93)90611-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Rawat, A., Rao, D., and Sahu, R. (2013). Effect of soybean inoculation with Bradyrhizobium and wheat inoculation with Azotobacter on their productivity and N turnover in a vertisol. Arch. Agron. Soil Sci. 59, 1559–1571. doi: 10.1080/03650340.2012.740555

CrossRef Full Text | Google Scholar

Reinprecht, Y., Schram, L., Marsolais, F., Smith, T., Hill, B., and Pauls, K. (2020). Effects of nitrogen application on nitrogen fixation in common bean production. Front. Plant Sci. 11:1172. doi: 10.3389/fpls.2020.01172

PubMed Abstract | CrossRef Full Text | Google Scholar

Ronner, E., Franke, A., Vanlauwe, B., Dianda, M., Edeh, E., Ukem, B., et al. (2016). Understanding variability in soybean yield and response to P-fertilizer and rhizobium inoculants on farmers’ fields in northern Nigeria. Field Crop Res. 186, 133–145. doi: 10.1016/j.fcr.2015.10.023

CrossRef Full Text | Google Scholar

Ruiz Diaz, D., Pedersen, P., and Sawyer, J. (2009). Soybean response to inoculation and nitrogen application following long-term grass pasture. Crop Sci. 49, 1058–1062. doi: 10.2135/cropsci2008.08.0510

CrossRef Full Text | Google Scholar

Sanyal, D., Osorno, J., and Chatterjee, A. (2020). Influence of rhizobium inoculation on dry bean yield and symbiotic nitrogen fixation potential. J. Plant Nutr. 43, 798–810. doi: 10.1080/01904167.2020.1711946

CrossRef Full Text | Google Scholar

Schroder, J., Zhang, H., Girma, K., Raun, W., Penn, C., and Payton, M. (2011). Soil acidification from long-term use of nitrogen fertilizers on winter wheat. Soil Sci. Soc. Am. J. 75, 957–964. doi: 10.2136/sssaj2010.0187

CrossRef Full Text | Google Scholar

Schütz, L., Gattinger, A., Meier, M., Müller, A., Boller, T., Mäder, P., et al. (2018). Improving crop yield and nutrient use efficiency via biofertilization – a global meta-analysis. Front. Plant Sci. 8:2204. doi: 10.3389/fpls.2017.02204

PubMed Abstract | CrossRef Full Text | Google Scholar

Semenov, M., Krasnov, G., Semenov, V., and van Bruggen, A. (2020). Long-term fertilization rather than plant species shapes rhizosphere and bulk soil prokaryotic communities in agroecosystems. Appl. Soil Ecol. 154:103641. doi: 10.1016/j.apsoil.2020.103641

CrossRef Full Text | Google Scholar

Shi, S., Nuccio, E., Herman, D., Rijkers, R., Estera, K., Li, J., et al. (2015). Successional trajectories of rhizosphere bacterial communities over consecutive seasons. MBio 6, e00746–e00715. doi: 10.1128/mbio.00746-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Smercina, D., Evans, S., Friesen, M., and Tiemann, L. (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.02546-18

CrossRef Full Text | Google Scholar

Sohn, S., Ahn, J., Pandian, S., Oh, Y., Shin, E., Kang, H., et al. (2021). Dynamics of bacterial community structure in the rhizosphere and root nodule of soybean: impacts of growth stages and varieties. Int. J. Mol. Sci. 22:5577. doi: 10.3390/ijms22115577

PubMed Abstract | CrossRef Full Text | Google Scholar

Soumare, A., Diedhiou, A., Thuita, M., Hafidi, M., Ouhdouch, Y., Gopalakrishnan, S., et al. (2020). Exploiting biological nitrogen fixation: a route towards a sustainable agriculture. Plan. Theory 9:1011. doi: 10.3390/plants9081011

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Y., Zhang, N., Wang, E., Yuan, H., Yang, J., and Chen, W. (2009). Influence of intercropping and intercropping plus rhizobial inoculation on microbial activity and community composition in rhizosphere of alfalfa (Medicago sativa L.) and Siberian wild rye (Elymus sibiricus L.). FEMS Microbiol. Ecol. 70, 218–226. doi: 10.1111/j.1574-6941.2009.00752.x

CrossRef Full Text | Google Scholar

Thilakarathna, M., and Raizada, M. (2017). A meta-analysis of the effectiveness of diverse rhizobia inoculants on soybean traits under field conditions. Soil. Boil. Biochem. 105, 177–196. doi: 10.1016/j.soilbio.2016.11.022

CrossRef Full Text | Google Scholar

Trabelsi, D., and Mhamdi, R. (2013). Microbial inoculants and their impact on soil microbial communities: a review. Biomed. Res. Int. 2013:11. doi: 10.1155/2013/863240

PubMed Abstract | CrossRef Full Text | Google Scholar

Trabelsi, D., Ammar, H., Mengoni, A., and Mhamdi, R. (2012). Appraisal of the crop-rotation effect of rhizobial inoculation on potato cropping systems in relation to soil bacterial communities. Soil. Boil. Biochem. 54, 1–6. doi: 10.1016/j.soilbio.2012.05.013

CrossRef Full Text | Google Scholar

Trabelsi, D., Mengoni, A., Ammar, H., and Mhamdi, R. (2011). Effect of on-field inoculation of Phaseolus vulgaris with rhizobia on soil bacterial communities. FEMS Microbiol. Ecol. 77, 211–222. doi: 10.1111/j.1574-6941.2011.01102.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Turan, V., Schröder, P., Bilen, S., Insam, H., and Juárez, M. (2019). Co-inoculation effect of Rhizobium and Achillea millefolium L. oil extracts on growth of common bean (Phaseolus vulgaris L.) and soil microbial-chemical properties. Sci. Rep. 9:15178. doi: 10.1038/s41598-019-51587-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ulzen, J., Abaidoo, R., Ewusi-Mensah, N., Osei, O., Masso, C., and Opoku, A. (2020). Organic manure improves soybean response to rhizobia inoculant and P-fertilizer in northern Ghana. Front. Agron. 2:9. doi: 10.3389/fagro.2020.00009

CrossRef Full Text | Google Scholar

Wang, Q., Jiang, X., Guan, D., Wei, D., Zhao, B., Ma, M., et al. (2018). Long-term fertilization changes bacterial diversity and bacterial communities in the maize rhizosphere of Chinese Mollisols. Appl. Soil Ecol. 125, 88–96. doi: 10.1016/j.apsoil.2017.12.007

CrossRef Full Text | Google Scholar

Watkin, E., O’Hara, G., Howieson, J., and Glenn, A. (2000). Identification of tolerance to soil acidity in inoculant strains of Rhizobium leguminosarum bv. Trifolii. Soil. Boil. Biochem. 32, 1393–1403. doi: 10.1016/s0038-0717(00)00057-2

CrossRef Full Text | Google Scholar

Wu, M., Qin, H., Chen, Z., Wu, J., and Wei, W. (2011). Effect of long-term fertilization on bacterial composition in rice paddy soil. Biol. Fertil. Soils 47, 397–405. doi: 10.1007/s00374-010-0535-z

CrossRef Full Text | Google Scholar

Wu, S., Cao, Z., Li, Z., Cheung, K., and Wong, M. (2005). Effects of biofertilizer containing N-fixer, P and K solubilizers and AM fungi on maize growth: a greenhouse trial. Geoderma 125, 155–166. doi: 10.1016/j.geoderma.2004.07.003

CrossRef Full Text | Google Scholar

Xiao, Q., Huang, Y., Wu, L., Tian, Y., Wang, Q., Wang, B., et al. (2021). Long-term manuring increases microbial carbon use efficiency and mitigates priming effect via alleviated soil acidification and resource limitation. Biol. Fertil. Soils 57, 925–934. doi: 10.1007/s00374-021-01583-z

CrossRef Full Text | Google Scholar

Xie, S., Yang, F., Feng, H., Yu, Z., Liu, C., Wei, C., et al. (2021). Organic fertilizer reduced carbon and nitrogen in runoff and buffered soil acidification in tea plantations: evidence in nutrient contents and isotope fractionations. Sci. Total Environ. 762:143059. doi: 10.1016/j.scitotenv.2020.143059

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Y., Wang, G., Jin, J., Liu, J., Zhang, Q., and Liu, X. (2009). Bacterial communities in soybean rhizosphere in response to soil type, soybean genotype, and their growth stage. Soil. Boil. Biochem. 41, 919–925. doi: 10.1016/j.soilbio.2008.10.027

CrossRef Full Text | Google Scholar

Yan, J., Han, X., Chen, X., Lu, X., Chen, W., Wang, E., et al. (2019). Effects of long-term fertilization strategies on soil productivity and soybean rhizobial diversity in a Chinese Mollisol. Pedosphere 29, 784–793. doi: 10.1016/S1002-0160(17)60470-3

CrossRef Full Text | Google Scholar

Yan, J., Han, X., Ji, Z., Li, Y., Wang, E., Xie, Z., et al. (2014). Abundance and diversity of soybean-nodulating rhizobia in black soil are impacted by land use and crop management. Appl. Environ. Microbiol. 80, 5394–5402. doi: 10.1128/AEM.01135-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., Ni, K., Shi, Y., Yi, X., Zhang, Q., Fang, L., et al. (2018). Effects of long-term nitrogen application on soil acidification and solution chemistry of a tea plantation in China. Agric. Ecosyst. Environ. 252, 74–82. doi: 10.1016/j.agee.2017.10.004

CrossRef Full Text | Google Scholar

Yanni, Y., Dazzo, F., Squartini, A., Zanardo, M., Zidan, M., and Elsadany, A. (2016). Assessment of the natural endophytic association between Rhizobium and wheat and its ability to increase wheat production in the Nile delta. Plant Soil 407, 367–383. doi: 10.1007/s11104-016-2895-0

CrossRef Full Text | Google Scholar

Yu, X., Kang, X., Li, Y., Cui, Y., Tu, W., Shen, T., et al. (2019). Rhizobia population was favoured during in situ phytoremediation of vanadium-titanium magnetite mine tailings dam using Pongamia pinnata. Environ. Pollut. 255:113167. doi: 10.1016/j.envpol.2019.113167

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, G., Song, K., Miao, X., Huang, Q., Ma, J., Gong, H., et al. (2021). Nitrous oxide emissions, ammonia volatilization, and grain-heavy metal levels during the wheat season: effect of partial organic substitution for chemical fertilizer. Agric. Ecosyst. Environ. 311:107340. doi: 10.1016/j.agee.2021.107340

CrossRef Full Text | Google Scholar

Zhang, K., Maltais-Landry, G., James, M., Mendez, V., Wright, D., George, S., et al. (2022). Absolute microbiome profiling highlights the links among microbial stability, soil health, and crop productivity under long-term sod-based rotation. Biol. Fertil. Soils 58, 883–901. doi: 10.1007/s00374-022-01675-4

CrossRef Full Text | Google Scholar

Zhang, X., Zhang, R., Gao, J., Wang, X., Fan, F., Ma, X., et al. (2017). Thirty-one years of rice-rice-green manure rotations shape the rhizosphere microbial community and enrich beneficial bacteria. Soil. Boil. Biochem. 104, 208–217. doi: 10.1016/j.soilbio.2016.10.023

CrossRef Full Text | Google Scholar

Zheng, M., Zhou, Z., Luo, Y., Zhao, P., and Mo, J. (2019). Global pattern and controls of biological nitrogen fixation under nutrient enrichment: a meta-analysis. Glob. Chang. Biol. 25, 3018–3030. doi: 10.1111/gcb.14705

CrossRef Full Text | Google Scholar

Zhong, Y., Yang, Y., Liu, P., Xu, R., Rensing, C., Fu, X., et al. (2019). Genotype and rhizobium inoculation modulate the assembly of soybean rhizobacterial communities. Plant Cell Environ. 42, 2028–2044. doi: 10.1111/pce.13519

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, J., Guan, D., Zhou, B., Zhao, B., Ma, M., Qin, J., et al. (2015). Influence of 34–years of fertilization on bacterial communities in an intensively cultivated black soil in Northeast China. Soil. Boil. Biochem. 90, 42–51. doi: 10.1016/j.soilbio.2015.07.005

CrossRef Full Text | Google Scholar

Zhou, N., Mu, M., Yang, M., Zhou, Y., and Ma, M. (2021). The effect of microbial fertilizer on the growth, rhizospheric environment and medicinal quality of fritillaria taipaiensis. Horticulturae 7:500. doi: 10.3390/horticulturae7110500

CrossRef Full Text | Google Scholar

Zhou, X., Liang, Y., Chen, H., Shen, S., and Jing, Y. (2006). Effects of rhizobia inoculation and nitrogen fertilization on photosynthetic physiology of soybean. Photosynthetica 44, 530–535. doi: 10.1007/s11099-006-0066-x

CrossRef Full Text | Google Scholar

Zhou, Z., Wei, W., Shi, X., Liu, Y., He, X., and Wang, M. (2019). Twenty-six years of chemical fertilization decreased soil RubisCO activity and changed the ecological characteristics of soil cbbL-carrying bacteria in an entisol. Appl. Soil Ecol. 141, 1–9. doi: 10.1016/j.apsoil.2019.05.005

CrossRef Full Text | Google Scholar

Zhu, L., Zhou, H., Xie, X., Li, X., Zhang, D., Jia, L., et al. (2018). Effects of floodgates operation on nitrogen transformation in a lake based on structural equation modeling analysis. Sci. Total Environ. 631-632, 1311–1320. doi: 10.1016/j.scitotenv.2018.03.134

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Z., Zhou, J., Shahbaz, M., Tang, H., Liu, S., Zhang, W., et al. (2021). Microorganisms maintain C:N stoichiometric balance by regulating the priming effect in long-term fertilized soils. Appl. Soil Ecol. 167:104033. doi: 10.1016/j.apsoil.2021.104033

CrossRef Full Text | Google Scholar

Zhuang, X., Chen, J., Shim, H., and Bai, Z. (2007). New advances in plant growth-promoting rhizobacteria for bioremediation. Environ. Int. 33, 406–413. doi: 10.1016/j.envint.2006.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Zilli, J., Pacheco, R., Gianluppi, V., Smiderle, O., Urquiaga, S., and Hungria, M. (2021). Biological N2 fixation and yield performance of soybean inoculated with Bradyrhizobium. Nutr. Cycl. Agroecosyst. 119, 323–336. doi: 10.1007/s10705-021-10128-7

CrossRef Full Text | Google Scholar

Keywords: fertilization, rhizobium inoculation, bacterial community, rhizosphere soil, soybean

Citation: Wei W, Guan D, Ma M, Jiang X, Fan F, Meng F, Li L, Zhao B, Zhao Y, Cao F, Chen H and Li J (2023) Long-term fertilization coupled with rhizobium inoculation promotes soybean yield and alters soil bacterial community composition. Front. Microbiol. 14:1161983. doi: 10.3389/fmicb.2023.1161983

Received: 09 February 2023; Accepted: 12 April 2023;
Published: 18 May 2023.

Edited by:

Bin Huang, Tobacco Research Institute, Chinese Academy of Agricultural Sciences, China

Reviewed by:

Zengtao Zhong, Nanjing Agricultural University, China
Jingjing Wang, Chinese Academy of Sciences (CAS), China

Copyright © 2023 Wei, Guan, Ma, Jiang, Fan, Meng, Li, Zhao, Zhao, Cao, Chen 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: Jun Li, lijun01@caas.cn

Disclaimer: 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.