Large-scale mitogenomic analysis of the phylogeography of the Late Pleistocene cave bear

Ursus spelaeus (Rosenmüller, 1794) during the MIS 3: new evidence from the Cioclovina Uscată Cave and radiocarbon age overview for the Carpathians

Ursus spelaeus, the Late Pleistocene a cave bear is known from numerous accumulations found in the fossil sector of caves situated in the Carpathian and Apuseni Mountains. In this study, we present new radiocarbon data along a profile of the Cioclovina Uscată Cave, which is situated in the South Carpathians. The data suggest that, during the entire Marine Isotope Stage 3 (MIS 3) interval, the cave was serving as a shelter for U. spelaeus, with the oldest dated bone indicating an age of > 47,710 and the youngest one, an age of 31,820 ± 400 years cal BP. Histogram plots of 110 radiocarbon data from different caves of the Carpathian and Apuseni Mountains as Cioclovina Uscată, Peștera (Cave) cu Oase, Peștera Muierii, or Peștera Urșilor, respectively, show a maximum expansion of the cave bear population between 50,000 and 40,000, a decline between 40,000 and 35,000 and a partial recovery from 35,000-30,000 years cal BP. Radiocarbon data of Homo sapiens remains, younger than 35,000 years cal BP, support the fact that H. sapiens accessed the same caves where the cave bear persisted to hibernate. Besides general cool conditions and restricted food sources, the presence of H. sapiens constituted an additional stress factor driving the cave bear to extinction.

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

Sequencing and assembling bear genomes: the bare necessities

Unique genetic adaptations are present in bears of every species across the world. From (nearly) shutting down important organs during hibernation to preventing harm from lifestyles that could easily cause metabolic diseases in humans, bears may hold the answer to various human ailments. However, only a few of these unique traits are currently being investigated at the molecular level, partly because of the lack of necessary tools. One of these tools is well-annotated genome assemblies from the different, extant bear species. These reference genomes are needed to allow us to identify differences in genetic variants, isoforms, gene expression, and genomic features such as transposons and identify those that are associated with biomedical-relevant traits. In this review we assess the current state of the genome assemblies of the eight different bear species, discuss current gaps, and the future benefits these reference genomes may have in informing human biomedical applications, while at the same time improving bear conservation efforts.

Ancient mitochondrial and modern whole genomes unravel massive genetic diversity loss during near extinction of Alpine ibex

Population bottlenecks can have dramatic consequences for the health and long-term survival of a species. Understanding of historic population size and standing genetic variation prior to a contraction allows estimating the impact of a bottleneck on the species genetic diversity. Although historic population sizes can be modelled based on extant genomics, uncertainty is high for the last 10-20 millenia. Hence, integrating ancient genomes provides a powerful complement to retrace the evolution of genetic diversity through population fluctuations. Here, we recover 15 high-quality mitogenomes of the once nearly extinct Alpine ibex spanning 8601 BP to 1919 CE and combine these with 60 published modern whole genomes. Coalescent demography simulations based on modern whole genomes indicate population fluctuations coinciding with the last major glaciation period. Using our ancient and historic mitogenomes, we investigate the more recent demographic history of the species and show that mitochondrial haplotype diversity was reduced to a fifth of the pre-bottleneck diversity with several highly differentiated mitochondrial lineages having co-existed historically. The main collapse of mitochondrial diversity coincides with elevated human population growth during the last 1-2 kya. After recovery, one lineage was spread and nearly fixed across the Alps due to recolonization efforts. Our study highlights that a combined approach integrating genomic data of ancient, historic and extant populations unravels major long-term population fluctuations from the emergence of a species through its near extinction up to the recent past.

Integrating multi-taxon palaeogenomes and sedimentary ancient DNA to study past ecosystem dynamics

