Genomic signatures of host‐associated divergence and adaptation in a coral‐eating snail,Coralliophila violacea (Kiener, 1836)

The fluid nature of the ocean, combined with planktonic dispersal of marine larvae, lowers physical barriers to gene flow. However, divergence can still occur despite gene flow if strong selection acts on populations occupying different ecological niches. Here, we examined the population genomics of an ectoparasitic snail, Coralliophila violacea (Kiener 1836), that specializes on Porites corals in the Indo‐Pacific. Previous genetic analyses revealed two sympatric lineages associated with different coral hosts. In this study, we examined the mechanisms promoting and maintaining the snails’ adaptation to their coral hosts. Genome‐wide single nucleotide polymorphism (SNP) data from type II restriction site‐associated DNA (2b‐RAD) sequencing revealed two differentiated clusters of C. violacea that were largely concordant with coral host, consistent with previous genetic results. However, the presence of some admixed genotypes indicates gene flow from one lineage to the other. Combined, these results suggest that differentiation between host‐associated lineages of C. violacea is occurring in the face of ongoing gene flow, requiring strong selection. Indeed, 2.7% of all SNP loci were outlier loci (73/2,718), indicative of divergence with gene flow, driven by adaptation of each C. violacea lineage to their specific coral hosts.     

Genetic evidence for ecological divergence in kokanee salmon

The evolution of locally adapted phenotypes among populations that experience divergent selective pressures is a central mechanism for generating and maintaining biodiversity. Recently, the advent of high‐throughput DNA sequencing technology has provided tools for investigating the genetic basis of this process in natural populations of nonmodel organisms. Kokanee, the freshwater form of sockeye salmon (Oncorhynchus nerka), occurs as two reproductive ecotypes, which differ in spawning habitat (tributaries vs. shorelines); however, outside of the spawning season the two ecotypes co‐occur in many lakes and lack diagnostic morphological characteristics. We used restriction site‐associated DNA (RAD) sequencing to identify 6145 SNPs and genotype kokanee from multiple spawning sites in Okanagan Lake (British Columbia, Canada). Outlier tests revealed 18 loci putatively under divergent selection between ecotypes, all of which exhibited temporally stable allele frequencies within ecotypes. Six outliers were annotated to sequences in the NCBI database, two of which matched genes associated with early development. There was no evidence for neutral genetic differentiation; however, outlier loci demonstrated significant structure with respect to ecotype and had high assignment accuracy in mixed composition simulations. The absence of neutral structure combined with a small number of highly divergent outlier loci is consistent with theoretical predictions for the early stages of ecological divergence. These outlier loci were then applied to a realistic fisheries scenario in which additional RAD sequencing was used to genotype kokanee collected by trawl in Okanagan Lake, providing preliminary evidence that this approach may be an effective tool for conservation and management.     

Genomic divergence between Spanish Littorina saxatilis ecotypes unravels limited admixture and extensive parallelism associated with population history

The rough periwinkle, Littorina saxatilis, is a model system for studying parallel ecological speciation in microparapatry. Phenotypically parallel wave‐adapted and crab‐adapted ecotypes that hybridize within the middle shore are replicated along the northwestern coast of Spain and have likely arisen from two separate glacial refugia. We tested whether greater geographic separation corresponding to reduced opportunity for contemporary or historical gene flow between parallel ecotypes resulted in less parallel genomic divergence. We sequenced double‐digested restriction‐associated DNA (ddRAD) libraries from individual snails from upper, mid, and low intertidal levels of three separate sites colonized from two separate refugia. Outlier analysis of 4256 SNP markers identified 34.4% sharing of divergent loci between two geographically close sites; however, these sites each shared only 9.9%–15.1% of their divergent loci with a third more‐distant site. STRUCTURE analysis revealed that genotypes from only three of 166 phenotypically intermediate mid‐shore individuals appeared to result from recent hybridization, suggesting that hybrids cannot be reliably identified using shell traits. Hierarchical AMOVA indicated that the primary source of genomic differentiation was geographic separation, but also revealed greater similarity of the same ecotype across the two geographically close sites than previously estimated with dominant markers. These results from a model system for ecological speciation suggest that genomic parallelism is affected by the opportunity for historical or contemporary gene flow between populations.     

