A unifying framework for rare variant association testing in family-based designs, including higher criticism approaches, SKATs, and burden tests

Julian Hecker; F William Townes; Priyadarshini Kachroo; Cecelia Laurie; Jessica Lasky-Su; John Ziniti; Michael H Cho; Scott T Weiss; Nan M Laird; Christoph Lange

doi:10.1093/bioinformatics/btaa1055

A unifying framework for rare variant association testing in family-based designs, including higher criticism approaches, SKATs, and burden tests

Bioinformatics. 2021 Apr 1;36(22-23):5432-5438. doi: 10.1093/bioinformatics/btaa1055.

Authors

Affiliations

¹ Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA 02115, USA.
² Department of Computer Science, Princeton University, Princeton, NJ 08540-5233, USA.
³ Department of Biostatistics, University of Washington, Seattle, WA 98195-1617, USA.
⁴ Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA.

Abstract

Motivation: Analysis of rare variants in family-based studies remains a challenge. Transmission-based approaches provide robustness against population stratification, but the evaluation of the significance of test statistics based on asymptotic theory can be imprecise. Also, power will depend heavily on the choice of the test statistic and on the underlying genetic architecture of the locus, which will be generally unknown.

Results: In our proposed framework, we utilize the FBAT haplotype algorithm to obtain the conditional offspring genotype distribution under the null hypothesis given the sufficient statistic. Based on this conditional offspring genotype distribution, the significance of virtually any association test statistic can be evaluated based on simulations or exact computations, without the need for asymptotic approximations. Besides standard linear burden-type statistics, this enables our approach to also evaluate other test statistics such as variance components statistics, higher criticism approaches, and maximum-single-variant-statistics, where asymptotic theory might be involved or does not provide accurate approximations for rare variant data. Based on these P-values, combined test statistics such as the aggregated Cauchy association test (ACAT) can also be utilized. In simulation studies, we show that our framework outperforms existing approaches for family-based studies in several scenarios. We also applied our methodology to a TOPMed whole-genome sequencing dataset with 897 asthmatic trios from Costa Rica.

Availability and implementation: FBAT software is available at https://sites.google.com/view/fbatwebpage. Simulation code is available at https://github.com/julianhecker/FBAT_rare_variant_test_simulations. Whole-genome sequencing data for 'NHLBI TOPMed: The Genetic Epidemiology of Asthma in Costa Rica' is available at https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000988.v4.p1.

Supplementary information: Supplementary data are available at Bioinformatics online.

Abstract

Grants and funding