Impact of life history traits on gene flow: A multispecies systematic review across oceanographic barriers in the Mediterranean Sea

Abstract

Background

Marine species can demonstrate strong genetic differentiation and population structure despite the hypothesis of open seas and high connectivity. Some suggested drivers causing the genetic breaks are oceanographic barriers and the species’ biology. We assessed the relevance of seven major oceanographic fronts on species connectivity while considering their dispersal capacity and life strategy.

Methods

We systematically reviewed the scientific articles reporting population genetic differentiation along the Mediterranean Sea and across the Atlantic-Mediterranean transition. We retained those considering at least one sampling locality at each side of an oceanographic front, and at least two localities with no-front between them to correctly assess the effect of the front. To estimate the impact of life history characteristics affecting connectivity we considered the planktonic larval duration (PLD) and adult life strategy.

Results

Oceanographic barriers in the Mediterranean Sea seem to reduce gene flow globally; however, this effect is not homogeneous considering the life history traits of the species. The effect of the oceanographic fronts reduces gene flow in highly mobile species with PLD larger than 2–4 weeks. Benthic sessile species and/or with short PLD (< 2 weeks) have more significant genetic breaks between localities than species with higher motility; however, genetic differentiation occurs independently of the presence of a front.

Conclusion

Genetic connectivity is important for populations to recover from anthropogenic or natural impacts. We show that species with low mobility, mostly habitat-formers, have high genetic differentiation but low gene flow reduction mediated by the front, therefore, considering the importance of these species, we emphasize the vulnerability of the Mediterranean ecosystems and the necessity of protection strategies based on the whole ecosystem.

About the Authors:

Marta Pascual

* E-mail: [email protected]

Affiliation: Dept Genetics, Microbiology and Statistics/IrBio, Universitat de Barcelona, Diagonal 643, Barcelona, Spain

ORCID http://orcid.org/0000-0002-6189-0612

Borja Rives

Affiliation: Dept Genetics, Microbiology and Statistics/IrBio, Universitat de Barcelona, Diagonal 643, Barcelona, Spain

Celia Schunter

Affiliation: KAUST Environmental Epigenetic Program (KEEP), Division of Biological and Environmental Sciences & Engineering and Division of Applied Mathematics and Computer Sciences, King Abdullah University of Science and Technology, Thuwal, Kingdom of Saudi Arabia

Enrique Macpherson

Affiliation: Centre d’Estudis Avançats de Blanes (CEAB-CSIC), Car. Acc. Cala St. Francesc 14, Blanes, Girona, Spain

Introduction

The Mediterranean Sea displays one of the world’s richest diversity [1]. This basin contains more than 20,000 species of fish, cetaceans, invertebrates, sea turtles, algae and seaweeds, with a large proportion of endemism (ca. 20%). For this reason the Mediterranean Sea has been defined as a biodiversity hotspot for conservation priorities. As a consequence, there is an increasing interest to protect this unique natural heritage, and several Congress and World Summits have called on countries to establish a system of networks of marine protected areas (MPAs) with the aim of covering 20 to 30% of the total area [2]. Furthermore, an emphasis is laid on commercial species, where the management structures have not been sufficiently enforced [3].

Modelling studies have indicated the great importance of spatial configurations of MPAs and stock identification units to promote population and ecosystem persistence [4,5]. Although these configurations are clear in modelling results, efforts to assess and design MPAs and stock identifications are hindered by the lack of knowledge of several important factors. One of the major issues is the uncertainty about propagule dispersal which is one of the essential processes connecting areas and populations [6–8]. Currents and other oceanographic processes, which are often not being considered in management policies, potentially influence this dispersal capability and connectivity [9,10].

The majority of marine animals show life histories characterized by a long planktonic larval phase potentially allowing long-distance dispersal by marine currents. Therefore, we would expect to see low genetic structuring between localities of many marine species [11,12]. This could be particularly true in species with high fecundity or very large population sizes that potentially can have long-distance dispersal of eggs, larvae or adults [13]. However, a large number of studies in the last decade have disproven the concept that the seas are ‘open’ and well–connected and it has been shown that a number of species present a spatial genetic differentiation which is higher than expected if we only consider their dispersive abilities [14–16]. Several mechanisms may cause genetic differentiation between populations such as vicariance processes, caused by historical barriers, oceanographic currents, habitat discontinuities, local adaptation, larval behaviour, isolation by distance and limited dispersal capabilities [17]. Therefore, the global level of genetic differentiation within species results from a complex equilibrium between structuring factors (e.g. oceanographic fronts, isolation by distance) and homogenising factors (e.g. long larval pelagic phase, migratory behaviour of adults) [11].

The length of the pelagic phase or planktonic larval duration (PLD) could be considered as proxy of the species’ dispersal potential, being one of the most important homogenising factors in the population structure [18]. In the marine environment, dispersal can be validated using genetic markers as a measure of connectivity among localities [19]. Some studies have demonstrated that PLD and genetic metrics typically reflect scales of dispersal [20] but see the review of Selkoe et al.[9], whereas others have shown no clear patterns between genetic connectivity and dispersal capabilities [21].

In the Mediterranean Sea, oceanographic processes, such as current patterns and oceanographic discontinuities are crucial factors influencing population genetic connectivity [17,21,22]. Hence, to understand connectivity patterns on a large-scale it is important to consider the physical processes possibly influencing gene flow between localities. The Mediterranean Sea is an ideal study area for a survey incorporating oceanographic features and gene flow. The circulation patterns within the Mediterranean Sea is well described [23,24]. Moreover, several oceanographic discontinuities (Fig 1), mostly on the Spanish coast, originated by the entry of less saline Atlantic waters throughout the Gibraltar Strait have been reported to act as barriers to gene flow for numerous species [21,25,26]. The best-studied discontinuity is the Almeria-Oran Front (AOF), which has been proposed to be the main point of genetic break between the Atlantic Ocean and the Mediterranean Sea [27]. Other fronts have also shown to be strong barriers for genetic exchange in some species, such as the Balearic Front (BF) [21] or the Ibiza channel (IC) [28]. Other oceanographic processes occurring within the Mediterranean Sea, e.g. along the Sicily Channel, the Otranto Channel or the southern margin of the Aegean Sea, can also act as barriers to gene flow, however, their importance is not well studied [29–31].

[Figure omitted. See PDF.]

Fig 1. Map of the Mediterranean Sea with the name of the sub-basins, main currents (white lines) and oceanographic fronts analysed (red lines).

The name of the fronts and the acronym used (in black) is as follows: GS (Gibraltar Strait), AOF (Almeria-Oran Front), IC (Ibiza Channel), BF (Balearic Front), SC (Sicily Channel), ADR (Otranto Channel), AEG (southern margin of the Aegean Sea).

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

