Monitoring of species' genetic diversity in Europe varies greatly and overlooks potential climate change impacts.

Genetic monitoring of populations currently attracts interest in the context of the Convention on Biological Diversity but needs long-term planning and investments. However, genetic diversity has been largely neglected in biodiversity monitoring, and when addressed, it is treated separately, detached from other conservation issues, such as habitat alteration due to climate change. We report an accounting of efforts to monitor population genetic diversity in Europe (genetic monitoring effort, GME), the evaluation of which can help guide future capacity building and collaboration towards areas most in need of expanded monitoring. Overlaying GME with areas where the ranges of selected species of conservation interest approach current and future climate niche limits helps identify whether GME coincides with anticipated climate change effects on biodiversity. Our analysis suggests that country area, financial resources and conservation policy influence GME, high values of which only partially match species' joint patterns of limits to suitable climatic conditions. Populations at trailing climatic niche margins probably hold genetic diversity that is important for adaptation to changing climate. Our results illuminate the need in Europe for expanded investment in genetic monitoring across climate gradients occupied by focal species, a need arguably greatest in southeastern European countries. This need could be met in part by expanding the European Union's Birds and Habitats Directives to fully address the conservation and monitoring of genetic diversity.

Genetic assignment predicts depth of benthic settlement for 0-group Atlantic cod.

Atlantic cod is a keystone species that remains among the most economically important demersal fish in the North Atlantic. Throughout its distribution range, Atlantic cod is composed of populations with varying environmental preferences and migratory propensities. This life-history variation is likely to have contributed to the niche width and large population sizes of Atlantic cod, and its relative resilience to environmental change and exploitation. The Icelandic cod stock is currently managed as a single unit, but early research indicates population variation by depth and temperature and distinct offshore and inshore spawning components. Pelagic 0-group juveniles from different spawning grounds coexist in nursery areas around Iceland, but their genetic composition or habitat partitioning had not been examined post benthic settlement. In the current study we examine the genetic composition of Atlantic cod juvenile aggregations at nearshore nursery grounds in NW-Iceland and report distinct segregation by the depth of offshore and inshore juvenile cod. The physiological mechanism of this segregation is not known, but the pattern demonstrates the need to consider population structure at nursery grounds in the application of marine spatial planning and other area-based conservation tools.

Holocene deglaciation drove rapid genetic diversification of Atlantic walrus.

Rapid global warming is severely impacting Arctic ecosystems and is predicted to transform the abundance, distribution and genetic diversity of Arctic species, though these linkages are poorly understood. We address this gap in knowledge using palaeogenomics to examine how earlier periods of global warming influenced the genetic diversity of Atlantic walrus (Odobenus rosmarus rosmarus), a species closely associated with sea ice and shallow-water habitats. We analysed 82 ancient and historical Atlantic walrus mitochondrial genomes (mitogenomes), including now-extinct populations in Iceland and the Canadian Maritimes, to reconstruct the Atlantic walrus' response to Arctic deglaciation. Our results demonstrate that the phylogeography and genetic diversity of Atlantic walrus populations was initially shaped by the last glacial maximum (LGM), surviving in distinct glacial refugia, and subsequently expanding rapidly in multiple migration waves during the late Pleistocene and early Holocene. The timing of diversification and establishment of distinct populations corresponds closely with the chronology of the glacial retreat, pointing to a strong link between walrus phylogeography and sea ice. Our results indicate that accelerated ice loss in the modern Arctic may trigger further dispersal events, likely increasing the connectivity of northern stocks while isolating more southerly stocks putatively caught in small pockets of suitable habitat.

Genomic diversity and differentiation between island and mainland populations of white-tailed eagles (Haliaeetus albicilla).

