- 1Reproductive Medicine Center, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
- 2Institute of Digestive Disease and Department of Medicine and Therapeutics, State Key Laboratory of Digestive Diseases, Li Ka Shing Institute of Health Sciences, The Chinese University of Hong Kong, Hong Kong, China
- 3Department of Gynecology, First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
Objective: We aim to identify the crucial genes or potential biomarkers associated with uterine fibroids (UFs), which may provide clinicians with evidence about the diagnostic biomarker of UFs and reveal the mechanism of its progression.
Methods: The gene expression and genome-wide DNA methylation profiles were obtained from Gene Expression Omnibus database (GEO). GSE45189, GSE31699, and GSE593 datasets were included. GEO2R and Venn diagrams were used to analyze the differentially expressed genes (DEGs) and extract the hub genes. Gene Ontology (GO) analysis was performed by the online tool Database for Annotation, Visualization, and Integrated Discovery (DAVID). The mRNA and protein expression of hub genes were validated by RT-qPCR, western blot, and immunohistochemistry. The receiver operating characteristic (ROC) curve was used to evaluate the diagnostic value.
Results: We detected 22 DEGs between UFs and normal myometrium, which were enriched in cell maturation, apoptotic process, hypoxia, protein binding, and cytoplasm for cell composition. By finding the intersection of the data between differentially expressed mRNA and DNA methylation profiles, 3 hub genes were identified, including transmembrane 4 L six family member 1 (TM4SF1), TNF superfamily member 10 (TNFSF10), and proteolipid protein 1 (PLP1). PLP1 was validated to be up-regulated significantly in UFs both at mRNA and protein levels. The area under the ROC curve (AUC) of PLP1 was 0.956, with a sensitivity of 79.2% and a specificity of 100%. Conclusion: Overall, our results indicate that PLP1 may be a potential diagnostic biomarker for uterine fibroids.
1 Introduction
Uterine fibroids (UFs) are one of the most common uterine benign neoplasms in women of reproductive age, with a morbidity of 77% (Stewart, 2005), and symptomatic lesions occur in 20%–40% of UFs patients (Leyland et al., 2022). The main clinical symptom includes menorrhagia, abnormal uterine bleeding, infertility, recurrent spontaneous abortion, and other pelvic disorder (Styer and Rueda, 2016; Dolmans et al., 2021). Moreover, UFs are the primary incidents of hysterectomy (Ciarmela et al., 2022) with a quantifiable economic and social burden (Cardozo et al., 2012). Ultrasound is the first-line imaging technique in the evaluation of UFs (Russo et al., 2022). It can provide information about some characteristics of morphology, such as cystic area, echogenicity, borders, and vascularization of the lesion. Nevertheless, it is difficult for clinicians to differentiate the benign myoma in the uterine from malignant leiomyosarcoma accurately. Recently, a novel diagnosis strategy has emerged that integrates the histological features and molecular biomarkers to provide a comprehensive assessment of UFs and determine whether a complete hysterectomy is required (Levy et al., 2013; Trovik et al., 2014; Croce and Chibon, 2021; Machado-Lopez et al., 2021). However, these potential biomarkers still lack reliable clinical utility (Levy et al., 2013), as the sensitivity or specificity of them is less than 75% or 99.6%, respectively (Anderson et al., 2010). Thus, more valuable biomarkers validated for the diagnosis of UFs are desperately required. It may also enable us better understand the mechanism of progression and some important features of UFs.
DNA methylation, one of the epigenetic modifications of DNA in mammalians, refers to the transfer of a methyl group to the fifth carbon of a cytosine residue on the DNA sequence to form 5-methylcytosine (Reik et al., 2001). It occurs in CpG dinucleotides that are clustered frequently in regions of about 1–2 kb in length, called CpG islands, in or near the promoter and first exon regions of genes (Jones, 2012; Schübeler, 2015; Dor and Cedar, 2018). The frequency of CpG in gene regulatory regions is different. It was demonstrated that in leiomyomas, CpG sites were hypomethylated in the distal region of the estrogen receptor-alpha (ER-alpha) promoter combined with the higher ER-alpha mRNA levels (Asada et al., 2008). Besides, the aberrant expression of methyltransferases (Li et al., 2003) and other existence of differently methylated genomic locus in fibroids were also reported to separate the UFs from myometrium (Croce and Chibon, 2015; Braný et al., 2019; Liu et al., 2019; Sato et al., 2019; Maekawa et al., 2022). Based on the specific hypomethylated/hypermethylated genes (Islam et al., 2013; Sato et al., 2016) and the genome-wide DNA methylation profiles of UFs (Navarro et al., 2012; Maekawa et al., 2013), DNA methylation is considered to be the mainstay epigenetic mechanism of UFs. It is involved in the developmental processes of UFs by silencing, switching, and stabilizing genes. Hence, genes associated with DNA methylation may offer us some useful clinical diagnostic biomarkers for UFs. Nevertheless, the hub gene is still unclear.
In the present study, three Gene Expression Omnibus (GEO) datasets were utilized for analyzing the key gene relevant to DNA methylation in UFs. The hub gene was further validated by RT-qPCR, western blot, and immunohistochemistry. Finally, the receiver operating characteristic (ROC) curve was used to evaluate the performance of this biomarker for diagnosing UFs.
2 Methods
2.1 Obtaining the gene expression profiles in UFs
All three gene expression profiles in leiomyoma and normal myometrium tissue (GSE45189, GSE31699, and GSE593) were obtained from the National Center of Biotechnology Information (NCBI) Gene Expression Omnibus (GEO). The retrieval strategy was present with several keywords: leiomyoma, myometrium, gene expression profiling, and DNA methylation genome-wide association study. The inclusion was as follows: 1) a case-control research design; 2) includes UFs and normal myometrium tissue; 3) The original profiles should contain a genome-wide assessment. The exclusion criteria were the following: 1) non-case–control research design; 2) Other tissue. The analysis of the GSE45189 data set was based on 3 frozen UFs and 3 normal myometrium tissue obtained from the uterus with leiomyoma. The GSE31699 data set includes the gene expression profile of 68 UFs and paired normal myometrium tissue. The GSE593 data set included 6 tissue samples for DEGs analysis only. All details of sample information and experiment type are shown in Table 1.
2.2 Identification of DEGs between uterine fibroids and normal myometrium
We analyzed the DEGs between UFs and normal myometrium tissue from the gene expression profiles of GSE45189, GSE593, and GSE31699 datasets respectively. The differential DNA methylation genes were analyzed from the genome-wide DNA methylation profiles in GSE45189 and GSE31699 datasets respectively. All differential genes were identified using the online analysis tool GEO2R (https://www.ncbi.nlm.nih.gov/geo/geo2r/). Benjamin-Hochberg was applied for the control of false discovery rate (FDR), and p < 0.05 was utilized as the database’s cut-off criteria. We draw the Venn diagram by the online tool (http://bioinformatics.psb.ugent.be/webtools/Venn/).
2.3 Gene Ontology analysis
In this study, Gene Ontology (GO) analysis was performed by the online tool, Database for Annotation, Visualization, and Integrated Discovery (DAVID version 2021, https://david.ncifcrf.gov/home.jsp) (Han et al., 2022). The 22 DEGs distracted from all 3 datasets were uploaded to DAVID, and p < 0.05 was identified as the critical threshold for significant enrichment. The GO term included the following three criteria: molecular function (MF), cell composition (CC), and biological process (BP).
2.4 PLP1 methylation analysis
The CpG islands around the PLP1 gene promoter were profiled by the UCSC Genome online tool (https://genome.ucsc.edu/). The DNA methylation data of the PLP1 gene was retrieved from the DiseaseMeth version 2.0 database (http://bio-bigdata.hrbmu.edu.cn/diseasemeth/) (Song et al., 2022). The RNA modification type of PLP1 was identified by m6A-atlas (Tang et al., 2021) (http://rnamd.org/m6a) and m5C-atlas (Ma et al., 2022) (http://rnamd.org/m5c-atlas/index.html). The possible m6a regulator of PLP1 was analyzed by the online tools m6a target (http://m6a2target.canceromics.org/).
2.5 Clinical data
A total of 48 patients of UFs were recruited for this study who underwent myomectomy or hysterectomy with a final histological diagnosis of uterine fibroids in Tongji Hospital from 2018–2020. 14 UFs-free individuals were considered a control group. The slices of normal myometrium tissue were difficult to obtain, especially for UFs patients who underwent myomectomy. We included the patients with single uterine prolapse (8/14) who underwent hysterectomy, or patients with cesarean section scar diverticulum (CDS) who were treated by hysteroscopy (6/14) as a comparable control. The normal myometrium tissue from the slices of the CDS patients was identified by pathologists and only the section of myometrium tissue was included for further IHC analysis. Exclusion criteria for all participants consisted of fibroid degeneration, leiomyosarcoma, adenomyosis, and other gynecologic or pelvic malignant disorders. Any women with complicated diseases, for example, metabolic disorders, hypertension, autoimmune diseases, and treated with hormones before surgery were excluded. The information of all patients was collected from electronic medical records in Tongji Hospital which contains age, myoma location (FIGO), the maximum diameter of fibroids, and previous history of pregnancies and surgery. All procedures of this study were approved by the Ethics Committee of Tongji Medical College, Huazhong University of Science and Technology (2022S068).
2.6 RT-qPCR
Total RNA was extracted from leiomyoma and normal myometrium tissue using RNA-easy Isolation Reagent (Vazyme, R701). Synthesis of cDNA was performed using the PrimeScript™ RT Master Mix (Takara, RR036A). Then, real-time PCR analyses (Vazyme, Q712-02) were carried out in triplicate for each sample. All gene expression was normalized to GAPDH. The expression levels were calculated using the 2-∆∆ Cq method (Livak and Schmittgen, 2001). The PCR primers were listed at supplemental Tabel1.
2.7 Western blot
All leiomyoma and normal myometrium tissue was lysed by RIPA buffer contended with 1% PMSG. Standard western blotting procedures were used (Liu C. et al., 2021). The primary antibody used PLP1 (Abcam, ab254363, 1:2000). Equal loading was confirmed using the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) antibody (CST, 5174S, 1:2000). The appropriate anti-Rabbit HRP-linked secondary antibody (CST, 7074, 1:3000) was used.
2.8 Immunohistochemistry and hematoxylin and eosin stain
The section (4-μm thick) of paraffin-embedded leiomyoma and normal myometrium were deparaffinized and rehydrated using a series of graded xylene and alcohol. All slices used EDTA for antigen retrieval. After 1 h cooled, 10% H2O2 was used to quench endogenous peroxidase activity. Blocking was performed using goat serum for 30 min, RT. The primary antibody used PLP1 (Abcam, ab254363, 1:2000). The HRP labeled anti-Rabbit secondary antibody was used the following day. Finally, the slices were mounted with the coverslips using Permount TM Mounting Medium. And all adjacent slices were stained with hematoxylin and eosin (H&E) based on the basic protocol. The percentage of positive stained was conculcated as follows (Karpathiou et al., 2021): 0 = 0%, 1 = 0–25%, 2 = 26–50%, 3 = 51–75%, 4 = 76–100%. The intensity scoring was conducted as follows: 0 = no staining, 1 = weak, 2 = moderate, 3 = strong. The final scores of all sections were based on multiplying the percentage by intensity. [0] = negative expression; [1–3] = low expression; [4–12] = high expression.
2.9 Statistical analysis
All data were presented as the mean ± SD, and data generated in vitro were compared using Student’s t-tests. We performed χ2 test to explore the relationship between UFs and normal myometrium for categorical data. Receiver operating characteristic (ROC) analysis based on the IHC score of all cases was performed to evaluate the diagnostic value of PLP1. The optimal cutoff value in the ROC curve was set to the value that maximizes the Youden index. Youden’s index was defined as sensitivity + specificity—1. The statistical significance threshold was set at a p-value of <0.05. SPSS v21.0 (IBM, United States) and GraphPad Prism 8.0 (GraphPad, United States) were used for statistical analysis and figures preparation.
3 Results
3.1 Identification of differentially expressed genes between uterine fibroids and normal myometrium
A total of 163 DEGs were identified between UFs and normal myometrium in GSE593 data set by GEO2R analysis, including 58 upregulating genes and 105 downregulating genes. DEGs were analyzed from the gene expression profiles in GSE45189 and GSE31699 datasets respectively. Thereinto, 189 upregulated and 309 downregulated genes were identified by GEO2R analysis in GSE45189 dataset. As for GSE31699 dataset, 2060 DEGs were identified, which included 1129 upregulating genes and 931 downregulated genes. The DEGs of UFs and normal myometrium for each dataset was visualized in the corresponding volcano plots (Figures 1A–C). The DEGs from all datasets were identified by the Venn diagram (Figure 1 D), 22 DEGs were shown in Table 2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1045395/fgene-13-1045395-HTML-r1/image_m/fgene-13-1045395-g001.jpg)
FIGURE 1. DEGs identified by GEO2R. Volcano plot of differentially expressed genes in GSE593 (A), GSE31699 (B), and GSE45189 (C) datasets. Up-regulation and down-regulation genes are marked with red and blue respectively. The criteria for a DEG are |log2FC|>1 and adjusted p-value < 0.05. (D) Venn diagram based on the DEGs from all 3 datasets, 22 genes were extracted. (E) Gene Ontology Enrichment analysis of 22 DEGs. P < 0.05 was identified as the critical threshold for significant enrichment. MF, molecular function, CC, cell composition, BP, and biological process.
3.2 Analysis of Gene Ontology Enrichment
The GO enriched terms were analyzed by DAVID database. The results showed that DEGs between UFs and normal myometrium of all 3 datasets were mainly enriched in cell maturation, regulation of the apoptotic process, cellular response to hypoxia, and response to testosterone in the biological process. Protein binding was the most enrichment term in the molecular function criterion. Considering the cell composition criterion, the results showed that DEGs mainly concentrated on the cytoplasm. All terms of GO analysis are presented in Figure 1E.
3.3 Identify the hub gene
Epigenomic aberrations, especially DNA methylation have been identified as one of the main mechanisms for UFs pathogenesis (Mlodawska et al., 2022). Over the years, literature has reported that aberrant DNA methylation occurs throughout the genome in UFs (Islam et al., 2013; Sato et al., 2014; Sato et al., 2016; Sato et al., 2019; Maekawa et al., 2022), accompanied by mRNA expression discrepancy, demonstrating that aberrant gene expression caused by aberrant DNA methylation plays a key role in the pathogenesis of UFs. The datasets included in this study, GSE31699 and GSE45189, provided data about genome-wide DNA methylation of UFs (Navarro et al., 2012; Maekawa et al., 2013). The differential DNA methylation genes were analyzed by GEO2R. The volcano plots of the hypermethylation/hypomethylation genes were presented in the supplemental figure. To investigate the accurate genes with aberrantly DNA promoter methylation, we drew the Venn diagram as followed (Figure 2A). The differential DNA methylation genes in both GSE31699 and GSE45189 datasets and the 22 DEGs extracted from all 3 datasets were analyzed, and three hub genes were found (TM4SF1, TNFSF10, PLP1). The RT-qPCR was implemented to verify the gene expression. Among them, PLP1 was overexpressed in UFs tissue (Figures 2B–D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1045395/fgene-13-1045395-HTML-r1/image_m/fgene-13-1045395-g002.jpg)
FIGURE 2. Identification of the hub gene. (A) The co-expressive differentially expressed genes in differently mRNA expression profiles (from GSE593, GSE31699, and GSE45189 datasets) and genome-wide DNA methylation genes (from GSE31699, GSE45189 datasets) by Venn diagram. M-GSE31699, M-GSE45189 indicated the differential DNA methylation genes. (B–D) The relative mRNA expression of TM4SF1, TNFSF10, and PLP1 by RT-qPCR. n = 14/groups. ***p < 0.001. ns no statistical significance.
3.4 The expression of PLP1 might be regulated by both DNA methylation and RNA modification
Multiple CpG islands were present in Supplemental Figure S2. Then, we found that the PLP1 transcripts were methylated to varying degrees in uterine carcinosarcoma and uterine corpus endometrial carcinoma (Supplemental Figure S2). Interestingly, we discovered that PLP1 was linked to N6-methyladenosine (m6A) modification. Metagene analysis of m6A indicating modification of PLP1 in 3′UTR gene region in the human embryonic stem cells (ESC). The RNA binding protein and binding region were present in Supplemental Table S2. We curated and analyzed a set of 5 acknowledged m6A regulators of PLP1 (4 readers and 1 writer). The detailed descriptions of m6A regulators were present in Supplemental Table S3.
3.5 Baseline characteristics of the patients
To verify the expression of PLP1 in UFs tissue, we collected tissue samples from UFs patients. A total of 48 UFs patients from Tongji Hospital who underwent abdominal surgery caused by uterine fibroids (UFs) and 14 UFs-free individuals were included in this study. The mean age of all participants enrolled in this study is approximately 45, of which the age of UFs patients is 42 and UFs-free participants is 44. There was no statistical difference in age, number of pregnancies, and whether previous abdominal surgery was performed between these two groups. The basic characteristics of all participants were shown in Table 3. The additional fibroid characteristics of UFs patients were summarized in Table 4. The mean maximum diameter of leiomyoma in this study is 6.4 cm. The mean size of UFs in this study was relatively large because all those fibroids were detected for surgical reasons. Myomectomy was operated on in 79.2% (38/48) of patients while the rest of the patients (20.8%, 10/48) underwent a hysterectomy. Notably, 12.5% (6/48) of patients had a previous myomectomy.
3.6 Elevated PLP1 expression in UFs tissue
The qRT-PCR and Western blot results showed PLP1 was overexpressed in UFs compared with that in normal myometrium (Figure 2B, Figures 3A,B). In addition, IHC staining of PLP1 was performed in UFs and paired normal myometrium (Figure 3C). The corresponding H&E-stained sections were also shown in Figure 3C to illustrate the histological characteristics of UFs. PLP1 was upregulated in UFs compared with that in normal myometrium tissue. The relative IHC score based on all cases demonstrated the expression of PLP1 significantly increased in UFs (Figure 3D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1045395/fgene-13-1045395-HTML-r1/image_m/fgene-13-1045395-g003.jpg)
FIGURE 3. PLP1 expression in UFs tissue (A) The PLP1 protein expression of tissue from fibroids and normal myometrium was assessed by Western blot. (B) The relative expression of PLP1 is based on the gray value of Western blot. N = 4/groups. *p < 0.05. (C) Immunohistochemical and corresponding hematoxylin and eosin stains results of fibroids and normal myometrium tissue (magnification 200). IHC, immunohistochemical, HE, hematoxylin and eosin stains, UFs, Uterine Fibroids, NM, Normal Myometrium. (D) Relative IHC score. UFs, n = 48. NM, n = 14. ***p < 0.001.
3.7 ROC curve analysis
The diagnostic value of PLP1 of uterine fibroids was determined by the ROC curve which was constructed by the IHC score (Figure 4). The area under the ROC curve (AUC) was 0.956 with p < 0.005. The results showed the cutoff value was 2.069 with a sensitivity of 79.2% and a specificity of 100%, suggesting PLP1 presented high diagnostic accuracy of UFs.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1045395/fgene-13-1045395-HTML-r1/image_m/fgene-13-1045395-g004.jpg)
FIGURE 4. ROC curve depicting the diagnostic value of PLP1. The values of AUC, optimum cutoff, sensitivity, and specificity are 0.956 (p < 0.001), 2.069, 79.2% and 100%, respectively.
4 Discussion
Uterine fibroids, formed by the proliferation of smooth muscle cells, are one of the most common benign tumors in women of reproductive age (Styer and Rueda, 2016). Numerous studies have demonstrated the potential biomarker for the diagnosis and surveillance of UFs, but the efficacies were still unclear (Levy et al., 2013). In this study, comprehensive bioinformatics methods were used to verify the biomarker as well as to investigate the possible molecular mechanism underlying the development of UFs.
We analyzed the DEGs of 3 datasets, including GSE593, GSE45189, and GSE31699 datasets, under the same criteria using GEO2R. 22 DEGs were found between UFs and normal myometrium tissue samples. The DEGs were mainly enriched in cell maturation and regulator of apoptotic. It is well known that apoptosis is one of the key regulators of fibroid growth, and the dysregulation of apoptotic pathways may contribute to the development of UFs (Okolo, 2008). Another main enrichment in the biological process was the cellular response to hypoxia. Leiomyoma grows in the hypoxia microenvironment, and such environment may lead to the formation of UFs (Zhou et al., 2011; Tal and Segars, 2014). The hypoxia-inducible factor-1protein was also overexpressed in UFs tissue compared with myometrium (Miyashita-Ishiwata et al., 2022a; b). Response to testosterone was also enriched in our analysis. Fujimoto et al. reported that testosterone increased after treatment with estradiol dipropionate in leiomyoma, while it not occurred in the myometrium, which indicated that testosterone might participate in the biological process of UFs (Fujimoto et al., 1994). Nonetheless, Ke LQ et al. failed to find reliable evidence to prove the effectiveness of danazol, a synthetic isoxazole derivative chemically related to 17-ethinyl testosterone, in UFs in clinical trials (Ke et al., 2009). Although there is still an uncertain conclusion on the response of testosterone in UFs, the aberrant activities of this process might impact the development of UFs, which is consistent with our enrichment analysis. Most of the DEGs were enriched in protein binding for molecular function (19/22) based on our results.
As one of the well-studied epigenomic processes in mammals (Bestor, 2000), DNA methylation was considered as a potential mechanism in the pathology of UFs. Many aberrantly hypermethylation/hypomethylation genes were detected in by genome-wide DNA methylation assays (Li et al., 2003; Yamagata et al., 2009; Navarro et al., 2012; Islam et al., 2013; Maekawa et al., 2013; Carbajo-García et al., 2022), and numerous genes were validated to participate in the developmental progress of the UFs in vitro experiments. SATB homeobox 2 and neuregulin 1 were proved to be the upregulated hypermethylated genes involved in the pathogenesis of uterine leiomyoma by activating the WNT/β-catenin and TGF-β pathways (Sato et al., 2019). Shimeng Liu et al. sorted cells from UFs tissue into stem cell-like cells and revealed that most of the stem cells in UFs were hypermethylated. Meanwhile, tumor growth was suppressed when administered the hypomethylating drug, 5′-Aza (Liu et al., 2020; Liu S. et al., 2021). The methylation condition of mediator complex subunit 12 (MED12), one of the most widely reported somatic mutation genes in UFs (Mäkinen et al., 2011), also could separate the UFs from myometrium on account of the aberrantly clustering molecular pathways based on the MED12 methylation-induced DEGs (Maekawa et al., 2022). All these facts combined with the genome-wide DNA methylation profile of UFs suggested that methylation was the vital epigenetic mechanism in UFs and the significant DEGs between UFs and myometrium might be induced by the local changes of DNA methylation at genome loci. To identify the key genes in this progress, we next combined the DEGs extracted above with the different DNA methylation conditions genes in both GSE45189 and GSE31699 datasets. Consequently, TM4SF1, TNFSF10, and PLP1 were identified, in which only PLP1 was significantly upregulated as verified by RT-qPCR.
PLP1 is the most abundant protein of myelination (Eng et al., 1968; Norton and Poduslo, 1973; Wight, 2017), and the mutation of PLP1 can lead to the X-chromosome-linked leukodystrophy Pelizaeus–Merzbacher disease (Elitt et al., 2020). PLP1 has been widely reported in the formation of the central nervous system while the aberrant expression of PLP1 in various malignant tumors was identified by the bioinformatic analysis (Li et al., 2017). Remarkably, the high level of PLP1 in primary colorectal cancer patients presents poorer overall survival times than those with low expression levels (Han et al., 2020). Although the underlying mechanism is unclear, PLP1 was still considered a potential biomarker in other diseases other than only in nervous system lesions (Khalaf et al., 2022). PLP1 was considered as the hypomethylation and transcriptionally upregulated genes in leiomyoma based on the genome-wide DNA methylation and mRNA expression analysis (Navarro et al., 2012). The results of this present study showed that the PLP1 was over-expressed in UFs tissue based on the comprehensive analysis of irregular methylation genes and DEGs in fibroids. Then the validation was conducted at both mRNA and protein levels based on the tissue samples from leiomyoma patients. To the best of our knowledge, it is the first time to validate the dysregulation of PLP1 in benign tumors. Our result suggested that the hypomethylation of PLP1 might be involved in the pathophysiology of UFs, but further experiments still need to implement. Moreover, overexpressed PLP1 exhibited major oxidative phosphorylation deficits (Wight, 2017) and the down-regulation of oxidative phosphorylation aggravates therapeutically adverse tumor hypoxia (Ashton et al., 2018). According to our result that the 22 DEGs including PLP1 were enriched in hypoxia, we constructed the posits that PLP1 might participate in the hypoxia program to onset the fibroids. However, the underlying mechanisms were speculated through different tissue types, which should be further investigated in UFs.
To evaluate the diagnostic meaning of PLP1 expression in UFs, we drew the ROC curve based on the calculation of the IHC score, indicating that PLP1 expression can be a convincing biomarker for UFs with the AUC, sensitivity, and specificity of 0.956, 79.2%, and 100%, respectively. However, the limitation of this ROC analysis is based on IHC score only, lacking clinical utility. The specificity analysis was limited by the fact that PLP1 expression tissue samples were derived only from UFs and normal myometrium, although it was credible.
5 Conclusion
In summary, we obtained 22 DEGs in UFs via bioinformatical analysis and identified PLP1 as the core gene by the combined analysis of genome-wide DNA methylation profiles and 22 DEGs. The over-expression of PLP1 in UFs tissue was validated in both mRNA and protein levels for the first time. Our findings indicated that PLP1 is a potential diagnostic biomarker of the UFs.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Tongji Medical College, Huazhong University of Science and Technology (No. 2022S068). The patients/participants provided their written informed consent to participate in this study.
Author contributions
LC collected data, organized the structure, prepared the figures, and drafted the manuscript. ZL and JL prepared the figures. ZL and SL checked the manuscript. RW and FR collected data. HZ, ZL, and SL participated in the discussion. All authors approved the final version to be published.
Funding
This work was supported by the Foundation of Tongji Hospital (No. 2020JZKT469) and Excellent Youth Funding of Henan Provincial Foundation Committee (No. 222300420091).
Acknowledgments
We appreciate Shitong Lin for his suggestions on the article.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary materials
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.1045395/full#supplementary-material
References
Anderson, G. L., Mcintosh, M., Wu, L., Barnett, M., Goodman, G., Thorpe, J. D., et al. (2010). Assessing lead time of selected ovarian cancer biomarkers: A nested case-control study. J. Natl. Cancer Inst. 102, 26–38. doi:10.1093/jnci/djp438
Asada, H., Yamagata, Y., Taketani, T., Matsuoka, A., Tamura, H., Hattori, N., et al. (2008). Potential link between estrogen receptor-alpha gene hypomethylation and uterine fibroid formation. Mol. Hum. Reprod. 14, 539–545. doi:10.1093/molehr/gan045
Ashton, T. M., Mckenna, W. G., Kunz-Schughart, L. A., and Higgins, G. S. (2018). Oxidative phosphorylation as an emerging target in cancer therapy. Clin. Cancer Res. 24, 2482–2490. doi:10.1158/1078-0432.CCR-17-3070
Bestor, T. H. (2000). The DNA methyltransferases of mammals. Hum. Mol. Genet. 9, 2395–2402. doi:10.1093/hmg/9.16.2395
Braný, D., Dvorská, D., Grendár, M., Ňachajová, M., Szépe, P., Lasabová, Z., et al. (2019). Different methylation levels in the KLF4, ATF3 and DLEC1 genes in the myometrium and in corpus uteri mesenchymal tumours as assessed by MS-HRM. Pathol. Res. Pract. 215, 152465. doi:10.1016/j.prp.2019.152465
Carbajo-García, M. C., Corachán, A., Juárez-Barber, E., Monleón, J., Payá, V., Trelis, A., et al. (2022). Integrative analysis of the DNA methylome and transcriptome in uterine leiomyoma shows altered regulation of genes involved in metabolism, proliferation, extracellular matrix, and vesicles. J. Pathol. 257, 663–673. doi:10.1002/path.5920
Cardozo, E. R., Clark, A. D., Banks, N. K., Henne, M. B., Stegmann, B. J., and Segars, J. H. (2012). The estimated annual cost of uterine leiomyomata in the United States. Am. J. Obstet. Gynecol. 206, 211.e211–e9. doi:10.1016/j.ajog.2011.12.002
Ciarmela, P., Delli Carpini, G., Greco, S., Zannotti, A., Montik, N., Giannella, L., et al. (2022). Uterine fibroid vascularization: From morphological evidence to clinical implications. Reprod. Biomed. Online 44, 281–294. doi:10.1016/j.rbmo.2021.09.005
Croce, S., and Chibon, F. (2015). MED12 and uterine smooth muscle oncogenesis: State of the art and perspectives. Eur. J. Cancer 51, 1603–1610. doi:10.1016/j.ejca.2015.04.023
Croce, S., and Chibon, F. (2021). Molecular prognostication of uterine smooth muscle neoplasms: From CGH array to CINSARC signature and beyond. Genes Chromosom. Cancer 60, 129–137. doi:10.1002/gcc.22906
Dolmans, M. M., Cacciottola, L., Donnez, J., and Poirot, C. (2021). Fertility preservation: How to preserve ovarian function in children, adolescents and adults. J. Clin. Med. 10, 5247. doi:10.3390/jcm10225247
Dor, Y., and Cedar, H. (2018). Principles of DNA methylation and their implications for biology and medicine. Lancet 392, 777–786. doi:10.1016/S0140-6736(18)31268-6
Elitt, M. S., Barbar, L., Shick, H. E., Powers, B. E., Maeno-Hikichi, Y., Madhavan, M., et al. (2020). Suppression of proteolipid protein rescues Pelizaeus-Merzbacher disease. Nature 585, 397–403. doi:10.1038/s41586-020-2494-3
Eng, L. F., Chao, F. C., Gerstl, B., Pratt, D., and Tavaststjerna, M. G. (1968). The maturation of human white matter myelin. Fractionation of the myelin membrane proteins. Biochemistry 7, 4455–4465. doi:10.1021/bi00852a042
Fujimoto, J., Nishigaki, M., Hori, M., Ichigo, S., Itoh, T., and Tamaya, T. (1994). The effect of estrogen and androgen on androgen receptors and mRNA levels in uterine leiomyoma, myometrium and endometrium of human subjects. J. Steroid Biochem. Mol. Biol. 50, 137–143. doi:10.1016/0960-0760(94)90020-5
Han, J., Zhang, X., Liu, Y., Jing, L., Liu, Y. B., and Feng, L. (2020). CLCA4 and MS4A12 as the significant gene biomarkers of primary colorectal cancer. Biosci. Rep. 40, BSR20200963. doi:10.1042/BSR20200963
Han, X., Yang, Y., Qi, J., Zhang, M., Xue, Y., Chu, X., et al. (2022). Protective effects and possible mechanism of 6-gingerol against arsenic trioxide-induced nephrotoxicity based on network pharmacological analysis and experimental validation. Int. Immunopharmacol. 110, 108926. doi:10.1016/j.intimp.2022.108926
Islam, M. S., Protic, O., Stortoni, P., Grechi, G., Lamanna, P., Petraglia, F., et al. (2013). Complex networks of multiple factors in the pathogenesis of uterine leiomyoma. Fertil. Steril. 100, 178–193. doi:10.1016/j.fertnstert.2013.03.007
Jones, P. A. (2012). Functions of DNA methylation: Islands, start sites, gene bodies and beyond. Nat. Rev. Genet. 13, 484–492. doi:10.1038/nrg3230
Karpathiou, G., Chauleur, C., Papoudou-Bai, A., Dagher, S., and Peoc'h, M. (2021). CD56 is expressed in uterine smooth muscle tumors. Int. J. Gynecol. Pathol. 40, 315–323. doi:10.1097/PGP.0000000000000696
Ke, L. Q., Yang, K., Li, J., and Li, C. M. (2009). Danazol for uterine fibroids. Cochrane Database Syst. Rev. 2009, CD007692. doi:10.1002/14651858.CD007692.pub2
Khalaf, G., Mattern, C., Begou, M., Boespflug-Tanguy, O., Massaad, C., and Massaad-Massade, L. (2022). Mutation of proteolipid protein 1 gene: From severe hypomyelinating leukodystrophy to inherited spastic paraplegia. Biomedicines 10, 1709. doi:10.3390/biomedicines10071709
Levy, G., Hill, M. J., Plowden, T. C., Catherino, W. H., and Armstrong, A. Y. (2013). Biomarkers in uterine leiomyoma. Fertil. Steril. 99, 1146–1152. doi:10.1016/j.fertnstert.2012.10.048
Leyland, N., Leonardi, M., Murji, A., Singh, S. S., Al-Hendy, A., and Bradley, L. (2022). A call-to-action for clinicians to implement evidence-based best practices when caring for women with uterine fibroids. Reprod. Sci. 29, 1188–1196. doi:10.1007/s43032-022-00877-3
Li, M., Sun, Q., and Wang, X. (2017). Transcriptional landscape of human cancers. Oncotarget 8, 34534–34551. doi:10.18632/oncotarget.15837
Li, S., Chiang, T. C., Richard-Davis, G., Barrett, J. C., and Mclachlan, J. A. (2003). DNA hypomethylation and imbalanced expression of DNA methyltransferases (DNMT1, 3A, and 3B) in human uterine leiomyoma. Gynecol. Oncol. 90, 123–130. doi:10.1016/s0090-8258(03)00194-x
Liu, C., Wang, M., Zhang, H., and Sui, C. (2021). Altered microRNA profiles of extracellular vesicles secreted by endometrial cells from women with recurrent implantation failure. Reprod. Sci. 28, 1945–1955. doi:10.1007/s43032-020-00440-y
Liu, S., Yin, P., Kujawa, S. A., Coon, J. S. T., Okeigwe, I., and Bulun, S. E. (2019). Progesterone receptor integrates the effects of mutated MED12 and altered DNA methylation to stimulate RANKL expression and stem cell proliferation in uterine leiomyoma. Oncogene 38, 2722–2735. doi:10.1038/s41388-018-0612-6
Liu, S., Yin, P., Xu, J., Dotts, A. J., Kujawa, S. A., Coon, V. J., et al. (2021). Progesterone receptor-DNA methylation crosstalk regulates depletion of uterine leiomyoma stem cells: A potential therapeutic target. Stem Cell Rep. 16, 2099–2106. doi:10.1016/j.stemcr.2021.07.013
Liu, S., Yin, P., Xu, J., Dotts, A. J., Kujawa, S. A., Coon, V. J., et al. (2020). Targeting DNA methylation depletes uterine leiomyoma stem cell-enriched population by stimulating their differentiation. Endocrinology 161, bqaa143. doi:10.1210/endocr/bqaa143
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25, 402–408. doi:10.1006/meth.2001.1262
Ma, J., Song, B., Wei, Z., Huang, D., Zhang, Y., Su, J., et al. (2022). m5C-Atlas: a comprehensive database for decoding and annotating the 5-methylcytosine (m5C) epitranscriptome. Nucleic Acids Res. 50, D196–D203. doi:10.1093/nar/gkab1075
Machado-Lopez, A., Simon, C., and Mas, A. (2021). Molecular and cellular insights into the development of uterine fibroids. Int. J. Mol. Sci. 22, 8483. doi:10.3390/ijms22168483
Maekawa, R., Sato, S., Tamehisa, T., Sakai, T., Kajimura, T., Sueoka, K., et al. (2022). Different DNA methylome, transcriptome and histological features in uterine fibroids with and without MED12 mutations. Sci. Rep. 12, 8912. doi:10.1038/s41598-022-12899-7
Maekawa, R., Sato, S., Yamagata, Y., Asada, H., Tamura, I., Lee, L., et al. (2013). Genome-wide DNA methylation analysis reveals a potential mechanism for the pathogenesis and development of uterine leiomyomas. PLoS One 8, e66632. doi:10.1371/journal.pone.0066632
Mäkinen, N., Mehine, M., Tolvanen, J., Kaasinen, E., Li, Y., Lehtonen, H. J., et al. (2011). MED12, the mediator complex subunit 12 gene, is mutated at high frequency in uterine leiomyomas. Science 334, 252–255. doi:10.1126/science.1208930
Miyashita-Ishiwata, M., El Sabeh, M., Reschke, L. D., Afrin, S., and Borahay, M. A. (2022a). Differential response to hypoxia in leiomyoma and myometrial cells. Life Sci. 290, 120238. doi:10.1016/j.lfs.2021.120238
Miyashita-Ishiwata, M., El Sabeh, M., Reschke, L. D., Afrin, S., and Borahay, M. A. (2022b). Hypoxia induces proliferation via NOX4-Mediated oxidative stress and TGF-β3 signaling in uterine leiomyoma cells. Free Radic. Res. 56, 163–172. doi:10.1080/10715762.2022.2061967
Mlodawska, O. W., Saini, P., Parker, J. B., Wei, J. J., Bulun, S. E., Simon, M. A., et al. (2022). Epigenomic and enhancer dysregulation in uterine leiomyomas. Hum. Reprod. Update 28, 518–547. doi:10.1093/humupd/dmac008
Navarro, A., Yin, P., Monsivais, D., Lin, S. M., Du, P., Wei, J. J., et al. (2012). Genome-wide DNA methylation indicates silencing of tumor suppressor genes in uterine leiomyoma. PLoS One 7, e33284. doi:10.1371/journal.pone.0033284
Norton, W. T., and Poduslo, S. E. (1973). Myelination in rat brain: Changes in myelin composition during brain maturation. J. Neurochem. 21, 759–773. doi:10.1111/j.1471-4159.1973.tb07520.x
Okolo, S. (2008). Incidence, aetiology and epidemiology of uterine fibroids. Best. Pract. Res. Clin. Obstet. Gynaecol. 22, 571–588. doi:10.1016/j.bpobgyn.2008.04.002
Reik, W., Dean, W., and Walter, J. (2001). Epigenetic reprogramming in mammalian development. Science 293, 1089–1093. doi:10.1126/science.1063443
Russo, C., Camilli, S., Martire, F. G., Di Giovanni, A., Lazzeri, L., Malzoni, M., et al. (2022). Ultrasound features of highly vascularized uterine myomas (uterine smooth muscle tumors) and correlation with histopathology. Ultrasound Obstet. Gynecol. 60, 269–276. doi:10.1002/uog.24855
Sato, S., Maekawa, R., Tamura, I., Shirafuta, Y., Shinagawa, M., Asada, H., et al. (2019). SATB2 and NGR1: Potential upstream regulatory factors in uterine leiomyomas. J. Assist. Reprod. Genet. 36, 2385–2397. doi:10.1007/s10815-019-01582-y
Sato, S., Maekawa, R., Yamagata, Y., Asada, H., Tamura, I., Lee, L., et al. (2014). Potential mechanisms of aberrant DNA hypomethylation on the x chromosome in uterine leiomyomas. J. Reprod. Dev. 60, 47–54. doi:10.1262/jrd.2013-095
Sato, S., Maekawa, R., Yamagata, Y., Tamura, I., Lee, L., Okada, M., et al. (2016). Identification of uterine leiomyoma-specific marker genes based on DNA methylation and their clinical application. Sci. Rep. 6, 30652. doi:10.1038/srep30652
Schübeler, D. (2015). Function and information content of DNA methylation. Nature 517, 321–326. doi:10.1038/nature14192
Song, B., Wang, X., Liang, Z., Ma, J., Huang, D., Wang, Y., et al. (2022). RMDisease V2.0: An updated database of genetic variants that affect RNA modifications with disease and trait implication. Nucleic Acids Res. 50, 10290–10310. doi:10.1093/nar/gkac830
Stewart, C. L. W. E. A., and Stewart, E. A. (2005). Uterine fibroids: The elephant in the room. science 308, 1589–1592. doi:10.1126/science.1112063
Styer, A. K., and Rueda, B. R. (2016). The epidemiology and Genetics of uterine leiomyoma. Best. Pract. Res. Clin. Obstet. Gynaecol. 34, 3–12. doi:10.1016/j.bpobgyn.2015.11.018
Tal, R., and Segars, J. H. (2014). The role of angiogenic factors in fibroid pathogenesis: Potential implications for future therapy. Hum. Reprod. Update 20, 194–216. doi:10.1093/humupd/dmt042
Tang, Y., Chen, K., Song, B., Ma, J., Wu, X., Xu, Q., et al. (2021). m6A-Atlas: a comprehensive knowledgebase for unraveling the N6-methyladenosine (m6A) epitranscriptome. Nucleic Acids Res. 49, D134–D143. doi:10.1093/nar/gkaa692
Trovik, J., Salvesen, H. B., Cuppens, T., Amant, F., and Staff, A. C. (2014). Growth differentiation factor-15 as biomarker in uterine sarcomas. Int. J. Gynecol. Cancer 24, 252–259. doi:10.1097/IGC.0000000000000037
Wight, P. A. (2017). Effects of intron 1 sequences on human PLP1 expression: Implications for PLP1-related disorders. ASN Neuro 9, 1759091417720583. doi:10.1177/1759091417720583
Yamagata, Y., Maekawa, R., Asada, H., Taketani, T., Tamura, I., Tamura, H., et al. (2009). Aberrant DNA methylation status in human uterine leiomyoma. Mol. Hum. Reprod. 15, 259–267. doi:10.1093/molehr/gap010
Keywords: uterine fibroids, bioinformatics analysis, DNA methylation, PLP1, biomarker
Citation: Cai L, Liao Z, Li S, Wu R, Li J, Ren F and Zhang H (2022) PLP1 may serve as a potential diagnostic biomarker of uterine fibroids. Front. Genet. 13:1045395. doi: 10.3389/fgene.2022.1045395
Received: 15 September 2022; Accepted: 14 October 2022;
Published: 31 October 2022.
Edited by:
Kunqi Chen, Fujian Medical University, ChinaReviewed by:
Qing Long, Kunming Medical University, ChinaHaolong Li, Peking Union Medical College Hospital, China
Ahmed Hammad, Egyptian Atomic Energy Authority, Egypt
Malihe Rastegarpanah, Breast Cancer Research Center, Motamed Cancer Institute, Iran
Haoran Shi, University of Giessen, Germany
Copyright © 2022 Cai, Liao, Li, Wu, Li, Ren and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hanwang Zhang, hwzhang605@126.com; Fang Ren, fccrenf@zzu.edu.cn