Seascape genomics of the pink abalone (Haliotis corrugata): An insight into a cross-border species in the northeast Pacific coast.

Seascape genomics gives insight into the geographic and environmental factors shaping local adaptations. It improves the understanding of the potential effects of climate change, which is relevant to provide the basis for the international management of fishery resources. The pink abalone (Haliotis corrugata) is distributed from California, United States to Baja California Sur, Mexico, exposed to a latitudinal environmental gradient in the California Current System. Management of the pink abalone contrasts between Mexico and the United States; Mexico has an active fishery organized in four administrative areas, while the United States has kept the fishery in permanent closure since 1996. However, the impact of environmental factors on genetic variation along the species distribution remains unknown, and understanding this relationship is crucial for effective spatial management strategies. This study aims to investigate the neutral and adaptive genomic structure of H. corrugata. A total of 203 samples from 13 locations were processed using ddRADseq, and covering the species' distribution. Overall, 2,231 neutral, nine potentially adaptive and three genomic-environmental association loci were detected. The neutral structure identified two groups: 1) California, United States and 2) Baja California Peninsula, México. In addition, the adaptive structure analysis also detected two groups with genetic divergence observed at Punta Eugenia. Notably, the seawater temperature significantly correlated with the northern group (temperate) and the southern (warmer) group. This study is a valuable foundation for future research and conservation initiatives, emphasizing the importance of considering neutral and adaptive genetic factors when developing management strategies for marine species.

Climate adaptation and vulnerability of foundation species in a global change hotspot.

Climate change is altering species ranges, and relative abundances within ranges, as populations become differentially adapted and vulnerable to the climates they face. Understanding present species ranges, whether species harbour and exchange adaptive variants, and how variants are distributed across landscapes undergoing rapid change, is therefore crucial to predicting responses to future climates and informing conservation strategies. Such insights are nonetheless lacking for most species of conservation concern. We assess genomic patterns of neutral variation, climate adaptation and climate vulnerability (offsets in predicted distributions of putatively adaptive variants across present and future landscapes) for sister foundation species, the marine tubeworms Galeolaria caespitosa and Galeolaria gemineoa, in a sentinel region for climate change impacts. We find that species are genetically isolated despite uncovering sympatry in their ranges, show parallel and nonparallel signals of thermal adaptation on spatial scales smaller than gene flow across their ranges, and are predicted to face different risks of maladaptation under future temperatures across their ranges. Our findings have implications for understanding local adaptation in the face of gene flow, and generate spatially explicit predictions for climatic disruption of adaptation and species distributions in coastal ecosystems that could guide experimental validation and conservation planning.

Clones on the run: The genomics of a recently expanded partially clonal species.

Why species that in their core areas mainly reproduce sexually become enriched with clones in marginal populations ("geographic parthenogenesis") remains unclear. Earlier hypotheses have emphasized that selection might promote clonality because it protects locally adapted genotypes. On the other hand, it also hampers recombination and adaptation to changing conditions. The aim of the present study was to investigate the early stages of range expansion in a partially clonal species and what drives an increase in cloning during such expansion. We used genome-wide sequencing to investigate the origin and evolution of large clones formed in a macroalgal species (Fucus vesiculosus) during a recent expansion into the postglacial Baltic Sea. We found low but persistent clonality in core populations, while at range margins, large dominant clonal lineages had evolved repeatedly from different sexual populations. A range expansion model showed that even when asexual recruitment is less favourable than sexual recruitment in core populations, repeated bottlenecks at the expansion front can establish a genetically eroded clonal wave that spreads ahead of a sexual wave into the new area. Genetic variation decreases by drift following repeated bottlenecks at the expansion front. This results in the emerging clones having low expected heterozygosity, which corroborated our empirical observations. We conclude that Baker's Law (clones being favoured by uniparental reproductive assurance in new areas) can play an important role during range expansion in partially clonal species, resulting in a complex spatiotemporal mosaic of clonal and sexual lineages that might persist during thousands of generations.

Balanced polymorphism fuels rapid selection in an invasive crab despite high gene flow and low genetic diversity.

