Bioluminescence capability and intensity in the dinoflagellate Alexandrium species

Article information

Algae. 2021;36(4):299-314
Publication date (electronic) : 2021 December 15
doi : https://doi.org/10.4490/algae.2021.36.12.6
1School of Earth and Environmental Sciences, College of Natural Sciences, Seoul National University, Seoul 08826, Korea
2Research Institute of Oceanography, Seoul National University, Seoul 08826, Korea
3Department of Marine Biotechnology, College of Ocean Sciences, Kunsan National University, Kunsan 54150, Korea
4Department of Marine Biotechnology, Anyang University, Incheon 23038, Korea
*Corresponding Author: E-mail: hjjeong@snu.ac.kr, Tel: +82-2-880-6746, Fax: +82-2-874-9695
Received 2021 October 30; Accepted 2021 December 6.

Abstract

Some species in the dinoflagellate genus Alexandrium are bioluminescent. Of the 33 formally described Alexandrium species, the bioluminescence capability of only nine species have been tested, and eight have been reported to be bioluminescent. The present study investigated the bioluminescence capability of seven Alexandrium species that had not been tested. Alexandrium mediterraneum, A. pohangense, and A. tamutum were bioluminescent, but A. andersonii, A. hiranoi, A. insuetum, and A. pseudogonyaulax were not. We also measured the bioluminescent intensity of A. affine, A. fraterculus, A. mediterraneum, A. ostenfeldii, A. pacificum, A. pohangense, A. tamarense, and A. tamutum. The mean 200-second-integrated bioluminescence intensity per cell ranged from 0.02 to 32.2 × 104 relative luminescence unit per cell (RLU cell−1), and the mean maximum bioluminescence intensity per cell per second (BLMax) ranged from 0.01 to 10.3 × 104 RLU cell−1 s−1. BLMax was significantly correlated with the maximum growth rates of Alexandrium species, except for A. tamarense. A phylogenetic tree based on large subunit ribosomal DNA (LSU rDNA) showed that the bioluminescent species A. affine, A. catenella, A. fraterculus, A. mediterraneum, A. pacificum, and A. tamarense formed a large clade. However, the toxicity or mixotrophic capability of these species was split. Thus, their bioluminescence capability in this clade was more consistent than their toxicity or mixotrophic capability. Phylogenetic trees based on LSU rDNA and the luciferase gene of Alexandrium were consistent except for A. pohangense. The results of the present study can provide a basis for understanding the interspecific diversity in bioluminescence of Alexandrium.

INTRODUCTION

Bioluminescence, the production of light by living organisms, is a widespread phenomenon in nature (Haddock et al. 2010, Valiadi and Iglesias-Rodriguez 2013). Light is produced when luciferase (enzyme) binding with luciferin causes an oxidation reaction of luciferin (substrate) (McCapra 1976, Wilson and Hastings 2013). A majority of the bioluminescent organisms reside in the ocean (Shimomura 2006, Widder 2010).

Dinoflagellates are ubiquitous protists and one of the most common bioluminescent producers (Haddock et al. 2010, Le Tortorec et al. 2016, Kang et al. 2019, Cusick and Widder 2020, Jang and Jeong 2020, Lee et al. 2020). In dinoflagellates, bioluminescence is a defense mechanism against predation and helps attract secondary predators (a function sometimes called a burglar alarm) (Esaias and Curl 1972, Haddock et al. 2010, Valiadi and Iglesias-Rodriguez 2013, Lindström et al. 2017). Several species of bioluminescent dinoflagellates form red tides, dense algal blooms (Valiadi et al. 2012, Cusick and Widder 2014, Jeong et al. 2015); due to their brightness, fish, mammals, and ships passing red-tide patches can be detected at night (Cram and Schülein 1974, Rohr et al. 1998, Miller et al. 2005). Thus, dinoflagellate bioluminescence has been studied in terms of ecophysiology and military uses. However, the bioluminescent capability of dinoflagellates has been tested in only a small portion of 3,450 formally described dinoflagellate species (e.g., Sweeney 1963, Kelly 1968, Swift et al. 1995, Latz and Jeong 1996, Jeong et al. 2021). Among the tested dinoflagellate species, only 70 (<3% of 3,450 formally described dinoflagellates) have been reported to be bioluminescent (Marcinko et al. 2013, Cusick and Widder 2014, 2020). Furthermore, the bioluminescence intensity of only 14 species (<1% of 3,450 formally described dinoflagellates) has been reported, including that of Alexandrium tamarense, A. catenella, Ceratocorys horrida, Lingulodinium polyedra, Noctiluca scintillans, Pyrodinium bahamense, Protoperidinium spp., Pyrocystis spp., and Tripos spp. (Seliger et al. 1969, Esaias and Curl 1972, Esaias et al. 1973, Swift et al. 1973, White 1979, Widder and Case 1981, Lapota et al. 1989, Latz and Lee 1995, Sullivan and Swift 1995, Latz et al. 2004, Cussatlegras and Le Gal 2007, Eckert 2015). Thus, the bioluminescent capability and intensity of more dinoflagellate species should be explored.

The genus Alexandrium is widely distributed across coastal waters and is a major organism causing red tides or harmful algal blooms (Anderson et al. 2012, Jeong et al. 2017). However, among the 33 known Alexandrium species (Guiry and Guiry 2021), only nine species have been tested for the bioluminescence capability (Esaias and Curl 1972, Esaias et al. 1973, Liu et al. 2004, Baker et al. 2008, Latz et al. 2008, Kremp et al. 2009, Martínez et al. 2016). Furthermore, bioluminescence intensity has been reported in only three Alexandrium species (Esaias and Curl 1972, White 1979), and it widely ranged from 2 × 104 to 7 × 107 photons per cell. As the bioluminescence intensity of Alexandrium spp. may be species-dependent, it should be measured in more Alexandrium species. In addition, bioluminescence intensity should be assessed using the same technique across species, in order to reveal the potential factors that affect the bioluminescence intensity.

In the present study, the bioluminescence capability of seven Alexandrium species (A. andersonii, A. hiranoi, A. insuetum, A. mediterraneum, A. pohangense, A. pseudogonyaulax, and A. tamutum) was investigated. The bioluminescence intensity of eight Alexandrium species (A. affine, A. fraterculus, A. mediterraneum, A. ostenfeldii, A. pacificum, A. pohangense, A. tamarense, and A. tamutum) was also explored. Moreover, the effects of cell size, mixotrophic ability, toxicity, maximum swimming speed, and maximum growth rate on bioluminescence intensity were investigated in eight Alexandrium species. In addition, the sequence of the luciferase gene (lcf) was analyzed in A. fraterculus, A. mediterraneum, A. pohangense, and A. tamutum, and the phylogenetic trees based on the lcf and large subunit ribosomal DNA (LSU rDNA) of Alexandrium spp. were compared. The results of the present study can provide a basis for understanding the interspecific diversity in the bioluminescence capability and intensity of Alexandrium.

MATERIALS AND METHODS

Cultures of experimental organisms

Cells of A. affine, A. fraterculus, A. mediterraneum, A. pacificum, A. pohangense, and A. tamutum were isolated from Korean coastal waters and then established using two serial single-cell isolations (Table 1). Cultures of the other six Alexandrium species (A. andersonii, A. minutum, A. insuetum, A. ostenfeldii, A. tamarense, and A. hiranoi) used in the present study were obtained from the National Center for Marine Algae and Microbiota (NCMA; USA), Roscoff Culture Collection (RCC; France), and the Microbial Culture Collection at the National Institute for Environmental Studies (NIES Collection, Japan). All experimental cultures were maintained at 20°C in a chamber under cool white fluorescent lights (20 μmol photons m−2 s−1) on a 14: 10 h light / dark cycle. The culture of A. pohangense was grown mixotrophically by feeding on Margalefidinium polykrikoides. Other cultures were grown autotrophically in F/2 (Guillard and Ryther 1962) or L1 (Guillard and Hargraves 1993) seawater medium without silicate (Table 1).

Information on the isolation and maintenance of the experimental organisms and also Alexandrium species in which bioluminescence capability was previously tested or newly tested

Bioluminescence capability

The bioluminescence capability of each Alexandrium species was tested when the cultures reached their late exponential to early stationary phases. Cultures of the seven Alexandrium species were mechanically and chemically stimulated (Table 1). For mechanical stimulation, dense Alexandrium cultures were placed in the dark for 3 h and then shaken by hand in a dark room. During mechanical stimulation, bioluminescence images were captured using a Canon EOS R with a Canon RF 50 mm F1.2L lens (Canon, Tokyo, Japan) at aperture F1.2, shutter speed 30 s, ISO4000.

