Hybridization between two high Arctic cetaceans confirmed by genomic analysis

Contrasting new and available reference genomes to highlight uncertainties in assemblies and areas for future improvement: an example with monodontid species

BACKGROUND

Reference genomes provide a foundational framework for evolutionary investigations, ecological analysis, and conservation science, yet uncertainties in the assembly of reference genomes are difficult to assess, and by extension rarely quantified. Reference genomes for monodontid cetaceans span a wide spectrum of data types and analytical approaches, providing the context to derive broader insights related to discrepancies and regions of uncertainty in reference genome assembly. We generated three beluga (Delphinapterus leucas) and one narwhal (Monodon monoceros) reference genomes and contrasted these with published chromosomal scale assemblies for each species to quantify discrepancies associated with genome assemblies.

RESULTS

The new reference genomes achieved chromosomal scale assembly using a combination of PacBio long reads, Illumina short reads, and Hi-C scaffolding data. For beluga, we identified discrepancies in the order and orientation of contigs in 2.2-3.7% of the total genome depending on the pairwise comparison of references. In addition, unsupported higher order scaffolding was identified in published reference genomes. In contrast, we estimated 8.2% of the compared narwhal genomes featured discrepancies, with inversions being notably abundant (5.3%). Discrepancies were linked to repetitive elements in both species.

CONCLUSIONS

We provide several new reference genomes for beluga (Delphinapterus leucas), while highlighting potential avenues for improvements. In particular, additional layers of data providing information on ultra-long genomic distances are needed to resolve persistent errors in reference genome construction. The comparative analyses of monodontid reference genomes suggested that the three new reference genomes for beluga are more accurate compared to the currently published reference genome, but that the new narwhal genome is less accurate than one published. We also present a conceptual summary for improving the accuracy of reference genomes with relevance to end-user needs and how they relate to levels of assembly quality and uncertainty.

Hybridization of Atlantic puffins in the Arctic coincides with 20th-century climate change

The Arctic is experiencing the fastest rates of global warming, leading to shifts in the distribution of its biota and increasing the potential for hybridization. However, genomic evidence of recent hybridization events in the Arctic remains unexpectedly rare. Here, we use whole-genome sequencing of contemporary and 122-year-old historical specimens to investigate the origin of an Arctic hybrid population of Atlantic puffins (Fratercula arctica) on Bjørnøya, Norway. We show that the hybridization between the High Arctic, large-bodied subspecies F. a. naumanni and the temperate, smaller-sized subspecies F. a. arctica began as recently as six generations ago due to an unexpected southward range expansion of F. a. naumanni. Moreover, we find a significant temporal loss of genetic diversity across Arctic and temperate puffin populations. Our observations provide compelling genomic evidence of the impacts of recent distributional shifts and loss of diversity in Arctic communities during the 20th century.

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.

Moving beyond the adaptationist paradigm for human evolution, and why it matters

The Journal of Human Evolution (JHE) was founded 50 years ago when much of the foundation for how we think about human evolution was in place or being put in place, providing the main framework for how we consider our origins today. Here, we will explore historical developments, including early JHE outputs, as they relate to our understanding of the relationship between phenotypic variation and evolutionary process, and use that as a springboard for considering our current understanding of these links as applied to human evolution. We will focus specifically on how the study of variation itself has shifted us away from taxonomic and adaptationist perspectives toward a richer understanding of the processes shaping human evolutionary history, using literature searches and specific test cases to highlight this. We argue that natural selection, gene exchange, genetic drift, and mutation should not be considered individually when considering the production of hominin diversity. In this context, we offer suggestions for future research directions and reflect on this more complex understanding of human evolution and its broader relevance to society. Finally, we end by considering authorship demographics and practices in the last 50 years within JHE and how a shift in these demographics has the potential to reshape the science of human evolution going forward.

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.

Mitochondrial diversity and inter-specific phylogeny among dolphins of the genus Stenella in the Southwest Atlantic Ocean

The genus Stenella is comprised of five species occurring in all oceans. Despite its wide distribution, genetic diversity information on these species is still scarce especially in the Southwest Atlantic Ocean. Some features of this genus can enhance opportunities for potential introgressive hybridization, e.g. sympatric distibution along the Brazilian coast, mixed known associations among species, karyotype uniformity and genome permeability. In this study we analyzed three genes of the mitochondrial genome to investigate the genetic diversity and occurrence of genetic mixture among eighty specimens of Stenella. All species exhibited moderate to high levels of genetic diversity (h = 0.833 to h = 1.000 and π = 0.006 to π = 0.015). Specimens of S. longirostris, S. attenuata and S. frontalis were clustered into differentiated haplogroups, in contrast, haplotypes of S. coeruleoalba and S. clymene were clustered together. We detected phylogenetic structure of mixed clades for S. clymene and S. coeruleoalba specimens, in the Southwest Atlantic Ocean, and also between S. frontalis and S. attenuata in the Northeast Atlantic Ocean, and between S. frontalis and S. longirostris in the Northwest Atlantic Ocean. These specimes were morphologically identified as one species but exhibited the maternal lineage of another species, by mitochondrial DNA. Our results demonstrate that ongoing gene flow is occurring among species of the genus Stenella reinforcing that this process could be one of the reasons for the confusing taxonomy and difficulties in elucidating phylogenetic relationships within this group.

Origin and expansion of the world's most widespread pinniped: range-wide population genomics of the harbour seal (Phoca vitulina)

