Partial genomic survival of cave bears in living brown bears

The rise of baobab trees in Madagascar

The baobab trees (genus Adansonia) have attracted tremendous attention because of their striking shape and distinctive relationships with fauna1. These spectacular trees have also influenced human culture, inspiring innumerable arts, folklore and traditions. Here we sequenced genomes of all eight extant baobab species and argue that Madagascar should be considered the centre of origin for the extant lineages, a key issue in their evolutionary history2,3. Integrated genomic and ecological analyses revealed the reticulate evolution of baobabs, which eventually led to the species diversity seen today. Past population dynamics of Malagasy baobabs may have been influenced by both interspecific competition and the geological history of the island, especially changes in local sea levels. We propose that further attention should be paid to the conservation status of Malagasy baobabs, especially of Adansonia suarezensis and Adansonia grandidieri, and that intensive monitoring of populations of Adansonia za is required, given its propensity for negatively impacting the critically endangered Adansonia perrieri.

Recent increase in species-wide diversity after interspecies introgression in the highly endangered Iberian lynx

Genetic diversity is lost in small and isolated populations, affecting many globally declining species. Interspecific admixture events can increase genetic variation in the recipient species' gene pool, but empirical examples of species-wide restoration of genetic diversity by admixture are lacking. Here we present multi-fold coverage genomic data from three ancient Iberian lynx (Lynx pardinus) approximately 2,000-4,000 years old and show a continuous or recurrent process of interspecies admixture with the Eurasian lynx (Lynx lynx) that increased modern Iberian lynx genetic diversity above that occurring millennia ago despite its recent demographic decline. Our results add to the accumulating evidence for natural admixture and introgression among closely related species and show that this can result in an increase of species-wide genetic diversity in highly genetically eroded species. The strict avoidance of interspecific sources in current genetic restoration measures needs to be carefully reconsidered, particularly in cases where no conspecific source population exists.

The origins and diversification of Holarctic brown bear populations inferred from genomes of past and present populations

The brown bear (Ursus arctos) is one of the survivors of the Late Quaternary megafauna extinctions. However, despite being widely distributed across the Holarctic, brown bears have experienced extensive range reductions, and even extirpations in some geographical regions. Previous research efforts using genetic data have provided valuable insights into their evolutionary history. However, most studies have been limited to contemporary individuals or mitochondrial DNA, limiting insights into population processes that preceded the present. Here, we present genomic data from two Late Pleistocene brown bears from Honshu, Japan and eastern Siberia, and combine them with published contemporary and ancient genomes from across the Holarctic range of brown bears to investigate the evolutionary relationships among brown bear populations through time and space. By including genomic data from Late Pleistocene and Holocene individuals sampled outside the current distribution range, we uncover diversity not present in contemporary populations. Notably, although contemporary individuals display geographically structured populations most likely driven by isolation-by-distance, this pattern varies among the ancient samples across different regions. The inclusion of ancient brown bears in our analysis provides novel insights into the evolutionary history of brown bears and contributes to understanding the populations and diversity lost during the Late Quaternary.

Evolution of the Mutation Spectrum Across a Mammalian Phylogeny

Although evolutionary biologists have long theorized that variation in DNA repair efficacy might explain some of the diversity of lifespan and cancer incidence across species, we have little data on the variability of normal germline mutagenesis outside of humans. Here, we shed light on the spectrum and etiology of mutagenesis across mammals by quantifying mutational sequence context biases using polymorphism data from thirteen species of mice, apes, bears, wolves, and cetaceans. After normalizing the mutation spectrum for reference genome accessibility and k-mer content, we use the Mantel test to deduce that mutation spectrum divergence is highly correlated with genetic divergence between species, whereas life history traits like reproductive age are weaker predictors of mutation spectrum divergence. Potential bioinformatic confounders are only weakly related to a small set of mutation spectrum features. We find that clock-like mutational signatures previously inferred from human cancers cannot explain the phylogenetic signal exhibited by the mammalian mutation spectrum, despite the ability of these signatures to fit each species' 3-mer spectrum with high cosine similarity. In contrast, parental aging signatures inferred from human de novo mutation data appear to explain much of the 1-mer spectrum's phylogenetic signal in combination with a novel mutational signature. We posit that future models purporting to explain the etiology of mammalian mutagenesis need to capture the fact that more closely related species have more similar mutation spectra; a model that fits each marginal spectrum with high cosine similarity is not guaranteed to capture this hierarchy of mutation spectrum variation among species.

Large mammals and trees in eastern monsoonal China: anthropogenic losses since the Late Pleistocene and restoration prospects in the Anthropocene

Massive human-induced declines of large-sized animals and trees (megabiota) from the Late Pleistocene to the Anthropocene have resulted in downsized ecosystems across the globe, in which components and functions have been greatly simplified. In response, active restoration projects of extant large-sized species or functional substitutes are needed at large scales to promote ecological processes that are important for ecosystem self-regulation and biodiversity maintenance. Despite the desired global scope of such projects, they have received little attention in East Asia. Here, we synthesise the biogeographical and ecological knowledge of megabiota in ancient and modern China, with relevant data mostly located in eastern monsoonal China (EMC), aiming to assess its potential for restoring functionally intact ecosystems modulated by megabiota. We found that during the Late Pleistocene, 12 mammalian megafaunal (carnivores ≥15 kg and herbivores ≥500 kg) species disappeared from EMC: one carnivore Crocuta ultima (East Asian spotted hyena) and 11 herbivores including six megaherbivores (≥1000 kg). The relative importance of climate change and humans in driving these losses remains debated, despite accumulating evidence in favour of the latter. Later massive depletion of megafauna and large-sized (45-500 kg) herbivores has been closely associated with agricultural expansion and societal development, especially during the late Holocene. While forests rich in large timber trees (33 taxa in written records) were common in the region 2000-3000 years ago, millennial-long logging has resulted in considerable range contractions and at least 39 threatened species. The wide distribution of C. ultima, which likely favoured open or semi-open habitats (like extant spotted hyenas), suggests the existence of mosaic open and closed vegetation in the Late Pleistocene across EMC, in line with a few pollen-based vegetation reconstructions and potentially, or at least partially, reflecting herbivory by herbivorous megafauna. The widespread loss of megaherbivores may have strongly compromised seed dispersal for both megafruit (fleshy fruits with widths ≥40 mm) and non-megafruit plant species in EMC, especially in terms of extra-long-distance (>10 km) dispersal, which is critical for plant species that rely on effective biotic agents to track rapid climate change. The former occurrence of large mammals and trees have translated into rich material and non-material heritages passed down across generations. Several reintroduction projects have been implemented or are under consideration, with the case of Elaphurus davidianus a notable success in recovering wild populations in the middle reaches of the Yangtze River, although trophic interactions with native carnivorous megafauna have not yet been restored. Lessons of dealing with human-wildlife conflicts are key to public support for maintaining landscapes shared with megafauna and large herbivores in the human-dominated Anthropocene. Meanwhile, potential human-wildlife conflicts, e.g. public health risks, need to be scientifically informed and effectively reduced. The Chinese government's strong commitment to improved policies of ecological protection and restoration (e.g. ecological redlines and national parks) provides a solid foundation for a scaling-up contribution to the global scope needed for solving the crisis of biotic downsizing and ecosystem degradation.

