Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSphere
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Applied and Environmental Science

Soil Health Management Enhances Microbial Nitrogen Cycling Capacity and Activity

Jialin Hu, Virginia L. Jin, Julie Y. M. Konkel, Sean M. Schaeffer, Liesel G. Schneider, Jennifer M. DeBruyn
Katherine McMahon, Editor
Jialin Hu
aDepartment of Biosystems Engineering & Soil Science, University of Tennessee, Knoxville, Tennessee, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Virginia L. Jin
bUSDA—Agricultural Research Service, Agroecosystem Management Research Unit, University of Nebraska-Lincoln, Lincoln, Nebraska, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Julie Y. M. Konkel
aDepartment of Biosystems Engineering & Soil Science, University of Tennessee, Knoxville, Tennessee, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sean M. Schaeffer
aDepartment of Biosystems Engineering & Soil Science, University of Tennessee, Knoxville, Tennessee, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Liesel G. Schneider
cDepartment of Animal Science, University of Tennessee, Knoxville, Tennessee, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jennifer M. DeBruyn
aDepartment of Biosystems Engineering & Soil Science, University of Tennessee, Knoxville, Tennessee, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jennifer M. DeBruyn
Katherine McMahon
University of Wisconsin-Madison
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSphere.01237-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Soil microbial transformations of nitrogen (N) can be affected by soil health management practices. Here, we report in situ seasonal dynamics of the population size (gene copy abundances) and functional activity (transcript copy abundances) of five bacterial genes involved in soil N cycling (ammonia-oxidizing bacteria [AOB] amoA, nifH, nirK, nirS, and nosZ) in a long-term continuous cotton production system under different management practices (cover crops, tillage, and inorganic N fertilization). Hairy vetch (Vicia villosa Roth), a leguminous cover crop, most effectively promoted the expression of N cycle genes, which persisted after cover crop termination throughout the growing season. Moreover, we observed similarly high or even higher N cycle gene transcript abundances under vetch with no fertilizer as no cover crop with N fertilization throughout the cover crop peak and cotton growing seasons (April, May, and October). Further, both the gene and transcript abundances of amoA and nosZ were positively correlated to soil nitrous oxide (N2O) emissions. We also found that the abundances of amoA genes and transcripts both positively correlated to field and incubated net nitrification rates. Together, our results revealed relationships between microbial functional capacity and activity and in situ soil N transformations under different agricultural seasons and soil management practices.

IMPORTANCE Conservation agriculture practices that promote soil health have distinct and lasting effects on microbial populations involved with soil nitrogen (N) cycling. In particular, using a leguminous winter cover crop (hairy vetch) promoted the expression of key functional genes involved in soil N cycling, equaling or exceeding the effects of inorganic N fertilizer. Hairy vetch also left a legacy on soil nutrient capacity by promoting the continued activity of N cycling microbes after cover crop termination and into the main growing season. By examining both genes and transcripts involved in soil N cycling, we showed different responses of functional capacity (i.e., gene abundances) and functional activity (i.e., transcript abundances) to agricultural seasons and management practices, adding to our understanding of the effects of soil health management practices on microbial ecology.

INTRODUCTION

Nitrogen (N) is a major nutrient for plant growth, and its transformation and availability in soil is a key factor in terrestrial ecosystem productivity (1). Nitrogen fertilizer has been widely used to improve crop yields in agroecosystems, but losses of excess N result in the eutrophication of waterbodies and exacerbate atmospheric concentrations of the greenhouse gas nitrous oxide (N2O) (2). Conservation agriculture and soil health practices, such as growing cover crops, reducing N fertilization, and minimizing tillage, ensure mineral nutrients can be properly used by crops and do not become surplus or deficient (3).

The C-to-N ratios of plant residues can drive the net balance of soil mineral N by affecting microbial mineralization and immobilization of N (4). Cover crops, in addition to the main crop, can take up soil nitrate (NO3−) to reduce N losses via leaching and denitrification (5). Further, leguminous cover crops can supplement crop N requirements through biological N fixation, thereby reducing exogenous N fertilizer requirements and associated N losses (6). The effects of tillage intensity on N losses through N2O emissions are cropping system specific, but reduced or no tillage generally mitigates N2O emissions compared to conventional tillage, especially in long-term (>10 years) production systems (7, 8).

Conservation agriculture practices can change soil properties, which are important factors influencing microbial functional groups involved in soil N cycling (9). For example, increased soil organic matter (SOM) and soil water content (SWC) under reduced tillage and winter cover cropping can create nutrient-rich anaerobic microsites that favor microbial denitrification (10, 11). N limitation caused by reduced N fertilization may decrease the population size of ammonia-oxidizing bacteria (AOB) (12) but promote the activity of diazotrophs for N fixation (13).

The abundance and expression of genes encoding key enzymes involved in soil N transformations have been widely used to evaluate the abundances and activity of microbial populations involved in specific N cycling processes (14). For example, different microbial assemblages are responsible for N fixation, nitrification, and denitrification (15). Biological N fixation is the reduction of atmospheric N2 to biologically available ammonium. Molecular analysis of N-fixing bacteria often targets nifH genes encoding a subunit of nitrogenase reductase, which catalyzes the reduction of dinitrogen to ammonia (16–18). Nitrification is the process of converting ammonia to nitrate and is one of the major processes contributing to N2O production (19–22) and nitrate leaching (23) in soils. The first and rate-limiting step of nitrification is microbial oxidation of ammonia to hydroxylamine catalyzed by ammonia monooxygenases (AMO). Thus, molecular analysis of ammonia-oxidizing bacteria (AOB) and archaea (AOA) commonly targets amoA genes encoding subunit A of AMO (24, 25). Although AOA may be more numerically abundant than AOB in soil (26), AOB have been shown to be generally more responsive than AOA to soil management practices, such as N fertilization (27–29) and tillage (30), and were shown to be functionally dominant for ammonia oxidation in agricultural topsoils (31). Denitrification is the stepwise reduction of nitrate to nitrite, nitric oxide, nitrous oxide, and dinitrogen gas under anaerobic conditions and is catalyzed by nitrate, nitrite, nitric oxide, and nitrous oxide reductases, respectively (32). Reduction of nitrite to nitric oxide is catalyzed by two structurally different nitrite reductases (NIR) carried by different taxa as follows: (i) Cu-containing NIR (CuNIR) encoded by a single gene nirK, sometimes accompanied by gene nirV, and (ii) cytochrome cd1 NIR (cdNIR) encoded by at least 10 genes (nirSECFDLGHJN) with nirS gene encoding the functional subunits of cdNIR (33). The reduction of N2O to N2 gas is catalyzed by nitrous oxide reductase (NOS), which is frequently studied by using the marker gene nosZ that encodes the multi-Cu NOS catalytic subunit (34, 35). The variation and the ratio between amoA, nirK, and nirS that are involved in N2O production and nosZ that is involved in N2O reduction have been used to predict system-level N2O emissions (36, 37).

The reported effects of conservation soil management practices on important functional groups involved in N cycling are inconsistent. For example, many studies indicate that reduced or no tillage can increase the abundance of both AOA and AOB amoA (38), nirK (39), nosZ (39, 40), and (nirK + nirS)/nosZ ratios (41). Other studies, however, reported a positive effect of no tillage on AOB amoA (30) but did not find significant effects of tillage on the abundances of AOA amoA (30), nirK (42), or nirS genes (39). Cover cropping may also affect some populations: nirK gene abundances were higher under cereal rye compared to hairy vetch while nirS gene abundances were unaffected by cover crop species (43). Nitrogen input can suppress nifH abundance when the available N is excessive (44, 45), but it can also improve abundances when available N is limited (46). These inconsistent and sometime contradictory results may be due to a lack of temporal resolution. For example, field trials may not have been conducted long enough to invoke a significant change in the soil system. Moreover, there is seasonal variability in microbial functional dynamics driving N transformations that are not captured if only a single time point sample is taken. Another critical gap in our current understanding of functional group dynamics is the relationship between functional capacity and activity. Many studies have only focused on N cycling functional groups at either the genetic level (i.e., using extracted DNA), which misses changes in activity or expression, or at the transcriptional level (i.e., using extracted RNA), which misses changes in population abundances.

The objectives of our study were to (i) determine the effects of tillage, cover cropping, and N fertilization on abundance and activity of N cycling microbes; (ii) investigate the seasonal dynamics of abundance and activity of N cycling functional groups under long-term conservation soil management; and (iii) identify relationships among the N cycle functional gene and transcript abundances, microbial-controlled N transformation rates, and soil physicochemical parameters. The study was conducted at a long-term (36 years) continuous cotton experiment with different tillage, cover cropping, and N fertilization regimes. To examine the dynamics of soil N cycling functional groups, we used quantitative PCR (qPCR) and real-time quantitative reverse transcription-PCR (qRT-PCR) to target the N cycle genes and transcripts of nifH, AOB amoA, nirK, nirS, and nosZ. We then explored the seasonal dynamics of soil N cycling populations under these different conservation agricultural management practices and their relationships to soil N pools and fluxes.

RESULTS

