Genome DNA Methylation Analysis Research Articles

DNA methylation analysis in patients with neurodevelopmental disorders improves variant interpretation and reveals complexity

Analysis of genomic DNA methylation by generating epigenetic signature profiles (“episignatures”) is increasingly being implemented in genetic diagnosis. Here we report our experience using episignature analysis to resolve both uncomplicated and complex cases of neurodevelopmental disorder (NDD). We analysed 97 NDDs divided into: (i) a validation cohort of 59 patients with likely pathogenic/pathogenic variants characterized by a known episignature and (ii) a test cohort of 38 patients harbouring variants of unknown significance (VUS) or unidentified variants. The expected episignature was obtained in most cases with likely pathogenic/pathogenic variants (53/59; 90%), a revealing exception being the overlapping profile of two SMARCB1 pathogenic variants with ARID1A/B:c.6200, confirmed by the overlapping clinical features. In the test cohort, five cases showed the expected episignature, including: (i) novel pathogenic variants in ARID1B and BRWD3; (ii) a deletion in ATRX causing MRXFH1 X-linked mental retardation and (iii) confirmed the clinical diagnosis of Cornelia de Lange (CdL) syndrome in mutation negative CdL patients. Episignatures analysis of the in BAF complex components revealed novel functional protein interactions and common episignatures affecting homologous residues in highly conserved paralogous proteins (SMARCA2 M856V and SMARCA4 M866V). Finally, we also found sex-dependent episignatures in X-linked disorders. Implementation of episignature profiling is still in its early days but with increasing utilization come increasing awareness of the capacity of this methodology to help resolve the complex challenges of genetic diagnoses.

Evaluating the utility of ZNF331 promoter methylation as a prognostic and predictive marker in stage III colon cancer: results from CALGB 89803 (Alliance)

ABSTRACT While epigenomic alterations are common in colorectal cancers (CRC), few epigenomic biomarkers that risk-stratify patients have been identified. We thus sought to determine the potential of ZNF331 promoter hypermethylation (mZNF331) as a prognostic and predictive marker in colon cancer. We examined the association of mZNF331 with clinicopathologic features, relapse, survival, and treatment efficacy in patients with stage III colon cancer treated within a randomized adjuvant chemotherapy trial (CALGB/Alliance89803). Residual tumour tissue was available for genomic DNA extraction and methylation analysis for 385 patients. ZNF331 promoter methylation status was determined by bisulphite conversion and fluorescence-based real-time polymerase chain reaction. Kaplan-Meier estimator and Cox proportional hazard models were used to assess the prognostic and predictive role of mZNF331 in this well-annotated dataset, adjusting for clinicopathologic features and standard molecular markers. mZNF331 was observed in 267/385 (69.4%) evaluable cases. Histopathologic features were largely similar between patients with mZNF331 compared to unmethylated ZNF331 (unmZNFF31). There was no significant difference in disease-free or overall survival between patients with mZNF331 versus unmZNF331 colon cancers, even when adjusting for clinicopathologic features and molecular marker status. Similarly, there was no difference in disease-free or overall survival across treatment arms when stratified by ZNF331 methylation status. While ZNF331 promoter hypermethylation is frequently observed in CRC, our current study of a small subset of patients with stage III colon cancer suggests limited applicability as a prognostic marker. Larger studies may provide more insight and clarity into the applicability of mZNF331 as a prognostic and predictive marker.

270 Inducing Neural Maturation and Increasing Survival in SHH Medulloblastoma Using EZH2 Inhibition in a Mouse Model

