Demo-Genetic Approach for the Conservation and Restoration of a Habitat-Forming Octocoral: The Case of Red Coral, Corallium rubrum, in the Réserve Naturelle de Scandola

Gazulla, Carlota R.; López-Sendino, Paula; Antunes, Agostinho; Aurelle, Didier; Montero-Serra, Ignasi; Dominici, Jean-Marie; Linares, Cristina; Garrabou, Joaquim; Ledoux, Jean-Baptiste

doi:10.3389/fmars.2021.633057

ORIGINAL RESEARCH article

Front. Mar. Sci., 17 June 2021

Sec. Marine Ecosystem Ecology

Volume 8 - 2021 | https://doi.org/10.3389/fmars.2021.633057

This article is part of the Research Topic Biogenic Reefs at Risk: Facing Globally Widespread Local Threats and Their Interaction with Climate Change View all 12 articles

Demo-Genetic Approach for the Conservation and Restoration of a Habitat-Forming Octocoral: The Case of Red Coral, Corallium rubrum, in the Réserve Naturelle de Scandola

$\r\nCarlota R. Gazulla,$ Carlota R. Gazulla^1,2

Paula López-Sendino¹

Agostinho Antunes^3,4

Didier Aurelle^5,6

Ignasi Montero-Serra⁷

Jean-Marie Dominici⁸

Cristina Linares⁷

Joaquim Garrabou^1,5*

Jean-Baptiste Ledoux^1,3*

¹Institut de Ciències del Mar, Consejo Superior de Investigaciones Científicas, Barcelona, Spain
²Departament de Genètica i Microbiologia, Universitat Autònoma de Barcelona, Barcelona, Spain
³CIIMAR/CIMAR, Centro Interdisciplinar de Investigação Marinha e Ambiental, Universidade do Porto, Porto, Portugal
⁴Department of Biology, Faculty of Sciences, University of Porto, Porto, Portugal
⁵Aix Marseille Univ, Université de Toulon, CNRS, IRD, MIO, Marseille, France
⁶CNRS, EPHE, Institut de Systématique, Evolution, Biodiversité, Muséum National d’Histoire Naturelle, Sorbonne Université, Paris, France
⁷Departament de Biologia Evolutiva, Ecologia i Ciències Ambientals, Institut de Recerca de la Biodiversitat, Universitat de Barcelona, Barcelona, Spain
⁸Réserve Naturelle de Scandola, Parc Régional de Corse, Galeria, France

Marine protected areas (MPAs) are one of the most efficient conservation tools to buffer marine biodiversity loss induced by human activities. Beside effective enforcement, an accurate understanding of the eco-evolutionary processes underlying the patterns of biodiversity is needed to reap the benefits of management policies. In this context, integrating population genetics with demographic data, the demo-genetic approach, is particularly relevant to shift from a “species-based pattern” toward an “eco-evolutionary-based processes” conservation. Here, targeting a key species in the Mediterranean coralligenous, the red coral, Corallium rubrum, in an emblematic Mediterranean MPA, the “Réserve Naturelle de Scandola” (France), we applied demo-genetic approaches at two contrasted spatial scales, among populations and within one population, to (i) infer the demographic connectivity among populations in the metapopulation network and (ii) shed new light on the genetic connectivity and on the demographic transitions underlying the dynamics of a near-pristine population. Integrating different spatial and temporal scales, we demonstrated (i) an apparent temporal stability in the pattern of genetic diversity and structure in the MPA in spite of a dramatic demographic decline and (ii) contrasted levels of genetic isolation but substantial demographic connectivity among populations. Focusing on the near-pristine population, we complemented the characterization of red coral demographic connectivity suggesting (iii) temporal variability and (iv) the occurrence of collective dispersal. In addition, we demonstrated (v) contrasted patterns of spatial genetic structure (SGS), depending on the considered stage-class (adults vs. juveniles), in the near-pristine population. This last result points out that the overall SGS resulted from a restricted dispersal of locally produced juveniles (SGS among adults and juveniles) combined to mortality during early life stages (decrease of SGS from juveniles to adults). Demonstrating the occurrence of two management units and the importance of two populations (CAVB and ALE) for the network of connectivity, we made recommendations for the management of the Réserve Naturelle de Scandola. Besides, we contributed to the implementation of scientifically driven restoration protocols in red coral by providing estimates for the size, density, and distances among patches of transplanted colonies.

Introduction

Global change is drastically modifying the networks of interactions among the different levels of marine biodiversity, from genes to ecosystems, impacting ecosystem functioning and related socioeconomic services. Marine protected areas (MPAs) are one of the most efficient tools to buffer this anthropogenic biodiversity loss. Besides an effective enforcement (Edgar et al., 2014; Costello and Ballantine, 2015), the benefits of MPAs rely on an accurate understanding of the eco-evolutionary processes shaping and maintaining the patterns of biodiversity (Sale et al., 2005). Connectivity, defined as the exchange of individuals among populations (demographic connectivity) and their successful reproduction (genetic connectivity) (Lowe and Allendorf, 2010), is of critical importance in the functioning of MPAs. It underlies the dynamics and genetics of populations (Palumbi, 2004; Gagnaire et al., 2015). However, to date, design and management of MPAs usually rely on “species-based pattern” rather than “eco-evolutionary-based processes” considerations, restraining potentially their benefits.

One way to achieve the required shift, from a pattern- to a process-focused management, relies on the integration of population genetics into conservation strategies. While this is acknowledged since the beginning of conservation biology (Soulé, 1985), population genetics and genetic diversity remain overlooked by biodiversity managers (Cook and Sgrò, 2018; Laikre et al., 2020). This is particularly detrimental considering the recent developments in demo-genetic approaches that, combining demographic information with population genetics, allow the inference of eco-evolutionary processes acting over contemporary timescale. At large spatial scale (i.e., among populations), demo-genetic approaches allow, for instance, for the estimation of contemporary connectivity (see Waples and Gaggiotti, 2006; Broquet and Petit, 2009). Populations can thus be classified in terms of their role in the functioning of the network (e.g., source vs. sink) helping, accordingly, to prioritize conservation efforts. Beside their inputs in conservation, demo-genetic approaches can be used to improve restoration practices, especially when applied at fine spatial scale (i.e., among individuals within population). In marine sessile invertebrates, the combination between the dispersive larval phase and complex reproductive strategies (Bierne et al., 2016) can limit our understanding of ecological processes driving the population dynamics. Characterizing the fine-scale SGS within population, the non-random spatial distribution of genotypes through the formation of local pedigree structure (Vekemans and Hardy, 2004) allows to infer demographic parameters such as the “neighborhood size” (Nb) and the “mean parent-offspring distance” (σ_g) (see Rousset, 1997, Rousset, 2000). These two parameters provide insight into the distance range at which genetic interactions occurred among individuals within population (Vekemans and Hardy, 2004). Interestingly, SGS is a dynamic pattern, which can increase (e.g., Pardini and Hamrick, 2008) or decrease (e.g., Hampe et al., 2010) from juveniles to adults. Demographic processes such as density-dependent or -independent mortality during early-life stages may decrease the strength of SGS from propagule to recruit to adult-stage classes (e.g., Chung et al., 2003, 2007). Accordingly, the decomposition of SGS among stage classes can shed new light on the population dynamics of the targeted species by characterizing the interactions between the demographic and evolutionary processes shaping the pool of reproductive adults (Kalisz et al., 2001; Jacquemyn et al., 2006; Hampe et al., 2010). The characterization of the dynamics of the SGS within population and the estimation of related demographic parameters (Nb and σ_g) are thus of primary importance to support the definition of scientifically driven transplantation protocols (e.g., size and density of restored patches) and to set the size of restoration actions (e.g., distance among restored patches).

