Antidepressant Exposure and DNA Methylation: Insights from a Methylome-Wide Association Study.

Importance: Understanding antidepressant mechanisms could help design more effective and tolerated treatments. Objective: Identify DNA methylation (DNAm) changes associated with antidepressant exposure. Design: Case-control methylome-wide association studies (MWAS) of antidepressant exposure were performed from blood samples collected between 2006-2011 in Generation Scotland (GS). The summary statistics were tested for enrichment in specific tissues, gene ontologies and an independent MWAS in the Netherlands Study of Depression and Anxiety (NESDA). A methylation profile score (MPS) was derived and tested for its association with antidepressant exposure in eight independent cohorts, alongside prospective data from GS. Setting: Cohorts; GS, NESDA, FTC, SHIP-Trend, FOR2107, LBC1936, MARS-UniDep, ALSPAC, E-Risk, and NTR. Participants: Participants with DNAm data and self-report/prescription derived antidepressant exposure. Main Outcomes and Measures: Whole-blood DNAm levels were assayed by the EPIC/450K Illumina array (9 studies, N exposed = 661, N unexposed = 9,575) alongside MBD-Seq in NESDA (N exposed = 398, N unexposed = 414). Antidepressant exposure was measured by self- report and/or antidepressant prescriptions. Results: The self-report MWAS (N = 16,536, N exposed = 1,508, mean age = 48, 59% female) and the prescription-derived MWAS (N = 7,951, N exposed = 861, mean age = 47, 59% female), found hypermethylation at seven and four DNAm sites (p < 9.42x10 -8 ), respectively. The top locus was cg26277237 ( KANK1, p self-report = 9.3x10 -13 , p prescription = 6.1x10 -3 ). The self-report MWAS found a differentially methylated region, mapping to DGUOK-AS1 ( p adj = 5.0x10 -3 ) alongside significant enrichment for genes expressed in the amygdala, the "synaptic vesicle membrane" gene ontology and the top 1% of CpGs from the NESDA MWAS (OR = 1.39, p < 0.042). The MPS was associated with antidepressant exposure in meta-analysed data from external cohorts (N studies = 9, N = 10,236, N exposed = 661, f3 = 0.196, p < 1x10 -4 ). Conclusions and Relevance: Antidepressant exposure is associated with changes in DNAm across different cohorts. Further investigation into these changes could inform on new targets for antidepressant treatments. 3 Key Points: Question: Is antidepressant exposure associated with differential whole blood DNA methylation?Findings: In this methylome-wide association study of 16,536 adults across Scotland, antidepressant exposure was significantly associated with hypermethylation at CpGs mapping to KANK1 and DGUOK-AS1. A methylation profile score trained on this sample was significantly associated with antidepressant exposure (pooled f3 [95%CI]=0.196 [0.105, 0.288], p < 1x10 -4 ) in a meta-analysis of external datasets. Meaning: Antidepressant exposure is associated with hypermethylation at KANK1 and DGUOK-AS1 , which have roles in mitochondrial metabolism and neurite outgrowth. If replicated in future studies, targeting these genes could inform the design of more effective and better tolerated treatments for depression.

A combined analysis of genetically correlated traits identifies 187 loci and a role for neurogenesis and myelination in intelligence.

Intelligence, or general cognitive function, is phenotypically and genetically correlated with many traits, including a wide range of physical, and mental health variables. Education is strongly genetically correlated with intelligence (rg = 0.70). We used these findings as foundations for our use of a novel approach-multi-trait analysis of genome-wide association studies (MTAG; Turley et al. 2017)-to combine two large genome-wide association studies (GWASs) of education and intelligence, increasing statistical power and resulting in the largest GWAS of intelligence yet reported. Our study had four goals: first, to facilitate the discovery of new genetic loci associated with intelligence; second, to add to our understanding of the biology of intelligence differences; third, to examine whether combining genetically correlated traits in this way produces results consistent with the primary phenotype of intelligence; and, finally, to test how well this new meta-analytic data sample on intelligence predicts phenotypic intelligence in an independent sample. By combining datasets using MTAG, our functional sample size increased from 199,242 participants to 248,482. We found 187 independent loci associated with intelligence, implicating 538 genes, using both SNP-based and gene-based GWAS. We found evidence that neurogenesis and myelination-as well as genes expressed in the synapse, and those involved in the regulation of the nervous system-may explain some of the biological differences in intelligence. The results of our combined analysis demonstrated the same pattern of genetic correlations as those from previous GWASs of intelligence, providing support for the meta-analysis of these genetically-related phenotypes.

