:: Algae

Jang and Jeong: Spatio-temporal distributions of the newly described mixotrophic dinoflagellate Yihiella yeosuensis (Suessiaceae) in Korean coastal waters and its grazing impact on prey populations

Research Article

Algae 2020; 35(1): 45-59.

Published online: March 15, 2020

DOI: https://doi.org/10.4490/algae.2020.35.2.24

Spatio-temporal distributions of the newly described mixotrophic dinoflagellate Yihiella yeosuensis (Suessiaceae) in Korean coastal waters and its grazing impact on prey populations

Se Hyeon Jang¹, Hae Jin Jeong^1,^2,^*

¹School of Earth and Environmental Sciences, College of Natural Sciences, Seoul National University, Seoul 08826, Korea

²Research Institute of Oceanography, Seoul National University, Seoul 08826, Korea

^*Corresponding Author: E-mail: hjjeong@snu.ac.kr,
Tel: +82-2-880-6746, Fax: +82-2-874-9695

Received January 10, 2020 Accepted February 24, 2020

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

ABSTRACT

To investigate the spatio-temporal distributions of the mixotrophic dinoflagellate Yihiella yeosuensis in Korean coastal waters and its grazing impact on prey populations, water samples were seasonally collected from 28 stations in the East, West, and South Seas of Korea and Jeju Island from April 2015 to October 2018. The abundances of Y. yeosuensis in the water samples were quantified using quantitative real-time polymerase chain reaction (qPCR). Simultaneously, the physical and chemical properties of water from all sampled stations were determined, and the abundances of the optimal prey species of Y. yeosuensis, the prasinophyte Pyramimonas sp. and the cryptophyte Teleaulax amphioxeia, were quantified using qPCR. Y. yeosuensis has a wide distribution, as is reflected by the detection of Y. yeosuensis cells at 23 sampling stations; however, this distribution has a strong seasonality, which is indicated by its detection at 22 stations in summer but only one station in winter. The abundance of Y. yeosuensis was significantly and positively correlated with those of Pyramimonas sp. and T. amphioxeia, as well as with water temperature. The highest abundance of Y. yeosuensis was 48.5 cells mL⁻¹ in Buan in July 2017, when the abundances of Pyramimonas sp. and T. amphioxeia were 917.6 and 210.4 cells mL⁻¹, respectively. The growth rate of Y. yeosuensis on Pyramimonas sp., calculated by interpolating the growth rates at the same abundance, was 0.49 d⁻¹, which is 37% of the maximum growth rate of Y. yeosuensis on Pyramimonas sp. obtained in the laboratory. Therefore, the field abundance of Pyramimonas sp. obtained in the present study can support a moderate positive growth of Y. yeosuensis. The maximum grazing coefficient for Y. yeosuensis on the co-occurring Pyramimonas sp. was 0.42 d⁻¹, indicating that 35% of the Pyramimonas sp. population were consumed in 1 d. Therefore, the spatio-temporal distribution of Y. yeosuensis in Korean coastal waters may be affected by those of the optimal prey species and water temperature. Moreover, Y. yeosuensis may potentially have considerable grazing impacts on populations of Pyramimonas sp.

Key words: food web; harmful algal bloom; mixotrophy; qPCR; red tide; temperature

INTRODUCTION

Dinoflagellates are a major component of marine ecosystems (Sherr and Sherr 2007, Taylor et al. 2008, Jeong et al. 2013). They play diverse roles in marine food webs (Hansen 1992, Coats 1999, Jeong 1999, 2010b, Yoo et al. 2009, Lee et al. 2016b, Stoecker et al. 2017, Kang et al. 2019a). Many autotrophic and mixotrophic dinoflagellates and a few heterotrophic dinoflagellates form red tides or harmful algal blooms that can alter the balance of marine food webs and cause large-scale mortalities of fish and shellfish (Glasgow et al. 1995, Park et al. 2013b, Jeong et al. 2017, Shumway et al. 2018). Thus, the distributions of dinoflagellates are a critical concern for scientists, government officials, aquaculture farmers, and the public (Burkholder et al. 1995, Jeong et al. 2006, Taylor et al. 2008, Kudela and Gobler 2012, Park et al. 2013a). In particular, when a dinoflagellate is newly described, its distribution should be explored.

For a long time, phototrophic dinoflagellates had been treated as obligate autotrophic dinoflagellates; thus, scientists had mainly studied their abiotic growth factors, such as light and nutrient availability (Eppley 1972, Egge and Aksnes 1992, Kim et al. 2004, Lee et al. 2017a, 2019a). However, in the last three decades, many phototrophic dinoflagellates have turned out to be mixotrophic dinoflagellates (Stoecker et al. 1997, Jeong et al. 1999, 2004, 2005a, 2005b, 2010a, 2010b, 2015, 2016, Adolf et al. 2008, Berge et al. 2008, Glibert et al. 2009, Lee et al. 2015, 2016a). Furthermore, several newly described dinoflagellates have been revealed to be mixotrophic (Yoo et al. 2010, Kang et al. 2011, Jeong et al. 2012, Lee et al. 2014a, 2014b, Lim et al. 2015). Thus, the biotic factor of prey availability should be considered to understand population dynamics of mixotrophic dinoflagellates (Hansen 2011, Jeong et al. 2012, 2013, 2015, Stoecker et al. 2017).

Identifying dinoflagellates under a light microscope is very difficult for scientists, especially when the sizes of these organisms are less than 20 μm and the morphological differences among dinoflagellate species within a genus are small. Molecular techniques, such as quantitative real-time polymerase chain reaction (qPCR), partially solve the difficulty of differentiating dinoflagellate species (Smith et al. 2016, Lee et al. 2017b, 2019b, Kang et al. 2019b). Thus, to quantify the abundance of a small dinoflagellate and subsequently understand its distribution, molecular techniques should be used.

The dinoflagellate Yihiella yeosuensis, belonging to the family Suessiaceae, was newly described in 2017 (Jang et al. 2017b). It is one of the smallest dinoflagellates, having a length of 10 μm. Recently, Y. yeosuensis was revealed to be mixotrophic, and this mixotrophy significantly increased its growth rate (Jang et al. 2017a). However, Y. yeosuensis could feed only on the prasinophyte Pyramimonas sp. and the cryptophyte Teleaulax amphioxeia among 19 tested algal prey species (Jang et al. 2017a). In addition, Y. yeosuensis exhibits an anti-predation behavior, fast-backward jumping, against common heterotrophic protists (Jeong et al. 2018). Thus, the distribution of Y. yeosuensis may be affected by prey availability and abiotic factors, and relationships between Y. yeosuensis distribution and environmental factors should be investigated.

