Histone Acetylome-wide Association Study of Autism Spectrum Disorder.

The association of histone modification changes with autism spectrum disorder (ASD) has not been systematically examined. We conducted a histone acetylome-wide association study (HAWAS) by performing H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq) on 257 postmortem samples from ASD and matched control brains. Despite etiological heterogeneity, ≥68% of syndromic and idiopathic ASD cases shared a common acetylome signature at >5,000 cis-regulatory elements in prefrontal and temporal cortex. Similarly, multiple genes associated with rare genetic mutations in ASD showed common "epimutations." Acetylome aberrations in ASD were not attributable to genetic differentiation at cis-SNPs and highlighted genes involved in synaptic transmission, ion transport, epilepsy, behavioral abnormality, chemokinesis, histone deacetylation, and immunity. By correlating histone acetylation with genotype, we discovered >2,000 histone acetylation quantitative trait loci (haQTLs) in human brain regions, including four candidate causal variants for psychiatric diseases. Due to the relative stability of histone modifications postmortem, we anticipate that the HAWAS approach will be applicable to multiple diseases.

Quantifying the proportion of different cell types in the human cortex using DNA methylation profiles.

BACKGROUND: Due to interindividual variation in the cellular composition of the human cortex, it is essential that covariates that capture these differences are included in epigenome-wide association studies using bulk tissue. As experimentally derived cell counts are often unavailable, computational solutions have been adopted to estimate the proportion of different cell types using DNA methylation data. Here, we validate and profile the use of an expanded reference DNA methylation dataset incorporating two neuronal and three glial cell subtypes for quantifying the cellular composition of the human cortex. RESULTS: We tested eight reference panels containing different combinations of neuronal- and glial cell types and characterised their performance in deconvoluting cell proportions from computationally reconstructed or empirically derived human cortex DNA methylation data. Our analyses demonstrate that while these novel brain deconvolution models produce accurate estimates of cellular proportions from profiles generated on postnatal human cortex samples, they are not appropriate for the use in prenatal cortex or cerebellum tissue samples. Applying our models to an extensive collection of empirical datasets, we show that glial cells are twice as abundant as neuronal cells in the human cortex and identify significant associations between increased Alzheimer's disease neuropathology and the proportion of specific cell types including a decrease in NeuNNeg/SOX10Neg nuclei and an increase of NeuNNeg/SOX10Pos nuclei. CONCLUSIONS: Our novel deconvolution models produce accurate estimates for cell proportions in the human cortex. These models are available as a resource to the community enabling the control of cellular heterogeneity in epigenetic studies of brain disorders performed on bulk cortex tissue.

Developmentally dynamic changes in DNA methylation in the human pancreas.

Development of the human pancreas requires the precise temporal control of gene expression via epigenetic mechanisms and the binding of key transcription factors. We quantified genome-wide patterns of DNA methylation in human fetal pancreatic samples from donors aged 6 to 21 post-conception weeks. We found dramatic changes in DNA methylation across pancreas development, with > 21% of sites characterized as developmental differentially methylated positions (dDMPs) including many annotated to genes associated with monogenic diabetes. An analysis of DNA methylation in postnatal pancreas tissue showed that the dramatic temporal changes in DNA methylation occurring in the developing pancreas are largely limited to the prenatal period. Significant differences in DNA methylation were observed between males and females at a number of autosomal sites, with a small proportion of sites showing sex-specific DNA methylation trajectories across pancreas development. Pancreas dDMPs were not distributed equally across the genome and were depleted in regulatory domains characterized by open chromatin and the binding of known pancreatic development transcription factors. Finally, we compared our pancreas dDMPs to previous findings from the human brain, identifying evidence for tissue-specific developmental changes in DNA methylation. This study represents the first systematic exploration of DNA methylation patterns during human fetal pancreas development and confirms the prenatal period as a time of major epigenomic plasticity.

Evaluation of nanopore sequencing for epigenetic epidemiology: a comparison with DNA methylation microarrays.

