Genome-wide study of major depression in 685,808 diverse individuals identifies 697 independent associations, infers causal neuronal subtypes and biological targets for novel pharmacotherapies.

In a genome-wide association study (GWAS) meta-analysis of 685,808 individuals with major depression (MD) and 4,364,225 controls from 29 countries and across diverse and admixed ancestries, we identify 697 independent associations at 636 loci, 293 of which are novel. Using fine-mapping and functional genomic tools, we find 308 high-confidence gene associations and enrichment of postsynaptic density and receptor clustering. Leveraging new single-cell gene expression data, we conducted a causal neural cell type enrichment analysis that implicates dysregulation of excitatory and inhibitory midbrain and forebrain neurons, peptidergic neurons, and medium spiny neurons in MD. Our findings are enriched for the targets of antidepressants and provide potential antidepressant repurposing opportunities (e.g., pregabalin and modafinil). Polygenic scores (PGS) trained using either European or multi-ancestry data significantly predicted MD status across all five diverse ancestries and explained a maximum of 5.8% of the variance in liability to MD in Europeans. These findings represent a major advance in our understanding of MD across global populations. MD GWAS reveals known and novel biological targets that may be used to target and develop pharmacotherapies addressing the considerable unmet need for effective treatment.

Genetic Architectures of Adolescent Depression Trajectories in 2 Longitudinal Population Cohorts.

Importance: Adolescent depression is characterized by diverse symptom trajectories over time and has a strong genetic influence. Research has determined genetic overlap between depression and other psychiatric conditions; investigating the shared genetic architecture of heterogeneous depression trajectories is crucial for understanding disease etiology, prediction, and early intervention. Objective: To investigate univariate and multivariate genetic risk for adolescent depression trajectories and assess generalizability across ancestries. Design, Setting, and Participants: This cohort study entailed longitudinal growth modeling followed by polygenic risk score (PRS) association testing for individual and multitrait genetic models. Two longitudinal cohorts from the US and UK were used: the Adolescent Brain and Cognitive Development (ABCD; N = 11 876) study and the Avon Longitudinal Study of Parents and Children (ALSPAC; N = 8787) study. Included were adolescents with genetic information and depression measures at up to 8 and 4 occasions, respectively. Study data were analyzed January to July 2023. Main Outcomes and Measures: Trajectories were derived from growth mixture modeling of longitudinal depression symptoms. PRSs were computed for depression, anxiety, neuroticism, bipolar disorder, schizophrenia, attention-deficit/hyperactivity disorder, and autism in European ancestry. Genomic structural equation modeling was used to build multitrait genetic models of psychopathology followed by multitrait PRS. Depression PRSs were computed in African, East Asian, and Hispanic ancestries in the ABCD cohort only. Association testing was performed between all PRSs and trajectories for both cohorts. Results: A total sample size of 14 112 adolescents (at baseline: mean [SD] age, 10.5 [0.5] years; 7269 male sex [52%]) from both cohorts were included in this analysis. Distinct depression trajectories (stable low, adolescent persistent, increasing, and decreasing) were replicated in the ALSPAC cohort (6096 participants; 3091 female [51%]) and ABCD cohort (8016 participants; 4274 male [53%]) between ages 10 and 17 years. Most univariate PRSs showed significant uniform associations with persistent trajectories, but fewer were significantly associated with intermediate (increasing and decreasing) trajectories. Multitrait PRSs-derived from a hierarchical factor model-showed the strongest associations for persistent trajectories (ABCD cohort: OR, 1.46; 95% CI, 1.26-1.68; ALSPAC cohort: OR, 1.34; 95% CI, 1.20-1.49), surpassing the effect size of univariate PRS in both cohorts. Multitrait PRSs were associated with intermediate trajectories but to a lesser extent (ABCD cohort: hierarchical increasing, OR, 1.27; 95% CI, 1.13-1.43; decreasing, OR, 1.23; 95% CI, 1.09-1.40; ALSPAC cohort: hierarchical increasing, OR, 1.16; 95% CI, 1.04-1.28; decreasing, OR, 1.32; 95% CI, 1.18-1.47). Transancestral genetic risk for depression showed no evidence for association with trajectories. Conclusions and Relevance: Results of this cohort study revealed a high multitrait genetic loading of persistent symptom trajectories, consistent across traits and cohorts. Variability in univariate genetic association with intermediate trajectories may stem from environmental factors. Multitrait genetics may strengthen depression prediction models, but more diverse data are needed for generalizability.

