Next Article in Journal
Impact of Goal Directed Therapy in Head and Neck Oncological Surgery with Microsurgical Reconstruction: Free Flap Viability and Complications
Next Article in Special Issue
Structural Bases for the Synergistic Inhibition of Human Thymidylate Synthase and Ovarian Cancer Cell Growth by Drug Combinations
Previous Article in Journal
Apoptosis-Inducing TNF Superfamily Ligands for Cancer Therapy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Prognostic Evidence of the miRNA-Based Ovarian Cancer Signature MiROvaR in Independent Datasets

1
Integrated Biology Platform, Department of Applied Research and Technology Development, Fondazione IRCCS Istituto Nazionale dei Tumori, 20133 Milan, Italy
2
Molecular Therapies Unit, Department of Research, Fondazione IRCCS Istituto Nazionale dei Tumori, 20133 Milan, Italy
3
Medical Statistics Unit, University of Campania “Luigi Vanvitelli”, 80138 Naples, Italy
4
Urogynaecological Medical Oncology Unit, Istituto Nazionale Tumori–IRCCS-“Fondazione G. Pascale”, 80131 Naples, Italy
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this paper.
Cancers 2021, 13(7), 1544; https://doi.org/10.3390/cancers13071544
Submission received: 2 February 2021 / Revised: 15 March 2021 / Accepted: 24 March 2021 / Published: 27 March 2021
(This article belongs to the Special Issue Preclinical and Translational Research in Gynecological Cancers)

Abstract

:

Simple Summary

Epithelial ovarian cancers (EOC) have an unpredictable frequent recurrence often associated with incurable chemo-resistant disease. Basing on the miRNA expression profile of 892 EOC patients, we previously developed a 35 miRNA-based classifier, MiROvaR, able to predict EOC risk of early relapse. Further independent analysis of prediction accuracy represents a crucial step in the test-validation phase. Here we exploited an external and independently collected, handled and profiled EOC cohort, to challenge MirovaR accuracy. Our analysis confirmed the MiROvaR prognostic power, thus opening the way to its prospective validation as a clinical grade assay entering into clinical practice to help in the refinement of therapeutic intervention for high risk EOC patients.

Abstract

Epithelial ovarian cancer (EOC) remains the second most common cause of gynecological cancer deaths. To improve patients’ outcomes, we still need reliable biomarkers of early relapse, of which external independent validation is a crucial process. Our previously established prognostic signature, MiROvaR, including 35 microRNAs (miRNA) able to stratify EOC patients for their risk of relapse, was challenged on a new independent cohort of 197 EOC patients included in the Pelvic Mass Study whose miRNA profile was made publically available, thus resulting in the only accessible database aside from the EOC TCGA collection. Following accurate data matrix adjustment to account for the use of different miRNA platforms, MiROvaR confirmed its ability to discriminate early relapsing patients. The model’s original cutoff separated 156 (79.2%) high- and 41 (20.8%) low-risk patients with median progression free survival (PFS) of 16.3 months and not yet reached (NYR), respectively (hazard ratio (HR): 2.42–95% Confidence Interval (CI) 1.49–3.93; Log-rank p = 0.00024). The MiROvaR predictive accuracy (area under the curve (AUC) = 0.68; 95% Cl 0.57–0.79) confirms its prognostic value. This external validation in a totally independently collected, handled and profiled EOC cohort suggests that MiROvaR is a strong and reliable biomarker of EOC early relapse, warranting prospective validation.

1. Introduction