Genomics of the divergence continuum in an African plant biodiversity hotspot,I: drivers of population divergence in Restio capensis (Restionaceae)

Understanding the drivers of population divergence, speciation and species persistence is of great interest to molecular ecology, especially for species‐rich radiations inhabiting the world's biodiversity hotspots. The toolbox of population genomics holds great promise for addressing these key issues, especially if genomic data are analysed within a spatially and ecologically explicit context. We have studied the earliest stages of the divergence continuum in the Restionaceae, a species‐rich and ecologically important plant family of the Cape Floristic Region (CFR) of South Africa, using the widespread CFR endemic Restio capensis (L.) H.P. Linder & C.R. Hardy as an example. We studied diverging populations of this morphotaxon for plastid DNA sequences and >14 400 nuclear DNA polymorphisms from Restriction site Associated DNA (RAD) sequencing and analysed the results jointly with spatial, climatic and phytogeographic data, using a Bayesian generalized linear mixed modelling (GLMM) approach. The results indicate that population divergence across the extreme environmental mosaic of the CFR is mostly driven by isolation by environment (IBE) rather than isolation by distance (IBD) for both neutral and non‐neutral markers, consistent with genome hitchhiking or coupling effects during early stages of divergence. Mixed modelling of plastid DNA and single divergent outlier loci from a Bayesian genome scan confirmed the predominant role of climate and pointed to additional drivers of divergence, such as drift and ecological agents of selection captured by phytogeographic zones. Our study demonstrates the usefulness of population genomics for disentangling the effects of IBD and IBE along the divergence continuum often found in species radiations across heterogeneous ecological landscapes.     

Phylogenetic evidence for multiple sympatric ecological diversification in a marine snail

Parallel speciation can occur when traits determining reproductive isolation evolve independently in different populations that experience a similar range of environments. However, a common problem in studies of parallel evolution is to distinguish this hypothesis from an alternative one in which different ecotypes arose only once in allopatry and now share a sympatric scenario with substantial gene flow between them. Here we show that the combination of a phylogenetic approach with life-history data is able to disentangle both hypotheses in the case of the intertidal marine snail Littorina saxatilis on the rocky shores of Galicia in northwestern Spain. In this system, numerous phenotypic and genetic differences have evolved between two sympatric ecotypes spanning a sharp ecological gradient, and as aside effect of the former have produced partial reproductive isolation. A mitochondrial phylogeny of these populations strongly suggests that the two sympatric ecotypes have originated independently several times. Building upon earlier work demonstrating size-based assortative mating as the main contributor to reproductive isolation among ecotypes, our analysis provides strong evidence that divergent selection across a sharp ecological gradient promoted the parallel divergence of body size and shape between two sympatric ecotypes. Thus, divergent selection occurring independently in different populations has produced the marine equivalent of host races, which may represent the first step in speciation.     

Genomic islands of divergence linked to ecotypic variation in sockeye salmon

Regions of the genome displaying elevated differentiation (genomic islands of divergence) are thought to play an important role in local adaptation, especially in populations experiencing high gene flow. However, the characteristics of these islands as well as the functional significance of genes located within them remain largely unknown. Here, we used data from thousands of SNPs aligned to a linkage map to investigate genomic islands of divergence in three ecotypes of sockeye salmon (Oncorhynchus nerka) from a single drainage in southwestern Alaska. We found ten islands displaying high differentiation among ecotypes. Conversely, neutral structure observed throughout the rest of the genome was low and not partitioned by ecotype. One island on linkage group So13 was particularly large and contained six SNPs with F_ST > 0.14 (average F_ST of neutral SNPs = 0.01). Functional annotation revealed that the peak of this island contained a nonsynonymous mutation in a gene involved in growth in other species (TULP4). The islands that we discovered were relatively small (80–402 Kb), loci found in islands did not show reduced levels of diversity, and loci in islands displayed slightly elevated linkage disequilibrium. These attributes suggest that the islands discovered here were likely generated by divergence hitchhiking; however, we cannot rule out the possibility that other mechanisms may have produced them. Our results suggest that islands of divergence serve an important role in local adaptation with gene flow and represent a significant advance towards understanding the genetic basis of ecotypic differentiation.     

