Urban rendezvous along the seashore: Ports as Darwinian field labs for studying marine evolution in the Anthropocene.

Humans have built ports on all the coasts of the world, allowing people to travel, exploit the sea, and develop trade. The proliferation of these artificial habitats and the associated maritime traffic is not predicted to fade in the coming decades. Ports share common characteristics: Species find themselves in novel singular environments, with particular abiotic properties-e.g., pollutants, shading, protection from wave action-within novel communities in a melting pot of invasive and native taxa. Here, we discuss how this drives evolution, including setting up of new connectivity hubs and gateways, adaptive responses to exposure to new chemicals or new biotic communities, and hybridization between lineages that would have never come into contact naturally. There are still important knowledge gaps, however, such as the lack of experimental tests to distinguish adaptation from acclimation processes, the lack of studies to understand the putative threats of port lineages to natural populations or to better understand the outcomes and fitness effects of anthropogenic hybridization. We thus call for further research examining "biological portuarization," defined as the repeated evolution of marine species in port ecosystems under human-altered selective pressures. Furthermore, we argue that ports act as giant mesocosms often isolated from the open sea by seawalls and locks and so provide replicated life-size evolutionary experiments essential to support predictive evolutionary sciences.

Genetic parallelism between European flat oyster populations at the edge of their natural range.

Although all marine ecosystems have experienced global-scale losses, oyster reefs have shown the greatest. Therefore, substantial efforts have been dedicated to restoration of such ecosystems during the last two decades. In Europe, several pilot projects for the restoration of the native European flat oyster, Ostrea edulis, recently begun and recommendations to preserve genetic diversity and to conduct monitoring protocols have been made. In particular, an initial step is to test for genetic differentiation against homogeneity among the oyster populations potentially involved in such programs. Therefore, we conducted a new sampling of wild populations at the European scale and a new genetic analysis with 203 markers to (1) confirm and study in more detail the pattern of genetic differentiation between Atlantic and Mediterranean populations, (2) identify potential translocations that could be due to aquaculture practices and (3) investigate the populations at the fringe of the geographical range, since they seemed related despite their geographic distance. Such information should be useful to enlighten the choice of the animals to be translocated or reproduced in hatcheries for further restocking. After the confirmation of the general geographical pattern of genetic structure and the identification of one potential case of aquaculture transfer at a large scale, we were able to detect genomic islands of differentiation mainly in the form of two groups of linked markers, which could indicate the presence of polymorphic chromosomal rearrangements. Furthermore, we observed a tendency for these two islands and the most differentiated loci to show a parallel pattern of differentiation, grouping the North Sea populations with the Eastern Mediterranean and Black Sea populations, against geography. We discussed the hypothesis that this genetic parallelism could be the sign of a shared evolutionary history of the two groups of populations despite them being at the border of the distribution nowadays.

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

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

Subtle limits to connectivity revealed by outlier loci within two divergent metapopulations of the deep-sea hydrothermal gastropod Ifremeria nautilei.

Hydrothermal vents form archipelagos of ephemeral deep-sea habitats that raise interesting questions about the evolution and dynamics of the associated endemic fauna, constantly subject to extinction-recolonization processes. These metal-rich environments are coveted for the mineral resources they harbour, thus raising recent conservation concerns. The evolutionary fate and demographic resilience of hydrothermal species strongly depend on the degree of connectivity among and within their fragmented metapopulations. In the deep sea, however, assessing connectivity is difficult and usually requires indirect genetic approaches. Improved detection of fine-scale genetic connectivity is now possible based on genome-wide screening for genetic differentiation. Here, we explored population connectivity in the hydrothermal vent snail Ifremeria nautilei across its species range encompassing five distinct back-arc basins in the Southwest Pacific. The global analysis, based on 10,570 single nucleotide polymorphism (SNP) markers derived from double digest restriction-site associated DNA sequencing (ddRAD-seq), depicted two semi-isolated and homogeneous genetic clusters. Demogenetic modeling suggests that these two groups began to diverge about 70,000 generations ago, but continue to exhibit weak and slightly asymmetrical gene flow. Furthermore, a careful analysis of outlier loci showed subtle limitations to connectivity between neighbouring basins within both groups. This finding indicates that migration is not strong enough to totally counterbalance drift or local selection, hence questioning the potential for demographic resilience at this latter geographical scale. These results illustrate the potential of large genomic data sets to understand fine-scale connectivity patterns in hydrothermal vents and the deep sea.

