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

Impact of land development along the western Mediterranean coast of Egypt regarding surfactant sources, interfering elements and ecotoxicity

The effect of recent land development of the western side of the Egyptian Mediterranean coast on the fates, behaviors, interactions, and ecotoxicology of surfactant (LAS), F, Br, B, Ca, Mg, and P was studied. Samples of seawater and sediments were collected from 15 stations at different depths representing, 5 perpendicular sectors. Elevated levels of LAS were identified in seawater columns in the El-Hamam (467.3 ± 220.8 μg/L) and El-Dabaa (480.0 ± 314.1 μg/L) stations. LAS homologue in sediment was in the range of 0.013-0.024, 0.042-0.184, 0.086-0.402, and 0.025-0.058 μg/g for C10, C11, C12, and C13, respectively. Studied parameters showed mixture risk characterization ratios RCRmix > 1 for algae, invertebrates, and fish in seawater and sediments, except for P, which showed low risk (RCRmix ≤ 1) in sediment. Acute relative contribution (RC) of LAS reflected that fish were the most sensitive species (RCFish = 48.5), followed by algae (RCAlgae = 44.4) and invertebrates (RCInvert = 7.1).

A comparative study of surfactant distribution and fate (western and eastern) Egyptian Mediterranean coasts focusing on its environmental toxicity

One of the most difficult-to-manage new contaminants constantly released into the environment is linear alkylbenzene sulphonate (LAS), an anionic surfactant. Significant volumes of LAS are received by the Mediterranean coast of Egypt. The current study is a comprehensive assessment of the environmental fate of the LAS 1505 km off the Mediterranean coast of Egypt in the fall of 2023 in order to track its geographic spread and eventual demise in the water column. Critical analysis of LAS revealed that it is vertically distributed in various ways according to sources, uses, production amounts, and salinity levels. The vertical variation of LAS can be explained by its amphiphilic structure. A significant increase in surfactant concentration (>300 μg/L) was recorded in 66% and 43% of the total samples, ranging from 301.128 to 455.36 and from 304.556 to 486.135 for the western and eastern sides along the Egyptian Mediterranean coast, respectively. Evaluation of the average acute and chronic risk quotient (RQ) along the investigated locations revealed that fish were the most susceptible to LAS in both long and short exposure periods. The presented results also indicated significant LAS toxicity to three trophic levels (RQ values > 1). LAS toxicity to marine organisms was greater in the western than in eastern coastal regions according to acute and chronic mixture risk characterization ratios (RCRmix). The three trophic levels in the study area had the following order of acute relative contribution (RC) to LAS toxicity: fish > invertebrates > algae. The ANOVA test results showed that in both the western and eastern regions, LAS varied significantly (p < 0.05) with salinity (1.04E-60 and 5.44E-42) and depth (6.02E-65 and 1.59E-47), respectively. In addition, a significant difference was observed using the ANOVA test between the eastern and western regions of the Egyptian Mediterranean coast.

Atlantic bonito (Sarda) genomic analysis reveals population differentiation across Northeast Atlantic and mediterranean locations: Implications for fishery management

The Atlantic bonito (Sarda, family Scombridae) is a pelagic species and one of the most exploited small tuna species. Despite its economic importance, biological information is scarce with no associated management directives. In this study, using a population genomic approach resulted in a lack of panmixia of two genetic pools with different effective population sizes (east central-tropical Atlantic and northeast Atlantic-Mediterranean) and an intermixing zone in Atlantic Morocco. The absence of genetic heterogeneity between the locations in Atlantic - Mediterranean transitional zone adds new evidence that challenges the Strait of Gibraltar as a phylogeographic barrier for marine pelagic species. These results are proposed to the ICCAT Commission to establish management areas for this species, although they are not consistent with the recently adopted pelagic ecoregions. Finally, a panel of highly informative SNPs was developed for efficient and low-cost population assignment and the analysis of unresolved population structures.

Population genetic insights into establishment, adaptation, and dispersal of the invasive quagga mussel across perialpine lakes

Human activities have facilitated the invasion of freshwater ecosystems by various organisms. Especially, invasive bivalves such as the quagga mussels, Dreissena bugensis, have the potential to alter ecosystem function as they heavily affect the food web. Quagga mussels occur in high abundance, have a high filtration rate, quickly spread within and between waterbodies via pelagic larvae, and colonize various substrates. They have invaded various waterbodies across the Northern Hemisphere. In Central Europe, they have invaded multiple large and deep perialpine lakes with first recordings in Lake Geneva in 2015 and 2016 in Lake Constance. In the deep perialpine lakes, quagga mussels quickly colonized the littoral zone but are also abundant deeper (>80 m), where they are often thinner and brighter shelled. We analysed 675 quagga mussels using ddRAD sequencing to gain in-depth insights into the genetic population structure of quagga mussels across Central European lakes and across various sites and depth habitats in Lake Constance. We revealed substantial genetic differentiation amongst quagga mussel populations from three unconnected lakes, and all populations showed high genetic diversity and effective population size. In Lake Constance, we detected no genetic differentiation amongst quagga mussels sampled across different sites and depth habitats. We also did not identify any convincing candidate loci evidential for adaptation along a depth gradient and a transplant experiment showed no indications of local adaptation to living in the deep based on investigating growth and survival. Hence, the shallow-water and the deep-water morphotypes seem to be a result of phenotypic plasticity rather than local adaptation to depth. In conclusion, our ddRAD approach revealed insight into the establishment of genetically distinct quagga mussel populations in three perialpine lakes and suggests that phenotypic plasticity and life history traits (broadcast spawner with high fecundity and dispersing pelagic larvae) facilitate the fast spread and colonization of various depth habitats by the quagga mussel.

Atlantic Oceanic Squids in the "Grey Speciation Zone"

Cryptic species complexes represent an important challenge for the adequate characterization of Earth's biodiversity. Oceanic organisms tend to have greater unrecognized cryptic biodiversity since the marine realm was often considered to lack hard barriers to genetic exchange. Here, we tested the effect of several Atlantic and Mediterranean oceanic barriers on 16 morphospecies of oceanic squids of the orders Oegopsida and Bathyteuthida using three mitochondrial and one nuclear molecular marker and five species delimitation methods. Number of species recognized within each morphospecies differed among different markers and analyses, but we found strong evidence of cryptic biodiversity in at least four of the studied species (Chtenopteryx sicula, Chtenopteryx canariensis, Ancistrocheirus lesueurii, and Galiteuthis armata). There were highly geographically structured units within Helicocranchia navossae that could either represent recently diverged species or population structure. Although the species studied here can be considered relatively passive with respect to oceanic currents, cryptic speciation patterns showed few signs of being related to oceanic currents. We hypothesize that the bathymetry of the egg masses and duration of the paralarval stage might influence the geographic distribution of oceanic squids. Because the results of different markers and different species delimitation methods are inconsistent and because molecular data encompassing broad geographic sampling areas for oceanic squids are scarce and finding morphological diagnostic characters for early life stages is difficult, it is challenging to assess the species boundaries for many of these species. Thus, we consider many to be in the "grey speciation zone." As many oceanic squids have cosmopolitan distributions, new studies combining genomic and morphological information from specimens collected worldwide are needed to correctly assess the actual oceanic squid biodiversity.

Insights on Pinna nobilis population genetic structure in the Aegean and Ionian Sea

