Abstract
Background/Aim: To investigate the role of the expression of hypoxia-related genes on the prognosis of ovarian clear cell carcinoma (OCCC). Materials and Methods: Basal mRNA levels of eight hypoxia-related genes were compared. Cell viability was assayed after treating ES-2 cells under hypoxic conditions. The mRNA and protein levels were evaluated after the induction of hypoxia and administration of increased doses of N-acetylcysteine (NAC). Finally, the prognostic role of their expression levels was evaluated in 61 patients with OCCC. Results: The mRNA and protein levels of maspin increased gradually with the induction of hypoxia. Maspin protein expression decreased after treatment with paclitaxel and NAC. High expression of maspin was related to poor progression-free and overall survival in patients with OCCC (adjusted hazard ratios, 3.97 and 7.47; 95% confidence intervals=1.34-11.81, and 1.98-28.13). Conclusion: High expression of maspin induced by hypoxia might be associated with poor prognosis of OCCC.
Ovarian cancer is a heterogeneous group of diseases and consists of different histologic types, demonstrating distinct biological properties (1). Even though type II ovarian cancers, including high grade serous carcinoma, have valuable biomarkers such as CA-125 and human epididymis protein 4, there are few diagnostic or prognostic biomarkers for type I ovarian cancers representing low-grade serous, mucinous, endometrioid, and clear cell carcinoma (2). Thus, it is difficult to predict the clinical prognosis of type I ovarian cancers due to the lack of biomarkers (3). Since previous studies demonstrated the correlation between hypoxia and ovarian cancer progression (4, 5), further investigations for identifying new biomarkers related to hypoxia are still meaningful. Among them, ovarian clear cell carcinoma (OCCC) shows rather unique entity in terms of clinical, genetic, and histological features (6). Notably, chronic inflammation and oxidative stress in precancerous lesions such as endometriosis have been suggested to trigger its malignant transformation to OCCC (7, 8). Although hypoxia-inducible factor-1α (HIF-1α) and vascular endothelial growth factor (VEGF) are related to carcinogenesis and progression of ovarian cancer by inducing chemo-resistance (9-11), other hypoxia-related genes such as HIF-1β, cAMP-response element-binding protein-binding protein (CBP), adenovirus early region 1A-binding protein P300 (P300), factor-inhibiting hypoxia-inducible factor (FIH), von Hippel-Lindau (VHL), and maspin have not been investigated sufficiently for their roles affecting drug resistance and prognosis. There is also a lack of relevant studies about the role of hypoxia-related genes as prognostic factors in ovarian cancer. Thus, this study aimed to investigate the role of the expression of hypoxia-related genes as potential prognostic biomarkers in type I ovarian cancer, especially OCCC.
Materials and Methods
Reagents, cell lines, and tissue samples. Cisplatin and N-acetylcysteine (NAC) were acquired from Sigma-Aldrich (St. Louis, MO, USA), and paclitaxel was donated by Samyang Biopharmaceuticals Corp. (Seongnam-si, Republic of Korea). For western blot and immunohistochemistry, anti-HIF-1α was obtained from Abcam (Cambridge, UK). Anti-CBP, anti-FIH, an-ti-HIF-1β, and anti-VHL antibodies were purchased from Cell Signaling Technology (Danvers, MA, USA). Moreover, anti-P300, anti-VEGF, and anti-maspin antibodies were obtained from Santa Cruz Biotechnology (Dallas, TX, USA).
ES-2, a clear cell carcinoma cell line of human ovarian origin (12), and five frozen tissues, each from either a normal ovary or OCCC, were purchased from the Cancer Tissue Bank of Seoul National University (Seoul, Republic of Korea).
Cell culture. ES-2 cells were cultured in a 1:1 mixture of MCDB 105/Medium 199 (Sigma-Aldrich, St. Louis, MO, USA) supplemented with 15% fetal bovine serum, 100 U/ml penicillin, and 100 μg/ml streptomycin. ES-2 cells were grown in the culture media until reaching up to 70% confluency in monolayers in 100-mm tissue culture dishes; subsequently, cells were serum-starved for 24 h before the assays.
Incubation under normoxic and hypoxic conditions. For normoxic conditions, ES-2 cells were grown in a standard humidified incubator containing 5% CO2 and 21% O2 at 37°C. For hypoxic conditions, cells were incubated in a commercially available hypoxia incubator chamber (STEMCELL Technologies Inc., Vancouver, Canada) containing 90% N2; 5% CO2; 5% O2, 92% N2; 5% CO2; 3% O2, and 94% N2; 5% CO2; 1% O2.
Cytotoxicity and cell viability assays. The cultured ES-2 cells were treated with gradually escalated doses of paclitaxel, and cisplatin from 0 to 100 μmol/l was applied to cultured ES-2 cells for determining the half-maximal inhibitory concentration (IC50) under normoxic conditions. The MTT [3-(4,5-dimethylthiazol-2-yl)-2,5- diphenyltetrazolium bromide] colorimetric assay was performed for assessing cell viability. Following serum starvation for 24 h, MTT solution at 5 mg/ml was applied to ES-2 cells for 4 h at 37°C. After that, the medium was substituted with 100 μl of 0.04 N HCl in isopropanol to dissolve the formazan crystals. Cell viability was quantified with a microplate reader (Bio-rad, Seoul, Republic of Korea) at 595 nm and 655 nm as a background. The percentage of cell viability was calculated according to the following equation:
Cell viability percentage (%)=(absorbance of hypoxia or chemical compound-treated cells divided by absorbance of control cells)×100.
Quantitative real time polymerase chain reaction assay. Total cellular RNA of ES-2 cells and each frozen tissue derived from either normal ovary or OCCC were extracted by utilizing Trizol reagent (Invitrogen, Carlsbad, CA, USA) based on the manufacturer’s guidelines. Complementary DNA was synthesized using total RNA extracted from ES-2 cells with the AccuPower RT PreMix (Bioneer, Daejeon, Republic of Korea). Gene expression was measured using SYBR Green (Sigma) and a StepOnePlus qRT-PCR system (Applied Biosystems, Foster City, CA, USA). For the control and normalization of variation, the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene expression was analyzed simultaneously in each experiment. The sequences of the forward and reverse primers of the eight hypoxia-related genes are shown in Supplementary Table I. The forward primer (5’-ACA CAG AAG ACG GTG GAT GG-3’) and the reverse primer (5’-GGC AGG TCA GGT CAA CAA CA-3’) for the GAPDH gene were constructed for the amplification of GAPDH, a 193 bp product. Subsequently, expression levels of the eight hypoxia-related genes were determined using the cycle threshold (CT) values and the standard curve method. The values were normalized based on the expression of the GAPDH gene.
The conditions of polymerase chain reaction (PCR) were the following: denaturation at 95°C for 3 min, followed by 40 cycles at 95°C for 30 s, annealing at 60°C for 30 s, and DNA synthesis at 72°C for 30 s. PCR was performed using a melting curve program, which increased the temperature at a rate of 0.5°C every 10 s, ranging from 55°C to 95°C with a continuous fluorescence measurement. ROX reference dye (Invitrogen, Carlsbad, CA, USA) was used as a negative control for the fluorescence measurements. Sequence-specific products were identified by generating a melting curve, where CT value represented the cycle number at which a fluorescent signal was statistically more significant than the background. Relative gene expression levels were quantified using the 2-ΔΔCTmethod (13).
Determination of the levels of intracellular reactive oxygen species. Intracellular reactive oxygen species (ROS) production was evaluated by using 2’,7’-dichlorofluorescein diacetate (DCFH-DA, Sigma), which was converted to fluorescent 2’,7’-dichlorofluorescein (DCF) in the presence of peroxides. ES-2 cells were detached with trypsin-ethylenediaminetetraacetic acid (EDTA), collected by centrifugation, and washed with phosphate-buffered saline (PBS). The cells were treated with 10 μM DCFH-DA for 30 min at 37°C. Then, the cells were washed with PBS twice and treated with NAC in a dose-dependent manner (0, 1, 2, 5, and 10 mM) for 1 h at 37°C in a CO2 incubator, with a combination of either paclitaxel or cisplatin. The treated cells were re-washed with PBS. Fluorescent DCF intensity was analyzed using a flow cytometer (BD Bioscience, Franklin Lakes, NJ, USA).
Western blot. Protein concentration of CBP, P300, HIF-1α, HIF-1β, FIH, VHL, VEGF, and maspin in whole-cell extracts was estimated using the Bradford protein assay (Bio-Rad, Hercules, CA, USA) with bovine serum albumin as the standard. Protein extract was denatured, separated using sodium dodecyl sulphate–polyacrylamide gel electrophoresis (SDS-PAGE), transferred to a nitrocellulose membrane and developed through enhanced chemiluminescence detection (Super-Signal, West Pico, IL, USA). The protein amount was quantified by measuring light intensity from correctly sized bands under ultraviolet light using a ChemiDoc EQ system and Quantity One software (Bio-Rad). Detection of immunore-active proteins was done using goat anti-rabbit polyclonal antibodies against phospho-proteins and total-proteins at a 1:1,000 dilution and 10% SDS/PAGE gel. Total protein levels were used as a loading control for the normalization of results from protein detection by western blotting. Multiple exposures of each western blot were performed to ensure the linearity of chemiluminescent signals.
Human population. The clinicopathological data of 61 patients with OCCC between July 2001 and November 2012 were gathered. The inclusion criteria were as follows: diagnosis of OCCC; primary cytoreductive surgery followed by adjuvant chemotherapy using paclitaxel (175 mg/m2) and cisplatin (50 mg/m2); Eastern Cooperative Oncology Group (ECOG) performance status ranging from 0-2; and no underlying diseases which may affect survival.
Data on the patients’ stage based on the International Federation of Gynecology and Obstetrics (FIGO) classification, age, the extent of cytoreductive surgery, platinum-resistance, progression-free survival (PFS), and overall survival (OS) were extracted. Optimal cytoreduction was defined as cytoreductive surgery with a maximal diameter of residual tumors measured ≤1 cm. Platinum resistance was defined as a treatment-free interval less than six months following completion of adjuvant chemotherapy. PFS was defined as the period between the date of surgery and the date of disease recurrence, and OS was defined as the period from the date of surgery to the date of death or end of this study.
Immunohistochemistry. We conducted immunohistochemistry (IHC) to evaluate the expression of the eight hypoxia-related genes in patients with OCCC. For the negative control, purified non-immune mouse IgG as primary antibody substitution was utilized at the same concentration. Representative sections of core tissue (2 mm diameter) were sliced from paraffin blocks and arranged in new tissue microarray blocks using the trephine apparatus (SuperBioChips Laboratories, Seoul, Republic of Korea). Subsequently, a pathologist who was unaware of the patients’ clinicopathological characteristics interpreted and scored the eight hypoxia-related genes’ expression levels semi-quantitatively. The scoring system was based on two parameters - the product of staining intensity, read as I score (1, weak; 2, moderate; and 3, strong) and the percentage of positive cells, read as P score (1, <25%; 2, 25-50%; 3, 51-75%; and 4, >75%). The immunoreactivity score was defined as I score multiplied by P score, ranging from 1 to 12. Protein expression was considered high when the score was six or more because the median score was six.
Statistical analysis. The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board of Seoul National University Hospital (No. 1311-050-533). Informed consent was waived due to the retrospective design of the study and there were no ethical risks.
All in vitro experiments were replicated three times for obtaining consistent results, and cell viability assays were analyzed by the analysis of variance (ANOVA) method based on the general linear model (PROC-GLM) of the SAS program (SAS Institute Inc., Cary, NC, USA). For identifying the roles of the eight hypoxia-related genes in influencing platinum resistance and PFS, the Chi-square test and Student’s t-test were performed. Kaplan-Meier analysis with the log-rank test, Cox’s proportional hazard, and logistic regression analyses were done to determine the hazard ratio (HR), odds ratio (OR), and 95% confidence interval (CI). Statistical analysis was performed using SPSS software version 21.0 (SPSS Inc., Chicago, IL, USA). A p-value <0.05 was considered statistically significant.
Results
Cell viability. IC50 values of the chemotherapeutic agents paclitaxel and cisplatin were 20 and 50 μM, respectively (Supplementary Figure 1). After determining IC50 values, ES-2 cells were treated with paclitaxel (20 μM) and cisplatin (50 μM) under different hypoxic conditions. Accordingly, ES-2 cell proliferation was less inhibited under normoxic conditions, and an increase in hypoxia had a minor influence on the inhibition of cell growth (Supplementary Figure 2). These results suggest that hypoxia might be associated with chemo-resistance of OCCC cells.
mRNA expression. We compared the mRNA levels of CBP, P300, HIF-1α, HIF-1β, FIH, VHL, VEGF, and maspin between normal ovarian and OCCC tissues under normoxic conditions. The mRNA levels of VEGF were similar between normal ovarian and OCCC tissues. In contrast, the mRNA levels of CBP, P300, HIF-1α, HIF-1β, FIH, VHL, and maspin were lower in OCCC than in normal ovarian tissues (Supplementary Figure 3). Next, we compared changes in the eight hypoxia-related genes’ mRNA levels after treating ES-2 cells with paclitaxel and cisplatin under hypoxic conditions. There were no changes in the mRNA levels of CBP, HIF-1α, and VEGF, whereas the mRNA levels of HIF-1β, FIH, and VHL showed inconsistent changes after administration of paclitaxel or cisplatin under hypoxia. Moreover, only the mRNA levels of maspin increased gradually after administration of paclitaxel or cisplatin under hypoxic conditions (Figure 1).
Quantitative comparison of relative mRNA expression of the eight hypoxia-related genes including cAMP-response element-binding protein-binding protein (CBP), adenovirus early region 1A-binding protein P300 (P300), hypoxia-inducible factor-1α (HIF-1α), hypoxia inducible factor-1β (HIF-1β), factor-inhibiting hypoxia-inducible factor (FIH), von Hippel-Lindau (VHL), vascular endothelial growth factor (VEGF) and maspin after treatment with (A) paclitaxel and (B) cisplatin under hypoxic conditions in ES-2 cells (*p<0.05) (horizontal axis: induction of hypoxia, vertical axis: relative mRNA expression).
Effects of N-acetylcysteine on cell proliferation and reactive oxygen species production. To evaluate the relationship between hypoxia and ROS generation in ES-2 cells, we estimated cellular proliferation and ROS production in ES-2 cells after treatment with NAC, a potent anti-oxidant. Cell proliferation was not affected after treatment with 0 to 10 mM of NAC, whereas it was decreased after treatment with 20 mM of NAC. Moreover, ROS production was gradually reduced after treatment with NAC. Then, we evaluated ROS production in ES-2 cells treated with either paclitaxel or cisplatin in the presence of increased concentrations of NAC. The addition of 0 to 5 mM of NAC after treatment with paclitaxel resulted in a gradual decrease in ROS production. However, there was no change in ROS production after administration of cisplatin in ES-2 cells despite the presence of NAC (Figure 2).
Cell viability and comparison of reactive oxygen species (ROS) production. (A) Cell proliferation following treatment with N-acetylcysteine (NAC); ROS production (B) after treatment with NAC, (C) after treatment with paclitaxel, and (D) after treatment with cisplatin in the presence of increased concentrations of NAC.
Protein expression. Administration of paclitaxel or cisplatin under hypoxic conditions had no effect on the protein expression of CBP, P300, HIF-1α, HIF-1β, FIH, and VEGF. However, VHL protein expression of decreased after administration of cisplatin, and that of maspin increased after administration of paclitaxel under hypoxic conditions. On the other hand, maspin protein expression decreased after administration of paclitaxel in the presence of increased concentrations of NAC in ES-2 cells (Figure 3).
Comparison of the expression of cAMP-response element-binding protein-binding protein (CBP), adenovirus early region 1A-binding protein P300 (P300), hypoxia-inducible factor-1α (HIF-1α), hypoxia-inducible factor-1β (HIF-1β), factor-inhibiting hypoxia-inducible factor (FIH), von Hippel-Lindau (VHL), vascular endothelial growth factor (VEGF) and maspin proteins after treatment with paclitaxel and cisplatin (A) under hypoxic conditions and (B) in the presence of increasing concentrations of N-acetylcysteine (NAC), in ES-2 cells by western blot.
Prognosis. The clinicopathological characteristics of 61 patients diagnosed with OCCC are summarized in Table I. The median age was 49 years, and 40 patients (65.6%) had early-stage OCCC (stage I-II disease). Moreover, 53 (86.9%) and 47 patients (77%) received optimal cytoreduction and six or more cycles of chemotherapy. Nuclear expression of CBP, HIF-1β, VHL, maspin, and cytoplasmic expression of P300, HIF-1α, FIH were observed by IHC (Supplementary Figure 4). Supplementary Table ⅠI depicts factors affecting platinum resistance in univariate and multivariate analyses. Only FIGO stage III-IV disease and suboptimal cytoreduction were risk factors for platinum resistance (adjusted ORs, 27.2 and 12.9; 95%CI=2.54-29.59, and 1.83-9.04; p<0.05).
Characteristics of 61 patients with clear cell carcinoma of the ovary.
In terms of survival, high expression of maspin was related to worse PFS and OS compared to low expression of maspin (median, PFS, 91.3 vs. 118.7 months; OS, 121.4 vs. 138.9 months; p<0.05; Figure 4). In multivariate analyses, high expression of maspin was a prognostic factor for poor PFS and OS (adjusted HRs, 3.97 and 7.47; 95%CI=1.34-11.81 and 1.98-28.13; Table II).
Comparison of (A) progression-free survival and (B) overall survival between low and high expression of the eight hypoxia-related genes including cAMP-response element-binding protein-binding protein (CBP), adenovirus early region 1A-binding protein P300 (P300), hypoxia-inducible factor-1α (HIF-1α), hypoxia-inducible factor-1β (HIF-1β), factor-inhibiting hypoxia-inducible factor (FIH), von Hippel-Lindau (VHL) and maspin in patients with ovarian clear cell carcinoma.
Factors affecting survival in 61 patients with ovarian clear cell carcinoma.
Discussion
OCCC has been suggested to develop by malignant transformation of precancerous lesions such as endometriosis under chronic inflammation and oxidative stress, unlike other histologic types that may develop by müllerian differentiation or exposure to unopposed estrogen (7). Thus, prognostic values of hypoxia-related genes can differ among various types of histology under hypoxic conditions that exacerbate chronic inflammation and oxidative stress.
This study showed that most hypoxia-related genes such as CBP, P300, HIF-1α, HIF-1β, FIH, and VEGF were not over-expressed under hypoxic conditions, which means that their roles in predicting the prognosis of OCCC may be limited. These findings are supported by a previous study that showed glutathione peroxidase 3 (GPX3) over-expression in OCCC compared with other histologic types of ovarian cancer (14). Since the major biochemical function of GPX3 is to catalyze the conversion of glutathione and hydrogen peroxide into glutathione disulfide and water as one of the most important cellular detoxifying mechanisms against oxidative damages (15), GPX3 over-expression in OCCC can suppress hypoxia-related gene expression, affecting prognosis.
Maspin, a protein encoded by the SERPIN5 gene in chromosome 18q21.3-q23, is classified under the mammalian serine protease inhibitor family and functions as a tumor suppressor gene (16). Although maspin has been suggested to play a role as an anti-angiogenic inhibitor in various types of malignancies (17), its role is still controversial in ovarian cancer because the relatively few reports have demonstrated conflicting results about the relationship between maspin expression and its prognosis (18-21).
This study demonstrated that maspin expression increased after administration of paclitaxel under hypoxic conditions or in the presence of increased concentrations of NAC. These findings are supported by a previous study, which reveal that maspin expression was high in cisplatin-sensitive CAOV-3 and OVCAR-3 cells and low in cisplatin-resistant ES-2 cells (20). Moreover, high expression of stathmin, a microtubule-destabilizing protein, and maspin have been reported to play a role in the progression of adenoid cystic carcinoma (22), and thereby high expression of stathmin can decrease microtubule stabilization, which might reduce the effect of paclitaxel, in particular, in cancers with wild-type p53 such as OCCC (23).
Moreover, high expression of maspin was associated with poor prognosis of OCCC despite the no relation to platinum resistance. These results suggest that high expression of maspin induced by hypoxia could lead to resistance to paclitaxel and a subsequent worse survival in patients with OCCC. This hypothesis is supported by a previous study where the reduction in the levels of heat shock protein 27 increased chemosensitivity to paclitaxel by increasing ROS production in ovarian cancer (24).
The current study has several limitations. First, clinical outcomes, such as platinum resistance and survival, were reviewed retrospectively. Second, we used only the ES-2 cell line as an OCCC cell line. Considering that there are considerable variations regarding mutation and activity of signaling pathways, the in vitro experiments should be extended to include additional OCCC cell lines. Third, the absence of comparison to normal ovarian cells and investigation of cell viability without chemotherapy under hypoxia as a control limit the interpretation of our results. Fourth, the detailed mechanisms of hypoxia-induced expression of maspin and its relation to the resistance to paclitaxel were not evaluated. Fifth, mRNA levels were inconsistent among some genes despite the repeated experiments, which should be further investigated considering the interaction of these genes in relation to hypoxia and chemotherapy. Sixth, the prognostic value of maspin should be validated in further studies including a large number of patients.
In conclusion, this study suggests that hypoxia-induced expression of maspin might be associated with poor prognosis of OCCC. Nevertheless, this hypothesis should be investigated by further basic research using various cell lines and large-scale clinical studies.
Acknowledgements
The Authors sincerely appreciate the assistance provided by Samyang Biopharmaceuticals Corp. in donating paclitaxel for this study.
Footnotes
Authors’ Contributions
HSK designed this study. EJL, SJP, and HSK collected the data. EJL, SJP, CL, GWY, JWK and HSK analyzed the results. EJL, GWY, JWK and HSK contributed to manuscript writing. All Authors read and approved the final manuscript.
Conflicts of Interest
All Authors declare no conflicts of interest in relation to this study.
Supplementary Material
Available at: https://docs.google.com/document/d/1FAIcBjl4__tVptMy-qd6i2sRzr0moa9eBbxmq60X0Qc/edit?usp=sharing
- Received September 7, 2021.
- Revision received October 11, 2021.
- Accepted October 12, 2021.
- Copyright © 2022 International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved