A 1 Ma sedimentary ancient DNA (sedaDNA) record of catchment vegetation changes and the developmental history of tropical Lake Towuti (Sulawesi, Indonesia).

Studying past ecosystems from ancient environmental DNA preserved in lake sediments (sedaDNA) is a rapidly expanding field. This research has mainly involved Holocene sediments from lakes in cool climates, with little known about the suitability of sedaDNA to reconstruct substantially older ecosystems in the warm tropics. Here, we report the successful recovery of chloroplast trnL (UAA) sequences (trnL-P6 loop) from the sedimentary record of Lake Towuti (Sulawesi, Indonesia) to elucidate changes in regional tropical vegetation assemblages during the lake's Late Quaternary paleodepositional history. After the stringent removal of contaminants and sequence artifacts, taxonomic assignment of the remaining genuine trnL-P6 reads showed that native nitrogen-fixing legumes, C3 grasses, and shallow wetland vegetation (Alocasia) were most strongly associated with >1-million-year-old (>1 Ma) peats and silts (114-98.8 m composite depth; mcd), which were deposited in a landscape of active river channels, shallow lakes, and peat-swamps. A statistically significant shift toward partly submerged shoreline vegetation that was likely rooted in anoxic muddy soils (i.e., peatland forest trees and wetland C3 grasses (Oryzaceae) and nutrient-demanding aquatic herbs (presumably Oenanthe javanica)) occurred at 76 mcd (~0.8 Ma), ~0.2 Ma after the transition into a permanent lake. This wetland vegetation was most strongly associated with diatom ooze (46-37 mcd), thought to be deposited during maximum nutrient availability and primary productivity. Herbs (Brassicaceae), trees/shrubs (Fabaceae and Theaceae), and C3 grasses correlated with inorganic parameters, indicating increased drainage of ultramafic sediments and laterite soils from the lakes' catchment, particularly at times of inferred drying. Downcore variability in trnL-P6 from tropical forest trees (Toona), shady ground cover herbs (Zingiberaceae), and tree orchids (Luisia) most strongly correlated with sediments of a predominantly felsic signature considered to be originating from the catchment of the Loeha River draining into Lake Towuti during wetter climate conditions. However, the co-correlation with dry climate-adapted trees (i.e., Castanopsis or Lithocarpus) plus C4 grasses suggests that increased precipitation seasonality also contributed to the increased drainage of felsic Loeha River sediments. This multiproxy approach shows that despite elevated in situ temperatures, tropical lake sediments potentially comprise long-term archives of ancient environmental DNA for reconstructing ecosystems, which warrants further exploration.

Phylogenetic insights into the genetic legacies of Hungarian-speaking communities in the Carpathian Basin.

This study focuses on exploring the uniparental genetic lineages of Hungarian-speaking minorities residing in rural villages of Baranja (Croatia) and the Zobor region (Slovakia). We aimed to identify ancestral lineages by examining genetic markers distributed across the entire mitogenome and on the Y-chromosome. This allowed us to discern disparities in regional genetic structures within these communities. By integrating our newly acquired genetic data from a total of 168 participants with pre-existing Eurasian and ancient DNA datasets, our goal was to enrich the understanding of the genetic history trajectories of Carpathian Basin populations. Our findings suggest that while population-based analyses may not be sufficiently robust to detect fine-scale uniparental genetic patterns with the sample sizes at hand, phylogenetic analysis of well-characterized Y-chromosomal Short Tandem Repeat (STR) data and entire mitogenome sequences did uncover multiple lineage ties to far-flung regions and eras. While the predominant portions of both paternal and maternal DNA align with the East-Central European spectrum, rarer subhaplogroups and lineages have unveiled ancient ties to both prehistoric and historic populations spanning Europe and Eastern Eurasia. This research augments the expansive field of phylogenetics, offering critical perspectives on the genetic constitution and heritage of the communities in East-Central Europe.

Ancient Mycobacterium leprae genome reveals medieval English red squirrels as animal leprosy host.

