Polygenic prediction of major depressive disorder and related traits in African ancestries UK Biobank participants

Genome-Wide Association Studies (GWAS) over-represent European ancestries, neglecting all other ancestry groups and low-income nations. Consequently, polygenic risk scores (PRS) more accurately predict complex traits in Europeans than African Ancestries groups. Very few studies have looked at the transferability of European-derived PRS for behavioural and mental health phenotypes to Africans. We assessed the comparative accuracy of depression PRS trained on European and African Ancestries GWAS studies to predict major depressive disorder (MDD) and related traits in African ancestry participants from the UK Biobank. UK Biobank participants were selected based on Principal component analysis clustering with an African genetic similarity reference population, MDD was assessed with the Composite International Diagnostic Interview (CIDI). PRS were computed using PRSice2 software using either European or African Ancestries GWAS summary statistics. PRS trained on European ancestry samples (246,363 cases) predicted case control status in Africans of the UK Biobank with similar accuracies (R2 = 2%, β = 0.32, empirical p-value = 0.002) to PRS trained on far much smaller samples of African Ancestries participants from 23andMe, Inc. (5045 cases, R² = 1.8%, β = 0.28, empirical p-value = 0.008). This suggests that prediction of MDD status from Africans to Africans had greater efficiency relative to discovery sample size than prediction of MDD from Europeans to Africans. Prediction of MDD status in African UK Biobank participants using GWAS findings of likely causal risk factors from European ancestries was non-significant. GWAS of MDD in European ancestries are inefficient for improving polygenic prediction in African samples; urgent MDD studies in Africa are needed.

Genome-wide association study of major anxiety disorders in 122,341 European-ancestry cases identifies 58 loci and highlights GABAergic signaling

The major anxiety disorders (ANX; including generalized anxiety disorder, panic disorder, and phobias) are highly prevalent, often onset early, persist throughout life, and cause substantial global disability. Although distinct in their clinical presentations, they likely represent differential expressions of a dysregulated threat-response system. Here we present a genome-wide association meta-analysis comprising 122,341 European ancestry ANX cases and 729,881 controls. We identified 58 independent genome-wide significant ANX risk variants and 66 genes with robust biological support. In an independent sample of 1,175,012 self-report ANX cases and 1,956,379 controls, 51 of the 58 associated variants were replicated. As predicted by twin studies, we found substantial genetic correlation between ANX and depression, neuroticism, and other internalizing phenotypes. Follow-up analyses demonstrated enrichment in all major brain regions and highlighted GABAergic signaling as one potential mechanism underlying ANX genetic risk. These results advance our understanding of the genetic architecture of ANX and prioritize genes for functional follow-up studies.

Genome-wide study of half a million individuals with major depression identifies 697 independent associations, infers causal neuronal subtypes and biological targets for novel pharmacotherapies

In a genome-wide association study (GWAS) 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 genetic loci, 293 of which are novel. Using fine-mapping and functional genomic datasets, 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 implicated excitatory and inhibitory midbrain and forebrain neurons, peptidergic neurons, and medium spiny neurons in MD. Critically, our findings are enriched for the targets of antidepressants and provide potential antidepressant repurposing opportunities (e.g., pregabalin and modafinil). Polygenic scores (PGS) from European ancestries explained up to 5.7% of the variance in liability to MD in European samples and PGS trained using either European or multi-ancestry data significantly predicted case control status across all four diverse ancestries. These findings represent a major advance in our understanding of MD across global populations. We provide evidence that 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.

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.

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.

Polygenic prediction of major depressive disorder and related traits in African ancestries UK Biobank participants

Introduction

Genome-Wide Association Studies (GWAS) over-represent European ancestries compared to the global population, neglecting all other ancestry groups and low-income nations. Consequently, polygenic risk scores (PRS) more accurately predict complex traits in Europeans than African Ancestries groups. Very few studies have looked at the transferability of European-derived PRS for behavioural and mental health phenotypes to non-Europeans. We assessed the comparative accuracy of PRS for Major Depressive Disorder (MDD) trained on European and African Ancestries GWAS studies to predict MDD and related traits in African Ancestries participants from the UK Biobank.

Methods

UK Biobank participants were selected based on Principal component analysis (PCA) clustering with an African genetic similarity reference population and MDD was assessed with the Composite International Diagnostic Interview (CIDI). Polygenic Risk Scores (PRS) were computed using PRSice2 using either European or African Ancestries GWAS summary statistics.

Results

PRS trained on European ancestry samples (246,363 cases) predicted case control status in Africans of the UK Biobank with similar accuracies (190 cases, R2=2%) to PRS trained on far much smaller samples of African Ancestries participants from 23andMe, Inc. (5045 cases, R2=1.8%). This suggests that prediction of MDD status from Africans to Africans had greater efficiency per unit increase in the discovery sample size than prediction of MDD from Europeans to Africans. Prediction of MDD status in African UK Biobank participants using GWAS findings of causal risk factors from European ancestries was non-significant.

