mCSEA: detecting subtle differentially methylated regions

Host-microbe multi-omics and succinotype profiling have prognostic value for future relapse in patients with inflammatory bowel disease

Crohn's disease (CD) and ulcerative colitis (UC) are chronic relapsing inflammatory bowel disorders (IBD), the pathogenesis of which is uncertain but includes genetic susceptibility factors, immune-mediated tissue injury and environmental influences, most of which appear to act via the gut microbiome. We hypothesized that host-microbe alterations could be used to prognostically stratify patients experiencing relapses up to four years after endoscopy. We therefore examined multiple omics data, including published and new datasets, generated from paired inflamed and non-inflamed mucosal biopsies from 142 patients with IBD (54 CD; 88 UC) and from 34 control (non-diseased) biopsies. The relapse-predictive potential of 16S rRNA gene and transcript amplicons (standing and active microbiota) were investigated along with host transcriptomics, epigenomics and genetics. While standard single-omics analysis could not distinguish between patients who relapsed and those that remained in remission within four years of colonoscopy, we did find an association between the number of flares and a patient's succinotype. Our multi-omics machine learning approach was also able to predict relapse when combining features from the microbiome and human host. Therefore multi-omics, rather than single omics, better predicts relapse within 4 years of colonoscopy, while a patient's succinotype is associated with a higher frequency of relapses.

DNA methylome biomarkers of rheumatoid arthritis-associated interstitial lung disease reflecting lung fibrosis pathways, an exploratory case-control study

Rheumatoid Arthritis-associated Interstitial Lung Disease (RA-ILD) significantly reduces life quality and survival, necessitating improvements in its understanding and clinical management. We addressed these goals using DNA methylation analysis, which has not been done in RA-ILD samples, by comparing 32 RA patients with ILD diagnosed less than one year before (cases) and 32 matched RA patients without ILD (controls). This analysis identified 6679 differentially methylated positions (DMPs) with Δβ ≥ 2% and FDR < 0.05, and 576 differentially methylated regions in RA-ILD. Some DMPs were near mucin, collagen, and telomere maintenance genes. Also, the most notably enriched gene set (up to padj = 1.9 × 10-38) included genes overexpressed in fibrosis by monocytes and alveolar macrophages. Other significantly enriched gene sets, known to be dysregulated in fibrosis, included the mitotic spindle and the Rho GTPases. Additionally, analysis of transcription factor binding sites around DMPs showed unique enrichment near the liver X receptor element (LXRE), which is associated with fibrosis in multiple tissues. These results were consistent and unaffected by stricter significance thresholds. They indicated that differential DNA methylation may serve as blood biomarkers for RA-ILD including some related to lung fibrosis pathways.

Genome-wide methylation association study in monozygotic twins discordant for curve severity of adolescent idiopathic scoliosis

BACKGROUND CONTEXT

Emerging evidence suggests that abnormal DNA methylation patterns may play a role in the progression of adolescent idiopathic scoliosis (AIS). However, the mechanisms underlying the influence of DNA methylation on curve severity remain largely unknown.

PURPOSE

To characterize the DNA methylation profiles associated with the curve severity of AIS.

STUDY DESIGN

Retrospective study with prospectively collected clinical data and blood samples.

METHODS

A total of 7 AIS monozygotic twin pairs discordant for curve severity were included. Genome-wide methylation profile from blood samples were quantified by Illumina Infinium MethylationEPIC BeadChip (850K chip). Cell type composition of the samples was estimated by RefbaseEWAS method. Differentially methylated CpG sites were identified through comparison between patients with low and high Cobb angle. We also performed a gene-based collapsing analysis using mCSEA by aggregating the CpG sites based on promoter region. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed using clusterProfiler package.

RESULTS

Genome-wide DNA methylation analysis identified multiple differentially methylated positions across the genome. Gene-based collapsing analysis identified 212 differentially methylated genes (FDR adjusted p<.05), most of which (186/212) were hypermethylated in the group with high Cobb angle. Some of the identified genes were well-documented in AIS literature, such as TBX1, PAX3 and LBX1. Functional enrichment analysis revealed that the differentially methylated genes (DMGs) were involved in pattern specification process, skeletal development, muscle function, neurotransmission and several signaling pathways (cAMP, Wnt and prolactin).

CONCLUSIONS

The study represents the largest systematic epigenomic analyses of monozygotic twins discordant for curve severity and supports the important role of altered DNA methylation in AIS.

CLINICAL SIGNIFICANCE

The identified CpG sites provide insight into novel biomarkers predicting curve progression of AIS. Furthermore, the differentially methylated genes and enriched pathways may serve as interventional therapy target for AIS patients.

DMRIntTk: Integrating different DMR sets based on density peak clustering

BACKGROUND

Identifying differentially methylated regions (DMRs) is a basic task in DNA methylation analysis. However, due to the different strategies adopted, different DMR sets will be predicted on the same dataset, which poses a challenge in selecting a reliable and comprehensive DMR set for downstream analysis.

RESULTS

Here, we develop DMRIntTk, a toolkit for integrating DMR sets predicted by different methods on a same dataset. In DMRIntTk, the genome is segmented into bins, and the reliability of each DMR set at different methylation thresholds is evaluated. Then, the bins are weighted based on the covered DMR sets and integrated into final DMRs using a density peak clustering algorithm. To demonstrate the practicality of DMRIntTk, it was applied to different scenarios, including tissues with relatively large methylation differences, cancer tissues versus normal tissues with medium methylation differences, and disease tissues versus normal tissues with subtle methylation differences. Our results show that DMRIntTk can effectively trim regions with small methylation differences from the original DMR sets and thereby enriching the proportion of DMRs with larger methylation differences. In addition, the overlap analysis suggests that the integrated DMR sets are quite comprehensive, and functional analyses indicate the integrated disease-related DMRs are significantly enriched in biological pathways associated with the pathological mechanisms of the diseases. A comparative analysis of the integrated DMR set versus each original DMR set further highlights the superiority of DMRIntTk, demonstrating the unique biological insights it can provide.

