Speciation is rarely observable directly. A way forward is to compare pairs of ecotypes that evolved in parallel in similar contexts but have reached different degrees of reproductive isolation. Such comparisons are possible in the marine snail Littorina saxatilis by contrasting barriers to gene flow between parallel ecotypes in Spain and Sweden. In both countries, divergent ecotypes have evolved to withstand either crab predation or wave action. Here, we explore transects spanning contact zones between the Crab and the Wave ecotypes using low-coverage whole-genome sequencing, morphological and behavioural traits. Despite parallel phenotypic divergence, distinct patterns of differentiation between the ecotypes emerged: a continuous cline in Sweden indicating a weak barrier to gene flow, but two highly genetically and phenotypically divergent, and partly spatially overlapping clusters in Spain suggesting a much stronger barrier to gene flow. The absence of Spanish early-generation hybrids supported strong isolation, but a low level of gene flow is evident from molecular data. In both countries, highly differentiated loci were located in both shared and country-specific chromosomal inversions but were also present in collinear regions. Despite being considered the same species and showing similar levels of phenotypic divergence, the Spanish ecotypes are much closer to full reproductive isolation than the Swedish ones. Barriers to gene flow of very different strengths between ecotypes within the same species might be explained by dissimilarities in the spatial arrangement of habitats, the selection gradients or the ages of the systems.

The symbiosis between clownfish and giant tropical sea anemones (Order Actiniaria) is one of the most iconic on the planet. Distributed on tropical reefs, 28 species of clownfishes form obligate mutualistic relationships with 10 nominal species of venomous sea anemones. Our understanding of the symbiosis is limited by the fact that most research has been focused on the clownfishes. Chromosome-scale reference genomes are available for all clownfish species, yet only short reads-based reference genomes are available for five species of host sea anemones. Recent studies have shown that the clownfish-hosting sea anemones belong to three distinct clades of sea anemones that have evolved symbiosis with clownfishes independently. Here we present the first high-quality long-read assemblies for three species of clownfish-hosting sea anemones belonging to each of these clades: Entacmaea quadricolor, Stichodactyla haddoni, and Radianthus doreensis. PacBio HiFi sequencing yielded 1,597,562, 3,101,773, and 1,918,148 million reads for E. quadricolor, S. haddoni, and R. doreensis, respectively. All three assemblies were highly contiguous and complete with N50 values above 4 Mb and BUSCO completeness above 95% on the Metazoa dataset. Genome structural annotation with BRAKER3 predicted 20,454, 18,948, and 17,056 protein-coding genes in E. quadricolor, S. haddoni, and R. doreensis genome, respectively. These new resources will form the basis of comparative genomic analyses that will allow us to deepen our understanding of this mutualism from the host perspective.

The symbiosis between clownfishes (or anemonefishes) and their host sea anemones ranks among the most recognizable animal interactions on the planet. Found on coral reef habitats across the Indian and Pacific Oceans, 28 recognized species of clownfishes adaptively radiated from a common ancestor to live obligately with only 10 nominal species of host sea anemones. Are the host sea anemones truly less diverse than clownfishes? Did the symbiosis with clownfishes trigger a reciprocal adaptive radiation in sea anemones, or minimally, a co-evolutionary response to the mutualism? To address these questions, we combined fine- and broad-scale biogeographic sampling with multiple independent genomic datasets for the bubble-tip sea anemone, Entacmaea quadricolor , the most common clownfish host anemone throughout the Indo-West Pacific. Fine-scale sampling and restriction site associated DNA sequencing (RADseq) throughout the Japanese Archipelago revealed three highly divergent cryptic species: two of which co-occur throughout the Ryukyu Islands and can be differentiated by the clownfish species they host. Remarkably, broader biogeographic sampling and bait-capture sequencing reveals that this pattern is not simply the result of local ecological processes unique to Japan, but part of a deeper evolutionary signal where some species of E. quadricolor serve as host to the generalist clownfish species Amphiprion clarkii and others serve as host to the specialist clownfish A. frenatus . In total, we delimit at least five cryptic species in E. quadricolor that have diversified within the last five million years. The rapid diversification of E. quadricolor combined with functional ecological and phenotypic differentiation supports the hypothesis that this may represent an adaptive radiation in response to mutualism with clownfishes. Our data indicate that clownfishes are not merely settling in locally available hosts but recruiting to specialized host lineages with which they have co-evolved. These findings have important implications for understanding how the clownfish-sea anemone symbiosis has evolved and will shape future research agendas on this iconic model system. ### Competing Interest Statement The authors have declared no competing interest.

