Species-specific responses of Late Quaternary megafauna to climate and humans

52,000 years of woolly rhinoceros population dynamics reveal extinction mechanisms

The extinction of the woolly rhinoceros (Coelodonta antiquitatis) at the onset of the Holocene remains an enigma, with conflicting evidence regarding its cause and spatiotemporal dynamics. This partly reflects challenges in determining demographic responses of late Quaternary megafauna to climatic and anthropogenic causal drivers with available genetic and paleontological techniques. Here, we show that elucidating mechanisms of ancient extinctions can benefit from a detailed understanding of fine-scale metapopulation dynamics, operating over many millennia. Using an abundant fossil record, ancient DNA, and high-resolution simulation models, we untangle the ecological mechanisms and causal drivers that are likely to have been integral in the decline and later extinction of the woolly rhinoceros. Our 52,000-y reconstruction of distribution-wide metapopulation dynamics supports a pathway to extinction that began long before the Holocene, when the combination of cooling temperatures and low but sustained hunting by humans trapped woolly rhinoceroses in suboptimal habitats along the southern edge of their range. Modeling indicates that this ecological trap intensified after the end of the last ice age, preventing colonization of newly formed suitable habitats, weakening stabilizing metapopulation processes, triggering the extinction of the woolly rhinoceros in the early Holocene. Our findings suggest that fragmentation and resultant metapopulation dynamics should be explicitly considered in explanations of late Quaternary megafauna extinctions, sending a clarion call to the fragility of the remaining large-bodied grazers restricted to disjunct fragments of poor-quality habitat due to anthropogenic environmental change.

Interactive effects of drought and deforestation on multitrophic communities and aquatic ecosystem functions in the Neotropics-a test using tank bromeliads

Background

Together with the intensification of dry seasons in Neotropical regions, increasing deforestation is expected to exacerbate species extinctions, something that could lead to dramatic shifts in multitrophic communities and ecosystem functions. Recent studies suggest that the effects of habitat loss are greater where precipitation has decreased. Yet, experimental studies of the pure and interactive effects of drought and deforestation at ecosystem level remain scarce.

Methods

Here, we used rainshelters and transplantation from rainforest to open areas of natural microcosms (the aquatic ecosystem and microbial-faunal food web found within the rainwater-filled leaves of tank bromeliads) to emulate drought and deforestation in a full factorial experimental design. We analysed the pure and interactive effects of our treatments on functional community structure (including microorganisms, detritivore and predatory invertebrates), and on leaf litter decomposition in tank bromeliad ecosystems.

Results

Drought or deforestation alone had a moderate impact on biomass at the various trophic level, but did not eliminate species. However, their interaction synergistically reduced the biomass of all invertebrate functional groups and bacteria. Predators were the most impacted trophic group as they were totally eliminated, while detritivore biomass was reduced by about 95%. Fungal biomass was either unaffected or boosted by our treatments. Decomposition was essentially driven by microbial activity, and did not change across treatments involving deforestation and/or drought.

Conclusions

Our results suggest that highly resistant microorganisms such as fungi (plus a few detritivores) maintain key ecosystem functions in the face of drought and habitat change. We conclude that habitat destruction compounds the problems of climate change, that the impacts of the two phenomena on food webs are mutually reinforcing, and that the stability of ecosystem functions depends on the resistance of a core group of organisms. Assuming that taking global action is more challenging than taking local-regional actions, policy-makers should be encouraged to implement environmental action plans that will halt habitat destruction, to dampen any detrimental interactive effect with the impacts of global climate change.

A palaeogenomic investigation of overharvest implications in an endemic wild reindeer subspecies

Overharvest can severely reduce the abundance and distribution of a species and thereby impact its genetic diversity and threaten its future viability. Overharvest remains an ongoing issue for Arctic mammals, which due to climate change now also confront one of the fastest changing environments on Earth. The high-arctic Svalbard reindeer (Rangifer tarandus platyrhynchus), endemic to Svalbard, experienced a harvest-induced demographic bottleneck that occurred during the 17-20th centuries. Here, we investigate changes in genetic diversity, population structure, and gene-specific differentiation during and after this overharvesting event. Using whole-genome shotgun sequencing, we generated the first ancient and historical nuclear (n = 11) and mitochondrial (n = 18) genomes from Svalbard reindeer (up to 4000 BP) and integrated these data with a large collection of modern genome sequences (n = 90) to infer temporal changes. We show that hunting resulted in major genetic changes and restructuring in reindeer populations. Near-extirpation followed by pronounced genetic drift has altered the allele frequencies of important genes contributing to diverse biological functions. Median heterozygosity was reduced by 26%, while the mitochondrial genetic diversity was reduced only to a limited extent, likely due to already low pre-harvest diversity and a complex post-harvest recolonization process. Such genomic erosion and genetic isolation of populations due to past anthropogenic disturbance will likely play a major role in metapopulation dynamics (i.e., extirpation, recolonization) under further climate change. Our results from a high-arctic case study therefore emphasize the need to understand the long-term interplay of past, current, and future stressors in wildlife conservation.

Genomic Analyses Capture the Human-Induced Demographic Collapse and Recovery in a Wide-Ranging Cervid

The glacial cycles of the Quaternary heavily impacted species through successions of population contractions and expansions. Similarly, populations have been intensely shaped by human pressures such as unregulated hunting and land use changes. White-tailed and mule deer survived in different refugia through the Last Glacial Maximum, and their populations were severely reduced after the European colonization. Here, we analyzed 73 resequenced deer genomes from across their North American range to understand the consequences of climatic and anthropogenic pressures on deer demographic and adaptive history. We found strong signals of climate-induced vicariance and demographic decline; notably, multiple sequentially Markovian coalescent recovers a severe decline in mainland white-tailed deer effective population size (Ne) at the end of the Last Glacial Maximum. We found robust evidence for colonial overharvest in the form of a recent and dramatic drop in Ne in all analyzed populations. Historical census size and restocking data show a clear parallel to historical Ne estimates, and temporal Ne/Nc ratio shows patterns of conservation concern for mule deer. Signatures of selection highlight genes related to temperature, including a cold receptor previously highlighted in woolly mammoth. We also detected immune genes that we surmise reflect the changing land use patterns in North America. Our study provides a detailed picture of anthropogenic and climatic-induced decline in deer diversity and clues to understanding the conservation concerns of mule deer and the successful demographic recovery of white-tailed deer.

Population genomics of the muskox' resilience in the near absence of genetic variation

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

Ancient mitogenomes reveal a high maternal genetic diversity of Pleistocene woolly rhinoceros in Northern China

BACKGROUND

