Co-benefits not trade-offs associated with heat tolerance in a reef building coral

doi:10.21203/rs.3.rs-1973987/v1

Download PDF

Article

Co-benefits not trade-offs associated with heat tolerance in a reef building coral

https://doi.org/10.21203/rs.3.rs-1973987/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 12 Apr, 2023

Read the published version in Communications Biology →

Version 1

posted

You are reading this latest preprint version

As marine species adapt to climate change, their heat tolerance will likely be under strong selection. Yet trade-offs between heat tolerance and other life history traits could compromise natural adaptation or restorative assisted evolution. This is particularly important for ecosystem engineers, such as reef-building corals, which support biodiversity yet are vulnerable to heatwave-induced mass bleaching and mortality. Here, we exposed 70 colonies of the reef-building coral Acropora digitifera to a long-term marine heatwave emulation experiment. We tested for trade-offs between heat tolerance and three traits – colony growth (3D photogrammetry models), fecundity (oocyte counts), and symbiont community composition (ITS2 sequencing). Despite observing remarkable within-population variability in heat tolerance, all coral colonies were dominated by Cladocopium C40 symbionts. We found no evidence for trade-offs between heat tolerance and fecundity or growth. Contrary to expectations, positive associations emerged with growth, such that faster-growing colonies tended to bleach and die at higher levels of heat stress. Collectively, our results suggest that these corals exist on an energetic continuum where some individuals have co-benefits across multiple traits. Within populations, trade-offs between heat tolerance and either growth or fecundity may not be major barriers to natural adaptation or the success of restorative assisted evolution interventions.

Biological sciences/Ecology/Climate-change ecology

Biological sciences/Ecology/Evolutionary ecology

Biological sciences/Ecology/Population dynamics

Bayesian statistics

climate change

coral bleaching and mortality

heat tolerance

ITS2 sequencing

life history traits

R-INLA

structure-from-motion.

Ocean warming is causing profound changes in marine ecosystems, and to keep pace and avoid extirpation, species must either migrate or adapt. The adaptive capacity of ecosystem engineers, such as reef-building corals, will play a disproportionately large role in the future biodiversity and function of marine ecosystems. Coral reefs continue to face unprecedented declines due to mass coral bleaching and mortality events caused by marine heatwaves ^1–3. These extended periods of anomalously high ocean temperatures are increasing in frequency and intensity under climate change ^4,5. The ability of individual corals to survive levels of heat stress sufficient to induce mass bleaching and mortality, hereafter "heat tolerance", will likely emerge as an important trait under natural selection during the coming decades. Survivorship rather than just bleaching is core component of heat tolerance as bleached corals can still recover and persist ^6,7.

Considerable variability in coral heat tolerance exists between individuals, even within a single coral population on a single reef ⁸. Currently, there are growing efforts to test novel restorative interventions, such as assisted evolution, which aims to bolster the heat tolerance of coral populations by seeding reefs with more tolerant coral colonies ^9,10. Understanding variation in heat tolerance is crucial to estimating the capacity for natural adaptation to climate change and the efficacy of assisted evolution interventions. Under both natural selection and assisted evolution, the relationships between coral heat tolerance and fitness traits (i.e., reproduction and survival – even in non-heat stress years) or other ecological traits (e.g., growth and fecundity which may or may not affect fitness) are of fundamental importance to future population persistence.

Organisms are limited by resource availability, forcing them to balance resource allocation between different physiological processes leading to “trade-offs” between resource-intensive traits. For instance, successful strategies to deal with drought stress are well known in long-lived birds, where less energy is allocated to reproduction in drought years to preserve cell maintenance and growth ¹¹. Such trade-offs will always occur between resource-limited processes. However, sometimes apparent positive associations can be found between resource-intensive traits, even when a trade-off might be expected. The variability of total resource budgets among individuals can explain this phenomenon ¹², and can be associated with co-benefits across multiple traits and co-tolerance to the impacts of multiple stressors among individuals ¹³, despite the presence of more nuanced resource trade-offs within individuals. For instance, in oysters there is genetic evidence for trade-offs (negative correlations) between reproductive effort and both survival and growth ¹⁴, likely due to resource allocation. However, these negative correlations turn positive when these traits are measured across numerous oysters under feeding treatments ¹⁴. This can occur due to variability in resource acquisition among individuals, which in turn can lead to differences in their total resource budgets. Subtle trait trade-offs are then easily masked by the broader population-scale energetic continuum (i.e., gradient of total resource budgets among individuals).

Heat tolerance in reef-building corals has been shown to have negative associations with growth at the colony level ^15,16, suggesting a resource trade-off. This can lead to considerable negative impacts on coral reefs at the ecosystem level ¹⁷. For instance, Acropora spp. and Pocillopora damicornis corals dominated by thermally tolerant Durusdinium spp. symbiotic microalgae show considerable reductions in vital cell processes, such as carbon storage ¹⁵, photosynthetic efficiency, and energetics ^15,16. Ultimately this results in reduced coral growth rates ^18,19. Notably, this growth disadvantage can be eliminated by extreme warming, as heat tolerance becomes more important than growth under such conditions ²⁰. The presence of mixed symbiont communities and symbiont shuffling post-bleaching can lead to flexibility in the magnitude of heat tolerance-growth trade-offs ^18,21. However, for other coral genera (e.g., Montipora sp.), even the expectedly heat-tolerant symbiont Durusdinium glynnii (rather than Cladocopium spp.) has been found to have no influence on coral host physiology and metabolic profiles ²².

