Blog

Apr 08, 2023

Photosynthetic usable energy explains vertical patterns of biodiversity in zooxanthellate corals

Scientific Reports volume 12,

Scientific Reports volume 12, Article number: 20821 (2022) Cite this article

1188 Accesses

2 Citations

78 Altmetric

Metrics details

The biodiversity in coral reef ecosystems is distributed heterogeneously across spatial and temporal scales, being commonly influenced by biogeographic factors, habitat area and disturbance frequency. A potential association between gradients of usable energy and biodiversity patterns has received little empirical support in these ecosystems. Here, we analyzed the productivity and biodiversity variation over depth gradients in symbiotic coral communities, whose members rely on the energy translocated by photosynthetic algal symbionts (zooxanthellae). Using a mechanistic model we explored the association between the depth-dependent variation in photosynthetic usable energy to corals and gradients of species diversity, comparing reefs with contrasting water clarity and biodiversity patterns across global hotspots of marine biodiversity. The productivity-biodiversity model explained between 64 and 95% of the depth-related variation in coral species richness, indicating that much of the variation in species richness with depth is driven by changes in the fractional contribution of photosynthetically fixed energy by the zooxanthellae. These results suggest a fundamental role of solar energy availability and photosynthetic production in explaining global-scale patterns of coral biodiversity and community structure along depth gradients. Accordingly, the maintenance of water optical quality in coral reefs is fundamental to protect coral biodiversity and prevent reef degradation.

Sunlight is the major source of energy for virtually the entire biochemical production of organic matter on Earth1. Among the different factors affecting primary producers both in terrestrial and aquatic ecosystems, sunlight is perhaps the most spatially and temporally heterogeneous, resulting in significant impacts on the energy transfer across trophic levels2,3,4. The variation in the supply of usable energy derived from primary production plays an important role on the spatial variation of species diversity in ecological communities, facilitating a decreased risk of species extinction from demographic and environmental stochasticity in highly productive environments5. Although a positive productivity-biodiversity relationship prevails in terrestrial and aquatic ecosystems6, significant debate remains over the strength of the linkage between these two features and whether it is predictable, considering the observed variations in the form of productivity-biodiversity relationship in some communities (e.g., positive, negative, unimodal and neutral)7,8,9,10. Since the amount of usable energy by a group of organisms in a community is difficult to measure directly, part of this variation has been associated with the surrogates used as roughly indicative of productivity and the confounding effects of environmental stressors, typically more accentuated in highly productive environments6,7.

Scleractinian corals are metazoans with the ability to calcify, responsible for building the most diverse and productive marine ecosystem from which humans obtain essential socio-economic benefits. The ecological and evolutionary success of these animals in oligotrophic environments since the Late Triassic is attributed to the nutritional endosymbiosis with photosynthetic unicellular algae (zooxanthellae)11, which confers on them the capacity to use sunlight as the major source of energy increasing their rates of calcification12. Due to the exponential attenuation of sunlight with depth13, the vertical distribution of symbiotic corals occurs along steep gradients governed by light in small spatial scales, while other limiting resources and physical factors (e.g., nitrogen, oxygen, temperature) vary less or remain nearly constant14,15. As benthic, sessile organisms prevalent in coastal waters, corals are also exposed to significant changes in light intensity in response to the optical properties of the water column (i.e., vertical attenuation coefficient for downwelling irradiance, Kd). This variation is determined by precipitation patterns and seaward fluxes of terrestrial nutrients and sediments13,16, having important implications on the primary production and energy balance of coral communities3,4. However, associations between gradients of productivity and biodiversity in coral communities have not yet been demonstrated, and patterns of coral biodiversity are primarily attributed to biogeographic factors, habitat area and disturbance frequency10,17,18,19,20,21,22.

The obligatory nature of the endosymbiosis between a primary producer and an animal, together with the vertical distribution of corals along strong gradients of sunlight energy, make coral communities particularly interesting for exploring the relationship between primary productivity and consumer-species diversity. Using a mechanistic model based on underwater optics and physiological principles, we estimated the variation in photosynthetic usable energy supplied by the symbiotic algae to coral hosts and examined its association with patterns of coral species diversity across depth gradients. We further tested the model using published datasets from reefs with contrasting evolutionary histories, environmental conditions, and patterns of coral species diversity, encompassing the three major hotspots of global marine biodiversity. Our analysis shows that much of the variation in coral diversity with depth is explained by changes in photosynthetic usable energy, providing support to a mechanistic link between productivity and biodiversity in coral communities.

The diurnal patterns of change of the photosystem II (PSII) effective quantum yield (ΔF/Fm’) measured with pulse amplitude modulated (PAM) fluorometry, followed an inverse association with light availability (Fig. 1a). The diurnal oscillation in ΔF/Fm’ results from the induction of photochemical (qP) and non-photochemical quenching (qN or NPQ) processes, which occur simultaneously in the photosynthetic apparatus in order to transform part of the solar energy absorbed into organic carbon (qP) and dissipate the energy absorbed in excess as heat (NPQ)23,24,25,26. On day one in Fig. 1a, a complete recovery of the maximum quantum yield of PSII (Fv/Fm) was observed by the end of the day when corals were exposed to similar light conditions (control) relative to previous days (9 mol quanta m−2 d−1 with a peak irradiance of 454 μmol quanta m−2 s−1). This observation indicates that no accumulation of photodamage occurred in these corals at the end of the day as NPQ was fully relaxed. On day two, light exposure was increased by nearly three times relative to control (26 mol quanta m−2 d−1 with a peak irradiance of 1130 μmol quanta m−2 s−1) and the amplitude of the ΔF/Fm’ oscillation was greater, reaching a minimum value at noon that was lower with respect to the previous day. This finding suggests a greater induction of NPQ processes to dissipate the increased amount of absorbed light energy. Additionally, the incomplete recovery of Fv/Fm observed at dusk (0.529) compared to former dawn levels (0.643) suggests the accumulation of photodamage in the photosynthetic apparatus, potentially associated with an incomplete relaxation of the photoinhibitory quenching (qI) component of NPQ which requires de novo synthesis of PSII proteins (e.g., D1) for its relaxation27,28,29. The accumulation of photodamage as light exposure was increased can be explained by the rate of damage to PSII exceeding that of repair. On day three, corals were exposed to nearly half of the light levels on day one (5 mol quanta m−2 d−1 with a peak of irradiance of 237 μmol quanta m−2 s−1), which resulted in reduced oscillation in ΔF/Fm’ and recovery in Fv/Fm at the end of the day (0.625) exceeding formed dawn levels (0.523). These responses to irradiance reduction indicate a lower induction of NPQ processes and a greater relaxation of qI at dusk as a result of greater rates of PSII repair compared to those of photodamage27,28,29 (Fig. 1a). The observed patterns in ΔF/Fm’ oscillation and daily Fv/Fm variation are associated with down-regulation responses of algal symbionts to protect their photosynthetic apparatus and optimize photosynthetic performance under variable light environments. Furthermore, they indicate that light exposure plays a critical role in the balance between the rates of PSII damage and subsequent repair, for which de novo synthesis of proteins for the reassembly of reaction centers is required.

