Einar Eg Nielsen

Application of SNPs for population genetics of nonmodel organisms: new opportunities and challenges

Recent improvements in the speed, cost and accuracy of next generation sequencing are revolutionizing the discovery of single nucleotide polymorphisms (SNPs). SNPs are increasingly being used as an addition to the molecular ecology toolkit in nonmodel organisms, but their efficient use remains challenging. Here, we discuss common issues when employing SNP markers, including the high numbers of markers typically employed, the effects of ascertainment bias and the inclusion of nonneutral loci in a marker panel. We provide a critique of considerations specifically associated with the application and population genetic analysis of SNPs in nonmodel taxa, focusing specifically on some of the most commonly applied methods.

Matrilinear phylogeography of Atlantic salmon ( Salmo salar L.) in Europe and postglacial colonization of the Baltic Sea area

Sixty‐four samples from 46 salmon populations totalling 2369 specimens were used for polymerase chain reaction–restriction fragment length polymorphism (PCR–RFLP) analysis of the mitochondrial ND1 region. The final analyses included 3095 specimens from 60 populations in Northern Europe. A subsample was analysed by RFLP of ND3/4/5/6. Representative RFLP haplotypes from different parts of the distribution area were sequenced and the phylogeny of European haplotypes and their relations to the North American lineage was described. The four common European haplotypes derive from the ancestral ND1‐BBBA (rooting the European clade to the North American) by one‐step substitutions: AAAA < AABA < BBBA > BBBB. The Swedish west‐coast populations differ from the geographically close southern Baltic, indicating absence of inward and limited outward gene flow through the Danish straits during the last 8000 years. Within the Baltic Sea, only three ND1 haplotypes were detected and there was no variation for ND3/4/5/6. In the whole southern Baltic and in lakes Vänern, Ladoga and Onega the haplotype AABA dominated. Proposed postglacial colonization routes to the Baltic Sea are discussed in relation to the haplotype distribution pattern.

Adaptive differences in gene expression in European flounder (Platichthys flesus)

Population structure was previously believed to be very limited or absent in classical marine fishes, but recently, evidence of weakly differentiated local populations has been accumulating using noncoding microsatellite markers. However, the evolutionary significance of such minute genetic differences remains unknown. Therefore, in order to elucidate the relationship between genetic markers and adaptive divergence among populations of marine fishes, we combined cDNA microarray and microsatellite analysis in European flounders (Platichthys flesus). We demonstrate that despite extremely low levels of neutral genetic divergence, a high number of genes were significantly differentially expressed between North Sea and Baltic Sea flounders maintained in a long‐term reciprocal transplantation experiment mimicking natural salinities. Several of the differentially regulated genes could be directly linked to fitness traits. These findings demonstrate that flounders, despite little neutral genetic divergence between populations, are differently adapted to local environmental conditions and imply that adaptation in gene expression could be common in other marine organisms with similar low levels of population subdivision.

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.

Genetic diversity and connectivity within Mytilus spp. in the subarctic and Arctic

Climate changes in the Arctic are predicted to alter distributions of marine species. However, such changes are difficult to quantify because information on present species distribution and the genetic variation within species is lacking or poorly examined. Blue mussels, Mytilus spp., are ecosystem engineers in the coastal zone globally. To improve knowledge of distribution and genetic structure of the Mytilus edulis complex in the Arctic, we analyzed 81 SNPs in 534 Mytilus spp. individuals sampled at 13 sites to provide baseline data for distribution and genetic variation of Mytilus mussels in the European Arctic. Mytilus edulis was the most abundant species found with a clear genetic split between populations in Greenland and the Eastern Atlantic. Surprisingly, analyses revealed the presence of Mytilus trossulus in high Arctic NW Greenland (77°N) and Mytilus galloprovincialis or their hybrids in SW Greenland, Svalbard, and the Pechora Sea. Furthermore, a high degree of hybridization and introgression between species was observed. Our study highlights the importance of distinguishing between congener species, which can display local adaptation and suggests that information on dispersal routes and barriers is essential for accurate predictions of regional susceptibility to range expansions or invasions of boreal species in the Arctic.

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.