Woolly rhinoceros (Coelodonta antiquitatis) is a typical indicator of cold-stage climate that was widely distributed in Northern Hemisphere during the Middle-Late Pleistocene. Although a plethora of fossils have been excavated from Northern China, their phylogenetic status, intraspecific diversity and phylogeographical structure are still vague.

RESULTS

In the present study, we generated four mitogenomes from Late Pleistocene woolly rhinoceros in Northern China and compared them with published data. Bayesian and network analyses indicate that the analyzed individuals contain at least four maternal haplogroups, and Chinese samples fall in three of them. One of our samples belongs to a previously unidentified early diverging clade (haplogroup D), which separated from other woolly rhinoceros around 0.57 Ma (95% CI: 0.76-0.41 Ma). The timing of this clade's origin coincides with the first occurrence of woolly rhinoceros, which are thought to have evolved in Europe. Our other three samples cluster in haplogroup C, previously only identified from one specimen from Wrangel Island (ND030) and initially considered to be an isolated clade. Herein, our findings suggest that ND030 is likely descended from a northward dispersal of the individuals carrying haplogroup C from Northern China. Additionally, Chinese woolly rhinoceros specimens exhibit higher nucleotide diversity than those from Siberia.

CONCLUSION

Our findings highlight Northern China as a possible refugium and a key evolution center of the Pleistocene woolly rhinoceros.

Was the Late Ordovician mass extinction truly exceptional?

The Late Ordovician mass extinction event is the oldest of the five great extinction events in the fossil record. It has long been regarded as an outlier among mass extinctions, primarily due to its association with a cooling climate. However, recent temporally better resolved fossil biodiversity estimates complicate this view, providing growing evidence for a prolonged but punctuated biodiversity decline modulated by changes in atmospheric composition, ocean chemistry, and viable habitat area. This evolving view invokes extinction drivers similar to those that occurred during other major extinctions; some are even factors in the current human-induced biodiversity crisis. Even this very ancient and, at first glance, exceptional event conveys important lessons about the intensifying 'sixth mass extinction'.

Whole-genome sequencing reveals adaptations of hairy-footed jerboas (Dipus, Dipodidae) to diverse desert environments

BACKGROUND

Environmental conditions vary among deserts across the world, spanning from hyper-arid to high-elevation deserts. However, prior genomic studies on desert adaptation have focused on desert and non-desert comparisons overlooking the complexity of conditions within deserts. Focusing on the adaptation mechanisms to diverse desert environments will advance our understanding of how species adapt to extreme desert environments. The hairy-footed jerboas are well adapted to diverse desert environments, inhabiting high-altitude arid regions, hyper-arid deserts, and semi-deserts, but the genetic basis of their adaptation to different deserts remains unknown.

RESULTS

Here, we sequenced the whole genome of 83 hairy-footed jerboas from distinct desert zones in China to assess how they responded under contrasting conditions. Population genomics analyses reveal the existence of three species in hairy-footed jerboas distributed in China: Dipus deasyi, Dipus sagitta, and Dipus sowerbyi. Analyses of selection between high-altitude desert (elevation ≥ 3000m) and low-altitude desert (< 500m) populations identified two strongly selected genes, ATR and HIF1AN, associated with intense UV radiation and hypoxia in high-altitude environments. A number of candidate genes involved in energy and water homeostasis were detected in the comparative genomic analyses of hyper-arid desert (average annual precipitation < 70mm) and arid desert (< 200mm) populations versus semi-desert (> 360mm) populations. Hyper-arid desert animals also exhibited stronger adaptive selection in energy homeostasis, suggesting water and resource scarcity may be the main drivers of desert adaptation in hairy-footed jerboas.

CONCLUSIONS

Our study challenges the view of deserts as homogeneous environments and shows that distinct genomic adaptations can be found among desert animals depending on their habitats.

Evolutionary history of the extinct wolf population from France in the context of global phylogeographic changes throughout the Holocene

Phylogeographic patterns in large mammals result from natural environmental factors and anthropogenic effects, which in some cases include domestication. The grey wolf was once widely distributed across the Holarctic, but experienced phylogeographic shifts and demographic declines during the Holocene. In the 19th-20th centuries, the species became extirpated from large parts of Europe due to direct extermination and habitat loss. We reconstructed the evolutionary history of the extinct Western European wolves based on the mitogenomic composition of 78 samples from France (Neolithic-20th century) in the context of other populations of wolves and dogs worldwide. We found a close genetic similarity of French wolves from ancient, medieval and recent populations, which suggests the long-term continuity of maternal lineages. MtDNA haplotypes of the French wolves showed large diversity and fell into two main haplogroups of modern Holarctic wolves. Our worldwide phylogeographic analysis indicated that haplogroup W1, which includes wolves from Eurasia and North America, originated in Northern Siberia. Haplogroup W2, which includes only European wolves, originated in Europe ~35 kya and its frequency was reduced during the Holocene due to an expansion of haplogroup W1 from the east. Moreover, we found that dog haplogroup D, currently restricted to Europe and the Middle East, was nested within the wolf haplogroup W2. This suggests European origin of haplogroup D, probably as a result of an ancient introgression from European wolves. Our results highlight the dynamic evolutionary history of European wolves during the Holocene, with a partial lineage replacement and introgressive hybridization with local dog populations.

Mitogenomes revealed the history of bison colonization of Northern Plains after the Last Glacial Maximum

American bison demonstrated differential patterns of extinction, survival, and expansion since the terminal Pleistocene. We determined population dynamics of the Northern Great Plains bison using 40 mitochondrial genomes from radiocarbon dated remains with the age ranging from 12,226 to 167 calibrated years before present. Population dynamics correlated with environmental and anthropogenic factors and was characterized by three primary periods: terminal Pleistocene population growth starting 14,000 years ago, mid Holocene demographic stability between 6700 and 2700 years ago, and late Holocene population decline in the last 2700 years. Most diversification of mtDNA haplotypes occurred in the early Holocene when bison colonized new territories opened by retreating ice sheets. Holocene mtDNA lineages were not found in modern bison and lacked association with archaeological sites and morphological forms.

Ecospaces of the Middle to Upper Paleolithic transition: The archaeofaunal record of the Iberian Peninsula

The rich archaeofaunal record of Iberia provides a means of exploring potential differences between Neanderthal and anatomically modern human interactions with the environment. In this article, we present an analysis of Iberian archaeofaunas dating between 60 and 30 ka to explore if, how, and why the faunal ecospaces of Neanderthals and anatomically modern humans differed. We test for impacts of chronology (as a proxy for Neanderthal and anatomically modern human exploitation) and environmental regionalization (using bioclimatic regions) on archaeofaunal composition, using a combination of cluster (unweighted pair-group method using arithmetic averages) and nonmetric multidimensional scaling. Our chronological analysis finds no significant compositional difference between Neanderthal and anatomically modern mammalian faunal assemblages; however, bioclimatic regionalization is stronger in anatomically modern human-affiliated assemblages than in Neanderthal ones, a finding that may indicate a difference in site occupation duration or foraging mobility between Neanderthals and anatomically modern humans.