Adaptation across environmental gradients has been demonstrated in numerous systems with extensive dispersal, despite high gene flow and consequently low genetic structure. The speed and mechanisms by which such adaptation occurs remain poorly resolved, but are critical to understanding species spread and persistence in a changing world. Here, we investigate these mechanisms in the European green crab Carcinus maenas, a globally distributed invader. We focus on a northwestern Pacific population that spread across >12 degrees of latitude in 10 years from a single source, following its introduction <35 years ago. Using six locations spanning >1500 km, we examine genetic structure using 9376 single nucleotide polymorphisms (SNPs). We find high connectivity among five locations, with significant structure between these locations and an enclosed lagoon with limited connectivity to the coast. Among the five highly connected locations, the only structure observed was a cline driven by a handful of SNPs strongly associated with latitude and winter temperature. These SNPs are almost exclusively found in a large cluster of genes in strong linkage disequilibrium that was previously identified as a candidate for cold tolerance adaptation in this species. This region may represent a balanced polymorphism that evolved to promote rapid adaptation in variable environments despite high gene flow, and which now contributes to successful invasion and spread in a novel environment. This research suggests an answer to the paradox of genetically depauperate yet successful invaders: populations may be able to adapt via a few variants of large effect despite low overall diversity.

Genetic patterns in Montipora capitata across an environmental mosaic in Kane'ohe Bay, O'ahu, Hawai'i.

Spatial genetic structure (SGS) is important to a population's ability to adapt to environmental change. For species that reproduce both sexually and asexually, the relative contribution of each reproductive mode has important ecological and evolutionary implications because asexual reproduction can have a strong effect on SGS. Reef-building corals reproduce sexually, but many species also propagate asexually under certain conditions. To understand SGS and the relative importance of reproductive mode across environmental gradients, we evaluated genetic relatedness in almost 600 colonies of Montipora capitata across 30 environmentally characterized sites in Kane'ohe Bay, O'ahu, Hawaii, using low-depth restriction digest-associated sequencing. Clonal colonies were relatively rare overall but influenced SGS. Clones were located significantly closer to one another spatially than average colonies and were more frequent on sites where wave energy was relatively high, suggesting a strong role of mechanical breakage in their formation. Excluding clones, we found no evidence of isolation by distance within sites or across the bay. Several environmental characteristics were significant predictors of the underlying genetic variation (including degree heating weeks, time spent above 30°C, depth, sedimentation rate and wave height); however, they only explained 5% of this genetic variation. Our results show that asexual fragmentation contributes to the ecology of branching corals at local scales and that genetic diversity is maintained despite strong environmental gradients in a highly impacted ecosystem, suggesting potential for broad adaptation or acclimatization in this population.

Climate differently influences the genomic patterns of two sympatric marine fish species.

Climate influences population genetic variation in marine species. Capturing these impacts remains challenging for marine fishes which disperse over large geographical scales spanning steep environmental gradients. It requires the extensive spatial sampling of individuals or populations, representative of seascape heterogeneity, combined with a set of highly informative molecular markers capable of revealing climatic-associated genetic variations. We explored how space, dispersal and environment shape the genomic patterns of two sympatric fish species in the Mediterranean Sea, which ranks among the oceanic basins most affected by climate change and human pressure. We hypothesized that the population structure and climate-associated genomic signatures of selection would be stronger in the less mobile species, as restricted gene flow tends to facilitate the fixation of locally adapted alleles. To test our hypothesis, we genotyped two species with contrasting dispersal abilities: the white seabream Diplodus sargus and the striped red mullet Mullus surmuletus. We collected 823 individuals and used genotyping by sequencing (GBS) to detect 8,206 single nucleotide polymorphisms (SNPs) for the seabream and 2,794 for the mullet. For each species, we identified highly differentiated genomic regions (i.e. outliers) and disentangled the relative contribution of space, dispersal and environmental variables (climate, marine primary productivity) on the outliers' genetic structure to test the prevalence of gene flow and local adaptation. We observed contrasting patterns of gene flow and adaptive genetic variation between the two species. The seabream showed a distinct Alboran sea population and panmixia across the Mediterranean Sea. The mullet revealed additional differentiation within the Mediterranean Sea that was significantly correlated to summer and winter temperatures, as well as marine primary productivity. Functional annotation of the climate-associated outlier SNPs then identified candidate genes involved in heat tolerance that could be examined to further predict species' responses to climate change. Our results illustrate the key steps of a comparative seascape genomics study aiming to unravel the evolutionary processes at play in marine species, to better anticipate their response to climate change. Defining population adaptation capacities and environmental niches can then serve to incorporate evolutionary processes into species conservation planning.

