# Coral Skeletal Cores as Windows Into Past Symbiodiniaceae Community Dynamics

Jose F Grillo, Vanessa Tirpitz, Jessica Reichert, Marine Canesi, Stéphanie Reynaud, Eric Douville, Maren Ziegler

► To cite this version:

Jose F Grillo, Vanessa Tirpitz, Jessica Reichert, Marine Canesi, Stéphanie Reynaud, et al., Coral Skeletal Cores as Windows Into Past Symbiodiniaceae Community Dynamics. Global Change Biology, 2025, 31 (11), 10.1111/gcb.70575. hal-05468862

HAL Id: hal-05468862

https://hal.science/hal-05468862v1

Submitted on 21 Jan 2026

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

RESEARCH ARTICLE

OPEN ACCESS

# Coral Skeletal Cores as Windows Into Past Symbiodiniaceae Community Dynamics

Jose F. Grillo<fnref n="1" /> ID
Eric Douville³ ID

Vanessa Tirpitz<fnref n="1" /> ID
Maren Ziegler<fnref n="1" /> ID

Jessica Reichert<fnref n="1" />,² ID
Marine Canesi³ ID

Stéphanie Reynaud⁴ ID

Marine Canesi³ ID
Stéphanie Reynaud⁴ ID

**Correspondence:** Maren Ziegler (maren.ziegler@bio.uni-giessen.de)

**Received:** 5 June 2025

**Revised:** 24 September 2025

**Accepted:** 26 September 2025

**Funding:** MZ acknowledges financial support from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project number 469364832 (M.Z.) – SPP 2299/Project number 441832482. The work was made possible with the support of a scholarship from the German Academic Exchange Service (DAAD) to JFG. X-ray analysis and growth modelling were carried out and funded by the Agence Nationale de la Recherche under the grant number ANR-22-CE02-0025 (project COR-Resilience).

**Keywords:** ancient DNA | climate change | community reconstruction | coral bleaching | ITS2-amplicon sequencing | shifting | shuffling

## ABSTRACT

The symbiosis between the dinoflagellate Symbiodiniaceae family and reef-building corals underpins the productivity of coral reefs. This relationship facilitates the deposition of calcium-carbonate skeletons that build the reef structure thanks to the energy derived from photosynthesis. The loss of Symbiodiniaceae from coral tissues—resulting in coral bleaching—impedes coral growth and can lead to mass mortality if the symbiosis fails to recover. Given that Symbiodiniaceae communities are dynamic and can shift in response to environmental stressors in the decades to centuries-long lifespan of coral colonies, understanding these changes is crucial. Although the reconstruction of Symbiodiniaceae communities from coral skeleton records has recently been demonstrated as feasible, no studies have yet assessed reconstructions across different species and locations. Here, we present an approach to use coral skeletons for reconstructing the Symbiodiniaceae community on decadal and centennial scales and for resolving dynamics related to coral species and the environmental history of sampling locations. For this, we used dated coral skeleton cores from *Porites lobata* and *Diploastrea heliopora*, species commonly used as climate archives, sampled in Palau and Papua New Guinea. We also examined the effect of various DNA extraction protocols on community reconstruction. Here we show that the reconstructed Symbiodiniaceae communities significantly varied across all cores and DNA extraction methods, with decalcification-based protocols enhancing the retrieval of skeletal-bound DNA. Moreover, we observed distinct community dynamics related to the specific coral host and sampling location. Notably, associations of Symbiodiniaceae dynamics with past heat stress events were apparent in cores of both species from Palau. Our findings enable a deeper understanding of the temporal and spatial variability in Symbiodiniaceae communities, offering insights that may refine the use of paleobiological proxies in climate studies and reveal broader ecological trends and microbially aided adaptation pathways in corals.

## 1 | Introduction

The health and productivity of coral reefs builds on the symbiosis of reef-building corals with photosynthetic dinoflagellates from

the diverse Symbiodiniaceae family that reside within their tissues (Nitschke et al. 2022; Helgoe et al. 2024). The translocation of carbohydrates produced by Symbiodiniaceae during photosynthesis is used by the coral to deposit their characteristic calcium

<footnote n="1">¹Marine Holobiomics Lab, Department of Animal Ecology and Systematics, Justus Liebig University Giessen, Giessen, Germany | ²Hawai'i Institute of Marine Biology, University of Hawai'i at Mānoa, Honolulu, Hawaii, USA | ³Laboratoire des Sciences du Climat et de l'Environnement, LSCE/IPSL, CEA-CNRS-UVSQ, Université Paris-Saclay, Gif-sur-Yvette, France | ⁴Centre Scientifique de Monaco, Principality of Monaco, Monaco</footnote>

<footnote>This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.</footnote>
carbonate (CaCO₃) skeleton, which builds the reef framework (Al-Horani 2015). At the same time, this symbiosis between corals and Symbiodiniaceae is also their Achilles' heel due to its susceptibility to heat stress (Helgoe et al. 2024). Increasing duration, frequency, and severity of marine heat waves cause coral bleaching, which is the loss of the Symbiodiniaceae from coral tissue (Frölicher et al. 2018; Helgoe et al. 2024). Coral bleaching results in starvation, reduced skeletal growth, and may lead to mortality (Helgoe et al. 2024). Marine heatwaves aggravated by climate change have caused mass bleaching events, which led to the loss of between 30% and 90% coral cover in reef ecosystems across the Pacific Ocean and the Arabian Seas (Riegl et al. 2018; Raymundo et al. 2019). In 2024, the Fourth Global Coral Bleaching Event was recorded, and while its impact is still being studied, future global-scale bleaching events are expected to continue with rising CO₂ emissions (Reimer et al. 2024).

Symbiodiniaceae are a diverse family, whose members are functionally distinct regarding environmental niche space and the traits they confer to the coral (LaJeunesse et al. 2018). These functions are assumed to be determined by one or a few abundant symbiont taxa that dominate the community along with several background symbionts that usually persist at low abundance (Quigley et al. 2014; Davies et al. 2023). Symbiodiniaceae species differ in the amount of carbon they translocate to the coral, influencing their host's growth rates (Pettay et al. 2015), and they notoriously differ in heat stress tolerance (Swain et al. 2017; Xiao et al. 2022). So-called adaptive bleaching (Buddemeier and Fautin 1993), a dynamic process in which the main symbiont taxon is replaced, may rescue the coral after bleaching and may subsequently result in higher heat stress tolerance (Baker et al. 2004; Boulotte et al. 2016). The bleached symbionts may be replaced from the extant community through shuffling or recruited from the environment through switching (Baker 2003; Quigley et al. 2022). As such, the study of Symbiodiniaceae through tissue or mucus sampling provides a snapshot of the community at particular time points that are limited by the availability of samples (Million et al. 2024). Extraction of historical DNA from dried museum specimens can extend surveys back in time (Rowe et al. 2011). A study conducted on up to 172-year-old soft corals in museum collections suggested that Symbiodiniaceae communities remained largely unchanged compared to recent samples (Baker et al. 2013). Yet, approaches that rely on tissue samples fail to reveal long-term dynamics following periods of environmental stress throughout the century-scale lifespan of reef-building corals.

