Pillars of Hercules: is the Atlantic-Mediterranean transition a phylogeographical break?

Exploring reported population differences in Norway lobster (Nephrops norvegicus) in the Pomo Pits region of the Adriatic Sea using genome-wide markers

The Norway lobster (Nephrops norvegicus) is one of the most important decapod crustacean seafood species in the Adriatic Sea. Previous research has identified significant differences in growth rates and maturation timing of Nephrops in the Pomo/Jabuka Pits area compared to other subpopulations in Adriatic fishing grounds. Here, we use 1,623 genome-wide single nucleotide polymorphisms (SNPs) to investigate whether the Pomo Pits subpopulation is genetically different from other sites in the Adriatic and neighbouring seas. We found no genetic differentiation among all sampled Adriatic sites, suggesting high gene flow between Pomo Pits Nephrops and those of surrounding areas. We also found genetic homogeneity between the Adriatic sites and single-site samples from the Aegean and Tyrrhenian Seas. However, we detected distinct genetic differentiation between all Mediterranean sites and an Atlantic site in western Scotland, which provides evidence for a phylogenetic break between the Atlantic and the Mediterranean. Our results indicate that Pomo Pits Nephrops are not genetically different from others sampled in the Adriatic and that key biological parameters in Pomo Pits Nephrops could be driven by spatial variation in fishing pressure and/or environmental factors rather than geographic isolation.

Monitoring the genetic diversity in holothurian populations from the Italian coasts with mitochondrial DNA sequences

Sea cucumbers are deposit feeding members of marine benthic communities. The over-exploitation of sea cucumber natural stocks, especially in the Mediterranean basin, is having negative impacts on the marine ecosystem. This concerns had led the Italian government to take legal actions to preserve these important marine resources. The aim of this study was to evaluate the level of genetic diversity and population structures within two sympatric Holothuria sea cucumber species, H. polii and H. tubulosa, across ten Mediterranean areas along Italy. A 323 bp portion of the COI mitochondrial gene was sequenced in a total of 441 holothurians (251 H. polii, 177 H. tubulosa and 13 from other species). Genetic diversity analyses and analysis of molecular variance (AMOVA) suggest that H. polii and H. tubulosa species are distinct, but within and among populations are homogeneous, indicating active gene flows across the Mediterranean areas investigated. H. polii showed a lower genetic diversity than H. tubulosa, probably related to differences in life history traits. Phylogenetic analyses showed a clear differentiation between the two species, even if six specimens morphologically assigned to a species clustered within the other species, indicating a possible occurrence of hybridization events. These data will be useful in implementing conservation actions for these holothurian genetic resources.

Seascape Genetics and Distinct Intraspecific Diversification of the Decapod Nephrops norvegicus in the Adriatic Sea

Norway lobster Nephrops norvegicus, a prized decapod crustacean species, is found at different depths across the East Atlantic Ocean and Mediterranean Sea. Despite management efforts, the stocks are globally characterised as overexploited. In the present study, the impact of biogeographical boundaries on the phylogeographical and demographic population status was investigated within the Adriatic Sea, addressing important genetic indices for decapod functional conservation management. Central Mediterranean, Adriatic Sea A total of 482 individuals of Nephrops divided into the 12 samples were collected across biogeographical range of the Adriatic Sea. Using the mtDNA D-loop and microsatellite markers, methods of phylogeography and seascape genetics were applied to infer offshore versus coastal population divergence, demography and structure. Significant findings include genetic differentiation between offshore and coastal samples, with higher diversity indices in open waters. The limited gene flow observed between these two areas emphasises the self-sustained nature of coastal populations. Recent demographic changes in coastal populations reflect geographical constraints, fishing pressures and fluctuations in self-recruitment success. Additionally, the study reveals historical biogeographic events shaping the Adriatic populations, with evidence suggesting lineage divergence during the upper Pleistocene and postglacial recolonisation from southern Adriatic refugia. The role of biogeographical conditions in shaping genetic structure and limited gene flow between inshore and offshore areas underscore the need for improved management strategies, emphasising the importance of marine protected areas in conserving coastal populations and maintaining overall genetic diversity of the Norway lobster in the Adriatic Sea. Genomic monitoring within current management practices is recommended.

