Ancient horse genomes reveal the timing and extent of dispersals across the Bering Land Bridge

Widespread horse-based mobility arose around 2200 BCE in Eurasia

Horses revolutionized human history with fast mobility1. However, the timeline between their domestication and their widespread integration as a means of transport remains contentious2-4. Here we assemble a collection of 475 ancient horse genomes to assess the period when these animals were first reshaped by human agency in Eurasia. We find that reproductive control of the modern domestic lineage emerged around 2200 BCE, through close-kin mating and shortened generation times. Reproductive control emerged following a severe domestication bottleneck starting no earlier than approximately 2700 BCE, and coincided with a sudden expansion across Eurasia that ultimately resulted in the replacement of nearly every local horse lineage. This expansion marked the rise of widespread horse-based mobility in human history, which refutes the commonly held narrative of large horse herds accompanying the massive migration of steppe peoples across Europe around 3000 BCE and earlier3,5. Finally, we detect significantly shortened generation times at Botai around 3500 BCE, a settlement from central Asia associated with corrals and a subsistence economy centred on horses6,7. This supports local horse husbandry before the rise of modern domestic bloodlines.

Comparison of Bayesian Coalescent Skyline Plot Models for Inferring Demographic Histories

Bayesian coalescent skyline plot models are widely used to infer demographic histories. The first (non-Bayesian) coalescent skyline plot model assumed a known genealogy as data, while subsequent models and implementations jointly inferred the genealogy and demographic history from sequence data, including heterochronous samples. Overall, there exist multiple different Bayesian coalescent skyline plot models which mainly differ in two key aspects: (i) how changes in population size are modeled through independent or autocorrelated prior distributions, and (ii) how many change-points in the demographic history are used, where they occur and if the number is pre-specified or inferred. The specific impact of each of these choices on the inferred demographic history is not known because of two reasons: first, not all models are implemented in the same software, and second, each model implementation makes specific choices that the biologist cannot influence. To facilitate a detailed evaluation of Bayesian coalescent skyline plot models, we implemented all currently described models in a flexible design into the software RevBayes. Furthermore, we evaluated models and choices on an empirical dataset of horses supplemented by a small simulation study. We find that estimated demographic histories can be grouped broadly into two groups depending on how change-points in the demographic history are specified (either independent of or at coalescent events). Our simulations suggest that models using change-points at coalescent events produce spurious variation near the present, while most models using independent change-points tend to over-smooth the inferred demographic history.

Characteristic of Przewalski horses population from Askania-Nova reserve based on genetic markers

Przewalski horses are considered the last living population of wild horses, however, they are secondarily feral offspring of herds domesticated ~ 5000 years ago by the Botai culture. After Przewalski horses were almost extinct at the beginning of the twentieth century, their population is about 2500 individuals worldwide, with one of the largest breeding centers in Askania-Nova Biosphere Reserve (Ukraine). The research aimed to establish the maternal variation of Przewalski horses population maintained in Askania-Nova Reserve based on mitochondrial DNA hypervariable 1 and hypervariable 2 regions profiling, as well as, analysis of Y chromosome single nucleotide polymorphism unique for Przewalski horses, and coat color markers: MC1R and TBX3. The mtDNA hypervariable regions analysis in 23 Przewalski horses allowed assigning them to three distinctly different haplotypes, showing the greatest similarity to the Equus caballus reference, the Equus przewalskii reference, and to extinct species-Haringtonhippus. The Y chromosome analysis using fluorescently labelled assays differentiated horses in terms of polymorphism (g731821T>C) characteristic of Equus przewalskii. All male individuals presented genotype C characteristics for Przewalski horses. The polymorphisms within the coat color genes indicated only native, wild genotypes. The Y chromosome and coat color analysis denied admixtures of the tested horses with other Equidae.

Ancient segmentally duplicated LCORL retrocopies in equids

LINE-1 is an active transposable element encoding proteins capable of inserting host gene retrocopies, resulting in retro-copy number variants (retroCNVs) between individuals. Here, we performed retroCNV discovery using 86 equids and identified 437 retrocopy insertions. Only 5 retroCNVs were shared between horses and other equids, indicating that the majority of retroCNVs inserted after the species diverged. A large number (17-35 copies) of segmentally duplicated Ligand Dependent Nuclear Receptor Corepressor Like (LCORL) retrocopies were present in all equids but absent from other extant perissodactyls. The majority of LCORL transcripts in horses and donkeys originate from the retrocopies. The initial LCORL retrotransposition occurred 18 million years ago (17-19 95% CI), which is coincident with the increase in body size, reduction in digit number, and changes in dentition that characterized equid evolution. Evolutionary conservation of the LCORL retrocopy segmental amplification in the Equidae family, high expression levels and the ancient timeline for LCORL retrotransposition support a functional role for this structural variant.

