Characterization of ancient and modern genomes by SNP detection and phylogenomic and metagenomic analysis using PALEOMIX

Ancient DNA suggests a historical demographic decline and genetic erosion in the Atlantic bluefin tuna

Overexploitation has depleted fish stocks during the past century; nonetheless, its genomic consequences remain poorly understood for most species. Characterizing the spatiotemporal patterns of these consequences may provide baseline estimates of past diversity and productivity to aid management targets, help predict future dynamics, and facilitate the identification of evolutionary factors limiting fish population recovery. Here, we evaluate human impacts on the evolution of the iconic Atlantic bluefin tuna (Thunnus thynnus), one of the longest and most intensely exploited marine fishes, with a tremendous cultural and economic importance. We sequenced whole genomes from modern (n = 49) and ancient (n = 41) specimens dating up to 5,000 y ago, uncovering several findings. First, we identify temporally stable patterns of population admixture, as bluefin tuna caught off Norway and in the eastern Mediterranean share a greater degree of ancestry with Gulf of Mexico bluefin tuna than western and central Mediterranean bluefin tuna. This suggests that Atlantic spawning areas are important mixing grounds for the genetic diversity of Mediterranean bluefin tuna. We model effective population size to show that Mediterranean bluefin tuna began to undergo a demographic decline by the year 1900 to an extent not observed across the previous millennia. Coinciding with this, we found that heterozygosity and nucleotide diversity were significantly lower in modern (2013 to 2020) than ancient (pre-1941) Mediterranean bluefin tuna, suggesting that bluefin tuna underwent a genetic bottleneck. With this work, we show how ancient DNA provides unique perspectives on ecological complexity with the potential to inform the management and conservation of fishes.

Complete mitochondrial genome data and phylogenetic analysis of the Plain Banded Awl, Hasora vitta (Lepidoptera: Hesperiidae: Coeliadinae) from Malaysia

We report the complete mitogenome of Hasora vitta (Plain Banded Awl) sampled from Malaysia. The mitogenome is 15,289 bp long, comprising of 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and a control region. All PCGs were initiated by the typical ATN codon, except for COX1 with a CGA start codon. Two PCGs (COX1 and COX2) were terminated with an incomplete stop codon T. Phylogenetic analysis highly supported the placement of Hasora vitta from Malaysia within Coeliadinae and is clustered with two Hasora vitta samples from China.

Ancient Microbiomes as Mirrored by DNA Extracted From Century-Old Herbarium Plants and Associated Soil

Numerous specimens stored in natural history collections have been involuntarily preserved together with their associated microbiomes. We propose exploiting century-old soils occasionally found on the roots of herbarium plants to assess the diversity of ancient soil microbial communities originally associated with these plants. We extracted total DNA and sequenced libraries produced from rhizospheric soils and roots of four plants preserved in herbaria for more than 120 years in order to characterise the preservation and taxonomic diversity that can be recovered in such contexts. Extracted DNA displayed typical features of ancient DNA, with cytosine deamination at the ends of fragments predominantly shorter than 50 bp. When compared to extant microbiomes, herbarium microbial communities clustered with soil communities and were distinct from communities from other environments. Herbarium communities also displayed biodiversity features and assembly rules typical of soil and plant-associated ones. Soil communities were richer than root-associated ones with which they shared most taxa. Regarding community turnover, we detected collection site, soil versus root and plant species effects. Eukaryotic taxa that displayed a higher abundance in roots were mostly plant pathogens that were not identified among soil-enriched ones. Conservation of these biodiversity features and assembly rules in herbarium-associated microbial communities indicates that herbarium-extracted DNA might reflect the composition of the original plant-associated microbial communities and that preservation in herbaria seemingly did not dramatically alter these characteristics. Using this approach, it should be possible to investigate historical soils and herbarium plant roots to explore the diversity and temporal dynamics of soil microbial communities.

The Evolutionary History of the Extinct Baltic Sea Harp Seal Population

The now-extinct harp seal population that inhabited the Baltic Sea from the Mesolithic to the Iron Age is an enigma. It occurred outside the species' contemporary Arctic range, likely deviated from typical harp seal migratory behaviour, and experienced body size reductions and dramatic population fluctuations leading up to its extinction. Here we use ancient DNA analyses to shed more light on the evolutionary history of the Baltic Sea harp seal population, including its origin, timing of colonisation, diversity and factors contributing to its demise. We generated 49 ancient Baltic and eight ancient Arctic harp seal mitogenomes, which we analysed together with 53 contemporary Arctic harp seal mitogenomes. We detected limited phylogeographic resolution among ancient and contemporary populations, which we interpret as a late Pleistocene range expansion from a common refugial population with subsequent gene flow. Ancient Baltic harp seals were significantly genetically differentiated from contemporary harp seal populations and retained their own genetic composition throughout time. The genetic diversity of Baltic harp seals decreased over time, yet was comparable to that of contemporary populations. This suggests that Baltic harp seals formed a distinct breeding population, which may occasionally have received immigrants from the Arctic but was itself confined in the Baltic Sea until the end. We hypothesise that loss of genetic diversity and the ultimate extinction of the Baltic harp seal population was a consequence of population fluctuations caused by climatic change, reduced salinity and biological productivity, and periodic intense human harvest.

Temporal Loss of Genome-Wide and Immunogenetic Diversity in a Near-Extinct Parrot

Loss of genetic diversity threatens a species' adaptive potential and long-term resilience. Predicted to be extinct by 2038, the orange-bellied parrot (Neophema chrysogaster) is a critically endangered migratory bird threatened by numerous viral, bacterial and fungal diseases. The species has undergone multiple population crashes, reaching a low of three wild-born females and 13 males in 2016, and is now represented by only a single wild population and individuals in the captive breeding program. Here we used our high-quality long-read reference genome, and contemporary (N = 19) and historical (N = 16) resequenced genomes from as early as 1829, to track the long-term genomic erosion and immunogenetic diversity decline in this species. 62% of genomic diversity was lost between historical (mean autosomal heterozygosity = 0.00149 ± 0.000699 SD) and contemporary (0.00057 ± 0.000026) parrots. A greater number and length of runs of homozygosity in contemporary samples were also observed. A temporal reduction in the number of alleles at Toll-like receptor genes was found (historical average alleles = 5.78 ± 2.73; contemporary = 3.89 ± 2.10), potentially exacerbating disease susceptibility in the contemporary population. Of particular concern is the new threat of avian influenza strain (HPAI) to Australia. We discuss the conservation implications of our findings and propose that hybridisation and synthetic biology may be required to address the catastrophic loss of genetic diversity that has occurred in this species in order to prevent extinction.

Genomes of critically endangered saola are shaped by population structure and purging

The saola is one of the most elusive large mammals, standing at the brink of extinction. We constructed a reference genome and resequenced 26 saola individuals, confirming the saola as a basal member of the Bovini. Despite its small geographic range, we found that the saola is partitioned into two populations with high genetic differentiation (FST = 0.49). We estimate that these populations diverged and started declining 5,000-20,000 years ago, possibly due to climate changes and exacerbated by increasing human activities. The saola has long tracts without genomic diversity; however, most of these tracts are not shared by the two populations. Saolas carry a high genetic load, yet their gradual decline resulted in the purging of the most deleterious genetic variation. Finally, we find that combining the two populations, e.g., in an eventual captive breeding program, would mitigate the genetic load and increase the odds of species survival.

Validating a Target-Enrichment Design for Capturing Uniparental Haplotypes in Ancient Domesticated Animals

