Search

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.

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.

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.

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

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

The identification of biomarkers that discriminate individual ageing trajectories is a principal target of ageing research. Some of the most promising predictors of biological ageing have been developed using DNA methylation. One recent candidate, which tracks age-related phenotypes in addition to chronological age, is ‘DNAm PhenoAge’. Here, we performed a phenome-wide association analysis of this biomarker in a cohort of older adults to assess its relationship with a comprehensive set of both historical, and contemporaneously-measured, phenotypes. Higher than expected DNAm PhenoAge compared with chronological age, known as epigenetic age acceleration, was found to associate with a number of blood, cognitive, physical fitness and lifestyle variables, and with mortality. Notably, DNAm PhenoAge, assessed at age 70, was associated with cognitive ability at age 11, and with educational attainment. Adjusting for age 11 cognitive ability attenuated the majority of the cross-sectional later-life associations between DNAm PhenoAge and health outcomes. These results highlight the importance of early life factors on healthy older ageing.

Background DNA methylation changes with age. Chronological age predictors built from DNA methylation are termed ‘epigenetic clocks’. The deviation of predicted age from the actual age (‘age acceleration residual’, AAR) has been reported to be associated with death. However, it is currently unclear how a better prediction of chronological age affects such association. Methods In this study, we build multiple predictors based on training DNA methylation samples selected from 13,661 samples (13,402 from blood and 259 from saliva). We use the Lothian Birth Cohorts of 1921 (LBC1921) and 1936 (LBC1936) to examine whether the association between AAR (from these predictors) and death is affected by (1) improving prediction accuracy of an age predictor as its training sample size increases (from 335 to 12,710) and (2) additionally correcting for confounders (i.e., cellular compositions). In addition, we investigated the performance of our predictor in non-blood tissues. Results We found that in principle, a near-perfect age predictor could be developed when the training sample size is sufficiently large. The association between AAR and mortality attenuates as prediction accuracy increases. AAR from our best predictor (based on Elastic Net, https://github.com/qzhang314/DNAm-based-age-predictor) exhibits no association with mortality in both LBC1921 (hazard ratio = 1.08, 95% CI 0.91–1.27) and LBC1936 (hazard ratio = 1.00, 95% CI 0.79–1.28). Predictors based on small sample size are prone to confounding by cellular compositions relative to those from large sample size. We observed comparable performance of our predictor in non-blood tissues with a multi-tissue-based predictor. Conclusions This study indicates that the epigenetic clock can be improved by increasing the training sample size and that its association with mortality attenuates with increased prediction of chronological age. Electronic supplementary material The online version of this article (10.1186/s13073-019-0667-1) contains supplementary material, which is available to authorized users.

Introduction: The gut microbiota has important roles in maintaining human health. The microbiota and its metabolic byproducts could play a role in the pathogenesis of neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS). Areas covered: The authors evaluate the methods of assessing the gut microbiota, and also review how the gut microbiota affects the various physiological functions of the gut. The authors then consider how gut dysbiosis could theoretically affect the pathogenesis of ALS. They present the current evidence regarding the composition of the gut microbiota in ALS and in rodent models of ALS. Finally, the authors review therapies that could improve gut dysbiosis in the context of ALS. Expert opinion: Currently reported studies suggest some instances of gut dysbiosis in ALS patients and mouse models; however, these studies are limited, and more information with well-controlled larger datasets is required to make a definitive judgment about the role of the gut microbiota in ALS pathogenesis. Overall this is an emerging field that is worthy of further investigation. The authors advocate for larger studies using modern metagenomic techniques to address the current knowledge gaps.

Despite extensive sex differences in human complex traits and disease, the male and female genomes differ only in the sex chromosomes. This implies that most sex-differentiated traits are the result of differences in the expression of genes that are common to both sexes. While sex differences in gene expression have been observed in a range of different tissues, the biological mechanisms for tissue-specific sex differences (TSSDs) in gene expression are not well understood. A total of 30 640 autosomal and 1021 X-linked transcripts were tested for heterogeneity in sex difference effect sizes in n = 617 individuals across 40 tissue types in Genotype–Tissue Expression (GTEx). This identified 65 autosomal and 66 X-linked TSSD transcripts (corresponding to unique genes) at a stringent significance threshold. Results for X-linked TSSD transcripts showed mainly concordant direction of sex differences across tissues and replicate previous findings. Autosomal TSSD transcripts had mainly discordant direction of sex differences across tissues. The top cis-expression quantitative trait loci (eQTLs) across tissues for autosomal TSSD transcripts are located a similar distance away from the nearest androgen and estrogen binding motifs and the nearest enhancer, as compared to cis-eQTLs for transcripts with stable sex differences in gene expression across tissue types. Enhancer regions that overlap top cis-eQTLs for TSSD transcripts, however, were found to be more dispersed across tissues. These observations suggest that androgen and estrogen regulatory elements in a cis region may play a common role in sex differences in gene expression, but TSSD in gene expression may additionally be due to causal variants located in tissue-specific enhancer regions.

