LSD1 Influences Both Histone and DNA Methylation During Terminal Erythroid Maturation

Blood ◽  
2013 ◽  
Vol 122 (21) ◽  
pp. 3671-3671
Author(s):  
Michael Getman ◽  
Jeffrey Malik ◽  
James Palis ◽  
Laurie A Steiner

Abstract The molecular mechanisms that drive the maturation of a committed erythroid progenitor to a functional red blood cell are incompletely understood. LSD1 (Lysine-Specific Histone Demethylase 1) is a widely expressed histone demethylase that plays an important role in erythroid maturation (Kereyni, elife, 2013). Although LSD1 is important for a number of biologic processes ranging from embryonic development to leukemogenesis, the molecular mechanisms underlying the influence of LSD1 on gene expression are incompletely understood. The goal of our study is to elucidate the molecular mechanisms by which LSD1 regulates erythroid gene expression and influences erythroid maturation. We hypothesize that LSD1 promotes specific patterns of histone and DNA methylation that facilitate gene expression changes necessary for normal erythroid maturation to occur. To address this hypothesis, the functional and molecular consequences of LSD1 knockdown were assessed in Extensively Self Renewing Erythroblasts (ESREs), a non-transformed, karyotypically normal model of terminal erythroid maturation (England, Blood, 2011). Primary fetal liver was cultured in the presence of EPO, SCF, IGF1 and dexamethasone to derive ESREs. The ESREs were capable of extensive ex-vivo expansion, doubling daily at the proerythroblast phase, however when matured, >90% of cells became benzidine positive and >65% enucleated within 3 days. Lentiviral-mediated shRNA was used to knock down LSD1 in expanding ESREs. Imaging flow cytometry done on maturation day 3 demonstrated that the knockdown cells had impairments in multiple facets of maturation, with larger cell and nuclear areas, higher kit expression, and lower rates of enucleation than the scramble control. LSD1 knockdown was also associated with impaired hemoglobin accumulation (78% vs. 95% benzidine positive; p<0.005). Treatment of ESREs with an inhibitor to LSD1 (Tranylcypromine; TCP) resulted in similar abnormalities in cell and nuclear size, kit expression, hemoglobin accumulation, and enucleation (40% vehicle vs.1% TCP). The functional deficits in maturation, including abnormal kit expression and low rates of enucleation, persisted on maturation day 4. To delineate the molecular mechanisms underlying this maturation impairment, RNA-seq was done in LSD1 knockdown and scramble control samples, and 230 differentially expressed genes (FDR<0.01) were identified using cuffdiff (Trapnell, Nat Biotech, 2013). Consistent with LSD1’s role in erythroid maturation, Ingenuity Pathway Analysis identified multiple networks involving hemoglobin synthesis, and GATA1, EPO, and KLF1 were all predicted as upstream regulators (p-values of 8.24e10-11, 7.25 e10-6, and 3.86e10-4, respectively). To better understand how LSD1 influences gene expression, chromatin immunoprecipitation coupled with high throughput sequencing was used to identify sites of H3K4me2 binding in the differentially expressed genes. 214/230 differentially expressed genes were associated with sites of H3K4me2 occupancy. Quantitative ChIP demonstrated that LSD1 inhibition was associated with increases in H3K4me2 levels at a subset of these sites, however consistent with previous studies, global levels of H3K4me2, determined by Enzyme Linked Immunosorbent Assay (ELIZA), did not change significantly. Although it is known that LSD1 demethylates and stabilizes the maintenance DNA methyltransferase DNMT1 (Wang, Nat Genet 2009), the consequences of LSD1 loss on DNA methylation (5-methyl cytosine; 5-mC) have yet to be investigated. To gain a comprehensive understanding of how LSD1 regulates erythroid gene expression, changes in the level of 5-mC were assessed after knockdown or inhibition of LSD1. Global 5-mC levels, determined by ELIZA assay, were ∼30% lower in TCP treated samples than vehicle treated control (p<0.02) and western blot demonstrated a 3-fold decrease in DNMT1 protein in the TCP treated samples. Both methyl binding domain pull-down coupled with quantitative PCR and genome-wide bisulfite sequencing were utilized to assess changes in 5-mC levels in the differentially expressed genes. Loss of LSD1 was associated with significantly lower levels of 5-mC at several differentially expressed, erythroid-specific genes, such as bh1. Taken together, these data support the hypothesis that LSD1 influences both histone and DNA methylation at genes important for erythroid maturation. Disclosures: No relevant conflicts of interest to declare.

Blood ◽  
2012 ◽  
Vol 120 (21) ◽  
pp. 3442-3442
Author(s):  
Michael Getman ◽  
Samantha J England ◽  
James Palis ◽  
Laurie A Steiner

Abstract Abstract 3442 The maturation of a committed erythroid progenitor to a functional red blood cell is a complex process involving significant changes in gene expression during a time of rapid cell division and nuclear condensation. LSD1 (Lysine-Specific Histone Demethylase 1) is a histone H3 lysine 4 (H3K4) and lysine 9 (H3K9) demethylase that plays pivotal role in this process. LSD1 participates in both enhancer and repressor complexes, and facilitates repression of γ-globin by participating in the Direct Repeat Erythroid Complex (Cui, MCB, 2011). LSD1 inhibitors Tranylcypromine (TCP) and Pargyline (PG) are being investigated as potential therapies for the β-globinopathies, however little is known about the broader functional or genomic consequences of LSD1 inhibition on terminal erythroid maturation. Both TCP and PG impair erythroid maturation in Extensively Self Renewing Erythroblasts (ESREs), a primary cell model of terminal erythroid maturation. ESREs are primary cells derived from fetal liver that proliferate extensively in culture, but retain the ability to appropriately mature and enucleate (England, Blood, 2011), making them ideal for functional and genomic studies of terminal erythroid maturation. In untreated or vehicle (DMSO) treated cultures >90% of cells are benzidine positive by day3 of maturation. In contrast, cultures treated with 400um PG, 1um TCP, or 2um TCP were 72, 42, and 33% benzidine positive by maturation day3, respectively. Cells in the TCP-and PG- treated cultures also had morphologic evidence of impaired maturation, with larger nuclei and more basophilic cytoplasm. In addition to its role as a histone demethylase, LSD1 stabilizes DNMT1 (DNA methyltransferase 1; Wang, Nat Genet 2009). We hypothesized that loss of DNA methylation contributes to the maturation impairment seen with LSD1 inhibitors, and that inhibition of DNMTs with decitabine would also impair terminal erythroid maturation. Consistent with this hypothesis, ESREs treated with decitabine demonstrated a dose-dependent impairment of maturation similar to that seen with PG and TCP. To elucidate the molecular mechanisms underlying the maturation impairment in TCP- and PG- treated cultures, levels of H3K4me2 and methylated DNA (5-methyl cytosine, 5-mC) were assessed both globally and at specific loci. An ELISA (Enzyme-linked Immunosorbent Assay) was used to assess global levels of H3K4me2 and 5-mC in vehicle-, PG-, and TCP-treated cultures after 24 hours of maturation. Global levels of H3K4me2 were significantly higher in PG- and TCP- treated samples than control. In maturing cells, there was no significant difference in the level of 5-mC in vehicle- and inhibitor- treated cultures. It is well established, however, that global DNA methylation decreases with erythroid maturation (Seashore, Science, 2011), and a significant decrease in 5-mC occurs in ESREs during the first 24hrs of maturation. As TCP- and vehicle- treated cultures mature differently, the effect of TCP on 5-mC levels was also assessed in self-renewing ESREs at the proerythroblast stage. Unlike maturing cells, TCP-treated proerythroblasts had a significant decrease in 5-mC levels compared to control. Chromatin immunoprecipitation (ChIP) was used to examine the local effects of LSD1 inhibition on H3K4me2 enrichment at erythroid-specific promoters. TCP-treated cultures had non-uniform changes in H3K4me2 enrichment, with levels increased at some promoters (e.g. protein 4.1,εy-globin), but unchanged at others (e.g. β-globin). To further study the relationship between LSD1 inhibition and H3K4me2 levels, ChIP-seq was used to identify LSD1 sites that co-localized with putative enhancers, defined as peaks of H3K4me2 binding > than 1kb from a transcription start site. ChIP-qPCR was used to compare the level of H3K4me2 at 5 validated enhancer-associated LSD1 sites in vehicle- and TCP-treated cells. The effect of TCP was variable, with only 2/5 enhancer-associated LSD1 sites having increased H3K4me2. Lastly, the local effects of inhibitors on 5-mC were examined using a methyl binding domain pulldown coupled with qPCR. In TCP-treated cells, 5-mC levels declined at several loci, most notably at the εy-globin promoter. Taken together, these results suggest that the impaired erythroid maturation associated with LSD1 inhibition results from the perturbation of multiple mechanisms of epigenetic regulation. Disclosures: No relevant conflicts of interest to declare.