Soil properties and ambient temperatures.All of the soil physicochemical parameters and N transformation rates measured in this study showed seasonal dynamics, and most of them were significantly affected by hairy vetch (see Table S2 and S3 in the supplemental material). In particular, increased nitrification and N mineralization rates in vetch plots were mainly observed in cover crop peak season (April) and spring crop transition (May) (Table S3). N fertilization significantly increased the concentration of NO3−-N, total soil C (TC), total soil N (TN), N2O-N, and the rate of incubated nitrification and N mineralization (Table S3). Compared to conventional tillage (CT), no tillage (NT) significantly increased the concentration of total extractable C (TEC), TC, and TN, which was particularly noticeable in April and October for TEC and in May and October for TC and TN (Table S3). Ambient temperatures during the cotton growing seasons (May and October, 16.1 and 17.2°C, respectively) were higher than during the cover crop growing seasons (November and April, 8.3 and 9.4°C, respectively).

N cycle genes and transcripts.The abundances of N cycle genes and transcripts were influenced differently by agricultural season and soil management practices. In general, for both relative (normalized to 16S rRNA genes and 16S rRNA) and absolute abundances (per gram dry weight soil), genes were mostly affected by only main effects or two-way interaction effects of agricultural season and soil management practices, while transcripts were affected by much more complex interaction effects among season and management practices (Table 1; see also Table S4 and S5 in the supplemental material).

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 1

Results of mixed model ANOVA testing effects of agricultural season and soil management practices on the normalized and absolute abundances of N cycle genes and transcriptsa

Both relative and absolute abundances of marker genes for N fixation (nifH) were significantly affected by agricultural season and tillage (Table 1), with highest abundances during cover crop harvest (April) and under CT (Fig. 1). Moreover, although no significant effect of N fertilization on the relative abundance of nifH genes was observed (Fig. 2A), absolute abundance of nifH was significantly higher under no N fertilization than N fertilization, especially in no cover (NC) and vetch (V) cropping systems (Fig. 2B). Expression of nifH had different and more complex dynamics. N fertilization (67 kg N ha−1 [67N]) increased nifH transcript relative abundance in NT but decreased it in CT (see Fig. S2A in the supplemental material). In general, nifH transcript abundances were higher during fall crop transition (November) than other seasons (Fig. 3 and 4). nifH transcript abundances were also significantly affected by the interaction of season, cover crop, and tillage (Table 1). The relative abundance of nifH transcripts was not significantly different regardless of cover crops or tillage in April but was significantly higher under CT than NT in spring crop transition (May) for NC plots, in cotton peak season (October) for NC and V plots, and in fall crop transition (November) for wheat (W) plots (Fig. 3A). The absolute abundance of nifH transcript was significantly higher under CT than NT in April and May for NC plots, in October for V plots, and in November for W plots (Fig. 3B). The interaction of season, cover crop, and N fertilization also significantly affected nifH transcripts (Table 1). Fertilizer generally had a positive effect on nifH transcripts in NC plots but a negative effect in W plots, especially in October (Fig. 4).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Seasonal dynamics of 16S rRNA gene normalized nifH, amoA, nirK, nirS, and nosZ gene abundances (A) and absolute abundances of nifH, amoA, nirK, nirS, nosZ, and 16S rRNA genes (B) in relation to tillage. Points represent the mean ± standard error (n = 24). NT, no tillage; CT, conventional tillage.

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Variation of 16S rRNA gene normalized nifH, amoA, nirK, nirS, and nosZ gene abundances (A) and absolute abundances of nifH, amoA, nirK, nirS, nosZ, and 16S rRNA genes (B) in relation to N fertilization rate (0N, no N fertilization; 67N, 67 kg N ha−1 fertilization) under different cover crop treatments (NC, no cover; V, vetch; W, wheat). Points represent the mean ± standard error (n = 32).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Seasonal dynamics of 16S rRNA normalized transcript abundances of nifH, amoA, nirK, nirS, and nosZ (A) and absolute transcript abundances of nifH, amoA, nirK, nirS, nosZ, and 16S rRNA gene (B) in relation to cover crops and tillage. Points represent the mean ± standard error (n = 8). NC, no cover; V, vetch; W, wheat; NT, no tillage; CT, conventional tillage.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Seasonal dynamics of 16S rRNA normalized transcript abundances of nifH, amoA, nirK, nirS, and nosZ (A) and absolute transcript abundances of nifH, amoA, nirK, nirS, nosZ, and 16S rRNA gene (B) in relation to cover crop treatments and N fertilization rate. Points represent the mean ± standard error (n = 8). NC, no cover; V, vetch; W, wheat; 0N, no N fertilization; 67N, 67 kg N ha−1 fertilization.

Both relative and absolute abundances of ammonia oxidation gene AOB amoA were significantly affected by season combined with tillage (Table 1), with the highest relative abundance in October under both CT and NT and lowest in November under CT (Fig. 1A). The highest absolute amoA abundance was in May under CT and in October under both CT and NT and lowest in November under both CT and NT (Fig. 1B). The amoA gene abundances were also significantly affected by the combination of cover crop and N fertilization (Table 1): N fertilization increased amoA gene relative abundance under V and NC but decreased both relative and absolute abundances under wheat (Fig. 2). Similar to nifH, addition of N fertilization increased relative abundance of amoA transcripts under NT but decreased under CT (Fig. S2A). Their absolute abundance was not affected by the interaction of N fertilization and tillage (Fig. S2B). In addition, amoA transcript abundances were significantly affected by the combination of season, cover crop, and tillage (Table 1). In April and May, amoA transcript abundances were highest under V in both NT and CT (Fig. 3). After a lag in April and May, NC plots had increased amoA transcripts in October under CT and in November under NT (Fig. 3). The interaction of season, cover crop, and N fertilization also significantly affected amoA transcripts (Table 1). Addition of fertilizer significantly increased amoA transcript relative abundance under W in April and under NC in May, though not to the same level as V (Fig. 4A). In October and November, 67N continued to exhibit increased relative abundance of amoA transcripts under NC; in contrast, 67N caused a decrease under W in October and November and V in November (Fig. 4A). The absolute abundance of amoA transcripts showed similar trends to their relative abundance (Fig. 4B).

The abundances of nitrite reduction gene nirK varied by agricultural season (Table 1), with lowest relative and absolute abundances in May and April, respectively (Fig. 1). CT increased absolute abundance rather than relative abundance of nirK genes (Fig. 1). N fertilization significantly increased their absolute abundance under V (Fig. 2B). The relative abundance of nirK transcripts increased with fertilization in NT plots but decreased under CT (Fig. S2A). In addition, the nirK transcript abundances were significantly affected by the combination of season, cover crop, and tillage (Table 1). In April, May, and October, nirK transcript abundances were generally highest under V, which was more obvious for their absolute abundance (Fig. 3). In November, CT plots had similar levels of nirK transcripts regardless of cover, while NT plots had much more variable abundances between cover crops (Fig. 3). The interaction of season, cover crop, and N fertilization rate also significantly affected both relative and absolute abundances of nirK transcripts (Table 1). Fertilization generally increased nirK transcripts under W in April and November but decreased in May and October (Fig. 4).

The abundances of nitrite reduction gene nirS were significantly affected by agricultural season and tillage (Table 1), with higher abundances in November and under CT (Fig. 1). nirS gene abundances were also significantly affected by the combination of cover crop and N fertilization rate (Table 1): N fertilization increased nirS gene relative abundance in NC and W plots but had no effect in V plots (Fig. 2A); N fertilization increased nirS gene absolute abundance in W plots but decreased in V plots (Fig. 2B). nirS transcript abundances were significantly affected by the interaction between agricultural season and all three soil management practices (Table 1). The highest transcript relative abundance of nirS was observed in NT-V-0 kg N ha−1 (0N), NT-V-67N, CT-V-0N, and NT-NC-67N in April, May, October, and November, respectively, while the highest absolute abundance was observed in NT-V-0N, CT-V-0N, CT-V-67N, and NT-NC-0N in April, May, October, and November, respectively (Table S4 and S5).

The abundances of nitrous oxide reduction gene nosZ were significantly affected by agricultural season (Table 1), with both relative and absolute abundances highest in October (Fig. 1). Their absolute abundance was higher in CT than NT (P < 0.0001) (Fig. 1B). nosZ gene relative abundance was highest under vetch (P = 0.0024) and 67N treatments (P < 0.0001), while absolute abundance was only increased by 67N under wheat (Fig. 2). As with nirS transcripts, there was a significant interaction between agricultural season and the three soil management practices on nosZ transcript abundances (Table 1). nosZ transcript relative abundance was highest in NT-NC-67N in April and November and CT-V-0N in May and October and lowest in NT-W-0N in April and November, NT-NC-0N in May, and CT-W-67N in October (Table S4). nosZ transcript absolute abundance was highest in CT-V-67N in April and October, CT-V-0N in May, and NT-V-67N in November and lowest in NT-NC-0N in April and May, CT-W-67N in October, and CT-W-0N in November (Table S5).

16S rRNA genes and 16S rRNA.The abundance of 16S rRNA genes showed significant seasonal dynamics (P < 0.0001). Highest abundances were observed in May and lowest in April (Table S5). Compared to NT, CT slightly but significantly promoted 16S rRNA gene abundances from a mean of 2.35 × 109 copies g−1 dry weight soil in NT to 3.02 × 109 copies g−1 dry weight soil in CT (Table S5). 16S rRNA gene copy abundances decreased from 2.92 × 109 copies g−1 dry weight soil in 0N to 2.43 × 109 copies g−1 dry weight soil in 67N (Table S5). 16S rRNA copy numbers were significantly affected by the combination of season and cover crop (P < 0.05), which was lowest in April under NC (9.01 × 1010 copies g−1 dry weight soil) and highest in October in V (5.74 × 1011 copies g−1 dry weight soil; see Table S5). Moreover, 67N increased 16S rRNA abundances in CT but decreased them in NT (Fig. S2). Unlike N cycle gene transcripts, the abundance of 16S rRNA was not significantly affected by the interaction of season, tillage, and cover crop/N fertilization (Table 1; Fig. 3 and 4).