INTRODUCTION: Medulloblastoma (MB) is the most common malignant pediatric brain tumor. Genetically, MB can be divided into four subgroups of which the SHH subtype histologically shows nodular architecture. Within this nodular architecture, there are islands of mature cells, with more abundant neuropil and a low proliferation rate, scattered among sheets of primitive cells. METHODS: We performed laser capture microdissection followed by whole transcriptome analysis, spatial transcriptomics using Digital Spatial Profiling, whole genome DNA methylation and ChIP-Seq analysis of mature and primitive areas from 8 medulloblastomas. We developed a genetically-engineered mouse model of SHH MB showing spontaneous maturation and lack of maturation with a conditional EZH2 genetic ablation or EZH2 overactivation respectively. Finally, we developed a fucoidan-based nanoparticle drug delivery across the blood brain barrier (BBB) for targeted molecular inhibition. RESULTS: Using whole transcriptome and DNA methylation analysis, we identified ∼120 differentially expressed genes between primitive and mature regions with enrichment for genes regulated by H3K4me3 and H3K27me3. ChIP-Seq analysis showed striking differences in H3K27me3 enrichment between primitive and mature medulloblastoma cells including at the EZH2 locus. Medulloblastoma specific EZH2 genetic ablation resulted in diffuse tumor cell differentiation and prolonged survival in mice (n = 10 per group, log-rank p = 0.01). Conversely, conditional EZH2 (Y641F) activation prevented medulloblastoma differentiation. A fucoidan-based nanoparticle successfully delivered the EZH2 inhibitor (EPZ-6438) across the murine BBB to achieve significant extension of mouse survival (median 70 days compared to 21 days in control mice; *p = 0.01, Mantel-Cox). CONCLUSIONS: Spontaneous maturation of medulloblastoma cells can be induced by inhibition of EZH2. Fucoidan-based nanoparticle delivery systems allows tumor specific targeted delivery across the BBB extending survival in mice.

Genetic and epigenetic features of bilateral Wilms tumor predisposition in patients from the Children’s Oncology Group AREN18B5-Q

Developing synchronous bilateral Wilms tumor suggests an underlying (epi)genetic predisposition. Here, we evaluate this predisposition in 68 patients using whole exome or genome sequencing (n = 85 tumors from 61 patients with matched germline blood DNA), RNA-seq (n = 99 tumors), and DNA methylation analysis (n = 61 peripheral blood, n = 29 non-diseased kidney, n = 99 tumors). We determine the predominant events for bilateral Wilms tumor predisposition: 1)pre-zygotic germline genetic variants readily detectable in blood DNA [WT1 (14.8%), NYNRIN (6.6%), TRIM28 (5%), and BRCA-related genes (5%)] or 2)post-zygotic epigenetic hypermethylation at 11p15.5 H19/ICR1 that may require analysis of multiple tissue types for diagnosis. Of 99 total tumor specimens, 16 (16.1%) have 11p15.5 normal retention of imprinting, 25 (25.2%) have 11p15.5 copy neutral loss of heterozygosity, and 58 (58.6%) have 11p15.5 H19/ICR1 epigenetic hypermethylation (loss of imprinting). Here, we ascertain the epigenetic and genetic modes of bilateral Wilms tumor predisposition.

Identification of intergenerational epigenetic inheritance by whole genome DNA methylation analysis in trios

Genome-wide association studies have identified thousands of loci associated with common diseases and traits. However, a large fraction of heritability remains unexplained. Epigenetic modifications, such as the observed in DNA methylation have been proposed as a mechanism of intergenerational inheritance. To investigate the potential contribution of DNA methylation to the missing heritability, we analysed the methylomes of four healthy trios (two parents and one offspring) using whole genome bisulphite sequencing. Of the 1.5 million CpGs (19%) with over 20% variability between parents in at least one family and compatible with a Mendelian inheritance pattern, only 3488 CpGs (0.2%) lacked correlation with any SNP in the genome, marking them as potential sites for intergenerational epigenetic inheritance. These markers were distributed genome-wide, with some preference to be located in promoters. They displayed a bimodal distribution, being either fully methylated or unmethylated, and were often found at the boundaries of genomic regions with high/low GC content. This analysis provides a starting point for future investigations into the missing heritability of simple and complex traits.

Comprehensive Analysis Reveals the Potential Roles of Transcription Factor Dp-1 in Lung Adenocarcinoma.

