Front. Microbiol. Frontiers in Microbiology Front. Microbiol. 1664-302X Frontiers Media S.A. 10.3389/fmicb.2023.1132016 Microbiology Original Research Effect of endophytic diazotroph Enterobacter roggenkampii ED5 on nitrogen-metabolism-related microecology in the sugarcane rhizosphere at different nitrogen levels Guo Dao-Jun 1 2 † Li Dong-Ping 3 † Yang Bin 1 Verma Krishan K. 2 Singh Rajesh Kumar 2 Singh Pratiksha 2 Khan Qaisar 4 Sharma Anjney 2 Qin Ying 4 Zhang Bao-Qing 2 Song Xiu-Peng 2 * Li Yang-Rui 2 * 1College of Life Sciences and Engineering, Hexi University, Zhangye, Gansu, China 2Key Laboratory of Sugarcane Biotechnology and Genetic Improvement (Guangxi), Ministry of Agriculture, Guangxi Key Laboratory of Sugarcane Genetic Improvement, Sugarcane Research Institute, Guangxi Academy of Agricultural Sciences, Nanning, Guangxi, China 3Microbiology Institute, Guangxi Academy of Agricultural Sciences, Nanning, Guangxi, China 4College of Agriculture, Guangxi University, Nanning, Guangxi, China

Edited by: Jesús Navas-Castillo, IHSM La Mayora, CSIC, Spain

Reviewed by: Anna Pytlak, Polish Academy of Sciences, Poland; César Marín, Santo Tomás University, Chile

 *Correspondence: Yang-Rui Li liyr@gxaas.net Xiu-Peng Song xiupengsong@163.com

†These authors have contributed equally to this work

 15 08 2023 2023 14 1132016 26 12 2022 26 07 2023 Copyright © 2023 Guo, Li, Yang, Verma, Singh, Singh, Khan, Sharma, Qin, Zhang, Song and Li. 2023 Guo, Li, Yang, Verma, Singh, Singh, Khan, Sharma, Qin, Zhang, Song 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.

Sugarcane is an important sugar and energy crop worldwide, requiring a large amount of nitrogen (N). However, excessive application of synthetic N fertilizer causes environmental pollution in farmland. Endophytic nitrogen-fixing bacteria (ENFB) provide N nutrition for plants through biological N fixation, thus reducing the need for chemical fertilizers. The present study investigated the effect of the N-fixing endophytic strain Enterobacter roggenkampii ED5 on phytohormone indole-3-acetic acid (IAA), N-metabolism enzyme activities, microbial community compositions, and N cycle genes in sugarcane rhizosphere soil at different N levels. Three levels of 15N-urea, such as low N (0 kg/ha), medium N (150 kg/ha), and high N (300 kg/ha), were applied. The results showed that, after inoculating strain ED5, the IAA content in sugarcane leaves was significantly increased by 68.82% under low N condition at the seedling stage (60 days). The nitrate reductase (NR) activity showed a downward trend. However, the glutamine synthase (GS) and NADH-glutamate dehydrogenase (NADH-GDH) activities were significantly enhanced compared to the control under the high N condition, and the GS and NR genes had the highest expression at 180 and 120 days, respectively, at the low N level. The total N content in the roots, stems, and leaves of sugarcane was higher than the control. The 15N atom % excess of sugarcane decreased significantly under medium N condition, indicating that the medium N level was conducive to N fixation in strain ED5. Metagenome analysis of sugarcane rhizosphere soil exhibited that the abundance of N-metabolizing microbial richness was increased under low and high N conditions after inoculation of strain ED5 at the genus level, while it was increased at the phylum level only under the low N condition. The LefSe (LDA > 2, p < 0.05) found that the N-metabolism-related differential microorganisms under the high N condition were higher than those under medium and low N conditions. It was also shown that the abundance of nifDHK genes was significantly increased after inoculation of ED5 at the medium N level, and other N cycle genes had high abundance at the high N level after inoculation of strain ED5. The results of this study provided a scientific reference for N fertilization in actual sugarcane production.

 Enterobacter roggenkampii ED5 sugarcane growth nitrogen fixation rhizosphere metagenomics nitrogen levels section-at-acceptance Microbe and Virus Interactions with Plants Introduction

Sugarcane (Saccharum spp.) is the leading sugar and energy crop in China, and its growth requires a large amount of nitrogen (N) (Li and Yang, 2015). N is a nutrient element of many important organic compounds in plants and significantly impacts plant life activities, crop yield, quality, and nutrient accumulation (Novoa and Loomis, 1981; Kraiser et al., 2011). Nowadays, the N nutrition of crops mainly comes from chemical N fertilizers. However, excessive N application also malnourishes soil at a global scale and causes leaching and eutrophication (Zhu and Chen, 2002; Savci, 2012; Qiao et al., 2018). Biological N fixation (BNF) by N-fixing bacteria converts atmospheric nitrogen into ammonia ions available to plants with the action of nitrogenase. The amount of N from BNF in sugarcane reached 40 kg N ha−1 annually in Brazil (Urquiaga et al., 2012). Therefore, utilizing BNF can ensure the ecologically sustainable development of sugarcane production. Some endophytic N-fixing strains have already been screened from sugarcane tissues, such as Beijerinckia, Novosphingobium sediminicola, Pantoea sp., Klebsiella, Enterobacter sacchari sp., Klebsiella variicola, Gluconacetobacter diazotrophicus, and Herbaspirillum sp. (Dobereiner, 1961; Loiret et al., 2004; Njoloma et al., 2006; Bertalan et al., 2009; Lin et al., 2012b; Zhu et al., 2013; Wei et al., 2014; Muangthong et al., 2015).

The 15N isotope dilution technique is commonly used to detect the effects of associated N-fixing bacteria and plants. It was reported that the inoculation of the N-fixing bacteria SC20 contributed 28% of the total N in sugarcane plants determined by the 15N isotope dilution technique (Sajjad Mirza et al., 2001). The strains Bacillus megaterium CY5 and Bacillus mycoides CA1 were found to have an apparent N fixation role in sugarcane by using the 15N isotope dilution technique (Singh et al., 2020b). The effect of N-fixing strains including Microbacterium sp. 16SH, Pseudomonas koreensi CY4, Kosakonia radicincitans BA1, and Stenotrophomonas maltophilia COA2 on sugarcane has been investigated in China using the 15N isotope dilution technique (Lin et al., 2012a; Li et al., 2017; Singh et al., 2020a).

For metagenomics technology, the DNA of microorganism samples taken from the environment is directly extracted to construct a metagenomic library. Genomics research methods are applied to study the genetic composition and community functions of microorganisms contained in environmental samples to reveal the micro-ecological mechanisms in the particular environment at the three levels of community, gene, and function (Hugenholtz and Tyson, 2008; Suenaga, 2012; Escobar-Zepeda et al., 2015; Awasthi et al., 2020). Using metagenomics technology, Pang et al. (2021) found that continuous cropping of sugarcane could change the composition of soil bacteria and fungi, resulting in a reduction of N, total sulfur, cane yield, and sugar content.

Application of an appropriate amount of N fertilizer can effectively regulate the structure and number of soil microbiota, which could support the stability and sustainability of the rhizosphere ecosystem (Dixon and Kahn, 2004; Geisseler and Scow, 2014; Zhao et al., 2019; Hu et al., 2022). In a study, four levels of N with urea, such as zero N (0 kg ha−1), low N (96 kg ha−1), medium N (482 kg ha−1), and high N (964 kg ha−1), were applied, but only the low N application significantly increased the richness and diversity of soil bacteria in sugarcane fields (Yang et al., 2021). In another study, medium N application (150 N kg ha−1) improved the soil bacterial community in general, and Acinetobacter and Proteobacteria phyla, in particular, were enhanced by 47 and 71%, respectively, compared to the control (Khan et al., 2022).

Our research team isolated the N-fixing strain Enterobacter roggenkampii ED5 (strain ED5) from sugarcane roots (Guo et al., 2020). It was found that, after inoculation with the endophytic diazotroph strain ED5, the physiological enzymes and plant hormone levels in sugarcane were significantly enhanced, and some growth-related key genes were enriched at the transcriptional level (Guo et al., 2021, 2022). The rational application of N fertilizer is important in sugarcane production. Using 15N isotope dilution technology, we preliminarily proved that strain ED5 has a strong N fixation ability associated with sugarcane. However, it is unknown whether different N levels affect the growth-promoting and N-fixing capabilities of strain ED5 and whether the intervention of ED5 inoculation under different levels of N affects the microecology of the sugarcane rhizosphere.

In this study, 15N-urea was used to set up three different N application levels, including low, medium, and high N. Strain ED5 was inoculated in sugarcane at three N levels. The study included (i) analyzing the indole-3-acetic acid (IAA) content in sugarcane, (ii) detecting the N-metabolism-related enzyme activities and gene expressions in sugarcane leaves, (iii) checking the biological N fixation characteristics of sugarcane, and (iv) studying the N-metabolism-related microbial community structure and functional composition and the genes in the N cycle process in sugarcane rhizosphere soil. This aim was to provide a reference for the application of N-fixing bacteria in sugarcane production.

 Materials and methods