The fan mussel Pinna nobilis Linnaeus, 1758 is an endemic species of the Mediterranean Sea, protected by international agreements. It is one of the largest bivalves in the world, playing an important role in the benthic communities; yet it has been recently characterized as Critically Endangered by the IUCN, due to mass mortality events. In this context, the assessment of the genetic variation of the remaining P. nobilis populations and the evaluation of connectivity among them are crucial elements for the conservation of the species. For this purpose, samples were collected from six regions of the Eastern Mediterranean Sea; the Islands of Karpathos, Lesvos and Crete; the Chalkidiki and Attica Peninsulas; and the Amvrakikos Gulf. Sampling was performed either by collecting tissue from the individuals or by using a non-invasive method, i.e., by scraping the inside of their shells aiming to collect their mucus and thus avoid stress induction to them. Conventional molecular techniques with the use of the COI and 16S rRNA mitochondrial markers were selected for the depiction of the intra-population genetic variability. The analyses included 104 samples from the present study and publicly available sequences of individuals across the whole Mediterranean Sea. The results of this work (a) suggest the use of eDNA as an efficient sampling method for protected bivalves and (b) shed light to the genetic structure of P. nobilis population in the Eastern Mediterranean; this latter knowledge might prove to be fundamental for the species conservation and hence the ecosystem resilience. The haplotype analyses reinforced the evidence that there is a certain degree of connectivity among the distinct regions of the Mediterranean; yet there is evidence of population distinction within the basin, namely between the Western and the Eastern basins. The combination of both genetic markers in the same analysis along with the inclusion of a large number of individuals produced more robust results, revealing a group of haplotypes being present only in the Eastern Mediterranean and providing insights for the species' most suitable conservation management.

Comparative population genomics of manta rays has global implications for management

Understanding population connectivity and genetic diversity is of fundamental importance to conservation. However, in globally threatened marine megafauna, challenges remain due to their elusive nature and wide-ranging distributions. As overexploitation continues to threaten biodiversity across the globe, such knowledge gaps compromise both the suitability and effectiveness of management actions. Here, we use a comparative framework to investigate genetic differentiation and diversity of manta rays, one of the most iconic yet vulnerable groups of elasmobranchs on the planet. Despite their recent divergence, we show how oceanic manta rays (Mobula birostris) display significantly higher heterozygosity than reef manta rays (Mobula alfredi) and that M. birostris populations display higher connectivity worldwide. Through inferring modes of colonization, we reveal how both contemporary and historical forces have likely influenced these patterns, with important implications for population management. Our findings highlight the potential for fisheries to disrupt population dynamics at both local and global scales and therefore have direct relevance for international conservation of marine species.

I Like This New Me: Unravelling Population Structure of Mediterranean Electric Rays and Taxonomic Uncertainties within Torpediniformes

The present study focused on the three species of electric rays known to occur in the Mediterranean Sea: Torpedo torpedo, Torpedo marmorata and Tetronarce nobiliana. Correct identification of specimens is needed to properly assess the impact of fisheries on populations and species. Unfortunately, torpedoes share high morphological similarities, boosting episodes of field misidentification. In this context, genetic data was used (1) to identify specimens caught during fishing operations, (2) to measure the diversity among and within these species, and (3) to shed light on the possible occurrence of additional hidden species in the investigated area. New and already published sequences of COI and NADH2 mitochondrial genes were analyzed, both at a small scale along the Sardinian coasts (Western Mediterranean) and at a large scale in the whole Mediterranean Sea. High levels of genetic diversity were found in Sardinian populations, being significantly different from other areas of the Eastern Mediterranean Sea due to the biotic and abiotic factors here discussed. Sardinian torpedoes can hence be indicated as priority populations/areas to be protected within the Mediterranean Sea. Moreover, sequence data confirmed that only the three species occur in the investigated area. The application of several 'species-delimitation' methods found evidence of cryptic species in the three species outside the Mediterranean Sea, as well as in other genera/families, suggesting the urgent need for future studies and a comprehensive revision of the order Torpediniformes for its effective conservation.

Matching an Old Marine Paradigm: Limitless Connectivity in a Deep-Water Fish over a Large Distance

Investigations of population structuring in wild species are fundamental to complete the bigger picture defining their ecological and biological roles in the marine realm, to estimate their recovery capacity triggered by human disturbance and implement more efficient management strategies for fishery resources. The Blackspot Seabream (Pagellus bogaraveo, Brünnich 1768) is a commercially valuable deep-water fish highly exploited over past decades. Considering its exploitation status, deepening the knowledge of intraspecific variability, genetic diversity, and differentiation using high-performing molecular markers is considered an important step for a more effective stock assessment and fishery management. With one of the largest efforts conceived of and completed by countries overlooking the Atlantic and Mediterranean coasts in recent years, a total of 320 individuals were collected from different fishing grounds in the Mediterranean Sea and Atlantic Ocean and analysed using 29 microsatellite loci. We applied multiple statistical approaches to investigate the species' connectivity and population structure across most of its described distribution area. Considering the incomplete knowledge regarding the migratory behaviour of adults, here we suggest the importance of egg and larval dispersal in sustaining the observed genetic connectivity on such a large geographical scale.

Exceptional population genomic homogeneity in the black brittle star Ophiocomina nigra (Ophiuroidea, Echinodermata) along the Atlantic-Mediterranean coast

The Atlantic-Mediterranean marine transition is characterised by strong oceanographic barriers and steep environmental gradients that generally result in connectivity breaks between populations from both basins and may lead to local adaptation. Here, we performed a population genomic study of the black brittle star, Ophiocomina nigra, covering most of its distribution range along the Atlantic-Mediterranean region. Interestingly, O. nigra is extremely variable in its coloration, with individuals ranging from black to yellow-orange, and different colour morphs inhabiting different depths and habitats. In this work, we used a fragment of the mitochondrial COI gene and 2,374 genome-wide ddRADseq-derived SNPs to explore: (a) whether the different colour morphs of O. nigra represent different evolutionary units; (b) the disruptive effects of major oceanographic fronts on its population structure; and (c) genomic signals of local adaptation to divergent environments. Our results revealed exceptional population homogeneity, barely affected by oceanographic fronts, with no signals of local adaptation nor genetic differentiation between colour morphs. This remarkable panmixia likely results from a long pelagic larval duration, a large effective population size and recent demographic expansions. Our study unveils an extraordinary phenotypic plasticity in O. nigra, opening further research questions on the ecological and molecular mechanisms underpinning coloration in Ophiuroidea.

Evaluating the taxonomic status of the large sized Tricolia Risso, 1826 in the Northeast Atlantic and Mediterranean Sea

Despite a long history of taxonomic studies on the genus Tricolia Risso, 1826, there is a shortfall on thorough systematic molecular reviews of the taxon from the NE Atlantic and Mediterranean coasts. Aiming to assess the genetic distinctness among morphospecies and the taxonomic status of currently accepted large sized species in these areas, we conducted a molecular phylogenetic analysis of the genus based on one mitochondrial (cox1) and two nuclear (28S and ITS2) markers. Seven Tricolia species were consistently retrieved in the analyses, including a new genetic lineage in the NE Atlantic designated as Tricolia sp. 1. Molecular analyses revealed that only one species, T. azorica, occurs in the NE Atlantic archipelagos. The sister taxa T. pullus (Mediterranean) and T. picta (NE Atlantic) should be classified as distinct species, instead of subspecies of the T. pullus group (sensu Gofas 1982). Tricolia miniata is also a complex of species in the Mediterranean and future studies across the distribution range are necessary to clarify its status.

Black Sea Turtle (Chelonia mydas agassizii) Life History in the Sanctuary of Colola Beach, Michoacan, Mexico