Seascape genomics reveals candidate molecular targets of heat stress adaptation in three coral species.

Anomalous heat waves are causing a major decline of hard corals around the world and threatening the persistence of coral reefs. There are, however, reefs that have been exposed to recurrent thermal stress over the years and whose corals appear to have been tolerant against heat. One of the mechanisms that could explain this phenomenon is local adaptation, but the underlying molecular mechanisms are poorly known. In this work, we applied a seascape genomics approach to study heat stress adaptation in three coral species of New Caledonia (southwestern Pacific) and to uncover the molecular actors potentially involved. We used remote sensing data to characterize the environmental trends across the reef system, and sampled corals living at the most contrasted sites. These samples underwent next generation sequencing to reveal single nucleotide polymorphisms (SNPs), frequencies of which were associated with heat stress gradients. As these SNPs might underpin an adaptive role, we characterized the functional roles of the genes located in their genomic region. In each of the studied species, we found heat stress-associated SNPs located in proximity of genes involved in pathways well known to contribute to the cellular responses against heat, such as protein folding, oxidative stress homeostasis, inflammatory and apoptotic pathways, and DNA damage-repair. In some cases, the same candidate molecular targets of heat stress adaptation recurred among species. Together, these results underline the relevance and the power of the seascape genomics approach for the discovery of adaptive traits that could allow corals to persist across wider thermal ranges.

Genomic vulnerability of a dominant seaweed points to future-proofing pathways for Australia's underwater forests.

Globally, critical habitats are in decline, threatening ecological, economic and social values and prompting calls for 'future proofing' efforts that enhance resilience to climate change. Such efforts rely on predicting how neutral and adaptive genomic patterns across a species' distribution will change under future climate scenarios, but data is scant for most species of conservation concern. Here, we use seascape genomics to characterise genetic diversity, structure and gene-environmental associations in a dominant forest-forming seaweed, Phyllospora comosa, along its entire latitudinal (12° latitude), and thermal (~14°C) range. Phyllospora showed high connectivity throughout its central range, with evidence of genetic structure and potential selection associated with sea surface temperatures (SSTs) at its rear and leading edges. Rear and leading-edge populations harboured only half the genetic diversity of central populations. By modelling genetic turnover as a function of SST, we assessed the genomic vulnerability across Phyllospora's distributional range under climate change scenarios. Despite low diversity, range-edge populations were predicted to harbour beneficial adaptations to marginal conditions and overall adaptability of the species may be compromised by their loss. Assisted gene flow from range edge populations may be required to enhance adaptation and increase resilience of central and leading-edge populations under warming oceans. Understanding genomic vulnerability can inform proactive restoration and future-proofing strategies for underwater forests and ensure their persistence in changing oceans.

Distinct interspecific and intraspecific vulnerability of coastal species to global change.

Characterising and predicting species responses to anthropogenic global change is one of the key challenges in contemporary ecology and conservation. The sensitivity of marine species to climate change is increasingly being described with forecasted species distributions, yet these rarely account for population level processes such as genomic variation and local adaptation. This study compares inter- and intraspecific patterns of biological composition to determine how vulnerability to climate change, and its environmental drivers, vary across species and populations. We compare species trajectories for three ecologically important southern African marine invertebrates at two time points in the future, both at the species level, with correlative species distribution models, and at the population level, with gradient forest models. Reported range shifts are species-specific and include both predicted range gains and losses. Forecasted species responses to climate change are strongly influenced by changes in a suite of environmental variables, from sea surface salinity and sea surface temperature, to minimum air temperature. Our results further suggest a mismatch between future habitat suitability (where species can remain in their ecological niche) and genomic vulnerability (where populations retain their genomic composition), highlighting the inter- and intraspecific variability in species' sensitivity to global change. Overall, this study demonstrates the importance of considering species and population level climatic vulnerability when proactively managing coastal marine ecosystems in the Anthropocene.