The harbour seal (Phoca vitulina) is the most widely distributed pinniped, occupying a wide variety of habitats and climatic zones across the Northern Hemisphere. Intriguingly, the harbour seal is also one of the most philopatric seals, raising questions as to how it colonised virtually the whole of the Northern Hemisphere. To shed light on the origin, remarkable range expansion, population structure and genetic diversity of this species, we used genotyping-by-sequencing to analyse ~13,500 biallelic SNPs from 286 individuals sampled from 22 localities across the species' range. Our results point to a Northeast Pacific origin, colonisation of the North Atlantic via the Canadian Arctic, and subsequent stepping-stone range expansions across the North Atlantic from North America to Europe, accompanied by a successive loss of genetic diversity. Our analyses further revealed a deep divergence between modern North Pacific and North Atlantic harbour seals, with finer-scale genetic structure at regional and local scales consistent with strong philopatry. The study provides new insights into the harbour seal's remarkable ability to colonise and adapt to a wide range of habitats. Furthermore, it has implications for current harbour seal subspecies delineations and highlights the need for international and national red lists and management plans to ensure the protection of genetically and demographically isolated populations.

The influence of a lost society, the Sadlermiut, on the environment in the Canadian Arctic

High latitude freshwater ecosystems are sentinels of human activity and environmental change. The lakes and ponds that characterize Arctic landscapes have a low resilience to buffer variability in climate, especially with increasing global anthropogenic stressors in recent decades. Here, we show that a small freshwater pond in proximity of the archaeological site “Native Point” on Southampton Island (Nunavut, Arctic Canada) is a highly sensitive environmental recorder. The sediment analyses allowed for pinpointing the first arrival of Sadlermiut culture at Native Point to ~ 1250 CE, followed by a dietary shift likely in response to the onset of cooling in the region ~ 1400 CE. The influence of the Sadlermiut on the environment persisted long after the last of their population perished in 1903. Presently, the pond remains a distorted ecosystem that has experienced fundamental shifts in the benthic invertebrate assemblages and accumulated anthropogenic metals in the sediment. Our multi-proxy paleolimnological investigation using geochemical and biological indicators emphasizes that direct and indirect anthropogenic impacts have long-term environmental implications on high latitude ecosystems.

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.

Integrating fish scale and bone isotopic compositions for 'deep time' retrospective studies

Isotopic studies of archived fish scales have tremendous potential to develop long-term retrospectives that provide important insights into how humans have altered aquatic ecosystems. However, fish specimens in museum archives and other repositories typically date to time periods when the impacts of industrial societies may have already caused profound environmental changes. Archaeological fish bones offer an opportunity to bridge this key temporal gap by providing samples spanning from the recent past to as far back as the Pleistocene. Collagen is the primary protein component of both fish scale and bone, but the comparability of isotopic compositions from these tissues has not been established experimentally. To lay the framework for integrating isotopic datasets from these tissues, we compare δ¹³C and δ¹⁵N of bone and scale collagen, as well as other tissues, from fish with life-time controlled diets. Results show that while there is no difference in δ¹³C between scale and bone collagen, there may be a very slight but meaningful inter-tissue offset in δ¹⁵N (<0.3‰). We discuss potential sources of δ¹⁵N variation in scale and bone collagen measurements. Because there is no difference in scale and bone δ¹³C, and the observed offset in δ¹⁵N is very small (less than analytical uncertainty in many studies), our findings demonstrate that collagen isotopic compositions from these tissues should be directly comparable when integrating datasets from modern and ancient samples to build more powerful, millennium-scale isotopic times series. In linking isotopic compositions of collagen from modern, historical (scales), and archaeological (bones) fish, our findings open the way for more nuanced contextualization of how ecosystems functioned prior to large-scale exploitation and how they have responded to mounting anthropogenic pressures in the intervening centuries.

Species assignment in forensics and the challenge of hybrids

Forensic identification of species is in growing demand, particularly from law enforcement authorities in the areas of wildlife, fisheries and hunting as well as food authentication. Within the non-human forensic genetics expanding applications' field, the major current difficulties result from the lack of standards and genetic databases as well as the poor or absent taxonomic definition of several groups. Here we focus on a forensically important and overlooked problem in species identification: the exclusive use of uniparental markers, a common practice in current genetic barcoding methodologies, may lead to incorrect or impossible assignment whenever hybrids can occur (frequently, not only in domesticates, but also in the wild). For example, if one of these cases involves a mammal, and mitochondrial DNA alone is used (which in instances may be the only type of DNA sequence available in databases), the sample will be wrongfully assigned to the female parental species, completely missing the detection of a possible hybrid animal. The importance of this issue in the forensic contributions to food authentication, wildlife and conservation genetics is analyzed. We present a cautionary guidance on the forensic reporting of results avoiding this error.

A mitochondrial genetic divergence proxy predicts the reproductive compatibility of mammalian hybrids

Numerous pairs of evolutionarily divergent mammalian species have been shown to produce hybrid offspring. In some cases, F1 hybrids are able to produce F2s through matings with F1s. In other instances, the hybrids are only able to produce offspring themselves through backcrosses with a parent species owing to unisexual sterility (Haldane's Rule). Here, we explicitly tested whether genetic distance, computed from mitochondrial and nuclear genes, can be used as a proxy to predict the relative fertility of the hybrid offspring resulting from matings between species of terrestrial mammals. We assessed the proxy's predictive power using a well-characterized felid hybrid system, and applied it to modern and ancient hominins. Our results revealed a small overlap in mitochondrial genetic distance values that distinguish species pairs whose calculated distances fall within two categories: those whose hybrid offspring follow Haldane's Rule, and those whose hybrid F1 offspring can produce F2s. The strong correlation between genetic distance and hybrid fertility demonstrated here suggests that this proxy can be employed to predict whether the hybrid offspring of two mammalian species will follow Haldane's Rule.