A comprehensive hierarchical comparison of structural connectomes in Major Depressive Disorder cases v. controls in two large population samples.

BACKGROUND: The brain can be represented as a network, with nodes as brain regions and edges as region-to-region connections. Nodes with the most connections (hubs) are central to efficient brain function. Current findings on structural differences in Major Depressive Disorder (MDD) identified using network approaches remain inconsistent, potentially due to small sample sizes. It is still uncertain at what level of the connectome hierarchy differences may exist, and whether they are concentrated in hubs, disrupting fundamental brain connectivity. METHODS: We utilized two large cohorts, UK Biobank (UKB, N = 5104) and Generation Scotland (GS, N = 725), to investigate MDD case-control differences in brain network properties. Network analysis was done across four hierarchical levels: (1) global, (2) tier (nodes grouped into four tiers based on degree) and rich club (between-hub connections), (3) nodal, and (4) connection. RESULTS: In UKB, reductions in network efficiency were observed in MDD cases globally (d = -0.076, pFDR = 0.033), across all tiers (d = -0.069 to -0.079, pFDR = 0.020), and in hubs (d = -0.080 to -0.113, pFDR = 0.013-0.035). No differences in rich club organization and region-to-region connections were identified. The effect sizes and direction for these associations were generally consistent in GS, albeit not significant in our lower-N replication sample. CONCLUSION: Our results suggest that the brain's fundamental rich club structure is similar in MDD cases and controls, but subtle topological differences exist across the brain. Consistent with recent large-scale neuroimaging findings, our findings offer a connectomic perspective on a similar scale and support the idea that minimal differences exist between MDD cases and controls.

Distinguishing different psychiatric disorders using DDx-PRS.

Despite great progress on methods for case-control polygenic prediction (e.g. schizophrenia vs. control), there remains an unmet need for a method that genetically distinguishes clinically related disorders (e.g. schizophrenia (SCZ) vs. bipolar disorder (BIP) vs. depression (MDD) vs. control); such a method could have important clinical value, especially at disorder onset when differential diagnosis can be challenging. Here, we introduce a method, Differential Diagnosis-Polygenic Risk Score (DDx-PRS), that jointly estimates posterior probabilities of each possible diagnostic category (e.g. SCZ=50%, BIP=25%, MDD=15%, control=10%) by modeling variance/covariance structure across disorders, leveraging case-control polygenic risk scores (PRS) for each disorder (computed using existing methods) and prior clinical probabilities for each diagnostic category. DDx-PRS uses only summary-level training data and does not use tuning data, facilitating implementation in clinical settings. In simulations, DDx-PRS was well-calibrated (whereas a simpler approach that analyzes each disorder marginally was poorly calibrated), and effective in distinguishing each diagnostic category vs. the rest. We then applied DDx-PRS to Psychiatric Genomics Consortium SCZ/BIP/MDD/control data, including summary-level training data from 3 case-control GWAS ( N =41,917-173,140 cases; total N =1,048,683) and held-out test data from different cohorts with equal numbers of each diagnostic category (total N =11,460). DDx-PRS was well-calibrated and well-powered relative to these training sample sizes, attaining AUCs of 0.66 for SCZ vs. rest, 0.64 for BIP vs. rest, 0.59 for MDD vs. rest, and 0.68 for control vs. rest. DDx-PRS produced comparable results to methods that leverage tuning data, confirming that DDx-PRS is an effective method. True diagnosis probabilities in top deciles of predicted diagnosis probabilities were considerably larger than prior baseline probabilities, particularly in projections to larger training sample sizes, implying considerable potential for clinical utility under certain circumstances. In conclusion, DDx-PRS is an effective method for distinguishing clinically related disorders.