New insights into population structure, demographic history, and effective population size of the critically endangered blue shark Prionace glauca in the Mediterranean Sea

The blue shark, Prionace glauca, is the most abundant pelagic shark in the open ocean but its vulnerability remains poorly understood while being one of the most fecund sharks. In the Mediterranean Sea, the blue shark is listed as Critically Endangered (CR) by the International Union for Conservation of Nature. The species is facing a strong decline due to fishing, and scientific data regarding its genetic structure and vulnerability are still lacking. Here, we investigated the genetic diversity, demographic history, and population structure of the blue shark within the Mediterranean Sea, from samples of the Gulf of Lion and Malta, using sequences of the mtDNA control region and 22 microsatellite markers. We also compared our mitochondrial data to previous studies to examine the Atlantic-Mediterranean population structure. We assessed the blue shark's genetic vulnerability in the Mediterranean basin by modelling its effective population size. Our results showed a genetic differentiation between the Atlantic and the Mediterranean basins, with limited gene flow between the two areas, and distinct demographic histories making the Mediterranean population an independent management unit. Within the Mediterranean Sea, no sign of population structure was detected, suggesting a single population across the Western and Central parts of the sea. The estimated effective population size was low and highlighted the high vulnerability of the Mediterranean blue shark population, as the estimated size we calculated might not be sufficient to ensure the long-term persistence of the population. Our data also provide additional evidence that the Gulf of Lion area acts as a nursery for P. glauca, where protection is essential for the conservation strategy of the species in the Mediterranean.

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.

Santamaría MT, Faria R, Machado AM, Fonseca MM, Froufe E, C. Castro LF. Population Genomics Reveals the Underlying Structure of the Small Pelagic European Sardine and Suggests Low Connectivity within Macaronesia

The European sardine (Sardina pilchardus, Walbaum 1792) is indisputably a commercially important species. Previous studies using uneven sampling or a limited number of makers have presented sometimes conflicting evidence of the genetic structure of S. pilchardus populations. Here, we show that whole genome data from 108 individuals from 16 sampling areas across 5000 km of the species' distribution range (from the Eastern Mediterranean to the archipelago of Azores) support at least three genetic clusters. One includes individuals from Azores and Madeira, with evidence of substructure separating these two archipelagos in the Atlantic. Another cluster broadly corresponds to the center of the distribution, including the sampling sites around Iberia, separated by the Almeria-Oran front from the third cluster that includes all of the Mediterranean samples, except those from the Alboran Sea. Individuals from the Canary Islands appear to belong to the Mediterranean cluster. This suggests at least two important geographical barriers to gene flow, even though these do not seem complete, with many individuals from around Iberia and the Mediterranean showing some patterns compatible with admixture with other genetic clusters. Genomic regions corresponding to the top outliers of genetic differentiation are located in areas of low recombination indicative that genetic architecture also has a role in shaping population structure. These regions include genes related to otolith formation, a calcium carbonate structure in the inner ear previously used to distinguish S. pilchardus populations. Our results provide a baseline for further characterization of physical and genetic barriers that divide European sardine populations, and information for transnational stock management of this highly exploited species towards sustainable fisheries.

Oceanic currents maintain the genetic structure of non-marine coastal taxa in the western Mediterranean Sea

