Previous Article | Next Article ![]()
Applied and Environmental Microbiology, July 2007, p. 4250-4258, Vol. 73, No. 13
0099-2240/07/$08.00+0 doi:10.1128/AEM.00081-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Agriculture and Agri-Food Canada, Potato Research Centre, Fredericton, New Brunswick E3B4Z7, Canada,1 Nova Scotia Agricultural College, Department of Environmental Sciences, Truro, Nova Scotia B2N5E3, Canada,2 University of Guelph, Department of Environmental Biology, Guelph, Ontario N1G2W1, Canada3
Received 12 January 2007/ Accepted 12 April 2007
|
|
|---|
|
|
|---|
Denitrifying bacteria are widespread in the soil environment and belong to diverse genera (6). Denitrification in bacteria consists of four steps for the reduction of NO3 to N2, catalyzed by nitrate, nitrite, nitric oxide, and nitrous oxide reductases. Denitrification functional genes have been isolated from cultured bacteria, and gene fragments have been amplified from cultured and environmental samples using a range of primers (5, 6, 12, 30). Studies of denitrifying bacteria tend to focus on the steps from nitrite reductase onward, because nitrate reductase is not always linked to complete denitrification but may be related to dissimilatory nitrate reduction to ammonia (27). Broad-range and specific primer sets for real-time quantification of denitrifiers have been developed based on the nitrite reductase genes nirS and nirK (11, 14) and, most recently, the nitrous oxide reductase gene nosZ (15). Nitric oxide reductases are responsible for the conversion of NO to N2O (35). Nitric oxide reductase has two forms in bacteria, cNOR (cytochrome c electron donor) and qNOR (quinol electron donor) (6). cNOR is most commonly associated with denitrifier populations, with qNOR found in some denitrifiers but also in nondenitrifying microorganisms with a detoxification function against NO (6).
Although several researchers have quantified denitrifier populations in soil samples and other environments (11, 14, 18), studies linking analysis of the denitrifying population abundance and actual denitrification rates have yet to be performed. Real-time PCR provides a simple, rapid method for quantifying target genes in complex samples, such as soil and other environmental matrices. The abundance of Pseudomonas stutzeri nirS was quantified using real-time PCR in a range of environmental sample matrices, including groundwater and agricultural soil (11). Broad-range primers to amplify nirK gene fragments from soil samples have been developed previously (14). Henry et al. reported that carbon addition increased the denitrifier population in soil as quantified using the nirK gene; however, the conditions under which this occurred are not well defined and the increased population was not linked to a measure of denitrification activity. Broad-range primers for analysis of the nosZ gene have also been developed (15) and used to quantify these targets and a selection of other denitrification genes using real-time PCR in a range of soil samples.
Literature reviews have identified the need to link the structure, abundance, and function of denitrifying populations to actual denitrification rates to determine the influence of the microbial population in this fundamental process (28, 29). The aims of this study were (i) to develop and validate primers for real-rime PCR quantification based on cnorB genes from cultured, soil-derived denitrifiers and (ii) to use the primers to quantify changes in denitrifier population densities under different carbon addition treatments, chosen to induce different levels of denitrification activity. Denitrification rates were monitored to determine whether population densities were related to actual denitrification in soil microcosms.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Specificity of PCR amplification of cnorB genes from control denitrifier and field isolatesa
|
|
View this table: [in a new window] |
TABLE 2. Primers/probes and conditions used for real-time PCR amplifications
|
Standard curves for DNA quantitation.
External standard curves were generated for each of the three target sequences, cnorBP, cnorBB, and 16S rRNA genes. cnorB PCR products obtained by amplification with primers cnorB2F and cnorB7R (6) were cloned into pGEM-T according to the manufacturer's instructions (Promega). 16S rRNA genes from a representative denitrifier strain were amplified using primers 27f and 1492r (20) and similarly cloned. Plasmid DNAs were extracted using plasmid Mini or Midi kits (QIAGEN, Inc., Mississauga, Ontario, Canada) and quantified using the fluorescent dye PicoGreen according to the manufacturer's specifications (Invitrogen). Prior to use in standard curves, plasmids were restriction digested with appropriate enzymes (single cut to linearize the vector) and heat treated to inactivate enzymes. The sizes of PCR product inserts used for the generation of standard curves were confirmed by sequencing. Copy numbers of plasmid standards were calculated directly from the concentration and length (base pairs) of the extracted plasmid DNA. Standard curves were generated for each target sequence, showing the relationship between cnorB or 16S rRNA gene copy numbers and threshold cycle values. Standard curves were run on each 96-well plate used for real-time PCR. Where possible, all samples from an experiment were run on a single plate; however, when this was not possible, replicate calibrator samples were run on each plate to adjust values obtained. Copy numbers of target sequences in unknown soil DNA extracts were determined from standard curves. 16S rRNA and cnorB copy numbers were not corrected for DNA extraction efficiency from soil.
Pure culture and soil DNA extraction procedures.
Genomic DNA was extracted from overnight cultures of denitrifiers by using the UltraClean microbial DNA kit according to the manufacturer's instructions (MO BIO Laboratories, Inc., Carlsbad, CA). Soil DNA was extracted using a modified procedure based on those in references 10 and 21. In brief, 0.25-g portions of freeze-dried soils were weighed into 2-ml screw-cap tubes containing 0.1 g each of washed and sterile 0.1-mm-diameter, 0.2- to 0.3-mm-diameter, and 0.7- to 1.2-mm-diameter-glass beads plus one 2.5-mm-diameter glass bead (Sigma-Aldrich Canada Ltd., Oakville, Ontario, Canada). Modified hexadecyltrimethylammonium bromide extraction buffer (0.5 ml; equal volumes of 10% [wt/vol] hexadecyltrimethylammonium bromide
in 0.7 M NaCl and 240 mM potassium phosphate buffer, pH 8.0, plus 2.5 mg.ml1 aurintricarboxylic acid
), aluminum ammonium sulfate (50 µl of 200 mM filter-sterilized solution), and phenol-chloroform-isoamyl alcohol (0.5 ml; 25:24:1) were added to tubes and mixed thoroughly. Tubes were shaken at 20 strokes/s for 10 min (MO BIO 96-well plate shaker with tube adapter set) to lyse cells and then centrifuged (16,000 x g, 10 min, 4°C) before removing the aqueous phase to a new tube. The aqueous phase was reextracted with chloroform-isoamyl alcohol (24:1) and centrifuged as above, and the aqueous phase was removed to a fresh tube. Total nucleic acids were precipitated by the addition of 2 volumes 30% (wt/vol) polyethylene glycol 6000 (Fluka BioChemika)-1.6 M NaCl for 2 h at room temperature and then centrifuged (18,000 x g, 30 min, 4°C). Pelleted nucleic acids were washed in ice-cold 70% (vol/vol) ethanol and air dried prior to resuspension in Tris-EDTA buffer (pH 8.0). Extracted DNA was visualized by agarose (0.8% [wt/vol]) gel electrophoresis. DNA was quantified spectrofluorometrically using the fluorescent dye PicoGreen (Invitrogen).
Detection limits and extraction efficiency.
Denitrifying isolates (Bosea sp. strain PD18 and Pseudomonas sp. strain PD21) were inoculated into sterile and nonsterile soil, and DNA was then extracted to determine the extraction efficiencies and detection limits of pure cultures. Different numbers of cells (approximately 108, 106, 104, and 102 cells/g dry weight soil) were added into the soil samples, and dilution plate counts were used to enumerate the added cell populations. Cell numbers inoculated were compared to copy numbers generated in real-time PCR assays. Gene copy numbers were calculated to be equivalent to cell numbers, as cnorB has only ever been found in single copy in the bacterial genome (35).
Soil microcosm experiments. (i) Experiment 1.
A preliminary experiment was conducted to determine whether real-time PCR primers for denitrifier populations and total bacterial population (16S rRNA genes) were capable of detecting differences in these populations in soil over time and to compare this method with biochemical methods for microbial biomass measurements. The experiment manipulated population densities by the addition of large amounts of readily available carbon to soil (500 mg glucose-C/kg soil per day). Fresh field soil (250 g) was placed in 1-liter glass jars and maintained at a moisture content of 29% (wt/wt) ± 1% for the duration of the experiment by the addition of nutrient or control treatments in distilled water (dH2O). Two nutrient treatments were used, a daily addition of 500 mg/kg glucose-C plus 100 mg/kg nitrate-N and a control with no nutrient addition (dH2O instead of nutrient solutions). Each treatment was manually stirred daily during additions of nutrients/water. Stock solutions of glucose and nitrate (in the form of KNO3) were prepared in dH2O for the addition of nutrients to soil microcosms. Sufficient jars were used to allow four replicates of each treatment to be sacrificed at 0, 1, 2, 3, 5, and 7 days. Soil from each jar was analyzed for denitrifier and total bacterial populations by real-time PCR methods as described above and by microbial biomass carbon as described below.
(ii) Experiment 2.
In contrast to the above experiment, where the primary objective was to evaluate the molecular measurement techniques over time in soil microcosms, this experiment was designed to evaluate the population dynamics of total bacteria and components of the denitrifier population under denitrifying conditions and to compare to actual denitrification rates when differing rates of glucose-C were applied. Soil was sampled as above and stored at 4°C prior to use. Soil was repacked into cores as previously described (9), and soil moisture was adjusted to 70% (wt/wt) water-filled pore space to ensure denitrifying conditions (3). Treatments were additions of 0, 125, 250, and 500 mg/kg of glucose-C. Glucose treatment levels were chosen to reflect amounts of available C equivalent to field application levels (9). Nitrate (KNO3-N) was added to all microcosms at 500 mg/kg to ensure nonlimiting concentrations for the duration of the experiment. Sufficient jars were used to allow four replicates of each treatment to be sacrificed at 0 and 2 days. The soil cores were placed in 1-liter glass jars and sealed, and the atmosphere was supplemented with acetylene to a final concentration of 10% (vol/vol). Acetylene blocks the reduction of N2O to N2. Soil treatments were prepared at 4°C and remained at this temperature overnight to allow for the diffusion of acetylene into the soil cores. Jars were then incubated at 25°C in the dark for up to 2 days. Time zero represents the time at which the treatments were transferred to incubation at 25°C. Samples were analyzed for DNA extraction, real-time PCR, and biochemical measurements as described for experiment 1.
Biochemical and analytical methods.
Microbial biomass carbon was measured using the CHCl3 fumigation-extraction method (34). Fumigated and nonfumigated soil samples (25 g) were extracted using 50 ml 0.5 M K2SO4 and analyzed for extractable organic carbon. Microbial biomass carbon was calculated using a factor of 0.35 (34). Nitrate was determined from K2SO4 extracts of nonfumigated samples. Segmented flow analysis (Technicon Industrial Systems, Tarrytown, MA) was used for the colorimetric determination of extractable organic carbon (Technicon Method 455-76 W/A) and NO3 concentrations (Technicon Method 100-70W). Gas analysis was performed for N2O and CO2 production by using a Varian Star 3800 gas chromatograph (Varian, Walnut Creek, CA) fitted with an electron capture detector and a Combi PAL autosampler (CTC Analytics, Zwingen, Switzerland). The electron capture detector was operated at 300°C with 90% Ar, 10% CH4 carrier gas at 20 ml/min in a HayeSep N 80/100 precolumn (0.32-cm diameter by 50-cm length) and HayeSep D 80/100 mesh analytical columns (0.32-cm diameter by 200-cm length) in a column oven operated at 70°C. The precolumn was used in combination with a four-port valve to remove water from samples.
Data analysis.
Analysis of variance was performed using the general linear model of SAS (version 8; SAS Institute Inc., Cary, NC). All nonnormal data were log transformed. Based on a factorial design, means comparisons were performed using the LSMeans test, although if the interaction of the factors showed no significant difference, treatment means were compared using single-degree-of-freedom contrasts. Treatment means and standard errors presented in the tables and figures are calculated from untransformed data. Significance was accepted at a level of probability (P) of <0.05.
|
|
|---|
Detection limits, extraction efficiency, and standard curves.
Nonsterile soil contained background numbers of cnorBP and cnorBB denitrifiers at 1.7 x 105 and 8.3 x 104 copy numbers/g soil, respectively (Table 3). Freshly grown Pseudomonas or Bosea pure cultures inoculated into this nonsterile soil at 106 and 108 cells/g soil could be detected above background levels (Table 3). Inoculated Pseudomonas cells were extracted and detectable at the levels at which they were inoculated (e.g., 108 gene copies were detected when
108 cells were added, or at close to 100% extraction efficiency) (Table 3). Bosea cells were not detected at the levels at which they were added into soil (
10 to 15% extraction efficiency of inoculated cells), indicating difficulties with cell lysis/DNA extraction from these cells in soil (Table 3). Modifications to the DNA extraction method (including size and amount of glass beads, time of bead beating, volume of soil, and buffers) to improve cell lysis were evaluated, but no improvements could be made to increase the extraction efficiency of these cells. Copy numbers are therefore likely to be underestimated when the cnorBB guild is analyzed in native soils. Autoclaved soil was also tested with inoculated cells, but background levels of amplifiable signal were obtained (data not shown). What the source of amplifiable signal was in those soil DNA extracts is not known. Soils were autoclaved three times over 3 consecutive days with incubation in between to remove viable bacteria, but it is possible that DNA from the background population of denitrifiers was still present in autoclaved soil and was subsequently extracted and amplified. Previous work (26) showed a background copy number of
103 for the detection of nahAc in autoclaved soil.
|
View this table: [in a new window] |
TABLE 3. Detection of inoculated denitrifier cells in nonsterile soila
|
![]() View larger version (13K): [in a new window] |
FIG. 1. Standard curves generated from restriction-digested plasmid standards for each of three primer sets, cnorBP, cnorBB, and 16S rRNA genes. Values represent means (n = 3) ± standard errors (error bars are too small to be seen).
|
![]() View larger version (15K): [in a new window] |
FIG. 2. Copy numbers of (A) total bacterial population (16S rRNA genes) and (B) P. mandelii (cnorBP) (squares) and Bosea-Bradyrhizobium-Ensifer (cnorBB) (triangles) denitrifier populations as measured by real-time PCR after the addition of excess glucose-C to aerobic soil microcosms. Solid symbols indicate the addition of glucose, and open symbols indicate controls (no glucose addition). Values are means (n = 4) ± standard errors (error bars).
|
The growth response of the Bosea-Bradyrhizobium-Ensifer (cnorBB) guild to the addition of glucose-C was different from that of the cnorBP population (Fig. 2B). cnorBB copy numbers in the glucose-treated microcosm increased from 2.1 x 105 copy numbers/g soil at day 0 to a maximum of 1.1 x 106 copy numbers/g soil at 7 days, whereas cnorBB copy numbers in the control microcosms increased from 1.39 x 105 at day 0 to a maximum of 4.1 x 105 copy numbers/g soil at day 7. In contrast to the cnorBP population, the cnorBB guild never increased its relative proportion (expressed as a percentage of total 16S rRNA gene copy numbers) above 0.1% in either the glucose-treated or the control microcosms (data not shown).
There was a positive correlation (r = +0.83) between the molecular method (16S rRNA gene real-time PCR) for estimating total bacterial population and a biochemical measure (microbial biomass carbon) of microbial biomass in the soil microcosms (Fig. 3). This indicates that the changes in total bacterial population measured using the molecular method are reflected in the biochemical measurement of microbial biomass carbon. It is recognized that the two measures are not equivalent, however, because microbial biomass carbon measures the total microbial biomass in soil, including bacterial, archaeal, and fungal components, whereas the 16S rRNA gene PCR targets only bacterial ribosomal genes. Previous studies examined the relationship between several measures of the microbial population in soil (such as levels of phospholipid fatty acid, substrate-induced respiration, and total DNA extracted [2, 22]), but to the best of our knowledge, this is the first example of this relationship between quantitative measures of 16S rRNA genes and microbial biomass carbon measurements.
![]() View larger version (13K): [in a new window] |
FIG. 3. Relationship between microbial biomass carbon and 16S rRNA gene copy number in soil microcosms. Values are taken from all microcosms at all time points sampled. The line of best fit indicates the linear relationship described by the following equation: y = 1.95 x 107x + 105 (r2 = 0.69).
|
|
View this table: [in a new window] |
TABLE 4. Denitrifier and total bacterial copy numbers and cumulative N2O-N and CO2-C emissions in response to glucose-C additions to soil microcosms
|
![]() View larger version (12K): [in a new window] |
FIG. 4. Relationship between rate of glucose-C added (mg/kg) to soil microcosms and the size of the cnorBP population as a percentage of total population measured by 16S gene copy number at time 2 days. The line of best fit indicates the linear relationship described by the following equation: y = 0.00116x + 0.0570. Values are means (n = 4) ± standard errors (error bars).
|
|
|
|---|
The nitric oxide reductase gene (cnorB) was used as a marker for the denitrification populations targeted in this study. Denitrification gene sequences were previously amplified and sequenced from a culture collection of denitrifiers isolated from potato field soil (7). cnorB was present in a majority of the denitrifiers isolated, and we were able to target specific bacterial denitrifier groups (Pseudomonas mandelii-like strains and the Bosea-Bradyrhizobium-Ensifer guild) using this gene. The primer sets developed were determined to be specific and sensitive enough for the detection of changes in denitrifying bacterial populations in soil microcosms. The sensitivity obtained was similar to that in a previous study (15), where a detection limit of 10 to 100 target molecules per assay was achieved, equivalent to 103 to 104/g dry soil for NTC samples. The use of a third step (80°C) in the PCR cycle was implemented with our cnorB primer sets. The fluorescence data acquired avoided any signal from primer dimers or lower-melting-temperature nonspecific products (15); however, in most cases this was simply a precaution, as few nonspecific products were observed.
The populations analyzed in these experiments were small components of the total bacterial population and small components of the total (as opposed to culturable) denitrifier population, but they are found in a range of soil samples obtained from our potato field sites (7). The response of these populations to denitrifying conditions provides knowledge of denitrifier population dynamics in this process in a natural soil. Depending on the availability and composition of nutrients and habitat in soil, the soil bacterial community may comprise different proportions of r- and K-strategists (32). The P. mandelii (cnorBP) population showed a classical r-strategy response, with rapid growth in response to the addition of glucose, a readily available carbon source. This population increased its relative proportion in the total population, indicating the ability of this population to successfully compete for available resources. In contrast, the Bosea-Bradyrhizobium-Ensifer (cnorBB) guild did not increase in response to glucose addition, indicating a different growth strategy (K-strategist) under these conditions. Glucose is often used as a substrate in soil studies since the majority of soil microorganisms can metabolize it (32). The response to other organic substrates, such as crop residues, may be substantially different for the denitrifying populations identified here.
Although real-time PCR provides a simple rapid method for the quantification of bacterial populations in soil, the values obtained may not be accurate or "absolute" for a number of reasons. The total population measured here was estimated using 16S rRNA gene primers and a TaqMan probe. The use of a probe and the limitations of primer design, while required for PCR specificity, may lead to an underestimation of the total population of bacteria if the probe and primer do not bind to all possible amplified target regions (33). The estimated abundances of denitrifier populations may not equal the true abundance of these groups in soil due to differences and biases in extraction efficiency. We showed that freshly grown Bosea cells inoculated into soil were not extracted efficiently. Whether this problem would also apply to native populations of this bacterium is unknown, but it is likely that extraction bias significantly affects the absolute numbers of targets extracted and subsequently amplified in real-time PCR. Although absolute numbers may not be achievable, gross differences and changes in population size are still detectable. The differences observed between the two denitrifier populations studied are then real differences in the responses of these populations to the conditions tested.
Denitrification was measured and observed to increase in response to increasing glucose additions, along with soil respiration. This increase in denitrification rate may result from the growth and/or revival (25) of denitrifier population (and hence increase in abundance of active denitrifiers), leading to an overall increase in the denitrification rate in soils, or from the growth of the total population, leading to increased respiration, lowered oxygen levels, and subsequent induction of denitrification genes in the denitrifier population. We have shown differences in the growth responses of two groups of the denitrifier population to the addition of an available carbon source. Although we measured overall denitrification activity, the specific contribution of each denitrifier population to denitrification was not measured. mRNA analysis (real-time reverse transcription-PCR) is the next step towards understanding the different factors that influence not only population density but also functional gene activity in response to conditions that influence denitrification and nitrous oxide emissions in soils.
The targeting of specific populations provided knowledge of the responses of subpopulations of denitrifiers to changing conditions and provides insights into factors that contribute to denitrifier population density in the field. An analysis of larger components of the denitrifier population, through the use of broad-range primer pairs for denitrification genes as described in references 14, 15, and 18, would provide a better overall understanding of the relationships of denitrifier communities to denitrifying conditions and carbon amendments. Absolute quantification of total denitrifier populations is still not technically feasible using current techniques due to the diversity of gene sequences for functional genes and to the diversity of phylogenetic genes (i.e., 16S rRNA genes) from which denitrifiers are obtained. Recent work (16) has highlighted the limitations of current broad-range primer sets for nitrite reductase genes (nirK and nirS) in cultured and uncultured denitrifiers. Two primer sets developed for nosZ, while both amplifying a broad range of target genes in soil samples, amplified different sets of sequences, verified when PCR products were sequenced (15). Primers have been developed based on a few conserved full-length genes and are unlikely to cover the entire diversity of denitrifiers in any given environment.
This study demonstrates the feasibility of using a molecular approach to understanding the effects of different treatments on the dynamics of denitrifier populations in incubation experiments. Such approaches should be applicable under field conditions. The application of broad-range primer sets as well as the more specific primer sets developed in this study will allow a comparison of the overall denitrifier population dynamics in response to environmental inputs as well as some of the internal dynamics of denitrifier population turnover in a mixed soil population.
Funding for this project was supplied by the GAPS program of Agriculture and Agri-Food Canada and an NSERC (Canada) team strategic grant. M.N.M. was the recipient of an NSERC postgraduate scholarship.
Published ahead of print on 20 April 2007. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»