Integrating genetic data and population viability analyses for the identification of harbour seal (Phoca vitulina) populations and management units.

Identification of populations and management units is an essential step in the study of natural systems. Still, there is limited consensus regarding how to define populations and management units, and whether genetic methods allow for inference at the relevant spatial and temporal scale. Here, we present a novel approach, integrating genetic, life history and demographic data to identify populations and management units in southern Scandinavian harbour seals. First, 15 microsatellite markers and model- and distance-based genetic clustering methods were used to determine the population genetic structure in harbour seals. Second, we used harbour seal demographic and life history data to conduct population viability analyses (PVAs) in the vortex simulation model in order to determine whether the inferred genetic units could be classified as management units according to Lowe and Allendorf's (Molecular Ecology, 19, 2010, 3038) 'population viability criterion' for demographic independence. The genetic analyses revealed fine-scale population structuring in southern Scandinavian harbour seals and pointed to the existence of several genetic units. The PVAs indicated that the census population size of each of these genetic units was sufficiently large for long-term population viability, and hence that the units could be classified as demographically independent management units. Our study suggests that population genetic inference can offer the same degree of temporal and spatial resolution as 'nongenetic' methods and that the combined use of genetic data and PVAs constitutes a promising approach for delineating populations and management units.

Herbarium specimens reveal a historical shift in phylogeographic structure of common ragweed during native range disturbance.

Invasive plants provide ample opportunity to study evolutionary shifts that occur after introduction to novel environments. However, although genetic characters pre-dating introduction can be important determinants of later success, large-scale investigations of historical genetic structure have not been feasible. Common ragweed (Ambrosia artemisiifolia L.) is an invasive weed native to North America that is known for its allergenic pollen. Palynological records from sediment cores indicate that this species was uncommon before European colonization of North America, and ragweed populations expanded rapidly as settlers deforested the landscape on a massive scale, later becoming an aggressive invasive with populations established globally. Towards a direct comparison of genetic structure now and during intense anthropogenic disturbance of the late 19th century, we sampled 45 natural populations of common ragweed across its native range as well as historical herbarium specimens collected up to 140 years ago. Bayesian clustering analyses of 453 modern and 473 historical samples genotyped at three chloroplast spacer regions and six nuclear microsatellite loci reveal that historical ragweed's spatial genetic structure mirrors both the palaeo-record of Ambrosia pollen deposition and the historical pattern of agricultural density across the landscape. Furthermore, for unknown reasons, this spatial genetic pattern has changed substantially in the intervening years. Following on previous work relating morphology and genetic expression between plants collected from eastern North America and Western Europe, we speculate that the cluster associated with humans' rapid transformation of the landscape is a likely source of these aggressive invasive populations.

Inferring recent historic abundance from current genetic diversity.

Recent historic abundance is an elusive parameter of great importance for conserving endangered species and understanding the pre-anthropogenic state of the biosphere. The number of studies that have used population genetic theory to estimate recent historic abundance from contemporary levels of genetic diversity has grown rapidly over the last two decades. Such assessments often yield unexpectedly large estimates of historic abundance. We review the underlying theory and common practices of estimating recent historic abundance from contemporary genetic diversity, and critically evaluate the potential issues at various estimation steps. A general issue of mismatched spatio-temporal scales between the estimation itself and the objective of the estimation emerged from our assessment; genetic diversity-based estimates of recent historic abundance represent long-term averages, whereas the objective typically is an estimate of recent abundance for a specific population. Currently, the most promising approach to estimate the difference between recent historic and contemporary abundance requires that genetic data be collected from samples of similar spatial and temporal duration. Novel genome-enabled inference methods may be able to utilize additional information of dense genome-wide distributions of markers, such as of identity-by-descent tracts, to infer recent historic abundance from contemporary samples only.

Using virtual reality for anatomical landmark annotation in geometric morphometrics.

To study the shape of objects using geometric morphometrics, landmarks are oftentimes collected digitally from a 3D scanned model. The expert may annotate landmarks using software that visualizes the 3D model on a flat screen, and interaction is achieved with a mouse and a keyboard. However, landmark annotation of a 3D model on a 2D display is a tedious process and potentially introduces error due to the perception and interaction limitations of the flat interface. In addition, digital landmark placement can be more time-consuming than direct annotation on the physical object using a tactile digitizer arm. Since virtual reality (VR) is designed to more closely resemble the real world, we present a VR prototype for annotating landmarks on 3D models. We study the impact of VR on annotation performance by comparing our VR prototype to Stratovan Checkpoint, a commonly used commercial desktop software. We use an experimental setup, where four operators placed six landmarks on six grey seal (Halichoerus grypus) skulls in six trials for both systems. This enables us to investigate multiple sources of measurement error. We analyse both for the configuration and for single landmarks. Our analysis shows that annotation in VR is a promising alternative to desktop annotation. We find that annotation precision is comparable between the two systems, with VR being significantly more precise for one of the landmarks. We do not find evidence that annotation in VR is faster than on the desktop, but it is accurate.