The endophytic diazotroph strain ED5 isolated and screened from sugarcane roots (Guo et al., 2020) was used in this study. The low N fixing sugarcane variety GT11 (more N fertilizer is needed for growth than other varieties) was selected. The pot size for the greenhouse experiment was 32 cm in upper diameter, 25 cm in lower diameter, and 27 cm in height. The N fertilizer 15N-urea with an abundance of 10.09% used in the experiment was produced by the Shanghai Research Institute of Chemical Industry Co., Ltd., China. The greenhouse experiment was done at the Sugarcane Research Institute of the Guangxi Academy of Agricultural Sciences (GXAAS) with Latitude 108.09834°, Longitude 22.90991°, and altitude 89 m in Nanning, China. The experiment soil was medium-fertility red soil from sugarcane fields in GXAAS.

Preparation of inoculum liquid and sugarcane seedling transplanting

Strain ED5 was inoculated in LB broth, incubated at 180 rpm and 32°C for 48 h, and centrifuged at 5,000 rpm for 5 min, and the supernatant was discarded. The bacteria were diluted and mixed with sterilized water so that the final concentration of bacteria reached 106 cfu mL−1 in the inoculum suspension. Sugarcane for the experiment was planted using the seedling tray method. The length and width of the seedling tray were 40 cm and 30 cm, respectively, and the culture substrate was sterilized by high-pressure steam sterilization. The seed cane was cut into bud segments of the same length, soaked in 1% carbendazim solution for 0.5 h, and then cultured in the tray. The seedlings were transplanted to the greenhouse when they grew to 2–3 leaves and were managed routinely, and those with the same growth momentum were selected for inoculation and transplantation. The root surface of the selected seedlings was washed with sterile water, soaked in the bacteria suspension for 1 h, and then transplanted to the experimental pots. Those soaked in sterile water were used as controls.

 Greenhouse experiment

In the greenhouse experiment, the pot-cultured plants were used to apply 15N-urea fertilizer at once after transplanting. Three N application levels were set at low N (0 g 15N-urea pot, equivalent to N 0 kg/ha in the field), medium N (1.52 g15N-urea per pot, equivalent to N 150 kg/ha in the field), and high N (3.05 g 15N-urea per pot, equivalent to N 300 kg/ha urea in the field). The plants were divided into treatments and control groups. The treatment groups were designed as follows: strain ED5 inoculated at low N (GLE), non-inoculated at low N (GL); strain ED5 inoculated at medium N (GME), non-inoculated at medium N (GM); and strain ED5 inoculated at high N (GHE), non-inoculated at high N (GH). A completely randomized block design was adopted with six pots for each treatment, and a total of 108 pots were used for all the treatments.

 IAA content

The IAA content in leaf +1 (Cheavegatti-Gianotto et al., 2011) was determined by the enzyme-linked immunosorbent assay (ELISA) kit produced by Jiangsu Enzymic Immune Industry Co., Ltd., China. The quantity of leaf +1 was weighted, PBS buffer (pH 7.4) was added, and it was ground evenly with a grinder three times, 60 s each time. The product was centrifuged at 4°C at 3,000 rpm for 20 min, and then the supernatant was collected for testing. For standard samples, 50 μl of standard products were injected into each standard well, and a standard curve was drawn. Samples and blank wells were set in the ELISA plate, and 40 μl of diluent and 10 μl of the sample were added for testing. After 100 μl of enzyme-labeled reagent was added to each sample well, the plate was sealed and incubated at 37°C for 60 min. Then the plate sealing film was removed, the liquid was discarded, the washing solution was added and kept for 30 s before the solution was removed, and the wash was repeated five times. The chromogenic reagent was added to each well and cultured in the dark at 37°C for 15 min, then 50 μl of stop solution was added to terminate the reaction, and the OD value was detected with a microplate reader (Spark) at a wavelength of 450 nm. The proteins were quantitated using the BCA kit (Solarbio, China). All the tests were performed with three biological replicates.

 Activities of N-metabolism-related enzymes

The samples of sugarcane leaf +1 were taken at 60, 120, 180, and 240 days after strain ED5 inoculation. The activities of three N-metabolism-related key enzymes, namely, nitrate reductase (NR), glutamine synthase (GS), and NADH-glutamate dehydrogenase (NADH-GDH), were analyzed by kits produced by Jiangsu Edison Biotechnology Co., Ltd., China. All the tests were performed with three biological replicates.

 Expression of N-metabolism-related genes

The expressions of nitrite reductase (NR) and glutamine synthase (GS) genes in sugarcane leaf +1 were analyzed by quantitative real-time PCR (qRT-PCR) at 60, 120, 180, and 240 days, respectively, after strain ED5 inoculation. The primers for gene amplification are shown in Supplementary Table S1 (Yang et al., 2019). The total RNA in sugarcane leaves was extracted with the plant whole RNA extraction kit (TSP412, Qingke Biology, China), and its integrity and purity were detected by 1.2% agarose gel electrophoresis (1.8 ≤ OD260/OD280 ≤ 2.1). The Goldenstar RT6 cDNA Synthesis Kit version 2, a reverse transcription kit, was used for RNA reverse transcription amplification. The cDNA products obtained from reverse transcription were diluted with qPCR templates appropriately. The 2×T5 Fast qPCR Mix (SYBR Green I, Qingke Biology, China) was used as the qPCR enzyme, and the expression was detected with a LightCycler 1480 II real-time fluorescence quantitative PCR instrument. The qRT-PCR reaction components are listed in Supplementary Table S2. The 2−ΔΔCt method was used to calculate the relative gene expression (Livak and Schmittgen, 2001).

 Total N content and <sup>15</sup>N atom %

After inoculation of strain ED5, the total N content and 15N atom % in sugarcane roots, stems, and leaves were tested at 60, 120, 180, and 240 days, respectively, after strain ED5 inoculation. The whole sugarcane plants were taken at each testing stage after transplanting, and the soil and impurities on the surface of the roots, stems, and leaves were removed by washing with tap water. The cleaned samples were treated at 105°C for 30 min and then dried at 65°C. The dried sample was crushed into powder for analysis. According to Sáez-Plaza et al. (2013), the Kjeldahl method was used for the total N content measurement. The 15N atom % (N enrichment, 15N atom) analysis was completed by the Institute of Genetic Physiology, Hebei Academy of Agriculture and Forestry Sciences, China. The 15N atom was calculated according to the method suggested by Oliveira et al. (2006). All the tests were performed with three biological replicates.

 Sugarcane rhizosphere soil DNA extraction

Sugarcane rhizosphere soil sample DNA was extracted using the kit method (FastDNA® Spin Kit, MP Biomedicals, United States) at 120 days after inoculating strain ED5. The DNA integrity was checked by 1% agarose gel electrophoresis, and its purity and concentration were analyzed using a NanoDrop2000 (Thermo Fisher Scientific) and TBS-380 (TurnerBioSystems).

 PE library construction and metagenomic sequencing

The NEXTFLEX Rapid DNA Seq Kit (Bio Scientific) was used to construct the PE library. The library molecules were amplified with bridge PCR to generate DNA clusters, and the Illumina NovaSeq second-generation sequencing platform was used for macrogenome sequencing. The raw metagenomic sequencing data have been submitted to the NCBI SRA database with the accession number PRJNA911313.

 Sequence quality control and genome assembly

The fastp software (version 0.20.0) was used to cut the paired-end adapter sequence of the original data, remove low-quality reads, and obtain high-quality reads (Chen et al., 2018). The results are shown in Supplementary Table S3. The optimized quality control sequences were spliced and assembled into contigs (single assembly) using the MEGAHIT (version 1.1.2) software. The final assembly results were screened for contigs ≥300 bp. The assembly results are listed in Supplementary Table S4.

 Gene prediction and construction of a non-redundant gene catalog

