Promoter DNA Hypermethylation and Paradoxical Gene Activation

DEP domain containing 1 as a biomarker for poor prognosis in lung adenocarcinoma

Objective

The DEP domain-containing 1 (DEPDC1) gene is essential in the development and advancement of different types of cancer. This study is to examine the levels of DEPDC1 in lung adenocarcinoma (LUAD), and to determine its relationship with clinical results and immune response. The goal is to assess its potential as a biomarker and therapeutic target for LUAD.

Methods

By comprehensively utilizing the Cancer Genome Atlas (TCGA), gene Expression Synthesis (GEO), UALCAN, cBioPortal, TISIDB databases and online platforms, we conducted a bioinformatics analysis to investigate DEPDC1 gene survival analysis, prognostic diagnosis, prognostic survival, immune cell infiltration, DNA methylation, and the correlation of genetic mutations in LUAD. The results were validated through cell assay and immunohistochemical staining.

Results

DEPDC1 shows high levels of expression in the majority of tumors, with its expression being notably elevated in LUAD compablue to normal tissues. The expression of DEPDC1 varies based on the clinical characteristics of patients with LUAD. DEPDC1 expression affects the survival prognosis and prognostic model construction of LUAD patients. In addition, the presence of DEPDC1 is linked to immune infiltration. Various chemokines and chemokine receptors, immunoinhibitors and immune-stimulators in LUAD are significantly correlated with DEPDC1 methylation levels. Cell experiments confirmed through qPCR that the mRNA expression of DEPDC1 in LUAD was markedly elevated in comparison to the normal population, and immunohistochemistry showed positive DEPDC1 expression in LUAD pathological sections.

Conclusion

Systematic analysis and experiments have verified that DEPDC1 serves as a biomarker for detecting early, prediction of survival, and evaluation of immune cell infiltration in LUAD.

The cuproptosis-related gene UBE2D2 functions as an immunotherapeutic and prognostic biomarker in pan-cancer

BACKGROUND

Cuproptosis, as a unique modality of regulated cell death, requires the involvement of ubiquitin-binding enzyme UBE2D2. However, the prognostic and immunotherapeutic values of UBE2D2 in pan-cancer remain largely unknown.

METHODS

Using UCSC Xena, TIMER, Clinical Proteomic Tumor Analysis Consortium (CPTAC), and Human Protein Atlas (HPA) databases, we aimed to explore the differential expression pattern of UBE2D2 across multiple cancer types and to evaluate its association with patient prognosis, clinical features, and genetic variations. The association between UBE2D2 and immunotherapy response was assessed by gene set enrichment analysis, tumor microenvironment, immune gene co-expression and drug half maximal inhibitory concentration (IC50) analysis.

RESULTS

The mRNA and protein levels of UBE2D2 were markedly elevated in most cancer types, and UBE2D2 exhibited prognostic significance in liver hepatocellular carcinoma (LIHC), kidney chromophobe (KICH), uveal melanomas (UVM), cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), and kidney renal papillary cell carcinoma (KIRP). UBE2D2 expression was correlated with clinical features, tumor mutation burden, microsatellite instability, and anti-tumor drug resistance in several tumor types. Gene enrichment analysis showed that UBE2D2 was significantly associated with immune-related pathways. The expression level of UBE2D2 was correlated with immune cell infiltration, including CD4 + T cells、Macrophages M2、CD8 + T cells in pan-cancer. PDCD1, CD274 and CTLA4 expression levels were positively correlated with UBE2D2 level in multiple cancers.

CONCLUSIONS

We comprehensively investigated the potential value of UBE2D2 as a prognostic and immunotherapeutic predictor for pan-cancer, providing a novel insight for cancer immunotherapy.

Molecular adaptations in response to exercise training are associated with tissue-specific transcriptomic and epigenomic signatures

Regular exercise has many physical and brain health benefits, yet the molecular mechanisms mediating exercise effects across tissues remain poorly understood. Here we analyzed 400 high-quality DNA methylation, ATAC-seq, and RNA-seq datasets from eight tissues from control and endurance exercise-trained (EET) rats. Integration of baseline datasets mapped the gene location dependence of epigenetic control features and identified differing regulatory landscapes in each tissue. The transcriptional responses to 8 weeks of EET showed little overlap across tissues and predominantly comprised tissue-type enriched genes. We identified sex differences in the transcriptomic and epigenomic changes induced by EET. However, the sex-biased gene responses were linked to shared signaling pathways. We found that many G protein-coupled receptor-encoding genes are regulated by EET, suggesting a role for these receptors in mediating the molecular adaptations to training across tissues. Our findings provide new insights into the mechanisms underlying EET-induced health benefits across organs.

Whole-genome bisulfite sequencing identifies stage- and subtype-specific DNA methylation signatures in pancreatic cancer

In pancreatic ductal adenocarcinoma (PDAC), no recurrent metastasis-specific mutation has been found, suggesting that epigenetic mechanisms, such as DNA methylation, are the major contributors of late-stage disease progression. Here, we performed the first whole-genome bisulfite sequencing (WGBS) on mouse and human PDAC organoid models to identify stage-specific and molecular subtype-specific DNA methylation signatures. With this approach, we identified thousands of differentially methylated regions (DMRs) that can distinguish between the stages and molecular subtypes of PDAC. Stage-specific DMRs are associated with genes related to nervous system development and cell-cell adhesions, and are enriched in promoters and bivalent enhancers. Subtype-specific DMRs showed hypermethylation of GATA6 foregut endoderm transcriptional networks in the squamous subtype and hypermethylation of EMT transcriptional networks in the progenitor subtype. These results indicate that aberrant DNA methylation contributes to both PDAC progression and subtype differentiation, resulting in significant and reoccurring DNA methylation patterns with diagnostic and prognostic potential.

Nitric oxide inhibits ten-eleven translocation DNA demethylases to regulate 5mC and 5hmC across the genome

DNA methylation at cytosine bases of eukaryotic DNA (5-methylcytosine, 5mC) is a heritable epigenetic mark that can regulate gene expression in health and disease. Enzymes that metabolize 5mC have been well-characterized, yet the discovery of endogenously produced signaling molecules that regulate DNA methyl-modifying machinery have not been described. Herein, we report that the free radical signaling molecule nitric oxide (NO) can directly inhibit the Fe(II)/2-OG-dependent DNA demethylases ten-eleven translocation (TET) and human AlkB homolog 2 (ALKBH2). Physiologic NO concentrations reversibly inhibited TET and ALKBH2 demethylase activity by binding to the mononuclear non-heme iron atom which formed a dinitrosyliron complex (DNIC) preventing cosubstrates (2-OG and O2) from binding. In cancer cells treated with exogenous NO, or cells endogenously synthesizing NO, there was a global increase in 5mC and 5-hydroxymethylcytosine (5hmC) in DNA, the substrates for TET, that could not be attributed to increased DNA methyltransferase activity. 5mC was also elevated in NO-producing cell-line-derived mouse xenograft and patient-derived xenograft tumors. Genome-wide DNA methylome analysis of cells chronically treated with NO (10 days) demonstrated enrichment of 5mC and 5hmC at gene-regulatory loci which correlated to changes in the expression of NO-regulated tumor-associated genes. Regulation of DNA methylation is distinctly different from canonical NO signaling and represents a novel epigenetic role for NO.

PTEN regulates expression of its pseudogene in glioblastoma cells in DNA methylation-dependent manner

Glioblastoma (GBM) is the most aggressive and frequent type of primary brain cancer in adult patients. One of the key molecular features associated with GBM pathogenesis is the dysfunction of PTEN oncosuppressor. In addition to PTEN gene, humans and several primates possess processed PTEN pseudogene (PTENP1) that gives rise to long non-coding RNA lncPTENP1-S. Regulation and functions of PTEN and PTENP1 are highly interconnected, however, the exact molecular mechanism of how these two genes affect each other remains unclear. Here, we analyzed the methylation level of the CpG islands (CpGIs) in the promoter regions of PTEN and PTENP1 in patient-derived GBM neurospheres. We found that increased PTEN methylation corelates with decreased PTEN mRNA level. Unexpectedly, we showed the opposite trend for PTENP1. Using targeted methylation and demethylation of PTENP1 CpGI, we demonstrated that DNA methylation increases lncPTENP1-S expression in the presence of wild type PTEN protein but decreases lncPTENP1-S expression if PTEN protein is absent. Further experiments revealed that PTEN protein binds to PTENP1 promoter region and inhibits lncPTENP1-S expression if its CpGI is demethylated. Interestingly, we did not detect any effect of lncPTENP1-S on the level of PTEN mRNA, indicating that in GBM cells PTENP1 is a downstream target of PTEN rather than its upstream regulator. Finally, we studied the functions of lncPTENP1-S and demonstrated that it plays a pro-oncogenic role in GBM cells by upregulating the expression of cancer stem cell markers and decreasing cell adhesion.

DNA methylation-altered genes in the rat hippocampal neurogenic niche after continuous exposure to amorphous curcumin

