Bisulfite Genomic Sequencing Research Articles

Integrative analysis of whole genome bisulfite and transcriptome sequencing reveals the effect of sodium butyrate on DNA methylation in the differentiation of bovine skeletal muscle satellite cells

Integrative analysis of whole genome bisulfite and transcriptome sequencing reveals the effect of sodium butyrate on DNA methylation in the differentiation of bovine skeletal muscle satellite cells

Methylome of circulating cell free DNA as a novel biomarker tool: from acute coronary syndrome to broader risk stratification of cardiovascular disease

Abstract Background Cell death during tissue injury leads to release of chromosomal DNA which is promptly degraded. However, short DNA fragments wrapped around nucleosomes survive longer in the bloodstream and can be isolated as circulating cell-free DNA (ccfDNA). Oncology field pioneered the use of cfDNA as clinical biomarker by identifying sequence modifications induced in cancer cells. Importantly, ccfDNA also retains the same epigenetic information as the tissue of origin. Despite major advances for disease diagnosis and therapy, coronary artery diseases (CAD) and associated conditions such acute myocardial infarction (AMI) or heart failure (HF) are responsible for leading number of deaths worldwide. Preventive and surveillance strategies with non-invasive testing are needed. Purpose Evaluation of ccfDNA methylation profiles as prognostic non-invasive biomarker for risk stratification of acute coronary syndromes and their potential application in preventive and surveillance strategies in other cardiovascular disorders. Methods CcfDNA isolated from citrate-plasma collected from our discovery cohort, comprising healthy controls, STEMI, NSTEMI and unstable angina (UA) patients, was analyzed applying Whole Genome Bisulfite Sequencing (WGBS) using a low-input BS-seq (PBAT) protocol. Computational analysis generated differentially methylated regions (DMRs) by tiling CpG areas and setting a threshold of minimal 25% methylation difference compared to the healthy group. To assess the validity of our findings, a selection of disease-specific DMRs were tested using targeted methylation sequencing protocol in a new set of patients. This procedure consisted in enzymatic conversion of ccfDNA methylated sites followed by probe-specific enrichment libraries and sequencing at a Novaseq 6000 platform. Results The study identified a total of 1941 DMRs and from those 486 were STEMI, 223 NSTEMI and 684 UA specific. From the identified DMRs, 74.13% from STEMI, 72.68% from NSTEMI and 61.62% from UA annotated with heart-related diseases using DisGeNET gene-disease associations. Interestingly, most of identified DMRs located within intronic regions (~60%) and also annotated as enhancers and lncRNAs. The independent validation with targeted methylation sequencing identified 566, 306 and 248 DMRs, respectively, and reduced the number of disease-specific DMRs to 171 in STEMI, 115 in NSTEMI and 82 in UA. Conclusion and Future Approaches Methylation profiles of ccfDNA showed potential as a biomarker to stratify the severity of acute coronary syndromes, a novel non-invasive detection methodology that could be extrapolated to risk stratification in other cardiovascular disease such as heart failure or coronary artery disease.

PCBS: an R package for fast and accurate analysis of bisulfite sequencing data.

Whole-genome bisulfite sequencing is a powerful tool for analyzing chromatin methylation genome-wide, but analysis of whole-genome bisulfite data is slow, inflexible, and often inaccurate. We developed PCBS (Principal Component BiSulfite), a computationally efficient R package for Whole Genome Bisulfite Sequencing analysis that demonstrates remarkable accuracy and flexibility compared to current tools. PCBS identifies differentially methylated loci, differentially methylated regions, and offers novel functionality that allows for more targeted methylation analyses. PCBS uses minimal computational resources; a complete pipeline in mouse can run on a local RStudio instance in a matter of minutes. PCBS is an R package available under a GNU GPLv3 license on GitHub: https://github.com/katlande/PCBS and CRAN: https://CRAN.R-project.org/package=PCBS. Instructions for use are available at: https://katlande.github.io/PCBS/.

PRKACB is a novel imprinted gene in marsupials

