Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2003, p. 6723-6730, Vol. 69, No. 11
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.11.6723-6730.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Institute for Limnology, Austrian Academy of Sciences, A-5310 Mondsee, Austria,1 Federal Environmental Agency, D-14195 Berlin, Germany2
Received 23 June 2003/ Accepted 29 August 2003
|
|
|---|
|
|
|---|
Recently, the TaqMan PCR, or the Taq nuclease assay (TNA), was introduced to quantify specific genotypes of picocyanobacteria (1) or microcystin-producing cyanobacteria in the field (10). This technique utilizes a sequence specific dually labeled fluorescent probe (TaqMan probe) and primers to quantify the level of DNA template initially present in a sample. The rate of exponential accumulation of the amplicon is monitored by the hydrolysis of the TaqMan probe, in which it generates a fluorescent signal during the amplification process. The threshold cycle (Ct) is the PCR cycle number at which the fluorescence passes a set threshold level and can be used to determine the starting DNA amount in the sample based on a standard curve (based on samples with a known concentration).
The usefulness of the TNA in microbial ecology could be constrained by two factors. Firstly, the DNA content of a specific organism is dependent upon the growth rate which has been shown to be maximal under the highest-growth-rate conditions (16). Consequently, physiological conditions may influence the result on the quantification of genotypes and calibration is usually done based on DNA concentration or genome copy number. However, for population ecology the quantification in terms of cell numbers would be most useful. Secondly, due to the linear-log calibration curves, the technique is very sensitive to variations in the slope induced by minor variations of the Cts. It is not known whether spatial or temporal variations in genotype abundance in nature can be reliably monitored in spite of the noise induced by the semilogarithmic calibration algorithm alone.
It was the aim of the study to develop a quantitative PCR approach useful for quantifying the total population of Microcystis sp. (as defined in reference 20) as well as the subpopulation comprising all microcystin genotypes in terms of cell numbers and for monitoring the seasonal development of microcystin genotypes in Lake Wannsee (Berlin, Germany). We developed two independent TNAs, one to quantify the total population of Microcystis sp. using the intergenic spacer region within the phycocyanin (PC) operon and the other assay to quantify microcystin genotypes using a region of mcyB part of the microcystin synthetase gene cluster (28). To validate the results of the TNA, (i) strains were quantified in batch culture during the growth cycle and compared in amplification efficiency in the laboratory, and (ii) cell numbers estimated by the TNA in the field were compared to cell numbers estimated by counting under an inverted microscope.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Microcystis species and strains used to test primer and TaqMan probe sensitivity and specificity (HUB524), to monitor DNA content over the growth cycle (HUB524), and to compare amplification efficiencies of individual strains
|
DNA extraction and TNA.
Extraction from frozen filters was performed as described previously (15). The TNA was used to quantify two specific gene regions, the intergenic spacer region within the PC operon and the mcyB region, which encodes one step in microcystin biosynthesis (28). A variable gene region of the PC operon was selected based on an alignment (ClustalW 1.8) of PC genes from several genera of cyanobacteria from GenBank, including Microcystis (AJ003168, AF385388; R. Kurmayer, unpublished data), Planktothrix (AJ132279 and AJ131820), Oscillatoria (AJ401186 and AJ401185), Nostoc (X05239), Synechocystis (AJ003180), Chroococcus (AJ003189), Synechococcus (AF223465, AF223464, and AF223463), and Lyngbya (AJ401187). The variability of the intergenic spacer region of the PC gene was sufficiently high to achieve specificity for amplification of the PC operon of Microcystis even in the presence of other cyanobacteria (e.g., Planktothrix; see below). The mcyB gene region was selected from Microcystis strain HUB524 (Z28338) and was located between core motifs A2 and A3 (17). From those gene regions, optimal primers and TaqMan probes were designed with Primer Express 2.0 software (Applied Biosystems, Vienna, Austria). The primers were 188F (5'-GCTACTTCGACCGCGCC-3') and 254R (5'-TCCTACGGTTTAATTGAGACTAGCC-3') for the PC operon and 30F 5'-CCTACCGAGCGCTTGGG-3' and 108R (5'-GAAAATCCCCTAAAGATTCCTGAGT-3') for mcyB, with amplicon sizes of 66 and 78 bp, respectively. The TaqMan probes were 5'-CCGCTGCTGTCGCCTAGTCCCTG-3' for the PC operon and 5'-CACCAAAGAAACACCCGAATCTGAGAGG-3' for mcyB. These probes each had a fluorescent reporter dye (6-carboxyfluorescein) covalently attached to the 5' end (5'-FAM) and a 3'-TAMRA fluorescent quencher dye (6-carboxytetramethylrhodamine). PCR was initiated with two holds, one for 2 min at 50°C (AmpErase UNG protection against carryover contamination) and one for 10 min at 95°C. Subsequently, a 95°C denaturation step for 15 s was followed by a 60°C annealing and extension step for 1 min, for 45 cycles. PCRs were performed with a GeneAmp 5700 sequence detection system (ABI, Vienna, Austria) using SDS 1.3 software. TNAs were performed with a volume of 25 µl, containing 12.5 µl of 2x TaqMan universal PCR master mix (ABI), a 300 nM (300 fmol µl-1) concentration of each primer, a 100 nM concentration of the TaqMan probe, and 5 µl of template containing various amounts of genomic DNA and filled up to 25 µl with sterile Millipore water. For mcyB, a 900 nM concentration of each primer and a 250 nM concentration of the TaqMan probe were used. Each measurement was performed in triplicate.
Specificity and sensitivity of the TNA.
For both genes (the PC operon and mcyB), a standard curve based on predetermined cell concentrations was established by relating the known DNA concentrations (in cell equivalents) to the Ct of the diluted samples. The threshold value for the fluorescence of all samples was set manually at 0.1 in accordance with the instruction manual of the GeneAmp 5700 sequence detection system. Five milliliters of a suspension of 2.08 x 107 cells of HUB524 ml-1 (determined by electronic particle counting) were filtered on a GF/C filter, and from the DNA extract, six dilutions ranging from 1:102 to 1:2 x 106 of template DNA (equivalent to 104,000 cells to 5.2 cells) were prepared. To test the specificity of both TNAs in the presence of a highly diversified DNA background from a natural phytoplankton community, water from Lake Wannsee was filtered on 4 May 2000 through a sieve (25-µm mesh size), selectively removing Microcystis but permitting other cyanobacteria to pass through the filter. Three hundred milliliters of the filtrate was filtered on a GF/C filter, and the DNA was extracted and added to the dilutions of DNA obtained from Microcystis strain HUB524. The counting of cyanobacteria under an inverted microscope (3.2-ml sedimentation chamber) revealed the dominance of filamentous cyanobacteria (Aphanizomenon spp., Limnothrix spp., Limnothrix redekei, and Planktothrix agardhii) but no Microcystis (see reference 15 for details on counting and biovolume determination). The dilutions 1:100 (DNA equivalent to a cyanobacterial biovolume of 6.6 x 10-4 mm3 in the template) and 1:1,000 (6.6 x 10-5 mm3) were used, and the maximum and minimum ratios of biovolume of cells of HUB524 to the background biovolume were 66 (at 100,000 cells, 1:1,000) and 6.6 x 10-4 (at 10 cells, 1:100), respectively (calculated with an average Microcystis cell biovolume of 42.2 µm3). To compare background effects, the number of cells estimated in the presence of a natural background was divided by the number estimated in the absence of a natural background, and the ratio of cells with background to cells without background was calculated.
Sample collection from Lake Wannsee.
Sampling was performed at Lake Wannsee (Berlin, Germany) from June 1999 to October 2000. The lake is shallow (mean depth, 5.5 m; maximum depth, 8.5 m; area, 2.82 km2), polymictic, hypertrophic, and regularly dominated by Microcystis sp. during the summer. The phytoplankton was sampled weekly (May to October) and fortnightly (November to April). Water samples were obtained at the deepest part of the lake and integrated into a 30-liter bucket by collecting 2 liters every meter from the surface to the sediment. Cells were harvested and extracted as described above. The cell number and the biovolume were determined by counting via the inverted-microscope technique. Details on counting procedures and ambient abiotic factors are provided in reference 15.
Statistical analyses.
The accuracy and the statistical significance of the measurements obtained by TNA were tested by regressing the TNA variables (Ct or cell numbers) against the respective cell numbers determined via the electronic particle counter or the microscope (as an independent variable). For linear regression analysis, the data were log transformed and tested for normal distribution (Kolmogorov-Smirnov test, P < 0.01) and for constant variance of the dependent variable regardless of the value of the independent variable (Spearman rank correlation, P < 0.01). The residuals were tested for their independence of each other (Durbin-Watson statistic; if the residuals are not correlated, the value will be 2 ± 0.5). All data passed those tests; however, the errors in field data were serially correlated (Durbin-Watson statistic = 1.1 [PC operon] and 1.07 [mcyB]). Consequently, the interval of field observations was increased and every third sampling date was omitted from statistical analysis to achieve the assumption that the error term is a random normal variable. The linear curves were fitted by using the least-square approximation and the associated statistical tests of Sigma Plot 2000 (version 6.10). For the field analysis, the two linear regressions between cell numbers as determined by microscopy and by TNA were compared in slope and intercept using a general factorial model of analysis of variance (ANOVA). The data were modeled as follows: yijz = µ +
j + ßz + (
ß)jz +
ijz, where yijz is the cell number measured by TNA, µ is the overall mean level,
j is the effect of the cell number in the microscope, ßz is the effect of the grouping factor as a covariate (PC operon and mcyB), (
ß)jz is the interaction between both factors, and
ijz is the unexplained part of the variance (27). An SPSS statistical package (release 6.0) was used for ANOVA.
|
|
|---|
![]() View larger version (24K): [in a new window] |
FIG. 1. (A) Standard curves (black symbols) for both the PC operon (circles) and mcyB (squares) based on predetermined cell concentrations of Microcystis HUB524 by relating the known DNA concentrations (in cell equivalents) to the Ct of the diluted samples. In addition, both TNAs were tested in the presence of a 1:100 dilution (6.6 x 10-4 mm3, white symbols) and a 1:1,000 dilution (6.6 x 10-5 mm3, dotted symbols) of a natural background (containing other cyanobacteria but no Microcystis) from Lake Wannsee. All data are means of three parallels ± 1 standard error (error bars that are not visible are hidden behind the symbols; however, error bars have been omitted for outliers shown at 105 cells). (B) Comparison between cell numbers estimated for both genes by TNA from the standard curve in the absence (white columns) or in the presence of 6.6 x 10-4 mm3 and 6.6 x 10-5 mm3 of background. (C) Number of cells estimated in the presence of a natural background divided by the cell number estimated in the absence of a natural background for both the PC operon (circles) and mcyB (squares) and both background dilutions (6.6 x 10-4 mm3 [black symbols] and 6.6 x 10-5 mm3 [white symbols]).
|
![]() View larger version (19K): [in a new window] |
FIG. 2. (A) Cell numbers of Microcystis HUB524 grown under high-light conditions in batch culture and quantified via electronic particle counting and TNA over 4 months (mean ± 1 standard error). Data are from two independent batch culture experiments. For reasons of clarity, only the cell numbers quantified with the PC gene are shown; for details on statistical correlation between cell numbers and both the PC gene and mcyB, see the text. The region within the dotted lines is considered the transition phase from logarithmic growth to the stationary phase. (B) Ratio of cell numbers determined by TNA to cell numbers determined by the particle counter during the same batch culture experiment (mean ± 1 standard error). The transition from the exponential growth phase to the stationary phase differed significantly (P < 0.001, Mann-Whitney rank sum test) from the earlier and later growth phases. (C) Regression of cell numbers (mean ± 1 standard error) determined by the TNA (PC operon) versus cell numbers determined by the electronic particle counter during the same batch culture experiment. Error bars that are not visible are hidden behind the symbols.
|
![]() View larger version (19K): [in a new window] |
FIG. 3. (A) Cell numbers of 10 Microcystis strains grown under high-light conditions in batch culture and quantified by the electronic particle counter and by the TNA (mean ± 1 standard error). Data are from two independent parallels. (B) Ratios of cell numbers determined by TNA to cell numbers determined by the particle counter observed among strains and during the stationary phase of strain HUB524. The boxes show the median (line) and the 25th and the 75th percentiles, and the whiskers indicate the 5th and the 95th percentiles.
|
Throughout the study period, Microcystis cells were detected during counting in the microscope. On 29 of 55 sampling dates, Microcystis made up more than two-thirds of the cyanobacterial biovolume. On both sampling dates in April 2000, the Cts were out of the range of the standard curve of the TNA for the PC operon, and on 22 August 2000, no result was obtained. Those values were omitted from further analysis. Over the entire study period a significant relationship between cell numbers estimated via the inverted microscope technique and TNA of the PC gene was found (Fig. 4A and D). For the PC operon the regression equation was y = -1.11 + 1.22x (R2 = 0.88, n = 34, P < 0.0001), where y is the log10 cell number as determined by the TNA and x is the log10 cell number counted in the microscope. In general, mcyB-containing cells were always found when the PC operon genotypes were detected (Fig. 4B). On six sampling dates (22 February, 12 and 25 April, 4 and 9 May, and 22 August 2000), the concentration of mcyB was lower than the detection limit (Ct > 45). The mean mcyB proportion over the study period was 11.4% ± 2.6% (95% CL, n = 46), with a minimum and maximum of 0.9 and 38.3%, respectively (one outlier on the 8 June 1999, 114%). The proportion of mcyB genotypes from June 1999 to September 1999 (11.9% ± 3.6% [95% CL]) did not differ significantly from that in the corresponding period in the summer of 2000 (12.4% ± 5.5%; Mann-Whitney rank sum test, n = 28, P = 0.88). It is concluded that microcystin genotypes constituted only the smaller part of the Microcystis population in Lake Wannsee. The cell number of mcyB genotypes was found to be one-to-one related to the cell number of the PC operon genotypes (Fig. 4C). The regression equation was y = -1.20 + 1.03x (R2 = 0.95, n = 50, P < 0.0001), where y is the log10 cell number as determined by the TNA (mcyB) and x is the log10 cell number determined by TNA (the PC operon). The linear regression between cell numbers determined by microscopy and by TNA for mcyB was y = -2.17 + 1.22x (R2 = 0.80, n = 34, P < 0.0001), where y is the log10 cell number determined by TNA and x is the log10 cell number counted in the microscope. When the regressions for both genes (Fig. 4D) in intercept and slope were compared via the ANOVA model, the influence of the microscopically determined cell numbers (
j) on the TNA estimation of cell numbers and the influence of the grouping variable (ßz) were found to be significant (P < 0.0001 [
j] and P < 0.1 [ßz], n = 69). In contrast, no significant interaction effect between both factors [(
ß)jz] was found (P = 0.98), indicating that both regressions differed significantly in intercept (ßz) but not in slope [(
ß)jz]. It is concluded that the mean proportion of microcystin genotypes was stable during the period of seasonal population growth from the lowest cell numbers in winter to the highest cell numbers in summer and did not depend on seasonal influences during the study period.
![]() View larger version (25K): [in a new window] |
FIG. 4. (A) Numbers of Microcystis cells in Lake Wannsee from June 1999 to October 2000, determined by counting under an inverted microscope or by TNA with the PC gene (mean ± 1 standard error). (B) Proportion of mcyB genotypes during the same period (one outlier on 8 June 1999, 114%, was omitted). The solid line indicates the mean water temperature integrated every meter over the total water column. (C) Dependence of cell numbers determined by TNA of the mcyB gene on cell numbers determined by TNA for PC. Error bars show the mean ± 1 standard error. (D) Comparison between cell number determined in the microscope and determined via TNA for PC (black) and mcyB (white) during the same study period. For details on statistical regression parameters, see the text.
|
|
|
|---|
The high correlation coefficients between cell numbers in the microscope and those determined by TNA in the field strongly suggest the robustness of both TNAs under natural conditions. Because the phytoplankton of Lake Wannsee typically consists of several cyanobacterial genera (Microcystis, Planktothrix, Aphanizomenon, Anabaena, Limnothrix, and Pseudanabaena) (9), it is expected that both TNAs will work with sufficient accuracy in any body of water. In addition, the population of Microcystis typically varies in abundance by 3 to 4 orders of magnitude over the year, and cell quantification on a logarithmic scale has sufficient accuracy to monitor the population cycle over time for most research and monitoring objectives. Compared with the time involved with the counting of cells of Microcystis under the inverted microscope, TNA can be considered less time-consuming and labor-intensive. The relatively high acquisition and running costs in addition to the costs for acquisition and maintenance of standard laboratory equipment (usually comprising an inverted microscope) may constitute a drawback in applying TNA on a broader scale.
There was a large significant increase in TNA fluorescence signal during the transition from the logarithmic phase to the stationary phase. The reason might be an uncoupling between DNA replication and subsequent cell division when resources become limiting. This result is in contrast to the data obtained by Mann and Carr (16) for Anacystis nidulans demonstrating that faster-growing cells have a higher DNA content per cell. However, Anacystis has a much shorter mean generation time (3 h) than Microcystis (in this study, the maximum generation time was 41 h), and for Microcystis the time needed for DNA replication is probably never a limiting factor during cell division as observed for Anacystis.
Seasonal stability of microcystin genotype proportion.
In this study the averaged proportion of microcystin genotypes was considered constant from 1999 to 2000, as shown by the parallel regression lines between cell numbers in the microscope and cells containing or lacking the microcystin gene as differentiated by the TNA. In addition, the correlation between PC genotypes (as the independent variable) and microcystin genotypes (the dependent variable) was close to 1, and 95% of the variance of mcyB was attributable to the influence of population growth (the increase in PC concentration). The results are in agreement with the conclusion obtained from the genetic analysis of different colony size classes (15) that shifts of microcystin genotype proportions between and within size classes are stable over the growing season.
The stability in microcystin genotype proportion between years may depend on inoculation and reinvasion events after strong winters when the total population was removed from the water column. Microcystis sp. are known to rely on reinvasion into the water column from the sediments (23). It has been reported elsewhere (15) that the Microcystis population of Lake Wannsee consists mainly of the morphospecies M. aeruginosa/Microcystis flos-aquae and M. ichthyoblabe. Those morphospecies also differ significantly in the percentage of microcystin genotypes; i.e., 73% of colonies assigned to M. aeruginosa contain mcyB, while only 16% of colonies assigned to M. ichthyoblabe contain it. M. aeruginosa has frequently been reported to form large and firm colonies, while the colonies of M. ichthyoblabe are typically small and fragile (31). It remains unknown whether there are differences in survival ability in the sediment between morphospecies (e.g., see references 2 and 3), and it cannot be excluded that the proportion of microcystin genotypes could be changed by reinvasion events from the sediment. In this study, Microcystis cells were always found during counting, even in winter, and the persisting population probably preserved the mcyB genotype composition from one year to the next.
Microcystin gene distribution and microcystin net production.
For the population of Lake Wannsee a close relationship between the occurrence of mcyB and microcystin net production has been observed (14, 15). In addition, 322 individual colonies sampled from numerous bodies of water in Europe were tested for mcyB gene distribution and microcystin net production by sensitive matrix-assisted laser desorption/ionization time-of-flight mass spectrometry in parallel, and only three individual colonies (1%) were found to contain mcyB but failed to show any detectable microcystin (L. Via-Odorika et al., unpublished data). Correspondingly, a significant relationship between the population growth rate and the microcystin net production rate for the same population in the summer of 2000 has been reported (15). Significant correlations between surrogate parameters such as chlorophyll a or algal biovolume and microcystin net production by Microcystis sp. have been reported by other authors as well (6, 12, 13, 18). Taken together, those results support the conclusion that it is possible to infer microcystin concentrations from surrogate parameters, for example, Microcystis cell numbers.
Over the entire study period the proportion of mcy genotypes made up the smaller part of the PC operon genotypes only. It is noteworthy that a number of field studies obtained similar results by the isolation of clonal strains from field samples. The percentage of hepatotoxic strains versus nonhepatotoxic strains isolated by Ohtake et al. (19) from Lake Kasumigaura (Japan) was 10% (n = 20), that obtained by Rohrlack et al. (24) from Lake Wannsee and Lake Pehlitzsee (Germany) was 45% (n = 22), that obtained by Shirai et al. (25) from Lake Kasumigaura was 34% (n = 68), and that obtained by Vezie et al. (30) from Lake Grand-Lieu (France) was 16% (n = 98). From the data obtained by PCR analysis of individually isolated colonies in a previous study (14), a higher mcyB genotype proportion for the total population could be expected. However, only the largest colonies of Microcystis (>100 µm) were sampled individually in the study of Kurmayer et al. (14), and it has also been found that the larger colonies (>340 µm) had the highest mcyB proportion but contributed the smallest fraction to the total cell number, <10% (15).
On the other hand, Microcystis populations have been found to vary in hepatotoxicity within a few weeks, e.g., in 50% lethal dose (expressed in milligrams [dry weight] per kilogram of body weight) (7) or in microcystin concentration (12, 21), and a mosaic structure for toxic cyanobacterial blooms has been suggested (4, 22, 30). In the present study, with one exception, short-term variation ranged from 1 to 38%, and the data on mcyB genotype proportions showed maxima on a weekly or biweekly scale (Fig. 4B). Both TNAs showed variability in cell number estimation compared to microscopy (Fig. 4D). While there was a statistically significant correlation between overestimation or underestimation of cell numbers determined in the microscope between both TNAs, the larger part of the variance remained unexplained (Kurmayer, unpublished). This unexplained part of the variance might be due to short-term shifts in mcyB genotype proportion. However, it cannot be excluded that a difference in precision between both TNAs may bias the mcyB genotype proportion in the short term, and the data on short-term variability need to be interpreted with caution. Further studies are needed to determine whether the data on short-term variability in mcyB genotype proportion are a real phenomenon or must be attributed to the noise induced by the real-time PCR approach.
This work was financed by the EU project TOPIC (Toxin Production in Cyanobacteria, CT 38-0246), the FWF project CYTOGENE (Linking Cyanotoxin Production to Genetic Diversity, P15709), and the EU project PEPCY (Toxic and Other Bioactive Peptides in Cyanobacteria, QLK4-CT-2002-02634).
|
|
|---|
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»