Brain age predicts mortality.

Age-associated disease and disability are placing a growing burden on society. However, ageing does not affect people uniformly. Hence, markers of the underlying biological ageing process are needed to help identify people at increased risk of age-associated physical and cognitive impairments and ultimately, death. Here, we present such a biomarker, 'brain-predicted age', derived using structural neuroimaging. Brain-predicted age was calculated using machine-learning analysis, trained on neuroimaging data from a large healthy reference sample (N=2001), then tested in the Lothian Birth Cohort 1936 (N=669), to determine relationships with age-associated functional measures and mortality. Having a brain-predicted age indicative of an older-appearing brain was associated with: weaker grip strength, poorer lung function, slower walking speed, lower fluid intelligence, higher allostatic load and increased mortality risk. Furthermore, while combining brain-predicted age with grey matter and cerebrospinal fluid volumes (themselves strong predictors) not did improve mortality risk prediction, the combination of brain-predicted age and DNA-methylation-predicted age did. This indicates that neuroimaging and epigenetics measures of ageing can provide complementary data regarding health outcomes. Our study introduces a clinically-relevant neuroimaging ageing biomarker and demonstrates that combining distinct measurements of biological ageing further helps to determine risk of age-related deterioration and death.

A DNA methylation biomarker of alcohol consumption.

The lack of reliable measures of alcohol intake is a major obstacle to the diagnosis and treatment of alcohol-related diseases. Epigenetic modifications such as DNA methylation may provide novel biomarkers of alcohol use. To examine this possibility, we performed an epigenome-wide association study of methylation of cytosine-phosphate-guanine dinucleotide (CpG) sites in relation to alcohol intake in 13 population-based cohorts (ntotal=13 317; 54% women; mean age across cohorts 42-76 years) using whole blood (9643 European and 2423 African ancestries) or monocyte-derived DNA (588 European, 263 African and 400 Hispanic ancestry) samples. We performed meta-analysis and variable selection in whole-blood samples of people of European ancestry (n=6926) and identified 144 CpGs that provided substantial discrimination (area under the curve=0.90-0.99) for current heavy alcohol intake (⩾42 g per day in men and ⩾28 g per day in women) in four replication cohorts. The ancestry-stratified meta-analysis in whole blood identified 328 (9643 European ancestry samples) and 165 (2423 African ancestry samples) alcohol-related CpGs at Bonferroni-adjusted P<1 × 10-7. Analysis of the monocyte-derived DNA (n=1251) identified 62 alcohol-related CpGs at P<1 × 10-7. In whole-blood samples of people of European ancestry, we detected differential methylation in two neurotransmitter receptor genes, the γ-Aminobutyric acid-A receptor delta and γ-aminobutyric acid B receptor subunit 1; their differential methylation was associated with expression levels of a number of genes involved in immune function. In conclusion, we have identified a robust alcohol-related DNA methylation signature and shown the potential utility of DNA methylation as a clinically useful diagnostic test to detect current heavy alcohol consumption.

An epigenome-wide association study meta-analysis of educational attainment.

The epigenome is associated with biological factors, such as disease status, and environmental factors, such as smoking, alcohol consumption and body mass index. Although there is a widespread perception that environmental influences on the epigenome are pervasive and profound, there has been little evidence to date in humans with respect to environmental factors that are biologically distal. Here we provide evidence on the associations between epigenetic modifications-in our case, CpG methylation-and educational attainment (EA), a biologically distal environmental factor that is arguably among the most important life-shaping experiences for individuals. Specifically, we report the results of an epigenome-wide association study meta-analysis of EA based on data from 27 cohort studies with a total of 10 767 individuals. We find nine CpG probes significantly associated with EA. However, robustness analyses show that all nine probes have previously been found to be associated with smoking. Only two associations remain when we perform a sensitivity analysis in the subset of never-smokers, and these two probes are known to be strongly associated with maternal smoking during pregnancy, and thus their association with EA could be due to correlation between EA and maternal smoking. Moreover, the effect sizes of the associations with EA are far smaller than the known associations with the biologically proximal environmental factors alcohol consumption, body mass index, smoking and maternal smoking during pregnancy. Follow-up analyses that combine the effects of many probes also point to small methylation associations with EA that are highly correlated with the combined effects of smoking. If our findings regarding EA can be generalized to other biologically distal environmental factors, then they cast doubt on the hypothesis that such factors have large effects on the epigenome.

Short telomere length is associated with impaired cognitive performance in European ancestry cohorts.

The association between telomere length (TL) dynamics on cognitive performance over the life-course is not well understood. This study meta-analyses observational and causal associations between TL and six cognitive traits, with stratifications on APOE genotype, in a Mendelian Randomization (MR) framework. Twelve European cohorts (N=17 052; mean age=59.2±8.8 years) provided results for associations between qPCR-measured TL (T/S-ratio scale) and general cognitive function, mini-mental state exam (MMSE), processing speed by digit symbol substitution test (DSST), visuospatial functioning, memory and executive functioning (STROOP). In addition, a genetic risk score (GRS) for TL including seven known genetic variants for TL was calculated, and used in associations with cognitive traits as outcomes in all cohorts. Observational analyses showed that longer telomeres were associated with better scores on DSST (ß=0.051 per s.d.-increase of TL; 95% confidence interval (CI): 0.024, 0.077; P=0.0002), and MMSE (ß=0.025; 95% CI: 0.002, 0.047; P=0.03), and faster STROOP (ß=-0.053; 95% CI: -0.087, -0.018; P=0.003). Effects for DSST were stronger in APOE ɛ4 non-carriers (ß=0.081; 95% CI: 0.045, 0.117; P=1.0 × 10-5), whereas carriers performed better in STROOP (ß=-0.074; 95% CI: -0.140, -0.009; P=0.03). Causal associations were found for STROOP only (ß=-0.598 per s.d.-increase of TL; 95% CI: -1.125, -0.072; P=0.026), with a larger effect in ɛ4-carriers (ß=-0.699; 95% CI: -1.330, -0.069; P=0.03). Two-sample replication analyses using CHARGE summary statistics showed causal effects between TL and general cognitive function and DSST, but not with STROOP. In conclusion, we suggest causal effects from longer TL on better cognitive performance, where APOE ɛ4-carriers might be at differential risk.

Genome-wide association study of cognitive functions and educational attainment in UK Biobank (N=112 151).

People's differences in cognitive functions are partly heritable and are associated with important life outcomes. Previous genome-wide association (GWA) studies of cognitive functions have found evidence for polygenic effects yet, to date, there are few replicated genetic associations. Here we use data from the UK Biobank sample to investigate the genetic contributions to variation in tests of three cognitive functions and in educational attainment. GWA analyses were performed for verbal-numerical reasoning (N=36 035), memory (N=112 067), reaction time (N=111 483) and for the attainment of a college or a university degree (N=111 114). We report genome-wide significant single-nucleotide polymorphism (SNP)-based associations in 20 genomic regions, and significant gene-based findings in 46 regions. These include findings in the ATXN2, CYP2DG, APBA1 and CADM2 genes. We report replication of these hits in published GWA studies of cognitive function, educational attainment and childhood intelligence. There is also replication, in UK Biobank, of SNP hits reported previously in GWA studies of educational attainment and cognitive function. GCTA-GREML analyses, using common SNPs (minor allele frequency>0.01), indicated significant SNP-based heritabilities of 31% (s.e.m.=1.8%) for verbal-numerical reasoning, 5% (s.e.m.=0.6%) for memory, 11% (s.e.m.=0.6%) for reaction time and 21% (s.e.m.=0.6%) for educational attainment. Polygenic score analyses indicate that up to 5% of the variance in cognitive test scores can be predicted in an independent cohort. The genomic regions identified include several novel loci, some of which have been associated with intracranial volume, neurodegeneration, Alzheimer's disease and schizophrenia.