BackgroundGenomic imprinting results in parent-of-origin-specific gene expression and, among vertebrates, is found only in therian mammals: marsupials and eutherians. A differentially methylated region (DMR), in which the methylation status of CpG dinucleotides differs between the two alleles, can mark the parental identity of imprinted genes. We developed a computational pipeline that detected CpG islands (CGIs) marked by both methylated and unmethylated signals in whole genome bisulfite sequencing data. This approach identified candidate marsupial DMRs in a publicly available koala methylome. One of these candidate DMRs was associated with PRKACB, a gene encoding the protein kinase A catalytic subunit beta. Nothing is known about the imprinting status of PRKACB in eutherian mammals although mutations of this gene are associated with endocrine neoplasia and other developmental disorders.ResultsIn the tammar wallaby and brushtail possum there was parent-of-origin-specific DNA methylation in the PRKACB DMR in which the maternal allele was methylated and the paternal allele was unmethylated. There were multiple RNAs transcribed from this locus. Allele-specific expression analysis identified paternal expression of a PRKACB lncRNA and an mRNA isoform. Comparison of the PRKACB gene start site between marsupials and eutherians demonstrated that the CGI is longer in marsupials. The PRKACB gene product functions in the same signalling pathway as the guanine nucleotide-binding protein alpha subunit encoded at the GNAS locus, a known eutherian imprinted gene. In a mouse methylome Gnas had three differentially methylated CGIs, while in the koala methylome the GNAS locus had two unmethylated CGIs.ConclusionsWe conclude that PRKACB is a novel, DMR-associated marsupial imprinted gene. Imprinting of PRKACB in marsupials and GNAS in eutherians may indicate a conserved selection pressure for imprinting of the protein kinase A signalling pathway in therians with the two lineages adapting by imprinting different genes.

DNA methylation correlates with transcriptional noise in response to elevated pCO2 in the eastern oyster (Crassostrea virginica)

Abstract Ocean acidification significantly affects marine calcifiers like oysters, warranting the study of molecular mechanisms like DNA methylation that contribute to adaptive plasticity in response to environmental change. However, a consensus has not been reached on the extent to which methylation modules gene expression, and in turn plasticity, in marine invertebrates. In this study, we investigated the impact of pCO2 on gene expression and DNA methylation in the eastern oyster, Crassostrea virginica. After a 30-day exposure to control (572 ppm) or elevated pCO2 (2,827 ppm), whole genome bisulfite sequencing (WGBS) and RNA-Seq data were generated from adult female gonad tissue and male sperm samples. Although differentially methylated loci (DML) were identified in females (89) and males (2,916), there were no differentially expressed genes, and only one differentially expressed transcript in females. However, gene body methylation impacted other forms of gene activity in sperm, such as the maximum number of transcripts expressed per gene and changes in the predominant transcript expressed. Elevated pCO2 exposure increased gene expression variability (transcriptional noise) in males but decreased noise in females, suggesting a sex-specific role of methylation in gene expression regulation. Functional annotation of genes with changes in transcript-level expression or containing DML revealed several enriched biological processes potentially involved in elevated pCO2 response, including apoptotic pathways and signal transduction, as well as reproductive functions. Taken together, these results suggest that DNA methylation may regulate gene expression variability to maintain homeostasis in elevated pCO2 conditions and could play a key role in environmental resilience in marine invertebrates.

Genome-wide DNA methylation sequencing reveals the involvement of ferroptosis in hepatotoxicity induced by dietary exposure to food-grade titanium dioxide