If mechanically stimulated bioluminescence was not visible using the digital camera, chemical stimulation was conducted. For the chemical stimulation, four 200 μL aliquots of dense cultures of each Alexandrium species were placed in four wells of Corning 96-well white plates (Corning Life Sciences, Amsterdam, Netherlands). Four culture media without Alexandrium cells were also placed in four wells as controls, and analysis of these wells confirmed no bioluminescence. The bioluminescence analysis of each Alexandrium species was performed 3 h after the beginning of the dark phase, when the cells are known to be fully dark-adapted and produce highly stimulated bioluminescence (Biggley et al. 1969, Krasnow et al. 1980). Chemically stimulated bioluminescence was measured after the addition of 50 μL of 1 M acetic acid to the well (Hastings and Sweeney 1957, Fogel and Hastings 1972, Sweeney 1986). Acidification by the addition of acetic acid, which reduces the intracellular pH to <6–7, is known to bind the luciferase and luciferin of bioluminescent species (Fogel and Hastings 1972, Sweeney 1986, Von Dassow and Latz 2002). The bioluminescence response was assayed for 200 s using a GloMax Navigator microplate luminometer (Promega, Madison, WI, USA).

Measurement of bioluminescence intensity

To measure the bioluminescence intensity of eight Alexandrium species, triplicate 1 mL aliquots were removed from a culture of each Alexandrium species and then fixed with 5% acid Lugol’s solution, and Alexandrium cells were enumerated under a light microscope to determine the density of each culture. Cells were added to a 50 mL culture flask to reach the target cell density of 1,000 cells mL−1, using media in which the cells were grown. A 200 μL aliquot of diluted culture was placed in a well in a Corning 96-well white plate (Corning Life Sciences), with four replicates. The analysis of bioluminescence was conducted using the method described above.

Bioluminescence data analysis

To obtain the bioluminescence intensity per cell per second of each Alexandrium species, the mean value obtained from the four control wells (containing seawater without bioluminescent Alexandrium cells) was subtracted from the bioluminescence intensity value of the cells in one well, and then the value was divided by the total number of cells in the well. Using the data on bioluminescence intensity per cell per second of each Alexandrium species, several parameters, such as the integrated bioluminescence intensity per cell for 200 s and the maximum bioluminescence intensity per cell per second, were calculated as described by Latz and Lee (1995). The mean 200-second-integrated bioluminescence intensity per cell (BLInt200s) was calculated by averaging the sums of the bioluminescence intensity per cell per 200 s in four wells. Furthermore, the mean maximum bioluminescence intensity per cell per second (BLMax) was calculated by averaging the highest bioluminescence intensities per cell per second in the four wells. In the present study, relative luminescence unit per cell (RLU cell−1) and RLU per cell per second (RLU cell−1 s−1) were used as the units of BLInt200s and BLMax, respectively (Lindström et al. 2017).

Sequencing large subunit ribosomal DNA and lcf

To obtain the LSU rDNA and lcf sequences, a 10 mL aliquot of each dense culture was harvested at 3,667 g (4,000 rpm) for 10 min in a 15 mL conical tube. The gDNA was extracted from the pellet using a Wizard SV Genomic DNA Purification System (Promega).

The amplification reaction mixtures (50 μL in total volume) were as follows: 5 μL of 10× F-StarTaq buffer, 1 μL of 10 mM of dNTP mix, 0.25 μL of 5 U μL−1 BioFACT F-Star Taq DNA polymerase (BioFACT Co., Ltd., Daejeon, Korea), 2 μL of 10 μM of each primer, and 1 μL of the extracted gDNA. The primers used to amplify regions of LSU rDNA were D1R (5-ACCCGCTGAATTTAAGCATA-3) (Scholin et al. 1994), 28–1483R (5-GCTACTACCACCAAGATCTGC-3) (Daugbjerg et al. 2000), and LSUB (5-ACGAACGATTTGCACGTCAG-3) (Litaker et al. 2003). The gDNA was amplified in an AllInOneCycler (Bioneer, Daejeon, Korea) under the following conditions: 2 min at 95°C for initial denaturation, 38 cycles of 40 s at 95°C, 30 s at the annealing temperature, and 1 min at 72°C, with a final extension of 5 min at 72°C. The universal lcf primers for amplifying the lcf of dinoflagellates were used: LcfUniCHF3 (5′-TCCAGGTTGCACGCTTCGA-3′) (Baker et al. 2008) and LcfUniCHR4 (5′-GGGTCTTGTCGCCGTACTCAAA-3′) (Baker et al. 2008). These primer sets (LcfCHF3 and LcfCHR4) targeted the lcf N-terminal region, which is more diverse than the central regions (Baker et al. 2008). Polymerase chain reactions (PCRs) were performed under the following conditions: 5 min at 95°C for initial denaturation, 35 cycles of 45 s at 95°C, 30 s at 62°C and 45 s at 68°C, and a final extension at 68°C for 5 min.

The PCR products were purified using an AccuPrep DNA Purification Kit (Bioneer), and sequencing was performed using an ABI 3730XL DNA Analyzer (Applied Biosystems, Foster City, CA, USA).

Phylogenetic analysis

To explore the phylogenetic relationships of the bioluminescent Alexandrium species, revealed in the present study or in previous studies, the sequences of the LSU rDNA region of the Alexandrium species were aligned using the software MEGA v4 (Tamura et al. 2007), including sequences obtained from GenBank. Maximum likelihood (ML) analysis of the region was conducted using the program RAxML 7.0.3, with the general tIme reversible plus gamma (GTR + GAMMA) model (Stamatakis 2006). Furthermore, 200 independent tree inferences were used to identify the best tree. The ML bootstrap values were determined using 1,000 replicates.

To investigate any differences in the phylogenetic trees based on the LSU rDNA region and lcf of the individual Alexandrium species, the lcf sequences of the six Alexandrium species obtained in the present study and the five Alexandrium species available in GenBank were aligned. Protoceratium reticulatum was added because it formed a clade with A. pohangense, and Gonyaulax spinifera was added because it was positioned as a base of an Alexandrium clade. Lingulodinium polyedra was added as an outgroup. ML analysis of the region was conducted as described previously.

Bayesian analysis was conducted using MrBayes v.3.1 (Huelsenbeck and Ronquist 2001, Ronquist and Huelsenbeck 2003) with the default GTR + G + I model to determine the best available model for the data. For all sequence regions, four independent Markov chain Monte Carlo runs were performed simultaneously until the average standard deviation of split frequencies reached <0.01. The trees were sampled every 1,000 generations. To ensure likelihood convergence, the first 1,000 trees were discarded as burn-ins. Five million generations of LSU rDNA and three million generations of the lcf of Alexandrium species were run in MrBayes, and at least two million generations were saved for each tree reconstruction. The parameters were as follows for assumed nucleotide frequencies of LSU rDNA: substitution rate matrix with A–C substitution = 0.106, A–G substitution = 0.236, A–T substitution = 0.103, C–G substitution = 0.043, C–T substitution = 0.416, G–T substitution = 0.095; the proportion of sites assumed to be invariable = 0.051; and the rate for variable sites assumed to follow a gamma distribution with shape parameter = 1.075. For the assumed nucleotide frequencies of the lcf, parameters were: substitution rate matrix with A–C substitution = 0.090, A–G substitution = 0.255, A–T substitution = 0.071, C–G substitution = 0.074, C–T substitution = 0.424, and G–T substitution = 0.086; the proportion of sites assumed to be invariable = 0.391; and the rate for variable sites assumed to follow a gamma distribution with shape parameter = 54.594. Moreover, the assumed nucleotide frequencies of LSU rDNA comprised the sequences of the Alexandrium species belonging to the lcf phylogeny, substitution rate matrix with A–C substitution = 0.120, A–G substitution = 0.228, A–T substitution = 0.089, C–G substitution = 0.066, C–T substitution = 0.395, and G–T substitution = 0.101; the proportion of sites assumed to be invariable = 0.041; and the rate for variable sites assumed to follow a gamma distribution with shape parameter = 1.353.

Statistical analysis

To test whether the BLInt200s or BLMax of eight Alexandrium species differ from one another, univariate analyses were performed. Normality and homogeneity of variance were tested using Shapiro-Wilk’s W and Levene’s tests, respectively. A parametric one-way analysis of variance (ANOVA) was performed. If the data did not satisfy the homogeneity assumption, Welch’s one-way ANOVA and the Games-Howell post-hoc test were performed (Welch 1947, Games and Howell 1976). When the data did not meet the normality assumption, a nonparametric Kruskal-Wallis test and Mann-Whitney U test with Bonferroni correction (p < 0.05) were conducted (Mann and Whitney 1947, Kruskal and Wallis 1952, Dunn 1961). Spearman’s rank correlation coefficient was used to examine relationships between the variables equivalent spherical diameter (ESD) of bioluminescence species, maximum swimming speed, maximum growth rate, BLInt200s, and BLMax. Statistical analyses were performed using SPSS ver. 25.0 (IBM-SPSS Inc., Armonk, NY, USA).

RESULTS

Bioluminescence capability of Alexandrium species

Among the seven Alexandrium species tested, A. mediterraneum AMYS1807, A. pohangense APPH1409, and A. tamutum ATSH1609 had bioluminescent capability, whereas A. andersonii CCMP2222, A. hiranoi NIES-3611, A. insuetum CCMP2082, and A. pseudogonyaulax KSUDinoB4 did not (Table 2, Fig. 1).

Mean 200-second-integrated bioluminescence intensity per cell (BLInt200s, ×104 RLU cell−1 ± standard error) and mean maximum bioluminescence intensity per cell per second (BLMax, ×104 RLU cell−1 s−1 ± standard error) of the 13 Alexandrium species tested in the present study

Fig. 1

Photographs of bioluminescence showing the bioluminescent capability of the Alexandrium species tested in the present study. A 250 mL culture flask containing a dense culture of each species in the late exponential stage was shaken by hand. (A) Alexandrium mediterraneum AMYS1807 (30,000 cells mL−1). (B) Alexandrium pohangense APPH1409 (2,000 cells mL−1). (C) Alexandrium tamutum ATSH1609 (10,000 cells mL−1). (D) Alexandrium andersonii CCMP2222 (7,000 cells mL−1). (E) Alexandrium hiranoi NIES-3611 (1,000 cells mL−1). (F) Alexandrium pseudogonyaulax KSUDi-noB4 (1,000 cells mL−1). (G) Alexandrium insuetum CCMP2082 (30,000 cells mL−1).

Position of bioluminescent Alexandrium species in a phylogenetic tree based on the LSU rDNA sequences

In a phylogenetic tree based on the LSU rDNA sequences, all 10 bioluminescent Alexandrium species were positioned in all four large clades (Fig. 2). Of the 10 species, six (A. affine, A. catenella, A. fraterculus, A. mediterraneum, A. pacificum, and A. tamarense) belonged to a large clade (called bioluminescence clade I, BLClade I), whereas two species (A. ostenfeldii and A. tamutum) belonged to another large clade (BLClade II). However, A. pohangense (BLClade III) and A. monilatum (BLClade IV) belonged to other large clades.

Fig. 2

Consensus Bayesian tree of nuclear large subunit ribosomal DNA region (based on 890 bp aligned position), using the GTR + G + I model. The bioluminescent dinoflagellates Lingulodinium polyedra, Protoceratium reticulatum, and Gonyaulax spinifera were used as the outgroups. The numbers above the branches indicate the Bayesian posterior probability (≥0.5, left) and maximum likelihood bootstrap values (≥50%, right), and solid black circles indicate full support (1.0 and 100%). Bioluminescent strains are indicated by a blue closed star, whereas nonbioluminescent strains are indicated by a purple open star. The strains tested in the present study are written in bold.

Phylogenetic tree based on the lcf sequences

In phylogenetic trees based on the lcf and LSU rDNA sequences, the positions of nine bioluminescent Alexandrium species were the same, with the exception of A. pohangense (Fig. 3). In the phylogenetic tree based on the lcf sequences, A. pohangense formed a separate clade with Protoceratium reticulatum (Fig. 3A); however, in the phylogenetic tree based on the LSU rDNA sequences, A. pohangense was positioned inside a large Alexandrium clade (Fig. 3B).

Fig. 3

Consensus Bayesian tree, based on 498 bp aligned position of nuclear luciferase gene (lcf) (A) and based on 890 bp aligned position of the large subunit ribosomal DNA (B), using the GTR + G + I model and Lingulodinium polyedra as the outgroup. The numbers above the branches indicate the Bayesian posterior probability (≥0.5, left) and maximum likelihood bootstrap values (≥50%, right), and solid black circles indicate full support (1.0 and 100%). The sequences obtained in the present study are in boldface.

Bioluminescence intensity of eight Alexandrium species

The BLInt200s of eight Alexandrium species ranged from 0.02 × 104 to 32.2 × 104 RLU cell−1 (Table 2, Fig. 4A); that of A. ostenfeldii was the highest, whereas that of A. tamarense was the lowest. The BLInt200s values of the eight Alexandrium species were significantly different (Welch’s ANOVA, F7, 9.06 = 67.15, p < 0.01) and were divided into three subsets (Games-Howell post-hoc test, p < 0.05) (Fig. 4A). The BLMax of the eight Alexandrium species ranged from 0.01 × 104 to 10.3 × 104 RLU cell−1 s−1 (Table 2, Fig. 4B). The BLMax values of the eight Alexandrium species were significantly different (Kruskal-Wallis test, H7 = 28.80, p < 0.01), and the Mann-Whitney U test with Bonferroni correction (p < 0.05) revealed that the BLMax values of the Alexandrium species were divided into three subsets (Fig. 4B).

Fig. 4

Bioluminescence intensity of each Alexandrium species. (A) The mean 200-second-integrated bioluminescence intensity per cell (BLInt200s). Significantly different subsets of the Alexandrium species based on Games-Howell post-hoc tests (p < 0.05) of Welch’s ANOVA: (a) Alexandrium tamarense (Ata) and A. mediterraneum (Ame); (b) A. tamutum (Atm) and A. pacificum (Apa); (abc) A. affine (Aaf) and A. fraterculus (Afr); (c) A. ostenfeldii (Aos) and A. pohangense (Apo). (B) The mean maximum bioluminescence intensity per cell per second (BLMax). Significantly different subsets of Alexandrium species based on the Bonferroni post-hoc test (p < 0.05) of Kruskal-Wallis test: (a’) Ata; (a’b’) Ame and Atm; (a’b’c’) Apa, Aaf, and Apo; (b’c’) Afr; (c’) Aos. RLU, relative luminescence unit. Values represent the treatment means ± standard error (n = 4).

DISCUSSION

Bioluminescence capability

The present study explored the bioluminescence capability of seven Alexandrium species (A. andersonii, A. hiranoi, A. insuetum, A. pseudogonyaulax, A. mediterraneum, A. pohangense, and A. tamutum). Previously, the bioluminescence capability of nine Alexandrium species was explored. Thus, the bioluminescence capability of a total of 16 Alexandrium species, almost half the 33 formally described Alexandrium species, has now been explored. It is worthwhile to investigate the bioluminescence capability of the remaining 17 Alexandrium species.

The present study is the first to examine bioluminescence in A. mediterraneum, A. pohangense, and A. tamutum. Previously, eight Alexandrium species were reported to be bioluminescent; thus, a total of 11 Alexandrium species, 33% of the 33 formally described Alexandrium species, have been revealed to be bioluminescent. Before the present study was conducted, 90% of the tested Alexandrium species (eight of nine tested species) were known to be bioluminescent. However, the addition of our data shows that only 69% of tested Alexandrium species were bioluminescent. Thus, it is incorrect to assert that as a rule, Alexandrium species are bioluminescent; in fact, bioluminescence capability of these organisms is species-dependent.

Position of bioluminescent Alexandrium species in a phylogenetic tree based on the LSU rDNA sequences

In a phylogenetic tree based on the LSU rDNA sequences, six of 10 bioluminescent Alexandrium species belonged to a large clade (BLClade I), although the 10 bioluminescent Alexandrium species spanned four large clades (BLClade I–IV). In BLClade I, all Alexandrium species are known to be bioluminescent, except for A. australiense, A. cohorticula, and A. tamiyavanichii, whose bioluminescence capability has not yet been tested. Thus, it is worthwhile to explore the bioluminescence capability of these three species to test the hypothesis that all Alexandrium species in BLClade I are bioluminescent.

