A data-driven approach to preprocessing Illumina 450K methylation array data

BACKGROUND

As the most stable and experimentally accessible epigenetic mark, DNA methylation is of great interest to the research community. The landscape of DNA methylation across tissues, through development and in disease pathogenesis is not yet well characterized. Thus there is a need for rapid and cost effective methods for assessing genome-wide levels of DNA methylation. The Illumina Infinium HumanMethylation450 (450K) BeadChip is a very useful addition to the available methods for DNA methylation analysis but its complex design, incorporating two different assay methods, requires careful consideration. Accordingly, several normalization schemes have been published. We have taken advantage of known DNA methylation patterns associated with genomic imprinting and X-chromosome inactivation (XCI), in addition to the performance of SNP genotyping assays present on the array, to derive three independent metrics which we use to test alternative schemes of correction and normalization. These metrics also have potential utility as quality scores for datasets.

RESULTS

The standard index of DNA methylation at any specific CpG site is β = M/(M + U + 100) where M and U are methylated and unmethylated signal intensities, respectively. Betas (βs) calculated from raw signal intensities (the default GenomeStudio behavior) perform well, but using 11 methylomic datasets we demonstrate that quantile normalization methods produce marked improvement, even in highly consistent data, by all three metrics. The commonly used procedure of normalizing betas is inferior to the separate normalization of M and U, and it is also advantageous to normalize Type I and Type II assays separately. More elaborate manipulation of quantiles proves to be counterproductive.

CONCLUSIONS

Careful selection of preprocessing steps can minimize variance and thus improve statistical power, especially for the detection of the small absolute DNA methylation changes likely associated with complex disease phenotypes. For the convenience of the research community we have created a user-friendly R software package called wateRmelon, downloadable from bioConductor, compatible with the existing methylumi, minfi and IMA packages, that allows others to utilize the same normalization methods and data quality tests on 450K data.

Alzheimer's disease: early alterations in brain DNA methylation at ANK1, BIN1, RHBDF2 and other loci

We used a collection of 708 prospectively collected autopsied brains to assess the methylation state of the brain's DNA in relation to Alzheimer's disease (AD). We found that the level of methylation at 71 of the 415,848 interrogated CpGs was significantly associated with the burden of AD pathology, including CpGs in the ABCA7 and BIN1 regions, which harbor known AD susceptibility variants. We validated 11 of the differentially methylated regions in an independent set of 117 subjects. Furthermore, we functionally validated these CpG associations and identified the nearby genes whose RNA expression was altered in AD: ANK1, CDH23, DIP2A, RHBDF2, RPL13, SERPINF1 and SERPINF2. Our analyses suggest that these DNA methylation changes may have a role in the onset of AD given that we observed them in presymptomatic subjects and that six of the validated genes connect to a known AD susceptibility gene network.

Central and systemic endotoxin challenges exacerbate the local inflammatory response and increase neuronal death during chronic neurodegeneration

The contribution of inflammation to the progression of neurodegenerative diseases such as Alzheimer's, Parkinson's, and prion diseases is poorly understood. Brain inflammation in animal models of these diseases is dominated by chronic microglial activation with minimal proinflammatory cytokine expression. However, these inflammatory cells are "primed" to produce exaggerated inflammatory responses to subsequent lipopolysaccharide (LPS) challenges. We show that, using the ME7 model of prion disease, intracerebral challenge with LPS results in dramatic interleukin-1beta (IL-1beta) expression, neutrophil infiltration, and inducible nitric oxide synthase expression in the brain parenchyma of prion-diseased mice compared with the same challenge in normal mice. Systemic inflammation evoked by LPS also produced greater increases in proinflammatory cytokines, pentraxin 3, and inducible nitric oxide synthase transcription in prion-diseased mice than in control mice and induced microglial expression of IL-1beta. These systemic challenges also increased neuronal apoptosis in the brains of ME7 animals. Thus, both central and peripheral inflammation can exacerbate local brain inflammation and neuronal death. The finding that a single acute systemic inflammatory event can induce neuronal death in the CNS has implications for therapy in neurodegenerative diseases.

