Correcting systematic inflation in genetic association tests that consider interaction effects: application to a genome-wide association study of posttraumatic stress disorder.

IMPORTANCE Genetic association studies of psychiatric outcomes often consider interactions with environmental exposures and, in particular, apply tests that jointly consider gene and gene-environment interaction effects for analysis. Using a genome-wide association study (GWAS) of posttraumatic stress disorder (PTSD), we report that heteroscedasticity (defined as variability in outcome that differs by the value of the environmental exposure) can invalidate traditional joint tests of gene and gene-environment interaction. OBJECTIVES To identify the cause of bias in traditional joint tests of gene and gene-environment interaction in a PTSD GWAS and determine whether proposed robust joint tests are insensitive to this problem. DESIGN, SETTING, AND PARTICIPANTS The PTSD GWAS data set consisted of 3359 individuals (978 men and 2381 women) from the Grady Trauma Project (GTP), a cohort study from Atlanta, Georgia. The GTP performed genome-wide genotyping of participants and collected environmental exposures using the Childhood Trauma Questionnaire and Trauma Experiences Inventory. MAIN OUTCOMES AND MEASURES We performed joint interaction testing of the Beck Depression Inventory and modified PTSD Symptom Scale in the GTP GWAS. We assessed systematic bias in our interaction analyses using quantile-quantile plots and genome-wide inflation factors. RESULTS Application of the traditional joint interaction test to the GTP GWAS yielded systematic inflation across different outcomes and environmental exposures (inflation-factor estimates ranging from 1.07 to 1.21), whereas application of the robust joint test to the same data set yielded no such inflation (inflation-factor estimates ranging from 1.01 to 1.02). Simulated data further revealed that the robust joint test is valid in different heteroscedasticity models, whereas the traditional joint test is invalid. The robust joint test also has power similar to the traditional joint test when heteroscedasticity is not an issue. CONCLUSIONS AND RELEVANCE We believe the robust joint test should be used in candidate-gene studies and GWASs of psychiatric outcomes that consider environmental interactions. To make the procedure useful for applied investigators, we created a software tool that can be called from the popular PLINK package for analysis.

[1]  Kerry J Ressler,et al.  How the neurocircuitry and genetics of fear inhibition may inform our understanding of PTSD. , 2010, The American journal of psychiatry.

[2]  B. Bradley,et al.  Impaired fear inhibition is a biomarker of PTSD but not depression , 2010, Depression and anxiety.

[3]  S. Sen,et al.  Gene × environment interaction models in psychiatric genetics. , 2012, Current topics in behavioral neurosciences.

[4]  Stacey J Winham,et al.  Gene-environment interactions in genome-wide association studies: current approaches and new directions. , 2013, Journal of child psychology and psychiatry, and allied disciplines.

[5]  H. H. Ku,et al.  Contributions to Probability and Statistics, Essays in Honor of Harold Hotelling. , 1961 .

[6]  M. Stein,et al.  Genetic and environmental influences on trauma exposure and posttraumatic stress disorder symptoms: a twin study. , 2002, The American journal of psychiatry.

[7]  John P.A. Ioannidis,et al.  Non-Replication and Inconsistency in the Genome-Wide Association Setting , 2007, Human Heredity.

[8]  N. Breslau.  The epidemiology of posttraumatic stress disorder: what is the extent of the problem? , 2001, The Journal of clinical psychiatry.

[9]  C. A. Morgan,et al.  Fear-potentiated startle in posttraumatic stress disorder , 1995, Biological Psychiatry.

[10]  G. King,et al.  Unifying Political Methodology: The Likelihood Theory of Statistical Inference , 1989 .

[11]  W R True,et al.  A twin study of genetic and environmental contributions to liability for posttraumatic stress symptoms. , 1993, Archives of general psychiatry.

[12]  H. Levene.  Robust tests for equality of variances , 1961 .

[13]  Claude Bouchard,et al.  A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance , 2012, Nature Genetics.

[14]  N. Chatterjee,et al.  Powerful multilocus tests of genetic association in the presence of gene-gene and gene-environment interactions. , 2006, American journal of human genetics.

[15]  Tao Wang,et al.  A partial least‐square approach for modeling gene‐gene and gene‐environment interactions when multiple markers are genotyped , 2009, Genetic epidemiology.

[16]  M. Stein,et al.  Interaction of FKBP5 with Childhood Adversity on Risk for Post-Traumatic Stress Disorder , 2010, Neuropsychopharmacology.

[17]  Patrick M O'Malley,et al.  Associations between access to food stores and adolescent body mass index. , 2007, American journal of preventive medicine.

[18]  J. van Os,et al.  Gene-environment-wide interaction studies in psychiatry. , 2009, The American journal of psychiatry.

[19]  Jung-Ying Tzeng,et al.  Studying gene and gene-environment effects of uncommon and common variants on continuous traits: a marker-set approach using gene-trait similarity regression. , 2011, American journal of human genetics.

[20]  A. Beck,et al.  An inventory for measuring depression. , 1961, Archives of general psychiatry.

[21]  Marilyn C Cornelis,et al.  Genetics of Post-Traumatic Stress Disorder: Review and Recommendations for Genome-Wide Association Studies , 2010, Current psychiatry reports.

[22]  M. J. van de Vijver,et al.  The pathology of familial breast cancer: predictive value of immunohistochemical markers estrogen receptor, progesterone receptor, HER-2, and p53 in patients with mutations in BRCA1 and BRCA2. , 2002, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[23]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[24]  W. Greene,et al.  计量经济分析 = Econometric analysis , 2009 .