In BLClade I, A. affine AAV1, A. catenella CCMP1719, A. pacificum CCMP1598, and A. tamarense UW2C are known to be toxic, whereas A. fraterculus AF0307MIE01 and A. mediterraneum SZN01 are known to be nontoxic (Higman et al. 2001, Nguyen-Ngoc 2004, Orr et al. 2011, John et al. 2014, Eckford-Soper et al. 2016, Subong et al. 2017, Blossom et al. 2019). Both bioluminescence and toxicity are tools for anti-predation (Esaias and Curl 1972, Colin and Dam 2003, Bergkvist et al. 2008, Wohlrab et al. 2010, Lindström et al. 2017). Although all six Alexandrium species in BLClade I were bioluminescent, the toxicity status of these six species was split (Table 3, Fig. 5). Thus, bioluminescence is likely to be a primary defense tool and toxicity as a supplementary tool in this clade. Furthermore, A. catenella, A. pacificum CCMP3434, and A. tamarense are known to be mixotrophic, whereas A. affine CCMP112, A. fraterculus AFYS1309, and A. mediterraneum CCMP3433 are not mixotrophic (Table 3, Fig. 5) (Jeong et al. 2005, Lee et al. 2016, Lim et al. 2019). Thus, the acquisition of a mixotrophic capability in the six Alexandrium species in BLClade I may have occurred later than that of bioluminescence capability. However, further studies are needed to confirm the timing of the acquisition of a mixotrophic and bioluminescence capability in Alexandrium species.

Comparison of bioluminescence capability (BC), mixotrophy, and toxicity of Alexandrium species

Fig. 5

Comparison of the presence (closed symbols) or absence (open symbols) of bioluminescent capacity (BL), toxicity (Toxi), and mixotrophic ability (Mixo) in the Consensus Bayesian tree of the nuclear large subunit ribosomal DNA region in Fig. 2. Asterisks indicate that some strains are toxic, whereas others are nontoxic (Prakash and Taylor 1966, Murakami et al. 1988, 1998, Jacobson and Anderson 1996, Satake et al. 1997, Draisci et al. 1999, Ciminiello et al. 2000, Higman et al. 2001, MacKenzie et al. 2004, Montresor et al. 2004, Nguyen-Ngoc 2004, Jeong et al. 2005, Lilly et al. 2005, Yoo et al. 2009, Orr et al. 2011, 2013, Blossom et al. 2012, 2017, John et al. 2014, Lim et al. 2015, 2019, Eckford-Soper et al. 2016, Lee et al. 2016, Subong et al. 2017, Murray et al. 2020).

Phylogenetic tree based on lcf sequences

The positions of nine bioluminescent Alexandrium species in the phylogenetic trees based on the lcf sequences were almost the same as those in the phylogenetic trees based on the LSU rDNA. This congruence is also shown in internal transcribed spacer and small subunit rDNA phylogenies (Baker et al. 2008, Valiadi et al. 2012, Le Tortorec et al. 2016). Thus, the lcf and ribosomal DNA may have coevolved. In the phylogenetic trees based on LSU rDNA sequences, A. pohangense was included in BLClade III. However, in the phylogenetic trees based on lcf genes, A. pohangense formed a clade with P. reticulatum, separate from the other Alexandrium species, unlike in the phylogenetic trees based on LSU rDNA sequences. The sequences in the lcf N-terminal region of A. pohangense differed from those of P. reticulatum by 41 bp (9.8%) (GenBank accession No. AY766386), but they differed from the other Alexandrium species by 43–81 bp (10.3–19.3%) (GenBank accession numbers listed in Fig. 3). Similarly, Valiadi et al. (2012) reported that the lcf in the catalytic domains of A. monilatum (included in BLClade IV of the phylogenetic tree based on LSU rDNA) was more similar to that of L. polyedra than that of Alexandrium species. Thus, it may not be surprising that A. pohangense forms a clade with P. reticulatum, rather than other Alexandrium species.

Bioluminescence intensity of eight Alexandrium species

The results of the present study clearly showed that there were interspecific variations in the bioluminescence intensity of Alexandrium species. However, none of the BLInt200s or BLMax values of the Alexandrium species investigated in the present study were significantly correlated with the ESD (Fig. 6). Thus, other factors may affect the BLInt200s or BLMax of Alexandrium species. In addition, none of the BLInt200s or BLMax values of the Alexandrium species investigated in the present study were significantly correlated with the maximum growth rate (Fig. 7). However, when the maximum growth rate of A. tamarense was excluded, the BLMax of Alexandrium species was negatively correlated with the maximum growth rate (Spearman’s rank correlation coefficient; rs = −1.00, p < 0.01). In general, the energy gained by a bioluminescent organism is used for reproduction or growth, respiration (swimming and metabolism), excretion, and bioluminescence. Thus, the bioluminescence intensity of bioluminescent organisms may be low if the growth rate is high.

Fig. 6

Mean 200-second-integrated bioluminescence intensity per cell (BLInt200s) (A) and mean maximum bioluminescence intensity per cell per second (BLMax) (B) of bioluminescent Alexandrium species as a function of cell size (equivalent spherical diameter [ESD], μm), as described in Table 2. The values represent the treatment means ± standard error (n = 4). The Spearman correlation coefficient and p-value were described in the figures. Aaf, A. affine; Afr, A. fraterculus; Ame, A. mediterraneum; Aos, A. ostenfeldii; Apa, A. pacificum; Apo, A. pohangense; Ata, A. tamarense; Atm, A. tamutum.

Fig. 7

Mean 200-second-integrated bioluminescence intensity per cell (BLInt200s) (A) and the mean maximum bioluminescence intensity per cell per second (BLMax) (B) of the bioluminescent Alexandrium species as a function of maximum growth rate (MGR; d−1). Values represent the treatment means ± standard error (n = 4). a: the Spearman correlation coefficient and p-value of all species. b: the Spearman correlation coefficient and p-value of all species except for A. tama-rense (red circle). Aaf, A. affine; Afr, A. fraterculus; Aos, A. ostenfeldii; Apa, A. pacificum; Apo, A. pohangense; Ata, A. tamarense. MGR data were obtained from Gu (2011), Laabir et al. (2011), Lim et al. (2015), Bill et al. (2016), and Lee et al. (2019).

None of the BLInt200s or BLMax values of the Alexandrium species investigated in the present study were significantly correlated with the maximum swimming speed (Fig. 8) (Spearman’s rank correlation coefficient; rs = −0.01 and −0.12, p > 0.05). However, the BLInt200s and BLMax values of A. ostenfeldii were considerably greater than those of the other Alexandrium species. Therefore, A. ostenfeldii may spend less energy on swimming and may demonstrate notably higher BLInt200s or BLMax values than other Alexandrium species.

Fig. 8

Mean 200-second-integrated bioluminescence intensity per cell (BLInt200s) (A) and the mean maximum bioluminescence intensity per cell per second (BLMax) (B) of the bioluminescent Alexandrium species as a function of maximum swimming speed (MSS; μm s−1). Values represent the treatment means ± standard error (n = 4). The Spearman correlation coefficient and p-value were described in the figures. Aaf, A. affine; Afr, A. fraterculus; Ame, A. mediterraneum; Aos, A. ostenfeldii; Apa, A. pacificum; Apo, A. pohangense; Ata, A. tamarense; Atm, A. tamutum. MSS data were obtained from Lewis et al. (2006), Lim et al. (2015, 2018), Lee et al. (2016), and Kang et al. (2018).

The present study enhanced our knowledge of bioluminescence in Alexandrium. It is the first to report bioluminescence capability in A. mediterraneum, A. pohangense, and A. tamutum and to report lack of bioluminescence capability in A. andersonii, A. hiranoi, A. insuetum, and A. pseudogonyaulax. Furthermore, this is the first study to report the bioluminescence intensity of more than three species in one genus and to explore factors affecting differences in bioluminescence intensity. Moreover, our data revealed that phylogenetic trees based on the LSU rDNA and lcf sequences of Alexandrium species were consistent except for A. pohangense. The bioluminescence capabilities of 17 Alexandrium species have not yet been explored and provide a rich opportunity for further studies of bioluminescence in Alexandrium species.

ACKNOWLEDGEMENTS

This research was supported by the National Research Foundation (NRF) funded by the Ministry of Science and ICT (NRF-2020M3F6A1110582; NRF-2021M3I6A1091272; 2021R1A2C1093379) and by the Useful Dinoflagellate program of Korea Institute of Marine Science and Technology Promotion (KIMST) funded by the Ministry of Oceans and Fisheries (MOF) award to HJJ.