The mutualism between clownfishes (or anemonefishes) and their giant host sea anemones are among the most immediately recognizable animal interactions on the planet and have attracted a great deal of popular and scientific attention [[1][1]-[5][2]]. However, our evolutionary understanding of this iconic symbiosis comes almost entirely from studies on clownfishes— a charismatic group of 28 described species in the genus Amphiprion [[2][3]]. Adaptation to venomous sea anemones (Anthozoa: Actiniaria) provided clownfishes with novel habitat space, ultimately triggering the adaptive radiation of the group [[2][3]]. Clownfishes diverged from their free-living ancestors 25-30 MYA with their adaptive radiation to sea anemones dating to 13.2 MYA [[2][3], [3][4]]. Far from being mere habitat space, the host sea anemones also receive substantial benefits from hosting clownfishes, making the mutualistic and co-dependent nature of the symbiosis well established [[4][5], [5][2]]. Yet the evolutionary consequences of mutualism with clownfishes have remained a mystery from the host perspective. Here we use bait-capture sequencing to fully resolve the evolutionary relationships among the 10 nominal species of clownfish-hosting sea anemones for the first time ([Figure 1][6]). Using time-calibrated divergence dating analyses we calculate divergence times of less than 25 MYA for each host species, with 9 of 10 host species having divergence times within the last 13 MYA ([Figure 1][6]). The clownfish-hosting sea anemones thus diversified coincidently with clownfishes, potentially facilitating the clownfish adaptive radiation, and providing the first strong evidence for co-evolutionary patterns in this iconic partnership.### Competing Interest StatementThe authors have declared no competing interest. [1]: #ref-1 [2]: #ref-5 [3]: #ref-2 [4]: #ref-3 [5]: #ref-4 [6]: #F1

The intertidal gastropod Littorina saxatilis is a model system to study speciation and local adaptation. The repeated occurrence of distinct ecotypes showing different levels of genetic divergence makes L. saxatilis particularly suited to study different stages of the speciation continuum in the same lineage. A major finding is the presence of several large chromosomal inversions associated with the divergence of ecotypes and, specifically, the species offers a system to study the role of inversions in this divergence. The genome of L. saxatilis is 1.35 Gb and composed of 17 chromosomes. The first reference genome of the species was assembled using Illumina data, was highly fragmented (N50 of 44 kb), and was quite incomplete, with a BUSCO completeness of 80.1% on the Metazoan dataset. A linkage map of one full-sibling family enabled the placement of 587 Mbp of the genome into 17 linkage groups corresponding to the haploid number of chromosomes, but the fragmented nature of this reference genome limited the understanding of the interplay between divergent selection and gene flow during ecotype formation. Here, we present a newly generated reference genome that is highly contiguous, with a N50 of 67 Mb and 90.4% of the total assembly length placed in 17 super-scaffolds. It is also highly complete with a BUSCO completeness of 94.1% of the Metazoa dataset. This new reference will allow for investigations into the genomic regions implicated in ecotype formation as well as better characterization of the inversions and their role in speciation.

Speciation typically occurs in a time frame too long to be observed directly. This issue can be over-come by studying pairs of populations at different points in the speciation continuum, ideally within clades so that patterns are not confounded by differences among taxa. Such comparisons are possible in the marine snail Littorina saxatilis because it shows repeated occurrence of ecotypes adapted to either crab predation or wave action that differ in age and environmental context. Here, we explored transects spanning hybrid zones between the crab and wave ecotypes to contrast barriers to gene flow in Spain and Sweden, using low coverage whole genome sequencing, shell features, and behavioural traits. The two countries showed parallel divergence but distinct patterns of differentiation between the ecotypes: a continuous cline in Sweden but two highly genetically and phenotypically divergent, and partly spatially-overlapping clusters in Spain. Spanish early-generation hybrids were not observed but a low level of gene flow still seems to occur. In both countries, highly differentiated loci are clustered in genomic regions covered by chromosomal inversions but also occur in collinear regions. Despite being the same species and showing similar levels of phenotypic divergence, the Spanish ecotypes are closer to full reproductive isolation than the Swedish ecotypes. We discuss potential mechanisms contributing to the evolution of these different levels of reproductive isolation, particularly the age of the population, the strength of selection, the spatial context, and the role of assortative mating.

Understanding population divergence that eventually leads to speciation is essential for evolutionary biology. High species diversity in the sea was regarded as a paradox when strict allopatry was considered necessary for most speciation events because geographical barriers seemed largely absent in the sea, and many marine species have high dispersal capacities. Combining genome-wide data with demographic modelling to infer the demographic history of divergence has introduced new ways to address this classical issue. These models assume an ancestral population that splits into two subpopulations diverging according to different scenarios that allow tests for periods of gene flow. Models can also test for heterogeneities in population sizes and migration rates along the genome to account, respectively, for background selection and selection against introgressed ancestry. To investigate how barriers to gene flow arise in the sea, we compiled studies modelling the demographic history of divergence in marine organisms and extracted preferred demographic scenarios together with estimates of demographic parameters. These studies show that geographical barriers to gene flow do exist in the sea but that divergence can also occur without strict isolation. Heterogeneity of gene flow was detected in most population pairs suggesting the predominance of semipermeable barriers during divergence. We found a weak positive relationship between the fraction of the genome experiencing reduced gene flow and levels of genome-wide differentiation. Furthermore, we found that the upper bound of the 'grey zone of speciation' for our dataset extended beyond that found before, implying that gene flow between diverging taxa is possible at higher levels of divergence than previously thought. Finally, we list recommendations for further strengthening the use of demographic modelling in speciation research. These include a more balanced representation of taxa, more consistent and comprehensive modelling, clear reporting of results and simulation studies to rule out nonbiological explanations for general results.