In the last three decades, DNA sequencing of ancient animal osteological assemblages has become an important tool complementing standard archaeozoological approaches to reconstruct the history of animal domestication. However, osteological assemblages of key archaeological contexts are not always available or do not necessarily preserve enough ancient DNA for a cost-effective genetic analysis. Here, we develop an in-solution target-enrichment approach, based on 80-mer species-specific RNA probes (ranging from 306 to 1686 per species) to characterise (in single experiments) the mitochondrial genetic variation from eight domesticated animal species of major economic interest: cattle, chickens, dogs, donkeys, goats, horses, pigs and sheep. We also illustrate how our design can be adapted to enrich DNA library content and map the Y-chromosomal diversity within Equus caballus. By applying our target-enrichment assay to an extensive panel of ancient osteological remains, farm soil, and cave sediments spanning the last 43 kyrs, we demonstrate that minimal sequencing efforts are necessary to exhaust the DNA library complexity and to characterise mitogenomes to an average depth-of-coverage of 19.4 to 2003.7-fold. Our assay further retrieved horse mitogenome and Y-chromosome data from Late Pleistocene coprolites, as well as bona fide mitochondrial sequences from species that were not part of the probe design, such as bison and cave hyena. Our methodology will prove especially useful to minimise costs related to the genetic analyses of maternal and paternal lineages of a wide range of domesticated and wild animal species, and for mapping their diversity changes over space and time, including from environmental samples.

Road to Extinction? Past and Present Population Structure and Genomic Diversity in the Koala

Koalas are arboreal herbivorous marsupials, endemic to Australia. During the late 1800s and early 1900s, the number of koalas declined dramatically due to hunting for their furs. In addition, anthropogenic activities have further decimated their available habitat, and decreased population numbers. Here, we utilize 37 historic and 25 modern genomes sampled from across their historic and present geographic range, to gain insights into how their population structure and genetic diversity have changed across time; assess the genetic consequences of the period of intense hunting, and the current genetic status of this iconic Australian species. Our analyses reveal how genome-wide heterozygosity has decreased through time and unveil previously uncharacterized mitochondrial haplotypes and nuclear genotypes in the historic dataset, which are absent from today's koala populations.

Five millennia of mitonuclear discordance in Atlantic bluefin tuna identified using ancient DNA

Mitonuclear discordance between species is readily documented in marine fishes. Such discordance may either be the result of past natural phenomena or the result of recent introgression from previously seperated species after shifts in their spatial distributions. Using ancient DNA spanning five millennia, we here investigate the long-term presence of Pacific bluefin tuna (Thunnus orientalis) and albacore (Thunnus alalunga) -like mitochondrial (MT) genomes in Atlantic bluefin tuna (Thunnus thynnus), a species with extensive exploitation history and observed shifts in abundance and age structure. Comparing ancient (n = 130) and modern (n = 78) Atlantic bluefin MT genomes from most of its range, we detect no significant spatial or temporal population structure, which implies ongoing gene flow between populations and large effective population sizes over millennia. Moreover, we identify discordant MT haplotypes in ancient specimens up to 5000 years old and find that the frequency of these haplotypes has remained similar through time. We therefore conclude that MT discordance in the Atlantic bluefin tuna is not driven by recent introgression. Our observations provide oldest example of directly observed MT discordance in the marine environment, highlighting the utility of ancient DNA to obtain insights in the long-term persistence of such phenomena.

Synteny Enabled Upgrade of the Galapagos Giant Tortoise Genome Improves Inferences of Runs of Homozygosity

The utility and importance of whole-genome sequences are recognized across various fields, including evolution and conservation. However, for some taxa, like extinct species, using methods to generate contiguous genomes that rely on high-quality DNA is impossible. In such cases, an alternative may be to employ synteny-based methods using a genome from a closely related taxon to generate more complete genomes. Here we update the reference genome for the Pinta Island Galapagos giant tortoise (Chelonoidis abingdonii) without conducting additional sequencing through rescaffolding against the most closely related chromosome-level genome assembly, the Aldabra giant tortoise (Aldabrachelys gigantea). This effort resulted in a much more contiguous genome, CheloAbing_2.0, with an N50 that is two orders of magnitude longer and large reductions in L50 and the number of gaps. We then examined the impact of the CheloAbing_2.0 genome on estimates of runs of homozygosity (ROH) using genome resequencing data from 37 individual Galapagos giant tortoises from the 13 extant lineages to test the mechanisms by which a fragmented assembly may over- or underestimate the number and extent of ROH. The use of CheloAbing_2.0 resulted in individual estimates of inbreeding, including ROH proportion (FROH), number (NROH), and cumulative length (SROH), that were statistically different from those derived from the earlier genome assembly. This improved genome will serve as a resource for future efforts focusing on the ecology, evolution, and conservation of this species group. More broadly, our results highlight that synteny-based scaffolding is promising for generating contiguous genomes without needing additional data types.

Chromosome-level genomics and historical museum collections reveal new insights into the population structure and chromosome evolution of waterbuck

Advances in the sequencing and assembly of chromosome-level genome assemblies has enabled the study of non-model animals, providing further insights into the evolution of genomes and chromosomes. Here, we present the waterbuck ( Kobus ellipsiprymnus ) as an emerging model antelope for studying population dynamics and chromosome evolution. Antelope evolutionary history has been shaped by Robertsonian (Rb) fusions, with waterbuck also showing variation in karyotype due to two polymorphic Rb fusions. These polymorphisms are variable between and within the two recognised subspecies, the common and defassa waterbuck. To provide new insights into waterbuck evolution, we firstly assembled a chromosome-level genome assembly for the defassa subspecies using PacBio HiFi and Hi-C sequencing. We then utilised museum collections to carry out whole genome sequencing (WGS) of 24 historical waterbuck skins from both subspecies. Combined with a previous WGS dataset (n = 119), this represents the largest study of waterbuck populations to date. We found novel population structure and gene flow between waterbuck populations and regions across the genome with high genomic differentiation between the two subspecies. Several of these regions were found around the centromeres of fixed and polymorphic Rb fusions, exhibiting signatures of low recombination and local population structure. Interestingly, these regions contain genes involved in development, fertility, and recombination. Our results highlight the importance of assembling genomes to the chromosome-level, the utility and value of historical collections in sampling a wide-ranging species to uncover fine-scale population structure, and the potential impacts of Rb fusions on genomic differentiation and the recombination landscape.

Ancient dog introgression into the Iberian wolf genome may have facilitated adaptation to human-dominated landscapes

Understanding how large carnivores respond to increasingly human-dominated landscapes will determine their future adaptive potential. The Iberian wolf (Canis lupus signatus), a gray wolf subspecies endemic to the Iberian Peninsula (Portugal and Spain), has uniquely persisted in human-dominated landscapes, unlike many other wolf populations that faced widespread extinction across Europe during the twentieth century. In this study, we conducted a comprehensive genome-wide analysis of 145 historical and contemporary Iberian wolf samples to investigate whether hybridization with domestic dogs resulted in genetic introgression. We identified a dog-derived block on Chromosome 2 in Iberian wolves, displaying signatures consistent with introgression and high nucleotide similarity among introgressed individuals. Additionally, our estimates place the average timing of introgression between 6100 and 3000 years ago, with low sequence divergence to dogs from the Iberian Peninsula suggesting a single local origin for the hybridization event. Using forward genetic simulations, we show that the introgressed haplotype is most likely being maintained in Iberian wolves by selection. The introgressed dog variants are located within the MAST4 gene, which has been linked to neurological disorders, including cognitive and motor developmental delays, hinting at a potential role in cognitive behavior in Iberian wolves. This study uncovers a case of putative adaptive introgression from domestic dogs into wolves, offering new insights into wild canids' adaptation to human-dominated landscapes.