Coastal habitats are amongst the most dynamic on Earth, due to their simultaneous exposure to terrestrial, oceanic and atmospheric processes. Coastal taxa are therefore often ecologically specialised and adapted to withstand frequent shifts in sea level, wave exposure, temperature or salinity. This specialisation often resulting in significant cryptic diversity. Previous molecular studies have suggested that genetic differentiation in non-marine coastal organisms may be influenced by oceanic currents and fronts, but the extent to which such processes affect dispersal and evolution of such taxa remains unclear. Here we explore whether population genetic structure in two supralittoral rockpool beetle species (genus Ochthebius) can be predicted from the general circulation pattern of the marine currents and associated oceanic fronts. We simulated dispersal using a Lagrangian particle tracking model and compared this with population genetic structure inferred from COI (mitochondrial) and wingless (nuclear) genes applying linear models and Mantel tests. We show that a biophysical model based on oceanic currents and fronts in the western Mediterranean Sea is a much better predictor of observed population genetic structure than isolation by distance in both species. Our results show that oceanic processes, besides shaping contemporary population connectivity in fully marine organisms, also exert a meaningful influence on terrestrially-derived coastal taxa such as supralittoral rockpool beetles - the first time this mode of dispersal has been demonstrated in an insect.

Mediterranean fish communities are struggling to adapt to global warming. Evidence from the western coast of Italy

Climate change has significant impacts on marine ecosystems, resulting in disruptions in biological interactions, shifts in community composition, and changes in the physiology of fish and other marine organisms. In this study conducted in the central Mediterranean Sea, the mean temperature of the catch (MTC) was employed as an indicator to investigate the climatological factors influencing the fish community. The MTC, which utilizes species-preferred temperatures, was calculated using bottom temperature (BT) data weighted against scientific catches. The estimated MTC increasing rates were 0.01 °C year-1 for the entire community, 0.017 °C year-1 for the shelf break, and 0.004 °C year-1 for the continental slope assemblage. We found that MTC is increasing at a lower rate compared to BT, suggesting a progressive under-adaptation of the fish community that seems not fully able to keep up with the ongoing pace of warming. The study identified sea surface temperature and bottom temperature as key drivers of changes in fish community composition. Notably, the fish community composition exhibited drastic changes over the studied period, and we suggest that the MTC can be a useful index to monitor such changes within the context of the EU's climate change adaptation strategy.

Insights on the Evolutionary History and Genetic Patterns of Octopus vulgaris Cuvier, 1797 in the Northeastern Atlantic Using Mitochondrial DNA

Octopus vulgaris is one of the most harvested octopus species in the world. In the Iberian Peninsula, there are several small-scale fisheries that have a long-term tradition of harvesting octopus. The Asturias fleet (in Northern Spain) has an internationally recognized MSC label for its exploitation. Of concern, genetic assessments of exploited stocks are currently scarce, which could prevent the implementation of adequate managing strategies. We use two mitochondrial regions (cytochrome oxidase subunit 1 and control region) to analyze the genetic status and evolutionary events that conditioned octopus populations' characteristics in the Northeastern Atlantic. A total of 90 individuals were sampled from three different localities in the Iberian Peninsula as well as a location in Macaronesia. Temporal genetic analyses on Asturias and Algarve populations were also performed. Results indicated the absence of fine spatial genetic structuring but showed the Canary Islands (in Macaronesia) as the most distinct population. Our analyses detected two distinct clades, already described in the literature, but, for the first time, we confirmed the presence of the α-southern haplogroup in the Northern Iberian Peninsula. This result indicates a more continuous cline for the distribution of these two haplogroups than previously reported. Temporal changes in the distribution of both haplogroups in contact zones were also detected.

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.

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.

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.

Heterogeneous Genomic Divergence Landscape in Two Commercially Important European Scallop Species

