Previous Article | Next Article ![]()
Applied and Environmental Microbiology, February 2007, p. 1166-1173, Vol. 73, No. 4
0099-2240/07/$08.00+0 doi:10.1128/AEM.01803-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Department of Biological Sciences, Stanford University, Stanford, California 94305
Received 31 July 2006/ Accepted 3 December 2006
|
|
|---|
|
|
|---|
Molecular techniques have been used to increase our understanding of the ecology of ammonia oxidizers, especially their distribution in response to environmental conditions. Both ribosomal and "functional" genes (i.e., genes that code for proteins involved directly in ammonia oxidation) have been used as molecular markers for studying the ammonia oxidizers. The gene coding for the
subunit of ammonia monooxygenase (amoA) has been used most frequently as a functional gene marker for studying the phylogeny and distribution of betaproteobacterial ammonia oxidizers (33; reviewed by Avrahami and Conrad [3]). Studies using these two markers have revealed shifts in the community composition of betaproteobacterial ammonia oxidizers in response to changes in the environment, including changes in ammonium concentrations (reviewed by Kowalchuk and Stephen [18]), pH (reviewed by Boer and Kowalchuck [10]), temperature (reviewed by Avrahami and Conrad [3]), and soil moisture (11). However, the interactions among these environmental factors are not well understood.
A number of phylogenetic groups of betaproteobacterial ammonia oxidizers have been defined based on 16S rRNA phylogeny (18, 23). These groups can also be identified based on AmoA phylogeny (4, 6), which is generally concordant with that of 16S rRNA (1, 23). Ecological studies have been performed on members of only a few of these groups. The least understood of these groups is the group "AmoA cluster 10." This group has a representative pure culture: Nitrosospira sp. strain AF, isolated from red sandy soil (pH 4.0) sampled in Kasama, Zambia (32). Nitrosospira sp. strain AF is distinct from the other studied species (i.e., strains 40KI, L115, and B6) in its extremely slow growth (a generation time of 138 h compared to those of three other strains ranging from 50 to 63 h [16]) and tolerance of low pH (as low as pH 5 [16]). In addition, Nitrosospira sp. strain AF showed a higher temperature optimum (31 to 33°C) than three Nitrosospira strains that had been isolated from Nordic soils (26 to 29°C). Environmental sequences closely related to Nitrosospira sp. strain AF have so far been detected only in soil from Tennessee (9, 30) and in Californian grassland soil (13). However, most of the past studies of betaproteobacrial ammonia oxidizers in nature used PCR primers (24) subsequently shown to be unable to detect amoA sequences related to Nitrosospira sp. strains 24C and AF if mixed with other sequences (6). Those studies that did detect cluster 10 sequences either used a degenerated forward primer modified by Stephen et al. (30) from the amoA-1F primer of Rotthauwe et al. (24) or used seminested-PCR amplification with A189 and amoA-2R primers (12, 13, 24). The failure to detect this cluster in more soils may be due to primer biases (3).
We studied the effect of changes in temperature, soil moisture, and fertilizer concentration on potential nitrification activity and the abundance and community structure of betaproteobacterial ammonia oxidizers by applying molecular methods and physiological measurements to Californian grassland soil manipulated in the laboratory. The high relative abundance of cluster 10 ammonia oxidizers in this soil provided a rare opportunity to study the ecology of this group. Our observations suggest that the response of cluster 10 ammonia-oxidizing bacteria to environmental change is different from that of the much-better-studied Nitrosospira multiformis-like ammonia-oxidizing bacteria.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Main characteristics of grassland soil from Jasper Ridge field
|
Potential nitrification activity.
Potential nitrification activity was measured after 7 weeks of incubation (26, 31; modified by Avrahami et al. [6]). Briefly, sterile Erlenmeyer flasks containing 45 ml phosphate buffer (1 mM; pH 7.4), 0.04 ml (NH4)2SO4 (0.25 M), and 5 g of soil were incubated at 25°C on a shaker for 46 h. Samples were taken at five time points, centrifuged for 5 min, and filtered through Acrodisc-CR 25-mm syringe filters (0.2 µm; Life Sciences). Samples were stored at 20°C until analysis of nitrite and nitrate in a DIONEX ion chromatograph system (Sunnyvale, CA) equipped with an AS11-HC analytical column. The sum of nitrite and nitrate concentrations increased linearly with time. Rates of potential nitrification activity were determined from the slope of a linear regression of nitrite plus nitrate concentration on time.
DNA extraction.
Samples for community analysis were taken at the beginning of the experiment and after 7 weeks of incubation. DNA was extracted from 0.5-g soil samples using the Fast DNA spin kit for soil (BIO 101, Carlsbad, CA) and the FastPrep instrument (twice for 30 s each at speed of 5.5) in accordance with the manufacturer's instructions. DNA was further purified using the Wizard DNA cleanup kit (Promega, Madison, WI).
PCR amplification of amoA.
The forward primer used for PCR amplification both for terminal restriction fragment length polymorphism (T-RFLP) and denaturing gradient gel electrophoresis (DGGE) analyses was amoA-1F*, a degenerated primer modified by Stephen et al. (30) from the amoA-1F primer of Rotthauwe et al. (24). This primer was 5' end labeled with 5-carboxyfluorescein for T-RFLP analysis, while for DGGE analysis, a GC clamp (21) was added to the 5' end of the same primer. The degenerated amoA-2R primer of Rotthauwe et al. (24) was used for T-RFLP fingerprinting analysis, while for DGGE, the reverse primer amoA-2R-GG (22) was used. In a previous study, this primer (termed amoAR1) showed the best results for soil samples (6). Amplification was performed by using 1 µM of each primer (QIAGEN, Valencia, CA), 0.48 µg/µl bovine serum albumin (Rosch, Penzberg, Germany), 0.625 units (for T-RFLP) or 1.25 units (for DGGE) AmpliTaq DNA polymerase (AmpliTaq, Applied Biosystems, Foster City, CA), and 12.5 µl (for T-RFLP) or 25 µl (for DGGE) of MasterAmp 2x PCR premix F containing 100 mM Tris-HCl (pH 8.3), 100 mM KCl, 7 mM MgCl2, 400 µM of each deoxynucleoside triphosphate, and PCR enhancer betaine (Epicentre Technologies, Madison, WI). DNA and water (Sigma-Aldrich, Deisenhofen, Germany) were added to a final volume of 25 (for T-RFLP) or 50 µl (for DGGE). All PCRs were performed in triplicate both for T-RFLP and for DGGE. Amplifications were always started by placing PCR tubes into a preheated (94°C) DNA Engine thermal cycler (MJ Research, San Francisco, CA). The thermal profile used for amplification included 5 min at 94°C, followed by 35 cycles of 45 s at 94°C, 60 s at 55°C, and 1 min at 72°C, and 10 min at 72°C for the last cycle.
DGGE fingerprinting analysis.
DGGE was performed as described previously (5) (modified from the method of Muyzer et al. [21]) with slight modifications. We used a DGGE system manufactured by C.B.S. Scientific (Solana Beach, CA). Gels were stained with SYBR green I (Invitrogen-Molecular Probes, Carlsbad, CA) and scanned with a UV scanner and a Kodak camera.
T-RFLP fingerprinting analysis.
T-RFLP fingerprinting analysis was conducted as described previously (13, 15) with some modification. Three PCRs from each sample were pooled and purified using the QIAquick PCR purification kit (QIAGEN, Valencia, CA). Purified PCR products (approximately 100 ng) were digested separately with 5 U of the restriction endonuclease Bsh12362I (MBI Fermentas, Hanover, MD), which has the recognition site GCGC. The digestions were carried out in a total volume of 10 µl for 4 h at 37°C, according to the instructions of the manufacturer. Enzyme inactivation was performed by incubation at 65°C for 20 min. The subsequent T-RFLP fingerprinting analysis was performed at the Genomics Technology Support Facility of Michigan State University (East Lansing, MI; http://genomics.msu.edu/), as described previously (14). The relative abundances of individual T-RFs in a given amoA PCR product were calculated based on the peak height of the individual T-RFs in relation to the total peak height of all T-RFs detected in the respective T-RFLP community fingerprint pattern. The peak heights were automatically quantified by the GeneScan software. To verify the assignments of T-RFs to our detected amoA gene types, we also tested individual clones by T-RFLP fingerprinting analysis.
We chose restriction endonuclease Bsh12362I for our T-RFLP analysis based on the analysis of three clone libraries from our soil. One library consisted of 77 clones (accession numbers AY369266 to AY369342) derived from an April 2000 soil sample (described by Horz et al. [13]), one library of 69 clones (accession numbers DQ393198 to DQ393266) derived from an April 2003 soil sample (S. Avrahami and B. Bohannan, unpublished data), and one library of 48 clones (accession numbers DQ396811 to DQ396858) derived from soil samples incubated for 16 to 20 weeks under laboratory conditions (Avrahami and Bohannan, unpublished). All three libraries contained clones from identical phylogenetic groups, although in different relative abundances. The phylogenetic groups which could be recognized via T-RFLP were (i) AmoA cluster 10 (phylogenetic tree in references 4 and 6), with Nitrosospira sp. strain AF as a representative pure culture (indicated by four T-RFs: 60, 133, 329, and 341 bp); (ii) AmoA cluster 3a (denoted previously as clade JR1 by Horz et al. [13]), with Nitrosospira multiformis as a representative pure culture (indicated by one T-RF of 434 bp); and (iii) AmoA cluster 9 (phylogenetic tree in references 4 and 6), with Nitrosospira sp. strain Nsp65 as a representative in pure culture (indicated by T-RFs 427 and 428 bp in length).
Real-time PCR amplification.
Real-time PCR amplification was modified from the method of Horz et al. (13). In contrast to Horz et al. (13), we used the same primer system (amoA-1F* and amoA-2R) for both rounds of PCR amplification. The first round of PCR amplification was performed as described above using 0.3 µl of DNA per 25-µl reaction mixture and the following thermal cycling protocol: 94°C for 5 min; 12 cycles of 94°C for 45s, 55°C for 60s, and 72°C for 180 s; and 72°C for 10 min for the last cycle. Each sample was amplified in triplicate; the resulting samples later were pooled for the second round of PCR amplification. The second round was performed in triplicate under the same conditions for 35 cycles. Each reaction contained 10 µl DyNAmo Hot Start SYBR green quantitative PCR (qPCR) kit (MJ Research, San Francisco, CA), 1 µM of each primer, and 0.3-µl aliquots of the first round in a final volume of 20 µl. The reaction mixtures were placed in hard-shell, thin-wall 96-well microplates (MJ Research, San Francisco, CA) using Ultra Clear caps (MJ Research, San Francisco, CA) and placed in the DNA Engine Opticon 2 system (MJ Research, San Francisco, CA). The primer annealing temperature was determined experimentally using the temperature gradient feature of the DNA Engine Opticon 2 system. In order to avoid any possible primer dimer interference, the fluorescence reading during each cycle was adjusted to 86°C, a temperature above the melting point of the primer dimers. In addition, to confirm nonspecific amplification, we performed melting curve analyses after each PCR amplification, ensuring the detection of a clear single peak at 90°C corresponding to the 491-bp amoA fragment, together with verification of the size of this fragment by agarose gel electrophoresis.
In order to estimate possible PCR interference and/or inhibition from soil substances, we amplified a plasmid with an insert from Nitrosospira multiformis (accession number AY177933) using the same amoA primer system described above. The plasmid was purified using a QIAprep spin miniprep kit (QIAGEN, Valencia, CA) by following the manufacturer's instructions. The dilution of the purified plasmid DNA template was mixed with 0.3 µl of DNA from the experimental sample and subjected to the real-time PCR protocol described above. The quantification analysis was performed as described previously (13).
Statistical analysis.
All statistical analyses were performed with SAS version 9.1; analysis of variance (ANOVA) was performed using both the GLM and MIXED procedures in SAS (SAS Institute, Inc., Cary, NC). Means were estimated as least square means, and the degrees of freedom were estimated using the Satterthwaite approximation. The data were transformed to achieve normality. For each data set, the transformation that resulted in the best fit of the data set to a normal distribution was chosen. For most data sets, a log transformation was best, with the exception of the nitrate concentration data (square root transformation), the relative abundance of cluster 10 data (arcsine followed by power of 2 transformation), the relative abundance of cluster 3a data (arcsine followed by power of 9.9 transformation), and the relative abundance of cluster 9 (arcsine followed by an inverse transformation).
|
|
|---|
The relative abundance of AmoA cluster 10 was significantly altered by both the temperature and fertilizer treatments (Table 2). There were also significant interactions between these two factors and between temperature and soil moisture (Table 2), such that the effects of either fertilizer or soil moisture were observed at only 30°C (Fig. 1a and b, respectively). At an elevated temperature, either increased soil moisture or increased fertilizer resulted in a decreased relative abundance of cluster 10. When the data from both the 20°C and the 30°C treatments were combined, the relative abundance of cluster 10 did not correlate with any other set of data from the experiment. However, when data from the 30°C treatment were analyzed separately, ammonium concentration was the strongest correlate of the relative abundance of cluster 10 (r = 0.67; P = 0.0023; n = 18), presumably due to the interaction between the fertilizer and temperature treatments. Cluster 3a exhibited the opposite response to temperature, fertilizer, and soil moisture. It increased in its relative abundance at a combination of 30°C and the HF treatment and of 30°C and high soil moisture. Fertilizer treatment was strongly correlated to ammonium concentration (r = 0.92; P < 0.0001). Therefore, the ammonium concentrations were assumed to be due to the fertilizer treatment we imposed. Regression values indicated that ammonium concentration could explain 21% of the variation in the relative abundance of cluster 3a and 58% of the variation when only data from samples incubated at 30°C were analyzed (P = 0.0026 and P = 0.0001, respectively). Cluster 9 also tended to increase in its relative abundance in response to interactions between temperature, soil moisture, and fertilizer (F4,18 = 3.86; P = 0.0195). However, this group could not be detected in the majority of the samples and its highest relative abundance was only 5.6%.
|
View this table: [in a new window] |
TABLE 2. ANOVA analysis for the relative abundance of AmoA cluster 10 dataa
|
![]() View larger version (18K): [in a new window] |
FIG. 1. Effect of (a) temperature and fertilizer concentration and (b) temperature and soil moisture on the relative abundance of AmoA cluster 10. The mean relative abundance of cluster 10 is depicted for all samples, grouped by (a) temperature and soil fertilizer treatments (including all moisture treatments; n = 36) and by (b) temperature and soil moisture treatments (including all fertilizer treatments; n = 36). (a) LF, MF, and HF treatments represent fertilizer concentrations of 0.05, 0.1, and 0.3%, respectively. (b) Soil moisture levels are 30, 45, and 60% of maximum water-holding capacity. Error bars are 95% confidence limits.
|
Incubation increased abundance substantially. The value before incubation was on average only 8.3% (±2.1%) of the value of the sample with the lowest abundance after incubation. Both temperature and soil moisture had significant effects on total abundance (Table 3; Fig. 2). Decreasing temperature as well as increasing soil moisture resulted in an increase in abundance (mean ratio values of 128 and 11 at 20 and 30°C, respectively, and of 10, 47, and 109 for 30, 45, and 60% of WHC, respectively). The highest abundances were found in samples incubated at 20°C, 60% of WHC, and with either MF or HF treatment (Fig. 2).
|
View this table: [in a new window] |
TABLE 3. ANOVA analysis for total abundance of amoA gene dataa
|
![]() View larger version (19K): [in a new window] |
FIG. 2. Effect of temperature, soil moisture, and fertilizer concentration on total abundance of betaproteobacterial ammonia oxidizers. The total abundance was normalized to a standard sample (see the text). The stacked bars represent the contribution of each phylogenetic group to the total abundance. Percent WHC represents the percent of the soil's maximum water-holding capacity. LF, MF, and HF treatments represent fertilizer concentrations of 0.05, 0.1, and 0.3%, respectively. Error bars are 95% confidence limits of the total abundance.
|
|
View this table: [in a new window] |
TABLE 4. ANOVA analysis for "scaled abundance" of AmoA cluster 10 dataa
|
|
View this table: [in a new window] |
TABLE 5. ANOVA analysis for potential nitrification activity dataa
|
![]() View larger version (23K): [in a new window] |
FIG. 3. Effect of (a) temperature and soil moisture, (b) temperature and fertilizer concentration, and (c) soil moisture and fertilizer concentration on potential nitrification activity. The mean potential nitrification activity is depicted for all samples grouped by (a) temperature and soil moisture treatments (including all fertilizer treatments; n = 36), (b) temperature and soil fertilizer treatments (including all moisture treatments; n = 36), and (c) soil moisture and soil fertilizer treatments (including all temperature treatments; n = 36). Soil moisture levels are 30, 45, and 60% of the soil's maximum water-holding capacity. LF, MF and HF represent fertilizer concentrations of 0.05, 0.1, and 0.3%, respectively). Rates were calculated per gram of dry soil. Error bars are 95% confidence limits.
|
![]() View larger version (21K): [in a new window] |
FIG. 4. Effect of (a) temperature and soil moisture and (b) soil moisture and fertilizer concentration on nitrate concentrations. The mean nitrate concentration is depicted for all samples grouped by (a) temperature and soil moisture treatments (including all fertilizer treatments; n = 36) and (b) soil moisture and soil fertilizer treatments (including all temperature treatments; n = 36). Temperatures of 20 and 30°C are displayed as black and gray bars, respectively; soil moisture levels are 30, 45, and 60% of the soil's maximum water-holding capacity. Nitrate concentrations were calculated per gram of dry soil. Error bars are 95% confidence limits.
|
|
|
|---|
Nitrosospira AmoA cluster 10 was the most common group in all samples both before and after incubation (representing 70.6 to 100% of the amoA genes detected). Consistent with previous studies of betaproteobacterial ammonia oxidizers, community shifts could not be detected by DGGE after comparable periods of incubation (e.g., 4 to 6.5 weeks [5, 6]). However, we did detect shifts in this study by T-RFLP, which is known to be more sensitive than DGGE (20). The relative abundance of AmoA cluster 10 was reduced through the synergistic effect of increased temperature, fertilizer, and soil moisture, with a parallel increase in the relative abundance of two other phylogenetic groups (AmoA clusters 3a and 9). It is not clear based on our data alone whether the decrease in relative abundance of AmoA cluster 10 is due to direct competition with these two other groups, a decrease in the fitness of AmoA cluster 10 due to suboptimal growth conditions, or both. However, as discussed above, Nitrosospira sp. strain AF has an optimal growth temperature of 31 to 33°C, which was determined at an ammonium concentration over two times higher than our MF treatment and only slightly lower than our HF treatment (16). If this isolate is characteristic of other cluster 10 organisms, these results indicate that high temperature and high ammonium concentrations are not suboptimal for these microorganisms, suggesting that competitive exclusion may underlie the decrease we observed. Further work with pure cultures would be required to test this hypothesis.
The predominance of AmoA cluster 10 in the other treatment combinations essentially provided an "enrichment culture" of cluster 10 organisms and a unique opportunity to study the response of cluster 10 organisms to environmental change.
The total abundance of members of cluster 10 in these treatments was significantly influenced by temperature and soil moisture. The samples at a combination of 20°C, 60% of WHC, and either MF or HF treatment had the highest total abundance, 3 orders of magnitude higher than that of samples from the other treatments. The soil samples before incubation were taken after the summer period, when the soil was exposed to high temperature (as high as 34°C at 2 cm below surface) and extremely low soil moisture (2.3 g/100 g dws, where dws indicates dry weight of soil, immediately after sampling in the field, equivalent to 9.1% of maximum water-holding capacity). The soil at that point had a very low total abundance of betaproteobacterial ammonia oxidizers, 1 order of magnitude lower than the lowest total abundance found in the samples after incubation. A significantly higher total abundance at 20°C than at 30°C must be due to cell growth and demonstrates the potential for proliferation by cluster 10 organisms under these conditions. This observation was consistent with the observation that increased temperature reduced the scaled abundance of AmoA cluster 10 (i.e., total abundance of AmoA cluster 10) in those treatments where this cluster declined over time.
Temperature had a significant effect not only on abundance but also on activity, although in a different direction. Mean potential nitrification activity was significantly higher at 30°C than at 20°C, although there was an antagonistic interaction between temperature and fertilizer level and a synergistic interaction between temperature and soil moisture. The positive effect of increasing soil moisture on nitrification activity is consistent with other studies (e.g., (29). Nitrate concentrations in the soil were highly correlated with potential nitrification activity, suggesting that in situ nitrification activity may respond to environmental change in a manner similar to potential nitrification activity. These similar trends could be explained by the linear increase in nitrate over the course of the experiment, indicating that the system did not reach a steady state after 7 weeks of incubation, presumably as a result of the fertilizer addition preventing the limitation of substrate.
We were surprised to observe that potential nitrification activity increased with temperature, despite decreases in the abundance of betaproteobacterial ammonia oxidizers. What might explain this observation? The diffusion rate of ammonia into the cell might potentially increase with increased temperature. However, high oxidation rates at such high temperatures have not been reported for other ammonia oxidizers. We can hypothesize that the enzymes of cluster 10 organisms are adapted to warm environments and are more active at higher temperatures than at lower temperatures, which is consistent with observations of the representative pure culture. However, based on the low abundance of cluster 10 organisms at 30°C, this high activity does not appear to provide a growth advantage, even after longer incubation (i.e., 16 to 20 weeks [data not shown]). Why might activity and abundance be decoupled? Being active under high temperatures might provide a short-term advantage through increased cell maintenance and survival under high temperatures and/or by decreasing the lag time before growth when conditions improve. This could be a major advantage in highly seasonal environments with hot, dry periods (such as in California's Mediterranean climate). Such a survival strategy has been proposed for Nitrosospira briensis under starvation conditions (8), but it has not previously been suggested for other stressful conditions or for other strains of Nitrosospira.
There is another potential explanation for our observation. It has been recently reported that archaeal ammonia oxidizers may be widespread and abundant in soils (19). Indeed, amoA genes clustering with those of archaeal ammonia oxidizers have recently been discovered in our soil (C. A. Francis et al., unpublished data) and the contribution of such organisms to potential nitrification activity in these soils is currently unknown. If the contribution of such organisms is substantial, especially at high temperatures, this could explain the decoupling of betaproteobacterial ammonia oxidizers' abundance and potential nitrification activity under these conditions. Studies of archaeal ammonia oxidizers in these soils are currently under way.
This work was supported by awards from the Mellon Foundation, the David and Lucile Packard Foundation, and the National Science Foundation (DEB-0221838 and DEB-0108556).
Published ahead of print on 8 December 2006. ![]()
Present address: Center for Ecology and Evolutionary Biology, University of Oregon, Eugene, OR 97403-5289. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»