Functional annotation of the human brain methylome identifies tissue-specific epigenetic variation across brain and blood

Background

Dynamic changes to the epigenome play a critical role in establishing and maintaining cellular phenotype during differentiation, but little is known about the normal methylomic differences that occur between functionally distinct areas of the brain. We characterized intra- and inter-individual methylomic variation across whole blood and multiple regions of the brain from multiple donors.

Results

Distinct tissue-specific patterns of DNA methylation were identified, with a highly significant over-representation of tissue-specific differentially methylated regions (TS-DMRs) observed at intragenic CpG islands and low CG density promoters. A large proportion of TS-DMRs were located near genes that are differentially expressed across brain regions. TS-DMRs were significantly enriched near genes involved in functional pathways related to neurodevelopment and neuronal differentiation, including BDNF, BMP4, CACNA1A, CACA1AF, EOMES, NGFR, NUMBL, PCDH9, SLIT1, SLITRK1 and SHANK3. Although between-tissue variation in DNA methylation was found to greatly exceed between-individual differences within any one tissue, we found that some inter-individual variation was reflected across brain and blood, indicating that peripheral tissues may have some utility in epidemiological studies of complex neurobiological phenotypes.

Conclusions

This study reinforces the importance of DNA methylation in regulating cellular phenotype across tissues, and highlights genomic patterns of epigenetic variation across functionally distinct regions of the brain, providing a resource for the epigenetics and neuroscience research communities.

Methylomic profiling implicates cortical deregulation of ANK1 in Alzheimer's disease

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.

Interindividual methylomic variation across blood, cortex, and cerebellum: implications for epigenetic studies of neurological and neuropsychiatric phenotypes

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/).

Association of plasma clusterin concentration with severity, pathology, and progression in Alzheimer disease

CONTEXT

Blood-based analytes may be indicators of pathological processes in Alzheimer disease (AD).

OBJECTIVE

To identify plasma proteins associated with AD pathology using a combined proteomic and neuroimaging approach.

DESIGN

Discovery-phase proteomics to identify plasma proteins associated with correlates of AD pathology. Confirmation and validation using immunodetection in a replication set and an animal model.

SETTING

A multicenter European study (AddNeuroMed) and the Baltimore Longitudinal Study of Aging.

PARTICIPANTS

Patients with AD, subjects with mild cognitive impairment, and healthy controls with standardized clinical assessments and structural neuroimaging.

MAIN OUTCOME MEASURES

Association of plasma proteins with brain atrophy, disease severity, and rate of clinical progression. Extension studies in humans and transgenic mice tested the association between plasma proteins and brain amyloid.

RESULTS

Clusterin/apolipoprotein J was associated with atrophy of the entorhinal cortex, baseline disease severity, and rapid clinical progression in AD. Increased plasma concentration of clusterin was predictive of greater fibrillar amyloid-beta burden in the medial temporal lobe. Subjects with AD had increased clusterin messenger RNA in blood, but there was no effect of single-nucleotide polymorphisms in the gene encoding clusterin with gene or protein expression. APP/PS1 transgenic mice showed increased plasma clusterin, age-dependent increase in brain clusterin, as well as amyloid and clusterin colocalization in plaques.

CONCLUSIONS

These results demonstrate an important role of clusterin in the pathogenesis of AD and suggest that alterations in amyloid chaperone proteins may be a biologically relevant peripheral signature of AD.

Clusterin regulates β-amyloid toxicity via Dickkopf-1-driven induction of the wnt-PCP-JNK pathway