Assessing the genetic overlap between BMI and cognitive function.

Obesity and low cognitive function are associated with multiple adverse health outcomes across the life course. They have a small phenotypic correlation (r=-0.11; high body mass index (BMI)-low cognitive function), but whether they have a shared genetic aetiology is unknown. We investigated the phenotypic and genetic correlations between the traits using data from 6815 unrelated, genotyped members of Generation Scotland, an ethnically homogeneous cohort from five sites across Scotland. Genetic correlations were estimated using the following: same-sample bivariate genome-wide complex trait analysis (GCTA)-GREML; independent samples bivariate GCTA-GREML using Generation Scotland for cognitive data and four other samples (n=20 806) for BMI; and bivariate LDSC analysis using the largest genome-wide association study (GWAS) summary data on cognitive function (n=48 462) and BMI (n=339 224) to date. The GWAS summary data were also used to create polygenic scores for the two traits, with within- and cross-trait prediction taking place in the independent Generation Scotland cohort. A large genetic correlation of -0.51 (s.e. 0.15) was observed using the same-sample GCTA-GREML approach compared with -0.10 (s.e. 0.08) from the independent-samples GCTA-GREML approach and -0.22 (s.e. 0.03) from the bivariate LDSC analysis. A genetic profile score using cognition-specific genetic variants accounts for 0.08% (P=0.020) of the variance in BMI and a genetic profile score using BMI-specific variants accounts for 0.42% (P=1.9 × 10(-7)) of the variance in cognitive function. Seven common genetic variants are significantly associated with both traits at P<5 × 10(-5), which is significantly more than expected by chance (P=0.007). All these results suggest there are shared genetic contributions to BMI and cognitive function.

Shared genetic aetiology between cognitive functions and physical and mental health in UK Biobank (N=112 151) and 24 GWAS consortia.

Causes of the well-documented association between low levels of cognitive functioning and many adverse neuropsychiatric outcomes, poorer physical health and earlier death remain unknown. We used linkage disequilibrium regression and polygenic profile scoring to test for shared genetic aetiology between cognitive functions and neuropsychiatric disorders and physical health. Using information provided by many published genome-wide association study consortia, we created polygenic profile scores for 24 vascular-metabolic, neuropsychiatric, physiological-anthropometric and cognitive traits in the participants of UK Biobank, a very large population-based sample (N=112 151). Pleiotropy between cognitive and health traits was quantified by deriving genetic correlations using summary genome-wide association study statistics and to the method of linkage disequilibrium score regression. Substantial and significant genetic correlations were observed between cognitive test scores in the UK Biobank sample and many of the mental and physical health-related traits and disorders assessed here. In addition, highly significant associations were observed between the cognitive test scores in the UK Biobank sample and many polygenic profile scores, including coronary artery disease, stroke, Alzheimer's disease, schizophrenia, autism, major depressive disorder, body mass index, intracranial volume, infant head circumference and childhood cognitive ability. Where disease diagnosis was available for UK Biobank participants, we were able to show that these results were not confounded by those who had the relevant disease. These findings indicate that a substantial level of pleiotropy exists between cognitive abilities and many human mental and physical health disorders and traits and that it can be used to predict phenotypic variance across samples.

Genetic contributions to variation in general cognitive function: a meta-analysis of genome-wide association studies in the CHARGE consortium (N=53949).