The Mediterranean Sea is a striking example of the challenges induced by global change on marine biodiversity. This biodiversity hotspot is submitted to direct anthropogenic pressures from overfishing to pollution and was recently recognized as a hotspot of climate change (Diffenbaugh et al., 2007; Cramer et al., 2018). In the last decades, mass mortality events (MMEs) linked to marine heat waves have been observed along thousands of kilometers of coastal habitats, impacting more than 90 species (Garrabou et al., 2019). Coralligenous biogenic reefs, which are among the richest Mediterranean communities (Ballesteros, 2006) and were recognized as priority habitats by the UNEP-MAP, have been particularly sensitive to MMEs (Garrabou et al., 2009, 2019). Populations of long-lived and habitat-forming macrobenthic species such as sponges, octocorals, or bryozoans were dramatically impacted by MMEs, with up to 80% of impacted individuals in some locations (Cerrano et al., 2000; Garrabou et al., 2001, 2009, 2019). The reported high incidence and mortality rates in these key species may induce cascading effects at the ecosystem level, questioning the future of coralligenous communities (Gómez-Gras et al., 2021).

Here, focusing on red coral, Corallium rubrum, a key species of the Mediterranean coralligenous, in an emblematic Mediterranean MPA, the “Réserve Naturelle de Scandola” (Parc Naturel Regional de Corse, France), our main objective is to provide insight into red coral connectivity among and within populations in order to refine the management of the MPA and to contribute to the definition of scientifically driven restoration protocols. The relevance of this case study relies on two main issues. First, the Scandola Marine Reserve is one of the most emblematic MPAs in the Mediterranean. Characterized by a well-conserved biodiversity, including coralligenous habitats, and a unique and near-pristine population of red coral (Cave-b; Garrabou et al., 2017), the MPA is a reference for managers and scientists. However, little is known regarding the eco-evolutionary processes underlying the functioning of the MPA, and particularly the pattern of connectivity among populations. Then, following almost three decades of multidisciplinary studies (Abbiati et al., 1993; Santangelo et al., 1993), red coral has been established as a model species in conservation biology. This habitat-forming octocoral, with low population dynamics and late sexual maturity (Garrabou and Harmelin, 2002; Marschal et al., 2004; Torrents et al., 2005), is found in contrasted habitats mainly in the Western Mediterranean and neighboring Atlantic (Zibrowius et al., 1984; Boavida et al., 2016), from 5- to more than 1000-m depth (Knittweis et al., 2016). Warming-induced mortalities (Garrabou et al., 2001) and overharvesting for its use in jewelry (Bruckner, 2009) induced a dramatic shift in the demographic structure of shallow populations, questioning the evolutionary trajectory of the species (Montero-Serra et al., 2019). Patterns of genetic structure and diversity in shallow red coral populations have been characterized from global, to regional (including the Scandola Marine Reserve; see below), to local scales. Overall, the species shows a significant genetic structure among populations at a scale of tens of meters (Costantini et al., 2007a), which increases with geographic distance (i.e., isolation by distance pattern, IBD; Ledoux et al., 2010a; Aurelle et al., 2011; Aurelle and Ledoux, 2013). While the IBD suggests a limited gene flow occurring mainly among close-by populations, it is noteworthy that the connectivity of the species remains poorly characterized. At regional scale, some heterogeneity in the patterns of structure has been revealed (Ledoux et al., 2010a). Analyzing populations from the Scandola Marine Reserve, Ledoux et al. (2010a) revealed a lack of IBD combined to a relatively low level of genetic diversity compared to other regions. At local scale, the significant SGSs reported among individuals within an unprotected population (Ledoux et al., 2010b) suggest that red coral populations may be considered as complex networks of genetically related individuals interacting locally (few centimeters). While these data improved our abilities to restore declining populations, the eco-evolutionary processes, such as genetic connectivity and density-dependent or -independent mortality underlying the establishment of the reproductive pool in red coral populations (Montero-Serra et al., 2015, 2018), remain to be characterized in order to provide scientifically driven restoration protocols.

Based on a sampling at the population and individual levels, we applied complementary demo-genetic approaches at two contrasting spatial scales: among populations and among individuals within one population. In a first step, we characterized the patterns of genetic diversity and structure among six red coral populations (i) to infer the demographic connectivity among these populations. In a second step, we used a geo-referenced sampling considering two stage-classes, adults vs. juveniles, (ii) to decompose the SGS in the Cave-b population, and (iii) to further characterize the potential impact of genetic connectivity and early life-stage mortality on the population dynamics of this near-pristine population. Our results suggested the occurrence of collective dispersal, improving our current knowledge regarding the ecology of C. rubrum. Besides, we used the genetic structure and demographic connectivity patterns to prioritize management efforts in the Scandola Marine Reserve. We translated the SGS results into restoration advices. In particular, we provided estimates for the size, density, and distance among patches to restore red coral populations.

Materials and Methods

Sampling and Spatial Data