Recently, some studies have provided interesting models defining hydrodynamic provinces, by coupling Lagrangian simulations of particles with oceanographic currents to modelate the transport of larvae in order to characterize marine connectivity [32,33]. The main aim of these models was to identify natural conservation units to be used in the establishment of an MPA network. These units would be defined as the Cells of Ecosystem Functioning (CEF), based on both oceanographic and ecological processes, arranged in space so as to account for both patterns (biodiversity distribution) and processes (ecosystem functioning) [34]. However, these CEFs do not include data on gene flow among localities hence lacking an important factor in connectivity processes. A review by Patarnello et al. [27] was the first attempt to summarize the genetic data regarding the biogeographical separation between the Mediterranean and Atlantic biota. This review (including 20 studies) showed steep changes of allele frequencies associated with the Almeria-Oran Front (AOF) but failed to relate biological traits with genetic differentiation. After this review numerous genetic studies have been carried out in different areas of the Mediterranean, most of them analysing one or more oceanographic fronts. Some of the studies, however, do not provide data comparing localities separated by fronts and control sites to correctly evaluate the influence of the front. Other studies contain genetic differences between sampling localities that include more than one front between them making it difficult to draw adequate conclusions on the influence of the oceanographic fronts on gene flow. Hence, there is a lack of a comprehensive analysis on the effect of oceanographic fronts on gene flow in the Mediterranean Sea, which can seriously compromise our knowledge on the connectivity patterns among localities and in the establishment of MPA networks [34].

In the present review we aim to assess the effect of life history characteristics on population differentiation considering the impact of oceanographic discontinuities. In order to correctly evaluate whether a front is the cause for genetic differences between localities we screened all population genetic studies based in the Mediterranean Sea considering at least one sampling locality at each side of the oceanographic front, and at least two localities with no-front between them. We focus our review on the seven most important Mediterranean oceanographic discontinuities: Gibraltar Strait (GS), Almeria-Oran Front (AOF), Ibiza Channel (IC), Balearic Front (BF), Sicily Channel (SC), Otranto Channel (ADR) and the southern margin of the Aegean Sea (AEG) [23,35,36]. Our main goal is to evaluate the impact of life history characteristics on connectivity reduction, in both the absence and presence of fronts, so we categorized all species based on their larval dispersal capabilities and their adult behaviour. Finally, we discuss the potential implications of life history traits and oceanographic discontinuities in the establishment of management units and networks of MPAs.

Materials and methods

We searched the ISI Web of Knowledge database for scientific articles dated until 2016 which evaluated population genetic differentiation along the Mediterranean Sea and across the Atlantic-Mediterranean transition. Keywords were selected to identify these studies: “Mediterranean Sea”, “gene flow”, “genetic/population structure”, “genetic differentiation”, “connectivity”, which resulted in a total of 718 studies (Fig 2). We assessed 440 papers but retained only 72 that could be used to correctly evaluate the impact of oceanographic fronts (Fig 2 and S1 Checklist). We retained only those studies meeting the following criteria: reporting the analysis of genetic differentiation between localities situated at both sides of one front and including genetic data from localities not separated by a front that could be used as control, and also based on molecular marker selection see below. A comparison was also excluded whenever more than one front was present between two localities. Control localities are necessary to correctly evaluate the effect of a front as they allow discriminating between reduced connectivity due to the species limited dispersal capabilities or to the effect of the oceanographic front [21]. Invasive species were not included in the analyses to avoid confounding effects due to genetic differentiation between localities driven by passive colonization, sometimes from different native areas and highly influenced by genetic drift during arrival of colonisers [37–39].

[Figure omitted. See PDF.]

Fig 2. PRISMA flow diagram for literature search.

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

Molecular markers selection

Different molecular markers have been suggested to be more adequate at identifying historical (e.g. mtDNA) or contemporary (e.g. microsatellites) processes [40]. We therefore classified marker types into two groups: Mitochondrial genes, including COI, 16S, cyt b were categorized as “MT”, and nuclear polymorphic genes including allozymes, microsatellites and SNPs under the category “NUC”. For some species more than one study or analysis with several different molecular markers and genetic distance measurements were available. To avoid biases due to overrepresentation of the same species in each front, only one analysis per species per front was considered. The following criteria were applied: (1) for studies applying nuclear and mitochondrial markers the former were preferentially chosen; (2) when several indices were reported we preferentially retained the data for F_ST values between localities against other less frequently used indices.

Oceanographic discontinuities

Between the Atlantic Ocean and Mediterranean Sea there are 7 major fronts frequently analysed in population genetic and biogeographic studies (see Fig 1): Gibraltar Strait (GS), Almeria-Oran (AOF), Ibiza Channel (IC), Balearic Front (BF), Sicily Channel (SC), Otranto Channel (ADR), and the southern margin of the Aegean Sea (AEG). The inflow of the North Atlantic Central Water (NACW) throughout the Gibraltar Strait (GS) is the most important oceanographic process between the Atlantic Ocean and Mediterranean Sea [36,41]. The GS discontinuity is located just before the entry of the Atlantic waters throughout the Gibraltar Strait, when the NACW recirculates near the Strait, in front of Cape Trafalgar towards the northwest along the coast of Cadiz [42]. Once in the Mediterranean Sea the NACW encounters the higher density Mediterranean water generating the Almeria-Oran front (AOF), a quasi-continuous front where its northern end detaches from the Spanish coast between Almeria and Cartagena and its southern end terminates around Oran on the North African coast [43]. The NACW is modified, increasing its salinity, into a mass usually called Modified Atlantic Water (MAW) that reaches the north-western Mediterranean basin and is deflected eastward by the cyclonic circulation around the Balearic Islands forming a well-defined second density front, the Balearic Front (BF). Just south of this front is the Ibiza Channel (IC), 80 Km width and 800 m depth, corresponding to the passage intersecting the Balearic topographic ridge between Ibiza and the Iberian Peninsula at Cape La Nao (Fig 1). The MAW also flows along the Algerian coast and crosses the Sicily Channel (SC), which divides the Western and Eastern basins of the Mediterranean Sea [35,36]. The next front considered going eastwards is the Otranto Channel (ADR) located at the entrance of the Adriatic Sea whose sill is 800 m deep [44]. Finally, around the southern margin of the Aegean Sea (AEG) another oceanographic discontinuity is observed formed by several cyclonic, anticyclonic gyres and eddies interconnected by currents and jets flowing at speeds of 20–30 cm s^-1 [45].

Life history data

To evaluate the impact of life history characteristics affecting connectivity we considered the planktonic larval duration (PLD), as a proxy of larval dispersal, and adult life behaviour. Each species was assigned to one of the following three groups according to their PLD: 1–15 days (S), 16–30 days (M), ≥31 days (L). Another categorical variable was constructed according to adult life strategy (LIFE): species were considered to be benthic sessile or with limited motility (BS), benthic vagile (BM), or pelagic (PEL). Finally, we also categorized the species according to major taxonomic groups: Angiosperm, Porifera, Cnidaria, Echinodermata, Mollusca, Crustacea, Tunicata and Pisces and grouped into a phylogenetic tree [46]. Biological information was obtained from the literature for angiosperm [47], porifera [48], cnidaria [49], echinodermata [50], mollusca [51–53], crustacea [54], tunicata [55] and pisces [56,57] (see also references cited in these articles).