Coral skeletons have been extensively used for the reconstruction of past environmental conditions through analyses of stable isotopes, trace elements, and biogeochemistry of the skeletal material (Thompson 2022; Canesi et al. 2024). The signals recorded in coral skeletons over centuries are proxies for past climates, ocean chemistry, and pollution (Ramos et al. 2009; Thompson 2022; Canesi et al. 2024). Recently, ancient DNA (aDNA) (DNA isolated from samples over 100 years old) applications probed sub-fossil deep-sea coral skeletons for species identification (Waller et al. 2007) and subsequently to identify genomic similarities between millennia-old and living colonies of *Acropora palmata* and their bacterial communities (Scott et al. 2022). The skeletons of tropical stony corals have successfully been used for reconstructions of host phylogenies and of prokaryotic and Symbiodiniaceae communities (Scott et al. 2022;

González-Pech et al. 2024; Rey et al. 2025). Nevertheless, the use of single cores is associated with inherent limitations that have prevented the study of Symbiodiniaceae community changes related to geographic and coral host-associated sources of variability that remain to be explored in future studies.

The overarching goal of this study was to evaluate the use of skeleton cores from massive coral species for the reconstruction of Symbiodiniaceae communities at decadal and centennial time scales across species and locations. For this, century-old coral skeleton cores were sampled from two massive species (*Porites lobata* and *Diploastrea heliopora*) from two locations in the Pacific (Papua New Guinea and Palau), and ancient Symbiodiniaceae community dynamics were analyzed using ITS2 amplicon sequencing. Our specific aims were to: (1) assess the effect of three different DNA extraction methods on the reconstructed Symbiodiniaceae communities, (2) evaluate changes in reconstructed communities obtained by species and sampling location, and (3) probe differences in the reconstruction dynamics between two massive-growing coral species in the context of past potential heat stress events. The results generated by this study establish the potential for long-term Symbiodiniaceae community reconstructions from multiple massive reef-building coral species across regions, filling a crucial information gap in current survivorship projections and coral restoration efforts.

## 2 | Materials and Methods

### 2.1 | Collection of Coral Cores

Four coral cores were drilled from colonies located in North Palau (7°47′N, 134°35′E) and Losuia Tabungora Island in Papua New Guinea (9°21′S, 152°02′E) during the Tara Pacific expedition (Planes et al. 2019) in 2017 and 2018, respectively (Figure 1). The cores were collected as described in Canesi et al. (2023, 2024) for geochemical reconstruction of past environmental conditions and stored at room temperature after sampling. Briefly, a 40–150 cm long skeletal core was taken with a hydraulic drill (Stanley Black & Decker, US) with a 7 cm diameter corer from one living colony of *Porites lobata* (Dana, 1846) and *Diploastrea heliopora* (Lamark, 1816) per site at 13–16 m. The cores were rinsed with fresh water at the time of collection, left to air dry, and stored at room temperature until processing. *Porites lobata* is one of the species most commonly used for reconstruction of climate pasts and has a relatively high linear extension rate of 7.5–16.3 mm per year (Klein and Loya 1991; Edinger et al. 2000). *Diploastrea heliopora* is a comparably slow-growing species (ca. 2–4 mm per year) which enables long-term time series reconstruction (Masuda et al. 2010; Canesi et al. 2024).

### 2.2 | Experimental Design

To test the effect of different DNA extraction methods on the reconstruction of Symbiodiniaceae communities from coral skeleton cores, three methods were compared in this study based on commercially available kits and the presence of decalcification steps prior to the DNA isolation (Figure 1).
## Sampling

<figure><img src="image_4.png" /></figure>

## Core processing

<figure><img src="image_5.png" /></figure>

**FIGURE 1** | Overview of study design. Coral cores of the species *P. lobata* and *D. heliopora* were sampled in Palau and Papua New Guinea. Fragments along the growth axis were drilled to (I) compare DNA extraction methods, (II) compare symbiont communities across species and sampling locations, and (III) link community dynamics to bleaching alerts of the past century. Map lines delineate study areas and do not necessarily depict accepted national boundaries.

Moreover, two time frames were chosen to reconstruct the Symbiodiniaceae community dynamics. The first 30 years before sampling were analyzed to study the recent dynamics of Symbiodiniaceae communities after potential bleaching events determined by combining sea surface temperature records and bleaching alerts at each of the sampling locations by the NOAA Coral Reef Watch dating back to 1985 (Skirving et al. 2020). Moreover, 80–90 year-old fragments of *P. lobata* and 100–110 year-old fragments of *D. heliopora* were sampled to assess the effectiveness of century-scale Symbiodiniaceae community reconstructions.

## 2.3 | Coral Core Sample Processing

Skeletal cores were dated through X-Ray along the longitudinal axis, and the fragment age, corresponding to the depth of the

core, was calculated based on the linear extension rate, where 1 year of growth was assigned as the combination of a high and low-density band. The original 7-cm diameter cylindrical core, including the tissue layer, was initially halved for dating and later sectioned into quarters, with one quarter being used in this study. The quarter of the core was further divided into 5-year fragments along the longitudinal axis, with the exception of the tissue layer fragment, which corresponds to 2 years of growth (Figure 1 and Figure S1). All cores were treated with UV light for 15 min for surface decontamination prior to sample processing. The coring equipment (drill bit, ceramic mortar, and pestle) and all working surfaces were sterilized. Approximately 1 g of skeleton was taken from the midpoint along the longitudinal axis of each 5-year fragment with a sterile metallic drill bit (diameter: 6 mm) and pulverized with a sterile ceramic mortar and pestle. Ten samples per coral core were used, resulting in a total of 40 skeletal samples from the two coral species from two locations.
## 2.4 | Molecular Sample Processing

A comparison of DNA extraction methods was performed with the 10 powder samples of the *P. lobata* core from Papua New Guinea to test the effect of the extraction method on the reconstruction of Symbiodiniaceae communities in ancient coral samples. These tests included an extraction blank for each processing batch. Three different types of extraction were compared: (1) DNA was extracted from a 250 mg stock sample using the PowerSoil Pro kit (Qiagen, Germany) following the manufacturer's instructions with a final elution volume of 50 µL. Methods 2 and 3 shared a decalcification step, which was not part of method 1. (2) A 250 mg stock sample was taken, first decalcified, and then the DNA was extracted with the same PowerSoil Pro kit (Qiagen, Germany). For decalcification, sterile-filtered EDTA solution 0.5 M (pH 7.5) was added to the skeletal powder and incubated at 37°C on a Thermomixer comfort (Eppendorf, Germany) at 500 RPM until all skeletal material was dissolved. (3) The ancient DNA (aDNA) extraction protocol proposed by Rohland et al. 2018) based on a decalcification step, the use of commercial silica spin columns and custom buffers that allow the recovery of short DNA fragments (as small as 35 bp) was tested with a 50 mg sample and the binding buffer D. Briefly, the coral skeleton powder is decalcified and lysed with a buffer containing EDTA 0.5 M (pH 8.0), 10 mg/mL proteinase K and Tween 20. The lysate supernatant is added to the binding buffer D (5 M guanidine hydrochloride, 40% (vol/vol) 2-propanol, 0.12 M sodium acetate and 0.05% (vol/vol) Tween 20) and transferred to the silica spin columns. The spin columns are washed with PE buffer (Qiagen, Germany) and the DNA is eluted with a custom-made buffer based on 1 M Tris-HCl (pH 8.0), EDTA 0.5 M (pH 8.0) and Tween 20.

