Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Correspondence
  • Published:

Reconciling S-LDSC and LDAK functional enrichment estimates

Nature Genetics volume 51, pages 1202–1204 (2019)Cite this article

Subjects

To the Editor — Recent work has highlighted the importance of accounting for linkage disequilibrium (LD)-dependent genetic architectures in analyses of heritability1,2,3,4,5. Two models incorporating LD-dependent architectures have been proposed for analyses of functional enrichment: the baseline-LD model4 used by stratified LD-score regression4,6 (S-LDSC) and the LDAK model1,3. Although both models include LD-dependent effects, they produce very different estimates of functional enrichment (for example, 9.35x ± 0.80 in ref. 4 and 1.34x ± 0.26 in ref. 3 for conserved regions), thus leading to different interpretations of the functional architecture of complex traits. To reconcile these findings, we performed a comprehensive set of formal model comparisons and empirical analyses. Each of these analyses supports the higher functional enrichment estimates of S-LDSC with the baseline-LD model; each paragraph below is detailed in a corresponding section of the Supplementary Note (and detailed analyses can be found in Supplementary Tables 1–10 and Supplementary Figs. 1–23).

We defined six heritability models, including the infinitesimal model that ref. 3 called the ‘GCTA model’7; the baseline-LD model4, combining functional annotations6 with LD-dependent and minor allele frequency (MAF)-dependent architectures; and the LDAK model3, combining LD-dependent and MAF-dependent architectures. Notably, the baseline-LD and LDAK models use very different LD-dependent architectures. For comparison purposes, we also defined the ‘α-model’1, comprising only the MAF-dependent part of the LDAK model; the ‘Gazal-LD model’, comprising only the LD-dependent and MAF-dependent parts of the baseline-LD model (analogous to the LDAK model); and the ‘baseline-LD + LDAK’ model, which flexibly combines baseline-LD and LDAK model annotations.

This is a preview of subscription content, access via your institution

Relevant articles

Open Access articles citing this article.

Access options

Buy this article

  • Purchase on Springer Link
  • Instant access to full article PDF
Buy now

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Likelihood comparison of different models of per-SNP heritability.
Fig. 2: Comparison of functional enrichment estimates in analyses of UK Biobank traits.

Data availability

Baseline-LD model version 1.1 can be found at https://data.broadinstitute.org/alkesgroup/LDSCORE/1000G_Phase3_baselineLD_v1.1_ldscores.tgz. UK Biobank association statistics, computed with BOLT-LMM v2.3, are available at http://data.broadinstitute.org/alkesgroup/UKBB/.

Code availability

LDSC software is available at https://github.com/bulik/ldsc. LDAK version 5 is available at http://dougspeed.com/downloads/.

References

  1. Speed, D., Hemani, G., Johnson, M. R. & Balding, D. J. Am. J. Hum. Genet. 91, 1011–1021 (2012).

    Article  CAS  Google Scholar 

  2. Yang, J. et al. Nat. Genet. 47, 1114–1120 (2015).

    Article  CAS  Google Scholar 

  3. Speed, D., Cai, N., Johnson, M. R., Nejentsev, S. & Balding, D. J. Nat. Genet. 49, 986–992 (2017).

    Article  CAS  Google Scholar 

  4. Gazal, S. et al. Nat. Genet. 49, 1421–1427 (2017).

    Article  CAS  Google Scholar 

  5. Speed, D. & Balding, D. J. Nat. Genet. 51, 277–284 (2019).

    Article  CAS  Google Scholar 

  6. Finucane, H. K. et al. Nat. Genet. 47, 1228–1235 (2015).

    Article  CAS  Google Scholar 

  7. Yang, J. et al. Nat. Genet. 42, 565–569 (2010).

    Article  CAS  Google Scholar 

  8. Gazal, S. et al. Nat. Genet. 50, 1600–1607 (2018).

    Article  CAS  Google Scholar 

  9. Farh, K. K.-H. et al. Nature 518, 337–343 (2015).

    Article  CAS  Google Scholar 

  10. Huang, H. et al. Nature 547, 173–178 (2017).

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We are grateful to P.-R. Loh for assistance with UK Biobank data and to L. O’Connor, D. Speed and D. Balding for helpful discussions. This research was conducted by using the UK Biobank Resource under Application 16549. A.L.P. is funded by NIH grants U01 HG009379, R01 MH101244 and R01 MH107649. S.G. is funded by NIH K99 HG010160-01. H.K.F. is funded by Eric and Wendy Schmidt and NIH DP5-OD024582. Computational analyses were performed on the Orchestra High-Performance Compute Cluster at Harvard Medical School.

Author information

Author notes

  1. These authors jointly supervised this work: Hilary K. Finucane, Alkes L. Price.

Authors and Affiliations

  1. Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA

    Steven Gazal & Alkes L. Price

  2. Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA, USA

    Steven Gazal, Carla Marquez-Luna, Hilary K. Finucane & Alkes L. Price

  3. Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA

    Carla Marquez-Luna & Alkes L. Price

Authors
  1. Steven Gazal
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Carla Marquez-Luna
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Hilary K. Finucane
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Alkes L. Price
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

S.G., H.K.F. and A.L.P. designed experiments. S.G. performed experiments. S.G., C.M.L. and H.K.F. analyzed data. S.G., H.K.F. and A.L.P., with assistance from C.M.L., wrote the manuscript.

Corresponding authors

Correspondence to Steven Gazal, Hilary K. Finucane or Alkes L. Price.

Ethics declarations

Competing interests

The authors declare no competing interests.

Supplementary information

Supplementary Information

Supplementary Notes, Supplementary Figs. 1–23 and Supplementary Tables 1 and 5

Reporting Summary

Supplementary Tables

Supplementary Tables 2–4 and 6–10

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gazal, S., Marquez-Luna, C., Finucane, H.K. et al. Reconciling S-LDSC and LDAK functional enrichment estimates. Nat Genet 51, 1202–1204 (2019). https://doi.org/10.1038/s41588-019-0464-1

Download citation

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41588-019-0464-1

This article is cited by

Search

Advanced search

Quick links

Nature Briefing: Translational Research

Sign up for the Nature Briefing: Translational Research newsletter — top stories in biotechnology, drug discovery and pharma.

Get what matters in translational research, free to your inbox weekly. Sign up for Nature Briefing: Translational Research