Blog

Apr 28, 2023

Acclimatization of a coral

Communications Biology volume

Communications Biology volume 6, Article number: 66 (2023) Cite this article

1080 Accesses

1 Citations

23 Altmetric

Metrics details

Ocean acidification caused by shifts in ocean carbonate chemistry resulting from increased atmospheric CO2 concentrations is threatening many calcifying organisms, including corals. Here we assessed autotrophy vs heterotrophy shifts in the Mediterranean zooxanthellate scleractinian coral Balanophyllia europaea acclimatized to low pH/high pCO2 conditions at a CO2 vent off Panarea Island (Italy). Dinoflagellate endosymbiont densities were higher at lowest pH Sites where changes in the distribution of distinct haplotypes of a host-specific symbiont species, Philozoon balanophyllum, were observed. An increase in symbiont C/N ratios was observed at low pH, likely as a result of increased C fixation by higher symbiont cell densities. δ13C values of the symbionts and host tissue reached similar values at the lowest pH Site, suggesting an increased influence of autotrophy with increasing acidification. Host tissue δ15N values of 0‰ strongly suggest that diazotroph N2 fixation is occurring within the coral tissue/mucus at the low pH Sites, likely explaining the decrease in host tissue C/N ratios with acidification. Overall, our findings show an acclimatization of this coral-dinoflagellate mutualism through trophic adjustment and symbiont haplotype differences with increasing acidification, highlighting that some corals are capable of acclimatizing to ocean acidification predicted under end-of-century scenarios.

The defining feature of the Anthropocene era1 is the emergence of human activities as a driving force of global change2, which is occurring at a rate that raises concerns whether organismal adaptation can keep pace with rapidly changing environmental conditions3. Stressors associated with ocean warming and acidification are among the most direct and pervasive anthropogenic changes for marine biota, including corals4. The decrease in pH from ca. 8.2 prior to the industrial revolution to ca. 8.1 with a doubling of CO2 is leading to a gradual decline in the saturation state of calcium carbonate in seawater5. This phenomenon has been projected to negatively impact the ability of corals to calcify6,7. However, empirical evidence, in terms of e.g., natural selection of tolerant bacterial and/or dinoflagellate symbionts, differential regulation of environmental stress response genes8,9, suggests an underappreciated capacity of corals to acclimate and genetically adapt to environmental change10. Indeed, these ancient organisms have survived, evolved, and adapted throughout hundreds of millions of years of global climate change11,12,13,14.

The symbiosis between scleractinian corals and their dinoflagellate microalgae (family Symbiodiniaceae), commonly referred to as zooxanthellae, has been extensively studied15,16,17. Zooxanthellae significantly contribute to the energy budget of the host by providing photosynthetically fixed carbon16 while recycling host respiration and excretion by-products18. In such symbiosis, both carbon and nitrogen can be obtained via heterotrophy and autotrophy and are recycled between the host and dinoflagellate symbionts19,20. Generally, the symbiosis provides most of the carbon needed for respiration21, while predation on zooplankton and particulate organic matter are still needed to meet nitrogen and phosphorus requirements16. Nonetheless, the relative contribution of heterotrophy versus autotrophy on host nutrition varies across species, populations, environments, and/or with ontogenesis22.

A critical limitation for many experimental studies has been replicating the rate (decades) and biological scales (ecosystems) at which ocean acidification operates. Natural CO2 vents acidify the surrounding seawater creating carbonate chemistry conditions that mimic future ocean acidification predictions23,24,25, even if with wide short-term variability26,27,28. By investigating natural populations living along transects radiating from the sources of CO2, these systems allow substituting time for space, providing valuable insights on acclimatization and adaptation to ocean acidification29. This study was conducted on natural populations of the zooxanthellate scleractinian coral Balanophyllia europaea living at a CO2 volcanic vent near Panarea Island (Italy). This underwater crater at 10 m depth releases persistent gaseous emissions (98–99% CO2 without instrumentally detectable toxic compounds), resulting in a stable pH gradient at ambient temperature30 with ocean acidification conditions projected for 2100 under conservative and worst-case IPCC scenarios31,32. With decreasing pH, B. europaea shows a decline in population density33 and net calcification rates, the latter as a result of increased skeletal porosity, whereas linear extension rate is preserved24, allowing the coral to reach size at sexual maturity28. Moreover, B. europaea maintains unchanged skeletal calcium carbonate polymorph, organic matrix content, aragonite fiber thickness and skeletal hardness, calcifying fluid pH, and gross calcification with decreasing pH34.

The aim of this study was to further investigate the physiological acclimatization of this coral-dinoflagellate mutualism with increasing acidification, by assessing the relative contribution of photosynthetically versus heterotrophically derived nutrition (autotrophic/heterotrophic ratio). This was compared with other key physiological parameters (i.e., photosynthetic efficiency, symbiont cell densities, and chlorophyll concentration) and with the distribution of genetically distinct haplotypes of the dinoflagellate symbiont Philozoon balanophyllum along the gradient. We expected to find an overall physiological acclimatization of the symbiont, allowing the coral to sustain energetically costly processes (e.g., maintaining elevated calcifying fluid pH, reproduction) that in previous studies were found to be constant along the natural pH gradient24,28,34.

Among measured seawater parameters (pH, total alkalinity, temperature, and salinity), only pH changed across Sites (Kruskal–Wallis test, H = 38, df = 2, p = 0.000; Supplementary Table 1). Even though fluctuations were observed for pH, the average value decreased from 7.97 (95% CI: 7.95–7.99) at Site 1, to 7.86 (95% CI: 7.83–7.90) at Site 2, to 7.64 (95% CI: 7.55–7.77) at Site 3 (Mann-Whitney U test, p = 0.000; Fig. 1; Supplementary Table 1). Average aragonite saturation (Ωarag) decreased by almost 30% from Site 1 (average: 3.2, 95%CI: 3.1–3.4) to Site 3 (average=2.3; 95% CI: 2.2–2.4) (Mann–Whitney U test, p = 0.000; Fig. 1; Supplementary Table 1). No significant difference in Ωarag was observed between Sites 2 and 3 (Mann–Whitney U test, p = 0.432; Supplementary Table 1).

The boxes indicate the 25th and 75th percentiles and the line within the boxes mark the medians. Whisker length is equal to 1.5 × interquartile range (IQR). Circles represent outliers. Different letters indicate statistical differences (p < 0.05; number of observations: pHTS = 185-198 per Site; Ωarag = 184–195 per Site).

Phosphate and inorganic nitrogen (Nitrate+Nitrite) were homogeneous across Sites (Kruskal–Wallis test, H = 1.192, df = 2, p = 0.551 and H = 3.082, df = 2, p = 0.214, respectively), while sulphate was significantly lower at Site 3 compared to Sites 1 and 2 (Mann–Whitney U test, p = 0.009 and p = 0.016, respectively; Supplementary Table 2). No significant difference in sulphate was observed between Sites 1 and 2 (Mann-Whitney U test, p = 0.293; Supplementary Table 2).