The samples from the remaining three cores were extracted using only the decalcification and Powersoil Pro approach following method 2 above to examine differences in the Symbiodiniaceae community dynamics between coral species and sampling locations. For the PCRs, 6.4 µL of DNA template was used with the Symbiodiniaceae-specific primers SYM_VAR_5.8S2 and SYM_VAR_REV (Hume et al. 2018) and the Multiplex PCR Master Mix (Qiagen, Germany) with a final volume of 40 µL. One negative control per PCR batch was included during the sample processing. PCRs were cleaned using the QIAquick PCR purification kit (Qiagen, Germany). The cleaned PCR products were sent for library preparation and sequencing on the DNBseq G400 sequencing platform (BGI, Hong Kong).

## 2.5 | Sequence Data Processing

Demultiplexed, quality-controlled reads were submitted to the SymPortal analytical workflow (Hume et al. 2019). SymPortal was designed for the analysis of Symbiodiniaceae communities through the recognition of co-occurring intragenomic sequences within each sample as so-called type profiles (Hume et al. 2019). Given the expected DNA degradation in the ancient DNA samples used here, we did not assume that all sequences necessary to distinguish SymPortal type profiles would be reliably recovered in the samples. In addition, the unknown mechanism of DNA or symbiont entrapment in the coral skeleton violates the SymPortal assumption of one Symbiodiniaceae species per

genus in a sample. Therefore, the resulting sequences from the SymPortal minimum entropy decomposition algorithm (post-MED sequences) were used for all analyses.

## 2.6 | Data Analyses

All analyses, visualizations, and statistical tests were performed in the R statistical software (version 4.3.2, R Core Team 2023). The absolute abundance count table of the post-MED sequences generated by the SymPortal workflow was pre-processed by removing all Symbiodiniaceae sequences present in the negative control samples for the DNA extraction and PCR negative control corresponding to each processing batch. Because of the known susceptibility of the aDNA approach to contamination (Yang and Watt 2005; Hofreiter et al. 2021), we chose this stringent approach. This resulted in an average of 39,298 post-MED sequences per sample (range: 16,940–110,234; Table S1). The relative abundance of post-MED sequences per sample was then calculated, and only sequences with a relative abundance > 0.1% over all samples were retained in the analysis. The relative abundance of Symbiodiniaceae sequences per sample was visualized with a stacked bar plot using tidyverse (version 2.0.0, Wickham et al. 2019) and the ggnested package for color palette generation (version 0.1.1). Bray-Curtis distances of the Symbiodiniaceae communities of each sample and DNA extraction (log-transformed) were visualized with an ordination plot based on non-metric Multidimensional Scaling (nMDS) with *metaMDS* from the vegan package (version 2.6-8, Oksanen et al. 2024). Finally, the number of shared Symbiodiniaceae sequences was visualized with an UpSet plot with the ComplexUpset package (version 1.3.3, Krassowski et al. 2022). Separate Permutational Multivariate Analysis of Variance (PERMANOVA, Anderson 2017) based on Bray-Curtis distances was performed to test for differences between the ITS2 communities obtained with different DNA extraction methods for the *P. lobata* core from Papua New Guinea and for differences between coral species and sampling locations in a fully crossed design, respectively. Both statistical tests were computed with 999 permutations, and multivariate homogeneity of variance was tested with the vegan package. Post hoc permutational pairwise comparisons were computed with the pairwise.adonis package (version 0.4.1, Martinez Arbizu 2017). Distance-decay plots were used to test the effect of fragment age on the similarity of the Symbiodiniaceae communities reconstructed from coral skeletons at decadal scales. Euclidean distances between fragments as a proxy for time and Bray-Curtis Symbiodiniaceae community similarity were fitted in a linear regression model. Mantel tests were performed to test for significant correlations between distance matrices based on 999 permutations.

The daily resolution sea surface temperature (SST) and bleaching alert data for both sampling locations obtained from the NOAA CRW (accessed on 18-02-2025, Skirving et al. 2020) were filtered to coincide with the 5-year intervals of the coral skeleton core and visualized as a bar plot. This metric distinguishes three bleaching alert levels, which correspond to the risk of reef-wide bleaching (level 1), the risk of reef-wide bleaching with mortality of heat-sensitive corals (level 2), and the risk of multi-species mortality (level 3).
<figure><img src="image_6.png" /><figcaption>**FIGURE 2** | Symbiodiniaceae community from one *P. lobata* skeleton core collected from Papua New Guinea reconstructed with three different DNA extraction methods. (a) Relative abundance of post-MED Symbiodiniaceae sequences per genus. The shades of color represent unique post-MED sequences. (b) Number of shared post-MED sequences between the DNA extraction methods. The intersection size (vertical bar plot) represents the number of shared sequences of a set of DNA methods represented by filled points. (c–e) Ordination of Symbiodiniaceae community composition with non-metric Multidimensional Scaling based on Bray-Curtis distances of post-MED sequences. Ellipses represent 95% CI.</figcaption></figure>

## 3 | Results

### 3.1 | Symbiodiniaceae Communities Differed According to DNA Extraction Protocol

The reconstruction of Symbiodiniaceae communities from coral skeleton cores was successful for the three DNA extraction protocols at decadal and centennial time scales, although differences were observed between the resulting communities (PERMANOVA, Pseudo-F = 1.7595, *p* < 0.001, Figure 2a,c–e; Table 1; Table S2). Of note, the two protocols including the decalcification step were more similar to each other than to the protocol lacking this step (Table S2). While they had similar multivariate homogeneity of variances

<table><caption>**TABLE 1** | Statistical test of differences in the reconstructed Symbiodiniaceae communities obtained with three DNA extraction protocols from a *P. lobata* core sampled in Papua New Guinea. PERMANOVA on Bray-Curtis distances with 999 permutations.</caption><thead><tr><th>Factor</th><th>DF</th><th>SS</th><th>R<sup>2</sup></th><th>Pseudo-F</th><th>p</th></tr></thead><tbody><tr><td>DNA extraction method</td><td>2</td><td>1.5086</td><td>0.1153</td><td>1.7597</td><td>0.0001</td></tr><tr><td>Residual</td><td>27</td><td>11.5739</td><td>0.8847</td><td></td><td></td></tr><tr><td>Total</td><td>29</td><td>13.0825</td><td>1.0000</td><td></td><td></td></tr></tbody></table>