Skull ecomorphological variation of narwhals (Monodon monoceros, Linnaeus 1758) and belugas (Delphinapterus leucas, Pallas 1776) reveals phenotype of their hybrids.

Narwhals and belugas are toothed whales belonging to the Monodontidae. Belugas have a circumpolar Arctic and sub-Artic distribution while narwhals are restricted to the Atlantic Arctic. Their geographical ranges overlap during winter migrations in the Baffin Bay area (Canada/West Greenland) and successful interbreeding may occur. Here, we employed geometric morphometrics on museum specimens to explore the cranium and mandible morphology of a known hybrid (NHMD MCE 1356) and the cranium morphology of a putative hybrid (NHMD 1963.44.1.4) relative to skull morphological variation in the parental species. Specifically, we used 3D models of skulls from 69 belugas, 86 narwhals, and the two known/putative hybrids and 2D left hemi-mandibles from 20 belugas, 64 narwhals and the known hybrid. Skull shape analyses allowed clear discrimination between species. Narwhals are characterised by a relatively short rostrum and wide neurocranium while belugas show a more elongated and narrower cranium. Sexual size dimorphism was detected in narwhals, with males larger than females, but no sexual shape dimorphism was detected in either species (excluding presence/absence of tusks in narwhals). Morphological skull variation was also dependent on different allometric slopes between species and sexes in narwhals. Our analyses showed that the cranium of the known hybrid was phenotypically close to belugas but its 2D hemi-mandible had a narwhal shape and size morphology. Both cranium and mandible were strongly correlated, with the pattern of covariation being similar to belugas. The putative hybrid was a pure male narwhal with extruded teeth. Comparison of genomic DNA supported this result, and stable carbon and nitrogen isotope values suggested that the putative hybrid had a more benthic foraging strategy compared to narwhals. This work demonstrates that although the known hybrid could be discriminated from narwhals and belugas, detection of its affinities with these parental species was dependent on the part of the skull analysed.

Measurement error using a SeeMaLab structured light 3D scanner against a Microscribe 3D digitizer.

BACKGROUND: Geometric morphometrics is a powerful approach to capture and quantify morphological shape variation. Both 3D digitizer arms and structured light surface scanners are portable, easy to use, and relatively cheap, which makes these two capturing devices obvious choices for geometric morphometrics. While digitizer arms have been the "gold standard", benefits of having full 3D models are manifold. We assessed the measurement error and investigate bias associated with the use of an open-source, high-resolution structured light scanner called SeeMaLab against the popular Microscribe 3D digitizer arm. METHODOLOGY: The analyses were based on 22 grey seal (Halichoerus grypus) skulls. 31 fixed anatomical landmarks were annotated both directly using a Microscribe 3D digitizer and on reconstructed 3D digital models created from structured light surface scans. Each skull was scanned twice. Two operators annotated the landmarks, each twice on all the skulls and 3D models, allowing for the investigation of multiple sources of measurement error. We performed multiple Procrustes ANOVAs to compare the two devices in terms of within- and between-operator error, to quantify the measurement error induced by device, to compare between-device error with other sources of variation, and to assess the level of scanning-related error. We investigated the presence of general shape bias due to device and operator. RESULTS: Similar precision was obtained with both devices. If landmarks that were identified as less clearly defined and thus harder to place were omitted, the scanner pipeline would achieve higher precision than the digitizer. Between-operator error was biased and seemed to be smaller when using the scanner pipeline. There were systematic differences between devices, which was mainly driven by landmarks less clearly defined. The factors device, operator and landmark replica were all statistically significant and of similar size, but were minor sources of total shape variation, compared to the biological variation among grey seal skulls. The scanning-related error was small compared to all other error sources. CONCLUSIONS: As the scanner showed precision similar to the digitizer, a scanner should be used if the advantages of obtaining detailed 3D models of a specimen are desired. To obtain high precision, a pre-study should be conducted to identify difficult landmarks. Due to the observed bias, data from different devices and/or operators should not be combined when the expected biological variation is small, without testing the landmarks for repeatability across platforms and operators. For any study necessitating the combination of landmark measurements from different operators, the scanner pipeline will be better suited. The small scanning-related error indicates that by following the same scanning protocol, different operators can be involved in the scanning process without introducing significant error.

Dental pulp exposure, periapical inflammation and suppurative osteomyelitis of the jaws in juvenile Baltic grey seals (Halichoerus grypus grypus) from the late 19th century.

