Analysis Of Differentially Methylated Regions Research Articles

Discovery and Validation of Methylation Signatures in Circulating Cell-Free DNA for the Detection of Colorectal Cancer.

This study was conducted with the primary objective of assessing the performance of cfDNA methylation in the detection of colorectal cancer (CRC). Five tumor tissue, 20 peripheral blood leucocyte, and 169 cfDNA samples were collected for whole-genome bisulfite sequencing (WGBS) analysis. Bioinformatic analysis was conducted to identify differentially methylated regions (DMRs) and their functional characteristics. Quantitative methylation-specific PCR (qMSP) was used to validate the methylation levels of DMRs in the tissues and leucocytes. cfDNA samples from CRC patients and healthy controls were used to evaluate the performance of the DMR analysis. WGBS analysis revealed a decrease in DNA methylation levels in the CpG context in CRC tumor tissues compared with adjacent normal tissues. A total of 132 DMRs in cfDNA were identified as potential markers for diagnosing CRC. In a cohort of 95 CRC patients and 74 healthy controls, a combination of the three DMRs (DAB1, PPP2R5C, and FAM19A5) yielded an AUC of 0.763, achieving 64.21% sensitivity and 78.38% specificity in discriminating CRC patients from healthy controls. This study provides insights into DNA methylation patterns in CRC and identifies a set of DMRs in cfDNA with potential diagnostic value for CRC. These DMRs hold promise as biomarkers for CRC detection, offering promise for non-invasive CRC diagnosis. Further research is warranted to validate these findings in larger cohorts.

Extensive differential DNA methylation between tuberculosis skin test positive and skin test negative cattle

