Search

Genetics have nominated many schizophrenia risk genes and identified convergent signals between schizophrenia and neurodevelopmental disorders. However, functional interpretation of the nominated genes in the relevant brain cell types is often lacking. We executed interaction proteomics for six schizophrenia risk genes that have also been implicated in neurodevelopment in human induced cortical neurons. The resulting protein network is enriched for common variant risk of schizophrenia in Europeans and East Asians, is down-regulated in layer 5/6 cortical neurons of individuals affected by schizophrenia, and can complement fine-mapping and eQTL data to prioritize additional genes in GWAS loci. A sub-network centered on HCN1 is enriched for common variant risk and contains proteins (HCN4 and AKAP11) enriched for rare protein-truncating mutations in individuals with schizophrenia and bipolar disorder. Our findings showcase brain cell-type-specific interactomes as an organizing framework to facilitate interpretation of genetic and transcriptomic data in schizophrenia and its related disorders.

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.

The Mood Disorder Questionnaire (MDQ) is a common screening tool for bipolar disorder that assesses manic symptoms. Its utility for genetic studies of mania or bipolar traits has not been fully examined. We psychometrically compared the MDQ to self-reported bipolar disorder in participants from the United Kingdom National Institute of Health and Care Research Mental Health BioResource. We conducted genome-wide association studies of manic symptom quantitative traits and symptom subgroups, derived from the MDQ items (N = 11,568–19,859). We calculated genetic correlations with bipolar disorder and other psychiatric and behavioral traits. The MDQ screener showed low positive predictive value (0.29) for self-reported bipolar disorder. Neither concurrent nor lifetime manic symptoms were genetically correlated with bipolar disorder. Lifetime manic symptoms had a highest genetic correlation (rg = 1.0) with posttraumatic stress disorder although this was not confirmed by within-cohort phenotypic correlations (rp = 0.41). Other significant genetic correlations included attention deficit hyperactivity disorder (rg = 0.69), insomnia (rg = 0.55), and major depressive disorder (rg = 0.42). Our study adds to existing literature questioning the MDQ's validity and suggests it may capture symptoms of general distress or psychopathology, rather than hypomania/mania specifically, in at-risk populations.

Contains fulltext : 291574.pdf (Publisher’s version ) (Open Access) Left-right asymmetry is an important organizing feature of the healthy brain that may be altered in schizophrenia, but most studies have used relatively small samples and heterogeneous approaches, resulting in equivocal findings. We carried out the largest case-control study of structural brain asymmetries in schizophrenia, with MRI data from 5,080 affected individuals and 6,015 controls across 46 datasets, using a single image analysis protocol. Asymmetry indexes were calculated for global and regional cortical thickness, surface area, and subcortical volume measures. Differences of asymmetry were calculated between affected individuals and controls per dataset, and effect sizes were meta-analyzed across datasets. Small average case-control differences were observed for thickness asymmetries of the rostral anterior cingulate and the middle temporal gyrus, both driven by thinner left-hemispheric cortices in schizophrenia. Analyses of these asymmetries with respect to the use of antipsychotic medication and other clinical variables did not show any significant associations. Assessment of age- and sex-specific effects revealed a stronger average leftward asymmetry of pallidum volume between older cases and controls. Case-control differences in a multivariate context were assessed in a subset of the data (N = 2,029), which revealed that 7% of the variance across all structural asymmetries was explained by case-control status. Subtle case-control differences of brain macrostructural asymmetry may reflect differences at the molecular, cytoarchitectonic, or circuit levels that have functional relevance for the disorder. Reduced left middle temporal cortical thickness is consistent with altered left-hemisphere language network organization in schizophrenia.

BackgroundCoded healthcare data may not capture all stroke cases and has limited accuracy for stroke subtypes. We sought to determine the incremental value of adding natural language processing (NLP) of free-text radiology reports to international classification of disease (ICD-10) codes to phenotype stroke, and stroke subtypes, in routinely collected healthcare datasets.MethodsWe linked participants in a community-based prospective cohort study, Generation Scotland, to clinical brain imaging reports (2008-2020) from five Scottish health boards. We used five combinations of NLP outputs and ICD-10 codes to define stroke phenotypes. With these phenotype models we measured the: stroke incidence standardised to a European Standardised Population; adjusted hazard ratio (aHR) of baseline hypertension for later stroke; and proportion of participants allocated stroke subtypes.ResultsOf 19,026 participants, over a mean follow-up of 10.2 years, 1938 had 3493 brain scans. Any stroke was identified in 534 participants: 319 with NLP alone, 59 with ICD-10 codes alone and 156 with both ICD-10 codes and an NLP report consistent with stroke. The stroke aHR for baseline hypertension was 1.47 (95%CI: 1.12-1.92) for NLP-defined stroke only; 1.57 (95%CI: 1.18-2.10) for ICD-10 defined stroke only; and 1.81 (95%CI: 1.20-2.72) for cases with ICD 10 stroke codes and NLP stroke phenotypes. The age-standardised incidence of stroke for these phenotype models was 1.35, 1.34, and 0.65 per 1000 person years, respectively. The proportion of strokes not subtyped was 26% (57/215) using only ICD-10, 9% (42/467) using only NLP, and 12% (65/534) using both NLP and ICD-10.ConclusionsAddition of NLP derived phenotypes to ICD-10 stroke codes identified approximately 2.5 times more stroke cases and greatly increased the proportion with subtyping. The phenotype model using ICD 10 stroke codes and NLP stroke phenotypes had the strongest association with baseline hypertension. This information is relevant to large cohort studies and clinical trials that use routine electronic health records for outcome ascertainment.