CONCLUSIONS

Conclusively, DMRIntTk can help researchers obtain a reliable and comprehensive DMR set from many prediction methods.

Multi-Omic Graph Diagnosis (MOGDx): a data integration tool to perform classification tasks for heterogeneous diseases

MOTIVATION

Heterogeneity in human diseases presents challenges in diagnosis and treatments due to the broad range of manifestations and symptoms. With the rapid development of labelled multi-omic data, integrative machine learning methods have achieved breakthroughs in treatments by redefining these diseases at a more granular level. These approaches often have limitations in scalability, oversimplification, and handling of missing data.

RESULTS

In this study, we introduce Multi-Omic Graph Diagnosis (MOGDx), a flexible command line tool for the integration of multi-omic data to perform classification tasks for heterogeneous diseases. MOGDx has a network taxonomy. It fuses patient similarity networks, augments this integrated network with a reduced vector representation of genomic data and performs classification using a graph convolutional network. MOGDx was evaluated on three datasets from the cancer genome atlas for breast invasive carcinoma, kidney cancer, and low grade glioma. MOGDx demonstrated state-of-the-art performance and an ability to identify relevant multi-omic markers in each task. It integrated more genomic measures with greater patient coverage compared to other network integrative methods. Overall, MOGDx is a promising tool for integrating multi-omic data, classifying heterogeneous diseases, and aiding interpretation of genomic marker data.

AVAILABILITY AND IMPLEMENTATION

MOGDx source code is available from https://github.com/biomedicalinformaticsgroup/MOGDx.

Blood methylation pattern reflects epigenetic remodelling in adipose tissue after bariatric surgery

BACKGROUND

Studies on DNA methylation following bariatric surgery have primarily focused on blood cells, while it is unclear to which extend it may reflect DNA methylation profiles in specific metabolically relevant organs such as adipose tissue. Here, we investigated whether adipose tissue depots specific methylation changes after bariatric surgery are mirrored in blood.

METHODS

Using Illumina 850K EPIC technology, we analysed genome-wide DNA methylation in paired blood, subcutaneous and omental visceral AT (SAT/OVAT) samples from nine individuals (N = 6 female) with severe obesity pre- and post-surgery.

FINDINGS

The numbers and effect sizes of differentially methylated regions (DMRs) post-bariatric surgery were more pronounced in AT (SAT: 12,865 DMRs from -11.5 to 10.8%; OVAT: 14,632 DMRs from -13.7 to 12.8%) than in blood (9267 DMRs from -8.8 to 7.7%). Cross-tissue DMRs implicated immune-related genes. Among them, 49 regions could be validated with similar methylation changes in blood from independent individuals. Fourteen DMRs correlated with differentially expressed genes in AT post bariatric surgery, including downregulation of PIK3AP1 in both SAT and OVAT. DNA methylation age acceleration was significantly higher in AT compared to blood, but remained unaffected after surgery.

INTERPRETATION

Concurrent methylation pattern changes in blood and AT, particularly in immune-related genes, suggest blood DNA methylation mirrors AT's inflammatory state post-bariatric surgery.

FUNDING

The funding sources are listed in the Acknowledgments section.

Cortico-striatal differences in the epigenome in attention-deficit/ hyperactivity disorder

While epigenetic modifications have been implicated in ADHD through studies of peripheral tissue, to date there has been no examination of the epigenome of the brain in the disorder. To address this gap, we mapped the methylome of the caudate nucleus and anterior cingulate cortex in post-mortem tissue from fifty-eight individuals with or without ADHD. While no single probe showed adjusted significance in differential methylation, several differentially methylated regions emerged. These regions implicated genes involved in developmental processes including neurogenesis and the differentiation of oligodendrocytes and glial cells. We demonstrate a significant association between differentially methylated genes in the caudate and genes implicated by GWAS not only in ADHD but also in autistic spectrum, obsessive compulsive and bipolar affective disorders through GWAS. Using transcriptomic data available on the same subjects, we found modest correlations between the methylation and expression of genes. In conclusion, this study of the cortico-striatal methylome points to gene and gene pathways involved in neurodevelopment, consistent with studies of common and rare genetic variation, as well as the post-mortem transcriptome in ADHD.

A randomized placebo-controlled trial of nicotinamide riboside in older adults with mild cognitive impairment