Many coral populations are dominated by a single symbiont taxon or a single symbiont community type ^23,24. In these cases, do trade-offs between heat tolerance and other traits persist? Recent genomic evidence based on corals from contrasting thermal environments suggests that the shift in allele frequencies associated with coral host-derived heat tolerance must be associated with a fitness cost ²³. However, it is yet to be tested whether trade-offs between heat tolerance and other ecological traits exist for corals that share the same Symbiodiniaceae community. Considering their prevalence in numerous other taxa, including Crustacea, Insecta, and Chordata ^25–27, it is likely that heat tolerance-related trait trade-offs also affect coral hosts. As coral adaptation occurs locally, not globally, and since endosymbiont communities are relatively uniform across local scales, it is important to resolve the extent of host-derived heat tolerance trade-offs to better predict coral adaptation to climate change.

Egg and sperm development are resource intensive processes and have been suggested as potential costs to growth and heat tolerance ²⁸. It is reasonable to expect trade-offs between heat tolerance and fecundity given the evidence for heat tolerance-growth trade-offs in corals, and the fact that growth and gamete production are both resource-intensive processes. Evidence has shown that temperature stress can reduce hard and soft coral fecundity (i.e., egg density and volume) ^29,30, and has suggested that corals may be able to reabsorb their oocytes to divert energy away from reproduction and into growth instead under certain types of stress such as fragmentation ^31,32. However, there has not yet been an assessment of the relationship between fecundity and heat tolerance to test for associations or trade-offs which could be crucial to understand population fitness and performance.

Here, we tested for ecologically relevant associations and trade-offs between heat tolerance and three ecological traits in a common species of Indo-Pacific coral: colony growth, fecundity, and symbiont community composition. To measure these traits we combined: (i) a long-term (sensu ³³) 5-week marine heatwave emulation experiment to measure heat tolerance; (ii) interannual comparisons of 3D models to measure growth; (iii) polyp counts and dissections to measure fecundity; and (iv) ITS2 sequencing to determine Symbiodiniaceae community composition.

2.1 | Heat tolerance variability

Within the studied coral population, we found substantial variability in coral heat tolerance measured as bleaching and mortality responses throughout a marine heatwave emulation experiment. At the beginning of the experiment all colonies had all replicate fragments healthy, corresponding to BSI (bleaching and survival index) vales of 1. The final heat stress exposure reached a DHW (degree heating weeks) of 10.7 ºC-weeks (Fig. 1A-B), a level which would likely induce a mass bleaching and mortality event in nature. By this final exposure, 47% of colonies had all replicate fragments dead (BSI = 0), while most remaining colony fragments were bleached (25% of all fragments), translating to BSI values < 0.56 (Fig. 1C). Meanwhile, in the unstressed procedural control tanks, all representative fragments from each colony remained alive, showing that the experimental setup (aquarium lights, flow through filtered water etc.) was suitable for keeping fragments alive and that heat stress was the ultimate cause of bleaching and mortality responses (i.e., rather than effects due to the tank setup). Coral heat tolerance throughout the experimental exposure, expressed as the average BSI, ranged from 0.47 to 0.91 among colonies. The critical level of accumulated heat stress (DHW – degree heating weeks) at onset of bleaching and mortality (BSI < 0.75) was highly variable across colonies, with an average critical DHW of 6.7 ºC-weeks (± 1.1 ºC-weeks SD) but ranging from 4.3 to 9.2°C-weeks in the least and most tolerant individuals. Notably, that corresponds to almost an entire categorical shift of the NOAA Coral Reef Watch bleaching alert system (e.g., Alert Level 1 to 2, moving from mass bleaching expected to mass mortality expected).

2.2 | Similar Symbiodiniaceae communities

Symbiodiniaceae community composition was consistent among colonies, whereby all individual colonies were dominated by C40 Cladocopium spp. (Fig. 1D) and 96% of colonies contained a single ITS2 type profile. For 70% of colonies, the dominant symbiont strain was C40-C3-C115-C40h. Only 2 colonies had mixed symbiont communities, with 29% relative abundance of C15h-C15hf-C15hg in colony A75 and 31% D1/D4-D4c-D1c-D1h-D2 in colony A96. There was no significant effect of symbiont ITS2 type on bleaching and mortality responses (Fig. 1D) with strongly overlapping confidence intervals across all ITS2 type profile groups (Fig. S1) (GLMM Tukey test, P > 0.05 for all pairwise comparisons, Table S2). This was also shown from principal component analysis based on ITS2 type profiles and grouping colony heat tolerance into broad equal-sized categories (average BSI: high ≥ 0.8 > medium ≥ 0.7 > low; Fig. S2), with strong overlap among each category.

2.3 | 3D colony growth

Coral colony growth was highly variable in terms of both live surface area growth (216 ± 722 cm²/yr, average ± SD) and volumetric growth (404 ± 647 cm³/yr, average ± SD). There were a range of relatively fast-growing colonies, stable-sized colonies, and others that either lost live tissue or volume through partial mortality or breakage (example colonies of each growth type shown in Fig. 1E, and all colony 3D models shown in Fig. S4). Once growth metrics were corrected by initial colony size, there was no significant trend between live surface area growth and colony size (linear regression, P > 0.05; Fig. S6). However, there was a significant negative trend between total volumetric growth and colony size, whereby larger corals had lower rates of volumetric change (linear regression, P = 0.042; Fig. S6). Notably, colony partial mortality had no effect on fecundity or symbiont community traits (ANOVA and Tukey tests or pairwise Bonferroni-corrected Wilcoxon tests, all P values > 0.05; Fig. S7), a weak negative effect on heat tolerance (average BSI; ANOVA, P = 0.023; Fig. S7), and a stronger negative effect on growth which was expected as tissue loss reduces size (pairwise Bonferroni-corrected Wilcoxon tests for live surface area growth, P = 0.006; ANOVA and Tukey tests for volumetric growth, P = 0.005; Fig. S7).

2.4 | Fecundity