Range-wide evolutionary relationships and historical demography of brown bears (Ursus arctos) revealed by whole-genome sequencing of isolated central Asian populations

Phylogeographic studies uncover hidden pathways of divergence and inform conservation. Brown bears (Ursus arctos) have one of the broadest distributions of all land mammals, ranging from Eurasia to North America, and are an important model for evolutionary studies. Although several whole genomes were available for individuals from North America, Europe and Asia, limited whole-genome data were available from Central Asia, including the highly imperilled brown bears in the Gobi Desert. To fill this knowledge gap, we sequenced whole genomes from nine Asian brown bears from the Gobi Desert of Mongolia, Northern Mongolia and the Himalayas of Pakistan. We combined these data with published brown bear sequences from Europe, Asia and North America, as well as other bear species. Our goals were to determine the evolutionary relationships among brown bear populations worldwide, their genetic diversity and their historical demography. Our analyses revealed five major lineages of brown bears based on a filtered set of 684,081 single nucleotide polymorphisms. We found distinct evolutionary lineages of brown bears in the Gobi, Himalayas, northern Mongolia, Europe and North America. The lowest level of genetic diversity and the highest level of inbreeding were found in Pakistan, the Gobi Desert and Central Italy. Furthermore, the effective population size (Ne ) for all brown bears decreased over the last 70,000 years. Our results confirm the genetic distinctiveness and ancient lineage of brown bear subspecies in the Gobi Desert of Mongolia and the Himalayas of Pakistan and highlight their importance for conservation.

Introgression between highly divergent fungal sister species

To understand how species evolve and adapt to changing environments, it is important to study gene flow and introgression due to their influence on speciation and radiation events. Here, we apply a novel experimental system for investigating these mechanisms using natural populations. The system is based on two fungal sister species with morphological and ecological similarities occurring in overlapping habitats. We examined introgression between these species by conducting whole genome sequencing of individuals from populations in North America and Europe. We assessed genome-wide nucleotide divergence and performed crossing experiments to study reproductive barriers. We further used ABBA-BABA statistics together with a network analysis to investigate introgression, and conducted demographic modelling to gain insight into divergence times and introgression events. The results revealed that the species are highly divergent and incompatible in vitro. Despite this, small regions of introgression were scattered throughout the genomes and one introgression event likely involves a ghost population (extant or extinct). This study demonstrates that introgression can be found among divergent species and that population histories can be studied without collections of all the populations involved. Moreover, the experimental system is shown to be a useful tool for research on reproductive isolation in natural populations.

Multiple introgressions shape mitochondrial evolutionary history in Drosophila paulistorum and the Drosophila willistoni group

Hybridization and the consequent introgression of genomic elements is an important source of genetic diversity for biological lineages. This is particularly evident in young clades in which hybrid incompatibilities are still incomplete and mixing between species is more likely to occur. Drosophila paulistorum, a representative of the Neotropical Drosophila willistoni subgroup, is a classic model of incipient speciation. The species is divided into six semispecies that show varying degrees of pre- and post-mating incompatibility with each other. In the present study, we investigate the mitochondrial evolutionary history of D. paulistorum and the willistoni subgroup. For that, we perform phylogenetic and comparative analyses of the complete mitochondrial genomes and draft nuclear assemblies of 25 Drosophila lines of the willistoni and saltans species groups. Our results show that the mitochondria of D. paulistorum are polyphyletic and form two non-sister clades that we name α and β. Identification and analyses of nuclear mitochondrial insertions further reveal that the willistoni subgroup has an α-like mitochondrial ancestor and strongly suggest that both the α and β mitochondria of D. paulistorum were acquired through introgression from unknown fly lineages of the willistoni subgroup. We also uncover multiple mitochondrial introgressions across D. paulistorum semispecies and generate novel insight into the evolution of the species.

Advancements and Challenges in Ancient DNA Research: Bridging the Global North-South Divide

Ancient DNA (aDNA) research first began in 1984 and ever since has greatly expanded our understanding of evolution and migration. Today, aDNA analysis is used to solve various puzzles about the origin of mankind, migration patterns, and the spread of infectious diseases. The incredible findings ranging from identifying the new branches within the human family to studying the genomes of extinct flora and fauna have caught the world by surprise in recent times. However, a closer look at these published results points out a clear Global North and Global South divide. Therefore, through this research, we aim to emphasize encouraging better collaborative opportunities and technology transfer to support researchers in the Global South. Further, the present research also focuses on expanding the scope of the ongoing conversation in the field of aDNA by reporting relevant literature published around the world and discussing the advancements and challenges in the field.

Range-wide whole-genome resequencing of the brown bear reveals drivers of intraspecies divergence

Population-genomic studies can shed new light on the effect of past demographic processes on contemporary population structure. We reassessed phylogeographical patterns of a classic model species of postglacial recolonisation, the brown bear (Ursus arctos), using a range-wide resequencing dataset of 128 nuclear genomes. In sharp contrast to the erratic geographical distribution of mtDNA and Y-chromosomal haplotypes, autosomal and X-chromosomal multi-locus datasets indicate that brown bear population structure is largely explained by recent population connectivity. Multispecies coalescent based analyses reveal cases where mtDNA haplotype sharing between distant populations, such as between Iberian and southern Scandinavian bears, likely results from incomplete lineage sorting, not from ancestral population structure (i.e., postglacial recolonisation). However, we also argue, using forward-in-time simulations, that gene flow and recombination can rapidly erase genomic evidence of former population structure (such as an ancestral population in Beringia), while this signal is retained by Y-chromosomal and mtDNA, albeit likely distorted. We further suggest that if gene flow is male-mediated, the information loss proceeds faster in autosomes than in X chromosomes. Our findings emphasise that contemporary autosomal genetic structure may reflect recent population dynamics rather than postglacial recolonisation routes, which could contribute to mtDNA and Y-chromosomal discordances.