Divergence in the face of high dispersal capabilities is a documented but poorly understood phenomenon. The white-tailed eagle (Haliaeetus albicilla) has a large geographic dispersal capability and should theoretically be able to maintain genetic homogeneity across its dispersal range. However, following analysis of the genomic variation of white-tailed eagles, from both historical and contemporary samples, clear signatures of ancient biogeographic substructure across Europe and the North-East Atlantic is observed. The greatest genomic differentiation was observed between island (Greenland and Iceland) and mainland (Denmark, Norway and Estonia) populations. The two island populations share a common ancestry from a single mainland population, distinct from the other sampled mainland populations, and despite the potential for high connectivity between Iceland and Greenland they are well separated from each other and are characterized by inbreeding and little variation. Temporal differences also highlight a pattern of regional populations persisting despite the potential for admixture. All sampled populations generally showed a decline in effective population size over time, which may have been shaped by four historical events: (1) Isolation of refugia during the last glacial period 110-115,000 years ago, (2) population divergence following the colonization of the deglaciated areas ~10,000 years ago, (3) human population expansion, which led to the settlement in Iceland ~1100 years ago, and (4) human persecution and exposure to toxic pollutants during the last two centuries.

Population structure and genetic variation of fragmented mountain birch forests in Iceland.

Betula pubescens Ehrh. (mountain birch) is the only forest-forming tree in Iceland. Since human settlement (874 AD), the continuous 25,000 to 30,000 km2 forest has shrunk to 1.200 km2 of fragmented patches, making it a good object to study population genetic consequences of habitat fragmentation and disturbance. Further, genetic studies have also shown that hybridization between the tetraploid (2n = 56) B. pubescens and the diploid (2n = 28) Betula nana L. (dwarf birch) occurs among Iceland's natural populations. This study assessed the genetic variation within and among 11 birch forests remaining across Iceland. Genotype-by-sequencing methodology provided a total of 24,585 single nucleotide polymorphisms (SNP´s), with a minor allele frequency >5% for genetic analyses. The analysis showed similar diversity within forests, suggesting that fragmentation and hybridization have had a limited effect on the genetic variation within sites. A clear genetic divergence is found among forests from the different regions of Iceland that may reflect historical isolation; the differentiation between forests increased with geographic distances reflecting isolation by distance. Information on the distribution of genetic variation of birch in Iceland is essential for its conservation and to establish genotype-phenotype associations to predict responses to new environmental conditions imposed by climate change and novel biotic/abiotic stressors.

Microbial Response to Increased Temperatures Within a Lava-Induced Hydrothermal System in Iceland: An Analogue for the Habitability of Volcanic Terrains on Mars.

Fossil hydrothermal systems on Mars are important exploration targets because they may have once been habitable and could still preserve evidence of microbial life. We investigated microbial communities within an active lava-induced hydrothermal system associated with the 2014-2015 eruption of Holuhraun in Iceland as a Mars analogue. In 2016, the microbial composition in the lava-heated water differed substantially from that of the glacial river and spring water sources that fed into the system. Several taxonomic and metabolic groups were confined to the water emerging from the lava and some showed the highest sequence similarities to subsurface ecosystems, including to the predicted thermophilic and deeply branching Candidatus Acetothermum autotrophicum. Measurements show that the communities were affected by temperature and other environmental factors. In particular, comparing glacial river water incubated in situ (5.7°C, control) with glacial water incubated within a lava-heated stream (17.5°C, warm) showed that microbial abundance, richness, and diversity increased in the warm treatment compared with the control, with the predicted major metabolism shifting from lithotrophy toward organotrophy and possibly phototrophy. In addition, thermophilic bacteria isolated from the lava-heated water and a nearby acidic hydrothermal system included the known endospore-formers Geobacillus stearothermophilus and Paenibacillus cisolokensis as well as a potentially novel taxon within the order Hyphomicrobiales. Similar lava-water interactions on Mars could therefore have generated habitable environments for microbial communities.

Invertebrate communities in springs across a gradient in thermal regimes.