Although the mechanism of Aβ action in the pathogenesis of Alzheimer's disease (AD) has remained elusive, it is known to increase the expression of the antagonist of canonical wnt signalling, Dickkopf-1 (Dkk1), whereas the silencing of Dkk1 blocks Aβ neurotoxicity. We asked if clusterin, known to be regulated by wnt, is part of an Aβ/Dkk1 neurotoxic pathway. Knockdown of clusterin in primary neurons reduced Aβ toxicity and DKK1 upregulation and, conversely, Aβ increased intracellular clusterin and decreased clusterin protein secretion, resulting in the p53-dependent induction of DKK1. To further elucidate how the clusterin-dependent induction of Dkk1 by Aβ mediates neurotoxicity, we measured the effects of Aβ and Dkk1 protein on whole-genome expression in primary neurons, finding a common pathway suggestive of activation of wnt-planar cell polarity (PCP)-c-Jun N-terminal kinase (JNK) signalling leading to the induction of genes including EGR1 (early growth response-1), NAB2 (Ngfi-A-binding protein-2) and KLF10 (Krüppel-like factor-10) that, when individually silenced, protected against Aβ neurotoxicity and/or tau phosphorylation. Neuronal overexpression of Dkk1 in transgenic mice mimicked this Aβ-induced pathway and resulted in age-dependent increases in tau phosphorylation in hippocampus and cognitive impairment. Furthermore, we show that this Dkk1/wnt-PCP-JNK pathway is active in an Aβ-based mouse model of AD and in AD brain, but not in a tau-based mouse model or in frontotemporal dementia brain. Thus, we have identified a pathway whereby Aβ induces a clusterin/p53/Dkk1/wnt-PCP-JNK pathway, which drives the upregulation of several genes that mediate the development of AD-like neuropathologies, thereby providing new mechanistic insights into the action of Aβ in neurodegenerative diseases.

A novel multi-tissue RNA diagnostic of healthy ageing relates to cognitive health status

BACKGROUND

Diagnostics of the human ageing process may help predict future healthcare needs or guide preventative measures for tackling diseases of older age. We take a transcriptomics approach to build the first reproducible multi-tissue RNA expression signature by gene-chip profiling tissue from sedentary normal subjects who reached 65 years of age in good health.

RESULTS

One hundred and fifty probe-sets form an accurate classifier of young versus older muscle tissue and this healthy ageing RNA classifier performed consistently in independent cohorts of human muscle, skin and brain tissue (n = 594, AUC = 0.83-0.96) and thus represents a biomarker for biological age. Using the Uppsala Longitudinal Study of Adult Men birth-cohort (n = 108) we demonstrate that the RNA classifier is insensitive to confounding lifestyle biomarkers, while greater gene score at age 70 years is independently associated with better renal function at age 82 years and longevity. The gene score is 'up-regulated' in healthy human hippocampus with age, and when applied to blood RNA profiles from two large independent age-matched dementia case-control data sets (n = 717) the healthy controls have significantly greater gene scores than those with cognitive impairment. Alone, or when combined with our previously described prototype Alzheimer disease (AD) RNA 'disease signature', the healthy ageing RNA classifier is diagnostic for AD.

CONCLUSIONS

We identify a novel and statistically robust multi-tissue RNA signature of human healthy ageing that can act as a diagnostic of future health, using only a peripheral blood sample. This RNA signature has great potential to assist research aimed at finding treatments for and/or management of AD and other ageing-related conditions.

A histone acetylome-wide association study of Alzheimer's disease identifies disease-associated H3K27ac differences in the entorhinal cortex

