Genetic Diversity and Local Connectivity in the Mediterranean Red Gorgonian Coral after Mass Mortality Events

Abstract

Estimating the patterns of connectivity in marine taxa with planktonic dispersive stages is a challenging but crucial task because of its conservation implications. The red gorgonian Paramuricea clavata is a habitat forming species, characterized by short larval dispersal and high reproductive output, but low recruitment. In the recent past, the species was impacted by mass mortality events caused by increased water temperatures in summer. In the present study, we used 9 microsatellites to investigate the genetic structure and connectivity in the highly threatened populations from the Ligurian Sea (NW Mediterranean). No evidence for a recent bottleneck neither decreased genetic diversity in sites impacted by mass mortality events were found. Significant IBD pattern and high global F_ST confirmed low larval dispersal capability in the red gorgonian. The maximum dispersal distance was estimated at 20–60 km. Larval exchange between sites separated by hundreds of meters and between different depths was detected at each site, supporting the hypothesis that deeper subpopulations unaffected by surface warming peaks may provide larvae for shallower ones, enabling recovery after climatically induced mortality events.

About the Authors:

Joanna Pilczynska

* E-mail: [email protected]

Affiliations Departamento de Biologia and CESAM - Centro de Estudos do Ambiente e do Mar, Universidade de Aveiro, Aveiro, Portugal, Department of Earth and Environmental Sciences, University of Pavia, Pavia, Italy

Silvia Cocito

Affiliation: ENEA, Marine Environment Research Centre, La Spezia, Italy

Joana Boavida

Affiliation: CCMAR - Centro de Ciências do Mar, Universidade do Algarve, Faro, Portugal

Ester Serrão

Affiliation: CCMAR - Centro de Ciências do Mar, Universidade do Algarve, Faro, Portugal

Henrique Queiroga

Affiliation: Departamento de Biologia and CESAM - Centro de Estudos do Ambiente e do Mar, Universidade de Aveiro, Aveiro, Portugal

Introduction

Extreme weather events, including floods, heat waves and droughts, are currently emerging as one of the most important facets of climate change, and a growing body of literature is focused on extreme events [1]. Anomalous and extreme events due to global warming have increased considerably during recent decades in temperate regions such as the Mediterranean Sea and an increase in the frequency of heat wave extremes of 200–500% is predicted at the end of the twenty-first century [2]. Extreme events, together with other sources of mortality caused by human impact, such as overfishing and environmental pollution, may cause significant impacts on genetic diversity as a result of population size decrease [3], as has been observed after mass mortality events [4; 5].

If natural populations consist of reduced numbers of individuals, loss of genetic variability may dramatically influence the populations themselves, since genetically depauperate populations might fail to adapt to future environmental changes, eventually causing their disappearance. Due to the predicted intensification of weather extremes, many species may not have sufficient genetic potential for the evolution of strategies able to mitigate their impact. Nevertheless, genetic diversity reduction after mass mortality events is not a rule in the marine realm [3; 6]. There is an urgent need to carry out dedicated research to acquire more extensive, sound data across a range of life-history and demographic features, with the ultimate aim of formulating better predictions about the role of catastrophic disturbances in determining genetic structure and genetic diversity.

Connectivity patterns and gene exchange among populations are major research topics in marine ecology and are essential for the planning of marine reserves [7], to function as interconnected networks that can supply recruits to sites that undergo population bottlenecks or local extinctions. Genetic recovery from events of mass mortality, either natural or human-induced, may be dependent on the possibility of dispersal from external sources or from small local populations at refugial pockets [8–10]. This can result in much slower recovery of genetic diversity than population density, because demographic recovery is possible from a few founder or bottleneck survivors but genetic recovery requires extensive levels of connectivity or a long period of time [4]. However, inferring levels of connectivity is particularly challenging in the marine environment, where many species disperse exclusively by means of planktonic propagules, such as larval or spore stages, in a 3-dimensional fluid medium. Direct estimation of dispersal by tracking marine propagules is often not feasible in most instances, because of their small size and the unbounded nature of the marine environment [11; 12], so indirect methods, such as the use of neutral genetic markers [13; 14], are commonly applied.

The red gorgonian Paramuricea clavata (Risso, 1826) is a key species of sublittoral rocky habitats [15], widespread in the western Mediterranean Sea and in the Adriatic Sea [16] and less common in the Aegean Sea [17]. The species is a surface brooder with a short larval dispersal phase [18; 19]. In situ observations indicate that the dispersive stage of the larvae may last only a few minutes and thus the larvae are expected to settle near the mother colony [18]. This mechanism does not favor dispersal, but possibly decreases larval mortality and wastage, contributing to the replenishment of local populations [20]. However, metamorphosis may be delayed in the laboratory for up to 25 days after egg collection, suggesting high dispersal capacity, at least under certain conditions [18]. In 1999 and 2003, two episodes of mass mortality, connected with increased water temperature reaching 24°C [21], affected several populations of benthic suspension feeders in the northwestern Mediterranean, causing a drastic decrease of P. clavata colony density [22–25]. The number of colonies affected by partial mortality decreased with depth [22; 23; 25], because mortality affected the community from the surface to approximately 20 m, the approximate depth of the thermocline. Hereafter we designate as deep the sites with colonies living below the thermocline, which were not exposed to abnormally long lasting summer conditions [26]. Therefore, the clear differences in the number of damaged colonies and the extent of injury were visible between shallow populations (dwelling above the thermocline) and the deep ones, even when the depth difference was only of a few meters [22; 25]. As hypothesized by Cerrano & Bavestrello [27; 28], the deep dwelling subpopulations may act as a reservoir, supplying larvae to shallower sites, especially because large colonies with high fecundity rates survived there. However, little is known about short distance dispersal of the species. Surveys based on microsatellite loci in the Mediterranean Sea, showed that P. clavata exhibits a high level of genetic differentiation at the small and large spatial scales, which is consistent with the short larval dispersal displayed by the species [29–31]. Connectivity among populations of P. clavata separated by less than 14 km was weak in the study of [31], since the mean immigration rate was below 9% and most of the immigrants came from neighboring populations located only hundreds of meters away.