In many respects, freshwater springs can be considered as unique ecosystems on the fringe of aquatic habitats. This integrates their uniqueness in terms of stability of environmental metrics. The main objective of our study was to evaluate how environmental variables may shape invertebrate diversity and community composition in different freshwater spring types and habitats within. In order to do so, we sampled invertebrates from 49 springs in Iceland, where we included both limnocrene and rheocrene springs. At each site, samples were taken from the benthic substrate of the spring ("surface") and the upwelling groundwater at the spring source ("source"). To collect invertebrates from the spring sources we used a modified method of "electrobugging" and Surber sampler for collecting invertebrates from the surface. In total, 54 invertebrate taxa were identified, mostly Chironomidae (Diptera). Chironomid larvae also dominated in terms of abundance (67%), followed by Ostracoda (12%) and Copepoda (9%). The species composition in the surface samples differed considerably between rheocrene and limnocrene springs and was characterised by several indicator species. Alpha diversity was greater at the surface of springs than at the source, but the beta diversity was higher at the source. Diversity, as summarized by taxa richness and Shannon diversity, was negatively correlated with temperature at the surface. At the source, on the other hand, Shannon diversity increased with temperature. The community assembly in springs appears to be greatly affected by water temperature, with the source community of hot springs being more niche-assembled (i.e., affected by mechanisms of tolerance and adaptation) than the source community of cold springs, which is more dispersal-assembled (i.e., by mechanisms of drift and colonization).

Evaluation of four methods to identify the homozygotic sex chromosome in small populations.

BACKGROUND: Whole genomes are commonly assembled into a collection of scaffolds and often lack annotations of autosomes, sex chromosomes, and organelle genomes (i.e., mitochondrial and chloroplast). As these chromosome types differ in effective population size and can have highly disparate evolutionary histories, it is imperative to take this information into account when analysing genomic variation. Here we assessed the accuracy of four methods for identifying the homogametic sex chromosome in a small population using two whole genome sequences (WGS) and 133 RAD sequences of white-tailed eagles (Haliaeetus albicilla): i) difference in read depth per scaffold in a male and a female, ii) heterozygosity per scaffold in a male and a female, iii) mapping to the reference genome of a related species (chicken) with annotated sex chromosomes, and iv) analysis of SNP-loadings from a principal components analysis (PCA), based on the low-depth RADseq data. RESULTS: The best performing approach was the reference mapping (method iii), which identified 98.12% of the expected homogametic sex chromosome (Z). Read depth per scaffold (method i) identified 86.41% of the homogametic sex chromosome with few false positives. SNP-loading scores (method iv) identified 78.6% of the Z-chromosome and had a false positive discovery rate of more than 10%. Heterozygosity per scaffold (method ii) did not provide clear results due to a lack of diversity in both the Z and autosomal chromosomes, and potential interference from the heterogametic sex chromosome (W). The evaluation of these methods also revealed 10 Mb of putative PAR and gametologous regions. CONCLUSION: Identification of the homogametic sex chromosome in a small population is best accomplished by reference mapping or examining differences in read depth between sexes.

Fishing for DNA? Designing baits for population genetics in target enrichment experiments: Guidelines, considerations and the new tool supeRbaits.

Targeted sequencing is an increasingly popular next-generation sequencing (NGS) approach for studying populations that involves focusing sequencing efforts on specific parts of the genome of a species of interest. Methodologies and tools for designing targeted baits are scarce but in high demand. Here, we present specific guidelines and considerations for designing capture sequencing experiments for population genetics for both neutral genomic regions and regions subject to selection. We describe the bait design process for three diverse fish species: Atlantic salmon, Atlantic cod and tiger shark, which was carried out in our research group, and provide an evaluation of the performance of our approach across both historical and modern samples. The workflow used for designing these three bait sets has been implemented in the R-package supeRbaits, which encompasses our considerations and guidelines for bait design for the benefit of researchers and practitioners. The supeRbaits R-package is user-friendly and versatile. It is written in C++ and implemented in R. supeRbaits and its manual are available from Github: https://github.com/BelenJM/supeRbaits.

A supergene underlies linked variation in color and morphology in a Holarctic songbird.