Genomic evidence for convergent evolution of a key trait underlying divergence in island birds

Reproductive isolation can be initiated by changes in one or a few key traits that prevent random mating among individuals in a population. During the early stages of speciation, when isolation is often incomplete, there will be a heterogeneous pattern of differentiation across regions of the genome between diverging populations, with loci controlling these key traits appearing the most distinct as a result of strong diversifying selection. In this study, we used Illumina‐sequenced ddRAD tags to identify genomewide patterns of differentiation in three recently diverged island populations of the Monarcha castaneiventris flycatcher of the Solomon Islands. Populations of this species have diverged in plumage colour, and these differences in plumage colour, in turn, are used in conspecific recognition and likely important in reproductive isolation. Previous candidate gene sequencing identified point mutations in MC1R and ASIP, both known pigmentation genes, to be associated with the difference in plumage colour between islands. Here, we show that background levels of genomic differentiation based on over 70,000 SNPs are extremely low between populations of distinct plumage colour, with no loci reaching the level of differentiation found in either candidate gene. Further, we found that a phylogenetic analysis based on these SNPs produced a taxonomy wherein the two melanic populations appear to have evolved convergently, rather than from a single common ancestor, in contrast to their original classification as a single subspecies. Finally, we found evidence that the pattern of low genomic differentiation is the result of both incomplete lineage sorting and gene flow between populations.     

Genomic divergence across ecological gradients in the Central African rainforest songbird (Andropadus virens)

The little greenbul, a common rainforest passerine from sub‐Saharan Africa, has been the subject of long‐term evolutionary studies to understand the mechanisms leading to rainforest speciation. Previous research found morphological and behavioural divergence across rainforest–savannah transition zones (ecotones), and a pattern of divergence with gene flow suggesting divergent natural selection has contributed to adaptive divergence and ecotones could be important areas for rainforests speciation. Recent advances in genomics and environmental modelling make it possible to examine patterns of genetic divergence in a more comprehensive fashion. To assess the extent to which natural selection may drive patterns of differentiation, here we investigate patterns of genomic differentiation among populations across environmental gradients and regions. We find compelling evidence that individuals form discrete genetic clusters corresponding to distinctive environmental characteristics and habitat types. Pairwise F_ST between populations in different habitats is significantly higher than within habitats, and this differentiation is greater than what is expected from geographic distance alone. Moreover, we identified 140 SNPs that showed extreme differentiation among populations through a genomewide selection scan. These outliers were significantly enriched in exonic and coding regions, suggesting their functional importance. Environmental association analysis of SNP variation indicates that several environmental variables, including temperature and elevation, play important roles in driving the pattern of genomic diversification. Results lend important new genomic evidence for environmental gradients being important in population differentiation.     

Proteomic evidence of a paedomorphic evolutionary process within a marine snail species: a strategy for adapting to extreme ecological conditions?

The exposed and sheltered ecotypes of the marine snail Littorina saxatilis from European rocky shores are considered a key model system to study adaptation and ecological speciation. Previous studies showed that two ecotypes (RB and SU) of this species in NW Spain have adapted differently to different shore levels and microhabitats. In order to understand how this divergent adaptive process has been accomplished, we followed a quantitative proteomic approach to investigate the proteome variation in a number of different biological factors, that is, ecotype, ontogeny and their interactions. This approach allowed testing the hypothesis that one of the ecotypes has evolved by paedomorphosis, and also whether or not the molecular mechanisms related to ecotype differentiation are set up in early developmental stages. Additionally, the identification of some candidate proteins using mass spectrometry provides some functional insights into these evolutionary processes. Results from this study provided evidence of higher ontogenetic differentiation at proteome level in the RB (metamorphic) than in SU (paedomorphic) ecotype that point to the possibility of juvenile stage retention in this latter ecotype. The level of protein expression (proteome) differences between ecotypes maintained nearly constant from late embryonic stages to adulthood, although some proteins involved in these changes considerably differed in embryonic compared to other ontogenetic stages. Paedomorphosis may be the evolutionary response of the SU ecotype of solving the trade‐off during sexually immaturity that is caused by the evolution of small size arising from adaptation to the wave‐exposed habitat. Some potential candidate genes of adaptation related to energetic metabolism have been identified, providing a promising baseline for future functional analyses.     