Global population genetic structure and lineage differentiation of the stable fly, Stomoxys calcitrans

The bloodsucking fly, Stomoxys calcitrans (Diptera: Muscidae), is a cosmopolitan pest that transmits potential pathogens mechanically. We conduct phylogeographic analyses of S. calcitrans to resolve its global population genetic structure for establishing baseline of molecular studies. Results from mitochondrial gene suggested that the major divergence of S. calcitrans predominantly occurred 0.32-0.47 million years ago (Mya) and the subsequent diversifications took place during 0.13-0.27 Mya. The Ethiopian region was deduced as the most likely origin of S. calcitrans and the Nearctic lineages were considered to have originated from Oriental or Palaearctic regions. Our results further revealed that each biogeographic region of S. calcitrans likely maintains its genetic specialty, and yet, those non-monophyletic relationships were possibly caused by ancestral retention, dispersal with mammals, long-distance migration, and the international livestock industries. Moreover, the three highly diverged Ethiopian lineages may be putative cryptic species that require clarification of their veterinary importance. Unravelling the genetic structure of stable fly and preventing gene flow among biogeographic regions through anthropogenic activities are thus pivotal in livestock industry administration, particularly genetic exchange among differentiated lineages that might lead to the consequence of ecological trait alterations.

Phylogenomics revealed migration routes and adaptive radiation timing of Holarctic malaria mosquito species of the Maculipennis Group

BACKGROUND

Phylogenetic analyses of closely related species of mosquitoes are important for better understanding the evolution of traits contributing to transmission of vector-borne diseases. Six out of 41 dominant malaria vectors of the genus Anopheles in the world belong to the Maculipennis Group, which is subdivided into two Nearctic subgroups (Freeborni and Quadrimaculatus) and one Palearctic (Maculipennis) subgroup. Although previous studies considered the Nearctic subgroups as ancestral, details about their relationship with the Palearctic subgroup, and their migration times and routes from North America to Eurasia remain controversial. The Palearctic species An. beklemishevi is currently included in the Nearctic Quadrimaculatus subgroup adding to the uncertainties in mosquito systematics.

RESULTS

To reconstruct historic relationships in the Maculipennis Group, we conducted a phylogenomic analysis of 11 Palearctic and 2 Nearctic species based on sequences of 1271 orthologous genes. The analysis indicated that the Palearctic species An. beklemishevi clusters together with other Eurasian species and represents a basal lineage among them. Also, An. beklemishevi is related more closely to An. freeborni, which inhabits the Western United States, rather than to An. quadrimaculatus, a species from the Eastern United States. The time-calibrated tree suggests a migration of mosquitoes in the Maculipennis Group from North America to Eurasia about 20-25 million years ago through the Bering Land Bridge. A Hybridcheck analysis demonstrated highly significant signatures of introgression events between allopatric species An. labranchiae and An. beklemishevi. The analysis also identified ancestral introgression events between An. sacharovi and its Nearctic relative An. freeborni despite their current geographic isolation. The reconstructed phylogeny suggests that vector competence and the ability to enter complete diapause during winter evolved independently in different lineages of the Maculipennis Group.

CONCLUSIONS

Our phylogenomic analyses reveal migration routes and adaptive radiation timing of Holarctic malaria vectors and strongly support the inclusion of An. beklemishevi into the Maculipennis Subgroup. Detailed knowledge of the evolutionary history of the Maculipennis Subgroup provides a framework for examining the genomic changes related to ecological adaptation and susceptibility to human pathogens. These genomic variations may inform researchers about similar changes in the future providing insights into the patterns of disease transmission in Eurasia.

DNA methylation-based profiling of horse archaeological remains for age-at-death and castration

Age profiling of archaeological bone assemblages can inform on past animal management practices, but is limited by the fragmentary nature of the fossil record and the lack of universal skeletal markers for age. DNA methylation clocks offer new, albeit challenging, alternatives for estimating the age-at-death of ancient individuals. Here, we take advantage of the availability of a DNA methylation clock based on 31,836 CpG sites and dental age markers in horses to assess age predictions in 84 ancient remains. We evaluate our approach using whole-genome sequencing data and develop a capture assay providing reliable estimates for only a fraction of the cost. We also leverage DNA methylation patterns to assess castration practice in the past. Our work opens for a deeper characterization of past husbandry and ritual practices and holds the potential to reveal age mortality profiles in ancient societies, once extended to human remains.