The genetic architecture of a phenotype can have considerable effects on the evolution of a trait or species. Characterizing genetic architecture provides insight into the complexity of a given phenotype and, potentially, the role of the phenotype in evolutionary processes like speciation. We use genome sequences to investigate the genetic basis of phenotypic variation in redpoll finches (Acanthis spp.). We demonstrate that variation in redpoll phenotype is broadly controlled by a ~55-Mb chromosomal inversion. Within this inversion, we find multiple candidate genes related to melanogenesis, carotenoid coloration, and bill shape, suggesting the inversion acts as a supergene controlling multiple linked traits. A latitudinal gradient in ecotype distribution suggests supergene driven variation in color and bill morphology are likely under environmental selection, maintaining supergene haplotypes as a balanced polymorphism. Our results provide a mechanism for the maintenance of ecotype variation in redpolls despite a genome largely homogenized by gene flow.

RAD sequencing of common whelk, Buccinum undatum, reveals fine-scale population structuring in Europe and cryptic speciation within the North Atlantic.

Buccinum undatum is a subtidal gastropod that exhibits clear spatial variation in several phenotypic shell traits (color, shape, and thickness) across its North Atlantic distribution. Studies of spatial phenotypic variation exist for the species; however, population genetic studies have thus far relied on a limited set of mitochondrial and microsatellite markers. Here, we greatly expand on previous work by characterizing population genetic structure in B. undatum across the North Atlantic from SNP variation obtained by RAD sequencing. There was a high degree of genetic differentiation between Canadian and European populations (Iceland, Faroe Islands, and England) consistent with the divergence of populations in allopatry (F ST > 0.57 for all pairwise comparisons). In addition, B. undatum populations within Iceland, the Faroe Islands, and England are typified by weak but significant genetic structuring following an isolation-by-distance model. Finally, we established a significant correlation between genetic structuring in Iceland and two phenotypic traits: shell shape and color frequency. The works detailed here enhance our understanding of genetic structuring in B. undatum and establish the species as an intriguing model for future genome-wide association studies.

Bacterial diversity in Icelandic cold spring sources and in relation to the groundwater amphipod Crangonyx islandicus.

Crangonyx islandicus is a groundwater amphipod endemic to Iceland, considered to have survived the Ice Ages in subglacial refugia. Currently the species is found in spring sources in lava fields along the tectonic plate boundary of the country. The discovery of a groundwater species in this inaccessible habitat indicates a hidden ecosystem possibly based on chemoautotrophic microorganisms as primary producers. To explore this spring ecosystem, we assessed its microbial diversity and analysed whether and how the diversity varied between the amphipods and the spring water, and if was dependent on environmental factors and geological settings. Isolated DNA from spring water and from amphipods was analysed using metabarcoding methods, targeting the 16S rRNA gene. Two genera of bacteria, Halomonas and Shewanella were dominating in the amphipod samples in terms of relative abundance, but not in the groundwater samples where Flavobacterium, Pseudomonas and Alkanindiges among others were dominating. The richness of the bacteria taxa in the microbial community of the groundwater spring sources was shaped by pH level and the beta diversity was shaped by geographic locations.

Disappearance of Icelandic Walruses Coincided with Norse Settlement.

There is a growing body of evidence demonstrating the impacts of human arrival in new "pristine" environments, including terrestrial habitat alterations and species extinctions. However, the effects of marine resource utilization prior to industrialized whaling, sealing, and fishing have largely remained understudied. The expansion of the Norse across the North Atlantic offers a rare opportunity to study the effects of human arrival and early exploitation of marine resources. Today, there is no local population of walruses on Iceland, however, skeletal remains, place names, and written sources suggest that walruses existed, and were hunted by the Norse during the Settlement and Commonwealth periods (870-1262 AD). This study investigates the timing, geographic distribution, and genetic identity of walruses in Iceland by combining historical information, place names, radiocarbon dating, and genomic analyses. The results support a genetically distinct, local population of walruses that went extinct shortly after Norse settlement. The high value of walrus products such as ivory on international markets likely led to intense hunting pressure, which-potentially exacerbated by a warming climate and volcanism-resulted in the extinction of walrus on Iceland. We show that commercial hunting, economic incentives, and trade networks as early as the Viking Age were of sufficient scale and intensity to result in significant, irreversible ecological impacts on the marine environment. This is to one of the earliest examples of local extinction of a marine species following human arrival, during the very beginning of commercial marine exploitation.

