A time-varying group sparse additive model for genome-wide association studies of dynamic complex traits

MOTIVATION Despite the widespread popularity of genome-wide association studies (GWAS) for genetic mapping of complex traits, most existing GWAS methodologies are still limited to the use of static phenotypes measured at a single time point. In this work, we propose a new method for association mapping that considers dynamic phenotypes measured at a sequence of time points. Our approach relies on the use of Time-Varying Group Sparse Additive Models (TV-GroupSpAM) for high-dimensional, functional regression. RESULTS This new model detects a sparse set of genomic loci that are associated with trait dynamics, and demonstrates increased statistical power over existing methods. We evaluate our method via experiments on synthetic data and perform a proof-of-concept analysis for detecting single nucleotide polymorphisms associated with two phenotypes used to assess asthma severity: forced vital capacity, a sensitive measure of airway obstruction and bronchodilator response, which measures lung response to bronchodilator drugs. AVAILABILITY AND IMPLEMENTATION Source code for TV-GroupSpAM freely available for download at http://www.cs.cmu.edu/~mmarchet/projects/tv_group_spam, implemented in MATLAB. CONTACT epxing@cs.cmu.edu SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

[1]  Guifang Fu,et al.  Dynamic semiparametric Bayesian models for genetic mapping of complex trait with irregular longitudinal data , 2013, Statistics in medicine.

[2]  Nicholas G Martin,et al.  Robust estimation of experimentwise P values applied to a genome scan of multiple asthma traits identifies a new region of significant linkage on chromosome 20q13. , 2005, American journal of human genetics.

[3]  H Yasue,et al.  Variants of NOS1, NOS2, and NOS3 genes in asthmatics. , 2000, Biochemical and biophysical research communications.

[4]  H. Timmerman,et al.  Mast cell tryptase and asthma , 1997, Mediators of inflammation.

[5]  D. Conrad,et al.  A potential new target for asthma therapy: A Disintegrin and Metalloprotease 10 (ADAM10) involvement in murine experimental asthma , 2011, Allergy.

[6]  C. Ober,et al.  Asthma genetics 2006: the long and winding road to gene discovery , 2006, Genes and Immunity.

[7]  John F Alcorn,et al.  A tale of two cytokines: IL-17 and IL-22 in asthma and infection , 2014, Expert review of respiratory medicine.

[8]  Byeong U. Park,et al.  Time-Varying Additive Models for Longitudinal Data , 2013 .

[9]  R. Tibshirani,et al.  Generalized Additive Models , 1991 .

[10]  P. Mahesh,et al.  An understanding of the genetic basis of asthma , 2011, The Indian journal of medical research.

[11]  Shuiwang Ji,et al.  SLEP: Sparse Learning with Efficient Projections , 2011 .

[12]  Runze Li,et al.  BAYESIAN GROUP LASSO FOR NONPARAMETRIC VARYING-COEFFICIENT MODELS WITH APPLICATION TO FUNCTIONAL GENOME-WIDE ASSOCIATION STUDIES. , 2015, The annals of applied statistics.

[13]  Andrew P Morris,et al.  Basic statistical analysis in genetic case-control studies , 2011, Nature Protocols.

[14]  Guifang Fu,et al.  The Bayesian lasso for genome-wide association studies , 2011, Bioinform..

[15]  Carole Ober,et al.  The genetics of asthma and allergic disease: a 21st century perspective , 2011, Immunological reviews.

[16]  Se-Min Park,et al.  Association analysis of thromboxane A synthase 1 gene polymorphisms with aspirin intolerance in asthmatic patients. , 2011, Pharmacogenomics.

[17]  Chunli Feng,et al.  Prostaglandin E2 deficiency uncovers a dominant role for thromboxane A2 in house dust mite-induced allergic pulmonary inflammation , 2012, Proceedings of the National Academy of Sciences.

[18]  Todd M Umstead,et al.  Essential regulation of lung surfactant homeostasis by the orphan G protein-coupled receptor GPR116. , 2013, Cell reports.

[19]  Y. Bossé,et al.  Updates on the COPD gene list , 2012, International journal of chronic obstructive pulmonary disease.

[20]  Rongling Wu,et al.  Nonparametric Functional Mapping of Quantitative Trait Loci , 2009, Biometrics.

[21]  Heather Eliassen,et al.  The Childhood Asthma Management Program (CAMP): design, rationale, and methods. Childhood Asthma Management Program Research Group. , 1999, Controlled clinical trials.

[22]  Shilpy Sharma,et al.  TGFbeta1 haplotypes and asthma in Indian populations. , 2005, The Journal of allergy and clinical immunology.

[23]  M. Yuan,et al.  Model selection and estimation in regression with grouped variables , 2006 .

[24]  Xi Chen,et al.  Group Sparse Additive Models , 2012, ICML.

[25]  D. Postma,et al.  Superoxide dismutases, lung function and bronchial responsiveness in a general population , 2009, European Respiratory Journal.

[26]  Steuart Rorke,et al.  Association of the ADAM33 gene with asthma and bronchial hyperresponsiveness , 2002, Nature.

[27]  Jyotsna Batra,et al.  Genetic contribution of chemokine receptor 2 (CCR2) polymorphisms towards increased serum total IgE levels in Indian asthmatics. , 2009, Genomics.

[28]  Jingyong Zhao,et al.  Death receptor-6 regulates the development of pulmonary eosinophilia and airway inflammation in a mouse model of asthma. , 2006, Immunology letters.

[29]  Runze Li,et al.  A dynamic model for genome-wide association studies , 2011, Human Genetics.

[30]  Trevor J. Hastie,et al.  Genome-wide association analysis by lasso penalized logistic regression , 2009, Bioinform..

[31]  Zitong Li,et al.  A Bayesian Nonparametric Approach for Mapping Dynamic Quantitative Traits , 2013, Genetics.

[32]  Wen-Miin Liang,et al.  Genetic variants in PLCB4/PLCB1 as susceptibility loci for coronary artery aneurysm formation in Kawasaki disease in Han Chinese in Taiwan , 2015, Scientific Reports.

[33]  Eleazar Eskin,et al.  Genome‐Wide Association Mapping With Longitudinal Data , 2012, Genetic epidemiology.

[34]  Qing-ling Huang,et al.  PLC/CAMK IV–NF-κB involved in the receptor for advanced glycation end products mediated signaling pathway in human endothelial cells , 2010, Molecular and Cellular Endocrinology.