Multi-model seascape genomics identifies distinct environmental drivers of selection among sympatric marine species.

BACKGROUND: As global change and anthropogenic pressures continue to increase, conservation and management increasingly needs to consider species' potential to adapt to novel environmental conditions. Therefore, it is imperative to characterise the main selective forces acting on ecosystems, and how these may influence the evolutionary potential of populations and species. Using a multi-model seascape genomics approach, we compare putative environmental drivers of selection in three sympatric southern African marine invertebrates with contrasting ecology and life histories: Cape urchin (Parechinus angulosus), Common shore crab (Cyclograpsus punctatus), and Granular limpet (Scutellastra granularis). RESULTS: Using pooled (Pool-seq), restriction-site associated DNA sequencing (RAD-seq), and seven outlier detection methods, we characterise genomic variation between populations along a strong biogeographical gradient. Of the three species, only S. granularis showed significant isolation-by-distance, and isolation-by-environment driven by sea surface temperatures (SST). In contrast, sea surface salinity (SSS) and range in air temperature correlated more strongly with genomic variation in C. punctatus and P. angulosus. Differences were also found in genomic structuring between the three species, with outlier loci contributing to two clusters in the East and West Coasts for S. granularis and P. angulosus, but not for C. punctatus. CONCLUSION: The findings illustrate distinct evolutionary potential across species, suggesting that species-specific habitat requirements and responses to environmental stresses may be better predictors of evolutionary patterns than the strong environmental gradients within the region. We also found large discrepancies between outlier detection methodologies, and thus offer a novel multi-model approach to identifying the principal environmental selection forces acting on species. Overall, this work highlights how adding a comparative approach to seascape genomics (both with multiple models and species) can elucidate the intricate evolutionary responses of ecosystems to global change.

Thermal selection as a driver of marine ecological speciation.

Intraspecific genetic structure in widely distributed marine species often mirrors the boundaries between temperature-defined bioregions. This suggests that the same thermal gradients that maintain distinct species assemblages also drive the evolution of new biodiversity. Ecological speciation scenarios are often invoked to explain such patterns, but the fact that adaptation is usually only identified when phylogenetic splits are already evident makes it impossible to rule out the alternative scenario of allopatric speciation with subsequent adaptation. We integrated large-scale genomic and environmental datasets along one of the world's best-defined marine thermal gradients (the South African coastline) to test the hypothesis that incipient ecological speciation is a result of divergence linked to the thermal environment. We identified temperature-associated gene regions in a coastal fish species that is spatially homogeneous throughout several temperature-defined biogeographic regions based on selectively neutral markers. Based on these gene regions, the species is divided into geographically distinct regional populations. Importantly, the ranges of these populations are delimited by the same ecological boundaries that define distinct infraspecific genetic lineages in co-distributed marine species, and biogeographic disjunctions in species assemblages. Our results indicate that temperature-mediated selection represents an early stage of marine ecological speciation in coastal regions that lack physical dispersal barriers.

Putatively adaptive genetic variation in the giant California sea cucumber (Parastichopus californicus) as revealed by environmental association analysis of restriction-site associated DNA sequencing data.