Nicotinamide riboside (NR) increases blood levels of NAD+, a cofactor central to energy metabolism, and improves brain function in some rodent models of neurodegeneration. We conducted a placebo-controlled randomized pilot study with the primary objective of determining safety of NR in older adults with mild cognitive impairment (MCI). Twenty subjects with MCI were randomized to receive placebo or NR using dose escalation to achieve, and maintain, a final dose of 1 g/day over a 10-week study duration. The primary outcome was post-treatment change from baseline measures of cognition (Montreal Cognitive Assessment, MoCA). Predefined secondary outcomes included post-treatment changes in cerebral blood flow (CBF); blood NAD+ levels; and additional neurocognitive, psychometric, and physical performance tests. DNA methylation was assessed in peripheral blood mononuclear cells (PBMCs) as an exploratory outcome. The target NR dose was safely achieved as evidenced by a 2.6-fold increase in blood NAD+ in the NR group (p < 0.001, 95% CI [17.77, 43.49]) with no between-group difference in adverse event reporting. MoCA and other neurocognitive and psychometric metrics remained stable throughout the study. NR reduced CBF in the default mode network (DMN) with greatest differences observed in the left inferior parietal lobe (IPL) (DMN p = 0.013, μ = 0.92, 95% CI [0.23, 1.62]; left IPL p = 0.009, μ = 1.66, 95% CI [0.5, 2.82]). Walking speed in the placebo group significantly improved across the study duration suggestive of a practice effect but did not change in the NR group (p = 0.0402 and p = 0.4698, respectively). Other secondary outcome measures remained stable. Global methylation analyses indicated a modest NR-associated increase in DNA methylation and concomitant reduction in epigenetic age as measured by PhenoAge and GrimAge epigenetic clock analyses. In summary, NR significantly increased blood NAD+ concentrations in older adults with MCI. NR was well tolerated and did not alter cognition. While CBF was reduced by NR treatment, statistical significance would not have withstood multiple comparisons correction. A larger trial of longer duration is needed to determine the potential of NR as a strategy to improve cognition and alter CBF in older adults with MCI. ClinicalTrials.gov NCT02942888.

Alterations in leukocyte DNA methylome are associated to immunosuppression in severe clinical phenotypes of septic patients

Introduction

Sepsis patients experience a complex interplay of host pro- and anti-inflammatory processes which compromise the clinical outcome. Despite considering the latest clinical and scientific research, our comprehension of the immunosuppressive events in septic episodes remains incomplete. Additionally, a lack of data exists regarding the role of epigenetics in modulating immunosuppression, subsequently impacting patient survival.

Methods

To advance the current understanding of the mechanisms underlying immunosuppression, in this study we explored the dynamics of DNA methylation using the Infinium Methylation EPIC v1.0 BeadChip Kit in leukocytes from patients suffering from sepsis, septic shock, and critically ill patients as controls, within the first 24 h after admission in the Intensive Care Unit of a tertiary hospital.

Results and discussion

Employing two distinct analysis approaches (DMRcate and mCSEA) in comparing septic shock and critically ill patients, we identified 1,256 differentially methylated regions (DMRs) intricately linked to critical immune system pathways. The examination of the top 100 differentially methylated positions (DMPs) between septic shock and critically ill patients facilitated a clear demarcation among the three patient groups. Notably, the top 6,657 DMPs exhibited associations with organ dysfunction and lactate levels. Among the individual genes displaying significant differential methylation, IL10, TREM1, IL1B, and TNFAIP8 emerged with the most pronounced methylation alterations across the diverse patient groups when subjected to DNA bisulfite pyrosequencing analysis. These findings underscore the dynamic nature of DNA methylation profiles, highlighting the most pronounced alterations in patients with septic shock, and revealing their close association with the disease.

A novel principal component based method for identifying differentially methylated regions in Illumina Infinium MethylationEPIC BeadChip data

Differentially methylated regions (DMRs) are genomic regions with methylation patterns across multiple CpG sites that are associated with a phenotype. In this study, we proposed a Principal Component (PC) based DMR analysis method for use with data generated using the Illumina Infinium MethylationEPIC BeadChip (EPIC) array. We obtained methylation residuals by regressing the M-values of CpGs within a region on covariates, extracted PCs of the residuals, and then combined association information across PCs to obtain regional significance. Simulation-based genome-wide false positive (GFP) rates and true positive rates were estimated under a variety of conditions before determining the final version of our method, which we have named DMRPC. Then, DMRPC and another DMR method, coMethDMR, were used to perform epigenome-wide analyses of several phenotypes known to have multiple associated methylation loci (age, sex, and smoking) in a discovery and a replication cohort. Among regions that were analysed by both methods, DMRPC identified 50% more genome-wide significant age-associated DMRs than coMethDMR. The replication rate for the loci that were identified by only DMRPC was higher than the rate for those that were identified by only coMethDMR (90% for DMRPC vs. 76% for coMethDMR). Furthermore, DMRPC identified replicable associations in regions of moderate between-CpG correlation which are typically not analysed by coMethDMR. For the analyses of sex and smoking, the advantage of DMRPC was less clear. In conclusion, DMRPC is a new powerful DMR discovery tool that retains power in genomic regions with moderate correlation across CpGs.

The specific DNA methylation landscape in focal cortical dysplasia ILAE type 3D

Focal Cortical Dysplasia (FCD) is a frequent cause of drug-resistant focal epilepsy in children and young adults. The international FCD classifications of 2011 and 2022 have identified several clinico-pathological subtypes, either occurring isolated, i.e., FCD ILAE Type 1 or 2, or in association with a principal cortical lesion, i.e., FCD Type 3. Here, we addressed the DNA methylation signature of a previously described new subtype of FCD 3D occurring in the occipital lobe of very young children and microscopically defined by neuronal cell loss in cortical layer 4. We studied the DNA methylation profile using 850 K BeadChip arrays in a retrospective cohort of 104 patients with FCD 1 A, 2 A, 2B, 3D, TLE without FCD, and 16 postmortem specimens without neurological disorders as controls, operated in China or Germany. DNA was extracted from formalin-fixed paraffin-embedded tissue blocks with microscopically confirmed lesions, and DNA methylation profiles were bioinformatically analyzed with a recently developed deep learning algorithm. Our results revealed a distinct position of FCD 3D in the DNA methylation map of common FCD subtypes, also different from non-FCD epilepsy surgery controls or non-epileptic postmortem controls. Within the FCD 3D cohort, the DNA methylation signature separated three histopathology subtypes, i.e., glial scarring around porencephalic cysts, loss of layer 4, and Rasmussen encephalitis. Differential methylation in FCD 3D with loss of layer 4 mapped explicitly to biological pathways related to neurodegeneration, biogenesis of the extracellular matrix (ECM) components, axon guidance, and regulation of the actin cytoskeleton. Our data suggest that DNA methylation signatures in cortical malformations are not only of diagnostic value but also phenotypically relevant, providing the molecular underpinnings of structural and histopathological features associated with epilepsy. Further studies will be necessary to confirm these results and clarify their functional relevance and epileptogenic potential in these difficult-to-treat children.

