Nina Overgaard Therkildsen

The simple fool's guide to population genomics via RNA‐Seq: an introduction to high‐throughput sequencing data analysis

High‐throughput sequencing technologies are currently revolutionizing the field of biology and medicine, yet bioinformatic challenges in analysing very large data sets have slowed the adoption of these technologies by the community of population biologists. We introduce the ‘Simple Fools Guide to Population Genomics via RNA‐seq’ (SFG), a document intended to serve as an easy‐to‐follow protocol, walking a user through one example of high‐throughput sequencing data analysis of nonmodel organisms. It is by no means an exhaustive protocol, but rather serves as an introduction to the bioinformatic methods used in population genomics, enabling a user to gain familiarity with basic analysis steps. The SFG consists of two parts. This document summarizes the steps needed and lays out the basic themes for each and a simple approach to follow. The second document is the full SFG, publicly available at http://sfg.stanford.edu, that includes detailed protocols for data processing and analysis, along with a repository of custom‐made scripts and sample files. Steps included in the SFG range from tissue collection to de novo assembly, blast annotation, alignment, gene expression, functional enrichment, SNP detection, principal components and FST outlier analyses. Although the technical aspects of population genomics are changing very quickly, our hope is that this document will help population biologists with little to no background in high‐throughput sequencing and bioinformatics to more quickly adopt these new techniques.

A genomic island linked to ecotype divergence in Atlantic cod

The genomic architecture underlying ecological divergence and ecological speciation with gene flow is still largely unknown for most organisms. One central question is whether divergence is genome‐wide or localized in ‘genomic mosaics’ during early stages when gene flow is still pronounced. Empirical work has so far been limited, and the relative impacts of gene flow and natural selection on genomic patterns have not been fully explored. Here, we use ecotypes of Atlantic cod to investigate genomic patterns of diversity and population differentiation in a natural system characterized by high gene flow and large effective population sizes, properties which theoretically could restrict divergence in local genomic regions. We identify a genomic region of strong population differentiation, extending over approximately 20 cM, between pairs of migratory and stationary ecotypes examined at two different localities. Furthermore, the region is characterized by markedly reduced levels of genetic diversity in migratory ecotype samples. The results highlight the genomic region, or ‘genomic island’, as potentially associated with ecological divergence and suggest the involvement of a selective sweep. Finally, we also confirm earlier findings of localized genomic differentiation in three other linkage groups associated with divergence among eastern Atlantic populations. Thus, although the underlying mechanisms are still unknown, the results suggest that ‘genomic mosaics’ of differentiation may even be found under high levels of gene flow and that marine fishes may provide insightful model systems for studying and identifying initial targets of selection during ecological divergence.

Evolutionary impact assessment: accounting for evolutionary consequences of fishing in an ecosystem approach to fisheries management

Managing fisheries resources to maintain healthy ecosystems is one of the main goals of the ecosystem approach to fisheries (EAF). While a number of international treaties call for the implementation of EAF, there are still gaps in the underlying methodology. One aspect that has received substantial scientific attention recently is fisheries-induced evolution (FIE). Increasing evidence indicates that intensive fishing has the potential to exert strong directional selection on life-history traits, behaviour, physiology, and morphology of exploited fish. Of particular concern is that reversing evolutionary responses to fishing can be much more difficult than reversing demographic or phenotypically plastic responses. Furthermore, like climate change, multiple agents cause FIE, with effects accumulating over time. Consequently, FIE may alter the utility derived from fish stocks, which in turn can modify the monetary value living aquatic resources provide to society. Quantifying and predicting the evolutionary effects of fishing is therefore important for both ecological and economic reasons. An important reason this is not happening is the lack of an appropriate assessment framework. We therefore describe the evolutionary impact assessment (EvoIA) as a structured approach for assessing the evolutionary consequences of fishing and evaluating the predicted evolutionary outcomes of alternative management options. EvoIA can contribute to EAF by clarifying how evolution may alter stock properties and ecological relations, support the precautionary approach to fisheries management by addressing a previously overlooked source of uncertainty and risk, and thus contribute to sustainable fisheries.

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.