An introgression breakthrough left by an anthropogenic contact between two ascidians.

Human-driven translocations of species have diverse evolutionary consequences such as promoting hybridization between previously geographically isolated taxa. This is well illustrated by the solitary tunicate, Ciona robusta, native to the North East Pacific and introduced in the North East Atlantic. It is now co-occurring with its congener Ciona intestinalis in the English Channel, and C. roulei in the Mediterranean Sea. Despite their long allopatric divergence, first and second generation crosses showed a high hybridization success between the introduced and native taxa in the laboratory. However, previous genetic studies failed to provide evidence of recent hybridization between C. robusta and C. intestinalis in the wild. Using SNPs obtained from ddRAD-sequencing of 397 individuals from 26 populations, we further explored the genome-wide population structure of the native Ciona taxa. We first confirmed results documented in previous studies, notably (i) a chaotic genetic structure at regional scale, and (ii) a high genetic similarity between C. roulei and C. intestinalis, which is calling for further taxonomic investigation. More importantly, and unexpectedly, we also observed a genomic hotspot of long introgressed C. robusta tracts into C. intestinalis genomes at several locations of their contact zone. Both the genomic architecture of introgression, restricted to a 1.5 Mb region of chromosome 5, and its absence in allopatric populations suggest introgression is recent and occurred after the introduction of the non-native species. Overall, our study shows that anthropogenic hybridization can be effective in promoting introgression breakthroughs between species at a late stage of the speciation continuum.

Getting out of crises: Environmental, social-ecological and evolutionary research is needed to avoid future risks of pandemics.

The implementation of One Health/EcoHealth/Planetary Health approaches has been identified as key (i) to address the strong interconnections between risk for pandemics, climate change and biodiversity loss and (ii) to develop and implement solutions to these interlinked crises. As a response to the multiple calls from scientists on that subject, we have here proposed seven long-term research questions regarding COVID-19 and emerging infectious diseases (EIDs) that are based on effective integration of environmental, ecological, evolutionary, and social sciences to better anticipate and mitigate EIDs. Research needs cover the social ecology of infectious disease agents, their evolution, the determinants of susceptibility of humans and animals to infections, and the human and ecological factors accelerating infectious disease emergence. For comprehensive investigation, they include the development of nature-based solutions to interlinked global planetary crises, addressing ethical and philosophical questions regarding the relationship of humans to nature and regarding transformative changes to safeguard the environment and human health. In support of this research, we propose the implementation of innovative multidisciplinary facilities embedded in social ecosystems locally: ecological health observatories and living laboratories. This work was carried out in the frame of the European Community project HERA (www.HERAresearchEU.eu), which aims to set priorities for an environment, climate and health research agenda in the European Union by adopting a systemic approach in the face of global environmental change.

Prevalence and polymorphism of a mussel transmissible cancer in Europe.

Transmissible cancers are parasitic malignant cell lineages that have acquired the ability to infect new hosts from the same species, or sometimes related species. First described in dogs and Tasmanian devils, transmissible cancers were later discovered in some marine bivalves affected by a leukaemia-like disease. In Mytilus mussels, two lineages of bivalve transmissible neoplasia (BTN) have been described to date (MtrBTN1 and MtrBTN2), both of which emerged in a Mytilus trossulus founder individual. Here, we performed extensive screening of genetic chimerism, a hallmark of transmissible cancer, by genotyping 106 single nucleotide polymorphisms of 5,907 European Mytilus mussels. Genetic analysis allowed us to simultaneously obtain the genotype of hosts - Mytilus edulis, M. galloprovincialis or hybrids - and the genotype of tumours of heavily infected individuals. In addition, a subset of 222 individuals were systematically genotyped and analysed by histology to screen for possible nontransmissible cancers. We detected MtrBTN2 at low prevalence in M. edulis, and also in M. galloprovincialis and hybrids although at a much lower prevalence. No MtrBTN1 or new BTN were found, but eight individuals with nontransmissible neoplasia were observed at a single polluted site on the same sampling date. We observed a diversity of MtrBTN2 genotypes that appeared more introgressed or more ancestral than MtrBTN1 and reference healthy M. trossulus individuals. The observed polymorphism is probably due to somatic null alleles caused by structural variations or point mutations in primer-binding sites leading to enhanced detection of the host alleles. Despite low prevalence, two sublineages divergent by 10% fixed somatic null alleles and one nonsynonymous mtCOI (mitochondrial cytochrome oxidase I) substitution are cospreading in the same geographical area, suggesting a complex diversification of MtrBTN2 since its emergence and host species shift.