Rat offspring who are exposed to an amorphous formula of curcumin (CUR) from the embryonic stage have anti-anxiety-like behaviors, enhanced fear extinction learning, and increased synaptic plasticity in the hippocampal dentate gyrus (DG). In the present study, we investigated the links between genes with altered methylation status in the neurogenic niche and enhanced neural functions after CUR exposure. We conducted methylation and RNA sequencing analyses of the DG of CUR-exposed rat offspring on day 77 after delivery. Methylation status and transcript levels of candidate genes were validated using methylation-sensitive high-resolution melting and real-time reverse-transcription PCR, respectively. In the CUR group, we confirmed the hypermethylation and downregulation of Gpr150, Mmp23, Rprml, and Pcdh8 as well as the hypomethylation and upregulation of Ppm1j, Fam222a, and Opn3. Immunohistochemically, reprimo-like+ hilar cells and protocadherin-8+ granule cells were decreased and opsin-3+ hilar cells were increased by CUR exposure. Both reprimo-like and opsin-3 were partially expressed on subpopulations of glutamic acid decarboxylase 67+ γ-aminobutyric acid-ergic interneurons. Furthermore, the transcript levels of genes involved in protocadherin-8-mediated N-cadherin endocytosis were altered with CUR exposure; this was accompanied by Ctnnb1 and Syp upregulation and Mapk14, Map2k3, and Grip1 downregulation, suggesting that CUR-induced enhanced synaptic plasticity is associated with cell adhesion. Together, our results indicate that functionally different genes have altered methylation and expression in different neuronal populations of the hippocampal neurogenic niche, thus enhancing synaptic plasticity after CUR exposure.

Pan-cancer analysis of NUP155 and validation of its role in breast cancer cell proliferation, migration, and apoptosis

NUP155 is reported to be correlated with tumor development. However, the role of NUP155 in tumor physiology and the tumor immune microenvironment (TIME) has not been previously examined. This study comprehensively investigated the expression, immunological function, and prognostic significance of NUP155 in different cancer types. Bioinformatics analysis revealed that NUP155 was upregulated in 26 types of cancer. Additionally, NUP155 upregulation was strongly correlated with advanced pathological or clinical stages and poor prognosis in several cancers. Furthermore, NUP155 was significantly and positively correlated with DNA methylation, tumor mutational burden, microsatellite instability, and stemness score in most cancers. Additionally, NUP155 was also found to be involved in TIME and closely associated with tumor infiltrating immune cells and immunoregulation-related genes. Functional enrichment analysis revealed a strong correlation between NUP155 and immunomodulatory pathways, especially antigen processing and presentation. The role of NUP155 in breast cancer has not been examined. This study, for the first time, demonstrated that NUP155 was upregulated in breast invasive carcinoma (BRCA) cells and revealed its oncogenic role in BRCA using molecular biology experiments. Thus, our study highlights the potential value of NUP155 as a biomarker in the assessment of prognostic prediction, tumor microenvironment and immunotherapeutic response in pan-cancer.

CDK12 is a potential biomarker for diagnosis, prognosis and immunomodulation in pan-cancer

Cell cycle-dependent protein kinase 12 (CDK12) plays a key role in a variety of carcinogenesis processes and represents a promising therapeutic target for cancer treatment. However, to date, there have been no systematic studies addressing its diagnostic, prognostic and immunological value across cancers. Here, we found that CDK12 was significantly upregulated in various types of cancers, and it expression increased with progression in ten cancer types, including breast cancer, cholangiocarcinoma and colon adenocarcinoma. Moreover, the ROC curves indicated that CDK12 showed diagnostic value in eight cancer types. High CDK12 expression was associated with poor prognosis in eight types of cancer, including low-grade glioma, mesothelioma, melanoma and pancreatic cancer. Furthermore, we conducted immunoassays to explore the exact mechanisms underlying CDK12-induced carcinogenesis, which revealed that increased expression of CDK12 allowed tumours to evade immune surveillance and upregulate immune checkpoint genes. Additionally, mutational studies have shown that amplification and missense mutations are the predominant mutational events affecting CDK12 across cancers. These findings establish CDK12 as a significant biological indicator of cancer diagnosis, prognosis, and immunotherapeutic targeting. Early surveillance and employment of CDK12 inhibitors, along with concomitant immunotherapy interventions, may enhance the clinical outcomes of cancer patients.

Development of machine learning-based malignant pericardial effusion-related model in breast cancer: Implications for clinical significance, tumor immune and drug-therapy

Background

Malignant pericardial effusion (MPE) is a common complication of advanced breast cancer (BRCA) and plays an important role in BRCA. This study is aims to construct a prognostic model based on MPE-related genes for predicting the prognosis of breast cancer.

Methods

The BRCA samples are analyzed based on the expression of MPE-related genes by using an unsupervised cluster analysis method. This study processes the data by least absolute shrinkage and selection operator and multivariate Cox analysis, and uses machine learning algorithms to construct BRCA prognostic model and develop web tool.

Results