Ancient DNA (aDNA) has played a major role in our understanding of the past. Important advances in the sequencing and analysis of aDNA from a range of organisms have enabled a detailed understanding of processes such as past demography, introgression, domestication, adaptation and speciation. However, to date and with the notable exception of microbiomes and sediments, most aDNA studies have focused on single taxa or taxonomic groups, making the study of changes at the community level challenging. This is rather surprising because current sequencing and analytical approaches allow us to obtain and analyse aDNA from multiple source materials. When combined, these data can enable the simultaneous study of multiple taxa through space and time, and could thus provide a more comprehensive understanding of ecosystem-wide changes. It is therefore timely to develop an integrative approach to aDNA studies by combining data from multiple taxa and substrates. In this review, we discuss the various applications, associated challenges and future prospects of such an approach.

Ancient Faunal History Revealed by Interdisciplinary Biomolecular Approaches

Starting four decades ago, studies have examined the ecology and evolutionary dynamics of populations and species using short mitochondrial DNA fragments and stable isotopes. Through technological and analytical advances, the methods and biomolecules at our disposal have increased significantly to now include lipids, whole genomes, proteomes, and even epigenomes. At an unprecedented resolution, the study of ancient biomolecules has made it possible for us to disentangle the complex processes that shaped the ancient faunal diversity across millennia, with the potential to aid in implicating probable causes of species extinction and how humans impacted the genetics and ecology of wild and domestic species. However, even now, few studies explore interdisciplinary biomolecular approaches to reveal ancient faunal diversity dynamics in relation to environmental and anthropogenic impact. This review will approach how biomolecules have been implemented in a broad variety of topics and species, from the extinct Pleistocene megafauna to ancient wild and domestic stocks, as well as how their future use has the potential to offer an enhanced understanding of drivers of past faunal diversity on Earth.

Early Alpine occupation backdates westward human migration in Late Glacial Europe

Before the end of the Last Glacial Maximum (LGM, ∼16.5 ka ago)¹ set in motion major shifts in human culture and population structure,² a consistent change in lithic technology, material culture, settlement pattern, and adaptive strategies is recorded in Southern Europe at ∼18-17 ka ago. In this time frame, the landscape of Northeastern Italy changed considerably, and the retreat of glaciers allowed hunter-gatherers to gradually recolonize the Alps.^3-6 Change within this renewed cultural frame (i.e., during the Late Epigravettian phase) is currently associated with migrations favored by warmer climate linked to the Bølling-Allerød onset (14.7 ka ago),^7-11 which replaced earlier genetic lineages with ancestry found in an individual who lived ∼14 ka ago at Riparo Villabruna, Italy, and shared among different contexts (Villabruna Cluster).⁹ Nevertheless, these dynamics and their chronology are still far from being disentangled due to fragmentary evidence for long-distance interactions across Europe.¹² Here, we generate new genomic data from a human mandible uncovered at Riparo Tagliente (Veneto, Italy), which we directly dated to 16,980-16,510 cal BP (2σ). This individual, affected by focal osseous dysplasia, is genetically affine to the Villabruna Cluster. Our results therefore backdate by at least 3 ka the diffusion in Southern Europe of a genetic component linked to Balkan/Anatolian refugia, previously believed to have spread during the later Bølling/Allerød event. In light of the new genetic evidence, this population replacement chronologically coincides with the very emergence of major cultural transitions in Southern and Western Europe.

Exploring the phylogeography and population dynamics of the giant deer (Megaloceros giganteus) using Late Quaternary mitogenomes

Late Quaternary climatic fluctuations in the Northern Hemisphere had drastic effects on large mammal species, leading to the extinction of a substantial number of them. The giant deer (Megaloceros giganteus) was one of the species that became extinct in the Holocene, around 7660 calendar years before present. In the Late Pleistocene, the species ranged from western Europe to central Asia. However, during the Holocene, its range contracted to eastern Europe and western Siberia, where the last populations of the species occurred. Here, we generated 35 Late Pleistocene and Holocene giant deer mitogenomes to explore the genetics of the demise of this iconic species. Bayesian phylogenetic analyses of the mitogenomes suggested five main clades for the species: three pre-Last Glacial Maximum clades that did not appear in the post-Last Glacial Maximum genetic pool, and two clades that showed continuity into the Holocene. Our study also identified a decrease in genetic diversity starting in Marine Isotope Stage 3 and accelerating during the Last Glacial Maximum. This reduction in genetic diversity during the Last Glacial Maximum, coupled with a major contraction of fossil occurrences, suggests that climate was a major driver in the dynamics of the giant deer.

