Isolation and characterization of eight polymorphic microsatellite loci for the Mediterranean gorgonian Paramuricea clavata

The gorgonian Paramuricea clavata is a benthic organism often included in conservation management plans since it creates complexity in the ecosystems and is extremely vulnerable to disturbances. Eight microsatellite markers isolated from an enriched genomic library were characterized in a total of 50 individuals from two north-western Mediterranean populations. All loci were polymorphic and the number of alleles per locus ranged between 2 and 14. Significant genetic differentiation was observed between populations. The polymorphic markers presently isolated will allow for the assessment of the spatial genetic structure between and within populations of this umbrella species of the Mediterranean Sea.     

Local adaptation and oceanographic connectivity patterns explain genetic differentiation of a marine diatom across the North Sea–Baltic Sea salinity gradient

Drivers of population genetic structure are still poorly understood in marine micro‐organisms. We exploited the North Sea–Baltic Sea transition for investigating the seascape genetics of a marine diatom, Skeletonema marinoi. Eight polymorphic microsatellite loci were analysed in 354 individuals from ten locations to analyse population structure of the species along a 1500‐km‐long salinity gradient ranging from 3 to 30 psu. To test for salinity adaptation, salinity reaction norms were determined for sets of strains originating from three different salinity regimes of the gradient. Modelled oceanographic connectivity was compared to directional relative migration by correlation analyses to examine oceanographic drivers. Population genetic analyses showed distinct genetic divergence of a low‐salinity Baltic Sea population and a high‐salinity North Sea population, coinciding with the most evident physical dispersal barrier in the area, the Danish Straits. Baltic Sea populations displayed reduced genetic diversity compared to North Sea populations. Growth optima of low salinity isolates were significantly lower than those of strains from higher native salinities, indicating local salinity adaptation. Although the North Sea–Baltic Sea transition was identified as a barrier to gene flow, migration between Baltic Sea and North Sea populations occurred. However, the presence of differentiated neutral markers on each side of the transition zone suggests that migrants are maladapted. It is concluded that local salinity adaptation, supported by oceanographic connectivity patterns creating an asymmetric migration pattern between the Baltic Sea and the North Sea, determines genetic differentiation patterns in the transition zone.     

Phylogeography of the Atlanto-Mediterranean sea cucumber Holothuria (Holothuria) mammata: the combined effects of historical processes and current oceanographical pattern

We assessed the genetic structure of populations of the widely distributed sea cucumber Holothuria (Holothuria) mammata Grube, 1840, and investigated the effects of marine barriers to gene flow and historical processes. Several potential genetic breaks were considered, which would separate the Atlantic and Mediterranean basins, the isolated Macaronesian Islands from the other locations analysed, and the Western Mediterranean and Aegean Sea (Eastern Mediterranean). We analysed mitochondrial 16S and COI gene sequences from 177 individuals from four Atlantic locations and four Mediterranean locations. Haplotype diversity was high (H=0.9307 for 16S and 0.9203 for COI), and the haplotypes were closely related (π=0.0058 for 16S and 0.0071 for COI). The lowest genetic diversities were found in the Aegean Sea population. Our results showed that the COI gene was more variable and more useful for the detection of population structure than the 16S gene. The distribution of mtDNA haplotypes, the pairwise F(ST) values and the results of exact tests and amova revealed: (i) a significant genetic break between the population in the Aegean Sea and those in the other locations, as supported by both mitochondrial genes, and (ii) weak differentiation of the Canary and Azores Islands from the other populations; however, the populations from the Macaronesian Islands, Algarve and West Mediterranean could be considered to be a panmictic metapopulation. Isolation by distance was not identified in H. (H.) mammata. Historical events behind the observed findings, together with the current oceanographic patterns, were proposed and discussed as the main factors that determine the population structure and genetic signature of H. (H.) mammata.     

High genetic differentiation of red gorgonian populations from the Atlantic Ocean and the Mediterranean Sea