Notes

CONFLICTS OF INTEREST

The authors declare that they have no potential conflicts of interest.

References

Anderson DM, Alpermann TJ, Cembella AD, Collos Y, Masseret E, Montresor M. 2012;The globally distributed genus Alexandrium: multifaceted roles in marine ecosystems and impacts on human health. Harmful Algae 14:10–35.
Baker A, Robbins I, Moline MA, Iglesias-Rodríguez MD. 2008;Oligonucleotide primers for the detection of bioluminescent dinoflagellates reveal novel luciferase sequences and information on the molecular evolution of this gene. J Phycol 44:419–428.
Bergkvist J, Selander E, Pavia H. 2008;Induction of toxin production in dinoflagellates: the grazer makes a difference. Oecologia 156:147–154.
Biggley WH, Swift E, Buchanan RJ, Seliger HH. 1969;Stimulable and spontaneous bioluminescence in the marine dinoflagellates, Pyrodinium bahumense, Gonyaulax polyedra, and Pyrocystis lunula . J Gen Physiol 54:96–122.
Bill BD, Moore SK, Hay LR, Anderson DM, Trainer VL. 2016;Effects of temperature and salinity on the growth of Alexandrium (Dinophyceae) isolates from the Salish Sea. J Phycol 52:230–238.
Blossom HE, Bædkel TD, Tillmann U, Hansen PJ. 2017;A search for mixotrophy and mucus trap production in Alexandrium spp. and the dynamics of mucus trap formation in Alexandrium pseudogonyaulax . Harmful Algae 64:51–62.
Blossom HE, Daugbjerg N, Hansen PJ. 2012;Toxic mucus traps: a novel mechanism that mediates prey uptake in the mixotrophic dinoflagellate Alexandrium pseudogonyaulax . Harmful Algae 17:40–53.
Blossom HE, Markussen B, Daugbjerg N, Krock B, Norlin A, Hansen PJ. 2019;The cost of toxicity in microalgae: direct evidence from the dinoflagellate Alexandrium . Front Microbiol 10:1065.
Ciminiello P, Fattorusso E, Forino M, Montresor M. 2000;Saxitoxin and neosaxitoxin as toxic principles of Alexandrium andersoni (Dinophyceae) from the Gulf of Naples, Italy. Toxicon 38:1871–1877.
Colin SP, Dam HG. 2003;Effects of the toxic dinoflagellate Alexandrium fundyense on the copepod Acartia hudsonica: a test of the mechanisms that reduce ingestion rates. Mar Ecol Prog Ser 248:55–65.
Cram DL, Schülein FH. 1974;Observations on surface-shoaling Cape hake off South West Africa. ICES J Mar Sci 35:272–275.
Cusick KD, Widder EA. 2014;Intensity differences in bioluminescent dinoflagellates impact foraging efficiency in a nocturnal predator. Bull Mar Sci 90:797–811.
Cusick KD, Widder EA. 2020;Bioluminescence and toxicity as driving factors in harmful algal blooms: ecological functions and genetic variability. Harmful Algae 98:101850.
Cussatlegras A-S, Le Gal P. 2007;Variability in the bioluminescence response of the dinoflagellate Pyrocystis lunula . J Exp Mar Biol Ecol 343:74–81.
Daugbjerg N, Hansen G, Larsen J, Moestrup Ø. 2000;Phylogeny of some of the major genera of dinoflagellates based on ultrastructure and partial LSU rDNA sequence data, including the erection of three new genera of unarmoured dinoflagellates. Phycologia 39:302–317.
Draisci R, Ferretti E, Palleschi L, Marchiafava C, Poletti R, Milandri A, Ceredi A, Pompei M. 1999;High levels of yessotoxin in mussels and presence of yessotoxin and homoyessotoxin in dinoflagellates of the Adriatic Sea. Toxicon 37:1187–1193.
Dunn OJ. 1961;Multiple comparisons among means. J Am Stat Assoc 56:52–64.
Eckert R. 2015. Excitation and luminescence in Noctiluca miliaris . In : Johnson FH, Haneda Y, eds. Bioluminescence in Progress Princeton University Press. Princeton, NJ: p. 269–300.
Eckford-Soper LK, Bresnan E, Lacaze J-P, Green DH, Davidson K. 2016;The competitive dynamics of toxic Alexandrium fundyense and non-toxic Alexandrium tamarense: the role of temperature. Harmful Algae 53:135–144.
Esaias WE, Curl HC Jr. 1972;Effect of dinoflagellate bioluminescence on copepod ingestion rates. Limnol Oceanogr 17:901–906.
Esaias WE, Curl HC Jr, Seliger HH. 1973;Action spectrum for a low intensity, rapid photoinhibition of mechanically stimulable bioluminescence in the marine dinoflagellates Gonyaulax catenella, G. acatenella, and G tamarensis . J Cell Physiol 82:363–372.
Fogel M, Hastings JW. 1972;Bioluminescence: mechanism and mode of control of scintillon activity. Proc Natl Acad Sci U S A 69:690–693.
Games PA, Howell JF. 1976;Pairwise multiple comparison procedures with unequal n’s and/or variances: a Monte Carlo study. J Educ Stat 1:113–125.
Gu H. 2011;Morphology, phylogenetic position, and ecophysiology of Alexandrium ostenfeldii (Dinophyceae) from the Bohai Sea, China. J Syst Evol 49:606–616.
Guillard RRL, Hargraves PE. 1993; Stichochrysis immobilis is a diatom, not a chrysophyte. Phycologia 32:234–236.
Guillard RRL, Ryther JH. 1962;Studies of marine planktonic diatoms. I. Cyclotella nana Hustedt and Detonula confervacea (Cleve) Gran. Can J Microbiol 8:229–239.
Guiry MD, Guiry GM. 2021. AlgaeBase Available from: https://www.algaebase.org . Accessed Oct 10 2021.
Haddock SHD, Moline MA, Case JF. 2010;Bioluminescence in the sea. Annu Rev Mar Sci 2:443–493.
Hastings JW, Sweeney BM. 1957;The luminescent reaction in extracts of the marine dinoflagellate, Gonyaulax polyedra . J Cell Comp Physiol 49:209–225.
Higman WA, Stone DM, Lewis JM. 2001;Sequence comparisons of toxic and non-toxic Alexandrium tamarense (Dinophyceae) isolates from UK waters. Phycologia 40:256–262.
Huelsenbeck JP, Ronquist F. 2001;MrBayes: Bayesian inference of phylogenetic trees. Bioinformatics 17:754–755.
Jacobson DM, Anderson DM. 1996;Widespread phagocytosis of ciliates and other protists by marine mixotrophic and heterotrophic thecate dinoflagellates. J Phycol 32:279–285.
Jang SH, Jeong HJ. 2020;Spatio-temporal distributions of the newly described mixotrophic dinoflagellate Yihiella yeosuensis (Suessiaceae) in Korean coastal waters and its grazing impact on prey populations. Algae 35:45–59.
Jeong HJ, Kang HC, Lim AS, Jang SH, Lee K, Lee SY, Ok JH, You JH, Kim JH, Lee KH, Park SA, Eom SH, Yoo YD, Kim KY. 2021;Feeding diverse prey as an excellent strategy of mixotrophic dinoflagellates for global dominance. Sci Adv 7eabe4214.
Jeong HJ, Lim AS, Franks PJS, Lee KH, Kim JH, Kang NS, Lee MJ, Jang SH, Lee SY, Yoon EY, Park JY, Yoo YD, Seong KA, Kwon JE, Jang TY. 2015;A hierarchy of conceptual models of red-tide generation: nutrition, behavior, and biological interactions. Harmful Algae 47:97–115.
Jeong HJ, Lim AS, Lee K, Lee MJ, Seong KA, Kang NS, Jang SH, Lee KH, Lee SY, Kim MO, Kim JH, Kwon JE, Kang HC, Kim JS, Yuh W, Shin K, Jang PK, Ryu J-H, Kim SY, Park JY, Kim KY. 2017;Ichthyotoxic Cochlodinium polykrikoides red tides offshore in the South Sea, Korea in 2014: I. Temporal variations in three-dimensional distributions of red-tide organisms and environmental factors. Algae 32:101–130.
Jeong HJ, Park JY, Nho JH, Park MO, Ha JH, Seong KA, Jeng C, Seong CN, Lee KY, Yih WH. 2005;Feeding by red-tide dinoflagellates on the cyanobacterium Synechococcus . Aquat Microb Ecol 41:131–143.
John U, Litaker RW, Montresor M, Murray S, Brosnahan ML, Anderson DM. 2014;Formal revision of the Alexandrium tamarense species complex (Dinophyceae) taxonomy: the introduction of five species with emphasis on molecular-based (rDNA) classification. Protist 165:779–804.
Kang HC, Jeong HJ, Kim SJ, You JH, Ok JH. 2018;Differential feeding by common heterotrophic protists on 12 different Alexandrium species. Harmful Algae 78:106–117.
Kang HC, Jeong HJ, Ok JH, You JH, Jang SH, Lee SY, Lee KH, Park JY, Rho J-R. 2019;Spatial and seasonal distributions of the phototrophic dinoflagellate Biecheleriopsis adriatica (Suessiaceae) in Korea: quantification using qPCR. Algae 34:111–126.
Kelly MG. 1968;The occurrence of dinoflagellate luminescence at Woods Hole. Biol Bull 135:279–295.
Krasnow R, Dunlap JC, Taylor W, Hastings JW, Vetterling W, Gooch V. 1980;Circadian spontaneous bioluminescent glow and flashing of Gonyau1ax polyedra . J Comp Phvsiol 138:19–26.
Kremp A, Lindholm T, Dreßler N, Erler K, Gerdts G, Eirtovaara S, Leskinen E. 2009;Bloom forming Alexandrium ostenfeldii (Dinophyceae) in shallow waters of the Åland archipelago, Northern Baltic Sea. Harmful Algae 8:318–328.
Kruskal WH, Wallis WA. 1952;Use of ranks in one-criterion variance analysis. J Am Stat Assoc 47:583–621.
Laabir M, Jauzein C, Genovesi B, Masseret E, Grzebyk D, Cecchi P, Vaquer A, Perrin Y, Collos Y. 2011;Influence of temperature, salinity and irradiance on the growth and cell yield of the harmful red tide dinoflagellate Alexandrium catenella colonizing Mediterranean waters. J Plankton Res 33:1550–1563.
Lapota D, Geiger ML, Stiffey AV, Rosenberger DE, Young DK. 1989;Correlations of planktonic bioluminescence with other oceanographic parameters from a Norwegian fjord. Mar Ecol Prog Ser 55:217–227.
Latz MI, Bovard M, VanDelinder V, Segre E, Rohr J, Groisman A. 2008;Bioluminescent response of individual dinoflagellate cells to hydrodynamic stress measured with millisecond resolution in a microfluidic device. J Exp Biol 211:2865–2875.
Latz MI, Jeong HJ. 1996;Effect of red tide dinoflagellate diet and cannibalism on the bioluminescence of the heterotrophic dinoflagellates Protoperidinium spp. Mar Ecol Prog Ser 132:275–285.
Latz MI, Lee AO. 1995;Spontaneous and stimulated bioluminescence of the dinoflagellate Ceratocorys horrzda (Peridiniales). J Phycol 31:120–132.
Latz MI, Nauen JC, Rohr J. 2004;Bioluminescence response of four species of dinoflagellates to fully developed pipe flow. J Plankton Res 26:1529–1546.
Le Tortorec AH, Tahvanainen P, Kremp A, Simis SGH. 2016;Diversity of luciferase sequences and bioluminescence production in Baltic Sea Alexandrium ostenfeldii . Eur J Phycol 51:317–327.
Lee KH, Jeong HJ, Kang HC, Ok JH, You JH, Park SA. 2019;Growth rates and nitrate uptake of co-occurring red-tide dinoflagellates Alexandrium affine and A. fraterculus as a function of nitrate concentration under light-dark and continuous light conditions. Algae 34:237–251.
Lee KH, Jeong HJ, Kwon JE, Kang HC, Kim JH, Jang SH, Park JY, Yoon EY, Kim JS. 2016;Mixotrophic ability of the phototrophic dinoflagellates Alexandrium andersonii, A. affine, and A. fraterculus . Harmful Algae 59:67–81.
Lee SY, Jeong HJ, Ok JH, Kang HC, You JH. 2020;Spatial-temporal distributions of the newly described mixotrophic dinoflagellate Gymnodinium smaydae in Korean coastal waters. Algae 35:225–236.
Lewis NI, Xu W, Jericho SK, Kreuzer HJ, Jericho MH, Cembella AD. 2006;Swimming speed of three species of Alexandrium (Dinophyceae) as determined by digital in-line holography. Phycologia 45:61–70.
Lilly EL, Halanych KM, Anderson DM. 2005;Phylogeny, biogeography, and species boundaries within the Alexandrium minutum group. Harmful Algae 4:1004–1020.
Lim AS, Jeong HJ, Kim JH, Jang SH, Lee MJ, Lee K. 2015;Mixotrophy in the newly described dinoflagellate Alexandrium pohangense: a specialist for feeding on the fast-swimming ichthyotoxic dinoflagellate Cochlodinium polykrikoides . Harmful Algae 49:10–18.
Lim AS, Jeong HJ, Ok JH. 2019;Five Alexandrium species lacking mixotrophic ability. Algae 34:289–301.
Lim AS, Jeong HJ, Ok JH, Kim SJ. 2018;Feeding by the harmful phototrophic dinoflagellate Takayama tasmanica (Family Kareniaceae). Harmful Algae 74:19–29.
Lindström J, Grebner W, Rigby K, Selander E. 2017;Effects of predator lipids on dinoflagellate defense mechanisms-increased bioluminescence capacity. Sci Rep 7:13104.
Litaker RW, Vandersea MW, Kibler SR, Reece KS, Stokes NA, Steidinger KA, Millie DF, Bendis BJ, Pigg RJ, Tester PA. 2003;Identification of Pfiesteria piscicida (Dinophyceae) and Pfiesteria-like organisms using internal transcribed spacer-specific PCR assays. J Phycol 39:754–761.
Liu L, Wilson T, Hastings JW. 2004;Molecular evolution of dinoflagellate luciferases, enzymes with three catalytic domains in a single polypeptide. Proc Natl Acad Sci U S A 101:16555–16560.
MacKenzie L, de Salas M, Adamson J, Beuzenberg V. 2004;The dinoflagellate genus Alexandrium (Halim) in New Zealand coastal waters: comparative morphology, toxicity and molecular genetics. Harmful Algae 3:71–92.
Mann HB, Whitney DR. 1947;On a test of whether one of two random variables is stochastically larger than the other. Ann Math Stat 18:50–60.
Marcinko CLJ, Painter SC, Martin AP, Allen JT. 2013;A review of the measurement and modelling of dinoflagellate bioluminescence. Prog Oceanogr 109:117–129.
Martínez A, Méndez S, Fabre A. 2016;First record of bioluminescence of Alexandrium fraterculus (dinoflagellate), in the Uruguayan coast, South Western Atlantic Ocean. Pan-Am. J Aquat Sci 11:356–360.
McCapra F. 1976;Chemical mechanisms in bioluminescence. Acc Chem Res 9:201–208.
Miller SD, Haddock SHD, Elvidge CD, Lee TF. 2005;Detection of a bioluminescent milky sea from space. Proc Natl Acad Sci U S A 102:14181–14184.
Montresor M, John U, Beran A, Medlin LK. 2004; Alexandrium tamutum sp. nov. (Dinophyceae): a new nontoxic species in the genus Alexandrium . J Phycol 40:398–411.
Murakami M, Makabe K, Yamaguchi K, Konosu S, Wälchli MR. 1988;Goniodomin A, a novel polyether macrolide from the dinoflagellate Goniodoma pseudogoniaulax . Tetrahedron Lett 29:1149–1152.
Murakami M, Okita Y, Matsuda H, Okino T, Yamaguchi K. 1998;From the dinoflagellate Alexandrium hiranoi . Phytochemistry 48:85–88.
Murray S, John U, Savela H, Kremp A. 2020. Alexandrium spp: genetic and ecological factors influencing saxitoxin production and proliferation. In : Botana LM, Louzao MC, Vilarino N, eds. Climate Change and Marine and Freshwater Toxins Walter de Gruyter. Berlin: p. 133–166.
Nguyen-Ngoc L. 2004;An autecological study of the potentially toxic dinoflagellate Alexandrium affine isolated from Vietnamese waters. Harmful Algae 3:117–129.
Orr RJS, Stüken A, Murray SA, Jakobsen KS. 2013;Evolutionary acquisition and loss of saxitoxin biosynthesis in dinoflagellates: the second “core” gene, sxtG. Appl Environ Microbiol 79:2128–2136.
Orr RJS, Stüken A, Rundberget T, Eikrem W, Jakobsen KS. 2011;Improved phylogenetic resolution of toxic and non-toxic Alexandrium strains using a concatenated rDNA approach. Harmful Algae 10:676–688.
Prakash A, Taylor FJR. 1966;A “red water” bloom of Gonyaulax acatenella in the Strait of Georgia and its relation to paralytic shellfish toxicity. J Fish Res Board Can 23:1265–1270.
Rohr J, Latz MI, Fallon S, Nauen JC, Hendricks E. 1998;Experimental approaches towards interpreting dolphin-stimulated bioluminescence. J Exp Biol 201:1447–1460.
Ronquist F, Huelsenbeck JP. 2003;MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19:1572–1574.
Satake M, MacKenzie L, Yasumoto T. 1997;Identification of Protoceratium reticulatum as the biogenetic origin of yessotoxin. Nat. Toxins 5:164–167.
Scholin CA, Herzog M, Sogin M, Anderson DM. 1994;Identification of group- and strain-specific genetic makers for globally distributed Alexandrium (Dinophyceae). II. Sequence analysis of a fragment of the LSU rRNA gene. J Phycol 30:999–1011.
Seliger HH, Biggley WH, Swift E. 1969;Absolute values of photon emission from the marine dinoflagellates Pyrodinium bahamense, Gonyaulax polyedra and Pyrocystis lunula. Photochem. Photobiol 10:227–232.
Shimomura O. 2006. Bioluminescence: chemical principles and methods World Scientific. Toh Tuck Link: p. 500.
Stamatakis A. 2006;RaxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22:2688–2690.
Subong BJJ, Benico GA, Sulit AKL, Mendoza CO, Cruz LJ, Azanza RV, Jimenez EC. 2017;Toxicity and protein expression of Alexandrium species collected in the Philippine waters. Philipp J Sci 146:425–436.
Sullivan JM, Swift E. 1995;Photoenhancement of bioluminescence capacity in natural and laboratory populations of the autotrophic dinoflagellate Ceratium fusus (Ehrenb.) Dujardin. J. Geophys. Res. Oceans 100:6565–6574.
Sweeney BM. 1963;Bioluminescent dinoflagellates. Biol Bull 125:177–181.
Sweeney BM. 1986;The loss of the circadian rhythm in photosynthesis in an old strain of Gonyaulax polyedra . Plant Physiol 80:978–981.
Swift E, Biggley WH, Seliger HH. 1973;Species of oceanic dinoflagellates in the genera Dissodinium and Pyrocystis: interclonal and interspecific comparisons of the color and photon yield of bioluminescence. J Phycol 9:420–426.
Swift E, Sullivan JM, Batchelder HP, Van Keuren J, Vaillancourt RD, Bidigare RR. 1995;Bioluminescent organisms and bioluminescence measurements in the North Atlantic Ocean near latitude 59.5°N, longitude 21°W. J Geophys Res Oceans 100:6527–6547.
Tamura K, Dudley J, Nei M, Kumar S. 2007;MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol 24:1596–1599.
Valiadi M, Iglesias-Rodriguez D. 2013;Understanding bioluminescence in dinoflagellates: how far have we come? Microorganisms 1:3–25.
Valiadi M, Iglesias-Rodriguez MD, Amorim A. 2012;Distribution and genetic diversity of the luciferase gene within marine dinoflagellates. J Phycol 48:826–836.
Von Dassow P, Latz MI. 2002;The role of Ca2+ in stimulated bioluminescence of the dinoflagellate Lingulodinium polyedrum . J Exp Biol 205:2971–2986.
Welch BL. 1947;The generalization of ‘student’s’ problem when several different population variances are involved. Biometrika 34:28–35.
White HH. 1979;Effects of dinoflagellate bioluminescence on the ingestion rates of herbivorous zooplankton. J Exp Mar Biol Ecol 36:217–224.
Widder EA. 2010;Bioluminescence in the ocean: origins of biological, chemical, and ecological diversity. Science 328:704–708.
Widder EA, Case JF. 1981;Two flash forms in the bioluminescent dinoflagellate, Pyrocystis fusiformis . J Comp Physiol 143:43–52.
Wilson T, Hastings JW. 2013. Bioluminescence: living lights, lights for living Harvard University Press. Cambridge, MA: p. 185.
Wohlrab S, Iversen MH, John U. 2010;A molecular and co-evolutionary context for grazer induced toxin production in Alexandrium tamarense . PLoS ONE 5:e15039.
Yoo YD, Jeong HJ, Kim MS, Kang NS, Song JY, Shin W, Kim KY, Lee K. 2009;Feeding by phototrophic red-tide dinoflagellates on the ubiquitous marine diatom Skeletonema costatum . J Eukaryot Microbiol 56:413–420.