Observational studies and randomized controlled trials presented inconsistent findings on the effects of cholesterol-lowering statins on depression. It therefore remains unclear whether statins have any beneficial effects on depression, and if so, what the underlying molecular mechanisms are. Here, we aimed to use genomic approaches to investigate this further. Using Connectivity Map (CMap), we first investigated whether statins and antidepressants shared pharmacological effects by interrogating gene expression responses to drug exposure in human cell lines. Second, using Mendelian randomization analysis, we investigated both on-target (through HMGCR inhibition) and potential off-target (through ITGAL and HDAC2 inhibition) causal effects of statins on depression risk and depressive symptoms, and traits related to the shared biological pathways identified from CMap analysis. Compounds inducing highly similar gene expression responses to statins in HA1E cells (indicated by an average connectivity score with statins > 90) were found to be enriched for antidepressants (12 out of 38 antidepressants; p = 9E-08). Genes perturbed in the same direction by both statins and antidepressants were significantly enriched for diverse cellular and metabolic pathways, and various immune activation, development and response processes. MR analysis did not identify any significant associations between statin exposure and depression risk or symptoms after multiple testing correction. However, genetically proxied HMGCR inhibition was strongly associated with alterations in platelets (a prominent serotonin reservoir) and monocyte percentage, which have previously been implicated in depression. Genetically proxied ITGAL inhibition was strongly associated with basophil, monocyte and neutrophil counts. We identified biological pathways that are commonly perturbed by both statins and antidepressants, and haematological biomarkers genetically associated with statin targets. Our findings warrant pre-clinical investigation of the causal role of these shared pathways in depression and potential as therapeutic targets, and investigation of whether blood biomarkers may be important considerations in clinical trials investigating effects of statins on depression.

We report results from the Bipolar Exome (BipEx) collaboration analysis of whole-exome sequencing of 13,933 patients with bipolar disorder (BD) matched with 14,422 controls. We find an excess of ultra-rare protein-truncating variants (PTVs) in patients with BD among genes under strong evolutionary constraint in both major BD subtypes. We find enrichment of ultra-rare PTVs within genes implicated from a recent schizophrenia exome meta-analysis (SCHEMA; 24,248 cases and 97,322 controls) and among binding targets of CHD8. Genes implicated from genome-wide association studies (GWASs) of BD, however, are not significantly enriched for ultra-rare PTVs. Combining gene-level results with SCHEMA, AKAP11 emerges as a definitive risk gene (odds ratio (OR) = 7.06, P = 2.83 × 10−9). At the protein level, AKAP-11 interacts with GSK3B, the hypothesized target of lithium, a primary treatment for BD. Our results lend support to BD’s polygenicity, demonstrating a role for rare coding variation as a significant risk factor in BD etiology.

Inflammation and ageing-related DNA methylation patterns in the blood have been linked to a variety of morbidities, including cognitive decline and neurodegenerative disease. However, it is unclear how these blood-based patterns relate to patterns within the brain, and how each associates with central cellular profiles. In this study, we profiled DNA methylation in both the blood and in five post-mortem brain regions (BA17, BA20/21, BA24, BA46 and hippocampus) in 14 individuals from the Lothian Birth Cohort 1936. Microglial burdens were additionally quantified in the same brain regions. DNA methylation signatures of five epigenetic ageing biomarkers (‘epigenetic clocks’), and two inflammatory biomarkers (DNA methylation proxies for C-reactive protein and interleukin-6) were compared across tissues and regions. Divergent correlations between the inflammation and ageing signatures in the blood and brain were identified, depending on region assessed. Four out of the five assessed epigenetic age acceleration measures were found to be highest in the hippocampus (β range=0.83-1.14, p≤0.02). The inflammation-related DNA methylation signatures showed no clear variation across brain regions. Reactive microglial burdens were found to be highest in the hippocampus (β=1.32, p=5×10-4); however, the only association identified between the blood- and brain-based methylation signatures and microglia was a significant positive association with acceleration of one epigenetic clock (termed DNAm PhenoAge) averaged over all five brain regions (β=0.40, p=0.002). This work highlights a potential vulnerability of the hippocampus to epigenetic ageing and provides preliminary evidence of a relationship between DNA methylation signatures in the brain and differences in microglial burdens.

Contains fulltext : 248376.pdf (Publisher’s version ) (Open Access) For many traits, males show greater variability than females, with possible implications for understanding sex differences in health and disease. Here, the ENIGMA (Enhancing Neuro Imaging Genetics through Meta-Analysis) Consortium presents the largest-ever mega-analysis of sex differences in variability of brain structure, based on international data spanning nine decades of life. Subcortical volumes, cortical surface area and cortical thickness were assessed in MRI data of 16,683 healthy individuals 1-90 years old (47% females). We observed significant patterns of greater male than female between-subject variance for all subcortical volumetric measures, all cortical surface area measures, and 60% of cortical thickness measures. This pattern was stable across the lifespan for 50% of the subcortical structures, 70% of the regional area measures, and nearly all regions for thickness. Our findings that these sex differences are present in childhood implicate early life genetic or gene-environment interaction mechanisms. The findings highlight the importance of individual differences within the sexes, that may underpin sex-specific vulnerability to disorders.

Special Issue: The ENIGMA Consortium: the first 10 years., Delineating the association of age and cortical thickness in healthy individuals is critical given the association of cortical thickness with cognition and behavior. Previous research has shown that robust estimates of the association between age and brain morphometry require large-scale studies. In response, we used cross-sectional data from 17,075 individuals aged 3–90 years from the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) Consortium to infer age-related changes in cortical thickness. We used fractional polynomial (FP) regression to quantify the association between age and cortical thickness, and we computed normalized growth centiles using the parametric Lambda, Mu, and Sigma method. Interindividual variability was estimated using meta-analysis and one-way analysis of variance. For most regions, their highest cortical thickness value was observed in childhood. Age and cortical thickness showed a negative association; the slope was steeper up to the third decade of life and more gradual thereafter; notable exceptions to this general pattern were entorhinal, temporopolar, and anterior cingulate cortices. Interindividual variability was largest in temporal and frontal regions across the lifespan. Age and its FP combinations explained up to 59% variance in cortical thickness. These results may form the basis of further investigation on normative deviation in cortical thickness and its significance for behavioral and cognitive outcomes., European Community's Seventh Framework Programme, Grant/Award Numbers: 278948, 602450, 603016, 602805; US National Institute of Child Health and Human Development, Grant/Award Numbers: RO1HD050735, 1009064, 496682; QIMR Berghofer Medical Research Institute and the Centre for Advanced Imaging, University of Queensland; ICTSI NIH/NCRR, Grant/Award Number: RR025761; European Community's Horizon 2020 Programme, Grant/Award Numbers: 667302, 643051; Vici Innovation Program, Grant/Award Numbers: #91619115, 016-130-669; NWO Brain & Cognition Excellence Program, Grant/Award Number: 433-09-229; Biobanking and Biomolecular Resources Research Infrastructure (Netherlands) (BBMRI-NL); Spinozapremie, Grant/Award Number: NWO-56-464-14192; Biobanking and Biomolecular Resources Research Infrastructure, Grant/Award Numbers: 184.033.111, 184.021.007; Netherlands Organization for Health Research and Development (ZonMW), Grant/Award Numbers: 480-15-001/674, 024.001.003, 911-09-032, 056-32-010, 481-08-011, 016-115-035, 31160008, 400-07-080, 400-05-717, 451-04-034, 463-06-001, 480-04-004, 904-61-193, 912-10-020, 985-10-002, 904-61-090; NIMH, Grant/Award Number: R01 MH090553; Geestkracht programme of the Dutch Health Research Council, Grant/Award Number: 10-000-1001; FP7 Ideas: European Research Council; Nederlandse Organisatie voor Wetenschappelijk Onderzoek, Grant/Award Numbers: NWO/SPI 56-464-14192, NWO-MagW 480-04-004, 433-09-220, NWO 51.02.062, NWO 51.02.061; National Center for Advancing Translational Sciences, National Institutes of Health, Grant/Award Number: UL1 TR000153; National Center for Research Resources; National Center for Research Resources at the National Institutes of Health, Grant/Award Numbers: NIH 1U24 RR025736-01, NIH 1U24 RR021992; NIH Institutes contributing to the Big Data to Knowledge; U.S. National Institutes of Health, Grant/Award Numbers: R01 CA101318, P30 AG10133, R01 AG19771; Medical Research Council, Grant/Award Numbers: U54EB020403, G0500092; National Institute of Mental Health, Grant/Award Numbers: R01MH117014, R01MH042191; Fundación Instituto de Investigación Marqués de Valdecilla, Grant/Award Numbers: API07/011, NCT02534363, NCT0235832; Instituto de Salud Carlos III, Grant/Award Numbers: PI14/00918, PI14/00639, PI060507, PI050427, PI020499; Swedish Research Council, Grant/Award Numbers: 523-2014-3467, 2017-00949, 521-2014-3487; South-Eastern Norway Health Authority; the Research Council of Norway, Grant/Award Number: 223273; South Eastern Norway Regional Health Authority, Grant/Award Numbers: 2017-112, 2019107; Icahn School of Medicine at Mount Sinai; Seventh Framework Programme (FP7/2007-2013), Grant/Award Number: 602450; National Institutes of Health, Grant/Award Numbers: R01 MH116147, R01 MH113619, R01 MH104284; South London and Maudsley NHS Foundation Trust; the National Institute for Health Research (NIHR)