Climate-driven habitat shifts of high-ranked prey species structure Late Upper Paleolithic hunting

Changing climates in the past affected both human and faunal population distributions, thereby structuring human diets, demography, and cultural evolution. Yet, separating the effects of climate-driven and human-induced changes in prey species abundances remains challenging, particularly during the Late Upper Paleolithic, a period marked by rapid climate change and marked ecosystem transformation. To disentangle the effects of climate and hunter-gatherer populations on animal prey species during the period, we synthesize disparate paleoclimate records, zooarchaeological data, and archaeological data using ecological methods and theory to test to what extent climate and anthropogenic impacts drove broad changes in human subsistence observed in the Late Upper Paleolithic zooarchaeological records. We find that the observed changes in faunal assemblages during the European Late Upper Paleolithic are consistent with climate-driven animal habitat shifts impacting the natural abundances of high-ranked prey species on the landscape rather than human-induced resource depression. The study has important implications for understanding how past climate change impacted and structured the diet and demography of human populations and can serve as a baseline for considerations of resilience and adaptation in the present.

Population Structure, Demographic History, and Adaptation of Giant Honeybees in China Revealed by Population Genomic Data

There have been many population-based genomic studies on human-managed honeybees (Apis mellifera and Apis cerana), but there has been a notable lack of analysis with regard to wild honeybees, particularly in relation to their evolutionary history. Nevertheless, giant honeybees have been found to occupy distinct habitats and display remarkable characteristics, which are attracting an increased amount of attention. In this study, we de novo sequenced and then assembled the draft genome sequence of the Himalayan giant honeybee, Apis laboriosa. Phylogenetic analysis based on genomic information indicated that A. laboriosa and its tropical sister species Apis dorsata diverged ∼2.61 Ma, which supports the speciation hypothesis that links A. laboriosa to geological changes throughout history. Furthermore, we re-sequenced A. laboriosa and A. dorsata samples from five and six regions, respectively, across their population ranges in China. These analyses highlighted major genetic differences for Tibetan A. laboriosa as well as the Hainan Island A. dorsata. The demographic history of most giant honeybee populations has mirrored glacial cycles. More importantly, contrary to what has occurred among human-managed honeybees, the demographic history of these two wild honeybee species indicates a rapid decline in effective population size in the recent past, reflecting their differences in evolutionary histories. Several genes were found to be subject to selection, which may help giant honeybees to adapt to specific local conditions. In summary, our study sheds light on the evolutionary and adaptational characteristics of two wild giant honeybee species, which was useful for giant honeybee conservation.

Ancient DNA of narrow-headed vole reveal common features of the Late Pleistocene population dynamics in cold-adapted small mammals

The narrow-headed vole, collared lemming and common vole were the most abundant small mammal species across the Eurasian Late Pleistocene steppe-tundra environment. Previous ancient DNA studies of the collared lemming and common vole have revealed dynamic population histories shaped by climatic fluctuations. To investigate the extent to which species with similar adaptations share common evolutionary histories, we generated a dataset comprised the mitochondrial genomes of 139 ancient and 6 modern narrow-headed voles from several sites across Europe and northwestern Asia covering approximately the last 100 thousand years (kyr). We inferred Bayesian time-aware phylogenies using 11 radiocarbon-dated samples to calibrate the molecular clock. Divergence of the main mtDNA lineages across the three species occurred during marine isotope stages (MIS) 7 and MIS 5, suggesting a common response of species adapted to open habitat during interglacials. We identified several time-structured mtDNA lineages in European narrow-headed vole, suggesting lineage turnover. The timing of some of these turnovers was synchronous across the three species, allowing us to identify the main drivers of the Late Pleistocene dynamics of steppe- and cold-adapted species.

The Exceptional Presence of Megaloceros giganteus in North-Eastern Iberia and Its Palaeoecological Implications: The Case of Teixoneres Cave (Moià, Barcelona, Spain)

In this article we announce the discovery of the first remains of Megaloceros giganteus found in Catalonia (north-eastern Iberia) from the Late Pleistocene: a fragment of maxillary. Dated between 35,000 and 37,000 cal BP, it is also among the youngest occurrence of this taxon in the Iberian Peninsula, while its last known occurrence is dated to the Neolithic period. Through a comparison with the giant deer of the northern Pyrenees, we analyzed the herbivore guilds in which this taxon was associated to understand the context in which it was able to enter the Iberian Peninsula. By comparing its diet with those of specimens from Northern Europe, we detail the ecological adaptations of this taxon in this new environment. We suggest that Megaloceros accompanied the migrations of cold-adapted species by taking advantage of the opening of corridors on both sides of the Pyrenees during the coldest periods of the Late Pleistocene. The diet of the Iberian individuals, which is oriented towards abrasive plants, suggests an adaptation to a different ecological niche than that found in Northern European individuals. The northern Iberian Peninsula may have been an extreme in the geographical expansion of M. giganteus. More specimens will be needed in the future to establish the variability of the southern Megaloceros populations.

The cumulative impacts of anthropogenic stressors vary markedly along environmental gradients

Understanding the cumulative effects of multiple stressors on biodiversity is key to managing their impacts. Stressor interactions are often studied using an additive/antagonistic/synergistic typology, aimed at identifying situations where individual stressor effects are reduced or amplified when they act in combination. Here, we analysed variation in the family richness of stream macroinvertebrates in the groups Ephemeroptera, Plecoptera and Trichoptera (EPT) at 4658 sites spanning a 32° latitudinal range in eastern Australia in relation to two largely human-induced stressors, salinity and turbidity, and two environmental gradients, temperature and slope. The cumulative and interactive effect of salinity and turbidity on EPT family richness varied across the landscape and by habitat (edge or riffle) such that we observed additive, antagonistic and synergistic outcomes depending on the environmental context. Our findings highlight the importance of understanding the consistency of multiple stressor impacts, which will involve higher-order interactions between multiple stressors and environmental factors.

A long-term perspective on Neanderthal environment and subsistence: Insights from the dental microwear texture analysis of hunted ungulates at Combe-Grenal (Dordogne, France)