Methodological Changes in the Field of Paleogenetics

Paleogenetics has significantly changed since its inception almost forty years ago. Initially, molecular techniques available to the researchers offered minimal possibilities for ancient DNA analysis. The subsequent expansion of the scientific tool cabinet allowed for more remarkable achievements, combined has with the newfound popularity of this budding field of science. Finally, a breakthrough was made with the development of next-generation sequencing (NGS) technologies and the update of DNA isolation protocols, through which even very fragmented aDNA samples could be used to sequence whole genomes. In this paper, we review the achievements made thus far and compare the methodologies utilized in this field of science, discussing their benefits and challenges.

Blue Turns to Gray: Paleogenomic Insights into the Evolutionary History and Extinction of the Blue Antelope (Hippotragus leucophaeus)

The blue antelope (Hippotragus leucophaeus) is the only large African mammal species to have become extinct in historical times, yet no nuclear genomic information is available for this species. A recent study showed that many alleged blue antelope museum specimens are either roan (Hippotragus equinus) or sable (Hippotragus niger) antelopes, further reducing the possibilities for obtaining genomic information for this extinct species. While the blue antelope has a rich fossil record from South Africa, climatic conditions in the region are generally unfavorable to the preservation of ancient DNA. Nevertheless, we recovered two blue antelope draft genomes, one at 3.4× mean coverage from a historical specimen (∼200 years old) and one at 2.1× mean coverage from a fossil specimen dating to 9,800-9,300 cal years BP, making it currently the oldest paleogenome from Africa. Phylogenomic analyses show that blue and sable antelope are sister species, confirming previous mitogenomic results, and demonstrate ancient gene flow from roan into blue antelope. We show that blue antelope genomic diversity was much lower than in roan and sable antelope, indicative of a low population size since at least the early Holocene. This supports observations from the fossil record documenting major decreases in the abundance of blue antelope after the Pleistocene-Holocene transition. Finally, the persistence of this species throughout the Holocene despite low population size suggests that colonial-era human impact was likely the decisive factor in the blue antelope's extinction.

Ancient DNA Reveals Maternal Philopatry of the Northeast Eurasian Brown Bear (Ursus arctos) Population during the Holocene

Significant palaeoecological and paleoclimatic changes that took place during Late Pleistocene-Early Holocene transition are considered important factors that led to megafauna extinctions. Unlike many other species, the brown bear (Ursus arctos) has survived this geological time. Despite the fact that several mitochondrial DNA clades of brown bears became extinct at the end of the Pleistocene, this species is still widely distributed in Northeast Eurasia. Here, using the ancient DNA analysis of a brown bear individual that inhabited Northeast Asia in the Middle Holocene (3460 ± 40 years BP) and comparative phylogenetic analysis, we show a significant mitochondrial DNA similarity of the studied specimen with modern brown bears inhabiting Yakutia and Chukotka. In this study, we clearly demonstrate the maternal philopatry of the Northeastern Eurasian U. arctos population during the several thousand years of the Holocene.

Ghost lineages can invalidate or even reverse findings regarding gene flow

Introgression, endosymbiosis, and gene transfer, i.e., horizontal gene flow (HGF), are primordial sources of innovation in all domains of life. Our knowledge on HGF relies on detection methods that exploit some of its signatures left on extant genomes. One of them is the effect of HGF on branch lengths of constructed phylogenies. This signature has been formalized in statistical tests for HGF detection and used for example to detect massive adaptive gene flows in malaria vectors or to order evolutionary events involved in eukaryogenesis. However, these studies rely on the assumption that ghost lineages (all unsampled extant and extinct taxa) have little influence. We demonstrate here with simulations and data reanalysis that when considering the more realistic condition that unsampled taxa are legion compared to sampled ones, the conclusion of these studies become unfounded or even reversed. This illustrates the necessity to recognize the existence of ghosts in evolutionary studies.

A beary good genome: Haplotype-resolved, chromosome-level assembly of the brown bear (Ursus arctos)

The brown bear (Ursus arctos) is the second largest and most widespread extant terrestrial carnivore on Earth and has recently emerged as a medical model for human metabolic diseases. Here, we report a fully phased chromosome-level assembly of a male North American brown bear built by combining Pacific Biosciences (PacBio) HiFi data and publicly available Hi-C data. The final genome size is 2.47 Gigabases (Gb) with a scaffold and contig N50 length of 70.08 and 43.94 Megabases (Mb), respectively. Benchmarking Universal Single-Copy Ortholog (BUSCO) analysis revealed that 94.5% of single copy orthologs from Mammalia were present in the genome (the highest of any ursid genome to date). Repetitive elements accounted for 44.48% of the genome and a total of 20,480 protein coding genes were identified. Based on whole genome alignment to the polar bear, the brown bear is highly syntenic with the polar bear, and our phylogenetic analysis of 7,246 single-copy orthologs supports the currently proposed species tree for Ursidae. This highly contiguous genome assembly will support future research on both the evolutionary history of the bear family and the physiological mechanisms behind hibernation, the latter of which has broad medical implications.

Uncovering the enigmatic evolution of bears in greater depth: The hybrid origin of the Asiatic black bear

Bears are fascinating mammals because of their complex pattern of speciation and rapid evolution of distinct phenotypes. Interspecific hybridization has been common and has shaped the complex evolutionary history of bears. In this study, based on the largest population-level genomic dataset to date involving all Ursinae species and recently developed methods for detecting hybrid speciation, we provide explicit evidence for the hybrid origin of Asiatic black bears, which arose through historical hybridization between the ancestor of polar bear/brown bear/American black bears and the ancestor of sun bear/sloth bears. This was inferred to have occurred soon after the divergence of the two parental lineages in Eurasia due to climate-driven population expansion and dispersal. In addition, we found that the intermediate body size of this hybrid species arose from its combination of relevant genes derived from two parental lineages of contrasting sizes. This and alternate fixation of numerous other loci that had diverged between parental lineages may have initiated the reproductive isolation of the Asiatic black bear from its two parents. Our study sheds further light on the evolutionary history of bears and documents the importance of hybridization in new species formation and phenotypic evolution in mammals.

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.