We quantified genome-wide patterns of lysine H3K27 acetylation (H3K27ac) in entorhinal cortex samples from Alzheimer's disease (AD) cases and matched controls using chromatin immunoprecipitation and highly parallel sequencing. We observed widespread acetylomic variation associated with AD neuropathology, identifying 4,162 differential peaks (false discovery rate < 0.05) between AD cases and controls. Differentially acetylated peaks were enriched in disease-related biological pathways and included regions annotated to genes involved in the progression of amyloid-β and tau pathology (for example, APP, PSEN1, PSEN2, and MAPT), as well as regions containing variants associated with sporadic late-onset AD. Partitioned heritability analysis highlighted a highly significant enrichment of AD risk variants in entorhinal cortex H3K27ac peak regions. AD-associated variable H3K27ac was associated with transcriptional variation at proximal genes including CR1, GPR22, KMO, PIM3, PSEN1, and RGCC. In addition to identifying molecular pathways associated with AD neuropathology, we present a framework for genome-wide studies of histone modifications in complex disease.

Plasma proteins predict conversion to dementia from prodromal disease

BACKGROUND

The study aimed to validate previously discovered plasma biomarkers associated with AD, using a design based on imaging measures as surrogate for disease severity and assess their prognostic value in predicting conversion to dementia.

METHODS

Three multicenter cohorts of cognitively healthy elderly, mild cognitive impairment (MCI), and AD participants with standardized clinical assessments and structural neuroimaging measures were used. Twenty-six candidate proteins were quantified in 1148 subjects using multiplex (xMAP) assays.

RESULTS

Sixteen proteins correlated with disease severity and cognitive decline. Strongest associations were in the MCI group with a panel of 10 proteins predicting progression to AD (accuracy 87%, sensitivity 85%, and specificity 88%).

CONCLUSIONS

We have identified 10 plasma proteins strongly associated with disease severity and disease progression. Such markers may be useful for patient selection for clinical trials and assessment of patients with predisease subjective memory complaints.

Mitochondrial genes are altered in blood early in Alzheimer's disease

Although mitochondrial dysfunction is a consistent feature of Alzheimer's disease in the brain and blood, the molecular mechanisms behind these phenomena are unknown. Here we have replicated our previous findings demonstrating reduced expression of nuclear-encoded oxidative phosphorylation (OXPHOS) subunits and subunits required for the translation of mitochondrial-encoded OXPHOS genes in blood from people with Alzheimer's disease and mild cognitive impairment. Interestingly this was accompanied by increased expression of some mitochondrial-encoded OXPHOS genes, namely those residing closest to the transcription start site of the polycistronic heavy chain mitochondrial transcript (MT-ND1, MT-ND2, MT-ATP6, MT-CO1, MT-CO2, MT-C03) and MT-ND6 transcribed from the light chain. Further we show that mitochondrial DNA copy number was unchanged suggesting no change in steady-state numbers of mitochondria. We suggest that an imbalance in nuclear and mitochondrial genome-encoded OXPHOS transcripts may drive a negative feedback loop reducing mitochondrial translation and compromising OXPHOS efficiency, which is likely to generate damaging reactive oxygen species.

Elevated DNA methylation across a 48-kb region spanning the HOXA gene cluster is associated with Alzheimer's disease neuropathology

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⁻¹⁸).

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.

Genome-wide association with MRI atrophy measures as a quantitative trait locus for Alzheimer's disease

Alzheimer's disease (AD) is a progressive neurodegenerative disorder with considerable evidence suggesting an initiation of disease in the entorhinal cortex and hippocampus and spreading thereafter to the rest of the brain. In this study, we combine genetics and imaging data obtained from the Alzheimer's Disease Neuroimaging Initiative and the AddNeuroMed study. To identify genetic susceptibility loci for AD, we conducted a genome-wide study of atrophy in regions associated with neurodegeneration in this condition. We identified one single-nucleotide polymorphism (SNP) with a disease-specific effect associated with entorhinal cortical volume in an intron of the ZNF292 gene (rs1925690; P-value=2.6 × 10(-8); corrected P-value for equivalent number of independent quantitative traits=7.7 × 10(-8)) and an intergenic SNP, flanking the ARPP-21 gene, with an overall effect on entorhinal cortical thickness (rs11129640; P-value=5.6 × 10(-8); corrected P-value=1.7 × 10(-7)). Gene-wide scoring also highlighted PICALM as the most significant gene associated with entorhinal cortical thickness (P-value=6.7 × 10(-6)).