Cover crop effects on N cycling microbial activity.In general, V plots tended to exhibit the highest N cycle gene transcript abundances across seasons. The positive effect of 67N addition on N cycle gene transcripts was primarily observed in plots with no cover crops (NC), especially during cotton peak season (October) (Fig. 4). It was notable that fertilization did not have a positive effect in V plots, with V-0N exhibiting similarly high or even higher N cycle gene transcript abundances as V-67N, especially during the cotton growing season (May and October) (Fig. 4). We also observed similarly high or even higher N cycle gene transcript abundances under vetch with no fertilizer (V-0N) as no cover crop with N fertilization (NC-67N) throughout the cover crop harvest and cotton growing seasons (April, May, and October) (Fig. 4).

The transcript/gene ratios of N cycle genes also generally remained higher in V plots under the treatment combination of cover crop with tillage and with N fertilization in April, May, and October (Fig. 5; see Table S6 in the supplemental material), reconfirming the enhancing effect of vetch on N cycle gene expressions. Moreover, the seasonal dynamics of transcript/gene ratio in relation to cover crop with tillage/N fertilization was generally similar to the trends in transcript absolute abundances (Fig. 5), indicating that season and soil management practices had stronger effects at the transcript level compared to the gene abundance level.

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Seasonal dynamics of transcript/gene ratio of nifH, amoA, nirK, nirS, and nosZ (B) in relation to cover crops and tillage (A) and in relation to cover crops and N fertilization (B). Points represent the mean ± standard error (n = 8). NC, no cover; V, vetch; W, wheat; NT, no tillage; CT, conventional tillage.

Variability of N cycling microbial functional potential and activity.Nonmetric multidimensional scaling (NMDS) analysis showed the variability of the nitrogen cycling bacterial populations in terms of functional potential (i.e., based on the five gene relative abundances) and functional activity (i.e., based on the five transcript relative abundances) (Fig. 6). As expected, the statistical dispersion, or variability, of the populations by activity was consistently and significantly higher than for functional potential (dispersion index = 0.3926 for transcripts and 0.2168 for genes; P < 0.001). Analysis of similarities (ANOSIM) testing showed overall significant differences in functional potential by season (R = 0.1748; P = 0.001), tillage (R = 0.0222; P = 0.020), and N fertilization (R = 0.0401; P = 0.003) but not cover crop. The population distribution by activity was significantly affected by season (R = 0.2656; P = 0.001), cover crop (R = 0.0495; P = 0.001), and tillage (R = 0.0205; P = 0.016) but not fertilization (see Table S7 in the supplemental material).

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

NMDS of Bray-Curtis distances between N cycling communities based on abundances of five N cycling genes (G) and transcripts (T) (stress = 0.015). Points on the ordination plots are colored by season (A), cover crop treatment (B), tillage (C), and fertilization rate (D). NC, no cover; V, vetch; W, wheat; NT, no tillage; CT, conventional tillage; 0N, no N fertilization; 67N, 67 kg N ha−1 fertilization. Ellipses represent 95% confidence interval of each group.

Relationships of soil N pools and processes and soil properties with functional gene abundance and expression.Correlation heatmaps showed significant correlations among absolute abundances of N cycle genes and their transcripts, transcript/gene ratio, and soil properties (Fig. 7). We observed significant positive correlations among the five gene transcripts and significant positive correlations between N cycle genes and their corresponding transcripts for all gene targets except nifH (Fig. 7). Moreover, the transcript/gene ratios of N cycle genes were more related to their transcript abundances than gene abundances, reconfirming that transcript variation depended more on changes in expression rather than population sizes.

FIG 7
  • Open in new tab
  • Download powerpoint
FIG 7

Heatmap showing correlation among genes, gene transcripts, transcript/gene ratios, and soil properties. SWC, soil water content; TC, total carbon; TN, total nitrogen; TEC, total extractable carbon; TEN, total extractable nitrogen; MBC, microbial biomass carbon; MBN, microbial biomass nitrogen.

N cycle genes and their transcripts, as well as their transcript/gene ratios, showed different patterns of correlations with soil properties (Fig. 7). The gene and transcript abundances of nifH, nirK, and nirS were all positively correlated with SWC (Fig. 7). amoA gene and transcript abundances were both significantly and positively correlated with N2O-N, NO3−-N, field nitrification rate, and incubated nitrification rate (Fig. 7). Notably, amoA gene was positively correlated with NH4+-N (R = 0.494) while amoA transcript had no relationship with NH4+-N, and the transcript/gene ratio of amoA was negatively correlated with NH4+-N (R = −0.260), indicating that although amoA gene abundance was promoted by NH4+, its expression was lower with higher NH4+-N concentration (Fig. 7). In addition, there were significant positive correlations between N2O-N and the gene and transcript abundances of amoA and nosZ (Fig. 7).

DISCUSSION

Seasonal dynamics in abundance and activity of nitrogen cycling populations.The primary objective of this study was to investigate the impact of soil health management practices on the seasonal population size and activity of soil bacterial communities involved in N transformations by using nifH, amoA, nirK, nirS, and nosZ as indicators. We found that the N cycling bacterial community activity, as indicated by gene and transcript abundances, respectively, showed seasonal changes.

Without considering soil management practices, the highest functional activity of diazotrophs, as indicated by nifH expression, was in November during the fall crop transition, which might be explained by lower N availability at this time. The inhibitory effect of increased N availability (NH4+-N and/or NO3−-N) on nifH abundance and expression that we observed has been reported in other systems as well (45–48). Diazotrophs will preferentially use easily available N (e.g., nitrate and ammonium), as this takes less energy to assimilate compared to fixing dinitrogen gas (4). Soil freeze-thaw cycles, warming, and N addition can stimulate net N mineralization, increasing N availability (49). Compared to November, the warmer temperatures and N fertilization after cotton planting in May resulted in higher soil N availability during cotton growing seasons, thereby possibly inhibiting the expression of nifH. The cotton and winter cover crops take up inorganic and dissolved organic N for growth, leading to increased soil N limitation, especially during the cover crop growing period, which might enhance the competitiveness of diazotrophs and result in their relatively high abundances in April and November. Moreover, N fixation can be affected by soil oxygen content because nitrogenase is an extremely oxygen-sensitive enzyme (4). The high soil water content (SWC) in November may have resulted in lower soil oxygen, contributing to a higher functional activity of diazotrophs. The positive correlation between SWC and nifH abundance in our study was also observed in other studies (50, 51). In addition, the relatively high soil C:N ratio and the decomposition of a large amount of plant residues and roots after cotton harvest season (October) may have led to soil N limitation, stimulating the expression of N fixation genes.

The amoA gene relative and absolute abundances were both highest in October and lowest in November, which corresponded to the highest soil NH4+-N concentrations in October and lowest in November. Moreover, a significantly positive correlation between amoA gene abundances and NH4+-N concentrations was observed in our study, indicating that concentrations of NH4+-N, the substrate for nitrification, influence the abundance of nitrifiers (28, 35). The decrease in amoA gene abundance in November may be due to higher SWC and lower soil oxygen availability because AOB are aerobic bacteria (35). The seasonal dynamics of amoA transcript relative and absolute abundances both depended on the interaction of cover crops with either tillage or N fertilization rate. However, contrary to amoA gene abundance, relative amoA expression (normalized to gene abundance) was higher in November than in October and was negatively correlated to NH4+-N concentrations, indicating that while population sizes of AOB were down, their functional activity was increased. One hypothesis to explain this unexpected result was that ammonia oxidizers increased expression of their functional genes under oxygen and/or nutrient stress. We observed a significant positive correlation between expression level of amoA and SWC, which was also reported in other studies. For example, Theodorakopoulos et al. found that >70% water-filled pore space favored transcription of the AOB amoA gene (37). A pure culture experiment by Yu and Chandran found that a low level of dissolved oxygen can promote amoA transcription in AOB Nitrosomonas europaea (52). In addition, Bollmann et al. reported that during short-term (2 weeks) starvation, the amoA mRNA of Nitrosospira briensis remained present, and the ammonia-oxidation functionality was maintained at a high level (53).

The abundances of nirS genes and transcripts were greater than those of nirK in all agricultural seasons, which suggested that nirS-harboring denitrifiers were more predominant at our field site. Both relative and absolute transcript abundances and expressions of nirK- and nirS-harboring denitrifiers were highest in November, which is likely due to low soil oxygen level resulting from high SWC at that time. The significantly positive correlations of SWC with the transcript abundances and expression levels of nirK and nirS also support this explanation. Both relative and absolute abundances of nosZ genes and transcripts were highest in October. However, the N2O-N emissions were highest in both May and October. The not fully consistent seasonal trend between nosZ and N2O-N emission suggests that the typical nosZ (clade I nosZ) that we targeted in this study may not independently predict the seasonal dynamics of N2O emission, and atypical nosZ (clade II nosZ) may be playing an important role (54, 55). In fact, several studies have reported that atypical nosZ is more abundant than clade I nosZ in agricultural soils (56–58).

