Strong population genetic structure and contrasting demographic histories for the small-spotted catshark (Scyliorhinus canicula) in the Mediterranean Sea

Coastal and demersal chondrichthyans, such as the small-spotted catshark, are expected to exhibit genetic differentiation in areas of complex geomorphology like the Mediterranean Basin because of their limited dispersal ability. To test this hypothesis, we used a fragment of the mitochondrial cytochrome c oxidase subunit I gene and 12 nuclear microsatellite loci in order to investigate the genetic structure and historical demography of this species, and to identify potential barriers to gene flow. Samples were collected from the Balearic Islands, the Algerian Basin, the Ionian Sea, the Corinthian Gulf and various locations across the Aegean Sea. Additional sequences from the Atlantic and the Levantine Basin retrieved from GenBank were included in the mitochondrial DNA analysis. Both mitochondrial and nuclear microsatellite DNA data revealed a strong genetic subdivision, mainly between the western and eastern Mediterranean, whereas the Levantine Basin shared haplotypes with both areas. The geographic isolation of the Mediterranean basins seems to enforce the population genetic differentiation of the species, with the deep sea acting as a strong barrier to its dispersal. Contrasting historical demographic patterns were also observed in different parts of the species'' distribution, most notably a population growth trend in the western Mediterranean/Atlantic area and a slight decreasing one in the Aegean Sea. The different effects of the Pleistocene glacial periods on the habitat availability may explain the contrasting demographic patterns observed. The current findings suggest that the small-spotted catshark exhibits several genetic stocks in the Mediterranean, although further study is needed.     

Modelling the presence of anchovy Engraulis encrasicolus in the Aegean Sea during early summer, based on satellite environmental data

Acoustic and satellite environmental data as well as bathymetry data were used to model the presence of anchovy, Engraulis encrasicolus during early summer in the northern Aegean Sea (Eastern Mediterranean). Generalized Additive Models (GAMs) were used for modelling and subsequently applied in a predictive mode to identify those areas in the Greek Seas and the entire Mediterranean basin that could support species’ presence. Model results were evaluated with the estimation of Receiver Operating Characteristic (ROC)-plots as well as qualitatively, based on (a) acoustic data from concurrent studies in certain areas of the northern Aegean Sea that were not included in the estimation of the GAM model and (b) historical acoustic data from the central Aegean and Ionian Seas. Mapping the estimated environmental conditions in the Mediterranean basin indicated areas that generally agree with the known distribution grounds of anchovy, such as the straits of Sicily and coastal waters of Tunisia, areas in the Tyrrhenian Sea, the Adriatic Sea, the Gulf of Lions and the Catalan Sea. Guest editor: V. D. Valavanis Essential Fish Habitat Mapping in the Mediterranean     

Genetic and morphological analyses of tub gurnard Chelidonichthys lucerna populations in Turkish marine waters

Genetic and morphological structure of tub gurnard Chelidonichthys lucerna populations in Turkish marine waters were investigated with mtDNA sequencing of 16S rRNA and morphological characters. C. lucerna samples were collected from the Black Sea, Marmara, Aegean and northeastern Mediterranean coasts of Turkey. The lowest genetic diversity was found in the northeastern Mediterranean (Iskenderun Bay) population, while the highest was in the Marmara population with overall average value of genetic diversity within populations. A total of 14 haplotypes was found, and the highest haplotype diversity was in the Black Sea whereas the lowest was in the northeastern Mediterranean population (Iskenderun Bay). The Black Sea and Iskenderun Bay populations showed the least genetic divergence (0.001081), while the highest was between the Marmara Sea and northeastern Mediterranean (Antalya Bay) populations (0.002067). Pairwise comparisons of genetic distance revealed statistically significant differences (P < 0.05) between the Marmara and both the Aegean and northeastern Mediterranean (Antalya Bay) samples. Neighbour joining tree analyses clustered the northeastern Mediterranean populations (Antalya Bay and Iskenderun Bay) as genetically more interrelated populations, whereas the Aegean Sea population was clustered as most isolated one. Discriminant function analysis of morphological characters showed that only the Black Sea population is differentiated from the other populations.     

