To Be, or Not to Be: That Is the Hamletic Question of Cryptic Evolution in the Eastern Atlantic and Mediterranean Raja miraletus Species Complex.

Despite a high species diversity, skates (Rajiformes) exhibit remarkably conservative morphology and ecology. Limited trait variations occur within and between species, and cryptic species have been reported among sister and non-sister taxa, suggesting that species complexes may be subject to stabilising selection. Three sibling species are currently recognised in the Raja miraletus complex: (i) R. miraletus occurring along the Portuguese and Mediterranean coasts, (ii) R. parva in the Central-Eastern Atlantic off West Africa and (iii) R. ocellifera in the Western Indian Ocean off South Africa. In the present study, the genetic variation at mitochondrial and nuclear markers was estimated in the species complex by analysing 323 individuals sampled across most of its geographical distribution area to test the hypothesis that restricted gene flow and genetic divergence within species reflect known climate and bio-oceanographic discontinuities. Our results support previous morphological studies and confirm the known taxonomic boundaries of the three recognised species. In addition, we identified multiple weakly differentiated clades in the Northeastern Atlantic Ocean and Mediterranean, at least two additional cryptic taxa off Senegal and Angola, a pronounced differentiation of ancient South African clades. The hidden genetic structure presented here may represent a valuable support to species' conservation action plans.

Oxygen and carbon isotope variations in Chamelea gallina shells: Environmental influences and vital effects.

Stable isotopes in mollusc shells, together with variable growth rates and other geochemical properties, can register different environmental clues, including seawater temperature, salinity and primary productivity. However, the strict biological control over the construction of biominerals exerted by many calcifying organisms can constrain the use of these organisms for paleoenvironmental reconstructions. Biologically controlled calcification is responsible for the so called vital effects that cause a departure from isotopic equilibrium during shell formation, resulting in lower shell oxygen and carbon compared to the equilibrium value. We investigated shell oxygen and carbon isotopic composition of the bivalve Chamelea gallina in six sites along with a latitudinal gradient on the Adriatic Sea (NE Mediterranean Sea). Seawater δ18 O and δ13 CDIC varied from North to South, reflecting variations in seawater temperature, salinity, and chlorophyll concentration among sites. Shell δ18 O and δ13 C differed among sites and exhibited a wide range of values along with the ~400 km latitudinal gradient, away from isotopic equilibrium for both isotopes. These results hampered the utilization of this bivalve as a proxy for environmental reconstructions, in spite of C. gallina showing promise as a warm temperature proxy. Rigorous calibration studies with a precise insight of environment and shell growth are crucial prior to considering this bivalve as a reliable paleoclimatic archive.

Multiscale analysis on otolith structural features reveals differences in ontogenesis and sex in Merluccius merluccius in the western Adriatic Sea.

Otolith biomineralization results from biochemical processes regulated by the interaction of internal (physiological) and external (environmental) factors which lead to morphological and ultrastructural variability at intra- and interspecific levels. The aim of this study was to conduct a multi-scale analysis of the sagittal otoliths of the Merlucius merlucius (European hake) from the western Adriatic Sea in order to correlate otolith features with fish ontogeny and sex. We show that otoliths of sexually undifferentiated (non-sexed) individuals having a fish body total length (TL) less than 15 cm had faster growth in length, width, area, perimeter, volume and weight and a higher amount of organic matrix compared with otoliths of sexually differentiated individuals (females and males) having a fish size range of 15-50 cm. Most importantly, with increasing fish TL, female saccular otoliths contained a higher number of protuberances and rougher surface compared with male specimens, which showed more uniform mean curvature density. The differences between females and males discovered in this study could be associated with fish hearing adaptation to reproductive behavioural strategies during the spawning season. The outcomes of this research provide insights on how size and sex-related variations in otolith features may be affected by fish ecological and behavioural patterns.

DNA metabarcoding suggests dietary niche partitioning in the Adriatic European hake.