The systematic analysis of museum collections can provide important insights into the dental and skeletal pathology of wild mammals. Here we present a previously unreported type of dental defect and related skull pathology in five juvenile Baltic grey seals that had been collected in the course of a seal culling program along the Danish coast in 1889 and 1890. All five skulls exhibited openings into the pulp cavities at the crown tips of all (four animals) or two (one animal) canines as well as several incisors and (in one animal) also some anterior premolars. The affected teeth showed wide pulp cavities and thin dentin. Pulp exposure had caused infection, inflammation, and finally necrosis of the pulp. As was evidenced by the extensive radiolucency around the roots of the affected teeth, the inflammation had extended from the pulp into the periapical space, leading to apical periodontitis with extensive bone resorption. Further spreading of the inflammation into the surrounding bone regions had then caused suppurative osteomyelitis of the jaws. The postcanine teeth of the pathological individuals typically had dentin of normal thickness and, except for one specimen, did not exhibit pulp exposure. The condition may have been caused by a late onset of secondary and tertiary dentin formation that led to pulp exposure in anterior teeth exposed to intense wear. Future investigations could address a possible genetic causation of the condition in the studied grey seals.

Genomics of host-pathogen interactions: challenges and opportunities across ecological and spatiotemporal scales.

Evolutionary genomics has recently entered a new era in the study of host-pathogen interactions. A variety of novel genomic techniques has transformed the identification, detection and classification of both hosts and pathogens, allowing a greater resolution that helps decipher their underlying dynamics and provides novel insights into their environmental context. Nevertheless, many challenges to a general understanding of host-pathogen interactions remain, in particular in the synthesis and integration of concepts and findings across a variety of systems and different spatiotemporal and ecological scales. In this perspective we aim to highlight some of the commonalities and complexities across diverse studies of host-pathogen interactions, with a focus on ecological, spatiotemporal variation, and the choice of genomic methods used. We performed a quantitative review of recent literature to investigate links, patterns and potential tradeoffs between the complexity of genomic, ecological and spatiotemporal scales undertaken in individual host-pathogen studies. We found that the majority of studies used whole genome resolution to address their research objectives across a broad range of ecological scales, especially when focusing on the pathogen side of the interaction. Nevertheless, genomic studies conducted in a complex spatiotemporal context are currently rare in the literature. Because processes of host-pathogen interactions can be understood at multiple scales, from molecular-, cellular-, and physiological-scales to the levels of populations and ecosystems, we conclude that a major obstacle for synthesis across diverse host-pathogen systems is that data are collected on widely diverging scales with different degrees of resolution. This disparity not only hampers effective infrastructural organization of the data but also data granularity and accessibility. Comprehensive metadata deposited in association with genomic data in easily accessible databases will allow greater inference across systems in the future, especially when combined with open data standards and practices. The standardization and comparability of such data will facilitate early detection of emerging infectious diseases as well as studies of the impact of anthropogenic stressors, such as climate change, on disease dynamics in humans and wildlife.

Dolphin Morbillivirus in a Fin Whale (Balaenoptera physalus) in Denmark, 2016.

We studied the etiology of encephalitis in a fin whale (Balaenoptera physalus) that stranded in 2016 on the coast of Denmark. Dolphin morbillivirus (DMV) was detected in the brain and other organs. Phylogenetics showed close relation to DMV isolated from a striped dolphin (Stenella coeruleoalba) from Spain in 2012.

What are the toxicological effects of mercury in Arctic biota?

This review critically evaluates the available mercury (Hg) data in Arctic marine biota and the Inuit population against toxicity threshold values. In particular marine top predators exhibit concentrations of mercury in their tissues and organs that are believed to exceed thresholds for biological effects. Species whose concentrations exceed threshold values include the polar bears (Ursus maritimus), beluga whale (Delphinapterus leucas), pilot whale (Globicephala melas), hooded seal (Cystophora cristata), a few seabird species, and landlocked Arctic char (Salvelinus alpinus). Toothed whales appear to be one of the most vulnerable groups, with high concentrations of mercury recorded in brain tissue with associated signs of neurochemical effects. Evidence of increasing concentrations in mercury in some biota in Arctic Canada and Greenland is therefore a concern with respect to ecosystem health.

Time trends of mercury in feathers of West Greenland birds of prey during 1851-2003.

Temporal trends of mercury (Hg) in West Greenland gyrfalcons, peregrine falcons, and white-tailed eagles were determined over 150 years from 1851 to 2003. Hg was measured in the fifth primary feather. Results showed that Hg increased in the order gyrfalcon (lowest) < peregrine falcon (intermediate) < white-tailed eagle (highest). All species showed significant age accumulations, which were taken into account in the temporal trend analysis. Of eight time trend analyses (three species and three age groups of which one was missing), seven showed an increase in primary feather concentrations. Of these, four were significant at the 5% level, two were close to being significant, and one was not significant. The linear regressions of which three out of four showed significant increases were for juvenile and immature gyrfalcon and juvenile peregrine falcon, which covered only periods prior to 1960, owing to limited data from the last half-century. The two sample comparisons of Hg 10-year medians for adult peregrine falcons and juvenile and adult white-tailed eagles indicated a continued increase during recent decades. However, low levels of Hg in a few recent collections among gyrfalcons and peregrines could indicate a change in the increasing trend.