2021 ◽  
Vol 14 (1) ◽  
Author(s):  
Xiao-Liang Xing ◽  
Zhi-Yong Yao ◽  
Chaoqun Xing ◽  
Zhi Huang ◽  
Jing Peng ◽  
...  

Abstract Background Colorectal cancer (CRC) is the second most prevalent cancer, as it accounts for approximately 10% of all annually diagnosed cancers. Studies have indicated that DNA methylation is involved in cancer genesis. The purpose of this study was to investigate the relationships among DNA methylation, gene expression and the tumor-immune microenvironment of CRC, and finally, to identify potential key genes related to immune cell infiltration in CRC. Methods In the present study, we used the ChAMP and DESeq2 packages, correlation analyses, and Cox regression analyses to identify immune-related differentially expressed genes (IR-DEGs) that were correlated with aberrant methylation and to construct a risk assessment model. Results Finally, we found that HSPA1A expression and CCRL2 expression were positively and negatively associated with the risk score of CRC, respectively. Patients in the high-risk group were more positively correlated with some types of tumor-infiltrating immune cells, whereas they were negatively correlated with other tumor-infiltrating immune cells. After the patients were regrouped according to the median risk score, we could more effectively distinguish them based on survival outcome, clinicopathological characteristics, specific tumor-immune infiltration status and highly expressed immune-related biomarkers. Conclusion This study suggested that the risk assessment model constructed by pairing immune-related differentially expressed genes correlated with aberrant DNA methylation could predict the outcome of CRC patients and might help to identify those patients who could benefit from antitumor immunotherapy.


2020 ◽  
Author(s):  
Na Li ◽  
Ru-feng Bai ◽  
Chun Li ◽  
Li-hong Dang ◽  
Qiu-xiang Du ◽  
...  

Abstract Background: Muscle trauma frequently occurs in daily life. However, the molecular mechanisms of muscle healing, which partly depend on the extent of the damage, are not well understood. This study aimed to investigate gene expression profiles following mild and severe muscle contusion, and to provide more information about the molecular mechanisms underlying the repair process.Methods: A total of 33 rats were divided randomly into control (n = 3), mild contusion (n = 15), and severe contusion (n = 15) groups; the contusion groups were further divided into five subgroups (1, 3, 24, 48, and 168 h post-injury; n = 3 per subgroup). Then full genome microarray of RNA isolated from muscle tissue was performed to access the gene expression changes during healing process.Results: A total of 2,844 and 2,298 differentially expressed genes were identified in the mild and severe contusion groups, respectively. The analysis of the overlapping differentially expressed genes showed that there are common mechanisms of transcriptomic repair of mild and severe contusion within 48 h post-contusion. This was supported by the results of principal component analysis, hierarchical clustering, and weighted gene co‐expression network analysis of the 1,620 coexpressed genes in mildly and severely contused muscle. From these analyses, we discovered that the gene profiles in functional modules and temporal clusters were similar between the mild and severe contusion groups; moreover, the genes showed time-dependent patterns of expression, which allowed us to identify useful markers of wound age. We then performed an analysis of the functions of genes (including Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway annotation, and protein–protein interaction network analysis) in the functional modules and temporal clusters, and the hub genes in each module–cluster pair were identified. Interestingly, we found that genes downregulated within 24−48 h of the healing process were largely associated with metabolic processes, especially oxidative phosphorylation of reduced nicotinamide adenine dinucleotide phosphate, which has been rarely reported. Conclusions: These results improve our understanding of the molecular mechanisms underlying muscle repair, and provide a basis for further studies of wound age estimation.


Blood ◽  
2011 ◽  
Vol 118 (21) ◽  
pp. 5023-5023
Author(s):  
Monika Belickova ◽  
Jaroslav Cermak ◽  
Jitka Vesela ◽  
Eliska Cechova ◽  
Zuzana Zemanova ◽  
...  