Principles and components of the coral productivity-biodiversity model. (a) Diurnal variation of ΔF/Fm’ in symbionts of Porites astreoides measured experimentally during a control-reference day (first day), and during days with increased (second day) and reduced light exposure (third day), emphasizing the consequences on the Fv/Fm recovery at the end of the diurnal cycle (arrows). The yellow shaded area represents the diurnal variation of irradiance (E). (b) Changes in the PSII t1/2 in response to contrasting light exposure (in % of sea surface), fitting the data to a power function (continuous line). The dotted line indicates the relative energy cost of repairing photodamaged PSII by the symbiotic algae (Ca) with light exposure, expected to be inversely proportional to the PSII t1/2. (c) Theoretical behavior of coral holobiont energetic performance (continuous line) in relation to Qm (discontinuous line) across a gradient of irradiance, highlighting the predicted effects of light stress (LS) and light limitation (LL) on coral holobiont energetic performance. (d) Schematic representation of the relative changes in light-induced processes affecting the energy balance of coral holobionts along a depth gradient. PUES: photosynthetic usable energy supply; Ca energy costs of repair from photodamage for the algae; and Pg: gross productivity.

To better understand the relationship between light exposure and the rates of PSII damage, a series of experiments were performed using chloramphenicol (CAP), an inhibitor of chloroplast protein synthesis that suppresses repair of damaged PSII reaction centers30. The difference between controls and samples exposed to CAP allowed us to determine the PSII half-time (t1/2) as a function of light exposure. The association found between PSII t1/2 and irradiance is non-linear and can be described using a power-function model (R2 = 0.97, p < 0.01) (Fig. 1b). Similar patterns have been observed in other primary producers31,32, suggesting that the non-linearity between PSII t1/2 and irradiance may be driven by common mechanisms across primary producers. A conceptual analysis of the changes in the maximum excitation pressure over PSII (Qm) as a function of irradiance indicates that both the shortage as well as the excess of light compromise the energetic performance of coral holobionts (Fig. 1c). This compromise results from two contrasting effects: a negligible photosynthetic activity of the zooxanthellae under light-limiting conditions (low Qm), and increased costs of repair from photodamage at high irradiance levels (high Qm)28,33.

The non-linearity between downwelling irradiance, photosynthetic production and energy expenditure in photorepair, leads to a unimodal, hump-shaped pattern of the photosynthetic usable energy to corals with depth (Fig. 1d). This pattern is associated with an energetic imbalance that results from the hyperbolic response of photosynthesis to light (i.e., saturation of photosynthesis above the light saturation point, Ek), and the linear relationship between light availability and light absorption34,35,36. The decline in usable energy toward both ends of the depth-mediated light gradient is explained by two different processes: (1) in shallow high-light environments, the energy expenditure in photorepair increases proportionally to the amount of light above Ek while photosynthesis is saturated (Pmax); and (2) in deep low-light environments, the light attenuation results in a gradual decline of the photosynthetic activity. At a particular intermediate depth, where the absorbed irradiance is close to that required for maximizing photosynthesis (Ek), the supply of photosynthetic usable energy and the holobiont energetic performance are predicted to reach maximum potential. The level of irradiance that is required to saturate the photosynthetic apparatus (Ek) determines the threshold for the escalation of the amount of solar energy absorbed in excess relative to the maximum photosynthetic capacity (Pmax) of corals, which remains constant above Ek.

The coral productivity-biodiversity model explained most of the variation in coral species number with depth in all reefs tested (Fig. 2, Supplementary Tables S1, S2). The explanatory power of the model was consistently high, between 80 and 95%, in sites with unimodal relationships between species richness and depth (e.g., atolls from the Chagos Archipelago in the Indian Ocean37, Discovery Bay in the Caribbean38, and Kimbe Bay in the Indo-Pacific10). The model also captured the overall pattern of change in species number with depth in sites with monotonic patterns, both positive (Gulf of Eilat, Red Sea39) and negative (Pulau Hantu, Singapore40). In these sites, the model respectively explained 64% and 75% of species number variation, and in the site with a positive, monotonic relationship, it predicted a reduction in species number toward non-sampled depths below ~ 40 m, suggesting that these monotonic patterns may represent only one-half of the unimodal productivity-biodiversity relationship depending on whether the sampled depth range was located above or below the mode (Fig. 2).

Coral species richness variation along depth gradients in reefs from major hotspots of biodiversity. The observed distribution of species richness (squares) is significantly explained by the coral productivity-biodiversity model at all sites. Continuous lines represent the trend in mean values and the discontinuous lines 95% confidence intervals. The goodness-of-fit (R2) and the statistical significance of the model (p-value), together with the local Kd are indicated for each site. Species richness was projected to depths at which the light intensity was estimated to be 2% of the incident irradiance at sea surface.

