- Open Access
DNA methylation-based profiling reveals distinct clusters with survival heterogeneity in high-grade serous ovarian cancer
Clinical Epigenetics volume 13, Article number: 190 (2021)
High-grade serous ovarian cancer (HGSOC) is the most common type of epigenetically heterogeneous ovarian cancer. Methylation typing has previously been used in many tumour types but not in HGSOC. Methylation typing in HGSOC may promote the development of personalized care. The present study used DNA methylation data from The Cancer Genome Atlas database and identified four unique methylation subtypes of HGSOC. With the poorest prognosis and high frequency of residual tumours, cluster 4 featured hypermethylation of a panel of genes, which indicates that demethylation agents may be tested in this group and that neoadjuvant chemotherapy may be used to reduce the possibility of residual lesions. Cluster 1 and cluster 2 were significantly associated with metastasis genes and metabolic disorders, respectively. Two feature CpG sites, cg24673765 and cg25574024, were obtained through Cox proportional hazards model analysis of the CpG sites. Based on the methylation level of the two CpG sites, the samples were classified into high- and low-risk groups to identify the prognostic information. Similar results were obtained in the validation set. Taken together, these results explain the epigenetic heterogeneity of HGSOC and provide guidance to clinicians for the prognosis of HGSOC based on DNA methylation sites.
Epithelial ovarian cancer (EOC) has the highest fatality rate of female reproductive cancers , primarily because more than 70% of EOCs are diagnosed at advanced stages and are associated with disseminated intraperitoneal disease . EOC is a heterogeneous disease that includes several subtypes with different clinical and molecular features . The most prevalent histotype of EOC is high-grade serous ovarian cancer (HGSOC), which is associated with a poor prognosis.
Epigenetic alterations have recently emerged as a common hallmark of human cancer [4,5,6]. DNA methylation is a core type of epigenetic alterations and a key epigenetic regulator of gene expression associated with different tumour types . DNA methylation also acts as a biomarker for the prognosis of cancer [8, 9]. Novel methods for tumour classification have been examined in recent years based on genome-wide DNA methylation data [10, 11]. Several DNA methylation pattern analyses have been reported for many cancer types using The Cancer Genome Atlas (TCGA) data [12,13,14,15,16]. There is also increasing interest in the role of DNA methylation in defining molecular subtypes to assist in elucidating the clinical characteristics and prognosis of EOC [17,18,19]. However, previous researches on the methylation of HGSOC are limited to appointed methylation sites or are established based on the relationship between multiple data integrations. A comprehensive and independent profile of the prognostic value of these aberrantly methylated biomarkers in the HGSOC subtype is not clear. Therefore, a comprehensive and independent assessment of DNA methylation is needed to better understand the methylation-based classification of HGSOC subtypes.
In the present population-based study, we aimed to investigate whether clinically relevant and different HGSOC subtypes could be distinguished using genome-wide DNA methylation pattern evaluation. We examined HGSOC classification strategies based on DNA methylation profiles from the TCGA database. Our methylation-based classification system improves the understanding of methylation heterogeneity and may lead to the identification of prognostic signatures and provide good practice guidelines for the clinical treatment of HGSOC patients.
Methylation data downloading and processing
HGSOC methylation data and clinical follow-up data of HGSOC patients were downloaded from TCGA. These data were collected and analysed in accordance with the TCGA Human Subjects Protection and Data Access Policies. The methylation level of each probe was represented by the β-value, which ranges from 0 (unmethylated) to 1 (fully methylated). To obtain more precise classification results, we deleted some methylation loci using the following exclusion criteria: (1) cytosines preceding guanosine sites (CpGs) in sex chromosomes and (2) missing data in more than 70% of the samples. The remaining missing probes were imputed using the k-nearest neighbours . DNA methylation in promoter regions (2 kb upstream to 0.5 kb downstream from transcription start sites) strongly influences gene expression. Therefore, we annotated these sites into promoter regions and reached a final matrix with 21,123 probes. By incorporating patient clinical information and integrating the DNA methylation array data, the ComBat algorithm in the sva R package was used to remove batch effects. From TCGA, we obtained clinical follow-up information for 587 serous ovarian cancer patients and excluded samples that did not meet the following inclusion criteria: (1) pathology limited to serous ovarian cancer, (2) disease was grade 3, (3) available survival time data, (4) primary tumour sample type, (5) available stage information, and (6) available age information. A total of 479 patients were ultimately included. The samples were randomly divided into a training group (n = 240) and a validation group (n = 239). Survival times less than 30 days were generally due to operative complications, and these samples were deleted. Samples without methylation information were also excluded from our study. Finally, there were 233 samples in the training group and 232 samples in the validation group (Fig. 1).
Determination of statistically significant methylation sites with Cox proportional risk regression models
Methylation sites associated with prognosis were regarded as classification characteristics. The methylation profiles and corresponding clinical characteristic data of the training set were used to construct the Cox regression model with the R package “survival”, and the methylation level of each site was analysed according to the following steps: (a) a univariate Cox proportional hazard regression model was constructed with p < 0.05; (b) the significant methylation sites obtained from univariate Cox regression models were used to establish the multivariate Cox proportional risk regression model based on age and stage as the covariants in survival analysis; (c) statistically significant sites (p < 0.05) in multivariate analyses were used as classification features, and the coxph function in survival package R was used.
Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis annotated by the CpG sites
The methylation sites that were statistically significant in multivariate analyses were annotated to the corresponding genes. The “DOSE”, “clusterProfiler” and “enrich plot” R packages were used to complete the GO and KEGG analyses. A maximum p value of 0.05 was chosen to select only significantly enriched categories.
Consensus clustering analysis based on methylation profiles to obtain methylation subtypes associated with HGSOC prognosis
Consensus Clustering Plus [21, 22] is an appropriate method to examine new subclasses of cancer. The present study performed clustering based on the previously obtained methylation sites that were significant in multivariate analyses with the Consensus Cluster Plus algorithm  of R software to determine the methylation subgroups of HGSOC. The similarity distance between samples was calculated by the Euclidean distance method , and κ-means  was used as the clustering algorithm to search for reliable and stable subgroup classifications. After executing Consensus Cluster Plus, consensus clustering results were obtained. The optimal cluster number was determined based on the following criteria: a relatively low-variation coefficient within the cluster, relatively high consistency, and no obviously increased area under the cumulative distribution function curve. The coefficient of variation was calculated with the following formula: CV = (SD/MN) * 100%, in which MN represents the mean of samples and SD represents the standard deviation.
Survival and clinical characteristic analyses
The log-rank test and Kaplan–Meier plots were used to determine the overall survival and significant differences between HGSOC methylation clusters. Associations between clinical characteristics and HGSOC methylation clusters were analysed with the chi-square test. All tests were two-sided, and differences with p < 0.05 were considered significant.
Construction of the prognostic model and verification using fivefold cross-validation and the validation data
The prognosis-related CpG sites were selected to establish the Cox proportional hazards model. Cox regression analysis was performed using multivariate Cox hazards regression model with stepwise method, which was implemented using survival package in R project. A risk score formula for predicting OS was developed based on a linear combination of the expression level multiplied regression coefficient derived from the cox regression model:
(where “Exp” denotes the methylation level of loci and “β” represents the regression coefficient from the multivariate Cox regression model [25, 26]). By utilizing the median risk score as the threshold, the HGSOC patients were stratified into high- and low-risk groups. Survival curves were estimated by the Kaplan–Meier and log-rank method. Fivefold cross-validation (CV) was used on the patients. The same processes were performed on the validation data to verify the stability and applicability of the model. R/Bioconductor tools and R functions were used for all analyses in R version 3.6.1.
Clinicopathological features and the initial screening of DNA methylation loci of HGSOC
Table 1 summarizes the clinicopathological characteristics of the training group and validation group. The mean age of all patients was 59 years. A total of 436 patients (93.8%) had advanced disease, and 58.1% achieved complete remission. 66.4% of the patients underwent optimal cytoreductive surgery. However, the integrity of information for venous invasion and lymphatic invasion was not satisfactory. Clinicopathological parameters were well balanced between the training and the validation groups (Table 1, p > 0.05). A univariate Cox proportional hazard regression model was used to distinguish the 1434 methylation sites that significantly correlated with survival (p < 0.05). Age and stage are universally acknowledged factors associated with HGSOC prognosis. Therefore, age, stage, and 1434 methylation sites were added to the multivariate Cox proportional hazard regression model, which revealed 780 sites that were significantly related to survival (Additional file 2: Table S1).
Enrichment analysis of differentially methylated genes
It is important to understand the molecular mechanisms of genes associated with methylation loci. To examine possible interactions between differentially methylated genes and their possible influence, we analysed the 780 sites that were significantly correlated with survival in the multivariate Cox proportional hazard regression model and examined their co-expressed genes. We performed GO (Fig. 2A–D) and KEGG pathway analyses (Fig. 2E, F). Biological significance (Fig. 2A, B, Additional file 3: Table S2) revealed enrichment in protein acylation, energy derivation by oxidation, and G1/S transition of the mitotic cell cycle. We also performed molecular function analysis and found that mitogen-activated protein kinase (MAPK) activity, growth factor activity, histone binding, and oxidoreductase activity were involved in the network (Fig. 2C, D, Additional file 4: Table S3). Pathways related to MAPK signalling, pyruvate metabolism, and glycolysis were enriched in the KEGG analysis (Fig. 2E, F, Additional file 5: Table S4). The methylation loci and corresponding genes involved in different pathways are shown in Additional file 6: Table S5. The results suggest that the methylation sites examined in this study may affect HGSOC tumorigenesis and development via MAPK signalling, pyruvate metabolism, and glycolysis.
Consensus clustering revealed distinct methylation subtypes associated with HGSOC prognosis
The Consensus Cluster Plus R software package was used to analyse the consensus clustering of 233 samples. The delta area under the curve had an appreciable increase from k = 2 to k = 4, but the increase rate was not obvious for k > 4 (Fig. 3A, B). Therefore, 4 was selected as a suitable cluster number for further analysis in this study. When k = 4, we named the resulting subgroups C1 (28.8%; 67 samples), C2 (22.7%; 53 samples), C3 (38.6%; 90 samples), and C4 (9.9%; 23 samples). By identifying DNA methylation panels that discriminated the four subgroups, we found that the methylation levels of the four subgroups were significantly related to different molecular features (Fig. 3C). C4 was the poorest prognostic group and exhibited hypermethylation of 54 methylation loci, which corresponded to 51 genes (Additional file 7: Table S6). Therefore, the 51 genes composed the hypermethylation panel, in which ANXA7, IGF2, and SLC5A8 harboured two hypermethylation sites in the samples from the C4 subgroup. C1 featured hypomethylation of cg03848675, cg12493906, and cg13055001, which were annotated as FOXF2, MMP26, and PPP1CA, respectively. FOXF2 and MMP26 were primarily related to tumour metastasis and invasion. PPP1CA was associated with the MAPK pathway. C3 exhibited hypermethylation of cg03848675, which was opposite to the pattern observed in C1, and exhibited hypomethylation of cg14290451 (RPL10A). C2 had the best prognosis and featured hypomethylation of cg13791131, cg25574024, cg24673765, and cg27239157, which were annotated as IGF2 (cg13791131, cg25574024), HSPB6 (cg24673765), and MCF2L2 (cg27239157). IGF2 plays a key role in glucose metabolism, HSPB6 is associated with insulin resistance , and MCF2L2 is related to type 1 diabetes  and polycystic ovary syndrome . All of these genes play a role in metabolic disorders. Therefore, the four subgroups classified based on the methylation levels may reflect changes in some molecular genetic features.
Survival and clinical characteristic analyses of methylation subtypes
We combined the four molecular subtypes with the clinical parameters and performed Kaplan–Meier analysis, and the log-rank tests showed that survival curves of these four molecular subtypes were significantly different (Table 2). C4 was obviously more malignant than the other clusters (Table 2) and exhibited hypermethylation at 54 methylation loci (Additional file 7: Table S6). However, the median overall survival time was significantly longer in C2 than the other clusters (64 months in C2 vs. 35 months in C1, 48 months in C3, 24 months in C4, p = 0.0001, Table 2). We also noted that the residual tumour and lymphatic invasion rates were significantly different between the four subgroups (p = 0.027 and 0.031, respectively, Table 2). C4 was highly interrelated with larger residual tumours (≥ 11 mm), and the lymphatic invasion rate was the lowest in C2, which partially explains why C4 had the worst prognosis and C2 had the best prognosis (Table 2). Median age did not differ between the four clusters (p = 0.774, Table 2).
Construction of the prognosis model and verification in validation data
Based on the four subtype classifications, we established a prognosis prediction model (Fig. 4). To optimize this model to contain only the most predictive genes, a stepwise Cox proportional hazards regression model was used, which identified two genes. Subsequently, a risk score was built: risk value = (1.207973 × cg25574024methylation value) − (1.35159 × cg24673765methylation value). The risk score for each patient was calculated using this formula. In the training cohort, patients were assigned into high-risk group and low-risk group based on the optimized risk value (Fig. 4A–C). Kaplan–Meier survival analyses showed that the rate of survival in the high-risk group was significantly higher than that in the low-risk group (Fig. 4D, p < 0.001). Furthermore, we used the validation cohort to validate the model. With the same formula, patients were divided into a high-risk group and low-risk group by their risk score (Fig. 4E–G). Additionally, results of OS showed that patients in the high-risk group had a significantly shorter survival than the counterparts (Fig. 4H, p = 0.014). Furthermore, the results of fivefold cross-validation show that the model established by the method in this paper has good performance in the prediction of prognosis, and the survival rate has significant differences between high-risk group and low-risk group in the fivefold cross-validation (Additional file 1: Fig. S1).
Classifications that are based solely on the pathological features of tissues have limitations. Remodelling of the epigenome is a fundamental factor in tumour prognosis prediction or risk stratification. DNA methylation is one of the most common epigenetic phenomena. It is universally acknowledged that DNA methylation alterations are potential molecular markers for cancer progression [30,31,32,33]. Several studies of ovarian cancer have reported that DNA methylation signatures play a pivotal role in the molecular classification, survival status, and adjuvant chemotherapy response [34,35,36,37,38,39,40,41]. Despite increasing knowledge of the most common type of ovarian cancer, HGSOC, there are no effective methylation-based molecular signatures . Therefore, we performed the present study to better understand the extent and heterogeneity of aberrant DNA methylation in HGSOC.
The present study analysed prognosis-related methylation sites and obtained the corresponding signalling pathways, which provided the research directions for the methylation of HGSOC in future. We identified four methylation subgroups with different prognoses. There were divergent biological characteristics in the four methylation subgroups, which confirmed the heterogeneity of HGSOC and the necessity of cautious classification. Molecular-targeted therapy is gaining traction, and we may identify the underlying mechanisms in different methylation subgroups to help design new strategies in future. We further established a prognosis prediction model based on the significantly methylated sites using multivariate Cox analysis for the convenience of clinical application. Validation was performed to illustrate the reliability of this prediction model. This study identified stable classification with methylation patterns and clinical meaning, which is informative for HGSOC biological characteristics and prognosis.
Molecular mechanistic studies based on bioinformatics analyses are highly significant to cancer research. Previous studies used methylated DNA Immunoprecipitation in A2780 and CaOV3 ovarian cancer cell lines to confirm a panel of six gene promoters that differentiated serous EOC from normal ovarian surface epithelial cells . Zeller et al.  identified loci at 4092 genes that were hypermethylated in chemoresistant A2780/cp70 cells compared to the parental-sensitive A2780 cell line. Further high-throughput DNA methylation profiling of 27 primary epithelial ovarian tumours and 15 ovarian cancer cell lines revealed significant differences in the DNA methylation profiles between ovarian cancer cell lines and tumours, which underscores the need for caution in the use of cell lines as tumour models for epigenetic molecular studies.
Based on the difference between ovarian tumours and cancer cell lines, preference is given to studies that directly examine the methylation profiles of cancer tissues. Notably, ovarian cancer is a heterogeneous disease that includes five major epithelial ovarian tumour subtypes (high- and low-grade serous, endometrioid, mucinous, and clear cell). Previous studies demonstrated that ovarian tumours of different histological types have distinct methylation profiles , which reinforces the need to treat different histotypes of ovarian cancer as separate diseases. For serous EOC, Keita et al.  found that widespread DNA hypermethylation occurs in tumours with low malignant potential (borderline) and significant DNA hypomethylation was observed only in grade 3 serous EOC tumours.
Although it is the most common histological type of EOC, the methylation status of HGSOC has not been studied in detail. Keita et al.  only included 10 cases of HGSOC, which could not be further studied as a separate group. Bodelon et al.  found three methylation subgroups in 61 HGSOCs. However, the small sample size is a limitation of the quality of data generated from these experiments, and the characteristics of methylation subgroups, which are a matter of cardinal importance for targeted therapy, were not elaborated in detail. Reyes et al.  focused on the difference in methylation between HGSOC and normal fallopian tube tissue and between primary and recurrent ovarian cancer. Montavon et al.  examined and compared the methylation patterns of 10 genes in a cohort of 80 primary HGSOC and 12 benign ovarian surface epithelium samples. Dai et al.  primarily profiled DNA methylation of genes in four pathways. Baranova et al.  identified the methylation of the CDH13, HNF1B, PCDH17, and GATA4 genes to distinguish HGSOC from normal samples. In summary, research on the methylation of HGSOC is limited to appointed methylation sites, differences between HGSOC and normal tissue, and differences between recurrent cancer and primary cancer. Previous TCGA studies  established a four-cluster system. However, the classification system was established based on multiple data integrations (mRNA and miRNA expression and DNA methylation) rather than single methylation data. The present study emphasized the importance of DNA methylation and considered DNA methylation as an independent system, which is different from TCGA articles.
Based on TCGA data, our study primarily focused on comprehensive methylation analyses of HGSOC. The corresponding genes of differentially methylated sites that were significant in the multivariate Cox analysis were analysed. Pathway enrichment analyses suggested that differentially methylated genes were primarily enriched in the MAPK signalling pathway, including the CACNB3, MAP3K12, PGF, MAPK13, and CSF1R genes. The MAPK signalling pathway plays a role in the regulation of gene expression, cellular growth, and survival, and it is implicated in most cancers. However, little research focused on the interaction between methylation modification and MAPK signalling pathway proteins. The glycolysis/gluconeogenesis pathway and pyruvate metabolism pathway were also enriched in KEGG analyses, and these pathways play a vital role in cancer progression [47,48,49]. However, the relationship between these metabolic pathways and DNA methylation is not known. DNA methylation may also be a direction in the study of proteoglycan pathways, including the VAV2, CCND1, DDX5, and ITPR2 genes.
The subtypes identified from the methylation profiling data classified HGSOC into four groups. Notably, the methylation levels of the different subgroups reflected different molecular features. C1 was associated with the hypomethylation of cg03848675, cg12493906, and cg13055001, which were annotated as FOXF2, MMP26, and PPP1CA, respectively. FOXF2 is a critical tumour suppressor in gastric carcinogenesis that mediates upregulation of the E3 ligase IRF2BPL to drive ubiquitylation and degradation of β-catenin, which blunts Wnt signalling and suppresses carcinogenesis . The function of FOXC2 in breast cancer is inconsistent. Some researchers reported that FOXC2 suppressed epithelial–mesenchymal transition and multidrug resistance in basal-like breast cancer . However, another study found that FOXF2 promoted the bone metastasis of breast cancer cells . lncRNA ADAMTS9-AS2 decreased tumour progression in ovarian cancer by regulating the miR-182-5p/FOXF2 axis . Matrix metalloproteases (MMPs) play a vital role in cancer metastasis. The immunostaining intensity of MMP-26 increased with ovarian tumour stage , which suggests a role for MMP-26 in ovarian cancer biological function. PPP1CA was associated with activation of the MAPK signalling pathway [55, 56]. Because the hypomethylation loci of the C1 subtype were closely related to tumour metastasis, C1 was defined as the metastasis subgroup. Targeted therapy inhibits metastasis may be more effective in this subgroup than the other subgroups.
The C2 subtype had the best prognosis and showed relative hypomethylation of cg13791131, cg25574024, cg24673765, and cg27239157, which were annotated as IGF2 (cg13791131, cg25574024), HSPB6 (cg24673765), and MCF2L2 (cg27239157), respectively. IGF2 plays a key role in glucose metabolism, HSPB6 might be associated with insulin resistance , and MCF2L2 might be one of the most important marker genes contributing to type 1 diabetes  and polycystic ovary syndrome . All of these diseases are metabolic disorders, and these three genes are associated with metabolism. Therefore, C2 was defined as the metabolism subtype. Future studies will investigate whether this subtype is related to metabolic disorders and examine the application value of metabolic drugs in this subtype.
C3 presented with hypermethylation associated with cg03848675, which was opposite to the patterns observed in C1, and featured hypomethylation of cg14290451(RPL10A). RpL10A stimulates cell proliferation via the insulin signalling pathway . C4 was the poorest prognostic group and exhibited hypermethylation at 54 methylation loci. The hypermethylation of tumour suppressor genes contributes to a more aggressive phenotype. Therefore, C4 was defined as the hypermethylation subtype, which suggests that demethylation agents could be preclinically tested for this group. It will be of great interest to clarify the underlying reasons for these unique subtypes and elucidate the relationship between different subtypes and their level of sensitivity to specific targeted agents. However, care must be taken. Intervening therapeutically to reverse a pattern seen in a cluster may have adverse effects. The adverse effects of therapeutic drugs require structural transformation by pharmaceutics, but this process is a long way off.
Clinical features, including survival outcome, residual tumours, and lymphatic invasion, were markedly different between the four subgroups. The metabolism subtype was associated with no macroscopic residual tumour and a high probability of negative lymphatic invasion status, which could explain the favourable prognosis. Notably, the frequency of residual tumours in the hypermethylation subtype was higher than that in the other subtypes, which indicates that this group could be treated with neoadjuvant chemotherapy to improve the quality of surgery and reduce the possibility of residual lesions.
We further developed a prediction model for prognostication and clinical application. The prognostic model distinguished the training data sets and the validation sets into different prognosis clusters. This clinically promising model may be used to predict the prognosis of HGSOC patients, and follow-up may be strengthened for high-risk patients.
The results of this exploratory study suggest four distinct HGSOC clusters that are distinguishable with DNA methylation profiling and highlight several important genetic characteristics. We showed that the metabolism subtype had a favourable prognosis and that the hypermethylation subtype had the worst prognosis. This result provides a more detailed explanation of HGSOC heterogeneity. Due to the high rates of residual sites, hypermethylation subtype tumours may be treated with neoadjuvant chemotherapy to improve the quality of surgery and reduce the residual rate. The sensitivity to demethylation agents in this subtype should be examined and elucidated. Agents regulating metabolism may be effective for the metabolism subtype, and agents disrupting tumour metastasis may have value for further exploration in the metastasis subtype. Our prediction model provides guidance for clinicians in decisions related to prognosis. Our findings lay the groundwork for an improved understanding of the methylation-based subtypes of HGSOC and provide a useful resource with clinical implications for further studies.
Availability of data and materials
Clinical information and DNA methylation data were retrieved from the TCGA data portal, which is a publicly available database.
Epithelial ovarian cancer
High-grade serous ovarian cancer
The Cancer Genome Atlas
Cytosines preceding guanosine sites
Kyoto Encyclopedia of Genes and Genomes
Mitogen-activated protein kinase
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin. 2018;68(1):7–30.
Barnholtz-Sloan JS, Schwartz AG, Qureshi F, Jacques S, Malone J, Munkarah AR. Ovarian cancer: changes in patterns at diagnosis and relative survival over the last three decades. Am J Obstet Gynecol. 2003;189(4):1120–7.
Bowtell DD. The genesis and evolution of high-grade serous ovarian cancer. Nat Rev Cancer. 2010;10(11):803–8.
Romani M, Pistillo MP, Banelli B. Epigenetic targeting of glioblastoma. Front Oncol. 2018;8:448.
Han TS, Ban HS, Hur K, Cho HS. The epigenetic regulation of HCC metastasis. Int J Mol Sci. 2018;19(12):3978.
Porten SP. Epigenetic alterations in bladder cancer. Curr Urol Rep. 2018;19(12):102.
Jones PA, Baylin SB. The epigenomics of cancer. Cell. 2007;128(4):683–92.
Li Y, Gong Y, Ning X, Peng D, Liu L, He S, et al. Downregulation of CLDN7 due to promoter hypermethylation is associated with human clear cell renal cell carcinoma progression and poor prognosis. J Exp Clin Cancer Res CR. 2018;37(1):276.
Hao X, Luo H, Krawczyk M, Wei W, Wang W, Wang J, et al. DNA methylation markers for diagnosis and prognosis of common cancers. Proc Natl Acad Sci USA. 2017;114(28):7414–9.
Sturm D, Orr BA, Toprak UH, Hovestadt V, Jones DTW, Capper D, et al. New brain tumor entities emerge from molecular classification of CNS-PNETs. Cell. 2016;164(5):1060–72.
Röhrich M, Koelsche C, Schrimpf D, Capper D, Sahm F, Kratz A, et al. Methylation-based classification of benign and malignant peripheral nerve sheath tumors. Acta Neuropathol. 2016;131(6):877–87.
Yang Z, Jones A, Widschwendter M, Teschendorff AE. An integrative pan-cancer-wide analysis of epigenetic enzymes reveals universal patterns of epigenomic deregulation in cancer. Genome Biol. 2015;16(1):140.
Noushmehr H, Weisenberger DJ, Diefes K, Phillips HS, Pujara K, Berman BP, et al. Identification of a CpG island methylator phenotype that defines a distinct subgroup of glioma. Cancer Cell. 2010;17(5):510–22.
Liao S, Hartmaier RJ, McGuire KP, Puhalla SL, Luthra S, Chandran UR, et al. The molecular landscape of premenopausal breast cancer. Breast Cancer Res BCR. 2015;17:104.
Aine M, Sjödahl G, Eriksson P, Veerla S, Lindgren D, Ringnér M, et al. Integrative epigenomic analysis of differential DNA methylation in urothelial carcinoma. Genome Med. 2015;7(1):23.
Dabrowski MJ, Draminski M, Diamanti K, Stepniak K, Mozolewska MA, Teisseyre P, et al. Unveiling new interdependencies between significant DNA methylation sites, gene expression profiles and glioma patients survival. Sci Rep. 2018;8(1):4390.
Shen H, Fridley BL, Song H, Lawrenson K, Cunningham JM, Ramus SJ, et al. Epigenetic analysis leads to identification of HNF1B as a subtype-specific susceptibility gene for ovarian cancer. Nat Commun. 2013;4:1628.
Cicek MS, Koestler DC, Fridley BL, Kalli KR, Armasu SM, Larson MC, et al. Epigenome-wide ovarian cancer analysis identifies a methylation profile differentiating clear-cell histology with epigenetic silencing of the HERG K+ channel. Hum Mol Genet. 2013;22(15):3038–47.
Cancer Genome Atlas Research Network. Integrated genomic analyses of ovarian carcinoma. Nature. 2011;474(7353):609–15.
Troyanskaya O, Cantor M, Sherlock G, Brown P, Hastie T, Tibshirani R, et al. Missing value estimation methods for DNA microarrays. Bioinformatics (Oxford, England). 2001;17(6):520–5.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Wilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics (Oxford, England). 2010;26(12):1572–3.
Ghosh A, Barman S. Application of Euclidean distance measurement and principal component analysis for gene identification. Gene. 2016;583(2):112–20.
Ahlqvist E, Storm P, Käräjämäki A, Martinell M, Dorkhan M, Carlsson A, et al. Novel subgroups of adult-onset diabetes and their association with outcomes: a data-driven cluster analysis of six variables. Lancet Diabetes Endocrinol. 2018;6(5):361–9.
Dave SS, Wright G, Tan B, Rosenwald A, Gascoyne RD, Chan WC, et al. Prediction of survival in follicular lymphoma based on molecular features of tumor-infiltrating immune cells. N Engl J Med. 2004;351(21):2159–69.
Hu J, Xu L, Shou T, Chen Q. Systematic analysis identifies three-lncRNA signature as a potentially prognostic biomarker for lung squamous cell carcinoma using bioinformatics strategy. Transl Lung Cancer Res. 2019;8(5):614–35.
Wang Y, Xu A, Ye J, Kraegen EW, Tse CA, Cooper GJ. Alteration in phosphorylation of P20 is associated with insulin resistance. Diabetes. 2001;50(8):1821–7.
Zhang D, Efendic S, Brismar K, Gu HF. Effects of MCF2L2, ADIPOQ and SOX2 genetic polymorphisms on the development of nephropathy in type 1 Diabetes Mellitus. BMC Med Genet. 2010;11:116.
Zheng Q, Shi Y, Yang Z, Xu X, Wang L, Xue F, et al. Family-based association study of the MCF2L2 gene and polycystic ovary syndrome. Gynecol Obstet Invest. 2009;68(3):171–3.
Dor Y, Cedar H. Principles of DNA methylation and their implications for biology and medicine. Lancet. 2018;392(10149):777–86.
Witt H, Gramatzki D, Hentschel B, Pajtler KW, Felsberg J, Schackert G, et al. DNA methylation-based classification of ependymomas in adulthood: implications for diagnosis and treatment. Neuro Oncol. 2018;20(12):1616–24.
Klutstein M, Nejman D, Greenfield R, Cedar H. DNA methylation in cancer and aging. Can Res. 2016;76(12):3446–50.
Koch A, Joosten SC, Feng Z, de Ruijter TC, Draht MX, Melotte V, et al. Analysis of DNA methylation in cancer: location revisited. Nat Rev Clin Oncol. 2018;15(7):459–66.
Keita M, Wang ZQ, Pelletier JF, Bachvarova M, Plante M, Gregoire J, et al. Global methylation profiling in serous ovarian cancer is indicative for distinct aberrant DNA methylation signatures associated with tumor aggressiveness and disease progression. Gynecol Oncol. 2013;128(2):356–63.
Bodelon C, Killian JK, Sampson JN, Anderson WF, Matsuno R, Brinton LA, et al. Molecular classification of epithelial ovarian cancer based on methylation profiling: evidence for survival heterogeneity. Clin Cancer Res Off J Am Assoc Cancer Res. 2019;25(19):5937–46.
Song F, Li L, Zhang B, Zhao Y, Zheng H, Yang M, et al. Tumor specific methylome in Chinese high-grade serous ovarian cancer characterized by gene expression profile and tumor genotype. Gynecol Oncol. 2020;158(1):178–87.
Dai W, Zeller C, Masrour N, Siddiqui N, Paul J, Brown R. Promoter CpG island methylation of genes in key cancer pathways associates with clinical outcome in high-grade serous ovarian cancer. Clin Cancer Res Off J Am Assoc Cancer Res. 2013;19(20):5788–97.
Montavon C, Gloss BS, Warton K, Barton CA, Statham AL, Scurry JP, et al. Prognostic and diagnostic significance of DNA methylation patterns in high grade serous ovarian cancer. Gynecol Oncol. 2012;124(3):582–8.
Baranova I, Kovarikova H, Laco J, Sedlakova I, Vrbacky F, Kovarik D, et al. Identification of a four-gene methylation biomarker panel in high-grade serous ovarian carcinoma. Clin Chem Lab Med. 2020;58(8):1332–40.
Zhang W, Barger CJ, Link PA, Mhawech-Fauceglia P, Miller A, Akers SN, et al. DNA hypomethylation-mediated activation of Cancer/Testis Antigen 45 (CT45) genes is associated with disease progression and reduced survival in epithelial ovarian cancer. Epigenetics. 2015;10(8):736–48.
Fang F, Cardenas H, Huang H, Jiang G, Perkins SM, Zhang C, et al. Genomic and epigenomic signatures in ovarian cancer associated with resensitization to platinum drugs. Can Res. 2018;78(3):631–44.
Hiss D. Optimizing molecular-targeted therapies in ovarian cancer: the renewed surge of interest in ovarian cancer biomarkers and cell signaling pathways. J Oncol. 2012;2012:737981.
Gloss BS, Patterson KI, Barton CA, Gonzalez M, Scurry JP, Hacker NF, et al. Integrative genome-wide expression and promoter DNA methylation profiling identifies a potential novel panel of ovarian cancer epigenetic biomarkers. Cancer Lett. 2012;318(1):76–85.
Zeller C, Dai W, Steele NL, Siddiq A, Walley AJ, Wilhelm-Benartzi CS, et al. Candidate DNA methylation drivers of acquired cisplatin resistance in ovarian cancer identified by methylome and expression profiling. Oncogene. 2012;31(42):4567–76.
Houshdaran S, Hawley S, Palmer C, Campan M, Olsen MN, Ventura AP, et al. DNA methylation profiles of ovarian epithelial carcinoma tumors and cell lines. PLoS ONE. 2010;5(2):e9359.
Reyes HD, Devor EJ, Warrier A, Newtson AM, Mattson J, Wagner V, et al. Differential DNA methylation in high-grade serous ovarian cancer (HGSOC) is associated with tumor behavior. Sci Rep. 2019;9(1):17996.
Ganapathy-Kanniappan S, Geschwind JF. Tumor glycolysis as a target for cancer therapy: progress and prospects. Mol Cancer. 2013;12:152.
Bergers G, Fendt SM. The metabolism of cancer cells during metastasis. Nat Rev Cancer. 2021;21(3):162–80.
Wang Z, Dong C. Gluconeogenesis in cancer: function and regulation of PEPCK, FBPase, and G6Pase. Trends Cancer. 2019;5(1):30–45.
Higashimori A, Dong Y, Zhang Y, Kang W, Nakatsu G, Ng SSM, et al. Forkhead box F2 suppresses gastric cancer through a novel FOXF2-IRF2BPL-β-catenin signaling axis. Can Res. 2018;78(7):1643–56.
Cai J, Tian AX, Wang QS, Kong PZ, Du X, Li XQ, et al. FOXF2 suppresses the FOXC2-mediated epithelial-mesenchymal transition and multidrug resistance of basal-like breast cancer. Cancer Lett. 2015;367(2):129–37.
Wang S, Li GX, Tan CC, He R, Kang LJ, Lu JT, et al. FOXF2 reprograms breast cancer cells into bone metastasis seeds. Nat Commun. 2019;10(1):2707.
Wang A, Jin C, Li H, Qin Q, Li L. LncRNA ADAMTS9-AS2 regulates ovarian cancer progression by targeting miR-182-5p/FOXF2 signaling pathway. Int J Biol Macromol. 2018;120(Pt B):1705–13.
Ripley D, Tunuguntla R, Susi L, Chegini N. Expression of matrix metalloproteinase-26 and tissue inhibitors of metalloproteinase-3 and -4 in normal ovary and ovarian carcinoma. Int J Gynecol Cancer Off J Int Gynecol Cancer Soc. 2006;16(5):1794–800.
Chen M, Wan L, Zhang J, Zhang J, Mendez L, Clohessy JG, et al. Deregulated PP1α phosphatase activity towards MAPK activation is antagonized by a tumor suppressive failsafe mechanism. Nat Commun. 2018;9(1):159.
Sun H, Ou B, Zhao S, Liu X, Song L, Liu X, et al. USP11 promotes growth and metastasis of colorectal cancer via PPP1CA-mediated activation of ERK/MAPK signaling pathway. EBioMedicine. 2019;48:236–47.
Chaichanit N, Wonglapsuwan M, Chotigeat W. Ribosomal protein L10A and signaling pathway. Gene. 2018;674:170–7.
We thank the patients and families who contributed to this study.
This study was funded by the Natural Science Foundation of Shanghai (Grant Nos. 17ZR1403500 and 20ZR1408900), Science and Technology Commission of Shanghai Municipality (Grant No. 19YF1404300), and National Natural Science Foundation of China (Grant Nos. 82002748 and 81802596 and 82072877).
Ethical approval and consent to participate
Consent for participation for all patients was obtained through The Cancer Genome Atlas Project.
The authors have no competing interests to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The results of fivefold cross-validation in patients of training group.
A multivariate Cox proportional hazard regression model revealed that 780 sites were significantly related to survival.
The top 20 statistically significant biological processes based on corresponding genes derived from the differentially methylated probe (Table S1)
The top 20 statistically significant molecular functions based on corresponding genes derived from the differentially methylated probe (Table S1).
The top 20 statistically significant pathways according to Kyoto Encyclopedia of Genes and Genomes pathway analysis(Table S1).
The methylation loci and corresponding genes involved in MAPK signalling, pyruvate metabolism, glycolysis, and proteoglycan pathways.
The molecular characteristics of cluster 1 to cluster 4.
About this article
Cite this article
Wang, J., Li, J., Chen, R. et al. DNA methylation-based profiling reveals distinct clusters with survival heterogeneity in high-grade serous ovarian cancer. Clin Epigenet 13, 190 (2021). https://doi.org/10.1186/s13148-021-01178-3
- Ovarian cancer
- High-grade serous ovarian cancer
- DNA profiling
- Methylation subtypes