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

Stable carbon (δ¹³C) and nitrogen (δ¹⁵N) 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 δ¹³C and δ¹⁵N 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. δ¹³C and δ¹⁵N were measured, and intra-individual bone and dentine isotopic compositions were used to calculate correction factors for each species. We detected differences in δ¹³C and δ¹⁵N. In both narwhals and belugas, we found lower average δ¹³C and δ¹⁵N 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.

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 km². 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.

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.

Ocean-wide genomic variation in Gray's beaked whales, Mesoplodon grayi

The deep oceans of the Southern Hemisphere are home to several elusive and poorly studied marine megafauna. In the absence of robust observational data for these species, genetic data can aid inferences on population connectivity, demography and ecology. A previous investigation of genetic diversity and population structure in Gray's beaked whale (Mesoplodon grayi) from Western Australia and New Zealand found high levels of mtDNA diversity, no geographic structure and stable demographic history. To further investigate phylogeographic and demographic patterns across their range, we generated complete mitochondrial and partial nuclear genomes of 16 of the individuals previously analysed and included additional samples from South Africa (n = 2) and South Australia (n = 4), greatly expanding the spatial range of genomic data for the species. Gray's beaked whales are highly elusive and rarely observed, and our data represents a unique and geographically broad dataset. We find relatively high levels of diversity in the mitochondrial genome, despite an absence of population structure at the mitochondrial and nuclear level. Demographic analyses suggest these whales existed at stable levels over at least the past 1.1 million years, with an approximately twofold increase in female effective population size approximately 250 thousand years ago, coinciding with a period of increased Southern Ocean productivity, sea surface temperature and a potential expansion of suitable habitat. Our results suggest that Gray's beaked whales are likely to be resilient to near-future ecosystem changes, facilitating their conservation. Our study demonstrates the utility of low-effort shotgun sequencing in providing ecological information on highly elusive species.

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 δ ¹⁵N 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.

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.

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.