In this study, to investigate the spatio-temporal distributions of the mixotrophic dinoflagellate Y. yeosuensis in Korean coastal waters and the grazing impact of Y. yeosuensis on prey populations, water samples were collected seasonally from 28 stations in the East, West, and South Seas of Korea and Jeju Island from April 2015 to October 2018, and the abundances of Y. yeosuensis and the optimal prey Pyramimonas sp. and T. amphioxeia were quantified using qPCR. Simultaneously, physical and chemical properties of the waters at all sampled stations, such as water temperature, salinity, nutrient concentrations, and dissolved oxygen concentration, were measured. Furthermore, the grazing impact of Y. yeosuensis on Pyramimonas sp. populations was calculated by combining field data on the abundances of the predator and prey with ingestion rates of Y. yeosuensis on Pyramimonas sp. reported by Jang et al. (2017a). The results of this study provide a basis for understanding the distribution of Y. yeosuensis, the critical environmental factors affecting this distribution, and also the grazing impact of Y. yeosuensis on prey populations.

MATERIALS AND METHODS

Preparation of experimental organisms

Cells of Y. yeosuensis YYYS1405 were isolated from surface waters off Yeosu, in the South Sea of Korea (34°32′ N, 127°55′ E), when the water temperature and salinity were 18.6°C and 33.8, respectively (Jang et al. 2017b). Cells of Pyramimonas sp. PSSH1204 were isolated from surface waters off Shiwha Bay, West Sea of Korea, when the water temperature and salinity were 15.3°C and 27.4, respectively. Cells of T. amphioxeia CR-MAL01 were isolated from surface waters off Gomso Bay, in the West Sea of Korea, when the water temperature and salinity were 7.8°C and 30.1, respectively (Yih et al. 2004). A clonal culture of each of these three species was used for qPCR, following two consecutive single-cell isolations. Pyramimonas sp. and T. amphioxeia were maintained in enriched f/2 seawater medium, while Y. yeosuensis was provided Pyramimonas sp. as prey. After the concentration of each species had increased sufficiently, the volume was increased to 32, 270, and 500 mL in polycarbonate bottles. The bottles were filled to capacity, capped, and placed on a shelf at 20°C, illuminated with an irradiance of 20 μmol photons m⁻² s⁻¹ provided by cool white fluorescent light, under a 14 : 10-h light-dark cycle.

Primer design and specificity testing

The internal transcribed spacer ribosomal (r) DNA sequences of Y. yeosuensis YYYS1405 and related dinoflagellate species obtained from GenBank were aligned using the program MEGA v.4 (Tamura et al. 2007). Additionally, the large subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase (rbcL) gene sequence of Pyramimonas sp. PSSH1204 and related prasinophyte species obtained from GenBank were aligned. Finally, the large subunit rDNA sequences of T. amphioxeia CR-MAL01 and related cryptophyte species obtained from GenBank were aligned. Manual searches of the alignments were conducted to determine the unique sequences of Y. yeosuensis, Pyramimonas sp., and T. amphioxeia and thereby develop target species-specific qPCR assays using primer-probe sets able to discriminate the target species from other related species. For determining the optimal melting temperatures and secondary structures, the primer and probe sequences for each target species were analyzed using Primer 3 (Whitehead Institute, Cambridge, MA, USA; Howard Hughes Medical Institute, Cambridge, MA, USA; Howard Hughes Medical Institute, Chevy Chase, MD, USA) and Oligo Calc: Oligonucleotide Properties Calculator (Kibbe 2007). Subsequently, primers and probes were synthesized by BioResearch Technologies (Petaluma, CA, USA). The probes were dual-labeled at the 5′ and 3′ ends with the fluorescent dyes FAM and BHQ1 (BioResearch Technologies), respectively (Table 1). A species specificity test for the primers and probe for each target species was conducted using related species (Supplementary Table S1).

Construction of standard curves and qPCR analysis

Standard curves were constructed using DNA obtained from dense cultures of the three species harvested by filtration through a 25-mm GF/C filter (Whatman Inc., Clifton, NJ, USA). Each culture in the exponential growth phase was used for qPCR. Specifically, DNA was extracted from cultures of Y. yeosuensis, Pyramimonas sp., and T. amphioxeia, each containing 100,000 cells in the final elution volume of 100 μL, using an AccuPrep Genomic DNA extraction kit (Bioneer, Daejeon, Korea), according to the manufacturer’s instructions. The extracted DNA was then serially diluted by adding predetermined volumes of deionized sterile water (DDW) (Bioneer) to 1.5-mL tubes to ultimately prepare a 6-fold dilution series for all species, targeting 100, 10, 1, 0.1, 0.01, and 0.001% of the initially extracted DNA. Subsequently, the DNA samples were stored at −20°C, and a qPCR amplification was conducted within a day. The qPCR assays used to prepare the standard curve were performed using 1 μL of the DNA template, 0.2 μM of primers (forward and reverse), 0.15 μM of the probe (final concentrations), and 5 μL of the qPCRBIO Probe Separate-ROX (Genepole, Gwangmyeong, Korea) prepared with DDW to achieve a final volume of 10 μL. The qPCR assays were performed using a Rotor-Gene Q (Qiagen, Hilden, Germany) and the following thermal cycling conditions: 3 min at 95°C, followed by 45 cycles for 10 s at 95°C and 40 s at 58°C.

Quantification of the abundances of Yihiella yeosuensis and its optimal prey species in Korean waters

For qPCR analysis, 50–300 mL of water samples collected from each station in each time interval were filtered through a 25-mm GF/C filter (Whatman Inc., Clifton, NJ, USA). The filter had been loosely rolled and placed into a 1.5-mL tube and frozen at −20°C until it was transported to the laboratory. For each sample, DNA was extracted from the cells captured on the filters using an AccuPrep Genomic DNA extraction kit (Bioneer), according to the manufacturer’s instructions.

To determine the abundances of Y. yeosuensis, Pyramimonas sp., and T. amphioxeia in field samples, the aforementioned qPCR assay conditions were used. The DNA from each sample was amplified four times to ensure the accuracy of the results. A sample containing DDW as the template was used as a negative control, whereas one used for standard curve construction was used as a positive and standard control.

Hydrographical properties

Using a clean bucket, surface water samples were collected from 28 stations located in the East, West, and South Seas of Korea and Jeju Island, Korea, in April, July, and October of 2015; January, March, July, October, and December of 2016; March, July, and October of 2017; and January, March, July, and October of 2018 (Fig. 1). As the water depth at most stations was less than 3 m, only surface samplings of the stations were performed.

The temperature, salinity, and dissolved oxygen (DO) of the sampled water were measured using a YSI Professional Plus instrument (YSI Inc., Yellow Springs, OH, USA). For nutrient concentration analysis, the water samples were gently filtered through GF/F filters (Whatman Inc., Florham Park, NJ, USA) and stored at −20°C until the concentrations of nitrate plus nitrite (NO₃ + NO₂, hereafter NO₃), phosphate (PO₄), and silicate (SiO₂) were measured using a nutrient auto-analyzer system (QuAAtro; Seal Analytical Gmbh, Norderstedt, Germany). The hydrographical profiles for the monitoring conducted during 2015–2018 were also described in our previous study (Kang et al. 2019b).

