DNA Methylation and Transcriptomic Features are Preserved Throughout Disease Recurrence and Chemoresistance in High Grade Serous Ovarian Cancers

Background Little is known about the role of global DNA methylation in recurrence and chemoresistance of high grade serous ovarian cancer (HGSOC). We performed whole genome bisulfite sequencing and transcriptome sequencing in 62 primary and recurrent tumors from 28 patients with stage III/IV HGSOC, of which 11 patients carried germline, pathogenic BRCA1 and/or BRCA2 mutations. Results Landscapes of genome-wide methylation (on average 24.2 million CpGs per tumor) and transcriptomes in primary and recurrent tumors showed extensive heterogeneity between patients but were highly preserved in tumors from the same patient. We identified significant differences in the burden of differentially methylated regions (DMRs) in tumors from BRCA1/2 compared to non-BRCA1/2 carriers (mean 659 DMRs and 388 DMRs in paired comparisons respectively). We identified overexpression of immune pathways in BRCA1/2 carriers compared to non-carriers, implicating an increased immune response in improved survival (P=0.006) in these BRCA1/2 carriers. Conclusions These findings indicate methylome and gene expression programs established in the primary tumor are conserved throughout disease progression, even extensive chemotherapy treatment, and that changes in methylation and gene expression are unlikely to serve as drivers for chemoresistance in HGSOC.

Single-Base Resolution Methylomes of Somatic Embryogenesis in Theobroma Cacao L. Reveal Epigenome Modifications Associated With Somatic Embryo Abnormalities

Propagation by somatic embryogenesis in Theobroma cacao has some issues to be solved, as many morphologically abnormal somatic embryos that do not germinate into plants are frequently observed, thus hampering plant production on a commercial scale. For the first time the methylome landscape of T. cacao somatic embryogenesis was examined, using whole-genome bisulfite sequencing technique, with the aim to understand the epigenetic basis of somatic embryo abnormalities. We identified 873 differentially methylated genes (DMGs) in the CpG context between zygotic embryos, normal and abnormal somatic embryos, with important roles in development, programmed cell death, oxidative stress, and hypoxia induction, which can help to explain the morphological abnormalities of somatic embryos. We also identified the role of ethylene and its precursor 1-aminocyclopropane-1-carboxylate in several biological processes, such as hypoxia induction, cell differentiation and cell polarity, that could be associated to the development of abnormal somatic embryos. The biological processes and the hypothesis of ethylene and its precursor involvement in the somatic embryo abnormalities in cacao are discussed.

Characterization of epigenetic alterations in esophageal cancer by whole-genome bisulfite sequencing

Esophageal carcinoma is a common and aggressive malignancy, and its patients have dismal clinical outcomes. The epigenetic dysregulation in both major subtypes, esophageal squamous cell carcinoma (ESCC) and adenocarcinoma (EAC), awaits further characterization. Here, we perform whole-genome bisulfite sequencing (WGBS) on a total of 43 esophageal cancer and normal samples, generating one of the largest WGBS datasets in this cancer to date. Focusing on hypomethylated regions in cancer, we show that they are associated with increased chromatin activity and enhancer RNA expression. Using this large collection of WGBS dataset, we reveal and validate novel clusters in both ESCC and EAC. We further identify specific molecular features in each cluster, with potential clinical implications. These data together advance our understanding of the epigenetic alterations in esophageal cancer and provide a rich resource for the research community of this disease.

Focal disruption of DNA methylation dynamics at enhancers in IDH-mutant AML cells

Recurrent mutations in IDH1 or IDH2 in acute myeloid leukemia (AML) are associated with increased DNA methylation, but the genome-wide patterns of this hypermethylation phenotype have not been comprehensively studied in AML samples. We analyzed whole-genome bisulfite sequencing data from 15 primary AML samples with IDH1 or IDH2 mutations, which identified ~4000 focal regions that were uniquely hypermethylated in IDHmut samples vs. normal CD34+ cells and other AMLs. These regions had modest hypermethylation in AMLs with biallelic TET2 mutations, and levels of 5-hydroxymethylation that were diminished in IDH and TET-mutant samples, indicating that this hypermethylation results from inhibition of TET-mediated demethylation. Focal hypermethylation in IDHmut AMLs occurred at regions with low methylation in CD34+ cells, implying that DNA methylation and demethylation are active at these loci. AML samples containing IDH and DNMT3AR882 mutations were significantly less hypermethylated, suggesting that IDHmut-associated hypermethylation is mediated by DNMT3A. IDHmut-specific hypermethylation was highly enriched for enhancers that form direct interactions with genes involved in normal hematopoiesis and AML, including MYC and ETV6. These results suggest that focal hypermethylation in IDH-mutant AML occurs by altering the balance between DNA methylation and demethylation, and that disruption of these pathways at enhancers may contribute to AML pathogenesis.