Patterns of genetic variation within a species may be used to infer past events in the evolutionary history of marine species. In the present study we aimed to compare the genetic diversity of the red gorgonian Paramuricea clavata in the Atlantic Ocean and the Mediterranean Sea. For genetic markers we used microsatellites and a mitochondrial gene fragment. Our results revealed a distinct genetic composition and diversity between the Mediterranean and the Atlantic. The Mediterranean samples had higher microsatellite heterozygosity, allelic richness and private allelic richness. The hypotheses that can explain these patterns are the isolation of Atlantic populations and/or a founder effect. Additionally, a clear difference was obtained from the mitochondrial locus, since sequences from Atlantic and Mediterranean samples diverged by 1%, which is high for soft corals.     

Phylogeography of the marbled crab Pachygrapsus marmoratus (Decapoda,Grapsidae) along part of the African Mediterranean coast reveals genetic homogeneity across the Siculo-Tunisian Strait versus heterogeneity across the Gibraltar Strait

We investigate the influence of previously postulated biogeographic barriers in the Mediterranean Sea on the population genetic structure of a highly dispersive and continuously distributed coastal species. In particular, we examine nuclear and mitochondrial genetic variation in the marbled crab, Pachygrapsus marmoratus, across part of the African Mediterranean coast in order to assess the influence of the Siculo-Tunisian Strait on its population genetic structure. Four polymorphic microsatellite loci were genotyped for 110 individuals, collected from eight locations covering parts of the Algerian, Tunisian and Libyan coasts. In addition, mtDNA corresponding to the Cox1 gene was sequenced for 80 samples. The corresponding results show contrasting patterns of genetic differentiation. While mtDNA results revealed a homogeneous haplotype composition in our study area, microsatellite data depicted genetic differentiation among populations, but not associated with any geographic barrier. This pattern, already recorded for this species from different geographic regions, may hint at the involvement of a complex series of abiotic and biotic factors in determining genetic structure. Demographic history reconstruction, inferred from mtDNA data, supports demographic and spatial expansion for the North African metapopulation dating back to the Mid-Pleistocene and following an historical bottleneck. Comparison of these African mitochondrial sequences with new sequences from a Turkish population and previously published sequences revealed a weak but significant separation of Atlantic and Mediterranean populations across the Gibraltar Strait, which was not recorded in previous studies of this grapsid species.     

Marked mitochondrial DNA differences between Mediterranean and Atlantic populations of the swordfish, Xiphias gladius

Restriction analysis of mitochondrial DNA (mtDNA) from 204 individuals of swordfish (Xiphias gladius) revealed no differentiation among samples from three sites in the Mediterranean Sea (Greece, Italy, Spain), but a high degree of differentiation between Mediterranean samples and a sample from the Gulf of Guinea. A fifth sample from the Atlantic side of the Straits of Gibraltar (Tarifa) consisted mostly of mitotypes that are common in the Mediterranean, but contained several of mtDNA types of the Guinea sample not found in the Mediterranean. We conclude that, in spite of free migration of swordfish across the Straits of Gibraltar, little genetic exchange occurs between the populations inhabiting the Mediterranean Sea and the tropical Atlantic ocean. This is the first evidence of genetic differentiation among geographic populations of this highly mobile species that supports a world-wide fishery.     

Genetic biodiversity in the Baltic Sea: species-specific patterns challenge management

Information on spatial and temporal patterns of genetic diversity is a prerequisite to understanding the demography of populations, and is fundamental to successful management and conservation of species. In the sea, it has been observed that oceanographic and other physical forces can constitute barriers to gene flow that may result in similar population genetic structures in different species. Such similarities among species would greatly simplify management of genetic biodiversity. Here, we tested for shared genetic patterns in a complex marine area, the Baltic Sea. We assessed spatial patterns of intraspecific genetic diversity and differentiation in seven ecologically important species of the Baltic ecosystem—Atlantic herring (Clupea harengus), northern pike (Esox lucius), European whitefish (Coregonus lavaretus), three-spined stickleback (Gasterosteus aculeatus), nine-spined stickleback (Pungitius pungitius), blue mussel (Mytilus spp.), and bladderwrack (Fucus vesiculosus). We used nuclear genetic data of putatively neutral microsatellite and SNP loci from samples collected from seven regions throughout the Baltic Sea, and reference samples from North Atlantic areas. Overall, patterns of genetic diversity and differentiation among sampling regions were unique for each species, although all six species with Atlantic samples indicated strong resistence to Atlantic-Baltic gene-flow. Major genetic barriers were not shared among species within the Baltic Sea; most species show genetic heterogeneity, but significant isolation by distance was only detected in pike and whitefish. These species-specific patterns of genetic structure preclude generalizations and emphasize the need to undertake genetic surveys for species separately, and to design management plans taking into consideration the specific structures of each species.     

