- Open Access
CAMK2N1/RUNX3 methylation is an independent prognostic biomarker for progression-free and overall survival of platinum-sensitive epithelial ovarian cancer patients
Clinical Epigenetics volume 13, Article number: 15 (2021)
To date, no predictive or prognostic molecular biomarkers except BRCA mutations are clinically established for epithelial ovarian cancer (EOC) despite being the deadliest gynecological malignancy. Aim of this biomarker study was the analysis of DNA methylation biomarkers for their prognostic value independent from clinical variables in a heterogeneous cohort of 203 EOC patients from two university medical centers.
The marker combination CAMK2N1/RUNX3 exhibited a significant prognostic value for progression-free (PFS) and overall survival (OS) of sporadic platinum-sensitive EOC (n = 188) both in univariate Kaplan–Meier (LogRank p < 0.05) and multivariate Cox regression analysis (p < 0.05; hazard ratio HR = 1.587). KRT86 methylation showed a prognostic value only in univariate analysis because of an association with FIGO staging (Fisher’s exact test p < 0.01). Thus, it may represent a marker for EOC staging. Dichotomous prognostic values were observed for KATNAL2 methylation depending on BRCA aberrations. KATNAL2 methylation exhibited a negative prognostic value for PFS in sporadic EOC patients without BRCA1 methylation (HR 1.591, p = 0.012) but positive prognostic value in sporadic EOC with BRCA1 methylation (HR 0.332, p = 0.04) or BRCA-mutated EOC (HR 0.620, n.s.).
The retrospective analysis of 188 sporadic platinum-sensitive EOC proved an independent prognostic value of the methylation marker combination CAMK2N1/RUNX3 for PFS and OS. If validated prospectively this combination may identify EOC patients with worse prognosis after standard therapy potentially benefiting from intensive follow-up, maintenance therapies or inclusion in therapeutic studies. The dichotomous prognostic value of KATNAL2 should be validated in larger sample sets of EOC.
Despite all medical efforts and advances, epithelial ovarian cancer (EOC) remains the most lethal gynecological malignancy. In Germany, this disease accounts for 3.2% of all newly diagnosed cancer cases and 5.3% of all cancer-related death of women . The lifetime risk of developing this cancer entity lies by 1–2%. There are several factors increasing or decreasing the lifetime risk. Demographic, environmental and hormonal facts possibly add to the risk. Additionally, various familiar factors influence the chances of developing EOC. Family cancer syndromes predispose patients for the development of this malignancy. A number of syndromes are associated with ovarian cancer, the hereditary breast and ovarian cancer (HBOC) syndrome being the most prominent one . Within this syndrome, the genes BRCA1 and BRCA2 are frequently mutated [3, 4]. Germline mutations account for 15–25% of all ovarian carcinoma [5, 6] while somatic mutations of the BRCA genes are additionally detected in 5–6% of EOC [7, 8]. The five-year survival rate of EOC lies between 35 and 45% [1, 9,10,11]. These low survival numbers can be explained by the lack of early detection methods resulting in the frequent diagnosis of progressed tumor stages (≥ FIGO III). Standard therapy comprises surgery to reduce the tumor burden combined with platin-taxane-based chemotherapy. A successful surgical intervention is crucial for improved survival. A macroscopic completely resected tumor is most preferable and improves the odds of survival and therapy success [12, 13]. Whereas approximately 20% of patients exhibit intrinsic resistant EOC most tumors of initially sensitive patients develop platin resistance and relapse within a median of 18 month [14, 15]. Neither prognostic markers for EOC nor predictive markers for standard chemotherapy are clinically established. Patients with BRCA mutations tend to respond better to platin and show improved progression-free survival but this does not necessarily translate into improved long-term overall survival [16, 17]. Additionally, BRCA mutations or homologous repair defects (HRD) can identify patients benefitting from PARPi treatment [18,19,20,21]. However, even patients with HRD-proficient tumors can be successfully treated with PARPi [19, 22] pointing to the potential requirement to establish better predictive marker. Moreover, no predictive marker is defined for anti-angiogenesis treatment with Bevacizumab. Thus, in the absence of predictive biomarker for any adjuvant treatment, prognostic biomarker may help to identify patients with worse outcome benefitting from maintenance treatment with these new treatment options. Both clinical parameters and different molecules or cell populations are targets of intense research to identify and validate prognostic marker. Clinical parameters can be analyzed posttreatment only (resection status, tumor stage, chemotherapeutic response) but are well established . Biomarker may be analyzed pre-therapeutically (blood-based liquid biopsy) or at least before start of adjuvant chemotherapy (tissue-based biomarker) but are far from being implemented in clinical routine. Biomarker classes for liquid biopsies can be circulating nucleic acids, immune parameters, serum proteins or extracellular vesicle carrying miRNA and other molecules [24,25,26,27,28]. Tissue-based marker may consist of aberrant DNA, coding- or non-coding RNA expression or the analysis of intratumoral immune cell populations [7, 29,30,31,32,33,34]. One DNA aberration which can be used as molecular biomarker is CpG hypermethylation belonging to the field of epigenetics.
Epigenetic changes influence gene expression in a reversible manner without changing the DNA sequence. Histone modification, DNA methylation and, micro-RNA associated silencing are modes of epigenetic alteration in both normal developmental processes and disease, e.g., cancer development [35, 36]. DNA methylation is a stable epigenetic modification catalyzed by DNA-Methyltransferases (DNMT’s) and recognized by methyl-CpG binding proteins (MBDs) resulting in changed chromatin states. Since 5-methyl-cytosines and the DNA itself are highly stable, efforts to use these as biomarker are the focus of current medical research [37,38,39]. Epigenetic changes are seen in various tumor entities, ovarian cancer being one of them . Currently, aberrant DNA methylation in EOC is detected using different approaches such as target-specific methylation-specific PCR (MSP)  or whole genome approaches [42,43,44]. A large number of genes are affected by DNA methylation in EOC . For instance, tumor suppressor genes such as BRCA1, RASSF1A  or MLH1 [41, 46]; cell adhesion genes as ICAM-1  and CDH1  and DNA repair genes PALB2  are affected by DNA methylation. BRCA1 promoter methylation is reported in 10–15% of sporadic OvCa [7, 50,51,52,53]. It appears that a reversal of this phenotype can take place over time leading to changed methylation level in relapsed disease . MLH1 methylation is potentially linked to chemotherapy response or resistance development . However, lack of candidate biomarker validation and low sensitivity or specificity prevents the clinical use of methylation-based biomarker .
In previous studies, we screened twelve primary high-grade serous ovarian carcinomas (HGSOC) using genome wide CpG microarrays for differentially methylated genes to identify potential prognostic biomarkers. 37 hypo-/hypermethylated regions were selected to verify promising candidates in a comparative cohort analysis of 36 samples (PFS < 3 years vs. > 3 years). The methylation marker genes ATL1, CAMK2N1, KATNAL2, KRT86 and RUNX3 were detected as most promising due to good discriminatory power and the combination CAMK2N1 and RUNX3 could identify EOC patients with significantly shortened PFS in univariate analysis . However, the small number of patients and missing multivariate analyzes for the prognostic value related to PFS and OS were limitations of this study. Nevertheless, subsequent functional analyses identified tumor suppressive functions for CAMK2N1 and a dichotomous role for RUNX3 transcript variants regulated by DNA methylation . Robust validation and proof of an independent prognostic value are crucial before biomarkers are tested prospectively in a clinical setting. Therefore, this study aimed to validate the previously identified and most promising epigenetic candidate genes—ATL1, CAMK2N1, KATNAL2, KRT86 and RUNX3—using a larger patient cohort from two different University hospitals enabling multivariate analyses. To provide further information for patients with a hereditary form of ovarian cancer, a small number of such patients were also included and the BRCA1 methylation status was additionally evaluated in all samples.
Methylation analysis of biomarker candidate genes
Patients from two university hospital centers were included in this retrospective analysis. Inclusion criteria were histologically confirmed primary epithelial ovarian cancer, available fresh frozen tissue and at least 36 months follow-up data. After screening of the accessible information and material, samples were excluded due to underrepresentation of tumor cells (< 10%) and tissue derived from relapse events. A total of 203 tumor samples existed for the validation of biomarker candidate genes. Samples from patients with known or identified BRCA1/2 mutations (7.5%; n = 15) were evaluated in a separate analysis (see below) to exclude any bias caused by the known prognostic value of BRCA mutations . Basic clinical data from the analyzed patients are shown in Table 1.
The cohort from Jena University Hospital (n = 101) consists of 65 new patients but also includes 36 samples from our previous study . To enable a correct validation, all samples used for initial marker identification were excluded. To validate the robustness of employed MSP assays and to exclude major bias by tumor heterogeneity, several experiments were done. First of all, a subset of samples analyzed in Häfner et al.  was re-tested. Since several years and different other analyses were performed using the very same tumor block, a deeper section of the tumor was analyzed. Secondly, from some patients, multiple biopsy tissues were available. Both analyses revealed a high consistency of the methylation pattern throughout the samples (95% concordance, Cohen’s Kappa = 0.881) proving the validity and robustness of methylation marker detection (Additional file 1: Fig. S1).
Altogether the cohort reflects the typical EOC patient and mainly comprises of late stage HGSOC patients who were treated by macroscopically complete resection and platin-taxane-based chemotherapy (Table 1). However, the cohort is more diverse than the establishing cohort in regard to tumor stage, histology, surgical resection status and adjuvant chemotherapy. Specifically, a subset of patients received maintenance treatment with Bevacizumab, currently standard of care for advanced stage EOC patients in Germany . Evaluating the clinical parameters of the 188 included sporadic EOC patients, we confirmed the prognostic value of the FIGO staging (FIGO I + II vs. FIGO III + IV) and surgery outcome (tumor residual 0 mm vs. < 10 mm/ > 10 mm) for overall survival (OS) and progression-free survival (PFS) (Table 2; Additional file 1: Fig. S2; LogRank p < 0.05). High-grade serous EOC tend to show a shorter PFS than other EOC (p = 0.08). The differentiation stage of the tumors (G1/2 vs. G3) did not show a prognostic value for PFS or OS.
Table 2 summarizes, in addition to the clinical parameters, the median and 95% confidence interval of PFS and OS as determined by univariate Kaplan–Meier analyzes for patients stratified by methylation biomarkers or combinations thereof. Combinations reflect the methylation status of at least one included marker (logical <OR> combination). Methylation of CAMK2N1, KRT86 and the combinations CAMK2N1/RUNX3 or CAMK2N1/RUNX3/KRT86 stratified patients into groups with significantly different PFS (Fig. 1, Table 2).
This confirmed our earlier results for these markers, whereas the prognostic value of the singular RUNX3 methylation could not be validated (Fig. 1, Table 2). Albeit KRT86 methylation showed the most significant discriminatory power in univariate analysis, it exhibited no independent prognostic value (p = 0.055, Table 3). This is mainly related to the exclusive methylation in late stage EOC and the resulting significant association of KRT86 methylation with FIGO staging (Fisher’s exact test, p < 0.01). The present large cohort with overall survival data enabled the additional determination of prognostic effects related to overall survival. Only the combination of methylation markers CAMK2N1/RUNX3 exhibited a significant prognostic value for OS (Fig. 1, Table 2). The other methylation marker could not significantly differentiate between patient groups with varying OS (Fig. 1, Table 2). However, patients differentiated by CAMK2N1 methylation as singular marker exhibited a trend to varying overall survival (p = 0.06). Additionally, to the marker validation, we included methylation analyses for BRCA1. We could not identify any prognostic value for this parameter (Table 2) confirming data from the TCGA cohort and a recent meta-analysis [7, 53].
All methylation markers showing a significant prognostic value in univariate analysis (LogRank test) were separately analyzed together with both significant clinical parameters (FIGO and resection status) in multivariate Cox regression analysis (Table 3). We also analyzed clinical parameters with non-significant results in univariate analysis (histology, grading, treatment type) but they did neither exhibit prognostic value in multivariate analysis nor change results for the other clinical or molecular marker (data not shown). The combinations CAMK2N1/RUNX3 and CAMK2N1/RUNX3/KRT86 retained their significant prognostic value for PFS in this analysis. The calculated hazard ratio of progression for late stage EOC, after suboptimal debulking, for CAMK2N1/RUNX3 methylation and for CAMK2N1/RUNX3/KRT86 methylation was 5.029, 1.865, 1.587 and 1.543, respectively (Table 3). Thus, the above combinations of methylation markers had a prognostic value for PFS independent from FIGO staging and resection status. In relation to overall survival, only the resection status and CAMK2N1/RUNX3 methylation showed an independent prognostic value (HR 2.45 and 1.60, respectively; Table 3). Late stage EOC exhibited a hazard ratio of 2.08 with p = 0.051.
Explorative analysis of patients with defects in homologous recombination DNA repair (HRD)
A total of 10 patients with a familiar predisposition for EOC and a known mutation in BRCA1, BRCA2 or RAD51C and available tumor tissue were included in the cohort. The information of the exact mutation was reported by genetic counseling (Additional file 1: Table S1). Additionally, we screened young EOC patients with available tissue for an exploratory somatic BRCA1 mutation analysis and could identify five patients with mutations and included them in this cohort (n = 15, named HRD-defective EOC afterward, Additional file 1: Table S1).
Because of the small number of samples with mutations, all analyses have explorative character. We did not detect any significant differences in clinical variables between BRCAwt and HRD-defective EOC (data not shown). This is likely caused by the fact, that our cohort is enriched for high-grade serous EOC resembling BRCA-mutated EOC. Only KRT86 showed significant different methylation frequencies in BRCAwt and HRD-defective EOC (22.9% and 0%, respectively; p = 0.044 Fisher’s exact Test). A trend of similar differences was seen for CAMK2N1 and BRCA1 (8.5% and 11.2% methylated vs. absent methylation in BRCAwt vs. HRD-defective EOC, respectively). However, no significant difference was detected in the combined methylation frequency for all genes (60% of samples in both groups). Additionally, the prognostic value for PFS did not strongly differ between BRCAwt and HRD-defective EOC for most methylation marker although the small sample size precluded significant results (data not shown). Interestingly, we observed a difference in prognosis of patients with KATNAL2-methylated BRCAwt and HRD-deficient tumors compared to the respective unmethylated group. Because HRD-deficiency can also be potentially caused by BRCA1 methylation, we analyzed the prognostic value of KATNAL2 methylation in BRCAwt EOC with and without BRCA1 methylation and HRD-deficient EOC (Fig. 2). Whereas KATNAL2 methylation showed no prognostic effect in the complete cohort (n = 188, Table 2) the stratification into BRCA aberration dependent groups resulted in significant prognostic effects (p < 0.05, LogRank test). KATNAL2 methylation correlated with a significant worse outcome in patients with BRCAwt, but with a prolonged PFS in HRD-deficient and BRCA1-methylated EOC patients compared to the patients with respective KATNAL2 unmethylated tumors (29 vs. 37 month (p = 0.012, HR 1.591), 35 vs. 31 month (n.s., HR 0.620) and 32 vs. 19 month (p = 0.04, HR 0.332), respectively; Fig. 2).
Aberrant DNA methylation can be used as molecular biomarker for disease diagnosis, prognosis or therapy stratification [58, 59]. Importantly, both the specific use and the target population must be defined and validated. Biomarker establishment and validation must balance between the power of the marker, the quality of the results and the homogeneity of the target population. In previous work, we identified 5 methylation markers that were prognostic for PFS by univariate analysis in a defined population—late stage type II EOC treated with platinum-taxane-based chemotherapy . The aim of the present study was to analyze the potential independent prognostic value of these markers for PFS and OS in a group of EOC patients more heterogeneous in regard to tumor stage, grading, tumor histology and treatment by multivariate analyses. To provide excellent sample quality and correct sample classification, pathological evaluation of tissues was done. Because samples originated from two different medical centers, handling differences cannot be excluded, but the processing (cryo-sectioning, DNA purification), bisulfite conversion and MSP was done in one laboratory only. Moreover, the robustness of methylation analyses was determined by repeated assays from selected patients using different tissue blocks of the same tumor or repeated sections from identical biopsies resulting in an almost perfect agreement (Cohen’s kappa 0.881). We included samples with tumor cell fraction > 10% and did not observe a strong influence of the tumor content on methylation frequency of candidate regions (Additional file 1: Fig. S3). This may relate to the low cut-off value used to assign a positive methylation status (1% relative methylation, see method section) and absent methylation in normal keratinocytes, PBMCs and whole tissue sections of normal ovaries (data not shown). However, several limitations prohibit the transfer of the presented results to the general EOC population. First of all, the analyzed cohort does not represent all histological subtypes. Secondly, both median progression-free and overall survival are longer for our cohort compared to the general EOC population in Germany (31 month vs. 18.2 month PFS, 73 month vs. 44.1 month OS; ). The 5-year survival rate of 56% in our cohort is higher than the population-based rate in Germany (41%) . Reasons for this are (1) the selection of platin-sensitive EOC and (2) the high proportion of optimal resected patients (macroscopically tumor-free, ~ 80%). A third limitation is that only a subset of patients had available information about germline mutation status (n = 10). An additional exploratory somatic BRCA1 mutational analysis was done in early onset patients (age at diagnosis < 40 years) and revealed 5 patients with BRCA1 mutations. Recent studies in Germany detected an approximately 20% BRCA germline mutation rate in unselected EOC patients (Harter ) suggesting that we underestimate the HRD-deficient fraction in our cohort (15/203, 7.4%). Nevertheless, this may only affect the prognostic value of KATNAL2 (see below) because other markers seem to be not affected by HRD-status.
Using this enlarged, more heterogeneous cohort of EOC patients (Table 1) we could confirm the prognostic value of singular markers (CAMK2N1, KRT86) and combinations (CAMK2N1/RUNX3; CAMK2N1/RUNX3/KRT86) for PFS in univariate analysis. Additionally, both combinations are prognostic independent from tumor stage and resection status (multivariate Cox regression analysis). The stronger effect of combinations is not caused by the pure combination of markers potentially resulting in the detection of EOC with generally higher methylation level because EOC unmethylated for all analyzed marker showed no significant different survival compared to the group of EOC with at least one methylated marker (Table 2). It is more likely that informative methylation markers detect specific but different subgroups of EOC with reduced survival. This is supported by the low methylation frequency of informative markers (CAMK2N1 7.9%, KRT86 22.3%, RUNX3 14.9%). Whereas CAMK2N1 and RUNX3 methylation is independent from FIGO staging, KRT86 methylation is not an independent prognostic factor but associated with late stage EOC. Thus, KRT86 methylation can potentially be used for pre-operative staging if it can be specifically detected in systemic samples, e.g., blood. Besides the prognostic effect for PFS we could also detect an association of the marker combination CAMK2N1/RUNX3 with overall survival both in univariate and multivariate analysis. EOC patients with methylation of CAMK2N1 or RUNX3 have a shorter overall survival than patients without methylation of these genes independently from FIGO staging or tumor resection status. The effect on overall survival is mainly caused by the shorter PFS in CAMK2N1/RUNX3 methylation positive EOC patients (Fig. 1, Table 2) as Kaplan–Meier analyses for the time interval between relapse and death or last information available did not show a prognostic value for the marker combination (subset of relapsed patients n = 133 (70.7%); data not shown).
If prospectively validated CAMK2N1/RUNX3 methylation may identify a subgroup of platin-sensitive patients that have worse prognosis. It can be speculated that this subgroup may benefit from maintenance treatment with Bevacizumab because both GOG-218 and ICON7 identified an overall survival benefit in a group of patients with high risk for progression differently defined by clinical variables [61, 62]. Although 29 patients of our cohort were treated with Bevacizumab, the number of methylation events (n = 2) is too low to draw any conclusion. However, published and unpublished data point to an influence of RUNX3 on the angiogenic potential of EOC by influencing the expression of thrombospondin-1 and other angiogenic factors ([63, 64], Heinze et al. in preparation). We also tried to identify an association of methylation marker status with BRCA status (BRCA mutation, BRCA1 methylation) but did not detect a largely different prognostic value of CAMK2N1/RUNX3 methylation in BRCAwt and potentially HRD-deficient EOC. Nevertheless, KATNAL2 methylation being not prognostic in the complete group of EOC showed opposite effects in BRCAwt EOC without BRCA1 methylation (HR 1.591, p = 0.012) and BRCAwt EOC with BRCA1 methylation (HR 0.332, p = 0.04) or BRCAmut EOC (HR 0.620, n.s.). Although the numbers of patients with either BRCA aberration are small, this suggests different consequences of KATNAL2 methylation in BRCAwt versus HRD-deficient EOC. Shared clinico-pathological features between BRCA1-mutated and BRCA1-methylated tumors may relate to the inactivation of BRCA1 albeit only BRCA1 mutations significantly influence patients survival . Likely, we did not identify a prognostic value of BRCA1 methylation (Table 2). For patients with BRCA1-mutated tumors, survival analyses are limited by the low number and the shorter follow-up period for these patients (median mutated vs. wt: 36 month vs. 55 month). The PFS was not significantly different in BRCA1-mutated and BRCA1wt patients in our cohort (data not shown). In preliminary experiments, we analyzed BRCA1 protein expression by immunohistochemistry in a subset of samples and could observe significantly different staining patterns between BRCA1wt and mutated or methylated tumors (Additional file 1: Table S3) supporting earlier results from methylated tumors . Both groups of tumors with BRCA1 aberration express lower levels of BRCA1. Thus, loss of BRCA1 may indeed result in the switch of the prognostic value of KATNAL2 methylation. KATNAL2 (Katanin Catalytic Subunit A1 Like 2) is a gene coding for a protein with microtubule binding and microtubule-severing ATPase activity whose detailed function in carcinogenesis is largely unknown. However, methylation of KATNAL2 is a potential marker for severe cervical intraepithelial neoplasia . The Katanin gene family regulates microtubule dynamics during mitosis, migration, ciliogenesis or cellular reorganization . Strong KATNAL2 downregulation leads to massive cell death, whereas mild downregulation causes microtubule stabilization and reduced turnover leading to mitotic defects and G2/M arrest in mouse cells . Therefore, it can be speculated that microtubule stabilization by methylation-dependent KATNAL2 downregulation induces mitotic defects, chromosome missegregation and chromosomal instability (CIN) potentially enabling tumor adaptation to chemotherapeutic intervention and worse outcome in BRCA1wt patients. However, BRCA proteins are not only involved in DNA repair but also microtubule organization—specifically in mitosis and cell polarization [68, 69]. Importantly, BRCA1 loss results in disturbed mitosis and CIN . Thus, in tumors with BRCA1 aberrations, the sum of mitotic defects may result in mitotic catastrophe, increased tumor cell death/better therapeutic response and improved survival. Besides the unknown basic mechanism for the different prognostic impact of KATNAL2 methylation, depending on HRD status, this points to the potential necessity to differentiate between these EOC groups for biomarker identification and validation studies.
New biomarker candidates for cancer are identified frequently, a small subset of them will be published in peer-reviewed journals and only a few are validated in independent studies. However, even a positive validation does not guarantee a successful biomarker candidate . Only ~ 0.8% of published methylation markers are commercially available . Additional knowledge of the underlying mechanism of the biomarker’s detection and function would be helpful to underline its significance. One validation step for DNA methylation marker can consist of immunohistochemical detection of the epigenetically regulated protein. However, DNA methylation may not result in largely different expression a priori but rather restricts the capacity for an increased expression. Such regulation of the expression potential of affected genes was observed for genes hypermethylated during platin resistance development in EOC cell lines (i.e., CAMK2N1 and TRIB2) in our laboratory [44, 73]. Both genes show gene expression differences between isogenic pairs of sensitive, unmethylated and resistant, hypermethylated cells under platin treatment only. Using the presently available chemotherapy-naïve EOC samples, we may not detect such changes in gene expression. Moreover, RUNX3 methylation affects only promoter P2 and the associated isoform that cannot be differentiated from other RUNX3 proteins by immunohistochemistry so far. Hence, we have not analyzed the protein expression of our candidate genes. However, first functional analyses of the biomarker candidate genes were conducted and a part of it is already published . Specifically, the genes of the identified best methylation marker combination CAMK2N1 and RUNX3 were shown to have tumor suppressive functions (CAMK2N1) or to be involved in platin resistance and migration (RUNX3). CAMK2N1, one of two endogenous CAMKII inhibitors, was firstly identified as tumor suppressor inducing cell cycle arrest by p27 stabilization . Accordingly, CAMK2N1 was not only described as tumor suppressor in EOC by our group but in multiple myeloma, oral squamous cell carcinoma, prostate cancer and thyroid cancer throughout the last years [75,76,77,78]. Therefore, a reduced PFS and OS of patients with CAMK2N1-methylated EOC vs. unmethylated EOC (p = 0.021 and p = 0.06, LogRank test) is in agreement with the gene function. RUNX3 is described both as oncogene and tumor suppressor gene depending on the analyzed tumor entity or cellular background and readout [79, 80]. A tumor suppressive role of RUNX3 in vitro is described for Wilm’s tumor, prostate cancer, hepatocellular carcinoma, lung cancer and glioma [81,82,83,84,85]. In addition, we have shown dichotomous functions of the RUNX3 isoforms . Importantly, methylation of RUNX3 affects only promotor P2 and the associated transcript variant 2. This transcript variant increases cisplatin sensitivity and reduces cell migration for A2780 and cisplatin sensitivity of OVCAR3 in vitro ( and unpublished data). In the present study, we could not confirm a significant prognostic value of singular RUNX3 methylation but a trend for decreased PFS for patients with tumors showing methylated RUNX3 (median PFS 21 month vs. 31 month; p = 0.1 LogRank test).
In conclusion, the retrospective analysis of 188 sporadic platinum-sensitive EOC proved an independent prognostic value of the methylation marker combination CAMK2N1/RUNX3 for PFS and OS. If validated prospectively this combination may identify EOC patients with worse prognosis after standard therapy potentially benefiting from intensive follow-up, maintenance therapies or inclusion in therapeutic studies.
A total of 203 fresh frozen (FF) and formalin-fixed paraffin-embedded (FFPE) tissue samples were used and divided into a validation and BRCA mutation cohort. The samples originated from the University Medical Center Hamburg-Eppendorf and the Jena University Hospital from 1996 to 2016. The use of patient’s material was approved by the Ethics committee Jena (#2582-6/09) and Hamburg (#190505). Clinical data were retrieved from the medical records. Progression-free survival was calculated from the date of primary surgery to first occurrence of relapse (confirmed by second-look surgery or non-invasive diagnostic tools) or last follow-up and overall survival was calculated from date of primary surgery until death or last follow-up. Surgical samples not needed for diagnostic procedures were snap-frozen in liquid nitrogen and stored at − 80 °C. The samples for the validation set were included based on the existence of three-years-follow-up data. Detailed and summary information on those two cohorts are depicted in Table 1 and Additional file 1: Table S1. Hematoxylin–Eosin stained sections of analyzed tissues were evaluated for histopathological classification and estimation of tumor cell fraction. Only tissues with > 10% tumor cells were used and 24.6%, 44.2% and 22.8% of samples contained 10–40%, 40–80% and > 80% tumor cells.
Sections from FF and FFPE tissue were used to perform analysis. The isolation was done according to the phenol-chloroform-isoamyl alcohol (PCI)-based DNA isolation protocol and using the QIAamp Kit according to the manufacture instruction.
RNA isolation and cDNA synthesis for BRCA mutation detection
Sections from selected FF tissue were used for RNA isolation and subsequent cDNA synthesis. The isolation was done with the NucleoSpin® RNA kit according to the manual (Macherey–Nagel). 500 ng of RNA was reverse transcribed according to established techniques following published protocols .
500 ng–1 µg of genomic DNA derived from tissue was converted using the EZ DNA Methylation-Gold™ according to the manual instruction (Zymo Research). The bsDNA concentration was determined and adjusted to 5 ng/µl. Identically treated gDNA from cervical swab of a healthy individual served as negative control sample. Positive control sample was created by in vitro methylation using the CpG-Methyltransferase (M.SssI), according the supplier’s instruction (New England BioLabs). A serial dilution using in vitro methylated and unmethylated control bsDNA was created to generate the standard of 0%, 1%, 10% and 100% methylation level.
The q-MSP was performed according to  using the RotorGeneQ thermocycler (Qiagen, Germany) and candidate gene-specific primer listed in Additional file 1: Table S2 [44, 52]. All MSP assays detect DNA methylation within CpG islands located in the promoter/first exon region of associated genes except KRT86. The analyzed CpG island for KRT86 is located at the 3’end (last exon) of the gene. Melting curve analyses confirmed the specificity of MSP amplification. The level of methylation was quantified stepwise by (1) calculating the relative amount (2∆Ct, relative methylation) of methylated target sequences in relation to a beta-actin fragment which was amplified by bisulfite-DNA-specific, but methylation-independent PCR (BS-PCR) and (2) comparing the relative methylation to an artificial dilution series (0–100% methylated DNA). Methylation data were evaluated qualitatively by application of a cutoff value and only samples showing a higher relative methylation as the 1% control were scored as methylated. Marker combinations reflect an <OR> combination of single marker results.
BRCA1 mutation detection
To detect BRCA1 exon mutations, a total of 6 PCRs were done using cDNA derived from patient’s RNA. Approximately 20 ng reverse transcribed nucleic acid were used. The reactions were conducted in 50 µl volume containing: dNTPs (240 µM each), forward and reverse primer (10 pmol each), DMSO (5%), MgCl2 (1.75 mM), AmpliTaqGold (1.25 U) and respective reaction buffer II (Applied Biosystems). The PCR steps were as follows: initial denaturation and activation at 95 °C for 10 min followed by 25 cycles of denaturation phase at 95 °C for 15 s, primer-specific annealing for 20 s at 56 °C and elongation at 72 °C for 30 s. The reaction products were size separated and sequenced in both directions. The readings were aligned to the BRCA1 sequence according to human reference sequence (hg19).
The statistical analysis was executed using SPSS and Microsoft Excel 2010/365. The unpaired student’s t-test was used for pairwise comparisons of continuous variables whereas chi-squared or Fishers Exact test was used for ordinal variables. The difference in survival between groups of patients was depicted by Kaplan–Meier plots and evaluated by LogRank test. Multivariate Cox regression analyses were used to calculated hazard ratios and determine the independent prognostic value of methylation markers.
Availability of data and materials
All data generated or analyzed during this study are included in this published article [and its supplementary information files]. Raw data are available from the corresponding author on reasonable request.
Epithelial ovarian cancer
International Federation of Gynecology and Obstetrics
Hereditary breast and ovarian cancer syndrome
High-grade serous ovarian cancer
Homologous DNA repair defects
Methyl-binding domain protein
RKI, GEKID. Krebs in Deutschland 2013/2014. Berlin: Robert-Koch-Institut; 2017. p. 88–91.
Lynch HT, Casey MJ, Snyder CL, et al. Hereditary ovarian carcinoma: heterogeneity, molecular genetics, pathology, and management. Mol Oncol. 2009;3(2):97–137.
Kobayashi H, Ohno S, Sasaki Y, et al. Hereditary breast and ovarian cancer susceptibility genes (review). Oncol Rep. 2013;30(3):1019–29.
Nielsen FC, van Overeem HT, Sorensen CS. Hereditary breast and ovarian cancer: new genes in confined pathways. Nat Rev Cancer. 2016;16(9):599–612.
Norquist BM, Harrell MI, Brady MF, et al. Inherited mutations in women with ovarian carcinoma. JAMA Oncol. 2016;2(4):482–90.
Harter P, Hauke J, Heitz F, et al. Prevalence of deleterious germline variants in risk genes including BRCA1/2 in consecutive ovarian cancer patients (AGO-TR-1). PLoS ONE. 2017;12(10):e0186043.
Cancer Genome Atlas Research Network. Integrated genomic analyses of ovarian carcinoma. Nature. 2011;474(7353):609–15.
Hauke J, Hahnen E, Schneider S, et al. Deleterious somatic variants in 473 consecutive individuals with ovarian cancer: results of the observational AGO-TR1 study (NCT02222883). J Med Genet. 2019;56(9):574–80.
Vaughan S, Coward JI, Bast RC Jr, et al. Rethinking ovarian cancer: recommendations for improving outcomes. Nat Rev Cancer. 2011;11(10):719–25.
Rescigno P, Cerillo I, Ruocco R, et al. New hypothesis on pathogenesis of ovarian cancer lead to future tailored approaches. Biomed Res Int. 2013;2013:852839.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2017. CA Cancer J Clin. 2017;67(1):7–30.
Sorensen SM, Schnack TH, Hogdall C. Impact of residual disease on overall survival in women with Federation of Gynecology and Obstetrics stage IIIB-IIIC vs stage IV epithelial ovarian cancer after primary surgery. Acta Obstet Gynecol Scand. 2019;98(1):34–43.
du Bois A, Reuss A, Pujade-Lauraine E, et al. Role of surgical outcome as prognostic factor in advanced epithelial ovarian cancer: a combined exploratory analysis of 3 prospectively randomized phase 3 multicenter trials: by the Arbeitsgemeinschaft Gynaekologische Onkologie Studiengruppe Ovarialkarzinom (AGO-OVAR) and the Groupe d’Investigateurs Nationaux Pour les Etudes des Cancers de l’Ovaire (GINECO). Cancer. 2009;115(6):1234–44.
Greenlee RT, Hill-Harmon MB, Murray T, et al. Cancer statistics, 2001. CA Cancer J Clin. 2001;51(1):15–36.
Davis A, Tinker AV, Friedlander M. “Platinum resistant” ovarian cancer: what is it, who to treat and how to measure benefit? Gynecol Oncol. 2014;133(3):624–31.
Bolton KL, Chenevix-Trench G, Goh C, et al. Association between BRCA1 and BRCA2 mutations and survival in women with invasive epithelial ovarian cancer. JAMA. 2012;307(4):382–90.
McLaughlin JR, Rosen B, Moody J, et al. Long-term ovarian cancer survival associated with mutation in BRCA1 or BRCA2. J Natl Cancer Inst. 2013;105(2):141–8.
Ledermann J, Harter P, Gourley C, et al. Olaparib maintenance therapy in patients with platinum-sensitive relapsed serous ovarian cancer: a preplanned retrospective analysis of outcomes by BRCA status in a randomised phase 2 trial. Lancet Oncol. 2014;15(8):852–61.
Mirza MR, Monk BJ, Herrstedt J, et al. Niraparib maintenance therapy in platinum-sensitive, recurrent ovarian cancer. N Engl J Med. 2016;375(22):2154–64.
Hill SJ, Decker B, Roberts EA, et al. Prediction of DNA repair inhibitor response in short-term patient-derived ovarian cancer organoids. Cancer Discov. 2018;8(11):1404–21.
Gentles L, Goranov B, Matheson E, et al. Exploring the frequency of homologous recombination DNA repair dysfunction in multiple cancer types. Cancers (Basel). 2019;11(3):354.
Gonzalez-Martin A, Pothuri B, Vergote I, et al. Niraparib in patients with newly diagnosed advanced ovarian cancer. N Engl J Med. 2019;381(25):2391–402.
Matulonis UA, Sood AK, Fallowfield L, et al. Ovarian cancer Nat Rev Dis Primers. 2016;2:16061.
Asante DB, Calapre L, Ziman M, et al. Liquid biopsy in ovarian cancer using circulating tumor DNA and cells: ready for prime time? Cancer Lett. 2020;468:59–71.
Mari R, Mamessier E, Lambaudie E, et al. Liquid biopsies for ovarian carcinoma: how blood tests may improve the clinical management of a deadly disease. Cancers (Basel). 2019;11(6):774.
Sheng R, Li X, Wang Z, et al. Circular RNAs and their emerging roles as diagnostic and prognostic biomarkers in ovarian cancer. Cancer Lett. 2020;473:139–47.
Zhu Y, Zhou S, Liu Y, et al. Prognostic value of systemic inflammatory markers in ovarian cancer: a PRISMA-compliant meta-analysis and systematic review. BMC Cancer. 2018;18(1):443.
Frawley T, Piskareva O. Extracellular vesicle dissemination of epidermal growth factor receptor and ligands and its role in cancer progression. Cancers (Basel). 2020;12(11):3200.
Drakes ML, Stiff PJ. Regulation of ovarian cancer prognosis by immune cells in the tumor microenvironment. Cancers (Basel). 2018;10(9):302.
Etemadmoghadam D, deFazio A, Beroukhim R, et al. Integrated genome-wide DNA copy number and expression analysis identifies distinct mechanisms of primary chemoresistance in ovarian carcinomas. Clin Cancer Res. 2009;15(4):1417–27.
Etemadmoghadam D, Weir BA, Au-Yeung G, et al. Synthetic lethality between CCNE1 amplification and loss of BRCA1. Proc Natl Acad Sci USA. 2013;110(48):19489–94.
Konecny GE, Winterhoff B, Wang C. Gene-expression signatures in ovarian cancer: promise and challenges for patient stratification. Gynecol Oncol. 2016;141(2):379–85.
Samuel P, Carter DR. The diagnostic and prognostic potential of microRNAs in epithelial ovarian carcinoma. Mol Diagn Ther. 2017;21(1):59–73.
Millstein J, Budden T, Goode EL, et al. Prognostic gene expression signature for high-grade serous ovarian cancer. Ann Oncol. 2020;31(9):1240–50.
Taby R, Issa JP. Cancer epigenetics. CA Cancer J Clin. 2010;60(6):376–92.
Jones PA, Baylin SB. The epigenomics of cancer. Cell. 2007;128(4):683–92.
Robertson KD. DNA methylation and human disease. Nat Rev Genet. 2005;6(8):597–610.
Jin Z, Liu Y. DNA methylation in human diseases. Genes Dis. 2018;5(1):1–8.
Duffy MJ, Napieralski R, Martens JW, et al. Methylated genes as new cancer biomarkers. Eur J Cancer. 2009;45(3):335–46.
Barton CA, Hacker NF, Clark SJ, et al. DNA methylation changes in ovarian cancer: implications for early diagnosis, prognosis and treatment. Gynecol Oncol. 2008;109(1):129–39.
Gifford G, Paul J, Vasey PA, et al. The acquisition of hMLH1 methylation in plasma DNA after chemotherapy predicts poor survival for ovarian cancer patients. Clin Cancer Res. 2004;10(13):4420–6.
Michaelson-Cohen R, Keshet I, Straussman R, et al. Genome-wide de novo methylation in epithelial ovarian cancer. Int J Gynecol Cancer. 2011;21(2):269–79.
Patch AM, Christie EL, Etemadmoghadam D, et al. Whole-genome characterization of chemoresistant ovarian cancer. Nature. 2015;521(7553):489–94.
Hafner N, Steinbach D, Jansen L, et al. RUNX3 and CAMK2N1 hypermethylation as prognostic marker for epithelial ovarian cancer. Int J Cancer. 2016;138(1):217–28.
Ibanez de Caceres I, Battagli C, Esteller M, et al. Tumor cell-specific BRCA1 and RASSF1A hypermethylation in serum, plasma, and peritoneal fluid from ovarian cancer patients. Cancer Res. 2004;64(18):6476–81.
Strathdee G, MacKean MJ, Illand M, et al. A role for methylation of the hMLH1 promoter in loss of hMLH1 expression and drug resistance in ovarian cancer. Oncogene. 1999;18(14):2335–41.
Arnold JM, Cummings M, Purdie D, et al. Reduced expression of intercellular adhesion molecule-1 in ovarian adenocarcinomas. Br J Cancer. 2001;85(9):1351–8.
Yuecheng Y, Hongmei L, Xiaoyan X. Clinical evaluation of E-cadherin expression and its regulation mechanism in epithelial ovarian cancer. Clin Exp Metastasis. 2006;23(1):65–74.
Potapova A, Hoffman AM, Godwin AK, et al. Promoter hypermethylation of the PALB2 susceptibility gene in inherited and sporadic breast and ovarian cancer. Cancer Res. 2008;68(4):998–1002.
Ruscito I, Dimitrova D, Vasconcelos I, et al. BRCA1 gene promoter methylation status in high-grade serous ovarian cancer patients—a study of the tumour Bank ovarian cancer (TOC) and ovarian cancer diagnosis consortium (OVCAD). Eur J Cancer. 2014;50(12):2090–8.
Baldwin RL, Nemeth E, Tran H, et al. BRCA1 promoter region hypermethylation in ovarian carcinoma: a population-based study. Cancer Res. 2000;60(19):5329–33.
Radosa MP, Hafner N, Camara O, et al. Loss of BRCA1 protein expression as indicator of the BRCAness phenotype is associated with favorable overall survival after complete resection of sporadic ovarian cancer. Int J Gynecol Cancer. 2011;21(8):1399–406.
Kalachand RD, Stordal B, Madden S, et al. BRCA1 promoter methylation and clinical outcomes in ovarian cancer: an individual patient data meta-analysis. J Natl Cancer Inst. 2020. https://doi.org/10.1093/jnci/djaa070.
Prieske K, Prieske S, Joosse SA, et al. Loss of BRCA1 promotor hypermethylation in recurrent high-grade ovarian cancer. Oncotarget. 2017;8(47):83063–74.
Natanzon Y, Goode EL, Cunningham JM. Epigenetics in ovarian cancer. Semin Cancer Biol. 2018;51:160–9.
Heinze K, Kritsch D, Mosig AS, et al. Functional analyses of RUNX3 and CaMKIINalpha in ovarian cancer cell lines reveal tumor-suppressive functions for CaMKIINalpha and dichotomous roles for RUNX3 transcript variants. Int J Mol Sci. 2018;19(1):253.
Leitlinienprogramm Onkologie: S3-Leitlinie Diagnostik, Therapie und Nachsorge maligner Ovarialtumoren AWMF-Registernummer: 032/035OL (2019).
Beltran-Garcia J, Osca-Verdegal R, Mena-Molla S, et al. Epigenetic IVD tests for personalized precision medicine in cancer. Front Genet. 2019;10:621.
Thomas ML, Marcato P. Epigenetic modifications as biomarkers of tumor development, therapy response, and recurrence across the cancer care continuum. Cancers (Basel). 2018;10(4):101.
Hanker LC, Loibl S, Burchardi N, et al. The impact of second to sixth line therapy on survival of relapsed ovarian cancer after primary taxane/platinum-based therapy. Ann Oncol. 2012;23(10):2605–12.
Oza AM, Cook AD, Pfisterer J, et al. Standard chemotherapy with or without bevacizumab for women with newly diagnosed ovarian cancer (ICON7): overall survival results of a phase 3 randomised trial. Lancet Oncol. 2015;16(8):928–36.
Tewari KS, Burger RA, Enserro D, et al. Final overall survival of a randomized trial of bevacizumab for primary treatment of ovarian cancer. J Clin Oncol. 2019;37(26):2317–28.
Shi X, Deepak V, Wang L, et al. Thrombospondin-1 is a putative target gene of Runx2 and Runx3. Int J Mol Sci. 2013;14(7):14321–32.
Peng Z, Wei D, Wang L, et al. RUNX3 inhibits the expression of vascular endothelial growth factor and reduces the angiogenesis, growth, and metastasis of human gastric cancer. Clin Cancer Res. 2006;12(21):6386–94.
Lendvai A, Johannes F, Grimm C, et al. Genome-wide methylation profiling identifies hypermethylated biomarkers in high-grade cervical intraepithelial neoplasia. Epigenetics Off J DNA Methylation Soc. 2012;7(11):1268–78.
Kuo YW, Howard J. Cutting, amplifying, and aligning microtubules with severing enzymes. Trends Cell Biol. 2020;31:50–61.
Ververis A, Christodoulou A, Christoforou M, et al. A novel family of katanin-like 2 protein isoforms (KATNAL2), interacting with nucleotide-binding proteins Nubp1 and Nubp2, are key regulators of different MT-based processes in mammalian cells. Cell Mol Life Sci. 2016;73(1):163–84.
Maxwell CA, Benitez J, Gomez-Baldo L, et al. Interplay between BRCA1 and RHAMM regulates epithelial apicobasal polarization and may influence risk of breast cancer. PLoS Biol. 2011;9(11):e1001199.
Petsalaki E, Zachos G. DNA damage response proteins regulating mitotic cell division: double agents preserving genome stability. FEBS J. 2020;287(9):1700–21.
Stolz A, Ertych N, Kienitz A, et al. The CHK2-BRCA1 tumour suppressor pathway ensures chromosomal stability in human somatic cells. Nat Cell Biol. 2010;12(5):492–9.
Hayes DF. Biomarker validation and testing. Mol Oncol. 2015;9(5):960–6.
Koch A, Joosten SC, Feng Z, et al. Analysis of DNA methylation in cancer: location revisited. Nat Rev Clin Oncol. 2018;15(7):459–66.
Kritsch D, Hoffmann F, Steinbach D, et al. Tribbles 2 mediates cisplatin sensitivity and DNA damage response in epithelial ovarian cancer. Int J Cancer. 2017;141(8):1600–14.
Wang C, Li N, Liu X, et al. A novel endogenous human CaMKII inhibitory protein suppresses tumor growth by inducing cell cycle arrest via p27 stabilization. J Biol Chem. 2008;283(17):11565–74.
Russo E, Salzano M, De Falco V, et al. Calcium/Calmodulin-dependent protein kinase II and its endogenous inhibitor alpha in medullary thyroid cancer. Clin Cancer Res. 2014;20(6):1513–20.
Wang T, Liu Z, Guo S, et al. The tumor suppressive role of CAMK2N1 in castration-resistant prostate cancer. Oncotarget. 2014;5(11):3611–21.
Xu K, Hu X, Sun L, et al. MicroRNA-532 exerts oncogenic functions in t(4;14) multiple myeloma by targeting CAMK2N1. Hum Cell. 2019;32(4):529–39.
Li N, Nan CC, Zhong XY, et al. miR-182-5p promotes growth in oral squamous cell carcinoma by inhibiting CAMK2N1. Cell Physiol Biochem. 2018;49(4):1329–41.
Ito Y, Bae SC, Chuang LS. The RUNX family: developmental regulators in cancer. Nat Rev Cancer. 2015;15(2):81–95.
Whittle MC, Izeradjene K, Rani PG, et al. RUNX3 controls a metastatic switch in pancreatic ductal adenocarcinoma. Cell. 2015;161(6):1345–60.
Li X, Zhong M, Wang J, et al. miR-301a promotes lung tumorigenesis by suppressing Runx3. Mol Cancer. 2019;18(1):99.
Steponaitis G, Kazlauskas A, Vaitkiene P, et al. Oncosuppressive role of RUNX3 in human astrocytomas. J Oncol. 2019;2019:1232434.
Gu H, Gu S, Zhang X, et al. miR-106b-5p promotes aggressive progression of hepatocellular carcinoma via targeting RUNX3. Cancer Med. 2019;8(15):6756–67.
Fan L, Wang Y, Huo W, et al. MicroRNA-301a-3p overexpression promotes cell invasion and proliferation by targeting runt-related transcription factor 3 in prostate cancer. Mol Med Rep. 2019;20(4):3755–63.
Zhao HZ, Zhao HL, Zhang YN, et al. MicroRNA-199b promotes cell proliferation and invasion in Wilms’ tumour by directly targeting Runt-related transcription factor 3. Mol Med Rep. 2018;18(2):1812–9.
This study was conducted as part of a doctoral thesis project and is published in a modified version at the medical faculty of Friedrich-Schiller-University Jena. This work was supported by a grant of the Deutsche Forschungsgemeinschaft to Norman Häfner.
Open Access funding enabled and organized by Projekt DEAL. This research was funded by Deutsche Forschungsgemeinschaft DFG, Grant Number HA5068/2-3. The funding body had no role in the design of the study or collection, analysis, and interpretation of data or in writing the manuscript.
Ethics approval and consent to participate
The use of patient’s material and clinical data was approved by the Ethics committee Jena (#2582-6/09) and Hamburg (#190505).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Graphical overview of methylations-specific PCR validation results. Figure S2. Kaplan–Meier plots for PFS and OS of EOC patients stratified by clinical variables. Figure S3. Comparison of methylation frequency in tumors with different tumor cell fraction. Table S1. Clinical data and mutation data of EOC with known mutations (n = 15). Table S2. Primer data for MSP and BRCA1 mutation analysis. Table S3. BRCA1 immunohistochemistry results for a subset of tumors depending on BRCA1 aberrations.
About this article
Cite this article
Heinze, K., Rengsberger, M., Gajda, M. et al. CAMK2N1/RUNX3 methylation is an independent prognostic biomarker for progression-free and overall survival of platinum-sensitive epithelial ovarian cancer patients. Clin Epigenet 13, 15 (2021). https://doi.org/10.1186/s13148-021-01006-8
- Epithelial ovarian cancer
- Prognostic marker
- DNA methylation
- CAMK2N1 (CaMKIINα)