Pleistocene parades of carnivores into North America

The distribution and movement of species, broadly known as biogeography, is one of the fundamental subfields of ecology and evolutionary biology. However, significant mysteries remain about the processes that gave rise to the modern distribution of biodiversity across the globe. Over the last several decades, the genetic study of ancient and subfossil specimens has started to shed light on past migrations of some species, with a particular focus on humans and megafauna. In this issue of Molecular Ecology, Salis et al. (2021) use ancient mitogenomes and a new phylogeographic method to add an important new piece of evidence to the mystery of megafaunal migrations into North America during the Pleistocene. They found a striking synchronicity of brown bear (Ursus arctos) and lion (Panthera spp.) migrations across the Bering Land Bridge at several time points during the late Pleistocene, which highlights the lasting impact of sea level change on the prehistoric and modern dispersal of terrestrial carnivores across continents.

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.

Genetic drift drives rapid speciation of an Arctic insular endemic shrew (Sorex pribilofensis)

Episodes of Quaternary environmental change shaped the genomes of extant species, influencing their response to contemporary environments, which are changing rapidly. Island endemics are among the most vulnerable to such change, accounting for a disproportionate number of recent extinctions. To prevent extinctions and conserve island biodiversity it is vital to combine knowledge of species' ecologies with their complex evolutionary histories. The Bering Sea has a history of cyclical island isolation and reconnection, coupled with modern rates of climate change that exceed global averages. The endangered Pribilof Island shrew (Sorex pribilofensis) is endemic to St. Paul Island, Alaska, which was isolated from mainland Beringia ~14,000 years ago by rising sea levels. Using ~11,000 single nucleotide polymorphisms, 17 microsatellites and mitochondrial sequence data, we test predictions about the evolutionary processes driving shrew speciation across Beringia. Our data show considerable differentiation of S. pribilofensis from mainland sibling species, relative to levels of divergence between mainland shrews. We also find a genome-wide loss of diversity and extremely low N_e for S. pribilofensis. We then show that intraspecific genetic diversity is significantly related to interspecific divergence, and that differentiation between S. pribilofensis and other Beringian shrews is highest across loci that are fixed in S. pribilofensis, indicating that strong drift has driven differentiation of this island species. Our findings show that drift as a consequence of Arctic climate cycling can rapidly reshape insular biodiversity. Arctic island species that lack genomic diversity and have evolved in response to past climate may have limited ability to respond to modern environmental changes.

Investigating the reliability of metapodials as taxonomic Indicators for Beringian horses

The metapodials of extinct horses have long been regarded as one of the most useful skeletal elements to determine taxonomic identity. However, recent research on both extant and extinct horses has revealed the possibility for plasticity in metapodial morphology, leading to notable variability within taxa. This calls into question the reliability of metapodials in species identification, particularly for species identified from fragmentary remains. Here, we use ten measurements of metapodials from 203 specimens of four Pleistocene horse species from eastern Beringia to test whether there are significant differences in metapodial morphology that support the presence of multiple species. We then reconstruct the body masses for every specimen to assess the range in body size within each species and determine whether species differ significantly from one another in mean body mass. We find that that taxonomic groups are based largely on the overall size of the metapodial, and that all metapodial measurements are highly autocorrelated. We also find that mean body mass differs significantly among most, but not all, species. We suggest that metapodial measurements are unreliable taxonomic indicators for Beringian horses given evidence for plasticity in metapodial morphology and their clear reflection of differences in body mass. We recommend future studies use more reliable indicators of taxonomy to identify Beringian horse species, particularly from localities from which fossils of several species have been recovered.

Evolution of the Family Equidae, Subfamily Equinae, in North, Central and South America, Eurasia and Africa during the Plio-Pleistocene

Studies of horse evolution arose during the middle of the 19th century, and several hypotheses have been proposed for their taxonomy, paleobiogeography, paleoecology and evolution. The present contribution represents a collaboration of 19 multinational experts with the goal of providing an updated summary of Pliocene and Pleistocene North, Central and South American, Eurasian and African horses. At the present time, we recognize 114 valid species across these continents, plus 4 North African species in need of further investigation. Our biochronology and biogeography sections integrate Equinae taxonomic records with their chronologic and geographic ranges recognizing regional biochronologic frameworks. The paleoecology section provides insights into paleobotany and diet utilizing both the mesowear and light microscopic methods, along with calculation of body masses. We provide a temporal sequence of maps that render paleoclimatic conditions across these continents integrated with Equinae occurrences. These records reveal a succession of extinctions of primitive lineages and the rise and diversification of more modern taxa. Two recent morphological-based cladistic analyses are presented here as competing hypotheses, with reference to molecular-based phylogenies. Our contribution represents a state-of-the art understanding of Plio-Pleistocene Equus evolution, their biochronologic and biogeographic background and paleoecological and paleoclimatic contexts.