Bovine tuberculosis (bTB), caused by Mycobacterium bovis (M. bovis), represents a significant problem for the agriculture industry as well as posing a risk for human health. Current diagnostic tests for bTB target the cell-mediated immune (CMI) response to infection with M. bovis, primarily through screening of animals with the tuberculin skin test. Epigenetic modifications have been shown to alter the course of the immune response and differentially methylated regions (DMRs) might also influence the outcome of the skin test in cattle. Whole Genome Bisulphite Sequencing (WGBS) was used to profile DNA methylation levels from peripheral blood of a group of cattle identified as test positive for M. bovis (positive for the single intradermal comparative tuberculin test (SICTT) and/or the interferon-γ release assay compared to a test negative control group [n = 8/group, total of 16 WGBS libraries]. Although global methylation profiles were similar for both groups across the genome, 223 DMRs and 159 Differentially Promoter Methylated Genes (DPMGs) were identified between groups with an excess of hypermethylated sites in SICTT positive cattle (threshold > 15% differential methylation). Genes located within these DMRs included the Interleukin 1 receptor (IL1R1) and MHC related genes (BOLA and BOLA-DQB). KEGG pathway analysis identified enrichment of genes involved in Calcium and MAPK signalling, as well as metabolism pathways. Analysis of DMRs in a subset of SICTT negative cattle that were IFN-γ positive showed differential methylation of genes including Interleukin 10 Receptor, alpha (IL10RA), Interleukin 17 F (IL17F) and host defence peptides (DEFB and BDEF109). This study has identified a number of immune gene loci at which differential methylation is associated with SICTT test results and the degree of methylation could influence effective host immune responses.

Genome-wide Methylation Dynamics and Context-dependent Gene Expression Variability in Differentiating Preadipocytes.

Obesity, characterized by the accumulation of excess fat, is a complex condition resulting from the combination of genetic and epigenetic factors. Recent studies have found correspondence between DNA methylation and cell differentiation, suggesting a role of the former in cell fate determination. There is a lack of comprehensive understanding concerning the underpinnings of preadipocyte differentiation, specifically when cells are undergoing terminal differentiation (TD). To gain insight into dynamic genome-wide methylation, 3T3 L1 preadipocyte cells were differentiated by a hormone cocktail. The genomic DNA was isolated from undifferentiated cells and 4 hours, 2 days postdifferentiated cells, and 15 days TD cells. We employed whole-genome bisulfite sequencing (WGBS) to ascertain global genomic DNA methylation alterations at single base resolution as preadipocyte cells differentiate. The genome-wide distribution of DNA methylation showed similar overall patterns in pre-, post-, and terminally differentiated adipocytes, according to WGBS analysis. DNA methylation decreases at 4 hours after differentiation initiation, followed by methylation gain as cells approach TD. Studies revealed novel differentially methylated regions (DMRs) associated with adipogenesis. DMR analysis suggested that though DNA methylation is global, noticeable changes are observed at specific sites known as "hotspots." Hotspots are genomic regions rich in transcription factor (TF) binding sites and exhibit methylation-dependent TF binding. Subsequent analysis indicated hotspots as part of DMRs. The gene expression profile of key adipogenic genes in differentiating adipocytes is context-dependent, as we found a direct and inverse relationship between promoter DNA methylation and gene expression.

Abstract 7004: DNA Methylation analysis of African American colorectal cancers reveal race-specific alterations

Abstract African Americans (AAs) have among the highest incidence and mortality rates of colorectal cancer (CRC) in the US. They present with more right-sided, microsatellite stable (MSS) tumors and are diagnosed at earlier ages than non-Hispanic Whites (NHW). While DNA methylation changes and their significance have been previously described in CRC, much less is known about the unique DNA methylation changes that occur in AA CRCs compared to NHW CRCs. In the current study, we analyzed DNA methylation changes in AA patients and compared those changes to data from NHWs. AA patients diagnosed with a MSS colorectal tumor at the Yale New Health System were included in the study. A pathologist reviewed archival slides, and selected areas from CRCs (n=160) and adjacent normal tissue (n=42) were cored (median age=59; 82 females, 78 males). DNA was bisulfite converted and analysis performed using EPIC arrays. The 450K methylation array from NHW CRCs was also analyzed for 217 CRCs and 20 normal tissues obtained from TCGA (median age=68; 103 females, 114 males). Differentially methylated regions (DMRs) with |Δβ|>0.1, FDR corrected p-value <0.05 between tumor and normal tissue were identified using SeSAMe. We identified 4881 DMRs (57% hyper- and 44% hypomethylated) between AA CRC and normal tissue. Of these, 46% were in CpG Islands, 15% in CpG Shores, 4% in CpG shelves, and 36% in open sea. Most DMRs (71%) were in promoter regions, with 18% in gene bodies, 5% in enhancer regions, and 6% in intergenic regions. Next, we performed additional DMR analysis using overlapping probes from both arrays to compare NHW with AA data. We identified 2239 DMRs between AA CRC and normal tissue compared to 4138 in NHW. Of the 2239 AA DMRs, 1419 (66%) were common with NHW DMRs, 790 (34%) were unique to AAs. There were many overlapping pathways and genes; however, the DMRs impacting these pathways differed between AA and NHW CRCs. Hypermethylation of LCK, a tyrosine kinase involved in T-cell receptor signaling, was found in 60% of AA CRCs compared to 5.9% of NHW CRCs. LCK hypermethylation leads to CD4/8 T-cell inactivation and linked to poor outcomes in lung cancer. We observed AA-specific hypermethylation of SSTR1 in 57.2% of AA CRCs compared to 27.6% of NHW CRCs. Somatostasin (SST) suppresses tumor growth by regulating certain factors involved in cell proliferation. Hypermethylation of SSTR1 has been associated with various malignancies including gastric cancer. Lastly, promoter hypermethylation of 3 protocadherin genes was found in 65.1% of AAs compared to 33.1% NHWs. Protocadherins regulate WNT signaling and hypermethylation is associated with poor outcomes in multiple cancers. Taken together, these results show AA-specific CRC methylation changes compared with NHW CRCs. Further work is warranted to investigate the potential functional and clinical implications of these methylation alterations, some of which may represent therapeutic targets and biomarkers. Citation Format: Seeta Rajpara, David N. Buckley, Rosa Xicola, Reger Mikaeel, Mary K. Yagle, Baris Kerimoglu, Curtis Thorne, Dante A. Bellomo, Megha Padi, Xavier Llor, Nathan Ellis, Bodour Salhia. DNA Methylation analysis of African American colorectal cancers reveal race-specific alterations [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7004.

Triterpenoid ursolic acid regulates the environmental carcinogen benzo[a]pyrene-driven epigenetic and metabolic alterations in SKH-1 hairless mice for skin cancer interception.

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental carcinogens accountable to developing skin cancers. Recently, we reported that exposure to benzo[a]pyrene (B[a]P), a common PAH, causes epigenetic and metabolic alterations in the initiation, promotion and progression of non-melanoma skin cancer (NMSC). As a follow-up investigation, this study examines how dietary triterpenoid ursolic acid (UA) regulates B[a]P-driven epigenetic and metabolic pathways in SKH-1 hairless mice. Our results show UA intercepts against B[a]P-induced tumorigenesis at different stages of NMSC. Epigenomic cytosines followed by guanine residues (CpG) methyl-seq data showed UA diminished B[a]P-mediated differentially methylated regions (DMRs) profiles. Transcriptomic RNA-seq revealed UA revoked B[a]P-induced differentially expressed genes (DEGs) of skin cancer-related genes, such as leucine-rich repeat LGI family member 2 (Lgi2) and kallikrein-related peptidase 13 (Klk13), indicating UA plays a vital role in B[a]P-mediated gene regulation and its potential consequences in NMSC interception. Association analysis of DEGs and DMRs found that the mRNA expression of KLK13 gene was correlated with the promoter CpG methylation status in the early-stage comparison group, indicating UA could regulate the KLK13 by modulating its promoter methylation at an early stage of NMSC. The metabolomic study showed UA alters B[a]P-regulated cancer-associated metabolisms like thiamin metabolism, ascorbate and aldarate metabolism during the initiation phase; pyruvate, citrate and thiamin metabolism during the promotion phase; and beta-alanine and pathothenate coenzyme A (CoA) biosynthesis during the late progression phase. Taken together, UA reverses B[a]P-driven epigenetic, transcriptomic and metabolic reprogramming, potentially contributing to the overall cancer interception against B[a]P-mediated NMSC.

Promoter DNA methylation patterns in oral, laryngeal and oropharyngeal anatomical regions are associated with tumor differentiation, nodal involvement and survival.

Differentially methylated regions (DMRs) can be used as head and neck squamous cell carcinoma (HNSCC) diagnostic, prognostic and therapeutic targets in precision medicine workflows. DNA from 21 HNSCC and 10 healthy oral tissue samples was hybridized to a genome-wide tiling array to identify DMRs in a discovery cohort. Downstream analyses identified differences in promoter DNA methylation patterns in oral, laryngeal and oropharyngeal anatomical regions associated with tumor differentiation, nodal involvement and survival. Genome-wide DMR analysis showed 2,565 DMRs common to the three subsites. A total of 738 DMRs were unique to laryngeal cancer (n=7), 889 DMRs were unique to oral cavity cancer (n=10) and 363 DMRs were unique to pharyngeal cancer (n=6). Based on the genome-wide analysis and a Gene Ontology analysis, 10 candidate genes were selected to test for prognostic value and association with clinicopathological features. TIMP3 was associated with tumor differentiation in oral cavity cancer (P=0.039), DAPK1 was associated with nodal involvement in pharyngeal cancer (P=0.017) and PAX1 was associated with tumor differentiation in laryngeal cancer (P=0.040). A total of five candidate genes were selected, DAPK1, CDH1, PAX1, CALCA and TIMP3, for a prevalence study in a larger validation cohort: Oral cavity cancer samples (n=42), pharyngeal cancer tissues (n=25) and laryngeal cancer samples (n=52). PAX1 hypermethylation differed across HNSCC anatomic subsites (P=0.029), and was predominantly detected in laryngeal cancer. Kaplan-Meier survival analysis (P=0.043) and Cox regression analysis of overall survival (P=0.001) showed that DAPK1 methylation is associated with better prognosis in HNSCC. The findings of the present study showed that the HNSCC subsites oral cavity, pharynx and larynx display substantial differences in aberrant DNA methylation patterns, which may serve as prognostic biomarkers and therapeutic targets.

Abstract C060: DNA Methylation analysis of African American colorectal cancers reveal race-specific alterations

Abstract African Americans (AAs) have among the highest incidence and mortality rates of colorectal cancer (CRC) in the US. They present with more right-sided, microsatellite stable tumors and are diagnosed at earlier ages than non-Hispanic Whites (NHW). While DNA methylation changes and their significance have been previously described in CRC, much less is known about the unique DNA methylation changes that occur in AA CRCs compared to NHW CRCs. In the current study, we analyzed DNA methylation changes in AA patients and compared those changes to data from NHWs. AA patients diagnosed with a microsatellite-stable colorectal tumor at the Yale New Health System were included in the study. A pathologist reviewed archival slides, and selected areas from tumor (n=160) and adjacent normal tissue (n=42) were cored. The median age of patients was 59 from 82 females and 78 males. Extracted DNA was bisulfite converted and tested using EPIC arrays. DNA methylation (450K) from NHW CRCs was also analyzed for 217 CRCs and 20 normal tissues obtained from TCGA (median age=70; 114 males, 103 females). Differentially methylated regions (DMRs) with |Db|>0.1, FDR corrected p-value <0.05 between tumor and normal tissue were identified using SeSAMe. We identified 4881 DMRs (57% hyper- and 44% hypomethylated) between AA CRC and normal tissue. Of these, 46% were in CpG Islands, 15% in CpG Shores, 4% in CpG shelves, and 36% in open sea. Most DMRs (71%) were in promoter regions, compared with 18% in gene bodies, 5% in enhancer regions, and 6% in intergenic regions. Next, we performed an additional DMR analysis using overlapping probes from the 450K and EPIC arrays to compare TCGA NHW with Yale AA data. We identified 2290 DMRs between AA CRC and normal tissue compared to 4657 in NHW. Of the 2290 AA DMRs, 1405 (61%) were common with NHW DMRs, 885 (39%) were unique to AAs. There were many overlapping pathways and genes involving metabolic functions, DNA damage response, immune response, and WNT signaling. However, the DMRs impacting these pathways differed between AA and NHW CRCs. Several AA-specific DMRs were associated with immune response and included LCK and GPSM3 genes. Hypermethylation of LCK, a tyrosine kinase involved in T-cell receptor signaling, was found in 60% of AA CRCs compared to 8% of NHW CRCs. LCK hypermethylation leads to CD4/8 T-cell inactivation and is linked to poor outcomes in lung cancer. We also observed AA-specific hypermethylation of GPSM3 in 60% of AA CRCs compared to 18% of NHW CRCs. GPSM3 is associated with CD4/8 T-cell infiltration. Lastly, promoter hypermethylation of 6 protocadherin genes was found in 59% of AAs compared to 20% NHWs. Protocadherins play a role in regulating WNT signaling and hypermethylation has been associated with poor outcomes in multiple cancers. Taken together, these results show AA-specific CRC methylation changes compared with NHW CRCs. Further work is warranted to investigate the potential functional and clinical implications of these methylation alterations, some of which may represent therapeutic targets and biomarkers. Citation Format: David N. Buckley, Rosa Xicola, Mikaeel Reger, Mary K. Yagle, Baris Kerimoglu, Curtis Thorne, Dante Bellomo, Megha Padi, Xavier Llor, Nathan Ellis, Bodour Salhia. DNA Methylation analysis of African American colorectal cancers reveal race-specific alterations [abstract]. In: Proceedings of the 16th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2023 Sep 29-Oct 2;Orlando, FL. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2023;32(12 Suppl):Abstract nr C060.

CfDNA Methylation-Based Pan-Hematologic Prognostic Classification

Background: The prognosis in hematologic (heme) malignancies remains challenging due to the complex and diverse biology of these diseases and invasive procedures. Although risk stratification/prognostic methods (e.g., cytogenetics, mutations, clinical parameters) exist for single heme indications, there is a lack of prognostic tools that simultaneously address the biological heterogeneity across multiple hematologic neoplasms. This gap underscores the need for a pan-heme prognostic classifier. Here, we developed a comprehensive prognostic classifier based on targeted methylation sequencing of cfDNA in peripheral blood. Methods: Our pan-heme classifier was trained using 621 heme participants (32 Hodgkin lymphoma, 265 Non-Hodgkin Lymphomas [NHL], 128 myeloproliferative neoplasms, 72 plasma cell neoplasms, 124 monoclonal gammopathy of undetermined significance [MGUS]) and evaluated with a subset of 320 heme patients (18 Hodgkin lymphomas, 188 NHL, 70 myeloproliferative neoplasms, 44 plasma cell neoplasms) from the Circulating Cell-free Genome Atlas substudy 2 using a six-fold nested cross-validation. The top K principal components (PCs) of methylation beta values with age, sex, body mass index, smoking status, stage, and heme subtype were used to fit a Cox model. Personalized risk scores were generated for each sample by exponentiating the linear predictor. The C-index was calculated using risk scores and survival outcomes. Participants were stratified into high- and low-risk groups based on the risk scores of corresponding indications, with scores above 66.7th percentile categorized as high risk and below this threshold categorized as low risk. The risk classification was assessed by Kaplan-Meier curves and log-rank tests. To better understand methylation biology driving the classification performance, we conducted differential methylation analysis between high-risk (N=213) and low-risk (N=107) groups, using beta-binomial regression with an arcsine link function to model the region-level count data. Hypothesis testing was performed by Wald-test for each region. The differentially methylated regions (DMRs) were considered significant if q-value was less than 0.05 and the absolute value of beta-binomial model delta coefficients was greater than 0.1. The closest genes were annotated to significant DMRs, and pathway over-representation analysis was conducted using the annotated genes from significant hypermethylated DMRs or hypomethylated DMRs. Results: Methylation features are significant predictors of survival outcome even after adjusting for clinical covariates. The pan-heme risk stratification achieved an overall C-index of 0.70 (95% CI: 0.62-0.78), indicating good discrimination. The survival of the high- and low-risk groups was significantly different (log-rank p-value<0.0001). Noteworthy results were observed for specific indications: Myelodysplastic syndromes (MDS; N=16) showed a robust discriminatory power with a C-index of 0.86 (95% CI: 0.67-1.00). Hodgkin lymphoma (N=18) achieved a C-index of 0.91 (95% CI: 0.76-1.00). Diffuse large B-cell lymphoma (DLBCL; N=33) yielded a C-index of 0.71 (95% CI: 0.52-0.89). Follicular lymphoma (N=50) reached a C-index of 0.63 (95% CI: 0.39-0.86). Plasma cell neoplasm (N=44) showed a C-index of 0.65 (95% CI: 0.37-0.90). Moreover, the differential methylation analysis identified 4,110 DMRs. Overlapped features in significant DMRs compared to background were enriched in genomic regions, including CpG islands, CpG shelves, and CpG shores. Pathway enrichment analysis of hypermethylated DMRs in high-risk groups revealed significant associations with calcium, cAMP, MAPK, Ras, Rap1, Wnt, and hippo signaling pathways. Conclusions: We developed a novel, minimally-invasive, pan-heme prognostic classifier utilizing cfDNA blood-based targeted methylation sequencing technology. This cfDNA-based approach demonstrated potential for risk stratification across a diverse range of heme malignancies. Further, an approach utilizing a simple, minimally-invasive blood draw instead of a combination of multiple testing modalities, offers a pragmatic approach to disease prognosis across large heterogeneous groups of hematologic malignancies.

Acupuncture increased the Number of Retrieved Oocytes in a Mouse Model of POR: The Involvement of DNA Methylation in the Oocytes

Background:: Poor ovarian response (POR) reduces the success rate of in vitro fertilization mainly because of fewer oocytes retrieved. Acupuncture (Ac) therapy can improve the number of retrieved oocytes in the controlled ovarian stimulation program. The role of Ac in the corresponding epigenetic mechanism of POR has not been studied. Objective:: This study was conducted to determine the effect of Ac on the number of retrieved oocytes and its role in DNA methylation in a mouse model of POR. Methods:: Forty C57BL/6N female mice with normal estrous cycles were randomly classified into 4 groups of 10 each: control (Con) group, Ac-Con group, POR group, and Ac-POR group. Mice in POR and Ac-POR groups received a gastric gavage of Tripterygium wilfordii polyglycoside suspension of 50 mg/kg-1 once a day for 14 consecutive days. Ac was applied at “Shenting” (DU 24), “Guanyuan” (CV 4), “Zusanli” (ST 36), and “Shenshu” (BL 23) in the Ac-POR group for 10 min per session, once a day for 14 consecutive days. All four groups were stimulated with pregnant mare serum gonadotropin and human chorionic gonadotropin, and the number of retrieved oocytes and proportion of mature oocytes were recorded. The DNA methylation level in a single mouse oocyte in each group was analyzed using single-cell genome-wide bisulfite sequencing (scBSseq), and key pathways were identified using GO and KEGG enrichment analyses. Results:: A dissecting microscope revealed that the Ac therapy improved the number of retrieved oocytes compared with the POR group (p < 0.05). ScBS-seq showed that there was no significant change in global DNA methylation levels between the POR model and control group mice. However, differences were primarily observed in the differentially methylated regions (DMRs) of each chromosome, and Ac decreased global DNA methylation. DMR analysis identified 13 genes that may be associated with Ac treatment. Cdk5rap2 and Igf1r, which mediate germ cell apoptosis, growth, and development, maybe most closely related to the Ac treatment of POR. KEGG analysis revealed that differentially expressed genes were mainly enriched in Wnt, GnRH, and calcium signaling pathways. The genes were closely related to the regulation of POR via Ac. Conclusion:: The results suggest that DNA methylation in oocytes is related to the development of POR and that the role of Ac in affecting DNA methylation in oocytes is associated with the Wnt, GnRH, and calcium signaling pathways as well as Cdk5rap2 and Igf1r in POR mice.

Methylome-wide analysis of milk somatic cells upon subclinical mastitis in dairy cattle

Better understanding of the molecular mechanisms behind bovine mastitis is fundamental for improving the management of this disease, which continues to be of major concern for the dairy industry, especially in its subclinical form. Disease severity and progression depend on numerous aspects, such as livestock genetics, and the interaction between the causative agent, the host, and the environment. In this context, epigenetic mechanisms have proven to have a role in controlling the response of the animal to inflammation. Therefore, in this study we aimed to explore genome wide DNA methylation of milk somatic cells (SC) in healthy cows (n = 15) and cows affected by naturally occurring subclinical mastitis by Streptococcus (Strep.) agalactiae (n = 12) and Prototheca spp. (n = 11), to better understand the role of somatic cells methylome in the host response to disease. Differentially methylated regions (DMR) were evaluated comparing: i) Strep. agalactiae-infected vs. healthy; ii) Prototheca-infected vs. healthy, and iii) mastitis vs. healthy and iv) Strep. agalactiae-infected vs. Prototheca-infected. The functional analysis was performed at 2 levels. To begin with, we extracted differentially methylated genes (DMG) from promoter - DMR, which were analyzed using the Cytoscape ClueGO plug-in. Coupled with this DMG-driven approach, all the genes associated with promoter - methylated regions were fed to the Pathifier algorithm. From the DMR analysis, we identified 1,081 hypermethylated and 361 hypomethylated promoter regions in Strep. agalactiae infected animals, while 1,514 hypermethylated and 358 hypomethylated promoter regions were identified in Prototheca infected animals, when compared with the healthy controls. When considering infected animals as a whole group (regardless of the pathogen), we found 1,576 hypermethylated and 460 hypomethylated promoter regions. Both pathogens were associated with methylation differences in genes involved in pathways related to meiosis, reproduction and tissue remodeling. Exploring the whole methylome, in subclinically infected cows we observed a strong deregulation of immune related pathways, such as nuclear factor kB and toll-like receptors signaling pathways, and of energy-related pathways such as the tricarboxylic acid cycle and unsaturated fatty acid biosynthesis. In conclusion, no evident pathogen-specific SC methylome signature was detected in the present study. Overall, we observed a clear regulation of host immune response driven by DNA methylation upon subclinical mastitis. Further studies on a larger cohort of animals are needed to validate our results and to possibly identify a unique SC methylome that signifies pathogen-specific alterations.

Genome-wide Epigenetic Study of Chronic Rhinosinusitis Tissues Reveals Dysregulated Inflammatory, Immunologic and Remodeling Pathways.

Epigenetics studies mechanisms such as DNA methylation, histone modifications, non-coding RNAs, and alternative polyadenylation that can modify gene activity without changing the underlying DNA nucleotide base-pair structure. Because these changes are reversible, they have potential in developing novel therapeutics. Currently, seven pharmaceutical agents targeting epigenetic changes are FDA approved and commercially available for treatment of certain cancers. However, studies investigating epigenetics in chronic rhinosinusitis (CRS) have not been undertaken previously in the United States. The goal of this study was to investigate sinonasal DNA methylation patterns in CRS versus controls, to discern environmentally-induced epigenetic changes impacting CRS subjects. Ethmoidal samples from CRS and inferior turbinate mucosal tissue samples from controls without CRS were studied. DNA methylation was studied by Reduced Representation Bisulfite Sequencing. RADMeth® biostatistical package was used to identify differentially methylated regions (DMRs) between CRS and controls. Ingenuity Pathway analysis of DMRs was performed to identify top upstream regulators and canonical pathways. Ninety-three samples from 64 CRS subjects (36 CRSwNP; 28 CRSsNP) and 29 controls were studied. CRS and control samples differed in 13 662 CpGs sites and 1381 DMRs. Top upstream regulators identified included: 1. CRS versus controls: TGFB1, TNF, TP53, DGCR8, and beta-estradiol. 2. CRSwNP versus controls: TGFB1, CTNNB1, lipopolysaccharide, ID2, and TCF7L2. 3. CRSsNP versus controls: MYOD1, acetone, ID2, ST8SIA4, and LEPR. Differential patterns of methylation were identified between controls and CRS, CRSwNP, and CRSsNP. Epigenetic, environmentally-induced changes related to novel, inflammatory, immunologic, and remodeling pathways appear to affect epithelial integrity, cell proliferation, homeostasis, vascular permeability, and other yet uncharacterized pathways and genes.

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