BRCA patients are classified into three clusters and a BRCA prognostic model is constructed containing 9 MPE-related genes. There are significant differences in signature pathways, immune infiltration, immunotherapy response and drug sensitivity testing between the high and low-risk groups. Of note, a web-based tool (http://wys.helyly.top/cox.html) is developed to predict overall survival as well as drug-therapy response of BRCA patients quickly and conveniently, which can provide a basis for clinicians to formulate individualized treatment plans.

Conclusion

The MPE-related prognostic model developed in this study can be used as an effective tool for predicting the prognosis of BRCA and provides new insights for the diagnosis and treatment of BRCA patients.

Regulation and tumor-suppressive function of the miR-379/miR-656 (C14MC) cluster in cervical cancer

Cervical cancer (CC) is a key contributor to cancer-related mortality in several countries. The identification of molecular markers and the underlying mechanism may help improve CC management. We studied the regulation and biological function of the chromosome 14 microRNA cluster (C14MC; miR-379/miR-656) in CC. Most C14MC members exhibited considerably lower expression in CC tissues and cell lines in The Cancer Genome Atlas (TCGA) cervical squamous cell carcinoma and endocervical adenocarcinoma patient cohorts. Bisulfite Sanger sequencing revealed hypermethylation of the C14MC promoter in CC tissues and cell lines. 5-aza-2 deoxy cytidine treatment reactivated expression of the C14MC members. We demonstrated that C14MC is a methylation-regulated miRNA cluster via artificial methylation and luciferase reporter assays. C14MC downregulation correlated with poor overall survival and may promote metastasis. C14MC activation via the lentiviral-based CRISPRa approach inhibited growth, proliferation, migration, and invasion; enhanced G2/M arrest; and induced senescence. Post-transcriptional regulatory network analysis of C14MC transcriptomic data revealed enrichment of key cancer-related pathways, such as metabolism, the cell cycle, and phosphatidylinositol 3-kinase (PI3K)-AKT signaling. Reduced cell proliferation, growth, migration, invasion, and senescence correlated with the downregulation of active AKT, MYC, and cyclin E1 (CCNE1) and the overexpression of p16, p21, and p27. We showed that C14MC miRNA activation increases reactive oxygen species (ROS) levels, intracellular Ca2+ levels, and lipid peroxidation rates, and inhibits epithelial-mesenchymal transition (EMT). C14MC targets pyruvate dehydrogenase kinase-3 (PDK3) according to the luciferase reporter assay. PDK3 is overexpressed in CC and is inversely correlated with C14MC. Both miR-494-mimic transfection and C14MC activation inhibited PDK3 expression. Reduced glucose uptake and lactate production, and upregulation of PDK3 upon C14MC activation suggest the potential role of these proteins in metabolic reprogramming. Finally, we showed that C14MC activation may inhibit EMT signaling. Thus, C14MC is a tumor-suppressive and methylation-regulated miRNA cluster in CC. Reactivation of C14MC can be useful in the management of CC.

Targeted DNA Methylation Editing Using an All-in-One System Establishes Paradoxical Activation of EBF3

Cutaneous melanoma is rapidly on the rise globally, surpassing the growth rate of other cancers, with metastasis being the primary cause of death in melanoma patients. Consequently, understanding the mechanisms behind this metastatic process and exploring innovative treatments is of paramount importance. Recent research has shown promise in unravelling the role of epigenetic factors in melanoma progression to metastasis. While DNA hypermethylation at gene promoters typically suppresses gene expression, we have contributed to establishing the newly understood mechanism of paradoxical activation of genes via DNA methylation, where high methylation coincides with increased gene activity. This mechanism challenges the conventional paradigm that promoter methylation solely silences genes, suggesting that, for specific genes, it might actually activate them. Traditionally, altering DNA methylation in vitro has involved using global demethylating agents, which is insufficient for studying the mechanism and testing the direct consequence of gene methylation changes. To investigate promoter hypermethylation and its association with gene activation, we employed a novel approach utilising a CRISPR-SunTag All-in-one system. Here, we focused on editing the DNA methylation of a specific gene promoter segment (EBF3) in melanoma cells using the All-in-one system. Using bisulfite sequencing and qPCR with RNA-Seq, we successfully demonstrated highly effective methylation and demethylation of the EBF3 promoter, with subsequent gene expression changes, to establish and validate the paradoxical role of DNA methylation. Further, our study provides novel insights into the function of the EBF3 gene, which remains largely unknown. Overall, this study challenges the conventional view of methylation as solely a gene-silencing mechanism and demonstrates a potential function of EBF3 in IFN pathway signalling, potentially uncovering new insights into epigenetic drivers of malignancy and metastasis.

Pan-cancer analysis identifies the IRF family as a biomarker for survival prognosis and immunotherapy

IRF family genes have been shown to be crucial in tumorigenesis and tumour immunity. However, information about the role of IRF in the systematic assessment of pan-cancer and in predicting the efficacy of tumour therapy is still unknown. In this work, we performed a systematic analysis of IRF family genes in 33 tumour samples, including expression profiles, genomics and clinical characteristics. We then applied Single-Sample Gene-Set Enrichment Analysis (ssGSEA) to calculate IRF-scores and analysed the impact of IRF-scores on tumour progression, immune infiltration and treatment efficacy. Our results showed that genomic alterations, including SNPs, CNVs and DNA methylation, can lead to dysregulation of IRFs expression in tumours and participate in regulating multiple tumorigenesis. IRF-score expression differed significantly between 12 normal and tumour samples and the impact on tumour prognosis and immune infiltration depended on tumour type. IRF expression was correlated to drug sensitivity and to the expression of immune checkpoints and immune cell infiltration, suggesting that dysregulation of IRF family expression may be a critical factor affecting tumour drug response. Our study comprehensively characterizes the genomic and clinical profile of IRFs in pan-cancer and highlights their reliability and potential value as predictive markers of oncology drug efficacy. This may provide new ideas for future personalized oncology treatment.

Unveiling the potential of HSPA4: a comprehensive pan-cancer analysis of HSPA4 in diagnosis, prognosis, and immunotherapy

With the global rise in cancer incidence and mortality rates, research on the topic has become increasingly urgent. Among the significant players in this field are heat shock proteins (HSPs), particularly HSPA4 from the HSP70 subfamily, which has recently garnered considerable interest for its role in cancer progression. However, despite numerous studies on HSPA4 in specific cancer types, a comprehensive analysis across all cancer types is lacking. This study employs various bioinformatics techniques to delve into the role of HSPA4 in pan-cancer. Our objective is to assess its potential in clinical diagnosis, prognosis, and as a future molecular target for therapy. The research findings reveal significant differences in HSPA4 expression across different cancer types, suggesting its diagnostic value and close association with cancer staging and patient survival rates. Furthermore, genetic variations and methylation status of HSPA4 play critical roles in tumorigenesis. Lastly, the interaction of HSPA4 with immune cells is linked to the tumor microenvironment (TME) and immunotherapy. In summary, HSPA4 emerges as a promising cancer biomarker and a vital member of the HSPs family, holding potential applications in diagnosis, prognosis, and immunotherapy.

DNMT1/DNMT3a-mediated promoter hypermethylation and transcription activation of ICAM5 augments thyroid carcinoma progression

Promoter methylation is one of the most studied epigenetic modifications and it is highly relevant to the onset and progression of thyroid carcinoma (THCA). This study investigates the promoter methylation and expression pattern of intercellular adhesion molecule 5 (ICAM5) in THCA. CpG islands with aberrant methylation pattern in THCA, and the expression profiles of the corresponding genes in THCA, were analyzed using bioinformatics. ICAM5 was suggested to have a hypermethylation status, and it was highly expressed in THCA tissues and cells. Its overexpression promoted proliferation, mobility, and tumorigenic activity of THCA cells. As for the downstream signaling, ICAM5 was found to activate the MAPK/ERK and MAPK/JNK signaling pathways. Either inhibition of ERK or JNK blocked the oncogenic effects of ICAM5. DNA methyltransferases 1 (DNMT1) and DNMT3a were found to induce promoter hypermethylation of ICAM5 in THCA cells. Knockdown of DNMT1 or DNMT3a decreased the ICAM5 expression and suppressed malignant properties of THCA cells in vitro and in vivo, which were, however, restored by further artificial ICAM5 overexpression. Collectively, this study reveals that DNMT1 and DNMT3a mediates promoter hypermethylation and transcription activation of ICAM5 in THCA, which promotes malignant progression of THCA through the MAPK signaling pathway.

Analysis of the potential biological value of pyruvate dehydrogenase E1 subunit β in human cancer

BACKGROUND

The pyruvate dehydrogenase E1 subunit β (PDHB) gene which regulates energy metabolism is located in mitochondria. However, few studies have elucidated the role and mechanism of PDHB in different cancers.

AIM

To comprehensive pan-cancer analysis of PDHB was performed based on bioinformatics approaches to explore its tumor diagnostic and prognostic value and tumor immune relevance in cancer. In vitro experiments were performed to examine the biological regulation of PDHB in liver cancer.

METHODS

Pan-cancer data related to PDHB were obtained from the Cancer Genome Atlas (TCGA) database. Analysis of the gene expression profiles of PDHB was based on TCGA and Genotype Tissue Expression Dataset databases. Cox regression analysis and Kaplan-Meier methods were used to assess the correlation between PDHB expression and survival prognosis in cancer patients. The correlation between PDHB and receiver operating characteristic diagnostic curve, clinicopathological staging, somatic mutation, tumor mutation burden (TMB), microsatellite instability (MSI), DNA methylation, and drug susceptibility in pan-cancer was also analyzed. Various algorithms were used to analyze the correlation between PDHB and immune cell infiltration and tumor chemotaxis environment, as well as the co-expression analysis of PDHB and immune checkpoint (ICP) genes. The expression and functional phenotype of PDHB in single tumor cells were studied by single-cell sequencing, and the functional enrichment analysis of PDHB-related genes was performed. The study also validated the level of mRNA or protein expression of PDHB in several cancers. Finally, in vitro experiments verified the regulatory effect of PDHB on the proliferation, migration, and invasion of liver cancer.

RESULTS

PDHB was significantly and differently expressed in most cancers. PDHB was significantly associated with prognosis in patients with a wide range of cancers, including kidney renal clear cell carcinoma, kidney renal papillary cell carcinoma, breast invasive carcinoma, and brain lower grade glioma. In some cancers, PDHB expression was clearly associated with gene mutations, clinicopathological stages, and expression of TMB, MSI, and ICP genes. The expression of PDHB was closely related to the infiltration of multiple immune cells in the immune microenvironment and the regulation of tumor chemotaxis environment. In addition, single-cell sequencing results showed that PDHB correlated with different biological phenotypes of multiple cancer single cells. This study further demonstrated that down-regulation of PDHB expression inhibited the proliferation, migration, and invasion functions of hepatoma cells.

CONCLUSION

As a member of pan-cancer, PDHB may be a novel cancer marker with potential value in diagnosing cancer, predicting prognosis, and in targeted therapy.

Colorectal cancer with low SLC35A3 is associated with immune infiltrates and poor prognosis

The expression level of SLC35A3 is associated with the prognosis of many cancers, but its role in colorectal cancer (CRC) is unclear. The purpose of our study was to elucidate the role of SLC35A3 in CRC. The expression levels of SLC35A3 in CRC were evaluated through tumor immune resource assessment (TIMER), The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), International Cancer Genome Consortium (ICGC), Human Protein Atlas (HPA), qRT-PCR, and immunohistochemical evaluation. TCGA, GEO, and ICGC databases were used to analyze the diagnostic and prognostic value of SLC35A3 in CRC. A overall survival (OS) model was constructed and validated based on the expression level of SLC35A3 and multivariable analysis results. The cBioPortal tool was used to analyze SLC35A3 mutation in CRC. The UALCAN tool was used to analyze the promoter methylation level of SLC35A3 in colorectal cancer. In addition, the role of SLC35A3 in CRC was determined through GO analysis, KEGG analysis, gene set enrichment analysis (GSEA), immune infiltration analysis, and immune checkpoint correlation analysis. In vitro experiments validated the function of SLC35A3 in colorectal cancer cells. Compared with adjacent normal tissues and colonic epithelial cells, the expression of SLC35A3 was decreased in CRC tissues and CRC cell lines. Low expression of SLC35A3 was associated with N stage, pathological stage, and lymphatic infiltration, and it was unfavorable for OS, disease-specific survival (DSS), recurrence-free survival (RFS), and post-progression survival (PPS). According to the Receiver Operating Characteristic (ROC) analysis, SLC35A3 is a potential important diagnostic biomarker for CRC patients. The nomograph based on the expression level of SLC35A3 showed a better predictive model for OS than single prognostic factors and TNM staging. SLC35A3 has multiple types of mutations in CRC, and its promoter methylation level is significantly decreased. GO and KEGG analysis indicated that SLC35A3 may be involved in transmembrane transport protein activity, cell communication, and interaction with neurotransmitter receptors. GSEA revealed that SLC35A3 may be involved in energy metabolism, DNA repair, and cancer pathways. In addition, SLC35A3 was closely related to immune cell infiltration and immune checkpoint expression. Immunohistochemistry confirmed the positive correlation between SLC35A3 and helper T cell infiltration. In vitro experiments showed that overexpression of SLC35A3 inhibited the proliferation and invasion capability of colorectal cancer cells and promoted apoptosis. The results of this study indicate that decreased expression of SLC35A3 is closely associated with poor prognosis and immune cell infiltration in colorectal cancer, and it can serve as a promising independent prognostic biomarker and potential therapeutic target.

Non-coding RNAs in hypertrophic scars and keloids: Current research and clinical relevance: A review

Hypertrophic scars (HS) and keloids (KD) are lesions that develop as a result of excessive fibroblast proliferation and collagen deposition in response to dermal injury, leading to dysregulation of the inflammatory, proliferative, and remodeling phases during wound healing. HS and KD affect up to 90 % of the population and are associated with lower quality of life, physical health, and mental status in patients. Efficient targeted treatment represents a significant challenge, primarily due to our limited understanding of their underlying pathogenesis. Non-coding RNAs (ncRNAs), which constitute a significant portion of the human transcriptome with minimal or no protein-coding capacity, have been implicated in various cellular physiologies and pathologies and may serve as diagnostic indicators or therapeutic targets. NcRNAs have been found to be aberrantly expressed and regulated in HS and KD. This review provides a summary of the expression profiles and molecular mechanisms of three common ncRNAs, including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs), in HS and KD. It also discusses their potential as biomarkers for the diagnosis and treatment of these diseases and provides novel insights into epigenetic-based diagnosis and treatment strategies for HS and KD.