Egg density averaged 4.9 eggs per polyp (± 1.4 eggs per polyp SD) across all coral colonies, with polyps containing between 1 and 11 eggs. The geometric mean diameter of eggs was in line with that shown for Acropora digitifera ³⁴, although slightly larger (Fig. 1.F, red dashed line). Average egg volume across all colonies was 0.11 ± 0.03 mm³ SD, ranging from 0.06 to 0.2 mm³, similar to that shown in previous work ²⁹. The estimates of total colony egg production (274,590 ± 236,624 eggs/colony, average ± SD) and total colony egg volume (32 ± 27 cm³/colony, average ± SD) were highly variable among colonies, likely due to variability in colony size.

2.5 | Co-benefits but not trade-offs with heat tolerance

There was no evidence for trade-offs between heat tolerance (average BSI per colony) and either growth (Fig. 2A), fecundity (Fig. 2B), or Symbiodiniaceae community composition (Fig. S1, Fig. S2). The lack of a relationship with Symbiodiniaceae relates to the high similarity of the symbiont community among colonies. Contrary to expectations, we found positive associations between heat tolerance and growth metrics (live surface area growth and volumetric growth). Despite having a weak association due to high levels of uncertainty (i.e., 95% credible intervals of the slope posterior distributions intersected zero, Fig. 2A), the probability that these slopes were positive, and that growth and heat tolerance act in concert was high; 90% probability for live surface area growth and 94% probability for volumetric growth (Fig. S8). In comparison, any effect of fecundity metrics (total colony egg production and total colony egg volume) on heat tolerance were unmeasurable given our sampling design (Fig. 2B), with close to 50:50 odds of the relationship being negative or positive (trade-off or co-benefit) (Fig. S8).

Throughout the course of the 5-week heatwave emulation experiment (Fig. 3), the relationships between coral colony growth (measured from 3D models on the reef) and instantaneous bleaching mortality responses were variable. The onset of bleaching responses occurred after approximately 3 weeks of elevated temperatures (Fig. 1B) at a DHW value of approximately 4°C-weeks (Fig. 3A). However, bleaching onset was delayed in fragments of positive-growth colonies compared to negative-growth colonies by approximately 1°C-week (Fig. 3A). Accordingly, during this period of differential bleaching onset (at approximately 5°C-weeks), the slope of the BSI-growth relationship (concept shown in Fig. 3B) was strongly positive showing a pronounced peak (e.g., 95% credible interval of the slope not intersecting zero, Fig. 3C). Faster growing colonies needed higher levels of heat stress to trigger the onset of bleaching responses.

The BSI-growth slope progressed throughout the experiment with a double peak pattern. Moving on to 7°C-weeks, as bleaching responses aligned across all colonies of different growth rates, the BSI-growth slope decreased from the first bleaching-associated peak toward zero (Fig. 3C). As the heat stress exposure reached 9°C-weeks, the BSI-growth slope again increased to a second peak, this time reflecting the delayed onset of mortality responses in positive-growth colonies compared to negative-growth colonies. Toward the end of the exposure (10.7°C-weeks), the BSI-growth slope moved back toward zero, reflecting critical bleaching and mortality responses across all colonies regardless of growth rates.

Notably, beyond 7°C-weeks, the uncertainty around the BSI-growth slope widened (Fig. 3C), intersecting with zero. Yet despite this, the probability of a positive BSI-growth slope (i.e., no trade-off) remained at 72–91% throughout the later stages of the heat stress exposure as bleaching and mortality responses progressed to critical levels.

Throughout the heat stress exposure, the progression of instantaneous BSI-growth relationships for volumetric growth (Fig. S9) showed a similar pattern to that of live surface area growth (Fig. 3). Slopes were consistently positive with credible intervals deviating from zero at bleaching onset (Fig. S10). However, the progression of instantaneous BSI-fecundity relationships (total colony egg production, and total colony egg volume) was static, centred at zero throughout the heat stress exposure (Fig. S9). These results show that growth has a positive association with heat tolerance, but that fecundity has no measurable relationship with heat tolerance (Fig. S10).

In coral reef ecology, theory and evidence suggest that high coral heat tolerance is associated with a trade-off with growth, especially for corals hosting certain symbiotic dinoflagellates (e.g., Durusdinium trenchii). Here, we investigated heat tolerance-associated trait trade-offs in a shallow outer reef crest coral population primarily hosting C40-dominated Symbiodiniaceae communities. We found no evidence for trade-offs between coral heat tolerance and either growth or fecundity.

Contrary to expectations, we found positive associations between heat tolerance and growth. Indeed, during a marine heatwave emulation experiment, fragments taken from faster growing coral colonies on the reef were able to withstand higher levels of experimental heat stress before the onset of bleaching and mortality. Previous work has identified considerable trade-offs between heat tolerance and growth caused by the presence of different symbionts ¹⁹, flagging this as a potential barrier to successful coral adaptation under climate change ¹⁷. However, our results show that heat tolerance and growth can be positively associated, offering a more optimistic outlook for coral populations. Particularly we found a double-peak pattern in BSI-growth regression slopes throughout a 5-week heat stress exposure (Fig. 3C). First the onset of bleaching was delayed in fragments of faster growing coral colonies (Fig. 3C, first peak), and then the onset of mortality was delayed (Fig. 3C, second peak). Our findings suggest that selecting corals for heat tolerance either through natural selection (i.e., selective mortality of heat sensitive individuals) or assisted evolution (i.e., propagating heat tolerant individuals) may not compromise growth or fecundity.

A key consideration here is separating the distribution of traits within a population from what happens to corals as temperatures continue to rise under climate change. In general, we found that more heat tolerant individuals also tended to have higher growth rates. This implies that post-bleaching coral populations may not necessarily have lower overall growth. However, the existence of thermal optima – as demonstrated for coral calcification ³⁵ and photosynthesis ³⁶ – still imply that long-lived corals may experience declines in their growth as temperatures rise, even if they are the more heat tolerant members of the earlier population. As we did not measure growth post-bleaching, there is a need for future research to understand the plasticity of trait associations after stress. As such, further work is also needed to understand whether selection for heat tolerance can also select for other beneficial traits. Our study focussed on corals from similar depths on a single reef to limit the influence of environment on organism physiology. Further research is needed to understand if co-benefits between heat tolerance and growth are also present in other coral species and over broader spatial scales.