The Kd's in sites with unimodal relationships were consistently low (between 0.06 and 0.082 m−1)10,38,41. In contrast, in sites with positive and negative monotonic relationships, the Kd's were respectively the lowest and highest of all sites considered in our analysis (0.047 m−1 and 0.303 m−1)40,41 (Supplementary Table S2). The deepest record of maximum coral species richness by location was associated with the water optical properties determined by the local Kd. A non-linear association between both parameters was observed and described using a power function (R2 = 0.98, p < 0.01, Fig. 3a), indicating a fast reduction of the depth of maximum coral richness with increasing Kd (Fig. 3a). The overall variation in the relative number of coral species by location along a similar light intensity gradient normalized to the percentage of surface irradiance mediated by the local Kd, followed a unimodal, bell-shaped pattern with an overall reduction in the relative number of species at both ends of the light intensity gradient (Fig. 3b).

Relationship between the distribution of coral species richness and the water optical properties. (a) Vertical attenuation coefficient for downwelling irradiance (Kd) versus the depth of maximum diversity by location. A power function depicted in the figure was used to fit a linear regression to the data. (b) Variation in coral species richness along a light intensity gradient determined by the local Kd. Relative light exposure is reported in % of sea surface. The overall data was fitted using a gaussian function (continuous line) with 95% confidence intervals (discontinuous lines).

The repair of the photodamaged PSII reaction centers is an essential component of the photoprotection capacity of the zooxanthellae to maintain an optimal photosynthetic performance of coral holobionts in contrasting light environments35,36. Although photodamage and photoprotection are processes intensively studied on intact corals and/or their isolated symbionts during last decades23,24,25,26,42,43,44, the understanding of their implications on coral energetic performance at the community and ecological levels is still limited. The induction of these processes can be inferred from variations in the effective (ΔF/Fm’) and maximum (Fv/Fm) quantum yield of PSII of algal symbionts in hospite, primarily studied by quenching analysis of the chlorophyll a fluorescence signal28,35. Two processes, photochemical (qP) and non-photochemical (qN or NPQ) regulate the quenching of the chlorophyll a fluorescence signal. The first pathway, qP, is directly associated with the photosynthetic conversion of sunlight energy to chemical energy, whereas qN or NPQ are two parameters developed in quenching analysis to quantify the photoprotective capacity of the photosynthetic apparatus to dissipate as heat the excessive excitation energy. Three main NPQ components with contrasting relaxation kinetics have been characterized: the energy-dependent quenching (qE), the quickest to relax associated with the operation of the xanthophyll cycle; the state-transition quenching (qT) associated with the redistribution of excitation energy between both photosystems, PSII and PSI, by physical modulation of their absorption cross-section; and the photoinhibition quenching (qI), the slowest to relax and related to the reversible photoinhibition of PSII after de novo synthesis of proteins. This last component reflects the capacity to repair inactive/photodamaged PSII27,28.

In this study, we parameterized the costs of repair of the photosynthetic apparatus on in hospite zooxanthellae as a function of irradiance. This parameterization was derived from the experimental determination of the variation of PSII half-time (t1/2) with light exposure, after suppressing the PSII repair capacity in the zooxanthellae. Differences in Fv/Fm at the beginning and end of the diurnal cycle when light exposure is increased or reduced, reflect changes in the accumulation of photoinactivated PSII reaction centers. These differences arise from disparities between the rates of photodamage and repair of the photosynthetic apparatus of in hospite zooxanthellae in particular light environments. The kinetics of the PSII t1/2 indicates that the rate of PSII damage is proportional to light exposure, and thus the energy cost of repair the damaged PSII reaction centers for optimizing the algal photosynthetic performance should also be variable and mediated by light exposure. This assumption is supported by the evidence that a constant synthesis and replacement of proteins (e.g., D1) is required during the photosynthetic activity for the re-assembly of PSII reaction centers, and that this protein turnover can be the largest single contributor to the costs of maintenance in primary producers45,46. The parameterization of the energy costs to maintain photosynthesis of in hospite zooxanthellae along depth-mediated light gradients is a key factor taken into account in our model because of the potential limiting effect on the translocation of photosynthetic usable energy to their coral hosts, expected to have relevant ecological implications.

Bearing in mind the consistency of the productivity-biodiversity model tested here, despite local differences in reef geological history, environmental conditions and patterns of diversity, our analysis indicates that much of the variation in coral species richness along depth gradients is driven by changes in the fractional contribution of photosynthetically fixed energy by the symbiotic algae to their coral hosts. These findings support the occurrence of a productivity-biodiversity association in reef-building coral communities, similar to that widely acknowledged in communities across terrestrial environments5,6,20. Furthermore, this linkage highlights the fundamental role of zooxanthellae productivity in the structure of coral reef communities. The lack of prior support for a productivity-biodiversity relationship in coral communities10 arise from unclear definitions of the actual energy physiologically available to the coral animal (e.g., solar energy), ignoring key physiological processes that constrain the photosynthetic activity of its algal symbionts. This contradiction illustrates the difficulty of actually measuring gradients of usable energy available to organisms in studies of productivity-biodiversity associations, and of choosing the energy-related variable that best explains species richness variation according to the system studied5,6,7.

Although coral richness followed a positive relationship with the photosynthetic usable energy supply, the productivity-biodiversity association is described by a unimodal, humpbacked curve with depth whose shape and mode localization are highly influenced by the water optical properties (i.e., local Kd). The overall humpbacked shape of the curve results from the non-linearity between depth, light availability, and the photosynthetic activity of the symbiotic algae. The increasing and decreasing phases of the humpbacked species-richness curve with depth result from two different processes that can limit the energetic output of the algae and, thus, the coral holobiont performance. In deep, low-light environments (the increasing phase of the richness curve), light deprivation results in a deficit of organic carbon derived from photosynthesis that can be translocated to the coral animal. In contrast, in shallow, high-light environments (the decreasing phase of the curve), the increased costs of maintenance of the photosynthetic activity once photosynthesis is fully saturated, limit the amount of photosynthetic usable energy that can be translocated to the coral host. This condition is probably coupled with a strong selective pressure exerted by intense light on the coral community in shallow habitats, including high levels of UVR47. At intermediate irradiance, the energetic output of the symbionts and the coral holobiont performance are predicted to reach their maximum potential, which can lead to reduced rates of species extinction and increased biodiversity5. The contrasting patterns observed in the productivity-biodiversity curves by location and the relationship between the depth of maximum richness and the local Kd illustrate the essential role of the water optical properties on the spatial organization and depth distribution of symbiotic coral communities. Thus, the relationship of the optical properties of the water is not only limited to the lower depth distribution of symbiotic corals41, but also, to patterns of coral biodiversity and community structure. It is worth noting that the water optical properties can be highly variable in some sites due to the influence of particular meteorologic and oceanographic conditions3,48, and that a single Kd may not represent the complex dynamics of the underwater light climate. Understanding this variability, although crucial for estimating the spatial–temporal variation in benthic primary production4,16, is beyond the scope of this study.