Detection and clinical significance of CEACAM5 methylation in colorectal cancer patients

Colorectal cancer (CRC) is a globally common cancer, and the serum carcinoembryonic antigen (sCEA) is widely applied as a diagnostic and prognostic tumor marker in CRC. This study aimed to elucidate the mechanism of CEA expression and corresponding clinical features to improve prognostic assessments. In CRC cells, hypomethylation of the CEACAM5 promoter enhanced CEA expression in HCT116 and HT29 cells with 5-aza-2'-deoxycytidine (5-Aza-dC) treatment. Our clinical data indicated that 64.7% (101/156) of CRC patients had an sCEA level above the normal range, and 76.2% (77/101) of those patients showed a lower average CpG methylation level of the CEACAM5 promoter. The methylation analysis showed that both CRC cell lines and patient samples shared the same critical methylation CpG regions at -200 to -500 and -1000 to -1400 bp of the CEACAM5 promoter. Patients with hypermethylation of the CEACAM5 promoter showed features of a BRAF mutation, TGFB2 mutation, microsatellite instability-high, and preference for right-sided colorectal cancer and peritoneal seeding presentation that had a similar clinical character to the consensus molecular subtype 1 (CMS1) of colorectal cancer. Additionally, hypermethylation of the CEACAM5 promoter combined with evaluated sCEA demonstrated the worst survival among the patients. Therefore, the methylation status of the CEACAM5 promoter also served as an effective biomarker for assessing disease prognosis. Results indicated that DNA methylation is a major regulatory mechanism for CEA expression in colorectal cancer. Moreover, our data also highlighted that patients in a subgroup who escaped from inactivation by DNA methylation had distinct clinical and pathological features and the worst survival.

Comprehensive analysis of the potential biological significance of cuproptosis-related gene LIPT2 in pan-cancer prognosis and immunotherapy

Lipoyltransferase 2 (LIPT2) acts as a key enzyme involved in fatty acid metabolism and cell membrane synthesis. However, the biological function of LIPT2 in various cancer types and its potential significance in prognosis continue to be unresolved. For this analysis, we evaluated the expression levels and the significance of prognosis of LIPT2 gene in all cancers by various bioinformatics methods. The results found that LIPT2 was dramatically overexpressed in the vast majority of cancers. The upregulated LIPT2 was related to bad prognosis in Brain Lower Grade Glioma (LGG), Glioma (GBMLGG), Glioblastoma multiforme (GBM), Kidney Chromophobe (KICH), and High-Risk Wilms Tumor (WT), while it had a favorable prognosis in Kidney renal clear cell carcinoma (KIRC), and Ovarian serous cystadenocarcinoma (OV), Pan-kidney cohort (KIPAN). Furthermore, we assessed the mutation status, methylation levels, and immune status of LIPT2 in pan-cancer. Single-cell sequencing results revealed the correlation of LIPT2 expression with various biological characteristics such as DNA lesion, tumor angiogenesis, cell apoptosis, metastasis, and invasion. Enrichment analysis unveiled potential molecular regulatory mechanisms. In conclusion, our research reveals a detailed key role of LIPT2 in the progression, prognosis, and immune efficacy of various forms of cancer. Therefore, we have reason to believe that LIPT2 has the potential to be a candidate biomarker for tumors.

High-resolution DNA methylation screening of the major histocompatibility complex in multiple sclerosis

Background

The HLA-DRB1 gene in the major histocompatibility complex (MHC) region in chromosome 6p21 is the strongest genetic factor identified as influencing multiple sclerosis (MS) susceptibility. DNA methylation changes associated with MS have been consistently detected at the MHC region. However, understanding the full scope of epigenetic regulations of the MHC remains incomplete, due in part to the limited coverage of this region by standard whole genome bisulfite sequencing or array-based methods.

Methods

We developed and validated an MHC capture protocol coupled with bisulfite sequencing and conducted a comprehensive analysis of the MHC methylation landscape in blood samples from 147 treatment naïve MS study participants and 129 healthy controls.

Results

We identified 132 differentially methylated region (DMRs) within MHC region associated with disease status. The DMRs overlapped with established MS risk loci. Integration of the MHC methylome with human leukocyte antigen (HLA) genetic data indicate that the methylation changes are significantly associated with HLA genotypes. Using DNA methylation quantitative trait loci (mQTL) mapping and the causal inference test (CIT), we identified 643 cis-mQTL-DMRs paired associations, including 71 DMRs possibly mediating causal relationships between 55 single nucleotide polymorphisms (SNPs) and MS risk.

Results

The results describe MS-associated methylation changes in MHC region and highlight the association between HLA genotypes and methylation changes. Results from the mQTL and CIT analyses provide evidence linking MHC region variations, methylation changes, and disease risk for MS.

pH-sensing G protein-coupled orphan receptor GPR68 is expressed in human cartilage and correlates with degradation of extracellular matrix during OA progression

Background

Osteoarthritis (OA) is a debilitating joints disease affecting millions of people worldwide. As OA progresses, chondrocytes experience heightened catabolic activity, often accompanied by alterations in the extracellular environment's osmolarity and acidity. Nevertheless, the precise mechanism by which chondrocytes perceive and respond to acidic stress remains unknown. Recently, there has been growing interest in pH-sensing G protein-coupled receptors (GPCRs), such as GPR68, within musculoskeletal tissues. However, function of GPR68 in cartilage during OA progression remains unknown. This study aims to identify the role of GPR68 in regulation of catabolic gene expression utilizing an in vitro model that simulates catabolic processes in OA.

Methods

We examined the expression of GPCR by analyzing high throughput RNA-Seq data in human cartilage isolated from healthy donors and OA patients. De-identified and discarded OA cartilage was obtained from joint arthroplasty and chondrocytes were prepared by enzymatic digestion. Chondrocytes were treated with GPR68 agonist, Ogerin and then stimulated IL1β and RNA isolation was performed using Trizol method. Reverse transcription was done using the cDNA synthesis kit and the expression of GPR68 and OA related catabolic genes was quantified using SYBR® green assays.

Results

The transcriptome analysis revealed that pH sensing GPCR were expressed in human cartilage with a notable increase in the expression of GPR68 in OA cartilage which suggest a potential role for GPR68 in the pathogenesis of OA. Immunohistochemical (IHC) and qPCR analyses in human cartilage representing various stages of OA indicated a progressive increase in GPR68 expression in cartilage associated with higher OA grades, underscoring a correlation between GPR68 expression and the severity of OA. Furthermore, IHC analysis of Gpr68 in murine cartilage subjected to surgically induced OA demonstrated elevated levels of GPR68 in knee cartilage and meniscus. Using IL1β stimulated in vitro model of OA catabolism, our qPCR analysis unveiled a time-dependent increase in GPR68 expression in response to IL1β stimulation, which correlates with the expression of matrix degrading proteases suggesting the role of GPR68 in chondrocytes catabolism and matrix degeneration. Using pharmacological activator of GPR68, our results further showed that GPR68 activation repressed the expression of MMPs in human chondrocytes.

Conclusions

Our results demonstrated that GPR68 was robustly expressed in human cartilage and mice and its expression correlates with matrix degeneration and severity of OA progression in human and surgical model. GPR68 activation in human chondrocytes further repressed the expression of MMPs under OA pathological condition. These results identify GPR68 as a possible therapeutic target in the regulation of matrix degradation during OA.

DNA Methylation Profiling of Ovarian Tissue of Climbing Perch (Anabas testudienus) in Response to Monocrotophos Exposure

Epigenetic modifications like DNA methylation can alter an organism's phenotype without changing its DNA sequence. Exposure to environmental toxicants has the potential to change the resilience of aquatic species. However, little information is available on the dynamics of DNA methylation in fish gonadal tissues in response to organophosphates. In the present work, reduced-representation bisulfite sequencing was performed to identify DNA methylation patterns in the ovarian tissues of Anabas testudienus exposed to organophosphates, specifically monocrotophos (MCP). Through sequencing, an average of 41,087 methylated cytosine sites were identified and distributed in different parts of genes, i.e., in transcription start sites (TSS), promoters, exons, etc. A total of 1058 and 1329 differentially methylated regions (DMRs) were detected as hyper-methylated and hypo-methylated in ovarian tissues, respectively. Utilizing whole-genome data of the climbing perch, the DMRs, and their associated overlapping genes revealed a total of 22 genes within exons, 45 genes at transcription start sites (TSS), and 218 genes in intergenic regions. Through gene ontology analysis, a total of 16 GO terms particularly involved in ovarian follicular development, response to oxidative stress, oocyte maturation, and multicellular organismal response to stress associated with reproductive biology were identified. After functional enrichment analysis, relevant DMGs such as steroid hormone biosynthesis (Cyp19a, 11-beta-HSD, 17-beta-HSD), hormone receptors (ar, esrrga), steroid metabolism (StAR), progesterone-mediated oocyte maturation (igf1ar, pgr), associated with ovarian development in climbing perch showed significant differential methylation patterns. The differentially methylated genes (DMGs) were subjected to analysis using real-time PCR, which demonstrated altered gene expression levels. This study revealed a molecular-level alteration in genes associated with ovarian development in response to chemical exposure. This work provides evidence for understanding the relationship between DNA methylation and gene regulation in response to chemicals that affect the reproductive fitness of aquatic animals.