Transcription factor Dp-1 (TFDP1) was overexpressed and interacted with other genes to impact multiple signaling pathways in various human cancers. However, there is less research about the TFDP1 specific roles in lung adenocarcinoma (LUAD). We first explored TFDP1 expression levels and relative diseases from a pan-cancer perspective using the ONCOMINE, TIMER, and Open Targets Platform databases. Then, we used UALCAN, GEPIA 2, TCGA-LUAD data, and Kaplan-Meier plotter to examine TFDP1 clinicopathological features and prognosis in LUAD patients. Genomic alterations and DNA methylation analysis were performed by cBioPortal and MethSurv, respectively. Then, we used a cancer single-cell state atlas (CancerSEA) to find TFDP1 functions at a single-cell resolution. LinkedOmics was used to find TFDP1 coexpressed genes, biological processes, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Then, Gene Set Cancer Analysis (GSCA) was used to examine the drug resistence of TFDP1 in LUAD. We found that TFDP1 was overexpressed in most human cancers and related to various diseases, including LUAD. Moreover, LUAD patients with high TFDP1 expression levels might be significantly associated with individual cancer stages and have a poor prognosis. Multivariate analysis revealed that the American Joint Committee on Cancer (AJCC) pathologic stage, AJCC stage T, and AJCC stage N were the independent prognostic factors. LUAD patients with TFDP1 alterations suggested poor overall survival (OS), and disease-free survival (DFS), while hypermethylation might lead to a good prognosis. TFDP1 and its coexpressed genes were enriched in multiple signaling pathways and biological processes involved in the cell cycle, spliceosome, and DNA replication. Furthermore, TFDP1 was strongly positively related to the half-maximal inhibitory concentration (IC50) values of multiple drugs. In summary, TFDP1 was a possible biomarker and potential therapeutic target for LUAD patients.

Molecular Effects of Silicon on Arabidopsis thaliana Seedlings under UV-B Stress.

Silicon-plant interaction studies have shown that silicon reduces the harmful effects of stress in plants. Ultraviolet-B (UV-B) radiation, one of the abiotic stress affecting plants, poses a severe problem due to global warming. In this context, it is crucial to examine silicon's effects on UV-B radiation stress at the molecular level. The experiments were carried out on 17 days old Arabidopsis seedlings that were treated with 800 μWatt cm-2 doses of UV-B for 60 min and harvested on the 28th day. 1 mM orthosilicic acid was applied to the in vitro plant tissue culture for experimental groups. According to the results of the osmolyte accumulation analyses, silicon has been shown to play a role in the osmotic stress response. Gene expression levels of DGK2, CHS, FLC, RAD51, and UVR8 were measured via qPCR, and it has been shown that silicon interacts with these genes under UV-B radiation stress. The result of genomic DNA methylation analysis demonstrated that silicon might affect DNA methylation levels by increasing the 5-mC percentage compared with the control group. This study focused on the molecular effects of silicon application. It supports silicon-plant interaction research by demonstrating that silicon might affect UV-B response at the molecular level.

Abstract 363: Irf8 And Securinine: Seeking Treatments For Abdominal Aortic Aneurysm

Objective: Abdominal aortic aneurysm (AAA) is defined as local weakening and dilatation of the abdominal aorta ≥50% and has high incidence in elderly. Current therapy is limited to surgical repair, only, no pharmacological therapy has shown effectiveness against AAA development and growth. Untreated AAA rupture is highly associated with death, and aged patients are at augmented risk for surgical complications. We therefore sought to identify novel agents that could be employed for medical AAA therapy. Methods: Using whole genome DNA methylation analysis (RRBS-Seq) in a transgenerational model of murine AAA, combined with transcriptional analysis of murine AAA tissue, we identified interferon regulatory factors (IRFs: a group of immune-system-related transcription factors (TF)), as major drivers of AAA risk and pathophysiology. We performed TF enrichment analysis that specifically suggested IRF8 as being crucial for AAA development. Drug Signatures Database (DSigDB) and Connectivity Map (CMap) analysis identified securinine, an alkaloid and traditional Chinese herbal medicine derived from Securinega suffruticosa , as a probable IRF8 inhibitor. Its vascular effects are relatively unknown. Results: We studied the effects of securinine using the porcine pancreatic elastase (PPE) model of AAA. Sonographic imaging data of 28-day securinine-treated mice (15 mg/kg) showed smaller aneurysm diameters and significantly downregulated gene expression of IRF8 and IFN-γ compared to control PPE groups. Conclusion: These data suggest that securinine may have therapeutic benefits on AAA development, potentially through regulation of IRF8.

A molecular approach integrating genomic and DNA methylation profiling for tissue of origin identification in lung-specific cancer of unknown primary

