Skip to main content


Suppressed N fixation and diazotrophs after four decades of fertilization



N fixation is one of the most important microbially driven ecosystem processes on Earth, allowing N to enter the soil from the atmosphere, and regulating plant productivity. A question that remains to be answered is whether such a fundamental process would still be that important in an over-fertilized world, as the long-term effects of fertilization on N fixation and associated diazotrophic communities remain to be tested. Here, we used a 35-year fertilization experiment, and investigated the changes in N fixation rates and the diazotrophic community in response to long-term inorganic and organic fertilization.


It was found that N fixation was drastically reduced (dropped by 50%) after almost four decades of fertilization. Our results further indicated that functionality losses were associated with reductions in the relative abundance of keystone and phylogenetically clustered N fixers such as Geobacter spp.


Our work suggests that long-term fertilization might have selected against N fixation and specific groups of N fixers. Our study provides solid evidence that N fixation and certain groups of diazotrophic taxa will be largely suppressed in a more and more fertilized world, with implications for soil biodiversity and ecosystem functions.


Biological nitrogen (N) fixation is one of the most important ecological processes on Earth, and is responsible for the fixation of up to 100 Tg N year−1 from the atmosphere globally [1,2,3]. However, N fixation and its associated microbial communities have been largely challenged by the industrial Haber process, and later by inorganic and organic fertilization [4], which provides 32 Tg N year−1 to global croplands [5]. Such a large amount of fertilization might relegate N fixers to a second place [6, 7] and could have long-term consequences for these important microbial communities and ecosystem processes in the future [8, 9]. Surprisingly, relatively little is known about the long-term effects of inorganic and organic fertilization on N fixation rates and their associated N fixers in terrestrial ecosystems.

Short-term additions of N fertilizers can result in an increase in the abundance of fast-growing diazotrophs [10]. These microbial communities may use resources from fertilizers to support their own vegetative growth, instead of fixing nitrogen [11], which is known to be an energy-expensive process [12]. Much less is known, however, about the long-term consequences (over decades) of soil fertilization (e.g., N fertilizer additions) on N fixation rates and their associated diazotrophic communities. We posit that, in fertilized environments, N fixation and fixers will become less and less important as time passes. However, experimental evidence supporting this hypothesis is lacking. Following natural selection theories [13,14,15], we hypothesized that fertilization should suppress N fixation and drastically change the community composition of N fixers, which may no longer be needed to fix N2 from the atmosphere. Fertilization could be particularly detrimental for oligotrophic microbial communities and for obligate N fixers that have limited ability to downregulate fixation. However, fertilization could benefit copiotrophic and facultative N fixers that are capable of downregulating fixation, such as Bradyrhizobium spp. [16, 17].

Here, we used soils from a 35-year fertilization experiment and the most advanced sequencing technology to target nifH genes that encode the reductase subunit of nitrogenase [18]. The role of fertilization in regulating N fixation and the phylogeny and community composition of N fixers were evaluated [19] by using contrasting fertilization management strategies: non-fertilization (control), chemical fertilization (NPK), chemical fertilization with wheat straw (NPK + WS), chemical fertilization with pig manure (NPK + PM), and chemical fertilization with cow manure (NPK + CM).


N fixation and N fixers under long-term fertilization scenarios