Statistical analyses

Pairwise genetic population distances (preferentially F_ST), comparing localities separated by fronts and one control pairwise comparison representative of no-fronts, were extracted from the selected scientific articles. Geographic distances between locations were approximated following the coastline with Google Earth.

For different life history categories (PLD and LIFE strategy) we tested differences in frequency of significant genetic distances between localities, separated by fronts or control sites, by chi-square tests [58].

In order to avoid differences in significance due to isolation by distance, for each selected article all pairwise genetic distances (F_ST) were also standardized by geographic distance in Km as in Galarza et al.[21]. We considered the existence of connectivity reduction to be mediated by the front when the standardized F_ST between localities separated by front was larger than between control sites. We compared the two standardized pairwise F_ST (with front and no-front) with a Wilcoxon matched pairs test for each variable (FRONT, PLD, LIFE strategy and MARKER) using the programme STATISTICA V 8 [59].

To analyse the interaction between the front effect and the other variables, we created a continuous variable ranging between 0 and 1 called NDIF. This variable was calculated as (x-xmin) / (xmax-xmin) were x represents the difference between front and no-front standardized F_ST values, xmax is the highest difference value and xmin the smallest one. We performed permutational multivariate analysis of variance (PERMANOVA) using the statistical package PRIMER-E v6 [60]. We used NDIF as dependent variable and PLD (S, M, L), LIFE strategy (BS, BM, PEL), MARKER (MT, NUC) and FRONT (Gibraltar Strait (GS), Almeria-Oran (AOF), Ibiza Channel (IC), Balearic Front (BF), Sicily Channel (SC), Otranto Channel (ADR), and the southern margin of the Aegean Sea (AEG)) as fix factors. PLD and LIFE strategy were not analysed together because not all levels were represented in both predictive variables. The assayed interactions included FRONT×PLD, FRONT×LIFE, and FRONT×MARKER. For each analysis 999 permutations were carried out.

Results

In total 72 scientific papers passed our requirements to assess the effect of the fronts along the Mediterranean Sea and across the Atlantic-Mediterranean transition. These papers measured genetic differentiation between localities influenced by a front as well as a control comparison between localities without no-front (S1 Table). The mean geographic distance (±SE) between the selected control locations was 305 Km (± 13.3 Km). Fronts were evaluated selecting two localities on each side of the front with a mean distance of 492 Km ± 25.2 Km. Overall, we had information of 176 datasets of 70 species analysing different fronts (S2 Table). Fishes were the most represented taxonomic group (41.1%), followed by crustaceans (21.0%) and molluscs (13.6%). The mean (±SE) number of fronts correctly analysed per species in each reference was 2.12 ± 0.13.

Genetic data and life history traits

We observed that 66% of comparisons involving species with short PLD (S) had significant differences between localities, while only 41% and 28% were significant for species with medium (M) and long (L) PLD, respectively (Fig 3A). This difference was observed in sites separated by a front (Chi: 16.53, p = 0.0003) and also in control sites not separated by any known front (Chi: 18.66, p = 0.0001). A similar trend was observed considering the LIFE strategy variable, where the percentage of significant comparisons reduced with increased mobility (Fig 4A). The percentage of benthic sessile (BS) species showing significant genetic differences was always higher than in species with higher motility capacity (BM and PEL). However, no significant differences between categories were observed for locations separated by a front (Chi: 4.33, p = 0.1146) or without front (Chi: 2.78, p = 0.2490).

[Figure omitted. See PDF.]

Fig 3. Influence of Planktonic larval duration (PLD) categories on genetic connectivity.

(A) Frequency of significant genetic differences (P<0.05) between localities separated by a front and without a front. (B) Frequency of comparisons showing genetic reduction mediated by the front. PLD categories are identified as S = 1–15 days, M = 16–30 days and L ≥ 31 days.

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

[Figure omitted. See PDF.]

Fig 4. Influence of adult LIFE strategy categories on genetic connectivity.

(A) Frequency of significant genetic differences (P<0.05) between localities separated by a front and without a front. (B) Frequency of comparisons showing genetic reduction mediated by the front. LIFE strategy categories are identified as BS = Benthic sessile or limited motility, BM = benthic vagile and PEL = pelagic.

https://doi.org/10.1371/journal.pone.0176419.g004

Effects of oceanographic fronts

We compared the standardized F_ST distances between localities separated by the front (F) and those without front (NF) and considered gene flow reduction when the former distance was larger. The frequency of comparisons with reduction increased with PLD (Fig 3B), being 49% for species with small PLD and 66% for species with large PLD. The same tendency was observed for adult LIFE strategy (Fig 4B), with benthic sessile or limited motility species showing the lowest number of comparisons with front reduction (47%) and pelagic species the largest (82%).

We also evaluated the impact of life history traits (PLD, LIFE strategy and the two of them combined to integrate the whole life-history) on the connectivity reduction mediated by fronts using Wilcoxon matched pairs tests (Table 1). Interestingly, oceanographic discontinuities significantly reduced gene flow in species with L PLD, with PEL LIFE strategy, as well as in species with BMM and PELL categories combining both adult and larval mobility capacities (Table 1). No significant reduction mediated by the front was observed with Wilcoxon matched pairs test for different types of markers (Table 1). Thus, although presenting higher genetic differentiation between localities (Figs 3A and 4A), species with low dispersal capabilities seemed less affected by the fronts (Figs 3B and 4B). On the contrary, the connectivity of species with larger motility capacity would be more affected by oceanographic discontinuities. We detected gene flow reduction for species with long PLDs (more than two weeks) in most Mediterranean fronts (Fig 5 and S3 Table). Furthermore, higher frequency of significant genetic differentiation in short PLD species and reduced adult mobility across discontinuities was observed across most Mediterranean oceanographic fronts (Fig 5 and S3 Table).

[Figure omitted. See PDF.]

Fig 5. Genetic differentiation patterns between localities separated by oceanographic fronts.

Frequency of species, for each PLD category (A) and LIFE strategy (B), showing significant genetic distances (black) and gene flow reduction (blue) between localities separated by oceanographic fronts. Each oceanographic front (red lines) is analysed separately: Gibraltar Strait (GS), Almeria-Oran Front (AOF), Ibiza Channel (IC), Balearic Front (BF), Sicily Channel (SC), Otranto Channel (ADR), and the southern margin of the Aegean Sea (AEG).

https://doi.org/10.1371/journal.pone.0176419.g005

[Figure omitted. See PDF.]

Table 1. Wilcoxon matched pairs test comparing the standardized F_ST distances between localities separated by fronts (F) and those without front (NF) for different larval and adult stage categories, and type of molecular marker.

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

Life history traits such as the length of PLD or the adult habitat show similar patterns across different taxa that seem to relate with genetic differentiation independent of phylogenetic relationships (Fig 6). For instance, species with low dispersal capabilities, in all taxa, frequently present significant differentiation between populations separated by fronts but not related to the front reducing connectivity (Fig 6).

