Evolving southwest African response to abrupt deglacial North Atlantic climate change events

Climate change during the last deglaciation was strongly influenced by the „bipolar seesaw‟, producing antiphase climate responses between the North and South Atlantic. However, mounting evidence demands refinements of this model, with the occurrence of abrupt events in southern low to mid latitudes occurring in-phase with North Atlantic climate. Improved constraints on the north-south phasing and spatial extent of these events are therefore critical to
understanding the mechanisms that propagate abrupt events within the climate system. We present a 19,400 year multi-proxy record of climate change obtained from a rock hyrax midden in southernmost Africa. Arid anomalies in phase with the Younger Dryas and 8.2 ka events are apparent, indicating a clear shift in the influence of the bipolar seesaw, which diminished as the Earth warmed, and was succeeded after ~14.6 ka by the emergence of a dominant interhemispheric atmospheric teleconnection.

North Atlantic surface ocean 14C variability between 30 – 40 kyr. BP

The Faroe-Shetland channel is situated in the main path of the inflow of warm North Atlantic surface water to the Nordic seas and further provides an escape route for the cold Norwegian Sea Deep Water. AMS 14C dates of planktonic foraminifera covering Marine Isotope Stage 3 from two cores in the Faroe-Shetland channel will be used to trace past variability of the Atlantic Meridional Overturning Circulation (AMOC). The reservoir age R shows considerable variability ranging between 50 to 2750 14C years. In particular high R values are observed during Heinrich event 4 (H4) with values around 1550 14C years and during the Laschamp magnetic excursion with R values as high as 2700 14C years. The period between Greenland interstadial 8 (GI8) and GI5 show highly variable R values with interstadial R values around 500 – 650 14C years, i.e. slightly higher than ‘normal’, whereas stadials show either significantly higher or lower R values. From GI5 towards the Last Glacial Maximum R values are generally around 1000 14C years or higher. Using magnetic susceptibility, IRD and δ13C and δ18O values measured on the planktic foraminifera species Neogloboquadrina pachyderma, we compare the observed R variability with reconstructed changes in the Atlantic Meridional Overturning Circulation (AMOC). Furthermore a climate model of intermediate complexity (GENIE) including 14C is used as conceptual tool for identifying oceanographic configuration explaining the observed R variability.

High-resolution genetic analysis reveals extensive gene flow within the jellyfish Pelagia noctiluca (Scyphozoa) in the North Atlantic and Mediterranean Sea

Despite the importance of gelatinous zooplankton as components of marine ecosystems, both ecologically and socio-economically, relatively little is known about population persistence or connectivity in jellyfish. In the present study, we employed a combination of nuclear microsatellite markers and sequence data from the mitochondrial cytochrome oxidase I (COI) gene to determine levels and patterns of population genetic structuring in the holoplanktonic jellyfish Pelagia noctiluca across the northeast Atlantic Ocean and Mediterranean Sea. Our results indicate a high degree of connectivity in P. noctiluca, with little evidence of geographical structuring of genetic variation. A small but significant differentiation of Atlantic Ocean and Mediterranean stocks was detected based on the microsatellite data, but no evidence of differentiation was observed with the mtDNA, probably due to the higher power of the microsatellites to detect low levels of genetic structuring. Two clearly distinct groups of genotypes were observed within the mtDNA COI, which probably diverged in the early Pleistocene, but with no evidence of geographical structuring. Palaeodistribution modelling of P. noctiluca at the Last Glacial Maximum (LGM; ca. 21 KYA) indicated large areas of suitable habitat south of the species’ current-day distribution, with little reduction in area. The congruent evidence for minimal genetic differentiation from the nuclear microsatellites and the mtDNA, coupled with the results of the palaeodistribution modelling, supports the idea of long-term population stability and connectivity, thus providing key insights into the population dynamics and demography of this important species

Rhodopsin gene polymorphism associated with divergent light environments in atlantic cod

The spectral sensitivity of visual pigments in vertebrate eyes is optimized for specific light conditions. One of such pigments, rhodopsin (RH1), mediates dim-light vision. Amino acid replacements at tuning sites may alter spectral sensitivity, providing a mechanism to adapt to ambient light conditions and depth of habitat in fish. Here we present a first investigation of RH1 gene polymorphism among two ecotypes of Atlantic cod in Icelandic waters, which experience divergent light environments throughout the year due to alternative foraging behaviour. We identified one synonymous single nucleotide polymorphism (SNP) in the RH1 protein coding region and one in the 3' untranslated region (3'-UTR) that are strongly divergent between these two ecotypes. Moreover, these polymorphisms coincided with the well-known panthophysin (Pan I) polymorphism that differentiates coastal and frontal (migratory) populations of Atlantic cod. While the RH1 SNPs do not provide direct inference for a specific molecular mechanism, their association with this dim-sensitive pigment indicates the involvement of the visual system in local adaptation of Atlantic cod.

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.

SNP discovery using Next Generation Transcriptomic Sequencing in Atlantic herring (Clupea harengus)