Two commercially important scallop species of the genus Pecten are found in Europe: the north Atlantic Pecten maximus and the Mediterranean Pecten jacobaeus whose distributions abut at the Almeria-Orán front. Whilst previous studies have quantified genetic divergence between these species, the pattern of differentiation along the Pecten genome is unknown. Here, we mapped RADseq data from 235 P. maximus and 27 P. jacobaeus to a chromosome-level reference genome, finding a heterogeneous landscape of genomic differentiation. Highly divergent genomic regions were identified across 14 chromosomes, while the remaining five showed little differentiation. Demographic and comparative genomics analyses suggest that this pattern resulted from an initial extended period of isolation, which promoted divergence, followed by differential gene flow across the genome during secondary contact. Single nucleotide polymorphisms present within highly divergent genomic regions were located in areas of low recombination and contrasting patterns of LD decay were found between the two species, hinting at the presence of chromosomal inversions in P. jacobaeus. Functional annotations revealed that highly differentiated regions were enriched for immune-related processes and mRNA modification. While future work is necessary to characterize structural differences, this study provides new insights into the speciation genomics of P. maximus and P. jacobaeus.

Three mitochondrial lineages and no Atlantic-Mediterranean barrier for the bogue Boops boops across its widespread distribution

Marine species exhibiting wide distributional ranges are frequently subdivided into discrete genetic units over limited spatial scales. This is often due to specific life-history traits or oceanographic barriers that prevent gene flow. Fine-scale sampling studies revealed distinct phylogeographic patterns in the northeastern Atlantic and the Mediterranean, ranging from panmixia to noticeable population genetic structure. Here, we used mitochondrial sequence data to analyse connectivity in the bogue Boops boops throughout most of its widespread distribution. Our results identified the existence of three clades, one comprising specimens from the Azores and eastern Atlantic/Mediterranean, another with individuals from the Canary Islands, Madeira and Cape Verde archipelagos, and the third with samples from Mauritania only. One of the branches of the northern subtropical gyre (Azores Current) that drifts towards the Gulf of Cádiz promotes a closer connection between the Azores, southern Portugal and the Mediterranean B. boops populations. The Almería-Oran Front, widely recognised as an oceanographic barrier for many organisms to cross the Atlantic-Mediterranean divide, does not seem to affect the dispersal of this benthopelagic species. The southward movement of the Cape Verde Frontal Zone during the winter, combined with the relatively short duration of the pelagic larval stage of B. boops, may be potential factors for preventing the connectivity between the Atlantic oceanic archipelagos and Mauritania shaping the genetic signature of this species.

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.

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.

Disentangling the Taxonomic Status of Caprella penantis sensu stricto (Amphipoda: Caprellidae) Using an Integrative Approach

Despite its importance in intertidal and shallow-water marine ecosystems, Caprella penantis continues to be one of the most taxonomically challenging amphipods in the world. A recent molecular study focusing on C. penantis sensu stricto pointed out the existence of three highly divergent lineages, indicating the possible existence of a process of ongoing speciation and, thus, casting doubt on the taxonomic status of this species. In the present study, we used an integrative approach to continue to shed light on the taxonomy and distribution of this caprellid. To this end, we combined morphological and genetic data (COI and 18S) and included, for the first time, populations from its type locality. Our analyses provide strong evidence of the existence of potentially three distinct species, genetically and geographically restricted, within C. penantis sensu stricto, with the distribution of the true C. penantis sensu stricto restricted to the UK (type locality), the northern coast of the Iberian Peninsula, and the Azores. Results show the co-occurrence of two of these species in a locality of northern Portugal and indicate the existence of distinct evolutionary and diversification patterns along the eastern Atlantic region. Overall, our study highlights the use of an integrative approach to properly assess species boundaries and unravel hidden biodiversity in amphipods.

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.

Natural experiments and long-term monitoring are critical to understand and predict marine host-microbe ecology and evolution

Marine multicellular organisms host a diverse collection of bacteria, archaea, microbial eukaryotes, and viruses that form their microbiome. Such host-associated microbes can significantly influence the host’s physiological capacities; however, the identity and functional role(s) of key members of the microbiome (“core microbiome”) in most marine hosts coexisting in natural settings remain obscure. Also unclear is how dynamic interactions between hosts and the immense standing pool of microbial genetic variation will affect marine ecosystems’ capacity to adjust to environmental changes. Here, we argue that significantly advancing our understanding of how host-associated microbes shape marine hosts’ plastic and adaptive responses to environmental change requires (i) recognizing that individual host–microbe systems do not exist in an ecological or evolutionary vacuum and (ii) expanding the field toward long-term, multidisciplinary research on entire communities of hosts and microbes. Natural experiments, such as time-calibrated geological events associated with well-characterized environmental gradients, provide unique ecological and evolutionary contexts to address this challenge. We focus here particularly on mutualistic interactions between hosts and microbes, but note that many of the same lessons and approaches would apply to other types of interactions.