[25]  B. Bradley,et al.  Childhood abuse is associated with increased startle reactivity in adulthood , 2009, Depression and anxiety.

[26]  Taylor J. Maxwell,et al.  APOE Modulates the Correlation Between Triglycerides, Cholesterol, and CHD Through Pleiotropy, and Gene-by-Gene Interactions , 2013, Genetics.

[27]  S. Cichon,et al.  Genomewide association studies: history, rationale, and prospects for psychiatric disorders. , 2009, The American journal of psychiatry.

[28]  M. Stein,et al.  Interactive effect of stressful life events and the serotonin transporter 5-HTTLPR genotype on posttraumatic stress disorder diagnosis in 2 independent populations. , 2009, Archives of general psychiatry.

[29]  D. Conti,et al.  Nicotinic acetylcholine receptor β2 subunit gene implicated in a systems-based candidate gene study of smoking cessation , 2008, Human molecular genetics.

[30]  Yurii S. Aulchenko,et al.  ProbABEL package for genome-wide association analysis of imputed data , 2010, BMC Bioinformatics.

[31]  A. Fitzpatrick,et al.  Modifications to the Healthy Eating Index and its ability to predict obesity: the Multi-Ethnic Study of Atherosclerosis. , 2008, The American journal of clinical nutrition.

[32]  U. Nöthlings,et al.  Genome-wide investigation of gene–environment interactions in colorectal cancer , 2012, Human Genetics.

[33]  Joseph E LeDoux,et al.  Response Variation following Trauma: A Translational Neuroscience Approach to Understanding PTSD , 2007, Neuron.

[34]  Peter Kraft,et al.  Genome-Wide Meta-Analysis of Joint Tests for Genetic and Gene-Environment Interaction Effects , 2011, Human Heredity.

[35]  Nilanjan Chatterjee,et al.  Semiparametric maximum likelihood estimation exploiting gene-environment independence in case-control studies , 2005 .

[36]  Thomas Lumley,et al.  Behavior of QQ-Plots and Genomic Control in Studies of Gene-Environment Interaction , 2011, PloS one.

[37]  Beate Ritz,et al.  Genome-Wide Gene-Environment Study Identifies Glutamate Receptor Gene GRIN2A as a Parkinson's Disease Modifier Gene via Interaction with Coffee , 2011, PLoS genetics.

[38]  A. Amstadter,et al.  Gene-environment interaction in posttraumatic stress disorder: an update. , 2009, Journal of traumatic stress.

[39]  Peter Kraft,et al.  The Impact of Gene-Environment Dependence and Misclassification in Genetic Association Studies Incorporating Gene-Environment Interactions , 2009, Human Heredity.

[40]  S. Galea,et al.  The serotonin transporter genotype and social support and moderation of posttraumatic stress disorder and depression in hurricane-exposed adults. , 2007, The American journal of psychiatry.

[41]  P. J. Huber.  The behavior of maximum likelihood estimates under nonstandard conditions , 1967 .

[42]  Lisa M. Shin,et al.  Neurocircuitry Models of Posttraumatic Stress Disorder and Extinction: Human Neuroimaging Research—Past, Present, and Future , 2006, Biological Psychiatry.

[43]  B. Bradley,et al.  Trauma exposure and stress-related disorders in inner city primary care patients. , 2009, General hospital psychiatry.

[44]  K. Ressler,et al.  Posttraumatic stress disorder among African Americans in an inner city mental health clinic. , 2005, Psychiatric services.

[45]  Peter Kraft,et al.  Exploiting Gene-Environment Interaction to Detect Genetic Associations , 2007, Human Heredity.

[46]  S. Galea,et al.  Modification of the association between serotonin transporter genotype and risk of posttraumatic stress disorder in adults by county-level social environment. , 2009, American journal of epidemiology.

[47]  L. Wasserman,et al.  Genomic control, a new approach to genetic-based association studies. , 2001, Theoretical population biology.

[48]  E. Foa,et al.  Comparison of the PTSD Symptom Scale–Interview Version and the Clinician-Administered PTSD Scale , 2000, Journal of traumatic stress.

[49]  Jason H. Moore,et al.  Missing heritability and strategies for finding the underlying causes of complex disease , 2010, Nature Reviews Genetics.

[50]  Eric J Tchetgen Tchetgen,et al.  On the robustness of tests of genetic associations incorporating gene-environment interaction when the environmental exposure is misspecified. , 2011, Epidemiology.

[51]  K. Brady,et al.  Screening for PTSD in a Substance Abuse Sample: Psychometric Properties of a Modified Version of the PTSD Symptom Scale Self-Report , 1998, Journal of traumatic stress.

[52]  H. White.  A Heteroskedasticity-Consistent Covariance Matrix Estimator and a Direct Test for Heteroskedasticity , 1980 .

[53]  K. Ressler,et al.  Association of FKBP5 polymorphisms and childhood abuse with risk of posttraumatic stress disorder symptoms in adults. , 2008, JAMA.

[54]  Kai Wang,et al.  Genetic association tests in the presence of epistasis or gene‐environment interaction , 2008, Genetic epidemiology.

[55]  K. Ressler,et al.  Pain medication use among patients with posttraumatic stress disorder. , 2006, Psychosomatics.

[56]  H. Resnick,et al.  The Modified PTSD Symptom Scale: A brief self-report measure of posttraumatic stress disorder. , 1993 .