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.

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.

Reference Ranges for NT-proBNP (N-Terminal Pro-B-Type Natriuretic Peptide) and Risk Factors for Higher NT-proBNP Concentrations in a Large General Population Cohort.

BACKGROUND: Demographic differences in expected NT-proBNP (N-terminal pro-B-type natriuretic peptide) concentration are not well established. We aimed to establish reference ranges for NT-proBNP and explore the determinants of moderately elevated NT-proBNP under the universal definition of heart failure criteria. METHODS: This is a cross-sectional study. NT-proBNP was measured in serum from 18 356 individuals without previous cardiovascular disease in the Generation Scotland Scottish Family Health Study. Age- and sex-stratified medians and 97.5th centiles were generated. Sex stratified risk factors for moderately elevated NT-proBNP (≥125 pg/mL) were investigated. RESULTS: In males, median (97.5th centile) NT-proBNP concentration at age <30 years was 21 (104) pg/mL, rising to 38 (195) pg/ml at 50 to 59 years, and 281 (6792) pg/mL at ≥80 years. In females, median NT-proBNP at age <30 years was 51 (196) pg/mL, 66 (299) pg/mL at 50 to 59 years, and 240 (2704) pg/mL at ≥80 years. At age <30 years, 9.8% of females and 1.4% of males had elevated NT-proBNP, rising to 76.5% and 81.0%, respectively, at age ≥80 years. After adjusting for risk factors, an NT-proBNP ≥125 pg/mL was more common in females than males (OR, 9.48 [95% CI, 5.60-16.1]). Older age and smoking were more strongly associated with elevated NT-proBNP in males than in females (Psex interaction <0.001, 0.07, respectively). Diabetes was inversely associated with odds of elevated NT-proBNP in females only (Psex interaction=0.007). CONCLUSIONS: An NT-proBNP ≥125 pg/mL is common in females without classical cardiovascular risk factors as well as older people. If NT-proBNP becomes widely used for screening in the general population, interpretation of NT-proBNP levels will require that age and sex-specific thresholds are used to identify patients with potential pathophysiology.

Reference ranges for GDF-15, and risk factors associated with GDF-15, in a large general population cohort.