Photosynthetic yield of the dinoflagellate symbionts, measured on 7–25 haphazardly chosen corals per Site (Supplementary Table 3), did not change along the gradient in any of the assessed time intervals. However, yield was significantly higher during the low light time intervals 18:00–20:00 and 20:00–22:00 than in 9:00-13:00 (Fig. 2; Supplementary Tables 4 and 5). Minimum (F) and maximum fluorescence (Fm’) were significantly higher at Sites 2 and 3 compared to Site 1 at all time intervals. Minimum fluorescence (F) increased from 20:00–22:00 to 9:00–13:00 to 18:00–20:00, while maximum fluorescence (Fm’) was significantly higher in the time interval 18:00–20:00 than in the other time intervals (Fig. 2; Supplementary Tables 4 and 5).

a Effective quantum yield (ΔF/Fm´), b minimum fluorescence (F), and c maximum fluorescence (Fm´). For each parameter, significant differences (p < 0.01) between time intervals are indicated by different colors (significantly higher mean value is ranked as white->light gray->dark grey). The boxes indicate the 25th and 75th percentiles and the line within the boxes mark the medians. Whisker length is equal to 1.5 × interquartile range (IQR). Circles represent outliers. Different letters indicate statistical differences (p < 0.05; number of specimens measured is reported in Supplementary Table 3).

Increasing numbers of symbiont cells from Site 1 to Site 3 were observed along the gastroderm (endoderm) of the mesenteries in histological sections (Fig. 3a–f). Indeed, the average number of symbiont cells per area was higher at Sites 2 and 3 compared to Site 1 (LSD post hoc test: p = 0.002; Fig. 3g). No significant difference in the number of symbiont cells was observed between Sites 2 and 3 (LSD post hoc test: p = 0.968; Fig. 3g). Chl-a concentration per symbiont cell was homogenous among studied Sites (One-way ANOVA, F2,9 = 0.359, p = 0.708; Fig. 3h). The average chl-a concentration, expressed as chl-a amount per polyp area (pg/mm2), was higher at Sites 2 and 3 compared to Site 1 (LSD post hoc test: p = 0.015 and p = 0.053, respectively; Fig. 3i). No significant difference in chl-a amount per polyp area was observed between Sites 2 and 3 (LSD post hoc test: p = 0.449; Fig. 3i).

a–c Histological transverse sections of representative polyps in Sites 1-3 (x 4 magnification). d–f Symbiont cells within the endoderm of the mesentery highlighted by x40 magnification. z: zooxanthellae; ec: ectoderm; m: mesoglea; en: endoderm. g Symbiont cell density and h chlorophyll-a (chl-a) concentration normalized over cell count or i polyp area. The boxes indicate the 25th and 75th percentiles and the line within the boxes mark the medians. Whisker length is equal to 1.5 × interquartile range (IQR). Circles represent outliers. Different letters indicate statistical differences (p < 0.05; number of corals = 4 per Site).

The δ13C values for both dinoflagellate symbionts and host tissue did not show significant statistical differences along the gradient (symbionts: Kruskal-Wallis test, H = 2.346, df = 2, p = 0.309; tissue: One-way ANOVA, F2,9 = 1.243, p = 0.334; Fig. 4; Supplementary Table 6). However, at Site 3 (pH=7.64, pCO2 = 1161) δ13Ch and δ13Cz reached similar values in some specimens (Fig. 4; Supplementary Fig. 1) and the proportion of carbon in the host supplied by the symbionts35 reached the highest values (Fig. 4).

The proportion of carbon in the host supplied by the symbionts was calculated using the formula: \({\delta }^{13}{{{{{{\rm{C}}}}}}}_{{{{{{\rm{tissue}}}}}}}=({\delta }^{13}{{{{{{\rm{C}}}}}}}_{{{{{{\rm{symbiont}}}}}}}){{{{{\rm{x}}}}}}+(1-{{{{{\rm{x}}}}}})({\delta }^{13}{{{{{{\rm{C}}}}}}}_{{{{{{\rm{zooplankton}}}}}}/{{{{{\rm{POC}}}}}}})\). We assume δ13Czooplankton/POC = −22‰37. The boxes indicate the 25th and 75th percentiles and the line within the boxes mark the medians. Whisker length is equal to 1.5 × interquartile range (IQR). Different letters indicate statistical differences (p < 0.05) between Sites (number of corals = 4 per Site).

δ15N values in symbionts were homogeneous among Sites (One-way ANOVA, F2,9 = 1.001, p = 0.405; Fig. 4; Supplementary Table 6), whereas δ15N in coral tissue showed heavier values in Site 1 compared to Sites 2 and 3 (LSD post hoc test: p = 0.000). No significant difference in host tissue δ15N was observed between Sites 2 and 3 (LSD post hoc test: p = 0.181).

Host tissue C/N ratio was significantly lower in Sites 2 and 3 compared to Site 1 (LSD post hoc test, p = 0.000; Fig. 4; Supplementary Table 6). No significant difference in host tissue C/N ratio was observed between Sites 2 and 3 (LSD post hoc test: p = 0.237). Symbiont C/N ratios increased from Site 1 to Sites 2 and 3 (LSD post hoc test, Site 1 vs Site 2: p = 0.005, Site 1 vs Site 3: p = 0.000, Site 2 vs Site 3, p = 0.002; Fig. 4; Supplementary Table 6).

Analysis of the psbAncr identified Philozoon balanophyllum in each sample of B. europaea analyzed (N = 4 per Site) (Supplementary Dataset 1). Moreover, each individual B. europaea harboured a single P. balanophyllum haplotype. When compared phylogenetically with sequences obtained from P. balanophyllum from other regions in the Tyrrhenian and Adriatic Seas, psbAncr haplotypes from the acidified Sites 2 and 3 grouped to a distinct genetic lineage statistically distinct from all other haplotypes representative of the species (Fig. 5: Supplementary Dataset 1). Haplotype sequences from the less acidified environment of Site 1 and from samples obtained around the Tyrrhenian and Adriatic Sea were more similar to each other (Fig. 5).

Philozoon balanophyllum halpotypes in B. europaea specimens from Sites 1-3 along the Panarea pH gradient and from around the Tyrrhenian and Adriatic Seas. Haplotypes of Philozoon actinarum, the symbiont of the sea anemone Anemonia viridis, common to the Mediterranean, are provided as the outgroup. (Photograph by Francesco Sesso).

As calcifying organisms, corals are threatened by increasing absorbance of atmospheric CO2 into Earth's oceans. Our research indicates that populations of B. europaea adjust to decreasing pH through changes in symbiont cell densities and their haplotype identities, which appear to affect the amount of photosynthetically fixed carbon relative to heterotrophically acquired carbon in the host. Therefore, certain coral-dinoflagellate mutualisms have the capacity to acclimatize to ocean acidification by various ways.