Zooxanthellate corals display contrasting photoacclimation responses, coral cover, colony morphologies and genetic richness along depth gradients33,49,50,51,52,53,54, which collectively suggest that coral species occupy different light niches. Particular colony geometries that optimize light capture and photosynthetic energy acquisition seem to be selected for maximizing the energy output at the colony level in specific light habitats49,54. Shallow-water corals, for example, adopt complex morphologies (e.g., branching and corymbose) to regulate within-colony light levels, potentially maximizing the light utilization efficiency and the photosynthetic output of the colony. In contrast, deep-water corals adopt flattened morphologies to maximize light capture and minimize self-shading50,51,52,53. These patterns suggest that a compromise between maximizing the photosynthetic production while minimizing the energy expenditure in photorepair may have been a driving force of colony morphology through coral evolution. Additionally, the association with symbionts with distinctive photoacclimatory potential can allow corals to cope with contrasting light climates along depth gradients and within colonies33,55,56. The location of maximum translocation of photosynthetic usable energy at intermediate irradiance may promote a more widely available specialization and evolutionary innovation space in that particular habitat with regard to variations in colony geometry and symbiotic associations, allowing more coral species to coexist8. Conversely, the variety of viable specializations that emerge from combinations between colony morphology and algal associations are predicted to be lower at both ends of the light intensity gradient due to reduced resource heterogeneity and stronger competition for resources7. Moreover, the reduced translocation of photosynthetic usable energy can lead to increased risk of species extinction and prevalence of few efficient competitors at exploiting the available light energy, either in excess or deficit.

Given their mixotrophic nature, heterotrophy is another aspect that has to be considered when analyzing the energy budget of zooxanthellate corals and its potential influence on biodiversity patterns across depths. Previous studies have demonstrated that some coral species are able to increase their metabolic reliance on heterotrophy to compensate for reduced photosynthetically derived energy acquisition in deep, low-light environments57,58, thereby potentially influencing biodiversity patterns across depths. Although we did not parameterize this aspect of the symbiosis in the documented productivity-biodiversity model, a potential effect of heterotrophic plasticity or trophic niche differentiation on coral richness cannot be discharged. However, the large component of the variability explained by the model and the contrasting features exhibited by zooxanthellate (non-facultative) corals suggest that coral adaptation and specialization to optimize solar energy utilization may have been a fundamental driving force in coral evolution, perhaps more important than enhancing heterotrophic feeding capabilities41,52. In this sense, if heterotrophy played a significant role in determining patterns of coral biodiversity with depth, a greater number of species with morphologies that facilitate suspension feeding (e.g., branching) would be expected to thrive with increasing depth. However, the empirical evidence indicates that flattened morphologies to maximize light capture prevail in the lower photic zone41,52.

Disturbance frequency and intensity have also been shown to affect the structure of coral communities, and it has been traditionally hypothesized that intermediate disturbance regimes lead to greater biodiversity17. The intermediate disturbance hypothesis, originally proposed as a conceptual model, has been supported and rejected in its capacity for explaining biodiversity patterns in ecological communities, both aquatic and terrestrial7,9. There are two major disturbances that can frequently affect shallow-water reefs, with little impact in deep-water counterparts: coral bleaching related to heat stress and high wave energy due to storms and hurricanes15,41. These disturbances can certainly influence both the local species diversity and the community composition of symbiotic corals18,19. However, the consistency of the explanatory power of the productivity-biodiversity model despite local environmental and ecological conditions among sites suggests that disturbances and other environmental factors may alter the location of the node as well as the slope of the increasing and/or decreasing phase of the unimodal curve, but not the overall pattern. Indeed, our results suggest that overlooking the role of productivity can obscure the underlying cause of biodiversity patterns in coral communities.

In summary, the results of this analysis indicate that solar energy and photosynthetic productivity are major driving forces of biodiversity patterns along depth gradients in symbiotic coral communities. The symbiosis with photosynthesizing dinoflagellates was a successful adaptive solution for corals to thrive in oligotrophic environments11 which, ultimately, led to the consolidation of one of the most biodiverse ecosystems on the planet. The gradient of downwelling irradiance mediated by the water optical properties, coupled with the metabolic and physiological constraints imposed by the obligate symbiosis with a primary producer, seems to be primary determinants for the establishment of global-scale patterns of biodiversity in scleractinian corals. As a major reef-building taxa, there patterns may potentially affect the associated communities that depend on corals for food and shelter. The increased worldwide degradation of the water optical properties in coastal environments associated with coastal development, nutrient enrichment, massive algal blooms, and terrestrial runoff59,60,61, may be an important underlying cause of biodiversity loss and change in the assemblage structure of coral reef communities. Local conservation actions seeking to maintain the water optical quality and the underwater light climate are essential to preserve coral biodiversity, while concerted global action to limit greenhouse emissions and slow global warming continues to move forward.

We assume that the energy cost of repairing the light-induced damage of the photosynthetic apparatus in the zooxanthellae is the main limiting factor for the photosynthetic usable energy supply by the algae to their coral hosts. This assumption is supported by evidence that the continual replacement of proteins required for the re-assembly of PSII reaction centers can be the largest single contributor to the costs of maintenance in primary producers45,46. The photosynthetic usable energy supply to corals is correlative to the "contribution of zooxanthellae carbon to the animal host respiration", a concept originally introduced by Muscatine et al.62. In their original formulation, however, the translocation of photosynthetically fixed carbon was calculated based on the oxygen production/consumption of the holobiont using a factor that was independent of light availability.

