Accelerated Pace of Aging in Schizophrenia: Five Case-Control Studies

BACKGROUND

Schizophrenia is associated with increased risk of developing multiple aging-related diseases, including metabolic, respiratory, and cardiovascular diseases, and Alzheimer's and related dementias, leading to the hypothesis that schizophrenia is accompanied by accelerated biological aging. This has been difficult to test because there is no widely accepted measure of biological aging. Epigenetic clocks are promising algorithms that are used to calculate biological age on the basis of information from combined cytosine-phosphate-guanine sites (CpGs) across the genome, but they have yielded inconsistent and often negative results about the association between schizophrenia and accelerated aging. Here, we tested the schizophrenia-aging hypothesis using a DNA methylation measure that is uniquely designed to predict an individual's rate of aging.

METHODS

We brought together 5 case-control datasets to calculate DunedinPACE (Pace of Aging Calculated from the Epigenome), a new measure trained on longitudinal data to detect differences between people in their pace of aging over time. Data were available from 1812 psychosis cases (schizophrenia or first-episode psychosis) and 1753 controls. Mean chronological age was 38.9 (SD = 13.6) years.

RESULTS

We observed consistent associations across datasets between schizophrenia and accelerated aging as measured by DunedinPACE. These associations were not attributable to tobacco smoking or clozapine medication.

CONCLUSIONS

Schizophrenia is accompanied by accelerated biological aging by midlife. This may explain the wide-ranging risk among people with schizophrenia for developing multiple different age-related physical diseases, including metabolic, respiratory, and cardiovascular diseases, and dementia. Measures of biological aging could prove valuable for assessing patients' risk for physical and cognitive decline and for evaluating intervention effectiveness.

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

Importance

Understanding antidepressant mechanisms could help design more effective and tolerated treatments.

Objective

Identify DNA methylation (DNAm) changes associated with antidepressant exposure.

Design

Case-control methylome-wide association studies (MWAS) of antidepressant exposure were performed from blood samples collected between 2006-2011 in Generation Scotland (GS). The summary statistics were tested for enrichment in specific tissues, gene ontologies and an independent MWAS in the Netherlands Study of Depression and Anxiety (NESDA). A methylation profile score (MPS) was derived and tested for its association with antidepressant exposure in eight independent cohorts, alongside prospective data from GS.

Setting

Cohorts; GS, NESDA, FTC, SHIP-Trend, FOR2107, LBC1936, MARS-UniDep, ALSPAC, E-Risk, and NTR.

Participants

Participants with DNAm data and self-report/prescription derived antidepressant exposure.

Main Outcomes and Measures

Whole-blood DNAm levels were assayed by the EPIC/450K Illumina array (9 studies, N exposed = 661, N unexposed = 9,575) alongside MBD-Seq in NESDA (N exposed = 398, N unexposed = 414). Antidepressant exposure was measured by self- report and/or antidepressant prescriptions.

Results

The self-report MWAS (N = 16,536, N exposed = 1,508, mean age = 48, 59% female) and the prescription-derived MWAS (N = 7,951, N exposed = 861, mean age = 47, 59% female), found hypermethylation at seven and four DNAm sites (p < 9.42x10 -8 ), respectively. The top locus was cg26277237 ( KANK1, p self-report = 9.3x10 -13 , p prescription = 6.1x10 -3 ). The self-report MWAS found a differentially methylated region, mapping to DGUOK-AS1 ( p adj = 5.0x10 -3 ) alongside significant enrichment for genes expressed in the amygdala, the "synaptic vesicle membrane" gene ontology and the top 1% of CpGs from the NESDA MWAS (OR = 1.39, p < 0.042). The MPS was associated with antidepressant exposure in meta-analysed data from external cohorts (N studies = 9, N = 10,236, N exposed = 661, f3 = 0.196, p < 1x10 -4 ).

Conclusions and Relevance

Antidepressant exposure is associated with changes in DNAm across different cohorts. Further investigation into these changes could inform on new targets for antidepressant treatments.

3 Key Points

Question: Is antidepressant exposure associated with differential whole blood DNA methylation?Findings: In this methylome-wide association study of 16,536 adults across Scotland, antidepressant exposure was significantly associated with hypermethylation at CpGs mapping to KANK1 and DGUOK-AS1. A methylation profile score trained on this sample was significantly associated with antidepressant exposure (pooled f3 [95%CI]=0.196 [0.105, 0.288], p < 1x10 -4 ) in a meta-analysis of external datasets. Meaning: Antidepressant exposure is associated with hypermethylation at KANK1 and DGUOK-AS1 , which have roles in mitochondrial metabolism and neurite outgrowth. If replicated in future studies, targeting these genes could inform the design of more effective and better tolerated treatments for depression.

Enrichment of the Local Synaptic Translatome for Genetic Risk Associated With Schizophrenia and Autism Spectrum Disorder

BACKGROUND

Genes that encode synaptic proteins or messenger RNA targets of the RNA-binding protein FMRP (fragile X messenger ribonucleoprotein) have been linked to schizophrenia and autism spectrum disorder (ASD) through the enrichment of genetic variants that confer risk for these disorders. FMRP binds many transcripts with synaptic functions and is thought to regulate their local translation, a process that enables rapid and compartmentalized protein synthesis required for development and plasticity.

METHODS