Insights into bear evolution from a Pleistocene polar bear genome

The polar bear (Ursus maritimus) has become a symbol of the threat to biodiversity from climate change. Understanding polar bear evolutionary history may provide insights into apex carnivore responses and prospects during periods of extreme environmental perturbations. In recent years, genomic studies have examined bear speciation and population history, including evidence for ancient admixture between polar bears and brown bears (Ursus arctos). Here, we extend our earlier studies of a 130,000- to 115,000-y-old polar bear from the Svalbard Archipelago using a 10× coverage genome sequence and 10 new genomes of polar and brown bears from contemporary zones of overlap in northern Alaska. We demonstrate a dramatic decline in effective population size for this ancient polar bear’s lineage, followed by a modest increase just before its demise. A slightly higher genetic diversity in the ancient polar bear suggests a severe genetic erosion over a prolonged bottleneck in modern polar bears. Statistical fitting of data to alternative admixture graph scenarios favors at least one ancient introgression event from brown bears into the ancestor of polar bears, possibly dating back over 150,000 y. Gene flow was likely bidirectional, but allelic transfer from brown into polar bear is the strongest detected signal, which contrasts with other published work. These findings may have implications for our understanding of climate change impacts: Polar bears, a specialist Arctic lineage, may not only have undergone severe genetic bottlenecks but also been the recipient of generalist, boreal genetic variants from brown bears during critical phases of Northern Hemisphere glacial oscillations.

Genomic variation from an extinct species is retained in the extant radiation following speciation reversal

Ecosystem degradation and biodiversity loss are major global challenges. When reproductive isolation between species is contingent on the interaction of intrinsic lineage traits with features of the environment, environmental change can weaken reproductive isolation and result in extinction through hybridization. By this process called speciation reversal, extinct species can leave traces in genomes of extant species through introgressive hybridization. Using historical and contemporary samples, we sequenced all four species of an Alpine whitefish radiation before and after anthropogenic lake eutrophication and the associated loss of one species through speciation reversal. Despite the extinction of this taxon, substantial fractions of its genome, including regions shaped by positive selection before eutrophication, persist within surviving species as a consequence of introgressive hybridization during eutrophication. Given the prevalence of environmental change, studying speciation reversal and its genomic consequences provides fundamental insights into evolutionary processes and informs biodiversity conservation.

Phylogeography of ancient and modern brown bears from eastern Eurasia

The brown bear (Ursus arctos) is an iconic carnivoran species of the Northern Hemisphere. Its population history has been studied extensively using mitochondrial markers, which demonstrated signatures of multiple waves of migration, arguably connected with glaciation periods. Among Eurasian brown bears, Siberian populations remain understudied. We have sequenced complete mitochondrial genomes of four ancient (~4.5-40 kya) bears from South Siberia and 19 modern bears from South Siberia and the Russian Far East. Reconstruction of phylogenetic relationships between haplotypes and evaluation of modern population structure have demonstrated that all the studied samples belong to the most widespread Eurasian clade 3. One of the ancient haplotypes takes a basal position relative to the whole of clade 3; the second is basal to the haplogroup 3a (the most common subclade), and two others belong to clades 3a1 and 3b. Modern Siberian bears retain at least some of this diversity; apart from the most common haplogroup 3a, we demonstrate the presence of clade 3b, which was previously found mainly in mainland Eurasia and Northern Japan. Our findings highlight the importance of South Siberia as a refugium for northern Eurasian brown bears and further corroborate the hypothesis of several waves of migration in the Pleistocene.

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.

OUP accepted manuscript

Most species are extinct, those that are not are often unknown. Sequenced and sampled species are often a minority of known ones. Past evolutionary events involving horizontal gene flow, such as horizontal gene transfer, hybridization, introgression, and admixture, are therefore likely to involve “ghosts,” that is extinct, unknown, or unsampled lineages. The existence of these ghost lineages is widely acknowledged, but their possible impact on the detection of gene flow and on the identification of the species involved is largely overlooked. It is generally considered as a possible source of error that, with reasonable approximation, can be ignored. We explore the possible influence of absent species on an evolutionary study by quantifying the effect of ghost lineages on introgression as detected by the popular D-statistic method. We show from simulated data that under certain frequently encountered conditions, the donors and recipients of horizontal gene flow can be wrongly identified if ghost lineages are not taken into account. In particular, having a distant outgroup, which is usually recommended, leads to an increase in the error probability and to false interpretations in most cases. We conclude that introgression from ghost lineages should be systematically considered as an alternative possible, even probable, scenario. [ABBA–BABA; D-statistic; gene flow; ghost lineage; introgression; simulation.]

Consequences of Hybridization in Mammals: A Systematic Review

Hybridization, defined as breeding between two distinct taxonomic units, can have an important effect on the evolutionary patterns in cross-breeding taxa. Although interspecific hybridization has frequently been considered as a maladaptive process, which threatens species genetic integrity and survival via genetic swamping and outbreeding depression, in some cases hybridization can introduce novel adaptive variation and increase fitness. Most studies to date focused on documenting hybridization events and analyzing their causes, while relatively little is known about the consequences of hybridization and its impact on the parental species. To address this knowledge gap, we conducted a systematic review of studies on hybridization in mammals published in 2010-2021, and identified 115 relevant studies. Of 13 categories of hybridization consequences described in these studies, the most common negative consequence (21% of studies) was genetic swamping and the most common positive consequence (8%) was the gain of novel adaptive variation. The total frequency of negative consequences (49%) was higher than positive (13%) and neutral (38%) consequences. These frequencies are biased by the detection possibilities of microsatellite loci, the most common genetic markers used in the papers assessed. As negative outcomes are typically easier to demonstrate than positive ones (e.g., extinction vs hybrid speciation), they may be over-represented in publications. Transition towards genomic studies involving both neutral and adaptive variation will provide a better insight into the real impacts of hybridization.

Ancient hybridization patterns between bighorn and thinhorn sheep

