An efficient pipeline for ancient DNA mapping and recovery of endogenous ancient DNA from whole-genome sequencing data

Searching for alternative high DNA-yielding bone types for DNA analysis of aged skeletal remains

The petrous bone contains significantly higher amounts of DNA than any other human bone. Because of highly destructive sampling and because it is not always part of the recovered remains, the need for alternative sources of DNA is important. To identify additional optimal bone types, petrous bones were compared to femurs, tali, and calcanei sampled from 66 adult skeletons from two distinct modern-era Christian cemeteries. An extraction method employing full demineralization was used to obtain DNA, real-time PCR quantification to ascertain DNA quantity and degradation, and a commercial forensic short tandem repeats (STR) PCR amplification kit to determine genetic profiles. Statistical analysis was performed to explore the differences in DNA yield, DNA degradation, and success of STR amplification. A systematic studies exploring intra-skeletal variability in DNA preservation including various excavation sites differing by time period and geographical position are rare, and the second part of the investigation was based on a comparison of both archaeological sites, which allowed us to compare the effect of different post-mortem intervals and environmental conditions on DNA preservation. The older burial site in Črnomelj was active between the 13th and 18th century, whereas the more recent Polje burial was in use from the 16th to 19th century, creating different temporal and geographical environments. Results for the Črnomelj burial site revealed that the petrous bone outperformed all other bone types studied, except the calcaneus. At the Polje archeological site calcanei, tali, and femurs yielded the same STR typing success as petrous bones. The results obtained highlight the importance of careful bone sample selection for DNA analysis of aged skeletal remains. In addition to petrous bones, calcanei were found to be an alternative source of DNA when older burial sites are investigated. When more recent burial sites are processed, calcanei, tali, and femurs should be sampled besides petrous bones, not only because they exhibited good performance, but also because of easier sampling and easier grinding in the case of trabecular bones. This study contributes valuable insights into the potential use of various skeletal types as a source of DNA for investigation of aged skeletal remains, and it offers practical implications for forensic and archaeological investigations.

Ancient genomes illuminate the demographic history of Shandong over the past two millennia

Shandong province, located in the Lower Yellow River, is one of the birthplaces of ancient Chinese civilization. However, the comprehensive genetic histories of this region have remained largely unknown until now due to a lack of ancient human genomes. Here, we present 21 ancient genomes from Shandong dating from the Warring States period to the Jin-Yuan Dynasties. Unlike the early Neolithic samples from Shandong, the historical samples are most closely related to post-Late Neolithic populations of the Middle Yellow River Basin, suggesting a population turnover in Shandong from the Neolithic Age to the Historical era. In addition, we detect a close genetic affinity between the historical samples in Shandong and present-day Han Chinese, showing long-term genetic stability in Han Chinese at least since the Warring States period.

Unravelling reference bias in ancient DNA datasets

MOTIVATION

The alignment of sequencing reads is a critical step in the characterization of ancient genomes. However, reference bias and spurious mappings pose a significant challenge, particularly as cutting-edge wet lab methods generate datasets that push the boundaries of alignment tools. Reference bias occurs when reference alleles are favoured over alternative alleles during mapping, whereas spurious mappings stem from either contamination or when endogenous reads fail to align to their correct position. Previous work has shown that these phenomena are correlated with read length but a more thorough investigation of reference bias and spurious mappings for ancient DNA has been lacking. Here, we use a range of empirical and simulated palaeogenomic datasets to investigate the impacts of mapping tools, quality thresholds, and reference genome on mismatch rates across read lengths.

RESULTS

For these analyses, we introduce AMBER, a new bioinformatics tool for assessing the quality of ancient DNA mapping directly from BAM-files and informing on reference bias, read length cut-offs and reference selection. AMBER rapidly and simultaneously computes the sequence read mapping bias in the form of the mismatch rates per read length, cytosine deamination profiles at both CpG and non-CpG sites, fragment length distributions, and genomic breadth and depth of coverage. Using AMBER, we find that mapping algorithms and quality threshold choices dictate reference bias and rates of spurious alignment at different read lengths in a predictable manner, suggesting that optimized mapping parameters for each read length will be a key step in alleviating reference bias and spurious mappings.

AVAILABILITY AND IMPLEMENTATION

AMBER is available for noncommercial use on GitHub (https://github.com/tvandervalk/AMBER.git). Scripts used to generate and analyse simulated datasets are available on Github (https://github.com/sdolenz/refbias_scripts).

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.

The genomic prehistory of the Indigenous peoples of Uruguay

The prehistory of the people of Uruguay is greatly complicated by the dramatic and severe effects of European contact, as with most of the Americas. After the series of military campaigns that exterminated the last remnants of nomadic peoples, Uruguayan official history masked and diluted the former Indigenous ethnic diversity into the narrative of a singular people that all but died out. Here, we present the first whole genome sequences of the Indigenous people of the region before the arrival of Europeans, from an archaeological site in eastern Uruguay that dates from 2,000 years before present. We find a surprising connection to ancient individuals from Panama and eastern Brazil, but not to modern Amazonians. This result may be indicative of a migration route into South America that may have occurred along the Atlantic coast. We also find a distinct ancestry previously undetected in South America. Though this work begins to piece together some of the demographic nuance of the region, the sequencing of ancient individuals from across Uruguay is needed to better understand the ancient prehistory and genetic diversity that existed before European contact, thereby helping to rebuild the history of the Indigenous population of what is now Uruguay.

Additional evaluations show that specific BWA-aln settings still outperform BWA-mem for ancient DNA data alignment

Xu et al. (2021) recently recommended a new parameterization of BWA-mem as a superior alternative to the widely-used BWA-aln algorithm to map ancient DNA sequencing data. Here, we compare the BWA-mem parameterization recommended by Xu et al. with the best-performing alignment methods determined in the recent benchmarks of Oliva and colleagues (2021), demonstrating that BWA-aln is still the gold-standard for ancient DNA read alignment .