Comparing geographical genetic differentiation between candidate and noncandidate loci for adaptation strengthens support for parallel ecological divergence in the marine snail Littorina saxatilis

The Galician sympatric ecotypes of Littorina saxatilis have been proposed as a model system for studying parallel ecological speciation. Such a model system makes a clear prediction: candidate loci (for divergent adaptation) should present a higher level of geographical differentiation than noncandidate (neutral) loci. We used 2356 amplified fragment length polymorphisms (AFLPs) and four microsatellite loci to identify candidate loci for ecological adaptation using the F _ST outlier method. Three per cent of the studied AFLP loci were identified as candidate loci associated with adaptation, after multitest adjustments, thus contributing to ecotype differentiation (candidate loci were not detected within ecotypes). Candidate and noncandidate loci were analysed separately at four different F _ST partitions: differences between ecotypes (overall and local), differences between localities and micro-geographical differences within ecotypes. The magnitude of F _ST differed between candidate and noncandidate loci for all partitions except in the case of microgeographical differentiation within ecotypes, and the microsatellites (putatively neutral) showed an identical pattern to noncandidate loci. Thus, variation in candidate loci is determined partially independent by divergent natural selection (in addition to stochastic forces) at each locality, while noncandidate loci are exclusively driven by stochastic forces. These results support the evolutionary history described for these particular populations, considered to be a clear example of incomplete sympatric ecological speciation.     

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.     

Evolution of adaptation through allometric shifts in a marine snail

Variation in ontogenetic development among individuals may be a major contributor to morphological variation within species. Evolution of different growth trajectories might, for example, evolve as a response to varying ecological contexts of individuals living in different environments, or by life-stage or gender differences. The intertidal periwinkle Littorina saxatilis is strongly polymorphic in shell shape. We compared ontogenetic trajectories between life stages, local populations, and sexes to understand how different morphological end points are reached during ontogeny and what might cause these differences. Applying landmark-based geometric morphometrics, we captured shell shape variation for four Swedish populations of this species. We also derived a method to visualize ontogenetic trajectories described by the relationship of size to the multivariate shape space. We found that growth trajectories differed between individuals living in different habitats, as well as between sexes and maturity stages. Males living on rocky cliffs grew isometrically throughout life, whereas females from the same habitat switched from isometric growth as juveniles to allometric growth as adults. In contrast, males and females living on boulders grew allometrically as juveniles but changed to isometric growth at sexual maturity. Thus, in this species, ontogenetic growth seems influenced by habitat-associated selection as well as by gender and age-specific selection. These differing selection regimes result in ontogenetic shifts in allometry in three of the four groups examined.     

Mitonuclear discordance in wolf spiders: Genomic evidence for species integrity and introgression

Systematists and taxonomists have benefited greatly from the emergence of molecular methods. Species identification has become straightforward through DNA barcoding and the rapid build‐up of massive DNA barcode reference libraries. In animals, mitonuclear discordance can significantly complicate the process of species identification and delimitation. The causes of mitonuclear discordance are either biological (e.g., introgression, incomplete lineage sorting, horizontal gene transfer androgenesis) or induced by operational factors (e.g., human error with specimen misidentification or incorrect species delimitation). Moreover, endosymbionts may play an important role in promoting fixation of mitochondrial genomes. Here, we study the mitonuclear discordance of wolf spiders species (Lycosidae) (independent cases from Alopecosa aculeata and Pardosa pullata groups) that share identical COI DNA barcodes. We approached the case utilizing double‐digest restriction site‐associated DNA sequencing (ddRADseq) to obtain and analyse genomic‐scale data. Our results suggest that the observed cases of mitonuclear discordance are not due to operational reasons but result from biological processes. Further analysis indicated introgression and that incomplete lineage sorting is unlikely to have been responsible for the observed discrepancy. Additional survey of endosymbionts provided ideas on further research and their role in shaping mitochondrial DNA distribution patterns. Thus, ddRADseq grants an efficient way to study the taxonomy of problematic groups with insight into underlying evolutionary processes.     