A Possible More Precise Management Unit Delineation Based on Epigenomic Differentiation of a Long-Distance-Migratory Marine Fish Scomberomorus niphonius

Understanding population structure and adaptive history is critical for designing appropriate management regulations for fisheries and conserving adaptive potential for the future. However, this is not easy for marine fish, especially those with long-distance migration abilities. In this study, we constructed a high-quality reference genome for Japanese Spanish mackerel (Scomberomorus niphonius) and explored its population structure using whole genomic and epigenomic data. Despite the high depth of the sequence data, we failed to identify geographical genetic differentiation of Japanese Spanish mackerel across Chinese coastal waters. However, whole-genome bisulphite sequencing can classify this species into the Bohai-Yellow Sea group and the East China Sea-South China Sea group. Genes involved in embryonic skeletal system development, limb morphogenesis functions, and adult locomotory behaviour were differentially methylated in the southern (Zhanjiang, ZJ) and northern (Western Dalian, WDL) populations and may play important roles as drivers of population structure in Japanese Spanish mackerel. Our study not only provides the first reference genome of the Japanese Spanish mackerel and sheds light on population differentiation at the epigenomic level, but also provides a methylome-based framework for population structure analyses of marine fish with long-distance migration ability. These findings are expected to facilitate the development of scientific programmes for the successful conservation of marine fishery resources.

Copy number variation contributes to parallel local adaptation in an invasive plant

Adaptation is a critical determinant of the diversification, persistence, and geographic range limits of species. Yet the genetic basis of adaptation is often unknown and potentially underpinned by a wide range of mutational types-from single nucleotide changes to large-scale alterations of chromosome structure. Copy number variation (CNV) is thought to be an important source of adaptive genetic variation, as indicated by decades of candidate gene studies that point to CNVs underlying rapid adaptation to strong selective pressures. Nevertheless, population-genomic studies of CNVs face unique logistical challenges not encountered by other forms of genetic variation. Consequently, few studies have systematically investigated the contributions of CNVs to adaptation at a genome-wide scale. We present a genome-wide analysis of CNV contributing to the adaptation of an invasive weed, Ambrosia artemisiifolia. CNVs show clear signatures of parallel local adaptation between North American (native) and European (invaded) ranges, implying widespread reuse of CNVs during adaptation to shared heterogeneous patterns of selection. We used a local principal component analysis (PCA) to genotype CNV regions in whole-genome sequences of samples collected over the last two centuries. We identified 16 large CNV regions of up to 11.85 megabases in length, eight of which show signals of rapid evolutionary change, with pronounced frequency shifts between historic and modern populations. Our results provide compelling genome-wide evidence that CNV underlies rapid adaptation over contemporary timescales of natural populations.

Time-lagged genomic erosion and future environmental risks in a bird on the brink of extinction

Global biodiversity is rapidly declining due to habitat degradation and genomic erosion, highlighting the urgent need to monitor endangered species and their genetic health. Temporal genomics and ecological modelling offer finer resolution than single-time-point measurements, providing a comprehensive view of species' recent and future trajectories. We investigated genomic erosion and environmental suitability in the critically endangered regent honeyeater (Anthochaera phrygia) by sequencing whole genomes of historical and modern specimens and building multi-temporal species distribution models (SDMs) across the last century. The species has declined from hundreds of thousands of individuals to fewer than 300 over the past 100 years. SDMs correctly predicted known patterns of local extinction in southeast Australia. Our demographic reconstructions revealed a gradual population decline from 2000 to 2500 years ago, sharply accelerating in the last 500 years due to climate variability and habitat loss. Despite this substantial demographic collapse, the regent honeyeater has lost only 9% of its genetic diversity, with no evidence of inbreeding or connectivity loss. Also, it exhibits higher diversity than many other threatened bird species. Forward-in-time genomic simulations indicate that this time lag between population decline and genetic diversity loss conceals the risk of ongoing genomic erosion into a future of rapidly degrading environmental suitability. Our work underscores the need for targeted conservation efforts and continuous genetic monitoring to prevent species extinction.

Rapid Parallel Adaptation in Distinct Invasions of Ambrosia Artemisiifolia Is Driven by Large-Effect Structural Variants

When introduced to multiple distinct ranges, invasive species provide a compelling natural experiment for understanding the repeatability of adaptation. Ambrosia artemisiifolia is an invasive, noxious weed, and chief cause of hay fever. Leveraging over 400 whole-genome sequences spanning the native-range in North America and 2 invasions in Europe and Australia, we inferred demographically distinct invasion histories on each continent. Despite substantial differences in genetic source and effective population size changes during introduction, scans of both local climate adaptation and divergence from the native-range revealed genomic signatures of parallel adaptation between invasions. Disproportionately represented among these parallel signatures are 37 large haploblocks-indicators of structural variation-that cover almost 20% of the genome and exist as standing genetic variation in the native-range. Many of these haploblocks are associated with traits important for adaptation to local climate, like size and the timing of flowering, and have rapidly reformed native-range clines in invaded ranges. Others show extreme frequency divergence between ranges, consistent with a response to divergent selection on different continents. Our results demonstrate the key role of large-effect standing variants in rapid adaptation during range expansion, a pattern that is robust to diverse invasion histories.

Beyond population size: Whole-genome data reveal bottleneck legacies in the peninsular Italian wolf

Preserving genetic diversity and adaptive potential while avoiding inbreeding depression is crucial for the long-term conservation of natural populations. Despite demographic increases, traces of past bottleneck events at the genomic level should be carefully considered for population management. From this perspective, the peninsular Italian wolf is a paradigmatic case. After being on the brink of extinction in the late 1960s, peninsular Italian wolves rebounded and recolonized most of the peninsula aided by conservation measures, including habitat and legal protection. Notwithstanding their demographic recovery, a comprehensive understanding of the genomic consequences of the historical bottleneck in Italian wolves is still lacking. To fill this gap, we sequenced whole genomes of 13 individuals sampled in the core historical range of the species in Central Italy to conduct population genomic analyses, including a comparison with wolves from two highly-inbred wolf populations (i.e. Scandinavia and Isle Royale). We found that peninsular Italian wolves, despite their recent recovery, still exhibit relatively low genetic diversity, a small effective population size, signatures of inbreeding, and a non-negligible genetic load. Our findings indicate that the peninsular Italian wolf population is still susceptible to bottleneck legacies, which could lead to local inbreeding depression in case of population reduction or fragmentations. This study emphasizes the importance of considering key genetic parameters to design appropriate long-term conservation management plans.

The Once and Future Fish: Assessing a Millennium of Atlantic Herring Exploitation Through Mixed-Stock Analysis and Ancient DNA

