Navigating the currents of seascape genomics: how spatial analyses can augment population genomic studies

Combining population genomics and biophysical modelling to assess connectivity patterns in an Antarctic fish

Connectivity is a fundamental process of population dynamics in marine ecosystems. In the last decade, with the emergence of new methods, combining different approaches to understand the patterns of connectivity among populations and their regulation has become increasingly feasible. The Western Antarctic Peninsula (WAP) is characterized by complex oceanographic dynamics, where local conditions could act as barriers to population connectivity. Here, the notothenioid fish Harpagifer antarcticus, a demersal species with a complex life cycle (adults with poor swim capabilities and pelagic larvae), was used to assess connectivity along the WAP by combining biophysical modelling and population genomics methods. Both approaches showed congruent patterns. Areas of larvae retention and low potential connectivity, observed in the biophysical model output, coincide with four genetic groups within the WAP: (1) South Shetland Islands, (2) Bransfield Strait, (3) the central and (4) the southern area of WAP (Marguerite Bay). These genetic groups exhibited limited gene flow between them, consistent with local oceanographic conditions, which would represent barriers to larval dispersal. The joint effect of geographic distance and larval dispersal by ocean currents had a greater influence on the observed population structure than each variable evaluated separately. The combined effect of geographic distance and a complex oceanographic dynamic would be generating limited levels of population connectivity in the fish H. antarcticus along the WAP. Based on this, population connectivity estimations and priority areas for conservation were discussed, considering the marine protected area proposed for this threatened region of the Southern Ocean.

Boundary Effects Cause False Signals of Range Expansions in Population Genomic Data

Studying range expansions is central for understanding genetic variation through space and time as well as for identifying refugia and biological invasions. Range expansions are characterized by serial founder events causing clines of decreasing genetic diversity away from the center of origin and asymmetries in the two-dimensional allele frequency spectra. These asymmetries, summarized by the directionality index (ψ), are sensitive to range expansions and persist for longer than clines in genetic diversity. In continuous and finite meta-populations, genetic drift tends to be stronger at the edges of the species distribution in equilibrium populations and populations undergoing range expansions alike. Such boundary effects are expected to affect geographic patterns in genetic diversity and ψ. Here we demonstrate that boundary effects cause high false positive rates in equilibrium meta-populations when testing for range expansions. In the simulations, the absolute value of ψ (|ψ|) in equilibrium data sets was proportional to the fixation index (FST). By fitting signatures of range expansions as a function of ɛ |ψ|/FST and geographic clines in ψ, strong evidence for range expansions could be detected in data from a recent rapid invasion of the cane toad, Rhinella marina, in Australia, but not in 28 previously published empirical data sets from Australian scincid lizards that were significant for the standard range expansion tests. Thus, while clinal variation in ψ is still the most sensitive statistic to range expansions, to detect true signatures of range expansions in natural populations, its magnitude needs to be considered in relation to the overall levels of genetic structuring in the data.

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.

Finding genes and pathways that underlie coral adaptation

Mass coral bleaching is one of the clearest threats of climate change to the persistence of marine biodiversity. Despite the negative impacts of bleaching on coral health and survival, some corals may be able to rapidly adapt to warming ocean temperatures. Thus, a significant focus in coral research is identifying the genes and pathways underlying coral heat adaptation. Here, we review state-of-the-art methods that may enable the discovery of heat-adaptive loci in corals and identify four main knowledge gaps. To fill these gaps, we describe an experimental approach combining seascape genomics with CRISPR/Cas9 gene editing to discover and validate heat-adaptive loci. Finally, we discuss how information on adaptive genotypes could be used in coral reef conservation and management strategies.

Genomic analyses indicate resilience of a commercially and culturally important marine gastropod snail to climate change

Genomic vulnerability analyses are being increasingly used to assess the adaptability of species to climate change and provide an opportunity for proactive management of harvested marine species in changing oceans. Southeastern Australia is a climate change hotspot where many marine species are shifting poleward. The turban snail, Turbo militaris is a commercially and culturally harvested marine gastropod snail from eastern Australia. The species has exhibited a climate-driven poleward range shift over the last two decades presenting an ongoing challenge for sustainable fisheries management. We investigate the impact of future climate change on T. militaris using genotype-by-sequencing to project patterns of gene flow and local adaptation across its range under climate change scenarios. A single admixed, and potentially panmictic, demographic unit was revealed with no evidence of genetic subdivision across the species range. Significant genotype associations with heterogeneous habitat features were observed, including associations with sea surface temperature, ocean currents, and nutrients, indicating possible adaptive genetic differentiation. These findings suggest that standing genetic variation may be available for selection to counter future environmental change, assisted by widespread gene flow, high fecundity and short generation time in this species. We discuss the findings of this study in the content of future fisheries management and conservation.

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.

Genetic Diversity, Population Structure, and Environmental Adaptation Signatures of Chinese Coastal Hard-Shell Mussel Mytilus coruscus Revealed by Whole-Genome Sequencing

The hard-shell mussel (Mytilus coruscus) is widespread in the temperate coastal areas of the northwest Pacific and holds a significant position in the shellfish aquaculture market in China. However, the natural resources of this species have been declining, and population genetic studies of M. coruscus are also lacking. In this study, we conducted whole-genome resequencing (WGR) of M. coruscus from eight different latitudes along the Chinese coast and identified a total of 25,859,986 single nucleotide polymorphism (SNP) markers. Our findings indicated that the genetic diversity of M. coruscus from the Zhoushan region was lower compared with populations from other regions. Furthermore, we observed that the evolutionary tree clustered into two primary branches, and the Zhangzhou (ZZ) population was in a separate branch. The ZZ population was partly isolated from populations in other regions, but the distribution of branches was not geographically homogeneous, and a nested pattern emerged, consistent with the population differentiation index (FST) results. To investigate the selection characteristics, we utilized the northern M. coruscus populations (Dalian and Qingdao) and the central populations (Zhoushan and Xiangshan) as reference populations and the southern ZZ population as the target population. Our selection scan analysis identified several genes associated with thermal responses, including Hsp70 and CYP450. These genes may play important roles in the adaptation of M. coruscus to different living environments. Overall, our study provides a comprehensive understanding of the genomic diversity of coastal M. coruscus in China and is a valuable resource for future studies on genetic breeding and the evolutionary adaptation of this species.

Holistic view of the seascape dynamics and environment impact on macro-scale genetic connectivity of marine plankton populations

BACKGROUND

Plankton seascape genomics studies have revealed different trends from large-scale weak differentiation to microscale structures. Previous studies have underlined the influence of the environment and seascape on species differentiation and adaptation. However, these studies have generally focused on a few single species, sparse molecular markers, or local scales. Here, we investigated the genomic differentiation of plankton at the macro-scale in a holistic approach using Tara Oceans metagenomic data together with a reference-free computational method.

RESULTS