Whole-genome sequencing has advanced the study of species evolution, including the detection of genealogical discordant events such as ancient hybridization and incomplete lineage sorting (ILS). The evolutionary history of bighorn (Ovis canadensis) and thinhorn (Ovis dalli) sheep present an ideal system to investigate evolutionary discordance due to their recent and rapid radiation and putative secondary contact between bighorn and thinhorn sheep subspecies, specifically the dark pelage Stone sheep (O. dalli stonei) and predominately white Dall sheep (O. dalli dalli), during the last ice age. Here, we used multiple genomes of bighorn and thinhorn sheep, together with snow (O. nivicola) and the domestic sheep (O. aries) as outgroups, to assess their phylogenomic history, potential introgression patterns and their adaptive consequences. Among the Pachyceriforms (snow, bighorn and thinhorn sheep) a consistent monophyletic species tree was retrieved; however, many genealogical discordance patterns were observed. Alternative phylogenies frequently placed Stone and bighorn as sister clades. This relationship occurred more often and was less divergent than that between Dall and bighorn. We also observed many blocks containing introgression signal between Stone and bighorn genomes in which coat colour genes were present. Introgression signals observed between Dall and bighorn were more random and less frequent, and therefore probably due to ILS or intermediary secondary contact. These results strongly suggest that Stone sheep originated from a complex series of events, characterized by multiple, ancient periods of secondary contact with bighorn sheep.

Finding unknown species in the genomes of extant species

Although many species have gone extinct, their genetic components might exist in extant species because of ancient hybridization. Via advances in genome sequencing and development of modern population genetics, one can find the legacy of unknown or extinct species in the context of available genomes from extant species. Such discovery can be used as a strategy to search for hidden species or fossils in conservation biology and archeology, gain novel insight into complex evolutionary history, and provide the new sources of genetic variation for breeding and trait improvement in agriculture.

Demographic History of the Brown Bear (Ursus arctos) on Hokkaido Island, Japan, Based on Whole-Genomic Sequence Analysis

Previous studies of the brown bear (Ursus arctos) on Hokkaido Island, Japan, have detected three geographically distinct subpopulations representing different mitochondrial lineages and shown that gene flow between subpopulations has occurred due to male-biased dispersal. In this study, we determined whole-genomic sequences for six Hokkaido brown bears and analyzed these data along with previously published genomic sequences of 17 brown bears from other parts of the world. We found that the Hokkaido population is genetically distinct from the other populations, keeping genetic diversity higher than the endangered populations in western Europe but lower than most populations on the continents. A reconstruction of historical demography showed no increase in population size for the Hokkaido population during the Eemian interglacial period (130,000–114,000 years ago). In a phylogenetic analysis of the autosomal data, the Hokkaido population formed a clade distinct from North American and European populations, showing that it has maintained genetic diversity independently from continental populations following geographical isolation on the island. This autosomal genetic similarity contrasts with the geographically separate mitochondrial lineages on Hokkaido and indicates the occurrence of male-driven gene flow between subpopulations.

Genomic introgression from a distant congener in the Levant fritillary butterfly, Melitaea acentria

Introgressive hybridization is more common in nature than previously thought, and its role and creative power in evolution is hotly discussed but not completely understood. Introgression occurs more frequently in sympatry between recently diverged taxa, or when the speciation process has not yet been completed. However, there are relatively few documented cases of hybridization that erodes reproductive barriers between distantly related species. Here, we use whole genome and mitochondrial data to examine how introgression from a distant congener affects pattern of genetic differentiation in the Levant fritillary butterfly Melitaea acentria. We show that this local taxon has evolved as a peripatric geographic isolate of the widespread Melitaea persea, and that there has been significant unidirectional gene flow from the sympatric, nonclosely related Melitaea didyma to M. acentria. We found direct evidence of ongoing sporadic hybridization between M. didyma and M. acentria, which are separated by at least 5 million years of independent evolution. Elevated differentiation and lower level of introgression on the sex Z chromosome compared to autosomes suggest that the Z chromosome has accumulated loci acting as intrinsic postzygotic barriers. Our results show that introgression from M. didyma has been an additional source of nucleotide diversity in the M. acentria population, providing material for drift and selection.

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

The Bering Land Bridge (BLB) last connected Eurasia and North America during the Late Pleistocene. Although the BLB would have enabled transfers of terrestrial biota in both directions, it also acted as an ecological filter whose permeability varied considerably over time. Here we explore the possible impacts of this ecological corridor on genetic diversity within, and connectivity among, populations of a once wide-ranging group, the caballine horses (Equus spp.). Using a panel of 187 mitochondrial and eight nuclear genomes recovered from present-day and extinct caballine horses sampled across the Holarctic, we found that Eurasian horse populations initially diverged from those in North America, their ancestral continent, around 1.0-0.8 million years ago. Subsequent to this split our mitochondrial DNA analysis identified two bidirectional long-range dispersals across the BLB ~875-625 and ~200-50 thousand years ago, during the Middle and Late Pleistocene. Whole genome analysis indicated low levels of gene flow between North American and Eurasian horse populations, which probably occurred as a result of these inferred dispersals. Nonetheless, mitochondrial and nuclear diversity of caballine horse populations retained strong phylogeographical structuring. Our results suggest that barriers to gene flow, currently unidentified but possibly related to habitat distribution across Beringia or ongoing evolutionary divergence, played an important role in shaping the early genetic history of caballine horses, including the ancestors of living horses within Equus ferus.

Predicting sample success for large-scale ancient DNA studies on marine mammals

In recent years, nonhuman ancient DNA studies have begun to focus on larger sample sizes and whole genomes, offering the potential to reveal exciting and hitherto unknown answers to ongoing biological and archaeological questions. However, one major limitation to such studies is the substantial financial and time investments still required during sample screening, due to uncertainty regarding successful sample selection. This study investigates the effect of a wide range of sample properties including latitude, sample age, skeletal element, collagen preservation, and context on endogenous content and DNA damage profiles for 317 ancient and historic pinniped samples collected from across the North Atlantic and surrounding regions. Using generalised linear and mixed-effect models, we found that a range of factors affected DNA preservation within each of the species under consideration. The most important findings were that endogenous content varied significantly within species according to context, the type of skeletal element, the collagen content and collection year. There also appears to be an effect of the sample's geographic origin, with samples from the Arctic generally showing higher endogenous content and lower damage rates. Both latitude and sample age were found to have significant relationships with damage levels, but only for walrus samples. Sex, ontogenetic age and extraction material preparation were not found to have any significant relationship with DNA preservation. Overall, skeletal element and sample context were found to be the most influential factors and should therefore be considered when selecting samples for large-scale ancient genome studies.

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

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

The Legacy of Recurrent Introgression during the Radiation of Hares

