Previous Article | Next Article ![]()
Applied and Environmental Microbiology, September 2002, p. 4486-4494, Vol. 68, No. 9
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.9.4486-4494.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Lehrstuhl für Physiologie und Biochemie der Pflanzen, Universität Konstanz, 78457 Konstanz, Germany,1 NIOO Centre for Estuarine and Coastal Ecology, 4400 AC Yerseke, The Netherlands2
Received 26 December 2001/ Accepted 4 June 2002
|
|
|---|
|
|
|---|
The APP of Lake Constance is dominated by phycoerythrin (PE)-rich Synechococcus spp., whereas phycocyanin (PC)-rich species and eukaryotic algae represent less than 5% of the pelagic community (10). Between 1988 and 1994, we isolated 26 PE- and PC-rich Synechococcus strains from the pelagic zone. Molecular fingerprinting by restriction fragment length polymorphism of psbA genes revealed 12 unique genotypes (29). 16S rRNA sequence-inferred phylogenetic analysis identified all Synechococcus spp. isolated from Lake Constance as members of three lineages of the picophytoplankton clade sensu Urbach et al. (36), which also comprises marine Synechococcus and Prochlorococcus spp. (A. Ernst, S. Becker, U. Wollenzien, and C. Postius, submitted for publication).
The 12 pelagic Synechococcus isolates from Lake Constance differed in pigmentation, surface structures, and physiological characteristics potentially influencing growth rates as well as trophic interactions (for a review, see reference 29). Differences in pigmentation and photosynthetic performance were also described for isolates of the phylogenetically related, marine Prochlorococcus spp. (22). Different genotypes of this population had been traced by fluorescence in situ hybridization targeting the 16S rRNA in combination with light microscopy or flow cytometry (42, 44). The depth-dependent distribution of genotypes confirmed the hypothesis deduced from physiological studies that they represent ecotypes adapted to different environmental conditions (8, 22, 36, 42). Based on the characteristics of a freshwater ecosystem, we had proposed that different ecotypes form subpopulations, which occur successively throughout the year (5). The succession of subpopulations could be driven by seasonal changes in the biotic and abiotic environment, as described for the succession of algae in spring, summer, and autumn populations (33). Alternatively, the subpopulations of the pelagic Synechococcus populations could exhibit a chaotic behavior in time as proposed by Huisman and Weissing (14) for mixed populations of algae.
To get insights into the internal dynamics of the freshwater APP, we wanted to trace ecophysiologically distinct, isolated strains in their natural environment. However, a low plating efficiency, the lack of morphological differences (7), and the close phylogenetic relation of these isolates (Ernst et al., submitted) precluded the in situ study of ecotype dynamics by traditional or modern methods (e.g., fluorescence in situ hybridization).
PE-rich Synechococcus spp. isolated from the APP of Lake Constance exhibit a sequence divergence of 0 to 8 nucleotides in 1,456 nucleotides of the 16S rRNA, with four clustered mutations that were not specific for the lineage (Ernst et al., submitted). However, as with many members of the picophytoplankton clade, they possess a long and highly variable internal transcribed spacer (ITS-1) separating 16S ribosomal DNA (rDNA) and 23S rDNA, which reflects the phylogenetic relations inferred from 16S rRNA sequences (17; Ernst et al., submitted). Therefore, we chose target sequences in the ITS-1 to trace the isolated Synechococcus strains in their natural habitat.
Unlike structural components of the ribosomal operon, the noncoding sections of the ribosomal ITS do not accumulate in cells. Hence, for the detection of different organisms, the ITS-1 has to be amplified by PCR. In this study, the PCR products were quantified directly in a Taq nuclease assay (TNA) or analyzed by denaturing gradient gel electrophoresis (DGGE). The latter method allows homologous PCR products to be distinguished on the basis of their melting behavior in a gradient of denaturants. As little as one base difference in the sequence can be sufficient to distinguish two amplicons of the same length (9). This technique has a high resolution and is now frequently used to demonstrate shifts in the genetic diversity of microbial communities (26), but quantitative aspects, which are frequently implicated in the interpretation of DGGE data, have been questioned. PCR is a nonlinear reaction in which small differences in the amplification efficiency and competitive effects between primer and amplicons can lead to a significant bias within and between reactions (3, 21, 24, 35, 38, 40). Quantification by PCR conducted over a fixed number of cycles requires analysis of serially diluted samples to which endogenous amplification standards have been added (32).
Alternatively, PCR-based quantification can be achieved by real-time PCR technology, in which the accumulation of the PCR product or the activity of Taq polymerase are continuously monitored by fluorescent labels. A particular form of real-time PCR is the TNA, also called a 5' nuclease assay or TaqMan PCR. In TNA, the activity of Taq polymerase is monitored by the hydrolysis of an oligonucleotide probe (TaqMan probe) labeled with two fluorescent dyes (18, 19). In the intact probes, the fluorescence of 6-carboxyfluorescein (FAM), acting as a reporter, is diminished by a quencher, 6-carboxy-tetramethylrhodamine (TAMRA). If this probe binds to a complementary target sequence during the annealing phase of PCR, the probe becomes hydrolyzed by the 5'
3' exonuclease activity of Taq polymerase (13). The release of the reporter during each round of amplification allows for rapid detection and quantification of target DNA without the need for post-PCR processing (18, 19). The reaction is calibrated by establishing a correlation between the concentration of target sequences (number of genomes) in the assay prior to amplification and the cycle at which fluorescence reaches a threshold value, the threshold cycle (11). Real-time PCR allows control of the amplification efficiency in every reaction and absolute determination of templates in an assay (3). As DGGE and TNA are not bound to particular genes or gene products, the application of these techniques can be tailored to the phylogenetic level at which microbial abundance and diversity are studied.
For the analysis of the Synechococcus-type cyanobacteria in the APP of Lake Constance, several PCR assays with different specificities were developed. For analysis of the genetic diversity in a lineage of PE-rich Synechococcus spp., a fragment of the ITS-1 was amplified by nested PCR and analyzed with DGGE. For quantitative PCR, two different primer pairs and TaqMan probes were used. One assay was designed to quantify all known (isolated) Synechococcus strains and all strains present in the environment that share target sequences for primers and probe in conserved sections in the ITS-1 with these isolates. A second TNA was designed for the specific detection of cells that share primer and probe sequences with a single isolated ecotype, Synechococcus sp. strain BO 8807 (3). Cultured PE-rich strain BO 8807 is distinguished from the seven closest related PE-rich isolates by a rod-shaped morphology with a longitudinal axis of 4.42 ± 2.62 µm (7), a highly glycosylated S-layer (6), and a lowered clearance rate in predation experiments with a chrysomonade and a nanoflagellate (1, 25). The feasibility of DGGE and TNA in tracing a particular ecotype in a natural population consisting of numerous closely related organisms and unknown diversity is discussed.
|
|
|---|
Sampling.
Every 2 or 3 weeks 10 liters of a mixed water sample (integration over 0 to 8 m depth) was collected at a fixed sampling site in the northwestern part (Überlinger See) of Lake Constance. Additionally, surface water (0-m depth) was collected from the pelagic and the littoral zone at the fixed sampling site near the maximum depth of Überlinger See and above a water column of 5 m, respectively. After a first filtration step through a 30-µm mesh, the water was filtered under low vacuum (0.3 x 105 Pa) through a 3-µm filter (cellulose), which was changed after filtration of 2 liters to avoid clogging of the filter. For the collection of picoplankton, 0.15-µm filters (cellulose acetate, 47-mm diameter) were used and stored at -20°C until DNA extraction. Filters were received from Schleicher & Schuell, Dassel, Germany.
Epifluorescence microscopy.
For cell counts of cyanobacterial picoplankton, water samples from Überlinger See or cultures of unialgal Synechococcus spp. were fixed with formalin (final concentration, 0.2% [vol/vol]) and concentrated under vacuum on black polycarbonate filters with a 0.2-µm pore size (Nuclepore, Tübingen, Germany). To achieve equal distribution of cells on filters, a second cellulose nitrate filter (0.45-µm pore size; Nuclepore) was placed beneath the filter for concentration of cells. Using a Labophot-2 microscope (Nikon) equipped with a 100/1.25 oil objective and an interference filter combination BA 590, the cell number of all single-cell coccoid cyanobacteria was determined, whereby colony-forming cells were excluded. No eukaryotic picoplankton was detected, since staining with 4',6-diamidino-2-phenylindole and excitation with filter BA 450 showed no presence of eukaryotic cells. Samples were analyzed in triplicate, and from each filter the cells of 10 different areas were counted.
Isolation of DNA.
DNA from cultivated strains was obtained by a phenol-chloroform extraction method described previously (3). The concentration and purity of genomic DNA was determined by measuring the absorption ratio A260/A280. For the estimation of genome copy numbers for pelagic Synechococcus spp., a genome size of 3 Mbp was assumed. Using an approximate molecular mass for a base pair of 650 Da, 1 ng of genomic DNA represented 3 x 105 copies of Synechococcus spp. genomes. For DNA extraction, filters with cells were cut in 16 equal pieces along lines imprinted by the filter support (Nalgene, Braunschweig, Germany) during the filtration process. Pieces were incubated with 400 µl of 5% (wt/vol) Chelex-100 (sodium form; 100 to 200 mesh; Bio-Rad, München, Germany) in a reaction tube for 30 min at 100°C and shaken several times while boiling (modified after reference 39). After vortexing at high speed for 10 s and centrifugation for 2 min, the supernatant was used as template in TNAs or in conventional PCR. From each filter the DNA of three different pieces was extracted and analyzed in duplicate.
DGGE.
For analysis of environmental samples with DGGE, PCR products were amplified in two successive PCRs, conducted as a nested PCR. In 25-µl volumes, 2 µl of supernatant of Chelex-100 extraction from filter pieces (see above) was mixed with 200 nM concentrations of the primers PITSANF and PITSEND (from MWG, Ebersberg, Germany, or Interactiva, Ulm, Germany), 2.5 mM Mg2+, 2.5 µl of 10x reaction buffer, and 0.625 U of Taq polymerase from Qiagen, Hilden, Germany. The PCR was conducted in a PTC-100 thermal cycler (MJ Research, Inc.) using a two-step cycling program as follows: after an initial denaturation at 95°C (3 min), the program comprised 30 cycles of a 2.5-min annealing-extension step at 60°C and a denaturation step (40 s) at 94°C. The reaction was terminated with a final polymerization step at 70°C for 5 min. For the second reaction, 2 µl of the first assay mixture was used as a template in 50-µl-volume assays for amplification. The reaction mixtures contained primers PITSGCANF and PITSGC (Table 1 and Fig. 1; from Interactiva or Genaxis, Spechbach, Germany) and other components as described above, and the PCR was performed after an initial incubation of 3 min at 95°C as follows: 35 cycles of 1.5 min at 65°C for annealing-polymerization and 30 s at 94°C for denaturation. The reactions were terminated with a final step at 70°C for 5 min. DGGE markers were generated by amplification of 10 ng of genomic DNA from isolated Synechococcus strains in a single reaction, using the primers PITSGCANF and PITSGC to produce GC-clamped fragments.
|
View this table: [in a new window] |
TABLE 1. Primers and probes used in this study
|
|
View larger version (10K): [in a new window] |
FIG. 1. Target positions of PCR primers and probes in the ribosomal operon of Synechococcus spp. used in this study. The target sequence and the arrows indicating the 5'-to-3' orientation of the oligonucleotides are not drawn to scale.
|
For confirmation of the identity of fragments, bands were excised from DGGE gels and transferred to 30 µl of elution buffer containing 0.5 M NH4-acetate and 1 mM EDTA, pH 8. The extraction of DNA followed the instructions of Ausubel et al. (2). In brief, the excised gel piece was shaken in the elution buffer overnight at 37°C, 800 rpm, and DNA was precipitated twice with cold ethanol. The extracted fragment was first dissolved in 20 µl of sterile water, diluted 20-fold, and used as template in PCR for reamplification (see above, second step of nested PCR). The migration behavior of these fragments was rechecked in 10 to 40% gradient gels before double-stranded sequencing was performed at GATC GmbH, Konstanz, Germany.
TNAs.
Primers and labeled probes for TNA (Table 1 and Fig. 1) were developed on the basis of ITS-1 sequences of Synechococcus strains isolated from the pelagic zone of Lake Constance as described previously (3). The TNA was conducted as reported previously, but the sample volume was reduced from 25 µl to 10 µl, of which 4 µl comprised supernatant from a DNA extraction with Chelex-100. The assays were pretreated for 2 min at 50°C and 10 min at 95°C before a two-step cycling program with 45 cycles of annealing and extension at 60°C for 60 s and 15-s denaturation at 94°C was carried out. Reactions were performed in an ABI PRISM 7700 sequence detection system (PE Biosystems, Foster City, Calif.). A normalized fluorescence of
RQ = 0.02 was used as the threshold to identify the reaction cycle CT for construction of standard curves (11). The constant amplification efficiency,
c, of TNAs achieved with different primer-probe combinations was calculated from the slope s of log-linear calibration curves using the equation
c = 10-1/s - 1 (15). The reaction efficiency at the threshold value (
RQ = 0.02) was calculated for every TNA using a nonlinear least-squares fit of the sigmoid product curve determined in the TNA, calculated with the equation T(i + 1) = Ti [1 + Km (Ti + Km)-1] (32), in which T0 and Km were varied. In TNAs, the template and amplicon concentration Ti is equivalent to the fluorescence
RQi of the reporter in the ith cycle. The reaction efficiency is calculated at the threshold value by using the equation
0.02 = Km (0.02+ Km)-1 (32; for details, see reference 3). For the least-squares fit, 5 to 13 successive cycles starting with the second cycle of consecutive positive
RQ values were used.
|
|
|---|
![]() View larger version (56K): [in a new window] |
FIG. 2. Separation of PE-rich Synechococcus spp. from Lake Constance with DGGE. Ten to 40% gradients were used; for details, see Materials and Methods. (A) Migration behavior of seven PE-rich Synechococcus isolates from Lake Constance, for establishment of a marker. (B) Band pattern determined for environmental samples 1 to 11 from 1999 (compare with Fig. 4A). (C) Band pattern with environmental samples 1 to 10 from 2000 (compare with Fig. 4B). M, marker, represents a mixture of fragments depicted in panel A. One microliter from PCR mixtures with isolated strains (panel A and each strain in the marker) and 5 µl from assays with environmental DNA were applied per lane.
|
![]() View larger version (20K): [in a new window] |
FIG. 3. Standard curves obtained by the threshold cycle method in real-time PCR. For TNAs, 10-ml-assay volumes contained approximately 101 to 106 copies of Synechococcus sp. strain BO 8807 genomes, 5 µl of 2x TaqMan Universal PCR master mix (5 mM Mg2+ final concentration), 300 nM concentrations of primers, and 200 nM concentrations of probe. Primers used with probe S8807A were P8807AP and P8807AM; primers used with probe S100A were P100PA and P3. For PCR conditions, see Materials and Methods. Fluorescence threshold RQ = 0.02; s = slope. Amplification efficiency was calculated as follows: c = 10-1/s - 1. Error bars represent the standard deviations of six experiments.
|
|
View this table: [in a new window] |
TABLE 2. Cell counts and genomes of autotrophic picoplankton detected by TNA probes S100A and S8807A
|
|
View this table: [in a new window] |
TABLE 3. Numbers of cells and genomes of APP passing filters with 30- and 3-µm pore size
|
![]() View larger version (22K): [in a new window] |
FIG. 4. Abundance of cyanobacterial picoplankton (Synechococcus spp.) and Synechococcus sp. strain BO 8807 in the pelagic zone of Lake Constance during the growth periods 1999 (A) and 2000 (B). Genomes per milliliter were determined by TNAs; the results were not corrected for prospective filtration losses. The percentage of BO 8807 is the relative genome numbers of Synechococcus sp. strain BO 8807 ( ) compared to the total genome number of Synechococcus spp. ( ), determined with probe S100A. Error bars represent standard deviations of three filter pieces analyzed in duplicate.
|
0.02 was determined at the fluorescence threshold value
RQ = 0.02 by applying the equations deduced by Schnell and Mendoza (32; see Materials and Methods). In the TNA conducted for calibration,
0.02 varied between 0.970 and 0.988 (data not shown), and similar values were achieved in the environmental samples (Table 4), indicating unbiased quantification of Synechococcus sp. BO 8807 genomes in all samples. |
View this table: [in a new window] |
TABLE 4. Amplification efficiency in TNA for quantification of Synechococcus sp. strain BO 8807 in environmental samples
|
|
|
|---|
Critical factors in environmental analysis are the extraction efficiency, purity, and stability of nucleic acids from environmental samples (27, 28). In this study we employed Chelex-100, a chelating agent with an ion-exchanging feature, for efficient extraction of DNA from filters. Serial dilution of this environmental DNA and application in 10-µl TNAs did not indicate any inhibitory effects of the 4 µl of the extract (data not shown). Furthermore, it was demonstrated that the reaction efficiency of TNAs observed at the threshold cycle is similar for DNA of the laboratory-grown strain used for calibration and for DNA from environmental samples (Table 4). This ensured that false quantification caused by PCR inhibitors or by competitive conditions in the PCR did not occur.
Tracing of isolated strains by TNA and DGGE.
In this study, real-time PCR with TaqMan probes and DGGE were used to trace phylogenetically closely related, PE-rich Synechococcus spp. that form subpopulations in the APP of Lake Constance. For the amplifications in nested PCR for DGGE and TNA for quantitative tracing on the ecotype level, target sequences of primers and probes were selected that are conserved in the ribosomal operon of 12 PE- and PC-rich strains isolated from Lake Constance and in eight additional Synechococcus spp. isolated from other fresh and brackish waters, comprising five different lineages in a 16S rRNA inferred phylogenetic tree (Ernst et al., submitted). The target sequences of the PCR used to produce a GC-clamped fragment for DGGE were selected to amplify all known PE-rich Synechococcus spp. of the pelagic habitat. Additionally, a strain-specific primer-probe combination was designed for quantitative detection by TNA of a single ecotype, Synechococcus sp. BO 8807. From the high amplification efficiency
0.02 at the cycle used for quantification (Table 4), which indicates the absence of homologous fragments that are amplified but not detected by the TaqMan probe, we can assume that this primer-probe combination was highly specific for tracing of Synechococcus sp. BO 8807 in the natural habitat. Hence, the detection of genome numbers of this strain over 4 orders of magnitude during the two periods of observation was possible (Fig. 4). Not only the absolute genome numbers but also the fractional contribution of strain BO 8807 to the TNA-quantified population of picocyanobacteria were highly variable. The discrepancies between the two quantitative determinations could not be explained without a complementary study of the population by DGGE, demonstrating frequent shifts in the genetic composition of the population throughout the observation period (Fig. 2). However, while TNA allowed tracing of 10 genomes in assays for the construction of standard curves (Fig. 3) or 16 BO 8807 genomes per ml in the habitat (Table 4), tracing of single ecotypes by DGGE was not always possible. The latter method is not only limited by the relatively low sensitivity of gel-based analysis methods of PCR products but also by a PCR bias, known from competitive PCR, which leads to suppression of a minor constituent in samples containing more than one template (3, 32). This explains why Synechococcus sp. strain BO 8807, which by quantitative analysis was shown to frequently contribute less than 1% to the TNA-quantified population (Fig. 4), was only sporadically detected by DGGE.
Another problem arose from the number of variable positions (22 in the seven isolated strains) in the 194-bp sequence of genuine ITS-1 in the GC-clamped fragments, which apparently can lead to compensatory effects in the melting behavior of different fragments. This problem precluded the visualization of shifts between genotypes that dominated the populations in both years, of which we have three cultivated strains, BO 8808, BO 8809, and BO 9402. This problem also obscured the detection of strain BO 8807 in the year 1999 (Fig. 2), but the absence of the comigrating fragment in the year 2000 demonstrated that even among codominant genotypes a year-to-year variability can be observed in the APP.
Cell and genome numbers.
Comparison of microscopic cell counts and determination of number of genomes by TNA resulted in 2.6 to 4.1 genomes per cell in the cultivated Synechococcus sp. strain BO 8807, depending on the age of the culture and the TNA used (Table 2). However, the genome number calculated per cell of the natural population was
1 (Tables 2 and 3). This low number could not be explained by a lowered DNA extraction efficiency from filters, because the addition of cultivated strain BO 8807 to a natural sample yielded the genome number expected from the consideration of the individual determinations (Table 2). We therefore have to consider that the APP is dominated by Synechococcus spp. that do not share characteristics of the cultivated strain BO 8807. The most obvious deviant characteristic of this strain are elongated cells containing several genomes (Table 2), which appear in stationary cultures, while this may not be the case with coccoid species. Another possibility is that genera of the APP are missing in our culture collection that were not amplified by the PCR primers used in the general Synechococcus TNA (see above). Thus, although an experimentally determined 1:1 relation of genomes and cells per milliliter was obtained (Tables 2 and 3), a general application of a conversion factor of 1 remains questionable.
Population dynamics in the APP.
In the years 1999 and 2000 we determined the number of genomes of Synechococcus spp. with known genetic profiles contributing to the cyanobacterial APP of Lake Constance. These known species in the APP showed two abundance patterns known from previous studies, in which cells were counted by epifluorescence microscopy (10). An intensive spring bloom with up to 2 x 105 genomes per ml was observed in May 1999 (Fig. 4A). It was terminated by a drop of 1 order of magnitude occurring at the start of the clear-water phase at the end of May (D. Straile, personal communication). The APP recovered in July and diminished again in October, when stratification of the water column was lost. In 2000, the known species in the APP were much less abundant in May (2 x 104 genomes per ml), but nevertheless genome number decreased by more than 1 order of magnitude during the clear-water phase in mid-June (D. Straile, personal communication) (Fig. 4B). Unusually though, the population recovered immediately after the clear-water phase to reach more than 105 genomes per ml until in mid-August the population started to decline (Fig. 4B). In both years, organisms with the ITS-1 sequence signature of Synechococcus sp. strain BO 8807 were detected in every sample examined, but in most cases they contributed less than 5% to the total APP. The development of this subpopulation seemed to be largely independent from that of the other members of the APP, as demonstrated by the wide range of relative abundances in both years. Interestingly though, the range of relative abundance, 0.02 and 5.7% in 1999 and 0.03 to 5.2% in 2000, and a high relative abundance during the clear-water phase were similar in both years. In contrast to other Synechococcus ecotypes we have in culture, Synechococcus sp. strain BO 8807 forms rods completely covered by a regularly structured glycosylated protein, forming an S-layer (6). We assumed that this surface structure was responsible for reduced predation in feeding experiments with a chrysomonad and a nanoflagellate from Lake Constance (1, 25). These experimental observations may relate to the high relative abundance of this strain during the clear-water phase. However, the data demonstrate that during summer, possibly due to succession in the heterotrophic nanoplankton, this protection diminished. In both years, a relative and absolute minimum of Synechococcus sp. strain BO 8807, about 20 genomes per ml, was observed at the end of the summer bloom. Production of and lysis by viruses were unlikely causes of this decline, because the ecotype had a low abundance throughout the summer population (up to several hundred genomes per milliliter), which is at or below the threshold concentration of hosts required for successful phage replication in pure cultures (16, 43).
Concluding remarks.
In this study we showed for the first time the dynamics of populations and subpopulations in the autotrophic picoplankton of a deep lake. The quantitative approach with TNA unveiled shortcomings in the sampling procedure and led to the conclusion that there must be a group of Synechococcus-type cyanobacteria not detected by our PCR-based assays. However, within the group of detectable organisms the wide dynamic range of TNA is an invaluable advantage in studying dynamics in microbial populations. On the other hand, diversity cannot be demonstrated by TNA and, thus, cannot replace DGGE. The study also showed that the use of ITS-1 as a target and the use of specific primers facilitated the recovery of signals from isolated strains in their natural habitat.
This work was supported by Deutsche Forschungsgemeinschaft through Sonderforschungsbereich 454 "Bodenseelitoral."
This is publication 2846 of NIOO, Centre for Estuarine and Coastal Ecology, Yerseke, The Netherlands. ![]()
|
|
|---|
3' exonuclease activity of Thermus aquaticus DNA polymerase. Proc. Natl. Acad. Sci. USA 88:7276-7280.
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»