We used summary statistics from large-scale genome-wide association studies of schizophrenia (74,776 cases, 101,023 controls) and ASD (18,381 cases, 27,969 controls) to test the hypothesis that the subset of synaptic genes that encode localized transcripts is more strongly associated with each disorder than nonlocalized transcripts. We also postulated that this subset of synaptic genes is responsible for associations attributed to FMRP targets.

RESULTS

Schizophrenia associations were enriched in genes encoding localized synaptic transcripts compared to the remaining synaptic genes or to the remaining localized transcripts; this also applied to ASD associations, although only for transcripts observed after stimulation by fear conditioning. The genetic associations with either disorder captured by these gene sets were independent of those derived from FMRP targets. Schizophrenia association was related to FMRP interactions with messenger RNAs in somata, but not in dendrites, while ASD association was related to FMRP binding in either compartment.

CONCLUSIONS

Our data suggest that synaptic transcripts capable of local translation are particularly relevant to the pathogenesis of schizophrenia and ASD, but they do not characterize the associations attributed to current sets of FMRP targets.

Meta-analysis of epigenetic aging in schizophrenia reveals multifaceted relationships with age, sex, illness duration, and polygenic risk

BACKGROUND

The study of biological age acceleration may help identify at-risk individuals and reduce the rising global burden of age-related diseases. Using DNA methylation (DNAm) clocks, we investigated biological aging in schizophrenia (SCZ), a mental illness that is associated with an increased prevalence of age-related disabilities and morbidities. In a whole blood DNAm sample of 1090 SCZ cases and 1206 controls across four European cohorts, we performed a meta-analysis of differential aging using three DNAm clocks (i.e., Hannum, Horvath, and Levine). To dissect how DNAm aging contributes to SCZ, we integrated information on duration of illness and SCZ polygenic risk, as well as stratified our analyses by chronological age and biological sex.

RESULTS

We found that blood-based DNAm aging is significantly altered in SCZ independent from duration of the illness since onset. We observed sex-specific and nonlinear age effects that differed between clocks and point to possible distinct age windows of altered aging in SCZ. Most notably, intrinsic cellular age (Horvath clock) is decelerated in SCZ cases in young adulthood, while phenotypic age (Levine clock) is accelerated in later adulthood compared to controls. Accelerated phenotypic aging was most pronounced in women with SCZ carrying a high polygenic burden with an age acceleration of + 3.82 years (CI 2.02-5.61, P = 1.1E-03). Phenotypic aging and SCZ polygenic risk contributed additively to the illness and together explained up to 14.38% of the variance in disease status.

CONCLUSIONS

Our study contributes to the growing body of evidence of altered DNAm aging in SCZ and points to intrinsic age deceleration in younger adulthood and phenotypic age acceleration in later adulthood in SCZ. Since 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. Our study did not find that DNAm aging could be explained by the duration of illness of patients, but we did observe age- and sex-specific effects that warrant further investigation. Finally, our results show that combining genetic and epigenetic predictors can improve predictions of disease outcomes and may help with disease management in schizophrenia.

Quantifying the proportion of different cell types in the human cortex using DNA methylation profiles

BACKGROUND

Due to interindividual variation in the cellular composition of the human cortex, it is essential that covariates that capture these differences are included in epigenome-wide association studies using bulk tissue. As experimentally derived cell counts are often unavailable, computational solutions have been adopted to estimate the proportion of different cell types using DNA methylation data. Here, we validate and profile the use of an expanded reference DNA methylation dataset incorporating two neuronal and three glial cell subtypes for quantifying the cellular composition of the human cortex.

RESULTS

We tested eight reference panels containing different combinations of neuronal- and glial cell types and characterised their performance in deconvoluting cell proportions from computationally reconstructed or empirically derived human cortex DNA methylation data. Our analyses demonstrate that while these novel brain deconvolution models produce accurate estimates of cellular proportions from profiles generated on postnatal human cortex samples, they are not appropriate for the use in prenatal cortex or cerebellum tissue samples. Applying our models to an extensive collection of empirical datasets, we show that glial cells are twice as abundant as neuronal cells in the human cortex and identify significant associations between increased Alzheimer's disease neuropathology and the proportion of specific cell types including a decrease in NeuNNeg/SOX10Neg nuclei and an increase of NeuNNeg/SOX10Pos nuclei.

CONCLUSIONS

Our novel deconvolution models produce accurate estimates for cell proportions in the human cortex. These models are available as a resource to the community enabling the control of cellular heterogeneity in epigenetic studies of brain disorders performed on bulk cortex tissue.

Global Energy Metabolism Deficit in Alzheimer Disease Brain

BACKGROUND

The understanding of Alzheimer's disease (AD) has been dominated by the amyloid hypothesis. However, therapies targeting beta-amyloid have largely failed, generating interest in other potential pathogenic factors including energy metabolism.

OBJECTIVES

To interrogate canonical energy metabolism pathways from human prefrontal cortical tissue samples obtained from necropsy comparing AD and control.

DESIGN, SETTING, AND PARTICIPANTS

Postmortem pre-frontal cortical tissue from 10 subjects histologically diagnosed with AD and 10 control (CTRL) subjects was subjected to untargeted metabolomics to interrogate energy metabolism pathways. The samples were matched by age, sex, and post-mortem interval. Metabolite Measurements: Untargeted metabolomics analyses were via Metabolon®.