Population genomics meet Lagrangian simulations: Oceanographic patterns and long larval duration ensure connectivity among Paracentrotus lividus populations in the Adriatic and Ionian seas

Connectivity between populations influences both their dynamics and the genetic structuring of species. In this study, we explored connectivity patterns of a marine species with long‐distance dispersal, the edible common sea urchin Paracentrotus lividus, focusing mainly on the Adriatic–Ionian basins (Central Mediterranean). We applied a multidisciplinary approach integrating population genomics, based on 1,122 single nucleotide polymorphisms (SNPs) obtained from 2b‐RAD in 275 samples, with Lagrangian simulations performed with a biophysical model of larval dispersal. We detected genetic homogeneity among eight population samples collected in the focal Adriatic–Ionian area, whereas weak but significant differentiation was found with respect to two samples from the Western Mediterranean (France and Tunisia). This result was not affected by the few putative outlier loci identified in our dataset. Lagrangian simulations found a significant potential for larval exchange among the eight Adriatic–Ionian locations, supporting the hypothesis of connectivity of P. lividus populations in this area. A peculiar pattern emerged from the comparison of our results with those obtained from published P. lividus cytochrome b (cytb) sequences, the latter revealing genetic differentiation in the same geographic area despite a smaller sample size and a lower power to detect differences. The comparison with studies conducted using nuclear markers on other species with similar pelagic larval durations in the same Adriatic–Ionian locations indicates species‐specific differences in genetic connectivity patterns and warns against generalizing single‐species results to the entire community of rocky shore habitats.     

Microsatellite DNA reveals genetically different populations of Atlantic bonito Sarda sarda in the Mediterranean Basin

Microsatellites were used to investigate population genetic structure of Atlantic bonito Sarda sarda from the Black Sea, Marmara Sea, Aegean Sea, north-eastern Mediterranean Sea and Adriatic Sea. Overall average observed heterozygosity was high (0.93). Average observed heterozygosity per locus ranged from 0.79 to 0.98. Pairwise F_ST estimates for all loci between populations ranged from 0 to 0.07626, and significant F_ST values (P < 0.001) were detected between populations; the Blacks Sea and Marmara Sea samples were not significantly different from each other, but significant different from the other samples, and Aegean Sea and north-eastern Mediterranean Sea samples were also not significantly different from each other, but significantly different from all other samples. The Adriatic Sea sample was significant different from all other samples. The Mantel test revealed a significant (P < 0.001, r = 0.68) isolation-by-distance for these 11 populations. Neighbour-joining analysis clustered the Black Sea and Marmara Sea samples together while collections from Aegean Sea and north-eastern Mediterranean Sea were clustered close to each other and far from the others. On the other hand, the Adriatic Sea collection presented very distinctive relationship from the others.     

Ecoregion-Based Conservation Planning in the Mediterranean: Dealing with Large-Scale Heterogeneity

Spatial priorities for the conservation of three key Mediterranean habitats, i.e. seagrass Posidonia oceanica meadows, coralligenous formations, and marine caves, were determined through a systematic planning approach. Available information on the distribution of these habitats across the entire Mediterranean Sea was compiled to produce basin-scale distribution maps. Conservation targets for each habitat type were set according to European Union guidelines. Surrogates were used to estimate the spatial variation of opportunity cost for commercial, non-commercial fishing, and aquaculture. Marxan conservation planning software was used to evaluate the comparative utility of two planning scenarios: (a) a whole-basin scenario, referring to selection of priority areas across the whole Mediterranean Sea, and (b) an ecoregional scenario, in which priority areas were selected within eight predefined ecoregions. Although both scenarios required approximately the same total area to be protected in order to achieve conservation targets, the opportunity cost differed between them. The whole-basin scenario yielded a lower opportunity cost, but the Alboran Sea ecoregion was not represented and priority areas were predominantly located in the Ionian, Aegean, and Adriatic Seas. In comparison, the ecoregional scenario resulted in a higher representation of ecoregions and a more even distribution of priority areas, albeit with a higher opportunity cost. We suggest that planning at the ecoregional level ensures better representativeness of the selected conservation features and adequate protection of species, functional, and genetic diversity across the basin. While there are several initiatives that identify priority areas in the Mediterranean Sea, our approach is novel as it combines three issues: (a) it is based on the distribution of habitats and not species, which was rarely the case in previous efforts, (b) it considers spatial variability of cost throughout this socioeconomically heterogeneous basin, and (c) it adopts ecoregions as the most appropriate level for large-scale planning.     

