Previous Article | Next Article 
Applied and Environmental Microbiology, February 2002, p. 989-994, Vol. 68, No. 2
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.2.989-994.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Detection and Quantification of Pfiesteria piscicida by Using the Mitochondrial Cytochrome b Gene
Huan Zhang and Senjie Lin*
Department of Marine Sciences, University of Connecticut, Groton, Connecticut 06340
Received 24 August 2001/
Accepted 13 November 2001

ABSTRACT
Mitochondrial cytochrome
b was isolated from the dinoflagellate
Pfiesteria piscicida, and the utility of the gene for species
identification was examined. One of the primer sets designed
was shown to be highly specific for
P. piscicida. A time step
PCR protocol was used to demonstrate the potential of this primer
set for quantification of this species.

INTRODUCTION
Pfiesteria piscicida is a heterotrophic dinoflagellate that
has been implicated in fish kill events in some estuaries on
the east coast of the United States (
4,
5). Accurate identification
and quantification of
P. piscicida and differentiation of this
species from morphologically similar dinoflagellates are essential
for verifying the association of
P. piscicida with fish kills
and for studying the population dynamics of this species. Molecular
techniques that can be used for rapid and accurate identification
and quantification of this species (e.g., quantitative PCR)
are highly desirable (
13), because the traditional method that
relies on thecal plate tabulation and electron microscopy is
not feasible for processing a large number of samples.
The mitochondrial cytochrome b (cob) gene has been found to be useful for phylogenetic analysis and identification of many organisms (1, 6, 7, 8, 14, 19) because of its relatively high mutation rate (for reviews see references 2 and 16). In this study, we isolated this gene from dinoflagellates and examined its utility for identifying and quantifying P. piscicida.

Algal cultures and sample collection.
The
Pfiesteria spp. and other algae used in this study are listed
in Table
1.
Pfiesteria spp.
, Cryptoperidiniopsis sp., and
Gyrodinium galatheanum were grown in 15-practical salinity unit (PSU) seawater
that was filtered (pore size, 0.45 µm) and autoclaved;
the cryptophyte
Rhodomonas sp. strain CCMP768 was provided as
food for these organisms. The
Rhodomonas strain was grown in
15-PSU seawater amended with f/2 nutrients (
9). Other algae
were grown in f/2 medium that was prepared with 28-PSU seawater.
Illumination was provided at a photon flux of approximately
100 microeinsteins m
-2 s
-1.
View this table:
[in this window]
[in a new window]
|
TABLE 1. Algal species that were tested with P. piscicida-specific mitochondrial cytochrome b (Ppcob) and 18S rRNA (Pp18S) primers
|
When
Pfiesteria cultures were in the exponential growth phase,
feeding was discontinued for 2 to 3 days and samples were harvested
by centrifugation at 3,000
x g for 20 min at 4°C. Samples
of all other algae were collected in a similar manner from actively
growing cultures.

DNA extraction, cloning, and sequencing.
DNA was extracted by using a previously described protocol (
11),
with a slight modification of the DNA extraction buffer (0.1
M EDTA, 1% sodium dodecyl sulfate, 200 µg of proteinase
K per ml). Partial
cob (
Ppcob) and
cox3 (subunit 3 of cytochrome
oxidase,
Ppcox3) fragments were obtained by screening a
P. piscicida cDNA library (Lin and Zhang, unpublished data). Based on the
sequences of these fragments,
Ppcob- and
Ppcox3-specific primers
were designed (PPCOB and PPCOX3, respectively) (Table
2),
and
these primers were used in PCR performed with genomic DNA as
templates to amplify the 5"-terminal and upstream regions of
Ppcob. Amplification was carried out by using the following
program: 1 min at 95°C, followed by 35 cycles of 20 s at
94°C, 30 s at 58°C, and 1 min at 72°C. PCR products
were subcloned into the pCR vector (Invitrogen) and were sequenced
by using a BigDye terminator cycle sequencing kit (Applied Biosystems).
Three structurally different DNA fragments were isolated, and
each of these fragments contained
cob,
cox3, and an unidentified
sequence (X region) flanking the 5" end of
Ppcob (Fig.
1). The
first clone, PpcytbG1 (3,404 bp; GenBank accession number
AF357520),
contained the majority of the coding region of
Ppcob and the
3"-terminal region of
cox3, interrupted by a 2,332-bp X region.
The second clone, PpcytbG2 (1,861 bp; GenBank accession number
AF357521), contained the same
Ppcob and
cox3 regions as PpcytbG1
but had a deletion (1,516 bp) in the X region. In the third
fragment (PpcytbG3; 977 bp; GenBank accession number
AF357522),
almost the whole X region and 282 bp of the
Ppcob 5"-terminal
coding region were deleted, suggesting that this may be a pseudogene.
The
Ppcob sequences in the three fragments were identical.