Leprosy, one of the oldest recorded diseases in human history, remains prevalent in Asia, Africa, and South America, with over 200,000 cases every year.1,2 Although ancient DNA (aDNA) approaches on the major causative agent, Mycobacterium leprae, have elucidated the disease's evolutionary history,3,4,5 the role of animal hosts and interspecies transmission in the past remains unexplored. Research has uncovered relationships between medieval strains isolated from archaeological human remains and modern animal hosts such as the red squirrel in England.6,7 However, the time frame, distribution, and direction of transmissions remains unknown. Here, we studied 25 human and 12 squirrel samples from two archaeological sites in Winchester, a medieval English city well known for its leprosarium and connections to the fur trade. We reconstructed four medieval M. leprae genomes, including one from a red squirrel, at a 2.2-fold average coverage. Our analysis revealed a phylogenetic placement of all strains on branch 3 as well as a close relationship between the squirrel strain and one newly reconstructed medieval human strain. In particular, the medieval squirrel strain is more closely related to some medieval human strains from Winchester than to modern red squirrel strains from England, indicating a yet-undetected circulation of M. leprae in non-human hosts in the Middle Ages. Our study represents the first One Health approach for M. leprae in archaeology, which is centered around a medieval animal host strain, and highlights the future capability of such approaches to understand the disease's zoonotic past and current potential.

New Canary Islands Roman mediated settlement hypothesis deduced from coalescence ages of curated maternal indigenous lineages.

Numerous genetic studies have contributed to reconstructing the human history of the Canary Islands population. The recent use of new ancient DNA targeted enrichment and next-generation sequencing techniques on new Canary Islands samples have greatly improved these molecular results. However, the bulk of the available data is still provided by the classic mitochondrial DNA phylogenetic and phylogeographic studies carried out on the indigenous, historical, and extant human populations of the Canary Islands. In the present study, making use of all the accumulated mitochondrial information, the existence of DNA contamination and archaeological sample misidentification in those samples is evidenced. Following a thorough review of these cases, the new phylogeographic analysis revealed the existence of a heterogeneous indigenous Canarian population, asymmetrically distributed across the various islands, which most likely descended from a unique mainland settlement. These new results and new proposed coalescent ages are compatible with a Roman-mediated arrival driven by the exploitation of the purple dye manufacture in the Canary Islands.

A High-Coverage Mesolithic Aurochs Genome and Effective Leveraging of Ancient Cattle Genomes Using Whole Genome Imputation.

Ancient genomic analyses are often restricted to utilizing pseudohaploid data due to low genome coverage. Leveraging low-coverage data by imputation to calculate phased diploid genotypes that enables haplotype-based interrogation and single nucleotide polymorphism (SNP) calling at unsequenced positions is highly desirable. This has not been investigated for ancient cattle genomes despite these being compelling subjects for archeological, evolutionary, and economic reasons. Here, we test this approach by sequencing a Mesolithic European aurochs (18.49×; 9,852 to 9,376 calBCE) and an Early Medieval European cow (18.69×; 427 to 580 calCE) and combine these with published individuals: two ancient and three modern. We downsample these genomes (0.25×, 0.5×, 1.0×, and 2.0×) and impute diploid genotypes, utilizing a reference panel of 171 published modern cattle genomes that we curated for 21.7 million (Mn) phased SNPs. We recover high densities of correct calls with an accuracy of >99.1% at variant sites for the lowest downsample depth of 0.25×, increasing to >99.5% for 2.0× (transversions only, minor allele frequency [MAF] ≥ 2.5%). The recovery of SNPs correlates with coverage; on average, 58% of sites are recovered for 0.25× increasing to 87% for 2.0×, utilizing an average of 3.5 million (Mn) transversions (MAF ≥2.5%), even in the aurochs, despite the highest temporal distance from the modern reference panel. Our imputed genomes behave similarly to directly called data in allele frequency-based analyses, for example consistently identifying runs of homozygosity >2 Mb, including a long homozygous region in the Mesolithic European aurochs.

Network of large pedigrees reveals social practices of Avar communities.