[Figure omitted. See PDF.]

Fig 6. Phylogenetic tree of analysed taxa representing life history traits and genetic differentiation patterns.

In the Life History column the number of species is given for each taxon, with the inner circle representing PLD categories (S = Orange, M = Brown, L = Blue) and the outer circle representing LIFE strategy categories (BS = Orange, BM = Brown, PEL = Blue). In the Genetic Differentiation column the number of front comparisons is given for each taxon, with the inner circle representing the frequency of comparisons showing genetic reduction patterns (YES = dark grey, NO = light grey), and the outer circle representing the frequency of significant genetic differences between localities separated by a front (YES = dark grey, NO = light grey).

https://doi.org/10.1371/journal.pone.0176419.g006

For each dataset we also compared with a Wilcoxon matched pairs test the standardized F_ST distances between localities separated by the front (F) and those without front (NF). No significant differences were observed for each front analysed separately or all together, despite having more species with higher standardized differentiation between populations when separated by a front (Table 2). The effect of different factors and their interactions (FRONT, PLD, LIFE and MARKER) on gene flow reduction was assessed with PERMANOVA using the standardized NDIF as dependent continuous variable. Only PLD showed a significant effect (Table 3). The differentiation between localities was not dependent on the front assayed nor the marker used (Table 3). Unfortunately, the variable LIFE_PLD combined could not be compared for each front since increasing the number of categories drastically reduces the number of comparisons to be assessed. More population genetic studies in a wider number of species with different dispersal capabilities are needed to statistically analyse this variable.

[Figure omitted. See PDF.]

Table 2. Wilcoxon matched pairs test for each oceanographic front comparing the standardized F_ST values between localities separated by the front (F) and those without front (NF).

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

[Figure omitted. See PDF.]

Table 3. PERMANOVA analyses considering the effect of different variables and their interaction on gene flow reduction.

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

Discussion

The results of the present study demonstrate that the oceanographic discontinuities in the Mediterranean Sea, in general, affect gene flow among localities. However, this effect is not homogeneous considering the life history trait of species (Fig 7), not restricted to phylogenetic groups. Species with short dispersal capabilities showed significant genetic differentiation between most locations, but the gene flow reduction was not mediated by the presence of a front. However, significant connectivity reduction due to an oceanographic front was detected in species with PLD longer than 2–4 weeks and adult mobility. The variability in this effect can have important consequences in the analysis of the effective connectivity processes in the Mediterranean Sea [32,33] and in the identification of natural conservation units for the establishment of an MPA network or in conservation strategies [8,34].

[Figure omitted. See PDF.]

Fig 7. Effect of life history traits on population genetic differentiation and their impact on connectivity reduction mediated by oceanographic discontinuities.

These effects influence the capacity of population recovery and species vulnerability.

https://doi.org/10.1371/journal.pone.0176419.g007

Life history traits and gene flow reduction

It is now more widely accepted that marine species can have strong genetic population structuring even on very small scales dismissing the previous hypothesis of open, homogeneous and interconnected seas. This strong population structure can be caused by oceanographic barriers such as the Almeria-Oran Front [21,27,61], or life history traits such as larval behaviour [62] and Isolation By Distance [63]. We show that genetic differences among localities are highly dependent on life history traits of species. The highest significant differences occur in sessile benthic species with short planktonic larval durations (PLD < two weeks) independently of the presence of a front.

PLD has a strong influence on the genetic structure but also on the geographical distribution of species [64] with adult traits being determinant as well [65]. The larvae of species with a short PLD in the Mediterranean Sea usually remain close to the coastline [57,66] and disperse only a few meters/kilometres from the parents, as seen for example in sponges or gorgonians [62,67]. Therefore, the population differentiation in these species seems mostly caused by the reduced dispersal capabilities of the larvae since their adults are sessile. For the sponge Crambe crambe the estimated mean dispersal distances per generation were only about 35 cm, suggesting that the observed fine-scale genetic structure may be common in invertebrates with lecitotrophic larvae [14]. A similar result was observed in the coral Astroides calycularis [68], where the low connectivity is explained by the negative buoyancy and demersal behaviour of the larvae. In the gorgonian Paramuricea clavata high levels of self-recruitment and parentage relationships were detected at a small scale [69] with genetic drift having a strong impact in populations [70]. No isolation by distance and high genetic differentiation between localities separated by dozens to hundreds of kilometres was detected in the seagrass Zostera noltei, where genetic and physical connectivity assessment also indicated that rare long distance dispersal was possible [71]. These characteristics can generate strong genetic differentiation among geographically closed localities and similarities between widely separated localities, independently of the presence or absence of an oceanographic discontinuity.

For most species with longer PLDs (> two weeks) larvae move along the continental shelf and slope [66]. The distribution of their larvae are strongly affected by currents [72,73], eddies [74] and oceanographic fronts [75]. In these species, genetic differentiation between control locations (no-front, NF) was generally smaller than between localities separated by oceanographic discontinuities. This difference indicates that the front has an additional effect reducing gene flow among localities, although genetic differentiation between localities was seldom significant. Therefore, most fish and crustacean species living on the shelf and slope (e.g. Diplodus vulgaris, Mullus spp., Serranus cabrilla, Liocarcinus depurator) present a genetic structure generated by the presence of fronts [17,21,25,28]. However, as several authors have pointed out, limited larval dispersal can also be observed in species with long PLD (>2 weeks) [6,76,77]. Habitat suitability may have important implication; in the fish Tripterygion delaisi the presence of continuous rocky habitat between localities prevents genetic differentiation, while large discontinuities of sand or deep-water channels seem to reduce gene flow [78]. The largest distance between recruitment locations in this species measured through kinship analysis was ca. 11.5 km [79]. These direct measurements on population connectivity can provide compelling evidence to estimate the size and distance among areas in order to link them into an ecologically coherent MPAs network.

Gene flow reduction by oceanographic front

The oceanographic discontinuities in the Mediterranean Sea restrict gene flow in numerous species. However, the dynamic behaviour of these fronts shows significant intra- and inter-annual variability [80,81]. When all fronts were analysed simultaneously, we observed gene flow reduction in 58.8% of the comparisons. This tendency seems to apply to all fronts with the smallest reduction observed in the Sicily Channel (SC) and the largest in the Gibraltar Strait (GS). Differences between fronts could be due to the species studied and to fluctuations in the front strength. The AOF has been described as the main oceanographic barrier causing genetic differentiation along the Atlantic-Mediterranean transition area [27,82]. However, this front presents inter-annual variability allowing some degree of homogenization between localities of the same species at certain years [22]. Therefore, although most fronts caused a clear restriction in gene flow, the effect of oceanographic barriers might be temporarily relaxed and different results obtained regarding the year of dispersal of the analysed samples [22,83]. The reproductive period of species may also determine whether fronts are the cause of the genetic structuring. Across the Ibiza Channel, localities will be genetically differentiated if the reproduction coincides with the largest intensity of the front, as shown for Serranus cabrilla [17], and not differentiated if the reproduction period coincides with the lowest intensity of the front, as observed in Epinephelus marginatus [84].