Hybridization may often be an important source of adaptive variation, but the extent and long-term impacts of introgression have seldom been evaluated in the phylogenetic context of a radiation. Hares (Lepus) represent a widespread mammalian radiation of 32 extant species characterized by striking ecological adaptations and recurrent admixture. To understand the relevance of introgressive hybridization during the diversification of Lepus, we analyzed whole exome sequences (61.7 Mb) from 15 species of hares (1-4 individuals per species), spanning the global distribution of the genus, and two outgroups. We used a coalescent framework to infer species relationships and divergence times, despite extensive genealogical discordance. We found high levels of allele sharing among species and show that this reflects extensive incomplete lineage sorting and temporally layered hybridization. Our results revealed recurrent introgression at all stages along the Lepus radiation, including recent gene flow between extant species since the last glacial maximum but also pervasive ancient introgression occurring since near the origin of the hare lineages. We show that ancient hybridization between northern hemisphere species has resulted in shared variation of potential adaptive relevance to highly seasonal environments, including genes involved in circadian rhythm regulation, pigmentation, and thermoregulation. Our results illustrate how the genetic legacy of ancestral hybridization may persist across a radiation, leaving a long-lasting signature of shared genetic variation that may contribute to adaptation. [Adaptation; ancient introgression; hybridization; Lepus; phylogenomics.].

African and Asian leopards are highly differentiated at the genomic level

Leopards are the only big cats still widely distributed across the continents of Africa and Asia. They occur in a wide range of habitats and are often found in close proximity to humans. But despite their ubiquity, leopard phylogeography and population history have not yet been studied with genomic tools. Here, we present population-genomic data from 26 modern and historical samples encompassing the vast geographical distribution of this species. We find that Asian leopards are broadly monophyletic with respect to African leopards across almost their entire nuclear genomes. This profound genetic pattern persists despite the animals' high potential mobility, and despite evidence of transfer of African alleles into Middle Eastern and Central Asian leopard populations within the last 100,000 years. Our results further suggest that Asian leopards originated from a single out-of-Africa dispersal event 500-600 thousand years ago and are characterized by higher population structuring, stronger isolation by distance, and lower heterozygosity than African leopards. Taxonomic categories do not take into account the variability in depth of divergence among subspecies. The deep divergence between the African subspecies and Asian populations contrasts with the much shallower divergence among putative Asian subspecies. Reconciling genomic variation and taxonomy is likely to be a growing challenge in the genomics era.

Ocean-wide genomic variation in Gray's beaked whales, Mesoplodon grayi

The deep oceans of the Southern Hemisphere are home to several elusive and poorly studied marine megafauna. In the absence of robust observational data for these species, genetic data can aid inferences on population connectivity, demography and ecology. A previous investigation of genetic diversity and population structure in Gray's beaked whale (Mesoplodon grayi) from Western Australia and New Zealand found high levels of mtDNA diversity, no geographic structure and stable demographic history. To further investigate phylogeographic and demographic patterns across their range, we generated complete mitochondrial and partial nuclear genomes of 16 of the individuals previously analysed and included additional samples from South Africa (n = 2) and South Australia (n = 4), greatly expanding the spatial range of genomic data for the species. Gray's beaked whales are highly elusive and rarely observed, and our data represents a unique and geographically broad dataset. We find relatively high levels of diversity in the mitochondrial genome, despite an absence of population structure at the mitochondrial and nuclear level. Demographic analyses suggest these whales existed at stable levels over at least the past 1.1 million years, with an approximately twofold increase in female effective population size approximately 250 thousand years ago, coinciding with a period of increased Southern Ocean productivity, sea surface temperature and a potential expansion of suitable habitat. Our results suggest that Gray's beaked whales are likely to be resilient to near-future ecosystem changes, facilitating their conservation. Our study demonstrates the utility of low-effort shotgun sequencing in providing ecological information on highly elusive species.

Extensive Genome-Wide Phylogenetic Discordance Is Due to Incomplete Lineage Sorting and Not Ongoing Introgression in a Rapidly Radiated Bryophyte Genus

The relative importance of introgression for diversification has long been a highly disputed topic in speciation research and remains an open question despite the great attention it has received over the past decade. Gene flow leaves traces in the genome similar to those created by incomplete lineage sorting (ILS), and identification and quantification of gene flow in the presence of ILS is challenging and requires knowledge about the true phylogenetic relationship among the species. We use whole nuclear, plastid, and organellar genomes from 12 species in the rapidly radiated, ecologically diverse, actively hybridizing genus of peatmoss (Sphagnum) to reconstruct the species phylogeny and quantify introgression using a suite of phylogenomic methods. We found extensive phylogenetic discordance among nuclear and organellar phylogenies, as well as across the nuclear genome and the nodes in the species tree, best explained by extensive ILS following the rapid radiation of the genus rather than by postspeciation introgression. Our analyses support the idea of ancient introgression among the ancestral lineages followed by ILS, whereas recent gene flow among the species is highly restricted despite widespread interspecific hybridization known in the group. Our results contribute to phylogenomic understanding of how speciation proceeds in rapidly radiated, actively hybridizing species groups, and demonstrate that employing a combination of diverse phylogenomic methods can facilitate untangling complex phylogenetic patterns created by ILS and introgression.

Examining Natural History through the Lens of Palaeogenomics

The many high-resolution tools that are uniquely applicable to specimens from the Quaternary period (the past ~2.5 Ma) provide an opportunity to cross-validate data and test hypotheses based on the morphology and distribution of fossils. Among these tools is palaeogenomics - the genome-scale sequencing of genetic material from ancient specimens - that can provide direct insight into ecology and evolution, potentially improving the accuracy of inferences about past ecological communities over longer timescales. Palaeogenomics has revealed instances of over- and underestimation of extinct diversity, detected cryptic faunal migration and turnover, allowed quantification of widespread sex biases and sexual dimorphism in the fossil record, revealed past hybridisation events and hybrid individuals, and has highlighted previously unrecognised routes of zoonotic disease transfer.

Phylotranscriptomic evidence for pervasive ancient hybridization among Old World salamanders

