AbstractAlzheimer’s disease (AD) is a chronic neurodegenerative disorder characterized by the progressive accumulation of amyloid-β (Aβ) plaques and neurofibrillary tangles in the neocortex. Recent studies have implicated a role for regulatory genomic variation in AD progression, finding widespread evidence for altered DNA methylation associated with neuropathology. To date, however, no study has systematically examined other types of regulatory genomic modifications in AD. In this study, we quantified genome-wide patterns of lysine H3K27 acetylation (H3K27ac) - a robust mark of active enhancers and promoters that is strongly correlated with gene expression and transcription factor binding - in entorhinal cortex samples from AD cases and matched controls (n = 47) using chromatin immunoprecipitation followed by highly parallel sequencing (ChIP-seq). Across ~182,000 robustly detected H3K27ac peak regions, we found widespread acetylomic variation associated with AD neuropathology, identifying 4,162 differential peaks (FDR < 0.05) between AD cases and controls. These differentially acetylated peaks are enriched in disease-specific biological pathways and include regions annotated to multiple genes directly involved in the progression of Aβ and tau pathology (e.g. APP, PSEN1, PSEN2, MAPT), as well as genomic regions containing variants associated with sporadic late-onset AD. This is the first study of variable H3K27ac yet undertaken in AD and the largest study investigating this modification in the entorhinal cortex. In addition to identifying molecular pathways associated with AD neuropathology, we present a framework for genome-wide studies of histone modifications in complex disease, integrating our data with results obtained from genome-wide association studies as well as other epigenetic marks profiled on the same samples.

ABSTRACTEpigenome-wide association studies of Alzheimer’s disease have highlighted neuropathology-associated DNA methylation differences, although existing studies have been limited in sample size and utilized different brain regions. Here, we combine data from six DNA methylomic studies of Alzheimer’s disease (N=1,453 unique individuals) to identify differential methylation associated with Braak stage in different brain regions and across cortex. We identify 236 CpGs in the prefrontal cortex, 95 CpGs in the temporal gyrus and ten CpGs in the entorhinal cortex at Bonferroni significance, with none in the cerebellum. Our cross-cortex meta-analysis (N=1,408 donors) identifies 220 CpGs associated with neuropathology, annotated to 121 genes, of which 84 genes have not been previously reported at this significance threshold. We have replicated our findings using two further DNA methylomic datasets consisting of a further > 600 unique donors. The meta-analysis summary statistics are available in our online data resource (www.epigenomicslab.com/ad-meta-analysis/).

Summary/AbstractBackgroundMost epigenome-wide association studies (EWAS) quantify DNA methylation (DNAm) in peripheral tissues such as whole blood to identify positions in the genome where variation is statistically associated with a trait or exposure. As whole blood comprises a mix of cell types, it is unclear whether trait-associated variation is specific to an individual cellular population.MethodsWe collected three peripheral tissues (whole blood, buccal and nasal epithelial cells) from thirty individuals. Whole blood samples were subsequently processed using fluorescence-activated cell sorting (FACS) to purify five constituent cell-types (monocytes, granulocytes, CD4+ T cells, CD8+ T cells, and B cells). DNAm was profiled in all eight sample-types from each individual using the Illumina EPIC array.ResultsWe identified significant differences in both the level and variability of DNAm between different tissues and cell types, and DNAm data-derived estimates of age and smoking were found to differ dramatically across sample types from the same individual. We found that for the majority of loci variation in DNAm in individual blood cell types was only weakly predictive of variance in DNAm measured in whole blood, however, the proportion of variance explained was greater than that explained by either buccal or nasal tissues. Instead we observe that DNAm variation in whole blood is additively influenced by a combination of the major blood cell types. For a subset of sites variable DNAm detected in whole blood can be attributed to variation in a single blood cell type providing potential mechanistic insight.ConclusionsWe identified major differences in DNAm between blood cell types and peripheral tissues, with each sample type being characterized by a unique DNAm signature across multiple loci. Our results suggest that associations between whole blood DNAm and traits or exposures reflect differences in multiple cell types and provide important insights for the interpretation of EWAS performed in whole blood.Key MessagesWe identified major differences in DNA methylation between blood cell types and peripheral tissues, with each sample type being characterized by a unique DNA methylation signature across multiple loci.Estimates of DNAmAge and tobacco smoking from DNA methylation data can be highly variable across different sample types collected from the same individual at the same time.While individual blood cell types did predict more of the variation in whole blood compared to buccal epithelial and nasal epithelial cells, the percentage of variance explained was still small.Instead our data indicate that at the majority of sites, variation in multiple blood cell types additively combines to drive variation in DNA methylation in whole blood.There are subset of sites where variable DNA methylation detected in whole blood can be attributed to variation in a single blood cell type.

AbstractBackgroundThe combination of sodium bisulfite treatment with highly-parallel sequencing is a common method for quantifying DNA methylation across the genome. The power to detect between-group differences in DNA methylation using bisulfite-sequencing approaches is influenced by both experimental (e.g. read depth, missing data and sample size) and biological (e.g. mean level of DNA methylation and difference between groups) parameters. There is, however, no consensus about the optimal thresholds for filtering bisulfite sequencing data with implications for the reproducibility of findings in epigenetic epidemiology.ResultsWe used a large reduced representation bisulfite sequencing (RRBS) dataset to assess the distribution of read depth across DNA methylation sites and the extent of missing data. To investigate how various study variables influence power to identify DNA methylation differences between groups, we developed a framework for simulating bisulfite sequencing data. As expected, sequencing read depth, group size, and the magnitude of DNA methylation difference between groups all impacted upon statistical power. The influence on power was not dependent on one specific parameter, but reflected the combination of study-specific variables. As a resource to the community, we have developed a tool, POWEREDBiSeq, which utilizes our simulation framework to predict study-specific power for the identification of DNAm differences between groups, taking into account user-defined read depth filtering parameters and the minimum sample size per group.ConclusionsOur data-driven approach highlights the importance of filtering bisulfite-sequencing data by minimum read depth and illustrates how the choice of threshold is influenced by the specific study design and the expected differences between groups being compared. The POWEREDBiSeq tool can help users identify the level of data filtering needed to optimize power and aims to improve the reproducibility of bisulfite sequencing studies.

AbstractDNA 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 2,425 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 peripherl 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.

AbstractSex is an important covariate of epigenome-wide association studies due to its strong influence on DNA methylation patterns across numerous genomic positions. Nevertheless, many samples on the Gene Expression Omnibus (GEO) frequently lack a sex annotation or are incorrectly labelled. Considering the influence that sex imposes on DNA methylation patterns, it is necessary to ensure that methods for filtering poor samples and checking of sex assignment are accurate and widely applicable. In this paper, we presented a novel method to predict sex using only DNA methylation density signals, which can be readily applied to almost all DNA methylation datasets of different formats (raw IDATs or text files with only density signals) uploaded to GEO. We identified 4345 significantly (p < 0.01) sex-associated CpG sites present on both 450K and EPIC arrays, and constructed a sex classifier based on the two first components of PCAs from the two sex chromosomes. The proposed method is constructed using whole blood samples and exhibits good performance across a wide range of tissues. We further demonstrated that our method can be used to identify samples with sex chromosome aneuploidy, this function is validated by five Turner syndrome cases and one Klinefelter syndrome case. The proposed method has been integrated into the wateRmelon Bioconductor package.

ABSTRACTMeasures to quantify changes in the pace of biological aging in response to intervention are needed to evaluate geroprotective interventions for humans. We used data from the Dunedin Study 1972-3 birth cohort tracking within-individual decline in 19 indicators of organ-system integrity across four time points spanning two decades to model Pace of Aging. We distilled this two-decade Pace of Aging into a single-time-point DNA-methylation blood-test using elastic-net regression and a DNA-methylation dataset restricted to exclude probes with low test-retest reliability. We evaluated the resulting measure, named DunedinPACE, in five additional datasets. DunedinPACE showed high test-retest reliability, was associated with morbidity, disability, and mortality, and indicated faster aging in young adults with childhood adversity. DunedinPACE effect-sizes were similar to GrimAge Clock effect-sizes. In analysis of morbidity, disability, and mortality, DunedinPACE and added incremental prediction beyond GrimAge. DunedinPACE is a novel blood biomarker of the pace of aging for gerontology and geroscience.

AbstractBackgroundThe study of biological age acceleration may help identify at-risk individuals and contribute to reduce the rising global burden of age-related diseases. Using DNA methylation (DNAm) clocks, we investigated biological aging in schizophrenia (SCZ), a severe mental illness that is associated with an increased prevalence of age-related disabilities and morbidities. In a multi-cohort whole blood sample consisting of 1,090 SCZ cases and 1,206 controls, we investigated differential aging using three DNAm clocks (i.e. Hannum, Horvath, Levine). These clocks are highly predictive of chronological age and are known to capture different processes of biological aging.ResultsWe found that blood-based DNAm aging is significantly altered in SCZ with age- and sex-specific effects that differ between clocks and map to distinct chronological age windows. Most notably, differential phenotypic age (Levine clock) was most pronounced in female SCZ patients in later adulthood compared to matched controls. Female patients with high SCZ polygenic risk scores (PRS) present the highest age acceleration in this age group with +4.30 years (CI: 2.40-6.20, P=1.3E-05). Phenotypic age and SCZ PRS contribute additively to the illness and together explain up to 22.4% of the variance in disease status in this study. This suggests that combining genetic and epigenetic predictors may improve predictions of disease outcomes.ConclusionsSince increased phenotypic age is associated with increased risk of all-cause mortality, our findings indicate that specific and identifiable patient groups are at increased mortality risk as measured by the Levine clock. These results provide new biological insights into the aging landscape of SCZ with age- and sex-specific effects and warrant further investigations into the potential of DNAm clocks as clinical biomarkers that may help with disease management in schizophrenia.

AbstractParkinson’s disease (PD) and dementia with Lewy bodies (DLB) are closely related progressive disorders with no available causal therapy, neuropathologically characterized by intraneuronal aggregates of misfolded α-synuclein. To explore the role of DNA methylation changes in PD and DLB pathogenesis, we performed an epigenome-wide association study (EWAS) of 322 postmortem frontal cortex samples and replicated results in an independent set of 219 donors. We report novel differentially methylated replicating loci associated with Braak Lewy body stage near SFMBT2, PHYHIP, BRF1/PACS2 and DGKG. The DGKG locus also showed evidence of DNA methylation changes in the earliest, preclinical stage of disease. Differentially methylated probes were independent of known PD genetic risk alleles. Meta-analysis provided suggestive evidence for a differentially methylated locus within the chromosomal region affected by the PD-associated 22q11.2 deletion. Our findings elucidate novel disease pathways in PD and DLB and generate hypotheses for future molecular studies of Lewy body pathology.

AbstractAlternative splicing is a post-transcriptional regulatory mechanism producing multiple distinct mRNA molecules from a single pre-mRNA. Alternative splicing has a prominent role in the central nervous system, impacting neurodevelopment and various neuronal functions as well as being increasingly implicated in brain disorders including autism, schizophrenia and Alzheimer’s disease. Standard short-read RNA-Seq approaches only sequence fragments of the mRNA molecule, making it difficult to accurately characterize the true nature of RNA isoform diversity. In this study, we used long-read isoform sequencing (Iso-Seq) to generate full-length cDNA sequences and map transcript diversity in the human and mouse cerebral cortex. We identify widespread RNA isoform diversity amongst expressed genes in the cortex, including many novel transcripts not present in existing genome annotations. Alternative splicing events were found to make a major contribution to RNA isoform diversity in the cortex, with intron retention being a relatively common event associated with nonsense-mediated decay and reduced transcript expression. of note, we found evidence for transcription from novel (unannotated genes) and fusion events between neighbouring genes. Although global patterns of RNA isoform diversity were found to be generally similar between human and mouse cortex, we identified some notable exceptions. We also identified striking developmental changes in transcript diversity, with differential transcript usage between human adult and fetal cerebral cortex. Finally, we found evidence for extensive isoform diversity in genes associated with autism, schizophrenia and Alzheimer’s disease. Our data confirm the importance of alternative splicing in the cerebral cortex, dramatically increasing transcriptional diversity and representing an important mechanism underpinning gene regulation in the brain. We provide this transcript level data as a resource to the scientific community.

AbstractAlzheimer’s disease is a highly heritable, common neurodegenerative disease characterised neuropathologically by the accumulation of β-amyloid plaques and tau-containing neurofibrillary tangles. In addition to the well-established risk associated with the APOE locus, there has been considerable success in identifying additional genetic variants associated with Alzheimer’s disease. Major challenges in understanding how genetic risk influences the development of Alzheimer’s disease are clinical and neuropathological heterogeneity, and the high level of accompanying comorbidities. We report a multimodal analysis integrating longitudinal clinical and cognitive assessment with neuropathological data collected as part of the Brains for Dementia Research (BDR) study to understand how genetic risk factors for Alzheimer’s disease influence the development of neuropathology and clinical performance. 693 donors in the BDR cohort with genetic data, semi-quantitative neuropathology measurements, cognitive assessments and established diagnostic criteria were included in this study. We tested the association of APOE genotype and Alzheimer’s disease polygenic risk score - a quantitative measure of genetic burden - with survival, four common neuropathological features in Alzheimer’s disease brains (neurofibrillary tangles, β-amyloid plaques, Lewy bodies and TDP-43 proteinopathy), clinical status (clinical dementia rating) and cognitive performance (Mini-Mental State Exam, Montreal Cognitive Assessment). The APOE ε4 allele was significantly associated with younger age of death in the BDR cohort. Our analyses of neuropathology highlighted two independent pathways from APOE ε4, one where β-amyloid accumulation mediates the development of tauopathy, and a second characterized by direct effects on tauopathy independent of β-amyloidosis. Although we also detected association between APOE ε4 and dementia status and cognitive performance, these were all mediated by tauopathy, highlighting that they are a consequence of the neuropathological changes. Analyses of polygenic risk score identified associations with tauopathy and β-amyloidosis, which appeared to have both shared and unique contributions, suggesting that different genetic variants associated with Alzheimer’s disease affect different features of neuropathology to different degrees. Taken together, our results provide insight into how genetic risk for Alzheimer’s disease influences both the clinical and pathological features of dementia, increasing our understanding about the interplay between APOE genotype and other genetic risk factors.

AbstractAmyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with an estimated heritability of around 50%. DNA methylation patterns can serve as biomarkers of (past) exposures and disease progression, as well as providing a potential mechanism that mediates genetic or environmental risk. Here, we present a blood-based epigenome-wide association study (EWAS) meta-analysis in 10,462 samples (7,344 ALS patients and 3,118 controls), representing the largest case-control study of DNA methylation for any disease to date. We identified a total of 45 differentially methylated positions (DMPs) annotated to 42 genes, which are enriched for pathways and traits related to metabolism, cholesterol biosynthesis, and immunity. We show that DNA-methylation-based proxies for HDL-cholesterol, BMI, white blood cell (WBC) proportions and alcohol intake were independently associated with ALS. Integration of these results with our latest GWAS showed that cholesterol biosynthesis was causally related to ALS. Finally, we found that DNA methylation levels at several DMPs and blood cell proportion estimates derived from DNA methylation data, are associated with survival rate in patients, and could represent indicators of underlying disease processes.

AbstractCharacterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. Here we describe results of DNA methylation-quantitative trait loci (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTL of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We reveal that the genetic architecture of DNAm levels is highly polygenic and DNAm exhibits signatures of negative and positive natural selection. Using shared genetic control between distal DNAm sites we construct networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic factors are associated with both blood DNAm levels and complex diseases but in most cases these associations do not reflect causal relationships from DNAm to trait or vice versa indicating a more complex genotype-phenotype map than has previously been hypothesised.

ABSTRACTObjectivePsychosis - a complex and heterogeneous neuropsychiatric condition characterized by hallucinations and delusions - is a common feature of schizophrenia. There is evidence for altered DNA methylation (DNAm) associated with schizophrenia in both brain and peripheral tissues. We aimed to undertake a systematic analysis of variable DNAm associated with psychosis, schizophrenia, and treatment-resistant schizophrenia, also exploring measures of biological ageing, smoking, and blood cell composition derived from DNAm data to identify molecular biomarkers of disease.MethodsWe quantified DNAm across the genome in blood samples from 4,483 participants from seven case-control cohorts including patients with schizophrenia or first-episode psychosis. Measures of biological age, cellular composition and smoking status were derived from DNAm data using established algorithms. DNAm and derived measures were analyzed within each cohort and the results combined by meta-analysis.ResultsPsychosis cases were characterized by significant differences in measures of blood cell proportions and elevated smoking exposure derived from the DNAm data, with the largest differences seen in treatment-resistant schizophrenia patients. DNAm at 95 CpG sites was significantly different between psychosis cases and controls, with 1,048 differentially methylated positions (DMPs) identified between schizophrenia cases and controls. Schizophrenia-associated DMPs colocalize to regions identified in genetic association studies, with genes annotated to these sites enriched for pathways relevant to disease. Finally, a number of the schizophrenia associated differences were only present in the treatment-resistant schizophrenia subgroup.ConclusionsWe show that DNAm data can be leveraged to derive measures of blood cell counts and smoking that are strongly associated with psychosis. Our DNAm meta-analysis identified multiple DMPs associated with both psychosis and a more refined diagnosis of schizophrenia, with evidence for differential methylation associated with treatment-resistant schizophrenia that potentially reflects exposure to clozapine.

ABSTRACTCharacterizing the complex relationship between genetic, epigenetic and transcriptomic variation has the potential to increase understanding about the mechanisms underpinning health and disease phenotypes. In this study, we describe the most comprehensive analysis of common genetic variation on DNA methylation (DNAm) to date, using the Illumina EPIC array to profile samples from the UK Household Longitudinal study. We identified 12,689,548 significant DNA methylation quantitative trait loci (mQTL) associations (P < 6.52x10-14) occurring between 2,907,234 genetic variants and 93,268 DNAm sites, including a large number not identified using previous DNAm-profiling methods. We demonstrate the utility of these data for interpreting the functional consequences of common genetic variation associated with > 60 human traits, using Summary data–based Mendelian Randomization (SMR) to identify 1,662 pleiotropic associations between 36 complex traits and 1,246 DNAm sites. We also use SMR to characterize the relationship between DNAm and gene expression, identifying 6,798 pleiotropic associations between 5,420 DNAm sites and the transcription of 1,702 genes. Our mQTL database and SMR results are available via a searchable online database (http://www.epigenomicslab.com/online-data-resources/) as a resource to the research community.

AbstractInduced pluripotent stem cells (iPSCs) and their differentiated neurons (iPSC-neurons) are a widely used cellular model in the research of the central nervous system. However, it is unknown how well they capture age-associated processes, particularly given that pluripotent cells are only present during the earliest stages of mammalian development. Epigenetic clocks utilize coordinated age-associated changes in DNA methylation to make predictions that correlate strongly with chronological age. It has been shown that the induction of pluripotency rejuvenates predicted epigenetic age. As existing clocks are not optimized for the study of brain development, we developed the fetal brain clock (FBC), a bespoke epigenetic clock trained in human prenatal brain samples in order to investigate more precisely the epigenetic age of iPSCs and iPSC-neurons. The FBC was tested in two independent validation cohorts across a total of 194 samples, confirming that the FBC outperforms other established epigenetic clocks in fetal brain cohorts. We applied the FBC to DNA methylation data from iPSCs and iPSC-derived neuronal precursor cells and neurons, finding that these cell types are epigenetically characterized as having an early fetal age. Furthermore, while differentiation from iPSCs to neurons significantly increases epigenetic age, iPSC-neurons are still predicted as being fetal. Together our findings reiterate the need to better understand the limitations of existing epigenetic clocks for answering biological research questions and highlight a limitation of iPSC-neurons as a cellular model of age-related diseases.

ABSTRACTBiological aging is the gradual, progressive decline in system integrity that occurs with advancing chronological age, causing morbidity and disability. Measurements of the pace of aging are needed to serve as surrogate endpoints in trials of therapies designed to prevent disease by slowing biological aging. We report a blood DNA-methylation measure that is sensitive to variation in the pace of biological aging among individuals born the same year. We first modeled longitudinal change in 18 biomarkers tracking organ-system integrity across 12 years of follow-up in the Dunedin birth cohort. Rates of change in each biomarker were composited to form a measure of aging-related decline, termed Pace of Aging. Elastic-net regression was used to develop a DNA-methylation predictor of Pace of Aging, called DunedinPoAm for Dunedin (P)ace (o)f (A)ging (m)ethylation. Validation analyses showed DunedinPoAm was associated with functional decline in the Dunedin Study and more advanced biological age in the Understanding Society Study, predicted chronic disease and mortality in the Normative Aging Study, was accelerated by early-life adversity in the E-risk Study, and DunedinPoAm prediction was disrupted by caloric restriction in the CALERIE trial. DunedinPoAm generally outperformed epigenetic clocks. Findings provide proof-of-principle for DunedinPoAm as a single-time-point measure of a person’s pace of biological aging.

AbstractHuman DNA-methylation data have been used to develop biomarkers of ageing - referred to ‘epigenetic clocks’ - that have been widely used to identify differences between chronological age and biological age in health and disease including neurodegeneration, dementia and other brain phenotypes. Existing DNA methylation clocks are highly accurate in blood but are less precise when used in older samples or on brain tissue. We aimed to develop a novel epigenetic clock that performs optimally in human cortex tissue and has the potential to identify phenotypes associated with biological ageing in the brain. We generated an extensive dataset of human cortex DNA methylation data spanning the life-course (n = 1,397, ages = 1 to 104 years). This dataset was split into ‘training’ and ‘testing’ samples (training: n = 1,047; testing: n = 350). DNA methylation age estimators were derived using a transformed version of chronological age on DNA methylation at specific sites using elastic net regression, a supervised machine learning method. The cortical clock was subsequently validated in a novel human cortex dataset (n = 1,221, ages = 41 to 104 years) and tested for specificity in a large whole blood dataset (n = 1,175, ages = 28 to 98 years). We identified a set of 347 DNA methylation sites that, in combination optimally predict age in the human cortex. The sum of DNA methylation levels at these sites weighted by their regression coefficients provide the cortical DNA methylation clock age estimate. The novel clock dramatically out-performed previously reported clocks in additional cortical datasets. Our findings suggest that previous associations between predicted DNA methylation age and neurodegenerative phenotypes might represent false positives resulting from clocks not robustly calibrated to the tissue being tested and for phenotypes that become manifest in older ages. The age distribution and tissue type of samples included in training datasets need to be considered when building and applying epigenetic clock algorithms to human epidemiological or disease cohorts.

ABSTRACTWidespread structural, chemical and molecular differences have been reported between the male and female human brain. Although several neurodevelopmental disorders are more commonly diagnosed in males, little is known regarding sex differences in early human brain development. Here, we used RNA sequencing data from a large collection of human brain samples from the second trimester of gestation (N = 120) to assess sex biases in gene expression within the human fetal brain. In addition to 43 genes (102 Ensembl transcripts) transcribed from the Y-chromosome in males, we detected sex differences in the expression of 2558 autosomal genes (2723 Ensembl transcripts) and 155 genes on the X-chromosome (207 Ensembl transcripts) at a false discovery rate (FDR) < 0.1. Genes exhibiting sex-biased expression in human fetal brain are enriched for high-confidence risk genes for autism and other developmental disorders. Male-biased genes are enriched for expression in neural progenitor cells, whereas female-biased genes are enriched for expression in Cajal-Retzius cells and glia. All gene- and transcript-level data are provided as an online resource (available at http://fgen.psycm.cf.ac.uk/FBSeq1) through which researchers can search, download and visualize data pertaining to sex biases in gene expression during early human brain development.

Microarray technology has been used to measure genome-wide DNA methylation in thousands of individuals. These studies typically test the associations between individual DNA methylation sites ("probes") and complex traits or diseases. The results can be used to generate methylation profile scores (MPS) to predict outcomes in independent data sets. Although there are many parallels between MPS and polygenic (risk) scores (PGS), there are key differences. Here, we review motivations, methods, and applications of DNA methylation-based trait prediction, with a focus on common diseases. We contrast MPS with PGS, highlighting where assumptions made in genetic modeling may not hold in epigenetic data.

AbstractAlzheimer’s disease (AD) is a chronic neurodegenerative disease characterized by the progressive accumulation of amyloid-beta and neurofibrillary tangles of tau in the neocortex. We profiled DNA methylation in two regions of the cortex from 631 donors, performing an epigenome-wide association study of multiple measures of AD neuropathology. We meta-analyzed our results with those from previous studies of DNA methylation in AD cortex (total n = 2013 donors), identifying 334 cortical differentially methylated positions (DMPs) associated with AD pathology including methylomic variation at loci not previously implicated in dementia. We subsequently profiled DNA methylation in NeuN+ (neuronal-enriched), SOX10+ (oligodendrocyte-enriched) and NeuN–/SOX10– (microglia- and astrocyte-enriched) nuclei, finding that the majority of DMPs identified in ‘bulk’ cortex tissue reflect DNA methylation differences occurring in non-neuronal cells. Our study highlights the power of utilizing multiple measures of neuropathology to identify epigenetic signatures of AD and the importance of characterizing disease-associated variation in purified cell-types.

BACKGROUND: Measures to quantify changes in the pace of biological aging in response to intervention are needed to evaluate geroprotective interventions for humans. Previously, we showed that quantification of the pace of biological aging from a DNA-methylation blood test was possible (Belsky et al. 2020). Here, we report a next-generation DNA-methylation biomarker of Pace of Aging, DunedinPACE (for Pace of Aging Calculated from the Epigenome). METHODS: We used data from the Dunedin Study 1972-1973 birth cohort tracking within-individual decline in 19 indicators of organ-system integrity across four time points spanning two decades to model Pace of Aging. We distilled this two-decade Pace of Aging into a single-time-point DNA-methylation blood-test using elastic-net regression and a DNA-methylation dataset restricted to exclude probes with low test-retest reliability. We evaluated the resulting measure, named DunedinPACE, in five additional datasets. RESULTS: DunedinPACE showed high test-retest reliability, was associated with morbidity, disability, and mortality, and indicated faster aging in young adults with childhood adversity. DunedinPACE effect-sizes were similar to GrimAge Clock effect-sizes. In analysis of incident morbidity, disability, and mortality, DunedinPACE and added incremental prediction beyond GrimAge. CONCLUSIONS: DunedinPACE is a novel blood biomarker of the pace of aging for gerontology and geroscience. FUNDING: This research was supported by US-National Institute on Aging grants AG032282, AG061378, AG066887, and UK Medical Research Council grant MR/P005918/1.

AbstractParkinson’s disease (PD) and dementia with Lewy bodies (DLB) are closely related progressive disorders with no available disease-modifying therapy, neuropathologically characterized by intraneuronal aggregates of misfolded α-synuclein. To explore the role of DNA methylation changes in PD and DLB pathogenesis, we performed an epigenome-wide association study (EWAS) of 322 postmortem frontal cortex samples and replicated results in an independent set of 200 donors. We report novel differentially methylated replicating loci associated with Braak Lewy body stage near TMCC2, SFMBT2, AKAP6 and PHYHIP. Differentially methylated probes were independent of known PD genetic risk alleles. Meta-analysis provided suggestive evidence for a differentially methylated locus within the chromosomal region affected by the PD-associated 22q11.2 deletion. Our findings elucidate novel disease pathways in PD and DLB and generate hypotheses for future molecular studies of Lewy body pathology.

Most epigenetic epidemiology to date has utilized microarrays to identify positions in the genome where variation in DNA methylation is associated with environmental exposures or disease. However, these profile less than 3% of DNA methylation sites in the human genome, potentially missing affected loci and preventing the discovery of disrupted biological pathways. Third generation sequencing technologies, including Nanopore sequencing, have the potential to revolutionize the generation of epigenetic data, not only by providing genuine genome-wide coverage but profiling epigenetic modifications direct from native DNA. Here we assess the viability of using Nanopore sequencing for epidemiology by performing a comparison with DNA methylation quantified using the most comprehensive microarray available, the Illumina EPIC array. We implemented a CRISPR-Cas9 targeted sequencing approach in concert with Nanopore sequencing to profile DNA methylation in three genomic regions to attempt to rediscover genomic positions that existing technologies have shown are differentially methylated in tobacco smokers. Using Nanopore sequencing reads, DNA methylation was quantified at 1779 CpGs across three regions, providing a finer resolution of DNA methylation patterns compared to the EPIC array. The correlation of estimated levels of DNA methylation between platforms was high. Furthermore, we identified 12 CpGs where hypomethylation was significantly associated with smoking status, including 10 within the AHRR gene. In summary, Nanopore sequencing is a valid option for identifying genomic loci where large differences in DNAm are associated with a phenotype and has the potential to advance our understanding of the role differential methylation plays in the etiology of complex disease.

BACKGROUND: There is growing interest in the role of DNA methylation in regulating the transcription of mitochondrial genes, particularly in brain disorders characterized by mitochondrial dysfunction. Here, we present a novel approach to interrogate the mitochondrial DNA methylome at single base resolution using targeted bisulfite sequencing. We applied this method to investigate mitochondrial DNA methylation patterns in post-mortem superior temporal gyrus and cerebellum brain tissue from seven human donors. RESULTS: We show that mitochondrial DNA methylation patterns are relatively low but conserved, with peaks in DNA methylation at several sites, such as within the D-LOOP and the genes MT-ND2, MT-ATP6, MT-ND4, MT-ND5 and MT-ND6, predominantly in a non-CpG context. The elevated DNA methylation we observe in the D-LOOP we validate using pyrosequencing. We identify loci that show differential DNA methylation patterns associated with age, sex and brain region. Finally, we replicate previously reported differentially methylated regions between brain regions from a methylated DNA immunoprecipitation sequencing study. CONCLUSIONS: We have annotated patterns of DNA methylation at single base resolution across the mitochondrial genome in human brain samples. Looking to the future this approach could be utilized to investigate the role of mitochondrial epigenetic mechanisms in disorders that display mitochondrial dysfunction.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with an estimated heritability between 40 and 50%. DNA methylation patterns can serve as proxies of (past) exposures and disease progression, as well as providing a potential mechanism that mediates genetic or environmental risk. Here, we present a blood-based epigenome-wide association study meta-analysis in 9706 samples passing stringent quality control (6763 patients, 2943 controls). We identified a total of 45 differentially methylated positions (DMPs) annotated to 42 genes, which are enriched for pathways and traits related to metabolism, cholesterol biosynthesis, and immunity. We then tested 39 DNA methylation–based proxies of putative ALS risk factors and found that high-density lipoprotein cholesterol, body mass index, white blood cell proportions, and alcohol intake were independently associated with ALS. Integration of these results with our latest genome-wide association study showed that cholesterol biosynthesis was potentially causally related to ALS. Last, DNA methylation at several DMPs and blood cell proportion estimates derived from DNA methylation data were associated with survival rate in patients, suggesting that they might represent indicators of underlying disease processes potentially amenable to therapeutic interventions.

Amyotrophic lateral sclerosis (ALS) is a fatal motoneuron disease with a monogenic cause in approximately 10% of cases. However, familial clustering of disease without inheritance in a Mendelian manner and the broad range of phenotypes suggest the presence of epigenetic mechanisms. Hence, we performed an epigenome-wide association study on sporadic, symptomatic and presymptomatic familial ALS cases with mutations in C9ORF72 and FUS and healthy controls studying DNA methylation in blood cells. We found differentially methylated DNA positions (DMPs) and regions embedding DMPs associated with either disease status, C9ORF72 or FUS mutation status. One DMP reached methylome-wide significance and is attributed to a region encoding a long non-coding RNA (LOC389247). Furthermore, we could demonstrate co-localization of DMPs with an ALS-associated GWAS region near the SCN7A/SCN9A and XIRP2 genes. Finally, a classifier model that predicts disease status (ALS, healthy) classified all but one presymptomatic mutation carrier as healthy, suggesting that the presence of ALS symptoms rather than the presence of ALS-associated genetic mutations is associated with blood cell DNA methylation.

Disadvantaged socio-economic position (SEP) is associated with greater biological age, relative to chronological age, measured by DNA methylation (positive 'age acceleration', AA). Social mobility has been proposed to ameliorate health inequalities. This study aimed to understand the association of social mobility with positive AA. Diagonal reference modelling and ordinary least square regression techniques were applied to explore social mobility and four measures of age acceleration (first-generation: 'Horvath', 'Hannum' and second-generation: 'Phenoage', DunedinPoAm) in n = 3140 participants of the UK Household Longitudinal Study. Disadvantaged SEP in early life is associated with positive AA for three (Hannum, Phenoage and DunedinPoAm) of the four measures examined while the second generation biomarkers are associated with SEP in adulthood (p 

AbstractEpigenetic dysregulation is thought to contribute to the etiology of schizophrenia (SZ), but the cell type-specificity of DNA methylation makes population-based epigenetic studies of SZ challenging. To train an SZ case–control classifier based on DNA methylation in blood, therefore, we focused on human genomic regions of systemic interindividual epigenetic variation (CoRSIVs), a subset of which are represented on the Illumina Human Methylation 450K (HM450) array. HM450 DNA methylation data on whole blood of 414 SZ cases and 433 non-psychiatric controls were used as training data for a classification algorithm with built-in feature selection, sparse partial least squares discriminate analysis (SPLS-DA); application of SPLS-DA to HM450 data has not been previously reported. Using the first two SPLS-DA dimensions we calculated a “risk distance” to identify individuals with the highest probability of SZ. The model was then evaluated on an independent HM450 data set on 353 SZ cases and 322 non-psychiatric controls. Our CoRSIV-based model classified 303 individuals as cases with a positive predictive value (PPV) of 80%, far surpassing the performance of a model based on polygenic risk score (PRS). Importantly, risk distance (based on CoRSIV methylation) was not associated with medication use, arguing against reverse causality. Risk distance and PRS were positively correlated (Pearson r = 0.28, P = 1.28 × 10−12), and mediational analysis suggested that genetic effects on SZ are partially mediated by altered methylation at CoRSIVs. Our results indicate two innate dimensions of SZ risk: one based on genetic, and the other on systemic epigenetic variants.

Epigenome-wide association studies of Alzheimer's disease have highlighted neuropathology-associated DNA methylation differences, although existing studies have been limited in sample size and utilized different brain regions. Here, we combine data from six DNA methylomic studies of Alzheimer's disease (N = 1453 unique individuals) to identify differential methylation associated with Braak stage in different brain regions and across cortex. We identify 236 CpGs in the prefrontal cortex, 95 CpGs in the temporal gyrus and ten CpGs in the entorhinal cortex at Bonferroni significance, with none in the cerebellum. Our cross-cortex meta-analysis (N = 1408 donors) identifies 220 CpGs associated with neuropathology, annotated to 121 genes, of which 84 genes have not been previously reported at this significance threshold. We have replicated our findings using two further DNA methylomic datasets consisting of a further >600 unique donors. The meta-analysis summary statistics are available in our online data resource ( www.epigenomicslab.com/ad-meta-analysis/ ).

Most epigenome-wide association studies (EWAS) quantify DNA methylation (DNAm) in peripheral tissues such as whole blood to identify positions in the genome where variation is statistically associated with a trait or exposure. As whole blood comprises a mix of cell types, it is unclear whether trait-associated DNAm variation is specific to an individual cellular population. We collected three peripheral tissues (whole blood, buccal epithelial and nasal epithelial cells) from thirty individuals. Whole blood samples were subsequently processed using fluorescence-activated cell sorting (FACS) to purify five constituent cell-types (monocytes, granulocytes, CD4+ T cells, CD8+ T cells, and B cells). DNAm was profiled in all eight sample-types from each individual using the Illumina EPIC array. We identified significant differences in both the level and variability of DNAm between different sample types, and DNAm data-derived estimates of age and smoking were found to differ dramatically across sample types from the same individual. We found that for the majority of loci variation in DNAm in individual blood cell types was only weakly predictive of variance in DNAm measured in whole blood, although the proportion of variance explained was greater than that explained by either buccal or nasal epithelial samples. Covariation across sample types was much higher for DNAm sites influenced by genetic factors. Overall, we observe that DNAm variation in whole blood is additively influenced by a combination of the major blood cell types. For a subset of sites, however, variable DNAm detected in whole blood can be attributed to variation in a single blood cell type providing potential mechanistic insight about EWAS findings. Our results suggest that associations between whole blood DNAm and traits or exposures reflect differences in multiple cell types and our data will facilitate the interpretation of findings in epigenetic epidemiology.

BACKGROUND: the combination of sodium bisulfite treatment with highly-parallel sequencing is a common method for quantifying DNA methylation across the genome. The power to detect between-group differences in DNA methylation using bisulfite-sequencing approaches is influenced by both experimental (e.g. read depth, missing data and sample size) and biological (e.g. mean level of DNA methylation and difference between groups) parameters. There is, however, no consensus about the optimal thresholds for filtering bisulfite sequencing data with implications for the reproducibility of findings in epigenetic epidemiology. RESULTS: We used a large reduced representation bisulfite sequencing (RRBS) dataset to assess the distribution of read depth across DNA methylation sites and the extent of missing data. To investigate how various study variables influence power to identify DNA methylation differences between groups, we developed a framework for simulating bisulfite sequencing data. As expected, sequencing read depth, group size, and the magnitude of DNA methylation difference between groups all impacted upon statistical power. The influence on power was not dependent on one specific parameter, but reflected the combination of study-specific variables. As a resource to the community, we have developed a tool, POWEREDBiSeq, which utilizes our simulation framework to predict study-specific power for the identification of DNAm differences between groups, taking into account user-defined read depth filtering parameters and the minimum sample size per group. CONCLUSIONS: Our data-driven approach highlights the importance of filtering bisulfite-sequencing data by minimum read depth and illustrates how the choice of threshold is influenced by the specific study design and the expected differences between groups being compared. The POWEREDBiSeq tool, which can be applied to different types of bisulfite sequencing data (e.g. RRBS, whole genome bisulfite sequencing (WGBS), targeted bisulfite sequencing and amplicon-based bisulfite sequencing), can help users identify the level of data filtering needed to optimize power and aims to improve the reproducibility of bisulfite sequencing studies.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with a lifetime risk of one in 350 people and an unmet need for disease-modifying therapies. We conducted a cross-ancestry genome-wide association study (GWAS) including 29,612 patients with ALS and 122,656 controls, which identified 15 risk loci. When combined with 8,953 individuals with whole-genome sequencing (6,538 patients, 2,415 controls) and a large cortex-derived expression quantitative trait locus (eQTL) dataset (MetaBrain), analyses revealed locus-specific genetic architectures in which we prioritized genes either through rare variants, short tandem repeats or regulatory effects. ALS-associated risk loci were shared with multiple traits within the neurodegenerative spectrum but with distinct enrichment patterns across brain regions and cell types. of the environmental and lifestyle risk factors obtained from the literature, Mendelian randomization analyses indicated a causal role for high cholesterol levels. The combination of all ALS-associated signals reveals a role for perturbations in vesicle-mediated transport and autophagy and provides evidence for cell-autonomous disease initiation in glutamatergic neurons.

We performed a systematic analysis of blood DNA methylation profiles from 4483 participants from seven independent cohorts identifying differentially methylated positions (DMPs) associated with psychosis, schizophrenia, and treatment-resistant schizophrenia. Psychosis cases were characterized by significant differences in measures of blood cell proportions and elevated smoking exposure derived from the DNA methylation data, with the largest differences seen in treatment-resistant schizophrenia patients. We implemented a stringent pipeline to meta-analyze epigenome-wide association study (EWAS) results across datasets, identifying 95 DMPs associated with psychosis and 1048 DMPs associated with schizophrenia, with evidence of colocalization to regions nominated by genetic association studies of disease. Many schizophrenia-associated DNA methylation differences were only present in patients with treatment-resistant schizophrenia, potentially reflecting exposure to the atypical antipsychotic clozapine. Our results highlight how DNA methylation data can be leveraged to identify physiological (e.g. differential cell counts) and environmental (e.g. smoking) factors associated with psychosis and molecular biomarkers of treatment-resistant schizophrenia.

Accumulating evidence suggests that individuals exposed to victimization at key developmental stages may have different epigenetic fingerprints compared to those exposed to no/minimal stressful events, however results are inconclusive. This study aimed to strengthen causal inference regarding the impact of adolescent victimization on the epigenome by controlling for genetic variation, age, gender, and shared environmental exposures. We conducted longitudinal epigenome-wide association analyses (EWAS) on DNA methylation (DNAm) profiles of 118 monozygotic (MZ) twin pairs from the Environmental Risk study with and without severe adolescent victimization generated using buccal DNA collected at ages 5, 10 and 18, and the Illumina EPIC array. Additionally, we performed cross-sectional EWAS on age-18 blood and buccal DNA from the same individuals to elucidate tissue-specific signatures of severe adolescent victimization. Our analyses identified 20 suggestive differentially methylated positions (DMPs) (P

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.

BACKGROUND: Sex is an important covariate of epigenome-wide association studies due to its strong influence on DNA methylation patterns across numerous genomic positions. Nevertheless, many samples on the Gene Expression Omnibus (GEO) frequently lack a sex annotation or are incorrectly labelled. Considering the influence that sex imposes on DNA methylation patterns, it is necessary to ensure that methods for filtering poor samples and checking of sex assignment are accurate and widely applicable. RESULTS: Here we presented a novel method to predict sex using only DNA methylation beta values, which can be readily applied to almost all DNA methylation datasets of different formats (raw IDATs or text files with only signal intensities) uploaded to GEO. We identified 4345 significantly (p

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype-phenotype map than previously anticipated.

AbstractMonozygotic (MZ) twins and higher-order multiples arise when a zygote splits during pre-implantation stages of development. The mechanisms underpinning this event have remained a mystery. Because MZ twinning rarely runs in families, the leading hypothesis is that it occurs at random. Here, we show that MZ twinning is strongly associated with a stable DNA methylation signature in adult somatic tissues. This signature spans regions near telomeres and centromeres, Polycomb-repressed regions and heterochromatin, genes involved in cell-adhesion, WNT signaling, cell fate, and putative human metastable epialleles. Our study also demonstrates a never-anticipated corollary: because identical twins keep a lifelong molecular signature, we can retrospectively diagnose if a person was conceived as monozygotic twin.

We evaluated whether the association between cigarette smoking and dementia risk is modified by genetic predisposition including apolipoprotein E (APOE) genotype and polygenic risk (excluding the APOE region). We included 193,198 UK Biobank participants aged 60-73 years without dementia at baseline. of non-APOE-ε4 carriers, 0.89% (95% CI 0.73-1.08%) current smokers developed dementia compared with 0.49% (95% CI 0.44-0.55%) of never smokers (adjusted HR 1.78; 95% CI 1.39-2.29). In contrast, of one APOE-ε4 allele carriers, 1.69% (95% CI 1.31-2.12%) current smokers developed dementia compared with 1.40% (95% CI 1.25-1.55%) of never smokers (adjusted HR 1.06; 95% CI 0.77-1.45); of two APOE-ε4 alleles carriers, 4.90% (95% CI 2.92-7.61%) current smokers developed dementia compared with 3.87% (95% CI 3.11-4.74%) of never smokers (adjusted HR 0.94; 95% CI 0.49-1.79). of participants with high polygenic risk, 1.77% (95% CI 1.35-2.27%) current smokers developed dementia compared with 1.05% (95% CI 0.91-1.21%) of never smokers (adjusted HR 1.63; 95% CI 1.16-2.28). A significant interaction was found between APOE genotype and smoking status (P = 0.002) while no significant interaction was identified between polygenic risk and smoking status (P = 0.25). APOE genotype but not polygenic risk modified the effect of smoking on dementia risk.

BACKGROUND: People with neurodegenerative disorders show diverse clinical syndromes, genetic heterogeneity, and distinct brain pathological changes, but studies report overlap between these features. DNA methylation (DNAm) provides a way to explore this overlap and heterogeneity as it is determined by the combined effects of genetic variation and the environment. In this study, we aim to identify shared blood DNAm differences between controls and people with Alzheimer's disease, amyotrophic lateral sclerosis, and Parkinson's disease. RESULTS: We use a mixed-linear model method (MOMENT) that accounts for the effect of (un)known confounders, to test for the association of each DNAm site with each disorder. While only three probes are found to be genome-wide significant in each MOMENT association analysis of amyotrophic lateral sclerosis and Parkinson's disease (and none with Alzheimer's disease), a fixed-effects meta-analysis of the three disorders results in 12 genome-wide significant differentially methylated positions. Predicted immune cell-type proportions are disrupted across all neurodegenerative disorders. Protein inflammatory markers are correlated with profile sum-scores derived from disease-associated immune cell-type proportions in a healthy aging cohort. In contrast, they are not correlated with MOMENT DNAm-derived profile sum-scores, calculated using effect sizes of the 12 differentially methylated positions as weights. CONCLUSIONS: We identify shared differentially methylated positions in whole blood between neurodegenerative disorders that point to shared pathogenic mechanisms. These shared differentially methylated positions may reflect causes or consequences of disease, but they are unlikely to reflect cell-type proportion differences.

Induced pluripotent stem cells (iPSCs) and their differentiated neurons (iPSC-neurons) are a widely used cellular model in the research of the central nervous system. However, it is unknown how well they capture age-associated processes, particularly given that pluripotent cells are only present during the earliest stages of mammalian development. Epigenetic clocks utilize coordinated age-associated changes in DNA methylation to make predictions that correlate strongly with chronological age. It has been shown that the induction of pluripotency rejuvenates predicted epigenetic age. As existing clocks are not optimized for the study of brain development, we developed the fetal brain clock (FBC), a bespoke epigenetic clock trained in human prenatal brain samples in order to investigate more precisely the epigenetic age of iPSCs and iPSC-neurons. The FBC was tested in two independent validation cohorts across a total of 194 samples, confirming that the FBC outperforms other established epigenetic clocks in fetal brain cohorts. We applied the FBC to DNA methylation data from iPSCs and embryonic stem cells and their derived neuronal precursor cells and neurons, finding that these cell types are epigenetically characterized as having an early fetal age. Furthermore, while differentiation from iPSCs to neurons significantly increases epigenetic age, iPSC-neurons are still predicted as being fetal. Together our findings reiterate the need to better understand the limitations of existing epigenetic clocks for answering biological research questions and highlight a limitation of iPSC-neurons as a cellular model of age-related diseases.

Common genetic variation appears to largely influence risk for neuropsychiatric disorders through effects on gene regulation. It is therefore possible to shed light on the biology of these conditions by testing for enrichment of associated genetic variation within regulatory genomic regions operating in specific tissues or cell types. Here, we have used the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-Seq) to map open chromatin (an index of active regulatory genomic regions) in bulk tissue, NeuN+ and NeuN- nuclei from the prenatal human frontal cortex, and tested enrichment of single-nucleotide polymorphism (SNP) heritability for five neuropsychiatric disorders (autism spectrum disorder, attention deficit hyperactivity disorder [ADHD], bipolar disorder, major depressive disorder, and schizophrenia) within these regions. We observed significant enrichment of SNP heritability for ADHD, major depressive disorder, and schizophrenia within open chromatin regions (OCRs) mapped in bulk fetal frontal cortex, and for all five tested neuropsychiatric conditions when we restricted these sites to those overlapping histone modifications indicative of enhancers (H3K4me1) or promoters (H3K4me3) in fetal brain. SNP heritability for neuropsychiatric disorders was significantly enriched in OCRs identified in fetal frontal cortex NeuN- as well as NeuN+ nuclei overlapping fetal brain H3K4me1 or H3K4me3 sites. We additionally demonstrate the utility of our mapped OCRs for prioritizing potentially functional SNPs at genome-wide significant risk loci for neuropsychiatric disorders. Our data provide evidence for an early neurodevelopmental component to a range of neuropsychiatric conditions and highlight an important role for regulatory genomic regions active within both NeuN+ and NeuN- cells of the prenatal brain.

A growing number of epigenome-wide association studies have demonstrated a role for DNA methylation in the brain in Alzheimer's disease. With the aim of exploring peripheral biomarker potential, we have examined DNA methylation patterns in whole blood collected from 284 individuals in the AddNeuroMed study, which included 89 nondemented controls, 86 patients with Alzheimer's disease, and 109 individuals with mild cognitive impairment, including 38 individuals who progressed to Alzheimer's disease within 1 year. We identified significant differentially methylated regions, including 12 adjacent hypermethylated probes in the HOXB6 gene in Alzheimer's disease, which we validated using pyrosequencing. Using weighted gene correlation network analysis, we identified comethylated modules of genes that were associated with key variables such as APOE genotype and diagnosis. In summary, this study represents the first large-scale epigenome-wide association study of Alzheimer's disease and mild cognitive impairment using blood. We highlight the differences in various loci and pathways in early disease, suggesting that these patterns relate to cognitive decline at an early stage.

Clozapine is an atypical antipsychotic medication that is used to treat schizophrenia patients who are resistant to other antipsychotic drugs. The molecular mechanisms mediating the effects of clozapine are not well understood and its use is often associated with severe side-effects. In this study, we exposed groups of wild-type zebrafish to two doses of clozapine ('low' (20 µg/L) and 'high' (70 µg/L)) over a 72-h period, observing dose-dependent effects on behaviour. Using RNA sequencing (RNA-seq) we identified multiple genes differentially expressed in the zebrafish brain following exposure to clozapine. Network analysis identified co-expression modules characterised by striking changes in module connectivity in response to clozapine, and these were enriched for regulatory pathways relevant to the etiology of schizophrenia. Our study highlights the utility of zebrafish as a model for assessing the molecular consequences of antipsychotic medications and identifies genomic networks potentially involved in schizophrenia.

Abstract
. Illumina DNA methylation arrays are a widely used tool for performing genome-wide DNA methylation analyses. However, measurements obtained from these arrays may be affected by technical artefacts that result in spurious associations if left unchecked. Cross-reactivity represents one of the major challenges, meaning that probes may map to multiple regions in the genome. Although several studies have reported on this issue, few studies have empirically examined the impact of cross-reactivity in an epigenome-wide association study (EWAS). In this paper, we report on cross-reactivity issues that we discovered in a large EWAS on the presence of the C9orf72 repeat expansion in ALS patients. Specifically, we found that that the majority of the significant probes inadvertently cross-hybridized to the C9orf72 locus. Importantly, these probes were not flagged as cross-reactive in previous studies, leading to novel insights into the extent to which cross-reactivity can impact EWAS. Our findings are particularly relevant for epigenetic studies into diseases associated with repeat expansions and other types of structural variation. More generally however, considering that most spurious associations were not excluded based on pre-defined sets of cross-reactive probes, we believe that the presented data-driven flag and consider approach is relevant for any type of EWAS.

AbstractAttention-deficit/hyperactivity disorder (ADHD) is a highly heritable neurodevelopmental disorder that often persists into adulthood. There is growing evidence that epigenetic dysregulation participates in ADHD. Given that only a limited number of epigenome-wide association studies (EWASs) of ADHD have been conducted so far and they have mainly focused on pediatric and population-based samples, we performed an EWAS in a clinical sample of adults with ADHD. We report one CpG site and four regions differentially methylated between patients and controls, which are located in or near genes previously involved in autoimmune diseases, cancer or neuroticism. Our sensitivity analyses indicate that smoking status is not responsible for these results and that polygenic risk burden for ADHD does not greatly impact the signatures identified. Additionally, we show an overlap of our EWAS findings with genetic signatures previously described for ADHD and with epigenetic signatures for smoking behavior and maternal smoking. These findings support a role of DNA methylation in ADHD and emphasize the need for additional efforts in larger samples to clarify the role of epigenetic mechanisms on ADHD across the lifespan.

Abstract
. Alzheimer’s disease is a highly heritable, common neurodegenerative disease characterized neuropathologically by the accumulation of β-amyloid plaques and tau-containing neurofibrillary tangles. In addition to the well-established risk associated with the APOE locus, there has been considerable success in identifying additional genetic variants associated with Alzheimer’s disease. Major challenges in understanding how genetic risk influences the development of Alzheimer’s disease are clinical and neuropathological heterogeneity, and the high level of accompanying comorbidities. We report a multimodal analysis integrating longitudinal clinical and cognitive assessment with neuropathological data collected as part of the Brains for Dementia Research study to understand how genetic risk factors for Alzheimer’s disease influence the development of neuropathology and clinical performance. Six hundred and ninety-three donors in the Brains for Dementia Research cohort with genetic data, semi-quantitative neuropathology measurements, cognitive assessments and established diagnostic criteria were included in this study. We tested the association of APOE genotype and Alzheimer’s disease polygenic risk score—a quantitative measure of genetic burden—with survival, four common neuropathological features in Alzheimer’s disease brains (neurofibrillary tangles, β-amyloid plaques, Lewy bodies and transactive response DNA-binding protein 43 proteinopathy), clinical status (clinical dementia rating) and cognitive performance (Mini-Mental State Exam, Montreal Cognitive Assessment). The APOE ε4 allele was significantly associated with younger age of death in the Brains for Dementia Research cohort. Our analyses of neuropathology highlighted two independent pathways from APOE ε4, one where β-amyloid accumulation co-occurs with the development of tauopathy, and a second characterized by direct effects on tauopathy independent of β-amyloidosis. Although we also detected association between APOE ε4 and dementia status and cognitive performance, these were all mediated by tauopathy, highlighting that they are a consequence of the neuropathological changes. Analyses of polygenic risk score identified associations with tauopathy and β-amyloidosis, which appeared to have both shared and unique contributions, suggesting that different genetic variants associated with Alzheimer’s disease affect different features of neuropathology to different degrees. Taken together, our results provide insight into how genetic risk for Alzheimer’s disease influences both the clinical and pathological features of dementia, increasing our understanding about the interplay between APOE genotype and other genetic risk factors.

Attention-deficit/hyperactivity disorder (ADHD) is a neurodevelopmental disorder that often persists into adulthood. ADHD and related personality traits, such as impulsivity and callousness, are caused by genetic and environmental factors and their interplay. Epigenetic modifications of DNA, including methylation, are thought to mediate between such factors and behavior and may behave as biomarkers for disorders. Here, we set out to study DNA methylation in persistent ADHD and related traits. We performed epigenome-wide association studies (EWASs) on peripheral whole blood from participants in the NeuroIMAGE study (age range 12-23 years). We compared participants with persistent ADHD (n = 35) with healthy controls (n = 19) and with participants with remittent ADHD (n = 19). Additionally, we performed EWASs of impulsive and callous traits derived from the Conners Parent Rating Scale and the Callous-Unemotional Inventory, respectively, across all participants. For every EWAS, the linear regression model analyzed included covariates for age, sex, smoking scores, and surrogate variables reflecting blood cell type composition and genetic background. We observed no epigenome-wide significant differences in single CpG site methylation between participants with persistent ADHD and healthy controls or participants with remittent ADHD. However, epigenome-wide analysis of differentially methylated regions provided significant findings showing that hypermethylated regions in the APOB and LPAR5 genes were associated with ADHD persistence compared to ADHD remittance (p = 1.68. 10-24 and p = 9.06. 10-7, respectively); both genes are involved in cholesterol signaling. Both findings appeared to be linked to genetic variation in cis. We found neither significant epigenome-wide single CpG sites nor regions associated with impulsive and callous traits; the top-hits from these analyses were annotated to genes involved in neurotransmitter release and the regulation of the biological clock. No link to genetic variation was observed for these findings, which thus might reflect environmental influences. In conclusion, in this pilot study with a small sample size, we observed several DNA-methylation-disorder/trait associations of potential significance for ADHD and the related behavioral traits. Although we do not wish to draw conclusions before replication in larger, independent samples, cholesterol signaling and metabolism may be of relevance for the onset and/or persistence of ADHD.

AbstractSuicide is the second leading cause of death globally among young people representing a significant global health burden. Although the molecular correlates of suicide remains poorly understood, it has been hypothesised that epigenomic processes may play a role. The objective of this study was to identify suicide-associated DNA methylation changes in the human brain by utilising previously published and unpublished methylomic datasets. We analysed prefrontal cortex (PFC, n = 211) and cerebellum (CER, n = 114) DNA methylation profiles from suicide completers and non-psychiatric, sudden-death controls, meta-analysing data from independent cohorts for each brain region separately. We report evidence for altered DNA methylation at several genetic loci in suicide cases compared to controls in both brain regions with suicide-associated differentially methylated positions enriched among functional pathways relevant to psychiatric phenotypes and suicidality, including nervous system development (PFC) and regulation of long-term synaptic depression (CER). In addition, we examined the functional consequences of variable DNA methylation within a PFC suicide-associated differentially methylated region (PSORS1C3 DMR) using a dual luciferase assay and examined expression of nearby genes. DNA methylation within this region was associated with decreased expression of firefly luciferase but was not associated with expression of nearby genes, PSORS1C3 and POU5F1. Our data suggest that suicide is associated with DNA methylation, offering novel insights into the molecular pathology associated with suicidality.

Methylome-wide association studies (MWASs) are promising complements to sequence variation studies. We used existing sequencing-based methylation data, which assayed the majority of all 28 million CpGs in the human genome, to perform an MWAS for schizophrenia in blood, while controlling for cell-type heterogeneity with a recently generated platform-specific reference panel. Next, we compared the MWAS results with findings from 3 existing large-scale array-based schizophrenia methylation studies in blood that assayed up to ~450 000 CpGs. Our MWAS identified 22 highly significant loci (P < 5 × 10-8) and 852 suggestively significant loci (P < 1 × 10-5). The top finding (P = 5.62 × 10-11, q = 0.001) was located in MFN2, which encodes mitofusin-2 that regulates Ca2+ transfer from the endoplasmic reticulum to mitochondria in cooperation with DISC1. The second-most significant site (P = 1.38 × 10-9, q = 0.013) was located in ALDH1A2, which encodes an enzyme for astrocyte-derived retinoic acid-a key neuronal morphogen with relevance for schizophrenia. Although the most significant MWAS findings were not assayed on the arrays, we observed significant enrichment of overlapping findings with 2 of the 3 array datasets (P = 0.0315, 0.0045, 0.1946). Overrepresentation analysis of Gene Ontology terms for the genes in the significant overlaps suggested high similarity in the biological functions detected by the different datasets. Top terms were related to immune and/or stress responses, cell adhesion and motility, and a broad range of processes essential for neurodevelopment.

Surgery is an invasive procedure evoking acute inflammatory and immune responses that can influence risk for postoperative complications including cognitive dysfunction and delirium. Although the specific mechanisms driving these responses have not been well-characterized, they are hypothesized to involve the epigenetic regulation of gene expression. We quantified genome-wide levels of DNA methylation in peripheral blood mononuclear cells (PBMCs) longitudinally collected from a cohort of elderly patients undergoing major surgery, comparing samples collected at baseline to those collected immediately post-operatively and at discharge from hospital. We identified acute changes in measured DNA methylation at sites annotated to immune system genes, paralleling changes in serum-levels of markers including C-reactive protein (CRP) and Interleukin 6 (IL-6) measured in the same individuals. Many of the observed changes in measured DNA methylation were consistent across different types of major surgery, although there was notable heterogeneity between surgery types at certain loci. The acute changes in measured DNA methylation induced by surgery are relatively stable in the post-operative period, generally persisting until discharge from hospital. Our results highlight the dramatic alterations in gene regulation induced by invasive surgery, primarily reflecting upregulation of the immune system in response to trauma, wound healing and anaesthesia.

DNA methylation plays an important role in both normal human development and risk of disease. The most utilized method of assessing DNA methylation uses BeadChips, generating an epigenome-wide "snapshot" of >450,000 observations (probe measurements) per assay. However, the reliability of each of these measurements is not equal, and little consideration is paid to consequences for research. We correlated repeat measurements of the same DNA samples using the Illumina HumanMethylation450K and the Infinium MethylationEPIC BeadChips in 350 blood DNA samples. Probes that were reliably measured were more heritable and showed consistent associations with environmental exposures, gene expression, and greater cross-tissue concordance. Unreliable probes were less replicable and generated an unknown volume of false negatives. This serves as a lesson for working with DNA methylation data, but the lessons are equally applicable to working with other data: as we advance toward generating increasingly greater volumes of data, failure to document reliability risks harming reproducibility.

Biological aging is the gradual, progressive decline in system integrity that occurs with advancing chronological age, causing morbidity and disability. Measurements of the pace of aging are needed as surrogate endpoints in trials of therapies designed to prevent disease by slowing biological aging. We report a blood-DNA-methylation measure that is sensitive to variation in pace of biological aging among individuals born the same year. We first modeled change-over-time in 18 biomarkers tracking organ-system integrity across 12 years of follow-up in n = 954 members of the Dunedin Study born in 1972–1973. Rates of change in each biomarker over ages 26–38 years were composited to form a measure of aging-related decline, termed Pace-of-Aging. Elastic-net regression was used to develop a DNA-methylation predictor of Pace-of-Aging, called DunedinPoAm for Dunedin(P)ace(o)f(A)ging(m)ethylation. Validation analysis in cohort studies and the CALERIE trial provide proof-of-principle for DunedinPoAm as a single-time-point measure of a person’s pace of biological aging.

Human DNA methylation data have been used to develop biomarkers of ageing, referred to as 'epigenetic clocks', which have been widely used to identify differences between chronological age and biological age in health and disease including neurodegeneration, dementia and other brain phenotypes. Existing DNA methylation clocks have been shown to be highly accurate in blood but are less precise when used in older samples or in tissue types not included in training the model, including brain. We aimed to develop a novel epigenetic clock that performs optimally in human cortex tissue and has the potential to identify phenotypes associated with biological ageing in the brain. We generated an extensive dataset of human cortex DNA methylation data spanning the life course (n = 1397, ages = 1 to 108 years). This dataset was split into 'training' and 'testing' samples (training: n = 1047; testing: n = 350). DNA methylation age estimators were derived using a transformed version of chronological age on DNA methylation at specific sites using elastic net regression, a supervised machine learning method. The cortical clock was subsequently validated in a novel independent human cortex dataset (n = 1221, ages = 41 to 104 years) and tested for specificity in a large whole blood dataset (n = 1175, ages = 28 to 98 years). We identified a set of 347 DNA methylation sites that, in combination, optimally predict age in the human cortex. The sum of DNA methylation levels at these sites weighted by their regression coefficients provide the cortical DNA methylation clock age estimate. The novel clock dramatically outperformed previously reported clocks in additional cortical datasets. Our findings suggest that previous associations between predicted DNA methylation age and neurodegenerative phenotypes might represent false positives resulting from clocks not robustly calibrated to the tissue being tested and for phenotypes that become manifest in older ages. The age distribution and tissue type of samples included in training datasets need to be considered when building and applying epigenetic clock algorithms to human epidemiological or disease cohorts.

We conducted DNA methylation association analyses using Illumina 450K data from whole blood for an Australian amyotrophic lateral sclerosis (ALS) case–control cohort (782 cases and 613 controls). Analyses used mixed linear models as implemented in the OSCA software. We found a significantly higher proportion of neutrophils in cases compared to controls which replicated in an independent cohort from the Netherlands (1159 cases and 637 controls). The OSCA MOMENT linear mixed model has been shown in simulations to best account for confounders. When combined in a methylation profile score, the 25 most-associated probes identified by MOMENT significantly classified case–control status in the Netherlands sample (area under the curve, AUC = 0.65, CI95% = [0.62–0.68], p = 8.3 × 10−22). The maximum AUC achieved was 0.69 (CI95% = [0.66–0.71], p = 4.3 × 10−34) when cell-type proportion was included in the predictor.

We conducted DNA methylation association analyses using Illumina 450K data from whole blood for an Australian amyotrophic lateral sclerosis (ALS) case-control cohort (782 cases and 613 controls). Analyses used mixed linear models as implemented in the OSCA software. We found a significantly higher proportion of neutrophils in cases compared to controls which replicated in an independent cohort from the Netherlands (1159 cases and 637 controls). The OSCA MOMENT linear mixed model has been shown in simulations to best account for confounders. When combined in a methylation profile score, the 25 most-associated probes identified by MOMENT significantly classified case-control status in the Netherlands sample (area under the curve, AUC = 0.65, CI95% = [0.62-0.68], p = 8.3 × 10-22). The maximum AUC achieved was 0.69 (CI95% = [0.66-0.71], p = 4.3 × 10-34) when cell-type proportion was included in the predictor.

Alzheimer's disease (AD) is associated with the intracellular aggregation of hyperphosphorylated tau and the accumulation of β-amyloid in the neocortex. We use transgenic mice harboring human tau (rTg4510) and amyloid precursor protein (J20) mutations to investigate transcriptional changes associated with the progression of tau and amyloid pathology. rTg4510 mice are characterized by widespread transcriptional differences in the entorhinal cortex with changes paralleling neuropathological burden across multiple brain regions. Differentially expressed transcripts overlap with genes identified in genetic studies of familial and sporadic AD. Systems-level analyses identify discrete co-expression networks associated with the progressive accumulation of tau that are enriched for genes and pathways previously implicated in AD pathology and overlap with co-expression networks identified in human AD cortex. Our data provide further evidence for an immune-response component in the accumulation of tau and reveal molecular pathways associated with the progression of AD neuropathology.

Abstract

Importance: Genetic factors increase risk of dementia, but the extent to which this can be offset by lifestyle factors is unknown.

Objective: to investigate whether a healthy lifestyle is associated with lower risk of dementia regardless of genetic risk.

Design, Setting, and Participants: a retrospective cohort study that included adults of European ancestry aged at least 60 years without cognitive impairment or dementia at baseline. Participants joined the UK Biobank study from 2006 to 2010 and were followed up until 2016 or 2017.

Exposures: a polygenic risk score for dementia with low (lowest quintile), intermediate (quintiles 2 to 4), and high (highest quintile) risk categories and a weighted healthy lifestyle score, including no current smoking, regular physical activity, healthy diet, and moderate alcohol consumption, categorized into favorable, intermediate, and unfavorable lifestyles.

Main Outcomes and Measures: Incident all-cause dementia, ascertained through hospital inpatient and death records.

Results: a total of 196 383 individuals (mean [SD] age, 64.1 [2.9] years; 52.7% were women) were followed up for 1 545 433 person-years (median [interquartile range] follow-up, 8.0 [7.4-8.6] years). Overall, 68.1% of participants followed a favorable lifestyle, 23.6% followed an intermediate lifestyle, and 8.2% followed an unfavorable lifestyle. Twenty percent had high polygenic risk scores, 60% had intermediate risk scores, and 20% had low risk scores. of the participants with high genetic risk, 1.23% (95% CI, 1.13%-1.35%) developed dementia compared with 0.63% (95% CI, 0.56%-0.71%) of the participants with low genetic risk (adjusted hazard ratio, 1.91 [95% CI, 1.64-2.23]). of the participants with a high genetic risk and unfavorable lifestyle, 1.78% (95% CI, 1.38%-2.28%) developed dementia compared with 0.56% (95% CI, 0.48%-0.66%) of participants with low genetic risk and favorable lifestyle (hazard ratio, 2.83 [95% CI, 2.09-3.83]). There was no significant interaction between genetic risk and lifestyle factors (P = .99). Among participants with high genetic risk, 1.13% (95% CI, 1.01%-1.26%) of those with a favorable lifestyle developed dementia compared with 1.78% (95% CI, 1.38%-2.28%) with an unfavorable lifestyle (hazard ratio, 0.68 [95% CI, 0.51-0.90]).

Conclusions and Relevance: Among older adults without cognitive impairment or dementia, both an unfavorable lifestyle and high genetic risk were significantly associated with higher dementia risk. A favorable lifestyle was associated with a lower dementia risk among participants with high genetic risk.

Common genetic variation contributes a substantial proportion of risk for both schizophrenia and bipolar disorder. Furthermore, there is evidence of significant, but not complete, overlap in genetic risk between the two disorders. It has been hypothesised that genetic variants conferring risk for these disorders do so by influencing brain development, leading to the later emergence of symptoms. The comparative profile of risk gene expression for schizophrenia and bipolar disorder across development over different brain regions however remains unclear. Using genotypes derived from genome-wide associations studies of the largest available cohorts of patients and control subjects, we investigated whether genes enriched for schizophrenia and bipolar disorder association show a bias for expression across any of 13 developmental stages in prefrontal cortical and subcortical brain regions. We show that genetic association with schizophrenia is positively correlated with expression in the prefrontal cortex during early midfetal development and early infancy, and negatively correlated with expression during late childhood, which stabilises in adolescence. In contrast, risk-associated genes for bipolar disorder did not exhibit a bias towards expression at any prenatal stage, although the pattern of postnatal expression was similar to that of schizophrenia. These results highlight the dynamic expression of genes harbouring risk for schizophrenia and bipolar disorder across prefrontal cortex development and support the hypothesis that prenatal neurodevelopmental events are more strongly associated with schizophrenia than bipolar disorder.

BACKGROUND: Previous studies have reported associations between attention-deficit/hyperactivity disorder symptoms and DNA methylation in children. We report the first epigenome-wide association study meta-analysis of adult attention-deficit/hyperactivity disorder symptoms, based on peripheral blood DNA methylation (Infinium HumanMethylation450K array) in three population-based adult cohorts. METHODS: an epigenome-wide association study was performed in the Netherlands Twin Register (N = 2258, mean age 37 years), Dunedin Multidisciplinary Health and Development Study (N = 800, age 38 years), and Environmental Risk Longitudinal Twin Study (N = 1631, age 18 years), and results were combined through meta-analysis (total sample size N = 4689). Region-based analyses accounting for the correlation between nearby methylation sites were also performed. RESULTS: One epigenome-wide significant differentially methylated position was detected in the Dunedin study, but meta-analysis did not detect differentially methylated positions that were robustly associated across cohorts. In region-based analyses, six significant differentially methylation regions (DMRs) were identified in the Netherlands Twin Register, 19 in the Dunedin study, and none in the Environmental Risk Longitudinal Twin Study. of these DMRs, 92% were associated with methylation quantitative trait loci, and 68% showed moderate to large blood-brain correlations for DNA methylation levels. DMRs included six nonoverlapping DMRs (three in the Netherlands Twin Register, three in the Dunedin study) in the major histocompatibility complex, which were associated with expression of genes in the major histocompatibility complex, including C4A and C4B, previously implicated in schizophrenia. CONCLUSIONS: Our findings point at new candidate loci involved in immune and neuronal functions that await further replication. Our work also illustrates the need for further research to examine to what extent epigenetic associations with psychiatric traits depend on characteristics such as age, comorbidities, exposures, and genetic background.

Large-scale epigenome-wide association meta-analyses have identified multiple 'signatures'' of smoking. Drawing on these findings, we describe the construction of a polyepigenetic DNA methylation score that indexes smoking behavior and that can be utilized for multiple purposes in population health research. To validate the score, we use data from two birth cohort studies: the Dunedin Longitudinal Study, followed to age-38 years, and the Environmental Risk Study, followed to age-18 years. Longitudinal data show that changes in DNA methylation accumulate with increased exposure to tobacco smoking and attenuate with quitting. Data from twins discordant for smoking behavior show that smoking influences DNA methylation independently of genetic and environmental risk factors. Physiological data show that changes in DNA methylation track smoking-related changes in lung function and gum health over time. Moreover, DNA methylation changes predict corresponding changes in gene expression in pathways related to inflammation, immune response, and cellular trafficking. Finally, we present prospective data about the link between adverse childhood experiences (ACEs) and epigenetic modifications; these findings document the importance of controlling for smoking-related DNA methylation changes when studying biological embedding of stress in life-course research. We introduce the polyepigenetic DNA methylation score as a tool both for discovery and theory-guided research in epigenetic epidemiology.

Most variants associated with complex phenotypes in genome-wide association studies (GWAS) do not directly index coding changes affecting protein structure. Instead they are hypothesized to influence gene regulation, with common variants associated with disease being enriched in regulatory domains including enhancers and regions of open chromatin. There is interest, therefore, in using epigenomic annotation data to identify the specific regulatory mechanisms involved and prioritize risk variants. We quantified lysine H3K27 acetylation (H3K27ac) - a robust mark of active enhancers and promoters that is strongly correlated with gene expression and transcription factor binding - across the genome in entorhinal cortex samples using chromatin immunoprecipitation followed by highly parallel sequencing (ChIP-seq). H3K27ac peaks were called using high quality reads combined across all samples and formed the basis of partitioned heritability analysis using LD score regression along with publicly-available GWAS results for seven psychiatric and neurodegenerative traits. Heritability for all seven brain traits was significantly enriched in these H3K27ac peaks (enrichment ranging from 1.09-2.13) compared to regions of the genome containing other active regulatory and functional elements across multiple cell types and tissues. The strongest enrichments were for amyotrophic lateral sclerosis (ALS) (enrichment = 2.19; 95% CI = 2.12-2.27), autism (enrichment = 2.11; 95% CI = 2.05-2.16) and major depressive disorder (enrichment = 2.04; 95% CI = 1.92-2.16). Much lower enrichments were observed for 14 non-brain disorders, although we identified enrichment in cortical H3K27ac domains for body mass index (enrichment = 1.16; 95% CI = 1.13-1.19), ever smoked (enrichment = 2.07; 95% CI = 2.04-2.10), HDL (enrichment = 1.53; 95% CI = 1.45-1.62) and trigylcerides (enrichment = 1.33; 95% CI = 1.24-1.42). These results indicate that risk alleles for brain disorders are preferentially located in regions of regulatory/enhancer function in the cortex, further supporting the hypothesis that genetic variants for these phenotypes influence gene regulation in the brain.

Autism spectrum disorder (ASD) encompasses a collection of complex neuropsychiatric disorders characterized by deficits in social functioning, communication and repetitive behaviour. Building on recent studies supporting a role for developmentally moderated regulatory genomic variation in the molecular aetiology of ASD, we quantified genome-wide patterns of DNA methylation in 223 post-mortem tissues samples isolated from three brain regions [prefrontal cortex, temporal cortex and cerebellum (CB)] dissected from 43 ASD patients and 38 non-psychiatric control donors. We identified widespread differences in DNA methylation associated with idiopathic ASD (iASD), with consistent signals in both cortical regions that were distinct to those observed in the CB. Individuals carrying a duplication on chromosome 15q (dup15q), representing a genetically defined subtype of ASD, were characterized by striking differences in DNA methylationacross a discrete domain spanning an imprinted gene cluster within the duplicated region. In addition to the dramatic cis-effects on DNA methylation observed in dup15q carriers, we identified convergent methylomic signatures associated with both iASD and dup15q, reflecting the findings from previous studies of gene expression and H3K27ac. Cortical co-methylation network analysis identified a number of co-methylated modules significantly associated with ASD that are enriched for genomic regions annotated to genes involved in the immune system, synaptic signalling and neuronal regulation. Our study represents the first systematic analysis of DNA methylation associated with ASD across multiple brain regions, providing novel evidence for convergent molecular signatures associated with both idiopathic and syndromic autism.

BACKGROUND: There has been a steady increase in the number of studies aiming to identify DNA methylation differences associated with complex phenotypes. Many of the challenges of epigenetic epidemiology regarding study design and interpretation have been discussed in detail, however there are analytical concerns that are outstanding and require further exploration. In this study we seek to address three analytical issues. First, we quantify the multiple testing burden and propose a standard statistical significance threshold for identifying DNA methylation sites that are associated with an outcome. Second, we establish whether linear regression, the chosen statistical tool for the majority of studies, is appropriate and whether it is biased by the underlying distribution of DNA methylation data. Finally, we assess the sample size required for adequately powered DNA methylation association studies. RESULTS: We quantified DNA methylation in the Understanding Society cohort (n = 1175), a large population based study, using the Illumina EPIC array to assess the statistical properties of DNA methylation association analyses. By simulating null DNA methylation studies, we generated the distribution of p-values expected by chance and calculated the 5% family-wise error for EPIC array studies to be 9 × 10- 8. Next, we tested whether the assumptions of linear regression are violated by DNA methylation data and found that the majority of sites do not satisfy the assumption of normal residuals. Nevertheless, we found no evidence that this bias influences analyses by increasing the likelihood of affected sites to be false positives. Finally, we performed power calculations for EPIC based DNA methylation studies, demonstrating that existing studies with data on ~ 1000 samples are adequately powered to detect small differences at the majority of sites. CONCLUSION: We propose that a significance threshold of P 

Childhood psychotic experiences (PEs), such as seeing or hearing things that others do not, or extreme paranoia, are relatively common with around 1 in 20 children reporting them at age 12. Childhood PEs are often distressing and can be predictive of schizophrenia, other psychiatric disorders, and suicide attempts in adulthood, particularly if they persist during adolescence. Previous research has demonstrated that methylomic signatures in blood could be potential biomarkers of psychotic phenomena. This study explores the association between DNA methylation (DNAm) and the emergence, persistence, and remission of PEs in childhood and adolescence. DNAm profiles were obtained from the ALSPAC cohort at birth, age 7, and age 15/17 (n = 901). PEs were assessed through interviews with participants at ages 12 and 18. We identified PE-associated probes (p 

Schizophrenia (SCZ) is associated with high mortality. DNA methylation levels vary over the life course, and pre-selected combinations of methylation array probes can be used to estimate "methylation age" (mAge). mAge correlates highly with chronological age but when it differs, termed mAge acceleration, it has been previously associated with all-cause mortality. We tested the association between mAge acceleration and mortality in SCZ and controls. We selected 190 SCZ cases and 190 controls from the Sweden Schizophrenia Study. Cases were identified from the Swedish Hospital Discharge Register with ≥5 specialist treatment contacts and ≥5 antipsychotic prescriptions. Controls had no psychotic disorder or antipsychotics. Subjects were selected if they had died or survived during follow-up (2:1 oversampling). Extracted DNA was assayed on the Illumina MethylationEPIC array. mAge was regressed on age at sampling to obtain mAge acceleration. Using Cox proportional hazards regression, the association between mAge acceleration and mortality was tested. After quality control, the following were available: n = 126 SCZ died, 63 SCZ alive, 127 controls died, 62 controls alive. In the primary analyses, we did not find a significant association between mAge acceleration and SCZ mortality (adjusted p > 0.005). Sensitivity analyses excluding SCZ cases with pre-existing cancer demonstrated a significant association between the Hannum mAge acceleration and mortality (hazard ratio = 1.13, 95% confidence interval = 1.04-1.22, p = 0.005). Per our pre-specified criteria, we did not confirm our primary hypothesis that mAge acceleration would predict subsequent mortality in people with SCZ, but we cannot rule out smaller effects or effects in patient subsets.

BACKGROUND: Alzheimer's disease is a progressive neurodegenerative disorder that is hypothesized to involve epigenetic dysfunction. Previous studies of DNA modifications in Alzheimer's disease have been unable to distinguish between DNA methylation and DNA hydroxymethylation. DNA hydroxymethylation has been shown to be enriched in the human brain, although its role in Alzheimer's disease has not yet been fully explored. Here, we utilize oxidative bisulfite conversion, in conjunction with the Illumina Infinium Human Methylation 450K microarray, to identify neuropathology-associated differential DNA methylation and DNA hydroxymethylation in the entorhinal cortex. RESULTS: We identified one experiment-wide significant differentially methylated position residing in the WNT5B gene. Next, we investigated pathology-associated regions consisting of multiple adjacent loci. We identified one significant differentially hydroxymethylated region consisting of four probes spanning 104 bases in the FBXL16 gene. We also identified two significant differentially methylated regions: one consisting of two probes in a 93 base-pair region in the ANK1 gene and the other consisting of six probes in a 99-base pair region in the ARID5B gene. We also highlighted three regions that show alterations in unmodified cytosine: two probes in a 39-base pair region of ALLC, two probes in a 69-base pair region in JAG2, and the same six probes in ARID5B that were differentially methylated. Finally, we replicated significant ANK1 disease-associated hypermethylation and hypohydroxymethylation patterns across eight CpG sites in an extended 118-base pair region in an independent cohort using oxidative-bisulfite pyrosequencing. CONCLUSIONS: Our study represents the first epigenome-wide association study of both DNA methylation and hydroxymethylation in Alzheimer's disease entorhinal cortex. We demonstrate that previous estimates of DNA hypermethylation in ANK1 in Alzheimer's disease were underestimates as it is confounded by hypohydroxymethylation.

BACKGROUND: the Horvath epigenetic clock is widely used. It predicts age quite well from 353 CpG sites in the DNA methylation profile in unknown samples and has been used to calculate "age acceleration" in various tissues and environments. RESULTS: the model systematically underestimates age in tissues from older people. This is seen in all examined tissues but most strongly in the cerebellum and is consistently observed in multiple datasets. Age acceleration is thus age-dependent, and this can lead to spurious associations. The current literature includes examples of association tests with age acceleration calculated in a wide variety of ways. CONCLUSIONS: the concept of an epigenetic clock is compelling, but caution should be taken in interpreting associations with age acceleration. Association tests of age acceleration should include age as a covariate.

X-chromosome inactivation (XCI), i.e. the inactivation of one of the female X chromosomes, restores equal expression of X-chromosomal genes between females and males. However, ~10% of genes show variable degrees of escape from XCI between females, although little is known about the causes of variable XCI. Using a discovery data-set of 1867 females and 1398 males and a replication sample of 3351 females, we show that genetic variation at three autosomal loci is associated with female-specific changes in X-chromosome methylation. Through cis-eQTL expression analysis, we map these loci to the genes SMCHD1/METTL4, TRIM6/HBG2, and ZSCAN9. Low-expression alleles of the loci are predominantly associated with mild hypomethylation of CpG islands near genes known to variably escape XCI, implicating the autosomal genes in variable XCI. Together, these results suggest a genetic basis for variable escape from XCI and highlight the potential of a population genomics approach to identify genes involved in XCI.

MotivationThe datasets generated by DNA methylation analyses are getting bigger. With the release of the HumanMethylationEPIC micro-array and datasets containing thousands of samples, analyses of these large datasets using R are becoming impractical due to large memory requirements. As a result there is an increasing need for computationally efficient methodologies to perform meaningful analysis on high dimensional data.ResultsHere we introduce the bigmelon R package, which provides a memory efficient workflow that enables users to perform the complex, large scale analyses required in epigenome wide association studies (EWAS) without the need for large RAM. Building on top of the CoreArray Genomic Data Structure file format and libraries packaged in the gdsfmt package, we provide a practical workflow that facilitates the reading-in, preprocessing, quality control and statistical analysis of DNA methylation data.We demonstrate the capabilities of the bigmelon package using a large dataset consisting of 1193 human blood samples from the Understanding Society: UK Household Longitudinal Study, assayed on the EPIC micro-array platform.Availability and implementationThe bigmelon package is available on Bioconductor (http://bioconductor.org/packages/bigmelon/). The Understanding Society dataset is available at https://www.understandingsociety.ac.uk/about/health/data upon request.Supplementary informationSupplementary data are available at Bioinformatics online.

INTRODUCTION: Alzheimer's disease is a neurodegenerative disorder that is hypothesized to involve epigenetic dysregulation of gene expression in the brain. METHODS: We performed an epigenome-wide association study to identify differential DNA methylation associated with neuropathology in prefrontal cortex and superior temporal gyrus samples from 147 individuals, replicating our findings in two independent data sets (N = 117 and 740). RESULTS: We identify elevated DNA methylation associated with neuropathology across a 48-kb region spanning 208 CpG sites within the HOXA gene cluster. A meta-analysis of the top-ranked probe within the HOXA3 gene (cg22962123) highlighted significant hypermethylation across all three cohorts (P = 3.11 × 10-18). DISCUSSION: We present robust evidence for elevated DNA methylation associated with Alzheimer's disease neuropathology spanning the HOXA gene cluster on chromosome 7. These data add to the growing evidence highlighting a role for epigenetic variation in Alzheimer's disease, implicating the HOX gene family as a target for future investigation.

Autism spectrum disorder (ASD) is a severe neurodevelopmental disorder characterized by deficits in social communication and restricted, repetitive behaviors, interests, or activities. The etiology of ASD involves both inherited and environmental risk factors, with epigenetic processes hypothesized as one mechanism by which both genetic and non-genetic variation influence gene regulation and pathogenesis. The aim of this study was to identify DNA methylation biomarkers of ASD detectable at birth.

BACKGROUND: Genetic influences on gene expression in the human fetal brain plausibly impact upon a variety of postnatal brain-related traits, including susceptibility to neuropsychiatric disorders. However, to date, there have been no studies that have mapped genome-wide expression quantitative trait loci (eQTL) specifically in the human prenatal brain. RESULTS: We performed deep RNA sequencing and genome-wide genotyping on a unique collection of 120 human brains from the second trimester of gestation to provide the first eQTL dataset derived exclusively from the human fetal brain. We identify high confidence cis-acting eQTL at the individual transcript as well as whole gene level, including many mapping to a common inversion polymorphism on chromosome 17q21. Fetal brain eQTL are enriched among risk variants for postnatal conditions including attention deficit hyperactivity disorder, schizophrenia, and bipolar disorder. We further identify changes in gene expression within the prenatal brain that potentially mediate risk for neuropsychiatric traits, including increased expression of C4A in association with genetic risk for schizophrenia, increased expression of LRRC57 in association with genetic risk for bipolar disorder, and altered expression of multiple genes within the chromosome 17q21 inversion in association with variants influencing the personality trait of neuroticism. CONCLUSIONS: We have mapped eQTL operating in the human fetal brain, providing evidence that these confer risk to certain neuropsychiatric disorders, and identifying gene expression changes that potentially mediate susceptibility to these conditions.

Major depressive disorder (MDD) is a common illness accompanied by considerable morbidity, mortality, costs, and heightened risk of suicide. We conducted a genome-wide association meta-analysis based in 135,458 cases and 344,901 controls and identified 44 independent and significant loci. The genetic findings were associated with clinical features of major depression and implicated brain regions exhibiting anatomical differences in cases. Targets of antidepressant medications and genes involved in gene splicing were enriched for smaller association signal. We found important relationships of genetic risk for major depression with educational attainment, body mass, and schizophrenia: lower educational attainment and higher body mass were putatively causal, whereas major depression and schizophrenia reflected a partly shared biological etiology. All humans carry lesser or greater numbers of genetic risk factors for major depression. These findings help refine the basis of major depression and imply that a continuous measure of risk underlies the clinical phenotype.

Heterozygous mutation of the transcription factor HNF1B is the most common cause of monogenetic developmental renal disease. Disease-associated mutations fall into two categories: HNF1B intragenic mutations and a 1.3 Mb deletion at chromosome 17q12. An increase in neurodevelopmental disorders has been observed in individuals harbouring the 17q12 deletion but not in patients with HNF1B coding mutations.Previous investigations have concentrated on identifying a genetic cause for the increase in behavioural problems seen in 17q12 deletion carriers. We have taken the alternative approach of investigating the DNA methylation profile of these two HNF1B genotype groups along with controls matched for age, gender and diabetes status using the Illumina 450K DNA methylation array (total sample n = 60).We identified a number of differentially methylated probes (DMPs) that were associated with HNF1B-associated disease and passed our stringent experiment-wide significance threshold. These associations were largely driven by the deletion patients and the majority of the significant probes mapped to the 17q12 deletion locus. The observed changes in DNA methylation at this locus were not randomly dispersed and occurred in clusters, suggesting a regulatory mechanism reacting to haploinsufficiency across the entire deleted region.Along with these deletion-specific changes in DNA methylation, we also identified a shared DNA methylation signature in both mutation and deletion patient groups indicating that haploinsufficiency of HNF1B impacts on the methylome of a number of genes, giving further insight to the role of HNF1B.

Characterizing the complex relationship between genetic, epigenetic, and transcriptomic variation has the potential to increase understanding about the mechanisms underpinning health and disease phenotypes. We undertook a comprehensive analysis of common genetic variation on DNA methylation (DNAm) by using the Illumina EPIC array to profile samples from the UK Household Longitudinal study. We identified 12,689,548 significant DNA methylation quantitative trait loci (mQTL) associations (p < 6.52 × 10−14) occurring between 2,907,234 genetic variants and 93,268 DNAm sites, including a large number not identified by previous DNAm-profiling methods. We demonstrate the utility of these data for interpreting the functional consequences of common genetic variation associated with > 60 human traits by using summary-data-based Mendelian randomization (SMR) to identify 1,662 pleiotropic associations between 36 complex traits and 1,246 DNAm sites. We also use SMR to characterize the relationship between DNAm and gene expression and thereby identify 6,798 pleiotropic associations between 5,420 DNAm sites and the transcription of 1,702 genes. Our mQTL database and SMR results are available via a searchable online database as a resource to the research community.

Cognitive functions are important correlates of health outcomes across the life-course. Individual differences in cognitive functions are partly heritable. Epigenetic modifications, such as DNA methylation, are susceptible to both genetic and environmental factors and may provide insights into individual differences in cognitive functions. Epigenome-wide meta-analyses for blood-based DNA methylation levels at ~420,000 CpG sites were performed for seven measures of cognitive functioning using data from 11 cohorts. CpGs that passed a Bonferroni correction, adjusting for the number of CpGs and cognitive tests, were assessed for: longitudinal change; being under genetic control (methylation QTLs); and associations with brain health (structural MRI), brain methylation and Alzheimer's disease pathology. Across the seven measures of cognitive functioning (meta-analysis n range: 2557-6809), there were epigenome-wide significant (P 

Accelerated DNA methylation age is linked to all-cause mortality and environmental factors, but studies of associations with socioeconomic position are limited. Researchers generally use small selected samples, and it is unclear how findings obtained with 2 commonly used methods for calculating methylation age (the Horvath method and the Hannum method) translate to general population samples including younger and older adults. Among 1,099 United Kingdom adults aged 28-98 years in 2011-2012, we assessed the relationship of Horvath and Hannum DNA methylation age acceleration with a range of social position measures: current income and employment, education, income and unemployment across a 12-year period, and childhood social class. Accounting for confounders, participants who had been less advantaged in childhood were epigenetically "older" as adults: in comparison with participants who had professional/managerial parents, Hannum age was 1.07 years higher (95% confidence interval: 0.20, 1.94) for participants with parents in semiskilled/unskilled occupations and 1.85 years higher (95% confidence interval: 0.67, 3.02) for those without a working parent at age 14 years. No other robust associations were seen. Results accord with research implicating early life circumstances as critical for DNA methylation age in adulthood. Since methylation age acceleration as measured by the Horvath and Hannum estimators appears strongly linked to chronological age, researchers examining associations with the social environment must take steps to avoid age-related confounding.

BACKGROUND: Numerous risk factors for dementia are well established, though the causal nature of these associations remains unclear. OBJECTIVE: to systematically review Mendelian randomization (MR) studies investigating causal relationships between risk factors and global cognitive function or dementia. METHODS: We searched five databases from inception to February 2017 and conducted citation searches including MR studies investigating the association between any risk factor and global cognitive function, all-cause dementia or dementia subtypes. Two reviewers independently assessed titles and abstracts, full-texts, and study quality. RESULTS: We included 18 MR studies investigating education, lifestyle factors, cardiovascular factors and related biomarkers, diabetes related and other endocrine factors, and telomere length. Studies were of predominantly good quality, however eight received low ratings for sample size and statistical power. The most convincing causal evidence was found for an association of shorter telomeres with increased risk of Alzheimer's disease (AD). Causal evidence was weaker for smoking quantity, vitamin D, homocysteine, systolic blood pressure, fasting glucose, insulin sensitivity, and high-density lipoprotein cholesterol. Well-replicated associations were not present for most exposures and we cannot fully discount survival and diagnostic bias, or the potential for pleiotropic effects. CONCLUSIONS: Genetic evidence supported a causal association between telomere length and AD, whereas limited evidence for other risk factors was largely inconclusive with tentative evidence for smoking quantity, vitamin D, homocysteine, and selected metabolic markers. The lack of stronger evidence for other risk factors may reflect insufficient statistical power. Larger well-designed MR studies would therefore help establish the causal status of these dementia risk factors.

Identifying genes regulating the pace of epigenetic ageing represents a new frontier in genome-wide association studies (GWASs). Here using 1,796 brain samples from 1,163 individuals, we carry out a GWAS of two DNA methylation-based biomarkers of brain age: the epigenetic ageing rate and estimated proportion of neurons. Locus 17q11.2 is significantly associated (P=4.5 × 10-9) with the ageing rate across five brain regions and harbours a cis-expression quantitative trait locus for EFCAB5 (P=3.4 × 10-20). Locus 1p36.12 is significantly associated (P=2.2 × 10-8) with epigenetic ageing of the prefrontal cortex, independent of the proportion of neurons. Our GWAS of the proportion of neurons identified two genome-wide significant loci (10q26 and 12p13.31) and resulted in a gene set that overlaps significantly with sets found by GWAS of age-related macular degeneration (P=1.4 × 10-12), ulcerative colitis (P

OBJECTIVE: Genome-wide association studies (GWAS) have been successful at identifying associations with stroke and stroke subtypes, but have not yet identified any associations solely with small vessel stroke (SVS). SVS comprises one quarter of all ischemic stroke and is a major manifestation of cerebral small vessel disease, the primary cause of vascular cognitive impairment. Studies across neurological traits have shown that younger-onset cases have an increased genetic burden. We leveraged this increased genetic burden by performing an age-at-onset informed GWAS meta-analysis, including a large younger-onset SVS population, to identify novel associations with stroke. METHODS: We used a three-stage age-at-onset informed GWAS to identify novel genetic variants associated with stroke. On identifying a novel locus associated with SVS, we assessed its influence on other small vessel disease phenotypes, as well as on messenger RNA (mRNA) expression of nearby genes, and on DNA methylation of nearby CpG sites in whole blood and in the fetal brain. RESULTS: We identified an association with SVS in 4,203 cases and 50,728 controls on chromosome 16q24.2 (odds ratio [OR; 95% confidence interval {CI}] = 1.16 [1.10-1.22]; p = 3.2 × 10-9 ). The lead single-nucleotide polymorphism (rs12445022) was also associated with cerebral white matter hyperintensities (OR [95% CI] = 1.10 [1.05-1.16]; p = 5.3 × 10-5 ; N = 3,670), but not intracerebral hemorrhage (OR [95% CI] = 0.97 [0.84-1.12]; p = 0.71; 1,545 cases, 1,481 controls). rs12445022 is associated with mRNA expression of ZCCHC14 in arterial tissues (p = 9.4 × 10-7 ) and DNA methylation at probe cg16596957 in whole blood (p = 5.3 × 10-6 ). INTERPRETATION: 16q24.2 is associated with SVS. Associations of the locus with expression of ZCCHC14 and DNA methylation suggest the locus acts through changes to regulatory elements. Ann Neurol 2017;81:383-394.

Major depressive disorder (MDD) represents a major social and economic health issue and constitutes a major risk factor for suicide. The molecular pathology of suicidal depression remains poorly understood, although it has been hypothesised that regulatory genomic processes are involved in the pathology of both MDD and suicidality. In this study, genome-wide patterns of DNA methylation were assessed in depressed suicide completers (n=20) and compared with non-psychiatric, sudden-death controls (n=20) using tissue from two cortical brain regions (Brodmann Area 11 (BA11) and Brodmann Area 25 (BA25)). Analyses focused on identifying differentially methylated regions (DMRs) associated with suicidal depression and epigenetic variation were explored in the context of polygenic risk scores for major depression and suicide. Weighted gene co-methylation network analysis was used to identify modules of co-methylated loci associated with depressed suicide completers and polygenic burden for MDD and suicide attempt. We identified a DMR upstream of the PSORS1C3 gene, subsequently validated using bisulfite pyrosequencing and replicated in a second set of suicide samples, which is characterised by significant hypomethylation in both cortical brain regions in MDD suicide cases. We also identified discrete modules of co-methylated loci associated with polygenic risk burden for suicide attempt, but not major depression. Suicide-associated co-methylation modules were enriched among gene networks implicating biological processes relevant to depression and suicidality, including nervous system development and mitochondria function. Our data suggest that there are coordinated changes in DNA methylation associated with suicide that may offer novel insights into the molecular pathology associated with depressed suicide completers.

Most genetic variants identified in genome-wide association studies (GWASs) of complex traits are thought to act by affecting gene regulation rather than directly altering the protein product. As a consequence, the actual genes involved in disease are not necessarily the most proximal to the associated variants. By integrating data from GWAS analyses with those from genetic studies of regulatory variation, it is possible to identify variants pleiotropically associated with both a complex trait and measures of gene regulation. In this study, we used summary-data-based Mendelian randomization (SMR), a method developed to identify variants pleiotropically associated with both complex traits and gene expression, to identify variants associated with complex traits and DNA methylation. We used large DNA methylation quantitative trait locus (mQTL) datasets generated from two different tissues (blood and fetal brain) to prioritize genes for >40 complex traits with robust GWAS data and found considerable overlap with the results of SMR analyses performed with expression QTL (eQTL) data. We identified multiple examples of variable DNA methylation associated with GWAS variants for a range of complex traits, demonstrating the utility of this approach for refining genetic association signals.

BACKGROUND: DNA methylation is an important epigenetic mechanism involved in gene regulation, with alterations in DNA methylation in the nuclear genome being linked to numerous complex diseases. Mitochondrial DNA methylation is a phenomenon that is receiving ever-increasing interest, particularly in diseases characterized by mitochondrial dysfunction; however, most studies have been limited to the investigation of specific target regions. Analyses spanning the entire mitochondrial genome have been limited, potentially due to the amount of input DNA required. Further, mitochondrial genetic studies have been previously confounded by nuclear-mitochondrial pseudogenes. Methylated DNA Immunoprecipitation Sequencing is a technique widely used to profile DNA methylation across the nuclear genome; however, reads mapped to mitochondrial DNA are often discarded. Here, we have developed an approach to control for nuclear-mitochondrial pseudogenes within Methylated DNA Immunoprecipitation Sequencing data. We highlight the utility of this approach in identifying differences in mitochondrial DNA methylation across regions of the human brain and pre-mortem blood. RESULTS: We were able to correlate mitochondrial DNA methylation patterns between the cortex, cerebellum and blood. We identified 74 nominally significant differentially methylated regions (p 

Genetic association studies provide evidence for a substantial polygenic component to schizophrenia, although the neurobiological mechanisms underlying the disorder remain largely undefined. Building on recent studies supporting a role for developmentally regulated epigenetic variation in the molecular aetiology of schizophrenia, this study aimed to identify epigenetic variation associated with both a diagnosis of schizophrenia and elevated polygenic risk burden for the disease across multiple brain regions. Genome-wide DNA methylation was quantified in 262 post-mortem brain samples, representing tissue from four brain regions (prefrontal cortex, striatum, hippocampus and cerebellum) from 41 schizophrenia patients and 47 controls. We identified multiple disease-associated and polygenic risk score-associated differentially methylated positions and regions, which are not enriched in genomic regions identified in genetic studies of schizophrenia and do not reflect direct genetic effects on DNA methylation. Our study represents the first analysis of epigenetic variation associated with schizophrenia across multiple brain regions and highlights the utility of polygenic risk scores for identifying molecular pathways associated with aetiological variation in complex disease.

Epigenetic disruption has been implicated in many diseases of aging, and age-associated DNA methylation changes at specific genomic loci in humans are strongly correlated with chronological age. The aim of this study was to explore the specificity of selected age-associated differentially methylated positions (aDMPs) identified in human epidemiological studies by quantifying DNA methylation across multiple tissues in homologous regions of the murine genome. We selected four high-confidence aDMPs (located in the vicinity of the ELOVL2, GLRA1, MYOD1 and PDE4C genes) and quantified DNA methylation across these regions in four tissues (blood, lung, cerebellum and hippocampus) from male and female C57BL/6J mice, ranging in age from fetal (embryonic day 17) to 630 days. We observed tissue-specific age-associated changes in DNA methylation that was directionally consistent with those observed in humans. These findings lend further support to the notion that changes in DNA methylation are associated with chronological age and suggest that these processes are often conserved across tissues and between mammalian species. Our data highlight the relevance of utilizing model systems, in which environmental and genetic influences can be carefully controlled, for the further study of these phenomena.

BACKGROUND: Schizophrenia is a highly heritable, neuropsychiatric disorder characterized by episodic psychosis and altered cognitive function. Despite success in identifying genetic variants associated with schizophrenia, there remains uncertainty about the causal genes involved in disease pathogenesis and how their function is regulated. RESULTS: We performed a multi-stage epigenome-wide association study, quantifying genome-wide patterns of DNA methylation in a total of 1714 individuals from three independent sample cohorts. We have identified multiple differentially methylated positions and regions consistently associated with schizophrenia across the three cohorts; these effects are independent of important confounders such as smoking. We also show that epigenetic variation at multiple loci across the genome contributes to the polygenic nature of schizophrenia. Finally, we show how DNA methylation quantitative trait loci in combination with Bayesian co-localization analyses can be used to annotate extended genomic regions nominated by studies of schizophrenia, and to identify potential regulatory variation causally involved in disease. CONCLUSIONS: This study represents the first systematic integrated analysis of genetic and epigenetic variation in schizophrenia, introducing a methodological approach that can be used to inform epigenome-wide association study analyses of other complex traits and diseases. We demonstrate the utility of using a polygenic risk score to identify molecular variation associated with etiological variation, and of using DNA methylation quantitative trait loci to refine the functional and regulatory variation associated with schizophrenia risk variants. Finally, we present strong evidence for the co-localization of genetic associations for schizophrenia and differential DNA methylation.

BACKGROUND: Single variant approaches have been successful in identifying DNA methylation quantitative trait loci (mQTL), although as with complex traits they lack the statistical power to identify the effects from rare genetic variants. We have undertaken extensive analyses to identify regions of low frequency and rare variants that are associated with DNA methylation levels. METHODS: We used repeated measurements of DNA methylation from five different life stages in human blood, taken from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. Variants were collapsed across CpG islands and their flanking regions to identify variants collectively associated with methylation, where no single variant was individually responsible for the observed signal. All analyses were undertaken using the sequence kernel association test. RESULTS: for loci where no individual variant mQTL was observed based on a single variant analysis, we identified 95 unique regions where the combined effect of low frequency variants (MAF ≤ 5%) provided strong evidence of association with methylation. For loci where there was previous evidence of an individual variant mQTL, a further 3 regions provided evidence of association between multiple low frequency variants and methylation levels. Effects were observed consistently across 5 different time points in the lifecourse and evidence of replication in the TwinsUK and Exeter cohorts was also identified. CONCLUSION: We have demonstrated the potential of this novel approach to mQTL analysis by analysing the combined effect of multiple low frequency or rare variants. Future studies should benefit from applying this approach as a complementary follow up to single variant analyses.

DNA methylation (DNAm) levels lend themselves for defining an epigenetic biomarker of aging known as the 'epigenetic clock'. Our genome-wide association study (GWAS) of cerebellar epigenetic age acceleration identifies five significant (P

Although mutations within the TREM2 gene have been robustly associated with Alzheimer's disease, it is not known whether alterations in the regulation of this gene are also involved in pathogenesis. Here, we present data demonstrating increased DNA methylation in the superior temporal gyrus in Alzheimer's disease brain at a CpG site located 289 bp upstream of the transcription start site of the TREM2 gene in 3 independent study cohorts using 2 different technologies (Illumina Infinium 450K methylation beadchip and pyrosequencing). A meta-analysis across all 3 cohorts reveals consistent AD-associated hypermethylation (p = 3.47E-08). This study highlights that extending genetic studies of TREM2 in AD to investigate epigenetic changes may nominate additional mechanisms by which disruption to this gene increases risk.

We characterized DNA methylation quantitative trait loci (mQTLs) in a large collection (n = 166) of human fetal brain samples spanning 56-166 d post-conception, identifying >16,000 fetal brain mQTLs. Fetal brain mQTLs were primarily cis-acting, enriched in regulatory chromatin domains and transcription factor binding sites, and showed substantial overlap with genetic variants that were also associated with gene expression in the brain. Using tissue from three distinct regions of the adult brain (prefrontal cortex, striatum and cerebellum), we found that most fetal brain mQTLs were developmentally stable, although a subset was characterized by fetal-specific effects. Fetal brain mQTLs were enriched amongst risk loci identified in a recent large-scale genome-wide association study (GWAS) of schizophrenia, a severe psychiatric disorder with a hypothesized neurodevelopmental component. Finally, we found that mQTLs can be used to refine GWAS loci through the identification of discrete sites of variable fetal brain methylation associated with schizophrenia risk variants.

BACKGROUND: the most widely utilized approaches for quantifying DNA methylation involve the treatment of genomic DNA with sodium bisulfite; however, this method cannot distinguish between 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC). Previous studies have shown that 5hmC is enriched in the brain, although little is known about its genomic distribution and how it differs between anatomical regions and individuals. In this study, we combine oxidative bisulfite (oxBS) treatment with the Illumina Infinium 450K BeadArray to quantify genome-wide patterns of 5hmC in two distinct anatomical regions of the brain from multiple individuals. RESULTS: We identify 37,145 and 65,563 sites passing our threshold for detectable 5hmC in the prefrontal cortex and cerebellum respectively, with 23,445 loci common across both brain regions. Distinct patterns of 5hmC are identified in each brain region, with notable differences in the genomic location of the most hydroxymethylated loci between these brain regions. Tissue-specific patterns of 5hmC are subsequently confirmed in an independent set of prefrontal cortex and cerebellum samples. CONCLUSIONS: This study represents the first systematic analysis of 5hmC in the human brain, identifying tissue-specific hydroxymethylated positions and genomic regions characterized by inter-individual variation in DNA hydroxymethylation. This study demonstrates the utility of combining oxBS-treatment with the Illumina 450k methylation array to systematically quantify 5hmC across the genome and the potential utility of this approach for epigenomic studies of brain disorders.

Given the tissue-specific nature of epigenetic processes, the assessment of disease-relevant tissue is an important consideration for epigenome-wide association studies (EWAS). Little is known about whether easily accessible tissues, such as whole blood, can be used to address questions about interindividual epigenomic variation in inaccessible tissues, such as the brain. We quantified DNA methylation in matched DNA samples isolated from whole blood and 4 brain regions (prefrontal cortex, entorhinal cortex, superior temporal gyrus, and cerebellum) from 122 individuals. We explored co-variation between tissues and the extent to which methylomic variation in blood is predictive of interindividual variation identified in the brain. For the majority of DNA methylation sites, interindividual variation in whole blood is not a strong predictor of interindividual variation in the brain, although the relationship with cortical regions is stronger than with the cerebellum. Variation at a subset of probes is strongly correlated across tissues, even in instances when the actual level of DNA methylation is significantly different between them. A substantial proportion of this co-variation, however, is likely to result from genetic influences. Our data suggest that for the majority of the genome, a blood-based EWAS for disorders where brain is presumed to be the primary tissue of interest will give limited information relating to underlying pathological processes. These results do not, however, discount the utility of using a blood-based EWAS to identify biomarkers of disease phenotypes manifest in the brain. We have generated a searchable database for the interpretation of data from blood-based EWAS analyses ( http://epigenetics.essex.ac.uk/bloodbrain/).

Childhood psychotic symptoms are associated with increased rates of schizophrenia, other psychiatric disorders, and suicide attempts in adulthood; thus, elucidating early risk indicators is crucial to target prevention efforts. There is considerable discordance for psychotic symptoms between monozygotic twins, indicating that child-specific non-genetic factors must be involved. Epigenetic processes may constitute one of these factors and have not yet been investigated in relation to childhood psychotic symptoms. Therefore, this study explored whether differences in DNA methylation at age 10 were associated with monozygotic twin discordance for psychotic symptoms at age 12. The Environmental Risk (E-Risk) Longitudinal Twin Study cohort of 2,232 children (1,116 twin pairs) was assessed for age-12 psychotic symptoms and 24 monozygotic twin pairs discordant for symptoms were identified for methylomic comparison. Children provided buccal samples at ages 5 and 10. DNA was bisulfite modified and DNA methylation was quantified using the Infinium HumanMethylation450 array. Differentially methylated positions (DMPs) associated with psychotic symptoms were subsequently tested in post-mortem prefrontal cortex tissue from adult schizophrenia patients and age-matched controls. Site-specific DNA methylation differences were observed at age 10 between monozygotic twins discordant for age-12 psychotic symptoms. Similar DMPs were not found at age 5. The top-ranked psychosis-associated DMP (cg23933044), located in the promoter of the C5ORF42 gene, was also hypomethylated in post-mortem prefrontal cortex brain tissue from schizophrenia patients compared to unaffected controls. These data tentatively suggest that epigenetic variation in peripheral tissue is associated with childhood psychotic symptoms and may indicate susceptibility to schizophrenia and other mental health problems.

BACKGROUND: Asthma is the most common chronic inflammatory disorder in children. The aetiology of asthma pathology is complex and highly heterogeneous, involving the interplay between genetic and environmental risk factors that is hypothesized to involve epigenetic processes. Our aim was to explore whether methylomic variation in early childhood is associated with discordance for asthma symptoms within monozygotic (MZ) twin pairs recruited from the Environmental Risk (E-Risk) longitudinal twin study. We also aimed to identify differences in DNA methylation that are associated with asthma that develops in childhood and persists into early adulthood as these may represent useful prognostic biomarkers. RESULTS: We examined genome-wide patterns of DNA methylation in buccal cell samples collected from 37 MZ twin pairs discordant for asthma at age 10. DNA methylation at individual CpG sites demonstrated significant variability within discordant MZ twin pairs with the top-ranked nominally significant differentially methylated position (DMP) located in the HGSNAT gene. We stratified our analysis by assessing DNA methylation differences in a sub-group of MZ twin pairs who remained persistently discordant for asthma at age 18. The top-ranked nominally significant DMP associated with persisting asthma is located in the vicinity of the HLX gene, which has been previously implicated in childhood asthma. CONCLUSIONS: We identified DNA methylation differences associated with childhood asthma in peripheral DNA samples from discordant MZ twin pairs. Our data suggest that differences in DNA methylation associated with childhood asthma which persists into early adulthood are distinct from those associated with asthma which remits.

Epigenetic processes play a key role in orchestrating transcriptional regulation during development. The importance of DNA methylation in fetal brain development is highlighted by the dynamic expression of de novo DNA methyltransferases during the perinatal period and neurodevelopmental deficits associated with mutations in the methyl-CpG binding protein 2 (MECP2) gene. However, our knowledge about the temporal changes to the epigenome during fetal brain development has, to date, been limited. We quantified genome-wide patterns of DNA methylation at ∼ 400,000 sites in 179 human fetal brain samples (100 male, 79 female) spanning 23 to 184 d post-conception. We identified highly significant changes in DNA methylation across fetal brain development at >7% of sites, with an enrichment of loci becoming hypomethylated with fetal age. Sites associated with developmental changes in DNA methylation during fetal brain development were significantly underrepresented in promoter regulatory regions but significantly overrepresented in regions flanking CpG islands (shores and shelves) and gene bodies. Highly significant differences in DNA methylation were observed between males and females at a number of autosomal sites, with a small number of regions showing sex-specific DNA methylation trajectories across brain development. Weighted gene comethylation network analysis (WGCNA) revealed discrete modules of comethylated loci associated with fetal age that are significantly enriched for genes involved in neurodevelopmental processes. This is, to our knowledge, the most extensive study of DNA methylation across human fetal brain development to date, confirming the prenatal period as a time of considerable epigenomic plasticity.

Inherited alleles account for most of the genetic risk for schizophrenia. However, new (de novo) mutations, in the form of large chromosomal copy number changes, occur in a small fraction of cases and disproportionally disrupt genes encoding postsynaptic proteins. Here we show that small de novo mutations, affecting one or a few nucleotides, are overrepresented among glutamatergic postsynaptic proteins comprising activity-regulated cytoskeleton-associated protein (ARC) and N-methyl-d-aspartate receptor (NMDAR) complexes. Mutations are additionally enriched in proteins that interact with these complexes to modulate synaptic strength, namely proteins regulating actin filament dynamics and those whose messenger RNAs are targets of fragile X mental retardation protein (FMRP). Genes affected by mutations in schizophrenia overlap those mutated in autism and intellectual disability, as do mutation-enriched synaptic pathways. Aligning our findings with a parallel case-control study, we demonstrate reproducible insights into aetiological mechanisms for schizophrenia and reveal pathophysiology shared with other neurodevelopmental disorders. © 2014 Macmillan Publishers Limited. All rights reserved.

Inherited alleles account for most of the genetic risk for schizophrenia. However, new (de novo) mutations, in the form of large chromosomal copy number changes, occur in a small fraction of cases and disproportionally disrupt genes encoding postsynaptic proteins. Here we show that small de novo mutations, affecting one or a few nucleotides, are overrepresented among glutamatergic postsynaptic proteins comprising activity-regulated cytoskeleton-associated protein (ARC) and N-methyl-d-aspartate receptor (NMDAR) complexes. Mutations are additionally enriched in proteins that interact with these complexes to modulate synaptic strength, namely proteins regulating actin filament dynamics and those whose messenger RNAs are targets of fragile X mental retardation protein (FMRP). Genes affected by mutations in schizophrenia overlap those mutated in autism and intellectual disability, as do mutation-enriched synaptic pathways. Aligning our findings with a parallel case-control study, we demonstrate reproducible insights into aetiological mechanisms for schizophrenia and reveal pathophysiology shared with other neurodevelopmental disorders.

Alzheimer's disease (AD) is a chronic neurodegenerative disorder that is characterized by progressive neuropathology and cognitive decline. We performed a cross-tissue analysis of methylomic variation in AD using samples from four independent human post-mortem brain cohorts. We identified a differentially methylated region in the ankyrin 1 (ANK1) gene that was associated with neuropathology in the entorhinal cortex, a primary site of AD manifestation. This region was confirmed as being substantially hypermethylated in two other cortical regions (superior temporal gyrus and prefrontal cortex), but not in the cerebellum, a region largely protected from neurodegeneration in AD, or whole blood obtained pre-mortem from the same individuals. Neuropathology-associated ANK1 hypermethylation was subsequently confirmed in cortical samples from three independent brain cohorts. This study represents, to the best of our knowledge, the first epigenome-wide association study of AD employing a sequential replication design across multiple tissues and highlights the power of this approach for identifying methylomic variation associated with complex disease.