The presence and impact of reference bias on population genomic studies of prehistoric human populations

PLoS Genet. 2019 Jul 26;15(7):e1008302. doi: 10.1371/journal.pgen.1008302. eCollection 2019 Jul.

Abstract

Haploid high quality reference genomes are an important resource in genomic research projects. A consequence is that DNA fragments carrying the reference allele will be more likely to map successfully, or receive higher quality scores. This reference bias can have effects on downstream population genomic analysis when heterozygous sites are falsely considered homozygous for the reference allele. In palaeogenomic studies of human populations, mapping against the human reference genome is used to identify endogenous human sequences. Ancient DNA studies usually operate with low sequencing coverages and fragmentation of DNA molecules causes a large proportion of the sequenced fragments to be shorter than 50 bp-reducing the amount of accepted mismatches, and increasing the probability of multiple matching sites in the genome. These ancient DNA specific properties are potentially exacerbating the impact of reference bias on downstream analyses, especially since most studies of ancient human populations use pseudo-haploid data, i.e. they randomly sample only one sequencing read per site. We show that reference bias is pervasive in published ancient DNA sequence data of prehistoric humans with some differences between individual genomic regions. We illustrate that the strength of reference bias is negatively correlated with fragment length. Most genomic regions we investigated show little to no mapping bias but even a small proportion of sites with bias can impact analyses of those particular loci or slightly skew genome-wide estimates. Therefore, reference bias has the potential to cause minor but significant differences in the results of downstream analyses such as population allele sharing, heterozygosity estimates and estimates of archaic ancestry. These spurious results highlight how important it is to be aware of these technical artifacts and that we need strategies to mitigate the effect. Therefore, we suggest some post-mapping filtering strategies to resolve reference bias which help to reduce its impact substantially.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Animals
  • Bias
  • DNA, Ancient / analysis*
  • Genome, Human
  • High-Throughput Nucleotide Sequencing / methods
  • Hominidae / genetics*
  • Humans
  • Metagenomics / methods*
  • Sequence Analysis, DNA / methods
  • Software

Substances

  • DNA, Ancient

Grants and funding

The computations were performed on resources provided by SNIC through Uppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX) under projects sllstore2017087, uppstore2018139, SNIC 2018/8-106 and 2018/8-239. TG was supported by grants from the Swedish Research Council Vetenskapsrådet (2017-05267, vr.se) and The Royal Physiographic Society of Lund (Nilsson-Ehle Endowments, fysiografen.se). CN was supported by a grant from Svenska Forskningsrådet Formas (Formas; 2017-00453, www.formas.se). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.