Mitochondrial DNA Sequence Variation Suggests the Lack of Genetic Heterogeneity in the Adriatic and Ionian Stocks of Sardina pilchardus

A genetic stock structure analysis of 11 sardine samples from the Adriatic Sea and Ionian neighboring area was carried out through sequence variation analysis of a 307-bp cytochrome b gene fragment in order to identify self-recruiting units in the Adriatic Sardina pilchardus stock. The overall lack of genetic subdivision among samples detected by analysis of molecular variance, pairwise Φ_st values, and the exact test of population differentiation indicates this sardine stock is part of a larger self-recruiting population whose boundaries are larger than the investigated area. This conclusion is in agreement with preliminary allozymic and mitochondrial DNA restriction fragment length polymorphism data, but contradicts the previous identification of 2 subpopulations of sardines in the Adriatic Sea argued on morphologic differences, which could be rather attributed to different hydrographic or ecologic conditions occurring in different areas of the Adriatic Sea. The reduced gene flow observed between Adriatic-Ionian and Spanish sardine geographic samples (P < 0.001) suggests that reproductively isolated populations of sardines may occur in the Mediterranean Sea.     

The Yellow Gorgonian Eunicella cavolini: Demography and Disturbance Levels across the Mediterranean Sea

The yellow octocoral Eunicella cavolini is one of the most common gorgonians thriving in Mediterranean hard-bottom communities. However, information regarding its distribution and ecology in several parts of the Mediterranean is lacking, while population trends and conservation status remain largely unknown. We investigated 19 populations of E. cavolini over three representative geographic regions: the NW Mediterranean, CE Adriatic, and N Aegean. Focusing on the upper bathymetric range of the species (<40 m), data were collected on the populations’ upper depth limit, density, colony height, and extent of injury. A three-level hierarchical sampling design was applied to assess the existence of spatial patterns, using: a) regions (located thousands of km apart), b) localities within regions (tens to hundreds of km apart), and c) sites within localities (hundreds of m to a few km apart). In the NW Mediterranean and CE Adriatic, the upper distribution limit was at depths ≤15 m, whereas in the N Aegean most populations were found deeper than 30 m. Population density ranged between 4.46-62 colonies per m², while mean colony height was 15.6±8.9 SD cm with a maximum of 62 cm. The NW Mediterranean sites were characterized by dense populations dominated by small colonies (<20 cm), periodic recruitment, and low proportion of large gorgonians (>30 cm). The CE Adriatic displayed intermediate densities, with well-structured populations, and continuous recruitment. In the N Aegean, most populations presented low densities, high proportion of large colonies, but low number of small colonies, signifying limited recruitment. Disturbance levels, as a function of extent and type of injury, are discussed in relation to past or present human-induced threats. This work represents geographically the most wide ranging demographic study of a Mediterranean octocoral to date. The quantitative information obtained provides a basis for future monitoring at a Mediterranean scale.     

Allozyme differentiation of bonito in the Mediterranean Sea

Bonito Sarda sarda samples collected in three areas from the Mediterranean Sea exhibited considerable genetic heterogeneity over all locations. This provided preliminary evidence for two different groups of bonito, one for the Ligurian and the Ionian Seas and the other for the Aegean Sea.     

Genetic and historic evidence for climate-driven population fragmentation in a top cetacean predator: the harbour porpoises in European water