Sea turtles present strategies that have allowed them to survive and reproduce. They spend most of their lives in the sea, except when they emerge as hatchlings from the nest and when the adult females return to nest. Those moments of their life cycle are vital for their reproductive success, conservation, and knowledge of their biology. This study reports the life history traits exhibited by female black sea turtles from Colola Beach, Mexico using morphometric and reproductive data obtained during 15 sampling seasons (1985-2000, n = 1500). The results indicate that nesting females have a mean body size of 85.7 cm and reach sexual maturity at 24 years old at a minimum size of 68 cm. Females deposit a mean of 69.3 eggs per clutch, and the mean fecundity was 196.4 eggs per female per season. The remigration intervals of 3 and 5 years were the most frequent registered. The life history traits found in the black sea turtle population present the lowest values reported with respect to studies conducted in the Atlantic and Indo-Pacific green turtle populations, which supports the hypothesis that this population is recovering, since morphometric and reproductive data represent young nesting turtles.

Limited population genetic variation but pronounced seascape genetic structuring in populations of the Mediterranean mussel (Mytilus galloprovincialis) from the eastern Adriatic Sea

Population genetic analysis of variation at five neutral microsatellite loci for Mediterranean mussels (Mytilus galloprovincialis) from 18 sites along the eastern Adriatic Sea revealed little or no spatial variation. In contrast, seascape genetics analysis revealed a pronounced locus-specific gradient in allelic and genotypic frequencies across the study region. At a sixth locus, MGE7, the frequencies of two alleles, MGE7²⁴³ and MGE7²⁴⁹, were strongly associated, negatively and positively, respectively, with a single environmental variable - minimum salinity (minSAL). The frequency of the MGE7^243/243 homozygous genotype was strongly negatively associated with minSAL, whereas the frequencies of the MGE7^246/249 and the MGE7^249/249 genotypes were strongly positively correlated with minSAL. Interpretation of these pronounced gradients is confounded by the fact that minSAL and another environmental variable, maximum sea surface temperature (maxSST), are highly correlated (R = -.911) and are therefore not necessarily acting independently. BLAST searches of the MGE7 locus against M. galloprovincialis whole genome shotgun sequence returned an alignment with contig mg10_S01094 (accession UYJE01010330.1) and 7 predicted M. galloprovincialis proteins VDI82194.1 - VDI82200.1. Conserved domain searches revealed a similar structure to the transcriptional regulator Msx2-interacting protein. The BLASTp search also returned significant alignments to Msx2-interacting proteins in Mytilus coruscus, Crassostrea virginica, and Haliotis rubra. The existence of the MGE7 gradient highlights the role that environmental variation may play in retarding gene flow among wild M. galloprovincialis populations, and also how the success of collection of young mussels (spat) from one site and their transfer to another site (the farm) may be influenced by a single factor such as minSAL or maxSST on a localized scale.

Horizontal and vertical segregation of polycyclic aromatic hydrocarbons in the Egyptian Mediterranean coast

Egyptian Mediterranean coast receives significant amounts of polycyclic aromatic hydrocarbons (PAHs) from industrial exhausts, riverine inputs, maritime shipping and fishers, and oil and natural gas production and exploration. The present study considers the first exhaustive assessment for the dissolved PAHs along the Egyptian Mediterranean coast (Alexandria to Manzallah) to monitor their spatial distribution and investigate the effect of the marine currents and the role of microorganisms in their distribution. Surface water levels ranged between 124.97 and 301.02 ng L^-1 with an average 223.68 ± 41.11 ng L^-1. The distribution increases from west to east based on the water circulation in the Mediterranean Sea. The levels in near shore stations were lower than those of middle and onshore stations. The intensive existence of micro-organisms near shore stations consumes great part of PAHs, while this bio-remediation process decreases gradually away from the shoreline leaving relative high concentrations of dissolved PAHs in the middle and onshore stations. Middle and deep-water levels ranged between 312.75 and 1042.95 ng L^-1 with an average 633.47 ± 225.53 ng L^-1. Deeper waters showed higher PAHs concentrations where the average concentrations of 50 m stations (868.12 ± 138.35 ng L^-1) ˃ 30 m stations (629.49 ± 143.85 ng L^-1) ˃ 10 m stations (402.79 ± 59.46 ng L^-1). The wind-induced waves re-suspend rich PAHs sediment particles to increase its concentration in the water column. Carcinogenic toxic equivalent quotient (TEQ) for total detected PAHs in the middle and deep water represented more than double (75.46 ng TEQ L^-1) the value in the surface water (34.76 ng TEQ L^-1). The diagnostic ratios and principal component analysis indicated mainly pyrogenic origin in surface, middle, and deep waters.

Morphospecies and molecular diversity of ‘lace corals’: the genus Reteporella (Bryozoa: Cheilostomatida) in the central North Atlantic Azores Archipelago

Background

As in most bryozoans, taxonomy and systematics of species in the genus Reteporella Busk, 1884 (family Phidoloporidae) has hitherto almost exclusively been based on morphological characters. From the central North Atlantic Azores Archipelago, nine Reteporella species have historically been reported, none of which have as yet been revised. Aiming to characterise the diversity and biogeographic distribution of Azorean Reteporella species, phylogenetic reconstructions were conducted on a dataset of 103 Azorean Reteporella specimens, based on the markers cytochrome C oxidase subunit 1, small and large ribosomal RNA subunits. Morphological identification was based on scanning electron microscopy and complemented the molecular inferences.

Results

Our results reveal two genetically distinct Azorean Reteporella clades, paraphyletic to eastern Atlantic and Mediterranean taxa. Moreover, an overall concordance between morphological and molecular species can be shown, and the actual bryozoan diversity in the Azores is greater than previously acknowledged as the dataset comprises three historically reported species and four putative new taxa, all of which are likely to be endemic. The inclusion of Mediterranean Reteporella specimens also revealed new species in the Adriatic and Ligurian Sea, whilst the inclusion of additional phidoloporid taxa hints at the non-monophyly of the genus Reteporella.

Conclusion

Being the first detailed genetic study on the genus Reteporella, the high divergence levels inferred within the genus Reteporella and family Phidoloporidae calls for the need of further revision. Nevertheless, the overall concordance between morphospecies and COI data suggest the potential adequacy of a 3% cut-off to distinguish Reteporella species. The discovery of new species in the remote Azores Archipelago as well as in the well-studied Mediterranean Sea indicates a general underestimation of bryozoan diversity. This study constitutes yet another example of the importance of integrative taxonomical approaches on understudied taxa, contributing to cataloguing genetic and morphological diversity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02080-z.

Genomic basis for early-life mortality in sharpsnout seabream

Mortality at early life stages of fishes is common in nature and can be shaped by stochastic and selective processes. Selective mortality has rarely been assessed in natural conditions but can now be studied by combining genomic data with information on different life stages that realates to fitness. Here we investigate selective mortality between settlers and six-month survivors of the sharpsnout seabream by genotype-phenotype/environmental association studies in three localities along a geographic gradient. We gathered information on 105 individuals at 85,031 SNPs, obtained from individual based 2b-RAD libraries, as well as 9 phenotypic and environmental variables derived from individual otolith readings. We found common signals across localities for potential selection, such as lower survival rates for individuals hatching earlier, growing faster and experiencing higher temperatures during their planktonic phase. We identified 122 loci with parallel significant association to phenotypic and environmental variables. Importantly, one of these loci mapped to the exonic region of the il20rb, a gene involved in immune response, in the phylogenetically closest reference genome, showing parallel frequency changes in non-synonymous mutations in the three studied populations. Further temporal assessments are needed to understand how polymorphisms that are key to selective mortality are maintained.

Spatial coalescent connectivity through multi-generation dispersal modelling predicts gene flow across marine phyla