BackgroundDetermining the tissue of origin (TOO) is essential for managing cancer of unknown primary (CUP). In this study, we evaluated the concordance between genome profiling and DNA methylation analysis in determining TOO for lung-specific CUP and assessed their performance by comparing the clinical responses and survival outcomes of patients predicted with multiple primary or with metastatic cancer.MethodsWe started by retrospectively screening for CUP patients who presented with both intra- and extrathoracic tumors. Tumor samples from included patients were analyzed with targeted sequencing with a 520-gene panel and targeted bisulfite sequencing. TOO inferences were made in parallel via an algorithm using genome profiles and time interval between tumors and via machine learning-based classification of DNA methylation profiles.ResultsFour hundred patients were screened retrospectively. Excluding patients definitively diagnosed with conventional diagnostic work-up or without available samples, 16 CUP patients were included. Both molecular approaches alone enabled inference of clonality for all analyzed patients. Genome profile enabled TOO inference for 43.8% (7/16) patients, and the percentage rose to 68.8% (11/16) after considering inter-tumor time lag. On the other hand, DNA methylation analysis was conclusive for TOO prediction for 100% (14/14) patients with available samples. The two approaches gave 100% (9/9) concordant inferences regarding clonality and TOO identity. Moreover, patients predicted with metastatic disease showed significantly shorter overall survival than those with multiple primary tumors.ConclusionsGenome and DNA methylation profiling have shown promise as individual analysis for TOO identification. This study demonstrated the feasibility of incorporating the two methods and proposes an integrative scheme to facilitate diagnosing and treating lung-specific CUPs.

Malignant melanotic nerve sheath tumor with PRKAR1A, KMT2C, and GNAQ mutations.

Malignant melanotic nerve sheath tumor (MMNST) is a rare and potentially aggressive lesion defined in the 2021 WHO Classification of Tumors of the Central Nervous System. MMNST demonstrate overlapping histologic and clinical features of schwannoma and melanoma. MMNST often harbor PRKAR1A mutations, especially within the Carney Complex. We present a case of aggressive MMNST of the sacral region in a 48-year-old woman. The tumor contained PRKAR1A frameshift pR352Hfs*89, KMT2C splice site c.7443-1G>T and GNAQ p.R183L missense mutations, as well as BRAF and MYC gains. Genomic DNA methylation analysis using the Illumina 850K EpicBead chip revealed that the lesion did not match an established methylation class; however, uniform manifold approximation and projection (UMAP) placed the tumor very near schwannomas. The tumor expressed PD-L1, and the patient was treated with radiation and immune checkpoint inhibitors following en bloc resection. Although she had symptomatic improvement, she suffered early disease progression with local recurrence, and distant metastases, and died 18 months after resection. It has been suggested that the presence of GNAQ mutations can differentiate leptomeningeal melanocytic neoplasms and uveal melanoma from MMNST. This case and others demonstrate that GNAQ mutations may exist in malignant nerve sheath tumors; that GNAQ and PRKAR1A mutations are not always mutually exclusive and that neither can be used to differentiate MMNST or MPNST from all melanocytic lesions.

Multiomics landscape of the autosomal dominant osteopetrosis type II disease-specific induced pluripotent stem cells

BackgroundAutosomal dominant osteopetrosis type II (ADO2) is a genetically and phenotypically metabolic bone disease, caused by osteoclast abnormalities. The pathways dysregulated in ADO2 could lead to the defects in osteoclast formation and function. However, the mechanism remains elusive.Materials and methodsTo systematically explore the molecular characterization of ADO2, we performed a multi-omics profiling from the autosomal dominant osteopetrosis type II iPSCs (ADO2-iPSCs) and healthy normal control iPSCs (NC-iPSCs) using whole genome re-sequencing, DNA methylation and N6-methyladenosine (m6A) analysis in this study.ResultsTotally, we detected 7,095,817 single nucleotide polymorphisms (SNPs) and 1,179,573 insertion and deletions (InDels), 1,001,943 differentially methylated regions (DMRs) and 2984 differential m6A peaks, and the comprehensive multi-omics profile was generated from the two cells. Interestingly, the ISG15 m6A level in ADO2-iPSCs is higher than NC-iPSCs by IGV software, and the differentially expressed m6A-modified genes (DEMGs) were highly enriched in the osteoclast differentiation and p53 signaling pathway, which associated with the development of osteopetrosis. In addition, combining our previously published transcriptome and proteome datasets, we found that the change in DNA methylation levels correlates inversely with some gene expression levels.ConclusionOur results indicate that the global multi-omics landscape not only provides a high-quality data resource but also reveals a dynamic pattern of gene expression, and found that the pathogenesis of ADO2 may begin early in life.