Recent climate change has triggered profound reorganization in northeast Atlantic ecosystems, with substantial impact on the distribution of marine assemblages from plankton to fishes. However, assessing the repercussions on apex marine predators remains a challenging issue, especially for pelagic species. In this study, we use Bayesian coalescent modelling of microsatellite variation to track the population demographic history of one of the smallest temperate cetaceans, the harbour porpoise (Phocoena phocoena) in European waters. Combining genetic inferences with palaeo-oceanographic and historical records provides strong evidence that populations of harbour porpoises have responded markedly to the recent climate-driven reorganization in the eastern North Atlantic food web. This response includes the isolation of porpoises in Iberian waters from those further north only approximately 300 years ago with a predominant northward migration, contemporaneous with the warming trend underway since the ‘Little Ice Age’ period and with the ongoing retreat of cold-water fishes from the Bay of Biscay. The extinction or exodus of harbour porpoises from the Mediterranean Sea (leaving an isolated relict population in the Black Sea) has lacked a coherent explanation. The present results suggest that the fragmentation of harbour distribution range in the Mediterranean Sea was triggered during the warm ‘Mid-Holocene Optimum’ period (approx. 5000 years ago), by the end of the post-glacial nutrient-rich ‘Sapropel’ conditions that prevailed before that time.     

The genetic population structure of Thunnus thynnus (Linnaeus, 1758) in the Mediterranean Sea,a controversial issue

The Atlantic bluefin tuna (ABFT), Thunnus thynnus (Linnaeus, 1758), is an important commercial species managed as two different stocks, western and eastern Atlantic, with their spawning grounds in the Gulf of Mexico and in the Mediterranean Sea, respectively. The eastern Atlantic stock has been overexploited in the last decades, leading to the application of specific management measures introduced by the International Commission for the Atlantic Tuna (ICCAT). A clear understanding of the genetic structure of ABFT Mediterranean population should be pursued in order to support management decisions. To date the genetic studies on the Mediterranean ABFT, carried out with different molecular markers and sampling procedures, have produced unclear results. Here, we analysed ABFT samples from central and western Mediterranean Sea with mitochondrial sequences and 11 microsatellite loci to investigate, among the others, the area of the Strait of Messina, where environmental conditions seem to support a resident population of ABFT. Furthermore, genetic analyses of mitochondrial sequences were carried out including nucleotide sequences of Adriatic ABFT wild larvae retrieved from GenBank. Among the investigated areas a genetic differentiation was detected between the Strait of Messina and the Tyrrhenian Sea with microsatellite loci according to the exact G test, but not to the Bayesian analyses carried out with STRUCTURE. The analyses with mitochondrial sequences do not reveal any differentiation among sampled areas, however, a highly significant genetic divergence was observed between the Adriatic mitochondrial sequences retrieved from GenBank and the central‐western Mediterranean sequences obtained in the present work. Our results provide some evidence of population structure of Mediterranean ABFT adding pieces to a still unclear picture.     

Into the depth of population genetics: pattern of structuring in mesophotic red coral populations

Deep-sea reef-building corals are among the most conspicuous invertebrates inhabiting the hard-bottom habitats worldwide and are particularly susceptible to human threats. The precious red coral (Corallium rubrum, L. 1758) has a wide bathymetric distribution, from shallow up to 800 m depth, and represents a key species in the Mediterranean mesophotic reefs. Several studies have investigated genetic variability in shallow-water red coral populations, while geographic patterns in mesophotic habitats are largely unknown. This study investigated genetic variability of C. rubrum populations dwelling between 55 and 120 m depth, from the Ligurian to the Ionian Sea along about 1500 km of coastline. A total of 18 deep rocky banks were sampled. Colonies were analyzed by means of a set of microsatellite loci and the putative control region of the mitochondrial DNA. Collected data were compared with previous studies. Both types of molecular markers showed high genetic similarity between populations within the northern (Ligurian Sea and Tuscan Archipelago) and the southern (Tyrrhenian and Ionian seas) study areas. Variability in habitat features between the sampling sites did not affect the genetic variability of the populations. Conversely, the patchy distribution of suitable habitats affected populations’ connectivity within and among deep coral banks. Based on these results and due to the emphasis on red coral protection in the Mediterranean Sea by international institutions, red coral could be promoted as a ‘focal species’ to develop management plans for the conservation of deep coralligenous reefs, a reservoir of marine biodiversity.     

The Mediterranean Sea under siege: spatial overlap between marine biodiversity,cumulative threats and marine reserves