Established in 1975 and listed as a world heritage site in 1983 by UNESCO, the Scandola Marine Reserve (Parc Naturel Regional de Corse, France) is particularly well enforced, belonging to the 0.23% of effectively protected Mediterranean surface (Claudet et al., 2020). In this MPA, individuals of C. rubrum from six populations were collected by SCUBA diving at depths from 15 to 25 m between 2015 and 2017 (Figure 1A). In five populations separated by distance ranging from 20 m to less than 2 km (DHM, ALI, ALE, PPL, and PLU), apical fragments (∼5 mm) from 30 to 51 mature individuals (i.e., colonies higher than 30 mm; Torrents et al., 2005) were randomly sampled. In Cave-b (CAVB), sampling was conducted following a 6.25-m-length horizontal transect at a depth of 23 m. Based on its demographic structure showing the highest biomass reported to date in the Mediterranean Sea (more than 100-fold) and the coexistence of a high number of juveniles, low size, and large centennial colonies in high density (201 colonies^∗m^–2), the Cave-b population is a unique and near-pristine population of red coral (Garrabou et al., 2017). The high density of this population combined to the morphological characteristics of the habitat prevents a georeferenced sampling at the individual level. Accordingly, 11 quadrats (25 × 25 cm) were sampled. The size of the quadrat was of the same order as the mean square parent–offspring distance (σ_g) estimated by Ledoux et al. (2010b). Nine quadrats were sampled every 50 cm along the transect. Two quadrats at the center of the transect were duplicated 25 cm above the transect. Within each quadrat, 15–20 adult colonies and a similar number of juvenile colonies were randomly collected (see Figure 1B and Table 1). In the red coral, sexual maturity was reported in colonies higher than 30 mm (Torrents et al., 2005). To avoid potential mis-assignment in the reproductive status of the samples, we considered adult and juvenile colonies as colonies higher than 100 mm and between 5 and 10 mm, respectively. For the adult colonies, small fragments (<10 mm) were sampled while, for the juveniles, the whole colonies were sampled.

FIGURE 1

Figure 1. (A) The « Réserve marine de Scandola » (Parc Naturel Régional de Corse, France) is located in the western coast of Corsica. Locations of the different populations sampled in the MPA are shown. The exact locations of CAVB and DHM are not shown on purpose. (B) Individual sampling conducted within the CAVB population. Adult and juvenile colonies were sampled within nine quadrats (25*25 cm) separated by 50 cm along a 6.25-m transect. Two quadrats were duplicated at the center of the transect, 25 cm under the main transect, resulting in 11 quadrats in total.

TABLE 1

Table 1. Location and genetic diversity estimators of the Corallium rubrum populations.

The resulting 634 apical fragments were conserved in 95% ethanol and stored at −20°C prior to DNA extraction.

DNA Extraction, Microsatellite Genotyping, and Quality Control

DNA extraction, microsatellite genotyping, and quality control are described in Supplementary Material. Briefly, all the individuals were genotyped using eight microsatellite loci amplified in three multiplexes: Multiplex 1 (Mic13/Mic20/Mic26), Multiplex 2 (Mic24/Mic25/Mic27), and Multiplex 3 (Mic22/Mic23) (Ledoux et al., 2010a). Following the quality control, statistical analyses were conducted with a dataset including 580 colonies, corresponding to the same number of unique multilocus genotypes, including 201 adults (CAVB-AD) and 166 juveniles (CAVB-JUV) from Cave-b and 213 individuals from five populations (Table 1).

Hardy–Weinberg Equilibrium and Genetic Diversity

Probability of identity (Pi; the average probability that two samples will share the same genotype), Hardy–Weinberg equilibrium, and genetic diversity analyses for each population, including total number of alleles, observed (H_o) and unbiased heterozygosity (H_e; Nei, 1973), f estimator of F_IS, and rarefied allelic richness (Ar₍₅₆₎; Petit et al., 1998), are described in Supplementary Material.

SGS Among Populations

A clustering analysis using STRUCTURE 2.2 (Pritchard et al., 2000) was conducted to evaluate the number of genetic clusters (K) from the individuals’ genotypes without assumptions on population boundaries and considering the whole dataset. Because of the unbalanced sampling, we considered the admixture model with correlated allele frequencies (Falush et al., 2003, 2007) and population-specific ancestry prior (α) with an α = 0.17 (=1/number of assumed populations; Wang, 2017). Ten independent runs were performed for each K with a burn-in period of 500,000 followed by 250,000 iterations. The K-value corresponding to “upper most hierarchical level of structure” was determined using Ln Pr(X| K) (Pritchard et al., 2000), the ΔK method (Evanno et al., 2005), and the MedMedK, MedMeaK, MaxMedK, and MaxMeaK statistics with three different thresholds (0.6, 0.7, and 0.8; Puechmaille, 2016), CLUMPP 1.1 (Jakobsson and Rosenberg, 2007), and DISTRUCT 1.1 (Rosenberg, 2003) were used for graphical outputs.

We conducted a discriminant analysis of principal components (DAPC, Jombart et al., 2010) in ADEGENET (Jombart, 2008). Contrary to STRUCTURE, DAPC describes patterns of diversity without assumption about the underlying population genetic model (Jombart et al., 2010). Data were transformed into principal components, and discriminant analyses were used to maximize variation among groups while minimizing variations within the group. We used the seven population locations (discriminating here among adults and juveniles from CAVB) as group prior. Based on the a-score method, the number of principal components was set to 46, while we retained three discriminant functions.

Global and pairwise differentiations were quantified using Weir and Cockerham’s (1984) estimators of F_ST in GENEPOP (Rousset, 2008). Genotypic differentiation was tested using an exact test based on the MC algorithm (Guo and Thompson, 1992) with default parameters in GENEPOP.

Isolation by distance among populations was tested using a linear regression of genetic distances computed as F_ST/(1 − F_ST) over logarithms of the geographic distances (ln(d)) (Rousset, 1997). The significance of the correlation between the genetic distances (F_ST/(1 − F_ST)) and the logarithms of geographic distances (Ln(d)) was tested by the Mantel test with 10,000 permutations in GENEPOP.

We used GESTE (Foll and Gaggiotti, 2006) to compute the population-specific F_ST. Measuring the genetic differentiation of each population, this method gives insight into the relative impact of genetic drift on the differentiation of the considered population (Gaggiotti and Foll, 2010).

Assignment Tests

Two sets of assignment tests in GeneClass2 (Piry et al., 2004) were conducted in order to (i) estimate the demographic connectivity among the six populations (not considering juveniles from CAVB) and (ii) detect migrants within the juveniles sampled in CAVB. First, we conducted a filtered assignment analysis following Lukoschek et al. (2016). We identified first-generation migrants (FGMs) in each population using the Bayesian criteria of Rannala and Mountain (1997) with 100,000 simulated genotypes and an alpha of 0.005. FGMs were removed from the dataset and reassigned to the reference dataset (i.e., dataset without FGMs). Migrants were assigned to a population if their assignment probability was >0.01. We allowed for multiple assignments (i.e., one individual originating potentially from different populations). When the assignment probability of an FGM was lower than 0.01 for all populations, we considered this FGM to come from an unsampled population.