BackgroundFollowing the announcement by the European Food Safety Authority that the food additive titanium dioxide (E 171) is unsafe for human consumption, and the subsequent ban by the European Commission, concerns have intensified over the potential risks E 171 poses to human vital organs. The liver is the main organ for food-grade nanoparticle metabolism. It is increasingly being found that epigenetic changes may play an important role in nanomaterial-induced hepatotoxicity. However, the profound effects of E 171 on the liver, especially at the epigenetic level, remain largely unknown.MethodsMice were exposed orally to human-relevant doses of two types of E 171 mixed in diet for 28 and/or 84 days. Conventional toxicology and global DNA methylation analyses were performed to assess E 171-induced hepatotoxicity and epigenetic changes. Whole genome bisulfite sequencing and further ferroptosis protein detection were used to reveal E 171-induced changes in liver methylation profiles and toxic mechanisms.ResultsExposed to E 171 for 28 and/or 84 days resulted in reduced global DNA methylation and hydroxymethylation in the liver of mice. E 171 exposure for 84 days elicited inflammation and damage in the mouse liver, whereas 28-day exposure did not. Whole-genome DNA methylation sequencing disclosed substantial methylation alterations at the CG and non-CG sites of the liver DNA in mice exposed to E 171 for 84 days. Mechanistic analysis of the DNA methylation alterations indicated that ferroptosis contributed to the liver toxicity induced by E 171. E 171-induced DNA methylation changes triggered NCOA4-mediated ferritinophagy, attenuated the protein levels of GPX4, FTH1, and FTL in the liver, and thereby caused ferroptosis.ConclusionsLong-term oral exposure to E 171 triggers hepatotoxicity and induces methylation changes in both CG and non-CG sites of liver DNA. These epigenetic alterations activate ferroptosis in the liver through NCOA4-mediated ferritinophagy, highlighting the role of DNA methylation and ferroptosis in the potential toxicity caused by E 171 in vivo.

DNA Methylation Carries Signatures of Sublethal Effects Under Thermal Stress in Loggerhead Sea Turtles.

To date, studies of the impacts of climate warming on individuals and populations have mostly focused on mortality and thermal tolerance. In contrast, much less is known about the consequences of sublethal effects, which are more challenging to detect, particularly in wild species with cryptic life histories. This necessitates the development of molecular tools to identify their signatures. In a split-clutch field experiment, we relocated clutches of wild, nesting loggerhead sea turtles (Caretta caretta) to an insitu hatchery. Eggs were then split into two sub-clutches and incubated under shallow or deep conditions, with those in the shallow treatment experiencing significantly higher temperatures in otherwise natural conditions. Although no difference in hatching success was observed between treatments, hatchlings from the shallow, warmer treatment had different length-mass relationships and were weaker at locomotion tests than their siblings incubated in the deep, cooler treatment. To characterise the molecular signatures of these thermal effects, we performed whole genome bisulfite sequencing on blood samples collected upon emergence. We identified 287 differentially methylated sites between hatchlings from different treatments, including on genes with neurodevelopmental, cytoskeletal, and lipid metabolism functions. Taken together, our results show that higher incubation temperatures induce sublethal effects in hatchlings, which are reflected in their DNA methylation status at identified sites. These sites could be used as biomarkers of thermal stress, especially if they are retained across life stages. Overall, this study suggests that global warming reduces hatchling fitness, which has implications for dispersal capacity and ultimately a population's adaptive potential. Conservation efforts for these endangered species and similar climate-threatened taxa will therefore benefit from strategies for monitoring and mitigating exposure to temperatures that induce sublethal effects.

Role of Epigenetic Factors in Determining the Biological Behavior and Prognosis of Hepatocellular Carcinoma.

The current study's objective is to evaluate the molecular genetic mechanisms influencing the biological behavior of hepatocellular carcinoma (HCC) by analyzing the transcriptomic and epigenetic signatures of the tumors. Transcriptomic data were downloaded from the NCBI GEO database. We investigated the expression differences between the GSE46444 (48 cirrhotic tissues versus 88 HCC tissues) and GSE63898 (168 cirrhotic tissues versus 228 HCC tissues) data sets using GEO2R. Differentially expressed genes were evaluated using GO and KEGG metabolic pathway analysis websites. Whole genome bisulfite sequencing (WGBS) and Methylated DNA Immunoprecipitation Sequencing (MeDIP-Seq) data sets (26 HCC tissues versus 26 adjacent non-tumoral tissues) were also downloaded from the NCBI SRA database. These data sets were analyzed using Bismark and QSEA, respectively. The methylation differences between the groups were assessed using functional enrichment analysis. In the GSE46444 data set, 80 genes were upregulated, and 315 genes were downregulated in the tumor tissue (HCC tissue) compared to the non-tumor cirrhotic tissue. In the GSE63898 data set, 1261 genes were upregulated, and 458 genes were downregulated in the cirrhotic tissue compared to the tumor tissues. WGBS revealed that 20 protein-coding loci were hypermethylated. while the hypomethylated regions were non-protein-coding. The methylated residues of the tumor tissue, non-tumorous cirrhotic tissue, and healthy tissue were comparable. MeDIP-Seq, conducted on tumoral and non-tumoral tissues, identified hypermethylated or hypomethylated areas as protein-coding regions. The functional enrichment analysis indicated that these genes were related to pathways including peroxisome, focal adhesion, mTOR, RAP1, Phospholipase D, Ras, and PI3K/AKT signal transduction. The investigation of transcriptomic and epigenetic mechanisms identified several genes significant in the biological behavior of HCC. These genes present potential targets for the development of targeted therapy.