Mitochondrial dysfunction and immune activation are detectable in early Alzheimer's disease blood

Alzheimer's disease (AD), like other dementias, is characterized by progressive neuronal loss and neuroinflammation in the brain. The peripheral leukocyte response occurring alongside these brain changes has not been extensively studied, but might inform therapeutic approaches and provide relevant disease biomarkers. Using microarrays, we assessed blood gene expression alterations occurring in people with AD and those with mild cognitive changes at increased risk of developing AD. Of the 2,908 differentially expressed probes identified between the three groups (p < 0.01), a quarter were altered in blood from mild cognitive impairment (MCI) and AD subjects, relative to controls, suggesting a peripheral response to pathology may occur very early. There was strong evidence for mitochondrial dysfunction with decreased expression of many of the respiratory complex I-V genes and subunits of the core mitochondrial ribosome complex. This mirrors changes previously observed in AD brain. A number of genes encoding cell adhesion molecules were increased, along with other immune-related genes. These changes are consistent with leukocyte activation and their increased the transition from circulation into the brain. In addition to expression changes, we also found increased numbers of basophils in people with MCI and AD, and increased monocytes in people with an AD diagnosis. Taken together this study provides both an insight into the functional response of circulating leukocytes during neurodegeneration and also identifies potential targets such as the respiratory chain for designing and monitoring future therapeutic interventions using blood.

Inflammatory proteins in plasma are associated with severity of Alzheimer's disease

Markers of Alzheimer’s disease (AD) are being widely sought with a number of studies suggesting blood measures of inflammatory proteins as putative biomarkers. Here we report findings from a panel of 27 cytokines and related proteins in over 350 subjects with AD, subjects with Mild Cognitive Impairment (MCI) and elderly normal controls where we also have measures of longitudinal change in cognition and baseline neuroimaging measures of atrophy. In this study, we identify five inflammatory proteins associated with evidence of atrophy on MR imaging data particularly in whole brain, ventricular and entorhinal cortex measures. In addition, we observed six analytes that showed significant change (over a period of one year) in people with fast cognitive decline compared to those with intermediate and slow decline. One of these (IL-10) was also associated with brain atrophy in AD. In conclusion, IL-10 was associated with both clinical and imaging evidence of severity of disease and might therefore have potential to act as biomarker of disease progression.

Telomere length and risk of idiopathic pulmonary fibrosis and chronic obstructive pulmonary disease: a mendelian randomisation study

BACKGROUND

Idiopathic pulmonary fibrosis (IPF) is a fatal lung disease accounting for 1% of UK deaths. In the familial form of pulmonary fibrosis, causal genes have been identified in about 30% of cases, and a majority of these causal genes are associated with telomere maintenance. Prematurely shortened leukocyte telomere length is associated with IPF and chronic obstructive pulmonary disease (COPD), a disease with similar demographics and shared risk factors. Using mendelian randomisation, we investigated evidence supporting a causal role for short telomeres in IPF and COPD.

METHODS

Mendelian randomisation inference of telomere length causality was done for IPF (up to 1369 cases) and COPD (13 538 cases) against 435 866 controls of European ancestry in UK Biobank. Polygenic risk scores were calculated and two-sample mendelian randomisation analyses were done using seven genetic variants previously associated with telomere length, with replication analysis in an IPF cohort (2668 cases vs 8591 controls) and COPD cohort (15 256 cases vs 47 936 controls).

FINDINGS