DILS: Demographic inferences with linked selection by using ABC.

We present DILS, a deployable statistical analysis platform for conducting demographic inferences with linked selection from population genomic data using an Approximate Bayesian Computation framework. DILS takes as input single-population or two-population data sets (multilocus fasta sequences) and performs three types of analyses in a hierarchical manner, identifying: (a) the best demographic model to study the importance of gene flow and population size change on the genetic patterns of polymorphism and divergence, (b) the best genomic model to determine whether the effective size Ne and migration rate N, m are heterogeneously distributed along the genome (implying linked selection) and (c) loci in genomic regions most associated with barriers to gene flow. Also available via a Web interface, an objective of DILS is to facilitate collaborative research in speciation genomics. Here, we show the performance and limitations of DILS by using simulations and finally apply the method to published data on a divergence continuum composed by 28 pairs of Mytilus mussel populations/species.

Pre-introduction introgression contributes to parallel differentiation and contrasting hybridization outcomes between invasive and native marine mussels.

Non-native species experience novel selection pressures in introduced environments and may interbreed with native lineages. Species introductions therefore provide opportunities to investigate repeated patterns of adaptation and introgression across replicated contact zones. Here, we investigate genetic parallelism between multiple introduced populations of the invasive marine mussel, Mytilus galloprovincialis, in the absence (South Africa and California) and presence of hybridization with a native congener (Mytilus planulatus in Batemans Bay and Sydney Harbour, Australia). Repeatability in post-introduction differentiation from native-range populations varied between genetically distinct Atlantic and Mediterranean lineages, with Atlantic-derived introductions displaying high differentiation (maxFST  > 0.4) and parallelism at outlier loci. Identification of long noncoding RNA transcripts (lncRNA) additionally allowed us to clarify that parallel responses are largely limited to protein-coding loci, with lncRNAs likely evolving under evolutionary constraints. Comparisons of independent hybrid zones revealed differential introgression most strongly in Batemans Bay, with an excess of M. galloprovincialis ancestry and resistance to introgression at loci differentiating parental lineages (M. planulatus and Atlantic M. galloprovincialis). Additionally, contigs putatively introgressed with divergent alleles from a closely related species, Mytilus edulis, showed stronger introgression asymmetries compared with genome-wide trends and also diverged in parallel in both Atlantic-derived introductions. These results suggest that divergent demographic histories experienced by introduced lineages, including pre-introduction introgression, influence contemporary admixture dynamics. Our findings build on previous investigations reporting contributions of historical introgression to intrinsic reproductive architectures shared between marine lineages and illustrate that interspecific introgression history can shape differentiation between colonizing populations and their hybridization with native congeners.

The geometry and genetics of hybridization.

When divergent populations form hybrids, hybrid fitness can vary with genome composition, current environmental conditions, and the divergence history of the populations. We develop analytical predictions for hybrid fitness, which incorporate all three factors. The predictions are based on Fisher's geometric model, and apply to a wide range of population genetic parameter regimes and divergence conditions, including allopatry and parapatry, local adaptation, and drift. Results show that hybrid fitness can be decomposed into intrinsic effects of admixture and heterozygosity, and extrinsic effects of the (local) adaptedness of the parental lines. Effect sizes are determined by a handful of geometric distances, which have a simple biological interpretation. These distances also reflect the mode and amount of divergence, such that there is convergence toward a characteristic pattern of intrinsic isolation. We next connect our results to the quantitative genetics of line crosses in variable or patchy environments. This means that the geometrical distances can be estimated from cross data, and provides a simple interpretation of the "composite effects." Finally, we develop extensions to the model, involving selectively induced disequilibria, and variable phenotypic dominance. The geometry of fitness landscapes provides a unifying framework for understanding speciation, and wider patterns of hybrid fitness.