Small pelagic fish support profitable fisheries and are important for food security around the world. Yet, their sustainable management can be hindered by the indiscriminate impacts of simultaneous exploitation of fish from multiple distinct biological populations over extended periods of time. The quantification of such impacts is greatly facilitated by recently developed molecular tools-including diagnostic single nucleotide polymorphism (SNP) panels for mixed-stock analysis (MSA)-that can accurately detect the population identity of individual fish. However, the biological relevance of such tools over longer periods of time remains unknown. Here, we demonstrate that diagnostic SNP panels designed for contemporary MSA in Atlantic herring have a millennium-long biological relevance and applicability. We assign the population identity of ancient Atlantic herring specimens-obtained through famously profitable historic fisheries-up to 1300 years old from eight archaeological sites across Europe. Analyzing contemporary and ancient whole-genome data, we obtain evidence for the long-term mixed-stock exploitation of Atlantic herring. Despite such mixed-stock exploitation, we exclusively identify autumn-spawning herring amongst these archaeological remains, indicative of a specific biological availability or cultural preference for certain herring ecotypes in the past. Moreover, our results show that herring demographic patterns were relatively stable until the dramatic disruptions and stock collapses during the 20th century. We find small but significant reductions in genetic diversity over time, indicating long-term evolutionary consequences from 20th-century stock declines. The long-term applicability of diagnostic SNP panels underscores their biological relevance and cost-effective application for the genetic monitoring of herring stocks and highlights the utility of ancient DNA to obtain insights in herring ecology and population dynamics.

Persistent Gene Flow Suggests an Absence of Reproductive Isolation in an African Antelope Speciation Model

African antelope diversity is a globally unique vestige of a much richer world-wide Pleistocene megafauna. Despite this, the evolutionary processes leading to the prolific radiation of African antelopes are not well understood. Here, we sequenced 145 whole genomes from both subspecies of the waterbuck (Kobus ellipsiprymnus), an African antelope believed to be in the process of speciation. We investigated genetic structure and population divergence and found evidence of a mid-Pleistocene separation on either side of the eastern Great Rift Valley, consistent with vicariance caused by a rain shadow along the so-called "Kingdon's Line." However, we also found pervasive evidence of both recent and widespread historical gene flow across the Rift Valley barrier. By inferring the genome-wide landscape of variation among subspecies, we found 14 genomic regions of elevated differentiation, including a locus that may be related to each subspecies' distinctive coat pigmentation pattern. We investigated these regions as candidate speciation islands. However, we observed no significant reduction in gene flow in these regions, nor any indications of selection against hybrids. Altogether, these results suggest a pattern whereby climatically driven vicariance is the most important process driving the African antelope radiation and suggest that reproductive isolation may not set in until very late in the divergence process. This has a significant impact on taxonomic inference, as many taxa will be in a gray area of ambiguous systematic status, possibly explaining why it has been hard to achieve consensus regarding the species status of many African antelopes. Our analyses demonstrate how population genetics based on low-depth whole genome sequencing can provide new insights that can help resolve how far lineages have gone along the path to speciation.

Human DNA from the oldest Eneolithic cemetery in Nalchik points the spread of farming from the Caucasus to the Eastern European steppes

The Darkveti-Meshoko culture (c.5000-3500/3300 BCE) is the earliest known farming community in the Northern Caucasus, but its contribution to the genetic profile of the neighboring steppe herders has remained unclear. We present analysis of human DNA from the Nalchik cemetery-the oldest Eneolithic site in the Northern Caucasus-which shows a link with the LowerVolga's first herders of the Khvalynsk culture. The Nalchik male genotype combines the genes of the Caucasus hunter-gatherers, the Eastern hunter-gatherers and the Pre-Pottery Neolithic (PPN) farmers of western Asia. Improved comparative analysis suggests that the genetic profile of certain Khvalynsk individuals shares the genetic ancestry of the Unakozovo-Nalchik type population of the Northern Caucasus' Eneolithic. Therefore, it seems that in the first half of the 5th millennium BCE, cultural and mating networks helped agriculture and pastoralism spread from West Asia across the Caucasian, into the steppes between the Don and the Volga in Eastern Europe.

Tracing 600 years of long-distance Atlantic cod trade in medieval and post-medieval Oslo using stable isotopes and ancient DNA

Marine resources have been important for the survival and economic development of coastal human communities across northern Europe for millennia. Knowledge of the origin of such historic resources can provide key insights into fishing practices and the spatial extent of trade networks. Here, we combine ancient DNA and stable isotopes (δ13C, δ15N, non-exchangeable δ2H and δ34S) to investigate the geographical origin of archaeological cod remains in Oslo from the eleventh to seventeenth centuries CE. Our findings provide genetic evidence that Atlantic cod was obtained from different geographical populations, including a variety of distant-water populations like northern Norway and possibly Iceland. Evidence for such long-distance cod trade is already observed from the eleventh century, contrasting with archaeological and historical evidence from Britain and other areas of Continental Europe around the North and Baltic Seas, where such trade increased during the thirteenth to fourteenth centuries. The genomic assignments of specimens to different populations coincide with significantly different δ13C values between those same specimens, indicating that multiple Atlantic cod populations living in different environments were exploited. This research provides novel information about the exploitation timeline of specific Atlantic cod stocks and highlights the utility of combining ancient DNA (aDNA) methods and stable isotope analysis to describe the development of medieval and post-medieval marine fisheries.

Genome-wide population affinities and signatures of adaptation in hydruntines, sussemiones and Asian wild asses

The extremely rich palaeontological record of the horse family, also known as equids, has provided many examples of macroevolutionary change over the last ~55 Mya. This family is also one of the most documented at the palaeogenomic level, with hundreds of ancient genomes sequenced. While these data have advanced understanding of the domestication history of horses and donkeys, the palaeogenomic record of other equids remains limited. In this study, we have generated genome-wide data for 25 ancient equid specimens spanning over 44 Ky and spread across Anatolia, the Caucasus, Central Asia and Mongolia. Our dataset includes the genomes from two extinct species, the European wild ass, Equus hydruntinus, and the sussemione Equus ovodovi. We document, for the first time, the presence of sussemiones in Mongolia and their survival around ~3.9 Kya, a finding that should be considered when discussing the timing of the first arrival of the domestic horse in the region. We also identify strong spatial differentiation within the historical ecological range of Asian wild asses, Equus hemionus, and incomplete reproductive isolation in several groups yet considered as different species. Finally, we find common selection signatures at ANTXR2 gene in European, Asian and African wild asses. This locus, which encodes a receptor for bacterial toxins, shows no selection signal in E. ovodovi, but a 5.4-kb deletion within intron 7. Whether such genetic modifications played any role in the sussemione extinction remains unknown.

Greenland Norse walrus exploitation deep into the Arctic

Walrus ivory was a prized commodity in medieval Europe and was supplied by Norse intermediaries who expanded across the North Atlantic, establishing settlements in Iceland and Greenland. However, the precise sources of the traded ivory have long remained unclear, raising important questions about the sustainability of commercial walrus harvesting, the extent to which Greenland Norse were able to continue mounting their own long-range hunting expeditions, and the degree to which they relied on trading ivory with the various Arctic Indigenous peoples that they were starting to encounter. We use high-resolution genomic sourcing methods to track walrus artifacts back to specific hunting grounds, demonstrating that Greenland Norse obtained ivory from High Arctic waters, especially the North Water Polynya, and possibly from the interior Canadian Arctic. These results substantially expand the assumed range of Greenland Norse ivory harvesting activities and support intriguing archaeological evidence for substantive interactions with Thule Inuit, plus possible encounters with Tuniit (Late Dorset Pre-Inuit).

Late Pleistocene polar bear genomes reveal the timing of allele fixation in key genes associated with Arctic adaptation