Background Epigenetic clocks can track both chronological age (cAge) and biological age (bAge). The latter is typically defined by physiological biomarkers and risk of adverse health outcomes, including all-cause mortality. As cohort sample sizes increase, estimates of cAge and bAge become more precise. Here, we aim to develop accurate epigenetic predictors of cAge and bAge, whilst improving our understanding of their epigenomic architecture. Methods First, we perform large-scale (N = 18,413) epigenome-wide association studies (EWAS) of chronological age and all-cause mortality. Next, to create a cAge predictor, we use methylation data from 24,674 participants from the Generation Scotland study, the Lothian Birth Cohorts (LBC) of 1921 and 1936, and 8 other cohorts with publicly available data. In addition, we train a predictor of time to all-cause mortality as a proxy for bAge using the Generation Scotland cohort (1214 observed deaths). For this purpose, we use epigenetic surrogates (EpiScores) for 109 plasma proteins and the 8 component parts of GrimAge, one of the current best epigenetic predictors of survival. We test this bAge predictor in four external cohorts (LBC1921, LBC1936, the Framingham Heart Study and the Women’s Health Initiative study). Results Through the inclusion of linear and non-linear age-CpG associations from the EWAS, feature pre-selection in advance of elastic net regression, and a leave-one-cohort-out (LOCO) cross-validation framework, we obtain cAge prediction with a median absolute error equal to 2.3 years. Our bAge predictor was found to slightly outperform GrimAge in terms of the strength of its association to survival (HRGrimAge = 1.47 [1.40, 1.54] with p = 1.08 × 10−52, and HRbAge = 1.52 [1.44, 1.59] with p = 2.20 × 10−60). Finally, we introduce MethylBrowsR, an online tool to visualise epigenome-wide CpG-age associations. Conclusions The integration of multiple large datasets, EpiScores, non-linear DNAm effects, and new approaches to feature selection has facilitated improvements to the blood-based epigenetic prediction of biological and chronological age.

BackgroundAnxiety and depressive disorders often co-occur and the order of their emergence may be associated with different clinical outcomes. However, minimal research has been conducted on anxiety-anxiety comorbidity. This study examined factors associated with anxiety comorbidity and anxiety-MDD temporal sequence.MethodsOnline, self-report data were collected from the UK-based GLAD and COPING NBR cohorts (N = 38,775). Logistic regression analyses compared differences in sociodemographic, trauma, and clinical factors between single anxiety, anxiety-anxiety comorbidity, anxiety-MDD (major depressive disorder) comorbidity, and MDD-only. Additionally, anxiety-first and MDD-first anxiety-MDD were compared. Differences in familial risk were assessed in those participants with self-reported family history or genotype data.ResultsAnxiety-anxiety and anxiety-MDD had higher rates of self-reported anxiety or depressive disorder diagnoses, younger age of onset, and higher recurrence than single anxiety. Anxiety-MDD displayed greater clinical severity/complexity than MDD only. Anxiety-anxiety had more severe current anxiety symptoms, less severe current depressive symptoms, and reduced likelihood of self-reporting an anxiety/depressive disorder diagnosis than anxiety-MDD. Anxiety-first anxiety-MDD had a younger age of onset, more severe anxiety symptoms, and less likelihood of self-reporting a diagnosis than MDD-first. Minimal differences in familial risk were found.LimitationsSelf-report, retrospective measures may introduce recall bias. The familial risk analyses were likely underpowered.ConclusionsAnxiety-anxiety comorbidity displayed a similarly severe and complex profile of symptoms as anxiety-MDD but distinct features. For anxiety-MDD, first-onset anxiety had an earlier age of onset and greater severity than MDD-first. Anxiety disorders and comorbidity warrant further investigation and attention in research and practice.

BACKGROUND: Metabolic differences have been reported between individuals with and without major depressive disorder (MDD), but their consistency and causal relevance have been unclear.METHODS: We conducted a metabolome-wide association study of MDD with 249 metabolomic measures available in the UK Biobank (n = 29,757). We then applied two-sample bidirectional Mendelian randomization and colocalization analysis to identify potentially causal relationships between each metabolite and MDD.RESULTS: A total of 191 metabolites tested were significantly associated with MDD (false discovery rate-corrected p < .05), which decreased to 129 after adjustment for likely confounders. Lower abundance of omega-3 fatty acid measures and a higher omega-6 to omega-3 ratio showed potentially causal effects on liability to MDD. There was no evidence of a causal effect of MDD on metabolite levels. Furthermore, genetic signals associated with docosahexaenoic acid colocalized with loci associated with MDD within the fatty acid desaturase gene cluster. Post hoc Mendelian randomization of gene-transcript abundance within the fatty acid desaturase cluster demonstrated a potentially causal association with MDD. In contrast, colocalization analysis did not suggest a single causal variant for both transcript abundance and MDD liability, but rather the likely existence of two variants in linkage disequilibrium with one another.CONCLUSIONS: Our findings suggest that decreased docosahexaenoic acid and increased omega-6 to omega-3 fatty acids ratio may be causally related to MDD. These findings provide further support for the causal involvement of fatty acids in MDD.

Introduction: Childhood trauma and adversity are common across societies and have strong associations with physical and psychiatric morbidity throughout the life-course. One mechanism through which childhood trauma may predispose individuals to poor psychiatric outcomes, such as raised risk of lifetime depression, could be via associations with brain structure. This study aimed to elucidate the associations between childhood trauma scores and brain structure across two large, independent community cohorts. Methods: The two samples comprised (i) a subsample of individuals from Generation Scotland with imaging and in-depth phenotyping, including the CTQ-28 (n=1,024); and (ii) individuals from UK Biobank with imaging and a modified summary CTQ measure (n=27,202). This comprised n=28,226 for mega-analysis. Scans were processed using FreeSurfer image processing software, providing cortical and subcortical as well as global brain metrics. Regression models were used to determine associations between these metrics and childhood trauma measures. Associations between childhood trauma measures and psychiatric phenotypes were also explored. Results: Childhood trauma measures associated with lifetime risk of depression diagnosis with similar ORs across cohorts (OR 1.06, 1.23 GS and UKB respectively), which also related to earlier onset and more recurrent course within both samples. There was also evidence for associations between childhood trauma measures and a range of brain structures. Replicated findings included reduced global brain volumes, reduced cortical surface area but not thickness, with highest effects at mega-analysis seen in the frontal (β=-0.0385, SE=0.0048, p(FDR)=5.43×10 -15) and parietal lobes (β=-36 0.0387, SE=0.005, p(FDR)=1.56×10 -14). At a regional level, one subcortical regional volume in particular - the ventral diencephalon (VDc) - displayed significant associations with childhood trauma measures across the two cohorts and at mega-analysis (β=-0.0232, SE=0.0039, p(FDR)=2.91×10 -8). There was also evidence for associations with reduced hippocampus, thalamus, and nucleus accumbens volumes, however these were not as consistent across cohorts. Discussion: There was strong evidence for associations between childhood trauma and reduced global and regional brain volumes across cohorts. In particular, the presence of an association between childhood trauma and the volume of the VDc (which includes the hypothalamic area), with replication, provides further evidence of the importance of neuroendocrine stress response pathways in links between early life stress and clinical outcomes.

Background Efforts to develop neuroimaging-based biomarkers in major depressive disorder (MDD), at the individual level, have been limited to date. As diagnostic criteria are currently symptom-based, MDD is conceptualized as a disorder rather than a disease with a known etiology; further, neural measures are often confounded by medication status and heterogeneous symptom states. Methods We describe a consortium to quantify neuroanatomical and neurofunctional heterogeneity via the dimensions of novel multivariate coordinate system (COORDINATE-MDD). Utilizing imaging harmonization and machine learning methods in a large cohort of medication-free, deeply phenotyped MDD participants, patterns of brain alteration are defined in replicable and neurobiologically-based dimensions and offer the potential to predict treatment response at the individual level. International datasets are being shared from multi-ethnic community populations, first episode and recurrent MDD, which are medication-free, in a current depressive episode with prospective longitudinal treatment outcomes and in remission. Neuroimaging data consist of de-identified, individual, structural MRI and resting-state functional MRI with additional positron emission tomography (PET) data at specific sites. State-of-the-art analytic methods include automated image processing for extraction of anatomical and functional imaging variables, statistical harmonization of imaging variables to account for site and scanner variations, and semi-supervised machine learning methods that identify dominant patterns associated with MDD from neural structure and function in healthy participants. Results We are applying an iterative process by defining the neural dimensions that characterise deeply phenotyped samples and then testing the dimensions in novel samples to assess specificity and reliability. Crucially, we aim to use machine learning methods to identify novel predictors of treatment response based on prospective longitudinal treatment outcome data, and we can externally validate the dimensions in fully independent sites. Conclusion We describe the consortium, imaging protocols and analytics using preliminary results. Our findings thus far demonstrate how datasets across many sites can be harmonized and constructively pooled to enable execution of this large-scale project.

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.

BackgroundBlood DNA methylation can inform us about the biological mechanisms that underlie common disease states. Previous epigenome-wide analyses of common diseases often focus solely on the prevalence or incidence of individual conditions and rely on small sample sizes, which may limit power to discover disease-associated loci.ResultsWe conduct blood-based epigenome-wide association studies on the prevalence of 14 common disease states in Generation Scotland (nindividuals≤18,413, nCpGs=752,722). We also utilise health record linkage to perform epigenome-wide analyses on the incidence of 19 disease states. We present a structured literature review on existing epigenome-wide analyses for all 19 disease states to assess the degree of replication within the existing literature and the novelty of the present findings.We identify 69 associations between CpGs and the prevalence of four disease states at baseline, of which 58 are novel. We also uncover 64 CpGs that associate with the incidence of two disease states (COPD and type 2 diabetes), of which 56 are novel. These associations were independent from common lifestyle risk factors. We highlight poor replication across the existing literature. Here, replication was defined by the reporting of at least one common gene in >2 studies examining the same disease state. Existing blood-based epigenome-wide analyses showed evidence of replication for only 4/19 disease states (with up-to-15% of unique genes replicated for lung cancer).ConclusionsOur summary data and structured review of the literature provide an important platform to guide future studies that examine the role of blood DNA methylation in complex disease states.

Objective: A major limitation of current suicide research is the lack of power to identify robust correlates of suicidal thoughts or behavior. Variation in suicide risk assessment instruments used across cohorts may represent a limitation to pooling data in international consortia. Method: Here, we examine this issue through two approaches: (a) an extensive literature search on the reliability and concurrent validity of the most commonly used instruments and (b) by pooling data (N ∼ 6,000 participants) from cohorts from the Enhancing NeuroImaging Genetics Through Meta-Analysis (ENIGMA) Major Depressive Disorder and ENIGMA–Suicidal Thoughts and Behaviour working groups, to assess the concurrent validity of instruments currently used for assessing suicidal thoughts or behavior. Results: We observed moderate-to-high correlations between measures, consistent with the wide range (κ range: 0.15–0.97; r range: 0.21–0.94) reported in the literature. Two common multi-item instruments, the Columbia Suicide Severity Rating Scale and the Beck Scale for Suicidal Ideation were highly correlated with each other (r = 0.83). Sensitivity analyses identified sources of heterogeneity such as the time frame of the instrument and whether it relies on self-report or a clinical interview. Finally, construct-specific analyses suggest that suicide ideation items from common psychiatric questionnaires are most concordant with the suicide ideation construct of multi-item instruments. Conclusions: Our findings suggest that multi-item instruments provide valuable information on different aspects of suicidal thoughts or behavior but share a modest core factor with single suicidal ideation items. Retrospective, multisite collaborations including distinct instruments should be feasible provided they harmonize across instruments or focus on specific constructs of suicidality.

Alcohol misuse is common in many societies worldwide and is associated with extensive morbidity and mortality, often leading to alcohol use disorders (AUD) and alcohol-related end-organ damage. The underlying mechanisms contributing to the development of AUD are largely unknown; however, growing evidence suggests that alcohol consumption is strongly associated with alterations in DNA methylation. Identification of alcohol-associated methylomic variation might provide novel insights into pathophysiology and novel treatment targets for AUD. Here we performed the largest single-cohort epigenome-wide association study (EWAS) of alcohol consumption to date (N = 8161) and cross-validated findings in AUD populations with relevant endophenotypes, as well as alcohol-related animal models. Results showed 2504 CpGs significantly associated with alcohol consumption (Bonferroni p value < 6.8 × 10−8) with the five leading probes located in SLC7A11 (p = 7.75 × 10−108), JDP2 (p = 1.44 × 10−56), GAS5 (p = 2.71 × 10−47), TRA2B (p = 3.54 × 10−42), and SLC43A1 (p = 1.18 × 10−40). Genes annotated to associated CpG sites are implicated in liver and brain function, the cellular response to alcohol and alcohol-associated diseases, including hypertension and Alzheimer’s disease. Two-sample Mendelian randomization confirmed the causal relationship of consumption on AUD risk (inverse variance weighted (IVW) p = 5.37 × 10−09). A methylation-based predictor of alcohol consumption was able to discriminate AUD cases in two independent cohorts (p = 6.32 × 10−38 and p = 5.41 × 10−14). The top EWAS probe cg06690548, located in the cystine/glutamate transporter SLC7A11, was replicated in an independent cohort of AUD and control participants (N = 615) and showed strong hypomethylation in AUD (p −17). Decreased CpG methylation at this probe was consistently associated with clinical measures including increased heavy drinking days (p −4), increased liver function enzymes (GGT (p = 1.03 × 10−21), ALT (p = 1.29 × 10−6), and AST (p = 1.97 × 10−8)) in individuals with AUD. Postmortem brain analyses documented increased SLC7A11 expression in the frontal cortex of individuals with AUD and animal models showed marked increased expression in liver, suggesting a mechanism by which alcohol leads to hypomethylation-induced overexpression of SLC7A11. Taken together, our EWAS discovery sample and subsequent validation of the top probe in AUD suggest a strong role of abnormal glutamate signaling mediated by methylomic variation in SLC7A11. Our data are intriguing given the prominent role of glutamate signaling in brain and liver and might provide an important target for therapeutic intervention.