The introduction of Next Generation Sequencing (NGS) has revolutionised population genetics, providing studies of non-model species with unprecedented genomic coverage, allowing evolutionary biologists to address questions previously far beyond the reach of available resources. Furthermore, the simple mutation model of Single Nucleotide Polymorphisms (SNPs) permits cost-effective high-throughput genotyping in thousands of individuals simultaneously. Genomic resources are scarce for the Atlantic herring (Clupea harengus), a small pelagic species that sustains high revenue fisheries. This paper details the development of 578 SNPs using a combined NGS and high-throughput genotyping approach. Eight individuals covering the species distribution in the eastern Atlantic were bar-coded and multiplexed into a single cDNA library and sequenced using the 454 GS FLX platform. SNP discovery was performed by de novo sequence clustering and contig assembly, followed by the mapping of reads against consensus contig sequences. Selection of candidate SNPs for genotyping was conducted using an in silico approach. SNP validation and genotyping were performed simultaneously using an Illumina 1,536 GoldenGate assay. Although the conversion rate of candidate SNPs in the genotyping assay cannot be predicted in advance, this approach has the potential to maximise cost and time efficiencies by avoiding expensive and time-consuming laboratory stages of SNP validation. Additionally, the in silico approach leads to lower ascertainment bias in the resulting SNP panel as marker selection is based only on the ability to design primers and the predicted presence of intron-exon boundaries. Consequently SNPs with a wider spectrum of minor allele frequencies (MAFs) will be genotyped in the final panel. The genomic resources presented here represent a valuable multi-purpose resource for developing informative marker panels for population discrimination, microarray development and for population genomic studies in the wild.

Environmental selection on transcriptome-derived SNPs in a high gene flow marine fish, the Atlantic herring (Clupea harengus)

High gene flow is considered the norm for most marine organisms and is expected to limit their ability to adapt to local environments. Few studies have directly compared the patterns of differentiation at neutral and selected gene loci in marine organisms. We analysed a transcriptome-derived panel of 281 SNPs in Atlantic herring (Clupea harengus), a highly migratory small pelagic fish, for elucidating neutral and selected genetic variation among populations and to identify candidate genes for environmental adaptation. We analysed 607 individuals from 18 spawning locations in the northeast Atlantic, including two temperature clines (5-12 °C) and two salinity clines (5-35‰). By combining genome scan and landscape genetic analyses, four genetically distinct groups of herring were identified: Baltic Sea, Baltic-North Sea transition area, North Sea/British Isles and North Atlantic; notably, samples exhibited divergent clustering patterns for neutral and selected loci. We found statistically strong evidence for divergent selection at 16 outlier loci on a global scale, and significant correlations with temperature and salinity at nine loci. On regional scales, we identified two outlier loci with parallel patterns across temperature clines and five loci associated with temperature in the North Sea/North Atlantic. Likewise, we found seven replicated outliers, of which five were significantly associated with low salinity across both salinity clines. Our results reveal a complex pattern of varying spatial genetic variation among outlier loci, likely reflecting adaptations to local environments. In addition to disclosing the fine scale of local adaptation in a highly vagile species, our data emphasize the need to preserve functionally important biodiversity.

Gene-associated markers can assign origin in a weakly structured fish, Atlantic herring

Regulations on the exploitation of populations of commercially important fish species and the ensuing consumer interest in sustainable products have increased the need to accurately identify the population of origin of fish and fish products. Although genomics-based tools have proven highly useful, there are relatively few examples in marine fish displaying accurate origin assignment. We synthesize data for 156 single-nucleotide polymorphisms typed in 1039 herring, Clupea harengus L., spanning the Northeast Atlantic to develop a tool that allows assignment of individual herring to their regional origin. We show the method's suitability to address specific biological questions, as well as management applications. We analyse temporally replicated collections from two areas, the Skagerrak (n = 81, 84, 66) and the western Baltic (n = 52, 52). Both areas harbour heavily fished mixed-origin stocks, complicating management issues. We report novel genetic evidence that herring from the Baltic Sea contribute to catches in the North Sea, and find support that western Baltic feeding aggregations mainly constitute herring from the western Baltic with contributions from the Eastern Baltic. Our study describes a general approach and outlines a database allowing individual assignment and traceability of herring across a large part of its East Atlantic distribution.

Stock structure of Atlantic herring Clupea harengus in the Norwegian Sea and adjacent waters

The genetic structure of Atlantic herring Clupea harengus L. was investigated in its north-easterly distribution in the Norwegian Sea and adjacent waters, using 23 neutral and one non-neutral (Cpa111) microsatellite loci. Fish from the suspected 2 main populations - the Norwegian spring-spawning herring (NSSH) and the Icelandic summer-spawning herring (ISSH) - were collected at spawning locations in their respective spawning seasons from 2009 to 2012. Samples were also collected from Norwegian autumn spawning locations, from different local Norwegian fjords such as the inner part of Trondheimsfjorden, Lindås pollene, Landvikvannet and Lusterfjorden, as well as from suspected Faroese spawning components. The observed level of genetic differentiation was significant but low (F_ST = 0.007) and mostly attributable to the differentiation of the local Norwegian fjord populations. The locus Cpa111, which was detected to putatively be under positive selection, exhibited the highest FST value (0.044). The observed genetic patterns were robust to exclusion of this locus. Landvikvannet herring was also genetically distinguishable from the 3 other fjord populations. In addition, the present study does not support genetic structuring among the ISSH and the NSSH.

Genetic structure of the lumpfish Cyclopterus lumpus across the North Atlantic

Lumpfish, or lumpsucker, Cyclopterus lumpus (Linnaeus, 1758) is widely distributed in the North Atlantic Ocean. It has a considerable economic value and substantial fisheries occur in several North Atlantic regions owing to the use of its fully ripe internal egg masses in the ovaries as an alternative to sturgeon caviar. Despite being intensively fished in several locations, biological knowledge is limited and no genetic structure information is available. In this study, the stock structure of C. lumpus was investigated across the North Atlantic using ten microsatellite loci. Out of ten loci, two exhibited higher level of differentiation but their inclusion/ exclusion from the analyses did not drastically change the observed genetic pattern. A total of three distinct genetic groups were detected: Maine-Canada-Greenland, Iceland-Norway and Baltic Sea. These results, discussed in terms of origin of differentiation, gene flow, and selection, showed that gene flow was rather limited among the detected groups, and also between Greenland and Maine-Canada.