Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2009, p. 54-63, Vol. 75, No. 1
0099-2240/09/$08.00+0 doi:10.1128/AEM.00818-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Ocean Science Department, University of California, Santa Cruz, 1156 High Street, Santa Cruz, California 95064,1 Moss Landing Marine Laboratories, 8272 Moss Landing Rd., Moss Landing, California 950392
Received 9 April 2008/ Accepted 6 November 2008
|
|
|---|
|
|
|---|
There are three confirmed producers of the YTX phycotoxin: the dinoflagellates Protoceratium reticulatum (Claparhde and Lachmann) Buetschli, Lingulodinium polyedrum (Stein) Dodge, and Gonyaulax spinifera Dodge. YTX production within and among dinoflagellate species tested to date is highly variable. Previous reports have identified both nontoxic and toxic isolates of L. polyedrum from Spain, the United Kingdom, and California (3, 18, 43, 47, 58). In addition, several YTX-positive isolates of L. polyedrum have been identified in Italy (9, 62) and Ireland (J. Silke, personal communication), yet only nontoxic isolates have been observed in Norway waters (44). Several nontoxic isolates of G. spinifera have been identified from the United Kingdom (21, 58) and one from the United States (18), yet several toxic isolates from New Zealand (45) and Italy (46) have been observed. Water samples from New Zealand yielded both nontoxic and toxic isolates of G. spinifera, including one isolate that produced the highest amount of YTX per cell for any genus to date (45; L. Rhodes, personal communication). YTX-producing isolates of P. reticulatum have been described from Canada (58), Italy (5), Norway (50), Japan (10, 52), New Zealand (30, 34, 51, 52), Spain (41, 42), South Africa (12, 21), and the United Kingdom (58). Together, these studies indicate that while YTX production by P. reticulatum appears to be globally distributed, a large variation in the toxicity (0.3 to 79 pg YTX cell–1) is observed between isolates of this species. To add to the challenges of establishing monitoring criteria for potential YTX-producing dinoflagellates, there have been conflicting studies published on the toxicity of P. reticulatum from U.S. coastal waters, with reports indicating both the presence and absence of YTX accumulation in cultured isolates from California, Florida, and Washington (41, 42; B. Paz, personal communication).
During environmental YTX events, many dinoflagellates (including species not previously identified to produce YTX) have been isolated, grown in culture, and tested for the production of YTX, particularly members of the Prorocentrum genus which frequently co-occur during these events. This process is labor intensive and dependent on the isolation of a number of cells sufficient for YTX analysis. Given the variation in the expression of toxicity among isolates of the same species and the limited morphological variability within species, we attempted to identify coarse scale genetic markers that would differentiate between spatially and temporally distinct isolates while simultaneously identifying YTX-accumulating strains of L. polyedrum, P. reticulatum, and G. spinifera.
Given the worldwide distribution of potentially toxic dinoflagellates and the global observations of YTX, we speculated that phylogenetic analysis of molecular sequence information from these YTX-producing dinoflagellates could be used as a tool to help identify other potential yet unidentified YTX-producing dinoflagellates. The resultant molecular database could also enhance the subsequent development of robust detection probes for monitoring YTX-producing species.
The use of molecular phylogenetics is now widely utilized in addition to traditional methodologies (such as morphology, ultrastructure, life-cycle information, and fossil record) to understand the evolutionary history of dinoflagellates and to evaluate the relationships within genera among morphologically indistinguishable species. The rRNA operon comprises the genome regions that code for the RNA components of the ribosomes and consists of the large subunit (LSU), the small subunit (SSU), and the 5.8S genes, the latter bound by internal transcribed spacer regions 1 and 2 (ITS1 and ITS2). The rRNA gene has been the target of many phylogenetic studies because ribosomes are universally present in living organisms and functional constraints have resulted in high sequence conservation within these domains (13, 26, 53). The LSU consists of many structural domains (D1 to D12), and the D1, D3, and D8 domains are particularly useful to evaluate the phylogenetic relationships of closely related species (27, 32, 54) since these variable regions flank more conservative regions. The ITS regions are commonly used to analyze closely related and geographically different species (6, 24, 28, 48, 53, 56). The rRNA operon has also been widely and successfully targeted for the design of primers and probes to detect or quantify harmful algal bloom species (2, 31, 32, 55).
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. List of the cultures analyzed, including the label of each, the GenBank accession number, the isolation location, the year of isolation, and the culture collection or person from which the culture was received
|
Duplicate samples of each culture were gently filtered onto Poretics, 5.0-µm polycarbonate filters (Osmonics, Inc.), at mid-exponential growth phase. The cells were resuspended in medium and centrifuged in a Fisher Scientific Marathon 8K centrifuge for 3 min at 1,073 x g (4,000 rpm). The supernatant was pipetted off, and the pellets were immediately frozen in liquid nitrogen and stored at –80°C until analysis.
DNA extraction and PCR amplification.
The genomic DNA was extracted from the pellets using NucleoSpin plant kits (BD Biosciences) according to the manufacturer's instructions. Approximately 50 ng of genomic DNA was amplified using 50-µl PCRs containing nuclease-free water (Fisher Scientific), 1x JumpStart PCR buffer (10 mM Tris-HCl [pH 8.3], 50 mM KCl; Sigma-Aldrich), 2.5 mM MgCl, 0.2 mM deoxynucleoside triphosphate (Advantage ultrapure PCR deoxynucleotide mix; BD Biosciences), 200 nM of gene-specific oligonucleotide primers, and 0.5 U JumpStart Taq DNA polymerase (Sigma-Aldrich). The ITS1, 5.8S, and ITS2 regions (approximately 800 bp) were amplified from each culture using several primers, each of which is listed in Table 2. The Lp1F and Lp2R oligonucleotide pairs were based on the L. polyedrum ribosomal DNA (rDNA) sequence AF377944 and designed to amplify the entire ITS1-5.8S-ITS2 region from the 3' end of the SSU rDNA and 5' end of the LSU rDNA. This dinoflagellate-specific primer set was used successfully to amplify this domain from isolates CCMP404, CCMP1889, and CTCC01 (all P. reticulatum) and CCMP1931, CCMP1936, 104A, and CCAP1121/2 (all L. polyedrum). The more generic oligonucleotide pair 1400F and 38R targeting the same domain (29) were needed to amplify the entire ITS1-5.8S-ITS2 domain from the G. spinifera isolate CCMP409. All amplifications were carried out in duplicate using the Perkin Elmer GeneAmp PCR system 2400 and the following cycling conditions: initial template denaturation for 5 min at 94°C, followed by 33 cycles of 94°C for 30 s, 58°C for 30 s, and extension at 72°C for 1 min and 15 s, followed by a final extension at 72°C for 7 min.
|
View this table: [in a new window] |
TABLE 2. List of primers used in the DNA amplification of cultures
|
The quality and specificity of all the amplification products were assessed by agarose gel electrophoresis. Individual product bands visualized by staining with ethidium bromide were purified using the Wizard SV gel and PCR clean-up system (Promega) according to the manufacturer's instructions. Products eluted from the solid-phase matrix were further cleaned by isopropanol precipitation. After verifying the purity and concentration by agarose gel electrophoresis, individual template samples were sequenced bidirectionally from the appropriate oligonucleotide primer using a commercial service (Northwoods DNA, Inc.). Consensus sequences obtained for this study were deposited in GenBank using the BankIT facility (Table 1).
Alignment and phylogenetic analysis.
The sequences were aligned using ClustalX, version 1.8 (60), and subsequently adjusted manually. The new ITS and LSU sequences from this study (Table 1) were aligned with existing ITS and LSU sequences available in GenBank (Tables 3 and 4). Phylogenetic relationships were inferred based on the simple pairwise differences of nucleotides (p-distance) derived from the alignment, ignoring gaps. The p-distance represents the number of nucleotide differences divided by the total number of nucleotides compared. The neighbor joining (NJ) distance method and maximum parsimony (MP) estimation of phylogenetic relationships were conducted using the Molecular Evolutionary Genetics Analysis (MEGA) program version 3.1 (23). Bootstrapping values of 1,000 replicates were used to assess the consistency of each derived topology. Due to differences in the availability of sequence datasets for the targeted domains, two outgroups were used separately in these analyses: the dinoflagellate Heterocapsa sp. (GenBank accession number AB084100) was used in the ITS analysis, and the dinoflagellate Symbiodinium microadriaticum (GenBank accession number AF060896) was used in the LSU analysis. These outgroups were selected based on their having the largest p-distance with respect to the YTX-producing taxa.
|
View this table: [in a new window] |
TABLE 3. List of sequences used in the ITS region rDNA phylogenetic analysis from GenBank
|
|
View this table: [in a new window] |
TABLE 4. List of sequences used in the LSU region rDNA phylogenetic analysis from GenBank
|
|
|
|---|
![]() View larger version (30K): [in a new window] |
FIG. 1. NJ analysis of the ITS rDNA for all species listed in Tables 1 and 3. Bootstrap values (1,000 replicates) are listed as percentages of 100, and only values greater than 50% are shown. Heterocapsa sp. (GenBank accession number AB084100) was the outgroup. Cultures of Lingulodinium polyedrum and Protoceratium reticulatum in bold are toxic, and cultures of Lingulodinium polyedrum and Gonyaulax spinifera in italics are nontoxic. Cultures of Protoceratium reticulatum for which there are conflicting studies on toxicity are indicated with double asterisks (**) after the name. The toxicity is unknown for all of the YTX-producing species in the regular font.
|
![]() View larger version (36K): [in a new window] |
FIG. 2. Maximum parsimony analysis of the ITS rDNA for all species listed in Tables 1 and 3. Bootstrap values (1,000 replicates) are listed as percentages of 100, and only values greater than 50% are shown. Heterocapsa sp. (GenBank accession number AB084100) was the outgroup. Cultures of Lingulodinium polyedrum and Protoceratium reticulatum in bold are toxic, and cultures of Lingulodinium polyedrum and Gonyaulax spinifera in italics are nontoxic. Cultures of Protoceratium reticulatum for which there are conflicting studies on toxicity are indicated with double asterisks (**) after the name. The toxicity is unknown for all of the YTX-producing species in the regular font.
|
LSU region.
In both the NJ (Fig. 3) and maximum parsimony (Fig. 4) analyses of the LSU domains D1 and D2, the Gonyaulacales order again formed a distinct clade, although with lower support (bootstrap values 80% and 56%, respectively) compared with the ITS analyses.
![]() View larger version (31K): [in a new window] |
FIG. 3. NJ analysis of the LSU rDNA for all species listed in Tables 1 and 4. Bootstrap values (1,000 replicates) are listed as percentages of 100, and only values greater than 50% are shown. Symbiodinium microadriaticum (GenBank accession number AF060896) was the outgroup. Cultures of Lingulodinium polyedrum, Protoceratium reticulatum, and Gonyaulax spinifera in bold are toxic and cultures of Lingulodinium polyedrum and Gonyaulax spinifera in italics are nontoxic. Cultures of Protoceratium reticulatum for which there are conflicting studies on toxicity are indicated with double asterisks (**) after the name. The toxicity is unknown for all of the YTX-producing species in the regular font.
|
![]() View larger version (33K): [in a new window] |
FIG. 4. Maximum parsimony analysis of the LSU rDNA for all species listed in Tables 1 and 4. Bootstrap values (1,000 replicates) are listed as percentages of 100, and only values greater than 50% are shown. Symbiodinium microadriaticum (GenBank accession number AF060896) was the outgroup. Cultures of Lingulodinium polyedrum, Protoceratium reticulatum, and Gonyaulax spinifera in bold are toxic, and cultures of Lingulodinium polyedrum and Gonyaulax spinifera in italics are nontoxic. Cultures of Protoceratium reticulatum for which there are conflicting studies on toxicity are indicated with double asterisks (**) after the name. The toxicity is unknown for all of the YTX-producing species in the regular font.
|
The sequence diversity was considerably higher among isolates designated as G. spinifera. Interestingly, the two toxin-producing strains of G. spinifera, DQ151557 and DQ151558 (45), were genetically identical to each other (0.000 p-distance) but distinct from the non-toxin-producing isolate EU532478 (0.330 p-distance) and GenBank strain AY154960 (0.3 p-distance), toxicity unknown. The toxic strains were less divergent from GenBank strain EF416284 (0.130 p-distance), toxicity unknown. The GenBank submission assigned as G. membranacea (AY154965) was found to be identical to the LSU sequence derived from the G. spinifera strain (CCMP409) used in this study, indicating an incorrect species annotation either in the original submission or in the CCMP culture collection.
|
|
|---|
There are a large number of studies that have used nuclear rRNA genes successfully to evaluate relationships among closely related taxa within genera and particularly within species that are morphologically indistinguishable using the variable domains within the LSU region (27, 29, 32, 54) and potentially more rapidly evolving ITS region (6, 24, 28, 48, 53, 56). While a consensus molecular phylogeny has not been accepted for the Dinophyceae as a whole, there are broad consistencies in the topological positions of some taxonomic orders within this group. For example, previous dinoflagellate phylogeny studies have found the Gonyaulacales order to be a monophyletic and more recently diverged group relative to other described orders within the Dinophyceae (37, 63).
In all phylogenetic analyses in this study, the Gonyaulacales order formed a distinct clade, consistent with previously published studies (7, 13, 37, 49, 59, 63). The molecular sequence analysis for the two loci examined in this study support the validity of the Gonyaulacales and firmly place all known YTX-producing species within this dinoflagellate order. The Prorocentrales order was included in these analyses because Prorocentrum micans was frequently present when YTX was detected in California mussels samples (18; M. Silver, unpublished data), and several isolates from the United Kingdom have been tested for YTX production (58). The results of this study, however, suggest that non-Gonyaulacaloid dinoflagellate species (such as Prorocentrum) are unlikely to produce YTX, and testing such species that fall outside of the Gonyaulacales order may be ineffective for identifying the biological origin of YTX during events.
The L. polyedrum species showed very low intraspecific diversity, and at this level of genetic assessment, there did not appear to be geographically distinct populations but rather a global distribution of ribotypes within this species. Consequently, rRNA operons provide no indirect genetic markers for YTX production, and lack of diversity at this level hints that variation in toxicity may be due to environmental conditions or genomic variability. Toxin production in other algae, such as Alexandrium spp., P. reticulatum, and Pseudo-nitzschia spp., has been shown to be influenced by environmental conditions, particularly nutrients (1, 4, 14, 19, 38, 39, 40), although under common culture conditions, some isolate-specific variation can be maintained (22, 57). Future studies should evaluate if there are specific genes that are "turned on" during toxin production (42), which would explain the differences between toxic and nontoxic isolates of the same, genetically similar species. The lack of genetic diversity indicates that overall quantification may be the most effective methodology in toxic bloom monitoring of L. polyedrum. A quantitative PCR method based on the SSU has been developed for the quantification of L. polyedrum from southern California (35). The results from this study of the LSU and ITS regions of L. polyedrum populations suggest that the SSU quantitative PCR method can be applied to other geographic regions to aid in monitoring L. polyedrum bloom dynamics.
There was a large sequence divergence observed between multiple isolates in the LSU analyses of G. spinifera, which is consistent with other studies of the LSU region of the Gonyaulax genus (11, 20, 45, 46). While a limited number of LSU sequences were used in this study, distinct ribotypes of toxic G. spinifera (DQ151557 and DQ151558) were revealed. In a similar phylogenetic analysis of the LSU region, toxic strains of G. spinifera from New Zealand (DQ151557 and DQ151558) (45) and Italy (46) grouped together and formed a clade distinct from other strains and had a high intraspecific variability compared with GenBank strains (46). The high levels of intraspecific genetic divergence detected in our analysis (30 to 40%) as well as from Riccardi et al. (46) suggest that G. spinifera is undergoing rapid diversification. Considerable morphological variability is observed for G. spinifera (8, 61), and it is possible that cryptic species may exist within the known isolates of G. spinifera. Notwithstanding current taxonomic assignments, the members of the genus Gonyaulax may provide fruitful targets for the development of molecular probes distinguishing isolates of differing toxicity and for the examination of how the expression of toxin production may drive the evolution or diversification of harmful algal species. Additional multilocus and ultrastructural studies using a large number of unique G. spinifera isolates will be needed to thoroughly describe the genetic variability within this species.
Due to the conflicting reports of toxicity for P. reticulatum cultures CCMP404 and CCMP1889, it is difficult to determine if there are genetically distinguishable isolates of this species based on toxicity. However, there are several published studies evaluating the influence of environmental conditions on the production of YTX in P. reticulatum. In isolates of P. reticulatum from Emilia-Romagna, Italy, YTX production increased within the cells when cultures were grown under higher salinity and temperature conditions and under both replete and severe phosphate-limited nutrient conditions (14). YTX released from cells into the medium was found to be higher in cultures grown under nitrogen limitation, lower temperature or lower salinity conditions. Those authors concluded that environmental conditions directly affect toxin production and that decreased temperature and salinity will decrease toxin production but will not terminate toxin production in cells (14). In culture, the stationary growth phase, as well as the addition of the micronutrient selenium but not iron or cobalt, significantly increased YTX production by P. reticulatum (33, 34). These published studies suggest that environmental conditions can influence toxicity and therefore that genetically distinct isolates based on toxicity may not exist. Future research on this species will need to evaluate the genetic diversity of strains for which a complete toxin profile has already been established, such as the isolated cultures from Spain (42).
We hypothesized that the rRNA genes may be useful for the evaluation of genetic variability among toxic and nontoxic isolates. However, the results of this study show that the constrained sequence variability of several rRNA operons do not provide robust markers of YTX toxicity among species for most genera in the Gonyaulacales order. While this study and the results of Riccardi et al. (46) suggest that this application might be possible for G. spinifera, these results are based on small sets of isolates, and therefore, a larger number of G. spinifera isolates needs to be analyzed to characterize the true genetic and toxicity associations within this species. If Gonyaulax is indeed undergoing rapid diversification, this genus may be suitable for genomic tools, but we suggest that future studies focus on both physiological and genomic assays for YTX production in Gonyaulacales. Our results do demonstrate, however, that confirmed YTX production is currently confined to the order Gonyaulacales within the Dinophyceae and that species within this taxonomic order should be given priority for future testing and field collections associated with monitoring for YTX contamination events. Interestingly, a previous study evaluating the origin of paralytic shellfish poisoning (PSP) toxins concluded that PSP toxin-producing species are randomly distributed throughout all dinoflagellate groups and, based on this widespread toxin distribution, speculated that PSP toxins had multiple independent origins in the Gymnodiniales and Gonyaulacales orders (63). While other species within the Gonyaulacales order have not been documented to produce YTX congeners (e.g., Alexandrium catenella), all YTX-producing species identified so far remain within this group. Our results contrast with this hypothesis of multiple origins for PSP toxins, as the occurrence of YTX production within several distinct genera in the Gonyaulacales support the hypothesis that YTX biosynthetic capacity arose early in the divergence of this order and consequently later in the evolutionary history of the Dinophyceae.
We thank Bethany Jenkins and Kendra Hayashi. We thank the Marine Pollution Studies Laboratory at Granite Canyon, University of California, Davis, for the seawater they provided, and Grant Pitcher and Peter Franks for providing algal strains used in this study. We also thank three anonymous reviewers who provided useful commentary on an earlier version of the manuscript.
Published ahead of print on 14 November 2008. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»