In the present study we used microsatellite markers to study the genetic structure of P. clavata populations in an area impacted by mass mortality events, in the Ligurian Sea, Mediterranean. Our objectives were twofold. The first objective was to investigate a possible bottleneck effect of past mortality events on genetic diversity in a region where some populations were highly affected by mass mortality events and others were not. In this case we are expecting to detect the decrease of genetic diversity at sites impacted by mortality. Our second objective was to understand the connectivity patterns and migration at the scale of a few tens of km since up to date very few studies have examined short distance migration in the species. We expect to detect migrations between non-impacted reefs and impacted ones, supporting recovery of damaged populations. Our research contributes to better understand the mechanisms that enable recovery of threatened populations by providing data about larval migrations between impacted and healthy populations of the red gorgonian. The results should prove to be particularly valuable for the conservation of soft corals communities and thus the overall marine biodiversity.

Material and Methods

Sampling

P. clavata colonies were sampled by scuba divers following a hierarchical sampling design. Samples were taken from three sites, Punta Mesco, La Spezia and Livorno, in the Ligurian Sea (NW Mediterranean). At each site two different reefs were chosen (Fig 1, Table 1). In the remaining of the text, each reef will be referred to by its code from Table 1. At Punta Mesco and Livorno reefs were sampled at different depths. At these two sites shallow reefs (<25 m depth) were impacted by mass mortality events, whereas colonies dwelling below 25 m (referred as deep reefs in the remaining of the text) remained non-impacted ([32]; Di Fiore M., pers com). At La Spezia both reefs were sampled at the same depth, since the rocky cliff ends on a muddy bottom at 24–25 m, and therefore both reefs were damaged during the mass mortality (Table 1). Populations impacted by mass mortality recovered significantly. Colonies density in La Spezia population declined from over 35 to nearly 8 colonies m^-2 shortly after the mass mortality and increased to 20 colonies m^-2 four years later [33]. In Portofino (Ligurian sea), located near Punta Mesco, colonies density decreased after the mortality event from nearly 20 to 5 colonies m^-2, but 3 years after the event the density recovered to pre-mortality levels [23]. We do not have data from Punta Mesco and Livorno but we may expect that the recovery pattern was similar and, therefore, all impacted populations consist of survivors and new colonies that settled after the mortality. Sites were separated by distances ranging from 20 to 60 km (Fig 1) whereas reefs within each site were separated by 200 to 300 m. Site 1, Punta Mesco, is located in the Cinque Terre Marine Protected Area. A small portion of colony branch (around 3–4 cm of branch tip) was taken from up to 30 different colonies randomly chosen from each reef. The branch tip from each colony was stored individually in a plastic tube underwater. Samples were placed on ice during transport and preserved in ethanol after arrival to the laboratory, no later than 3 hours after collection.

[Figure omitted. See PDF.]

Fig 1. Sampling sites in the Mediterranean.

Reefs impacted by mass mortality events (in black) and healthy reefs (in white).

https://doi.org/10.1371/journal.pone.0150590.g001

[Figure omitted. See PDF.]

Table 1. Sampling sites characteristics.

https://doi.org/10.1371/journal.pone.0150590.t001

The Cinque Terre Marine Protected Area authorized ENEA to conduct fieldwork and sampling at Punta Mesco (reef 1 and 2). La Spezia Islands (reef 3 and 4) are included in a recently established marine conservation area where ENEA is authorized by the Regional Natural Park of Porto Venere to conduct fieldwork and sampling. No specific permission was required for sampling at Livorno. P. clavata is not endangered nor protected.

Microsatellite analysis

Coral DNA was extracted using E.Z.N.A. Mollusc DNA Kit according to the manufacturer handbook. We analyzed 10 microsatellite, developed by Agell et al. [35] and Mokhtar-Jamaï et al. [36] following the protocols published by the authors. Loci Par_a, Par_b, Par_d, Par_f and Par_m were amplified from total genomic DNA in 10 μl solution of dNTPs (0.25 mM each), selected primers (0.25 μM each), 4 mM of MgCl₂, 1x manufacturer-supplied buffer and 0.25 u DFS-Taq DNA Polymerase (Bioron). The PCR program was: 2 min 94°C, (10 sec 94°C, 20 sec AT, 1 min 72°C)x30, 5 min 72°C. Annealing temperatures (AT): Par_a: 59°C, Par_b: 47°C, Par_d: 51°C, Par_f, Par_m: 52°C. To amplify loci Parcla_9, Parcla_10, Parcla_12, Parcla_14 and Parcla_17, total genomic DNA was dissolved in 10 μl solution of dNTPs (125 μM each), selected primers (0.5 μM each), 0.25 u GoTaq^® DNA Polymerase (Promega) and 1x manufacturer-supplied PCR buffer. The PCR program was: 3 min 94°C, (1 min 94°C, 1 min 60°C, 1 min 72°C)x30, 5 min 72°C. The length of amplified fragments was analyzed on an ABI 3730XL Genetic Analyzer using an internal size standard (GeneScan 500 LIZ). The analysis of DNA fragment length was performed with STRand [37]. Scored microsatellite fragment sizes were then visualized in R environment using the MsatAllele_1.02 package to track and reanalyze scoring errors.

Genetic diversity

MICRO-CHECKER v.2.2.3 [38] was used to estimate null allele frequency and to check for scoring errors owing to stutters and large allele dropout. Linkage disequilibrium among all pairs of loci was tested in GENEPOP 4.2. [39; 40] with significance levels determined by the Markov chain method (dememorization = 5000, batches = 500, iterations = 10 000).

Observed (H_o) and Nei’s [41] unbiased expected heterozygosity (H_e) were computed in GENETIX v.4.05 [42]. The rarefaction procedure implemented in HP-RARE software [43] was used to estimate allelic richness (A_r) and private allelic richness (A_p). The minimum number of genes was set to 25 (the minimum sample size). Differences in heterozygosity and allelic richness between impacted and healthy reefs were tested using Kruskal-Wallis. The power of the test and minimum sample size to achieve 90% power were calculated in PASS 14 [44]. Single and multi-locus Weir & Cockerham’s [45] f estimator of F_IS were calculated using GENEPOP 4.2. Departures from Hardy–Weinberg (HW) equilibrium within sample for each locus and over all loci were tested in GENEPOP 4.2. The level of significance was determined by the Markov chain method using the default parameters (dememorization = 1000, batches = 100, iterations = 1000).