Under climate change, coral heat tolerance will likely be one of the most important fitness-related traits in determining population persistence or collapse ^37,38. However, heat tolerance can come at a cost to other traits, like growth. This premise is typically based on a Symbiodiniaceae-derived coral heat tolerance, where certain dominant symbiont taxa (e.g., Durusdinium trenchii) confer higher tolerance at the expense of photosynthetic energetics and ultimately growth ³⁹. Many coral populations, particularly in the Indo-Pacific, host Symbiodiniaceae communities that are either dominated by a single taxon of symbiotic algae or by a single community type with similar relative abundance of different symbiont taxa ^23,40. Determining whether heat tolerance trade-offs persist for corals hosting broadly similar symbiont communities can improve our understanding of coral population functioning and the potential of adaptation to climate change.

Trade-off theory suggests that corals have contrasting strategies, either being resistant to high temperatures or showing enhanced growth rates, with concurrent disadvantages being reported for both strategies ¹⁹. However, our results show that heat tolerance can be positively associated with colony growth, where some individuals would expectedly have higher fitness due to excelling in multiple ecological traits simultaneously. Organisms must partition resources among costly physiological processes, suggesting that one should find trade-offs between these processes. However, total resource budgets can be highly variable among individuals, especially in wild populations, due to processes such as the efficiency of resource acquisition from the environment ¹². As such, even though trade-offs must occur at some levels, they can be masked at ecological scales due to high variability of resource acquisition among individuals. In line with trade-off theory, such energetic variability among individuals can result in apparent positive associations among resource intensive traits ^12,14. It is likely that the physiological processes underpinning high growth rates are also linked to high resilience to heat stress, as suggested by the concept of co-tolerance ¹³. None of the colonies with high live surface area growth rates underwent partial mortality. Those with even low levels of partial mortality were likely to have been affected by other stressors, such as disease, competition, physical damage, or predation, and were also then associated with lower levels of heat tolerance. The physiological cost of tissue repair, fighting infection or regrowth after breakage could deplete energy reserves ^41,42 rendering corals more susceptible to bleaching and mortality under acute heat stress.

Tolerance to extreme temperature stress is a vital trait for corals in the weeks or months during marine heatwaves. However, heatwaves currently do not happen every year and generally occur only in the warmest months. As such, heat tolerance is unlikely to directly benefit corals during cooler months or years, without considering associations between heat tolerance and other traits. Comparatively, other traits like growth or fecundity are of importance throughout the year in sexually mature adult corals (i.e., over 3 years old for Acropora spp. ⁴³). Colonies grow year-round and typically spawn during one season per year ⁴⁴, whilst developing eggs for the rest of the annual gametogenic cycle. Together, these results suggest that corals exist along an energetic continuum, where co-benefits may be derived from underlying physiological drivers like immunity ⁴², feeding efficiency ⁴⁵, or energy storage ⁴⁶. Energetic variability could then result in higher levels of fitness and better performance across suites of different traits. Weak positive associations with heat tolerance occurred with growth but not fecundity, suggesting that the drivers of energy allocation to fecundity may act independently of heat tolerance.

A trade-off between coral heat tolerance and key ecological traits like growth would have considerable negative implications for natural evolution under climate change ¹⁷. This would also apply to restoration efforts involving assisted evolution that aim to boost population resistance to heat stress by propagating more heat tolerant coral individuals via selective breeding or assisted gene flow. If coral heat tolerance was associated with lower growth or fecundity for instance, then out-planting large numbers of corals with these traits would have potentially damaging effects on natural population fitness. However, we found no evidence for such trade-offs between heat tolerance and either colony growth or fecundity, for a coral population associated with the same Symbiodiniaceae community. Although further work will be needed to understand whether these trends persist across larger spatial scales and for other species, our results suggest that selecting corals for heat tolerance through either natural selection or assisted evolution is unlikely to come at a cost to growth or fecundity. Under climate change, coral heat tolerance will be increasingly important to the persistence of coral populations. Contrary to expectations, selection for heat tolerance may not necessarily compromise other important parts of the coral life-history.

4.1 | Model System

To minimise the influence of environmental and interspecific drivers, we measured variation in traits on a single coral population at Mascherchur reef, Palau, Micronesia (7° 17' 29.3'' N, 134° 31' 8.0'' E), a semi-sheltered outer reef crest. The shallow-water Indo-Pacific reef builder, Acropora digitifera, was chosen as a model species due to its high local abundance, broad geographic distribution, and corymbose growth form that allows sub-sampling of branches (fragments) to provide intra-organism statistical replication without sacrificing the colony. Seventy coral colonies were tagged at 2–3 m depth and surveyed repeatedly for different traits between 2017 and 2019. Large adult colonies of similar diameter (24 ± 8 cm, average ± SD) were chosen to limit the size-related variability in total resource budgets among individuals which could obscure trade-off relationships ¹².

4.2 | Marine heatwave emulation experiment

