Blood-based epigenome-wide analyses of chronic low-grade inflammation across diverse population cohorts

Chronic inflammation is a hallmark of age-related disease states. The effectiveness of inflammatory proteins including C-reactive protein (CRP) in assessing long-term inflammation is hindered by their phasic nature. DNA methylation (DNAm) signatures of CRP may act as more reliable markers of chronic inflammation. We show that inter-individual differences in DNAm capture 50% of the variance in circulating CRP (N = 17,936, Generation Scotland). We develop a series of DNAm predictors of CRP using state-of-the-art algorithms. An elastic-net-regression-based predictor outperformed competing methods and explained 18% of phenotypic variance in the Lothian Birth Cohort of 1936 (LBC1936) cohort, doubling that of existing DNAm predictors. DNAm predictors performed comparably in four additional test cohorts (Avon Longitudinal Study of Parents and Children, Health for Life in Singapore, Southall and Brent Revisited, and LBC1921), including for individuals of diverse genetic ancestry and different age groups. The best-performing predictor surpassed assay-measured CRP and a genetic score in its associations with 26 health outcomes. Our findings forge new avenues for assessing chronic low-grade inflammation in diverse populations.

Revealing uncertainty in the status of biodiversity change

Biodiversity faces unprecedented threats from rapid global change1. Signals of biodiversity change come from time-series abundance datasets for thousands of species over large geographic and temporal scales. Analyses of these biodiversity datasets have pointed to varied trends in abundance, including increases and decreases. However, these analyses have not fully accounted for spatial, temporal and phylogenetic structures in the data. Here, using a new statistical framework, we show across ten high-profile biodiversity datasets2-11 that increases and decreases under existing approaches vanish once spatial, temporal and phylogenetic structures are accounted for. This is a consequence of existing approaches severely underestimating trend uncertainty and sometimes misestimating the trend direction. Under our revised average abundance trends that appropriately recognize uncertainty, we failed to observe a single increasing or decreasing trend at 95% credible intervals in our ten datasets. This emphasizes how little is known about biodiversity change across vast spatial and taxonomic scales. Despite this uncertainty at vast scales, we reveal improved local-scale prediction accuracy by accounting for spatial, temporal and phylogenetic structures. Improved prediction offers hope of estimating biodiversity change at policy-relevant scales, guiding adaptive conservation responses.

Morphological profiling in human neural progenitor cells classifies hits in a pilot drug screen for Alzheimer's disease

Alzheimer's disease accounts for 60-70% of dementia cases. Current treatments are inadequate and there is a need to develop new approaches to drug discovery. Recently, in cancer, morphological profiling has been used in combination with high-throughput screening of small-molecule libraries in human cells in vitro. To test feasibility of this approach for Alzheimer's disease, we developed a cell morphology-based drug screen centred on the risk gene, SORL1 (which encodes the protein SORLA). Increased Alzheimer's disease risk has been repeatedly linked to variants in SORL1, particularly those conferring loss or decreased expression of SORLA, and lower SORL1 levels are observed in post-mortem brain samples from individuals with Alzheimer's disease. Consistent with its role in the endolysosomal pathway, SORL1 deletion is associated with enlarged endosomes in neural progenitor cells and neurons. We, therefore, hypothesized that multi-parametric, image-based cell phenotyping would identify features characteristic of SORL1 deletion. An automated morphological profiling method (Cell Painting) was adapted to neural progenitor cells and used to determine the phenotypic response of SORL1-/- neural progenitor cells to treatment with compounds from a small internationally approved drug library (TargetMol, 330 compounds). We detected distinct phenotypic signatures for SORL1-/- neural progenitor cells compared to isogenic wild-type controls. Furthermore, we identified 16 compounds (representing 14 drugs) that reversed the mutant morphological signatures in neural progenitor cells derived from three SORL1-/- induced pluripotent stem cell sub-clones. Network pharmacology analysis revealed the 16 compounds belonged to five mechanistic groups: 20S proteasome, aldehyde dehydrogenase, topoisomerase I and II, and DNA synthesis inhibitors. Enrichment analysis identified DNA synthesis/damage/repair, proteases/proteasome and metabolism as key pathways/biological processes. Prediction of novel targets revealed enrichment in pathways associated with neural cell function and Alzheimer's disease. Overall, this work suggests that (i) a quantitative phenotypic metric can distinguish induced pluripotent stem cell-derived SORL1-/- neural progenitor cells from isogenic wild-type controls and (ii) phenotypic screening combined with multi-parametric high-content image analysis is a viable option for drug repurposing and discovery in this human neural cell model of Alzheimer's disease.

Epigenetic Contributions to Clinical Risk Prediction of Cardiovascular Disease

BACKGROUND

Cardiovascular disease (CVD) is among the leading causes of death worldwide. The discovery of new omics biomarkers could help to improve risk stratification algorithms and expand our understanding of molecular pathways contributing to the disease. Here, ASSIGN-a cardiovascular risk prediction tool recommended for use in Scotland-was examined in tandem with epigenetic and proteomic features in risk prediction models in ≥12 657 participants from the Generation Scotland cohort.

METHODS

Previously generated DNA methylation-derived epigenetic scores (EpiScores) for 109 protein levels were considered, in addition to both measured levels and an EpiScore for cTnI (cardiac troponin I). The associations between individual protein EpiScores and the CVD risk were examined using Cox regression (ncases≥1274; ncontrols≥11 383) and visualized in a tailored R application. Splitting the cohort into independent training (n=6880) and test (n=3659) subsets, a composite CVD EpiScore was then developed.

RESULTS

Sixty-five protein EpiScores were associated with incident CVD independently of ASSIGN and the measured concentration of cTnI (P<0.05), over a follow-up of up to 16 years of electronic health record linkage. The most significant EpiScores were for proteins involved in metabolic, immune response, and tissue development/regeneration pathways. A composite CVD EpiScore (based on 45 protein EpiScores) was a significant predictor of CVD risk independent of ASSIGN and the concentration of cTnI (hazard ratio, 1.32; P=3.7×10-3; 0.3% increase in C-statistic).

CONCLUSIONS

EpiScores for circulating protein levels are associated with CVD risk independent of traditional risk factors and may increase our understanding of the etiology of the disease.

Comprehensive evaluation of smoking exposures and their interactions on DNA methylation

BACKGROUND

Smoking impacts DNA methylation, but data are lacking on smoking-related differential methylation by sex or dietary intake, recent smoking cessation (<1 year), persistence of differential methylation from in utero smoking exposure, and effects of environmental tobacco smoke (ETS).

METHODS

We meta-analysed data from up to 15,014 adults across 5 cohorts with DNA methylation measured in blood using Illumina's EPIC array for current smoking (2560 exposed), quit < 1 year (500 exposed), in utero (286 exposed), and ETS exposure (676 exposed). We also evaluated the interaction of current smoking with sex or diet (fibre, folate, and vitamin C).

FINDINGS

Using false discovery rate (FDR < 0.05), 65,857 CpGs were differentially methylated in relation to current smoking, 4025 with recent quitting, 594 with in utero exposure, and 6 with ETS. Most current smoking CpGs attenuated within a year of quitting. CpGs related to in utero exposure in adults were enriched for those previously observed in newborns. Differential methylation by current smoking at 4-71 CpGs may be modified by sex or dietary intake. Nearly half (35-50%) of differentially methylated CpGs on the 450 K array were associated with blood gene expression. Current smoking and in utero smoking CpGs implicated 3049 and 1067 druggable targets, including chemotherapy drugs.

INTERPRETATION

Many smoking-related methylation sites were identified with Illumina's EPIC array. Most signals revert to levels observed in never smokers within a year of cessation. Many in utero smoking CpGs persist into adulthood. Smoking-related druggable targets may provide insights into cancer treatment response and shared mechanisms across smoking-related diseases.

FUNDING

Intramural Research Program of the National Institutes of Health, Norwegian Ministry of Health and Care Services and the Ministry of Education and Research, Chief Scientist Office of the Scottish Government Health Directorates and the Scottish Funding Council, Medical Research Council UK and the Wellcome Trust.

Blood-based epigenome-wide analyses of 19 common disease states: A longitudinal, population-based linked cohort study of 18,413 Scottish individuals

BACKGROUND

DNA methylation is a dynamic epigenetic mechanism that occurs at cytosine-phosphate-guanine dinucleotide (CpG) sites. Epigenome-wide association studies (EWAS) investigate the strength of association between methylation at individual CpG sites and health outcomes. Although blood methylation may act as a peripheral marker of common disease states, previous EWAS have typically focused only on individual conditions and have had limited power to discover disease-associated loci. This study examined the association of blood DNA methylation with the prevalence of 14 disease states and the incidence of 19 disease states in a single population of over 18,000 Scottish individuals.