We reconstructed the FST-based genomic differentiation of 113 marine planktonic taxa occurring in the North and South Atlantic Oceans, Southern Ocean, and Mediterranean Sea. These taxa belong to various taxonomic clades spanning Metazoa, Chromista, Chlorophyta, Bacteria, and viruses. Globally, population genetic connectivity was significantly higher within oceanic basins and lower in bacteria and unicellular eukaryotes than in zooplankton. Using mixed linear models, we tested six abiotic factors influencing connectivity, including Lagrangian travel time, as proxies of oceanic current effects. We found that oceanic currents were the main population genetic connectivity drivers, together with temperature and salinity. Finally, we classified the 113 taxa into parameter-driven groups and showed that plankton taxa belonging to the same taxonomic rank such as phylum, class or order presented genomic differentiation driven by different environmental factors.

CONCLUSION

Our results validate the isolation-by-current hypothesis for a non-negligible proportion of taxa and highlight the role of other physicochemical parameters in large-scale plankton genetic connectivity. The reference-free approach used in this study offers a new systematic framework to analyse the population genomics of non-model and undocumented marine organisms from a large-scale and holistic point of view.

Large effective population size masks population genetic structure in Hirondellea amphipods within the deepest marine ecosystem, the Mariana Trench

The examination of genetic structure in the deep-ocean hadal zone has focused on divergence between tectonic trenches to understand how environment and geography may drive species divergence and promote endemism. There has been little attempt to examine localized genetic structure within trenches, partly because of logistical challenges associated with sampling at an appropriate scale, and the large effective population sizes of species that can be sampled adequately may mask underlying genetic structure. Here we examine genetic structure in the superabundant amphipod Hirondellea gigas in the Mariana Trench at depths of 8126-10,545 m. RAD sequencing was used to identify 3182 loci containing 43,408 single nucleotide polymorphisms (SNPs) across individuals after stringent pruning of loci to prevent paralogous multicopy genomic regions being erroneously merged. Principal components analysis of SNP genotypes resolved no genetic structure between sampling locations, consistent with a signature of panmixia. However, discriminant analysis of principal components identified divergence between all sites driven by 301 outlier SNPs in 169 loci and significantly associated with latitude and depth. Functional annotation of loci identified differences between singleton loci used in analysis and paralogous loci pruned from the data set and also between outlier and nonoutlier loci, all consistent with hypotheses explaining the role of transposable elements driving genome dynamics. This study challenges the traditional perspective that highly abundant amphipods within a trench form a single panmictic population. We discuss the findings in relation to eco-evolutionary and ontogenetic processes operating in the deep sea, and highlight key challenges associated with population genetic analysis in nonmodel systems with inherent large effective population sizes and genomes.

A decade of progress in marine evolutionary biology

This article summarizes the Evolutionary Applications Special Issue, "A decade of progress in Marine Evolutionary Biology." The globally connected ocean, from its pelagic depths to its highly varied coastlines, inspired Charles Darwin to develop the theory of evolution during the voyage of the Beagle. As technology has developed, there has been a dramatic increase in our knowledge about life on our blue planet. This Special Issue, composed of 19 original papers and seven reviews, represents a small contribution to the larger picture of recent research in evolutionary biology, and how such advancements come about through the connection of researchers, their fields, and their knowledge. The first European network for marine evolutionary biology, the Linnaeus Centre for Marine Evolutionary Biology (CeMEB), was developed to study evolutionary processes in the marine environment under global change. Though hosted by the University of Gothenburg in Sweden, the network quickly grew to encompass researchers throughout Europe and beyond. Today, more than a decade after its foundation, CeMEB's focus on the evolutionary consequences of global change is more relevant than ever, and knowledge gained from marine evolution research is urgently needed in management and conservation. This Special Issue, organized and developed through the CeMEB network, contains contributions from all over the world and provides a snapshot of the current state of the field, thus forming an important basis for future research directions.

Population structure of mud flounder Paralichthys orbignyanus from the south-western Atlantic Ocean

The mud flounder Paralichthys orbignyanus (Pleuronectiformes, Paralichthyidae) inhabits shallow waters of low salinities and mud bottoms in the temperate marine coastal regions of the Bonaerensean Ecoregion of the Argentinean Biogeographic Province in the south-western Atlantic Ocean. Specimens of P. orbignyanus were collected from Lagoa dos Patos (LDP) (southern Brazil), Mar Chiquita (MCH) and Marisol (MAR) both located in Buenos Aires (Argentina), and San Antonio Oeste (SAO) in the San Matías Gulf, Rio Negro (Argentina). A fragment of the mitochondrial DNA of the Control Region and seven microsatellite loci were characterized. In the Control Region, P. orbignyanus showed high variability, low nucleotide diversity, mild population expansion and a coalescence time of 35,000 years before the present. Flounders provided evidence of a genetic structure between the sampling sites LDP, MCH, MAR vs. SAO. On the other hand, P. orbignyanus displayed a lower to moderate contemporary genetic structure among all samples except between LDP and MCH. With no evidence of isolation by distance, this analysis supports a model of limited gene flow that is likely to be associated with a consistent larvae retention in all sampling sites. In addition, the present connectivity is ascribed to a lower migration process from SAO in the San Matías Gulf congruent with the prevailing littoral drift.

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.

Sardines at a junction: Seascape genomics reveals ecological and oceanographic drivers of variation in the NW Mediterranean Sea

By evaluating genetic variation across the entire genome, one can address existing questions in a novel way while raising new ones. The latter includes how different local environments influence adaptive and neutral genomic variation within and among populations, providing insights into local adaptation of natural populations and their responses to global change. Here, under a seascape genomic approach, ddRAD data of 4609 single nucleotide polymorphisms (SNPs) from 398 sardines (Sardina pilchardus) collected in 11 Mediterranean and one Atlantic site were generated. These were used along with oceanographic and ecological information to detect signals of adaptive divergence with gene flow across environmental gradients. The studied sardines constitute two clusters (F_ST  = 0.07), a pattern attributed to outlier loci, highlighting putative local adaptation. The trend in the number of days with sea surface temperature above 19°C, a critical threshold for successful sardine spawning, was crucial at all levels of population structuring with implications on the species' key biological processes. Outliers link candidate SNPs to the region's environmental heterogeneity. Our findings provide evidence for a dynamic equilibrium in which population structure is maintained by physical and ecological factors under the opposing influences of migration and selection. This dynamic in a natural system warrants continuous monitoring under a seascape genomic approach that might benefit from a temporal and more detailed spatial dimension. Our results may contribute to complementary studies aimed at providing deeper insights into the mechanistic processes underlying population structuring. Those are key to understanding and predicting future changes and responses of this highly exploited species in the face of climate change.

Seascape genomics of common dolphins (Delphinus delphis) reveals adaptive diversity linked to regional and local oceanography

Background

High levels of standing genomic variation in wide-ranging marine species may enhance prospects for their long-term persistence. Patterns of connectivity and adaptation in such species are often thought to be influenced by spatial factors, environmental heterogeneity, and oceanographic and geomorphological features. Population-level studies that analytically integrate genome-wide data with environmental information (i.e., seascape genomics) have the potential to inform the spatial distribution of adaptive diversity in wide-ranging marine species, such as many marine mammals. We assessed genotype-environment associations (GEAs) in 214 common dolphins (Delphinus delphis) along > 3000 km of the southern coast of Australia.