Large bovids and cervids constituted major components of the European Middle Palaeolithic faunas and hence a key resource for Neanderthal populations. In paleoenvironmental reconstructions, red deer (Cervus elaphus) occurrence is classically considered as a tree-cover indicator while Bovinae (Bison priscus and Bos primigenius) and reindeer (Rangifer tarandus) occurrences are typically associated with open landscapes. However, insights into the ecology of extant ungulate populations show a more complex reality. Exploring the diet of past ungulates allows to better comprehend the hunting strategies of Palaeolithic populations and to reconstruct the modifications through time of past landscapes. By reflecting what animals have eaten during the last days or weeks of their life, dental microwear textures of herbivores link a population and its environment. Here we analyzed, via Dental Microwear Texture Analysis (DMTA), the diet of 50 Bos/Bison, 202 R. tarandus and 116 C. elaphus preyed upon by the Neanderthals that occupied Combe-Grenal rock-shelter, one of the most important Mousterian archaeo-sequences in southwestern France considering its long stratigraphy, abundance of faunal remains and the variations perceptible in Palaeolithic material culture. Grazers and mixed-feeders are the most represented dietary categories among Combe-Grenal's guild of herbivores, highlighting the availability, along the sequence, of open landscapes. The absence of clear changes in the use of plant resources by hunted ungulates through time, even though palaeoenvironmental changes were well-documented by previous studies along the sequence, is interpreted as resulting from the hunting of non-randomly selected prey by Neanderthals, preferentially in open environments. Thus, these results provide further insight into the hunting strategies of Neanderthals and modify our perception of potential links between subsistence and material culture. Combe-Grenal hunters "stayed in the open" through millennia, and were not forced to switch to hunting tactics and material technology adapted to close encounters in forested environments.

Lions and brown bears colonized North America in multiple synchronous waves of dispersal across the Bering Land Bridge

The Bering Land Bridge connecting North America and Eurasia was periodically exposed and inundated by oscillating sea levels during the Pleistocene glacial cycles. This land connection allowed the intermittent dispersal of animals, including humans, between Western Beringia (far northeast Asia) and Eastern Beringia (northwest North America), changing the faunal community composition of both continents. The Pleistocene glacial cycles also had profound impacts on temperature, precipitation and vegetation, impacting faunal community structure and demography. While these palaeoenvironmental impacts have been studied in many large herbivores from Beringia (e.g., bison, mammoths, horses), the Pleistocene population dynamics of the diverse guild of carnivorans present in the region are less well understood, due to their lower abundances. In this study, we analyse mitochondrial genome data from ancient brown bears (Ursus arctos; n = 103) and lions (Panthera spp.; n = 39), two megafaunal carnivorans that dispersed into North America during the Pleistocene. Our results reveal striking synchronicity in the population dynamics of Beringian lions and brown bears, with multiple waves of dispersal across the Bering Land Bridge coinciding with glacial periods of low sea levels, as well as synchronous local extinctions in Eastern Beringia during Marine Isotope Stage 3. The evolutionary histories of these two taxa underline the crucial biogeographical role of the Bering Land Bridge in the distribution, turnover and maintenance of megafaunal populations in North America.

Population dynamics and demographic history of Eurasian collared lemmings

BACKGROUND

Ancient DNA studies suggest that Late Pleistocene climatic changes had a significant effect on population dynamics in Arctic species. The Eurasian collared lemming (Dicrostonyx torquatus) is a keystone species in the Arctic ecosystem. Earlier studies have indicated that past climatic fluctuations were important drivers of past population dynamics in this species.

RESULTS

Here, we analysed 59 ancient and 54 modern mitogenomes from across Eurasia, along with one modern nuclear genome. Our results suggest population growth and genetic diversification during the early Late Pleistocene, implying that collared lemmings may have experienced a genetic bottleneck during the warm Eemian interglacial. Furthermore, we find multiple temporally structured mitogenome clades during the Late Pleistocene, consistent with earlier results suggesting a dynamic late glacial population history. Finally, we identify a population in northeastern Siberia that maintained genetic diversity and a constant population size at the end of the Pleistocene, suggesting suitable conditions for collared lemmings in this region during the increasing temperatures associated with the onset of the Holocene.

CONCLUSIONS

This study highlights an influence of past warming, in particular the Eemian interglacial, on the evolutionary history of the collared lemming, along with spatiotemporal population structuring throughout the Late Pleistocene.

Spatiotemporal influences of climate and humans on muskox range dynamics over multiple millennia

Processes leading to range contractions and population declines of Arctic megafauna during the late Pleistocene and early Holocene are uncertain, with intense debate on the roles of human hunting, climatic change, and their synergy. Obstacles to a resolution have included an overreliance on correlative rather than process-explicit approaches for inferring drivers of distributional and demographic change. Here, we disentangle the ecological mechanisms and threats that were integral in the decline and extinction of the muskox (Ovibos moschatus) in Eurasia and in its expansion in North America using process-explicit macroecological models. The approach integrates modern and fossil occurrence records, ancient DNA, spatiotemporal reconstructions of past climatic change, species-specific population ecology, and the growth and spread of anatomically modern humans. We show that accurately reconstructing inferences of past demographic changes for muskox over the last 21,000 years require high dispersal abilities, large maximum densities, and a small Allee effect. Analyses of validated process-explicit projections indicate that climatic change was the primary driver of muskox distribution shifts and demographic changes across its previously extensive (circumpolar) range, with populations responding negatively to rapid warming events. Regional analyses show that the range collapse and extinction of the muskox in Europe (~13,000 years ago) was likely caused by humans operating in synergy with climatic warming. In Canada and Greenland, climatic change and human activities probably combined to drive recent population sizes. The impact of past climatic change on the range and extinction dynamics of muskox during the Pleistocene-Holocene transition signals a vulnerability of this species to future increased warming. By better establishing the ecological processes that shaped the distribution of the muskox through space and time, we show that process-explicit macroecological models have important applications for the future conservation and management of this iconic species in a warming Arctic.

Following the niche: the differential impact of the last glacial maximum on four European ungulates

Predicting the effects of future global changes on species requires a better understanding of the ecological niche dynamics in response to climate; the large climatic fluctuations of the last 50,000 years can be used as a natural experiment to that aim. Here we test whether the realized niche of horse, aurochs, red deer, and wild boar changed between 47,000 and 7500 years ago using paleoecological modelling over an extensive archaeological database. We show that they all changed their niche, with species-specific responses to climate fluctuations. We also suggest that they survived the climatic turnovers thanks to their flexibility and by expanding their niche in response to the extinction of competitors and predators. Irrespective of the mechanism behind such processes, the fact that species with long generation times can change their niche over thousands of years cautions against assuming it to stay constant both when reconstructing the past and predicting the future.

European megafaunal ungulates living in open habitats over the last 50,000 years showed evidence for niche change, possibly driven by climatic change and extinction of competitors and predators

Ecosystem productivity affected the spatiotemporal disappearance of Neanderthals in Iberia