ABSTRACT 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.

Abstract 6698: Combinatorial genomic and epigenomic cell-free DNA analysis of high-risk metastatic castration resistant prostate cancer reveals prognostic liquid biopsy signatures

Abstract Background: Androgen receptor signaling inhibitors (ARSI) such abiraterone and enzalutamide have significantly improved clinical outcomes in metastatic castrate-resistant prostate cancer (mCRPC) patients. However, patients with genomic alterations in the androgen receptor (AR) and its enhancer region do not respond well and acquire resistance to these inhibitors. Here, we applied our previously developed cell-free DNA (cfDNA) liquid biopsy assay (EnhanceAR-Seq) to detect these high-risk mCRPC patients prior to the administration of first-line AR-directed therapy and correlated with survival. We also interrogated the plasma methylome to identify differentially methylated regions (DMRs) across these high-risk patients. Methodology: We applied EnhanceAR-Seq to plasma cfDNA isolated from 99 mCRPC patients enrolled from two institutions (n=52 Tulane; n=47 WashU). Plasma samples were collected prior to ARSI initiation (n=63) or during treatment (n=36). We also performed Enzymatic Methyl-seq in pre-treatment plasma from 43 patients. We split these 43 patients into cell-free genomically high-risk and low-risk groups using our EnhanceAR-seq results, and conducted DMR analysis using metilene. To identify significant DMRs, we performed multiple hypothesis testing and required q<0.05 using the Benjamini-Hochberg procedure, and further required at least 5 CpGs per DMR. Results: EnhanceAR-Seq detected AR/enhancer alterations in 35% of all plasma samples. Cell-free AR/enhancer detection was highly prognostic (PFS HR=2.80, p=0.0002; OS HR=2.6, p=0.01). When considering only pre-treatment plasma, AR/enhancer alterations detected in 44% (28/63) of samples correlated significantly with worse PFS (HR=2.21, p=0.009) and OS (HR=2.60, p=0.02). AR/enhancer alterations detected in 19% (7/36) of samples collected during ARSI were also associated with worse PFS (HR=15.6, p=0.0002) and OS (HR=8.09, p=0.05). Plasma methylome analysis revealed that for cell-free genomically high-risk mCRPC patients (based on AR/enhancer alterations detected in cfDNA), significantly hypomethylated DMRs were found in the AR promoter, upstream AR enhancer, and in AR-associated genes including FOXP1 and FOLH1. Significantly hypomethylated DMRs were also observed in DNA damage repair and cell cycle genes including MSH6, MSH3, FANCD2, CDK12 and RAD51B. Hypermethylated DMRs were seen in tumor suppressor genes including ZBTB16, BRCA2, WT1 and GNAS. Conclusions: AR/enhancer alterations detected in plasma cfDNA predicted inferior survival in mCRPC patients. Cell-free genomically high-risk mCRPC patients could be distinguished from low-risk patients based on distinct methylation signatures. Citation Format: Irfan Alahi, Pradeep S. Chauhan, Alexander L. Shiang, Jace Webster, Ha X. Dang, Lilli Greiner, Breanna yang, Elisa M. Ledet, Ramandeep K. Babbra, Wenjia Feng, Peter K. Harris, Ellen B. Jaeger, Patrick J. Miller, Sydney A. Caputo, Giordano Cittolin Santos, Oliver Sartor, Russell K. Pachynski, Christopher A. Maher, Aadel A. Chaudhuri. Combinatorial genomic and epigenomic cell-free DNA analysis of high-risk metastatic castration resistant prostate cancer reveals prognostic liquid biopsy signatures [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6698.

Genome-wide analysis of DNA methylation reveals significant impact of two typical high temperature stress on biological functions in Apostichopus japonicus

Sea cucumber Apostichopus japonicus, an important marine benthic animal, is particularly sensitive to environmental temperature stress, which affected the growth, metabolism, and reproduction even survival. In this study, we used whole genome bisulfite sequencing (WGBS) to perform a comprehensive comparative analysis of the DNA methylation profiles of A. japonicus at three temperatures: normal (20 °C), aestivation (26 °C), and lethal (32 °C). The results showed that the total methylation level was highest in group T26 and lowest in group T32. Most of the DNA methylation occurred on the C bases of the CG type, and the degree of DNA methylation on the CG-type C bases was 13.32–14.99%, reflecting a moderate level of DNA methylation. Comparison of differentially methylated regions (DMRs) in the CG sequence context, we identified 17,174 DMRs in group T26_vs_T20 and 16,216 DMRs in group T32_vs_T26. There were more hypermethylated regions than hypomethylated regions, and this difference was more obvious in group T32_vs_T26. DMRs analysis of the promoter regions and the gene body regions identified more DMR-related genes in the gene body region than in the promoter region in both comparison groups (T26_vs_T20 and T32_vs_T26). GO and KEGG analysis showed that pathways such as protein processing in the endoplasmic reticulum, ubiquitin mediated proteolysis, and endocytosis, etc., were highly enriched in the DMR-related genes in the comparison groups. Some of these pathways that were significantly enriched in the hypermethylated genes in group T26_vs_T20 but were enriched in the hypomethylated genes in group T32_vs_T26. Enrichment analysis of genes both differentially methylated and differentially expressed showed that downregulated genes in group T32_vs_T26 that were also hypomethylated in the gene body region were significantly enriched in the metabolic and the antioxidant-related glutathione metabolism pathway. Our results suggested that DNA methylation develop a lot of variation in the regulation of the response of A. japonicus to two levels of heat stress: at aestivation temperature and at lethal temperature, and there are many differences between the two. The results of this study provide important insights into the epigenetic mechanisms underlying the response of A. japonicus to the pressure of environmental temperature.

A Race-Specific, DNA Methylation Analysis of Aging in Normal Rectum: Implications for the Biology of Aging and Its Relationship to Rectal Cancer.

Approximately 90% of colorectal cancer (CRC) develop over the age of 50, highlighting the important role of aging in CRC risk. African Americans (AAs) shoulder a greater CRC burden than European Americans (EA) and are more likely to develop CRC at a younger age. The effects of aging in AA and EA normal rectal tissue have yet to be defined. Here, we performed epigenome-wide DNA methylation analysis in the first, large-scale biracial cohort of normal rectum (n = 140 samples). We identified increased epigenetic age acceleration in EA than AA rectum (p = 3.91 × 10-4) using linear regression. We also identified differentially methylated regions (DMRs) associated with chronological aging in AA and EA, separately using DMRcate. Next, a consensus set of regions associated with cancer was identified through DMR analysis of two rectal cancer cohorts. The vast majority of AA DMRs were present in our analysis of aging in rectum of EA subjects, though rates of epigenetic drift were significantly greater in AA (p = 1.94 × 10-45). However, 3.66-fold more DMRs were associated with aging in rectum of EA subjects, many of which were also associated with rectal cancer. Our findings reveal a novel relationship between race, age, DNA methylation and rectal cancer risk that warrants further investigation.

Epigenome‐wide association study and integrative cross‐omics analyses of cerebral white matter hyperintensities implicate novel gene networks related to Alzheimer’s disease

AbstractBackgroundCerebral white matter hyperintensities (WMH) on magnetic resonance imaging are markers of cerebral small vessel disease (cSVD), a major risk factor for Alzheimer’s disease and related dementia (ADRD). Despite the successful identification of multiple genetic variants associated with this highly heritable condition, its genetic architecture remains incompletely understood. More specifically, the role of DNA methylation (DNAm) has received little attention.MethodWe investigated the association of blood DNAm at approximately 450,000 CpG sites with WMH burden in 9,732 middle to late age adults from 14 cohorts. Integrative cross‐omics analyses followed to elucidate the role of identified DNAm in cSVD.ResultWe identified 12 single‐CpG and 46 region‐based DNAm associations with WMH burden. Our top discovery single CpG, cg24202936 (P = 7.6×10−8), was associated with F2 expression in blood (P = 6.4×10−5), and colocalized with FOLH1 expression in brain (posterior probability = 0.75). Some single‐CpG loci (cg17417856 and cg06809326) with suggestive evidence (P<1×10−5) were also identified as differentially methylated regions (DMRs) in PRMT1 and in CCDC144NL‐AS1. Mendelian randomization analyses showed that cg06809326 is causally associated with WMH burden (OR [95% CI]: 1.39 [1.03‐1.87], P = 0.03) and that expression of CCDC144NL‐AS1 mediates this association. DMR analysis, joint epigenetic association analysis, and multi‐omics colocalization analysis consistently identified a role of DNAm near SH3PXD2A, a locus previously identified in genome‐wide association studies of WMH. Gene set enrichment analyses revealed functions of the WMH‐associated DNAm loci in the blood‐brain barrier (BBB) and in the immune response. Integrative cross‐omics analysis identified 19 key regulatory genes in two WMH‐associated networks related to extracellular matrix (ECM) organization, and lipid and lipoprotein metabolism. A drug repositioning analysis indicated peroxisome proliferator‐activated receptor alpha (PPAR‐α) agonist as a target drug for WMH, which was also proposed as a promising target for ADRD.ConclusionOur epigenome‐wide association study and integrative cross‐omics analyses implicate novel genes influencing WMH burden, which converged on pathways related to the immune response and to a compromised BBB possibly due to disrupted cell‐cell and cell‐ECM interactions. The results also suggest that PPAR‐α may contribute to lowering risk for WMH possibly through protection against BBB disruption, which may potentially modify the ADRD risk.

PTEN-knockout regulates metabolic rewiring and epigenetic reprogramming in prostate cancer and chemoprevention by triterpenoid ursolic acid.

PTEN (phosphatase and tensin homolog deleted on chromosome 10) is one of the most frequently mutated/deleted tumor suppressor genes in many human cancers. Ursolic acid (UA) is a natural triterpenoid possessing antioxidant, anti-inflammatory, and anticancer effects. However, how PTEN impacts metabolic rewiring and how UA modifies PTEN-driven metabolic and epigenetic reprogramming in prostate cancer (PCa) remains unknown. In the current study, we found that UA protects against PTEN knockout (KO)-induced tumorigenesis at different stages of PCa. Epigenomic CpG methyl-seq revealed UA attenuated PTEN KO-induced differentially methylated regions (DMRs) profiles. Transcriptomic RNA-seq showed UA abrogated PTEN KO-induced differentially expressed genes (DEGs) of PCa-related oncogenes' Has3, Cfh, and Msx1 overexpression, indicating UA plays a crucial role in PTEN KO-mediated gene regulation and its potential consequences on cancer interception. Association analysis of DEGs and DMRs identified that the mRNA expression of tumor suppressor gene BDH2, and oncogenes Ephas, Isg15, and Nos2 were correlated with the promoter CpG methylation status in the early-stage comparison groups indicating UA could regulate the oncogenes or tumor suppressor genes by modulating their promoter methylation at an early stage of prostate tumorigenesis. The metabolomic study showed UA attenuated PTEN KO-regulated cancer-associated metabolisms like purine metabolism/metabolites correlating with RNAseq findings, glycolysis/gluconeogenesis metabolism, as well as epigenetic-related metabolites pyruvate and lactate indicating UA plays a critical role in PTEN KO-mediated metabolic and epigenetic reprogramming and its consequences on cancer development. In this context, UA impacts metabolic rewiring causing epigenetic and transcriptomic reprogramming potentially contributing to the overall protection against prostate-specific PTEN KO-mediated PCa.

Long‐term myocardial DNA methylation and metabolomic responses to exertional heat stroke in mice

Epidemiological studies indicate that exposure to exertional heat stroke (EHS) can increase the risk of all‐cause cardiovascular and systemic disease and promote early mortality. Using a preclinical mouse model of EHS, we previously characterized profound DNA methylation changes in immune cells at 30 d of recovery and severe cardiometabolic disturbances in the ventricular myocardium (VM) at 9‐14 d. These observations indicate that a single exposure to EHS can result in long‐term epigenetic and metabolic alterations. However, whether a long‐term epigenetic responses occur in the VM and or whether cardiometabolic perturbations persist is unknown.PURPOSETo determine if long‐term alterations to (1) DNA methylation and (2) the metabolome occur within the VM after 30 d of recovery from EHS and (3) whether epigenetic alterations predict imbalances within the metabolome.METHODSSixteen C57BL6 female mice underwent a standardized EHS protocol using a forced running wheel (environmental temp: 37.5°C, 40% humidity) or a matched exercise control trial (EXC) (22.5°C). The EHS mice achieved peak core temps of ~42.4°C, accompanied by transient loss of consciousness. Mice were sacrificed after 30 d of recovery; left ventricles were isolated and frozen in liquid nitrogen for DNA methylation sequencing and untargeted metabolomics.RESULTSFollowing 30 d of recovery from EHS, more than 6,000 differentially methylated cytosines (DMCs) were observed within promoter regions; 900 differentially methylated regions (DMRs) were also observed (i.e ≥ 5 differentially DMCs within 300 base pairs) at (Q < 0.05) compared to EXC. Pathway analysis of promoter DMRs indicated changes in genes that modulate basic cell functions, DNA binding, transcription, and metabolism. Cardiometabolic perturbations were found to largely resolve compared to observations at 9‐14 d. Some metabolic disturbances remained, e.g. reductions in α‐ketoglutarate and increases in several lipid intermediates. These remnant alterations are known to be common in early onset heart failure and chronic kidney disease. Repressive DMCs (P < 0.05) were observed within transcriptional regulatory regions of enzymes that regulate α‐ketoglutarate, e.g. isocitrate dehydrogenase 1 & 2 (IDH1 & IDH2), suggesting a potential link with metabolomic findings.CONCLUSIONSA single EHS episode results in a robust epigenetic response in VM when observed at 30 d of recovery, whereas all but a subset of metabolites within the VM had fully recovered. These long‐term molecular changes may induce altered or maladaptive physiological outcomes in response to future environmental or medical challenges.

Exposure to 3,3',4,4',5-Pentachlorobiphenyl (PCB126) Causes Widespread DNA Hypomethylation in Adult Zebrafish Testis.

Exposure to environmental toxicants during preconception has been shown to affect offspring health and epigenetic mechanisms such as DNA methylation are hypothesized to be involved in adverse outcomes. However, studies addressing the effects of exposure to environmental toxicants during preconception on epigenetic changes in gametes are limited. The objective of this study is to determine the effect of preconceptional exposure to a dioxin-like polychlorinated biphenyl (3,3',4,4',5-pentachlorobiphenyl [PCB126]) on DNA methylation and gene expression in testis. Adult zebrafish were exposed to 3 and 10 nM PCB126 for 24 h and testis tissue was sampled at 7 days postexposure for histology, DNA methylation, and gene expression profiling. Reduced representation bisulfite sequencing revealed 37 and 92 differentially methylated regions (DMRs) in response to 3 and 10 nM PCB126 exposures, respectively. Among them, 19 DMRs were found to be common between both PCB126 treatment groups. Gene ontology (GO) analysis of DMRs revealed that enrichment of terms such as RNA processing, iron-sulfur cluster assembly, and gluconeogenesis. Gene expression profiling showed differential expression of 40 and 1621 genes in response to 3 and 10 nM PCB126 exposures, respectively. GO analysis of differentially expressed genes revealed enrichment of terms related to xenobiotic metabolism, oxidative stress, and immune function. There is no overlap in the GO terms or individual genes between DNA methylation and RNA sequencing results, but functionally many of the altered pathways have been shown to cause spermatogenic defects.

Novel Zinc-Related Differentially Methylated Regions in Leukocytes of Women With and Without Obesity.

IntroductionNutriepigenetic markers are predictive responses associated with changes in “surrounding” environmental conditions of humans, which may influence metabolic diseases. Although rich in calories, Western diets could be linked with the deficiency of micronutrients, resulting in the downstream of epigenetic and metabolic effects and consequently in obesity. Zinc (Zn) is an essential nutrient associated with distinct biological roles in human health. Despite the importance of Zn in metabolic processes, little is known about the relationship between Zn and epigenetic. Thus, the present study aimed to identify the epigenetic variables associated with Zn daily ingestion (ZnDI) and serum Zinc (ZnS) levels in women with and without obesity.Materials and MethodsThis is a case-control, non-randomized, single-center study conducted with 21 women allocated into two groups: control group (CG), composed of 11 women without obesity, and study group (SG), composed of 10 women with obesity. Anthropometric measurements, ZnDI, and ZnS levels were evaluated. Also, leukocyte DNA was extracted for DNA methylation analysis using 450 k Illumina BeadChips. The epigenetic clock was calculated by Horvath method. The chip analysis methylation pipeline (ChAMP) package selected the differentially methylated regions (DMRs).ResultsThe SG had lower ZnS levels than the CG. Moreover, in SG, the ZnS levels were negatively associated with the epigenetic age acceleration. The DMR analysis revealed 37 DMRs associated with ZnDI and ZnS levels. The DMR of PM20D1 gene was commonly associated with ZnDI and ZnS levels and was hypomethylated in the SG.ConclusionOur findings provide new information on Zn's modulation of DNA methylation patterns and bring new perspectives for understanding the nutriepigenetic mechanisms in obesity.