Colony heat tolerance was determined in August 2018 at the Palau International Coral Reef Center by subjecting replicate fragments of each colony to a long-term 5-week marine heatwave emulation experiment (all tank experiment details are given in Table S1). In comparison to short heat-shock experiments that typically last 1–2 days, this experimental temperature profile was designed to match more closely the duration of natural marine heatwaves ^8,33, with the assumption that the phenotypic bleaching and mortality responses would be more ecologically relevant. After fragments were collected and given a 7–10-day acclimatisation period, temperature was increased gradually over the time course of the experiment (35 days), reaching a final bleaching-level temperature of approximately 33 ºC, or 3.5 ºC above the local climatological baseline (Fig. 1A, Table S1). The use of flow-through tank systems allowed an element of natural diel temperature variability ⁴⁷ in all tanks (4 heat stress and 2 procedural control tanks), while aquarium lights provided light conditions close to those measured on the reef ⁸. Calibrated HOBO loggers were placed in each tank and recorded temperatures at 10-minute intervals. To relate coral bleaching and mortality responses to accumulated heat stress, not instantaneous temperature, we calculated heat stress using the Degree Heating Weeks (DHW) metric. DHW was developed by the National Oceanic and Atmospheric Administration’s Coral Reef Watch to provide a real-time coral bleaching alert system based on satellite-derived sea surface temperatures. DHW is a daily measure of both the intensity and duration of heat stress, calculated by accumulating temperature anomalies > 1 ºC relative to a climatological baseline (MMM – maximum of monthly means) over a 12-week (84-day) rolling window ⁴⁸. To allow for comparisons between the DHW within our experiment and the NOAA bleaching forecasts, we adjusted the MMM from the 5 km grid cell encompassing the collection site based on the relationship between daily time series of 5 km sea surface temperature (CoralTemp v3.1) and daily averaged in situ temperature (recorded from additional HOBO loggers at the collection site; ⁸), producing the local climatological baseline (MMM_adj – adjusted MMM). The final accumulated heat stress exposure reached in this experiment was 10.7 ºC-weeks.

4.3 | Heat tolerance

Fragments (6 per colony) were dispersed in random locations among the four heat stress tanks (4 fragments/colony) and two procedural control tanks (2 fragments/colony). If fragments from a colony died in the procedural control tank which was under non-stressful ambient temperature conditions, it was an indication of handling effects for that colony, so all remaining fragments from the colony were removed from the experiment and the colony was not assigned a heat tolerance score. The health status of each fragment was scored visually at intervals of between 1 and 3 days. We used a bleaching survival index (BSI), the inverse of the commonly used bleaching and mortality index (BMI), to categorise coral bleaching and survivorship responses ⁴⁹. This was done in order to have a positive correlation between the bleaching survival index and heat tolerance. For BSI, c₁ to c₅ are the proportion of replicate fragments (per colony) recorded as healthy (c₁), half bleached (c₂), bleached (c₃), partial mortality (c₄), or dead (c₅), and N is the total number of categories (here N = 5). For example, a colony whose replicate fragments are either all healthy or all dead, will have a BSI value of 1 or 0, respectively.

$$BSI=1-BMI$$

$$BMI=\frac{{0c}_{1}+{1c}_{2}+{2c}_{3}+{3c}_{4}+{4c}_{5}}{N-1}$$

Here, we define the onset of the bleaching and mortality response to occur when BSI declines below 0.75. This BSI score was chosen as it represents a colony with either all replicate fragments in the ‘partially bleached’ category, or with some replicate fragments healthy while others from the same colony were fully bleached or dead.

BSI was standardised to the amount of accumulated heat stress (i.e., DHW) as in ⁸. In this way, the bleaching stress response can be compared to publicly available bleaching forecasts to facilitate the interpretation of our results. Notably the NOAA CRW bleaching risk forecast considers DHW of 4 and 8 ºC-weeks as Alert Level 1 (significant bleaching expected) and Alert Level 2 (significant bleaching and mortality expected), respectively ⁴⁸. As the BSI of a particular colony is representative of only a single timepoint and will change throughout the heat stress exposure, the colony’s overall heat tolerance was calculated as the average BSI through time.

4.4 | Symbiont Identification

The composition of the symbiont community was identified from one tissue scraping (< 1 cm) per colony sampled in March 2018 (September 2018 for two colonies not found in March), and is assumed to be stable since no major disturbances occurred during this time period ²⁴. Due to some colonies not being relocated during specific surveys, DNA samples were available for 51 colonies. Polymerase chain reaction (PCR) was used to amplify DNA extracted from coral tissue and then sent for ITS2 sequencing. The hyper conserved ITS2 region was chosen to facilitate integration with the SymPortal database which assigns symbiont ITS2 Type profiles or distinct intragenomic variants (DIVs) that represent different symbiont taxa ⁵⁰. A full description of the DNA extraction, sequencing protocol and Symportal analysis are provided in the supplementary materials (S1).

4.5 | Size and growth metrics

Coral colony size was measured from 3D models at successive timepoints throughout the study (November 2017, May 2018, and February 2019). Overall, 114 3D reconstructions were built for 45 coral colonies. Following Ferrari et al. (2017), Metashape Professional (v 1.5, Agisoft) was used to construct 3D surface meshes using structure from motion photogrammetry, which had > 93% photo alignment on average and an average (± SD) resolution (distance between mesh vertices) of 1.1 (0.4) mm and average (± SD) scaling error of 0.36 (0.38) mm. Geomagic Control (v 2015, 3D Systems) was used to align successive coral colony models to give a common base, and compute three size metrics: live surface area (LiveSA, excluding dead regions without live tissue), total surface area (SA) and total volume (V). The initial maximum diameter (D) was recorded in situ in 2017 for each coral colony using a measurement tape (i.e., independently of the 3D models). Further details of photogrammetry methods can be found in the supplementary materials (S2).

Annual growth rates, in terms of live surface area and volume, were estimated as the change in size relative to the time interval between successive photogrammetry surveys. For those colonies with two size comparisons (i.e., 2017–2018 and 2018–2019), the average of both was taken. Due to some colonies not being relocated during specific surveys or model building issues, photogrammetry model comparisons were available for 45 colonies. However, variation in initial colony size can introduce bias to raw areal or volumetric growth measurements (i.e., for identical linear extension rates, large colonies appear to grow faster than small colonies in terms of surface area and volume). Therefore, growth metrics must be corrected for colony size. One option is to calculate percentage growth, by dividing the growth rate by the corresponding initial colony size. However, this results in growth overestimation for small colonies that may easily double in size from one year to the next (i.e., growth rate of 100%). Another option, which we adopted here, is to multiply the growth rate by a dimensionless adjustment factor (AF) that represents the difference between the population mean colony size and the size of an individual colony. AFs for each colony (i) were calculated as the ratio between the average initial size across all colonies relative to the initial size of colony i, such that

