Efficient variant set mixed model association tests for continuous and binary traits in large-scale whole genome sequencing studies
Jennifer E. Huffman,
Jennifer A Brody,
Stephanie M Gogarten,
Lawrence F Bielak,
Joshua C Bis,
Russell P. Bowler,
Brian E. Cade,
Michael H. Cho,
Joanne E. Curran,
Paul S. de Vries,
David C Glahn,
Andrew D. Johnson,
Joshua P. Lewis,
Rasika A. Mathias,
Braxton D Mitchell,
Jeffrey R. O’Connell,
Patricia A. Peyser,
Wendy S. Post,
Alex P Reiner,
Stephen S Rich,
Jerome I. Rotter,
Edwin K Silverman,
Jennifer A. Smith,
Ramachandran S Vasan,
James G Wilson,
Lisa R. Yanek,
NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium, TOPMed Hematology and Hemostasis Working Group,
Nicholas L. Smith,
Ingrid B Borecki,
L. Adrienne Cupples,
Cathy C Laurie,
Alanna C. Morrison,
Kenneth M Rice,
Posted 20 Aug 2018
bioRxiv DOI: 10.1101/395046 (published DOI: 10.1016/j.ajhg.2018.12.012)
Posted 20 Aug 2018
With advances in Whole Genome Sequencing (WGS) technology, more advanced statistical methods for testing genetic association with rare variants are being developed. Methods in which variants are grouped for analysis are also known as variant-set, gene-based, and aggregate unit tests. The burden test and Sequence Kernel Association Test (SKAT) are two widely used variant-set tests, which were originally developed for samples of unrelated individuals and later have been extended to family data with known pedigree structures. However, computationally-efficient and powerful variant-set tests are needed to make analyses tractable in large-scale WGS studies with complex study samples. In this paper, we propose the variant-Set Mixed Model Association Tests (SMMAT) for continuous and binary traits using the generalized linear mixed model framework. These tests can be applied to large-scale WGS studies involving samples with population structure and relatedness, such as in the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine (TOPMed) program. SMMAT tests share the same null model for different variant sets, and a virtue of this null model, which includes covariates only, is that it needs to be only fit once for all tests in each genome-wide analysis. Simulation studies show that all the proposed SMMAT tests correctly control type I error rates for both continuous and binary traits in the presence of population structure and relatedness. We also illustrate our tests in a real data example of analysis of plasma fibrinogen levels in the TOPMed program (n = 23,763), using the Analysis Commons, a cloud-based computing platform.
- Downloaded 802 times
- Download rankings, all-time:
- Site-wide: 19,062 out of 101,401
- In genetics: 1,182 out of 5,038
- Year to date:
- Site-wide: 44,988 out of 101,401
- Since beginning of last month:
- Site-wide: 71,278 out of 101,401
Downloads over time
Distribution of downloads per paper, site-wide
- 20 Oct 2020: Support for sorting preprints using Twitter activity has been removed, at least temporarily, until a new source of social media activity data becomes available.
- 18 Dec 2019: We're pleased to announce PanLingua, a new tool that enables you to search for machine-translated bioRxiv preprints using more than 100 different languages.
- 21 May 2019: PLOS Biology has published a community page about Rxivist.org and its design.
- 10 May 2019: The paper analyzing the Rxivist dataset has been published at eLife.
- 1 Mar 2019: We now have summary statistics about bioRxiv downloads and submissions.
- 8 Feb 2019: Data from Altmetric is now available on the Rxivist details page for every preprint. Look for the "donut" under the download metrics.
- 30 Jan 2019: preLights has featured the Rxivist preprint and written about our findings.
- 22 Jan 2019: Nature just published an article about Rxivist and our data.
- 13 Jan 2019: The Rxivist preprint is live!