This Essay argues that in order to truly understand how marine hosts benefit from the immense diversity of microbes, we need to expand towards long-term, multi-disciplinary research focussing on few areas of the world’s ocean that we refer to as “natural experiments,” where processes can be studied at scales that far exceed those captured in laboratory experiments.

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.

Comparative genetics of scyphozoan species reveals the geological history and contemporary processes of the Mediterranean Sea

Jellyfish are useful genetic indicators for aquatic ecosystems as they have limited mobility and are highly exposed to the water column. By using comparative genomics and the molecular clock (timetree) of Rhizostoma pulmo, we revealed a divergence point between the East and West Mediterranean Sea (MS) populations that occurred 4.59 million years ago (mya). It is suggested that the two distinct ecological environments we know today were formed at this time. We propose that before this divergence, the highly mixed Atlantic and Mediterranean waters led to the wide dispersal of different species including R. pulmo. At 4.59 mya, the Western and Eastern MS were formed, indicating the possibility of a dramatic environmental event. For the first time, we find that for the jellyfish we examined, the division of the MS in east and west is not at the Straits of Sicily as generally thought, but significantly to the east. Using genomics of the Aurelia species, we examined contemporary anthropogenic impacts with a focus on migration of scyphozoa across the Suez Canal (Lessepsian migration). Aurelia sp. is among the few scyphozoa we find in both the MS and the Red Sea, but our DNA analysis revealed that the Red Sea Aurelia sp. did not migrate or mix with MS species. Phyllorhiza punctata results showed that this species was only recently introduced to the MS as a result of anthropogenic transportation activity, such as ballast water discharge, and revealed a migration vector from Australia to the MS. Our findings demonstrate that jellyfish genomes can be used as a phylogeographic molecular tool to trace past events across large temporal scales and reveal invasive species introduction due to human activity.

A star is torn-molecular analysis divides the Mediterranean population of Poli's stellate barnacle, Chthamalus stellatus (Cirripedia, Chtamalidae)

Poli's stellate barnacle, Chthamalus stellatus Poli, populates the Mediterranean Sea, the North-Eastern Atlantic coasts, and the offshore Eastern Atlantic islands. Previous studies have found apparent genetic differences between the Atlantic and the Mediterranean populations of C. stellatus, suggesting possible geological and oceanographic explanations for these differences. We have studied the genetic diversity of 14 populations spanning from the Eastern Atlantic to the Eastern Mediterranean, using two nuclear genes sequences revealing a total of 63 polymorphic sites. Both genotype-based, haplotype-based and the novel SNP distribution population-based methods have found that these populations represent a geographic cline along the west to east localities. The differences in SNP distribution among populations further separates a major western cluster into two smaller clusters, the Eastern Atlantic and the Western Mediterranean. It also separates the major eastern cluster into two smaller clusters, the Mid-Mediterranean and Eastern Mediterranean. We suggested here environmental conditions like surface currents, water salinity and temperature as probable factors that have formed the population structure. We demonstrate that C. stellatus is a suitable model organism for studying how geological events and hydrographic conditions shape the fauna in the Mediterranean Sea.

Towards omics-based predictions of planktonic functional composition from environmental data