The SEQC2 epigenomics quality control (EpiQC) study

Background Cytosine modifications in DNA such as 5-methylcytosine (5mC) underlie a broad range of developmental processes, maintain cellular lineage specification, and can define or stratify types of cancer and other diseases. However, the wide variety of approaches available to interrogate these modifications has created a need for harmonized materials, methods, and rigorous benchmarking to improve genome-wide methylome sequencing applications in clinical and basic research. Here, we present a multi-platform assessment and cross-validated resource for epigenetics research from the FDA’s Epigenomics Quality Control Group. Results Each sample is processed in multiple replicates by three whole-genome bisulfite sequencing (WGBS) protocols (TruSeq DNA methylation, Accel-NGS MethylSeq, and SPLAT), oxidative bisulfite sequencing (TrueMethyl), enzymatic deamination method (EMSeq), targeted methylation sequencing (Illumina Methyl Capture EPIC), single-molecule long-read nanopore sequencing from Oxford Nanopore Technologies, and 850k Illumina methylation arrays. After rigorous quality assessment and comparison to Illumina EPIC methylation microarrays and testing on a range of algorithms (Bismark, BitmapperBS, bwa-meth, and BitMapperBS), we find overall high concordance between assays, but also differences in efficiency of read mapping, CpG capture, coverage, and platform performance, and variable performance across 26 microarray normalization algorithms. Conclusions The data provided herein can guide the use of these DNA reference materials in epigenomics research, as well as provide best practices for experimental design in future studies. By leveraging seven human cell lines that are designated as publicly available reference materials, these data can be used as a baseline to advance epigenomics research.

Clinical response to nivolumab in an INI1-deficient pediatric chordoma correlates with immunogenic recognition of brachyury

Poorly differentiated chordoma (PDC) is a recently recognized subtype of chordoma characterized by expression of the embryonic transcription factor, brachyury, and loss of INI1. PDC primarily affects children and is associated with a poor prognosis and limited treatment options. Here we describe the molecular and immune tumour microenvironment profiles of two paediatric PDCs produced using whole-genome, transcriptome and whole-genome bisulfite sequencing (WGBS) and multiplex immunohistochemistry. Our analyses revealed the presence of tumour-associated immune cells, including CD8+ T cells, and expression of the immune checkpoint protein, PD-L1, in both patient samples. Molecular profiling provided the rationale for immune checkpoint inhibitor (ICI) therapy, which resulted in a clinical and radiographic response. A dominant T cell receptor (TCR) clone specific for a brachyury peptide–MHC complex was identified from bulk RNA sequencing, suggesting that targeting of the brachyury tumour antigen by tumour-associated T cells may underlie this clinical response to ICI. Correlative analysis with rhabdoid tumours, another INI1-deficient paediatric malignancy, suggests that a subset of tumours may share common immune phenotypes, indicating the potential for a therapeutically targetable subgroup of challenging paediatric cancers.

The DNMT1 inhibitor GSK-3484862 mediates global demethylation in murine embryonic stem cells

Background DNA methylation plays an important role in regulating gene expression in mammals. The covalent DNMT1 inhibitors 5-azacytidine and decitabine are widely used in research to reduce DNA methylation levels, but they impart severe cytotoxicity which limits their demethylation capability and confounds interpretation of experiments. Recently, a non-covalent inhibitor of DNMT1 called GSK-3484862 was developed by GlaxoSmithKline. We sought to determine whether GSK-3484862 can induce demethylation more effectively than 5-azanucleosides. Murine embryonic stem cells (mESCs) are an ideal cell type in which to conduct such experiments, as they have a high degree of DNA methylation but tolerate dramatic methylation loss. Results We determined the cytotoxicity and optimal concentration of GSK-3484862 by treating wild-type (WT) or Dnmt1/3a/3b triple knockout (TKO) mESC with different concentrations of the compound, which was obtained from two commercial sources. Concentrations of 10 µM or below were readily tolerated for 14 days of culture. Known DNA methylation targets such as germline genes and GLN-family transposons were upregulated within 2 days of the start of GSK-3484862 treatment. By contrast, 5-azacytidine and decitabine induced weaker upregulation of methylated genes and extensive cell death. Whole-genome bisulfite sequencing showed that treatment with GSK-3484862 induced dramatic DNA methylation loss, with global CpG methylation levels falling from near 70% in WT mESC to less than 18% after 6 days of treatment with GSK-3484862. The treated cells showed a methylation level and pattern similar to that observed in Dnmt1-deficient mESCs. Conclusions GSK-3484862 mediates striking demethylation in mESCs with minimal non-specific toxicity.

