Evolutionary history and seascape genomics of Harbour porpoises (Phocoena phocoena) across environmental gradients in the North Atlantic and adjacent waters.

The Harbour porpoise (Phocoena phocoena) is a highly mobile cetacean species primarily occurring in coastal and shelf waters across the Northern hemisphere. It inhabits heterogeneous seascapes broadly varying in salinity and temperature. Here, we produced 74 whole genomes at intermediate coverage to study Harbour porpoise's evolutionary history and investigate the role of local adaptation in the diversification into subspecies and populations. We identified ~6 million high quality SNPs sampled at eight localities across the North Atlantic and adjacent waters, which we used for population structure, demographic and genotype-environment association analyses. Our results suggest a genetic differentiation between three subspecies (P.p. relicta, P.p. phocoena and P.p. meridionalis), and three distinct populations within P.p. phocoena: Atlantic, Belt Sea and Proper Baltic Sea. Effective population size and Tajima's D suggest population contraction in Black Sea and Iberian porpoises, but expansion in the P.p. phocoena populations. Phylogenetic trees indicate post-glacial colonization from a southern refugium. Genotype-environment association analysis identified salinity as major driver in genomic variation and we identified candidate genes putatively underlying adaptation to different salinity. Our study highlights the value of whole genome resequencing to unravel subtle population structure in highly mobile species, shows how strong environmental gradients and local adaptation may lead to population differentiation, and how neutral and adaptive markers can give different perspectives on population subdivision. The results have great conservation implications as we found inbreeding and low genetic diversity in the endangered Black Sea subspecies and identified the critically endangered Proper Baltic Sea porpoises as a separate population.

Southern introgression increases adaptive immune gene variability in northern range margin populations of Fire-bellied toad.

Northern range margin populations of the European fire-bellied toad (Bombina bombina) have rapidly declined during recent decades. Extensive agricultural land use has fragmented the landscape, leading to habitat disruption and loss, as well as eutrophication of ponds. In Northern Germany (Schleswig-Holstein) and Southern Sweden (Skåne), this population decline resulted in decreased gene flow from surrounding populations, low genetic diversity, and a putative reduction in adaptive potential, leaving populations vulnerable to future environmental and climatic changes. Previous studies using mitochondrial control region and nuclear transcriptome-wide SNP data detected introgressive hybridization in multiple northern B. bombina populations after unreported release of toads from Austria. Here, we determine the impact of this introgression by comparing the body conditions (proxy for fitness) of introgressed and nonintrogressed populations and the genetic consequences in two candidate genes for putative local adaptation (the MHC II gene as part of the adaptive immune system and the stress response gene HSP70 kDa). We detected regional differences in body condition and observed significantly elevated levels of within individual MHC allele counts in introgressed Swedish populations, associated with a tendency toward higher body weight, relative to regional nonintrogressed populations. These differences were not observed among introgressed and nonintrogressed German populations. Genetic diversity in both MHC and HSP was generally lower in northern than Austrian populations. Our study sheds light on the potential benefits of translocations of more distantly related conspecifics as a means to increase adaptive genetic variability and fitness of genetically depauperate range margin populations without distortion of local adaptation.

Supporting evidence for PCB pollution threatening global killer whale population.

A recent Science report predicted the global killer whale population to collapse due to PCB pollution. Here we present empirical evidence, which supports and extends the reports' statement. In 2016, a neonate male killer whale stranded on the German island of Sylt. Neonatal attributes indicated an age of at least 3 days. The stomach contained ∼20 mL milk residue and no pathologies explaining the cause of death could be detected. Blubber samples presenting low lipid concentrations were analysed for persistent organic pollutants. Skin samples were collected for genotyping of the mitochondrial control region. The blubber PCB concentrations were very high [SPCBs, 225 mg/kg lipid weight (lw)], largely exceeding the PCB toxicity thresholds reported for the onset of immunosuppression [9 mg/kg lw ∑PCB] and for severe reproductive impairment [41 mg/kg lw ∑PCB] reported for marine mammals. Additionally, this individual showed equally high concentrations in p,p'-DDE [226 mg/kg lw], PBDEs [5 mg/kg lw] and liver mercury levels [1.1 µg/g dry weight dw]. These results suggest a high placental transfer of pollutants from mother to foetus. Consequently, blubber and plasma PCB concentrations and calf mortality rates are both high in primiparous females. With such high pollutant levels, this neonate had poor prerequisites for survival. The neonate belonged to Ecotype I (generalist feeder) and carried the mitochondrial haplotype 35 present in about 16% of the North Atlantic killer whale from or close to the North Sea. The relevance of this data becomes apparent in the UK West Coast Community, the UK's only residentorca population, which is currently composed of only eight individuals (each four males and females) and no calves have been reported over the last 19 years.Despite worldwide regulations, PCBs persist in the environment and remain a severe concern for killer whale populations, placing calves at high risk due to the mother-offspring PCB-transfer resulting in a high toxicological burden of the neonates.

