PhenoSpD: an integrated toolkit for phenotypic correlation estimation and multiple testing correction using GWAS summary statistics

Research output: Contribution to journal/Conference contribution in journal/Contribution to newspaperJournal articleResearchpeer-review

Documents

  • giy090

    Final published version, 1 MB, PDF document

DOI

  • Jie Zheng, Bristol University
  • ,
  • Tom G. Richardson, Bristol University
  • ,
  • Louise A.C. Millard, Bristol University
  • ,
  • Gibran Hemani, Bristol University
  • ,
  • Benjamin L. Elsworth, Bristol University
  • ,
  • Christopher A. Raistrick, Bristol University
  • ,
  • Bjarni Vilhjalmsson
  • Benjamin M. Neale, Broad Institute, Massachusetts General Hospital and Harvard Medical School
  • ,
  • Philip C. Haycock, Bristol University
  • ,
  • George Davey Smith, Bristol University
  • ,
  • Tom R. Gaunt, Bristol University

Background: Identifying phenotypic correlations between complex traits and diseases can provide useful etiological insights. Restricted access to much individual-level phenotype data makes it difficult to estimate large-scale phenotypic correlation across the human phenome. Two state-of-the-art methods, metaCCA and LD score regression, provide an alternative approach to estimate phenotypic correlation using only genome-wide association study (GWAS) summary results.

Results: Here, we present an integrated R toolkit, PhenoSpD, to use LD score regression to estimate phenotypic correlations using GWAS summary statistics and to utilize the estimated phenotypic correlations to inform correction of multiple testing for complex human traits using the spectral decomposition of matrices (SpD). The simulations suggest that it is possible to identify nonindependence of phenotypes using samples with partial overlap; as overlap decreases, the estimated phenotypic correlations will attenuate toward zero and multiple testing correction will be more stringent than in perfectly overlapping samples. Also, in contrast to LD score regression, metaCCA will provide approximate genetic correlations rather than phenotypic correlation, which limits its application for multiple testing correction. In a case study, PhenoSpD using UK Biobank GWAS results suggested 399.6 independent tests among 487 human traits, which is close to the 352.4 independent tests estimated using true phenotypic correlation. We further applied PhenoSpD to an estimated 5,618 pair-wise phenotypic correlations among 107 metabolites using GWAS summary statistics from Kettunen's publication and PhenoSpD suggested the equivalent of 33.5 independent tests for these metabolites.

Conclusions: PhenoSpD extends the use of summary-level results, providing a simple and conservative way to reduce dimensionality for complex human traits using GWAS summary statistics. This is particularly valuable in the age of large-scale biobank and consortia studies, where GWAS results are much more accessible than individual-level data.

Original languageEnglish
JournalGigaScience
Volume7
Issue8
Pages (from-to)1-10
Number of pages10
ISSN2047-217X
DOIs
Publication statusPublished - 1 Aug 2018

See relations at Aarhus University Citationformats

Download statistics

No data available

ID: 135344243