What role did fluctuations play in biomass availability for secondary consumers in the disappearance of Neanderthals and the survival of modern humans? To answer this, we quantify the effects of stadial and interstadial conditions on ecosystem productivity and human spatiotemporal distribution patterns during the Middle to Upper Palaeolithic transition (50,000-30,000 calibrated years before the present) in Iberia. First, we used summed probability distribution, optimal linear estimation and Bayesian age modelling to reconstruct an updated timescale for the transition. Next, we executed a generalized dynamic vegetation model to estimate the net primary productivity. Finally, we developed a macroecological model validated with present-day observations to calculate herbivore abundance. The results indicate that, in the Eurosiberian region, the disappearance of Neanderthal groups was contemporaneous with a significant decrease in the available biomass for secondary consumers, and the arrival of the first Homo sapiens populations coincided with an increase in herbivore carrying capacity. During stadials, the Mediterranean region had the most stable conditions and the highest biomass of medium and medium-large herbivores. These outcomes support an ecological cause for the hiatus between the Mousterian and Aurignacian technocomplexes in Northern Iberia and the longer persistence of Neanderthals in southern latitudes.

The recent Asian elephant range expansion in Yunnan, China, is associated with climate change and enforced protection efforts in human-dominated landscapes

Recently, the northward movement of Asian elephants (Elephas maximus) in Yunnan, China, has attracted international attention. Climate change or human disturbances have been proposed to be the key drivers, but these hypotheses have not been rigorously tested. In this study, we quantified the relationship between climate change and human impacts on the recent range expansion of Asian elephants in southwest China. We found that the first observation probability of this species in a new place during 1959–2021 had a significant and positive association with change in air temperature and human density, resulting in a movement toward a high-latitude region with a warmer climate and higher human density; however, its association with precipitation was scale-dependent in time: positive or negative during the past 10 or 5 years, respectively. Under the enforced protection policy, human-dominated areas became preferred habitats for elephants. Our results indicate that climate change and enforced protection efforts in human-dominated landscapes in the last few decades are significant drivers of the recent range expansion of Asian elephants in Yunnan, China. It is necessary to expand the current protected areas or habitat corridors toward the north or set up new reserves in the north and set up barriers between human settlements and elephant habitats to facilitate elephant movements and minimize human-elephant conflicts under accelerated global change.

The untapped potential of macrofossils in ancient plant DNA research

The rapid development of ancient DNA analysis in the last decades has induced a paradigm shift in ecology and evolution. Driven by a combination of breakthroughs in DNA isolation techniques, high-throughput sequencing, and bioinformatics, ancient genome-scale data for a rapidly growing variety of taxa are now available, allowing researchers to directly observe demographic and evolutionary processes over time. However, the vast majority of paleogenomic studies still focus on human or animal remains. In this article, we make the case for a vast untapped resource of ancient plant material that is ideally suited for paleogenomic analyses: plant remains, such as needles, leaves, wood, seeds, or fruits, that are deposited in natural archives, such as lake sediments, permafrost, or even ice caves. Such plant remains are commonly found in large numbers and in stratigraphic sequence through time and have so far been used primarily to reconstruct past local species presences and abundances. However, they are also unique repositories of genetic information with the potential to revolutionize the fields of ecology and evolution by directly studying microevolutionary processes over time. Here, we give an overview of the current state-of-the-art, address important challenges, and highlight new research avenues to inspire future research.

Birds adapted to cold conditions show greater changes in range size related to past climatic oscillations than temperate birds

Investigation of ecological responses of species to past climate oscillations provides crucial information to understand the effects of global warming. In this work, we investigated how past climate changes affected the distribution of six bird species with different climatic requirements and migratory behaviours in the Western Palearctic and in Africa. Species Distribution Models and Marine Isotopic Stage (MIS) 2 fossil occurrences of selected species were employed to evaluate the relation between changes in range size and species climatic tolerances. The Last Glacial Maximum (LGM) range predictions, generally well supported by the MIS 2 fossil occurrences, suggest that cold-dwelling species considerably expanded their distribution in the LGM, experiencing more pronounced net changes in range size compared to temperate species. Overall, the thermal niche proves to be a key ecological trait for explaining the impact of climate change in species distributions. Thermal niche is linked to range size variations due to climatic oscillations, with cold-adapted species currently suffering a more striking range reduction compared to temperate species. This work also supports the persistence of Afro-Palearctic migrations during the LGM due to the presence of climatically suitable wintering areas in Africa even during glacial maxima.

Climate Warming Since the Holocene Accelerates West-East Communication for the Eurasian Temperate Water Strider Species Aquarius paludum

Holocene climate warming has dramatically altered biological diversity and distributions. Recent human-induced emissions of greenhouse gases will exacerbate global warming and thus induce threats to cold-adapted taxa. However, the impacts of this major climate change on transcontinental temperate species are still poorly understood. Here, we generated extensive genomic datasets for a water strider, Aquarius paludum, which was sampled across its entire distribution in Eurasia and used these datasets in combination with ecological niche modeling (ENM) to elucidate the influence of the Holocene and future climate warming on its population structure and demographic history. We found that A. paludum consisted of two phylogeographic lineages that diverged in the middle Pleistocene, which resulted in a “west–east component” genetic pattern that was probably triggered by Central Asia-Mongoxin aridification and Pleistocene glaciations. The diverged western and eastern lineages had a second contact in the Holocene, which shaped a temporary hybrid zone located at the boundary of the arid–semiarid regions of China. Future predictions detected a potentially novel northern corridor to connect the western and eastern populations, indicating west–east gene flow would possibly continue to intensify under future warming climate conditions. Further integrating phylogeographic and ENM analyses of multiple Eurasian temperate taxa based on published studies reinforced our findings on the “west–east component” genetic pattern and the predicted future northern corridor for A. paludum. Our study provided a detailed paradigm from a phylogeographic perspective of how transcontinental temperate species differ from cold-adapted taxa in their response to climate warming.

Mammal extinction facilitated biome shift and human population change during the last glacial termination in East-Central Europe