Aim A large body of knowledge exists on individual anthropogenic threats that have an impact on marine biodiversity in the Mediterranean Sea, although we know little about how these threats accumulate and interact to affect marine species and ecosystems. In this context, we aimed to identify the main areas where the interaction between marine biodiversity and threats is more pronounced and to assess their spatial overlap with current marine protected areas in the Mediterranean. Location Mediterranean Sea. Methods We first identified areas of high biodiversity of marine mammals, marine turtles, seabirds, fishes and commercial or well‐documented invertebrates. We mapped potential areas of high threat where multiple threats are occurring simultaneously. Finally we quantified the areas of conservation concern for biodiversity by looking at the spatial overlap between high biodiversity and high cumulative threats, and we assessed the overlap with protected areas. Results Our results show that areas with high marine biodiversity in the Mediterranean Sea are mainly located along the central and north shores, with lower values in the south‐eastern regions. Areas of potential high cumulative threats are widespread in both the western and eastern basins, with fewer areas located in the south‐eastern region. The interaction between areas of high biodiversity and threats for invertebrates, fishes and large animals in general (including large fishes, marine mammals, marine turtles and seabirds) is concentrated in the coastal areas of Spain, Gulf of Lions, north‐eastern Ligurian Sea, Adriatic Sea, Aegean Sea, south‐eastern Turkey and regions surrounding the Nile Delta and north‐west African coasts. Areas of concern are larger for marine mammal and seabird species. Main conclusions These areas may represent good candidates for further research, management and protection activities, since there is only a maximum 2% overlap between existing marine protected areas (which cover 5% of the Mediterranean Sea) and our predicted areas of conservation concern for biodiversity.     

Assessing the geographic scale of genetic population management with microsatellites and introns in the clam Ruditapes decussatus

The clam Ruditapes decussatus is commercially important in southwestern Europe, suffering from population decline and hybridization with exotic Manila clam (R. philippinarum). Previous studies with intronic markers showed a genetic subdivision of the species in three races (Atlantic, West Mediterranean, and Adriatic‐Aegean). However, detailed population genetic studies to help management of the main production areas in the southwest of Europe are missing. We have analyzed eight Atlantic and two Mediterranean populations from the Spanish coasts using 14 microsatellites and six intronic markers. Microsatellites confirmed the Atlantic and West Mediterranean races detected with introns and showed that genetic variability was higher in Mediterranean than in Atlantic populations. Both marker types showed that genetic differentiation of Atlantic populations was low and indicated that populations could be managed at the regional level in the case of Cantabrian and Gulf of Cadiz areas, but not in the case of Rias Baixas and the Mediterranean. This study shows the interest of including different types of markers in studies of genetic population structure of marine organisms.     

Demosponge distribution in the eastern Mediterranean: a NW–SE gradient

The purpose of this paper was to investigate patterns of demosponge distribution along gradients of environmental conditions in the biogeographical subzones of the eastern Mediterranean (Aegean and Levantine Sea). The Aegean Sea was divided into six major areas on the basis of its geomorphology and bathymetry. Two areas of the Levantine Sea were additionally considered. All available data on demosponge species numbers and abundance in each area, as well as their vertical and general geographical distribution were taken from the literature. Multivariate analysis revealed a NW–SE faunal gradient, showing an apparent dissimilarity among the North Aegean, the South Aegean and the Levantine Sea, which agrees with the differences in the geographical, physicochemical and biological characteristics of the three areas. The majority of demosponge species has been recorded in the North Aegean, while the South Aegean is closer, in terms of demosponge diversity, to the oligotrophic Levantine Sea. The number of studies in the Aegean and Levantine subareas was positively correlated with the number of species recorded within each Aegean subarea. Demosponge species with an Altanto-Mediterranean distribution prevailed in the Aegean and the Levantine. The reduced contribution of the endemic component, as compared to the western Mediterranean, is consistent with the general NW–SE decrease in the number of endemic species in the Mediterranean. Demosponge distribution at the order level showed also a NW–SE gradient, similar to that observed in the entire Atlantic–Mediterranean region, suggesting a warm water affinity. Sublittoral, circalittoral and bathyal zones were clearly distinguishable in the Aegean Sea on the basis of their sponge fauna. The total number of species was an exponential negative function of depth.     