Abbreviations: DF, degrees of freedom; Pseudo-F, pseudo-F statistic; SS, sum of squares.
<figure><img src="image_7.png" /><figcaption>**FIGURE 3** | Symbiodiniaceae community composition based on non-metric Multidimensional Scaling of Bray-Curtis distances of post-MED sequences. Ellipses represent 95% CI. Symbiodiniaceae communities (a) from both coral species and sampling locations, (b) from *P. lobata* and *D. heliopora* sampled in Palau, (c) in Papua New Guinea, and (d) from *P. lobata*, and (e) *D. heliopora* at both locations. Samples Anc1-3 are 80–90 years old for *P. lobata* and 100–110 years old for *D. heliopora*.</figcaption></figure>

(Betadisper, *p* > 0.05, Table S3), the overall patterns were also reflected in the high number of 421 to 435 unique post-MED sequences per DNA extraction method (Figure 2b). The protocols with a decalcification step shared a higher number of sequences (67) than each of them with the third protocol (51 and 63) or all protocols (37) (Figure 2b).

The Symbiodiniaceae genera *Symbiodinium*, *Breviolum*, *Cladocopium*, and *Durusdinium* were recovered from all coral core samples and all DNA extraction methods, although community composition differed between protocols (Figure 2a,c-e). Post-MED sequences in the genus *Symbiodinium* were most diverse and they also represented the most abundant fraction in all skeletal samples, but not in the top skeletal sample that included the tissue layer, where *Cladocopium* tended to be more abundant for the protocols employing a decalcification step (Figure 2a).

## 3.2 | Reconstructed Symbiodiniaceae Dynamics Differed by Species and Location

Symbiodiniaceae communities differed between coral species (PERMANOVA, Pseudo-*F*=1.63, *p*<0.001), sampling locations (PERMANOVA, Pseudo-*F*=1.69, *p*<0.001), and between all cores

across species and locations (PERMANOVA, Pseudo-*F*=2.2304, *p*<0.001, Figure 3; Tables S4 and S5). Symbiodiniaceae communities between the *P. lobata* cores differed in multivariate variance (Betadisper, *p*<0.05, Table S6), which was driven by the similarity of the core top fragments with tissue, while the older skeletal samples had distinct communities between cores (Figure 3d). No significant correlations were found among Symbiodiniaceae community similarity and fragment age in either species and sampling location (Figure S2a-d and Table S7).

Reconstructed Symbiodiniaceae communities from Papua New Guinea were more stable over decadal and centennial timeframes (Figure 4a). With the exception of the most recent sample from the top of the cores, no shifts in dominance between Symbiodiniaceae genera were observed, though sequences within genera varied. In contrast, shifts in the dominant symbiont genus were recorded in the cores from Palau (Figure 4a). While the skeletons of both coral species were largely dominated by sequences belonging to the genus *Symbiodinium*, the 10-year old sample from *P. lobata* was dominated by *Cladocopium* sequences. The 30-year old sample of *D. heliopora* was dominated by *Durusdinium* sequences, and the sequences recovered from the 25-year old sample indicate a trend where *Durusdinium* dominance slowly reversed back to *Symbiodinium* dominance that had reestablished in the 20-year
<figure><img src="image_8.png" /><figcaption>**FIGURE 4** | Symbiodiniaceae communities recovered from coral skeleton cores in relation to historical bleaching alerts. (a) Relative abundance of Symbiodiniaceae genera of *P. lobata* and *D. heliopora* sampled in Palau and Papua New Guinea. Color shades represent unique post-MED sequences. (b) Daily bleaching alert reconstructions from sea surface water temperatures for the time intervals corresponding with the coral core fragment ages. Samples Anc1-3 are 80–90 years old for *P. lobata* and 100–110 years for *D. heliopora*.</figcaption></figure>

old sample (Figure 4a). These site-specific patterns were also reflected in a higher number of shared post-MED sequences for the cores from Papua New Guinea (80) when compared to Palau (23) (Figure S3).

## 3.3 | Symbiodiniaceae Dynamics Partially Coincided with Historical Heat Stress Events

The historical SST data obtained from the NOAA Coral Reef Watch revealed an increasing frequency and severity of SST

anomalies associated with bleaching alerts in both locations over time, which peaked at over 220 and 150 bleaching-alert days for the sampling year of the cores for Palau and Papua New Guinea, respectively (Figure 4b and Figure S4). In general, a higher number of annual bleaching-alert days was recorded at Palau than at Papua New Guinea. The occurrence of a severe alert level 3 bleaching episode in Palau coincided with a change in dominant symbiont in *P. lobata* (Figure 4), while the changes observed in *D. heliopora* could not be matched to a distinct event based on the NOAA bleaching alert data.
## 4.1 | aDNA Extraction Protocol and Sample Decalcification Affect Symbiodiniaceae Community Reconstructions

Reconstruction of Symbiodiniaceae communities was successful at decadal and centennial scales from *P. lobata* skeleton cores with all DNA extraction methods. The main differences between the tested protocols included the use of a decalcification step with EDTA prior to the extraction in two out of three methods and the use of a commercial extraction kit versus a bespoke protocol for aDNA. The protocols that included the decalcification step shared the dominance of *Cladocopium* sequences in the core top sample as well as the higher number of shared post-MED sequences. These results suggest that the dissolution of the aragonite (CaCO₃) matrix yielded more comparable results regardless of the subsequent extraction protocol. The use of a decalcification step is common for aDNA extractions of coral skeletons, other marine calcifying invertebrates including mollusk shells and foraminifera, and mammal bones to release DNA bound to the skeletal architecture (Waller et al. 2007; Lendvay et al. 2020; Scott et al. 2022; González-Pech et al. 2024; Martin-Roy et al. 2024). Nevertheless, aDNA extractions employing coral skeletons have been successfully performed without the addition of a decalcification step, which is also supported by our data (Connelly et al. 2024; Rey et al. 2025). Another difference between protocols was the lower starting material in the aDNA approach of 50 mg versus 250 mg in the other protocols. However, these subsamples were all taken from the same powdered stock material, minimizing batch effects. In addition, despite differences in input mass, the two protocols that included a decalcification step yielded more similar Symbiodiniaceae communities than those of the two protocols with the same input mass. Together, these data highlight that a fraction of Symbiodiniaceae DNA is bound within the aragonite skeleton and that a decalcification step before DNA extraction enables the capture of this diversity regardless of which protocol is subsequently used for extraction.

Commercially available DNA extraction kits have previously been used for aDNA extraction of invertebrate aragonite skeletons. The PowerSoil Pro kit was selected for testing due to the inclusion of guanidine hydrochloride and silica spin columns that prevent the co-purification of PCR inhibitors common in environmental samples (Martin et al. 2021; Connelly et al. 2024; Martin-Roy et al. 2024). While lower DNA yields were reported with this kit when compared to the Blood and Tissue DNA kit (Qiagen, Germany) (Connelly et al. 2024), the combination with the decalcification step provided a higher amplification success in our early testing. This suggests that commercial kits may provide a simple alternative to the use of specialized aDNA extraction protocols with user-made buffers without a loss of reconstructed Symbiodiniaceae communities at century scales. This result is supported by the successful reconstruction of communities from 170-year-old coral skeleton core fragments employing the QIAamp DNA Micro Kit (Qiagen, Germany) (Rey et al. 2025). Moreover, the use of the aDNA protocol did not allow for the recovery of sufficient quality Symbiodiniaceae sequences from millennia-old coral

(Scott et al. 2022), suggesting that even optimized protocols for the recovery of shorter DNA fragments common in old and poorly preserved samples may not provide additional benefits for Symbiodiniaceae reconstructions.

Taken together, our results show that the DNA extraction protocol has an effect on the recovered Symbiodiniaceae communities, which has implications for comparability across studies. As such, we recommend the use of one consistent DNA extraction method with a decalcification step for all samples of a study to capture the additional diversity. Moreover, additional DNA extraction optimizations that improved yields in forensic mammalian bones could be evaluated in further studies. These include testing DNA concentrations in the EDTA decalcification solution and DNA enzymatic repair techniques to test their effects in the reconstructed Symbiodiniaceae communities from coral skeletons (Hofreiter et al. 2021).

## 4.2 | aDNA from Coral Cores Recovers Diverse Symbiodiniaceae Communities

Different Symbiodiniaceae communities were reconstructed from all cores, highlighting the community specificity of each coral species and between geographic locations. Coral species are known to associate with specialized Symbiodiniaceae communities (Davies et al. 2023), and within species, these communities may vary along spatial, temporal, or environmental gradients (Terraneo et al. 2019; de Souza et al. 2022; Lewis et al. 2024; Ng et al. 2024). Taken together, our results indicate that the reconstructions are sensitive to changes associated with taxonomic and geographical factors but not with fragment age, a novel finding that previous Symbiodiniaceae community reconstructions have not evaluated (González-Pech et al. 2024; Rey et al. 2025).

The recovered Symbiodiniaceae genera partially correspond to modern coral tissue samples. The dominance of *Cladocopium* in the tissue fragments of both *P. lobata* cores and the high relative abundance of *Durusdinium* in the tissue fragment of *D. heliopora* sampled in Palau is consistent with samples from other reefs of the Indo-Pacific, including Palau and Papua New Guinea, during the Tara Pacific expedition as well as long-term community monitoring studies (Hume et al. 2020; Lewis et al. 2024; Ng et al. 2024). Nevertheless, *Symbiodinium* was found in high abundances (> 50%) in most samples below the core top, regardless of coral species and sampling location, followed by a lower (~ 10%) but consistent abundance of *Breviolum*. The genus *Symbiodinium* has been detected to be present mostly as a low-abundant background and rarely as a dominant symbiont in multiple coral species, including *P. lobata* and *P. lutea* sampled in Palau and Papua New Guinea (Smith et al. 2017; Hume et al. 2020; Ng et al. 2024), partially explaining its occurrence in the skeleton samples. Notably, the genus *Breviolum* is more frequently reported from Atlantic coral reefs, where it has also been reconstructed from coral cores of *Orbicella faveolata* (Garrido et al. 2022; González-Pech et al. 2024). In the Pacific, it has only occasionally been associated with corals, reef water, and sediments in the Great Barrier Reef, with other cnidarians such as hydrocorals and anemones in the eastern and central Pacific, and with the
coral *Plesiastrea versipora* in the northern temperate Pacific (LaJeunesse et al. 2018; Ziegler et al. 2018, 2019; Nitschke et al. 2022). Similar to finding *Symbiodinium* in the coral skeletons, while not implausible, the detection of *Breviolum* raises questions about the process of Symbiodiniaceae deposition in coral skeletons.

There is currently no established model of Symbiodiniaceae DNA deposition in coral skeletons. The coral skeletal matrix is a low-light environment that experiences periodic anoxic conditions (Iha et al. 2021). Shallow coral skeletal bands may support a distinct endolithic community often composed of green algae from the genus *Ostreobium*, fungi, and bacteria (Tribollet et al. 2006; Iha et al. 2021). These endolithic organisms have adaptations to the challenging conditions within the skeleton, which Symbiodiniaceae cells do not seem to possess (Iha et al. 2021). One possible mechanism involves the accidental entrapment of the native Symbiodiniaceae communities from the coral tissue in the skeleton, e.g., when horizontal disseпipments are deposited during calcification or after injury. Whether and for how long these symbiont cells remain viable is unknown. A complementary explanation involves the entrapment of Symbiodiniaceae DNA in the aragonite matrix of the coral skeleton similar to the preferential DNA preservation in the columnar aragonite CaCO₃ layer of mollusk shells (Martin et al. 2021). Additionally, mollusk shells capture and preserve environmental DNA (eDNA) (der Sarkissian et al. 2017). The transport of seawater minerals to the calcifying fluid of corals might provide an explanation for the preservation of exogenous DNA in skeletons, although possible transport mechanisms need to be investigated (Tambutté et al. 2011; der Sarkissian et al. 2017; Thompson 2022). Our results suggest the deposition of multiple Symbiodiniaceae genera in coral skeletons. Further studies are needed to elucidate the underlying mechanisms that lead to the mismatch in relative abundances between genera between coral tissue and the underlying skeleton in aDNA-derived communities. For instance, the Symbiodiniaceae cells may use the skeleton as a habitat, or the cells or only their DNA may be trapped. An example of entrapment of exogenous organisms is the record of the dinoflagellate *Azadinium spinosum* in coral skeletons from aDNA-based reconstructions in *P. lobata*. *Azadinium spinosum* is not commonly associated with corals, and it has been hypothesized that its traces were incorporated in the skeleton from the surrounding seawater by ingestion during algae bloom events (Rey et al. 2025). The presence of diverse Symbiodiniaceae communities in the seawater (Claar et al. 2020; Million et al. 2024) may thus explain their detection in coral skeletons, and the record of these communities in coral skeletons highlights the potential to explore the use of coral skeletons for eDNA-based reconstructions in time.

**4.3 | Reconstructions from Coral Cores Capture Past Symbiodiniaceae Community Dynamics**

Symbiodiniaceae communities in cores from Papua New Guinea were relatively stable at decadal and centennial time scales, and only the tissue-associated core tops were distinct from the older skeletal samples. In contrast, both cores sampled in Palau suggest changes in the dominant symbiont community during the last decades. Notably, the large shifts in symbiont community

composition among time points were characterized by dynamics of sequences within the extant community and novel sequences, indicating a combination of symbiont shuffling and shifting as the underlying mechanism (Quigley et al. 2022).

Notably, overall higher SST and frequency of days with severe bleaching alerts were observed in Palau when compared to Papua New Guinea, suggesting a larger influence of historical heat stress events that could be associated with Symbiodiniaceae community dynamics (Lewis et al. 2024; Ng et al. 2024). The higher frequency of days with severe bleaching alerts coincided with the dominance of *Cladocopium* over *Symbiodinium* in *P. lobata* within the first 10 years of the reconstruction, which is similar to previously reported patterns in *O. faveolata* reconstructions (González-Pech et al. 2024). While *Porites* spp. are often tightly associated with the C15-type complex of *Cladocopium*, *P. lobata* was previously shown to undergo a lab-induced symbiont switch (Scharfenstein et al. 2022), as well as a high spatio-temporal structure in the reef (Ziegler et al. 2015), deeming a symbiont change plausible.

