Combining δ13C and δ15N from bone and dentine in marine mammal palaeoecological research: insights from toothed whales.

Stable carbon (δ13C) and nitrogen (δ15N) isotopic compositions of bone and dentine collagen extracted from museum specimens have been widely used to study the paleoecology of past populations. Due to possible systematic differences in stable isotope values between bone and dentine, dentine values need to be transformed into bone-collagen equivalent using a correction factor to allow comparisons between the two collagen sources. Here, we provide correction factors to transform dentine δ13C and δ15N values into bone-collagen equivalent for two toothed whales: narwhal and beluga. We sampled bone and dentine from the skulls of 11 narwhals and 26 belugas. In narwhals, dentine was sampled from tusk and embedded tooth; in belugas, dentine was sampled from tooth. δ13C and δ15N were measured, and intra-individual bone and dentine isotopic compositions were used to calculate correction factors for each species. We detected differences in δ13C and δ15N. In both narwhals and belugas, we found lower average δ13C and δ15N in bone compared with dentine. The correction factors provided by the study enable the combined analysis of stable isotope data from bone and dentine in these species.

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.

How low can you go? Introducing SeXY: sex identification from low-quantity sequencing data despite lacking assembled sex chromosomes.

Accurate sex identification is crucial for elucidating the biology of a species. In the absence of directly observable sexual characteristics, sex identification of wild fauna can be challenging, if not impossible. Molecular sexing offers a powerful alternative to morphological sexing approaches. Here, we present SeXY, a novel sex-identification pipeline, for very low-coverage shotgun sequencing data from a single individual. SeXY was designed to utilize low-effort screening data for sex identification and does not require a conspecific sex-chromosome assembly as reference. We assess the accuracy of our pipeline to data quantity by downsampling sequencing data from 100,000 to 1000 mapped reads and to reference genome selection by mapping to a variety of reference genomes of various qualities and phylogenetic distance. We show that our method is 100% accurate when mapping to a high-quality (highly contiguous N50 > 30 Mb) conspecific genome, even down to 1000 mapped reads. For lower-quality reference assemblies (N50 < 30 Mb), our method is 100% accurate with 50,000 mapped reads, regardless of reference assembly quality or phylogenetic distance. The SeXY pipeline provides several advantages over previously implemented methods; SeXY (i) requires sequencing data from only a single individual, (ii) does not require assembled conspecific sex chromosomes, or even a conspecific reference assembly, (iii) takes into account variation in coverage across the genome, and (iv) is accurate with only 1000 mapped reads in many cases.

Complete mitochondrial genome of the giant root-rat (Tachyoryctes macrocephalus).

The endangered giant root-rat (Tachyoryctes macrocephalus, also known as giant mole rat) is a fossorial rodent endemic to the afro-alpine grasslands of the Bale Mountains in Ethiopia. The species is an important ecosystem engineer with the majority of the global population found within 1000 km2. Here, we present the first complete mitochondrial genome of the giant root-rat and the genus Tachyoryctes, recovered using shotgun sequencing and iterative mapping. A phylogenetic analysis including 15 other representatives of the family Spalacidae placed Tachyoryctes as sister genus to Rhizomys with high support. This position is in accordance with a recent study revealing the topology of the Spalacidae family. The full mitochondrial genome of the giant root-rat presents an important resource for further population genetic studies.

Population-specific sex and size variation in long-term foraging ecology of belugas and narwhals.

Intraspecific variation in resource use by individuals of different age, sex or size may reflect differing energetic requirements and physiological constraints. Males and females often show differences in diet owing to sexual size dimorphism, different life histories and/or habitat use. Here, we investigate how sex and size influence the long-term foraging ecology of belugas and narwhals in Greenland, using stable isotopes of carbon and nitrogen from bone collagen. We show that males have a higher trophic level and a larger ecological niche than females in West Greenland belugas and in East Greenland narwhals. In addition, for these two populations, we find that δ 15N increases with size, particularly in males. We hypothesize that sexual size dimorphism together with strong maternal investment drive these differences. By contrast, we find no differences in foraging ecology between sexes in West Greenland narwhals and observe no influence of size on trophic level. This may reflect the influence of interspecific competition in West Greenland, where the distributions of belugas and narwhals overlap, and/or geographical resource partitioning among different summer aggregations of narwhals. Our results suggest that sex and size variations in diet are population dependent, and probably the result of varying ecological interactions.

Circumpolar phylogeography and demographic history of beluga whales reflect past climatic fluctuations.