Populations that experienced a recent mass mortality are predicted to lose allelic diversity faster than heterozygosity [46] and thus exhibit a heterozygosity excess relative to the heterozygosity expected from the observed number of alleles. To establish whether there is a heterozygosity excess or deficit, the BOTTLENECK software [47] computes a distribution of the expected heterozygosity under the assumption of mutation-drift equilibrium, calculated from the observed number of alleles, and compares it to the heterozygosity expected under Hardy-Weinberg equilibrium. The presence of a possible bottleneck effect was tested using 9999 simulations under the Two-Phase Model (TPM), since it generally fits microsatellite evolution better than either pure stepwise or infinite allele models [46, 47]. TPM assumes that the majority of mutations are single steps, when alleles increase or decrease by one repeat unit. The mutation sizes for remaining mutations are drawn from a geometric distribution and larger mutations are rare, but do occur. The frequency of step mutations was set to 0.9 (ps), the variance of mutations to 12 and the Wilcoxon test was used to test the null hypothesis of no significant heterozygosity excess.

Population differentiation and connectivity

Global and pairwise Weir & Cockerham’s [45] estimator of F_ST was estimated in GENEPOP. The genotypic differentiation between all pairs of reefs was tested in GENEPOP with default parameters.

The isolation by distance pattern was tested in GENEPOP. The shortest possible distance over sea between each reef pair was measured using Google Earth 7.1.2.2041. The relationship between genetic distance [F_ST/(1-F_ST)] and spatial distance [km] was tested using a Mantel test (n = 2000).

In order to quantify genetic variation within reefs, among reefs within a site and among sites a hierarchical analysis of molecular variance (AMOVA) was performed in ARLEQUIN 3.5 [48]. The significance of these variance components was tested using 50000 permutations.

The Bayesian approach implemented in STRUCTURE v.2.2 was used to investigate population structure. The recessive allele option was used to deal with null alleles [49]. The number of clusters (K) in the data set was evaluated under the admixture model with correlated allele frequencies. First run of 10 iterations, burnin of 10000 and MCMC = 50000 was computed for K from 1 to 10. The value of K that captures the major structure in the data was selected based on the plot of logarithm of the likelihood of observing the data [LnP(D)] as a function of K [50]. STRUCTURE was then run 30 times for K values from 2 to 6. The results were merged in CLUMPP [51] and graphically displayed in DISTRUCT [52]. The analysis was then repeated in the groups defined in the first run of STRUCTURE, when K = 3, to search for substructure within the groups.

A Bayesian assignment method [53] implemented in GENECLASS2 [54] was used to detect putative first generation migrants (F0). A Monte Carlo resampling method, as described in Paetkau et al. [55], was performed to evaluate each individual’s probability of belonging to a population from each reef.

Whenever multiple tests were conducted (linkage disequilibrium, HW equilibrium, genetic differentiation), the level of significance was adjusted using a false discovery rate (FDR) [56].

Results

Genetic diversity

No large allele dropout was detected by MICRO-CHECKER at any locus, but evidence of scoring errors due to stuttering was found in Parcla_12 and this locus was excluded from further analysis. The mean null allele frequency across all reefs varied from 0 for Parcla_10, Parcla_17 and Par_f to 0.18 for Par_m. No significant linkage disequilibrium was observed between any pair of loci (all p>0.05 after FDR correction), thus all loci were considered as genetically independent.

Eight of the loci were polymorphic at all sites, whereas one locus, Par_b, was monomorphic at LaSp3 and LaSp4 according to the 0.95 frequency criterium. The total number of alleles ranged from 4 for Par_b, Par_d and Par_f to 27 for Par_m. Unbiased expected heterozygosity (H_e) varied between 0.57 at PMes1 to 0.66 at Liv6, with a mean value of 0.62 (0.03). H_o ranged from 0.42 for LaSp4 to 0.59 at Liv5, with a mean value of 0.51 (0.07) (Table 2). The lowest allelic richness (A_r) was found at PMes1 (4.77) whereas the highest value was found at PMes2 (5.84). Private allelic richness (A_p) varied from 0.34 at PMes1 to 1.17 at PMes2 (Table 2). There was no evidence that reefs affected by mass mortality events had lower genetic diversity than healthy ones (Kruskal-Wallis, all p>0.05). The power of the test was low for He (20%) and Ar (20%), and moderate for Ho (70%) and Ap (66%). In our case, to achieve test power reaching 90% for the 0.05 confidence level, sample size would need to be at least 166 colonies per reef. Multilocus F_IS ranged from 0.08 at Liv5 to 0.31 at LaSp4 (Table 2). When all loci were examined separately, the lowest value was -0.26, for Par_b, and the highest one, 0.70, for Par_a.

[Figure omitted. See PDF.]

Table 2. Measures of genetic diversity.

https://doi.org/10.1371/journal.pone.0150590.t002

Significant heterozygote deficiency was detected at all reefs (Table 2). However, the departures from HW equilibrium were not evident for all loci at all sites. Populations from all reefs revealed a departure from HW equilibrium at locus Par_m, but the highest frequency of null alleles was found in this locus and evidence for null alleles was also detected at all reefs for this locus. When Par_m was excluded from the analyses, heterozygote deficiency was still significant at all reefs except Liv5 after FDR correction.

No evidence for a recent genetic bottleneck was detected for any of the investigated reefs. A heterozygote excess (an indicator of recent bottleneck) was not found at any reef (Wilcoxon, all p>0.58) neither was a heterozygote deficit (i.e., a sign of expansion) (all p>0.08).

Population differentiation and connectivity

Global F_ST value was 0.118, whereas pairwise comparisons between all pair of reefs varied from -0.003 between Liv5 and Liv6 to 0.189 between PMes1 and Liv6 (Table 3). All comparisons were significant (p<0.05), except Liv5 and Liv6, which revealed also no statistical differences in genotypic composition (Chi2 = 18.5, df = 18, p = 0.42).