The study of local extinction times, together with the associated environmental and human population changes in the last glacial termination, provides insights into the causes of mega- and microfauna extinctions. In East-Central (EC) Europe, groups of Palaeolithic humans were present throughout the last glacial maximum, but disappeared suddenly around 15,200 cal BP. In this study cave sediment profiles dated using radiocarbon techniques and a large set of mammal bones dated directly by AMS ¹⁴C were used to determine local extinction times. These were, in turn, compared to changes in the total megafauna population of EC Europe derived from coprophilous fungi, the Epigravettian population decline, quantitative climate models, pollen and plant macrofossil inferred climate, as well as to biome reconstructions. The results suggest that the population size of large herbivores decreased in the area after 17,700 cal BP, when temperate tree abundance and warm continental steppe cover both increased in the lowlands. Boreal forest expansion started around 16,200 cal BP. Cave sediments show the decline of narrow-headed vole and arctic lemming populations specifically associated with a tundra environment at the same time and the expansion of the common vole, an inhabitant of steppes. The last dated appearance of arctic lemming was at ~ 16,640 cal BP, while that of the narrow-headed vole at ~ 13,340, and the estimated extinction time of woolly mammoth was either at 13,830 (GRIWM) or 15,210 (PHASE), and reindeer at 11,860 (GRIWM) or 12,550 cal BP (PHASE). The population decline of the large herbivore fauna slightly preceded changes in terrestrial vegetation, and likely facilitated it via a reduction in the intensity of grazing and the concomitant accumulation of plant biomass. Furthermore, it is possible to conclude that the Late Epigravettian population had high degree of quarry-fidelity; they left the basin when these mammals vanished.

Strong and lasting impacts of past global warming on baleen whales and their prey

Global warming is affecting the population dynamics and trophic interactions across a wide range of ecosystems and habitats. Translating these real-time effects into their long-term consequences remains a challenge. The rapid and extreme warming period that occurred after the Last Glacial Maximum (LGM) during the Pleistocene-Holocene transition (7-12 thousand years ago) provides an opportunity to gain insights into the long-term responses of natural populations to periods with global warming. The effects of this post-LGM warming period have been assessed in many terrestrial taxa, whereas insights into the impacts of rapid global warming on marine taxa remain limited, especially for megafauna. In order to understand how large-scale climate fluctuations during the post-LGM affected baleen whales and their prey, we conducted an extensive, large-scale analysis of the long-term effects of the post-LGM warming on abundance and inter-ocean connectivity in eight baleen whale and seven prey (fish and invertebrates) species across the Southern and the North Atlantic Ocean; two ocean basins that differ in key oceanographic features. The analysis was based upon 7032 mitochondrial DNA sequences as well as genome-wide DNA sequence variation in 100 individuals. The estimated temporal changes in genetic diversity during the last 30,000 years indicated that most baleen whale populations underwent post-LGM expansions in both ocean basins. The increase in baleen whale abundance during the Holocene was associated with simultaneous changes in their prey and climate. Highly correlated, synchronized and exponential increases in abundance in both baleen whales and their prey in the Southern Ocean were indicative of a dramatic increase in ocean productivity. In contrast, the demographic fluctuations observed in baleen whales and their prey in the North Atlantic Ocean were subtle, varying across taxa and time. Perhaps most important was the observation that the ocean-wide expansions and decreases in abundance that were initiated by the post-LGM global warming, continued for millennia after global temperatures stabilized, reflecting persistent, long-lasting impacts of global warming on marine fauna.

An intermediate crocodylian linking two extant gharials from the Bronze Age of China and its human-induced extinction

A solid phylogenetic framework is the basis of biological studies, yet higher level relationships are still unresolved in some major vertebrate lineages. One such group is Crocodylia, where the branching pattern of three major families (Alligatoridae, Crocodylidae and Gavialidae) has been disputed over decades due to the uncertain relationship of two slender-snouted lineages, gavialines and tomistomines. Here, we report a bizarre crocodylian from the Bronze Age of China, which shows a mosaic of gavialine and tomistomine features across the skeleton, rendering support to their sister taxon relationship as molecular works have consistently postulated. Gavialine characters of the new Chinese crocodylian include a novel configuration of the pterygoid bulla, a vocal structure known in mature male Indian gharials. Extinct gavialines have repeatedly evolved potentially male-only acoustic apparatus of various shapes, illuminating the deep history of sexual selection on acoustic signalling in a slender-snouted group of crocodylians. Lastly, a cutmark analysis combined with accelerator mass spectrometry (AMS) radiocarbon dating of bone remains demonstrated that two individuals from Shang and Zhou dynasties in Guangdong, China, suffered head injuries and decapitation. Archaeological evidence together with historical accounts suggests the human-induced extinction of this unique crocodylian only a few hundred years ago.

Response of Iranian lizards to future climate change by poleward expansion, southern contraction, and elevation shifts

This study explores the relationships between recent Iranian lizard species distributions and the observed climate, as well as potential future distributions of species. For this purpose, an ensemble of seven algorithms was used to forecast the distributions of 30 species for the recent and future (2070) based on the averages of 14 global climate models under optimistic (RCP2.6) and pessimistic (RCP8.5) scenarios. Annual precipitation (n = 16) and annual mean temperature (n = 7) were identified as the most important variables in determining the distribution of 76.66% (23 out of 30) of the species. The consensus model predicts that the ranges of 83.33% of species (n = 25) have the potential to expand poleward at higher latitudes while preserving the majority of their recent distributions (except for four species). Furthermore, the ranges of the remaining species (n = 5) will be preserved at higher latitudes. However, they (n = 22) may contract slightly (n = 13) or excessively (n = 9) in the south of their distribution range at lower latitudes. These results indicate that species (N = 19) situated in mountainous areas such as the Zagros, Alborz, and Kopet Dagh may move or maintain their range at higher elevations as a result of future climate change. Finally, this study suggests that 30% of species (n = 9) may be threatened by future climate change and that they should be prioritized in conservation efforts.

The pre-eminent role of directional selection in generating extreme morphological change in glyptodonts (Cingulata; Xenarthra)

The prevalence of stasis on macroevolution has been classically taken as evidence of the strong role of stabilizing selection in constraining morphological change. Rates of evolution calculated over longer timescales tend to fall below the expected under genetic drift, suggesting that directional selection signals are erased at longer timescales. Here, we investigated the rates of morphological evolution of the skull in a fossil lineage that underwent extreme morphological modification, the glyptodonts. Contrary to what was expected, we show here that directional selection was the primary process during the evolution of glyptodonts. Furthermore, the reconstruction of selection patterns shows that traits selected to generate a glyptodont morphology are markedly different from those operating on extant armadillos. Changes in both direction and magnitude of selection are probably tied to glyptodonts' invasion of a specialist-herbivore adaptive zone. These results suggest that directional selection might have played a more critical role in the evolution of extreme morphologies than previously imagined.

Collapse of the mammoth-steppe in central Yukon as revealed by ancient environmental DNA

The temporal and spatial coarseness of megafaunal fossil records complicates attempts to to disentangle the relative impacts of climate change, ecosystem restructuring, and human activities associated with the Late Quaternary extinctions. Advances in the extraction and identification of ancient DNA that was shed into the environment and preserved for millennia in sediment now provides a way to augment discontinuous palaeontological assemblages. Here, we present a 30,000-year sedimentary ancient DNA (sedaDNA) record derived from loessal permafrost silts in the Klondike region of Yukon, Canada. We observe a substantial turnover in ecosystem composition between 13,500 and 10,000 calendar years ago with the rise of woody shrubs and the disappearance of the mammoth-steppe (steppe-tundra) ecosystem. We also identify a lingering signal of Equus sp. (North American horse) and Mammuthus primigenius (woolly mammoth) at multiple sites persisting thousands of years after their supposed extinction from the fossil record.

