Seasonal alternation of the ontogenetic development of the moon jellyfish Aurelia coerulea in Maizuru Bay, Japan
Authors:
Kentaro S. Suzuki aff001; Keita W. Suzuki aff002; Emi Kumakura aff003; Kana Sato aff003; Yutaro Oe aff004; Tasuku Sato aff004; Hideki Sawada aff002; Reiji Masuda aff002; Yasuyuki Nogata aff001
Authors place of work:
Environmental Science Research Laboratory, Central Research Institute of Electric Power Industry, Abiko, Chiba, Japan
aff001; Maizuru Fisheries Research Station, Field Science Education and Research Center, Kyoto University, Maizuru, Kyoto, Japan
aff002; CERES Inc., Abiko, Chiba, Japan
aff003; College of Bioresource Sciences, Nihon University, Fujisawa, Kanagawa, Japan
aff004
Published in the journal:
PLoS ONE 14(11)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0225513
Summary
Outbreaks of moon jellyfish Aurelia spp. are frequently reported from many parts of the world’s coastal areas. Aurelia spp. canonically show a metagenetic life cycle in which planulae transform into sessile polyps, which can drastically increase in number through asexual reproduction. Therefore, their asexual reproduction has been recognized as one of the major causes of the outbreaks. Aurelia spp. also show direct development that lacks asexual reproduction during the polyp stage, which prevents us from understanding the mechanisms of its outbreaks. To clarify the seasonality of the metagenetic and direct-development life cycles of Aurelia sp. in Maizuru Bay, Japan, we conducted field observations and laboratory experiments throughout the year. Additionally, the two life cycle types were genetically analyzed to confirm that they belong to the single species Aurelia coerulea, which dominates in coastal waters in Japan. From July until October, Aurelia coerulea produced smaller eggs and planulae all of which developed into polyps. However, from December until May, larger eggs and planulae were produced and 90% of the planulae developed into planktonic ephyrae bypassing the sessile polyp stage. Our results demonstrated that a single species, A. coerulea, seasonally shifts between their two life cycle types at a water temperature threshold of 20°C in Maizuru Bay. The higher energy storage of larger planulae was suggested to enable the planulae to develop into ephyrae without external energy input through feeding during the polyp stage. The adaptive significances of the two life cycle types were also discussed.
Keywords:
Ovaries – Copepods – Predation – Spawning – Life cycles – Salinity – Oocytes – Asexual reproduction
Introduction
Outbreaks of scyphozoan jellyfish have been frequently reported from many parts of the world’s coastal areas [1–3]. Among the scyphozoans, moon jellyfish Aurelia spp. are most common in temperate coastal waters, and outbreaks of Aurelia spp. can have severe effects on coastal ecosystems and human enterprises, such as fisheries and coastal power plant operations [1,4–6]. Therefore, mechanisms that cause the outbreaks of Aurelia spp. are of ecological and socio-economic interest.
Population dynamics have been studied for many species of marine organisms based on their life cycles (e.g., [7–9]). Similar to other scyphozoans, Aurelia spp. canonically have both asexual/sessile and sexual/planktonic stages in the life cycle (Fig 1, reviewed in [10]). In the metagenetic life cycle, a fertilized egg develops into a ciliated non-feeding planula larva, which settles down on a suitable substrate and then forms an asexual polyp. Environmental stimuli, such as temperature decline [11,12], cause a polyp to metamorphose into a strobila with several transverse segments in the body column, giving rise to multiple planktonic ephyrae, which grow up into sexually reproductive adult medusae. Aurelia spp. can significantly increase in number through asexual reproduction (e.g. budding and strobilation [12,13]) during the sessile period. Therefore, asexual reproduction has been recognized as a key process for outbreaks (e.g., [1,14]), and thus environmental factors that control asexual reproduction, such as temperature, prey abundance, and light, have been well studied [15–18].
Besides the above-mentioned canonical life cycle (hereafter described as “metagenetic life cycle”), it is reported from Japan and Germany that Aurelia spp. have another type of life cycle which lacks an asexual polyp stage (Fig 1, [12,19–22]). Contrary to holoplanktonic scyphozoans, such as Pelagia noctiluca and Periphylla periphylla, Aurelia spp. planulae of this alternative life cycle temporarily settle on substrate, and one settled planula produces one ephyra [12,19]. Though there is such a difference of Aurelia spp. from the holoplanktonic scyphozoans, this alternative life cycle of Aurelia spp. has been described as “direct development” because of the lack of an asexual polyp stage [10,19,23,24]. In the present study, we utilize the term “direct development” to be consistent with previous studies. In the direct-development type life cycle, Aurelia spp. do not increase in number through asexual reproduction during the polyp stage, and thus the mechanisms of outbreaks should be different between the two life cycle types. However, little is known about the direct development of Aurelia spp. Planulae from the direct-development type, which were collected in Tsuruga Bay, Japan, were larger (500–700 μm in the major axis, [22]) than planulae that transformed into polyps in other regions (220–375 μm, [25,26]). In addition, the medusae of Aurelia sp. sexually reproduce mainly from winter to spring in Tsuruga Bay [19,27], which is contrary to the summer spawning of Aurelia spp. reported in many other regions (reviewed in [23]). Therefore, larger planulae and spawning in the cold season are hypothesized to be key characteristics of the direct development of Aurelia spp. However, the two life cycle types of Aurelia spp. have never been studied to test this hypothesis or elucidate the population dynamics.
To test the hypothesis of seasonal alternation between the two life cycle types in Maizuru Bay, Japan, field observations and laboratory experiments were conducted throughout the year. In Maizuru Bay, Aurelia sp. medusae occur throughout the year and produce planulae even in winter (R Masuda, personal communication). Additionally, Aurelia sp. medusae representing both life cycle types were genetically analyzed to confirm that they belonged to a single species, because the genus Aurelia contains many cryptic species [28–31]. The two life cycle types were discussed in terms of controlling factors in the environment and ecological significance in their adaptations.
Materials and methods
Ethics statement
Field observations in the present study were approved by the harbormaster of Maizuru Bay. Aurelia sp. is not an endangered or protected species, and no permits were needed to sample Aurelia sp. medusae in Maizuru Bay. Samplings of medusae and zooplankton did not involve endangered or protected species.
Maizuru Bay
Maizuru Bay is a typical semi-enclosed water area with a surface area of approximately 23 km2, located in the western part of Wakasa Bay, Japan (Fig 2). The mean sea bottom depth is approximately 20 m, while the deepest area is 30 m deep. The tidal range is small (maximum: approximately 30 cm), because the Sea of Japan is connected to the Pacific Ocean through the narrow and shallow Korea, Tsugaru, and Soya Straits. Several rivers discharge into the bay. In this area, precipitation is rather high from summer to winter due to the rainy season, typhoon passages, and heavy snow on the Sea of Japan side.
Aurelia sp. collections
From November 2016 to January 2018, Aurelia sp. medusae were collected monthly from aggregations using a ship, the RV Ryokuyo Maru of Kyoto University, or from the pier of Maizuru Fisheries Research Station, Kyoto University. Aggregations of Aurelia sp. were found by visual surveys of the surface water or on the display of a fishfinder HDS-12 Gen2 Touch with a transducer of 200 kHz (Lowrance Electronics Inc., Tulsa, OK, USA). In each month, at least 21 medusae were collected with hand nets for the surface aggregations or with an obliquely towed conical net for the sub-surface aggregations (mouth diameter of 0.8 m, length of 2.3 m and mesh size of 5 mm). The bell diameters of collected medusae were measured and it was visually observed whether they bore planulae in brood sacs on the oral arms. In principle, samples of the planula, ovary, and tentacle at the bell margin were taken from five planula-bearing females in each month. If no medusa bore planulae, five females with bell diameters of > 14 cm, which is large enough for spawning in Maizuru Bay, were sampled. The planulae were rinsed thoroughly from the selected females with seawater. Collected planulae were rinsed with filtered seawater (0.22 μm) on a 60 μm mesh to remove mucus, and were then stored in plastic bags with the filtered seawater for transportation to the laboratory. The plastic bags were stored in thermos bottles, which were installed in a portable cooler to maintain a constant temperature during transportation. The mean water temperature change was 1.2 ± 0.8°C (mean value ± SD, the same applies hereinafter) during the transportation of approximately 16 h. In the laboratory, planulae were stored in an incubator at the same temperature as the in situ water temperature for several hours until further observations and developmental experiments. The collected ovaries were preserved in 5% buffered formalin solution for further observations of the oocytes. Tentacles were preserved in 99% ethanol for further DNA analysis. The bell diameter data were classified into two groups based on results of the developmental experiments: samples from December 2016 until May 2017 when the direct development occurred (direct-development period) and from June 2017 until October 2017 when all planulae developed into polyps (metagenetic period). Then, the bell diameters were compared between the periods of the two life cycle types with Welch’s t-test due to unequal sample numbers and unequal variances. In the present study, all statistical tests were carried out using R 3.4.3 [32].
Planula and oocyte observations
To estimate the total number of planulae per female medusa, which has been utilized as an indicator of egg production [33–35], the numbers of planulae in three replicates of 1-ml subsamples were counted for each sample. Twenty planulae for each female were measured for the length of the major and minor axes under a dissecting microscope. The volume of the planula was calculated from the length of the major and minor axes with an assumption in which the planula was in a spheroid shape. One of four ovaries was selected for each female. The selected ovaries were weighed, and then the weight was multiplied by four to estimate the total ovary weight of each female, as the weight of each ovary accounted for approximately 25%, ranging from 20 to 30%, of the total ovary weight (n = 5). The diameter of the oocytes (up to 200 in number), which are in a sphere-like shape, were measured under a dissecting microscope. The numbers of planulae per medusa and per ovary weight, the lengths of the planula, and the diameters of the top 5% of oocytes were compared between the direct-development and metagenetic periods using Welch’s t-test due to unequal sample numbers and unequal variances. The volume of the planulae and the top 5% of oocytes, i.e., most matured oocytes, were also compared in each month.
Environmental parameters
To relate the seasonality of the two life cycle types to environmental parameters, global horizontal irradiance, water temperature, salinity, chlorophyll fluorescence, and copepod abundance data were obtained. As weekly averaged global horizontal irradiance data have been unavailable at Maizuru Marine Observatory since March 2013, they were obtained from Fukui Local Meteorological Observatory, which is approximately 100 km away from Maizuru Bay (Fig 2). Weekly averaged global horizontal irradiance at the two observatories from January 2000 to March 2013 were highly correlated (r = 0.95), and thus the data from the Fukui Local Meteorological Observatory can be applicable to Maizuru Bay. The water temperature, salinity, and chlorophyll fluorescence were measured semi-weekly with the AAQ-RINKO (JFE Advantech Co., Ltd., Nishinomiya, Japan) at 0.1 m depth intervals from the observational pier of Maizuru Fisheries Research Station, Kyoto University (bottom depth: approximately 7 m, Fig 2). For water temperature and salinity, the mean surface and bottom layer values were calculated from data of 0.1–1.0 m below the surface and those of 1.0–2.0 m above the bottom, respectively. Chlorophyll fluorescence data were averaged throughout the water column.
Zooplankton were collected with vertical tows of a cylindrical-conical net (mouth diameter of 0.3 m, length of 1 m and mesh size of 100 μm) from near the bottom to the surface at eight stations monthly from November 2016 to January 2018, except for August and September 2017, when the ship was unavailable (Fig 2). The volume of water filtered through the net was calculated from the length of the wire and the area of the net mouth. Zooplankton samples were preserved in 5% buffered formalin solution for later identification under a dissecting microscope. In the laboratory, zooplankton samples were divided into subsamples using a wide-bore pipette immediately after thorough mixing by pipetting. The subsamples were observed until the number of zooplankton exceeded 200 individuals. Then, zooplankton were identified to the family or genus level and the sizes were measured individually by length. In the present study, copepod biomass was used as a proxy of the availability of prey zooplankton for Aurelia sp., because copepods are a major prey of Aurelia spp. (e.g., [36,37]). The carbon amount of each copepod taxon was calculated from the size according to Uye [38], and then data from the eight stations were averaged monthly.
Observation of ontogenetic development
To observe the development from planulae to the further stages, we incubated the collected planulae under darkness at the same temperature as the in situ water temperature. In the incubation, for each collected female, we first put approximately 100 planulae into a 24-well plate filled with 2 ml of filtered seawater (0.22 μm). Then, 20 settled planulae were randomly selected, and the others were removed to equalize the experimental settings among the wells. Ontogenetic development was observed daily under a dissecting microscope, and was categorized into three stages: settled individuals without tentacles, polyps, and liberated ephyrae. The weighted mean duration required for direct development from planula to ephyra (Dpla–ephy) was calculated for each experiment with the following equation:
DNA sequencing
The mitochondrial cytochrome c oxidase subunit I (COI) gene has been utilized to identify sibling species of Aurelia (e.g., [31,40]). To compare haplotypes of Aurelia sp. in Maizuru Bay between the two life cycle types and also to haplotypes of Aurelia coerulea, which has been treated as Aurelia sp. 1 and is distributed along the coast of Japan [31], a partial region of the COI gene was chosen for analysis. Total DNA was extracted from the ethanol-preserved tentacles of female medusa using the DNeasy Blood & Tissue Kit (Qiagen Inc., Hilden, Germany) according to the manufacturer’s manual. A partial region of the COI gene was amplified using the forward primer of HCO2198 (5′-taaacttcagggtgaccaaaaaatca-3′, [41]) and the reverse primer of AaCOIi-L (5′-gcccgtyytaataggrgggtttgg-3′, [40]) in PCR according to Matsumura et al. [42]. Direct sequencing was performed on the ABI 3500 Genetic Analyzer (Thermo Fisher Scientific, Waltham, MA, USA) using the HCO2198 primer and BigDye Terminator v3.1 Cycle Sequencing Kit (Thermo Fisher Scientific Inc., Waltham, MA, USA). The nucleotide sequence data were checked visually in Sequencher v5.4.6 (Gene Codes Corporation, Ann Arbor, MI, USA), and then misreads were corrected. The COI sequences from Maizuru Bay were classified into two groups based on the life cycle type (direct-development or metagenetic period). We obtained COI sequences of A. coerulea, which were collected in the North Pacific, from GenBank according to Scorrano et al. [31]: AY903080, AY903168, EU010386. Multiple COI sequences were aligned using Clustal W [43]. A COI phylogenetic tree was constructed by the maximum likelihood method and Kimura 2-parameter model using MEGA X [44,45]. Statistical support for tree nodes was determined by bootstrap analysis of 1000 replicates. A COI sequence of Aurelia limbata (AY903189), which is the most closely related species to A. coerulea [40], was used as an outgroup.
Results
Environmental parameters
The global horizontal irradiance was low (< 10 MJ m-2) from December until March and peaked in June (23 MJ m-2, Fig 3). The water temperature at the bottom layer showed the minimum (approximately 11°C) and maximum values (approximately 28°C) in February and August, respectively (Fig 3). The water temperature at the surface layer showed a sharp decline in winter due to snow. The salinity at the bottom layer was almost constant (32 to 33) throughout the observation, while the surface salinity often decreased in the rainy and snowy seasons. The surface salinity dropped in response to heavy precipitation incidental to typhoons (August 7, September 17, and October 22 in 2017; according to Japan Meteorological Agency). Chlorophyll fluorescence exhibited peaks repeatedly through the year without a clear seasonal trend (Fig 3). Copepod biomass was higher from January until June than in other months (Fig 3).
Development from planulae
A high proportion of planulae (approximately 90%) developed directly into ephyrae from December until May (direct-development period), while all planulae developed into polyps from July until October (metagenetic period, Figs 4 and 5). In June and November, the ontogenetic development of planulae was not observed, because no medusae bore planulae. Dpla–ephy declined exponentially with temperature increase, and the mean values at 10.0°C and 18.3°C were 16.7 days and 6.6 days, respectively (Fig 6). The mean total bell diameter (TBD) and central disc diameter (CDD) of the newly liberated ephyrae were 2.5 ± 0.4 mm and 0.9 ± 0.1 mm, respectively.
Seasonal variations in medusae, oocytes, and planulae
The bell diameters of medusae were almost constant throughout the study period, being approximately 15 cm on average (Fig 7). Proportions of small medusae with bell diameters of < 10 cm were higher in May and June (> 10%) than in the other months (< 5%). The bell diameters of medusae were not significantly different between the direct-development and metagenetic periods (Fig 7; P = 0.25, Welch’s t-test). The proportion of planula-bearing medusa was 13% to 64% in the direct-development and metagenetic periods, while no medusae bore planulae in June and November (Fig 7). The number of planulae per medusa was not significantly different between the direct-development and metagenetic periods (Fig 7; direct-development period: 736 ± 731 planulae medusae-1; metagenetic period: 1028 ± 810 planulae medusae-1; P = 0.50, Welch’s t-test). In contrast, the number of planulae per ovary weight was significantly higher in the metagenetic period than that in the direct-development period (Fig 6; direct-development period: 10.6 ± 10.0 planulae mg-1; metagenetic period: 44.2 ± 24.5 planulae mg-1; P < 0.01, Welch’s t-test).
During the metagenetic period, the monthly maximum of oocyte diameter ranged from 250 to 408 μm, while larger oocytes with a diameter of 504–571 μm were observed every month during the direct-development period (Fig 8). Although no planulae were found in November and June, oocytes in November and June were similar in size to those in December and July, respectively. The diameters of the top 5% of oocytes were significantly longer in the direct-development period (424 ± 37 μm) than in the metagenetic period (213 ± 17 μm, P < 0.01, Welch’s t-test; Fig 9). Similarly, the length of the major axis of the planula was significantly longer in the direct-development period (679 ± 142 μm) than in the metagenetic period (305 ± 75 μm, P < 0.01, Welch’s t-test). In both the planulae and the top 5% of oocytes, the volume was approximately 10-fold larger in the direct-development period than in the metagenetic period. Volumes of planulae and the top 5% of oocytes were not significantly different in each month (P > 0.05, Welch’s t-test), except for in May 2017 when the planulae were significantly larger than the top 5% of oocytes (Fig 9; P < 0.05, Welch’s t-test).
COI sequences
From the COI sequence data in 374 bp, 13 haplotypes were detected from 57 individuals (Fig 10), most of which matched those of Aurelia coerulea previously reported from other Japanese regions [42,46,47]. In the COI phylogenetic tree, specimens from the two life cycle types did not show clear differences. The sequences of Aurelia from Maizuru Bay were much more similar to those of A. coerulea than to those of A. limbata.
Discussion
Seasonality of the two life cycle types of A. coerulea in Maizuru Bay
In Maizuru Bay, Aurelia sp. exhibited seasonal alternation of ontogenetic development, which corresponded closely to seasonal changes in the size of oocytes and planulae (Figs 5, 8 and 9). In addition, from the COI haplotypes, both types of Aurelia sp. in Maizuru Bay were suggested to belong to A. coerulea (Fig 10, hereafter described as A. coerulea). Previous studies have demonstrated that larger egg sizes are generally associated with direct development in marine invertebrates (reviewed in [48]). However, intraspecific variations in egg size and developmental type have rarely been reported, especially from a single population [49,50].
The oocytes of the direct-development life cycle were reproduced by overwintering medusae (Figs 5 and 7). In Lake Nakaumi, a brackish lake in Japan, A. coerulea medusae are suggested to be able to overwinter when the water temperature in winter is higher than 5°C [51,52]. In the present study, the lowest water temperature in winter was higher than 5°C (Fig 3), and thus the water temperature was sufficient for the medusae to overwinter. In general, the longevity of Aurelia spp. medusa is from four months to less than two years [23]. Therefore, ephyrae, which are liberated through direct development from December until May, are suggested to develop into medusae and then reproduce in the following summer (metagenetic life cycle) and the next winter to spring (direct-development life cycle; Fig 11).
The canonical reproduction of the metagenetic life cycle was observed mainly in summer, similarly to the other regions (Fig 5, [23]). In Japan, ephyrae are liberated through strobilation mainly from December until March and then develop into young medusae in spring to early summer [53–55]. Assuming that the small medusae observed in May and June were new recruits, Maizuru Bay is similar in the timing of strobilation to the other regions in Japan. Therefore, the ephyrae liberated through strobilation are suggested to reproduce in the following summer (metagenetic life cycle) and the next winter to spring (direct-development life cycle) in the same way as in the direct-development life cycle (Fig 11).
Environmental factors in relation to the two life cycle types
In the present study, A. coerulea clearly switched the two life cycle types. The planula size changed discontinuously in June and November when no medusae bore planulae (Figs 7 and 9). Although oocytes were not histologically analyzed in the present study, the top 5% of oocytes were considered to be fully matured (i.e., eggs), because there was no significant difference in volume between the top 5% of oocytes and planulae in most months of the year (Fig 9). Intraspecific variations in egg size are well reported from a wide variety of taxa, and mother size and environmental stresses are suggested to be the major causes (reviewed in [58]). The mother sizes of A. coerulea in the present study were almost constant throughout the year. Thus, in Maizuru Bay, A. coerulea is suggested to switch life cycle type due not to the mother size, but to the environmental changes in June and November (hereafter described as switching periods).
The environmental stimuli that cause the switch in egg size should be different between the direct-development and metagenetic periods and should also show changes in the switching periods. Salinity, prey availability and water temperature are reported as the major environmental stresses affecting the egg size of marine animals [58]. In the present study, neither the salinity nor copepod density showed significant changes in the switching periods (Fig 3), and thus salinity and prey availability are unlikely to be the causes. Negative correlation between water temperature and egg size is a well-known phenomenon for a number of fish and invertebrate species (e.g., [48,59,60]). In the bottom layer, where A. coerulea usually aggregates in Japanese coastal waters [61–63], water temperatures of 19.8–22.0°C and 17.1–21.1°C were observed in June and November, respectively (Fig 3). As the water temperature was different between the direct-development and metagenetic periods and showed similar values (approximately 20°C) in the switching periods, a water temperature of 20°C is suggested to be the threshold for switching the egg size and ontogenetic development of A. coerulea in Maizuru Bay.
Egg size in relation to the two life cycle types
During the direct-development period, A. coerulea exhibited larger eggs than A. coerulea in metagenetic period (Fig 9). The carbon amount of the planula of Aurelia sp. with a major axis length of 200–320 μm was reported by [64] to be 0.36 μg, which is similar to the size of the metagenetic planulae in the present study (Fig 9). The carbon amount of the metagenetic planula is much smaller than that of the newly liberated ephyra in the present study (1.0 μg), estimated from the CDD with the allometric equation of Kamiyama [65]. This indicates that the metagenetic planulae do not store enough energy to transform directly to ephyrae. Assuming that the carbon amount of a planula is correlated to the volume, the planula during the direct-development period is estimated to contain approximately 4 μg C, which is even higher than the carbon amount of a newly liberated ephyra. Therefore, the high energy storage of the large eggs and planulae during the direct-development period is suggested to enable the planulae to develop into ephyrae without external energy input through feeding. From many studies on marine invertebrates, larger eggs are suggested to be necessary for direct development because of higher requirements for energy in early development (reviewed in [48]). Berrill [66] hypothesized that among scyphozoan jellyfish, larger eggs would be associated with direct development. Egg sizes of the holoplanktonic jellyfish P. noctiluca and P. periphylla, whose planulae directly develop into ephyrae or medusae, are larger than those of metagenetic scyphozoans (Table 1). This finding that egg sizes are closely related to seasonal alternation of the ontogenetic development of A. coerulea is in accordance with the hypothesis of Berrill. The number of planulae per medusa did not show significant difference between the direct-development and metagenetic periods, suggesting that the relationship between the size and number of eggs cannot be explained by the well-known trade-off [67–69]. However, the rate of egg production is more important than the number of existing eggs in population dynamics. Information on the turnover rate of eggs is necessary for a strict comparison of egg production rate between the direct-development and metagenetic periods.
Adaptive significances of the two life cycle types
Different life cycles generally have different adaptive significances. To discuss the adaptive significances of the two life cycle types of A. coerulea in Maizuru Bay, we focus on the temporal overlapping of ephyra liberation, the seasonal separation of spawning and difference in the duration of the sessile stages (sessile duration) between the two life cycle types.
Both life cycle types are suggested to produce ephyrae from winter to spring (Fig 11). In addition, the ephyrae produced from the direct-development type life cycle (2.5 mm in TBB) were similar in size to ephyrae produced from the strobilation of the metagenetic A. coerulea (2.1–3.7 mm, [54,84]). Thus, the performances of newly liberated ephyrae in survival are probably comparable between the two life cycle types. In winter to spring, the copepod biomass was high in Maizuru Bay (Fig 3), similarly to other regions [85]. The ephyrae of A. coerulea show higher growth under higher prey availability [86], and the rapid growth during the ephyra stage is suggested to decrease mortality from predation [84,87]. Therefore, A. coerulea may produce ephyrae under the abundant prey condition to decrease mortality in the ephyra stage with the rapid growth.
In the metagenetic life cycle, A. coerulea spawns in the warmer season, and then the polyps attach themselves to substrates for several months to a half year until winter strobilation (Fig 11). Warmer conditions are favorable for polyps to increase their number through asexual reproduction, while, under colder conditions, polyps tend to allocate ingested energy for somatic growth to prepare for the following strobilation [15]. Thus, the summer spawning of metagenetic A. coerulea would be adaptive in terms of increasing their population size. On the other hand, the metagenetic life cycle would face higher risks of mortality than the direct-development life cycle. During the polyp stage, the major causes of mortality are predation and spatial competition. Numerous marine animals are known to prey on Aurelia spp. polyps [88–92]. The polyps also face spatial competition against other sessile animals, such as the mussel Mytilus galloprovincialis, which can eliminate polyps of A. coerulea [57]. The risks of predation and spatial competition during the polyp stage increase as the sessile duration increases, while the sessile duration of A. coerulea in the metagenetic life cycle (several months to a half year) is much longer than that of A. coerulea in the direct-development life cycle (less than half a month; Fig 6). Therefore, the summer spawning of metagenetic A. coerulea is suggested to have an advantage in increasing their population size through asexual reproduction and have risks of predation and spatial competition (Fig 11).
In comparison to the metagenetic life cycle, the direct-development life cycle has a disadvantage in increasing their population size in the short term because of a lack of asexual reproduction during the polyp stage. The winter spawning of the direct-development life cycle may benefit from a lower mortality. As planulae of A. coerulea are lecithotrophic, they die naturally unless they settle down to suitable substrates before they run out of their energy storage. In other words, a longer longevity of planula leads a higher potential to encounter suitable substrates, although there are some risks of predation and dispersion during the planktonic period. The sizes of planulae of three pocilloporid coral species show intraspecific positive correlation with longevity [93]. In addition, the planktonic larval duration of marine animals generally increases as the temperature decreases [94]. Thus, planulae in the direct-development life cycle, which are large and produced in winter, can have a longer longevity and thus a higher survival rate than those in the metagenetic life cycle. In addition, the shorter sessile duration of A. coerulea in the direct-development life cycle until ephyra liberation plays an important role in producing ephyrae in an appropriate season. Under the cold condition of 15°C, A. coerulea polyps do not transform into strobilae within five months from the settlement of planulae [95]. Thus, in the case of the winter spawning, a direct-development life cycle with a shorter sessile duration would be necessary for A. coerulea to produce ephyrae in winter to spring with high prey availability. Therefore, the direct-development life cycle plays an important role in producing ephyrae in the appropriate season, and potentially decreases the mortality of planulae both from predation and natural death. However, the direct-development life cycle also has a disadvantage in increasing the population size in the short term (Fig 11).
Conclusions
While asexual reproduction in the metagenetic life cycle of A. coerulea has been recognized as a key reason for their outbreaks, our results prompt a rethink of its reproductive strategies. Our study clarified that in Maizuru Bay, A. coerulea is capable to alternate two strategies of reproduction in response to environmental conditions. These alternative life cycle strategies may play important roles in adapting to various environmental conditions, as have been suggested for other organisms (reviewed in [49]). From the present study, the overwintering of mature medusae is suggested to be an indicator of the direct-development life cycle. The direct development of Aurelia spp. has been reported from only a few regions [12,19–21]; this may be because most studies concerning the medusa stage have been conducted in warmer seasons. To clarify the importance of the direct-development life cycle in the population dynamics of Aurelia spp., further investigations of the mechanisms of overwintering and the generality of the direct development among overwintering populations are required. As Helm et al. [24] suggested, the alternation of two types of life cycles of A. coerulea in Maizuru Bay could be a suitable model for studying the evolution of direct-development life cycles in scyphozoan jellyfish.
Supporting information
S1 File [xlsx]
Raw data for each figure.
Zdroje
1. Purcell JE, Uye S, Lo W-T. Anthropogenic causes of jellyfish blooms and their direct consequences for humans: A review. Mar Ecol Prog Ser. 2007;350: 153–174.
2. Richardson AJ, Bakun A, Hays GC, Gibbons MJ. The jellyfish joyride: causes, consequences and management responses to a more gelatinous future. Trends Ecol Evol. 2009;24: 312–322. doi: 10.1016/j.tree.2009.01.010 19324452
3. Robinson KL, Graham WM. Long-term change in the abundances of northern Gulf of Mexico scyphomedusae Chrysaora sp. and Aurelia spp. with links to climate variability. Limnol Oceanogr. 2013;58: 235–253.
4. Möller H. Reduction of a larval herring population by jellyfish predator. Science. 1984;224: 621–622. doi: 10.1126/science.224.4649.621 17838355
5. Uye S, Shimauchi H. Population biomass, feeding, respiration and growth rates, and carbon budget of the scyphomedusa Aurelia aurita in the Inland Sea of Japan. J Plankton Res. 2005;27: 237–248.
6. Dong Z, Liu D, Keesing JK. Jellyfish blooms in China: Dominant species, causes and consequences. Mar Pollut Bull. 2010;60: 954–963. doi: 10.1016/j.marpolbul.2010.04.022 20553695
7. Roughgarden J, Gaines S, Possingham H. Recruitment dynamics in complex life cycles. Science. 1988;241: 1460–1466. doi: 10.1126/science.11538249 11538249
8. Eckman JE. Closing the larval loop: Linking larval ecology to the population dynamics of marine benthic invertebrates. J Exp Mar Biol Ecol. 1996;200: 207–237.
9. Carrier TJ, Reitzel AM, Heyland A. Evolutionary Ecology of Marine Invertebrate Larvae. Oxford: Oxford University Press; 2018.
10. Arai MN. A functional biology of Scyphozoa. London: Chapman & Hall; 1997.
11. Kakinuma Y. On some factors for the differentiations of Cladonema uchidai and Aurelia aurita. Bull Mar Biol Stn Asamushi. 1962;11: 81–85.
12. Kakinuma Y. An experimental study of the life cycle and organ differentiation of Aurelia aurita Lamarck. Bull Mar Biol Stn Asamushi. 1975;15: 101–112.
13. Schiariti A, Morandini AC, Jarms G, von Glehn Paes R, Franke S, Mianzan H. Asexual reproduction strategies and blooming potential in Scyphozoa. Mar Ecol Prog Ser. 2014;510: 241–253.
14. Lucas CH, Graham WM, Widmer C. Jellyfish Life Histories: Role of Polyps in Forming and Maintaining Scyphomedusa Populations. Advances in Marine Biology. Amsterdam: Academic Press; 2012. p. 133–196. doi: 10.1016/B978-0-12-394282-1.00003-X 22877612
15. Han C-H, Uye S. Combined effects of food supply and temperature on asexual reproduction and somatic growth of polyps of the common jellyfish Aurelia aurita s.l. Plankton Benthos Res. 2010;5: 98–105.
16. Willcox S, Moltschaniwskyj NA, Crawford C. Asexual reproduction in scyphistomae of Aurelia sp.: Effects of temperature and salinity in an experimental study. J Exp Mar Biol Ecol. 2007;353: 107–114.
17. Liu W-C, Lo W-T, Purcell JE, Chang H-H. Effects of temperature and light intensity on asexual reproduction of the scyphozoan, Aurelia aurita (L.) in Taiwan. Hydrobiologia. 2009;616: 247–258.
18. Thein H, Ikeda H, Uye S. The potential role of podocysts in perpetuation of the common jellyfish Aurelia aurita s.l. (Cnidaria: Scyphozoa) in anthropogenically perturbed coastal waters. Hydrobiologia. 2012;690: 157–167.
19. Yasuda T. Ecological studies on the jelly-fish, Aurelia aurita (Linné), in Urazoko Bay, Fukui Prefecture-XI. An observation on ephyra formation. Publ Seto Mar Biol Lab. 1975;22: 75–80.
20. Haeckel EHPA. Metagenesis und hypogenesis von Aurelia aurita. Jena: G. Fischer; 1881.
21. Hirai E. On the developmental cycles of Aurelia aurita and Dactylometra pacifica. Bull Mar Biol Stn Asamushi. 1958;9: 81.
22. Yasuda T. Jellyfish: UFO in the sea. Tokyo: Koseisha-koseikaku; 2003. (in Japanese)
23. Lucas CH. Reproduction and life history strategies of the common jellyfish, Aurelia aurita, in relation to its ambient environment. Hydrobiologia. 2001;451: 229–246.
24. Helm RR, Tiozzo S, Lilley MKS, Lombard F, Dunn CW. Comparative muscle development of scyphozoan jellyfish with simple and complex life cycles. EvoDevo. 2015;6: 11. doi: 10.1186/s13227-015-0005-7 25932322
25. Gröndahl F. Interactions between polyps of Aurelia aurita and planktonic larvae of scyphozoans: an experimental study. Mar Ecol Prog Ser. 1988;45: 87–93.
26. Conley K, Uye S. Effects of hyposalinity on survival and settlement of moon jelly fish (Aurelia aurita) planulae. J Exp Mar Biol Ecol. 2015;462: 14–19.
27. Yasuda T. Ecological studies on the jelly-fish, Aurelia aurita in Urazoko Bay, Fukui Prefecture–IV. Monthly change in the bell-length composition and breeding season. Bull. Japanese Soc. Sci. Fisheries. 1971;37: 364–370. (in Japanese with English abstract)
28. Dawson MN, Martin LE. Geographic variation and ecological adaptation in Aurelia (Scyphozoa, Semaeostomeae): Some implications from molecular phylogenetics. Hydrobiologia. 2001;451: 259–273.
29. Ki J-S, Hwang D-S, Shin K, Yoon WD, Lim D, Kang YS, et al. Recent moon jelly (Aurelia sp.1) blooms in Korean coastal waters suggest global expansion: examples inferred from mitochondrial COI and nuclear ITS-5.8S rDNA sequences. ICES J Mar Sci. 2008;65: 443–452.
30. Chiaverano LM, Bayha KW, Graham WM. Local versus generalized phenotypes in two sympatric Aurelia species: Understanding jellyfish ecology using genetics and morphometrics. PLoS One. 2016;11: e0156588. doi: 10.1371/journal.pone.0156588 27332545
31. Scorrano S, Aglieri G, Boero F, Dawson MN, Piraino S. Unmasking Aurelia species in the Mediterranean Sea: An integrative morphometric and molecular approach. Zool J Linn Soc. 2017;180: 243–267.
32. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2016. Available from: https://www.R-project.org/.
33. Schneider G. Larvae production of the common jellyfish Aurelia aurita in the western Baltic 1982–1984. Kiel Meeresforsch. 1988;6: 295–300.
34. Lucas CH. Population dynamics of Aurelia aurita (Scyphozoa) from an isolated brackish lake, with particular reference to sexual reproduction. J Plankton Res. 1996;18: 987–1007.
35. Lucas CH, Lawes S. Sexual reproduction of the scyphomedusa Aurelia aurita in relation to temperature and variable food supply. Mar Biol. 1998;131: 629–638.
36. D’Ambra I, Graham WM, Carmichael RH, Hernandez FJ. Dietary overlap between jellyfish and forage fish in the northern Gulf of Mexico. Mar Ecol Prog Ser. 2018;587: 31–40.
37. Ishii H, Tanaka F. Food and feeding of Aurelia aurita in Tokyo Bay with an analysis of stomach contents and a measurement of digestion times. Hydrobiologia. 2001;451: 311–320.
38. Uye S. Length-weight relationships of important zooplankton from the Inland Sea of Japan. J Oceanogr Soc Jpn. 1982;38: 149–158.
39. Straehler-Pohl I, Jarms G. Identification key for young ephyrae: A first step for early detection of jellyfish blooms. Hydrobiologia. 2010;645: 3–21.
40. Dawson MN, Jacobs DK. Molecular evidence for cryptic species of Aurelia aurita (Cnidaria, Scyphozoa). Biol Bull. 2001;200: 92–96. doi: 10.2307/1543089 11249217
41. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3: 294–299. 7881515
42. Matsumura K, Kamiya K, Yamashita K, Hayashi F, Watanabe I, Murao Y, et al. Genetic polymorphism of the adult medusae invading an electric power station and wild polyps of Aurelia aurita in Wakasa Bay, Japan. J Mar Biol Assoc U K. 2005;85: 563–568.
43. Larkin MA, Blackshields G, Brown NP, Chenna R, Mcgettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23: 2947–2948. doi: 10.1093/bioinformatics/btm404 17846036
44. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35: 1547–1549. doi: 10.1093/molbev/msy096 29722887
45. Kimura M. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16: 111–120. doi: 10.1007/bf01731581 7463489
46. Endo N, Nogata Y, Sato K, Yoshimura E, Matsumura K, Sakaguchi I, et al. Genetic relationship of moon jellyfish Aurelia aurita s.l. in Tokyo Bay and Ise Bay based on nucleotide sequences of mitochondrial DNA. Bull Plankton Soc Japan. 2011;58: 1–11. (in Japanese with English abstract)
47. Dawson MN, Sen Gupta A, England MH. Coupled biophysical global ocean model and molecular genetic analyses identify multiple introductions of cryptogenic species. Proc Natl Acad Sci U S A. 2005;102: 11968–11973. doi: 10.1073/pnas.0503811102 16103373
48. Marshall DJ, Krug PJ, Kupriyanova EK, Byrne M, Emlet RB. The biogeography of marine invertebrate life histories. Annu Rev Ecol Evol Syst. 2012;43: 97–114.
49. Chia F-S, Gibson G, Qian P-Y. Poecilogony as a reproductive strategy of marine invertebrates. Oceanol Acta. 1996;19: 203–208.
50. Vendetti JE, Trowbridge CD, Krug PJ. Poecilogony and population genetic structure in Elysia pusilla (Heterobranchia: Sacoglossa), and reproductive data for five sacoglossans that express dimorphisms in larval development. Integr Comp Biol. 2012;52: 138–150. doi: 10.1093/icb/ics077 22659202
51. Han C-H, Kawahara M, Uye S. Seasonal variations in the trophic relationship between the scyphomedusa Aurelia aurita s.l. and mesozooplankton in a eutrophic brackish-water lake, Japan. Plankton Benthos Res. 2009;4: 14–22.
52. Makabe R, Kurihara T, Uye S. Spatio-temporal distribution and seasonal population dynamics of the jellyfish Aurelia aurita s.l. studied with Dual-frequency IDentification SONar (DIDSON). J Plankton Res. 2012;34: 936–950.
53. Miyake H, Iwao K, Kakinuma Y. Life history and environment of Aurelia aurita. South Pac Study. 1997;17: 273–285.
54. Makabe R, Furukawa R, Takao M, Uye S. Marine artificial structures as amplifiers of Aurelia aurita s.l. blooms: A case study of a newly installed floating pier. J Oceanogr. 2014;70: 447–455.
55. Watanabe T, Ishii H. In situ estimation of ephyrae liberated from polyps of Aurelia aurita using settling plates in Tokyo Bay, Japan. Hydrobiologia. 2001;451: 247–258.
56. Kim YS, Moon TS. Filtering rate with effect of water temperature and size of two farming ascidians Styela clava and S. plicata, and a farming mussel Mytilus edulis. J Korean Fish Soc. 1998;31: 272–277. (in Korean with English abstract)
57. Ishii H, Katsukoshi K. Seasonal and vertical distribution of Aurelia aurita polyps on a pylon in the innermost part of Tokyo Bay. J Oceanogr. 2010;66: 329–336.
58. Marshall DJ, Keough MJ. The evolutionary ecology of offspring size in marine invertebrates. Adv Mar Biol. 2007;53: 1–60. doi: 10.1016/S0065-2881(07)53001-4 17936135
59. Thatje S, Hall S. The effect of temperature on the evolution of per offspring investment in a globally distributed family of marine invertebrates (Crustacea: Decapoda: Lithodidae). Mar Biol. 2016;163: 1–9.
60. Bagenal TB. The interrelation of the size of fish eggs, the date of spawning and the production cycle. J Fish Biol. 1971;3: 207–219.
61. Suzuki KS, Yasuda A, Murata Y, Kumakura E, Yamada S, Endo N, et al. Quantitative effects of pycnocline and dissolved oxygen on vertical distribution of moon jellyfish Aurelia aurita s.l.: A case study of Mikawa Bay, Japan. Hydrobiologia. 2016;766: 151–163.
62. Suzuki KS, Kumakura E, Endo N, Ishii H, Nogata Y. Effects of water mass structure on vertical distribution of moon jellyfish Aurelia aurita s.l. within aggregations in four Japanese coastal areas. Bull Plankton Soc Japan. 2017;64: 114–123. (in Japanese with English abstract)
63. Itoh H, Tachibana A, Nomura H, Tanaka Y, Furota T, Ishimaru T. Vertical distribution of planktonic copepods in Tokyo Bay in summer. Plankton Benthos Res. 2011;6: 129–134.
64. Schneider G, Weisse T. Metabolism measurements of Aurelia aurita planulae larvae, and calculation of maximal survival period of the free swimming stage. Helgol Meeresunters. 1985;39: 43–47.
65. Kamiyama T. Planktonic ciliates as food for the scyphozoan Aurelia coerulea: feeding and growth responses of ephyra and metephyra stages. J Oceanogr. 2017; 1–11.
66. Berrill NJ. Developmental analysis of scyphomedusae. Biol Rev. 1949;24: 393–409. doi: 10.1111/j.1469-185x.1949.tb00581.x 24536313
67. Parker GA, Begon M. Optimal egg size and clutch size: Effects of environment and maternal phenotype. Am Nat. 1986;128: 573–592.
68. Smith CC, Fretwell SD. The optimal balance between size and number of offspring. Am Nat. 1974;108: 499–506.
69. Thorson G. Reproductive and larval ecology of marine bottom invertebrates. Biol Rev. 1950;25: 1–45. doi: 10.1111/j.1469-185x.1950.tb00585.x 24537188
70. Komai T. On Stephanoscyphus and Nausithoe. Mem Coll Sci Kyoto Imp Univ Ser B. 1935;10: 290–339.
71. Morandini AC, Da Silveira FL. Sexual reproduction of Nausithoe aurea (Scyphozoa, Coronatae). Gametogenesis, egg release, embryonic development, and gastrulation. Sci Mar. 2001;65: 139–149.
72. Metschnikoff E. Embryologische Studien an Medusen, Ein Beitrag zur Genealogie der Primitiv-Organe. Hölder A, editor. Vienna; 1886.
73. Conklin EG. The habits and early development of Linerges mercurius. Publs Carnegie Instn. 1908;103: 155–170.
74. Jarms G, Båmstedt U, Tiemann H, Martinussen MB, Fosså JH. The holopelagic life cycle of the deep-sea medusa Periphylla periphylla (Scyphozoa, Coronatae). Sarsia. 1999;84: 55–65.
75. Uchida T. The anatomy and development of rhizostome medusa, Mastigias papua L. Agassiz with observation on the phylogeny of Rhizostomae. J Fac Sci Imp Univ Tokyo, Sect IV, Zool. 1926;1: 45–95.
76. Goette A. Entwicklungsgeschichte der Aurelia aurita und Cotylorhiza tuberculata. Hamburg und Leipzig: Leopold Voss; 1887.
77. Ikeda H, Ohtsu K, Uye S. Structural changes of gonads during artificially induced gametogenesis and spawning in the giant jellyfish Nemopilema nomurai (Scyphozoa: Rhizostomeae). J Mar Biol Assoc U K. 2011;91: 215–227.
78. Hyde IH. Entwickelungsgeschichte einiger Scyphomedusen. Zeitschrift für Wissenschaftliche Zool. 1894;58: 531–564.
79. Goette A. Vergleichende Entwicklungsgeschichte von Pelagia noctiluca Pér. Z wiss Zool. 1893;55: 645–695.
80. Hadzi J. Über die Podozysten der Scyphopolypen. Biol Zentralb. 1912;32: 52–60.
81. Teissier G. La croissance embryonnaire de Chrysaora hysoscella (L.). Arch Zool Exp Gen. 1929;69: 137–178.
82. Littleford RA. The life cycle of Dactylometra quinquecirrha, L. Agassiz in the Chesapeake Bay. Biol Bull. 1939;77: 368–381.
83. Hargitt CW, Hargitt GT. Studies on the development of scyphomedusae. J Morphol. 1910;21: 217–262.
84. Suzuki KS, Kumakura E, Nogata Y. Incidental consumption of ephyrae of moon jellyfish Aurelia aurita s.l. by three filter-feeding sessile organisms: laboratory experiments. Fish Sci. 2016;82: 923–930.
85. Iguchi N. Spatial/temporal variations in zooplankton biomass and ecological characteristics of major species in the southern part of the Japan Sea: A review. Prog Oceanogr. 2004;61: 213–225.
86. Ishii H, Kojima S, Tanaka Y. Survivorship and production of Aurelia aurita ephyrae in the innermost part of Tokyo Bay, Japan. Plankton Biol Ecol. 2004;51: 26–35.
87. Lamb PD, Hunter E, Pinnegar JK, van der Kooij J, Creer S, Taylor MI. Cryptic diets of forage fish: jellyfish consumption observed in the Celtic Sea and western English Channel. J Fish Biol. 2019;94: 1026–1032. doi: 10.1111/jfb.13926 30746684
88. Miyajima-Taga Y, Masuda R, Morimitsu R, Ishii H, Nakajima K, Yamashita Y. Ontogenetic changes in the predator–prey interactions between threadsail filefish and moon jellyfish. Hydrobiologia. 2016;772: 175–187.
89. Feng S, Wang S, Zhang G, Sun S, Zhang F. Selective suppression of in situ proliferation of scyphozoan polyps by biofouling. Mar Pollut Bull. 2017;114: 1046–1056. doi: 10.1016/j.marpolbul.2016.10.062 27890264
90. Hoover RA, Armour R, Dow I, Purcell JE. Nudibranch predation and dietary preference for the polyps of Aurelia labiata (Cnidaria: Scyphozoa). Hydrobiologia. 2012;690: 199–213.
91. Takao M, Okawachi H, Uye S. Natural predators of polyps of Aurelia aurita s.l. (Cnidaria: Scyphozoa: Semaeostomeae) and their predation rates. Plankton Benthos Res. 2014;9: 105–113.
92. Hernroth L, Gröndahl F. On the biology of Aurelia aurita (L.) 3. Predation by Coryphella verrucosa (Gastropoda, Opisthobranchia), a major factor regulationg the development of Aurelia populations in the Gullmar Fjord, western Sweden. Ophelia. 1985;24: 37–45.
93. Isomura N, Nishihira M. Size variation of planulae and its effect on the lifetime of planulae in three pocilloporid corals. Coral Reefs. 2001;20: 309–315.
94. O’Connor MI, Bruno JF, Gaines SD, Halpern BS, Lester SE, Kinlan BP, et al. Temperature control of larval dispersal and the implications for marine ecology, evolution, and conservation. Proc Natl Acad Sci U S A. 2007;104: 1266–1271. doi: 10.1073/pnas.0603422104 17213327
95. Ishii H, Watanabe T. Experimental study of growth and asexual reproduction in Aurelia aurita polyps. Sessile Org. 2003;20: 69–73.
Článek vyšel v časopise
PLOS One
2019 Číslo 11
- Řekněte mi, co jíte, a já vám řeknu, jaká je úroveň vaší kognice
- Duševní onemocnění se u teenagerů přenáší v rámci školní třídy
- Od kostiček k lidské kůži. Jak si postavit 3D biotiskárnu?
- Pomohou nám s detekcí karcinomu plic i včely?
- Původní český test ALBA pro odhalení mírných kognitivních deficitů získal certifikaci MZ ČR
Nejčtenější v tomto čísle
- A daily diary study on maladaptive daydreaming, mind wandering, and sleep disturbances: Examining within-person and between-persons relations
- A 3’ UTR SNP rs885863, a cis-eQTL for the circadian gene VIPR2 and lincRNA 689, is associated with opioid addiction
- A substitution mutation in a conserved domain of mammalian acetate-dependent acetyl CoA synthetase 2 results in destabilized protein and impaired HIF-2 signaling
- Molecular validation of clinical Pantoea isolates identified by MALDI-TOF