Results

We identified 747 candidate adaptive SNPs out of a filtered panel of 17,327 SNPs, and five putatively locally-adapted populations with high levels of standing genomic variation were disclosed along environmentally heterogeneous coasts. Current velocity, sea surface temperature, salinity, and primary productivity were the key environmental variables associated with genomic variation. These environmental variables are in turn related to three main oceanographic phenomena that are likely affecting the dispersal of common dolphins: (1) regional oceanographic circulation, (2) localised and seasonal upwellings, and (3) seasonal on-shelf circulation in protected coastal habitats. Signals of selection at exonic gene regions suggest that adaptive divergence is related to important metabolic traits.

Conclusion

To the best of our knowledge, this represents the first seascape genomics study for common dolphins (genus Delphinus). Information from the associations between populations and their environment can assist population management in forecasting the adaptive capacity of common dolphins to climate change and other anthropogenic impacts.

Genomic differentiation of three pico-phytoplankton species in the Mediterranean Sea

For more than a decade, high-throughput sequencing has transformed the study of marine planktonic communities and has highlighted the extent of protist diversity in these ecosystems. Nevertheless, little is known relative to their genomic diversity at the species-scale as well as their major speciation mechanisms. An increasing number of data obtained from global scale sampling campaigns is becoming publicly available, and we postulate that metagenomic data could contribute to deciphering the processes shaping protist genomic differentiation in the marine realm. As a proof of concept, we developed a findable, accessible, interoperable and reusable (FAIR) pipeline and focused on the Mediterranean Sea to study three a priori abundant protist species: Bathycoccus prasinos, Pelagomonas calceolata and Phaeocystis cordata. We compared the genomic differentiation of each species in light of geographic, environmental and oceanographic distances. We highlighted that isolation-by-environment shapes the genomic differentiation of B. prasinos, whereas P. cordata is impacted by geographic distance (i.e. isolation-by-distance). At present time, the use of metagenomics to accurately estimate the genomic differentiation of protists remains challenging since coverages are lower compared to traditional population surveys. However, our approach sheds light on ecological and evolutionary processes occurring within natural marine populations and paves the way for future protist population metagenomic studies.

Population connectivity and genetic offset in the spawning coral Acropora digitifera in Western Australia

Anthropogenic climate change has caused widespread loss of species biodiversity and ecosystem productivity across the globe, particularly on tropical coral reefs. Predicting the future vulnerability of reef-building corals, the foundation species of coral reef ecosystems, is crucial for cost-effective conservation planning in the Anthropocene. In this study, we combine regional population genetic connectivity and seascape analyses to explore patterns of genetic offset (the mismatch of gene-environmental associations under future climate conditions) in Acropora digitifera across 12 degrees of latitude in Western Australia. Our data revealed a pattern of restricted gene flow and limited genetic connectivity among geographically distant reef systems. Environmental association analyses identified a suite of loci strongly associated with the regional temperature variation. These loci helped forecast future genetic offset in gradient forest and generalised dissimilarity models. These analyses predicted pronounced differences in the response of different reef systems in Western Australia to rising temperatures. Under the most optimistic future warming scenario (RCP 2.6), we predicted a general pattern of increasing genetic offset with latitude. Under the extreme climate scenario (RCP 8.5 in 2090-2100), coral populations at the Ningaloo World Heritage Area were predicted to experience a higher mismatch between current allele frequencies and those required to cope with local environmental change, compared to populations in the inshore Kimberley region. The study suggests complex and spatially heterogeneous patterns of climate-change vulnerability in coral populations across Western Australia, reinforcing the notion that regionally tailored conservation efforts will be most effective at managing coral reef resilience into the future.

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)'.

Combining population genomics with demographic analyses highlights habitat patchiness and larval dispersal as determinants of connectivity in coastal fish species

Gene flow shapes spatial genetic structure and the potential for local adaptation. Among marine animals with nonmigratory adults, the presence or absence of a pelagic larval stage is thought to be a key determinant in shaping gene flow and the genetic structure of populations. In addition, the spatial distribution of suitable habitats is expected to influence the distribution of biological populations and their connectivity patterns. We used whole genome sequencing to study demographic history and reduced representation (double‐digest restriction associated DNA) sequencing data to analyse spatial genetic structure in broadnosed pipefish (Syngnathus typhle). Its main habitat is eelgrass beds, which are patchily distributed along the study area in southern Norway. Demographic connectivity among populations was inferred from long‐term (~30‐year) population counts that uncovered a rapid decline in spatial correlations in abundance with distance as short as ~2 km. These findings were contrasted with data for two other fish species that have a pelagic larval stage (corkwing wrasse, Symphodus melops; black goby, Gobius niger). For these latter species, we found wider spatial scales of connectivity and weaker genetic isolation‐by‐distance patterns, except where both species experienced a strong barrier to gene flow, seemingly due to lack of suitable habitat. Our findings verify expectations that a fragmented habitat and absence of a pelagic larval stage promote genetic structure, while presence of a pelagic larvae stage increases demographic connectivity and gene flow, except perhaps over extensive habitat gaps.

Seascape genomics of coastal bottlenose dolphins along strong gradients of temperature and salinity

Heterogeneous seascapes and strong environmental gradients in coastal waters are expected to influence adaptive divergence, particularly in species with large population sizes where selection is expected to be highly efficient. However, these influences might also extend to species characterized by strong social structure, natal philopatry and small home ranges. We implemented a seascape genomic study to test this hypothesis in Indo-Pacific bottlenose dolphins (Tursiops aduncus) distributed along the environmentally heterogeneous coast of southern Australia. The datasets included oceanographic and environmental variables thought to be good predictors of local adaptation in dolphins and 8,081 filtered single nucleotide polymorphisms (SNPs) genotyped for individuals sampled from seven different bioregions. From a neutral perspective, population structure and connectivity of the dolphins were generally influenced by habitat type and social structuring. Genotype-environment association analysis identified 241 candidate adaptive loci and revealed that sea surface temperature and salinity gradients influenced adaptive divergence in these animals at both large- (1,000s km) and fine-scales (<100 km). Enrichment analysis and annotation of candidate genes revealed functions related to sodium-activated ion transport, kidney development, adipogenesis and thermogenesis. The findings of spatial adaptive divergence and inferences of putative physiological adaptations challenge previous suggestions that marine megafauna is most likely to be affected by environmental and climatic changes via indirect, trophic effects. Our work contributes to conservation management of coastal bottlenose dolphins subjected to anthropogenic disturbance and to efforts of clarifying how seascape heterogeneity influences adaptive diversity and evolution in small cetaceans.

Seascape genomics identify adaptive barriers correlated to tidal amplitude in the shore crab Carcinus maenas