Effects of cover crops on abundance and activity of nitrogen cycling populations.Different responses were observed between the two types of cover crops. In general, compared to winter wheat, hairy vetch had significantly higher transcript abundance and transcript/gene ratios for all N cycle genes. The elevated activity of N cycling microbes under vetch was not limited to just the cover crop period (April) but persisted past cover crop termination through the cotton growing season (October). This supports the idea that leguminous cover crops have a lasting stimulation of overall N cycling microbes in agricultural soil. This is consistent with our previous findings from this site, where we documented significantly increased microbial biomass (both total and microbial biomass N) and enzyme activities in the vetch treatments (59). It may be that the biological N fixation rate stimulated by legumes and the low C:N ratio of legume residues supplied more N to the soil, providing a competitive advantage for N cycling populations (6, 60). Winter wheat is a nonleguminous cover crop with high C:N residues that takes up additional N from soil for growth, which may decrease soil N availability and inhibit the functional gene expression of N cycling microbes (except nifH). Some studies observed that leguminous cover crops can promote the release of mineral N and increase denitrification and N2O emissions relative to other cover crops or no cover (43, 61, 62). The relatively high and generally similar level of N cycle gene expression under combinations of hairy vetch with no fertilization and no cover crop with N fertilization supports the idea that, in some cases, N-fixing cover crops can replace N fertilization to maintain the N cycling rate.

Effects of tillage and fertilization on abundance and activity of nitrogen cycling populations.Overall, tillage only affected the relative abundances of genes nifH and nirS, but it affected the absolute abundances of all five functional genes and 16S rRNA genes, with higher abundances under conventional tillage than under no tillage. This may be because no-till treatment can decrease surface soil bulk density, increase soil porosity, and result in higher oxygen content in surface soils (63). The aerobic conditions resulting from conventional tillage may suppress the growth of anaerobic N cycling bacteria, such as N fixers and denitrifiers. Another possible explanation is that conventional tillage makes substrates more available for soil microbes by destroying soil aggregates (64, 65), increasing soil microbial populations, including the N cycling populations.

Compared to no fertilization, N fertilization with 67 kg N ha−1 only significantly increased the relative abundance of genes nirS and nosZ but only when no cover crops and wheat were used. This effect was observed only in wheat plots when absolute abundance was taken into consideration. This is partly in accordance with a meta-analysis by Ouyang et al., which showed that N fertilization had no effect on nifH abundance, but significantly increased the abundance of amoA, nirK, nirS, and nosZ genes, respectively (66). Although some studies using NH4NO3 as N fertilizer also reported that N fertilization can increase AOB diversity, abundance, and activity (30, 67, 68), we did not observe this in our study, likely because similar concentrations of NH4+-N and field nitrification rates were observed under both fertilizer treatments. Increased abundance of denitrification genes nirS and nosZ with N fertilization corresponded to high NO3−-N concentration and high N2O-N emission in fertilized soil.

Linking the nitrogen cycling groups with nitrogen cycling processes.In this study, we assessed the relationships between abundance and expression of functional genes and soil properties. We were particularly interested in determining whether N2O emissions, a potent greenhouse gas, could be correlated to the population size or activity of bacteria involved in N cycling processes responsible for N2O production and transformation. Both the absolute abundances of genes and transcripts of amoA and nosZ were positively correlated to N2O emissions. This suggests that soil AOB populations might be a main contributor to N2O production through nitrification, nitrifier denitrification, or both (69, 70). Likewise, Soares et al. also observed that AOB, rather than denitrifiers, were the main contributors to N2O emissions in tropical agricultural soils (71). Moreover, research by Theodorakopoulos et al. suggested that AOB amoA transcript abundances, rather than gene abundances, were a good indicator of N2O emissions, and ammonia oxidation was a major process contributing to N2O emissions (37). However, some studies have found no correlations between N2O emissions and abundance of related N cycle genes (AOB amoA, nirK, nirS, and nosZ) (36, 42). Therefore, the contribution of AOB and denitrifiers to N2O emissions might be site-specific and/or influenced by local soil properties and environmental factors. In addition, although we found significant positive correlation between field and incubated net nitrification rates with absolute abundances of both amoA genes and transcripts, a disconnection between AOB amoA gene abundances and nitrification rates has been observed in some studies (38, 67). This inconsistent result may be site-specific and might be influenced by the proportion of dormant microbes in soils, substrate limitation, microsite heterogeneity, and/or the fact that measured process rates reflect net, rather than gross, rates.

Conclusions.In summary, our results showed that the distribution patterns, abundances, and expressions of N cycle functional genes generally changed significantly between seasons and were affected by soil management practices. We demonstrated that leguminous cover crops like hairy vetch can promote functional activity of N cycling populations through the growing season, equaling or exceeding the effects of inorganic N fertilizer. Our findings highlight the fact that both genetic and transcriptional assessments can provide important information about microbial populations and functions. Transcriptional-level study has higher resolution than the genetic level in examining the effects of soil management practices on N cycling communities, with larger variation under different season and management practices observed. In addition, the interaction effects of season and soil management practices on N cycle gene expressions based on relative (normalized by 16S rRNA copies) and absolute (copies per gram soil) transcript abundances showed generally similar trends, indicating that active soil N cycling communities were more responsive than total active soil bacterial populations to agricultural management practices. Together, our work revealed insights into important relationships between microbial functional capacity and activity and associated nitrogen pools and processes in the field, adding to our understanding of the effects of soil health management practices on soil microbial ecology.

MATERIALS AND METHODS

Study site, experimental design, and sampling.This study was conducted at a long-term continuous cotton conservation agriculture field experiment, established in 1981, at the University of Tennessee West Tennessee Research and Education Center (WTREC) in Jackson, Tennessee. Soils at this site are classified as well-drained Lexington silt loams (fine-silty, mixed, thermic, Ultic Hapludalfs) with a 0 to 2% slope. The main crop was rainfed cotton (Gossypium hirsutum L.). The field study implemented a randomized complete block design with split-split plot treatment arrangements in four replicates such that N fertilization rate was the main plot treatment, cover crop was the subplot treatment, and tillage was the sub-subplot treatment; a subset of 12 of the 32 treatment combinations was used here (see Fig. S1 in the supplemental material). The 12 cropping system management combinations included 2 types of tillage (conventional tillage [CT], no tillage [NT]); 3 types of winter cover crops (hairy vetch [Vicia villosa Roth] [V], winter wheat [Triticum aestivum L.] [W], no cover crop [NC]); and 2 levels of NH4NO3 fertilization rates (0, 67 kg N ha−1 [0N, 67N, respectively]) (Table 2). For tilled plots, tillage was done twice before planting by a standard disc harrow followed by smoothing and breaking up of clods by a triple-K harrow. Field operations for these plots, including cover crop termination, cotton planting, and phosphorus and potassium fertilizer application are described in more detail in previous studies of this site (59, 72).

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 2

Summary of 12 cropping system treatments used in this study

FIG S1

Schematic diagram of the experimental cotton plots located at the University of Tennessee West Tennessee Research and Education Center in Jackson, Tennessee, USA. Download FIG S1, TIF file, 0.2 MB.
Copyright © 2021 Hu et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S2

16S rRNA normalized transcript abundances of nifH, amoA, nirK, nirS, and nosZ (A) and absolute transcript abundances of nifH, amoA, nirK, nirS, nosZ, and 16S rRNA (B) in relation to the combination of tillage and N fertilization (0 or 67 kg N ha−1 fertilization). Each value represents the mean ± standard error (n = 48). Download FIG S2, TIF file, 0.08 MB.
Copyright © 2021 Hu et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Subsamples (2.5-cm diameter by 7.5-cm-long cores) of soil were collected from 10 random locations in each experimental sub-subplot (12 cropping systems × 4 replicates = 48 experimental plots), composited, sieved through a 2-mm mesh, and then refrigerated until analysis. Collections were obtained four times corresponding to different stages of crop growth in 2017, including cover crop harvest (April, immediately following cover crop harvest), spring crop transition (May, after cotton planting and before fertilizer application), cotton crop peak (October, before cotton harvest), and fall crop transition (November, after cotton harvest and before winter cover crop planting). Soil properties measured included N2O flux (N2O-N), soil water content (SWC), soil pH, NO3−-N, NH4+-N, total extractable N (TEN), total extractable C (TEC), total soil C (TC), total soil N (TN), soil ratio of C to N (C:N ratio), microbial biomass N (MBN), microbial biomass C (MBC), field net nitrification rate (field nitrification), field net N mineralization rate (field N mineralization), incubated net nitrification rate (incubated nitrification), and incubated net N mineralization rate (incubated N mineralization). Detailed methods on the measurement of these soil properties are provided in the supplemental material (see Text S1, Methods S1 to S5).

TEXT S1

Description of methods used for soil property measurements. Download Text S1, DOCX file, 0.02 MB.
Copyright © 2021 Hu et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Soil DNA and RNA extraction.Soil genomic DNA was extracted from 0.25 g of soil using the DNeasy PowerSoil Kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. DNA concentrations were quantified using NanoDrop One spectrophotometry (NanoDrop Technologies, Wilmington, DE). DNA was stored at −20°C.