Gene flow governs the contemporary spatial structure and dynamic of populations as well as their long-term evolution. For species that disperse using atmospheric or oceanic flows, biophysical models allow predicting the migratory component of gene flow, which facilitates the interpretation of broad-scale spatial structure inferred from observed allele frequencies among populations. However, frequent mismatches between dispersal estimates and observed genetic diversity prevent an operational synthesis for eco-evolutionary projections. Here we use an extensive compilation of 58 population genetic studies of 47 phylogenetically divergent marine sedentary species over the Mediterranean basin to assess how genetic differentiation is predicted by Isolation-By-Distance, single-generation dispersal and multi-generation dispersal models. Unlike previous approaches, the latter unveil explicit parents-to-offspring links (filial connectivity) and implicit links among siblings from a common ancestor (coalescent connectivity). We find that almost 70 % of observed variance in genetic differentiation is explained by coalescent connectivity over multiple generations, significantly outperforming other models. Our results offer great promises to untangle the eco-evolutionary forces that shape sedentary population structure and to anticipate climate-driven redistributions, altogether improving spatial conservation planning.

Genetic Structure of the Endangered Coral Cladocora caespitosa Matches the Main Bioregions of the Mediterranean Sea

Population connectivity studies are a useful tool for species management and conservation planning, particular of highly threatened or endangered species. Here, we evaluated the genetic structure and connectivity pattern of the endangered coral Cladocora caespitosa across its entire distribution range in the Mediterranean Sea. Additionally, we examined the relative importance of sexual and asexual reproduction in the studied populations and their genetic diversity. A total of 541 individuals from 20 localities were sampled and analysed with 19 polymorphic microsatellite markers. Of the genotyped individuals, 482 (89%) had unique multilocus genotypes. Clonality percentages of the populations varied from 0% (in eight populations) to nearly 69% (in one population from Crete). A heterozygosity deficit and a high degree of inbreeding was the general trend in our data set. Population differentiation in C. caespitosa was characterised by significant pairwise FST values with lower ones observed at an intraregional scale and higher ones, between populations from different biogeographic regions. Genetic structure analyses showed that the populations are divided according to the three main sub-basins of the Mediterranean Sea: the Western (Balearic, Ligurian and Tyrrhenian seas), the Central (Adriatic and Ionian seas) and the Eastern (Levantine and Aegean seas), coinciding with previously described gene flow barriers. However, the three easternmost populations were also clearly separated from one another, and a substructure was observed for the other studied areas. An isolation-by-distance pattern was found among, but not within, the three main population groups. This substructure is mediated mainly by dispersal along the coastline and some resistance to larval movement through the open sea. Despite the low dispersal ability and high self-recruitment rate of C. caespitosa, casual dispersive events between regions seem to be enough to maintain the species’ considerable genetic diversity. Understanding the population connectivity and structure of this endangered scleractinian coral allows for more informed conservation decision making.

Differences in Air and Sea Surface Temperatures in the Northern and Southern Part of the Adriatic Sea

The paper compares air and sea surface temperatures in recent years on two islands in the Adriatic Sea. The data measured at the climatological station Krk on the island of Krk and the main meteorological station Lastovo on the island of Lastovo are used. The island of Krk is located in the north of the Adriatic Sea and Lastovo in the south. Since a significant increase in air and sea surface temperatures has been observed over the last thirty years, the goal is to establish how they reflect at these two stations, 313 km apart. The goal of the analysis is to monitor the changes in these two islands to reduce the negative impacts they may cause. The analysis of sea temperatures showed that global warming has a greater impact in the northern Adriatic than in the southern Adriatic. Air and sea surface temperatures have a faster upward trend on Krk than on Lastovo. Similar to the Mediterranean Sea, a positive trend was observed in the Adriatic Sea for both sea surface temperature and air temperature.

Development of 27 new microsatellite markers for the shanny Lipophrys pholis

BACKGROUND

The shanny Lipophrys pholis is an intertidal fish that is widely distributed throughout the Northeast Atlantic. Characterized by limited adult mobility and a long pelagic larval duration, the shanny stands as an ideal model to better understand larval dispersal and connectivity dynamics, which are critical parameters with implications for marine conservation and management.

METHODS AND RESULTS

To this aim, we developed 27 highly polymorphic microsatellite markers and characterized a population of 42 individuals, presenting an average allelic diversity of 20.1 alleles per locus and heterozygosity ranging from 0.619 to 1.

CONCLUSIONS

This set of newly developed microsatellite markers will be useful in providing critical insights into the processes which shape L. pholis gene flow and connectivity patterns and can be used to investigate local parentage lineages.

Interannual fluctuations in connectivity among crab populations (Liocarcinus depurator) along the Atlantic-Mediterranean transition

An interesting evolutionary question that still remains open is the connectivity between marine populations. Marine currents can favour the dispersal of larvae or adults, but they can also produce eddies and gyres generating oceanographic fronts, thus limiting gene flow. To address this subject, we selected the Atlantic-Mediterranean transition, where several fronts are located: Gibraltar Strait (GS), Almeria-Oran Front (AOF) and Ibiza Channel (IC). Seven populations of the marine crab Liocarcinus depurator (Cadiz, West and East Alboran, Alacant, Valencia, Ebro Delta and North Catalonia) located along this transition were analysed in six consecutive years (2014–2019) using a fragment of the COI (Cytochrome Oxidase subunit I) gene. All sequences (966) belonged to two well defined haplogroups: ATL (most abundant in Atlantic waters) and MED (predominant in Mediterranean waters). Following a geographic variation, the frequency of ATL decreased significantly from Cadiz to North Catalonia. However, this variation presented steps due to the effect of oceanographic restrictions/fronts. Significant effects were recorded for GS (2015, 2017, 2018 and 2019), AOF (all years except 2018) and IC (2016). The intensity and precise location of these fronts changed over time. Multivariate analyses distinguished three main population groups: Cadiz, Alboran Sea and the remaining Mediterranean populations. These findings could be relevant to properly define Marine Protected Areas and for conservation and fisheries policies.

Clarifying the taxonomy of some cryptic blennies (Blenniidae) in their native and introduced range

Omobranchus punctatus is native to the Indo-Pacific region and invasive in the Atlantic region, currently being considered one of the most widely distributed blenny species. However, recent molecular studies indicated that O. punctatus is a complex of species, with three divergent mtDNA lineages identified to date, stressing the need for a taxonomic revision. In this study, we used an integrative approach, combining morphological and genetic data, to shed light on the taxonomy and distribution of O. punctatus. Moreover, we provide the first genetic records of introduced populations in Brazil and discuss the introduction pattern of this species in this region. Morphological data shows that O. punctatus consists of at least five distinct and geographically restricted species: O. punctatus sensu stricto, O. dispar, O. sewalli, O. cf. kochi, and O. cf. japonicus. Species delimitation analyses performed using the mtDNA data available confirmed that O. punctatus sensu stricto, O. dispar and O. sewalli correspond to different species that started to diverge about 2.6 Mya. Furthermore, O. sewalli was identified as the invasive species colonizing Atlantic shores. The existence of historical oceanographic barriers, such as the emergence of the Sunda Shelf in the Eastern Indian Ocean during the Pleistocene, and the biological traits of these blennies are the most likely factors responsible for their genetic differentiation and subsequent speciation.

Genomic distances reveal relationships of wild and cultivated beets