RESULTS

Glucose-derived energy metabolites in the glycolytic and pentose phosphate pathway and the ketone body β-hydroxybutyrate were uniformly decreased in AD brain vs. CTRL brain.

CONCLUSION

This pilot study aimed to identify energy metabolism abnormalities using untargeted brain metabolomics in two independent subject cohorts. Our study revealed a pattern of global energy deficit in AD brain, supporting a growing body of evidence of deficient energy metabolism in AD.

Uncertainty quantification of reference-based cellular deconvolution algorithms

The majority of epigenetic epidemiology studies to date have generated genome-wide profiles from bulk tissues (e.g., whole blood) however these are vulnerable to confounding from variation in cellular composition. Proxies for cellular composition can be mathematically derived from the bulk tissue profiles using a deconvolution algorithm; however, there is no method to assess the validity of these estimates for a dataset where the true cellular proportions are unknown. In this study, we describe, validate and characterize a sample level accuracy metric for derived cellular heterogeneity variables. The CETYGO score captures the deviation between a sample's DNA methylation profile and its expected profile given the estimated cellular proportions and cell type reference profiles. We demonstrate that the CETYGO score consistently distinguishes inaccurate and incomplete deconvolutions when applied to reconstructed whole blood profiles. By applying our novel metric to >6,300 empirical whole blood profiles, we find that estimating accurate cellular composition is influenced by both technical and biological variation. In particular, we show that when using a common reference panel for whole blood, less accurate estimates are generated for females, neonates, older individuals and smokers. Our results highlight the utility of a metric to assess the accuracy of cellular deconvolution, and describe how it can enhance studies of DNA methylation that are reliant on statistical proxies for cellular heterogeneity. To facilitate incorporating our methodology into existing pipelines, we have made it freely available as an R package (https://github.com/ds420/CETYGO).

Whole blood DNA methylation changes are associated with anti-TNF drug concentration in patients with Crohn's disease

BACKGROUND AND AIMS

Anti-TNF treatment failure in patients with inflammatory bowel disease (IBD) is common and frequently related to low drug concentrations. In order to identify patients who may benefit from dose optimisation at the outset of anti-TNF therapy, we sought to define epigenetic biomarkers in whole blood at baseline associated with anti-TNF drug concentrations at week 14.

METHODS

DNA methylation from 1,104 whole blood samples from 385 patients in the Personalised Anti-TNF Therapy in Crohn's disease (PANTS) study were assessed using the Illumina EPIC Beadchip (v1.0) at baseline, weeks 14, 30 and 54. We compared DNA methylation profiles in anti-TNF-treated patients who experienced primary non-response at week 14 and if they were assessed at subsequent time points, were not in remission at week 30 or 54 (infliximab n = 99, adalimumab n = 94), with patients who responded at week 14 and when assessed at subsequent time points, were in remission at week 30 or 54 (infliximab n = 99, adalimumab n = 93).

RESULTS

