# Genetic diversity increases with depth in red gorgonian populations of the Mediterranean Sea and the Atlantic Ocean

Joanna Pilczynska, Silvia Cocito, Joana Ruela Heimbürger Boavida, Ester A Serrão, Jorge Assis, Eliza Fragkopoulou, Henrique Queiroga

► To cite this version:

Joanna Pilczynska, Silvia Cocito, Joana Ruela Heimbürger Boavida, Ester A Serrão, Jorge Assis, et al.. Genetic diversity increases with depth in red gorgonian populations of the Mediterranean Sea and the Atlantic Ocean. PeerJ, 2019, 7, pp.e6794. 10.7717/peerj.e6794. hal-03043794

HAL Id: hal-03043794

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

Submitted on 26 Jun 2025

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 Attribution 4.0 International License
# Genetic diversity increases with depth in red gorgonian populations of the Mediterranean Sea and the Atlantic Ocean

Joanna Pilczynska<fnref n="1" />,<fnref n="2" />, Silvia Cocito<fnref n="3" />, Joana Boavida<fnref n="4" />,<fnref n="5" />, Ester A. Serrão<fnref n="4" />, Jorge Assis<fnref n="4" />, Eliza Fragkopoulou<fnref n="4" /> and Henrique Queiroga<fnref n="1" />

## ABSTRACT

**Background.** In the ocean, the variability of environmental conditions found along depth gradients exposes populations to contrasting levels of perturbation, which can be reflected in the overall patterns of species genetic diversity. At shallow sites, resource availability may structure large, persistent and well-connected populations with higher levels of diversity. In contrast, the more extreme conditions, such as thermal stress during heat waves, can lead to population bottlenecks and genetic erosion, inverting the natural expectation. Here we examine how genetic diversity varies along depth for a long-lived, important ecosystem-structuring species, the red gorgonian, *Paramuricea clavata*.

**Methods.** We used five polymorphic microsatellite markers to infer differences in genetic diversity and differentiation, and to detect bottleneck signs between shallow and deeper populations across the Atlantic Ocean and the Mediterranean Sea. We further explored the potential relationship between depth and environmental gradients (temperature, ocean currents, productivity and slope) on the observed patterns of diversity by means of generalized linear mixed models.

**Results.** An overall pattern of higher genetic diversity was found in the deeper sites of the Atlantic Ocean and the Mediterranean Sea. This pattern was largely explained by bottom temperatures, with a linear pattern of decreasing genetic diversity with increasing thermal stress. Genetic differentiation patterns showed higher gene flow within sites (i.e., shallow vs. deeper populations) than between sites. Recent genetic bottlenecks were found in two populations of shallow depths.

**Discussion.** Our results highlight the role of deep refugial populations safeguarding higher and unique genetic diversity for marine structuring species. Theoretical regression modelling demonstrated how thermal stress alone may reduce population sizes and diversity levels of shallow water populations. In fact, the examination of time series on a daily basis showed the upper water masses repeatedly reaching lethal temperatures for *P. clavata*. Differentiation patterns showed that the deep richer populations are

<footnote n="1">¹ Departamento de Biologia and CESAM—Centro de Estudos do Ambiente e do Mar, Universidade de Aveiro, Aveiro, Portugal</footnote>

<footnote n="2">² Department of Earth and Environmental Sciences, University of Pavia, Pavia, Italy</footnote>

<footnote n="3">³ Italian Agency for New Technologies, Energy and Sustainable Economic Development—ENEA, Marine Environment Research Centre, La Spezia, Italy</footnote>

<footnote n="4">⁴ CCMAR—Centro de Ciências do Mar, Universidade do Algarve, Faro, Portugal</footnote>

<footnote n="5">⁵ Aix Marseille Université, CNRS/INSU, Université de Toulon, IRD, Mediterranean Institute of Oceanography (MIO) UM 110, Marseille, France</footnote>
isolated. Gene flow was also inferred across different depths; however, not in sufficient levels to offset the detrimental effects of surface environmental conditions on genetic diversity. The identification of deep isolated areas with high conservation value for the red gorgonian represents an important step in the face of ongoing and future climate changes.

**Subjects** Biodiversity, Ecology, Genetics, Marine Biology, Climate Change Biology

**Keywords** Genetic diversity, Depth refugia, Mediterranean Sea, Atlantic Ocean, Paramuricea clavata, Climate change

## INTRODUCTION

Extreme environmental conditions may change the distribution of intra-specific biodiversity (Provan & Bennett, 2008). However, responses may differ significantly between environments and ecological groups. The trends in genetic diversity of mountain plants and vertebrates along altitudinal gradients are a well-known example, varying from decreased diversity with altitude due to drift and bottlenecks during vertical range expansion, to increased diversity associated with selective pressures at higher altitudes (Giordano, Ridenhour & Storfer, 2007; Ohsawa & Ide, 2008).

In the marine environment, environmental gradients are known to affect the genetic diversity levels of the populations (e.g., Costantini et al., 2011; Costantini et al., 2016; Johannesson & André, 2006). Particularly, the exposure to limiting niche conditions may reduce population sizes, leading to genetic erosion through bottlenecks and drift (Eckert, Samis & Lougheed, 2008). Conversely, where conditions are stable for long time, populations may persist and retain ancient genetic diversity (i.e., climatic refugia; Maggs et al., 2008; Provan & Bennett, 2008). Even small populations of species with reduced dispersal potential may harbour distinct ancient genetic diversity if persisting in refugial areas (e.g., Diekmann & Serrão, 2012; Assis et al., 2016).

In the coastal zone, depth gradients are associated with increased environmental variability due to the stratification of wind-induced turbulence, light attenuation, nutrient availability, sedimentation and the presence of thermoclines and haloclines (Garrabou, Ballesteros & Zabala, 2002; Assis et al., 2017). Together, they underpin well-structured gradients of species occurrences, abundances and genetic diversity. While extreme conditions experienced in shallow waters can be detrimental to genetic diversity, deep waters are relatively more stable and can provide genetic refugia (e.g., Smith et al., 2014; Assis et al., 2016). The hypothesis of deep populations harbouring higher and unique genetic diversity is of great conservational, biogeographical and evolutionary relevance. The loss of such genetically rich populations poses disproportionate risks for the species as a whole, considering the loss of adaptive variation for selection (Hampe & Petit, 2005). While there is increased evidence of deep populations not being directly affected by events of extreme environmental conditions (e.g., Cerrano et al., 2005; Linares et al., 2005; Smith et al., 2014), these are not immune to disturbances (Bavestrello et al., 2014; Frade et al., 2018). Furthermore, the hypothesis that deeper reefs may be more fecund (Holstein et al.,
2015), serving as a source of recruits for the recovery of shallower reefs is controversial between species (*Bongaerts et al., 2010; Bongaerts et al., 2017*), but also within the same species (*Mokhtar-Jamaï et al., 2011; Van Oppen et al., 2011; Pilczynska et al., 2016*).