P04.07 The molecular evolution of oligodendrogliomas

Abstract BACKGROUND Oligodendroglioma (OD) is defined by the presence of both IDH1/2 mutation and 1p/19q codeletion. Although prognosis of OD patients is relatively favorable, tumours usually relapse and often evolve to a higher malignancy grade, with some acquiring the treatment induced hypermutated phenotype. To better understand how these tumours evolve in time, we examined the molecular differences between matched primary and recurrent ODs. MATERIAL AND METHODS We identified 21 patients who underwent surgery at least twice [male: 11, female: 10, median age: 44 years (31–66)]. Clinical data were available for 14/21 patients: 5/14 received a treatment between resections [4 radiotherapy, 1 radiotherapy followed by PCV chemotherapy]; median time from the first to the second surgery was 71.5 months (12–158). Whole genome DNA-methylation analysis was performed using Illumina’s MethylationEPIC ‘850K’ BeadChip. Results were evaluated using the Molecularneuropathology.org platform (version 3.1.5) and in R. RESULTS Most samples were WHO grade 2 ODs [14, 10 and 1 tumours in first, second and third resection group, respectively]; WHO grade 3 was found in 6, 10 and 3 tumours in first, second and third resection, respectively; in 4 patients the tumour showed malignant progression from grade 2 to 3. Most ODs exhibited an IDH1 R132H mutation [17/21 patients]; in no cases was IDH1/2 mutation lost during progression. DNA methylation analysis was successfully performed in 41/45 cases [primary OD: 17, recurrent OD: 24] for a total of 18 matched pairs. 37 samples were assigned to the “IDH mutant glioma, subclass 1p/19q codeleted OD”; the remaining 4 were assigned to various other methylation classes but CNV (copy number variation) analysis confirmed the 1p19q codeletion in all samples. Recurrent tumours exhibited de novo loss of chromosome 4 in 3/24 cases (12.5%) and loss of chromosome 13 in 3/24 cases (12.5%). In unsupervised analysis of the 1000 most variable CpG sites, samples from the same patient clustered together. This indicates that the inter-tumour variability is greater than the intra-, temporal- or grading variability between tumours. There were no overt differences in DNA methylation levels between the primary and matched recurrent OD. However, lower genome wide DNA methylation levels were observed in tumours that dedifferentiated to grade 3 ODs compared to those of grade 2, indicating that DNA demethylation is associated to higher malignancy grade. CONCLUSION DNA methylation analysis in a cohort of primary and recurrent oligodendrogliomas highlights the genomic and epi-genetic changes that are acquired at tumour progression. We are currently expanding the cohort and collecting/integrating the clinical data to better explore the evolution of recurrent ODs.

Genomic DNA methylation analysis reveals that BLNK is a key potential gene in the regulation of autophagy-related thyroid cancer progression

The purpose of this study was to explore the relationship between autophagy and DNA methylation, and to identify key genes for autophagy-regulated thyroid cancer progression. We divided patients with thyroid cancer into high-autophagy score (AS) group and low-AS group based on their AS values. The results found that AS was associated with the distant metastasis of thyroid cancer, and adversely affected prognosis. Then, we screened 359 differently expressed genes (DEGs) with DNA methylation status consistent with gene expression change. Functional classification analysis demonstrated that the 359 DEGs consistent with DNA methylation status were significantly involved in adhesion, migration and differentiation of immune cells. To further screen the key genes in the autophagy-related thyroid cancer progression, we constructed a protein-protein interactions (PPI) network and performed prognostic analysis. B cell linker (BLNK) was identified as the key potential gene affecting autophagy-related thyroid cancer progression. Finally, we verified that BLNK promoted the proliferation of thyroid cancer cells, and BLNK expression was regulated by DNA methylation. Our research provides a new perspective for exploring the relationship between autophagy and DNA methylation during the progression of thyroid cancer, and provides a new target for the treatment of metastatic thyroid cancer.