Potential effect in the establishment of a MPA network and stock identifications in the Mediterranean Sea

The identification of population units and the definition of population boundaries are one of the highest priorities in the management of marine ecosystems [85,86]. Studies on the connectivity assessment between MPA networks, with an integrated multispecies approach encompassing oceanography, dispersal capabilities and genetics, are strongly needed [34,87]. In order to implement functional networks of conservation a number of umbrella species representatives of the ecosystem should be selected and analysed with genetic tools [88]. Therefore, we recommend that species with different life histories to be selected within the umbrella species group as they provide different, but complementary information on connectivity of the whole ecosystem. The role of oceanographic fronts in the restriction of gene flow in numerous species supports the eco-regionalization of the Mediterranean Sea proposed by different studies [32,33]. This regionalization based on dispersal patterns and connectivity mostly affect species with longer PLDs, a common characteristic of most commercial species (e.g. fishes and crustaceans), and any protection measures, including the establishment of MPAs networks or management plans, should consider them.

However, the establishment of MPA networks must contemplate the whole ecosystem, with special caution to those species with low dispersal capabilities. Oceanographic discontinuities did not generally affect the connectivity patterns in species with low dispersal larval capabilities (PLD < 2 weeks). These species include numerous sponges, gorgonians, angiosperms that are habitat formers [47,89] and essential in providing refuges or food for settlers of fishes [90]. The so-called “ecosystem engineers”, i.e. species that generate new habitat for other species, need an adequate protection since their elimination will be hard to recover given their low dispersal potential capabilities [91]. In the Mediterranean Sea destructive fishing practices, eutrophication and coastal development are among the major impacts responsible for habitat change [92] and in numerous ecosystems the “engineers species” have disappeared [93,94]. Furthermore, reducing cumulative local human impacts cannot reverse the loss of natural capital and the recovery of these species only occurs by transplanting propagules or juvenile stages [95] since most algae, sponges and other “engineers species” disperse only few meters and, in general, less than 1 km [91]. Moreover, local adaptation processes should also be considered in restoration actions given the importance of genetic adaptation to environmental pressures [96,97].

An optimal MPA network should protect areas of adequate size, including different home ranges, to help engineer species survive and flourish, thus providing habitat and food for larger species such as commercial fish species. Our results confirm the study of Shanks et al.[18] that suggested for marine reserve networks to be designed large enough (ca. 4–6 km in diameter) to contain short-distance dispersing larvae and be spaced far enough (ca. 10–20 km) for long-distance dispersing propagules released from one reserve to be able to settle in adjacent reserves [8]. These protected areas should be further interconnected providing suitable settlement habitat for far travelling larvae [98]. Finally, these networks need be created on either side of oceanographic fronts, since they generally are barriers to gene flow and often determine boundaries between hydrodynamic provinces [33], ensuring resilience of the whole ecosystem.

Conclusions

Genetic connectivity is important for recovery from anthropogenic of natural impacts. We show that both larval and adult life characteristics impact connectivity among localities and across oceanographic discontinuities. We observed that species with lower mobility potential have higher frequency of pairwise localities showing significant genetic differentiation than species with higher motility capacity, but independent of the existence of fronts and more subject to genetic drift. Moreover, oceanographic discontinuities reduce gene flow in species with medium to high dispersal abilities (Fig 7). We encourage, for the resilience of the whole ecosystem, to consider oceanographic discontinuities and umbrella species with different life characteristics when identifying management units and for designing networks of Marine Protected Areas.

Supporting information

[Figure omitted. See PDF.]

S1 Checklist. PRISMA Checklist for current study.

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

(PDF)

S1 Table. List of studies and species considered in the present study for the seven Mediterranean oceanographic discontinuities showing significant genetic distances and gene flow reduction by oceanographic front.

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

(PDF)

S2 Table. Genetic and geographic distances between localities separated by fronts and no-fronts for the 177 datasets.

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

(PDF)

S3 Table. Number of species, for each PLD category and LIFE strategy, showing Non significant/Significant genetic distances and No reduction/ Reduction gene flow between localities separated by each oceanographic front.

https://doi.org/10.1371/journal.pone.0176419.s004

(PDF)

Acknowledgments

We thank Maria Garcia and Javier Macpherson for help on drawing Fig 7, Enric Ballesteros for the photographs, and Ferran Palero for his help in the design of Fig 6.

Author Contributions

1. Conceptualization: MP CS EM.

2. Data curation: EM BR.

3. Formal analysis: MP BR EM.

4. Funding acquisition: MP EM.

5. Investigation: MP BR CS EM.

6. Methodology: MP CS EM.

7. Resources: EM CS BR.

8. Software: MP BR EM.

9. Supervision: MP EM.

10. Writing – original draft: MP BR CS EM.

11. Writing – review & editing: MP BR CS EM.

Citation: Pascual M, Rives B, Schunter C, Macpherson E (2017) Impact of life history traits on gene flow: A multispecies systematic review across oceanographic barriers in the Mediterranean Sea. PLoS ONE 12(5): e0176419. https://doi.org/10.1371/journal.pone.0176419

References

1. Coll M, Piroddi C, Steenbeek J, Kaschner K, Ben Rais Lasram F, Aguzzi J, et al. The Biodiversity of the Mediterranean Sea: Estimates, Patterns, and Threats. PLoS One. Public Library of Science; 2010;5: e11842. pmid:20689844

2. Olsen EM, Fonseca L, Goñi R, Johnson D, Katsanevakis S, Macpherson E, et al. Achieving Ecologically Coherent MPA Networks in Europe: Science Needs and Priorities. Marine Board Position Paper 18. Larkin K, McDonough N, editors. Ostend, Belgium: European Marine Board; 2013.

3. Lleonart J, Maynou F. Fish stock assessments in the Mediterranean: state of the art. Sci Mar. 2003;67: 37–49.

4. Kaplan DM. Fish life histories and marine protected areas: an odd couple? Mar Ecol Prog Ser. 2009;377: 213–225.

5. Bell JJ. Connectivity between island Marine Protected Areas and the mainland. Biol Conserv. 2008;141: 2807–2820.

6. Almany GR, Hamilton RJ, Bode M, Matawai M, Potuku T, Saenz-Agudelo P, et al. Dispersal of Grouper Larvae Drives Local Resource Sharing in a Coral Reef Fishery. Curr Biol. 2013;23: 626–630. pmid:23541728

7. Magris RA, Pressey RL, Weeks R, Ban NC. Integrating connectivity and climate change into marine conservation planning. Biological Conservation. 2014. pp. 207–221.

8. Dubois M, Rossi V, Ser-Giacomi E, Arnaud-Haond S, López C, Hernández-García E. Linking basin-scale connectivity, oceanography and population dynamics for the conservation and management of marine ecosystems. Glob Ecol Biogeogr. 2016;25: 503–515.