Large effective population size and temporal genetic stability in Atlantic cod (Gadus morhua) in the southern Gulf of St. Lawrence

Worldwide, many commercial fish stocks have experienced dramatic declines due to overfishing. Such fisheries-induced population reductions could potentially erode the genetic diversity of marine fish populations. Based on analyses of DNA extracted from archived and contemporary samples, this paper compares the genetic variability at nine microsatellite loci in a Canadian population of Atlantic cod (Gadus morhua) over 80 years, spanning from before the fish- ery intensified to now when the population is at historically low abundance. Extensively validated genetic data from the temporally spaced samples were used to estimate the effective population size. Over the period, we observed no loss of ei- ther heterozygosity or allelic diversity. Several of the estimation methods applied could not distinguish the effective popu- lation size from infinity, and the lower 95% confidence limit on estimates was generally >500, suggesting that the effective population size is probably considerably larger than this. Hence, it appears that the southern Gulf of St. Lawrence cod stock has maintained genetic variability to sustain future evolution despite a dramatic population decline. Resume ´ : Alechelle globale, plusieurs stocks commerciaux de poissons ont connu des declins spectaculaires a cause de la surpeche. De telles reductions demographiques dues a la peche pourraient potentiellement diminuer la diversitegene ´- tique des populations de poissons marins. Sappuyant sur des analyses dADN faites sur des echantillons archivese t contemporains, notre etude compare la variabilitegenetique a neuf locus microsatellites dans une population canadienne de morues franches (Gadus morhua) sur une periode de 80 ans couvrant le temps davant lintensification de la peche jus- quamaintenant alors que la population connaoˆt une abondance historiquement basse. Des donnees genetiques solidement validees provenant dechantillons etales dans le temps ont servi a estimer la taille effective de la population. Durant cette periode, nous nobservons de perte ni dheterozygotie, ni de diversiteallelique. Plusieurs des methodes destimation utili- sees ne peuvent pas distinguer la taille effective de la population de linfini et la limite inferieure de lintervalle de confiance de 95 % des estimations est generalement >500, ce qui laisse croire que la taille effective est probablement beaucoup plus elevee que ce nombre. Il apparaoˆt donc que le stock de morues du sud du golfe du Saint-Laurent a conserve ´ la variabilitegenetique necessaire pour assurer son evolution future, malgreson declin demographique spectaculaire. (Traduit par la Redaction)

Microevolution in time and space: SNP analysis of historical DNA reveals dynamic signatures of selection in Atlantic cod.

Little is known about how quickly natural populations adapt to changes in their environment and how temporal and spatial variation in selection pressures interact to shape patterns of genetic diversity. We here address these issues with a series of genome scans in four overfished populations of Atlantic cod (Gadus morhua) studied over an 80‐year period. Screening of >1000 gene‐associated single‐nucleotide polymorphisms (SNPs) identified 77 loci that showed highly elevated levels of differentiation, likely as an effect of directional selection, in either time, space or both. Exploratory analysis suggested that temporal allele frequency shifts at certain loci may correlate with local temperature variation and with life history changes suggested to be fisheries induced. Interestingly, however, largely nonoverlapping sets of loci were temporal outliers in the different populations and outliers from the 1928 to 1960 period showed almost complete stability during later decades. The contrasting microevolutionary trajectories among populations resulted in sequential shifts in spatial outliers, with no locus maintaining elevated spatial differentiation throughout the study period. Simulations of migration coupled with observations of temporally stable spatial structure at neutral loci suggest that population replacement or gene flow alone could not explain all the observed allele frequency variation. Thus, the genetic changes are likely to at least partly be driven by highly dynamic temporally and spatially varying selection. These findings have important implications for our understanding of local adaptation and evolutionary potential in high gene flow organisms and underscore the need to carefully consider all dimensions of biocomplexity for evolutionarily sustainable management.

Population Genomics of Marine Fishes: Next-Generation Prospects and Challenges