Comprehensive Analysis of Methylome and Transcriptome to Identify Potential Genes Regulating Porcine Testis Development.

DNA methylation plays a critical role in regulating gene expression during testicular development. However, few studies report on candidate genes related to the DNA methylation regulation of porcine testicular development. This study examined the differentially expressed genes (DEGs) and their methylation levels in testicular tissues from pigs at 60 days of age (60 d) and 180 days of age (180 d) using RNA-Seq and whole genome bisulfite sequencing (WGBS). It was determined that DNA methylation primarily occurs in the cytosine-guanine (CG) context, and the analysis identified 106,282 differentially methylated regions (DMRs) corresponding to 12,385 differentially methylated genes (DMGs). Further integrated analysis of RNA-Seq and WGBS data revealed 1083 DMGs negatively correlated with the expression of DEGs. GO analysis showed that these genes were significantly enriched in spermatogenesis, germ cell development, and spermatid differentiation. The screening of enriched genes revealed that hyper-methylation repressed ADAM30, ADAM3A, DPY19L2, H2BC1, MAK, RPL10L, SPATA16, and YBX2, while hypo-methylation elevated CACNA1I, CADM1, CTNNB1, JAM2, and PAFAH1B3 expression. Additionally, the methylation status of the key genes ADAM3A, ADAM30, YBX2, JAM2, PAFAH1B3, and CTNNB1 was detected by bisulfite sequencing PCR (BSP). This study offers insights into the epigenetic regulation mechanisms underlying porcine testicular development.

DNA methylation analysis reveals local changes in resistant and susceptible soybean lines in response to Phytophthora sansomeana.

Phytophthora sansomeana is an emerging oomycete pathogen causing root rot in many agricultural species including soybean. However, as of now, only one potential resistance gene has been identified in soybean, and our understanding of how genetic and epigenetic regulation in soybean contributes to responses against this pathogen remains largely unknown. In this study, we performed whole genome bisulfite sequencing (WGBS) on two soybean lines, Colfax (resistant) and Williams 82 (susceptible), in response to P. sansomeana at two time points: 4 and 16 hours post-inoculation to compare their methylation changes. Our findings revealed that there were no significant changes in genome-wide CG, CHG (H = A, T, or C), and CHH methylation. However, we observed local methylation changes, specially an increase in CHH methylation around genes and transposable elements (TEs) after inoculation, which occurred earlier in the susceptible line and later in the resistant line. After inoculation, we identified differentially methylated regions (DMRs) in both Colfax and Williams 82, with a predominant presence in TEs. Notably, our data also indicated that more TEs exhibited changes in their methylomes in the susceptible line compared to the resistant line. Furthermore, we discovered 837 DMRs within or flanking 772 differentially expressed genes (DEGs) in Colfax and 166 DMRs within or flanking 138 DEGs in Williams 82. These DEGs had diverse functions, with Colfax primarily showing involvement in metabolic process, defense response, plant and pathogen interaction, anion and nucleotide binding, and catalytic activity, while Williams 82 exhibited a significant association with photosynthesis. These findings suggest distinct molecular responses to P. sansomeana infection in the resistant and susceptible soybean lines.

Females with autism spectrum disorders show stronger DNA methylation signatures than males in perinatal tissues.