The red gorgonian *Paramuricea clavata* (Risso) lives on shadowed rocky substrates down to 120 m in the Mediterranean Sea (*Salomidi et al., 2009; Bo et al., 2012*) and to 100 m in the Atlantic Ocean (*Boavida et al., 2016a*). Shallow populations in the Mediterranean have recently suffered massive mortality events caused by thermal stress (*Perez et al., 2000; Romano et al., 2000; Garrabou et al., 2009*). Damage intensity has been reported to decrease with depth, with communities dwelling below the thermocline of 25 to 30 m being less affected (*Cerrano et al., 2005; Linares et al., 2005*). In the Atlantic Ocean, temperature-driven mortality events are presumed to be less frequent as water masses are usually mixed due to summer upwelling (*Relvas et al., 2007*). However, there are sporadic events of upwelling relaxation, with raising temperatures persisting for several days, responsible for the mortality of shallow water gastropod populations (*Lima et al., 2006*) and marine forests (*Araújo et al., 2016*), which eventually may reach levels that could become limiting for the red gorgonian. Here we investigate changes in the genetic diversity levels of *P. clavata* across populations in different depth ranges from the Mediterranean Sea and the Atlantic Ocean. Considering the different potential effect of temperature extremes between different depths, we hypothesize that shallow populations have lower genetic diversity when compared to deeper ones. We further explore the observed patterns of diversity using generalized linear mixed models with environmental information (temperature, current velocity, primary productivity and slope).

## MATERIALS & METHODS

We sampled two populations at different depths (shallow vs. deeper) in three sites of the Atlantic Ocean (Portuguese coast) and two sites of the Mediterranean Sea (Italy coast; Fig. 1; Table S1) by means of SCUBA diving. Atlantic samples were located between 12 m and 60 m depth, while those in the Mediterranean were between 20 m and 30 m. At each population, 10 cm apical branches were collected haphazardly from well-separated colonies to avoid clones, as the species can have asexual reproduction, though negligible (*Coma, Zabala & Gili, 1995; Pilczynska et al., 2017*). The branch tip from each colony was stored individually in plastic tube underwater and, after transportation to the lab, the samples were preserved in 96% ethanol until DNA extraction. Sampling permissions for the Mediterranean sites were authorised by Cinque Terre Marine Protected Area (*Pilczynska et al., 2016*). For the Atlantic sites, permission was authorized by the Institution for Conservation of Nature and Forest (ICNF, Portuguese governmental body responsible for the management of Protected Areas).

Genomic DNA was extracted with a CTAB protocol (*Winnepenninckx, Backeljau & De Wachter, 1993*), using proteinase K, with purification by standard chloroform:isoamyl alcohol (24:1) followed by DNA precipitation. Samples were genotyped at five microsatellite loci: Parcla 09, Parcla 10, Parcla 12, Parcla 14, Parcla 17 (*Mokhtar-Jamaï et al., 2011*). PCR conditions were as described in *Mokhtar-Jamaï et al. (2011)* with minor modifications (1.5
<figure><img src="image_2.png" /><figcaption>**Figure 1** Sampling sites of the red gorgonian (*Paramuricea clavata*) populations in the Atlantic Ocean and the Mediterranean Sea. The exact coordinates of each sampling site are available in the supplementary information (S3). Map source "©OpenStreetMap contributors", available under the Open Database License. This figure is published under CC BY SA: https://www.openstreetmap.org/copyright.</figcaption></figure>

Full-size 🗠 DOI: 10.7717/peerj.6794/fig-1

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

mM MgCl2, cycle: 95 °C 3 min, 94 °C 20 s, 45 °C 20 s, 72 °C 20 s for 40 cycles, final extension 72 °C 10 min). PCR products were analysed on an ABI 3730XL Genetic Analyser using an internal size standard (GeneScan 600 LIZ; Applied Biosystems). STRand version 2.2.30 was used to score alleles (*Locke, Baack & Toonen, 2000*) and the R package MsatAllel_1.02 (*Alberto, 2009*) allowed to visualise, track and reanalyse putative scoring errors.

Additional genetic data for the same microsatellite loci were compiled from Mokhtar-Jamaï et al. (2011). Extra sites were chosen when comprising two populations at different depths (shallow vs. deeper). The new data (five sites) were located in the Spanish and French Mediterranean coasts, between 15 m and 40 m depth (Fig. 1; Table S1).

The shallow and deeper populations sampled in the Atlantic and Mediterranean do not coincide with the same exact depths, as vertical distribution limits varied between sites. Thus, comparisons of genetic diversity levels between shallow and deeper populations were made at the site scale, and not between sites. Allelic richness, number of private alleles and gene diversity (expected heterozygosity) per population (shallow and deeper) were standardized to the smallest sample size found within sites using 10⁴ randomizations. Significant differences in mean diversity levels within each site were tested using a non-parametric Wilcoxon signed-rank test with 10⁴ randomizations (Assis et al., 2018).