Interactive DNA Methylation Array Analysis with ShinyÉPICo

Arrays provide a cost-effective platform for the analysis of human DNA methylation. ShinyÉPICo is an interactive, web-based, and graphical tool that allows the user to analyze Illumina DNA methylation arrays (450 k and EPIC), from the user's own computer or from a server. This tool covers the analysis entirely, from the raw data input to the final list of differentially methylated positions or regions. Here, we describe the steps of the analysis, the different parameters available, and useful information to understand and select the best options in each step.

Circulating cell-free DNA methylation mirrors alterations in cerebral patterns in epilepsy

BACKGROUND

DNA methylation profiling of circulating cell-free DNA (cfDNA) has rapidly become a promising strategy for biomarker identification and development. The cell-type-specific nature of DNA methylation patterns and the direct relationship between cfDNA and apoptosis can potentially be used non-invasively to predict local alterations. In addition, direct detection of altered DNA methylation patterns performs well as a biomarker. In a previous study, we demonstrated marked DNA methylation alterations in brain tissue from patients with mesial temporal lobe epilepsy with hippocampal sclerosis (MTLE-HS).

RESULTS

We performed DNA methylation profiling in cfDNA isolated from the serum of MTLE patients and healthy controls using BeadChip arrays followed by systematic bioinformatic analysis including deconvolution analysis and integration with DNase accessibility data sets. Differential cfDNA methylation analysis showed an overrepresentation of gene ontology terms and transcription factors related to central nervous system function and regulation. Deconvolution analysis of the DNA methylation data sets ruled out the possibility that the observed differences were due to changes in the proportional contribution of cortical neurons in cfDNA. Moreover, we found no overrepresentation of neuron- or glia-specific patterns in the described cfDNA methylation patterns. However, the MTLE-HS cfDNA methylation patterns featured a significant overrepresentation of the epileptic DNA methylation alterations previously observed in the hippocampus.

CONCLUSIONS

Our results support the use of cfDNA methylation profiling as a rational approach to seeking non-invasive and reproducible epilepsy biomarkers.

An evaluation of the genome-wide false positive rates of common methods for identifying differentially methylated regions using illumina methylation arrays

Differentially methylated regions (DMRs) are genomic regions with specific methylation patterns across multiple loci that are associated with a phenotype. We examined the genome-wide false positive (GFP) rates of five widely used DMR methods: comb-p, Bumphunter, DMRcate, mCSEA and coMethDMR using both Illumina HumanMethylation450 (450 K) and MethylationEPIC (EPIC) data and simulated continuous and dichotomous null phenotypes (i.e., generated independently of methylation data). coMethDMR provided well-controlled GFP rates (~5%) except when analysing skewed continuous phenotypes. DMRcate generally had well-controlled GFP rates when applied to 450 K data except for the skewed continuous phenotype and EPIC data only for the normally distributed continuous phenotype. GFP rates for mCSEA were at least 0.096 and comb-p yielded GFP rates above 0.34. Bumphunter had high GFP rates of at least 0.35 across conditions, reaching as high as 0.95. Analysis of the performance of these methods in specific regions of the genome found that regions with higher correlation across loci had higher regional false positive rates on average across methods. Based on the false positive rates, coMethDMR is the most recommended analysis method, and DMRcate had acceptable performance when analysing 450 K data. However, as both could display higher levels of FPs for skewed continuous distributions, a normalizing transformation of skewed continuous phenotypes is suggested. This study highlights the importance of genome-wide simulations when evaluating the performance of DMR-analysis methods.

GPU-accelerated and pipelined methylation calling

Motivation

The third-generation DNA sequencing technologies, such as Nanopore Sequencing, can operate at very high speeds and produce longer reads, which in turn results in a challenge for the computational analysis of such massive data. Nanopolish is a software package for signal-level analysis of Oxford Nanopore sequencing data. Call-methylation module of Nanopolish can detect methylation based on Hidden Markov Model (HMM). However, Nanopolish is limited by the long running time of some serial and computationally expensive processes. Among these, Adaptive Banded Event Alignment (ABEA) is the most time-consuming step, and the prior work, f5c, has already parallelized and optimized ABEA on GPU. As a result, the remaining methylation score calculation part, which uses HMM to identify if a given base is methylated or not, has become the new performance bottleneck.

Results

This article focuses on the call-methylation module that resides in the Nanopolish package. We propose Galaxy-methyl, which parallelizes and optimizes the methylation score calculation step on GPU and then pipelines the four steps of the call-methylation module. Galaxy-methyl increases the execution concurrency across CPUs and GPUs as well as hardware resource utilization for both. The experimental results collected indicate that Galaxy-methyl can achieve 3×-5× speedup compared with Nanopolish, and reduce the total execution time by 35% compared with f5c, on average.

Availability and implementation

The source code of Galaxy-methyl is available at https://github.com/fengyilin118/.

