Marine populations are typically characterized by weak genetic differentiation due to the potential for long‐distance dispersal favouring high levels of gene flow. However, strong directional ...advection of water masses or retentive hydrodynamic forces can influence the degree of genetic exchange among marine populations. To determine the oceanographic drivers of genetic structure in a highly dispersive marine invertebrate, the giant California sea cucumber (Parastichopus californicus), we first tested for the presence of genetic discontinuities along the coast of North America in the northeastern Pacific Ocean. Then, we tested two hypotheses regarding spatial processes influencing population structure: (i) isolation by distance (IBD: genetic structure is explained by geographic distance) and (ii) isolation by resistance (IBR: genetic structure is driven by ocean circulation). Using RADseq, we genotyped 717 individuals from 24 sampling locations across 2,719 neutral SNPs to assess the degree of population differentiation and integrated estimates of genetic variation with inferred connectivity probabilities from a biophysical model of larval dispersal mediated by ocean currents. We identified two clusters separating north and south regions, as well as significant, albeit weak, substructure within regions (FST = 0.002, p = .001). After modelling the asymmetric nature of ocean currents, we demonstrated that local oceanography (IBR) was a better predictor of genetic variation (R2 = .49) than geographic distance (IBD) (R2 = .18), and directional processes played an important role in shaping fine‐scale structure. Our study contributes to the growing body of literature identifying significant population structure in marine systems and has important implications for the spatial management of P. californicus and other exploited marine species.
Habitat loss and fragmentation have important consequences for population persistence by reducing dispersal, connectivity, and population size. Characterizing patterns of population genetic structure ...and gene flow, and testing for the genetic signature of population declines and inbreeding can provide key insights into conservation for at-risk species. We use 12 microsatellite DNA markers to quantify the spatial structuring across two regional populations of the threatened eastern hog-nosed snake (Heterodon platirhinos) in Ontario, Canada. We further test for inbreeding and recent, rapid population declines, for the latter using both a ‘heterozygosity excess’ approach and Approximate Bayesian computation. Spatial Bayesian assignment tests revealed four distinct genetic clusters that were not apparent from species occurrence data, suggesting that the current designation of a single management unit for this species is invalid. We found some evidence for inbreeding in all studied locales. However, heterozygosity excess tests provided only weak evidence for recent bottlenecks and varied according to mutational model specified. In contrast, ABC provided strong support for a recent and rapid reduction in population size within one sample locale (Wasaga Beach), where individuals are largely confined to a small area circumscribed by expanding urban areas. For three other locales investigated, we found equal support for either recent rapid or gradual population declines. Immediate action is required for Wasaga Beach if this population is to persist.
Interactions between environmental factors and complex life‐history characteristics of marine organisms produce the genetic diversity and structure observed within species. Our main goal was to test ...for genetic differentiation among eastern oyster populations from the coastal region of Canadian Maritimes against expected genetic homogeneity caused by historical events, taking into account spatial and environmental (temperature, salinity, turbidity) variation. This was achieved by genotyping 486 individuals originating from 13 locations using RADSeq. A total of 11,321 filtered SNPs were used in a combination of population genomics and environmental association analyses. We revealed significant neutral genetic differentiation (mean FST = 0.009) between sampling locations, and the occurrence of six major genetic clusters within the studied system. Redundancy analyses (RDAs) revealed that spatial and environmental variables explained 3.1% and 4.9% of the neutral genetic variation and 38.6% and 12.2% of the putatively adaptive genetic variation, respectively. These results indicate that these environmental factors play a role in the distribution of both neutral and putatively adaptive genetic diversity in the system. Moreover, polygenic selection was suggested by genotype–environment association analysis and significant correlations between additive polygenic scores and temperature and salinity. We discuss our results in the context of their conservation and management implications for the eastern oyster.
Conservation units (CUs) are important tools for supporting the implementation of standardized management practices for exploited species. Following the adoption of the Wild Salmon Policy in Canada, ...CUs were defined for Pacific salmon based on characteristics related to ecotype, life history and genetic variation using microsatellite markers as indirect measures of local adaptation. Genomic data sets have the potential to improve the definition of CUs by reducing variance around estimates of population genetic parameters, thereby increasing the power to detect more subtle patterns of population genetic structure and by providing an opportunity to incorporate adaptive information more directly with the identification of variants putatively under selection. We used one of the largest genomic data sets recently published for a nonmodel species, comprising 5662 individual Coho salmon (Oncorhynchus kisutch) from 149 sampling locations and a total of 24,542 high‐quality SNPs obtained using genotyping‐by‐sequencing and mapped to the Coho salmon reference genome to (1) evaluate the current delineation of CUs for Coho in Canada and (2) compare patterns of population structure observed using neutral and outlier loci from genotype–environment association analyses to determine whether separate CUs that capture adaptive diversity are needed. Our results reflected CU boundaries on the whole, with the majority of sampling locations managed in the same CU clustering together within genetic groups. However, additional groups that are not currently represented by CUs were also uncovered. We observed considerable overlap in the genetic clusters identified using neutral or candidate loci, indicating a general congruence in patterns of genetic variation driven by local adaptation and gene flow in this species. Consequently, we suggest that the current CU boundaries for Coho salmon are largely well‐suited for meeting the Canadian Wild Salmon Policy's objective of defining biologically distinct groups, but we highlight specific areas where CU boundaries may be refined.
Abstract
Dense single nucleotide polymorphism (SNP) arrays are essential tools for rapid high-throughput genotyping for many genetic analyses, including genomic selection and high-resolution ...population genomic assessments. We present a high-density (200 K) SNP array developed for the Eastern oyster (Crassostrea virginica), which is a species of significant aquaculture production and restoration efforts throughout its native range. SNP discovery was performed using low-coverage whole-genome sequencing of 435 F1 oysters from families from 11 founder populations in New Brunswick, Canada. An Affymetrix Axiom Custom array was created with 219,447 SNPs meeting stringent selection criteria and validated by genotyping more than 4,000 oysters across 2 generations. In total, 144,570 SNPs had a call rate >90%, most of which (96%) were polymorphic and were distributed across the Eastern oyster reference genome, with similar levels of genetic diversity observed in both generations. Linkage disequilibrium was low (maximum r2 ∼0.32) and decayed moderately with increasing distance between SNP pairs. Taking advantage of our intergenerational data set, we quantified Mendelian inheritance errors to validate SNP selection. Although most of SNPs exhibited low Mendelian inheritance error rates overall, with 72% of called SNPs having an error rate of <1%, many loci had elevated Mendelian inheritance error rates, potentially indicating the presence of null alleles. This SNP panel provides a necessary tool to enable routine application of genomic approaches, including genomic selection, in C. virginica selective breeding programs. As demand for production increases, this resource will be essential for accelerating production and sustaining the Canadian oyster aquaculture industry.
The complexity of climate change impacts on ecological processes necessitates flexible and adaptive conservation strategies that cross traditional disciplines. Current strategies involving protected ...areas are predominantly fixed in space, and may on their own be inadequate under climate change. Here, we propose a novel approach to climate adaptation that combines permanent protected areas with temporary conservation areas to create flexible networks. Previous work has tended to consider permanent and dynamic protection as separate actions, but their integration could draw on the strengths of both approaches to improve biodiversity conservation and help manage for ecological uncertainty in the coming decades. As there are often time lags in the establishment of new permanent protected areas, the inclusion of dynamic conservation areas within permanent networks could provide critical transient protection to mitigate land-use changes and biodiversity redistributions. This integrated approach may be particularly useful in highly human-modified and fragmented landscapes where areas of conservation value are limited and long-term place-based protection is unfeasible. To determine when such an approach may be feasible, we propose the use of a decision framework. Under certain scenarios, these coupled networks have the potential to increase spatio-temporal network connectivity and help maintain biodiversity and ecological processes under climate change. Implementing these networks would require multidisciplinary scientific evidence, new policies, creative funding solutions, and broader acceptance of a dynamic approach to biodiversity conservation.
Coastal marine environments are subject to a variety of anthropogenic pressures that can negatively impact habitats and the biodiversity they harbor. Conservation actions such as marine protected ...areas, marine reserves, and other effective area-based conservation measures, are pivotal tools for protecting coastal biodiversity. However, to be effective, conservation area networks must be planned through a systematic conservation planning (SCP) approach. Recently, such approaches have begun to orient their goals toward the conservation of different biodiversity facets and to integrate different types of data. In this review, we illustrate how genetic data and molecular techniques can bring useful knowledge for SCP approaches that are both more comprehensive (sampling the full range of biodiversity) and more adequate (ensuring the long-term persistence of biodiversity). With an emphasis on coastal organisms and habitats, we focus on phylogenetic analysis, the estimation of neutral and adaptive intraspecific genetic diversity at different spatial levels (alpha, beta, and gamma), the study of connectivity and dispersal, and the information obtainable from environmental DNA techniques. For each of these applications, we discuss the benefits of its integration into SCP for coastal systems, its strengths and weaknesses, and the aspects yet to be developed.
The eastern hog-nosed snake (Heterodon platirhinos) is a federally listed species of conservation concern in Canada. Microsatellite markers are desired for research on the spatial structure, degree ...of isolation and neutral genetic status of the two regional isolated populations in the province of Ontario. We developed 14 microsatellite loci and genotyped 27 individuals from one of the two regional populations. The total number of alleles and the level of expected heterozygosity ranged from 4 to 8 and from 0.35 to 0.84, respectively. These microsatellite loci will be useful for assessing neutral genetic diversity and demarcating functionally independent groups for effective management practices.
As climate change threatens species' persistence, predicting the potential for species to adapt to rapidly changing environments is imperative for the development of effective conservation ...strategies. Eco-evolutionary individual-based models (IBMs) can be useful tools for achieving this objective. We performed a literature review to identify studies that apply these tools in marine systems. Our survey suggested that this is an emerging area of research fuelled in part by developments in modelling frameworks that allow simulation of increasingly complex ecological, genetic and demographic processes. The studies we identified illustrate the promise of this approach and advance our understanding of the capacity for adaptation to outpace climate change. These studies also identify limitations of current models and opportunities for further development. We discuss three main topics that emerged across studies: (i) effects of genetic architecture and non-genetic responses on adaptive potential; (ii) capacity for gene flow to facilitate rapid adaptation; and (iii) impacts of multiple stressors on persistence. Finally, we demonstrate the approach using simple simulations and provide a framework for users to explore eco-evolutionary IBMs as tools for understanding adaptation in changing seas.
The availability of genomic data for an increasing number of species makes it possible to incorporate evolutionary processes into conservation plans. Recent studies show how genetic data can inform ...spatial conservation prioritization (SCP), but they focus on metrics of diversity and distinctness derived primarily from neutral genetic data sets. Identifying adaptive genetic markers can provide important information regarding the capacity for populations to adapt to environmental change. Yet, the effect of including metrics based on adaptive genomic data into SCP in comparison to more widely used neutral genetic metrics has not been explored. We used existing genomic data on a commercially exploited species, the giant California sea cucumber (Parastichopus californicus), to perform SCP for the coastal region of British Columbia (BC), Canada. Using a RAD‐seq data set for 717 P. californicus individuals across 24 sampling locations, we identified putatively adaptive (i.e., candidate) single nucleotide polymorphisms (SNPs) based on genotype–environment associations with seafloor temperature. We calculated various metrics for both neutral and candidate SNPs and compared SCP outcomes with independent metrics and combinations of metrics. Priority areas varied depending on whether neutral or candidate SNPs were used and on the specific metric used. For example, targeting sites with a high frequency of warm‐temperature‐associated alleles to support persistence under future warming prioritized areas in the southern coastal region. In contrast, targeting sites with high expected heterozygosity at candidate loci to support persistence under future environmental uncertainty prioritized areas in the north. When combining metrics, all scenarios generated intermediate solutions, protecting sites that span latitudinal and thermal gradients. Our results demonstrate that distinguishing between neutral and adaptive markers can affect conservation solutions and emphasize the importance of defining objectives when choosing among various genomic metrics for SCP.
Incorporación de Datos Genómicos Putativamente Neutros y Adaptativos dentro de la Planeación de la Conservación Marina
Resumen
La disponibilidad de los datos genómicos para un número creciente de especies posibilita la incorporación de los procesos evolutivos dentro de los planes de conservación. Los estudios recientes muestran cómo los datos genéticos pueden informar a la priorización de la conservación espacial (PCE) pero tienden a enfocarse más en las medidas de la diversidad y la distinción derivadas principalmente de los conjuntos de datos genéticos neutrales. La identificación de los marcadores genéticos adaptativos puede proporcionar información importante con respecto a la capacidad de las poblaciones para adaptarse al cambio ambiental. Aun así, no se ha explorado el efecto de la inclusión de las medidas basadas en los datos genéticos adaptativos dentro de la PCE y cómo se comparan con las medidas genéticas neutrales de uso más amplio. Usamos datos genómicos existentes sobre una especie de explotación comercial, el pepino de mar gigante de California (Parastichopus californicus), para realizar la PCE para la región costera de la Columbia Británica (BC) en Canadá. Usamos un conjunto de datos RAD‐seq para 717 individuos de la especie P. californicus en 24 localidades de muestreo para identificar los polimorfismos de un solo nucleótido (PSNs) putativamente adaptativos (es decir, candidatos) con base en las asociaciones genotipo‐ambiente manifestadas con la temperatura del fondo marino. Calculamos varias medidas para los PSNs neutrales y los PSNs candidatos y comparamos los resultados de la PCE con medidas independientes y con combinaciones de medidas. Las áreas prioritarias variaron dependiendo de si se usaron los SNP neutrales o los candidatos y de la medida específica que se utilizó. Por ejemplo, enfocarse en sitios con una frecuencia alta de alelos asociados con agua cálida para fortalecer la persistencia frente al futuro calentamiento prioriza las áreas en la región del sur de la costa. Al contrario, enfocarse en sitios con una alta heterocigosidad esperada en los loci de los candidatos para fortalecer la persistencia frente a la incertidumbre ambiental prioriza las áreas en la parte norte de la costa. Cuando combinamos las medidas, todos los escenarios generaron soluciones intermedias, protegiendo así los sitios que abarcan gradientes latitudinales y gradientes térmicos. Nuestros resultados demuestran que la distinción entre los marcadores neutrales y los adaptativos puede afectar las soluciones de conservación y también enfatizan la importancia de la definición de los objetivos cuando se elige entre varias medidas genómicas para la PCE.
Article impact statement: Genomic metrics derived from neutral and adaptive genomic data result in different spatial conservation prioritization solutions.
PDF
