A data-driven approach to preprocessing Illumina 450K methylation array data

As the most stable and experimentally accessible epigenetic mark, DNA methylation is of great interest to the research community. The landscape of DNA methylation across tissues, through development and in disease pathogenesis is not yet well characterized. Thus there is a need for rapid and cost effective methods for assessing genome-wide levels of DNA methylation. The Illumina Infinium HumanMethylation450 (450K) BeadChip is a very useful addition to the available methods for DNA methylation analysis but its complex design, incorporating two different assay methods, requires careful consideration. Accordingly, several normalization schemes have been published. We have taken advantage of known DNA methylation patterns associated with genomic imprinting and X-chromosome inactivation (XCI), in addition to the performance of SNP genotyping assays present on the array, to derive three independent metrics which we use to test alternative schemes of correction and normalization. These metrics also have potential utility as quality scores for datasets.|The standard index of DNA methylation at any specific CpG site is β = M/(M + U + 100) where M and U are methylated and unmethylated signal intensities, respectively. Betas (βs) calculated from raw signal intensities (the default GenomeStudio behavior) perform well, but using 11 methylomic datasets we demonstrate that quantile normalization methods produce marked improvement, even in highly consistent data, by all three metrics. The commonly used procedure of normalizing betas is inferior to the separate normalization of M and U, and it is also advantageous to normalize Type I and Type II assays separately. More elaborate manipulation of quantiles proves to be counterproductive.|Careful selection of preprocessing steps can minimize variance and thus improve statistical power, especially for the detection of the small absolute DNA methylation changes likely associated with complex disease phenotypes. For the convenience of the research community we have created a user-friendly R software package called wateRmelon, downloadable from bioConductor, compatible with the existing methylumi, minfi and IMA packages, that allows others to utilize the same normalization methods and data quality tests on 450K data

DNA methylation at the Igf2/H19 imprinting control region is associated with cerebellum mass in outbred mice

Background: Insulin-like growth factor 2 (Igf2) is a paternally expressed imprinted gene regulating fetal growth, playing an integral role in the development of many tissues including the brain. The parent-of-origin specific expression of Igf2 is largely controlled by allele-specific DNA methylation at CTCF-binding sites in the imprinting control region (ICR), located immediately upstream of the neighboring H19 gene. Previously we reported evidence of a negative correlation between DNA methylation in this region and cerebellum weight in humans. Results: We quantified cerebellar DNA methylation across all four CTCF binding sites spanning the murine Igf2/H19 ICR in an outbred population of Heterogeneous Stock (HS) mice (n = 48). We observe that DNA methylation at the second and third CTCF binding sites in the Igf2/H19 ICR shows a negative relationship with cerebellar mass, reflecting the association observed in human post-mortem cerebellum tissue. Conclusions: Given the important role of the cerebellum in motor control and cognition, and the link between structural cerebellar abnormalities and neuropsychiatric phenotypes, the identification of epigenetic factors associated with cerebellum growth and development may provide important insights about the etiology of psychiatric disorders

Schizophrenia-associated methylomic variation: molecular signatures of disease and polygenic risk burden across multiple brain regions

Genetic association studies provide evidence for a substantial polygenic component to schizophrenia, although the neurobiological mechanisms underlying the disorder remain largely undefined. Building on recent studies supporting a role for developmentally regulated epigenetic variation in the molecular aetiology of schizophrenia, this study aimed to identify epigenetic variation associated with both a diagnosis of schizophrenia and elevated polygenic risk burden for the disease across multiple brain regions. Genome-wide DNA methylation was quantified in 262 post-mortem brain samples, represent- ing tissue from four brain regions (prefrontal cortex, striatum, hippocampus and cerebellum) from 41 schizophrenia patients and 47 controls. We identified multiple disease-associated and polygenic risk score-associated differentially methylated posi- tions and regions, which are not enriched in genomic regions identified in genetic studies of schizophrenia and do not reflect direct genetic effects on DNA methylation. Our study represents the first analysis of epigenetic variation associated with schizophrenia across multiple brain regions and highlights the utility of polygenic risk scores for identifying molecular path- ways associated with aetiological variation in complex disease

Functional annotation of the human brain methylome identifies tissue-specific epigenetic variation across brain and blood

notes: PMCID: PMC3446315© 2012 Davies et al.; licensee BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.Dynamic changes to the epigenome play a critical role in establishing and maintaining cellular phenotype during differentiation, but little is known about the normal methylomic differences that occur between functionally distinct areas of the brain. We characterized intra- and inter-individual methylomic variation across whole blood and multiple regions of the brain from multiple donors

Methylome of human skeletal muscle after acute & chronic resistance exercise training, detraining & retraining

DNA methylation is an important epigenetic modification that can regulate gene expression following environmental encounters without changes to the genetic code. Using Infinium MethylationEPIC BeadChip Arrays (850,000 CpG sites) we analysed for the first time, DNA isolated from untrained human skeletal muscle biopsies (vastus lateralis) at baseline (rest) and immediately following an acute (single) bout of resistance exercise. In the same participants, we also analysed the methylome following a period of muscle growth (hypertrophy) evoked via chronic (repeated bouts-3 sessions/wk) resistance exercise (RE) (training) over 7-weeks, followed by complete exercise cessation for 7-weeks returning muscle back to baseline levels (detraining), and finally followed by a subsequent 7-week period of RE-induced hypertrophy (retraining). These valuable methylome data sets described in the present manuscript and deposited in an open-access repository can now be shared and re-used to enable the identification of epigenetically regulated genes/ networks that are modified after acute anabolic stimuli and hypertrophy, and further investigate the phenomenon of epigenetic memory in skeletal muscle