Over the past few years, technological advances have facilitated giant leaps forward in our ability to generate genome-wide molecular data, offering exciting opportunities for gaining new insights into the ecology and evolution of species where genomic information is still limited. Marine fishes are valuable organisms for advancing our understanding of evolution on historical and contemporary time scales, and here we highlight areas in which research on these species is likely to be particularly important in the near future. These include possibilities for gaining insights into processes on ecological time scales, identifying genomic signatures associated with population divergence under gene flow, and determining the genetic basis of phenotypic traits. We also consider future challenges pertaining to the implementation of genome-wide coverage through next-generation sequencing and genotyping methods in marine fishes. Complications associated with fast decay of linkage disequilibrium, as expected for species with large effective population sizes, and the possibility that adaptation is associated with both soft selective sweeps and polygenic selection, leaving complex genomic signatures in natural populations, are likely to challenge future studies. However, the combination of high genome coverage and new statistical developments offers promising solutions. Thus, the next generation of studies is likely to truly facilitate the transition from population genetics to population genomics in marine fishes. This transition will advance our understanding of basic evolutionary processes and will offer new possibilities for conservation and management of valuable marine resources.

Identifying genetic signatures of selection in a non-model species, alpine gentian ( Gentiana nivalis L.), using a landscape genetic approach

It is generally accepted that most plant populations are locally adapted. Yet, understanding how environmental forces give rise to adaptive genetic variation is a challenge in conservation genetics and crucial to the preservation of species under rapidly changing climatic conditions. Environmental variation, phylogeographic history, and population demographic processes all contribute to spatially structured genetic variation, however few current models attempt to separate these confounding effects. To illustrate the benefits of using a spatially-explicit model for identifying potentially adaptive loci, we compared outlier locus detection methods with a recently-developed landscape genetic approach. We analyzed 157 loci from samples of the alpine herb Gentiana nivalis collected across the European Alps. Principle coordinates of neighbor matrices (PCNM), eigenvectors that quantify multi-scale spatial variation present in a data set, were incorporated into a landscape genetic approach relating AFLP frequencies with 23 environmental variables. Four major findings emerged. 1) Fifteen loci were significantly correlated with at least one predictor variable (Radj2xa0>xa00.5). 2) Models including PCNM variables identified eight more potentially adaptive loci than models run without spatial variables. 3) When compared to outlier detection methods, the landscape genetic approach detected four of the same loci plus 11 additional loci. 4) Temperature, precipitation, and solar radiation were the three major environmental factors driving potentially adaptive genetic variation in G. nivalis. Techniques presented in this paper offer an efficient method for identifying potentially adaptive genetic variation and associated environmental forces of selection, providing an important step forward for the conservation of non-model species under global change.

Conserving marine biodiversity: insights from life-history trait candidate genes in Atlantic cod (Gadus morhua)

Recent technological developments have facilitated an increased focus on identifying genomic regions underlying adaptive trait variation in natural populations, and it has been advocated that this information should be important for designating population units for conservation. In marine fishes, phenotypic studies have suggested adaptation through divergence of life-history traits among natural populations, but the distribution of adaptive genetic variation in these species is still relatively poorly known. In this study, we extract information about the geographical distribution of genetic variation for 33 single nucleotide polymorphisms (SNPs) associated with life-history trait candidate genes, and compare this to variation in 70 putatively neutral SNPs in Atlantic cod (Gadus morhua). We analyse samples covering the major population complexes in the eastern Atlantic and find strong evidence for non-neutral levels and patterns of population structuring for several of the candidate gene-associated markers, including two SNPs in the growth hormone 1 gene. Thus, this study aligns with findings from phenotypic studies, providing molecular data strongly suggesting that these or closely linked genes are under selection in natural populations of Atlantic cod. Furthermore, we find that patterns of variation in outlier markers do not align with those observed at selectively neutral markers, and that outlier markers identify conservation units on finer geographical scales than those revealed when analysing only neutral markers. Accordingly, results also suggest that information about adaptive genetic variation will be useful for targeted conservation and management in this and other marine species.

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.