Multi-ancestry genome-wide association study of major depression aids locus discovery, fine mapping, gene prioritization and causal inference.

Most genome-wide association studies (GWAS) of major depression (MD) have been conducted in samples of European ancestry. Here we report a multi-ancestry GWAS of MD, adding data from 21 cohorts with 88,316 MD cases and 902,757 controls to previously reported data. This analysis used a range of measures to define MD and included samples of African (36% of effective sample size), East Asian (26%) and South Asian (6%) ancestry and Hispanic/Latin American participants (32%). The multi-ancestry GWAS identified 53 significantly associated novel loci. For loci from GWAS in European ancestry samples, fewer than expected were transferable to other ancestry groups. Fine mapping benefited from additional sample diversity. A transcriptome-wide association study identified 205 significantly associated novel genes. These findings suggest that, for MD, increasing ancestral and global diversity in genetic studies may be particularly important to ensure discovery of core genes and inform about transferability of findings.

Comprehensive assessment of sleep duration, insomnia, and brain structure within the UK Biobank cohort.

STUDY OBJECTIVES: To assess for associations between sleeping more than or less than recommended by the National Sleep Foundation (NSF), and self-reported insomnia, with brain structure. METHODS: Data from the UK Biobank cohort were analyzed (N between 9K and 32K, dependent on availability, aged 44 to 82 years). Sleep measures included self-reported adherence to NSF guidelines on sleep duration (sleeping between 7 and 9 hours per night), and self-reported difficulty falling or staying asleep (insomnia). Brain structural measures included global and regional cortical or subcortical morphometry (thickness, surface area, volume), global and tract-related white matter microstructure, brain age gap (difference between chronological age and age estimated from brain scan), and total volume of white matter lesions. RESULTS: Longer-than-recommended sleep duration was associated with lower overall grey and white matter volumes, lower global and regional cortical thickness and volume measures, higher brain age gap, higher volume of white matter lesions, higher mean diffusivity globally and in thalamic and association fibers, and lower volume of the hippocampus. Shorter-than-recommended sleep duration was related to higher global and cerebellar white matter volumes, lower global and regional cortical surface areas, and lower fractional anisotropy in projection fibers. Self-reported insomnia was associated with higher global gray and white matter volumes, and with higher volumes of the amygdala, hippocampus, and putamen. CONCLUSIONS: Sleeping longer than recommended by the NSF is associated with a wide range of differences in brain structure, potentially indicative of poorer brain health. Sleeping less than recommended is distinctly associated with lower cortical surface areas. Future studies should assess the potential mechanisms of these differences and investigate long sleep duration as a putative marker of brain health.

Meta-Analyses of Genome-Wide Association Studies for Postpartum Depression.