How do species barriers decay? Concordance and local introgression in mosaic hybrid zones of mussels.

The Mytilus complex of marine mussel species forms a mosaic of hybrid zones, found across temperate regions of the globe. This allows us to study 'replicated' instances of secondary contact between closely related species. Previous work on this complex has shown that local introgression is both widespread and highly heterogeneous, and has identified SNPs that are outliers of differentiation between lineages. Here, we developed an ancestry-informative panel of such SNPs. We then compared their frequencies in newly sampled populations, including samples from within the hybrid zones, and parental populations at different distances from the contact. Results show that close to the hybrid zones, some outlier loci are near to fixation for the heterospecific allele, suggesting enhanced local introgression, or the local sweep of a shared ancestral allele. Conversely, genomic cline analyses, treating local parental populations as the reference, reveal a globally high concordance among loci, albeit with a few signals of asymmetric introgression. Enhanced local introgression at specific loci is consistent with the early transfer of adaptive variants after contact, possibly including asymmetric bi-stable variants (Dobzhansky-Muller incompatibilities), or haplotypes loaded with fewer deleterious mutations. Having escaped one barrier, however, these variants can be trapped or delayed at the next barrier, confining the introgression locally. These results shed light on the decay of species barriers during phases of contact.

Anthropogenic hybridization at sea: three evolutionary questions relevant to invasive species management.

Species introductions promote secondary contacts between taxa with long histories of allopatric divergence. Anthropogenic contact zones thus offer valuable contrasts to speciation studies in natural systems where past spatial isolations may have been brief or intermittent. Investigations of anthropogenic hybridization are rare for marine animals, which have high fecundity and high dispersal ability, characteristics that contrast to most terrestrial animals. Genomic studies indicate that gene flow can still occur after millions of years of divergence, as illustrated by invasive mussels and tunicates. In this context, we highlight three issues: (i) the effects of high propagule pressure and demographic asymmetries on introgression directionality, (ii) the role of hybridization in preventing introduced species spread, and (iii) the importance of postzygotic barriers in maintaining reproductive isolation. Anthropogenic contact zones offer evolutionary biologists unprecedented large scale hybridization experiments. In addition to breaking the highly effective reproductive isolating barrier of spatial segregation, they allow researchers to explore unusual demographic contexts with strong asymmetries. The outcomes are diverse, from introgression swamping to strong barriers to gene flow, and lead to local containment or widespread invasion. These outcomes should not be neglected in management policies of marine invasive species. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Within-Generation Polygenic Selection Shapes Fitness-Related Traits across Environments in Juvenile Sea Bream.

Understanding the genetic underpinnings of fitness trade-offs across spatially variable environments remains a major challenge in evolutionary biology. In Mediterranean gilthead sea bream, first-year juveniles use various marine and brackish lagoon nursery habitats characterized by a trade-off between food availability and environmental disturbance. Phenotypic differences among juveniles foraging in different habitats rapidly appear after larval settlement, but the relative role of local selection and plasticity in phenotypic variation remains unclear. Here, we combine phenotypic and genetic data to address this question. We first report correlations of opposite signs between growth and condition depending on juvenile habitat type. Then, we use single nucleotide polymorphism (SNP) data obtained by Restriction Associated DNA (RAD) sequencing to search for allele frequency changes caused by a single generation of spatially varying selection between habitats. We found evidence for moderate selection operating at multiple loci showing subtle allele frequency shifts between groups of marine and brackish juveniles. We identified subsets of candidate outlier SNPs that, in interaction with habitat type, additively explain up to 3.8% of the variance in juvenile growth and 8.7% in juvenile condition; these SNPs also explained significant fraction of growth rate in an independent larval sample. Our results indicate that selective mortality across environments during early-life stages involves complex trade-offs between alternative growth strategies.

A single clonal lineage of transmissible cancer identified in two marine mussel species in South America and Europe.

Transmissible cancers, in which cancer cells themselves act as an infectious agent, have been identified in Tasmanian devils, dogs, and four bivalves. We investigated a disseminated neoplasia affecting geographically distant populations of two species of mussels (Mytilus chilensis in South America and M. edulis in Europe). Sequencing alleles from four loci (two nuclear and two mitochondrial) provided evidence of transmissible cancer in both species. Phylogenetic analysis of cancer-associated alleles and analysis of diagnostic SNPs showed that cancers in both species likely arose in a third species of mussel (M. trossulus), but these cancer cells are independent from the previously identified transmissible cancer in M. trossulus from Canada. Unexpectedly, cancers from M. chilensis and M. edulis are nearly identical, showing that the same cancer lineage affects both. Thus, a single transmissible cancer lineage has crossed into two new host species and has been transferred across the Atlantic and Pacific Oceans and between the Northern and Southern hemispheres.

