Search

Systemic low-grade inflammation is a feature of chronic disease. C-reactive protein (CRP) is a common biomarker of inflammation and used as an indicator of disease risk; however, the role of inflammation in disease is not completely understood. Methylation is an epigenetic modification in the DNA which plays a pivotal role in gene expression. In this study we evaluated differential DNA methylation patterns associated with blood CRP level to elucidate biological pathways and genetic regulatory mechanisms to improve the understanding of chronic inflammation. The racially and ethnically diverse participants in this study were included as 50% White, 41% Black or African American, 7% Hispanic or Latino/a, and 2% Native Hawaiian, Asian American, American Indian, or Alaska Native (total n = 13,433) individuals. We replicated 113 CpG sites from 87 unique loci, of which five were novel (CADM3, NALCN, NLRC5, ZNF792, and cg03282312), across a discovery set of 1,150 CpG sites associated with CRP level (p

Worldwide trends to delay childbearing have increased parental ages at birth. Older parental age may harm offspring health, but mechanisms remain unclear. Alterations in offspring DNA methylation (DNAm) patterns could play a role as aging has been associated with methylation changes in gametes of older individuals. We meta-analyzed epigenome-wide associations of parental age with offspring blood DNAm of over 9500 newborns and 2000 children (5-10 years old) from the Pregnancy and Childhood Epigenetics consortium. In newborns, we identified 33 CpG sites in 13 loci with DNAm associated with maternal age (PFDR

Whole-genome bisulfite sequencing (BS-Seq) measures cytosine methylation changes at single-base resolution and can be used to profile cell-free DNA (cfDNA). In plasma, ultrashort single-stranded cfDNA (uscfDNA, ∼50 nt) has been identified together with 167 bp double-stranded mononucleosomal cell-free DNA (mncfDNA). However, the methylation profile of uscfDNA has not been described. Conventional BS-Seq workflows may not be helpful because bisulfite conversion degrades larger DNA into smaller fragments, leading to erroneous categorization as uscfDNA. We describe the '5mCAdpBS-Seq' workflow in which pre-methylated 5mC (5-methylcytosine) single-stranded adapters are ligated to heat-denatured cfDNA before bisulfite conversion. This method retains only DNA fragments that are unaltered by bisulfite treatment, resulting in less biased uscfDNA methylation analysis. Using 5mCAdpBS-Seq, uscfDNA had lower levels of DNA methylation (∼15%) compared to mncfDNA and was enriched in promoters and CpG islands. Hypomethylated uscfDNA fragments were enriched in upstream transcription start sites (TSSs), and the intensity of enrichment was correlated with expressed genes of hemopoietic cells. Using tissue-of-origin deconvolution, we inferred that uscfDNA is derived primarily from eosinophils, neutrophils, and monocytes. As proof-of-principle, we show that characteristics of the methylation profile of uscfDNA can distinguish non-small cell lung carcinoma from non-cancer samples. The 5mCAdpBS-Seq workflow is recommended for any cfDNA methylation-based investigations.

BackgroundSystemic lupus erythematosus (SLE) has numerous symptoms across organs and an unpredictable flare-remittance pattern. This has made it challenging to understand drivers of long-term SLE outcomes. Our objective was to identify whether changes in DNA methylation over time, in an actively flaring SLE cohort, were associated with remission and whether these changes meaningfully subtype SLE patients.MethodsFifty-nine multi-ethnic SLE patients had clinical visits and DNA methylation profiles at a flare and approximately 3 months later. Methylation was measured using the Illumina EPIC array. We identified sites where methylation change between visits was associated with remission at the follow-up visit using limma package and a time x remission interaction term. Models adjusted for batch, age at diagnosis, time between visits, age at flare, sex, medications, and cell-type proportions. Separately, a paired T-test identified Bonferroni significant methylation sites with ≥ 3% change between visits (n = 546). Methylation changes at these sites were used for unsupervised consensus hierarchical clustering. Associations between clusters and patient features were assessed.ResultsNineteen patients fully remitted at the follow-up visit. For 1,953 CpG sites, methylation changed differently for remitters vs. non-remitters (Bonferroni p

Maternal use of e-cigarette (e-cig) aerosols poses significant risks to fetal brain development, potentially increasing susceptibility to neurodevelopmental disorders in later life. However, the underlying mechanisms remain incompletely understood. This study aimed to understand the effects of fetal e-cig exposure on DNA methylome and transcriptomic changes in the neonatal brain. Pregnant rats were exposed to e-cig aerosols, and neonatal brains (5 males and 5 females/group) from both control and e-cig-exposed groups were used for experimental analysis. Results indicated that prenatal e-cig exposure altered site-specific DNA methylation patterns at both CpG and CH (non-CpG) sites, predominantly in intergenic and intronic regions, with sex dimorphism in methylation and gene expression changes. Gene ontology analysis revealed that e-cig exposure not only affected neuron projection development and axonogenesis but also altered pathways related to neurodegeneration and long-term depression. These findings provide novel insights into the dynamic changes of CpG and CH methylation induced by e-cig exposure, underscoring the susceptibility of the developing brain to maternal e-cig exposure and its potential implications for developmental disorders and neurodegenerative diseases later in life.

BackgroundEpigenome-wide association studies have identified multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases. This study aimed to test the hypothesis that an alcohol consumption epigenetic risk score (ERS) is associated with blood pressure (BP) traits.ResultsWe implemented an ERS based on a previously reported epigenetic signature of 144 alcohol-associated CpGs in meta-analysis of participants of European ancestry. We found a one-unit increment of ERS was associated with eleven drinks of alcohol consumed per day, on average, across several cohorts (p

Chronic kidney disease (CKD) impacts about 1 in 7 adults in the United States, but African Americans (AAs) carry a disproportionately higher burden of disease. Epigenetic modifications, such as DNA methylation at cytosine-phosphate-guanine (CpG) sites, have been linked to kidney function and may have clinical utility in predicting the risk of CKD. Given the dynamic relationship between the epigenome, environment, and disease, AAs may be especially sensitive to environment-driven methylation alterations. Moreover, risk models incorporating CpG methylation have been shown to predict disease across multiple racial groups. In this study, we developed a methylation risk score (MRS) for CKD in cohorts of AAs. We selected nine CpG sites that were previously reported to be associated with estimated glomerular filtration rate (eGFR) in epigenome-wide association studies to construct a MRS in the Hypertension Genetic Epidemiology Network (HyperGEN). In logistic mixed models, the MRS was significantly associated with prevalent CKD and was robust to multiple sensitivity analyses, including CKD risk factors. There was modest replication in validation cohorts. In summary, we demonstrated that an eGFR-based CpG score is an independent predictor of prevalent CKD, suggesting that MRS should be further investigated for clinical utility in evaluating CKD risk and progression.

Age-associated changes in DNA methylation have been characterized across various animals, but not yet in amphibians, which are of particular interest because they include widely studied model organisms. In this study, we present clear evidence that the aquatic vertebrate species Xenopus tropicalis displays patterns of age-associated changes in DNA methylation. We have generated whole-genome bisulfite sequencing (WGBS) profiles from skin samples of nine frogs representing young, mature, and old adults and characterized the gene- and chromosome-scale DNA methylation changes with age. Many of the methylation features and changes we observe are consistent with what is known in mammalian species, suggesting that the mechanism of age-related changes is conserved. Moreover, we selected a few thousand age-associated CpG sites to build an assay based on targeted DNA methylation analysis (TBSseq) to expand our findings in future studies involving larger cohorts of individuals. Preliminary results of a pilot TBSeq experiment recapitulate the findings obtained with WGBS setting the basis for the development of an epigenetic clock assay. The results of this study will allow us to leverage the unique resources available for Xenopus to study how DNA methylation relates to other hallmarks of ageing.

Single-nucleotide variants (SNVs) in segmental duplications (SDs) have not been systematically assessed because of the limitations of mapping short-read sequencing data1,2. Here we constructed 1:1 unambiguous alignments spanning high-identity SDs across 102 human haplotypes and compared the pattern of SNVs between unique and duplicated regions3,4. We find that human SNVs are elevated 60% in SDs compared to unique regions and estimate that at least 23% of this increase is due to interlocus gene conversion (IGC) with up to 4.3 megabase pairs of SD sequence converted on average per human haplotype. We develop a genome-wide map of IGC donors and acceptors, including 498 acceptor and 454 donor hotspots affecting the exons of about 800 protein-coding genes. These include 171 genes that have 'relocated' on average 1.61 megabase pairs in a subset of human haplotypes. Using a coalescent framework, we show that SD regions are slightly evolutionarily older when compared to unique sequences, probably owing to IGC. SNVs in SDs, however, show a distinct mutational spectrum: a 27.1% increase in transversions that convert cytosine to guanine or the reverse across all triplet contexts and a 7.6% reduction in the frequency of CpG-associated mutations when compared to unique DNA. We reason that these distinct mutational properties help to maintain an overall higher GC content of SD DNA compared to that of unique DNA, probably driven by GC-biased conversion between paralogous sequences5,6.

Despite the prognostic value of arterial stiffness (AS) and pulsatile hemodynamics (PH) for cardiovascular morbidity and mortality, epigenetic modifications that contribute to AS/PH remain unknown. To gain a better understanding of the link between epigenetics (DNA methylation) and AS/PH, we examined the relationship of eight measures of AS/PH with CpG sites and co-methylated regions using multi-ancestry participants from Trans-Omics for Precision Medicine (TOPMed) Multi-Ethnic Study of Atherosclerosis (MESA) with sample sizes ranging from 438 to 874. Epigenome-wide association analysis identified one genome-wide significant CpG (cg20711926-CYP1B1) associated with aortic augmentation index (AIx). Follow-up analyses, including gene set enrichment analysis, expression quantitative trait methylation analysis, and functional enrichment analysis on differentially methylated positions and regions, further prioritized three CpGs and their annotated genes (cg23800023-ETS1, cg08426368-TGFB3, and cg17350632-HLA-DPB1) for AIx. Among these, ETS1 and TGFB3 have been previously prioritized as candidate genes. Furthermore, both ETS1 and HLA-DPB1 have significant tissue correlations between Whole Blood and Aorta in GTEx, which suggests ETS1 and HLA-DPB1 could be potential biomarkers in understanding pathophysiology of AS/PH. Overall, our findings support the possible role of epigenetic regulation via DNA methylation of specific genes associated with AIx as well as identifying potential targets for regulation of AS/PH.

BackgroundCongenital heart defects (CHDs) affect approximately half of individuals with Down syndrome (DS), but the molecular reasons for incomplete penetrance are unknown. Previous studies have largely focused on identifying genetic risk factors associated with CHDs in individuals with DS, but comprehensive studies of the contribution of epigenetic marks are lacking. We aimed to identify and characterize DNA methylation differences from newborn dried blood spots (NDBS) of DS individuals with major CHDs compared to DS individuals without CHDs.MethodsWe used the Illumina EPIC array and whole-genome bisulfite sequencing (WGBS) to quantitate DNA methylation for 86 NDBS samples from the California Biobank Program: (1) 45 DS-CHD (27 female, 18 male) and (2) 41 DS non-CHD (27 female, 14 male). We analyzed global CpG methylation and identified differentially methylated regions (DMRs) in DS-CHD versus DS non-CHD comparisons (both sex-combined and sex-stratified) corrected for sex, age of blood collection, and cell-type proportions. CHD DMRs were analyzed for enrichment in CpG and genic contexts, chromatin states, and histone modifications by genomic coordinates and for gene ontology enrichment by gene mapping. DMRs were also tested in a replication dataset and compared to methylation levels in DS versus typical development (TD) WGBS NDBS samples.ResultsWe found global CpG hypomethylation in DS-CHD males compared to DS non-CHD males, which was attributable to elevated levels of nucleated red blood cells and not seen in females. At a regional level, we identified 58, 341, and 3938 CHD-associated DMRs in the Sex Combined, Females Only, and Males Only groups, respectively, and used machine learning algorithms to select 19 Males Only loci that could distinguish CHD from non-CHD. DMRs in all comparisons were enriched for gene exons, CpG islands, and bivalent chromatin and mapped to genes enriched for terms related to cardiac and immune functions. Lastly, a greater percentage of CHD-associated DMRs than background regions were differentially methylated in DS versus TD samples.ConclusionsA sex-specific signature of DNA methylation was detected in NDBS of DS-CHD compared to DS non-CHD individuals. This supports the hypothesis that epigenetics can reflect the variability of phenotypes in DS, particularly CHDs.

BackgroundGenetic variants can modulate phenotypic outcomes via epigenetic intermediates, for example at methylation quantitative trait loci (mQTL). We present the first large-scale assessment of mQTL at human genomic regions selected for interindividual variation in CpG methylation, which we call correlated regions of systemic interindividual variation (CoRSIVs). These can be assayed in blood DNA and do not reflect interindividual variation in cellular composition.ResultsWe use target-capture bisulfite sequencing to assess DNA methylation at 4086 CoRSIVs in multiple tissues from each of 188 donors in the NIH Gene-Tissue Expression (GTEx) program. At CoRSIVs, DNA methylation in peripheral blood correlates with methylation and gene expression in internal organs. We also discover unprecedented mQTL at these regions. Genetic influences on CoRSIV methylation are extremely strong (median R2=0.76), cumulatively comprising over 70-fold more human mQTL than detected in the most powerful previous study. Moreover, mQTL beta coefficients at CoRSIVs are highly skewed (i.e., the major allele predicts higher methylation). Both surprising findings are independently validated in a cohort of 47 non-GTEx individuals. Genomic regions flanking CoRSIVs show long-range enrichments for LINE-1 and LTR transposable elements; the skewed beta coefficients may therefore reflect evolutionary selection of genetic variants that promote their methylation and silencing. Analyses of GWAS summary statistics show that mQTL polymorphisms at CoRSIVs are associated with metabolic and other classes of disease.ConclusionsA focus on systemic interindividual epigenetic variants, clearly enhanced in mQTL content, should likewise benefit studies attempting to link human epigenetic variation to the risk of disease.

Sufficient evidence supports a relationship between certain myeloid neoplasms and exposure to benzene or formaldehyde. DNA methylation could underlie benzene- and formaldehyde-induced health outcomes, but data in exposed human populations are limited. We conducted two cross-sectional epigenome-wide association studies (EWAS), one in workers exposed to benzene and another in workers exposed to formaldehyde. Using HumanMethylation450 BeadChips, we investigated differences in blood cell DNA methylation among 50 benzene-exposed subjects and 48 controls, and among 31 formaldehyde-exposed subjects and 40 controls. We performed CpG-level and regional-level analyses. In the benzene EWAS, we found genome-wide significant alterations, i.e., FWER-controlled P-values

Epigenome-wide association studies (EWAS) are widely implemented in epidemiology, and the Illumina HumanMethylationEPIC BeadChip (EPIC) DNA microarray is the most-used technology. Recently, next-generation sequencing (NGS)-based methods, which assess DNA methylation at single-base resolution, have become more affordable and technically feasible. While the content of microarray technology is fixed, NGS-based approaches, such as the Roche Nimblegen, SeqCap Epi Enrichment System (SeqCap), offer the flexibility of targeting most CpGs in a gene. With the current usage of microarrays and emerging NGS-based technologies, it is important to establish whether data generated from the two platforms are comparable. We harnessed 112 samples from the Center for the Health Assessment of Mothers and Children of Salinas (CHAMACOS) birth cohort study and compared DNA methylation between the EPIC microarray and SeqCap for PON1 and nine additional candidate genes, by evaluating epigenomic coverage and correlations. We conducted multivariable linear regression and principal component analyses to assess the ability of the EPIC array and SeqCap to detect biological differences in gene methylation by the PON1-108 single nucleotide polymorphism. We found an overall high concordance (r = 0.84) between SeqCap and EPIC DNA methylation, among highly methylated and minimally methylated regions. However, substantial disagreement was present between the two methods in moderately methylated regions, with SeqCap measurements exhibiting greater within-site variation. Additionally, SeqCap did not capture PON1 SNP associated differences in DNA methylation that were evident with the EPIC array. Our findings indicate that microarrays perform well for analysing DNA methylation in large cohort studies but with limited coverage.

BackgroundIdentifying blood-based DNA methylation patterns is a minimally invasive way to detect biomarkers in predicting age, characteristics of certain diseases and conditions, as well as responses to immunotherapies. As microarray platforms continue to evolve and increase the scope of CpGs measured, new discoveries based on the most recent platform version and how they compare to available data from the previous versions of the platform are unknown. The neutrophil dexamethasone methylation index (NDMI 850) is a blood-based DNA methylation biomarker built on the Illumina MethylationEPIC (850K) array that measures epigenetic responses to dexamethasone (DEX), a synthetic glucocorticoid often administered for inflammation. Here, we compare the NDMI 850 to one we built using data from the Illumina Methylation 450K (NDMI 450).ResultsThe NDMI 450 consisted of 22 loci, 15 of which were present on the NDMI 850. In adult whole blood samples, the linear composite scores from NDMI 450 and NDMI 850 were highly correlated and had equivalent predictive accuracy for detecting DEX exposure among adult glioma patients and non-glioma adult controls. However, the NDMI 450 scores of newborn cord blood were significantly lower than NDMI 850 in samples measured with both assays.ConclusionsWe developed an algorithm that reproduces the DNA methylation glucocorticoid response score using 450K data, increasing the accessibility for researchers to assess this biomarker in archived or publicly available datasets that use the 450K version of the Illumina BeadChip array. However, the NDMI850 and NDMI450 do not give similar results in cord blood, and due to data availability limitations, results from sample types of newborn cord blood should be interpreted with care.

BackgroundMyeloid-derived suppressor cells (MDSCs), which include monocytic (mMDSCs) and granulocytic (gMDSCs) cells, are an immunosuppressive, heterogeneous population of cells upregulated in cancer and other pathologic conditions, in addition to normal conditions of stress. The origin of MDSCs is debated, and the regulatory pattern responsible for gMDSC differentiation remains unknown. Since DNA methylation (DNAm) contributes to lineage differentiation, we have investigated whether it contributes to the acquisition of the gMDSC phenotype.ResultsUsing the Illumina EPIC array to measure DNAm of gMDSCs and neutrophils from diverse neonatal and adult blood sources, we found 189 differentially methylated CpGs between gMDSCs and neutrophils with a core of ten differentially methylated CpGs that were consistent across both sources of cells. Genes associated with these loci that are involved in immune responses include VCL, FATS, YAP1, KREMEN2, UBTF, MCC-1, and EFCC1. In two cancer patient groups that reflected those used to develop the methylation markers (head and neck squamous cell carcinoma (HNSCC) and glioma), all of the CpG loci were differentially methylated, reaching statistical significance in glioma cases and controls, while one was significantly different in the smaller HNSCC group.ConclusionsOur findings indicate that gMDSCs have a core of distinct DNAm alterations, informing future research on gMDSC differentiation and function.

ObjectiveFindings from cross-sectional studies have revealed associations between DNA methylation and systemic lupus erythematosus (SLE) outcomes. This study was undertaken to investigate the dynamics of DNA methylation by examining participants from an SLE longitudinal cohort using samples collected at 2 time points.MethodsA total of 101 participants from the California Lupus Epidemiology Study were included in our analysis. DNA was extracted from blood samples collected at the time of enrolment in the cohort and samples collected after 2 years and was analyzed using Illumina EPIC BeadChip kit. Paired t-tests were used to identify genome-wide changes which included 256 CpG sites previously found to be associated with SLE subtypes. Linear mixed models were developed to understand the relationship between DNA methylation and disease activity, medication use, and sample cell-type proportions, adjusted for age, sex, and genetic principal components.ResultsThe majority of CpGs that were previously determined to be associated with SLE subtypes remained stable over 2 years (185 CpGs [72.3%]; t-test false discovery rate >0.05). Compared to background genome-wide methylation, there was an enrichment of SLE subtype-associated CpGs that changed over time (27.7% versus 0.34%). Changes in cell-type proportions were associated with changes at 67 CpGs (P

BackgroundLRIG1, the founding member of the LRIG (leucine-rich repeat and immunoglobulin-like domain) family of transmembrane proteins, is a negative regulator of receptor tyrosine kinases and a tumour suppressor. Decreased LRIG1 expression is consistently observed in cancer, across diverse tumour types, and is linked to poor patient prognosis. However, mechanisms by which LRIG1 is repressed are not fully understood. Silencing of LRIG1 through promoter CpG island methylation has been reported in colorectal and cervical cancer but studies in breast cancer remain limited.MethodsIn silico analysis of human breast cancer patient data were used to demonstrate a correlation between DNA methylation and LRIG1 silencing in basal/triple-negative breast cancer, and its impact on patient survival. LRIG1 gene expression, protein abundance, and methylation enrichment were examined by quantitative reverse-transcription PCR, immunoblotting, and methylation immunoprecipitation, respectively, in breast cancer cell lines in vitro. We examined the impact of global demethylation on LRIG1 expression and methylation enrichment using 5-aza-2'-deoxycytidine. We also examined the effects of targeted demethylation of the LRIG1 CpG island, and transcriptional activation of LRIG1 expression, using the RNA guided deadCas9 transactivation system.ResultsAcross breast cancer subtypes, LRIG1 expression is lowest in the basal/triple-negative subtype so we investigated whether differential methylation may contribute to this. Indeed, we find that LRIG1 CpG island methylation is most prominent in basal/triple-negative cell lines and patient samples. Use of the global demethylating agent 5-aza-2'-deoxycytidine decreases methylation leading to increased LRIG1 transcript expression in basal/triple-negative cell lines, while having no effect on LRIG1 expression in luminal/ER-positive cell lines. Using a CRISPR/deadCas9 (dCas9)-based targeting approach, we demonstrate that TET1-mediated demethylation (Tet1-dCas9) along with VP64-mediated transcriptional activation (VP64-dCas9) at the CpG island, increased endogenous LRIG1 expression in basal/triple-negative breast cancer cells, without transcriptional upregulation at predicted off-target sites. Activation of LRIG1 by the dCas9 transactivation system significantly increased LRIG1 protein abundance, reduced site-specific methylation, and reduced cancer cell viability. Our findings suggest that CRISPR-mediated targeted activation may be a feasible way to restore LRIG1 expression in cancer.ConclusionsOur study contributes novel insight into mechanisms which repress LRIG1 in triple-negative breast cancer and demonstrates for the first time that targeted de-repression of LRIG1 in cancer cells is possible. Understanding the epigenetic mechanisms associated with repression of tumour suppressor genes holds potential for the advancement of therapeutic approaches.

Although male-female differences in placental structure and function have been observed, little is understood about their molecular underpinnings. Here, we present a mega-analysis of 14 publicly available placenta DNA methylation (DNAm) microarray datasets to identify individual CpGs and regions associated with fetal sex. In the discovery dataset of placentas from full term pregnancies (N = 532 samples), 5212 CpGs met genome-wide significance (p

BackgroundAltered DNA methylation may be an intermediate phenotype between breast cancer risk factors and disease. Mammographic density is a strong risk factor for breast cancer. However, no studies to date have identified an epigenetic signature of mammographic density. We performed an epigenome-wide association study of mammographic density.MethodsWhite blood cell DNA methylation was measured for 385 postmenopausal women using the Illumina Infinium MethylationEPIC BeadChip array. Differential methylation was assessed using genome-wide, probe-level, and regional analyses. We implemented a resampling-based approach to improve the stability of our findings.ResultsOn average, women with elevated mammographic density exhibited DNA hypermethylation within CpG islands and gene promoters compared to women with lower mammographic density. We identified 250 CpG sites for which DNA methylation was significantly associated with mammographic density. The top sites were located within genes associated with cancer, including HDLBP, TGFB2, CCT4, and PAX8, and were more likely to be located in regulatory regions of the genome. We also identified differential DNA methylation in 37 regions, including within the promoters of PAX8 and PF4, a gene involved in the regulation of angiogenesis. Overall, our results paint a picture of epigenetic dysregulation associated with mammographic density.ConclusionMammographic density is associated with differential DNA methylation throughout the genome, including within genes associated with cancer. Our results suggest the potential involvement of several genes in the biological mechanisms behind differences in breast density between women. Further studies are warranted to explore these potential mechanisms and potential links to breast cancer risk.

DNA methylation is a labile modification associated with gene expression control and environmental adaptations. High throughput, scalable and quantitative assessments of specific DNA methylation modifications in complex genomic regions for use in large population studies are needed. The performance of Droplet Digital™ PCR (ddPCR™) was investigated for DNA methylation detection against next-generation bisulfite sequencing (NGS) to demonstrate the ability of ddPCR to detect and validate DNA methylation levels and complex patterns among neighboring CpGs in regions associated with prenatal tobacco exposure. While both techniques are reproducible, ddPCR demonstrates a unique advantage for high-throughput DNA methylation analysis in large-scale population studies and provides the specificity to accurately measure DNA methylation of target CpGs in complex regions.

The completion of a telomere-to-telomere human reference genome, T2T-CHM13, has resolved complex regions of the genome, including repetitive and homologous regions. Here, we present a high-resolution epigenetic study of previously unresolved sequences, representing entire acrocentric chromosome short arms, gene family expansions, and a diverse collection of repeat classes. This resource precisely maps CpG methylation (32.28 million CpGs), DNA accessibility, and short-read datasets (166,058 previously unresolved chromatin immunoprecipitation sequencing peaks) to provide evidence of activity across previously unidentified or corrected genes and reveals clinically relevant paralog-specific regulation. Probing CpG methylation across human centromeres from six diverse individuals generated an estimate of variability in kinetochore localization. This analysis provides a framework with which to investigate the most elusive regions of the human genome, granting insights into epigenetic regulation.

Cognitive skills are a strong predictor of a wide range of later life outcomes. Genetic and epigenetic associations across the genome explain some of the variation in general cognitive abilities in the general population and it is plausible that epigenetic associations might arise from prenatal environmental exposures and/or genetic variation early in life. We investigated the association between cord blood DNA methylation at birth and cognitive skills assessed in children from eight pregnancy cohorts within the Pregnancy And Childhood Epigenetics (PACE) Consortium across overall (total N = 2196), verbal (total N = 2206) and non-verbal cognitive scores (total N = 3300). The associations at single CpG sites were weak for all of the cognitive domains investigated. One region near DUSP22 on chromosome 6 was associated with non-verbal cognition in a model adjusted for maternal IQ. We conclude that there is little evidence to support the idea that variation in cord blood DNA methylation at single CpG sites is associated with cognitive skills and further studies are needed to confirm the association at DUSP22.

Health outcomes are frequently shaped by difficult to dissect inter-relationships between biological, behavioral, social and environmental factors. DNA methylation patterns reflect such multivariate intersections, providing a rich source of novel biomarkers and insight into disease etiologies. Recent advances in whole-genome bisulfite sequencing enable investigation of DNA methylation over all genomic CpGs, but existing bioinformatic approaches lack accessible system-level tools. Here, we develop the R package Comethyl, for weighted gene correlation network analysis of user-defined genomic regions that generates modules of comethylated regions, which are then tested for correlations with multivariate sample traits. First, regions are defined by CpG genomic location or regulatory annotation and filtered based on CpG count, sequencing depth and variability. Next, correlation networks are used to find modules of interconnected nodes using methylation values within the selected regions. Each module containing multiple comethylated regions is reduced in complexity to a single eigennode value, which is then tested for correlations with experimental metadata. Comethyl has the ability to cover the noncoding regulatory regions of the genome with high relevance to interpretation of genome-wide association studies and integration with other types of epigenomic data. We demonstrate the utility of Comethyl on a dataset of male cord blood samples from newborns later diagnosed with autism spectrum disorder (ASD) versus typical development. Comethyl successfully identified an ASD-associated module containing regions mapped to genes enriched for brain glial functions. Comethyl is expected to be useful in uncovering the multivariate nature of health disparities for a variety of common disorders. Comethyl is available at github.com/cemordaunt/comethyl with complete documentation and example analyses.

Type 1 diabetes affects over nine million individuals globally, with approximately 40% developing diabetic kidney disease. Emerging evidence suggests that epigenetic alterations, such as DNA methylation, are involved in diabetic kidney disease. Here we assess differences in blood-derived genome-wide DNA methylation associated with diabetic kidney disease in 1304 carefully characterised individuals with type 1 diabetes and known renal status from two cohorts in the United Kingdom-Republic of Ireland and Finland. In the meta-analysis, we identify 32 differentially methylated CpGs in diabetic kidney disease in type 1 diabetes, 18 of which are located within genes differentially expressed in kidneys or correlated with pathological traits in diabetic kidney disease. We show that methylation at 21 of the 32 CpGs predict the development of kidney failure, extending the knowledge and potentially identifying individuals at greater risk for diabetic kidney disease in type 1 diabetes.

BackgroundOvarian clear cell carcinoma (OCCC) is a rare ovarian cancer histotype that tends to be resistant to standard platinum-based chemotherapeutics. We sought to better understand the role of DNA methylation in clinical and biological subclassification of OCCC.MethodsWe interrogated genome-wide methylation using DNA from fresh frozen tumors from 271 cases, applied nonsmooth nonnegative matrix factorization (nsNMF) clustering, and evaluated clinical associations and biological pathways.ResultsTwo approximately equally sized clusters that associated with several clinical features were identified. Compared with Cluster 2 (N = 137), Cluster 1 cases (N = 134) presented at a more advanced stage, were less likely to be of Asian ancestry, and tended to have poorer outcomes including macroscopic residual disease following primary debulking surgery (P < 0.10). Subset analyses of targeted tumor sequencing and IHC data revealed that Cluster 1 tumors showed TP53 mutation and abnormal p53 expression, and Cluster 2 tumors showed aneuploidy and ARID1A/PIK3CA mutation (P < 0.05). Cluster-defining CpGs included 1,388 CpGs residing within 200 bp of the transcription start sites of 977 genes; 38% of these genes (N = 369 genes) were differentially expressed across cluster in transcriptomic subset analysis (P < 10-4). Differentially expressed genes were enriched for six immune-related pathways, including IFNα and IFNγ responses (P < 10-6).ConclusionsDNA methylation clusters in OCCC correlate with disease features and gene expression patterns among immune pathways.ImpactThis work serves as a foundation for integrative analyses that better understand the complex biology of OCCC in an effort to improve potential for development of targeted therapeutics.

BackgroundBeef tenderness is a complex trait of economic importance for the beef industry. Understanding the epigenetic mechanisms underlying this trait may help improve the accuracy of breeding programs. However, little is known about epigenetic effects on Bos taurus muscle and their implications in tenderness, and no studies have been conducted in Bos indicus.ResultsComparing methylation profile of Bos indicus skeletal muscle with contrasting beef tenderness at 14 days after slaughter, we identified differentially methylated cytosines and regions associated with this trait. Interestingly, muscle that became tender beef had higher levels of hypermethylation compared to the tough group. Enrichment analysis of predicted target genes suggested that differences in methylation between tender and tough beef may affect signal transduction pathways, among which G protein signaling was a key pathway. In addition, different methylation levels were found associated with expression levels of GNAS, PDE4B, EPCAM and EBF3 genes. The differentially methylated elements correlated with EBF3 and GNAS genes overlapped CpG islands and regulatory elements. GNAS, a complex imprinted gene, has a key role on G protein signaling pathways. Moreover, both G protein signaling pathway and the EBF3 gene regulate muscle homeostasis, relaxation, and muscle cell-specificity.ConclusionsWe present differentially methylated loci that may be of interest to decipher the epigenetic mechanisms affecting tenderness. Supported by the previous knowledge about regulatory elements and gene function, the methylation data suggests EBF3 and GNAS as potential candidate genes and G protein signaling as potential candidate pathway associated with beef tenderness via methylation.

DNA methylation microarrays can be employed to interrogate cell-type composition in complex tissues. Here, we expand reference-based deconvolution of blood DNA methylation to include 12 leukocyte subtypes (neutrophils, eosinophils, basophils, monocytes, naïve and memory B cells, naïve and memory CD4 + and CD8 + T cells, natural killer, and T regulatory cells). Including derived variables, our method provides 56 immune profile variables. The IDOL (IDentifying Optimal Libraries) algorithm was used to identify libraries for deconvolution of DNA methylation data for current and previous platforms. The accuracy of deconvolution estimates obtained using our enhanced libraries was validated using artificial mixtures and whole-blood DNA methylation with known cellular composition from flow cytometry. We applied our libraries to deconvolve cancer, aging, and autoimmune disease datasets. In conclusion, these libraries enable a detailed representation of immune-cell profiles in blood using only DNA and facilitate a standardized, thorough investigation of immune profiles in human health and disease.

Infinium methylation arrays are not available for the vast majority of non-human mammals. Moreover, even if species-specific arrays were available, probe differences between them would confound cross-species comparisons. To address these challenges, we developed the mammalian methylation array, a single custom array that measures up to 36k CpGs per species that are well conserved across many mammalian species. We designed a set of probes that can tolerate specific cross-species mutations. We annotate the array in over 200 species and report CpG island status and chromatin states in select species. Calibration experiments demonstrate the high fidelity in humans, rats, and mice. The mammalian methylation array has several strengths: it applies to all mammalian species even those that have not yet been sequenced, it provides deep coverage of conserved cytosines facilitating the development of epigenetic biomarkers, and it increases the probability that biological insights gained in one species will translate to others.

BackgroundFew studies have examined prognostic outcomes-associated molecular signatures other than overall survival (OS) for gastric cancer (GC). We aimed to identify DNA methylation biomarkers associated with multiple prognostic outcomes of GC in an epigenome-wide association study.MethodsBased on the Cancer Genome Atlas (TCGA), DNA methylation loci associated with OS (n = 381), disease-specific survival (DSS, n = 372), and progression-free interval (PFI, n = 383) were discovered in training set subjects (false discovery rates

BackgroundObsessive-compulsive disorder (OCD) is characterized by intrusive thoughts and repetitive actions, that presents the involvement of the cortico-striatal areas. The contribution of environmental risk factors to OCD development suggests that epigenetic mechanisms may contribute to its pathophysiology. DNA methylation changes and gene expression were evaluated in post-mortem brain tissues of the cortical (anterior cingulate gyrus and orbitofrontal cortex) and ventral striatum (nucleus accumbens, caudate nucleus and putamen) areas from eight OCD patients and eight matched controls.ResultsThere were no differentially methylated CpG (cytosine-phosphate-guanine) sites (DMSs) in any brain area, nevertheless gene modules generated from CpG sites and protein-protein-interaction (PPI) showed enriched gene modules for all brain areas between OCD cases and controls. All brain areas but nucleus accumbens presented a predominantly hypomethylation pattern for the differentially methylated regions (DMRs). Although there were common transcriptional factors that targeted these DMRs, their targeted differentially expressed genes were different among all brain areas. The protein-protein interaction network based on methylation and gene expression data reported that all brain areas were enriched for G-protein signaling pathway, immune response, apoptosis and synapse biological processes but each brain area also presented enrichment of specific signaling pathways. Finally, OCD patients and controls did not present significant DNA methylation age differences.ConclusionsDNA methylation changes in brain areas involved with OCD, especially those involved with genes related to synaptic plasticity and the immune system could mediate the action of genetic and environmental factors associated with OCD.

BackgroundCardiovascular health (CVH) has been defined by the American Heart Association (AHA) as the presence of the "Life's Simple 7" ideal lifestyle and clinical factors. CVH is known to predict longevity and freedom from cardiovascular disease, the leading cause of death for women in the United States. DNA methylation markers of aging have been aggregated into a composite epigenetic age score, which is associated with cardiovascular morbidity and mortality. However, it is unknown whether poor CVH is associated with acceleration of aging as measured by DNA methylation markers in epigenetic age.Methods and resultsWe performed a cross-sectional analysis of racially/ethnically diverse post-menopausal women enrolled in the Women's Health Initiative cohort recruited between 1993 and 1998. Epigenetic age acceleration (EAA) was calculated using DNA methylation data on a subset of participants and the published Horvath and Hannum methods for intrinsic and extrinsic EAA. CVH was calculated using the AHA measures of CVH contributing to a 7-point score. We examined the association between CVH score and EAA using linear regression modeling adjusting for self-reported race/ethnicity and education. Among the 2,170 participants analyzed, 50% were white and mean age was 64 (7 SD) years. Higher or more favorable CVH scores were associated with lower extrinsic EAA (~ 6 months younger age per 1 point higher CVH score, p

BackgroundThe neuropeptide oxytocin regulates mammalian social behavior. Disruptions in oxytocin signaling are a feature of many psychopathologies. One commonly studied biomarker for oxytocin involvement in psychiatric diseases is DNA methylation at the oxytocin receptor gene (OXTR). Such studies focus on DNA methylation in two regions of OXTR, exon 3 and a region termed MT2 which overlaps exon 1 and intron 1. However, the relative contribution of exon 3 and MT2 in regulating OXTR gene expression in the brain is currently unknown.ResultsHere, we use the prairie vole as a translational animal model to investigate genetic, epigenetic, and environmental factors affecting Oxtr gene expression in a region of the brain that has been shown to drive Oxtr related behavior in the vole, the nucleus accumbens. We show that the genetic structure of Oxtr in prairie voles resembles human OXTR. We then studied the effects of early life experience on DNA methylation in two regions of a CpG island surrounding the Oxtr promoter: MT2 and exon 3. We show that early nurture in the form of parental care results in DNA hypomethylation of Oxtr in both MT2 and exon 3, but only DNA methylation in MT2 is associated with Oxtr gene expression. Network analyses indicate that CpG sites in the 3' portion of MT2 are most highly associated with Oxtr gene expression. We also identify two novel SNPs in exon 3 of Oxtr in prairie voles and a novel alternative transcript originating from the third intron of the gene. Expression of the novel alternative transcript is associated with genotype at SNP KLW2.ConclusionsThese results identify putative regulatory features of Oxtr in prairie voles which inform future studies examining OXTR in human social behaviors and disorders. These studies indicate that in prairie voles, DNA methylation in MT2, particularly in the 3' portion, is more predictive of Oxtr gene expression than DNA methylation in exon 3. Similarly, in human temporal cortex, we find that DNA methylation in the 3' portion of MT2 is associated with OXTR expression. Together, these results suggest that among the CpG sites studied, DNA methylation of MT2 may be the most reliable indicator of OXTR gene expression. We also identify novel features of prairie vole Oxtr, including SNPs and an alternative transcript, which further develop the prairie vole as a translational model for studies of OXTR.

BackgroundDNA methylation (DNAm) is associated with gene regulation and estimated glomerular filtration rate (eGFR), a measure of kidney function. Decreased eGFR is more common among US Hispanics and African Americans. The causes for this are poorly understood. We aimed to identify trans-ethnic and ethnic-specific differentially methylated positions (DMPs) associated with eGFR using an agnostic, genome-wide approach.MethodsThe study included up to 5428 participants from multi-ethnic studies for discovery and 8109 participants for replication. We tested the associations between whole blood DNAm and eGFR using beta values from Illumina 450K or EPIC arrays. Ethnicity-stratified analyses were performed using linear mixed models adjusting for age, sex, smoking, and study-specific and technical variables. Summary results were meta-analyzed within and across ethnicities. Findings were assessed using integrative epigenomics methods and pathway analyses.ResultsWe identified 93 DMPs associated with eGFR at an FDR of 0.05 and replicated 13 and 1 DMPs across independent samples in trans-ethnic and African American meta-analyses, respectively. The study also validated 6 previously published DMPs. Identified DMPs showed significant overlap enrichment with DNase I hypersensitive sites in kidney tissue, sites associated with the expression of proximal genes, and transcription factor motifs and pathways associated with kidney tissue and kidney development.ConclusionsWe uncovered trans-ethnic and ethnic-specific DMPs associated with eGFR, including DMPs enriched in regulatory elements in kidney tissue and pathways related to kidney development. These findings shed light on epigenetic mechanisms associated with kidney function, bridging the gap between population-specific eGFR-associated DNAm and tissue-specific regulatory context.

BackgroundBiological aging estimators derived from DNA methylation data are heritable and correlate with morbidity and mortality. Consequently, identification of genetic and environmental contributors to the variation in these measures in populations has become a major goal in the field.ResultsLeveraging DNA methylation and SNP data from more than 40,000 individuals, we identify 137 genome-wide significant loci, of which 113 are novel, from genome-wide association study (GWAS) meta-analyses of four epigenetic clocks and epigenetic surrogate markers for granulocyte proportions and plasminogen activator inhibitor 1 levels, respectively. We find evidence for shared genetic loci associated with the Horvath clock and expression of transcripts encoding genes linked to lipid metabolism and immune function. Notably, these loci are independent of those reported to regulate DNA methylation levels at constituent clock CpGs. A polygenic score for GrimAge acceleration showed strong associations with adiposity-related traits, educational attainment, parental longevity, and C-reactive protein levels.ConclusionThis study illuminates the genetic architecture underlying epigenetic aging and its shared genetic contributions with lifestyle factors and longevity.

DNA methylation may be regulated by genetic variants within a genomic region, referred to as methylation quantitative trait loci (mQTLs). The changes of methylation levels can further lead to alterations of gene expression, and influence the risk of various complex human diseases. Detecting mQTLs may provide insights into the underlying mechanism of how genotypic variations may influence the disease risk. In this article, we propose a methylation random field (MRF) method to detect mQTLs by testing the association between the methylation level of a CpG site and a set of genetic variants within a genomic region. The proposed MRF has two major advantages over existing approaches. First, it uses a beta distribution to characterize the bimodal and interval properties of the methylation trait at a CpG site. Second, it considers multiple common and rare genetic variants within a genomic region to identify mQTLs. Through simulations, we demonstrated that the MRF had improved power over other existing methods in detecting rare variants of relatively large effect, especially when the sample size is small. We further applied our method to a study of congenital heart defects with 83 cardiac tissue samples and identified two mQTL regions, MRPS10 and PSORS1C1, which were colocalized with expression QTL in cardiac tissue. In conclusion, the proposed MRF is a useful tool to identify novel mQTLs, especially for studies with limited sample sizes.

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.

Biological mechanisms underlying human germline mutations remain largely unknown. We statistically decompose variation in the rate and spectra of mutations along the genome using volume-regularized nonnegative matrix factorization. The analysis of a sequencing dataset (TOPMed) reveals nine processes that explain the variation in mutation properties between loci. We provide a biological interpretation for seven of these processes. We associate one process with bulky DNA lesions that are resolved asymmetrically with respect to transcription and replication. Two processes track direction of replication fork and replication timing, respectively. We identify a mutagenic effect of active demethylation primarily acting in regulatory regions and a mutagenic effect of long interspersed nuclear elements. We localize a mutagenic process specific to oocytes from population sequencing data. This process appears transcriptionally asymmetric.

BackgroundParental smoking is implicated in the etiology of acute lymphoblastic leukemia (ALL), the most common childhood cancer. We recently reported an association between an epigenetic biomarker of early-life tobacco smoke exposure at the AHRR gene and increased frequency of somatic gene deletions among ALL cases.MethodsHere, we further assess this association using two epigenetic biomarkers for maternal smoking during pregnancy-DNA methylation at AHRR CpG cg05575921 and a recently established polyepigenetic smoking score-in an expanded set of 482 B-cell ALL (B-ALL) cases in the California Childhood Leukemia Study with available Illumina 450K or MethylationEPIC array data. Multivariable Poisson regression models were used to test the associations between the epigenetic biomarkers and gene deletion numbers.ResultsWe found an association between DNA methylation at AHRR CpG cg05575921 and deletion number among 284 childhood B-ALL cases with MethylationEPIC array data, with a ratio of means (RM) of 1.31 [95% confidence interval (CI), 1.02-1.69] for each 0.1 β value reduction in DNA methylation, an effect size similar to our previous report in an independent set of 198 B-ALL cases with 450K array data [meta-analysis summary RM (sRM) = 1.32; 95% CI, 1.10-1.57]. The polyepigenetic smoking score was positively associated with gene deletion frequency among all 482 B-ALL cases (sRM = 1.31 for each 4-unit increase in score; 95% CI, 1.09-1.57).ConclusionsWe provide further evidence that prenatal tobacco-smoke exposure may influence the generation of somatic copy-number deletions in childhood B-ALL.ImpactAnalyses of deletion breakpoint sequences are required to further understand the mutagenic effects of tobacco smoke in childhood ALL.

BackgroundEstrogens are thought to contribute to breast cancer risk through cell cycling and accelerated breast aging. We hypothesize that lifetime estrogen exposure drives early epigenetic breast aging observed in healthy women. In this study, we examined associations between hormonal factors and epigenetic aging measures in healthy breast tissues.MethodsWe extracted DNA from breast tissue specimens from 192 healthy female donors to the Susan G. Komen Tissue Bank at the Indiana University Simon Cancer Center. Methylation experiments were performed using the Illumina EPIC 850K array platform. Age-adjusted regression models were used to examine for associations between factors related to estrogen exposure and five DNA methylation-based estimates: Grim age, pan-tissue age, Hannum age, phenotypic age, and skin and blood clock age.ResultsWomen were aged 19-90 years, with 95 premenopausal, and 97 nulliparous women. The age difference (Grim age - chronologic age) was higher at earlier ages close to menarche. We found significant associations between earlier age at menarche and age-adjusted accelerations according to the Grim clock, the skin and blood clock, and between higher body mass index (BMI) and age-adjusted accelerations in the Grim clock, Hannum clock, phenotypic clock, and skin and blood clock.ConclusionsEarlier age at menarche and higher BMI are associated with elevations in DNA methylation-based age estimates in healthy breast tissues, suggesting that cumulative estrogen exposure drives breast epigenetic aging.ImpactEpigenetic clock measures may help advance inquiry into the relationship between accelerated breast tissue aging and an elevated incidence of breast cancer in younger women.

BackgroundDiet quality is a risk factor for chronic disease and mortality. Differential DNA methylation across the epigenome has been associated with chronic disease risk. Whether diet quality is associated with differential methylation is unknown. This study assessed whether diet quality was associated with differential DNA methylation measured across 445 548 loci in the Women's Health Initiative (WHI) and the TwinsUK cohort.DesignThe discovery cohort consisted of 4355 women from the WHI. The replication cohort consisted of 571 mono- and dizygotic twins from the TwinsUK cohort. DNA methylation was measured in whole blood using the Illumina Infinium HumanMethylation450 Beadchip. Diet quality was assessed using the Alternative Healthy Eating Index 2010 (AHEI-2010). A meta-analysis, stratified by study cohort, was performed using generalized linear models that regressed methylation on AHEI-2010, adjusting for cell composition, chip number and location, study characteristics, principal components of genetic relatedness, age, smoking status, race/ethnicity and body mass index (BMI). Statistical significance was defined as a false discovery rate < 0.05. Significant sites were tested for replication in the TwinsUK cohort, with significant replication defined by P

CpG islands often exhibit low DNA methylation, high histone H3 lysine 4 trimethylation, low nucleosome density, and high DNase I hypersensitivity, yet the rules by which CpG islands are sensed remain poorly understood. In this study, we first evaluated the relationships between the DNA and the chromatin properties of CpG islands in embryonic stem cells using modified bacterial artificial chromosomes. Then, using a bioinformatic approach, we identified strict CpG-island density and length thresholds in mouse embryonic stem and differentiated cells that consistently specify low DNA methylation levels. Surprisingly, the human genome exhibited a dramatically different relationship between DNA properties and DNA methylation levels of CpG islands. Further analysis allowed speculation that this difference is accommodated in part by evolutionary changes in the nucleotide composition of orthologous promoters. Thus, a change in the rules by which CpG-island properties are sensed may have co-evolved with compensatory genome adaptation events during mammalian evolution.

A general approach for heritably altering gene expression has the potential to enable many discovery and therapeutic efforts. Here, we present CRISPRoff-a programmable epigenetic memory writer consisting of a single dead Cas9 fusion protein that establishes DNA methylation and repressive histone modifications. Transient CRISPRoff expression initiates highly specific DNA methylation and gene repression that is maintained through cell division and differentiation of stem cells to neurons. Pairing CRISPRoff with genome-wide screens and analysis of chromatin marks establishes rules for heritable gene silencing. We identify single guide RNAs (sgRNAs) capable of silencing the large majority of genes including those lacking canonical CpG islands (CGIs) and reveal a wide targeting window extending beyond annotated CGIs. The broad ability of CRISPRoff to initiate heritable gene silencing even outside of CGIs expands the canonical model of methylation-based silencing and enables diverse applications including genome-wide screens, multiplexed cell engineering, enhancer silencing, and mechanistic exploration of epigenetic inheritance.

Neonatal dried blood spots (NDBS) are a widely banked sample source that enables retrospective investigation into early life molecular events. Here, we performed low-pass whole genome bisulfite sequencing (WGBS) of 86 NDBS DNA to examine early life Down syndrome (DS) DNA methylation profiles. DS represents an example of genetics shaping epigenetics, as multiple array-based studies have demonstrated that trisomy 21 is characterized by genome-wide alterations to DNA methylation. By assaying over 24 million CpG sites, thousands of genome-wide significant (q

Increased expression of vascular cell adhesion molecule (VCAM)-1 on the activated arterial endothelial cell (EC) surface critically contributes to atherosclerosis which may in part be regulated by epigenetic mechanisms. This study investigated whether and how the clinically available histone deacetylases 1 and 2 (HDAC1/2) inhibitor drug Romidepsin epigenetically modulates VCAM-1 expression to suppress atherosclerosis. Methods: VCAM-1 expression was analyzed in primary human aortic EC (HAEC) treated with Romidepsin or transfected with HDAC1/2-targeting siRNA. Methylation of GATA6 promoter region was examined with methylation-specific PCR assay. Enrichment of STAT3 to GATA6 promoter was detected with chromatin immunoprecipitation. Lys685Arg mutation was constructed to block STAT3 acetylation. The potential therapeutic effect of Romidepsin on atherosclerosis was evaluated in Apoe -/- mice fed with a high-fat diet. Results: Romidepsin significantly attenuated TNFα-induced VCAM-1 expression on HAEC surface and monocyte adhesion through simultaneous inhibition of HDAC1/2. This downregulation of VCAM-1 was attributable to reduced expression of transcription factor GATA6. Romidepsin enhanced STAT3 acetylation and its binding to DNA methyltransferase 1 (DNMT1), leading to hypermethylation of the GATA6 promoter CpG-rich region at +140/+255. Blocking STAT3 acetylation at Lys685 disrupted DNMT1-STAT3 interaction, decreased GATA6 promoter methylation, and reversed the suppressive effects of HDAC1/2 inhibition on GATA6 and VCAM-1 expression. Finally, intraperitoneal administration of Romidepsin reduced diet-induced atherosclerotic lesion development in Apoe -/- mice, accompanied by a reduction in GATA6/VCAM-1 expression in the aorta. Conclusions: HDAC1/2 contributes to VCAM-1 expression and atherosclerosis by suppressing STAT3 acetylation-dependent GATA6 promoter methylation. These findings may provide a rationale for HDAC1/2-targeting therapy in atherosclerotic heart disease.

Most human cancers converge to a deregulated methylome with reduced global levels and elevated methylation at select CpG islands. To investigate the emergence and dynamics of the cancer methylome, we characterized genome-wide DNA methylation in pre-neoplastic monoclonal B cell lymphocytosis (MBL) and chronic lymphocytic leukemia (CLL), including serial samples collected across disease course. We detected the aberrant tumor-associated methylation landscape at CLL diagnosis and found no significantly differentially methylated regions in the high-count MBL-to-CLL transition. Patient methylomes showed remarkable stability with natural disease and post-therapy progression. Single CLL cells were consistently aberrantly methylated, indicating a homogeneous transition to the altered epigenetic state, and a distinct expression profile together with MBL cells compared to normal B cells. Our longitudinal analysis reveals the cancer methylome to emerge early, which may provide a platform for subsequent genetically-driven growth dynamics and together with its persistent presence suggests a central role in the normal-to-cancer transition.

PurposeKnown clinical and genetic markers have limitations in predicting disease course and outcome in juvenile myelomonocytic leukemia (JMML). DNA methylation patterns in JMML have correlated with outcome across multiple studies, suggesting it as a biomarker to improve patient stratification. However, standardized approaches to classify JMML on the basis of DNA methylation patterns are lacking. We, therefore, sought to define an international consensus for DNA methylation subgroups in JMML and develop classification methods for clinical implementation.Experimental designPublished DNA methylation data from 255 patients with JMML were used to develop and internally validate a classifier model. Accuracy across platforms (EPIC-arrays and MethylSeq) was tested using a technical validation cohort (32 patients). The suitability of both methods for single-patient classification was demonstrated using an independent cohort (47 patients).ResultsAnalysis of pooled, published data established three DNA methylation subgroups as a de facto standard. Unfavorable prognostic parameters (PTPN11 mutation, elevated fetal hemoglobin, and older age) were significantly enriched in the high methylation (HM) subgroup. A classifier was then developed that predicted subgroups with 98% accuracy across different technological platforms. Applying the classifier to an independent validation cohort confirmed an association of HM with secondary mutations, high relapse incidence, and inferior overall survival (OS), while the low methylation subgroup was associated with a favorable disease course. Multivariable analysis established DNA methylation subgroups as the only significant factor predicting OS.ConclusionsThis study provides an international consensus definition for DNA methylation subgroups in JMML. We developed and validated methods which will facilitate the design of risk-stratified clinical trials in JMML.

The mammalian DNA methylome is formed by two antagonizing processes, methylation by DNA methyltransferases (DNMT) and demethylation by ten-eleven translocation (TET) dioxygenases. Although the dynamics of either methylation or demethylation have been intensively studied in the past decade, the direct effects of their interaction on gene expression remain elusive. Here, we quantify the concurrence of DNA methylation and demethylation by the percentage of unmethylated CpGs within a partially methylated read from bisulfite sequencing. After verifying 'methylation concurrence' by its strong association with the co-localization of DNMT and TET enzymes, we observe that methylation concurrence is strongly correlated with gene expression. Notably, elevated methylation concurrence in tumors is associated with the repression of 40~60% of tumor suppressor genes, which cannot be explained by promoter hypermethylation alone. Furthermore, methylation concurrence can be used to stratify large undermethylated regions with negligible differences in average methylation into two subgroups with distinct chromatin accessibility and gene regulation patterns. Together, methylation concurrence represents a unique methylation metric important for transcription regulation and is distinct from conventional metrics, such as average methylation and methylation variation.

Down syndrome is associated with genome-wide perturbation of gene expression, which may be mediated by epigenetic changes. We perform an epigenome-wide association study on neonatal bloodspots comparing 196 newborns with Down syndrome and 439 newborns without Down syndrome, adjusting for cell-type heterogeneity, which identifies 652 epigenome-wide significant CpGs (P

BackgroundDown syndrome (DS) is characterized by a genome-wide profile of differential DNA methylation that is skewed towards hypermethylation in most tissues, including brain, and includes pan-tissue differential methylation. The molecular mechanisms involve the overexpression of genes related to DNA methylation on chromosome 21. Here, we stably overexpressed the chromosome 21 gene DNA methyltransferase 3L (DNMT3L) in the human SH-SY5Y neuroblastoma cell line and assayed DNA methylation at over 26 million CpGs by whole genome bisulfite sequencing (WGBS) at three different developmental phases (undifferentiated, differentiating, and differentiated).ResultsDNMT3L overexpression resulted in global CpG and CpG island hypermethylation as well as thousands of differentially methylated regions (DMRs). The DNMT3L DMRs were skewed towards hypermethylation and mapped to genes involved in neurodevelopment, cellular signaling, and gene regulation. Consensus DNMT3L DMRs showed that cell lines clustered by genotype and then differentiation phase, demonstrating sets of common genes affected across neuronal differentiation. The hypermethylated DNMT3L DMRs from all pairwise comparisons were enriched for regions of bivalent chromatin marked by H3K4me3 as well as differentially methylated sites from previous DS studies of diverse tissues. In contrast, the hypomethylated DNMT3L DMRs from all pairwise comparisons displayed a tissue-specific profile enriched for regions of heterochromatin marked by H3K9me3 during embryonic development.ConclusionsTaken together, these results support a mechanism whereby regions of bivalent chromatin that lose H3K4me3 during neuronal differentiation are targeted by excess DNMT3L and become hypermethylated. Overall, these findings demonstrate that DNMT3L overexpression during neurodevelopment recreates a facet of the genome-wide DS DNA methylation signature by targeting known genes and gene clusters that display pan-tissue differential methylation in DS.