Epithelial ovarian cancer (EOC) is a life-threatening disease characterized by late-stage presentation, high pathological and molecular heterogeneity and frequent progression to an incurable state of chemo-resistant recurrent disease [1].
In the last 10 years, great efforts have been made to better characterize EOC to improve patients’ stratification. Broad international, collaborative studies have been performed primarily on high-grade serous EOC (HGSOC) to derive transcriptome-based molecular subtypes to guide patients’ management. However, the great genomic and spatial heterogeneity of HGSOCs posed a serious limitation to their clinical application. Yet, the recently defined ProTYPE [2] HGSOC molecular classifier, although it has reached the clinical grade level, still needs prospective validation without clear prognostic potential. In fact, with the exception of BRCA1/2 mutations guiding the use of PARP inhibitors [3], no transcriptome-based molecular classifiers have entered clinical practice.
The identification of high-risk EOC patients still remains an urgent clinical need to improve the design of tailored therapy. By relying on the master layer of regulation for gene expression provided by microRNAs (miRNAs) [4], we identified a 35 miRNA-based molecular predictor—MiROvaR—able to stratify EOC patients for their risk of relapse independently from the two strongest prognostic clinical variables so far available for EOC: the International federation of Gynecology and Obstetrics (FIGO) stage and residual disease after primary surgery [5]. The molecular predictor was developed by profiling 179 patients from the Multicenter Italian Trials in Ovarian Cancer (MITO)2 clinical trial [6], and validated in two independent cohorts: the TCGA collection of EOC [7] and a set of patients from our laboratory. We then demonstrated its ability to encompass the biological and molecular differences among the histological subtypes of EOC.
Here, we challenge MiROvaR performance in external and independent EOC Danish case material profiled for the identification of prognostic/predictive miRNA [8].

2. Materials and Methods

2.1. Data Processing

A total of 197 formalin-fixed, paraffin-embedded specimens from patients with EOC belonging to the “Pelvic Mass Study,” recruited at the Gynecologic Department, Rigshospitalet (Copenhagen University Hospital, Denmark), were profiled for microRNA expression [8] using GeneChip 1.0 miRNA microarrays (Affymetrix) that allowed the detection of 847 different human miRNAs. Normalized data were retrieved from NCBI’s Gene Expression Omnibus (GEO) repository under the accession number GSE94320. To account for the different miRNA platforms used, Affymetrix in the “Pelvic Mass Study” [8] and Agilent/Illumina in our paper [5], Prahm’s data [8] were processed using the justRMA function [9] to apply the RMA procedure. The intensity measures were background adjusted and normalized by the quantile method, with resulting expression levels on a log2-based scale [10]. Based on the 385 miRNAs shared among platforms in our previous work [5], 382 (99.2%) were present in the Affymetrix arrays. To assess MiROvaR on GSE94320, this data matrix and those used by Bagnoli et al. [5] were adjusted to reduce the likelihood of systemic, non-biological, technical experimental biases by ComBat algorithm [11].
The rescaled data distribution, resulting in skewness of −0.756 and kurtosis of 2.93, enabled the application of the MiROvaR model with threshold for patients’ stratification in risk classes, as defined in Bagnoli et al. [5]. Skewness is a measure of the symmetry of the data distribution, while kurtosis measures the tail-heaviness of the distribution. Skewness and kurtosis were computed by D’Agostino and Anscombe-Glynn tests, respectively [12].

2.2. Statistical Analysis