Age-related promoter-switch regulates Runx1 expression in adult rat hearts

BACKGROUND

Runt-related transcription factor-1 (RUNX1), a key member of the core-binding factor family of transcription factors, has emerged as a novel therapeutic target for cardiovascular disease. There is an urgent need to fully understand the expression pattern of Runx1 in the heart and the mechanisms by which it is controlled under normal conditions and in response to disease. The expression of Runx1 is regulated at the transcriptional level by two promoters designated P1 and P2. Alternative usage of these two promoters creates differential mRNA transcripts diversified in distribution and translational potential. While the significance of P1/P2 promoter-switch in the transcriptional control of Runx1 has been highlighted in the embryogenic process, very little is known about the level of P1- and P2-specific transcripts in adult hearts, and the underlying mechanisms controlling the promoter-switch.

METHODS

To amplify P1/P2 specific sequences in the heart, we used two different sense primers complementary to either P1 or P2 5'-regions to monitor the expression of P1/P2 transcripts. DNA methylation levels were assessed at the Runx1 promoter regions. Rats were grouped by age.

RESULTS

The expression levels of both P1- and P2-derived Runx1 transcripts were decreased in older rats when compared with that in young adults, paralleled with an age-dependent decline in Runx1 protein level. Furthermore, older rats demonstrated a higher degree of DNA methylation at Runx1 promoter regions. Alternative promoter usage was observed in hearts with increased age, as reflected by altered P1:P2 mRNA ratio.

CONCLUSION

Our data demonstrate that the expression of Runx1 in the heart is age-dependent and underscore the importance of gene methylation in the promoter-mediated transcriptional control of Runx1, thereby providing new insights to the role of epigenetic regulation in the heart.

Integrated molecular analysis reveals hypermethylation and overexpression of HOX genes to be poor prognosticators in isocitrate dehydrogenase mutant glioma

BACKGROUND

Diffuse gliomas represent over 80% of malignant brain tumors ranging from low-grade to aggressive high-grade lesions. Within isocitrate dehydrogenase (IDH)-mutant gliomas, there is a high variability in survival and a need to more accurately predict outcome.

METHODS

To identify and characterize a predictive signature of outcome in gliomas, we utilized an integrative molecular analysis (using methylation, mRNA, copy number variation (CNV), and mutation data), analyzing a total of 729 IDH-mutant samples including a test set of 99 from University Health Network (UHN) and 2 validation cohorts including the German Cancer Research Center (DKFZ) and The Cancer Genome Atlas (TCGA).

RESULTS

Cox regression analysis of methylation data from the UHN cohort identified CpG-based signatures that split the glioma cohort into 2 prognostic groups strongly predicting survival that were validated using 2 independent cohorts from TCGA and DKFZ (all P-values < .0001). The methylation signatures that predicted poor outcomes also exhibited high CNV instability and hypermethylation of HOX gene probes. Integrated multi-platform analyses using mRNA and methylation (iRM) showed that parallel HOX gene overexpression and simultaneous hypermethylation were significantly associated with increased mutational load, high aneuploidy, and worse survival (P-value < .0001). A 7-HOX gene signature was developed and validated using the most significantly associated HOX genes with patient outcome in both 1p/19q codeleted and non-codeleted IDHmut gliomas.

CONCLUSIONS

HOX gene methylation and expression provide important prognostic information in IDH-mutant gliomas that are not captured by current molecular diagnostics. A 7-HOX gene signature of outcome shows significant survival differences in both 1p/19q codeleted and non-codeleted IDH-mutant gliomas.

The comprehensive analysis of the prognostic and functional role of N-terminal methyltransferases 1 in pan-cancer

Background

NTMT1, a transfer methylase that adds methyl groups to the N-terminus of proteins, has been identified as a critical player in tumor development and progression. However, its precise function in pan-cancer is still unclear. To gain a more comprehensive understanding of its role in cancer, we performed a thorough bioinformatics analysis.

Methods

To conduct our analysis, we gathered data from multiple sources, including RNA sequencing and clinical data from the TCGA database, protein expression data from the UALCAN and HPA databases, and single-cell expression data from the CancerSEA database. Additionally, we utilized TISIDB to investigate the interaction between the tumor and the immune system. To assess the impact of NTMT1 on the proliferation of SNU1076 cells, we performed a CCK8 assay. We also employed cellular immunofluorescence to detect DNA damage and used flow cytometry to measure tumor cell apoptosis.

Results

Our analysis revealed that NTMT1 was significantly overexpressed in various types of tumors and that high levels of NTMT1 were associated with poor survival outcomes. Functional enrichment analysis indicated that NTMT1 may contribute to tumor development and progression by regulating pathways involved in cell proliferation and immune response. In addition, we found that knockdown of NTMT1 expression led to reduced cell proliferation, increased DNA damage, and enhanced apoptosis in HNSCC cells.

Conclusion

High expression of NTMT1 in tumors is associated with poor prognosis. The underlying regulatory mechanism of NTMT1 in cancer is complex, and it may be involved in both the promotion of tumor development and the inhibition of the tumor immune microenvironment.

Insight into the Epigenetics of Kaposi's Sarcoma-Associated Herpesvirus

Epigenetic reprogramming represents a series of essential events during many cellular processes including oncogenesis. The genome of Kaposi's sarcoma-associated herpesvirus (KSHV), an oncogenic herpesvirus, is predetermined for a well-orchestrated epigenetic reprogramming once it enters into the host cell. The initial epigenetic reprogramming of the KSHV genome allows restricted expression of encoded genes and helps to hide from host immune recognition. Infection with KSHV is associated with Kaposi's sarcoma, multicentric Castleman's disease, KSHV inflammatory cytokine syndrome, and primary effusion lymphoma. The major epigenetic modifications associated with KSHV can be labeled under three broad categories: DNA methylation, histone modifications, and the role of noncoding RNAs. These epigenetic modifications significantly contribute toward the latent-lytic switch of the KSHV lifecycle. This review gives a brief account of the major epigenetic modifications affiliated with the KSHV genome in infected cells and their impact on pathogenesis.

Comprehensive analysis of TLX2 in pan cancer as a prognostic and immunologic biomarker and validation in ovarian cancer

T cell leukemia homeobox 2 (TLX2) plays an important role in some tumors. Bioinformatics and experimental validation represent a useful way to explore the mechanisms and functions of TLX2 gene in the cancer disease process from a pan cancer perspective. TLX2 was aberrantly expressed in pan cancer and cell lines and correlated with clinical stage. High TLX2 expression was significantly associated with poor overall survival in COAD, KIRC, OC, and UCS. The greatest frequency of TLX2 alterations in pan cancer was amplification. Alterations of NXF2B, MSLNL, PCGF1, INO80B-WBP1, LBX2-AS1, MRPL53, LBX2, TTC31, WDR54, and WBP1 co-occurred in the TLX2 alteration group. PFS was significantly shorter in the TLX2-altered group (n = 6) compared to the TLX2-unaltered group (n = 400). Methylation levels of TLX2 were high in 17 tumors. TLX2 expression was associated with MSI in seven tumors and TMB in five tumors. TLX2 expression was associated with immune infiltration and immune checkpoint genes. TLX2 may be associated with some pathways and chemoresistance. We constructed a possible competing endogenous RNA (ceRNA) network of LINC01010/miR-146a-5p/TLX2 in OC. TLX2 expression was significantly upregulated in ovarian cancer cell lines compared to ovarian epithelial cell lines. Aberrant expression of TLX2 in pan cancer may promote tumorigenesis and progression through different mechanisms. TLX2 may represent an important therapeutic target for human cancers.

Association between DNA Methylation in the Core Promoter Region of the CUT-like Homeobox 1 (CUX1) Gene and Lambskin Pattern in Hu Sheep

CUT-like homeobox 1 (CUX1) has been proven to be a key regulator in sheep hair follicle development. In our previous study, CUX1 was identified as a differential expressed gene between Hu sheep lambskin with small wave patterns (SM) and straight wool patterns (ST); however, the exact molecular mechanism of CUX1 expression has been obscure. As DNA methylation can regulate the gene expression, the potential association between CUX1 core promotor region methylation and lambskin pattern in Hu sheep was explored in the present study. The results show that the core promoter region of CUX1 was present at (-1601-(-1) bp) upstream of the transcription start site. A repressive region (-1151-(-751) bp) was also detected, which had a strong inhibitory effect on CUX1 promoter activity. Bisulfite amplicon sequencing revealed that no significant difference was detected between the methylation levels of CUX1 core promoter region in SM tissues and ST tissues. Although the data demonstrated the differential expression of CUX1 between SM and ST probably has no association with DNA methylation, the identification of the core region and a potential repressive region of CUX1 promoter can enrich the role of CUX1 in Hu sheep hair follicle development.

The role of salinity on genome-wide DNA methylation dynamics in European sea bass gills