Balanophyllia europaea's relative dependence on autotrophy and heterotrophy for nutrient acquisition appears in part contingent on ocean pH. The δ13C values in host tissues and dinoflagellate symbionts are indicative of greater photosynthetic inputs at the low pH/high pCO2 Sites36,37. The higher seawater pCO2 values at Site 3 probably account for lighter δ13C values in the symbionts22. The latter result is consistent with previous findings reporting a greater net autotrophic input to the carbon budget under low pH/high pCO2 conditions in the temperate sea anemone Anemonia viridis investigated at a similar CO2 vent system at Vulcano Island (Italy)37. Moreover, greater autotrophy through increased symbiont abundances, enhanced C fixation, and less heterotrophically derived nitrogen may also explain the increase in carbon relative to nitrogen concentration in dinoflagellate symbiont cells from Sites 2 and 3 compared to Site 138. Symbiont cell densities and host tissue thickness are known to vary in corals photo-acclimating to seasonal changes in solar irradiance39,40,41. The increase in host tissue (endoderm) thickness from ~15 µm in Site 1 (pH 8.0) to ~40 µm in Site 3 (pH 7.6) and the corresponding increase in symbiont cell densities is similar to winter phenotypes observed in tropical corals40. Tissue thickening has been observed in other coral species experiencing low seawater pH conditions42,43 and may therefore constitute a general phenotypic response to ocean acidification. Moreover, enhanced translocation of photosynthetically fixed carbon to the host as a result of high symbiont densities could contribute to tissue growth and increased energy reserves. Indeed, previous findings have shown that energy derived from photosynthesis may have up to 20-fold greater effect on tissue production compared to skeleton deposition44.

Zooxanthellate corals typically associate with dinoflagellate species that raise the functional tolerance of the mutualism when existing in thermally stressful environments45,46. Physiological constraints related to the pH gradient across the study transect may explain the non-random distribution of distinct haplotypes of Philozoon balanophyllum (formerly temperate Clade A) in B. europaea. The small but fixed differences in psbAncr between populations over a distance of 30 meters or less is greater than the genetic differences in this symbiont distributed over many hundreds of kilometers across the Mediterranean basin (Fig. 5). It is compelling to assume the related haplotypes at Sites 2 and 3 may have attributes (allelic combinations) better fit to the lower pH conditions, possibly contributing to the increase in symbiont cell density observed at the low pH Sites compared to Site 1. However, this possibility will require additional physiological characterization47.

Tissue C/N ratios in the low pH/high pCO2 Sites were unusually low compared to ratios from the control Site (Fig. 4). The latter were consistent with values typically found in marine organisms (from 6.5 to 8.7)48. C/N ratios can reach values around 2 in marine particulate organic matter49,50 and even drop below 1 in sediments characterized by significant amounts of fixed nitrogen51. Surprisingly, we found that the C/N ratio of the symbionts followed an opposite trend compared to the tissue, increasing with decreasing pH. Indeed, assuming that the elemental composition of an animal reflects that of its diet52, this observation may, at first, appear contradictory: If the algae feed the host, and the symbiont's C/N ratio rises, should not the host's C/N ratio follow suit? However, this logic negates the potential for the host to receive nutritional input from multiple sources. Many corals harbor bacteria capable of fixing N2 gas (diazotrophs)53 and a growing body of literature indicates that ocean acidification may promote the enrichment of dinitrogen fixing bacterial communities in natural microbiomes8,54,55. Indeed, recent work at our study site found that the prevalence of genes associated with N2 fixation, as well as the production of N storage molecules, increases in the B. europaea microbiome under low pH conditions56. Increased assimilation of diazotroph-derived N by the coral is supported by our isotopic data (Fig. 4). At low pH Sites, the δ15N of host tissue was near 0‰. This value is significantly lower than typically reported for zooxanthellate coral species57 and strongly suggests that diazotroph N2 fixation is occurring within the coral tissue/mucus58,59. Transcriptomic and proteomic studies aimed at identifying expressed proteins involved in nitrogen fixation (e.g., nitrogenase)60 as well as isotope tracer experiments (e.g., 15N-labeled dinitrogen) aimed at quantifying N2 fixation61 are needed to have a better understanding of this process in B. europaea along the Panarea pH gradient.

Figure 6 summarizes the main skeletal and physiological changes displayed by B. europaea acclimatized to low pH/high pCO2 conditions at the Panarea CO2 vent. While exposure to low pH creates more porous skeletons, linear extension rates are maintained24, allowing corals to reach size at sexual maturity and reproduce normally28, despite possessing a more fragile skeleton24. The acquisition of symbiont strains potentially better adapted to acidified conditions could further contribute to increased nutrient cycling and animal growth. Bacteria involved in nutrient-cycling and N2 fixation may play a crucial role in supplementing coral-Symbiodiniaceae symbiosis with additional nitrogen under ocean acidification, helping sustain the higher photosynthetic rates expected under acidified conditions62,63,64. Likewise, it has been hypothesized that photosynthetically fixed carbon supplied by the dinoflagellates may serve as an energy source for N2 fixation in cyanobacterial symbionts in corals60, suggesting that in zooxanthellate corals, photosynthesis and N2 fixation may be, in some cases, interdependent. Moreover, enhanced N2 fixation by diazotrophs living in association with the coral tissue/mucus56 may partially explain the observed increase in dinoflagellate symbiont cell densities in B. europaea at the low pH Sites65, and viceversa. Enhanced translocation of photosynthetically fixed carbon as a result of high dinoflagellate symbiont cell densities may help sustain the high costs of diazotroph N2 fixation under increased acidification. Indeed, dinoflagellates supply glycerol as an energy source for cyanobacterial symbionts found in Montipora cavernosa, providing a steady supply of reductant (e.g., NADPH) and ATP for dinitrogen fixation in the cyanobacteria60. Additionally, N2 fixation by coral-associated diazotrophs, coral calcification, reproduction, and tissue thickening are energy-intensive processes that probably compete for energy deriving from photosynthesis within the coral holobiont66. Thus, massive N2 fixation, together with maintenance of coral reproduction and thickening of coral tissue occurring at the low pH Sites could be absorbing most of the photosynthetic energy from the symbiotic algae, creating an energy deficit that could possibly explain the previously observed decline in net calcification rates of this species along the same pH gradient24,34. Taken together, current and previous findings highlight the importance of the interactions among all the components of the holobiont to unveil how and to what degree corals will endure ocean acidification predicted for the end of the century.