Hybridization can leave genealogical signatures in an organism's genome, originating from the parental lineages and persisting over time. This potentially confounds phylogenetic inference methods that aim to represent evolution as a strictly bifurcating tree. We apply a phylotranscriptomic approach to study the evolutionary history of, and test for inter-lineage introgression in the Salamandridae, a Holarctic salamanders group of interest in studies of toxicity and aposematism, courtship behavior, and molecular evolution. Although the relationships between the 21 currently recognized salamandrid genera have been the subject of numerous molecular phylogenetic studies, some branches have remained controversial and sometimes affected by discordances between mitochondrial vs. nuclear trees. To resolve the phylogeny of this family, and understand the source of mito-nuclear discordance, we generated new transcriptomic (RNAseq) data for 20 salamandrids and used these along with published data, including 28 mitochondrial genomes, to obtain a comprehensive nuclear and mitochondrial perspective on salamandrid evolution. Our final phylotranscriptomic data set included 5455 gene alignments for 40 species representing 17 of the 21 salamandrid genera. Using concatenation and species-tree phylogenetic methods, we find (1) Salamandrina sister to the clade of the "True Salamanders" (consisting of Chioglossa, Mertensiella, Lyciasalamandra, and Salamandra), (2) Ichthyosaura sister to the Near Eastern genera Neurergus and Ommatotriton, (3) Triturus sister to Lissotriton, and (4) Cynops paraphyletic with respect to Paramesotriton and Pachytriton. Combining introgression tests and phylogenetic networks, we find evidence for introgression among taxa within the clades of "Modern Asian Newts" and "Modern European Newts". However, we could not unambiguously identify the number, position, and direction of introgressive events. Combining evidence from nuclear gene analysis with the observed mito-nuclear phylogenetic discordances, we hypothesize a scenario with hybridization and mitochondrial capture among ancestral lineages of (1) Lissotriton into Ichthyosaura and (2) Triturus into Calotriton, plus introgression of nuclear genes from Triturus into Lissotriton. Furthermore, both mitochondrial capture and nuclear introgression may have occurred among lineages assigned to Cynops. More comprehensive genomic data will, in the future, allow testing this against alternative scenarios involving hybridization with other, extinct lineages of newts.

Dental Calculus as a Tool to Study the Evolution of the Mammalian Oral Microbiome

Dental calculus, the calcified form of the mammalian oral microbial plaque biofilm, is a rich source of oral microbiome, host, and dietary biomolecules and is well preserved in museum and archaeological specimens. Despite its wide presence in mammals, to date, dental calculus has primarily been used to study primate microbiome evolution. We establish dental calculus as a valuable tool for the study of nonhuman host microbiome evolution, by using shotgun metagenomics to characterize the taxonomic and functional composition of the oral microbiome in species as diverse as gorillas, bears, and reindeer. We detect oral pathogens in individuals with evidence of oral disease, assemble near-complete bacterial genomes from historical specimens, characterize antibiotic resistance genes, reconstruct components of the host diet, and recover host genetic profiles. Our work demonstrates that metagenomic analyses of dental calculus can be performed on a diverse range of mammalian species, which will allow the study of oral microbiome and pathogen evolution from a comparative perspective. As dental calculus is readily preserved through time, it can also facilitate the quantification of the impact of anthropogenic changes on wildlife and the environment.

AuthentiCT: a model of ancient DNA damage to estimate the proportion of present-day DNA contamination

Contamination from present-day DNA is a fundamental issue when studying ancient DNA from historical or archaeological material, and quantifying the amount of contamination is essential for downstream analyses. We present AuthentiCT, a command-line tool to estimate the proportion of present-day DNA contamination in ancient DNA datasets generated from single-stranded DNA libraries. The prediction is based solely on the patterns of post-mortem damage observed on ancient DNA sequences. The method has the power to quantify contamination from as few as 10,000 mapped sequences, making it particularly useful for analysing specimens that are poorly preserved or for which little data is available.

Analyses of key genes involved in Arctic adaptation in polar bears suggest selection on both standing variation and de novo mutations played an important role

BACKGROUND

Polar bears are uniquely adapted to an Arctic existence. Since their relatively recent divergence from their closest living relative, brown bears, less than 500,000 years ago, the species has evolved an array of novel traits suited to its Arctic lifestyle. Previous studies sought to uncover the genomic underpinnings of these unique characteristics, and disclosed the genes showing the strongest signal of positive selection in the polar bear lineage. Here, we survey a comprehensive dataset of 109 polar bear and 33 brown bear genomes to investigate the genomic variants within these top genes present in each species. Specifically, we investigate whether fixed homozygous variants in polar bears derived from selection on standing variation in the ancestral gene pool or on de novo mutation in the polar bear lineage.

RESULTS

We find that a large number of sites fixed in polar bears are biallelic in brown bears, suggesting selection on standing variation. Moreover, we uncover sites in which polar bears are fixed for a derived allele while brown bears are fixed for the ancestral allele, which we suggest may be a signal of de novo mutation in the polar bear lineage.

CONCLUSIONS

Our findings suggest that, among other mechanisms, natural selection acting on changes in genes derived from a combination of variation already in the ancestral gene pool, and from de novo missense mutations in the polar bear lineage, may have enabled the rapid adaptation of polar bears to their new Arctic environment.

Recent introgression between Taiga Bean Goose and Tundra Bean Goose results in a largely homogeneous landscape of genetic differentiation

Several studies have uncovered a highly heterogeneous landscape of genetic differentiation across the genomes of closely related species. Specifically, genetic differentiation is often concentrated in particular genomic regions (“islands of differentiation”) that might contain barrier loci contributing to reproductive isolation, whereas the rest of the genome is homogenized by introgression. Alternatively, linked selection can produce differentiation islands in allopatry without introgression. We explored the influence of introgression on the landscape of genetic differentiation in two hybridizing goose taxa: the Taiga Bean Goose (Anser fabalis) and the Tundra Bean Goose (A. serrirostris). We re-sequenced the whole genomes of 18 individuals (9 of each taxon) and, using a combination of population genomic summary statistics and demographic modeling, we reconstructed the evolutionary history of these birds. Next, we quantified the impact of introgression on the build-up and maintenance of genetic differentiation. We found evidence for a scenario of allopatric divergence (about 2.5 million years ago) followed by recent secondary contact (about 60,000 years ago). Subsequent introgression events led to high levels of gene flow, mainly from the Tundra Bean Goose into the Taiga Bean Goose. This scenario resulted in a largely undifferentiated genomic landscape (genome-wide F_ST = 0.033) with a few notable differentiation peaks that were scattered across chromosomes. The summary statistics indicated that some peaks might contain barrier loci while others arose in allopatry through linked selection. Finally, based on the low genetic differentiation, considerable morphological variation and incomplete reproductive isolation, we argue that the Taiga and the Tundra Bean Goose should be treated as subspecies.

Heavy reliance on plants for Romanian cave bears evidenced by amino acid nitrogen isotope analysis