Most epigenetic epidemiology to date has utilized microarrays to identify positions in the genome where variation in DNA methylation is associated with environmental exposures or disease. However, these profile less than 3% of DNA methylation sites in the human genome, potentially missing affected loci and preventing the discovery of disrupted biological pathways. Third generation sequencing technologies, including Nanopore sequencing, have the potential to revolutionize the generation of epigenetic data, not only by providing genuine genome-wide coverage but profiling epigenetic modifications direct from native DNA. Here we assess the viability of using Nanopore sequencing for epidemiology by performing a comparison with DNA methylation quantified using the most comprehensive microarray available, the Illumina EPIC array. We implemented a CRISPR-Cas9 targeted sequencing approach in concert with Nanopore sequencing to profile DNA methylation in three genomic regions to attempt to rediscover genomic positions that existing technologies have shown are differentially methylated in tobacco smokers. Using Nanopore sequencing reads, DNA methylation was quantified at 1779 CpGs across three regions, providing a finer resolution of DNA methylation patterns compared to the EPIC array. The correlation of estimated levels of DNA methylation between platforms was high. Furthermore, we identified 12 CpGs where hypomethylation was significantly associated with smoking status, including 10 within the AHRR gene. In summary, Nanopore sequencing is a valid option for identifying genomic loci where large differences in DNAm are associated with a phenotype and has the potential to advance our understanding of the role differential methylation plays in the etiology of complex disease.

Exploring the mediation of DNA methylation across the epigenome between childhood adversity and First Episode of Psychosis-findings from the EU-GEI study.

ABTRACT: Studies conducted in psychotic disorders have shown that DNA-methylation (DNAm) is sensitive to the impact of Childhood Adversity (CA). However, whether it mediates the association between CA and psychosis is yet to be explored. Epigenome wide association studies (EWAS) using the Illumina Infinium-Methylation EPIC array in peripheral blood tissue from 366 First-episode of psychosis and 517 healthy controls was performed. Adversity scores were created for abuse, neglect and composite adversity with the Childhood Trauma Questionnaire (CTQ). Regressions examining (I) CTQ scores with psychosis; (II) with DNAm EWAS level and (III) between DNAm and caseness, adjusted for a variety of confounders were conducted. Divide-Aggregate Composite-null Test for the composite null-hypothesis of no mediation effect was conducted. Enrichment analyses were conducted with missMethyl package and the KEGG database. Our results show that CA was associated with psychosis (Composite: OR = 1.68; p = <0.001; abuse: OR = 2.16; p < 0.001; neglect: OR = 2.27; p = <0.001). None of the CpG sites significantly mediated the adversity-psychosis association after Bonferroni correction (p < 8.1 × 10-8). However, 28, 34 and 29 differentially methylated probes associated with 21, 27, 20 genes passed a less stringent discovery threshold (p < 5 × 10-5) for composite, abuse and neglect respectively, with a lack of overlap between abuse and neglect. These included genes previously associated to psychosis in EWAS studies, such as PANK1, SPEG TBKBP1, TSNARE1 or H2R. Downstream gene ontology analyses did not reveal any biological pathways that survived false discovery rate correction. Although at a non-significant level, DNAm changes in genes previously associated with schizophrenia in EWAS studies may mediate the CA-psychosis association. These results and associated involved processes such as mitochondrial or histaminergic disfunction, immunity or neural signalling requires replication in well powered samples. The lack of overlap between mediating genes associated with abuse and neglect suggests differential biological trajectories linking CA subtypes and psychosis.

Assessing the co-variability of DNA methylation across peripheral cells and tissues: Implications for the interpretation of findings in epigenetic epidemiology.

