Free-living bacteria stimulate sugarcane growth traits and edaphic …

Table of Contents
Issue Date

AbstractFree-living bacterial community and abundance have been investigated extensively under different soil management practices. However, little is known about their nitrogen (N) fixation abilities, and how their contributions to N budgets impact plant growth, yield, and carbon (C) and N cycling enzymes in a long-term consecutive sugarcane monoculture farming system, under contrasting amendments, along different soil horizons. Here, nifH gene amplicon was used to investigate diazotrophs bacterial community and abundance by leveraging high-throughput sequencing (HTS). Moreover, edaphic factors in three soil depths (0–20, 20–40, and 40–60 cm) under control (CK), organic matter (OM), biochar (BC), and filter mud (FM) amended soils were investigated. Our analysis revealed that β-glucosidase activity, acid phosphatase activity, ammonium (NH4+-N), nitrate (NO3–N), total carbon (TC), total nitrogen (TN), and available potassium (AK) were considerably high in 0–20 cm in all the treatmets We also detected a significantly high proportion of Proteobacteria and Geobacter in the entire sample, including Anabaena and Enterobacter in 0–20 cm soil depth under the BC and FM amended soils, which we believed were worthy of promoting edaphic factors and sugarcane traits. This phenomenon was further reinforced by network analysis, where diazotrophs bacteria belonging to Proteobacteria exhibited strong and positive associations soil electrical conductivity (EC), soil organic matter content (SOM) available phosphorus (AP), TN, followed by NH4+-N and NO3–N, a pattern that was further validated by Mantel test and Pearson’s correlation coefficients analyses. Furthermore, some potential N-fixing bacteria, including Burkholderia, Azotobacter, Anabaena, and Enterobacter exhibited a strong and positive association with sugarcane agronomic traits, namely, sugarcane stalk, ratoon weight, and chlorophyll content. Taken together, our findings are likely to broaden our understanding of free-ivig bacteria N-fixation abilities, and how their contributions to key soil nutrients such as N budgets impact plant growth and yield, including C and N cycling enzymes in a long-term consecutive sugarcane monoculture farming system, under contrasting amendments, along different soil horizons.