OBJECTIVE: Postpartum depression (PPD) is a common subtype of major depressive disorder (MDD) that is more heritable, yet is understudied in psychiatric genetics. The authors conducted meta-analyses of genome-wide association studies (GWASs) to investigate the genetic architecture of PPD. METHOD: Meta-analyses were conducted on 18 cohorts of European ancestry (17,339 PPD cases and 53,426 controls), one cohort of East Asian ancestry (975 cases and 3,780 controls), and one cohort of African ancestry (456 cases and 1,255 controls), totaling 18,770 PPD cases and 58,461 controls. Post-GWAS analyses included 1) single-nucleotide polymorphism (SNP)-based heritability ([Formula: see text]), 2) genetic correlations between PPD and other phenotypes, and 3) enrichment of the PPD GWAS findings in 27 human tissues and 265 cell types from the mouse central and peripheral nervous system. RESULTS: No SNP achieved genome-wide significance in the European or the trans-ancestry meta-analyses. The [Formula: see text] of PPD was 0.14 (SE=0.02). Significant genetic correlations were estimated for PPD with MDD, bipolar disorder, anxiety disorders, posttraumatic stress disorder, insomnia, age at menarche, and polycystic ovary syndrome. Cell-type enrichment analyses implicate inhibitory neurons in the thalamus and cholinergic neurons within septal nuclei of the hypothalamus, a pattern that differs from MDD. CONCLUSIONS: While more samples are needed to reach genome-wide levels of significance, the results presented confirm PPD as a polygenic and heritable phenotype. There is also evidence that despite a high correlation with MDD, PPD may have unique genetic components. Cell enrichment results suggest GABAergic neurons, which converge on a common mechanism with the only medication approved by the U.S. Food and Drug Administration for PPD (brexanolone).

Pathway-Based Polygenic Risk Scores for Schizophrenia and Associations With Reported Psychotic-like Experiences and Neuroimaging Phenotypes in the UK Biobank.