Conclusion

GWAS studies of MDD in European ancestries are an inefficient means of improving polygenic prediction accuracy in African samples.

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.

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.

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.

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.

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.

Complex trait methylation scores in the prediction of major depressive disorder

BACKGROUND

DNA methylation (DNAm) is associated with time-varying environmental factors that contribute to major depressive disorder (MDD) risk. We sought to test whether DNAm signatures of lifestyle and biochemical factors were associated with MDD to reveal dynamic biomarkers of MDD risk that may be amenable to lifestyle interventions.

METHODS

Here, we calculated methylation scores (MS) at multiple p-value thresholds for lifestyle (BMI, smoking, alcohol consumption, and educational attainment) and biochemical (high-density lipoprotein (HDL) and total cholesterol) factors in Generation Scotland (GS) (N=9,502) and in a replication cohort (ALSPAC_adults, N=565), using CpG sites reported in previous well-powered methylome-wide association studies. We also compared their predictive accuracy for MDD to a MDD MS in an independent GS sub-sample (N=4,432).

FINDINGS

Each trait MS was significantly associated with its corresponding phenotype in GS (β_range=0.089-1.457) and in ALSPAC (β_range=0.078-2.533). Each MS was also significantly associated with MDD before and after adjustment for its corresponding phenotype in GS (β_range=0.053-0.145). After accounting for relevant lifestyle factors, MS for educational attainment (β=0.094) and alcohol consumption (MS_p-value<0.01-0.5; β_range=-0.069-0.083) remained significantly associated with MDD in GS. Smoking (AUC=0.569) and educational attainment (AUC=0.585) MSs could discriminate MDD from controls better than the MDD MS (AUC=0.553) in the independent GS sub-sample. Analyses implicating MDD did not replicate across ALSPAC, although the direction of effect was consistent for all traits when adjusting for the MS corresponding phenotypes.

INTERPRETATION

We showed that lifestyle and biochemical MS were associated with MDD before and after adjustment for their corresponding phenotypes (p_nominal<0.05), but not when smoking, alcohol consumption, and BMI were also included as covariates. MDD results did not replicate in the smaller, female-only independent ALSPAC cohort (N_ALSPAC=565; N_GS=9,502), potentially due to demographic differences or low statistical power, but effect sizes were consistent with the direction reported in GS. DNAm scores for modifiable MDD risk factors may contribute to disease vulnerability and, in some cases, explain additional variance to their observed phenotypes.

FUNDING

Wellcome Trust.

Genetic architecture of 11 major psychiatric disorders at biobehavioral, functional genomic and molecular genetic levels of analysis

We interrogate the joint genetic architecture of 11 major psychiatric disorders at biobehavioral, functional genomic and molecular genetic levels of analysis. We identify four broad factors (neurodevelopmental, compulsive, psychotic and internalizing) that underlie genetic correlations among the disorders and test whether these factors adequately explain their genetic correlations with biobehavioral traits. We introduce stratified genomic structural equation modeling, which we use to identify gene sets that disproportionately contribute to genetic risk sharing. This includes protein-truncating variant-intolerant genes expressed in excitatory and GABAergic brain cells that are enriched for genetic overlap across disorders with psychotic features. Multivariate association analyses detect 152 (20 new) independent loci that act on the individual factors and identify nine loci that act heterogeneously across disorders within a factor. Despite moderate-to-high genetic correlations across all 11 disorders, we find little utility of a single dimension of genetic risk across psychiatric disorders either at the level of biobehavioral correlates or at the level of individual variants.

Structural neuroimaging measures and lifetime depression across levels of phenotyping in UK biobank

Depression is assessed in various ways in research, with large population studies often relying on minimal phenotyping. Genetic results suggest clinical diagnoses and self-report measures of depression show some core similarities, but also important differences. It is not yet clear how neuroimaging associations depend on levels of phenotyping. We studied 39,300 UK Biobank imaging participants (20,701 female; aged 44.6 to 82.3 years, M = 64.1, SD = 7.5) with structural neuroimaging and lifetime depression data. Past depression phenotypes included a single-item self-report measure, an intermediate measure of 'probable' lifetime depression, derived from multiple questionnaire items relevant to a history of depression, and a retrospective clinical diagnosis according to DSM-IV criteria. We tested (i) associations between brain structural measures and each depression phenotype, and (ii) effects of phenotype on these associations. Depression-brain structure associations were small (β < 0.1) for all phenotypes, but still significant after FDR correction for many regional metrics. Lifetime depression was consistently associated with reduced white matter integrity across phenotypes. Cortical thickness showed negative associations with Self-reported Depression in particular. Phenotype effects were small across most metrics, but significant for cortical thickness in most regions. We report consistent effects of lifetime depression in brain structural measures, including reduced integrity of thalamic radiations and association fibres. We also observed significant differences in associations with cortical thickness across depression phenotypes. Although these results did not relate to level of phenotyping as expected, effects of phenotype definition are still an important consideration for future depression research.