Marine microbes play a crucial role in climate regulation, biogeochemical cycles, and trophic networks. Unprecedented amounts of data on planktonic communities were recently collected, sparking a need for innovative data-driven methodologies to quantify and predict their ecosystemic functions. We reanalyze 885 marine metagenome-assembled genomes through a network-based approach and detect 233,756 protein functional clusters, from which 15% are functionally unannotated. We investigate all clusters' distributions across the global ocean through machine learning, identifying biogeographical provinces as the best predictors of protein functional clusters' abundance. The abundances of 14,585 clusters are predictable from the environmental context, including 1347 functionally unannotated clusters. We analyze the biogeography of these 14,585 clusters, identifying the Mediterranean Sea as an outlier in terms of protein functional clusters composition. Applicable to any set of sequences, our approach constitutes a step towards quantitative predictions of functional composition from the environmental context.

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.

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.

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.

Molecular analysis reveals a cryptic species of Chthamalus (Crustacea: Cirripedia) in the Cape Verde Islands

The intertidal barnacle Chthamalus stellatus has a broad distribution, occurring in the Mediterranean, the east Atlantic shores and east Atlantic Macaronesian Islands (Madeira, the Canaries and the Azores). Traditionally, based on morphological characters, Chthamalus of the Cape Verde Islands were also regarded as C. stellatus. However, using a mitochondrial gene and two nuclear genes, we found that although Chthamalus from Cape Verde is morphologically similar to C. stellatus, there are genetic differences between the two that are larger than those found between different species of Chthamalus. We thus claim that these genetic differences justify the assignment of the Cape Verde populations as an evolutionarily significant unit and a sister clade to C. stellatus. We also show that the connection between taxonomic units that are close to each other lies not only in the resemblance between DNA sequences. We have found that numerous point mutations characterizing the Cape Verde Chthamalus are present as infrequent alleles in C. stellatus, indicating that two close taxonomic units can also share polymorphisms present in their common ancestor.

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.

Multiple Non-Species-Specific Pathogens Possibly Triggered the Mass Mortality in Pinna nobilis

The fan mussel, Pinna nobilis, represents the largest bivalve endemic to the Mediterranean Sea. Since 2016, dramatic mass mortality of this species has been observed in several areas. The first surveys suggested that Haplosporidium pinnae (currently considered species-specific) was the main etiological agent, but recent studies have indicated that a multifactorial disease may be responsible for this phenomenon. In this study, we performed molecular diagnostic analyses on P. nobilis, P. rudis, and bivalve heterologous host species from the island of Sardinia to shed further light on the pathogens involved in the mass mortality. The results support the occurrence of a multifactorial disease and that Mycobacterium spp. and H. pinnae are not necessarily associated with the illness. Indeed, our analyses revealed that H. pinnae is not species-specific for P. nobilis, as it was present in other bivalves at least three years before the mass mortality began, and species of Mycobacterium were also found in healthy individuals of P. nobilis and P. rudis. We also detected the species Rhodococcus erythropolis, representing the first report in fan mussels of a bacterium other than Mycobacterium spp. and Vibrio spp. These results depict a complicated scenario, further demonstrating how the P. nobilis mass mortality event is far from being fully understood.

Mitogenome phylogenetics in the genus Palaemon (Crustacea: Decapoda) sheds light on species crypticism in the rockpool shrimp P. elegans

The genus Palaemon comprises worldwide marine and freshwater shrimps and prawns, and some of them are ecologically or commercially important species. Palaemon is not currently a monophyletic group, so phylogenetics and systematics are constantly changing. Species crypticism has been pointed out in several Palaemon species, being the clearest evidence in the European rockpool shrimp P. elegans. Here we sequenced and described seven European Palaemon mitochondrial genomes. The mitochondrial protein-coding genes were used, along with those of three other Palaemon species, to perform mitogenome phylogenetic analyses to clarify the evolutionary relationships within the genus, and particularly to shed light on the cryptic species found within P. elegans. The Messinian Salinity Crisis (5.3-5.9 Ma, late Miocene) was proposed to be the origin of this cryptic species and it was used as aged constraint for calibration analysis. We provide the largest and the first time-calibrated mitogenome phylogeny of the genus Palaemon and mitogenome substitution rate was estimated (1.59% per million years) in Decapoda for the first time. Our results highlighted the need for future systematics changes in Palaemon and crypticism in P. elegans was confirmed. Mitochondrial genome and cox1 (1.41%) substitution rate estimates matched those published elsewhere, arguing that the Messinian Salinity Crisis was a plausible event driving the split between P. elegans and its cryptic species. Molecular dating suggested that Pleistocene glaciations were likely involved in the differentiation between the Atlantic and Mediterranean populations of P. elegans. On the contrary, the divergence between the Atlantic and Mediterranean populations of the common littoral shrimp P. serratus was greater and dated to be much older (4.5-12.3 Ma, Plio-Miocene), so we considered that they could represent two separated species. Therefore, species crypticism in the genus Palaemon seems to be a common phenomenon.