Most marine invertebrates disperse during a planktonic larval stage that may drift for weeks with ocean currents. A challenge for larvae of coastal species is to return to coastal nursery habitats. Shore crab (Carcinus maenas L.) larvae are known to show tidal rhythmicity in vertical migration in tidal areas and circadian rhythmicity in microtidal areas, which seems to increase successful coastal settlement. We studied genome‐wide differentiation based on 24,000 single nucleotide polymorphisms of 12 native populations of shore crab sampled from a large tidal amplitude gradient from macrotidal (~8 m) to microtidal (~0.2 m). Dispersal and recruitment success of larvae was assessed with a Lagrangian biophysical model, which showed a strong effect of larval behaviour on long‐term connectivity, and dispersal barriers that partly coincided with different tidal environments. The genetic population structure showed a subdivision of the samples into three clusters, which represent micro‐, meso‐ and macrotidal areas. The genetic differentiation was mostly driven by 0.5% outlier loci, which showed strong allelic clines located at the limits between the three tidal areas. Demographic modelling suggested that the two genetic barriers have different origins. Differential gene expression of two clock genes (cyc and pdp1) further highlighted phenotypic differences among genetic clusters that are potentially linked to the differences in larval behaviour. Taken together, our seascape genomic study suggests that tidal regime acts as a strong selection force on shore crab population structure, consistent with larval behaviour affecting dispersal and recruitment success.

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.

Applying genomics in assisted migration under climate change: Framework, empirical applications, and case studies

The rate of global climate change is projected to outpace the ability of many natural populations and species to adapt. Assisted migration (AM), which is defined as the managed movement of climate-adapted individuals within or outside the species ranges, is a conservation option to improve species' adaptive capacity and facilitate persistence. Although conservation biologists have long been using genetic tools to increase or maintain diversity of natural populations, genomic techniques could add extra benefit in AM that include selectively neutral and adaptive regions of the genome. In this review, we first propose a framework along with detailed procedures to aid collaboration among scientists, agencies, and local and regional managers during the decision-making process of genomics-guided AM. We then summarize the genomic approaches for applying AM, followed by a literature search of existing incorporation of genomics in AM across taxa. Our literature search initially identified 729 publications, but after filtering returned only 50 empirical studies that were either directly applied or considered genomics in AM related to climate change across taxa of plants, terrestrial animals, and aquatic animals; 42 studies were in plants. This demonstrated limited application of genomic methods in AM in organisms other than plants, so we provide further case studies as two examples to demonstrate the negative impact of climate change on non-model species and how genomics could be applied in AM. With the rapidly developing sequencing technology and accumulating genomic data, we expect to see more successful applications of genomics in AM, and more broadly, in the conservation of biodiversity.

Genomic survey of edible cockle ( Cerastoderma edule ) in the Northeast Atlantic: a baseline for sustainable management of its wild resources

Knowledge on correlations between environmental factors and genome divergence between populations of marine species is crucial for sustainable management of fisheries and wild populations. The edible cockle (Cerastoderma edule) is a marine bivalve distributed along the Northeast Atlantic coast of Europe and is an important resource from both commercial and ecological perspectives. We performed a population genomics screening using 2b‐RAD genotyping on 9309 SNPs localized in the cockle's genome on a sample of 536 specimens pertaining to 14 beds in the Northeast Atlantic Ocean to analyse the genetic structure with regard to environmental variables. Larval dispersal modelling considering species behaviour and interannual/interseasonal variation in ocean conditions was carried out as an essential background to which compare genetic information. Cockle populations in the Northeast Atlantic displayed low but significant geographical differentiation between populations (F_ST = 0.0240; p < 0.001), albeit not across generations. We identified 742 and 36 outlier SNPs related to divergent and balancing selection in all the geographical scenarios inspected, and sea temperature and salinity were the main environmental correlates suggested. Highly significant linkage disequilibrium was detected at specific genomic regions against the very low values observed across the whole genome. Two main genetic groups were identified, northwards and southwards of French Brittany. Larval dispersal modelling suggested a barrier for larval dispersal linked to the Ushant front that could explain these two genetic clusters. Further genetic subdivision was observed using outlier loci and considering larval advection. The northern group was divided into the Irish/Celtic Seas and the English Channel/North Sea, while the southern group was divided into three subgroups. This information represents the baseline for the management of cockles, designing conservation strategies, founding broodstock for depleted beds and producing suitable seed for aquaculture production.

Comparative cytogenetic patterns in Carangidae fishes in association with their distribution range

Carangidae are an important and widespreaded family of pelagic predatory fishes that inhabit reef regions or open ocean areas, some species occupying a vast circumglobal distribution. Cytogenetic comparisons among representatives of its different tribes help to understand the process of karyotype divergence in marine ecosystems due to the variable migratory ability of species. In this sense, conventional cytogenetic investigations (Giemsa staining, Ag-NORs, and C-banding), GC base-specific fluorochrome staining and FISH mapping of ribosomal DNAs were performed. Four species, Elagatisbipinnulata (Quoy et Gaimard, 1825) and Seriolarivoliana (Valenciennes, 1883) (Naucratini), with circumtropical distributions, Gnathanodonspeciosus (Forsskål, 1775) (Carangini), widely distributed in the tropical and subtropical waters of the Indian and Pacific oceans, and Trachinotuscarolinus (Linnaeus, 1766) (Trachinotini), distributed along the western Atlantic Ocean, were analyzed, thus encompassing representatives of three out its four tribes. All species have diploid chromosome number 2n = 48, with karyotypes composed mainly by acrocentric chromosomes (NF = 50-56). The 18S rDNA/Ag-NORs/GC+ and 5S rDNA loci were located on chromosomes likely homeologs. Karyotypes showed a pattern considered basal for the family or with small variations in their structures, apparently due to pericentric inversions. The migratory capacity of large pelagic swimmers, in large distribution areas, likely restricts the fixation of chromosome changes in Carangidae responsible for a low level of karyotype diversification.

Optimizing sampling design and sequencing strategy for the genomic analysis of quantitative traits in natural populations

Mapping the genes underlying ecologically relevant traits in natural populations is fundamental to develop a molecular understanding of species adaptation. Current sequencing technologies enable the characterization of a species' genetic diversity across the landscape or even over its whole range. The relevant capture of the genetic diversity across the landscape is critical for a successful genetic mapping of traits and there are no clear guidelines on how to achieve an optimal sampling and which sequencing strategy to implement. Here we determine, through simulation, the sampling scheme that maximizes the power to map the genetic basis of a complex trait in an outbreeding species across an idealized landscape and draw genomic predictions for the trait, comparing individual and pool sequencing strategies. Our results show that quantitative trait locus detection power and prediction accuracy are higher when more populations over the landscape are sampled and this is more cost-effectively done with pool sequencing than with individual sequencing. Additionally, we recommend sampling populations from areas of high genetic diversity. As progress in sequencing enables the integration of trait-based functional ecology into landscape genomics studies, these findings will guide study designs allowing direct measures of genetic effects in natural populations across the environment.

Genetic differentiation and signatures of local adaptation revealed by RADseq for a highly dispersive mud crab Scylla olivacea (Herbst, 1796) in the Sulu Sea