High-quality whole-genome sequence of an abundant Holarctic odontocete, the harbour porpoise (Phocoena phocoena).

The harbour porpoise (Phocoena phocoena) is a highly mobile cetacean found across the Northern hemisphere. It occurs in coastal waters and inhabits basins that vary broadly in salinity, temperature and food availability. These diverse habitats could drive subtle differentiation among populations, but examination of this would be best conducted with a robust reference genome. Here, we report the first harbour porpoise genome, assembled de novo from an individual originating in the Kattegat Sea (Sweden). The genome is one of the most complete cetacean genomes currently available, with a total size of 2.39 Gb and 50% of the total length found in just 34 scaffolds. Using 122 of the longest scaffolds, we were able to show high levels of synteny with the genome of the domestic cattle (Bos taurus). Our draft annotation comprises 22,154 predicted genes, which we further annotated through matches to the NCBI nucleotide database, GO categorization and motif prediction. Within the predicted genes, we have confirmed the presence of >20 genes or gene families that have been associated with adaptive evolution in other cetaceans. Overall, this genome assembly and draft annotation represent a crucial addition to the genomic resources currently available for the study of porpoises and Phocoenidae evolution, phylogeny and conservation.

Targeted resequencing of coding DNA sequences for SNP discovery in nonmodel species.

Targeted capture coupled with high-throughput sequencing can be used to gain information about nuclear sequence variation at hundreds to thousands of loci. Divergent reference capture makes use of molecular data of one species to enrich target loci in other (related) species. This is particularly valuable for nonmodel organisms, for which often no a priori knowledge exists regarding these loci. Here, we have used targeted capture to obtain data for 809 nuclear coding DNA sequences (CDS) in a nonmodel organism, the Eurasian lynx Lynx lynx, using baits designed with the help of the published genome of a related model organism (the domestic cat Felis catus). Using this approach, we were able to survey intraspecific variation at hundreds of nuclear loci in L. lynx across the species' European range. A large set of biallelic candidate SNPs was then evaluated using a high-throughput SNP genotyping platform (Fluidigm), which we then reduced to a final 96 SNP-panel based on assay performance and reliability; validation was carried out with 100 additional Eurasian lynx samples not included in the SNP discovery phase. The 96 SNP-panel developed from CDS performed very successfully in the identification of individuals and in population genetic structure inference (including the assignment of individuals to their source population). In keeping with recent studies, our results show that genic SNPs can be valuable for genetic monitoring of wildlife species.

Inter-individual differences in contamination profiles as tracer of social group association in stranded sperm whales.

Ecological and physiological factors lead to different contamination patterns in individual marine mammals. The objective of the present study was to assess whether variations in contamination profiles are indicative of social structures of young male sperm whales as they might reflect a variation in feeding preferences and/or in utilized feeding grounds. We used a total of 61 variables associated with organic compounds and trace element concentrations measured in muscle, liver, kidney and blubber gained from 24 sperm whales that stranded in the North Sea in January and February 2016. Combining contaminant and genetic data, there is evidence for at least two cohorts with different origin among these stranded sperm whales; one from the Canary Island region and one from the northern part of the Atlantic. While genetic data unravel relatedness and kinship, contamination data integrate over areas, where animals occured during their lifetime. Especially in long-lived animals with a large migratory potential, as sperm whales, contamination data may carry highly relevant information about aggregation through time and space.