The three-spined stickleback (Gasterosteus aculeatus) is an important model system for the study of parallel evolution in the wild, having repeatedly colonized and adapted to freshwater from the sea throughout the northern hemisphere. Previous studies identified numerous genomic regions showing consistent genetic differentiation between freshwater and marine ecotypes, but these had typically limited geographic sampling and mostly focused on the Eastern Pacific region. We analysed population genomic data from the three-spined stickleback marine and freshwater ecotypes covering the entire species' range to detect loci involved in parallel evolution at different geographic scales. Most signatures of parallel evolution were unique to the Eastern Pacific and trans-oceanic 2 marine-freshwater differentiation was restricted to a limited number of shared genomic regions, including three chromosomal inversions. Based on simulations and empirical data, we demonstrate that this could result from the stochastic loss of freshwater-adapted alleles during the invasion of the Atlantic basin and selection against freshwater-adapted variants in the sea, both of which can reduce standing genetic variation available for freshwater adaptation outside the Eastern Pacific region. Moreover, the elevated linkage disequilibrium associated with marine-freshwater differentiation in the Eastern Pacific is consistent with secondary contact between marine and freshwater populations that evolved in isolation from each other during past glacial periods. Thus, contrary to what earlier studies from the Eastern Pacific region have led us to believe, parallel marinefreshwater differentiation in sticklebacks is far less prevalent and pronounced in all other parts of the species global distribution range.

Divergent selection may initiate ecological speciation extremely rapidly. How often and at what pace ecological speciation proceeds to yield strong reproductive isolation is more uncertain. Here, we document a case of extraordinarily rapid speciation associated with ecological selection in the postglacial Baltic Sea. European flounders (Platichthys flesus) in the Baltic exhibit two contrasting reproductive behaviors: pelagic and demersal spawning. Demersal spawning enables flounders to thrive in the low salinity of the Northern Baltic, where eggs cannot achieve neutral buoyancy. We show that demersal and pelagic flounders are a species pair arising from a recent event of speciation. Despite having a parapatric distribution with extensive overlap, the two species are reciprocally monophyletic and show strongly bimodal genotypic clustering and no evidence of contemporary migration, suggesting strong reproductive isolation. Divergence across the genome is weak but shows strong signatures of selection, a pattern suggestive of a recent ecological speciation event. We propose that spawning behavior in Baltic flounders is the trait under ecologically based selection causing reproductive isolation, directly implicating a process of ecological speciation. We evaluated different possible evolutionary scenarios under the approximate Bayesian computation framework and estimate that the speciation process started in allopatry ∼2,400 generations ago, following the colonization of the Baltic by the demersal lineage. This is faster than most known cases of ecological speciation and represents the most rapid event of speciation ever reported for any marine vertebrate.ecological speciation | genomics | evolution | rapid speciation | Baltic Sea

Stickleback fishes in the family Gasterosteidae have become model organisms in ecology and evolutionary biology. However, even in the case of the most widely studied species in this family - the three-spined stickleback (Gasterosteus aculeatus) - the worldwide phylogenetic relationships and colonization history of the different populations and lineages remain poorly resolved. Using a large collection of samples covering most parts of the species distribution range, we subjected thousands of single nucleotide polymorphisms to coalescent analyses in order to reconstruct a robust worldwide phylogeny of extant G. aculeatus populations, as well as their ancestral geographic distributions using Statistical-Dispersal Vicariance and Bayesian Binary MCMC analyses. The results suggest that contemporary populations originated from the Pacific Ocean in the Late Pleistocene, and the Atlantic was colonized through the Arctic Ocean by a lineage that diverged from Pacific sticklebacks ca 44.6 Kya. This lineage contains two branches: one that is distributed in the Mediterranean area, from the Iberian Peninsula to the Black Sea ('Southern European Clade'), and another that is comprised of populations from northern Europe and the east coast of North America ('Trans-Atlantic Clade'). Hence, the results suggest that the North American East Coast was colonized by trans-Atlantic migration. Coalescence-based divergence time estimates suggest that divergence among major clades is much more recent than previously estimated.

Hybridization and convergent evolution are phenomena of broad interest in evolutionary biology, but their occurrence poses challenges for reconstructing evolutionary affinities among affected taxa. Sticklebacks in the genus Pungitius are a case in point: evolutionary relationships and taxonomic validity of different species and populations in this circumpolarly distributed species complex remain contentious due to convergent evolution of traits regarded as diagnostic in their taxonomy, and possibly also due to frequent hybridization among taxa. To clarify the evolutionary relationships among different Pungitius species and populations globally, as well as to study the prevalence and extent of introgression among recognized species, genomic data sets of both reference genome‐anchored single nucleotide polymorphisms and de novo assembled RAD‐tag loci were constructed with RAD‐seq data. Both data sets yielded topologically identical and well‐supported species trees. Incongruence between nuclear and mitochondrial DNA‐based trees was found and suggested possibly frequent hybridization and mitogenome capture during the evolution of Pungitius sticklebacks. Further analyses revealed evidence for frequent nuclear genetic introgression among Pungitius species, although the estimated proportions of autosomal introgression were low. Apart from providing evidence for frequent hybridization, the results challenge earlier mitochondrial and morphology‐based hypotheses regarding the number of species and their affinities in this genus: at least seven extant species can be recognized on the basis of genetic data. The results also shed new light on the biogeographical history of the Pungitius‐complex, including suggestion of several trans‐Arctic invasions of Europe from the Northern Pacific. The well‐resolved phylogeny should facilitate the utility of this genus as a model system for future comparative evolutionary studies.

Population genetic theory predicts that small effective population sizes (Ne) and restricted gene flow limit the potential for local adaptation. In particular, the probability of evolving similar phenotypes based on shared genetic mechanisms (i.e. parallel evolution), is expected to be reduced. We tested these predictions in a comparative genomic study of two ecologically similar and geographically co-distributed stickleback species (viz. Gasterosteus aculeatus and Pungitius pungitius). We found that P. pungitius harbours less genetic diversity and exhibits higher levels of genetic differentiation and isolation-by-distance than G. aculeatus. Conversely, G. aculeatus exhibits a stronger degree of genetic parallelism across freshwater populations than P. pungitius: 2996 vs. 379 SNPs located within 26 vs. nine genomic regions show evidence of selection in multiple freshwater populations of G. aculeatus and P. pungitius, respectively. Most regions involved in parallel evolution in G. aculeatus showed increased levels of divergence, suggestive of selection on ancient haplotypes. In contrast, haplotypes involved in freshwater adaptation in P. pungitius were younger, and often associated with reduced diversity. In accordance with theory, the results suggest that connectivity and genetic drift play crucial roles in determining the levels and geographic distribution of standing genetic variation, providing evidence that population subdivision limits local adaptation and therefore also the likelihood of parallel evolution.