Risperidone-induced changes in DNA methylation in peripheral blood from first-episode schizophrenia patients parallel changes in neuroimaging and cognitive phenotypes

BACKGROUND

Second generation antipsychotics such as risperidone are first-line pharmacotherapy treatment choices for schizophrenia. However, our ability to reliably predict and monitor treatment reaction is impeded by the lack of relevant biomarkers. As a biomarker for the susceptibility of schizophrenia and clozapine treatment response, DNA methylation (DNAm) has been studied, but the impact of antipsychotics on DNAm has not been explored in drug-naïve patients.

OBJECTIVE

The aim of the present study was to examine changes of DNAm after short-term antipsychotic therapy in first-episode drug-naïve schizophrenia (FES) to identify the beneficial and adverse effect of risperidone on DNAm and their relation to treatment outcome.

METHODS

Thirty-eight never treated schizophrenia patients and 38 demographically matched individuals (healthy controls) were assessed at baseline and at 8-week follow-up with symptom ratings, and cognitive and functional imaging procedures, at which time a blood draw for DNAm studies was performed. During the 8-week period, patients received treatment with risperidone monotherapy. An independent data set was used as replication.

RESULTS

We identified brain related pathways enriched in 4,888 FES-associated CpG sites relative to controls. Risperidone administration in patients altered DNAm in 5,979 CpG sites relative to baseline. Significant group differences in DNAm at follow-up were seen in FES patients at 6,760 CpG sites versus healthy controls. Through comparison of effect size, we found 87.54% out of the risperidone-associated changes in DNAm showed possible beneficial effect, while only 12.46% showed potential adverse effect. There were 580 DNAm sites in which changes shifted methylation levels to be indistinguishable from controls after risperidone treatment. The DNAm changes of some sites that normalized after treatment were correlated with treatment-related changes in symptom severity, spontaneous neurophysiological activity, and cognition. We replicated our results in an independent data set.

CONCLUSION

The normalizing effect of risperidone monotherapy on gene DNAm, and its correlation with clinically relevant phenotypes, indicates that risperidone therapy is associated with DNAm changes that are related to changes in brain physiology, cognition and symptom severity.

Lineage-coupled clonal capture identifies clonal evolution mechanisms and vulnerabilities of BRAFV600E inhibition resistance in melanoma

Targeted cancer therapies have revolutionized treatment but their efficacies are limited by the development of resistance driven by clonal evolution within tumors. We developed "CAPTURE", a single-cell barcoding approach to comprehensively trace clonal dynamics and capture live lineage-coupled resistant cells for in-depth multi-omics analysis and functional exploration. We demonstrate that heterogeneous clones, either preexisting or emerging from drug-tolerant persister cells, dominated resistance to vemurafenib in BRAFV600E melanoma. Further integrative studies uncovered diverse resistance mechanisms. This includes a previously unrecognized and clinically relevant mechanism, chromosome 18q21 gain, which leads to vulnerability of the cells to BCL2 inhibitor. We also identified targetable common dependencies of captured resistant clones, such as oxidative phosphorylation and E2F pathways. Our study provides new therapeutic insights into overcoming therapy resistance in BRAFV600E melanoma and presents a platform for exploring clonal evolution dynamics and vulnerabilities that can be applied to study treatment resistance in other cancers.

Pioneer transcription factors are associated with the modulation of DNA methylation patterns across cancers

Methylation of cytosines on DNA is a prominent modification associated with gene expression regulation. Aberrant DNA methylation patterns have recurrently been linked to dysregulation of the regulatory program in cancer cells. To shed light on the underlying molecular mechanism driving this process, we hypothesised that aberrant methylation patterns could be controlled by the binding of specific transcription factors (TFs) across cancer types. By combining DNA methylation arrays and gene expression data with TF binding sites (TFBSs), we explored the interplay between TF binding and DNA methylation in 19 cancer types. We performed emQTL (expression-methylation quantitative trait loci) analyses independently in each cancer type and identified 13 TFs whose expression levels are correlated with local DNA methylation patterns around their binding sites in at least 2 cancer types. The 13 TFs are mainly associated with local demethylation and are enriched for pioneer function, suggesting a specific role for these TFs in modulating chromatin structure and transcription in cancer patients. Furthermore, we confirmed that de novo methylation is precluded across cancers at CpGs lying in genomic regions enriched for TF binding signatures associated with SP1, CTCF, NRF1, GABPA, KLF9, and/or YY1. The modulation of DNA methylation associated with TF binding was observed at cis-regulatory regions controlling immune- and cancer-associated pathways, corroborating that the emQTL signals were derived from both cancer and tumor-infiltrating cells. As a case example, we experimentally confirmed that FOXA1 knock-down is associated with higher methylation in regions bound by FOXA1 in breast cancer MCF-7 cells. Finally, we reported physical interactions between FOXA1 with TET1 and TET2 both in an in vitro setup and in vivo at physiological levels in MCF-7 cells, adding further support for FOXA1 attracting TET1 and TET2 to induce local demethylation in cancer cells.

Astroblastomas exhibit radial glia stem cell lineages and differential expression of imprinted and X-inactivation escape genes

