Changes in N-Transforming Archaea and Bacteria in Soil during the Establishment of Bioenergy Crops
2011

Y Mao
A C Yannarell
R I Mackie
Changes in N-Transforming Archaea and Bacteria in Soil during the Establishment of Bioenergy Crops

PLoS ONE
6924750201110.1371/journal.pone.0024750Received May 24, 2011; Accepted August 16, 2011;Editor: Jack Anthony Gilbert, Argonne National Laboratory, United States of America Funding: This work was funded by the Energy Biosciences Institute, Environmental Impact and Sustainability of Feedstock Production Program at the University of Illinois, Urbana. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist. *
Widespread adaptation of biomass production for bioenergy may influence important biogeochemical functions in the landscape, which are mainly carried out by soil microbes. Here we explore the impact of four potential bioenergy feedstock crops (maize, switchgrass, Miscanthus X giganteus, and mixed tallgrass prairie) on nitrogen cycling microorganisms in the soil by monitoring the changes in the quantity (real-time PCR) and diversity (barcoded pyrosequencing) of key functional genes (nifH, bacterial/archaeal amoA and nosZ) and 16S rRNA genes over two years after bioenergy crop establishment. The quantities of these N-cycling genes were relatively stable in all four crops, except maize (the only fertilized crop), in which the population size of AOB doubled in less than 3 months. The nitrification rate was significantly correlated with the quantity of ammonia-oxidizing archaea (AOA) not bacteria (AOB), indicating that archaea were the major ammonia oxidizers. Deep sequencing revealed high diversity of nifH, archaeal amoA, bacterial amoA, nosZ and 16S rRNA genes, with 229, 309, 330, 331 and 8989 OTUs observed, respectively. Rarefaction analysis revealed the diversity of archaeal amoA in maize markedly decreased in the second year. Ordination analysis of T-RFLP and pyrosequencing results showed that the Ntransforming microbial community structures in the soil under these crops gradually differentiated. Thus far, our two-year study has shown that specific N-transforming microbial communities develop in the soil in response to planting different bioenergy crops, and each functional group responded in a different way. Our results also suggest that cultivation of maize with N-fertilization increases the abundance of AOB and denitrifiers, reduces the diversity of AOA, and results in significant changes in the structure of denitrification community.

Introduction

Bioenergy derived from cellulosic ethanol is a potential sustainable alternative to fossil fuel-based energy, since the energy from green plants is renewable and largely carbon neutral in comparison to fossil fuel combustion. Perennial grasses, such as switchgrass (Panicum virgatum) and Miscanthus6giganteus, with large annual biomass production potential, are proposed as biofuel feedstocks that can maximize ethanol production without adversely affecting the market for food crops (e.g. maize). However, our knowledge of the impacts of various bioenergy feedstock production systems on the soil microbial ecosystem is still very limited. The chemistry of perennial crop residues and plant root exudates may stimulate or inhibit the growth and activity of different fractions of the soil microbial community, and thus the planting of different crops can result in distinct microbial communities [1,2,3]. Differences in management techniques between traditional row-crop agriculture and perennial biomass feedstocks represent different soil disturbance regimes, altered water use, differing rates of fertilizer application, etc., and these should have a direct impact on soil microbial dynamics, subsequently influencing the terrestrial biogeochemical cycles. In particular, we predict that the cultivation of high nitrogen-use efficiency perennial grasses will result in altered nitrogentransforming microbial communities in comparison to those found under N-fertilized maize.

The biological nitrogen cycle is one of the most important nutrient cycles in the terrestrial ecosystem. It includes four major processes: nitrogen fixation, mineralization (decay), nitrification and denitrification. Because many of the microorganisms responsible for these processes are recalcitrant to laboratory cultivation, previous studies of the distribution and diversity of nitrogentransforming microorganisms have employed cultivation-independent techniques targeting functional genes: nifH, amoA and nosZ genes, which encode the key enzymes in nitrogen fixation, ammonia oxidization and complete denitrification, respectively [4,5,6,7,8,9].

Biological nitrogen fixation, which converts atmospheric N 2 into ammonium that is available to organisms, is an important natural input of available nitrogen in many terrestrial habitats [10]. Although nitrogen fixation in terrestrial ecosystems is thought to be mainly carried out by the symbiotic bacteria in association with plants, free-living diazotrophs in soils can play important roles in N cycling in a number of ecosystems [11,12]. In average, 2-3 kg N ha 21 year 21 could be imported by free living N-fixers [13]. Various field experiments have shown that the biomass yield of one candidate biofuel feedstock crop, Miscanthus6giganteus, is not significantly increased by the addition of mineral N fertilizer [14]. The lack of response to nitrogen fertilization and the high biomass production suggest that biological nitrogen fixation may play an important role in supplying the nitrogen needs of Miscanthus [15]. Plant species have previously been shown to have a significant effect on the composition of diazotrophs in the field; for example, diazotroph diversity is higher in soil under Acacia tortilis ssp. raddiana (a leguminous tree) than Balanites aegyptiaca (a non-leguminous tree) [16]. Plant genotype also has a strong effect on the rhizosphere diazotrophs of rice [17]. Agronomic practices can also influence soil diazotrophs, e.g. application of N-fertilization can reduce the diversity of diazotrophs [17]. Therefore, we hypothesize that the cultivation of maize with inorganic N-fertilizer will reduce the abundance and diversity of diazotrophs in the soil ecosystem, while biofuel feedstocks receiving little or no N-fertilizer (e.g. Miscanthus) will encourage the development of active diazotrophic communities.

Nitrification, which converts ammonium to nitrate, includes two steps: ammonia oxidation to nitrite, and nitrite oxidation to nitrate. The production of nitrate in soil not only supplies nutrition for plants, but it can also mobilize nitrogen to groundwater through nitrate leaching. Ammonia oxidation is the first and rate-limiting step of nitrification [18]. It is typically thought to be carried out by a few groups in band c-Proteobacteria, referred to as ammoniaoxidizing bacteria (AOB) [18]. However, recent environmental metagenomic analyses revealed that ammonia monooxygenase asubunit (amoA) genes are also present in archaea (AOA) [19], and archaeal amoA has been shown to be widespread in many environments, e.g. soils, hot springs and marine water [6,19,20,21,22]. Recent work has found that AOA can be up to 3000 times more abundant in soil than AOB [22,23,24], meaning that AOA are the most abundant ammonia oxidizing organisms in soil ecosystems [25]. The soil ammonia oxidizing community is known to be influenced by plant types and management, but different segments of this community respond differently [24,26,27]. For example, the abundance and composition of AOB is significantly altered by long-term fertilization, but AOA are rarely affected [24,27]. The nitrification activity in soil ecosystems is known to be correlated with the abundances and structures of ammonia oxidizers [24,28,29]. We therefore hypothesize that different biofuel cropping systems, especially those that rely on Nfertilization, will influence the composition of ammonium oxidizers in soil, with potential consequences for nitrification rates.

Denitrification, which reduces nitrate to N 2 gas, is carried out by a diverse group of microorganisms belonging to more than 60 genera of bacteria, archaea, and some eukaryotes [30]. Complete denitrification involves four steps: NO 3 2 RNO 2 2 RNOR N 2 ORN 2 . The enzyme nitrous oxide reductase (encoded by nosZ) that reduces N 2 O to N 2 is essential for complete conversion of NO 3 2 to N 2 . Approximately 17 Tg N is estimated to be lost from the land surface through denitrification every year [31]. It is known that the structure and activity of denitrifiers in the terrestrial ecosystem could be significantly influenced by the plant species [7,29,32]. In a study of a maize-cropped field, it was found that organic or mineral fertilizer applications could affect both the structure and activity of the denitrifying community in the long term, with changes persisting for at least 14 months [33]. The potential denitrification rate was found to be significantly correlated to the denitrifier density, as estimated by the quantification of nosZ gene copy numbers [34]. Denitrification releases mineralized nitrogen in the soil ecosystem to the atmosphere, and thus, the balance between denitrification and N-fixation, can determine the biologically available N for the biosphere (Arp, 2000).

It is known that plant species can change the soil microbial community [1]. However, while much previous work has examined the microbial community differences between the established crops [7,29,34,35,36], less is known about how microbial communities in the agricultural soils develop during the transition from one cropping system to another (e.g. annual row crops to perennial biofuel feedstocks). Thus, to improve our knowledge of the effects of bioenergy feedstock production on the complex N-cycling microbial communities of terrestrial ecosystems, we followed the changes in soil microbial communities during a two-year establishment period of maize, switchgrass, Miscanthus6giganteus, and mixed tallgrass prairie. We monitored the abundance of key genes for nitrogen fixation, ammonia oxidation and complete denitrification (nifH, bacterial/archaeal amoA and nosZ) as well as the structural changes of these N-cycling genes and bacterial/archaeal 16S rRNA genes using real-time PCR and barcoded pyrosequencing methods respectively.

Materials and Methods

Study site and sampling

The experiment was conducted at the Energy Biosciences Institute's Energy Farm located southwest of Urbana, Illinois, USA. Miscanthus (Miscanthus6giganteus, MG), switchgrass (Panicum virgatum, PV), maize (Zea mays, ZM) and restored tallgrass prairie (used as control, NP) were planted in the spring of 2008. Miscanthus was replanted in the spring of 2009 due to its poor growth in 2008. Each crop was planted in a randomized block design, with a 0.7-ha plot of each crop randomly positioned within four blocks (n = 4 for each crop). Samples were collected in April 2008, before planting of these crops, in order to characterize the background soil microbial communities. Bulk soil samples were collected at monthly intervals during the growing seasons (June-September 2008 and 2009). 10 soil cores (0-10 cm depth, 1.8 cm diameter) were collected from each plot and homogenized in an ethanol-sanitized, plastic bucket. About 60 g of the well-mixed soil was then subsampled into a 50 mL tube for each plot, and kept on ice until brought to the lab and stored in a 280uC freezer. In total, 112 soil samples were collected. Following standard agricultural practices, only maize was fertilized (17 g/m 2 in 2008, 20 g/m 2 in 2009) with a mixture of urea, ammonia and nitrate (28% UAN). Herbicides were applied in these crops except the restored tallgrass prairie: 1.56 l/ha of Roundup (only applied in 2008) and 4.13 l/ ha of Lumax for maize; 1.37 l/ha of 2,4-Dichlorophenoxyacetic acid for Miscanthus; 1.37 l/ha of 2,4-Dichlorophenoxyacetic acid for switchgrass in 2008.

