Feeding by common heterotrophic protists on the phototrophic dinoflagellate Biecheleriopsis adriatica (Suessiaceae) compared to that of other suessioid dinoflagellates
Article information
Abstract
The species in the dinoflagellate order Suessiales have 5–24 latitudinal paraplate series and include many fossil and extant species. There have been a few studies on the ecophysiology of the phototrophic species Biecheleriopsis adriatica, and no study on its predators. Thus, we explored the feeding occurrence by common heterotrophic protists on B. adriatica and the growth and ingestion rates of the heterotrophic dinoflagellate Oxyrrhis marina on B. adriatica BATY06 as a function of prey concentration. The common heterotrophic dinoflagellates Aduncodinium glandula, O. marina, Gyrodinium dominans, Gyrodinium moestrupii, Luciella masanensis, Pfiesteria piscicida, and Oblea rotunda and two naked ciliates Strombidinopsis sp. and Pelagostrobilidium sp. were able to feed on B. adriatica, but the heterotrophic dinoflagellate Polykrikos kofoidii was not. However, B. adriatica supported the positive growth of O. marina, but did not support that of G. dominans and O. rotunda. With increasing prey concentrations, the growth and ingestion rates of O. marina on B. adriatica increased and became saturated. The maximum growth rate of O. marina on B. adriatica was 0.162 d−1. Furthermore, the maximum ingestion rate of O. marina on B. adriatica was 0.2 ng C predator−1 d−1 (2.0 cells predator−1 d−1). In the order Suessiales, the feeding occurrence by common heterotrophic protists on B. adriatica is similar to that on Effrenium voratum and Biecheleria cincta, but different from that on Yihiella yeosuensis. However, the growth and ingestion rates of O. marina on B. adriatica are considerably lower than those on E. voratum and B. cincta, but higher than those on Y. yeosuensis. Therefore, B. adriatica may be less preferred prey for O. marina than E. voratum and B. cincta, but more preferred prey than Y. yeosuensis.
INTRODUCTION
Dinoflagellates are ubiquitous and a major component of marine ecosystems (Lessard 1984, Jeong 1999, Lim et al. 2017b). They have diverse trophic modes, such as exclusively autotrophic, mixotrophic, and heterotrophic (Stoecker 1999, Jeong et al. 2010). Due to their trophic mode diversity, they can play diverse roles in marine planktonic food webs as primary producers, predators on diverse prey items (including bacteria, microalgae, and metazoans), and prey for heterotrophic protists and metazoans (Coats 1999, Tillmann 2004, Hansen 2011, Johnson 2015, Stoecker et al. 2017). Therefore, dinoflagellates are largely involved in the cycling of materials and energy flow in marine ecosystems (Eppley et al. 1973, Carlsson et al. 1999, Calbet and Landry 2004). Many dinoflagellates species are known to form red tides or harmful algal blooms (HABs), which often cause the large-scale mortality of marine organisms and subsequently great economic loss in diverse industries (Adolf et al. 2015, Menden-Deuer and Montalbano 2015). Thus, to preserve marine organisms and also reduce the economic losses due to red tides or HABs caused by a dinoflagellate, the population dynamics of the dinoflagellate should be well understood. In population dynamics models, the growth and mortality rates of the dinoflagellate are two critical parameters (Jeong et al. 2015). Thus, to understand the roles of a dinoflagellate in marine ecosystems, interactions between the dinoflagellate species and other related organisms (as a potential prey or predator) should be well documented.
Species in the dinoflagellate order Suessiales are known as dinoflagellates, whose cell surface is covered with many thin plates (Lindberg et al. 2005). There are many species in this order (Algaebase, http://www.algaebase.org). However, compared to species in other orders, such as Peridiniales, Gonyaulacales, and Gymnodiniales, there have been fewer species whose predators and mortality rates are reported (Jeong et al. 1997, 2010, Kang et al. 2018). In the order Suessiales, the predators and mortality rates of Effrenium (Symbiodinium) voratum, Biecheleria cincta, Protodinium simplex, and Yihiella yeosuensis have been reported (Strom and Morello 1998, Montagnes and Lessard 1999, Yoo et al. 2013c, Jeong et al. 2014, 2018a). Biecheleriopsis adriatica was reported as a new genus and species of the order Suessiales in 2009 (Moestrup et al. 2009). This species is characterized from other groups of woloszynskioid dinoflagellates by having a nuclear fibrous connective, and a distinct tongue-like process in the left ventral corner of the asymmetric hyposome (Moestrup et al. 2009, Jang et al. 2015). The vegetative cells of B. adriatica were widely distributed in Korean waters (Kang et al. 2019). To understand its ecological roles in marine food webs, its predators and mortality rates should be explored.
Heterotrophic protists are also a major component of marine ecosystems (Sherr and Sherr 2002, Jeong et al. 2010, Lim et al. 2017b). Their grazing impact on phytoplankton populations are usually greater than that of metazoans because the abundance of heterotrophic protist predators is much greater than that of metazoan predators (Turner and Borkman 2005, Lim et al. 2017b). Heterotrophic dinoflagellates (HTDs) and ciliates are major heterotrophic protist groups (Jeong et al. 1999, Levinsen and Nielsen 2002). They are often effective predators of phototrophic dinoflagellates and sometimes control prey populations (Jeong et al. 2003, Lim et al. 2017a). Aduncodinium glandula, Gyrodinium dominans, Gyrodinium moestrupii, Luciella masanensis, Oblea rotunda, Oxyrrhis marina, Polykrikos kofoidii, and Pfiesteria piscicida are common HTDs and Pelagostrobilidium sp. and Strombidinopsis sp. are also common naked ciliates in many marine environments (Strom and Buskey 1993, Claessens et al. 2008, Taylor et al. 2008, Watts et al. 2010, Calbet et al. 2013, Tillmann and Hoppenrath 2013, Yoo et al. 2013b). There are usually large variations in feeding occurrence and growth and ingestion rates of heterotrophic protist predator species when diverse prey species are provided (Hansen 1992, Menden-Deuer et al. 2005, Jeong et al. 2018a, 2018b, Kang et al. 2018). Thus, it is worthwhile to explore the interactions between B. adriatica and these potential heterotrophic protist predators.
In the present study, we investigated the types of predators that are able to feed on a Korean strain of B. adriatica. Furthermore, the growth and ingestion rates of O. marina on B. adriatica were measured as a function of prey concentration. The growth and ingestion rates were compared to the rates of the same predator on different species in the order Suessiales. The results of this study provide a basis for understanding the interactions between B. adriatica and common heterotrophic protists, and their ecological roles in the marine planktonic community.
MATERIALS AND METHODS
Preparation of experimental organisms
Cells of B. adriatica were isolated from plankton samples that were collected from surface waters off the coast of Tongyoung, Korea using plankton samplers in August 2006, when the water temperature and salinity were 28.0°C and 31.0, respectively (Table 1). The collected samples were screened softly by using a 154-μm Nitex mesh. The clonal culture of B. adriatica BATY06 was established using two consecutive single-cell isolations (Jang et al. 2015). When the B. adriatica concentration had increased sufficiently, the volume of the bottles increased to 32, 270, and 500-mL PC bottles containing fresh f/2 medium. The bottles were placed on a shelf at 20°C, illuminated with an irradiance of 20 μmol photons m−2 s−1 provided by cool white fluorescent lights, under a 14 : 10 h light : dark cycle. Only cultures in the exponential growth phase were used.
The HTDs A. glandula, G. dominans, G. moestrupii, L. masanensis, O. rotunda, O. marina, and P. kofoidii, isolated from plankton samples collected from the coastal waters off Masan, Shiwha, Saemankeum, Jinhae, and Jangheung in 2007–2016 were used in this study (Table 1). A clonal culture of each HTD species was established using two consecutive single-cell isolations, except for that of P. piscicida, which was obtained from the National Center for Marine Algae and Microbiota (NCMA), USA. The naked ciliates Pelagostrobilidium sp. and Strombidinopsis sp. were isolated from plankton samples, collected using 20- and 10-μm mesh nets from the coastal waters off Tongyoung and Yeosu in August 2017 and July 2018, respectively (Table 1). A clonal culture for each of Pelagostrobilidium sp. and Strombidinopsis sp. was also established using two serial single-cell isolations.
The carbon content of B. adriatica BATY06 (0.1 ng C per cell) was estimated from the cell volume, according to the equation suggested by Menden-Deuer and Lessard (2000). Furthermore, the carbon contents of all predator species used in the present study were also estimated from the cell volume as described above. The cell volumes of the HTD predators were estimated using the methods of Jang et al. (2016) for A. glandula; Kim and Jeong (2004) and Yoo et al. (2013b) for G. dominans and G. moestrupii, respectively; Jeong et al. (2007) for L. masanensis and P. piscicida; Ok et al. (2017) for O. rotunda; Jeong et al. (2008b) for O. marina; Jeong et al. (2001b) for P. kofoidii; Kim et al. (2019) for Strombidinopsis sp.; and Jeong et al. (2018b) for Pelagostrobilidium sp. (Table 1).
Interactions between Bicheleriopsis adriatica and heterotrophic protists
Experiment (Expt) 1 was designed to investigate feeding by each of the HTDs and ciliates on B. adriatica BATY06, after mixing B. adriatica with potential predator species (Table 2). In this experiment, whether the target heterotrophic protist was able to feed on B. adriatica and/or other interactions were observed.
Dense cultures of B. adriatica (ca. 40,000 cells mL−1) and each of the HTDs and ciliates were added to each 42-mL PC bottle (Table 2). For each experiment, one experiment (with prey and predator), one prey control (without predator), and one predator control (without prey) bottles were set up. The bottles for the A. glandula predator were placed on a shelf, but the bottles for the rest of the predators were placed on a rotating wheel at 0.00017 g (0.9 rpm). All bottles were incubated at 20°C, under an illumination of 20 μmol photons m−2 s−1 and a 14 : 10 h light : dark cycle.
After 2, 24, and 48 h of incubation, 3 mL aliquots were taken from each bottle and transferred into the wells of 6-well plate chambers. More than 30 cells of each predator in one plate chamber were tracked for 2 min under a dissecting microscope at 20–63× magnification to determine whether the predator was able to feed on B. adriatica. The cells of predators having ingested cells of B. adriatica in their body were photographed on slides with cover-glasses at a 200–400× magnification using a digital camera (Zeiss-AxioCam MRc5; Carl Zeiss Ltd., Göttingen, Germany) attached to an inverted light microscope (Zeiss-Axiovert 200 M; Carl Zeiss Ltd.). The feeding process of each of these heterotrophic protists on B. adriatica was recorded using a video analyzing system (Sony DXC-C33; Sony Co., Tokyo, Japan) and captured using the digital camera.
Growth and ingestion rates of Oxyrrhis marina feeding on Biecheleriopsis adriatica as a function of prey concentration
Expt 2 was designed to measure the growth and ingestion rates of O. marina feeding on B. adriatica BATY06 as a function of prey concentration (Table 2). In preliminary tests, B. adriatica BATY06 supported the positive growth of O. marina, but did not support the positive growth of G. dominans, O. rotunda, and Strombidinopsis sp. Whether B. adriatica supports the positive growth of A. glandula, P. piscicida, and L. masanensis was not tested because these dinoflagellates were not distinguishable from B. adriatica when fixed with Lugol’s solution and formalin.
Dense cultures of O. marina grown with Amphidinium carterae as prey were transferred into 250-mL PC bottles when cells of A. carterae were not detectable for 24 h. The bottles were filled to capacity with freshly filtered seawater, capped, placed on a plankton wheel rotating at 0.00017 g (0.9 rpm), and incubated at 20°C under an illumination of 20 μmol photons m−2 s−1 and a 14 : 10 h light : dark cycle. This was done to minimize the possible residual growth resulting from the ingestion of prey during the batch culture. After one day, the cells in three 1 mL aliquots from each bottle were counted using a compound microscope to determine the concentration of predator cells, and the cultures were used in further experiments.
The initial concentrations of O. marina and B. adriatica were established in eight different combinations. Triplicate 42-mL PC experimental bottles (mixtures of predator and prey) and triplicate control bottles (prey only) were set up at each predator–prey combination. In addition, triplicate control bottles (predator only) were established at a single predator concentration. Predetermined volumes of O. marina and B. adriatica were added into each bottle using an autopipette. To obtain similar water conditions, the water of the predator culture was filtered through a 0.7-μm GF/F filter, and then added to the prey control bottles in the same amount as the volume of the predator culture added into the experiment bottles for each predator–prey combination. Similarly, the water of the prey culture was filtered through a 0.7-μm GF/F filter and then added to the predator control bottles in the same amount as the volume of the prey culture added into the experimental bottles. To all the bottles, 5 mL of f/2 medium was added to provide sufficient nutrients to B. adriatica grown autotrophically, and then filled to capacity with freshly filtered seawater and capped. To determine the actual initial predator and prey densities at the beginning of the experiment, a 5 mL aliquot was sampled from each bottle, fixed with 5% Lugol’s solution, and enumerated in three 1-mL Sedgewick Rafter chambers (SRCs). The bottles were then refilled to capacity with freshly filtered seawater, capped, and placed on rotating wheels at 0.00017 g under the conditions described above. Dilution of the cultures from the refilling of bottles was considered when calculating the growth and ingestion rates. A 10 mL aliquot was taken from each bottle at 48 h and fixed with 5% Lugol’s solution, and the abundances of predators and prey were then determined by counting all or >200 cells in three 1-mL SRCs.
The specific growth rate of a heterotrophic protist predator, μ (d−1) was calculated according to following formula:
, where P0 and Pt = the concentrations of the predator at 0 and 2 d.
Data for O. marina growth rates were fitted to a modified Michaelis–Menten equation:
, where μmax = the maximum growth rate (d−1); x = prey concentration (cells mL−1 or ng C mL−1), x′ = threshold of prey concentration (where the prey concentration μ = 0), and KGR = the prey concentration sustaining 1/2 μmax. Data were iteratively fitted to the model using DeltaGraph (Red Rock Software Inc., Salt Lake, UT, USA).
The ingestion rate and mean prey concentration were calculated using the equations of Frost (1972) and Heinbokel (1978), respectively. The incubation time for calculating ingestion was the same as that for estimating growth rate. Data for O. marina ingestion rates (IR, cells predator−1 d−1 or ng C predator−1 d−1) were fitted into a modified Michaelis–Menten equation:
, where Imax = the maximum ingestion rate (cells predator−1 d−1 or ng C predator−1 d−1); x = prey concentration (cells mL−1 or ng C mL−1), and KIR = the prey concentration sustaining 1/2 Imax.
RESULTS
Interactions between Biecheleriopsis adriatica and heterotrophic protists
The HTDs A. glandula, O. marina, G. dominans, G. moestrupii, L. masanensis, P. piscicida, and O. rotunda and the ciliates Strombidinopsis sp. and Pelagostrobilidium sp. were able to feed on B. adriatica BATY06, but the HTD P. kofoidii was not (Table 1, Fig. 1). Cells of O. marina, G. dominans, G. moestrupii, Strombidinopsis sp., and Pelagostrobilidium sp. engulfed a whole B. adriatica cell. Furthermore, L. masanensis and P. piscicida fed on a B. adriatica cell using a peduncle. Moreover, O. rotunda fed on a B. adriatica cell using a pallium after anchoring the prey cell using a tow filament. However, P. kofoidii did not deploy the nematocyst–taeniocyst complex to B. adriatica cells.
O. marina engulfed whole B. adriatica cells (Fig. 2A–F). The time for a B. adriatica cell to be completely engulfed by O. marina was ca. 40–50 s. Meanwhile, A. glandula deployed the peduncle to the body of B. adriatica and then sucked the body materials of B. adriatica (Fig. 2G–L). The A. glandula sometimes left some residual materials of the B. adriatica body. The time for a B. adriatica cell to be ingested by A. glandula was ca. 120–180 s (Fig. 2G–L).
Growth and ingestion rates of Oxyrrhis marina feeding on Biecheleriopsis adriatica as a function of prey concentration
With increasing mean prey concentrations, the specific growth rate of O. marina feeding on B. adriatica BATY06 increased at B. adriatica concentrations <493 ng C mL−1 (4,930 cells mL−1), but became saturated at higher mean prey concentrations (Fig. 3). When the data were fitted to Eq. (2), the calculated maximum growth rate (μmax) of O. marina on B. adriatica was 0.162 d−1.
With increasing mean prey concentrations, the ingestion rate of O. marina feeding on B. adriatica increased at B. adriatica concentrations <493 ng C mL−1 (4,930 cells mL−1), but became almost saturated at higher mean prey concentrations (Fig. 4). The highest ingestion rate was 0.35 ng C predator−1 d−1 (3.5 cells predator−1 d−1). However, when the data were fitted to Eq. (3), the calculated maximum ingestion rate (Imax) of O. marina on B. adriatica was 0.2 ng C predator−1 d−1 (2.0 cells predator−1 d−1).
DISCUSSION
Prior to the present study, among the species belonging to the order Suessiales, only a few have been tested regarding whether common heterotrophic protist predators are able to feed on the prey species or not. The results of the present study clearly showed that among the common heterotrophic protist predators tested, all of the predators except for P. kofoidii were able to feed on B. adriatica BATY06. The minimum prey size that P. kofoidii is able to feed on is known to be approximately 10 μm (Jeong et al. 2001b). The average size of B. adriatica tested in this study was 10.1 μm. Thus, the small size of B. adriatica may be partially responsible for it not been fed upon by P. kofoidii, which deploys the nematocyst–taeniocyst complex to anchor prey cells. The types of common heterotrophic protist predators that are able to feed on B. adriatica were identical or very similar to E. voratum and B. cincta, but largely different from Y. yeosuensis, which only O. marina, A. glandula, and Strombidinopsis sp. were able to feed on (Table 3). Thus, B. adriatica is vulnerable to common heterotrophic protist predators as much as E. voratum and B. cincta; however, it is more vulnerable than Y. yeosuensis. Furthermore, the μmax and Imax of O. marina on B. adriatica was lower than those of E. voratum and B. cincta, but greater than those of Y. yeosuensis (Table 4). Thus, it is expected that in marine environments, O. marina is less abundant when B. adriatica is abundant than when E. voratum or B. cincta is abundant, but more abundant than when Y. yeosuensis prey is abundant. Meanwhile, either the μmax or Imax of O. marina feeding on B. adriatica, B. cincta, E. voratum, and Y. yeosuensis was not significantly related to the size of the prey species (Fig. 5A & B). Thus, factors other than prey size may affect the μmax or Imax of O. marina. The maximum swimming speeds of B. adriatica and Y. yeosuensis were much greater than those of B. cincta and E. voratum (Table 4). Furthermore, B. adriatica and Y. yeosuensis jump backward when a predator attacks, whereas B. cincta and E. voratum do not jump (Kang et al. 2011, Yoo et al. 2013c, Jeong et al. 2014, 2018a, Jang et al. 2015, 2017, personal observation). Thus, the much higher maximum swimming speeds and jumping behaviors of B. adriatica and Y. yeosuensis may be partially responsible for the lower μmax or Imax of O. marina on these prey species compared to those on B. cincta and E. voratum. The μmax of O. marina feeding on B. adriatica, B. cincta, E. voratum, and Y. yeosuensis was also not significantly related to Imax (Fig. 5C). Thus, the gross growth efficiencies (GGEs) of O. marina on these four dinoflagellate prey species are different (Table 4). The GGE of O. marina on B. adriatica (18%) is lower than that on B. cincta (27%), but higher than that on E. voratum (10%) or Y. yeosuensis (0%). Thus, the GGEs of O. marina on these suessioid dinoflagellate prey species are wide.
When the μmax and Imax of O. marina on B. adriatica BATY06 were compared with those on the dinoflagellate species belonging to the other orders, the μmax of O. marina on B. adriatica were higher than that on the HTD L. masanensis and the mixotrophic dinoflagellates Paragymnodinium shiwhaense and Takayama helix, but lower than that on many species belonging to other orders, such as Amphidiniales, Gymnodiniales, Thoracosphaerales, Chlamydomonadales, Isochrysidales, Chattonellales, and Eutreptiales (Table 5). Furthermore, the Imax of O. marina on B. adriatica was higher than that on L. masanensis, Stoeckeria algicida, P. shiwhaense, and T. helix, but lower than that on many species belonging to other orders Amphidiniales, Gymnodiniales, Thoracosphaerales, Chlamydomonadales, Isochrysidales, Chattonellales, and Eutreptiales (Table 5). Therefore, B. adriatica may not be preferred prey for O. marina.
Vegetative cells of B. adriatica was present in the waters of 20 stations when surface water samples were collected from 28 stations along the Korean Peninsula and Jeju Island from April 2015 to October 2018 (Kang et al. 2019). However, the maximum abundance of B. adriatica in that study was 41.7 cells mL−1 (4.2 ng C mL−1), which is lower than the threshold prey concentration for the growth of O. marina feeding on B. adriatica (729 cells mL−1, 72.9 ng C mL−1) obtained in the present study. Furthermore, the highest reported abundance of B. adriatica in the world’s oceans was 305 cells mL−1 (30.5 ng C mL−1), which was obtained from Bolinao in the Philippines (Benico et al. 2019). This abundance is also lower than the threshold prey concentration for the growth of O. marina feeding on B. adriatica. Vegetative cells of B. adriatica have been found in the East, West, and South Seas of Korea (Fig. 6). Cells of O. marina have also been found in the East, West, and South Seas of Korea (Fig. 6). Thus, there is a high possibility that O. marina encounter and feed on B. adriatica in these waters. However, O. marina may not be abundant when B. adriatica co-occurs in Korean waters due to its abundances being lower than the threshold prey concentration for growth. Furthermore, G. dominans and P. kofoidii have been found in Korean waters in all four seasons (Jeong et al. 2002, Kang et al. 2013, Yoo et al. 2013a, Lim et al. 2017b). Thus, there is also a high possibility that these predators encounter and feed on B. adriatica. However, G. dominans and P. kofoidii did not grow on B. adriatica in this study; these common heterotrophic protists may not be abundant as the same time that B. adriatica co-occurs in these waters.
In conclusion, the combination of the results, such as the lower μmax or Imax of O. marina on B. adriatica BATY06 than those on other prey species in the order Suessiales and most phototrophic prey species in the other orders; higher threshold prey concentration for the growth of O. marina feeding on B. adriatica than the highest reported abundance of B. adriatica; and no growth of G. dominans, O. rotunda, P. kofoidii, and Strombidinopsis sp. on B. adriatica, suggests that B. adriatica may have an advantage over other competing prey species regarding their survival in marine environments.
ACKNOWLEDGEMENTS
We thank Sung Yeon Lee, An Suk Lim, and Sang Ah Park for technical support. This research 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 Development of the methods for controlling and managing marine ecological disturbance causing and harmful organisms (MEDHO) funded by the Ministry of Oceans and Fisheries, the National Research Foundation (NRF) funded by the Ministry of Science and ICT (NRF-2015M1A5A1041806; NRF-2017R1E1A1A01074419) award to HJJ.