Understanding the spatial scale of local adaptation and the factors associated with adaptive diversity are important objectives for ecology and evolutionary biology, and have significant implications for effective conservation and management of wild populations and natural resources. In this study, we used an environmental association analysis to identify important bioclimatic variables correlated with putatively adaptive genetic variation in a benthic marine invertebrate-the giant California sea cucumber (Parastichopus californicus)-spanning coastal British Columbia and southeastern Alaska. We used a redundancy analysis (RDA) with 3,699 single nucleotide polymorphisms (SNPs) obtained using RAD sequencing to detect candidate markers associated with 11 bioclimatic variables, including sea bottom and surface conditions, across two spatial scales (entire study area and within subregions). At the broadest scale, RDA revealed 59 candidate SNPs, 86% of which were associated with mean bottom temperature. Similar patterns were identified when population structure was accounted for. Additive polygenic scores, which provide a measure of the cumulative signal across all candidate SNPs, were strongly correlated with mean bottom temperature, consistent with spatially varying selection across a thermal gradient. At a finer scale, 23 candidate SNPs were detected, primarily associated with surface salinity (26%) and bottom current velocity (17%). Our findings suggest that environmental variables may play a role as drivers of spatially varying selection for P. californicus. These results provide context for future studies to evaluate the genetic basis of local adaptation in P. californicus and help inform the relevant scales and environmental variables for in situ field studies of putative adaptive variation in marine invertebrates.

Seascapes Shaped the Local Adaptation and Population Structure of South China Coast Yellowfin Seabream (Acanthopagrus latus).

Understanding the genetic composition and regional adaptation of marine species under environmental heterogeneity and fishing pressure is crucial for responsible management. In order to understand the genetic diversity and adaptability of yellowfin seabream (Acanthopagrus latus) along southern China coast, this study was conducted a seascape genome analysis on yellowfin seabream from the ecologically diverse coast, spanning over 1600 km. A total of 92 yellowfin seabream individuals from 15 sites were performed whole-genome resequencing, and 4,383,564 high-quality single nucleotide polymorphisms (SNPs) were called. By conducting a genotype-environment association analysis, 29,951 adaptive and 4,328,299 neutral SNPs were identified. The yellowfin seabream exhibited two distinct population structures, despite high gene flow between sites. The seascape genome analysis revealed that genetic structure was influenced by a variety of factors including salinity gradients, habitat distance, and ocean currents. The frequency of allelic variation at the candidate loci changed with the salinity gradient. Annotation of these loci revealed that most of the genes are associated with osmoregulation, such as kcnab2a, kcnk5a, and slc47a1. These genes are significantly enriched in pathways associated with ion transport including G protein-coupled receptor activity, transmembrane signaling receptor activity, and transporter activity. Overall, our findings provide insights into how seascape heterogeneity affects adaptive evolution, while providing important information for regional management in yellowfin seabream populations.

Local adaptation with gene flow in a highly dispersive shark.

Adaptive divergence in response to environmental clines are expected to be common in species occupying heterogeneous environments. Despite numerous advances in techniques appropriate for non-model species, gene-environment association studies in elasmobranchs are still scarce. The bronze whaler or copper shark (Carcharhinus brachyurus) is a large coastal shark with a wide distribution and one of the most exploited elasmobranchs in southern Africa. Here, we assessed the distribution of neutral and adaptive genomic diversity in C. brachyurus across a highly heterogeneous environment in southern Africa based on genome-wide SNPs obtained through a restriction site-associated DNA method (3RAD). A combination of differentiation-based genome-scan (outflank) and genotype-environment analyses (redundancy analysis, latent factor mixed models) identified a total of 234 differentiation-based outlier and candidate SNPs associated with bioclimatic variables. Analysis of 26,299 putatively neutral SNPs revealed moderate and evenly distributed levels of genomic diversity across sites from the east coast of South Africa to Angola. Multivariate and clustering analyses demonstrated a high degree of gene flow with no significant population structuring among or within ocean basins. In contrast, the putatively adaptive SNPs demonstrated the presence of two clusters and deep divergence between Angola and all other individuals from Namibia and South Africa. These results provide evidence for adaptive divergence in response to a heterogeneous seascape in a large, mobile shark despite high levels of gene flow. These results are expected to inform management strategies and policy at the national and regional level for conservation of C. brachyurus populations.

Seascape genomics reveals limited dispersal and suggests spatially varying selection among European populations of sea lamprey (Petromyzon marinus).

