A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies

Xihao Li,Han Chen,Margaret Sunitha Selvaraj,Eric Van Buren,Hufeng Zhou,Yuxuan Wang,Ryan Sun,Zachary R McCaw,Zhi Yu,Donna K Arnett,Joshua C Bis,John Blangero,Eric Boerwinkle,Donald W Bowden,Jennifer A Brody,Brian E Cade,April P Carson,Jenna C Carlson,Nathalie Chami,Yii-Der Ida Chen,Joanne E Curran,Paul S de Vries,Myriam Fornage,Nora Franceschini,Barry I Freedman,Charles Gu,Nancy L Heard-Costa,Jiang He,Lifang Hou,Yi-Jen Hung,Marguerite R Irvin,Robert C Kaplan,Sharon L R Kardia,Tanika Kelly,Iain Konigsberg,Charles Kooperberg,Brian G Kral,Changwei Li,Ruth J F Loos,Michael C Mahaney,Lisa W Martin,Rasika A Mathias,Ryan L Minster,Braxton D Mitchell,May E Montasser,Alanna C Morrison,Nicholette D Palmer,Patricia A Peyser,Bruce M Psaty,Laura M Raffield,Susan Redline,Alexander P Reiner,Stephen S Rich,Colleen M Sitlani,Jennifer A Smith,Kent D Taylor,Hemant Tiwari,Ramachandran S Vasan,Zhe Wang,Lisa R Yanek,Bing Yu,NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium,Kenneth M Rice,Jerome I Rotter,Gina M Peloso,Pradeep Natarajan,Zilin Li,Zhonghua Liu,Xihong Lin

DOI: https://doi.org/10.1101/2023.10.30.564764

2023-11-02

bioRxiv

Abstract:Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.