The Northern Adriatic Sea (FAO Geographical Sub-Area 17) is one of the most productive fishing areas of the Mediterranean Sea and it includes a broad diversity of habitats. In the Northern Adriatic basin, the Pomo Pit (200-273 m of depth) is one of the most important areas of aggregation for some demersal stocks shared in the Adriatic Sea and it is an important spawning/nursery area of the European hake (Merluccius merluccius). Through a metabarcoding approach we investigated the feeding habits of European hake, both inside and outside the Pomo Pit, and their temporal variability comparing samples collected in 2016 and 2014. Our analyses proved the presence of an ontogenetic shift from a diet based mainly on crustaceans in juveniles to a more piscivorous feeding behaviour in adult hakes and suggested the presence of a specific niche partitioning and food preferences between hakes living inside and outside the Pomo Pit. The main differences among adult hakes refer to the presence of molluscs in the stomachs of hakes collected within the Pomo Pit and the presence of high depth prey species (i.e., Micromesistius poutassou). Metabarcoding revealed the relevant ecological role played by the Pomo Pit in M. merluccius feeding behaviour and ontogenetic development, promoting a careful ecosystem-based management of fisheries in this area through focused conservation measures.

Climate variation during the Holocene influenced the skeletal properties of Chamelea gallina shells in the North Adriatic Sea (Italy).

Understanding how marine taxa will respond to near-future climate changes is one of the main challenges for management of coastal ecosystem services. Ecological studies that investigate relationships between the environment and shell properties of commercially important marine species are commonly restricted to latitudinal gradients or small-scale laboratory experiments. This paper aimed to explore the variations in shell features and growth of the edible bivalve Chamelea gallina from the Holocene sedimentary succession to present-day thanatocoenosis of the Po Plain-Adriatic Sea system (Italy). Comparing the Holocene sub-fossil record to modern thanatocoenoses allowed obtaining an insight of shell variations dynamics on a millennial temporal scale. Five shoreface-related assemblages rich in C. gallina were considered: two from the Middle Holocene, when regional sea surface temperatures were higher than today, representing a possible analogue for the near-future global warming, one from the Late Holocene and two from the present-day. We investigated shell biometry and skeletal properties in relation to the valve length of C. gallina. Juveniles were found to be more porous than adults in all horizons. This suggested that C. gallina promoted an accelerated shell accretion with a higher porosity and lower density at the expense of mechanically fragile shells. A positive correlation between sea surface temperature and both micro-density and bulk density were found, with modern specimens being less dense, likely due to lower aragonite saturation state at lower temperature, which could ultimately increase the energetic costs of shell formation. Since no variation was observed in shell CaCO3 polymorphism (100% aragonite) or in compositional parameters among the analyzed horizons, the observed dynamics in skeletal parameters are likely not driven by a diagenetic recrystallization of the shell mineral phase. This study contributes to understand the response of C. gallina to climate-driven environmental shifts and offers insights for assessing anthropogenic impacts on this economic relevant species.

Environmental influence on calcification of the bivalve Chamelea gallina along a latitudinal gradient in the Adriatic Sea.

Environmental factors are encoded in shells of marine bivalves in the form of geochemical properties, shell microstructure and shell growth rate. Few studies have investigated how shell growth is affected by habitat conditions in natural populations of the commercial clam Chamelea gallina. Here, skeletal parameters (micro-density and apparent porosity) and growth parameters (bulk density, linear extension and net calcification rates) were investigated in relation to shell sizes and environmental parameters along a latitudinal gradient in the Adriatic Sea (400 km). Net calcification rates increased with increasing solar radiation, sea surface temperature and salinity and decreasing Chlorophyll concentration in immature and mature shells. In immature shells, which are generally more porous than mature shells, enhanced calcification was due to an increase in bulk density, while in mature shells was due to an increase in linear extension rates. The presence of the Po river in the Northern Adriatic Sea was likely the main driver of the fluctuations observed in environmental parameters, especially salinity and Chlorophyll concentration, and seemed to negatively affect the growth of C. gallina.

A metabarcoding approach for the feeding habits of European hake in the Adriatic Sea.

European hake (Merluccius merluccius) is one of the most economically important fish for the Mediterranean Sea. It is an important predator of deep upper shelf slope communities currently characterized by growth overexploitation: the understanding of hake's diet might support next generation management tools. However, all current European hake diet studies depend on the morphological identification of prey remains in stomach content, with consequent limitations. In this study, we set up a metabarcoding approach based on cytochrome oxidase I PCR amplification and Miseq Illumina paired-end sequencing of M. merluccius stomach content remains and compared the results to classic morphological analyses. A total of 95 stomach contents of M. merluccius sampled in the North-Central Adriatic Sea were analyzed with both the metabarcoding and morphological approaches. Metabarcoding clearly outperformed the morphological method in the taxonomic identification of prey describing more complex trophic relationships even when considering the morphological identification of 200 stomach contents. Statistical analysis of diet composition revealed a weak differentiation among the hake's size classes, confirming an opportunistic feeding behavior. All the analyses performed showed the presence of a core of shared prey among the size classes and a cloud of size-specific prey. Our study highlights the exceptional potential of metabarcoding as an approach to provide unprecedented taxonomic resolution in the diet of M. merluccius and potentially of other marine predators, due to the broad-spectrum of detection of the primers used. A thorough description of these complex trophic relationships is fundamental for the implementation of an ecosystem approach to fisheries.