Cultivated beets (Beta vulgaris ssp. vulgaris), including sugar beet, rank among the most important crops. The wild ancestor of beet crops is the sea beet Beta vulgaris ssp. maritima. Species and subspecies of wild beets are readily crossable with cultivated beets and are thus available for crop improvement. To study genomic relationships in the genus Beta, we sequence and analyse 606 beet genomes, encompassing sugar beet, sea beet, B. v. adanensis, B. macrocarpa, and B. patula. We observe two genetically distinct groups of sea beets, one from the Atlantic coast and the other from the Mediterranean area. Genomic comparisons based on k-mers identify sea beets from Greece as the closest wild relatives of sugar beet, suggesting that domestication of the ancestors of sugar beet may be traced to this area. Our work provides comprehensive insight into the phylogeny of wild and cultivated beets and establishes a framework for classification of further accessions of unknown (sub-)species assignment.

Mitonuclear genetic patterns of divergence in the marbled crab, Pachygrapsus marmoratus (Fabricius, 1787) along the Turkish seas

Biogeographical transition zones present good opportunities for studying the effect of the past ice ages on genetic structure of species because secondary contact zones of post-glacial lineages can be formed. In this study, we investigated the population genetic structure of the marbled rock crab, Pachygrapsus marmoratus along the coasts of Turkey. We genotyped 334 individuals from the Black Sea, the Turkish Straits System (TSS), the Aegean, and the Eastern Mediterranean basins. In order to reveal its evolutionary history and its population connectivity, we used mitochondrial CO1 region and five microsatellite loci. CO1 analyzes also included 610 additional samples from Genbank, which covered most of its distribution range. Both microsatellites and mtDNA showed decreased diversity in sampling sites of the TSS and the Black Sea as compared to those along the Aegean and the Levantine coasts. There is an especially strong geographical pattern in distributions of haplotypes in mtDNA, most probably as a result of genetic drift in the Black Sea and the Sea of Marmara (SoM). Microsatellite data analyses revealed two genetically distinct clusters of P. marmoratus (clusters C and M). While individuals belonging to cluster C are present in all the sampling locations, those belonging to cluster M are only detected along the Mediterranean coasts including the Aegean and the Levantine basins. These clusters shared similar haplotypes in the Mediterranean. Haplotypes of two sympatric clusters could be similar due to incomplete lineage sorting of ancestral polymorphisms. In order to retrieve the complex demographic history and to investigate evolutionary processes resulting in sympatric clusters in the Aegean Sea and the Levantine basin, mitochondrial markers with faster mutation rates than CO1 and/or SNP data will be useful.

Climate differently influences the genomic patterns of two sympatric marine fish species

Climate influences population genetic variation in marine species. Capturing these impacts remains challenging for marine fishes which disperse over large geographical scales spanning steep environmental gradients. It requires the extensive spatial sampling of individuals or populations, representative of seascape heterogeneity, combined with a set of highly informative molecular markers capable of revealing climatic-associated genetic variations. We explored how space, dispersal and environment shape the genomic patterns of two sympatric fish species in the Mediterranean Sea, which ranks among the oceanic basins most affected by climate change and human pressure. We hypothesized that the population structure and climate-associated genomic signatures of selection would be stronger in the less mobile species, as restricted gene flow tends to facilitate the fixation of locally adapted alleles. To test our hypothesis, we genotyped two species with contrasting dispersal abilities: the white seabream Diplodus sargus and the striped red mullet Mullus surmuletus. We collected 823 individuals and used genotyping by sequencing (GBS) to detect 8,206 single nucleotide polymorphisms (SNPs) for the seabream and 2,794 for the mullet. For each species, we identified highly differentiated genomic regions (i.e. outliers) and disentangled the relative contribution of space, dispersal and environmental variables (climate, marine primary productivity) on the outliers' genetic structure to test the prevalence of gene flow and local adaptation. We observed contrasting patterns of gene flow and adaptive genetic variation between the two species. The seabream showed a distinct Alboran sea population and panmixia across the Mediterranean Sea. The mullet revealed additional differentiation within the Mediterranean Sea that was significantly correlated to summer and winter temperatures, as well as marine primary productivity. Functional annotation of the climate-associated outlier SNPs then identified candidate genes involved in heat tolerance that could be examined to further predict species' responses to climate change. Our results illustrate the key steps of a comparative seascape genomics study aiming to unravel the evolutionary processes at play in marine species, to better anticipate their response to climate change. Defining population adaptation capacities and environmental niches can then serve to incorporate evolutionary processes into species conservation planning.

Regional patterns of δ 13C and δ 15N for European common cuttlefish (Sepia officinalis) throughout the Northeast Atlantic Ocean and Mediterranean Sea

Stable isotope compositions of carbon and nitrogen (expressed as δ 13C and δ 15N) from the European common cuttlefish (Sepia officinalis) were measured in order to evaluate the utility of using these natural tracers throughout the Northeast Atlantic Ocean and Mediterranean Sea (NEAO-MS). Mantle tissue was obtained from S. officinalis collected from 11 sampling locations spanning a wide geographical coverage in the NEAO-MS. Significant differences of both δ 13C and δ 15N values were found among S. officinalis samples relative to sampling location. δ 13C values did not show any discernable spatial trends; however, a distinct pattern of lower δ 15N values in the Mediterranean Sea relative to the NEAO existed. Mean δ 15N values of S. officinalis in the Mediterranean Sea averaged 2.5‰ lower than conspecifics collected in the NEAO and showed a decreasing eastward trend within the Mediterranean Sea with the lowest values in the most eastern sampling locations. Results suggest δ 15N may serve as a useful natural tracer for studies on the population structure of S. officinalis as well as other marine organisms throughout the NEAO-MS.

Genetic homogeneity, lack of larvae recruitment, and clonality in absence of females across western Mediterranean populations of the starfish Coscinasterias tenuispina

We here analysed the populations' genetic structure of Coscinasterias tenuispina, an Atlantic-Mediterranean fissiparous starfish, focusing on the western Mediterranean, to investigate: the distribution and prevalence of genetic variants, the relative importance of asexual reproduction, connectivity across the Atlantic-Mediterranean transition, and the potential recent colonisation of the Mediterranean Sea. Individuals from 11 Atlantic-Mediterranean populations of a previous study added to 172 new samples from five new W Mediterranean sites. Individuals were genotyped at 12 microsatellite loci and their gonads histologically analysed for sex determination. Additionally, four populations were genotyped at two-time points. Results demonstrated genetic homogeneity and low clonal richness within the W Mediterranean, due to the dominance of a superclone, but large genetic divergence with adjacent areas. The lack of new genotypes recruitment over time, and the absence of females, confirmed that W Mediterranean populations were exclusively maintained by fission and reinforced the idea of its recent colonization. The existence of different environmental conditions among basins and/or density-depend processes could explain this lack of recruitment from distant areas. The positive correlation between clonal richness and heterozygote excess suggests that most genetic diversity is retained within individuals in the form of heterozygosity in clonal populations, which might increase their resilience.

Genetic connectivity of two marine gastropods in the Mediterranean Sea: seascape genetics reveals species-specific oceanographic drivers of gene flow

Oceanographic features such as currents, waves, temperature and salinity, together with life history traits, control patterns and rates of gene flow and contribute to shaping the population genetic structure of marine organisms. Seascape genetics is an emerging discipline that adopts a spatially explicit approach to examine biotic and abiotic factors that drive gene flow in marine environments. In this study, we examined factors that contribute to genetic differentiation in two coastal Mediterranean gastropods whose geographical ranges overlap but which inhabit different environments. The two species differ in several life history traits and in their dispersal capabilities. Genetic differentiation was relatively low for the trochid species Gibbula divaricata (F_ST  =0.059), and high for the vermetid species Dendropoma lebeche (F_ST  =0.410). Salinity emerged as the most important variable explaining the genetic structure of both species; sea surface temperature was also important for G. divaricata. For the more sessile D. lebeche, the coastline was predicted to provide important pathways for stepping-stone connectivity and gene flow. Our results provide a greater understanding of the factors influencing marine population connectivity, which may be useful to guide marine conservation and management in the Mediterranean.

