Search

Testing the effect of rare variants on phenotypic variation is difficult due to the need for extremely large cohorts to identify associated variants given expected effect sizes. An alternative approach is to investigate the effect of rare genetic variants on low-level genomic traits, such as gene expression or DNA methylation (DNAm), as effect sizes are expected to be larger for low-level compared to higher-order complex traits. Here, we investigate DNAm in healthy ageing populations - the Lothian Birth cohorts of 1921 and 1936 and identify both transient and stable outlying DNAm levels across the genome. We find an enrichment of rare genetic variants within 1kb of DNAm sites in individuals with stable outlying DNAm, implying genetic control of this extreme variation. Using a family-based cohort, the Brisbane Systems Genetics Study, we observed increased sharing of DNAm outliers among more closely related individuals, consistent with these outliers being driven by rare genetic variation. We demonstrated that outlying DNAm levels have a functional consequence on gene expression levels, with extreme levels of DNAm being associated with gene expression levels towards the tails of the population distribution. Overall, this study demonstrates the role of rare variants in the phenotypic variation of low-level genomic traits, and the effect of extreme levels of DNAm on gene expression.

BackgroundAmyotrophic lateral sclerosis (ALS), the most predominant form of Motor Neuron Disease (MND), is a progressive and fatal neurodegenerative condition that spreads throughout the neuromotor system by afflicting upper and lower motor neurons. Lower motor neurons project from the central nervous system and innervate muscle fibres at motor endplates, which degrade over the course of the disease leading to muscle weakness. The direction of neurodegeration from or to the point of neuromuscular junctions and the role of muscle itself in pathogenesis has continued to be a topic of debate in ALS research.MethodsTo assess the variation in gene expression between affected and nonaffected muscle tissue that might lead to this local degeneration of motor units, we generated RNA-seq skeletal muscle transcriptomes from 28 MND cases and 18 healthy controls and conducted differential expression analyses on gene-level counts, as well as an isoform switching analysis on isoform-level counts.ResultsWe identified 52 differentially-expressed genes (Benjamini-Hochberg-adjustedp< 0.05) within this comparison, including 38 protein coding, 9 long non-coding RNA, and 5 pseudogenes. Of protein-coding genes, 31 were upregulated in cases including with notable genes including the collagenicCOL25A1(p= 3.1 × 10−10),SAA1which is released in response to tissue injury (p= 3.6 × 10−5) as well as others of the SAA family, and the actin-encodingACTC1(p= 2.3 × 10−5). Additionally, we identified 17 genes which exhibited a functional isoform switch with likely functional consequences between cases and controls.ConclusionsOur analyses provide evidence of increased tissue generation in MND cases, which likely serve to compensate for the degeneration of motor units and skeletal muscle.