Under low pH conditions, coral population density decreases33, net calcification is depressed, while linear extension rate is maintained constant, allowing the coral to reach critical size at sexual maturity and reproduce. Moreover, an overall rearrangement of the coral holobiont is documented by: (i) an increase in symbiont cell density, triggered by thickening of the coral tissue and the establishment of novel dinoflagellate haplotypes possibly better adapted to lower pH conditions, and (ii) an increase, within the coral tissue/mucus in microbial communities capable of dinitrogen fixation as well as N storage and mobilization. Variations displayed by corals living at average pH 7.6 compared to corals at average pH 8.0 are shown with grey symbols and are listed in the following order: skeleton (from micro to macro), symbionts (from micro to macro), trophic strategy, reproduction, tissue/mucus microbial community. Darker and lighter shades of green in the coral sketches represent higher and lower symbiont cell density. The color scale bar highlighting the seawater pH change across the gradient does not match the color scale of pH test strips. Image was assembled using Adobe Photoshop CC (19.1.6).

Temperature, salinity, and pH (NBS scale) were measured at three Sites respectively 34 m (Site 1), 13 m (Site 2) and 9 m (Site 3) away from the center of the crater with a multi-parametric probe (600 R, YSI Incorporated) powered from a small boat and operated by SCUBA divers. Total alkalinity was determined by Gran titration (888 Titrando) from bottom-water samples collected at the three Sites and poisoned by adding 1% of saturated HgCl2 shortly after collection33. Certified reference materials (Batch 187) provided by Andrew Dickson (Scripps Institution of Oceanography, La Jolla, CA) were used to ascertain the quality of results obtained. Environmental data were collected during several expeditions between 2010–201324,32,33 and 2019-2020 (this study). Measured pHNBS were converted to the total scale using CO2SYS software67. Mean pHTS (back-transformed hydrogen ion concentrations) was calculated for each Site and used with total alkalinity, salinity and temperature to calculate carbonate chemistry parameters using the software CO2SYS with referenced dissociation constants68,69,70. Bottom water samples for dissolved inorganic nutrients were collected at the four Sites using 100 ml plastic bottles (two replicates for each Site) and frozen at −20 °C. Inorganic nitrogen (nitrite-NO2 + nitrate-NO3) and phosphate (orthophosphate-PO4) were determined using a colorimetric method71,72. Absorbances were measured with an AxFlow quAAtro AutoAnalyzers (N-detection limit/sensitivity: 0.006/0.001 µM; P-detection limit/sensitivity: 0.006/0.001 µM). Sulphate was measured by ion chromatography using a Metrohm 761 compact IC.

Photosynthetic yield was estimated by measuring PSII (photosystem II) fluorescence using a Diving PAM (Pulse Amplitude Modulator) (Heinz Walz GmbH). Measurements were taken in July 2019 and February 2020 on a total of 353 B. europaea specimens naturally living at Sites 1-3 at the Bottaro CO2 vents off Panarea Island24,33. During each dive, 7–25 corals per Site (Supplementary Table 3) were haphazardly selected for photosynthetic yield measurement. Divers used very dim, red-filtered diving flashlights to avoid biases during the measurements in the time interval 20:00-22:00 (i.e. in the dark). During the day, measurements were performed between 9:00–13:00 and 18:00–20:00. For each coral, the probe of the Diving PAM was pointed vertically to the coral mouth at a distance of 1 cm and the PSII fluorescence measurement was taken after holding the sensor steady until fluorescence values were steady. The diving PAM produces weak red flashes (measuring light) and detects the initial core fluorescence returns at 780–800 nm (F), upon measurement a strong white light saturates PSII reaction centers thus most of the measuring light is dissipated via fluorescence (maximal fluorescence, Fm′). From these parameters, maximal (night time) and effective quantum yield (day time) (ΔF/Fm′) were calculated using the following equation73:

Coral samples (N = 8 corals per Site) were collected in July 2019 using a hammer and a chisel by SCUBA diving and placed in label plastic containers. All samples were transported in ice to the laboratory and frozen at −20 °C until further processing. Coral tissue was removed with an airbrush connected to a reservoir of phosphate buffer saline (PBS) solution filtered through a 0.22 µm filter, and the skeleton was kept for further analysis. Separation of symbiont cells and coral tissue was performed by mechanical homogenization and centrifugation of the homogenate (5000 x g for 5 min at 4 °C) using a protocol adapted from previous studies, that accounts for the validity of the separation method from a qualitative (no tissue/debris contamination in the symbiont pellet) and quantitative (complete symbiont cell precipitation) point of view36,74,75,76. Following separation, aliquots of symbiont cell suspensions were processed for cell count and chlorophyll concentration assays. Tissue and symbiont suspensions were frozen at −20 °C, then lyophilized separately and stored at −20 °C prior to δ13C and δ15N stable isotope analyses.

Symbiont cell density in the homogenate was determined by fluorescent microscopic counts (Nikon Eclipse Ti, Japan) using a hemocytometer (BOECO, Germany) and 5 replicate (1 mm2 each) cell counts per sample (N = 4 corals per Site). Each replica was photographed both in brightfield and in fluorescent light using 440 nm emission to identify chlorophyll. Cell counting was performed using NIS-Elements Advanced Research (version 4.50.00, Nikon, Japan) with 0.5< Circularity <1, and the typical diameter parameter was to set between 5 and 15 µm. Chlorophyll-a (chl-a) concentration was measured in 2 ml of resuspended symbiont cell homogenate that was filtered onto a Whatman GF/C filter and incubated overnight with 1 ml 90% acetone at 4 °C. After incubation, the filter was manually homogenized, and the solution was filtered through a 0.22 µm syringe filter. A NanoDrop (Thermo-Fisher, United States) was used for spectrophotometric measurements at 630, 647, 664 and 697 nm wavelengths and absorbance values were used to calculate the chl-a concentration based on the following equation77:

Polyps collected for histology were immediately fixed in a formalin solution (10% formaldehyde in 37% seawater saturated with calcium carbonate) before being transferred to the lab. Samples were then postfixed in Bouin solution (composed of 15 ml picric acid saturated aqueous solution, 5 ml formaldehyde 37%, and 1 ml glacial acetic acid). After decalcification in EDTA and dehydration in increasing concentration of ethanol (from 80% to 100%), polyps were embedded in paraffin. Sections were cut at 7 μm intervals along the oral-aboral axis. Tissues were stained with Mayer's hematoxylin (Carlo Erba) and eosin (Sigma-Aldrich®)78.

Histological observations of the coral tissue were performed using a light microscope NIKON Eclipse 80i coupled with a Nikon NIS-Elements D high-resolution digital camera. Sections of the three polyps were photographed at the same distance from the oral pole at 4x magnification. Subsequently, a detail of the mesentery was photographed at 40x magnification for the observation of the symbiont cells within the endoderm.