9. Selkoe KA, Gaggiotti OE, ToBo Laboratory, Bowen BW, Toonen RJ. Emergent patterns of population genetic structure for a coral reef community. Mol Ecol. 2014;23: 3064–3079. pmid:24866831

10. White C, Selkoe KA, Watson J, Siegel DA, Zacherl DC, Toonen RJ. Ocean currents help explain population genetic structure. Proc R Soc B. 2010;277: 1685–94. pmid:20133354

11. Charrier G, Chenel T, Durand JD, Girard M, Quiniou L, Laroche J. Discrepancies in phylogeographical patterns of two European anglerfishes (Lophius budegassa and Lophius piscatorius). Mol Phylogenet Evol. 2006;38: 742–754. pmid:16309924

12. Cordero D, Peña JB, Saavedra C. Phylogeographic analysis of introns and mitochondrial DNA in the clam Ruditapes decussatus uncovers the effects of Pleistocene glaciations and endogenous barriers to gene flow. Mol Phylogenet Evol. 2014;71: 274–287. pmid:24269315

13. Palumbi SR. Genetic divergence, reproductive isolation, and marine speciation. Annu Rev Ecol Syst. Annual Reviews 4139 El Camino Way, P.O. Box 10139, Palo Alto, CA 94303–0139, USA; 1994;25: 547–572.

14. Calderón I, Ortega N, Duran S, Becerro M, Pascual M, Turon X. Finding the relevant scale: clonality and genetic structure in a marine invertebrate (Crambe crambe, Porifera). Mol Ecol. 2007;16: 1799–1810. pmid:17444893

15. Palero F, Abelló P, Macpherson E, Gristina M, Pascual M, Abello P, et al. Phylogeography of the European spiny lobster (Palinurus elephas): Influence of current oceanographical features and historical processes. Mol Phylogenet Evol. 2008;48: 708–717. pmid:18515152

16. Pérez-Ruzafa Á, González-Wangüemert M, Lenfant P, Marcos C, García-Charton JA. Effects of fishing protection on the genetic structure of fish populations. Biol Conserv. 2006;129: 244–255.

17. Schunter C, Carreras-Carbonell J, MacPherson E, Tintoré J, Vidal-Vijande E, Pascual A, et al. Matching genetics with oceanography: Directional gene flow in a Mediterranean fish species. Mol Ecol. 2011;20: 5167–5181. pmid:22097887

18. Shanks AL, Grantham BA, Carr MH. Propagule dispersal distance and the size and spacing of marine reserves. Ecol Appl. Ecological Society of America; 2003;13: 159–169.

19. Selkoe KA, D’Aloia CC, Crandall ED, Iacchei M, Liggins L, Puritz JB, et al. A decade of seascape genetics: Contributions to basic and applied marine connectivity. Marine Ecology Progress Series. 2016. pp. 1–19.

20. 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.

21. Galarza JA, Carreras-Carbonell J, Macpherson E, Pascual M, Roques S, Turner GF, et al. The influence of oceanographic fronts and early-life-history traits on connectivity among littoral fish species. Proc Natl Acad Sci U S A. 2009;106: 1473–1478. pmid:19164518

22. Pascual M, Palero F, García-Merchán VH, Macpherson E, Robainas-Barcia A, Mestres F, et al. Temporal and spatial genetic differentiation in the crab Liocarcinus depurator across the Atlantic-Mediterranean transition. Sci Rep. Nature Publishing Group; 2016;6: 29892. pmid:27431989

23. Fernández V, Dietrich DE, Haney RL, Tintoré J. Mesoscale, seasonal and interannual variability in the Mediterranean Sea using a numerical ocean model. Prog Oceanogr. 2005;66: 321–340.

24. Rio MH, Poulain PM, Pascual A, Mauri E, Larnicol G, Santoleri R. A Mean Dynamic Topography of the Mediterranean Sea computed from altimetric data, in-situ measurements and a general circulation model. J Mar Syst. Liege, BELGIUM; 2007;65: 484–508.

25. Galarza JA, Turner GF, Macpherson E, Rico C. Patterns of genetic differentiation between two co-occurring demersal species: the red mullet (Mullus barbatus) and the striped red mullet (Mullus surmuletus). Can J Fish Aquat Sci. 2009;66: 1478–1490.

26. Marie AD, Lejeusne C, Karapatsiou E, Cuesta JA, Drake P, Macpherson E, et al. Implications for management and conservation of the population genetic structure of the wedge clam Donax trunculus across two biogeographic boundaries. Sci Rep. 2016;6.

27. Patarnello T, Volckaert F, Castilho R. Pillars of Hercules: is the Atlantic-Mediterranean transition a phylogeographical break? Mol Ecol. 2007;16: 4426–4444. pmid:17908222

28. García-Merchán VH, Robainas-Barcia A, Abelló P, Macpherson E, Palero F, García-Rodríguez M, et al. Phylogeographic patterns of decapod crustaceans at the Atlantic–Mediterranean transition. Mol Phylogenet Evol. 2012;62: 664–672. pmid:22138160

29. De Innocentiis S, Lesti A, Livi S, Rossi AR, Crosetti D, Sola L. Microsatellite markers reveal population structure in gilthead sea bream Sparus auratus from the Atlantic Ocean and Mediterranean Sea. Fish Sci. 2004;70: 852–859.

30. Zitari-Chatti R, Chatti N, Fulgione D, Caiazza I, Aprea G, Elouaer A, et al. Mitochondrial DNA variation in the caramote prawn Penaeus (Melicertus) kerathurus across a transition zone in the Mediterranean Sea. Genetica. 2009;136: 439–447. pmid:19109695

31. Serra IA, Innocenti A, Di Maida G, Calvo S, Migliaccio M, Zambianchi E, et al. Genetic structure in the Mediterranean seagrass Posidonia oceanica: disentangling past vicariance events from contemporary patterns of gene flow. Mol Ecol. 2010;19: 557–568. pmid:20051010

32. Berline L, Rammou A-M, Doglioli A, Molcard A, Petrenko A. A Connectivity-Based Eco-Regionalization Method of the Mediterranean Sea. MacKenzie BR, editor. PLoS One. Public Library of Science; 2014;9: e111978. pmid:25375212

33. Rossi V, Ser-Giacomi E, Lõpez C, Hernández-García E. Hydrodynamic provinces and oceanic connectivity from a transport network help designing marine reserves. Geophys Res Lett. 2014;41: 2883–2891.

34. Boero F. The future of the Mediterranean Sea Ecosystem: towards a different tomorrow. Rend Lincei. Springer Milan; 2015;26: 3–12.

35. Astraldi M, Balopoulos S, Candela J, Font J, Gacic M, Gasparini G., et al. The role of straits and channels in understanding the characteristics of Mediterranean circulation. Prog Oceanogr. 1999;44: 65–108.

36. Millot C, Taupier-Letage I. Circulation in the Mediterranean Sea. Hdb Environ Chem. 2005;5: 29–66.