Middle Pleistocene genome calibrates a revised evolutionary history of extinct cave bears

Palaeogenomes provide the potential to study evolutionary processes in real time, but this potential is limited by our ability to recover genetic data over extended timescales.¹ As a consequence, most studies so far have focused on samples of Late Pleistocene or Holocene age, which covers only a small part of the history of many clades and species. Here, we report the recovery of a low coverage palaeogenome from the petrous bone of a ∼360,000 year old cave bear from Kudaro 1 cave in the Caucasus Mountains. Analysis of this genome alongside those of several Late Pleistocene cave bears reveals widespread mito-nuclear discordance in this group. Using the time interval between Middle and Late Pleistocene cave bear genomes, we directly estimate ursid nuclear and mitochondrial substitution rates to calibrate their respective phylogenies. This reveals post-divergence mitochondrial transfer as the dominant factor explaining their mito-nuclear discordance. Interestingly, these transfer events were not accompanied by large-scale nuclear introgression. However, we do detect additional instances of nuclear admixture among other cave bear lineages, and between cave bears and brown bears, which are not associated with mitochondrial exchange. Genomic data obtained from the Middle Pleistocene cave bear petrous bone has thus facilitated a revised evolutionary history of this extinct megafaunal group. Moreover, it suggests that petrous bones may provide a means of extending both the magnitude and time depth of palaeogenome retrieval over substantial portions of the evolutionary histories of many mammalian clades.

Cost-effective straightforward method for captured whole mitogenome sequencing of ancient DNA

Working with mitochondrial DNA from highly degraded samples is challenging. We present a whole mitogenome Illumina-based sequencing method suitable for highly degraded samples. The method makes use of double-stranded library preparation with hybridization-based target enrichment. The aim of the study was to implement a new user-friendly method for analysing many ancient DNA samples at low cost. The method combines the Swift 2S™ Turbo library preparation kit and xGen® panel for mitogenome enrichment. Swift allows to use low input of aDNA and own adapters and primers, handles inhibitors well, and has only two purification steps. xGen is straightforward to use and is able to leverage already pooled libraries. Given the ancient DNA is more challenging to work with, the protocol was developed with several improvements, especially multiplying DNA input in case of low concentration DNA extractions followed by AMPure® beads size selection and real-time pre-capture PCR monitoring in order to avoid cycle-optimization step. Nine out of eleven analysed samples successfully retrieved mitogenomes. Hence, our method provides an effective analysis of whole mtDNA, and has proven to be fast, cost-effective, straightforward, with utilisation in population-wide research of burial sites.

Genomic patterns in the widespread Eurasian lynx shaped by Late Quaternary climatic fluctuations and anthropogenic impacts

Disentangling the contribution of long-term evolutionary processes and recent anthropogenic impacts to current genetic patterns of wildlife species is key to assessing genetic risks and designing conservation strategies. Here, we used 80 whole nuclear genomes and 96 mitogenomes from populations of the Eurasian lynx covering a range of conservation statuses, climatic zones and subspecies across Eurasia to infer the demographic history, reconstruct genetic patterns, and discuss the influence of long-term isolation and/or more recent human-driven changes. Our results show that Eurasian lynx populations shared a common history until 100,000 years ago, when Asian and European populations started to diverge and both entered a period of continuous and widespread decline, with western populations, except Kirov, maintaining lower effective sizes than eastern populations. Population declines and increased isolation in more recent times probably drove the genetic differentiation between geographically and ecologically close westernmost European populations. By contrast, and despite the wide range of habitats covered, populations are quite homogeneous genetically across the Asian range, showing a pattern of isolation by distance and providing little genetic support for the several proposed subspecies. Mitogenomic and nuclear divergences and population declines starting during the Late Pleistocene can be mostly attributed to climatic fluctuations and early human influence, but the widespread and sustained decline since the Holocene is more probably the consequence of anthropogenic impacts which intensified in recent centuries, especially in western Europe. Genetic erosion in isolated European populations and lack of evidence for long-term isolation argue for the restoration of lost population connectivity.