Understanding the genetic targets of natural selection is one of the most challenging goals of population genetics. Some of the earliest candidate genes were identified from associations between allozyme allele frequencies and environmental variation. One such example is the clinal polymorphism in the arginine kinase (Ak) gene in the marine snail Littorina fabalis. While other enzyme loci do not show differences in allozyme frequencies among populations, the Ak alleles are near differential fixation across repeated wave exposure gradients in Europe. Here, we use this case to illustrate how a new sequencing toolbox can be employed to characterize the genomic architecture associated with historical candidate genes. We found that the Ak alleles differ by nine nonsynonymous substitutions, which perfectly explain the different migration patterns of the allozymes during electrophoresis. Moreover, by exploring the genomic context of the Ak gene, we found that the three main Ak alleles are located on different arrangements of a putative chromosomal inversion that reaches near fixation at the opposing ends of two transects covering a wave exposure gradient. This shows Ak is part of a large (3/4 of the chromosome) genomic block of differentiation, in which Ak is unlikely to be the only target of divergent selection. Nevertheless, the nonsynonymous substitutions among Ak alleles and the complete association of one allele with one inversion arrangement suggest that the Ak gene is a strong candidate to contribute to the adaptive significance of the inversion.

Potential for, and limits to, adaptation to environmental changes are critical for resilience and risk mitigation. The Mediterranean basin is a mosaic of biodiversity-rich ecosystems long affected by human influence, whose resilience is now questioned by climate change. After reviewing the different components of biological adaptation, we present the main characteristics of marine and terrestrial biodiversity in the Mediterranean basin and of the pressures they face. Taking climatic trends into consideration, we discuss the adaptive potential of a range of ecosystems dominated by species without active dispersal. We argue that the high heterogeneity of Mediterranean landscapes and seascapes constitutes a laboratory for the study of adaptation when environmental conditions change rapidly and may provide opportunities for adaptation and adaptability of species and ecosystems. Adaptive management in the Mediterranean can and should harness the nature-based solutions offered by both ecological and evolutionary processes for increasing the resilience of ecosystems to climate change.

Habitat connectivity is central to several key evolutionary and ecological processes, having implications for the spatial structuring of marine populations. For benthic species that have no or little mobility as adults, connectivity is evaluated by analyzing the dispersal of propagules across the seascape. We investigate fine-scale connectivity of coralligenous habitats in Marseille Bay (located in the north-western Mediterranean Sea) using high-resolution cartography and a particle-tracking model ran over a range of parameters derived from the biological traits of potential target species. We present annual and seasonal means of various connectivity diagnostics measuring the retention and exchange of propagules among coralligenous patches. A synthetic view is obtained by clustering individual patches into "coralligenous provinces" (ensemble of patches tightly connected by ocean currents). When discussing our results against historical observations, we highlight some genetic breaks, sharp community changes and ecological clusters that fit well our simulated connectivity patterns. Consistent findings include reduced biodiversity along Cote Bleue, high biodiversity at Planier Island and the presence of various dispersal barriers, which evolve with dispersal durations and provide the backbone of habitat connectivity. Our results help to apprehend and test hypothesis on marine population structures, providing useful information for ecologists and conservationists.

We investigated the validity of Autonomous Reef Monitoring Structures (ARMS) as monitoring tools for hard bottoms across a wide geographic and environmental range. We deployed 36 ARMS in the northeast Atlantic, northwest Mediterranean, Adriatic and Red Sea at 7-17 m depth. After 12-16 months, community composition was inferred from photographs, in six plate-faces for each ARMS. Overall, we found a highly significant effect of sea region, site (within seas), and plate-face on community composition. Plate-faces thus represent distinct micro-habitats and provide pseudo-replicates, increasing statistical power. Within each sea region taken individually, there was also a highly significant effect of site and plate-face. Because strong effects were obtained despite the fusion of taxonomic categories at high taxonomic ranks (to ensure comparability among biogeographic provinces), ARMS photo-analysis appears a promising monitoring tool for each sea region. We recommend keeping three ARMS per site and analyzing more numerous sites within a sea region to investigate environmental effects.