The polar bear (Ursus maritimus) occupies a relatively narrow ecological niche, with many traits adapted for cold temperatures, movement across snow, ice and open water, and for consuming highly lipid-dense prey species. The divergence of polar bears from brown bears (Ursus arctos) and their adaptation to their Arctic lifestyle is a well-known example of rapid evolution. Previous research investigating whole genomes uncovered twelve key genes that are highly differentiated between polar and brown bears, show signatures of selection in the polar bear lineage, and are associated with polar bear adaptation to the Arctic environment. Further research suggested fixed derived alleles in these genes arose from selection on both standing variation and de novo mutations in the evolution of polar bears. Here, we reevaluate these findings based on a larger and geographically more representative dataset of 119 polar bears and 135 brown bears, and assess the timing of derived allele fixation in polar bears by incorporating the genomes of two Late Pleistocene individuals (aged 130-100,000 years old and 100-70,000 years old). In contrast with previous results, we found no evidence of derived alleles fixed in present-day polar bears within the key genes arising from de novo mutation. Most derived alleles fixed in present-day polar bears were also fixed in the Late Pleistocene polar bears, suggesting selection occurred prior to 70,000 years ago. However, some derived alleles fixed in present-day polar bears were not fixed in the two Late Pleistocene polar bears, including at sites within APOB, LYST, and TTN. These three genes are associated with cardiovascular function, metabolism, and pigmentation, suggesting selection may have acted on different loci at different times.

Long genetic and social isolation in Neanderthals before their extinction

Neanderthal genomes have been recovered from sites across Eurasia, painting an increasingly complex picture of their populations' structure that mostly indicates that late European Neanderthals belonged to a single metapopulation with no significant evidence of population structure. Here, we report the discovery of a late Neanderthal individual, nicknamed "Thorin," from Grotte Mandrin in Mediterranean France, and his genome. These dentognathic fossils, including a rare example of distomolars, are associated with a rich archeological record of Neanderthal final technological traditions in this region ∼50-42 thousand years ago. Thorin's genome reveals a relatively early divergence of ∼105 ka with other late Neanderthals. Thorin belonged to a population with a small group size that showed no genetic introgression with other known late European Neanderthals, revealing some 50 ka of genetic isolation of his lineage despite them living in neighboring regions. These results have important implications for resolving competing hypotheses about causes of the disappearance of the Neanderthals.

Elucidating the sustainability of 700 y of Inuvialuit beluga whale hunting in the Mackenzie River Delta, Northwest Territories, Canada

Beluga whales play a critical role in the subsistence economies and cultural heritage of Indigenous communities across the Arctic, yet the effects of Indigenous hunting on beluga whales remain unknown. Here, we integrate paleogenomics, genetic simulations, and stable δ13C and δ15N isotope analysis to investigate 700 y of beluga subsistence hunting in the Mackenzie Delta area of northwestern Canada. Genetic identification of the zooarchaeological remains, which is based on radiocarbon dating, span three time periods (1290 to 1440 CE; 1450 to 1650 CE; 1800 to 1870 CE), indicates shifts across time in the sex ratio of the harvested belugas. The equal number of females and males harvested in 1450 to 1650 CE versus more males harvested in the two other time periods may reflect changes in hunting practices or temporal shifts in beluga availability. We find temporal shifts and sex-based differences in δ13C of the harvested belugas across time, suggesting historical adaptability in the foraging ecology of the whales. We uncovered distinct mitochondrial diversity unique to the Mackenzie Delta belugas, but found no changes in nuclear genomic diversity nor any substructuring across time. Our findings indicate the genomic stability and continuity of the Mackenzie Delta beluga population across the 700 y surveyed, indicating the impact of Inuvialuit subsistence harvests on the genetic diversity of contemporary beluga individuals has been negligible.

Morphometrics and Phylogenomics of Coca (Erythroxylum spp.) Illuminate Its Reticulate Evolution, With Implications for Taxonomy

South American coca (Erythroxylum coca and E. novogranatense) has been a keystone crop for many Andean and Amazonian communities for at least 8,000 years. However, over the last half-century, global demand for its alkaloid cocaine has driven intensive agriculture of this plant and placed it in the center of armed conflict and deforestation. To monitor the changing landscape of coca plantations, the United Nations Office on Drugs and Crime collects annual data on their areas of cultivation. However, attempts to delineate areas in which different varieties are grown have failed due to limitations around identification. In the absence of flowers, identification relies on leaf morphology, yet the extent to which this is reflected in taxonomy is uncertain. Here, we analyze the consistency of the current naming system of coca and its four closest wild relatives (the "coca clade"), using morphometrics, phylogenomics, molecular clocks, and population genomics. We include name-bearing type specimens of coca's closest wild relatives E. gracilipes and E. cataractarum. Morphometrics of 342 digitized herbarium specimens show that leaf shape and size fail to reliably discriminate between species and varieties. However, the statistical analyses illuminate that rounder and more obovate leaves of certain varieties could be associated with the subtle domestication syndrome of coca. Our phylogenomic data indicate extensive gene flow involving E. gracilipes which, combined with morphometrics, supports E. gracilipes being retained as a single species. Establishing a robust evolutionary-taxonomic framework for the coca clade will facilitate the development of cost-effective genotyping methods to support reliable identification.

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.

When birds of a feather flock together: Severe genomic erosion and the implications for genetic rescue in an endangered island passerine

The Seychelles magpie-robin's (SMR) five island populations exhibit some of the lowest recorded levels of genetic diversity among endangered birds, and high levels of inbreeding. These populations collapsed during the 20th century, and the species was listed as Critically Endangered in the IUCN Red List in 1994. An assisted translocation-for-recovery program initiated in the 1990s increased the number of mature individuals, resulting in its downlisting to Endangered in 2005. Here, we explore the temporal genomic erosion of the SMR based on a dataset of 201 re-sequenced whole genomes that span the past ~150 years. Our sample set includes individuals that predate the bottleneck by up to 100 years, as well as individuals from contemporary populations established during the species recovery program. Despite the SMR's recent demographic recovery, our data reveal a marked increase in both the genetic load and realized load in the extant populations when compared to the historical samples. Conservation management may have reduced the intensity of selection by increasing juvenile survival and relaxing intraspecific competition between individuals, resulting in the accumulation of loss-of-function mutations (i.e. severely deleterious variants) in the rapidly recovering population. In addition, we found a 3-fold decrease in genetic diversity between temporal samples. While the low genetic diversity in modern populations may limit the species' adaptability to future environmental changes, future conservation efforts (including IUCN assessments) may also need to assess the threats posed by their high genetic load. Our computer simulations highlight the value of translocations for genetic rescue and show how this could halt genomic erosion in threatened species such as the SMR.

Ancient Mitochondrial Genomes Provide New Clues in the History of the Akhal-Teke Horse in China

This study analyzed ancient DNA from the remains of horses unearthed from the Shihuyao tombs. These were found to date from the Han and Tang Dynasties in Xinjiang (approximately 2200 to 1100 years ago). Two high-quality mitochondrial genomes were acquired and analyzed using next-generation sequencing. The genomes were split into two maternal haplogroups, B and D, according to a study that included ancient and contemporary samples from Eurasia. A close genetic affinity was observed between the horse of the Tang Dynasty and Akhal-Teke horses according to the primitive horse haplotype G1. Historical evidence suggests that the ancient Silk Road had a vital role in their dissemination. Additionally, the matrilineal history of the Akhal-Teke horse was accessed and suggested that the early domestication of the breed was for military purposes.

Minimally destructive hDNA extraction method for retrospective genetics of pinned historical Lepidoptera specimens