Astroblastomas (ABs) are rare brain tumors of unknown origin. We performed an integrative genetic and epigenetic analysis of AB-like tumors. Here, we show that tumors traceable to neural stem/progenitor cells (radial glia) that emerge during early to later brain development occur in children and young adults, respectively. Tumors with MN1-BEND2 fusion appear to present exclusively in females and exhibit overexpression of genes expressed prior to 25 post-conception weeks (pcw), including genes enriched in early ventricular zone radial glia and ependymal tumors. Other, histologically classic ABs overexpress or harbor mutations of mitogen-activated protein kinase pathway genes, outer and truncated radial glia genes, and genes expressed after 25 pcw, including neuronal and astrocyte markers. Findings support that AB-like tumors arise in the context of epigenetic and genetic changes in neural progenitors. Selective gene fusion, variable imprinting and/or chromosome X-inactivation escape resulting in biallelic overexpression may contribute to female predominance of AB molecular subtypes.

Leukocyte methylomic imprints of exposure to the genocide against the Tutsi in Rwanda: a pilot epigenome-wide analysis

Aim & methods: We conducted a pilot epigenome-wide association study of women from Tutsi ethnicity exposed to the genocide while pregnant and their resulting offspring, and a comparison group of women who were pregnant at the time of the genocide but living outside of Rwanda.Results: Fifty-nine leukocyte-derived DNA samples survived quality control: 33 mothers (20 exposed, 13 unexposed) and 26 offspring (16 exposed, 10 unexposed). Twenty-four significant differentially methylated regions (DMRs) were identified in mothers and 16 in children. Conclusions:In utero genocide exposure was associated with CpGs in three of the 24 DMRs: BCOR, PRDM8 and VWDE, with higher DNA methylation in exposed versus unexposed offspring. Of note, BCOR and VWDE show significant correlation between brain and blood DNA methylation within individuals, suggesting these peripherally derived signals of genocide exposure may have relevance to the brain.

Coordinated glucocorticoid receptor and MAFB action induces tolerogenesis and epigenome remodeling in dendritic cells

Glucocorticoids (GCs) exert potent anti-inflammatory effects in immune cells through the glucocorticoid receptor (GR). Dendritic cells (DCs), central actors for coordinating immune responses, acquire tolerogenic properties in response to GCs. Tolerogenic DCs (tolDCs) have emerged as a potential treatment for various inflammatory diseases. To date, the underlying cell type-specific regulatory mechanisms orchestrating GC-mediated acquisition of immunosuppressive properties remain poorly understood. In this study, we investigated the transcriptomic and epigenomic remodeling associated with differentiation to DCs in the presence of GCs. Our analysis demonstrates a major role of MAFB in this process, in synergy with GR. GR and MAFB both interact with methylcytosine dioxygenase TET2 and bind to genomic loci that undergo specific demethylation in tolDCs. We also show that the role of MAFB is more extensive, binding to thousands of genomic loci in tolDCs. Finally, MAFB knockdown erases the tolerogenic properties of tolDCs and reverts the specific DNA demethylation and gene upregulation. The preeminent role of MAFB is also demonstrated in vivo for myeloid cells from synovium in rheumatoid arthritis following GC treatment. Our results imply that, once directly activated by GR, MAFB plays a critical role in orchestrating the epigenomic and transcriptomic remodeling that define the tolerogenic phenotype.

Integrative epigenomics in Sjögren´s syndrome reveals novel pathways and a strong interaction between the HLA, autoantibodies and the interferon signature

Primary Sjögren's syndrome (SS) is a systemic autoimmune disease characterized by lymphocytic infiltration and damage of exocrine salivary and lacrimal glands. The etiology of SS is complex with environmental triggers and genetic factors involved. By conducting an integrated multi-omics study, we confirmed a vast coordinated hypomethylation and overexpression effects in IFN-related genes, what is known as the IFN signature. Stratified and conditional analyses suggest a strong interaction between SS-associated HLA genetic variation and the presence of Anti-Ro/SSA autoantibodies in driving the IFN epigenetic signature and determining SS. We report a novel epigenetic signature characterized by increased DNA methylation levels in a large number of genes enriched in pathways such as collagen metabolism and extracellular matrix organization. We identified potential new genetic variants associated with SS that might mediate their risk by altering DNA methylation or gene expression patterns, as well as disease-interacting genetic variants that exhibit regulatory function only in the SS population. Our study sheds new light on the interaction between genetics, autoantibody profiles, DNA methylation and gene expression in SS, and contributes to elucidate the genetic architecture of gene regulation in an autoimmune population.

Severity of Idiopathic Scoliosis Is Associated with Differential Methylation: An Epigenome-Wide Association Study of Monozygotic Twins with Idiopathic Scoliosis

Epigenetic mechanisms may contribute to idiopathic scoliosis (IS). We identified 8 monozygotic twin pairs with IS, 6 discordant (Cobb angle difference > 10°) and 2 concordant (Cobb angle difference ≤ 2°). Genome-wide methylation in blood was measured with the Infinium HumanMethylation EPIC Beadchip. We tested for differences in methylation and methylation variability between discordant twins and tested the association between methylation and curve severity in all twins. Differentially methylated region (DMR) analyses identified gene promoter regions. Methylation at cg12959265 (chr. 7 DPY19L1) was less variable in cases (false discovery rate (FDR) = 0.0791). We identified four probes (false discovery rate, FDR < 0.10); cg02477677 (chr. 17, RARA gene), cg12922161 (chr. 2 LOC150622 gene), cg08826461 (chr. 2), and cg16382077 (chr. 7) associated with curve severity. We identified 57 DMRs where hyper- or hypo-methylation was consistent across the region and 28 DMRs with a consistent association with curve severity. Among DMRs, 21 were correlated with bone methylation. Prioritization of regions based on methylation concordance in bone identified promoter regions for WNT10A (WNT signaling), NPY (regulator of bone and energy homeostasis), and others predicted to be relevant for bone formation/remodeling. These regions may aid in understanding the complex interplay between genetics, environment, and IS.