Soil RNA was extracted using the RNeasy PowerSoil total RNA kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. RNA was eluted in RNase-free water and stored at −80°C. RNA quality was checked by 1% agarose gel electrophoresis and NanoDrop One spectrophotometry (NanoDrop Technologies, Wilmington, DE). Samples were checked for DNA contamination by attempting to amplify the N cycle functional genes by PCR and subsequent agarose gel electrophoresis. Only samples that gave a negative electrophoresis result (i.e., no DNA contamination) were used for reverse transcription. cDNA was produced from RNA using SuperScript IV reverse transcriptase (Invitrogen, Paisley, UK) according to the manufacturer’s protocol. Random hexamer primers (Invitrogen) were used at a final concentration of 2.5 μM per reaction. cDNA samples were diluted to a final concentration of 15 ng μl−1 and stored at −20°C.

Molecular cloning.Genes nifH, amoA, nirK, nirS, and nosZ were amplified from extracted soil DNA with previously designed primer pairs and reaction conditions (see Table S1 in the supplemental material). PCR products were purified with QIAquick PCR purification kit (Qiagen), ligated into pCR 4-TOPO vectors (Invitrogen), and transformed into Transform One Shot TOP10 competent cells (Invitrogen) following the manufacturer’s protocol. After shaking at 225 rpm at 37°C for 1 h, cells were spread on LB solid medium containing 50 μg ml−1 ampicillin and cultured overnight. Positive clones were grown in liquid LB medium amended with ampicillin medium overnight, and plasmid DNA was extracted using the PureLink quick plasmid miniprep kit (Invitrogen) according to the manufacturer’s instructions. Plasmid DNA concentrations were determined on a Hoefer DQ 200 fluorometer (Hoefer Inc., San Francisco, CA), which uses a Hoechst dye specific for double-stranded DNA that compares DNA fluorescence with a known standard, and then sequenced with primers M13 Forward (−20)/M13 Reverse to confirm that the correct target gene was cloned. The copy numbers of nifH, amoA, nirK, nirS, and nosZ genes were calculated using the concentration of the extracted plasmid DNA, size (base pairs) of the plasmids containing cloned products, and the molecular weight of DNA. Ten-fold serial dilutions of the plasmid DNA (102 to 108 copies μl−1) were used during qPCR assays to generate standard curves for quantifying N cycle genes and transcripts in soil samples.

TABLE S1

Primers and PCR conditions used for qPCR and qRT-PCR. Download Table S1, DOCX file, 0.01 MB.
Copyright © 2021 Hu et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S2

Results of mixed model ANOVA (based on GLIMMIX procedure in SAS) testing effects of agricultural season and soil management practices on soil parameters. The significances were shown with bold fonts and asterisks. F values are reported. Download Table S2, DOCX file, 0.02 MB.
Copyright © 2021 Hu et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S3

Mean values of soil parameters in relation to agricultural season and soil management practices. Only treatments or treatment combinations that significantly affected at least one parameter were shown. Different lowercase letters indicate significant differences between treatment levels within groups. The letters were only shown in groups that have significant effects. Means were compared using least significant difference tests (α = 0.05). Download Table S3, DOCX file, 0.04 MB.
Copyright © 2021 Hu et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S4

Mean values of log-transformed relative abundances of N cycle genes and their transcripts (normalized by 16S rRNA genes and 16S rRNA) in relation to agricultural season and soil management practices. Download Table S4, DOCX file, 0.05 MB.
Copyright © 2021 Hu et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S5

Mean values of log-transformed mean absolute abundances of N cycle genes and their transcripts in relation to agricultural season and soil management practices. Download Table S5, DOCX file, 0.06 MB.
Copyright © 2021 Hu et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S6

Mean values of log-transformed transcript/gene ratios of N cycle genes in relation to agricultural season and soil management practices. Download Table S6, DOCX file, 0.03 MB.
Copyright © 2021 Hu et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S7

ANOSIM analysis of the distribution patterns of N cycling functional groups in four agricultural seasons and under different soil management practices. Download Table S7, DOCX file, 0.02 MB.
Copyright © 2021 Hu et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Quantitative analysis of N cycle genes and gene transcripts.For both DNA and cDNA samples, qPCR of nifH, amoA, nirK, nirS, and nosZ were performed on a CFX96 Optical real-time detection system (Bio-Rad, Laboratories Inc., Hercules, CA, USA). The 20-μl qPCR reaction mixture contained 10 μl Maxima SYBR green qPCR master mix (2×) (Thermo Scientific, USA), 1 μl PCR forward and reverse primer (both 10 μM), 2.5 μl DNA or cDNA template, and 5.5 μl nuclease-free water. Primer sets and reaction parameters used are listed in the supplemental material (see Table S1). Following the reactions, a melt curve analysis was performed to confirm the specificity of the PCR product for each real-time PCR amplification. Quantitative PCR of 16S rRNA genes and 16S rRNA was done using Femto bacterial DNA quantification kit (Zymo Research Corp., CA, USA) following the manufacturer’s protocols.

Both absolute abundances (copies per gram dry weight soil) and relative abundances (normalized to 16S rRNA genes or to 16S rRNA) of functional genes and transcripts were analyzed and compared to ensure the reliability of the result. Absolute abundance helped to achieve better understanding of the variation of N cycling bacteria with seasons and management practices, while relative abundance mitigated the impact of DNA and RNA extraction efficiency on data accuracy and revealed the relative changes of N cycling population within the total soil bacteria. The relative abundance of N cycle genes and transcripts used in this study were calculated as follows. NGA =N cycle gene copies/μl DNA extract16S rRNA gene copies/μl DNA extract NTA = N cycle transcript copies/μl RNA extract16S rRNA copies/μl RNA extract

where NGA is the relative abundance of the gene, and NTA is the relative abundance of the gene transcript.

The transcript-to-gene copy number ratio (transcript/gene ratio) for each gene was also calculated to provide a means to determine whether transcript variation was due to changes in population sizes (i.e., genes) or changes in expression (i.e., transcripts) (73) and was calculated as follows. Transcript/gene ratio =transcript copies/g dry weight soilgene copies/g dry weight soil

Statistical analysis.Statistical analyses to test the effects of agricultural season and management practices on the abundance and expression of N cycle genes were performed in SAS 9.4 (SAS Inst., Cary, NC, USA). The relative and absolute abundances of functional genes and transcripts as well as their transcript/gene ratio were both log transformed to achieve normal distributions of residuals. These log-transformed data were analyzed using a repeated-measures mixed model analysis of variance (ANOVA) within the GLIMMIX procedure. The fixed effects included season, tillage system, cover crop, and N fertilizer rate, as well as all interactions. Random effects included the block, block by nitrogen within cover crop, and the block by tillage by nitrogen within cover crop. Means were compared using a least significant difference (LSD) method. The effects of agriculture practices on N cycling microbes were visualized using scatterplots made by graphing software Origin Pro 2019b.

Pearson correlation analyses were used to determine correlation among genes, gene transcripts, transcript/gene ratios, and soil properties by rcorr function of package Hmisc in R 3.6.1 (the R Foundation for Statistical Computing). Correlation coefficient P values were then adjusted with a correction for multiple comparisons using the BH method of the Hmisc package. A heatmap based on Pearson correlations was used to visualize the significant positive and negative correlations among measured soil properties, N cycle gene abundance, gene transcript abundance, and transcript/gene ratios, and was performed by the pheatmap package in R 3.6.1.

Nonmetric multidimensional scaling (NMDS) based on Bray-Curtis distances was performed in R 3.6.1 with packages vegan 2.5-5 and phyloseq to reveal the distribution patterns of N cycling functional groups and their activity among seasons and management practices. Analysis of similarities (ANOSIM), a rank-based nonparametric statistical test, was also performed by R with the function anosim in package vegan to compare groups and test the null hypothesis that the similarity between groups is greater than or equal to the similarity within the groups. Dispersion indices were calculated in R with the function betadisper in package vegan to assess the multivariate homogeneity of group dispersions. The function TukeyHSD.betadisper in package vegan was used to calculate Tukey’s honest significant differences between groups. In this study, statistically significant difference was accepted at a P value of <0.05 unless otherwise noted.

ACKNOWLEDGMENTS

This research was supported by a USDA award (USDA-NIFA 2015-67019-23604) to S.M.S. and J.M.D. J.H. received financial support from the China Scholarship Council.

The authors are grateful to Brian Kozlowski and the staff of the West Tennessee Research and Education Center who manage and maintain the field trials; Megan Davis, Lydia Turpin, Lidong Li, and Kelly Cobaugh for field sampling and data collection; and Fei Yao for assistance in the laboratory. Finally, we gratefully acknowledge the late Don Tyler who, decades ago, had the vision to establish conservation agriculture trials that have provided invaluable information about long-term soil health management.

We declare no competing interests.