[Figure omitted. See PDF.]

Table 3. Pairwise F_ST values.

https://doi.org/10.1371/journal.pone.0150590.t003

The correlation between F_ST/(1-F_ST) and distance was significant (p = 0.01), supporting an isolation by distance model of gene flow in P. clavata at the local scale (Fig 2).

[Figure omitted. See PDF.]

Fig 2. The isolation by distance pattern for P. clavata.

Linear regression of the genetic distance measured as F_ST / (1—F_ST) over the geographic distance (m).

https://doi.org/10.1371/journal.pone.0150590.g002

The AMOVA (Table 4) revealed that a highly significant percentage (11.19%) of the total genetic variation occurred among sites, whereas a smaller, but still significant, percentage of variation was explained by differences among reefs within sites (1.17%). Indeed, most of the variance was explained by differences within reefs (87.64%) and this was highly significant.

[Figure omitted. See PDF.]

Table 4. AMOVA.

https://doi.org/10.1371/journal.pone.0150590.t004

In the STRUCTURE analysis the plot of LnP(D) as a function of K revealed a plateau for K≥3 (S2 Fig). Samples were divided into 3 clearly separated clusters, each of them grouping the two reefs from the same site (Fig 3). Samples from Livorno (Liv5 and Liv6) displayed a high coefficient of population membership while reefs from Punta Mesco (PMes1 and PMes2) and La Spezia (LaSp3 and LaSp4) showed a low level of admixture. A second run of STRUCTURE did not reveal any genetic structure within the groups (data not shown). When K was set to 2, reefs from Livorno were separated from the others, whereas a K value over 3 did not reveal any additional structure.

[Figure omitted. See PDF.]

Fig 3. Clustering analysis.

Population structure revealed by the first run of clustering analyses in STRUCTURE. Each individual is represented by a vertical line, divided into segments representing the proportion of the genome of the individual that is assigned to each cluster. The number of clusters was set to 2–6. Reefs are separated by a black vertical line.

https://doi.org/10.1371/journal.pone.0150590.g003

According to GENECLASS2, 37.9% of the colonies were first generation migrants (F0). Most of them (31.6%) were exchanged between reefs separated by 200–300 m within the same site. The contribution of F0 migrants from the same site varied from 23.3 to 31.0% for PMes1 and LaSp3, respectively (Table 5). Between-sites migration was detected only between La Spezia and Punta Mesco, separated by 20km. The data suggest that migration occurred in both senses, with 10.0% of the colonies from PMes (percentage from both reefs pooled together) assumed as migrants from LaSp and 8.5% of the colonies form LaSp contributed from PMes. Migrants exchanged between Liv5 and Liv6 were not evaluated because these sites do not differ significantly (F_ST = -0.00332).

[Figure omitted. See PDF.]

Table 5. First generation migrants.

https://doi.org/10.1371/journal.pone.0150590.t005

Discussion

Genetic diversity

Our results did not support the hypothesis of genetic effects of past mass mortality events in P. clavata from the Ligurian Sea. No differences in genetic diversity between reefs affected by mass mortality events and healthy ones were detected. Additionally, the levels of genetic diversity found in P. clavata from the Ligurian Sea were not much lower than the range of values reported previously for the species [29] and other Mediterranean [57] and tropical corals [58–60]. Numerous gorgonian populations in the Mediterranean were affected by mass mortality in the recent past [22–25] and there is the concern that these may have gone through genetic bottlenecks as a result of decreased population densities. However, confounding effects may obscure bottleneck results. For example, in Corallium rubrum (Linnaeus, 1758), 16 of 40 Mediterranean shallow populations showed a sign of recent expansion after a bottleneck [61], although a Wahlund effect could not be discarded. In a study of P. clavata from Ibiza (Balearic Islands, Spain), partial mortality, measured as the proportion of colony tissue damaged, was negatively correlated with effective population size, mean number of alleles per population and proportion of recent migration rates [31]. These results indicated that populations with colonies that are partially affected by mortality are less diverse, undergo a larger effect of drift and receive less immigrants than healthy gorgonian populations. Additionally, in our study, allelic richness was slightly higher in populations non-impacted by mass mortality, which is in agreement with published results [31]. Yet, in spite of the extensive density reductions that mass mortality events had on the Ligurian red gorgonian (ca. four-fold; see methods section), our results show no difference in genetic diversity between healthy and impacted sites. This may suggest that the adaptive potential of surviving populations was not reduced. Nonetheless, this interpretation should be cautious given that our sampling design was limited to six reefs, among which only two were not impacted by mass mortality. Consequently the test presents low to moderate power (20–70%) to detect differences between genetic diversity of healthy and impacted sites (high probability of making a type II error, i.e., failing to detect existing differences). The need to increase power in such tests with higher number of samples and loci has been highlighted in previous studies [62]. To achieve a power of 90%, we would need to sample 166 colonies per reef (ca. 1000 colonies in total), which is unrealistic, because of the logistic effort and the limited number of colonies at some reefs. It is also noteworthy that heterozygosity excess tests have been shown to have a limited power to detect mild-bottlenecks of 10-1000-fold population declines [63]. Even if past mass mortality events have affected genetic diversity, it may not be detected by our study, because of the lack of power of the statistical tests and other confounding factors such as recent expansion or Wahlund effect.

The departures from Hardy-Weinberg equilibrium found in the present research, indicated by high and significant F_IS values, are typical for species exhibiting low larval dispersal, and were already reported in P. clavata populations [29]. Heterozygote deficit as a result of inbreeding was previously reported in coral species characterized by having a short larval dispersal [58; 60; 64]. Our findings confirm previous studies and are in accordance with the reproductive biology and larval ecology of P. clavata. Planulae exhibit negative phototaxis and negative buoyancy and settle near the mother colony [18], increasing the probability of subsequent mating with closely related individuals, associated with reduced dispersal of the male gametes also. High F_IS values may also result from a high number of F0 migrants, analogously to Wahlund effect, since migrants come from genetically distinct population. The presence of null alleles may also partially explain the result, but heterozygote deficit was still prevalent in the absence of the main locus suspected to have null alleles.