Adaptive Preferential Sampling in Phylodynamics With an Application to SARS-CoV-2

Longitudinal molecular data of rapidly evolving viruses and pathogens provide information about disease spread and complement traditional surveillance approaches based on case count data. The coalescent is used to model the genealogy that represents the sample ancestral relationships. The basic assumption is that coalescent events occur at a rate inversely proportional to the effective population size N_e(t), a time-varying measure of genetic diversity. When the sampling process (collection of samples over time) depends on N_e(t), the coalescent and the sampling processes can be jointly modeled to improve estimation of N_e(t). Failing to do so can lead to bias due to model misspecification. However, the way that the sampling process depends on the effective population size may vary over time. We introduce an approach where the sampling process is modeled as an inhomogeneous Poisson process with rate equal to the product of N_e(t) and a time-varying coefficient, making minimal assumptions on their functional shapes via Markov random field priors. We provide efficient algorithms for inference, show the model performance vis-a-vis alternative methods in a simulation study, and apply our model to SARS-CoV-2 sequences from Los Angeles and Santa Clara counties. The methodology is implemented and available in the R package adapref. Supplementary files for this article are available online.

Museum Genomics

Natural history collections are invaluable repositories of biological information that provide an unrivaled record of Earth's biodiversity. Museum genomics-genomics research using traditional museum and cryogenic collections and the infrastructure supporting these investigations-has particularly enhanced research in ecology and evolutionary biology, the study of extinct organisms, and the impact of anthropogenic activity on biodiversity. However, leveraging genomics in biological collections has exposed challenges, such as digitizing, integrating, and sharing collections data; updating practices to ensure broadly optimal data extraction from existing and new collections; and modernizing collections practices, infrastructure, and policies to ensure fair, sustainable, and genomically manifold uses of museum collections by increasingly diverse stakeholders. Museum genomics collections are poised to address these challenges and, with increasingly sensitive genomics approaches, will catalyze a future era of reproducibility, innovation, and insight made possible through integrating museum and genome sciences.

Comment on "A global environmental crisis 42,000 years ago"

Cooper et al. (Research Articles, 19 February 2021, p. 811) propose that the Laschamps geomagnetic inversion ~42 ka BP drove global climatic shifts, causing major behavioural changes within prehistoric groups, and events of human and megafaunal extinction. Other scientific studies indicate that this proposition is unproven from the current archaeological, paleoanthropological, and genetic records.

The ghost of hosts past: impacts of host extinction on parasite specificity

A growing body of research is focused on the extinction of parasite species in response to host endangerment and declines. Beyond the loss of parasite species richness, host extinction can impact apparent parasite host specificity, as measured by host richness or the phylogenetic distances among hosts. Such impacts on the distribution of parasites across the host phylogeny can have knock-on effects that may reshape the adaptation of both hosts and parasites, ultimately shifting the evolutionary landscape underlying the potential for emergence and the evolution of virulence across hosts. Here, we examine how the reshaping of host phylogenies through extinction may impact the host specificity of parasites, and offer examples from historical extinctions, present-day endangerment, and future projections of biodiversity loss. We suggest that an improved understanding of the impact of host extinction on contemporary host-parasite interactions may shed light on core aspects of disease ecology, including comparative studies of host specificity, virulence evolution in multi-host parasite systems, and future trajectories for host and parasite biodiversity. This article is part of the theme issue 'Infectious disease macroecology: parasite diversity and dynamics across the globe'.

Process-explicit models reveal pathway to extinction for woolly mammoth using pattern-oriented validation

Pathways to extinction start long before the death of the last individual. However, causes of early stage population declines and the susceptibility of small residual populations to extirpation are typically studied in isolation. Using validated process-explicit models, we disentangle the ecological mechanisms and threats that were integral in the initial decline and later extinction of the woolly mammoth. We show that reconciling ancient DNA data on woolly mammoth population decline with fossil evidence of location and timing of extinction requires process-explicit models with specific demographic and niche constraints, and a constrained synergy of climatic change and human impacts. Validated models needed humans to hasten climate-driven population declines by many millennia, and to allow woolly mammoths to persist in mainland Arctic refugia until the mid-Holocene. Our results show that the role of humans in the extinction dynamics of woolly mammoth began well before the Holocene, exerting lasting effects on the spatial pattern and timing of its range-wide extinction.

The Tien Shan vole (Microtus ilaeus; Rodentia: Cricetidae) as a new species in the Late Pleistocene of Europe

Grey voles (subgenus Microtus) represent a complex of at least seven closely related and partly cryptic species. The range of these species extends from the Atlantic to the Altai Mountains, but most of them occur east of the Black Sea. Using ancient DNA analyses of the Late Pleistocene specimens, we identified a new mtDNA lineage of grey voles in Europe. Phylogenetic analysis of mitochondrial DNA cytochrome b sequences from 23 voles from three caves, namely, Emine-Bair-Khosar (Crimea, Ukraine), Cave 16 (Bulgaria), and Bacho Kiro (Bulgaria), showed that 14 specimens form a previously unrecognized lineage, sister to the Tien Shan vole. The average sequence divergence of this lineage and the extant Tien Shan vole was 4.8%, which is similar to the divergence of grey vole forms, which are considered distinct species or being on the verge of speciation; M. arvalis and M. obscurus or M. mystacinus and M. rossiaemeridionalis. We estimated the time to the most recent common ancestor of the grey voles to be 0.66 Ma, which is over twice the recent estimates, while the divergence of the extant Tien Shan vole and the new lineage to be 0.29 Ma. Our discovery suggests that grey voles may have been more diversified in the past and that their ranges may have differed substantially from current ones. It also underlines the utility of ancient DNA to decipher the evolutionary history of voles.

Life-history traits and habitat availability shape genomic diversity in birds: implications for conservation

More than 25% of species assessed by the International Union for Conservation of Nature (IUCN) are threatened with extinction. Understanding how environmental and biological processes have shaped genomic diversity may inform management practices. Using 68 extant avian species, we parsed the effects of habitat availability and life-history traits on genomic diversity over time to provide a baseline for conservation efforts. We used published whole-genome sequence data to estimate overall genomic diversity as indicated by historical long-term effective population sizes (Ne) and current genomic variability (H), then used environmental niche modelling to estimate Pleistocene habitat dynamics for each species. We found that Ne and H were positively correlated with habitat availability and related to key life-history traits (body mass and diet), suggesting the latter contribute to the overall genomic variation. We found that H decreased with increasing species extinction risk, suggesting that H may serve as a leading indicator of demographic trends related to formal IUCN conservation status in birds. Our analyses illustrate that genome-wide summary statistics estimated from sequence data reflect meaningful ecological attributes relevant to species conservation.