Autism spectrum disorder (ASD) comprises a group of neurodevelopmental conditions currently diagnosed by behavioral assessment in childhood, with reported underdiagnosis in females. Though diagnosis in early life is linked to improved outcomes, we currently lack objective screening tools for newborns. To address this gap, we sought to identify a sex-specific DNA methylation signature for ASD using perinatal tissues that reflect dysregulation in the brain. DNA methylation was assayed from ASD and typically developing (TD) newborn blood, umbilical cord blood, placenta, and post-mortem cortex samples using whole genome bisulfite sequencing (WGBS) in a total of 511 samples. We found that methylation levels of differentially methylated regions (DMRs) differentiated samples by ASD diagnosis in females more than males across the perinatal tissues. We tested three theories for ASD sex differences in newborn blood, finding epigenetic support for an X chromosome-related female protective effect, as well as a high replication rate of DMRs (48.1%) in females across two independent cohorts. In our pan-tissue analysis, three genes (X-linked BCOR, GALNT9, OPCML) mapped to ASD DMRs replicated in all four female tissues. ASD DMRs from all tissues were enriched for neuro-related processes (females) and SFARI ASD-risk genes (females and males). Overall, we found a highly replicated methylation signature of ASD in females across perinatal tissues that reflected dysregulation in the brain and involvement of X chromosome epigenetics. This comparative study of perinatal tissues shows the promise of newborn blood DNA methylation biomarkers for early detection of females at risk for ASD and emphasizes the importance of sex-stratification in ASD studies.

Creation and validation of the first infinium DNA methylation array for the human imprintome

BackgroundDifferentially methylated imprint control regions (ICRs) regulate the monoallelic expression of imprinted genes. Their epigenetic dysregulation by environmental exposures throughout life results in the formation of common chronic diseases. Unfortunately, existing Infinium methylation arrays lack the ability to profile these regions adequately. Whole genome bisulfite sequencing (WGBS) is the unique method able to profile the ICRs. However, it is very expensive and it requires not only a high coverage, but it is also computationally intensive to assess these regions.FindingsTo address this deficiency, we developed a custom methylation array containing 22,819 probes. Among them, 10,438 are CG probes targeting unique CpG sites, with 9,757 probes successfully mapping to 1,088 out of the 1,488 candidate ICRs recently described. To assess the performance of the array, we created matched samples processed with the Human Imprintome array and WGBS, which is the current standard method for assessing the methylation of the Human Imprintome. We compared the methylation levels from the shared CpG sites, and obtained a mean R2 = 0.569. We also created matched samples processed with the Human Imprintome array and the Infinium Methylation EPIC v2 array, and obtained a mean R2 = 0.796. Furthermore, replication experiments demonstrated high reliability (ICC: 0.799–0.945).ConclusionsOur custom array will be useful for replicable and accurate assessment, mechanistic insight, and targeted investigation of ICRs. This tool should accelerate the discovery of ICRs associated with a wide range of diseases and exposures, and advance our understanding of genomic imprinting and its relevance in development and disease formation throughout the life course.

SETDB1 regulates short interspersed nuclear elements and chromatin loop organization in mouse neural precursor cells

BackgroundTransposable elements play a critical role in maintaining genome architecture during neurodevelopment. Short Interspersed Nuclear Elements (SINEs), a major subtype of transposable elements, are known to harbor binding sites for the CCCTC-binding factor (CTCF) and pivotal in orchestrating chromatin organization. However, the regulatory mechanisms controlling the activity of SINEs in the developing brain remains elusive.ResultsIn our study, we conduct a comprehensive genome-wide epigenetic analysis in mouse neural precursor cells using ATAC-seq, ChIP-seq, whole genome bisulfite sequencing, in situ Hi-C, and RNA-seq. Our findings reveal that the SET domain bifurcated histone lysine methyltransferase 1 (SETDB1)-mediated H3K9me3, in conjunction with DNA methylation, restricts chromatin accessibility on a selective subset of SINEs in neural precursor cells. Mechanistically, loss of Setdb1 increases CTCF access to these SINE elements and contributes to chromatin loop reorganization. Moreover, de novo loop formation contributes to differential gene expression, including the dysregulation of genes enriched in mitotic pathways. This leads to the disruptions of cell proliferation in the embryonic brain after genetic ablation of Setdb1 both in vitro and in vivo.ConclusionsIn summary, our study sheds light on the epigenetic regulation of SINEs in mouse neural precursor cells, suggesting their role in maintaining chromatin organization and cell proliferation during neurodevelopment.