Insight Between the Epigenetics and Transcription Responding of Cotton Hypocotyl Cellular Elongation Under Salt-Alkaline Stress

Gossypium barbadense is a cultivated cotton not only known for producing superior fiber but also for its salt and alkaline resistance. Here, we used Whole Genome Bisulfite Sequencing (WGBS) technology to map the cytosine methylation of the whole genome of the G. barbadense hypocotyl at single base resolution. The methylation sequencing results showed that the mapping rates of the three samples were 75.32, 77.54, and 77.94%, respectively. In addition, the Bisulfite Sequence (BS) conversion rate was 99.78%. Approximately 71.03, 53.87, and 6.26% of the cytosine were methylated at CG, CHG, and CHH sequence contexts, respectively. A comprehensive analysis of DNA methylation and transcriptome data showed that the methylation level of the promoter region was a positive correlation in the CHH context. Saline-alkaline stress was related to the methylation changes of many genes, transcription factors (TFs) and transposable elements (TEs), respectively. We explored the regulatory mechanism of DNA methylation in response to salt and alkaline stress during cotton hypocotyl elongation. Our data shed light into the relationship of methylation regulation at the germination stage of G. barbadense hypocotyl cell elongation and salt-alkali treatment. The results of this research help understand the early growth regulation mechanism of G. barbadense in response to abiotic stress.

Integrative Analysis of Epigenome and Transcriptome Data Reveals Aberrantly Methylated Promoters and Enhancers in Hepatocellular Carcinoma

DNA methylation is a key transcription regulator, whose aberration was ubiquitous and important in most cancers including hepatocellular carcinoma (HCC). Whole-genome bisulfite sequencing (WGBS) was conducted for comparison of DNA methylation in tumor and adjacent tissues from 33 HCC patients, accompanying RNA-seq to determine differentially methylated region-associated, differentially expressed genes (DMR-DEGs), which were independently replicated in the TCGA-LIHC cohort and experimentally validated via 5-aza-2-deoxycytidine (5-azadC) demethylation. A total of 9,867,700 CpG sites showed significantly differential methylation in HCC. Integrations of mRNA-seq, histone ChIP-seq, and WGBS data identified 611 high-confidence DMR-DEGs. Enrichment analysis demonstrated activation of multiple molecular pathways related to cell cycle and DNA repair, accompanying repression of several critical metabolism pathways such as tyrosine and monocarboxylic acid metabolism. In TCGA-LIHC, we replicated about 53% of identified DMR-DEGs and highlighted the prognostic significance of combinations of methylation and expression of nine DMR-DEGs, which were more efficient prognostic biomarkers than considering either type of data alone. Finally, we validated 22/23 (95.7%) DMR-DEGs in 5-azadC-treated LO2 and/or HepG2 cells. In conclusion, integration of epigenome and transcriptome data depicted activation of multiple pivotal cell cycle-related pathways and repression of several metabolic pathways triggered by aberrant DNA methylation of promoters and enhancers in HCC.

25-Hydroxycholesterol 3-Sulfate Recovers Acetaminophen Induced Acute Liver Injury via Stabilizing Mitochondria in Mouse Models

Acetaminophen (APAP) overdose is one of the most frequent causes of acute liver failure (ALF). N-acetylcysteine (NAC) is currently being used as part of the standard care in the clinic but its usage has been limited in severe cases, in which liver transplantation becomes the only treatment option. Therefore, there still is a need for a specific and effective therapy for APAP induced ALF. In the current study, we have demonstrated that treatment with 25-Hydroxycholesterol 3-Sulfate (25HC3S) not only significantly reduced mortality but also decreased the plasma levels of liver injury markers, including LDH, AST, and ALT, in APAP overdosed mouse models. 25HC3S also decreased the expression of those genes involved in cell apoptosis, stabilized mitochondrial polarization, and significantly decreased the levels of oxidants, malondialdehyde (MDA), and reactive oxygen species (ROS). Whole genome bisulfite sequencing analysis showed that 25HC3S increased demethylation of 5mCpG in key promoter regions and thereby increased the expression of those genes involved in MAPK-ERK and PI3K-Akt signaling pathways. We concluded that 25HC3S may alleviate APAP induced liver injury via up-regulating the master signaling pathways and maintaining mitochondrial membrane polarization. The results suggest that 25HC3S treatment facilitates the recovery and significantly decreases the mortality of APAP induced acute liver injury and has a synergistic effect with NAC in propylene glycol (PG) for the injury.