Connectivity of marine populations is shaped by complex interactions between biological and physical processes across the seascape. The influence of environmental features on the genetic structure of populations has key implications for the dynamics and persistence of populations, and an understanding of spatial scales and patterns of connectivity is crucial for management and conservation. This study employed a seascape genomics approach combining larval dispersal modeling and population genomic analysis using single nucleotide polymorphisms (SNPs) obtained from RADseq to examine environmental factors influencing patterns of genetic structure and connectivity for a highly dispersive mud crab Scylla olivacea (Herbst, 1796) in the Sulu Sea. Dispersal simulations reveal widespread but asymmetric larval dispersal influenced by persistent southward and westward surface circulation features in the Sulu Sea. Despite potential for widespread dispersal across the Sulu Sea, significant genetic differentiation was detected among eight populations based on 1,655 SNPs (FST  = 0.0057, p < .001) and a subset of 1,643 putatively neutral SNP markers (FST  = 0.0042, p < .001). Oceanography influences genetic structure, with redundancy analysis (RDA) indicating significant contribution of asymmetric ocean currents to neutral genetic variation ( R adj 2  = 0.133, p = .035). Genetic structure may also reflect demographic factors, with divergent populations characterized by low effective population sizes (N e < 50). Pronounced latitudinal genetic structure was recovered for loci putatively under selection (FST  = 0.2390, p < .001), significantly correlated with sea surface temperature variabilities during peak spawning months for S. olivacea ( R adj 2  = 0.692-0.763; p < .050), suggesting putative signatures of selection and local adaptation to thermal clines. While oceanography and dispersal ability likely shape patterns of gene flow and genetic structure of S. olivacea across the Sulu Sea, the impacts of genetic drift and natural selection influenced by sea surface temperature also appear as likely drivers of population genetic structure. This study contributes to the growing body of literature documenting population genetic structure and local adaptation for highly dispersive marine species, and provides information useful for spatial management of the fishery resource.

Charting a course for genetic diversity in the UN Decade of Ocean Science

The health of the world's oceans is intrinsically linked to the biodiversity of the ecosystems they sustain. The importance of protecting and maintaining ocean biodiversity has been affirmed through the setting of the UN Sustainable Development Goal 14 to conserve and sustainably use the ocean for society's continuing needs. The decade beginning 2021-2030 has additionally been declared as the UN Decade of Ocean Science for Sustainable Development. This program aims to maximize the benefits of ocean science to the management, conservation, and sustainable development of the marine environment by facilitating communication and cooperation at the science-policy interface. A central principle of the program is the conservation of species and ecosystem components of biodiversity. However, a significant omission from the draft version of the Decade of Ocean Science Implementation Plan is the acknowledgment of the importance of monitoring and maintaining genetic biodiversity within species. In this paper, we emphasize the importance of genetic diversity to adaptive capacity, evolutionary potential, community function, and resilience within populations, as well as highlighting some of the major threats to genetic diversity in the marine environment from direct human impacts and the effects of global climate change. We then highlight the significance of ocean genetic diversity to a diverse range of socioeconomic factors in the marine environment, including marine industries, welfare and leisure pursuits, coastal communities, and wider society. Genetic biodiversity in the ocean, and its monitoring and maintenance, is then discussed with respect to its integral role in the successful realization of the 2030 vision for the Decade of Ocean Science. Finally, we suggest how ocean genetic diversity might be better integrated into biodiversity management practices through the continued interaction between environmental managers and scientists, as well as through key leverage points in industry requirements for Blue Capital financing and social responsibility.

Restricted dispersal in a sea of gene flow

How far do marine larvae disperse in the ocean? Decades of population genetic studies have revealed generally low levels of genetic structure at large spatial scales (hundreds of kilometres). Yet this result, typically based on discrete sampling designs, does not necessarily imply extensive dispersal. Here, we adopt a continuous sampling strategy along 950 km of coast in the northwestern Mediterranean Sea to address this question in four species. In line with expectations, we observe weak genetic structure at a large spatial scale. Nevertheless, our continuous sampling strategy uncovers a pattern of isolation by distance at small spatial scales (few tens of kilometres) in two species. Individual-based simulations indicate that this signal is an expected signature of restricted dispersal. At the other extreme of the connectivity spectrum, two pairs of individuals that are closely related genetically were found more than 290 km apart, indicating long-distance dispersal. Such a combination of restricted dispersal with rare long-distance dispersal events is supported by a high-resolution biophysical model of larval dispersal in the study area, and we posit that it may be common in marine species. Our results bridge population genetic studies with direct dispersal studies and have implications for the design of marine reserve networks.

Genomic signatures of divergent selection are associated with social behaviour for spinner dolphin ecotypes

Understanding the genomic basis of adaptation is critical for understanding evolutionary processes and predicting how species will respond to environmental change. Spinner dolphins in the eastern tropical Pacific (ETP) present a unique system for studying adaptation. Within this large geographical region are four spinner dolphin ecotypes with weak neutral genetic divergence and no obvious barriers to gene flow, but strong spatial variation in morphology, behaviour and habitat. These ecotypes have large population sizes, which could reduce the effects of drift and facilitate selection. To identify genomic regions putatively under divergent selective pressures between ecotypes, we used genome scans with 8994 RADseq single nucleotide polymorphisms (SNPs) to identify population differentiation outliers and genotype-environment association outliers. Gene ontology enrichment analyses indicated that outlier SNPs from both types of analyses were associated with multiple genes involved in social behaviour and hippocampus development, including 15 genes associated with the human social disorder autism. Evidence for divergent selection on social behaviour is supported by previous evidence that these spinner dolphin ecotypes differ in mating systems and associated social behaviours. In particular, three of the ETP ecotypes probably have a polygynous mating system characterized by strong premating competition among males, whereas the fourth ecotype probably has a polygynandrous mating system characterized by strong postmating competition such as sperm competition. Our results provide evidence that selection for social behaviour may be an evolutionary force driving diversification of spinner dolphins in the ETP, potentially as a result of divergent sexual selection associated with different mating systems. Future studies should further investigate the potential adaptive role of the candidate genes identified here, and could probably find further signatures of selection using whole genome sequence data.

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.

Population structure in a continuously distributed coastal marine species, the harbor porpoise, based on microhaplotypes derived from poor-quality samples

Harbor porpoise in the North Pacific are found in coastal waters from southern California to Japan, but population structure is poorly known outside of a few local areas. We used multiplexed amplicon sequencing of 292 loci and genotyped clusters of single nucleotide polymoirphisms as microhaplotypes (N = 271 samples) in addition to mitochondrial (mtDNA) sequence data (N = 413 samples) to examine the genetic structure from samples collected along the Pacific coast and inland waterways from California to southern British Columbia. We confirmed an overall pattern of strong isolation-by-distance, suggesting that individual dispersal is restricted. We also found evidence of regions where genetic differences are larger than expected based on geographical distance alone, implying current or historical barriers to gene flow. In particular, the southernmost population in California is genetically distinct (F_ST  = 0.02 [microhaplotypes]; 0.31 [mtDNA]), with both reduced genetic variability and high frequency of an otherwise rare mtDNA haplotype. At the northern end of our study range, we found significant genetic differentiation of samples from the Strait of Georgia, previously identified as a potential biogeographical boundary or secondary contact zone between harbor porpoise populations. Association of microhaplotypes with remotely sensed environmental variables indicated potential local adaptation, especially at the southern end of the species' range. These results inform conservation and management for this nearshore species, illustrate the value of genomic methods for detecting patterns of genetic structure within a continuously distributed marine species, and highlight the power of microhaplotype genotyping for detecting genetic structure in harbor porpoises despite reliance on poor-quality samples.