From AD 567-568, at the onset of the Avar period, populations from the Eurasian Steppe settled in the Carpathian Basin for approximately 250 years1. Extensive sampling for archaeogenomics (424 individuals) and isotopes, combined with archaeological, anthropological and historical contextualization of four Avar-period cemeteries, allowed for a detailed description of the genomic structure of these communities and their kinship and social practices. We present a set of large pedigrees, reconstructed using ancient DNA, spanning nine generations and comprising around 300 individuals. We uncover a strict patrilineal kinship system, in which patrilocality and female exogamy were the norm and multiple reproductive partnering and levirate unions were common. The absence of consanguinity indicates that this society maintained a detailed memory of ancestry over generations. These kinship practices correspond with previous evidence from historical sources and anthropological research on Eurasian Steppe societies2. Network analyses of identity-by-descent DNA connections suggest that social cohesion between communities was maintained via female exogamy. Finally, despite the absence of major ancestry shifts, the level of resolution of our analyses allowed us to detect genetic discontinuity caused by the replacement of a community at one of the sites. This was paralleled with changes in the archaeological record and was probably a result of local political realignment.

Ancient DNA indicates a century of overhunting did not reduce genetic diversity in Pacific Walruses (Odobenus rosmarus divergens).

Pacific Walruses (Odobenus rosmarus divergens [Illiger 1815]) are gregarious marine mammals considered to be sentinels of the Arctic because of their dependence on sea ice for feeding, molting, and parturition. Like many other marine mammal species, their population sizes were decimated by historical overhunting in the nineteenth and twentieth centuries. Although they have since been protected from nearly all commercial hunting pressure, they now face rapidly accelerating habitat loss as global warming reduces the extent of summer sea ice in the Arctic. To investigate how genetic variation was impacted by overhunting, we obtained mitochondrial DNA sequences from historic Pacific Walrus samples in Alaska that predate the period of overhunting, as well as from extant populations. We found that genetic variation was unchanged over this period, suggesting Pacific Walruses are resilient to genetic attrition in response to reduced population size, and that this may be related to their high vagility and lack of population structure. Although Pacific Walruses will almost certainly continue to decline in number as the planet warms and summer sea ice is further reduced, they may be less susceptible to the ratcheting effects of inbreeding that typically accompany shrinking populations.

Palaeogenomic insights into the origins of early settlers on the island of Cyprus.

Archaeological evidence supports sporadic seafaring visits to the Eastern Mediterranean island of Cyprus by Epipaleolithic hunter-gatherers over 12,000 years ago, followed by permanent settlements during the early Neolithic. The geographical origins of these early seafarers have so far remained elusive. By systematically analysing all available genomes from the late Pleistocene to early Holocene Near East (c. 14,000-7000 cal BCE), we provide a comprehensive overview of the genetic landscape of the early Neolithic Fertile Crescent and Anatolia and infer the likely origins of three recently published genomes from Kissonerga-Mylouthkia (Cypriot Late Pre-Pottery Neolithic B, c. 7600-6800 cal BCE). These appear to derive roughly 80% of their ancestry from Aceramic Neolithic Central Anatolians residing in or near the Konya plain, and the remainder from a genetically basal Levantine population. Based on genome-wide weighted ancestry covariance analysis, we infer that this admixture event took place roughly between 14,000 and 10,000 BCE, coinciding with the transition from the Cypriot late Epipaleolithic to the Pre-Pottery Neolithic A (PPNA). Additionally, we identify strong genetic affinities between the examined Cypro-LPPNB individuals and later northwestern Anatolians and the earliest European Neolithic farmers. Our results inform archaeological evidence on prehistoric demographic processes in the Eastern Mediterranean, providing important insights into early seafaring, maritime connections, and insular settlement.

Using ancient sedimentary DNA to forecast ecosystem trajectories under climate change.

Ecosystem response to climate change is complex. In order to forecast ecosystem dynamics, we need high-quality data on changes in past species abundance that can inform process-based models. Sedimentary ancient DNA (sedaDNA) has revolutionised our ability to document past ecosystems' dynamics. It provides time series of increased taxonomic resolution compared to microfossils (pollen, spores), and can often give species-level information, especially for past vascular plant and mammal abundances. Time series are much richer in information than contemporary spatial distribution information, which have been traditionally used to train models for predicting biodiversity and ecosystem responses to climate change. Here, we outline the potential contribution of sedaDNA to forecast ecosystem changes. We showcase how species-level time series may allow quantification of the effect of biotic interactions in ecosystem dynamics, and be used to estimate dispersal rates when a dense network of sites is available. By combining palaeo-time series, process-based models, and inverse modelling, we can recover the biotic and abiotic processes underlying ecosystem dynamics, which are traditionally very challenging to characterise. Dynamic models informed by sedaDNA can further be used to extrapolate beyond current dynamics and provide robust forecasts of ecosystem responses to future climate change. This article is part of the theme issue 'Ecological novelty and planetary stewardship: biodiversity dynamics in a transforming biosphere'.

PNPLA3 fatty liver allele was fixed in Neanderthals and segregates neutrally in humans.

OBJECTIVE: Fat deposition is modulated by environmental factors and genetic predisposition. Genome-wide association studies identified PNPLA3 p.I148M (rs738409) as a common variant that increases risk of developing liver steatosis. When and how this variant evolved in humans has not been studied to date. DESIGN: Here we analyse ancient DNA to track the history of this allele throughout human history. In total, 6444 published ancient (modern humans, Neanderthal, Denisovan) and 3943 published present day genomes were used for analysis after extracting genotype calls for PNPLA3 p.I148M. To quantify changes through time, logistic and, by grouping individuals according to geography and age, linear regression analyses were performed. RESULTS: We find that archaic human individuals (Neanderthal, Denisovan) exclusively carried a fixed PNPLA3 risk allele, whereas allele frequencies in modern human populations range from very low in Africa to >50% in Mesoamerica. Over the last 15 000 years, distributions of ancestral and derived alleles roughly match the present day distribution. Logistic regression analyses did not yield signals of natural selection during the last 10 000 years. CONCLUSION: Archaic human individuals exclusively carried a fixed PNPLA3 allele associated with fatty liver, whereas allele frequencies in modern human populations are variable even in the oldest samples. Our observation might underscore the advantage of fat storage in cold climate and particularly for Neanderthal under ice age conditions. The absent signals of natural selection during modern human history does not support the thrifty gene hypothesis in case of PNPLA3 p.I148M.

Comparison and optimization of protocols and whole-genome capture conditions for ancient DNA samples.

Ancient DNA (aDNA) obtained from human remains is typically fragmented and present in relatively low amounts. Here we investigate a set of optimal methods for producing aDNA data by comparing silica-based DNA extraction and aDNA library preparation protocols. We also test the efficiency of whole-genome enrichment (WGC) on ancient human samples by modifying a number of parameter combinations. We find that the Dabney extraction protocol performs significantly better than alternatives. We further observed a positive trend with the BEST library protocol indicating lower clonality. Notably, our results suggest that WGC is effective at retrieving endogenous DNA, particularly from poorly-preserved human samples, by increasing human endogenous proportions by 5x. Thus, aDNA studies will be most likely to benefit from our results.

Ancient genomes and the evolutionary path of modern humans.

Growing evidence from archaic and early modern human genomes brings new insights to the emergence of modern humans. We recount recent information collected from ancient DNA studies that inform us about the evolutionary pathway to modern humanity. These findings point to both individual- and population-level advantages underlying modern human expansion.

Benchmarking a targeted 16S ribosomal RNA gene enrichment approach to reconstruct ancient microbial communities.

