Abstract
Introduction:
Understanding how organisms adapt to complex and variable environments requires in situ analysis of ecological performance in different environments.
Methods:
In this multi-year study, we investigate the photophysiological performance of dinoflagellate endosymbionts (Symbiodiniaceae) of the golden jellyfish Mastigias papua across a natural environmental gradient of marine lakes and semienclosed coves in the Palau Archipelago (Western Caroline Islands, Pacific Ocean).
Results:
These sites vary in terms of age, distance from the nearby ocean (lagoon) and underwater light conditions, all of which influence the ecology of the host and that of its symbionts. Jellyfish endosymbionts from isolated lakes ‐ remote from the ocean and characterized by more turbid conditions – were found to exhibit significantly lower photosynthetic performance (quantum yield and PSII performance index) than endosymbionts from sites more connected to the nearby ocean and with clearer waters. The pigment composition of endosymbionts varied little between those collected from inland lakes and ocean-connected coves, and showed little sign of adaptation to light spectra dominated by reddish wavelengths.
Discussion:
This suggests a mismatch between the optical characteristics of waters in more turbid lake environments and those of pigments that appear to remain optimized for blue light use in marine dinoflagellates. While published studies show evidence of morphological and behavioral adaptations in Mastigias jellyfish living in Palau’s marine lakes, this study suggests that their symbionts are limited by reduced pigment plasticity and lower photosynthetic performance. These results highlight the importance of symbiont characteristics and physiological compatibility in driving host adaptation, and the potential for evolutionary mismatch in tightly integrated symbioses, under novel environmental conditions.
1 Introduction
In recent years, important progresses have been made in understanding the molecular biology of cells, mapping genomes and identifying gene functions (
;
;
). However, a substantial gap remains in translating molecular-level knowledges into an understanding of large-scale ecological and evolutionary processes, such as for example how organisms living in ecosystems adapt to complex and changing environments (
;
;
). While molecular biology tools allow to make important functional predictions based on gene content and functions, these predictions are often confounded by several factors. Notably, the flexibility of organismal behavior and physiology introduces significant variability, and many genes still remain poorly characterized or have an unknown function (
;
). This may partly limit the predictive power of purely genomic approaches, when it comes to studying adaptation directly at ecosystem level. In this context, empirical field studies remain invaluable. It is by studying phenotypes
in situ, in a wide range of environmental contexts, that hypotheses concerning organism’s adaptation can be tested and refined (
). Such fieldwork remains crucial for a better understanding of how molecular and genetic mechanisms scale up to influence organism’s fitness, population dynamics, or interactions under contrasted natural environmental conditions. Importantly, organisms do not adapt in isolation. All life forms interact with other co-occurring species - through predation, competition, mutualism, etc., - and these biotic interactions are critical components of their adaptive landscapes (
). For example, many species depend on obligate symbioses with species-specific microbial symbionts for their energy and nutrients requirements (
). A typical example is animals that harbor photosynthetic endosymbionts, such as cyanobacteria or algae, which provide their hosts with photosynthetically fixed carbon and other metabolic benefits (
;
). These symbioses add a layer of complexity to adaptation, because the evolutionary rates of hosts and their symbionts may differ significantly. Species with short generation times, such as unicellular algae or bacteria, can evolve and adapt rapidly. In contrast, their multicellular hosts, which often have much longer life cycles, may adapt over far longer timescales. This mismatch raises questions about how such symbiotic association respond to environmental changes: Do for example the faster-evolving symbionts buffer the host against environmental stress? Or does this disparity introduce instability into the interaction? As
pointed out, adaptation in interacting organisms may occur on different timescales, with implications for their co-evolutionary dynamics. The phylum Cnidaria offers several examples of such tightly integrated symbiotic associations. While reef-building corals are perhaps the best-known cnidarian hosts of photosynthetic symbionts, many jellyfish species also depend on endosymbiotic algae, commonly dinoflagellates of the family Symbiodiniaceae, often referred to as zooxanthellae (
). These symbioses are central to the biology and ecology of the hosts. One of the examples is found in the golden jellyfish
Mastigias papua, which form stable holobionts with their dinoflagellate endosymbionts (
). These jellyfish are widely distributed across tropical Pacific marine environments but are especially well-known for forming large and isolated populations in Palau’s marine lakes (
;
). In the archipelago of Palau, located in the western Pacific Ocean, ∼70 marine lakes–small and saline water masses within the karstic limestone topography of the islands–have been described (
). These lakes are partially isolated from the nearby ocean, with more or less important connections via tunnels and/or narrow subterranean fissures. As a result, the exchange of water - and thus of marine organisms and nutrients - is partly limited. Over thousands of years, these geographic barriers have allowed each lake to evolve as unique ecosystems, with distinct physicochemical properties and biological communities (
). The most famous of these lakes is Jellyfish Lake (Ongeim’l Tketau), which is inhabited by a population of the golden jellyfish
Mastigias papuathat undertake daily horizontal migrations in response to sunlight (
). However, this lake is just one among many in Palau Islands hosting populations of different subspecies of
M. papua(
;
). These lakes are characterized by a gradient of environmental conditions, ranging from semi-enclosed coves, closely connected to the neighboring ocean lagoon, to fully land-locked lakes. This environmental gradient provides a set of natural laboratories for studying how isolation, water biogeochemistry, and light availability influence the ecological and evolutionary processes of marine organisms. These marine lakes were formed approximately 12,000 years ago, at the end of the last ice age, when rising sea levels filled several low-lying basins in the limestone islands (
). The timing and extent of lake formation varied depending on their location and depth. Deeper basins filled earlier, while shallower basins filled later on (
). Each population of
M. papuashows evidences of variation in morphology, behavior, and genetic markers, despite being derived from a common ancestor (
;
). Remarkably, all populations continue to rely on the same genus of photosynthetic endosymbionts (Dinophyceae, Suessiales, Symbiodiniaceae). It is likely even the same clade, based on current phylogenetic evidence (
), which was recently named
Cladocopium(
). However, the environmental conditions to which these holobionts are subjected differ from site to site, particularly in terms of depth, nutrient content, light availability, and trophic status (oligotrophic, mesotrophic, eutrophic). In this respect, light is an essential resource for endosymbionts. Coastal coves and lakes with strong connections to the nearby ocean are generally oligotrophic, with clearer waters and light regimes similar to those found in open marine environments. By contrast, lakes further inland from the lagoon often receive humic substances and organic matter from the surrounding tropical forests, via surface runoff. These inputs of dissolved organic matter, combined with the development of phytoplankton, result in the browning and opacification of surface waters, and a concomitant reduction in light penetration and spectral quality (
). These significant differences in luminosity between lakes and coves are then likely to influence the physiology of endosymbionts, possibly leading to plastic and evolutionary responses. To explore these dynamics, we used the natural variations between Palau’s marine lakes mentioned above to study the interaction between environmental conditions and the photosynthetic characteristics and performance of symbionts. Specifically, we addressed the following questions:
Does the pigment composition of endosymbiotic algae vary in response to different light environments between open and closed marine systems?
Do endosymbiotic algae exhibit different photosynthetic performance between open marine systems and closed marine lakes along a shore distance gradient?
Are there specific photo physiological characteristics that suggest adaptation in endosymbiotic algae living in different lacustrine environments?
By examining these questions, using samples of jellyfish collected over several years in these different marine lake environments, we aim to shed light on the mechanisms underlying host-symbiont interactions and their ability to adapt - or acclimatize - to environmental change.
2 Materials and methods
2.1 Study area and sample collection
Field surveys were conducted in the Palau Archipelago, Micronesia (7°–8° N, 133°–135° E), during multiple expeditions between 2010 and 2015 (Figure 1; Table 1). Sampling sites included coastal marine lakes and coves inhabited by symbiotic jellyfish of the genera Mastigias (Cnidaria: Scyphozoa: Rhizostomeae; ; ).
FIGURE 1
TABLE 1
| Site name(a) | Habitat type | Jellyfish sampling | Geography | BGC | Optical properties | |||||
|---|---|---|---|---|---|---|---|---|---|---|
| Years (i)(b) | ni(j) | ntot(j) | SD (m) | MD(c) (m) | Age (d) (years) | TP(e) (µg L-1) | KPAR(f) (m-1) | Secchi depth(g) (m) | ||
| Clear (CLM) | Lake | 2013, 2015 | 3, 5 | 8 | 244 | 34 | ? | 37.74 | 0.34 | 8.7 (1.0) |
| Goby (GLK) | Lake | 2010, 2011, 2015 | 9(5), 10(5), 5 | 24(10) | 89 | 15 | 9,000 | 68.21 | 0.42 | 5.5 (0.3) |
| Ngermeuangel (BJLK) | Lake | 2010, 2011, 2013 | 16(5), 5(5), 2 | 23(10) | 184 | 38 | 15,000 | 28.66 | 0.10 | 12.7 (0.4) |
| Ongael (OLO) | Lake | 2010, 2013 | 5(5), 3 | 8(5) | 119 | 4 | 6,000 | 39.99 | 0.66 | 3.2-Bottom |
| Ongeim'1 tketau (OTM) | Lake | 2010, 2011, 2015 | 8(6), 5(5), 5 | 18(11) | 96 | 33 | 10,000 | 29.32 | 0.39 | 8.4 (0.2) |
| Tarzan | Lake | 2015 | 5 | 5 | 36 | ∼5–10 | ? | 3.40 | (0.03)(e) | n.d |
| Jurassic (ulong island) | Cove/Tidal sill | 2010 | 10(5) | 10(5) | 23 | ∼5–10 | ? | 9.98 | (0.09)(e) | Bottom |
| Secret | Cove/Tidal sill | 2015 | 3 | 3 | 18 | >10 | ? | 14.72 | (0.14)(e) | n.d |
| Ulebsechel (RCA, Risong bay) | Cove/Tidal sill | 2015 | 3 | 3 | 15 | 20 | ? | 24.09 | 0.17 | 10.8 (0.3) |
| Unnamed (ngerchaol island) | Cove | 2011 | 3(3) | 3(3) | 0 | ∼5–10 | ? | 22.34 | (0.21)(e) | n.d |
| Unnamed (Risong bay, site2) | Cove | 2015 | 3 | 3 | 0 | ∼5–10 | ? | n.d | n.d | n.d |
| | | | Total = | 108(44) | | | | | | |
Summary of site characteristics and jellyfish sampling information from 2010 to 2015.
Sampling dates were: 2010 (19 August–09 September), 2011 (27 July–09 August), 2013 (06–27 August), and 2015 (18 April–09 May). Columns indicate site name, habitat type, and sampling year (i); with total number of individuals collected per year, ni(j), and total number of individuals sampled per site, ntot(j), with (j) the number of individuals subsampled for photosynthetic pigment analysis. Reported parameters include geographical and biogeochemical (BGC) characteristics—shore distance (SD), maximum depth (MD), lake age (Age) and total phosphorus (TP) — and optical properties: mean extinction coefficient (KPAR; PAR: photosynthetically active radiation from 400 to 700 nm), and Secchi depth.
Site names as recorded in .
Underlined Years: years (2010 and 2011) for which photosynthetic pigment samples were able to be sampled. For example, in Goby Lake, a total of 9 jellyfish individuals were collected in 2010. Among these 9, j = 5 individuals were subsampled for subsequent analysis of photosynthetic pigments.
Adapted from from . For Tarzan, Jurassic, Secret, Ngerchaol, and Risong bay (site 2) MD, is not known precisely.
Based on , and .
From .
KPAR values in parentheses: Estimations based on the established relationship between TP and KPAR including all other sites (without jellyfish) for which TP and KPAR data were both available (from , for TP), with: ln (KPAR) = 1.10 (SE: ±0.20) ln (TP) – 4.95(SE: ±0.69) (p < 0.0001) (Supplementary Material S1).
Mean Secchi depth (± standard error) calculated from the monitoring program for CLM (n = 7), GLK (n = 30), BJLK (n = 30), OLO (n = 8), OTM (n = 48), RCA (n = 48), and Jurassic (n = 1). For Lake OLO, the Secchi depth reached the bottom (4 m) 5 times out of 8. See Supplementary Material S2.
Individual medusae were collected on site by snorkeling and transported in containers filled with ambient seawater to the Coral Reef Research Foundation (CRRF, Koror, Palau). During transport from the sampling site to CRRF (travel time of ≤1 h), individual jellyfish were kept in the dark in their containers. Once at CRRF, animals were maintained in clean, aerated seawater for 2–3 h in darkness prior to processing. Bell diameter and wet weight were measured using a graduated ruler (±0.1 cm) and a laboratory bench balance (±0.1 g) (data not shown here). For bell diameter measurements, individuals were gently placed on a laboratory dissecting pan (PVC). They were then weighed after being placed in a pre-tared weighing dish. Subsequently, fragments of oral arms were excised from each specimen for (i) photophysiological measurements (2010, 2011, 2013, and 2015) and (ii) pigment composition analyses (2010 and 2011).
2.2 Environmental parameters
2.2.1 Distance from the seashore, maximum depth, site age, and total phosphorus concentrations
The shortest distance from each site to the seashore (SD, m) was estimated using satellite imagery from Google Earth. The maximum depth was initially sourced from , and updated, if possible, during the field surveys. Bathymetry data were recorded at 2-s intervals using a Furuno single-beam echosounder paired with a Garmin GPS. The resulting bathymetric data were then analyzed, interpolated, and gridded in RStudio. The maximum depth (MD) for each site was then extracted from the gridded bathymetry model. For some other sites, MD is not precisely known, and only approximate estimates are reported in Table 1. Site ages were available for some sites as well (based on ). Total phosphorus (TP) concentrations were sourced from .
2.2.2 Irradiance
Vertical profiles of spectral irradiance (320–950 nm, Δλ = 3 nm) were measured using a RAMSES ACC-VIS hyperspectral radiometer (TriOS, Germany), from the surface to 10 m depth at 0.5–1 m intervals. Photosynthetically available radiation (PAR, 400–700 nm) was then calculated, and the extinction coefficient (KPAR, m-1) was estimated from the exponential decay of light intensity with depth, using the equation PAR(z) = PAR(z=0) × exp(−KPAR× Z), with Z = depth (m). In instances where the RAMSES ACC-VIS hyperspectral radiometer was unavailable (for specific years or sites), the diffuse attenuation coefficient for photosynthetic active radiation (KPAR) was estimated using the empirical relationship between total phosphorus (TP) and KPAR. This relationship was derived from data collected at all other sites where both TP and KPAR measurements were available, regardless of jellyfish presence. The relationship is described by the equation: ln (KPAR) = 1.10 (±0.20) × ln (TP) – 4.95 (±0.69) (pvalue < 0.0001; see Supplementary Material S1 for details). In addition, Secchi depth (m) is recorded as part of a monitoring program established by the CRRF (Palau) for 7 of the 11 sites studied here (Table 1, Supplementary Material S2).
2.3 Jellyfish
2.3.1 Photosynthetic pigment analysis
Jellyfish samples were collected for chlorophyll a and accessory pigments analysis in 2010 and 2011 (for a total of n = 44 jellyfish individual samples, Table 1). Tissue samples from the oral arm of jellyfish were rinsed in clean seawater, weighed, and stored at −20 °C in cryogenic tubes. For both years, photosynthetic pigment samples were processed and analyzed immediately upon return from the field expedition. Pigments were extracted in 3 mL of 95% ethanol at −20 °C for 24 h and quantified by HPLC (Agilent 1200) following , using certified pigment standards (DHI, Denmark).
2.3.2 Photosynthetic efficiency of endosymbionts
Photophysiological measurements of symbiotic dinoflagellates (Symbiodiniaceae, genus: Cladocopium) were performed on a total of 108 individuals in 2010, 2011, 2013 and 2015 (Table 1), using customized AquaPen-C AP-C 100 fluorometers (Photon Systems Instruments, Czech Republic). Jellyfish oral arm fragments were dark-adapted for ∼30 min, and chlorophyll a fluorescence transients were recorded using LEDs at 455, 470, 505, 530, 590, 620, or 630 nm and actinic light intensity of 60 µmol photons m-2 s-1. A white-light-equipped fluorometer (400–700 nm) was used to normalize across sensors.
Fluorescence parameters were extracted from chlorophyll fluorescence induction curves (OJIP transients), where each letter - O, J, I, and P - corresponds to a distinct step associated with the redox states of photosystems II. The OJIP transient curve included the following parameters (see
):
F0: Minimum fluorescence intensity at 50 µs; measured after dark adaptation
Fm: Maximum fluorescence intensity
Fv/Fm: Maximum quantum yield of PSII photochemistry (ϕpo), calculated as
ψ0: defined as 1 - (Fj - F0)/(Fm - F0), with Fj = fluorescence intensity at J-step (at 2 ms)
PIabs: Performance index on absorption basis, calculated as:
where RC is the flux trapped by reaction centers and ABS is the flux of photons absorbed by antenna pigments. PI
absis a well-established photosynthetic index that provides a more sensitive indicator of environmental stress conditions on PS II than quantum yield alone (
;
). More details are supplied in
Supplementary Material S3.
For each jellyfish sample, baseline fluorescence measurements at different illumination wavelengths, F
0(λ), were normalized by the maximum fluorescence value, F
0max, observed at λ = 455 or 470 nm. The normalized fluorescence, F
0n(λ), was calculated as:
F0n(λ) = F0(λ)/F0max
This normalization step was applied prior to any comparative or statistical analysis to account for potential biases arising from variations in jellyfish sub-sample sizes or year-to-year changes in instrument gain of the AquaPen-C AP-C 100 fluorometers.
2.4 Statistical analysis
All statistical analyses and figures were generated using RStudio (version 2025.05.1). Non-metric multidimensional scaling (NMDS) analysis was employed to investigate differences in photosynthetic pigment composition among jellyfish from different lakes and coves. Relationships between photosynthetic performance and shore distance were assessed using linear correlations, and non-parametric Pearson correlation coefficients.
3 Results
3.1 Environmental parameters
3.1.1 Distance from the shore, maximum depth and site age
In total, 11 sites were analyzed from 2010 to 2015 (Table 1); including six lakes, two open coves directly connected to the lagoon, and three semi-enclosed coves partially separated from the lagoon by a tidal sill. The minimum distance (SD) separating each site from the lagoon ranges from 0 m (Ngerchaol and Risong coves) to a maximum of 244 m (Clear Lake, CLM). Maximum depths (m) range from a minimum of ∼5–10 m (in coves) to a maximum of 38 m (Ngermeuangel, BJLK). Estimated lake ages reported in the literature range from ∼6,000 years for Lake Ongael (OLO) to ∼15,000 years for Lake Ngermeuangel (BJLK) (Table 1), and deeper lakes are presumed to have filled earlier, during the last deglaciation (). Thus, an approximate age estimate can be derived from the Age–Maximum Depth relationship - using GLK, BJLK, OLO and OTM Age and MD data–and expressed as: Age (years) ∼ 5000(SE: ±2000) + 200(SE: ±70) × MD (R2adj = 0.74, p = 0.092). Based on this relationship, Clear Lake (CLM; max depth: 34 m) filled ∼12,000 years ago, with shallower coves filling less than ∼5,000–10,000 years ago. These estimates, however, remain speculative owing to the age-depth model’s limited robustness.
3.1.2 Optical site characteristics and biogeochemical context
The extinction coefficients (KPAR, m-1) and mean Secchi depth (m), are given in Table 1. We used an established relationship between total phosphorus concentration (TP, µg·L-1) and KPAR–from data collected across all sites surveyed from 2010 to 2015 (irrespective of jellyfish presence) (Supplementary Material S1) - to estimate KPAR for sites lacking RAMSES ACC-VIS hyperspectral radiometer deployments (Table 1):
Given this approach, in lakes and coves inhabited by Mastigias jellyfish, photosynthetically active radiation (PAR) extinction coefficients (KPAR) exhibit variations ranging from 0.03 m-1 in Tarzan Lake, to 0.66 m-1 in Ongael Lake (OLO). As expected, mean Secchi depth and KPAR are correlated (Pearson correlation coefficient, r = −0.95, p-value = 0.0133), with Secchi depth ranging from a minimum of ∼3–4 m in the more turbid waters of Ongael Lake (OLO), to a maximum of ∼12–13 m in the clearer waters of Ngermeuangel Lake (BJLK). Regular monitoring of Secchi depth - conducted as part of the monitoring program established by CRRF (Palau) - revealed consistent year-to-year differences between sites (Supplementary Material S2). However, intra-site interannual variability is also observed, depending in particular on tidal and weather conditions. For example, in Goby Lake (GLK), Secchi depth readings dropped to 1.0–2.5 m in May–June 2015, following a lake overturning event in April–May 2015 (Supplementary Material S2).
3.1.3 Jellyfish photosynthetic pigment content
Irrespective of sites (lakes or coves) and sampling year (2010 and 2011), endosymbionts contained a common pool of accessory pigments. A total of 10 photosynthetic pigments were detected by HPLC (on n = 44 individual samples, Table 1). For each jellyfish sample, the relative photosynthetic accessory pigment composition of endosymbiotic zooxanthellae was expressed using a ratio calculated as follow:
The most abundant pigments were, in descending order of relative mass concentrations, peridin (Piratio = 2.71 ± 0.92), chlorophyll c2 (Piratio = 1.58 ± 0.45), monovinyl (MV) chlorophyll a (Piratio = 0.82 ± 0.07), diadinoxanthin (Piratio = 0.33 ± 0.05), chlorophyllide a (Piratio = 0.18 ± 0.07), and pheophytin a (Piratio = 0.15 ± 0.09). β-carotene, and the sum of other xanthophylls (zeaxanthin, diatoxanthin, and violaxanthin) were also detected, but in smaller quantities (Piratio < 0.03).
Alongside this common pool of accessory pigments, NMDS (Non-metric MultiDimensional Scaling) analysis revealed differences in Piratios, both between years, and within and between sites (Figure 2; NMDS stress = 0.042). Axis 1 of the NMDS shows lower MV chlorophyll a Piratios, and higher peridinin, chlorophyllide a, chlrophyll c2, and pheoporbide Piratios, in some of the individuals sampled in 2011, compared to those sampled in 2010, in lake OTM, lake BJLK, and Ngelchaol Cove. Similarly, NMDS axis 2 shows that some individuals collected in 2011 in lake GLK and OTM had higher diadinoxanthin, and other xanthophylls Piratios. Conversely, individuals collected from lake OLO in 2010 had higher β-carotene and pheophytin a Piratios.
FIGURE 2
3.1.4 Photosynthetic dynamics of the photosystem II (PSII)
Normalized baseline fluorescence (F0n) in all jellyfish samples exhibited peak values at λ = 455 nm and 470 nm, then declined along the wavelength gradient to a minimum at λ = 590 nm (Figure 3). Additionally, when data from all years were pooled, F0n at all wavelengths exhibited no significant variation along the shore distance gradient for individuals sampled in lakes and coves (Figure 3).
FIGURE 3
Conversely, both the quantum yield (Fv/Fm) and the photosynthetic performance index (PIabs) were negatively correlated with increasing shore distance across all wavelength bands, except at 590 nm (yellow) and 620 nm (orange) (Figures 4, 5).
FIGURE 4