Genetic structure and connectivity

The high value of global F_ST indicates strong genetic differentiation among reefs and additionally supports the low dispersal capability of the larvae. The value obtained in the present study was nearly equal to the one previously found in the species (0.116 versus 0.118 in the present study) [29]. Our results indicated that the majority of variance in the population can be explained by the within reefs variation, indicating large variability among individuals. Although the variation among reefs within sites is smaller, reefs differ significantly in relation to the variability present at their respective site. Finally, variation among sites is high and significant as well. Therefore, the high diversity within each population should not be interpreted to conclude that populations differ only slightly. Higher differences between sites than between reefs within sites are consistent with a significant IBD pattern. Distance clearly acts as a barrier to gene flow in P. clavata, which is typical for species with short larval dispersal [33]. Moreover, isolation by distance is a phenomenon that occurs at local scales [65] and therefore it is evident among the closely located populations studied here. At larger spatial scales (hundreds to thousands of km), P. clavata also displayed a significant IBD pattern, similarly to another Mediterranean coral, C. rubrum [66]. However, at a fine spatial scale (cm to m), red gorgonian colonies did not show significant IBD, in contrast to C. rubrum [66].

Short distance migration (hundreds of meters) is likely to be the dominant scale of dispersal in P. clavata, but our results also suggest that migration from close undisturbed sites may be a significant source of recruits for disturbed areas to recover, since the maximum larval dispersal was between 20 to 60 km, with migrants detected in reefs separated by 20 km, but not by 60 km. This strengthens evidence from studies showing that coral population recovery after catastrophic mortality may be mainly supported by local migration from undamaged sites. A study of genetic connectivity in coral populations recovering after catastrophic bleaching revealed that the majority of detected migrants originated from the only site that was not decimated by a recent mortality event. Most of these immigrants were received by the site which reached pre-bleaching diversity, highlighting their role in population recovery [8]. Our results demonstrate that larvae may disperse from the Cinque Terre Marine Protected Area (Punta Mesco) to the adjacent populations, supporting their recovery after the disturbance. Additionally, the larval transport from deep, healthy reefs to the shallow ones, impacted by past mass mortality events, may have great importance for recovery after climatically induced population collapses.

Our findings indicate higher larval dispersal potential than reported in [19], which observed that larvae settle immediately on the substrate surrounding their mother colony. We found that transport of the larvae was not only common over distances of hundreds of meters, but also for tens of kilometers. The maximum migration distance detected in the present study reached 20 km, i.e. the distance that separates the two sites of La Spezia and Punta Mesco. Ten percent of the colonies from Punta Mesco were estimated as migrants from La Spezia. In this case, the larval dispersal was consistent with predominant currents, with the large scale Ligurian circulation being characterized by a cyclonic, east-to-west flow, active all year round and modulated by seasonality and wind forcing [67]. High contribution of migrants from Punta Mesco (8.5%) found in La Spezia may be caused by the predominantly southwesterly summer wind (Libeccio), which occasionally reverses coastal current direction [68]. In contrast to the previous two sites, Livorno was genetically homogeneous. This site is separated from La Spezia by 60 km of sandy bottom, unsuitable for the red gorgonian. The nearest known population to the south is located around 25 km from Livorno, at the deep rocky shoal below 40–50 m. It seems that both distance and depth differences may isolate Livorno reefs from other populations.

Theoretically, the high number of migrants exchanged between Punta Mesco and La Spezia should lead to panmixia in a few generations [69], but the genetic structure between these sites remains relatively high. These observations may be explained by an increased recruitment rate after the mortality. In stable populations the number of recruits is low, but after mass mortality the amount of available substrate increased and was occupied by new settlers, including migrants. When colonies recruited after the mass mortality events reached maturity and started reproducing, most of the available substrate was already occupied and only a minor proportion of their offspring could settle and survive, effectively reducing the opportunity for successful out-crossing. Additionally, P. clavata reaches maturity when colonies are from 3 years old (data from La Spezia population, [70]) to 7 years old (data from Medes Islands, Spain, [71]) and therefore only a few generations may have appeared after the mortality.

Our dispersal estimate (IBD regression slope of 0.037, R² = 0.78) was similar to the one reported previously for populations from the whole Mediterranean, when IBD slope was 0.033 (R² = 0.51) [29]. However, when comparing our results with a study that investigated a small spatial scale (IBD slope of 0.012, R² = 0.85 [31]), the isolation by distance appears stronger in our study, implying a lower dispersal. Variability in dispersal estimates may arise form larval responses to environmental conditions, modifying planktonic larval duration (PLD), or from differential establishment success after dispersal. Indeed, PLD of up to 25 days, estimated from laboratory experiments [18], indicate that larvae can delay their metamorphosis when lacking the necessary settlement stimuli [65] and, therefore, may be capable of long distance dispersal. The demographic state of the local receiver populations may strongly condition the success of establishment of long distance migrants. In mature populations, the rare arrival of few migrants may pass unnoticed, lost among the high recruitment mortality bottleneck that is typical of most populations with overlapping generations. Therefore large-scale dispersal is rarely detected as effective gene flow. However, where mass mortality has occurred, the opportunity of expansion in freely available habitat magnifies the probabilities of success of long distance dispersers. Moreover, the study of [31] examined not only populations that experienced mortality, but also one population that was recently founded. These processes, which have been named density-barrier effects or prior colonization effects, are most strikingly demonstrated for species which have been able to rapidly expand for thousands of km along novel available habitat (e.g., in invasive expansions or postglacial recolonizations) whereas in the ancient native ranges they remain highly structured even across a few tens of km [72; 73]. Therefore, inferring probability of connectivity from genetic data alone may be misleading, where the demographic conditions of the populations may prevail over dispersal capabilities in determining the possibilities of population recovery by migration in a metapopulation.