Grazing impacts

The grazing coefficients attributable to Y. yeosuensis feeding on the co-occurring Pyramimonas sp. were calculated by combining field data on the abundances of Y. yeosuensis and Pyramimonas sp. obtained in this study with ingestion rates of Y. yeosuensis on Pyramimonas sp. obtained by Jang et al. (2017a). The grazing coefficient g (d⁻¹) was calculated using the following equation:

(1)

g=CR×PC×24

, where CR is the clearance rate (mL predator⁻¹ h⁻¹) of Y. yeosuensis for Pyramimonas sp. at a given prey concentration and PC is the predator concentration (cells mL⁻¹). CR values were calculated using the following equation:

(2)

CR=IR/x

, where IR is the ingestion rate (cells eaten predator⁻¹ h⁻¹) of Y. yeosuensis for Pyramimonas sp. and x is the prey concentration (cells mL⁻¹). CR values were corrected using Q₁₀ = 2.8 (Hansen et al. 1997) because the in situ water temperatures and the temperatures used in the laboratory for the experiment were different.

Data analysis

Statistical analyses were performed on a personal computer using the statistical software package SPSS ver. 23.0 (IBM Corp., Armonk, NY, USA). To investigate a one by one relationship between cell abundances and hydrographic factors, the Pearson correlation analysis test was used.

RESULTS

Spatial and temporal distributions of Yihiella yeosuensis and its optimal prey species in Korean waters

Y. yeosuensis showed a wide distribution in Korean coastal waters through the qPCR method. Y. yeosuensis was detected at 23 of 28 sampling stations in the study period (Fig. 2, Supplementary Table S2). Y. yeosuensis was not detected in the waters of the Seocheon (West Sea), Sokcho (East Sea), Wimi, Seongsan, and Gimnyeong stations (Jeju Island). The highest abundance of Y. yeosuensis (48.5 cells mL⁻¹) was recorded at the Buan station (West Sea) in July 2017. One of the optimal prey species for Y. yeosuensis, Pyramimonas sp. was detected at all the stations (Fig. 3). The highest abundance of Pyramimonas sp. (917.6 cells mL⁻¹) was obtained at the Buan station (West Sea) in July 2017. The other optimal prey species, T. amphioxeia, was also found at all the sampling stations (Fig. 4). The highest abundance of T. amphioxeia (666.5 cells mL⁻¹) was found at the Seogwipo station (Jeju Island) in July 2018. Furthermore, the highest combined abundance of Pyramimonas sp. and T. amphioxeia (1,127.9 cells mL⁻¹) was also found at the Buan station (West Sea) in July 2017.

Y. yeosuensis was identified at 22 stations in the summer but only at three stations in the spring and autumn and only one station in winter, showing a strong seasonality (Fig. 2). Furthermore, four of the five highest abundances of Y. yeosuensis were found in the summer (Supplementary Table S2).

Cells of Pyramimonas sp. were found in all four seasons (Fig. 3). However, the five highest abundances of Pyramimonas sp. were found in the summer (Supplementary Table S2). Cells of T. amphioxeia were also found in all four seasons, but four of the five highest abundances of T. amphioxeia were observed in the summer (Fig. 4, Supplementary Table S2).

Hydrographical properties of assessed stations during the study period

During the study period, the water temperatures at all the stations ranged from 0.2 to 28.0°C, while the salinity at all the stations ranged from 0.1 to 35.6 (Table 2). The concentrations of NO₃ and PO₄ at all the stations ranged from not detectable (ND) to 149.0 μM and from ND to 6.3 μM, respectively, while the concentration of SiO₂ ranged from ND to 453.4 μM. The concentration of DO at all the stations ranged from 0.2 to 14.8 mg L⁻¹ (Table 2).

During the study period, Y. yeosuensis, Pyramimonas sp., and T. amphioxeia were present in Korean coastal waters when water temperatures were 2.3–28.0, 3.5–27.0, and 1.1–28.0°C, respectively; salinities were 9.9–34.4, 9.9–34.6, and 0.1–35.6, respectively (Fig. 5A–C); concentrations of NO₃ were 0.8–91.9, ND-149.0, and ND-118.3 μM, respectively; concentrations of PO₄ were 0.1–2.7, ND-6.3, and ND-6.3 μM, respectively; concentrations of SiO₂ were 2.1–316.6, ND-453.4, and ND-448.4 μM, respectively (Fig. 5D–F); and concentrations of DO were 1.7–11.3, 0.4–12.5, and 0.2–14.8 mg L⁻¹, respectively.

Additionally, during the study period, the highest abundances of Y. yeosuensis, Pyramimonas sp., and T. amphioxeia were found when water temperatures were 25.0, 25.0, and 23.3°C, respectively (Fig. 5A–C); salinities were 29.2, 29.2, and 22.4, respectively (Fig. 5A–C); concentrations of NO₃ were 15.2, 15.2, and 53.4 μM, respectively; concentrations of PO₄ were 1.8, 1.8, and 0.7 μM, respectively; concentrations of SiO₂ were 53.5, 53.5, and 300.2 μM, respectively (Fig. 5D–F); and concentrations of DO were 6.1, 6.1, and 7.0 mg L⁻¹, respectively.

Relationships between the abundance of Yihiella yeosuensis and environmental factors

Pearson’s correlation analysis showed that the abundance of Y. yeosuensis was significantly and positively correlated with those of Pyramimonas sp. (p < 0.01, r = 0.968) and T. amphioxeia (p < 0.01, r = 0.352) (Table 3, Fig. 6). Furthermore, the abundance of Y. yeosuensis was significantly and positively correlated with the water temperature (p < 0.05, r = 0.1) (Table 3). However, there were no significant correlations between the abundance of Y. yeosuensis and water salinity, concentrations of nutrients, or DO concentration (Table 3).

The abundance of Y. yeosuensis was affected mainly by the abundances of Pyramimonas sp. and T. amphioxeia. We further investigated the relationships between environmental factors and the distributions of Pyramimonas sp. and T. amphioxeia. The abundance of T. amphioxeia was significantly and positively correlated with the water temperature and the concentrations of NO₃ and SiO₂ but was significantly and negatively correlated with salinity (Table 3). However, the abundance of Pyramimonas was not significantly correlated with any of the investigated environmental factors (Table 3).

Grazing impacts of Yihiella yeosuensis on populations of Pyramimonas sp

When the abundances of Pyramimonas sp. and Y. yeosuensis were 0.0–917.6 and 0.0–48.5 cells mL⁻¹, respectively, the calculated grazing coefficients (g) attributable to Y. yeosuensis on the co-occurring Pyramimonas sp. were up to 0.42 d⁻¹ (mean ± standard error = 0.03 ± 0.01 d⁻¹, n = 31) (Fig. 7).