Then, focusing on CAVB, we assigned CAVB-JUV considering the remaining samples (i.e., five populations + CAVB-AD) as reference. We used the Bayesian assignment method of Rannala and Mountain (1997) simulating 100,000 individuals and an alpha of 0.005. Following this analysis, CAVB-JUV was divided in two samples, CAVB-JUV-L and CAVB-JUV-M, distinguishing the juveniles originating from CAVB (locally produced) from the immigrant juveniles, respectively.

Fine-Scale SGS and Neighborhood Size

Focusing on CAVB, we analyzed the SGS by regressing the Nason’s kinship coefficients (F_ij; Loiselle et al., 1995) among pairs of individuals on the logarithms of pairwise geographic distances, in SPAGEDI (Hardy and Vekemans, 2002). We first conducted this analysis considering a “global dataset” including CAVB-AD and CAVB-JUV-L. Then, we decomposed the SGS by conducting the same analyses considering three categories: CAVB-AD, CAVB-JUV-L, and a “between-generation” (i.e., each pair is composed by one juvenile and one adult), referring to allele frequencies within each group. In each case, significance of the slope of the linear regression (b) was tested using 10,000 permutations. The pattern of SGS was represented with a correlogram using eight distance classes ([0–25]; [25–75]; [75–150]; [150–231]; [231–300]; [300–376]; [376–451]; [451–625]). Approximate confidence intervals (95%) for average kinship values for each spatial distance were obtained by permuting individual locations 10,000 times. Differences in SGS among the different categories were assessed using the Sp statistics with Sp = −b/(1 − F̂₍₁₎), with F̂₍₁₎ being the mean F_ij among individuals belonging to the first distance interval (Vekemans and Hardy, 2004).

For the global dataset (CAVB-AD and CAVB-JUV-L), we computed the neighborhood size (N_b) as N_b = 1/Sp and we deduced the moment estimate of the mean parent–offspring distance, σ_g, from N_b = 4^∗π^∗D_e^∗σ_g² (see Rousset, 1997) with the effective density, D_e = 201 ind^∗m^–2 (all adult colonies take part in reproductive events; Garrabou et al., 2017) and D_e = 100 ind^∗m^–2 (half of the adult colonies take part in reproductive events).

Sibship and Parentage Analysis

We characterized the pattern of relationships combining three different methods and considering two sets of individuals: (i) CAVB-JUV-L and CAVB-AD and (ii) CAVB-JUV-M.

First, we used the F_ij Nason’s kinship coefficients. Theoretically, F_ij discriminates among different types of relationships: F_ij = 0.25 or F_ij = 0.125 is expected for full-sibs and half-sibs, respectively. However, depending on allele frequencies, F_ij is influenced by the polymorphism of the molecular markers. We conducted simulations based on the allele frequencies accounting for the six populations to estimate F_ij values linked to four levels of relationships: unrelated (UN), half-sibs (HS), full-sibs (FS), and parent-offspring (PO). RELATED R-package (Pew et al., 2015) was used to simulate 1,000 pairs of genotypes for each of the four categories. This allowed us to define threshold values for the different categories of relationships (see section “Results”). We then computed pairwise F_ijs within CAVB-JUV-L, within CAVB-JUV-M, and between CAVB-JUV-L and CAVB-AD, referring to the allele frequencies of the whole dataset, and we used the threshold values to characterize the pattern of relationships.

The maximum-likelihood method implemented in COLONY 2.0.6.1 (Jones and Wang, 2010) was used to reconstruct the parental relationships among individuals based on their multilocus genotypes. COLONY reconstructs sibship and assigns paternity jointly considering the likelihood over the entire pedigree. Analyses were conducted considering two sets of individuals (i.e., CAVB-JUV-L/CAVB-AD and CAVB-JUV-M). For each set, we conducted three runs with different seed numbers considering females and males as polygamous and using the full-likelihood method with very high likelihood precision and very long length of run. For the first set, individuals from CAVB-AD were considered as potential fathers or mothers, while the second set was conducted without potential parents.

The two methods showed contrasted results with a low number of concordant relationships, which may be explained by different factors such as the presence of null alleles. Accordingly, we contrasted the results obtained with the two methods in a third step, using the statistical approach implemented in ML-Relate (Kalinowski et al., 2006). The analysis based on F_ij involving a high number of pairs (>60,000), we thus focused on juveniles and, more particularly, on juveniles with F_ij > 0.2625 and on adult/juvenile pairs with F_ij > 0.4228 (see section “Results”). Regarding COLONY, we accounted for all the inferred HS, FS, and PO relationships. Overall, the first two methods identified 28 PO, 193 HS, and 16 FS relationships (see section “Results”). For each of the relationships identified with one method but not with the other, we considered the identified relationship as the putative relationship and the UN relationship as the alternative relationship in ML-Relate. The putative relationship was tested based on the sampling distribution of a test statistic obtained after simulations of 1,000 genotypes under the alternative hypothesis (see Kalinowski et al., 2006 for details).

Results

Hardy–Weinberg Equilibrium and Genetic Diversity Estimators

Significant linkage disequilibrium after FDR correction was observed between Mic13 and Mic25 when considering all populations. However, no significant disequilibrium was observed when considering each population separately. The estimated frequency of null alleles was between 0.09 for PPL and 0.13 for ALI (mean over populations ± SE = 0.10 ± 0.01). The probability of identity (Pi) was 2.6 ^∗ 10^–11, supporting the validity of the set of microsatellites to infer relationships among individuals. Observed heterozygosity values were between 0.50 for DHM and 0.59 for PLU (mean over populations ± SE = 0.54 ± 0.04). The gene diversity H_e ranged between 0.67 (DHM) and 0.74 (ALE) (mean over populations ± SE = 0.71 ± 0.02). The f estimator of F_IS varied between 0.19 (PPL) and 0.23 (CAVB-AD) (mean over populations ± SE = 0.23 ± 0.04). Significant departure from panmixia was observed in all populations (significant p-Value after FDR correction). The lowest and highest values of Ar₍₅₆₎ were observed for PLU (9.71) and DHM (11.36) (mean Ar₍₅₆₎ over populations ± SE = 10.63 ± 0.65) (see Table 1). None of the genetic diversity estimators was significantly different among CAVB-AD vs. CAVB-JUV (all p-Values > 0.05).

SGS Among Populations