Genetic hypervariability of a Northeastern Atlantic venomous rockfish

BACKGROUND

Understanding the interplay between climate and current and historical factors shaping genetic diversity is pivotal to infer changes in marine species range and communities' composition. A phylogeographical break between the Atlantic and the Mediterranean has been documented for several marine organisms, translating into limited dispersal between the two basins.

METHODS

In this study, we screened the intraspecific diversity of 150 individuals of the Madeira rockfish (Scorpaena maderensis) across its distributional range (seven sampling locations in the Atlantic and Mediterranean basins) using the mitochondrial control region and the nuclear S7 first intron.

RESULTS

The present work is the most comprehensive study done for this species, yielding no genetic structure across sampled locations and no detectable Atlantic-Mediterranean break in connectivity. Our results reveal deep and hyper-diverse bush-like genealogies with large numbers of singletons and very few shared haplotypes. The genetic hyper-diversity found for the Madeira rockfish is relatively uncommon in rocky coastal species, whose dispersal capability is limited by local oceanographic patterns. The effect of climate warming on the distribution of the species is discussed.

Weak biodiversity connectivity in the European network of no-take marine protected areas

The need for international cooperation in marine resource management and conservation has been reflected in the increasing number of agreements aiming for effective and well-connected networks of Marine Protected Areas (MPAs). However, the extent to which individual MPAs are connected remains mostly unknown. Here, we use a biophysical model tuned with empirical data on species dispersal ecology to predict connectivity of a vast spectrum of biodiversity in the European network of marine reserves (i.e., no-take MPAs). Our results highlight the correlation between empirical propagule duration data and connectivity potential and show weak network connectivity and strong isolation for major ecological groups, resulting from the lack of direct connectivity corridors between reserves over vast regions. The particularly high isolation predicted for ecosystem structuring species (e.g., corals, sponges, macroalgae and seagrass) might potentially undermine biodiversity conservation efforts if local retention is insufficient and unmanaged populations are at risk. Isolation might also be problematic for populations' persistence in the light of climate change and expected species range shifts. Our findings provide novel insights for management directives, highlighting the location of regions requiring additional marine reserves to function as stepping-stone connectivity corridors.

Genetic landscape clustering of a large DNA barcoding data set reveals shared patterns of genetic divergence among freshwater fishes of the Maroni Basin

The Maroni is one of the most speciose basins of the Guianas and hosts a megadiverse freshwater fish community. Although taxonomic references based on morphological identification exist for both the Surinamese and Guianese parts of the basin, there are still taxonomic uncertainties concerning the status of several species. We used COI sequences of 1284 fish in conjunction with morphological and biogeographical evidence to assist with species delineation and discovery in order to validate and standardize the current taxonomy. This resulted in a final DNA barcode data set of 199 fish species (125 genera, 36 families and eight orders; 68.86% of strictly freshwater fishes from the basin), among which 25 are new putative candidate species flagged as requiring taxonomic update. DNA barcoding delineation through Barcode Index Numbers (BINs) revealed further cryptic diversity (230 BINs in total). To explore global genetic patterns across the basin, genetic divergence landscapes were computed for 128 species, showing a global trend of high genetic divergence between the Surinamese southwest (Tapanahony and Paloemeu), the Guianese southeast (Marouini, Litany, Tampok, etc.), and the river outlet in the north. This could be explained by lower levels of connectivity between these three main areas and/or the exchange of individuals between these areas and the neighbouring basins. A new method of ordination of genetic landscapes successfully assigned species into cluster groups based on their respective pattern of genetic divergence across the Maroni Basin: genetically homogeneous species were effectively discriminated from species showing high spatial genetic fragmentation and possible lower capacity for dispersal.

Molecular taxonomy confirms that the northeastern Atlantic and Mediterranean Sea harbor a single lancelet, Branchiostoma lanceolatum (Pallas, 1774) (Cephalochordata: Leptocardii: Branchiostomatidae)

Branchiostomatidae (lancelets or amphioxus) comprises about 30 species, several of which are well-established models in evolutionary development. Our zoological and ecological knowledge of the family is nonetheless limited. Despite evident differences can be found among known populations, the taxonomy of Branchiostoma lanceolatum (type species of the genus Branchiostoma) has never been investigated with modern methods through its range in the northeastern Atlantic and Mediterranean Sea. We address this via a multilocus molecular approach and comparing specimens collected from different European populations. Results obtained here confirm the presence of a single species inhabiting the range between the topotypical localities of B. lanceolatum (Atlantic Ocean) and of its junior synonym B. lubricum (Mediterranean Sea), without evincing geographical structure between populations. This suggests that environment most likely drives the characteristics observed in different geographic areas. The long larval phase and the slow mutation rate in lancelets may have played a key role in the evolutionary history of this iconic species.

Drivers for genetic structure at different geographic scales for Pacific red snapper (Lutjanus peru) and yellow snapper (Lutjanus argentiventris) in the tropical eastern Pacific

The tropical eastern Pacific (TEP) is a highly dynamic region and a model system to study how habitat discontinuities affect the distribution of shorefishes, particularly for species that display ontogenetic habitat shifts, including snappers (Lutjanidae). To evaluate the genetic structure of the Pacific red snapper (Lutjanus peru) and the yellow snapper (Lutjanus argentiventris) throughout their distribution range along the TEP, 13 and 11 microsatellite loci were analysed, respectively. The genetic diversity of L. peru (N = 446) and L. argentiventris (N = 170) was evaluated in 10 and 5 localities, respectively, showing slightly higher but non-significant values in the Gulf of California for both species. The genetic structure analysis identified the presence of significant genetic structure in both species, but the locations of the identified barriers for the gene flow differed between species. The principal driver for the genetic structure at large scales >2500 km was isolation by distance. At smaller scales (<250 km), the habitat discontinuity for juveniles and adults and the environmental differences throughout the distribution range represented potential barriers to gene flow between populations for both species.

Warming intensifies the interaction between the temperate seagrass Posidonia oceanica and its dominant fish herbivore Sarpa salpa

Apart from directly influencing individual life histories of species, climate change is altering key biotic interactions as well, causing community processes to unravel. With rising temperatures, disruptions to producer-consumer relationships can have major knock-on effects, particularly when the producer is a habitat-forming species. We studied how sea surface temperature (SST) modifies multiple pathways influencing the interaction between the foundational seagrass species, Posidonia oceanica, and its main consumer, the fish Sarpa salpa in the Mediterranean Sea. We used a combination of a field-based temperature gradient approaches and experimental manipulations to assess the effect of temperature on seagrass performance (growth) and fish early life history (larval development) as well as on the interaction itself (seagrass palatability and fish foraging activity). Within the range of temperatures assessed, S. salpa larvae grew slightly faster at warmer conditions but maintained their settlement size, resulting in a relatively small reduction in pelagic larval duration (PLD) and potentially reducing dispersion. Under warmer conditions (>24 °C), P. oceanica reduced its growth rate considerably and seemed to display fewer deterring mechanisms as indicated by a disproportionate consumption in choice experiments. However, our field-based observations along the temperature gradient showed no change in fish foraging time, or in other aspects of feeding behaviour. As oceans warm, our results indicate that, while S. salpa may show little change in early life history, its preference towards P. oceanica might increase, which, together with reduced seagrass growth, could considerably intensify the strength of herbivory. It is unclear if P. oceanica meadows can sustain such an intensification, but it will clearly add to the raft of pressures this threatened ecosystem already faces from global and local environmental change.

Reconciling seascape genetics and fisheries science in three codistributed flatfishes