Phylogenetic analysis.
A comparison of the
Ppcob genomic sequence with the corresponding
sequences of other organisms (Fig.
2) indicated that the
Ppcob sequence exhibited about 94% identity to the partial
cob sequence
of the dinoflagellate
Prorocentrum minimum CCMP696, 60% identity
to protozoan sequences, and <55% identity to the sequences
of other eukaryotes. In contrast, the 18S rRNA sequence of
P. piscicida exhibited 95% identity to the
P. minimum sequence
and >70% identity to the sequences of other eukaryotes.
Phylogenetic trees were constructed by using the nucleotide
sequences of both
cob and 18S rRNA and the neighbor-joining
method (
17) and were corrected by the Kimura method (
10); the
trees were rooted with the green alga
Scenedesmus obliquus (cob accession number,
X17375; 18S rRNA accession number,
AB037088).
Sequences were aligned by using the CLUSTAL W (version 1.8)
server at the DNA Data Bank of Japan (http://www.ddbj.nig.ac.jp/Welcome-e.html)
and the default values, and the alignments were adjusted by
using the Sea View program to maintain codon integrity prior
to analysis with the Phylo_win package. The reliability of the
tree topology was evaluated by bootstrap analysis with 1,000
replicates. The tree constructed by using the
cob genes of dinoflagellates
and other organisms showed clear clustering of
P. piscicida with
P. minimum (the bootstrap values for this cluster were
100%) (Fig.
2A). This tree also showed that there is a close
relationship to the apicomplexan
Plasmodium. The 18S rRNA tree
and the
cob gene tree indicated similar phylogenetic relationships
for all of the organisms (Fig.
2B). However, the branch lengths
in the
cob tree were more even than those in the 18S rRNA tree,
in agreement with the recent finding that mitochondrial genes
behave more like a molecular clock and can be ideal genes for
phylogenetic analysis (
16).

Design and test of P. piscidida-specific primers.
All of the primers used in this study are listed in Table
2.
Based on a multiple-sequence alignment of
cob nucleotide sequences
of
P. piscicida and other organisms, the following species-specific
primers were designed for
P. piscicida: PPCOBF1 (forward primer)
and PPCOBR1 (reverse primer). The other primers (forward primers
PPCOBF2 and PPCOBF3 and reverse primers PPCOBR2 and PPCOBR3)
were designed based on the sequences of the 5"-terminal and
upstream regions of
Ppcob. For comparison, the following
P. piscicida-specific primers for 18S rRNA (Pp 185) were also designed
based on the reported sequence (
13,
15): PP18SF1 and PP18SR1.
The forward primer (PP18SF1) was
Pfiesteria specific and identical
to primer 65 For designed by Rublee et al. (
15), while the reverse
primer (PP18SR1) was universal for eukaryotes. As a control
for PCR, a set of universal primers for eukaryotic 18S rRNA
(18SCOMF and 18SCOMR) was also designed (Table
2).
PCR were performed by using the specific primers and approximately 50 ng of genomic DNA from each species or strain (Table 1). Amplification was carried out by using the following program: 1 min at 95°C, followed by 35 cycles of 20 s at 94°C, 30 s at 60°C, and 40 s at 72°C. PCR products obtained with both cob and 18S rRNA were sequenced directly.
The specificities and sensitivities of the primer sets tested varied slightly. Primers PPCOB3F and PPCOB3R appeared to be the most specific (Fig. 3),
while primers PPCOB2F and PPCOB2R appeared to be the most sensitive (data not shown). Since specificity is more important, the PPCOB3F-PPCOB3R primer pair was used for further tests. The Ppcob PCR results were consistently positive when DNA templates were derived from cultures previously identified as P. piscicida (CCMP1830, CCMP1831, CCMP1834, CCMP1921, CCMP1928, and NCSU113-3) (Fig. 3A). The nucleotide sequences of all of the PCR fragments were identical, in complete agreement with the results of P. piscicida-specific 18S rRNA PCR (Fig. 3B).
The
P. piscicida-specific
cob primers never produced positive
results in cases in which the Pp18S primers gave negative results.
The
Ppcob and Pp18S primers produced no DNA product with a blank
control (H
2O) and with samples of
Rhodomonas, the
Pfiesteria-like
dinoflagellates
Cryptoperidiniopsis sp. strains CCMP1827 and
CCMP1828, and other dinoflagellates. The absence of a PCR product
in these cases was not due to poor quality of the DNA, as the
same DNA templates yielded abundant PCR products when a universal
18S rRNA primer set was used (Fig.
3C). On the contrary, the
Pp18S primers sometimes generated weak DNA bands (e.g., with
Pfiesteria shumwayae, CCMP1828, CCMP1829, and CCMP1835) (Fig.
3B), while the
Ppcob primers did not produce positive results
(Fig.
3A) and sequencing results showed that the weak bands
were nonspecific products.
The clear absence of PCR products when the Ppcob primers were used and the sequencing results obtained with the faint Pp18S PCR bands suggested that CCMP1829 and CCMP1835, two unclassified Pfiesteria-like dinoflagellates, were unlikely to be P. piscicida. In contrast, when CCMP1929, another unclassified Pfiesteria-like dinoflagellate, was examined, PCR with both Ppcob and Pp18S primers generated P. piscicida type bands (Fig. 3A and B), and the nucleotide sequences of the DNA fragments were identical to the nucleotide sequences of the DNA fragments of P. piscicida strains (GenBank accession numbers AF357519 [cob] and AF149793 [18S rRNA]). These results suggest that CCMP1929 is a P. piscicida strain.

Time step PCR to determine sensitivity and quantitation capability.
To test the ability of the
Ppcob primers to quantify
P. piscicida,
genomic DNA was isolated from 2
x 10
7 P. piscicida cells, dissolved
in 500 µl of TE buffer, and used as a stock preparation.
Serial dilutions of this DNA stock preparation equivalent to
0.2, 1, 5, 25, 125, 625, 3,125, and 15,625
P. piscicida cells
were PCR amplified by using primers PPCOB3F and PPCOB3R. Each
reaction mixture was overlaid with 1 drop of mineral oil. A
5-µl aliquot was removed from each reaction mixture at
cycles 25 and 30 during the 35-cycle PCR (called time step PCR
in this study).
The results of the time step PCR when the DNA serial dilutions were used showed that there was a good linear correlation between cell number and the quantity of the PCR product (Fig. 4).
After 35 cycles, as little as 0.2 cell was detected, and the PCR product reached a plateau when DNA from more than 125 cells was used. A linear increase in the amount of the PCR product was obtained for 5 to 625 cells with 30 cycles. When fewer cycles (25 cycles) were used, no saturation problem was encountered even for the highest number of cells (15,625 cells), although the detection limit increased to 125 cells. A wide range of cell numbers (0.2 to 15,625 cells) can be quantified by this method.

Tests with utermohls solution-preserved and field-collected samples.
We examined the potential effects on PCR of Utermohl's solution
(
18), a fixative commonly used to preserve plankton samples
collected in the field. Two
P. piscicida culture samples were
collected; one was fixed in Utermohl's solution for 2 months
at room temperature before DNA extraction, and the other was
processed immediately. The
Ppcob PCR results obtained with the
fresh samples and the samples preserved with Utermohl's solution
were similar (Fig.
5A).
Field samples were collected in December 2000 and June 2001
from the dock at the Avery Point Campus of the University of
Connecticut, which is located at the eastern end of Long Island
Sound, where no
P. piscicida has been found. Examination of
the samples by light microscopy revealed that diatoms were dominant
and that few dinoflagellates were present. Samples were fixed
with Utermohl's solution, and 50-ml subsamples were centrifuged
to obtain cell pellets. DNA was extracted and dissolved in 50
µl of TE buffer, and 1 µl of each DNA solution (equivalent
to 1 ml of a field sample) was used in a PCR performed as described
above. In parallel, the same amounts of the field DNA samples
were spiked with
P. piscicida stock DNA equivalent to the DNA
from 10 cells and analyzed. While PCR with both
Ppcob and Pp18S
primers yielded negative results for the field samples (Fig.
5B), samples spiked with
P. piscicida cells yielded positive
results, and the sensitivity was as expected.

Utility of Ppcob AND Pp18S in ecological studies.
High levels of accuracy and sensitivity are required to detect
P. piscicida in the field, especially because
P. piscicida normally
is only a minor component of the plankton community. Results
of this study showed that the
Ppcob-specific primer set consisting
of primers PPCOB3F and PPCOB3R meets these requirements. The
accuracy of this primer set was manifested by the PCR results
obtained with several
P. piscicida strains and the unclassified
P. piscicida-like dinoflagellates. The high level of sensitivity
was shown by the low detection limit, 0.2 cell ml
-1 based on
our time step PCR protocol, which is far lower than the minimum
P. piscicida concentration that causes ichthyotoxicity (more
than 250 cells ml
-1) (
3,
4) and concentrations often found in
affected estuaries (>40 cells ml
-1) (
3,
4,
12).
In the natural environment, there are numerous unknown plankton species that may sometimes produce spurious positive results in a survey based on a single gene sequence. As mentioned above, occasional nonspecific positive results were obtained in our Pp18S rRNA experiments. Although no nonspecific positive results have been found for Ppcob so far, the possibility that some nonspecific positive results may occur cannot be eliminated. Use of both Ppcob and 18S rRNA simultaneously, as demonstrated in our tests, can provide more reliable identification of P. piscicida.
In addition, often field-collected samples may need to be preserved and stored for some period of time before they can be processed. The results of this study showed that storage of samples fixed with Utermohl's solution for up to 2 months did not have noticeable effects on the PCR results compared with the results obtained with fresh unfixed samples. Finally, unknown organic materials present in field samples often coprecipitate with DNA and inhibit PCR, and variations may occur between PCR and between gels. Thus, we propose that for each field sample two subsamples should be analyzed and that to one of the subsamples a known amount of P. piscicida DNA should be added. A standard PCR in which serial dilutions of a known amount of DNA are used should also be performed simultaneously with the PCR performed with samples, and standard curves, such as those shown in Fig. 4B, should be constructed. With these precautions, the use of Ppcob and Pp18S to identify P. piscicida and the use of Ppcob to quantify P. piscicida are feasible with the simple time step PCR.

ACKNOWLEDGMENTS
We thank Robert Andersen of Bigelow Laboratory and JoAnne Burkholder
of North Carolina State University for
P. piscicida and
Pfiesteria-like
cultures, Donald M. Anderson of WHOI for the
Alexandrium tamarense culture, Patricia Tester of National Ocean Services for the
P. shumwayae culture, and Gary Wikfors of National Marine Fisheries
Service, and Patricia M. Glibert and Diane K. Stoecker of Horn
Point Environmental Laboratory for
Prorocentrum and
Gyrodinium cultures. We also thank Timothy Feinstein and Keri Costa for
assistance with growing the cultures and collecting field samples.
This research was supported by NOAA-ECOHAB grant NA86OP0491.

FOOTNOTES
* Corresponding author. Mailing address: Department of Marine Sciences, University of Connecticut, Groton, CT 06340. Phone: (860) 405-9168. Fax: (860) 405-9153. E-mail:
senjie.lin{at}uconn.edu.

ECOHAB publication number 34. 

REFERENCES
1
- Akihito, I. A., T. Kobayashi, K. Ikeo, T. Imanishi, H. Ono, Y. Umehara, C. Hamamatsu, K. Sugiyama, Y. Ikeda, Y. Sakemoto, A. Fumihito, S. Ohno, and T. Gojobori. 2000. Evolutionary aspects of gobioid fishes based upon a phylogenetic analysis of mitochondrial cytochrome b genes. Gene 259:5-15.[CrossRef][Medline]
2
- Avise, J. C. 1994. Molecular markers, natural history and evolution. Chapman & Hall, New York, N.Y.
3
- Burkholder, J. M., and H. B. Glasgow, Jr. 1997. Pfiesteria piscicida and other Pfiesteria-like dinoflagellates: behavior, impacts, and environmental controls. Limnol. Oceanogr. 42:1052-1055.
4
- Burkholder, J. M., H. B. Glasgow, and N. Deamer-Melia. 2001. Overview and present status of the toxic Pfiesteria piscicida (Dinophyceae). Phycologia 40:186-214.
5
- Burkholder, J. M., E. J. Noga, C. H. Hobbs, and H. B. Glasgow, Jr. 1992. New "phantom" dinoflagellate is the causative agent of major estuarine fish kills. Nature 358:407-410.[CrossRef][Medline]
6
- Callejas, C., and M. D. Ochando. 2000. Recent radiation of Iberian barbel fish (Teleostei, Cyprinidae) inferred from cytochrome b genes. J. Hered. 91:283-288.[Abstract/Free Full Text]
7
- Conroy, C. J., and J. A. Cook. 1999. MtDNA evidence for repeated pulses of speciation within arvicoline and murid rodents. J. Mammal. Evol. 6:221-245.[CrossRef]
8
- Conway, D. J., C. Fanello, J. M. Lloyd, B. M. A. Al-Joubori, A. H. Baloch, S. D. Somanath, C. Roper, A. M. J. Oduola, B. Mulder, M. M. Povoa, B. Singh, and A. W. Thomas. 2000. Origin of Plasmodium falciparum malaria is traced by mitochondrial DNA. Mol. Biochem. Parasitol. 111:163-171.[CrossRef][Medline]
9
- Guillard, R. R., and J. H. Ryther. 1962. Studies of marine planktonic diatoms. I. Cyclotella nana Hustedt and Detonula confervacea (Cleve.) Gran. Can. J. Microbiol. 18:229-239.
10
- Kimura, M. 1980. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. E vol. 16:111-120.
11
- Lin, S., and E. J. Carpenter. 1999. A PSTTLRE-form of cdc2-like gene in the marine microalga Dunaliella tertiolecta. Gene 239:39-48.[CrossRef][Medline]
12
- Marshall, H. G., D. Seaborn, and J. Wolny. 1999. Monitoring results for Pfiesteria piscicida and Pfiesteria-like organisms from Virginia waters in 1998. Va. J. Sci. 50:287-298.
13
- Oldach, D., C. F. Delwiche, K. S. Jakobsen, T. Tengs, E. G. Brown, J. W. Kempton, E. F. Schaefer, H. A. Bowers, H. B. Glasgow, Jr., J. M. Burkholder, K. A. Steidinger, and P. A. Rublee. 2000. Heteroduplex mobility assay-guided sequence discovery: elucidation of the small subunit (18S) rRNA sequences of Pfiesteria piscicida and related dinoflagellates from complex algal culture and environmental sample DNA pools. Proc. Natl. Acad. Sci. USA 97:4303-4308.[Abstract/Free Full Text]
14
- Pearson, O. P., and M. F. Smith. 1999. Genetic similarity between Akodon olivaceus and Akodon xanthorhinus (Rodentia: Muridae) in Argentina. J. Zool. 247:43-52.[CrossRef]
15
- Rublee, P. A., J. Kempton, E. Schaefer, J. M. Burkholder, H. B. Glasgow, and D. W. Oldach. 1999. PCR and FISH detection extends the range of Pfiesteria piscicida in estuarine waters. Va. J. Sci. 50:325-326.
16
- Saccone, C., C. Gissi, C. Lanave, A. Larizza, G. Pesole, and A. Reyes. 2000. Evolution of the mitochondrial genetic system: an overview. Gene 261:153-159.[CrossRef][Medline]
17
- Saitou, N., and M. Nei. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. E vol. 4:406-425.
18
- Utermohl, H. 1958. Zur Vervollkommnung der quantitativen Phytoplankton-Methodik. Mitt. Int. Ver. Theor. Angew. Limnol. 9:1-38.
19
- Zhang, H., N. Mikawa, Y. Yamada, N. Horie, A. Okamura, T. Utoh, S. Tanaka, and T. Motonobu. 1999. Detection of foreign eels in the natural waters of Japan by PCR. Fish. Sci. 65:684-686.
Applied and Environmental Microbiology, February 2002, p. 989-994, Vol. 68, No. 2
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.2.989-994.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
This article has been cited by other articles:
-
Zhang, H., Litaker, W., Vandersea, M. W., Tester, P., Lin, S.
(2008). Geographic distribution of Karlodinium veneficum in the US east coast as detected by ITS-ferredoxin real-time PCR assay. J PLANKTON RES
30: 905-922
[Abstract]
[Full Text]
-
Zhang, H., Bhattacharya, D., Maranda, L., Lin, S.
(2008). Mitochondrial cob and cox1 Genes and Editing of the Corresponding mRNAs in Dinophysis acuminata from Narragansett Bay, with Special Reference to the Phylogenetic Position of the Genus Dinophysis. Appl. Environ. Microbiol.
74: 1546-1554
[Abstract]
[Full Text]
-
Saito, K., Drgon, T., Krupatkina, D. N., Drgonova, J., Terlizzi, D. E., Mercer, N., Vasta, G. R.
(2007). Effect of Biotic and Abiotic Factors on In Vitro Proliferation, Encystment, and Excystment of Pfiesteria piscicida. Appl. Environ. Microbiol.
73: 6410-6420
[Abstract]
[Full Text]
-
Nash, E. A., Barbrook, A. C., Edwards-Stuart, R. K., Bernhardt, K., Howe, C. J., Nisbet, R. E. R.
(2007). Organization of the Mitochondrial Genome in the Dinoflagellate Amphidinium carterae. Mol Biol Evol
24: 1528-1536
[Abstract]
[Full Text]
-
Lin, S., Zhang, H., Dubois, A.
(2006). Low abundance distribution of Pfiesteria piscicida in Pacific and Western Atlantic as detected by mtDNA-18S rDNA real-time polymerase chain reaction. J PLANKTON RES
28: 667-681
[Abstract]
[Full Text]
-
Zhang, H., Lin, S.
(2005). Development of a cob-18S rRNA Gene Real-Time PCR Assay for Quantifying Pfiesteria shumwayae in the Natural Environment. Appl. Environ. Microbiol.
71: 7053-7063
[Abstract]
[Full Text]
-
Coyne, K. J., Burkholder, J. M., Feldman, R. A., Hutchins, D. A., Cary, S. C.
(2004). Modified Serial Analysis of Gene Expression Method for Construction of Gene Expression Profiles of Microbial Eukaryotic Species. Appl. Environ. Microbiol.
70: 5298-5304
[Abstract]
[Full Text]
-
Lin, S., Zhang, H.
(2003). Mitogen-Activated Protein Kinase in Pfiesteria piscicida and Its Growth Rate-Related Expression. Appl. Environ. Microbiol.
69: 343-349
[Abstract]
[Full Text]
-
Saito, K., Drgon, T., Robledo, J. A. F., Krupatkina, D. N., Vasta, G. R.
(2002). Characterization of the rRNA Locus of Pfiesteria piscicida and Development of Standard and Quantitative PCR-Based Detection Assays Targeted to the Nontranscribed Spacer. Appl. Environ. Microbiol.
68: 5394-5407
[Abstract]
[Full Text]