Most epigenome-wide association studies (EWAS) quantify DNA methylation (DNAm) in peripheral tissues such as whole blood to identify positions in the genome where variation is statistically associated with a trait or exposure. As whole blood comprises a mix of cell types, it is unclear whether trait-associated DNAm variation is specific to an individual cellular population. We collected three peripheral tissues (whole blood, buccal epithelial and nasal epithelial cells) from thirty individuals. Whole blood samples were subsequently processed using fluorescence-activated cell sorting (FACS) to purify five constituent cell-types (monocytes, granulocytes, CD4+ T cells, CD8+ T cells, and B cells). DNAm was profiled in all eight sample-types from each individual using the Illumina EPIC array. We identified significant differences in both the level and variability of DNAm between different sample types, and DNAm data-derived estimates of age and smoking were found to differ dramatically across sample types from the same individual. We found that for the majority of loci variation in DNAm in individual blood cell types was only weakly predictive of variance in DNAm measured in whole blood, although the proportion of variance explained was greater than that explained by either buccal or nasal epithelial samples. Covariation across sample types was much higher for DNAm sites influenced by genetic factors. Overall, we observe that DNAm variation in whole blood is additively influenced by a combination of the major blood cell types. For a subset of sites, however, variable DNAm detected in whole blood can be attributed to variation in a single blood cell type providing potential mechanistic insight about EWAS findings. Our results suggest that associations between whole blood DNAm and traits or exposures reflect differences in multiple cell types and our data will facilitate the interpretation of findings in epigenetic epidemiology.

A comparison of feature selection methodologies and learning algorithms in the development of a DNA methylation-based telomere length estimator.

BACKGROUND: The field of epigenomics holds great promise in understanding and treating disease with advances in machine learning (ML) and artificial intelligence being vitally important in this pursuit. Increasingly, research now utilises DNA methylation measures at cytosine-guanine dinucleotides (CpG) to detect disease and estimate biological traits such as aging. Given the challenge of high dimensionality of DNA methylation data, feature-selection techniques are commonly employed to reduce dimensionality and identify the most important subset of features. In this study, our aim was to test and compare a range of feature-selection methods and ML algorithms in the development of a novel DNA methylation-based telomere length (TL) estimator. We utilised both nested cross-validation and two independent test sets for the comparisons. RESULTS: We found that principal component analysis in advance of elastic net regression led to the overall best performing estimator when evaluated using a nested cross-validation analysis and two independent test cohorts. This approach achieved a correlation between estimated and actual TL of 0.295 (83.4% CI [0.201, 0.384]) on the EXTEND test data set. Contrastingly, the baseline model of elastic net regression with no prior feature reduction stage performed less well in general-suggesting a prior feature-selection stage may have important utility. A previously developed TL estimator, DNAmTL, achieved a correlation of 0.216 (83.4% CI [0.118, 0.310]) on the EXTEND data. Additionally, we observed that different DNA methylation-based TL estimators, which have few common CpGs, are associated with many of the same biological entities. CONCLUSIONS: The variance in performance across tested approaches shows that estimators are sensitive to data set heterogeneity and the development of an optimal DNA methylation-based estimator should benefit from the robust methodological approach used in this study. Moreover, our methodology which utilises a range of feature-selection approaches and ML algorithms could be applied to other biological markers and disease phenotypes, to examine their relationship with DNA methylation and predictive value.

Epigenetic age acceleration is associated with oligodendrocyte proportions in MSA and control brain tissue.

AIMS: Epigenetic clocks are widely applied as surrogates for biological age in different tissues and/or diseases, including several neurodegenerative diseases. Despite white matter (WM) changes often being observed in neurodegenerative diseases, no study has investigated epigenetic ageing in white matter. METHODS: We analysed the performances of two DNA methylation-based clocks, DNAmClockMulti and DNAmClockCortical , in post-mortem WM tissue from multiple subcortical regions and the cerebellum, and in oligodendrocyte-enriched nuclei. We also examined epigenetic ageing in control and multiple system atrophy (MSA) (WM and mixed WM and grey matter), as MSA is a neurodegenerative disease comprising pronounced WM changes and α-synuclein aggregates in oligodendrocytes. RESULTS: Estimated DNA methylation (DNAm) ages showed strong correlations with chronological ages, even in WM (e.g., DNAmClockCortical , r = [0.80-0.97], p < 0.05). However, performances and DNAm age estimates differed between clocks and brain regions. DNAmClockMulti significantly underestimated ages in all cohorts except in the MSA prefrontal cortex mixed tissue, whereas DNAmClockCortical tended towards age overestimations. Pronounced age overestimations in the oligodendrocyte-enriched cohorts (e.g., oligodendrocyte-enriched nuclei, p = 6.1 × 10-5 ) suggested that this cell type ages faster. Indeed, significant positive correlations were observed between estimated oligodendrocyte proportions and DNAm age acceleration estimated by DNAmClockCortical (r > 0.31, p < 0.05), and similar trends were obtained with DNAmClockMulti . Although increased age acceleration was observed in MSA compared with controls, no significant differences were detected upon adjustment for possible confounders (e.g., cell-type proportions). CONCLUSIONS: Our findings show that oligodendrocyte proportions positively influence epigenetic age acceleration across brain regions and highlight the need to further investigate this in ageing and neurodegeneration.