Analyses were performed at the Godwin laboratory for Paleoclimate Research, Dept. of Earth Sciences, Cambridge University (UK). Tissue and symbiont samples (N = 4 corals per Site) were analysed for percentage carbon and nitrogen, 12C/13C and 14N/15N using a Costech Elemental Analyser attached to a Thermo DELTA V mass spectrometer in continuous flow mode. Reference standards from IAEA in Vienna were analysed along with the samples. The dried sample/standard was carefully weighed into a tin capsule, sealed and loaded into the auto-sampler. Reference standards were run at intervals throughout the sequence and these values were used to calibrate to the international standards for 14N/15N (δ15N air) and 12C/13C (δ13C VPDB). Precision of analyses is +/−0.05 % for C and N, better than 0.1 % for 12C/13C and better than 0.1 % for 14N/15N.

Frozen polyps of B. europaea (N = 4 corals per Site) were powdered with a mortar pestle using liquid nitrogen79. DNA was isolated with the Wizard Genomic DNA Purification kit (Promega) according to the manufacturer's instructions. Quality and quantity of extracted DNA was double checked using electrophoresis (0.8% agarose gel) and spectrophotometric measurements (λ=260 nm/280 nm).

The high resolution psbA non-coding region (psbAncr) from the chloroplast mini-circle genes of dinoflagellates80,81 was amplified and then directly sequenced. The ‘universal’ primers psbAFor_1 (5´GCA GCT CAT GGT TAT TTT GGT AGA C 3´) and psbARev_1 (5´AAT TCC CAT TCT CTA CCC ATC C 3´), designed to amplify the psbAncr for most Symbiodiniaceae81, were used with the following PCR conditions: 94 °C for 2 min; then 40 cycles of 94 °C 10 s, 55 °C for 30 s and 72 °C for 2 min; and a final extension at 72 °C for 10 min. The internal primers Philozoon-psbAF (5´ATT TGG TTC ACA GCG CTT GG 3´) and Philozoon-psbAR (5´CCA TTT GAC TCC CAC ACT GGA) were also used for nucleotide sequencing through the middle region of the amplified fragment82. Direct Sanger sequencing on PCR amplified DNA was performed using Big Dye 3.1 reagents (Life Sciences) and the Applied Biosystems 3730XL instrument.

Data were checked for normality using a Kolmogorov–Smirnov test (N > 50) and Shapiro-Wilk test (N < 50) and for homogeneity using Levene's Test. One-way analysis of variance (ANOVA) and the non-parametric Kruskal–Wallis equality-of-populations rank were used to assess differences in environmental parameters (number of observations for each environmental parameter are reported in Supplementary Table 1), symbiont cell density, chlorophyll-a concentration, δ13C, δ15N, C/N ratios, and amount of carbon translocated from the symbiont to the host among Sites (N = 4 corals per Site for all listed biological parameters). Where significant, pairwise comparisons between species were performed via LSD or Mann–Whitney U test. Data analyses were performed using SPSS Statistics 26.0 and GraphPad Prism 9 software. Due to the heteroskedastic dataset, mean minimum fluorescence (F), maximum fluorescence (Fm’) and effective quantum yield (ΔF/Fm′) were compared among Sites and time intervals with a permutation multivariate analysis of variance (PERMANOVA)83 based on Euclidean distances, using a crossed design with two fixed factors (factor "Site" with 3 levels: Site 1, Site 2, Site 3; factor "time interval" with 3 levels: 9:00–13:00, 18:00–20:00, 20:00–22:00) and 999 permutations (number of corals analyzed are reported in Supplementary Table 3). PERMANOVA analyses were performed with software Primer 6 (Primer-e Ltd).