The taxonomic characterization of ancient microbiomes is a key step in the rapidly growing field of paleomicrobiology. While PCR amplification of the 16S ribosomal RNA (rRNA) gene is a widely used technique in modern microbiota studies, this method has systematic biases when applied to ancient microbial DNA. Shotgun metagenomic sequencing has proven to be the most effective method in reconstructing taxonomic profiles of ancient dental calculus samples. Nevertheless, shotgun sequencing approaches come with inherent limitations that could be addressed through hybridization enrichment capture. When employed together, shotgun sequencing and hybridization capture have the potential to enhance the characterization of ancient microbial communities. Here, we develop, test, and apply a hybridization enrichment capture technique to selectively target 16S rRNA gene fragments from the libraries of ancient dental calculus samples generated with shotgun techniques. We simulated data sets generated from hybridization enrichment capture, indicating that taxonomic identification of fragmented and damaged 16S rRNA gene sequences was feasible. Applying this enrichment approach to 15 previously published ancient calculus samples, we observed a 334-fold increase of ancient 16S rRNA gene fragments in the enriched samples when compared to unenriched libraries. Our results suggest that 16S hybridization capture is less prone to the effects of background contamination than 16S rRNA amplification, yielding a higher percentage of on-target recovery. While our enrichment technique detected low abundant and rare taxa within a given sample, these assignments may not achieve the same level of specificity as those achieved by unenriched methods.

Ancient DNA and osteological analyses of a unique paleo-archive reveal Early Holocene faunal expansion into the Scandinavian Arctic.

Paleo-archives are essential for our understanding of species responses to climate warming, yet such archives are extremely rare in the Arctic. Here, we combine morphological analyses and bulk-bone metabarcoding to investigate a unique chronology of bone deposits sealed in the high-latitude Storsteinhola cave system (68°50' N 16°22' E) in Norway. This deposit dates to a period of climate warming from the end of the Late Glacial [~13 thousand calibrated years before the present (ka cal B.P.)] to the Holocene thermal maximum (~5.6 ka cal B.P.). Paleogenetic analyses allow us to exploit the 1000s of morphologically unidentifiable bone fragments resulting in a high-resolution sequence with 40 different taxa, including species not previously found here. Our record reveals borealization in both the marine and terrestrial environments above the Arctic Circle as a naturally recurring phenomenon in past periods of warming, providing fundamental insights into the ecosystem-wide responses that are ongoing today.

Assessing the impact of post-mortem damage and contamination on imputation performance in ancient DNA.

Low-coverage imputation is becoming ever more present in ancient DNA (aDNA) studies. Imputation pipelines commonly used for present-day genomes have been shown to yield accurate results when applied to ancient genomes. However, post-mortem damage (PMD), in the form of C-to-T substitutions at the reads termini, and contamination with DNA from closely related species can potentially affect imputation performance in aDNA. In this study, we evaluated imputation performance (i) when using a genotype caller designed for aDNA, ATLAS, compared to bcftools, and (ii) when contamination is present. We evaluated imputation performance with principal component analyses and by calculating imputation error rates. With a particular focus on differently imputed sites, we found that using ATLAS prior to imputation substantially improved imputed genotypes for a very damaged ancient genome (42% PMD). Trimming the ends of the sequencing reads led to similar improvements in imputation accuracy. For the remaining genomes, ATLAS brought limited gains. Finally, to examine the effect of contamination on imputation, we added various amounts of reads from two present-day genomes to a previously downsampled high-coverage ancient genome. We observed that imputation accuracy drastically decreased for contamination rates above 5%. In conclusion, we recommend (i) accounting for PMD by either trimming sequencing reads or using a genotype caller such as ATLAS before imputing highly damaged genomes and (ii) only imputing genomes containing up to 5% of contamination.

Ancient mitogenomes suggest complex maternal history of one of the oldest settlements of western India.