METHODS AND FINDINGS

DNA methylation was assayed at 752,722 CpG sites in whole-blood samples from 18,413 volunteers in the family-structured, population-based cohort study Generation Scotland (age range 18 to 99 years). EWAS tested for cross-sectional associations between baseline CpG methylation and 14 prevalent disease states, and for longitudinal associations between baseline CpG methylation and 19 incident disease states. Prevalent cases were self-reported on health questionnaires at the baseline. Incident cases were identified using linkage to Scottish primary (Read 2) and secondary (ICD-10) care records, and the censoring date was set to October 2020. The mean time-to-diagnosis ranged from 5.0 years (for chronic pain) to 11.7 years (for Coronavirus Disease 2019 (COVID-19) hospitalisation). The 19 disease states considered in this study were selected if they were present on the World Health Organisation's 10 leading causes of death and disease burden or included in baseline self-report questionnaires. EWAS models were adjusted for age at methylation typing, sex, estimated white blood cell composition, population structure, and 5 common lifestyle risk factors. A structured literature review was also conducted to identify existing EWAS for all 19 disease states tested. The MEDLINE, Embase, Web of Science, and preprint servers were searched to retrieve relevant articles indexed as of March 27, 2023. Fifty-four of approximately 2,000 indexed articles met our inclusion criteria: assayed blood-based DNA methylation, had >20 individuals in each comparison group, and examined one of the 19 conditions considered. First, we assessed whether the associations identified in our study were reported in previous studies. We identified 69 associations between CpGs and the prevalence of 4 conditions, of which 58 were newly described. The conditions were breast cancer, chronic kidney disease, ischemic heart disease, and type 2 diabetes mellitus. We also uncovered 64 CpGs that associated with the incidence of 2 disease states (COPD and type 2 diabetes), of which 56 were not reported in the surveyed literature. Second, we assessed replication across existing studies, which was defined as the reporting of at least 1 common site in >2 studies that examined the same condition. Only 6/19 disease states had evidence of such replication. The limitations of this study include the nonconsideration of medication data and a potential lack of generalizability to individuals that are not of Scottish and European ancestry.

CONCLUSIONS

We discovered over 100 associations between blood methylation sites and common disease states, independently of major confounding risk factors, and a need for greater standardisation among EWAS on human disease.

Phenome-wide analyses identify an association between the parent-of-origin effects dependent methylome and the rate of aging in humans

BACKGROUND

The variation in the rate at which humans age may be rooted in early events acting through the genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions enriched for genetically controlled imprinting effects (the typical type of POE) and regions influenced by environmental effects associated with parents (the atypical POE). This part of the methylome is heavily influenced by early events, making it a potential route connecting early exposures, the epigenome, and aging. We aim to test the association of POE-CpGs with early and later exposures and subsequently with health-related phenotypes and adult aging.

RESULTS

We perform a phenome-wide association analysis for the POE-influenced methylome using GS:SFHS (Ndiscovery = 5087, Nreplication = 4450). We identify and replicate 92 POE-CpG-phenotype associations. Most of the associations are contributed by the POE-CpGs belonging to the atypical class where the most strongly enriched associations are with aging (DNAmTL acceleration), intelligence, and parental (maternal) smoking exposure phenotypes. A proportion of the atypical POE-CpGs form co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased within-module methylation connectivity with age. The atypical POE-CpGs also display high levels of methylation heterogeneity, fast information loss with age, and a strong correlation with CpGs contained within epigenetic clocks.

CONCLUSIONS

These results identify the association between the atypical POE-influenced methylome and aging and provide new evidence for the "early development of origin" hypothesis for aging in humans.

Development and validation of DNA methylation scores in two European cohorts augment 10-year risk prediction of type 2 diabetes

Type 2 diabetes mellitus (T2D) presents a major health and economic burden that could be alleviated with improved early prediction and intervention. While standard risk factors have shown good predictive performance, we show that the use of blood-based DNA methylation information leads to a significant improvement in the prediction of 10-year T2D incidence risk. Previous studies have been largely constrained by linear assumptions, the use of cytosine-guanine pairs one-at-a-time and binary outcomes. We present a flexible approach (via an R package, MethylPipeR) based on a range of linear and tree-ensemble models that incorporate time-to-event data for prediction. Using the Generation Scotland cohort (training set n_cases = 374, n_controls = 9,461; test set n_cases = 252, n_controls = 4,526) our best-performing model (area under the receiver operating characteristic curve (AUC) = 0.872, area under the precision-recall curve (PRAUC) = 0.302) showed notable improvement in 10-year onset prediction beyond standard risk factors (AUC = 0.839, precision-recall AUC = 0.227). Replication was observed in the German-based KORA study (n = 1,451, n_cases = 142, P = 1.6 × 10^-5).

Refining epigenetic prediction of chronological and biological age

BACKGROUND

Epigenetic clocks can track both chronological age (cAge) and biological age (bAge). The latter is typically defined by physiological biomarkers and risk of adverse health outcomes, including all-cause mortality. As cohort sample sizes increase, estimates of cAge and bAge become more precise. Here, we aim to develop accurate epigenetic predictors of cAge and bAge, whilst improving our understanding of their epigenomic architecture.

METHODS

