Previous Article | Next Article ![]()
Applied and Environmental Microbiology, July 2008, p. 4440-4453, Vol. 74, No. 14
0099-2240/08/$08.00+0 doi:10.1128/AEM.02692-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Department of Oceanography, Florida State University, Tallahassee, Florida 32306,1 Department of Oceanography, Texas A&M University, College Station, Texas 77843,2 Institute of Marine and Coastal Sciences, Rutgers University, New Brunswick, New Jersey 089013
Received 29 November 2007/ Accepted 12 May 2008
|
|
|---|
|
|
|---|
Nitrogen often limits primary production in marine environments (30), and the predominant sink or loss of nitrogen in the coastal ocean is nitrification-denitrification. Nitrification and denitrification are bacterially mediated processes that, when coupled, link the mineralization of nitrogenous compounds to the liberation of gaseous nitrogen that can be released to the atmosphere (15, 42, 70). Understanding the diversity of the bacterial groups associated with nitrification and denitrification is critical to understanding the factors that may influence this important part of the nitrogen cycle in the marine environment. In order to elucidate the nitrifier communities, the functional gene, amoA, was assayed which encodes the first subunit of ammonia monooxygenase, a protein involved in the first step of nitrification. This gene has been used to identify nitrifiers, and a previous study by Hunter et al. (33) reported a low diversity of amoA genes in sandy marine sediments. However, a limited number of samples were examined, and very little is known about the diversity of bacterial amoA genes in permeable marine sediments across environmental gradients. For monitoring denitrifying bacteria in marine environments, three functional genes (nirS, nirK, and nosZ) have been used to (9, 10, 43, 66, 67) due to their broad phylogenetic diversity. While nirS and nirK encode enzymes active early in the denitrification pathway, nosZ encodes the enzyme nitrous oxide reductase, which catalyzes the final step in denitrification, thus representing the process that leads to the loss of biologically available nitrogen from the sediment. Although reports of denitrifier diversity in continental shelf sediments have been published (66-68), the database for marine sediments remains small, and past studies have often focused on methodological development.
The objective of the present study was to expand the characterization of bacterial community diversity in understudied permeable shelf sediments across spatial, temporal, and hydrodynamic gradients. Simultaneous purification and analysis of community DNA and RNA compared total and metabolically active bacterial populations. High-throughput techniques were developed and coupled with clonal analysis to facilitate the determination of overall, denitrifier, and nitrifier bacterial community composition at increased spatial resolution. Whereas clonal analysis revealed no substantial changes across environmental gradients, terminal restriction fragment length polymorphism (TRFLP) profiles showed distinct differences in bacterial community composition according to site, depth, and time period sampled.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Biogeochemical characteristics of gulf and bay sediments sampled in March and May
|
Bacterial SSU rRNA, nosZ, and amoA genes were amplified by PCR. General reaction conditions included a standard PCR mix (1x PCR buffer containing 1.5 mM MgCl2 [Takara Bio, Inc., Japan], 250 µM concentrations of each deoxynucleoside triphosphate [Takara Bio], 1 to 15 pmol each of forward and reverse primers, 1 to 2 U of Taq polymerase [Takara Bio], deionized H2O, and 10 to 20 ng of DNA). The primer sequences and specific PCR conditions for DNA amplification are listed in Table 2.
|
View this table: [in a new window] |
TABLE 2. Oligonucleotide primer pairs and PCR conditions used in this study
|
PCR and reverse transcription-PCR products were cleaned by using a Qiagen PCR purification kit and cloned by using a TOPO TA cloning kit (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions. A total of 15 libraries (4 SSU rRNA gene [DNA derived], 4 SSU rRNA [RNA derived], 4 nosZ, and 3 amoA) were constructed from bay and gulf sediment samples collected from the 0- to 2-cm and 18- to 20-cm depth intervals. No amoA amplification was detected in DNA extracted from either of the three individual or combined March bay cores at 18- to 20-cm depth; thus, three amoA clone libraries were constructed. Cloned inserts were amplified by using the vector-specific primers M13F and M13R (50) and corresponding PCR conditions (Table 2). After amplicon size verification, 1 µg of the clone PCR product was digested by using HaeIII (New England Biolabs, Inc., Beverly, MA) and MspI (Promega) restriction enzymes for 2.5 h at 37°C according to the manufacturer's instructions for RFLP analysis. Digested DNA fragments were size fractionated by a 2% agarose gel stained with ethidium bromide and visualized by using a Bio-Rad GelDoc XR system (Bio-Rad Laboratories, Inc., Hercules, CA). DNA fragment sizes were estimated by comparison to molecular mass standards (1-kb and 50-bp DNA ladder; Promega). Clones were grouped into phylotypes according to banding patterns with representatives of select SSU rRNA and SSU rRNA gene phylotypes (i.e., phylotypes containing more than one clone), and all nosZ and amoA gene phylotypes were bidirectionally sequenced by using an Applied Biosystems 3100 genetic analyzer (ABI, Foster City, CA) at the Florida State University sequencing facility. Select phylotypes with banding patterns that matched previously sequenced phylotypes from Hunter et al. (33) were sequenced to verify similarity between studies. Clones with designations beginning LC or SC are from Hunter et al. (33) and represent phylotypes detected in both studies.
TRFLP analysis.
For TRFLP, the bacterial SSU rRNA and nosZ genes were PCR amplified from the same community DNA used for clone library construction. Similar PCR amplification conditions were used as described above, except the forward primers were labeled with 6-carboxyfluorescien (ABI). Fluorescent amplification of amoA was unsuccessful due to gene template concentrations being below our detection limit. TRFLP analysis of SSU rRNA and nosZ genes was performed according to the method of Kerkhof et al. (40) using MnlI as the restriction enzyme and the forward primers listed in Table 2. TRFLP fingerprinting was performed on an ABI 310 genetic analyzer (ABI) using Genescan software. To identify individual peaks, the 5' terminal restriction fragment (TRF) length for each sequenced clone was determined by an in silico digest using the MnlI recognition site and matched to the peak sizes from the TRFLP analysis. Ten phylotypes were chosen at random to have in silico digests verified by individual TRFLP analysis. For each phylotype, the in silico predictions were within one base of the size observed on TRFLP electropherograms with a majority being exact matches. In silico digests from Hunter et al. (33) were also included in the analysis of TRFLP profiles. Similar TRF lengths that were obtained from multiple unrelated clone sequences were discarded, and associated electropherogram peaks were not included in statistical calculations. TRFs from representative phylotype sequences were matched to peaks from each of the 28 sediment sample profiles (i.e., seven depths, two sites, and two months). Further, comparative analysis of the TRFLP profiles was performed based on the Sorenson's similarity index (47, 56) and UPGMA (for unweighted pair-group method with arithmetic averages) analysis using the COMbinatorial Polythetic Agglomerative Hierarchical clustering package (COMPAH96; http://www.es.umb.edu/edgwebp.htm). Dendrograms were constructed by inserting distance matrices into Clustering Calculator (www2.biology.ualberta.ca/jbrzusto/cluster.php#ClusterCalc) using Canberra Distance for distance measure and Saitou and Nei neighbor-joining for the clustering method. Output trees were visualized on TreeView (taxonomy.zoology.gla.ac.uk/rod/treeview.html).
Phylogenetic and statistical analyses.
All clone sequences were checked for chimeras using Chimera Check from the Ribosomal Database Project II Release 9 (48). Sequences from the present study and reference sequences, as determined by highest sequence similarity during BLAST analysis, were subsequently aligned by using the Fast Aligner algorithm in the ARB package (71). All SSU rRNA alignments were visually verified and manually adjusted according to Escherichia coli SSU rRNA secondary structure. According to the method of Hunter et al. (33), neighbor-joining trees incorporating a Jukes-Cantor distance correction were created from the alignments using the ARB software package (71). An average of 500 (i.e., amoA) to 1,000 (i.e., SSU rRNA and nosZ clones) nucleotides were included in the phylogenetic analyses. Bootstrap data represented 1,000 samplings. RNA derived SSU rRNA clones were not included in tree construction due to the reduced lengths of cloned sequences (
500 bp) available for comparison. Rarefaction analysis was performed using equations described by Heck et al. (27). The percent coverage (C) of the clone libraries was calculated according to the following equation: C = [1 – (ni/N)] x 100 (24, 54), where ni is the number of unique clones as determined by RFLP analysis, and N is the total number of clones in the library. Sorensen's index, the Shannon-Wiener index, and species richness were determined by EstimateS (14, 16, 17).
Nucleotide sequence accession numbers.
The 52 total nucleotide SSU rRNA, nosZ, and amoA sequences reported here were submitted to the GenBank database under accession numbers DQ431855 to DQ431906.
|
|
|---|
RFLP and statistical analyses of SSU rRNA, nosZ, and amoA clone libraries.
Fifteen clone libraries constructed by using nucleic acids extracted from March bay 0- to 2-cm (3B02) and 18- to 20-cm (3B1820) and gulf 0- to 2-cm (3G02) and 18- to 20-cm (3G1820) sediment core sections resulted in a total of 249 Bacteria SSU rRNA gene clones (DNA-derived), 136 Bacteria SSU rRNA clones (RNA-derived), 231 nosZ clones, and 64 amoA clones. The sample abbreviations given in parentheses will henceforth be used to identify the clone sequences retrieved. Clone sequences with designations beginning with LC or SC are from Hunter et al. (33) and represent phylotypes detected in both studies. The percent coverage for all four SSU rRNA gene clone libraries was 66% (Table 3), with individual library coverage ranging from 32 to 42%. RNA-derived SSU rRNA clone libraries had higher percent coverages for the individual libraries (54 to 87%) and the combined, total library (88%). All diversity indices suggested lower diversity in the RNA-derived libraries compared to the corresponding DNA-derived libraries (Table 3). The total nosZ clone library, integrating all samples, had a calculated total percent coverage of 74% with lower individual clone library percent coverages (39 to 65%) (Table 3). Statistical estimators suggested the nosZ gene diversity was similar to that of the SSU rRNA (Table 3), but much higher than that of the amoA gene. Two phylotypes were detected in the three amoA clone libraries (see Table S1 in the supplemental material) with a percent coverage of 100% for each (data not shown). Rarefaction curves for the total SSU rRNA (DNA- and RNA-derived) and nosZ libraries and the total and individual amoA gene libraries suggested a sufficient number of clones were sampled to represent library diversity (see Fig. S1 in the supplemental material). Although rarefaction curves for the remaining libraries did not indicate sample saturation, sufficient clones were screened to indicate trends within the bacterial communities.
|
View this table: [in a new window] |
TABLE 3. Statistical analyses of SSU rRNA and nosZ gene clone libraries using standard ecological and molecular estimates of sequence diversity
|
![]() View larger version (26K): [in a new window] |
FIG. 1. Phylum Proteobacteria neighbor-joining phylogenetic tree, incorporating a Jukes-Cantor distance correction, of SSU rRNA gene sequences from 3B02, 3B1820, 3G02, and 3G1820 samples. Sequences from the present study and close relatives were aligned by using the Fast Aligner algorithm, verified by hand, and compared to the E. coli SSU rRNA secondary structure using the ARB software package. Sequence names in boldface type represent both cloned sequences and sequences from Hunter et al. (33) used in TRFLP peak analysis. Bootstrap analyses were conducted on 1,000 samples, and percentages greater than 50% are indicated at the nodes. Methanococcus maripaludis was used as the outgroup. Scale bar, 0.10 change per nucleotide position. Abbreviations: Desulfobu, Desulfobulbaceae; Desulfoba.(F), Desulfobacteraceae; Desulfoba.(O), Desulfobacterales; Rhodo.(F), Rhodobacteraceae; Rhodo.(O), Rhodobacterales; Hyphomicr., Hyphomicrobiaceae; Ectothio, Ectothiorhodospiraceae; Chroma., Chromatiales.
|
![]() View larger version (28K): [in a new window] |
FIG. 2. Non-Proteobacteria neighbor-joining phylogenetic tree, incorporating a Jukes-Cantor distance correction, of SSU rRNA gene sequences from 3B02, 3B1820, 3G02, and 3G1820 samples. Sequences from the present study and close relatives were aligned by using the Fast Aligner algorithm, verified by hand, and compared to the E. coli SSU rRNA secondary structure using the ARB software package. Sequence names in boldface type represent both cloned sequences and sequences from Hunter et al. (33) used in TRFLP peak analysis. Bootstrap analyses were conducted on 1,000 samples, and percentages greater than 50% are indicated at the nodes. Methanococcus maripaludid was used as the outgroup. Scale bar, 0.10 change per nucleotide position. Abbreviations: Plancto.(F), Planctomycetaceae; Plancto.(O), Planctomycetales; Chrooc., Chroocales; Chlorob.(F), Chlorobiaceae; Chlorob.(O), Chlorobiales; Sphingo., Sphingobacteriales; Flavo.(F), Flavobacteriaceae; Flavo.(O); Flavobacteriales; Gemma.(F), Gemmatimonadaceae; Gemma.(O), Gemmatimonadales; Acido.(F), Acidobacteriaceae, Acido.(O), Acidobacteriales; Paeni., Paenibacillaceae; Turici., Turicibacteraceae; Anaero.(F), Anaerolinaceae; Anaero.(O), Anaerolinales.
|
Non-Proteobacteria lineages.
Eight non-Proteobacteria phyla were identified by phylogenetic analysis. Of the 26 non-Proteobacteria-related phylotypes, 5 grouped within the phylum Planctomycetes, which contains the Pirellula-associated anammox functional group (Fig. 2). Although two phylotypes branch within the Pirellula clade, the other three phylotypes, represented by clones 3B1820-36, LC1-9, and LC1-1, diverged to form a second clade that was supported by strong bootstrap values (Fig. 2). The Bacteroidetes/Chlorobi superphylum contained six diverse phylotypes, accounting for 14% of the total SSU rRNA gene clones (Table 4), and was related to the orders Sphingobacteriales and Chlorobiales and the families Saprospiraceae, Flexibacteraceae, and Flavobacteriaceae. All Acidobacteria- and Actinobacteria-related phylotypes had low sequence similarity to identified cultured isolates; thus, classification beyond the phyla level was not possible (Fig. 2). Five of the six phylotypes in the phylum Actinobacteria, clustered within their own clade and apart from any cultured isolates (Fig. 2), with only SC3-41 (33) being affiliated with the cultured genus "Candidatus Microthrix parvicella" (89%). A single phylotype, 3B1820-43, clustered within the phylum Cyanobacteria and was 99% similar to the chloroplast DNA sequence from Alnus incana (Fig. 2 and see Table S2 in the supplemental material). Three phylotypes—LC1-24, 3G02-01, and 3G1820-58—clustered within the family Anaerolinaceae in the phylum Chloroflexi (Fig. 2). Phylotypes 3G02-02 and 3B1820-01 clustered with the families Paenibacillaceae and Turicibacteraceae, respectively, within the phylum Firmicutes (see Table S2 in the supplemental material and Fig. 2). The final two phylotypes, 3G1820-03 and 3G1820-56, branched into the aerobic heterotrophic phylum Gemmatimonadetes (Fig. 2).
SSU rRNA (RNA-derived) phylogenetic analysis.
Sequence analysis of 23 of the 42 SSU rRNA phylotypes (RNA-derived) indicated six distinct phyla. Although Firmicutes, Bacteriodetes/Chlorobi, and Gemmatimonadetes were detected in the DNA-derived SSU rRNA gene libraries, no RNA-derived clones grouped into these phyla. Similar to the DNA-derived clones, nearly half of the RNA-derived clone sequences (41%) grouped into 12 Proteobacteria-related phylotypes; however, Gammaproteobacteria-related clones represented 59% of the total Proteobacteria (see Table S3 in the supplemental material).
Twelve RNA-derived Proteobacteria-related phylotypes grouped into three classes: Alpha-, Delta-, and Gammaproteobacteria. A majority of the DNA-derived phylotypes (n = 9) was most related to the Alphaproteobacteria (see Table S2 in the supplemental material), while only one RNA-derived phylotype grouped in this lineage. Six of the twelve Proteobacteria-related RNA-derived phylotypes grouped within the class Gammaproteobacteria. Although no Pseudomonadaceae-related clones were detected in the DNA-derived libraries, they were the most frequently detected Gammaproteobacteria lineage in the RNA-derived library (see Table S3 in the supplemental material) and were widely identified in the TRFLP profiles. Deltaproteobacteria-related phylotypes were more depth specific in RNA-derived libraries compared to the DNA-derived libraries.
A total of 11 RNA-derived non-Proteobacteria-related phylotypes grouped into five different phyla, all represented in the DNA-derived libraries. Actinobacteria, Planctomycete, and Acidobacteria-related phylotypes, incorporating 13, 12, and 5% of the total RNA-derived clones, respectively, were detected in all depths at each site (see Table S3 in the supplemental material). All three phylotypes related to phototrophic Cyanobacteria (11% of the total RNA-derived clones) were detected in the 18- to 20-cm libraries.
Phylogenetic analysis of the nosZ gene.
Sequence analysis of the 42 nosZ gene phylotypes indicated that all sequences were most closely related to sequences of uncultured organisms that putatively clustered with Alphaproteobacteria-related nosZ genes (Fig. 3) However, due to low sequence similarity to cultured strains and the presence of nosZ on mobile genetic elements (69), characterization of these genes as being Alphaproteobacteria-related is reported with caution. The most frequently detected phylotype (17% of the total nosZ clones), 3G02N-02, was 81% similar to the environmental clone CSS bacterium 696M, retrieved from permeable sediments off the coast of New Jersey (66) (see Table S4 in the supplemental material). Several deep branching clades were comprised solely of clones from the South Atlantic Bight and NEGOM (Fig. 3) (2, 66).
![]() View larger version (21K): [in a new window] |
FIG. 3. Neighbor-joining phylogenetic tree, incorporating a Jukes-Cantor distance correction, of the nosZ gene sequences from 3B02, 3B1820, 3G02, and 3G1820 samples. Sequences from the present study and close relatives were aligned by using the Fast Aligner algorithm and verified by hand by using the ARB software package. Sequence names in boldface type represent both cloned sequences and sequences from Hunter et al. (33) used in TRFLP peak analysis. Bootstrap analyses were conducted on 1,000 samples, and percentages greater than 50% are indicated at the nodes. Ralstonia eutropha strain H16 was used as the outgroup. Scale bar, 0.10 change per nucleotide position.
|
![]() View larger version (12K): [in a new window] |
FIG. 4. Replicate TRFLP fingerprints of May bay and gulf 6- to 8-cm samples. Replicates from this depth demonstrate the reproducibility of TRFLP profiles from different cores collected at the same time and from the same site. Difference in peak presence and height can be observed between sites.
|
![]() View larger version (32K): [in a new window] |
FIG. 5. Sorensen's index-based dendrograms of SSU rRNA gene (A) and nosZ gene (B). TRFLP profiles from each depth, site, and time were compared for peak absence or presence to produce a distance matrix pairwise comparison. Neighbor-joining trees were then constructed from the distance matrices.
|
![]() View larger version (28K): [in a new window] |
FIG. 6. Relative abundance of microbial populations across sediment depth, site and time. Phylogenetic groups were determined by comparing in silico digests of cloned sequences to electropherogram peak sizes during TRFLP analysis. The percentage of the total peak area was then reported for each clone and combined with similar lineages to produce abundance profiles. TRFLP peak area profiles for Gemmatimonadetes were omitted due to a single clone detected in March Gulf samples from 4 to 10 cm. Please note the differences in x-axis scales.
|
Multiple non-Proteobacteria-related phylotypes were detected and showed similar site, depth, and time specificity. Pirellula-related phylotypes were detected at mid-depths in bay and gulf samples (Fig. 6), whereas Actinobacteria-related phylotypes from bay sediments were observed higher in the sediment column compared to those detected in gulf sediments (Fig. 6). Firmicutes-related phylotypes 3G02-02 and 3B1820-01 clustered with the families Paenibacillaceae and Turicibacteraceae, respectively (see Table S2 in the supplemental material and Fig. 2) and were only detected in the shallow sediment depth TRFLP profiles (Fig. 6).
Similar to the community fingerprints of the SSU rRNA gene, approximately half of the TRFLP peaks from the nosZ fingerprints were identified and given a corresponding phylotype. Although numerous phylotypes were detected in both the bay and gulf samples including several that were detected at nearly all depths (i.e., SC3N-21/SC1N-18, SC3N-27, SC3N-19; see Fig. S3 in the supplemental material), some phylotypes were exclusive to one site or detected at limited depths within the cores. Peaks associated with phylotypes 3G1820N-06, 3G1820N-48, and SC3N-10 were detected on both sample dates but only in gulf sediments (see Fig. S3 in the supplemental material). A total of 87% of the clones detected in phylotype 3G02N-02 were isolated from gulf samples and 78% of these were in 3G02N (see Table S4 in the supplemental material). This is in agreement with TRFLP peak areas associated with this phylotype being much higher, with wider distribution through the gulf samples (see Fig. S3 in the supplemental material). More than 30% of the TRFLP peak area in the shallow bay sediment samples was associated with the phylotype 3G1820N-14 (10% of the total nosZ clones). Interestingly, this phylotype was not detected at any depth or time point in gulf sediments.
|
|
|---|
Phylogenetic analysis of the total and metabolically active bacterial communities in marine sands.
Similar to previous studies of marine sediments (12, 33, 44, 45, 55), as well as of marine bacterioplanktonic communities (72), we showed that the Gammaproteobacteria, Planctomycetes, and Bacteroidetes groups predominate over the sequence diversity of marine sands. In contrast to previous work (45, 55), we observed that the Alphaproteobacteria showed the highest relative abundance in overall and denitrifying communities. The Deltaproteobacteria and gram-positive groups (Actinobacteria and Firmicutes) also comprised a substantial fraction of the overall diversity observed in permeable sediments, in agreement with previous work (12, 33, 45, 55). The relative abundance of Alphaproteobacteria phylotypes in particular suggests that this group is well adapted to sediments undergoing continuous and rapid pore water advection as indicated by the high permeability and rapid rates of organic matter turnover that we observed in the sands of the NEGOM. The fact that a predominance of Alphaproteobacteria sequences is also observed in marine bacterioplankton communities (72) may be explained by the recent observations showing enhanced exchange of particles between permeable sediments and the overlying water column (5, 18, 31, 36). In comparison to muddy marine sediments, bacteria could more readily switch between attached and planktonic lifestyles in marine sands.
A surprising observation of our study was that Gammaproteobacteria-related phylotypes were more frequently detected in RNA-derived clone libraries, suggesting high levels of metabolic activity within this group. Alternatively, Gammaproteobacteria found in the present study may maintain a higher ribosome content per cell under nongrowth conditions, as has been demonstrated in some pure cultures (22) and seawater enrichments (21). Regardless of the cause, increased rRNA levels may permit certain Gammaproteobacteria to respond better to the rapid changes in oxygen and organic substrates that occur in marine sands. One genus with the capacity for denitrification, i.e., Marinobacter (64), was detected within both gulf and bay sediment RNA-derived libraries, and previous cultivation work suggests that this group is among the most abundant of denitrifiers in marine sediments (4, 25). Specifically targeting this lineage with quantitative PCR methods would provide valuable insight into its ecological impact.
The results from the denitrifying communities are consistent with the rRNA findings with a large sequence diversity detected. However, nosZ-based phylogenetic classification beyond the class level, i.e., Alphaproteobacteria, is tentative at best. Linking nosZ diversity with SSU rRNA diversity is difficult due to primarily partial nosZ sequences from environmental samples being predominant in the database. Although additional marine denitrifiers have been cultured (62), gaps exist in the sequence database between SSU rRNA and nosZ phylogenies. For example, while Hyphomicrobiaceae, a family with known denitrifying lineages, was frequently detected in the DNA libraries, no nosZ genes and no full genome sequences for this group are available in public databases. Thus, many nosZ sequences from the present study remain tightly clustered with previously reported clone sequences from environmental samples (33, 66) and not to cultivated species. These results underscore the need for cultivation of marine denitrifying model microorganisms that better represent the indigenous populations.
Analysis of the amoA gene showed a low diversity (n = 2) in NEGOM sediments in comparison with previous studies of marine muds (3, 13, 23, 58, 60). Bacterial amoA sequences showed a high sequence identity (98%) to sequences retrieved from the water column and muddy sediments of other coastal marine ecosystems. However, two ubiquitous phylotypes at both NEGOM sites represent a comparable level of diversity to a previous sandy sediment study (33). Therefore, while ammonia oxidation has been shown to play a critical role in the coastal ecosystem in other marine habitats (1, 8, 29, 49, 51, 57, 73), results presented here suggest that ammonia oxidizing-related taxa may be under-represented in these clone libraries, potentially due to targeting only the bacterial contribution to the ammonia oxidizing community (22a).
Profiling of permeable marine sediment communities over expanded scales using DNA fingerprinting.
An understanding of bacterial community structure and function in marine sediments has been hampered by the excessive time and cost of sample processing. Previous studies of bacterial diversity in marine sediments have relied mainly on clonal analysis and have often lacked replication or based interpretation of bacterial community structure on relatively few samples. In the present study, DNA fingerprinting provided improved replication and spatial coverage that was not achievable with clonal analysis. Consistent with past work in marine sands (55), we observed community composition at the phylum/subphylum level to be relatively stable in clone libraries, whereas we observed substantial trends with depth and between sites at the family to species level using our community fingerprinting approach. Such differences in community interpretation highlight the need for more robust replication in sampling and characterization. The results presented here are indicative of in situ conditions, not total community potential as shown by the microcosm experiments of Hunter et al. (33).
Sorensen's index-based analysis of TRFLP profiles provided evidence for bacterial niche specialization and suggested that the bacterial community diversity is affected by advective pore water flow in permeable sands. The gulf site was exposed to higher flushing rates due to increased exposure to larger waves compared to bay site sediments. Lower sediment flushing rates and a well-defined redox gradient, as determined by a clearly visible zone of black iron sulfides at a depth of
4 cm, were observed in the bay sediments. These differences may thus lead to at least some of the variations in the community structure that were observed. By correlating specific TRFLP peaks to cloned sequences, changes in community structure at each depth interval could be analyzed for potential ecological impacts. Two Alphaproteobacteria phylotypes were related to a known denitrifier, i.e., Hyphomicrobium denitrificans strain DSM 1869 (62), and were detected at all sites, depths, and time points, thus supporting the broad distribution observed in a previous lab-based sandy sediment community structure determination (33). Hyphomicrobiaceae lineages have been detected in a wide variety of environments (i.e., wastewater [19], deep-sea sediment [53, 74], and Antarctic sediment 52]) and linked not only to nitrogen cycling but also to growth on methylated sulfur compounds (52). Such metabolically diverse lineages with the capacity for denitrification may prove to be a significant contributor to geochemical processes in dynamic environments such as coastal sediments that are subjected to variable pulses of carbon and oxygen. At deeper depths, Deltaproteobacteria (i.e., Desulfobulbaceae and Desulfobacteraceae) became more relatively abundant, supporting the location of the redox boundary as observed by the presence of iron sulfides in the cores. Due to pore water flow within these sediments, it is interesting to hypothesize that the depth-specific taxa identified are surface adherent and thus less affected by water influx.
Large differences in community structure with depth were expected, and the differences in the distribution of bacterial taxa between sites at comparable depth are consistent with prior findings in coastal sands from the Mid-Atlantic Bight (67). The Sorensen's index-based analysis of the TRFLP profiles indicated as low as 50% similarity for SSU rRNA gene and 38% similarity for nosZ between the bay and gulf bacterial communities. Thus, differences in flushing rates or redox zone stability may allow different populations to become more numerous between sites. Free-living, non-surface-adherent populations (26) may be forced by pore water flow through the sediment and thus detected at multiple locations within the core. Such movement within the sediment may explain why some groups (i.e., denitrifiers and phototrophs) were detected throughout the core despite having limited geochemical ranges.
Vertical and lateral movement of sediment particles may also contribute to site-specific similarities in the bacterial community structure. Although wave action and bottom water currents can suspend and mix surface sediments (61, 65), homogenization may be limited to only the uppermost sediment layers. Surface-adherent populations would then be less mixed throughout the core compared to the effect of advective flow on free-living populations. Such mixing would explain the broad depth range of some detected lineages within the present study, including the aforementioned Deltaproteobacteria taxa. In addition, members of the Rhodobacteraceae and Cyanobacteria are known phototrophs; however, their profiles within the sediment were very different. While Rhodobacteraceae were detected in the shallow depths exclusively, Cyanobacteria-related lineages were observed at all depths.
Published ahead of print on 16 May 2008. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»