The occurrence of a severe alert level 3 bleaching episode in Palau coincided with a change in dominant symbiont in *P. lobata*. While the *Durusdinium* dominance in *D. heliopora* in Palau could not be linked to a NOAA bleaching alert, this sample was taken following the El Niño-associated bleaching in 1982/1983 (Coffroth et al. 1990). Given the low spatial and temporal resolution in environmental history and proof-of-concept sampling design, the causal link between environmental history and community dynamics remains to be confirmed. The subsequent sample from 5 years later indicates a lower but still dominant *Durusdinium* presence, which had reversed back to *Symbiodinium* dominance in the skeleton another 5 years later. This could reflect the reversal in symbiont communities to their original state, as reported from the field (Smith et al. 2017; Quigley et al. 2022). The coral skeleton-based Symbiodiniaceae community reconstruction cannot unequivocally confirm past changes of tissue-associated communities due to a lack of understanding of DNA entrapment mechanisms. However, as suggested by our data, reconstructions reveal patterns that may indicate changes in dominance of the community that were previously inaccessible due to a lack of long-term time-series samples. Moreover, by combining past ocean chemistry and climate reconstructions from the same cores employed for aDNA sampling, high-resolution correlations between Symbiodiniaceae dynamics and past environmental conditions may be established. The information derived from such multidisciplinary studies may be used to better inform projections of coral health and status under future ocean conditions.

**4.4 | Perspective**

Ancient DNA approaches to reconstruct microbial dynamics or host adaptation processes from coral skeletons are challenged by the environmental conditions during the lifetime of the coral that favor DNA degradation. In this study, aDNA as far back as 115 years was successfully recovered, validating the use of coral skeleton records for comparative analyses of long-term Symbiodiniaceae community dynamics. Similar to the reconstruction of past ocean conditions from geochemical proxies,
the reconstruction of Symbiodiniaceae communities is limited
by the available skeletal material, which often lacks biological
replicates per species and site, hindering the study of biological
patterns at the population level. Further research on the mech-
anisms of aDNA deposition in the coral skeleton will help to
interpret reconstruction dynamics and to understand natural
long-term adaptations to environmental stressors.

Author Contributions

**Jose F. Grillo:** formal analysis, investigation, writing – original draft, writing – review and editing. **Vanessa Tirpitz:** investigation. **Jessica Reichert:** investigation. **Marine Canesi:** investigation. **Stéphanie Reynaud:** resources. **Eric Douville:** resources. **Maren Ziegler:** conceptualization, investigation, methodology, project administration, resources, supervision, writing – original draft, writing – review and editing.

Acknowledgements

Special thanks to the Tara Ocean Foundation, the R/V Tara crew,
and the Tara Pacific Expedition Participants (https://doi.org/10.
5281/zenodo.3777760). We are keen to thank the commitment of
the following institutions for their financial and scientific support
that made this unique Tara Pacific Expedition possible: CNRS, PSL,
CSM, EPHE, Genoscope, CEA, Inserm, Université Côte d'Azur,
ANR, agnès b., UNESCO-IOC, the Veolia Foundation, the Prince
Albert II de Monaco Foundation, Région Bretagne, Billerudkorsnas,
AmerisourceBergen Company, Lorient Agglomération, Oceans by
Disney, L'Oréal, Biotherm, France Collectivités, Fonds Français
pour l'Environnement Mondial (FFEM), Etienne Bourgois, and the
Tara Ocean Foundation teams. Tara Pacific would not exist without
the continuous support of the participating institutes. The authors
also particularly thank Serge Planes, Denis Allemand, and the Tara
Pacific consortium. The authors wish to thank Guillaume Iwankow
for the sampling of the coral cores, Silvia Nachtigall for support with
lab work, and the Marine Holobiomics Lab for discussions of the
data. We are further keen to thank the local people and authorities
who made this study possible. Specifically, we thank the Ministry
of Natural Resources, Environment and Tourism of the Republic
of Palau, and in Papua New Guinea we thank the Marine Science
Research Committee (MSRC), the Department of Foreign Affairs. and
the University of Papua New Guinea. This is publication number #40
of the Tara Pacific Consortium. Open Access funding enabled and or-
ganized by Projekt DEAL.

Conflicts of Interest

The authors declare no conflicts of interest.

Data Availability Statement

The data and code that support the findings of this study are openly
available in Zenodo at https://doi.org/10.5281/zenodo.17246431. The
daily bleaching alert data were obtained from the NOAA Coral Reef
Watch at https://pae-paha.pacioos.hawaii.edu/erddap/griddap/dhw_
5km.html and are archived with the raw data for this study. All se-
quences are accessible under the NCBI bioproject ID PRJNA1330243.

References

Al-Horani, F. A. 2015. "Physiology of Skeletogenesis in Scleractinian Coral." In *Diseases of Coral*, edited by C. M. Woodley, C. A. Downs, A. W. Bruckner, J. W. Porter, and S. B. Galloway, 192–205. Wiley.

Anderson, M. J. 2017. "Permutational Multivariate Analysis of Variance (PERMANOVA)." In *Wiley StatsRef: Statistics Reference Online*, edited

by N. Balakrishnan, T. Colton, B. Everitt, W. Piegorsch, F. Ruggeri, and
J. L. Teugels, 1–15. Wiley.

Baker, A. C. 2003. "Flexibility and Specificity in Coral-Algal Symbiosis: Diversity, Ecology, and Biogeography of *Symbiodinium*." *Annual Review of Ecology, Evolution, and Systematics* 34: 661–689.

Baker, A. C., C. J. Starger, T. R. McClanahan, and P. W. Glynn. 2004. "Corals' Adaptive Response to Climate Change." *Nature* 430: 741.

Baker, D. M., L. Weigt, M. Fogel, and N. Knowlton. 2013. "Ancient DNA From Coral-Hosted Symbiodinium Reveal a Static Mutualism Over the Last 172 Years." *PLoS One* 8: e55057.

Boulotte, N. M., S. J. Dalton, A. G. Carroll, et al. 2016. "Exploring the *Symbiodinium* Rare Biosphere Provides Evidence for Symbiont Switching in Reef-Building Corals." *ISME Journal* 10: 2693–2701.

Buddemeier, R. W., and D. G. Fautin. 1993. "Coral Bleaching as an Adaptive Mechanism." *Bioscience* 43: 320–326.

Canesi, M., E. Douville, P. Montagna, et al. 2023. "Differences in Carbonate Chemistry Up-Regulation of Long-Lived Reef-Building Corals." *Scientific Reports* 13: 11589.

Canesi, M., E. Douville, P. Montagna, et al. 2024. "Sea Surface Temperature Reconstruction in the Pacific Ocean Using Multi-Elemental Proxy in Porites and Diploastrea Corals: Application to Palau Archipelago." *Chemical Geology* 645: 121884.

Claar, D. C., K. L. Tietjen, K. D. Cox, R. D. Gates, and J. K. Baum. 2020. "Chronic Disturbance Modulates Symbiont (Symbiodiniaceae) Beta Diversity on a Coral Reef." *Scientific Reports* 10: 4492.