The clinical endpoint of the study is progression-free survival (PFS), defined as the time from primary surgery until relapse, progressive disease or death of any cause, whichever occurred first. All clinical data were retrieved from Prahm’s paper [8].
For biomarker optimal cutoff determination, we used the Cutoff Finder R package available at http://molpath.charite.de/cutoff (accessed on 15 January 2021) [13], with optimization of the correlation with PFS data; in this case, the survival analysis was performed using the coxph and survfit functions of the survival R package [14]. A Cox proportional hazard model was fitted to the dichotomized variable and the survival data. The point with the most significant (log-rank test) split was considered as being the optimal cutoff. Hazard ratios (HRs) including 95% confidence intervals were calculated and plotted in the range of biomarker expression.
To assess the overall performance of our model, the prediction error [15] was evaluated, defined through Brier’s score as a function of survival time and computed using the R package pec [16]. The prediction error established the relative worth of prediction risk obtained by our signature over time, and was computed by including our miRNA model and fitting a likelihood Cox-proportional hazard model. With this tool, the predicted risks of progression are grouped according to the deciles of their distribution, and for each decile, the observed proportion of an event is plotted against the mean value predicted by the miRNA model. The null model corresponds to a model wherein no data on the covariates are used. To evaluate the prediction error curves, the miRNA model and the null model are plotted over time. The benchmark value of 0.25 reached by the null model corresponds to a 50% risk prediction probability for every sample. The leave-one-out cross validation estimate of the prediction error was calculated for all event times, using a time-dependent adaption of the Brier score.
To evaluate a prediction model that gives a continuous range of probabilities, discrimination and calibration are recommended [17]. Discrimination was evaluated by time-dependent receiver operating curve (ROC) [18], computed using the timeROC R package [19]. With the information of follow-up and PFS, the AUC can be computed at several time points, and MirOvaR prognostic ability can be evaluated over time.
All tests were two-sided with a p-value considered significant at <0.05, and were performed using R software, version 3.6.0 (http://www.r-project.org/ accessed on 15 January 2021).

3. Results and Discussion

The Danish case material described by Prahm et al. [8] includes 197 EOC patients, representative of EOC demography and incidence, with a median follow-up time of 88 months, a median progression-free survival (PFS) of 19 months and a median overall survival of 49 months (clinical-pathological characteristics have been already detailed [8] and summarized in Table S1). As reported by Prahm et al. [8], the patients’ classification according to MiROvaR prognostic risk index, using the original cutoff defined in our paper [5], was inapplicable due to differences in miRNA dynamic ranges across the microarray platforms used. However, as suggested by the authors [8], when the median MiROvaR value was imposed as a cutoff, a significant prognostic stratification of EOC according to their risk of relapse was obtained (Figure S1).
To challenge the MiROvaR model on this independent case material, data had to be adjusted taking into account the different methodologies used in Prahm’s [8] and our [5] studies, which were based on Affymetrix and Agilent/Illumina platforms, respectively. Essentially, following background subtraction and normalization, 382 (99.2%) of the 385 miRNAs shared among our Agilent/Illumina platforms were also detected in the Affymetrix arrays. The data were further adjusted to reduce the likelihood of systemic, non-biological, technical experimental biases by the ComBat algorithm [9]. All the 35 miRNAs belonging to MiROvaR signature were detected on the Affymetrix arrays. The rescaled data distribution allowed for the assessment of the prognostic potential of MiROvaR by imposing the original model cutoff (0.07359) that stratified the patients into 156 (79.2%) high- and 41 (20.8%) low-risk (Figure 1A) categories. The model separated high- and low-risk patients with a median PFS of 16.3 months and not yet reached (NYR), respectively (HR: 2.42–95% CI 1.49–3.93; Log-rank p = 0.00024) (Figure 1B). Moreover, when the MiROvaR index was plotted against HR, the optimal cutoff value (the cutoff point with the most significant HR as determined by applying the Cutoff Finder tool) [10] corresponded to 0.1085, close to the published model cutoff (Figure 1C). With the optimal cutoff, only a marginal improvement was observed as compared to the original model cutoff (see Figure 1B); in fact, the best cutoff separated 152 (78.2%) high- and 45 (22.8%) low-risk patients (Figure S2A), with median progression-free survival (PFS) of 16 months and not yet reached (NYR), respectively (HR: 2.65–95% Cl 1.65–4.27; Log-rank p = 0.00003) (Figure S2B).
The traditional approach to assess the performance of a prognostic model in independent external datasets is to quantify how close predictions are to the actual outcome [20], assessing measures of overall performance and measures of discrimination (sensitivity and specificity). For overall performance, we evaluated the prediction error through the Brier score. As shown in Figure 1D, the score for MiROvaR (red line) is lower than the null model, denoting superior model performance over the period of follow-up [21]. The accuracy of the model was assessed by generation of the ROC curve and evaluation of AUC that reaches a value of 0.68; 95% CI 0.57–0.79 (Figure 1E). For the sake of completeness, a time-dependent ROC curve was computed measuring the performance of our model over time (Figure S3).
Our primary aim was to validate the performance of our miRNA-based molecular predictor of EOC early relapse. Extending our analysis to clinically relevant covariates available in Prahm’s dataset, a significant association of MiROvaR high risk was observed with serous histotype, high grade, and the presence of residual disease at primary surgery (Table 1). MiROvaR is intended to be a widely useful tool that could encompass the biological/molecular differences among the histological sub-types of EOC. However, we confirmed its validity in a homogeneous sub-group of HGSOC (from the TCGA dataset and from the second validation set used in our original paper [5]). When considering HGSOC only in the Danish case material, the proportion of low-risk HGSOC was too low to have a statistically significant impact; however, the median PFS for the MiROvaR low-risk group was 36 months (50% recurrence rate) versus 15 months (83% recurrence rate) for the MiROvaR high-risk group.
Although in the Danish case material the categorization of the patient population for FIGO stage was described (Table S1), this clinical parameter was not associated with the single-patient molecular profile, and therefore it cannot be included into the association analysis (see Table 1).
The lack of annotation of FIGO stage among the characteristics included into the matrix of molecular data effected also the possibility of assessing the prognostic independence of MiROvaR, adjusting for the clinical parameters that we used in our former study [5], i.e., residual disease and FIGO stage. Eventually, in a bivariate analysis adjusting for residual disease (the only prognostic clinical covariate here available), MiROvaR confirmed its independent prognostic value (Table 2).
The proposed categorization for residual disease (suboptimal vs. optimal) has been included into the analysis, since the training set used for MiROvaR development was derived from a MITO trial [6] that used this categorization. However, we clearly showed that in the three independent multicenter case materials used for its development and totally accounting for over 890 patients, MiROvaR also performed with the residual disease categorized as NED (not evident disease) vs. <1 cm vs. >1 cm [5]. A post-hoc analysis performed categorizing those patients according to current guidelines for residual disease (presence of residual vs. NED) again confirmed MiROvaR independency. In the Prahm case material, when adopting this classification, we observed a significant interaction between MiROvaR and residual disease (p = 0.016). However, it is worth noting that in patients with no residual disease, MiROvaR significantly and efficiently discriminated those with worse prognosis (HR 2.70; 95%CI 1.26–5.79; p-value 0.011). This can clearly have clinical relevance, since those patients—according to a prognostic stratification based on residual disease and without a further molecular classification—would all have been considered as exhibiting good prognosis.
The management of ovarian cancer is evolving from a “one-size-fits-all” approach to more precise interventions that take into account the tumors’ molecular characteristics [22]. These tailored approaches need the discovery of prognostic/predictive biomarkers guiding the selection of patients who are most likely to benefit from further therapeutic interventions, or a de-escalation to reduce unnecessary toxicity. While a clinical-grade assay for the transcriptome-based classification of HGSOC histotype is ready for “the prime time” [2], the new knowledge derived by single cell sequencing approaches adds a further level of complexity that suggests a re-interpretation, and actually limits the efficacy of transcriptome-based discrete HGSOC subtyping [23]. We have shown the potential role of miRNAs in predicting disease progression [5,24], and the Prahm et al. study [8] enabled us to test our MiROvaR model in a totally independently collected, handled and profiled EOC patient dataset, thus following the rules for biomarkers validation [25], and confirming MiROvaR reliability as a biomarker of EOC early relapse, regardless of their histotypes. A biomarker such as MiROvaR, based on a tumor’s molecular characteristics at diagnosis, once validated also on liquid biopsy, might help in patients’ selection before any therapeutic treatment. Furthermore, by giving information on the biological characteristics of the tumor, MiROvaR may help in identifying new actionable targets. From a clinical point of view, the subgroup of patients with unfavorable prognosis identified by MiROvaR might be a candidate for more aggressive treatment modalities. In the late stage, setting MiROvaR may help in classifying HR-proficient patients at high risk of relapse who could benefit from Bevacizubam treatment/maintenance rather than PARP inhibitors. In the early stage setting, MiROvaR could help in identifying a group of patients at higher risk of relapse, who would really benefit from the adjuvant chemotherapy, and a group of patients at lower risk for whom unnecessary treatments can be avoided.

4. Conclusions

At a methodological level, we should highlight the importance of publically available genomics data with associated, well-controlled clinical data to enable their wise reuse with new bioinformatic tools. Furthermore, at the clinical level, we believe that the herein presented data represent an important step forward in establishing the role of miRNAs as biomarkers in ovarian cancer, thus warranting a prospective validation hopefully resulting in MiROvaR entering into clinical practice.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers13071544/s1. Figure S1: Progression-free survival of EOC patients in Prahm’s dataset1 (GSE94320) stratified by risk according to MiROvaR median cutoff, Figure S2: Progression-free survival of EOC patients in Prahm’s dataset 1 (GSE94320) stratified by risk according to MiROvaR best cutoff, Figure S3: MiROvaR predictive accuracy over-time, Table S1: Clinical and pathological characteristics of patients from the Danish case material.

Author Contributions

All the authors have participated to the preparation of this manuscript. L.D.C., D.M. designed the study; L.D.C., P.C., M.B., S.P. and D.M. analyzed and interpreted the data; D.M., M.B., S.P. provided funding; L.D.C., D.M., wrote the first version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially funded by Fondazione Italiana per la Ricerca sul Cancro (AIRC IG 17475 to DM; AIRC IG 19243 to DM; AIRC IG 5776 to SP); CARIPLO Foundation 2013-0865 to DM; Italian Ministry of Health (RF 2016-02363995 to DM and 5 × 1000 Funds–2013 to MB). The funders had no role in: the design of the study; the collection, analysis, or interpretation of the data; the writing of the manuscript; or the decision to submit the manuscript for publication.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

This is an in-silico analysis of publically available data; please refer to the original publication [8].

Data Availability Statement

Normalized data for the reported analysis were retrieved from NCBI’s Gene Expression Omnibus (GEO) repository under the accession number GSE94320.

Acknowledgments

The authors would like to thank the network of the Danish Pelvic Mass Study for making the microarray and clinical data freely available to the scientific community allowing the independent evaluation of MiROvaR. The authors are grateful to the MITO translational group for the fruitful discussion and support on prospective validation of MiROvaR and, eventually, for its entering into clinical practice. The authors would like to thank Silvana Canevari for her critical revision of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The authors hadn’t any personal, economic, or scientific involvement with the team of the Danish Pelvic Mass Study, which generated the data used in the present paper.

References

  1. Lheureux, S.; Gourley, C.; Vergote, I.; Oza, A.M. Epithelial ovarian cancer. Lancet 2019, 393, 1240–1253. [Google Scholar] [CrossRef] [Green Version]
  2. Talhouk, A.; George, J.; Wang, C.; Budden, T.; Tan, T.Z.; Chiu, D.S.; Kommoss, S.; Leong, H.S.; Chen, S.; Intermaggio, M.P.; et al. Development and Validation of the Gene Expression Predictor of High-grade Serous Ovarian Carcinoma Molecular SubTYPE (PrOTYPE). Clin. Cancer Res. 2020, 26, 5411–5423. [Google Scholar] [CrossRef]
  3. Dizon, D.S. PARP inhibitors for targeted treatment in ovarian cancer. Lancet 2017, 390, 1929–1930. [Google Scholar] [CrossRef]
  4. Krol, J.; Loedige, I.; Filipowicz, W. The widespread regulation of microRNA biogenesis, function and decay. Nat. Rev. Genet. 2010, 11, 597–610. [Google Scholar] [CrossRef]
  5. Bagnoli, M.; Canevari, S.; Califano, D.; Losito, S.; Maio, M.D.; Raspagliesi, F.; Carcangiu, M.L.; Toffoli, G.; Cecchin, E.; Sorio, R.; et al. Development and Validation of a microRNA-Based Signature (MiROvaR) to Predict Early Relapse or Progression of Epithelial Ovarian Cancer: A Cohort Study. Lancet Oncol. 2016, 17, 1137–1146. [Google Scholar] [CrossRef]
  6. Pignata, S.; Scambia, G.; Ferrandina, G.; Savarese, A.; Sorio, R.; Breda, E.; Gebbia, V.; Musso, P.; Frigerio, L.; Del Medico, P.; et al. Carboplatin Plus Paclitaxel Versus Carboplatin Plus Pegylated Liposomal Doxorubicin as First-Line Treatment for Patients With Ovarian Cancer: The MITO-2 Randomized Phase III Trial. J. Clin. Oncol. 2011, 29, 3628–3635. [Google Scholar] [CrossRef] [PubMed]
  7. Kang, J.; D’Andrea, A.D.; Kozono, D. A DNA Repair Pathway–Focused Score for Prediction of Outcomes in Ovarian Cancer Treated With Platinum-Based Chemotherapy. J. Natl. Cancer Inst. 2012, 104, 670–681. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Prahm, K.P.; Høgdall, C.; Karlsen, M.A.; Christensen, I.J.; Novotny, G.W.; Høgdall, E. Identification and validation of potential prognostic and predictive miRNAs of epithelial ovarian cancer. PLoS ONE 2018, 13, e0207319. [Google Scholar] [CrossRef] [Green Version]
  9. Gautier, L.; Cope, L.; Bolstad, B.M.; Irizarry, R.A. affy--analysis of Affymetrix GeneChip data at the probe level. Bioinformatics 2004, 20, 307–315. [Google Scholar] [CrossRef]
  10. Huber, W.; Carey, V.J.; Gentleman, R.; Anders, S.; Carlson, M.; Carvalho, B.S.; Bravo, H.C.; Davis, S.; Gatto, L.; Girke, T.; et al. Orchestrating high-throughput genomic analysis with Bioconductor. Nat. Methods 2015, 12, 115–121. [Google Scholar] [CrossRef] [PubMed]
  11. Johnson, W.E.; Li, C.; Rabinovic, A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 2006, 8, 118–127. [Google Scholar] [CrossRef]
  12. Wessa, P. Free Statistics Software, Office for Research Development and Education, Version 1.2.1, URL. 2020. Available online: https://www.wessa.wet/ (accessed on 15 January 2021).
  13. Budczies, J.; Klauschen, F.; Sinn, B.V.; Győrffy, B.; Schmitt, W.D.; Darb-Esfahani, S.; Denkert, C. Cutoff Finder: A Comprehensive and Straightforward Web Application Enabling Rapid Biomarker Cutoff Optimization. PLoS ONE 2012, 7, e51862. [Google Scholar] [CrossRef] [Green Version]
  14. Therneau, T.; Lumley, T. Survival: Survival Analysis, Including Penalized Likelihood. 2011. Available online: Http://CRAN.R-Project.Org/Package=survival (accessed on 15 January 2021).
  15. Gerds, T.A.; Schumacher, M. Efron-Type Measures of Prediction Error for Survival Analysis. Biometrics 2007, 63, 1283–1287. [Google Scholar] [CrossRef]
  16. Mogensen, U.B.; Ishwaran, H.; Gerds, T.A. Evaluating Random Forests for Survival Analysis Using Prediction Error Curves. J. Stat. Softw. 2012, 50, 1–23. [Google Scholar] [CrossRef] [PubMed]
  17. Collins, G.S.; Reitsma, J.B.; Altman, D.G.; Moons, K.G. Transparent Reporting of a multivariable prediction model for Individual Prognosis or Diagnosis (TRIPOD). Ann. Intern. Med. 2015, 162, 735–736. [Google Scholar] [CrossRef] [PubMed]
  18. Heagerty, P.J.; Lumley, T.; Pepe, M.S. Time-Dependent ROC Curves for Censored Survival Data and a Diagnostic Marker. Biometrics 2000, 56, 337–344. [Google Scholar] [CrossRef] [PubMed]
  19. Blanche, P.; Dartigues, J.-F.; Jacqmin-Gadda, H. Estimating and comparing time-dependent areas under receiver operating characteristic curves for censored event times with competing risks. Stat. Med. 2013, 32, 5381–5397. [Google Scholar] [CrossRef]
  20. Steyerberg, E.W.; Vickers, A.J.; Cook, N.R.; Gerds, T.; Gonen, M.; Obuchowski, N.; Pencina, M.J.; Kattan, M.W. Assessing the Performance of Prediction Models: A Framework for Traditional and Novel Measures. Epidemiology 2010, 21, 128–138. [Google Scholar] [CrossRef] [Green Version]
  21. Schumacher, M.; Graf, E.; Gerds, T. How to Assess Prognostic Models for Survival Data: A Case Study in Oncology. Methods Inf. Med. 2003, 42, 564–571. [Google Scholar] [CrossRef] [PubMed]
  22. Lheureux, S.; Msc, M.B.; Oza, A.M. Epithelial ovarian cancer: Evolution of management in the era of precision medicine. CA A Cancer J. Clin. 2019, 69, 280–304. [Google Scholar] [CrossRef] [Green Version]
  23. Geistlinger, L.; Oh, S.; Ramos, M.; Schiffer, L.; LaRue, R.S.; Henzler, C.M.; Munro, S.A.; Daughters, C.; Nelson, A.C.; Winterhoff, B.J.; et al. Multiomic Analysis of Subtype Evolution and Heterogeneity in High-Grade Serous Ovarian Carcinoma. Cancer Res. 2020, 80, 4335–4345. [Google Scholar] [CrossRef] [PubMed]
  24. Cramer, D.W.; Elias, K.M. A prognostically relevant miRNA signature for epithelial ovarian cancer. Lancet Oncol. 2016, 17, 1032–1033. [Google Scholar] [CrossRef]
  25. Hentze, J.L.; Høgdall, C.; Kjær, S.K.; Blaakær, J.; Høgdall, E. Searching for new biomarkers in ovarian cancer patients: Rationale and design of a retrospective study under the Mermaid III project. Contemp. Clin. Trials Commun. 2017, 8, 167–174. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Progression-free survival (PFS) of epithelial ovarian cancers (EOC) patients in Prahm’s dataset (GSE94320) stratified by risk according to MiROvaR model. (A). MiROvaR index. Affymetrix microarray data were retrieved from Gene Expression Omnibus (GEO) repository and MiROvaR index was computed in GSE94320 after adjustment to account for the microarray platforms. The bar-plot depicts the MiROvaR index showing skewness = −0.756 and kurtosis = 2.93. The red bar shows the model cutoff value (=0.07359) as determined in our original paper and used in the present analysis. (B). Kaplan–Meier curves according to the MiROvaR value as cutoff: blue and red lines indicate low- and high-risk patients reaching HR = 1.42 (CI 1.49–3.93), p = 0.00024. High- and low-risk curves were compared with the long-rank test. HR = hazard ratio. Shadows indicate upper and lower 95% confidence intervals. (C). Hazard ratio assessed with PFS as the endpoint and independent of the cutoff point for the MiROvaR index. The red bar shows the model cutoff value (=0.07359), while the green bar designates the optimal cutoff point (=0.1085) with the most significant HR split. Solid and broken lines indicate the HR and the 95% confidence intervals. (D) Curves of time-dependent prediction errors by Brier scores to evaluate MiROvaR performance for predicting PFS in EOC. The Brier score for our model (red line) is computed along with the reference (i.e., the marginal Kaplan–Meier estimator, ignoring the predictors). (E). Time-dependent receiver operating curve (ROC) in Prahm’s dataset for MirOvaR predicting 72-month time point. The AUC is = 0.68 (95% CI 0.57–0.79).
Figure 1. Progression-free survival (PFS) of epithelial ovarian cancers (EOC) patients in Prahm’s dataset (GSE94320) stratified by risk according to MiROvaR model. (A). MiROvaR index. Affymetrix microarray data were retrieved from Gene Expression Omnibus (GEO) repository and MiROvaR index was computed in GSE94320 after adjustment to account for the microarray platforms. The bar-plot depicts the MiROvaR index showing skewness = −0.756 and kurtosis = 2.93. The red bar shows the model cutoff value (=0.07359) as determined in our original paper and used in the present analysis. (B). Kaplan–Meier curves according to the MiROvaR value as cutoff: blue and red lines indicate low- and high-risk patients reaching HR = 1.42 (CI 1.49–3.93), p = 0.00024. High- and low-risk curves were compared with the long-rank test. HR = hazard ratio. Shadows indicate upper and lower 95% confidence intervals. (C). Hazard ratio assessed with PFS as the endpoint and independent of the cutoff point for the MiROvaR index. The red bar shows the model cutoff value (=0.07359), while the green bar designates the optimal cutoff point (=0.1085) with the most significant HR split. Solid and broken lines indicate the HR and the 95% confidence intervals. (D) Curves of time-dependent prediction errors by Brier scores to evaluate MiROvaR performance for predicting PFS in EOC. The Brier score for our model (red line) is computed along with the reference (i.e., the marginal Kaplan–Meier estimator, ignoring the predictors). (E). Time-dependent receiver operating curve (ROC) in Prahm’s dataset for MirOvaR predicting 72-month time point. The AUC is = 0.68 (95% CI 0.57–0.79).
Cancers 13 01544 g001
Table 1. Distribution of MiROvaR high- and low-risk patients from the Danish case material in relation to clinical and pathological variables.
Table 1. Distribution of MiROvaR high- and low-risk patients from the Danish case material in relation to clinical and pathological variables.
Clinical CharacteristicsLow RiskHigh Riskp Value
N%N%
Age, years
Median6066
Range31–8131–89
Histology <0.0001
Serous221414086
Endometroid747853
Mucinous873327
Clear Cells444556
Grade <0.0001
1, well differentiated1155945
2, moderately differentiated22228078
3, poorly differentiated8116689
Missing information 1
Surgical debulking 0.002
Optimal (<1 cm)35289172
Suboptimal (>1 cm)686592
Table 2. Cox proportional hazard regression analysis.
Table 2. Cox proportional hazard regression analysis.
CovariatesUnivariate AnalysisMultivariate Analysis
HR (95% CI)p-ValueHR (95% CI)p-Value
MiROvaR (high- vs. low-risk)2.42 (1.49–3.93)0.0003671.75 (1.1–2.89)0.0282
Residual disease (suboptimal vs. optimal)4.28 (3–6.1)<0.00013.82 (2.65–5.49)<0.0001
HR, Hazard ratio; CI, confidence interval.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

De Cecco, L.; Bagnoli, M.; Chiodini, P.; Pignata, S.; Mezzanzanica, D. Prognostic Evidence of the miRNA-Based Ovarian Cancer Signature MiROvaR in Independent Datasets. Cancers 2021, 13, 1544. https://doi.org/10.3390/cancers13071544

AMA Style

De Cecco L, Bagnoli M, Chiodini P, Pignata S, Mezzanzanica D. Prognostic Evidence of the miRNA-Based Ovarian Cancer Signature MiROvaR in Independent Datasets. Cancers. 2021; 13(7):1544. https://doi.org/10.3390/cancers13071544

Chicago/Turabian Style

De Cecco, Loris, Marina Bagnoli, Paolo Chiodini, Sandro Pignata, and Delia Mezzanzanica. 2021. "Prognostic Evidence of the miRNA-Based Ovarian Cancer Signature MiROvaR in Independent Datasets" Cancers 13, no. 7: 1544. https://doi.org/10.3390/cancers13071544

APA Style

De Cecco, L., Bagnoli, M., Chiodini, P., Pignata, S., & Mezzanzanica, D. (2021). Prognostic Evidence of the miRNA-Based Ovarian Cancer Signature MiROvaR in Independent Datasets. Cancers, 13(7), 1544. https://doi.org/10.3390/cancers13071544

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop