Search

Abstract Fixation of atmospheric N2 by free-living diazotrophs accounts for an important proportion of nitrogen naturally introduced to temperate grasslands. The effect of plants or fertilization on the general microbial community has been extensively studied, yet an understanding of the potential combinatorial effects on the community structure and activity of free-living diazotrophs is lacking. In this study we provide a multilevel assessment of the single and interactive effects of different long-term fertilization treatments, plant species and vicinity to roots on the free-living diazotroph community in relation to the general microbial community in grassland soils. We sequenced the dinitrogenase reductase (nifH) and the 16S rRNA genes of bulk soil and root-associated compartments (rhizosphere soil, rhizoplane and root) of two grass species (Arrhenatherum elatius and Anthoxanthum odoratum) and two herb species (Galium album and Plantago lanceolata) growing in Austrian grassland soils treated with different fertilizers (N, P, NPK) since 1960. Overall, fertilization has the strongest effect on the diazotroph and general microbial community structure, however with vicinity to the root, the plant effect increases. Despite the long-term fertilization, plants strongly influence the diazotroph communities emphasizing the complexity of soil microbial communities’ responses to changing nutrient conditions in temperate grasslands.

Diazotrophic microorganisms introduce biologically available nitrogen (N) to the global N cycle through the activity of the nitrogenase enzyme. The genetically conserved dinitrogenase reductase (nifH) gene is phylogenetically distributed across four clusters (I–IV) and is widely used as a marker gene for N2 fixation, permitting investigators to study the genetic diversity of diazotrophs in nature and target potential participants in N2 fixation. To date there have been limited, standardized pipelines for analyzing the nifH functional gene, which is in stark contrast to the 16S rRNA gene. Here we present a bioinformatics pipeline for processing nifH amplicon datasets – NifMAP (“NifH MiSeq Illumina Amplicon Analysis Pipeline”), which as a novel aspect uses Hidden-Markov Models to filter out homologous genes to nifH. By using this pipeline, we evaluated the broadly inclusive primer pairs (Ueda19F–R6, IGK3–DVV, and F2–R6) that target the nifH gene. To evaluate any systematic biases, the nifH gene was amplified with the aforementioned primer pairs in a diverse collection of environmental samples (soils, rhizosphere and roots samples, biological soil crusts and estuarine samples), in addition to a nifH mock community consisting of six phylogenetically diverse members. We noted that all primer pairs co-amplified nifH homologs to varying degrees; up to 90% of the amplicons were nifH homologs with IGK3–DVV in some samples (rhizosphere and roots from tall oat-grass). In regards to specificity, we observed some degree of bias across the primer pairs. For example, primer pair F2–R6 discriminated against cyanobacteria (amongst others), yet captured many sequences from subclusters IIIE and IIIL-N. These aforementioned subclusters were largely missing by the primer pair IGK3–DVV, which also tended to discriminate against Alphaproteobacteria, but amplified sequences within clusters IIIC (affiliated with Clostridia) and clusters IVB and IVC. Primer pair Ueda19F–R6 exhibited the least bias and successfully captured diazotrophs in cluster I and subclusters IIIE, IIIL, IIIM, and IIIN, but tended to discriminate against Firmicutes and subcluster IIIC. Taken together, our newly established bioinformatics pipeline, NifMAP, along with our systematic evaluations of nifH primer pairs permit more robust, high-throughput investigations of diazotrophs in diverse environments.