LnP(D) increased slightly for the whole range of considered K-values (i.e., no plateau) with the strongest increases observed for K = 2 and K = 3. Evanno’s method (Evanno et al., 2005) identified two different genetic clusters, while Peuchmaille statistics supported K = 2 (median/mean with threshold 0.7 and 0.8) and K = 3 (remaining statistics) (Figure 2A; Supplementary Material). For K = 2, the first cluster encompassed all the individuals from CAVB and DHM (mean membership coefficient = 0.93). The remaining individuals were grouped in a second cluster with a relatively high mean membership coefficient (0.87). When considering K = 3, this second cluster was divided in two clusters segregating the individuals from PLU (mean membership coefficient = 0.91). Those results were consistent considering null alleles as recessive alleles or missing data (data not shown).

FIGURE 2

Figure 2. Result of the clustering analysis of Corallium rubrum from the Scandola Marine Reserve conducted with STRUCTURE considering two and three clusters (K). (A) Each individual is represented by a vertical line partitioned in K–colored segments, which represent the individual membership fraction in K clusters. Thin and thick black vertical lines delineate the different populations and the different clusters, respectively. Sample names are shown above the assignment plots, while the mean membership coefficients for each cluster are shown below the assignment plots. (B) Scatter plot of the discriminant analysis of principal components (DAPC) based on a a–score of 46. Each dot corresponds to one individual (n = 580) from the seven populations (here CAVB-AD and CAVB-JUV were considered separately), which are represented by different colors. Inertia ellipses center on the mean for each population and include 67% of the individuals. The two axes of the scatter plot explained 52.1% and 27.3% of the total genetic variation in the dataset.

In the DAPC analyses, individuals from CAVB and DHM were separated from the remaining individuals along the first axis, while the second axis separated the individuals from PLU from the individuals from ALI, ALE, and PPL. Those two axes represented 79.4% of the total variation in the data (Figure 2B).

The global F_ST was 0.069. The pairwise F_STs ranged from 0.022 for ALI vs. ALE to 0.124 for CAVB vs. PLU. The genetic distance (F_ST/(1 − F_ST)) and the geographic distance (Ln(d)) were not correlated rejecting the IBD among the six populations (p-Value = 0.22). The exact tests for genotypic differentiation were significant at the global level and for all pairwise comparisons except between ALI and ALE, separated by 20 m (see Supplementary Material). Regarding the temporal differentiation, the pairwise F_ST between CAVB-AD and CAVB-JUV was low (0.001), but the genotypic differentiation was significant (p-Value < 0.01).

The lowest population-specific F_ST was observed for ALE (0.050; 95% CI: 0.036–0.067) whereas the highest value was observed for PLU (0.097; 95% CI: 0.070–0.123). Based on 95% CI, significantly different population-specific F_STs were observed for PLU vs. CAVB and PLU vs. ALE (Table 1).

Assignment Tests

Considering the whole dataset, 69 of the 419 individuals (16.4%) were identified as FGM, with eight (PPL) to 15 (DHM) FGMs per population. Considering assignment probability > 0.01 and allowing for multiple assignments, we assigned seven FGMs (10.1%) to three populations (CAVB, ALE, and PPL) indicating that 89.9% of the FGMs came from unsampled populations. Three populations (DHM, ALI, and PLU) did not produce any FGM. From the assigned FGMs, 50% came from CAVB.

When focusing on CAVB, 17 juveniles were not assigned to CAVB-AD (i.e., assignment probability < 0.005), while 10 juveniles were assigned to one of five other populations, with a probability higher than the assignment probability to CAVB-AD. Accordingly, those 27 juveniles were considered as migrants. Considering only assignment probability higher than 0.01, nine of these 27 (33.3%) migrants were assigned to DHM. The 18 remaining migrants were considered as coming from unsampled population(s) (assignment probability < 0.01 in all populations).

Fine-Scale SGS

We demonstrated an overall SGS based on the marginally significant and negative correlation between Nason’s kinship coefficients (F_ij) and ln(d) observed when considering the “global dataset” (b = −1.7 ^∗ 10^–3; p-Value = 0.04). From the neighborhood size (Nb = 575.3 individuals), we estimated the mean parent–offspring distance σ_g = 15.09 cm and 21.39 cm for D_e = 201 individuals^∗m^–2 and D_e = 100 individuals^∗m^–2, respectively. The decrease in kinship with geographic distance remained significant considering the CAVB-JUV-L dataset (b = −2.4^∗10^–3; p-Value = 0.03) and the CAVB-AD/CAVB-JUV-L dataset (b = −2.2^∗10^–3; p-Value = 0.04). No significant correlation was observed (b = −8.2^∗10^–4; p-Value = 0.29) when considering CAVB-AD. In spite of a trend toward stronger SGS for CAVB-JUV-L and CAVB-AD/CAVB-JUV-L (Sp = 0.0024 and 0.0022, respectively) compared to the whole dataset (Sp = 0.0017), no significant difference in Sp was observed (Figure 3 and Table 2).

FIGURE 3

Figure 3. Decomposing the spatial genetic structure among red coral colonies within the CAVB population. (A) Mean kinship coefficient (Fij) considering all the individuals (adults and juveniles) for each of distance class (d). Vertical lines correspond to standard errors. Dashed lines correspond to 95% confidence interval computed after 10,000 permutations of individual locations. The correlation among the genetic (F_ij) and geographic (ln(d)) distances among individuals was significant (p-Value = 0.04). (B) The global spatial genetic structure considering all the individuals (black line) is decomposed considering (i) only the pairs of adults (dotted line); (ii) only the pairs of juveniles (dark gray line); and (iii) only the pairs involving one adult and one juvenile (light gray line). The correlation among genetic (F_ij) and geographic (ln(d)) distances among individuals was significant when considering the juvenile pairs (p-Value = 0.04) and the adult–juvenile pairs (p-Value = 0.03) but not the adult pairs.

TABLE 2

Table 2. Decomposing the spatial genetic structure among individuals within Cave-b population: Summary statistics of the regression analyses between the geographic (ln(d)) and genetic (F_ij) distances among pairs of individuals.

Sibship and Parentage Analysis