Article information Continued

Fig. 1

Photographs of bioluminescence showing the bioluminescent capability of the Alexandrium species tested in the present study. A 250 mL culture flask containing a dense culture of each species in the late exponential stage was shaken by hand. (A) Alexandrium mediterraneum AMYS1807 (30,000 cells mL−1). (B) Alexandrium pohangense APPH1409 (2,000 cells mL−1). (C) Alexandrium tamutum ATSH1609 (10,000 cells mL−1). (D) Alexandrium andersonii CCMP2222 (7,000 cells mL−1). (E) Alexandrium hiranoi NIES-3611 (1,000 cells mL−1). (F) Alexandrium pseudogonyaulax KSUDi-noB4 (1,000 cells mL−1). (G) Alexandrium insuetum CCMP2082 (30,000 cells mL−1).

Fig. 2

Consensus Bayesian tree of nuclear large subunit ribosomal DNA region (based on 890 bp aligned position), using the GTR + G + I model. The bioluminescent dinoflagellates Lingulodinium polyedra, Protoceratium reticulatum, and Gonyaulax spinifera were used as the outgroups. The numbers above the branches indicate the Bayesian posterior probability (≥0.5, left) and maximum likelihood bootstrap values (≥50%, right), and solid black circles indicate full support (1.0 and 100%). Bioluminescent strains are indicated by a blue closed star, whereas nonbioluminescent strains are indicated by a purple open star. The strains tested in the present study are written in bold.