DISCUSSION

This study showed that Y. yeosuensis cells are widely distributed in Korean coastal waters. Prior to our study, the presence of Y. yeosuensis had been reported at only one location (34°32′ N, 127°55′ E) (Jang et al. 2017b). Thus, this study extends the latitudinal range in which Y. yeosuensis is present from 34°32′ N to 33°17′-38°12′ N.

During the study period, the ranges of water temperature, salinity, and the concentrations of NO₃, PO₄, and SiO₂ at the stations in which Y. yeosuensis was present were wide. However, the highest abundance of Y. yeosuensis was found in summer when the water temperature was 25.0°C. Meanwhile, the water temperature, salinity, and concentrations of NO₃, PO₄, and SiO₂ at the stations in which Pyramimonas sp. and T. amphioxeia were found also tended to vary widely. However, the highest abundances of these two species were also found in the summer when water temperatures were 23.3–25.0°C. The maximum growth rates of Y. yeosuensis on Pyramimonas sp. and T. amphioxeia (i.e., mixotrophic growth rate) were 1.32 and 1.29 d⁻¹, respectively, but the growth rate of Y. yeosuensis without prey was negligible (Jang et al. 2017a). Thus, Y. yeosuensis was expected to be abundant only when its prey is abundant (Jang et al. 2017a). Significant positive correlations between Y. yeosuensis and Pyramimonas sp. or T. amphioxeia found in this study support this hypothesis. During the study period, the maximum abundance of Pyramimonas sp. was 917.6 cells mL⁻¹ (36.7 ng C mL⁻¹). Simultaneously, the abundance of Y. yeosuensis was also the highest at 48.5 cells mL⁻¹. When the equation of the growth rates of Y. yeosuensis on Pyramimonas sp. as a function of prey concentration given by Jang et al. (2017a) was used, the calculated growth rate of Y. yeosuensis on Pyramimonas sp. at the prey concentration of 36.7 ng C mL⁻¹ was 0.49 d⁻¹. This calculated growth rate is 37% of the maximum growth rate of Y. yeosuensis on Pyramimonas sp. Therefore, the abundance of Pyramimonas sp. determined in this study can support moderate positive growth of Y. yeosuensis. Meanwhile, Y. yeosuensis is known to have an anti-predation behavior, fast-backward jumping, against common heterotrophic protists (Jeong et al. 2018). Thus, in the population dynamics of Y. yeosuensis, its abundance may be mainly affected by prey availability rather than mortality due to predation. In addition, Y. yeosuensis is known to encyst under prey-depletion conditions (Jang et al. 2019). Therefore, to predict the population dynamics of Y. yeosuensis in natural environments, prey abundance should be determined.

This study is the first report regarding the abundance of Pyramimonas sp., belonging to the subgenus Trichocystis, and T. amphioxeia in Korean waters. It is challenging to distinguish Pyramimonas sp. and T. amphioxeia from other species in each genus through light microscopy because their cells are approximately 5–10 μm in length. However, the abundances of these species were quantified in this study by a qPCR method employing species-specific primers and probes. In the global oceans, cells of Pyramimonas spp. have been found at temperatures and salinities of −1.5 to 28°C and 3.8 to 34.7, respectively (Van Valkenburg et al. 1978, Parry et al. 1989, Domingos and Menezes 1998, Andreoli et al. 2000, Bergesch et al. 2008, Stonik 2009, Tas and Yilmaz 2015), and cells of T. amphioxeia have been found at temperatures and salinities of 2.9–25°C and 0.5–34.8, respectively (Schmidt et al. 1998, Hajdu et al. 2007, Johnson et al. 2016, Herfort et al. 2017, Kim et al. 2018). This study extends the ranges of water temperature and salinity in which cells of T. amphioxeia can exist to 1.1–28.0°C and 0.1–35.6, respectively. Thus, there is a distinct possibility that Y. yeosuensis encounters Pyramimonas sp. and T. amphioxeia in the global oceans. However, the threshold Pyramimonas sp. concentration for the positive growth of Y. yeosuensis is 3.18 ng C mL⁻¹ (79.5 cells mL⁻¹). The maximum growth and ingestion rates of Y. yeosuensis for T. amphioxeia were similar to those for Pyramimonas sp. (Jang et al. 2017a). This indicated that the threshold T. amphioxeia concentration for the positive growth of Y. yeosuensis is likely to be similar to the threshold Pyramimonas sp. concentration. In our study, there were three water samples in which the biomass of Pyramimonas sp. exceeded 3.18 ng C mL⁻¹ and three water samples in which the biomass of T. amphioxeia exceeded 3.18 ng C mL⁻¹. In addition, there were five water samples in which the combined biomass of Pyramimonas sp. and T. amphioxeia exceeded 3.18 ng C mL⁻¹. Thus, Y. yeosuensis may have difficulty in encountering Pyramimonas sp. or T. amphioxeia with a biomass exceeding 3.18 ng C mL⁻¹ in Korean waters. There are yet to be studies on red tides or blooms due to Y. yeosuensis. Low biomasses of the only prey species of Y. yeosuensis, Pyramimonas sp. and T. amphioxeia, are likely to limit the population of Y. yeosuensis.

Our molecular analysis showed that the strain of Pyramimonas sp., PSSH1204, used in this study belongs to the subgenus Trichocystis. The maximum abundances of Pyramimonas spp. belonging to the subgenus Trichocystis in the brackish Patterson Lakes of Australia (300,000 cells mL⁻¹) (McFadden et al. 1986) and the estuary of the Marmara Sea, Turkey (1,600 cells mL⁻¹) (Tas and Yilmaz 2015) are greater than that of Pyramimonas sp. recorded in this study (918 cells mL⁻¹). Thus, if Y. yeosuensis is present in these waters, its maximum abundance may be higher than that found in the present study.

The maximum grazing coefficients (g) of Y. yeosuensis on the co-occurring Pyramimonas sp. was 0.42 d⁻¹ (i.e., up to 35% of the Pyramimonas sp. population was consumed in 1 d). Therefore, Y. yeosuensis may sometimes have considerable grazing impacts on populations of Pyramimonas sp. Thus, we conclude that the spatio-temporal distributions of Y. yeosuensis in the Korean coastal waters may affect and also be affected by those of its optimal prey species.

SUPPLEMENTARY MATERIALS

Supplementary Table S1