Epigenetic modifications, like DNA methylation, generate phenotypic diversity in fish and ultimately lead to adaptive evolutionary processes. Euryhaline marine species that migrate between salinity-contrasted habitats have received little attention regarding the role of salinity on whole-genome DNA methylation. Investigation of salinity-induced DNA methylation in fish will help to better understand the potential role of this process in salinity acclimation. Using whole-genome bisulfite sequencing, we compared DNA methylation patterns in European sea bass (Dicentrarchus labrax) juveniles in seawater and after freshwater transfer. We targeted the gill as a crucial organ involved in plastic responses to environmental changes. To investigate the function of DNA methylation in gills, we performed RNAseq and assessed DNA methylome-transcriptome correlations. We showed a negative correlation between gene expression levels and DNA methylation levels in promoters, first introns and first exons. A significant effect of salinity on DNA methylation dynamics with an overall DNA hypomethylation in freshwater-transferred fish compared to seawater controls was demonstrated. This suggests a role of DNA methylation changes in salinity acclimation. Genes involved in key functions as metabolism, ion transport and transepithelial permeability (junctional complexes) were differentially methylated and expressed between salinity conditions. Expression of genes involved in mitochondrial metabolism (tricarboxylic acid cycle) was increased, whereas the expression of DNA methyltransferases 3a was repressed. This study reveals novel links between DNA methylation, mainly in promoters and first exons/introns, and gene expression patterns following salinity change.

High temperature influences DNA methylation and transcriptional profiles in sea urchins (Strongylocentrotus intermedius)

BACKGROUND

DNA methylation plays an important role in life processes by affecting gene expression, but it is still unclear how DNA methylation is controlled and how it regulates gene transcription under high temperature stress conditions in Strongylocentrotus intermedius. The potential link between DNA methylation variation and gene expression changes in response to heat stress in S. intermedius was investigated by MethylRAD-seq and RNA-seq analysis. We screened DNA methylation driver genes in order to comprehensively elucidate the regulatory mechanism of its high temperature adaptation at the DNA/RNA level.

RESULTS

The results revealed that high temperature stress significantly affected not only the DNA methylation and transcriptome levels of S. intermedius (P < 0.05), but also growth. MethylRAD-seq analysis revealed 12,129 CG differential methylation sites and 966 CWG differential methylation sites, and identified a total of 189 differentially CG methylated genes and 148 differentially CWG methylated genes. Based on KEGG enrichment analysis, differentially expressed genes (DEGs) are mostly enriched in energy and cell division, immune, and neurological damage pathways. Further RNA-seq analysis identified a total of 1968 DEGs, of which 813 genes were upregulated and 1155 genes were downregulated. Based on the joint MethylRAD-seq and RNA-seq analysis, metabolic processes such as glycosaminoglycan degradation, oxidative phosphorylation, apoptosis, glutathione metabolism, thermogenesis, and lysosomes are regulated by DNA methylation.

CONCLUSIONS

High temperature affected the DNA methylation and expression levels of genes such as MOAP-1, GGT1 and RDH8, which in turn affects the metabolism of HPSE, Cox, glutathione, and retinol, thereby suppressing the immune, energy metabolism, and antioxidant functions of the organism and finally manifesting as stunted growth. In summary, the observations in the present study improve our understanding of the molecular mechanism of the response to high temperature stress in sea urchin.

The pharmacoepigenomic landscape of cancer cell lines reveals the epigenetic component of drug sensitivity

Aberrant DNA methylation accompanies genetic alterations during oncogenesis and tumour homeostasis and contributes to the transcriptional deregulation of key signalling pathways in cancer. Despite increasing efforts in DNA methylation profiling of cancer patients, there is still a lack of epigenetic biomarkers to predict treatment efficacy. To address this, we analyse 721 cancer cell lines across 22 cancer types treated with 453 anti-cancer compounds. We systematically detect the predictive component of DNA methylation in the context of transcriptional and mutational patterns, i.e., in total 19 DNA methylation biomarkers across 17 drugs and five cancer types. DNA methylation constitutes drug sensitivity biomarkers by mediating the expression of proximal genes, thereby enhancing biological signals across multi-omics data modalities. Our method reproduces anticipated associations, and in addition, we find that the NEK9 promoter hypermethylation may confer sensitivity to the NEDD8-activating enzyme (NAE) inhibitor pevonedistat in melanoma through downregulation of NEK9. In summary, we envision that epigenomics will refine existing patient stratification, thus empowering the next generation of precision oncology.

Reverse-ChIP Techniques for Identifying Locus-Specific Proteomes: A Key Tool in Unlocking the Cancer Regulome

A phenotypic hallmark of cancer is aberrant transcriptional regulation. Transcriptional regulation is controlled by a complicated array of molecular factors, including the presence of transcription factors, the deposition of histone post-translational modifications, and long-range DNA interactions. Determining the molecular identity and function of these various factors is necessary to understand specific aspects of cancer biology and reveal potential therapeutic targets. Regulation of the genome by specific factors is typically studied using chromatin immunoprecipitation followed by sequencing (ChIP-Seq) that identifies genome-wide binding interactions through the use of factor-specific antibodies. A long-standing goal in many laboratories has been the development of a 'reverse-ChIP' approach to identify unknown binding partners at loci of interest. A variety of strategies have been employed to enable the selective biochemical purification of sequence-defined chromatin regions, including single-copy loci, and the subsequent analytical detection of associated proteins. This review covers mass spectrometry techniques that enable quantitative proteomics before providing a survey of approaches toward the development of strategies for the purification of sequence-specific chromatin as a 'reverse-ChIP' technique. A fully realized reverse-ChIP technique holds great potential for identifying cancer-specific targets and the development of personalized therapeutic regimens.

Genetic and epigenetic features of neuroendocrine prostate cancer and their emerging applications

Prostate cancer is the second most prevalent cancer in men globally. De novo neuroendocrine prostate cancer (NEPC) is uncommon at initial diagnosis, however, (treatment-induced) t-NEPC emerges in up to 25% of prostate adenocarcinoma (PRAD) cases treated with androgen deprivation, carrying a drastically poor prognosis. The transition from PRAD to t-NEPC is underpinned by several key genetic mutations; TP53, RB1, and MYCN are the main genes implicated, bearing similarities to other neuroendocrine tumours. A broad range of epigenetic alterations, such as aberrations in DNA methylation, histone post-translational modifications, and non-coding RNAs, may drive lineage plasticity from PRAD to t-NEPC. The clinical diagnosis of NEPC is hampered by a lack of accessible biomarkers; recent advances in liquid biopsy techniques assessing circulating tumour cells and ctDNA in NEPC suggest that the advent of non-invasive means of monitoring progression to NEPC is on the horizon. Such techniques are vital for NEPC management; diagnosis of t-NEPC is crucial for implementing effective treatment, and precision medicine will be integral to providing the best outcomes for patients.

A Novel Repetition Frequency-Based DNA Encoding Scheme to Predict Human and Mouse DNA Enhancers with Deep Learning

Recent studies have shown that DNA enhancers have an important role in the regulation of gene expression. They are responsible for different important biological elements and processes such as development, homeostasis, and embryogenesis. However, experimental prediction of these DNA enhancers is time-consuming and costly as it requires laboratory work. Therefore, researchers started to look for alternative ways and started to apply computation-based deep learning algorithms to this field. Yet, the inconsistency and unsuccessful prediction performance of computational-based approaches among various cell lines led to the investigation of these approaches as well. Therefore, in this study, a novel DNA encoding scheme was proposed, and solutions were sought to the problems mentioned and DNA enhancers were predicted with BiLSTM. The study consisted of four different stages for two scenarios. In the first stage, DNA enhancer data were obtained. In the second stage, DNA sequences were converted to numerical representations by both the proposed encoding scheme and various DNA encoding schemes including EIIP, integer number, and atomic number. In the third stage, the BiLSTM model was designed, and the data were classified. In the final stage, the performance of DNA encoding schemes was determined by accuracy, precision, recall, F1-score, CSI, MCC, G-mean, Kappa coefficient, and AUC scores. In the first scenario, it was determined whether the DNA enhancers belonged to humans or mice. As a result of the prediction process, the highest performance was achieved with the proposed DNA encoding scheme, and an accuracy of 92.16% and an AUC score of 0.85 were calculated, respectively. The closest accuracy score to the proposed scheme was obtained with the EIIP DNA encoding scheme and the result was observed as 89.14%. The AUC score of this scheme was measured as 0.87. Among the remaining DNA encoding schemes, the atomic number showed an accuracy score of 86.61%, while this rate decreased to 76.96% with the integer scheme. The AUC values of these schemes were 0.84 and 0.82, respectively. In the second scenario, it was determined whether there was a DNA enhancer and, if so, it was decided to which species this enhancer belonged. In this scenario, the highest accuracy score was obtained with the proposed DNA encoding scheme and the result was 84.59%. Moreover, the AUC score of the proposed scheme was determined as 0.92. EIIP and integer DNA encoding schemes showed accuracy scores of 77.80% and 73.68%, respectively, while their AUC scores were close to 0.90. The most ineffective prediction was performed with the atomic number and the accuracy score of this scheme was calculated as 68.27%. Finally, the AUC score of this scheme was 0.81. At the end of the study, it was observed that the proposed DNA encoding scheme was successful and effective in predicting DNA enhancers.

Profiling the immune epigenome across global cattle breeds

BACKGROUND

Understanding the variation between well and poorly adapted cattle breeds to local environments and pathogens is essential for breeding cattle with improved climate and disease-resistant phenotypes. Although considerable progress has been made towards identifying genetic differences between breeds, variation at the epigenetic and chromatin levels remains poorly characterized. Here, we generate, sequence and analyse over 150 libraries at base-pair resolution to explore the dynamics of DNA methylation and chromatin accessibility of the bovine immune system across three distinct cattle lineages.

RESULTS