Shell morphology and color of the subtidal whelk Buccinum undatum exhibit fine-scaled spatial patterns.

Geographical patterns in morphology can be the result of divergence among populations due to neutral or selective changes and/or phenotypic plasticity in response to different environments. Marine gastropods are ideal subjects on which to explore these patterns, by virtue of the remarkable intraspecific variation in life-history traits and morphology often observed across relatively small spatial scales. The ubiquitous N-Atlantic common whelk (Buccinum undatum) is well known for spatial variation in life-history traits and morphology. Previous studies on genetic population structure have revealed that it exhibits significant differentiation across geographic distances. Within Breiðafjörður Bay, a large and shallow bay in W-Iceland, genetic differentiation was demonstrated between whelks from sites separated by just 20 km. Here, we extended our previous studies on the common whelk in Breiðafjörður Bay by quantifying phenotypic variation in shell morphology and color throughout the Bay. We sought to test whether trait differentiation is dependent on geographic distance and/or environmental variability. Whelk in Breiðafjörður Bay displayed fine-scale patterns of spatial variation in shape, thickness, and color diversity. Differentiation increased with increasing distance between populations, indicating that population connectivity is limited. Both shape and color varied along a gradient from the inner part of the bay in the east to the outer part in the west. Whelk shells in the innermost part of Breiðafjörður Bay were thick with an elongate shell, round aperture, and low color diversity, whereas in the outer part of the bay the shells were thinner, rounder, with a more elongate aperture and richer color diversity. Significant site-specific difference in shell traits of the common whelk in correlation with environmental variables indicates the presence of local ecotypes and limited demographic connectivity.

Historical demography and genetic differentiation of the giant freshwater prawn Macrobrachium rosenbergii in Bangladesh based on mitochondrial and ddRAD sequence variation.

Macrobrachium rosenbergii, the giant freshwater prawn, is an important source of high quality protein and occurs naturally in rivers as well as commercial farms in South and South-East Asia, including Bangladesh. This study investigated the genetic variation and population structure of M. rosenbergii sampled from four rivers in Bangladesh (sample size ranged from 19 to 20), assessing sequence variation, both in the mitochondrial cytochrome oxidase subunit 1 (CO1) gene and in 106 single nucleotide polymorphisms (SNPs) sampled randomly from the genome with double digest RAD sequencing (ddRADseq). The mitochondrial variation presented a shallow genealogy with high haplotype diversity (h = 0.95), reflecting an expansion in population size for the last ~82 kyr. Based on the CO1 variation the current effective population size (Ne) was 9.7 × 106 (CI: 1.33 × 106 - 35.84 × 106) individuals. A significant population differentiation was observed with the mitochondrial CO1 sequence variation and based on the ddRADseq variation, which could be traced to the divergence of the population in the Naf River in the South-East border with Myanmar from the other populations. A differentiation in mtDNA haplotype frequencies was also observed between the Biskhali River and the Karnaphuli Rivers in eastern Bangladesh. This study demonstrated the use of high-throughput genotyping based on the ddRADseq method to reveal population structure at a small geographical scale for an important freshwater prawn. The information from this study can be utilized for management and conservation of this species in Bangladesh.

Population structure of Purple Sandpipers (Calidris maritima) as revealed by mitochondrial DNA and microsatellites.

The Purple Sandpiper (Calidris maritima) is a medium-sized shorebird that breeds in the Arctic and winters along northern Atlantic coastlines. Migration routes and affiliations between breeding grounds and wintering grounds are incompletely understood. Some populations appear to be declining, and future management policies for this species will benefit from understanding their migration patterns. This study used two mitochondrial DNA markers and 10 microsatellite loci to analyze current population structure and historical demographic trends. Samples were obtained from breeding locations in Nunavut (Canada), Iceland, and Svalbard (Norway) and from wintering locations along the coast of Maine (USA), Nova Scotia, New Brunswick, and Newfoundland (Canada), and Scotland (UK). Mitochondrial haplotypes displayed low genetic diversity, and a shallow phylogeny indicating recent divergence. With the exception of the two Canadian breeding populations from Nunavut, there was significant genetic differentiation among samples from all breeding locations; however, none of the breeding populations was a monophyletic group. We also found differentiation between both Iceland and Svalbard breeding populations and North American wintering populations. This pattern of divergence is consistent with a previously proposed migratory pathway between Canadian breeding locations and wintering grounds in the United Kingdom, but argues against migration between breeding grounds in Iceland and Svalbard and wintering grounds in North America. Breeding birds from Svalbard also showed a genetic signature intermediate between Canadian breeders and Icelandic breeders. Our results extend current knowledge of Purple Sandpiper population genetic structure and present new information regarding migration routes to wintering grounds in North America.

