Search

Stable carbon (

Environmental variation across the distribution of wild species can lead to local adaptations. The Baltic Sea was formed when the Fenno-Scandian ice sheet retreated around 12 thousand years ago, creating a new brackish water habitat colonised by both marine and freshwater fish species. The European perch (Perca fluviatilis) is a predatory freshwater fish with a large geographical distribution across Eurasia, where it inhabits a wide range of environmental niches. In the Baltic Sea region it has even developed a specialised brackish water phenotype that can tolerate environmental salinity levels, which are lethal to the ancestral freshwater phenotype. However, very little is known about the colonisation history and underlying genomic mechanisms facilitating the colonisation and adaptation of perch to the Baltic Sea. Here, we use Genotyping-By-Sequencing data from six freshwater and six brackish water localities to disclose the evolutionary relationship between the freshwater and brackish water phenotype. Our results show that the brackish water perch phenotype occurs in multiple distinct genetic clusters. We find that gene flow between brackish water phenotypes with full access to the sea likely led to lower levels of differentiation and higher diversity than in freshwater phenotypes. Selection analyses suggest that genomic adaptation played a role in the colonisation of the Baltic Sea and that the top three regions under selection harbour salinity tolerance genes. We also find a link between the historic salinity of the Baltic Sea and the demographic history of the brackish water phenotypes and go on to discuss the implications of our findings for management of brackish water perch in the western Baltic sea.HighlightsGBS data from 12 perch populations, six with brackish and six with freshwater originColonisation history and differentiated gene flow shaped the current population structureThe brackish water ecotype was found in all three major genetic cladesTop three regions under selection harboured salinity tolerance genesSalinity influenced Ne during the formation of the Baltic SeaGraphical abstract

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.

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(2). 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.

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., Narwhals avoid areas with warm temperatures. Both their prey capturing and diving behavior show their clear preference for cold habitats. They preferentially select habitats where the temperature ranges between 0 and 1°C. Photo: Carsten Egevang

Stable carbon (δ13C) and nitrogen (δ15N) isotope 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 tooth dentine from the skulls of 11 narwhals and 26 belugas. In narwhals, dentine was sampled from tusk and embedded tooth; in beluga, dentine was sampled from tooth. δ13C and δ15N were measured and intraindividual bone and dentine isotopic compositions were used to calculate correction factors for each species.We detected differences in δ13C and δ15N. In narwhals, we found (i) lower average δ13C and δ15N in bone compared with dentine; (ii) no difference in dentine δ13C between tusk and embedded tooth. For belugas, we also detected lower δ13C and δ15N in bone compared with tooth dentine. The correction factors provided by the study enable the combined analysis of stable isotope data from bone and dentine in these species.

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 utilise 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 1,000 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 1,000 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 1,000 mapped reads in many cases.

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.

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 δ 15 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.

The fungal pathogenNeonectria neomacrosporais of increasing concern in Europe where, within the last decade, it has caused substantial damage to forest stands and ornamental trees of the genusAbies(Mill.). Using whole-genome sequencing of a comprehensive collection of isolates, we show the extent of three major clades withinN. neomacrospora, which most likely diverged around the end of the last Ice Age. We find it likely that the current European epidemic ofN. neomacrosporawas founded from a population belonging to the east North American clade. All European isolates (1957-2019) had a common evolutionary history, but substantial and asymmetrical gene flow from the larger American source population could be detected. The European population shows multiple signs of having gone through a bottleneck and subsequent population expansion.

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.

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.

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).

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.

In the 1970s, a flood-protection system dramatically changed a large part of the coastal environment of Køge Bugt, a bay in the western Baltic Sea, from open coast to a brackish lagoon habitat. An anadromous stock of European perch Perca fluviatilis seems to have benefitted from this change, but details about their spawning behavior remain unknown. We used oviductinserted acoustic transmitters to reveal the pre-spawning behavior and spawning sites of this population. Thirteen female perch were caught in the lower stream basin of St. Vejle Å, and were tagged with acoustic transmitters inserted through the oviduct. The fish were tracked from March 2 to May 24, 2012 with both passive and active telemetry systems. The pre-spawning behavior involved short trips between the stream and adjacent lagoons. Twelve of the 13 transmitters (92%) were expulsed during spawning, providing for the first time a strong proof of concept of oviduct-inserted acoustic transmitters in brackish and marine fish spawning studies. The transmitter expulsions were validated using an egg map, which was based on visual observations of perch egg-strands, and 11 of the 12 expulsed transmitters (92%) were located in areas with eggs. Many fish spawned in the brackish water with salinities up to 9.6 PSU. These salinities are higher than those previously observed for European perch spawning in the wild, and call for further investigations of salinity tolerance in perch eggs

Estuaries are important spawning areas for many freshwater and marine fishes. This is also the case for Baltic populations of European perch Perca fluviatilis, an important commercial species in most of Europe. There is, however, inconsistency in the literature between the maximum salinity tolerance of the eggs in vivo, and the salinities at which some populations spawn in situ (7 vs. 9.6‰). In the present study, hatching success of perch was determined in vivo for a Danish, western Baltic, brackish water population at salinities of 4, 7, 10 and 12‰. Furthermore, in order to place the population genetically among other European perch populations, individual egg samples were sequenced for a 390 base pair fragment of the mtDNA D-loop region. Hatching occurred at all four salinities, with no statistical differences among treatments. Successful hatching at 12‰ is well above salinities of 7‰, which has previously been the highest reported from in vivo studies. This discrepancy is likely to be a result of methodological differences (e.g. different temperature) or perhaps interspecific variability in egg hatching abilities among perch populations. The perch from the present study consisted of a mixture of haplotypes similar to the haplotypes known to dominate in the Central and Western Europe and the Baltic Sea regions. Our results highlight the potential for recruitment of perch in coastal waters in the western Baltic Sea and have implications for both coastal management and aquaculture industry.