Analysis of Genomic DNA Methylation in Six Tissues of the Greater Horseshoe Bat (Rhinolophus ferrumequinum)

Analysis of Genomic DNA Methylation in Six Tissues of the Greater Horseshoe Bat (Rhinolophus ferrumequinum)

Comparative genome-wide methylation analysis of longissimus dorsi muscles in Yorkshire and Wannanhua pigs.

DNA methylation was one of the earliest discovered epigenetic modifications in vertebrates, and is an important epigenetic mechanism involved in the expression of genes in many biological processes, including muscle growth and development. Its effects on economically important traits are evidenced in reported differences in meat quality traits between Chinese indigenous pig breeds (Wannanhua pig) and Western commercial pig breeds (Yorkshire pig), and this presents a unique model for analyzing the effects of DNA methylation on these traits. In the present study, a whole genome DNA methylation analysis was performed on the two breeds using methylated DNA immunoprecipitation. GO functional enrichment and pathway enrichment analyses identified differentially methylated genes primarily associated with fatty acid metabolism, biological processes of muscle development and signaling pathways related to muscle development and pork quality. Differentially methylated genes were verified by sodium pyrosequencing, and the results were consistent with the sequencing results. The results of the integrative analysis between DNA methylation and gene expression revealed that the DNA methylation levels showed a significantly negative correlation with gene expression levels around the transcription start site of genes. In total, 41 genes were both differentially expressed and methylated; these genes were related to fat metabolism, lipid metabolism and skeletal muscle development. This study could help further explore the molecular mechanisms and phenotypic differences in pig growth and development among different breeds.

Diabetes changes gene expression but not DNA methylation in cardiac cells

BackgroundDiabetes mellitus is a worldwide epidemic that causes high mortality due to cardiovascular complications, in particular heart failure. Diabetes is associated with profound pathophysiological changes in the heart. The aim of this study was to investigate the impact of diabetes on gene expression and DNA methylation in cardiac cells. Methods and resultsTranscriptome analysis of heart tissue from mice with streptozotocin-induced diabetes revealed only 39 genes regulated, whereas cell type-specific analysis of the diabetic heart was more sensitive and more specific than heart tissue analysis and revealed a total of 3205 differentially regulated genes in five cell types. Whole genome DNA methylation analysis with basepair resolution of distinct cardiac cell types identified highly specific DNA methylation signatures of genic and regulatory regions. Interestingly, despite marked changes in gene expression, DNA methylation remained stable in streptozotocin-induced diabetes. Integrated analysis of cell type-specific gene expression enabled us to assign the particular contribution of single cell types to the pathophysiology of the diabetic heart. Finally, analysis of gene regulation revealed ligand-receptor pairs as potential mediators of heterocellular interaction in the diabetic heart, with fibroblasts and monocytes showing the highest degree of interaction. ConclusionIn summary, cell type-specific analysis reveals differentially regulated gene programs that are associated with distinct biological processes in diabetes. Interestingly, despite these changes in gene expression, cell type-specific DNA methylation signatures of genic and regulatory regions remain stable in diabetes. Analysis of heterocellular interactions in the diabetic heart suggest that the interplay between fibroblasts and monocytes is of pivotal importance.

Whole Genome DNA Methylation Analysis of Active Pulmonary Tuberculosis Disease Identifies Novel Epigenotypes: PARP9/miR-505/RASGRP4/GNG12 Gene Methylation and Clinical Phenotypes.

