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.

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.

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.

IgG fractions from patients with antiphospholipid syndrome and systemic lupus erythematosus bind to platelets, but do not affect collagen-induced platelet activation

Anti-beta-2 glycoprotein 1 (anti-β2GP1) is an antiphospholipid antibody found in patients with antiphospholipid syndrome (APS) and systemic lupus erythematosus (SLE). Its presence commonly is associated with thrombosis; however, the mechanisms of interaction of anti-β2GP1 antibodies and platelets remain unclear. We investigated the effects of APS and SLE patient-derived IgG fractions on collagen-mediated platelet aggregation and examined the binding of patient-derived IgG to platelets before and after activation by collagen. IgG fractions, 150, 200, 300 or 350 µg/ml, isolated from 11 patients with APS and SLE were incubated with two sets of platelet-rich plasma (PRP) in the incubation wells of an aggregometer. The first set was activated by collagen and the other set was incubated for an additional 10 min. All platelets were collected by centrifugation and fixed in cell blocks. We assessed binding of IgG to platelets using immunocytochemistry (ICC). Patient-derived IgG fractions did not affect collagen-induced platelet aggregation. ICC staining using anti-human IgG antibodies demonstrated that patient-derived IgG fractions had greater affinity for non-activated platelets than those activated by 0.75 µg/ml collagen. Patient-derived IgG fractions bound to the surface of platelets and potentially could be internalized by platelets. IgG fractions from APS and SLE patients may sensitize non-activated platelets, which could increase platelet reactivity and thrombotic risk in patients. We did not detect secondary effects of patient-derived IgG fractions.

Identifying the Common Genetic Basis of Antidepressant Response

Background

Antidepressants are a first-line treatment for depression. However, only a third of individuals experience remission after the first treatment. Common genetic variation, in part, likely regulates antidepressant response, yet the success of previous genome-wide association studies has been limited by sample size. This study performs the largest genetic analysis of prospectively assessed antidepressant response in major depressive disorder to gain insight into the underlying biology and enable out-of-sample prediction.

Methods

Genome-wide analysis of remission (n remit = 1852, n nonremit = 3299) and percentage improvement (n = 5218) was performed. Single nucleotide polymorphism-based heritability was estimated using genome-wide complex trait analysis. Genetic covariance with eight mental health phenotypes was estimated using polygenic scores/AVENGEME. Out-of-sample prediction of antidepressant response polygenic scores was assessed. Gene-level association analysis was performed using MAGMA and transcriptome-wide association study. Tissue, pathway, and drug binding enrichment were estimated using MAGMA.

Results

Neither genome-wide association study identified genome-wide significant associations. Single nucleotide polymorphism-based heritability was significantly different from zero for remission (h 2 = 0.132, SE = 0.056) but not for percentage improvement (h 2 = -0.018, SE = 0.032). Better antidepressant response was negatively associated with genetic risk for schizophrenia and positively associated with genetic propensity for educational attainment. Leave-one-out validation of antidepressant response polygenic scores demonstrated significant evidence of out-of-sample prediction, though results varied in external cohorts. Gene-based analyses identified ETV4 and DHX8 as significantly associated with antidepressant response.

Conclusions

This study demonstrates that antidepressant response is influenced by common genetic variation, has a genetic overlap schizophrenia and educational attainment, and provides a useful resource for future research. Larger sample sizes are required to attain the potential of genetics for understanding and predicting antidepressant response.

DNA methylome-wide association study of genetic risk for depression implicates antigen processing and immune responses

BACKGROUND

Depression is a disabling and highly prevalent condition where genetic and epigenetic, such as DNA methylation (DNAm), differences contribute to disease risk. DNA methylation is influenced by genetic variation but the association between polygenic risk of depression and DNA methylation is unknown.

METHODS

We investigated the association between polygenic risk scores (PRS) for depression and DNAm by conducting a methylome-wide association study (MWAS) in Generation Scotland (N = 8898, mean age = 49.8 years) with replication in the Lothian Birth Cohorts of 1921 and 1936 and adults in the Avon Longitudinal Study of Parents and Children (ALSPAC) (Ncombined = 2049, mean age = 79.1, 69.6 and 47.2 years, respectively). We also conducted a replication MWAS in the ALSPAC children (N = 423, mean age = 17.1 years). Gene ontology analysis was conducted for the cytosine-guanine dinucleotide (CpG) probes significantly associated with depression PRS, followed by Mendelian randomisation (MR) analysis to infer the causal relationship between depression and DNAm.

RESULTS

Widespread associations (NCpG = 71, pBonferroni < 0.05, p < 6.3 × 10-8) were found between PRS constructed using genetic risk variants for depression and DNAm in CpG probes that localised to genes involved in immune responses and neural development. The effect sizes for the significant associations were highly correlated between the discovery and replication samples in adults (r = 0.79) and in adolescents (r = 0.82). Gene Ontology analysis showed that significant CpG probes are enriched in immunological processes in the human leukocyte antigen system. Additional MWAS was conducted for each lead genetic risk variant. Over 47.9% of the independent genetic risk variants included in the PRS showed associations with DNAm in CpG probes located in both the same (cis) and distal (trans) locations to the genetic loci (pBonferroni < 0.045). Subsequent MR analysis showed that there are a greater number of causal effects found from DNAm to depression than vice versa (DNAm to depression: pFDR ranged from 0.024 to 7.45 × 10-30; depression to DNAm: pFDR ranged from 0.028 to 0.003).

CONCLUSIONS

PRS for depression, especially those constructed from genome-wide significant genetic risk variants, showed methylome-wide differences associated with immune responses. Findings from MR analysis provided evidence for causal effect of DNAm to depression.

Cohort profile for the STratifying Resilience and Depression Longitudinally (STRADL) study: A depression-focused investigation of Generation Scotland, using detailed clinical, cognitive, and neuroimaging assessments

STratifying Resilience and Depression Longitudinally (STRADL) is a population-based study built on the Generation Scotland: Scottish Family Health Study (GS:SFHS) resource. The aim of STRADL is to subtype major depressive disorder (MDD) on the basis of its aetiology, using detailed clinical, cognitive, and brain imaging assessments. The GS:SFHS provides an important opportunity to study complex gene-environment interactions, incorporating linkage to existing datasets and inclusion of early-life variables for two longitudinal birth cohorts. Specifically, data collection in STRADL included: socio-economic and lifestyle variables; physical measures; questionnaire data that assesses resilience, early-life adversity, personality, psychological health, and lifetime history of mood disorder; laboratory samples; cognitive tests; and brain magnetic resonance imaging. Some of the questionnaire and cognitive data were first assessed at the GS:SFHS baseline assessment between 2006-2011, thus providing longitudinal measures relevant to the study of depression, psychological resilience, and cognition. In addition, routinely collected historic NHS data and early-life variables are linked to STRADL data, further providing opportunities for longitudinal analysis. Recruitment has been completed and we consented and tested 1,188 participants.

Methylome-wide association study of antidepressant use in Generation Scotland and the Netherlands Twin Register implicates the innate immune system

Antidepressants are an effective treatment for major depressive disorder (MDD), although individual response is unpredictable and highly variable. Whilst the mode of action of antidepressants is incompletely understood, many medications are associated with changes in DNA methylation in genes that are plausibly linked to their mechanisms. Studies of DNA methylation may therefore reveal the biological processes underpinning the efficacy and side effects of antidepressants. We performed a methylome-wide association study (MWAS) of self-reported antidepressant use accounting for lifestyle factors and MDD in Generation Scotland (GS:SFHS, N = 6428, EPIC array) and the Netherlands Twin Register (NTR, N = 2449, 450 K array) and ran a meta-analysis of antidepressant use across these two cohorts. We found ten CpG sites significantly associated with self-reported antidepressant use in GS:SFHS, with the top CpG located within a gene previously associated with mental health disorders, ATP6V1B2 (β = -0.055, pcorrected = 0.005). Other top loci were annotated to genes including CASP10, TMBIM1, MAPKAPK3, and HEBP2, which have previously been implicated in the innate immune response. Next, using penalised regression, we trained a methylation-based score of self-reported antidepressant use in a subset of 3799 GS:SFHS individuals that predicted antidepressant use in a second subset of GS:SFHS (N = 3360, β = 0.377, p = 3.12 × 10-11, R2 = 2.12%). In an MWAS analysis of prescribed selective serotonin reuptake inhibitors, we showed convergent findings with those based on self-report. In NTR, we did not find any CpGs significantly associated with antidepressant use. The meta-analysis identified the two CpGs of the ten above that were common to the two arrays used as being significantly associated with antidepressant use, although the effect was in the opposite direction for one of them. Antidepressants were associated with epigenetic alterations in loci previously associated with mental health disorders and the innate immune system. These changes predicted self-reported antidepressant use in a subset of GS:SFHS and identified processes that may be relevant to our mechanistic understanding of clinically relevant antidepressant drug actions and side effects.

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.].

Sex-Dependent Shared and Nonshared Genetic Architecture Across Mood and Psychotic Disorders

BACKGROUND

Sex differences in incidence and/or presentation of schizophrenia (SCZ), major depressive disorder (MDD), and bipolar disorder (BIP) are pervasive. Previous evidence for shared genetic risk and sex differences in brain abnormalities across disorders suggest possible shared sex-dependent genetic risk.

METHODS

We conducted the largest to date genome-wide genotype-by-sex (G×S) interaction of risk for these disorders using 85,735 cases (33,403 SCZ, 19,924 BIP, and 32,408 MDD) and 109,946 controls from the PGC (Psychiatric Genomics Consortium) and iPSYCH.

RESULTS

Across disorders, genome-wide significant single nucleotide polymorphism-by-sex interaction was detected for a locus encompassing NKAIN2 (rs117780815, p = 3.2 × 10-8), which interacts with sodium/potassium-transporting ATPase (adenosine triphosphatase) enzymes, implicating neuronal excitability. Three additional loci showed evidence (p < 1 × 10-6) for cross-disorder G×S interaction (rs7302529, p = 1.6 × 10-7; rs73033497, p = 8.8 × 10-7; rs7914279, p = 6.4 × 10-7), implicating various functions. Gene-based analyses identified G×S interaction across disorders (p = 8.97 × 10-7) with transcriptional inhibitor SLTM. Most significant in SCZ was a MOCOS gene locus (rs11665282, p = 1.5 × 10-7), implicating vascular endothelial cells. Secondary analysis of the PGC-SCZ dataset detected an interaction (rs13265509, p = 1.1 × 10-7) in a locus containing IDO2, a kynurenine pathway enzyme with immunoregulatory functions implicated in SCZ, BIP, and MDD. Pathway enrichment analysis detected significant G×S interaction of genes regulating vascular endothelial growth factor receptor signaling in MDD (false discovery rate-corrected p < .05).

CONCLUSIONS

In the largest genome-wide G×S analysis of mood and psychotic disorders to date, there was substantial genetic overlap between the sexes. However, significant sex-dependent effects were enriched for genes related to neuronal development and immune and vascular functions across and within SCZ, BIP, and MDD at the variant, gene, and pathway levels.

Lifestyle and Genetic Factors Modify Parent-of-Origin Effects on the Human Methylome

BACKGROUND

parent-of-origin effects (POE) play important roles in complex disease and thus understanding their regulation and associated molecular and phenotypic variation are warranted. Previous studies mainly focused on the detection of genomic regions or phenotypes regulated by POE. Understanding whether POE may be modified by environmental or genetic exposures is important for understanding of the source of POE-associated variation, but only a few case studies addressing modifiable POE exist.

METHODS

in order to understand this high order of POE regulation, we screened 101 genetic and environmental factors such as 'predicted mRNA expression levels' of DNA methylation/imprinting machinery genes and environmental exposures. POE-mQTL-modifier interaction models were proposed to test the potential of these factors to modify POE at DNA methylation using data from Generation Scotland: The Scottish Family Health Study(N=2315).

FINDINGS

a set of vulnerable/modifiable POE-CpGs were identified (modifiable-POE-regulated CpGs, N=3). Four factors, 'lifetime smoking status' and 'predicted mRNA expression levels' of TET2, SIRT1 and KDM1A, were found to significantly modify the POE on the three CpGs in both discovery and replication datasets. We further identified plasma protein and health-related phenotypes associated with the methylation level of one of the identified CpGs.

INTERPRETATION

the modifiable POE identified here revealed an important yet indirect path through which genetic background and environmental exposures introduce their effect on DNA methylation, motivating future comprehensive evaluation of the role of these modifiers in complex diseases.

FUNDING

NSFC (81971270),H2020-MSCA-ITN(721815), Wellcome (204979/Z/16/Z,104036/Z/14/Z), MRC (MC_UU_00007/10, MC_PC_U127592696), CSO (CZD/16/6,CZB/4/276, CZB/4/710), SFC (HR03006), EUROSPAN (LSHG-CT-2006-018947), BBSRC (BBS/E/D/30002276), SYSU, Arthritis Research UK, NHLBI, NIH.

Characterisation of age and polarity at onset in bipolar disorder

BACKGROUND

Studying phenotypic and genetic characteristics of age at onset (AAO) and polarity at onset (PAO) in bipolar disorder can provide new insights into disease pathology and facilitate the development of screening tools.

AIMS

To examine the genetic architecture of AAO and PAO and their association with bipolar disorder disease characteristics.

METHOD

Genome-wide association studies (GWASs) and polygenic score (PGS) analyses of AAO (n = 12 977) and PAO (n = 6773) were conducted in patients with bipolar disorder from 34 cohorts and a replication sample (n = 2237). The association of onset with disease characteristics was investigated in two of these cohorts.

RESULTS

Earlier AAO was associated with a higher probability of psychotic symptoms, suicidality, lower educational attainment, not living together and fewer episodes. Depressive onset correlated with suicidality and manic onset correlated with delusions and manic episodes. Systematic differences in AAO between cohorts and continents of origin were observed. This was also reflected in single-nucleotide variant-based heritability estimates, with higher heritabilities for stricter onset definitions. Increased PGS for autism spectrum disorder (β = -0.34 years, s.e. = 0.08), major depression (β = -0.34 years, s.e. = 0.08), schizophrenia (β = -0.39 years, s.e. = 0.08), and educational attainment (β = -0.31 years, s.e. = 0.08) were associated with an earlier AAO. The AAO GWAS identified one significant locus, but this finding did not replicate. Neither GWAS nor PGS analyses yielded significant associations with PAO.

CONCLUSIONS

AAO and PAO are associated with indicators of bipolar disorder severity. Individuals with an earlier onset show an increased polygenic liability for a broad spectrum of psychiatric traits. Systematic differences in AAO across cohorts, continents and phenotype definitions introduce significant heterogeneity, affecting analyses.

Brain structural associations with depression in a large early adolescent sample (the ABCD study®)

BACKGROUND

Depression is the leading cause of disability worldwide with > 50% of cases emerging before the age of 25 years. Large-scale neuroimaging studies in depression implicate robust structural brain differences in the disorder. However, most studies have been conducted in adults and therefore, the temporal origins of depression-related imaging features remain largely unknown. This has important implications for understanding aetiology and informing timings of potential intervention.

METHODS

Here, we examine associations between brain structure (cortical metrics and white matter microstructural integrity) and depression ratings (from caregiver and child), in a large sample (N = 8634) of early adolescents (9 to 11 years old) from the US-based, Adolescent Brain and Cognitive Development (ABCD) Study®. Data was collected from 2016 to 2018.

FINDINGS

We report significantly decreased global cortical and white matter metrics, and regionally in frontal, limbic and temporal areas in adolescent depression (Cohen's d = -0⋅018 to -0⋅041, β = -0·019 to -0⋅057). Further, we report consistently stronger imaging associations for caregiver-reported compared to child-reported depression ratings. Divergences between reports (caregiver vs child) were found to significantly relate to negative socio-environmental factors (e.g., family conflict, absolute β = 0⋅048 to 0⋅169).

INTERPRETATION

Depression ratings in early adolescence were associated with similar imaging findings to those seen in adult depression samples, suggesting neuroanatomical abnormalities may be present early in the disease course, arguing for the importance of early intervention. Associations between socio-environmental factors and reporter discrepancy warrant further consideration, both in the wider context of the assessment of adolescent psychopathology, and in relation to their role in aetiology.

FUNDING

Wellcome Trust (References: 104036/Z/14/Z and 220857/Z/20/Z) and the Medical Research Council (MRC, Reference: MC_PC_17209).

Associations between major psychiatric disorder polygenic risk scores and blood-based markers in UK biobank

Major depressive disorder (MDD), schizophrenia (SCZ), and bipolar disorder (BD) have both shared and discrete genetic risk factors, and are associated with peripheral abnormalities. The relationships between such genetic architectures and blood-based markers are, however, unclear. We investigated relationships between polygenic risk scores (PRS) for these disorders and peripheral markers in the UK Biobank cohort. We calculated polygenic risk scores for n = 367,329 (MDD PRS), n = 366,465 (SCZ PRS), and n = 366,383 (BD PRS) UK Biobank cohort subjects. We then examined associations between disorder PRS and 58 inflammatory/immune, hematological, bone, cardiovascular, hormone, liver, renal and diabetes-associated blood markers using two generalized linear regression models: 'minimally adjusted' controlling for variables such as age and sex, and 'fully adjusted' including additional lifestyle covariates: BMI, alcohol and smoking status, and medication intake. There were 38/58 MDD PRS, 32/58 SCZ PRS, and 20/58 BD PRS-blood marker associations detected for our minimally adjusted model. Of these, 13/38 (MDD PRS), 14/32 (SCZ PRS), and 10/20 (BD PRS) associations remained significant after controlling for lifestyle factors. Many were disorder-specific, with 8/13 unique MDD PRS associations identified. Several disorder-specific associations for MDD and SCZ were immune-related, with mostly positive and negative associations identified for MDD and SCZ PRS respectively. This study suggests that MDD, SCZ and BD have both shared and distinct peripheral markers associated with disorder-specific genetic risk. The results also implicate inflammatory dysfunction in MDD and SCZ, albeit with differences in patterns between the two conditions, and enrich our understanding of potential underlying pathophysiological mechanisms in major psychiatric disorders.

Methylome-wide association study of early life stressors and adult mental health

The environment and events that we are exposed to in utero, during birth and in early childhood influence our future physical and mental health. The underlying mechanisms that lead to these outcomes are unclear, but long-term changes in epigenetic marks, such as DNA methylation, could act as a mediating factor or biomarker. DNA methylation data were assayed at 713 522 CpG sites from 9537 participants of the Generation Scotland: Scottish Family Health Study, a family-based cohort with extensive genetic, medical, family history and lifestyle information. Methylome-wide association studies of eight early life environment phenotypes and two adult mental health phenotypes (major depressive disorder and brief resilience scale) were conducted using DNA methylation data collected from adult whole blood samples. Two genes involved with different developmental pathways (PRICKLE2, Prickle Planar Cell Polarity Protein 2 and ABI1, Abl-Interactor-1) were annotated to CpG sites associated with preterm birth (P < 1.27 × 10⁻⁹). A further two genes important to the development of sensory pathways (SOBP, Sine Oculis Binding Protein Homolog and RPGRIP1, Retinitis Pigmentosa GTPase Regulator Interacting Protein) were annotated to sites associated with low birth weight (P < 4.35 × 10⁻⁸). The examination of methylation profile scores and genes and gene-sets annotated from associated CpGs sites found no evidence of overlap between the early life environment and mental health conditions. Birth date was associated with a significant difference in estimated lymphocyte and neutrophil counts. Previous studies have shown that early life environments influence the risk of developing mental health disorders later in life; however, this study found no evidence that this is mediated by stable changes to the methylome detectable in peripheral blood.

Brain Correlates of Suicide Attempt in 18,925 Participants Across 18 International Cohorts

BACKGROUND

Neuroimaging studies of suicidal behavior have so far been conducted in small samples, prone to biases and false-positive associations, yielding inconsistent results. The ENIGMA-MDD Working Group aims to address the issues of poor replicability and comparability by coordinating harmonized analyses across neuroimaging studies of major depressive disorder and related phenotypes, including suicidal behavior.

METHODS

Here, we pooled data from 18 international cohorts with neuroimaging and clinical measurements in 18,925 participants (12,477 healthy control subjects and 6448 people with depression, of whom 694 had attempted suicide). We compared regional cortical thickness and surface area and measures of subcortical, lateral ventricular, and intracranial volumes between suicide attempters, clinical control subjects (nonattempters with depression), and healthy control subjects.

RESULTS

We identified 25 regions of interest with statistically significant (false discovery rate < .05) differences between groups. Post hoc examinations identified neuroimaging markers associated with suicide attempt including smaller volumes of the left and right thalamus and the right pallidum and lower surface area of the left inferior parietal lobe.

CONCLUSIONS

This study addresses the lack of replicability and consistency in several previously published neuroimaging studies of suicide attempt and further demonstrates the need for well-powered samples and collaborative efforts. Our results highlight the potential involvement of the thalamus, a structure viewed historically as a passive gateway in the brain, and the pallidum, a region linked to reward response and positive affect. Future functional and connectivity studies of suicidal behaviors may focus on understanding how these regions relate to the neurobiological mechanisms of suicide attempt risk.

A genome-wide association study finds genetic variants associated with neck or shoulder pain in UK Biobank

Background

Common types of musculoskeletal conditions include pain in the neck and shoulder areas. This study seeks to identify the genetic variants associated with neck or shoulder pain based on a genome-wide association approach using 203 309 subjects from the UK Biobank cohort and look for replication evidence from the Generation Scotland: Scottish Family Health Study (GS:SFHS) and TwinsUK.

Methods

A genome-wide association study was performed adjusting for age, sex, BMI and nine population principal components. Significant and independent genetic variants were then sent to GS:SFHS and TwinsUK for replication.

Results

We identified three genetic loci that were associated with neck or shoulder pain in the UK Biobank samples. The most significant locus was in an intergenic region in chromosome 17, rs12453010, having P = 1.66 × 10⁻¹¹. The second most significant locus was located in the FOXP2 gene in chromosome 7 with P = 2.38 × 10⁻¹⁰ for rs34291892. The third locus was located in the LINC01572 gene in chromosome 16 with P = 4.50 × 10⁻⁸ for rs62053992. In the replication stage, among four significant and independent genetic variants, rs2049604 in the FOXP2 gene and rs62053992 in the LINC01572 gene were weakly replicated in GS:SFHS (P = 0.0240 and P = 0.0202, respectively).

Conclusions

We have identified three loci associated with neck or shoulder pain in the UK Biobank cohort, two of which were weakly supported in a replication cohort. Further evidence is needed to confirm their roles in neck or shoulder pain.

Genetic and shared couple environmental contributions to smoking and alcohol use in the UK population

Alcohol use and smoking are leading causes of death and disability worldwide. Both genetic and environmental factors have been shown to influence individual differences in the use of these substances. In the present study we tested whether genetic factors, modelled alongside common family environment, explained phenotypic variance in alcohol use and smoking behaviour in the Generation Scotland (GS) family sample of up to 19,377 individuals. SNP and pedigree-associated effects combined explained between 18 and 41% of the variance in substance use. Shared couple effects explained a significant amount of variance across all substance use traits, particularly alcohol intake, for which 38% of the phenotypic variance was explained. We tested whether the within-couple substance use associations were due to assortative mating by testing the association between partner polygenic risk scores in 34,987 couple pairs from the UK Biobank (UKB). No significant association between partner polygenic risk scores were observed. Associations between an individual's alcohol PRS (b = 0.05, S.E. = 0.006, p < 2 × 10-16) and smoking status PRS (b = 0.05, S.E. = 0.005, p < 2 × 10-16) were found with their partner's phenotype. In support of this, G carriers of a functional ADH1B polymorphism (rs1229984), known to be associated with greater alcohol intake, were found to consume less alcohol if they had a partner who carried an A allele at this SNP. Together these results show that the shared couple environment contributes significantly to patterns of substance use. It is unclear whether this is due to shared environmental factors, assortative mating, or indirect genetic effects. Future studies would benefit from longitudinal data and larger sample sizes to assess this further.

Cohort profile for the STratifying Resilience and Depression Longitudinally (STRADL) study: A depression-focused investigation of Generation Scotland, using detailed clinical, cognitive, and neuroimaging assessments

STratifying Resilience and Depression Longitudinally (STRADL) is a population-based study built on the Generation Scotland: Scottish Family Health Study (GS:SFHS) resource. The aim of STRADL is to subtype major depressive disorder (MDD) on the basis of its aetiology, using detailed clinical, cognitive, and brain imaging assessments. The GS:SFHS provides an important opportunity to study complex gene-environment interactions, incorporating linkage to existing datasets and inclusion of early-life variables for two longitudinal birth cohorts. Specifically, data collection in STRADL included: socio-economic and lifestyle variables; physical measures; questionnaire data that assesses resilience, early-life adversity, personality, psychological health, and lifetime history of mood disorder; laboratory samples; cognitive tests; and brain magnetic resonance imaging. Some of the questionnaire and cognitive data were first assessed at the GS:SFHS baseline assessment between 2006-2011, thus providing longitudinal measures relevant to the study of depression, psychological resilience, and cognition. In addition, routinely collected historic NHS data and early-life variables are linked to STRADL data, further providing opportunities for longitudinal analysis. Recruitment has been completed and we consented and tested 1,188 participants.

Epigenome-wide association study and multi-tissue replication of individuals with alcohol use disorder: evidence for abnormal glucocorticoid signaling pathway gene regulation

Alcohol use disorder (AUD) is a chronic debilitating disorder with limited treatment options and poorly defined pathophysiology. There are substantial genetic and epigenetic components; however, the underlying mechanisms contributing to AUD remain largely unknown. We conducted the largest DNA methylation epigenome-wide association study (EWAS) analyses currently available for AUD (total N = 625) and employed a top hit replication (N = 4798) using a cross-tissue/cross-phenotypic approach with the goal of identifying novel epigenetic targets relevant to AUD. Results show that a network of differentially methylated regions in glucocorticoid signaling and inflammation-related genes were associated with alcohol use behaviors. A top probe consistently associated across all cohorts was located in the long non-coding RNA growth arrest specific five gene (GAS5) (p < 10-24). GAS5 has been implicated in regulating transcriptional activity of the glucocorticoid receptor and has multiple functions related to apoptosis, immune function and various cancers. Endophenotypic analyses using peripheral cortisol levels and neuroimaging paradigms showed that methylomic variation in GAS5 network-related probes were associated with stress phenotypes. Postmortem brain analyses documented increased GAS5 expression in the amygdala of individuals with AUD. Our data suggest that alcohol use is associated with differential methylation in the glucocorticoid system that might influence stress and inflammatory reactivity and subsequently risk for AUD.

DNA methylation signatures of aggression and closely related constructs: A meta-analysis of epigenome-wide studies across the lifespan

DNA methylation profiles of aggressive behavior may capture lifetime cumulative effects of genetic, stochastic, and environmental influences associated with aggression. Here, we report the first large meta-analysis of epigenome-wide association studies (EWAS) of aggressive behavior (N = 15,324 participants). In peripheral blood samples of 14,434 participants from 18 cohorts with mean ages ranging from 7 to 68 years, 13 methylation sites were significantly associated with aggression (alpha = 1.2 × 10-7; Bonferroni correction). In cord blood samples of 2425 children from five cohorts with aggression assessed at mean ages ranging from 4 to 7 years, 83% of these sites showed the same direction of association with childhood aggression (r = 0.74, p = 0.006) but no epigenome-wide significant sites were found. Top-sites (48 at a false discovery rate of 5% in the peripheral blood meta-analysis or in a combined meta-analysis of peripheral blood and cord blood) have been associated with chemical exposures, smoking, cognition, metabolic traits, and genetic variation (mQTLs). Three genes whose expression levels were associated with top-sites were previously linked to schizophrenia and general risk tolerance. At six CpGs, DNA methylation variation in blood mirrors variation in the brain. On average 44% (range = 3-82%) of the aggression-methylation association was explained by current and former smoking and BMI. These findings point at loci that are sensitive to chemical exposures with potential implications for neuronal functions. We hope these results to be a starting point for studies leading to applications as peripheral biomarkers and to reveal causal relationships with aggression and related traits.

CRISPR disruption and UK Biobank analysis of a highly conserved polymorphic enhancer suggests a role in male anxiety and ethanol intake

Excessive alcohol intake is associated with 5.9% of global deaths. However, this figure is especially acute in men such that 7.6% of deaths can be attributed to alcohol intake. Previous studies identified a significant interaction between genotypes of the galanin (GAL) gene with anxiety and alcohol abuse in different male populations but were unable to define a mechanism. To address these issues the current study analysed the human UK Biobank cohort and identified a significant interaction (n = 115,865; p = 0.0007) between allelic variation (GG or CA genotypes) in the highly conserved human GAL5.1 enhancer, alcohol intake (AUDIT questionnaire scores) and anxiety in men. Critically, disruption of GAL5.1 in mice using CRISPR genome editing significantly reduced GAL expression in the amygdala and hypothalamus whilst producing a corresponding reduction in ethanol intake in KO mice. Intriguingly, we also found the evidence of reduced anxiety-like behaviour in male GAL5.1KO animals mirroring that seen in humans from our UK Biobank studies. Using bioinformatic analysis and co-transfection studies we further identified the EGR1 transcription factor, that is co-expressed with GAL in amygdala and hypothalamus, as being important in the protein kinase C (PKC) supported activity of the GG genotype of GAL5.1 but less so in the CA genotype. Our unique study uses a novel combination of human association analysis, CRISPR genome editing in mice, animal behavioural analysis and cell culture studies to identify a highly conserved regulatory mechanism linking anxiety and alcohol intake that might contribute to increased susceptibility to anxiety and alcohol abuse in men.

Mental health before and during the COVID-19 pandemic in two longitudinal UK population cohorts

BACKGROUND

The COVID-19 pandemic and mitigation measures are likely to have a marked effect on mental health. It is important to use longitudinal data to improve inferences.

AIMS

To quantify the prevalence of depression, anxiety and mental well-being before and during the COVID-19 pandemic. Also, to identify groups at risk of depression and/or anxiety during the pandemic.

METHOD

Data were from the Avon Longitudinal Study of Parents and Children (ALSPAC) index generation (n = 2850, mean age 28 years) and parent generation (n = 3720, mean age 59 years), and Generation Scotland (n = 4233, mean age 59 years). Depression was measured with the Short Mood and Feelings Questionnaire in ALSPAC and the Patient Health Questionnaire-9 in Generation Scotland. Anxiety and mental well-being were measured with the Generalised Anxiety Disorder Assessment-7 and the Short Warwick Edinburgh Mental Wellbeing Scale.

RESULTS

Depression during the pandemic was similar to pre-pandemic levels in the ALSPAC index generation, but those experiencing anxiety had almost doubled, at 24% (95% CI 23-26%) compared with a pre-pandemic level of 13% (95% CI 12-14%). In both studies, anxiety and depression during the pandemic was greater in younger members, women, those with pre-existing mental/physical health conditions and individuals in socioeconomic adversity, even when controlling for pre-pandemic anxiety and depression.

CONCLUSIONS

These results provide evidence for increased anxiety in young people that is coincident with the pandemic. Specific groups are at elevated risk of depression and anxiety during the COVID-19 pandemic. This is important for planning current mental health provisions and for long-term impact beyond this pandemic.

Polygenic contributions to alcohol use and alcohol use disorders across population-based and clinically ascertained samples

BACKGROUND

Studies suggest that alcohol consumption and alcohol use disorders have distinct genetic backgrounds.

METHODS

We examined whether polygenic risk scores (PRS) for consumption and problem subscales of the Alcohol Use Disorders Identification Test (AUDIT-C, AUDIT-P) in the UK Biobank (UKB; N = 121 630) correlate with alcohol outcomes in four independent samples: an ascertained cohort, the Collaborative Study on the Genetics of Alcoholism (COGA; N = 6850), and population-based cohorts: Avon Longitudinal Study of Parents and Children (ALSPAC; N = 5911), Generation Scotland (GS; N = 17 461), and an independent subset of UKB (N = 245 947). Regression models and survival analyses tested whether the PRS were associated with the alcohol-related outcomes.

RESULTS

In COGA, AUDIT-P PRS was associated with alcohol dependence, AUD symptom count, maximum drinks (R2 = 0.47-0.68%, p = 2.0 × 10-8-1.0 × 10-10), and increased likelihood of onset of alcohol dependence (hazard ratio = 1.15, p = 4.7 × 10-8); AUDIT-C PRS was not an independent predictor of any phenotype. In ALSPAC, the AUDIT-C PRS was associated with alcohol dependence (R2 = 0.96%, p = 4.8 × 10-6). In GS, AUDIT-C PRS was a better predictor of weekly alcohol use (R2 = 0.27%, p = 5.5 × 10-11), while AUDIT-P PRS was more associated with problem drinking (R2 = 0.40%, p = 9.0 × 10-7). Lastly, AUDIT-P PRS was associated with ICD-based alcohol-related disorders in the UKB subset (R2 = 0.18%, p < 2.0 × 10-16).

CONCLUSIONS

AUDIT-P PRS was associated with a range of alcohol-related phenotypes across population-based and ascertained cohorts, while AUDIT-C PRS showed less utility in the ascertained cohort. We show that AUDIT-P is genetically correlated with both use and misuse and demonstrate the influence of ascertainment schemes on PRS analyses.

Sex-stratified genome-wide association study of multisite chronic pain in UK Biobank

Chronic pain is highly prevalent worldwide and imparts a significant socioeconomic and public health burden. Factors influencing susceptibility to, and mechanisms of, chronic pain development, are not fully understood, but sex is thought to play a significant role, and chronic pain is more prevalent in women than in men. To investigate sex differences in chronic pain, we carried out a sex-stratified genome-wide association study of Multisite Chronic Pain (MCP), a derived chronic pain phenotype, in UK Biobank on 178,556 men and 209,093 women, as well as investigating sex-specific genetic correlations with a range of psychiatric, autoimmune and anthropometric phenotypes and the relationship between sex-specific polygenic risk scores for MCP and chronic widespread pain. We also assessed whether MCP-associated genes showed expression pattern enrichment across tissues. A total of 123 SNPs at five independent loci were significantly associated with MCP in men. In women, a total of 286 genome-wide significant SNPs at ten independent loci were discovered. Meta-analysis of sex-stratified GWAS outputs revealed a further 87 independent associated SNPs. Gene-level analyses revealed sex-specific MCP associations, with 31 genes significantly associated in females, 37 genes associated in males, and a single gene, DCC, associated in both sexes. We found evidence for sex-specific pleiotropy and risk for MCP was found to be associated with chronic widespread pain in a sex-differential manner. Male and female MCP were highly genetically correlated, but at an rg of significantly less than 1 (0.92). All 37 male MCP-associated genes and all but one of 31 female MCP-associated genes were found to be expressed in the dorsal root ganglion, and there was a degree of enrichment for expression in sex-specific tissues. Overall, the findings indicate that sex differences in chronic pain exist at the SNP, gene and transcript abundance level, and highlight possible sex-specific pleiotropy for MCP. Results support the proposition of a strong central nervous-system component to chronic pain in both sexes, additionally highlighting a potential role for the DRG and nociception.

Factors associated with sharing e-mail information and mental health survey participation in large population cohorts

BACKGROUND

People who opt to participate in scientific studies tend to be healthier, wealthier and more educated than the broader population. Although selection bias does not always pose a problem for analysing the relationships between exposures and diseases or other outcomes, it can lead to biased effect size estimates. Biased estimates may weaken the utility of genetic findings because the goal is often to make inferences in a new sample (such as in polygenic risk score analysis).

METHODS

We used data from UK Biobank, Generation Scotland and Partners Biobank and conducted phenotypic and genome-wide association analyses on two phenotypes that reflected mental health data availability: (i) whether participants were contactable by e-mail for follow-up; and (ii) whether participants responded to follow-up surveys of mental health.

RESULTS

In UK Biobank, we identified nine genetic loci associated (P <5 × 10-8) with e-mail contact and 25 loci associated with mental health survey completion. Both phenotypes were positively genetically correlated with higher educational attainment and better health and negatively genetically correlated with psychological distress and schizophrenia. One single nucleotide polymorphism association replicated along with the overall direction of effect of all association results.

CONCLUSIONS

Re-contact availability and follow-up participation can act as further genetic filters for data on mental health phenotypes.

Structural brain correlates of serum and epigenetic markers of inflammation in major depressive disorder

Inflammatory processes are implicated in the aetiology of Major Depressive Disorder (MDD); however, the relationship between peripheral inflammation, brain structure and depression remains unclear, partly due to complexities around the use of acute/phasic inflammatory biomarkers. Here, we report the first large-scale study of both serological and methylomic signatures of CRP (considered to represent acute and chronic measures of inflammation respectively) and their associations with depression status/symptoms, and structural neuroimaging phenotypes (T1 and diffusion MRI) in a large community-based sample (Generation Scotland; NMDD cases = 271, Ncontrols = 609). Serum CRP was associated with overall MDD severity, and specifically with current somatic symptoms- general interest (β = 0.145, PFDR = 6 × 10-4) and energy levels (β = 0.101, PFDR = 0.027), along with reduced entorhinal cortex thickness (β = -0.095, PFDR = 0.037). DNAm CRP was significantly associated with reduced global grey matter/cortical volume and widespread reductions in integrity of 16/24 white matter tracts (with greatest regional effects in the external and internal capsules, βFA= -0.12 to -0.14). In general, the methylation-based measures showed stronger associations with imaging metrics than serum-based CRP measures (βaverage = -0.15 versus βaverage = 0.01 respectively). These findings provide evidence for central effects of peripheral inflammation from both serological and epigenetic markers of inflammation, including in brain regions previously implicated in depression. This suggests that these imaging measures may be involved in the relationship between peripheral inflammation and somatic/depressive symptoms. Notably, greater effects on brain morphology were seen for methylation-based rather than serum-based measures of inflammation, indicating the importance of such measures for future studies.

Grey and white matter associations of psychotic-like experiences in a general population sample (UK Biobank)

There has been a substantial amount of research reporting the neuroanatomical associations of psychotic symptoms in people with schizophrenia. Comparatively little attention has been paid to the neuroimaging correlates of subclinical psychotic symptoms, so-called "psychotic-like experiences" (PLEs), within large healthy populations. PLEs are relatively common in the general population (7-13%), can be distressing and negatively affect health. This study therefore examined gray and white matter associations of four different PLEs (auditory or visual PLEs, and delusional ideas about conspiracies or communications) in subjects of the UK Biobank study with neuroimaging data (N = 21,390, mean age = 63 years). We tested for associations between any PLE (N = 768) and individual PLEs with gray and white matter brain structures, controlling for sex, age, intracranial volume, scanning site, and position in the scanner. Individuals that reported having experienced auditory hallucinations (N = 272) were found to have smaller volumes of the caudate, putamen, and accumbens (β = -0.115-0.134, pcorrected = 0.048-0.036), and reduced temporal lobe volume (β = -0.017, pcorrected = 0.047) compared to those that did not. People who indicated that they had ever believed in unreal conspiracies (N = 111) had a larger volume of the left amygdala (β = 0.023, pcorrected = 0.038). Individuals that reported a history of visual PLEs (N = 435) were found to have reduced white matter microstructure of the forceps major (β = -0.029, pcorrected = 0.009), an effect that was more marked in participants who reported PLEs as distressing. These associations were not accounted for by diagnoses of psychotic or depressive illness, nor the known risk factors for psychotic symptoms of childhood adversity or cannabis use. These findings suggest altered regional gray matter volumes and white matter microstructure in association with PLEs in the general population. They further suggest that these alterations may appear more frequently with the presentation of different psychotic symptoms in the absence of clinically diagnosed psychotic disorders.

Potential Genetic Overlap Between Insomnia and Sleep Symptoms in Major Depressive Disorder: A Polygenic Risk Score Analysis

Background: The prevalence of insomnia and hypersomnia in depressed individuals is substantially higher than that found in the general population. Unfortunately, these concurrent sleep problems can have profound effects on the disease course. Although the full biology of sleep remains to be elucidated, a recent genome-wide association (GWAS) of insomnia, and other sleep traits in over 1 million individuals was recently published and provides many promising hits for genetics of insomnia in a population-based sample. Methods: Using data from the largest available GWAS of insomnia and other sleep traits, we sought to test if sleep variable PRS scores derived from population-based studies predicted sleep variables in samples of depressed cases [Psychiatric Genomics Consortium - Major Depressive Disorder subjects (PGC MDD)]. A leave-one-out analysis was performed to determine the effects that each individual study had on our results. Results: The only significant finding was for insomnia, where p-value threshold, p = 0.05 was associated with insomnia in our PGC MDD sample (R ² = 1.75^-3, p = 0.006). Conclusion: Our results reveal that <1% of variance is explained by the variants that cover the two significant p-value thresholds, which is in line with the fact that depression and insomnia are both polygenic disorders. To the best of our knowledge, this is the first study to investigate genetic overlap between the general population and a depression sample for insomnia, which has important treatment implications, such as leading to novel drug targets in future research efforts.

Educational attainment impacts drinking behaviors and risk for alcohol dependence: results from a two-sample Mendelian randomization study with ~780,000 participants

Observational studies suggest that lower educational attainment (EA) may be associated with risky alcohol use behaviors; however, these findings may be biased by confounding and reverse causality. We performed two-sample Mendelian randomization (MR) using summary statistics from recent genome-wide association studies (GWAS) with >780,000 participants to assess the causal effects of EA on alcohol use behaviors and alcohol dependence (AD). Fifty-three independent genome-wide significant SNPs previously associated with EA were tested for association with alcohol use behaviors. We show that while genetic instruments associated with increased EA are not associated with total amount of weekly drinks, they are associated with reduced frequency of binge drinking ≥6 drinks (ß_IVW = -0.198, 95% CI, -0.297 to -0.099, P_IVW = 9.14 × 10^-5), reduced total drinks consumed per drinking day (ß_IVW = -0.207, 95% CI, -0.293 to -0.120, P_IVW = 2.87 × 10^-6), as well as lower weekly distilled spirits intake (ß_IVW = -0.148, 95% CI, -0.188 to -0.107, P_IVW = 6.24 × 10^-13). Conversely, genetic instruments for increased EA were associated with increased alcohol intake frequency (ß_IVW = 0.331, 95% CI, 0.267-0.396, P_IVW = 4.62 × 10^-24), and increased weekly white wine (ß_IVW = 0.199, 95% CI, 0.159-0.238, P_IVW = 7.96 × 10^-23) and red wine intake (ß_IVW = 0.204, 95% CI, 0.161-0.248, P_IVW = 6.67 × 10^-20). Genetic instruments associated with increased EA reduced AD risk: an additional 3.61 years schooling reduced the risk by ~50% (OR_IVW = 0.508, 95% CI, 0.315-0.819, P_IVW = 5.52 × 10^-3). Consistency of results across complementary MR methods accommodating different assumptions about genetic pleiotropy strengthened causal inference. Our findings suggest EA may have important effects on alcohol consumption patterns and may provide potential mechanisms explaining reported associations between EA and adverse health outcomes.