Parallel pattern of differentiation at a genomic island shared between clinal and mosaic hybrid zones in a complex of cryptic seahorse lineages.

Diverging semi-isolated lineages either meet in narrow clinal hybrid zones, or have a mosaic distribution associated with environmental variation. Intrinsic reproductive isolation is often emphasized in the former and local adaptation in the latter, although both reduce gene flow between groups. Rarely are these two patterns of spatial distribution reported in the same study system. Here, we report that the long-snouted seahorse Hippocampus guttulatus is subdivided into discrete panmictic entities by both types of hybrid zones. Along the European Atlantic coasts, a northern and a southern lineage meet in the southwest of France where they coexist in sympatry-i.e., in the same geographical zone-with little hybridization. In the Mediterranean Sea, two lineages have a mosaic distribution, associated with lagoon-like and marine habitats. A fifth lineage was identified in the Black Sea. Genetic homogeneity over large spatial scales contrasts with isolation maintained in sympatry or close parapatry at a fine scale. A high variation in locus-specific introgression rates provides additional evidence that partial reproductive isolation must be maintaining the divergence. We find that fixed differences between lagoon and marine populations in the Mediterranean Sea belong to the most differentiated SNPs between the two Atlantic lineages, against the genome-wide pattern of structure that mostly follow geography. These parallel outlier SNPs cluster on a single chromosome-wide island of differentiation. Since Atlantic lineages do not map to lagoon-sea habitat variation, genetic parallelism at the genomic island suggests a shared genetic barrier contributes to reproductive isolation in contrasting contexts-i.e., spatial versus ecological. We discuss how a genomic hotspot of parallel differentiation could have evolved and become associated both with space and with a patchy environment in a single study system.

The hidden side of a major marine biogeographic boundary: a wide mosaic hybrid zone at the Atlantic-Mediterranean divide reveals the complex interaction between natural and genetic barriers in mussels.

The Almeria-Oran Front (AOF) is a recognised hotspot of genetic differentiation in the sea, with genetic discontinuities reported in more than 50 species. The AOF is a barrier to dispersal and an ecological boundary; both can determine the position of these genetic breaks. However, the maintenance of genetic differentiation is likely reinforced by genetic barriers. A general drawback of previous studies is an insufficient density of sampling sites at the transition zone, with a conspicuous lack of samples from the southern coastline. We analysed the fine-scale genetic structure in the mussel Mytilus galloprovincialis using a few ancestry-informative loci previously identified from genome scans. We discovered a 600-km-wide mosaic hybrid zone eastward of the AOF along the Algerian coasts. This mosaic zone provides a new twist to our understanding of the Atlantic-Mediterranean transition because it demonstrates that the two lineages can live in sympatry with ample opportunities to interbreed in a large area, but they hardly do so. This implies that some form of reproductive isolation must exist to maintain the two genetic backgrounds locally cohesive. The mosaic zone ends with an abrupt genetic shift at a barrier to dispersal in the Gulf of Bejaia, Eastern Algeria. Simulations of endogenous or exogenous selection in models that account for the geography and hydrodynamic features of the region support the hypothesis that sister hybrid zones could have been differentially trapped at two alternative barriers to dispersal and/or environmental boundaries, at Almeria in the north and Bejaia in the south. A preponderantly unidirectional north-south gene flow next to the AOF can also maintain a patch of intrinsically maintained genetic background in the south and the mosaic structure, even in the absence of local adaptation. Our results concur with the coupling hypothesis that suggests that natural barriers can explain the position of genetic breaks, while their maintenance depends on genetic barriers.

Coadapted genomes and selection on hybrids: Fisher's geometric model explains a variety of empirical patterns.