Reconciling seascape genetics and fisheries science in three codistributed flatfishes

Uncertainty hampers innovative mixed-fisheries management by the scales at which connectivity dynamics are relevant to management objectives. The spatial scale of sustainable stock management is species-specific and depends on ecology, life history and population connectivity. One valuable approach to understand these spatial scales is to determine to what extent population genetic structure correlates with the oceanographic environment. Here, we compare the level of genetic connectivity in three codistributed and commercially exploited demersal flatfish species living in the North East Atlantic Ocean. Population genetic structure was analysed based on 14, 14 and 10 neutral DNA microsatellite markers for turbot, brill and sole, respectively. We then used redundancy analysis (RDA) to attribute the genetic variation to spatial (geographical location), temporal (sampling year) and oceanographic (water column characteristics) components. The genetic structure of turbot was composed of three clusters and correlated with variation in the depth of the pycnocline, in addition to spatial factors. The genetic structure of brill was homogenous, but correlated with average annual stratification and spatial factors. In sole, the genetic structure was composed of three clusters, but was only linked to a temporal factor. We explored whether the management of data poor commercial fisheries, such as in brill and turbot, might benefit from population-specific information. We conclude that the management of fish stocks has to consider species-specific genetic structures and may benefit from the documentation of the genetic seascape and life-history traits.

Population genomics and history of speciation reveal fishery management gaps in two related redfish species (Sebastes mentella and Sebastes fasciatus)

Understanding the processes shaping population structure and reproductive isolation of marine organisms can improve their management and conservation. Using genomic markers combined with estimation of individual ancestries, assignment tests, spatial ecology, and demographic modeling, we (i) characterized the contemporary population structure, (ii) assessed the influence of space, fishing depth, and sampling years on contemporary distribution, and (iii) reconstructed the speciation history of two cryptic redfish species, Sebastes mentella and S. fasciatus. We genotyped 860 individuals in the Northwest Atlantic Ocean using 24,603 filtered single nucleotide polymorphisms (SNPs). Our results confirmed the clear genetic distinctiveness of the two species and identified three ecotypes within S. mentella and five populations in S. fasciatus. Multivariate analyses highlighted the influence of spatial distribution and depth on the overall genomic variation, while demographic modeling revealed that secondary contact models best explained inter- and intragenomic divergence. These species, ecotypes, and populations can be considered as a rare and wide continuum of genomic divergence in the marine environment. This acquired knowledge pertaining to the evolutionary processes driving population divergence and reproductive isolation will help optimizing the assessment of demographic units and possibly to refine fishery management units.

Coral cover surveys corroborate predictions on reef adaptive potential to thermal stress

As anomalous heat waves are causing the widespread decline of coral reefs worldwide, there is an urgent need to identify coral populations tolerant to thermal stress. Heat stress adaptive potential is the degree of tolerance expected from evolutionary processes and, for a given reef, depends on the arrival of propagules from reefs exposed to recurrent thermal stress. For this reason, assessing spatial patterns of thermal adaptation and reef connectivity is of paramount importance to inform conservation strategies. In this work, we applied a seascape genomics framework to characterize the spatial patterns of thermal adaptation and connectivity for coral reefs of New Caledonia (Southern Pacific). In this approach, remote sensing of seascape conditions was combined with genomic data from three coral species. For every reef of the region, we computed a probability of heat stress adaptation, and two indices forecasting inbound and outbound connectivity. We then compared our indicators to field survey data, and observed that decrease of coral cover after heat stress was lower at reefs predicted with high probability of adaptation and inbound connectivity. Last, we discussed how these indicators can be used to inform local conservation strategies and preserve the adaptive potential of New Caledonian reefs.

Extreme seascape drives local recruitment and genetic divergence in brooding and spawning corals in remote north-west Australia

Management strategies designed to conserve coral reefs threatened by climate change need to incorporate knowledge of the spatial distribution of inter- and intra-specific genetic diversity. We characterized patterns of genetic diversity and connectivity using single nucleotide polymorphisms (SNPs) in two reef-building corals to explore the eco-evolutionary processes that sustain populations in north-west Australia. Our sampling focused on the unique reefs of the Kimberley; we collected the broadcast spawning coral Acropora aspera (n = 534) and the brooding coral Isopora brueggemanni (n = 612) across inter-archipelago (tens to hundreds of kilometres), inter-reef (kilometres to tens of kilometres) and within-reef (tens of metres to a few kilometres) scales. Initial analysis of A. aspera identified four highly divergent lineages that were co-occurring but morphologically similar. Subsequent population analyses focused on the most abundant and widespread lineage, Acropora asp-c. Although the overall level of geographic subdivision was greater in the brooder than in the spawner, fundamental similarities in patterns of genetic structure were evident. Most notably, limits to gene flow were observed at scales <35 kilometres. Further, we observed four discrete clusters and a semi-permeable barrier to dispersal that were geographically consistent between species. Finally, sites experiencing bigger tides were more connected to the metapopulation and had greater gene diversity than those experiencing smaller tides. Our data indicate that the inshore reefs of the Kimberley are genetically isolated from neighbouring oceanic bioregions, but occasional dispersal between inshore archipelagos is important for the redistribution of evolutionarily important genetic diversity. Additionally, these results suggest that networks of marine reserves that effectively protect reefs from local pressures should be spaced within a few tens of kilometres to conserve the existing patterns of demographic and genetic connectivity.

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.

Seascape genomics as a new tool to empower coral reef conservation strategies: An example on north-western Pacific Acropora digitifera

Coral reefs are suffering a major decline due to the environmental constraints imposed by climate change. Over the last 20 years, three major coral bleaching events occurred in concomitance with anomalous heatwaves, provoking a severe loss of coral cover worldwide. The conservation strategies for preserving reefs, as they are implemented now, cannot cope with global climatic shifts. Consequently, researchers are advocating for preservation networks to be set-up to reinforce coral adaptive potential. However, the main obstacle to this implementation is that studies on coral adaption are usually hard to generalize at the scale of a reef system. Here, we study the relationships between genotype frequencies and environmental characteristics of the sea (seascape genomics), in combination with connectivity analysis, to investigate the adaptive potential of a flagship coral species of the Ryukyu Archipelago (Japan). By associating genotype frequencies with descriptors of historical environmental conditions, we discovered six genomic regions hosting polymorphisms that might promote resistance against heat stress. Remarkably, annotations of genes in these regions were consistent with molecular roles associated with heat responses. Furthermore, we combined information on genetic and spatial distances between reefs to predict connectivity at a regional scale. The combination of these results portrayed the adaptive potential of this population: we were able to identify reefs carrying potential heat stress adapted genotypes and to understand how they disperse to neighbouring reefs. This information was summarized by objective, quantifiable and mappable indices covering the whole region, which can be extremely useful for future prioritization of reefs in conservation planning. This framework is transferable to any coral species on any reef system and therefore represents a valuable tool for empowering preservation efforts dedicated to the protection of coral reefs in warming oceans.