$${Areal\_AF}_{i}= \frac{\stackrel{-}{D }}{{D}_{i}}$$

$${Volumetric\_AF}_{i}= \frac{\stackrel{-}{A }}{{A}_{i}}$$

To avoid further bias, the colony size metrics used for the AF calculation were derived from D independently of the 3D models and were one dimension lower than their corresponding growth metric (e.g., correct the SA growth rate based on the diameter). For three colonies that lacked empirical measurements of D but had 3D photogrammetry models available, D was predicted from the initial SA of the photogrammetry model, such that $\text{l}\text{o}\text{g}\left(D\right)=-0.738+0.537\times \text{l}\text{o}\text{g}\left(initial SA\right)$ (F_40/1=239.1, R² = 0.86, P < 0.001) (Fig. S5). The initial area (A) for use in the volumetric AF was given by the product of D times initial height (assumed to be D/2). Size-dependencies of each metric were tested using linear regressions. Notably, partial mortality was recorded from field surveys categorically, such that a colony without partial mortality must have been healthy continuously from the first until last 3D photogrammetry survey.

4.6 | Fecundity

Fecundity was measured as egg density (number of eggs per polyp) and egg volume (average volume per polyp). Two fragments from each tagged coral colony in the reef were removed prior to spawning in March 2018 or 2019, decalcified using 10% hydrochloric acid and stored in ethanol. Due to some colonies not being relocated during specific surveys or sample preservation issues, polyp fecundity measurements were available for 47 colonies. Ten polyps from each fragment were dissected using a dissection scope and each egg was photographed alongside a scale using an attached digital camera. All image analysis was performed using the semi-automated SizeExtractR workflow ⁵² to annotate polyp images and compute egg density and volume from the 20 polyps per colony. Egg volume (EV) was estimated from the geometric mean diameter (GMD) ⁵². We then estimated total colony egg production (TEP) and total colony egg volume (TEV) by combining dissection-derived fecundity data (EC – egg counts per polyp – and EV, respectively) with per colony estimates of polyp density per unit area (PD) and 3D model-derived estimates of total live surface area (LiveSA).

$$EV= \frac{4}{3} \times \pi \times {\left(\frac{GMD}{2}\right)}^{3}$$

$$TEP=EC\times PD\times LiveSA$$

$$TEV=EV \times TEP$$

4.7 | Statistical approach

Traits hypothesised to influence coral heat tolerance (average BSI) were colony growth (in terms of both LiveSA and V), colony fecundity (total egg production and total egg volume), and symbiont community. Trade-off analyses were conducted by regressing predictor variables against colony heat tolerance using Bayesian general linear models. In these regressions a negative slope represents a trait trade-off while a positive slope shows co-benefits among traits. Bayesian models were used to facilitate a more intuitive interpretation of model uncertainty ^53,54 and fit using integrated nested LaPlace approximation in the statistical package R-INLA ⁵⁴. The effect of multivariate symbiont community type on sigmoidal BSI-DHW responses was tested using generalised linear models with a binomial response distribution and post-hoc Tukey tests. In addition, confounding effects of partial mortality on all traits were tested using ANOVA paired to post-hoc Tukey tests or non-parametric equivalents.

During a long-term heat stress exposure, the bleaching and mortality response of sampled fragments may vary through the time course of the experiment. Therefore, we were interested to test not only the relationship between heat tolerance (BSI averaged across the whole heat stress exposure) and other traits, but also the relationship at individual timepoints between the instantaneous BSI based on replicate fragments (response variable) and the known colony-level traits (predictor variable, e.g., fecundity or colony growth). To link bleaching responses to their ultimate driver – heat stress (not time) – we evaluated BSI-trait relationships at specific DHW levels instead of at survey time points. For each DHW level, we computed the linear regression between the instantaneous BSI of each colony against their corresponding colony-level trait of interest. This allowed us to evaluate how the trait relationships change throughout the duration of the heat stress exposure, shown as changes in BSI-trait regression slopes. All LMs were fit using Bayesian inference in R-INLA ⁵⁴, and uncertainty in regression slope estimations were quantified using Bayesian credible intervals (95%, 75%, 50%).

Acknowledgments

We thank the numerous staff at the Palau International Coral Reef Centre (PICRC) who supported this research, Dr. Jamie Craggs for contributions to the experimental aquaria system, and Faith Paysinger for work on fecundity measurement. This research was funded by Natural Environment Research Council’s ONE Planet Doctoral Training Partnership (NE/S007512/1) to LL and the European Research Council Horizon 2020 project CORALASSIST (725848) awarded to JRG.

Author Contributions

The first draft of the manuscript, 3D model comparisons, and statistical analyses were conducted by LL. The marine heatwave emulation experiment was conducted by AH, HM, JB, and JRG, AJE, ES. 3D photogrammetry surveys were conducted by RF, WFF, DRP, BS, and JRG, and photogrammetry models were constructed by DRP, RF, WFF, and BS. Symbiont DNA extractions were conducted by JCB and AH. All authors contributed to the development of ideas and writing the final manuscript.

Data Availability Statement