Gene set enrichment analysis for genome-wide DNA methylation data

DNA methylation is one of the most commonly studied epigenetic marks, due to its role in disease and development. Illumina methylation arrays have been extensively used to measure methylation across the human genome. Methylation array analysis has primarily focused on preprocessing, normalization, and identification of differentially methylated CpGs and regions. GOmeth and GOregion are new methods for performing unbiased gene set testing following differential methylation analysis. Benchmarking analyses demonstrate GOmeth outperforms other approaches, and GOregion is the first method for gene set testing of differentially methylated regions. Both methods are publicly available in the missMethyl Bioconductor R package.

Identification of two methylated fragments of an SDC2 CpG island using a sliding window technique for early detection of colorectal cancer

Colorectal cancer (CRC) is one of the most common cancer types globally with a 5‐year survival rate of < 50% in China. Aberrant DNA methylation is one of the hallmarks of tumor initiation, progression, and metastasis. Here, we investigated the clinical performance of two differentially methylated regions (DMRs) in SDC2 CpG islands for the detection of CRC. A sliding window technique was used to identify the DMRs, and methylation‐specific PCR assay was used to assess the DMRs in 198 CRC samples and 54 normal controls. Two DMRs (DMR2 and DMR5) were identified using The Cancer Genome Atlas (TCGA) data, and the hypermethylation of DMR2 and DMR5 was detected in 90.91% (180/198) and 89.90% (178/198) of CRC samples, respectively. When combining DMR2 and DMR5, the sensitivity for CRC detection was 94.4% higher than that of DMR2 or DMR5 alone. Based on the above results, we propose using DMR2 and DMR5 as a sensitive biomarker to detect CRC.

Detecting Differentially Methylated Promoters in Genes Related to Disease Phenotypes Using R

DNA methylation in gene promoters plays a major role in gene expression regulation, and alterations in methylation patterns have been associated with several diseases. In this context, different software suites and statistical methods have been proposed to analyze differentially methylated positions and regions. Among them, the novel statistical method implemented in the mCSEA R package proposed a new framework to detect subtle, but consistent, methylation differences. Here, we provide an easy-to-use pipeline covering all the necessary steps to detect differentially methylated promoters with mCSEA from Illumina 450K and EPIC methylation BeadChips data. This protocol covers the download of data from public repositories, quality control, data filtering and normalization, estimation of cell type proportions, and statistical analysis. In addition, we show the procedure to compare disease vs. normal phenotypes, obtaining differentially methylated regions including promoters or CpG Islands. The entire protocol is based on R programming language, which can be used in any operating system and does not require advanced programming skills.

Wavelet Screening identifies regions highly enriched for differentially methylated loci for orofacial clefts

DNA methylation is the most widely studied epigenetic mark in humans and plays an essential role in normal biological processes as well as in disease development. More focus has recently been placed on understanding functional aspects of methylation, prompting the development of methods to investigate the relationship between heterogeneity in methylation patterns and disease risk. However, most of these methods are limited in that they use simplified models that may rely on arbitrarily chosen parameters, they can only detect differentially methylated regions (DMRs) one at a time, or they are computationally intensive. To address these shortcomings, we present a wavelet-based method called 'Wavelet Screening' (WS) that can perform an epigenome-wide association study (EWAS) of thousands of individuals on a single CPU in only a matter of hours. By detecting multiple DMRs located near each other, WS identifies more complex patterns that can differentiate between different methylation profiles. We performed an extensive set of simulations to demonstrate the robustness and high power of WS, before applying it to a previously published EWAS dataset of orofacial clefts (OFCs). WS identified 82 associated regions containing several known genes and loci for OFCs, while other findings are novel and warrant replication in other OFCs cohorts.

DNA Methylation Analysis to Unravel Altered Genetic Pathways Underlying Early Onset and Late Onset Neonatal Sepsis. A Pilot Study

Background: Neonatal sepsis is a systemic condition widely affecting preterm infants and characterized by pro-inflammatory and anti-inflammatory responses. However, its pathophysiology is not yet fully understood. Epigenetics regulates the immune system, and its alteration leads to the impaired immune response underlying sepsis. DNA methylation may contribute to sepsis-induced immunosuppression which, if persistent, will cause long-term adverse effects in neonates.

Objective: To analyze the methylome of preterm infants in order to determine whether there are DNA methylation marks that may shed light on the pathophysiology of neonatal sepsis.

Design: Prospective observational cohort study performed in the neonatal intensive care unit (NICU) of a tertiary care center.

Patients: Eligible infants were premature ≤32 weeks admitted to the NICU with clinical suspicion of sepsis. The methylome analysis was performed in DNA from blood using Infinium Human Methylation EPIC microarrays to uncover methylation marks.

Results: Methylation differential analysis revealed an alteration of methylation levels in genomic regions involved in inflammatory pathways which participate in both the innate and the adaptive immune response. Moreover, differences between early and late onset sepsis as compared to normal controls were assessed.

Conclusions: DNA methylation marks can serve as a biomarker for neonatal sepsis and even contribute to differentiating between early and late onset sepsis.

shinyÉPICo: A graphical pipeline to analyze Illumina DNA methylation arrays

SUMMARY

Illumina DNA methylation bead arrays provide a cost-effective platform for the simultaneous analysis of a high number of human samples. However, the analysis can be time-demanding and requires some computational expertise. shinyÉPICo is an interactive, web-based, and graphical tool that allows the user to analyze Illumina DNA methylation arrays (450k and EPIC), from the user's own computer or from a server. The tool covers the entire analysis, from the raw data to the final list of differentially methylated positions and differentially methylated regions between sample groups. It allows the user to test several normalization methods, linear model parameters, including covariates, and differentially methylated CpGs filters, in a quick and easy manner, with interactive graphics helping to select the options in each step. shinyÉPICo represents a comprehensive tool for standardizing and accelerating DNA methylation analysis, as well as optimizing computational resources in laboratories studying DNA methylation.

