Over the last few centuries, millions of plant specimens have accumulated within herbaria and biocultural collections. These include type specimens, used to define species, and populations that are rare, extinct or difficult to access. They therefore represent a considerable resource for a broad range of scientific uses. However, collections degrade over time, and become increasingly difficult to characterise their genetic signatures, even considering exponential advancements in sequencing technologies. Here, we tested the genotyping performance on highly degraded samples using a commonly used high-throughput sequencing (HtS) technique, genome skimming, against a recent alternative target capture kit, the universal set Angiosperm-353. We performed phylogenomic analyses of modern leaf and historical barks of Cinchona, including 23 historical barks and six fresh leaf specimens. DNA within historical barks is highly degraded, therefore a customised DNA extraction method was developed before library preparation. We show that sample degradation over time directly impacted the quantity and quality of the data produced by both methodologies (in terms of reads mapped to the references). However, we find that both approaches generate enough data to infer phylogenetic relationships, even between highly degraded specimens that are over 230 years old. Within historical barks, the target capture kit is more advantageous than genome skimming in profiling Cinchona species since it was possible to retrieve nuclear and plastid data to infer phylogenies. This study showcases the value of historical samples in genetic studies and paves the way for further experiments across different taxonomic groups with varying levels of genetic variation or hybridisation.

Using whole mitochondrial DNA sequences from 89 white-tailed eagles (Haliaeetus albicilla) sampled from Iceland, Greenland, Norway, Denmark and Estonia between 1990-2018, we investigate the mitogenomic variation within and between countries. We show that there is a substantial population differentiation between the countries, reflecting similar major phylogeographic patterns obtained previously for the control region of the mitochondria, which suggested two main refugia during the last glacial period of Ice Age. Distinct mitogenomic lineages are observed within countries which divergence times exceeds the time since last glacial period of Ice Age ended. The lineages appear to have been maintained by natural selection. An excess of segregating amino acids in comparison with number of fixations, as revealed by the neutrality index suggests a load of deleterious mutations. The maintenance of mtDNA lineages within countries inflates our estimates of effective national population sizes and the times of their divergence.

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: I) isolation of refugia during the last glacial period 110-115,000 years ago, II) population divergence following the colonization of the deglaciated areas ~10,000 years ago, III) human population expansion, which led to the settlement in Iceland ~1,100 years ago, and IV) human persecution and exposure to toxic pollutants during the last two centuries.

Aim: We have studied population genetic change through time in the Northern dragonhead, Dracocephalum ruyschiana (Lamiaceae); a plant species that has experienced a drastic population decline and habitat loss in Europe. We aimed at adding a historic level to the monitoring of dragonhead by testing a microfluidic SNP array approach on herbarium specimens up to 200 years old and comparing the genomic results with that of modern populations in Norway. We also aimed to gain a more holistic species knowledge to guide monitoring efforts by combining herbarium genomics with ecological niche modelling (ENM). Location: Europe (mainly Norway) Methods: We have applied a microfluidic array consisting of 96 SNP markers on 130 herbarium specimens collected from 1820 to 2008. Obtained genotype data were compared with SNP data from modern samples using various population genetic analyses. We used sample metadata and observational records to model the species’ environmental niche. Results: The SNP array successfully genotyped all included herbarium specimens but was less capable of capturing diversity outside of Norway, which was genetically highly divergent from the Norwegian dragonheads. The historic-modern comparison revealed similar genetic structure in space and limited change through time in Norway. The ENM suggests that dragonhead has not fully achieved its potential distribution in Norway, which is anchored in warmer and drier regions, including areas where it does not occur today. Main conclusions: With the appropriate design procedures, the SNP array technology is promising for genotyping old herbarium specimens; an invaluable source of information from the past. We found no signs of the severe reduction in population size in our temporal genomic data of Norwegian dragonhead. Regardless, the regional populations in Norway are genetically divergent, both from each other and more so from populations outside of Norway, rendering continued protection of all existing populations of the species relevant.

Using whole genome shotgun sequences from 92 white-tailed eagles (Haliaeetus albicilla) sampled from Greenland, Iceland, Norway, Denmark, Estonia, and Turkey between 1885-1950 and after 1990, we investigate the genomic variation within countries over time, and between countries. Clear genetic differentiation is observed between samples from the different countries, with the largest differences between the island and mainland populations, and indications that the island populations share the most recent ancestry with the Norwegian population. We find signs of strong inbreeding in the island populations. Further, temporal differences are observed in some populations, for example, replacement of the Danish gene pool following its population’s extinction in the early 20th century, as well as a change in the genetic diversity of the Icelandic population following a severe bottleneck during the last century, all of which could warrant a further conservation effort in Iceland. More generally, all populations show a decline in effective population size, which may have been shaped by I) distinct refugia during the last glacial period, II) population divergence following the colonization of the deglaciated areas ~10,000 years ago, III) human population expansion and e.g., settlement in Iceland ~1,100 years ago, and IV) human persecution and toxic pollutants during the last two centuries.

Runs of homozygosity (ROH) occur when offspring inherit haplotypes that are identical by descent from each parent. Length distributions of ROH are informative about population history; specifically the probability of inbreeding mediated by mating system and/or population demography. Here, we investigate whether variation in killer whale (Orcinus orca) demographic history is reflected in genome-wide heterozygosity and ROH length distributions, using a global dataset of 26 genomes representative of geographic and ecotypic variation in this species, and two F1 admixed individuals with Pacific-Atlantic parentage. We first reconstruct demographic history for each population as changes in effective population size through time using the pairwise sequential Markovian coalescent (PSMC) method. We find a subset of populations declined in effective population size during the Late Pleistocene, while others had more stable demography. Genomes inferred to have undergone ancestral declines in effective population size, were autozygous at hundreds of short ROH (<1Mb), reflecting high background relatedness due to coalescence of haplotypes deep within the pedigree. In contrast, longer and therefore younger ROH (>1.5 Mb) were found in low latitude populations and populations of known conservation concern, including a Scottish population, for which 37.8% of the autosomes comprised of ROH >1.5 Mb in length.