Individual-based population genomics reveal different drivers of adaptation in sympatric fish

Connectivity and local adaptation are two contrasting evolutionary forces highly influencing population structure. To evaluate the impact of early-life traits and environmental conditions on genetic structuring and adaptation, we studied two sympatric fish species in the Western Mediterranean Sea: Symphodus tinca and S. ocellatus. We followed an individual-based approach and measured early-life history traits from otolith readings, gathered information on environmental variables and obtained genome-wide markers from genotyping-by-sequencing (GBS). The two species presented contrasting population structure across the same geographic gradient, with high and significant population differentiation in S. ocellatus, mostly determined by oceanographic fronts, and low differentiation and no front effect in S. tinca. Despite their different levels of genetic differentiation, we identified in both species candidate regions for local adaptation by combining outlier analysis with environmental and phenotypic association analyses. Most candidate loci were associated to temperature and productivity in S. ocellatus and to temperature and turbulence in S. tinca suggesting that different drivers may determine genomic diversity and differentiation in each species. Globally, our study highlights that individual-based approach combining genomic, environmental and phenotypic information is key to identify signals of selection and the processes mediating them.

Population genetic structure of the short-beaked common dolphin from the Black Sea and the Turkish Straits System

Our study aims to assess the population connectivity, evolutionary history, and conservation status of the short-beaked common dolphin in the Black Sea and Turkish Straits System (TSS). We also include DNA sequences from the Atlantic Ocean and the Mediterranean Sea to provide a regional perspective to our localized study. Analysis of 366 base pairs of mitochondrial DNA D-loop fragments from 37 samples collected from short-beaked common dolphins in the Black Sea, TSS, and Aegean Sea revealed 13 haplotypes, eight of which have not been previously reported. While analysis of samples archived on GenBank revealed 89 different haplotypes across the region. The haplotype network contains two main peripheral groups that include individuals from all locations. Haplotypes from the Atlantic Ocean are scattered across the network and no obvious population separation was detected. Some shared haplotypes potentially indicate multi-directional colonization events of the Mediterranean Sea from the eastern Atlantic Ocean. Moreover, some less widely distributed haplotypes suggest some level of more recent genetic connectivity through the Strait of Gibraltar and the TSS and point out the importance of these straits in the dispersal of short-beaked common dolphins. The haplotype and nucleotide diversity values were lower in the Black Sea, TSS, and western Mediterranean Sea when compared to the Atlantic Ocean, supporting the expansion of Atlantic populations into the Mediterranean and the Black Seas. Differentiation was observed between the Atlantic Ocean, and the Mediterranean Sea, TSS and the Black Sea based on Фst but not between Mediterranean and the Black Seas. For common dolphins, which have high dispersal potential, the protection of open seas and narrow seaways to enhance connectivity may be crucial.

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.

Deep phylogeographic structure may indicate cryptic species within the Sparid genus Spondyliosoma

