Functional diversity and abundance of nitrogen cycle-related genes in paddy soil

The nitrogen cycle and the associated microbes play an important role in natural ecosystems, including terrestrial habitats; they also have a major effect on climate change. The aim of this study was to explore microbial communities in rice paddy soil by detecting and quantifying some key functional genes involved in the nitrogen cycle using molecular techniques such as conventional polymerase chain reaction (PCR), clone library construction, sequencing, phylogenetic analysis, and real-time PCR. The genes analyzed were as follows: nitrogenase reductase gene (nifH), hydrazine synthase gene (hzsA), nitrous oxide reductase gene (nosZ), copper-containing (nirK) and cytochrome cd1-containing (nirS) nitrite reductase genes, nitrite oxidoreductase gene (nxrB), and ammonium monooxygenase gene (amoA). The sequence assessment using the clone library targeting these genes revealed high diversity and dominance of bacterial communities. Furthermore, real-time PCR using SYBR green dye and some primers specific for each gene revealed the high abundance of nxrB (4.1 × 109 ± 0.4 × 109 copies g−1 soil) and low abundance of hzsA (4.0 × 105 ± 1.1 × 105 copies g−1 soil). The findings of our study will be useful to explore microbial communities in terrestrial habitats, such as agricultural paddy fields.


Introduction
Soil is a mixture of organic matter, minerals, gases, water, and microorganisms, which together support life on earth. Soil is considered a major component of the earth's ecosystem as it plays various roles owing to its physiochemical heterogeneity; for example, it acts as a carbon reservoir, water storage system, modifier of earth's atmosphere, and habitat for plants and microorganisms including bacteria and fungi [1]. The high heterogeneity of soil may lead to the heterogeneity of habitats. A small amount of soil can contain up to one billion bacteria [2]. Furthermore, the diversity and composition of soil communities determine the functionality of an ecosystem [3]. Moreover, ecosystems are affected by processes occurring in the soil, especially those driven by microbes, such as the nitrogen cycle. The nitrogen cycle and the associated microbial communities play a crucial role in terrestrial habitats; they also have a major effect on the environment. They have been a subject of interest for researchers and the general public because the nitrogen cycle and the associated microbial communities determine the fate and availability of nitrogen, which is essential to support life [4].
The nitrogen cycle involves the transformation of nitrogen from one form to other. Nitrogen fixation, nitrification or ammonia oxidation, and denitrification are the major processes of the nitrogen cycle. Most of these processes are mediated by microbial communities, which obtain energy or accumulate nitrogen via these processes according to their growth requirements. In nitrogen fixation, nitrogen gas is transformed to ammonium by nitrogen-fixing bacteria. This process involves the nitrogenase enzyme, whose multiple subunits are encoded by nifH, nifD, and nifK (Rubio Lm 2002). Among these three genes, nifH has been the most sequenced and extensively used marker gene for studying phylogeny, diversity, and abundance of nitrogen-fixing bacteria [5,6]. In nitrification, ammonia or ammonium is oxidized to nitrite, and this step is catalyzed by the ammonium monooxygenase enzyme encoded by amoA [7]. Nitrite is then oxidized to nitrate, which is catalyzed by nitrite oxidoreductases encoded by genes such as nxrB [8]. Mostly ammonia-oxidizing bacteria (AOB) and ammonia-oxidizing archaea (AOA) mediate the nitrification process [7]. Furthermore, with the discovery of anammox bacteria, it has been reported that ammonia oxidation is also anaerobically mediated by hydrazine synthase encoded by hzsA [9]. During denitrification, nitrate is reduced ultimately to nitrogen gas via a series of enzymatic reactions by diverse nitrogen oxide reductases encoded by nirS, nirK, and nosZ. Any disturbance in this process due to an imbalance in diverse nitrogen compounds might lead to negative consequences such as global warming and climate change. Therefore, understanding, investigating, and obtaining detailed information about the N cycle-related microorganisms are major requisites for improving the current state of the environment.
Although the nitrogen cycle has been investigated and reported in diverse types of soils such as agricultural field soils [10], acidic forest soils [11], paddy soils [12], intertidal soils [13], and wetland soils [14], information about these processes and probes for detecting specific genes is limited. The most widely used phylogenetic marker for studying microbial communities is the 16S rRNA gene. However, some microorganisms performing a common function are not numerically dominant and are hardly detectable using 16S rRNA gene-based probes. Thus, the use of functional gene markers has been the alternative approach for studying environmental samples. Furthermore, this method provides resolution below species level, because of the high evolutionary rates of the less conserved functional molecules [15,16]. Considering the advantages of functional genes over the 16S rRNA gene, in the present study, we detected and quantified nitrogen cyclerelated genes in a rice paddy field using functional gene markers.

Rice paddy soil sampling
A rice paddy field site, located in Iksan, South Korea (latitude 35° 49′ 42.4" N and longitude 127° 02′ 38.4" E), was selected for soil sampling. Soil samples were collected at 20 cm depth from the surface of the field using a small shovel. Soil slurry with water was prepared at a ratio of 1:5 by shaking for 30 min in a shaker to measure pH and electrical conductivity using portable probes (Thermo Scientific, Pittsburgh, PA, USA).