Uncertainty hampers innovative mixed-fisheries management by the scales at which connectivity dynamics are relevant to management objectives. The spatial scale of sustainable stock management is species-specific and depends on ecology, life history and population connectivity. One valuable approach to understand these spatial scales is to determine to what extent population genetic structure correlates with the oceanographic environment. Here, we compare the level of genetic connectivity in three codistributed and commercially exploited demersal flatfish species living in the North East Atlantic Ocean. Population genetic structure was analysed based on 14, 14 and 10 neutral DNA microsatellite markers for turbot, brill and sole, respectively. We then used redundancy analysis (RDA) to attribute the genetic variation to spatial (geographical location), temporal (sampling year) and oceanographic (water column characteristics) components. The genetic structure of turbot was composed of three clusters and correlated with variation in the depth of the pycnocline, in addition to spatial factors. The genetic structure of brill was homogenous, but correlated with average annual stratification and spatial factors. In sole, the genetic structure was composed of three clusters, but was only linked to a temporal factor. We explored whether the management of data poor commercial fisheries, such as in brill and turbot, might benefit from population-specific information. We conclude that the management of fish stocks has to consider species-specific genetic structures and may benefit from the documentation of the genetic seascape and life-history traits.

Temperature reduces fish dispersal as larvae grow faster to their settlement size

As species struggle to cope with rising ocean temperatures, temperate marine assemblages are facing major reorganization. Many benthic species have a brief but critical period dispersing through the plankton, when they are particularly susceptible to variations in temperature. Impacts of rising temperatures can thus ripple through the population with community-wide consequences. However, responses are highly species-specific, making it difficult to discern assemblage-wide patterns in the life histories of different fish species. Here, we evaluate the responses to temperature in the early life histories of several fish species using otolith reconstructive techniques. We also assess the consequences of future warming scenarios to this assemblage. We sampled recent settlers of nine common species across a temperature gradient in the Mediterranean Sea and obtained environmental data for each individual. Using otolith microstructure, we measured early life traits including pelagic larval duration (PLD), growth rate, settlement size, hatching and settlement dates. We used a GLM framework to examine how environmental variables influenced early life-history parameters. We show that increasing temperature results in considerable reduction in the dispersal potential of temperate fish. We find a nearly universal, assemblage-wide decline in pelagic larval duration (PLD) of between 10% and 25%. This was because, with increasing temperature, larvae grew quicker to their settlement size. Settlement size itself was less affected by temperature and appears to be an ontogenetically fixed process. Given current estimates of ocean warming, there could be an assemblage-wide reduction in larval dispersal of up to 50 km across the Mediterranean, reducing connectivity and potentially isolating populations as waters warm.

Chaotic genetic structure and past demographic expansion of the invasive gastropod Tritia neritea in its native range, the Mediterranean Sea

To better predict population evolution of invasive species in introduced areas it is critical to identify and understand the mechanisms driving genetic diversity and structure in their native range. Here, we combined analyses of the mitochondrial COI gene and 11 microsatellite markers to investigate both past demographic history and contemporaneous genetic structure in the native area of the gastropod Tritia neritea, using Bayesian skyline plots (BSP), multivariate analyses and Bayesian clustering. The BSP framework revealed population expansions, dated after the last glacial maximum. The haplotype network revealed a strong geographic clustering. Multivariate analyses and Bayesian clustering highlighted the strong genetic structure at all scales, between the Black Sea and the Adriatic Sea, but also within basins. Within basins, a random pattern of genetic patchiness was observed, suggesting a superimposition of processes involving natural biological effects (no larval phase and thus limited larval dispersal) and putative anthropogenic transport of specimens. Contrary to the introduced area, no isolation-by-distance patterns were recovered in the Mediterranean or the Black Seas, highlighting different mechanisms at play on both native and introduced areas, triggering unknown consequences for species’ evolutionary trajectories. These results of Tritia neritea populations on its native range highlight a mixture of ancient and recent processes, with the effects of paleoclimates and life history traits likely tangled with the effects of human-mediated dispersal.

Population Genetic Diversity of Two Marine Gobies (Gobiiformes: Gobiidae) from the North-Eastern Atlantic and the Mediterranean Sea

Gobies (Gobiiformes: Gobiidae) are the most species-rich family of fishes in general, and the most abundant fish group in the European seas. Nonetheless, our knowledge on many aspects of their biology, including the population genetic diversity, is poor. Although barriers to gene flow are less apparent in the marine environment, the ocean is not a continuous habitat, as has been shown by studies on population genetics of various marine biota. For the first time, European marine goby species which cannot be collected by common fishery techniques were studied. The population genetic structure of two epibenthic species, Gobius geniporus and Gobius cruentatus, from seven localities across their distribution ranges was assessed, using one mitochondrial (cytochrome b) and one nuclear gene (first intron of ribosomal protein gene S7). Our results showed that there is a great diversity of haplotypes of mitochondrial gene cytochrome b in both species at all localities. Global fixation indices (FST) indicated a great differentiation of populations in both studied gobies. Our results did not show a geographic subdivision to individual populations. Instead, the data correspond with the model of migration which allow divergence and recurrent migration from the ancestral population. The estimated migration routes coincide with the main currents in the studied area. This matches well the biology of the studied species, with adults exhibiting only short-distance movements and planktonic larval stages.

Investigation of mechanisms underlying chaotic genetic patchiness in the intertidal marbled crab Pachygrapsus marmoratus (Brachyura: Grapsidae) across the Ligurian Sea

Background

Studies on marine community dynamics and population structures are limited by the lack of exhaustive knowledge on the larval dispersal component of connectivity. Genetic data represents a powerful tool in understanding such processes in the marine realm. When dealing with dispersion and connectivity in marine ecosystems, many evidences show patterns of genetic structure that cannot be explained by any clear geographic trend and may show temporal instability. This scenario is usually referred to as chaotic genetic patchiness, whose driving mechanisms are recognized to be selection, temporal shifts in local population dynamics, sweepstakes reproductive success and collective dispersal.

In this study we focused on the marbled crab Pachygrapsus marmoratus that inhabits the rocky shores of the Mediterranean Sea, Black Sea and East Atlantic Ocean, and disperses through planktonic larvae for about 1 month. P. marmoratus exhibits unexpectedly low connectivity levels at local scale, although well-defined phylogeographic patterns across the species’ distribution range were described. This has been explained as an effect of subtle geographic barriers or due to sweepstake reproductive success. In order to verify a chaotic genetic patchiness scenario, and to explore mechanisms underlying it, we planned our investigation within the Ligurian Sea, an isolated basin of the western Mediterranean Sea, and we genotyped 321 individuals at 11 microsatellite loci.

Results

We recorded genetic heterogeneity among our Ligurian Sea samples with the occurrence of genetic clusters not matching the original populations and a slight inter-population divergence, with the geographically most distant populations being the genetically most similar ones. Moreover, individuals from each site were assigned to all the genetic clusters. We also recorded evidences of self-recruitment and a higher than expected within-site kinship.

Conclusions

Overall, our results suggest that the chaotic genetic patchiness we found in P. marmoratus Ligurian Sea populations is the result of a combination of differences in reproductive success, en masse larval dispersion and local larval retention. This study defines P. marmoratus as an example of marine spawner whose genetic pool is not homogenous at population level, but rather split in a chaotic mosaic of slightly differentiated genetic patches derived from complex and dynamic ecological processes.

Global patterns of population genetic differentiation in seed plants

Evaluating the factors that drive patterns of population differentiation in plants is critical for understanding several biological processes such as local adaptation and incipient speciation. Previous studies have given conflicting results regarding the significance of pollination mode, seed dispersal mode, mating system, growth form and latitudinal region in shaping patterns of genetic structure, as estimated by F_ST values, and no study to date has tested their relative importance together across a broad scale. Here, we assembled a 337-species data set for seed plants from publications with data on F_ST from nuclear markers and species traits, including variables pertaining to the sampling scheme of each study. We used species traits, while accounting for sampling variables, to perform phylogenetic multiple regressions. Results demonstrated that F_ST values were higher for tropical, mixed-mating, non-woody species pollinated by small insects, indicating greater population differentiation, and lower for temperate, outcrossing trees pollinated by wind. Among the factors we tested, latitudinal region explained the largest portion of variance, followed by pollination mode, mating system and growth form, while seed dispersal mode did not significantly relate to F_ST . Our analyses provide the most robust and comprehensive evaluation to date of the main ecological factors predicted to drive population differentiation in seed plants, with important implications for understanding the basis of their genetic divergence. Our study supports previous findings showing greater population differentiation in tropical regions and is the first that we are aware of to robustly demonstrate greater population differentiation in species pollinated by small insects.

Genetic diversity and population structure of the rockpool shrimp Palaemon elegans based on microsatellites: evidence for a cryptic species and differentiation across the Atlantic-Mediterranean transition

The rockpool shrimp Palaemon elegans is an ecologically important crustacean species within the European coastline fauna. In the present study, genetic diversity and population structure and connectivity were assessed by examining 21 polymorphic microsatellite loci at 13 sampling sites located along the natural distribution range. All localities showed similar levels of genetic variability. Significant deficits of heterozygosity were recorded, most likely due to the presence of null alleles. Genetic structure analyses revealed two clearly genetically distinct groups within P. elegans but without following any geographical or oceanographic basis. Thus, our results provided nuclear evidence for the existence of a Mediterranean cryptic species within P. elegans, highlighting the need to revise its taxonomic status. Regarding P. elegans sensu stricto, population structuring was reported across the Atlantic-Mediterranean transition area, where the Almería-Orán Front restricts the gene flow between the Atlantic and the Mediterranean population. Moreover, while population connectivity was suggested between all Mediterranean localities, some substructure was found within the Atlantic group. Canary Islands exhibited a weak but significant genetic differentiation from all Atlantic mainland localities, consistent with the isolation-by-distance pattern detected throughout the Atlantic population. Overall, all these findings provided new insights into the population biology of P. elegans complex.

Towards an Ecosystem-Based Marine Spatial Planning in the deep Mediterranean Sea

The deep sea covers about 79% of the Mediterranean basin, including habitats potentially able to deliver multiple ecosystem services and numerous resources of high economic value. Thus, the deep Mediterranean Sea represents an important frontier for marine resources exploitation, which is embedded within the European Blue Growth Strategy goals and agendas. The deep sea is crucial for the ecological functioning of the entire basin. For this reason, the deep Mediterranean deserves protection from the potential cumulative impacts derived from existent and developing human activities. Marine Spatial Planning (MSP) has been identified as key instrument for spatially allocating maritime uses in the sea space avoiding spatial conflicts between activities, and between activities and the environment. Indeed, MSP incorporates the ecosystem-based approach (EB-MSP) to balance both socio-economic and environmental objectives, in line with the Maritime Spatial Planning Directive and the Marine Strategy Framework Directive. Despite MSP is under implementation in Europe, the Directive is not applied yet for the managing and monitoring of the environmental status of the deep sea. In the Mediterranean, deep areas fall both in internal and territorial waters, and in High Seas, and its management framework turns out to be complicated. Moreover, a certain level of cumulative impacts in the deep Mediterranean has been already identified and likely underestimated because of paucity of knowledge related with deep-sea ecosystems. Thus, the implementation of scientific knowledge and the establishment of a sustainable management regime of deep-sea resources and space are urgent. This study aims at reflecting on the best available ecological knowledge on the deep Mediterranean to incorporate conservation objectives in EB-MSP. We propose a framework to include key ecological principles in the relevant phases of any EB-MSP processes taking in consideration existing socio-economic and conservation scenarios in the region. We add the uncertainty principle to reflect on the still unexplored and missing knowledge related to the deep Mediterranean. Here, we resume some guidelines to overcome limits and bottlenecks while ensuring protection of deep-sea ecosystems and resources in the Mediterranean Sea.

Comparative genomics approach to evolutionary process connectivity

The influence of species life history traits and historical demography on contemporary connectivity is still poorly understood. However, these factors partly determine the evolutionary responses of species to anthropogenic landscape alterations. Genetic connectivity and its evolutionary outcomes depend on a variety of spatially dependent evolutionary processes, such as population structure, local adaptation, genetic admixture, and speciation. Over the last years, population genomic studies have been interrogating these processes with increasing resolution, revealing a large diversity of species responses to spatially structured landscapes. In parallel, multispecies meta-analyses usually based on low-genome coverage data have provided fundamental insights into the ecological determinants of genetic connectivity, such as the influence of key life history traits on population structure. However, comparative studies still lack a thorough integration of macro- and micro-evolutionary scales to fully realize their potential. Here, I present how a comparative genomics framework may provide a deeper understanding of evolutionary process connectivity. This framework relies on coupling the inference of long-term demographic and selective history with an assessment of the contemporary consequences of genetic connectivity. Standardizing this approach across several species occupying the same landscape should help understand how spatial environmental heterogeneity has shaped the diversity of historical and contemporary connectivity patterns in different taxa with contrasted life history traits. I will argue that a reasonable amount of genome sequence data can be sufficient to resolve and connect complex macro- and micro-evolutionary histories. Ultimately, implementing this framework in varied taxonomic groups is expected to improve scientific guidelines for conservation and management policies.

From metabarcoding to metaphylogeography: separating the wheat from the chaff

Metabarcoding is by now a well-established method for biodiversity assessment in terrestrial, freshwater, and marine environments. Metabarcoding data sets are usually used for α- and β-diversity estimates, that is, interspecies (or inter-MOTU [molecular operational taxonomic unit]) patterns. However, the use of hypervariable metabarcoding markers may provide an enormous amount of intraspecies (intra-MOTU) information-mostly untapped so far. The use of cytochrome oxidase (COI) amplicons is gaining momentum in metabarcoding studies targeting eukaryote richness. COI has been for a long time the marker of choice in population genetics and phylogeographic studies. Therefore, COI metabarcoding data sets may be used to study intraspecies patterns and phylogeographic features for hundreds of species simultaneously, opening a new field that we suggest to name metaphylogeography. The main challenge for the implementation of this approach is the separation of erroneous sequences from true intra-MOTU variation. Here, we develop a cleaning protocol based on changes in entropy of the different codon positions of the COI sequence, together with co-occurrence patterns of sequences. Using a data set of community DNA from several benthic littoral communities in the Mediterranean and Atlantic seas, we first tested by simulation on a subset of sequences a two-step cleaning approach consisting of a denoising step followed by a minimal abundance filtering. The procedure was then applied to the whole data set. We obtained a total of 563 MOTUs that were usable for phylogeographic inference. We used semiquantitative rank data instead of read abundances to perform AMOVAs and haplotype networks. Genetic variability was mainly concentrated within samples, but with an important between seas component as well. There were intergroup differences in the amount of variability between and within communities in each sea. For two species, the results could be compared with traditional Sanger sequence data available for the same zones, giving similar patterns. Our study shows that metabarcoding data can be used to infer intra- and interpopulation genetic variability of many species at a time, providing a new method with great potential for basic biogeography, connectivity and dispersal studies, and for the more applied fields of conservation genetics, invasion genetics, and design of protected areas.