FOOTNOTES

    • Received December 2, 2020.
    • Accepted December 18, 2020.
  • Copyright © 2021 Hu et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. LeBauer DS,
    2. Treseder KK
    . 2008. Nitrogen limitation of net primary productivity in terrestrial ecosystems is globally distributed. Ecology 89:371–379. doi:10.1890/06-2057.1.
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    1. Galloway JN,
    2. Townsend AR,
    3. Erisman JW,
    4. Bekunda M,
    5. Cai Z,
    6. Freney JR,
    7. Martinelli LA,
    8. Seitzinger SP,
    9. Sutton MA
    . 2008. Transformation of the nitrogen cycle: recent trends, questions, and potential solutions. Science 320:889–892. doi:10.1126/science.1136674.
    OpenUrlAbstract/FREE Full Text
  3. 3.↵
    1. White PJ,
    2. Crawford JW,
    3. Díaz Álvarez MC,
    4. García Moreno R
    . 2012. Soil management for sustainable agriculture. Appl Environ Soil Sci 2012:850739. doi:10.1155/2012/850739.
    OpenUrlCrossRef
  4. 4.↵
    1. Chapin FS,
    2. Matson PA,
    3. Mooney HA
    . 2002. Principles of terrestrial ecosystem ecology. Springer, New York, NY.
  5. 5.↵
    1. Parkin T,
    2. Kaspar T
    . 2006. Nitrous oxide emissions from corn-soybean systems in the Midwest. J Environ Qual 35:1496–1506. doi:10.2134/jeq2005.0183.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    1. Carsky RJ,
    2. Singh BB,
    3. Oyewole B
    . 2001. Contribution of early season cowpea to late season maize in the savanna zone of West Africa. Biol Agric Hortic 18:303–315. doi:10.1080/01448765.2001.9754894.
    OpenUrlCrossRef
  7. 7.↵
    1. van Kessel C,
    2. Venterea R,
    3. Six J,
    4. Adviento-Borbe MA,
    5. Linquist B,
    6. van Groenigen KJ
    . 2013. Climate, duration, and N placement determine N2O emissions in reduced tillage systems: a meta-analysis. Glob Change Biol 19:33–44. doi:10.1111/j.1365-2486.2012.02779.x.
    OpenUrlCrossRefPubMedWeb of Science
  8. 8.↵
    1. Feng J,
    2. Li F,
    3. Zhou X,
    4. Xu C,
    5. Ji L,
    6. Chen Z,
    7. Fang F
    . 2018. Impact of agronomy practices on the effects of reduced tillage systems on CH4 and N2O emissions from agricultural fields: a global meta-analysis. PLoS One 13:e0196703. doi:10.1371/journal.pone.0196703.
    OpenUrlCrossRef
  9. 9.↵
    1. Forbes MS,
    2. Broos K,
    3. Baldock JA,
    4. Gregg AL,
    5. Wakelin SA
    . 2009. Environmental and edaphic drivers of bacterial communities involved in soil N-cycling. Soil Res 47:380–388. doi:10.1071/SR08126.
    OpenUrlCrossRef
  10. 10.↵
    1. Cambardella CA,
    2. Johnson JM,
    3. Varvel GE
    . 2012. Soil carbon sequestration in central USA agroecosystems, p 41–58. In Liebig MA, Franzluebbers AJ, Follett RF (ed), Managing agricultural greenhouse gases: coordinated agricultural research through GRACEnet to address our changing climate. Elsevier, San Diego, CA.
  11. 11.↵
    1. Linn DM,
    2. Doran JW
    . 1984. Aerobic and anaerobic microbial populations in no-till and plowed soils. Soil Sci Soc Am J 48:794–799. doi:10.2136/sssaj1984.03615995004800040019x.
    OpenUrlCrossRefWeb of Science
  12. 12.↵
    1. Yang D,
    2. Xiao X,
    3. He N,
    4. Zhu W,
    5. Liu M,
    6. Xie G
    . 2020. Effects of reducing chemical fertilizer combined with organic amendments on ammonia-oxidizing bacteria and archaea communities in a low-fertility red paddy field. Environ Sci Pollut Res Int 27:29422–29432. doi:10.1007/s11356-020-09120-5.
    OpenUrlCrossRef
  13. 13.↵
    1. Bahulikar RA,
    2. Chaluvadi SR,
    3. Torres-Jerez I,
    4. Mosali J,
    5. Bennetzen JL,
    6. Udvardi M
    . 22 January 2020. Nitrogen fertilization reduces nitrogen fixation activity of diverse diazotrophs in switchgrass roots. Phytobiomes J doi:10.1094/pbiomes-09-19-0050-fi.
    OpenUrlCrossRef
  14. 14.↵
    1. Wang J,
    2. Chadwick DR,
    3. Cheng Y,
    4. Yan X
    . 2018. Global analysis of agricultural soil denitrification in response to fertilizer nitrogen. Sci Total Environ 616–617:908–917. doi:10.1016/j.scitotenv.2017.10.229.
    OpenUrlCrossRef
  15. 15.↵
    1. Kuypers MMM,
    2. Marchant HK,
    3. Kartal B
    . 2018. The microbial nitrogen-cycling network. Nat Rev Microbiol 16:263–276. doi:10.1038/nrmicro.2018.9.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Orr CH,
    2. James A,
    3. Leifert C,
    4. Cooper JM,
    5. Cummings SP
    . 2011. Diversity and activity of free-living nitrogen-fixing bacteria and total bacteria in organic and conventionally managed soils. Appl Environ Microbiol 77:911–919. doi:10.1128/AEM.01250-10.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. Akter Z,
    2. Pageni BB,
    3. Lupwayi NZ,
    4. Balasubramanian PM
    . 2014. Biological nitrogen fixation and nifH gene expression in dry beans (Phaseolus vulgaris L.). Can J Plant Sci 94:203–212. doi:10.4141/cjps2013-200.
    OpenUrlCrossRef
  18. 18.↵
    1. Gaby JC,
    2. Buckley DH
    . 2014. A comprehensive aligned nifH gene database: a multipurpose tool for studies of nitrogen-fixing bacteria. Database (Oxford) 2014:bau001. doi:10.1093/database/bau001.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Liu R,
    2. Hu H,
    3. Suter H,
    4. Hayden HL,
    5. He J,
    6. Mele P,
    7. Chen D
    . 2016. Nitrification Is a primary driver of nitrous oxide production in laboratory microcosms from different land-use soils. Front Microbiol 7:1373. doi:10.3389/fmicb.2016.01373.
    OpenUrlCrossRef
  20. 20.↵
    1. Inatomi M,
    2. Hajima T,
    3. Ito A
    . 2019. Fraction of nitrous oxide production in nitrification and its effect on total soil emission: a meta-analysis and global-scale sensitivity analysis using a process-based model. PLoS One 14:e0219159. doi:10.1371/journal.pone.0219159.
    OpenUrlCrossRef
  21. 21.↵
    1. Khalil K,
    2. Mary B,
    3. Renault P
    . 2004. Nitrous oxide production by nitrification and denitrification in soil aggregates as affected by O2 concentration. Soil Biol Biochem 36:687–699. doi:10.1016/j.soilbio.2004.01.004.
    OpenUrlCrossRefWeb of Science
  22. 22.↵
    1. White CJ,
    2. Lehnert N
    . 2016. Is there a pathway for N2O production from hydroxylamine oxidoreductase in ammonia-oxidizing bacteria? Proc Natl Acad Sci U S A 113:14474–14476. doi:10.1073/pnas.1617953114.
    OpenUrlFREE Full Text
  23. 23.↵
    1. Bárta J,
    2. Tahovská K,
    3. Šantrůčková H,
    4. Oulehle F
    . 2017. Microbial communities with distinct denitrification potential in spruce and beech soils differing in nitrate leaching. Sci Rep 7:9738. doi:10.1038/s41598-017-08554-1.
    OpenUrlCrossRef
  24. 24.↵
    1. Ward BB,
    2. Bouskill NJ
    . 2011. The utility of functional gene arrays for assessing community composition, relative abundance, and distribution of ammonia-oxidizing bacteria and archaea. Methods Enzymol 496:373–396. doi:10.1016/B978-0-12-386489-5.00015-4.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Norton JM,
    2. Low JM,
    3. Klotz MG
    . 1996. The gene encoding ammonia monooxygenase subunit A exists in three nearly identical copies in Nitrosospira sp. NpAV. FEMS Microbiol Lett 139:181–188. doi:10.1111/j.1574-6968.1996.tb08200.x.
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    1. Leininger S,
    2. Urich T,
    3. Schloter M,
    4. Schwark L,
    5. Qi J,
    6. Nicol GW,
    7. Prosser JI,
    8. Schuster SC,
    9. Schleper C
    . 2006. Archaea predominate among ammonia-oxidizing prokaryotes in soils. Nature 442:806–809. doi:10.1038/nature04983.
    OpenUrlCrossRefPubMedWeb of Science
  27. 27.↵
    1. Ouyang Y,
    2. Norton JM,
    3. Stark JM,
    4. Reeve JR,
    5. Habteselassie MY
    . 2016. Ammonia-oxidizing bacteria are more responsive than archaea to nitrogen source in an agricultural soil. Soil Biol Biochem 96:4–15. doi:10.1016/j.soilbio.2016.01.012.
    OpenUrlCrossRef
  28. 28.↵
    1. Ouyang Y,
    2. Norton JM,
    3. Stark JM
    . 2017. Ammonium availability and temperature control contributions of ammonia oxidizing bacteria and archaea to nitrification in an agricultural soil. Soil Biol Biochem 113:161–172. doi:10.1016/j.soilbio.2017.06.010.
    OpenUrlCrossRef
  29. 29.↵
    1. Ying J,
    2. Li X,
    3. Wang N,
    4. Lan Z,
    5. He J,
    6. Bai Y
    . 2017. Contrasting effects of nitrogen forms and soil pH on ammonia oxidizing microorganisms and their responses to long-term nitrogen fertilization in a typical steppe ecosystem. Soil Biol Biochem 107:10–18. doi:10.1016/j.soilbio.2016.12.023.
    OpenUrlCrossRef
  30. 30.↵
    1. Segal LM,
    2. Miller DN,
    3. McGhee RP,
    4. Loecke TD,
    5. Cook KL,
    6. Shapiro CA,
    7. Drijber RA
    . 2017. Bacterial and archaeal ammonia oxidizers respond differently to long-term tillage and fertilizer management at a continuous maize site. Soil Till Res 168:110–117. doi:10.1016/j.still.2016.12.014.
    OpenUrlCrossRef
  31. 31.↵
    1. Jia Z,
    2. Conrad R
    . 2009. Bacteria rather than archaea dominate microbial ammonia oxidation in an agricultural soil. Environ Microbiol 11:1658–1671. doi:10.1111/j.1462-2920.2009.01891.x.
    OpenUrlCrossRefPubMedWeb of Science
  32. 32.↵
    1. Colliver BB,
    2. Stephenson T
    . 2000. Production of nitrogen oxide and dinitrogen oxide by autotrophic nitrifiers. Biotechnol Adv 18:219–232. doi:10.1016/s0734-9750(00)00035-5.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    1. Sanchez C,
    2. Minamisawa K
    . 2018. Redundant roles of Bradyrhizobium oligotrophicum Cu-type (NirK) and cd1-type (NirS) nitrite reductase genes under denitrifying conditions. FEMS Microbiol Lett 365:fny015. doi:10.1093/femsle/fny015.
    OpenUrlCrossRef
  34. 34.↵
    1. Chan Y-K,
    2. McCormic WA,
    3. Watson RJ
    . 1997. A new nos gene downstream from nosDFY is essential for dissimilatory reduction of nitrous oxide by Rhizobium (Sinorhizobium) meliloti. Microbiology 143:2817–2824. doi:10.1099/00221287-143-8-2817.
    OpenUrlCrossRefPubMedWeb of Science
  35. 35.↵
    1. Levy-Booth DJ,
    2. Prescott CE,
    3. Grayston SJ
    . 2014. Microbial functional genes involved in nitrogen fixation, nitrification and denitrification in forest ecosystems. Soil Biol Biochem 75:11–25. doi:10.1016/j.soilbio.2014.03.021.
    OpenUrlCrossRef
  36. 36.↵
    1. Pereira EIP,
    2. Suddick EC,
    3. Mansour I,
    4. Mukome FND,
    5. Parikh SJ,
    6. Scow K,
    7. Six J
    . 2015. Biochar alters nitrogen transformations but has minimal effects on nitrous oxide emissions in an organically managed lettuce mesocosm. Biol Fertil Soils 51:573–582. doi:10.1007/s00374-015-1004-5.
    OpenUrlCrossRef
  37. 37.↵
    1. Theodorakopoulos N,
    2. Lognoul M,
    3. Degrune F,
    4. Broux F,
    5. Regaert D,
    6. Muys C,
    7. Heinesch B,
    8. Bodson B,
    9. Aubinet M,
    10. Vandenbol M
    . 2017. Increased expression of bacterial amoA during an N2O emission peak in an agricultural field. Agric Ecosyst Environ 236:212–220. doi:10.1016/j.agee.2016.12.002.
    OpenUrlCrossRef
  38. 38.↵
    1. Li S,
    2. Jiang X,
    3. Wang X,
    4. Wright AL
    . 2015. Tillage effects on soil nitrification and the dynamic changes in nitrifying microorganisms in a subtropical rice-based ecosystem: a long-term field study. Soil Till Res 150:132–138. doi:10.1016/j.still.2015.02.005.
    OpenUrlCrossRef
  39. 39.↵
    1. Tatti E,
    2. Goyer C,
    3. Burton D,
    4. Wertz S,
    5. Zebarth B,
    6. Chantigny M,
    7. Filion M
    . 2015. Tillage management and seasonal effects on denitrifier community abundance, gene expression and structure over winter. Microb Ecol 70:795–808. doi:10.1007/s00248-015-0591-x.
    OpenUrlCrossRef
  40. 40.↵
    1. Badagliacca G,
    2. Benitez E,
    3. Amato G,
    4. Badalucco L,
    5. Giambalvo D,
    6. Laudicina VA,
    7. Ruisi P
    . 2018. Long-term no-tillage application increases soil organic carbon, nitrous oxide emissions and faba bean (Vicia faba L.) yields under rain-fed Mediterranean conditions. Sci Total Environ 639:350–359. doi:10.1016/j.scitotenv.2018.05.157.
    OpenUrlCrossRef
  41. 41.↵
    1. Wang J,
    2. Zou J
    . 2020. No-till increases soil denitrification via its positive effects on the activity and abundance of the denitrifying community. Soil Biol Biochem 142:107706. doi:10.1016/j.soilbio.2020.107706.
    OpenUrlCrossRef
  42. 42.↵
    1. Maul JE,
    2. Cavigelli MA,
    3. Vinyard B,
    4. Buyer JS
    . 2019. Cropping system history and crop rotation phase drive the abundance of soil denitrification genes nirK, nirS and nosZ in conventional and organic grain agroecosystems. Agric Ecosyst Environ 273:95–106. doi:10.1016/j.agee.2018.11.022.
    OpenUrlCrossRef
  43. 43.↵
    1. Bowen H,
    2. Maul JE,
    3. Poffenbarger H,
    4. Mirsky S,
    5. Cavigelli M,
    6. Yarwood S
    . 2018. Spatial patterns of microbial denitrification genes change in response to poultry litter placement and cover crop species in an agricultural soil. Biol Fertil Soils 54:769–781. doi:10.1007/s00374-018-1301-x.
    OpenUrlCrossRef
  44. 44.↵
    1. Wang C,
    2. Zheng M,
    3. Song W,
    4. Wen S,
    5. Wang B,
    6. Zhu C,
    7. Shen R
    . 2017. Impact of 25 years of inorganic fertilization on diazotrophic abundance and community structure in an acidic soil in southern China. Soil Biol Biochem 113:240–249. doi:10.1016/j.soilbio.2017.06.019.
    OpenUrlCrossRef
  45. 45.↵
    1. Cusack DF,
    2. Silver W,
    3. McDowell WH
    . 2009. Biological nitrogen fixation in two tropical forests: ecosystem-level patterns and effects of nitrogen fertilization. Ecosystems 12:1299–1315. doi:10.1007/s10021-009-9290-0.
    OpenUrlCrossRef
  46. 46.↵
    1. Lindsay EA,
    2. Colloff MJ,
    3. Gibb NL,
    4. Wakelin SA
    . 2010. The abundance of microbial functional genes in grassy woodlands is influenced more by soil nutrient enrichment than by recent weed invasion or livestock exclusion. Appl Environ Microbiol 76:5547–5555. doi:10.1128/AEM.03054-09.
    OpenUrlAbstract/FREE Full Text
  47. 47.↵
    1. Pereira ESMC,
    2. Schloter-Hai B,
    3. Schloter M,
    4. van Elsas JD,
    5. Salles JF
    . 2013. Temporal dynamics of abundance and composition of nitrogen-fixing communities across agricultural soils. PLoS One 8:e74500. doi:10.1371/journal.pone.0074500.
    OpenUrlCrossRef
  48. 48.↵
    1. Reardon CL,
    2. Gollany HT,
    3. Wuest SB
    . 2014. Diazotroph community structure and abundance in wheat-fallow and wheat-pea crop rotations. Soil Biol Biochem 69:406–412. doi:10.1016/j.soilbio.2013.10.038.
    OpenUrlCrossRef
  49. 49.↵
    1. Contosta AR,
    2. Frey SD,
    3. Cooper AB
    . 2011. Seasonal dynamics of soil respiration and N mineralization in chronically warmed and fertilized soils. Ecosphere 2:art36. doi:10.1890/ES10-00133.1.
    OpenUrlCrossRef
  50. 50.↵
    1. Chen J,
    2. Shen W,
    3. Xu H,
    4. Li Y,
    5. Luo T
    . 2019. The composition of nitrogen-fixing microorganisms correlates with soil nitrogen content during reforestation: a comparison between legume and non-legume plantations. Front Microbiol 10:508. doi:10.3389/fmicb.2019.00508.
    OpenUrlCrossRef
  51. 51.↵
    1. Zhang X,
    2. Jia X,
    3. Wu H,
    4. Li J,
    5. Yan L,
    6. Wang J,
    7. Li Y,
    8. Kang X
    . 2020. Depression of soil nitrogen fixation by drying soil in a degraded alpine peatland. Sci Total Environ 747:141084. doi:10.1016/j.scitotenv.2020.141084.
    OpenUrlCrossRef
  52. 52.↵
    1. Yu R,
    2. Chandran K
    . 2010. Strategies of Nitrosomonas europaea 19718 to counter low dissolved oxygen and high nitrite concentrations. BMC Microbiol 10:70. doi:10.1186/1471-2180-10-70.
    OpenUrlCrossRefPubMed
  53. 53.↵
    1. Bollmann A,
    2. Schmidt I,
    3. Saunders AM,
    4. Nicolaisen MH
    . 2005. Influence of starvation on potential ammonia-oxidizing activity and amoA mRNA levels of Nitrosospira briensis. Appl Environ Microbiol 71:1276–1282. doi:10.1128/AEM.71.3.1276-1282.2005.
    OpenUrlAbstract/FREE Full Text
  54. 54.↵
    1. Sanford RA,
    2. Wagner DD,
    3. Wu Q,
    4. Chee-Sanford JC,
    5. Thomas SH,
    6. Cruz-Garcia C,
    7. Rodriguez G,
    8. Massol-Deya A,
    9. Krishnani KK,
    10. Ritalahti KM,
    11. Nissen S,
    12. Konstantinidis KT,
    13. Loffler FE
    . 2012. Unexpected nondenitrifier nitrous oxide reductase gene diversity and abundance in soils. Proc Natl Acad Sci U S A 109:19709–19714. doi:10.1073/pnas.1211238109.
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    1. Orellana LH,
    2. Rodriguez RL,
    3. Higgins S,
    4. Chee-Sanford JC,
    5. Sanford RA,
    6. Ritalahti KM,
    7. Loffler FE,
    8. Konstantinidis KT
    . 2014. Detecting nitrous oxide reductase (nosZ) genes in soil metagenomes: method development and implications for the nitrogen cycle. mBio 5:e01193-14. doi:10.1128/mBio.01193-14.
    OpenUrlAbstract/FREE Full Text
  56. 56.↵
    1. Jones CM,
    2. Graf DR,
    3. Bru D,
    4. Philippot L,
    5. Hallin S
    . 2013. The unaccounted yet abundant nitrous oxide-reducing microbial community: a potential nitrous oxide sink. ISME J 7:417–426. doi:10.1038/ismej.2012.125.
    OpenUrlCrossRefPubMed
  57. 57.↵
    1. Tsiknia M,
    2. Paranychianakis NV,
    3. Varouchakis EA,
    4. Nikolaidis NP,
    5. Sessitsch A
    . 2015. Environmental drivers of the distribution of nitrogen functional genes at a watershed scale. FEMS Microbiol Ecol 91:fiv052. doi:10.1093/femsec/fiv052.
    OpenUrlCrossRefPubMed
  58. 58.↵
    1. Juhanson J,
    2. Hallin S,
    3. Söderström M,
    4. Stenberg M,
    5. Jones CM
    . 2017. Spatial and phyloecological analyses of nosZ genes underscore niche differentiation amongst terrestrial N2O reducing communities. Soil Biol Biochem 115:82–91. doi:10.1016/j.soilbio.2017.08.013.
    OpenUrlCrossRef
  59. 59.↵
    1. Mbuthia LW,
    2. Acosta-Martínez V,
    3. DeBruyn J,
    4. Schaeffer S,
    5. Tyler D,
    6. Odoi E,
    7. Mpheshea M,
    8. Walker F,
    9. Eash N
    . 2015. Long term tillage, cover crop, and fertilization effects on microbial community structure, activity: implications for soil quality. Soil Biol Biochem 89:24–34. doi:10.1016/j.soilbio.2015.06.016.
    OpenUrlCrossRef
  60. 60.↵
    1. Wang H,
    2. Li X,
    3. Li X,
    4. Li X,
    5. Wang J,
    6. Zhang H
    . 2017. Changes of microbial population and N-cycling function genes with depth in three Chinese paddy soils. PLoS One 12:e0189506. doi:10.1371/journal.pone.0189506.
    OpenUrlCrossRef
  61. 61.↵
    1. Abdalla M,
    2. Rueangritsarakul K,
    3. Jones M,
    4. Osborne B,
    5. Helmy M,
    6. Roth B,
    7. Burke J,
    8. Nolan P,
    9. Smith P,
    10. Williams M
    . 2012. How effective is reduced tillage-cover crop management in reducing N2O fluxes from arable crop soils? Water Air Soil Pollut 223:5155–5174. doi:10.1007/s11270-012-1268-4.
    OpenUrlCrossRef
  62. 62.↵
    1. Hauggaard-Nielsen H,
    2. Lachouani P,
    3. Knudsen MT,
    4. Ambus P,
    5. Boelt B,
    6. Gislum R
    . 2016. Productivity and carbon footprint of perennial grass-forage legume intercropping strategies with high or low nitrogen fertilizer input. Sci Total Environ 541:1339–1347. doi:10.1016/j.scitotenv.2015.10.013.
    OpenUrlCrossRef
  63. 63.↵
    1. López-Vázquez A,
    2. Cadena-Zapata M,
    3. Campos-Magaña S,
    4. Zermeño-Gonzalez A,
    5. Mendez-Dorado M
    . 2019. Comparison of energy used and effects on bulk density and yield by tillage systems in a semiarid condition of Mexico. Agronomy 9:189. doi:10.3390/agronomy9040189.
    OpenUrlCrossRef
  64. 64.↵
    1. Elliott ET
    . 1986. Aggregate structure and carbon, nitrogen, and phosphorus in native and cultivated soils. Soil Sci Soc Am J 50:627–633. doi:10.2136/sssaj1986.03615995005000030017x.
    OpenUrlCrossRefWeb of Science
  65. 65.↵
    1. Ashman MR,
    2. Hallett PD,
    3. Brookes PC
    . 2003. Are the links between soil aggregate size class, soil organic matter and respiration rate artefacts of the fractionation procedure? Soil Biol Biochem 35:435–444. doi:10.1016/S0038-0717(02)00295-X.
    OpenUrlCrossRef
  66. 66.↵
    1. Ouyang Y,
    2. Evans SE,
    3. Friesen ML,
    4. Tiemann LK
    . 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. doi:10.1016/j.soilbio.2018.08.024.
    OpenUrlCrossRef
  67. 67.↵
    1. Liu S,
    2. Coyne MS,
    3. Grove JH
    . 2017. Long-term tillage and nitrogen fertilization: consequences for nitrifier density and activity. Appl Soil Ecol 120:121–127. doi:10.1016/j.apsoil.2017.07.034.
    OpenUrlCrossRef
  68. 68.↵
    1. Liu S,
    2. Coyne MS,
    3. Grove JH,
    4. Flythe MD
    . 2018. Nitrogen, season, and tillage management influence ammonia oxidizing bacterial communities in long-term maize. Appl Soil Ecol 129:98–106. doi:10.1016/j.apsoil.2018.05.002.
    OpenUrlCrossRef
  69. 69.↵
    1. Zhu X,
    2. Burger M,
    3. Doane TA,
    4. Horwath WR
    . 2013. Ammonia oxidation pathways and nitrifier denitrification are significant sources of N2O and NO under low oxygen availability. Proc Natl Acad Sci U S A 110:6328–6333. doi:10.1073/pnas.1219993110.
    OpenUrlAbstract/FREE Full Text
  70. 70.↵
    1. Huang T,
    2. Gao B,
    3. Hu XK,
    4. Lu X,
    5. Well R,
    6. Christie P,
    7. Bakken LR,
    8. Ju XT
    . 2014. Ammonia-oxidation as an engine to generate nitrous oxide in an intensively managed calcareous fluvo-aquic soil. Sci Rep 4:3950. doi:10.1038/srep03950.
    OpenUrlCrossRefPubMed
  71. 71.↵
    1. Soares JR,
    2. Cassman NA,
    3. Kielak AM,
    4. Pijl A,
    5. Carmo JB,
    6. Lourenco KS,
    7. Laanbroek HJ,
    8. Cantarella H,
    9. Kuramae EE
    . 2016. Nitrous oxide emission related to ammonia-oxidizing bacteria and mitigation options from N fertilization in a tropical soil. Sci Rep 6:30349. doi:10.1038/srep30349.
    OpenUrlCrossRef
  72. 72.↵
    1. Nouri A,
    2. Lee J,
    3. Yin X,
    4. Tyler DD,
    5. Saxton AM
    . 2019. Thirty-four years of no-tillage and cover crops improve soil quality and increase cotton yield in Alfisols, Southeastern USA. Geoderma 337:998–1008. doi:10.1016/j.geoderma.2018.10.016.
    OpenUrlCrossRef
  73. 73.↵
    1. Cusick KD,
    2. Fitzgerald LA,
    3. Cockrell AL,
    4. Biffinger JC
    . 2015. Selection and evaluation of reference genes for reverse transcription-quantitative PCR expression studies in a thermophilic bacterium grown under different culture conditions. PLoS One 10:e0131015. doi:10.1371/journal.pone.0131015.
    OpenUrlCrossRef
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
Soil Health Management Enhances Microbial Nitrogen Cycling Capacity and Activity
Jialin Hu, Virginia L. Jin, Julie Y. M. Konkel, Sean M. Schaeffer, Liesel G. Schneider, Jennifer M. DeBruyn
mSphere Jan 2021, 6 (1) e01237-20; DOI: 10.1128/mSphere.01237-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSphere article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Soil Health Management Enhances Microbial Nitrogen Cycling Capacity and Activity
(Your Name) has forwarded a page to you from mSphere
(Your Name) thought you would be interested in this article in mSphere.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Soil Health Management Enhances Microbial Nitrogen Cycling Capacity and Activity
Jialin Hu, Virginia L. Jin, Julie Y. M. Konkel, Sean M. Schaeffer, Liesel G. Schneider, Jennifer M. DeBruyn
mSphere Jan 2021, 6 (1) e01237-20; DOI: 10.1128/mSphere.01237-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

soil health
conservation agriculture
soil nitrogen cycle
soil microbes
quantitative PCR (qPCR)
quantitative reverse transcription PCR (qRT-PCR)
agroecosystems
denitrification
nitrification
nitrogen fixation
soil microbiology

Related Articles

Cited By...

About

  • About mSphere
  • Board of Editors
  • Policies
  • For Reviewers
  • For the Media
  • Embargo Policy
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSphereJ

@ASMicrobiology

       

 

Website feedback

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5042