We hypothesized that DNA methylation patterns may contribute to the development of active pulmonary tuberculosis (TB). Illumina’s DNA methylation 450 K assay was used to identify differentially methylated loci (DML) in a discovery cohort of 12 active pulmonary TB patients and 6 healthy subjects (HS). DNA methylation levels were validated in an independent cohort of 64 TB patients and 24 HS. Microarray analysis identified 1028 DMLs in TB patients versus HS, and 3747 DMLs in TB patients after versus before anti-TB treatment, while autophagy was the most enriched signaling pathway. In the validation cohort, PARP9 and miR505 genes were hypomethylated in the TB patients versus HS, while RASGRP4 and GNG12 genes were hypermethylated, with the former two further hypomethylated in those with delayed sputum conversion, systemic symptoms, or far advanced lesions. MRPS18B and RPTOR genes were hypomethylated in TB patients with pleural involvement. RASGRP4 gene hypermethylation and RPTOR gene down-regulation were associated with high mycobacterial burden. TB patients with WIPI2/GNG12 hypermethylation or MRPS18B/FOXO3 hypomethylation had lower one-year survival. In vitro ESAT6 and CFP10 stimuli of THP-1 cells resulted in DNA de-methylation changes of the PARP9, RASGRP4, WIPI2, and FOXO3 genes. In conclusions, aberrant DNA methylation over the PARP9/miR505/RASGRP4/GNG12 genes may contribute to the development of active pulmonary TB disease and its clinical phenotypes, while aberrant DNA methylation over the WIPI2/GNG12/MARPS18B/FOXO3 genes may constitute a determinant of long-term outcomes.

Analysis of shading on DNA methylation by MSAP in Pinellia ternata

Pinellia ternata is a medicinal herb of Araceae, and its tubers are used as medicines. It is a common Chinese herbal medicine in China and has a large market demand. When exposing to strong light intensity and high temperature during the growth process, P. ternata withers in a phenomenon known as "sprout tumble", which largely limits tuber production. Shade can effectively delay sprout tumble formation and increase its yield, however the relevant regulation mechanism is unclear. DNA methylation, as a self-modifying response to environmental changes, is often involved in the regulation of plant growth and development. In this study, P. ternata grown under natural light and 90% shading were selected as the control group and the experimental group for genomic DNA methylation analysis by using methylate sensitive amplification polymorphism(MSAP). The results showed that a total of 617 loci were detected with 20 pairs of primers, of which 311 were in the natural light group and 306 in the shading group. The methylation sites in the light and shading groups accounted for 58.2% and 71.57%, respectively, and the methylation ratios in the methylation sites were 27.65% and 29.41%, respectively, indicating that shading significantly induced the genome DNA methylation of P. ternata. Compared to the natural light group, shading promoted 32.51% of the genes methylation, while inducing 16.25% gene demethylation. This study reveals the DNA methylation variation of P. ternata under shading conditions, which lays a preliminary theoretical foundation for further analysis of the mechanism of shading regulation of P. ternata growth from epigenetic level.

Abstract 839: Whole genome DNA methylation analysis of multiple myeloma identifies pervasive hypomethylation and biomarkers of survival