Several Arctic marine mammal species are predicted to be negatively impacted by rapid sea ice loss associated with ongoing ocean warming. However, consequences for Arctic whales remain uncertain. To investigate how Arctic whales responded to past climatic fluctuations, we analysed 206 mitochondrial genomes from beluga whales (Delphinapterus leucas) sampled across their circumpolar range, and four nuclear genomes, covering both the Atlantic and the Pacific Arctic region. We found four well-differentiated mitochondrial lineages, which were established before the onset of the last glacial expansion ~110 thousand years ago. Our findings suggested these lineages diverged in allopatry, reflecting isolation of populations during glacial periods when the Arctic sea-shelf was covered by multiyear sea ice. Subsequent population expansion and secondary contact between the Atlantic and Pacific Oceans shaped the current geographic distribution of lineages, and may have facilitated mitochondrial introgression. Our demographic reconstructions based on both mitochondrial and nuclear genomes showed markedly lower population sizes during the Last Glacial Maximum (LGM) compared to the preceding Eemian and current Holocene interglacial periods. Habitat modelling similarly revealed less suitable habitat during the LGM (glacial) than at present (interglacial). Together, our findings suggested the association between climate, population size, and available habitat in belugas. Forecasts for year 2100 showed that beluga habitat will decrease and shift northwards as oceans continue to warm, putatively leading to population declines in some beluga populations. Finally, we identified vulnerable populations which, if extirpated as a consequence of ocean warming, will lead to a substantial decline of species-wide haplotype diversity.

Some like it cold: Temperature-dependent habitat selection by narwhals.

The narwhal (Monodon monoceros) is a high-Arctic species inhabiting areas that are experiencing increases in sea temperatures, which together with reduction in sea ice are expected to modify the niches of several Arctic marine apex predators. The Scoresby Sound fjord complex in East Greenland is the summer residence for an isolated population of narwhals. The movements of 12 whales instrumented with Fastloc-GPS transmitters were studied during summer in Scoresby Sound and at their offshore winter ground in 2017-2019. An additional four narwhals provided detailed hydrographic profiles on both summer and winter grounds. Data on diving of the whales were obtained from 20 satellite-linked time-depth recorders and 16 Acousonde™ recorders that also provided information on the temperature and depth of buzzes. In summer, the foraging whales targeted depths between 300 and 850 m where the preferred areas visited by the whales had temperatures ranging between 0.6 and 1.5°C (mean = 1.1°C, SD = 0.22). The highest probability of buzzing activity during summer was at a temperature of 0.7°C and at depths > 300 m. The whales targeted similar depths at their offshore winter ground where the temperature was slightly higher (range: 0.7-1.7°C, mean = 1.3°C, SD = 0.29). Both the probability of buzzing events and the spatial distribution of the whales in both seasons demonstrated a preferential selection of cold water. This was particularly pronounced in winter where cold coastal water was selected and warm Atlantic water farther offshore was avoided. It is unknown if the small temperature niche of whales while feeding is because prey is concentrated at these temperature gradients and is easier to capture at low temperatures, or because there are limitations in the thermoregulation of the whales. In any case, the small niche requirements together with their strong site fidelity emphasize the sensitivity of narwhals to changes in the thermal characteristics of their habitats.

Influence of past climate change on phylogeography and demographic history of narwhals, Monodon monoceros.

The Arctic is warming at an unprecedented rate, with unknown consequences for endemic fauna. However, Earth has experienced severe climatic oscillations in the past, and understanding how species responded to them might provide insight into their resilience to near-future climatic predictions. Little is known about the responses of Arctic marine mammals to past climatic shifts, but narwhals (Monodon monoceros) are considered one of the endemic Arctic species most vulnerable to environmental change. Here, we analyse 121 complete mitochondrial genomes from narwhals sampled across their range and use them in combination with species distribution models to elucidate the influence of past and ongoing climatic shifts on their population structure and demographic history. We find low levels of genetic diversity and limited geographic structuring of genetic clades. We show that narwhals experienced a long-term low effective population size, which increased after the Last Glacial Maximum, when the amount of suitable habitat expanded. Similar post-glacial habitat release has been a key driver of population size expansion of other polar marine predators. Our analyses indicate that habitat availability has been critical to the success of narwhals, raising concerns for their fate in an increasingly warming Arctic.

Hybridization between two high Arctic cetaceans confirmed by genomic analysis.