The ancient township of Vadnagar tells a story of a long chain of cultural diversity and exchange. Vadnagar has been continuously habituated and shows a presence of rich cultural amalgamation and continuous momentary sequences between the 2th century BCE and present-day. Seven cultural periods developed a complex and enriched settlement at Vadnagar in spatio-temporality. Although archaeological studies done on this oldest settlement suggested a rich cultural heritage, the genetic studies to pinpoint the genetic ancestry was lacking till date. In our current study we have for the first time reconstructed the complete mitogenomes of medieval individuals of the Vadnagar archaeological site in Gujarat. The study aimed to investigate the cosmopolitan nature of the present population as well as the migratory pattern and the inflow of different groups through trade, cultural and religious practices. Our analysis suggests heterogeneous nature of the medieval population of Vadnagar with presence of deeply rooted local ancestral components as well as central Asian genetic ancestry. This Central Asian component associated with mitochondrial haplotype U2e was not shared with any individual from India, but rather with individuals from the Bronze Age of Tajikistan and with an earlier age of coalescence. In summary, we propose that the medieval site of Vadnagar in western India was rich in cultural and genetic aspects, with both local and western Eurasian components.

Mitochondrial genomes of Pleistocene megafauna retrieved from recent sediment layers of two Siberian lakes.

Ancient environmental DNA (aeDNA) from lake sediments has yielded remarkable insights for the reconstruction of past ecosystems, including suggestions of late survival of extinct species. However, translocation and lateral inflow of DNA in sediments can potentially distort the stratigraphic signal of the DNA. Using three different approaches on two short lake sediment cores of the Yamal peninsula, West Siberia, with ages spanning only the past hundreds of years, we detect DNA and identified mitochondrial genomes of multiple mammoth and woolly rhinoceros individuals-both species that have been extinct for thousands of years on the mainland. The occurrence of clearly identifiable aeDNA of extinct Pleistocene megafauna (e.g. >400 K reads in one core) throughout these two short subsurface cores, along with specificities of sedimentology and dating, confirm that processes acting on regional scales, such as extensive permafrost thawing, can influence the aeDNA record and should be accounted for in aeDNA paleoecology.

Predicting Functional Consequences of Recent Natural Selection in Britain.

Ancient DNA can directly reveal the contribution of natural selection to human genomic variation. However, while the analysis of ancient DNA has been successful at identifying genomic signals of selection, inferring the phenotypic consequences of that selection has been more difficult. Most trait-associated variants are noncoding, so we expect that a large proportion of the phenotypic effects of selection will also act through noncoding variation. Since we cannot measure gene expression directly in ancient individuals, we used an approach (Joint-Tissue Imputation [JTI]) developed to predict gene expression from genotype data. We tested for changes in the predicted expression of 17,384 protein coding genes over a time transect of 4,500 years using 91 present-day and 616 ancient individuals from Britain. We identified 28 genes at seven genomic loci with significant (false discovery rate [FDR] < 0.05) changes in predicted expression levels in this time period. We compared the results from our transcriptome-wide scan to a genome-wide scan based on estimating per-single nucleotide polymorphism (SNP) selection coefficients from time series data. At five previously identified loci, our approach allowed us to highlight small numbers of genes with evidence for significant shifts in expression from peaks that in some cases span tens of genes. At two novel loci (SLC44A5 and NUP85), we identify selection on gene expression not captured by scans based on genomic signatures of selection. Finally, we show how classical selection statistics (iHS and SDS) can be combined with JTI models to incorporate functional information into scans that use present-day data alone. These results demonstrate the potential of this type of information to explore both the causes and consequences of natural selection.

Ancient genomes illuminate Eastern Arabian population history and adaptation against malaria.

The harsh climate of Arabia has posed challenges in generating ancient DNA from the region, hindering the direct examination of ancient genomes for understanding the demographic processes that shaped Arabian populations. In this study, we report whole-genome sequence data obtained from four Tylos-period individuals from Bahrain. Their genetic ancestry can be modeled as a mixture of sources from ancient Anatolia, Levant, and Iran/Caucasus, with variation between individuals suggesting population heterogeneity in Bahrain before the onset of Islam. We identify the G6PD Mediterranean mutation associated with malaria resistance in three out of four ancient Bahraini samples and estimate that it rose in frequency in Eastern Arabia from 5 to 6 kya onward, around the time agriculture appeared in the region. Our study characterizes the genetic composition of ancient Arabians, shedding light on the population history of Bahrain and demonstrating the feasibility of studies of ancient DNA in the region.