ENIGMA-CNV working group., Low-frequency 1q21.1 distal deletion and duplication copy number variant (CNV) carriers are predisposed to multiple neurodevelopmental disorders, including schizophrenia, autism and intellectual disability. Human carriers display a high prevalence of micro- and macrocephaly in deletion and duplication carriers, respectively. The underlying brain structural diversity remains largely unknown. We systematically called CNVs in 38 cohorts from the large-scale ENIGMA-CNV collaboration and the UK Biobank and identified 28 1q21.1 distal deletion and 22 duplication carriers and 37,088 non-carriers (48% male) derived from 15 distinct magnetic resonance imaging scanner sites. With standardized methods, we compared subcortical and cortical brain measures (all) and cognitive performance (UK Biobank only) between carrier groups also testing for mediation of brain structure on cognition. We identified positive dosage effects of copy number on intracranial volume (ICV) and total cortical surface area, with the largest effects in frontal and cingulate cortices, and negative dosage effects on caudate and hippocampal volumes. The carriers displayed distinct cognitive deficit profiles in cognitive tasks from the UK Biobank with intermediate decreases in duplication carriers and somewhat larger in deletion carriers—the latter potentially mediated by ICV or cortical surface area. These results shed light on pathobiological mechanisms of neurodevelopmental disorders, by demonstrating gene dose effect on specific brain structures and effect on cognitive function., 1000BRAINS: The 1000BRAINS study was funded by the Institute of Neuroscience and Medicine, Research Center Juelich, Germany. We thank the Heinz Nixdorf Foundation (Germany) for the generous support of the Heinz Nixdorf Recall Study on which 1000BRAINS is based. We also thank the scientists and the study staff of the Heinz Nixdorf Recall Study and 1000BRAINS. Funding was also granted by the Initiative and Networking Fund of the Helmholtz Association (Caspers) and the European Union’s Horizon 2020 Research and Innovation Programme under Grant Agreement 945539 (Human Brain Project SGA3; Amunts, Caspers, Cichon). Brainscale: The Brainscale study was supported by the Netherlands Organization for Scientific Research MagW 480-04-004 (Dorret I. Boomsma), 51.02.060 (Hilleke E. Hulshoff Pol), 668.772 (Dorret I. Boomsma and Hilleke E. Hulshoff Pol); NWO/SPI 56-464-14192 (Dorret I. Boomsma), the European Research Council (ERC-230374) (Dorret I. Boomsma), High Potential Grant Utrecht University (Hilleke E.Hulshoff Pol) and NWO Brain and Cognition 433-09-220 (Hilleke E.Hulshoff Pol). Betula: The Betula study was funded by the Knut and Alice Wallenberg (KAW) foundation (Nyberg). The Freesurfer segmentations on the Betula sample were performed on resources provided by the Swedish National Infrastructure for Computing (SNIC) at HPC2N (in Umeå, Sweden), partially funded by the Swedish Research Council through grant agreement no. 2018-05973. 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 from 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 support from the participating departments and centres and from external grants, that is, the Biobanking and Biomolecular Resources Research Infrastructure (Netherlands) (BBMRI-NL), the Hersenstichting Nederland and the Netherlands Organization 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 agreement nos. 602450 (IMAGEMEND), 278948 (TACTICS) and 602805 (Aggressotype), as well as from the European Community’s Horizon 2020 programme under grant agreement no. 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. deCODE genetics: deCODE genetics acknowledges support from the Innovative Medicines Initiative Joint Undertaking under grant agreement nos. 115008 (NEWMEDS) and 115300 (EUAIMS), of which resources are composed of EFPIA in-kind contribution and financial contribution from the European Union’s Seventh Framework Programme (EU-FP7/2007-2013), EU-FP7-funded grant agreement no. 602450 (IMAGEMEND) and EU-funded FP7-People-2011-IAPP grant agreement no. 286213 (PsychDPC). Dublin: This work was supported by Science Foundation Ireland (SFI grant 12/IP/1359 to Gary Donohoe and grant SFI08/IN.1/B1916-Corvin to Aidan C. Corvin). ECHO-DEFINE: The ECHO study acknowledges funding from a Medical Research Council (MRC) Centre Grant to Michael J. Owen (G0801418), the Wellcome Trust (Institutional Strategic Support Fund (ISSF) to van den Bree and Clinical Research Training Fellowship to Joanne L. Doherty), the Waterloo Foundation (WF 918-1234 to van den Bree), the Baily Thomas Charitable Fund (2315/1 to van den Bree), National Institute of Mental Health (NIMH 5UO1MH101724 to van den Bree and Michael J. Owen), the IMAGINE-2 study (funded by the MRC (MR/T033045/1) to van den Bree, Jeremy Hall and Michael J. Owen), the IMAGINE-ID study (funded by MRC (MR/N022572/1) to Jeremy Hall, van den Bree and Owen). The DEFINE study was supported by a Wellcome Trust Strategic Award (100202/Z/12/Z) to Michael J. Owen. ENIGMA: ENIGMA is supported in part by NIH grants U54 EB20403, R01MH116147 and R56AG058854. NIA T32AG058507; NIH/NIMH 5T32MH073526. EPIGEN-Dublin: The EPIGEN-Dublin cohort was supported by a Science Foundation Ireland Research Frontiers Programme award (08/RFP/GEN1538). EPIGEN-UK (Sisodiya): 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. GAP: This work was supported by the National Institute for Health Research (NIHR) Mental Health Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and the Institute of Psychiatry, Psychology and Neuroscience, King’s College London. The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health. GOBS: The GOBS study data collection was supported in part by the National Institutes of Health (NIH) grants: R01 MH078143, R01 MH078111 and R01 MH083824, with work conducted in part in facilities constructed under the support of NIH grant C06 RR020547. GSP: Data were in part provided by the Brain Genomics Superstruct Project (GSP) of Harvard University and Massachusetts General Hospital (MGH) (Principal Investigators: Randy Buckner, Jordan Smoller and Joshua Roffman), with support from the Center for Brain Science Neuroinformatics Research Group, Athinoula A. Martinos Center for Biomedical Imaging, Center for Genomic Medicine and Stanley Center for Psychiatric Research. Twenty individual investigators at Harvard and MGH generously contributed data to the overall project. We would like to thank Randy Buckner for insightful comments and feedback on this work. HUBIN: The HUBIN study was financed by the Swedish Research Council (K2010-62X-15078-07-2, K2012-61X-15078-09-3, 521-2014-3487 K2015-62X-15077-12-3, 2017-00949), the regional agreement on medical training and clinical research between Stockholm County Council and the Karolinska Institutet. HUNT: The HUNT study is a collaboration between HUNT Research Centre (Faculty of Medicine and Movement Sciences, NTNU—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: This work received support from the following sources: the European Union-funded FP6 Integrated 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), Human Brain Project (HBP SGA 2, 785907),the FP7 projects IMAGEMEND(602450; IMAging GEnetics for MENtal Disorders) and MATRICS (603016), the Innovative Medicine Initiative Project EUAIMS (115300-2), 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 01EE1406A, 01EE1406B), the Deutsche Forschungsgemeinschaft (DFG grants, SM 80/7-2, SFB 940/2), the Medical Research Foundation and Medical Research Council (grants MR/R00465X/1 and MR/S020306/1). Further support was provided by grants from: ANR (project AF12-NEUR0008-01—WM2NA, ANR-12-SAMA-0004), the Eranet Neuron (ANR-18-NEUR00002-01), the Fondation de France (00081242), the Fondation pour la Recherche Médicale (DPA20140629802), 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 Fondation de l’Avenir (grant AP-RM-17-013), the Fédération pour la Recherche sur le Cerveau; 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. Lifespan: The study is funded by the Research Council of Norway (230345, 288083 and 223273). NCNG: NCNG sample collection was supported by grants from the Bergen Research Foundation and the University of Bergen, the Dr Einar Martens Fund, the Research Council of Norway, to le Hellard, Steen and Espeseth. The Bergen group was supported by grants from the Western Norway Regional Health Authority (Grant 911593 to Arvid Lundervold, Grant 911397 and 911687 to Astri Johansen Lundervold). NTR: The NTR cohort was supported by the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organisation for Health Research and Development (ZonMW) grants 904-61-090, 985-10-002, 912-10-020, 904-61-193, 480-04-004,463-06-001, 451-04-034, 400-05-717, Addiction-31160008, 016-115-035, 481-08-011, 056-32-010, Middelgroot-911-09-032, OCW_NWO Gravity programme—024.001.003, NWO-Groot 480-15-001/674, Center for Medical Systems Biology (CSMB, NWO Genomics), NBIC/BioAssist/RK(2008.024), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL, 184.021.007 and 184.033.111); Spinozapremie (NWO-56-464-14192), KNAW Academy Professor Award (PAH/6635) and University Research Fellow grant (URF) to Dorret I. Boomsma; Amsterdam Public Health research institute (former EMGO+), Neuroscience Amsterdam research institute (former NCA); the European Science Foundation (ESF, EU/QLRT-2001-01254), the European Community’s Seventh Framework Programme (FP7- HEALTH-F4-2007-2013, grant 01413: ENGAGE and grant 602768: ACTION); the European Research Council (ERC Starting 284167, ERC Consolidator 771057, ERC Advanced 230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the National Institutes of Health (NIH, R01D0042157-01A1, R01MH58799-03, MH081802, DA018673, R01 DK092127-04, Grand Opportunity grants 1RC2 MH089951 and 1RC2 MH089995); the Avera Institute for Human Genetics, Sioux Falls, South Dakota (USA). 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 NWO through grant 2018/EW/00408559, BiG Grid, the Dutch e-Science Grid and SURFSARA. OATS: The OATS study has been funded by a National Health & Medical Research Council (NHMRC) and Australian Research Council (ARC) Strategic Award Grant of the Ageing Well, Ageing Productively Programme (ID No. 401162) and NHMRC Project Grants (ID Nos. 1045325 and 1085606). This research was facilitated through Twins Research Australia, a national resource in part supported by an NHMRC Centre for Research Excellence Grant (ID No.: 1079102). We thank the participants for their time and generosity in contributing to this research. We acknowledge the contribution of the OATS research team (https://cheba.unsw.edu.au/project/older-australian-twins-study) to this study. OATS genotyping was partly funded by a Commonwealth Scientific and Industrial Research Organization Flagship Collaboration Fund Grant. Osaka: Osaka study was supported by the Brain Mapping by Integrated Neurotechnologies for Disease Studies (Brain/MINDS: Grant Number JP18dm0207006), Brain/MINDS& beyond studies (Grant Number JP20dm0307002) and Health and Labour Sciences Research Grants for Comprehensive Research on Persons with Disabilities (Grant Number JP20dk0307081) from the Japan Agency for Medical Research and Development (AMED), Grants-in-Aid for Scientific Research (KAKENHI; Grant Numbers JP25293250 and JP16H05375). Some computations were performed at the Research Center for Computational Science, Okazaki, Japan. PAFIP: The PAFIP study was supported by Instituto de Salud Carlos III, FIS 00/3095, 01/3129, PI020499, PI060507, PI10/00183, the SENY Fundació Research Grant CI2005-0308007 and the FundaciónMarqués de Valdecilla API07/011. Biological samples from our cohort were stored at the Valdecilla Biobank and genotyping services were conducted at the Spanish ‘Centro Nacional de Genotipado’ (CEGEN-ISCIII). MCIC/COBRE: The study is funded by the National Institutes of Health studies R01EB006841, P20GM103472 and P30GM122734 and Department of Energy DE-FG02-99ER62764. PING: Data collection and sharing for the Paediatric Imaging, Neurocognition and Genetics (PING) Study (National Institutes of Health Grant RC2DA029475) were funded by the National Institute on Drug Abuse and the Eunice Kennedy Shriver National Institute of Child Health & Human Development. A full list of PING investigators is at http://pingstudy.ucsd.edu/investigators.html. QTIM: The QTIM study was supported by the National Institute of Child Health and Human Development (R01 HD050735) and the National Health and Medical Research Council (NHMRC 486682, 1009064), Australia. Genotyping was supported by NHMRC (389875). Medland is supported in part by an NHMRC fellowship (APP1103623). 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 (grant nos. 01ZZ9603, 01ZZ0103 and 01ZZ0403), the Ministry of Cultural Affairs and the Social Ministry of the Federal State of Mecklenburg-West Pomerania. Genome-wide single-nucleotide polymorphism typing in SHIP and MRI scans in SHIP and SHIP-TREND have been supported by a joint grant from Siemens Healthineers, Erlangen, Germany and the Federal State of Mecklenburg-West Pomerania. StrokeMRI: StrokeMRI was supported by the Norwegian ExtraFoundation for Health and Rehabilitation(2015/FO5146), the Research Council of Norway (249795, 262372), the South-Eastern Norway Regional Health Authority (2014097, 2015044, 2015073) and the Department of Psychology, University of Oslo. Sydney MAS: The Sydney Memory and Aging Study (Sydney MAS) is funded by National and HealthMedical Research Council (NHMRC) Programme and Project Grants (ID350833, ID568969 and ID109308). We also thank the Sydney MAS participants and the Research Team. SYS: The SYS Study is supported by Canadian Institutes of Health Research. TOP: Centre of Excellence: RCN #23273 and RCN #226971. Part of this work was performed on the TSD (Tjeneste for Sensitive Data) facilities, owned by the University of Oslo, operated and developed by the TSD service group at the University of Oslo, IT-Department (USIT) (tsd-drift@usit.uio.no). The research leading to these results has received funding from the European Union Seventh Framework Programme (FP7-PEOPLE-2013-COFUND) under grant agreement no. 609020—Scientia Fellows; the Research Council of Norway (RCN) #276082—A lifespan perspective on mental illness: toward precision medicine using multimodal brain imaging and genetics. Ida E. Sønderby and Rune Bøen are supported by South-Eastern Norway Regional Health Authority (#2020060). Ida E. Sønderby and Ole A. Andreassen have received funding from the European Union’s Horizon 2020 Research and Innovation Programme under Grant agreement no. 847776 (CoMorMent project) and the KG Jebsen Foundation (SKGJ-MED-021). UCLA_UMCU: The UCLA_UMCU cohort comprises of six studies which were supported by National Alliance for Research in Schizophrenia and Affective Disorders (NARSAD) (20244 to Prof. Hillegers), The Netherlands Organisation for Health Research and Development (ZonMw) (908-02-123 to Prof. Hulshoff Pol), and Netherlands Organisation for Scientific Research (NWO 9120818 and NWO-VIDI 917-46-370 to Prof. Hulshoff Pol). The GROUP study was funded through the Geestkracht programme of the Dutch Health Research Council (ZonMw, grant number 10-000-1001), and matching funds from participating pharmaceutical companies (Lundbeck, AstraZeneca, Eli Lilly and Janssen Cilag) and universities and mental health care organizations (Amsterdam: Academic Psychiatric Centre of the Academic Medical Center and the mental health institutions: GGZ inGeest, Arkin, Dijk en Duin, GGZ Rivierduinen, Erasmus Medical Centre, GGZ Noord-Holland-Noord. Groningen: University Medical Center Groningen and the mental health institutions: Lentis, GGZ Friesland, GGZ Drenthe, Dimence, Mediant, GGNet Warnsveld, Yulius Dordrecht and Parnassia Psycho-medical Center, The Hague. Maastricht: Maastricht University Medical Centre and the mental health institutions: GGzE, GGZ Breburg, GGZ Oost-Brabant, Vincent van Gogh, voor Geestelijke Gezondheid, Mondriaan, Virenzeriagg, Zuyderland GGZ, MET ggz, Universitair Centrum Sint-JozefKortenberg, CAPRI University of Antwerp, PC Ziekeren Sint-Truiden, PZ Sancta Maria Sint-Truiden, GGZ Overpelt, OPZ Rekem. Utrecht: University Medical Center Utrecht and the mental health institutions: Altrecht, GGZ Centraal and Delta.). UK Biobank: This work made use of data sharing from UK Biobank (under project code 27412). Others: Work by Pierre Vanderhaeghen was funded by Grants of the European Research Council (ERC Adv Grant GENDEVOCORTEX), the EOS Programme, the Belgian FWO, the AXA Research Fund and the Belgian Queen Elizabeth Foundation. Ikuo K. Suzuki was supported by a postdoctoral fellowship of the FRS/FNRS.

DNA methylation quantitative trait locus (mQTL) analyses on 32,851 participants identify genetic variants associated with DNA methylation at 420,509 sites in blood, resulting in a database of >270,000 independent mQTLs. 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.

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. Here we describe results of DNA methylation-quantitative trait loci (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 mQTL of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We reveal that the genetic architecture of DNAm levels is highly polygenic and DNAm exhibits signatures of negative and positive natural selection. Using shared genetic control between distal DNAm sites we construct networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic factors are associated with both blood DNAm levels and complex diseases but in most cases these associations do not reflect causal relationships from DNAm to trait or vice versa indicating a more complex genotype-phenotype map than has previously been hypothesised.

Context Circulating concentrations of free triiodothyronine (fT3), free thyroxine (fT4), and thyrotropin (TSH) are partly heritable traits. Recent studies have advanced knowledge of their genetic architecture. Epigenetic modifications, such as DNA methylation (DNAm), may be important in pituitary-thyroid axis regulation and action, but data are limited. Objective To identify novel associations between fT3, fT4, and TSH and differentially methylated positions (DMPs) in the genome in subjects from 2 Australian cohorts. Method We performed an epigenome-wide association study (EWAS) of thyroid function parameters and DNAm using participants from: Brisbane Systems Genetics Study (median age 14.2 years, n = 563) and the Raine Study (median age 17.0 years, n = 863). Plasma fT3, fT4, and TSH were measured by immunoassay. DNAm levels in blood were assessed using Illumina HumanMethylation450 BeadChip arrays. Analyses employed generalized linear mixed models to test association between DNAm and thyroid function parameters. Data from the 2 cohorts were meta-analyzed. Results We identified 2 DMPs with epigenome-wide significant (P Conclusions This study has demonstrated associations between blood-based DNAm and both fT3 and TSH. This may provide insight into mechanisms underlying thyroid hormone action and/or pituitary-thyroid axis function.

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; R2 = 4.4%, p = 2.1 × 10−5). 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−7; DNAm IL-6 score: N = 7028, β = 0.02, SE = 0.0009, p < 2 × 10−16). 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, pFDR < 2 × 10−16). 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.

Promoter-anchored chromatin interactions (PAIs) play a pivotal role in transcriptional regulation. Current high-throughput technologies for detecting PAIs, such as promoter capture Hi-C, are not scalable to large cohorts. Here, we present an analytical approach that uses summary-level data from cohort-based DNA methylation (DNAm) quantitative trait locus (mQTL) studies to predict PAIs. Using mQTL data from human peripheral blood (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n = 1980$$\end{document}n=1980), we predict 34,797 PAIs which show strong overlap with the chromatin contacts identified by previous experimental assays. The promoter-interacting DNAm sites are enriched in enhancers or near expression QTLs. Genes whose promoters are involved in PAIs are more actively expressed, and gene pairs with promoter-promoter interactions are enriched for co-expression. Integration of the predicted PAIs with GWAS data highlight interactions among 601 DNAm sites associated with 15 complex traits. This study demonstrates the use of mQTL data to predict PAIs and provides insights into the role of PAIs in complex trait variation., Promoter-anchored chromatin interaction (PAI) is a mechanism by which gene expression is regulated, but methods to measure PAIs are costly and currently not scalable. Here, the authors develop an approach by which PAIs can be predicted using summary-level data from methylation QTL studies.

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

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. Therefore, 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). We identified 13 SNPs associated with 13 proteins (n = 1 SNP each) concordant across OLS and Bayesian methods. We identified three CpG sites spread across three 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 one 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. Our data may aid understanding of the molecular regulation of the circulating inflammatory proteome as well as causal relationships between inflammatory mediators and disease.