A multi-proxy approach to exploring Homo sapiens' arrival, environments and adaptations in Southeast Asia

The capability of Pleistocene hominins to successfully adapt to different types of tropical forested environments has long been debated. In order to investigate environmental changes in Southeast Asia during a critical period for the turnover of hominin species, we analysed palaeoenvironmental proxies from five late Middle to Late Pleistocene faunas. Human teeth discoveries have been reported at Duoi U'Oi, Vietnam (70-60 ka) and Nam Lot, Laos (86-72 ka). However, the use of palaeoproteomics allowed us to discard the latter, and, to date, no human remains older than ~ 70 ka are documented in the area. Our findings indicate that tropical rainforests were highly sensitive to climatic changes over that period, with significant fluctuations of the canopy forests. Locally, large-bodied faunas were resilient to these fluctuations until the cooling period of the Marine Isotope Stage 4 (MIS 4; 74-59 ka) that transformed the overall biotope. Then, under strong selective pressures, populations with new phenotypic characteristics emerged while some other species disappeared. We argue that this climate-driven shift offered new foraging opportunities for hominins in a novel rainforest environment and was most likely a key factor in the settlement and dispersal of our species during MIS 4 in SE Asia.

Humans, megafauna and landscape structure – Rock engravings from Namibia encourage a comparative approach to central Europe and southern Africa

This paper deals with reflections that arose after observing prehistoric rock engravings at different locations in Namibia. These observations stimulated comparative considerations with focus on southern Africa and central Europe. Similar to the Aurignacian rock art of European origin, the most common motifs in the Namibian rock engravings are large animals. While in Europe, the species that served as a blueprint for the illustration of Aurignacian rock art have mostly disappeared, the megafauna illustrated on the rock engravings in Namibia can still be found in the immediate vicinity of the rock art. Against this background, we discuss and further develop a comparative regional approach. We reconstruct and evaluate the suitability of African savannas and still-existing megafauna communities as an appropriate reference-frame for natural European grassland systems and extinct associated warm-adapted megafauna (Eemian Interglacial megafauna). Special attention is laid on the unique situation in Africa in the light of a global extinction wave of megafauna following increasing human activity in the Late Quaternary. This leads us to discuss the use of domesticated ungulates as surrogate taxa to fulfill ecosystem functions in Europe as part of concepts termed “rewilding” or “naturalistic grazing”. After critically examining these concepts, we conclude that using domesticated forms as representatives of extinct or locally disappeared species in Europe has its justification to some extent. If, however, the naturally occurring megaherbivore community still exists (Africa), these naturally occurring species should be given priority due to their organismic abilities and limitations adapted to the harsh conditions in their specific environment. Finally, we discuss the application of (transboundary) protected areas as effective instruments to mitigate human-wildlife conflicts. A holistic approach, including nature conservation and preservation of cultural achievements (domesticated forms, grazing systems), appears promising for the effective protection of the natural African savanna ecosystems with their unique fauna elements, as illustrated in rock engravings that inspired us to write this paper.

Quantifying effects of snow depth on caribou winter range selection and movement in Arctic Alaska

BACKGROUND

Caribou and reindeer across the Arctic spend more than two thirds of their lives moving in snow. Yet snow-specific mechanisms driving their winter ecology and potentially influencing herd health and movement patterns are not well known. Integrative research coupling snow and wildlife sciences using observations, models, and wildlife tracking technologies can help fill this knowledge void.

METHODS

Here, we quantified the effects of snow depth on caribou winter range selection and movement. We used location data of Central Arctic Herd (CAH) caribou in Arctic Alaska collected from 2014 to 2020 and spatially distributed and temporally evolving snow depth data produced by SnowModel. These landscape-scale (90 m), daily snow depth data reproduced the observed spatial snow-depth variability across typical areal extents occupied by a wintering caribou during a 24-h period.

RESULTS

We found that fall snow depths encountered by the herd north of the Brooks Range exerted a strong influence on selection of two distinct winter range locations. In winters with relatively shallow fall snow depth (2016/17, 2018/19, and 2019/20), the majority of the CAH wintered on the tundra north of the Brooks Range mountains. In contrast, during the winters with relatively deep fall snow depth (2014/15, 2015/16, and 2017/18), the majority of the CAH caribou wintered in the mountainous boreal forest south of the Brooks Range. Long-term (19 winters; 2001-2020) monitoring of CAH caribou winter distributions confirmed this relationship. Additionally, snow depth affected movement and selection differently within these two habitats: in the mountainous boreal forest, caribou avoided areas with deeper snow, but when on the tundra, snow depth did not trigger significant deep-snow avoidance. In both wintering habitats, CAH caribou selected areas with higher lichen abundance, and they moved significantly slower when encountering deeper snow.

CONCLUSIONS

In general, our findings indicate that regional-scale selection of winter range is influenced by snow depth at or prior to fall migration. During winter, daily decision-making within the winter range is driven largely by snow depth. This integrative approach of coupling snow and wildlife observations with snow-evolution and caribou-movement modeling to quantify the multi-facetted effects of snow on wildlife ecology is applicable to caribou and reindeer herds throughout the Arctic.

Population dynamics of caribou shaped by glacial cycles before the last glacial maximum

Pleistocene glacial cycles influenced the diversification of high‐latitude wildlife species through recurrent periods of range contraction, isolation, divergence, and expansion from refugia and subsequent admixture of refugial populations. We investigate population size changes and the introgressive history of caribou (Rangifer tarandus) in western Canada using 33 whole genome sequences coupled with larger‐scale mitochondrial data. We found that a major population expansion of caribou occurred starting around 110,000 years ago (kya), the start of the last glacial period. Additionally, we found effective population sizes of some caribou reaching ~700,000 to 1,000,000 individuals, one of the highest recorded historical effective population sizes for any mammal species thus far. Mitochondrial analyses dated introgression events prior to the LGM dating to 20–30 kya and even more ancient at 60 kya, coinciding with colder periods with extensive ice coverage, further demonstrating the importance of glacial cycles and events prior to the LGM in shaping demographic history. Reconstructing the origins and differential introgressive history has implications for predictions on species responses under climate change. Our results have implications for other whole genome analyses using pairwise sequentially Markovian coalescent (PSMC) analyses, as well as highlighting the need to investigate pre‐LGM demographic patterns to fully reconstruct the origin of species diversity, especially for high‐latitude species.

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.