Previous Article | Next Article ![]()
Applied and Environmental Microbiology, March 2003, p. 1359-1371, Vol. 69, No. 3
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.3.1359-1371.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Department of Molecular and Cell Biology, Institute of Medical Sciences, University of Aberdeen, Foresterhill, Aberdeen AB25 2ZD, United Kingdom
Received 17 July 2002/ Accepted 26 November 2002
|
|
|---|
-proteobacterial ammonia-oxidizing bacteria were characterized by analysis of 16S rRNA genes amplified using group-specific primers by PCR- and reverse transcription-PCR amplification of 16S rDNA and RNA, respectively. Amplification products were analyzed by sequencing of clones and by denaturant gradient gel electrophoresis (DGGE). Amplification of primers specific for Planctomycetales-Verrucomicrobia and ß-proteobacterial ammonia-oxidizing bacteria generated products at all sampling sites and depths, but no product was generated using primers specific for
-proteobacterial ammonia-oxidizing bacteria. 16S rDNA DGGE banding patterns indicated complex communities of ß-proteobacterial ammonia-oxidizing bacteria in anoxic marine sediments. Phylogenetic analysis of sequences from clones and those excised from DGGE gels suggests dominance of Nitrosospira cluster 1-like organisms and of strains belonging to a novel cluster represented in dominant bands in 16S rRNA DGGE banding patterns. Their presence indicates a group of organisms closely related to recognized ß-proteobacterial ammonia-oxidizing bacteria that may be selected in anoxic environments and may be capable of anoxic ammonia oxidation. Sequence analysis of planctomycete clone libraries and sequences excised from DGGE gels also demonstrated a diverse microbial community and suggested the presence of new subdivisions, but no sequence related to recognized Anammox organisms was detected. |
|
|---|
Physiological studies demonstrate complete anoxic conversion of ammonia by Nitrosomonas europaea and Nitrosomonas eutropha with nitrite as electron acceptor (4) and by N. eutropha with NO2 as electron acceptor (48, 49). Autotrophic AOB have also been detected in temporarily and permanently anoxic rhizosphere habitats (5). Further potential for ammonia oxidation in aquatic sediments may lie in the activity of anaerobic AOB within the Planctomycetales, utilizing nitrite as the electron acceptor, as demonstrated in wastewater treatment processes (the Anammox process [22, 54]).
The application of molecular techniques, in particular analysis of 16S rRNA genes, provides new opportunities for the assessment of ammonia-oxidizing populations in aquatic sediments. Phylogenetic analysis of 16S rRNA genes of pure and mixed cultures places AOB in three groups. Organisms involved in the Anammox process belong to the order Planctomycetales (47, 54). Although planctomycetes are present in anoxic marine sediments (31, 34), no sequence with a close phylogenetic association with those of recognized Anammox organisms has yet been reported. The second group forms a deep branch within the
-proteobacteria of the class Proteobacteria and consists of only three recognized, closely related Nitrosococcus species (44, 62). The third group, including the majority of cultured strains, forms a monophyletic group within the ß-proteobacteria and comprises two clades, Nitrosomonas and Nitrosospira. Phylogenetic analysis of 16S rDNA sequences (44, 53) indicates the existence of clusters within these two clades, and their distribution has been investigated in a number of environments (3, 26, 27, 28, 36, 37, 42, 43, 52). Pure-culture representatives have been isolated for all groups (44), with the exception of Nitrosospira cluster 1 and Nitrosomonas cluster 5, for which only clone sequences are available from marine environments (3, 37, 42).
Molecular analysis of ammonia oxidizer 16S rDNA fragments, amplified from environmental DNA, by denaturing gradient gel electrophoresis (DGGE) characterizes community structure, enables rapid analysis of clone libraries, and, through excision, reamplification, and sequencing of bands from gels, provides information on species composition. This approach has demonstrated differences in populations of Nitrosomonas-like ß-proteobacterial AOB associated with polluted marine fish farm sediments (37) and marine aggregates (43) and domination by Nitrosospira-like organisms in Arctic Ocean waters (3). Analysis of 16S rDNA does not discriminate between metabolically active and quiescent cells. Analysis of environmental 16S rRNA, using reverse transcription-PCR (RT-PCR), provides greater sensitivity, because of the higher target copy number, and may indicate which members of the community are more metabolically active, if active cells contain larger numbers of ribosomes (13, 27, 58). The aim of this study was to determine the distribution and community composition of all three groups of AOB within an anoxic marine sediment system showing strong evidence of anoxic nitrification (40). As a first step in assessing their potential contribution to anoxic ammonia oxidation, 16S rRNA genes were amplified from DNA and RNA and their products analyzed by DGGE analysis and sequencing, with detailed phylogenetic analysis of almost complete 16S rRNA gene sequences from clone libraries.
|
|
|---|
DNA and RNA extraction.
The water content of sediments was reduced to ca. 5% (wt/vol) of the total water content by centrifugation at 16,000 x g for 10 min. Nucleic acids were extracted from 0.5 g of compacted sediment by disruption with a Ribolyser cell disruptor (Hybaid, Ltd., Ashford, United Kingdom) (speed 4, two 10-s bursts) in 2-ml screw cap Blue Matrix Ribolyser tubes (Hybaid) containing a mixture of ceramic and silica beads, 0.5 ml of hexadecyltrimethylammonium bromide (CTAB; Sigma, Poole, United Kingdom), extraction buffer, and 0.5 ml of phenol-chloroform-isoamyl alcohol (25:24:1) (pH 8.0). Preparation of CTAB buffer and further extraction and precipitation of nucleic acids were carried out as described by Griffiths et al. (17). Extracted nucleic acids were resuspended in 50 µl of RNase-free sterile water and divided into two aliquots for the preparation of DNA or RNA templates. DNA templates were further purified using 0.5-ml of Vivaspin concentrator columns (Vivascience, Sartorius, Germany).
PCR amplification of environmental 16S rDNA and rRNA fragments.
Extracted DNA and cDNA from RT reactions (see below) were amplified by nested PCR for DGGE analysis using 16S rDNA primer sets specific for ß-proteobacterial AOB (CTO189f-GC and CTO654r [28]), the NOC1-45f and NOC2-1168r primer set specific for
-proteobacterial AOB (57, 61), and the Planctomycetales-Verrucomicrobia assay (PV assay [8]) using the PLA40f primer in combination with the general eubacterial primer pf1053r (11). Alignments of PLA40f and pf1053r with sequences of recognized Anammox strains revealed no mismatches. However, mismatches in the NOC primer set (62), following alignments of
-proteobacterial AOB from the GenBank database, led to modification of primer NOC1-45f to CGTYGGAATCTGGCCTCTAGA and of NOC2-1168r to AGATTAGCTCCGCATCGCTG (modifications are shown in bold). PCR products were reamplified using the general eubacterial 357f-GC and 518r primers (41). Prior to secondary amplification, PCR products were diluted 1:50 to prevent the amplification of nontarget bacteria without affecting DGGE community profiles. All specific primer sets were tested for stringency against Planctomyces maris AL (NCIMB 2232), Planctomyces brasiliensis VL32 (NCIMB 13185), Nitrosococcus oceanus (NCIMB 11848), Nitrosomonas europaea (NCIMB 11850), Nitrosospira multiformis (NCIMB 11849), and Escherichia coli (NCIMB 86).
RT-PCR was done by the method of Griffiths et al. (17) with some modifications. DNA-free RNA for RT-PCR analysis was obtained by treating 8 µl of extracted RNA with 2 U of RQ1 RNase-free DNase (Promega Corp., Madison, Wis.) for 30 min at 37°C as specified by the manufacturer. Prior to RT, the RNA secondary structure was melted by incubating the RNA samples with reverse 16S rRNA PCR primer (see above) at 0.9 µM at 70°C for 10 min. Samples of annealed primer-template were snapped on ice, and 8 µl of RT reaction mixture was added as specified by the manufacturer (SuperScript RT RNase H-Reverse Transcriptase; GibcoBRL). RT was carried out at 42°C for 50 min, and the enzyme was subsequently heat-inactivated for 15 min at 70°C. Performing PCR on RNA samples after DNase treatment confirmed the amplification of RNA templates free of DNA contamination. Following transcription of 16S rRNA templates into cDNA, PCR was carried out as described above for DNA templates. PCR-amplified fragments were resolved by electrophoresis on 1% (wt/vol) horizontal agarose gels with visualization of DNA by ethidium bromide fluorescence and with Hyperladder DNA fragment size markers (Bioline). PCR amplification was also carried out on pure cultures or clones representative of each recognized cluster of ß-proteobacterial AOB (53) for use as reference sequences in CTO189f-GC plus CTO654r 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).
Cloning inserts of 1.1 and 1.4 kb specific for ß-proteobacterial AOB and Planctomycetales-Verrucomicrobia, respectively, were created by direct PCR amplification of DNA extracts of sample "core 8/depth 20" with the ßAMO161f and ßAMO1301r primer set, which is selective but not completely specific for ß-proteobacterial AOB (38), or by the PV assay, in conjunction with the general eubacterial primer 1492r (8, 11). PCR fragments were purified by standard agarose electrophoresis, and excised DNA bands were cleaned using QIAquick columns (Qiagen, Hilden, Germany). Primer pairs and PCR conditions for all amplifications used are summarized in Table 1. Individual reagents and their concentrations were as follows: 1x PCR buffer (Bioline), 1 U of BioTaq DNA polymerase (Bioline), 1.5 mM MgCl2, each primer at a concentration of 0.2 µM, and each deoxynucleoside triphosphate dNTP (Bioline) at a concentration of 250 µM. Amplification was performed with a total volume of 50 µl using an Omn-E thermal cycler (Hybaid), applying the thermal cycle conditions shown in Table 1.
|
View this table: [in a new window] |
TABLE 1. Primers and PCR amplification conditions
|
Cloning and sequence analysis.
Three purified PCR products resulting from amplification with primers ßAMOf plus ßAMOr or PLA40f plus 1492r were pooled to minimize PCR drift (59), ligated into the pGEM T-vector system (Promega Ltd., Southampton, United Kingdom), and transformed into XL1-Blue MRF Kan supercompetent E. coli cells (Stratagene Inc., Maidstone, United Kingdom) as specified by the manufacturer. Transformed colonies were screened for inserts of the correct size by PCR amplification with the specific primers described above. For each PCR amplicon, 150 clones with the correct fragment size insert were subjected to nested PCR with the primer sets described above and the 357f-GC plus 518r PCR products were screened by DGGE for sequence differences within the hypervariable V3 region. For ß-proteobacterial AOB, clone libraries selected for full-length sequencing of the entire insert represent all recognized clones with different DGGE migration patterns. Ten replicate clones of the recognized main ß-proteobacterial sequences showing identical migration patterns on the DGGE gels were also selected for sequencing. Clones from libraries generated by the PV assay were selected at random for full-length sequencing of the entire insert. The selected clones were amplified with vector primers M13f and M13r (Promega Ltd.), and M13 PCR products were purified by standard preparative agarose gel electrophoresis as described above. To ensure complete full-length sequence read without any ambiguities of sequences of potentially closely related ß-proteobacterial AOB, the M13 vector products were sequenced in both reading directions using the SP6 and T7 vector primers (Promega Ltd.), the internal 16S rDNA primer CTO654r, and 518f (reverse complement to the above described 518r). M13 vector products derived from the PV assay were sequenced using vector primers SP6 and T7 and the internal 16S rDNA primer pf1053r.
Individual DGGE bands were excised, and acrylamide slices were crushed and resuspended overnight at 4°C in 30 µl of sterile water to elute the DNA. Reamplified PCR products were again subjected to DGGE analysis to ensure purity and correct migration within the gels. Reamplified PCR products showing single bands during DGGE analysis were purified as described above and sequenced with the 518r primer. Sequencing reactions were performed using the BigDye Terminator cycle-sequencing kit (PE Biosystems, Warrington, United Kingdom), and the cycle-sequencing products were analyzed with a model ABI377 automated sequencer (PE Biosystems).
Sequence assembly and manual refinement of alignments were carried out using the Sequencer 4.1 program (Genes Codes Corp., Ann Arbor, Mich.) or the multiple-alignment algorithms as implemented in Clustalx 1.81 (EMBL, Heidelberg, Germany). For phylogenetic analysis, 1.1- and 1.4-kb clone sequences for ß-proteobacterial AOB and Planctomycetales-Verrucomicrobia, respectively, were aligned with closely related 16S rDNA sequences retrieved from the GenBank database at the National Center for Biotechnology Information (35) using the Basic Local Alignment Search Tool (BLASTn [2]). Mismatches between clone sequences, sequences derived from DGGE bands, and sequences retrieved from GenBank were analyzed using the B12Seq algorithm as implemented in BLASTn, and DNA-DNA sequence similarity percentages were calculated using the sequence similarity matrix function as implemented in BioEdit (freeware by Tom Hall, Department of Microbiology, University of North Carolina). Sequence inserts that were present in a minority of aligned sequences were omitted from further analysis. Phylogenetic relationships between pairs of 16S rDNA clone sequences were calculated using maximum parsimony (Phylip 3.62) and distance analysis as implemented in PAUP 4.1 (55). Maximum-parsimony analysis was carried out using the global rearrangement criterion and randomized input of sequence order with 1,000 times bootstrap resampling. A maximum-likelihood analysis of the most parsimonious tree, constructed from a subset of sequences representing the major sequence groups, was used to estimate the proportion of variable sites needed for the analysis of LogDet/Paralinear distances (29) of variable sites (32). LogDet/Paralinear distances were calculated using the minimal evolution criterion, and data were bootstrapped 1,000 times. Phylogenetic trees were constructed from distance analysis using the neighbor-joining method (46). Additional distance matrix corrections (Kimura2; Jukes and Cantor) were used to infer established tree topologies. To exclude chimeric 16S rDNA primary structures of clone sequences prior to phylogenetic analysis, sequences were analyzed using CHIMERA CHECK of the Ribosomal Database Project (RDP [35]), and the terminal 500 to 600 nucleotide sequence positions of the 5' and 3' ends were used in separate phylogenetic analyses, conducted as described above.
Nucleotide sequence accession numbers.
The partial clone sequences determined in this study have been deposited in the GenBank database under accession no. AY114311 to AY114349.
|
|
|---|
16S rDNA and rRNA DGGE analysis.
DGGE banding patterns of CTO189f-GC plus CTO654r 16S rDNA PCR products, specific for ß-proteobacterial AOB, were similar for all replicates and sampling depths. Two major bands were present in all samples (Fig. 1), with the upper band migrating to a position between Nitrosomonas cluster 5 and 6 controls and the lower band migrating to a position between Nitrosospira cluster 1 and 4 controls. However, PCR products generated using CTO primers did not comigrate with any of the control sequences representing recognized ß-proteobacterial AOB clusters. In addition, the presence of multiple bands, resulting from degeneracy of the CTO primer set, hampered reliable identification of closely migrating bands (28, 63). DGGE analysis was then carried out on nested 160-bp 16S rDNA PCR products, generated using primers 357f-GC plus 518r and spanning only the hypervariable V3 region. These PCR products were shorter, and their denaturing characteristics were more suitable for DGGE analysis. Banding patterns proved to be more discriminative and demonstrated diverse communities of ß-proteobacterial AOB. More than 15 distinct DGGE bands were obtained per sample, five of which were dominant (Fig. 2), comprising three closely migrating bands of high intensity (AI, AII, and AIII) and two bands of lower intensity (AIV and AV). DGGE banding profiles provided no evidence of major changes in community structure with sediment depth. Dominant bands were present in all samples at all depths, and variation of the banding pattern was apparent only through the presence or absence of less prominent bands. Variation in the intensity of major DGGE bands, indicating changes in relative abundance, were less pronounced than in CTO-derived profiles, probably due to saturation of PCR templates in the nested PCR approach.
![]() View larger version (37K): [in a new window] |
FIG. 1. DGGE analysis of direct-amplification 465-bp ß-proteobacterial AOB 16S rDNA sequences amplified from DNA extracted from Loch Duich marine sediments and representative cluster controls. Lanes: CL1, EnvB1-8 (Nitrosospira cluster 1); CL2, pH4.2A/27 (Nitrosospira cluster 2); CL3, pH4.2A/4 (Nitrosospira cluster 3); CL4, pH7B/C3 (Nitrosospira cluster 4); CL5, EnvA1-21 (Nitrosomonas cluster 5); CL6, EnvC1-19 (Nitrosomonas cluster 6); CL7, N. europaea (Nitrosomonas cluster 7). The right-hand lanes show a depth profile of Loch Duich sediment samples.
|
![]() View larger version (86K): [in a new window] |
FIG. 2. DGGE analysis of nested 160-bp ß-proteobacterial ammonia oxidizer 16S rDNA sequences amplified from DNA extracted from Loch Duich marine sediments and reference sequences retrieved from Loch Duich sediments by cloning and by sequencing of excised DGGE bands. Lanes: EL, environmental library (DGGE sequences excised and reamplified from 357f-GC plus 518r 16S rDNA PCR products); CL, clone library (ßAMO161f plus ßAMO1301 clone sequences retrieved from Loch Duich sediments and reamplified with 357f-GC plus 518r for DGGE analysis); 1 to 40, depth profile of Loch Duich 16S rDNA 357f-GC plus 518r PCR products. Clone sequences shown in lane CL are (from top to bottom) LD1-B6, LD1-B37, LD1-A40, LD1-A3, LD1-A1, LD1-A2, LD1-B20, LD1-B10, LD1-B28, LD1-A10, LD1-A15, LD1-B9, and LD1-A20; DGGE excised sequences shown in lane EL are (from top to bottom) LD1-env-A12, LD1-env-B12, LD1-env-B11, LD1-env-B5, LD1-env-A8, LD1-env-A5, LD1-env-A1, and LD1-env-A2. Arrows (right) mark the prominent bands as referred to in the text. No sequence was generated from band X when it was excised and reamplified.
|
![]() View larger version (90K): [in a new window] |
FIG. 3. DGGE analysis of nested 160-bp RT-PCR products derived from ß-proteobacterial ammonia oxidizer 16S rRNA extracted from Loch Duich marine sediments and reference sequences retrieved from Loch Duich sediments by cloning and excising bands from DGGE bands. Lane designations are as in the legend to Fig. 2. Arrows (right) mark the prominent bands as referred to in text. No sequence was generated from band X when it was excised and reamplified.
|
![]() View larger version (90K): [in a new window] |
FIG. 4. DGGE analysis of nested 160-bp Planctomycetales-Verrucomicrobia 16S rDNA sequences amplified from DNA extracted from Loch Duich marine sediments and reference sequences retrieved from Loch Duich sediments by cloning and sequencing of excised bands from DGGE gels. Lanes EL, environmental library (DGGE sequences excised and reamplified from 357f-GC plus 518r 16S rDNA PCR products); CL, clone library (PV clone sequences retrieved from Loch Duich sediments, reamplified with 357f-GC plus 518r for DGGE analysis); 1 to 36, depth profile 16S rDNA PCR products of Loch Duich sediment samples. Clone sequences in lane CL are (from top to bottom) LD1-PA42, LD1-PA15, LD1-PB1, LD1-PA13, LD1-PA26, LD1-PB2, LD1-PA40, LD1-PA32, LD1-PA16, LD1-PB20, LD1-PA30, LD1-PA20, LD1-PA21, LD1-PB12, LD1-PA11, LD1-PB3, LD1-PA50, LD1-PA34, and LD1-27; DGGE excised sequences in lane EL are (from top to bottom) LD1-env-P8, LD1-env-P7, LD1-env-P6, LD1-env-P1, LD1-env-P4, LD1-env-P3, LD1-env-P2, and LD1-env-P5. No sequence was generated from the band marked with an arrow when it was excised and reamplified.
|
![]() View larger version (101K): [in a new window] |
FIG. 5. DGGE analysis of nested 160-bp Planctomycetales-Verrucomicrobia 16S rDNA products amplified from RNA extracted from Loch Duich marine sediments and reference sequences retrieved from Loch Duich sediments by cloning and sequencing excised bands from DGGE gels. Lane designations are as in the legend to Fig. 4. No sequence was generated from the band marked with an arrow when it was excised and reamplified.
|
The sequence similarity of all ß-proteobacterial AOB clones generated in this study was high, on average 98.6%, indicating closely related clones. Phylogenetic and 16S rDNA sequence similarity analysis showed that the clones fell within two major groups. Group 1 contained closely related clones with sequence similarities between 99.7 and 97.9%, divided into five subgroups (Fig. 6). Sequence similarities within subgroups containing two or more clones were greater than 99.4%, and differences may have resulted from sequencing or cloning errors. However, the maximum average similarity among the five subgroups was 98.1 to 98.5%, indicating independent sequences. Group 1 clones fell within Nitrosospira cluster 1, representatives of which have previously been amplified from marine environments (28, 37, 43). The phylogenetic discrimination of cluster 1 within the nitrosospiras was better supported in all tree topologies than that between all other Nitrosospira subclusters. The closest recognized AOB relative to cluster 1 clone LD1-A1 was Nitrosospira sp. strain 04 (97.5% sequence similarity).
![]() View larger version (41K): [in a new window] |
FIG. 6. Evolutionary distance dendrogram showing the positions of environmental 16S rDNA 1.1-kb clone sequences (labels in bold), recovered from anoxic marine sediments, in relation to representative members within the ß-proteobacterial ammonia oxidizers. The tree is based on results of neighbor-joining analysis of LogDet/Paralinear distances. The multifurcation connects branches for which a relative order could not be determined unambiguously in the majority of bootstrap resamplings (1,000 replicates). Conservative bootstrap values of maximum parsimony or LogDet/Paralinear distances are indicated. The scale bar indicates an estimated 0.05 change per nucleotide position. Cluster designation are as suggested by Purkhold et al. (44) and Stephen et al. (53). GenBank accession numbers of closely related clone sequences are shown in italics.
|
|
View this table: [in a new window] |
TABLE 2. Similarities of 16S rDNA sequences from excised DGGE bands representative of ß-proteobacterial AOB and Planctomycetales-Verrucomicrobia 16S rDNA sequences (160 bp) to sequences retrieved from GenBank or clone sequences generated in this study
|
|
View this table: [in a new window] |
TABLE 3. Closest-match similarities of 1.1-kb 16S rDNA clone sequences of proposed novel ß-proteobacterial ammonia oxidizer cluster with closest-match sequences retrieved from GenBank.
|
Phylogenetic analysis of all sequences necessitated the restriction of alignments with 1.0-kb sequence fragments to the closest matches from GenBank. Phylogenetic analysis of 1.4- and 1.0-kb alignments yielded stable lineage or division association of clone sequences in all tree topologies inferred, even when using only partial sequences (the first or last 0.5 kb). The majority of PV clone sequences were associated with subdivision 5 of the Verrucomicrobia division or with bacterial candidate division WC1 and proposed bacterial candidate division BRC1 (Fig. 7). Clones LD1-PA32 and LD1-PA38 (99.3% sequence similarity to each other) may represent a new subdivision within the Verrucomicrobia since the closest match (clone wb1 A08) had only 80.2% sequence identity and phylogenetic analysis demonstrated no association with previously recognized subdivisions. Clones LD1-PA25 and LD1-PA42 may represent a previously unrecognized subdivision within Chlamydia, since both clones (99.3% sequence similarity to each other) were equally distant from recognized Chlamydia and Parachlamydia (87% sequence similarity). Of all analyzed clones, only LD1-P40 was associated with the Planctomycetales division. The closest match to this clone was marine seep clone BPC016 (89.0% sequence similarity), and the closest match to a recognized planctomycete was 87.8% to Planctomycetales sp. strain 116. Anammox-related sequences were grouped as a stable cluster in all tree topologies inferred (even if the helix 9 insertion was excluded for phylogenetic analyses), but in all tree topologies, none of the clones generated in this study was associated with the Anammox cluster (closest similarity, 77.5% for clone LD1-PA13).
![]() View larger version (43K): [in a new window] |
FIG.7. Evolutionary distance dendrogram showing the positions of environmental 16S rDNA 1.4- and 1.0-kb clone sequences (labels in bold) recovered from anoxic marine sediments, in relation to recognized and associated Anammox sequences (47), representative members of the order Planctomycetales, the divisions Verrucomicrobia and Chlamydia, and closely related bacterial lineages. The tree is based on results of neighbor-joining analysis of LogDet/Paralinear distances. The multifurcation connects branches for which a relative order could not be determined unambiguously in the majority of bootstrap resamplings (1,000 replicates). Conservative bootstrap values of maximum parsimony or LogDet/Paralinear distances are indicated. Subdivision (SD5) and bacterial lineage candidate designations (BRC1, WS3) are suggested by Hugenholtz et al. (20), Doika et al. (9), and Derakshani et al. (8). GenBank accession numbers of closely related clone sequences are shown in italics. The scale bar indicates an estimated 0.05 change per nucleotide position.
|
|
|
|---|
-proteobacterial ammonia oxidizers, but they are known to be present in marine environments (6, 62). Primers designed to be specific for these groups were used to amplify 16S rRNA genes from sediment samples, with subsequent analysis by cloning, sequencing, and DGGE. No
-proteobacterial sequence was obtained. This group is represented by a small number of pure cultures isolated from marine environments, and few studies have attempted to retrieve representative sequences from natural environments (62). The inability to detect sequences from this group may be due to inadequacy of the primer set used, low population levels, or biases associated with the DNA extraction procedure. However, since
-proteobacterial ammonia oxidizers are assumed to be ubiquitous in marine planktonic systems (62), their absence from the investigated sediments may indicate that communities of ß-proteobacterial AOB are adapted to the geochemical conditions within anoxic marine sediments rather than representing relics of marine sedimentation processes. Nevertheless, this study provided no evidence for the potential role of
-proteobacterial AOBs in ammonia oxidation in anoxic sediments. Furthermore, Nold et al. (42), analyzing amoA gene sequences, which discriminate between ß- and
-proteobacterial AOB, found no evidence for the presence of
-proteobacterial AOB in northwest Pacific marine sediments. Sequences representative of ß-proteobacterial AOB and planctomycetes were obtained from all samples and are discussed below.
ß-Proteobacterial ammonia oxidizers.
DGGE analysis, sequencing, and phylogenetic analysis of cloned 16S rRNA genes demonstrated the presence of two groups of ß-proteobacterial ammonia oxidizer communities, each with significant sequence diversity. Sequences belonging to the first group, Nitrosospira cluster 1, have previously been retrieved from other marine environments (3, 28, 43, 53), but no laboratory isolate representative of this group has yet been obtained. The second group contained sequences of particular interest (LD1-B6 and LD1-B37), which formed a novel cluster, not previously demonstrated, within the ß-proteobacterial ammonia oxidizers together with sequences obtained from ammonia-oxidizing isolates from seawater (TT140-098-2 and TT140-098-A) and marine sediments (C17) (57, 60). Isolates TT140-098-2, TT140-098-A, and C17 have been described as rod shaped and showed immunological reactivity with both Nitrosococcus oceani and Nitrosomonas marina antisera (C17) or no pronounced reactivity with either serum (TT140-098-2 and TT140-098-A) (60). Additionally, ammonia oxidation has recently been demonstrated in an estuarine isolate with high sequence affinity to the clone sequences falling within this clade (U. Purkhold, M. Wagner, and H.-P. Koops, personal communication). The novel cluster has low sequence similarity to currently recognized AOB clusters, forms a deeply branching clade based on 1.1 kb of 16S rDNA sequence, and may represent a new species. The association of LD1-B6 and LD1-B37 with isolates for which nitrifying activity has been demonstrated supports the assumption that the organisms from which these sequences are derived are ammonia oxidizers. Additionally, amoA gene sequences from marine sediments that showed no clear association with Nitrosomonas or Nitrosospira lineages (42) may belong to representatives of this novel cluster.
Inferring the AOB tree topology, based on 1.1-kb 16S rDNA sequences (382 informative base pairs), with the Log/Det Paralinear distance function yielded the most conservative bootstrap values of all methods applied, and an unambiguous association of the novel cluster, clones of Nitrosomonas cluster 5, and N. cryotolerans (retrieved only from marine environments) with the Nitrosomonas or Nitrosospira clusters was not possible. The uncertain association of N. cryotolerans within the Nitrosomonas clade is supported by phylogenetic analysis of full-length 16S rDNA and amoA gene sequences (25, 44) and its characteristic morphological features (23). Our experience indicates that the association of the N. cryotolerans cluster with Nitrosomonas is strongly influenced by the inclusion or exclusion of the new cluster sequences within the analysis. We suggest a review of the designation of N. cryotolerans within the genus Nitrosomonas once the available cultured representatives of the novel cluster have been characterized in detail. Furthermore, our results indicate that Nitrosomonas cluster 5 may represent a new species and is not closely associated with other Nitrosomonas clusters (36). However, support for this hypothesis requires analysis of full-length gene sequences from cultured representatives of all proposed clusters and additional information from DNA-DNA hybridization and physiological studies.
Environmental relevance of clone sequences.
One source of error in 16S rRNA-based analyses of microbial communities is that of sequence errors introduced during cloning and PCR (51). To assess the significance of this problem, DGGE analysis was carried out both on PCR products of environmental DNA and on clones generated by a different primer approach. In general, there was good correspondence between profiles from the two approaches, indicating that clone libraries represented the natural population well. Almost all clone sequences generated in this study have been confirmed by the two different primer approaches and by their presence as DGGE bands in PCR amplicons from environmental samples. In this sense, we question the relevance of closely related clone sequences that were not represented by DGGE migration patterns in any of the environmental samples (LD1-A15, LD1-B9, and LD1-A20).
Planctomycete sequences.
DGGE profiles and sequence analysis of amplified planctomycete 16S rRNA genes indicated high diversity within this group, but no sequence fell within the currently recognized Anammox group and all sequences showed low similarity to Anammox organism sequences and Anammox-associated clone sequences (47) in the database. There was therefore no evidence for their role in anaerobic ammonia oxidation in these sediments. The inability to detect Anammox sequences may be due to biases associated with DNA extraction and PCR amplification or to a low relative abundance within planctomycete sequences amplified using the PV assay. Additionally, only 26 of >40 different clone types indicated in DGGE profiles and only 8 of the environmental DGGE bands were sequenced. A nested PCR approach would provide the sensitivity to detect low copy numbers of Anammox-related sequences in natural environments in future studies. The majority of studies on Anammox organisms have involved laboratory enrichments or systems related to wastewater treatment plants (12, 22, 54). The relatively few environments from which enrichment cultures have been obtained, the phylogenetic depth within these lineages, and the assumption that potential Anammox planctomycetes adapted to marine sediments may be evolutionarily distant, limit the assessment of the true extent of the phylogenetic distribution of Anammox organisms within the order Planctomycetales. This is the first study targeting planctomycete clone sequences from marine environments. The PV assay was designed by Derakshani et al. (8) on the basis of signature nucleotides described as indicative of Planctomycetales- and Verrucomicrobia-like 16S rDNA sequences for which a common ancestry had been assumed (30). Using this assay, Derakshani et al. (8) isolated three clones from anoxic soil systems that were associated with the Anammox cluster but also showed a great evolutionary distance (PBS-114, PBS-III-2, and PBS-II-37 [Fig. 7]). The majority of sequences retrieved in this study were associated with the Verrucomicrobia and the deep-branching bacterial lineages WS3 and BRC1, supporting them as candidate divisions. However, it should be emphasized that the PV assay enabled the detection of previously unrecognized sequences that had evaded detection with general bacterial primers in studies of comparable environments (31, 34). None of the clone or DGGE sequences could be associated with high sequence similarity to previously described clone sequences or sequences of cultured representatives, and almost all clone sequences retrieved as closest-match sequences from GenBank have been isolated from anoxic or marine environments, suggesting their ecological relevance in these habitats.
Presence of ß-proteobacterial AOB and potential activity.
The detection of 16S rRNA gene sequences amplified from DNA provides strong evidence for the presence of ß-proteobacterial AOB in anoxic marine sediments but does not necessarily demonstrate their activity. However, sediments were sampled to a depth of 40 cm, and bioturbation and mixing are believed to be minimal in these systems (18, 40). The sedimentation rate has been estimated to be 4 mm year-1, and sediment at 40 cm deep is estimated to be approximately 100 years old (18). The presence of ammonia oxidizers at this depth indicates either their survival for prolonged periods under anoxic conditions or their ability to grow during this period. Analysis of 16S rRNA gene sequences derived from 16S rRNA may provide information about which components of the population were growing or active, if ribosome levels increase with growth rate and if ribosome turnover is rapid following cessation of growth (10). There are indications that rRNA turnover in ammonia oxidizers is slow following substrate exhaustion (58), and it has been suggested that rRNA-based activity studies might not be suitable for slow-growing Anammox- and Nitrosomonas-like AOB (47, 50, 58). Nevertheless, the AOB in the anoxic Loch Duich sediments apparently maintain high 16S rRNA levels for decades, and significant differences were found between DGGE profiles of DNA- and RNA-derived PCR products. In particular, DNA- and RNA-derived PCR products of sequences associated with the new cluster differed, suggesting either differences in growth and activity of closely related AOB or alternative strategies for rRNA synthesis and turnover. This, together with the diversity of both DNA- and RNA-derived gene sequences, suggests spatial and temporal heterogeneity in growth and activity within the community, presumably resulting from heterogeneity in physicochemical characteristics within the sediments and the existence of microenvironments.
In conclusion, this study demonstrates the potential for ammonia oxidation in anoxic marine sediments. The presence of both ß-proteobacterial AOB 16S rDNA and rRNA suggests that anoxic metabolism may enable their growth, activity, and survival in such environments. There was no evidence for the presence or activity of
-proteobacterial ammonia oxidizers or Anammox organisms in these sediments, although this may indicate limitations of the methodological approaches used.
Established pathways for anaerobic oxidation by ß-proteobacterial ammonia oxidizers and planctomycetes require substrates (NO2 or NO2-) derived from aerobic processes, which are unlikely to penetrate to the depths sampled in this study. Luther et al. (33) and Hulth et al. (21), however, provided indirect evidence for the potential for coupled anoxic nitrification-denitrification and manganese oxidation-reduction in marine sediments whereby Mn2+ was oxidized by nitrate and ammonia was oxidized by manganous oxides. Mortimer et al. (40) suggested that these processes occur throughout the iron reduction zone of the Loch Duich sediments and that the observed pronounced nitrate peaks may result from locally uncoupled nitrification-denitrification processes or locally increased concentrations of manganous oxides. The demonstrated presence and activity of ammonia oxidizers in these environments add further evidence for coupled anoxic nitrification-manganese reduction processes and have significant implications for descriptions and predictions of nitrogen turnover and geochemical cycling in marine sediments.
We thank P. Hayes and I. Davies (Marine Laboratory, Aberdeen), R. Mortimer, M. Krom, and S. Harris (University of Leeds), and J. Barnes (University of Newcastle) for helpful discussions and collection of samples.
|
|
|---|
-proteobacterium, detected by PCR and sequencing of 16S rRNA and amoA genes. Appl. Environ. Microbiol. 68:4153-4157.
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»