In the UK Biobank, a genetically instrumented one-SD shorter telomere length was associated with higher odds of IPF (odds ratio [OR] 4·19, 95% CI 2·33-7·55; p=0·0031) but not COPD (1·07, 0·88-1·30; p=0·51). Similarly, an association was found in the IPF replication cohort (12·3, 5·05-30·1; p=0·0015) and not in the COPD replication cohort (1·04, 0·71-1·53; p=0·83). Meta-analysis of the two-sample mendelian randomisation results provided evidence inferring that shorter telomeres cause IPF (5·81 higher odds of IPF, 95% CI 3·56-9·50; p=2·19 × 10-12). There was no evidence to infer that telomere length caused COPD (OR 1·07, 95% CI 0·90-1·27; p=0·46).

INTERPRETATION

Cellular senescence is hypothesised as a major driving force in IPF and COPD; telomere shortening might be a contributory factor in IPF, suggesting divergent mechanisms in COPD. Defining a key role for telomere shortening enables greater focus in telomere-related diagnostics, treatments, and the search for a cure in IPF. Investigation of therapies that improve telomere length is warranted.

FUNDING

Medical Research Council.

Systemic inflammation modulates Fc receptor expression on microglia during chronic neurodegeneration

Chronic neurodegeneration is a major worldwide health problem, and it has been suggested that systemic inflammation can accelerate the onset and progression of clinical symptoms. A possible explanation is that systemic inflammation "switches" the phenotype of microglia from a relatively benign to a highly aggressive and tissue-damaging phenotype. The current study investigated the molecular mechanism underlying this microglia phenotype "switching." We show in mice with chronic neurodegeneration (ME7 prion model) that there is increased expression of receptors that have a key role in macrophage activation and associated signaling pathways, including TREM-2, Siglec-F, CD200R, and FcγRs. Systemic inflammation induced by LPS further increased protein levels of the activating FcγRIII and FcγRIV, but not of other microglial receptors, including the inhibitory FcγRII. In addition to these changes in receptor expression, IgG levels in the brain parenchyma were increased during chronic neurodegeneration, and these IgG levels further increased after systemic inflammation. γ-Chain-deficient mice show modified proinflammatory cytokine expression in the brain after systemic inflammation. We conclude that systemic inflammation during chronic neurodegeneration increases the expression levels of activating FcγR on microglia and thereby lowers the signaling threshold for Ab-mediated cell activation. At the same time, IgG influx into the brain could provide a cross-linking ligand resulting in excessive microglia activation that is detrimental to neurons already under threat by misfolded protein.

Recalibrating the epigenetic clock: implications for assessing biological age in the human cortex

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.

Cross-region reduction in 5-hydroxymethylcytosine in Alzheimer's disease brain

Epigenetic processes play a key role in the central nervous system and altered levels of 5-methylcytosine have been associated with a number of neurologic phenotypes, including Alzheimer's disease (AD). Recently, 3 additional cytosine modifications have been identified (5-hydroxymethylcytosine, 5-formylcytosine, and 5-carboxylcytosine), which are thought to be intermediate steps in the demethylation of 5-methylcytosine to unmodified cytosine. Little is known about the frequency of these modifications in the human brain during health or disease. In this study, we used immunofluorescence to confirm the presence of each modification in human brain and investigate their cross-tissue abundance in AD patients and elderly control samples. We identify a significant AD-associated decrease in global 5-hydroxymethylcytosine in entorhinal cortex and cerebellum, and differences in 5-formylcytosine levels between brain regions. Our study further implicates a role for epigenetic alterations in AD.

Longitudinal analyses of the DNA methylome in deployed military servicemen identify susceptibility loci for post-traumatic stress disorder