Fig. 3

Consensus Bayesian tree, based on 498 bp aligned position of nuclear luciferase gene (lcf) (A) and based on 890 bp aligned position of the large subunit ribosomal DNA (B), using the GTR + G + I model and Lingulodinium polyedra as the outgroup. The numbers above the branches indicate the Bayesian posterior probability (≥0.5, left) and maximum likelihood bootstrap values (≥50%, right), and solid black circles indicate full support (1.0 and 100%). The sequences obtained in the present study are in boldface.

Fig. 4

Bioluminescence intensity of each Alexandrium species. (A) The mean 200-second-integrated bioluminescence intensity per cell (BLInt200s). Significantly different subsets of the Alexandrium species based on Games-Howell post-hoc tests (p < 0.05) of Welch’s ANOVA: (a) Alexandrium tamarense (Ata) and A. mediterraneum (Ame); (b) A. tamutum (Atm) and A. pacificum (Apa); (abc) A. affine (Aaf) and A. fraterculus (Afr); (c) A. ostenfeldii (Aos) and A. pohangense (Apo). (B) The mean maximum bioluminescence intensity per cell per second (BLMax). Significantly different subsets of Alexandrium species based on the Bonferroni post-hoc test (p < 0.05) of Kruskal-Wallis test: (a’) Ata; (a’b’) Ame and Atm; (a’b’c’) Apa, Aaf, and Apo; (b’c’) Afr; (c’) Aos. RLU, relative luminescence unit. Values represent the treatment means ± standard error (n = 4).

Fig. 5

Comparison of the presence (closed symbols) or absence (open symbols) of bioluminescent capacity (BL), toxicity (Toxi), and mixotrophic ability (Mixo) in the Consensus Bayesian tree of the nuclear large subunit ribosomal DNA region in Fig. 2. Asterisks indicate that some strains are toxic, whereas others are nontoxic (Prakash and Taylor 1966, Murakami et al. 1988, 1998, Jacobson and Anderson 1996, Satake et al. 1997, Draisci et al. 1999, Ciminiello et al. 2000, Higman et al. 2001, MacKenzie et al. 2004, Montresor et al. 2004, Nguyen-Ngoc 2004, Jeong et al. 2005, Lilly et al. 2005, Yoo et al. 2009, Orr et al. 2011, 2013, Blossom et al. 2012, 2017, John et al. 2014, Lim et al. 2015, 2019, Eckford-Soper et al. 2016, Lee et al. 2016, Subong et al. 2017, Murray et al. 2020).