Natural history and molecular evolution of demersal Mediterranean sharks and skates inferred by comparative phylogeographic and demographic analyses.

BACKGROUND: The unique and complex paleoclimatic and paleogeographic events which affected the Mediterranean Sea since late Miocene deeply influenced the distribution and evolution of marine organisms and shaped their genetic structure. Following the Messinian salinity crisis and the sea-level fluctuations during the Pleistocene, several Mediterranean marine species developed deep genetic differentiation, and some underwent rapid radiation. Here, we consider two of the most prioritized groups for conservation in the light of their evolutionary history: sharks and rays (elasmobranchs). This paper deals with a comparative multispecies analysis of phylogeographic structure and historical demography in two pairs of sympatric, phylogenetically- and ecologically-related elasmobranchs, two scyliorhinid catsharks (Galeus melastomus, Scyliorhinus canicula) and two rajid skates (Raja clavata, Raja miraletus). Sampling and experimental analyses were designed to primarily test if the Sicilian Channel can be considered as effective eco-physiological barrier for Mediterranean demersal sympatric elasmobranchs. METHODS: The phylogeography and the historical demography of target species were inferred by analysing the nucleotide variation of three mitochondrial DNA markers (i.e., partial sequence of COI, NADH2 and CR) obtained from a total of 248 individuals sampled in the Western and Eastern Mediterranean Sea as well as in the adjacent northeastern Atlantic Ocean. Phylogeographic analysis was performed by haplotype networking and testing spatial genetic differentiation of samples (i.e., analysis of molecular variance and of principal components). Demographic history of Mediterranean populations was reconstructed using mismatch distribution and Bayesian Skyline Plot analyses. RESULTS: No spatial genetic differentiation was identified in either catshark species, while phylogeographic structure of lineages was identified in both skates, with R. miraletus more structured than R. clavata. However, such structuring of skate lineages was not consistent with the separation between Western and Eastern Mediterranean. Sudden demographic expansions occurred synchronously during the upper Pleistocene (40,000-60,000 years ago) in both skates and G. melastomus, likely related to optimal environmental conditions. In contrast, S. canicula experienced a slow and constant increase in population size over the last 350,000 years. DISCUSSION: The comparative analysis of phylogeographic and historical demographic patterns for the Mediterranean populations of these elasmobranchs reveals that historical phylogeographic breaks have not had a large impact on their microevolution. We hypothesize that interactions between environmental and ecological/physiological traits may have been the driving force in the microevolution of these demersal elasmobranch species in the Mediterranean rather than oceanographic barriers.

Genetic differentiation and phylogeography of Mediterranean-North Eastern Atlantic blue shark (Prionace glauca, L. 1758) using mitochondrial DNA: panmixia or complex stock structure?

BACKGROUND: The blue shark (Prionace glauca, Linnaeus 1758) is one of the most abundant epipelagic shark inhabiting all the oceans except the poles, including the Mediterranean Sea, but its genetic structure has not been confirmed at basin and interoceanic distances. Past tagging programs in the Atlantic Ocean failed to find evidence of migration of blue sharks between the Mediterranean and the adjacent Atlantic, despite the extreme vagility of the species. Although the high rate of by-catch in the Mediterranean basin, to date no genetic study on Mediterranean blue shark was carried out, which constitutes a significant knowledge gap, considering that this population is classified as "Critically Endangered", unlike its open-ocean counterpart. METHODS: Blue shark phylogeography and demography in the Mediterranean Sea and North-Eastern Atlantic Ocean were inferred using two mitochondrial genes (Cytb and control region) amplified from 207 and 170 individuals respectively, collected from six localities across the Mediterranean and two from the North-Eastern Atlantic. RESULTS: Although no obvious pattern of geographical differentiation was apparent from the haplotype network, Φst analyses indicated significant genetic structure among four geographical groups. Demographic analyses suggest that these populations have experienced a constant population expansion in the last 0.4-0.1 million of years. DISCUSSION: The weak, but significant, differences in Mediterranean and adjacent North-eastern Atlantic blue sharks revealed a complex phylogeographic structure, which appears to reject the assumption of panmixia across the study area, but also supports a certain degree of population connectivity across the Strait of Gibraltar, despite the lack of evidence of migratory movements observed by tagging data. Analyses of spatial genetic structure in relation to sex-ratio and size could indicate some level of sex/stage biased migratory behaviour.

