Previous Article | Next Article ![]()
Applied and Environmental Microbiology, April 2004, p. 2028-2037, Vol. 70, No. 4
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.4.2028-2037.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Biology Department, Woods Hole Oceanographic Institution, Woods Hole, Massachusetts 02543,1 Department of Biological Sciences, University of Southern California, Los Angeles, California 90089-03712
Received 26 September 2003/ Accepted 1 January 2004
|
|
|---|
|
|
|---|
Molecular genetic methods have become fairly well established as alternatives to morphological identification of organisms in natural microbial assemblages during the last 15 years (2, 11, 24). To date, these approaches have been developed and applied primarily for assessing bacterial and archaeal diversity, because many prokaryotic species are currently unculturable and morphological criteria for their identification are virtually lacking. More recently, genetic and immunological methods have been applied in studies of protistan communities to identify and enumerate individual species of ecological interest and to examine the diversity of these assemblages (14, 17, 19, 29, 37).
Denaturing gradient gel electrophoresis (DGGE) is a PCR-based tool that has been applied extensively and effectively for analyzing the phylotype diversity of bacterial and archaeal assemblages in different environments (for examples, see references 3, 13, and 21 to 23). Application of this method to eukaryotic microbes has been accomplished in only a few studies (5, 28, 36). Our work represents some of the first to utilize the method to analyze protistan population relationships in different microhabitats (ice, seawater, and slush) in the Antarctic marine environment.
In our analysis of the genetic diversity of Ross Sea protists, we have targeted the small subunit ribosomal DNA (srDNA) with both general eukaryote-specific primers as well as primers selective for groups within the microbial eukaryotic realm (e.g., diatoms) to generate products for DGGE analysis. Banding patterns from different samples were compared, and similarities were determined for both sets of primers. The most intense bands were isolated from the gels, sequenced, and examined for their taxonomic affinities by BLAST analysis (1). This work has shown that DGGE can be an effective alternative means of assaying the genetic structure of protistan communities in the Antarctic ecosystem. This approach permitted the efficient comparison of samples from different habitats over large spatial and temporal scales and allowed the genetic identification of members of the assemblage.
|
|
|---|
![]() View larger version (26K): [in a new window] |
FIG. 1. Cruise track of cruise NBP 99-01. Triangles indicate the stations sampled (see the text for samples analyzed in this report). W, water; I, ice; SL, slush. Numbers are individual sample indicators. The map was generated by Antarctic Support Associates as part of the cruise data.
|
Primer design.
Primers for amplification of DGGE fragments were derived from primers that we have previously used for the amplification and sequencing of eukaryotic small subunit rRNA genes, 960F and 1200R (Table 1). A GC clamp (20) was added to the 5' end of 960F, and the primer sequence was modified to remove the three bases at the 3' end. This primer pair amplifies a region of the srDNA known as V7 and is just prior to the region amplified by van Hannen et al. (36).
|
View this table: [in a new window] |
TABLE 1. Primers used in this study
|
PCR amplification.
PCR amplification (27) was carried out to generate products for DGGE. One microliter of extracted nucleic acid was used as the template in 50-µl (with drop of oil overlay) touchdown amplification reaction mixtures with an initial annealing temperature of 65°C for two cycles (45 s at 95°C, 45 s at 65°C, and 45 s at 72°C), followed by a decrease in annealing temperature of 2°C every two cycles down to 55°C. Twenty-five cycles were then carried out at the 55°C annealing temperature (45 s at 95°C, 45 s at 55°C, and 45 s at 72°C). All amplification reactions were carried out using a Perkin-Elmer model 480 thermal cycler. Triplicate DGGE PCRs were carried out, visualized on 1% agarose gels stained with ethidium bromide, and then combined and ethanol precipitated for analysis. The precipitated products were resuspended in 5 µl of sterile distilled water and 5 µl of DGGE loading dye (40% Ficoll 400, 10 mM Tris [pH 7.8], 1 mM EDTA, 0.1% bromphenol blue), and 3 to 5 µl was loaded per lane on the gel.
PCR amplifications using group-selective primers were carried out in standard 50-µl reaction mixtures using 1 µl of extracted nucleic acid as template and 35 cycles of 45 s at 95°C, 45 s at 42°C, and 3 min at 72°C. Each group-selective primer was used in conjunction with a general eukaryotic small subunit ribosomal gene end primer (18). Primer A was paired with reverse group-selective primers, and primer B was paired with forward group-selective primers (Table 1). One microliter of this reaction mixture was used as the template for amplification with DGGE primers as described previously.
DGGE standard markers.
Standard markers to enable the comparison of different DGGE gels to each other were generated from protistan cells that were micropipetted from seawater or slush samples. These cells included an oligotrich ciliate, two different tintinnid ciliates, a copepod, the foraminifer Neogloboquadrina pachyderma, and solitary and colonial forms of Phaeocystis. Each of the cell nucleic acid extracts for these species yielded a single band by amplification and in DGGE analysis, and each of the bands had a different position or melting point in the gel. The DGGE amplification products were cloned, and plasmid preps were generated as stocks for future amplification of markers. To generate marker mixes, each plasmid was amplified separately, products were precipitated separately, and the presence of a single DGGE band was confirmed. All seven products were then combined in approximately equal proportions.
DGGE gels and analysis.
Our DGGE gels were run with a denaturing gradient of 35 to 85% urea at 60°C at 95 V overnight (16 h) using the CBS Scientific model DGGE-2000 gel apparatus. Band patterns were analyzed within and between gels using the BioImage whole-band analyzer software (Millipore Corporation). Bands were automatically detected by the software, but final confirmation of their presence (or absence) was accomplished manually. We added bands that we believed were overlooked, or we deleted ones that appeared to be miscalled by the program. The use of markers permitted the standardization of gels and allowed us to compare lanes within and between gels. Bands were identified as being the same using a 1% deviation in band position to account for slight differences in gel runs. Density differences between bands were not established, as the amount of PCR product loaded for each sample was not quantified. Dendrograms based upon the presence or absence of bands at each position were generated using the unweighted pair group method with averages (UPGMA) cluster analysis parameters in the BioImage software.
Sequence analysis of isolated DGGE bands.
Individual bands were cut from the DGGE gel using new razor blades, placed in 25 µl of sterile Milli-Q water, and allowed to elute overnight at 4°C. Isolated bands were identified by a gel identification number and a band number (e.g., in Fig. 2A below, band 1 is number 12500.1). The water containing the eluted DGGE band was recovered and placed in a new sterile microcentrifuge tube. Two or three microliters of the eluted band was reamplified with non-GC-clamped 960F and 1200R primers. These products were ethanol precipitated to remove excess primers and then resuspended in 10 µl of sterile Milli-Q water. Five microliters of product was sequenced directly using an ABI Prism Big Dye Terminator cycle sequencing ready reaction kit and the 960F or 1200R primer. Sequencing reactions were run on an ABI 377 apparatus (PE Applied Biosystems) and analyzed using Sequencher version 4.1 (Gene Codes Corporation). BLAST (1) searches were accomplished to determine relative taxonomic affinities of the DGGE band sequences. Bands that gave what appeared to be mixed sequences were cloned into the pGEM-T Easy vector system I (Promega). Four clones were picked for each band, and individual clones were sequenced as described previously.
![]() View larger version (98K): [in a new window] |
FIG. 2. Denaturing gradient gels for samples collected from the middle transect line of the cruise (Fig. 1). (A) Gel number 12500. Lane 1, W126; lane 2, W137; lane 3, I141; lane 4, W145; lane 5, SL148; lane 6, I 146; lane 7, W150; lane 8, SL154; lane 9, I155a; lane 10, I155b. (B) Gel number 12500. Lane 1, W159; lane 2, SL163; lane 3, I164; lane 4, W168; lane 5, SL172; lane 6, I 173; lane 7, W176; lane 8, SL180; lane 9, I181a; lane 10, I181b. Small circles with numbers indicate bands recovered and sequenced. a and b designate different melting methods for duplicate ice samples.
|
|
|
|---|
The greatest yield in number of bands from a single sample was nearly 30, but most samples yielded between 10 and 20 bands (Fig. 2). The gels in Fig. 2 are examples of several that were run and analyzed for this project. Some lanes had failed or faint reactions (e.g., Fig. 2A, lane 2), and these were reamplified and run on another gel for analysis. When samples were compared, usually more than half of the bands appearing on a gel were unique for each sample, with relatively few bands shared between samples. Ice samples shared the most bands when analyzed by DGGE, and these samples were also generally observed to have the lowest total number of bands. In the instances where we melted duplicate ice samples in 1 liter of sterile seawater or without seawater, the same major bands were usually present (Fig. 2A and B, lanes 9 and 10). Some minor bands differed between these replicates, indicating that the melting conditions may affect recovery of rare or delicate organisms. However, the dominant bands were still present, and thus the overall genetic pattern being detected often did not change appreciably. We did not assess whether the filtered sterile seawater had amplifiable DNA present, because we believed that any potentially contaminating DNA would have been in very low abundance compared to the target organisms present in colored ice or slush. This view was supported by the fact that we saw very little difference between duplicate ice sample DGGE patterns that were melted with or without the sterile seawater.
Dendrograms generated by comparisons of the banding patterns from all samples indicated that patterns obtained from the same microhabitat (i.e., water, ice, or slush) generally shared more bands than patterns from different microhabitat types (i.e., water versus ice, ice versus slush, or water versus slush), even when comparing samples collected from different locations (Fig. 3). Furthermore, assemblages collected from coastal samples of water and ice tended to group separately from those collected from the pack ice region of the Ross Sea when a single transect was examined (Fig. 3B). It should be noted that even when sample types were more similar to each other, they usually did not share a large number of bands (Fig. 3 [percent similarity scale]).
![]() View larger version (14K): [in a new window] |
FIG. 3. Dendrogram generated by a UPGMA cluster analysis comparison of DGGE patterns from direct amplification of ice, water, and slush samples collected in the Ross Sea, Antarctica. (A) All samples; (B) samples from the second transect only. Water samples 126, 137, and 148 and ice samples 141 and 146 were all coastal samples. a and b designate different melting methods for the same ice sample.
|
Sequences of recovered DGGE bands.
Partial nuclear srDNA sequences for microalgae (phototrophic protists), protozoa (heterotrophic protists), copepods, and metazoa were identified from the recovered and reamplified DGGE bands (Tables 2 and 3). It was feasible to obtain sequence information directly from some of the reamplified bands, but more often than not it was difficult to obtain unambiguous sequences. It seemed likely that more than one band was being recovered, either due to the presence of multiple bands that were not separated on the gels or because the method of cutting out bands did not resolve very closely spaced bands. Therefore, DNA from bands that did not give clear sequencing results were cloned. We rarely recovered more than one sequence type when analyzing cloned bands. This result could have been due to either the low number of clones that we sampled and sequenced (usually four) or to the fact that only the dominant sequence type was recovered. A small amount of contaminating sequence may have caused a problem for direct sequence analysis but yet not have contributed significantly in the cloning process.
|
View this table: [in a new window] |
TABLE 2. Primary BLAST matches of sequences obtained from DGGE bands from direct amplification products
|
|
View this table: [in a new window] |
TABLE 3. Sequences of DGGE bands from group-selective amplification products
|
Bands isolated from similar positions in the denaturing gradient (but different samples) generally yield the same sequence, but in several instances bands from different positions yielded similar sequences. For example, bands 12500.2, 12500.9, and 12500.21 (Fig. 2; Table 1) share the same position and the same diatom sequence affinity (Thalassiosira antarctica), while band 120600.7 shares the same sequence affinity but not the same gel position. The unidentified dinoflagellate bands 120600.34 and 120600.39 occurred at the same gel position, but 120600.25 did not. These similarities in taxonomic affinity, but differences in gel position, indicate that DGGE was able to detect genetic differences that may relate to species- or strain-level variation that would normally be impossible to identify using traditional methods. It is possible that in our evaluation of the DGGE sequences we have missed some of the sequence information that could account for the difference in gel position, because the length of sequence that we were able to recover from the bands varied (from 105 to 237 bases out of approximately 250 total bases) and even a single base difference can result in the distinct separation of two bands by DGGE (30).
Some of the bands recovered from the direct DGGE amplification of seawater samples yielded copepod sequences despite our prefiltration through 200-µm Nitex screening (5 out of 22 bands sequenced). Juvenile specimens, and/or appendages lost from adult copepods, may have gone through the screening, with the multicellular structures contributing significantly in the amplification of eukaryotic ribosomal genes in some of our samples. Copepod sequences were also recovered from ice samples, but these samples weren't prefiltered prior to extraction. It is of interest that the copepod sequences yielded the same type of BLAST search results that our protistan sequences didthey were similar but generally not identical to the sequences in the database.
Group-selective amplification.
The recovery of sequences from some protistan populations that were visually observed to be abundant in the samples (e.g., small flagellates) was less than anticipated (Table 2). This result may have been due to the inherent bias of PCR to amplify the most abundant sequences present in the sample, inefficient extraction (or damage) of nucleic acids from some protistan species, and our inability to recover all bands from the gel. These issues led to our development and application of group-selective amplification primers. Group-selective oligonucleotides were utilized in a primary amplification from three water DNA samples. DGGE primers were then used to reamplify the group-selective products, and these fragments were separated by DGGE.
Using group-selective primers for the primary amplification, it was possible to acquire bands not detected in the original amplification (Fig. 4). The band profiles shown for each of the group-selective amplifications contained many bands not represented in either the original direct amplification (compare lane 1 to the others) or in the other selective amplifications from the same sample. In Fig. 4A, 29 of the bands observed using the various selective amplifications were not found in the direct amplification. In Fig. 4B, 22 of the bands were not observed in the direct amplification. This result indicates the efficacy of the group-selective primers to recover previously undetected sequences and illustrates that direct amplification DGGE examines only a portion of what appears to be a highly diverse community.
![]() View larger version (71K): [in a new window] |
FIG. 4. Denaturing gradient gels for group-selective amplifications of W145 (gel ID number 120600) (A) and W150 (gel ID number 120600) (B). Lanes correspond to different amplification primers. Lane 1, direct DGGE amplification; lane 2, choanoflagellate-selective primer amplification; lane 3, bicosecid-selective primer amplification; lane 4, diatom-selective primer amplification; lane 5, A/B primer amplification; lane 6, prymnesiophyte-selective primer amplification; lane 7, chrysophyte/synurophyte-selective primer amplification. M, marker. Small circles with numbers indicate bands recovered and sequenced.
|
![]() View larger version (30K): [in a new window] |
FIG. 5. Dendrogram generated by UPGMA cluster analysis comparison of DGGE patterns from group-selective amplifications of three water samples. Labels indicate the sample number (W145, W150, and W176) followed by the group-selective amplification primer used.
|
General eukaryotic primers A and B (18) were also used to amplify full-length nuclear srDNA, the products of which were reamplified with DGGE primers. After examining the gels in Fig. 4, it was apparent that some of the change in band patterns was likely due to the change in template abundances as a result of the reamplification, although several of the bands that were detected in the original gels were also present. We view this general approach as being useful for expanding the amount of genetic information recovered from a sample and providing a more complete taxonomic profile of the community.
|
|
|---|
Phylotypes that we recovered from within ice cores were typically quite different than those recovered from seawater. This finding is contrary to studies that have suggested that species found commonly in the plankton are also common in sea ice (8, 9, 25). The uniqueness of the banding patterns from ice and seawater samples does not necessarily mean that the same species were not present in both microhabitats, but rather it implies that the dominant species in the two types of microhabitats differed markedly (or that morphologically similar but genetically distinct strains were present).
Despite the low degree of commonality among all samples in this study, samples from the same microhabitat but from different sampling sites showed more similarity than samples from the same sampling site but different microhabitats. This finding implies that environmental conditions within each of the microhabitat types enrich the abundances of habitat-specific types of protists from within the complex mixture of species present in the Ross Sea ecosystem. Some degree of similarity among the microflora and microfauna is not surprising given the results of previous microscopic studies, but the horizontal extent of this trend is striking given the temporal and spatial sampling scales involved in this study. Each station on a transect is 1° in latitude away from its closest neighbor, which is 60 nautical miles. Samples along the same transect were collected within 7 days of each other, whereas almost a month elapsed between the start of the first transect and the end of the third. It is highly noteworthy that we have been able to discern a general pattern of genetic similarity in different microhabitats despite the significantly large scales involved.
Analyses of samples from a single transect enabled discrimination between the samples that we considered coastal (open water within the polynya on the southern side of the pack ice and next to the continent, including land-fast ice) versus those that were collected from within the pack ice (Fig. 3B). Coastal water samples and ice samples formed two separate groups within the pack ice water, slush, and ice samples. These groupings indicate subtle differences in the protistan assemblages of the coastal and pack ice environments. Land-fast ice (I141 and I146) may contain multiyear growth of microbial communities and thus might be expected to contain different assemblages than ice in the pack ice region. Similarly, succession in coastal water within the polynya could explain the divergence of these planktonic communities from the water within the pack ice. The latter is more ephemeral temporally and may differ significantly from water in the polynya in physical features and nutrient chemistry. Trends related to coastal versus pack ice noted along individual transects were not evident when all samples were considered together (Fig. 3A). This inability to distinguish clearly between coastal and pack ice samples when all three transects were combined was likely due to the smaller number of bands actually shared between the samples as the communities changed over the month-long sampling scale of the project.
In DGGE analyses, the total number of bands in a lane is often taken to be a measurement of the diversity of that sample, but it is perhaps more appropriate to view the number of bands as representing the composition of the dominant populations in the sample (7). This view is supported by both the gel results (Fig. 2 and 4) and the sequence results (Tables 2 and 3). On our DGGE gels, we were able to observe up to 27 bands in a single sample by the direct amplification of that sample (Fig. 2B, lane 5). When we used the group-selective amplification strategy, we increased the total number of bands observed per sample by almost an equal amount (20 to 30). This finding indicates a high degree of diversity in the protistan assemblages from the Ross Sea, implying that our direct amplification results potentially represent only the more abundant organisms in the samples. This finding also implies that further work and new priming strategies will be necessary to fully characterize the diversity of these assemblages using this molecular approach.
Our sequencing results for DGGE bands further indicated that we had characterized only some of the protistan diversity in our samples. We were successful in recovering sequences from some taxa that have been described as conspicuous components of Antarctic protistan assemblages, such as Phaeocystis and dinoflagellates (Tables 2 and 3), but we had only limited success for others, particularly diatoms. Numerous diatom species have been reported in samples from the Antarctic environment, yet our results here indicated only three. We speculate that nucleic acid extraction efficiencies were low for some protistan species in our samples and, although the combination of bead beating and hot detergent was fairly rigorous, we may have had limited success at breaking open some cells.
A practical consideration explaining the relatively low number of bands recovered and sequenced in this study was that we were not able to exhaustively sample all of our DGGE bands. Some bands were simply too faint to be isolated reliably from the gels, and so we were constrained to analyzing the brightest and most intense bands. The brightness of a band is indicative of its concentration in the total PCR product, and that can reflect the abundance of the original template, although there are other factors that can affect this as well (26, 34). By examining the brighter bands, we were presumably observing a limited population of the most dominant or abundant organisms present in the original sample. Better methods for band isolation and removal and cloning of reamplified products will likely help to resolve this issue, although some bands will remain difficult to recover.
Morphology-based studies of Antarctic protistan diversity have described the presence of diverse and dynamic assemblages of protists in the water surrounding the continent (as examples, see references 9, 25, and 31 to 33). Our molecular approach to community composition generally is in agreement with those published reports in the sense that we have shown that the ice, slush, and seawater assemblages are genetically quite diverse and harbor protistan assemblages that can be quite heterogeneous on both the macroscale (Ross Sea) and microscale (ice and slush). In particular, our genetic results compare very well with the morphological analyses of Garrison and Buck (8). As in that study, we were able to document the distinct assemblages of organisms formed by ice and water samples, along with spatial heterogeneity among the water, ice, and slush samples. Interestingly, our DGGE analysis indicated that there was little genetic overlap between different sample types. Even samples that were considered similar still often shared fewer than 50% of their bands. Our finding of high diversity (low similarity) among samples in the same group is in contrast to the 54% similarity of the organisms from different water samples in Garrison and Buck's morphological analysis (8). This difference may be due to the fact that we examined the total microbial eukaryotic community rather than restricting our analysis to the phytoplankton. Nevertheless, the qualitative similarities between morphology-based and culture-independent approaches for assessing protistan diversity are encouraging. Given the relatively youthful state of development of DGGE and the improvements in resolution and speed of analysis expected in the future, this molecular approach should prove highly useful for ecological studies of these species. Continued development should provide insightful, culture-independent techniques for assessing the diversity of protistan assemblages in Antarctic microhabitats.
We thank everyone involved in making the NBP99-01 cruise a success, especially the captain and crew of the RVIB Nathaniel B Palmer and Martin Jeffries (chief scientist). Special thanks are given to Alison Murray for her assistance with the DGGE methods. We also thank Peter Countway, Dawn Moran, David Beaudoin, and three anonymous reviewers for their insightful and productive comments regarding the manuscript.
WHOI contribution number 11028. ![]()
|
|
|---|
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»