Mapping Bias Overestimates Reference Allele Frequencies at the HLA Genes in the 1000 Genomes Project Phase I Data

Abstract

Next-generation sequencing (NGS) technologies have become the standard for data generation in studies of population genomics, as the 1000 Genomes Project (1000G). However, these techniques are known to be problematic when applied to highly polymorphic genomic regions, such as the human leukocyte antigen (HLA) genes. Because accurate genotype calls and allele frequency estimations are crucial to population genomics analyses, it is important to assess the reliability of NGS data. Here, we evaluate the reliability of genotype calls and allele frequency estimates of the single-nucleotide polymorphisms (SNPs) reported by 1000G (phase I) at five HLA genes (HLA-A, -B, -C, -DRB1, and -DQB1). We take advantage of the availability of HLA Sanger sequencing of 930 of the 1092 1000G samples and use this as a gold standard to benchmark the 1000G data. We document that 18.6% of SNP genotype calls in HLA genes are incorrect and that allele frequencies are estimated with an error greater than ±0.1 at approximately 25% of the SNPs in HLA genes. We found a bias toward overestimation of reference allele frequency for the 1000G data, indicating mapping bias is an important cause of error in frequency estimation in this dataset. We provide a list of sites that have poor allele frequency estimates and discuss the outcomes of including those sites in different kinds of analyses. Because the HLA region is the most polymorphic in the human genome, our results provide insights into the challenges of using of NGS data at other genomic regions of high diversity.

Shown for reading. Its publisher allows quoting it, not republishing it.

References 27

1. 1.

   Andersen K. G. Shylakhter I. Tabrizi S. Grossman S. R. Happi C. T., 2012 Genome-wide scans provide evidence for positive selection of genes implicated in Lassa fever. Philos. Trans. R. Soc. Lond. B Biol. Sci. 367: 868–877.

   Not in our library10.1098/rstb.2011.0299

2. 2.

   Bjorkman P. J. Saper M. A. Samraoui B. Bennett W. S. Strominger J. L., 1987 Structure of the human class I histocompatibility antigen, HLA-A2. Nature 329: 506–512.

   Not in our library10.1038/329506a0

3. 3.
   HLA typing from RNA-Seq sequence reads · 2012

   No notice
4. 4.

   Brown J. H. Jardetzky T. S. Gorga J. C. Stern L. J. Urban R. G., 1993 Three-dimensional structure of the human class II histocompatibility antigen HLA-DR1. Nature 364: 33–39.

   Not in our library10.1038/364033a0

5. 5.

   Chapman S. J. Hill A. V. S., 2012 Human genetic susceptibility to infectious disease. Nat. Rev. Genet. 13: 175–188.

   Not in our library10.1038/nrg3114

6. 6.

   Danecek P. Auton A. Abecasis G. R. Albers C. a. Banks E., 2011 The variant call format and VCFtools. Bioinformatics 27: 2156–2158.

   Not in our library10.1093/bioinformatics/btr330

7. 7.

   De Santis D. Dinauer D. Duke J. Erlich H. A. Holcomb C. L., 2013 16(th) IHIW: review of HLA typing by NGS. Int. J. Immunogenet. 40: 72–76.

   Not in our library10.1111/iji.12024

8. 8.

   Degner J. F. Marioni J. C. Pai A. A. Pickrell J. K. Nkadori E., 2009 Effect of read-mapping biases on detecting allele-specific expression from RNA-sequencing data. Bioinformatics 25: 3207–3212.

   Not in our library10.1093/bioinformatics/btp579

9. 9.

   Dilthey, A., C. Cox, Z. Iqbal, M. R. Nelson, and G. McVean, 2014 Improved genome inference in the MHC using a population reference graph. bioRxiv. Available from: http://biorxiv.org/content/early/2014/07/08/006973. Accessed March 20, 2015.

   Not in our library10.1038/ng.3257

10. 10.

    Erlich R. L. Jia X. Anderson S. Banks E. Gao X., 2011 Next-generation sequencing for HLA typing of class I loci. BMC Genomics 12: 42.

    Not in our library10.1186/1471-2164-12-42

11. 11.
    HLA Diversity in the 1000 Genomes Dataset · 2014

    No notice
12. 12.

    Helmberg W. Feolo M. Dunivin R. Hoffman D., 2014 dbMHC.

    Not in our library