Sea lamprey Petromyzon marinus is an anadromous and semelparous fish without homing behaviors. Despite being a freshwater, free-living organism for a large part of their life cycle, its adulthood is spent as a parasite of marine vertebrates. In their native European range, while it is well-established that sea lampreys comprise a single nearly-panmictic population, few studies have further explored the evolutionary history of natural populations. Here, we performed the first genome-wide characterization of sea lamprey's genetic diversity in their European natural range. The objectives were to investigate the connectivity among river basins and explore evolutionary processes mediating dispersal during the marine phase, with the sequencing of 186 individuals from 8 locations spanning the North Eastern Atlantic coast and the North Sea with double-digest RAD-sequencing, obtaining a total of 30,910 bi-allelic SNPs. Population genetic analyses reinforced the existence of a single metapopulation encompassing freshwater spawning sites within the North Eastern Atlantic and the North Sea, though the prevalence of private alleles at northern latitudes suggested some limits to the species' dispersal. Seascape genomics suggested a scenario where oxygen concentration and river runoffs impose spatially varying selection across their distribution range. Exploring associations with the abundance of potential hosts further suggested that hake and cod could also impose selective pressures, although the nature of such putative biotic interactions was unresolved. Overall, the identification of adaptive seascapes in a panmictic anadromous species could contribute to conservation practices by providing information for restoration activities to mitigate local extinctions on freshwater sites.

Invader at the edge - Genomic origins and physiological differences of round gobies across a steep urban salinity gradient.

Species invasions are a global problem of increasing concern, especially in highly connected aquatic environments. Despite this, salinity conditions can pose physiological barriers to their spread, and understanding them is important for management. In Scandinavia's largest cargo port, the invasive round goby (Neogobius melanostomus) is established across a steep salinity gradient. We used 12,937 SNPs to identify the genetic origin and diversity of three sites along the salinity gradient and round goby from western, central and northern Baltic Sea, as well as north European rivers. Fish from two sites from the extreme ends of the gradient were also acclimated to freshwater and seawater, and tested for respiratory and osmoregulatory physiology. Fish from the high-salinity environment in the outer port showed higher genetic diversity, and closer relatedness to the other regions, compared to fish from lower salinity upstream the river. Fish from the high-salinity site also had higher maximum metabolic rate, fewer blood cells and lower blood Ca2+. Despite these genotypic and phenotypic differences, salinity acclimation affected fish from both sites in the same way: seawater increased the blood osmolality and Na+ levels, and freshwater increased the levels of the stress hormone cortisol. Our results show genotypic and phenotypic differences over short spatial scales across this steep salinity gradient. These patterns of the physiologically robust round goby are likely driven by multiple introductions into the high-salinity site, and a process of sorting, likely based on behaviour or selection, along the gradient. This euryhaline fish risks spreading from this area, and seascape genomics and phenotypic characterization can inform management strategies even within an area as small as a coastal harbour inlet.

Sponge diversification in marine lakes: Implications for phylogeography and population genomic studies on sponges.

The relative influence of geography, currents, and environment on gene flow within sessile marine species remains an open question. Detecting subtle genetic differentiation at small scales is challenging in benthic populations due to large effective population sizes, general lack of resolution in genetic markers, and because barriers to dispersal often remain elusive. Marine lakes can circumvent confounding factors by providing discrete and replicated ecosystems. Using high-resolution double digest restriction-site-associated DNA sequencing (4826 Single Nucleotide Polymorphisms, SNPs), we genotyped populations of the sponge Suberites diversicolor (n = 125) to test the relative importance of spatial scales (1-1400 km), local environmental conditions, and permeability of seascape barriers in shaping population genomic structure. With the SNP dataset, we show strong intralineage population structure, even at scales <10 km (average F ST = 0.63), which was not detected previously using single markers. Most variation was explained by differentiation between populations (AMOVA: 48.8%) with signatures of population size declines and bottlenecks per lake. Although the populations were strongly structured, we did not detect significant effects of geographic distance, local environments, or degree of connection to the sea on population structure, suggesting mechanisms such as founder events with subsequent priority effects may be at play. We show that the inclusion of morphologically cryptic lineages that can be detected with the COI marker can reduce the obtained SNP set by around 90%. Future work on sponge genomics should confirm that only one lineage is included. Our results call for a reassessment of poorly dispersing benthic organisms that were previously assumed to be highly connected based on low-resolution markers.