General cognitive function is substantially heritable across the human life course from adolescence to old age. We investigated the genetic contribution to variation in this important, health- and well-being-related trait in middle-aged and older adults. We conducted a meta-analysis of genome-wide association studies of 31 cohorts (N=53,949) in which the participants had undertaken multiple, diverse cognitive tests. A general cognitive function phenotype was tested for, and created in each cohort by principal component analysis. We report 13 genome-wide significant single-nucleotide polymorphism (SNP) associations in three genomic regions, 6q16.1, 14q12 and 19q13.32 (best SNP and closest gene, respectively: rs10457441, P=3.93 × 10(-9), MIR2113; rs17522122, P=2.55 × 10(-8), AKAP6; rs10119, P=5.67 × 10(-9), APOE/TOMM40). We report one gene-based significant association with the HMGN1 gene located on chromosome 21 (P=1 × 10(-6)). These genes have previously been associated with neuropsychiatric phenotypes. Meta-analysis results are consistent with a polygenic model of inheritance. To estimate SNP-based heritability, the genome-wide complex trait analysis procedure was applied to two large cohorts, the Atherosclerosis Risk in Communities Study (N=6617) and the Health and Retirement Study (N=5976). The proportion of phenotypic variation accounted for by all genotyped common SNPs was 29% (s.e.=5%) and 28% (s.e.=7%), respectively. Using polygenic prediction analysis, ~1.2% of the variance in general cognitive function was predicted in the Generation Scotland cohort (N=5487; P=1.5 × 10(-17)). In hypothesis-driven tests, there was significant association between general cognitive function and four genes previously associated with Alzheimer's disease: TOMM40, APOE, ABCG1 and MEF2C.

Human cognitive ability is influenced by genetic variation in components of postsynaptic signalling complexes assembled by NMDA receptors and MAGUK proteins.

Differences in general cognitive ability (intelligence) account for approximately half of the variation in any large battery of cognitive tests and are predictive of important life events including health. Genome-wide analyses of common single-nucleotide polymorphisms indicate that they jointly tag between a quarter and a half of the variance in intelligence. However, no single polymorphism has been reliably associated with variation in intelligence. It remains possible that these many small effects might be aggregated in networks of functionally linked genes. Here, we tested a network of 1461 genes in the postsynaptic density and associated complexes for an enriched association with intelligence. These were ascertained in 3511 individuals (the Cognitive Ageing Genetics in England and Scotland (CAGES) consortium) phenotyped for general cognitive ability, fluid cognitive ability, crystallised cognitive ability, memory and speed of processing. By analysing the results of a genome wide association study (GWAS) using Gene Set Enrichment Analysis, a significant enrichment was found for fluid cognitive ability for the proteins found in the complexes of N-methyl-D-aspartate receptor complex; P=0.002. Replication was sought in two additional cohorts (N=670 and 2062). A meta-analytic P-value of 0.003 was found when these were combined with the CAGES consortium. The results suggest that genetic variation in the macromolecular machines formed by membrane-associated guanylate kinase (MAGUK) scaffold proteins and their interaction partners contributes to variation in intelligence.

Associations between polymorphisms in five inflammation-related genes and cognitive ability in older persons.

Several studies have found associations between inflammatory biomarker levels and cognitive ability. This study tested the relationship between polymorphisms in genes that are associated with or encode the biomarkers and cognitive ability and estimated lifetime cognitive change. Data came from the aspirin for asymptomatic atherosclerosis trial (n = 2091, mean age = 67.2 years ). Twelve single nucleotide polymorphisms (SNPs) were genotyped from five genes (IL-1alpha, IL-1beta, IL-6, HNF-1A and F13A1). Cognition was assessed via administration of a five-test battery of psychometric tests, which were used to derive a general intelligence factor, g. A vocabulary-based cognitive test was also administered and adjusted for in the analysis to enable an estimation of lifetime cognitive change. Age- and sex-adjusted analyses yielded one weakly significant association between the IL-1alpha rs2856838 SNP and a measure of mental flexibility/processing speed (P = 0.044). Adjustment for the vocabulary-based scores resulted in a single, significant association between the IL-1alpha rs3783546 SNP and a measure of processing speed (P = 0.048). There is little evidence to suggest an association between SNPs in the inflammation-related genes IL-1alpha, IL-1beta, IL-6, TCF-1 and F13A1 and cognition in an elderly population of community-dwelling Scottish citizens.