Matching genetics with oceanography: directional gene flow in a Mediterranean fish species

Genetic connectivity and geographic fragmentation are two opposing mechanisms determining the population structure of species. While the first homogenizes the genetic background across populations the second one allows their differentiation. Therefore, knowledge of processes affecting dispersal of marine organisms is crucial to understand their genetic distribution patterns and for the effective management of their populations. In this study, we use genetic analyses of eleven microsatellites in combination with oceanographic satellite and dispersal simulation data to determine distribution patterns for Serranus cabrilla, a ubiquitous demersal broadcast spawner, in the Mediterranean Sea. Pairwise population F_ST values ranged between ?0.003 and 0.135. Two genetically distinct clusters were identified, with a clear division located between the oceanographic discontinuities at the Ibiza Channel (IC) and the Almeria‐Oran Front (AOF), revealing an admixed population in between. The Balearic Front (BF) also appeared to dictate population structure. Directional gene flow on the Spanish coast was observed as S. cabrilla dispersed from west to east over the AOF, from north to south on the IC and from south of the IC towards the Balearic Islands. Correlations between genetic and oceanographic data were highly significant. Seasonal changes in current patterns and the relationship between ocean circulation patterns and spawning season may also play an important role in population structure around oceanographic fronts.     

Genetic Structuring across Marine Biogeographic Boundaries in Rocky Shore Invertebrates

Biogeography investigates spatial patterns of species distribution. Discontinuities in species distribution are identified as boundaries between biogeographic areas. Do these boundaries affect genetic connectivity? To address this question, a multifactorial hierarchical sampling design, across three of the major marine biogeographic boundaries in the central Mediterranean Sea (Ligurian-Tyrrhenian, Tyrrhenian-Ionian and Ionian-Adriatic) was carried out. Mitochondrial COI sequence polymorphism of seven species of Mediterranean benthic invertebrates was analysed. Two species showed significant genetic structure across the Tyrrhenian-Ionian boundary, as well as two other species across the Ionian Sea, a previously unknown phylogeographic barrier. The hypothesized barrier in the Ligurian-Tyrrhenian cannot be detected in the genetic structure of the investigated species. Connectivity patterns across species at distances up to 800 km apart confirmed that estimates of pelagic larval dispersal were poor predictors of the genetic structure. The detected genetic discontinuities seem more related to the effect of past historical events, though maintained by present day oceanographic processes. Multivariate statistical tools were used to test the consistency of the patterns across species, providing a conceptual framework for across-species barrier locations and strengths. Additional sequences retrieved from public databases supported our findings. Heterogeneity of phylogeographic patterns shown by the 7 investigated species is relevant to the understanding of the genetic diversity, and carry implications for conservation biology.     

Population differentiation in the open sea: insights from the pelagic copepod Pleuromamma xiphias

Although a number of recent studies of marine holoplankton have reported significant genetic structure among populations, little is currently known about the biological and oceanographic processes that influence population connectivity in oceanic plankton. In order to examine how depth preferences influence dispersal in oceanic plankton, I characterized the genetic structure of a copepod with diel vertical migration (DVM) (Pleuromamma xiphias), throughout its global distribution, and compared these results to those expected given the interaction of this species' habitat depth with ocean circulation and bathymetry. Mitochondrial COI sequences from 651 individuals from 28 sites in the Indian, Pacific, and Atlantic Oceans revealed highly significant genetic differentiation both within and among ocean basins. Limited dispersal among distinct pelagic provinces seems to have played a major role in population differentiation in this species, with strong genetic breaks observed across known oceanographic fronts or current systems in all three ocean basins. The Indo-West Pacific (IWP) holds a highly distinct genetic population of this species that was sampled in both the western Pacific and eastern Indian Oceans. This suggests that the IWP does not act as a strong barrier to gene flow between basins, as expected, despite the relatively shallow water depth (<200 m) and vertically extensive (>400 m) diel migration of this species. A pattern of isolation by distance was observed in the Indian Ocean with genetic differentiation among samples down to spatial scales of ～800 km, indicating that realized dispersal in P. xiphias occurs over much smaller spatial scales than in previously reported oceanic holoplankton. Given its highly regionalized population genetic structure, P. xiphias may have some capacity to adapt to local oceanographic conditions, and it should not be assumed that populations of this species in distinct pelagic biomes will respond in the same way to shared physical or climatic forcing.     

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

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