New insights into the population genetics of partially clonal organisms: When seagrass data meet theoretical expectations

Seagrass meadows are among the most important coastal ecosystems in terms of both spatial extent and ecosystem services, but they are also declining worldwide. Understanding the drivers of seagrass meadow dynamics is essential for designing sound management, conservation and restoration strategies. However, poor knowledge of the effect of clonality on the population genetics of natural populations severely limits our understanding of the dynamics and connectivity of meadows. Recent modelling approaches have described the expected distributions of genotypic and genetic descriptors under increasing clonal rates, which may help us better understand and interpret population genetics data obtained for partial asexuals. Here, in the light of these recent theoretical developments, we revisited population genetics data for 165 meadows of four seagrass species. Contrasting shoot lifespan and rhizome turnover led to the prediction that the influence of asexual reproduction would increase along a gradient from Zostera noltii to Zostera marina, Cymodocea nodosa and Posidonia oceanica, with increasing departure from Hardy-Weinberg equilibrium (F_is ), mostly towards heterozygote excess, and decreasing genotypic richness (R). This meta-analysis provides a nested validation of this hypothesis at both the species and meadow scales through a significant relationship between F_is and R within each species. By empirically demonstrating the theoretical expectations derived from recent modelling approaches, this work calls for the use of Hardy-Weinberg equilibrium (F_is ) rather than only the strongly sampling-sensitive R to assess the importance of clonal reproduction (c), at least when the impact of selfing on F_is can be neglected. The results also emphasize the need to revise our appraisal of the extent of clonality and its influence on the dynamics, connectivity and evolutionary trajectory of partial asexuals in general, including in seagrass meadows, to develop the most accurate management strategies.

Signatures of local adaptation in the spatial genetic structure of the ascidian Pyura chilensis along the southeast Pacific coast

The highly heterogeneous Humboldt Current System (HCS) and the 30°S transition zone on the southeast Pacific coast, represent an ideal scenario to test the influence of the environment on the spatial genomic structure in marine near-shore benthic organisms. In this study, we used seascape genomic tools to evaluate the genetic structure of the commercially important ascidian Pyura chilensis, a species that exhibits a low larval transport potential but high anthropogenic dispersal. A recent study in this species recorded significant genetic differentiation across a transition zone around 30°S in putatively adaptive SNPs, but not in neutral ones, suggesting an important role of environmental heterogeneity in driving genetic structure. Here, we aim to understand genomic-oceanographic associations in P. chilensis along the Southeastern Pacific coast using two combined seascape genomic approaches. Using 149 individuals from five locations along the HCS, a total of 2,902 SNPs were obtained by Genotyping-By-Sequencing, of which 29–585 were putatively adaptive loci, depending on the method used for detection. In adaptive loci, spatial genetic structure was better correlated with environmental differences along the study area (mainly to Sea Surface Temperature, upwelling-associated variables and productivity) than to the geographic distance between sites. Additionally, results consistently showed the presence of two groups, located north and south of 30°S, which suggest that local adaptation processes seem to allow the maintenance of genomic differentiation and the spatial genomic structure of the species across the 30°S biogeographic transition zone of the Humboldt Current System, overriding the homogenizing effects of gene flow.

Investigating population-scale allelic differential expression in wild populations of Oithona similis (Cyclopoida, Claus, 1866)

Acclimation allowed by variation in gene or allele expression in natural populations is increasingly understood as a decisive mechanism, as much as adaptation, for species evolution. However, for small eukaryotic organisms, as species from zooplankton, classical methods face numerous challenges. Here, we propose the concept of allelic differential expression at the population-scale (psADE) to investigate the variation in allele expression in natural populations. We developed a novel approach to detect psADE based on metagenomic and metatranscriptomic data from environmental samples. This approach was applied on the widespread marine copepod, Oithona similis, by combining samples collected during the Tara Oceans expedition (2009-2013) and de novo transcriptome assemblies. Among a total of 25,768 single nucleotide variants (SNVs) of O. similis, 572 (2.2%) were affected by psADE in at least one population (FDR < 0.05). The distribution of SNVs under psADE in different populations is significantly shaped by population genomic differentiation (Pearson r = 0.87, p = 5.6 × 10-30), supporting a partial genetic control of psADE. Moreover, a significant amount of SNVs (0.6%) were under both selection and psADE (p < .05), supporting the hypothesis that natural selection and psADE tends to impact common loci. Population-scale allelic differential expression offers new insights into the gene regulation control in populations and its link with natural selection.

Genomic Signatures of Local Adaptation in Clam Shrimp (Eulimnadia texana) from Natural Vernal Pools

Vernal pools are unique in their isolation and the strong selection acting on their resident species. Vernal pool clam shrimp (Eulimnadia texana) are a promising model due to ease of culturing, short generation time, small genomes, and obligate desiccated diapaused eggs. Clam shrimp are also androdioecious (sexes include males and hermaphrodites), and here we use population-scaled recombination rates to support the hypothesis that the heterogametic sex is recombination free in these shrimp. We collected short-read sequence data from pooled samples from different vernal pools to gain insights into local adaptation. We identify genomic regions in which some populations have allele frequencies that differ significantly from the metapopulation. BayPass (Gautier M. 2015. Genome-wide scan for adaptive divergence and association with population-specific covariates. Genetics 201(4):1555-1579.) detected 19 such genomic regions showing an excess of population subdivision. These regions on average are 550 bp in size and had 2.5 genes within 5 kb of them. Genes located near these regions are involved in Malpighian tubule function and osmoregulation, an essential function in vernal pools. It is likely that salinity profiles vary between pools and over time, and variants at these genes are adapted to local salinity conditions.

Asymmetrical gene flow in five co-distributed syngnathids explained by ocean currents and rafting propensity

Ocean circulation driving macro-algal rafting is believed to serve as an important mode of dispersal for many marine organisms, leading to predictions on population-level genetic connectivity and the directionality of effective dispersal. Here, we use genome-wide single nucleotide polymorphism data to investigate whether gene flow directionality in two seahorses (Hippocampus) and three pipefishes (Syngnathus) follows the predominant ocean circulation patterns in the Gulf of Mexico and northwestern Atlantic. In addition, we explore whether gene flow magnitudes are predicted by traits related to active dispersal ability and habitat preference. We inferred demographic histories of these co-distributed syngnathid species, and coalescent model-based estimates indicate that gene flow directionality is in agreement with ocean circulation data that predicts eastward and northward macro-algal transport. However, the magnitude to which ocean currents influence this pattern appears strongly dependent on the species-specific traits related to rafting propensity and habitat preferences. Higher levels of gene flow and stronger directionality are observed in Hippocampus erectus, Syngnathus floridae and Syngnathus louisianae, which closely associated with the pelagic macro-algae Sargassum spp., compared to Hippocampus zosterae and the Syngnathus scovelli/Syngnathus fuscus sister-species pair, which prefer near shore habitats and are weakly associated with pelagic Sargassum. This study highlights how the combination of population genomic inference together with ocean circulation data can help explain patterns of population structure and diversity in marine ecosystems.