Background: Thyroid hormones play a key role in differentiation and metabolism and are known regulators of gene expression through both genomic and epigenetic processes including DNA methylation. The aim of this study was to examine associations between thyroid hormones and DNA methylation. Methods: We carried out a fixed-effect meta-analysis of epigenome-wide association study (EWAS) of blood DNA methylation sites from 8 cohorts from the ThyroidOmics Consortium, incorporating up to 7073 participants of both European and African ancestry, implementing a discovery and replication stage. Statistical analyses were conducted using normalized beta CpG values as dependent and log-transformed thyrotropin (TSH), free thyroxine, and free triiodothyronine levels, respectively, as independent variable in a linear model. The replicated findings were correlated with gene expression levels in whole blood and tested for causal influence of TSH and free thyroxine by two-sample Mendelian randomization (MR). Results: Epigenome-wide significant associations (p-value

The recent increase in obesity levels across many countries is likely to be driven by nongenetic factors. The epigenetic modification DNA methylation (DNAm) may help to explore this as it is sensitive to both genetic and environmental exposures. While the relationship between DNAm and body fat traits has been extensively studied [1–9], there is limited literature on the shared associations of DNAm variation across such traits. Akin to genetic correlation estimates, which measure the degree of common genetic control between two traits, here we introduce an approach to evaluate the similarities in DNAm associations between traits, DNAm correlations. As DNAm can be both a cause and consequence of complex traits [5, 10, 11], DNAm correlations have the potential to provide novel insights into trait relationships above that currently obtained from genetic and phenotypic correlations. Utilising 7,519 unrelated individuals from Generation Scotland (GS), we calculated DNAm correlations using the bivariate OREML framework in the OSCA software [12] to investigate the shared associations of DNAm variation between traits. For each trait we also estimated the shared contribution of DNAm between sexes. We identified strong, positive DNAm correlations between each of the body fat traits (BMI, body fat % and waist to hip ratio; ranging from 0.96 to 1.00), finding larger associations than those identified by genetic and phenotypic correlations. We identified a significant deviation from 1 in the rDNAmfor BMI between males and females, with sex-specific DNAm changes associated with BMI identified at eight DNAm probes. Employing genome-wide DNAm correlations to evaluate the similarities in the associations of DNAm with complex traits has provided novel insight into obesity related traits beyond that provided by genetic correlations.

Endometriosis is a leading cause of pain and infertility affecting millions of women globally. Identifying biologic and genetic effects on DNA methylation (DNAm) in endometrium increases understanding of mechanisms that influence gene regulation predisposing to endometriosis and offers an opportunity for novel therapeutic target discovery. Herein, we characterize variation in endometrial DNAm and its association with menstrual cycle phase, endometriosis, and genetic variants through analysis of genome-wide genotype data and methylation at 759,345 DNAm sites in endometrial samples from 984 deeply-phenotyped participants. We identify significant differences in DNAm profiles between menstrual cycle phases and at four DNAm sites between stage III/IV endometriosis and controls. We estimate that 15.4% of the variation in endometriosis is captured by DNAm, and identify DNAm networks associated with endometriosis. DNAm quantitative trait locus (mQTL) analysis identified 118,185 independentcis-mQTL including some tissue-specific effects. We find significant differences in DNAm profiles between endometriosis sub- phenotypes and a significant association between genetic regulation of methylation in endometrium and disease risk, providing functional evidence for genomic targets contributing to endometriosis risk and pathogenesis.

The complement system, including complement components 3 and 4 (C3, C4), traditionally has been linked to innate immunity. More recently, complement components have also been implicated in brain development and the risk of schizophrenia. Based on a large, population-based case-cohort study, we measured the blood concentrations of C3 and C4 in 68,768 neonates. We found a strong correlation between the concentrations of C3 and C4 (phenotypic correlation = 0.65,P-value < 1.0×10−100, genetic correlation = 0.38,P-value = 1.9×10−35). A genome-wide association study (GWAS) for C4 protein concentration identified 36 independent loci, 30 of which were in or near the major histocompatibility complex on chromosome 6 (which includes theC4gene), while six loci were found on six other chromosomes. A GWAS for C3 identified 15 independent loci, seven of which were located in theC3gene on chromosome 19, and eight loci on five other chromosomes. We found no association between (a) measured neonatal C3 and C4 concentrations, imputed C4 haplotypes, or predictedC4gene expression, with (b) schizophrenia (SCZ), bipolar disorder (BIP), depression (DEP), autism spectrum disorder, attention deficit hyperactivity disorder or anorexia nervosa diagnosed in later life. Mendelian randomisation (MR) suggested a small positive association between higher C4 protein concentration and an increased risk of SCZ, BIP, and DEP, but these findings did not persist in more stringent analyses. Evidence from MR supported causal relationships between C4 concentration and several autoimmune disorders: systemic lupus erythematosus (SLE, OR and 95% confidence interval, 0.37, 0.34 – 0.42); type-1 diabetes (T1D, 0.54, 0.50 - 0.58); multiple sclerosis (MS, 0.68, 0.63 - 0.74); rheumatoid arthritis (0.85, 0.80 - 0.91); and Crohn’s disease (1.26, 1.19 - 1.34). A phenome-wide association study (PheWAS) in UK Biobank confirmed that the genetic correlates of C4 concentration were associated a range of autoimmune disorders including coeliac disease, thyrotoxicosis, hypothyroidism, T1D, sarcoidosis, psoriasis, SLE and ankylosing spondylitis. We found no evidence of associations between C3 versus mental or autoimmune disorders based on either MR or PheWAS. In general, our results do not support the hypothesis that C4 is causally associated with the risk of SCZ (nor several other mental disorders). We provide new evidence to support the hypothesis that higher C4 concentration is associated with lower risks of autoimmune disorders.

Strong associations between neural tube defects (NTDs) and monozygotic (MZ) twinning have long been noted, and it has been suggested that NTD cases who do not present as MZ twins may be the survivors of MZ twinning events. We have recently shown that MZ twins carry a strong, distinctive DNA methylation signature and have developed an algorithm based on genomewide DNA methylation array data that distinguishes MZ twins from dizygotic twins and other relatives at well above chance level. We have applied this algorithm to published methylation data from five fetal tissues (placental chorionic villi, kidney, spinal cord, brain and muscle) collected from spina bifida cases (n = 22), anencephalic cases (n = 15) and controls (n = 19). We see no difference in signature between cases and controls, providing no support for a common etiological role of MZ twinning in NTDs. The strong associations therefore continue to await elucidation.

Amyotrophic Lateral Sclerosis (ALS) is recognised to be a complex neurodegenerative disease involving both genetic and non-genetic risk factors. The underlying causes and risk factors for the majority of cases remain unknown; however, ever-larger genetic data studies and methodologies promise an enhanced understanding. Recent analyses using published summary statistics from the largest ALS genome-wide association study (GWAS) (20,806 ALS cases and 59,804 healthy controls) identified that schizophrenia (SCZ), cognitive performance (CP) and educational attainment (EA) related traits were genetically correlated with ALS. To provide additional evidence for these correlations, we built single and multi-trait genetic predictors using GWAS summary statistics for ALS and these traits, (SCZ, CP, EA) in an independent Australian cohort (846 ALS cases and 665 healthy controls). We compared methods for generating the risk predictors and found that the combination of traits improved the prediction (Nagelkerke-R-2) of the case-control logistic regression. The combination of ALS, SCZ, CP, and EA, using the SBayesR predictor method gave the highest prediction (Nagelkerke-R-2) of 0.027 (P value = 4.6 x 10(-8)), with the odds-ratio for estimated disease risk between the highest and lowest deciles of individuals being 3.15 (95% CI 1.96-5.05). These results support the genetic correlation between ALS, SCZ, CP and EA providing a better understanding of the complexity of ALS. Refereed/Peer-reviewed

Genetic association studies have identified 44 common genome-wide significant risk loci for late-onset Alzheimer’s disease (LOAD). However, LOAD genetic architecture and prediction are unclear. Here we estimate the optimal P-threshold (Poptimal) of a genetic risk score (GRS) for prediction of LOAD in three independent datasets comprising 676 cases and 35,675 family history proxy cases. We show that the discriminative ability of GRS in LOAD prediction is maximised when selecting a small number of SNPs. Both simulation results and direct estimation indicate that the number of causal common SNPs for LOAD may be less than 100, suggesting LOAD is more oligogenic than polygenic. The best GRS explains approximately 75% of SNP-heritability, and individuals in the top decile of GRS have ten-fold increased odds when compared to those in the bottom decile. In addition, 14 variants are identified that contribute to both LOAD risk and age at onset of LOAD., Despite the identification of genetic risk loci for late-onset Alzheimer’s disease (LOAD), the genetic architecture and prediction remains unclear. Here, the authors use genetic risk scores for prediction of LOAD across three datasets and show evidence suggesting oligogenic variant architecture for this disease.

DNA methylation has fundamental roles in gene programming and aging that may help predict mortality. However, no large-scale study has investigated whether site-specific DNA methylation predicts all-cause mortality. We used the Illumina-HumanMethylation450-BeadChip to identify blood DNA methylation sites associated with all-cause mortality for 12, 300 participants in 12 Cohorts of the Heart and Aging Research in Genetic Epidemiology (CHARGE) Consortium. Over an average 10-year follow-up, there were 2,561 deaths across the cohorts. Nine sites mapping to three intergenic and six gene-specific regions were associated with mortality (P < 9.3x10-7) independently of age and other mortality predictors. Six sites (cg14866069, cg23666362, cg20045320, cg07839457, cg07677157, cg09615688)—mapping respectively to BMPR1B, MIR1973, IFITM3, NLRC5, and two intergenic regions—were associated with reduced mortality risk. The remaining three sites (cg17086398, cg12619262, cg18424841)—mapping respectively to SERINC2, CHST12, and an intergenic region—were associated with increased mortality risk. DNA methylation at each site predicted 5%¬¬–15% of all deaths. We also assessed the causal association of those sites to age-related chronic diseases by using Mendelian randomization, identifying weak causal relationship between cg18424841 and cg09615688 with coronary heart disease. Of the nine sites, three (cg20045320, cg07839457, cg07677157) were associated with lower incidence of heart disease risk and two (cg20045320, cg07839457) with smoking and inflammation in prior CHARGE analyses. Methylation of cg20045320, cg07839457, and cg17086398 was associated with decreased expression of nearby genes (IFITM3, IRF, NLRC5, MT1, MT2, MARCKSL1) linked to immune responses and cardiometabolic diseases. These sites may serve as useful clinical tools for mortality risk assessment and preventative care.

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.

ENIGMA-CNV Working Group: van der Meer, Dennis; Sonderby, Ida E; Kaufmann, Tobias; Walters, G Bragi; Abdellaoui, Abdel; Ames, David; Amunts, Katrin; Andersson, Micael; Armstrong, Nicola J; Bernard, Manon; Blackburn, Nicholas B; Blangero, John; Boomsma, Dorret I; Brodaty, Henry; Brouwer, Rachel M; Bulow, Robin; Cahn, Wiepke; Calhoun, Vince D; Caspers, Svenja; Cavalleri, Gianpiero L; Ching, Christopher R K; Cichon, Sven; Ciufolini, Simone; Corvin, Aiden; Crespo-Facorro, Benedicto; Curran, Joanne E; Dalvie, Shareefa; Dazzan, Paola; de Geus, Eco J C; de Zubicaray, Greig I; de Zwarte, Sonja M C; Delanty, Norman; den Braber, Anouk; Desrivieres, Sylvane; Di Forti, Marta; Doherty, Joanne L; Donohoe, Gary; Ehrlich, Stefan; Eising, Else; Espeseth, Thomas; Fisher, Simon E; Fladby, Tormod; Frei, Oleksandr; Frouin, Vincent; Fukunaga, Masaki; Gareau, Thomas; Glahn, David C; Grabe, Hans J; Groenewold, Nynke A; Gustafsson, Omar; Haavik, Jan; Haberg, Asta K; Hashimoto, Ryota; Hehir-Kwa, Jayne Y; Hibar, Derrek P; Hillegers, Manon H J; Hoffmann, Per; Holleran, Laurena; Hottenga, Jouke-Jan; Hulshoff Pol, Hilleke E; Ikeda, Masashi; Jacquemont, Sebastien; Jahanshad, Neda; Jockwitz, Christiane; Johansson, Stefan; Jonsson, Erik G; Kikuchi, Masataka; Knowles, Emma E M; Kwok, John B; Le Hellard, Stephanie; Linden, David E J; Liu, Jingyu; Lundervold, Arvid; Lundervold, Astri J; Martin, Nicholas G; Mather, Karen A; Mathias, Samuel R; McMahon, Katie L; McRae, Allan F; Medland, Sarah E; Moberget, Torgeir; Moreau, Clara; Morris, Derek W; Muhleisen, Thomas W; Murray, Robin M; Nordvik, Jan E; Nyberg, Lars; Olde Loohuis, Loes M; Ophoff, Roel A; Owen, Michael J; Paus, Tomas; Pausova, Zdenka; Peralta, Juan M; Pike, Bruce; Prieto, Carlos; Quinlan, Erin Burke; Reinbold, Celine S; Reis Marques, Tiago; Rucker, James J H; Sachdev, Perminder S; Sando, Sigrid B; Schofield, Peter R; Schork, Andrew J; Schumann, Gunter; Shin, Jean; Shumskaya, Elena; Silva, Ana I; Sisodiya, Sanjay M; Steen, Vidar M; Stein, Dan J; Strike, Lachlan T; Tamnes, Christian K; Teumer, Alexander; Thalamuthu, Anbupalam; Tordesillas-Gutierrez, Diana; Uhlmann, Anne; Ulfarsson, Magnus O; van 't Ent, Dennis; van den Bree, Marianne B M; Vassos, Evangelos; Wen, Wei; Wittfeld, Katharina; Wright, Margaret J; Zayats, Tetyana; Dale, Anders M; Djurovic, Srdjan; Agartz, Ingrid; Westlye, Lars T; Stefansson, Hreinn; Stefansson, Kari; Thompson, Paul M; Andreassen, Ole A., [Importance] Recurrent microdeletions and duplications in the genomic region 15q11.2 between breakpoints 1 (BP1) and 2 (BP2) are associated with neurodevelopmental disorders. These structural variants are present in 0.5% to 1.0% of the population, making 15q11.2 BP1-BP2 the site of the most prevalent known pathogenic copy number variation (CNV). It is unknown to what extent this CNV influences brain structure and affects cognitive abilities., [Objective] To determine the association of the 15q11.2 BP1-BP2 deletion and duplication CNVs with cortical and subcortical brain morphology and cognitive task performance., [Design, Setting, and Participants] In this genetic association study, T1-weighted brain magnetic resonance imaging were combined with genetic data from the ENIGMA-CNV consortium and the UK Biobank, with a replication cohort from Iceland. In total, 203 deletion carriers, 45 247 noncarriers, and 306 duplication carriers were included. Data were collected from August 2015 to April 2019, and data were analyzed from September 2018 to September 2019., [Main Outcomes and Measures] The associations of the CNV with global and regional measures of surface area and cortical thickness as well as subcortical volumes were investigated, correcting for age, age2, sex, scanner, and intracranial volume. Additionally, measures of cognitive ability were analyzed in the full UK Biobank cohort., [Results] Of 45 756 included individuals, the mean (SD) age was 55.8 (18.3) years, and 23 754 (51.9%) were female. Compared with noncarriers, deletion carriers had a lower surface area (Cohen d = −0.41; SE, 0.08; P = 4.9 × 10−8), thicker cortex (Cohen d = 0.36; SE, 0.07; P = 1.3 × 10−7), and a smaller nucleus accumbens (Cohen d = −0.27; SE, 0.07; P = 7.3 × 10−5). There was also a significant negative dose response on cortical thickness (β = −0.24; SE, 0.05; P = 6.8 × 10−7). Regional cortical analyses showed a localization of the effects to the frontal, cingulate, and parietal lobes. Further, cognitive ability was lower for deletion carriers compared with noncarriers on 5 of 7 tasks., [Conclusions and Relevance] These findings, from the largest CNV neuroimaging study to date, provide evidence that 15q11.2 BP1-BP2 structural variation is associated with brain morphology and cognition, with deletion carriers being particularly affected. The pattern of results fits with known molecular functions of genes in the 15q11.2 BP1-BP2 region and suggests involvement of these genes in neuronal plasticity. These neurobiological effects likely contribute to the association of this CNV with neurodevelopmental disorders., This study is supported in part by grants U54 EB20403, R01MH116147, and R56AG058854 from the National Institutes of Health, grant 609020 from the European Union Seventh Framework Programme, and grants 223273 and 276082 from the Research Council Norway. Part of this work was performed using the Service for Sensitive Data (TSD), which is developed and operated by the TSD Service Groupand owned by the University of Oslo.

An improved understanding of etiological mechanisms in Parkinson’s disease (PD) is urgently needed because the number of affected individuals is projected to increase rapidly as populations age. We present results from a blood-based methylome-wide association study of PD involving meta-analysis of 229 K CpG probes in 1,132 cases and 999 controls from two independent cohorts. We identify two previously unreported epigenome-wide significant associations with PD, including cg06690548 on chromosome 4. We demonstrate that cg06690548 hypermethylation in PD is associated with down-regulation of the SLC7A11 gene and show this is consistent with an environmental exposure, as opposed to medications or genetic factors with effects on DNA methylation or gene expression. These findings are notable because SLC7A11 codes for a cysteine-glutamate anti-porter regulating levels of the antioxidant glutathione, and it is a known target of the environmental neurotoxin β-methylamino-L-alanine (BMAA). Our study identifies the SLC7A11 gene as a plausible biological target in PD., Parkinson’s disease (PD) is a common neurodegenerative disorder with a complex etiology involving genetics and the environment. Here, Vallerga et al. identify two CpG probes associated with PD in a blood cell type-corrected epigenome-wide meta-analysis, implicating the SLC7A11 gene as a plausible biological target.

Publisher's version (útgefin grein), Carriers of large recurrent copy number variants (CNVs) have a higher risk of developing neurodevelopmental disorders. The 16p11.2 distal CNV predisposes carriers to e.g., autism spectrum disorder and schizophrenia. We compared subcortical brain volumes of 12 16p11.2 distal deletion and 12 duplication carriers to 6882 non-carriers from the large-scale brain Magnetic Resonance Imaging collaboration, ENIGMA-CNV. After stringent CNV calling procedures, and standardized FreeSurfer image analysis, we found negative dose-response associations with copy number on intracranial volume and on regional caudate, pallidum and putamen volumes (β = −0.71 to −1.37; P < 0.0005). In an independent sample, consistent results were obtained, with significant effects in the pallidum (β = −0.95, P = 0.0042). The two data sets combined showed significant negative dose-response for the accumbens, caudate, pallidum, putamen and ICV (P = 0.0032, 8.9 × 10−6, 1.7 × 10− 9, 3.5 × 10−12 and 1.0 × 10−4, respectively). Full scale IQ was lower in both deletion and duplication carriers compared to non-carriers. This is the first brain MRI study of the impact of the 16p11.2 distal CNV, and we demonstrate a specific effect on subcortical brain structures, suggesting a neuropathological pattern underlying the neurodevelopmental syndromes., 1000BRAINS: 1000BRAINS is a population-based cohort based on the Heinz-Nixdorf Recall Study and is supported in part by the German National Cohort. We thank the Heinz Nixdorf Foundation (Germany) for their generous support in terms of the Heinz Nixdorf Study. The HNR study is also supported by the German Ministry of Education and Science (FKZ 01EG940), and the German Research Council (DFG, ER 155/6-1). The authors are supported by the Initiative and Networking Fund of the Helmholtz Association (Svenja Caspers) and the European Union’s Horizon 2020 Research and Innovation Programme under Grant Agreement 7202070 (Human Brain Project SGA1; Katrin Amunts, Sven Cichon). This work was further supported by the German Federal Ministry of Education and Research (BMBF) through the Integrated Network IntegraMent (Integrated Understanding of Causes and Mechanisms in Mental Disorders) under the auspices of the e:Med Program (grant 01ZX1314A to M.M.N. and S.C.), and by the Swiss National Science Foundation (SNSF, grant 156791 to S.C.). 16p.11.2 European Consortium: B.D. is supported by the Swiss National Science Foundation (NCCR Synapsy, project grant Nr 32003B_159780) and Foundation Synapsis. LREN is very grateful to the Roger De Spoelberch and Partridge Foundations for their generous financial support. This work was supported by grants from the Simons Foundation (SFARI274424) and the Swiss National Science Foundation (31003A_160203) to A.R. and S.J. Betula: The relevant Betula data collection and analyses were supported by a grant from the Knut & Alice Wallenberg (KAW) to L. Nyberg. Brainscale: the Brainscale study was supported by the Netherlands Organization for Scientific Research MagW 480-04-004 (Dorret Boomsma), 51.02.060 (Hilleke Hulshoff Pol), 668.772 (Dorret Boomsma & Hilleke Hulshoff Pol); NWO/SPI 56-464-14192 (Dorret Boomsma), the European Research Council (ERC-230374) (Dorret Boomsma), High Potential Grant Utrecht University (Hilleke Hulshoff Pol), NWO Brain and Cognition 433-09-220 (Hilleke Hulshoff Pol). Brain Imaging Genetics (BIG): This work makes use of the BIG database, first established in Nijmegen, The Netherlands, in 2007. This resource is now part of Cognomics (www.cognomics.nl), a joint initiative by researchers of the Donders Centre for Cognitive Neuroimaging, the Human Genetics and Cognitive Neuroscience departments of the Radboud university medical centre and the Max Planck Institute for Psycholinguistics in Nijmegen. The Cognomics Initiative has received supported from the participating departments and centres and from external grants, i.e., the Biobanking and Biomolecular Resources Research Infrastructure (the Netherlands) (BBMRI-NL), the Hersenstichting Nederland, and the Netherlands Organisation for Scientific Research (NWO). The research leading to these results also receives funding from the NWO Gravitation grant ‘Language in Interaction’, the European Community’s Seventh Framework Programme (FP7/2007–2013) under grant agreements n° 602450 (IMAGEMEND), n°278948 (TACTICS), and n°602805 (Aggressotype) as well as from the European Community’s Horizon 2020 programme under grant agreement n° 643051 (MiND) and from ERC-2010-AdG 268800-NEUROSCHEMA. In addition, the work was supported by a grant for the ENIGMA Consortium (grant number U54 EB020403) from the BD2K Initiative of a cross-NIH partnership. COBRE: This work was supported by a NIH COBRE Phase I grant (1P20RR021938, Lauriello, PI and 2P20GM103472, Calhoun, PI) awarded to the Mind Research Network. We wish to express our gratitude to numerous investigators who were either external consultants to the Cores and projects, mentors on the projects, members of the external advisory committee and members of the internal advisory committee. Decode: The research leading to these results has received financial contribution from the European Union’s Seventh Framework Programme (EU-FP7/2007–2013), EU-FP7 funded grant no. 602450 (IMAGEMEND) as well as support from the Innovative Medicines Initiative Joint Undertaking under grant agreement no.115300 (EUAIMS). DemGene: Norwegian Health Association and Research Council of Norway. Dublin: Work was supported by Science Foundation Ireland (SFI grant 12/IP/1359 to Gary Donohoe and SFI08/IN.1/B1916-Corvin to Aidan C Corvin) and the European Research Council (ERC-StG-2015-677467). EPIGEN-UK (SMS, CL): The work was partly undertaken at UCLH/UCL, which received a proportion of funding from the UK Department of Health’s NIHR Biomedical Research Centres funding scheme. We are grateful to the Wolfson Trust and the Epilepsy Society for supporting the Epilepsy Society MRI scanner, and the Epilepsy Society for supporting CL. Haavik: The work at the K.G.Jebsen center for neuropsychiatric disorders at the University of Bergen, Norway, was supported by Stiftelsen K.G. Jebsen, European Community’s Seventh Framework Program under grant agreement no 602805 and the H2020 Research and Innovation Program under grant agreement numbers 643051 and 667302. HUNT: The HUNT Study is a collaboration between HUNT Research Centre (Faculty of Medicine, Norwegian University of Science and Technology), Nord-Trøndelag County Council, Central Norway Health Authority, and the Norwegian Institute of Public Health. HUNT-MRI was funded by the Liaison Committee between the Central Norway Regional Health Authority and the Norwegian University of Science and Technology, and the Norwegian National Advisory Unit for functional MRI. IMAGEN: The work received support from the European Union-funded FP6Integrated Project IMAGEN (Reinforcement-related behaviour in normal brain function and psychopathology) (LSHM-CT- 2007-037286), the Horizon 2020 funded ERC Advanced Grant ‘STRATIFY’ (Brain network based stratification of reinforcement-related disorders) (695313), ERANID (Understanding the Interplay between Cultural, Biological and Subjective Factors in Drug Use Pathways) (PR-ST-0416-10004), BRIDGET (JPND: BRain Imaging, cognition Dementia and next generation GEnomics) (MR/N027558/1), the FP7 projects IMAGEMEND (602450; IMAging GEnetics for MENtal Disorders) and MATRICS (603016), the Innovative Medicine Initiative Project EU-AIMS (115300), the Medical Research Council Grant ‘c-VEDA’ (Consortium on Vulnerability to Externalizing Disorders and Addictions) (MR/N000390/1), the Swedish Research Council FORMAS, the Medical Research Council, the National Institute for Health Research (NIHR) Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King’s College London, the Bundesministeriumfür Bildung und Forschung (BMBF grants 01GS08152; 01EV0711; eMED SysAlc01ZX1311A; Forschungsnetz AERIAL), the Deutsche Forschungsgemeinschaft (DFG grants SM 80/7-1, SM 80/7-2, SFB 940/1). Further support was provided by grants from: ANR (project AF12-NEUR0008-01—WM2NA, and ANR-12-SAMA-0004), the Fondation de France, the Fondation pour la Recherche Médicale, the Mission Interministérielle de Lutte-contre-les-Drogues-et-les-Conduites-Addictives (MILDECA), the Assistance-Publique-Hôpitaux-de-Paris and INSERM (interface grant), Paris Sud University IDEX 2012; the National Institutes of Health, Science Foundation Ireland (16/ERCD/3797), USA (Axon, Testosterone and Mental Health during Adolescence; RO1 MH085772-01A1), and by NIH Consortium grant U54 EB020403, supported by a cross-NIH alliance that funds Big Data to Knowledge Centres of Excellence. MCIC: This work was supported primarily by the Department of Energy DE-FG02-99ER62764 through its support of the Mind Research Network and the consortium as well as by the National Association for Research in Schizophrenia and Affective Disorders (NARSAD) Young Investigator Award (to SE) as well as through the Blowitz-Ridgeway and Essel Foundations, and through NWO ZonMw TOP 91211021, the DFG research fellowship (to SE), the Mind Research Network, National Institutes of Health through NCRR 5 month-RR001066 (MGH General Clinical Research Center), NIMH K08 MH068540, the Biomedical Informatics Research Network with NCRR Supplements to P41 RR14075 (MGH), M01 RR 01066 (MGH), NIBIB R01EB006841 (MRN), R01EB005846 (MRN), 2R01 EB000840 (MRN), 1RC1MH089257 (MRN), as well as grant U24 RR021992. NCNG: this sample collection was supported by grants from the Bergen Research Foundation and the University of Bergen, the Dr Einar Martens Fund, the K.G. Jebsen Foundation, the Research Council of Norway, to SLH, VMS and TE. The Bergen group was supported by grants from the Western Norway Regional Health Authority (Grant 911593 to AL, Grant 911397 and 911687 to AJL). NESDA: Funding for NESDA was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10-000-1002); the Center for Medical Systems Biology (CSMB, NWO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University’s Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, R01D0042157-01A, MH081802, Grand Opportunity grants 1RC2 MH089951 and 1RC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health.Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. NTR: The NTR study was supported by the Netherlands Organization for Scientific Research (NWO), MW904-61-193 (Eco de Geus & Dorret Boomsma), MaGW-nr: 400-07- 080 (Dennis van ‘t Ent), MagW 480-04-004 (Dorret Boomsma), NWO/SPI 56-464-14192 (Dorret Boomsma), the European Research Council, ERC-230374 (Dorret Boomsma), and Amsterdam Neuroscience. OATS: OATS (Older Australian Twins Study) was facilitated by access to Twins Research Australia, which is funded by a National Health and Medical Research Council (NHMRC) Enabling Grant 310667. OATS is also supported via a NHMRC/Australian Research Council Strategic Award (401162) and a NHMRC Project Grant (1045325). DNA extraction was performed by Genetic Repositories Australia, which was funded by a NHMRC Enabling Grant (401184). OATS genotyping was partly funded by a Commonwealth Scientific and Industrial Research Organisation Flagship Collaboration Fund Grant. PAFIP: PAFIP data were collected at the Hospital Universitario Marqués de Valdecilla, University of Cantabria, Santander, Spain, under the following grant support: Carlos III Health Institute PIE14/00031 and SAF2013-46292-R and SAF2015-71526-REDT. We wish to acknowledge IDIVAL Neuroimaging Unit for imaging acquirement and analysis.We want to particularly acknowledge the patients and the BioBankValdecilla (PT13/0010/0024) integrated in the Spanish National Biobanks Network for its collaboration. QTIM: The QTIM study was supported by grants from the US National Institute of Child Health and Human Development (R01 HD050735) and the Australian National Health and Medical Research Council (NHMRC) (486682, 1009064). Genotyping was supported by NHMRC (389875). Lachlan Strike is supported by an Australian Postgraduate Award (APA). AFM is supported by NHMRC CDF 1083656. We thank the twins and siblings for their participation, the many research assistants, as well as the radiographers, for their contribution to data collection and processing of the samples. SHIP: SHIP is part of the Community Medicine Research net of the University of Greifswald, Germany, which is funded by the Federal Ministry of Education and Research (grants no. 01ZZ9603, 01ZZ0103, 01ZZ0403 and 01ZZ0701), the Ministry of Cultural Affairs as well as the Social Ministry of the Federal State of Mecklenburg-West Pomerania, and the network ‘Greifswald Approach to Individualized Medicine (GANI_MED)’ funded by the Federal Ministry of Education and Research (grant 03IS2061A). Genome-wide data have been supported by the Federal Ministry of Education and Research (grant no. 03ZIK012) and a joint grant from Siemens Healthineers, Erlangen, Germany and the Federal State of Mecklenburg- West Pomerania. Whole-body MR imaging was supported by a joint grant from Siemens Healthineers, Erlangen, Germany and the Federal State of Mecklenburg West Pomerania. The University of Greifswald is a member of the Caché Campus program of the InterSystems GmbH. StrokeMRI: StrokeMRI has been supported by the Research Council of Norway (249795), the South-Eastern Norway Regional Health Authority (2014097, 2015044, 2015073) and the Norwegian ExtraFoundation for Health and Rehabilitation. TOP: TOP is supported by the Research Council of Norway (223273, 213837, 249711), the South East Norway Health Authority (2017-112), the Kristian Gerhard Jebsen Stiftelsen (SKGJ‐MED‐008) and the European Community’s Seventh Framework Programme (FP7/2007–2013), grant agreement no. 602450 (IMAGEMEND). We acknowledge the technical support and service from the Genomics Core Facility at the Department of Clinical Science, the University of Bergen

BackgroundCpG methylation levels can help to explain inter-individual differences in phenotypic traits. Few studies have explored whether identifying CpG subsets based on biological and statistical properties can maximise predictions while minimising array content.MethodsVariance component analyses and penalised regression (epigenetic predictors) were used to test the influence of (i) the number of CpGs considered, (ii) mean CpG methylation variability and (iii) methylation QTL status on the variance captured in eighteen traits by blood DNA methylation. Training and test sets comprised ≤4,450 and ≤2,578 unrelated individuals from Generation Scotland, respectively.ResultsAs the number of CpG sites under consideration decreased, so too did the estimates from the variance components and prediction analyses. Methylation QTL status and mean CpG variability did not influence variance components. However, relative effect sizes were 15% larger for epigenetic predictors based on CpGs with methylation QTLs compared to sites without methylation QTLs. Relative effect sizes were 45% larger for predictors based on CpGs with mean beta-values between 10%-90% compared to those using hypo- or hypermethylated CpGs (beta-value ≤10% or ≥90%).ConclusionArrays with fewer CpGs 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 CpGs in array content selection.

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 137 connections between EpiScores for proteins and a future diagnosis of common adverse health outcomes. These included diabetes, stroke, depression, Alzheimer’s dementia, 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.

Ankylosing Spondylitis (AS) is a highly heritable inflammatory arthritis which occurs more frequently in men than women. In their recent publication examining sex differences in the genetic aetiology of common complex traits and diseases, Bernabeu et al. (2021) observe differences in heritability of AS between sexes, and a genome-wide significant genotype by sex interaction in risk of AS at the major histocompatability (MHC) locus1. The authors then present evidence suggesting that this genotype by sex interaction arises primarily as a result of differential expression of the gene MICA across the sexes in skeletal muscle tissue. Through a series of conditional association analyses in the UK Biobank, reanalysis of the GTEx gene expression resource and RNASeq experiments on peripheral blood cells from AS cases and controls, we show that the genotype by sex interaction the authors’ report is unlikely to be a result of variation in MICA, but probably reflects a known interaction between the HLA-B gene, sex and risk of AS. We demonstrate that the diagnostic accuracy of AS in the UK Biobank is low, particularly amongst women, likely explaining some of the observed differences in heritability across the sexes and the difficulty in precisely locating association signals in the cohort.

The mechanisms underlying how monozygotic (or identical) twins arise are yet to be determined. Here, the authors investigate this in an epigenome-wide association study, showing that monozygotic twinning has a characteristic DNA methylation signature in adult somatic tissues. Monozygotic (MZ) twins and higher-order multiples arise when a zygote splits during pre-implantation stages of development. The mechanisms underpinning this event have remained a mystery. Because MZ twinning rarely runs in families, the leading hypothesis is that it occurs at random. Here, we show that MZ twinning is strongly associated with a stable DNA methylation signature in adult somatic tissues. This signature spans regions near telomeres and centromeres, Polycomb-repressed regions and heterochromatin, genes involved in cell-adhesion, WNT signaling, cell fate, and putative human metastable epialleles. Our study also demonstrates a never-anticipated corollary: because identical twins keep a lifelong molecular signature, we can retrospectively diagnose if a person was conceived as monozygotic twin., Netherlands Organization for Scientific Research (NWO); Biobanking and Biomolecular Research Infrastructure (BBMRI–NL); NWO Large Scale infrastructures; X-Omics