Multi-omic characterization of allele-specific regulatory variation in hybrid pigs

Hybrid mapping is a powerful approach to efficiently identify and characterize genes regulated through mechanisms in cis. In this study, using reciprocal crosses of the phenotypically divergent Duroc and Lulai pig breeds, we perform a comprehensive multi-omic characterization of regulatory variation across the brain, liver, muscle, and placenta through four developmental stages. We produce one of the largest multi-omic datasets in pigs to date, including 16 whole genome sequenced individuals, as well as 48 whole genome bisulfite sequencing, 168 ATAC-Seq and 168 RNA-Seq samples. We develop a read count-based method to reliably assess allele-specific methylation, chromatin accessibility, and RNA expression. We show that tissue specificity was much stronger than developmental stage specificity in all of DNA methylation, chromatin accessibility, and gene expression. We identify 573 genes showing allele specific expression, including those influenced by parent-of-origin as well as allele genotype effects. We integrate methylation, chromatin accessibility, and gene expression data to show that allele specific expression can be explained in great part by allele specific methylation and/or chromatin accessibility. This study provides a comprehensive characterization of regulatory variation across multiple tissues and developmental stages in pigs.

Genetic and epigenetic differentiation in response to genomic selection for avian lay date.

Anthropogenic climate change has led to globally increasing temperatures at an unprecedented pace and, to persist, wild species have to adapt to their changing world. We, however, often fail to derive reliable predictions of species' adaptive potential. Genomic selection represents a powerful tool to investigate the adaptive potential of a species, but constitutes a 'blind process' with regard to the underlying genomic architecture of the relevant phenotypes. Here, we used great tit (Parus major) females from a genomic selection experiment for avian lay date to zoom into this blind process. We aimed to identify the genetic variants that responded to genomic selection and epigenetic variants that accompanied this response and, this way, might reflect heritable genetic variation at the epigenetic level. We applied whole genome bisulfite sequencing to blood samples of individual great tit females from the third generation of bidirectional genomic selection lines for early and late lay date. Genomic selection resulted in differences at both the genetic and epigenetic level. Genetic variants that showed signatures of selection were located within genes mostly linked to brain development and functioning, including LOC107203824 (SOX3-like). SOX3 is a transcription factor that is required for normal hypothalamo-pituitary axis development and functioning, an essential part of the reproductive axis. As for epigenetic differentiation, the early selection line showed hypomethylation relative to the late selection line. Sites with differential DNA methylation were located in genes important for various biological processes, including gonadal functioning (e.g., MSTN and PIK3CB). Overall, genomic selection for avian lay date provided insights into where within the genome the heritable genetic variation for lay date, on which selection can operate, resides and indicates that some of this variation might be reflected by epigenetic variants.

Oral pimonidazole unveils clinicopathologic and epigenetic features of hypoxic tumour aggressiveness in localized prostate cancer

BackgroundTumor hypoxia is associated with prostate cancer (PCa) treatment resistance and poor prognosis. Pimonidazole (PIMO) is an investigational hypoxia probe used in clinical trials. A better understanding of the clinical significance and molecular alterations underpinning PIMO-labeled tumor hypoxia is needed for future clinical application. Here, we investigated the clinical significance and molecular alterations underpinning PIMO-labeled tumor hypoxia in patients with localized PCa, in order to apply PIMO as a prognostic tool and to identify potential biomarkers for future clinical translation.MethodsA total of 39 patients with localized PCa were recruited and administered oral PIMO before undergoing radical prostatectomy (RadP). Immunohistochemical staining for PIMO was performed on 37 prostatectomy specimens with staining patterns evaluated and clinical association analyzed. Whole genome bisulfite sequencing was performed using laser-capture of microdissected specimen sections comparing PIMO positive and negative tumor areas. A hypoxia related methylation molecular signature was generated by integrating the differentially methylated regions with previously established RNA-seq datasets.ResultsThree PIMO staining patterns were distinguished: diffuse, focal, and comedo-like. The comedo-like staining pattern was more commonly associated with adverse pathology. PIMO-defined hypoxia intensity was positively correlated with advanced pathologic stage, tumor invasion, and cribriform and intraductal carcinoma morphology. The generated DNA methylation signature was found to be a robust hypoxia biomarker, which could risk-stratify PCa patients across multiple clinical datasets, as well as be applicable in other cancer types.ConclusionsOral PIMO unveiled clinicopathologic features of disease aggressiveness in localized PCa. The generated DNA methylation signature is a novel and robust hypoxia biomarker that has the potential for future clinical translation.

ATRT-16. MAPPING THE DISORDERED METHYLOME IN ATRT

Abstract BACKGROUND Atypical Teratoid Rhabdoid Tumor (ATRT) is an aggressive pediatric brain tumor driven by inactivation of the chromatin regulator SMARCB1, with remarkably few other mutations. SMARCB1 loss leads to epigenetic dysregulation which plays a pivotal role in ATRT development and progression, and DNA methylation profiling has distinguished molecular subgroups of this disease. However, studies that provide a static view of the epigenome may not capture epigenetic variability that drives cellular plasticity and heterogeneity. Mapping epigenetic variability is crucial to understanding the phenotypic plasticity of this tumor and its ability to evade therapeutic interventions in the clinic. METHODS We performed whole genome bisulfite sequencing (WGBS) on 22 ATRT samples, 4 normal fetal brain controls, and 1 hindbrain neural stem cell control, and applied a novel analytic pipeline (InformME) to generate DNA methylation potential energy landscapes encapsulating epigenetic stochasticity. This characterizes mean methylation as well as methylation entropy. We identified key genes and regulatory regions with the most stochastic methylation signatures. RESULTS We find that ATRT exhibits a globally destabilized epigenome, with markedly elevated methylation entropy (a measure of variability/disorder). Methylation discordance mapped to key regulatory regions including bivalent promoters, quiescent chromatin regions, and enhancers, and specific genes, with an enrichment for targets of the Polycomb repressive complex 2 and regulators of pluripotency and differentiation. CONCLUSIONS The disordered methylome of ATRT could play a critical role in the tumor’s ability to evade therapies. Mapping the targets of epigenetic instability can identify targets for intervention. Moreover, strategies to constrain epigenetic variability may have therapeutic benefit by limiting the number of epigenetic states available to cancer cells.

DIPG-58. TARGETING DISORDERED DNA METHYLATION IN DIPG TO CONSTRAIN VARIABILITY AND INDUCE IMMUNE SIGNALING

Abstract BACKGROUND Diffuse Intrinsic Pontine Glioma (DIPG) is a universally fatal pediatric brain cancer characterized by the histone H3 K27M mutation. The known downstream consequences of this mutation include loss of repressive chromatin marks, global DNA hypomethylation, and altered gene expression. METHODS We sought to investigate the role of epigenetic variability as a basis of DIPG cellular heterogeneity and plasticity by using Whole Genome Bisulfite Sequencing (WGBS). Existing methods of analysis of WGBS are limited in that they average/smooth methylation data and fail to capture intrinsic variability. We performed (WGBS) on 23 primary patient samples of DIPG and applied an analysis that captures the variability of DNA methylation reads to characterize methylation stochasticity. We then sought to modulate the epigenome pharmacologically using a DNA methyltransferase inhibitor and assessed changes through RNA-Seq and immunopeptidomics. RESULTS We find that DIPG has a marked increase in methylation stochasticity (quantified as methylation entropy). Additionally, while DIPG is globally hypomethylated, it retains DNA methylation at specific regulatory regions, such as bivalent promoters, and key genes, such as FOXG1, CXCR4, and KLF4. We sought to reverse these changes using the DNA methyltransferase inhibitor, decitabine, in 4 DIPG neurosphere cell lines. Treatment with decitabine reversed the entropy seen in DIPG and induced innate immune signaling and interferon signaling in DIPG cells. Using RNA-seq, we found that decitabine treatment dramatically alters gene expression, including activation of endogenous retroviral elements, activation of interferon signaling, de-repression of hypermethylated targets, and expression of putative neoantigens, such as PRAME and DAZL. Immunopeptidomics profiling of MHC Class I bound canonical and non-canonical peptides of DIPG cell lines after decitabine treatment also revealed decitabine-induced differences. CONCLUSIONS Taken together, our study finds that the methylome of DIPG is highly disordered but is also highly responsive to pharmacologic modulation by hypomethylating agents. Furthermore, hypomethylating drugs can increase the immunogenicity of DIPG, thus offering the potential for future combination with immunotherapy.

PATH-10. ACCELERATING COMPREHENSIVE CNS TUMOR MOLECULAR DIAGNOSTICS WITH RAPID-CNS2 AND MNP-FLEX: A PROSPECTIVE MULTI-CENTER VALIDATION

Abstract BACKGROUND The 2021 WHO classification update highlights the necessity of integrating molecular alterations for precise central nervous system (CNS) tumor diagnoses. However, current molecular reporting methods are hindered by significant initial investment, labor-intensive protocols, and lengthy turnaround times. Methylation-based classification has emerged as a pivotal diagnostic tool but is currently limited to array-based techniques. This necessitates exploration of novel technologies to streamline molecular analysis. METHODS We implemented Rapid-CNS2 - our adaptive sampling-based nanopore sequencing workflow- on 190 adult and pediatric samples at University Hospital Heidelberg and University of Nottingham. Intraoperative potential was assessed through real-time analysis followed by 24-hour sequencing for comprehensive genomic insights. Additionally, we developed MNP-Flex, a platform-agnostic version of the Heidelberg methylation classifier covering 184 CNS tumor classes. We evaluated MNP-flex on a global cohort of over 78,000 samples from methylation arrays, whole genome bisulfite sequencing, nanopore whole genome sequencing, methylation panels and Rapid-CNS2. RESULTS Rapid-CNS2 validation yielded accurate integrated diagnoses in all 190 samples. Within a crucial 30-minute timeframe, we reported accurate methylation families and arm-level copy number profiles followed by next-day reporting of fine-grained methylation classification, SNVs, focal CNVs, MGMT status, fusions and novel structural variants. Moreover, MNP-Flex achieved 92% accuracy over the validation dataset spanning over 78,000 samples from five different technologies. CONCLUSIONS The adoption of Rapid-CNS2 and MNP-Flex enables rapid intraoperative broad methylation classification and copy number alteration reporting within 30 minutes, with additional clinically relevant, fine-grained molecular insights available the following day. It offers clinicians rapid access to comprehensive molecular information critical for treatment decisions. Furthermore, MNP-Flex extends the utility of the Heidelberg methylation classifier to diverse sequencing-based data. By overcoming the limitations of currently available methods, our workflow represents a paradigm shift in the field, promising improved management of CNS tumor patients. Rapid-CNS2 can be executed with a single command, while MNP-Flex is publicly available as a web service, enhancing accessibility and usability for clinical applications.

Azacitidine as epigenetic priming for chemotherapy is safe and well-tolerated in infants with newly diagnosed KMT2A-rearranged acute lymphoblastic leukemia: Children's Oncology Group trial AALL15P1.

Infants less than 1 year old diagnosed with KMT2A-rearranged (KMT2A-r) acute lymphoblastic leukemia (ALL) are at high risk of remission failure, relapse, and death due to leukemia, despite intensive therapies. Infant KMT2A-r ALL blasts are characterized by DNA hypermethylation. Epigenetic priming with DNA methyltransferase inhibitors increases the cytotoxicity of chemotherapy in preclinical studies. The Children's Oncology Group trial AALL15P1 tested the safety and tolerability of five days of azacitidine immediately prior to the start of chemotherapy on day six, in four post-induction chemotherapy courses for infants with newly diagnosed KMT2A-r ALL. The treatment was welltolerated, with only two of 31 evaluable patients (6.5%) experiencing dose-limiting toxicity. Whole genome bisulfite sequencing of peripheral blood mononuclear cells (PBMCs) demonstrated decreased DNA methylation in 87% of samples tested following five days of azacitidine. Event-free survival was similar to prior studies of newly diagnosed infant ALL. Azacitidine is safe and results in decreased DNA methylation of PBMCs in infants with KMT2A-r ALL, but the incorporation of azacitidine to enhance cytotoxicity did not impact survival. Clinicaltrials.gov identifier: NCT02828358.