Global determinants of freshwater and marine fish genetic diversity

Genetic diversity is estimated to be declining faster than species diversity under escalating threats, but its spatial distribution remains poorly documented at the global scale. Theory predicts that similar processes should foster congruent spatial patterns of genetic and species diversity, but empirical studies are scarce. Using a mined database of 50,588 georeferenced mitochondrial DNA barcode sequences (COI) for 3,815 marine and 1,611 freshwater fish species respectively, we examined the correlation between genetic diversity and species diversity and their global distributions in relation to climate and geography. Genetic diversity showed a clear spatial organisation, but a weak association with species diversity for both marine and freshwater species. We found a predominantly positive relationship between genetic diversity and sea surface temperature for marine species. Genetic diversity of freshwater species varied primarily across the regional basins and was negatively correlated with average river slope. The detection of genetic diversity patterns suggests that conservation measures should consider mismatching spatial signals across multiple facets of biodiversity.

Biogeographic patterns of genetic diversity are poorly documented, especially for fish species. Here the authors show that (mitochondrial) genetic diversity has global spatial organization patterns with different environmental drivers for marine and freshwater fishes, where genetic diversity is only partly congruent with species richness.

Ocean currents and the population genetic signature of fish migrations

Animal migrations are a fascinating and global phenomenon, yet they are often difficult to study and sometimes poorly understood. Here, we build on classic ecological theory by hypothesizing that some enigmatic spawning migrations across coastal marine habitats can be inferred from the population genetic signature of larval dispersal by ocean currents. We test this assumption by integrating spatially realistic simulations of alternative spawning migration routes, associated patterns of larval dispersal, and associated variation in the population genetic structure of eastern Australian sea mullet (Mugil cephalus). We then use simulation results to assess the implications of alternative spawning destinations for larval replenishment, and we contrast simulated against measured population genetic variation. Both analyses suggest that the spawning migrations of M. cephalus in eastern Australia are likely to be localized (approximately 100 km along the shore), and that spawning is likely to occur in inshore waters. Our conclusions are supported by multiple lines of evidence available through independent studies, but they challenge the more traditional assumption of a single, long-distance migration event with subsequent offshore spawning in the East Australian Current. More generally, our study operationalizes classic theory on the relationship between fish migrations, ocean currents, and reproductive success. However, rather than confirming the traditionally assumed adaptation of migratory behavior to dominant ocean current flow, our findings support the concept of a genetically measurable link between fish migrations and local oceanographic conditions, specifically water temperature and coastal retention of larvae. We believe that future studies using similar approaches for high resolution and spatially realistic ecological-genetic scenario testing can help rapidly advance our understanding of key ecological processes in many other marine species.

Sex matters: Otolith shape and genomic variation in deacon rockfish (Sebastes diaconus)

Little is known about intraspecific variation within the deacon rockfish (Sebastes diaconus), a recently described species found in the northeast Pacific Ocean. We investigated population structure among fish sampled from two nearshore reefs (Siletz Reef and Seal Rock) and one offshore site (Stonewall Bank) within a <50-km2 area off the Oregon coast. Fish from the three sample sites exhibited small but statistically significant differences based on genetic variation at >15,000 neutral loci, whether analyzed independently or classified into nearshore and offshore groups. Male and females were readily distinguished using genetic data and 92 outlier loci were associated with sex, potentially indicating differential selection between males and females. Morphometric results indicated that there was significant secondary sexual dimorphism in otolith shape, but further sampling is required to disentangle potential confounding influence of age. This study is the first step toward understanding intraspecific variation within the deacon rockfish and the potential management implications. Since differentiation among the three sample sites was small, we consider the results to be suggestive of a single stock. However, future studies should evaluate how the stock is affected by differences in sex, age, and gene flow between the nearshore and offshore environments.

Sampling strategy optimization to increase statistical power in landscape genomics: A simulation-based approach

An increasing number of studies are using landscape genomics to investigate local adaptation in wild and domestic populations. Implementation of this approach requires the sampling phase to consider the complexity of environmental settings and the burden of logistical constraints. These important aspects are often underestimated in the literature dedicated to sampling strategies. In this study, we computed simulated genomic data sets to run against actual environmental data in order to trial landscape genomics experiments under distinct sampling strategies. These strategies differed by design approach (to enhance environmental and/or geographical representativeness at study sites), number of sampling locations and sample sizes. We then evaluated how these elements affected statistical performances (power and false discoveries) under two antithetical demographic scenarios. Our results highlight the importance of selecting an appropriate sample size, which should be modified based on the demographic characteristics of the studied population. For species with limited dispersal, sample sizes above 200 units are generally sufficient to detect most adaptive signals, while in random mating populations this threshold should be increased to 400 units. Furthermore, we describe a design approach that maximizes both environmental and geographical representativeness of sampling sites and show how it systematically outperforms random or regular sampling schemes. Finally, we show that although having more sampling locations (between 40 and 50 sites) increase statistical power and reduce false discovery rate, similar results can be achieved with a moderate number of sites (20 sites). Overall, this study provides valuable guidelines for optimizing sampling strategies for landscape genomics experiments.

Conservation genetics of elasmobranchs of the Mexican Pacific Coast, trends and perspectives

One of the most critical threats to biodiversity is the high extinction rate driven by human activities. Reducing extinction rates requires the implementation of conservation programmes based on robust scientific data. Elasmobranchs are important ecological components of the ocean, and several species sustain substantial economic activities. Unfortunately, elasmobranchs are one of the most threatened and understudied animal taxa. The Mexican Pacific Coast (MPC) is a region with high elasmobranch diversity and is the seat of major elasmobranch fisheries. But it is also a developing region with several conservation and management challenges which require national and international attention. Here, we review the conservation genetics literature of elasmobranchs from the MPC. We present a synthesis of the works using samples from the region and emphasize the main gaps and biases in these data. In addition, we discuss the benefits and challenges of generating genomic information to improve the management and conservation of an elasmobranch biodiversity hotspot in a developing country. We found 47 elasmobranch genetic articles that cover <30% of the elasmobranch diversity in the region. These studies mainly used mitochondrial DNA sequences to analyse the genetic structure of commercially important and abundant species of the order Carcharhiniformes. Some of these papers also assessed mating systems, demographic parameters, and taxonomic uncertainties, all of which are important topics for efficient management decisions. In terms of conservation genetics, elasmobranchs from the MPC remain understudied. However, high-throughput sequencing technologies have increased the power and accessibility of genomic tools, even in developing countries such as Mexico. The tools described here provide information relevant for biodiversity conservation. Therefore, we strongly suggest that investment in genomic research will assist implementation of efficient management strategies. In time, this will reduce the extinction risk of the unique elasmobranch biodiversity from the MPC.