GABA is an inhibitory neurotransmitter in the neural circuit regulating metamorphosis in a marine snail

The marine mud snail, Tritia (=Ilyanassa) obsoleta, displays a biphasic life cycle. During the initial phase of early development, embryos hatch from benthic egg capsules to become weakly swimming veliger larvae. In the second phase, adult T. obsoleta are facultative carnivores and major agents of community disturbance. Metamorphosis is the irreversible developmental event that links these two life history stages. When physiologically competent, larvae can respond to appropriate environmental cues by settling onto their mudflat habitat and transforming themselves into miniature adult snails. Two neurotransmitters—serotonin and nitric oxide—have opposing effects on the metamorphic process in this species. In multiple other species of gastropod and bivalve molluscs, a third neurotransmitter, the classically inhibitory compound γ‐aminobutyric acid (GABA), can induce settlement or metamorphosis upon external application to competent larvae. In this situation, GABA is presumed to mimic the action of ligands from the juvenile environment that bind to larval chemosensory receptors and activate the metamorphic pathway. Results of our experiments contradict this commonly reported action of GABA on molluscan larvae. External application of GABA to competent larvae of T. obsoleta elicited no response, but instead attenuated the action of serotonin (5‐HT), a metamorphic inducer. Our investigations into the responses of larval T. obsoleta to multiple GABAergic reagents support our hypothesis that GABA functions internally as a neurotransmitter in the pathway that controls the initiation of metamorphosis. Our results also suggest that GABA acts directly on or downstream from serotonergic neurons to regulate the metamorphosis‐inducing effects of this neurotransmitter. © 2018 Wiley Periodicals, Inc. Develop Neurobiol 78: 736–753, 2018     

Patterns of genetic diversity and candidate genes for ecological divergence in a homoploid hybrid sunflower, Helianthus anomalus

Natural hybridization accompanied by a shift in niche preference by hybrid genotypes can lead to hybrid speciation. Natural selection may cause the fixation of advantageous alleles in the ecologically diverged hybrids, and the loci experiencing selection should exhibit a reduction in allelic diversity relative to neutral loci. Here, we analyzed patterns of genetic diversity at 59 microsatellite loci associated with expressed sequence tags (ESTs) in a homoploid hybrid sunflower species, Helianthus anomalus. We used two indices, ln RV and ln RH, to compare variation and heterozygosity (respectively) at each locus between the hybrid species and its two parental species, H. annuus and H. petiolaris. Mean values of ln RV and ln RH were significantly lower than zero, which implies that H. anomalus experienced a population bottleneck during its recent evolutionary history. After correcting for the apparent bottleneck, we found six loci with a significant reduction in variation or with heterozygosity in the hybrid species, compared to one or both of the parental species. These loci should be viewed as a ranked list of candidate loci, pending further sequencing and functional analyses. Sequence data were generated for two of the candidate loci, but population genetics tests failed to detect deviations from neutral evolution at either locus. Nonetheless, a greater than eight-fold excess of nonsynonymous substitutions was found near a putative N-myristoylation motif at the second locus (HT998), and likelihood-based models indicated that the protein has been under selection in H. anomalus in the past and, perhaps, in one or both parental species. Finally, our data suggest that selective sweeps may have united populations of H. anomalus isolated by a mountain range, indicating that even low gene-flow species may be held together by the spread of advantageous alleles.     

Genomic microsatellite adaptive divergence of wild barley by microclimatic stress in 'Evolution Canyon', Israel