Heavy reliance on plants is rare in Carnivora and mostly limited to relatively small species in subtropical settings. The feeding behaviors of extinct cave bears living during Pleistocene cold periods at middle latitudes have been intensely studied using various approaches including isotopic analyses of fossil collagen. In contrast to cave bears from all other regions in Europe, some individuals from Romania show exceptionally high δ¹⁵N values that might be indicative of meat consumption. Herbivory on plants with high δ¹⁵N values cannot be ruled out based on this method, however. Here we apply an approach using the δ¹⁵N values of individual amino acids from collagen that offsets the baseline δ¹⁵N variation among environments. The analysis yielded strong signals of reliance on plants for Romanian cave bears based on the δ¹⁵N values of glutamate and phenylalanine. These results could suggest that the high variability in bulk collagen δ¹⁵N values observed among cave bears in Romania reflects niche partitioning but in a general trophic context of herbivory.

The contribution of ancient admixture to reproductive isolation between European sea bass lineages

Understanding how new species arise through the progressive establishment of reproductive isolation (RI) barriers between diverging populations is a major goal in Evolutionary Biology. An important result of speciation genomics studies is that genomic regions involved in RI frequently harbor anciently diverged haplotypes that predate the reconstructed history of species divergence. The possible origins of these old alleles remain much debated, as they relate to contrasting mechanisms of speciation that are not yet fully understood. In the European sea bass (Dicentrarchus labrax), the genomic regions involved in RI between Atlantic and Mediterranean lineages are enriched for anciently diverged alleles of unknown origin. Here, we used haplotype-resolved whole-genome sequences to test whether divergent haplotypes could have originated from a closely related species, the spotted sea bass (Dicentrarchus punctatus). We found that an ancient admixture event between D. labrax and D. punctatus is responsible for the presence of shared derived alleles that segregate at low frequencies in both lineages of D. labrax. An exception to this was found within regions involved in RI between the two D. labrax lineages. In those regions, archaic tracts originating from D. punctatus locally reached high frequencies or even fixation in Atlantic genomes but were almost absent in the Mediterranean. We showed that the ancient admixture event most likely occurred between D. punctatus and the D. labrax Atlantic lineage, while Atlantic and Mediterranean D. labrax lineages were experiencing allopatric isolation. Our results suggest that local adaptive introgression and/or the resolution of genomic conflicts provoked by ancient admixture have probably contributed to the establishment of RI between the two D. labrax lineages.

Hyena paleogenomes reveal a complex evolutionary history of cross-continental gene flow between spotted and cave hyena

The genus Crocuta (African spotted and Eurasian cave hyenas) includes several closely related extinct and extant lineages. The relationships among these lineages, however, are contentious. Through the generation of population-level paleogenomes from late Pleistocene Eurasian cave hyena and genomes from modern African spotted hyena, we reveal the cross-continental evolutionary relationships between these enigmatic hyena lineages. We find a deep divergence (~2.5 Ma) between African and Eurasian Crocuta populations, suggesting that ancestral Crocuta left Africa around the same time as early Homo. Moreover, we find discordance between nuclear and mitochondrial phylogenies and evidence for bidirectional gene flow between African and Eurasian Crocuta after the lineages split, which may have complicated prior taxonomic classifications. Last, we find a number of introgressed loci that attained high frequencies within the recipient lineage, suggesting some level of adaptive advantage from admixture.

Biomolecular analyses reveal the age, sex and species identity of a near-intact Pleistocene bird carcass

Ancient remains found in permafrost represent a rare opportunity to study past ecosystems. Here, we present an exceptionally well-preserved ancient bird carcass found in the Siberian permafrost, along with a radiocarbon date and a reconstruction of its complete mitochondrial genome. The carcass was radiocarbon dated to approximately 44–49 ka BP, and was genetically identified as a female horned lark. This is a species that usually inhabits open habitat, such as the steppe environment that existed in Siberia at the time. This near-intact carcass highlights the potential of permafrost remains for evolutionary studies that combine both morphology and ancient nucleic acids.

Nicolas Dussex et al. identify a 44,000–49,000 year old bird found in Siberian permafrost as a female horned lark using ancient DNA. This exceptionally well-preserved specimen illustrates the potential contribution to science of permafrost deposits, such as the study of ecology and evolution of ancient ecosystems, calibration of molecular clocks, and furthering our understanding of processes such as biological regulation and gene expression in relation to climate change.

The problem with petrous? A consideration of the potential biases in the utilization of pars petrosa for ancient DNA analysis

Advances in NGS sequencing technologies, improved laboratory protocols and new bioinformatic workflows have seen huge increases in ancient DNA (aDNA) research on archaeological materials. A large proportion of aDNA work now utilizes the petrous portion of the temporal bone (pars petrosa), which is recognized as an excellent skeletal element for long-term ancient endogenous (host) DNA survival. This has been significant due to the often low endogenous content of other skeletal elements, meaning that large amounts of sequencing are frequently required to obtain sufficient genetic coverage. However, exclusive sampling of the petrous for aDNA analysis introduces a new set of potential biases into our scientific studies - and these issues are yet to be considered by ancient DNA researchers. This paper aims to outline the possible biases of utilizing petrous bones to undertake aDNA analyses and highlight how these complications may potentially be overcome in future research.

Consensify: A Method for Generating Pseudohaploid Genome Sequences from Palaeogenomic Datasets with Reduced Error Rates

A standard practise in palaeogenome analysis is the conversion of mapped short read data into pseudohaploid sequences, frequently by selecting a single high-quality nucleotide at random from the stack of mapped reads. This controls for biases due to differential sequencing coverage, but it does not control for differential rates and types of sequencing error, which are frequently large and variable in datasets obtained from ancient samples. These errors have the potential to distort phylogenetic and population clustering analyses, and to mislead tests of admixture using D statistics. We introduce Consensify, a method for generating pseudohaploid sequences, which controls for biases resulting from differential sequencing coverage while greatly reducing error rates. The error correction is derived directly from the data itself, without the requirement for additional genomic resources or simplifying assumptions such as contemporaneous sampling. For phylogenetic and population clustering analysis, we find that Consensify is less affected by artefacts than methods based on single read sampling. For D statistics, Consensify is more resistant to false positives and appears to be less affected by biases resulting from different laboratory protocols than other frequently used methods. Although Consensify is developed with palaeogenomic data in mind, it is applicable for any low to medium coverage short read datasets. We predict that Consensify will be a useful tool for future studies of palaeogenomes.