13. 13.

    Hernandez R. D. Kelley J. L. Elyashiv E. Melton S. C. Auton A., 2011 Classic selective sweeps were rare in recent human evolution. Science 331: 920–924.

    Not in our library10.1126/science.1198878

14. 14.

    Hill-Burns E. M. Factor S. A. Zabetian C. P. Thomson G. Payami H., 2011 Evidence for more than one Parkinson’s disease-associated variant within the HLA region. PLoS One 6: e27109.

    Not in our library10.1371/journal.pone.0027109

15. 15.

    Kitts A. Feolo M. Helmberg W., 2003 The major histocompatibility complex database, dbMHC. In: National Center for Biotechnology Information NIH, ed. The NCBI Handbook. Bethesda: National Center for Biotechnology Information NIH, p.1–29.

    Not in our library

16. 16.

    Lappalainen T. Sammeth M. Friedländer M. R. ’t Hoen P. C. Monlong J., 2013 Transcriptome and genome sequencing uncovers functional variation in humans. Nature 501: 506–511.

    Not in our library10.1038/nature12531

17. 17.

    Leffler E. M. Gao Z. Pfeifer S. Ségurel L. Auton A., 2013 Multiple instances of ancient balancing selection shared between humans and chimpanzees. Science 339: 1578–1582.

    Not in our library10.1126/science.1234070

18. 18.
    HLA typing from 1000 genomes whole genome and whole exome illumina data. · 2013

    Not checked
19. 19.
    Nomenclature for factors of the HLA system, 2010 · 2010

    No notice
20. 20.

    Meyer D. Thomson G., 2001 How selection shapes variation of the human major histocompatibility complex: a review. Ann. Hum. Genet. 65: 1–26.

    Not in our library10.1046/j.1469-1809.2001.6510001.x

21. 21.

    Nielsen R. Paul J. S. Albrechtsen A. Song Y. S., 2011 Genotype and SNP calling from next-generation sequencing data. Nat. Rev. Genet. 12: 443–451.

    Not in our library10.1038/nrg2986

22. 22.

    Quinlan A. R. Hall I. M., 2010 BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26: 841–842.

    Not in our library10.1093/bioinformatics/btq033

23. 23.

    Robinson J. Halliwell J. A. McWilliam H. Lopez R. Parham P., 2013 The IMGT/HLA database. Nucleic Acids Res. 41: D1222–D1227.

    Not in our library10.1093/nar/gks949

24. 24.

    Sollid L. M. Pos W. Wucherpfennig K. W., 2014 Molecular mechanisms for contribution of MHC molecules to autoimmune diseases. Curr. Opin. Immunol. 31C: 24–30.

    Not in our library10.1016/j.coi.2014.08.005

25. 25.

    The 1000 Genomes Project Consortium, 2012 An integrated map of genetic variation from 1,092 human genomes. Nature 491: 56–65.

    Not in our library10.1038/nature11632

26. 26.

    Treangen T. J. Salzberg S. L., 2012 Repetitive DNA and next-generation sequencing: computational challenges and solutions. Nat. Rev. Genet. 13: 36–46.

    Not in our library10.1038/nrg3117

27. 27.

    Ward L. D. Kellis M., 2012 Evidence of abundant purifying selection in humans for recently acquired regulatory functions. Science 337: 1675–1678.

    Not in our library10.1126/science.1225057

Cited by papers we hold 9

• Reference flow: reducing reference bias using multiple population genomes · 2021
  No notice
• Nanopore sequencing and the Shasta toolkit enable efficient de novo assembly of eleven human genomes · 2020
  No notice
• TGS-GapCloser: A fast and accurate gap closer for large genomes with low coverage of error-prone long reads · 2020
  No notice
• Developing a High-Quality Linkage Map for the Atlantic Killifish Fundulus heteroclitus · 2019
  Not checked
• Paragraph: a graph-based structural variant genotyper for short-read sequence data · 2019
  No notice
• The presence and impact of reference bias on population genomic studies of prehistoric human populations · 2019
  No notice
• Long reads: their purpose and place · 2018
  No notice
• A genomic perspective on HLA evolution · 2017
  No notice
• Genome graphs and the evolution of genome inference · 2017
  No notice

Reads

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.