The pairwise F_ij values for simulated unrelated (UN), halfsibs (HS), fullsibs (FS), and parent–offspring (PO) ranged from −0.1984 to 0.2625 (mean ± SD = −0.001 ± 0.071), −0.1160 to 0.4228 (0.126 ± 0.087), −0.0246 to 0.6393 (0.248 ± 0.106), and 0.0622 to 0.6556 (0.252 ± 0.077), respectively. Considering the overlaps of simulated F_ij distributions for the four categories of relationships (see Supplementary Material 4), we distinguished UN (F_ij < −0.1160) from related (HS, FS, and PO) pairs (F_ij > 0.2625). Within the related pairs, juvenile pairs with 0.2625 < F_ij < 0.4228 were considered as HS or FS. Finally, pairs with F_ij > 0.4228 were assigned to PO or FS depending on whether they involved one adult/one juvenile or two juveniles, respectively. Note that some pairs (e.g., all the pairs with −0.1160 < F_ij < 0.2625) were not further considered since we were unable to assign a relationship category.

For the 351 pairs involving the 27 migrant juveniles (CAVB-JUV-M), the mean F_ij (±SD) was 0.0084 (0.1124). We identified eight related pairs (F_ij > 0.2625; 2.2%), from which seven were considered as HS or FS (0.2625 < F_ij < 0.4228; 1.9%) and one as FS (F_ij > 0.4228; 0.3%). Forty pairs were considered as UN (F_ij < −0.1159; 11.4%). UN and HS pairs occurred among juveniles belonging to the same or to different quadrats, while the FS involved individuals sampled within the same quadrat.

When considering 139 local juveniles (CAVB-JUV-L), the mean F_ij (±SD) was 0.0166 (0.0981). Among the 9,591 pairs, 689 (7.18%) were considered as UN (F_ij < −0.1159), while 153 HS or FS (0.2625 < F_ij < 0.4228; 1.6%) and six FS (F_ij > 0.4228; 0.06%) pairs were revealed. Two FS pairs involved juveniles sampled at the two extremes of the transect while 66 UN pairs involved juveniles sampled in the same quadrat.

The mean F_ij (±SD) among adults (CAVB-AD) was 0.0086 (0.1007). Among the 20,100 pairs, we identified 1761 (8.76%) UN (F_ij < −0.1159) and 329 (1.63%) related (F_ij > 0.2625) pairs. From the related pairs and considering the CAVB-AD samples involving different generations, 306 (1.52%) were considered as HS or similar second-degree relatives (e.g., grand-parent/grand-offspring, avunculate) (0.2625 < F_ij < 0.4228), while 23 (0.01%) were considered as FS or first-degree relatives (F_ij > 0.4228).

The mean F_ij (±SD) for CAVB-AD vs. CAVB-JUV-L was 0.0104 (0.0990). Among the 30,351 pairs, we identified 2,629 (8.66%) UN (F_ij < −0.1159) and 412 (1.35%) related (F_ij > 0.2625) pairs. From the related pairs, 387 (1.52%) were considered as HS or similar degree relationships (0.2625 < F_ij < 0.4228), while 25 (0.08%) were considered as PO (F_ij > 0.4228).

Regarding COLONY analyses, no HS and two (0.6%) FS pairs (p > 0.6) were observed in migrant juveniles (CAVB-JUV-M). Considering local juveniles (CAVB-JUV-L) and adults (CAVB-AD), 10 (0.1%) FS and 40 (0.4%) HS pairs were identified (p > 0.6) in addition to three PO pairs. All the PO pairs involved the same parent.

Filtering these results using the statistical test implemented in ML-Relate, we retained 114 HS, 7 FS among the local juveniles (CAVB-JUV-L), and 11 PO pairs when considering the adults. Within the migrant juveniles (CAVB-JUV-M), two HS and one FS pairs were retained.

Discussion

Genetic Diversity and Structure in the Scandola Marine Reserve: Apparent Stability of the Diversity Pattern in Spite of a Dramatic Demographic Decline

We used complementary approaches to characterize the SGS among six red coral populations in the Scandola Marine Reserve. When considering the clustering approach, the first level of structure (K = 2) supports the occurrence of a geographic break among the northern (PLU, PPL, ALI, and ALE) and southern (CAVB and DHM) populations with high membership coefficients (0.87 and 0.93, respectively). For K = 3, individuals from PLU were segregated from the individuals from PPL, ALI, and ALE forming two subclusters, also characterized by high membership coefficients (0.91 and 0.81, respectively). This pattern was fully supported by the DAPC analyses, which discriminated among northern and southern populations along the first axis (52.1%), while individuals from PLU were separated from PPL, ALI, and ALE along the second axis (27.3%). When shifting at the population level, we found significant genotypic differentiation among all population pairs, with the exception of ALE vs. ALI, which belong to the same location (only apart by 15 m). Worthy of note is that the genetic distances between populations were not correlated with the geographic distances, rejecting the isolation-by-distance pattern (IBD). Computations of the population-specific F_STs identified PLU as genetically isolated, compared to the remaining populations, which showed relatively homogenous values. Accordingly, a geographic imprint explains part of the global pattern (i.e., northern vs. southern clusters). Yet, the differentiation of PLU points out the importance of local genetic drift effects in the survey area, which may potentially explain the lack of IBD usually reported in this species (e.g., Ledoux et al., 2010a; Aurelle et al., 2011). Non-exclusive hypotheses can explain this genetic differentiation of PLU. Located in a cave, PLU can be isolated from neighboring populations due to hydrodynamic factors (e.g., Costantini et al., 2018). PLU is also impacted by a dramatic demographic decline with more than 90% of biomass lost during the last 15 years due to warming-induced MMEs (Gómez-Gras et al., 2021). This decline can potentially increase the effect of genetic drift in this particular population.

Overall, the pattern reported here, combining different genetic clusters with a lack of IBD, matches the pattern previously reported in the same region almost 10 years ago (Ledoux et al., 2010a). The same happens for the genetic diversity estimators, with values estimated here similar to previously reported values (Ledoux et al., 2010a). Changes in heterozygosity values are difficult to detect (Riquet et al., 2017), and a delay in the detection of genetic erosion is expected in long-lived species with overlapping generations such as red coral (e.g., Hailer et al., 2006). However, this apparent stability in the pattern of genetic structure and diversity remains surprising considering the dramatic demographic decline reported in red coral in this area (Gómez-Gras et al., 2021). While deserving cautious interpretation, this result highlights the need to consider complementary metrics such as effective population sizes to characterize temporal trend in genetic diversity.

Integrating “Among Populations” and “Within Population” Demo-Genetic Approaches: Demographic Connectivity of Red Coral in the Scandola Marine Reserve

