Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2002, p. 5685-5692, Vol. 68, No. 11
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.11.5685-5692.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Max-Planck-Institute for Terrestrial Microbiology, Karl-von-Frisch-Strasse, 35043 Marburg, Germany
Received 30 May 2002/ Accepted 6 August 2002
|
|
|---|
|
|
|---|
The aerobic ammonia-oxidizing nitrifiers in soils have mainly been studied by targeting the 16S rRNA genes. The known ammonia oxidizers that exist in pure culture comprise two monophyletic groups (31). One group belongs to the
subdivision of the class Proteobacteria, with Nitrosococus oceani and Nitrosococus halophilus as the only known species. The other group belongs to the ß subdivision of the class Proteobacteria and includes two genera: Nitrosomonas and Nitrosospira. Ammonia-oxidizing bacteria and clones belonging to the ß-Proteobacteria were divided into at least seven clusters based on 16S rRNA gene sequences (31, 39). Soils were dominated by Nitrosospira species of clusters 2, 3, and 4 (8, 17, 29, 38). The community patterns of ammonia oxidizer 16S rRNA clones were found to depend on the soil acidity, with Nitrosospira cluster 2 being established at low soil pH (21, 38). Another selective factor was the concentration of ammonium, with Nitrosospira cluster 3 dominating in early successional soils with relatively high ammonium concentrations while Nitrosospira clusters 2 and 4 dominated in old successional soils with low ammonium concentrations (8, 20, 21). Fields to which manure or compost had been added contained in addition a significant population of Nitrosomonas species (14, 18). Management of grassland soil also affected the ammonia oxidizer populations, which were more diverse in unimproved than in improved soil (40). Soil improved by the addition of N fertilizer was dominated by Nitrosospira cluster 3 and Nitrosomonas cluster 7 (40).
The key enzyme for all aerobic ammonia oxidizers is ammonia monooxygenase. The gene coding for a subunit of this enzyme (amoA) can serve as a useful target for environmental studies, since it reflects the phylogeny of the ammonia oxidizers very well (2, 31, 33). Although the amoA gene is both specific and universal for all ammonia oxidizers, so far it has not been used to the same extent as the 16S rRNA gene in studying soil ammonia oxidizers (19).
For comparison, we also investigated the structure of the denitrifier community and that of the total bacterial community by targeting the nirK genes (5) and the 16S rRNA genes, respectively. nirK encodes the copper-containing nitrite reductase. NirK and the cytochrome cd1-containing nitrite reductase (NirS), which are functionally equivalent although structurally different, are the key enzymes of the denitrification process. Nitrite reductase genes are functional marker genes of denitrifying bacteria, since this physiological group is widespread among phylogenetically unrelated groups. Both genes have been used to detect denitrifiers in environmental samples. In a recent study (30), PCR-based assays detected nirK more readily in soils whereas nirS was preferentially found in marine sediments (6). Community structures of denitrifiers were analyzed by cloning of both genes (6, 30) and by terminal restriction fragment length polymorphism (T-RFLP) analysis of nirS genes (4) from environmental samples. In this study, we introduce community analysis by T-RFLP based on nirK genes.
Our objective was to test how a change in soil ammonium concentration affects the activity and structure of the ammonia oxidizers in comparison to other soil bacteria, denitrifiers in particular.
|
|
|---|
Experimental setup for N2O release.
The experiment described by Gödde and Conrad (13) was slightly modified. Ammonium was added to the soil, since the original ammonium concentration was very low (0.19 µg of NH4+-N g [dry weight] of soil-1). To investigate the effect of the addition of ammonium to the soil, the ammonium content was adjusted to three different concentrations, i.e., 6.5, 58, and 395 µg of NH4+-N g (dry weight) of soil-1, followed by incubation at 4°C. The ammonium concentrations were measured after 5 days of incubation. The release of N2O was measured using a Carlo Erba 8000 gas chromatograph equipped with an electron capture detector as described previously (13). The contribution of nitrification to N2O release was measured after preincubation with and without 10 Pa of acetylene as described previously (13).
Molecular analysis.
Samples were incubated for 4 weeks at 4°C at low, medium, and high ammonium (LA, MA, and HA, respectively) concentrations. The 120-ml serum bottles containing 5 g of soil were opened for aeration every 3 days for 10 min each time to maintain aerobic conditions and stable water content. Approximately 500 mg (wet weight) of soil was transferred from each treatment to a 2-ml screw-cap tube. In addition, samples (ZT) were taken before incubation as a reference.
DNA extraction.
DNA was extracted from these samples following the protocol of Lüdemann et al. (22) with slight modifications. Each sample was mixed with 750 µl of Na-phosphate buffer (120 mM; pH 8) and 187 µl of sodium dodecyl sulfate solution (10% sodium dodecyl sulfate, 0.5 M Tris-HCl [pH 8], 0.1 M NaCl). The samples were resuspended homogenously by vortexing them. After incubation of the solution for 10 min at 60°C, 0.5 g of glass beads (0.17- to 0.18-mm diameter) was added and the suspension was shaken for 1 min at maximum speed in a bead beater (Dismembrator-S; B. Braun Biotech, Melsungen, Germany). After centrifugation (10 min; 20,817 x g at 4°C), the supernatant was collected and extracted three times with phenol (pH 8), phenol-chloroform-isoamyl alcohol (25:24:1; pH 8) and chloroform-isoamyl alcohol (25:24). The addition of ammonium acetate (to a final concentration of 2.5 M) and an equal volume of cold isopropanol followed by 30 min at -20°C and centrifugation (20 min; 20,817 x g at 4°C) allowed the precipitation of DNA. Subsequently, the DNA pellet was washed with 80% ethanol, centrifuged (5 min; 20,817 x g at 4°C), and dried under vacuum. Finally, the DNA was resuspended with 100 µl of TE buffer (10 mM Tris base, 1 mM EDTA, pH 8). DNA was cleaned from humic acids by using polyvinylpolypyrrolidone columns as described by Henckel et al. (15). The purity of the DNA was high as judged by the ratios of absorption at 260 and 280 nm of 1.49 to 1.82 (n = 6).
PCR amplification of amoA.
The primers used for PCR amplification were an amoA set which was described previously by Rotthauwe et al. (33). For denaturant gradient gel electrophoresis (DGGE) analysis, a guanosine-cytosine clamp, as described by Muyzer et al. (26), was added to the 5' end of the forward primer. Amplification was performed by using 0.5 µM each primer, 1 U of AmpliTaq DNA polymerase (Perkin-Elmer Applied Biosystems, Weiterstadt, Germany), and 25 µl of MasterAmp 2x PCR premix F containing 100 mM Tris-HCl (pH 8.3), 100 mM KCl, 7 mM MgCl2, 400 µM each deoxynucleoside triphosphate, and the PCR enhancer betaine (Epicentre Technologies, Madison, Wis.). The DNA was diluted 10 times, due to humic acid inhibition, and 1 µl was added to a final volume of 50 µl. Amplifications were always started by placing the PCR tubes into the preheated (94°C) thermal block of a Mastercycler Gradient thermocycler (Eppendorf, Hamburg, Germany). The thermal profile which was used for amplification was modified after that of Rotthauwe et al. (33) and included 5 min at 94°C followed by 35 cycles of 45 s at 94°C, 30 s at 60°C, and 1 min at 72°C, and 6 min at 72°C for the last cycle.
DGGE and cloning of amoA.
DGGE was performed as described previously with slight modifications (26). PCR products were separated on polyacrylamide gels, using a gradient of 45% (6% [wt/vol] acrylamide-bisacrylamide [37.5:1; Bio-Rad Laboratories GmbH, Munich, Germany], 18% deionized formamide, 3.1 M urea) to 65% (6% [wt/vol] acrylamide-bisacrylamide [37.5:1; Bio-Rad], 26% deionized formamide, 4.5 M urea). The gels were electrophoresed by using the D GENE system (Bio-Rad) with 0.5x Tris-acetate-EDTA at 60°C and 100 V for 17 h. The gels were stained with 1:50,000 (vol/vol) SYBR Green (Biozym, Hessisch-Oldendorf, Germany) and scanned with a Storm 860 PhosphorImager (Molecular Dynamics, Sunnyvale, Calif.).
Bands were excised from DGGE gels using a Dark Reader transilluminator (Clare Chemical Research, Ross on Wye, United Kingdom). The excised bands were resuspended in 200 µl of PCR water, reamplified, and then separated on DGGE again. Due to the appearance of multiple bands in each sample, the reamplified PCR products from the excised bands were ligated to pGEM T-Easy vectors, which were used to transform Escherichia coli JM109 competent cells (Promega, Madison, Wis.) according to the manufacturer's instructions. Clones containing a correct insert were reamplified using DGGE-specific primers and screened by DGGE by comparing the clones to environmental and experimental samples. Different clone types reamplified with amoA primers were sequenced as described below.
Community analysis of bacteria and denitrifiers.
For amplification of 16S ribosomal DNA (rDNA), DNA from soil samples was extracted according to the modified protocol of Lüdemann et al. (22) as described previously for amoA. DNA (1 µl) from a 1:10-diluted extract was used for PCRs (50 µl) as described elsewhere (4) with the bacterial 16S rDNA primers 8-27F [5'-AGAGTTTGATC(A/C)TGGCTCAG-3'] and 1392-1407R (5'-ACGGGCGGTGTGTACA-3') described by Amann et al. (3) and modified by C. L. Moyer (School of Ocean and Earth Science and Technology, University of Hawaii) (unpublished results). For nirK, DNA was extracted using the FastDNA Soil Spin kit (BIO101, Carlsbad, Calif.). PCRs (50 µl) were performed with 1 µl of a 1:2-diluted DNA extract as described elsewhere (6). For T-RFLP, the bacterial forward 16S rDNA primer and the reverse nirK primer were 5'-end labeled with 6-carboxyfluorescein, and a 2 mM concentration of the labeled nirK primer was used. Three PCRs each were pooled and purified using the QIAquick PCR purification kit (Qiagen, Hilden, Germany). The purified PCR products (100 ng) were cleaved with 3 U of HhaI (Promega, Mannheim, Germany) for 16S rDNA and with HaeIII (New England Biolabs, Frankfurt, Germany) for nirK for 3 h at 37°C. The restriction fragments were separated on an ABI 377 automated sequencer for 15 h with limits of 2,500 V and 40 mA. After electrophoresis, the lengths of the fluorescently labeled terminal restriction fragments (T-RFs) were determined by comparison to the internal standards Genescan 2500-ROX (Applied Biosystems, Warrington, United Kingdom) for terminal 16S rDNA fragments and Genescan 1000-ROX (Applied Biosystems) for terminal nirK fragments by using GeneScan 2.1 software (Applied Biosystems). Peaks with a fluorescence of 100 U over the background fluorescence and larger than 50 bp were analyzed by peak height. The relative abundances of T-RFs were determined by calculating the ratio between the height of each peak and the total peak height of all peaks within one sample. Ratios were converted to percentages, and the results are displayed as histograms.
Cloning of nirK genes from soil.
Amplified nirK genes from sample MA were cloned using the TOPO TA cloning kit (Invitrogen, Paisley, United Kingdom) according to the manufacturer's instructions. White colonies picked at random were screened for inserts of the proper size by PCR amplification of the inserts with nirK primers used for T-RFLP under the conditions described previously. The clones were also screened by T-RFLP, and those showing correspondence to T-RFs from the environmental DNA extract were chosen for sequencing. The inserts were PCR amplified using the T7 promoter primer and the M13 reverse primer as described elsewhere (6), except that the annealing temperature was lowered to 55°C.
Sequencing.
PCR products from clones were purified with the QIAquick PCR purification kit (Qiagen). Both strands were sequenced directly from 70 ng of PCR product with vector-specific primers and the ABI Big Dye terminator kit version 2.0 (Applied Biosystems) according to the manufacturer's instructions. Subsequently, excess primers and dye terminators were removed with Autoseq G-50 columns (Amersham-Pharmacia Biotech, Freiburg, Germany), and cycle-sequencing reactions were analyzed with an ABI 377 DNA sequencer.
Phylogenetic analysis.
The sequences were aligned to sequences from the EMBL database using the ARB Fast aligner feature (http://www.arb-home.de). For phylogenetic analysis, 150 deduced amino acid positions for amoA and 102 for nirK were included, omitting insertions or deletions due to ambiguous positional homology. Phylogenetic analyses were performed using ARB and the PHYLIP software package version 3.6a2.1 (12). Trees were reconstructed using the distance matrix-based methods neighbor joining (ARB and PHYLIP) and FITCH (PHYLIP), parsimony (PROTPARS [PHYLIP]), and maximum likelihood (MOLPHY [Institute Pasteur, Paris, France; http://bioweb.pasteur.fr/sequanal/interfaces/molphy.html]).
Nucleotide sequence accession numbers.
The amoA and nirK gene sequences have been deposited in the EMBL nucleotide sequence database under accession no. AY098867 through AY098890 and AJ487545 through AJ487555, respectively.
|
|
|---|
![]() View larger version (20K): [in a new window] |
FIG. 1. Effect of incubation at LA, MA, and HA concentrations on ammonium transformation. The stacked bars indicate percent contributions of nitrification (open) and denitrification (shaded) to total N2O emission, and the small squares indicate the rates of total N2O emission. Means ± standard errors are shown (n = 3).
|
![]() View larger version (86K): [in a new window] |
FIG. 2. DGGE analysis of amoA fragments from soil samples taken before treatment (ZT) and after incubation with LA, MA, and HA concentrations at 4°C for 4 weeks.
|
![]() View larger version (32K): [in a new window] |
FIG. 3. Phylogenetic Fitch-Margoliash tree (using global rearrangement and randomized input order [three jumbles]) based on partial amoA sequences (150 amino acids). Clones obtained from this experiment are shown in boldface. The sources of sequences were soil samples taken before treatment (ZT) and after incubation with LA, MA, and HA concentrations at 4°C for 4 weeks. The scale bar indicates 10 mutations per 100 sequence positions.
|
![]() View larger version (27K): [in a new window] |
FIG. 4. T-RFs of amplified nirK fragments (514 bp) with major differences in relative abundance according to soil treatment. Soil samples were taken before treatment (ZT) and after incubation with LA, MA, and HA concentrations at 4°C for 4 weeks. The key indicates the lengths of the T-RFs.
|
![]() View larger version (48K): [in a new window] |
FIG. 5. Neighbor-joining tree of partial nirK genes based on 102 amino acids. The consensus tree was reconstructed based on neighbor-joining, parsimony, FITCH, and maximum-likelihood (MOLPHY) analyses. Unresolved nodes are displayed as multifurcations and are indicated by dashed lines. Clones obtained from this experiment are shown in boldface, and the calculated sizes of T-RFs are indicated in parentheses. The scale bar represents 10 mutations per 100 sequence positions.
|
|
|
|---|
On the other hand, our study showed that the soil nitrifier population represented by DGGE patterns of the amoA gene did not change significantly during 4 weeks of incubation at 4°C, even when the soil was amended with HA concentrations. Therefore, the response of the activity was probably not the result of a major change in the ammonia-oxidizing population but of a physiological shift. Mendum et al. (23) arrived at a similar conclusion after studying ammonia-oxidizing populations in fertilized and unfertilized soil over a period of 6 weeks. This conclusion was further supported by analysis of the phylogenetic structure of amoA sequences retrieved from the soil ammonia-oxidizing community, which did not reveal any clustering of genes with respect to ammonium treatment. In contrast, our analysis revealed an amazingly high diversity of amoA sequences from the different ammonium treatments, affiliated with five Nitrosospira clusters and one Nitrosomonas cluster. No clones were affiliated with Nitrosospira sp. strain AP. AHB1 from cluster 2, which was isolated from an acidic soil (21, 38). It should also be noted that all amoA sequences had a similarity of >85% (amino acid sequence level) to those of known ammonia oxidizer species. Hence, all the diversity was displayed within a relatively narrow range of species. Generally, the diversity of ammonia oxidizers in tilled agricultural soils is expected to be much lower than in native soils (8). However, our soil had been treated briefly with composted organic material before sampling, which may have led to an input of organisms to the native ammonia oxidizer populations. Although amoA exists in two or three copies in the genome of ammonia oxidizers (7, 27), it is unlikely that part of the diversity is due to multiple copies of amoA. Direct sequencing of amoA genes from 31 pure cultures indicated that partial sequences obtained with the primer set used in our study were unambiguous and thus were not contributing to diversity (1).
The observed invariability of the ammonia oxidizer community upon incubation at 4°C does not support the hypothesis that the pattern of N2O release, which had been observed in earlier studies upon adaptation to particular temperatures (13), might be the result of changes in the community structure. However, other bacterial populations, such as denitrifiers, contribute to N2O release from soils and presumably are affected by the addition of ammonium to soils. Enhanced nitrification rates yield an increase in soil nitrate, which can then be reduced via denitrification. To determine whether changes in the total bacterial community or in the denitrifying population had an impact on the physiological status of the soil, we also explored changes in the community of bacteria and denitrifiers by T-RFLP based on 16S rDNA and nirK genes, respectively. nirS genes were not detected in our soil. However, Priemé et al. (30) successfully amplified nirS genes from a mixed deciduous marsh in Michigan, using the nirS primer sets we tried to use in our study. Similar to the results for the ammonia-oxidizing population, only minor changes were detected for the bacterial community from the untreated soil and after treatment with different ammonium additions. Ammonia oxidizers and denitrifiers comprise only a subfraction of the total soil community. Consequently, in this approach, changes within these groups are supposedly hidden by an overall stable total bacterial community. However, a shift was detected for the denitrifying population based on nirK. Consistent with the shift toward a higher contribution of nitrification to the overall release of N2O upon the addition of medium and large amounts of ammonium, a shift was observed for T-RFLP. Patterns were similar for these two samples but differed from those for the untreated soil and soil treated with a small amount of ammonium. Although T-RFLP has a slightly higher level of resolution than DGGE (24), community shifts for our incubation conditions generally are probably more likely to occur for denitrifiers than for ammonia oxidizers, since the latter are slow-growing organisms (35). Consequently, an adaptation period of 4 weeks might be too limited.
To confirm the specificities of T-RFs, amplified nirK genes from soil sample MA were cloned and screened by T-RFLP, and representative clones were sequenced. Screening revealed that the clones corresponded to major T-RFs, but not all T-RFs detected in the T-RFLP of the soil sample were found in the clone library. The screening results were confirmed by in silico analysis of restriction sites. Generally, the sizes of T-RFs determined by gel electrophoresis seemed to be below the calculated size by 1 or 2 bp. Thus, clones AgMA6 and AgMA36 would represent the T-RF of 59 bp rather than that of 61 bp, although the relative abundance of the 61-bp T-RF exceeded that of the 59-bp T-RF by severalfold and thus it was expected to be cloned preferentially. However, cloning is known to be a factor that possibly introduces bias. Phylogenetic analysis placed nirK genes into soil cluster I, consisting of nirK genes obtained from marsh wetland soil and forested upland soil in Michigan (30). One clone clustered with clones from soil cluster III, consisting of other marsh wetland clones. The origins of these genes remain obscure, since isolates with this genotype are lacking. Thus, the genes seem to be derived from novel, presumably yet-uncultivated denitrifiers present in temperate world soils. Similar to ammonia oxidizers, no nirK clones closely related to those found in a highly acidic soil were found (32), suggesting that pH is indeed a selective factor. However, in their study, Rösch et al. used a different set of nirK primers (32). Thus, soil cluster II was either not detected in our soil due to bias introduced by our PCR primers or it was indeed absent. Surprisingly, in our soil sample, nirK genes were found which show close relationship to genes from denitrifying strains of Ochrobactrum sp. and Alcaligenes faecalis. Genes were also found within the Hyphomicrobium-Rhizobium radiation, although they were not as closely related. Strains of Ochrobactrum spp., Hyphomicrobium spp., and Rhizobium spp. are readily isolated from soils. This suggests that these cloned nirK genes which cluster with the soil bacteria originate from similar soil organisms. However, the nirK genes also might be derived from organisms with a different phylotype, since the phylogeny of nirK is not consistent with the phylogeny of these organisms generally, e.g., Ochrobactrum spp. and A. faecalis, belonging to the
- and ß-Proteobacteria, respectively.
Our results show that rates of N2O emission from soil are positively correlated to soil ammonium concentrations. The contribution of ammonia oxidation to N2O release increases and thus provides the electron acceptor for denitrification. Consequently, both processes lead to higher overall emission rates. Community shifts were negligible for ammonia oxidizers and the total bacterial community but were detected on the level of the faster-growing denitrifiers. Thus, adaptation occurs on the level of physiological status rather than by community shifts.
This study was supported financially by the German Federal Ministry for Education and Research within the BIOLOG Biodiversity Program (01LC0021).
|
|
|---|
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»