To infer the drivers shaping genetic variation among sites, genetic diversity estimates were standardized to the smallest size of all samples using 10⁴ randomizations. The estimates of diversity were modelled with linear regression against depth (null model) and a set of important environmental predictors affecting the physiology of *P. clavata* (Maximum bottom temperature; *Boavida et al., 2016b*) and proxies of essential resources (minimum bottom productivity, minimum bottom current velocity and slope, for food intake; *Boavida et al., 2016b*). Other important predictors were not considered since (1) they did
not vary between sampling sites (e.g., salinity) or (2) were correlated with productivity (e.g., inorganic nutrients such as nitrates and phosphates). The predictors were developed with three-dimensional profiles of monthly data compiled from the Global Ocean Physics Reanalysis (ORAP; http://www.marine.copernicus.eu/) and the Biogeochemistry Non-Assimilative Hindcast Simulation (PISCES; http://www.marine.copernicus.eu/). Bottom environmental data for each population was obtained using trilinear interpolation (Assis et al., 2017) weighting location (longitude and latitude) and depth. Long-term minimum (productivity and currents) and maximum (temperature) extremes were averaged for the years 2000–2014. Slope was computed using the “terrain” function of the R package “raster” (3.5.2 version; R Core Team, 2018) in bathymetry.

Given the sampling design (no precise coincidence of depths between sites) we adopted the Generalized Linear Mixed Models (GLMM; Bolker et al., 2009) framework. This is suitable for modelling unbalanced designs (Zhang & Chen, 2013) and complex spatial and temporal correlation structures since it accounts for dependencies within hierarchical groups by introducing random effects (Bolker et al., 2009). While these effects generally comprise the blocking in experimental treatments, they can be used to encompass variation among geographical regions (Dormann et al., 2007). Accordingly, “site” was included in the models as a random effect term (Zuur, Ieno & Smith, 2007; Ludwig et al., 2012). The best-fit selection between the null and the environmental models followed the relative Akaike Information Criterion (AIC), while goodness-of-fit was inferred with the R2GLMM algorithm from Edwards et al. (2008). Partial dependence plots were also produced to illustrate the effect of predictors on the response of models, by accounting for the mean effect of all other predictors (Elith, Leathwick & Hastie, 2008). All models were performed in the R environment (3.5.2 version; R Core Team, 2018).

To infer larval dispersal potential, pairwise genetic differentiation was estimated between and within sampled populations (i.e., shallow vs. deeper populations) using $F_{ST}$ estimator and an analysis of molecular variance (AMOVA) based on allele frequencies, computed with Genodive (Meirmans & Van Tienderen, 2004) under $10^4$ randomizations (Assis et al., 2013).

Evidence of genetic bottlenecks were inferred by testing for heterozygosity excess (Piry, Luikart & Cornuet, 1999). This is rooted on the assumption that populations that have recently experienced a bottleneck event are predicted to temporarily reduce allelic diversity at a faster rate than heterozygosity. Such an excess in heterozygosity rate was tested for each population with the software Bottleneck (Piry, Luikart & Cornuet, 1999) using $10^4$ simulations. We used the suggested and more appropriate (realistic) parameters for microsatellites: the Two-Phase Model (TPM) (Luikart & Cornuet, 1998; Piry, Luikart & Cornuet, 1999) with a step in mutations of 0.9 (ps) and a variance in mutations of 12 (Piry, Luikart & Cornuet, 1999; Busch, Waser & DeWoody, 2007). Because the genetic dataset has less than 20 loci, we used the Wilcoxon test to address the null hypothesis of no heterozygosity excess (on average) across loci (Cornuet & Luikart, 1996; Luikart & Cornuet, 1998).

## RESULTS

Genetic diversity as allelic richness (A) and number of private alleles (PA) was consistently higher in deeper populations, with the exception of only one site (Petit Congloue; Figs. 2A,
<figure><img src="image_4.png" /><figcaption>**Figure 2** Genetic diversity (mean ± SD) as (A) allelic richness, (B) private alleles and (C) expected heterozygosity of *Paramuricea clavata* in shallow (light grey) and deeper (dark grey) populations. Asterisks indicate significant differences in diversity levels (*P* < 0.05). Sampling site name, depth (*m*) and number of samples (*n*) are described for each population.</figcaption></figure>

Full-size DOI: 10.7717/peerj.6794/fig-2

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

2B). The greatest difference in A and PA between shallow and deeper populations was observed in Pharillions and Sagres (A: 2.7 and 2.07; PA: 13.52 and 10.33, respectively), while the smallest was observed in Livorno (A: 0.09 and PA: 0.43). The variation in expected heterozygosity (He) did not follow a clear pattern, with only one sampling site (Grotte Peres) having a higher He in its deeper population (Fig. 2C).

The null linear regression models using depth alone were outperformed by the models using environmental predictors (environmental models showed lower AIC and higher R²; Table 1). Both approaches found significant relationships while fitting the predictors against standardized allelic richness and number of private alleles (Table 1), however, they failed to explain the variability in expected heterozygosity. The environmental models showing better fitting were mostly explained by temperature alone (Table 1). This predictor
<table><caption>**Table 1** Summary of linear regression models testing genetic diversity indices against depth and a set of environmental predictors. Akaike Information Criterion, R-squared and *p*-value scores are shown for each model. Bold represents higher values when comparisons were made.</caption><thead><tr><th rowspan="2">Model type</th><th rowspan="2">Predictors</th><th colspan="3">Allelic richness</th><th colspan="3">Private alleles</th><th colspan="3">Expected heterozygosity</th></tr><tr><th>AIC</th><th>R<sup>2</sup></th><th>p-value</th><th>AIC</th><th>R<sup>2</sup></th><th>p-value</th><th>AIC</th><th>R<sup>2</sup></th><th>p-value</th></tr></thead><tbody><tr><td>Null</td><td>Depth</td><td>71.05</td><td>0.22</td><td>0.01</td><td>92.84</td><td>0.60</td><td>0.01</td><td>-35.26</td><td>0.01</td><td>0.77</td></tr><tr><td rowspan="4">Environmental</td><td>Model effect</td><td><strong>46.83</strong></td><td><strong>0.63</strong></td><td></td><td><strong>74.31</strong></td><td><strong>0.62</strong></td><td></td><td><strong>-39.33</strong></td><td><strong>0.43</strong></td><td></td></tr><tr><td>Min. current velocity</td><td></td><td></td><td>0.98</td><td></td><td></td><td>0.45</td><td></td><td></td><td>0.29</td></tr><tr><td>Min. productivity</td><td></td><td></td><td>0.99</td><td></td><td></td><td>0.95</td><td></td><td></td><td>0.81</td></tr><tr><td>Slope</td><td></td><td></td><td>0.72</td><td></td><td></td><td>0.56</td><td></td><td></td><td>0.13</td></tr><tr><td></td><td>Max. temperature</td><td></td><td></td><td>0.01</td><td></td><td></td><td>0.01</td><td></td><td></td><td>0.53</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Figure 3** Partial dependence functions depicting the effect of bottom temperature on (A) allelic richness, (B) private alleles and (C) expected heterozygosity.</figcaption></figure>

Full-size DOI: 10.7717/peerj.6794/fig-3

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

produced a negative response on the models, with lower allelic richness and private alleles with increasing bottom temperatures (Fig. 3).

Genetic differentiation ($F_{ST}$) was higher between than within sampling sites, with two exceptions for Sagres and Riou Sud (Fig. 4). The former site (i.e., Sagres) showed the highest differentiation of all pairwise comparisons at the site scale (i.e., shallow vs. deeper sample; $F_{ST} \approx 0.4$; Table S2). No significant differentiation was found within 7 Mediterranean sites (Punta Mesco, Livorno, Pota Del Llop, Grotte Peres and Petit Congloue; Table S2) and between the two populations of Grotte Peres and the deeper population of Petit Congloue (Table S2), which are approximately 2 km far apart.

Bottleneck events (heterozygosity excess) were detected in the shallow populations of Berlengas and Punta Mesco (Table 1; S3). All remaining populations showed no heterozygosity excess (on average) across loci (Table 1; S3).
<figure><img src="image_8.png" /><figcaption>**Figure 4** Pairwise genetic differentiation *F*<sub>ST</sub> (mean ± SD) within sampling site (i.e., shallow vs. deeper populations; light grey) and between sampling sites (dark grey). Asterisks indicate significant differences in genetic differentiation levels (*P* < 0.05).</figcaption></figure>

Full-size

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

DOI: 10.7717/peerj.6794/fig-4

DISCUSSION

The distribution of genetic diversity corroborated the expectations associating extreme
environmental changes with the reduction of effective population sizes and genetic diversity
levels (Garrabou et al., 2009). Our model species *P. clavata* consistently showed less allelic
richness, a smaller number of private alleles and bottleneck signs (for two sites) in shallow
waters populations. This pattern found across the Atlantic Ocean and Mediterranean Sea
was explained by the negative relationship between thermal stress and depth, and not
by the used proxies of essential resources (i.e., productivity, current velocity and slope).
In fact, the examination of daily time series of bottom temperatures for each site and
depth range (S4) shows shallow water masses recurrently surpassing the physiological
threshold of 24-25 °C from which mortality occurs (Previati et al., 2010; Pairaud et al.,
2014; Boavida et al., 2016b), particularly in the Mediterranean Sea. This is in line with
previous studies linking mass-mortality events of *P. clavata* with marine heatwaves, which
reported increased mortality in more shallow environments (Cerrano et al., 2005; Linares
et al., 2005). Still, the observed differences in genetic diversity may have resulted from
additional factors beyond temperature changes, not considered in our study. For instance,
the scarcity of favourable rocky bottoms may strongly reduce population sizes of *P. clavata*
(Linares et al., 2008; Boavida et al., 2016b). This could be a plausible explanation for the
one exception found in the diversity levels of deep populations. The lower diversity in the
deeper population of Petite Congloue may result from a smaller and patchier distribution
of individuals, if the availability of rocky bottoms is limited. In a particular survey, Gori
et al., 2011 found the density of *P. clavata* decreasing with depth, as populations were
highly dependent on the presence of rocky vertical walls in shallower waters. In our study,
the availability of rocky bottoms was not considered, although slope may be considered
as a proxy for such environmental factor at biogeographical scales (e.g., Boavida et al.,
2016b), it might not have reflected site scale habitat changes between our sampling depths.
Further, competition with other species, such as photosynthetic algae (Zabala & Ballesteros,
1989), together with higher hydrodynamics regimes and human induced disturbances (e.g.,
anchoring, scuba diving), may further help explaining the overall pattern found, restricting
large, continuous and richer populations to deeper sites (Harmelin & Marinopoulos, 1994; Cúrdia et al., 2013).

Contrarily to allelic richness and number of private alleles, gene diversity (i.e., expected heterozygosity) did not correlated with depth, nor with the additional predictors considered in regression analyses. Gene diversity is not as sensitive as allelic richness to detecting historical population changes (Leberg, 1992; Petit, El Mousadik & Pons, 1998; Spencer, Neigel & Leberg, 2000) since drift resulting from population size reductions theoretically affects more the rare alleles than the frequent ones (Nei, Maruyama & Chakraborty, 1975). Accordingly, the diversity pattern found in deeper populations implies persistence without significant population reduction, raising the hypothesis of deep refugia for *P. clavata*. These cryptic populations may play an important role in buffering the loss of diversity in shallow waters, as reported for other corals (Bongaerts et al., 2010; Bongaerts et al., 2017; Smith et al., 2014). Similar to terrestrial elevational refugia, the stability of deep environments may be an important mechanism safeguarding regional genetic diversity for the species as a whole (Epps et al., 2006; Assis et al., 2016; Lourenço et al., 2016).

While diversity levels of *P. clavata* followed the patterns of thermal stress, the bottleneck tests did not provide a strong empirical support of recent bottlenecks in shallow populations (exceptions for Berlengas e Punta Mesco). Although recent environmental changes may account for genetic diversity losses, the fluctuations of past climate changes may further shape the genetic structure of marine species, both latitudinally and along depth gradients (Maggs et al., 2008; Provan, 2013; Assis et al., 2014; Assis et al., 2016; Assis et al., 2018; Neiva et al., 2016). The lower diversity found in shallow populations could have resulted from past climate extremes, an historical effect that cannot be detected with tests based on heterozygosity excess (Peery et al., 2012). Also, one cannot discard the limited statistical power of these tests when using small sample sizes, as in our case (small number of individuals and loci; (Piry, Luikart & Cornuet, 1999; Peery et al., 2012). Our approach could have overlooked actual population declines. For instance, previous studies on *P. clavata*, with similar sample sizes, drew unambiguous conclusions by failing to detect bottlenecks in populations that actually suffered mass-mortality events (Pilczynska et al., 2016; Padrón et al., 2018).

Population differentiation was found higher between than within sampled sites. This is in line with the low dispersal potential of the species, posing highly structured gene pools throughout the Atlantic Ocean and the Mediterranean Sea (Linares et al., 2007; Mokhtar-Jamaï et al., 2011; Pilczynska et al., 2016; Padrón et al., 2018). At the site scale, the non-significant differentiation levels (Punta Mesco, Livorno, Pota Del Llop, Grotte Peres and Petit Congloue) suggest that gene flow and admixture occurs along depth ranges. While summer thermoclines were thought to be a major force structuring genetic diversity of gorgonian species (Costantini et al., 2011; Costantini et al., 2016), our results do not support such hypothesis. The shallow populations of *P. clavata* harbouring lower diversity levels are not in complete isolation from those in deep richer areas, as previously suggested for the species (Cerrano & Bavestrello, 2008). However, contrasting patterns of coral gene flow can occur across different environments (e.g., Van Oppen et al., 2011; Mokhtar-Jamaï et al., 2011; Pilczynska et al., 2016). For instance, the higher within site differentiation in
Sagres can result from the alongshore coastal circulation regime in the Gulf of Cadiz that extends to the Sagres region (Garel et al., 2016), promoting high connectivity levels along deep neighbouring populations, or from the fact that the shallow population of Sagres was located within a cave, possibly, further hindering connectivity. This was an exception though, as the overall connectivity pattern found is one of restricted connectivity between deeper populations. Thus, their diversity levels may have resulted from persistence alone, and not by admixture processes, as also suggested by the number of private alleles. The differentiation pattern further allows concluding that connectivity along depth ranges (within site) does not seem to completely offset the more detrimental effects occurring in shallow water populations, reducing genetic diversity of *P. clavata*, as observed.

Our findings of deep isolated richer populations have significant conservation value for the species as a whole. Future environmental changes are predicted to produce a major redistribution of marine biodiversity in the Atlantic Ocean and the Mediterranean Sea (Albouy et al., 2013; Assis, Araújo & Serrão, 2017), particularly in the business as usual climate scenario (RCP8.5; Assis, Araújo & Serrão, 2017; Assis et al., 2017). The reduction of biodiversity through the loss of rich and distinct gene pools found in deep populations may contribute to the loss of genes responsible for the species adaptation and evolution (Hughes & Stachowicz, 2004; Hampe & Petit, 2005; Reusch et al., 2005). The identified genetically rich populations of *P. clavata* represent an important baseline for Climate Change Integrated Conservation Strategies for phylogeographic lineages (Hannah, Midgley & Millar, 2002).

## CONCLUSIONS

This study demonstrates that deeper marine populations can consistently harbour higher genetic diversity than those in shallower environments. Theoretical regression modelling using environmental data showed that such genetic pattern may arise from thermal stress alone, reducing population sizes and diversity levels. These findings raised the hypothesis of deep refugia for *P. clavata*, with deeper populations in more stable environments safeguarding the species biodiversity across the Mediterranean and Atlantic populations. In fact, the analyses of temperature time-series showed shallow populations exposed to temperatures surpassing the species thermal tolerance, often for as long as over 2 months, while deep populations were almost never exposed to such conditions. The patterns of population differentiation revealed that gene flow occurs between shallow and deeper populations, although not in sufficient levels to homogenise depth differences in diversity levels. Empirical evidence is provided for deep persistent populations of gorgonians with the potential to safeguard richer and unique gene pools. While additional genetic data could better support the major findings, particularly by increasing the power of bottleneck tests (more molecular markers and sampling sites), this study represents a timely baseline to conserve populations with higher conservation value in risk of disappearing.
# ADDITIONAL INFORMATION AND DECLARATIONS

## Funding

This study was supported by financial support to CESAM (UID/AMB/50017/2019), to FCT/MCTES through national funds, and the co-funding by the FEDER, within the PT2020 Partnership Agreement and Compete 2020. This study was also supported by a Pew Marine Fellowship (USA), the National Geographic Channel through project Deep Reefs, a National Geographic/Waitt grant (no. W153-11), the InAqua Conservation Fund (Oceanário de Lisboa), the European Regional Development Fund (COMPETE program) and the Foundation for Science and Technology (FCT) of Portugal through postdoctoral fellowship SFRH/BPD/111003/2015 and programs CCMAR/Multi/04326/2013, UID/MAR/04292/2013 and PTDC/BIA-BIC/114526/2009 (DiverseShores—Testing associations between genetic and community diversity in European rocky shore environments). Joanna Pilczynska was supported by a MARES Grant. MARES is a Joint Doctorate programme selected under Erasmus Mundus coordinated by Ghent University (FPA 2011-0016). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## Grant Disclosures

The following grant information was disclosed by the authors:

CESAM: UID/AMB/50017/2019.

FCT/MCTES through national funds.

FEDER.

Pew Marine Fellowship (USA), the National Geographic Channel.

National Geographic/Waitt: W153-11.

InAqua Conservation Fund (Oceanário de Lisboa).

European Regional Development Fund (COMPETE program).

Foundation for Science and Technology (FCT) of Portugal through postdoctoral fellowship: SFRH/BPD/111003/2015, CCMAR/Multi/04326/2013, UID/MAR/04292/2013, PTDC/BIA-BIC/114526/2009.

DiverseShores—Testing associations.

Erasmus Mundus coordinated by Ghent University: FPA 2011-0016.

## Competing Interests

The authors declare there are no competing interests.

## Author Contributions

* Joanna Pilczynska and Joana Boavida conceived and designed the experiments, performed the experiments, analyzed the data, contributed reagents/materials/analysis tools, prepared figures and/or tables, authored or reviewed drafts of the paper, approved the final draft.

* Silvia Cocito, Ester A Serrão and Henrique Queiroga conceived and designed the experiments, analyzed the data, contributed reagents/materials/analysis tools, authored or reviewed drafts of the paper, approved the final draft.
* Jorge Assis performed the experiments, analyzed the data, contributed reagents/materials/analysis tools, prepared figures and/or tables, authored or reviewed drafts of the paper, approved the final draft.

* Eliza Fragkopoulou analyzed the data, contributed reagents/materials/analysis tools, prepared figures and/or tables, authored or reviewed drafts of the paper, approved the final draft.

## Field Study Permissions

The following information was supplied relating to field study approvals (i.e., approving body and any reference numbers):

1. Parco Nazionale delle Cinque Terre Area Marina Protetta delle Cinque Terre.

2. Instituto de Conservação da Natureza e da Biodiversidade

## Data Availability

The following information was supplied regarding data availability:

Microsatellite data for Paramuricea clavata samples are provided in Dataset S1. The data are available in GEN format and can be used in GENEPOP software package (http://genepop.curtin.edu.au). Each sample is named after the initials of the sampling location, where the last letter “S” stands for shallow populations and “D” stands for deep populations.

## Supplemental Information

Supplemental information for this article can be found online at http://dx.doi.org/10.7717/peerj.6794#supplemental-information.

## REFERENCES

* **Alberto F. 2009.** MsatAllele-1.0: an R package to visualize the binning of microsatellite alleles. *Journal of Heredity* **100**:394–397 DOI 10.1093/jhered/esn110.

* **Albouy C, Guilhaumon F, Leprieur F, Lasram FBR, Somot S, Aznar R, Velez L, Le Loch F, Mouillot D. 2013.** Projected climate change and the changing biogeography of coastal Mediterranean fishes. *Journal of Biogeography* **40**(3):534–547.

* **Araújo RM, Assis J, Aguillar R, Airoldi L, Bárbara I, Bartsch I, Bekkby T, Christie H, Davoult D, Derrien-Courtel S, Fernandez C, Fredriksen S, Gevaert F, Gundersen H, Le Gal A, Lévêque L, Mieszkowska N, Norderhaug KM, Oliveira P, Puente A, Rico JM, Rinde E, Schubert H, Strain EM, Valero M, Viard F, Sousa-Pinto I. 2016.** Status, trends and drivers of kelp forests in Europe: an expert assessment. *Biodiversity and Conservation* **25**:1319–1348 DOI 10.1007/s10531-016-1141-7.

* **Assis J, Álvares E, Nelson S, Coelho C, Tempera F, Valero M, Alberto F. 2018.** Past climate changes and strong oceanographic barriers structured low- latitude genetic relics for the golden kelp *Laminaria ochroleuca*. *Journal of Biogeography* **45**:2326–2336 DOI 10.1111/jbi.13425.

* **Assis J, Araújo BM, Serrão EA. 2017.** Projected climate changes threaten ancient refugia of kelp forests in the North Atlantic. *Global Change Biology* **24**:e55–e66.
**Assis J, Castilho Coelho N, Alberto F, Valero M, Raimondi P, Reed D, Alvares Serrão E. 2013.** High and distinct range-edge genetic diversity despite local bottlenecks. *PLOS ONE* 8(7):e68646 DOI 10.1371/journal.pone.0068646.

**Assis J, Coelho NC, Lamy T, Valero M, Alberto F, Serrão EA. 2016.** Deep reefs are climatic refugia for genetic diversity of marine forests. *Journal of Biogeography* 43:833–844 DOI 10.1111/jbi.12677.

**Assis J, Serrão EA, Claro B, Perrin C, Pearson GA. 2014.** Climate-driven range shifts explain the distribution of extant gene pools and predict future loss of unique lineages in a marine brown alga. *Molecular Ecology* 23:2797–2810 DOI 10.1111/mec.12772.

**Assis J, Tyberghein L, Bosch S, Verbruggen H, Serrão EA, De Clerck O. 2017.** BioORACLE v2.0: extending marine data layers for bioclimatic modelling. *Global Ecology and Biogeography* 27(3):277–284 DOI 10.1111/geb.12693.

**Bavestrello G, Bo M, Canese S, Sandulli R, Cattaneo-Vietti R. 2014.** The red coral populations of the gulfs of Naples and Salerno: human impact and deep mass mortalities. *Italian Journal of Zoology* 81:552–563 DOI 10.1080/11250003.2014.950349.

**Bo M, Canese S, Spaggiari C, Pusceddu A, Bertolino M, Angiolillo M, Giusti M, Loreto MF, Salvati E, Greco S, Bavestrello G. 2012.** Deep coral oases in the South Tyrrhenian Sea. *PLOS ONE* 7(11):e49870 DOI 10.1371/journal.pone.0049870.

**Boavida J, Paulo D, Aurelle D, Arnaud-Haound S, Marschal C, Gonçalves JMS, Serrão EA. 2016a.** A well-kept treasure at depth: precious red coral rediscovered in Atlantic deep coral gardens (SW Portugal) after 300 years. *PLOS ONE* 11(1):e0147228 DOI 10.1371/journal.pone.0147228.

**Boavida J, Silva I, Assis J, Serrão EA. 2016b.** Overlooked habitat of a vulnerable gorgonian revealed in the Atlantic and Mediterranean by niche modelling. *Scientific Reports* 6:36460 DOI 10.1038/srep36460.

**Bolker BM, Brooks ME, Clark CJ, Geange SW, Poulsen JR, Stevens MHH, White JSS. 2009.** Generalized linear mixed models: a practical guide for ecology and evolution. *Trends in Ecology and Evolution* 24:127–135 DOI 10.1016/j.tree.2008.10.008.

**Bongaerts P, Ridgway T, Sampayo EM, Hoegh-Guldberg O. 2010.** Assessing the deep reef refugia hypothesis: focus on Caribbean reefs. *Coral Reefs* 29:1–19 DOI 10.1007/s00338-009-0580-y.

**Bongaerts P, Riginos C, Brunner R, Englebert N, Smith SR, Hoegh-Guldberg O. 2017.** Deep reefs are not universal refuges: reseeding potential varies among coral species. *Science Advances* 3(2):e1602373 DOI 10.1126/sciadv.1602373.

**Busch JD, Waser PM, DeWoody JA. 2007.** Recent demographic bottlenecks are not accompanied by a genetic signature in banner-tailed kangaroo rats (Dipodomys spectabilis). *Molecular Ecology* 16:2450–2462 DOI 10.1111/j.1365-294X.2007.03283.x.

**Cerrano C, Arillo A, Azzini F, Calcinai B, Castellano L, Muti C, Valisano L, Zega G, Bavestrello G. 2005.** Gorgonian population recovery after a mass mortality event. *Aquatic Conservation* 15:147–157 DOI 10.1002/aqc.661.
**Cerrano C, Bavestrello G. 2008.** Medium-term effects of die-off of rocky benthos in the Ligurian SeaWhat can we learn from gorgonians? *Chemistry and Ecology* **24**:73–82
DOI 10.1080/02757540801979648.

**Coma R, Zabala M, Gili JM. 1995.** Sexual reproductive effort in the Mediterranean gorgonian *Paramuricea clavata*. *Marine Ecology Progress Series* **117**:185–192
DOI 10.3354/meps117185.

**Cornuet JM, Luikart G. 1996.** Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. *Genetics* **144**:2001–2004.

**Costantini F, Rossi S, Pintus E, Cerrano C, Gili JM, Abbiati M. 2011.** Low connectivity and declining genetic variability along a depth gradient in *Corallium rubrum* populations. *Coral Reefs* **30**:991–1003 DOI 10.1007/s00338-011-0771-1.

**Costantini F, Gori A, Lopez-González P, Bramanti L, Rossi S, Gili JM, Abbiati M. 2016.** Limited genetic connectivity between gorgonian morphotypes along a depth gradient. *PLOS ONE* **11**:1–20.

**Cúrdia J, Monteiro P, Afonso CML, Santos MN, Cunha MR, Gonçalves JMS. 2013.** Spatial and depth-associated distribution patterns of shallow gorgonians in the Algarve coast (Portugal, NE Atlantic). *Helgoland Marine Research* **67**:521–534
DOI 10.1007/s10152-012-0340-1.

**Diekmann OE, Serrão EA. 2012.** Range-edge genetic diversity: locally poor extant south southern patches maintain a regionally diverse hotspot in the seagrass *Zostera marina*. *Molecular Ecology* **21**:1647–1657 DOI 10.1111/j.1365-294X.2012.05500.x.

**Dormann CF, McPherson JM, Araújo MB, Bivand R, Bolliger J, Carl G Davies, RG, Hirzel A, Jetz W, Kissling DW, Kühn I, Ohlemüller R, Peres-Neto PR, Reineking B, Schröder B, Schurr FM, Wilson R. 2007.** Methods to account for spatial autocorrelation in the analysis of species distributional data: a review. *Ecography* **30**:609–628
DOI 10.1111/j.2007.0906-7590.05171.x.

**Eckert CG, Samis KE, Lougheed SC. 2008.** Genetic variation across species’ geographical ranges: the central-marginal hypothesis and beyond. *Molecular Ecology* **17**:1170–1188 DOI 10.1111/j.1365-294X.2007.03659.x.

**Edwards LJ, Muller KE, Wolfinger RD, Qaqish BF, Schabenberger O. 2008.** An R2 statistic for fixed effects in the linear mixed model. *Statistics in Medicine* **27**(29):6137–6157 DOI 10.1002/sim.3429.

**Elith J, Leathwick JR, Hastie T. 2008.** A working guide to boosted regression trees. *Journal of Animal Ecology* **77**:802–813 DOI 10.1111/j.1365-2656.2008.01390.x.

**Epps CW, Palsboll PJ, Wehausen JD, Roderick GK, McCullough DR. 2006.** Elevation and connectivity define genetic refugia for mountain sheep as climate warms. *Molecular Ecology* **15**:4295–4302 DOI 10.1111/j.1365-294X.2006.03103.x.

**Frade PR, Bongaerts P, Englebert N, Rogers A, Gonzalez-Rivero M, Hoegh-Guldberg O. 2018.** Deep reefs of the Great Barrier Reef offer limited thermal refuge during mass coral bleaching. *Nature Communications* **9**(1):3447
DOI 10.1038/s41467-018-05741-0.
**Garel E, Laiz I, Drago T, Relvas P. 2016.** Characterisation of coastal counter-currents on the inner shelf of the Gulf of Cadiz. *Journal of Marine Systems* **155**:19–34 DOI 10.1016/j.jmarsys.2015.11.001.

**Garrabou J, Ballesteros E, Zabala M. 2002.** Structure and dynamics of north-western Mediterranean rocky benthic communities along a depth gradient. *Estuarine, Coastal and Shelf Science* **55**:493–508 DOI 10.1006/ecss.2001.0920.

**Garrabou J, Coma R, Bensoussan N, Chevaldonné P, Cigliano M, Diaz D, Harmelin LG, Gambi MC, Kerstingk DK, Ledoux JB, Lejeusne C, Linares C, Marschal C, Pérez T, Ribes M, Romano JC, Serrano E, Teixido N, Torrents O, Zabala M, Zuberer F, Ferrano C. 2009.** Mass mortality in NW Mediterranean rocky benthic communities: effects of the 2003 heat wave. *Global Change Biology* **15**:1090–1103 DOI 10.1111/j.1365-2486.2008.01823.x.

**Giordano AR, Ridenhour BJ, Storfer A. 2007.** The influence of altitude and topography on genetic structure in the long-toed salamander (*Ambystoma macrodactylym*). *Molecular Ecology* **16**:1625–1637 DOI 10.1111/j.1365-294X.2006.03223.x.

**Hampe A, Petit RJ. 2005.** Conserving biodiversity under climate change: the rear edge matters. *Ecology Letters* **8**:461–467 DOI 10.1111/j.1461-0248.2005.00739.x.

**Hannah L, Midgley GF, Millar D. 2002.** Climate change-integrated conservation strategies. *Global Ecology and Biogeography* **11**:485–495 DOI 10.1046/j.1466-822X.2002.00306.x.

**Harmelin J-G, Marinopoulos J. 1994.** Population structure and partial mortality of the gorgonian *Paramuricea clavata* (Risso) in the north-western Mediterranean (France, Port-Cros Island). *Marine Life* **4**:5–13.

**Holstein DM, Smith TB, Gyory J, Paris CB. 2015.** Fertile fathoms: deep reproductive refugia for threatened shallow corals. *Scientific Reports* **5**:12407 DOI 10.1038/srep12407.

**Hughes AR, Stachowicz JJ. 2004.** Genetic diversity enhances the resistance of a seagrass ecosystem to disturbance. *Proceedings of the National Academy of Sciences of the United States of America* **101**(24):8998–9002 DOI 10.1073/pnas.0402642101.

**Johannesson K, André C. 2006.** Life on the margin: genetic isolation and diversity loss in a peripheral marine ecosystem, the Baltic Sea. *Molecular Ecology* **15**:2013–2029 DOI 10.1111/j.1365-294X.2006.02919.x.

**Leberg PL. 1992.** Effects of population bottlenecks on genetic diversity as measured by allozyme electrophoresis. *Evolution* **46**:477–494 DOI 10.1111/j.1558-5646.1992.tb02053.x.

**Lima FP, Queiroz N, Ribeiro PA, Hawkins SJ, Santos AM. 2006.** Recent changes in the distribution of a marine gastropod *Patella rustica* Linnaeus, 1758, and their relationship to unusual climatic events. *Journal of Biogeography* **33**:812–822 DOI 10.1111/j.1365-2699.2006.01457.x.

**Linares C, Coma R, Diaz D, Zabala M, Hereu B, Dantart L. 2005.** Immediate and delayed effects of a mass mortality event on gorgonian population dynamics and benthic community structure in the NW Mediterranean Sea. *Marine Ecology Progress Series* **305**:127–137 DOI 10.3354/meps305127.
**Linares C, Coma R, Garrabou J, Díaz D, Zabala M. 2008.** Size distribution, density and disturbance in two Mediterranean gorgonians: *Paramuricea clavata* and *Eunicella singularis*. *Journal of Applied Ecology* **45**:688–699.

**Linares C, Doak DF, Coma R, Díaz D, Zabala M. 2007.** Life history and viability of a long-lived marine invertebrate: the octocoral *Paramuricea clavata*. *Ecology* **88**:918–928 DOI 10.1890/05-1931.

**Locke M, Baack E, Toonen R. 2000.** STRand. Version 2.2.30. Veterinary Genetics Lab. Davis: University of California.

**Lourenço CR, Zardi GI, McQuaid CD, Serrão EA, Pearson GA, Jacinto R, Nicastro KR. 2016.** Upwelling areas as climate change refugia for the distribution and genetic diversity of a marine macroalga. *Journal of Biogeography* **43**:1595–1607 DOI 10.1111/jbi.12744.

**Ludwig M, Schlinkert H, Holzschuh A, Fischer C, Scherber C, Trnka A, Tscharntke T, Batáry P. 2012.** Landscape-moderated bird nest predation in hedges and forest edges. *Acta Oecologica* **45**:50–56 DOI 10.1016/j.actao.2012.08.008.

**Luikart G, Cornuet JM. 1998.** Empirical evaluation of a test for identifying recently bot-tlencked populations from allele frequency data. *Conservation Biology* **12**:228–237 DOI 10.1046/j.1523-1739.1998.96388.x.

**Maggs CA, Castilho R, Foltz D, Henzler C, Jolly MT, Kelly J, Olsen J, Perez K, Stam W, Vainola R, Viard F, Wares J. 2008.** Evaluating signatures of glacial refugia for North Atlantic Marine Organisms. *Ecology* **89**:S108–S122 DOI 10.1890/08-0257.1.

**Meirmans PG, Van Tienderen PH. 2004.** GENOTYPE and GENODIVE: two programs for the analysis of genetic diversity of asexual organisms. *Molecular Ecology Notes* **4**:792–794 DOI 10.1111/j.1471-8286.2004.00770.x.

**Mokhtar-Jamaï K, Pascual M, Jb Ledoux, Coma R, Féral JP, Garrabou J, Aurelle D. 2011.** From global to local genetic structuring in the red gorgonian *Paramuricea clavata*: the interplay between oceanographic conditions and limited larval dispersal. *Molecular Ecology* **20**:3291–3305 DOI 10.1111/j.1365-294X.2011.05176.x.

**Nei M, Maruyama T, Chakraborty R. 1975.** The bottleneck effect and genetic variability in populations. *Evolution* **29**:1–10 DOI 10.1111/j.1558-5646.1975.tb00807.x.

**Neiva J, Serrão E, Assis J, Pearson G, Coyer J, Olsen J, Valero M et al.. 2016.** Climate oscillations, range shifts and phylogeographic pat-terns of North Atlantic Fucaceae. In: Hu Z-M, Hu & Fraser C, eds. *Sea-weed phylogeography*. Dordrecht: Springer, 279–308.

**Ohsawa T, Ide Y. 2008.** Global patterns of genetic variation in plant species along vertical and horizontal gradients on mountains. *Global Ecology and Biogeography* **17**:152–163 DOI 10.1111/j.1466-8238.2007.00357.x.

**Padrón M, Costantini F, Bramanti L, Guizien K, Abbiati M. 2018.** Genetic connectivity supports recovery of gorgonian populations affected by climate change. *Aquatic Conservation: Marine and Freshwater Ecosystems* **28**:776–787 DOI 10.1002/aqc.2912.

**Pairaud IL, Bensoussan N, Garreau P, Faure V, Garrabou J. 2014.** Impacts of climate change on coastal benthic ecosystems: assessing the current risk of mortality
outbreaks associated with thermal stress in NW Mediterranean coastal areas. *Ocean Dynamics* **64**:103–115 DOI 10.1007/s10236-013-0661-x.

Peery ZM, Kirby R, Reid BN, Stoelting R, Doucet-Bëer E, Robinson S, Vásquez-Carrillo C, Pauli JN, Palsboll PJ. 2012. Reliability of genetic bottleneck tests for detecting recent population declines. *Molecular Ecology* **21**:3403–3418 DOI 10.1111/j.1365-294X.2012.05635.x.

Perez T, Garrabou J, Sartoretto S, Harmelin JG, Francour P, Vacelet J. 2000. Mass mortality of marine invertebrates: an unprecedented event in the Northwestern Mediterranean. *Comptes Rendus de l' Académie des Sciences Paris* **323**:853–865.

Petit RJ, El Mousadik A, Pons O. 1998. Identifying populations for conservation on the basis of genetic markers. *Conservation Biology* **12**:844–855 DOI 10.1046/j.1523-1739.1998.96489.x.

Pilczynska J, Boavida J, Cocito S, Lombardi C, Peirano A, Queiroga H. 2017. Low clonal propagation in Atlantic and Mediterranean populations of the red gorgonian *Paramuricea clavata* (Octocorallia). *Scientia Marina* **81**:103–110 DOI 10.3989/scimar.04473.05A.

Pilczynska J, Cocito S, Boavida J, Serrão E, Queiroga H. 2016. Genetic diversity and local connectivity in the mediterranean red gorgonian coral after mass mortality events. *PLOS ONE* **11**(3):e0150590 DOI 10.1371/journal.pone.0150590.

Piry S, Luikart G, Cornuet JM. 1999. BOTTLENECK: a computer program for detecting recent reductions in the effective population size using allele frequency data. *Journal of Heredity* **90**:502–503 DOI 10.1093/jhered/90.4.502.

Previati M, Scinto A, Cerrano C, Osinga R. 2010. Oxygen consumption in Mediterranean octocorals under different temperatures. *Journal of Experimental Marine Biology and Ecology* **390**:39–48 DOI 10.1016/j.jembe.2010.04.025.

Provan J. 2013. The effects of past, present and future climate changeon range-wide genetic diversity in northern North Atlantic marinespecies. *Frontiers of Biogeography* **5**:60–66.

Provan J, Bennett KD. 2008. Phylogeographic insights into cryptic glacial refugia. *Trends in Ecology and Evolution* **23**:564–571 DOI 10.1016/j.tree.2008.06.010.

R Core Team. 2018. R: a language and environment for statistical computing. Version 3.5.2. Vienna: R Foundation for Statistical Computing. Available at https://www.R-project.org/.

Relvas P, Barton ED, Dubert J, Oliveira PB, Peliz A, Da Silva JCB, Santos AMP. 2007. Physical oceanography of the western Iberia ecosystem: latest views and challenges. *Progress Oceanography* **74**:149–173 DOI 10.1016/j-pocean.2007.04.021.

Reusch TB, Ehlers A, Hämmerli A, Worm B. 2005. Ecosystem recovery after climatic extremes enhanced by genotypic diversity. *Proceedings of the National Academy of Sciences of the United States of America* **102**:2826–2831 DOI 10.1073/pnas.0500008102.

Romano JC, Bensoussan N, Younes WA, Arlhac D. 2000. Anomalie thermique dans les eaux du golfe de Marseille durant l'été 1999. Une explication partielle de la mortalité d'invertébrés fixés? *Comptes Rendus de l'Académie des Sciences-Series III-Sciences de la Vie* **323**:415–427 DOI 10.1016/S0764-4469(00)00141-4.
**Salomidi M, Smith C, Katsanevakis S, Panayotidis P, Papathanassiou V. 2009.** Some observations on the structure and distribution of Gorgonian assemblages in the Eastern Mediterranean Sea. In: *Proceedings of the 1st symposium on conservation of the coralligenous bio-concretions (Tabarka 16–19 January)*. 242–245.

**Smith TB, Glynn PW, Maté JL, Toth LT, Gyory J. 2014.** A depth refugium from catastrophic coral bleaching prevents regional extinction. *Ecology* **95**:1663–1673
DOI 10.1890/13-0468.1.

**Spencer CC, Neigel JE, Leberg PL. 2000.** Experimental evaluation of the usefulness of microsatellite DNA for detecting demographic bottlenecks. *Molecular Ecology* **9**:1517–1528 DOI 10.1046/j.1365-294x.2000.01031.x.

**Van Oppen MJH, Bongaerts P, Underwood JN, Peplow LM, Cooper TF. 2011.** The role of deep reefs in shallow reef recovery: an assessment of vertical connectivity in a brooding coral from west and east Australia. *Molecular Ecology* **20**:1647–1660 DOI 10.1111/j.1365-294X.2011.05050.x.

**Winnepenninckx B, Backeljau T, De Wachter R. 1993.** Extraction of high molecular weight DNA from molluscs. *Trends in Genetics* **9**:407
DOI 10.1016/0168-9525(93)90102-N.

**Zabala M, Ballesteros E. 1989.** Surface-dependent strategies and energy flux in benthic marine communities or, why corals do not exist in the Mediterranean. *Scientia Marina* **53**:3–17.

**Zhang G, Chen JJ. 2013.** Adaptive Fitting of Linear Mixed-Effects Models with Correlated Random-effects. *Journal of Statistical Computation and Simulation* **83**:2291–2314 DOI 10.1080/00949655.2012.690763.

**Zuur AF, Ieno EN, Smith GM. 2007.** *Analyzing ecological data*. New York: Springer.