Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2002, p. 20-30, Vol. 68, No. 1
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.1.20-30.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Department of Molecular and Cell Biology, Institute of Medical Sciences, University of Aberdeen, Foresterhill, Aberdeen, AB25 2ZD,1 Department of Zoology, Natural History Museum, London, SW7 5BD, United Kingdom2
Received 18 April 2001/ Accepted 9 October 2001
|
|
|---|
|
|
|---|
Investigation of microbial diversity in natural communities has been made possible by the development of molecular fingerprinting techniques, in particular those based on analysis of 16S rRNA gene sequences, but the majority of studies have assessed "total" diversity, using primers which amplify 16S rRNA genes from bacterial and archaeal domains. The considerable diversity within these groups makes detection of treatment effects difficult (23) and analysis of narrower groups has been proposed (30).
Nitrification, the conversion of ammonium to nitrate via nitrite, is critical to the cycling of nitrogen in terrestrial environments, increasing nitrogen losses through leaching and denitrification of nitrate (34). This process is carried out in most ecosystems by autotrophic bacteria and is often rate limited by the activities of ammonia-oxidizing bacteria that are responsible for the oxidation of ammonia to nitrite by the enzyme ammonia monooxygenase. Ammonia oxidizers are therefore likely to be influenced by traditional grassland management regimes, such as N fertilization and soil liming, and populations of these organisms increase within managed environments (2, 20, 26).
Despite their significance in a range of environments, the slow growth and low yield of ammonia oxidizers observed during laboratory cultivation have restricted physiological and ecological studies on these organisms (34, 38). However, the application of 16S ribosomal DNA (rDNA)-based techniques circumvents many of these problems and has enabled ecological studies of ammonia oxidizer populations within a range of environments (22, 24, 27, 39, 40, 41). Comparative 16S ribosomal gene sequencing analysis has placed autotrophic ammonia oxidizers into two distinct monophyletic groups (8, 44). The first group, characterized by a limited number of marine nitrosococci (15, 49), belongs to the
-subclass of the Proteobacteria. The second group, within the ß -proteobacteria, includes most of the known isolates and contains the genera Nitrosomonas and Nitrosospira (8). Further molecular analysis, based on 16S rDNA sequences of natural populations of ß-proteobacterial ammonia oxidizers, from a range of different environments, has demonstrated the existence of at least seven sequence clusters (Nitrosospira clusters 1 to 4; Nitrosomonas clusters 5 to 7) (39). The presence and relative abundance of different clusters appear to be related to environmental factors (16, 24, 32, 39, 40), but the physiological basis for the observed differences is unknown. For example, sequences representative of Nitrosospira cluster 2 have been found in greater relative abundance in acid agricultural soils, while cluster 3 sequences appear to be more common in neutral agricultural soil, suggesting the possibility that cluster 2 organisms are, in some way, adapted to growth at low pH. Such differences in community structure may explain the occurrence of nitrification in acid soils (3, 4, 47), while autotrophic nitrification in liquid culture rarely occurs below pH 7 (1).
The availability of techniques for analysis of natural populations of ammonia-oxidizing bacteria, their importance in biogeochemical cycling, and the potential effects of N fertilizer applications on community structure provide an ideal system for assessing the effects of anthropogenic activities on microbial diversity and heterogeneity. This was facilitated by the availability at the chosen study site of soils from three different types of British upland grass pasture: improved, semi-improved, and unimproved grasslands (37). These grasslands differ with respect to different management practices, including sowing of different grass species, N fertilization, and sheep grazing, operated over a period of 30 years. N fertilization, in particular, is likely to influence the diversity and community structure of ammonia-oxidizing bacteria. Ammonia oxidizer populations within the different soils were therefore investigated by denaturing gradient gel electrophoresis (DGGE), probe hybridization, and phylogenetic clustering of 16S rRNA and amoA genes (encoding the
-subunit of ammonia monooxygenase) amplified directly from extracted DNA.
|
|
|---|
At each sampling site, three 5-by-5-m quadrats (plots) were randomly located, and five soil cores (diameter, 10 cm; depth, 5 cm) from each replicate plot were collected, pooled, and sieved (mesh size, <2 mm) to remove stones and plant material. Fresh, sieved soil was used for measurement of pH, nitrification potential, and other soil characteristics. Subsamples (5 to 10 g) of sieved soil were stored at70°C until required for molecular analyses.
Chemical analyses.
pH was measured in 1:2 (wt/vol) suspensions of soil in either 0.01 M CaCl2 or distilled water. The soil water content was determined as weight loss after drying for 24 h at 105°C. Soil organic matter was determined as loss on ignition at 430°C for 24 h. Total N and total C measurements were undertaken on oven-dried soil (48) by using an automated Fison NA 1500 NCS Analyser mass spectrophotometer (Elemental Microanalysis Ltd., Oakhampton, United Kingdom). Exchangeable ammonium and nitrate were determined after 1:4 (wt/vol) extraction of soil with 2 M KCl. Ammonia and nitrate concentrations were measured colorimetrically by using an Alpkem RFA Autoanalyser (Alpkem Corporation, Clackamas, Oreg.). Soil nitrification potential was measured over 7 days by the shaken soil-slurry method (7) with an initial ammonium sulfate concentration of 50 mM and regular sampling for analysis of the nitrite and nitrate concentrations. Statistical analysis by one-way analysis of variance was undertaken on triplicate samples of each replicate plot by using the Fisher exact test.
DNA extraction.
DNA was extracted from 0.5 g of soil, and samples were disrupted with a Hybaid Ribolyser Cell Disrupter (Hybaid, Ltd., Ashford, United Kingdom) (speed 4, 10 s; two times). Extractions were carried out in 2-ml screw-cap Blue Matrix Ribolyser tubes containing a mixture of ceramic and silica beads (Hybaid), 300 µl of 0.12 M sodium phosphate buffer (pH 8.0), 500 µl of Tris-buffered phenol (Fisher Scientific U.K., Loughborough, United Kingdom), and 200 µl of 1 M Tris-HCl (pH 8). After centrifugation at 13,000 rpm for 10 min, the supernatant was transferred to a 1.5-ml Eppendorf tube and reextracted with an equal volume of Tris-buffered phenol, followed by an equal volume of chloroform. The upper aqueous phase was dialyzed, to concentrate the DNA and remove humic acids by using a Microcon YM-100 centrifugal filter device (Millipore Corporation, Bedford, Mass.) with four changes of 100 µl of TE buffer (10 mM Tris, 1 mM EDTA; pH 8). Concentrated DNA was resuspended in 50µl of TE buffer before being further purified by standard low-melting-point agarose (1.0% [wt/vol]) gel electrophoresis with
DNA EcoRI/HindIII markers (NBL Gene Sciences, Cramlington, United Kingdom). High-molecular-weight DNA (>10 kb) was removed and subsequently cleaned by using a Bio-Rad Freeze and Squeeze spin column (Bio-Rad Laboratories, Hemel Hempstead, United Kingdom) according to the manufacturers instructions.
PCR amplification of 16S rRNA gene fragments.
Extracted DNA was amplified by nested PCR by using the ammonia oxidizer specific, 16S rDNA primer sets ßAMOf-ßAMOr (22) and CTO189f-CTO654r (16) and the amoA specific primer sets amoA-2F (AARGCGGCSAAGATGCCGCC)-amoA-5R (TTATTTGATCCCCTC) and amoA-3F (ACCTACCACATGCACTT)-amoA-4R (GGGTAGTGYGACCACCAGTA), supplied by Karl-Paul Witzel, Max-Planck-Institut für Limnologie, Plön, Germany. The amoA primers 2F and 5R bind at positions 279 to 298 and positions 1065 to 1079 of the amoA gene sequence of Nitrosomonas europaea (25) amplifying a 800-bp fragment (positions 31 to 50 and 817 to 831 of the open reading frame). The 3F and 4R primers bind at positions 369 to 385 and 588 to 607 (positions 121 to 137 and 339 to 358 of the open reading frame), respectively, generating a 238-bp fragment. The primers and the PCR conditions for each primer pair are summarized in Table 1. PCR buffer (Bioline, London, United Kingdom), 1 U of Biotaq DNA polymerase (Bioline), 1 µl of each primer (final concentrations of 0.1 [amoA] and 0.4 pmol [16S rDNA] µl1), each deoxynucleoside triphosphate (Bioline) at a concentration of 250µM, and 1 µl of template DNA were added to each 50-µl (total volume) reaction mixture in a 0.5-ml test tube stored on ice. Amplifications were started by placing cooled tubes into a preheated (94°C) Omn-E thermal cycler (Hybaid Ltd.), followed by the thermal cycle conditions shown in Table 1. PCR products (5 µl) were checked by standard agarose (1.2% [wt/vol]) gel electrophoresis with visualization of DNA by ethidium bromide fluorescence and Bioline HyperLadder DNA markers (Bioline) prior to analysis by DGGE. PCR amplification was also undertaken on representative 16S rDNA clones or pure cultures for each currently recognized cluster of ß-subgroup ammonia oxidizers (39), for use as reference sequences for 16S rRNA gene DGGE analysis, as follows: EnvB1-8 (Nitrosospira cluster 1), pH4.2A/27 (Nitrosospira cluster 2), pH4.2A/4 (Nitrosospira cluster 3), pH7B/C3 (Nitrosospira cluster 4), EnvA1-21 (Nitrosomonas cluster 5), EnvC1-19 (Nitrosomonas cluster 6), and N. europaea (Nitrosomonas cluster 7). For amoA gene DGGE analysis, the following pure cultures of ammonia oxidizing bacteria were used as references:Nitrosospira sp. strain B6 (cluster 2), Nitrosospira sp. strain AV (cluster 3),Nitrosospira sp. strain 40KI (cluster 4), N. europaea (cluster 7), and N. eutropha (cluster 7).
|
View this table: [in a new window] |
TABLE 1. Summary of phylogenetic primers, functional primers, and PCR conditions used in this study
|
Recovery and sequence analysis of bands from DGGE gels.
Samples of individual DGGE bands were reamplified with 1 µl as a template from a suspension of the excised band resuspended in 20 µl of sterile distilled water. Reamplified PCR products of excised DGGE bands were purified by dialysis with sterile distilled water by using a Microcon YM-100 centrifugal filter device (Millipore). The purified PCR product was then resuspended in about 20 µl of sterile distilled water to give a final concentration of DNA of between 15 and 20 ngµl1 as determined by standard agarose (1.2% [wt/vol]) gel electrophoresis and Bioline HyperLadder DNA markers (Bioline). Sequencing reactions were performed by using the BigDye Terminator Cycle Sequencing Kit (PE Biosystems, Warrington, United Kingdom) with either CTO189f or amoA-2F primer, and the cycle sequencing products were analyzed with a Model ABI377 automated sequencer (PE Biosystems).
Sequence chromatographs were analyzed by using the Chromas software package (version 1.42; C. McCarthy, Griffith University, Brisbane, Queensland, Australia). Partial nucleotide sequences (272 bp for 16S rDNA and 186 bp for amoA) were aligned manually by using GDE version 2.2 running in ARB (Technical University of Munich, Munich, Germany) with sequences retrieved from the GenBank database (amoA sequences) or from ammonia oxidizer 16S ribosomal sequence databases obtained from previous studies (24, 32, 39).
Phylogenetic relationships between pairs of 16S rDNA sequences were calculated by using both distance and maximum parsimony. For distance matrix analyses, LogDet/Paralinear distances (17) of variable sites (18) was used and implemented in PAUP version 4.0b4a (43). The maximum-likelihood method was used to estimate the proportion of variable sites (10). All LogDet/Paralinear distances trees were constructed by using the minimum evolution criterion, and the data were bootstrapped 100 times.
Ammonia monooxygenase (amoA) nucleotide sequences were analyzed as described above by using LogDet/Paralinear distances of variable sites. In addition, LogDet/Paralinear distances were carried out with codon position 3 excluded from the phylogenetic analysis in order to compensate for potential mutational saturation and nucleotide base composition heterogeneity (11). Trees were constructed as described above, and a 50% majority rule consensus tree was calculated. Base composition data were calculated by using the function in PAUP version 4.0b4a.
Probe labeling, DNA transfer, and hybridization.
Cluster-specific and genus-specific ammonia oxidizer oligonucleotide probes (40) were 3' end labeled by using terminal transferase (50 U) in the supplied buffer with 0.05 mM digoxigenin-11-ddUTP (DIG label) at 37°C according to the manufacturers instructions (Roche Diagnostics Ltd., Lewes, UK). DIG-labeled probes were reoptimized for specificity by using membrane-bound target ammonia oxidizer sequences and checking probe hybridization or cross-reactivity over a range of hybridization and washing temperatures. The optimized hybridization temperatures and number of mismatches to target and nontarget sequences for the DIG-labeled probes are summarized in Table 2.
|
View this table: [in a new window] |
TABLE 2. Oligonucleotide probes, hybridization oven temperatures, and the number of mismatches to both target and nontarget sequences for the DIG-labeled probes for analysis of ammonia oxidizer 16S rDNA fragments
|
Prehybridization and hybridization were carried out in standard hybridization buffer (5x SSC, 0.1% [wt/vol]N-lauroylsarcosine, 0.1% [wt/vol] sodium dodecyl sulfate [SDS], 1% [wt/vol] blocking reagent [Roche Diagnostics, Ltd.]) at the temperatures shown in Table 2 in a Hybaid Shake n Stack hybridization oven (Hybaid, Ltd.). Membranes were washed free of unbound and nonspecifically bound probe by washing them in 2x SSC with 0.1% (wt/vol) SDS at room temperature (2 times; 5 min), followed by two stringency washes for 20 min each at the hybridization temperature with 0.5x SSC and 0.1% (wt/vol) SDS.
Chemiluminescent detection of the hybridised probe was carried out by first equilibrating the membrane for 1 min in maleic acid buffer (100 mM maleic acid; 150 mM NaCl; pH 7.5, 0.3% [vol/vol] Tween 20) before blocking for 30 to 60 min in 1% (wt/vol) blocking reagent in maleic acid buffer. After blocking, the membrane was incubated in a 1:10,000 (vol/vol) dilution of anti-digoxigenin Fab fragments conjugated to alkaline phosphatase (Roche Diagnostics, Ltd.) for 30 min. Membranes were washed free of unbound antibody by washing the membrane twice, 15 min per wash, in washing buffer before incubation with a 1:100 (vol/vol) dilution of CSPD chemiluminescent substrate (disodium-3-(4-methoxyspiro(1,2-dioxetane-3,2'-(5' -chloro)tricyclo[3.3.1.13,7]decan)-4-yl)phenyl phosphate) as recommended by the manufacturer. Detection of the chemiluminescent signal was undertaken by exposure of the membrane to X-ray film (Fuji-RX X-ray film) for between 10 min and 1 h at room temperature. Membranes were stripped prior to reprobing by two washes in 0.2 M NaOH with 0.1% (wt/vol) SDS for 10 min each at 37°C, followed by rinsing in 2x SSC.
Nucleotide sequence accession numbers.
Accession numbers from the GenBank database of previously determined ammonia oxidizer 16S rRNA and amoA gene sequences used in the analyses are as follows. For 16S rDNA the sequences are N. europaea (AF037106), N. eutropha (M96402), N. ureae (Z46993), Nitrosococcus mobilis (M96403), N. marina (M96400), Nitrosospira sp. strain 40KI (X84656),Nitrosospira sp. strain D11 (X84660), Nitrosospira multiformis (M96401), Nitrosospira sp. strain AV (Y10127) Nitrosospira sp. strain C-141 (M96397), N. tenuis (M96404), Nitrosospira sp. strain AHB1 (X90820),Nitrosospira sp. strain B6 (X84657), clone EnvA2-13 (Z69097), clone EnvA2-17 (Z69098), clone EnvB2-6 (Z69106), clone EnvC1-18 (Z69117), clone EnvC2-4 (Z69121), clone EnvC2-6 (Z69122), clone pH4.2A/3E (Z69155), clone pH4.2A/F2 (Z69163), clone pH4.2A/G2 (Z69164), clone pH4.2A/K2 (Z69167), clone pH4.2A/L2 (Z69168), clone pH7C/37 (Z69192), clone pH7C/53 (Z69195), enrichment culture E20_pH4/5/UP (AF178104), and enrichment culture E29_pH7/50/UL (AF178112). For amoA the sequences are N. europaea (L08050), N. eutropha amoA1 (U51630), N. eutropha amoA2 (U72670), Nitrosospira sp. strain AV amoA1 (AF032438), Nitrosospira sp. strain AV amoA2 (AF016003), Nitrosospira sp. strain AV amoA3 (U92432), Nitrosospira tenuis (U76552),Nitrosospira multiformis amoA1 (U91603), N. multiformis amoA2 (U15733), N. multiformis amoA3 (U89833), Nitrosospira sp. strain AHB1 (X90821),Nitrosospira sp. strain 39-19 amoA1 (AF042170),Nitrosospira sp. strain 39-19 amoA2 (AF016002), and Nitrosospira sp. strain 39-19 amoA3 (AF006692).
Sequences of partial amoA and 16S rRNA gene fragments of excised DGGE bands have been deposited under accession numbers AJ308604 to AJ308622.
|
|
|---|
|
View this table: [in a new window] |
TABLE 3. Characteristics of grassland soils at the time of sampling (August 1998)
|
PCR amplification.
DNA was thought to be sufficiently pure for molecular analysis since all DNA samples could be amplified readily with the 16S rDNA eubacterial primer set, 357f/518r (reference 29 and data not shown) However, PCR amplification performed with a single set of primers for either ammonia oxidizer-specific 16S rRNA or amoA genes did not consistently yield detectable PCR products. It was therefore necessary to undertake nested PCR, with primer set ßAMOf-ßAMOr, followed by CTO189f-CTO654r, for 16S rRNA genes and primer set amoA-2F-amoA-5R, followed by amoA-3F-amoA-4R, for amoA genes (Table 1). The amoA primer sets yielded noticeably lower PCR product concentrations with all DNA templates (ca. 12 ng µl1), assessed by ethidium bromide fluorescence, compared to yields obtained by the 16S rDNA primer sets (20 to 30 ng µl1).
16S rDNA DGGE, probe hybridization, and sequence analysis.
Initial analysis was carried out on individual samples of soil obtained from replicate plots of the three grassland soil sites at Sourhope. DGGE banding patterns for ammonia oxidizer 16S rRNA genes from these sites were similar (Fig. 1), with the exception of a single U4a replicate [U4a(1)]. PCR products from all sites comigrated with products from cloned standards representing Nitrosomonas clusters 6 or 7 (band I) and Nitrosospira clusters 2 and 3 or Nitrosomonas cluster 5 (band II). U4a(1) generated a single band (band III) that comigrated with the Nitrosospira cluster 1 standard. However, comigration of certain bands from different clusters and the presence of multiple bands, which are due to the degeneracy of the reverse PCR primer (16, 24) and possibly the forward primer, prevents unambiguous identification of ammonia oxidizer clusters based on migration patterns alone. It was therefore necessary to excise and reamplify specific DGGE bands of interest for direct sequencing and phylogenetic analysis (Fig. 2). Previous 16S rRNA gene studies have shown that ß-proteobacterial ammonia oxidizers are monophyletic with two genera, which can be split into seven distinct phylogenetic clusters (16, 39, 40). LogDet/Paralinear distances of variable sites produced essentially the same phylogenetic tree topology and cluster designations as previously reported (39), although the Nitrosospira cluster 4 in our analysis appeared to be paraphyletic and weakly supported, presumably due to the relatively short sequence used. Phylogenetic analysis confirmed the environmental sequences found in this study to be 16S rDNA ofß -proteobacterial ammonia oxidizers and assigned them to Nitrosomonas cluster 7 [band I, sequences U4b(3)/I and MG6(2)/I], which is closely related to that of N. europaea, Nitrosospira cluster 3 [band II, sequences U4a(2)/II and MG6(2)/II], and Nitrosospira cluster 1 from the unimproved grassland site, U4a(1) [band III, sequence U4a(1)/III] (Fig. 2). In addition to sequence analysis, band identity was confirmed by DNA hybridization with specific ammonia oxidizer genus and/or cluster specific probes (40). Samples of the PCR products shown in Fig. 1 were analyzed by dot blot DNA hybridization with ammonia oxidizer cluster-specific probes NspCL2_458 and NspCl3_454 (Table 2). Positive DNA hybridization was only observed with the Nitrosospira cluster 3 probe (data not shown) demonstrating the absence of cluster 2 organisms in these samples.
![]() View larger version (36K): [in a new window] |
FIG. 1. DGGE analysis of 16S rDNA sequences from Sourhope grassland soil sites with representative cluster controls. Lanes 1 to 7 show ß-ammonia oxidizer cluster controls: lane 1, EnvB1-8 (Nitrosospira cluster 1); lane 2, pH4.2A/27 (Nitrosospira cluster 2); lane 3, pH4.2A/4 (Nitrosospira cluster 3); lane 4, pH7B/C3 (Nitrosospira cluster 4); lane 5, EnvA1-21 (Nitrosomonas cluster 5); lane 6, EnvC1-19 (Nitrosomonas cluster 6); and lane 7, N. europaea (Nitrosomonas cluster 7). Lanes 8 to 16 show Sourhope grassland samples from triplicate plots: unimproved, U4a (lanes 8 to 10); semi-improved, U4b (lanes 11 to 13); and improved, MG6 (14 to 16). The sequences of the individual bands were obtained from the following lanes: band I, lanes 13 and 15, U4b(3) and MG6(2); band II, lanes 9 and 15, U4a(2) and MG6(2); and band III, lane 8, U4a(1)
|
![]() View larger version (30K): [in a new window] |
FIG. 2. Phylogenetic tree showing the relationship of the sequences derived from excised DGGE bands to environmental and pure culture ammonia oxidizer sequences based on 272 bases of aligned 16S rRNA gene sequences. Minimum evolution tree derived from LogDet/Paralinear distances of variable sites (estimated value of proportion of invariable sites = 0.659). The tree shown is one of six equally valid trees that showed minor rearrangements within clusters. Bootstrap support values over 50% (100 replicates) are shown in boxes. NS denotes support below 50%. Upper value, bootstrap derived by maximum parsimony; lower value, derived by LogDet/Paralinear distances of variable sites.
|
![]() View larger version (74K): [in a new window] |
FIG. 3. Repeat DGGE analysis of 16S rDNA sequences from Sourhope grassland soil sites with representative cluster controls and DIG-labeled probing. Lanes 1 to 7 in panels a to d, ß-ammonia oxidizer cluster controls (see legend for Fig. 1); lanes 8 to 16, Sourhope grassland samples from triplicate plots: unimproved, U4a (lanes 8 to 10); semi-improved, U4b (lanes 11 to 13); and improved, MG6 (lanes 14 to 16). (a) DGGE gel before membrane transfer and DNA hybridization. (b to d) Membrane hybridized with NspCL1_249 (b), NspCL3_454 (c), and Nmo254a (d).
|
![]() View larger version (81K): [in a new window] |
FIG. 4. DGGE analysis of 16S rDNA sequences derived from replicate 0.5-g soil samples taken from the same unimproved U4a(1) (A) and improved MG6(1) (B) grassland sites, with representative cluster controls. Lanes 1 to 7, ß-ammonia oxidizer cluster controls (see legend for Fig. 1); lanes 8 to 12, replicate grassland soil samples.
|
![]() View larger version (27K): [in a new window] |
FIG. 5. DGGE analysis of 16S rDNA sequences derived from replicate PCR reactions of the same 0.5-g soil DNA extract. Lanes 1 to 5, replicate PCR reactions from soil DNA sample U4a(1) replicate E; lanes 6 to 10, replicate PCR reactions from soil DNA sample MG6(1) replicate E.
|
|
View this table: [in a new window] |
TABLE 4. Identification by probe hybridization of DGGE banding patterns of replicate 0.5-g samples from the same field plot
|
![]() View larger version (58K): [in a new window] |
FIG. 6. DGGE analysis of amoA gene sequences derived from Sourhope grassland soil sites and ammonia-oxidizing bacteria pure cultures. B6, Nitrosospira sp. strain B6; 40KI,Nitrosospira sp. strain 40KI; NpAV, Nitrosospira sp. strain NpAV; eut, N. eutropha; eur, N. europaea. Arrows represent bands excised for sequence and phylogenetic analysis.
|
![]() View larger version (39K): [in a new window] |
FIG. 7. Phylogenetic tree showing the relationship of the sequences derived from excised DGGE bands to pure culture ammonia oxidizer sequences based on 186 bases of aligned amoA gene sequences. A 50% majority-rule consensus tree (of 47 trees) derived by LogDet/Paralinear distances of variable sites (estimated value of proportion of invariable sites = 0.681) is shown. The values in parentheses represent the % G+C content of the variable sites analyzed. Bootstrap support values of >50% (100 replicates) are shown in boxes; NS denotes support below 50%. First value, bootstrap derived by maximum parsimony; second value, derived by LogDet/Paralinear distances of variable sites; third value, bootstrap derived by LogDet/Paralinear distances of variable sites with amoA third codon position removed.
|
|
View this table: [in a new window] |
TABLE 5. Heterogeneity in soil pH and ammonium concentration between different 0.5-g soil samples from the same replicate field plot bulk soil samplea
|
|
|
|---|
Unimproved grassland soils had the lowest pH and the lowest net nitrification potential, in comparison with improved soils, which received greater fertilizer input. Nitrate concentrations and most-probable-number estimates of ammonia oxidizers have been shown to be lower in undisturbed soils than in agricultural soils (38). Nitrification may be undetectable in climax communities, such as permanent grasslands (19, 35) and forest soils (12), where nutrient cycling is highly effective. Ammonia oxidizer populations can vary from <500 g1 in acidic soils to ca. 105 g1 in agricultural soils (2, 13, 20), and the higher nitrification potential in improved pastures may therefore result from increases in active ammonia oxidizer populations through greater N input.
Amplification of ammonia oxidizer 16S rRNA gene sequences from DNA extracts of Sourhope soil required the use of a nested PCR system, which appears to be essential for samples with low numbers of ammonia oxidizers (9, 32). Nested PCR was also required for amplification of amoA gene sequences, and yields of PCR amplified 16S rRNA genes were generally lower in unimproved (20 ngµl1) than in improved grassland soils (30 ngµl1), possibly reflecting ammonia oxidizer cell numbers. It is also possible that differences in 16S rDNA yields may be due to differences in the levels of PCR-inhibiting compounds, such as humic acids (31, 45), although PCR amplification with the eubacterial primer set (29) resulted in the same yield of PCR product for both unimproved and improved sites (data not shown).
DGGE analysis (28) of the PCR products of both the 16S rRNA gene and the functional gene, amoA, were used to assess ammonia oxidizer species composition. Ammonia oxidizer 16S rDNA sequences from unimproved soils were more diverse, at the cluster level, than those from the improved site, and Nitrosospira cluster 1 sequences were only observed within samples from U4a.Nitrosospira cluster 1 sequences had previously only been detected in marine environments (24, 32, 39). However, cluster 1 sequences have recently been detected in acid forest soils from Finland and Scotland (T. Aarnio, C. J. Phillips, and J. I. Prosser, unpublished data; G. R. Campbell, unpublished data), although their frequency of occurrence within soil samples appears to be less than in marine environments. Nitrosospira cluster 3 and Nitrosomonas cluster 7 were present in all grassland types. Pure culture representatives of these clusters have been found in soil enrichments, and cluster 3 sequences often dominate in molecular studies of ammonia oxidizer diversity in soil (2, 9, 26, 33, 39, 40). Nitrosospira cluster 3 contains the majority of cultured ammonia oxidizers, while Nitrosomonas cluster 7 includes the most-studied species, N. europaea. Although, N. europaea was isolated from soil, no cluster 7 sequences have been detected in soil samples. However, sequences closely related to the amoA sequence of N. europaea were found in an agricultural soil that regularly received high levels of N fertilizer (50). We recovered sequences that were closely related to the amoA sequence of N. europaea, in addition to those related to Nitrosospira species. DGGE analysis of amoA sequences demonstrated differences in amoA gene diversity between unimproved and improved grass pastures, while the semi-improved site shared similarities with both.
We were able to clearly resolve differences between amoA sequences at the genus level. Species level relationships were poorly resolved due to a low level of sequence variation, most of which occurred at codon position three. There were marked differences between the base compositions of amoA sequences assigned to Nitrosomonas (43 to 57% G+C content) and Nitrosospira (67 to 83% G+C content), suggesting that they evolved under different selection pressures (14).
Although reported to be dominant in acid agricultural soils (39, 40), Nitrosospira cluster 2 sequences were not detected in initial experiments (Fig. 1, 2, and 3), possibly due to low abundance or spatial distribution of the organisms within these soils. However, probe hybridization (Table 4) identified cluster 2 sequences in DGGE bands from repeat replicate samples of the improved site, MG6(1) but not in repeated samples of U4a(1). Nitrosospira cluster 2 strains have been isolated from a sewage treatment works in Norway (46), suggesting physiological diversity within this group and possible selection by high ammonium concentrations.
Repeated DNA extractions, PCR amplifications of 16S rDNA genes, and subsequent analysis by DGGE demonstrated variation in migration patterns between replicate subsamples of unimproved soil, although consistently similar banding patterns were obtained from similar subsamples from the improved site. Nitrosospira cluster 4 sequences, not detected in earlier experiments, were also identified in one replicate sample of U4a. These results suggest that heterogeneity and diversity of ammonia oxidizer populations are greater in unimproved soils than in managed soils, possibly through selection for particular groups, as observed for
-Proteobacteria at this site (23). The greater diversity of ammonia oxidizer 16S rDNA sequences from unimproved grassland soils may reflect the physical and chemical heterogeneity of these undisturbed soils. Strong et al. (42), using 0.8-g-samples of soil, found substantial heterogeneity between samples, which was not compromised by the sampling scale. We observed, at a similar spatial scale (0.5 g), spatial heterogeneity in both soil pH and NH4+N that was lower in improved than in unimproved soil due to regular additions of nitrogen fertilizer. Both pH and ammonium concentration are important for ammonia oxidizer growth and activity, and heterogeneity in these factors may have driven heterogeneity in ammonia oxidizer populations. Reduced diversity of ammonia oxidizer populations in fertilized soils has also been observed previously by Bruns et al (2). The present study therefore demonstrates a significant impact of soil management on the heterogeneity of soil physicochemical characteristics of relevance to ammonia oxidizer growth and on the diversity and heterogeneity of ammonia oxidiser populations.
We thank Allison McCaig and Arjen Speksnijder for help with the phylogenetic analysis and Karl-Paul Witzel for the amoA primer sets.
|
|
|---|
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»