Genetic population structure of the endemic fourline wrasse (Larabicus quadrilineatus) suggests limited larval dispersal distances in the Red Sea

The connectivity among marine populations is determined by the dispersal capabilities of adults as well as their eggs and larvae. Dispersal distances and directions have a profound effect on gene flow and genetic differentiation within species. Genetic homogeneity over large areas is a common feature of coral reef fishes and can reflect high dispersal capability resulting in high levels of gene flow. If fish larvae return to their parental reef, gene flow would be restricted and genetic differentiation could occur. Larabicus quadrilineatus (Labridae) is considered as an endemic fish species of the Red Sea and Gulf of Aden. The juveniles of this species are cleaner fish that feed on ectoparasites of other fishes. Here, we investigated the genetic population structure and gene flow in L. quadrilineatus among five locations in the Red Sea to infer connectivity among them. To estimate genetic diversity, we analysed 369 bp of 237 mitochondrial DNA control region sequences. Haplotype and nucleotide diversities were higher in the southern than in the northern Red Sea. Analysis of molecular variance (amova) detected the highest significant genetic variation between northern and central/southern populations (Phi(CT) = 0.01; P < 0.001). Migration analysis revealed a several fold higher northward than southward migration, which could be explained by oceanographic conditions and spawning season. Even though the Phi(ST) value of 0.01 is rather low and implies a long larval dispersal distance, estimates based on the isolation-by-distance model show a very low mean larval dispersal distance (0.44-5.1 km) compared to other studies. In order to enable a sustainable ornamental fishery on the fourline wrasse, the results of this study suggest that populations in the northern and southern Red Sea should be managed separately as two different stocks. The rather low larval dispersal distance of about 5 km needs to be considered in the design of marine protected areas to enable connectivity and self-seeding.     

Habitat continuity and stepping‐stone oceanographic distances explain population genetic connectivity of the brown alga Cystoseira amentacea

Effective predictive and management approaches for species occurring in a metapopulation structure require good understanding of interpopulation connectivity. In this study, we ask whether population genetic structure of marine species with fragmented distributions can be predicted by stepping‐stone oceanographic transport and habitat continuity, using as model an ecosystem‐structuring brown alga, Cystoseira amentacea var. stricta. To answer this question, we analysed the genetic structure and estimated the connectivity of populations along discontinuous rocky habitat patches in southern Italy, using microsatellite markers at multiple scales. In addition, we modelled the effect of rocky habitat continuity and ocean circulation on gene flow by simulating Lagrangian particle dispersal based on ocean surface currents allowing multigenerational stepping‐stone dynamics. Populations were highly differentiated, at scales from few metres up to thousands of kilometres. The best possible model fit to explain the genetic results combined current direction, rocky habitat extension and distance along the coast among rocky sites. We conclude that a combination of variable suitable habitat and oceanographic transport is a useful predictor of genetic structure. This relationship provides insight into the mechanisms of dispersal and the role of life‐history traits. Our results highlight the importance of spatially explicit modelling of stepping‐stone dynamics and oceanographic directional transport coupled with habitat suitability, to better describe and predict marine population structure and differentiation. This study also suggests the appropriate spatial scales for the conservation, restoration and management of species that are increasingly affected by habitat modifications.     

Phylogeography of the red coral (<Emphasis Type="Italic">Corallium rubrum</Emphasis>): inferences on the evolutionary history of a temperate gorgonian