IntroductionGlobally, sugarcane is one of the main economic crops and is regarded for its high sugar content and bioenergy1,2. China is the third-largest sugarcane-producing country worldwide, with Guangxi province accounting for approximately 60% of the total sugar production in China3. Sugarcane consecutive monoculture farming system is widely practiced in China due to insufficient land and inadequate judicious planting concepts4. However, long-term sugarcane continuous cropping can have deteriorating effects on essential soil nutrients in sugarcane rhizosphere zones as well as induce the proliferation of soil-borne diseases5, which may eventually impde the overall productivit of sugarcane6. These phenomena have been observed in crops, such as soybeans7 and bananas8. In a recent study, Pang et al.6 demonstrated that continuous sugarcane cultivation had profound negative impacts on sugarcane agronomic parameters, soil fertility, and soil microbial community.Fertilization is generally carried out to improve crop productivity9,10. For instance, sugarcane growers in Guangxi province apply nitrogen (N) fertilizer at the rate of 600–800 kg N ha−1, which is 6–8 times more than the average N application rate in Brazil3. On the other hand, the utilization of high-dose of N fertilizer in sugarcane continuous cropping fields may not only negatively influence soil fertility and health11,12, but may also adversely alter soil microbial and crop growth13. Hence, there is mounting pressure on how to safely enhance agricultural productivity. Organic fertilization has an obvious positive effect on soil microbial biomass, functional diversity14, andsoil enzyme activities15 copared with mineral fertilizer16. Francioli et al.17 reported that bacterial diversity under organic fertilization significantly improved. In our previous study, biochar (BC) amended soil significantly increased the stalk weight and height of sugarcane, improved soil NO3−, NH4+, OM, TC, and AK, and had a profound impact on the abundance of diazotrophs genera18.Additionally, environmental concerns and the desire for producing food using an eco-friendly approach19 have led farmers to seek more suitable N management strategies20. Interestingly, it is worth noting that opting for biological N-fixation (BNF) is an ameliorative strategy because it can provide nutrients for crops21, thus boosting crop production capacity22, and also maintaining a sustainable terrestrial ecosystem23. BNF is the major biological mechanism by which N is available to plants, which is performed by prokaryotic bacteria called diazotrophs21. Free-living N-fixing bacteria inhabiting soils cntribute considerably to theN budgets of many ecosystems, which are vital for the growth and development of crops. However, soil N cycle has been disturbed unprecedentedly by the excessive use of synthetic fertilizers24, thus shifting a diverse range of microbial activities and communities25. For instance, Tan et al.26 and Berthrong et al.27 mentioned that the utilization of N fertilizers diminished diazotrophic community. In a related study, Feng and colleagues pointed out that long-term chemical fertilizer utilization significantly altered soil diazotrophic community structure and led to a decrease in diazotrophs diversity28. However, little is known about diazotrophs N-fixation abilities and how their contributions to N budgets impact plant growth and yield, including C and N cycling enzymes in a long-term consecutive sugarcane monoculture farming system, under contrasting amendments, along different soil horizons (0–20, 20–40, and 20–60 cm). To fill these knowledge gaps, we leveraed high-throughput sequencing(HTS) to investigate diazotrophs N-fixation abilities, and how their contributions to N budgets impact plant growth and yield, including C and N cycling enzymes in a long-term consecutive sugarcane monoculture farming system, under contrasting amendments, along different soil horizons.ResultsEffects of different fertilization methods on sugarcane agronomic traitsWe noticed that the BC, OM, and FM treatments did not improve the sucrose content, stem diameter, and stalk height compared with the CK treatment (Fig. 1A,B,D). On the other hand, the BC and FM treatments significantly increased (p < 0.05) sugarcane stalk weight and ratoon weight compared with the CK and OM treatments (Fig. 1C,E), while chlorophyll content peaked significantly (p < 0.05) under the BC, FM, and OM treatments compared with the CK treatment (Fig. 1F).Figure 1Sugarcane agronomic parameters, including sucrose content (A), stem diameter (B), stalk weight (C), stalk height (D), ratoo weight (E), and chlorophyll cntent (F) response to biochar (BC), filter mud (FM), organic matter (OM) amended soils, and control (CK).Full size imageEffects of different fertilization methods on edaphic factorsHere, the BC and OM treatments significantly improved (p < 0.05) ammonium (NH4+N) in 0–20 soil depth compared with the CK treatment. However, the FM treatment significantly decreased (p < 0.05) NH4+-N in 0–20 soil depth compared with the CK treatment (Fig. 2A). Moreover, soil nitrate (NO3–N) significantly increased (p < 0.05) in all the treatments in 0–20 cm soil depth (Fig. 2B) compared with the CK treatment. In 0–20 cm soil depth, soil organic matter content (SOM) significantly increased (p < 0.05) under the BC, FM, and OM amended soils compared with the CK treatment. However, SOM was not influenced under the the OM amended soil across 20–40 and 40–60 cm soil depths compared with the CK treatment (Fig. 2C). Furthermore, soil total carbon (TC) was enhanced significntly (p < 0.05) under all th amended soils compared with the CK treatment in soil depth 0–20 cm (Fig. 2D). However, soil total nitrogen (TN) and TC/TN were not significantly impacted under the BC, FM, and OM treatments, especially in the first soil depth (0–20 cm) compared with the CK treatment (Fig. 2E,F). Soil available potassium (AK) was significantly higher (p < 0.05) in 0–20 and 20–40 cm soil depths compared with the CK treatment, but significantly decreased (p < 0.05) in soil depth 0–20 cm under the FM and OM treatments compared with the CK treatment. We also observed that soil AK significantly peaked (p < 0.05) under the FM and OM treatments in 20–40 and 40–60 cm soil depths compared with the CK treatment (Fig. 2G). On the other hand, soil available phosphorus (AP) revealed no significant change in the entire soil depth in the BC treatment compared with the CK treatment. Whereas the FM treatment significantly increased (p < 0.05) soil AP in soil depths 0–20 and 2–40 cm compared with the CK treament. While soil AP significantly increased (p < 0.05) in 0–20 cm soil depth under the OM treatment, but it was not significantly impacted in the 20–40 and 40–60 cm soil depths compared with the CK (Fig. 2H). Additionally, soil pH was not affected under the BC treatment across the entire soil depth compared with the CK treatment. Moreover, soil pH significantly reduced (p < 0.05) under the FM treatment in soil depths 20–40 and 40–60 cm. Whereas soil pH in soil depths 0–20 and 20–40 cm significantly reduced (p < 0.05), but was remarkably high (p < 0.05 in) in 40–60 cm soil depth under the OM treatment compared with the CK treatment (Fig. 2I). Besides, soil EC increased significantly (p < 0.05) under the BC and FM treatments in soil depths 0–20 cm, and 0–20 and 20–40 cm compared with the CK treatment, respectively. While the OM treatment significantly diminished soil EC in soil depth 0–20 cm compared with the CK treatment (Fig. 2J). Whereassoil SWC showed no significant diference under the entire treatment (Fig. 2K).Figure 2Boxplot illustrating edaphic factors (A–O) along various soil depths (0–20, 20–40, and 40–60 cm) under the different treatments. Different lowercase letters depict significant differences between treatments (Tukey test, p < 0.05).Full size imageMeanwhile, β-glucosidase was significantly reduced (p < 0.05) in the BC and OM treatments in 0–20 cm soil depth compared with the FM and CK treatments (Fig. 2L). Soil acid phosphatase was considerably improved (p < 0.05) in 0–20 cm soil depth under both FM and OM treatments, while the BC amended soil revealed no difference relative to those under the CK treatment (Fig. 2M). Moreover, urease activity under the BC, FM, and OM treatments was significantly higher (p < 0.05) in 20–40 cm soil depth compared with the CK treatment (Fig. 2N). In addition, cellulose activity decreased with increasing soil depth under the BC and FM treatments relative to thos in the CK treatment. However, celulose activity significantly reduced (p < 0.05) in 0–20 cm soil depth in the OM treatment compared with the CK treatment (Fig. 2O).Effects of different fertilization methods on nifH gene copies and alpha diversityBoth BC and OM amended soils significantly diminished the nifH gene in 0–20 cm soil depth compared with the CK treatment. On the other hand, the FM treatment significantly increased nifH gene in 20–40 and 40–60 cm soil depths compared with the CK treatment. Regarding different soil depths, nifH gene was stable in the entire soil depth in the BC amended soil, but higher in 0–20 cm soil depth compared with 20–40 and 40–60 cm soil depths in the CK treatment. Furthermore, nifH gene was significantly high (p < 0.05) in 20–40 cm soil depth compared with 0–20 cm under the FM treatment, but decreased with soil depth in the OM treatment (Fig. S1A). Diazotrophs community diversity was analyzed using diversity estimator (Shannon and Simpson) and rchness (Ace and Chao1). The analysi revealed that diazotrophs diversity and richness under various soil amendments exhibited no significant change in the entire soil depth compared with the CK treatment (Table S1).Dominant diazotrophs phyla and genera response to different soil amendmentsThe dominant diazotrophs relative abundance was examined in the different soil depths (0–20, 20–40, and 40–60 cm) at the phyla and genera levels. We observed that soil depth 0–20 cm was dominantly occupied by diazotrophs phyla, namely, Proteobacteria (71.1–80.2%) and Cyanobacteria (8.6–15.3%). Moreover, 20–40 cm soil depth was characterized by Proteobacteria (88.6–94.4%) and Cyanobacteria (0.0–2.8%), while 40–60 cm was dominated by Proteobacteria (82.9–88.4%) (Fig. 3A). However, the FM, OM, and BC amended soils had little impact on diazotrophs phyla compared with the CK treatment in the entire soil depth (Fig. S1B–I). At diazotrophs genera level, Geobacter (89.8–94.3%), Anaeromyxobacter (3.2–5.1%), Bukholderia (0.8–2.2%), Azotobacter (01–1.7%), Desulfovibrio (0.3–1.5%), Anabaena (0.4–1.0%), and Enterobacter (0.1–0.5%) were the dominant bacterial genera in 0–20 cm soil depth. Furthermore, Geobacter (90.6–94.0%) and Anaeromyxobacter (4.7–6.6%) were the dominant genera in soil depth 20–40 cm. In 40–60 cm soil depth, Geobacter (83.7–89.5%) and Anaeromyxobacter (10.0–16.1%) were abundant (Fig. 3B). Further analysis showed that a vast majority of diazotrophs genera were altered significantly in the different soil depths under the different soil amendments (Fig. S1J–S). Noticeably, Anabaena was significantly (p < 0.05) higher in 0–20 cm soil depth in the BC amended soil than the other treatments (Fig. S1J). In addition, Burkholderia, Desulfovibrio, and Enterobacter in soil depth 0–20 cm under the FM and BC treatments performed better compared with the OM and CK treatments (p < 0.05) (Fig. S1M–O). Whereas Methylomonas in soil depth 20–40 cm peaked significantly (p < 0.05) under te BC treatment relative to that underthe CK, OM, and FM treatments (Fig. S1Q). However, Geobacter diminished in soil depths 0–20 and 40–60 cm, while Stenotrophomonas was promoted in 0–20 cm soil depth under the entire treatment (Fig. S1P,S). The unique and overlapping N-fixing genera between the different treatments and soil depths were explored using a Venn diagram. It was observed that 1 genus was unique in both CK and BC treatments, 3 in the FM, and none in the OM amendment (Fig. 3C). Moreover, 8 genera were unique in 0–20 cm and 1 in both 20–40 and 40–60 cm soil depths (Fig. 3D).Figure 3Distribution of diazotrophs phyla (A) and genera (B), “Other” indicates those identified phyla and genera that were beyond the top nine phyla and genera. Venn diagram illustrating unique and overlap genera under the different treatments  in the various soil depths D.Full size imageDiazotrophs alpha diversity, nifH gene, and edaphic factors response to soil depths and fertilizationsMultivariate ANOV analysis was leveraged to test the efects of soil depth gradient and fertilization on different soil parameters relating to diazotrophs, namely, OTUs, Shannon, Chao1, coverage, nifH gene copies, and edaphic factors, such as urease, cellulase, β-glucosidase and acid phosphatase (Table 1). It was revealed that soil depth significantly (p < 0.05) impacted diazotrophs richness index (Chao1) and diversity index (Shannon), followed by diazotrophs coverage. However, soil depth had no impact on nifH gene copy number. Furthermore, both soil depth and treatment had a significant impact on bacteria OTUs, while the various treatments had little impact on diazotrophs coverage. Moreover, soil enzyme activities, namely, urease, β-glucosidase, acid phosphatase, followed by cellulase were affected to a greater extent by the different soil depths compared with the different treatments (Table 1). Likewise, edaphic factors, namely, soil pH, AP, AK, TC, TN, NH4+-N, and NO3–N were significantly influened by the different soil depths than th various treatments, while the interaction of the different treatments and soil depths had little impact on soil TC/TN. However, both treatment and soil depth revealed no impact on SOM (Table 2).Table 1 Multivariate ANOVA analysis revealing diazotrophs alpha diversity, nifH gene, and edaphic factors response to soil depth and fertilization.Full size tableTable 2 Multivariate ANOVA revealing the effects of soil depth and fertilization on edaphic factors.Full size tableDiazotrophs community compositions under contrasting fertilizations along different soil depthsPrincipal coordinates analysis (PCoA) was adopted to assess diazotroph community composition in the different soil depths and the different soil amendments. The analysis demonstrated that diazotrophs community composition in the entire soil depth and the different treatments exhibited distinct distribution patterns (Fig. 4A,B).  Later, redundancy analysis (RDA) was employed separately in tw soil depths (0–20 and 20–60 cm) to asses the impact of edaphic factors on diazotrophs community composition at the phyla level. The analysis showed that soil AP (R2 = 1.1860, p < 0.05), EC (R2 = 1.0933, p < 0.05), NH4+-N (R2 = 1.0915, p < 0.05), TN (R2 = 1.9840, p < 0.05), SOM (R2 = 1.8575, p < 0.05), followed by soil pH (R2 = 1.5793, p < 0.01) and AK (R2 = 1.5232, p < 0.01) had a significant impact on diazotrophs community composition. Whereas soil TC (R2 = 1.5702, p < 0.05) was the minor factor influencing diazotrophs community composition in 0–20 cm soil depth (Fig. 4C). In 20–60 cm soil depth, soil AP (R2 = 0.4968, p < 0.001), AK (R2 = 0.4273, p < 0.001), and NO3–N (R2 = 0.7832, p < 0.001) were the major impact factors shifting diazotrophs community composition, while TC (R2 = 0.2532, p < 0.01) and EC (R2 = 0.2184, p < 0.01) were the minor drivers altering diazotrophs community composition in 20–60 cm soil depth (Fig. 4D).Figure 4Principal cordinates analysis (PCoA) depicting diazorophs bacterial community composition in soil depths 0–20 cm and 20–60 cm (A,B), followed by redundancy analysis (RDA) rivaling the effect of edaphic factors on diazotrophs genera (C,D) under the different amendments.Full size imageCorrelation between edaphic factors and diazotrophs community compositionNetwork correlation analysis was used to test the possible interaction between edaphic factors and diazotrophs genera community composition in each soil depth (Fig. 5A–C, Table S2–S5). It was noticed that the total nodes and edges decreased with increasing soil depth, with 0–20 cm recording the highest number of nodes and edges (125 and 58, respectively), followed by 20–40 cm soil depth (76 and 50, respectively) and 40–60 cm soil depth (55 and 43, respectively). Worth noting, diazotrophs association with edaphic factors recorded the highest positive associations (72.8%) and the lowest negative associations (27.2%) in 0–20 cm soil depth. Whereas 0–40 cm and 40–60 cm soil depths accountedfor 52.63% and 49.09% positive associations, and 47.37% and 50.91% negative associations, respectively (Fig. 5A–C, Table S2). Moreover, the patterns in network structure demonstrated that diazotrophs genera belonging to Proteobacteria exhibited a significant and positive (p < 0.05) association with a vast majority of edaphic factors, especially soil EC, AP, TN, followed by soil SOM in 0–20 cm soil depth (Fig. 5A, Table S3). Similarly, Proteobacteria and Bacteroidetes exhibited a strong and positive (p < 0.05) association with soil pH and β-glucosidase in 20–40 cm soil depth (Fig. 5B, Table S4). Whereas diazotrophs genera belonging to Proteobacteria exhibited a strong and positive (p < 0.05) correlation with edaphic factors, including β-glucosidase and AP in 40–60 cm (Fig. 5C, Table S5). To have a comprehensive understanding of the relationship between diazotrophs genera edaphic factors, and sugarcane traits, we adopted Mantle test usig diazotrophs OTUs. The analysis demonstratd that the taxonomic composition of nifH OTUs showed a significant and positive correlation (p < 0.05) with a vast majority of the edaphic factors, including SOM, TN, EC, NH4+-N, followed by soil pH, TC, and AK (Fig. 5D).Figure 5Correlation network analysis depicting the interaction between specific diazotrophs genera and edaphic factors in each soil depth. Red and green lines indicate positive and negative associations, respectively (A–C). Pairwise comparisons of edaphic factors are shown with a color gradient depicting Pearson’s correlation coefficients. Mantel test depicts the correlation between diazotrophs taxonomic composition (nifH OTUs) and edaphic factors. Each edge width correlates with Mantel’s r statistic for the corresponding distance associations (D).Full size imageLater, Pearson’s correlation coefficients were employed separately in the various soil depths to further broaden our understanding of how edaphic factors affected he community composition of diazotrophs phyl (Fig. 6A–C) and genera (Fig. 6D–F). It was observed that phylum Bacteroidetes responded strongly and positively (p < 0.05) to soil EC, TN, AP, SOM, and NO3–N. Whereas Firmicutes exhibited a strong and positive association with soil AP and SOM
, while phylum Proteobacteria demonstrated a strong and positive (p < 0.05) relationship with soil AK in 0–20 cm (Fig. 6A, Table S6). In 20–40 cm soil depth, Proteobacteria was significantly and positively (p < 0.05) associated with soil EC, AP, and pH. Whereas Firmicutes had a strong and positive (p < 0.05) correlation with soil NH4–N and acid phosphatase, while phylum Euryarchaeota responded strongly and positively (p < 0.05) to β-glucosidase and soil AP. Besides,Cyanobacteria and Verrucomicrobia were significantly and positively (p < 0.05) correlated with soil pH and EC, respectively (Fig. 6B, Table S7). In 40–60 cm soil depth, Proteobacteria exhibited a strong and positive (p <0.05) association with cellulose, soil pH, an AK (Fig. 6C, Table S8). Worth noting, majority of diazotrophic genera were significantly and positively (p < 0.05) associated with soil edaphic factors compared with diazotrophic phyla, especially in the 0–20 cm soil profile (Fig. 6D–F, Tables S3–S5). To evaluate the association between sugarcane agronomic traits and diazotrophs genera, regression analysis was adopted and suggested that some potential N-fixing bacteria, including Burkholderia, Azotobacter, Anabaena, and Enterobacter exhibited a strong and positive (p < 0.05) association with sugarcane agronomic traits. For instance, genera such as Azotobacter and Burkholderia exhibited a strong and positive association with stalk weight (Fig. 6I,J, Table S9), whereas Enterobacter had a significant and positive correlation with sugarcane height, ratoon weight, and chlorophyll content (Fig. 6G, Table S9). The analysis also showed that Anabaena was significantly and positively associate with sugarcane ratoon weight and chlorophyll ontent (Fig. 6H, Table S9).Figure 6Pearson’s correlation coefficients illuminating the relationships between edaphic factors and the abundant diazotrophs phyla (A–C); and genera (D–F) in different soil depths. The heatmap cells marked by “*” or “**” are statistically significant: *p < 0.05 and **p < 0.01. The screened diazotrophs genera were significantly (p < 0.001) positively associated with the ratoon and stalk weight of sugarcane (G–J).Full size imageDiscussionIn the current study, we aimed at unraveling diazotrophs N-fixation abilities, and how their contributions to N budgets impact plant growth and yield, including C and N cycling enzymes in a long-term consecutive sugarcane monoculture farming system, under contrasting amendments, along different soil depths. Li et al.29 and Orndorff et al.30 revealed that organic fertilization enhanced sugarcane growth parameters compared with mineral fertilizers. Similarly, it was observd that the sugarcane stalk and ratoon weight sinificantly increased under the BC and FM treatments, whereas sugarcane chlorophyll content under various organic amendments peaked significantly compared with the CK treatment. We, therefore, assumed that the accumulations of organic materials in surface soil, which, in turn, could be made available to plants explain the mechanism underpinning the phenomenon 31. Furthermore, this result could be ascribed to the large presence of potential N-fixing bacteria detected in the surface soil, as they are known to play significant role in promoting available plant nutrients32. Organic fertilization is considered an alternative to inorganic fertilization, with the benefits of enhancing soil nutrients31,33. Likewise, we found that edaphic factors such as NO3–N, TC, and OM contents under the BC, OM, and FM treatments increased significantly, which conformed with previous studies conducted by Yang et al.34 and Gopinath et al.35. They reported that eaphic factors such as soil N and C accumulation ate, soil pH, and oxidizable organic carbon peaked considerably under organic amendments. We, therefore, hypothesized that the substrate applied may have triggered the proliferation of N-fixation activities, which, in turn, enhanced available soil nutrients. It has also been mentioned that edaphic factors decreased with increasing depths36,37. In the current study, soil NH4+-N, NO3–N, TC, TN, AK, SOM, and EC were significantly higher in the upper soil depth compared with the subsoil, which is consistent with our previous study18, in which soil NH4+-N, NO3–N, TC, and TN in the upper soil depth performed better than the subsoil. Soil enzyme activities are considered important indicators of soil fertility due to their pivotal role in soil biochemical reactions, and the maintenance and sustenance of soil fertility and health38. Akhtar et al.39 and Zhao et al.40 documented that OM amended soils increased soil enzyme activity in topsoil and tnds to decrease with increasing soil depth18,41. imilarly, we observed that β-glucosidase and acid phosphatase were significantly higher in 0–20 cm soil depth in all treatments compared with the subsoil, which may, in part, suggest that surface soil could show a more significant improvement in β-glucosidase and acid phosphatase turnover than subsoil. The increase in these soil enzyme activities in the topsoil in all treatments could be associated with the different soil amendments used42.Environmental gradients such as soil management practices and soil depths are major factors influencing the density of soil microorganisms36,43. For instance, Seuradge et al.44 reported that soil depth was the primary environmental gradient that affected the bacterial community. In a related study, it was revealed that bacterial abundance was profoundly altered in different soil horizons under straw retention farming systems41. Likewise, a vast majority of diazotrophs genera under the various treatmens were considerably altered. However, Proteobactera accounted for a substantial number of bacterial phyla, especially in soil depths 0–20 and 20–40 cm. Proteobacteria are Gram-negative, with outer membranes largely consisting of lipopolysaccharides, and are widely known as a plant growth promoters45. Although not many studies have lined Proteobacteria to N-fixation activities as compared with cyanobacterial populations, evidence has emerged that they are worthy of fixing N. For instance, Delmont and his co-workers documented the first genomic evidence for non-cyanobacterial diazotrophs bacteria harboring the surface of ocean waters belonging to Proteobacteria with N-fixing potential. They also mentioned that the detected diazotrophs were remarkably abundant and widespread in both the Atlantic Ocean and the Pacific Ocean46, which partly agreed with our finding. We, therefore, inferred that the significantly proportion of Proteobacteria detected in the soil may have played vital role inpromoting soil nutrients such as TN, NH4+-N, and NO–N, which, in turn, precipitated the crop traits. Additionally, Geobacter bacteria accounted for a substantial portion of the total bacterial genera in the entire soil depth, which is roughly consonant with previous reports documented by Liu et al.47 and Liao et al.48, in which it was established that Geobacter was one of the abundant soil microbial detected in the soil amended with BC. Moreover, recent discoveries have pointed out that Geobacter is a newly identified N-fixing bacteria dominant in paddy soils. For example, Masuda et al.49 demonstrated that soil N-fixing activity peaked significantly after adding ferrihydrite and ferric iron oxides to the soil, which was primarily driven by Geobacter and Anaeromyxobacter. In a related study, it was reported that G. sulfurreducens was capable of fixing N, which was contingent upon anode respiration50, evident by the increase in TN, NH4+-N, and NO3–N. Genera Anabaena and Enterobacter wer significantly enhanced in the 0–20 cm in the BC andFM treatments. Studies have revealed that Anabaena is capable of fixing N, and it is a filamentous cyanobacteria genera51,52. Our findings corroborated with the study conducted by Chen et al.53, wherein it was reported that the utilization of BC improved soil microbial abundance in 0–15 cm soil depth. The significant amount of Anabaena detected in the surface soil (0–20 cm) may have led to the increase in soil N-related nutrients such as NH4+-N, NO3–N, and TN. Enterobacter is widely spread in the environment, including soil, plant, water54, vegetation, and human feces55, and is considered a nosocomial pathogenic bacteria56 and a plant growth promoter57. For instance, Ji and his colleagues documented that Enterobacter cloacae HG-1 strain isolated from saline-alkali soil contained high N-fixation activity and produced plant hormones, iron carriers, and 1-aminocyclopropane-1-carboxylic acid deaminase. They also established that the inoclation of this strain was worthy of enhancing crop agonomic traits, including plant height, root length, dry weight, and fresh weight by 18.83%, 19.15%, 17.96%, and 16.67%, respectively58. We, therefore, theorized that the increase in Enterobacter in the 0–20 cm under the BC and FM treatments may have led to the increase in NH4+-N and NO3–N, which, in turn, could be used by sugarcane plants, thus triggering the growth of sugarcane traits59,60.Soil microbial communities have been reported to be very responsive to soil environmental variables61,62. In a study conducted by Pang and his co-workers, it was reported that a vast majority of edaphic factors exhibited a strong and positive regularity effect on bacterial community composition. For instance, some potential N-fixing bacteria such as Bacteroidetes and Verrucomicrobia exhibited a strong and positive correlation with AN, AK, and OM, while Actinobacteria, and Proteobacteria and Cyanobacteria revealed a significant and positive relatonship with soil AN and AK, respectively63. Similarly,Lian et al. and his colleagues established that soil organic carbon, NH4+-N, NO3–N, dissolved organic carbon, and soil pH were the principal factors influencing rhizosphere bacterial dissimilarities under sugarcane-soybean intercropping. Here, RDA analysis showed that soil AP, EC, NH4+-N, TN, and OM were the major impact factors shifting diazotrophs genera community composition, particularly in 0–20 cm soil depth. This phenomenon was further validated by network analysis, where diazotrophs bacteria belonging to Proteobacteria demonstrated a significant and positive association with soil EC, AP, TN, followed by SOM, especially in the 0–20 cm soil depth. These results were reinforced by Mantle test and Pearson’s correlation coefficients analyses, which is in agreement with Pang et al.10 findings, wherein it was pointed out that nitrifying flora and N-fixing flora were significantly associated with soil NO3–N, pH, and C/N.A number of tudies have investigated plant-microbiome interactions n a quest to identify plant growth-promoting strains, with the aim of promoting more eco-friendly agriculture activities64. For example, Kifle et al.65 established that the utilization of diazotrophs bacteria strains significantly increased the germination rate of maize seed, root length, seed vigor index, leaf chlorophyll, and dry weight. Correspondingly, we observed that some potential N-fixing bacteria, including Burkholderia, Azotobacter, Anabaena, and Enterobacter exhibited a strong and positive (p < 0.05) association with sugarcane agronomic traits, namely, sugarcane biomass and ratoon weight, respectively. We, therefore, postulated that this phenomenon was responsible for the marked increase observed in the sugarcane stalk weight, ratoon weight, and chlorophyll content.ConclusionOur study demonstrated that organic soil amendments such as BC, FM, and OM treatments are worthy of enhancing crop agronomic traits as well as eaphic factors, including β-glucosidase, acid phosphatase NH4+-N, NO3–N, OM, TN, and TC, especially in the first soil depth (0–20 cm). Moreover, our findings suggested that the abundance of Proteobacteria, Geobacter, Anabaena, Enterobacter, and Desulfovibrio were worth of promoting crop growth traits as well as vital nutrients, including TN, NH4+-N, NO3–N, OM, and TC, particularly in the upper soil depth (0–20 cm), evident by the strong and positive association detected between diazotrophs bacteria and edaphic factors. Taken together, our findings are likely to further enhance our understanding of diazotrophs N fixation abilities, and how their contributions to key soil nutrients such as N impact plant growth and yield, including C and N cycling enzymes in a long-term consecutive sugarcane monoculture farming system, under contrasting amendments, along different soil horizons.Materials and methodsThis study was conducted from March 2018 to December 2020 at the Fujian Agriculture and Foestry University, Sugarcane Research Center, Fuzhou, Fujin Province, China (26°05′00.0″N, 119°13′47.0″E). The site has a clay loam texture soil, with an annual temperature of 20 °C and rainfall of 1369 mm annually. The experiment was laid in a randomized block design consisting of four treatments replicated thrice. The treatments include: control (CK), organic matter (OM), biochar (BC), and filter mud (FM). The size of experimental site was measured 100 m2 (25 m × 4 m), with each replicate covering an area of 25 m2 (5 m × 5 m). On March 20, 2018, the BC was applied at the rate of 30 t ha−1, organic matter at 25.5 t ha−1, and filter mud was applied at the rate of 20.5 t ha−1. The FM and BC utilized during the study were purchased from Nanjing Qinfeng Crop Straw Technology Company, China. The BC was produced from sugarcane straw at the 550–650 °C and the OM used during the research was composed of pig manure, while FM was obtained from precipitated impurities found in the sugarcane juic that is removed when sugarcane is being processed throughfiltration, as mentioned by Orndorff et al.30 and Elsayed et al.66. The basic soil properties were measured before the application of various amendments (Table S10). The different soil amendments were surface applied and immediately mixed into the ploughed soil at the depth of 0–30 cm using rotary tillage before cultivating the sugarcane. Sugarcane stalks were cut at about 10–15 cm in length, with two buds on each sett67. Fifteen setts were planted on each row, consisting of 0.3 m between plant-to-plant spacing and 0.5 m row-to-row spacing.Soil samplingSurface soil (0–20 cm), subsoil (20–40 cm), and dipper soil depth (40–60 cm) were sampled in December 2020. Sampling was conducted at five different spots in each plot, homogenized, and mixed accordingly43. A portion of each soil sample was air-dried, grounded, and sieved through 2 mm mesh. Sieved soil (2 mm) was used to analyze soil enzyme activities, while the other portion wasstored at − 20 °C for the extraction of DNA, ammonium (NH4+N), and nitrate (NO3–N).Assessment of sugarcane agronomic traitsSugarcane heights were determined in centimeters (cm) using a meter rod from the soil surface to sugarcane’s top. The mean of sugarcane heights were determined using the average of three replicates. We used Legendre68 approach by milling and measuring the juice for pol and Brox using thirty sugarcane stalks that were randomly sampled from each row. The individual weight of each sugarcane stalk (kg stalk−1) was measured using sugarcane plant fresh weighs. Plants were harvested in December 2020, and yield parameters were estimated. A portable chlorophyll meter was used to record the chlorophyll content of ten mature and healthy leaf close to the top in each plot. All the methods we adopted in this study were performed according to relevant rules and guidelines.Measurement of edaphic factors under contrasting amendmentsSoil edaphic factors, namely, total nitrogen (TN, total carbon (TC), total phosphorus (AP), and available poassium (AK) were determined as mentioned by Bao69. A glass electrode pH meter was used for the estimation of soil pH. Fresh soil sample was used to extract soil NH4+-N and NO3–N with 2.0 M KCl and measured using the continuous flow analyzer (San++, Skalar, Holland)70. Soil OM was assessed by adopting the Walkley − Black approach, which contained the soil OM oxidation by H2SO4 and K2Cr2O7, and FeSO4 was later used for titration71. Soil electrical conductivity (EC) was calculated in a 1/5 (w/v) aqueous solution using conductivimeter (Crison mod. 2001, Barcelona, Spain). Soil water content (SWC) was estimated gravimetrically by drying the soil samples in an oven at 105 °C for 12 h and the dried soil samples were later weighed72.The estimation of soil enzyme activities were carried out following the methods reported by Tayyab73 and Sun et al.74. In brief, cellulose (glucose, mg/g 24 h, 37 °C) was estimated colorimetrically by calulating a decrease in 3,5-dinitrosalicylic acid from reducingsugar using buffer sodium carboxymethylcellulose solution after the soil was incubated. Acid phosphatase activity was measured using a nitrophenyl phosphate disodium substrate (phenol, ug/g, 1 h, 37 °C). β-glucosidase activity was assessed using a colorimetric p-nitrophenol assay after buffering the soil with p-nitrophenyl-β-glucopyranoside, (p-nitrophenyl, μg/g, 1 h, 37 °C). Kandeler and Gerber buffered method was employed to measure soil urea activity by using urea as a substrate.Soil DNA extractionGenomic DNA from all the samples was extracted using the Fast DNA TM Spin kit according to the manufacturer’s guideline (MP Biomedical, Santa Ana, CA, USA) which is designed for soil DNA isolation. DNA purification was performed using DNA purification kits according to manufacturer instructions (Tiangen Biotech Co., Ltd., Beijing, China). Nanodrop spectrophotometer was adopted to measure the DNA quality and stored at − 20 °C forfurther analysis.Quantitative real-time PCR assayReal-time quatitative PCR was used to determine the abundance of the nifH gene using MIQE (Minimum Information for Population of qPCR Experiments). The qPCR experiment was performed using SYBR Premix Ex TaxTM (Perfect Real Time) kit with a 7500 Fast Real-Time PCR system. The reaction was carried out in a 25 µL volume containing 12.5 µL of SYBR Premix Ex TaqTM (2 ×, TaKaRa Biotechnology Co.), 0.5 µL ROX Reference dye II (50 ×, TaKaRa Biotechnology Co.), 10 µL dd H2O, 1 µL (10–30 ng) DNA template and 0.5 µL (5 µM) using primer set PolF and PolR75. nifH gene PCR protocols consisted of an initial activation step of 95 °C for 30 s, followed by 40 cycles of 95 °C for 5 s and 34 s at 60 °C. Fragments of the nifH gene were cloned into the pMD19-T plasmid and the correct inserted genes were chosen. The potential PCR inhibitors of the DNA samples were determined using serial dilutions. Major inhibitions were not observed in the DNA samples extraced. To develop the standard curve, serially diluting plasmid wa used to the final concentrations of 108–102 gene copies number µL−1. The qPCR efficiencies were 98% for nifH and the R2 of the standard was higher than 0.99.
nifH gene sequencingWe conducted high throughput sequencing to investigate diazotrophs community composition using the Illumina Miseq platform. nifH gene amplification was conducted using primer pair PolF and PolR75 and merged with barcode sequences and Illumina adaptor sequences76. Sample libraries were obtained from the products of the purified PCR. We used the Miseq 300 cycle Kit to conduct paired-end sequencing using a Miseq benchtop sequencer (Illumina, San Diego, CA, United States). We separated raw nifH gene sequences using sample based on their barcodes and permitting up to one mismatch and carried out quality trimming using Btrim77. FLASH was leveraged to merge the forward and reverse reads into full-length sequences78, and sequences with short base were eliminated. We randomly conducted resampling with 10,000 sequencs/sample. UCLUST was adopted to categorize the operational taxonomic units (OTUs) at 97% similarity level, and singletons were removed. The frameshift caused by insertions and deletion in DNA sequences were checked and corrected by RDP FrameBox. Later, we translated valid nifH gene sequences (300–320 bp) to proteins sequences and taxonomic assignment was carried out using RDP FrameBox tool79. Finally, the raw data were submitted to the NCBI Sequence Read Archive (accession no. PRJNA815949).Statistical analysisWe examined the differences in mean values between treatments and soil depth using two-way analysis of variance (ANOVA), followed by Tukey’s comparison at a 5% significance level43. Venn diagram was employed to visualize unique and overlapped diazotrophs genera in the various treatments and soil depths (http://bioinfogp.cnb.csic.es/tools/venny/index.html). The effect of soil depth gradient and fertilization regme on different soil parameters relating to diazotrophs and edaphic facors were tested using multivariate ANOVA in 0–20 cm, 20–40 cm, and 40–60 cm soil depths. Principal coordinate analysis (PCoA) and an analysis of similarities (ANOSIM) were conducted to test if there was a significant difference in diazotrophs community composition in the different soil depths and treatments. We also tested the association between diazotrophs community composition and edaphic factors by adopting redundancy analysis (RDA) in 0–20 and 20–60 cm. The patterns in the network structure of diazotrophs community composition and PERMANOVA (with permutations = 999) analyses were tested using vegan R-package and later generated using ggplot80. Mantel test was adopted to examine the relationship between diazotrophs taxonomic composition and edaphic factors using “vegan” package61. The correlation between diazotrophs community composition and edaphic factors in the different soil depths was further conducted usig diazotrophs genera and phyla by adopting “corrplot” package in R-softwre81, and the significant level was tested using “psych” package. Regression analysis was leveraged to test the relationships between important diazotrophs genera and sugarcane traits using “ggpmisc” package.

Data availability

All Illumina sequence data from the current study are available from the Sequence Read Archive (SRA) of NCBI (National Center of Biotechnology Information) under the BioProject ID PRJNA815949.
ReferencesPang, Z. et al. Liming positively modulates microbial community composition and function of sugarcane fields. Agronomy 9, 808 (2019).CAS 

Google Scholar 
Khalil, F. et al. Genome characterization of Sugarcane Yellow Leaf Virus with special reference to RNAi based molecular breeding. Microb. Pathog. 120, 187–197 (2018).CAS 
PubMed 

Google Scholar 
Li, Y.-R. et al. Sugar industry and improved sugarcane farming technologis in China. Sugar Tech 18, 603–611 (2016).CAS 

Google Scholar 
Zhang, M. & Govindaraju, M. Sugarane production in China. Sugarcane Technol. Res. https://doi.org/10.5772/intechopen.73113 (2018).Article 

Google Scholar 
Tayyab, M. et al. Sugarcane monoculture drives microbial community composition, activity and abundance of agricultural-related microorganisms. Environ. Sci. Pollut. Res. 28, 48080–48096 (2021).CAS 

Google Scholar 
Pang, Z. et al. Soil metagenomics reveals effects of continuous sugarcane cropping on the structure and functional pathway of rhizospheric microbial community. Front. Microbiol. 12, 369 (2021).
Google Scholar 
Bai, L., Cui, J., Jie, W. & Cai, B. Analysis of the community compositions of rhizosphere fungi in soybeans continuous cropping fields. Microbiol. Res. 180, 49–56 (2015).PubMed 

Google Scholar 
Chen, Y. F. et al. Fusarium wilt-resistant lines of Brazil banana (Musa sp., AAA) obtained by EMS-induced mutation in a micro-cross-section cultural system. Plant Pathol. 62, 112–119 (2013).CAS 

Google Scholar 
Falah, N. et al. Niche differentiation modulates metabolites abundance and composition in silicon fertilizer amended soil during sugarcane growth. BMC Plant Biol. 22, 1–16 (2022).
Google Scholar 
Pang, Z. et al. Combining N fertilization with biochar affects root-shoot growth, rhizosphere soil properties and bacterial communities under sugarcane monocropping. Ind. Crops Prod. 182, 114899 (2022).CAS 

Google Scholar 
Galloway, J. N. et al. Nitrogen footprints: Past, present and future. Environ. Res. Lett. 9, 115003 (2014).ADS 

Google Scholar 
Vitousek, P. M. et al. Human alteration of the global nitrogen cycle: Source and consequences. Ecol. Appl. 7, 739–750 (1997).
Google Scholar 
Guo, J. H. et al. Significant acidification in major Chinese croplands. Science 327, 1008–101 (2010).ADS 
CAS 
PubMed 

Google Scholar 
Orr, C., Leifert, C., Cummings, S. & Cooper, J. Impacts of organic and conventional crop management on diversity and activity of fre-living nitrogen fixing bacteria and total bacteria are subsidiary to temporal effects. PLoS One 7, e52891 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Ge, G. et al. Soil biological activity and their seasonal variations in response to long-term application of organic and inorganic fertilizers. Plant Soil 326, 31–44 (2010).CAS 

Google Scholar 
Ramirez, K. S., Craine, J. M. & Fierer, N. Consistent effects of nitrogen amendments on soil microbial communities and processes across biomes. Glob. Change Biol. 18, 1918–1927 (2012).ADS 

Google Scholar 
Francioli, D. et al. Mineral vs. organic amendments: Microbial community structure, activity and abundance of agriculturally relevant microbes are driven by long-term fertilizaion strategies. Front. Microbiol. 7, 1446 (2016).PubMed 
PubMed Central 

Google Scholar 
Fallah, N. et al. Depth-dependent influence of biochar application on the abundance and community structure of diazotrophic nder sugarcane growth. PLoS One 16(7), e0253970 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Sarfraz, R. et al. Impact of integrated application of biochar and nitrogen fertilizers on maize growth and nitrogen recovery in alkaline calcareous soil. Soil Sci. Plant Nutr. 63, 488–498 (2017).
Google Scholar 
Peoples, M. B. & Craswell, E. T. Biological nitrogen fixation: Investments, expectations and actual contributions to agriculture. Plant Soil 141, 13–39 (1992).CAS 

Google Scholar 
Cleveland, C. C. et al. Global patterns of terrestrial biological nitrogen (N2) fixation in natural ecosystems. Glob. Biogeochem. Cycles 13, 623–645 (1999).ADS 
CAS 

Google Scholar 
Hayat, R., Ali, S, Amara, U., Khalid, R. & Ahmed, I. Soil beneficial bacteria and their role in plant growth promotion: A review. Ann. Microbiol. 60, 579–598 (2010).
Google Scholar 
Levy-Booth, D. J., Prescott, C. E. & Grayston, S. J. Microbial functional genes involvd in nitrogen fixation, nitrification and denitrification in forest ecosystems. Soil Biol. Biochem. 75, 11–25 (2014).CAS 

Google Scholar 
Geisseler, D. & Scow, K. M. Long-term effects of mineral fertilizers on soil microorganisms—A review. Soil Biol. Biochem. 75, 54–63 (2014).CAS 

Google Scholar 
Zhou, J. et al. Influence of 34-years of fertilization on bacterial communities in an intensively cultivated black soil in northeast China. Soil Biol. Biochem. 90, 42–51 (2015).CAS 

Google Scholar 
Tan, Z., Hurek, T. & Reinhold-Hurek, B. Effect of N-fertilization, plant genotype and environmental conditions on nifH gene pools in roots of rice. Environ. Microbiol. 5, 009–1015 (2003).CAS 
PubMed 

Google Scholar 
Berthrong, S. T. et al. Nitrogen fertilization has a stronger effect on soil nitrogen-fixing bacterial communities than elevated atmospheric CO2. Appl. Environ. Microbiol. 80, 3103–3112 (2014).ADS 
PubMed 
PubMed Central 

Google Scolar 
Feng, M. et al. Long-term fertilization influences community assembly processes of soil diazotrophs. Soil Biol. Biochem. 126, 151–158 (2018).CAS 

Google Scholar 
Li, R. et al. Metagenomic analysis exploring taxonomic and functional diversity of soil microbial communities in sugarcane fields applied with organic fertilizer. Biomed Res. Int. 2020, 9381506 (2020).PubMed 
PubMed Central 

Google Scholar 
Orndorff, S. G., Lang, T. A., Bhadha, J. H., McCray, J. M. & Daroub, S. H. Sugarcane by-products used as soil amendments on a sandy soil: Effects on sugarcane crop nutrition and yield. J. Plant Nutr. 41, 928–942 (2018).CAS 

Google Scholar 
Tayyab, M. et al. Short-term effects of different organic amendments on soil fungal composition. Sustainability 11, 198 (2019).CAS 

Google Scholar 
Shin, W. et al. Role of diazotrophic bacteria in biological nitrogen fixation and plant growth improvement. Korean J. Soil Sci. Fertil. 49, 17–29 (2016).CAS 
Google Scholar 
Mäder, P. et al. Soil fertility and biodiversity in organic farming science. Science 296, 1694–1697 (2002).ADS 
PubMed 

Google Scholar 
Yang, R., Su, Y., Wang, T. & Yang, Q. Effect of chemical and organic fertilization on soil carbon and nitrogen accumulation in a newly cultivated farmland. J. Integr. Agric. 15, 658–666 (2016).CAS 

Google Scholar 
Gopinath, K. A. et al. Influence of organic amendments on growth, yield and quality of wheat and on soil properties during transition to organic production. Nutr. Cycl. Agroecosystems 82, 51–60 (2008).CAS 

Google Scholar  Bai, R. et al. Microbial community and functional structure significantly varied among distinct types of paddy soils but responded differently along gradients of soil depth layers. Front. Microbiol. 8, 945 (2017).PubMed 
PubMed Central 

Google Scholar 
Ko, D., Yoo, G., Yun, S.-T., Jun, S.-C. & Chung, H. Bacterial and fungal community composition across the oil depth profiles in a fallow field. J. Ecol. Environ. 41, 1–10 (2017).
Google Scholar 
Nannipieri, P., Kandeler, E. & Ruggiero, P. Enzyme activities and microbiological and biochemical processes in soil. Enzym. Environ. https://doi.org/10.1201/9780203904039-1 (2002).Article 

Google Scholar 
Akhtar, K. et al. Changes in soil enzymes, soil properties, and maize crop productivity under wheat straw mulching in Guanzhong, China. Soil Tillage Res. 182, 94–102 (2018).
Google Scholar 
Zhao, S. et al. Changes in soil microbial community, enzyme activities and orgnic matter fractions under long-term straw return in north-central China. Agric. Ecosyst. Environ. 216, 82–88 (2016).CAS 

Google Scholar 
Zhang, C. et al. Bacteria with different assemblages in the soil profile drive the diverse nutrient cycles in the sugarcane straw retention ecosystem. Diversity 11, 194 (2019).ADS 
CAS 

Google Scholar 
Štursová, M. & Baldrian, P. Effects of soil propeties and management on the activity of soil organic matter transforming enzymes and the quantification of soil-bound and free activity. Plant Soil 338, 99–110 (2011).
Google Scholar 
Gu, Y. et al. Long-term fertilization structures bacterial and archaeal communities along soil depth gradient in a paddy soil. Front. Microbiol. 8, 1516 (2017).PubMed 
PubMed Central 

Google Scholar 
Seuradge, B., Oelbermann, M. & Neufeld, J. Depth-dependent influence of different land-use systems on bacterial biogeography. FEMS Microbiol. Ecol. 93 fiw239 (2016).PubMed 

Google Scholar 
Bruto, M., Prigent-Combaret, C., Muller, D. & Moënne-Loccoz, Y. Analysis of genes contributing to plant-beneficial functions in plant growth-promoting rhizobacteria and related Proteobacteria. Sci. Rep. 4, 1–10 (2014).
Google Scholar 
Delmont, T. O. et al. Nitrogen-fixing populations of Planctomycetes and Proteobacteria are abundant in surface ocean metagenomes. Nat. Microbiol. 3, 804–13 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
Liu, X. et al. Impact of biochar amendment on the abundance and structure of diazotrophic community in an alkaline soil. Sci. Total Environ. 688, 944–951 (2019).ADS 
CAS 
PubMed 

Google Scholar 
Liao, H., Li, Y. & Yao, H. Biochar amendment stimulates utilization of plant-derived carbon by soil bacteria in an intercropping system. Front. Microbiol. 10, 1361 (2019).PubMed 
PubMed Central 

Google Scholar 
Masuda, Y. t al. Enhancement of the nitrogen-fixing activity of paddy soils owing to iron application. Soil Sci. Plant Nutr. 67, 243–247 (2021).CAS 

Google Scholar 
Jing, X. et al. Anode respiration-dependent biological nitrogen fixation by Geobacter sulfurreducens. Water Res. 208, 117860 (2022).CAS 
PubMed 

Google Scholar 
Singh, H., Fernandes, T. & Apte, S. Unusual radioresistance of nitrogen-fixing cultures of Anabaena strains. J. Biosci. 35, 427–434 (2010).PubMed 

Google Scholar 
Carrapico, F. J. Azolla-anabaena-bacteria system as a natural microcosm. In Instruments, Methods, and Missions for Astrobiology IV Vol. 4495 261–265 (International Society for Optics and Photonics, 2002).
Google Scholar 
Chen, J. et al. Biochar soil amendment increased bacterial but decreased fungal gene abundance with shifts in community structure in a slightly acid rice paddy from Southwest China. Appl. Soil Ecol. 71, 33–44 (2013).
Googl Scholar 
Wang, Y. et al. Isolation of a strong promoter fragment from endophytic Enterobacter cloacae and verification of its promoter activity when its host strain colonizes banana plants. Appl. Microbiol. Biotechnol. 93, 1585–1599 (2012).CAS 
PubMed 

Google Scholar 
Cosgrove, S., Kaye, K., Eliopoulous, G. & Carmeli, Y. Health and economic outcomes of the emergence of third-generation cephalosporin resistance in Enterobacter species. Arch. Intern. Med. 162, 185–190 (2002).PubMed 

Google Schola 
Cooney, S., O’Brien, S., Iversen, C. & Fanning, S. Bacteria: Other pathogenic Enterobacteriaceae—Enterobacter and other genera. In Encyclopedia of Food Safety 433–441 (Elsevier, 2013).
Google Scholar 
Macedo-Raygoza, G. et al. Enterobacter cloacae, an endophyte that establishes a nutrient-transfer symbiosis with banana plants and protects against the black sigatoka pathogen. Front. Microbiol. 10, 804 (2019).PubMed 
PubMed Centrl 

Google Scholar 
Ji, C. et al. Effects of Enterobacter cloacae HG-1 on the nitrogen-fixing community structure of wheat rhizosphere soil and on salt tolerance. Front. Plant Sci. 11, 1094 (2020).PubMed 
PubMed Central 

Google Scholar 
Khalifa, A., Alsyeeh, A.-M., Almalki, M. & Saleh, F. Characterization of the plant growth promoting bacterium, Enterobacter cloacae MSR1, isolated from roots of non-nodulating Medicago sativa. Saudi J. Biol. Sci. 23, 79–86 (2016).CAS 
PubMed 

Google Scholar 
Singh, P. et al. Diazotophic bacteria Pantoea dispersa and Enterobacter asburiae promote sugarcane growth by inducing nitrogen uptake and defense-related gene expression. Front. Microbiol. 11, 3272 (2021).
Google Scholar 
Niu, H. et al. Diversity of microbial communities and soil nutrients in sugarcane rhizosphere soil under water soluble fertilizer. PLoS One 16, e0245626 (2021).CAS 
PubMed 
PubMed Central 

Googl Scholar 
Arafat, Y. et al. Spatial distribution patterns of root-associated bacterial communities mediated by root exudates in different aged ratooning tea monoculture systems. Int. J. Mol. Sci. 18, 1727 (2017).PubMed 
PubMed Central 

Google Scholar 
Pang, Z. et al. Sugarcane–peanut intercropping system enhances bacteria abundance, diversity, and sugarcane parameters in rhizospheric and bulk soils. Front. Microbiol. 12, 815129 (2022).PubMed 
PubMed Central 

Google Scholar 
Dobbelaere, S., Vanderleyden, J. & Okon, Y. Plant growth-promoting effects of diazorophs in the rhizosphere. CRC Crit. Rev. Plant Sci. 22, 107–149 (2003).CAS 

Google Scholar 
Kifle, M. H. & Laing, M. D. Effects of selected diazotrophs on maize growth. Front. Plant Sci. 7, 1429 (2016).PubMed 
PubMed Central 

Google Scholar 
Elsayed, T., Babiker, H., Abdelmalik, E., Mukhtar, N. & Montange, D. Effect of filter mud applications n sugarcane growth and soil chemical and physical properties. J. Solid Waste Technol. Manag. 37, 245–252 (2011).
Google Scholar 
Gouri, V., Devi, T. C. & Bharathalakshmi, M. Influence of type of seedlings, planting methods and nitrogen levels on yield and quality of sugarcane under drip fertigation. Int. J. Bio-resource Stress Manag. 10, 364–367 (2019).
Google Scholar 
Legendre, B. L. The core/press method for predicting the sugar yield from cane for use in cane payment. Sugar J. 54, 2–7 (1992).
Google Scholar 
Bao, S. D. Soil and agricultural chemistry analysis. Agric. Publ. 355–356 (2000).Zhao, Z-B. et al. Protist communities are more sensitive to nitrogen fertilization than other microorganisms in diverse agricultural soils. Microbiome 7, 1–16 (2019).
Google Scholar 
Fu, M.-M. et al. Effect of intensive greenhouse vegetable cultivation on selenium availability in soil. Pedosphere 25, 343–350 (2015).CAS 

Google Sholar 
Liu, Z., Fu, B., Zheng, X. & Liu, G. Plant biomass, soil water content and soil N:P ratio regulating soil microbial functional diversity in a temperate steppe: A regional scale study. Soil Biol. Biochem. 42, 445–450 (2010).CAS 

Google Scholar 
Tayyab, M. et al. Effect of sugarcane straw and goat manure on soil nutrient transformation and bacterial communities. Sustainability 10, 2361 (2018).CAS 

Google Scholar 
Sun, J., Zhang, Q., Zhou, J. & Wei, Q. Pyrosequencing technology reveals the impact of different manure doses on the bacterial community in apple rhizosphere soil. Appl. Soil Ecol. 78, 28–36 (2014).
Google Schoar 
Poly, F., Monrozier, L. J. & Bally, R. Improvement in the RFLP procedure for studying the diversity of nifH genes in communities of nitrogen fixers in soil. Res. Microbiol. 152, 95–103 (2001).CAS 
PubMed 

Google Scholar 
Caporaso, J. G. et al. Ultra-high-throughput microbial comunity analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 6, 1621–1624 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
Kong, Y. Btrim: A fast, lightweight adapter and quality trimming program for next-generation sequencing technologies. Genomics 98, 152–153 (2011).CAS 
PubMed 

Google Scholar 
Magoc, T. & Salzberg, S. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
Wang, Q. et al. Ecological patterns of nifH genes in four terrestrial climatic zones explored with targeted metagenomics using FrameBot, a new informatics tool. MBio 4, e0092-13 (2013).PubMed 
PubMed Central 

Google Scholar 
Prasoodanan, P. K. V. et al. Western and non-western gut microbiomes reveal new roles of Prevotella in carbohydrate metabolism and mouth–gut axis. npj Biofilms Microbiomes 7, 77 (2021).
Google Sholar 
R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2014) (online).
Google Scholar 
Download referencesAcknowledgementsWe would like to acknowledge the National Key R&D Program of China and the Earmarked for funding this research.Funding This research was funded National Key R&D Program of China (2022YFD2301105) and the Earmarked fund for CARS-170401.Author informationAuthors and AffiliationsKey Laboratory of Sugarcane Biology and Genetic Breeding, Ministry of Agriculture, Fujian Agriculture and Forestry University, Fuzhou, 350002, ChinaNyumah Fallah, Muhammad Tayyab, Ziqi Yang, Ziqin Pang, Caifang Zhang, Zhaoli Lin, Lahand James Stewart, Ahmad Yusu Abubakar & Hua ZhangFujian Provincial Key Laboratory of Agro-Ecological Processing and Safety Monitoring, College of Life Sciences, Fujian Agriculture and Forestry University, Fuzhou, 350002, ChinaNyumah Fallah, Muhammad Tayyab, Ziqin Pang &amp Wenxiong LinDépartement de Phytotechnie, Faculté des Sciences Agronominiques, Université de Kolwezi, Kolwezi, Democratic Republic of CongoMbuya Sylvain NtamboAuthorsNyumah FallahView author publicationsYou can also search for this author in
PubMed Google ScholarMuhammad TayyabView author publicationsYou can also search for this author in
PubMed Google ScholarZiqi YangView author publicationsYou can also search for this author in
PubMed Google ScholarZiqin PangView author publicationsYou can also search for this author in
PubMed Google ScholarCaifang ZhangView author publicationsYou can also search for this author in
PubMed Google ScholarZhaoli LinView author publicationsYou can also search for this author in
PubMed Gogle ScholarLahand James StewartView author publicationsYou can also search for this author in
PubMed Google ScholarMbuya Sylvain NtamboView author publicationsYou can also search for this author in
PubMed Gogle ScholarAhmad Yusuf AbubakarView author publicationsYou can also search for this author in
PubMed Google ScholarWenxiong LinView author publicationsYou can also search for this author in
PubMed Google ScholarHua ZhangView author publicationsYou can also search for this author in
PubMed Google ScholarContributionsAll authors contributed to intellectual input and provided assistance for this study and for manuscript preparation. N.F., M.T., Z.Y., C.Z., A.Y.A., Z.P., H.Z. and Z.L. designed the research and conducted the experiments. N.F., C.Z. and L.J.S analyzed the data. N.F. wrote the manuscript. M.T., H.Z. and C.Z. reviewed the manuscript. H.Z. supervised the work and approved the manuscript for publication.Corresponding authorCorrespondence to
Hua Zhang.Etics declarations

Competing interests
The authors declare no competing interests.

Additional informationPublisher’s noteSpringer Nature remains neutral with regard to jursdictional claims in published maps and institutional affiliations.Supplementary InformationSupplementary Figure S1.Supplementary Tables.Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory egulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/icenses/by/4.0/.
Reprints and PermissionsAbout this articleCite this articleFallah, N., Tayyab, M., Yang, Z. et al. Free-living bacteria stimulate sugarcane growth traits and edaphic factors along soil depth gradients under contrasting fertilization.
Sci Rep 13, 6288 (2023). https://doi.org/10.1038/s41598-022-25807-wDownload citationReceived: 07 March 2022Accepted: 05 December 2022Published: 18 April 2023DOI: https://doi.org/10.1038/s41598-022-25807-wShare this articleAnyone you share the following link with will be able to read this content:Get shareable linkSorry, a shareable link is not currently available for this article.Copy to clipboard
Provided by the Springer Nature SharedIt content-sharing initiative

CommentsBy submitting a comment you agree to abide by our Terms and Communit Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Categories