The larval exchange between deep and shallow reefs, such as that inferred in Punta Mesco and Livorno, has significant importance for the species' conservation. Colonies dwelling below 25–30 m were not affected by mass mortality events, in contrast to shallow subpopulations [22; 25]. Our research supports the hypothesis of [27; 28] that deeper subpopulations may supply larvae to shallower sites. However, depth differences in the present study were relatively small, aimed at representing nearby areas above and below the thermocline warming effects. Further studies on larval dispersal, including deeper subpopulations, are needed to develop a more complete picture of larval exchange between different depths. Deeper reefs may have an even greater larval contribution for the shallower and other deep reefs than what we estimated here, especially after disturbances. In the STRUCTURE analysis, some of the migrants exchanged between Punta Mesco and La Spezia were grouped into a separate cluster when k = 4. They might have originated from a different population, possibly from a greater depth at Punta Mesco, where P. clavata occurs down to 60 m. The evidence for longer larval dispersal, between Punta Mesco and La Spezia, suggests that damaged reefs may benefit from external larval sources. Unexpectedly high numbers of recruits in the La Spezia population were observed shortly after the mass mortality event [24]. In the subsequent years recruits density was five times higher than in the pre mortality period [33]. Our findings indicate that not only increased reproductive output might be responsible for the recovery of the damaged population, but also larval migration.

Conclusions

Our study failed to detect any genetic diversity loss in the P. clavata populations affected by mass mortality events. This may be due to the lack of test power and other confounding factors, including recent expansion or Wahlund effects. Our research confirmed low larval dispersal capability in the red gorgonian, since the maximum dispersal distance inferred from our data was between 20 and 60 km. However, this reduced ability for dispersal during the larval phase may still be ecologically significant for population replenishment and persistence, enabling migration between local populations. Population recovery after mortality events may be dependent on the possibility of propagule immigration from external sources, such as Protected Areas. Additionally, migration between reefs located at different depths implies that deeper refugia may provide larvae for shallow subpopulation recovery after climatically-induced mortality events affecting mostly shallow sites.

Supporting Information

[Figure omitted. See PDF.]

S1 Fig. The plot of the LnP(D).

https://doi.org/10.1371/journal.pone.0150590.s001

(PDF)

S2 Fig. The plot of Evanno deltaK.

https://doi.org/10.1371/journal.pone.0150590.s002

(PDF)

S1 Table. Microsatellite data set.

The allele sizes in 10 loci studied in 6 P. clavata populations in the Ligurian Sea.

https://doi.org/10.1371/journal.pone.0150590.s003

(XLSX)

Acknowledgments

The authors wish to thank C. Lombardi, A. Peirano, E. Mancuso, R. Albuquerque, F. Fernandes and J. Rodrigues for help during the fieldwork.

Author Contributions

Conceived and designed the experiments: SC HQ. Performed the experiments: JP. Analyzed the data: JP. Contributed reagents/materials/analysis tools: HQ. Wrote the paper: JP HQ SC. Contributed to sampling and data analysis: JB. Contributed to data analysis and reviewed the manuscript: ES.

Citation: 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. https://doi.org/10.1371/journal.pone.0150590

References

1. Jentsch A, Kreyling J, Beierkuhnlein C. A new generation of climate-change experiments: events, not trends. Front Ecol Environ. 2007;5(7): 365–374.

2. Giorgi F, Lionello P. Climate change projections for the Mediterranean region. Glob Planet Change. 2008; 63: 90–104.

3. Pujolar JM, Bevacqua D, Capoccioni F, Ciccotti E, De Leo GA, Zane L. No apparent genetic bottleneck in the demographically declining European eel using molecular genetics and forward-time simulations. Conserv Genet. 2011;12(3): 813–825.

4. Arnaud-Haond S, Duarte CM, Teixeira S, Massa SI, Terrados J, Tri NH, et al. Genetic recolonization of mangrove: genetic diversity still increasing in the Mekong Delta 30 years after Agent Orange. Mar Ecol Prog Ser. 2009;390: 129–135.

5. Pujolar JM, Vincenzi S, Zane L, Jesensek D, De Leo GA, Crivelli AJ. The effect of recurrent floods on genetic composition of marble trout populations. PLoS One. 2011;6(9): e23822. pmid:21931617

6. Colson I, Huges RN. Rapid recovery of genetic diversity of dogwhelk (Nucella lapillus L.) populations after local extinction and recolonization contradicts predictions from life-history characteristics. Mol Ecol. 2004;13(8): 2223–2233. pmid:15245396

7. Palumbi SR. Population genetics, demographic connectivity, and the design of marine reserves. Ecol Appl. 2003;13(1): 146–158.

8. Underwood JN, Smith LD, Van Oppen MJH, Gilmour JP. Multiple scales of genetic connectivity in a brooding coral on isolated reefs following catastrophic bleaching. Mol Ecol. 2007;16(4): 771–784. pmid:17284210

9. Hughes JM. Constraints on recovery: using molecular methods to study connectivity of aquatic biota in rivers and streams. Freshw Biol. 2007;52(4): 616–631.

10. Maggs CA, Castilho R, Foltz D, Henzler C, Jolly MT, Kelly J, et al. Evaluating signatures of glacial refugia for North Atlantic benthic marine taxa. Ecology. 2008;89: 108–122.

11. Bullock JM, Shea K, Skarpaas O. Measuring plant dispersal: an introduction to field methods and experimental design. Plant Ecol. 2006;186: 217–234.

12. Levin LA. Recent progress in understanding larval dispersal: new directions and digressions. Integr Comp Biol. 2006;46: 282–297. pmid:21672742

13. Gilg MR, Hilbish TJ. The geography of marine larval dispersal: coupling genetics with fine-scale physical oceanography. Ecology. 2003;84(11): 2989–2998.

14. Planes S, Jones GP, Thorrold SR. Larval dispersal connects fish populations in a network of marine protected areas. Proc Natl Acad Sci U S A. 2009;106(14): 5693–5697. pmid:19307588

15. Ballesteros E. Mediterranean coralligenous assemblages: A synthesis of the present knowledge. Ann Rev Mar Sci. 2006;44: 123–195.

16. Carpine C, Grasshof M. Les gorgonaires de la Méditerranée. Bull. Inst. Oceanogr. 1975;71: 1–140.