We find extensive epigenetic divergence between the taurine and indicine cattle breeds across immune cell types, which is linked to the levels of local DNA sequence divergence between the two cattle sub-species. The unique cell type profiles enable the deconvolution of complex cellular mixtures using digital cytometry approaches. Finally, we show distinct sub-categories of CpG islands based on their chromatin and methylation profiles that discriminate between classes of distal and gene proximal islands linked to discrete transcriptional states.

CONCLUSIONS

Our study provides a comprehensive resource of DNA methylation, chromatin accessibility and RNA expression profiles of three diverse cattle populations. The findings have important implications, from understanding how genetic editing across breeds, and consequently regulatory backgrounds, may have distinct impacts to designing effective cattle epigenome-wide association studies in non-European breeds.

Partial erosion on under-methylated regions and chromatin reprogramming contribute to oncogene activation in IDH mutant gliomas

BACKGROUND

IDH1/2 hotspot mutations are well known to drive oncogenic mutations in gliomas and are well-defined in the WHO 2021 classification of central nervous system tumors. Specifically, IDH mutations lead to aberrant hypermethylation of under-methylated regions (UMRs) in normal tissues through the disruption of TET enzymes. However, the chromatin reprogramming and transcriptional changes induced by IDH-related hypermethylation in gliomas remain unclear.

RESULTS

Here, we have developed a precise computational framework based on Hidden Markov Model to identify altered methylation states of UMRs at single-base resolution. By applying this framework to whole-genome bisulfite sequencing data from 75 normal brain tissues and 15 IDH mutant glioma tissues, we identified two distinct types of hypermethylated UMRs in IDH mutant gliomas. We named them partially hypermethylated UMRs (phUMRs) and fully hypermethylated UMRs (fhUMRs), respectively. We found that the phUMRs and fhUMRs exhibit distinct genomic features and chromatin states. Genes related to fhUMRs were more likely to be repressed in IDH mutant gliomas. In contrast, genes related to phUMRs were prone to be up-regulated in IDH mutant gliomas. Such activation of phUMR genes is associated with the accumulation of active H3K4me3 and the loss of H3K27me3, as well as H3K36me3 accumulation in gene bodies to maintain gene expression stability. In summary, partial erosion on UMRs was accompanied by locus-specific changes in key chromatin marks, which may contribute to oncogene activation.

CONCLUSIONS

Our study provides a computational strategy for precise decoding of methylation encroachment patterns in IDH mutant gliomas, revealing potential mechanistic insights into chromatin reprogramming that contribute to oncogenesis.

Emerging role of non-invasive and liquid biopsy biomarkers in pancreatic cancer

A global increase in the incidence of pancreatic cancer (PanCa) presents a major concern and health burden. The traditional tissue-based diagnostic techniques provided a major way forward for molecular diagnostics; however, they face limitations based on diagnosis-associated difficulties and concerns surrounding tissue availability in the clinical setting. Late disease development with asymptomatic behavior is a drawback in the case of existing diagnostic procedures. The capability of cell free markers in discriminating PanCa from autoimmune pancreatitis and chronic pancreatitis along with other precancerous lesions can be a boon to clinicians. Early-stage diagnosis of PanCa can be achieved only if these biomarkers specifically discriminate the non-carcinogenic disease stage from malignancy with respect to tumor stages. In this review, we comprehensively described the non-invasive disease detection approaches and why these approaches are gaining popularity for their early-stage diagnostic capability and associated clinical feasibility.

Epigenetics Approaches toward Precision Medicine for Idiopathic Pulmonary Fibrosis: Focus on DNA Methylation

Genetic information is not transmitted solely by DNA but by the epigenetics process. Epigenetics describes molecular missing link pathways that could bridge the gap between the genetic background and environmental risk factors that contribute to the pathogenesis of pulmonary fibrosis. Specific epigenetic patterns, especially DNA methylation, histone modifications, long non-coding, and microRNA (miRNAs), affect the endophenotypes underlying the development of idiopathic pulmonary fibrosis (IPF). Among all the epigenetic marks, DNA methylation modifications have been the most widely studied in IPF. This review summarizes the current knowledge concerning DNA methylation changes in pulmonary fibrosis and demonstrates a promising novel epigenetics-based precision medicine.

HIG1 domain family member 1A is a crucial regulator of disorders associated with hypoxia

Mitochondria play a central role in cellular energy conversion, metabolism, and cell proliferation. The regulation of mitochondrial function by HIGD1A, which is located on the inner membrane of the mitochondria, is essential to maintain cell survival under hypoxic conditions. In recent years, there have been shown other cellular pathways and mechanisms involving HIGD1A diametrically or through its interaction. As a novel regulator, HIGD1A maintains mitochondrial integrity and enhances cell viability under hypoxic conditions, increasing cell resistance to hypoxia. HIGD1A mainly targets cytochrome c oxidase by regulating downstream signaling pathways, which affects the ATP generation system and subsequently alters mitochondrial respiratory function. In addition, HIGD1A plays a dual role in cell survival in distinct degree hypoxia regions of the tumor. Under mild and moderate anoxic areas, HIGD1A acts as a positive regulator to promote cell growth. However, HIGD1A plays a role in inhibiting cell growth but retaining cellular activity under severe anoxic areas. We speculate that HIGD1A engages in tumor recurrence and drug resistance mechanisms. This review will focus on data concerning how HIGD1A regulates cell viability under hypoxic conditions. Therefore, HIGD1A could be a potential therapeutic target for hypoxia-related diseases.

Aberrant DNA Methylation Patterns of Deleted in Liver Cancer 1 Isoforms in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC), a common primary liver cancer, is the third leading cause of death worldwide. DNA methylation changes are common in HCC and have been studied to be associated with hepatocarcinogenesis. In our study, we used the MassARRAY^® EpiTYPER technology to investigate the methylation differences of deleted in liver cancer 1 (DLC1) (isoform 1 and 3) promoter between HCC tissues and corresponding adjacent noncancerous tissues and the association between methylation levels and clinicopathological features. In addition, the modified CRISPR-Cas9 system and the DNA methyltransferase inhibitor (DNMTi) were utilized to explore the functional correlation of epigenetic modifications and DLC1 gene regulation. The methylation levels of the DLC1 isoforms in HCC samples were found significantly lower than those in the adjacent noncancerous tissues (all p < 0.0001). Also, we found that the expression of DLC1 could be bidirectionally regulated by the modified CRISPR-Cas9 system and the DNMTi. Moreover, the hypomethylation of DLC1 in HCC samples was connected with the presence of satellite lesions (p = 0.0305) and incomplete tumor capsule (p = 0.0204). Receiver operator characteristic curve analysis demonstrated that the methylation levels of DLC1 could be applied to discriminate HCC patients (area under the curve = 0.728, p < 0.0001). The hypomethylation status was a key regulatory mechanism of DLC1 expression and might serve as a potential biomarker for HCC.

Development of epigenetic biomarkers for the identification of sex and thermal stress in fish using DNA methylation analysis and machine learning procedures

The sex ratio is a key ecological demographic parameter crucial for population viability. However, the epigenetic mechanisms operating during gonadal development regulating gene expression and the sex ratio remain poorly understood. Moreover, there is interest in the development of epigenetic markers associated with a particular phenotype or as sentinels of environmental effects. Here, we profiled DNA methylation and gene expression of 10 key genes related to sex development and stress, including steroidogenic enzymes, and growth and transcription factors. We provide novel information on the sex-related differences and on the influence of elevated temperature on these genes in zebrafish, a species with mixed genetic and environmental influences on sex ratios. We identified both positive (e.g., amh, cyp11c and hsd11b2) and negative (e.g., cyp11a1 and dmrt1) correlations in unexposed males, and negative correlation (amh) in exposed females between DNA methylation and gene expression levels. Further, we combined DNA methylation analysis with machine learning procedures and found a series of informative CpGs capable not only of correctly identifying sex (based on cyp19a1a DNA methylation levels) but also of identifying whether males and females had been exposed to abnormally elevated temperature when young (based on amh and foxl2a DNA methylation levels, respectively). This was achieved in the absence of conspicuous morphological alterations of the gonads. These DNA methylation-based epigenetic biomarkers represent molecular resources that can correctly recapitulate past thermal history and pave the way for similar findings in other species to assess potential ecological effects of environmental disturbances in the context of climate change.

Risk Stratification of Pancreatic Neuroendocrine Neoplasms Based on Clinical, Pathological, and Molecular Characteristics

Pancreatic neuroendocrine neoplasms consist of heterogeneous diseases. Depending on the novel features detected by various modern technologies, their classification and related prognosis predictions continue to change and develop. The role of traditional clinicopathological prognostic factors, including classification systems, is also being refined, and several attempts have been made to predict a more accurate prognosis through novel serum biomarkers, genetic factors, and epigenetic factors that have been identified through various state-of-the-art molecular techniques with multiomics sequencing. In this review article, the latest research results including the traditional approach to prognostic factors and recent advanced strategies for risk stratification of pancreatic neuroendocrine neoplasms based on clinical, pathological, and molecular characteristics are summarized. Predicting prognosis through multi-factorial assessments seems to be more efficacious, and prognostic factors through noninvasive methods are expected to develop further advances in liquid biopsy in the future.

The DNA methylation landscape of the root-knot nematode-induced pseudo-organ, the gall, in Arabidopsis, is dynamic, contrasting over time, and critically important for successful parasitism