Two geographically nonoverlapping species are currently described within the sparid genus Spondyliosoma: Spondyliosoma cantharus (Black Seabream) occurring across Mediterranean and eastern Atlantic waters from NW Europe to Angola and S. emarginatum (Steentjie) considered endemic to southern Africa. To address prominent knowledge gaps this study investigated range-wide phylogeographic structure across both species. Mitochondrial DNA sequences revealed deep phylogeographic structuring with four regionally partitioned reciprocally monophyletic clades, a Mediterranean clade and three more closely related Atlantic clades [NE Atlantic, Angola and South Africa (corresponding to S. emarginatum)]. Divergence and distribution of the lineages reflects survival in, and expansion from, disjunct glacial refuge areas. Cytonuclear differentiation of S. emarginatum supports its validity as a distinct species endemic to South African waters. However, the results also indicate that S. cantharus may be a cryptic species complex wherein the various regional lineages represent established/incipient species. A robust multilocus genetic assessment combining morphological data and detailing interactions among lineages is needed to determine the full diversity within Spondyliosoma and the most adequate biological and taxonomic status.

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.

Highly regional population structure of Spondyliosoma cantharus depicted by nuclear and mitochondrial DNA data

Resolution of population structure represents an effective way to define biological stocks and inform efficient fisheries management. In the present study, the phylogeography of the protogynous sparid Spondyliosoma cantharus, in the East Atlantic and Mediterranean Sea, was investigated with nuclear (S7) and mitochondrial (cytochrome b) DNA markers. Significant divergence of four regional genetic groups was observed: North Eastern Atlantic, Mediterranean Sea, Western African Transition (Cape Verde) and Gulf of Guinea (Angola). The two southern populations (Cape Verde and Angola) each comprised reciprocally monophyletic mtDNA lineages, revealed low levels of diversity in Cape Verde and high diversity for Angola despite being represented by only 14 individuals. A complete divergence between North Atlantic and Mediterranean populations was depicted by the mitochondrial marker, but a highly shared nuclear haplotype revealed an incomplete lineage sorting between these regions. Bayesian skyline plots and associated statistics revealed different dynamics among the four regions. Cape Verde showed no expansion and the expansion time estimated for Angola was much older than for the other regions. Mediterranean region seems to have experienced an early population growth but has remained with a stable population size for the last 30000 years while the North Atlantic population has been steadily growing. The lack of genetic structuring within these regions should not be taken as evidence of demographic panmixia in light of potential resolution thresholds and previous evidence of intra-regional phenotypic heterogeneity.

Speciation history of European (Anguilla anguilla) and American eel (A. rostrata), analysed using genomic data

Speciation in the ocean could differ from terrestrial environments due to fewer barriers to gene flow. Hence, sympatric speciation might be common, with American and European eel being candidates for exemplifying this. They show disjunct continental distributions on both sides of the Atlantic, but spawn in overlapping regions of the Sargasso Sea from where juveniles are advected to North American, European and North African coasts. Hybridization and introgression are known to occur, with hybrids almost exclusively observed in Iceland. Different speciation scenarios have been suggested, involving either vicariance or sympatric ecological speciation. Using RAD sequencing and whole-genome sequencing data from parental species and F1 hybrids, we analysed speciation history based on the joint allele frequency spectrum (JAFS) and pairwise sequentially Markovian coalescent (PSMC) plots. JAFS supported a model involving a split without gene flow 150,000-160,000 generations ago, followed by secondary contact 87,000-92,000 generations ago, with 64% of the genome experiencing restricted gene flow. This supports vicariance rather than sympatric speciation, likely associated with Pleistocene glaciation cycles and ocean current changes. Whole-genome PSMC analysis of F1 hybrids from Iceland suggested divergence 200,000 generations ago and indicated subsequent gene flow rather than strict isolation. Finally, simulations showed that results from both approaches (JAFS and PSMC) were congruent. Hence, there is strong evidence against sympatric speciation in North Atlantic eels. These results reiterate the need for careful consideration of cases of possible sympatric speciation, as even in seemingly barrier-free oceanic environments palaeoceanographic factors may have promoted vicariance and allopatric speciation.