Brain DNA methylomic analysis of frontotemporal lobar degeneration reveals OTUD4 in shared dysregulated signatures across pathological subtypes.

Frontotemporal lobar degeneration (FTLD) is an umbrella term describing the neuropathology of a clinically, genetically and pathologically heterogeneous group of diseases, including frontotemporal dementia (FTD) and progressive supranuclear palsy (PSP). Among the major FTLD pathological subgroups, FTLD with TDP-43 positive inclusions (FTLD-TDP) and FTLD with tau-positive inclusions (FTLD-tau) are the most common, representing about 90% of the cases. Although alterations in DNA methylation have been consistently associated with neurodegenerative diseases, including Alzheimer's disease and Parkinson's disease, little is known for FTLD and its heterogeneous subgroups and subtypes. The main goal of this study was to investigate DNA methylation variation in FTLD-TDP and FTLD-tau. We used frontal cortex genome-wide DNA methylation profiles from three FTLD cohorts (142 FTLD cases and 92 controls), generated using the Illumina 450K or EPIC microarrays. We performed epigenome-wide association studies (EWAS) for each cohort followed by meta-analysis to identify shared differentially methylated loci across FTLD subgroups/subtypes. In addition, we used weighted gene correlation network analysis to identify co-methylation signatures associated with FTLD and other disease-related traits. Wherever possible, we also incorporated relevant gene/protein expression data. After accounting for a conservative Bonferroni multiple testing correction, the EWAS meta-analysis revealed two differentially methylated loci in FTLD, one annotated to OTUD4 (5'UTR-shore) and the other to NFATC1 (gene body-island). Of these loci, OTUD4 showed consistent upregulation of mRNA and protein expression in FTLD. In addition, in the three independent co-methylation networks, OTUD4-containing modules were enriched for EWAS meta-analysis top loci and were strongly associated with the FTLD status. These co-methylation modules were enriched for genes implicated in the ubiquitin system, RNA/stress granule formation and glutamatergic synaptic signalling. Altogether, our findings identified novel FTLD-associated loci, and support a role for DNA methylation as a mechanism involved in the dysregulation of biological processes relevant to FTLD, highlighting novel potential avenues for therapeutic development.

Uncovering the genetic architecture of broad antisocial behavior through a genome-wide association study meta-analysis.

Despite the substantial heritability of antisocial behavior (ASB), specific genetic variants robustly associated with the trait have not been identified. The present study by the Broad Antisocial Behavior Consortium (BroadABC) meta-analyzed data from 28 discovery samples (N = 85,359) and five independent replication samples (N = 8058) with genotypic data and broad measures of ASB. We identified the first significant genetic associations with broad ASB, involving common intronic variants in the forkhead box protein P2 (FOXP2) gene (lead SNP rs12536335, p = 6.32 × 10-10). Furthermore, we observed intronic variation in Foxp2 and one of its targets (Cntnap2) distinguishing a mouse model of pathological aggression (BALB/cJ strain) from controls (BALB/cByJ strain). Polygenic risk score (PRS) analyses in independent samples revealed that the genetic risk for ASB was associated with several antisocial outcomes across the lifespan, including diagnosis of conduct disorder, official criminal convictions, and trajectories of antisocial development. We found substantial genetic correlations of ASB with mental health (depression rg = 0.63, insomnia rg = 0.47), physical health (overweight rg = 0.19, waist-to-hip ratio rg = 0.32), smoking (rg = 0.54), cognitive ability (intelligence rg = -0.40), educational attainment (years of schooling rg = -0.46) and reproductive traits (age at first birth rg = -0.58, father's age at death rg = -0.54). Our findings provide a starting point toward identifying critical biosocial risk mechanisms for the development of ASB.

Genome-wide DNA methylation analysis of aggressive behaviour: a longitudinal population-based study.

BACKGROUND: Human aggression is influenced by an interplay between genetic predisposition and experience across the life span. This interaction is thought to occur through epigenetic mechanisms, inducing differential gene expression, thereby moderating neuronal cell and circuit function, and thus shaping aggressive behaviour. METHODS: Genome-wide DNA methylation (DNAm) levels were measured in peripheral blood obtained from 95 individuals participating in the Estonian Children Personality Behaviours and Health Study (ECPBHS) at 15 and 25 years of age. We examined the association between aggressive behaviour, as measured by Life History of Aggression (LHA) total score and DNAm levels both assessed at age 25. We further examined the pleiotropic effect of genetic variants regulating LHA-associated differentially methylated positions (DMPs) and multiple traits related to aggressive behaviours. Lastly, we tested whether the DNA methylomic loci identified in association with LHA at age 25 were also present at age 15. RESULTS: We found one differentially methylated position (DMP) (cg17815886; p = 1.12 × 10-8 ) and five differentially methylated regions (DMRs) associated with LHA after multiple testing adjustments. The DMP annotated to the PDLIM5 gene, and DMRs resided in the vicinity of four protein-encoding genes (TRIM10, GTF2H4, SLC45A4, B3GALT4) and a long intergenic non-coding RNA (LINC02068). We observed evidence for the colocalization of genetic variants associated with top DMPs and general cognitive function, educational attainment and cholesterol levels. Notably, a subset of the DMPs associated with LHA at age 25 also displayed altered DNAm patterns at age 15 with high accuracy in predicting aggression. CONCLUSIONS: Our findings highlight the potential role of DNAm in the development of aggressive behaviours. We observed pleiotropic genetic variants associated with identified DMPs, and various traits previously established to be relevant in shaping aggression in humans. The concordance of DNAm signatures in adolescents and young adults may have predictive value for inappropriate and maladaptive aggression later in life.

A comparison of blood and brain-derived ageing and inflammation-related DNA methylation signatures and their association with microglial burdens.

Inflammation and ageing-related DNA methylation patterns in the blood have been linked to a variety of morbidities, including cognitive decline and neurodegenerative disease. However, it is unclear how these blood-based patterns relate to patterns within the brain and how each associates with central cellular profiles. In this study, we profiled DNA methylation in both the blood and in five post mortem brain regions (BA17, BA20/21, BA24, BA46 and hippocampus) in 14 individuals from the Lothian Birth Cohort 1936. Microglial burdens were additionally quantified in the same brain regions. DNA methylation signatures of five epigenetic ageing biomarkers ('epigenetic clocks'), and two inflammatory biomarkers (methylation proxies for C-reactive protein and interleukin-6) were compared across tissues and regions. Divergent associations between the inflammation and ageing signatures in the blood and brain were identified, depending on region assessed. Four out of the five assessed epigenetic age acceleration measures were found to be highest in the hippocampus (ß range = 0.83-1.14, p ≤ 0.02). The inflammation-related DNA methylation signatures showed no clear variation across brain regions. Reactive microglial burdens were found to be highest in the hippocampus (ß = 1.32, p = 5 × 10-4 ); however, the only association identified between the blood- and brain-based methylation signatures and microglia was a significant positive association with acceleration of one epigenetic clock (termed DNAm PhenoAge) averaged over all five brain regions (ß = 0.40, p = 0.002). This work highlights a potential vulnerability of the hippocampus to epigenetic ageing and provides preliminary evidence of a relationship between DNA methylation signatures in the brain and differences in microglial burdens.

Growth disrupting mutations in epigenetic regulatory molecules are associated with abnormalities of epigenetic aging.

Germline mutations in fundamental epigenetic regulatory molecules including DNA methyltransferase 3 alpha (DNMT3A) are commonly associated with growth disorders, whereas somatic mutations are often associated with malignancy. We profiled genome-wide DNA methylation patterns in DNMT3A c.2312G > A; p.(Arg771Gln) carriers in a large Amish sibship with Tatton-Brown-Rahman syndrome (TBRS), their mosaic father, and 15 TBRS patients with distinct pathogenic de novo DNMT3A variants. This defined widespread DNA hypomethylation at specific genomic sites enriched at locations annotated as genes involved in morphogenesis, development, differentiation, and malignancy predisposition pathways. TBRS patients also displayed highly accelerated DNA methylation aging. These findings were most marked in a carrier of the AML-associated driver mutation p.Arg882Cys. Our studies additionally defined phenotype-related accelerated and decelerated epigenetic aging in two histone methyltransferase disorders: NSD1 Sotos syndrome overgrowth disorder and KMT2D Kabuki syndrome growth impairment. Together, our findings provide fundamental new insights into aberrant epigenetic mechanisms, the role of epigenetic machinery maintenance, and determinants of biological aging in these growth disorders.

DNA methylation signatures of aggression and closely related constructs: A meta-analysis of epigenome-wide studies across the lifespan.

DNA methylation profiles of aggressive behavior may capture lifetime cumulative effects of genetic, stochastic, and environmental influences associated with aggression. Here, we report the first large meta-analysis of epigenome-wide association studies (EWAS) of aggressive behavior (N = 15,324 participants). In peripheral blood samples of 14,434 participants from 18 cohorts with mean ages ranging from 7 to 68 years, 13 methylation sites were significantly associated with aggression (alpha = 1.2 × 10-7; Bonferroni correction). In cord blood samples of 2425 children from five cohorts with aggression assessed at mean ages ranging from 4 to 7 years, 83% of these sites showed the same direction of association with childhood aggression (r = 0.74, p = 0.006) but no epigenome-wide significant sites were found. Top-sites (48 at a false discovery rate of 5% in the peripheral blood meta-analysis or in a combined meta-analysis of peripheral blood and cord blood) have been associated with chemical exposures, smoking, cognition, metabolic traits, and genetic variation (mQTLs). Three genes whose expression levels were associated with top-sites were previously linked to schizophrenia and general risk tolerance. At six CpGs, DNA methylation variation in blood mirrors variation in the brain. On average 44% (range = 3-82%) of the aggression-methylation association was explained by current and former smoking and BMI. These findings point at loci that are sensitive to chemical exposures with potential implications for neuronal functions. We hope these results to be a starting point for studies leading to applications as peripheral biomarkers and to reveal causal relationships with aggression and related traits.

Functional characterization of the schizophrenia associated gene AS3MT identifies a role in neuronal development.

Genome-wide association studies (GWAS) have identified multiple genomic regions associated with schizophrenia, although many variants reside in noncoding regions characterized by high linkage disequilibrium (LD) making the elucidation of molecular mechanisms challenging. A genomic region on chromosome 10q24 has been consistently associated with schizophrenia with risk attributed to the AS3MT gene. Although AS3MT is hypothesized to play a role in neuronal development and differentiation, work to fully understand the function of this gene has been limited. In this study we explored the function of AS3MT using a neuronal cell line (SH-SY5Y). We confirm previous findings of isoform specific expression of AS3MT during SH-SY5Y differentiation toward neuronal fates. Using CRISPR-Cas9 gene editing we generated AS3MT knockout SH-SY5Y cell lines and used RNA-seq to identify significant changes in gene expression in pathways associated with neuronal development, inflammation, extracellular matrix formation, and RNA processing, including dysregulation of other genes strongly implicated in schizophrenia. We did not observe any morphological changes in cell size and neurite length following neuronal differentiation and MAP2 immunocytochemistry. These results provide novel insights into the potential role of AS3MT in brain development and identify pathways through which genetic variation in this region may confer risk for schizophrenia.

Characterizing the properties of bisulfite sequencing data: maximizing power and sensitivity to identify between-group differences in DNA methylation.

BACKGROUND: The combination of sodium bisulfite treatment with highly-parallel sequencing is a common method for quantifying DNA methylation across the genome. The power to detect between-group differences in DNA methylation using bisulfite-sequencing approaches is influenced by both experimental (e.g. read depth, missing data and sample size) and biological (e.g. mean level of DNA methylation and difference between groups) parameters. There is, however, no consensus about the optimal thresholds for filtering bisulfite sequencing data with implications for the reproducibility of findings in epigenetic epidemiology. RESULTS: We used a large reduced representation bisulfite sequencing (RRBS) dataset to assess the distribution of read depth across DNA methylation sites and the extent of missing data. To investigate how various study variables influence power to identify DNA methylation differences between groups, we developed a framework for simulating bisulfite sequencing data. As expected, sequencing read depth, group size, and the magnitude of DNA methylation difference between groups all impacted upon statistical power. The influence on power was not dependent on one specific parameter, but reflected the combination of study-specific variables. As a resource to the community, we have developed a tool, POWEREDBiSeq, which utilizes our simulation framework to predict study-specific power for the identification of DNAm differences between groups, taking into account user-defined read depth filtering parameters and the minimum sample size per group. CONCLUSIONS: Our data-driven approach highlights the importance of filtering bisulfite-sequencing data by minimum read depth and illustrates how the choice of threshold is influenced by the specific study design and the expected differences between groups being compared. The POWEREDBiSeq tool, which can be applied to different types of bisulfite sequencing data (e.g. RRBS, whole genome bisulfite sequencing (WGBS), targeted bisulfite sequencing and amplicon-based bisulfite sequencing), can help users identify the level of data filtering needed to optimize power and aims to improve the reproducibility of bisulfite sequencing studies.

DNA methylation-based sex classifier to predict sex and identify sex chromosome aneuploidy.

BACKGROUND: Sex is an important covariate of epigenome-wide association studies due to its strong influence on DNA methylation patterns across numerous genomic positions. Nevertheless, many samples on the Gene Expression Omnibus (GEO) frequently lack a sex annotation or are incorrectly labelled. Considering the influence that sex imposes on DNA methylation patterns, it is necessary to ensure that methods for filtering poor samples and checking of sex assignment are accurate and widely applicable. RESULTS: Here we presented a novel method to predict sex using only DNA methylation beta values, which can be readily applied to almost all DNA methylation datasets of different formats (raw IDATs or text files with only signal intensities) uploaded to GEO. We identified 4345 significantly (p<0.01) sex-associated CpG sites present on both 450K and EPIC arrays, and constructed a sex classifier based on the two first principal components of the DNA methylation data of sex-associated probes mapped on sex chromosomes. The proposed method is constructed using whole blood samples and exhibits good performance across a wide range of tissues. We further demonstrated that our method can be used to identify samples with sex chromosome aneuploidy, this function is validated by five Turner syndrome cases and one Klinefelter syndrome case. CONCLUSIONS: This proposed sex classifier not only can be used for sex predictions but also applied to identify samples with sex chromosome aneuploidy, and it is freely and easily accessible by calling the 'estimateSex' function from the newest wateRmelon Bioconductor package ( https://github.com/schalkwyk/wateRmelon ).

The association of epigenetic clocks in brain tissue with brain pathologies and common aging phenotypes.

Epigenetic clocks are calculated by combining DNA methylation states across select CpG sites to estimate biologic age, and have been noted as the most successful markers of biologic aging to date. Yet, limited research has considered epigenetic clocks calculated in brain tissue. We used DNA methylation states in dorsolateral prefrontal cortex specimens from 721 older participants of the Religious Orders Study and Rush Memory and Aging Project, to calculate DNA methylation age using four established epigenetic clocks: Hannum, Horvath, PhenoAge, GrimAge, and a new Cortical clock. The four established clocks were trained in blood samples (Hannum, PhenoAge, GrimAge) or using 51 human tissue and cell types (Horvath); the recent Cortical clock is the first trained in postmortem cortical tissue. Participants were recruited beginning in 1994 (Religious Orders Study) and 1997 (Memory and Aging Project), and followed annually with questionnaires and clinical evaluations; brain specimens were obtained for 80-90% of participants. Mean age at death was 88.0 (SD 6.7) years. We used linear regression, logistic regression, and linear mixed models, to examine relations of epigenetic clock ages to neuropathologic and clinical aging phenotypes, controlling for chronologic age, sex, education, and depressive symptomatology. Hannum, Horvath, PhenoAge and Cortical clock ages were related to pathologic diagnosis of Alzheimer's disease (AD), as well as to Aß load (a hallmark pathology of Alzheimer's disease). However, associations were substantially stronger for the Cortical than other clocks; for example, each standard deviation (SD) increase in Hannum, Horvath, and PhenoAge clock age was related to approximately 30% greater likelihood of pathologic AD (all p < 0.05), while each SD increase in Cortical age was related to 90% greater likelihood of pathologic AD (odds ratio = 1.91, 95% confidence interval 1.38, 2.62). Moreover, Cortical age was significantly related to other AD pathology (eg, mean tau tangle density, p = 0.003), and to odds of neocortical Lewy body pathology (for each SD increase in Cortical age, odds ratio = 2.00, 95% confidence 1.27, 3.17), although no clocks were related to cerebrovascular neuropathology. Cortical age was the only epigenetic clock significantly associated with the clinical phenotypes examined, from dementia to cognitive decline (5 specific cognitive systems, and a global cognitive measure averaging 17 tasks) to Parkinsonian signs. Overall, our findings provide evidence of the critical necessity for bespoke clocks of brain aging for advancing research to understand, and eventually prevent, neurodegenerative diseases of aging.

Genome-wide DNA methylation profiling identifies convergent molecular signatures associated with idiopathic and syndromic autism in post-mortem human brain tissue.

Autism spectrum disorder (ASD) encompasses a collection of complex neuropsychiatric disorders characterized by deficits in social functioning, communication and repetitive behaviour. Building on recent studies supporting a role for developmentally moderated regulatory genomic variation in the molecular aetiology of ASD, we quantified genome-wide patterns of DNA methylation in 223 post-mortem tissues samples isolated from three brain regions [prefrontal cortex, temporal cortex and cerebellum (CB)] dissected from 43 ASD patients and 38 non-psychiatric control donors. We identified widespread differences in DNA methylation associated with idiopathic ASD (iASD), with consistent signals in both cortical regions that were distinct to those observed in the CB. Individuals carrying a duplication on chromosome 15q (dup15q), representing a genetically defined subtype of ASD, were characterized by striking differences in DNA methylationacross a discrete domain spanning an imprinted gene cluster within the duplicated region. In addition to the dramatic cis-effects on DNA methylation observed in dup15q carriers, we identified convergent methylomic signatures associated with both iASD and dup15q, reflecting the findings from previous studies of gene expression and H3K27ac. Cortical co-methylation network analysis identified a number of co-methylated modules significantly associated with ASD that are enriched for genomic regions annotated to genes involved in the immune system, synaptic signalling and neuronal regulation. Our study represents the first systematic analysis of DNA methylation associated with ASD across multiple brain regions, providing novel evidence for convergent molecular signatures associated with both idiopathic and syndromic autism.

Leveraging DNA-Methylation Quantitative-Trait Loci to Characterize the Relationship between Methylomic Variation, Gene Expression, and Complex Traits.

Characterizing the complex relationship between genetic, epigenetic, and transcriptomic variation has the potential to increase understanding about the mechanisms underpinning health and disease phenotypes. We undertook a comprehensive analysis of common genetic variation on DNA methylation (DNAm) by using the Illumina EPIC array to profile samples from the UK Household Longitudinal study. We identified 12,689,548 significant DNA methylation quantitative trait loci (mQTL) associations (p < 6.52 × 10-14) occurring between 2,907,234 genetic variants and 93,268 DNAm sites, including a large number not identified by previous DNAm-profiling methods. We demonstrate the utility of these data for interpreting the functional consequences of common genetic variation associated with > 60 human traits by using summary-data-based Mendelian randomization (SMR) to identify 1,662 pleiotropic associations between 36 complex traits and 1,246 DNAm sites. We also use SMR to characterize the relationship between DNAm and gene expression and thereby identify 6,798 pleiotropic associations between 5,420 DNAm sites and the transcription of 1,702 genes. Our mQTL database and SMR results are available via a searchable online database as a resource to the research community.