Home LiteratureArticle Details
PMID: 24990607 Published · ppublish English Journal Article Research Support, N.I.H., Extramural Research Support, Non-U.S. Gov't

Fast and accurate imputation of summary statistics enhances evidence of functional enrichment.

Bioinformatics (Oxford, England) ·Vol. 30 ·No. 20 ·2014-10-15 ·Pages 2906-14

Pasaniuc B, Zaitlen N, Shi H, Bhatia G, Gusev A, Pickrell J, Hirschhorn J, Strachan DP, Patterson N, Price AL

Abstract

Imputation using external reference panels (e.g. 1000 Genomes) is a widely used approach for increasing power in genome-wide association studies and meta-analysis. Existing hidden Markov models (HMM)-based imputation approaches require individual-level genotypes. Here, we develop a new method for Gaussian imputation from summary association statistics, a type of data that is becoming widely available. In simulations using 1000 Genomes (1000G) data, this method recovers 84% (54%) of the effective sample size for common (>5%) and low-frequency (1-5%) variants [increasing to 87% (60%) when summary linkage disequilibrium information is available from target samples] versus the gold standard of 89% (67%) for HMM-based imputation, which cannot be applied to summary statistics. Our approach accounts for the limited sample size of the reference panel, a crucial step to eliminate false-positive associations, and it is computationally very fast. As an empirical demonstration, we apply our method to seven case-control phenotypes from the Wellcome Trust Case Control Consortium (WTCCC) data and a study of height in the British 1958 birth cohort (1958BC). Gaussian imputation from summary statistics recovers 95% (105%) of the effective sample size (as quantified by the ratio of [Formula: see text] association statistics) compared with HMM-based imputation from individual-level genotypes at the 227 (176) published single nucleotide polymorphisms (SNPs) in the WTCCC (1958BC height) data. In addition, for publicly available summary statistics from large meta-analyses of four lipid traits, we publicly release imputed summary statistics at 1000G SNPs, which could not have been obtained using previously published methods, and demonstrate their accuracy by masking subsets of the data. We show that 1000G imputation using our approach increases the magnitude and statistical evidence of enrichment at genic versus non-genic loci for these traits, as compared with an analysis without 1000G imputation. Thus, imputation of summary statistics will be a valuable tool in future functional enrichment analyses. Publicly available software package available at http://bogdan.bioinformatics.ucla.edu/software/. [email protected] or [email protected] Supplementary materials are available at Bioinformatics online.