Fig. 6

Mean 200-second-integrated bioluminescence intensity per cell (BLInt200s) (A) and mean maximum bioluminescence intensity per cell per second (BLMax) (B) of bioluminescent Alexandrium species as a function of cell size (equivalent spherical diameter [ESD], μm), as described in Table 2. The values represent the treatment means ± standard error (n = 4). The Spearman correlation coefficient and p-value were described in the figures. Aaf, A. affine; Afr, A. fraterculus; Ame, A. mediterraneum; Aos, A. ostenfeldii; Apa, A. pacificum; Apo, A. pohangense; Ata, A. tamarense; Atm, A. tamutum.

Fig. 7

Mean 200-second-integrated bioluminescence intensity per cell (BLInt200s) (A) and the mean maximum bioluminescence intensity per cell per second (BLMax) (B) of the bioluminescent Alexandrium species as a function of maximum growth rate (MGR; d−1). Values represent the treatment means ± standard error (n = 4). a: the Spearman correlation coefficient and p-value of all species. b: the Spearman correlation coefficient and p-value of all species except for A. tama-rense (red circle). Aaf, A. affine; Afr, A. fraterculus; Aos, A. ostenfeldii; Apa, A. pacificum; Apo, A. pohangense; Ata, A. tamarense. MGR data were obtained from Gu (2011), Laabir et al. (2011), Lim et al. (2015), Bill et al. (2016), and Lee et al. (2019).

Fig. 8

Mean 200-second-integrated bioluminescence intensity per cell (BLInt200s) (A) and the mean maximum bioluminescence intensity per cell per second (BLMax) (B) of the bioluminescent Alexandrium species as a function of maximum swimming speed (MSS; μm s−1). Values represent the treatment means ± standard error (n = 4). The Spearman correlation coefficient and p-value were described in the figures. Aaf, A. affine; Afr, A. fraterculus; Ame, A. mediterraneum; Aos, A. ostenfeldii; Apa, A. pacificum; Apo, A. pohangense; Ata, A. tamarense; Atm, A. tamutum. MSS data were obtained from Lewis et al. (2006), Lim et al. (2015, 2018), Lee et al. (2016), and Kang et al. (2018).

Table 1

Information on the isolation and maintenance of the experimental organisms and also Alexandrium species in which bioluminescence capability was previously tested or newly tested

Organisms Strain name Location Date T S Media / Prey Biolu
Alexandrium andersonii CCMP2222 Gulf of Naples, Italy Aug 1, 1996 - - L1 NewBL
A. minutum CCMP1888 Laguna Obidos, Portugal - - - L1 PNB
A. tamutum ATSH1609 Shiwha, Korea Sep 9, 2016 25.9 32.6 F/2 NewBL
A. insuetum CCMP2082 Uchiumi Bay, Japan Jun 6, 1985 - - F/2 NewBL
A. ostenfeldii NIES-4274 Lake Koyama, Japan Nov 29, 2016 - - F/2 PB
A. pacificum KSUDinoE4 Kunsan, Korea Sep 2020 23.5 29.2 F/2 PB
A. affine AATA1308 Taean, Korea Aug 14, 2013 21.5 32.2 F/2 PB
A. pohangense APPH1409 Pohang, Korea Sep 2014 23.3 31.1 Mp NewBL
A. fraterculus AFYS1309 Yeosu, Korea Sep 6, 2013 23.4 32.8 L1 PB
A. mediterraneum AMYS1807 Yeosu, Korea Jul 2018 24.2 30.0 F/2 NewBL
A. tamarense CCMP115 English Channel, UK Jun 24, 1957 - - F/2 PB
A. pseudogonyaulax KSUDinoB4 Kunsan, Korea Sep 2020 23.5 29.2 F/2 NewBL
A. hiranoi NIES-3611 Muroto Cape, Japan Apr 29, 2008 - - F/2 NewBL
Protoceratium reticulatum RCC4104 Pacific Ocean, Japan Aug 30, 2013 - - L1 PB
Lingulodinium polyedra CCMP1931 La Jolla, California, USA May 1998 - - L1 PB

Culture information was obtained from the National Center for Marine Algae and Microbiota (NCMA; USA), Roscoff Culture Collection (RCC; France), and Microbial Culture Collection at the National Institute for Environmental Studies (NIES; Japan).

T, temperature (°C); S, salinity; NewBL, newly tested for bioluminescence in the present study; PNB, previously known to be not bioluminescent; PB, previously known to be bioluminescent; Mp, Margalefidinium polykrikoides.

Table 2

Mean 200-second-integrated bioluminescence intensity per cell (BLInt200s, ×104 RLU cell−1 ± standard error) and mean maximum bioluminescence intensity per cell per second (BLMax, ×104 RLU cell−1 s−1 ± standard error) of the 13 Alexandrium species tested in the present study

Species Strain name ESD BC BLInt200s BLMax
Alexandrium andersonii CCMP2222 14.9 × 0.00 ± 0.0
A. minutum CCMP1888 20.4 × 0.00 ± 0.0
A. insuetum CCMP2082 26.4 × 0.00 ± 0.0
A. pseudogonyaulax KSUDinoB4 35.4 × 0.00 ± 0.0
A. hiranoi NIES-3611 40.0 × 0.00 ± 0.0
A. tamutum ATSH1609 22.5 1.54 ± 0.1 0.19 ± 0.1
A. ostenfeldii NIES-4274 29.8 32.2 ± 4.5 10.3 ± 1.3
A. pacificum KSUDinoE4 30.3 1.55 ± 0.1 0.34 ± 0.1
A. affine AATA1308 31.4 9.23 ± 2.1 2.80 ± 0.7
A. pohangense APPH1409 32.0 19.2 ± 2.0 2.75 ± 0.5
A. fraterculus AFYS1309 32.3 12.7 ± 2.1 4.53 ± 0.9
A. mediterraneum AMYS1807 33.0 0.70 ± 0.1 0.14 ± 0.0
A. tamarense CCMP115 34.6 0.02 ± 0.0 0.01 ± 0.0

RLU, relative luminescence unit; ESD, equivalent spherical diameter (μm); BC, bioluminescence capability; ○, bioluminescent; ×, not bioluminescent.

Table 3

Comparison of bioluminescence capability (BC), mixotrophy, and toxicity of Alexandrium species

Species Strain BC Mixotrophy Toxicity Reference
Alexandrium insuetum CCMP2082 × × × Orr et al. (2013), Lim et al. (2019)
A. minutum CCMP1888 × Jeong et al. (2005), Orr et al. (2011)
A. andersonii CCMP2222 × × Orr et al. (2011), Lee et al. (2016)
A. pseudogonyaulax KSUDinoB4 × a a Murakami et al. (1988), Blossom et al. (2012, 2017)
A. hiranoi NIES-3611 × ×a Murakami et al. (1998), Blossom et al. (2012)
A. acatenella - a Prakash and Taylor (1966), Esaias and Curl (1972)
A. affine CCMP112 × × Orr et al. (2011), Lee et al. (2016)
A. fraterculus AFYS1309 × ×a MacKenzie et al. (2004), Lee et al. (2016)
A. catenella CCMP1719 a Jeong et al. (2005), Orr et al. (2013), Lim et al. (2019)
A. mediterraneum AMYS1807 ×a ×a John et al. (2014), Lim et al. (2019)
A. monilatum - a Sweeney (1963), MacKenzie et al. (2004)
A. tamarense CCMP115 a × Jeong et al. (2005), Orr et al. (2013), Lim et al. (2019)
A. pacificum KSUDinoE4 ×a a John et al. (2014), Lim et al. (2019), Murray et al. (2020)
A. tamutum ATSH1609 ×a ×a Montresor et al. (2004), Lim et al. (2019)
A. ostenfeldii NIES-4274 a a Jacobson and Anderson (1996), Lilly et al. (2005)
A. pohangense APPH1409 - Lim et al. (2015)

The data were obtained from the present study and from the literature.

○, presence of bioluminescent capability, mixotrophic ability, or toxicity; ×, absence of bioluminescent capability, mixotrophic ability, or toxicity.

a

The strains used for testing mixotrophic ability and toxicity differed from the strain used for testing bioluminescence.