The millions of specimens stored in entomological collections provide a unique opportunity to study historical insect diversity. Current technologies allow to sequence entire genomes of historical specimens and estimate past genetic diversity of present-day endangered species, advancing our understanding of anthropogenic impact on genetic diversity and enabling the implementation of conservation strategies. A limiting challenge is the extraction of historical DNA (hDNA) of adequate quality for sequencing platforms. We tested four hDNA extraction protocols on five body parts of pinned false heath fritillary butterflies, Melitaea diamina, aiming to minimise specimen damage, preserve their scientific value to the collections, and maximise DNA quality and yield for whole-genome re-sequencing. We developed a very effective approach that successfully recovers hDNA appropriate for short-read sequencing from a single leg of pinned specimens using silica-based DNA extraction columns and an extraction buffer that includes SDS, Tris, Proteinase K, EDTA, NaCl, PTB, and DTT. We observed substantial variation in the ratio of nuclear to mitochondrial DNA in extractions from different tissues, indicating that optimal tissue choice depends on project aims and anticipated downstream analyses. We found that sufficient DNA for whole genome re-sequencing can reliably be extracted from a single leg, opening the possibility to monitor changes in genetic diversity maintaining the scientific value of specimens while supporting current and future conservation strategies.

An endothelial regulatory module links blood pressure regulation with elite athletic performance

The control of transcription is crucial for homeostasis in mammals. A previous selective sweep analysis of horse racing performance revealed a 19.6 kb candidate regulatory region 50 kb downstream of the Endothelin3 (EDN3) gene. Here, the region was narrowed to a 5.5 kb span of 14 SNVs, with elite and sub-elite haplotypes analyzed for association to racing performance, blood pressure and plasma levels of EDN3 in Coldblooded trotters and Standardbreds. Comparative analysis of human HiCap data identified the span as an enhancer cluster active in endothelial cells, interacting with genes relevant to blood pressure regulation. Coldblooded trotters with the sub-elite haplotype had significantly higher blood pressure compared to horses with the elite performing haplotype during exercise. Alleles within the elite haplotype were part of the standing variation in pre-domestication horses, and have risen in frequency during the era of breed development and selection. These results advance our understanding of the molecular genetics of athletic performance and vascular traits in both horses and humans.

Facilitating accessible, rapid, and appropriate processing of ancient metagenomic data with AMDirT

Background