Ecosystem engineering species alter the physical structure of their environment and can create or modify habitats, having a massive impact on local biodiversity. Coralligenous reefs are highly diverse habitats endemic to the Mediterranean Sea built by calcareous benthic organisms among which Crustose Coralline Algae are the main engineering species. We analyzed the diversity of Lithophyllum stictiforme or L. cabiochiae in coralligenous habitats combining a multiple barcode and a population genomics approach with seascape features. Population genomics allowed disentangling pure spatial effects from environmental effects. We found that these taxa form a complex of eight highly divergent cryptic species that are easily identifiable using classic barcode markers (psbA, LSU, COI). Three factors have a significant effect on the relative abundances of these cryptic species: the location along the French Mediterranean coast, depth and Photosynthetic Active Radiation (PAR). The analysis of around 5000 SNPs for the most abundant species revealed genetic differentiation among localities in the Bay of Marseille but no differentiation between depths within locality. Thus, the effect of depth and PAR on cryptic species communities is not a consequence of restricted connectivity but rather due to differential settlement or survival among cryptic species. This differential is more likely driven by irradiance levels rather than by pressure or temperature. Both the genetic and species diversity patterns are congruent with the main patterns of currents in the Bay. Ecological differentiation among these engineering cryptic species, sensitive to ocean warming and acidification, could have important consequences on the diversity and structure of the coralligenous communities.

Genetic diversity is crucial for species' maintenance and persistence, yet is often overlooked in conservation studies. Species diversity is more often reported due to practical constraints, but it is unknown if these measures of diversity are correlated. In marine invertebrates, adults are often sessile or sedentary and populations exchange genes via dispersal of gametes and larvae. Species with a larval period are expected to have more connected populations than those without larval dispersal. We assessed the relationship between measures of species and genetic diversity, and between dispersal ability and connectivity. We compiled data on genetic patterns and life history traits in nine species across five phyla. Sampling sites spanned 600km in the northwest Mediterranean Sea and focused on a 50-km area near Marseilles, France. Comparative population genetic approaches yielded three main results. (i) Species without larvae showed higher levels of genetic structure than species with free-living larvae, but the role of larval type (lecithotrophic or planktotrophic) was negligible. (ii) A narrow area around Marseilles, subject to offshore advection, limited genetic connectivity in most species. (iii) We identified sites with significant positive contributions to overall genetic diversity across all species, corresponding with areas near low human population densities. In contrast, high levels of human activity corresponded with a negative contribution to overall genetic diversity. Genetic diversity within species was positively and significantly linearly related to local species diversity. Our study suggests that local contribution to overall genetic diversity should be taken into account for future conservation strategies.

The great phenotypic variability and the lack of diagnostic characters in the genus Pylaiella render the systematic study of this genus problematic. In this study, we investigated the diversity of Pylaiella littoralis along the Brittany (France) coast using a DNA barcoding multilocus approach with mitochondrial (cox1, nad1, and atp9) and chloroplastic (rbcL and atpB) markers associated with a population genetics approach using 10 microsatellite markers. In addition, spatio‐temporal sampling was conducted along the Brittany coast. We sampled 140 individuals from four sites located between Saint‐Malo and Concarneau (380 km) from April to October. Mitochondrial sequence data revealed the occurrence of two sibling species, with a minimum of 2.4% divergence between them. Microsatellite genotypic data congruently revealed two well‐supported clusters matching the two mitochondrial clades of Pylaiella. Although gene flow is limited between species, occurrence of genetic admixtures in some populations suggested that reproductive isolation is not complete. Our study highlighted the complementarity of barcoding and population genetics approaches to shed light on the evolutionary processes that lead to speciation.

Marine phytoplankton produce half of the oxygen we breathe and their astounding diversity is just starting to be unraveled. Many microbial phytoplankton are thought to be phototrophic, depending solely on inorganic sources of carbon and minerals for growth rather than preying on other planktonic cells. However, there is increasing evidence that symbiotic associations, to a large extent with bacteria, are required for vitamin or nutrient uptake for many eukaryotic microalgae. Here, we use in silico approaches to look for putative symbiotic interactions by analysing the gene content of microbial communities associated with 13 different Ostreococcus tauri (Chlorophyta, Mamilleophyceae) cultures sampled from the Mediterranean Sea. While we find evidence for bacteria in all cultures, there is no ubiquitous bacterial group, and the most prevalent group, Flavobacteria, is present in 10 out of 13 cultures. Among seven of the microbiomes, we detected genes predicted to encode type 3 secretion systems (T3SS, in 6/7 microbiomes) and/or putative type 6 secretion systems (T6SS, in 4/7 microbiomes). Phylogenetic analyses show that the corresponding genes are closely related to genes of systems identified in bacterial-plant interactions, suggesting that these T3SS might be involved in cell-to-cell interactions with O. tauri.