The red coral Corallium rubrum (Cnidaria, Octocorallia) is an exploited, long-lived sessile species from the Mediterranean Sea and the adjacent coastline in the Atlantic Ocean. Surveys of genetic variation using microsatellites have shown that populations of C. rubrum are characterized by strong differentiation at the local scale but a study of the phylogeography of this species was still lacking. Here, we used seven polymorphic microsatellite loci, together with sequence data from an intron of the elongation factor 1 (EF1) gene, to investigate the genetic structure of C. rubrum across its geographical range in the western Mediterranean Sea and in the Adriatic Sea. The EF1 sequences were also used to analyse the consequences of demographic fluctuations linked with past environmental change. Clustering analysis with microsatellite loci highlighted three to seven genetic groups with the distinction of North African and Adriatic populations; this distinction appeared significant with AMOVA and differentiation tests. Microsatellite and EF1 data extended the isolation by distance pattern previously observed for this species at the western Mediterranean scale. EF1 sequences confirmed the genetic differentiation observed between most samples with microsatellites. A statistical parsimony network of EF1 haplotypes provided no evidence of high sequence divergence among regions, suggesting no long-term isolation. Selective neutrality tests on microsatellites and EF1 were not significant but should be interpreted with caution in the case of EF1 because of the low sample sizes for this locus. Our results suggest that recent Quaternary environmental fluctuations had a limited impact on the genetic structure of C. rubrum.     

Contrasted levels of genetic diversity in a benthic Mediterranean octocoral: Consequences of different demographic histories?

Understanding the factors explaining the observed patterns of genetic diversity is an important question in evolutionary biology. We provide the first data on the genetic structure of a Mediterranean octocoral, the yellow gorgonian Eunicella cavolini, along with insights into the demographic history of this species. We sampled populations in four areas of the Mediterranean Sea: continental France, Algeria, Turkey, and the Balearic and Corsica islands. Along French coasts, three sites were sampled at two depths (20 and 40 m). We demonstrated a high genetic structure in this species (overall F_ST = 0.13), and most pairwise differentiation tests were significant. We did not detect any difference between depths at the same site. Clustering analyses revealed four differentiated groups corresponding to the main geographical areas. The levels of allelic richness and heterozygosity were significantly different between regions, with highest diversity in Algeria and lowest levels in Turkey. The highest levels of private allelic richness were observed in Algeria followed by Turkey. Such contrasted patterns of genetic diversity were not observed in other Mediterranean octocorals and could be the result of different evolutionary histories. We also provide new empirical evidence of contrasting results between tests and model‐based studies of demographic history. Our results have important consequences for the management of this species.     

Genetic population structure of turbot (Scophthalmus maximus L.) supports the presence of multiple hybrid zones for marine fishes in the transition zone between the Baltic Sea and the North Sea

Genetic population structure of turbot (Scophthalmus maximus L.) in the Northeast Atlantic was investigated using eight highly variable microsatellite loci. In total 706 individuals from eight locations with temporal replicates were assayed, covering an area from the French Bay of Biscay to the Aaland archipelago in the Baltic Sea. In contrast to previous genetic studies of turbot, we found significant genetic differentiation among samples with a maximum pairwise FST of 0.032. Limited or no genetic differentiation was found among samples within the Atlantic/North Sea area and within the Baltic Sea, suggesting high gene flow among populations in these areas. In contrast, there was a sharp cline in genetic differentiation going from the low saline Baltic Sea to the high saline North Sea. The data were explained best by two divergent populations connected by a hybrid zone; however, a mechanical mixing model could not be ruled out. A significant part of the genetic variance could be ascribed to variation among years within locality. Nevertheless, the population structure was relatively stable over time, suggesting that the observed pattern of genetic differentiation is biologically significant. This study suggests that hybrid zones are a common phenomenon for marine fishes in the transition area between the North Sea and the Baltic Sea and highlights the importance of using interspecific comparisons for inferring population structure in high gene flow species such as most marine fishes.     

Parallel genetic divergence among coastal–marine ecotype pairs of European anchovy explained by differential introgression after secondary contact