Coffroth, M. A., H. R. Lasker, and J. K. Oliver. 1990. "Coral Mortality Outside of the Eastern Pacific During 1982-1983: Relationship to El Niño." In *Elsevier Oceanography Series*, 141-182. Elsevier.

Connelly, M. T., M. G. Catapang, and A. M. Quattrini. 2024. "Unlocking the Treasure Trove: Leveraging Dry Coral Specimens for Museum Genomics." *Coral Reefs* 43: 1153-1159.

Davies, S. W., M. H. Gamache, L. I. Howe-Kerr, et al. 2023. "Building Consensus Around the Assessment and Interpretation of Symbiodiniaceae Diversity." *PeerJ* 11: e15023.

de Souza, M. R., C. Caruso, L. Ruiz-Jones, C. Drury, R. Gates, and R. J. Toonen. 2022. "Community Composition of Coral-Associated Symbiodiniaceae Differs Across Fine-Scale Environmental Gradients in Kāne Ohe Bay." *Royal Society Open Science* 9: 212042.

der Sarkissian, C., V. Pichereau, C. Dupont, et al. 2017. "Ancient DNA Analysis Identifies Marine Mollusc Shells as New Metagenomic Archives of the Past." *Molecular Ecology Resources* 17: 835–853.

Edinger, E. N., G. V. Limmon, J. Jompa, W. Widjatmoko, J. M. Heikoop, and M. J. Risk. 2000. "Normal Coral Growth Rates on Dying Reefs: Are Coral Growth Rates Good Indicators of Reef Health?" *Marine Pollution Bulletin* 40: 404–425.

Frölicher, T. L., E. M. Fischer, and N. Gruber. 2018. "Marine Heatwaves Under Global Warming." *Nature* 560: 360–364.

Garrido, A. G., D. C. de Assis Leite, L. F. Machado, R. S. Peixoto, and C. Zilberberg. 2022. "Highly Diverse and Geographically Differentiated Symbiodiniaceae Communities Associated With the Hydrocoral *Millepora alcicornis* in the Atlantic Ocean." *Coral Reefs* 41: 937–950.

González-Pech, R. A., C. Howe, L. D. Lizcano-Sandoval, et al. 2024. *A Paleogenomic Approach to Reconstruct Historical Responses of Coral Reefs to Anthropogenic Change*. bioRxiv.

Helgoe, J., S. K. Davy, V. M. Weis, and M. Rodriguez-Lanetty. 2024. "Triggers, Cascades, and Endpoints: Connecting the Dots of Coral Bleaching Mechanisms." *Biological Reviews* 99: 715–752.

Hofreiter, M., J. Sneberger, M. Pospisek, and D. Vanek. 2021. "Progress in Forensic Bone DNA Analysis: Lessons Learned From Ancient DNA." *Forensic Science International: Genetics* 54: 102538.
Hume, B. C. C., J. Poulain, S. Pesant, et al. 2020. Tara Pacific ITS2 Symbiodiniaceae Data Release Version 1.

Hume, B. C. C., E. G. Smith, M. Ziegler, et al. 2019. "SymPortal: A Novel Analytical Framework and Platform for Coral Algal Symbiont Next-Generation Sequencing ITS2 Profiling." *Molecular Ecology Resources* 19: 1063–1080.

Hume, B. C. C., M. Ziegler, J. Poulain, et al. 2018. "An Improved Primer Set and Amplification Protocol With Increased Specificity and Sensitivity Targeting the *Symbiodinium* ITS2 Region." *PeerJ* 6: e4816.

Iha, C., K. E. Dougan, J. A. Varela, et al. 2021. "Genomic Adaptations to an Endolithic Lifestyle in the Coral-Associated Alga Ostreobium." *Current Biology* 31: 1393–1402.

Klein, R., and Y. Loya. 1991. "Skeletal Growth and Density Patterns of Two Pontes Corals From the Gulf of Eilat, Red Sea." *Marine Ecology Progress Series* 77: 253–259.

Krassowski, M., M. Arts, and C. Lagger. 2022. Krassowski/Complex-Upset: v1.3.5.

LaJeunesse, T. C., J. E. Parkinson, P. W. Gabrielson, et al. 2018. "Systematic Revision of Symbiodiniaceae Highlights the Antiquity and Diversity of Coral Endosymbionts." *Current Biology* 28: 2570–2580.

Lendvay, B., L. E. Cartier, M. Gysi, et al. 2020. "DNA Fingerprinting: An Effective Tool for Taxonomic Identification of Precious Corals in Jewelry." *Scientific Reports* 10: 8287.

Lewis, A. M., C. C. Butler, K. E. Turnham, et al. 2024. "The Diversity, Distribution, and Temporal Stability of Coral 'Zooxanthellae' on a Pacific Reef: From the Scale of Individual Colonies to Across the Host Community." *Coral Reefs* 43: 841–856.

Martin, K. R., L. P. Waits, and C. E. Parent. 2021. "Teaching an Old Shell New Tricks: Extracting DNA From Current, Historical, and Ancient Mollusk Shells." *Bioscience* 71: 235–248.

Martinez Arbizu, P. 2017. *pairwiseAdonis: Pairwise Multilevel Comparison Using Adonis*.

Martin-Roy, R., J. Thyrring, X. Mata, et al. 2024. "Advancing Responsible Genomic Analyses of Ancient Mollusc Shells." *PLoS One* 19: e0302646.

Masuda, S., T. Awaji, N. Sugiura, et al. 2010. "Simulated Rapid Warming of Abyssal North Pacific Waters." *Science* 329: 319–322.

Million, W. C., C. R. Voolstra, G. Perna, G. Puntin, K. Rowe, and M. Ziegler. 2024. Resolving Symbiodiniaceae Diversity Across Coral Microhabitats and Reef Niches.

Ng, M. S., N. Soon, L. Afiq-Rosli, et al. 2024. "Highly Diverse Symbiodiniaceae Types Hosted by Corals in a Global Hotspot of Marine Biodiversity." *Microbial Ecology* 87: 92.

Nitschke, M. R., S. L. Rosset, C. A. Oakley, et al. 2022. "The Diversity and Ecology of Symbiodiniaceae: A Traits-Based Review." *Advances in Marine Biology* 55: 127.

Oksanen, J., G. L. Simpson, F. G. Blanchet, et al. 2024. *Vegan: Community Ecology Package*.

Pettay, D. T., D. C. Wham, R. T. Smith, R. Iglesias-Prieto, and T. C. LaJeunesse. 2015. "Microbial Invasion of the Caribbean by an Indo-Pacific Coral Zooxanthella." *Proceedings of the National Academy of Sciences of the United States of America* 112: 7513–7518.

Planes, S., D. Allemand, S. Agostini, et al. 2019. "The Tara Pacific expe-
dition—A pan-ecosystemic approach of the “-omics” complexity of coral
reef holobionts across the Pacific Ocean." *PLOS Biology* 17: e3000483.

Quigley, K. M., S. W. Davies, C. D. Kenkel, B. L. Willis, M. V. Matz, and L. K. Bay. 2014. "Deep-Sequencing Method for Quantifying Background Abundances of Symbiodinium Types: Exploring the Rare Symbiodinium Biosphere in Reef-Building Corals." *PLoS One* 9: e94297.