37. Rius M, Turon X, Ordóñez V, Pascual M. Tracking invasion histories in the sea: Facing complex scenarios using multilocus data. PLoS One. 2012;7: e35815. pmid:22545140

38. Pineda M-C, Lorente B, López-Legentil S, Palacín C, Turon X. Stochasticity in space, persistence in time: genetic heterogeneity in harbour populations of the introduced ascidian Styela plicata. PeerJ. 2016;4: e2158. pmid:27366653

39. Maltagliati F, Lupi L, Castelli A, Pannacciulli FG. The genetic structure of the exotic ascidian Styela plicata (Tunicata) from Italian ports, with a re-appraisal of its worldwide genetic pattern. Mar Ecol. 2016;37: 492–502.

40. Wang IJ. Recognizing the temporal distinctions between landscape genetics and phylogeography. Mol Ecol. 2010;19: 2605–2608. pmid:20561197

41. Salat J. Review of hydrographic environmental factors that may influence anchovy habitats in northwestern Mediterranean. Sci Mar. 1996;60 (S2): 21–32.

42. Garcia-Lafuente J, Delgado J, Criado-Aldeanueva F, Bruno M, del Rio J, Vargas JM. Water mass circulation on the continental shelf of the Gulf of Cadiz. Deep Res Part Ii-Topical Stud Oceanogr. 2006;53: 1182–1197.

43. Tintore J, Laviolette PE, Blade I, Cruzado A. A Study of an Intense Density Front in the Eastern Alboran Sea: The Almeria-Oran Front. J Phys Oceanogr. 1988;18: 1384–1397.

44. Artegiani A, Paschini E, Russo A, Bregant D, Raicich F, Pinardi N, et al. The Adriatic Sea General Circulation. Part I: Air–Sea Interactions and Water Mass Structure. J Phys Oceanogr. 1997;27: 1492–1514.

45. Theocharis A, Balopoulos E, Kioroglou S, Kontoyiannis H, Iona A. A synthesis of the circulation and hydrography of the South Aegean Sea and the Straits of the Cretan Arc (March 1994 January 1995). Prog Oceanogr. 1999;44: 469–509.

46. Simion P, Philippe H, Baurain D, Jager M, Richter DJ, Di Franco A, et al. A Large and Consistent Phylogenomic Dataset Supports Sponges as the Sister Group to All Other Animals. Curr Biol. Ruprecht-Karls-University; 2017;432: 618–622.

47. Ballesteros E. Mediterranean coralligenous assemblages: A synthesis of present knowledge. Oceanogr Mar Biol. 2006;44: 123–195.

48. Uriz MJ, Maldonado M, Turon X, Martí R. How do reproductive output, larval behaviour, and recruitment contribute to adult spatial patterns in Mediterranean encrusting sponges? Mar Ecol Prog Ser. 1998;167: 137–148.

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

50. López S, Turon X, Montero E, Palacín C, Duarte C, Tarjuelo I. Larval abundance, recruitment and early mortality in Paracentrotus lividus (Echinoidea). Interannual variability and plankton-benthos coupling. Mar Ecol Prog Ser. 1998;172: 239–251.

51. Villanueva R, Vidal EAG, Fernández-Álvarez FA, Nabhitabhata J. Early Mode of Life and Hatchling Size in Cephalopod Molluscs: Influence on the Species Distributional Ranges. Kiel S, editor. PLoS One. Public Library of Science; 2016;11: e0165334. pmid:27829039

52. Cabanellas-Reboredo M, Deudero S, Alós J, Valencia JM, March D, Hendriks IE, et al. Recruitment of Pinna nobilis (Mollusca: Bivalvia) on artificial structures. Mar Biodivers Rec. Cambridge University Press; 2009;2: e126.

53. Ribeiro PMA. Dispersal and connectivity of Northeastern Atlantic Patellid limpets: a multidisciplinary approach. University of Southampton. 2008.

54. González-Gordillo JI, Dos Santos A, Rodríguez A. Checklist and annotated bibliography of decapod crustacean larvae from the Southwestern European coast (Gibraltar Strait area). Sci Mar. 2001;65: 275–305.

55. Tarjuelo I, Posada D, Crandall KA, Pascual M, Turon X. Phylogeography and speciation of colour morphs in the colonial ascidian Pseudodistoma crucigaster. Mol Ecol. 2004;13: 3125–3136. pmid:15367125

56. Raventos N, Macpherson E. Planktonic larval duration and settlement marks on the otoliths of Mediterranean littoral fishes. Mar Biol. 2001;138: 1115–1120.

57. Macpherson E, Raventos N. Relationship between pelagic larval duration and geographic distribution of Mediterranean littoral fishes. Mar Ecol Ser. 2006;327: 257–265.

58. Preacher KJ. Calculation for the chi-square test: An interactive calculation tool for chi-square tests of goodness of fit and independence [Internet]. 2001. Available: http://quantpsy.org

59. Stat Soft. Inc. STATISTICA (data analysis software system), version 8.0. www.statsoft.com. 2007.

60. Anderson MJ, Gorley RN, Clarke KR. PERMANOVA+ for PRIMER: Guide to Software and Statistical Methods. Primer-E. Plymouth, RU; 2008. p. 214.

61. Babbucci M, Buccoli S, Cau A, Cannas R, Goñi R, Díaz D, et al. Population structure, demographic history, and selective processes: Contrasting evidences from mitochondrial and nuclear markers in the European spiny lobster Palinurus elephas (Fabricius, 1787). Mol Phylogenet Evol. 2010;56: 1040–1050. pmid:20510378

62. Blanquer A, Uriz MJ. Population genetics at three spatial scales of a rare sponge living in fragmented habitats. Bmc Evol Biol. 2010;10.

63. Zulliger DE, Tanner S, Ruch M, Ribi G. Genetic structure of the high dispersal Atlanto-Mediterreanean sea star Astropecten aranciacus revealed by mitochondrial DNA sequences and microsatellite loci. Mar Biol. 2009;156: 597–610.

64. Riginos C, Douglas KE, Jin Y, Shanahan DF, Treml EA. Effects of geography and life history traits on genetic differentiation in benthic marine fishes. Ecography (Cop). Blackwell Publishing Ltd; 2011;34: 566–575.

65. Luiz OJ, Allen AP, Robertson DR, Floeter SR, Kulbicki M, Vigliola L, et al. Adult and larval traits as determinants of geographic range size among tropical reef fishes. Proc Natl Acad Sci U S A. National Academy of Sciences; 2013;110: 16498–502. pmid:24065830

66. Sabatés A. Distribution pattern of larval fish populations in the northwestern Mediterranean. Mar Ecol Prog Ser. 1990;59: 75–82.

67. 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: 3291–3305. pmid:21762434

68. Casado-Amezúa P, Goffredo S, Templado J, Machordom A. Genetic assessment of population structure and connectivity in the threatened mediterranean coral Astroides calycularis (Scleractinia, Dendrophylliidae) at different spatial scales. Mol Ecol. 2012;21: 3671–3685. pmid:22646530