Relationships between maximum quantum yield (FV/Fm) and shore distance (m) under different wavelength illuminations. Symbols represent means ± standard error (vertical bars). Pearson correlation coefficients were calculated across all samples. Only significant correlations are highlighted with a dotted line (fitted to the data using the “geom_smooth ()” function from the R ggplot2 package). Lake names follow
FIGURE 5

Relationships between the relative performance index on an absorption basis (PIabs) and shore distance (m) under different wavelength illuminations. Symbols represent means ± SE (vertical bars). Pearson correlation coefficients were calculated across all samples. Only significant correlations are highlighted with a dotted line (fitted to the data using the “geom_smooth ()” function from the R ggplot2 package). Lake names follow
Plotting the residuals from the regressions of FV/Fm and PIabs versus shore distance (Figures 4, 5) against Maximum depth (a proxy of lake age) revealed no significant relationship (Figure 6; Pearson correlation coefficient, p > 0.05).
FIGURE 6

Relationships between the residuals of (A) the maximum quantum yield (FV/Fm), and (B) the relative performance index on an absorption basis (PIabs)—both calculated from the regression of FV/Fm or PIabs against shore distance (m) under white light illumination (see Figures 4, 5)—and maximum depth (m). Sites with unknown maximum depth (Tarzan, Jurassic, Secret, and two unnamed coves: Ngerchaol and Risong Bay, site 2; see Table 1) were excluded from this analysis.
4 Discussion
Many marine organisms form endosymbiosis with autotrophic dinoflagellates, so called zooxanthellae. Zooxanthella could be important for supporting nutrition but also for life cycle dynamics such as for example facilitating strobilation in some jellyfish species (
Pertaining to marine lakes, their underwater light environment is more similar to typical freshwater lakes than to open marine sites. Hence, it seems that for the jellyfish living in marine lakes the pigments in the endosymbionts and the underwater light climate in the lakes do not match as well as for jellyfish living in lagoons or open ocean areas. Our multi-annual measurements of the photosynthetic index, or the photosynthetic performance of PSII, of jellyfish endosymbionts supports this mismatch. The photosynthetic performance of Mastigias in the lakes was negatively correlated with the distance from the surrounding lagoon. One would expect a larger input of humic substances affecting light climate towards a higher contribution of red wavelengths with larger terrestrial catchment areas indicated by larger distance from the lagoon.
While Mastigias evolved differently in lakes towards smaller sizes, smaller tentacles, and complex behaviors such as diel vertical migration (
This mismatch between the spectral qualities of light in marine lakes and the pigment composition of Cladocopium sp. may have far-reaching consequences for the energy budgets and overall fitness of Mastigias inhabiting these environments. Reduced photosynthetic efficiency, as demonstrated by the lower PSII performance in lake-dwelling jellyfish, implies lower net carbon fixation by the symbionts, which in turn may reduce the autotrophic contribution to the host’s metabolic demands. This metabolic shortfall might necessitate compensatory strategies such as increased heterotrophic feeding, a behavior observed in several mixotrophic cnidarians (
Interestingly, the ecological flexibility of the host appears to exceed that of its symbiont. As noted, Mastigias in marine lakes have evolved substantial morphological and behavioral adaptations to cope with their isolated and often stratified habitats (
One might also consider whether the persistence of Cladocopium spp. in marine lake jellyfish is a consequence of historical constraints rather than present-day optimality. Given the isolated nature of marine lakes, opportunities for acquiring novel endosymbionts may be limited. Horizontal transmission of symbionts is common in Mastigias spp. (
Interestingly a recent study found evidence for significant adaptive differentiation in Cladocopium associated with corals across reef environments in Palau (
In summary, the reduced photosynthetic performance of jellyfish endosymbionts in marine lakes likely reflects a fundamental ecological mismatch driven by the light environment. While Mastigias papua spp. have adapted morphologically and behaviorally to these habitats, their dinoflagellate endosymbionts appear to face significant physiological constraints. This highlights the critical role of symbiont availability and compatibility in shaping host success and opens avenues for exploring co-evolution, symbiont plasticity, and the ecological limits of mutualistic associations in novel environments.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The manuscript presents research on animals that do not require ethical approval for their study.
Author contributions
PP: Writing – original draft, Conceptualization, Formal Analysis, Investigation, Supervision, Validation, Data curation, Funding acquisition, Project administration, Resources, Visualization, Writing – review and editing. HS: Writing – original draft, Conceptualization, Formal Analysis, Funding acquisition, Investigation, Project administration, Supervision, Validation, Data curation, Resources, Visualization, Writing – review and editing. JY-TH: Writing – original draft, Conceptualization, Formal Analysis, Investigation, Validation, Data curation, Visualization. MS: Writing – original draft, Conceptualization, Formal Analysis, Investigation, Validation, Data curation, Project administration, Resources, Supervision, Visualization, Writing – review and editing. SB: Writing – original draft, Conceptualization, Formal Analysis, Investigation, Validation, Data curation, Visualization. TS: Writing – original draft, Formal Analysis, Investigation, Methodology, Data curation, Visualization. SP: Writing – original draft, Investigation. GU: Writing – original draft, Investigation.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This research received base funding from the Ludwig-Maximilians-Universität (LMU, München, Germany) and the University of Bretagne Occidentale (Institute Universitaire Européen de la Mer, IUEM, UBO).
Acknowledgments
We thank Pat Colin, Lori J. Bell, Michael N. Dawson, Mathew Mesubed, and Emilio Basilius for their support and logistical assistance, both at the Coral Reef Research Foundation (CRRF) and in the field. We particularly wish to honour the memory of Emilio Basilius (1953–2019), whose contributions greatly supported field research in Palau’s marine lakes. Finally, we thank the Palau National Bureau of Marine Resources and the Koror State Government for authorizing our research in Palau.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphbi.2026.1670101/full#supplementary-material
References
1
AraiY.GotohR. O.YokoyamaJ.SatoC.OkuizumiK.HanzawaN. (2017). Phylogenetic relationships and morphological variations of upside-down jellyfishes, Cassiopea spp. inhabiting Palau Islands. Biogeography19, 133–141. 10.11358/biogeo.19.133
2
ArmstrongK. C.LippertM.HansonE.NestorV.CornwellB.WalkerN. S.et al (2024). Fine-scale geographic variation of cladocopium in Acropora hyacinthus across the Palauan archipelago. Ecol. Evol.14, e70650. 10.1002/ece3.70650
3
BakerA. C. (2003). Flexibility and specificity in coral–algal symbiosis: diversity, ecology, and biogeography of Symbiodinium. Ann. Rev. Ecol., Evol., Syst.34 (1), 661–689. 10.1146/annurev.ecolsys.34.011802.132417
4
CeustersN.ValckeR.FransM.ClaesJ. E.Van den EndeW.CeustersJ. (2019). Performance index and PSII connectivity under drought and contrasting light regimes in the CAM orchid Phalaenopsis. Front. Plant Sci.10, 1–15. 10.3389/fpls.2019.01012
5
ColinP. L. (2009). “Marine environments of Palau,” in Indo-Pacific. Press, San Diego, CA, USA.
6
DavyS. K.AllemandD.WeisV. M. (2012). Cell biology of cnidarian–dinoflagellate symbiosis. Microb. Mol. Biol. Rev.76 (2), 229–261. 10.1128/mmbr.05014-11
7
DawsonM. N.HamnerW. M. (2003). Geographic variation and behavioral evolution in marine plankton: the case of Mastigias (Scyphozoa, Rhizostomeae). Mar. Biol.143 (6), 1161–1174. 10.1007/s00227-003-1130-6
8
DawsonM. N.HamnerW. M. (2005). Rapid evolutionary radiation of marine zooplankton in peripheral environments. Proc. Natl. Acad. Sci. U. S. A.102 (26), 9235–9240. 10.1073/pnas.0503635102
9
DjeghriN.PondavenP.StiborH.DawsonM. N. (2019). Review of the diversity, traits, and ecology of zooxanthellate jellyfishes. Mar. Biol.166 (11), 147. 10.1007/s00227-019-3581-6
10
DjeghriN.PondavenP.StockenreiterM.BehlS.HuangJ. Y. T.HansenT.et al (2020). Isotopic and elemental compositions reveal density-dependent nutrition pathways in a population of mixotrophic jellyfish. Ecosphere11 (11), e03295. 10.1002/ecs2.3295
11
ErikssonS. Y.MatzM. V.VizeP. D.RosserN. L. (2025). Whole Genome sequencing reveals genetic differences between Symbiodiniaceae populations among reproductively and geographically isolated Acropora colonies in Western Australia. Ecol. Evol.15 (1), e70771. 10.1002/ece3.70771
12
FalkowskiP. G.LaRocheJ. (1991). Acclimation to spectral irradiance in algae. J. Phycol.27 (1), 8–14. 10.1111/j.0022-3646.1991.00008.x
13
FalkowskiP. G.KatzM. E.KnollA. H.QuiggA.RavenJ. A.SchofieldO.et al (2004). The evolution of modern eukaryotic phytoplankton. Science305 (5682), 354–360. 10.1126/science.1095964
14
FederM. E.Mitchell-OldsT. (2003). Evolutionary and ecological functional genomics. Nat. Rev. Genet.4 (8), 651–657. 10.1038/nrg1128
15
FerreiraM. R.ClearyD. F.BatN. K.PolóniaA. R.GomesN. C. (2022). Assessing the core microbial symbionts of jellyfish in Indonesian and Vietnamese marine lakes. J. Mar. Biol. Ass. U. K.102 (7), 486–495. 10.1017/S0025315422000777
16
GoodsonM. S.WhiteheadL. F.DouglasA. E. (2001). Symbiotic dinoflagellates in marine Cnidaria: diversity and function. Hydrobiologia461 (1), 79–82. 10.1023/A:1012766609448
17
GoodwinS.McPhersonJ. D.McCombieW. R. (2016). Coming of age: ten years of next-generation sequencing technologies. Nat. Rev. Genet.17 (6), 333–351. 10.1038/nrg.2016.49
18
HamnerW. M.HamnerP. P. (1998). Stratified marine lakes of Palau (Western Caroline islands). Phys. Geogr.19 (3), 175–220. 10.1080/02723646.1998.10642647
19
HamnerW. M.GilmerR. W.HamnerP. P. (1982). The physical, chemical, and biological dynamics of a stratified, saline lake in Palau. Limnol. Oceanogr.27 (4), 896–909. 10.4319/lo.1982.27.5.0896
20
KristensenT. N.KetolaT.KronholmI. (2018). Adaptation to environmental stress at different timescales. Ann. N. Y. Acad. Sci.1476 (1), 5–22. 10.1111/nyas.13974
21
LaJeunesseT. C.ParkinsonJ. E.GabrielsonP. W.JeongH. J.ReimerJ. D.VoolstraC. R.et al (2018). Systematic revision of Symbiodiniaceae highlights the antiquity and diversity of coral endosymbionts. Curr. Biol.28 (16), 2570–2580.e6. 10.1016/j.cub.2018.07.008
22
LewinH. A.RobinsonG. E.KressW. J.BakerW. J.CoddingtonJ.CrandallK. A.et al (2018). Earth BioGenome Project: sequencing life for the future of life. Proc. Natl. Acad. Sci. USA.115 (17), 4325–4333. 10.1073/pnas.1720115115
23
MatthewsJ. L.BlackallL. L.van OppenM. J. H.BournD. G.SuggettD. J.CampE. F.et al (2023). Symbiodiniaceae photophysiology and stress resilience is enhanced by microbial associations. Sci. Rep.13, 20724. 10.1038/s41598-023-47900-3
24
McFall-NgaiM.HadfieldM. G.BoschT. C. G.CareyH. V.Domazet-LošoT.DouglasA. E.et al (2013). Animals in a bacterial world: a new imperative for the life sciences. Proc. Natl. Acad. Sci. USA.110 (9), 3229–3236. 10.1073/pnas.1218525110
25
MuscatineL. (1990). “The role of symbiotic algae in carbon and energy flux in reef corals,” in Ecosystems of the World 25: coral reefs. Editor DubinskyZ. (Amsterdam: Elsevier), 75–87.
26
PageT. M.LawleyJ. W. (2022). The next generation is here: a review of transcriptomic approaches in marine ecology. Front. Mar. Sci.9, 757921. 10.3389/fmars.2022.757921
27
Rodríguez Del RíoÁ.Giner-LamiaJ.CantalapiedraC. P.BotasJ.DengZ.Hernández-PlazaA.et al (2023). Functional and evolutionary significance of unknown genes from uncultivated taxa. Nature626, 377–384. 10.1038/s41586-023-06609-1
28
SavolainenO.LascouxM.MeriläJ. (2013). Ecological genomics of local adaptation. Nat. Rev. Genet.14 (11), 807–820. 10.1038/nrg3522
29
SelosseM. A.CharpinM.NotF. (2017). Mixotrophy everywhere on land and in water: the “grand écart” hypothesis. Ecol. Lett.20 (2), 246–263. 10.1111/ele.12714
30
StanleyJ. G. D.Van De SchootbruggeB. (2009). “The evolution of the coral–algal symbiosis,” in Coral bleaching: patterns, processes, causes and consequences. Berlin, Heidelberg: Springer Berlin Heidelberg, 7–19.
31
StinchcombeJ. R.HoekstraH. E. (2008). Combining population genomics and quantitative genetics: finding the genes underlying ecologically important traits. Heredity100 (2), 158–170. 10.1038/sj.hdy.6800937
32
StompM.HuismanJ.StalL. J.MatthijsH. C. P. (2007). Colorful niches of phototrophic microorganisms shaped by vibrations of the water column. ISME J.1 (4), 271–282. 10.1038/ismej.2007.59
33
StrasserR. J.SrivastavaA.Tsimilli-MichaelM. (2000). The fluorescence transient as a tool to characterize and screen photosynthetic samples. Probing Photosynthesis Mechanisms, Regulation Adaptation25, 445–483.
34
SugiuraY. (1964). On the life-history of Rhizostome medusae: II. Indispensability of zooxanthellae for strobilation in Mastigias papua. Embryologia8 (3), 223–233. 10.1111/j.1440-169X.1964.tb00200.x
35
SummererM.SonntagB.SommarugaR. (2008). Ciliate-symbiont specifity of freshwater endosymbiotic Chlorella. J. Phycol.44, 77–84. 10.1111/j.1529-8817.2007.00455.x
36
Van HeukelemL.ThomasC. S. (2001). Computer-assisted high-performance liquid chromatography method development with applications to the isolation and analysis of phytoplankton pigments. J. Chromatogr. A910 (1), 31–49. 10.1016/S0378-4347(00)00603-4
37
Vega de LunaF.DangK. V.CardolM.RobertyS.CardolP. (2019). Photosynthetic capacity of the endosymbiotic dinoflagellate Cladocopium sp. is preserved during digestion of its jellyfish host Mastigias papua by the anemone Entacmaea medusivora. FEMS Microbiol. Ecol.1 (10), fiz141. 10.1093/femsec/fiz141
38
VennA. A.LoramJ. E.DouglasA. E. (2008). Photosynthetic symbioses in animals. J. Exp. Bot.59 (5), 1069–1080. 10.1093/jxb/erm328
39
West-EberhardM. J. (2003). Developmental plasticity and evolution. Oxford University Press.
40
WollrabS.PondavenP.BehlS.BekerB.StiborH. (2021). Differences in size distribution of marine phytoplankton in presence versus absence of jellyfish support theoretical predictions on top-down control patterns along alternative energy pathways. Mar. Biol.167, 9. 10.1007/s00227-019-3621-2
41
ZuoG.MeiW.FengN.ZhengD. (2025). Photosynthetic performance index (PIabs) and malondialdehyde (MDA) content determine rice biomass under combined salt stress and prohexadione-calcium treatment. BMC Plant Biol.25, 823. 10.1186/s12870-025-06826-x
Summary
Keywords
endosymbiotic dinoflagellates, environmental gradient analysis, marine lakes, Mastigias papua, Palau
Citation
Pondaven P, Stibor H, Yu-Tung Huang J, Stockenreiter M, Behl S, Stieglitz TC, Patris S and Ucharm G (2026) Photophysiological performance of zooxanthellate endosymbionts (Symbiodiniaceae) in the golden jellyfish Mastigias papua across a natural environmental gradient in marine lakes of Palau (Micronesia). Front. Photobiol. 4:1670101. doi: 10.3389/fphbi.2026.1670101
Received
21 July 2025
Revised
18 December 2025
Accepted
29 January 2026
Published
05 March 2026
Volume
4 - 2026
Edited by
Caroline Eve Dubé, Laval University, Canada
Reviewed by
Todd LaJeunesse, The Pennsylvania State University (PSU), United States
Ruben X. G. Silva, University of Aveiro, Portugal
Updates

Check for updates
Copyright
© 2026 Pondaven, Stibor, Yu-Tung Huang, Stockenreiter, Behl, Stieglitz, Patris and Ucharm.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Philippe Pondaven, philippe.pondaven@univ-brest.fr; Herwig Stibor, stibor@zi.biologie.uni-muenchen.de
ORCID: Sharon Patris, orcid.org/0000-0001-8860-560X; Philippe Pondaven, orcid.org/0000-0002-3427-7767; Herwig Stibor, orcid.org/0000-0001-5918-6020; Maria Stockenreiter, orcid.org/0000-0001-7380-071X; Thomas C. Stieglitz, orcid.org/0000-0002-2946-391X
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.