Genetic architecture of the marbled goby Pomatoschistus marmoratus (Perciformes, Gobiidae) in the Mediterranean Sea

The marbled goby Pomatoschistus marmoratus, a species inhabiting coastal Mediterranean lagoons, has been studied by measuring its mitochondrial DNA variation. This analysis revealed a Mediterranean west vs east split and, subsequently, an eastern differentiation among the Libyan-Tunisian Gulf, the Adriatic Sea and the Aegean Sea. The high cohesion between the samples collected in the vast area of western Mediterranean contrasts with the genetic mosaic of the more sub-structured eastern Mediterranean. This western homogeneity can not yet be fully explained even if a human-mediated migratory flow, due to a maritime traffic, has been posited. The pattern in the eastern basin revealed a genetic architecture possibly due to the non-migratory habit of the gobid. Within this perspective, the role of the Mediterranean lagoon habitat should be related to how much it amplifies the effects of historical (e.g. past sea-level changes) and environmental (e.g. present-day hydrographic regime) processes as regards the genetic structure of the inhabiting species.     

Assessing Dispersal Patterns of Fish Propagules from an Effective Mediterranean Marine Protected Area

Successfully enforced marine protected areas (MPAs) have been widely demonstrated to allow, within their boundaries, the recovery of exploited species and beyond their boundaries, the spillover of juvenile and adult fish. Little evidence is available about the so-called ‘recruitment subsidy’, the augmented production of propagules (i.e. eggs and larvae) due to the increased abundance of large-sized spawners hosted within effective MPAs. Once emitted, propagules can be locally retained and/or exported elsewhere. Patterns of propagule retention and/or export from MPAs have been little investigated, especially in the Mediterranean. This study investigated the potential for propagule production and retention/export from a Mediterranean MPA (Torre Guaceto, SW Adriatic Sea) using the white sea bream, Diplodus sargus sargus, as a model species. A multidisciplinary approach was used combining 1) spatial distribution patterns of individuals (post-settlers and adults) assessed through visual census within Torre Guaceto MPA and in northern and southern unprotected areas, 2) Lagrangian simulations of dispersal based on an oceanographic model of the region and data on early life-history traits of the species (spawning date, pelagic larval duration) and 3) a preliminary genetic study using microsatellite loci. Results show that the MPA hosts higher densities of larger-sized spawners than outside areas, potentially guaranteeing higher propagule production. Model simulations and field observation suggest that larval retention within and long-distance dispersal across MPA boundaries allow the replenishment of the MPA and of exploited populations up to 100 km down-current (southward) from the MPA. This pattern partially agrees with the high genetic homogeneity found in the entire study area (no differences in genetic composition and diversity indices), suggesting a high gene flow. By contributing to a better understanding of propagule dispersal patterns, these findings provide crucial information for the design of MPAs and MPA networks effective to replenish fish stocks and enhance fisheries in unprotected areas.     

Notes on the distribution and biology of the deep-sea crab Bathynectes maravigna (Brachyura: Portunidae) in the Mediterranean Sea

Data on the distribution and biology of the deep-sea portunid crab Bathynectes maravigna are reported for the Mediterranean Sea, based on several fisheries research surveys. Densities are low and, therefore, biological data are scarce. In the western Mediterranean, the species is much commoner in Alborán Sea than in the Catalano-Balearic Sea. Occurrences are also scarce in the southern Adriatic and northwestern Thyrrenian Sea, as well as in the Ionian Sea. The Alborán Sea and the seas surrounding the southern Italian peninsula are the areas where densities are the highest. The occurrence depth range was found to be 245–786 m, but most of the occurrences took place deeper than 500 m. Sizes ranged between 9 and 51 mm carapace length in males and between 12 and 51 mm in females. Ovigerous females have been only reported in October–December and March–May. Eighty three percent of both males and females are right-handed. Sexual dimorphism was present in cheliped length with males having longer chelae than females. The species appears to be much commoner in those areas where Atlantic influence is the highest.