We examined diversity levels and patterns of 19 nuclear microsatellites and four chloroplast microsatellites in 275 genotypes of wild barley Hordeum spontaneum, in seven stations at the ‘Evolution Canyon’ (EC) microsite, Lower Nahal Oren, Mt. Carmel, Israel. EC is sharply subdivided ecologically into a tropical savannoid, ‘African’, xeric, south‐facing slope (SFS) abutting the temperate, dense, liveoak, brushwood, ‘European’, mesic, north‐facing slope (NFS). We found the following. (i) 17 of 19 (89.5%) nuDNA simple sequence repeats (SSRs) were polymorphic across all seven subpopulations and three chDNA SSRs were polymorphic. (ii) A total of 216 nuDNA SSR alleles, with a maximum of 23 alleles in a nuclear locus, and ten chDNA SSRs, with a maximum of four alleles in a locus, were registered. (iii) There were striking and significant inter‐ and intraslope diversities, based on the 19 nuDNA SSRs, climaxing with a remarkable genetic distance between the mid‐slope stations on opposite slopes (D_A = 0.481), across a distance of 200 m. This genetic distance is as large as that between the H. spontaneum populations of Jerusalem and Sede Boqer, which are separated by 100 km (× 500 larger in transect length). (iv) Slope‐unique alleles (103 = 45.6%) were higher on the ‘European’ than on the ‘African’ slope. Slope‐specific (predominant) alleles (17) were equal on opposite slopes. (v) nuDNA SSR gene diversity was higher on the ‘European’ slope and the opposite was found for the chDNA SSR. (vi) nuDNA SSR genic differentiation was very high between opposite slopes, with Gst = 0.187; for chDNA SSR this value was 0.127. Our results are inexplicable by stochastic processes and suggest that: (i) microclimatic diversifying selection is the major evolutionary, fast‐acting, interslope force, overriding migration and drift, and (ii) ecological stress can generate local, regional and global adaptive patterns, suggesting that natural selection is a major differentiating force of both coding and noncoding SSRs linking micro‐ and macroevolutionary processes. © 2005 The Linnean Society of London, Biological Journal of the Linnean Society, 2005, 84 , 205–224.     

The close relationship between estimated divergent selection and observed differentiation supports the selective origin of a marine snail hybrid zone

To study the role of divergent selection in the differentiation of the two morphs in a hybrid zone of the intertidal snail Littorina saxatilis, we compared the strength of the divergent selection acting on a series of shell characters (as estimated by the viability of snails in a reciprocal transplant experiment) with the contribution of these characters to the phenotypic differences between the morphs. We found a close correlation between selection and differentiation, which suggests a cause-effect relationship, i.e. that all present differentiation is the result of past divergent selection. In addition, divergent selection was a very important component of the total natural selection acting on shell measures. These novel results support previous evidence, based on allozyme analysis, of a parapatric origin for this hybrid zone. We discuss possible limitations of this interpretation and the circumstances under which allopatric differentiation would produce the same results. Phenotypic analysis of divergent selection may be a useful method of investigating the evolutionary mechanisms involved in differentiation processes.     

Genomic divergence between the migratory and stationary ecotypes of Atlantic cod

Atlantic cod displays a range of phenotypic and genotypic variations, which includes the differentiation into coastal stationary and offshore migratory types of cod that co‐occur in several parts of its distribution range and are often sympatric on the spawning grounds. Differentiation of these ecotypes may involve both historical separation and adaptation to ecologically distinct environments, the genetic basis of which is now beginning to be unravelled. Genomic analyses based on recent sequencing advances are able to document genomic divergence in more detail and may facilitate the exploration of causes and consequences of genome‐wide patterns. We examined genomic divergence between the stationary and migratory types of cod in the Northeast Atlantic, using next‐generation sequencing of pooled DNA from each of two population samples. Sequence data was mapped to the published cod genome sequence, arranged in more than 6000 scaffolds (611 Mb). We identified 25 divergent scaffolds (26 Mb) with a higher than average gene density, against a backdrop of overall moderate genomic differentiation. Previous findings of localized genomic divergence in three linkage groups were confirmed, including a large (15 Mb) genomic region, which seems to be uniquely involved in the divergence of migratory and stationary cod. The results of the pooled sequencing approach support and extend recent findings based on single‐nucleotide polymorphism markers and suggest a high degree of reproductive isolation between stationary and migratory cod in the North‐east Atlantic.