A reactive oxygen species scoring system predicts cisplatin sensitivity and prognosis in ovarian cancer patients
BMC Cancer volume 19, Article number: 1061 (2019)
To reveal roles of reactive oxygen species (ROS) status in chemotherapy resistance and to develop a ROS scoring system for prognosis prediction in ovarian cancer.
We tested the sensitizing effects of ROS elevating drugs to cisplatin (cDDP) in ovarian cancer both in vitro and in vivo. A ROS scoring system was developed using The Cancer Genome Atlas (TCGA) database of ovarian cancer. The associations between ROS scores and overall survival (OS) were analyzed in TCGA, Tothill dataset, and our in-house dataset (TJ dataset).
ROS-inducing drugs increased cisplatin-induced ovarian cancer cell injury in vitro and in vivo. ROS scoring system was established using 25 ROS-related genes. Patients were divided into low (scores 0–12) and high (scores 13–25) score groups. Improved patient survival was associated with higher scores (TCGA dataset hazard ratio (HR) = 0.43, P < 0.001; Tothill dataset HR = 0.65, P = 0.022; TJ dataset HR = 0.40, P = 0.003). The score was also significantly associated with OS in multiple datasets (TCGA dataset r2 = 0.574, P = 0.032; Thothill dataset r2 = 0.266, P = 0.049; TJ dataset r2 = 0.632, P = 0.001) and with cisplatin sensitivity in ovarian cancer cell lines (r2 = 0.799, P = 0.016) when used as a continuous variable. The scoring system showed better prognostic performance than other clinical factors by receiver operating characteristic (ROC) curves (TCGA dataset area under the curve (AUC) = 0.71 v.s. 0.65, Tothill dataset AUC = 0.73 v.s. 0.67, TJ dataset AUC = 0.74 v.s. 0.66).
ROS status is associated with chemotherapy resistance. ROS score system might be a prognostic biomarker in predicting the survival benefit from ovarian cancer patients.
Ovarian cancer is the second common diagnosed and the most lethal of the various gynecologic malignancies . Although ovarian cancer patients are sensitive to platinum- and taxane-based chemotherapy during initial treatment, a significant proportion of patients relapse and develop platinum resistance [1, 2].
Reactive oxygen species (ROS) are oxygen-containing reactive chemical molecules generated during metabolic processes. ROS play an essential role in signal transduction pathways , cell cycle progression [3,4,5], gene transcription [3, 6], cell differentiation [7, 8], and cell death . Elevated oxidative stress and delicate redox balance were detected in cancer cells due to activation of oncogene, high metabolic activity, and mitochondrial malfunction [6, 9, 10]. Deprivation of the redox balance through an increase in ROS levels or a decrease in the cellular antioxidant capacity to induce cellular ROS burst have shown therapeutic benefits in cancer cells [11,12,13]. Most chemotherapeutics, including platinum and taxanes, exert anti-cancer effects by inducing ROS-mediated cell damage in cancer cells [14,15,16,17]. Since new therapeutic approaches combining chemotherapeutics with ROS-elevating drugs have exhibited improvement of the cytotoxicity and reduction of the resistance [18,19,20]. Moreover, some studies have shown that part of patients with drug resistance attribute to the lower level of tumor cell oxidative stress and stronger antioxidant ability [21, 22].
Therefore, it is of increasing interest to develop a prognostic method to predict patients who will benefit from chemotherapy or additive ROS inducer based on the quantifiable criteria of ROS status. In this paper, we reflected that ROS is involved in the drug resistance and chemo-sensitivity in vitro and in vivo. We then established a comprehensive scoring system through analyzing the relationship between the expressions of ROS pathway genes, including genes involved in oxidative stress, oxidation reaction, antioxidant response, and prognosis of patients in TCGA database. We validated the effect of this system in Tothill database and ovarian cancer patients from our hospital. The result showed that this scoring system might be clinically applied to predict the outcome of chemotherapy in ovarian cancer patients.
SKOV3 (HTB-77), Caov3 (HTB-75), OVCAR3 (HTB-161), and OV-90 (CRL-11732) ovarian cancer cell lines were purchased from the American Type Culture Collection (ATCC, Manassas, VA, US) and cultured as recommended. Caov3 was cultured in DMEM with 10% FBS (Invitrogen) and OVCAR3 was cultured in RPMI 1640 with 10% FBS (Invitrogen) and 10μg/ml insulin (Bovine). SKOV3 was cultured in McCoy 5A with 10% FBS (Invitrogen). OV-90 was grown in a 1:1 mixture of MCDB 105 medium containing a final concentration of 1.5 g/L sodium bicarbonate and Medium 199 containing a final concentration of 2.2 g/L sodium bicarbonate with 15% FBS (Thermo Scientific). Cisplatin-sensitive ovarian cancers cell line (OV2008) and its resistant variant (C13*) were gifts from Prof. Benjamin K. Tsang in the Ottawa Health Research Institute, Ottawa, Canada  and cultured in RPMI 1640 medium with 10% FBS (Invitrogen). All cells were free from mycoplasma and were used between 3 and 5 passages after thawing. All cells were authenticated by china center for type culture collection (CTCC, Wuhan, China) using short tandem repeat (STR) DNA profiling. Primary cell lines were isolated from ovarian cancer tissue specimens of patients undergoing surgical resection as previously described  and cultured in DMEM/F12 medium (Invitrogen) with 20% FBS (Invitrogen). The passage number of primary ovarian cancer cells during the experiment was ranged from 3 to 8 generations. All cell lines were cultured in a 37 °C humidified atmosphere containing 5% CO2.
Assessment of cell viability
Viability of cells were assessed by Cell Counting Kit-8 reagent (CCK8, Dojindo, Tokyo, Japan). 5 × 103 cancer cells were seed in 96-well plate and treated with cDDP at different concentrations for 48 h with or without ROS-elevating (PLX4032 (1 μM), Piperlongumine (PIPER, 10 μM) and β-phenylethyl isothiocyanate (PEITC, 10 μM)) or ROS-scavenging drugs (glutathione (GSH, 2 mM), N-acetyl cysteine (NAC, 1 mM) and Vitamin C (VitC, 1 mM)). Supernatants were removed and 100 μl of CCK8 solution (1:10 dilution) were added to the cancer cells. After 2 h incubation at 37 °C in dark, optical density (OD) at 450 nm was measured by a microplate reader. The IC50 value for each cell line was determined by nonlinear regression analysis using GraphPad Prism (GraphPad Software Inc., San Diego, CA). The results were tested by three independent experiments.
ROS measurement was assayed using dichloro-dihydro-fluorescein diacetate (DCFDA; Beyotime, Shanghai, China), according to the manufacturers’ instructions. Briefly, cells were loaded with DCFH-DA, washed with ice-cold HBSS. Then, the fluorescence intensity of the cells was measured at 488 nm by flow cytometry.
Tumor Xenograft studies
The study was approved by the Ethical Committee of the Medical Faculty of Tongji Medical College (Wuhan, China), and performed in accordance with the relevant guidelines and regulations. Six-week-old athymic female homozygous BALB/c nude mice (SPF) were bought from Beijing Hua Fukang biological Polytron Technologies Inc., and reared in accordance with the relevant guidelines and regulations. C13* cells (107 cells/0.1 ml PBS/mice) intraperitoneally injected into the right flank of BALB/C nude mice, under isoflurane-induced anesthesia. Treatment began when the tumor volume reached between 70 and 100mm3. The mice were randomly divided into 4 groups (n = 6) including PBS group, cDDP treatment group (cDDP 2.5 mg/kg, i.p., every 4 days for 28 days), PIPER treatment group (PIPER 2 mg/kg, i.p., daily for 28 consecutive days), and cDDP combination with PIPER treatment group (same dose as used in the single-agent groups). Following the initial treatment, the tumor sizes were measured every 2 days. Tumor volumes (V) were calculated by the following formula: V (mm3) = length × (square of width)/2. The mice were euthanized by cervical dislocation. Tumors were excised, weighed, paraformaldehyde-fixed paraffin embedding and used for ex vivo immunohistochemical staining.
Immunohistochemistry (IHC) and scoring
IHC staining was performed as described previously . Briefly, tissue sections were incubated with antibody γ-H2AX (Abcam, Cat: ab2893, dilution 1:200), Ki67 (Abcam, Cat: ab15580, dilution 1:200), CD34 Mouse monoclonal (Abcam, Cat: ab198395, 1:1000), and cleaved Caspase-3 (Cell signaling Technology, Cat: 9661,1:200) overnight at 4 °C and stained by 3,3′-diaminobenzidine (DAB). Tumor-cell staining was assigned a score as described previously . All specimens were evaluated by two independent experts simultaneously.
Study design, patients, and sample processing
This study was designed using a discovery stage and validation phase. In the discovery stage, 511 SOC patients with level 3 mRNA data were obtained from the TCGA database  to establish a scoring system. In the validation phase, the scores were validated using the largest outside independent dataset-Tothill dataset (GSE9899, n = 285). Patients lacking the serous pathologic type (n = 45) were excluded from the Tothill dataset. To further validate the scoring system, 105 blocks of formalin-fixed, paraffin-embedded (FFPE) tissues from primary epithelial ovarian cancer were obtained. The study was approved by the Ethical Committee of the Medical Faculty of Tongji Medical College. All patients written informed consents. The surgical staging was assessed in accordance with the International Federation of Gynecology and Obstetrics (FIGO) classification. Optimal debulking was defined as ≤1 cm residual disease. All clinicopathological characteristics are reported in Additional file 1: Table S1.
Quantitative real-time PCR (qRT-PCR)
RNAs from 105 FFPE cases were extracted from four 10-μm-thick FFPE sections using the miRNeasy FFPE kit (Qiagen, Valencia, CA, USA). The cDNA was synthesized by the SuperScript® IV First-Strand Synthesis System (Thermo Fisher Scientific, China). Real-time PCR amplification was performed on an CFX Connect™ Real-Time PCR Detection System with SYBR reagent (Bio-Rad, China). GAPDH was used as an internal control.
Student’s t-test was performed to compare the statistical difference between two groups. Multiple comparisons were accessed using the one-way analysis of variance (ANOVA). Survival was analyzed by the Kaplan–Meier method with the log-rank test. Univariate and multivariable Cox regression analyses were used to test for statistical independence between the score, pathological, and clinical variables. The relationship between the score and median OS was measured by pearson correlation coefficient. Area under the curve (AUC) values were calculated from the ROC curves. All tests were two-sided, and P-values < 0.05 were considered to indicate a statistically significant difference. All calculations were performed with SPSS (Version 25.0).
ROS levels are associated with cDDP sensitivity of ovarian cancer both in vitro and in vivo
To determine whether the ROS levels of ovarian cancer cells play a role in cisplatin (cDDP) resistance, 6 cell lines, including SKOV3, Caov3, OVCAR3, OV-90, OV2008, and C13* were exposed to different concentrations of cDDP (Fig. 1a and Additional file 2: Figure S1A). The IC50 of CDDP-treated cell lines are displayed in Table S2. We also measured the intrinsic ROS levels in each of the cell lines using flow cytometry (Additional file 2: Figure S1B). However, intrinsic ROS levels had no statistically significant impact on cDDP sensitivity. Then, the cancer cells were treated with different cDDP concentrations in addition to fixed concentrations of ROS-elevating (PLX4032 (1 μM), Piperlongumine (PIPER, 10 μM) and β-phenylethyl isothiocyanate (PEITC, 10 μM)) or ROS-scavenging drugs (glutathione (GSH, 2 mM), N-acetyl cysteine (NAC, 1 mM) and Vitamin C (VitC, 1 mM)). The doses of ROS-elevating or ROS-scavenging were preselected while mono-therapy has no effect on cancer cell proliferation or apoptosis (Additional file 2: Figure S1C). Interestingly, most ROS-elevating drugs increased cDDP cytotoxicity in all ovarian cancer cell lines, especially in C13*, OV2008 and SKOV3 (Fig. 1a, Additional file 2: Figure S1A, and Additional file 3: Table S2). By contrast, the cytotoxicity of cDDP in cell lines such as OV90 and OVCAR3 was reduced when combined with ROS-scavenging drugs (Additional file 2: Figure S1A and Additional file 3: Table S2). While Caov3 exhibited a moderate response in the presence of both ROS-scavenging and ROS-elevating drugs (Additional file 2: Figure S1A and Additional file 3: Table S2). Furthermore, the sensitizing effects of ROS-elevating drugs were confirmed in 6 primary cells derived from ovarian cancer patients (Fig. 1b, Additional file 2: Figure S1D and Additional file 3: Table S2), supporting its clinical relevance. PIPER showed a better sensitizing effect of cDDP cytotoxicity than other ROS-scavenging drugs, we further verified its sensitizing effect in another 5 primary cells. Five primary cell lines were derived from 3 patients with recurrent ovarian cancer and 2 patients with primary ovarian cancer. As expected, no matter what kind of ovarian cancer patients, PIPER increased cDDP cytotoxicity (Fig. 1c and Additional file 2: Figure S1E).
On the basis of the sensitizing effects of ROS-elevating drugs on cDDP in cell lines and primary ovarian cancer cells, we explored cDDP and PIPER combinations in C13* (cDDP resistant) xenograft tumors. As expected, C13* tumors are highly resistant to cDDP mono-therapy. Combination of PIPER and cDDP markedly delay tumor growth, while PIPER mono-therapy showed a minimal effect on tumor growth (Fig. 1d and e). IHC analysis showed that combination therapy diminished blood vessels (CD34), suppressed proliferation (Ki67) and increased DNA damage (H2ax) and apoptosis (cleaved caspase-3) compared to cDDP or PIPER mono-therapy (Fig. 1f).
These results revealed that baseline ROS levels in ovarian cancer cells measured by DCFDA do not accurately predict their sensitivity to cDDP. ROS-elevating drugs increased ovarian cancer cell sensitivity to cDDP in varying degrees. So, it is necessary to build a scoring system to assess ovarian cancer patients who may benefit from the combination of cDDP and ROS-elevating drugs.
Establishment of the ROS scoring system in ovarian cancer patients
ROS-elevating drugs sensitized ovarian cancer cells to cDDP, which to some extent depended on intrinsic ROS levels in the cells. However, there is no reliable and convenient methods to quantified ROS status in tumors. So, we developed a ROS scoring system based on expressions of ROS related genes. To identify ROS pathway genes, we selected 179 ROS related genes according to our knowledge and works of literature. Ovarian cancer patients in TCGA (n = 511) were divided into two groups according to the median expression values of the 179 ROS related genes individually. 25 of 179 ROS related genes were selected whose expression levels were associated (P < 0.15) with the OS of ovarian cancer patients in TCGA dataset (n = 511) (Table 1 and Additional file 4: Figure S2). For each patient, if high gene expression was associated with a good prognosis, all patients with a higher expression value than the median expression obtained one point, and vice versa. The scores of all candidate genes of each patient were added to obtain a total score, which was called the ROS score.
The score divided patients into two categories (“low” or “high”) with equal ranges. Kaplan-Meier and univariate Cox proportional hazard regression analyses revealed patients with high scores (scores 13–25) had better OS compared to patients with low scores (scores 0–12) (high v.s. low scores, median OS = 4.66 years v.s. 2.92 years, HR = 0.43; 95% CI = 0.34 to 0.55, P < 0.001) (Fig. 2a). Moreover, the ROS scoring system performed well independent of tumor stages, treatments, and different molecular subtypes. Patients with high ROS scores are associated with better prognosis in stages III/IV patients (HR = 0.44; 95% CI = 0.34 to 0.56, P < 0.001) (Fig. 2b), in stage III/IV patients receiving first-line chemotherapy with platinum and taxane regimens (HR = 0.38; 95% CI = 0.27 to 0.53, P < 0.001) (Fig. 2c), or in different ovarian cancer molecular subtypes (immunoreactive, differentiated, proliferative, mesenchymal) (Additional file 5: Figure S3).
Prognosis prediction value of the ROS scoring system in ovarian cancer patients
The association between ROS score and OS was confirmed in an independent outside dataset (Tothill, GSE9899), the largest ovarian cancer datasets available (high v.s. low scores, median OS = 4.77 years v.s. 3.04 years; HR = 0.65; 95% CI = 0.44 to 0.94, P = 0.022) (Fig. 3a). To further validate the prognosis prediction values of this ROS scoring system and broaden the scope of applications. We quantified these 25 survival-related ROS genes using a widely-applicable and convenient method, qRT-PCR analysis from FFPE tissues in 105 ovarian tumors. Consistent with afore-mentioned associations based on microarray data, patients with high ROS scores have longer survival times (high v.s. low, median OS = 3.76 years v.s. 2.36 years; HR = 0.40; 95% CI = 0.22 to 0.74, P = 0.002) (Fig. 3b). Univariate and multivariate Cox proportional hazard regression analyses were applied to estimate the relationship between OS and the score (high v.s. low) compared with other clinical factors, including age (≤59 v.s. ≥60 years), FIGO stage (IV v.s. III), histological grade (3 v.s. 1–2), and extent of surgical debulking (0–10 v.s. ≥11 mm residual tumor). Only the score is an independent prognostic factor in TCGA, Tothill dataset and our in house patients datasets (Fig. 3c, Additional file 5: Figure S4).
The contribution of the ROS scoring system as a continuous variable toward prediction of OS in all datasets and cDDP sensitivity in 6 ovarian cancer cell lines
To further verify the associations between scores and patients’ survival, we performed correlation analysis between each score and the median survival times in patients with the same scores. Interestingly, there was a positive correlation between scores and the median survival times (r = 0.758, P = 0.032 in TCGA dataset; r = 0.516, P = 0.049 in Tothill dataset; and r = 0.795, P = 0.001in TJ dataset) (Fig. 4).
Univariate and multivariate analyses were repeated when the ROS score was assessed as a continuous variable. Again, the ROS score outperformed other clinical covariates as an independent prognostic factor in TCGA dataset, Tothill dataset, and in TJ dataset (Table 2).
Predictive accuracy of the ROS scoring system in ovarian Cancer patients
To further evaluate the contribution of the score to OS prediction, ROC curve analysis was performed using the following variables: clinical covariates (age, grade, stage, and residual tumor (AGSR)); ROS score (Score); and clinical covariates plus ROS score (AGSR + Score).
For these analyses, patients were divided into two groups with a survival time higher or lower than the median OS (TCGA dataset median OS =3.79 years, Tothill dataset median OS = 3.62 years, TJ dataset median OS = 2.81 years), and those with a survival time shorter than the median OS at last follow-up were excluded. Interestedly, the score alone have higher predictive values than clinical covariates in all datasets (TCGA dataset AUC = 0.71 v.s. 0.65, Tothill dataset AUC = 0.73 v.s. 0.67, TJ dataset AUC = 0.74 v.s. 0.66) (Fig. 5a). Moreover, combined with clinical covariates and ROS score could further improve predictive performance in all dataset (TCGA dataset AUC = 0.73, Tothill dataset AUC = 0.81, TJ dataset AUC = 0.78) (Fig. 5b and c).
Compared with other types of cancers, one unique feature of ovarian cancer is that over 50% of ovarian cancers contain p53 mutation . Specifically, p53 mutation was identified in 96% of all serous ovarian tumors . Suppression of p53 led to significant decreases in the expression of SESN1, SESN2, and GPX1, suggesting that p53 is involved in cellular metabolism and antioxidant response [29,30,31]. So, p53 mutation could increase ROS levels and oxidative damage of DNA in ovarian cancer cells. Thus, alterations in the expression of ROS genes that affect ROS production or scavenging may be closely linked to the resistance of ovarian cancer cells to chemotherapy.
An increasing number of studies have identified relationships between ROS related genes (such as SOD, CAT, GLS2 and so on) and drug resistant [32, 33]. We found that ROS pathway function or activity plays a crucial role in chemotherapy responses in ovarian cancer cells and transplanted mouse models. Our results suggest that ROS related gene expression changes are important mechanisms by which ovarian cancer cells acquire resistance to anticancer drugs, and these changes result in different outcomes and prognoses in ovarian cancer patients.
In this study, we established a quantifiable ROS scoring system able to predict ovarian cancer patient prognosis based on the expression levels of ROS related genes in the TCGA dataset (n = 511). Moreover, we validated this system in another published dataset (Tothill dataset, GSE9899, n = 241). We further validated the scoring system in our in-house patient dataset (TJ dataset, n = 105). We indicated that the scoring system accurately determined the prognosis of ovarian cancer patients. The use of FFPE sections and qPCR also extended the use of the scoring system. Both datasets demonstrated that the system is prognostic for survival.
A number of gene profile-based prognosis techniques, used in combination with microarrays or PCR, have been previously developed to predict survival in patients with ovarian cancer [34,35,36], but the results have not been satisfactory. We demonstrated that our scoring system is superior to other known clinical factors in predicting OS, not only in the TCGA dataset but also in our dataset and another online validation set. TCGA divided ovarian cancer into four molecular subtypes (immunoreactive, differentiated, proliferative, mesenchymal) based on gene clustering, but these clusters did not associate with OS. However, statistical significance was observed when the score was applied to all subtypes except the proliferative subtype. Our score extends application of the TCGA classification model. Our system was also able to predict outcomes to first-line platinum and taxane chemotherapy in ovarian cancer. This feature has profound clinical significance because there are no other good clinical factors to predict the response to platinum-based standard chemotherapy. Most patients with advanced serous ovarian cancer will relapse after a few years even after standard therapies like thorough operation and chemotherapy are used. In addition, about 30% of patients with primary platinum resistance undergo multiple cycles of useless and potentially toxic treatment before second-line drug treatments are used. Moreover, agents that increase the ROS levels of cancer cells could be used as a standard treatment to improve chemotherapeutic responses for patients with ovarian cancers with low ROS levels.
In this study, we just demonstrate the prognostic value of the ROS scoring system, which leads to the possibility of clinical application. For individual ovarian cancer patient with poor prognosis predicted by the ROS scoring system, if possible, we can combine ROS-inducing drugs with platinum- and taxane-based chemotherapies to improve outcomes. We know there are lots of problems, such as the in vivo stability of ROS-elevating drugs, targeted property, and safety need to be resolved before ROS becoming a therapeutic target. However, further study based on patient-derived tumor xenograft (PDX) animal models for intraperitoneal administration of ROS-elevating drugs may lead to the possibility of clinical transformation.
This study has several limitations. First, although we reproduced our findings in two other datasets, this study is a retrospective analysis, and sample selection bias may exist. Of course, we hope that this score will be tested prospectively in a clinical trial, and we believe that the score is ready for such testing. Second, our gene expression profiling and analysis is only limited to ROS related pathways. However, other gene expression pathways that may be important in survival predictions were neglected. Third, this study is limited in its gene expression profiling. Other mechanisms of gene regulation, including microRNAs, DNA methylation, and CNV region changes were not considered. We are looking forward to future studies of this type and the development of more comprehensive prediction models.
We established a ROS scoring system that could predict the outcomes of ovarian cancer patients. This system offers considerable improvement over existing methods for prognostic classification and has the potential to provide clinicians with useful, readily available information for personalized chemotherapy in the future.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Reactive oxygen species
The Cancer Genome Atlas
International Federation of Gynecology and Obstetrics
Area under the curve.
Torre LA, Trabert B, DeSantis CE, Miller KD, Samimi G, Runowicz CD, Gaudet MM, Jemal A, Siegel RL. Ovarian cancer statistics, 2018. CA Cancer J Clin. 2018;68(4):284–96.
Cortez AJ, Tudrej P, Kujawa KA, Lisowska KM. Advances in ovarian cancer therapy. Cancer Chemother Pharmacol. 2018;81(1):17–38.
Zhang J, Wang X, Vikash V, Ye Q, Wu D, Liu Y, Dong W. ROS and ROS-mediated cellular signaling. Oxidative Med Cell Longev. 2016;2016:4350965.
Boonstra J, Post JA. Molecular events associated with reactive oxygen species and cell cycle progression in mammalian cells. Gene. 2004;337:1–13.
Patterson JC, Joughin BA, van de Kooij B, Lim DC, Lauffenburger DA, Yaffe MB. ROS and Oxidative stress are elevated in mitosis during asynchronous cell cycle progression and are exacerbated by mitotic arrest. Cell Syst. 2019;8(2):163–7 e162.
Kumari S, Badana AK, G MM, G S, Malla R: Reactive Oxygen Species: A Key Constituent in Cancer Survival Biomark Insights 2018, 13:1177271918755391.
Zhou D, Shao L, Spitz DR. Reactive oxygen species in normal and tumor stem cells. Adv Cancer Res. 2014;122:1–67.
Wang J, Luo B, Li X, Lu W, Yang J, Hu Y, Huang P, Wen S. Inhibition of cancer growth in vitro and in vivo by a novel ROS-modulating agent with ability to eliminate stem-like cancer cells. Cell Death Dis. 2017;8(6):e2887.
Liou GY, Storz P. Reactive oxygen species in cancer. Free Radic Res. 2010;44(5):479–96.
Sosa V, Moline T, Somoza R, Paciucci R, Kondoh H, ME LL: Oxidative stress and cancer: an overview. Ageing Res Rev 2013, 12(1):376–390.
Liu J, Wang Z. Increased oxidative stress as a selective anticancer therapy. Oxidative Med Cell Longev. 2015;2015:294303.
Postovit L, Widmann C, Huang P, Gibson SB. Harnessing oxidative stress as an innovative target for Cancer therapy. Oxidative Med Cell Longev. 2018;2018:6135739.
Leone A, Roca MS, Ciardiello C, Costantini S, Budillon A. Oxidative stress gene expression profile correlates with Cancer patient poor prognosis: identification of crucial pathways might select novel therapeutic approaches. Oxidative Med Cell Longev. 2017;2017:2597581.
Zhang J, Lei W, Chen X, Wang S, Qian W. Oxidative stress response induced by chemotherapy in leukemia treatment. Mol Clin Oncol. 2018;8(3):391–9.
Huang CY, Ju DT, Chang CF, Muralidhar Reddy P, Velmurugan BK. A review on the effects of current chemotherapy drugs and natural agents in treating non-small cell lung cancer. Biomedicine (Taipei). 2017;7(4):23.
Il'yasova D, Kennedy K, Spasojevic I, Wang F, Tolun AA, Base K, Young SP, Kelly Marcom P, Marks J, Millington DS, et al. Individual responses to chemotherapy-induced oxidative stress. Breast Cancer Res Treat. 2011;125(2):583–9.
Conklin KA. Chemotherapy-associated oxidative stress: impact on chemotherapeutic effectiveness. Integr Cancer Ther. 2004;3(4):294–300.
Zhang W, Hu X, Shen Q, Xing D. Mitochondria-specific drug release and reactive oxygen species burst induced by polyprodrug nanoreactors can enhance chemotherapy. Nat Commun. 2019;10(1):1704.
Molavian HR, Goldman A, Phipps CJ, Kohandel M, Wouters BG, Sengupta S, Sivaloganathan S. Drug-induced reactive oxygen species (ROS) rely on cell membrane properties to exert anticancer effects. Sci Rep. 2016;6:27439.
Wang J, Yi J. Cancer cell killing via ROS: to increase or decrease, that is the question. Cancer Biol Ther. 2008;7(12):1875–84.
Filippova M, Filippov V, Williams VM, Zhang K, Kokoza A, Bashkirova S, Duerksen-Hughes P. Cellular levels of oxidative stress affect the response of cervical cancer cells to chemotherapeutic agents. Biomed Res Int. 2014;2014:574659.
Liu Y, Li Q, Zhou L, Xie N, Nice EC, Zhang H, Huang C, Lei Y. Cancer drug resistance: redox resetting renders a way. Oncotarget. 2016;7(27):42740–61.
Asselin E, Mills GB, Tsang BK. XIAP regulates Akt activity and caspase-3-dependent cleavage during cisplatin-induced apoptosis in human ovarian epithelial cancer cells. Cancer Res. 2001;61(5):1862–8.
Li H, Gao Q, Guo L, Lu SH. The PTEN/PI3K/Akt pathway regulates stem-like cells in primary esophageal carcinoma cells. Cancer Biol Ther. 2011;11(11):950–8.
Sun CY, Su TF, Li N, Zhou B, Guo ES, Yang ZY, Liao J, Ding D, Xu Q, Lu H, et al. A chemotherapy response classifier based on support vector machines for high-grade serous ovarian carcinoma. Oncotarget. 2016;7(3):3245–54.
Goldman M, Craft B, Hastie M, Repečka K, McDade F, Kamath A, Banerjee A, Luo Y, Rogers D, Brooks AN et al: The UCSC Xena platform for public and private cancer genomics data visualization and interpretation 2019:326470.
Kmet LM, Cook LS, Magliocco AM. A review of p53 expression and mutation in human benign, low malignant potential, and invasive epithelial ovarian tumors. Cancer. 2003;97(2):389–404.
Cancer Genome Atlas Research N. Integrated genomic analyses of ovarian carcinoma. Nature. 2011;474(7353):609–15.
Wu Y, Lee S, Bobadilla S, Duan SZ, Liu X. High glucose-induced p53 phosphorylation contributes to impairment of endothelial antioxidant system. Biochim Biophys Acta Mol basis Dis. 2017;1863(9):2355–62.
Berkers CR, Maddocks OD, Cheung EC, Mor I, Vousden KH. Metabolic regulation by p53 family members. Cell Metab. 2013;18(5):617–33.
Mantovani F, Collavin L, Del Sal G. Mutant p53 as a guardian of the cancer cell. Cell Death Differ. 2019;26(2):199–212.
Stiewe T, Haran TE. How mutations shape p53 interactions with the genome to promote tumorigenesis and drug resistance. Drug Resist Updat. 2018;38:27–43.
Alam F, Mezhal F, El Hasasna H, Nair VA, Aravind SR, Saber Ayad M, El-Serafi A, Abdel-Rahman WM. The role of p53-microRNA 200-Moesin axis in invasion and drug resistance of breast cancer cells. Tumour Biol. 2017;39(9):1010428317714634.
Chiang AJ, Chen J, Chung YC, Huang HJ, Liou WS, Chang C. A longitudinal analysis with CA-125 to predict overall survival in patients with ovarian cancer. J Gynecol Oncol. 2014;25(1):51–7.
Paik ES, Lee JW, Park JY, Kim JH, Kim M, Kim TJ, Choi CH, Kim BG, Bae DS, Seo SW. Prediction of survival outcomes in patients with epithelial ovarian cancer using machine learning methods. J Gynecol Oncol. 2019;30(4):e65.
Prahm KP, Hogdall C, Karlsen MA, Christensen IJ, Novotny GW, Hogdall E. Identification and validation of potential prognostic and predictive miRNAs of epithelial ovarian cancer. PLoS One. 2018;13(11):e0207319.
This study is supported by the National Key Research and Development Program (2016YFC1303100), the National Basic Research Program of China (973 Program, 2015CB553903), the National Science-technology Supporting Plan Projects (2015BAI13B05), the Chinese National Key Plan of Precision Medicine Research (2016YFC0902901), and Nature and Science Foundation of China (81402163, 81402164, 81472783, 81572569, 81501530, 81671394, 81370469), the International S&T Cooperation Program of China (No. 2013DFA31400), and the Research Project of Health and Family Planning Commission of Hubei Province (WJ2015MA001). The funding bodies had no influence on the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Primary cell lines and animal studies were approved by the Ethical Committee of the Medical Faculty of Tongji Medical College (Wuhan, China), and were performed according to the relevant guidelines and regulations. All procedures performed in studies were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. Written informed consent was obtained from all individual participants included in the study.
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.
Clinicopathologic characteristics of ovarian cancer patients in 3 datasets. (docx 16.8 kb) (DOCX 16 kb)
ROS level is related to the survival of ovarian cancer cells treated with cDDP. (jpg 2.83mb) (A) Cell viability of ovarian cancer cell lines Caov3, OV90 and OVCAR3 was measured after treatment with gradient concentrations of cDDP with or without ROS-elevating (PLX4032, 1 μM, Piperlongumine (PIPER, 10 μM) and β-phenylethyl isothiocyanate (PEITC, 10 μM)) or scavenging drugs (glutathione (GSH, 2 mM), N-acetyl cysteine (NAC, 1 mM) and Vitamin C (VitC, 1 mM)) for 48 h by CCK-8. (B) Intracellular ROS concentrations of SKOV3, Caov3, OVCAR3, OV-90, OV2008, and C13* were measured by DCF-DA staining. (C) CCK8 detected cell viability of ovarian cancer cell lines after treatment with ROS-elevating. (D) cDDP IC50 curves for 3 strains of primary cancer cells with or without ROS-elevating or scavenging drugs for 48 h by CCK-8. (E) Cell viability of primary cancer cells derived from patients with recurrent or primary ovarian cancer was measured after treatment with gradient concentrations of cDDP with or without PIPER for 48 h by CCK-8. The two-tailed P-values < 0.05 were considered to indicate statistically significant differences. The results were tested by three independent experiments. (JPG 2907 kb)
The IC50 of ovarian cancer cell lines and primary cancer cells treated with different combinations of drugs. (docx 18.8 kb) (DOCX 18 kb)
84 of 179 ROS related genes with Kaplan-Meier log-rank P values < 0.5. (tif 328 kb) Genes with P < 0.15 are highlighted in red. (TIF 328 kb)
Prognostic Value of the ROS Scoring System in four ovarian cancer molecular subtypes of TCGA. (tif 495 kb) A Kaplan-Meier analysis of overall survival (OS) were performed in the differentiated, immunoreactive, mesenchymal, and proliferative for ovarian cancer patients in TCGA dataset with the ROS scoring system (high [scores 13–25], the green line v.s. low [scores 0–12], the blue line) is shown (P < 0.001, log-rank test). All statistical tests were two-sided. (TIF 495 kb)
Univariate analyses were performed in TCGA and Tothill dataset. (tif 1.21mb) Univariate analyses incorporating the score and known prognostic clinical factors, including age at diagnosis (≤59 v.s. ≥60 years), stage (III v.s. IV), grade (1–2 v.s. 3), and surgical debulking (0–10 v.s. ≥11 mm residual tumor). Solid squares represent the hazard ratio (HR) of death and open-ended horizontal lines represent the 95% confidence intervals (CIs). All P values were calculated using Cox proportional hazards analysis. (TIF 1241 kb)
ARRIVE checklist. (PDF 397 kb)
About this article
Cite this article
Sun, C., Guo, E., Zhou, B. et al. A reactive oxygen species scoring system predicts cisplatin sensitivity and prognosis in ovarian cancer patients. BMC Cancer 19, 1061 (2019). https://doi.org/10.1186/s12885-019-6288-7
- Serous ovarian Cancer
- Scoring system