Root-knot nematodes (RKNs) induce giant cells (GCs) within galls which are characterized by large-scale gene repression at early stages. However, the epigenetic mechanism(s) underlying gene silencing is (are) still poorly characterized. DNA methylation in Arabidopsis galls induced by Meloidogyne javanica was studied at crucial infection stages (3 d post-infection (dpi) and 14 dpi) using enzymatic, cytological, and sequencing approaches. DNA methyltransferase mutants (met1, cmt2, cmt3, cmt2/3, drm1/2, ddc) and a DNA demethylase mutant (ros1), were analyzed for RKN resistance/tolerance, and galls were characterized by confocal microscopy and RNA-seq. Early galls were hypermethylated, and the GCs were found to be the major contributors to this hypermethylation, consistent with the very high degree of gene repression they exhibit. By contrast, medium/late galls showed no global increase in DNA methylation compared to uninfected roots, but exhibited large-scale redistribution of differentially methylated regions (DMRs). In line with these findings, it was also shown that DNA methylation and demethylation mutants showed impaired nematode reproduction and gall/GC-development. Moreover, siRNAs that were exclusively present in early galls accumulated at hypermethylated DMRs, overlapping mostly with retrotransposons in the CHG/CG contexts that might be involved in their repression, contributing to their stability/genome integrity. Promoter/gene methylation correlated with differentially expressed genes encoding proteins with basic cell functions. Both mechanisms are consistent with reprogramming host tissues for gall/GC formation. In conclusion, RNA-directed DNA methylation (RdDM; DRM2/1) pathways, maintenance methyltransferases (MET1/CMT3) and demethylation (ROS1) appear to be prominent mechanisms driving a dynamic regulation of the epigenetic landscape during RKN infection.

Genome-wide DNA methylation reveals potential epigenetic mechanism of age-dependent viral susceptibility in grass carp

Background

Grass carp are an important farmed fish in China that are infected by many pathogens, especially grass carp reovirus (GCRV). Notably, grass carp showed age-dependent susceptibility to GCRV; that is, grass carp not older than one year were sensitive to GCRV, while those over three years old were resistant to this virus. However, the underlying mechanism remains unclear. Herein, whole genome-wide DNA methylation and gene expression variations between susceptible five-month-old (FMO) and resistant three-year-old (TYO) grass carp were investigated aiming to uncover potential epigenetic mechanisms.

Results

Colorimetric quantification revealed that the global methylation level in TYO fish was higher than that in FMO fish. Whole-genome bisulfite sequencing (WGBS) of the two groups revealed 6214 differentially methylated regions (DMRs) and 4052 differentially methylated genes (DMGs), with most DMRs and DMGs showing hypermethylation patterns in TYO fish. Correlation analysis revealed that DNA hypomethylation in promoter regions and DNA hypermethylation in gene body regions were associated with gene expression. Enrichment analysis revealed that promoter hypo-DMGs in TYO fish were significantly enriched in typical immune response pathways, whereas gene body hyper-DMGs in TYO fish were significantly enriched in terms related to RNA transcription, biosynthesis, and energy production. RNA-seq analysis of the corresponding samples indicated that most of the genes in the above terms were upregulated in TYO fish. Moreover, gene function analysis revealed that the two genes involved in energy metabolism displayed antiviral effects.

Conclusions

Collectively, these results revealed genome-wide variations in DNA methylation between grass carp of different ages. DNA methylation and gene expression variations in genes involved in immune response, biosynthesis, and energy production may contribute to age-dependent susceptibility to GCRV in grass carp. Our results provide important information for disease-resistant breeding programs for grass carp and may also benefit research on age-dependent diseases in humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12979-022-00285-w.

MAPK8IP2 is a potential prognostic biomarker and promote tumor progression in prostate cancer

Background

MAPK8IP2 is one of the JNK-interacting proteins (JIPs) family members, and is involved in the regulation of the JNK and P38 MAPK signaling pathways. MAPK8IP2 has been reported to be closely associated with several cancers. However, the biological function of MAPK8IP2 in prostate cancer (PCa) remains unclear.

Methods

MAPK8IP2 expression in PCa and subgroups of PCa was analyzed by public databases. The prognostic role of MAPK8IP2 in prostate cancer was analyzed using the Cox regression method. The potential mechanism by which MAPK8IP2 affects PCa progression was investigated by utilizing public data, including genetic alteration, DNA methylation, m6A methylation, and immune infiltration data. We further performed in vitro assays to validate the effect of MAPK8IP2 on PCa cell proliferation, migration and invasion.

Results

MAPK8IP2 is highly expressed in PCa tissues. Overexpression of MAPK8IP2 is associated with adverse clinicopathological factors and a poor prognosis in PCa. Receiver operating curve analysis showed that MAPK8IP2 can distinguish PCa tissues from non-PCa tissues with a certain accuracy (AUC = 0.814). The MAPK8IP2 genetic alteration rate was 2.6% and MAPK8IP2 alterations correlated with a poor prognosis. We also found that CDK12 and TP53 mutations were associated with MAPK8IP2 expression. The DNA methylation level of MAPK8IP2 was higher in primary tumors than in normal tissues, and the high MAPK8IP2 DNA methylation group of PCa patients had poor survival. Enrichment analysis indicated that MAPK8IP2 was involved in the MAPK signaling pathway. In vitro, knockdown of MAPK8IP2 inhibited PCa cell proliferation, migration and invasion.

Conclusion

MAPK8IP2 is a potential target for PCa treatment and can serve as a novel biomarker for PCa diagnosis and prognosis evaluation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-10259-2.

The critical role of epigenetic mechanisms involved in nanotoxicology

Over the past decades, nanomaterials (NMs) have been widely applied in the cosmetic, food, engineering, and medical fields. Along with the prevalence of NMs, the toxicological characteristics exhibited by these materials on health and the environment have gradually attracted attentions. A growing number of evidences have indicated that epigenetics holds an essential role in the onset and development of various diseases. NMs could cause epigenetic alterations such as DNA methylation, noncoding RNA (ncRNA) expression, and histone modifications. NMs might alternate either global DNA methylation or the methylation of specific genes to affect the biological function. Abnormal upregulation or downregulation of ncRNAs might also be a potential mechanism for the toxic effects caused by NMs. In parallel, the phosphorylation, acetylation, and methylation of histones also take an important part in the process of NMs-induced toxicity. As the adverse effects of NMs continue to be explored, mechanisms such as chromosomal remodeling, genomic imprinting, and m⁶ A modification are also gradually coming into the limelight. Since the epigenetic alterations often occur in the early development of diseases, thus the relevant studies not only provide insight into the pathogenesis of diseases, but also screen for the prospective biomarkers for early diagnosis and prevention. This review summarizes the epigenetic alterations elicited by NMs, hoping to provide a clue for nanotoxicity studies and security evaluation of NMs. This article is categorized under: Toxicology and Regulatory Issues in Nanomedicine > Toxicology of Nanomaterials.

A pan-cancer analysis of copper homeostasis-related gene lipoyltransferase 1: Its potential biological functions and prognosis values

As a key copper homeostasis-related molecule, lipoyltransferase 1 (LIPT1) is an essential enzyme for the activation of mitochondrial 2-ketoacid dehydrogenase, participating in fatty acylation. However, the biological significances of LIPT1 in the pan-cancer are unclear. Here, we comprehensively analyzed the functional characteristics of LIPT1 in human cancers and its roles in immune response. We found that LIPT1 was down-regulated in some cancers. And LIPT1 overexpression is associated with favorable prognosis in these patients, such as breast cancer, clear cell renal cell carcinoma, ovarian cancer and gastric cancer. We also explored the mutational status and methylation levels of LIPT1 in human cancers. Gene enrichment analysis indicated that abnormally expressed LIPT1 was significantly associated with immune cells infiltration, such as B cells, CD8⁺ T cells and cancer-associated fibroblast cells. The result from single cell sequencing reflected the important roles of LIPT1 in the regulation of several biological behaviors of cancer cells, such as DNA damage response and cell apoptosis. Taken together, our research could provide a comprehensive overview about the significances of LIPT1 in human pan-cancer progression, prognosis and immune.

Prediction of the Effect of Methylation in the Promoter Region of ZP2 Gene on Egg Production in Jinghai Yellow Chickens

Egg production in chickens is a quantitative trait. The aim of this study was to investigate the effect of promoter methylation of the Zona pellucida 2 (ZP2) gene on egg production. Real-time fluorescence quantification showed that the expression of the ZP2 gene in the ovaries of 300-day-old Jinghai yellow chickens in the high-laying group was significantly higher than that in the low-laying group (p < 0.01). A series of deletion fragments of the ZP2 gene promoter in Jinghai yellow chickens had different promoter activities in DF-1 cells, and the core region of the ZP2 gene promoter was found to be between −1552 and −1348. Four CpG islands in the promoter region of the ZP2 gene were detected by software prediction. The overall degree of methylation of the ZP2-1 amplified fragment was negatively correlated with mRNA expression to some extent (R = −0.197); the overall degree of methylation of the ZP2-2 amplified fragment was also negatively correlated with mRNA expression to some extent (R = −0.264), in which the methylation of methylcytosine (mC)-9, mC-20, and mC-21 sites was significantly negatively correlated with mRNA expression (p < 0.05). In addition, the mC-20 and mC-21 sites are located on the Sp1 transcription factor binding site, and it is speculated that these two sites may be the main sites for regulating transcription. In summary, the methylation sites mC-20 and mC-21 of the ZP2 gene may inhibit the binding of Sp1 and DNA, affect the transcription of the ZP2 gene, and then affect the number of eggs produced by the Jinghai yellow chickens.