DNA extraction and purification

Soil samples were freeze-dried overnight until completely dry and then manually homogenized with a sterile screwdriver. DNA was extracted from 0.3 g soil using the FastDNA SPIN Kit For Soil (MP Biomedicals) according to manufacturer's protocol. Extracted DNA was then purified using CTAB. DNA concentrations were determined using the Qubit quantification platform with Quant-iT TM dsDNA BR Assay Kit (Invitrogen). DNA was diluted to 10 ng/mL and stored in 280uC freezer for the following molecular applications.

Real-time PCR

The abundances of nifH, archaeal amoA, bacterial amoA and nosZ genes in all the soil samples were quantified using real-time PCR. Quantitative real-time PCR was performed according to the methods modified from previous studies: nifH (as a measure of Nfixing bacteria) used primers PolF (59-TGC GAY CCS AAR GCB GAC TC-39) and PolR (59-ATS GCC ATC ATY TCR CCG GA-39) [5]; archaeal amoA (as a measure of ammonia-oxidizing archaea) used primers Arch-amoAF (59-STA ATG GTC TGG CTT AGA CG-39) and Arch-amoAR (59-GCG GCC ATC CAT CTG TAT GT-39) [6]; bacterial amoA (as a measure of ammoniaoxidizing bacteria) used primers amoA-1F (59-GGG GTT TCT ACT GGT GGT-39) and amoA-2R (59-CCC CTC KGS AAA GCC TTC TTC-39) [4]; and nosZ (as a measure of denitrification bacteria) used primers nosZ-F (59-CGY TGT TCM TCG ACA GCC AG-39) [37] and nosZ 1622R (59-CGS ACC TTS TTG CCS TYG CG-39) [7]. Purified PCR products from a common DNA mixture (equal amounts of DNA from all samples collected in August of 2008 and 2009) were used to prepare sample-derived quantification standards as previously described [38]. The copy number of gene in each standard was calculated by DNA concentration (ng/mL, measured by Qubit) divided by the average molecular weight (estimated based on the barcoded-pyrosequencing results) of that gene. In comparison to using a clone (plasmid) as standard, this method avoids the difference of PCR amplification efficiency between standards and samples caused by the different sequence composition in the PCR templates (single sequence in a plasmid for the standard vs. mixture of thousands of sequences in a soil sample). The 10 mL reaction mixture contained 5 mL of 26Power SYBR Green Master Mix (Applied Biosystems), 0.5 mL of BSA (10 mg/mL, New England Biolabs), 0.4 mL of each primer (10 mM) and 5 ng of DNA template. Real-time amplification was performed in an ABI Prism 7700 Sequence Detector with MicroAmp Optical 384-Well Reaction Plate and Optical Adhesive Film (Applied Biosystems) using the following program: 94uC for 5 min; 40 cycles of 94uC for 45 s, 56uC for 1 min (54uC for nifH gene), 72uC for 1 min. A dissociation step was added at the end of the qPCR to assess amplification quality. The specificity of the PCR was further evaluated by running twenty randomly selected samples (for each gene) on a 1% (w/v) agarose gel. The corresponding real-time PCR efficiency for each of these genes was estimated based on a two-fold serial dilution of the common DNA mixture described above. The qPCR efficiency (E) was calculated according to the equation E = 10 [21/slope] [39]. Triplicate qPCR repetitions were performed for each of the gene for all the samples. The real-time PCR amplification efficiency of nifH, archaeal amoA, bacterial amoA and nosZ genes was 1.9060.01, 1.9060.06, 1.7660.01 and 1.8260.01, respectively. The R 2 of all these standard curves was higher than 0.99. The detection limit of this real-time PCR assay was 10 copies/mL. The copy numbers of these genes per gram of dry soil was calculated by the copy numbers of the gene per ng of DNA multiplied by the amount of DNA contained in each gram of dry soil. The quantities of these genes were corrected assuming a DNA extraction efficiency of 30% [40,41].

T-RFLP

The soil samples collected from four replicated (blocks) plots of the four crops prior to planting, and then August of establishment years 1 and 2 (2008 and 2009; 48 samples in total) were analyzed by terminal restriction fragment length polymorphism (T-RFLP). The nifH, archaeal amoA, bacterial amoA, nosZ gene were amplified from these samples with FAM-labeled (on forward primer) primers PolF/R, Arch-amoAF/R, amoA-1F/2R and nosZ-F/1622R (see Table S1). The 16S rRNA gene was amplified with 8F (59-FAM-AGAGTTTGATCMTGGCTCAG-39) and 1492R (59-GGTT-ACCTTGTTACGACTT-39) [42,43]. The PCR reaction mixture (25 ml) contained 5 ml GoTaq Flexi Buffer (56), 2 ml MgCl 2 (25 mM), 0.25 ml DNA Polymerase (5 U/ml, Promega, Madison, Wis.), 1.25 ml BSA (10 mg/ml), 1 ml forward primer (10 mM), 1 ml reverse primer (10 mM), 1.25 ml dNTP Mix (10 mM), 2 ml DNA template (10 ng/ml). The PCR reaction was performed in a thermo cycler (BioRad, Hercules, CA) using a 5-min heating step at 94uC followed by 30 cycles of denaturing at 94uC for 1 min, annealing at 60uC (54uC for nifH) for 45 s, and extension at 72uC for 1 min, with a final extension step of 5 min at 72uC. The PCR products were purified by QIAquick PCR Purification Kit (QIAGEN) and digested at 37uC overnight in a 20-ml mixture containing 2 ml NEB Buffer (106), 0.5 ml AluI/HhaI (20 U/ml) and 5 ml PCR product. 5 ml of the digested product was sent to the Core Sequencing Facility (University of Illinois at Urbana-Champaign) for fragment analysis. ROX1000 was used as inner standard. T-RFLP profiles were analyzed by GeneMarker (v 1.85) according to the manufacturer's instruction. Fragments with sizes between 50 bp and the length of the PCR products and peak area .500 were selected for T-RFLP profile statistical analysis. The profile data were normalized by calculating the relative abundance (percentage) of each fragment (individual peak area divided by the total peak area).

Barcoded pyrosequencing

The same samples used in T-RFLP were also used in barcoded pyrosequencing. The four replicated samples collected from the same crop at the same time were combined in to one composite sample for the construction of sequence libraries. Altogether, nine samples (one sample for background soil [mixed from the soils from all the plots before planting these crops], four samples from each of the different crops at Aug 2008, and four from Aug 2009) were obtained. Furthermore, all of these N-transforming genes of the sample collected from MG at the end of the second year (MG2) were sequenced twice in two different lanes to estimate the variation of the sequencing method. The 16S rRNA gene of ZM2 was also sequenced twice. Details of primers and PCR conditions used in the study are listed in Table S1.

The nifH, archaeal amoA, bacterial amoA, nosZ and 16S rRNA genes (V4-V5 region) were amplified using the barcode primers PolF/R, Arch-amoAF/R, amoA-1F/2R, nosZ-F/1622R and U519F/U926R, respectively (primers are shown in Table S1). The primers (HPLC purified, Integrated DNA Technologies) were designed as 59-Fusion Primer+barcode+gene specific primer-39 (Fusion Primer A, 59-CGTATCGCCTCCCTCGCGCCAT-CAG-39; Fusion Primer B, 59-CTATGCGCCTTGCCAGC-CCGCTCAG-39). The PCR conditions were optimized and primers with appropriate barcodes (10 bp) were selected. The barcodes used for each primer are described in NCBI SRA, with accession number SRA023700. The 50 mL PCR mixture contained 10 mL of 56 Phusion HF Buffer (Phusion GC Buffer was used for bacterial amoA gene amplification, both buffer contains 7.5 mM MgCl 2 ), 1 mL of 10 mM dNTPs, 2.5 mL of 10 mg/mL BSA, 0.5 mL of 2 U/mL Phusion Hot Start DNA Polymerase (FINNZYMES), 4 mL of 10 mM forward/reverse primer mixture and 4 mL of 10 ng/mL DNA templates. 1 mL of 100% DMSO was supplemented into the PCR mixture in bacteria amoA gene (GC rich) amplification. The PCR amplification was performed in a thermal cycler (BioRad, Hercules, CA) using the program 98uC for 2 min; 30 cycles of 98uC for 10 s, 60uC (54uC for nifH gene and 56uC for bacterial amoA gene) for 30 s, 72uC for 20 s; 72uC for 5 min. The PCR product was first checked on a 1.2% w/v agarose gel, and then purified by QIAquick PCR Purification Kit (QIAGEN). The DNA concentration of the purified PCR product was measured by Qubit Fluorometer using the Quant-iT TM dsDNA BR Assay Kit (Invitrogen) according to the manual. PCR products of the same gene, to be run together in the same lane (1/16 plate) in 454 sequencing, were mixed in equal mole amounts and run on a 2% w/v agarose gel. The target bands were cut from the gel and purified by QIAquick Gel Extraction Kit (QIAGEN). The DNA concentrations of the purified PCR products were measured by Qubit Fluorometer and adjusted to 50 nM. The nifH, archaeal amoA, bacterial amoA and nosZ genes PCR products were then mixed in equal mole amounts and sequenced on a Genome Sequencer FLX Instrument (Roche) using GS FLX Titanium series reagents. The 16S rRNA gene was run in a separate lane.

Sequence analysis