We conducted DNA methylation association analyses using Illumina 450K data from whole blood for an Australian amyotrophic lateral sclerosis (ALS) case–control cohort (782 cases and 613 controls). Analyses used mixed linear models as implemented in the OSCA software. We found a significantly higher proportion of neutrophils in cases compared to controls which replicated in an independent cohort from the Netherlands (1159 cases and 637 controls). The OSCA MOMENT linear mixed model has been shown in simulations to best account for confounders. When combined in a methylation profile score, the 25 most-associated probes identified by MOMENT significantly classified case–control status in the Netherlands sample (area under the curve, AUC = 0.65, CI95% = [0.62–0.68], p = 8.3 × 10−22). The maximum AUC achieved was 0.69 (CI95% = [0.66–0.71], p = 4.3 × 10−34) when cell-type proportion was included in the predictor.

Background Advanced age is associated with cognitive and physical decline and is a major risk factor for a multitude of disorders. There is also a gap in life expectancy between males and females. DNA methylation differences have been shown to be associated with both age and sex. Here, we investigate age-by-sex differences in blood-based DNA methylation in an unrelated cohort of 2586 individuals between the ages of 18 and 87 years, with replication in a further 4450 individuals between the ages of 18 and 93 years. Methods Linear regression models were applied, with stringent genome-wide significance thresholds (p −8) used in both the discovery and replication data. A second, highly conservative mixed linear model method that better controls the false-positive rate was also applied, using the same genome-wide significance thresholds. Results Using the linear regression method, 52 autosomal and 597 X-linked CpG sites, mapping to 251 unique genes, replicated with concordant effect size directions in the age-by-sex interaction analysis. The site with the greatest difference mapped to GAGE10, an X-linked gene. Here, DNA methylation levels remained stable across the male adult age range (DNA methylation by age r = 0.02) but decreased across female adult age range (DNA methylation by age r = − 0.61). One site (cg23722529) with a significant age-by-sex interaction also had a quantitative trait locus (rs17321482) that is a genome-wide significant variant for prostate cancer. The mixed linear model method identified 11 CpG sites associated with the age-by-sex interaction. Conclusion The majority of differences in age-associated DNA methylation trajectories between sexes are present on the X chromosome. Several of these differences occur within genes that have been implicated in sexually dimorphic traits.

