Advanced search
Start date
(Reference retrieved automatically from Web of Science through information on FAPESP grant and its corresponding number as mentioned in the publication by the authors.)

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

Full text
Brandt, Debora Y. C. [1] ; Aguiar, Vitor R. C. [1] ; Bitarello, Barbara D. [1] ; Nunes, Kelly [1] ; Goudet, Jerome [2] ; Meyer, Diogo [1]
Total Authors: 6
[1] Univ Sao Paulo, Dept Genet & Evolutionary Biol, BR-05508090 Sao Paulo, SP - Brazil
[2] Univ Lausanne, Dept Ecol & Evolut, Biophore, CH-1015 Lausanne - Switzerland
Total Affiliations: 2
Document type: Journal article
Source: G3-GENES, GENOMES, GENETICS; v. 5, n. 5, p. 931-941, MAY 1 2015.
Web of Science Citations: 31

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. (AU)

FAPESP's process: 14/12123-2 - Expression and eQTL mapping of HLA genes: analyses based on large-scale RNAseq assays
Grantee:Vitor Rezende da Costa Aguiar
Support type: Scholarships in Brazil - Post-Doctorate
FAPESP's process: 11/12500-2 - Maladaptation as a byproduct of adaptation: a genomic scale study
Grantee:Bárbara Domingues Bitarello
Support type: Scholarships in Brazil - Doctorate
FAPESP's process: 12/22796-9 - Population differentiation on genes under strong balancing selection: a case study on the HLA genes
Grantee:Débora Yoshihara Caldeira Brandt
Support type: Scholarships in Brazil - Master
FAPESP's process: 12/09950-9 - Evolution of HLA genes: population differentiation and signatures of recent selection in native and admixed populations from Brazil
Grantee:Kelly Nunes
Support type: Scholarships in Brazil - Post-Doctorate
FAPESP's process: 13/12162-5 - Balancing selection and population differentiation on HLA genes
Grantee:Débora Yoshihara Caldeira Brandt
Support type: Scholarships abroad - Research Internship - Master's degree
FAPESP's process: 12/18010-0 - Balancing selection in the human genome: detection, causes and consequences
Grantee:Diogo Meyer
Support type: Regular Research Grants