Previous Article | Next Article ![]()
Applied and Environmental Microbiology, February 2004, p. 781-789, Vol. 70, No. 2
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.2.781-789.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Department of Marine Sciences, University of Georgia, Athens, Georgia 30602-3636
Received 2 July 2003/ Accepted 17 November 2003
|
|
|---|
|
|
|---|
The physical environments of the two polar oceans differ in a number of important ways. The Arctic Ocean is perennially ice covered; is surrounded by continents with wide, shallow continental shelves; and receives
10% of the freshwater flowing into the world ocean. Freshwater inflow, brine exclusion during sea ice formation, the annual cycle of ice formation and melting, and the effective absence of wind mixing result in a very stable, highly stratified water column consisting of three major water masses (1, 2). The upper layer (to a depth of
60 m) consists of relatively well-mixed, low-salinity (<34
) water with a seasonal, shallow halocline. The intermediate layer (between 75 and 150 m) contains high-salinity (>34
) cold (<-1.5°C) water separated from the surface layer by an extremely stable, essentially permanent halocline. Relatively warm (2 to 4°C) Atlantic Ocean water is found below this intermediate layer. In contrast, the Southern Ocean surrounds a continent (Antarctica) with narrow, deep continental shelves and is bounded by an oceanic front. The water column is not strongly stratified. Coastal waters stratify during periods of meltwater formation when the pack ice retreats during spring and summer. The Southern Ocean receives virtually no freshwater inflow, no terrigenous organic matter, and only small amounts of terrestrially derived micronutrients (e.g., Fe) (20, 23, 42).
Archaea is divided into two kingdoms, the Euryarchaeota and the Crenarchaeota, with the possibility of a third kingdom, the Korarchaeota (6, 39). Cultivation-based analysis of the distribution of Archaea led to the belief that Archaea were restricted to extreme environments. Recent studies based on analyses of 16S rRNA gene sequences have demonstrated that Archaea are much more diverse and widespread than previously suspected (11, 18, 24). Karner et al. (19) concluded that there are 1.3 x 1028 archaeal cells (of which
20% are Crenarchaeota) and 3.1 x 1028 bacterial cells in the world ocean. Most Crenarchaeota 16S rRNA gene sequences from the marine environment are closely related and are identified as group I Crenarchaeota. The other predominantly marine group is affiliated with the Euryarchaeota and is identified as marine group II Euryarchaeota. This group has greater phylogenetic diversity among its sequenced representatives than the marine group I Crenarchaeota.
We have investigated the spatial and temporal distribution of Archaea in the Arctic Ocean to complement previous work on the Bacteria (5) and to obtain a more robust assessment of the diversity and richness of prokaryotes in the Arctic Ocean. We were also interested in comparing the results of our assessment of Archaea community composition with similar studies of Archaea in antarctic waters to identify differences, if they exist, in the compositions of these assemblages, with the ultimate goal of relating them to differences in the environmental conditions of the two polar oceans.
|
|
|---|
![]() View larger version (55K): [in a new window] |
FIG. 1. Locations of Arctic Ocean stations where the samples used in this study were collected. The locations of samples used to generate clone libraries are shown in boldface and larger.
|
|
View this table: [in a new window] |
TABLE 1. Environmental data for Arctic Ocean samplesa
|
DGGE.
The v3 region of the Archaea 16S rRNA gene was amplified from genomic DNA using the Archaea-specific (35) forward primer at position 344 (ARC344f; 5'-ACGGGGCGCAGCAGGCGCGA-3') and a universal reverse primer at position 517 (517r; 5'-ATTACCGCGGCTGCTGG-3'). All nucleotide positions given are in reference to the Escherichia coli 16S rRNA gene (9). A 40-bp GC clamp (5'-CGCCCGCCGCGCCCCGCGCCCGTCCCGCCGCCCCCGCCCC-3') (31) was added to the 5' end of the ARC344f primer, and fluorescein was attached to the 5' end of primer 517r. Each PCR mixture (100 µl) contained 1x PCR buffer (Promega), 2.5 mM MgCl2, 200 µM (each) deoxyribonucleoside triphosphate (dATP, dCTP, dGTP, and dTTP), 1.0 µM (each) primer, and 20 to 100 ng of template DNA. PCR amplification was performed under the following conditions: initial denaturation of the template DNA at 95°C for 5 min, a pause at 82°C to add Taq DNA polymerase (5 U; Promega), and then 35 cycles of denaturation, annealing, and extension using a touchdown program as described previously (16). The concentration of the resulting PCR product was estimated by the Hoechst dye assay (32), and then the mixed-template product was resolved by DGGE as described previously (4, 5). Bands of interest were excised from the gel, DNA was eluted from them into 50 µl of water by incubation at 55°C for 1 h, and then the eluted DNA was amplified as described above. The PCR products of the second amplification were resolved by DGGE, bands were excised, and DNA was eluted and amplified as described above. The PCR product was purified using Wizard PCR preps (Promega) and sequenced on an automated sequencer at the University of Georgia Molecular Genetics Instrument Facility with primer ARC344f to yield
120 bp of sequence information.
DGGE banding patterns were compared using Molecular Analyst software version 1.12 (Bio-Rad) essentially as described previously (16). The gel image was converted to a densitometry scan, standards run in the outside lanes of gels were used to normalize the gel lanes, the similarities of pairs of densitometry curves were quantified using Pearson's correlation coefficients, and cluster analysis was performed using the unweighted pair group method with averaging with fine alignment. Standard lanes from all gels were included in this analysis as controls; they clustered together at >90% similarity.
Clone libraries.
Because the sequences from DGGE bands were too short to provide adequate phylogenetic resolution, we sequenced a total of 198 cloned fragments of Archaea 16S rRNA genes amplified from nine samples covering three seasons and the range of depths sampled (Table 2). Archaea 16S rRNA genes were amplified using the Archaea-specific primers 21f (5'-TTCCGGTTGATCCYGCCGGA-3') and 958r (5'-YCCGGCGTTGAMTCCAATT-3'). PCR mixtures were prepared as described above. PCR was performed under the following conditions: initial denaturation at 94°C for 5 min, followed by 25 cycles consisting of 94°C for 30 s, 55°C for 30 s, and 72°C for 60 s, and a final elongation step at 72°C for 10 min. To minimize the potential effect of PCR biases in single reactions (34), three or four independent PCRs were pooled. The combined PCR product was purified with QiaQuick PCR purification columns in accordance with the manufacturer's instructions. The amplicons were cloned using a PGEM T-Easy Vector cloning kit (Promega) as recommended by the manufacturer. Clones from each library were selected randomly and sequenced from either end at the Molecular Genetics Instrument Facility using plasmid-specific primers (or from both ends for a subset of representative cloned amplicons), yielding
600 bp of sequence per reading, or on an ABI 310 in the Department of Marine Sciences, yielding
425 bp of sequence per reading.
|
View this table: [in a new window] |
TABLE 2. Number and percent of Crenarchaeota and Euryarchaeota found in clone libraries constructed from samples collected during the SCICEX 95, 96, and 97 cruises
|
Nucleotide sequence accession numbers.
The sequences of cloned fragments have been deposited in GenBank (see Fig. 5 for the accession numbers).
![]() ![]() View larger version (131K): [in a new window] |
FIG. 5. Neighbor-joining trees showing phylogenetic affiliations of representative partial 16S rRNA gene sequences retrieved from Arctic Ocean samples to closely related database sequences. Only one representative from each clone library of sequences that are 99% similar is shown; the total number of clones represented by a sequence is given in parentheses. Clones from this study are indicated in boldface. Bootstrap values higher than 50% are shown. The tree is unrooted, and the bar indicates a Jukes-Cantor distance of 0.1. (A) Marine group I Crenarchaeota (positions 536 to 941); clone SB95-72 is used as an outgroup; clusters are named as in reference 24. (B) Marine group II and IV Euryarchaeota (positions 512 to 937); clone SB95-57 is used as an outgroup.
|
|
|
|---|
![]() View larger version (97K): [in a new window] |
FIG. 2. DGGE fingerprints of Archaea amplicons. Clone libraries were constructed from the samples marked by arrows and from two other SCICEX 97 samples; however, these samples were not analyzed by DGGE. The numbered boxes indicate bands that were excised for sequencing.
|
![]() View larger version (157K): [in a new window] |
FIG. 3. Denaturing gradient gel comparing antarctic and arctic samples. The numbered boxes indicate bands that were excised for sequencing.
|
![]() View larger version (26K): [in a new window] |
FIG. 4. Dendrogram showing similarities of DGGE banding patterns for Arctic Ocean samples shown in Fig. 2. Std, standards.
|
We also found that clones containing inserts affiliated with group I Crenarchaeota and group II Euryarchaeota were equally abundant below the surface mixed layer in Arctic Ocean samples. Several reports have indicated a distinct partitioning of the marine water column between Crenarchaeota and Euryarchaeota, with Euryarchaeota dominant in the photic zone and Crenarchaeota dominant in deep water, for example, in the Santa Barbara Channel and Monterey Bay, California (13, 25, 33); the Mediterranean Sea (24); the North Atlantic Ocean and North Sea (24, 33); and west of the Antarctic Peninsula (10). Massana et al. (24) concluded that group I Crenarchaeota are ubiquitous and abundant and are often the dominant prokaryotes in marine waters, whereas group II Euryarchaeota are dominant only at the surface in temperate waters. In the Southern Ocean, group I Crenarchaeota phylotypes were dominant throughout the water column and group II Euryarchaeota were almost absent from deep-water samples (12, 24). Pernthaler et al. (33) demonstrated seasonal variation in the relative abundance of Euryarchaeota, with the largest populations occurring in surface waters during spring and summer. Karner et al. (19) conducted a quantitative analysis of the vertical distribution of Archaea in the Pacific Ocean using FISH. In contrast to previous reports of the dominance of Euryarchaeota in surface waters, they found the same low abundances of Euryarchaeota throughout the water column (surface to 4,750 m). In accordance with other studies, they also found that Crenarchaeota were dominant in deep water. In another study, only Euryarchaeota sequences were retrieved from 3,000 m in the antarctic Polar Front (21).
Crenarchaeota.
The majority of the group I Crenarchaeota sequences we obtained (104 of 111) form a closely related, monophyletic clade with sequence similarities of 98 to 100% (Fig. 5A). This clade is identified as cluster 1-
, following the designation given by Massana et al. (24). Sequences in this clade are closely related (99 to 100% similarity) to sequences cloned from various environments, such as the antarctic (7), the North Atlantic at a depth of 500 m (27), marine geothermal fluids sampled near the Cascadia Margin (8), a deep-sea hydrothermal vent (38), sediments of the Aegean Sea (GenBank description), and coastal marine environments (11). All of the Crenarchaeota clones we obtained from surface mixed-layer samples belonged to cluster 1-
, which is consistent with previous reports (24).
Cluster 1-
contained sequences from deep water that were closely related (99%) to cloned sequences obtained from an East Pacific Rise deep-sea hydrothermal-vent sample (GenBank description) and particles from the North Sea (40), the subtropical Atlantic Ocean (GenBank description), and the Santa Barbara Channel (25). Our results are in agreement with previous reports and suggest that members of cluster 1-
are adapted to live in the aphotic zone (24).
Despite the overall dominance (56%) of group I Crenarchaeota sequences in our clone libraries, none of the sequences we obtained from DGGE bands were related to this group. This may be explained by a mismatch at the 3' end of the putatively universal 517r primer. There are two other mismatches in the middle of the primer against the sequences of both group I Crenarchaeota and group II Euryarchaeota. Despite these mismatches, this primer set amplified a diverse Euryarchaeota population (41). It has been shown that a single mismatch at or near the terminal 3' base of a primer affects PCR more dramatically than single mismatches located internally or at the 5' end of a primer (37). Vetriani et al. (41) also used primer 517r for amplifying Archaea and identified only three Crenarchaeota bands in their DGGE profiles, even though they detected much greater Crenarchaeota diversity in clone libraries from the same samples. They also concluded that this discrepancy is due to the mismatch at the 3' end of primer 517r.
Euryarchaeota.
Marine group II Euryarchaeota sequences were more diverse than group I Crenarchaeota sequences (Fig. 5B), consistent with previous studies (24). Sequences in this group were assigned to nine different clusters, depending on sequence similarities and tree topology. Another cluster (Fig. 5B, cluster 10) contained sequences affiliated with marine group IV Euryarchaeota. Shorter sequences from DGGE bands were contained in sequences representing five of these clusters. The discussion below is ordered by the number of the DGGE band as identified in Fig. 2 and 3, with clusters that did not contain sequences from DGGE bands discussed last.
Sequences in cluster 3 (Fig. 5B) contained the sequence obtained from DGGE band 1 (100% similarity). Figures 2 and 3 show that band 1 was intense in fingerprints of samples collected below the halocline (235 m) on SCICEX 96 and 97 and in the antarctic samples GER 500 m and DB 150 m but was faint or absent in fingerprints of samples from shallower depths. In contrast, this band was intense in halocline samples from SCICEX 95. Sequences in cluster 3 were obtained from 131- and 235-m Arctic Ocean samples and were 100% identical over 520 bp to a sequence cloned from a sample collected at 1,000 m in the subtropical Atlantic near Bermuda (GenBank description). These sequences were also 98 to 99% similar to sequences cloned from deep-sea hydrothermal-vent fluids (38) and from flounder feces (40).
Cluster 6 was dominated by sequences from the halocline. Sequences in cluster 6 contained (100% similarity) the shorter sequence obtained from DGGE band 2. Band 2 was intense in fingerprints of Arctic Ocean halocline samples (Fig. 2) and was found in fingerprints of antarctic samples (GER 125 m and below; DB 0 and 150 m). Sequences in cluster 6 were distantly related to sequences in GenBank. The most closely related sequence (96% similarity) was cloned from a coastal marine environment (11). Cluster 5 contained only sequences obtained from mixed-layer samples. Sequences in cluster 5 contained (100% similarity) the shorter sequence obtained from DGGE band 3. Band 3 was intense in DGGE fingerprints from Arctic Ocean surface mixed-layer samples (Fig. 2). This band was faint or absent in fingerprints from antarctic samples (Fig. 3). Sequences in cluster 5 were 97% similar to sequences retrieved from 10 m in the southern Aegean Sea (28).
Cluster 1 contained six sequences from the halocline and one from the mixed layer. Sequences in cluster 1 contained (100% similarity) the shorter sequence obtained from DGGE band 4. Band 4 was intense in fingerprints from samples of the Arctic Ocean halocline (Fig. 2). In antarctic samples, band 4 was intense in DGGE fingerprints from deep-water samples (GER 500), but it was faint or absent in fingerprints from surface water samples (Fig. 3). Although we did not retrieve this sequence from clone libraries constructed from SCICEX 95 and 96 samples, DGGE analysis indicated that the ribotype represented by band 4 was present in all three seasons. Cluster 1 sequences are 97% similar to a sequence cloned from samples collected in hypersaline, anoxic basins in the Mediterranean Sea (GenBank description).
All (except one) of the sequences in cluster 9 were from the mixed-layer samples, and they formed two subclusters (9A and 9B) that were 98% similar to each other. Sequences in subclusters 9A and 9B contained (100% similarity) the shorter sequence obtained from DGGE bands 6B and 6A, respectively, suggesting that the small differences between them are not due to cloning or sequencing errors. DGGE analysis showed that bands 6A and 6B were intense in fingerprints from Arctic Ocean mixed-layer samples. None of the sequences cloned from SCICEX 95 and 97 samples fell into this cluster; however, bands 6A and 6B were present in DGGE fingerprints of samples from these cruises (Fig. 2, verified by sequencing DNA extracted from bands), indicating the year-round presence of these ribotypes. Bands 6A and 6B were also intense in DGGE fingerprints of antarctic surface water samples, but they were absent from fingerprints of deep-water samples (Fig. 3). Bands 6A and 6B had similar electrophoretic mobilities in DGGE and were not always distinct. Sequences in clusters 9A and 9B were 99% similar to a sequence cloned from continental shelf sediments collected off Antarctica (GenBank description) and from marine surface waters of Antarctica (14), respectively. Massana et al. (24) reported that 25 clones from surface water and only 1 clone from deep water fell into their cluster II-
(the same as our cluster 9a and b) and concluded that organisms represented by sequences in this cluster were adapted to surface waters.
We were not able to establish correspondence between the remaining cloned sequences and sequences from DGGE bands. Sequences in clusters 2 and 4 were retrieved from halocline and deep-water samples, respectively, and were 97 and 99% similar, respectively, to sequences retrieved from samples collected in the southern Aegean Sea at 10-m depth (28) and from the subtropical Atlantic (GenBank description). Sequences that corresponded to cluster 7 were retrieved from all depths but dominated libraries from Arctic Ocean mixed-layer samples. A previous study by Massana et al. (24) of the distribution of Archaea in different oceanic regions also found that sequences in this cluster (our cluster 7 corresponds to a subcluster of their group II-ß) were retrieved primarily from surface samples, suggesting that this organism is adapted to surface conditions. We were unable to identify a DGGE band corresponding to these sequences in our Arctic Ocean samples. However, they were 99% similar to a sequence retrieved from an antarctic sample by DeLong et al. (14). We also found a DGGE band sequence (band 7) in our antarctic samples that was 97% similar to cluster 7 sequences. Band 7 was intense in DGGE fingerprints of samples from all depths in the antarctic (Fig. 3). Sequences in cluster 8 were 98 to 99% similar to sequences obtained from deep water of the northeast Atlantic (GenBank description) and the Santa Barbara Channel (25).
One clone from the surface mixed layer was assigned to cluster 10. This sequence was 96% similar to the sequence of clone DH148-Y16, a newly identified lineage of Archaea, marine group IV Euryarchaeota, related to Haloarchaeon (22). Group IV Euryarchaeota sequences were detected in antarctic samples and in the Atlantic and Mediterranean using group-specific primers (22). Prior to our study, group IV Euryarchaeota sequences had been retrieved solely from samples collected at bathypelagic depths and had not been reported in surface waters.
The sequence from DGGE band 5 was not found in any of the clones we sequenced. This band was present in fingerprints of Arctic Ocean samples from both the mixed layer and the halocline but was faint or absent from 235-m samples. In antarctic samples, the band was intense in fingerprints from deep-water samples (GER 500 m) and faint or absent at shallower depths (GER 250 m and shallower). The sequence we obtained from it was 95% similar to sequences of group III Euryarchaeota environmental clones: AEGEAN_62, retrieved from the southern Aegean Sea at 10-m depth (28), and ME-450_P9, retrieved from the Mediterranean Sea at 450 m (24). Group III Euryarchaeota have not previously been reported in antarctic coastal water samples; however, a sequence that was 88% similar to our band 5 sequence was retrieved from a water sample collected at 3,000 m from the polar front just north of the antarctic peninsula (21). Previous studies retrieved marine group III Euryarchaeota sequences primarily from mesopelagic samples (17, 24). Massana et al. (24) screened 2,328 Archaea clones from different oceanic provinces and found only one group III Euryarchaeota clone (retrieved from 450 m in the Mediterranean Sea). In contrast, Moeseneder et al. (28) reported that 12% of the clones retrieved from the photic zone at 10-m depth in the eastern Mediterranean Sea belonged to marine group III Euryarchaeota.
Our study shows that although there appears to be significant overlap in the compositions of deep-water Archaea assemblages of the two polar oceans, there appear to be significant differences among assemblages in surface water samples. Previous studies have reported that antarctic planktonic Archaea assemblages were dominated by group I Crenarchaeota and that group II Euryarchaeota were a minor fraction of the Archaea assemblage (29, 30). In contrast, analysis of our Arctic Ocean clone libraries suggests that the abundances of the two groups are similar there. In addition, most of the group II Euryarchaeota sequences we retrieved were <97% similar to previously reported antarctic sequences. Only clusters 7 and 9 contained Arctic Ocean sequences that were similar (>99%) to previously reported antarctic sequences.
The biogeography of polar Archaea appears to differ from that of polar Bacteria in two significant ways. The first is that the richness of Archaea sequences encountered in this study appears to be lower than the richness of Arctic Ocean Bacteria (5) reported previously. The Bacteria study (5) used a similar combination of DGGE and cloning to examine the distribution of Bacteria in the Arctic Ocean, so the difference in apparent richness is not likely to be due to differences in methodology. However, differences in the robustness of primers for Bacteria versus Archaea, as discussed above, may contribute to the apparent discrepancy. Murray et al. (29) reported a similar discrepancy in the antarctic, and Massana et al. (24), who screened >2,000 clones from all over the world, concluded that the diversity of marine planktonic Archaea is low and is comparable to that detected in previous studies based on fewer samples. Alternatively, differences in richness, and presumably also diversity, may represent a real difference in the population genetics of these two groups. Secondly, in contrast to Arctic Ocean Bacteria sequences, which in general were related to polar or deep-sea clones or isolates (5), more of the Archaea we detected in this study seem to be affiliated with ribotypes found in extreme environments, warm seas, and temperate surface watershabitats that are very unlike cold polar oceans. This may simply be a function of the relative abundance of Bacteria versus Archaea sequences in GenBank and of the distribution of sampling efforts.
This work was supported by NSF awards OPP 95-00237, OPP 96-25131 (reissued as OPP 97-96261), and OPP 98-09971 to J.T.H. S.R. was supported by a CURO Internship from the UGA Honors Program, and B.R. was supported by REU supplements to NSF grant MCB 99-77886.
|
|
|---|
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»