MeSH Terms
Algorithms Biostatistics/methods Case-Control Studies Cohort Studies Genome-Wide Association Study/methods Genotype Humans Linkage Disequilibrium Phenotype Polymorphism, Single Nucleotide Software Time Factors
Authors & Affiliations
10 authors, click to expand affiliations / ORCID
Pasaniuc Bogdan
Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK.
Zaitlen Noah
Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK.
Shi Huwenbo
Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK.
Bhatia Gaurav
Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Har
Gusev Alexander
Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Har
Pickrell Joseph
Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK.
Hirschhorn Joel
Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK.
Strachan David P
Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK.
Patterson Nick
Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK.
Price Alkes L
Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Harvard School of Public Health, Boston, 02115, Departments of Epidemiology and Biostatistics, Harvard School of Public Health, Boston, MA, 02115, Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, 02142, Department of Genetics Harvard Medical School, Boston, MA, 02115 and Division of Population Health Sciences and Education, St George's, University of London, UK Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, 90024, Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, 90024, Department of Medicine, Lung Biology Center, University of California San Francisco, San Francisco, 94143, Program in Genetic Epidemiology and Statistical Genetics, Har
References (39)
39 references, click to expand
  1. An integrated map of genetic variation from 1,092 human genomes.
    Nature. 2012 Nov 1;491(7422):56-65 PMID: 23128226
  2. DIST: direct imputation of summary statistics for unmeasured SNPs.
    Bioinformatics. 2013 Nov 15;29(22):2925-7 PMID: 23990413
  3. Asking for more.
    Nat Genet. 2012 Jun 27;44(7):733 PMID: 22735581
  4. Fast and accurate genotype imputation in genome-wide association studies through pre-phasing.
    Nat Genet. 2012 Jul 22;44(8):955-9 PMID: 22820512
  5. Phasing of many thousands of genotyped samples.
    Am J Hum Genet. 2012 Aug 10;91(2):238-51 PMID: 22883141
  6. All SNPs are not created equal: genome-wide association studies reveal a consistent pattern of enrichment among functionally annotated SNPs.
    PLoS Genet. 2013 Apr;9(4):e1003449 PMID: 23637621
  7. MaCH: using sequence and genotype data to estimate haplotypes and unobserved genotypes.
    Genet Epidemiol. 2010 Dec;34(8):816-34 PMID: 21058334
  8. Low-coverage sequencing: implications for design of complex trait association studies.
    Genome Res. 2011 Jun;21(6):940-51 PMID: 21460063
  9. ProbABEL package for genome-wide association analysis of imputed data.
    BMC Bioinformatics. 2010 Mar 16;11:134 PMID: 20233392
  10. Genotype imputation for genome-wide association studies.
    Nat Rev Genet. 2010 Jul;11(7):499-511 PMID: 20517342
  11. Hundreds of variants clustered in genomic loci and biological pathways affect human height.
    Nature. 2010 Oct 14;467(7317):832-8 PMID: 20881960
  12. Lifecourse influences on health among British adults: effects of region of residence in childhood and adulthood.
    Int J Epidemiol. 2007 Jun;36(3):522-31 PMID: 17255346
  13. Rapid and accurate multiple testing correction and power estimation for millions of correlated markers.
    PLoS Genet. 2009 Apr;5(4):e1000456 PMID: 19381255
  14. BLUP genotype imputation for case-control association testing with related individuals and missing data.
    J Comput Biol. 2012 Jun;19(6):756-65 PMID: 22697245
  15. Efficiency and power as a function of sequence coverage, SNP array density, and imputation.
    PLoS Comput Biol. 2012;8(7):e1002604 PMID: 22807667
  16. Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture.
    Nat Genet. 2012 Apr 15;44(5):491-501 PMID: 22504420
  17. Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes.
    Nat Genet. 2009 Jun;41(6):703-7 PMID: 19430480
  18. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering.
    Am J Hum Genet. 2007 Nov;81(5):1084-97 PMID: 17924348
  19. General framework for meta-analysis of rare variants in sequencing association studies.
    Am J Hum Genet. 2013 Jul 11;93(1):42-53 PMID: 23768515
  20. Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes.
    Nat Genet. 2012 Sep;44(9):981-90 PMID: 22885922
  21. Genome-wide meta-analyses identify three loci associated with primary biliary cirrhosis.
    Nat Genet. 2010 Aug;42(8):658-60 PMID: 20639880
  22. Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits.
    Nat Genet. 2012 Mar 18;44(4):369-75, S1-3 PMID: 22426310
  23. A large-scale, consortium-based genomewide association study of asthma.
    N Engl J Med. 2010 Sep 23;363(13):1211-1221 PMID: 20860503
  24. Potential etiologic and functional implications of genome-wide association loci for human diseases and traits.
    Proc Natl Acad Sci U S A. 2009 Jun 9;106(23):9362-7 PMID: 19474294
  25. USING LINEAR PREDICTORS TO IMPUTE ALLELE FREQUENCIES FROM SUMMARY OR POOLED GENOTYPE DATA.
    Ann Appl Stat. 2010 Sep;4(3):1158-1182 PMID: 21479081
  26. Genomic inflation factors under polygenic inheritance.
    Eur J Hum Genet. 2011 Jul;19(7):807-12 PMID: 21407268
  27. Genotype and SNP calling from next-generation sequencing data.
    Nat Rev Genet. 2011 Jun;12(6):443-51 PMID: 21587300
  28. HAPGEN2: simulation of multiple disease SNPs.
    Bioinformatics. 2011 Aug 15;27(16):2304-5 PMID: 21653516
  29. Genotype imputation via matrix completion.
    Genome Res. 2013 Mar;23(3):509-18 PMID: 23233546
  30. Genomic privacy and limits of individual detection in a pool.
    Nat Genet. 2009 Sep;41(9):965-7 PMID: 19701190
  31. A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals.
    Am J Hum Genet. 2009 Feb;84(2):210-23 PMID: 19200528
  32. Resolving individuals contributing trace amounts of DNA to highly complex mixtures using high-density SNP genotyping microarrays.
    PLoS Genet. 2008 Aug 29;4(8):e1000167 PMID: 18769715
  33. Biological, clinical and population relevance of 95 loci for blood lipids.
    Nature. 2010 Aug 5;466(7307):707-13 PMID: 20686565
  34. Meta-analysis of gene-level associations for rare variants based on single-variant statistics.
    Am J Hum Genet. 2013 Aug 8;93(2):236-48 PMID: 23891470
  35. Genomic control for association studies.
    Biometrics. 1999 Dec;55(4):997-1004 PMID: 11315092
  36. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls.
    Nature. 2007 Jun 7;447(7145):661-78 PMID: 17554300
  37. Leveraging genetic variability across populations for the identification of causal variants.
    Am J Hum Genet. 2010 Jan;86(1):23-33 PMID: 20085711
  38. So many correlated tests, so little time! Rapid adjustment of P values for multiple correlated tests.
    Am J Hum Genet. 2007 Dec;81(6):1158-68 PMID: 17966093
  39. Extremely low-coverage sequencing and imputation increases power for genome-wide association studies.
    Nat Genet. 2012 May 20;44(6):631-5 PMID: 22610117
Article Info
Journal
Bioinformatics (Oxford, England)
Abbr.
Bioinformatics
ISSN
1367-4811
DOI
10.1093/bioinformatics/btu416
Published
2014-10-15
Epub
2014-00-01
Pages
2906-14
Language
English
Region
England
NLM ID
9808944
PMCID
PMC4184260
Subset
IM
Grants
Medical Research Council · G1001799 · United Kingdom
NIGMS NIH HHS · F32 GM106584 · United States
NIEHS NIH HHS · T32 ES007142 · United States
Wellcome Trust · 068545/Z/02 · United Kingdom
NHGRI NIH HHS · R01 HG006399 · United States
NCI NIH HHS · R03 CA162200 · United States
Medical Research Council · G0000934 · United Kingdom
NIGMS NIH HHS · R01 GM053275 · United States
External Links
PubMed source DOI: 10.1093/bioinformatics/btu416 PMC Full Text View journal details
Analysis Services
Analysis Services

Contact

No. 2 Wenbo Road, Zhangqiu District, Jinan, Shandong

Qilu Normal University · Genelibs Bioinformatics Lab

750 Shunhua Rd, Jinan

2F, Bldg F, University Science Park

Tel: 0531-88819269

WeChat Official Account

Follow our WeChat subscription account for real-time updates and the latest in medical and biological research.


Business Email

E-mail: [email protected]

Copyright © 2026.Genelibs Jinan Honghui Biotechnology Co., Ltd. 鲁ICP备20007159号-2