Spatiotemporal SNP analysis reveals pronounced biocomplexity at the northern range margin of Atlantic cod Gadus morhua

Accurate prediction of species distribution shifts in the face of climate change requires a sound understanding of population diversity and local adaptations. Previous modeling has suggested that global warming will lead to increased abundance of Atlantic cod (Gadus morhua) in the ocean around Greenland, but the dynamics of earlier abundance fluctuations are not well understood. We applied a retrospective spatiotemporal population genomics approach to examine the temporal stability of cod population structure in this region and to search for signatures of divergent selection over a 78-year period spanning major demographic changes. Analyzing >900 gene-associated single nucleotide polymorphisms in 847 individuals, we identified four genetically distinct groups that exhibited varying spatial distributions with considerable overlap and mixture. The genetic composition had remained stable over decades at some spawning grounds, whereas complete population replacement was evident at others. Observations of elevated differentiation in certain genomic regions are consistent with adaptive divergence between the groups, indicating that they may respond differently to environmental variation. Significantly increased temporal changes at a subset of loci also suggest that adaptation may be ongoing. These findings illustrate the power of spatiotemporal population genomics for revealing biocomplexity in both space and time and for informing future fisheries management and conservation efforts

Archived DNA reveals fisheries and climate induced collapse of a major fishery

Fishing and climate change impact the demography of marine fishes, but it is generally ignored that many species are made up of genetically distinct locally adapted populations that may show idiosyncratic responses to environmental and anthropogenic pressures. Here, we track 80 years of Atlantic cod (Gadus morhua) population dynamics in West Greenland using DNA from archived otoliths in combination with fish population and niche based modeling. We document how the interacting effects of climate change and high fishing pressure lead to dramatic spatiotemporal changes in the proportions and abundance of different genetic populations, and eventually drove the cod fishery to a collapse in the early 1970s. Our results highlight the relevance of fisheries management at the level of genetic populations under future scenarios of climate change

Data from: Practical low-coverage genomewide sequencing of hundreds of individually barcoded samples for population and evolutionary genomics in nonmodel species

Today most population genomic studies of nonmodel organisms either sequence a subset of the genome deeply in each individual or sequence pools of unlabelled individuals. With a step-by-step workflow, we illustrate how low-coverage whole-genome sequencing of hundreds of individually barcoded samples is now a practical alternative strategy for obtaining genomewide data on a population scale. We used a highly efficient protocol to generate high-quality libraries for ~6.5 USD from each of 876 Atlantic silversides (a teleost fish with a genome size ~730 Mb) that we sequenced to 1–4× genome coverage. In the absence of a reference genome, we developed a bioinformatic pipeline for mapping the genomic reads to a de novo assembled reference transcriptome. This provides an ‘in silico’ method for exome capture that avoids the complexities and expenses of using wet chemistry for target isolation. Using novel tools for analysis of low-coverage data, we extracted population allele frequencies, individual genotype likelihoods and polymorphism data for 2 504 335 SNPs across the exome for the 876 fish. To illustrate the use of the resulting data, we present a preliminary analysis of geographical patterns in the exome data and a comparison of complete mitochondrial genome sequences for each individual (constructed from the low-coverage data) that show population colonization patterns along the US east coast. With a total cost per sample of less than 50 USD (including sequencing) and ability to prepare 96 libraries in only 5 h, our approach adds a viable new option to the population genomics toolbox

QC_Stats_Therkildsen_and_Palumbi

This Excel file contains the raw data used to evaluate the quality of sequencing libraries. Please see the tab "Column title descriptions" for explanations of the column headers

Commands for running the pipeline

This file contains a list of the programs and parameter values we used for the bioinformatic analysis pipeline described in the paper. Note that the sequence of steps described in this file was designed to carefully evaluate the quality of sequencing data generated with a new library preparation method. For routine applications (where it may not be necessary to quantify the loss of data at each individual filtering step), several steps can be skipped, optimized or combined for a more streamlined pipeline