Abstract Multiple myeloma is a malignancy of terminally differentiated, antibody secreting B cells known as plasma cells. Normal plasma cell differentiation and cell fate are coupled to epigenetic and transcriptional reprogramming, including a proliferation-dependent global loss of DNA methylation. However, relatively little is known about the epigenetic changes that underlie myelomagenesis and how these contribute to disease etiology. To this end, we have analyzed the DNA methylome of 119 myeloma specimens from the CoMMpass study (NCT01454297) by whole genome bisulfite sequencing (WGBS) and more than 90% of these same specimens were also characterized for structural variants (long-insert whole genome sequencing) and gene expression (RNA-seq). Unsupervised hierarchical clustering grouped together specimens with t(4;14) translocations that upregulate the H3K36 dimethyltransferase NSD2 (also known as MMSET and WHSC1), which likely impacts the DNA methylation state through epigenetic cross-talk. These data also revealed a dramatic genome-wide hypomethylation where myeloma samples had a median global CpG methylation level of 41% as compared to 71% and 89% in normal plasma cells and B cells, respectively. Demethylation of the myeloma methylome occurred in large megabase domains encompassing genes that were devoid of gene expression. Conversely, DNA methylation remained mostly unchanged in the bodies of genes that were highly expressed. Although the majority of these hypomethylated domains were common across myelomas, many regions of variable methylation exist and these differences corresponded with proximal gene expression differences. These variably methylated regions were compared to PFS and OS and this identified 6,314 CpG loci where the level of DNA methylation was prognostic of outcome (P≤0.00001). These loci were clustered into discrete regions and in the majority of cases (79%), reduced DNA methylation at these loci corresponded with poor outcome. For example, several loci in the gene bodies of PRKCE, MGMT, FHIT, WWOX were prognostic of poor survival. Interestingly, myeloma t(14;16) translocations disrupt the tumor suppressor WWOX and induce the oncogene MAF. These data suggest that WWOX expression may also be lost by epigenetic mechanisms. Finally, we analyzed the DNA methylome of primary and relapsed samples for 22 patients, including multiple relapsed samples for 2 patients. These data identified genome-wide DNA methylation remodeling with a median of 1.9 million differential methylated CpGs between the newly diagnosed and relapsed specimens. These relapse differentially methylated loci coincided at the same regions in several patients and significantly overlapped loci where the DNA methylation level was prognostic of outcome. Integrative genetic, epigenetic, and transcriptional analyses for 120 myeloma samples will be presented. Citation Format: Benjamin G. Barwick, Doris R. Powell, Daniel Penaherrera, Sheri Skerget, Jonathan J. Keats, Daniel Auclair, Sagar Lonial, Lawrence H. Boise, Paula M. Vertino. Whole genome DNA methylation analysis of multiple myeloma identifies pervasive hypomethylation and biomarkers of survival [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 839.

S137 DNA METHYLATION DYNAMICS IDENTIFY LINEAGE‐SPECIFIC REGULATION PATTERNS OF HEMATOPOIETIC DIFFERENTIATION

Background:The classical hierarchical model of hematopoiesis has been revised based on single‐cell RNA sequencing (scRNAseq) data, which suggest a continuous rather than a step‐wise differentiation process. However, while differentiation trajectories can be inferred from scRNAseq data, it is impossible to discern whether discreet lineage commitment decisions occur at specific points along these trajectories. The existence of such commitment points could still be compatible with a step‐wise differentiation model and their characterization at the molecular level would be essential to accurately model the hematopoietic hierarchy.Aims:This study investigates DNA methylation programming in normal hematopoiesis to identify novel molecular commitment marks during hematopoietic differentiation. We hypothesized that whole genome DNA methylation analysis facilitates the identification of such molecular commitment marks due to the progressive nature of how this mark is programmed during differentiation.Methods:Using tagmentation‐based whole‐genome bisulfite sequencing, we generated a genome‐wide DNA methylation map of murine hematopoiesis. This map encompasses 26 murine hematopoietic cell populations throughout the hematopoietic hierarchy. Per cell population we sequenced at least three independent biological replicates. Differentially methylated regions (DMRs) were identified across all populations using DSS without smoothing. Multi‐tier single cell scRNAseq data were generated from murine bone marrow using the 10X Genomics Chromium Single Cell 3’ kit.Results:Across all cell populations studied, we identified 147,232 DMRs. Hierarchical clustering of these DMRs revealed co‐regulated epigenetic regions that show progressive DNA methylation programming during hematopoietic differentiation. These dynamically regulated regions can be interpreted as specific DNA methylation programs: pan‐hematopoietic‐, lineage‐ and cell type‐specific DNA methylation programs. The DNA methylation programs were strongly enriched in hematopoietic transcription factor binding sites as well as in lineage and cell type‐specific enhancer programs. To gain further insight into how the DNA methylation programming relates to regulation of gene expression, we generated a comprehensive multi‐tier single cell gene expression map of murine bone‐marrow cells. Integration of DNA methylome programs with single cell gene expression profiles revealed a strong anti‐correlation of promoter DNA methylation dynamics and cell‐state specific gene expression patterns. The progressive and lineage‐specific nature of DNA methylation programming occurring during hematopoietic differentiation allowed us to infer a phylogenetic tree of the hematopoietic system, which is solely based on DNA methylation dynamics. We demonstrated that previously defined subpopulations of common myeloid progenitors (CMPs) show distinct DNA methylation patterns: CD55+ CMPs showed megakaryocyte/erythrocyte patterns, while CD55‐ CMPs showed DNA methylation patterns compatible with granulocyte and monocyte differentiation. Furthermore, we showed that the identified DNA methylation programs are conserved across different mouse strains which underlines the broad applicability of this resource.Summary/Conclusion:We provide a novel view on the regulation of hematopoiesis based on DNA methylation dynamics which complements recent findings from scRNAseq studies. Our work provides a rich resource to investigate DNA methylation in normal as well as in abnormal hematopoiesis across a broad range of conditions and mouse strains.