Abstract Abstract 5023 A direct effects of lenalidomide on gene expression in 5q- patients was studied using HumanRef-8 v2 Expression BeadChips (Illumina). Expression profiles of 6 patients (before treatment and at the time of the first erytroid response) and 6 healthy controls were investigated from CD14+ monocytes of peripheral blood. Differentially expressed genes were identified by Significance Analysis of Microarrays (SAM). Simultaneously, selected genes (TNF, JUN, IL1) were monitored in the course of treatment using Real-Time PCR with Taqman Gene Expression Assays. A comparison of gene expression levels before and during lenalidomide treatment revealed 97 differentially expressed genes (FC >2; p<0.05) related to following biological processes: immune response (16 genes), inflammatory response (11 genes), response to bacteria (8 genes), anti-apoptosis (7 genes), regulation of MAP kinase activity (5 genes), oxygen transport (4 genes), and regulation of cell proliferation (11 genes). An overexpression of a number of cytokines (e.g. TNF, IL8, IL1B, CCL3L, CXCL2, and TNFAIP3) was detected in patients before treatment, after lenalidomide administration expression of the majority of the up-regulated cytokine genes decreased to the control baseline level. Detected overproduction of the cytokines in 5q- syndrome may lead to an increased apoptosis of hematopoietic progenitor cells and together with excessive oxidative stress may contribute to the damage the hematopoietic niche. In the same manner, untreated patients showed suppressed expression of two genes (CXCR4, CRTAP) which play an important role in the stem cell niche. After treatment, we detected increased expression of these genes. Both the observations might explain favorable effects of lenalidomide on the bone marrow stroma defect seen in 5q- syndrome. On the other hand, a substantial increase of the ARPC1B gene (an activator and a substrate of Aurora A) expression was detected after lenalidomide treatment. Since overexpression of Aurora A leads to polyploidy and chromosomal instability, ARPC1B might play a role in the disease progression observed in some patients treated with lenalidomide. To conclude, described changes in genes expression may contribute to identification of the pathways affected by lenalidomide and to the explanation of some effects of this drug that have not been fully understood yet. Supported by grants NS/9634 MZCR, UHKT2005 00023736, MSM0021620808 and COST EUGESMA Disclosures: No relevant conflicts of interest to declare.


2018 ◽  
Vol 1 (3) ◽  
Author(s):  
Li Gao ◽  
Yong Jie Yang ◽  
En Qi Li ◽  
Jia Ning Mao

Objective Evidence indicates that physical activity influence bone health. However, the molecular mechanisms mediating the beneficial adaptations to exercise are not well understood. The purpose of this study was to examine the differentially expressed genes in PBMC between athletes and healthy controls, and to analyze the important functional genes and signal pathways that cause increased bone mineral density in athletes, in order to further reveal the molecular mechanisms of exercise promoting bone health. Methods Five professional trampoline athletes and five age-matched untrained college students participated in this study. Used the human expression Microarray V4.0 expression profiling chip to detect differentially expressed genes in the two groups, and performed KEGG Pathway analysis and application of STRING database to construct protein interaction Network; Real-Time PCR technology was used to verify the expression of some differential genes.  Results Compared with healthy controls, there were significant improvement in lumbar spine bone mineral density, and 236 up-regulated as well as 265 down-regulated in serum samples of athletes. The differentially expressed genes involved 28 signal pathways, such as cell adhesion molecules. Protein interaction network showed that MYC was at the core node position. Real-time PCR results showed that the expression levels of CD40 and ITGα6 genes in the athletes were up-regulated compared with the healthy controls, the detection results were consistent with that of the gene chip. Conclusions The findings highlight that long-term high-intensity trampoline training could induce transcriptional changes in PBMC of the athletes. These data suggest that gene expression fingerprints can serve as a powerful research tool to design novel strategies for monitoring exercise. The findings of the study also provide support for the notion that PBMC could be used as a substitute to study exercise training that affects bone health.


2021 ◽  
Author(s):  
Yue Zhao ◽  
Chen Wang ◽  
Wangxia Li ◽  
Bingyu Jin ◽  
Yang Xiang ◽  
...  

