CAMK2N1/RUNX3 methylation is an independent prognostic biomarker for progression-free and overall survival of platinum-sensitive epithelial ovarian cancer patients

Background 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. Results 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.). 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. The dichotomous prognostic value of KATNAL2 should be validated in larger sample sets of EOC.

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 [2]. 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 [23]. Biomarker may be analyzed pre-therapeutically (blood-based liquid biopsy) or at least before start of adjuvant chemotherapy (tissuebased 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 [40]. Currently, aberrant DNA methylation in EOC is detected using different approaches such as target-specific methylationspecific PCR (MSP) [41] or whole genome approaches [42][43][44]. A large number of genes are affected by DNA methylation in EOC [40]. For instance, tumor suppressor genes such as BRCA1, RASSF1A [45] or MLH1 [41,46]; cell adhesion genes as ICAM-1 [47] and CDH1 [48] and DNA repair genes PALB2 [49] 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 [54]. MLH1 methylation is potentially linked to chemotherapy response or resistance development [46]. However, lack of candidate biomarker validation and low sensitivity or specificity prevents the clinical use of methylation-based biomarker [55].
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 [44]. 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 [56]. 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 [16]. 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 [44]. 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. [44] 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 [57]. 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 BRCA wt 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 BRCA wt 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 BRCA wt 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 BRCA wt 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 BRCA wt 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  Fig. 2).

Discussion
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 [44]. 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  [60]). The 5-year survival rate of 56% in our cohort is higher than the population-based rate in Germany (41%) [1]. 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 [6]) 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 BRCA wt and potentially HRD-deficient EOC. Nevertheless, KATNAL2 methylation being not prognostic in the complete group of EOC showed opposite effects in BRCA wt EOC without BRCA1 methylation (HR 1.591, p = 0.012) and BRCA wt EOC with BRCA1 methylation (HR 0.332, p = 0.04) or BRCA mut 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 BRCA wt 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 [53]. 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 [52]. 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 [65]. The Katanin gene family regulates microtubule dynamics during mitosis, migration, ciliogenesis or cellular reorganization [66]. 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 [67]. 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 [70]. 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 [71]. Only ~ 0.8% of published methylation markers are commercially available [72]. 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 [56]. 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 [74]. 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 [56]. 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 ( [56] 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).

Conclusion
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.

Patients material
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.

DNA isolation
Sections from FF and FFPE tissue were used to perform analysis. The isolation was done according to the phenolchloroform-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 [56].

Bisulfite conversion
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.

Methylation-specific PCR
The q-MSP was performed according to [44] 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 methylationindependent 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%), MgCl 2 (1.75 mM), Ampli-TaqGold (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).

Statistical analysis
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.
Additional file 1: Figure S1. 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.