Improving the Conservation of Mediterranean Chondrichthyans: The ELASMOMED DNA Barcode Reference Library.

Cartilaginous fish are particularly vulnerable to anthropogenic stressors and environmental change because of their K-selected reproductive strategy. Accurate data from scientific surveys and landings are essential to assess conservation status and to develop robust protection and management plans. Currently available data are often incomplete or incorrect as a result of inaccurate species identifications, due to a high level of morphological stasis, especially among closely related taxa. Moreover, several diagnostic characters clearly visible in adult specimens are less evident in juveniles. Here we present results generated by the ELASMOMED Consortium, a regional network aiming to sample and DNA-barcode the Mediterranean Chondrichthyans with the ultimate goal to provide a comprehensive DNA barcode reference library. This library will support and improve the molecular taxonomy of this group and the effectiveness of management and conservation measures. We successfully barcoded 882 individuals belonging to 42 species (17 sharks, 24 batoids and one chimaera), including four endemic and several threatened ones. Morphological misidentifications were found across most orders, further confirming the need for a comprehensive DNA barcoding library as a valuable tool for the reliable identification of specimens in support of taxonomist who are reviewing current identification keys. Despite low intraspecific variation among their barcode sequences and reduced samples size, five species showed preliminary evidence of phylogeographic structure. Overall, the ELASMOMED initiative further emphasizes the key role accurate DNA barcoding libraries play in establishing reliable diagnostic species specific features in otherwise taxonomically problematic groups for biodiversity management and conservation actions.

Shell properties of commercial clam Chamelea gallina are influenced by temperature and solar radiation along a wide latitudinal gradient.

Phenotype can express different morphologies in response to biotic or abiotic environmental influences. Mollusks are particularly sensitive to different environmental parameters, showing macroscale shell morphology variations in response to environmental parameters. Few studies concern shell variations at the different scale levels along environmental gradients. Here, we investigate shell features at the macro, micro and nanoscale, in populations of the commercially important clam Chamelea gallina along a latitudinal gradient (~400 km) of temperature and solar radiation in the Adriatic Sea (Italian cost). Six populations of clams with shells of the same length were analyzed. Shells from the warmest and the most irradiated population were thinner, with more oval shape, more porous and lighter, showing lower load fracture. However, no variation was observed in shell CaCO3 polymorphism (100% aragonite) or in compositional and textural shell parameters, indicating no effect of the environmental parameters on the basic processes of biomineralization. Because of the importance of this species as commercial resource in the Adriatic Sea, the experimentally quantified and significant variations of mass and fracture load in C. gallina shells along the latitudinal gradient may have economic implications for fisheries producing different economical yield for fishermen and consumers along the Adriatic coastline.

Outlier SNP markers reveal fine-scale genetic structuring across European hake populations (Merluccius merluccius).

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome-wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large- and fine-scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (F(CT) = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (F(CT) range 0.275-0.705) and fine-scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.

Genetic structure of bluefin tuna in the mediterranean sea correlates with environmental variables.

BACKGROUND: Atlantic Bluefin Tuna (ABFT) shows complex demography and ecological variation in the Mediterranean Sea. Genetic surveys have detected significant, although weak, signals of population structuring; catch series analyses and tagging programs identified complex ABFT spatial dynamics and migration patterns. Here, we tested the hypothesis that the genetic structure of the ABFT in the Mediterranean is correlated with mean surface temperature and salinity. METHODOLOGY: We used six samples collected from Western and Central Mediterranean integrated with a new sample collected from the recently identified easternmost reproductive area of Levantine Sea. To assess population structure in the Mediterranean we used a multidisciplinary framework combining classical population genetics, spatial and Bayesian clustering methods and a multivariate approach based on factor analysis. CONCLUSIONS: FST analysis and Bayesian clustering methods detected several subpopulations in the Mediterranean, a result also supported by multivariate analyses. In addition, we identified significant correlations of genetic diversity with mean salinity and surface temperature values revealing that ABFT is genetically structured along two environmental gradients. These results suggest that a preference for some spawning habitat conditions could contribute to shape ABFT genetic structuring in the Mediterranean. However, further studies should be performed to assess to what extent ABFT spawning behaviour in the Mediterranean Sea can be affected by environmental variation.