The MetaGene software (http://metagene.cb.k.u-tokyo.ac.jp) was used to predict the open reading frames (ORFs) of the contigs in the spliced results (Noguchi et al., 2006), and the gene prediction results are shown in Supplementary Table S5. The genes were selected and translated into amino acid sequences, and the nucleic acid length was equal to or >100 bp. The predicted gene sequences of all the samples (90%, identity; 90%, coverage) were clustered by CD-HIT software (http://www.bioinformatics.org/cd-hit/, version 4.6.1), and the longest gene of each class was taken as the representative sequence to construct a non-redundant gene catalog (Fu et al., 2012). The gene number and length after removing redundancy are listed in Supplementary Table S6.

 Bioinformatics and statistical analysis

Diamond software (http://www.diamondsearch.org/index.php, version 0.8.35) was used to compare the amino acid sequence of the non-redundant gene catalog in the NR, eggNOG (version 4.5.1), COG, and KEGG (version 94.2) databases (BLASTP E-value ≤ 1e-05) (Buchfink et al., 2015). The number of genes belonging to bacteria was used to calculate the bacterial abundance at different taxonomic levels such as domain, kingdom, phylum, class, order, family, genus, and species, and the table of taxonomic abundance levels was constructed accordingly. The gene catalog of N metabolism was constructed according to the KEGG database pathway (ko00910), and the total gene abundance corresponding to KO, EC, and Module was used to calculate the abundance of the corresponding functional categories.

All the data on IAA content, N-metabolism-related enzyme activities, and qRT-PCR were analyzed using Microsoft Excel 2013. The SPSS 22.0 software was used for data processing and statistical analysis, and Duncan's multiple range test was used to analyze variance and multiple comparisons (p = 0.05).

 Results Total N in sugarcane

The total N content was significantly increased in sugarcane roots, stems, and leaves inoculated with strain ED5 at all stages compared with the non-inoculated plants (Figure 1; Supplementary Table S7). At the medium N level, it showed no significant difference at 180 days but was significantly enhanced at other stages. However, the total N content in sugarcane tissue was significantly decreased at the high N level, except for the stems and leaves at 180 and 240 days (Supplementary Table S7). Overall, the total N content in sugarcane showed an increasing trend compared to the control.

Heatmap analysis of the trend of total nitrogen content changes in sugarcane GT11 inoculated with E. roggenkampii ED5. (A) Root (N%); (B) Stem (N%); (C) Leaf (N%). The color intensity in each panel of the heatmap showed the total nitrogen content of sugarcane GT11.

 <sup>15</sup>N atom %

The results of the 15N atom % in sugarcane are listed in Table 1. The 15N atom % in all the organs of sugarcane was significantly decreased compared to the control after strain ED5 inoculation at the medium N level, and the highest loss was observed at 60 days with 24.71% in roots, 26.69% in stems, and 22.01% in leaves, respectively (Table 1). At the high N level, there was no significant difference in 15N atom % in roots among different stages. However, it increased by 8.25 and 11.73% in stems and leaves, respectively, at 60 days, and there was no significant change at 240 days (Table 1). The results showed that the low-dose N fertilizer was more conducive to the BNF effects of strain ED5 on sugarcane, and the high dose of N fertilizer may inhibit the function of BNF.

Effect of E. roggenkampii ED5 on 15N atom % of sugarcane at different nitrogen application levels.

Sugarcane organs Time (days) Medium nitrogen High nitrogen
GM GME GH GHE
Root (15N atom %) 60 3.367 ± 0.234b 2.535 ± 0.184c 4.591 ± 0.249a 4.578 ± 0.210a
120 2.580 ± 0.099b 2.136 ± 0.158c 3.910 ± 0.015a 3.719 ± 0.065a
180 2.311 ± 0.503b 1.599 ± 0.072c 3.290 ± 0.011a 3.559 ± 0.045a
240 1.822 ± 0.040b 1.489 ± 0.104c 2.941 ± 0.049a 2.881 ± 0.040a
Stem (15N atom %) 60 3.774 ± 0.248c 2.88 ± 0.094d 4.826 ± 0.088b 5.224 ± 0.016a
120 2.709 ± 0.058b 2.559 ± 0.064c 4.304 ± 0.007a 4.277 ± 0.007a
180 2.278 ± 0.020c 1.788 ± 0.028d 4.166 ± 0.053a 3.828 ± 0.076b
240 1.892 ± 0.047b 1.656 ± 0.072c 2.838 ± 0.041a 2.676 ± 0.109a
Leaf (15N atom %) 60 3.539 ± 0.099c 2.760 ± 0.136d 4.526 ± 0.081b 5.057 ± 0.057a
120 2.657 ± 0.063b 2.445 ± 0.105c 4.183 ± 0.003a 4.182 ± 0.008a
180 2.156 ± 0.062c 1.650 ± 0.087d 3.878 ± 0.030a 3.637 ± 0.082b
240 1.938 ± 0.053b 1.717 ± 0.170c 3.114 ± 0.023a 3.081 ± 0.046a

The data are presented as mean ± SE (n = 3, Duncan's multiple range test). The means with the same alphabets in the same line are not significantly different at p < 0.05.

 IAA content in sugarcane leaves

The IAA content in sugarcane leaves after inoculation with strain ED5 showed a diverse trend (Figure 2). At the low N level, the IAA content increased by 69.82 and 23.91% at 60 days and 180 days, respectively, but decreased dramatically at 120 and 240 days. At the medium N level, the IAA content increased by 11.78 and 136.59% at 60 and 120 days but decreased at 180 and 240 days. The IAA content at low and medium N application levels has the same trend (Figure 2A). The results of this study indicated that the IAA content in sugarcane at the seedling stage (60 days) increases at different N levels, so it is suggested that the sugarcane seedling stage may be the key stage for strain ED5 to promote plant growth.

Effect of E. roggenkampii ED5 on the IAA contents and activities of N-metabolism-related enzymes in leaves of sugarcane GT11 at different nitrogen application levels. (A) IAA; (B) NR; (C) GS; (D) NADH-GDH. The same letter indicates no significant difference was detected at Duncan's multiple range test, P ≤ 0.05 (n = 3).

 N-metabolism enzyme activities

This study assessed the activities of N-metabolism-related enzymes, NR, GS, and NADH-GDH (Figure 2). It was found that the NR activity showed an overall downward trend under different N levels at all the growth stages (Figure 2B). Similarly, the GS enzyme activity exhibited a downward trend at both low and medium N levels but an upward trend at the high N level compared to the control. At the high N level, the GS activity increased significantly by 42.14% at 60 days (Figure 2C). The NADH-GDH activity was significantly higher than that in the control, except for that at 240 days. The highest increase in NADH-GDH activity was 45.22% at 120 days at the low N level (Figure 2D). The NADH-GDH activity exhibited an opposite trend at medium and high N levels compared to the control, which was decreased at the low N level but enhanced at the high N level at 60 days (Figure 2D). Therefore, the high N condition significantly improved the GS activity in sugarcane after inoculating strain ED5. Therefore, it is important to improve the NADH-GDH activity and decrease the NR activity at low N conditions, especially at the seedling stage (60 days).

 Expression of NR and GS genes

The expressions of qRT-PCR analysis of GS and NR coding genes in sugarcane after inoculation with strain ED5 at different N levels are presented in Figure 3. It was found that the gene GS had the highest expression at 180 days. However, the expression of the gene NR was significantly increased at 120 days with the highest expression at the low N level (Figure 3). This showed that the strain ED5 might significantly impact sugarcane N metabolism under low N conditions.

Expression of NR and GS genes in leaves of sugarcane GT11 at different nitrogen application levels. (A) GS; (B) NR. The same letter indicates no significant difference was detected at Duncan's multiple range test, P ≤ 0.05 (n = 3).

 N-metabolism-related microbial community structure

The N-metabolism-related gene catalog was constructed according to the KEGG database pathway (ko00910), and the microbial annotation was performed in the NR database, with an abundance of < 0.01 classified as others. The results (Figure 4) showed that, in total, the microorganisms in the GEL and GHE groups were annotated in 59 phyla, including one unique phylum at the phylum level, while 61 phyla containing three unique phyla were annotated in the GH group and 58 phyla with no unique phylum were annotated in the GL, GM, and GME groups (Figure 4A). There were nine phyla in the rhizosphere soil of sugarcane at different N levels, which were Actinobacteria (35.89–37.93%), Proteobacteria (26.47–28.69%), Acidobacteria (10.84–11.26%), Chloroflexi (7.07–7.42%), Rokubacteria (3.23–3.76%), Gemmatimonadetes (2.92–3.33%), Planctomycetes (1.58–1.74%), Verrucomicrobia (1.10–1.27%), and Nitrospirae (0.71–1.02%) (Figure 4B). At the genus level, there were 17, 15, 14, 12, 12, and 21 unique genera in the GL, GLE, GM, GME, GH, and GHE groups, respectively (Figure 4C). It means there was higher microbial richness after inoculation with strain ED5 at a high N level. Additionally, 22 different genera related to N metabolism were detected, of which 14 belonged to unclassified genera, and the other eight genera were Nocardioides (4.53–5.25%), Solirubrobacter (2.01–2.44%), Gaiella (1.82–2.02%), Streptomyces (1.79–1.94%), Phycicoccus (1.42–1.52%), Pseudonocardia (0.78–1.03%), Bradyrhizobium (0.71–1.01%), and Steroidobacter (0.58–1.09%) (Figure 4D).

N-metabolism-related microbial community structures in sugarcane variety GT11 rhizosphere soil at different nitrogen application levels. (A, B) On the phylum level; (C, D), on the genus level.

 Differences in abundance of N-metabolism-related microorganisms

This study used the LefSe (LDA >2, p < 0.05) bar chart to analyze the N metabolism in various microorganisms from phylum to genus. A total of 135 differentially abundant taxonomic clades were identified, resulting in divergences in different treatments (Figure 5, Supplementary Table S8). The representative differentially abundant microbial taxa included Proteobacteria (LDA = 4.57, p = 0.048), Betaproteobacteria (LDA = 4.04, p = 0.038), Gammaproteobacteria (LDA = 3.82, p = 0.039), and Deltaproteobacteria (LDA = 3.66, p = 0.037), which were enriched in the GH group. In the GHE group, Nocardioides (LDA = 3.82, p = 0.028), Sphingomonadaceae (LDA = 3.26, p = 0.041), Verrucomicrobiales (LDA = 3.07, p = 0.033), and Mizugakiibacter (LDA = 2.93, p = 0.017) were found as the representative differentially microbial taxa. Differently abundant microorganisms, including Acidobacteria (LDA = 4.12, p = 0.041), Rhizobiales (LDA = 3.91, p = 0.048), Thermoleophilia (LDA = 3.80, p = 0.016), and Solirubrobacterales (LDA = 3.75, p = 0.018) were present in the GL group. Differentially abundant microbial taxa of Bradyrhizobiaceae (LDA = 3.45, p = 0.034) and Phycicoccus (LDA = 3.26, p = 0.043) were found in the GLE group. Solirubrobacter (LDA = 3.60, p = 0.030) and Nevskiales (LDA = 3.49, p = 0.035) were enriched in the GM group, and Burkholderiales (LDA = 3.42, p = 0.050) and Bdellovibrionales (LDA = 2.70, p = 0.017) were found in the GME group. These results indicated that there were more N-metabolism-related microbial abundant taxonomic clades in the rhizosphere soil of sugarcane under the high N condition than in medium and low N conditions, and those in the control group were higher than those in the strain ED5 inoculation group.

The LEfSe bar chart of N-metabolism-related differences in taxon abundance in the rhizosphere soil of sugarcane GT11 at different nitrogen levels.

 N-metabolism-related functional composition

The top 30 classes of KEGG and KO hierarchical N-metabolism-related functional composition in the rhizosphere soil of sugarcane were analyzed by heatmap. The results (Figure 6) showed that the abundance of glutamine synthetase (K01915), glutamate synthase (GOGAT) large chain (K00265), NR (K00370), and GDH (K15371) encoding genes was higher than that of other genes (Figure 6). However, the KO function of carbonic anhydrase (K01673) had lower gene abundance in the GM and GH groups, and hydroxylamine reductase (K15864) had the lowest gene abundance in all the treatment groups. Overall, the KO functional composition for N metabolism in the rhizosphere soil of sugarcane was similar in different treatment groups.

Heatmap analysis showing functions of N-metabolism-related microorganisms in the rhizosphere soil of sugarcane GT11 at different nitrogen application levels. The color intensity in each panel of the heatmap showed the functional abundance of KO.

 Comparison of N-cycle process genes

The abundance of the N-fixing genes nifDKH was analyzed in the rhizosphere soil of sugarcane at different N levels. The results (Figure 7) showed that the abundance of nifDKH genes in the group inoculated with strain ED5 was significantly higher than the control group by 571.41% at medium N level. However, other treatment groups had no significant difference in N-fixing gene abundance (Figure 7A). The abundance of ammonia monooxygenase-encoding genes AmoABC was significantly lower by 21.93% than that in the control group at the low N level. However, it was substantially higher by 15.35%, although it was not statistically significant at the medium N level (Figure 7B). The nitrate reductase gene NarGHI in sugarcane rhizosphere soil was significantly enhanced at medium and high N levels by 5.37 and 15.21%, respectively, and significantly reduced by 9.66% compared to the control at the low N level. The abundance of nitrate reductase (cytochrome)-encoding genes NapAB decreased considerably by 7.09% at medium N level and increased by 24.67% compared to control at low N level. The abundance of the nitrite reductase (NO−)-coding gene NiKS was significantly increased by 19.99% at the high N level. In contrast, there was no significant difference in low and medium N conditions compared to the control. Additionally, the abundance of the NO reductase subunit-encoding gene NorBC and the N2O reductase-encoding gene NosZ was significantly increased under the high N condition (Figure 7C). The assimilative nitrate reductase catalytic subunit encoding genes NasAB and nitrite ferredoxin reductase-encoding gene NirA were increased significantly, with increments of 10.18 and 8.55% at medium level and 22.99 and 38.58% at the high N level. The abundance of the ferredoxin-nitrate reductase-encoding gene NarB was significantly increased by 56.75% at the high N level as compared to the control (Figure 7D). Based on these results, the medium N condition might improve the N fixation activity in sugarcane rhizosphere soil microorganisms after strain ED5 inoculation. The high N condition had a significant impact on the processes of nitrification, denitrification, and nitrate assimilation and reduction to ammonia in sugarcane rhizosphere soil after inoculation of strain ED5.

The abundance of nitrogen fixation functional genes in the rhizosphere of sugarcane GT11. (A) Nitrogenase genes; (B) nitrification process genes; (C) denitrification process genes; (D) nitrate assimilation and reduction to ammonia-related genes. The same letter indicates no significant difference was detected at Duncan's multiple range test, P ≤ 0.05 (n = 3).

 Discussion

Nitrogen (N) nutrition is critical to improve sugarcane productivity. However, excessive application of chemical N fertilizer pollutes the environment and increases the cost of sugarcane production (Bordonal et al., 2018; Sanches and Otto, 2022). Utilizing the N-fixing function of microorganisms is an ecological way to reduce the amount of chemical N fertilizer (Shridhar, 2012). It was reported that the combined application of N fertilizer and Bacillus pumilus TUAT-1 improved rice yield and growth parameters and enhanced the ability of rice to withstand salinity stress (Win et al., 2022). Under different doses of N fertilizer, PGPB strains Azospirillum brasilense and Bacillus circulans were inoculated with stevia, which improved soil quality and nutrient uptake and significantly increased the content of steviol glycosides (Alotaibi et al., 2022). Combining N fertilizer and azotobacter can improve wheat yields more than azotobacter or N fertilizer alone (El-Sorady et al., 2022).

At present, the biological N fixation in rice, sugarcane, maize, soybean, wheat, and other crops has been studied extensively (Baldani et al., 2002; Ladha and Reddy, 2003; Ciampitti and Salvagiotti, 2018; Rosenblueth et al., 2018; Sheoran et al., 2021). The 15N isotope dilution technique is a common research method for biological N fixation in plants (Carranca et al., 1999). Five mungbean genotypes were inoculated with Bradyrhizobium spp. under greenhouse conditions, and it was found that different soil properties significantly affected the N fixation effect and plant growth (Diatta et al., 2020). Reinprecht et al. (2020) reported that applying N fertilizer reduced the symbiotic N fixation effect of Rhizobium leguminosarum and soybean. Using the 15N isotope dilution method, Hardarson et al. (1991) found that the N2 fixation rate was slightly reduced to 75% at the 100 kg N ha−1 application. In contrast, it will be further reduced to 60 and 43% at 200 and 400 kg N ha−1 application conditions, respectively. The N fixation characteristics of Lupinus angustifolius cv. Iliyarrie were evaluated on different N application levels of ammonium nitrate (60, 90, 100, 120, and 150 kg N ha−1), and it was shown that mineral N replaced biologically fixed N, thereby reducing the soil N use efficiency of Lupinus angustifolius cv. Iliyarrie N fixation with an average loss of 33 kg N ha−1 (Evans et al., 1987). The biomass and N use efficiency of maize and sorghum were significantly improved by inoculating PGPB under zero N fertilization; however, the dry stem quality and N fixation efficiency were enhanced under N fertilizer application (Aquino et al., 2021). In the present study, the percentage of total N in the roots, stems, and leaves of sugarcane inoculated with strain ED5 was significantly increased, and the 15N atom % was considerably lower than the control at medium N level, indicating that strain ED5 has a good function of N fixation in sugarcane.

The secretion of endogenous plant hormones by the PGPB is a direct cause of plant growth-promoting effects. The IAA content in sugarcane was significantly enhanced after inoculation with strains Streptomyces chartreusis WZS02, Pseudomonas aeruginosa B18, and Gluconacetobacter diazotrophicus PAL5, and plant growth was promoted (Ullah et al., 2019; Wang et al., 2019; Singh et al., 2021). The strain PS3 Rhodopseudomonas palustris was inoculated into Brassica rapa var. chinensis at hydroponic conditions with nitrate as the primary N source, which significantly increased the endogenous IAA level in the plant leaves (Hsu et al., 2021). Our previous studies found that strain ED5 can efficiently secrete IAA in vitro and reach 732.93 μg mL−1 under 5% tryptophan condition (Guo et al., 2020). In this study, the N participation significantly affected the changes in IAA content in sugarcane compared to the control, indicating that strain ED5 may regulate sugarcane growth by secreting IAA.

NR, GS, and NADH-GDH enzymes are vital for N metabolism in plants. In the present study, the activities of N-metabolism-related enzymes NR, NADH-GDH, and GS in sugarcane were significantly affected by inoculating strain ED5, among which the NR activity was increased considerably under medium N condition compared to the control at the late growth stage of sugarcane. In contrast, the GS activity was significantly higher than that in the control at different stages. dos Santos et al. (2020) found that the N-assimilating enzymes NR and GS activities were affected by inoculation with five kinds of N-fixing bacteria. The activity of GS was higher than that in the control. In contrast, the NR activity showed variations. The activities of enzymes related to cucumber N metabolism (NR, GS, and GDH) were increased after PGPB strain Bacillus tequilensis SX31 inoculation (Wang et al., 2022). Similar to this study, diazotrophic Paenibacillus beijingensis BJ-18 was inoculated in maize, cucumber, and wheat with different N fertilization levels, and the results showed that N uptake, nitrate transporter (NRT), and NiR (sub-NR), NR, GS, and GOGAT were upregulated by 1.5–91.9-fold (Li et al., 2019). The relative expression of NR and GS genes in sugarcane leaves analyzed using qRT-PCR showed that the effects of strain ED5 on gene expression were different at different N levels.

The soil N-metabolism cycle plays a vital role in the plant ecosystem, in which a large number of microorganisms and plants interact to catalyze the redox reaction of N compounds, and microbes are the main drivers for maintaining the N cycle in the rhizosphere of plants (Knops et al., 2002; Chapman et al., 2006; Isobe and Ohte, 2014; Coskun et al., 2017; Moreau et al., 2019; Ramm et al., 2022). The application of N fertilizer can change the microbial community in the plant rhizosphere, directly affecting the N cycle (Nemergut et al., 2008; Dincá et al., 2022). Some studies reported that microbial fertilizer preparations could effectively change the structure of soil microbial communities, improve the ecological environment of plant rhizosphere microorganisms, and indirectly promote plant growth (Chen et al., 2011; Elnahal et al., 2022). N addition was the main driving force of soil microecological change. The microbial diversity will be reduced by long-term N addition. The soil acidification caused by N application is the crucial reason for the decline of soil microbial diversity because bacteria are suitable for surviving in a neutral environment, and N addition reduces bacterial diversity (Dai et al., 2018; Yang et al., 2020). In this study, the microbial richness of N-metabolizing bacteria in the sugarcane rhizosphere was improved under the high N fertilization condition, indicating that the high N condition is more conducive to the survival of N-metabolism-related bacteria. In an experiment, tomatoes were inoculated with Bacillus pumilus (PGPB) and planted with native soil without N fertilization and 150 mg N kg−1 urea soil in pots, which increased plant growth, N uptake, soil NH4+ concentration, rhizosphere bacterial population, soil bacterial gene expression, and soil nitrogenase activity (Masood et al., 2020). According to previous reports, the contents of physiological and chemical factors such as available phosphorus, N, and pH in soil are related to the abundance of N-fixing bacteria (Orr et al., 2012; Che et al., 2018). The associative N-fixation process of N-fixing bacteria and plants requires more energy, and organic carbon is the main energy source for N-fixing bacteria in soil (Herbert, 1999). However, applying N fertilizer can increase the organic carbon content, so the abundance of N-fixing genes and the diversity of N-fixing bacteria in the soil could be increased (Rodrigues Coelho et al., 2008). N fertilizer addition increased the number and diversity of endophytic diazotrophic bacterial communities in maize, while there was no significant difference in the rhizosphere soil (Rodríguez-Blanco et al., 2015). After diazotrophic Paenibacillus triticisoli BJ-18 was inoculated in maize on different N fertilization levels, the abundance and diversity of the bacterial community and diazotrophic community were significantly higher in the rhizosphere than in roots and stems through 16S rRNA and nifH gene sequencing (Li et al., 2021b).

Excessive application of N fertilizer can reduce the abundance of N-fixing genes and microbial richness in crops (Wu, 2011) and the abundance of the nitrogenase gene nifH in upland sorghum (Rodrigues Coelho et al., 2008). It was also found that N fertilization did not affect the abundance of the nitrogenase nifH gene in the farmland soil (Ouyang et al., 2018). The effects of diazotrophic microbial taxa on the abundance, activity, diversity, and community composition of grass taxa (switchgrass and big bluestem) rhizosphere soil microbes were evaluated under different N application conditions (0, 67, and 202 kg N ha −1), and the results showed that the diazotrophic microbial taxa were more vital in grass green-up, initial harvest, and second harvest agricultural seasons. The abundance of nifH transcripts of diazotrophic microbial taxa was reduced at excessive fertilization (202 kg N ha−1) compared to moderate fertilization (67 kg N ha−1), which reflected that high N fertilization has an inhibitory effect on soil diazotrophic microbial richness (Hu et al., 2021). The diversity of diazobacteria in red soil with long-term N fertilization was analyzed by PCR-RFLP, and the results showed that the diversity of nirK and nosZ genes was lower than that of nifH, but the sensitivity of the nirK and nosZ genes was higher than that of the nifH gene under different N fertilization conditions (Zhang et al., 2021). Similarly, long-term N addition significantly reduced the abundance of nifH and 16S rDNA genes, and high doses of N and phosphorus fertilizers decreased the abundance of diazotrophic bacteria populations, but low amounts of N and phosphorus fertilizers did not reduce their abundance (Zhou et al., 2021). The structure and abundance of the endophytic bacterial community in sweet sorghum were found to be affected by the amplification of rrs and nifH genes based on PCR-DGGE, qPCR, and high-throughput sequencing, but the diversity was not affected (Mareque et al., 2018). The diazotrophic bacterium Paenibacillus triticisoli BJ-18 was continuously used as an inoculum in winter wheat planted for three consecutive years, and the results showed that the abundance of some nitrogenase genes was increased significantly by using metagenome sequencing technology (Li et al., 2021a). In this study, we found that the abundance of N-fixing functional genes nifDHK in sugarcane inoculated with strain ED5 was significantly increased at medium N levels, whereas there was no significant difference under low and high N conditions, indicating that appropriate application of N fertilizer is beneficial to improve the associative N-fixing bacteria in sugarcane.

Nitrification oxidizes ammonia into nitrate, which can weaken the N retention capacity of the soil and accelerate the loss of N in agricultural production through nitrate leaching. This process impacts the N uptake and N use efficiency in plants and affects the N cycle in the soil environment (Subbarao et al., 2015). Nitrification is a two-step process involving the oxidation of ammonia and nitrite by bacteria to eventually produce nitrate (Van Kessel et al., 2015). Ammonia monooxygenase plays a key catalytic role in the nitrification process, and its encoding genes are AmoABC. It has been reported that the net nitrification rate of N in vegetable rhizosphere soil was significantly increased under single N fertilization, and ammonia-oxidizing bacteria were correlated with the net nitrification rate (Bi et al., 2017). In this study, the abundance of AmoABC genes was significantly decreased after inoculation with strain ED5 in the rhizosphere of sugarcane under the low N condition, whereas it increased dramatically at the high N level. Rehr and Klemme (1989) reported that applying chemical fertilizers may lead to increased denitrification and reduced NO3--N to N2, cause N loss, and accumulate a large amount of nitrous oxide in the soil. In this study, the abundance of NarGHI, NapAB, NiKS, and NorBC genes related to the denitrification process in sugarcane rhizosphere was significantly increased after inoculation with strain ED5 under the high N condition, indicating that inoculation with strain ED5 might have induced the nitrification or denitrification process in sugarcane rhizosphere under the high N condition, which was not conducive to N fixation. The dissimilatory reduction of nitrate to ammonium (DNRA) uses two enzymes (NaR and NiR) to reduce nitrate to ammonium salts, producing cationic ammonium, which retains N in the soil due to its weak mobility. The bacteria catabolize the carbon source and have a higher utilization rate of electron acceptors (Kraft et al., 2014). In the present study, the abundance of nitrite reductase (NADH)-encoding genes NirBD increased significantly in the sugarcane rhizosphere in the high and medium N treatments inoculated with strain ED5, indicating that the role of strain ED5 can effectively alleviate the loss of N in the soil. When ammonium is limited in soil, some bacteria grow with nitrate or nitrite as the only N source and synthesize nitrate into ammonia with the nitrate/nitrite assimilation reductase system (NAS system), which is controlled by NAS gene cluster regulation (Kraft et al., 2014). In this study, the abundance of nitrate reductase catalytic subunit genes NasAB and nitrite ferredoxin reductase gene NirA were significantly increased after inoculation with strain ED5 at the high N level, indicating that strain ED5 accelerated the assimilation of nitrate under N fertilization condition, which is beneficial to improve the N use efficiency in sugarcane.

 Conclusion

In this study, the sugarcane variety GT11 was inoculated with E. roggenkampii ED5 at different N levels. The results showed that the IAA content in sugarcane significantly increased at the seedling stage, and the activities of N-metabolism-related enzymes showed a downward trend, but the activities of NADH and GS significantly increased under high N condition after ED5 inoculation. The total N content and N fixation efficiency in sugarcane were significantly improved after inoculating strain ED5 under low N condition. There were differences in N-metabolism-related microbial community compositions in sugarcane rhizosphere soil under low and high N application conditions after inoculation of strain ED5 at the genus level, and the N-metabolism-related microorganisms under high N condition were more different than those at medium and low N levels. The enzyme activity and abundance of N-metabolism-related genes in the sugarcane rhizosphere can be enhanced by inoculation with strain ED5, especially nitrogenase genes. It is concluded that the medium N condition was good for the biological N fixation of sugarcane with endophytic diazotroph strain ED5, which could be suggested for sugarcane production.

 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 below: https://www.ncbi.nlm.nih.gov/, PRJNA911313.

 Author contributions

D-JG, D-PL, and Y-RL: planned the proposal and experiments. Y-RL and X-PS: study and resources. D-JG, D-PL, YQ, BY, PS, and RS: accomplished the experiments. AS, QK, KV, B-QZ, and X-PS: data examination. D-JG, QK, and KV: writing the original manuscript. Y-RL: review and editing. All authors contributed to the article and approved the submitted version.

 Funding

This study was supported by the Guangxi Key Laboratory of Sugarcane Genetic Improvement (21-238-16-K-03-03), the Fund for Guangxi Innovation Teams of Modern Agriculture Technology (nycytxgxcxtd-2021-03-01), the Guangxi Science and Technology Major Project (GKAA2211 7002-7), the Guangxi Academy of Agricultural Sciences (GNKB2017028, 2021YT11, GNKB2018034, and GNK2021JM87), and the New Agriculture Research and Reform Practice (2020103).

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

 References Alotaibi M. O. Alotibi M. M. Eissa M. A. Ghoneim A. M. (2022). Compost and plant growth-promoting bacteria enhanced steviol glycoside synthesis in stevia (Stevia rebaudiana Bertoni) plants by improving soil quality and regulating nitrogen uptake. S. Afr. J. Bot. 151, 306–314. 10.1016/j.sajb.2022.10.010 Aquino J. P. A. Antunes J. E. L. Bonifácio A. Rocha S. M. B. Amorim M. R. Alcântara Neto F. . (2021). Plant growth-promoting bacteria improve growth and nitrogen metabolism in maize and sorghum. Theor. Exp. Plant phys. 33, 249–260. 10.1007/s40626-021-00209-x Awasthi M. K. Ravindran B. Sarsaiya S. Chen H. Wainaina S. Singh E. . (2020). Metagenomics for taxonomy profiling: tools and approaches. Bioengineered. 11, 356–374. 10.1080/21655979.2020.173623836154550 Baldani J. I. Reis V. M. Baldani V. L. Döbereiner J. (2002). A brief story of nitrogen fixation in sugarcane—reasons for success in Brazil. Funct. Plant Biol. 29, 417–423. 10.1071/PP0108332689486 Bertalan M. Albano R. de Pádua V. Rouws L. Rojas C. Hemerly A. . (2009). Complete genome sequence of the sugarcane nitrogen-fixing endophyte Gluconacetobacter diazotrophicus Pal5. BMC Genomics. 10, 450. 10.1186/1471-2164-10-45019775431 Bi Q.-F. Chen Q.-H. Yang X.-R. Li H. Zheng B.-X. Zhou W.-W. . (2017). Effects of combined application of nitrogen fertilizer and biochar on the nitrification and ammonia oxidizers in an intensive vegetable soil. AMB Express. 7, 198. 10.1186/s13568-017-0498-729116481 Bordonal R. d. O Carvalho J. L. N. Lal R. de Figueiredo E. B. de Oliveira B. G. . (2018). Sustainability of sugarcane production in Brazil. A review. Agron. Sustain. Dev. 38, 13. 10.1007/s13593-018-0490-x Buchfink B. Xie C. Huson D. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. 10.1038/nmeth.317625402007 Carranca C. De Varennes A. Rolston D. (1999). Biological nitrogen fixation by fababean, pea and chickpea, under field conditions, estimated by the 15N isotope dilution technique. Eur. J. Agron. 10, 49–56. 10.1016/S1161-0301(98)00049-5 Chapman S. K. Langley J. A. Hart S. C. Koch G. W. (2006). Plants actively control nitrogen cycling: uncorking the microbial bottleneck. New Phytol. 169, 27–34. 10.1111/j.1469-8137.2005.01571.x16390416 Che R. Deng Y. Wang F. Wang W. Xu Z. Hao Y. . (2018). Autotrophic and symbiotic diazotrophs dominate nitrogen-fixing communities in Tibetan grassland soils. Sci. Total Environ. 639, 997–1006. 10.1016/j.scitotenv.2018.05.23829929338 Cheavegatti-Gianotto A. De Abreu H. M. C. Arruda P. Bespalhok Filho J. C. Burnquist W. L. Creste . (2011). Sugarcane (Saccharum X officinarum): a reference study for the regulation of genetically modified cultivars in Brazil. Trop. Plant Biol. 4, 62–89. 10.1007/s12042-011-9068-321614128 Chen S. Gu J. Gao H. Qin Q. (2011). Effect of microbial fertilizer on microbial activity and microbial community diversity in the rhizosphere of wheat growing on the Loess Plateau. Afr. J. Microbiol. Res. 5, 137–143. Chen S. Zhou Y. Chen Y. Gu J. (2018). fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 34, i884–i890. 10.1093/bioinformatics/bty56030423086 Ciampitti I. A. Salvagiotti F. (2018). New insights into soybean biological nitrogen fixation. Agron. J. 110, 1185–1196. 10.2134/agronj2017.06.0348 Coskun D. Britto D. T. Shi W. Kronzucker H. J. (2017). How plant root exudates shape the nitrogen cycle. Trends Plant Sci. 661–673. 10.1016/j.tplants.2017.05.00428601419 Dai Z. Su W. Chen H. Barberán A. Zhao H. Yu M. . (2018). Long-term nitrogen fertilization decreases bacterial diversity and favors the growth of Actinobacteria and Proteobacteria in agro-ecosystems across the globe. Global Change Biol. 24, 3452–3461. 10.1111/gcb.1416329645398 Diatta A. A. Thomason W. E. Abaye O. Thompson T. L. Battaglia M. L. Vaughan L. J. . (2020). Assessment of nitrogen fixation by mungbean genotypes in different soil textures using 15N natural abundance method. J. Soil Sci. Plant Nut. 20, 2230–2240. 10.1007/s42729-020-00290-2 Dincá L. C. Grenni P. Onet C. Onet A. (2022). Fertilization and soil microbial community: a review. Appl. Sci-Basel. 12, 1198. 10.3390/app12031198 Dixon R. Kahn D. (2004). Genetic regulation of biological nitrogen fixation. Nat. Rev. Microbiol. 2, 621–631. 10.1038/nrmicro95415263897 Dobereiner J. (1961). Nitrogen-fixing bacteria of the genusBeijerinckia Derx in the rhizosphere of sugar cane. Plant Soil. 15, 211–216. 10.1007/BF01400455 dos Santos S. G. da Silva Ribeiro F. Alves G. C. Santos L. A. Reis V. M. (2020). Inoculation with five diazotrophs alters nitrogen metabolism during the initial growth of sugarcane varieties with contrasting responses to added nitrogen. Plant Soil. 451, 25–44. 10.1007/s11104-019-04101-1 Elnahal A. S. El-Saadony M. T. Saad A. M. Desoky E.-S. M. El-Tahan A. M. Rady M. M. . (2022). The use of microbial inoculants for biological control, plant growth promotion, and sustainable agriculture: A review. Eur. J. Plant pathol. 162, 759–792. 10.1007/s10658-021-02393-7 El-Sorady G. A. El-Banna A. A. Abdelghany A. M. Salama E. A. Ali H. M. Siddiqui M. H. . (2022). Response of bread wheat cultivars inoculated with azotobacter species under different nitrogen application rates. Sustainability-basel. 14, 8394. 10.3390/su14148394 Escobar-Zepeda A. Vera-Ponce de León A. Sanchez-Flores A. (2015). The road to metagenomics: from microbiology to DNA sequencing technologies and bioinformatics. Front. Genet. 6, 00348. 10.3389/fgene.2015.0034826734060 Evans J. O'Connor G. Turner G. Bergersen F. (1987). Influence of mineral nitrogen on nitrogen fixation by lupin (Lupinus angustifolius) as assessed by 15N isotope dilution methods. Field Crop. Res. 17, 109–120. 10.1016/0378-4290(87)90086-4 Fu L. Niu B. Zhu Z. Wu S. Li W. (2012). CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 28, 3150–3152. 10.1093/bioinformatics/bts56523060610 Geisseler D. Scow K. M. J. S. B. (2014). Long-term effects of mineral fertilizers on soil microorganisms–A review. Soil Biol. Biochem. 75, 54–63. 10.1016/j.soilbio.2014.03.023 Guo D.-J. Li D.-P. Singh R. K. Singh P. Sharma A. Verma K. K. . (2021). Differential protein expression analysis of two sugarcane varieties in response to diazotrophic plant growth-promoting endophyte Enterobacter roggenkampii ED5. Front. Plant Sci. 12, 727741. 10.3389/fpls.2021.72774134887881 Guo D.-J. Li D.-P. Singh R. K. Singh P. Verma K. K. Sharma A. . (2022). Comparative transcriptome analysis of two sugarcane varieties in response to diazotrophic plant growth promoting endophyte Enterobacter roggenkampii ED5. J. Plant Interact. 17, 75–84. 10.1080/17429145.2021.2012608 Guo D.-J. Singh R. K. Singh P. Li D.-P. Sharma A. Xing Y.-X. . (2020). Complete genome sequence of Enterobacter roggenkampii ED5, a nitrogen fixing plant growth promoting endophytic bacterium with biocontrol and stress tolerance properties, isolated from sugarcane root. Front. Microbiol. 11, 580081. 10.3389/fmicb.2020.58008133072048 Hardarson G. Danso S. Zapata F. Reichardt K. (1991). Measurements of nitrogen fixation in fababean at different N fertilizer rates using the 15N isotope dilution and “A-value” methods. Plant Soil. 131, 161–168. 10.1007/BF00009445 Herbert R. (1999). Nitrogen cycling in coastal marine ecosystems. Fems Microbiol. Rev. 23, 563–590. 10.1111/j.1574-6976.1999.tb00414.x10525167 Hsu S.-H. Shen M.-W. Chen J.-C. Lur H.-S. Liu C.-T. (2021). The photosynthetic bacterium Rhodopseudomonas palustris strain PS3 exerts plant growth-promoting effects by stimulating nitrogen uptake and elevating auxin levels in expanding leaves. Front. Plant Sci. 12, 573634. 10.3389/fpls.2021.57363433613595 Hu J. Richwine J. D. Keyser P. D. Li L. Yao F. Jagadamma S. . (2021). Nitrogen fertilization and native C4 grass species alter abundance, activity, and diversity of soil diazotrophic communities. Front. Microbiol. 12, 675693. 10.3389/fmicb.2021.67569334305840 Hu Q. Liu T. Ding H. Guo L. Li C. Jiang Y. . (2022). Application rates of nitrogen fertilizers change the pattern of soil organic carbon fractions in a rice-wheat rotation system in China. Agr. Ecosyst. Environ. 338, 108081. 10.1016/j.agee.2022.108081 Hugenholtz P. Tyson G. W. (2008). Metagenomics. Nature. 455, 481–483. 10.1038/455481a18818648 Isobe K. Ohte N. (2014). Ecological perspectives on microbes involved in N-cycling. Microbes Environ. 29, 4–16. 10.1264/jsme2.ME1315924621510 Khan A. Jiang H. Bu J. Adnan M. Gillani S. W. Hussain M. A. . (2022). Untangling the Rhizosphere bacterial community composition and response of soil physiochemical properties to different nitrogen applications in sugarcane field. Front. Microbiol. 13, 856078. 10.3389/fmicb.2022.85607835369493 Knops J. Bradley K. Wedin D. (2002). Mechanisms of plant species impacts on ecosystem nitrogen cycling. Ecol. Lett. 5, 454–466. 10.1046/j.1461-0248.2002.00332.x28308347 Kraft B. Tegetmeyer H. E. Sharma R. Klotz M. G. Ferdelman T. G. Hettich R. L. . (2014). The environmental controls that govern the end product of bacterial nitrate respiration. Science. 345, 676–679. 10.1126/science.125407025104387 Kraiser T. Gras D. E. Gutiérrez A. G. González B. Gutiérrez R. A. (2011). A holistic view of nitrogen acquisition in plants. J. Exp. Bot. 62, 1455–1466. 10.1093/jxb/erq42521239377 Ladha J. Reddy P. (2003). Nitrogen fixation in rice systems: state of knowledge and future prospects. Plant Soil. 252, 151–167. 10.1023/A:1024175307238 Li H.-B. Singh R. K. Singh P. Song Q.-Q. Xing Y.-X. Yang L.-T. . (2017). Genetic diversity of nitrogen-fixing and plant growth promoting Pseudomonas species isolated from sugarcane rhizosphere. Front. Microbiol. 8, 1268. 10.3389/fmicb.2017.0126828769881 Li Y. Li Q. Chen S. (2021b). Diazotroph Paenibacillus triticisoli BJ-18 drives the variation in bacterial, diazotrophic and fungal communities in the rhizosphere and root/shoot endosphere of maize. Int. J. Mol. Sci. 22, 1460. 10.3390/ijms2203146033540521 Li Y. Li Y. Zhang H. Wang M. Chen S. (2019). Diazotrophic Paenibacillus beijingensis BJ-18 provides nitrogen for plant and promotes plant growth, nitrogen uptake and metabolism. Front. Microbiol. 10, 1119. 10.3389/fmicb.2019.0111931191471 Li Y. Wang M. Chen S. (2021a). Application of N2-fixing Paenibacillus triticisoli BJ-18 changes the compositions and functions of the bacterial, diazotrophic, and fungal microbiomes in the rhizosphere and root/shoot endosphere of wheat under field conditions. Biol. Fert. Soils. 57, 347–362. 10.1007/s00374-020-01528-y Li Y.-R. Yang L.-T. (2015). Sugarcane agriculture and sugar industry in China. Sugar Tech. 17, 1–8. 10.1007/s12355-014-0342-1 Lin L. Guo W. Xing Y. Zhang X. Li Z. Hu C. . (2012b). The actinobacterium Microbacterium sp. 16SH accepts pBBR1-based pPROBE vectors, forms biofilms, invades roots, and fixes N2 associated with micropropagated sugarcane plants. Appl. Microbiol. Biot. 93, 1185–1195. 10.1007/s00253-011-3618-322002067 Lin L. Li Z. Hu C. Zhang X. Chang S. Yang L. . (2012a). Plant growth-promoting nitrogen-fixing enterobacteria are in association with sugarcane plants growing in Guangxi, China. Microbes Environ. 27, 391–398. 10.1264/jsme2.ME1127522510648 Livak K. J. Schmittgen T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 25, 402–408. 10.1006/meth.2001.126225558171 Loiret F. Ortega E. Kleiner D. Ortega-Rodés P. Rodes R. Dong Z. (2004). A putative new endophytic nitrogen-fixing bacterium Pantoea sp. from sugarcane. J. Appl. Microbiol. 97, 504–511. 10.1111/j.1365-2672.2004.02329.x15281930 Mareque C. da J appl microbiol Silva T. F. Voll,ú R. E. Beracochea M. Seldin L. Battistoni F. (2018). The endophytic bacterial microbiota associated with sweet sorghum (Sorghum bicolor) is modulated by the application of chemical N fertilizer to the field. Int. J. Genomics. 2018, 7403670. 10.1155/2018/740367030363992 Masood S. Zhao X. Q. Shen R. F. (2020). Bacillus pumilus promotes the growth and nitrogen uptake of tomato plants under nitrogen fertilization. Sci. Hortic-Amsterdam. 272, 109581. 10.1016/j.scienta.2020.109581 Moreau D. Bardgett R. D. Finlay R. D. Jones D. L. Philippot L. (2019). A plant perspective on nitrogen cycling in the rhizosphere. Funct. Ecol. 33, 540–552. 10.1111/1365-2435.13303 Muangthong A. Youpensuk S. Rerkasem B. (2015). Isolation and characterisation of endophytic nitrogen fixing bacteria in sugarcane. Trop. Life Sci. Res. 26, 41–51.26868592 Nemergut D. R. Townsend A. R. Sattin S. R. Freeman K. R. Fierer N. Neff J. C. . (2008). The effects of chronic nitrogen fertilization on alpine tundra soil microbial communities: implications for carbon and nitrogen cycling. Environ. Microbiol. 10, 3093–3105. 10.1111/j.1462-2920.2008.01735.x18764871 Njoloma J. Tanaka K. Shimizu T. Nishiguchi T. Zakria M. Akashi R. . (2006). Infection and colonization of aseptically micropropagated sugarcane seedlings by nitrogen-fixing endophytic bacterium, Herbaspirillum sp. B501gfp1. Biol. Fert. Soils. 43, 137–143. 10.1007/s00374-006-0078-5 Noguchi H. Park J. Takagi T. (2006). MetaGene: prokaryotic gene finding from environmental genome shotgun sequences. Nucleic Acids Res. 34, 5623–5630. 10.1093/nar/gkl72317028096 Novoa R. Loomis R. (1981). Nitrogen and plant production. Plant Soil. 58, 177–204. 10.1007/BF02180053 Oliveira A. L. M. de Canuto E. L. Urquiaga S. Reis V. M. Baldani J. I. (2006). Yield of micropropagated sugarcane varieties in different soil types following inoculation with diazotrophic bacteria. Plant Soil. 284, 23–32. 10.1007/s11104-006-0025-0 Orr C. H. Leifert C. Cummings S. P. Cooper J. M. (2012). Impacts of organic and conventional crop management on diversity and activity of free-living nitrogen fixing bacteria and total bacteria are subsidiary to temporal effects. PLoS ONE. 7, e52891. 10.1371/journal.pone.005289123285218 Ouyang Y. Evans S. E. Friesen M. L. Tiemann L. K. (2018). Effect of nitrogen fertilization on the abundance of nitrogen cycling genes in agricultural soils: a meta-analysis of field studies. Soil Biol. Biochem. 127, 71–78. 10.1016/j.soilbio.2018.08.024 Pang Z. Dong F. Liu Q. Lin W. Hu C. Yuan Z. (2021). Soil metagenomics reveals effects of continuous sugarcane cropping on the structure and functional pathway of rhizospheric microbial community. Front. Microbiol. 12, 627569. 10.3389/fmicb.2021.62756933746921 Qiao C. Xu B. Han Y. Wang J. Wang X. Liu L. . (2018). Synthetic nitrogen fertilizers alter the soil chemistry, production and quality of tea. A meta-analysis. Agron. Sustain. Dev. 38, 1–10 10.1007/s13593-017-0485-z Ramm E. Liu C. Ambus P. Butterbach-Bahl K. Hu B. Martikainen P. J. . (2022). A review of the importance of mineral nitrogen cycling in the plant-soil-microbe system of permafrost-affected soils—changing the paradigm. Environ. Res. Lett. 17, 013004. 10.1088/1748-9326/ac417e Rehr B. Klemme J.-H. (1989). Competition for nitrate between denitrifying Pseudomonas stutzeri and nitrate ammonifying enterobacteria. Fems Microbiol. Ecol. 5, 51–57. 10.1111/j.1574-6968.1989.tb03657.x2619287 Reinprecht Y. Schram L. Marsolais F. Smith T. H. Hill B. Pauls K. P. (2020). Effects of nitrogen application on nitrogen fixation in common bean production. Front. Plant Sci. 11, 1172. 10.3389/fpls.2020.0117232849727 Rodrigues Coelho M. R. De Vos M. Carneiro N. P. Marriel I. E. Paiva E. Seldin L. (2008). Diversity of nifH gene pools in the rhizosphere of two cultivars of sorghum (Sorghum bicolor) treated with contrasting levels of nitrogen fertilizer. Fems Microbiol. Lett. 279, 15–22. 10.1111/j.1574-6968.2007.00975.x18070072 Rodríguez-Blanco A. Sicardi M. Frioni L. (2015). Plant genotype and nitrogen fertilization effects on abundance and diversity of diazotrophic bacteria associated with maize (Zea mays L.). Biol. Fert. Soils. 51, 391–402. 10.1007/s00374-014-0986-8 Rosenblueth M. Ormeño-Orrillo E. López-López A. Rogel M. A. Reyes-Hernández B. J. Martínez-Romero J. C. . (2018). Nitrogen fixation in cereals. Front. Microbiol. 9, 1794. 10.3389/fmicb.2018.0179430140262 Sáez-Plaza P. Michałowski T. Navas M. J. Asuero A. G. Wybraniec S. (2013). An overview of the Kjeldahl method of nitrogen determination. Part I. Early history, chemistry of the procedure, and titrimetric finish. Crit. Rev. Anal. Chem. 43, 178–223. 10.1080/10408347.2012.751786 Sajjad Mirza M. Ahmad W. Latif F. Haurat J. Bally R. Normand P. . (2001). Isolation, partial characterization, and the effect of plant growth-promoting bacteria (PGPB) on micro-propagated sugarcane in vitro. Plant Soil. 237, 47–54 10.1023/A:1013388619231 Sanches G. M. Otto R. (2022). A novel approach for determining nitrogen requirement based on a new agronomic principle—sugarcane as a crop model. Plant Soil. 472, 29–43. 10.1007/s11104-021-05263-7 Savci S. (2012). Investigation of effect of chemical fertilizers on environment. Apcbee Procedia 1, 287–292. 10.1016/j.apcbee.2012.03.047 Sheoran S. Kumar S. Kumar P. Meena R. S. Rakshit S. (2021). Nitrogen fixation in maize: breeding opportunities. Theor. Appl. Genet. 134, 1263–1280. 10.1007/s00122-021-03791-533677701 Shridhar B. S. (2012). Nitrogen fixing microorganisms. Int. J. Microbiol. Res. 3, 46–52. Singh P. Singh R. K. Guo D.-J. Sharma A. Singh R. N. Li D.-P. . (2021). Whole genome analysis of sugarcane root-associated endophyte Pseudomonas aeruginosa B18—A plant growth-promoting bacterium with antagonistic potential against Sporisorium scitamineum. Front. Microbiol. 12, 628376. 10.3389/fmicb.2021.62837633613496 Singh R. K. Singh P. Li H.-B. Guo D.-J. Song Q.-Q. Yang L.-T. . (2020a). Plant-PGPR interaction study of plant growth-promoting diazotrophs Kosakonia radicincitans BA1 and Stenotrophomonas maltophilia COA2 to enhance growth and stress-related gene expression in Saccharum spp. J. Plant Interact. 15, 427–445. 10.1080/17429145.2020.1857857 Singh R. K. Singh P. Li H.-B. Song Q.-Q. Guo D.-J. Solanki M. K. . (2020b). Diversity of nitrogen-fixing rhizobacteria associated with sugarcane: a comprehensive study of plant-microbe interactions for growth enhancement in Saccharum spp. BMC Plant Biol. 20, 1–21. 10.1186/s12870-020-02400-932423383 Subbarao G. V. Yoshihashi T. Worthington M. Nakahara K. Ando Y. Sahrawat K. L. . (2015). Suppression of soil nitrification by plants. Plant Sci. 233, 155–164. 10.1016/j.plantsci.2015.01.01225711823 Suenaga H. (2012). Targeted metagenomics: a high-resolution metagenomics approach for specific gene clusters in complex microbial communities. Environ. Microbiol. 14, 13–22. 10.1111/j.1462-2920.2011.02438.x21366818 Ullah S. Ashraf M. Asghar H. N. Iqbal Z. Ali R. (2019). Review Plant growth promoting rhizobacteria-mediated amelioration of drought in crop plants. Soil Envron. 38, 1–20. 10.25252/SE/19/71760 Urquiaga S. Xavier R. P. de Morais R. F. Batista R. B. Schultz N. Leite J. M. . (2012). Evidence from field nitrogen balance and 15N natural abundance data for the contribution of biological N2 fixation to Brazilian sugarcane varieties. Plant Soil. 356, 5–21. 10.1007/s11104-011-1016-3 Van Kessel M. A. Speth D. R. Albertsen M. Nielsen P. H. Op den Camp H. J. Kartal B. . (2015). Complete nitrification by a single microorganism. Nature. 528, 555–559. 10.1038/nature1645926610025 Wang J. Qu F. Liang J. Yang M. Hu X. (2022). Bacillus velezensis SX13 promoted cucumber growth and production by accelerating the absorption of nutrients and increasing plant photosynthetic metabolism. Sci. Hortic-Amsterdam. 301, 111151. 10.1016/j.scienta.2022.111151 Wang Z. Solanki M. K. Yu Z.-X. Yang L.-T. An Q.-L. Dong D.-F. . (2019). Draft genome analysis offers insights into the mechanism by which Streptomyces chartreusis WZS021 increases drought tolerance in sugarcane. Front. Microbiol. 9, 3262. 10.3389/fmicb.2018.0326230687260 Wei C.-Y. Lin L. Luo L.-J. Xing Y.-X. Hu C.-J. Yang L.-T. . (2014). Endophytic nitrogen-fixing Klebsiella variicola strain DX120E promotes sugarcane growth. Biol. Fert. Soils. 50, 657–666. 10.1007/s00374-013-0878-326203334 Win K. T. Oo A. Z. Yokoyama T. (2022). Plant Growth and Yield Response to Salinity Stress of Rice Grown under the Application of Different Nitrogen Levels and Bacillus pumilus Strain TUAT-1. Crops. 2, 435–444. 10.3390/crops2040031 Wu T. (2011). Can ectomycorrhizal fungi circumvent the nitrogen mineralization for plant nutrition in temperate forest ecosystems? Soil Biol. Biochem. 43, 1109–1117. 10.1016/j.soilbio.2011.02.003 Yang S. Xiao J. Liang T. He W. Tan H. (2021). Response of soil biological properties and bacterial diversity to different levels of nitrogen application in sugarcane fields. AMB Express. 11, 1–11. 10.1186/s13568-021-01331-434919198 Yang Y. Cheng H. Gao H. An S. (2020). Response and driving factors of soil microbial diversity related to global nitrogen addition. Land Degrad. Dev. 31, 190–204. 10.1002/ldr.3439 Yang Y. Gao S. Su Y. Lin Z. Guo J. Li M. . (2019). Transcripts and low nitrogen tolerance: Regulatory and metabolic pathways in sugarcane under low nitrogen stress. Environ. Exp. Bot. 163, 97–111. 10.1016/j.envexpbot.2019.04.010 Zhang C. Cao K. Li Y. Zhao J. Peng W. Cao H. . (2021). Long-term nitrogen fertilization shaped the nifH, nirK, and nosZ gene community patterns in red paddy soil in south China. Can. J. Microbiol. 67, 310–322. 10.1139/cjm-2020-030033022188 Zhao Z.-B. He J.-Z. Geisen S. Han L.-L. Wang J.-T. Shen J.-P. . (2019). Protist communities are more sensitive to nitrogen fertilization than other microorganisms in diverse agricultural soils. Microbiome. 7, 1–16. 10.1186/s40168-019-0647-030813951 Zhou J. Ma M. Guan D. Jiang X. Zhang N. Shu F. . (2021). Nitrogen has a greater influence than phosphorus on the diazotrophic community in two successive crop seasons in Northeast China. Sci. Rep. 11, 6303. 10.1038/s41598-021-85829-833737649 Zhu B. Zhou Q. Lin L. Hu C. Shen P. Yang L. . (2013). Enterobacter sacchari sp. nov., a nitrogen-fixing bacterium associated with sugar cane (Saccharum officinarum L.). Int. J. Syst. Evol. Micr. 63, 2577–2582. 10.1099/ijs.0.045500-023291881 Zhu Z. L. Chen D. (2002). Nitrogen fertilizer use in China–Contributions to food production, impacts on the environment and best management strategies. Nutr. Cycl. Agroecosys. 63, 117–127. 10.1023/A:1021107026067