Natural selection plays a variety of roles in hybridization, speciation, and admixture. Most research has focused on two extreme cases: crosses between closely related inbred lines, where hybrids are fitter than their parents, or crosses between effectively isolated species, where hybrids suffer severe breakdown. But many natural populations must fall into intermediate regimes, with multiple types of gene interaction, and these are more difficult to study. Here, we develop a simple fitness landscape model, and show that it naturally interpolates between previous modeling approaches, which were designed for the extreme cases, and invoke either mildly deleterious recessives, or discrete hybrid incompatibilities. Our model yields several new predictions, which we test with genomic data from Mytilus mussels, and published data from plants (Zea, Populus, and Senecio) and animals (Mus, Teleogryllus, and Drosophila). The predictions are generally supported, and the model explains a number of surprising empirical patterns. Our approach enables novel and complementary uses of genome-wide datasets, which do not depend on identifying outlier loci, or "speciation genes" with anomalous effects. Given its simplicity and flexibility, and its predictive successes with a wide range of data, the approach should be readily extendable to other outstanding questions in the study of hybridization.

Analysis of Genome-Wide Differentiation between Native and Introduced Populations of the Cupped Oysters Crassostrea gigas and Crassostrea angulata.

The Pacific cupped oyster is genetically subdivided into two sister taxa, Crassostrea gigas and Crassostrea angulata, which are in contact in the north-western Pacific. The nature and origin of their genetic and taxonomic differentiation remains controversial due the lack of known reproductive barriers and the high degree of morphologic similarity. In particular, whether the presence of ecological and/or intrinsic isolating mechanisms contributes to species divergence is unknown. The recent co-introduction of both taxa into Europe offers a unique opportunity to test how genetic differentiation is maintained under new environmental and demographic conditions. We generated a pseudochromosome assembly of the Pacific oyster genome using a combination of BAC-end sequencing and scaffold anchoring to a new high-density linkage map. We characterized genome-wide differentiation between C. angulata and C. gigas in both their native and introduced ranges, and showed that gene flow between species has been facilitated by their recent co-introductions in Europe. Nevertheless, patterns of genomic divergence between species remain highly similar in Asia and Europe, suggesting that the environmental transition caused by the co-introduction of the two species did not affect the genomic architecture of their partial reproductive isolation. Increased genetic differentiation was preferentially found in regions of low recombination. Using historical demographic inference, we show that the heterogeneity of differentiation across the genome is well explained by a scenario whereby recent gene flow has eroded past differentiation at different rates across the genome after a period of geographical isolation. Our results thus support the view that low-recombining regions help in maintaining intrinsic genetic differences between the two species.

The divergence history of European blue mussel species reconstructed from Approximate Bayesian Computation: the effects of sequencing techniques and sampling strategies.

Genome-scale diversity data are increasingly available in a variety of biological systems, and can be used to reconstruct the past evolutionary history of species divergence. However, extracting the full demographic information from these data is not trivial, and requires inferential methods that account for the diversity of coalescent histories throughout the genome. Here, we evaluate the potential and limitations of one such approach. We reexamine a well-known system of mussel sister species, using the joint site frequency spectrum (jSFS) of synonymous mutations computed either from exome capture or RNA-seq, in an Approximate Bayesian Computation (ABC) framework. We first assess the best sampling strategy (number of: individuals, loci, and bins in the jSFS), and show that model selection is robust to variation in the number of individuals and loci. In contrast, different binning choices when summarizing the jSFS, strongly affect the results: including classes of low and high frequency shared polymorphisms can more effectively reveal recent migration events. We then take advantage of the flexibility of ABC to compare more realistic models of speciation, including variation in migration rates through time (i.e., periodic connectivity) and across genes (i.e., genome-wide heterogeneity in migration rates). We show that these models were consistently selected as the most probable, suggesting that mussels have experienced a complex history of gene flow during divergence and that the species boundary is semi-permeable. Our work provides a comprehensive evaluation of ABC demographic inference in mussels based on the coding jSFS, and supplies guidelines for employing different sequencing techniques and sampling strategies. We emphasize, perhaps surprisingly, that inferences are less limited by the volume of data, than by the way in which they are analyzed.

Publisher Correction: The origin and remolding of genomic islands of differentiation in the European sea bass.

The originally published version of this Article contained errors in Figure 5, whereby the sign for Spearman's rho was incorrect in panels b and c. These errors have now been corrected in both the PDF and HTML versions of the Article.