Genomic Signatures of Local Adaptation under High Gene Flow in Lumpfish-Implications for Broodstock Provenance Sourcing and Larval Production.

Aquaculture of the lumpfish (Cyclopterus lumpus L.) has become a large, lucrative industry owing to the escalating demand for "cleaner fish" to minimise sea lice infestations in Atlantic salmon mariculture farms. We used over 10K genome-wide single nucleotide polymorphisms (SNPs) to investigate the spatial patterns of genomic variation in the lumpfish along the coast of Norway and across the North Atlantic. Moreover, we applied three genome scans for outliers and two genotype-environment association tests to assess the signatures and patterns of local adaptation under extensive gene flow. With our 'global' sampling regime, we found two major genetic groups of lumpfish, i.e., the western and eastern Atlantic. Regionally in Norway, we found marginal evidence of population structure, where the population genomic analysis revealed a small portion of individuals with a different genetic ancestry. Nevertheless, we found strong support for local adaption under high gene flow in the Norwegian lumpfish and identified over 380 high-confidence environment-associated loci linked to gene sets with a key role in biological processes associated with environmental pressures and embryonic development. Our results bridge population genetic/genomics studies with seascape genomics studies and will facilitate genome-enabled monitoring of the genetic impacts of escapees and allow for genetic-informed broodstock selection and management in Norway.

Biophysical models of dispersal contribute to seascape genetic analyses.

Dispersal is generally difficult to directly observe. Instead, dispersal is often inferred from genetic markers and biophysical modelling where a correspondence indicates that dispersal routes and barriers explain a significant part of population genetic differentiation. Biophysical models are used for wind-driven dispersal in terrestrial environments and for propagules drifting with ocean currents in the sea. In the ocean, such seascape genetic or seascape genomic studies provide promising tools in applied sciences, as actions within management and conservation rely on an understanding of population structure, genetic diversity and presence of local adaptations, all dependent on dispersal within the metapopulation. Here, we surveyed 87 studies that combine population genetics and biophysical models of dispersal. Our aim was to understand if biophysical dispersal models can generally explain genetic differentiation. Our analysis shows that genetic differentiation and lack of genetic differentiation can often be explained by dispersal, but the realism of the biophysical model, as well as local geomorphology and species biology also play a role. The review supports the use of a combination of both methods, and we discuss our findings in terms of recommendations for future studies and pinpoint areas where further development is necessary, particularly on how to compare both approaches. This article is part of the theme issue 'Species' ranges in the face of changing environments (part I)'.

Population and seascape genomics of a critically endangered benthic elasmobranch, the blue skate Dipturus batis.

The blue skate (Dipturus batis) has a patchy distribution across the North-East Atlantic Ocean, largely restricted to occidental seas around the British Isles following fisheries-induced population declines and extirpations. The viability of remnant populations remains uncertain and could be impacted by continued fishing and by-catch pressure, and the projected impacts of climate change. We genotyped 503 samples of D. batis, obtained opportunistically from the widest available geographic range, across 6 350 single nucleotide polymorphisms (SNPs) using a reduced-representation sequencing approach. Genotypes were used to assess the species' contemporary population structure, estimate effective population sizes and identify putative signals of selection in relation to environmental variables using a seascape genomics approach. We identified genetic discontinuities between inshore (British Isles) and offshore (Rockall and Faroe Island) populations, with differentiation most pronounced across the deep waters of the Rockall Trough. Effective population sizes were largest in the Celtic Sea and Rockall, but low enough to be of potential conservation concern among Scottish and Faroese sites. Among the 21 candidate SNPs under positive selection was one significantly correlated with environmental variables predicted to be affected by climate change, including bottom temperature, salinity and pH. The paucity of well-annotated elasmobranch genomes precluded us from identifying a putative function for this SNP. Nevertheless, our findings suggest that climate change could inflict a strong selective force upon remnant populations of D. batis, further constraining its already-restricted habitat. Furthermore, the results provide fundamental insights on the distribution, behaviour and evolutionary biology of D. batis in the North-East Atlantic that will be useful for the establishment of conservation actions for this and other critically endangered elasmobranchs.