1
|
Ko S, German CA, Jensen A, Shen J, Wang A, Mehrotra DV, Sun YV, Sinsheimer JS, Zhou H, Zhou JJ. GWAS of longitudinal trajectories at biobank scale. Am J Hum Genet 2022; 109:433-445. [PMID: 35196515 PMCID: PMC8948167 DOI: 10.1016/j.ajhg.2022.01.018] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2021] [Accepted: 01/25/2022] [Indexed: 12/12/2022] Open
Abstract
Biobanks linked to massive, longitudinal electronic health record (EHR) data make numerous new genetic research questions feasible. One among these is the study of biomarker trajectories. For example, high blood pressure measurements over visits strongly predict stroke onset, and consistently high fasting glucose and Hb1Ac levels define diabetes. Recent research reveals that not only the mean level of biomarker trajectories but also their fluctuations, or within-subject (WS) variability, are risk factors for many diseases. Glycemic variation, for instance, is recently considered an important clinical metric in diabetes management. It is crucial to identify the genetic factors that shift the mean or alter the WS variability of a biomarker trajectory. Compared to traditional cross-sectional studies, trajectory analysis utilizes more data points and captures a complete picture of the impact of time-varying factors, including medication history and lifestyle. Currently, there are no efficient tools for genome-wide association studies (GWASs) of biomarker trajectories at the biobank scale, even for just mean effects. We propose TrajGWAS, a linear mixed effect model-based method for testing genetic effects that shift the mean or alter the WS variability of a biomarker trajectory. It is scalable to biobank data with 100,000 to 1,000,000 individuals and many longitudinal measurements and robust to distributional assumptions. Simulation studies corroborate that TrajGWAS controls the type I error rate and is powerful. Analysis of eleven biomarkers measured longitudinally and extracted from UK Biobank primary care data for more than 150,000 participants with 1,800,000 observations reveals loci that significantly alter the mean or WS variability.
Collapse
Affiliation(s)
- Seyoon Ko
- Department of Computational Medicine, University of California, Los Angeles, Los Angeles, CA 90095, USA,Department of Biostatistics, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Christopher A. German
- Department of Biostatistics, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Aubrey Jensen
- Department of Biostatistics, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Judong Shen
- Biostatistics and Research Decision Sciences, Merck & Co., Inc., Kenilworth, NJ 07033, USA
| | - Anran Wang
- Biostatistics and Research Decision Sciences, Merck & Co., Inc., Kenilworth, NJ 07033, USA
| | - Devan V. Mehrotra
- Biostatistics and Research Decision Sciences, Merck & Co., Inc., Kenilworth, NJ 07033, USA
| | - Yan V. Sun
- Department of Epidemiology, Emory University, Atlanta, GA 30322, USA
| | - Janet S. Sinsheimer
- Department of Computational Medicine, University of California, Los Angeles, Los Angeles, CA 90095, USA,Department of Biostatistics, University of California, Los Angeles, Los Angeles, CA 90095, USA,Department of Human Genetics, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Hua Zhou
- Department of Computational Medicine, University of California, Los Angeles, Los Angeles, CA 90095, USA; Department of Biostatistics, University of California, Los Angeles, Los Angeles, CA 90095, USA.
| | - Jin J. Zhou
- Department of Biostatistics, University of California, Los Angeles, Los Angeles, CA 90095, USA,Department of Medicine, University of California, Los Angeles, Los Angeles, CA 90095, USA,Department of Epidemiology and Biostatistics, University of Arizona, Tucson, AZ 85721, USA,Corresponding author
| |
Collapse
|