Extracting DNA from ‘jaws’: high yield and quality from archived tiger shark (Galeocerdo cuvier) skeletal material

Archived specimens are highly valuable sources of DNA for retrospective genetic/genomic analysis. However, often limited effort has been made to evaluate and optimize extraction methods, which may be crucial for downstream applications. Here, we assessed and optimized the usefulness of abundant archived skeletal material from sharks as a source of DNA for temporal genomic studies. Six different methods for DNA extraction, encompassing two different commercial kits and three different protocols, were applied to material, so‐called bio‐swarf, from contemporary and archived jaws and vertebrae of tiger sharks (Galeocerdo cuvier). Protocols were compared for DNA yield and quality using a qPCR approach. For jaw swarf, all methods provided relatively high DNA yield and quality, while large differences in yield between protocols were observed for vertebrae. Similar results were obtained from samples of white shark (Carcharodon carcharias). Application of the optimized methods to 38 museum and private angler trophy specimens dating back to 1912 yielded sufficient DNA for downstream genomic analysis for 68% of the samples. No clear relationships between age of samples, DNA quality and quantity were observed, likely reflecting different preparation and storage methods for the trophies. Trial sequencing of DNA capture genomic libraries using 20 000 baits revealed that a significant proportion of captured sequences were derived from tiger sharks. This study demonstrates that archived shark jaws and vertebrae are potential high‐yield sources of DNA for genomic‐scale analysis. It also highlights that even for similar tissue types, a careful evaluation of extraction protocols can vastly improve DNA yield.

Geographic extent of introgression in Sebastes mentella and its effect on genetic population structure

Genetic population structure is often used to identify management units in exploited species, but the extent of genetic differentiation may be inflated by geographic variation in the level of hybridization between species. We identify the genetic population structure of Sebastes mentella and investigate possible introgression within the genus by analyzing 13 microsatellites in 2,562 redfish specimens sampled throughout the North Atlantic. The data support an historical divergence between the “shallow” and “deep” groups, beyond the Irminger Sea where they were described previously. A third group, “slope,” has an extended distribution on the East Greenland Shelf, in addition to earlier findings on the Icelandic slope. Furthermore, S. mentella from the Northeast Arctic and Northwest Atlantic waters are genetically different populations. In both areas, interspecific introgression may influence allele frequency differences among populations. Evidence of introgression was found for almost all the identified Sebastes gene pools, but to a much lower extent than suggested earlier. Greenland waters appear to be a sympatric zone for many of the genetically independent Sebastes groups. This study illustrates that the identified groups maintain their genetic integrity in this region despite introgression.

A microsatellite baseline for genetic stock identification of European Atlantic salmon (Salmo salar L.)

This work forms part of the SALSEA-Merge research project (Project No. 212529) and was funded by the European Union under theme six of the Seventh Framework programme. It was also co-sponsored by the Atlantic Salmon Trust and the Total Foundation, who we thank for financial support. PMcG and JC were partly supported by the Beaufort Marine Research Award in Fish Population Genetics funded by the Irish Government under the Sea Change Programme. The work was also supported under financial support of the program of fundamental research of Presidium of RAS “Searching fundamental scientific investigations in the interests of development of the Arctic zone of Russian Federation.”

Population structure and connectivity of tiger sharks (Galeocerdo cuvier) across the indo-pacific ocean basin

Population genetic structure using nine polymorphic nuclear microsatellite loci was assessed for the tiger shark (Galeocerdo cuvier) at seven locations across the Indo-Pacific, and one location in the southern Atlantic. Genetic analyses revealed considerable genetic structuring (FST > 0.14, p < 0.001) between all Indo-Pacific locations and Brazil. By contrast, no significant genetic differences were observed between locations from within the Pacific or Indian Oceans, identifying an apparent large, single Indo-Pacific population. A lack of differentiation between tiger sharks sampled in Hawaii and other Indo-Pacific locations identified herein is in contrast to an earlier global tiger shark nDNA study. The results of our power analysis provide evidence to suggest that the larger sample sizes used here negated any weak population subdivision observed previously. These results further highlight the need for cross-jurisdictional efforts to manage the sustainable exploitation of large migratory sharks like G. cuvier.