Our analyses suggested substantial demographic connectivity among the studied populations. Indeed, 16.4% of the mature individuals were identified as FGMs, in the six populations, including the drifting population of PLU. This result may appear counterintuitive in the context of the significant genetic structure observed even between populations separated by tens of meters (e.g., PPL vs. ALE). However, it is supported by recent studies in other Mediterranean octocorals reporting some levels of contemporary connectivity at least during recolonization steps (e.g., Arizmendi-Mejía et al., 2015; Padrón et al., 2018; Aurelle et al., 2020). Conciliating these, a priori, contradictory results imply further studies such as a detailed characterization of the balance between genetic drift and gene flow, during the course of population succession (i.e., from colonization to senescence). We hypothesize that the genetic imprint of migrants on a population will decrease from colonization to population maturation due to an increase of the genetic imprint of locally produced recruits. Our current understanding of the colonization processes (see below) combined to a potential underestimation of the influence of genetic drift on the patterns of spatial structure in Corallium rubrum, and other octocorals (Ledoux et al., 2015; Crisci et al., 2017; Pratlong et al., 2018) are in line with this hypothesis. Low reproductive success of migrants (i.e., limited genetic connectivity), due for instance to microenvironmental conditions favoring locally produced genotypes (e.g., Gorospe and Karl, 2013), may be another hypothesis to explore to solve this apparent paradox.

Following the filtering assignment, most of the FGMs identified (89.9%) came from unsampled populations. When considering the assigned migrants, two populations (CAVB and ALE) were identified as the main sources in the population network. In spite of a lack of isolation by distance, most of these assigned migrants came from neighboring populations separated by tens of meters (e.g., from CAVB to DHM), suggesting a low level of demographic connectivity among the Northern and Southern main genetic clusters. These analyses also support that the pool of migrants identified in each population came from (at least two) different sources. Accordingly, the multiple origins of migrants demonstrated in recently colonized and mature populations of other Mediterranean octocorals (Paramuricea clavata and Eunicella cavolini; see Mokhtar-Jamaï et al., 2011; Arizmendi-Mejía et al., 2015; Ledoux et al., 2018; Padrón et al., 2018; Aurelle et al., 2020) can be expanded to C. rubrum.

Focusing on the CAVB population, we refined the spatial features of the red coral demographic connectivity. Over the 166 juveniles analyzed here, 27 (16.2%) were considered as migrants, suggesting that, in spite of a substantial demographic connectivity, the vast majority of the juveniles was locally produced. A high level of self-recruitment was previously suggested as a central process in red coral population dynamics (Costantini et al., 2007b; Ledoux et al., 2010b). Considering the migrant juveniles sampled in CAVB (CAVB-JUV-M), the sibship analyses revealed the occurrence sibship aggregation. We identified two pairs of halfsibs and one pair of fullsibs in these juveniles, the two individuals of the three pairs being sampled in the same quadrat. Accordingly, collective dispersal, the dispersal of at least two sibs originating from the same reproductive event in a distant population (Broquet et al., 2013), may occur in red coral. The factors driving these events are still largely unknown (but see Ottmann et al., 2016). Beside abiotic factors (e.g., no mixing of mass water), biotic factors both at the larval (behavior) and adult (synchronous larval release) levels, combined to species reproductive strategies, have been called on to explain collective dispersal. In particular, most of the collective dispersals involving marine invertebrate larvae were reported in brooding or partially brooding species (e.g., Riquet et al., 2017, but see Dubé et al., 2020). The internal fertilization (Santangelo et al., 2003), brooding period (Bramanti et al., 2005), and larvae active swimming behavior (Martínez-Quintana et al., 2015) may drive collective dispersal in red coral.

Finally, we gained a first insight into the temporal features of red coral demographic connectivity. The assignment analyses suggest a lack of constant pattern, at least for the CAVB population. In the juveniles sampled in CAVB (CAVB-JUV), most of the migrants came from DHM. This DHM to CAVB migration was not reported when considering the pool of mature colonies. Temporal variability in the genetic make-up of larvae and juveniles has been previously reported in various marine species from fishes (e.g., Salles et al., 2016) to invertebrates (e.g., David et al., 1997; Riquet et al., 2017). This result is also in line with the recruitment variability reported in red coral (Garrabou and Harmelin, 2002). Considering the overlapping generations in this species, identifying the origin of migrants in different cohorts is a challenging task that was recently successfully achieved by analyzing different recruitment events, using settlement plates (Costantini et al., 2018). In this experimental framework, genetic heterogeneity among juveniles coming from two consecutive reproductive events was revealed (Costantini et al., 2018), consistently with the temporal variations reported here.

While the link between demographic and genetic connectivity remained to be fully characterized (but see below), the characterization of red coral demographic connectivity opens new perspectives for the management of the Scandola Marine Reserve.

Spatiotemporal Genetic Structure Among Individuals in a Near-Pristine Population: Complementing the Eco-Evolutionary Feedbacks Underlying Red Coral Population Dynamics

We demonstrated a statistically significant SGS when considering the adults and the local juvenile-stage classes in CAVB, a population with a near-pristine demographic structure (high size colonies at high density = 201 colonies^∗m^–2; Garrabou et al., 2017). A statistically significant SGS within population has been previously demonstrated in red coral (Ledoux et al., 2010b), albeit in a declining population characterized by different demographic features (high density of small size colonies; 160 colonies^∗m^–2). This shared pattern in such contrasted demographic situations suggests that SGS among colonies may be a common feature of red coral populations. Interestingly, the moment estimator of the neighborhood size (Nb) estimated in CAVB was one order of magnitude higher than the Nb previously reported in the declining population (575.3 vs. 75 individuals), resulting in a slightly lower effective dispersal (σ_g; from 15.1 to 21.4 cm vs. from 22.6 to 32.1 cm). While this comparison should be taken with caution owing to the different sampling strategies (quadrats vs. individual sampling), a decrease of σ_g with a higher effective density is expected (see Vekemans and Hardy, 2004). This feedback between effective density and dispersal (i.e., higher dispersal in low-density population) was proposed as a process buffering, to some point, the increase of genetic drift expected in declining populations and may be critical for the maintenance of red coral populations (Ledoux et al., 2020).

Such restricted σ_g seems conflicting with the pattern of demographic connectivity previously drawn based on assignment method. One should keep in mind the different dispersal components (demographic vs. genetic) estimated by the two methods (Broquet and Petit, 2009). Moreover, estimation of σ_g based on SGS among individuals may underestimate long-distance dispersal (Leblois et al., 2004). Consistent with a potential underestimation of effective dispersal, the sibship analyses revealed for instance one parent–offspring pair separated by 300 cm, suggesting that longer-distance dispersal may happen.