69. Mokhtar-Jamaï K, Coma R, Wang J, Zuberer F, Féral J-P, 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: 1765–1779. pmid:23789084

70. 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: e0119585. pmid:25774522

71. Jahnke M, Christensen A, Micu D, Milchakova N, Sezgin M, Todorova V, et al. Patterns and mechanisms of dispersal in a keystone seagrass species. Mar Environ Res. 2016;117: 54–62. pmid:27085058

72. Sabatés A, Masó M. Effect of a shelf-slope front on the spatial distribution of mesopelagic fish larvae in the western Mediterranean. Deep Sea Res Part A, Oceanogr Res Pap. Elsevier; 1990;37: 1085–1098.

73. Sabatés A, Salat J, Olivar MP. Advection of continental water as an export mechanism for anchovy, Engraulis encrasicolus, larvae. Sci Mar. 2001;65: 77–87.

74. Sabatés A, Olivar MP, Salat J, Palomera I, Alemany F. Physical and biological processes controlling the distribution of fish larvae in the NW Mediterranean. Prog Oceanogr. 2007;74: 355–376.

75. Siokou-Frangou I, Christaki U, Mazzocchi MG, Montresor M, Ribera D’Alcala M, Vaque D, et al. Plankton in the open mediterranean Sea: A review. Biogeosciences. 2010;7: 1543–1586.

76. Shanks AL. Pelagic Larval Duration and Dispersal Distance Revisited. Biol Bull. 2009;216: 373–385. pmid:19556601

77. Harrison HB, Williamson DH, Evans RD, Almany GR, Thorrold SR, Russ GR, et al. Larval Export from Marine Reserves and the Recruitment Benefit for Fish and Fisheries. Curr Biol. 2012;22: 1023–1028. pmid:22633811

78. Carreras-Carbonell J, Macpherson E, Pascual M. Population structure within and between subspecies of the Mediterranean triplefin fish Tripterygion delaisi revealed by highly polymorphic microsatellite loci. Mol Ecol. 2006;15: 3527–3539. pmid:17032255

79. Schunter C, Pascual M, Garza JC, Raventos N, Macpherson E. Kinship analyses identify fish dispersal events on a temperate coastline. Proc R Soc B. 2014;281: 20140556. pmid:24812064

80. Balbín R, López-Jurado JL, Flexas MM, Reglero P, Vélez-Velchí P, González-Pola C, et al. Interannual variability of the early summer circulation around the Balearic Islands: Driving factors and potential effects on the marine ecosystem. J Mar Syst. 2014;138: 70–81.

81. Renault L, Oguz T, Pascual A, Vizoso G, Tintore J. Surface circulation in the Alborn Sea (western Mediterranean) inferred from remotely sensed data. J Geophys Res Ocean. 2012;117: 1–12.

82. Reuschel S, Cuesta JA, Schubart CD. Marine biogeographic boundaries and human introduction along the European coast revealed by phylogeography of the prawn Palaemon elegans. Mol Phylogenet Evol. 2010;55: 765–775. pmid:20307676

83. Calderón I, Pita L, Brusciotti S, Palacín C, Turon X. Time and space: genetic structure of the cohorts of the common sea urchin Paracentrotus lividus in Western Mediterranean. Mar Biol. 2012;159: 187–197.

84. Schunter C, Carreras-Carbonell J, Planes S, Sala E, Ballesteros E, Zabala M, et al. Genetic connectivity patterns in an endangered species: The dusky grouper (Epinephelus marginatus). J Exp Mar Bio Ecol. 2011;401: 126–133.

85. Hauser L, Carvalho GR. Paradigm shifts in marine fisheries genetics: ugly hypotheses slain by beautiful facts. Fish Fish. 2008;9: 333–362.

86. Sahyoun R, Guidetti P, Di Franco A, Planes S. Patterns of fish connectivity between a marine protected area and surrounding fished areas. PLoS One. Public Library of Science; 2016;11: e0167441. pmid:27907100

87. García-Rubies A, Cebrian E, Schembri PJ, Evans J, Macpherson E. Ecological Effects and Benefits of Mediterranean Marine Protected Areas: Management Implications. In: Goriup PD, editor. Management of Marine Protected Areas: A Network Perspective. John Wiley & Sons Ltd.; 2017.

88. Baguette M, Blanchet S, Legrand D, Stevens VM, Turlure C. Individual dispersal, landscape connectivity and ecological networks. Biol Rev. 2013;88: 310–326. pmid:23176626

89. Cerrano C, Danovaro R, Gambi C, Pusceddu A, Riva A, Schiaparelli S. Gold coral (Savalia savaglia) and gorgonian forests enhance benthic biodiversity and ecosystem functioning in the mesophotic zone. Biodivers Conserv. 2009;19: 153–167.

90. García-Rubies A, Macpherson E. Substrate use and temporal pattern of recruitment in juvenile fishes of the Mediterranean littoral. Mar Biol. 1995;124: 35–42.

91. Thibaut T, Bottin L, Aurelle D, Boudouresque C-F, Blanfuné A, Verlaque M, et al. Connectivity of Populations of the Seaweed Cystoseira amentacea within the Bay of Marseille (Mediterranean Sea): Genetic Structure and Hydrodynamic Connections. Cryptogam Algol. 2016;37: 233–255.

92. Coll M, Piroddi C, Albouy C, Ben Rais Lasram F, Cheung WWL, Christensen V, et al. The Mediterranean Sea under siege: spatial overlap between marine biodiversity, cumulative threats and marine reserves. Glob Ecol Biogeogr. 2012;21: 465–480.

93. Cebrian E, Martí R, Uriz JM, Turon X. Sublethal effects of contamination on the Mediterranean sponge Crambe crambe: metal accumulation and biological responses. Mar Pollut Bull. 2003;46: 1273–1284. pmid:14550340

94. Thibaut T, Pinedo S, Torras X, Ballesteros E. Long-term decline of the populations of Fucales (Cystoseira spp. and Sargassum spp.) in the Albères coast (France, North-western Mediterranean). Mar Pollut Bull. 2005;50: 1472–1489. pmid:16026805

95. Sales M, Ballesteros E, Vidal E, Tomas F, Moranta J, Cebrian E. New Method for restoring Cystoseira forest. 6th European Phycological Congress. London; 2015.

96. Bradbury IR, Hubert S, Higgins B, Borza T, Bowman S, Paterson IG, et al. Parallel adaptive evolution of Atlantic cod on both sides of the Atlantic Ocean in response to temperature. Proc R Soc London B Biol Sci. 2010;277: 3725–34.

97. Carreras C, Ordóñez V, Zane L, Kruschel C, Nasto I, Macpherson E, et al. Population genomics of an endemic Mediterranean fish: differentiation by fine scale dispersal and adaptation. Sci Rep. 2017;7.

98. Buonomo R, Assis J, Fernandes F, Engelen A, Airoldi L, Serrão E. Habitat continuity and stepping-stone oceanographic distances explain population genetic connectivity of the brown alga Cystoseira amentacea. Mol Ecol. 2016;26: 766–780.