17. Öztürk B, Aktan Y, Topaloğlu B, Keskin Ç, Öztürk AA, Dede A, et al. Marine life of Turkey in the Aegean & Mediterranean Seas. Turkish Marine Research Foundation (TUDAV) Publications. Marine Education Series. Number 10. Istanbul, Turkey; 2004.

18. Linares C, Coma R, Mariani S, Diaz D, Hereu B, Zabala M. Early life history of the Mediterranean Gorgonian Paramuricea clavata: implications for population dynamics. Invertebr Biol. 2007;127(1): 1–11.

19. Coma R, Ribes M, Zabala M, Gili J-M. Reproduction and cycle of gonadal development in the Mediterranean gorgonian Paramuricea clavata. Mar Ecol Prog Ser. 1995;117: 173–183.

20. Linares C, Doak DF, Coma R, Díaz D, Zabala M. Life history and viability of a long-lived marine invertebrate: the octocoral Paramuricea clavata. Ecology. 2007;88(4): 918–928. pmid:17536708

21. Romano JC, Bensoussan N, Younes WA, Arlhac D. 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. 2000;323(4): 415–427.

22. Linares C, Coma R, Diaz D, Zabala M, Hereu B, Dantart L. Immediate and delayed effects of a mass mortality event on gorgonian population dynamics and benthic community structure in the NW Mediterranean Sea. Mar Ecol Prog Ser. 2005;305: 127–137.

23. Cerrano C, Arillo A, Azzini F, Calcinai B, et al. Gorgonian population recovery after a mass mortality event. Aquat Conserv. 2005;15: 147–157.

24. Cupido R, Cocito S, Sgorbini S, Bordone A, Santangelo G. Response of a gorgonian (Paramuricea clavata) population to mortality events: recovery or loss? Aquat Conserv. 2008;18 (6): 984–992.

25. Huete-Stauffer C, Vielmini I, Palma M, Navone A, Panzalis P, Vezzulli L, et al. Paramuricea clavata (Anthozoa, Octocorallia) loss in the Marine Protected Area of Tavolara (Sardinia, Italy) due to a mass mortality event. Mar Ecol. 2011;32: 107–116.

26. Coma R, Ribes M, Serrano E, Jiménez E, Salat J, Pascual J. Global warming-enhanced stratification and mass mortality events in the Mediterranean. Proc Natl Acad Sci U S A. 2009;106(15): 6176–6181. pmid:19332777

27. Cerrano C, Bavestrello G. Medium-term effects of die-off of rocky benthos in the Ligurian Sea. What can we learn from gorgonians? Chem Ecol. 2008;24 (S1): 73–82.

28. Cerrano C, Bavestrello G. Massive mortalities and extinctions. In: Marine Hard Bottom Communities. Patterns, Dynamics, Diversity, and Change. (Wahl M ed). Cap. 21. Ecological Studies 206: 295–307. Springer-Verlag Berlin Heidelberg; 2009

29. Mokhtar-Jamaï K, Pascual M, Ledoux JB, Coma R, Féral JP, Garrabou J, et al. From global to local genetic structuring in the red gorgonian Paramuricea clavata: the interplay between oceanographic conditions and limited larval dispersal. Mol Ecol. 2011;20(16): 3291–3385. pmid:21762434

30. Mokhtar-Jamaï K, Coma R, Wang J, Zuberer F, Féral JP, Aurelle D. Role of evolutionary and ecological factors in the reproductive success and the spatial genetic structure of the temperate gorgonian Paramuricea clavata. Ecol Evol. 2013;3(6): 1765–1779. pmid:23789084

31. Arizmendi-Mejía R, Linares C, Garrabou J, Antunes A, Ballesteros E, Cebrian E, et al. Combining Genetic and Demographic Data for the Conservation of a Mediterranean Marine Habitat-Forming Species. PLoS One. 2015;10(3): e0119585. pmid:25774522

32. Peirano A, Sgorbini S, Cupido R, Lombardi C, Cocito S. Decadal monitoring of corralligenous and bioconstruction organisms in the eastern Ligurian sea (NW Mediterraean). Proceedings of the 1st Mediterranean Symposium on the Coralligenous and other calcareous bio-concretions of the Mediterranean Sea 2009; 233–235.

33. Santangelo G, Cupido R, Cocito S, Bramanti L, Priori C, Erra F, Iannelli M. Effects of increased mortality on gorgonian corals (Cnidaria, Octocorallia): different demographic features may lead affected populations to unexpected recovery and new equilibrium points. Hydrobiologia. 2015; 1–17.

34. Cerrano C, Bavestrello G, Bianchi CN, Cattaneo-Vietti R, Bava S, Morganti C et al. A catastrophic mass-mortality episode of gorgonians and other organisms in the Ligurian Sea (North-western Mediterranean), summer 1999. Ecol Lett. 2000;3(4): 284–293.

35. Agell G, Rius M, Pascual M. Isolation and characterization of eight polymorphic microsatellite loci for the Mediterranean gorgonian Paramuricea clavata. Conserv Genet. 2009;10: 2025–2027.

36. Molecular Ecology Resources Primer Development Consortium, Aurelle D, Baker A et al. Permanent genetic resources added to the molecular ecology resources database 1 February 2010–31 March 2010. Mol Ecol Resour. 2010;10: 751–754. pmid:21565086

37. Toonen RJ, Hughes S. Increased Throughput for Fragment Analysis on ABI Prism 377 Automated Sequencer Using a Membrane Comb and STRand Software. Biotechniques. 2001; 31: 1320–1324. pmid:11768661

38. Van Oosterhout CV, Hutchinson WF, Wills DPM, Shipley P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4: 535–538.

39. Raymond M, Rousset F. GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered. 1995;86: 248–249.

40. Rousset F. Genepop'007: a complete reimplementation of the Genepop software for Windows and Linux. Mol Ecol Resour. 2008;8(1): 103–106. pmid:21585727

41. Nei M. Analysis of gene diversity in subdivided populations. Proc Natl Acad Sci U S A. 1973;70: 3321–3323. pmid:4519626

42. Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX 4.05, Logiciel Sous Windows TM Pour La Geʹneʹtique Des Populations. Laboratoire, Geʹnome, Populations, Interactions, CNRS UMR 5000, Universiteʹ de Montpellier II, Montpellier, France; 2004.