Sequences were first extracted from the raw data according the Genome Sequencer Data Analysis Software Manual (Software Version 2.0.00, October 2008) by the sequencing center (Roy J. Carver Biotechnology Center, University of Illinois at Urbana-Champaign). The sequences with low quality (length ,50 bp, which ambiguous base 'N', and average base quality score ,20, for detail see manual) were removed. The sequences that fully matched with the barcodes were selected and distributed to separate files for each of the different genes, after removal of the barcode, using RDP Pipeline Initial Process (http://pyro.cme. msu.edu/). For each gene, the sequences that didn't match with the gene specific primers or had a read length shorter than 350 bp were removed. The sequences that matched with the reverse primer were converted to their reverse complement counterparts using BioEdit to make all the sequences forward-oriented.

The 16S rRNA gene sequences were aligned by NAST (Greengenes). The sequences with significant matched minimum length ,300 and identity ,75% were removed. The aligned 16S rRNA gene sequences were used for chimera check using Bellerophon method in Mothur [44]. Distance matrices were calculated by ARB using the neighbor joining method [45]. A lane mask was used in calculating the 16S rRNA gene sequences to filter out the hyper variable regions. Operational Taxonomic Units (OTUs) were then classified using a 97% nucleotide sequence similarity cutoff and rarefaction curves were constructed based on the distance matrices (both of nucleotide and amino acid sequences) using DOTUR [46]. The phylogenetic affiliation of each 16S rRNA gene sequence was analyzed by RDP CLASSI-FIER (http://rdp.cme.msu.edu/) using confidence level of 80%.

The 16S rRNA gene sequences were also processed by QIIME pipeline and denoised by Denoiser V0.91 according to the manual [47,48]. The results were compared to that obtained by RDP pipeline. In total, 26,431 valid sequences were obtained after denoising using QIIME, which is 12.2% less than that obtained by RDP pipeline (without denoising). Using the 97% similarity cutoff, 8,568 OTUs were obtained, which is 4.7% lower than that observed by RDP pipeline. After random re-sampling to the same sequence depth (1789 sequences per sample) using Daisy_chopper (http://www.genomics.ceh.ac.uk/GeneSwytch/Tools.html), the number of OTUs for each sample obtained by two different processing methods (QIIME, denoised and RDP, non-denoised) was compared (Fig. S1). The estimated number of OTUs after denoising was similar to that obtained by RDP pyrosequencing pipeline (without denoising), showing that the denoising process had a very limited influence on our diversity analysis. The data reported in this paper was analyzed using RPD pipeline described in the previous paragraph.

The nifH, archaeal amoA, bacterial amoA and nosZ genes sequences were blasted against a non-redundant protein sequence database (download from NCBI) using BLASTX with an E-value cutoff of 0.001. The top 10 closest matches of each sequence were estimated using a custom made Perl script to remove possible chimeras and sequences with sequencing errors causing frameshifts. Sequences with different regions matching the same sequence in the database but with different frame positions were considered to be frameshifts. Sequences that matched two or more different origin sequences were classified as chimeras. The nucleotide sequences of nifH, archaeal amoA, bacterial amoA and nosZ genes were translated into amino acid by Geneious (http://www.geneious.com/) based on the frame positions obtained from BLASTX. The redundant sequences (identical sequences) were removed using CD-Hit [49], and the representatives with longest length were selected for following phylogenetic analysis. Both of the nucleotide and amino acid sequences of these N-cycling genes were aligned by MUSCLE 3.7 [50] using program default settings. Operational Taxonomic Units (OTUs) were then classified and rarefaction curves were constructed based on the distance matrices (both of nucleotide and amino acid sequences) using DOTUR [46]. Previous studies showed that the amino acid sequences of AmoA and NosZ similarity around 90% is generally relevant to 97% similarity of 16S rRNA gene [51,52]. Thus, all these N-cycling gene sequences were classified into OTUs using a 90% amino acid sequence similarity cutoff, and phylogenetic trees were built in ARB using the neighbor-joining method. Sequences of all the samples and genes were also randomly resampled to identical sequencing depth (the smallest sequencing effort) using Daisy_chopper (http://www.genomics.ceh.ac.uk/Gen-eSwytch/Tools.html) to avoid the potential bias caused by sequencing effort difference [53].

The 454-pyrosequencing data were deposited in NCBI SRA under accession number SRA023700.

Statistical analysis

ANOVA combined with post hoc Tukey B test was used to estimate the difference of archaeal/bacterial amoA, nifH and nosZ genes abundances under different crops based on the quantitative PCR results from the replicated plots. The T-RFLP data from the replicated plots were used to follow the structural changes of soil microbial communities by plant types, and significance tests for these changes were conducted using Analysis of Similarity (ANOSIM) based on Bray-Curtis similarity coefficients. Correspondence analysis (CA) and Canonical correspondence analysis (CCA) were also used to visualize the predominant microbial community changes of archaeal/bacterial amoA, nifH, nosZ and 16S rRNA genes after planting bioenergy crops based on the T-RFLP data. These statistical analyses were done using the free software PAST (http://folk.uio.no/ohammer/past/). Based on our extensive pyrosequencing library, the OTUs/genera that showed monotonic (i.e. continuously increasing or decreasing) trends for each crop treatment over the two year establishment were presumed to be particularly noteworthy in terms of crop impact. The populations with continuously increased or decreased abundance in the two-year period after planting these bioenergy crops were selected using a custom Perl script.

Results

Quantification of nifH, archaeal amoA, bacterial amoA and nosZ genes Quantities of AOA in all of crops fluctuated over the two growing seasons in a similar pattern (Fig. 1), but the abundance of this gene was always higher in MG than NP and PV. The quantity of bacterial amoA genes in ZM significantly increased from 1.4760.61610 8 to 3.2660.94610 8 during the first three months of establishment and thereafter remained higher in ZM than in the other cropping systems. The nitrification rates under these crops were analyzed in the second year by estimating the accumulation of the nitrate in buried soil bags, and linear regression revealed that the nitrification rates were significantly related to the quantities (log) of archaeal amoA genes (R 2 = 0.61, P = 0.03, n = 12), but not to the quantities of bacterial amoA genes (Fig. 2). The abundance of nifH genes remained stable for all the crops, ranging from 7610 7 to 9610 7 copies per gram of dry soil (Fig. 1). No significant differences were observed among the different crops in the first year. In the second year, the population sizes of diazotrophs in PV and NP had significantly increased in comparison to the first year (P = 0.0001 and 0.0002). The population size of denitrifiers was less variable in comparison to the other N-cycling populations; however, the copy number of nosZ increased in ZM during the second year of the study and remained higher than in MG (Fig. 1).

Structural changes of N-cycling genes and microbial communities after planting of bioenergy crops

The community structural differentiation of nifH, archaeal amoA, bacterial amoA, nosZ and 16S rRNA genes under different bioenergy crops were analyzed by T-RFLP. These analyses used the fully replicated sample set from the randomized block design.  Correspondence analysis (Fig. 3) showed that the soil microbial communities in the initial plots did not show any relationship with the crop treatments being applied. During the establishment of these bioenergy crops, the community composition of denitrification bacteria (nosZ) under ZM was completely separated (ANOSIM, P,0.05, Table S2) from those under the other crops by the end of the second year along the first-axis, which explained 51.7% of the total variance. None of the other groups showed significant clustering by plant, although the community composition of AOB (bacteria amoA) under ZM appeared to be separated from that of MG (ANOSIM, P = 0.17) along the second axis, which explained 14.2% of the total variance. In addition to plant species, the changes of soil microbial communities also could be caused by the variation of environmental conditions. To compare the magnitude of the changes of soil microbial communities related only to plant species, a constrained ordination method was also used. Canonical Correspondence analysis (CCA, Fig. S2) revealed that, at the end of the second year, the microbial communities under ZM were most different from the three cropping systems for bacterial amoA, nosZ, and 16S rRNA. Archaeal amoA was most distinct under MG, followed by ZM, and nifH was equally separated under all cropping systems (Fig. S2).

Diversity of nifH, archaeal amoA, bacterial amoA, nosZ and 16S rRNA genes To further understand the composition of microbial community in the field, the nifH, archaeal amoA, bacterial amoA, nosZ and 16S rRNA genes were deeply sequenced using the pyrosequencing approach. In total, 143,487 reads were obtained for these genes. The numbers and qualities of these sequences are described in Table 1, Table S3 and Text S1. The reproducibility of the pyrosequencing result was estimated by comparing the observed microbial composition between repeat sequencing runs for all these genes (Fig. 4). Linear regression analysis indicated a high reproducibility (R 2 = 0.95) of our pyrosequencing.

High diversity of nifH, archaeal amoA, bacterial amoA and nosZ genes were observed with 10899, 3187, 3945 and 11242 unique nucleotide sequences and 2286, 2246, 3633 and 4208 unique deduced amino acid sequences respectively (Fig. S3). These sequences were then translated to amino acid sequence according to the BLASTX report. The amino acid sequences of nifH, archaeal amoA, bacterial amoA and nosZ genes were classified into 229, 309, 330 and 331 OTUs, respectively, with a similarity cutoff of 90%. After random re-sampling to the same sequencing depth (697 sequences for each sample), the adjusted total number of OTUs for these genes were 217, 303, 319 and 278, respectively (Table S4). Rarefaction analysis of these genes showed that the diversity of archaeal amoA gene in ZM2 (second year ZM) and MG1 (first year MG) was markedly lower than the others (Fig.  S4). The diversity of nifH and nosZ genes slightly decreased in ZM2.

The diversity of bacterial and archaeal 16S rRNA genes was much higher than these N-cycling genes. In total, 19,824 unique 16S rRNA gene sequences and 8,989 species (OTUs classified at 97% similarity cutoff) were observed. RDP classification showed that these sequences covered 16 bacterial and 1 archaeal phyla, including 201 genera (Fig. S5). Proteobacteria and Acidobacteria were the most predominant phyla in the soil (.20%). The sequences belonging to Proteobacteria were distributed over 86 different genera, while 94.1% of the sequences in Acidobacteria belonged to Family Gemmatimonadaceae, with GP1 as the most predominant genus (accounted 35.9% of the Acidobacteria).

To understand which phylotypes were impacted by vegetation type, the OTUs of nifH, archaeal amoA, bacterial amoA, nosZ and 16S rRNA (genus for 16S rRNA) genes that continuously increased or decreased over the two-year establishment of these bioenergy crops were identified (Fig. 5 and Fig. S6, S7, S8, S9). After planting of these bioenergy crops, 27.5%, 15.4%, 22.7% and 14.5% of the total archaeal amoA, bacterial amoA, nifH, and nosZ phylotypes, respectively, were found to be continuously increasing or decreasing ( Table 2). Details of these continuously changed Ncycling OTUs are described in Text S1 and Figure S6, S7, S8, S9.

Pyrosequencing of 16S rRNA gene revealed 19.9% of the bacterial genera (39), spanning six phyla, continuously changed after planting of these bioenergy crops (Fig. 5). Only genus Methylibium was changed in all the crops, with decreased abundance in MG and increased abundance in the other crops. Rhodanobacter only appeared after planting of ZM (7 sequences for both of ZM1 and ZM2), and it was undetectable either in the background soil or in the soil under other crops. Consistent with the changes of Nitrosospira-like bacterial amoA OTU (see above), the abundance of genus Nitrosospira in the 16S rRNA library also increased in ZM. The abundance of genus Nitrospira, which is known as a nitrite-oxidizing bacteria, also increased in ZM. Most of the changed genera in MG decreased or even disappeared, except Terrabacter and Herbaspirillum. All of them were found at low abundance (,1%). Although many genera in Proteobacteria were changed, the total abundance of this most predominant phylum was quite stable under all of the crops (Fig. S5).

Discussion

In this study, we monitored the structural and quantitative changes of the key genes involved in N-cycling as well as the overall bacterial/archaeal community during two-year establishment of four different bioenergy feedstock crops, and analyzed the shifts of specific soil microbial populations in response to different types of crops. We were able to detect significant changes in the abundance of many of these microbial functional groups within 2 years of initial crop establishment. We also found that traditional row-crop agriculture of maize has a larger Valid sequences of nifH, archaeal amoA, bacterial amoA and nosZ genes were defined as high quality sequences with correct barcode and primer (at 59-end), length .350 bp and that did not have frameshifts and chimeric structure. The possible sequencing errors causing frameshifts and chimeras were removed based on the BLASTX result. Sequences with different regions matching the same sequence in the database but with different frame positions were considered to be frameshifts. Sequences that matched two or more different origin sequences were classified as chimeras. Valid sequences of 16S rRNA gene were sequences with correct barcode and primer, length .350 bp and passed the chimeric check program in Greengenes with the Bellerophon method. The sequence numbers for each sample are listed in Table S3. doi:10.1371/journal.pone.0024750.t001 impact on the soil N-cycling community than any of the perennial bioenergy feedstock crops (Figs. 1, S2), while the perennial crops were associated with overall community shifts in the phyla Planctomyces, Firmicutes, and Actinobacteria (Fig. 5).

Traditional maize cultivation significantly increased the total abundance of ammonia-oxidizing bacteria and denitrifying bacteria (Fig. 1), altered the community composition of denitrifying bacteria (Fig. 3) and decreased the diversity of ammonia-oxidizing archaea, denitrifying bacteria, and diazotrophs (Fig. S4). This may be due to the application of Nfertilizer, which occurred only in ZM plots. Ammonia oxidizers are sensitive to N-fertilizer [24,27], and these responses were manifested in the increased population size of AOB and the high number of markedly changed AOA species. The nitrification rate was significantly correlated with the quantity of archaeal amoA, but not bacterial amoA, indicating AOA was the major ammonia-oxidizer. Deep understanding of the structural shifts of key functional genes can help us to better understand changes in microbial activity in the environment. From the present database, we know that the global diversity of the nifH, archaeal amoA, bacterial amoA and nosZ genes, as well as the other functional genes of microorganisms, is high. The traditional approaches (e.g. clone library, DGGE and T-RFLP) used in previous studies may largely underestimate the diversity of microbial communities involved in soil nitrogen cycling. Mounier et al. (2004) revealed that even a large library with 713 clones was insufficient to enumerate the diversity nosZ gene in maize rhizosphere, showing the high complexity of N-cycling genes [54]. Thus, high-  OTUs of N-cycling genes were classified based on a cutoff of 90% amino acid sequences similarity. 16S rRNA genes were classified into genus level by RDP Classifier. Details of the abundance changes of the OTUs or genera are shown in Fig. 5  throughput deep sequencing approaches are essential to improve our knowledge of the diversity of these functional genes. In the present study, using barcoded 454-pyrosequencing approach, we found high diversity (ranging from ,3100 to ,11200 unique nucleotide sequences) of nifH, archaeal amoA, bacterial amoA and nosZ genes in the soil ecosystem, which far surpasses the diversity of the N-cycling genes observed in previous studies [54,55,56,57,58,59,60,61]. The rarefaction curves of these genes were close to saturation after sequencing ,1000 for each sample, indicating that such a sequencing effort is sufficient to elucidate the diversity and structure of the complex soil N-cycling communities. The high similarity between repeat runs of these genes (Fig. 4) demonstrates the high reproducibility and reliability of this barcoded pyrosequencing method. This result also indicates that the variation of pyrosequencing, resulting from random sampling of gene targets during emulsion PCR [62], can be greatly reduced by increasing the sequencing depth and library coverage. During the establishment of these bioenergy crops, about 15%-30% of N-cycling genes and the detected bacterial/ archaeal genera were continuously changed, indicating that a large proportion of soil microbes were affected by the transition to different bioenergy feedstock systems. Most of these phylotypes changed uniquely in one of the crops, indicating that the changes were mainly caused by the particular experimental crop treatment (specific plant species or management, such as fertilization) and not due to the environmental conditions that fluctuated in all treatments (e.g. temperature and moisture). Contrary changes of certain populations in different crops also support this conclusion; for example, the abundance of Methylibium (belonging to b-Proteobacteria) decreased in MG but increased in the other crops. The abundance of Bacteroidetes was previously found to be lower in the soil of Miscanthus-dominated grasslands (4%) in comparison to forest soil (6%) [63]. We found that the decrease of Bacteroidetes in MG was mainly due to the disappearance of the genus Ferruginibacter ( Fig. 5 and S5).

The structure of nifH gene was completely separated according to vegetation by the end of the second year (Fig. S2), which suggests that the structure of N-fixation bacterial population was particularly sensitive to plant genotype. Tan et al. (2003) has revealed that the structure of diazotrophs was not only different among rice species, but also changed rapidly with fertilization. The diversity of the nifH gene was obviously reduced within 15 days after fertilization [17]. Thus, the decreased diversity of nifH in ZM may be also due to the application of fertilizer not the presence of maize. However, the population size of N-fixing bacteria did not change under ZM, indicating N-fertilization may not change the quantity of soil diazotroph [64]. The N-fixing activity was expected to increase in MG [14,15]. Although the population size of total free-living soil N-fixing bacteria was not significantly increased by growth of MG in the two year period, the abundance of genus Herbaspirillum increased. Herbaspirillum species are known as endophytic diazotrophs that are enriched by C4-prennial grasses including Miscanthus [65,66,67]. Thus, we speculate that the abundance increase of Herbaspirillum in the bulk soil was likely due to the root exudates (e.g. organic carbons) released by Miscanthus, which favored the growth of this population. Our results also suggest that Miscanthus may only selectively enhance the activity of specific diazotrophs, not the whole N-fixing microbial community.

The AOA are thought to be more stable and less responsive to environmental differences than AOB, as revealed by previous quantitative studies [6,68]. However, in the present study we found that, while the population size of AOA was relatively stable, the structure of the AOA community was sensitive to the different cropping systems. The diversity of AOA markedly decreased after planting of maize, with 41 of the AOA phylotypes disappearing (Fig. S6). In contrast, the population size of AOB significantly increased after planting of maize in both the qPCR of bacterial amoA and the 16S rRNA pyrosequencing results. In addition to the increase of genus Nitrosospira (AOB) [69], the abundance of Nitrospira (nitrite-oxidizing bacteria) [70] also increased in N-fertilized maize (Fig. 5). However, the number of changed bacterial amoA phylotypes in ZM was much less than the other crops. Therefore, these two different groups of ammonia-oxidizers respond to the N-fertilization in a very different way. The population size of AOB increased immediately in the first growing season, thus, we hypothesize that the increased AOB abundance in maize was likely due not to the growth of maize plants, but to the application of N-fertilizer, which increased the ammonia content in the soil [24,27]. It has been found that AOB population size increased in seven days after applying of N-fertilizer, and it was still significantly greater than unfertilized soil 8 months after the last application of ammonia [41]. Consistently, we found that the population size of AOB doubled in less than three months and maintained a relatively high level over the two-year study even though measurable NH 4 + in the soil declined over this time period to levels close to that of the unfertilized plots ( Fig. 1 and S10). These results indicate the AOB population size can be quickly increased by N-fertilization and can remain for a long period even after the measurable ammonia has been consumed.

The nitrification rate was found to be significantly correlated with the quantity of archaeal amoA rather than bacterial amoA, indicating AOA rather than AOB may be the major active ammonia-oxidizer in these soils. Contradictory conclusions on the relative importance of AOA and AOB in soil nitrification have been previously reported where nitrification was found to be associated with the changes of archaeal amoA abundance or higher archaeal transcriptional activity in some of the soils [58,71,72]. In contrast, the nitrification kinetics in the other soils were correlated with the growth of AOB [73,74]. It has been found that the ammonia affinity of ''Candidatus Nitrosopumilus maritimus'' (a marine AOA) is much higher than AOBs, and its growth may be enhanced by relatively low ammonia concentration [75]. Thus, the contradictory conclusions from these studies may be due to the different soil ammonia concentration used in these experiments [58,73]. These results hint that AOB may be more active in soils amended with ammonia, while AOA are more active in soils with low ammonia concentration [58]. The nitrification rate outlier (ZM Aug; Fig. 2) had highest population size of AOB, suggesting that the activity of AOB was enhanced by N-fertilization and also supported the above speculation. In support of this speculation, a recent publication shows that recovery of nitrification potential after disruption was dominated by AOB in cropped soils while AOA were responsible RNP in pasture soils [76].

It is known that nitrogen fertilization can change the structure and activity of denitrifying community, and subsequently affect the N 2 O emission [33,34,77,78]. Large amounts (1.3%) of the applied N-fertilizer in maize fields (north Colorado) are converted to N 2 O by the combination of nitrification and denitrification [79]. However, it is still unclear how fertilization changes the microbial community, since most of the previous studies are based on the already established fields [7,29,33]. Our study revealed that the structure of denitrifying bacteria in maize soil was significantly differentiated from the other crops at the second year (Fig. 3).

However, the population size of denitrifiers was relatively stable in all the crops in comparison to other N-cycling microbial communities, which only slightly increased in maize at the second year. The high stability of denitrifying population abundance could be explained by the high diversity and functional redundancy of denitrification community [30,80].

In conclusion, our two-year study of transitional agriculture shows that specific N-transforming microbial communities develop in the soil in response to different bioenergy crops. Each N-cycling microbial group responded in a different way after planting with different bioenergy crops. In general, planting of maize has a larger impact on the soil N-cycling community than the other bioenergy crops. Our results also indicate that application of Nfertilizer may not only cause short-term environmental problems, e.g. water contamination, but also can have long-term influence on the global biogeochemical cycles through changing the soil microbial community structure and abundance. Since soil types and other environmental factors may also impact the N-cycling microbial community, the universality of our findings needs to be confirmed by additional study at different sites. Figure S1 Number of OTUs obtained by two different processing methods: QIIME (denoised) and RDP pyrosequencing pipeline (non-denoised). (TIF) Figure S2 Structural changes of archaeal amoA, bacterial amoA, nifH, nosZ and 16S rRNA genes after planting Miscanthus6giganteus (MG), Panicum virgatum (PV), restored prairie (NP) and Zea mays (ZM) revealed by T-RFLP and Canonical correspondence analysis (CCA). The number on each axis shows the explained total variation. The soil samples were collected from four replicated plots for each plant at each time point. * Correspondence analysis was used for the samples collected before planting bioenergy crops. (TIF) Figure S3 OTU classification of valid sequences at different distance levels based on nucleotide and deduced amino acid sequences. (TIF) Figure S4 Rarefaction analysis of the diversities of nifH, archaeal amoA, bacterial amoA, nosZ and 16S rRNA genes in the soil underneath different bioenergy crops. The OTUs of nifH, archaeal amoA, bacterial amoA and nosZ genes were classified at 90% similarity cutoff based on amino acid sequences, and 16S rRNA gene was classified at 97% similarity cutoff on nucleotide sequences. BG0 represents the samples collected before planting bioenergy crops. MG, PV, NP, and ZM represent Miscanthus6giganteus, Panicum virgatum, restored prairie and Zea mays respectively. 1 and 2 represent samples collected in the first and second growing seasons.    Text S1

Supporting Information

Figure 1 .
1Changes in abundance of nifH, archaeal amoA, bacterial amoA and nosZ genes in plots after planting Miscanthus6giganteus (MG), Panicum virgatum (PV), restored prairie (NP) and Zea mays (ZM). The copy number of genes in each gram of dry soil was estimated based on the results of real-time PCR (copy number in each ng DNA) and the average amount of extracted DNA (6.23 mg per dry soil) and assuming DNA extraction efficiency was 30%[40]. The R 2 of the standard curve of all these genes was higher than 0.99. The real-time PCR amplification efficiency of nifH, archaeal amoA, bacterial amoA and nosZ genes was 1.9060.01, 1.9060.06, 1.7660.01 and 1.8260.01 respectively. *Represents values that are significantly different (P,0.01). doi:10.1371/journal.pone.0024750.g001

Figure 2 .
2Relationship between the concentration of ammonia-oxidizing archaea/bacteria and nitrification rate. The nitrification rate was determined over the same time as our sample collection in 2009. Nitrification rate was calculated based on the accumulation of nitrate in soil bags incubated in the field (0-10 cm depth) for 15 to 32 days. doi:10.1371/journal.pone.0024750.g002

Figure 3 .
3Structural changes of archaeal amoA, bacterial amoA, nifH, nosZ and 16S rRNA genes after planting Miscanthus6giganteus (MG), Panicum virgatum (PV), restored prairie (NP) and Zea mays (ZM) revealed by T-RFLP and Correspondence analysis (CA). The number on each axis shows the percentage of total variation explained. The soil samples were collected from four replicate plots for each plant at each time point. doi:10.1371/journal.pone.0024750.g003

Figure 4 .
4Reproducibility of the pyrosequencing replicates. OTUs of nifH, archaeal amoA, bacterial amoA and nosZ genes were classified at a nucleotide similarity cutoff 90%. The 16S rRNA gene sequences were classified to genus level by RDP classifier. One of the samples was duplicated for each gene. doi:10.1371/journal.pone.0024750.g004

Figure 5 .
5Microbial genera that changed after planting of Miscanthus6giganteus (MG), Panicum virgatum (PV), restored prairie (NP) and Zea mays (ZM). Sequences were classified by RDP Classifier project. 2/+ represents the genus continuously decreased/increased after planting the crops; 22/++ represents the genus disappeared/appeared after planting the crops. *sequences belonging to Crenarchaeota, which could not be classified to genus level. doi:10.1371/journal.pone.0024750.g005

Figure S5
S5Phylum level microbial community composition in the soil under different plants before and for two years after transition to bioenergy cropping. * represent significantly changed phylum. MG, Miscanthus6giganteus; PV, Panicum virgatum; NP, restored prairie; ZM, Zea mays. (TIF) Figure S6 (a) Phylogenetic tree of and (b) abundance of archaeal amoA OTUs that continuously changed after planting Miscanthus6 giganteus (MG), Panicum virgatum (PV), restored prairie (NP) and Zea mays (ZM). OTUs were classified based on a cutoff of 90% amino acid sequence similarity. (TIF) Figure S7 (a) Phylogenetic tree of and (b) abundance of bacterial amoA OTUs that continuously changed after planting Miscanthus6 giganteus (MG), Panicum virgatum (PV), restored prairie (NP) and Zea mays (ZM). OTUs were classified based on a cutoff of 90% amino acid sequence similarity. (TIF) Figure S8 (a) Phylogenetic tree of and (b) abundance of nifH OTUs that continuously changed after planting Miscanthus6giganteus (MG), Panicum virgatum (PV), restored prairie (NP) and Zea mays (ZM). OTUs were classified based on a cutoff of 90% amino acid sequence similarity. (TIF) Figure S9 (a) Phylogenetic tree of and (b) abundance of nosZ OTUs that continuously changed after planting Miscanthus6giganteus (MG), Panicum virgatum (PV), restored prairie (NP) and Zea mays (ZM). OTUs of were classified based on a cutoff of 90% amino acid sequence similarity. (TIF) Figure S10 Nitrate and ammonia concentration in bulk soil. Soil samples for chemical and microbiological analysis were collected in the same week for each time point, except Sep 2008 when the nitrate and ammonia concentrations were not measured. MG, Miscanthus6giganteus; PV, Panicum virgatum; NP, restored prairie; ZM, Zea mays. These data were measured over the same period of our sample collection by the Biogeochemistry laboratory (C. Smith and M. David). (TIF) Figure S11 Pyrosequencing read length based on the raw sequence reads. (TIF)

Table 1 .
1Quality of barcoded pyrosequencing reads.Number of sequences

Genes
Correct barcode and primer
Length .350 bp
a Valid
Each sample (range)

nifH
28,334
27,781
21,111
1,312-2,956

Archaeal amoA
16,978
16,226
14,025
697-1,792

Bacterial amoA
28,254
27,874
21,817
1,726-2,569

nosZ
33,838
28,819
22,590
1,600-2,951

16S rRNA
30,487
30,175
30,101
2,034-3,488

Total
137,891 (96.1%)
130,875 (91.2%)
109,644 (76.4%)
697-3,488

Total number of raw reads was 143,487.

a

Table 2 .
2Number of OTUs or genera that changed continuously after planting Miscanthus6giganteus (MG), Panicum virgatum (PV), restored prairie (NP) and Zea mays (ZM).Genes
MG
PV
ZM
NP
*Total

Archaeal amoA
4
1 6
6 1
2 3
8 5

Bacterial amoA
18
18
7
20
51

nifH
8
2 1
1 4
2 3
5 2

nosZ
12
18
24
16
48

16S rRNA
17
15
14
21
40

Table S1
S1Primers and annealing temperature for nifH, archaeal amoA, bacterial amoA, nosZ and 16S rRNA genes. (DOC)Table S2 Comparsion of the microbial community structures between different crops. (DOC)

Table S3
S3Number of valid sequences for each gene in each sample. (DOCX)

Table S4
S4Number of OTUs observed after random re-sampling to the identical sequencing depth (697 sequences/sample. (DOCX)
PLoS ONE | www.plosone.org
September 2011 | Volume 6 | Issue 9 | e24750
AcknowledgmentsThe authors gratefully acknowledge help from A. Duong, D. Fieckert, A. Groll, N. Peld, and M. Masters with field sampling and laboratory work. C. Smith and M. David kindly provided denitrification rate estimates for sampled plots. P.Y. Hong provided helpful comments on this manuscript.Author Contributions
Microbial diversity in soil: selection microbial populations by plant and soil type and implications for disease suppressiveness. P Garbeva, J A Van Veen, J D Van Elsas, Annu Rev Phytopathol. 42Garbeva P, van Veen JA, van Elsas JD (2004) Microbial diversity in soil: selection microbial populations by plant and soil type and implications for disease suppressiveness. Annu Rev Phytopathol 42: 243-270.

Effects of above-ground plant species composition and diversity on the diversity of soil-borne microorganisms. G A Kowalchuk, D S Buma, W De Boer, Pgl Klinkhamer, J A Van Veen, Antonie Van Leeuwenhoek International Journal of General and Molecular Microbiology. 81Kowalchuk GA, Buma DS, de Boer W, Klinkhamer PGL, van Veen JA (2002) Effects of above-ground plant species composition and diversity on the diversity of soil-borne microorganisms. Antonie Van Leeuwenhoek International Journal of General and Molecular Microbiology 81: 509-520.

Ecological linkages between aboveground and belowground biota. D A Wardle, R D Bardgett, J N Klironomos, H Setala, W H Van Der Putten, Science. 304Wardle DA, Bardgett RD, Klironomos JN, Setala H, van der Putten WH, et al. (2004) Ecological linkages between aboveground and belowground biota. Science 304: 1629-1633.

The ammonia monooxygenase structural gene amoA as a functional marker: molecular fine-scale analysis of natural ammonia-oxidizing populations. J H Rotthauwe, K P Witzel, W Liesack, Appl Environ Microbiol. 63Rotthauwe JH, Witzel KP, Liesack W (1997) The ammonia monooxygenase structural gene amoA as a functional marker: molecular fine-scale analysis of natural ammonia-oxidizing populations. Appl Environ Microbiol 63: 4704-4712.

Comparison of nifH gene pools in soils and soil microenvironments with contrasting properties. F Poly, L Ranjard, S Nazaret, F Gourbiere, L J Monrozier, Appl Environ Microbiol. 67Poly F, Ranjard L, Nazaret S, Gourbiere F, Monrozier LJ (2001) Comparison of nifH gene pools in soils and soil microenvironments with contrasting properties. Appl Environ Microbiol 67: 2255-2262.

Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. C A Francis, K J Roberts, J M Beman, A E Santoro, B B Oakley, Proc Natl Acad Sci U S A. 102Francis CA, Roberts KJ, Beman JM, Santoro AE, Oakley BB (2005) Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. Proc Natl Acad Sci U S A 102: 14683-14688.

Structure and function of denitrifying and nitrifying bacterial communities in relation to the plant species in a constructed wetland. O Ruiz-Rueda, S Hallin, L Baneras, FEMS Microbiol Ecol. 67Ruiz-Rueda O, Hallin S, Baneras L (2009) Structure and function of denitrifying and nitrifying bacterial communities in relation to the plant species in a constructed wetland. FEMS Microbiol Ecol 67: 308-319.

Determinants of the distribution of nitrogen-cycling microbial communities at the landscape scale. D Bru, A Ramette, N P Saby, S Dequiedt, L Ranjard, ISME J. 5Bru D, Ramette A, Saby NP, Dequiedt S, Ranjard L, et al. (2011) Determinants of the distribution of nitrogen-cycling microbial communities at the landscape scale. ISME J 5: 532-542.

Diversity and activity of free-living nitrogen-fixing bacteria and total bacteria in organic and conventionally managed soils. C H Orr, A James, C Leifert, J M Cooper, S P Cummings, Appl Environ Microbiol. 77Orr CH, James A, Leifert C, Cooper JM, 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.

Nitrogenase gene diversity and microbial community structure: a cross-system comparison. J P Zehr, B D Jenkins, S M Short, G F Steward, Environ Microbiol. 5Zehr JP, Jenkins BD, Short SM, Steward GF (2003) Nitrogenase gene diversity and microbial community structure: a cross-system comparison. Environ Microbiol 5: 539-554.

Global patterns of terrestrial biological nitrogen (N-2) fixation in natural ecosystems. C C Cleveland, A R Townsend, D S Schimel, H Fisher, R W Howarth, Global Biogeochemical Cycles. 13Cleveland CC, Townsend AR, Schimel DS, Fisher H, Howarth RW, et al. (1999) Global patterns of terrestrial biological nitrogen (N-2) fixation in natural ecosystems. Global Biogeochemical Cycles 13: 623-645.

Evidence for the functional significance of diazotroph community structure in soil. S F Hsu, D H Buckley, Isme Journal. 3Hsu SF, Buckley DH (2009) Evidence for the functional significance of diazotroph community structure in soil. Isme Journal 3: 124-136.

Non-symbiotic nitrogen fixation in forest ecosystems. Y Son, Ecological Research. 16Son Y (2001) Non-symbiotic nitrogen fixation in forest ecosystems. Ecological Research 16: 183-196.

The effect of fertilization on yield and quality of Miscanthus sinensis. H Schwarz, P Liebhard, K Ehrendorfer, P Ruckenbauer, Giganteus'. Ind Crop Prod. 2Schwarz H, Liebhard P, Ehrendorfer K, Ruckenbauer P (1994) The effect of fertilization on yield and quality of Miscanthus sinensis 'Giganteus'. Ind Crop Prod 2: 153-159.

Comparative biogeochemical cycles of bioenergy crops reveal nitrogen-fixation and low greenhouse gas emissions in a Miscanthus6giganteus agro-ecosystem ecosystems. S Davis, W Parton, F Dohleman, C Smith, S Grosso, DOI10.1007/s10021-009-9306-9Davis S, Parton W, Dohleman F, Smith C, Grosso S, et al. (2010) Comparative biogeochemical cycles of bioenergy crops reveal nitrogen-fixation and low greenhouse gas emissions in a Miscanthus6giganteus agro-ecosystem ecosystems DOI 10.1007/s10021-009-9306-9.

Polymerase chain reaction denaturing gradient gel electrophoresis analysis of the N 2 -fixing bacterial diversity in soil under Acacia tortilis ssp. raddiana and Balanites aegyptiaca in the dryland part of Senegal. Demba Diallo, M Willems, A Vloemans, N Cousin, S Vandekerckhove, T T , Environ Microbiol. 6Demba Diallo M, Willems A, Vloemans N, Cousin S, Vandekerckhove TT, et al. (2004) Polymerase chain reaction denaturing gradient gel electrophoresis analysis of the N 2 -fixing bacterial diversity in soil under Acacia tortilis ssp. raddiana and Balanites aegyptiaca in the dryland part of Senegal. Environ Microbiol 6: 400-415.

Effect of N-fertilization, plant genotype and environmental conditions on nifH gene pools in roots of rice. Z Tan, T Hurek, B Reinhold-Hurek, Environ Microbiol. 5Tan Z, Hurek T, Reinhold-Hurek B (2003) Effect of N-fertilization, plant genotype and environmental conditions on nifH gene pools in roots of rice. Environ Microbiol 5: 1009-1015.

Ammonia-oxidizing bacteria: a model for molecular microbial ecology. G A Kowalchuk, J R Stephen, Annu Rev Microbiol. 55Kowalchuk GA, Stephen JR (2001) Ammonia-oxidizing bacteria: a model for molecular microbial ecology. Annu Rev Microbiol 55: 485-529.

Genomic studies of uncultivated archaea. C Schleper, G Jurgens, M Jonuscheit, Nat Rev Microbiol. 3Schleper C, Jurgens G, Jonuscheit M (2005) Genomic studies of uncultivated archaea. Nat Rev Microbiol 3: 479-488.

Isolation of an autotrophic ammonia-oxidizing marine archaeon. M Konneke, A E Bernhard, J R De La Torre, C B Walker, J B Waterbury, Nature. 437Konneke M, Bernhard AE, de la Torre JR, Walker CB, Waterbury JB, et al. (2005) Isolation of an autotrophic ammonia-oxidizing marine archaeon. Nature 437: 543-546.

A moderately thermophilic ammonia-oxidizing crenarchaeote from a hot spring. R Hatzenpichler, E V Lebedeva, E Spieck, K Stoecker, A Richter, Proceedings of the National Academy of Sciences of the United States of America. 105Hatzenpichler R, Lebedeva EV, Spieck E, Stoecker K, Richter A, et al. (2008) A moderately thermophilic ammonia-oxidizing crenarchaeote from a hot spring. Proceedings of the National Academy of Sciences of the United States of America 105: 2134-2139.

Archaea predominate among ammonia-oxidizing prokaryotes in soils. S Leininger, T Urich, M Schloter, L Schwark, J Qi, Nature. 442Leininger S, Urich T, Schloter M, Schwark L, Qi J, et al. (2006) Archaea predominate among ammonia-oxidizing prokaryotes in soils. Nature 442: 806-809.

Quantitative analyses of the abundance and composition of ammonia-oxidizing bacteria and ammonia-oxidizing archaea of a Chinese upland red soil under long-term fertilization practices. J Z He, J P Shen, L M Zhang, Y G Zhu, Y M Zheng, Environ Microbiol. 9He JZ, Shen JP, Zhang LM, Zhu YG, Zheng YM, et al. (2007) Quantitative analyses of the abundance and composition of ammonia-oxidizing bacteria and ammonia-oxidizing archaea of a Chinese upland red soil under long-term fertilization practices. Environ Microbiol 9: 2364-2374.

Abundance and composition of ammonia-oxidizing bacteria and ammonia-oxidizing archaea communities of an alkaline sandy loam. J P Shen, L M Zhang, Y G Zhu, J B Zhang, J Z He, Environ Microbiol. 10Shen JP, Zhang LM, Zhu YG, Zhang JB, He JZ (2008) Abundance and composition of ammonia-oxidizing bacteria and ammonia-oxidizing archaea communities of an alkaline sandy loam. Environ Microbiol 10: 1601-1611.

Spatial distribution of ammonia-oxidizing bacteria and archaea across a 44-hectare farm related to ecosystem functioning. E Wessen, M Soderstrom, M Stenberg, D Bru, M Hellman, ISME J. Wessen E, Soderstrom M, Stenberg M, Bru D, Hellman M, et al. (2011) Spatial distribution of ammonia-oxidizing bacteria and archaea across a 44-hectare farm related to ecosystem functioning. ISME J.

Influence of different cultivars on populations of ammonia-oxidizing bacteria in the root environment of rice. A M Briones, S Okabe, Y Umemiya, N B Ramsing, W Reichardt, Appl Environ Microbiol. 68Briones AM, Okabe S, Umemiya Y, Ramsing NB, Reichardt W, et al. (2002) Influence of different cultivars on populations of ammonia-oxidizing bacteria in the root environment of rice. Appl Environ Microbiol 68: 3067-3075.

Community composition of ammoniaoxidizing bacteria and archaea in rice field soil as affected by nitrogen fertilization. Y Wang, X Ke, L Wu, Y Lu, Syst Appl Microbiol. 32Wang Y, Ke X, Wu L, Lu Y (2009) Community composition of ammonia- oxidizing bacteria and archaea in rice field soil as affected by nitrogen fertilization. Syst Appl Microbiol 32: 27-36.

Effects of aboveground grazing on coupling among nitrifier activity, abundance and community structure. Le Roux, X Poly, F Currey, P Commeaux, C Hai, B , Isme J. 2Le Roux X, Poly F, Currey P, Commeaux C, Hai B, et al. (2008) Effects of aboveground grazing on coupling among nitrifier activity, abundance and community structure. Isme J 2: 221-232.

Effects of management regime and plant species on the enzyme activity and genetic structure of N-fixing, denitrifying and nitrifying bacterial communities in grassland soils. A K Patra, L Abbadie, A Clays-Josserand, V Degrange, S J Grayston, Environ Microbiol. 8Patra AK, Abbadie L, Clays-Josserand A, Degrange V, Grayston SJ, et al. (2006) Effects of management regime and plant species on the enzyme activity and genetic structure of N-fixing, denitrifying and nitrifying bacterial communities in grassland soils. Environ Microbiol 8: 1005-1016.

Ecology of denitrifying prokaryotes in agricultural soil. L Philippot, S Hallin, M Schloter, Advances in Agronomy. 96Philippot L, Hallin S, Schloter M (2007) Ecology of denitrifying prokaryotes in agricultural soil. Advances in Agronomy, Vol 96 96: 249-305.

On the fate of anthropogenic nitrogen. W H Schlesinger, Proc Natl Acad Sci U S A. 106Schlesinger WH (2009) On the fate of anthropogenic nitrogen. Proc Natl Acad Sci U S A 106: 203-208.

Community composition and functioning of denitrifying bacteria from adjacent meadow and forest soils. J J Rich, R S Heichen, P J Bottomley, K Cromack, Jr Myrold, D D , Appl Environ Microbiol. 69Rich JJ, Heichen RS, Bottomley PJ, Cromack K, Jr., Myrold DD (2003) Community composition and functioning of denitrifying bacteria from adjacent meadow and forest soils. Appl Environ Microbiol 69: 5974-5982.

Structure and activity of the denitrifying community in a maize-cropped field fertilized with composted pig manure or ammonium nitrate. C Dambreville, S Hallet, C Nguyen, T Morvan, J C Germon, FEMS Microbiol Ecol. 56Dambreville C, Hallet S, Nguyen C, Morvan T, Germon JC, et al. (2006) Structure and activity of the denitrifying community in a maize-cropped field fertilized with composted pig manure or ammonium nitrate. FEMS Microbiol Ecol 56: 119-131.

Relationship between Ncycling communities and ecosystem functioning in a 50-year-old fertilization experiment. S Hallin, C M Jones, M Schloter, L Philippot, Isme J. 3Hallin S, Jones CM, Schloter M, Philippot L (2009) Relationship between N- cycling communities and ecosystem functioning in a 50-year-old fertilization experiment. Isme J 3: 597-605.

Mechanisms of plant species impacts on ecosystem nitrogen cycling. Jmh Knops, K L Bradley, D A Wedin, Ecology Letters. 5Knops JMH, Bradley KL, Wedin DA (2002) Mechanisms of plant species impacts on ecosystem nitrogen cycling. Ecology Letters 5: 454-466.

Microbial activities related to C and N cycling and microbial community structure in the rhizospheres of Pinus sylvestris, Picea abies and Betula pendula seedlings in an organic and mineral soil. O Priha, S J Grayston, T Pennanen, A Smolander, FEMS Microbiol Ecol. 30Priha O, Grayston SJ, Pennanen T, Smolander A (1999) Microbial activities related to C and N cycling and microbial community structure in the rhizospheres of Pinus sylvestris, Picea abies and Betula pendula seedlings in an organic and mineral soil. FEMS Microbiol Ecol 30: 187-199.

Denitrification within the genus Azospirillum and other associative bacteria. K Kloos, A Mergel, C Rosch, H Bothe, Australian Journal of Plant Physiology. 28Kloos K, Mergel A, Rosch C, Bothe H (2001) Denitrification within the genus Azospirillum and other associative bacteria. Australian Journal of Plant Physiology 28: 991-998.

Development and application of real-time PCR assays for quantification of erm genes conferring resistance to macrolides-lincosamides-streptogramin B in livestock manure and manure management systems. J Chen, Z T Yu, F C Michel, T Wittum, M Morrison, Applied and Environmental Microbiology. 73Chen J, Yu ZT, Michel FC, Wittum T, Morrison M (2007) Development and application of real-time PCR assays for quantification of erm genes conferring resistance to macrolides-lincosamides-streptogramin B in livestock manure and manure management systems. Applied and Environmental Microbiology 73: 4407-4416.

Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. S A Bustin, Journal of Molecular Endocrinology. 25Bustin SA (2000) Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. Journal of Molecular Endocri- nology 25: 169-193.

Convenient determination of DNA extraction efficiency using an external DNA recovery standard and quantitativecompetitive PCR. K L Mumy, R H Findlay, Journal of Microbiological Methods. 57Mumy KL, Findlay RH (2004) Convenient determination of DNA extraction efficiency using an external DNA recovery standard and quantitative- competitive PCR. Journal of Microbiological Methods 57: 259-268.

Application of real-time PCR to study effects of ammonium on population size of ammonia-oxidizing bacteria in soil. Y Okano, K R Hristova, C M Leutenegger, L E Jackson, R F Denison, Appl Environ Microbiol. 70Okano Y, Hristova KR, Leutenegger CM, Jackson LE, Denison RF, et al. (2004) Application of real-time PCR to study effects of ammonium on population size of ammonia-oxidizing bacteria in soil. Appl Environ Microbiol 70: 1008-1016.

Analysis of actinomycete communities by specific amplification of genes encoding 16S rRNA and gel-electrophoretic separation in denaturing gradients. H Heuer, M Krsek, P Baker, K Smalla, E M Wellington, Appl Environ Microbiol. 63Heuer H, Krsek M, Baker P, Smalla K, Wellington EM (1997) Analysis of actinomycete communities by specific amplification of genes encoding 16S rRNA and gel-electrophoretic separation in denaturing gradients. Appl Environ Microbiol 63: 3233-3241.

16S/23S rRNA sequencing. D J Lane, Nucleic acid techniques in bacterial systematics. InIn E. Stackebrandt, M. GoodfellowLane DJ (1991) 16S/23S rRNA sequencing. InIn E. Stackebrandt, M. Goodfellow, eds. Nucleic acid techniques in bacterial systematics. pp 115-175.

Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. P D Schloss, S L Westcott, T Ryabin, J R Hall, M Hartmann, Appl Environ Microbiol. 75Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, et al. (2009) Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol 75: 7537-7541.

ARB: a software environment for sequence data. W Ludwig, O Strunk, R Westram, L Richter, H Meier, Nucleic Acids Research. 32Ludwig W, Strunk O, Westram R, Richter L, Meier H, et al. (2004) ARB: a software environment for sequence data. Nucleic Acids Research 32: 1363-1371.

Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness. P D Schloss, J Handelsman, Applied and Environmental Microbiology. 71Schloss PD, Handelsman J (2005) Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness. Applied and Environmental Microbiology 71: 1501-1506.

QIIME allows analysis of high-throughput community sequencing data. J G Caporaso, J Kuczynski, J Stombaugh, K Bittinger, F D Bushman, Nat Methods. 7Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, et al. (2010) QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7: 335-336.

Rapidly denoising pyrosequencing amplicon reads by exploiting rank-abundance distributions. J Reeder, R Knight, Nat Methods. 7Reeder J, Knight R (2010) Rapidly denoising pyrosequencing amplicon reads by exploiting rank-abundance distributions. Nat Methods 7: 668-669.

Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. W Li, A Godzik, Bioinformatics. 22Li W, Godzik A (2006) Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22: 1658-1659.

MUSCLE: multiple sequence alignment with high accuracy and high throughput. R C Edgar, Nucleic Acids Research. 32Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research 32: 1792-1797.

Genome-derived criteria for assigning environmental narG and nosZ sequences to operational taxonomic units of nitrate reducers. K Palmer, H L Drake, M A Horn, Appl Environ Microbiol. 75Palmer K, Drake HL, Horn MA (2009) Genome-derived criteria for assigning environmental narG and nosZ sequences to operational taxonomic units of nitrate reducers. Appl Environ Microbiol 75: 5170-5174.

The lithoautotrophic ammonia-oxidizing bacteria. H-P Koops, U Purkhold, A Pommerening-Röser, G Timmermann, M Wagner, Prokaryotes. 5Koops H-P, Purkhold U, Pommerening-Röser A, Timmermann G, Wagner M (2006) The lithoautotrophic ammonia-oxidizing bacteria. Prokaryotes 5: 788-811.

The taxonomic and functional diversity of microbes at a temperate coastal site: a 'multi-omic' study of seasonal and diel temporal variation. J A Gilbert, D Field, P Swift, S Thomas, D Cummings, PLoS One. 515545Gilbert JA, Field D, Swift P, Thomas S, Cummings D, et al. (2010) The taxonomic and functional diversity of microbes at a temperate coastal site: a 'multi-omic' study of seasonal and diel temporal variation. PLoS One 5: e15545.

Influence of maize mucilage on the diversity and activity of the denitrifying community. E Mounier, S Hallet, D Cheneby, E Benizri, Y Gruet, Environ Microbiol. 6Mounier E, Hallet S, Cheneby D, Benizri E, Gruet Y, et al. (2004) Influence of maize mucilage on the diversity and activity of the denitrifying community. Environ Microbiol 6: 301-312.

Analysis of nifH gene diversity in red soil amended with manure in Jiangxi, South China. Q Teng, B Sun, X Fu, S Li, Z Cui, J Microbiol. 47Teng Q, Sun B, Fu X, Li S, Cui Z, et al. (2009) Analysis of nifH gene diversity in red soil amended with manure in Jiangxi, South China. J Microbiol 47: 135-141.

High diversity of diazotrophs in the forefield of a receding alpine glacier. L Duc, M Noll, B E Meier, H Burgmann, J Zeyer, Microb Ecol. 57Duc L, Noll M, Meier BE, Burgmann H, Zeyer J (2009) High diversity of diazotrophs in the forefield of a receding alpine glacier. Microb Ecol 57: 179-190.

Association of novel and highly diverse acid-tolerant denitrifiers with N 2 O fluxes of an acidic fen. K Palmer, H L Drake, M A Horn, Appl Environ Microbiol. 76Palmer K, Drake HL, Horn MA (2010) Association of novel and highly diverse acid-tolerant denitrifiers with N 2 O fluxes of an acidic fen. Appl Environ Microbiol 76: 1125-1134.

Autotrophic ammonia oxidation by soil thaumarchaea. L M Zhang, P R Offre, J Z He, D T Verhamme, G W Nicol, Proc Natl Acad Sci U S A. 107Zhang LM, Offre PR, He JZ, Verhamme DT, Nicol GW, et al. (2010) Autotrophic ammonia oxidation by soil thaumarchaea. Proc Natl Acad Sci U S A 107: 17240-17245.

Responses of ammoniaoxidizing bacterial and archaeal populations to organic nitrogen amendments in low-nutrient groundwater. D W Reed, J M Smith, C A Francis, Y Fujita, Appl Environ Microbiol. 76Reed DW, Smith JM, Francis CA, Fujita Y (2010) Responses of ammonia- oxidizing bacterial and archaeal populations to organic nitrogen amendments in low-nutrient groundwater. Appl Environ Microbiol 76: 2517-2523.

Distribution and diversity of archaeal and bacterial ammonia oxidizers in salt marsh sediments. N S Moin, K A Nelson, A Bush, A E Bernhard, Appl Environ Microbiol. 75Moin NS, Nelson KA, Bush A, Bernhard AE (2009) Distribution and diversity of archaeal and bacterial ammonia oxidizers in salt marsh sediments. Appl Environ Microbiol 75: 7461-7468.

The influence of soil pH on the diversity, abundance and transcriptional activity of ammonia oxidizing archaea and bacteria. G W Nicol, S Leininger, C Schleper, J I Prosser, Environ Microbiol. 10Nicol GW, Leininger S, Schleper C, Prosser JI (2008) The influence of soil pH on the diversity, abundance and transcriptional activity of ammonia oxidizing archaea and bacteria. Environ Microbiol 10: 2966-2978.

Reproducibility and quantitation of amplicon sequencing-based detection. J Zhou, L Wu, Y Deng, X Zhi, Y H Jiang, ISME J. Zhou J, Wu L, Deng Y, Zhi X, Jiang YH, et al. (2011) Reproducibility and quantitation of amplicon sequencing-based detection. ISME J.

Bacterial community of very wet and acidic subalpine forest and fire-induced grassland soils Plant and Soil. Y T Lin, C P Lin, S M Chaw, W B Whitman, D C Coleman, 10.1007/s11104-010-0308-3Lin YT, Lin CP, Chaw SM, Whitman WB, Coleman DC, et al. (2010) Bacterial community of very wet and acidic subalpine forest and fire-induced grassland soils Plant and Soil DOI:10.1007/s11104-010-0308-3.

The effects of stubble retention and nitrogen application on soil microbial community structure and functional gene abundance under irrigated maize. S A Wakelin, M J Colloff, P R Harvey, P Marschner, A L Gregg, FEMS Microbiol Ecol. 59Wakelin SA, Colloff MJ, Harvey PR, Marschner P, Gregg AL, et al. (2007) The effects of stubble retention and nitrogen application on soil microbial community structure and functional gene abundance under irrigated maize. FEMS Microbiol Ecol 59: 661-670.

Endophytic root colonization of gramineous plants by Herbaspirillum frisingense. M Rothballer, B Eckert, M Schmid, A Fekete, M Schloter, FEMS Microbiol Ecol. 66Rothballer M, Eckert B, Schmid M, Fekete A, Schloter M, et al. (2008) Endophytic root colonization of gramineous plants by Herbaspirillum frisingense. FEMS Microbiol Ecol 66: 85-95.

Novel endophytic nitrogenfixing clostridia from the grass Miscanthus sinensis as revealed by terminal restriction fragment length polymorphism analysis. T Miyamoto, M Kawahara, K Minamisawa, Appl Environ Microbiol. 70Miyamoto T, Kawahara M, Minamisawa K (2004) Novel endophytic nitrogen- fixing clostridia from the grass Miscanthus sinensis as revealed by terminal restriction fragment length polymorphism analysis. Appl Environ Microbiol 70: 6580-6586.

Herbaspirillum frisingense sp. nov., a new nitrogen-fixing bacterial species that occurs in C4-fibre plants. G Kirchhof, B Eckert, M Stoffels, J I Baldani, V M Reis, Int J Syst Evol Microbiol. 51Kirchhof G, Eckert B, Stoffels M, Baldani JI, Reis VM, et al. (2001) Herbaspirillum frisingense sp. nov., a new nitrogen-fixing bacterial species that occurs in C4-fibre plants. Int J Syst Evol Microbiol 51: 157-168.

Quantification of Key Genes Steering the microbial nitrogen cycle in the rhizosphere of sorghum cultivars in tropical agroecosystems. B Hai, N H Diallo, S Sall, F Haesler, K Schauss, Applied and Environmental Microbiology. 75Hai B, Diallo NH, Sall S, Haesler F, Schauss K, et al. (2009) Quantification of Key Genes Steering the microbial nitrogen cycle in the rhizosphere of sorghum cultivars in tropical agroecosystems. Applied and Environmental Microbiology 75: 4993-5000.

Analysis of ammonia-oxidizing bacteria of the beta subdivision of the class Proteobacteria in coastal sand dunes by denaturing gradient gel electrophoresis and sequencing of PCR-amplified 16S ribosomal DNA fragments. G A Kowalchuk, J R Stephen, De Boer, W Prosser, J I Embley, T M , Appl Environ Microbiol. 63Kowalchuk GA, Stephen JR, De Boer W, Prosser JI, Embley TM, et al. (1997) Analysis of ammonia-oxidizing bacteria of the beta subdivision of the class Proteobacteria in coastal sand dunes by denaturing gradient gel electrophoresis and sequencing of PCR-amplified 16S ribosomal DNA fragments. Appl Environ Microbiol 63: 1489-1497.

A new obligately chemolithoautotrophic, nitrite-oxidizing bacterium, Nitrospira moscoviensis sp. nov. and its phylogenetic relationship. S Ehrich, D Behrens, E Lebedeva, W Ludwig, E Bock, Arch Microbiol. 164Ehrich S, Behrens D, Lebedeva E, Ludwig W, Bock E (1995) A new obligately chemolithoautotrophic, nitrite-oxidizing bacterium, Nitrospira moscoviensis sp. nov. and its phylogenetic relationship. Arch Microbiol 164: 16-23.

Growth of ammonia-oxidizing archaea in soil microcosms is inhibited by acetylene. P Offre, J I Prosser, G W Nicol, Fems Microbiology Ecology. 70Offre P, Prosser JI, Nicol GW (2009) Growth of ammonia-oxidizing archaea in soil microcosms is inhibited by acetylene. Fems Microbiology Ecology 70: 99-108.

Growth, activity and temperature responses of ammonia-oxidizing archaea and bacteria in soil microcosms. M Tourna, T E Freitag, G W Nicol, J I Prosser, Environmental Microbiology. 10Tourna M, Freitag TE, Nicol GW, Prosser JI (2008) Growth, activity and temperature responses of ammonia-oxidizing archaea and bacteria in soil microcosms. Environmental Microbiology 10: 1357-1364.

Bacteria rather than Archaea dominate microbial ammonia oxidation in an agricultural soil. Z Jia, R Conrad, Environ Microbiol. Jia Z, Conrad R (2009) Bacteria rather than Archaea dominate microbial ammonia oxidation in an agricultural soil. Environ Microbiol.

Nitrification driven by bacteria and not archaea in nitrogen-rich grassland soils. H J Di, K C Cameron, J P Shen, C S Winefield, M O'callaghan, Nature Geoscience. 2Di HJ, Cameron KC, Shen JP, Winefield CS, O'Callaghan M, et al. (2009) Nitrification driven by bacteria and not archaea in nitrogen-rich grassland soils. Nature Geoscience 2: 621-624.

Ammonia oxidation kinetics determine niche separation of nitrifying Archaea and Bacteria. W Martens-Habbena, P M Berube, H Urakawa, J R De La Torre, D A Stahl, Nature. 461Martens-Habbena W, Berube PM, Urakawa H, de la Torre JR, Stahl DA (2009) Ammonia oxidation kinetics determine niche separation of nitrifying Archaea and Bacteria. Nature 461: 976-979.

Evidence for different contributions of archaea and bacteria to the ammonia-oxidizing potential of diverse Oregon soils. A E Taylor, L H Zeglin, S Dooley, D D Myrold, P J Bottomley, Appl Environ Microbiol. 76Taylor AE, Zeglin LH, Dooley S, Myrold DD, Bottomley PJ (2010) Evidence for different contributions of archaea and bacteria to the ammonia-oxidizing potential of diverse Oregon soils. Appl Environ Microbiol 76: 7691-7698.

Reduced nitrate leaching and enhanced denitrifier activity and efficiency in organically fertilized soils. S B Kramer, J P Reganold, J D Glover, B J Bohannan, H A Mooney, Proc Natl Acad Sci U S A. 103Kramer SB, Reganold JP, Glover JD, Bohannan BJ, Mooney HA (2006) Reduced nitrate leaching and enhanced denitrifier activity and efficiency in organically fertilized soils. Proc Natl Acad Sci U S A 103: 4522-4527.

N fertilization effects on denitrification and N cycling in an aggrading forest. M D Wallenstein, W T Peterjohn, W H Schlesinger, Ecol Appl. 16Wallenstein MD, Peterjohn WT, Schlesinger WH (2006) N fertilization effects on denitrification and N cycling in an aggrading forest. Ecol Appl 16: 2168-2176.

Nitrous oxide emissions from an irrigated cornfield. G L Hutchinson, A R Mosier, Science. 205Hutchinson GL, Mosier AR (1979) Nitrous oxide emissions from an irrigated cornfield. Science 205: 1125-1127.

Maintenance of soil functioning following erosion of microbial diversity. S Wertz, V Degrange, J I Prosser, F Poly, C Commeaux, Environ Microbiol. 8Wertz S, Degrange V, Prosser JI, Poly F, Commeaux C, et al. (2006) Maintenance of soil functioning following erosion of microbial diversity. Environ Microbiol 8: 2162-2169.