Background: Schizophrenia is a heritable psychiatric disorder with a polygenic architecture. Genome-wide association studies have reported that an increasing number of risk-associated variants and polygenic risk scores (PRSs) explain 17% of the variance in the disorder. Substantial heterogeneity exists in the effect of these variants, and aggregating them based on biologically relevant functions may provide mechanistic insight into the disorder. Methods: Using the largest schizophrenia genome-wide association study conducted to date, we associated PRSs based on 5 gene sets previously found to contribute to schizophrenia pathophysiology-postsynaptic density of excitatory synapses, postsynaptic membrane, dendritic spine, axon, and histone H3-K4 methylation-along with respective whole-genome PRSs, with neuroimaging (n > 29,000) and reported psychotic-like experiences (n > 119,000) variables in healthy UK Biobank subjects. Results: Several variables were significantly associated with the axon gene-set (psychotic-like communications, parahippocampal gyrus volume, fractional anisotropy thalamic radiations, and fractional anisotropy posterior thalamic radiations (ß range -0.016 to 0.0916, false discovery rate-corrected p [pFDR] ≤ .05), postsynaptic density gene-set (psychotic-like experiences distress, global surface area, and cingulate lobe surface area [ß range -0.014 to 0.0588, pFDR ≤ .05]), and histone gene set (entorhinal surface area: ß = -0.016, pFDR = .035). From these, whole-genome PRSs were significantly associated with psychotic-like communications (ß = 0.2218, pFDR = 1.34 × 10-7), distress (ß = 0.1943, pFDR = 7.28 × 10-16), and fractional anisotropy thalamic radiations (ß = -0.0143, pFDR = .036). Permutation analysis revealed that these associations were not due to chance. Conclusions: Our results indicate that genetic variation in 3 gene sets relevant to schizophrenia may confer risk for the disorder through effects on previously implicated neuroimaging variables. Because associations were stronger overall for whole-genome PRSs, findings here highlight that selection of biologically relevant variants is not yet sufficient to address the heterogeneity of the disorder.

Genetic structure of major depression symptoms across clinical and community cohorts.

Diagnostic criteria for major depressive disorder allow for heterogeneous symptom profiles but genetic analysis of major depressive symptoms has the potential to identify clinical and aetiological subtypes. There are several challenges to integrating symptom data from genetically-informative cohorts, such as sample size differences between clinical and community cohorts and various patterns of missing data. We conducted genome-wide association studies of major depressive symptoms in three clinical cohorts that were enriched for affected participants (Psychiatric Genomics Consortium, Australian Genetics of Depression Study, Generation Scotland) and three community cohorts (Avon Longitudinal Study of Parents and Children, Estonian Biobank, and UK Biobank). We fit a series of confirmatory factor models with factors that accounted for how symptom data was sampled and then compared alternative models with different symptom factors. The best fitting model had a distinct factor for Appetite/Weight symptoms and an additional measurement factor that accounted for missing data patterns in the community cohorts (use of Depression and Anhedonia as gating symptoms). The results show the importance of assessing the directionality of symptoms (such as hypersomnia versus insomnia) and of accounting for study and measurement design when meta-analysing genetic association data.

Polygenic risk prediction: why and when out-of-sample prediction R2 can exceed SNP-based heritability.

In polygenic score (PGS) analysis, the coefficient of determination (R2) is a key statistic to evaluate efficacy. R2 is the proportion of phenotypic variance explained by the PGS, calculated in a cohort that is independent of the genome-wide association study (GWAS) that provided estimates of allelic effect sizes. The SNP-based heritability (hSNP2, the proportion of total phenotypic variances attributable to all common SNPs) is the theoretical upper limit of the out-of-sample prediction R2. However, in real data analyses R2 has been reported to exceed hSNP2, which occurs in parallel with the observation that hSNP2 estimates tend to decline as the number of cohorts being meta-analyzed increases. Here, we quantify why and when these observations are expected. Using theory and simulation, we show that if heterogeneities in cohort-specific hSNP2 exist, or if genetic correlations between cohorts are less than one, hSNP2 estimates can decrease as the number of cohorts being meta-analyzed increases. We derive conditions when the out-of-sample prediction R2 will be greater than hSNP2 and show the validity of our derivations with real data from a binary trait (major depression) and a continuous trait (educational attainment). Our research calls for a better approach to integrating information from multiple cohorts to address issues of between-cohort heterogeneity.

Phenome-wide analyses identify an association between the parent-of-origin effects dependent methylome and the rate of aging in humans.

BACKGROUND: The variation in the rate at which humans age may be rooted in early events acting through the genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions enriched for genetically controlled imprinting effects (the typical type of POE) and regions influenced by environmental effects associated with parents (the atypical POE). This part of the methylome is heavily influenced by early events, making it a potential route connecting early exposures, the epigenome, and aging. We aim to test the association of POE-CpGs with early and later exposures and subsequently with health-related phenotypes and adult aging. RESULTS: We perform a phenome-wide association analysis for the POE-influenced methylome using GS:SFHS (Ndiscovery = 5087, Nreplication = 4450). We identify and replicate 92 POE-CpG-phenotype associations. Most of the associations are contributed by the POE-CpGs belonging to the atypical class where the most strongly enriched associations are with aging (DNAmTL acceleration), intelligence, and parental (maternal) smoking exposure phenotypes. A proportion of the atypical POE-CpGs form co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased within-module methylation connectivity with age. The atypical POE-CpGs also display high levels of methylation heterogeneity, fast information loss with age, and a strong correlation with CpGs contained within epigenetic clocks. CONCLUSIONS: These results identify the association between the atypical POE-influenced methylome and aging and provide new evidence for the "early development of origin" hypothesis for aging in humans.

A Meta-Analysis of the Genome-Wide Association Studies on Two Genetically Correlated Phenotypes Suggests Four New Risk Loci for Headaches.

Headache is one of the commonest complaints that doctors need to address in clinical settings. The genetic mechanisms of different types of headache are not well understood while it has been suggested that self-reported headache and self-reported migraine were genetically correlated. In this study, we performed a meta-analysis of genome-wide association studies (GWAS) on the self-reported headache phenotype from the UK Biobank and the self-reported migraine phenotype from the 23andMe using the Unified Score-based Association Test (metaUSAT) software for genetically correlated phenotypes (N = 397,385). We identified 38 loci for headaches, of which 34 loci have been reported before and four loci were newly suggested. The LDL receptor related protein 1 (LRP1)-Signal Transducer and Activator of Transcription 6 (STAT6)-S hort chain D ehydrogenase/R eductase family 9C member 7 (SDR9C7) region in chromosome 12 was the most significantly associated locus with a leading p value of 1.24 × 10-62 of rs11172113. The One Cut homeobox 2 (ONECUT2) gene locus in chromosome 18 was the strongest signal among the four new loci with a p value of 1.29 × 10-9 of rs673939. Our study demonstrated that the genetically correlated phenotypes of self-reported headache and self-reported migraine can be meta-analysed together in theory and in practice to boost study power to identify more variants for headaches. This study has paved way for a large GWAS meta-analysis involving cohorts of different while genetically correlated headache phenotypes. Supplementary Information: The online version contains supplementary material available at 10.1007/s43657-022-00078-7.

Phenome-wide analysis identifies parent-of-origin effects on the human methylome associated with changes in the rate of aging.

Variation in the rate at which humans age may be rooted in early life events acting through genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions either enriched for genetically controlled imprinting effects (the typical type of POE) or atypical POE introduced by environmental effects associated with parents. This part of the methylome is heavily influenced by early life events, making it a potential route connecting early environmental exposures, the epigenome and the rate of aging. Here, we aim to test the association of POE-influenced methylation of CpG dinucleotides (POE-CpG sites) with early and later environmental exposures and subsequently with health-related phenotypes and adult aging phenotypes. We do this by performing phenome-wide association analyses of the POE-influenced methylome using a large family-based population cohort (GS:SFHS, Ndiscovery=5,087, Nreplication=4,450). At the single CpG level, 92 associations of POE-CpGs with phenotypic variation were identified and replicated. Most of the associations were contributed by POE-CpGs belonging to the atypical class and the most strongly enriched associations were with aging (DNAmTL acceleration), intelligence and parental (maternal) smoking exposure phenotypes. We further found that a proportion of the atypical-POE-CpGs formed co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased internal module connectivity (strength of methylation correlation across constituent CpGs) with age. Atypical POE-CpGs also displayed high levels of methylation heterogeneity and epigenetic drift (i.e. information loss with age) and a strong correlation with CpGs contained within epigenetic clocks. These results identified associations between the atypical-POE-influenced methylome and aging and provided new evidence for the "early development of origin" hypothesis for aging in humans.

Associations of negative affective biases and depressive symptoms in a community-based sample.

BACKGROUND: Major depressive disorder (MDD) was previously associated with negative affective biases. Evidence from larger population-based studies, however, is lacking, including whether biases normalise with remission. We investigated associations between affective bias measures and depressive symptom severity across a large community-based sample, followed by examining differences between remitted individuals and controls. METHODS: Participants from Generation Scotland (N = 1109) completed the: (i) Bristol Emotion Recognition Task (BERT), (ii) Face Affective Go/No-go (FAGN), and (iii) Cambridge Gambling Task (CGT). Individuals were classified as MDD-current (n = 43), MDD-remitted (n = 282), or controls (n = 784). Analyses included using affective bias summary measures (primary analyses), followed by detailed emotion/condition analyses of BERT and FAGN (secondary analyses). RESULTS: For summary measures, the only significant finding was an association between greater symptoms and lower risk adjustment for CGT across the sample (individuals with greater symptoms were less likely to bet more, despite increasingly favourable conditions). This was no longer significant when controlling for non-affective cognition. No differences were found for remitted-MDD v. controls. Detailed analysis of BERT and FAGN indicated subtle negative biases across multiple measures of affective cognition with increasing symptom severity, that were independent of non-effective cognition [e.g. greater tendency to rate faces as angry (BERT), and lower accuracy for happy/neutral conditions (FAGN)]. Results for remitted-MDD were inconsistent. CONCLUSIONS: This suggests the presence of subtle negative affective biases at the level of emotion/condition in association with depressive symptoms across the sample, over and above those accounted for by non-affective cognition, with no evidence for affective biases in remitted individuals.

Early-life inflammatory markers and subsequent psychotic and depressive episodes between 10 to 28 years of age.

Inflammation is implicated in depression and psychosis, including association of childhood inflammatory markers on the subsequent risk of developing symptoms. However, it is unknown whether early-life inflammatory markers are associated with the number of depressive and psychotic symptoms from childhood to adulthood. Using the prospective Avon Longitudinal Study of Children and Parents birth cohort (N = up-to 6401), we have examined longitudinal associations of early-life inflammation [exposures: interleukin-6 (IL-6), C-reactive protein (CRP) levels at age 9y; IL-6 and CRP DNA-methylation (DNAm) scores at birth and age 7y; and IL-6 and CRP polygenic risk scores (PRSs)] with the number of depressive episodes and psychotic experiences (PEs) between ages 10-28 years. Psychiatric outcomes were assessed using the Short Mood and Feelings Questionnaire and Psychotic Like Symptoms Questionnaires, respectively. Exposure-outcome associations were tested using negative binomial models, which were adjusted for metabolic and sociodemographic factors. Serum IL-6 levels at age 9y were associated with the total number of depressive episodes between 10 and 28y in the base model (n = 4835; ß = 0.066; 95%CI:0.020-0.113; pFDR = 0.041) which was weaker when adjusting for metabolic and sociodemographic factors. Weak associations were observed between inflammatory markers (serum IL-6 and CRP DNAm scores) and total number of PEs. Other inflammatory markers were not associated with depression or PEs. Early-life inflammatory markers are associated with the burden of depressive episodes and of PEs subsequently from childhood to adulthood. These findings support a potential role of early-life inflammation in the aetiology of depression and psychosis and highlight inflammation as a potential target for treatment and prevention.

Genome-by-Trauma Exposure Interactions in Adults With Depression in the UK Biobank.

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.

Corrigendum: Potential genetic overlap between insomnia and sleep symptoms in major depressive disorder: A polygenic risk score analysis.

[This corrects the article DOI: 10.3389/fpsyt.2021.734077.].

Comparing personalized brain-based and genetic risk scores for major depressive disorder in large population samples of adults and adolescents.

BACKGROUND: Major depressive disorder (MDD) is a polygenic disorder associated with brain alterations but until recently, there have been no brain-based metrics to quantify individual-level variation in brain morphology. Here, we evaluated and compared the performance of a new brain-based 'Regional Vulnerability Index' (RVI) with polygenic risk scores (PRS), in the context of MDD. We assessed associations with syndromal MDD in an adult sample (N = 702, age = 59 ± 10) and with subclinical depressive symptoms in a longitudinal adolescent sample (baseline N = 3,825, age = 10 ± 1; 2-year follow-up N = 2,081, age = 12 ± 1). METHODS: MDD-RVIs quantify the correlation of the individual's corresponding brain metric with the expected pattern for MDD derived in an independent sample. Using the same methodology across samples, subject-specific MDD-PRS and six MDD-RVIs based on different brain modalities (subcortical volume, cortical thickness, cortical surface area, mean diffusivity, fractional anisotropy, and multimodal) were computed. RESULTS: In adults, MDD-RVIs (based on white matter and multimodal measures) were more strongly associated with MDD (ß = 0.099-0.281, PFDR = 0.001-0.043) than MDD-PRS (ß = 0.056-0.152, PFDR = 0.140-0.140). In adolescents, depressive symptoms were associated with MDD-PRS at baseline and follow-up (ß = 0.084-0.086, p = 1.38 × 10-4-4.77 × 10-4) but not with any MDD-RVIs (ß < 0.05, p > 0.05). CONCLUSIONS: Our results potentially indicate the ability of brain-based risk scores to capture a broader range of risk exposures than genetic risk scores in adults and are also useful in helping us to understand the temporal origins of depression-related brain features. Longitudinal data, specific to the developmental period and on white matter measures, will be useful in informing risk for subsequent psychiatric illness.