Ecophenotypic differentiation among replicate ecotype pairs within a species complex is often attributed to independent outcomes of parallel divergence driven by adaptation to similar environmental contrasts. However, the extent to which parallel phenotypic and genetic divergence patterns have emerged independently is increasingly questioned by population genomic studies. Here, we document the extent of genetic differentiation within and among two geographic replicates of the coastal and marine ecotypes of the European anchovy (Engraulis encrasicolus) gathered from Atlantic and Mediterranean locations. Using a genome‐wide data set of RAD‐derived SNPs, we show that habitat type (marine vs. coastal) is the most important component of genetic differentiation among populations of anchovy. By analysing the joint allele frequency spectrum of each coastal–marine ecotype pair, we show that genomic divergence patterns between ecotypes can be explained by a postglacial secondary contact following a long period of allopatric isolation (c. 300 kyrs). We found strong support for a model including heterogeneous migration among loci, suggesting that secondary gene flow has eroded past differentiation at different rates across the genome. Markers experiencing reduced introgression exhibited strongly correlated differentiation levels among Atlantic and Mediterranean regions. These results support that partial reproductive isolation and parallel genetic differentiation among replicate pairs of anchovy ecotypes are largely due to a common divergence history prior to secondary contact. They moreover provide comprehensive insights into the origin of a surprisingly strong fine‐scale genetic structuring in a high gene flow marine fish, which should improve stock management and conservation actions.     

Genetic Diversity and Local Connectivity in the Mediterranean Red Gorgonian Coral after Mass Mortality Events

Estimating the patterns of connectivity in marine taxa with planktonic dispersive stages is a challenging but crucial task because of its conservation implications. The red gorgonian Paramuricea clavata is a habitat forming species, characterized by short larval dispersal and high reproductive output, but low recruitment. In the recent past, the species was impacted by mass mortality events caused by increased water temperatures in summer. In the present study, we used 9 microsatellites to investigate the genetic structure and connectivity in the highly threatened populations from the Ligurian Sea (NW Mediterranean). No evidence for a recent bottleneck neither decreased genetic diversity in sites impacted by mass mortality events were found. Significant IBD pattern and high global F_ST confirmed low larval dispersal capability in the red gorgonian. The maximum dispersal distance was estimated at 20–60 km. Larval exchange between sites separated by hundreds of meters and between different depths was detected at each site, supporting the hypothesis that deeper subpopulations unaffected by surface warming peaks may provide larvae for shallower ones, enabling recovery after climatically induced mortality events.     

Cryptic habitats and cryptic diversity: unexpected patterns of connectivity and phylogeographical breaks in a Mediterranean endemic marine cave mysid

The marine cave‐dwelling mysid Hemimysis margalefi is distributed over the whole Mediterranean Sea, which contrasts with the poor dispersal capabilities of this brooding species. In addition, underwater marine caves are a highly fragmented habitat which further promotes strong genetic structuring, therefore providing highly informative data on the levels of marine population connectivity across biogeographical regions. This study investigates how habitat and geography have shaped the connectivity network of this poor disperser over the entire Mediterranean Sea through the use of several mitochondrial and nuclear markers. Five deeply divergent lineages were observed among H. margalefi populations resulting from deep phylogeographical breaks, some dating back to the Oligo‐Miocene. Whether looking at the intralineage or interlineage levels, H. margalefi populations present a high genetic diversity and population structuring. This study suggests that the five distinct lineages observed in H. margalefi actually correspond to as many separate cryptic taxa. The nominal species, H. margalefi sensu stricto, corresponds to the westernmost lineage here surveyed from the Alboran Sea to southeastern Italy. Typical genetic breaks such as the Almeria‐Oran Front or the Siculo‐Tunisian Strait do not appear to be influential on the studied loci in H. margalefi sensu stricto. Instead, population structuring appears more complex and subtle than usually found for model species with a pelagic dispersal phase. The remaining four cryptic taxa are all found in the eastern basin, but incomplete lineage sorting is suspected and speciation might still be in process. Present‐day population structure of the different H. margalefi cryptic species appears to result from past vicariance events started in the Oligo‐Miocene and maintained by present‐day coastal topography, water circulation and habitat fragmentation.