Previous Article | Next Article ![]()
Applied and Environmental Microbiology, April 2007, p. 2552-2560, Vol. 73, No. 8
0099-2240/07/$08.00+0 doi:10.1128/AEM.02389-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

School of Plant Science, University of Tasmania, Private Bag 55, Hobart, Tasmania 7001, Australia,1 School of Aquaculture, University of Tasmania, Locked Bag 1-370, Launceston, Tasmania 7250, Australia2
Received 10 October 2006/ Accepted 12 February 2007
|
|
|---|
|
|
|---|
Pfiesteria and PLDs are morphologically defined by their distinct Kofoidian thecal plate formulae (37). The morphological characterization of C. brodyi and Pfiesteria species requires thecal plate analysis by scanning electron microscopy (SEM) on cultured strains in combination with a membrane stripping process (24) that is not suitable for rapid sample processing. While fluorescent in situ hybridization might be a useful alternative, this method is labor intensive and time consuming and suffers from interference by nonspecific binding in complex mixtures (32).
Real-time quantitative PCR that incorporates fluorogenic 5' nuclease chemistry (TaqMan) is a technique that offers highly sensitive, rapid, and quantitative analysis. The TaqMan system uses sequence-based fluorogenic probes to detect amplified products during PCR cycling, and the detection of PCR amplicons is monitored by measuring the increase in fluorescence. This approach has been used successfully to detect small-subunit ribosomal DNA (SSU rDNA) of P. piscicida in mixed algal cultures; however, there is a major concern with regard to reliable quantification of real-time PCR assays from environmental water samples. Environmental samples often contain organic and inorganic substances with the potential to inhibit PCR assays (18, 42). Thus, the effective removal of PCR inhibitors from samples has to be considered for any field-based survey using real-time PCR assays.
The present study developed a C. brodyi-specific real-time PCR using a TaqMan probe based on internal transcribed spacer 2 (ITS2) of the rDNA. For detection of P. piscicida and P. shumwayae, existing P. piscicida-and P. shumwayae-specific real-time PCR assays (3, 42) were used and tested against different genotypes from Australia. Methods for overcoming PCR-inhibitory substances in environmental samples were examined for application of the real-time PCR protocol with water samples. Temporal occurrences of P. piscicida, P. shumwayae, and C. brodyi were characterized during an 18-month field survey using species-specific real-time PCR assays in a eutrophic Tasmanian estuary.
|
|
|---|
f/2 medium (11a) at 20°C in the dark. Rhodomonas salina (CS-24; obtained from CSIRO Marine Laboratories, Hobart, Australia) was supplied to C. brodyi and Pfiesteria species as food. All autotrophic cultures listed in Table 1 were obtained from the University of Tasmania's collection of microalgae. These were grown in 28
f/2 medium at 20°C with cool white fluorescent lamps on a 12:12-h light:dark cycle. An additional eight heterotrophic strains (CBCM1 to CBCM8, unidentified Pfiesteria look-alike cultures) were collected by Wasele Hosja from the Chapman River, Western Australia, during September 2005, and their identity was assessed by P. piscicida-, P. shumwayae-, and C. brodyi-specific real-time PCR assays. |
View this table: [in a new window] |
TABLE 1. Cultures used in this study
|
![]() View larger version (16K): [in a new window] |
FIG. 1. Locations where water samples were collected for this study.
|
PCR inhibitor removal from environmental samples.
Three different methods were compared for their capacity to remove PCR inhibitors from environmental samples: PVPP (polyvinylpolypyrrolidone), PVP360 (polyvinylpyrrolidone), and a 100-fold dilution of the DNA extract, using nontreated samples as a control. Water samples collected on 15 February 2005 were tested with real-time PCR to confirm the absence of C. brodyi and then used to detect PCR inhibitors. The following treatments were applied to water samples without natural C. brodyi, spiked with 1 µl of CBDE10 DNA (1.2 ng µl1): (i) serial DNA dilutions of water sample spiked with CBDE10 DNA, (ii) water sample spiked with CBDE10 DNA for treatment with PVPP, and (iii) water sample spiked with CBDE10 DNA for treatment with PVP360. DNA was extracted from the cell pellets using a DNeasy tissue kit. ATL (180 µl) and proteinase K (20 µl) were added and thoroughly mixed by vortexing. Following incubation of the cell lysates for 1 h at 55°C, PVPP or PVP360 was added to a final concentration of 2% and the cell lysates were incubated for 30 min at 55°C. The cell lysates were then processed according to the manufacturer's instructions. C. brodyi-specific real-time PCR was performed with DNA from these treatments. The absence of PCR inhibitors in the 100-fold DNA dilution of all three treatments was confirmed. The absence of PCR inhibitors in the 100-fold DNA dilution was also confirmed with 10 other water samples that were randomly chosen among the field samples. The results of the 100-fold dilution of the DNA were calculated into cell estimates for the quantitative component of the assay. When large volumes of cell pellets (over 10 mg) were processed for DNA extraction, the DNeasy kit reagent volume was increased twice.
Design of TaqMan probe and primer set for detection of C. brodyi.
C. brodyi and related dinoflagellate rDNA sequences available from GenBank were aligned using the program ClustalX (40). Manual searches of the alignments were conducted to determine unique sequences and to develop a C. brodyi-specific real-time PCR assay. The sequences for the primer-probe set were conserved within C. brodyi strains, but discriminated between C. brodyi and other dinoflagellates. To minimize the variation in assay sensitivity caused by intraspecies variations within C. brodyi strains, mixed-base degenerate primer codes were used (Table 2). The primer and probe sequences of the target species were analyzed with Primer 3 (Whitehead Institute and Howard Hughes Medical Institute, Maryland) and Oligo analyzer 3.0 (Integrated DNA Technologies, Inc., Iowa) software for optimal melting temperature and secondary structure, and subsequently, the primers and probe were synthesized by Geneworks Pty. Ltd. (Hindmarsh, South Australia, Australia). The probe was dual-labeled with the fluorescent dyes FAM and 6-carboxytetramethyl-rhodamine (TAMRA) at the 5' and 3' ends, respectively.
|
View this table: [in a new window] |
TABLE 2. Primers and TaqMan probes for species-specific real-time PCR
|
Optimization of primer and probe concentrations.
To achieve optimal performance, a dilution series of primers and probe concentrations was tested by real-time PCR assay. The testing primer and probe concentrations ranged from 100 to 500 nM and 100 to 300 nM, respectively. The primer and probe concentrations that provided the lowest CT value and the highest total increase of fluorescence compared to the background were selected. The primer and probe concentrations were 0.2 µM and 0.15 µM, respectively.
Real-time PCR assay conditions.
PCR assays using TaqMan probes were performed on a Rotor-Gene RG-3000A (Corbett Research, New South Wales, Australia). The following reagents were added for a 20-µl reaction mixture: 10 µl of platinum quantitative PCR supermix-UDG (Invitrogen, New South Wales, Australia), forward and reverse primers each at a final concentration of 0.2 µM, fluorogenic probe at a final concentration of 0.15 µM, 2 µl of template DNA, and PCR-grade water to a final volume of 20 µl. For PCR assays using SYBR green, the following components were added for a 20-µl reaction mixture: 10 µl of platinum SYBR green real-time PCR supermix-UDG (Invitrogen) and concentrations of forward and reverse primers, template DNA, and PCR-grade water as described above. The thermal cycling conditions consisted of 2 min at 50°C and 2 min at 95°C followed by 45 cycles (50 cycles for SYBR green) of 10 s at 95°C and 45 s (40 s for SYBR green) at 60°C. Fluorescence data were collected at the end of each cycle, and determination of the cycle threshold line was carried out automatically by the instrument. PCR amplicons that were positive with SYBR green real-time PCR were further analyzed by gel electrophoresis through 2% agarose gels to confirm the amplicon size.
Construction of standard curves.
Standard curves using cell numbers were constructed from C. brodyi (CBWA12) and P. piscicida (PPSB25) cultures. Cell numbers (55,000 cells for C. brodyi and 11,000 cells for P. piscicida; n = 1) were estimated by LM in a hemocytometer (Blau Brand, Germany) before harvesting the cells, and genomic DNA was extracted. Tenfold serial dilutions of the DNA extracts were used to construct the standard curve (duplicate measurements by real-time PCR). The cell number of the target species in field samples was determined from CT values from comparison with the standard curve. To evaluate the accuracy and precision of the standard curves, known concentrations of C. brodyi cells (CBWA12; 480, 4,800, and 48,000 cells; n = 3) were spiked with sterile, filtered, environmental samples (0.2-µm membrane filter; MFS, California) followed by DNA extraction, 100-fold dilution of the DNA extracts, and real-time PCR. The results of the spiked experiments were then compared to the standard curves.
Species verification of PCR amplicons from environmental samples by sequence analysis.
PCR amplicons obtained in November 2004 to May 2005 that gave positive real-time PCR results (seven PCR products with C. brodyi-specific primers and one product with P. piscicida-specific primers) were visualized on an agarose gel stained with ethidium bromide. The band on the agarose gel was reamplified by use of the "band stab" approach (2). The band was stabbed with a stainless steel syringe needle, and the needle was dipped into the PCR mixture. PCR assays were carried out immediately, and the PCR products were purified. The purified PCR amplicons were cloned using the pGEM-T vector system (Promega, Madison, WI), and then one clone of each PCR product was sequenced with pUC/M13 primers (Promega, Madison, WI) on a Beckman-Coulter CEQ2000 capillary electrophoresis sequencer.
Phylogenetic analysis.
Partial ITS and SSU rDNA sequences of environmental samples were aligned with sequences of closely related organisms available on GenBank by using ClustalX (40). The dinoflagellate Karlodinium veneficum was chosen as the outgroup taxon, and the final alignment was refined manually. Bayesian analysis (likelihood) was carried out using the computer program MrBayes v. 3.12 (13). Nucleotides were treated as independent unordered characters of equal weight, and gaps were interpreted as missing data. The Bayesian analysis was run for 150,000 generations, using a 4-by-4 nucleotide-substitution, general time-reversible model with gamma-distributed among-site rate variation and a proportion of invariable sites (GTR+G+I). The resulting tree was sampled every 100 generations.
Microscopic observation.
The presence of Pfiesteria or Pfiesteria-like dinoflagellates in environmental samples collected every 2 weeks was checked by LM before the samples were stored at 70°C. Microscopic observations were compared to the real-time PCR results.
Nucleotide sequence accession numbers.
The sequences of field samples DE250505, DE281104, DE190105, DE140305, DE280305, DE120405, DE270405, and DE160505 have been deposited in GenBank under accession numbers DQ991364 to DQ991371, respectively. Sample codes: DE, Derwent River in Tasmania; the code for the sampling station (DE) is followed by a number indicating the day, month, and year of the sampling event (e.g., DE250505, Derwent River sampled on 25 May 2005).
|
|
|---|
![]() View larger version (72K): [in a new window] |
FIG. 2. Agarose gel analysis showing C. brodyi-selective real-time PCR products. The real-time PCR was carried out using a SYBR green-based system with the primers CBITSF and CBITSR. The arrow indicates a positive real-time PCR product of 128 bp. Lanes: 1, 2-kbp-ladder molecular size marker; 2, C. brodyi; 3, P. piscicida; 4, P. shumwayae; 5, Karlodinium veneficum (=micrum); 6, Karenia mikimotoi; 7, Prorocentrum minimum; 8, no-template control.
|
The specificities of the P. piscicida and P. shumwayae assays (SSU rDNA-based primers and mitochondrial cytochrome b-based primers, respectively) were tested with Australian genotypes of C. brodyi. No nonspecific reactions were produced by these primers. An environmental sample from the Derwent River samples assayed with P. piscicida-specific primers was also sequenced and was identical to the documented P. piscicida sequence.
Phylogenetic analysis.
PCR amplicons generated from 7 environmental samples with C. brodyi-specific primers were identical, or differed by only 1 bp from each other, and differed from C. brodyi CBWA12 by 2.3 to 3.1%. The sequence comparison between C. brodyi and one of the field samples is shown in Fig. 3. Bayesian analysis of environmental samples and other related species showed that the sequences of the PCR amplicons clustered in a clade with C. brodyi at the 51 to 99% bootstrap values (Fig. 4). C. brodyi, P. piscicida, and other related dinoflagellates were rooted as distinct clades.
|
View larger version (4K): [in a new window] |
FIG. 3. Alignment of partial ITS2 rDNA sequences from environmental sample and C. brodyi (CBWA12). The field sample DE160505 was sampled from the Derwent River on 16 May 2005 and was detected by C. brodyi-specific real-time PCR. Dots indicate nucleotides identical to those of C. brodyi.
|
![]() View larger version (31K): [in a new window] |
FIG. 4. Phylogram based on Bayesian analysis of partial ITS2 rDNA sequences (128 bp) of environmental samples and closely related organisms. The environmental samples DE281104 to DE160505 were sampled from the Derwent River from 28 November 2004 to 16 May 2005. These samples were detected by C. brodyi-specific real-time PCR and were sequenced. The sample DE250505 was collected on 25 May 2005 and amplified with P. piscicida-specific real-time PCR. The numbers adjacent to each node represent probabilities which were derived from 150,000 generations using a general time-reversible evolution model with gamma-distributed (0.248) among-site rate variation. The strains sequenced in the present study are shown in bold. Consistency index, 0.837; rescaled consistency index, 0.696; retention index, 0.831; homoplasy index, 0.162.
|
0.997 (Fig. 5). These assays could detect less than one cell of the target species in a reaction. When 100-fold dilutions of DNA extracts from C. brodyi cultures (480, 4,800, and 48,000 cells) spiked with environmental samples were compared to the standard curves, concentrations of the spiked cells (4.0 ± 0.9, 41 ± 10, and 564 ± 99 cells, respectively; n = 3; P > 0.05 by Student's t test) did not significantly differ from those of standard curves (4.8, 48, and 480 cells). The DNA content in C. brodyi cells was also measured by using a PicoGreen dsDNA quantification kit, and the amount of DNA was 7,274 ± 1,254 fg (mean ± standard deviation) per cell.
![]() View larger version (33K): [in a new window] |
FIG. 5. Linear relationship between the CT values and the cell numbers for C. brodyi (r2 = 0.997) (A) and for P. piscicida (r2 = 0.998) (B). The standard errors from two measurements are shown as error bars.
|
|
View this table: [in a new window] |
TABLE 3. Real-time PCR efficiencies for Derwent River samples (15 February 2005) with PVPP or PVP360 and different dilutions of the samples
|
![]() View larger version (16K): [in a new window] |
FIG. 6. Temporal variations in C. brodyi, P. piscicida, and P. shumwayae abundances in the Derwent River, Tasmania, during 2004 to 2006, quantified by species-specific real-time PCR. The values are the means of triplicate wells.
|
|
View this table: [in a new window] |
TABLE 4. C. brodyi abundance in Derwent River from November 2004 to April 2006a
|
|
|
|---|
When environmental samples are examined to obtain quantitative results, one main problem is false-negative results caused by real-time PCR inhibitors coextracted with target DNA. Surface water contains phenolic compounds, humic acids, and heavy metals, all known to be PCR inhibitors (41). The efficiency of removal of PCR inhibitors is generally dependent on the sample purity and volume. PVPP and PVP are known to remove humic compounds from DNA (10). In the present study, a dilution step and chemical treatments to remove these real-time PCR inhibitors were tested. The chemical treatment might not be sufficient to remove impurities from 10-liter Derwent River samples that might interfere with the real-time PCR assay. Processing low sample volumes and increasing PCR volumes may help to increase the DNA extraction efficiency and reduce concentrations of PCR inhibitors, so that overall sensitivity may be increased. However, where target organisms are present in low abundance in environmental water samples, stochastic effects on cell quantification become significant, by reducing the amount of the water sample being processed. Therefore, it may be a better strategy to sacrifice assay sensitivity than to process smaller volumes of water. However, an acceptable compromise should be decided by the objective of the study and potential inhibition should be checked by parallel analysis of natural samples with and without spiked C. brodyi cells.
Specificity of real-time PCR assay.
Previous studies have demonstrated that a primer set used in the SYBR green system can be used in a standard PCR assay to determine the presence or absence of target DNA (1). The primer set used in this study may also be applied in standard PCR assay format. The SYBR green and TaqMan probe systems have been known to have similar assay sensitivities, but with different levels of specificity (12, 16). In this study, no non-C. brodyi strain produced a signal in either the SYBR green-based or the TaqMan probe-based real-time PCR. However, when real-time PCR is applied to field samples, uncertainty in the assay specificity can result from the potential presence of unknown related species with similar sequences to the primer sites. Therefore, sequencing analysis of positive real-time PCR amplicons is required as further confirmation of assay specificity. More-variable DNA sequences, such as ITS rDNA regions, are useful genetic markers for the development of highly accurate species-specific real-time PCR. Homology searches of a nucleotide database revealed that the ITS2 rDNA region contains sequences unique to the target molecule, to which a primer-probe set of the real-time PCR was designed. Given the higher variability of the ITS regions compared to the SSU, LSU, and 5.8S genes, ITS sequence data have been used to analyze other dinoflagellates at the interspecies level (8, 21, 27). The higher divergence means that designing real-time PCR primer-probes based on ITS regions may reduce false positives.
Sequence and phylogenetic analyses.
One copy of an ITS2 rDNA fragment from each water sample was sequenced for confirmation of assay specificity. Three nucleotides differed from C. brodyi (CBWA12) DNA, and the sequenced samples were all identified as identical to C. brodyi by phylogenetic analysis. The sequence variation in Fig. 3 could be a result of intraclonal variation between rDNA loci in single C. brodyi cells, because only one clone was sequenced and the sequences were not obtained by consensus sequencing. Phylogenetic analyses of C. brodyi and other related dinoflagellates using ITS-5.8S rDNA sequences were previously described (23). When the partial ITS-based phylogram from this study was compared with the entire ITS-5.8S-based phylogram, the order of divergences among this dinoflagellate grouping differed slightly. Although the partial ITS-based phylogram clearly clustered the Derwent River samples (GenBank no. DQ991365 to DQ991371) with C. brodyi, the phylogram from this study may not resolve the phylogeny of C. brodyi and related dinoflagellates, because sort sequences (128 bp) were used to generate this phylogenetic tree.
Seasonal changes in C. brodyi, P. piscicida, and P. shumwayae in Derwent River, Tasmania.
Quantification of cells by real-time PCR could be conducted either by using the absolute number of gene copies in the target species or by relative concentrations of the target species. A single P. piscicida cell has been reported to have 100 to 200 rDNA copies (34). However, the number of rDNA copies in C. brodyi and P. shumwayae is not yet known. Therefore, the abundances of these species in environmental samples were calculated by comparisons with a standard curve generated from serial dilutions of DNA. The rDNA copy number in the natural populations was assumed to be the same as in the culture used as the standard.
The 18-month survey from the Derwent River showed that vegetative cells of C. brodyi were regularly present at very low abundance (mostly below 25 cells liter1) but were slightly more abundant (55 to 112 cells liter1) in April and May. Accurate quantification of C. brodyi in field samples has not been previously reported, due to lack of a quantitative molecular tool. Similarly low abundances were also shown for Pfiesteria species, which were detected by real-time PCR or PCR-fluorescent fragment detection techniques (9, 18, 42).
Zoospore sizes of C. brodyi range from 6 to 11 µm in length and from 5 to 10 µm in width (38). P. piscicida has been reported to have vegetative cell sizes ranging from 8.5 to 12 µm (20). In this study, a phytoplankton net (5-µm open mesh area) was used to concentrate large volumes (10 liters) of environmental samples for rapid sample processing. The mesh size was close to the size of the zoospores, and thus, it may allow some C. brodyi and P. piscicida cells to pass through the mesh. Centrifugation of the cells at 4,000 x g for 10 min also might be somewhat insufficient to settle all the cells. This could have affected the seasonal abundance results presented in Fig. 6 and underestimated the cell numbers, although this is unlikely to alter the conclusion of low abundance of C. brodyi and P. piscicida. Processing water samples through filtration with glass microfiber filters could minimize the potential loss of the organisms during sample collection.
Prey-dependent growth has been shown for C. brodyi and related dinoflagellates (6, 17, 26, 35). These dinoflagellates are known to feed on various microalgae, bacteria, protozoan ciliates, rotifers, tissues of fish, and planktonic shellfish larvae (5, 7, 35, 36). When preferred prey, such as cryptophytes, including R. salina, are supplied at high rates, higher cell densities can be obtained for these dinoflagellates (17, 26). Size, shape, and/or chemical cues of prey may affect food preferences of heterotrophic dinoflagellates (14, 25). Field experiments showed that the relationship between P. piscicida and total chlorophyll a concentrations may be influenced by appropriate prey rather than the abundance of total phytoplankton (9, 18). The seasonal changes in C. brodyi in the Derwent River thus may be associated with availability of suitable prey at times, implying that prey co-occurred in all seasons and prey cell densities increased during April and May, when Tasmanian coastal seawater reaches its annual temperature maximum. A positive correlation between a heterotrophic dinoflagellate and prey abundances has been shown in a previous study (15). When cell densities of Stoeckeria algicida (morphologically similar to Shepherd's Crook) peaked at 17,400 cells ml1 in July, being the warmer season in Korea, cell densities of Heterosigma akashiwo (algal prey for S. algicida) decreased, indicating that S. algicida populations may have considerable grazing impact on H. akashiwo and contribute to the decline of H. akashiwo blooms (15). However, the microscopic counts used in that study were inadequate for discriminating morphologically similar species and may erroneously have included cells that resemble S. algicida, such as C. brodyi, Pfiesteria spp., and Karlodinium veneficum. As the protocol developed in this study allowed for the estimation of C. brodyi abundances from environmental samples, the use of real-time PCR will provide higher detection accuracies of the marine dinoflagellates in field-based studies than traditional methods.
Eight heterotrophic cultures identified by C. brodyi-specific real-time PCR.
Identification of eight cultures from the Chapman River, Western Australia, was confirmed by species-specific real-time PCR. This assay was completed in approximately 1.5 h. Therefore, this assay would be useful for the rapid and sensitive identification of clonal cultures, as well as environmental samples. However, for confirmation of the morphological and genetic characters of these cultures, SEM and sequencing analyses are required. The finding of C. brodyi in the Chapman River confirms that this species is widely distributed and commonly present in Australia. P. shumwayae has been detected previously on every continent by using molecular probes based on SSU rDNA, suggesting that P. shumwayae is a cosmopolitan species (4, 28, 29, 30, 31, 32, 33, 39). In Australia, this species was detected in various locations in Queensland, South Australia, Victoria, Western Australia, and Tasmania (P. A. Rublee, personal communication). However, no P. shumwayae was detected by mitochondrial cytochrome b-based probes during this study. The discrepancy between previous findings and this study raises a question about whether this is due to the limited sampling area of this study or false-positive reactions by the SSU rDNA-based probes. Companion work testing the specificity of the P. shumwayae-selective probes on C. brodyi and related species will be published elsewhere (25a).
In conclusion, the real-time PCR assay developed in this study demonstrated high sensitivity and specificity to detect and quantify C. brodyi from environmental samples. A real-time PCR assay was successfully implemented for the ecological study of this dinoflagellate. Our real-time PCR approach readily detected low abundances of C. brodyi (1 to 112 cells liter1), revealing that C. brodyi is generally present as a minor component of phytoplankton in the Derwent River and that Pfiesteria species are rare in Tasmania. The detection of C. brodyi and P. piscicida at such low concentrations also suggests that real-time PCR, with appropriate adaptation of sample extraction, is an effective tool for investigating the distributions and seasonal dynamics of dinoflagellates in the marine environment.
This work was supported by a grant from the Internal Research Grant Scheme (IRGS), University of Tasmania.
Published ahead of print on 23 February 2007. ![]()
|
|
|---|
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»