Browse/search for people

Publication - Professor Tom Gaunt

    PhenoSpD

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

    Citation

    Zheng, J, Richardson, T, Millard, L, Hemani, G, Elsworth, B, Raistrick, C, Vilhjalmsson, B, Neale, BM, Haycock, P, Smith, GD & Gaunt, T, 2018, ‘PhenoSpD: an integrated toolkit for phenotypic correlation estimation and multiple testing correction using GWAS summary statistics’. GigaScience, vol 7.

    Abstract

    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 1) use LD score regression to estimate phenotypic cor-relations using GWAS summary statistics; and 2) 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 1) it is possible to identify non-independence of phenotypes using samples with partial overlap, as overlap decreases the estimated phenotypic correlations will attenuate towards zero and multiple testing correction will be more stringent than in perfectly overlapping samples; 2) in contrast to LD score regression, metaCCA will provide approximate genetic corre-lations 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 estimat-ed 5618 pair-wise phenotypic correlations among 107 metabolites using GWAS summary statistics from Kettunen et al. and PhenoSpD suggested the equivalent of 33.5 independent tests for theses metabolites.
    Conclusion: 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.
    Availability: R code and documentation for PhenoSpD V1.0.0 is available online https://github.com/MRCIEU/PhenoSpD.

    Full details in the University publications repository