DNA extraction from soil samples and conventional polymerase chain reaction
Genomic DNA was extracted from 0.25 mg of paddy soil using the MoBio Powersoil DNA kit (Carlsbad, CA, USA), following the manufacturer's protocol. The DNA quantity was measured by fluorometry (Qubit 3.0 Fluorometer, Invitrogen, Pittsburgh, PA, USA). Genomic DNA was diluted with ultra-pure water to 1-5 ng µL −1 for further use. All the primers used in this study are summarized in Additional file 1: Table S1. The PCR mixture of total volume 50 µL contained 1-2 µL of 1-5 ng µL −1 genomic DNA and 1 µL each of 10 µM forward and reverse primers. The conditions for the PCR cycles for each target gene are summarized in Additional file 1: Table S2. The quality and quantity of the purified DNA were determined by gel electrophoresis and fluorometer, respectively.

Construction of clone libraries and confirmation of positive clones
Clone libraries were constructed by using DNA fragments of specific band sizes after amplification by PCR and purification. DNA amplicons were ligated into pTOP TA V2 (vector) using the TOP Cloner TA Kit (Enzynomics, Daejon, South Korea) or pLUG-Prime TA Cloning Kit II (Intron Biotechnology, Seongnam, South Korea) and transformed into Escherichia coli DH5α ™ competent cells (Enzynomics, Daejon, South Korea). Blue-white screening was performed using ampicillin (0.1 mg mL −1 ), IPTG (0.1 mM), and X-gal (40 µg mL −1 ) in lysogeny broth (LB) agar plates. White colonies were selected and incubated overnight in LB broth supplemented with 100 µg mL −1 ampicillin. Plasmid DNA extraction was carried out using the Inclone Mini Plasmid Preparation Kit (Inclone biotech, Yongin, South Korea). The plasmid DNAs isolated from the white colonies were subjected to PCR amplification using the multiple cloning site (MCS) primer sets from the respective vector kits, M13F (-20) and M13R (-40) for pTOP-TA-V2 vector and M13F and M13R for pLUG-prime vector. The PCR product was checked on a 1% agarose gel to confirm positive clones. DNA isolated from the clone library was sequenced at Genotech (Daejeon, South Korea).

Phylogenetic analysis
The forward and reverse nucleotide sequences of cloned functional fragments of each gene were assembled together. The obtained sequences were compared with those in the NCBI GenBank by a BLAST search (https ://blast .ncbi. nlm.nih.gov) [17] in order to identify known sequences, for each gene. Phylogenetic trees of the cloned sequences of each gene were constructed along with similar sequences retrieved from the GenBank using Mega 7 [18]. The evolutionary history was inferred using the neighbor-joining method [19,20]. All the sequences were deposited in the GenBank with accession numbers as follows, nifH:

Quantitative PCR analysis
qPCR was used to quantify nitrogen-fixing bacterial gene (nifH), nitrifying bacterial genes (nxrB and amoA), denitrifying bacterial genes (nirK, nirS, and nosZ), and anaerobic ammonia oxidizing bacteria-specific gene (hzsA). It was performed using the primer sets listed in Additional file 1: Table S1 and Greenstar Accupower Mastermix (Bioneer, Daejeon, South Korea), according to the manufacturer's protocol, on the Biorad CFX Connect Real-Time System (Bio-Rad Laboratories, Inc., Hercules, CA, US). To calculate the gene copy numbers in a known concentration of DNA, Eq. 1 was used.
Each reaction was carried out with a mixture of total volume 20 μL, containing 10 μL of 2 × Greenstar qPCR master mix, 1 µL of DNA template, and 1 μL each of forward and reverse primer specific for each gene. To prepare the standard, the confirmed clones of each gene after sequencing were used as shown in Additional file 1: Table S3. Serial dilution was performed from 10 −1 to 10 −7 . The PCR of the standard and unknown samples of each gene was performed in triplicate. The gene copy number of unknown samples was determined using the standard regression curves of gene copy number.

Physico-chemical properties of soil
The soil was found to be slightly alkaline with a pH of 7.86. Soil electrical conductivity (EC 1:5 ) was 12.7 dS m −1 indicating high salinity, as reported for saturated paste equivalent in all soils [21].

Sequencing and phylogenetic analysis
The cloned fragments after the conventional PCR were sequenced, and the results revealed successful cloning of all the genes. After comparing and analyzing the sequences of the cloned genes with the known sequences of bacterial genes retrieved from the NCBI GenBank, the phylogenetic trees were constructed using the retrieved bacterial genes with high sequence similarities.
Here, we explored the diversity of microbial communities in a rice paddy soil by detecting and quantifying some key functional genes involved in the nitrogen cycle by using molecular techniques. Our study will help better understand the nitrogen cycle and the diversity of associated microbial communities in terrestrial habitats, such as paddy fields. Although we used primers to determine the nitrogen cycle-related genes based on the nitrogen-transforming reactions mediated by bacteria, there are several undiscovered nitrogen-transforming reactions that are thermodynamically feasible and the microorganisms that catalyze these reactions.
Additional file 1: Table S1. List of primers used in this study. Table S2. PCR conditions used for each gene. Table S3. Positive clones from each gene used to prepare standard for real-time PCR. Figure S1. PCR amplification of functional genes: nifH, amoA, nxrB, hzsA, nosZ, nirK, nirS genes, respectively. AK: performed the experiments; wrote the manuscript. J-HL: designed the experiments; wrote the manuscript. Both authors read and approved the final manuscript.