We determined the pattern of change of the zooxanthellae energy expenditure in photorepair as a function of light exposure using small fragments (~ 5 × 5 cm, n = 50) of the Caribbean coral Porites astreoides collected from La Bocana Reef at 5 m depth in Puerto Morelos, México. Two weeks after acclimation in a running seawater aquarium with neutral density filters simulating the light intensity at the collection depth and controlling temperature at 28 °C, corals were exposed to contrasting light intensities to measure changes in: (1) the photosystem II (PSII) effective (ΔF/Fm’) and maximum (Fv/Fm) photochemical quantum yields over diurnal cycles using pulse amplitude modulated technique23,33, and (2) the PSII half-time (t1/2) by analyzing the diurnal oscillations of ΔF/Fm’ of coral fragments exposed to 100 µg of chloramphenicol (CAP), a chloroplast protein-synthesis inhibitor30, relative to controls without CAP.

The PSII t1/2 was calculated as the amount of time required to inactivate 50% of the PSII. We explored the optimum relationship between PSII t1/2 and light exposure (as % of sea surface) using power-function modeling. The following power function explained the variation of the PSII t1/2 with light exposure:

where M is the coefficient of the power function representing the maximum theoretical value of PSII t1/2, E is irradiance (used here to denote variation in light exposure), and Δ is the rate of change of PSII t1/2 with respect to available light.

The energy expenditure in photorepair by the symbiotic algae over a diurnal cycle (Ca) was calculated through the relation:

where R corresponds to the relative energy cost of protein turnover for the reassembly of PSII reaction centers over a diurnal cycle (12 h of daylight).

To analyze the relationship between photosynthetic usable energy and coral biodiversity across depths, we used a numerical model with a lumped-parameter approach assuming that a specific set of parameters represent the average response of the whole community63. The model acknowledged three physical and physiological principles that govern the vertical distribution of light in the water column and the photosynthetic activity of primary producers: (1) solar radiation diminishes exponentially with depth according to the water optical properties (i.e., local Kd); (2) algal photosynthetic production describes a hyperbolic tangent response to light availability; and (3) energy expenditure for maintaining the algal photosynthetic activity is mediated by the light exposure following Eq. (2) and varies with depth according to Eq. (3). This approach allows us to draw general conclusions about light-driven processes at the community level at the expense of simplifying the complex dynamics of the water optical properties3,16,48, and acclimation strategies of corals to maximize performance under particular light and temperature regimes51,64,65.

The available solar radiation across depths, as the energy source for the symbiotic algae to fix inorganic carbon via photosynthesis, is calculated following Kirk13 as:

where Ez and E0 are the irradiance (E) at depth z and just below the sea surface, respectively.

The hyperbolic response of algal photosynthesis (Pg) to irradiance is calculated according to Jassby and Platt34, as:

where α is the photosynthetic efficiency corresponding to the slope of the linear increase of photosynthesis at subsaturating irradiance and Pmax is the maximum photosynthetic rates at light saturation. We assume no temporal reduction in light-saturated photosynthesis in response to photoinhibition at excessive irradiance considering the apparent absence of this phenomenon in symbiotic dinoflagellates in hospite24.

To estimate the variation of Ca as a function of depth, Ca(z), we combined Eqs. (1), (2) and (3) as:

Finally, the photosynthetic usable energy supply (PUES) is parameterized with the following equation:

We tested the explanatory power of the productivity-biodiversity model in a wide range of coral reef habitats along the world's major hotspots of marine biodiversity (Indo-Pacific, Caribbean and Indian Ocean66). We analyzed published datasets from studies that reported the change in number of symbiotic coral species with depth in seven reefs with contrasting geological history, environmental conditions and biodiversity patterns (unimodal, monotonic increase and monotonic decrease) (Supplementary Table S1)10,37,38,39,40. Data were extracted directly from the text, tables, or original datasets. When the data was only presented in graphical format, raw values were extracted using the software WebPlotDigitizer version 4.4 (https://automeris.io/WebPlotDigitizer). When the local Kd was not reported in the manuscript, it was calculated from the light change across depths following Eq. (3). When neither of these options was possible, it was extracted from the literature for the same or comparable sites (Supplementary Table S2).

For each site, a two-step optimization procedure was used to calculate those parameters that could not be estimated independently with empirical data (M, Δ, Pmax, α, R) (Supplementary Table S2). The optimization of these parameters produced the best possible fit between the model output and the mean number of observed species across depths. First, an iterative algorithm that allows the variation of each parameter within a suspected range while fulfilling the model equations exactly was used for fitting a linear model with coral richness as a function of the photosynthetic usable energy. Using the Akaike Information Criterion (AIC), we selected the set of values that minimized the misfit between the productivity-biodiversity model output and the empirical data. Then, the values of the targeted parameters that minimized the AIC were used as the starting values in a bound-constrained minimization algorithm (optim function) based on the Nelder-Mead method67. The lower and upper bounds of each parameter were determined based on empiric data and physiological limits and were constant among sites. A diurnal cycle of solar radiation describing a 12 h sinusoidal curve peaking at 1800 μmol quanta m−2 s−1 at noon was used as forcing function (chosen maximum is a random value close to naturally occurring maximum irradiance at ocean surface in tropical areas13). All analyses were conducted using R version 3.6.168.