A complex interplay of genetic and environmental risk factors influence global brain structural alterations associated with brain health and disease. Epigenome-wide association studies (EWAS) of global brain imaging phenotypes have the potential to reveal the mechanisms of brain health and disease and can lead to better predictive analytics through the development of risk scores.We perform an EWAS of global brain volumes in Generation Scotland using peripherally measured whole blood DNA methylation (DNAm) from two assessments, (i) at baseline recruitment, ~6 years prior to MRI assessment (N = 672) and (ii) concurrent with MRI assessment (N=565). Four CpGs at baseline were associated with global cerebral white matter, total grey matter, and whole-brain volume (Bonferroni p≤7.41×10-8, βrange = -1.46x10-6 to 9.59 × 10-7). These CpGs were annotated to genes implicated in brain-related traits, including psychiatric disorders, development, and ageing. We did not find significant associations in the meta-analysis of the EWAS of the two sets concurrent with imaging at the corrected level.These findings reveal global brain structural changes associated with DNAm measured ~6 years previously, indicating a potential role of early DNAm modifications in brain structure. Although concurrent DNAm was not associated with global brain structure, the nominally significant findings identified here present a rationale for future investigation of associations between DNA methylation and structural brain phenotypes in larger population-based samples.

Hypothalamic–pituitary–adrenal (HPA) axis dysregulation has been commonly reported in major depressive disorder (MDD), but with considerable heterogeneity of results; potentially due to the predominant use of acute measures of an inherently variable/phasic system. Chronic longer-term measures of HPA-axis activity have yet to be systematically examined in MDD, particularly in relation to brain phenotypes, and in the context of early-life/contemporaneous stress. Here, we utilise a temporally stable measure of cumulative HPA-axis function (hair glucocorticoids) to investigate associations between cortisol, cortisone and total glucocorticoids with concurrent measures of (i) lifetime-MDD case/control status and current symptom severity, (ii) early/current-life stress and (iii) structural neuroimaging phenotypes, in N = 993 individuals from Generation Scotland (mean age = 59.1 yrs). Increased levels of hair cortisol were significantly associated with reduced global and lobar brain volumes with reductions in the frontal, temporal and cingulate regions (βrange = −0.057 to −0.104, all PFDR βrange = 0.071 to 0.115, all PFDR = β = 0.083, P = 0.017), and with reduced global and regional brain volumes (global: β = −0.059, P = 0.043; nucleus accumbens: β = −0.075, PFDR = 0.044). Associations with total glucocorticoids followed a similar pattern to the cortisol findings. In this large community-based sample, elevated glucocorticoids were significantly associated with MDD, with early, but not later-life stress, and with reduced global and regional brain phenotypes. These findings provide important foundations for future mechanistic studies to formally explore causal relationships between early adversity, chronic rather than acute measures of glucocorticoids, and neurobiological associations relevant to the aetiology of MDD.

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

