Original Research ARTICLE
Evaluation of primers targeting the diazotroph functional gene and development of NifMAP – a bioinformatics pipeline for analyzing nifH amplicon data
- 1Department of Microbiology and Ecosystem Science / Division of Microbial Ecology, Universität Wien, Austria
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 the nifH functional gene, which is in stark contrast to the 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, 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 discriminated 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.
Keywords: Nitrogen Fixation, Primer evaluation, nifH gene, NifMAP, Illumina amplicon sequencing
Received: 28 Jan 2018;
Accepted: 27 Mar 2018.
Edited by:Marcus A. Horn, Leibniz University of Hanover, Germany
Reviewed by:Petri Penttinen, Zhejiang Agriculture and Forestry University, China
Oliver Schmidt, University of Bayreuth, Germany
Copyright: © 2018 Angel, Nepel, Panhölzl, Schmidt, Herbold, Eichorst and Woebken. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
* Correspondence: Dr. Dagmar Woebken, Universität Wien, Department of Microbiology and Ecosystem Science / Division of Microbial Ecology, Althanstrasse 14, Vienna, 1090, Austria, email@example.com