AVAILABILITY AND IMPLEMENTATION

shinyÉPICo is freely available as an R package at the Bioconductor project (http://bioconductor.org/packages/shinyepico/) and GitHub (https://github.com/omorante/shinyepico) under an AGPL3 license.

PTSD is associated with increased DNA methylation across regions of HLA-DPB1 and SPATC1L

Posttraumatic stress disorder (PTSD) is characterized by intrusive thoughts, avoidance, negative alterations in cognitions and mood, and arousal symptoms that adversely affect mental and physical health. Recent evidence links changes in DNA methylation of CpG cites to PTSD. Since clusters of proximal CpGs share similar methylation signatures, identification of PTSD-associated differentially methylated regions (DMRs) may elucidate the pathways defining differential risk and resilience of PTSD. Here we aimed to identify epigenetic differences associated with PTSD. DNA methylation data profiled from blood samples using the MethylationEPIC BeadChip were used to perform a DMR analysis in 187 PTSD cases and 367 trauma-exposed controls from the Grady Trauma Project (GTP). DMRs were assessed with R package bumphunter. We identified two regions that associate with PTSD after multiple test correction. These regions were in the gene body of HLA-DPB1 and in the promoter of SPATC1L. The DMR in HLA-DPB1 was associated with PTSD in an independent cohort. Both DMRs included CpGs whose methylation associated with nearby sequence variation (meQTL) and that associated with expression of their respective genes (eQTM). This study supports an emerging literature linking PTSD risk to genetic and epigenetic variation in the HLA region.

DECO: decompose heterogeneous population cohorts for patient stratification and discovery of sample biomarkers using omic data profiling

Motivation

Patient and sample diversity is one of the main challenges when dealing with clinical cohorts in biomedical genomics studies. During last decade, several methods have been developed to identify biomarkers assigned to specific individuals or subtypes of samples. However, current methods still fail to discover markers in complex scenarios where heterogeneity or hidden phenotypical factors are present. Here, we propose a method to analyze and understand heterogeneous data avoiding classical normalization approaches of reducing or removing variation.

Results

DEcomposing heterogeneous Cohorts using Omic data profiling (DECO) is a method to find significant association among biological features (biomarkers) and samples (individuals) analyzing large-scale omic data. The method identifies and categorizes biomarkers of specific phenotypic conditions based on a recurrent differential analysis integrated with a non-symmetrical correspondence analysis. DECO integrates both omic data dispersion and predictor–response relationship from non-symmetrical correspondence analysis in a unique statistic (called h-statistic), allowing the identification of closely related sample categories within complex cohorts. The performance is demonstrated using simulated data and five experimental transcriptomic datasets, and comparing to seven other methods. We show DECO greatly enhances the discovery and subtle identification of biomarkers, making it especially suited for deep and accurate patient stratification.

Availability and implementation

DECO is freely available as an R package (including a practical vignette) at Bioconductor repository (http://bioconductor.org/packages/deco/).

Supplementary information

Supplementary data are available at Bioinformatics online.

Longitudinal epigenome-wide association studies of three male military cohorts reveal multiple CpG sites associated with post-traumatic stress disorder

BACKGROUND

Epigenetic mechanisms have been suggested to play a role in the development of post-traumatic stress disorder (PTSD). Here, blood-derived DNA methylation data (HumanMethylation450 BeadChip) collected prior to and following combat exposure in three cohorts of male military members were analyzed to assess whether DNA methylation profiles are associated with the development of PTSD. A total of 123 PTSD cases and 143 trauma-exposed controls were included in the analyses. The Psychiatric Genomics Consortium (PGC) PTSD EWAS QC pipeline was used on all cohorts, and results were combined using a sample size weighted meta-analysis in a two-stage design. In stage one, we jointly analyzed data of two new cohorts (N = 126 and 78) for gene discovery, and sought to replicate significant findings in a third, previously published cohort (N = 62) to assess the robustness of our results. In stage 2, we aimed at maximizing power for gene discovery by combining all three cohorts in a meta-analysis.

RESULTS

Stage 1 analyses identified four CpG sites in which, conditional on pre-deployment DNA methylation, post-deployment DNA methylation was significantly associated with PTSD status after epigenome-wide adjustment for multiple comparisons. The most significant (intergenic) CpG cg05656210 (p = 1.0 × 10-08) was located on 5q31 and significantly replicated in the third cohort. In addition, 19 differentially methylated regions (DMRs) were identified, but failed replication. Stage 2 analyses identified three epigenome-wide significant CpGs, the intergenic CpG cg05656210 and two additional CpGs located in MAD1L1 (cg12169700) and HEXDC (cg20756026). Interestingly, cg12169700 had an underlying single nucleotide polymorphism (SNP) which was located within the same LD block as a recently identified PTSD-associated SNP in MAD1L1. Stage 2 analyses further identified 12 significant differential methylated regions (DMRs), 1 of which was located in MAD1L1 and 4 were situated in the human leukocyte antigen (HLA) region.

CONCLUSIONS

This study suggests that the development of combat-related PTSD is associated with distinct methylation patterns in several genomic positions and regions. Our most prominent findings suggest the involvement of the immune system through the HLA region and HEXDC, and MAD1L1 which was previously associated with PTSD.