Data generated and analyzed during the current study, as well as the code used to run these analyses, are openly available on Figshare (https://figshare.com/s/b73708c98e8108ec0352).

Field, C. B., Behrenfeld, M. J., Randerson, J. T. & Falkowski, P. Primary production of the biosphere: Integrating terrestrial and oceanic components. Science 281, 237–240. https://doi.org/10.1126/science.281.5374.237 (1998).

Article CAS PubMed Google Scholar

Valladares, F. In Progress in Botany Vol. 64 (eds Esser, K. et al.) 439–471 (Springer, 2003).

Chapter Google Scholar

Anthony, K. R. N., Ridd, P. V., Orpin, A. R., Larcombe, P. & Lough, J. Temporal variation of light availability in coastal benthic habitats: Effects of clouds, turbidity, and tides. Limnol. Oceanogr. 49, 2201–2211. https://doi.org/10.4319/lo.2004.49.6.2201 (2004).

Article Google Scholar

Gattuso, J. P. et al. Light availability in the coastal ocean: Impact on the distribution of benthic photosynthetic organisms and their contribution to primary production. Biogeosciences 3, 489–513. https://doi.org/10.5194/bg-3-489-2006 (2006).

Article Google Scholar

Wright, D. H. Species-energy theory: An extension of species-area theory. Oikos 41, 496–506 (1983).

Article Google Scholar

Cusens, J., Wright, S. D., McBride, P. D. & Gillman, L. N. What is the form of the productivity–animal-species-richness relationship? A critical review and meta-analysis. Ecology 93, 2241–2252. https://doi.org/10.1890/11-1861.1 (2012).

Article PubMed Google Scholar

Rosenzweig, M. L. & Abramsky, Z. in Species Diversity in Ecological Communities. Historical and Geographical Perspectives (eds Ricklefs, R. E. & Schluter, D.) Ch. 5, 52–65 (The University of Chicago Press, 1993).

Abrams, P. A. Monotonic or unimodal diversity-productivity gradients: What does competition theory predict?. Ecology 76, 2019–2027 (1995).

Article Google Scholar

Huston, M. A. Disturbance, productivity, and species diversity: Empiricism vs. logic in ecological theory. Ecology 95, 2382–2396 (2014).

Article Google Scholar

Roberts, T. E. et al. Testing biodiversity theory using species richness of reef-building corals across a depth gradient. Biol. Lett. 15, 20190493. https://doi.org/10.1098/rsbl.2019.0493 (2019).

Article PubMed PubMed Central Google Scholar

Frankowiak, K. et al. Photosymbiosis and the expansion of shallow-water corals. Sci. Adv. 2, e1601122. https://doi.org/10.1126/sciadv.1601122 (2016).

Article CAS PubMed PubMed Central Google Scholar

Goreau, T. F. & Goreau, N. I. The physiology of skeleton formation in corals. II. Calcium deposition by hermatypic corals under various conditions in the reef. Biol. Bull. 117, 239–250. https://doi.org/10.2307/1538903 (1959).

Article CAS Google Scholar

Kirk, J. T. O. Light and Photosynthesis in Aquatic Ecosystems 3rd edn. (Cambridge University Press, 2011).

Google Scholar

Stoddart, D. R. Ecology and morphology of recent coral reefs. Biol. Rev. 44, 433–498. https://doi.org/10.1111/j.1469-185X.1969.tb00609.x (1969).

Article Google Scholar

Lesser, M. P., Slattery, M. & Leichter, J. J. Ecology of mesophotic coral reefs. J. Exp. Mar. Biol. Ecol. 375, 1–8 (2009).

Article Google Scholar

Ackleson, S. G. Light in shallow waters: A brief research review. Limnol. Oceanogr. 48, 323–328. https://doi.org/10.4319/lo.2003.48.1_part_2.0323 (2003).

Article Google Scholar

Connell, J. H. Diversity in tropical rain forests and coral reefs. High diversity of trees and corals is maintained only in a nonequilibrium state. Science 199, 1302–1310. https://doi.org/10.1126/science.199.4335.1302 (1978).

Article CAS PubMed Google Scholar

Dollar, S. J. Wave stress and coral community structure in Hawaii. Coral Reefs 1, 71–81. https://doi.org/10.1007/BF00301688 (1982).

Article Google Scholar

Hughes, T. P. Community structure and diversity of coral reefs: The role of history. Ecology 70, 275–279. https://doi.org/10.2307/1938434 (1989).

Article Google Scholar

Fraser, R. H. & Currie, D. J. The species richness-energy hypothesis in a system where historical factors are thought to prevail: Coral reefs. Am. Nat. 148, 138–159 (1996).

Article Google Scholar

Cornell, H. V. & Karlson, R. H. Coral species richness: Ecological versus biogeographical influences. Coral Reefs 19, 37–49 (2000).

Article Google Scholar

Bellwood, D. R., Hughes, T., Connolly, S. & Tanner, J. Environmental and geometric constraints on Indo-Pacific coral reef biodiversity. Ecol. Lett. 8, 643–651. https://doi.org/10.1111/j.1461-0248.2005.00763.x (2005).

Article Google Scholar

Brown, B. E. et al. Diurnal changes in photochemical efficiency and xanthophyll concentrations in shallow water reef corals: Evidence for photoinhibition and photoprotection. Coral Reefs 18, 99–105 (1999).

Article Google Scholar

Hoegh-Guldberg, O. & Jones, R. J. Photoinhibition and photoprotection in symbiotic dinoflagellates from reef-building corals. Mar. Ecol. Prog. Ser. 183, 73–86. https://doi.org/10.3354/meps183073 (1999).

Article Google Scholar

Lesser, M. P. & Gorbunov, M. Y. Diurnal and bathymetric changes in chlorophyll fluorescence yields of reef corals measured in situ with a fast repetition rate fluorometer. Mar. Ecol. Prog. Ser. 212, 69–77. https://doi.org/10.3354/meps212069 (2001).

Article CAS Google Scholar

Hoogenboom, M. O., Anthony, K. R. N. & Connolly, S. R. Energetic cost of photoinhibition in corals. Mar. Ecol. Prog. Ser. 313, 1–12. https://doi.org/10.3354/meps313001 (2006).

Article CAS Google Scholar

Huot, Y. & Babin, M. Chlorophyll a Fluorescence in Aquatic Sciences: Methods and Applications 31–74 (Springer, 2010).

Book Google Scholar

Warner, M. E., Lesser, M. P. & Ralph, P. J. Chlorophyll a Fluorescence in Aquatic Sciences: Methods and Applications Ch. Chapter 10, 209–222 (Springer Science+Business Media B.V., 2010).

Skirving, W. et al. Remote sensing of coral bleaching using temperature and light: Progress towards an operational algorithm. Remote Sens. 10, 18 (2018).

Article Google Scholar

Enríquez, S., Merino, M. & Iglesias-Prieto, R. Variations in the photosynthetic performance along the leaves of the tropical seagrass Thalassia testudinum. Mar. Biol. 140, 891–900. https://doi.org/10.1007/s00227-001-0760-y (2002).

Article CAS Google Scholar

Sundby, C., McCaffery, S. & Anderson, J. M. Turnover of the photosystem II D1 protein in higher plants under photoinhibitory and nonphotoinhibitory irradiance. J. Biol. Chem. 268, 25476–25482 (1993).

Article CAS PubMed Google Scholar

Tyystjärvi, E. & Aro, E. M. The rate constant of photoinhibition, measured in lincomycin-treated leaves, is directly proportional to light intensity. Proc. Natl. Acad. Sci. U. S. A. 93, 2213–2218. https://doi.org/10.1073/pnas.93.5.2213 (1996).

Article PubMed PubMed Central Google Scholar

Iglesias-Prieto, R., Beltrán, V. H., LaJeunesse, T. C., Reyes-Bonilla, H. & Thomé, P. E. Different algal symbionts explain the vertical distribution of dominant reef corals in the eastern Pacific. Proc. R. Soc. Lond. B 271, 1757–1763. https://doi.org/10.1098/rspb.2004.2757 (2004).

Article CAS Google Scholar

Jassby, A. D. & Platt, T. Mathematical formulation of the relationship between photosynthesis and light for phytoplankton. Limnol. Oceanogr. 21, 540–547 (1976).

Article CAS Google Scholar

Long, S. P., Humphries, S. & Falkowski, P. G. Photoinhibition of photosynthesis in nature. Annu. Rev. Plant Physiol. Plant Mol. Biol. 45, 633–662. https://doi.org/10.1146/annurev.pp.45.060194.003221 (1994).

Article CAS Google Scholar

Huner, N. P. A., Öuist, G. & Sarhan, F. Energy balance and acclimation to light and cold. Trends Plant Sci. 3, 224–230 (1998).

Article Google Scholar

Sheppard, C. R. C. Coral cover, zonation and diversity on reef slopes of Chagos Atolls, and population structures of the major species. Mar. Ecol. Prog. Ser. 2, 193–205 (1980).

Article Google Scholar

Huston, M. A. Patterns of species diversity in relation to depth at Discovery Bay, Jamaica. Bull. Mar. Sci. 37, 928–935 (1985).

Google Scholar

Loya, Y. Community structure and species diversity of hermatypic corals at Eilat, Red Sea. Mar. Biol. 13, 100–123. https://doi.org/10.1007/BF00366561 (1972).

Article Google Scholar

Chow, G. S. E., Chan, Y. K. S., Jain, S. S. & Huang, D. Light limitation selects for depth generalists in urbanised reef coral communities. Mar. Environ. Res. 147, 101–112. https://doi.org/10.1016/j.marenvres.2019.04.010 (2019).

Article CAS PubMed Google Scholar

Kahng, S. E. et al. Community ecology of mesophotic coral reef ecosystems. Coral Reefs 29, 255–275. https://doi.org/10.1007/s00338-010-0593-6 (2010).

Article Google Scholar

Iglesias-Prieto, R. Temperature-dependent inactivation of Photosystem II in symbiotic dinoflagellates. in Proc. 8th Int. Coral Reef Sym, 1313–1318 (1997).

Jones, R. J., Hoegh-Guldberg, O., Larkum, A. W. D. & Schreiber, U. Temperature-induced bleaching of corals begins with impairment of the CO2 fixation mechanism in zooxanthellae. Plant Cell Environ. 21, 1219–1230. https://doi.org/10.1046/j.1365-3040.1998.00345.x (1998).

Article CAS Google Scholar

Hennige, S. J., Suggett, D. J., Warner, M. E., McDougall, K. E. & Smith, D. J. Photobiology of Symbiodinium revisited: Bio-physical and bio-optical signatures. Coral Reefs 28, 179–195. https://doi.org/10.1007/s00338-008-0444-x (2008).

Article Google Scholar

Quigg, A. & Beardall, J. Protein turnover in relation to maintenance metabolism at low photon flux in two marine microalgae. Plant Cell Environ. 26, 693–703. https://doi.org/10.1046/j.1365-3040.2003.01004.x (2003).

Article CAS Google Scholar

Järvi, S., Suorsa, M. & Aro, E. M. Photosystem II repair in plant chloroplasts—Regulation, assisting proteins and shared components with photosystem II biogenesis. Biochim. Biophys. Acta Bioenerg. 900–909, 2015. https://doi.org/10.1016/j.bbabio.2015.01.006 (1847).

Article CAS Google Scholar

Jokiel, P. L. Solar ultraviolet radiation and coral reef epifauna. Science 207, 1069–1071 (1980).

Article CAS PubMed Google Scholar

López-Londoño, T. et al. Physiological and ecological consequences of the water optical properties degradation on reef corals. Coral Reefs 40, 1243–1256. https://doi.org/10.1007/s00338-021-02133-7 (2021).

Article Google Scholar

Vermeij, M. J. A. & Bak, R. P. M. How are coral populations structured by light? Marine light regimes and the distribution of Madracis. Mar. Ecol. Prog. Ser. 233, 105–116. https://doi.org/10.3354/meps233105 (2002).

Article Google Scholar

Hoogenboom, M. O., Connolly, S. R. & Anthony, K. R. N. Interactions between morphological and physiological plasticity optimize energy acquisition in corals. Ecology 89, 1144–1154. https://doi.org/10.1890/07-1272.1 (2008).

Article PubMed Google Scholar

Kaniewska, P., Anthony, K., Sampayo, E., Campbell, P. & Hoegh-Guldberg, O. Implications of geometric plasticity for maximizing photosynthesis in branching corals. Mar. Biol. 161, 313–328 (2014).

Article CAS Google Scholar

Kramer, N., Tamir, R., Eyal, G. & Loya, Y. Coral morphology portrays the spatial distribution and population size-structure along a 5–100 m depth gradient. Front. Mar. Sci. https://doi.org/10.3389/fmars.2020.00615 (2020).

Article Google Scholar

Lesser, M. P., Mobley, C. D., Hedley, J. D. & Slattery, M. Incident light on mesophotic corals is constrained by reef topography and colony morphology. Mar. Ecol. Prog. Ser. 670, 49–60. https://doi.org/10.3354/meps13756 (2021).

Article Google Scholar

Prada, C. et al. Linking photoacclimation responses and microbiome shifts between depth-segregated sibling species of reef corals. R. Soc. Open Sci. 9, 211591. https://doi.org/10.1098/rsos.211591 (2022).

Article PubMed PubMed Central Google Scholar

Rowan, R., Knowlton, N., Baker, A. & Jara, J. Landscape ecology of algal symbionts creates variation in episodes of coral bleaching. Nature 388, 265–269. https://doi.org/10.1038/40843 (1997).

Article CAS PubMed Google Scholar

Warner, M. E., LaJeunesse, T. C., Robison, J. D. & Thur, R. M. The ecological distribution and comparative photobiology of symbiotic dinoflagellates from reef corals in Belize: Potential implications for coral bleaching. Limnol. Oceanogr. 51, 1887–1897. https://doi.org/10.4319/lo.2006.51.4.1887 (2006).

Article Google Scholar

Anthony, K. R. N. & Fabricius, K. E. Shifting roles of heterotrophy and autotrophy in coral energetics under varying turbidity. J. Exp. Mar. Biol. Ecol. 252, 221–253 (2000).

Article CAS PubMed Google Scholar

Hoogenboom, M., Rodolfo-Metalpa, R. & Ferrier-Pagès, C. Co-variation between autotrophy and heterotrophy in the Mediterranean coral Cladocora caespitosa. J. Exp. Biol. 213, 2399–2409 (2010).

Article PubMed Google Scholar

Carlson, R. R., Foo, S. A. & Asner, G. P. Land use impacts on coral reef health: A ridge-to-reef perspective. Front. Mar. Sci 6, 562. https://doi.org/10.3389/fmars.2019.00562 (2019).

Article Google Scholar

Wang, M. et al. The great Atlantic Sargassum belt. Science 365, 83–87. https://doi.org/10.1126/science.aaw7912 (2019).

Article CAS PubMed Google Scholar

Alvarez-Filip, L., González-Barrios, F. J., Pérez-Cervantes, E., Molina-Hernández, A. & Estrada-Saldívar, N. Stony coral tissue loss disease decimated Caribbean coral populations and reshaped reef functionality. Commun. Biol. 5, 440. https://doi.org/10.1038/s42003-022-03398-6 (2022).

Article PubMed PubMed Central Google Scholar

Muscatine, L., McCloskey, L. R. & Marian, R. E. Estimating the daily contribution of carbon from zooxanthellae to coral animal respiration. Limnol. Oceanogr. 26, 601–611. https://doi.org/10.4319/lo.1981.26.4.0601 (1981).

Article CAS Google Scholar

Jørgensen, S. E. & Bendoricchio, G. Fundamentals of Ecological Modelling 3rd edn, Vol. 21 (Elsevier Sceince B. V., 2001).

Google Scholar

Hennige, S. J. et al. Acclimation and adaptation of scleractinian coral communities along environmental gradients within an Indonesian reef system. J. Exp. Mar. Biol. Ecol. 391, 143–152. https://doi.org/10.1016/j.jembe.2010.06.019 (2010).

Article Google Scholar

Scheufen, T., Iglesias-Prieto, R. & Enríquez, S. Changes in the number of symbionts and Symbiodinium cell pigmentation modulate differentially coral light absorption and photosynthetic performance. Front. Mar. Sci 4, 309. https://doi.org/10.3389/fmars.2017.00309 (2017).

Article Google Scholar

Veron, J. E. N. Corals in Space and Time. The Biogeography and Evolution of the Scleractinia 321 (Cornell University Press, 1995).

Google Scholar

Nelder, J. A. & Mead, R. A simplex method for function minimization. J. Comput. 7, 308–313. https://doi.org/10.1093/comjnl/7.4.308 (1965).

Article MathSciNet MATH Google Scholar

R: A languate and environment for statistical computing. Retrieved from http://www.R-project.org (R Foundation for Statistical Computing, Vienna, Austria, 2010).

Download references

We dedicate this manuscript to the memory of Professor Robert (Bob) K. Trench, a seminal figure for our current understanding of algal invertebrate symbioses. The authors thank Todd LaJeunesse, Mitzy F. Porras, Carlos Prada, Teresa Morán-López and M. Isabel Martínez-Rugerio for their valuable contributions; the National Commission on Aquaculture and Fisheries (CONAPESCA) for research permit (DGOPA 08606.251011.3021). A PASPA fellowship from the DGAPA-UNAM supported the visit of SE to the Biology Department at Pennsylvania State University (PSU). This work was supported by Pennsylvania State University startup funds to RI-P.

Department of Biology, Pennsylvania State University, University Park, PA, 16802, USA

Tomás López-Londoño, Kelly Gómez-Campo & Roberto Iglesias-Prieto

Unidad Académica de Sistemas Arrecifales, Instituto de Ciencias del Mar y Limnología, Universidad Nacional Autónoma de México, 77500, Cancún, Quintana Roo, México

Xavier Hernández-Pech & Susana Enríquez

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

T.L.-L., K.G.-C., S.E. and R.I.-P. designed the study. X.H.-P. and K.G.-C. conducted coral experiments. T.L.-L. conducted meta-analysis and eco-physiological modelling under the supervision of S.E. and R.I.-P. T.L.-L. and R.I.-P. wrote the manuscript with contributions and feedback from all authors. All authors read and approved the final version of the manuscript.

Correspondence to Tomás López-Londoño or Roberto Iglesias-Prieto.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

López-Londoño, T., Gómez-Campo, K., Hernández-Pech, X. et al. Photosynthetic usable energy explains vertical patterns of biodiversity in zooxanthellate corals. Sci Rep 12, 20821 (2022). https://doi.org/10.1038/s41598-022-25094-5

Download citation

Received: 10 June 2022

Accepted: 24 November 2022

Published: 02 December 2022

DOI: https://doi.org/10.1038/s41598-022-25094-5

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.