We then decomposed the SGS accounting for the stage class of the individuals (i.e., adults and juveniles). We first revealed a significant SGS among adults and juveniles and showed that the significant SGS was retained in the juveniles but not in the adults. The significant SGS among stage classes corroborates the central role of restricted dispersal in the structure and in the dynamics of red coral populations. In spite of the demographic connectivity among populations previously discussed and of the occurrence of long-distance dispersal events, larvae recruit preferentially in the close vicinity of their parents (Ledoux et al., 2010b, this article). This restricted dispersal builds up the SGS within local juveniles. While potential methodological limitations (i.e., overlapping generations in the adult sample; see Berens et al., 2014) may take part in the random pattern observed among adult colonies, the decrease in SGS from juveniles to adults can be linked to an imprint of mortality during the early life stage (e.g., Chung et al., 2003; Hampe et al., 2010), corroborating the high mortality rates in the juveniles compared to adults reported by demographic surveys (Garrabou and Harmelin, 2002; Bramanti et al., 2007; Santangelo et al., 2012). The process underlying this demographic transition (e.g., self-thinning vs. density-independent mortality) remains a matter of discussion in red coral (Cau et al., 2016; Garrabou et al., 2017).

Shifting From Patterns to Processes to Improve the Conservation of Red Coral in the Scandola Marine Reserve and the Restoration of Declining Populations

Implemented more than 40 years ago, the Scandola Marine Reserve belongs to the 0.23% of the Mediterranean that is fully protected (Claudet et al., 2020). Beside harboring a unique and near-pristine shallow population of red coral (Garrabou et al., 2017), the benefits of the management conducted in this MPA on red coral are noteworthy (Linares et al., 2010, 2012). We revealed here contrasted connectivity among populations forming the red coral network. We highlighted the restricted connectivity, at least on a contemporary timescale, among the Northern and Southern parts of the MPA and, accordingly, the need to consider these localities as distinct management units. The potential for genetic isolation at local scale, as observed for PLU, is also of management interest. Indeed, isolated populations may be targeted by restoration to increase their effective population size and to buffer potential drift and related negative effects (e.g., expression of deleterious alleles and decrease population fitness; see Frankham, 2005; Garner et al., 2020). Still on the management prioritization, we identified two particular populations, namely, CAVB and ALE, as the main sources of migrants within the network. Additional populations must be sampled to refine this connectivity pattern. Moreover, while our results suggested relatively homogeneous levels of genetic diversity among the different populations, this pattern does not account for adaptive genetic variations (i.e., variations that produce advantages in fitness), which may be limited in this particular area (Ledoux et al., 2015). This is a crucial knowledge gap for an effective management of the MPA, which should be filled using available genomic resources (Pratlong et al., 2015, 2018).

Restoration actions are a promising tool in red coral at local scale (Montero-Serra et al., 2017). Restoration protocols usually rely on different components including the number of individuals to transplant, the transplanting design, and the choice of source populations (Weeks et al., 2011; Mijangos et al., 2015; Breed et al., 2019). Ledoux et al. (2010b) suggest to create small size and dense patches (tens to a hundred of colonies over half square meter) of reproductive colonies to mimic the patterns of structure characterizing red coral populations. Our results support this preliminary statement. We confirmed here the importance of restricted dispersal and of interactions among genetically and spatially related individuals in the dynamics of red coral populations. The size of a functional restored patch, i.e., promoting genetic interactions among colonies, depends on the density of transplanted colonies. Interestingly, in the two contrasted density situations analyzed so far (i.e., a declining population vs. near-pristine population), the mean effective dispersal remains within the order of tens of cm, with occasional more distant dispersal events. This restricted range of dispersal distances sets the scale for restoration actions (i.e., tens of cm). Yet, the parentage analyses revealed the potential for long-distance effective dispersal within population, beyond the mean parent–offspring distance (σ_g). Accordingly, it may be more efficient to restore a particular location creating small and dense patches of colonies (tens to a hundred of colonies in 0.5 to 1 m² area), separated by 2 to 3 m, in order to “capture” these long-distance dispersers, rather than a low-density and continuous patch. Finally, sourcing is a critical aspect of restoration. The identification of FGMs from different origins supports the need to consider different populations for sourcing. However, one important gap here relies on the ability of red coral colonies to deal with the ongoing warming (Garrabou et al., 2019). In this context, the identification of the genomic factors and eco-evolutionary processes driving the response to thermal stress in red coral (Pratlong et al., 2015, 2018) is the next critical step to improve restoration actions in the warming context.

Overall, the spatially integrated demo-genetic approach presented here revealed poorly known aspects of red coral eco-evolution, with strong implications for the species conservation and restoration. While targeting this species is particularly relevant considering its key ecological role in the coralligenous community, we call for the development of comparative studies among species to increase the robustness of the management inputs coming from eco-evolution studies.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author Contributions

J-BL, J-MD, and JG designed the study. J-BL, IM-S, J-MD, CL, and JG collected the samples. CG, PL-S, and J-BL performed the molecular analyses. CG and J-BL conducted the statistical analyses and led the writing of the manuscript. All authors contributed critically to the drafts and gave the final approval for publication.

Funding

This research was supported by the Strategic Funding UIDB/04423/2020 and UIDP/04423/2020 through national funds provided by the FCT – Foundation for Science and Technology and European Regional Development Fund (ERDF), in the framework of the program PT2020, the Spanish MINECO (CGL2012-32194), the TOTAL Foundation PERFECT project, the MIMOSA project funded by the foundation Prince Albert II de Monaco, and the European Union’s Horizon 2020 Research and Innovation Program under grant agreement N° 689518 (MERCES). This output reflects only the authors’ view, and the European Union cannot be held responsible for any use that may be made of the information contained therein. J-BL was funded by an assistant researcher contract framework of the RD Unit—UID/Multi/04423/2019 – Interdisciplinary Centre of Marine and Environmental Research—financed by the European Regional Development Fund (ERDF) through COMPETE2020 – Operational Program for Competitiveness and Internationalization (POCI) and national funds through FCT/MCTES (PIDDAC). Genotyping was performed at the Genome Transcriptome Facility of Bordeaux (grants from the Conseil Régional d’Aquitaine n 20030304002FA and 20040305003FA, from the European Union FEDER n 2003227 and from Investissements d’Avenir ANR-10-EQPX-16-01). This work acknowledges the “Severo Ochoa Centre of Excellence” accreditation (CEX2019-000928-S). We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2021.633057/full#supplementary-material

References

Abbiati, M., Santangelo, G., and Novelli, S. (1993). Genetic variation within and between two Tyrrhenian populations of the Mediterranean alcyonarian Corallium rubrum. Mar. Ecol. Prog. Ser. 95, 245–250. doi: 10.3354/meps095245