Quigley, K. M., B. Ramsby, P. Laffy, J. Harris, V. J. L. Mocellin, and L. K. Bay. 2022. "Symbioses Are Restructured by Repeated Mass Coral Bleaching." *Science Advances* 8: eabq8349.

R Core Team. 2023. R: A Language and Environment for Statistical Computing.

Ramos, R., R. Cipriani, H. M. Guzman, and E. Garcia. 2009. "Chronology of Mercury Enrichment Factors in Reef Corals From Western Venezuela." *Marine Pollution Bulletin* 58: 222–229.

Raymundo, L. J., D. Burdick, W. C. Hoot, et al. 2019. "Successive Bleaching Events Cause Mass Coral Mortality in Guam, Micronesia." *Coral Reefs* 38: 677–700.

Reimer, J. D., R. S. Peixoto, S. W. Davies, et al. 2024. "The Fourth Global Coral Bleaching Event: Where Do We Go From Here?" *Coral Reefs* 43: 1121–1125.

Rey, O., D. Dissard, E. Toulza, et al. 2025. "Reconstructing the Dynamics of Past Coral Endosymbiotic Algae Communities Using Coral Ancient DNA (coraDNA)." *Coral Reefs* 44: 529–539.

Riegl, B., M. Johnston, S. Purkis, et al. 2018. "Population Collapse Dynamics in *Acropora downingi*, an Arabian/Persian Gulf Ecosystem-Engineering Coral, Linked to Rising Temperature." *Global Change Biology* 24: 2447–2462.

Rohland, N., I. Glocke, A. Aximu-Petri, and M. Meyer. 2018. "Extraction of Highly Degraded DNA From Ancient Bones, Teeth and Sediments for High-Throughput Sequencing." *Nature Protocols* 13: 2447–2461.

Rowe, K. C., S. Singhal, M. D. Macmanes, et al. 2011. "Museum Genomics: Low-Cost and High-Accuracy Genetic Data From Historical Specimens." *Molecular Ecology Resources* 11: 1082–1092.

Scharfenstein, H. J., W. Y. Chan, P. Buerger, C. Humphrey, and M. J. H. van Oppen. 2022. "Evidence for de Novo Acquisition of Microalgal Symbionts by Bleached Adult Corals." *ISME Journal* 16: 1676–1679.

Scott, C. B., A. Cárdenas, M. Mah, et al. 2022. "Millennia-Old Coral Holobiont DNA Provides Insight Into Future Adaptive Trajectories." *Molecular Ecology* 31: 4979–4990.

Skirving, W., B. Marsh, J. de la Cour, et al. 2020. "CoralTemp and the Coral Reef Watch Coral Bleaching Heat Stress Product Suite Version 3.1." *Remote Sensing* 12: 3856.

Smith, E. G., G. O. Vaughan, R. N. Ketchum, D. McParland, and J. A. Burt. 2017. "Symbiont Community Stability Through Severe Coral Bleaching in a Thermally Extreme Lagoon." *Scientific Reports* 7: 2428.

Swain, T. D., J. Chandler, V. Backman, and L. Marcelino. 2017. "Consensus Thermotolerance Ranking for 110 *Symbiodinium* Phylotypes: An Exemplar Utilization of a Novel Iterative Partial-Rank Aggregation Tool With Broad Application Potential." *Functional Ecology* 31: 172–183.

Tambutté, S., M. Holcomb, C. Ferrier-Pagès, et al. 2011. "Coral Biomineralization: From the Gene to the Environment." *Journal of Experimental Marine Biology and Ecology* 408: 58–78.

Terraneo, T. I., M. Fusi, B. C. C. Hume, et al. 2019. "Environmental Latitudinal Gradients and Host-Specificity Shape Symbiodiniaceae Distribution in Red Sea *Porites* Corals." *Journal of Biogeography* 46: 2323–2335.

Thompson, D. M. 2022. "Environmental Records From Coral Skeletons: A Decade of Novel Insights and Innovation." *WIREs Climate Change* 13: e745.

Tribollet, A., C. Langdon, S. Golubic, and M. Atkinson. 2006. "Endolithic Microflora Are Major Primary Producers in Dead Carbonate Substrates of Hawaiian Coral Reefs." *Journal of Phycology* 42: 292–303.

Waller, R. G., J. F. Adkins, L. F. Robinson, and T. M. Shank. 2007. "Ancient DNA Techniques: Applications for Deep-Water Corals." *Bulletin of Marine Science* 81: 351–360.
Wickham, H., M. Averick, J. Bryan, et al. 2019. "Welcome to the Tidyverse." *Journal of Open Source Software* 4: 1686.

Xiao, L., S. Johansson, S. Rughöft, et al. 2022. "Photophysiological Response of Symbiodiniaceae Single Cells to Temperature Stress." *ISME Journal* 16: 2060–2064.

Yang, D. Y., and K. Watt. 2005. "Contamination Controls When Preparing Archaeological Remains for Ancient DNA Analysis." *Journal of Archaeological Science* 32: 331–336.

Ziegler, M., C. Arif, and C. R. Voolstra. 2019. "Symbiodiniaceae Diversity in Red Sea Coral Reefs & Coral Bleaching." In *Coral Reefs of the Red Sea*, edited by C. R. Voolstra and M. L. Berumen, 69–89. Springer International Publishing.

Ziegler, M., C. Roder, C. Büchel, and C. Voolstra. 2015. "Niche Acclimatization in Red Sea Corals Is Dependent on Flexibility of Host-Symbiont Association." *Marine Ecology Progress Series* 533: 149–161.

Ziegler, M., E. Stone, D. Colman, C. Takacs-Vesbach, and U. Shepherd. 2018. "Patterns of *Symbiodinium* (Dinophyceae) Diversity and Assemblages Among Diverse Hosts and the Coral Reef Environment of Lizard Island, Australia." *Journal of Phycology* 54: 447–460.

## Supporting Information

Additional supporting information can be found online in the Supporting Information section. **Table S1:** Post-MED sequence number of each core fragment for *P. lobata* and *D. heliopora*. **Table S2:** Pairwise statistical test of differences in the reconstructed Symbiodiniaceae. **Table S3:** Statistical test of multivariate dispersion in the reconstructed Symbiodiniaceae. **Table S4:** Statistical test of differences in the reconstructed Symbiodiniaceae communities of. **Table S5:** Pairwise statistical test of differences in the reconstructed Symbiodiniaceae. **Table S6:** Statistical test of multivariate dispersion in the reconstructed Symbiodiniaceae. **Table S7:** Mantel test performed on Symbiodiniaceae community similarities calculated with. **Figure S1:** Photographs of the coral skeleton cores. (A) Halved coral core of *Porites lobata*. **Figure S2:** Distance-decay linear regression fitted on community similarity Bray-Curtis. **Figure S3:** Number of shared post-MED sequences between *P. lobata* and *D. heliopora*. **Figure S4:** Daily sea surface temperature and days with bleaching alerts for the sampling.