In 1990, a skull from a morphologically unusual Monodontid was found in West Greenland and collected for the Natural History Museum of Denmark, University of Copenhagen. From its intermediate morphology, the skull was hypothesized to be a beluga/narwhal hybrid. If confirmed, the specimen would, to our knowledge, represent the sole evidence of hybridization between the only two toothed whale species endemic to the Arctic. Here we present genome-wide DNA sequence data from the specimen and investigate its origin using a genomic reference panel of eight belugas and eight narwhals. Our analyses reveal that the specimen is a male, first-generation hybrid between a female narwhal and a male beluga. We use stable carbon and nitrogen isotope analysis to investigate the dietary niche of the hybrid and find a higher δ13C value than in both belugas and narwhals, suggesting a foraging strategy unlike either parental species. These results further our understanding of the interaction between belugas and narwhals, and underscore the importance of natural history collections in monitoring changes in biodiversity. In addition, our study exemplifies how recent major advances in population genomic analyses using genotype likelihoods can provide key biological and ecological insights from low-coverage data (down to 0.05x).

Mind the gut: genomic insights to population divergence and gut microbial composition of two marine keystone species.

BACKGROUND: Deciphering the mechanisms governing population genetic divergence and local adaptation across heterogeneous environments is a central theme in marine ecology and conservation. While population divergence and ecological adaptive potential are classically viewed at the genetic level, it has recently been argued that their microbiomes may also contribute to population genetic divergence. We explored whether this might be plausible along the well-described environmental gradient of the Baltic Sea in two species of sand lance (Ammodytes tobianus and Hyperoplus lanceolatus). Specifically, we assessed both their population genetic and gut microbial composition variation and investigated not only which environmental parameters correlate with the observed variation, but whether host genome also correlates with microbiome variation. RESULTS: We found a clear genetic structure separating the high-salinity North Sea from the low-salinity Baltic Sea sand lances. The observed genetic divergence was not simply a function of isolation by distance, but correlated with environmental parameters, such as salinity, sea surface temperature, and, in the case of A. tobianus, possibly water microbiota. Furthermore, we detected two distinct genetic groups in Baltic A. tobianus that might represent sympatric spawning types. Investigation of possible drivers of gut microbiome composition variation revealed that host species identity was significantly correlated with the microbial community composition of the gut. A potential influence of host genetic factors on gut microbiome composition was further confirmed by the results of a constrained analysis of principal coordinates. The host genetic component was among the parameters that best explain observed variation in gut microbiome composition. CONCLUSIONS: Our findings have relevance for the population structure of two commercial species but also provide insights into potentially relevant genomic and microbial factors with regards to sand lance adaptation across the North Sea-Baltic Sea environmental gradient. Furthermore, our findings support the hypothesis that host genetics may play a role in regulating the gut microbiome at both the interspecific and intraspecific levels. As sequencing costs continue to drop, we anticipate that future studies that include full genome and microbiome sequencing will be able to explore the full relationship and its potential adaptive implications for these species.

Mitochondrial genome divergence between beluga whales in Baffin Bay and the Sea of Okhotsk.

The beluga whale is one of three endemic Arctic whales. The species is philopatric, and its migration patterns are passed from mother to calf. Management of the species is informed by the levels of genetic structuring among summer aggregation sites based on mitochondrial D-Loop data. To assess the levels of differentiation across the entire mitochondrial genome within belugas, we present a comparison between the first two complete mitochondrial genomes from opposite sides of their distribution range: Baffin Bay and the Russian Far East. Our analyses reveal that additional phylogenetic insights can be gained from expanding the genetic region analyzed. Further, we estimate the divergence time between the two mitochondrial genomes to be 0.469 MYA.

Genomic population structure of freshwater-resident and anadromous ide (Leuciscus idus) in north-western Europe.

Climate change experts largely agree that future climate change and associated rises in oceanic water levels over the upcoming decades, will affect marine salinity levels. The subsequent effects on fish communities in estuarine ecosystems however, are less clear. One species that is likely to become increasingly affected by changes in salinity is the ide (Leuciscus idus). The ide is a stenohaline freshwater fish that primarily inhabits rivers, with frequent anadromous behavior when sea salinity does not exceed 15%. Unlike most other anadromous Baltic Sea fish species, the ide has yet to be subjected to large-scale stocking programs, and thus provides an excellent opportunity for studying the natural population structure across the current salinity gradient in the Danish Belts. To explore this, we used Genotyping-by-Sequencing to determine genomic population structure of both freshwater resident and anadromous ide populations in the western Baltic Sea region, and relate the results to the current salinity gradient and the demographic history of ide in the region. The sample sites separate into four clusters, with all anadromous populations in one cluster and the freshwater resident populations in the remaining three. Results demonstrate high level of differentiation between sites hosting freshwater resident populations, but little differentiation among anadromous populations. Thus ide exhibit the genomic population structure of both a typical freshwater species, and a typical anadromous species. In addition to providing a first insight into the population structure of north-western European ide, our data also (1) provide indications of a single illegal introduction by man; (2) suggest limited genetic effects of heavy pollution in the past; and (3) indicate possible historical anadromous behavior in a now isolated freshwater population.