First, we perform large-scale (N = 18,413) epigenome-wide association studies (EWAS) of chronological age and all-cause mortality. Next, to create a cAge predictor, we use methylation data from 24,674 participants from the Generation Scotland study, the Lothian Birth Cohorts (LBC) of 1921 and 1936, and 8 other cohorts with publicly available data. In addition, we train a predictor of time to all-cause mortality as a proxy for bAge using the Generation Scotland cohort (1214 observed deaths). For this purpose, we use epigenetic surrogates (EpiScores) for 109 plasma proteins and the 8 component parts of GrimAge, one of the current best epigenetic predictors of survival. We test this bAge predictor in four external cohorts (LBC1921, LBC1936, the Framingham Heart Study and the Women's Health Initiative study).

RESULTS

Through the inclusion of linear and non-linear age-CpG associations from the EWAS, feature pre-selection in advance of elastic net regression, and a leave-one-cohort-out (LOCO) cross-validation framework, we obtain cAge prediction with a median absolute error equal to 2.3 years. Our bAge predictor was found to slightly outperform GrimAge in terms of the strength of its association to survival (HRGrimAge = 1.47 [1.40, 1.54] with p = 1.08 × 10-52, and HRbAge = 1.52 [1.44, 1.59] with p = 2.20 × 10-60). Finally, we introduce MethylBrowsR, an online tool to visualise epigenome-wide CpG-age associations.

CONCLUSIONS

The integration of multiple large datasets, EpiScores, non-linear DNAm effects, and new approaches to feature selection has facilitated improvements to the blood-based epigenetic prediction of biological and chronological age.

Additive Effects of Stress and Alcohol Exposure on Accelerated Epigenetic Aging in Alcohol Use Disorder

BACKGROUND

Stress contributes to premature aging and susceptibility to alcohol use disorder (AUD), and AUD itself is a factor in premature aging; however, the interrelationships of stress, AUD, and premature aging are poorly understood.

METHODS

We constructed a composite score of stress from 13 stress-related outcomes in a discovery cohort of 317 individuals with AUD and control subjects. We then developed a novel methylation score of stress (MS stress) as a proxy of composite score of stress comprising 211 CpGs selected using a penalized regression model. The effects of MS stress on health outcomes and epigenetic aging were assessed in a sample of 615 patients with AUD and control subjects using epigenetic clocks and DNA methylation-based telomere length. Statistical analysis with an additive model using MS stress and a MS for alcohol consumption (MS alcohol) was conducted. Results were replicated in 2 independent cohorts (Generation Scotland, N = 7028 and the Grady Trauma Project, N = 795).

RESULTS

Composite score of stress and MS stress were strongly associated with heavy alcohol consumption, trauma experience, epigenetic age acceleration (EAA), and shortened DNA methylation-based telomere length in AUD. Together, MS stress and MS alcohol additively showed strong stepwise increases in EAA. Replication analyses showed robust association between MS stress and EAA in the Generation Scotland and Grady Trauma Project cohorts.

CONCLUSIONS

A methylation-derived score tracking stress exposure is associated with various stress-related phenotypes and EAA. Stress and alcohol have additive effects on aging, offering new insights into the pathophysiology of premature aging in AUD and, potentially, other aspects of gene dysregulation in this disorder.

Phenome-wide analysis identifies parent-of-origin effects on the human methylome associated with changes in the rate of aging

Variation in the rate at which humans age may be rooted in early life events acting through genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions either enriched for genetically controlled imprinting effects (the typical type of POE) or atypical POE introduced by environmental effects associated with parents. This part of the methylome is heavily influenced by early life events, making it a potential route connecting early environmental exposures, the epigenome and the rate of aging. Here, we aim to test the association of POE-influenced methylation of CpG dinucleotides (POE-CpG sites) with early and later environmental exposures and subsequently with health-related phenotypes and adult aging phenotypes. We do this by performing phenome-wide association analyses of the POE-influenced methylome using a large family-based population cohort (GS:SFHS, Ndiscovery=5,087, Nreplication=4,450). At the single CpG level, 92 associations of POE-CpGs with phenotypic variation were identified and replicated. Most of the associations were contributed by POE-CpGs belonging to the atypical class and the most strongly enriched associations were with aging (DNAmTL acceleration), intelligence and parental (maternal) smoking exposure phenotypes. We further found that a proportion of the atypical-POE-CpGs formed co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased internal module connectivity (strength of methylation correlation across constituent CpGs) with age. Atypical POE-CpGs also displayed high levels of methylation heterogeneity and epigenetic drift (i.e. information loss with age) and a strong correlation with CpGs contained within epigenetic clocks. These results identified associations between the atypical-POE-influenced methylome and aging and provided new evidence for the "early development of origin" hypothesis for aging in humans.

Loss of SORCS2 is Associated with Neuronal DNA Double-Strand Breaks

SORCS2 is one of five proteins that constitute the Vps10p-domain receptor family. Members of this family play important roles in cellular processes linked to neuronal survival, differentiation and function. Genetic and functional studies implicate SORCS2 in cognitive function, as well as in neurodegenerative and psychiatric disorders. DNA damage and DNA repair deficits are linked to ageing and neurodegeneration, and transient neuronal DNA double-strand breaks (DSBs) also occur as a result of neuronal activity. Here, we report a novel role for SORCS2 in DSB formation. We show that SorCS2 loss is associated with elevated DSB levels in the mouse dentate gyrus and that knocking out SORCS2 in a human neuronal cell line increased Topoisomerase IIβ-dependent DSB formation and reduced neuronal viability. Neuronal stimulation had no impact on levels of DNA breaks in vitro, suggesting that the observed differences may not be the result of aberrant neuronal activity in these cells. Our findings are consistent with studies linking the VPS10 receptors and DNA damage to neurodegenerative conditions.

Opioid medication use and blood DNA methylation: epigenome-wide association meta-analysis

Aim: To identify differential methylation related to prescribed opioid use. Methods: This study examined whether blood DNA methylation, measured using Illumina arrays, differs by recent opioid medication use in four population-based cohorts. We meta-analyzed results (282 users; 10,560 nonusers) using inverse-variance weighting. Results: Differential methylation (false discovery rate <0.05) was observed at six CpGs annotated to the following genes: KIAA0226, CPLX2, TDRP, RNF38, TTC23 and GPR179. Integrative epigenomic analyses linked implicated loci to regulatory elements in blood and/or brain. Additionally, 74 CpGs were differentially methylated in males or females. Methylation at significant CpGs correlated with gene expression in blood and/or brain. Conclusion: This study identified DNA methylation related to opioid medication use in general populations. The results could inform the development of blood methylation biomarkers of opioid use.

Local CpG density affects the trajectory and variance of age-associated DNA methylation changes

BACKGROUND

DNA methylation is an epigenetic mark associated with the repression of gene promoters. Its pattern in the genome is disrupted with age and these changes can be used to statistically predict age with epigenetic clocks. Altered rates of aging inferred from these clocks are observed in human disease. However, the molecular mechanisms underpinning age-associated DNA methylation changes remain unknown. Local DNA sequence can program steady-state DNA methylation levels, but how it influences age-associated methylation changes is unknown.

RESULTS

We analyze longitudinal human DNA methylation trajectories at 345,895 CpGs from 600 individuals aged between 67 and 80 to understand the factors responsible for age-associated epigenetic changes at individual CpGs. We show that changes in methylation with age occur at 182,760 loci largely independently of variation in cell type proportions. These changes are especially apparent at 8322 low CpG density loci. Using SNP data from the same individuals, we demonstrate that methylation trajectories are affected by local sequence polymorphisms at 1487 low CpG density loci. More generally, we find that low CpG density regions are particularly prone to change and do so variably between individuals in people aged over 65. This differs from the behavior of these regions in younger individuals where they predominantly lose methylation.

CONCLUSIONS

Our results, which we reproduce in two independent groups of individuals, demonstrate that local DNA sequence influences age-associated DNA methylation changes in humans in vivo. We suggest that this occurs because interactions between CpGs reinforce maintenance of methylation patterns in CpG dense regions.

Alcohol use disorder is associated with DNA methylation-based shortening of telomere length and regulated by TESPA1: implications for aging

Chronic heavy alcohol consumption is associated with increased mortality and morbidity and often leads to premature aging; however, the mechanisms of alcohol-associated cellular aging are not well understood. In this study, we used DNA methylation derived telomere length (DNAmTL) as a novel approach to investigate the role of alcohol use on the aging process. DNAmTL was estimated by 140 cytosine phosphate guanines (CpG) sites in 372 individuals with alcohol use disorder (AUD) and 243 healthy controls (HC) and assessed using various endophenotypes and clinical biomarkers. Validation in an independent sample of DNAmTL on alcohol consumption was performed (N = 4219). Exploratory genome-wide association studies (GWAS) on DNAmTL were also performed to identify genetic variants contributing to DNAmTL shortening. Top GWAS findings were analyzed using in-silico expression quantitative trait loci analyses and related to structural MRI hippocampus volumes of individuals with AUD. DNAmTL was 0.11-kilobases shorter per year in AUD compared to HC after adjustment for age, sex, race, and blood cell composition (p = 4.0 × 10-12). This association was partially attenuated but remained significant after additionally adjusting for BMI, and smoking status (0.06 kilobases shorter per year, p = 0.002). DNAmTL shortening was strongly associated with chronic heavy alcohol use (ps < 0.001), elevated gamma-glutamyl transferase (GGT), and aspartate aminotransferase (AST) (ps < 0.004). Comparison of DNAmTL with PCR-based methods of assessing TL revealed positive correlations (R = 0.3, p = 2.2 × 10-5), highlighting the accuracy of DNAmTL as a biomarker. The GWAS meta-analysis identified a single nucleotide polymorphism (SNP), rs4374022 and 18 imputed ones in Thymocyte Expressed, Positive Selection Associated 1(TESPA1), at the genome-wide level (p = 3.75 × 10-8). The allele C of rs4374022 was associated with DNAmTL shortening, lower hippocampus volume (p < 0.01), and decreased mRNA expression in hippocampus tissue (p = 0.04). Our study demonstrates DNAmTL-related aging acceleration in AUD and suggests a functional role for TESPA1 in regulating DNAmTL length, possibly via the immune system with subsequent biological effects on brain regions negatively affected by alcohol and implicated in aging.

Identification of influential probe types in epigenetic predictions of human traits: implications for microarray design

BACKGROUND

CpG methylation levels can help to explain inter-individual differences in phenotypic traits. Few studies have explored whether identifying probe subsets based on their biological and statistical properties can maximise predictions whilst minimising array content. Variance component analyses and penalised regression (epigenetic predictors) were used to test the influence of (i) the number of probes considered, (ii) mean probe variability and (iii) methylation QTL status on the variance captured in eighteen traits by blood DNA methylation. Training and test samples comprised ≤ 4450 and ≤ 2578 unrelated individuals from Generation Scotland, respectively.

RESULTS

As the number of probes under consideration decreased, so too did the estimates from variance components and prediction analyses. Methylation QTL status and mean probe variability did not influence variance components. However, relative effect sizes were 15% larger for epigenetic predictors based on probes with known or reported methylation QTLs compared to probes without reported methylation QTLs. Relative effect sizes were 45% larger for predictors based on probes with mean Beta-values between 10 and 90% compared to those based on hypo- or hypermethylated probes (Beta-value ≤ 10% or ≥ 90%).

CONCLUSIONS

Arrays with fewer probes could reduce costs, leading to increased sample sizes for analyses. Our results show that reducing array content can restrict prediction metrics and careful attention must be given to the biological and distribution properties of CpG probes in array content selection.

Integrated methylome and phenome study of the circulating proteome reveals markers pertinent to brain health

Characterising associations between the methylome, proteome and phenome may provide insight into biological pathways governing brain health. Here, we report an integrated DNA methylation and phenotypic study of the circulating proteome in relation to brain health. Methylome-wide association studies of 4058 plasma proteins are performed (N = 774), identifying 2928 CpG-protein associations after adjustment for multiple testing. These are independent of known genetic protein quantitative trait loci (pQTLs) and common lifestyle effects. Phenome-wide association studies of each protein are then performed in relation to 15 neurological traits (N = 1,065), identifying 405 associations between the levels of 191 proteins and cognitive scores, brain imaging measures or APOE e4 status. We uncover 35 previously unreported DNA methylation signatures for 17 protein markers of brain health. The epigenetic and proteomic markers we identify are pertinent to understanding and stratifying brain health.

Pulmonary Function and Blood DNA Methylation: A Multiancestry Epigenome-Wide Association Meta-analysis

Rationale: Methylation integrates factors present at birth and modifiable across the lifespan that can influence pulmonary function. Studies are limited in scope and replication. Objectives: To conduct large-scale epigenome-wide meta-analyses of blood DNA methylation and pulmonary function. Methods: Twelve cohorts analyzed associations of methylation at cytosine-phosphate-guanine probes (CpGs), using Illumina 450K or EPIC/850K arrays, with FEV₁, FVC, and FEV₁/FVC. We performed multiancestry epigenome-wide meta-analyses (total of 17,503 individuals; 14,761 European, 2,549 African, and 193 Hispanic/Latino ancestries) and interpreted results using integrative epigenomics. Measurements and Main Results: We identified 1,267 CpGs (1,042 genes) differentially methylated (false discovery rate, <0.025) in relation to FEV₁, FVC, or FEV₁/FVC, including 1,240 novel and 73 also related to chronic obstructive pulmonary disease (1,787 cases). We found 294 CpGs unique to European or African ancestry and 395 CpGs unique to never or ever smokers. The majority of significant CpGs correlated with nearby gene expression in blood. Findings were enriched in key regulatory elements for gene function, including accessible chromatin elements, in both blood and lung. Sixty-nine implicated genes are targets of investigational or approved drugs. One example novel gene highlighted by integrative epigenomic and druggable target analysis is TNFRSF4. Mendelian randomization and colocalization analyses suggest that epigenome-wide association study signals capture causal regulatory genomic loci. Conclusions: We identified numerous novel loci differentially methylated in relation to pulmonary function; few were detected in large genome-wide association studies. Integrative analyses highlight functional relevance and potential therapeutic targets. This comprehensive discovery of potentially modifiable, novel lung function loci expands knowledge gained from genetic studies, providing insights into lung pathogenesis.

Complex trait methylation scores in the prediction of major depressive disorder

BACKGROUND

DNA methylation (DNAm) is associated with time-varying environmental factors that contribute to major depressive disorder (MDD) risk. We sought to test whether DNAm signatures of lifestyle and biochemical factors were associated with MDD to reveal dynamic biomarkers of MDD risk that may be amenable to lifestyle interventions.

METHODS

Here, we calculated methylation scores (MS) at multiple p-value thresholds for lifestyle (BMI, smoking, alcohol consumption, and educational attainment) and biochemical (high-density lipoprotein (HDL) and total cholesterol) factors in Generation Scotland (GS) (N=9,502) and in a replication cohort (ALSPAC_adults, N=565), using CpG sites reported in previous well-powered methylome-wide association studies. We also compared their predictive accuracy for MDD to a MDD MS in an independent GS sub-sample (N=4,432).

FINDINGS

Each trait MS was significantly associated with its corresponding phenotype in GS (β_range=0.089-1.457) and in ALSPAC (β_range=0.078-2.533). Each MS was also significantly associated with MDD before and after adjustment for its corresponding phenotype in GS (β_range=0.053-0.145). After accounting for relevant lifestyle factors, MS for educational attainment (β=0.094) and alcohol consumption (MS_p-value<0.01-0.5; β_range=-0.069-0.083) remained significantly associated with MDD in GS. Smoking (AUC=0.569) and educational attainment (AUC=0.585) MSs could discriminate MDD from controls better than the MDD MS (AUC=0.553) in the independent GS sub-sample. Analyses implicating MDD did not replicate across ALSPAC, although the direction of effect was consistent for all traits when adjusting for the MS corresponding phenotypes.

INTERPRETATION

We showed that lifestyle and biochemical MS were associated with MDD before and after adjustment for their corresponding phenotypes (p_nominal<0.05), but not when smoking, alcohol consumption, and BMI were also included as covariates. MDD results did not replicate in the smaller, female-only independent ALSPAC cohort (N_ALSPAC=565; N_GS=9,502), potentially due to demographic differences or low statistical power, but effect sizes were consistent with the direction reported in GS. DNAm scores for modifiable MDD risk factors may contribute to disease vulnerability and, in some cases, explain additional variance to their observed phenotypes.

FUNDING

Wellcome Trust.

Genome- and epigenome-wide studies of plasma protein biomarkers for Alzheimer's disease implicate TBCA and TREM2 in disease risk

Introduction

The levels of many blood proteins are associated with Alzheimer's disease (AD) or its pathological hallmarks. Elucidating the molecular factors that control circulating levels of these proteins may help to identify proteins associated with disease risk mechanisms.

Methods

Genome-wide and epigenome-wide studies (nindividuals ≤1064) were performed on plasma levels of 282 AD-associated proteins, identified by a structured literature review. Bayesian penalized regression estimated contributions of genetic and epigenetic variation toward inter-individual differences in plasma protein levels. Mendelian randomization (MR) and co-localization tested associations between proteins and disease-related phenotypes.

Results

Sixty-four independent genetic and 26 epigenetic loci were associated with 45 proteins. Novel findings included an association between plasma triggering receptor expressed on myeloid cells 2 (TREM2) levels and a polymorphism and cytosine-phosphate-guanine (CpG) site within the MS4A4A locus. Higher plasma tubulin-specific chaperone A (TBCA) and TREM2 levels were significantly associated with lower AD risk.

Discussion

Our data inform the regulation of biomarker levels and their relationships with AD.

DNA methylome-wide association study of genetic risk for depression implicates antigen processing and immune responses

BACKGROUND

Depression is a disabling and highly prevalent condition where genetic and epigenetic, such as DNA methylation (DNAm), differences contribute to disease risk. DNA methylation is influenced by genetic variation but the association between polygenic risk of depression and DNA methylation is unknown.

METHODS

We investigated the association between polygenic risk scores (PRS) for depression and DNAm by conducting a methylome-wide association study (MWAS) in Generation Scotland (N = 8898, mean age = 49.8 years) with replication in the Lothian Birth Cohorts of 1921 and 1936 and adults in the Avon Longitudinal Study of Parents and Children (ALSPAC) (Ncombined = 2049, mean age = 79.1, 69.6 and 47.2 years, respectively). We also conducted a replication MWAS in the ALSPAC children (N = 423, mean age = 17.1 years). Gene ontology analysis was conducted for the cytosine-guanine dinucleotide (CpG) probes significantly associated with depression PRS, followed by Mendelian randomisation (MR) analysis to infer the causal relationship between depression and DNAm.

RESULTS

Widespread associations (NCpG = 71, pBonferroni < 0.05, p < 6.3 × 10-8) were found between PRS constructed using genetic risk variants for depression and DNAm in CpG probes that localised to genes involved in immune responses and neural development. The effect sizes for the significant associations were highly correlated between the discovery and replication samples in adults (r = 0.79) and in adolescents (r = 0.82). Gene Ontology analysis showed that significant CpG probes are enriched in immunological processes in the human leukocyte antigen system. Additional MWAS was conducted for each lead genetic risk variant. Over 47.9% of the independent genetic risk variants included in the PRS showed associations with DNAm in CpG probes located in both the same (cis) and distal (trans) locations to the genetic loci (pBonferroni < 0.045). Subsequent MR analysis showed that there are a greater number of causal effects found from DNAm to depression than vice versa (DNAm to depression: pFDR ranged from 0.024 to 7.45 × 10-30; depression to DNAm: pFDR ranged from 0.028 to 0.003).

CONCLUSIONS

PRS for depression, especially those constructed from genome-wide significant genetic risk variants, showed methylome-wide differences associated with immune responses. Findings from MR analysis provided evidence for causal effect of DNAm to depression.

Cohort profile for the STratifying Resilience and Depression Longitudinally (STRADL) study: A depression-focused investigation of Generation Scotland, using detailed clinical, cognitive, and neuroimaging assessments

STratifying Resilience and Depression Longitudinally (STRADL) is a population-based study built on the Generation Scotland: Scottish Family Health Study (GS:SFHS) resource. The aim of STRADL is to subtype major depressive disorder (MDD) on the basis of its aetiology, using detailed clinical, cognitive, and brain imaging assessments. The GS:SFHS provides an important opportunity to study complex gene-environment interactions, incorporating linkage to existing datasets and inclusion of early-life variables for two longitudinal birth cohorts. Specifically, data collection in STRADL included: socio-economic and lifestyle variables; physical measures; questionnaire data that assesses resilience, early-life adversity, personality, psychological health, and lifetime history of mood disorder; laboratory samples; cognitive tests; and brain magnetic resonance imaging. Some of the questionnaire and cognitive data were first assessed at the GS:SFHS baseline assessment between 2006-2011, thus providing longitudinal measures relevant to the study of depression, psychological resilience, and cognition. In addition, routinely collected historic NHS data and early-life variables are linked to STRADL data, further providing opportunities for longitudinal analysis. Recruitment has been completed and we consented and tested 1,188 participants.

Epigenome-wide association study of alcohol consumption in N = 8161 individuals and relevance to alcohol use disorder pathophysiology: identification of the cystine/glutamate transporter SLC7A11 as a top target

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

Methylome-wide association study of antidepressant use in Generation Scotland and the Netherlands Twin Register implicates the innate immune system

Antidepressants are an effective treatment for major depressive disorder (MDD), although individual response is unpredictable and highly variable. Whilst the mode of action of antidepressants is incompletely understood, many medications are associated with changes in DNA methylation in genes that are plausibly linked to their mechanisms. Studies of DNA methylation may therefore reveal the biological processes underpinning the efficacy and side effects of antidepressants. We performed a methylome-wide association study (MWAS) of self-reported antidepressant use accounting for lifestyle factors and MDD in Generation Scotland (GS:SFHS, N = 6428, EPIC array) and the Netherlands Twin Register (NTR, N = 2449, 450 K array) and ran a meta-analysis of antidepressant use across these two cohorts. We found ten CpG sites significantly associated with self-reported antidepressant use in GS:SFHS, with the top CpG located within a gene previously associated with mental health disorders, ATP6V1B2 (β = -0.055, pcorrected = 0.005). Other top loci were annotated to genes including CASP10, TMBIM1, MAPKAPK3, and HEBP2, which have previously been implicated in the innate immune response. Next, using penalised regression, we trained a methylation-based score of self-reported antidepressant use in a subset of 3799 GS:SFHS individuals that predicted antidepressant use in a second subset of GS:SFHS (N = 3360, β = 0.377, p = 3.12 × 10-11, R2 = 2.12%). In an MWAS analysis of prescribed selective serotonin reuptake inhibitors, we showed convergent findings with those based on self-report. In NTR, we did not find any CpGs significantly associated with antidepressant use. The meta-analysis identified the two CpGs of the ten above that were common to the two arrays used as being significantly associated with antidepressant use, although the effect was in the opposite direction for one of them. Antidepressants were associated with epigenetic alterations in loci previously associated with mental health disorders and the innate immune system. These changes predicted self-reported antidepressant use in a subset of GS:SFHS and identified processes that may be relevant to our mechanistic understanding of clinically relevant antidepressant drug actions and side effects.

Blood-based epigenome-wide analyses of cognitive abilities

BACKGROUND

Blood-based markers of cognitive functioning might provide an accessible way to track neurodegeneration years prior to clinical manifestation of cognitive impairment and dementia.

RESULTS

Using blood-based epigenome-wide analyses of general cognitive function, we show that individual differences in DNA methylation (DNAm) explain 35.0% of the variance in general cognitive function (g). A DNAm predictor explains ~4% of the variance, independently of a polygenic score, in two external cohorts. It also associates with circulating levels of neurology- and inflammation-related proteins, global brain imaging metrics, and regional cortical volumes.

CONCLUSIONS

As sample sizes increase, the ability to assess cognitive function from DNAm data may be informative in settings where cognitive testing is unreliable or unavailable.

Epigenetic scores for the circulating proteome as tools for disease prediction

Protein biomarkers have been identified across many age-related morbidities. However, characterising epigenetic influences could further inform disease predictions. Here, we leverage epigenome-wide data to study links between the DNA methylation (DNAm) signatures of the circulating proteome and incident diseases. Using data from four cohorts, we trained and tested epigenetic scores (EpiScores) for 953 plasma proteins, identifying 109 scores that explained between 1% and 58% of the variance in protein levels after adjusting for known protein quantitative trait loci (pQTL) genetic effects. By projecting these EpiScores into an independent sample (Generation Scotland; n = 9537) and relating them to incident morbidities over a follow-up of 14 years, we uncovered 137 EpiScore-disease associations. These associations were largely independent of immune cell proportions, common lifestyle and health factors, and biological aging. Notably, we found that our diabetes-associated EpiScores highlighted previous top biomarker associations from proteome-wide assessments of diabetes. These EpiScores for protein levels can therefore be a valuable resource for disease prediction and risk stratification.

Associations between alcohol use and accelerated biological ageing

Harmful alcohol use is a leading cause of premature death and is associated with age-related disease. Biological ageing is highly variable between individuals and may deviate from chronological ageing, suggesting that biomarkers of biological ageing (derived from DNA methylation or brain structural measures) may be clinically relevant. Here, we investigated the relationships between alcohol phenotypes and both brain and DNA methylation age estimates. First, using data from UK Biobank and Generation Scotland, we tested the association between alcohol consumption (units/week) or hazardous use (Alcohol Use Disorders Identification Test [AUDIT] scores) and accelerated brain and epigenetic ageing in 20,258 and 8051 individuals, respectively. Second, we used Mendelian randomisation (MR) to test for a causal effect of alcohol consumption levels and alcohol use disorder (AUD) on biological ageing. Alcohol use showed a consistent positive association with higher predicted brain age (AUDIT-C: β = 0.053, p = 3.16 × 10-13 ; AUDIT-P: β = 0.052, p = 1.6 × 10-13 ; total AUDIT score: β = 0.062, p = 5.52 × 10-16 ; units/week: β = 0.078, p = 2.20 × 10-16 ), and two DNA methylation-based estimates of ageing, GrimAge (units/week: β = 0.053, p = 1.48 × 10-7 ) and PhenoAge (units/week: β = 0.077, p = 2.18x10-10 ). MR analyses revealed limited evidence for a causal effect of AUD on accelerated brain ageing (β = 0.118, p = 0.044). However, this result should be interpreted cautiously as the significant effect was driven by a single genetic variant. We found no evidence for a causal effect of alcohol consumption levels on accelerated biological ageing. Future studies investigating the mechanisms associating alcohol use with accelerated biological ageing are warranted.

Meta-analyses identify DNA methylation associated with kidney function and damage

Chronic kidney disease is a major public health burden. Elevated urinary albumin-to-creatinine ratio is a measure of kidney damage, and used to diagnose and stage chronic kidney disease. To extend the knowledge on regulatory mechanisms related to kidney function and disease, we conducted a blood-based epigenome-wide association study for estimated glomerular filtration rate (n = 33,605) and urinary albumin-to-creatinine ratio (n = 15,068) and detected 69 and seven CpG sites where DNA methylation was associated with the respective trait. The majority of these findings showed directionally consistent associations with the respective clinical outcomes chronic kidney disease and moderately increased albuminuria. Associations of DNA methylation with kidney function, such as CpGs at JAZF1, PELI1 and CHD2 were validated in kidney tissue. Methylation at PHRF1, LDB2, CSRNP1 and IRF5 indicated causal effects on kidney function. Enrichment analyses revealed pathways related to hemostasis and blood cell migration for estimated glomerular filtration rate, and immune cell activation and response for urinary albumin-to-creatinineratio-associated CpGs.

Lifestyle and Genetic Factors Modify Parent-of-Origin Effects on the Human Methylome

BACKGROUND

parent-of-origin effects (POE) play important roles in complex disease and thus understanding their regulation and associated molecular and phenotypic variation are warranted. Previous studies mainly focused on the detection of genomic regions or phenotypes regulated by POE. Understanding whether POE may be modified by environmental or genetic exposures is important for understanding of the source of POE-associated variation, but only a few case studies addressing modifiable POE exist.

METHODS

in order to understand this high order of POE regulation, we screened 101 genetic and environmental factors such as 'predicted mRNA expression levels' of DNA methylation/imprinting machinery genes and environmental exposures. POE-mQTL-modifier interaction models were proposed to test the potential of these factors to modify POE at DNA methylation using data from Generation Scotland: The Scottish Family Health Study(N=2315).

FINDINGS

a set of vulnerable/modifiable POE-CpGs were identified (modifiable-POE-regulated CpGs, N=3). Four factors, 'lifetime smoking status' and 'predicted mRNA expression levels' of TET2, SIRT1 and KDM1A, were found to significantly modify the POE on the three CpGs in both discovery and replication datasets. We further identified plasma protein and health-related phenotypes associated with the methylation level of one of the identified CpGs.

INTERPRETATION

the modifiable POE identified here revealed an important yet indirect path through which genetic background and environmental exposures introduce their effect on DNA methylation, motivating future comprehensive evaluation of the role of these modifiers in complex diseases.

FUNDING

NSFC (81971270),H2020-MSCA-ITN(721815), Wellcome (204979/Z/16/Z,104036/Z/14/Z), MRC (MC_UU_00007/10, MC_PC_U127592696), CSO (CZD/16/6,CZB/4/276, CZB/4/710), SFC (HR03006), EUROSPAN (LSHG-CT-2006-018947), BBSRC (BBS/E/D/30002276), SYSU, Arthritis Research UK, NHLBI, NIH.

Methylome-wide association study of early life stressors and adult mental health

The environment and events that we are exposed to in utero, during birth and in early childhood influence our future physical and mental health. The underlying mechanisms that lead to these outcomes are unclear, but long-term changes in epigenetic marks, such as DNA methylation, could act as a mediating factor or biomarker. DNA methylation data were assayed at 713 522 CpG sites from 9537 participants of the Generation Scotland: Scottish Family Health Study, a family-based cohort with extensive genetic, medical, family history and lifestyle information. Methylome-wide association studies of eight early life environment phenotypes and two adult mental health phenotypes (major depressive disorder and brief resilience scale) were conducted using DNA methylation data collected from adult whole blood samples. Two genes involved with different developmental pathways (PRICKLE2, Prickle Planar Cell Polarity Protein 2 and ABI1, Abl-Interactor-1) were annotated to CpG sites associated with preterm birth (P < 1.27 × 10⁻⁹). A further two genes important to the development of sensory pathways (SOBP, Sine Oculis Binding Protein Homolog and RPGRIP1, Retinitis Pigmentosa GTPase Regulator Interacting Protein) were annotated to sites associated with low birth weight (P < 4.35 × 10⁻⁸). The examination of methylation profile scores and genes and gene-sets annotated from associated CpGs sites found no evidence of overlap between the early life environment and mental health conditions. Birth date was associated with a significant difference in estimated lymphocyte and neutrophil counts. Previous studies have shown that early life environments influence the risk of developing mental health disorders later in life; however, this study found no evidence that this is mediated by stable changes to the methylome detectable in peripheral blood.

Genome-wide methylation data improves dissection of the effect of smoking on body mass index

Variation in obesity-related traits has a genetic basis with heritabilities between 40 and 70%. While the global obesity pandemic is usually associated with environmental changes related to lifestyle and socioeconomic changes, most genetic studies do not include all relevant environmental covariates, so the genetic contribution to variation in obesity-related traits cannot be accurately assessed. Some studies have described interactions between a few individual genes linked to obesity and environmental variables but there is no agreement on their total contribution to differences between individuals. Here we compared self-reported smoking data and a methylation-based proxy to explore the effect of smoking and genome-by-smoking interactions on obesity related traits from a genome-wide perspective to estimate the amount of variance they explain. Our results indicate that exploiting omic measures can improve models for complex traits such as obesity and can be used as a substitute for, or jointly with, environmental records to better understand causes of disease.

Epigenetic prediction of major depressive disorder

Variation in DNA methylation (DNAm) is associated with lifestyle factors such as smoking and body mass index (BMI) but there has been little research exploring its ability to identify individuals with major depressive disorder (MDD). Using penalised regression on genome-wide CpG methylation, we tested whether DNAm risk scores (MRS), trained on 1223 MDD cases and 1824 controls, could discriminate between cases (n = 363) and controls (n = 1417) in an independent sample, comparing their predictive accuracy to polygenic risk scores (PRS). The MRS explained 1.75% of the variance in MDD (β = 0.338, p = 1.17 × 10-7) and remained associated after adjustment for lifestyle factors (β = 0.219, p = 0.001, R2 = 0.68%). When modelled alongside PRS (β = 0.384, p = 4.69 × 10-9) the MRS remained associated with MDD (β = 0.327, p = 5.66 × 10-7). The MRS was also associated with incident cases of MDD who were well at recruitment but went on to develop MDD at a later assessment (β = 0.193, p = 0.016, R2 = 0.52%). Heritability analyses found additive genetic effects explained 22% of variance in the MRS, with a further 19% explained by pedigree-associated genetic effects and 16% by the shared couple environment. Smoking status was also strongly associated with MRS (β = 0.440, p ≤ 2 × 10-16). After removing smokers from the training set, the MRS strongly associated with BMI (β = 0.053, p = 0.021). We tested the association of MRS with 61 behavioural phenotypes and found that whilst PRS were associated with psychosocial and mental health phenotypes, MRS were more strongly associated with lifestyle and sociodemographic factors. DNAm-based risk scores of MDD significantly discriminated MDD cases from controls in an independent dataset and may represent an archive of exposures to lifestyle factors that are relevant to the prediction of MDD.

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.

An epigenetic predictor of death captures multi-modal measures of brain health

Individuals of the same chronological age exhibit disparate rates of biological ageing. Consequently, a number of methodologies have been proposed to determine biological age and primarily exploit variation at the level of DNA methylation (DNAm). A novel epigenetic clock, termed 'DNAm GrimAge' has outperformed its predecessors in predicting the risk of mortality as well as many age-related morbidities. However, the association between DNAm GrimAge and cognitive or neuroimaging phenotypes remains unknown. We explore these associations in the Lothian Birth Cohort 1936 (n = 709, mean age 73 years). Higher DNAm GrimAge was strongly associated with all-cause mortality over the eighth decade (Hazard Ratio per standard deviation increase in GrimAge: 1.81, P < 2.0 × 10-16). Higher DNAm GrimAge was associated with lower age 11 IQ (β = -0.11), lower age 73 general cognitive ability (β = -0.18), decreased brain volume (β = -0.25) and increased brain white matter hyperintensities (β = 0.17). There was tentative evidence for a longitudinal association between DNAm GrimAge and cognitive decline from age 70 to 79. Sixty-nine of 137 health- and brain-related phenotypes tested were significantly associated with GrimAge. Adjusting all models for childhood intelligence attenuated to non-significance a small number of associations (12/69 associations; 6 of which were cognitive traits), but not the association with general cognitive ability (33.9% attenuation). Higher DNAm GrimAge associates with lower cognitive ability and brain vascular lesions in older age, independently of early-life cognitive ability. This epigenetic predictor of mortality associates with different measures of brain health and may aid in the prediction of age-related cognitive decline.

Cohort profile for the STratifying Resilience and Depression Longitudinally (STRADL) study: A depression-focused investigation of Generation Scotland, using detailed clinical, cognitive, and neuroimaging assessments

STratifying Resilience and Depression Longitudinally (STRADL) is a population-based study built on the Generation Scotland: Scottish Family Health Study (GS:SFHS) resource. The aim of STRADL is to subtype major depressive disorder (MDD) on the basis of its aetiology, using detailed clinical, cognitive, and brain imaging assessments. The GS:SFHS provides an important opportunity to study complex gene-environment interactions, incorporating linkage to existing datasets and inclusion of early-life variables for two longitudinal birth cohorts. Specifically, data collection in STRADL included: socio-economic and lifestyle variables; physical measures; questionnaire data that assesses resilience, early-life adversity, personality, psychological health, and lifetime history of mood disorder; laboratory samples; cognitive tests; and brain magnetic resonance imaging. Some of the questionnaire and cognitive data were first assessed at the GS:SFHS baseline assessment between 2006-2011, thus providing longitudinal measures relevant to the study of depression, psychological resilience, and cognition. In addition, routinely collected historic NHS data and early-life variables are linked to STRADL data, further providing opportunities for longitudinal analysis. Recruitment has been completed and we consented and tested 1,188 participants.

Epigenome-wide association study and multi-tissue replication of individuals with alcohol use disorder: evidence for abnormal glucocorticoid signaling pathway gene regulation

Alcohol use disorder (AUD) is a chronic debilitating disorder with limited treatment options and poorly defined pathophysiology. There are substantial genetic and epigenetic components; however, the underlying mechanisms contributing to AUD remain largely unknown. We conducted the largest DNA methylation epigenome-wide association study (EWAS) analyses currently available for AUD (total N = 625) and employed a top hit replication (N = 4798) using a cross-tissue/cross-phenotypic approach with the goal of identifying novel epigenetic targets relevant to AUD. Results show that a network of differentially methylated regions in glucocorticoid signaling and inflammation-related genes were associated with alcohol use behaviors. A top probe consistently associated across all cohorts was located in the long non-coding RNA growth arrest specific five gene (GAS5) (p < 10-24). GAS5 has been implicated in regulating transcriptional activity of the glucocorticoid receptor and has multiple functions related to apoptosis, immune function and various cancers. Endophenotypic analyses using peripheral cortisol levels and neuroimaging paradigms showed that methylomic variation in GAS5 network-related probes were associated with stress phenotypes. Postmortem brain analyses documented increased GAS5 expression in the amygdala of individuals with AUD. Our data suggest that alcohol use is associated with differential methylation in the glucocorticoid system that might influence stress and inflammatory reactivity and subsequently risk for AUD.

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.

Birth weight associations with DNA methylation differences in an adult population

The Developmental Origins of Health and Disease (DOHaD) theory predicts that prenatal and early life events shape adult health outcomes. Birth weight is a useful indicator of the foetal experience and has been associated with multiple adult health outcomes. DNA methylation (DNAm) is one plausible mechanism behind the relationship of birth weight to adult health. Through data linkage between Generation Scotland and historic Scottish birth cohorts, and birth records held through the NHS Information and Statistics Division, a sample of 1,757 individuals with available birth weight and DNAm data was derived. Epigenome-wide association studies (EWAS) were performed in two independently generated DNAm subgroups (nSet1 = 1,395, nSet2 = 362), relating adult DNAm from whole blood to birth weight. Meta-analysis yielded one genome-wide significant CpG site (p = 5.97x10-9), cg00966482. There was minimal evidence for attenuation of the effect sizes for the lead loci upon adjustment for numerous potential confounder variables (body mass index, educational attainment, and socioeconomic status). Associations between birth weight and epigenetic measures of biological age were also assessed. Associations between lower birth weight and higher Grim Age acceleration (p(FDR) = 3.6x10-3) and shorter DNAm-derived telomere length (p(FDR) = 1.7x10-3) are described, although results for three other epigenetic clocks were null. Our results provide support for an association between birth weight and DNAm both locally at one CpG site, and globally via biological ageing estimates.

Epigenome-wide association study of kidney function identifies trans-ethnic and ethnic-specific loci

BACKGROUND

DNA methylation (DNAm) is associated with gene regulation and estimated glomerular filtration rate (eGFR), a measure of kidney function. Decreased eGFR is more common among US Hispanics and African Americans. The causes for this are poorly understood. We aimed to identify trans-ethnic and ethnic-specific differentially methylated positions (DMPs) associated with eGFR using an agnostic, genome-wide approach.

METHODS

The study included up to 5428 participants from multi-ethnic studies for discovery and 8109 participants for replication. We tested the associations between whole blood DNAm and eGFR using beta values from Illumina 450K or EPIC arrays. Ethnicity-stratified analyses were performed using linear mixed models adjusting for age, sex, smoking, and study-specific and technical variables. Summary results were meta-analyzed within and across ethnicities. Findings were assessed using integrative epigenomics methods and pathway analyses.

RESULTS

We identified 93 DMPs associated with eGFR at an FDR of 0.05 and replicated 13 and 1 DMPs across independent samples in trans-ethnic and African American meta-analyses, respectively. The study also validated 6 previously published DMPs. Identified DMPs showed significant overlap enrichment with DNase I hypersensitive sites in kidney tissue, sites associated with the expression of proximal genes, and transcription factor motifs and pathways associated with kidney tissue and kidney development.

CONCLUSIONS

We uncovered trans-ethnic and ethnic-specific DMPs associated with eGFR, including DMPs enriched in regulatory elements in kidney tissue and pathways related to kidney development. These findings shed light on epigenetic mechanisms associated with kidney function, bridging the gap between population-specific eGFR-associated DNAm and tissue-specific regulatory context.

Creating and validating a DNA methylation-based proxy for interleukin-6

Background

Studies evaluating the relationship between chronic inflammation and cognitive functioning have produced heterogeneous results. A potential reason for this is the variability of inflammatory mediators which could lead to misclassifications of individuals’ persisting levels of inflammation. DNA methylation (DNAm) has shown utility in indexing environmental exposures and could be leveraged to provide proxy signatures of chronic inflammation.

Method

We conducted an elastic net regression of interleukin-6 (IL-6) in a cohort of 875 older adults (Lothian Birth Cohort 1936; mean age: 70 years) to develop a DNAm-based predictor. The predictor was tested in an independent cohort (Generation Scotland; N = 7028 [417 with measured IL-6], mean age: 51 years).

Results

A weighted score from 35 CpG sites optimally predicted IL-6 in the independent test set (Generation Scotland; R² = 4.4%, p = 2.1 × 10⁻⁵). In the independent test cohort, both measured IL-6 and the DNAm proxy increased with age (serum IL-6: n = 417, β = 0.02, SE = 0.004, p = 1.3 × 10⁻⁷; DNAm IL-6 score: N = 7028, β = 0.02, SE = 0.0009, p < 2 × 10⁻¹⁶). Serum IL-6 did not associate with cognitive ability (n = 417, β = −0.06, SE = 0.05, p = .19); however, an inverse association was identified between the DNAm score and cognitive functioning (N = 7028, β = −0.16, SE = 0.02, p_FDR < 2 × 10⁻¹⁶).

Conclusions

These results suggest methylation-based predictors can be used as proxies for inflammatory markers, potentially allowing for further insight into the relationship between inflammation and pertinent health outcomes.

Structural brain correlates of serum and epigenetic markers of inflammation in major depressive disorder

Inflammatory processes are implicated in the aetiology of Major Depressive Disorder (MDD); however, the relationship between peripheral inflammation, brain structure and depression remains unclear, partly due to complexities around the use of acute/phasic inflammatory biomarkers. Here, we report the first large-scale study of both serological and methylomic signatures of CRP (considered to represent acute and chronic measures of inflammation respectively) and their associations with depression status/symptoms, and structural neuroimaging phenotypes (T1 and diffusion MRI) in a large community-based sample (Generation Scotland; NMDD cases = 271, Ncontrols = 609). Serum CRP was associated with overall MDD severity, and specifically with current somatic symptoms- general interest (β = 0.145, PFDR = 6 × 10-4) and energy levels (β = 0.101, PFDR = 0.027), along with reduced entorhinal cortex thickness (β = -0.095, PFDR = 0.037). DNAm CRP was significantly associated with reduced global grey matter/cortical volume and widespread reductions in integrity of 16/24 white matter tracts (with greatest regional effects in the external and internal capsules, βFA= -0.12 to -0.14). In general, the methylation-based measures showed stronger associations with imaging metrics than serum-based CRP measures (βaverage = -0.15 versus βaverage = 0.01 respectively). These findings provide evidence for central effects of peripheral inflammation from both serological and epigenetic markers of inflammation, including in brain regions previously implicated in depression. This suggests that these imaging measures may be involved in the relationship between peripheral inflammation and somatic/depressive symptoms. Notably, greater effects on brain morphology were seen for methylation-based rather than serum-based measures of inflammation, indicating the importance of such measures for future studies.

Identification of epigenome-wide DNA methylation differences between carriers of APOE ε4 and APOE ε2 alleles

BACKGROUND

The apolipoprotein E (APOE) ε4 allele is the strongest genetic risk factor for late onset Alzheimer's disease, whilst the ε2 allele confers protection. Previous studies report differential DNA methylation of APOE between ε4 and ε2 carriers, but associations with epigenome-wide methylation have not previously been characterised.

METHODS

Using the EPIC array, we investigated epigenome-wide differences in whole blood DNA methylation patterns between Alzheimer's disease-free APOE ε4 (n = 2469) and ε2 (n = 1118) carriers from the two largest single-cohort DNA methylation samples profiled to date. Using a discovery, replication and meta-analysis study design, methylation differences were identified using epigenome-wide association analysis and differentially methylated region (DMR) approaches. Results were explored using pathway and methylation quantitative trait loci (meQTL) analyses.

RESULTS

We obtained replicated evidence for DNA methylation differences in a ~ 169 kb region, which encompasses part of APOE and several upstream genes. Meta-analytic approaches identified DNA methylation differences outside of APOE: differentially methylated positions were identified in DHCR24, LDLR and ABCG1 (2.59 × 10-100 ≤ P ≤ 2.44 × 10-8) and DMRs were identified in SREBF2 and LDLR (1.63 × 10-4 ≤ P ≤ 3.01 × 10-2). Pathway and meQTL analyses implicated lipid-related processes and high-density lipoprotein cholesterol was identified as a partial mediator of the methylation differences in ABCG1 and DHCR24.

CONCLUSIONS

APOE ε4 vs. ε2 carrier status is associated with epigenome-wide methylation differences in the blood. The loci identified are located in trans as well as cis to APOE and implicate genes involved in lipid homeostasis.

Epigenetic Age Acceleration and Cognitive Function in African American Adults in Midlife: The Atherosclerosis Risk in Communities Study

Methylation levels measured at defined sites across the genome have recently been shown to be correlated with an individual's chronological age. Age acceleration, or the difference between age estimated from DNA methylation status and chronological age, has been proposed as a novel biomarker of aging. In this study, the cross-sectional association between two different measures of age acceleration and cognitive function was investigated using whole blood samples from 2,157 African American participants 47-70 years of age in the population-based Atherosclerosis Risk in Communities (ARIC) Study. Cognition was evaluated using three domain-specific tests. A significant inverse association between a 1-year increase in age acceleration calculated using a blood-based age predictor and scores on the Word Fluency Test was found using a general linear model adjusted for chronological age, gender, and years of education (β = -0.140 words; p = .001) and after adding other potential confounding variables (β = -0.104 words, p = .023). The results were replicated in 1,670 European participants in the Generation Scotland: Scottish Family Health Study (fully adjusted model: β = -0.199 words; p = .034). A significant association was also identified in a trans-ethnic meta-analysis across cohorts that included an additional 708 European American ARIC study participants (fully adjusted model: β = -0.110 words, p = .003). There were no associations found using an estimate of age acceleration derived from multiple tissues. These findings provide evidence that age acceleration is a correlate of performance on a test of verbal fluency in middle-aged adults.

Butterfly richness and abundance along a gradient of imperviousness and the importance of matrix quality

Heterogeneity in quantity and quality of resources provided in the urban matrix may mitigate adverse effects of urbanization intensity on the structure of biotic communities. To assess this we quantified the spatial variation in butterfly richness and abundance along an impervious surface gradient using three measures of urban matrix quality: floral resource availability and origin (native vs. exotic plants), tree cover, and the occurrence of remnant habitat patches. Butterfly richness and abundance were surveyed in 100 cells (500 × 500 m), selected using a random-stratified sampling design, across a continuous gradient of imperviousness in Melbourne, Australia. Sampling occurred twice during the butterfly flight season. Occurrence data were analyzed using generalized linear models at local and mesoscales. Despite high sampling completeness, we did not detect 75% of species from the regional species pool in the urban area, suggesting that urbanization has caused a large proportion of the region's butterflies to become absent or extremely rare within Melbourne's metro-area. Those species that do remain are largely very generalist in their choice of larval host plants. Butterfly species richness and abundance declined with increasing impervious surface cover and, contrary to evidence for other taxa, there was no evidence that richness peaked at intermediate levels of urbanization. Declines in abundance appeared to be more noticeable when impervious surface cover exceeded 25%, while richness declined linearly with increasing impervious surface cover. We find evidence that the quality of the urban matrix (floral resources and remnant vegetation) influenced butterfly richness and abundance although the effects were small. Total butterfly abundance responded negatively to exotic floral abundance early in the sampling season and positively to total floral abundance later in the sampling season. Butterfly species richness increased with tree cover. Negative impacts of increased urbanization intensity on butterfly species richness and abundance may be mitigated to some extent by improving the quality of the urban matrix by enhancing tree cover and the provision of floral resources, with some evidence that native plants are more effective.

Epigenome-wide analyses identify DNA methylation signatures of dementia risk

INTRODUCTION

Dementia pathogenesis begins years before clinical symptom onset, necessitating the understanding of premorbid risk mechanisms. Here we investigated potential pathogenic mechanisms by assessing DNA methylation associations with dementia risk factors in Alzheimer's disease (AD)-free participants.

METHODS

Associations between dementia risk measures (family history, AD genetic risk score [GRS], and dementia risk scores [combining lifestyle, demographic, and genetic factors]) and whole-blood DNA methylation were assessed in discovery and replication samples (n = ~400 to ~5000) from Generation Scotland.

RESULTS

AD genetic risk and two dementia risk scores were associated with differential methylation. The GRS associated predominantly with methylation differences in cis but also identified a genomic region implicated in Parkinson disease. Loci associated with dementia risk scores were enriched for those previously associated with body mass index and alcohol consumption.

DISCUSSION

Dementia risk measures show widespread association with blood-based methylation, generating several hypotheses for assessment by future studies.

Epigenetic measures of ageing predict the prevalence and incidence of leading causes of death and disease burden

BACKGROUND

Individuals of the same chronological age display different rates of biological ageing. A number of measures of biological age have been proposed which harness age-related changes in DNA methylation profiles. These measures include five 'epigenetic clocks' which provide an index of how much an individual's biological age differs from their chronological age at the time of measurement. The five clocks encompass methylation-based predictors of chronological age (HorvathAge, HannumAge), all-cause mortality (DNAm PhenoAge, DNAm GrimAge) and telomere length (DNAm Telomere Length). A sixth epigenetic measure of ageing differs from these clocks in that it acts as a speedometer providing a single time-point measurement of the pace of an individual's biological ageing. This measure of ageing is termed DunedinPoAm. In this study, we test the association between these six epigenetic measures of ageing and the prevalence and incidence of the leading causes of disease burden and mortality in high-income countries (n ≤ 9537, Generation Scotland: Scottish Family Health Study).

RESULTS

DNAm GrimAge predicted incidence of clinically diagnosed chronic obstructive pulmonary disease (COPD), type 2 diabetes and ischemic heart disease after 13 years of follow-up (hazard ratios = 2.22, 1.52 and 1.41, respectively). DunedinPoAm predicted the incidence of COPD and lung cancer (hazard ratios = 2.02 and 1.45, respectively). DNAm PhenoAge predicted incidence of type 2 diabetes (hazard ratio = 1.54). DNAm Telomere Length associated with the incidence of ischemic heart disease (hazard ratio = 0.80). DNAm GrimAge associated with all-cause mortality, the prevalence of COPD and spirometry measures at the study baseline. These associations were present after adjusting for possible confounding risk factors including alcohol consumption, body mass index, deprivation, education and tobacco smoking and surpassed stringent Bonferroni-corrected significance thresholds.

CONCLUSIONS

Our data suggest that epigenetic measures of ageing may have utility in clinical settings to complement gold-standard methods for disease assessment and management.

Characterisation of an inflammation-related epigenetic score and its association with cognitive ability

BACKGROUND

Chronic systemic inflammation has been associated with incident dementia, but its association with age-related cognitive decline is less clear. The acute responses of many inflammatory biomarkers mean they may provide an unreliable picture of the chronicity of inflammation. Recently, a large-scale epigenome-wide association study identified DNA methylation correlates of C-reactive protein (CRP)-a widely used acute-phase inflammatory biomarker. DNA methylation is thought to be relatively stable in the short term, marking it as a potentially useful signature of exposure.

METHODS

We utilise a DNA methylation-based score for CRP and investigate its trajectories with age, and associations with cognitive ability in comparison with serum CRP and a genetic CRP score in a longitudinal study of older adults (n = 889) and a large, cross-sectional cohort (n = 7028).

RESULTS

We identified no homogeneous trajectories of serum CRP with age across the cohorts, whereas the epigenetic CRP score was consistently found to increase with age (standardised β = 0.07 and 0.01) and to do so more rapidly in males compared to females. Additionally, the epigenetic CRP score had higher test-retest reliability compared to serum CRP, indicating its enhanced temporal stability. Higher serum CRP was not found to be associated with poorer cognitive ability (standardised β = - 0.08 and - 0.05); however, a consistent negative association was identified between cognitive ability and the epigenetic CRP score in both cohorts (standardised β = - 0.15 and - 0.08).

CONCLUSIONS

An epigenetic proxy of CRP may provide a more reliable signature of chronic inflammation, allowing for more accurate stratification of individuals, and thus clearer inference of associations with incident health outcomes.

Multi-method genome- and epigenome-wide studies of inflammatory protein levels in healthy older adults

BACKGROUND

The molecular factors which control circulating levels of inflammatory proteins are not well understood. Furthermore, association studies between molecular probes and human traits are often performed by linear model-based methods which may fail to account for complex structure and interrelationships within molecular datasets.

METHODS

In this study, we perform genome- and epigenome-wide association studies (GWAS/EWAS) on the levels of 70 plasma-derived inflammatory protein biomarkers in healthy older adults (Lothian Birth Cohort 1936; n = 876; Olink® inflammation panel). We employ a Bayesian framework (BayesR+) which can account for issues pertaining to data structure and unknown confounding variables (with sensitivity analyses using ordinary least squares- (OLS) and mixed model-based approaches).

RESULTS

We identified 13 SNPs associated with 13 proteins (n = 1 SNP each) concordant across OLS and Bayesian methods. We identified 3 CpG sites spread across 3 proteins (n = 1 CpG each) that were concordant across OLS, mixed-model and Bayesian analyses. Tagged genetic variants accounted for up to 45% of variance in protein levels (for MCP2, 36% of variance alone attributable to 1 polymorphism). Methylation data accounted for up to 46% of variation in protein levels (for CXCL10). Up to 66% of variation in protein levels (for VEGFA) was explained using genetic and epigenetic data combined. We demonstrated putative causal relationships between CD6 and IL18R1 with inflammatory bowel disease and between IL12B and Crohn's disease.

CONCLUSIONS

Our data may aid understanding of the molecular regulation of the circulating inflammatory proteome as well as causal relationships between inflammatory mediators and disease.