OBJECTIVES: Growth differentiation factor (GDF)-15 is attracting interest as a biomarker in several areas of medicine. We aimed to evaluate the reference range for GDF-15 in a general population, and to explore demographics, classical cardiovascular disease risk factors, and other cardiac biomarkers associated with GDF-15. METHODS: GDF-15 was measured in serum from 19,462 individuals in the Generation Scotland Scottish Family Health Study. Associations of cardiometabolic risk factors with GDF-15 were tested using adjusted linear regression. Among 18,507 participants with no heart disease, heart failure, or stroke, and not pregnant, reference ranges (median and 97.5th centiles) were derived by decade age bands and sex. RESULTS: Among males in the reference range population, median (97.5th centile) GDF-15 concentration at age <30 years was 537 (1,135) pg/mL, rising to 931 (2,492) pg/mL at 50-59 years, and 2,152 (5,972) pg/mL at ≥80 years. In females, median GDF-15 at age <30 years was 628 (2,195) pg/mL, 881 (2,323) pg/mL at 50-59 years, and 1847 (6,830) pg/mL at ≥80 years. Among those known to be pregnant, median GDF-15 was 19,311 pg/mL. After adjustment, GDF-15 was higher in participants with adverse cardiovascular risk factors, including current smoking (+26.1%), those with previous heart disease (+12.7%), stroke (+17.1%), heart failure (+25.3%), and particularly diabetes (+60.2%). GDF-15 had positive associations with cardiac biomarkers cardiac troponin I, cardiac troponin T, and N-terminal pro B-type natriuretic peptide (NT-proBNP). CONCLUSIONS: These data define reference ranges for GDF-15 for comparison in future studies, and identify potentially confounding risk factors and mediators to be considered in interpreting GDF-15 concentrations.

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.

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 (ALSPACadults, 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 (MSp-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 (pnominal<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 (NALSPAC=565; NGS=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.

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.

Although our genetic code does not change throughout our lives, our genes can be turned on and off as a result of epigenetics. Epigenetics can track how the environment and even certain behaviors add or remove small chemical markers to the DNA that makes up the genome. The type and location of these markers may affect whether genes are active or silent, this is, whether the protein coded for by that gene is being produced or not. One common epigenetic marker is known as DNA methylation. DNA methylation has been linked to the levels of a range of proteins in our cells and the risk people have of developing chronic diseases. Blood samples can be used to determine the epigenetic markers a person has on their genome and to study the abundance of many proteins. Gadd, Hillary, McCartney, Zaghlool et al. studied the relationships between DNA methylation and the abundance of 953 different proteins in blood samples from individuals in the German KORA cohort and the Scottish Lothian Birth Cohort 1936. They then used machine learning to analyze the relationship between epigenetic markers found in people's blood and the abundance of proteins, obtaining epigenetic scores or 'EpiScores' for each protein. They found 109 proteins for which DNA methylation patterns explained between at least 1% and up to 58% of the variation in protein levels. Integrating the 'EpiScores' with 14 years of medical records for more than 9000 individuals from the Generation Scotland study revealed 130 connections between EpiScores for proteins and a future diagnosis of common adverse health outcomes. These included diabetes, stroke, depression, various cancers, and inflammatory conditions such as rheumatoid arthritis and inflammatory bowel disease. Age-related chronic diseases are a growing issue worldwide and place pressure on healthcare systems. They also severely reduce quality of life for individuals over many years. This work shows how epigenetic scores based on protein levels in the blood could predict a person's risk of several of these diseases. In the case of type 2 diabetes, the EpiScore results replicated previous research linking protein levels in the blood to future diagnosis of diabetes. Protein EpiScores could therefore allow researchers to identify people with the highest risk of disease, making it possible to intervene early and prevent these people from developing chronic conditions as they age.

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.

Genome-wide interaction study with major depression identifies novel variants associated with cognitive function.

Major Depressive Disorder (MDD) often is associated with significant cognitive dysfunction. We conducted a meta-analysis of genome-wide interaction of MDD and cognitive function using data from four large European cohorts in a total of 3510 MDD cases and 6057 controls. In addition, we conducted analyses using polygenic risk scores (PRS) based on data from the Psychiatric Genomics Consortium (PGC) on the traits of MDD, Bipolar disorder (BD), Schizophrenia (SCZ), and mood instability (MIN). Functional exploration contained gene expression analyses and Ingenuity Pathway Analysis (IPA®). We identified a set of significantly interacting single nucleotide polymorphisms (SNPs) between MDD and the genome-wide association study (GWAS) of cognitive domains of executive function, processing speed, and global cognition. Several of these SNPs are located in genes expressed in brain, with important roles such as neuronal development (REST), oligodendrocyte maturation (TNFRSF21), and myelination (ARFGEF1). IPA® identified a set of core genes from our dataset that mapped to a wide range of canonical pathways and biological functions (MPO, FOXO1, PDE3A, TSLP, NLRP9, ADAMTS5, ROBO1, REST). Furthermore, IPA® identified upstream regulator molecules and causal networks impacting on the expression of dataset genes, providing a genetic basis for further clinical exploration (vitamin D receptor, beta-estradiol, tadalafil). PRS of MIN and meta-PRS of MDD, MIN and SCZ were significantly associated with all cognitive domains. Our results suggest several genes involved in physiological processes for the development and maintenance of cognition in MDD, as well as potential novel therapeutic agents that could be explored in patients with MDD associated cognitive dysfunction.

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.

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.

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.

Genetic and shared couple environmental contributions to smoking and alcohol use in the UK population.

Alcohol use and smoking are leading causes of death and disability worldwide. Both genetic and environmental factors have been shown to influence individual differences in the use of these substances. In the present study we tested whether genetic factors, modelled alongside common family environment, explained phenotypic variance in alcohol use and smoking behaviour in the Generation Scotland (GS) family sample of up to 19,377 individuals. SNP and pedigree-associated effects combined explained between 18 and 41% of the variance in substance use. Shared couple effects explained a significant amount of variance across all substance use traits, particularly alcohol intake, for which 38% of the phenotypic variance was explained. We tested whether the within-couple substance use associations were due to assortative mating by testing the association between partner polygenic risk scores in 34,987 couple pairs from the UK Biobank (UKB). No significant association between partner polygenic risk scores were observed. Associations between an individual's alcohol PRS (b = 0.05, S.E. = 0.006, p < 2 × 10-16) and smoking status PRS (b = 0.05, S.E. = 0.005, p < 2 × 10-16) were found with their partner's phenotype. In support of this, G carriers of a functional ADH1B polymorphism (rs1229984), known to be associated with greater alcohol intake, were found to consume less alcohol if they had a partner who carried an A allele at this SNP. Together these results show that the shared couple environment contributes significantly to patterns of substance use. It is unclear whether this is due to shared environmental factors, assortative mating, or indirect genetic effects. Future studies would benefit from longitudinal data and larger sample sizes to assess this further.

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.

Bayesian reassessment of the epigenetic architecture of complex traits.

Linking epigenetic marks to clinical outcomes improves insight into molecular processes, disease prediction, and therapeutic target identification. Here, a statistical approach is presented to infer the epigenetic architecture of complex disease, determine the variation captured by epigenetic effects, and estimate phenotype-epigenetic probe associations jointly. Implicitly adjusting for probe correlations, data structure (cell-count or relatedness), and single-nucleotide polymorphism (SNP) marker effects, improves association estimates and in 9,448 individuals, 75.7% (95% CI 71.70-79.3) of body mass index (BMI) variation and 45.6% (95% CI 37.3-51.9) of cigarette consumption variation was captured by whole blood methylation array data. Pathway-linked probes of blood cholesterol, lipid transport and sterol metabolism for BMI, and xenobiotic stimuli response for smoking, showed >1.5 times larger associations with >95% posterior inclusion probability. Prediction accuracy improved by 28.7% for BMI and 10.2% for smoking over a LASSO model, with age-, and tissue-specificity, implying associations are a phenotypic consequence rather than causal.

Variants associated with HHIP expression have sex-differential effects on lung function.

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P<5x10 -8) interactions with sex on lung function, and 21 showed suggestive interactions (P<1x10 -6). The strongest signal, from rs7697189 (chr4:145436894) on forced expiratory volume in 1 second (FEV 1) (P=3.15x10 -15), was replicated (P=0.016) in SpiroMeta. The C allele increased FEV 1 more in males (untransformed FEV 1 ß=0.028 [SE 0.0022] litres) than females (ß=0.009 [SE 0.0014] litres), and this effect was not accounted for by differential effects on height, smoking or pubertal age. rs7697189 resides upstream of the hedgehog-interacting protein ( HHIP) gene and was previously associated with lung function and HHIP lung expression. We found HHIP expression was significantly different between the sexes (P=6.90x10 -6), but we could not detect sex differential effects of rs7697189 on expression. Conclusions: We identified a novel genotype-by-sex interaction at a putative enhancer region upstream of the HHIP gene. Establishing the mechanism by which HHIP SNPs have different effects on lung function in males and females will be important for our understanding of lung health and diseases in both sexes.

A meta-analysis of genome-wide association studies of epigenetic age acceleration.

'Epigenetic age acceleration' is a valuable biomarker of ageing, predictive of morbidity and mortality, but for which the underlying biological mechanisms are not well established. Two commonly used measures, derived from DNA methylation, are Horvath-based (Horvath-EAA) and Hannum-based (Hannum-EAA) epigenetic age acceleration. We conducted genome-wide association studies of Horvath-EAA and Hannum-EAA in 13,493 unrelated individuals of European ancestry, to elucidate genetic determinants of differential epigenetic ageing. We identified ten independent SNPs associated with Horvath-EAA, five of which are novel. We also report 21 Horvath-EAA-associated genes including several involved in metabolism (NHLRC, TPMT) and immune system pathways (TRIM59, EDARADD). GWAS of Hannum-EAA identified one associated variant (rs1005277), and implicated 12 genes including several involved in innate immune system pathways (UBE2D3, MANBA, TRIM46), with metabolic functions (UBE2D3, MANBA), or linked to lifespan regulation (CISD2). Both measures had nominal inverse genetic correlations with father's age at death, a rough proxy for lifespan. Nominally significant genetic correlations between Hannum-EAA and lifestyle factors including smoking behaviours and education support the hypothesis that Hannum-based epigenetic ageing is sensitive to variations in environment, whereas Horvath-EAA is a more stable cellular ageing process. We identified novel SNPs and genes associated with epigenetic age acceleration, and highlighted differences in the genetic architecture of Horvath-based and Hannum-based epigenetic ageing measures. Understanding the biological mechanisms underlying individual differences in the rate of epigenetic ageing could help explain different trajectories of age-related decline.

Identification of novel differentially methylated sites with potential as clinical predictors of impaired respiratory function and COPD.

BACKGROUND: The causes of poor respiratory function and COPD are incompletely understood, but it is clear that genes and the environment play a role. As DNA methylation is under both genetic and environmental control, we hypothesised that investigation of differential methylation associated with these phenotypes would permit mechanistic insights, and improve prediction of COPD. We investigated genome-wide differential DNA methylation patterns using the recently released 850 K Illumina EPIC array. This is the largest single population, whole-genome epigenetic study to date. METHODS: Epigenome-wide association studies (EWASs) of respiratory function and COPD were performed in peripheral blood samples from the Generation Scotland: Scottish Family Health Study (GS:SFHS) cohort (n = 3781; 274 COPD cases and 2919 controls). In independent COPD incidence data (n = 149), significantly differentially methylated sites (DMSs; p < 3.6 × 10-8) were evaluated for their added predictive power when added to a model including clinical variables, age, sex, height and smoking history using receiver operating characteristic analysis. The Lothian Birth Cohort 1936 (LBC1936) was used to replicate association (n = 895) and prediction (n = 178) results. FINDINGS: We identified 28 respiratory function and/or COPD associated DMSs, which mapped to genes involved in alternative splicing, JAK-STAT signalling, and axon guidance. In prediction analyses, we observed significant improvement in discrimination between COPD cases and controls (p < .05) in independent GS:SFHS (p = .016) and LBC1936 (p = .010) datasets by adding DMSs to a clinical model. INTERPRETATION: Identification of novel DMSs has provided insight into the molecular mechanisms regulating respiratory function and aided prediction of COPD risk. Further studies are needed to assess the causality and clinical utility of identified associations. FUND: Wellcome Trust Strategic Award 10436/Z/14/Z.

Assessment of dried blood spots for DNA methylation profiling.

Background: DNA methylation reflects health-related environmental exposures and genetic risk, providing insights into aetiological mechanisms and potentially predicting disease onset, progression and treatment response. An increasingly recognised need for large-scale, longitudinally-profiled samples collected world-wide has made the development of efficient and straightforward sample collection and storage procedures a pressing issue. An alternative to the low-temperature storage of EDTA tubes of venous blood samples, which are frequently the source of the DNA used in such studies, is to collect and store at room temperature blood samples using purpose built filter paper, such as Whatman FTA® cards. Our goal was to determine whether DNA stored in this manner can be used to generate DNA methylation profiles comparable to those generated using blood samples frozen in EDTA tubes. Methods: DNA methylation profiles were obtained from matched EDTA tube and Whatman FTA® card whole-blood samples from 62 Generation Scotland: Scottish Family Health Study participants using the Infinium HumanMethylation450 BeadChip. Multiple quality control procedures were implemented, the relationship between the two sample types assessed, and epigenome-wide association studies (EWASs) performed for smoking status, age and the interaction between these variables and sample storage method. Results: Dried blood spot (DBS) DNA methylation profiles were of good quality and DNA methylation profiles from matched DBS and EDTA tube samples were highly correlated (mean r = 0.991) and could distinguish between participants. EWASs replicated established associations for smoking and age, with no evidence for moderation by storage method. Conclusions: Our results support the use of Whatman FTA® cards for collecting and storing blood samples for DNA methylation profiling. This approach is likely to be particularly beneficial for large-scale studies and those carried out in areas where freezer access is limited. Furthermore, our results will inform consideration of the use of newborn heel prick DBSs for research use.