List of species used to determine primer specificity, and quantitative real-time PCR (qPCR) results (https://e-algae.org).

algae-2020-35-1-45-suppl1.pdf

Supplementary Table S2

Abundance (cells mL⁻¹) of Yihiella yeosuensis (Yy), Teleaulax amphioxeia (Ta), and Pyramimonas sp. (Pyra) at each sampling station from April 2015 to October 2018 (https://e-algae.org).

algae-2020-35-1-45-suppl2.xlsx

ACKNOWLEDGEMENTS

We thank An Suk Lim, Sung Yeon Lee, Jin Hee Ok, Hee Chang Kang, Ji Hyun You, Ji Eun Kwon, and So Jin Kim for technical supports. This work was supported by the Useful Dinoflagellate Program of Korea Institute of Marine Science and Technology Promotion (KIMST) funded by the Ministry of Oceans and Fisheries (MOF) and the National Research Foundation (NRF) funded by the Ministry of Science and ICT (NRF-2015M1A5A1041806; NRF-2017R1E1A1A01074419) award to HJJ.

Fig. 1

Map showing the sampling stations in Korean coastal waters during the study period. SC, Sokcho; JMJ, Jumunjin; DH, Donghae; UJ, Uljin; PH, Pohang; US, Ulsan; BS, Busan; DDP, Dadaepo; MS, Masan; JH, Jinhae; TY, Tongyoung; YS, Yeosu; KY, Kwangyang; GH, Goheung; JAH, Jangheung; AS, Ansan; DAJ, Dangjin; MGP, Mageompo; TA, Taean; SCN, Seocheon; KS, Kunsan; BA, Buan; MP, Mokpo; AW, Aewol; SGP, Seogwipo; WM, Wimi; SS, Seongsan; GN, Gimnyeong.

Fig. 2

Map of the sampling stations, indicating the presence and absence of Yihiella yeosuensis in spring (March or April) (A), summer (June or July) (B), autumn (September or October) (C), and winter (December or January) (D) from 2015 to 2018. The green closed circles indicate the stations at which Y. yeosuensis cells were detected, whereas the black open circles indicate the stations at which Y. yeosuensis cells were not detected. SC, Sokcho; JMJ, Jumunjin; DH, Donghae; UJ, Uljin; PH, Pohang; US, Ulsan; BS, Busan; DDP, Dadaepo; MS, Masan; JH, Jinhae; TY, Tongyoung; YS, Yeosu; KY, Kwangyang; GH, Goheung; JAH, Jangheung; AS, Ansan; DAJ, Dangjin; MGP, Mageompo; TA, Taean; SCN, Seocheon; KS, Kunsan; BA, Buan; MP, Mokpo; AW, Aewol; SGP, Seogwipo; WM, Wimi; SS, Seongsan; GN, Gimnyeong.

Fig. 3

Map of the sampling stations, indicating the presence and absence of Pyramimonas sp. in spring (March or April) (A), summer (June or July) (B), autumn (September or October) (C), and winter (December or January) (D) from 2015 to 2018. The red closed circles indicate the stations at which Pyramimonas sp. cells were detected, whereas the black open circles indicate the stations at which Pyramimonas sp. cells were not detected. SC, Sokcho; JMJ, Jumunjin; DH, Donghae; UJ, Uljin; PH, Pohang; US, Ulsan; BS, Busan; DDP, Dadaepo; MS, Masan; JH, Jinhae; TY, Tongyoung; YS, Yeosu; KY, Kwangyang; GH, Goheung; JAH, Jangheung; AS, Ansan; DAJ, Dangjin; MGP, Mageompo; TA, Taean; SCN, Seocheon; KS, Kunsan; BA, Buan; MP, Mokpo; AW, Aewol; SGP, Seogwipo; WM, Wimi; SS, Seongsan; GN, Gimnyeong.

Fig. 4

Map of the sampling stations, indicating the presence and absence of Teleaulax amphioxeia in spring (March or April) (A), summer (June or July) (B), autumn (September or October) (C), and winter (December or January) (D) from 2015 to 2018. The blue closed circles indicate the stations at which T. amphioxeia cells were detected, whereas the black open circles indicate the stations at which T. amphioxeia cells were not detected. SC, Sokcho; JMJ, Jumunjin; DH, Donghae; UJ, Uljin; PH, Pohang; US, Ulsan; BS, Busan; DDP, Dadaepo; MS, Masan; JH, Jinhae; TY, Tongyoung; YS, Yeosu; KY, Kwangyang; GH, Goheung; JAH, Jangheung; AS, Ansan; DAJ, Dangjin; MGP, Mageompo; TA, Taean; SCN, Seocheon; KS, Kunsan; BA, Buan; MP, Mokpo; AW, Aewol; SGP, Seogwipo; WM, Wimi; SS, Seongsan; GN, Gimnyeong.

Fig. 5

The abundances (cells mL⁻¹) of Yihiella yeosuensis (Yy) (A), Pyramimonas sp. (Pyra) (B), and Teleaulax amphioxeia (Ta) (C) as a function of water temperature (T) and salinity (S) at all stations from 2015–2018 and those of Yy (D), Pyra (E), and Ta (F) as a function of NO₃ and PO₄ concentrations. The scale of the circles is the abundance of each species (cells mL⁻¹).

Fig. 6

(A) The abundance (cells mL⁻¹) of Yihiella yeosuensis as a function of the abundances of Pyramimonas sp. and Teleaulax amphioxeia at all the stations from 2015–2018. (B–D) The correlations between the abundance of Y. yeosuensis and those of Pyramimonas sp. (B), T. amphioxeia (C), and T. amphioxeia + Pyramimonas sp. (D) (n = 45). The equations of the curves are as follows: (B) Y. yeosuensis (cells mL⁻¹) = 0.052 (Pyramimonas sp.; cells mL⁻¹) + 0.065, r² = 0.937; (C) Y. yeosuensis (cells mL⁻¹) = 0.020 (T. amphioxeia; cells mL⁻¹) + 0.112, r² = 0.124; (D) Y. yeosuensis (cells mL⁻¹) = 0.029 (T. amphioxeia + Pyramimonas sp.; cells mL⁻¹) - 0.048, r² = 0.696. The scale of the circles in (A) is cells mL⁻¹.

Fig. 7

The calculated grazing coefficients (g, n = 31) of Yihiella yeosuensis in relation to the concentration of co-occurring Pyramimonas sp. (see text for calculation). Clearance rates, measured under the conditions indicated by Jang et al. (2017a), were corrected using Q₁₀ = 2.8 (Hansen et al. 1997) because in situ water temperatures obtained in this study and the temperature used in the laboratory by Jang et al. (2017a) (20°C) were occasionally different. The scale of the circles is g (d⁻¹).

Table 1

Sequences of the primer sets and probes for Yihiella yeosuensis (Yy), Pyramimonas sp. (Pyra), and Teleaulax amphioxeia (Ta) used in this study

Target species	Target gene	Analysis	Primer name		Primer sequence (5′-3′)	Reference
Yy	ITS rDNA	PCR	ITSF2	Forward	TACGTCCCTGCCCTTTGTAC	Litaker et al. (2003)
			LSU500R	Reverse	CCCTCATGGTACTTGTTTGC	Litaker et al. (2003)
		qPCR	Yyeosuensis_F	Forward	TCTCATGCAGGTCTTTTCTCAA	This study
			Yyeosuensis_R	Reverse	CCAGTCACAAGTTGCAGATGA	This study
			Yyeosuensis_P	Probe	TTGTCTGTGTCTTTGTGTGCCATAGTGC	This study
Pyra	rbcL	PCR	RH-1S	Forward	ATGTCACCACAAACAGAAACT	Daugbjerg et al. (1994)
			Ce1161R	Reverse	CATGTGCAATACGTGAATACC	Daugbjerg et al. (1994)
		qPCR	Pyramimonas_F	Forward	GCTTCAAGGCTCTACGTGCT	This study
			Pyramimonas_R	Reverse	GACGGTCACGCTCTACCTG	This study
			Pyramimonas_P	Probe	AAGACTTTCCTTGGTCCACCACATGGT	This study
Ta	LSU rDNA	PCR	Cryptophyte_F	Forward	CYTGCTTGGGAATGCAGGYC	This study
			LSUB	Reverse	ACGAACGATTTGCACGTCAG	Litaker et al. (2003)
		qPCR	Tamphioxeia_F	Forward	GGACGTGAGAAACCCTGTTG	This study
			Tamphioxeia_R	Reverse	CTGCGTCAGAAACACACGAG	This study
			Tamphioxeia_P	Probe	GGGCATGCGAGTGAATGATCCATTG	This study

ITS, internally transcribed spacer; rbcL, large subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase; LSU rDNA, large subunit ribosomal DNA;

PCR, polymerase chain reaction; qPCR, quantitative real-time PCR.

Table 2

The ranges of the measures of abiotic environmental factors during the study period and when each of Yihiella yeosuensis, Pyramimonas sp., and Teleaulax amphioxeia was detected in the Korean coastal waters from April 2015 to October 2018

Species	Temperature (°C)	Salinity	NO₃ (μM)	PO₄ (μM)	SiO₂ (μM)	DO (mg L⁻¹)
In this study period	0.2–28.0	0.1–35.6	ND-149.0	ND-6.3	ND-453.4	0.2–14.8
Yihiella yeosuensis	2.3–28.0	9.9–34.4	0.8–91.9	0.1–2.7	2.1–316.6	1.7–11.3
Pyramimonas sp.	3.5–27.0	9.9–34.6	ND-149.0	ND-6.3	ND-453.4	0.4–12.5
Teleaulax amphioxeia	1.1–28.0	0.1–35.6	ND-118.3	ND-6.3	ND-448.4	0.2–14.8

NO₃, concentrations of nitrate + nitrite; PO₄, phosphate; SiO₂, silicate; DO, dissolved oxygen; ND, not detectable.

Table 3

Pearson’s correlation coefficients for the relationships among the abundances of Yihiella yeosuensis (cells mL⁻¹) and its prey species and the relationships between each of these abundances and the abiotic factors during the study period

Species	Biotic environmental factors			Abiotic environmental factors

	Pyra	Ta	Pyra + Ta	T	S	NO₃	PO₄	SiO₂	DO
Y. yeosuensis	0.968**	0.352**	0.834**	0.100*	−0.039	0.008	0.045	0.014	−0.070
Pyramimonas sp.	-	-	-	0.096	−0.032	0.008	0.053	0.014	−0.048
T. amphioxeia	-	-	-	0.135**	−0.107*	0.112*	0.017	0.182**	−0.049

Pyra, Pyramimonas sp. (cells mL⁻¹); Ta, Teleaulax amphioxeia (cells mL⁻¹); T, temperature (°C); S, salinity; NO₃, concentrations of nitrate + nitrite (μM); PO₄, phosphate (μM); SiO₂, silicate (μM); DO, dissolved oxygen (mg L⁻¹).

^* p < 0.05,

^** p < 0.01 (n = 390).

REFERENCES

Adolf, JE, Bachvaroff, T & Place, AR 2008. Can cryptophyte abundance trigger toxic Karlodinium veneficum blooms in eutrophic estuaries? Harmful Algae. 8:119–128.

Andreoli, C, Moro, I, La Rocca, N, Valle, LD, Masiero, L, Rascio, N & Vecchia, FD 2000. Ecological, physiological, and biomolecular surveys on microalgae from Ross Sea (Antarctica). Ital J Zool. 67(S1):147–156.

Berge, T, Hansen, PJ & Moestrup, Ø 2008. Feeding mechanism, prey specificity and growth in light and dark of the plastidic dinoflagellate Karlodinium armiger . Aquat Microb Ecol. 50:279–288.

Bergesch, M, Odebrecht, C & Moestrup,, Ø 2008. Nanoflagellates from coastal waters of southern Brazil (32°S). Bot Mar. 51:35–50.

Burkholder, JM, Glasgow, HB Jr & Hobbs, CW 1995. Fish kills linked to a toxic ambush-predator dinoflagellate: distribution and environmental conditions. Mar Ecol Prog Ser. 124:43–61.

Coats, DW 1999. Parasitic life styles of marine dinoflagellates. J Eukaryot Microbiol. 46:402–409.

Daugbjerg, N, Moestrup, Ø & Arctander, P 1994. Phylogeny of the genus Pyramimonas (Prasinophyceae, Chlorophyta) inferred from the rbc L gene. J Phycol. 30:991–999.

Domingos, P & Menezes, M 1998. Taxonomic remarks on planktonic phytoflagellates in a hypertrophic tropical lagoon (Brazil). Hydrobiologia. 369:297–313.

Egge, JK & Aksnes, DL 1992. Silicate as regulating nutrient in phytoplankton competition. Mar Ecol Prog Ser. 83:281–289.

Eppley, RW 1972. Temperature and phytoplankton growth in the sea. Fish Bull. 70:1063–1085.

Glasgow, HB Jr, Burkholder, JM, Schmechel, DE, Tester, PA & Rublee, PA 1995. Insidious effects of a toxic estuarine dinoflagellate on fish survival and human health. J Toxicol Environ Health. 46:501–522.

Glibert, PM, Burkholder, JM, Kana, TM, Alexander, J, Skelton, H & Shilling, C 2009. Grazing by Karenia brevis on Synechococcus enhances its growth rate and may help to sustain blooms. Aquat Microb Ecol. 55:17–30.

Hajdu, S, Höglander, H & Larsson, U 2007. Phytoplankton vertical distributions and composition in Baltic Sea cyanobacterial blooms. Harmful Algae. 6:189–205.

Hansen, PJ 1992. Prey size selection, feeding rates and growth dynamics of heterotrophic dinoflagellates with special emphasis on Gyrodinium spirale . Mar Biol. 114:327–334.

Hansen, PJ 2011. The role of photosynthesis and food uptake for the growth of marine mixotrophic dinoflagellates. J Eukaryot Microbiol. 58:203–214.

Hansen, PJ, Bjørnsen, PK & Hansen, BW 1997. Zooplankton grazing and growth: scaling within the 2–2,000-μm body size range. Limnol Oceanogr. 42:687–704.

Herfort, L, Maxey, K, Voorhees, I, Simon, HM, Grobler, K, Peterson, TD & Zuber, P 2017. Use of highly specific molecular markers reveals positive correlation between abundances of Mesodinium cf. major and its preferred prey, Teleaulax amphioxeia, during red water blooms in the Columbia river estuary. J Eukaryot Microbiol. 64:740–755.

Jang, SH, Jeong, HJ & Chon, JK 2019. De novo transcriptome of the newly described phototrophic dinoflagellate Yihiella yeosuensis: comparison between vegetative cells and cysts. Mar Biol. 166:104 pp.

Jang, SH, Jeong, HJ, Kwon, JE & Lee, KH 2017a. Mixotrophy in the newly described dinoflagellate Yihiella yeosuensis: a small, fast dinoflagellate predator that grows mixotrophically, but not autotrophically. Harmful Algae. 62:94–103.

Jang, SH, Jeong, HJ, Moestrup, Ø, Kang, NS, Lee, SY, Lee, KH & Seong, KA 2017. b. Yihiella yeosuensis gen. et sp. nov. (Suessiaceae, Dinophyceae), a novel dinoflagellate isolated from the coastal waters of Korea. J Phycol. 53:131–145.

Jeong, HJ 1999. The ecological roles of heterotrophic dinoflagellates in marine planktonic community. J Eukaryotic Microbiol. 46:390–396.

Jeong, HJ, Ha, JH, Park, JY, Kim, JH, Kang, NS, Kim, S, Kim, JS, Yoo, YD & Yih, WH 2006. Distribution of the heterotrophic dinoflagellate Pfiesteria piscicida in Korean waters and its consumption of mixotrophic dinoflagellates, raphidophytes and fish blood cells. Aquat Microb Ecol. 44:263–278.

Jeong, HJ, Kang, HC, You, JH & Jang, SH 2018. Interactions between the newly described small- and fast- swimming mixotrophic dinoflagellate Yihiella yeosuensis and common heterotrophic protists. J Eukaryotic Microbiol. 65:612–626.

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, Yih, 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, Ok, JH, Lim, AS, Kwon, JE, Kim, SJ & Lee, SY 2016. Mixotrophy in the phototrophic dinoflagellate Takayama helix (family Kareniaceae): predator of diverse toxic and harmful dinoflagellates. Harmful Algae. 60:92–106.

Jeong, HJ, Shim, JH, Kim, JS, Park, JY, Lee, CW & Lee, Y 1999. Feeding by the mixotrophic thecate dinoflagellate Fragilidium cf. mexicanum on red-tide and toxic dinoflagellates. Mar Ecol Prog Ser. 176:263–277.

Jeong, HJ, Yoo, YD, Kang, NS, Lim, AS, Seong, KA, Lee, SY, Lee, MJ, Lee, KH, Kim, HS, Shin, W, Nam, SW, Yih, W & Lee, K 2012. Heterotrophic feeding as a newly identified survival strategy of the dinoflagellate Symbiodinium . Proc Natl Acad Sci U S A. 109:12604–12609.

Jeong, HJ, Yoo, YD, Kang, NS, Rho, JR, Seong, KA, Park, JW, Nam, GS & Yih, W 2010. a. Ecology of Gymnodinium aureolum. I. Feeding in western Korean waters. Aquat Microb Ecol. 59:239–255.

Jeong, HJ, Yoo, YD, Kim, JS, Kim, TH, Kim, JH, Kang, NS & Yih, W 2004. Mixotrophy in the phototrophic harmful alga Cochlodinium polykrikoides (Dinophycean): prey species, the effects of prey concentration, and grazing impact. J Eukaryotic Microbiol. 51:563–569.

Jeong, HJ, Yoo, YD, Kim, JS, Seong, KA, Kang, NS & Kim, TH 2010. b. Growth, feeding and ecological roles of the mixotrophic and heterotrophic dinoflagellates in marine planktonic food webs. Ocean Sci J. 45:65–91.

Jeong, HJ, Yoo, YD, Lee, KH, Kim, TH, Seong, KA, Kang, NS, Lee, SY, Kim, JS, Kim, S & Yih, WH 2013. Red tides in Masan Bay, Korea in 2004–2005: I. Daily variations in the abundance of red-tide organisms and environmental factors. Harmful Algae. 30(Suppl 1):S75–S88.

Jeong, HJ, Yoo, YD, Park, JY, Song, JY, Kim, ST, Lee, SH, Kim, KY & Yih, WH 2005. a. Feeding by phototrophic red-tide dinoflagellates: five species newly revealed and six species previously known to be mixotrophic. Aquat Microb Ecol. 40:133–150.

Jeong, HJ, Yoo, YD, Seong, KA, Kim, JH, Park, JY, Kim, S, Lee, SH, Ha, JH & Yih, WH 2005. b. Feeding by the mixotrophic red-tide dinoflagellate Gonyaulax polygramma: mechanisms, prey species, effects of prey concentration, and grazing impact. Aquat Microb Ecol. 38:249–257.

Johnson, MD, Beaudoin, DJ, Laza-Martinez, A, Dyhrman, ST, Fensin, E, Lin, S, Merculief, A, Nagai, S, Pompeu, M, Setälä, O & Stoecker, DK 2016. The genetic diversity of Mesodinium and associated cryptophytes. Front Microbiol. 7:2017.

Kang, HC, Jeong, HJ, Jang, SH & Lee, KH 2019a. Feeding by common heterotrophic protists on the phototrophic dinoflagellate Biecheleriopsis adriatica (Suessiaceae) compared to that of other suessioid dinoflagellates. Algae. 34:127–140.

Kang, HC, Jeong, HJ, Ok, JH, You, JH, Jang, SH, Lee, SY, Lee, KH, Park, JY & Rho, J-R 2019b. Spatial and seasonal distributions of the phototrophic dinoflagellate Biecheleriopsis adriatica (Suessiaceae) in Korea: quantification using qPCR. Algae. 34:111–126.

Kang, NS, Jeong, HJ, Yoo, YD, Yoon, EY, Lee, KH, Lee, K & Kim, G 2011. Mixotrophy in the newly described phototrophic dinoflagellate Woloszynskia cincta from western Korean waters: feeding mechanism, prey species and effect of prey concentration. J Eukaryotic Microbiol. 58:152–170.

Kibbe, WA 2007. OligoCalc: an online oligonucleotide properties calculator. Nucleic Acids Res. 35(Suppl 2):W43–W46.

Kim, D, Matsuyama, Y, Nagasoe, S, Yamaguchi, M, Yoon, Y, Oshima, Y, Imada, N & Honjo, T 2004. Effects of temperature, salinity and irradiance on the growth of the harmful red tide dinoflagellate Cochlodinium polykrikoides Margalef (Dinophyceae). J Plankton Res. 26:61–66.

Kim, M, Ji, H, Cho, J & Cho, S 2018. Identification of red tide-causing organism and characteristics of red tide occurrence in the Oncheon Down Stream, Busan. J Korean Soc Water Environ. 34:285–292, (in Korean).

Kudela, RM & Gobler, CJ 2012. Harmful dinoflagellate blooms caused by Cochlodinium sp.: global expansion and ecological strategies facilitating bloom formation. Harmful Algae. 14:71–86.

Lee, KH, Jeong, HJ, Jang, TY, Lim, AS, Kang, NS, Kim, J, Kim, KY, Park, K & Lee, K 2014. a. Feeding by the newly described mixotrophic dinoflagellate Gymnodinium smaydae: feeding mechanism, prey species, and effect of prey concentration. J Exp Mar Biol Ecol. 459:114–125.

Lee, KH, Jeong, HJ, Kang, HC, Ok, JH, You, JH & Park, SA 2019a. 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, Kim, HJ & Lim, AS 2017a. Nitrate uptake of the red tide dinoflagellate Prorocentrum micans measured using a nutrient repletion method: effect of light intensity. Algae. 32:139–153.

Lee, KH, Jeong, HJ, Kwon, JE, Kang, HC, Kim, JH, Jang, SH, Park, JY, Yoon, EY & Kim, JS 2016a. Mixotrophic ability of the phototrophic dinoflagellates Alexandrium andersonii, A. affine, and A. fraterculus . Harmful Algae. 59:67–81.

Lee, MJ, Jeong, HJ, Jang, SH, Lee, SY, Kang, NS, Lee, KH, Kim, HS, Wham, DC & LaJeunesse, TC 2016. b. Most low-abundance “background” Symbiodinium spp. are transitory and have minimal functional significance for symbiotic corals. Microb Ecol. 71:771–783.

Lee, MJ, Jeong, HJ, Lee, KH, Jang, SH, Kim, JH & Kim, KY 2015. Mixotrophy in the nematocyst-taeniocyst complex-bearing phototrophic dinoflagellate Polykrikos hartmannii . Harmful Algae. 49:124–134.

Lee, SK, Jeong, HJ, Jang, SH, Lee, KH, Kang, NS, Lee, MJ & Potvin, É 2014b. Mixotrophy in the newly described dinoflagellate Ansanella granifera: feeding mechanism, prey species, and effect of prey concentration. Algae. 29:137–152.

Lee, SY, Jeong, HJ, Kwon, JE, You, JH, Kim, SJ, Ok, JH, Kang, HC & Park, JY 2019b. First report of the photosynthetic dinoflagellate Heterocapsa minima in the Pacific Ocean: morphological and genetic characterizations and the nationwide distribution in Korea. Algae. 34:7–21.

Lee, SY, Jeong, HJ, Seong, KA, Lim, AS, Kim, JH, Lee, KH, Lee, MJ & Jang, SH 2017b. Improved real-time PCR method for quantification of the abundance of all known ribotypes of the ichthyotoxic dinoflagellate Cochlodinium polykrikoides by comparing 4 different preparation methods. Harmful Algae. 63:23–31.

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.

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.

McFadden, GI, Hill, DRA & Wetherbee, R 1986. A study of the genus Pyramimonas (Prasinophyceae) from southeastern Australia. Nord J Bot. 6:209–234.

Park, J, Jeong, HJ, Yoo, YD & Yoon, EY 2013a. Mixotrophic dinoflagellate red tides in Korean waters: distribution and ecophysiology. Harmful Algae. 30(Suppl 1):S28–S40.

Park, TG, Lim, WA, Park, YT, Lee, CK & Jeong, HJ 2013b. Economic impact, management and mitigation of red tides in Korea. Harmful Algae. 30(Suppl 1):S131–S143.

Parry, GD, Langdon, JS & Huisman, JM 1989. Toxic effects of a bloom of the diatom Rhizosolenia chunii on shellfish in Port Phillip Bay, southeastern Australia. Mar Biol. 102:25–41.

Schmidt, K, Kähler, P & von Bodungen, B 1998. Copepod egg production rates in the Pomeranian Bay (Southern Baltic Sea) as a function of phytoplankton abundance and taxonomic composition. Mar Ecol Prog Ser. 174:183–195.

Sherr, EB & Sherr, BF 2007. Heterotrophic dinoflagellates: a significant component of microzooplankton biomass and major grazers of diatoms in the sea. Mar Ecol Prog Ser. 352:187–197.

Shumway, SE, Burkholder, JM & Morton, SL 2018. Harmful algal blooms: a compendium desk reference. John Wiley & Sons, Hoboken, NJ, 667 pp.

Smith, KF, Rhodes, L, Harwood, DT, Adamson, J, Moisan, C, Munday, R & Tillmann, U 2016. Detection of Azadinium poporum in New Zealand: the use of molecular tools to assist with species isolations. J Appl Phycol. 28:1125–1132.

Stoecker, DK, Hansen, PJ, Caron, DA & Mitra, A 2017. Mixotrophy in the marine plankton. Annu Rev Mar Sci. 9:311–335.

Stoecker, DK, Li, A, Coats, DW, Gustafson, DE & Nannen, MK 1997. Mixotrophy in the dinoflagellate Prorocentrum minimum . Mar Ecol Prog Ser. 152:1–12.

Stonik, IV 2009. Species composition and seasonal dynamics of the population density and biomass of the genus Pyramimonas (Prasinophyceae) from the Russian waters of East/Japan Sea. Ocean Sci J. 44:21–26.

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.

Tas, S & Yilmaz, IN 2015. Potentially harmful microalgae and algal blooms in a eutrophic estuary in Turkey. Mediterr Mar Sci. 16:432–443.

Taylor, FJR, Hoppenrath, M & Saldarriaga, JF 2008. Dinoflagellate diversity and distribution. Biodivers Conserv. 17:407–418.

Van Valkenburg, SD, Jones, JK & Heinle, DR 1978. A comparison by size class and volume of detritus versus phytoplankton in Chesapeake Bay. Estuar Coast Mar Sci. 6:569–582.

Yih, W, Kim, HS, Jeong, HJ, Myung, G & Kim, YG 2004. Ingestion of cryptophyte cells by the marine photosynthetic ciliate Mesodinium rubrum . Aquat Microb Ecol. 36:165–170.

Yoo, YD, Jeong, HJ, Kang, NS, Song, JY, Kim, KY, Lee, G & Kim, J 2010. Feeding by the newly described mixotrophic dinoflagellate Paragymnodinium shiwhaense: feeding mechanism, prey species, and effect of prey concentration. J Eukaryot Microbiol. 57:145–158.

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.