The burden of large and rare copy number genetic variants (CNVs) as well as certain specific CNVs increase the risk of developing schizophrenia. Several cognitive measures are purported schizophrenia endophenotypes and may represent an intermediate point between genetics and the illness. This paper investigates the influence of CNVs on cognition. We conducted a systematic review and meta-analysis of the literature exploring the effect of CNV burden on general intelligence. We included ten primary studies with a total of 18,847 participants and found no evidence of association. In a new psychosis family study, we investigated the effects of CNVs on specific cognitive abilities. We examined the burden of large and rare CNVs (>200 kb, This work was supported by the Medical Research Council (G0901310) and the Wellcome Trust (grants 085475/B/08/Z, 085475/Z/08/Z). This study was supported by the NIHR Biomedical Research Centre at University College London Hospitals NHS Foundation Trust and University College London and by the NIHR Biomedical Research Centre for Mental Health at the South London and Maudsley NHS Foundation Trust at King’s College London. Further support to EB: Mental Health Research UK’s John Grace QC award, BMA Margaret Temple grants 2016 and 2006, MRC—Korean Health Industry Development Institute Partnering Award (MC_PC_16014), MRC New Investigator Award and a MRC Centenary Award (G0901310), National Institute of Health Research UK post-doctoral fellowship, the Psychiatry Research Trust, the Schizophrenia Research Fund, the Brain and Behaviour Research foundation’s NARSAD Young Investigator Awards 2005, 2008, Wellcome Trust Research Training Fellowship, the NIHR Biomedical Research Centre at UCLH, and the NIHR Biomedical Research Centre for Mental Health at the South London and Maudsley NHS Foundation Trust and Institute of Psychiatry King’s College London. Further support to co-authors: The Brain and Behaviour Research foundation’s (NARSAD’s) Young Investigator Award (Grant 22604, awarded to CI). The BMA Margaret Temple grant 2016 to JT. A 2014 European Research Council Marie Curie award to A Díez-Revuelta. HI has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 747429. A Medical Research Council doctoral studentship to JH-S, IA-Z and AB. A Mental Health Research UK studentship to RM. VB is supported by a Wellcome Trust Seed Award in Science (200589/Z/16/Z). FWO Senior Clinical Fellowship to RvW. The infrastructure for the GROUP consortium is funded through the Geestkracht programme of the Dutch Health Research Council (ZON-MW, grant number 10-000-1001), and matching funds from participating pharmaceutical companies (Lundbeck, AstraZeneca, Eli Lilly, Janssen Cilag) and universities and mental health care organisations (Amsterdam: Academic Psychiatric Centre of the Academic Medical Centre and the mental health institutions: GGZ Ingeest, Arkin, Dijk en Duin, GGZ Rivierduinen, Erasmus Medical Centre, GGZ Noord Holland Noord. Groningen: University Medical Centre Groningen and the mental health institutions: Lentis, GGZ Friesland, GGZ Drenthe, Dimence, Mediant, GGNet Warnsveld, Yulius Dordrecht and Parnassia psycho-medical centre The Hague. Maastricht: Maastricht University Medical Centre and the mental health institutions: GGZ Eindhoven en De Kempen, GGZ Breburg, GGZ Oost-Brabant, Vincent van Gogh voor Geestelijke Gezondheid, Mondriaan, Virenze riagg, Zuyderland GGZ, MET ggz, Universitair Centrum Sint-Jozef Kortenberg, CAPRI University of Antwerp, PC Ziekeren Sint-Truiden, PZ Sancta Maria Sint-Truiden, GGZ Overpelt, OPZ Rekem. Utrecht: University Medical Centre Utrecht and the mental health institutions Altrecht, GGZ Centraal and Delta). The Santander cohort was supported by Instituto de Salud Carlos III (PI020499, PI050427, PI060507), SENY Fundació (CI 2005-0308007), Fundacion Ramón Areces and Fundacion Marqués de Valdecilla (API07/011, API10/13). We thank Valdecilla Biobank for providing the biological PAFIP samples and associated data included in this study and for its help in the technical execution of this work; we also thank IDIVAL Neuroimaging Unit for its help in the acquisition and processing of imaging PAFIP data.

Emerging evidence suggests that obesity impacts brain physiology at multiple levels. Here we aimed to clarify the relationship between obesity and brain structure using structural MRI (n = 6420) and genetic data (n = 3907) from the ENIGMA Major Depressive Disorder (MDD) working group. Obesity (BMI > 30) was significantly associated with cortical and subcortical abnormalities in both mass-univariate and multivariate pattern recognition analyses independent of MDD diagnosis. The most pronounced effects were found for associations between obesity and lower temporo-frontal cortical thickness (maximum Cohen´s d (left fusiform gyrus) = −0.33). The observed regional distribution and effect size of cortical thickness reductions in obesity revealed considerable similarities with corresponding patterns of lower cortical thickness in previously published studies of neuropsychiatric disorders. A higher polygenic risk score for obesity significantly correlated with lower occipital surface area. In addition, a significant age-by-obesity interaction on cortical thickness emerged driven by lower thickness in older participants. Our findings suggest a neurobiological interaction between obesity and brain structure under physiological and pathological brain conditions.

There is increasing interest in using data-driven unsupervised methods to identify structural underpinnings of common mental illnesses, including Major Depressive Disorder (MDD) and associated traits such as cognition. However, studies are often limited to severe clinical cases with small sample sizes and most do not include replication. Here, we examine two relatively large samples with structural magnetic resonance imaging (MRI), measures of lifetime MDD and cognitive variables: Generation Scotland (GS subsample, N = 980), and UK Biobank (UKB, N = 8900); for discovery and replication, using an exploratory approach. Regional measures of FreeSurfer derived cortical thickness (CT), cortical surface area (CSA), cortical volume (CV) and subcortical volume (subCV) were input into a clustering process, controlling for common covariates. The main analysis steps involved constructing participant K-nearest neighbour graphs, and graph partitioning with Markov Stability to determine optimal clustering of participants. Resultant clusters were i) checked whether they were replicated in an independent cohort, and ii) tested for associations with depression status, and cognitive measures. Participants separated into two clusters based on structural brain measurements in GS subsample, with large Cohen's d effect sizes between clusters in higher order cortical regions, commonly associated with executive function and decision making. Clustering was replicated in the UKB sample, with high correlations of cluster effect sizes for CT, CSA, CV and subCV between cohorts across regions. The identified clusters were not significantly different with respect to MDD case-control status in either cohort (GS subsample: pFDR = 0.2239 - 0.6585; UKB: pFDR = 0.2003 - 0.7690). Significant differences in general cognitive ability were, however, found between the clusters for both datasets; for CSA, CV, subCV (GS subsample: d = 0.2529 - 0.3490, pFDR < 0.005; UKB: d = 0.0868 - 0.1070, pFDR < 0.005). Our results suggest that there are replicable natural groupings of participants based on cortical and subcortical brain measures, which may be related to differences in cognitive performance, but not to the MDD case-control status.