All original data and R code produced in this study are publicly available on Figshare at 10.25405/data.ncl.20411589 (until publication, the repository is link-access only and can be found at https://figshare.com/s/f1016b972bd2f159faa2). ITS2 sequences have been archived publicly at NCBI under BioProject 864615 and processed symbiont community composition can be explored publicly at https://symportal.org. Any additional information required to reanalyse the data reported in this paper is available from the lead contact upon request.

Conflicts of interest

The authors have no conflicts of interest.

Hughes, T.P., et al.: Global warming and recurrent mass bleaching of corals. Nature. 543, 373–377 (2017)
Cheung, M.W.M., Hock, K., Skirving, W., Mumby, P.J.: Cumulative bleaching undermines systemic resilience of the Great Barrier Reef. Curr. Biol. 1–8 (2021). doi:10.1016/j.cub.2021.09.078
Lachs, L., et al.: Linking population size structure, heat stress and bleaching responses in a subtropical endemic coral. Coral Reefs. 40, 777–790 (2021)
Li, X., Donner, S.D.: Lengthening of warm periods increased the intensity of warm-season marine heatwaves over the past 4 decades. Clim. Dyn. (2022). doi:10.1007/s00382-022-06227-y
Oliver, E.C.J., et al.: Longer and more frequent marine heatwaves over the past century. Nat. Commun. 9, 1–12 (2018)
Claar, D.C., et al.: Dynamic symbioses reveal pathways to coral survival through prolonged heatwaves. Nat. Commun. 11, 1–10 (2020)
Buddemeier, R.W., Fautin, D.G.: Coral Bleaching as an Adaptive Mechanism. Bioscience. 43, 320–326 (1993)
Humanes, A., et al.: Within-population variability in coral heat tolerance indicates climate adaptation potential. Proc. R. Soc. B Biol. Sci. 20220872 doi: (2022). 10.1098/rspb.2022.0872
van Oppen, M.J.H., et al.: Shifting paradigms in restoration of the world’s coral reefs. Glob Chang. Biol. 23, 3437–3448 (2017)
Humanes, A., et al.: An Experimental Framework for Selectively Breeding Corals for Assisted Evolution. Front. Mar. Sci. 8, 1–18 (2021)
Martin, T.E., Mouton, J.C.: Longer-lived tropical songbirds reduce breeding activity as they buffer impacts of drought. Nat. Clim. Chang. 10, 953–958 (2020)
Smallegange, I.M.: Life History Trade-offs. Encycl Evol. Biol. 2, 390–393 (2016)
Vinebrooke, R.D., et al.: Impacts of multiple stressors on biodiversity and ecosystem functioning: The role of species co-tolerance. Oikos. 104, 451–457 (2004)
Ernande, B., Boudry, P., Clobert, J., Haure, J.: Plasticity in resource allocation based life history traits in the Pacific oyster, Crassostrea gigas. I. Spatial variation in food abundance. J. Evol. Biol. 17, 342–356 (2004)
Cantin, N.E., Van Oppen, M.J.H., Willis, B.L., Mieog, J.C., Negri, A.P.: Juvenile corals can acquire more carbon from high-performance algal symbionts.Coral Reefsvol. 28 (2009)
Jones, A., Berkelmans, R.: The photokinetics of thermo-tolerance in Symbiodinium. Mar. Ecol. 33, 490–498 (2012)
Ortiz, J.C., González-Rivero, M., Mumby, P.J.: Can a thermally tolerant symbiont improve the future of Caribbean coral reefs? Glob Chang. Biol. 19, 273–281 (2013)
Little, A.F., Van Oppen, M.J.H., Willis, B.L.: Flexibility in algal endosymbioses shapes growth in reef corals. Sci. (80-). 304, 1492–1494 (2004)
Jones, A., Berkelmans, R.: Potential costs of acclimatization to a warmer climate: growth of a reef coral with heat tolerant vs. sensitive symbiont types.PLoS One5, (2010)
Cunning, R., Gillette, P., Capo, T., Galvez, K., Baker, A.C.: Growth tradeoffs associated with thermotolerant symbionts in the coral Pocillopora damicornis are lost in warmer oceans. Coral Reefs. 34, 155–160 (2015)
Silverstein, R.N., Cunning, R., Baker, A.C.: Change in algal symbiont communities after bleaching, not prior heat exposure, increases heat tolerance of reef corals. Glob Chang. Biol. 21, 236–249 (2015)
Matthews, J.L., et al.: Metabolite pools of the reef building coral Montipora capitata are unaffected by Symbiodiniaceae community composition. Coral Reefs. (2020). doi:10.1007/s00338-020-01999-3
Thomas, L., et al.: Spatially varying selection between habitats drives physiological shifts and local adaptation in a broadcast spawning coral on a remote atoll in Western Australia. 9185,1–12(2022)
Jung, E.M.U., Stat, M., Thomas, L., Schoepf, V.: Coral host physiology and symbiont dynamics associated with differential recovery from mass bleaching in an extreme, macro-tidal reef environment in northwest Australia. Coral Reefs. (2021). doi:10.1007/s00338-021-02094-x
Anderson, A.R., Collinge, J.E., Hoffmann, A.A., Kellett, M., McKechnie, S.W.: Thermal tolerance trade-offs associated with the right arm of chromosome 3 and marked by the hsr-omega gene in Drosophila melanogaster. Heredity (Edinb). 90, 195–202 (2003)
Kelly, M.W., Grosberg, R.K., Sanford, E.: Trade-offs, geography, and limits to thermal adaptation in a tide pool copepod. Am. Nat. 181, 846–854 (2013)
Roze, T., Christen, F., Amerand, A., Claireaux, G.: Trade-off between thermal sensitivity, hypoxia tolerance and growth in fish. J. Therm. Biol. 38, 98–106 (2013)
Howells, E.J., Berkelmans, R., Van Oppen, M.J.H., Willis, B.L., Bay, L.K.: Historical thermal regimes define limits to coral acclimatization. Ecology. 94, 1078–1088 (2013)
Paxton, C.W., Baria, M.V.B., Weis, V.M., Harii, S.: Effect of elevated temperature on fecundity and reproductive timing in the coral Acropora digitifera. Zygote. 24, 511–516 (2016)
Michalek-Wagner, K., Willis, B.L.: Impacts of bleaching on the soft coral Lobophytum compactum. I. Fecundity, fertilization and offspring viability. Coral Reefs. 19, 231–239 (2001)
Okubo, N., Taniguchi, H., Motokawa, T.: Successful methods for transplanting fragments of Acropora formosa and Acropora hyacinthus. Coral Reefs. 24, 333–342 (2005)
Okubo, N., Motokawa, T., Omori, M.: When fragmented coral spawn? Effect of size and timing on survivorship and fecundity of fragmentation in Acropora formosa. Mar. Biol. 151, 353–363 (2007)
McLachlan, R.H., Price, J.T., Solomon, S.L., Grottoli, A.: Thirty years of coral heat-stress experiments: a review of methods. Coral Reefs. (2020). doi:10.1007/s00338-020-01931-9
Ayalon, I., et al.: Report Coral Gametogenesis Collapse under Artificial Light Pollution ll ll Report Coral Gametogenesis Collapse under Artificial Light Pollution. Curr. Biol. 31, 413–419e3 (2021)
Marshall, A.T., Clode, P.: Calcification rate and the effect of temperature in a zooxanthellate and an azooxanthellate scleractinian reef coral. Coral Reefs. 23, 218–224 (2004)
Silbiger, N.J., Goodbody-Gringley, G., Bruno, J.F., Putnam, H.M.: Comparative thermal performance of the reef-building coral Orbicella franksi at its latitudinal range limits.Mar. Biol.166, (2019)
Donner, S.D., Skirving, W.J., Little, C.M., Oppenheimer, M., Hoegh-Gulberg, O.: Global assessment of coral bleaching and required rates of adaptation under climate change. Glob Chang. Biol. 11, 2251–2265 (2005)
Quigley, K.M., Bay, L.K., Oppen, M.J.H.: The active spread of adaptive variation for reef resilience. Ecol. Evol. 1–14 (2019). doi:10.1002/ece3.5616
Jones, A.M., Berkelmans, R.: Tradeoffs to Thermal Acclimation: Energetics and Reproduction of a Reef Coral with Heat Tolerant Symbiodinium Type-D. J. Mar. Biol. (2011). (2011)
Ros, M., et al.: Symbiont shuffling across environmental gradients aligns with changes in carbon uptake and translocation in the reef-building coral Pocillopora acuta. Coral Reefs. 40, 595–607 (2021)
Kaufman, M.L., Watkins, E., van Hooidonk, R., Baker, A.C., Lirman, D.: Thermal history influences lesion recovery of the threatened Caribbean staghorn coral Acropora cervicornis under heat stress. Coral Reefs. 40, 289–293 (2021)
Palmer, C.V., Traylor-Knowles, N.G., Willis, B.L., Bythell, J.C.: Corals use similar immune cells and wound-healing processes as those of higher organisms. PLoS One. 6, 1–11 (2011)
Baria, M.V.B., Cruz, D.W., Villanueva, R.D., Guest, J.R.: Spawning of three-year-old Acropora millepora corals reared from larvae in Northwestern Phillipines. Bull. Mar. Sci. 88, 61–62 (2012)
Guest, J.R., Baird, A.H., Goh, B.P., Chou, L.M.: Seasonal reproduction in equatorial reef corals. Invertebr Reprod. Dev. 48, 207–218 (2005)
Anthony, K.R.N.: Enhanced particle-feeding capacity of corals on turbid reefs (Great Barrier Reef, Australia). Coral Reefs. 19, 59–67 (2000)
Hinrichs, S., et al.: Seasonal variations in energy levels and metabolic processes of two dominant Acropora species (A. spicifera and A. digitifera) at Ningaloo Reef. Coral Reefs. 32, 623–635 (2013)
Craggs, J., et al.: Inducing broadcast coral spawning ex situ: Closed system mesocosm design and husbandry protocol. Ecol. Evol. 7, 11066–11078 (2017)
Skirving, W., et al.: Coraltemp and the coral reef watch coral bleaching heat stress product suite version 3.1. Remote Sens. 12, 1–10 (2020)
McClanahan, T.R., Baird, A.H., Marshall, P.A., Toscano, M.A.: Comparing bleaching and mortality responses of hard corals between southern Kenya and the Great Barrier Reef, Australia. Mar. Pollut Bull. 48, 327–335 (2004)
Hume, B.C.C., et al.: SymPortal: A novel analytical framework and platform for coral algal symbiont next-generation sequencing ITS2 profiling. Mol. Ecol. Resour. 19, 1063–1080 (2019)
Ferrari, R., et al.: 3D photogrammetry quantifies growth and external erosion of individual coral colonies and skeletons. Sci. Rep. 7, 1–9 (2017)
Lachs, L., et al.: SizeExtractR: a workflow for rapid reproducible extraction of object size metrics from scaled images.Ecol. Evol.12, (2022)
Lachs, L., et al.: Fine-tuning heat stress algorithms to optimise global predictions of mass coral bleaching.Remote Sens.13, (2021)
Bakka, H., et al.: Spatial modeling with R-INLA: A review. Wiley Interdiscip. Rev. Comput. Stat. 10, 1–24 (2018)

There is NO Competing Interest.

LachsSupplementaryMaterial.docx

Download PDF

Journal Publication

published 12 Apr, 2023

Read the published version in Communications Biology →

Version 1

posted

You are reading this latest preprint version

Co-benefits not trade-offs associated with heat tolerance in a reef building coral

Status:

Journal Publication

Version 1

Abstract

Figures

1 | Introduction

2 | Results

2.1 | Heat tolerance variability

2.2 | Similar Symbiodiniaceae communities

2.3 | 3D colony growth

2.4 | Fecundity

2.5 | Co-benefits but not trade-offs with heat tolerance

3 | Discussion

4 | Materials And Methods

4.1 | Model System

4.2 | Marine heatwave emulation experiment

4.3 | Heat tolerance

4.4 | Symbiont Identification

4.5 | Size and growth metrics

4.6 | Fecundity

4.7 | Statistical approach

Declarations

References

Additional Declarations

Supplementary Files

Status:

Journal Publication

Version 1