43. Kalinowski ST. HP-RARE 1.0: a computer program for performing rarefaction on measures of allelic richness. Mol Ecol Notes. 2005;5: 187–189.

44. PASS 14 Power Analysis and Sample Size Software. NCSS, LLC. Kaysville, Utah, USA, 2015. ncss.com/software/pass.

45. Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution (N Y). 1984;38: 1358–1370.

46. Luikart G, Cornuet JM. Empirical evaluation of a test for identifying recently bottlenecked populations from allele frequency data. Conserv Biol. 1998;12(1): 228–237.

47. Piry S, Luikart G, Cornuet JM. Bottleneck: a computer program for detecting recent reductions in the effective population size using allele frequency data. J Hered. 1999;90: 502–503.

48. Excoffier L, Laval G, Schneider S. Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol Bioinform Online. 2005;1: 47–50.

49. Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: dominant markers and null alleles. Mol Ecol Notes. 2007;7: 574–578. pmid:18784791

50. Pritchard J, Wen X, Falush D. Documentation for the STRUCTURE software, Version 2, Chicago 2007. http://pritch.bds.uchicago.edu.

51. Jakobsson M, Rosenberg NA. CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007;23: 1801–1806. pmid:17485429

52. Rosenberg NA. DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes. 2004:4: 137–138.

53. Rannala B, Mountain JL. Detecting immigration by using multilocus genotypes. Proc Natl Acad Sci U S A. 1997;94(17): 9197–9201. pmid:9256459

54. Piry S, Alapetite A, Cornuet J et al. GENECLASS2: a software for genetic assignment and first-generation migrant detection. J Hered. 2004;95: 536–539. pmid:15475402

55. Paetkau D, Slade R, Burden M, Estoup A. Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power. Mol Ecol. 2004;13: 55–65. pmid:14653788

56. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995;57: 289–300.

57. Costantini F, Fauvelot C, Abbiati M. Genetic structuring of the temperate gorgonian coral (Corallium rubrum) across the western Mediterranean Sea revealed by microsatellites and nuclear sequences. Mol Ecol. 2007;16(24): 5168–5182. pmid:17986195

58. Gutiérrez-Rodríguez C, Lasker HR. Microsatellite variation reveals high levels of genetic variability and population structure in the gorgonian coral Pseudopterogorgia elisabethae across the Bahamas. Mol Ecol. 2004;13: 2211–2221. pmid:15245395

59. Maier E, Tollrian R, Rinkevich B, Nurnberger B. Isolation by distance in the scleractinian coral Seriatopora hystrix from the Red Sea. Mar Biol. 2005;147: 1109–1120.

60. Magalon H, Adjeroud M, Veuille M. Patterns of genetic variation do not correlate with geographical distance in the reef-building coral Pocillopora meandrina in the South Pacific. Mol Ecol. 2005;14: 1861–1868. pmid:15910311

61. Ledoux JB, Mokhtar-Jamaï K, Roby C, Feral JP, Garrabou J, Aurelle D. Genetic survey of shallow populations of the Mediterranean red coral [Corallium rubrum (Linnaeus, 1758)]: new insights into evolutionary processes shaping nuclear diversity and implications for conservation. Mol Ecol. 2010;19(4): 675–690. pmid:20074314

62. Peery MZ, Kirby R, Reid BN, Stoelting R, Doucet-Beer E., Robinson S, et al. Reliability of genetic bottleneck tests for detecting recent population declines. Mol Ecol 2012;21(14): 3403–3418. pmid:22646281

63. Girod C, Vitalis R, Lebois R, Freville H. Inferring population decline and expansion from microsatellite data: a simulation-based evaluation of the MSvar methods. Genetics 2011;188: 165–179. pmid:21385729

64. Costantini F, Fauvelot C, Abbiati M. Fine-scale genetic structuring in Corallium rubrum: evidence of inbreeding and limited effective larval dispersal. Mar Ecol Prog Ser. 2007;340: 109–119.

65. Selkoe KA, Toonen RJ. Marine connectivity: a new look at pelagic larval duration and genetic metrics of dispersal. Mar Ecol Prog Ser. 2011;436: 291–305.

66. Ledoux JB, Aurelle D, Arizmendi-Mejía R, Frleta-Valić M, Linares C, Mokhtar-Jamaï K, Pralong M, Antunes A, Garrabou, J. Population genetics studies over contrasted spatial scales of two structural gorgonian species belonging to the coralligenous: state-of-the-art and conservation implications. 2ème Symposium Méditerranéen sur la conservation du Coralligène et autres Bio-Concrétions; 2014.

67. Astraldi M, Gasparini GP. The seasonal characteristics of the circulation in the North Mediterranean Basin and their relationship with the atmospheric-climatic conditions. J Geophys Res. 1992;97(C6): 9531–9540.

68. Haza AC, Özgökmen TM, Griffa A, Molcard A, Poulain PM, Peggion G. Transport properties in small-scale coastal flows: relative dispersion from VHF radar measurements in the Gulf of La Spezia. Ocean Dyn. 2010;60(4): 861–882.

69. Slatkin M. Gene flow in natural populations. Annu Rev Ecol Evol Syst. 1985;16: 393–430.

70. Cupido R, Cocito S, Manno V, Ferrando S, Peirano A, Iannelli M, et al. Sexual structure of a highly reproductive, recovering gorgonian population: quantifying reproductive output. Mar Ecol Prog Ser. 2012;469: 25.

71. Coma R, Zabala M, Gili J-M. Sexual reproductive effort in the Mediterranean gorgonian Paramuricea clavata. Mar Ecol Prog Ser. 1995;117: 185–192.

72. Neiva J, Pearson GA, Valero M, Serrao EA. Fine scale genetic breaks driven by historical range dynamics and ongoing density barriers-effects in the estuarine seaweed Fucus ceranoides L. BMC Evol Biol. 2012;12: 78. pmid:22672720

73. De Meester L, Gomez A, Okamura B, Schwenk K: The monopolization hypothesis and the dispersal-gene flow paradox in aquatic organisms. Acta Oecol. 2002;23(3): 121–135.