Our results indicated that N fixation rates were significantly suppressed by a wide range of fertilizers after almost four decades of fertilization (Fig. 1a). We found that N fixation rates dropped by 50%, which was more noticeable in bulk soils than in the rhizosphere (Additional file 1: Table S6). We then evaluated the effects of long-term fertilization on N fixers. To do so, we built a correlation network incorporating the detected dominant diazotrophic phylotypes and found three ecological clusters of N fixers strongly co-occurring with each other (modules #1, #2, and #3; Fig. 1b). Each ecological cluster consisted of multiple diazotrophic species attributing to different genera (Fig. 1c; Fig. 2a). Bradyrhizobium and Burkholderia were the most dominant genera of N fixers in module #1 and module #2; Geobacter and Anaeromyxobacter dominated module #3 (Fig. 1c). Long-term fertilization resulted in drastic changes in the relative abundance of ecological clusters; the relative abundance of module #3 was strongly reduced and that of modules #1 and #2 was somewhat increased—particularly under NPK + CM additions (Fig. 1d). Similar results were found for bulk and rhizosphere soil (Additional file 1: Table S8).

Fig. 1

a Nitrogen fixation rates in different fertilization treatments. b Network diagram with nodes colored according to each of the three main ecological clusters (modules #1–3). c Operational taxonomic unit (OTU) number properties of the dominant diazotrophic genera in the three main ecological clusters. d Relative abundance of the ecological clusters in different fertilization treatments. Different letters indicate the values that differ significantly among treatments at P < 0.05 (Duncan’s test) in bulk soil (A, B, C) and rhizosphere soil (a, b, c). BS: bulk soil, RS: rhizosphere soil. Fertilization treatments: control, non-fertilization; NPK fertilization, NPK (urea, superphosphate, and potassium chloride); NPK + WS, NPK with wheat straw; NPK + PM, NPK with pig manure; NPK + CM, NPK with cow manure (NPK + CM)

Fig. 2

a Phylogenetic tree displaying the taxonomic information on dominant soil diazotrophic phylotypes in three main ecological clusters (modules #1–3). b The standardized difference, in units of standard deviation (z-score), between observed and expected phylogenetic diversity assuming random sampling for each module. The dotted brown line represents the expected phylogenetic diversity for each treatment and the blue dashed lines represent the 95% confidence intervals. The points either above (2) or below (− 2) the blue lines represent the degree to which those ecological clusters are phylogenetically over-dispersed or clustered, respectively. Red, module #1; blue, module #2; cyan, module #3. BS: bulk soil, RS: rhizosphere soil. Fertilization treatments: control, non-fertilization; NPK fertilization, NPK (urea, superphosphate, and potassium chloride); NPK + WS, NPK with wheat straw; NPK + PM, NPK with pig manure; NPK + CM, NPK with cow manure (NPK + CM)

We also constructed a phylogenetic tree with the dominant diazotrophic phylotypes, and found that the N fixers of the Geobacter genus were highly clustered in module #3. On the contrary, N fixers within the dominant genera Bradyrhizobium and Burkholderia were randomly distributed in module #1 and module #2 (Fig. 2a). We then calculated the phylogenetic diversity of each ecological cluster and further compared these observed values with the expected random values for each cluster [20]. We found that the observed phylogenetic diversity for module #1 and module #2 was consistent with the random predictions (within the 95% confidence interval) across different fertilization treatments (Fig. 2b). However, the observed phylogenetic diversity for module #3 deviated significantly below the random prediction under long-term non-fertilization scenarios; this is indicative of phylogenetic clustering. Meanwhile, under long-term application of fertilizers, the trends for module #3 were indicative of phylogenetic randomness (Fig. 2b). These results suggest that long-term fertilization may have selected against the N fixers associated with module #3 (mostly Geobacter spp).

Linking N fixers to N fixation under long-term fertilization scenarios

A strong and significant positive association between the relative abundance of module #3 and N fixation rates was found (Fig. 3). However, no significant association was detected between the relative abundance of module #1, module #2, and N fixation rates (Fig. 3). Fifty diazotrophic phylotypes were detected to be highly positively correlated with nitrogen fixation rates using random forest regression (Additional file 2: Figure S2 and S3). These phylotypes were mostly included in module #3 (20/50) when compared to module #1 (3/50) and module #2 (4/50) (Additional file 1: Table S10).

Fig. 3

Regressions between the nitrogen fixation rates and the relative abundance of the main diazotrophic ecological clusters. From left to right are modules #1, #2, and #3, represented by red, blue, and cyan plots, respectively

Structural equation modeling (SEM) was then used to evaluate the potential associations between ecological clusters of N fixers and N fixation rates under different fertilization scenarios, and to gain a deeper knowledge of the indirect and direct effects of fertilization on N fixation when considering multiple environmental factors simultaneously. Our SEM explained 85% of the variation in N fixation rates (Fig. 4a). The relative abundance of module #1 and module #2 had direct negative effects on N fixation rates. However, a positive and significant association was still found between the relative abundance of module #3 and N fixation rates (Fig. 4). Thus, it indicated that long-term fertilization indirectly reduced N fixation by decreasing the relative abundance of diazotrophs within module #3. From a management perspective, the long-term negative effects of fertilization on N fixers in the module #3 seemed to be minimized when using NPK + CM as a fertilizer (Fig. 4a; box 2).

Fig. 4

a A structural equation model describing the effects of relative abundance of main ecological clusters on the nitrogen fixation rates. Soil properties include soil pH, total phosphorus (TP), total carbon (TC), and total nitrogen (TN). Numbers labeling the arrow lines are indicative of the correlations. R2 denotes the proportion of variance explained. Significance levels of each predictor are *P < 0.05 and **P < 0.01. b Standardized total effects (STE) from the SEM. This is the sum of direct and indirect effects from each variable on nitrogen fixation rate. Fertilization treatments: control, non-fertilization; NPK fertilization, NPK (urea, superphosphate, and potassium chloride); NPK + WS, NPK with wheat straw; NPK + PM, NPK with pig manure; NPK + CM, NPK with cow manure (NPK + CM)


Our work provides solid evidence that, after almost four decades of experiment, fertilization largely suppressed N fixation (about 50% decrease), and the relative abundance of specific N fixers (e.g., Geobacter spp.) that were reported to be positively associated with N fixation rates [21]. Our SEM provided further evidence that long-term fertilization indirectly reduced N fixation rates by decreasing the relative abundance of keystone and phylogenetically clustered N fixers within module #3 (those were positively associated with N fixation rates). In addition, fertilizations resulted in a change from phylogenetic clustering to phylogenetic randomness for the N fixers within module #3. These results suggest that long-term fertilization selects against N fixation and their associated N fixers. Our work has unveiled the basic mechanisms explaining the long-term effects of fertilization on N fixation and its associated microbial communities, and further suggest that the fundamental process of N fixation, and some keystone diazotrophs, will become increasingly suppressed as we continue to increase soil fertilization.

Our results identified a subset of positive (winners) and negative (losers) associations between N fixers and long-term fertilizations. For instance, long-term fertilization was positively associated with the relative abundance of modules #1 and #2 (winners under fertilization scenarios). Taxa within the dominant genera Bradyrhizobium and Burkholderia are known to be facultative N fixers [22, 23], which can consume soil resources from fertilization to support vegetative growth instead of fixing nitrogen [24]. Moreover, these genera are found within the classes Alphaproteobacteria and Betaproteobacteria, respectively, which are often classified as copiotrophs [25]. These dominant copiotrophic soil organisms may benefit from the resources coming from fertilizers and using them to support their own growth [10]. This idea is supported by the lack of correlation between the relative abundance of these taxa and N fixation rates. Our results suggest that taxa within modules #1 and #2 could benefit from long-term fertilization scenarios, providing a list of “winner” taxa.

On the other hand, long-term fertilization was negatively correlated with the relative abundance of diazotrophs within module #3 (losers under fertilization scenarios). Interestingly, fertilization also appeared to select against taxa within this ecological cluster, which were found to be phylogenetically clustered in the unfertilized field, but became phylogenetically randomness under long-term application of fertilizations. Reductions in the relative abundance of these taxa could also negatively influence N fixation rates. In fact, long-term fertilizations were found to be indirectly negatively associated with N fixation, by decreasing the relative abundance of N fixers within module #3. This indicates that taxa within module #3 will be suppressed under long-term fertilization with negative consequences for N fixation rates. In this respect, our study identified a list of “loser” taxa including Geobacter and Anaeromyxobacter, under long-term fertilization scenarios. There are a couple of explanations for this result. The genera Geobacter and Anaeromyxobacter, which belong to the class Deltaproteobacteria, are often classified as oligotrophic taxa [25]. Therefore, the abundance of this type of soil organism is expected to be inhibited in high-nutrient environments [26]. Moreover, taxa within this ecological cluster could have a lower capacity to downregulate N fixation, as they are outcompeted by other taxa under high N conditions but are highly competitive under low-nutrient conditions. Supporting this idea, some taxa within module #3, including Geobacter, have been reported to be particularly successful in fixing N in nutrient-poor environments [27].


Overall, our results suggest that long-term fertilization dramatically suppressed N fixation rates and the relative abundance of keystone and phylogenetically clustered N fixers. These findings have deepened our understanding on the linkage between N fixation and its associated N fixers under long-term fertilization scenarios. Moreover, our work provide solid evidence that the fundamental process of N fixation, and its associated microbial communities, will become more and more suppressed as terrestrial fertilization continues to increase.


Experimental design and sample collection

The experiment was set up in 1982 in Mengcheng county, Anhui province, China (33° 13′ N, 116° 35′ E, 42 m elevation) with typical lime concretion black soil. The annual temperature is 14.8 °C and the annual precipitation is 872 mm. Five fertilization treatments with wheat-soybean crop rotation were compared in a completely randomized block design with four replicates (each plot is 70 m2) [28, 29]: (1) control, non-fertilization; (2) NPK, NPK chemical fertilizers comprising urea (180 kg N ha−1 year−1), superphosphate (90 kg P2O5 ha−1 year−1) and potassium chloride (86 kg K2O ha− 1 y−1); (3) NPK + WS, NPK chemical fertilizers plus 7500 kg wheat straw ha−1 year−1; (4) NPK + PM, NPK chemical fertilizers plus 15,000 kg fresh pig manure ha−1 year−1; (5) NPK + CM, NPK chemical fertilizers plus 30,000 kg fresh cow manure ha−1 year−1. In the NPK + WS treatment, all the wheat straw were returned to the field, the pig manure in the NPK + PM treatment and the cow manure in the NPK + CM had the similar amount of organic carbon with the added wheat straw. Moreover, these contrasting types of fertilizers included in our fertilization treatments have different levels of availability for plants and microbes, e.g., from more labile (pig manure) to more recalcitrant (wheat straw and cow manure). We used a wide range of fertilization treatments aiming to make our results representative and applicable to contrasting management practices.

We dug around the wheat group (containing 30 to 40 wheat plants during the wheat filling stage on the 20th of April, 2017) to keep the root systems as intact as possible. The rhizosphere soil that was tightly adhered to the roots was then brushed. At the same time, the topsoil (0–15 cm deep) was collected as bulk soil using an auger corer (approximately 20 cm away from the plants). The collected soil was sieved through a 2 mm mesh to remove the impurities such as roots and stones. Some of the soil was stored at 4 °C for chemical analyses, and the rest was stored at − 40 °C for DNA extraction.

Soil physicochemical analysis

A pH meter (FE20 FiveEasy™, Mettler Toledo, Germany) was used to measure the soil pH at a soil to distilled water ratio of 1:5 (weight/volume). The soil moisture was determined gravimetrically by drying 5 g of fresh soil at approximately 105–108 °C to reach a constant weight and then calculating the weight ratio (evaporated water to dried soil). The total carbon (TC) and total nitrogen (TN) contents of the soil were determined by combustion of air-dried soil using a CNS-2000 analyzer (LECO, St. Joseph, MI, USA), after sieving the soil through a 0.15 mm mesh. The total phosphorus (TP) and total potassium (TK) contents of the soil were extracted after HF-HClO4 digestion and measured using the molybdenum blue method and flame spectrophotometry method (FP640, INASA, China), respectively. Dissolved organic carbon (DOC) was extracted by adding 50 mL distilled water to 5 g fresh soil, shaking for 1 h, and vacuum filtering through a G4 glass fiber filter with a pore space of 1.2 μm (Fisher), and then, the carbon contents in the extracts were determined by a total organic carbon analyzer (multi N/C 3000, Analytik Jena, Germany). Nitrate (NO3-N), ammonium (NH4+-N), and dissolved total nitrogen (DTN) were extracted at a ratio of 5 g fresh soil to 50 mL 2 M KCl. After shacking for 1 h, the extracts were filtered through a G4 glass fiber filter with a pore space of 1.2 μm (Fisher), and then, a continuous flow analytical system (San++ system, Skalar, Holland) was used to analyze the content of NO3-N, NH4+-N, and DTN. Dissolved organic nitrogen (DON) was calculated using the following formula: DON = DTN − NH4+-N − NO3-N. The available phosphorus (AP) in the soil was extracted by 0.5 M NaHCO3 and determined by using molybdenum blue method. Available potassium (AK) was extracted by 1 M ammonium acetate and determined by flame photometer (FP640, INASA, China) (Additional file 3: Appendix 1).

Determination of nitrogen fixation rates

The 15N2-labeling method is one of the most common and widely applied methods used for measuring N fixation rates [30, 31]. Five grams of soil were placed into 18 × 150 mm Balch tubes, and the headspace was replaced with synthetic air containing 80% 15N2 and 20% O2. The controls were filled with unlabeled N2 gas and processed in parallel. The tubes were incubated horizontally in the dark at room temperature for 22 days. The atom % 15N of soil samples was determined using a stable isotope ratio mass spectrometer (Flash 2000 HT/Conflo IV/Delta V, Thermo Fisher Scientific, Germany). Then, we calculated the net potential N fixation rate by comparing the difference of total 15N in soils receiving 15N2 relative to control.

High-throughput sequencing and bioinformatics analysis

For the DNA extraction, 0.5 g of fresh soil was used with the Fast DNA SPIN Kit (MP Biomedicals, Santa Ana, CA, USA). The nifH genes were amplified using primer pairs nifH-F/nifH-R (5′-AAAGGYGGWATCGGYAARTCCACCAC-3′)/(5′-TTGTTSGCSGCRTACATSGCCATCAT-3′) [32]. PCR reactions were performed in a 20 μL reaction containing 4 μL of 5 × FastPfu buffer, 2 μL of 2.5 mM dNTPs, 0.8 μL of 5 μM forward primer, 0.8 μL of 5 μM reverse primer, 0.4 μL of fastPfu Polymerase, 10 ng of template DNA, and double distilled water (ddH2O). Amplification was performed at 95 °C for 3 min, with 35 cycles of 95 °C for 30 s, 55 °C for 30 s, and 72 °C for 45 s, and extension at 72 °C for 10 min. PCR amplicons were purified by Agarose Gel DNA purification kit (TaKaRa Bio), and triplicate PCR amplifications for each sample were conducted and pooled as a PCR product and then sequenced on the platform of Illumina MiSeq PE300 (Majorbio Company in Shanghai, China). After sequencing, the nifH nucleotide sequences were analyzed using the QIIME-1.9.1 pipeline ( [33]. Firstly, the low-quality sequences (those with a quality score < 20, containing ambiguous nucleotides, or not matching the primer and barcode) were removed and the remaining sequences were further converted to amino acid sequences using the FunGene Pipeline of the Ribosomal Database Project [34]. The sequences encoding proteins that did not match the nifH protein sequence or that contained termination codons were discarded. The remaining sequences were aligned against the nifH gene database [35], removing both the failed and chimeric sequences. The remaining high-quality sequences were clustered into operational taxonomic units (OTUs) at 95% similarity with UCLUST [36] running in de novo mode, and all singleton OTUs were deleted.

Co-occurrence network analysis

We constructed a co-occurrence network with all the samples (rhizosphere and bulk soil) and identified the main ecological clusters of strongly associated OTUs. The top OTUs, accounting for more than 80% of the relative abundance in the total community, were chosen [37]. All pair-wised Spearman correlations between OTUs were calculated, and the correlations with a Spearman’s coefficient of less than 0.65 and a P value of more than 0.01, were removed. This allowed us to focus only on the OTUs that strongly co-occurred and were more likely to interact with each other. The main modules (ecological clusters) in the network were visualized using Gephi ( The relative abundance of each ecological cluster was calculated by averaging the standardized relative abundances (z-score) of the species that belonged to it (Additional file 3: Appendix 3).

Statistical analysis

ANOVA and pairwise t tests were used to compare the soil variables, dominant microbial taxa, and the microbial alpha diversity between different fertilization treatments (Additional file 3: Appendix 1). These tests were implemented using SPSS 21. Mantel test was used to analyze the correlations between the diazotrophic community and physicochemical properties (Additional file 3: Appendix 2). This was performed using the “vegan” package in R × 32 (3.2.2). A principal coordinate analysis (PCoA) was used to find significant differences in diazotrophic communities between sampling groups (Additional file 3: Appendix 2). The PCoA was carried out using the “labdsv” package R × 32 (3.2.2) (

Phylogenetic analyses

The nifH gene provides sufficient phylogenetic resolution [38] in ecological studies. The phylogenetic tree for the 481 dominant diazotrophic phylotypes in the ecological clusters was built using FastTree [39], and visualized using GraPhlAn [40]. Phylogenetic sampling theory can be analytically employed (assuming the random sampling from the phylogenetic tree as the predicted phylogenetic diversity in a local community and then comparing the observed phylogenetic diversity with those predictions) [20] to determine the degree to which diazotrophic community appear random (between − 2 and + 2), over-dispersed (above + 2), or clustered (below − 2). Phylogenetic sampling theory was performed using the R package “picante” [41]. An advantage of randomly sampling the regional phylogenetic tree is that it can be used to compare samples of unequal sizes based on the binomial sampling model [42]. The differences between observed and expected phylogenetic diversity were determined by calculating and comparing z-scores for each ecological cluster. When the observed phylogenetic diversity is less than the expected diversity (below − 2), the microbial community in the ecological cluster is considered to be phylogenetically clustered, which means that closely related taxa are more likely to be sampled [43] and actively selected by the environment [42].

Structural equation modeling analysis

The SEM [44] was conducted using IBM SPSS Amos 21 (Chicago, IL: Amos Development Corporation). It was used to evaluate the direct and indirect effects of the soil physicochemical properties and the relative abundance of the main ecological clusters on the N fixation rates. The physiochemical properties of the soil included soil pH, total carbon, total nitrogen, and total phosphorus. In the model, the treatments (control, NPK, NPK + WS, NPK + PM, and NPK + CM) were categorical variables with two levels: 1 (a particular treatment) and 0 (remaining considered treatments). In addition, bootstrapping was used to test the probability that path coefficients differed from zero, as a few of the variables were not normally distributed. We also calculated the standardized total effects (STEs) of the soil properties, fertilization treatments, and rhizosphere effect on the N fixation rate to aid interpretation of the SEM.

Random Forest modeling analysis

Random Forest regression (R package “randomForest”) was used to regress the normalized OTUs in different treatments. The 10-fold cross-validation method was used to determine the optimal set of OTUs correlated to the N fixation rates [45]. Ranked lists of OTUs in order of Random Forests reported feature importance scores were achieved based on the increase in mean-square error of nitrogen fixation rates predicted over 100 iterations of the algorithm. The 50 marker OTUs were chosen based on the minimum average cross-validation mean-squared errors, which were obtained from five trials of the 10-fold cross-validation.

Availability of data and materials

The obtained sequences were submitted to the NCBI Sequence Read Archive (SRA) with accession number SRP149667 ( Other data and result supporting the findings of the study are available in this article and its supplementary information files.


  1. 1.

    Vitousek PM, Menge DNL, Reed SC, Cleveland CC. Biological nitrogen fixation: rates, patterns and ecological controls in terrestrial ecosystems. Philos Trans R Soc B Biol Sci. 2013;368:20130119.

  2. 2.

    Nadelhoffer KJ. Chapter 14 - the impacts of nitrogen deposition on forest ecosystems. In: Hatfield JL, Follett RF, editors. Nitrogen in the Environment (Second Edition). Edited by. San Diego: Academic Press; 2008. p. 463–82.

  3. 3.

    Burns RC, Hardy RW. Nitrogen fixation in bacteria and higher plants. Mol Biol Biochem Biophys. 1975;21:1-189.

  4. 4.

    Wang C, Zheng M, Song W, Wen S, Wang B, Zhu C, Shen R. Impact of 25 years of inorganic fertilization on diazotrophic abundance and community structure in an acidic soil in southern China. Soil Biol Biochem. 2017;113:240–9.

  5. 5.

    Rockström J, Steffen W, Noone K, Persson Å, Chapin FS III, Lambin E, Lenton TM, Scheffer M, Folke C, Schellnhuber HJ. Planetary boundaries: exploring the safe operating space for humanity. Ecol Soc. 2009;14(2):32.

  6. 6.

    van Kessel C, Hartley C. Agricultural management of grain legumes: has it led to an increase in nitrogen fixation? Field Crops Res. 2000;65:165–81.

  7. 7.

    Wang J, Zhang D, Zhang L, Li J, Raza W, Huang Q, Shen Q. Temporal variation of diazotrophic community abundance and structure in surface and subsoil under four fertilization regimes during a wheat growing season. Agric Ecosyst Environ. 2016;216:116–24.

  8. 8.

    Chase Jonathan M, Myers JA. Disentangling the importance of ecological niches from stochastic processes across scales. Philos Trans R Soc B Biol Sci. 2011;366:2351–63.

  9. 9.

    Feng M, Adams JM, Fan K, Shi Y, Sun R, Wang D, Guo X, Chu H. Long-term fertilization influences community assembly processes of soil diazotrophs. Soil Biol Biochem. 2018;126:151–8.

  10. 10.

    Fierer N, Bradford MA, Jackson RB. Toward an ecological classification of soil bacteria. Ecology. 2007;88:1354–64.

  11. 11.

    Postgate JR. New advances and future potential in biological nitrogen fixation. J Appl Microbiol. 1974;37:185–202.

  12. 12.

    Mancinelli RL, McKay CP. The evolution of nitrogen cycling. Orig Life Evol Biosph. 1988;18:311–25.

  13. 13.

    Vitti JJ, Grossman SR, Sabeti PC. Detecting natural selection in genomic data. Annu Rev Genet. 2013;47:97–120.

  14. 14.

    Jerne NK. The natural-selection theory of antibody formation. Proc Natl Acad Sci U S A. 1955;41:849–57.

  15. 15.

    Stinchcombe JR, Rutter MT, Burdick DS, Tiffin P, Rausher MD, Mauricio R. Testing for environmentally induced bias in phenotypic estimates of natural selection: theory and practice. Am Nat. 2002;160:511–23.

  16. 16.

    Duncan N, Menge L, Levin SA, Hedin LO. Facultative versus obligate nitrogen fixation strategies and their ecosystem consequences. Am Nat. 2009;174:465–77.

  17. 17.

    Sheffer E, Batterman SA, Levin SA, Hedin LO. Biome-scale nitrogen fixation strategies selected by climatic constraints on nitrogen cycle. Nat Plants. 2015;1:15182.

  18. 18.

    Collavino MM, Tripp HJ, Frank IE, Vidoz ML, Calderoli PA, Donato M, Zehr JP, Aguilar OM. nifH pyrosequencing reveals the potential for location-specific soil chemistry to influence N 2-fixing community dynamics. Environ Microbiol. 2014;16:3211–23.

  19. 19.

    Frank IE, Turk-Kubo KA, Zehr JP. Rapid annotation of nifH gene sequences using classification and regression trees facilitates environmental functional gene analysis. Environ Microbiol Rep. 2016;8:905–16.

  20. 20.

    O'Dwyer JP, Kembel SW, Green JL. Phylogenetic diversity theory sheds light on the structure of microbial communities. PLoS Comput Biol. 2012;8:e1002832.

  21. 21.

    Lovley DR, Ueki T, Zhang T, Malvankar NS, Shrestha PM, Flanagan KA, Aklujkar M, Butler JE, Giloteaux L, Rotaru A-E, et al. Geobacter: the microbe electric’s physiology, ecology, and practical applications. Adv Microb Physiol. 2011;59:1–100.

  22. 22.

    Coenye T, Vandamme P. Diversity and significance of Burkholderia species occupying diverse ecological niches. Environ Microbiol. 2003;5:719–29.

  23. 23.

    Rivas R, Martens M, De Lajudie P, Willems A. Multilocus sequence analysis of the genus Bradyrhizobium. Syst Appl Microbiol. 2009;32:101–10.

  24. 24.

    Weyens N, van der Lelie D, Taghavi S, Newman L, Vangronsveld J. Exploiting plant–microbe partnerships to improve biomass production and remediation. Trends Biotechnol. 2009;27:591–8.

  25. 25.

    Trivedi P, Anderson IC, Singh BK. Microbial modulators of soil carbon storage: integrating genomic and metabolic knowledge for global prediction. Trends Microbiol. 2013;21:641–51.

  26. 26.

    Koch AL. Oligotrophs versus copiotrophs. Bioessays. 2001;23:657–61.

  27. 27.

    Bazylinski DA, Dean AJ, Schüler D, Phillips EJ, Lovley DR. N2-dependent growth and nitrogenase activity in the metal-metabolizing bacteria, Geobacter and Magnetospirillum species. Environ Microbiol. 2000;2:266–73.

  28. 28.

    Sun R, Zhang X-X, Guo X, Wang D, Chu H. Bacterial diversity in soils subjected to long-term chemical fertilization can be more stably maintained with the addition of livestock manure than wheat straw. Soil Biol and Biochem. 2015;88:9–18.

  29. 29.

    Sun R, Dsouza M, Gilbert JA, Guo X, Wang D, Guo Z, Ni Y, Chu H. Fungal community composition in soils subjected to long-term chemical fertilization is most influenced by the type of organic matter. Environ Microbiol. 2016;18:5137–50.

  30. 30.

    Hsu S-F, Buckley DH. Evidence for the functional significance of diazotroph community structure in soil. ISME J. 2009;3:124.

  31. 31.

    Montoya JP, Voss M, Kahler P, Capone DG. A simple, high-precision, high-sensitivity tracer assay for N2 fixation. Appl Environ Microbiol. 1996;62:986–93.

  32. 32.

    Rösch C, Mergel A, Bothe H. Biodiversity of denitrifying and dinitrogen-fixing bacteria in an acid forest soil. Appl Environ Microbiol. 2002;68:3818–29.

  33. 33.

    Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.

  34. 34.

    Wang Q, Quensen JF, Fish JA, Lee TK, Sun Y, Tiedje JM, Cole JR. Ecological patterns of nifH genes in four terrestrial climatic zones explored with targeted metagenomics using FrameBot, a new informatics tool. MBio. 2013;4:e00592–13.

  35. 35.

    Gaby JC, Buckley DH. A comprehensive aligned nifH gene database: a multipurpose tool for studies of nitrogen-fixing bacteria. Database. 2014;2014:bau001.

  36. 36.

    Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.

  37. 37.

    Delgado-Baquerizo M, Reith F, Dennis PG, Hamonts K, Powell JR, Young A, Singh BK, Bissett A. Ecological drivers of soil microbial diversity and soil biological networks in the Southern Hemisphere. Ecology. 2018;99:583–96.

  38. 38.

    Pinto A, Evangelista T, Carvalho M, Alves M, Luis MS. Respiratory assistance with a non-invasive ventilator (Bipap) in MND/ALS patients: survival rates in a controlled trial. J Neurol Sci. 1995;129:19–26.

  39. 39.

    Price MN, Dehal PS, Arkin AP. FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5:e9490.

  40. 40.

    Asnicar F, Weingart G, Tickle TL, Huttenhower C, Segata N. Compact graphical representation of phylogenetic data and metadata with GraPhlAn. PeerJ. 2015;3:e1029.

  41. 41.

    Kembel SW, Cowan PD, Helmus MR, Cornwell WK, Morlon H, Ackerly DD, Blomberg SP, Webb CO. Picante: R tools for integrating phylogenies and ecology. Bioinformatics. 2010;26:1463–4.

  42. 42.

    Burns AR, Stephens WZ, Stagaman K, Wong S, Rawls JF, Guillemin K, Bohannan BJ. Contribution of neutral processes to the assembly of gut microbial communities in the zebrafish over host development. ISME J. 2016;10:655–64.

  43. 43.

    Horner-Devine MC, Bohannan BJ. Phylogenetic clustering and overdispersion in bacterial communities. Ecology. 2006;87:S100–8.

  44. 44.

    Grace JB. Structural equation modeling and natural systems. Cambridge: Cambridge University Press; 2006.

  45. 45.

    Subramanian S, Huq S, Yatsunenko T, Haque R, Mahfuz M, Alam MA, Benezra A, DeStefano J, Meier MF, Muegge BD, et al. Persistent gut microbiota immaturity in malnourished Bangladeshi children. Nature. 2014;510:417.

Download references


We thank Zhibin Guo, Keke Hua, Yingying Ni, Yuying Ma, Hongfei Wang, Liang Chen, Yu Shi, Yuntao Li, Ruibo Sun, and Xiangtian Meng for their assistance in field management, soil sampling, and lab analyses.


This work was funded by the Natural Science Foundation of China (31870480), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB15010101), and the China Biodiversity Observation Networks (Sino BON). M.D-B. also acknowledges the support from the Marie Sklodowska-Curie Actions of the Horizon 2020 Framework Program H2020-MSCA-IF-2016 under REA grant agreement no. 702057.

Author information

HC, MDB, and YZ designed research. KF, XG, DW, YW, WY, MZ, and HY performed the research. KF, MDB, and HC analyzed data. KF, MDB, YZ, and HC wrote the paper. All authors read and approved the final manuscript.

Correspondence to Haiyan Chu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1. Physicochemical soil properties in bulk soil and rhizosphere soil among different fertilization treatments. Table S2. The relative abundance of dominant diazotrophic genera among different treatments. Table S3. Pairwise t-tests for the relative abundance of dominant diazotrophic genera between bulk soil and rhizosphere soil in different fertilization treatments. Table S4. Comparison of diazotrophic alpha-diversity indexes among different fertilization treatments. Table S5. Pairwise t-testing of diazotrophic alpha-diversity in bulk soil and rhizosphere soil. Table S6. Pairwise t-tests of nitrogen fixation rates in bulk soil and rhizosphere soil in different fertilization treatments. Table S7. ADONIS double factor analysis for the diazotrophic community. Table S8. Pairwise t-tests of relative abundance of main ecological clusters within the diazotrophic community in bulk soil and rhizosphere soil under different fertilization treatments. Table S9. The network properties for the main ecological clusters of the diazotrophic community. Table S10. Operational taxonomic unit (OTU) count properties of important species for nitrogen fixation rates found by the Random Forest model in the three main ecological clusters. Table S11. The correlations (r) and significance (p) were determined using a Mantel test between the diazotrophic community and environmental variables in bulk soil and rhizosphere soil. Table S12. Spearman correlation between physicochemical soil properties and diazotrophic alpha-diversity. Table S13. Spearman correlation between physicochemical soil properties and relative abundance of the main diazotrophic ecological clusters.

Additional file 2: Figure S1. Relative abundance of the dominant diazotrophic genera in different fertilization treatments. Figure S2. A random forest model was applied to regress the diazotrophic OTU profiling in bulk soil and rhizosphere soil against the nitrogen fixation rates. Figure S3. Correlations between the relative abundance of important species for nitrogen fixation rates found by the Random Forest model and their Importance Index in different fertilization treatments. Figure S4. Diazotrophic community variations in different fertilization samples; and diazotrophic community composition variations which were based on Bray-Curtis distances by principal coordinate analysis.

Additional file 3. Supplementary Results. Appendix 1. Soil properties and diazotrophic community under long-term fertilization scenarios. Appendix 2. Edaphic factors associated with the soil diazotrophic community under long-term fertilization scenarios. Appendix 3. Diazotrophic ecological clusters and associated edaphic factors.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Fan, K., Delgado-Baquerizo, M., Guo, X. et al. Suppressed N fixation and diazotrophs after four decades of fertilization. Microbiome 7, 143 (2019).

Download citation


  • Diazotrophs
  • Nitrogen fixation rates
  • Ecological clusters
  • Long-term fertilization