DNA methylation, which is modulated by both genetic factors and environmental exposures, may offer a unique opportunity to discover novel biomarkers of disease-related brain phenotypes, even when measured in other tissues than brain, such as blood. A few studies of small sample sizes have revealed associations between blood DNA methylation and neuropsychopathology, however, large-scale epigenome-wide association studies (EWAS) are needed to investigate the utility of DNA methylation profiling as a peripheral marker for the brain. Here, in an analysis of eleven international cohorts, totalling 3337 individuals, we report epigenome-wide meta-analyses of blood DNA methylation with volumes of the hippocampus, thalamus and nucleus accumbens (NAcc)—three subcortical regions selected for their associations with disease and heritability and volumetric variability. Analyses of individual CpGs revealed genome-wide significant associations with hippocampal volume at two loci. No significant associations were found for analyses of thalamus and nucleus accumbens volumes. Cluster-based analyses revealed additional differentially methylated regions (DMRs) associated with hippocampal volume. DNA methylation at these loci affected expression of proximal genes involved in learning and memory, stem cell maintenance and differentiation, fatty acid metabolism and type-2 diabetes. These DNA methylation marks, their interaction with genetic variants and their impact on gene expression offer new insights into the relationship between epigenetic variation and brain structure and may provide the basis for biomarker discovery in neurodegeneration and neuropsychiatric conditions.

Importance: Depressive disorders arise from a combination of genetic and environmental risk factors. Epigenetic disruption provides a plausible mechanism through which gene-environment interactions lead to depression. Large-scale, epigenome-wide studies on depression are missing, hampering the identification of potentially modifiable biomarkers.Objective: To identify epigenetic mechanisms underlying depression in middle-aged and elderly persons, using DNA methylation in blood.Design, Setting, and Participants: To date, the first cross-ethnic meta-analysis of epigenome-wide association studies (EWAS) within the framework of the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium was conducted. The discovery EWAS included 7948 individuals of European origin from 9 population-based cohorts. Participants who were assessed for both depressive symptoms and whole-blood DNA methylation were included in the study. Results of EWAS were pooled using sample-size weighted meta-analysis. Replication of the top epigenetic sites was performed in 3308 individuals of African American and European origin from 2 population-based cohorts.Main Outcomes and Measures: Whole-blood DNA methylation levels were assayed with Illumina-Infinium Human Methylation 450K BeadChip and depressive symptoms were assessed by questionnaire.Results: The discovery cohorts consisted of 7948 individuals (4104 [51.6%] women) with a mean (SD) age of 65.4 (5.8) years. The replication cohort consisted of 3308 individuals (2456 [74.2%] women) with a mean (SD) age of 60.3 (6.4) years. The EWAS identified methylation of 3 CpG sites to be significantly associated with increased depressive symptoms: cg04987734 (P = 1.57 × 10−08; n = 11 256; CDC42BPB gene), cg12325605 (P = 5.24 × 10−09; n = 11 256; ARHGEF3 gene), and an intergenic CpG site cg14023999 (P = 5.99 × 10−08; n = 11 256; chromosome = 15q26.1). The predicted expression of the CDC42BPB gene in the brain (basal ganglia) (effect, 0.14; P = 2.7 × 10−03) and of ARHGEF3 in fibroblasts (effect, −0.48; P = 9.8 × 10−04) was associated with major depression.Conclusions and Relevance: This study identifies 3 methylated sites associated with depressive symptoms. All 3 findings point toward axon guidance as the common disrupted pathway in depression. The findings provide new insights into the molecular mechanisms underlying the complex pathophysiology of depression. Further research is warranted to determine the utility of these findings as biomarkers of depression and evaluate any potential role in the pathophysiology of depression and their downstream clinical effects.

