Abstract
IMPORTANCE: Self-reported trauma exposure has consistently been found to be a risk factor for major depressive disorder (MDD), and several studies have reported interactions with genetic liability. To date, most studies have examined gene-environment interactions with trauma exposure using genome-wide variants (single-nucleotide variations [SNVs]) or polygenic scores, both typically capturing less than 3% of phenotypic risk variance. OBJECTIVE: To reexamine genome-by-trauma interaction associations using genetic measures using all available genotyped data and thus, maximizing accounted variance. DESIGN, SETTING, AND PARTICIPANTS: The UK Biobank study was conducted from April 2007 to May 1, 2016 (follow-up mental health questionnaire). The current study used available cross-sectional genomic and trauma exposure data from UK Biobank. Participants who completed the mental health questionnaire and had available genetic, trauma experience, depressive symptoms, and/or neuroticism information were included. Data were analyzed from April 1 to August 30, 2021. EXPOSURES: Trauma and genome-by-trauma exposure interactions. MAIN OUTCOMES AND MEASURES: Measures of self-reported depression, neuroticism, and trauma exposure with whole-genome SNV data are available from the UK Biobank study. Here, a mixed-model statistical approach using genetic, trauma exposure, and genome-by-trauma exposure interaction similarity matrices was used to explore sources of variation in depression and neuroticism. RESULTS: Analyses were conducted on 148 129 participants (mean [SD] age, 56 [7] years) of which 76 995 were female (52.0%). The study approach estimated the heritability (SE) of MDD to be approximately 0.160 (0.016). Subtypes of self-reported trauma exposure (catastrophic, adult, childhood, and full trauma) accounted for a significant proportion of the variance of MDD, with heritability (SE) ranging from 0.056 (0.013) to 0.176 (0.025). The proportion of MDD risk variance accounted for by significant genome-by-trauma interaction revealed estimates (SD) ranging from 0.074 (0.006) to 0.201 (0.009). Results from sex-specific analyses found genome-by-trauma interaction variance estimates approximately 5-fold greater for MDD in male participants (0.441 [0.018]) than in female participants (0.086 [0.009]). CONCLUSIONS AND RELEVANCE: This cross-sectional study used an approach combining all genome-wide SNV data when exploring genome-by-trauma interactions in individuals with MDD; findings suggest that such interactions were associated with depression manifestation. Genome-by-trauma interaction accounts for greater trait variance in male individuals, which points to potential differences in depression etiology between the sexes. The methodology used in this study can be extrapolated to other environmental factors to identify modifiable risk environments and at-risk groups to target with interventions.