Base calling on chromatograms was visually inspected for accuracy (Geneious v. 11.0.3) and the edited sequences aligned initially using the online application of ClustalW2 (http://www.ebi.ac.uk/Tools/msa/clustalw2/) (N = 4 corals per Site). Further adjustments to alignments were made upon visual inspection of the output file. Final edited sequences were deposited in GenBank. Phylogenetic analyses using Maximum Parsimony, confirmed with Maximum Likelihood, were conducted using the software PAUP (v. 4.0a136; Swofford, 2014) on aligned sequences. One-thousand bootstrap replicates were used to assess statistical significance of internal branching.

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

All source data underlying the graphs presented in the main figures are reported in Supplementary Data 1 and 2. All data and materials produced by this study are available from the corresponding author upon request.

Steffen, W. Introducing the Anthropocene: The human epoch. Ambio 50, 1784–1787 (2021).

Article Google Scholar

Keys, P. W. et al. Anthropocene risk. Nat. Sustain. 2, 667–673 (2019).

Article Google Scholar

Bell, G. Evolutionary rescue and the limits of adaptation. Philos. Trans. R. Soc. B Biol. Sci. 368, 20120080 (2013).

Article Google Scholar

Byrne, M. & Przeslawski, R. Multistressor impacts of warming and acidification of the ocean on marine invertebrates’ life histories. Integr. Comp. Biol. 53, 582–596 (2013).

Article CAS Google Scholar

Feely, R. A. et al. Impact of anthropogenic CO2 on the CaCO3 system in the oceans. Science 305, 362–366 (2004).

Article CAS Google Scholar

Doney, S. C., Fabry, V. J., Feely, R. A. & Kleypas, J. A. Ocean acidification: the other CO2 problem. Ann. Rev. Mar. Sci. 1, 169–192 (2009).

Article Google Scholar

Hill, T. S. & Hoogenboom, M. O. The indirect effects of ocean acidification on corals and coral communities. Coral Reefs https://doi.org/10.1007/s00338-022-02286-z (2022).

Biagi, E. et al. Patterns in microbiome composition differ with ocean acidification in anatomic compartments of the Mediterranean coral Astroides calycularis living at CO2 vents. Sci. Total Environ. 724, 138048 (2020).

Article CAS Google Scholar

Chen, B. et al. Microbiome community and complexity indicate environmental gradient acclimatisation and potential microbial interaction of endemic coral holobionts in the South China Sea. Sci. Total Environ. 765, 142690 (2021).

Article CAS Google Scholar

Palumbi, S. R., Barshis, D. J., Traylor-Knowles, N. & Bay, R. A. Mechanisms of reef coral resistance to future climate change. Science 344, 895–898 (2014).

Article CAS Google Scholar

Wood, R. The ecological evolution of reefs. Annu. Rev. Ecol. Syst. 29, 179–206 (1998).

Article Google Scholar

Drake, J. L. et al. How corals made rocks through the ages. Glob. Chang. Biol. 26, 31–53 (2020).

Article Google Scholar

Stanley, G. D. Photosymbiosis and the evolution of modern coral reefs. Science 312, 857–858 (2006).

Article CAS Google Scholar

Kitahara, M. V., Cairns, S. D., Stolarski, J., Blair, D. & Miller, D. J. A comprehensive phylogenetic analysis of the scleractinia (Cnidaria, Anthozoa) based on mitochondrial CO1 sequence data. PLoS One. 5, e11490 (2010).

Article Google Scholar

Dubinsky, Z. & Jokiel, P. Ratio of energy and nutrient fluxes regulates symbiosis between zooxanthellae and corals. Pac. Sci. 48, 313–324 (1994).

Google Scholar

Falkowski, P. G., Dubinsky, Z., Muscatine, L. & Porter, J. W. Light and the bioenergetics of a symbiotic coral. Bioscience 34, 705–709 (1984).

Article CAS Google Scholar

Frankowiak, K., Roniewicz, E. & Stolarski, J. Photosymbiosis in Late Triassic scleractinian corals from the Italian Dolomites. PeerJ 9, e11062 (2021).

Article Google Scholar

Davy, S. K., Allemand, D. & Weis, V. M. Cell biology of cnidarian-dinoflagellate symbiosis. Microbiol. Mol. Biol. Rev. 76, 229–261 (2012).

Article CAS Google Scholar

Kremer, P. Ingestion and elemental budgets for Linuche unguiculata, a scyphomedusa with zooxanthellae. J. Mar. Biol. Assoc. UK. 85, 613–625 (2005).

Article Google Scholar

Welsh, D. T., Dunn, R. J. K. & Meziane, T. Oxygen and nutrient dynamics of the upside down jellyfish (Cassiopea sp.) and its influence on benthic nutrient exchanges and primary production. Hydrobiologia 635, 351–362 (2009).

Article CAS 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 (1981).

Article CAS Google Scholar

Ferrier‐Pagès, C. & Leal, M. C. Stable isotopes as tracers of trophic interactions in marine mutualistic symbioses. Ecol. Evol. 9, 723–740 (2019).

Article Google Scholar

Teixidó, N. et al. Ocean acidification causes variable trait shifts in a coral species. Glob. Chang. Biol. 26, 6813–6830 (2020).

Article Google Scholar

Fantazzini, P. et al. Gains and losses of coral skeletal porosity changes with ocean acidification acclimation. Nat. Commun. 6, 7785 (2015).

Article CAS Google Scholar

Prada, F. et al. Coral micro- and macro-morphological skeletal properties in response to life-long acclimatization at CO2 vents in Papua New Guinea. Sci. Rep. 11, 19927 (2021).

Article CAS Google Scholar

Kerrison, P., Hall-Spencer, J. M., Suggett, D. J., Hepburn, L. J. & Steinke, M. Assessment of pH variability at a coastal CO2 vent for ocean acidification studies. Estuar. Coast. Shelf Sci. 94, 129–137 (2011).

Article CAS Google Scholar

Johnson, V. R., Russell, B. D., Fabricius, K. E., Brownlee, C. & Hall-Spencer, J. M. Temperate and tropical brown macroalgae thrive, despite decalcification, along natural CO2 gradients. Glob. Chang. Biol. 18, 2792–2803 (2012).

Article Google Scholar

Caroselli, E. et al. Low and variable pH decreases recruitment efficiency in populations of a temperate coral naturally present at a CO2 vent. Limnol. Oceanogr. 64, 1059–1069 (2019).

Article CAS Google Scholar

González-Delgado, S. & Hernández, J. C. The importance of natural acidified systems in the study of ocean acidification: what have we learned? Adv. Mar. Biol. 80, 57–99 (2018).

Article Google Scholar

Capaccioni, B., Tassi, F., Vaselli, O., Tedesco, D. & Poreda, R. Submarine gas burst at Panarea Island (southern Italy) on 3 November 2002: A magmatic versus hydrothermal episode. J. Geophys. Res. 112, B05201 (2007).

Google Scholar

Reggi, M. et al. Biomineralization in mediterranean corals: The role of the intraskeletal organic matrix. Cryst. Growth Des. 14, 4310–4320 (2014).

Article CAS Google Scholar

Prada, F. et al. Ocean warming and acidification synergistically increase coral mortality. Sci. Rep. 7, 1–10 (2017).

Article Google Scholar

Goffredo, S. et al. Biomineralization control related to population density under ocean acidification. Nat. Clim. Chang. 4, 593–597 (2014).

Article CAS Google Scholar

Wall, M. et al. Linking internal carbonate chemistry regulation and calcification in corals growing at a Mediterranean CO2 vent. Front. Mar. Sci. 6, 699 (2019).

Article Google Scholar

Zohary, T., Erez, J., Gophen, M., Berman-Frank, I. & Stiller, M. Seasonality of stable carbon isotopes within the pelagic food web of Lake Kinneret. Limnol. Oceanogr. 39, 1030–1043 (1994).

Article CAS Google Scholar

Xu, S. et al. Spatial variations in the trophic status of Favia palauensis corals in the South China Sea: Insights into their different adaptabilities under contrasting environmental conditions. Sci. China Earth Sci. 64, 839–852 (2021).

Article Google Scholar

Horwitz, R., Borell, E. M., Yam, R., Shemesh, A. & Fine, M. Natural high pCO2 increases autotrophy in Anemonia viridis (Anthozoa) as revealed from stable isotope (C, N) analysis. Sci. Rep. 5, 1–9 (2015).

Article Google Scholar

Chen, B., Zou, D., Zhu, M. & Yang, Y. Effects of CO2 levels and light intensities on growth and amino acid contents in red seaweed Gracilaria lemaneiformis. Aquac. Res. 48, 2683–2690 (2017).

Article CAS Google Scholar

Winters, G., Beer, S., Zvi, B., Brickner, I. & Loya, Y. Spatial and temporal photoacclimation of Stylophora pistillata: zooxanthella size, pigmentation, location and clade. Mar. Ecol. Prog. Ser. 384, 107–119 (2009).

Article Google Scholar

Fitt, W. K., McFarland, F. K., Warner, M. E. & Chilcoat, G. C. Seasonal patterns of tissue biomass and densities of symbiotic dinoflagellates in reef corals and relation to coral bleaching. Limnol. Oceanogr. 45, 677–685 (2000).

Article CAS Google Scholar

Wangpraseurt, D., Larkum, A. W. D., Ralph, P. J. & Kühl, M. Light gradients and optical microniches in coral tissues. Front. Microbiol. 3, 1–9 (2012).

Article Google Scholar

Krief, S. et al. Physiological and isotopic responses of scleractinian corals to ocean acidification. Geochim. Cosmochim. Acta. 74, 4988–5001 (2010).

Article CAS Google Scholar

Scucchia, F., Malik, A., Zaslansky, P., Putnam, H. M. & Mass, T. Combined responses of primary coral polyps and their algal endosymbionts to decreasing seawater pH. Proc. R. Soc. B Biol. Sci. 288, 20210328 (2021).

Article CAS Google Scholar

Anthony, K. R. N., Connolly, S. R. & Willis, B. L. Comparative analysis of energy allocation to tissue and skeletal growth in corals. Limnol. Oceanogr. 47, 1417–1429 (2002).

Article Google Scholar

LaJeunesse, T. C. et al. Systematic revision of Symbiodiniaceae highlights the antiquity and diversity of coral endosymbionts. Curr. Biol. 28, 2570–2580.e6 (2018).

Article CAS Google Scholar

Howells, E. J. et al. Coral thermal tolerance shaped by local adaptation of photosymbionts. Nat. Clim. Chang. 2, 116–120 (2012).

Article Google Scholar

Brading, P. et al. Differential effects of ocean acidification on growth and photosynthesis among phylotypes of Symbiodinium (Dinophyceae). Limnol. Oceanogr. 56, 927–938 (2011).

Article CAS Google Scholar

Takahashi, T., Broecker, W. S. & Langer, S. Redfield ratio based on chemical data from isopycnal surfaces. J. Geophys. Res. 90, 6907 (1985).

Article CAS Google Scholar

Xu, Z. et al. Changes of carbon to nitrogen ratio in particulate organic matter in the marine mesopelagic zone: A case from the South China Sea. Mar. Chem. 231, 103930 (2021).

Article CAS Google Scholar

Crawford, D. W. et al. Low particulate carbon to nitrogen ratios in marine surface waters of the Arctic. Glob. Biogeochem. Cycles. 29, 2021–2033 (2015).

Article CAS Google Scholar

Kikumoto, R. et al. Nitrogen isotope chemostratigraphy of the Ediacaran and Early Cambrian platform sequence at Three Gorges, South China. Gondwana Res. 25, 1057–1069 (2014).

Article CAS Google Scholar

DeNiro, M. J. & Epstein, S. Influence of diet on the distribution of carbon isotopes in animals. Geochim. Cosmochim. Acta 42, 495–506 (1978).

Article CAS Google Scholar

Benavides, M., Bednarz, V. N. & Ferrier-Pagès, C. Diazotrophs: Overlooked key players within the coral symbiosis and tropical reef ecosystems? Front. Mar. Sci. 4, 10 (2017).

Article Google Scholar

Wannicke, N., Frey, C., Law, C. S. & Voss, M. The response of the marine nitrogen cycle to ocean acidification. Glob. Chang. Biol. 24, 5031–5043 (2018).

Article Google Scholar

Bourne, D. G., Morrow, K. M. & Webster, N. S. Insights into the coral microbiome: underpinning the health and resilience of reef ecosystems. Annu. Rev. Microbiol. 70, 317–340 (2016).

Article CAS Google Scholar

Palladino, G. et al. Metagenomic shifts in mucus, tissue and skeleton of the coral Balanophyllia europaea living along a natural CO2 gradient. ISME Commun. 2, 65 (2022).

Article Google Scholar

Muscatine, L. et al. Stable isotopes (δ13C and δ15N) of organic matrix from coral skeleton. Proc. Natl Acad. Sci. 102, 1525–1530 (2005).

Article CAS Google Scholar

Lesser, M. P. et al. Nitrogen fixation by symbiotic cyanobacteria provides a source of nitrogen for the scleractinian coral Montastraea cavernosa. Mar. Ecol. Prog. Ser. 346, 143–152 (2007).

Article CAS Google Scholar

Alamaru, A., Loya, Y., Brokovich, E., Yam, R. & Shemesh, A. Carbon and nitrogen utilization in two species of Red Sea corals along a depth gradient: Insights from stable isotope analysis of total organic material and lipids. Geochim. Cosmochim. Acta. 73, 5333–5342 (2009).

Article CAS Google Scholar

Lesser, M. P., Mazel, C. H., Gorbunov, M. Y. & Falkowski, P. G. Discovery of symbiotic nitrogen-fixing cyanobacteria in corals. Science 305, 997–1000 (2004).

Article CAS Google Scholar

Lesser, M. P., Morrow, K. M., Pankey, S. M. & Noonan, S. H. C. Diazotroph diversity and nitrogen fixation in the coral Stylophora pistillata from the Great Barrier Reef. ISME J. 12, 813–824 (2018).

Article CAS Google Scholar

Marcelino, V. R., Morrow, K. M., Oppen, M. J. H., Bourne, D. G. & Verbruggen, H. Diversity and stability of coral endolithic microbial communities at a naturally high pCO2 reef. Mol. Ecol. 26, 5344–5357 (2017).

Article CAS Google Scholar

Rädecker, N., Pogoreutz, C., Voolstra, C. R., Wiedenmann, J. & Wild, C. Nitrogen cycling in corals: the key to understanding holobiont functioning? Trends Microbiol. 23, 490–497 (2015).

Article Google Scholar

Santos, H. F. et al. Climate change affects key nitrogen-fixing bacterial populations on coral reefs. ISME J. 8, 2272–2279 (2014).

Article Google Scholar

Olson, N. D., Ainsworth, T. D., Gates, R. D. & Takabayashi, M. Diazotrophic bacteria associated with Hawaiian Montipora corals: Diversity and abundance in correlation with symbiotic dinoflagellates. J. Exp. Mar. Bio. Ecol. 371, 140–146 (2009).

Article CAS Google Scholar

Zheng, X. et al. Effects of ocean acidification on carbon and nitrogen fixation in the hermatypic coral Galaxea fascicularis. Front. Mar. Sci. 8, 644965 (2021).

Article Google Scholar

Lewis, E. & Wallace, D. Program developed for CO2 system calculations. Ornl/Cdiac-105 1–21 (1998).

Dickson, A. G. & Millero, F. J. A comparison of the equilibrium constants for the dissociation of carbonic acid in seawater media. Deep Sea Res. Part A. Oceanogr. Res. Pap. 34, 1733–1743 (1987).

Article CAS Google Scholar

Dickson, A. G. Thermodynamics of the dissociation of boric acid in potassium chloride solutions from 273.15 to 318.15 K. J. Chem. Eng. Data. 35, 253–257 (1990).

Article CAS Google Scholar

Mehrbach, C., Culberson, C. H., Hawley, J. E. & Pytkowicx, R. M. Measurement of the apparent dissociation constants of carbonic acid in seawater at atmospheric pressure. Limnol. Oceanogr. 18, 897–907 (1973).

Article CAS Google Scholar

Ivancic, I. & Degobbis, D. An optimal manual procedure for ammonia analysis in natural waters by the indophenol blue method. Water Res. 18, 1143–1147 (1984).

Article CAS Google Scholar

Parson, T. R., Maita, Y. & Llli, C. M. A manual of chemical & biological methods for seawater analysis. (Elsevier, 1984). https://doi.org/10.1016/C2009-0-07774-5

Schreiber, U. Pulse-Amplitude-Modulation (PAM) fluorometry and saturation pulse method: an overview. in Chlorophyll a Fluorescence 1367, 279–319 (Springer Netherlands, 2004).

Grover, R., Maguer, J. F., Reynaud-Vaganay, S. & Ferrier-Pagès, C. Uptake of ammonium by the scleractinian coral Stylophora pistillata: Effect of feeding, light, and ammonium concentrations. Limnol. Oceanogr. 47, 782–790 (2002).

Article Google Scholar

Tremblay, P., Grover, R., Maguer, J. F., Hoogenboom, M. & Ferrier-Pagès, C. Carbon translocation from symbiont to host depends on irradiance and food availability in the tropical coral Stylophora pistillata. Coral Reefs. 33, 1–13 (2014).

Article Google Scholar

Pupier, C. A. et al. Productivity and carbon fluxes depend on species and symbiont density in soft coral symbioses. Sci. Rep. 9, 17819 (2019).

Article Google Scholar

Ritchie, R. J. Universal chlorophyll equations for estimating chlorophylls a, b, c, and d and total chlorophylls in natural assemblages of photosynthetic organisms using acetone, methanol, or ethanol solvents. Photosynthetica 46, 115–126 (2008).

Article CAS Google Scholar

Goffredo, S., Arnone, S. & Zaccanti, F. Sexual reproduction in the Mediterranean solitary coral Balanophyllia europaea (Scleractinia, Dendrophylliidae). Mar. Ecol. Prog. Ser. 229, 83–94 (2002).

Article Google Scholar

Barshis, D. J. et al. Genomic basis for coral resilience to climate change. Proc. Natl Acad. Sci. 110, 1387–1392 (2013).

Article CAS Google Scholar

Moore, R. B. Highly organized structure in the non-coding region of the psbA minicircle from clade C Symbiodinium. Int. J. Syst. Evol. Microbiol. 53, 1725–1734 (2003).

Article CAS Google Scholar

LaJeunesse, T. C. & Thornhill, D. J. Improved resolution of reef-coral endosymbiont (Symbiodinium) species diversity, ecology, and evolution through psbA non-coding region genotyping. PLoS One. 6, e29013 (2011).

Article CAS Google Scholar

LaJeunesse, T. C. et al. Revival of Philozoon Geddes for host-specialized dinoflagellates, ‘zooxanthellae’, in animals from coastal temperate zones of northern and southern hemispheres. Eur. J. Phycol. 57, 166–180 (2022).

Article Google Scholar

Anderson, M. J. PERMANOVA: a FORTRAN computer program for permutational multivariate analysis of variance. Wiley StatsRef: Statistics Reference Online (2005).

Download references

The research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Programme (FP7/2007–2013)/ERC grant agreement n° 249930—CoralWarm: Corals and global warming: the Mediterranean versus the Red Sea. Research project implemented under the National Recovery and Resilience Plan (NRRP), Mission 4 Component 2 Investment 1.4 - Call for tender No. 3138 of 16 December 2021, rectified by Decree n.3175 of 18 December 2021 of Italian Ministry of University and Research funded by the European Union – NextGenerationEU. Project code CN_00000033, Concession Decree No. 1034 of 17 June 2022 adopted by the Italian Ministry of University and Research, CUP J33C22001190001, Project title "National Biodiversity Future Center - NBFC". LaJeunesse was supported by funding from the USA National Science Foundation (grant OCE-1636022). Bartolo Basile, Francesco Sesso and Eolo Sub diving center assisted in the field. The Scientific Diving School collaborated with the underwater activities. We are grateful to Francesco Sesso for the image of B. europaea. We are also thankful to Prof. Costantino Vetriani and Dr. Corday Selden from Rutgers University whose comments were instrumental in developing certain parts of the Discussion.

Fiorella Prada

Present address: Environmental Biophysics and Molecular Ecology Program, Department of Marine and Coastal Sciences, Rutgers University, New Brunswick, NJ, 08901, USA

These authors contributed equally: Fiorella Prada, Silvia Franzellitti.

Marine Science Group, Department of Biological, Geological and Environmental Sciences, University of Bologna, Via F. Selmi 3, 40126, Bologna, Italy

Fiorella Prada, Erik Caroselli, Arianna Mancuso, Chiara Marchini & Stefano Goffredo

Fano Marine Center, The Inter-Institute Center for Research on Marine Biodiversity, Resources and Biotechnologies, Viale Adriatico 1/N, 61032, Fano, Italy

Fiorella Prada, Silvia Franzellitti, Erik Caroselli, Mauro Marini, Arianna Mancuso, Chiara Marchini, Marco Candela, Giuseppe Falini & Stefano Goffredo

Animal and Environmental Physiology Laboratory, Department of Biological, Geological and Environmental Sciences, University of Bologna, via S. Alberto 163, 48123, Ravenna, Italy

Silvia Franzellitti, Lorenzo Sana & Alessia Puglisi

The Interuniversity Institute for Marine Sciences in Eilat, PO Box 469, Eilat, 88103, Israel

Itay Cohen

Institute of Biological Resources and Marine Biotechnology, National Research Council (CNR), Largo Fiera della Pesca 2, 60125, Ancona, Italy

Mauro Marini & Alessandra Campanelli

Unit of Microbiome Science and Biotechnology, Department of Pharmacy and Biotechnology, University of Bologna, 40126, Bologna, Italy

Marco Candela

Department of Marine Biology, The Leon H. Charney School of Marine Sciences, University of Haifa, Haifa, Israel

Tali Mass

Department of Earth Sciences, University of Florence, via la Pira 4, Firenze, Italy

Franco Tassi

Institute of Geosciences and Earth Resources (IGG), National Research Council of Italy (CNR), via la Pira 4, Firenze, Italy

Franco Tassi

Department of Biology, The Pennsylvania State University, 208 Mueller Laboratory, University Park, PA, 16802, USA

Todd C. LaJeunesse

The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat-Gan, 52900, Israel

Zvy Dubinsky

Department of Chemistry "Giacomo Ciamician", University of Bologna, 40126, Bologna, Italy

Giuseppe Falini

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

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

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

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

S.G., G.F., and Z.D. conceived and designed the research. F.P., E.C., A.M., and S.G. collected the samples and performed the diving field work. F.P., S.F., E.C., M.M., A.C., L.S., C.M., A.P., T.M., F.T., T.C.L. performed the laboratory analyses. F.P., E.C., S.F., I.C., M.C., T.C.L. and S.G. analysed the data. F.P., S.F., and E.C. wrote the first draft. All authors contributed to writing the manuscript and participated in the scientific discussion.

Correspondence to Erik Caroselli, Todd C. LaJeunesse or Stefano Goffredo.

The authors declare no competing interests.

Communications Biology thanks Chris Langdon and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editors: Linn Hoffmann and Luke R. Grinham. Peer reviewer reports are available.

Publisher's note 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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Prada, F., Franzellitti, S., Caroselli, E. et al. Acclimatization of a coral-dinoflagellate mutualism at a CO2 vent. Commun Biol 6, 66 (2023). https://doi.org/10.1038/s42003-022-04327-3

Download citation

Received: 29 May 2022

Accepted: 01 December 2022

Published: 18 January 2023

DOI: https://doi.org/10.1038/s42003-022-04327-3

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.