Abstract BackgroundThe mobidity and mortality of coronary artery disease (CAD) is increasing year by year. Hence it is urgent to probe into the molecular mechanism of CAD and seek new therapeutic strategies. The purpose of our study was to screen genes associated with the development of CAD by using bioinformatics tools and clinical samples. MethodsMicroarray datasets from the Gene Expression Omnibus (GEO) database of peripheral blood cells (PBLs) were chosen for this study, and candidate differentially expressed microRNAs (DEMs) were screened using the limma and weighted co-expression network analysis (WGCNA) packages in R (v4.0). Subsequently, we construct a competitive endogenous RNAs (ceRNA) network and perform enrichment analysis of genes in the network. Meanwhile, differentially methylated genes (DMGs) in PBLs were identified using the "ChAMP" package in a DNA methylation chip. We then constructed the methylation-associated ceRNA network in CAD. Eventually, the methylation levels of genes and the relationship with the expression of genes in ceRNA were validated in PBLs samples using the Illumina Methylation 850K chip and transcriptome sequencing, while gene expressions were verified by qRT-PCR. And the regulation of DNA methylation on gene expression was verified in the THP-1 cells treated with 5-Aza-2'-deoxycytidine (5-AZA). ResultsA total of 71 differentially expressed miRNAs were screened by both WGCNA and limma. Then the ceRNA network in CAD was constructed with 269 nodes and 705 edges, which were significantly enriched in the chemokine-mediated signaling pathway and so on. Furthermore, from 4354 identified DMGs in a methylation data, 34 methylation-associated differentially expressed genes (DEGs) and 1 differentially expressed lncRNA (DEL) were obtained. After verification of methylation experiments in study population A, three genes were found to have altered methylation consistent with the bioinformatics results. And these genes were correlated in terms of methylation and expression levels. Corresponding with the bioinformatics results, qRT-PCR results in validation set B also showed that the expression of AGPAT4 and FAM169A were significantly lower in CAD. In addition, 5-AZA treatment could increase the expression of AGPAT4 and FAM169A in THP-1 cells. ConclusionsOur study deepens the understanding of the molecular mechanisms underlying the pathogenesis of CAD and provides new ideas for its treatment.


2019 ◽  
Vol 19 (1) ◽  
Author(s):  
Cheng Zhang ◽  
Bingye Zhang ◽  
Di Meng ◽  
Chunlin Ge

Abstract Background The incidence of cholangiocarcinoma (CCA) has risen in recent years, and it has become a significant health burden worldwide. However, the mechanisms underlying tumorigenesis and progression of this disease remain largely unknown. An increasing number of studies have demonstrated crucial biological functions of epigenetic modifications, especially DNA methylation, in CCA. The present study aimed to identify and analyze methylation-regulated differentially expressed genes (MeDEGs) involved in CCA tumorigenesis and progression by bioinformatics analysis. Methods The gene expression profiling dataset (GSE119336) and gene methylation profiling dataset (GSE38860) were obtained from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) and differentially methylated genes (DMGs) were identified using the limma packages of R and GEO2R, respectively. The MeDEGs were obtained by overlapping the DEGs and DMGs. Functional enrichment analyses of these genes were then carried out. Protein–protein interaction (PPI) networks were constructed using STRING and visualized in Cytoscape to determine hub genes. Finally, the results were verified based on The Cancer Genome Atlas (TCGA) database. Results We identified 98 hypermethylated, downregulated genes and 93 hypomethylated, upregulated genes after overlapping the DEGs and DMGs. These genes were mainly enriched in the biological processes of the cell cycle, nuclear division, xenobiotic metabolism, drug catabolism, and negative regulation of proteolysis. The top nine hub genes of the PPI network were F2, AHSG, RRM2, AURKB, CCNA2, TOP2A, BIRC5, PLK1, and ASPM. Moreover, the expression and methylation status of the hub genes were significantly altered in TCGA. Conclusions Our study identified novel methylation-regulated differentially expressed genes (MeDEGs) and explored their related pathways and functions in CCA, which may provide novel insights into a further understanding of methylation-mediated regulatory mechanisms in CCA.


2015 ◽  
Vol 17 (4) ◽  
pp. 373-383 ◽  
Author(s):  
Nada Lukkahatai ◽  
Brian Walitt ◽  
Alexandra Espina ◽  
Dan Wang ◽  
Leorey N. Saligan