Mutation of the diamond-blackfan anemia gene Rps7 in mouse results in morphological and neuroanatomical phenotypes.

The ribosome is an evolutionarily conserved organelle essential for cellular function. Ribosome construction requires assembly of approximately 80 different ribosomal proteins (RPs) and four different species of rRNA. As RPs co-assemble into one multi-subunit complex, mutation of the genes that encode RPs might be expected to give rise to phenocopies, in which the same phenotype is associated with loss-of-function of each individual gene. However, a more complex picture is emerging in which, in addition to a group of shared phenotypes, diverse RP gene-specific phenotypes are observed. Here we report the first two mouse mutations (Rps7(Mtu) and Rps7(Zma)) of ribosomal protein S7 (Rps7), a gene that has been implicated in Diamond-Blackfan anemia. Rps7 disruption results in decreased body size, abnormal skeletal morphology, mid-ventral white spotting, and eye malformations. These phenotypes are reported in other murine RP mutants and, as demonstrated for some other RP mutations, are ameliorated by Trp53 deficiency. Interestingly, Rps7 mutants have additional overt malformations of the developing central nervous system and deficits in working memory, phenotypes that are not reported in murine or human RP gene mutants. Conversely, Rps7 mouse mutants show no anemia or hyperpigmentation, phenotypes associated with mutation of human RPS7 and other murine RPs, respectively. We provide two novel RP mouse models and expand the repertoire of potential phenotypes that should be examined in RP mutants to further explore the concept of RP gene-specific phenotypes

Methylation QTLs in the developing brain and their enrichment in schizophrenia risk loci

We characterized DNA methylation quantitative trait loci (mQTLs) in a large collection (n = 166) of human fetal brain samples spanning 56-166 d post-conception, identifying >16,000 fetal brain mQTLs. Fetal brain mQTLs were primarily cis-acting, enriched in regulatory chromatin domains and transcription factor binding sites, and showed substantial overlap with genetic variants that were also associated with gene expression in the brain. Using tissue from three distinct regions of the adult brain (prefrontal cortex, striatum and cerebellum), we found that most fetal brain mQTLs were developmentally stable, although a subset was characterized by fetal-specific effects. Fetal brain mQTLs were enriched amongst risk loci identified in a recent large-scale genome-wide association study (GWAS) of schizophrenia, a severe psychiatric disorder with a hypothesized neurodevelopmental component. Finally, we found that mQTLs can be used to refine GWAS loci through the identification of discrete sites of variable fetal brain methylation associated with schizophrenia risk variants

Variation in 5-hydroxymethylcytosine across human cortex and cerebellum

Background: The most widely utilized approaches for quantifying DNA methylation involve the treatment of genomic DNA with sodium bisulfite; however, this method cannot distinguish between 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC). Previous studies have shown that 5hmC is enriched in the brain, although little is known about its genomic distribution and how it differs between anatomical regions and individuals. In this study, we combine oxidative bisulfite (oxBS) treatment with the Illumina Infinium 450K BeadArray to quantify genome-wide patterns of 5hmC in two distinct anatomical regions of the brain from multiple individuals. Results: We identify 37,145 and 65,563 sites passing our threshold for detectable 5hmC in the prefrontal cortex and cerebellum respectively, with 23,445 loci common across both brain regions. Distinct patterns of 5hmC are identified in each brain region, with notable differences in the genomic location of the most hydroxymethylated loci between these brain regions. Tissue-specific patterns of 5hmC are subsequently confirmed in an independent set of prefrontal cortex and cerebellum samples. Conclusions: This study represents the first systematic analysis of 5hmC in the human brain, identifying tissue-specific hydroxymethylated positions and genomic regions characterized by inter-individual variation in DNA hydroxymethylation. This study demonstrates the utility of combining oxBS-treatment with the Illumina 450k methylation array to systematically quantify 5hmC across the genome and the potential utility of this approach for epigenomic studies of brain disorders

Epigenetic signatures of starting and stopping smoking

Acknowledgements: This work was supported by Alzheimer's Research UK Major Project Grant [ARUK–PG2017B–10]. Generation Scotland received core funding from the Chief Scientist Office of the Scottish Government Health Directorates [CZD/16/6] and the Scottish Funding Council [HR03006]. We are grateful to all the families who took part, the general practitioners and the Scottish School of Primary Care for their help in recruiting them, and the whole Generation Scotland team, which includes interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists, healthcare assistants, and nurses. Genotyping of the GS:SFHS samples was carried out by the Genetics Core Laboratory at the Wellcome Trust Clinical Research Facility, Edinburgh, Scotland and was funded by the Medical Research Council UK and the Wellcome Trust (Wellcome Trust Strategic Award “STratifying Resilience and Depression Longitudinally” (STRADL) [104036/Z/14/Z]. DNA methylation data collection was funded by the Wellcome Trust Strategic Award [10436/Z/14/Z]. The research was conducted in The University of Edinburgh Centre for Cognitive Ageing and Cognitive Epidemiology (CCACE), part of the cross–council Lifelong Health and Wellbeing Initiative [MR/K026992/1]; funding from the Biotechnology and Biological Sciences Research Council (BBSRC) and Medical Research Council (MRC) is gratefully acknowledged. CCACE supports Ian Deary, with some additional support from Dementias Platform UK [MR/L015382/1]. HCW is supported by a JMAS SIM fellowship from the Royal College of Physicians of Edinburgh. AMM and HCW have received support from the Sackler InstitutePeer reviewedPublisher PD