Demography and evolutionary history of grey wolf populations around the Bering Strait

Glacial and interglacial periods throughout the Pleistocene have been substantial drivers of change in species distributions. Earlier analyses suggested that modern grey wolves (Canis lupus) trace their origin to a single Late Pleistocene Beringian population that expanded east and westwards, starting c. 25,000 years ago (ya). Here, we examined the demographic and phylogeographic histories of extant populations around the Bering Strait with wolves from two inland regions of the Russian Far East (RFE) and one coastal and two inland regions of North‐western North America (NNA), genotyped for 91,327 single nucleotide polymorphisms. Our results indicated that RFE and NNA wolves had a common ancestry until c. 34,400 ya, suggesting that these populations started to diverge before the previously proposed expansion out of Beringia. Coastal and inland NNA populations diverged c. 16,000 ya, concordant with the minimum proposed date for the ecological viability of the migration route along the Pacific Northwest coast. Demographic reconstructions for inland RFE and NNA populations reveal spatial and temporal synchrony, with large historical effective population sizes that declined throughout the Pleistocene, possibly reflecting the influence of broadscale climatic changes across continents. In contrast, coastal NNA wolves displayed a consistently lower effective population size than the inland populations. Differences between the demographic history of inland and coastal wolves may have been driven by multiple ecological factors, including historical gene flow patterns, natural landscape fragmentation, and more recent anthropogenic disturbance.

Grey wolf genomic history reveals a dual ancestry of dogs

The grey wolf (Canis lupus) was the first species to give rise to a domestic population, and they remained widespread throughout the last Ice Age when many other large mammal species went extinct. Little is known, however, about the history and possible extinction of past wolf populations or when and where the wolf progenitors of the present-day dog lineage (Canis familiaris) lived^1–8. Here we analysed 72 ancient wolf genomes spanning the last 100,000 years from Europe, Siberia and North America. We found that wolf populations were highly connected throughout the Late Pleistocene, with levels of differentiation an order of magnitude lower than they are today. This population connectivity allowed us to detect natural selection across the time series, including rapid fixation of mutations in the gene IFT88 40,000–30,000 years ago. We show that dogs are overall more closely related to ancient wolves from eastern Eurasia than to those from western Eurasia, suggesting a domestication process in the east. However, we also found that dogs in the Near East and Africa derive up to half of their ancestry from a distinct population related to modern southwest Eurasian wolves, reflecting either an independent domestication process or admixture from local wolves. None of the analysed ancient wolf genomes is a direct match for either of these dog ancestries, meaning that the exact progenitor populations remain to be located.

DNA from ancient wolves spanning 100,000 years sheds light on wolves’ evolutionary history and the genomic origin of dogs.

A polar bear paleogenome reveals extensive ancient gene flow from polar bears into brown bears

Polar bears (Ursus maritimus) and brown bears (Ursus arctos) are sister species possessing distinct physiological and behavioural adaptations that evolved over the last 500,000 years. However, comparative and population genomics analyses have revealed that several extant and extinct brown bear populations have relatively recent polar bear ancestry, probably as the result of geographically localized instances of gene flow from polar bears into brown bears. Here, we generate and analyse an approximate 20X paleogenome from an approximately 100,000-year-old polar bear that reveals a massive prehistoric admixture event, which is evident in the genomes of all living brown bears. This ancient admixture event was not visible from genomic data derived from living polar bears. Like more recent events, this massive admixture event mainly involved unidirectional gene flow from polar bears into brown bears and occurred as climate changes caused overlap in the ranges of the two species. These findings highlight the complex reticulate paths that evolution can take within a regime of radically shifting climate.

Radiocarbon and genomic evidence for the survival of Equus Sussemionus until the late Holocene