Background: Fibromyalgia syndrome (FMS), a chronic musculoskeletal condition characterized by diffuse pain, fatigue, sleep impairment, and cognitive dysfunction, is associated with significant functional disability. Its underlying biological mechanisms are unknown. This study investigated differentially expressed genes between women with FMS and healthy volunteers. Methods: Women who met the 1990 or 2010 American College of Rheumatology fibromyalgia criteria were compared to age- and race-matched pain-free healthy women. Peripheral blood samples were collected, and a full genome microarray gene expression analysis was performed. One-way analysis of variance was used to identify differentially expressed genes using the filtering criterion of 1% false discovery rate. Analysis of canonical pathways associated with these genes was performed. Confirmatory quantitative real-time polymerase chain reaction and enzyme-linked immunosorbent assay verified microarray results. Independent t-tests compared gene and protein expression between groups. Result: Participants were 54 women with FMS and 25 controls. Expression arrays from a subset of women with FMS ( n = 29) and controls ( n = 20) showed upregulation of 12 genes (>1.8-fold change, p < .05) in the FMS sample. Differentially expressed genes were related to B-cell development, primary immunodeficiency signaling, and mitotic roles of polo-like kinase. CENPK and HSP90AA1 were the most differentially expressed genes ( p < .01). Conclusion: Activity of interrelated pathways related to immune response, and homeostasis appears to be relevant to the experience of FMS. Replication and exploration of the relationship between gene expression and symptom severity will help determine clinical relevance of these findings.


Diagnostics ◽  
2020 ◽  
Vol 10 (6) ◽  
pp. 416
Author(s):  
Zeenat Mirza ◽  
Umama A. Abdel-dayem

Endometriosis is characterized by ectopic endometrial tissue implantation, mostly within the peritoneum, and affects women in their reproductive age. Studies have been done to clarify its etiology, but the precise molecular mechanisms and pathophysiology remain unclear. We downloaded genome-wide mRNA expression and clinicopathological data of endometriosis patients and controls from NCBI’s Gene Expression Omnibus, after a systematic search of multiple independent studies comprising 156 endometriosis patients and 118 controls to identify causative genes, risk factors, and potential diagnostic/therapeutic biomarkers. Comprehensive gene expression meta-analysis, pathway analysis, and gene ontology analysis was done using a bioinformatics-based approach. We identified 1590 unique differentially expressed genes (129 upregulated and 1461 downregulated) mapped by IPA as biologically relevant. The top upregulated genes were FOS, EGR1, ZFP36, JUNB, APOD, CST1, GPX3, and PER1, and the top downregulated ones were DIO2, CPM, OLFM4, PALLD, BAG5, TOP2A, PKP4, CDC20B, and SNTN. The most perturbed canonical pathways were mitotic roles of Polo-like kinase, role of Checkpoint kinase proteins in cell cycle checkpoint control, and ATM signaling. Protein–protein interaction analysis showed a strong network association among FOS, EGR1, ZFP36, and JUNB. These findings provide a thorough understanding of the molecular mechanism of endometriosis, identified biomarkers, and represent a step towards the future development of novel diagnostic and therapeutic options.


2020 ◽  
Author(s):  
Shijie Gao ◽  
Guowang Li ◽  
Hao Yu ◽  
Shiyang Yuan ◽  
Wenxiang Li ◽  
...  

Abstract Background DNA methylation is a common epigenetic regulatory way, and it plays a critical role in various human diseases. However, the potential role of how DNA methylation impacts Ewing’s sarcoma (ES) is not clear. This study aimed to explore the regulatory role of DNA methylation in ES. Methods The microarray data of gene expression and methylation were downloaded from Gene Expression Omnibus (GEO) database, and analyzed via GEO2R. Venn analysis was then applied to identify aberrantly methylated differentially expressed genes (DEGs). Subsequently, Function and pathway enrichment analysis was conducted. Protein-protein interaction (PPI) network was constructed. Hub genes were determined. Besides, a connectivity map (CMap) analysis was performed to screen bioactive compounds for ES treatment. Results A total of 135 hypomethylated high expression genes and 523 hypermethylated low expression genes were identified. The hypomethylated high expression genes were enriched in signal transduction and the apoptosis process. Meanwhile, hypermethylated low expression genes were related to DNA replication and transcription regulation. We next determined 10 hub genes through PPI analysis, among them, C3, TF, and TCEB1 might serve as diagnostic and therapeutic targets. Furthermore, CMap analysis revealed 6 chemicals as potential options for ES treatment. Conclusions For the first time, we jointly analyzed gene profiling and methylation data about ES. The introduction of DNA methylation characteristics over DEGs is helpful to understand the pathogenesis of ES. The identified hub aberrantly methylated DEGs and chemicals might provide some novel insights on ES treatment.


Sign in / Sign up

Export Citation Format

Share Document