In order to determine the impact of the epigenetic response to traumatic stress on post-traumatic stress disorder (PTSD), this study examined longitudinal changes of genome-wide blood DNA methylation profiles in relation to the development of PTSD symptoms in two prospective military cohorts (one discovery and one replication data set). In the first cohort consisting of male Dutch military servicemen (n=93), the emergence of PTSD symptoms over a deployment period to a combat zone was significantly associated with alterations in DNA methylation levels at 17 genomic positions and 12 genomic regions. Evidence for mediation of the relation between combat trauma and PTSD symptoms by longitudinal changes in DNA methylation was observed at several positions and regions. Bioinformatic analyses of the reported associations identified significant enrichment in several pathways relevant for symptoms of PTSD. Targeted analyses of the significant findings from the discovery sample in an independent prospective cohort of male US marines (n=98) replicated the observed relation between decreases in DNA methylation levels and PTSD symptoms at genomic regions in ZFP57, RNF39 and HIST1H2APS2. Together, our study pinpoints three novel genomic regions where longitudinal decreases in DNA methylation across the period of exposure to combat trauma marks susceptibility for PTSD.

Alzheimer's disease-associated (hydroxy)methylomic changes in the brain and blood

BACKGROUND

Late-onset Alzheimer's disease (AD) is a complex multifactorial affliction, the pathogenesis of which is thought to involve gene-environment interactions that might be captured in the epigenome. The present study investigated epigenome-wide patterns of DNA methylation (5-methylcytosine, 5mC) and hydroxymethylation (5-hydroxymethylcytosine, 5hmC), as well as the abundance of unmodified cytosine (UC), in relation to AD.

RESULTS

We identified epigenetic differences in AD patients (n = 45) as compared to age-matched controls (n = 35) in the middle temporal gyrus, pertaining to genomic regions close to or overlapping with genes such as OXT (- 3.76% 5mC, pŠidák = 1.07E-06), CHRNB1 (+ 1.46% 5hmC, pŠidák = 4.01E-04), RHBDF2 (- 3.45% UC, pŠidák = 4.85E-06), and C3 (- 1.20% UC, pŠidák = 1.57E-03). In parallel, in an independent cohort, we compared the blood methylome of converters to AD dementia (n = 54) and non-converters (n = 42), at a preclinical stage. DNA methylation in the same region of the OXT promoter as found in the brain was found to be associated with subsequent conversion to AD dementia in the blood of elderly, non-demented individuals (+ 3.43% 5mC, pŠidák = 7.14E-04).

CONCLUSIONS

The implication of genome-wide significant differential methylation of OXT, encoding oxytocin, in two independent cohorts indicates it is a promising target for future studies on early biomarkers and novel therapeutic strategies in AD.

Parallel profiling of DNA methylation and hydroxymethylation highlights neuropathology-associated epigenetic variation in Alzheimer's disease

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.

Identification of cis-regulatory variation influencing protein abundance levels in human plasma

Proteins are central to almost all cellular processes, and dysregulation of expression and function is associated with a range of disorders. A number of studies in human have recently shown that genetic factors significantly contribute gene expression variation. In contrast, very little is known about the genetic basis of variation in protein abundance in man. Here, we assayed the abundance levels of proteins in plasma from 96 elderly Europeans using a new aptamer-based proteomic technology and performed genome-wide local (cis-) regulatory association analysis to identify protein quantitative trait loci (pQTL). We detected robust cis-associations for 60 proteins at a false discovery rate of 5%. The most highly significant single nucleotide polymorphism detected was rs7021589 (false discovery rate, 2.5 × 10(-12)), mapped within the gene coding sequence of Tenascin C (TNC). Importantly, we identified evidence of cis-regulatory variation for 20 previously disease-associated genes encoding protein, including variants with strong evidence of disease association show significant association with protein abundance levels. These results demonstrate that common genetic variants contribute to the differences in protein abundance levels in human plasma. Identification of pQTLs will significantly enhance our ability to discover and comprehend the biological and functional consequences of loci identified from genome-wide association study of complex traits. This is the first large-scale genetic association study of proteins in plasma measured using a novel, highly multiplexed slow off-rate modified aptamer (SOMAmer) proteomic platform.

A meta-analysis of epigenome-wide association studies in Alzheimer's disease highlights novel differentially methylated loci across cortex

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/ ).