Access to sample-level metadata is important when selecting public metagenomic sequencing datasets for reuse in new biological analyses. The Standards, Precautions, and Advances in Ancient Metagenomics community (SPAAM, https://spaam-community.org) has previously published AncientMetagenomeDir, a collection of curated and standardised sample metadata tables for metagenomic and microbial genome datasets generated from ancient samples. However, while sample-level information is useful for identifying relevant samples for inclusion in new projects, Next Generation Sequencing (NGS) library construction and sequencing metadata are also essential for appropriately reprocessing ancient metagenomic data. Currently, recovering information for downloading and preparing such data is difficult when laboratory and bioinformatic metadata is heterogeneously recorded in prose-based publications.

Methods

Through a series of community-based hackathon events, AncientMetagenomeDir was updated to provide standardised library-level metadata of existing and new ancient metagenomic samples. In tandem, the companion tool 'AMDirT' was developed to facilitate rapid data filtering and downloading of ancient metagenomic data, as well as improving automated metadata curation and validation for AncientMetagenomeDir.

Results

AncientMetagenomeDir was extended to include standardised metadata of over 6000 ancient metagenomic libraries. The companion tool 'AMDirT' provides both graphical- and command-line interface based access to such metadata for users from a wide range of computational backgrounds. We also report on errors with metadata reporting that appear to commonly occur during data upload and provide suggestions on how to improve the quality of data sharing by the community.

Conclusions

Together, both standardised metadata reporting and tooling will help towards easier incorporation and reuse of public ancient metagenomic datasets into future analyses.

Advancing responsible genomic analyses of ancient mollusc shells

The analysis of the DNA entrapped in ancient shells of molluscs has the potential to shed light on the evolution and ecology of this very diverse phylum. Ancient genomics could help reconstruct the responses of molluscs to past climate change, pollution, and human subsistence practices at unprecedented temporal resolutions. Applications are however still in their infancy, partly due to our limited knowledge of DNA preservation in calcium carbonate shells and the need for optimized methods for responsible genomic data generation. To improve ancient shell genomic analyses, we applied high-throughput DNA sequencing to 27 Mytilus mussel shells dated to ~111-6500 years Before Present, and investigated the impact, on DNA recovery, of shell imaging, DNA extraction protocols and shell sub-sampling strategies. First, we detected no quantitative or qualitative deleterious effect of micro-computed tomography for recording shell 3D morphological information prior to sub-sampling. Then, we showed that double-digestion and bleach treatment of shell powder prior to silica-based DNA extraction improves shell DNA recovery, also suggesting that DNA is protected in preservation niches within ancient shells. Finally, all layers that compose Mytilus shells, i.e., the nacreous (aragonite) and prismatic (calcite) carbonate layers, with or without the outer organic layer (periostracum) proved to be valuable DNA reservoirs, with aragonite appearing as the best substrate for genomic analyses. Our work contributes to the understanding of long-term molecular preservation in biominerals and we anticipate that resulting recommendations will be helpful for future efficient and responsible genomic analyses of ancient mollusc shells.

Colonial-driven extinction of the blue antelope despite genomic adaptation to low population size

Low genomic diversity is generally indicative of small population size and is considered detrimental by decreasing long-term adaptability.1,2,3,4,5,6 Moreover, small population size may promote gene flow with congeners and outbreeding depression.7,8,9,10,11,12,13 Here, we examine the connection between habitat availability, effective population size (Ne), and extinction by generating a 40× nuclear genome from the extinct blue antelope (Hippotragus leucophaeus). Historically endemic to the relatively small Cape Floristic Region in southernmost Africa,14,15 populations were thought to have expanded and contracted across glacial-interglacial cycles, tracking suitable habitat.16,17,18 However, we found long-term low Ne, unaffected by glacial cycles, suggesting persistence with low genomic diversity for many millennia prior to extinction in ∼AD 1800. A lack of inbreeding, alongside high levels of genetic purging, suggests adaptation to this long-term low Ne and that human impacts during the colonial era (e.g., hunting and landscape transformation), rather than longer-term ecological processes, were central to its extinction. Phylogenomic analyses uncovered gene flow between roan (H. equinus) and blue antelope, as well as between roan and sable antelope (H. niger), approximately at the time of divergence of blue and sable antelope (∼1.9 Ma). Finally, we identified the LYST and ASIP genes as candidates for the eponymous bluish pelt color of the blue antelope. Our results revise numerous aspects of our understanding of the interplay between genomic diversity and evolutionary history and provide the resources for uncovering the genetic basis of this extinct species' unique traits.

Introgression and disruption of migration routes have shaped the genetic integrity of wildebeest populations

The blue wildebeest (Connochaetes taurinus) is a keystone species in savanna ecosystems from southern to eastern Africa, and is well known for its spectacular migrations and locally extreme abundance. In contrast, the black wildebeest (C. gnou) is endemic to southern Africa, barely escaped extinction in the 1900s and is feared to be in danger of genetic swamping from the blue wildebeest. Despite the ecological importance of the wildebeest, there is a lack of understanding of how its unique migratory ecology has affected its gene flow, genetic structure and phylogeography. Here, we analyze whole genomes from 121 blue and 22 black wildebeest across the genus' range. We find discrete genetic structure consistent with the morphologically defined subspecies. Unexpectedly, our analyses reveal no signs of recent interspecific admixture, but rather a late Pleistocene introgression of black wildebeest into the southern blue wildebeest populations. Finally, we find that migratory blue wildebeest populations exhibit a combination of long-range panmixia, higher genetic diversity and lower inbreeding levels compared to neighboring populations whose migration has recently been disrupted. These findings provide crucial insights into the evolutionary history of the wildebeest, and tangible genetic evidence for the negative effects of anthropogenic activities on highly migratory ungulates.

Giraffe lineages are shaped by major ancient admixture events

Strong genetic structure has prompted discussion regarding giraffe taxonomy,1,2,3 including a suggestion to split the giraffe into four species: Northern (Giraffa c. camelopardalis), Reticulated (G. c. reticulata), Masai (G. c. tippelskirchi), and Southern giraffes (G. c. giraffa).4,5,6 However, their evolutionary history is not yet fully resolved, as previous studies used a simple bifurcating model and did not explore the presence or extent of gene flow between lineages. We therefore inferred a model that incorporates various evolutionary processes to assess the drivers of contemporary giraffe diversity. We analyzed whole-genome sequencing data from 90 wild giraffes from 29 localities across their current distribution. The most basal divergence was dated to 280 kya. Genetic differentiation, FST, among major lineages ranged between 0.28 and 0.62, and we found significant levels of ancient gene flow between them. In particular, several analyses suggested that the Reticulated lineage evolved through admixture, with almost equal contribution from the Northern lineage and an ancestral lineage related to Masai and Southern giraffes. These new results highlight a scenario of strong differentiation despite gene flow, providing further context for the interpretation of giraffe diversity and the process of speciation in general. They also illustrate that conservation measures need to target various lineages and sublineages and that separate management strategies are needed to conserve giraffe diversity effectively. Given local extinctions and recent dramatic declines in many giraffe populations, this improved understanding of giraffe evolutionary history is relevant for conservation interventions, including reintroductions and reinforcements of existing populations.

The complete mitochondrial genome of Meller's mongoose (Rhynchogale melleri)

Meller's mongoose (Rhynchogale melleri) is a member of the family Herpestidae (Mammalia: Carnivora) and the sole species in the genus Rhynchogale. It is primarily found in savannas and open woodlands of eastern sub-Saharan Africa. Here, we report the first complete mitochondrial genome for a female Meller's mongoose collected in Tanzania, generated using a genome-skimming approach. The mitogenome had a final length of 16,644 bp and a total of 37 annotated genes. Phylogenetic analysis validated the placement of this species in the herpestid subfamily Herpestinae. Ultimately, the outcomes of this research offer a genetic foundation for future studies of Meller's mongoose.

High-throughput sequencing of insect specimens with sub-optimal DNA preservation using a practical, plate-based Illumina-compatible Tn5 transposase library preparation method

Entomological sampling and storage conditions often prioritise efficiency, practicality and conservation of morphological characteristics, and may therefore be suboptimal for DNA preservation. This practice can impact downstream molecular applications, such as the generation of high-throughput genomic libraries, which often requires substantial DNA input amounts. Here, we use a practical Tn5 transposase tagmentation-based library preparation method optimised for 96-well plates and low yield DNA extracts from insect legs that were stored under sub-optimal conditions for DNA preservation. The samples were kept in field vehicles for extended periods of time, before long-term storage in ethanol in the freezer, or dry at room temperature. By reducing DNA input to 6ng, more samples with sub-optimal DNA yields could be processed. We matched this low DNA input with a 6-fold dilution of a commercially available tagmentation enzyme, significantly reducing library preparation costs. Costs and workload were further suppressed by direct post-amplification pooling of individual libraries. We generated medium coverage (>3-fold) genomes for 88 out of 90 specimens, with an average of approximately 10-fold coverage. While samples stored in ethanol yielded significantly less DNA compared to those which were stored dry, these samples had superior sequencing statistics, with longer sequencing reads and higher rates of endogenous DNA. Furthermore, we find that the efficiency of tagmentation-based library preparation can be improved by a thorough post-amplification bead clean-up which selects against both short and large DNA fragments. By opening opportunities for the use of sub-optimally preserved, low yield DNA extracts, we broaden the scope of whole genome studies of insect specimens. We therefore expect these results and this protocol to be valuable for a range of applications in the field of entomology.

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

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

Retrieving complete plastid genomes of endangered Guibourtia timber using hybridization capture for forensic identification and phylogenetic analysis

The high economic value and increased demand for timber have led to illegal logging and overexploitation, threatening wild populations. In this context, there is an urgent need to develop effective and accurate forensic tools for identifying endangered Guibourtia timber species to protect forest ecosystem resources and regulate their trade. In this study, a hybridization capture method was developed and applied to explore the feasibility of retrieving complete plastid genomes from Guibourtia sapwood and heartwood specimens stored in a xylarium (wood collection). We then carried out forensic identification and phylogenetic analyses of Guibourtia within the subfamily Detarioideae. This study is the first to successfully retrieve high-quality plastid genomes from xylarium specimens, with 76.95-99.97% coverage. The enrichment efficiency, sequence depth, and coverage of plastid genomes from sapwood were 16.73 times, 70.47 times and 1.14 times higher, respectively, than those from heartwood. Although the DNA capture efficiency of heartwood was lower than that of sapwood, the hybridization capture method used in this study is still suitable for heartwood DNA analysis. Based on the complete plastid genome, we identified six endangered or commonly traded Guibourtia woods at the species level. This technique also has the potential for geographical traceability, especially for Guibourtia demeusei and Guibourtia ehie. Meanwhile, Bayesian phylogenetic analysis suggested that these six Guibourtia species diverged from closely related species within the subfamily Detarioideae ca. 18 Ma during the Miocene. The DNA reference database established based on the xylarium specimens provides admissible evidence for diversity conservation and evolutionary analyses of endangered Guibourtia species.

Ancient mitochondrial genomes unveil the origins and evolutionary history of New Zealand's enigmatic takahē and moho

Many avian species endemic to Aotearoa New Zealand were driven to extinction or reduced to relict populations following successive waves of human arrival, due to hunting, habitat destruction and the introduction of mammalian predators. Among the affected species were the large flightless South Island takahē (Porphyrio hochstetteri) and the moho (North Island takahē; P. mantelli), with the latter rendered extinct and the former reduced to a single relictual population. Little is known about the evolutionary history of these species prior to their decline and/or extinction. Here we sequenced mitochondrial genomes from takahē and moho subfossils (12 takahē and 4 moho) and retrieved comparable sequence data from takahē museum skins (n = 5) and contemporary individuals (n = 17) to examine the phylogeny and recent evolutionary history of these species. Our analyses suggest that prehistoric takahē populations lacked deep phylogeographic structure, in contrast to moho, which exhibited significant spatial genetic structure, albeit based on limited sample sizes (n = 4). Temporal genetic comparisons show that takahē have lost much of their mitochondrial genetic diversity, likely due to a sudden demographic decline soon after human arrival (~750 years ago). Time-calibrated phylogenetic analyses strongly support a sister species relationship between takahē and moho, suggesting these flightless taxa diverged around 1.5 million years ago, following a single colonisation of New Zealand by a flighted Porphyrio ancestor approximately 4 million years ago. This study highlights the utility of palaeogenetic approaches for informing the conservation and systematic understanding of endangered species whose ranges have been severely restricted by anthropogenic impacts.

2-Hydroxyglutarate modulates histone methylation at specific loci and alters gene expression via Rph1 inhibition

2-Hydroxyglutarate (2-HG) is an oncometabolite that accumulates in certain cancers. Gain-of-function mutations in isocitrate dehydrogenase lead to 2-HG accumulation at the expense of alpha-ketoglutarate. Elevated 2-HG levels inhibit histone and DNA demethylases, causing chromatin structure and gene regulation changes with tumorigenic consequences. We investigated the effects of elevated 2-HG levels in Saccharomyces cerevisiae, a yeast devoid of DNA methylation and heterochromatin-associated histone methylation. Our results demonstrate genetic background-dependent gene expression changes and altered H3K4 and H3K36 methylation at specific loci. Analysis of histone demethylase deletion strains indicated that 2-HG inhibits Rph1 sufficiently to induce extensive gene expression changes. Rph1 is the yeast homolog of human KDM4 demethylases and, among the yeast histone demethylases, was the most sensitive to the inhibitory effect of 2-HG in vitro. Interestingly, Rph1 deficiency favors gene repression and leads to further down-regulation of already silenced genes marked by low H3K4 and H3K36 trimethylation, but abundant in H3K36 dimethylation. Our results provide novel insights into the genome-wide effects of 2-HG and highlight Rph1 as its preferential demethylase target.

Tracking population structure and phenology through time using ancient genomes from waterlogged white oak wood

Whole genome characterizations of crop plants based on ancient DNA have provided unique keys for a better understanding of the evolutionary origins of modern cultivars, the pace and mode of selection underlying their adaptation to new environments and the production of phenotypes of interest. Although forests are among the most biologically rich ecosystems on earth and represent a fundamental resource for human societies, no ancient genome sequences have been generated for trees. This contrasts with the generation of multiple ancient reference genomes for important crops. Here, we sequenced the first ancient tree genomes using two white oak wood remains from Germany dating to the Last Little Ice Age (15th century CE, 7.3× and 4.0×) and one from France dating to the Bronze Age (1700 BCE, 3.4×). We assessed the underlying species and identified one medieval remains as a hybrid between two common oak species (Quercus robur and Q. petraea) and the other two remains as Q. robur. We found that diversity at the global genome level had not changed over time. However, exploratory analyses suggested that a reduction of diversity took place at different time periods. Finally, we determined the timing of leaf unfolding for ancient trees for the first time. The study extends the application of ancient wood beyond the classical proxies of dendroclimatology, dendrochronology, dendroarchaeology and dendroecology, thereby enhancing resolution of inferences on the responses of forest ecosystems to past environmental changes, epidemics and silvicultural practices.

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.

Equine herpesvirus 4 infected domestic horses associated with Sintashta spoke-wheeled chariots around 4,000 years ago

Equine viral outbreaks have disrupted the socio-economic life of past human societies up until the late 19th century and continue to be of major concern to the horse industry today. With a seroprevalence of 60-80 per cent, equine herpesvirus 4 (EHV-4) is the most common horse pathogen on the planet. Yet, its evolutionary history remains understudied. Here, we screen the sequenced data of 264 archaeological horse remains to detect the presence of EHV-4. We recover the first ancient EHV-4 genome with 4.2× average depth-of-coverage from a specimen excavated in the Southeastern Urals and dated to the Early Bronze Age period, approximately 3,900 years ago. The recovery of an EHV-4 virus outside the upper respiratory tract not only points to an animal particularly infected but also highlights the importance of post-cranial bones in pathogen characterisation. Bayesian phylogenetic reconstruction provides a minimal time estimate for EHV-4 diversification to around 4,000 years ago, a time when modern domestic horses spread across the Central Asian steppes together with spoke-wheeled Sintashta chariots, or earlier. The analyses also considerably revise the diversification time of the two EHV-4 subclades from the 16th century based solely on modern data to nearly a thousand years ago. Our study paves the way for a robust reconstruction of the history of non-human pathogens and their impact on animal health.

A modular dCas9-based recruitment platform for combinatorial epigenome editing

Targeted epigenome editing tools allow precise manipulation and investigation of genome modifications, however they often display high context dependency and variable efficacy between target genes and cell types. While systems that simultaneously recruit multiple distinct 'effector' chromatin regulators can improve efficacy, they generally lack control over effector composition and spatial organisation. To overcome this we have created a modular combinatorial epigenome editing platform, called SSSavi. This system is an interchangeable and reconfigurable docking platform fused to dCas9 that enables simultaneous recruitment of up to four different effectors, allowing precise control of effector composition and spatial ordering. We demonstrate the activity and specificity of the SSSavi system and, by testing it against existing multi-effector targeting systems, demonstrate its comparable efficacy. Furthermore, we demonstrate the importance of the spatial ordering of the recruited effectors for effective transcriptional regulation. Together, the SSSavi system enables exploration of combinatorial effector co-recruitment to enhance manipulation of chromatin contexts previously resistant to targeted editing.

The Natural and Clinical History of Plague: From the Ancient Pandemics to Modern Insights

The human pathogen Yersinia pestis is responsible for bubonic, septicemic, and pneumonic plague. A deeply comprehensive overview of its historical context, bacteriological characteristics, genomic analysis based on ancient DNA (aDNA) and modern strains, and its impact on historical and actual human populations, is explored. The results from multiple studies have been synthesized to investigate the origins of plague, its transmission, and effects on different populations. Additionally, molecular interactions of Y. pestis, from its evolutionary origins to its adaptation to flea-born transmission, and its impact on human and wild populations are considered. The characteristic combinations of aDNA patterns, which plays a decisive role in the reconstruction and analysis of ancient genomes, are reviewed. Bioinformatics is fundamental in identifying specific Y. pestis lineages, and automated pipelines are among the valuable tools in implementing such studies. Plague, which remains among human history's most lethal infectious diseases, but also other zoonotic diseases, requires the continuous investigation of plague topics. This can be achieved by improving molecular and genetic screening of animal populations, identifying ecological and social determinants of outbreaks, increasing interdisciplinary collaborations among scientists and public healthcare providers, and continued research into the characterization, diagnosis, and treatment of these diseases.

African bushpigs exhibit porous species boundaries and appeared in Madagascar concurrently with human arrival

Several African mammals exhibit a phylogeographic pattern where closely related taxa are split between West/Central and East/Southern Africa, but their evolutionary relationships and histories remain controversial. Bushpigs (Potamochoerus larvatus) and red river hogs (P. porcus) are recognised as separate species due to morphological distinctions, a perceived lack of interbreeding at contact, and putatively old divergence times, but historically, they were considered conspecific. Moreover, the presence of Malagasy bushpigs as the sole large terrestrial mammal shared with the African mainland raises intriguing questions about its origin and arrival in Madagascar. Analyses of 67 whole genomes revealed a genetic continuum between the two species, with putative signatures of historical gene flow, variable FST values, and a recent divergence time (<500,000 years). Thus, our study challenges key arguments for splitting Potamochoerus into two species and suggests their speciation might be incomplete. Our findings also indicate that Malagasy bushpigs diverged from southern African populations and underwent a limited bottleneck 1000-5000 years ago, concurrent with human arrival in Madagascar. These results shed light on the evolutionary history of an iconic and widespread African mammal and provide insight into the longstanding biogeographic puzzle surrounding the bushpig's presence in Madagascar.