The exceptionally rich fossil record available for the equid family has provided textbook examples of macroevolutionary changes. Horses, asses, and zebras represent three extant subgenera of Equus lineage, while the Sussemionus subgenus is another remarkable Equus lineage ranging from North America to Ethiopia in the Pleistocene. We sequenced 26 archaeological specimens from Northern China in the Holocene that could be assigned morphologically and genetically to Equus ovodovi, a species representative of Sussemionus. We present the first high-quality complete genome of the Sussemionus lineage, which was sequenced to 13.4× depth of coverage. Radiocarbon dating demonstrates that this lineage survived until ~3500 years ago, despite continued demographic collapse during the Last Glacial Maximum and the great human expansion in East Asia. We also confirmed the Equus phylogenetic tree and found that Sussemionus diverged from the ancestor of non-caballine equids ~2.3–2.7 million years ago and possibly remained affected by secondary gene flow post-divergence. We found that the small genetic diversity, rather than enhanced inbreeding, limited the species’ chances of survival. Our work adds to the growing literature illustrating how ancient DNA can inform on extinction dynamics and the long-term resilience of species surviving in cryptic population pockets.

Pleistocene mitogenomes reconstructed from the environmental DNA of permafrost sediments

Traditionally, paleontologists have relied on the morphological features of bones and teeth to reconstruct the evolutionary relationships of extinct animals.¹ In recent decades, the analysis of ancient DNA recovered from macrofossils has provided a powerful means to evaluate these hypotheses and develop novel phylogenetic models.² Although a great deal of life history data can be extracted from bones, their scarcity and associated biases limit their information potential. The paleontological record of Beringia³-the unglaciated areas and former land bridge between northeast Eurasia and northwest North America-is relatively robust thanks to its perennially frozen ground favoring fossil preservation.^4,5 However, even here, the macrofossil record is significantly lacking in small-bodied fauna (e.g., rodents and birds), whereas questions related to migration and extirpation, even among well-studied taxa, remain crudely resolved. The growing sophistication of ancient environmental DNA (eDNA) methods have allowed for the identification of species within terrestrial/aquatic ecosystems,^6-12 in paleodietary reconstructions,^13-19 and facilitated genomic reconstructions from cave contexts.^8,20-22 Murchie et al.^6,23 used a capture enrichment approach to sequence a diverse range of faunal and floral DNA from permafrost silts deposited during the Pleistocene-Holocene transition.²⁴ Here, we expand on their work with the mitogenomic assembly and phylogenetic placement of Equus caballus (caballine horse), Bison priscus (steppe bison), Mammuthus primigenius (woolly mammoth), and Lagopus lagopus (willow ptarmigan) eDNA from multiple permafrost cores spanning the last 30,000 years. We identify a diverse metagenomic spectra of Pleistocene fauna and identify the eDNA co-occurrence of distinct Eurasian and American mitogenomic lineages.

The Sequence Analysis of Mitochondrial DNA Revealed Some Major Centers of Horse Domestications: The Archaeologist's Cut

The question about the time and the place of horse domestication, a process which had a profound impact on the progress of mankind, is disputable. According to the most widely accepted hypothesis, the earliest domestication of the horse happened in the western parts of the Eurasian steppes, between the Northern Black Sea region and present-day Kazakhstan and Turkmenistan. It seems that it occurred not earlier than the first half and most probably during the middle (even the last third) of the fourth millennium BC (from ∼ 5.5 kya). The next steps of large-scale horse breeding occurred almost simultaneously in Eurasia and North Africa due to the development of the social structure of human communities. On the other hand, the morphological differences between wild and domestic animals are rather vague and the genetic introgression between them is speculative. In this review, we have tried to gather all available scientific data on the existing possible hypotheses for the earliest domestication of the horse, as well as to highlight some data on the most plausible ones. This is due to the frequency of some significant data on the frequency of strictly defined mitotypes in different historical periods of human civilizations existing in the same periods.

Silver Linings at the Dawn of a “Golden Age”

Nearly four decades after the first applications of strontium isotope analyses in archaeology and paleoecology research, it could be said that we are entering a “Golden Age”. Here, we reflect on major past developments and current strengths in strontium isotope research, as well as speculate on future directions. We review (1) the currently limited number of (but much needed) controlled feeding experiments, (2) recent advances in isoscape mapping and spatial assignment, and (3) the strength of multi-proxy approaches (including both the integration of strontium isotopes with other isotope systems and complementary techniques such as ancient DNA analyses). We also explore the integration of strontium isotope research with other types of paleoecological or archaeology data, as well as with evidence and interpretative frameworks from other fields (such as conservation ecology, conservation paleobiology or history). This blending is critical as we seek to advance the field beyond simply distinguishing local or relatively sedentary individuals from those that were non-local or highly mobile. We finish with a call for future research centered on balancing methodological developments and novel applications with critical self-reflection, deeper theoretical considerations and cross-disciplinarity.

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.