Overall, between baseline and week 14, we observed 4,999 differentially methylated probes (DMPs) annotated to 2376 genes following anti-TNF treatment. Pathway analysis identified 108 significant gene ontology terms enriched in biological processes related to immune system processes and responses.Epigenome-wide association (EWAS) analysis identified 323 DMPs annotated to 210 genes at baseline associated with higher anti-TNF drug concentrations at week 14. Of these, 125 DMPs demonstrated shared associations with other common traits (proportion of shared CpGs compared to DMPs) including body mass index (23.2%), followed by CRP (11.5%), smoking (7.4%), alcohol consumption per day (7.1%) and IBD type (6.8%). EWAS of primary non-response to anti-TNF identified 20 DMPs that were associated with both anti-TNF drug concentration and primary non-response to anti-TNF with a strong correlation of the coefficients (Spearman's rho = -0.94, p < 0.001).

CONCLUSION

Baseline DNA methylation profiles may be used as a predictor for anti-TNF drug concentration at week 14 to identify patients who may benefit from dose optimisation at the outset of anti-TNF therapy.

Epigenetic Age Acceleration in Frontotemporal Lobar Degeneration: A Comprehensive Analysis in the Blood and Brain

Frontotemporal lobar degeneration (FTLD) includes a heterogeneous group of disorders pathologically characterized by the degeneration of the frontal and temporal lobes. In addition to major genetic contributors of FTLD such as mutations in MAPT, GRN, and C9orf72, recent work has identified several epigenetic modifications including significant differential DNA methylation in DLX1, and OTUD4 loci. As aging remains one of the major risk factors for FTLD, we investigated the presence of accelerated epigenetic aging in FTLD compared to controls. We calculated epigenetic age in both peripheral blood and brain tissues of multiple FTLD subtypes using several DNA methylation clocks, i.e., DNAmClockMulti, DNAmClockHannum, DNAmClockCortical, GrimAge, and PhenoAge, and determined age acceleration and its association with different cellular proportions and clinical traits. Significant epigenetic age acceleration was observed in the peripheral blood of both frontotemporal dementia (FTD) and progressive supranuclear palsy (PSP) patients compared to controls with DNAmClockHannum, even after accounting for confounding factors. A similar trend was observed with both DNAmClockMulti and DNAmClockCortical in post-mortem frontal cortex tissue of PSP patients and in FTLD cases harboring GRN mutations. Our findings support that increased epigenetic age acceleration in the peripheral blood could be an indicator for PSP and to a smaller extent, FTD.

Brain DNA methylomic analysis of frontotemporal lobar degeneration reveals OTUD4 in shared dysregulated signatures across pathological subtypes

Frontotemporal lobar degeneration (FTLD) is an umbrella term describing the neuropathology of a clinically, genetically and pathologically heterogeneous group of diseases, including frontotemporal dementia (FTD) and progressive supranuclear palsy (PSP). Among the major FTLD pathological subgroups, FTLD with TDP-43 positive inclusions (FTLD-TDP) and FTLD with tau-positive inclusions (FTLD-tau) are the most common, representing about 90% of the cases. Although alterations in DNA methylation have been consistently associated with neurodegenerative diseases, including Alzheimer's disease and Parkinson's disease, little is known for FTLD and its heterogeneous subgroups and subtypes. The main goal of this study was to investigate DNA methylation variation in FTLD-TDP and FTLD-tau. We used frontal cortex genome-wide DNA methylation profiles from three FTLD cohorts (142 FTLD cases and 92 controls), generated using the Illumina 450K or EPIC microarrays. We performed epigenome-wide association studies (EWAS) for each cohort followed by meta-analysis to identify shared differentially methylated loci across FTLD subgroups/subtypes. In addition, we used weighted gene correlation network analysis to identify co-methylation signatures associated with FTLD and other disease-related traits. Wherever possible, we also incorporated relevant gene/protein expression data. After accounting for a conservative Bonferroni multiple testing correction, the EWAS meta-analysis revealed two differentially methylated loci in FTLD, one annotated to OTUD4 (5'UTR-shore) and the other to NFATC1 (gene body-island). Of these loci, OTUD4 showed consistent upregulation of mRNA and protein expression in FTLD. In addition, in the three independent co-methylation networks, OTUD4-containing modules were enriched for EWAS meta-analysis top loci and were strongly associated with the FTLD status. These co-methylation modules were enriched for genes implicated in the ubiquitin system, RNA/stress granule formation and glutamatergic synaptic signalling. Altogether, our findings identified novel FTLD-associated loci, and support a role for DNA methylation as a mechanism involved in the dysregulation of biological processes relevant to FTLD, highlighting novel potential avenues for therapeutic development.

Exploring the mediation of DNA methylation across the epigenome between childhood adversity and First Episode of Psychosis-findings from the EU-GEI study

ABTRACT

Studies conducted in psychotic disorders have shown that DNA-methylation (DNAm) is sensitive to the impact of Childhood Adversity (CA). However, whether it mediates the association between CA and psychosis is yet to be explored. Epigenome wide association studies (EWAS) using the Illumina Infinium-Methylation EPIC array in peripheral blood tissue from 366 First-episode of psychosis and 517 healthy controls was performed. Adversity scores were created for abuse, neglect and composite adversity with the Childhood Trauma Questionnaire (CTQ). Regressions examining (I) CTQ scores with psychosis; (II) with DNAm EWAS level and (III) between DNAm and caseness, adjusted for a variety of confounders were conducted. Divide-Aggregate Composite-null Test for the composite null-hypothesis of no mediation effect was conducted. Enrichment analyses were conducted with missMethyl package and the KEGG database. Our results show that CA was associated with psychosis (Composite: OR = 1.68; p = <0.001; abuse: OR = 2.16; p < 0.001; neglect: OR = 2.27; p = <0.001). None of the CpG sites significantly mediated the adversity-psychosis association after Bonferroni correction (p < 8.1 × 10-8). However, 28, 34 and 29 differentially methylated probes associated with 21, 27, 20 genes passed a less stringent discovery threshold (p < 5 × 10-5) for composite, abuse and neglect respectively, with a lack of overlap between abuse and neglect. These included genes previously associated to psychosis in EWAS studies, such as PANK1, SPEG TBKBP1, TSNARE1 or H2R. Downstream gene ontology analyses did not reveal any biological pathways that survived false discovery rate correction. Although at a non-significant level, DNAm changes in genes previously associated with schizophrenia in EWAS studies may mediate the CA-psychosis association. These results and associated involved processes such as mitochondrial or histaminergic disfunction, immunity or neural signalling requires replication in well powered samples. The lack of overlap between mediating genes associated with abuse and neglect suggests differential biological trajectories linking CA subtypes and psychosis.

A comparison of feature selection methodologies and learning algorithms in the development of a DNA methylation-based telomere length estimator

BACKGROUND

The field of epigenomics holds great promise in understanding and treating disease with advances in machine learning (ML) and artificial intelligence being vitally important in this pursuit. Increasingly, research now utilises DNA methylation measures at cytosine-guanine dinucleotides (CpG) to detect disease and estimate biological traits such as aging. Given the challenge of high dimensionality of DNA methylation data, feature-selection techniques are commonly employed to reduce dimensionality and identify the most important subset of features. In this study, our aim was to test and compare a range of feature-selection methods and ML algorithms in the development of a novel DNA methylation-based telomere length (TL) estimator. We utilised both nested cross-validation and two independent test sets for the comparisons.

RESULTS

We found that principal component analysis in advance of elastic net regression led to the overall best performing estimator when evaluated using a nested cross-validation analysis and two independent test cohorts. This approach achieved a correlation between estimated and actual TL of 0.295 (83.4% CI [0.201, 0.384]) on the EXTEND test data set. Contrastingly, the baseline model of elastic net regression with no prior feature reduction stage performed less well in general-suggesting a prior feature-selection stage may have important utility. A previously developed TL estimator, DNAmTL, achieved a correlation of 0.216 (83.4% CI [0.118, 0.310]) on the EXTEND data. Additionally, we observed that different DNA methylation-based TL estimators, which have few common CpGs, are associated with many of the same biological entities.

CONCLUSIONS

The variance in performance across tested approaches shows that estimators are sensitive to data set heterogeneity and the development of an optimal DNA methylation-based estimator should benefit from the robust methodological approach used in this study. Moreover, our methodology which utilises a range of feature-selection approaches and ML algorithms could be applied to other biological markers and disease phenotypes, to examine their relationship with DNA methylation and predictive value.

Novel CYP1B1-RMDN2 Alzheimer's disease locus identified by genome-wide association analysis of cerebral tau deposition on PET

Determining the genetic architecture of Alzheimer's disease (AD) pathologies can enhance mechanistic understanding and inform precision medicine strategies. Here, we performed a genome-wide association study of cortical tau quantified by positron emission tomography in 3,136 participants from 12 independent studies. The CYP1B1-RMDN2 locus was associated with tau deposition. The most significant signal was at rs2113389, which explained 4.3% of the variation in cortical tau, while APOE4 rs429358 accounted for 3.6%. rs2113389 was associated with higher tau and faster cognitive decline. Additive effects, but no interactions, were observed between rs2113389 and diagnosis, APOE4 , and Aβ positivity. CYP1B1 expression was upregulated in AD. rs2113389 was associated with higher CYP1B1 expression and methylation levels. Mouse model studies provided additional functional evidence for a relationship between CYP1B1 and tau deposition but not Aβ. These results may provide insight into the genetic basis of cerebral tau and novel pathways for therapeutic development in AD.

Genome-wide DNA methylation analysis of aggressive behaviour: a longitudinal population-based study

BACKGROUND

Human aggression is influenced by an interplay between genetic predisposition and experience across the life span. This interaction is thought to occur through epigenetic mechanisms, inducing differential gene expression, thereby moderating neuronal cell and circuit function, and thus shaping aggressive behaviour.

METHODS

Genome-wide DNA methylation (DNAm) levels were measured in peripheral blood obtained from 95 individuals participating in the Estonian Children Personality Behaviours and Health Study (ECPBHS) at 15 and 25 years of age. We examined the association between aggressive behaviour, as measured by Life History of Aggression (LHA) total score and DNAm levels both assessed at age 25. We further examined the pleiotropic effect of genetic variants regulating LHA-associated differentially methylated positions (DMPs) and multiple traits related to aggressive behaviours. Lastly, we tested whether the DNA methylomic loci identified in association with LHA at age 25 were also present at age 15.

RESULTS

We found one differentially methylated position (DMP) (cg17815886; p = 1.12 × 10^-8 ) and five differentially methylated regions (DMRs) associated with LHA after multiple testing adjustments. The DMP annotated to the PDLIM5 gene, and DMRs resided in the vicinity of four protein-encoding genes (TRIM10, GTF2H4, SLC45A4, B3GALT4) and a long intergenic non-coding RNA (LINC02068). We observed evidence for the colocalization of genetic variants associated with top DMPs and general cognitive function, educational attainment and cholesterol levels. Notably, a subset of the DMPs associated with LHA at age 25 also displayed altered DNAm patterns at age 15 with high accuracy in predicting aggression.

CONCLUSIONS

Our findings highlight the potential role of DNAm in the development of aggressive behaviours. We observed pleiotropic genetic variants associated with identified DMPs, and various traits previously established to be relevant in shaping aggression in humans. The concordance of DNAm signatures in adolescents and young adults may have predictive value for inappropriate and maladaptive aggression later in life.

An overview of DNA methylation-derived trait score methods and applications

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.

Epigenetic age acceleration is associated with oligodendrocyte proportions in MSA and control brain tissue

AIMS

Epigenetic clocks are widely applied as surrogates for biological age in different tissues and/or diseases, including several neurodegenerative diseases. Despite white matter (WM) changes often being observed in neurodegenerative diseases, no study has investigated epigenetic ageing in white matter.

METHODS

We analysed the performances of two DNA methylation-based clocks, DNAmClockMulti and DNAmClockCortical , in post-mortem WM tissue from multiple subcortical regions and the cerebellum, and in oligodendrocyte-enriched nuclei. We also examined epigenetic ageing in control and multiple system atrophy (MSA) (WM and mixed WM and grey matter), as MSA is a neurodegenerative disease comprising pronounced WM changes and α-synuclein aggregates in oligodendrocytes.

RESULTS

Estimated DNA methylation (DNAm) ages showed strong correlations with chronological ages, even in WM (e.g., DNAmClockCortical , r = [0.80-0.97], p < 0.05). However, performances and DNAm age estimates differed between clocks and brain regions. DNAmClockMulti significantly underestimated ages in all cohorts except in the MSA prefrontal cortex mixed tissue, whereas DNAmClockCortical tended towards age overestimations. Pronounced age overestimations in the oligodendrocyte-enriched cohorts (e.g., oligodendrocyte-enriched nuclei, p = 6.1 × 10-5 ) suggested that this cell type ages faster. Indeed, significant positive correlations were observed between estimated oligodendrocyte proportions and DNAm age acceleration estimated by DNAmClockCortical (r > 0.31, p < 0.05), and similar trends were obtained with DNAmClockMulti . Although increased age acceleration was observed in MSA compared with controls, no significant differences were detected upon adjustment for possible confounders (e.g., cell-type proportions).

CONCLUSIONS

Our findings show that oligodendrocyte proportions positively influence epigenetic age acceleration across brain regions and highlight the need to further investigate this in ageing and neurodegeneration.

Genome-wide characterization of mitochondrial DNA methylation in human brain

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.

A comparison of blood and brain-derived ageing and inflammation-related DNA methylation signatures and their association with microglial burdens

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

Uncovering the genetic architecture of broad antisocial behavior through a genome-wide association study meta-analysis

Despite the substantial heritability of antisocial behavior (ASB), specific genetic variants robustly associated with the trait have not been identified. The present study by the Broad Antisocial Behavior Consortium (BroadABC) meta-analyzed data from 28 discovery samples (N = 85,359) and five independent replication samples (N = 8058) with genotypic data and broad measures of ASB. We identified the first significant genetic associations with broad ASB, involving common intronic variants in the forkhead box protein P2 (FOXP2) gene (lead SNP rs12536335, p = 6.32 × 10-10). Furthermore, we observed intronic variation in Foxp2 and one of its targets (Cntnap2) distinguishing a mouse model of pathological aggression (BALB/cJ strain) from controls (BALB/cByJ strain). Polygenic risk score (PRS) analyses in independent samples revealed that the genetic risk for ASB was associated with several antisocial outcomes across the lifespan, including diagnosis of conduct disorder, official criminal convictions, and trajectories of antisocial development. We found substantial genetic correlations of ASB with mental health (depression rg = 0.63, insomnia rg = 0.47), physical health (overweight rg = 0.19, waist-to-hip ratio rg = 0.32), smoking (rg = 0.54), cognitive ability (intelligence rg = -0.40), educational attainment (years of schooling rg = -0.46) and reproductive traits (age at first birth rg = -0.58, father's age at death rg = -0.54). Our findings provide a starting point toward identifying critical biosocial risk mechanisms for the development of ASB.

DNA methylation signatures of Alzheimer's disease neuropathology in the cortex are primarily driven by variation in non-neuronal cell-types

Alzheimer'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.

Epigenome-wide association study of human frontal cortex identifies differential methylation in Lewy body pathology

Parkinson'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.

Functional characterization of the schizophrenia associated gene AS3MT identifies a role in neuronal development

Genome-wide association studies (GWAS) have identified multiple genomic regions associated with schizophrenia, although many variants reside in noncoding regions characterized by high linkage disequilibrium (LD) making the elucidation of molecular mechanisms challenging. A genomic region on chromosome 10q24 has been consistently associated with schizophrenia with risk attributed to the AS3MT gene. Although AS3MT is hypothesized to play a role in neuronal development and differentiation, work to fully understand the function of this gene has been limited. In this study we explored the function of AS3MT using a neuronal cell line (SH-SY5Y). We confirm previous findings of isoform specific expression of AS3MT during SH-SY5Y differentiation toward neuronal fates. Using CRISPR-Cas9 gene editing we generated AS3MT knockout SH-SY5Y cell lines and used RNA-seq to identify significant changes in gene expression in pathways associated with neuronal development, inflammation, extracellular matrix formation, and RNA processing, including dysregulation of other genes strongly implicated in schizophrenia. We did not observe any morphological changes in cell size and neurite length following neuronal differentiation and MAP2 immunocytochemistry. These results provide novel insights into the potential role of AS3MT in brain development and identify pathways through which genetic variation in this region may confer risk for schizophrenia.

Evaluation of nanopore sequencing for epigenetic epidemiology: a comparison with DNA methylation microarrays

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 revolutionise 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 aetiology of complex disease.

Genome-wide study of DNA methylation shows alterations in metabolic, inflammatory, and cholesterol pathways in ALS

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.

DunedinPACE, a DNA methylation biomarker of the pace of aging

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.

Common and rare variant association analyses in amyotrophic lateral sclerosis identify 15 risk loci with distinct genetic architectures and neuron-specific biology

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.

Full-length transcript sequencing of human and mouse cerebral cortex identifies widespread isoform diversity and alternative splicing

Alternative splicing is a post-transcriptional regulatory mechanism producing distinct mRNA molecules from a single pre-mRNA with a prominent role in the development and function of the central nervous system. We used long-read isoform sequencing to generate full-length transcript sequences in the human and mouse cortex. We identify novel transcripts not present in existing genome annotations, including transcripts mapping to putative novel (unannotated) genes and fusion transcripts incorporating exons from multiple genes. Global patterns of transcript diversity are similar between human and mouse cortex, although certain genes are characterized by striking differences between species. We also identify developmental changes in alternative splicing, with differential transcript usage between human fetal and adult cortex. Our data confirm the importance of alternative splicing in the cortex, dramatically increasing transcriptional diversity and representing an important mechanism underpinning gene regulation in the brain. We provide transcript-level data for human and mouse cortex as a resource to the scientific community.

•

There is widespread transcript diversity in the cortex and many novel transcripts

•

Some genes display big differences in isoform number between human and mouse cortex

•

There is evidence of differential transcript usage between human fetal and adult cortex

•

There are many novel isoforms of genes associated with human brain disease

Identical twins carry a persistent epigenetic signature of early genome programming

Monozygotic (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.

The association of epigenetic clocks in brain tissue with brain pathologies and common aging phenotypes

Epigenetic clocks are calculated by combining DNA methylation states across select CpG sites to estimate biologic age, and have been noted as the most successful markers of biologic aging to date. Yet, limited research has considered epigenetic clocks calculated in brain tissue. We used DNA methylation states in dorsolateral prefrontal cortex specimens from 721 older participants of the Religious Orders Study and Rush Memory and Aging Project, to calculate DNA methylation age using four established epigenetic clocks: Hannum, Horvath, PhenoAge, GrimAge, and a new Cortical clock. The four established clocks were trained in blood samples (Hannum, PhenoAge, GrimAge) or using 51 human tissue and cell types (Horvath); the recent Cortical clock is the first trained in postmortem cortical tissue. Participants were recruited beginning in 1994 (Religious Orders Study) and 1997 (Memory and Aging Project), and followed annually with questionnaires and clinical evaluations; brain specimens were obtained for 80-90% of participants. Mean age at death was 88.0 (SD 6.7) years. We used linear regression, logistic regression, and linear mixed models, to examine relations of epigenetic clock ages to neuropathologic and clinical aging phenotypes, controlling for chronologic age, sex, education, and depressive symptomatology. Hannum, Horvath, PhenoAge and Cortical clock ages were related to pathologic diagnosis of Alzheimer's disease (AD), as well as to Aβ load (a hallmark pathology of Alzheimer's disease). However, associations were substantially stronger for the Cortical than other clocks; for example, each standard deviation (SD) increase in Hannum, Horvath, and PhenoAge clock age was related to approximately 30% greater likelihood of pathologic AD (all p < 0.05), while each SD increase in Cortical age was related to 90% greater likelihood of pathologic AD (odds ratio = 1.91, 95% confidence interval 1.38, 2.62). Moreover, Cortical age was significantly related to other AD pathology (eg, mean tau tangle density, p = 0.003), and to odds of neocortical Lewy body pathology (for each SD increase in Cortical age, odds ratio = 2.00, 95% confidence 1.27, 3.17), although no clocks were related to cerebrovascular neuropathology. Cortical age was the only epigenetic clock significantly associated with the clinical phenotypes examined, from dementia to cognitive decline (5 specific cognitive systems, and a global cognitive measure averaging 17 tasks) to Parkinsonian signs. Overall, our findings provide evidence of the critical necessity for bespoke clocks of brain aging for advancing research to understand, and eventually prevent, neurodegenerative diseases of aging.

Sites of active gene regulation in the prenatal frontal cortex and their role in neuropsychiatric disorders

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.

Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation

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.

The molecular etiology of Alzheimer's disease

Alzheimer's disease (AD) is a growing global healthcare epidemic. Owing to advances in technology, genome-scale studies of various layers of molecular information have been undertaken in recent years and robust variation in key loci have now been published and reproduced by others. This mini-symposium highlights four key areas of current research in the field of molecular biology in AD, including articles focused on large-scale genomic profiling, epigenetic research, integrative multi-omic approaches and how these can be appropriately modeled to address reverse causality. This mini-symposium provides a timely update on research focused on elucidating the molecular etiology of AD to date and highlights new methodological advances that could enable neuroscientists to identify novel therapeutic targets.

A machine learning case-control classifier for schizophrenia based on DNA methylation in blood

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

DNA methylation-based sex classifier to predict sex and identify sex chromosome aneuploidy

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<0.01) sex-associated CpG sites present on both 450K and EPIC arrays, and constructed a sex classifier based on the two first principal components of the DNA methylation data of sex-associated probes mapped on 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.

Conclusions

This proposed sex classifier not only can be used for sex predictions but also applied to identify samples with sex chromosome aneuploidy, and it is freely and easily accessible by calling the ‘estimateSex’ function from the newest wateRmelon Bioconductor package (https://github.com/schalkwyk/wateRmelon).

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-07675-2).

Novel epigenetic clock for fetal brain development predicts prenatal age for cellular stem cell models and derived neurons

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.

Characterizing the properties of bisulfite sequencing data: maximizing power and sensitivity to identify between-group differences in DNA methylation

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.

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

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

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

A central role for anterior cingulate cortex in the control of pathological aggression

Controlling aggression is a crucial skill in social species like rodents and humans and has been associated with anterior cingulate cortex (ACC). Here, we directly link the failed regulation of aggression in BALB/cJ mice to ACC hypofunction. We first show that ACC in BALB/cJ mice is structurally degraded: neuron density is decreased, with pervasive neuron death and reactive astroglia. Gene-set enrichment analysis suggested that this process is driven by neuronal degeneration, which then triggers toxic astrogliosis. cFos expression across ACC indicated functional consequences: during aggressive encounters, ACC was engaged in control mice, but not BALB/cJ mice. Chemogenetically activating ACC during aggressive encounters drastically suppressed pathological aggression but left species-typical aggression intact. The network effects of our chemogenetic perturbation suggest that this behavioral rescue is mediated by suppression of amygdala and hypothalamus and activation of mediodorsal thalamus. Together, these findings highlight the central role of ACC in curbing pathological aggression.

Meta-analysis of genome-wide DNA methylation identifies shared associations across neurodegenerative disorders

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.

Assessing the co-variability of DNA methylation across peripheral cells and tissues: Implications for the interpretation of findings in epigenetic epidemiology

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.

As epigenetic variation is cell-type specific, an ongoing challenge in epigenetic epidemiology is how to interpret studies performed using bulk tissue (for example, whole blood) which comprises a mix of different cell types. In this study, we identified major differences in DNA methylation (DNAm) across multiple peripheral tissues and different blood cell types, with each sample type being characterized by a unique signature across multiple genomic loci. We demonstrate how these differences influence commonly used prediction scores derived from DNAm data for age and tobacco smoking, with estimates for the same individual being highly variable across tissues and cell types. Our results enabled us to assess the extent to which variable DNAm in each individual blood cell type relates to variation measured in whole blood. We found that although individual blood cell types predict more of the variation in DNAm in whole blood compared to buccal and nasal epithelial cells, the actual proportion of variance explained is relatively small, except for at sites where DNAm is under genetic control. Our data indicate that for most sites variation in multiple blood cell types additively combines to drive variation in DNAm measured in whole blood. Of note, for a subset of sites, variation in DNAm detected in whole blood can be attributed to a specific blood cell type, potentially facilitating the interpretation of EWAS findings.

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.

DNA methylation meta-analysis reveals cellular alterations in psychosis and markers of treatment-resistant schizophrenia

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.

The histone modification H3K4me3 is altered at the ANK1 locus in Alzheimer's disease brain

Several epigenome-wide association studies of DNA methylation have highlighted altered DNA methylation in the ANK1 gene in Alzheimer's disease (AD) brain samples. However, no study has specifically examined ANK1 histone modifications in the disease. We use chromatin immunoprecipitation-qPCR to quantify tri-methylation at histone 3 lysine 4 (H3K4me3) and 27 (H3K27me3) in the ANK1 gene in entorhinal cortex from donors with high (n = 59) or low (n = 29) Alzheimer's disease pathology. We demonstrate decreased levels of H3K4me3, a marker of active gene transcription, with no change in H3K27me3, a marker of inactive genes. H3K4me3 is negatively correlated with DNA methylation in specific regions of the ANK1 gene. Our study suggests that the ANK1 gene shows altered epigenetic marks indicative of reduced gene activation in Alzheimer's disease.

4D flow magnetic resonance imaging to assess left atrial haemodynamics in healthy and hypertrophic subjects

Funding Acknowledgements

Type of funding sources: Public grant(s) – National budget only. Main funding source(s): - University, research centre and hospital foundation grants for the contracting of new research staff (FI 2020) - Spanish Ministry of Economy and Competitiveness Retos investigacion project

Introduction

The assessment of the left atrium (LA) haemodynamics is key to better understand the development of LA-related pathological processes. In this regard 4D flow magnetic resonance imaging (MRI) can provide complementary information to standard Doppler echocardiographic studies and identify complex blood flow patterns. Yet, until recently, the left atrium (LA) has been largely left aside in 4D flow MRI studies.

Purpose

We aimed at assessing the LA haemodynamics of healthy and hypertrophic cardiomyopathy (HCM) subjects with a qualitative visualization of flow patterns and deriving quantitative indices related to ventricular dysfunction from pulmonary veins (PV)  and mitral valve (MV) velocity profiles.

Methods

Segmentation was performed directly over 4D flow angiograms. A total of 20 cases were processed, 11 healthy and 9 HCM subjects. 4D velocity matrices were masked with the segmented mask to isolate LA haemodynamics. Velocity profiles were then obtained in the PV and MV and integrated over planes perpendicular to the lumen of the vessels to create velocity spectrograms. Fourier spectral analysis was applied to the velocity curves to highlight differences that might go unnoticed in the time domain. In addition, the Q-Criterion was computed for vortex identification, visually inspecting both cohorts across the whole cardiac cycle.

Results

Fourier spectral analysis of the velocity curves suggested that overall, healthy patients have higher dynamic range of the velocity curves. It can be observed in Figure 1, that the usual E/A MV velocity pattern is preserved in 10 of the 11 healthy subjects while 5 of the HCM patients present significant alterations of said curve. In fact, patients 4, 6, 7 and 8 seem to present a 3 peaked MV velocity curve. The vortex analysis identified 3 main types of vortices in healthy subjects: a ‘filling’ systolic vortex (10/11) arising near the most dominant PV (usually the left superior PV) as seen in Figure 2; a conduit phase vortex (7/11), similar in nature to the preceding systolic vortex; and an E-wave vortex (9/11) attached to the LA ostium. Four of the HCM patients (out of the five with altered MV velocity profile) also showed a systolic vortex, but with more complex blood flow patterns and emerging far from the PVs. One of such vortices is shown in Figure 2, composed of two distinct eddies near the MV. The E-wave vortex was also observed but was less predominant than in healthy subjects (3/9).

Conclusions

4D Flow analysis of the LA is feasible and might hold promise in the understanding of the complex haemodynamics in ventricular dysfunction.

Abstract Figure. Velocity Spectrograms and Vortices