Population limitation in a non-cyclic arctic fox population in a changing climate.

Arctic foxes Vulpes lagopus (L.) display a sharp 3- to 5-year fluctuation in population size where lemmings are their main prey. In areas devoid of lemmings, such as Iceland, they do not experience short-term fluctuations. This study focusses on the population dynamics of the arctic fox in Iceland and how it is shaped by its main prey populations. Hunting statistics from 1958-2003 show that the population size of the arctic fox was at a maximum in the 1950s, declined to a minimum in the 1970s, and increased steadily until 2003. Analysis of the arctic fox population size and their prey populations suggests that fox numbers were limited by rock ptarmigan numbers during the decline period. The recovery of the arctic fox population was traced mostly to an increase in goose populations, and favourable climatic conditions as reflected by the Subpolar Gyre. These results underscore the flexibility of a generalist predator and its responses to shifting food resources and climate changes.

Latitudinal Gradient in Otolith Shape among Local Populations of Atlantic Herring (Clupea harengus L.) in Norway.

Otolith shape analysis of Atlantic herring (Clupea harengus) in Norwegian waters shows significant differentiation among fjords and a latitudinal gradient along the coast where neighbouring populations are more similar to each other than to those sampled at larger distances. The otolith shape was obtained using quantitative shape analysis, the outlines were transformed with Wavelet and analysed with multivariate methods. The observed morphological differences are likely to reflect environmental differences but indicate low dispersal among the local herring populations. Otolith shape variation suggests also limited exchange between the local populations and their oceanic counterparts, which could be due to differences in spawning behaviour. Herring from the most northerly location (69°N) in Balsfjord, which is genetically more similar to Pacific herring (Clupea pallasii), differed in otolith shape from all the other populations. Our results suggest that the semi-enclosed systems, where the local populations live and breed, are efficient barriers for dispersal. Otolith shape can thus serve as a marker to identify the origin of herring along the coast of Norway.

ShapeR: an R package to study otolith shape variation among fish populations.

ShapeR is an open source software package that runs on the R platform and is specifically designed to study otolith shape variation among fish populations. The package extends previously described software used for otolith shape analysis by allowing the user to automatically extract closed contour outlines from a large number of images, perform smoothing to eliminate pixel noise, choose from conducting either a Fourier or Wavelet transform to the outlines and visualize the mean shape. The output of the package are independent Fourier or Wavelet coefficients which can be directly imported into a wide range of statistical packages in R. The package might prove useful in studies of any two dimensional objects.

Historical DNA reveals the demographic history of Atlantic cod (Gadus morhua) in medieval and early modern Iceland.

Atlantic cod (Gadus morhua) vertebrae from archaeological sites were used to study the history of the Icelandic Atlantic cod population in the time period of 1500-1990. Specifically, we used coalescence modelling to estimate population size and fluctuations from the sequence diversity at the cytochrome b (cytb) and Pantophysin I (PanI) loci. The models are consistent with an expanding population during the warm medieval period, large historical effective population size (NE), a marked bottleneck event at 1400-1500 and a decrease in NE in early modern times. The model results are corroborated by the reduction of haplotype and nucleotide variation over time and pairwise population distance as a significant portion of nucleotide variation partitioned across the 1550 time mark. The mean age of the historical fished stock is high in medieval times with a truncation in age in early modern times. The population size crash coincides with a period of known cooling in the North Atlantic, and we conclude that the collapse may be related to climate or climate-induced ecosystem change.