Background Migraine is a common heritable neurovascular disorder typically characterised by episodic attacks of severe pulsating headache and nausea, often accompanied by visual, auditory or other sensory symptoms. Although genome-wide association studies have identified over 40 single nucleotide polymorphisms associated with migraine, there remains uncertainty about the casual genes involved in disease pathogenesis and how their function is regulated. Results We performed an epigenome-wide association study, quantifying genome-wide patterns of DNA methylation in 67 migraine cases and 67 controls with a matching age and sex distribution. Association analyses between migraine and methylation probe expression, after adjustment for cell type proportions, indicated an excess of small P values, but there was no significant single-probe association after correction for multiple testing (P

We show that genotype-by-environment interaction can be inferred from an analysis without environmental data in a large sample., Genotype-by-environment interaction (GEI) is a fundamental component in understanding complex trait variation. However, it remains challenging to identify genetic variants with GEI effects in humans largely because of the small effect sizes and the difficulty of monitoring environmental fluctuations. Here, we demonstrate that GEI can be inferred from genetic variants associated with phenotypic variability in a large sample without the need of measuring environmental factors. We performed a genome-wide variance quantitative trait locus (vQTL) analysis of ~5.6 million variants on 348,501 unrelated individuals of European ancestry for 13 quantitative traits in the UK Biobank and identified 75 significant vQTLs with P < 2.0 × 10−9 for 9 traits, especially for those related to obesity. Direct GEI analysis with five environmental factors showed that the vQTLs were strongly enriched with GEI effects. Our results indicate pervasive GEI effects for obesity-related traits and demonstrate the detection of GEI without environmental data.

The rapid increase of omic data in the past decades has greatly facilitated the investigation of associations between omic profiles such as DNA methylation (DNAm) and complex traits in large cohorts. Here, we proposed a mixed-linear-model-based method (called MOMENT) that tests for association between a DNAm probe and trait with all other distal probes fitted in multiple random-effect components to account for the effects of unobserved confounders as well as the correlations between distal probes induced by the confounders. We demonstrated by simulations that MOMENT showed a lower false positive rate and more robustness than existing methods. MOMENT has been implemented in a versatile software package (called OSCA) together with a number of other implementations for omic-data-based analysis including the estimation of variance in a trait captured by all measures of multiple omic profiles, omic-data-based quantitative trait locus (xQTL) analysis, and meta-analysis of xQTL data.

Although plasma proteins may serve as important markers of disease risk in neurological conditions, the molecular mechanisms responsible for inter-individual variation in plasma protein levels are poorly understood. In this study, we conducted genome- and epigenome-wide association studies on the levels of 92 neurological proteins to identify genetic and epigenetic loci associated with their plasma concentrations (n = 750). We identified 62 independent genome-wide significant loci for 37 proteins (P < 5.4 × 10−10) and 68 epigenome-wide significant sites associated with the levels of 7 proteins (P < 3.9 × 10−10). Using this information, we identified biological pathways in which putative neurological biomarkers are implicated as well as molecular mechanisms through which genetic variation may perturb plasma protein levels. Additionally, we found evidence that poliovirus receptor is causally associated with Alzheimer’s disease. In conclusion, we identified many novel genetic and epigenetic factors that are associated with neurological protein levels which may inform disease biology and establish causal relationships between biomarkers and neurological diseases.

Contains fulltext : 157507.pdf (Publisher’s version ) (Closed access) Spontaneous dizygotic (DZ) twinning occurs in 1%–4% of women, with familial clustering and unknown physiological pathways and genetic origin. DZ twinning might index increased fertility and has distinct health implications for mother and child. We performed a GWAS in 1,980 mothers of spontaneous DZ twins and 12,953 control subjects. Findings were replicated in a large Icelandic cohort and tested for association across a broad range of fertility traits in women. Two SNPs were identified (rs11031006 near FSHB, p = 1.54 × 10−9, and rs17293443 in SMAD3, p = 1.57 × 10−8) and replicated (p = 3 × 10−3 and p = 1.44 × 10−4, respectively). Based on ∼90,000 births in Iceland, the risk of a mother delivering twins increased by 18% for each copy of allele rs11031006-G and 9% for rs17293443-C. A higher polygenic risk score (PRS) for DZ twinning, calculated based on the results of the DZ twinning GWAS, was significantly associated with DZ twinning in Iceland (p = 0.001). A higher PRS was also associated with having children (p = 0.01), greater lifetime parity (p = 0.03), and earlier age at first child (p = 0.02). Allele rs11031006-G was associated with higher serum FSH levels, earlier age at menarche, earlier age at first child, higher lifetime parity, lower PCOS risk, and earlier age at menopause. Conversely, rs17293443-C was associated with later age at last child. We identified robust genetic risk variants for DZ twinning: one near FSHB and a second within SMAD3, the product of which plays an important role in gonadal responsiveness to FSH. These loci contribute to crucial aspects of reproductive capacity and health. 11 p.

DNA methylation, which is modulated by both genetic factors and environmental exposures, may offer a unique opportunity to discover novel biomarkers of disease-related brain phenotypes, even when measured in other tissues than brain, such as blood. A few studies of small sample sizes have revealed associations between blood DNA methylation and neuropsychopathology, however, large-scale epigenome-wide association studies (EWAS) are needed to investigate the utility of DNA methylation profiling as a peripheral marker for the brain. Here, in an analysis of eleven international cohorts, totalling 3,337 individuals, we report epigenome-wide meta-analyses of blood DNA methylation with volumes of the hippocampus, thalamus and nucleus accumbens (NAcc) –three subcortical regions selected for their associations with disease and heritability and volumetric variability. Analyses of individual CpGs revealed genome-wide significant associations with hippocampal volume at two loci. No significant associations were found for analyses of thalamus and nucleus accumbens volumes. CpG sites associated with hippocampus volume were significantly enriched within cancer-related genes and within regulatory elements containing the transcriptionally repressive histone H3K27 tri-methylation mark that is vital for stem cell fate specification. Cluster-based analyses revealed additional differentially methylated regions (DMRs) associated with hippocampal volume. DNA methylation at these loci affected expression of proximal genes involved in in learning and memory, stem cell maintenance and differentiation, fatty acid metabolism and type-2 diabetes. These DNA methylation marks, their interaction with genetic variants and their impact on gene expression offer new insights into the relationship between epigenetic variation and brain structure and may provide the basis for biomarker discovery in neurodegeneration and neuropsychiatric conditions.

Background Major challenges in understanding the functional consequences of genetic risk factors for human disease are which tissues and cell types are affected and the limited availability of suitable tissue. The aim of this study was to evaluate tissue-specific genotype-epigenetic characteristics in DNA samples from both endometrium and blood collected from women at different stages of the menstrual cycle and relate results to genetic risk factors for reproductive traits and diseases. Results We analysed DNA methylation (DNAm) data from endometrium and blood samples from 66 European women. Methylation profiles were compared between stages of the menstrual cycle, and changes in methylation overlaid with changes in transcription and genotypes. We observed large changes in methylation (27,262 DNAm probes) across the menstrual cycle in endometrium that were not observed in blood. Individual genotype data was tested for association with methylation at 443,016 and 443,101 DNAm probes in endometrium and blood respectively to identify methylation quantitative trait loci (mQTLs). A total of 4546 sentinel cis-mQTLs (P

Background DNA methylation levels change along with age, but few studies have examined the variation in the rate of such changes between individuals. Methods We performed a longitudinal analysis to quantify the variation in the rate of change of DNA methylation between individuals using whole blood DNA methylation array profiles collected at 2–4 time points (N = 2894) in 954 individuals (67–90 years). Results After stringent quality control, we identified 1507 DNA methylation CpG sites (rsCpGs) with statistically significant variation in the rate of change (random slope) of DNA methylation among individuals in a mixed linear model analysis. Genes in the vicinity of these rsCpGs were found to be enriched in Homeobox transcription factors and the Wnt signalling pathway, both of which are related to ageing processes. Furthermore, we investigated the SNP effect on the random slope. We found that 4 out of 1507 rsCpGs had one significant (P

BACKGROUNDEpigenetic age acceleration (an older methylation age compared to chronological age) correlates strongly with various age-related morbidities and mortality. Chronic systemic inflammation is thought to be a hallmark of ageing but the relationship between an increased epigenetic age and this likely key phenotype of ageing has not yet been extensively investigated.METHODSWe modelled the trajectories of the inflammatory biomarkers C-reactive protein (CRP; measured using both a high- and low-sensitivity assay), and interleukin-6 (IL-6) over the 8th decade in the Lothian Birth Cohort 1936. We additionally investigated the association between CRP and imputed leukocyte counts. Using linear mixed models we examined the cross-sectional and longitudinal association between the inflammatory biomarkers and two measures of epigenetic age acceleration, derived from the Horvath and Hannum epigenetic clocks.RESULTSLow-sensitivity CRP declined, high-sensitivity CRP did not change, and IL-6 increased over time. CRP levels inversely associated with total counts of CD8+T cells and CD4+T cells, and positively associated with senescent CD8+T cells, plasmablasts and granulocytes. Cross-sectionally, the Hannum, but not the Horvath, measure of age acceleration was positively associated with low-sensitivity CRP, high-sensitivity CRP, IL-6 and a restricted measure of CRP (≤10mg/L) likely reflecting levels relevant to chronic inflammation.CONCLUSIONSWe found a divergent relationship between inflammation and immune system parameters in older age. We additionally report the Hannum measure of epigenetic age acceleration associated with an elevated inflammatory profile cross-sectionally, but not longitudinally.

DNA methylation is associated with age. The deviation of age predicted from DNA methylation from actual age has been proposed as a biomarker for ageing. However, a better prediction of chronological age implies less opportunity for biological age. Here we used 13,661 samples (from blood and saliva) in the age range of 2 to 104 years from 14 cohorts measured on Illumina HumanMethylation450/EPIC arrays to perform prediction analyses. We show that increasing the sample size achieves a smaller prediction error and higher correlations in test datasets. We demonstrate that smaller prediction errors provide a limit to how much variation in biological ageing can be captured by methylation and provide evidence that age predictors from small samples are prone to confounding by cell composition. Our predictor shows a similar or better performance in non-blood tissues including saliva, endometrium, breast, liver, adipose and muscle, compared with Horvath’s across-tissue age predictor.

BackgroundGenome-wide DNA methylation (DNAm) profiling has allowed for the development of molecular predictors for a multitude of traits and diseases. Such predictors may be more accurate than the self-reported phenotypes, and could have clinical applications. Here, penalised regression models were used to develop DNAm predictors for body mass index (BMI), smoking status, alcohol consumption, and educational attainment in a cohort of 5,100 individuals. Using an independent test cohort comprising 906 individuals, the proportion of phenotypic variance explained in each trait was examined for DNAm-based and genetic predictors. Receiver operator characteristic curves were generated to investigate the predictive performance of DNAm-based predictors, using dichotomised phenotypes. The relationship between DNAm scores and all-cause mortality (n = 214 events) was assessed via Cox proportional-hazards models.ResultsThe DNAm-based predictors explained different proportions of the phenotypic variance for BMI (12%), smoking (60%), alcohol consumption (12%) and education (3%). The combined genetic and DNAm predictors explained 20% of the variance in BMI, 61% in smoking, 13% in alcohol consumption, and 6% in education. DNAm predictors for smoking, alcohol, and education but not BMI predicted mortality in univariate models. The predictors showed moderate discrimination of obesity (AUC=0.67) and alcohol consumption (AUC=0.75), and excellent discrimination of current smoking status (AUC=0.98). There was poorer discrimination of college-educated individuals (AUC=0.59).ConclusionsDNAm predictors correlate with lifestyle factors that are associated with health and mortality. They may supplement DNAm-based predictors of age to identify the lifestyle profiles of individuals and predict disease risk.List of abbreviationsDNAmDNA methylationBMIBody mass indexAUCArea under the curveCpGCytosine phosphate Guanine dinucleotideEWASEpigenome-wide association studyGS:SFHSGeneration Scotland: The Scottish family health studyLBC1936Lothian birth cohort 1936LASSOLeast absolute shrinkage and selector operatorHRHazard ratioCIConfidence intervalSTRADLStratifying resilience and depression longitudinally

BackgroundThe relationship between obesity and adverse health is well established, but little is known about the contribution of DNA methylation to obesity-related health outcomes. Additionally, it is of interest whether such contributions are independent of those attributed by the most widely used clinical measure of body mass – the Body Mass Index (BMI).MethodWe tested whether an epigenetic BMI score accounts for inter-individual variation in health-related, cognitive, psychosocial and lifestyle outcomes in the Lothian Birth Cohort 1936 (n=903). Weights for the epigenetic BMI score were derived using penalised regression on methylation data from unrelated Generation Scotland participants (n=2566).ResultsThe Epigenetic BMI score was associated with variables related to poor physical health (R2 ranges from 0.02-0.10), metabolic syndrome (R2 ranges from 0.01-0.09), lower crystallised intelligence (R2=0.01), lower health-related quality of life (R2=0.02), physical inactivity (R2=0.02), and social deprivation (R2=0.02). The epigenetic BMI score (per SD) was also associated with self-reported type 2 diabetes (OR 2.25, 95 % CI 1.74, 2.94), cardiovascular disease (OR 1.44, 95 % CI 1.23, 1.69) and high blood pressure (OR 1.21, 95% CI 1.13, 1.48; all at pConclusionsOur results show that regression models with epigenetic and phenotypic BMI scores as predictors account for a greater proportion of all outcome variables than either predictor alone, demonstrating independent and additive effects of epigenetic and phenotypic BMI scores.

Alzheimer’s disease (AD) is a public health priority for the 21st century. Risk reduction currently revolves around lifestyle changes with much research trying to elucidate the biological underpinnings. Using self-report of parental history of Alzheimer’s dementia for case ascertainment in a genome-wide association study of over 300,000 participants from UK Biobank (32,222 maternal cases, 16,613 paternal cases) and meta-analysing with published consortium data (n=74,046 with 25,580 cases across the discovery and replication analyses), six new AD-associated loci (P−8) are identified. Three contain genes relevant for AD and neurodegeneration: ADAM10, ADAMTS4, and ACE. Suggestive loci include drug targets such as VKORC1 (warfarin dose) and BZRAP1 (benzodiazepine receptor). We report evidence that association of SNPs and AD at the PVR gene is potentially mediated by both gene expression and DNA methylation in the prefrontal cortex. Our discovered loci may help to elucidate the biological mechanisms underlying AD and, given that many are existing drug targets for other diseases and disorders, warrant further exploration for potential precision medicine applications.

Age-related changes in DNA methylation have been observed in many cross-sectional studies, but longitudinal evidence is still very limited. Here, we aimed to characterize longitudinal age-related methylation patterns (Illumina HumanMethylation450 array) using 1011 blood samples collected from 385 old Swedish twins (mean age of 69 at baseline) up to five times over 20 years. We identified 1316 age-associated methylation sites (p−7) using a longitudinal epigenome-wide association study design. We measured how estimated cellular compositions changed with age and how much they confounded the age effect. We validated the results in two independent longitudinal cohorts, where 118 CpGs were replicated in PIVUS (p−5) and 594 were replicated in LBC (p−5). Functional annotation of age-associated CpGs showed enrichment in CCCTC-binding factor (CTCF) and other unannotated transcription factor binding sites. We further investigated genetic influences on methylation (methylation quantitative trait loci) and found no interaction between age and genetic effects in the 1316 age-associated CpGs. Moreover, in the same CpGs, methylation differences within twin pairs increased over time, where monozygotic twins had smaller intra-pair differences than dizygotic twins. We show that age-related methylation changes persist in a longitudinal perspective, and are fairly stable across cohorts. Moreover, the changes are under genetic influence, although this effect is independent of age. In addition, inter-individual methylation variations increase over time, especially in age-associated CpGs, indicating the increase of environmental contributions on DNA methylation with age.

Genome-wide association studies have identified hundreds of genetic variants associated with blood pressure (BP), but sequence variation accounts for a small fraction of the phenotypic variance. Epigenetic changes may alter the expression of genes involved in BP regulation and explain part of the missing heritability. We therefore conducted a two-stage meta-analysis of the cross-sectional associations of systolic and diastolic BP with blood-derived genome-wide DNA methylation measured on the Infinium HumanMethylation450 BeadChip in 17,010 individuals of European, African American, and Hispanic ancestry. Of 31 discovery-stage cytosine-phosphate-guanine (CpG) dinucleotides, 13 replicated after Bonferroni correction (discovery: N = 9,828, p < 1.0× 10-7; replication: N = 7,182, p 30%) and independent of known BP genetic variants, explaining an additional 1.4% and 2.0% of the interindividual variation in systolic and diastolic BP, respectively. Bidirectional Mendelian randomization among up to 4,513 individuals of European ancestry from 4 cohorts suggested that methylation at cg08035323 (TAF1B-YWHAQ) influences BP, while BP influences methylation at cg00533891 (ZMIZ1), cg00574958 (CPT1A), and cg02711608 (SLC1A5). Gene expression analyses further identified six genes (TSPAN2, SLC7A11, UNC93B1, CPT1A, PTMS, and LPCAT3) with evidence of triangular associations between methylation, gene expression, and BP. Additional integrative Mendelian randomization analyses of gene expression and DNA methylation suggested that the expression of TSPAN2 is a putative mediator of association between DNA methylation at cg23999170 and BP. These findings suggest that heritable DNA methylation plays a role in regulating BP independently of previously known genetic variants.

DNA methylation age is an accurate biomarker of chronological age and predicts lifespan, but its underlying molecular mechanisms are unknown. In this genome-wide association study of 9,907 individuals, we found gene variants mapping to five loci associated with intrinsic epigenetic age acceleration (IEAA) and gene variants in 3 loci associated extrinsic epigenetic age acceleration (EEAA). Mendelian randomization analysis suggested causal influences of menarche and menopause on IEAA and lipid levels on IEAA and EEAA. Variants associated with longer leukocyte telomere length (LTL) in the telomerase reverse transcriptase gene (TERT) locus at 5p15.33 confer higher IEAA (P-11). Causal modelling indicatesTERT-specific and independent effects on LTL and IEAA. Experimental hTERT expression in primary human fibroblasts engenders a linear increase in DNA methylation age with cell population doubling number. Together, these findings indicate a critical role for hTERT in regulating the DNA methylation clock, in addition to its established role of compensating for cell replication-dependent telomere shortening.

Estimation of the joint distribution of effect size and minor allele frequency (MAF) for genetic variants is important for understanding the genetic basis of complex trait variation and can be used to detect signature of natural selection. We develop a Bayesian mixed linear model that simultaneously estimates SNP-based heritability, polygenicity (i.e. the proportion of SNPs with nonzero effects) and the relationship between effect size and MAF for complex traits in conventionally unrelated individuals using genome-wide SNP data. We apply the method to 28 complex traits in the UK Biobank data (N = 126,752), and show that on average across 28 traits, 6% of SNPs have nonzero effects, which in total explain 22% of phenotypic variance. We detect significant (p < 0.05/28 =1.8×10−3) signatures of natural selection for 23 out of 28 traits including reproductive, cardiovascular, and anthropometric traits, as well as educational attainment. We further apply the method to 27,869 gene expression traits (N = 1,748), and identify 30 genes that show significant (p < 2.3×10−6) evidence of natural selection. All the significant estimates of the relationship between effect size and MAF in either complex traits or gene expression traits are consistent with a model of negative selection, as confirmed by forward simulation. We conclude that natural selection acts pervasively on human complex traits shaping genetic variation in the form of negative selection.

Background: Gene discovery has provided remarkable biological insights into amyotrophic lateral sclerosis (ALS). One challenge for clinical application of genetic testing is critical evaluation of the significance of reported variants. Methods: We use whole exome sequencing (WES) to develop a clinically relevant approach to identify a subset of ALS patients harboring likely pathogenic mutations. In parallel, we assess if DNA methylation can be used to screen for pathogenicity of novel variants since a methylation signature has been shown to associate with the pathogenic C9orf72 expansion, but has not been explored for other ALS mutations. Australian patients identified with ALS-relevant variants were cross-checked with population databases and case reports to critically assess whether they were “likely causal,” “uncertain significance,” or “unlikely causal.”. Results: Published ALS variants were identified in >10% of patients; however, in only 3% of patients (4/120) could these be confidently considered pathogenic (in SOD1 and TARDBP). We found no evidence for a differential DNA methylation signature in these mutation carriers. Conclusions: The use of WES in a typical ALS clinic demonstrates a critical approach to variant assessment with the capability to combine cohorts to enhance the largely unknown genetic basis of ALS Refereed/Peer-reviewed

The epigenome has been shown to be influenced by biological factors, such as disease status, and environmental factors, such as smoking, alcohol consumption, and body mass index. Although there is a widespread perception that environmental influences on the epigenome are pervasive and profound, there has been little evidence to date in humans with respect to environmental factors that are biologically distal. Here, we provide evidence on the associations between epigenetic modifications—in our case, CpG methylation—and educational attainment (EA), a biologically distal environmental factor that is arguably among of the most important life-shaping experiences for individuals. Specifically, we report the results of an epigenome-wide association study meta-analysis of EA based on data from 27 cohort studies with a total of 10,767 individuals. While we find that 9 CpG probes are significantly associated with EA, only two remain associated when we restrict the sample to never-smokers. These two are known to be strongly associated with maternal smoking during pregnancy, and thus their association with EA could be due to correlation between EA and maternal smoking. Moreover, their effect sizes on EA are far smaller than the known associations between CpG probes and biologically proximal environmental factors. Two analyses that combine the effects of many probes—polygenic methylation score and epigenetic-clock analyses—both suggest small associations with EA. If our findings regarding EA can be generalized to other biologically distal environmental factors, then they cast doubt on the hypothesis that such factors have large effects on the epigenome.

Background The link between DNA methylation, obesity, and adiposity-related diseases in the general population remains uncertain. Methods and Findings We conducted an association study of body mass index (BMI) and differential methylation for over 400,000 CpGs assayed by microarray in whole-blood-derived DNA from 3,743 participants in the Framingham Heart Study and the Lothian Birth Cohorts, with independent replication in three external cohorts of 4,055 participants. We examined variations in whole blood gene expression and conducted Mendelian randomization analyses to investigate the functional and clinical relevance of the findings. We identified novel and previously reported BMI-related differential methylation at 83 CpGs that replicated across cohorts; BMI-related differential methylation was associated with concurrent changes in the expression of genes in lipid metabolism pathways. Genetic instrumental variable analysis of alterations in methylation at one of the 83 replicated CpGs, cg11024682 (intronic to sterol regulatory element binding transcription factor 1 [SREBF1]), demonstrated links to BMI, adiposity-related traits, and coronary artery disease. Independent genetic instruments for expression of SREBF1 supported the findings linking methylation to adiposity and cardiometabolic disease. Methylation at a substantial proportion (16 of 83) of the identified loci was found to be secondary to differences in BMI. However, the cross-sectional nature of the data limits definitive causal determination. Conclusions We present robust associations of BMI with differential DNA methylation at numerous loci in blood cells. BMI-related DNA methylation and gene expression provide mechanistic insights into the relationship between DNA methylation, obesity, and adiposity-related diseases., Daniel Levy and colleagues examine mechanisms that could link body mass index to cardiovascular diseases via gene expression., Author Summary Why Was This Study Done? Genetic sequence variants explain only a modest proportion of the variation in body mass index (BMI) and cardiometabolic disease in the general population. There is limited understanding of the link of DNA methylation—a well-characterized epigenetic modification—with BMI and cardiometabolic disease in the general population. What Did the Researchers Do and Find? We conducted a cross-sectional analysis of the association of BMI with leukocyte DNA methylation at over 400,000 sites in the genome among 7,798 community-dwelling adults. We identified associations between BMI and methylation at 83 replicated sites (including 50 novel sites) and concurrent differences in expression in whole blood of genes overrepresented in lipid metabolism pathways. Using genetic sequence variants to model exposure to differential DNA methylation and tissue-specific gene expression, we found differential methylation and expression of SREBF1 to be implicated in BMI, adiposity-related traits, and coronary artery disease. Using genetic sequence variants to model exposure to differences in BMI, we found a substantial proportion of the differentially methylated sites (16 of 83) to be downstream of BMI. What Do These Findings Mean? Evidence is accumulating that epigenetic modifications, such as DNA methylation, are related to obesity-related diseases in the general population. We provide support for a role of genomic regulation of a lipid metabolism transcription factor, SREBF1, in adiposity and coronary artery disease. Mendelian randomization approaches can help prioritize relevant loci for future functional studies, but the cross-sectional observational nature of our study limits definitive causal inference.

Supplemental Digital Content is available in the text., Background— Genome-wide association studies have identified loci influencing circulating lipid concentrations in humans; further information on novel contributing genes, pathways, and biology may be gained through studies of epigenetic modifications. Methods and Results— To identify epigenetic changes associated with lipid concentrations, we assayed genome-wide DNA methylation at cytosine–guanine dinucleotides (CpGs) in whole blood from 2306 individuals from 2 population-based cohorts, with replication of findings in 2025 additional individuals. We identified 193 CpGs associated with lipid levels in the discovery stage (P

Epistasis has rarely been shown among natural polymorphisms in human traits; this research using advanced computation and gene expression data reveals many instances of epistasis between common single nucleotide polymorphisms in humans, with epistasis and the direction of its effect replicating in independent cohorts. Although frequently demonstrated in model organisms and domesticated species, very few examples of epistasis — where the effect of one polymorphism on a trait depends on other polymorphisms present in the genome — have been demonstrated in humans. Using advanced computation and a gene expression study design, these authors reveal hundreds of instances of epistasis between common single nucleotide polymorphisms (SNPs) in humans. They show that epistasis and direction of its effect replicate in independent cohorts. Epistatic networks of three or more SNPs are shown to influence the expression levels of many genes, whereby one cis-acting SNP is modulated by several trans-acting SNPs. Epistasis is the phenomenon whereby one polymorphism’s effect on a trait depends on other polymorphisms present in the genome. The extent to which epistasis influences complex traits1 and contributes to their variation2,3 is a fundamental question in evolution and human genetics. Although often demonstrated in artificial gene manipulation studies in model organisms4,5, and some examples have been reported in other species6, few examples exist for epistasis among natural polymorphisms in human traits7,8. Its absence from empirical findings may simply be due to low incidence in the genetic control of complex traits2,3, but an alternative view is that it has previously been too technically challenging to detect owing to statistical and computational issues9. Here we show, using advanced computation10 and a gene expression study design, that many instances of epistasis are found between common single nucleotide polymorphisms (SNPs). In a cohort of 846 individuals with 7,339 gene expression levels measured in peripheral blood, we found 501 significant pairwise interactions between common SNPs influencing the expression of 238 genes (P

Importance: Genome-wide association studies (GWASs) in European populations have identified more than 100 schizophrenia-associated loci. A schizophrenia GWAS in a unique Indian population offers novel findings. Objective: To discover and functionally evaluate genetic loci for schizophrenia in a GWAS of a unique Indian population. Design, Setting, and Participants: This GWAS included a sample of affected individuals, family members, and unrelated cases and controls. Three thousand ninety-two individuals were recruited and diagnostically ascertained via medical records, hospitals, clinics, and clinical networks in Chennai and surrounding regions. Affected participants fulfilled DSM-IV diagnostic criteria for schizophrenia. Unrelated control participants had no personal or family history of psychotic disorder. Recruitment, genotyping, and analysis occurred in consecutive phases beginning January 1, 2001. Recruitment was completed on February 28, 2018, and genotyping and analysis are ongoing. Main Outcomes and Measures: Associations of single-nucleotide polymorphisms and gene expression with schizophrenia. Results: The study population included 1321 participants with schizophrenia, 885 family controls, and 886 unrelated controls. Among participants with schizophrenia, mean (SD) age was 39.1 (11.4) years, and 52.7% were male. This sample demonstrated uniform ethnicity, a degree of inbreeding, and negligible rates of substance abuse. A novel genome-wide significant association was observed between schizophrenia and a chromosome 8q24.3 locus (rs10866912, allele A; odds ratio [OR], 1.27 [95% CI, 1.17-1.38]; P = 4.35 × 10) that attracted support in the schizophrenia Psychiatric Genomics Consortium 2 data (rs10866912, allele A; OR, 1.04 [95% CI, 1.02-1.06]; P = 7.56 × 10). This locus has undergone natural selection, with the risk allele A declining in frequency from India (approximately 72%) to Europe (approximately 43%). rs10866912 directly modifies the abundance of the nicotinate phosphoribosyltransferase gene (NAPRT1) transcript in brain cortex (normalized effect size, 0.79; 95% CI, 0.6-1.0; P = 5.8 × 10). NAPRT1 encodes a key enzyme for niacin metabolism. In Indian lymphoblastoid cell lines, (risk) allele A of rs10866912 was associated with NAPRT1 downregulation (AA: 0.74, n = 21; CC: 1.56, n = 17; P =.004). Preliminary zebrafish data further suggest that partial loss of function of NAPRT1 leads to abnormal brain development. Conclusions and Relevance: Bioinformatic analyses and cellular and zebrafish gene expression studies implicate NAPRT1 as a novel susceptibility gene. Given this gene's role in niacin metabolism and the evidence for niacin deficiency provoking schizophrenialike symptoms in neuropsychiatric diseases such as pellagra and Hartnup disease, these results suggest that the rs10866912 genotype and niacin status may have implications for schizophrenia susceptibility and treatment.

The mitochondrial and nuclear genomes coordinate and co-evolve in eukaryotes in order to adapt to environmental changes. Variation in the mitochondrial genome is capable of affecting expression of genes on the nuclear genome. Sex-specific mitochondrial genetic control of gene expression has been demonstrated in Drosophila melanogaster, where males were found to drive most of the total variation in gene expression. This has potential implications for male-related health and disease resulting from variation in mtDNA solely inherited from the mother. We used a family-based study comprised of 47,323 gene expression probes and 78 mitochondrial SNPs (mtSNPs) from n = 846 individuals to examine the extent of mitochondrial genetic control of gene expression in humans. This identified 15 significant probe-mtSNP associations (P

We propose a method (fastBAT) that performs a fast set-based association analysis for human complex traits using summary-level data from genome-wide association studies (GWAS) and linkage disequilibrium (LD) data from a reference sample with individual-level genotypes. We demonstrate using simulations and analyses of real datasets that fastBAT is more accurate and orders of magnitude faster than the prevailing methods. Using fastBAT, we analyze summary data from the latest meta-analyses of GWAS on 150,064–339,224 individuals for height, body mass index (BMI), and schizophrenia. We identify 6 novel gene loci for height, 2 for BMI, and 3 for schizophrenia at PfastBAT −8. The gain of power is due to multiple small independent association signals at these loci (e.g. the THRB and FOXP1 loci for schizophrenia). The method is general and can be applied to GWAS data for all complex traits and diseases in humans and to such data in other species.

Estimates of biological age based on DNA methylation patterns, often referred to as "epigenetic age", "DNAm age", have been shown to be robust biomarkers of age in humans. We previously demonstrated that independent of chronological age, epigenetic age assessed in blood predicted all-cause mortality in four human cohorts. Here, we expanded our original observation to 13 different cohorts for a total sample size of 13,089 individuals, including three racial/ethnic groups. In addition, we examined whether incorporating information on blood cell composition into the epigenetic age metrics improves their predictive power for mortality. All considered measures of epigenetic age acceleration were predictive of mortality (p ≤ 8.2 x 10-9), independent of chronological age, even after adjusting for additional risk factors (p -4), and within the racial/ethnic groups that we examined (non-Hispanic whites, Hispanics, African Americans). Epigenetic age estimates that incorporated information on blood cell composition led to the smallest p-values for time to death (p≤ 7.5 x 10-43). Overall, this study a) strengthens the evidence that epigenetic age predicts all-cause mortality above and beyond chronological age and traditional risk factors, and b) demonstrates that epigenetic age estimates that incorporate information on blood cell counts lead to highly significant associations with all-cause mortality.

Background Ankylosing spondylitis (AS) is an immune-mediated arthritis particularly targeting the spine and pelvis and is characterised by inflammation, osteoproliferation and frequently ankylosis. Current treatments that predominately target inflammatory pathways have disappointing efficacy in slowing disease progression. Thus, a better understanding of the causal association and pathological progression from inflammation to bone formation, particularly whether inflammation directly initiates osteoproliferation, is required. Methods The proteoglycan-induced spondylitis (PGISp) mouse model of AS was used to histopathologically map the progressive axial disease events, assess molecular changes during disease progression and define disease progression using unbiased clustering of semi-quantitative histology. PGISp mice were followed over a 24-week time course. Spinal disease was assessed using a novel semi-quantitative histological scoring system that independently evaluated the breadth of pathological features associated with PGISp axial disease, including inflammation, joint destruction and excessive tissue formation (osteoproliferation). Matrix components were identified using immunohistochemistry. Results Disease initiated with inflammation at the periphery of the intervertebral disc (IVD) adjacent to the longitudinal ligament, reminiscent of enthesitis, and was associated with upregulated tumor necrosis factor and metalloproteinases. After a lag phase, established inflammation was temporospatially associated with destruction of IVDs, cartilage and bone. At later time points, advanced disease was characterised by substantially reduced inflammation, excessive tissue formation and ectopic chondrocyte expansion. These distinct features differentiated affected mice into early, intermediate and advanced disease stages. Excessive tissue formation was observed in vertebral joints only if the IVD was destroyed as a consequence of the early inflammation. Ectopic excessive tissue was predominantly chondroidal with chondrocyte-like cells embedded within collagen type II- and X-rich matrix. This corresponded with upregulation of mRNA for cartilage markers Col2a1, sox9 and Comp. Osteophytes, though infrequent, were more prevalent in later disease. Conclusions The inflammation-driven IVD destruction was shown to be a prerequisite for axial disease progression to osteoproliferation in the PGISp mouse. Osteoproliferation led to vertebral body deformity and fusion but was never seen concurrent with persistent inflammation, suggesting a sequential process. The findings support that early intervention with anti-inflammatory therapies will be needed to limit destructive processes and consequently prevent progression of AS. Electronic supplementary material The online version of this article (doi:10.1186/s13075-015-0805-0) contains supplementary material, which is available to authorized users.

Human iris patterns are highly variable. The origins of this variation are of interest in the study of iris-related eye diseases and forensics, as well as from an embryological developmental perspective, with regard to their possible relationship to fundamental processes of neurodevelopment. We have performed genome-wide association scans on four iris characteristics (crypt frequency, furrow contractions, presence of peripupillary pigmented ring, and number of nevi) in three Australian samples of European descent. Both the discovery (n = 2121) and replication (n = 499 and 73) samples showed evidence for association between (1) crypt frequency and variants in the axonal guidance gene SEMA3A (p = 6.6 × 10(-11)), (2) furrow contractions and variants within the cytoskeleton gene TRAF3IP1 (p = 2.3 × 10(-12)), and (3) the pigmented ring and variants in the well-known pigmentation gene SLC24A4 (p = 7.6 × 10(-21)). These replicated findings individually accounted for around 1.5%-3% of the variance for these iris characteristics. Because both SEMA3A and TRAFIP1 are implicated in pathways that control neurogenesis, neural migration, and synaptogenesis, we also examined the evidence of enhancement among such genes, finding enrichment for crypts and furrows. These findings suggest that genes involved in normal neuronal pattern development may also influence tissue structures in the human iris.