Cross-validation of genes potentially associated with overall survival and drug resistance in ovarian cancer

  • Authors:
    • Yutao Gao
    • Xia Liu
    • Ting Li
    • Luwei Wei
    • Antai Yang
    • Yi Lu
    • Jian Zhang
    • Li Li
    • Sumei Wang
    • Fuqiang Yin
  • View Affiliations

  • Published online on: March 27, 2017     https://doi.org/10.3892/or.2017.5534
  • Pages: 3084-3092
Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

Ovarian cancer is the leading cause of death among malignancies of the female reproductive system. The 5-year survival rates of ovarian cancer (OC) patients are very poor as a result of recurrent disease and emergence of drug resistance; thus, studies to find predictive markers and factors for drug resistance are ongoing. In the present study, based on the microarrays from The Cancer Genome Atlas (TCGA), and the Gene Expression Omnibus (GEO) profiles covering 1648 OC patients, 11 out of 136 genes that were found to be significantly dysregulated in OC were associated with overall survival (OS) in 489 OC patients of the TCGA cohort. Of these genes, CRISP3, LYVE1, OVGP1 and BCHE were identified as independent prognostic factors, with decreased expression of the first three genes predicting shorter OS, and decreased BCHE predicting longer OS. OVGP1, BCHE and further two genes, CKAP2 and CLDN10, were consistently and remarkably associated with OS when the number of patients increased from 489 to 1583, with increased CKAP2 and decreased CLDN10 predicted shorter OS; combining the four genes provided better predictions. Associations among the four genes with OS in subgroups of OC were further verified. Downregulation of OVGP1 was significantly associated with shorter OS in all subgroups of OC patients, including subgroups of 752 patients treated with chemotherapy regimens containing taxol, 763 with both platin and taxol, 1364 with platin, 371 patients with grade 1-2 disease, 968 with grade 3 disease, 1148 with stage III-IV disease, and 439 with TP53 mutations. In addition, CKAP2 expression was significantly associated with shorter OS in 515 OC patients who had low CA125 levels. Furthermore, comprehensive analyses that including RT-qPCR, bioinformatics analysis and clinical data revealed an association of CKAP2, BCHE, CLDN10 and OVGP1 with drug resistance in OC. The genes identified in the present study might be prognostic factors as well as potential therapeutic targets in the treatment of OC.

Introduction

Ovarian cancer is the most lethal cancer among gynecologic malignancies (1). Early-stage ovarian cancer is frequently asymptomatic and difficult to detect, and thus diagnosis usually occurs after the disease has disseminated beyond the ovaries (2). The current standard treatment for advanced ovarian cancer is surgical debulking followed by platinum/taxane-based chemotherapy (3). Although this standard treatment significantly reduces the mortality rates and prolongs the survival time, the majority of patients will eventually relapse (4). The main obstacle to a successful treatment for ovarian cancer is the development of drug resistance that finally leads to fatal disease (5), and 5-year survival rates of ovarian cancer are less than 40%, with only modest improvement over the past 40 years (2). Therefore, there is a consistent and urgent need to understand the mechanism of drug resistance and to identify useful biomarkers for overall survival in ovarian cancer (6).

The present study was based on profiles of 1648 ovarian cancer patients, in microarrays from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). We identified several genes significantly associated with overall survival (OS); thus, these genes are potential biomarkers for prognosis in OC. We also found that these genes were potentially involved in regulation of drug resistance in OC.

Materials and methods

Cell culture

The human epithelial OC cell lines SKOV3 and A2780 were maintained in our laboratory and propagated in vitro by serial passage in RPMI-1640 medium supplemented with 10% fetal bovine serum (FBS). The cisplatin-resistant cell line SKOV3-DDP and the carboplatin-resistant cell line A2780-CBP were established by sequential exposure of cells to increasing concentrations of cisplatin and carboplatin, respectively (7).

Real-time quantitative polymerase chain reaction analysis

Total RNA was isolated from the cell lines SKOV3, SKOV3-DDP, A2780 and A2780-CBP, using TRIzol reagent (Life Technologies, Grand Island, NY, USA). The quantity and quality of the RNA were measured using a Thermo Scientific NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA). cDNA was synthesized from 2 µg of RNA using the SuperScript III First-Strand synthesis system (Life Technologies). mRNA expression was measured using real-time quantitative polymerase chain reaction (RT-qPCR) and the Power SYBR-Green PCR Master Mix (Applied Biosystems, Waltham, MA, USA). Data were collected with the Applied Biosystems StepOne RT-PCR system in accordance with the manufacturers instructions. The RT-qPCR gene-specific primers were: CAKP2: forward primer, 5-AA GCCACAAAACCTCAGCCT-3 and reverse primer, 5-CA TGAGGCCCTTTCCGGATT-3. BCHE: forward primer, 5-GGCCTGTCTTCAAAAGCACTG-3 and reverse primer, 5-TCCTGCTTTCCACTCCCATTC-3. CLDN10: forward primer, 5-GGAGGCTCCGATAAAGCCAA-3 and reverse primer, 5-GTGGCCCCGTTGTATGTGTA-3. OVGP1: forward primer, 5-AGCGAAGAAGCACTGGATTGA-3 and reverse primer, 5-ATTCACAGCAGATGACAGCCA-3. For GAPDH, used as the control, the forward primer was 5-GA GGTGAAGGTCGGAGT-3 and the reverse primer 5-GAA GATGGTGATGGGATTT-3.

Gene expression profiles

Gene expression data and survival information of the 1648 OC patients were deposited in the KM plotter, which was established based on the microarrays from GEO and TCGA (8). Twelve independent microarrays from GEO profiles covered 1083 patients, including GSE14764 (n=80), GSE15622 (n=36), GSE18520 (n=63), GSE19829 (n=28), GSE23554 (n=28), GSE26193 (n=107), GSE26712 (n=195), GSE27651 (n=49), GSE30161 (n=58), GSE3149 (n=166), GSE51373 (n=28) and GSE9891 (n=285); the TCGA ovarian cohort covered 565 patients. Detailed information of gene expression and clinical characteristics of the 565 OC patients in TCGA cohort were retrieved from the cBioPortal for cancer genomics (http://cbioportal.org) (9,10).

Microarray GDS3754 (Human U133 plus 2.0 GeneChips; Affymetrix, Inc., Santa Clara, CA, USA) was deposited in GEO Profiles (11), and used for the analyses of gene expression in drug-resistant OC cells. In this analysis, in the case of two probe sets that target one gene, only the probe set with significant variability was retained. In the case of more than three probe sets that target one gene, the set exhibiting the most divergent expression was excluded and the set with significant variability was retained (12).

The TCGA ovarian cohort and GEO profiles are well known publicly available cohorts that can be downloaded by researchers for further analysis. No changes were made to mRNA expression values used in the analysis.

Bioinformatics analysis

Protein/gene interaction network was generated using the GeneMANIA (http://www.genemania.org/), which is a web-based database and a tool for the prediction of gene functions on the basis of multiple networks derived from different genomic or proteomic data/sources (13). It is fast enough to predict gene functions with great accuracy using GeneMANIA because hundreds of datasets from GEO, BioGRID, Pathway Commons and I2D, as well as organismspecific functional genomics data sets have been collected in this software (14).

Biological process annotation was performed using the Coremine Medical (http://www.coremine.com/medical/), which is a gene/protein database and a web-based tool for text mining. It carries out automated extraction of experimental and theoretical biomedical knowledge from publicly available gene and text databases to create a gene-to-gene co-citation network for millions of named human genes by automated analysis of titles and abstracts in over 10 million MEDLINE records (15).

Data analysis

The data were analyzed using SPSS 20.0 software. The probability of survival and significance were calculated using the Kaplan-Meier method and KM plotter (8). The correlation between gene expression and the clinicopathological characteristics was evaluated by Pearsons χ2 test (two-sided). Gene expression was dichotomized into high and low values using the median as a cut-off in all above analyses, in accordance with previous studies (16). Gene mRNA expression levels are shown as mean ± SD. Homogeneity of the variance was analyzed using the t-test. P<0.05 was considered significant. The desired Affymetrix ID for CKAP2, BCHE, CLDN10, and OVGP1 is 218252_at, 205433_at, 205328_at, and 205432_at, respectively.

Results

Determination of genes associated with OS in OC

We previously reported that 60 genes were upregulated, and 126 genes were downregulated by at least 4-fold in 594 ovarian serous cystadenocarcinomas compared with 8 normal ovaries, respectively, according to TCGA ovarian statistics data deposited in the Oncomine database (17). Of these, we selected 136 genes (38 upregulated and 98 downregulated genes) for which we found ≤2 articles in PubMed using a search with (ovarian[Title]) and ‘gene’[Title/Abstract], for further analysis of their relationships with OS. The mRNA expression data of the 136 genes and related clinical data for 489 OC patients in the TCGA cohort were retrieved from cBioPortal for Cancer Genomics (9,10). Gene expression was dichotomized into high and low values using the median as a cut-off in all above analyses, in accordance with a previous study (16). As shown in Table I, among all 136 genes, 11 were closely associated with OS, as determined by Kaplan-Meier analysis. CKAP2 is upregulated, and ATXN10, BCHE, CLDN10, CRISP3, FCGBP, LYVE1, NDNF, OVGP1, PTGIS and REEP1 are downregulated in OC (17). We found that the upregulation of CKAP2, and downregulation of CLDN10, CRISP3 and OVGP1 were associated with shorter OS, whereas the downregulation of ATXN10, BCHE, FCGBP, LYVE1, NDNF, PTGIS and REEP1 were associated with longer OS. To elucidate whether any of the above genes was an independent factor for predicting OS, we performed multivariate analyses of histological grade, tumor stage, residual tumor and platinum status for the 11 genes, using a Cox proportional hazards model (Table II). The results indicated that BCHE (P=0.026), CRISP3 (P=0.031), LYVE1 (P=0.014) and OVGP1 (P=0.001) were independent prognostic factors for OS in 489 OC patients.

Table I.

Eleven genes were determined by Kaplan-Meier analysis to be closely associated with overall survival in ovarian cancer.

Table I.

Eleven genes were determined by Kaplan-Meier analysis to be closely associated with overall survival in ovarian cancer.

95% CI

Dysregulation of genes EstimateSth. ErrorLowerUpperP-value
Upregulated genes
  CKAP2H39.9002.48235.03544.7650.024
L47.5003.06041.50353.497
Overall43.8002.10539.67547.925
Downregulated genes
  ATXN10H41.0002.78235.54746.4530.03
L45.3002.67140.06650.534
Overall43.8002.10539.67547.925
  BCHEH36.3002.28331.82640.7740.000
L49.0003.20642.71755.283
Overall43.8002.10539.67547.925
  CLDN10H48.3001.79744.77751.8230.036
L38.2002.07534.13442.266
Overall43.8002.10539.67547.925
  CRISP3H47.7001.79544.18151.2190.043
L38.3002.02334.33642.264
Overall43.8002.10539.67547.925
  FCGBPH41.0002.37136.35245.6480.025
L44.9002.86639.28350.517
Overall43.8002.10539.67547.925
  LYVE1H41.5002.22137.14845.8520.032
L44.9003.64337.76152.039
Overall43.8002.10539.67547.925
  NDNFH41.0002.07036.94445.0560.034
L45.1002.08641.01249.188
Overall43.8002.10539.67547.925
  OVGP1H49.0003.58141.98156.0190.003
L37.9001.94734.08441.716
Overall43.8002.10539.67547.925
  PTGISH40.4002.89034.73546.0650.004
L47.5001.97543.62951.371
Overall43.8002.10539.67547.925
  REEP1H38.4002.05334.37742.4230.007
L47.7003.02141.77953.621
Overall43.8002.10539.67547.925

[i] H, high expression; L, low expression. Gene expression and survival data of 489 ovarian cancer patients in a TCGA cohort were used for the analysis. Expression values of a gene were dichotomized into high and low expression using the median as a cut-off.

Table II.

Multivariate analysis of overall survival of ovarian cancer patients using Cox proportional hazard model.

Table II.

Multivariate analysis of overall survival of ovarian cancer patients using Cox proportional hazard model.

95% CI

FactorsBSEWalddfSig.Exp(B)LowerUpper
CKAP20.1720.1351.62610.2021.1870.9121.546
ATXN100.1560.1261.53810.2151.1690.9131.496
BCHE0.3100.1394.95510.0261.3631.0381.790
CLDN10−0.1840.1301.98810.1590.8320.6441.074
CRISP3−0.2820.1314.64710.0310.7540.5840.975
FCGBP0.1960.1282.33010.1271.2170.9461.565
LYVE10.3240.1326.00410.0141.3821.0671.791
NDNF0.0570.1320.18410.6681.0580.8181.370
OVGP1−0.4390.13510.57710.0010.6450.4950.840
PTGIS0.0300.1420.04510.8321.0300.7811.360
REEP10.0750.1320.32510.5681.0780.8331.395
Grade
  (G2/G3)−0.2790.1972.01010.1560.7560.5141.113
Platinum status
  (R/S)1.2300.16853.69010.0003.4202.4624.753
Residual (mm)
  (≤10/>10)−0.1760.1451.48210.2230.8380.6311.113
Stage
  (II/III–IV)−0.6360.4062.44910.1180.5300.2391.174

[i] R/S, platinum resistance/platinum sensitive. Gene expression and survival data of 489 ovarian cancer patients in a TCGA cohort were used for the analysis. Gene expression was dichotomized into high and low values using the median as a cut-off.

The associations of the 11 genes with OS were further validated using KM plotter which included gene expression data and survival information of 1287 OC patients downloaded from GEO and TCGA ovarian data in 2012 (8), and updated to 1648 patients for the 2015 version. As shown in Fig. 1, among the 11 genes, CKAP2, BCHE, CLDN10 and OVGP1 were consistently associated with OS when the number of patients increased from 489 in TCGA cohort to 1583 in KM plotter (Table I and Fig. 1). Of these, upregulation of CKAP2, and downregulation of CLDN10 and OVGP1 were associated with shorter OS, whereas downregulation of BCHE was associated with longer OS. When the four genes were combined, we found that their mean expression in combination was greatly more associated with OS than any of the genes used separately (Fig. 2).

Verification of genes associated with OS in OC patients treated with different chemotherapy regimens

As shown in Table III, in subgroups of patients treated with different regimens, we found that the high expression of CKAP2 was associated with shorter OS in 752 patients treated with taxol-based regimens, and in 763 treated with regimens containing both platin and taxol, but the association was not significant in 1364 patient treated with only platins; downregulation of BCHE predicts a longer OS in patients treated with platin, but not in patients treated with taxols or taxols with platin; downregulation of CLDN10 was associated with shorter OS in patients treated with platin or platin with taxol, but not in the taxol-only group, and downregulation of OVGP1 predicted a shorter OS in all three subgroups of ovarian patients treated with different regimens.

Table III.

Four genes were associated with overall survival in subgroups of ovarian patients treated with different chemotherapy regiments, as determined by KM plotter.

Table III.

Four genes were associated with overall survival in subgroups of ovarian patients treated with different chemotherapy regiments, as determined by KM plotter.

Chemotherapy containing (no. of patients)

Platin (1364)Taxol (780)Platin and taxol (763)



GeneHR (95% CI)P-valueHR (95% CI)P-valueHR (95% CI)P-value
CKAP21.09 (0.94–1.26)0.271.24 (1.02–1.5)0.0331.24 (1.01–1.51)0.037
BCHE1.26 (1.09–1.46)0.00171.19 (0.97–1.44)0.0881.19 (0.97–1.45)0.091
CLDN100.76 (0.66–0.88)0.000290.82 (0.68–1)0.0520.81 (0.67–0.99)0.042
OVGP10.77 (0.66–0.89)0.000370.74 (0.61–0.91)0.00310.74 (0.61–0.91)0.0035

[i] Gene expression was dichotomized into high and low values using the median as a cut-off. HR, hazard ratio.

Verification of genes associated with OS in subgroups of OC patients

The relationships between dysregulation of the four genes (CKAP2, BCHE, CLDN10 and OVGP1) with OS in subgroups of OC patients with different CA125 levels, grades, stages and TP53 mutation status were further investigated. We revealed that the high CKAP2 expression significantly associated with shorter OS in 515 OC patients having low CA125 levels (the average CA125 below lower quartile) (Fig. 3). Besides, as shown in Table IV, all four genes were significantly associated with OS in 968 OC patients with grade 3 disease, but only BCHE and OVGP1 were associated with OS in 371 patients with grade 1–2 disease. CLDN10 and OVGP1 were associated with OS in 1148 patients with stage III–IV OC, but none of the genes was associated with OS in 133 patients with stage I–II OC. In the subgroup of 439 patients with TP53 mutation, downregulation of OVGP1 predicted shorter OS, whereas downregulation of CLDN10 predicted shorter OS in the TP53 wild-type group. In addition, downregulation of BCHE was associated with longer PFS and PPS in cohorts of 1307 and 708 patients, respectively; whereas downregulation of CLDN10 and OVGP1 were associated with shorter PFS and PPS (Fig. 4).

Table IV.

The four genes were associated with overall survival in subgroups of ovarian patients with different grades, stages and TP53 mutation status, as determined by KM plotter.

Table IV.

The four genes were associated with overall survival in subgroups of ovarian patients with different grades, stages and TP53 mutation status, as determined by KM plotter.

CKAP2BCHECLDN10OVGP1




SubgroupNo. of patientsHR (95% CI)P-valueHR (95% CI)P-valueHR (95%CI)P-valueHR (95%CI)P-value
Grade
  I–II3711.22 (0.91–1.64)0.181.73 (1.28–2.34)0.000260.78 (0.58–1.04)0.0910.58 (0.43–0.78)0.00027
  III9681.21 (1.02–1.43)0.0321.20 (1.01–1.42)0.040.75 (0.63–0.89)0.000880.80 (0.67–0.95)0.011
Stage
  I–II1331.36 (0.61–3.01)0.451.04 (0.47–2.31)0.920.90 (0.41–1.97)0.780.70 (0.32–1.55)0.38
  III–IV11481.1 (0.94–1.29)0.221.13 (0.97–1.32)0.120.84 (0.72–0.98)0.0240.73 (0.62–0.85)8.0e-5
TP53
  Wild860.99 (0.55–1.78)0.972.0 (1.09–3.66)0.0230.78 (0.43–1.39)0.40.8 (0.44–1.43)0.45
  Mutation4391.02 (0.79–1.31)0.91.12 (0.87–1.44)0.370.85 (0.66–1.09)0.20.75 (0.58–0.96)0.024

[i] Gene expression was dichotomized into high and low values using the median as a cut-off. HR, hazard ratio.

Clinical importance of the analysis

The four OS-related genes, CKAP2, BCHE, CLDN10 and OVGP1, with clinical characteristics such as tumor stage, histological grade, and primary therapy outcome were checked in the TCGA cohort (n=489). As shown in Table V, BCHE and OVGP1 were closely associated with primary therapy outcomes. Downregulation of BCHE was associated with good complete response, whereas downregulation of OVGP1 was associated with poor complete response. Besides, downregulation of BCHE was associated with lower tumor stage, in particular, stage II. Only CLDN10 was associated with platinum status; its downregulation might predict drug resistance. The 4 genes were not significantly associated with histological grade in this analysis.

Table V.

Correlation between gene expression and the clinicopathological characteristics of 489 patients with ovarian cancer in TCGA cohort.

Table V.

Correlation between gene expression and the clinicopathological characteristics of 489 patients with ovarian cancer in TCGA cohort.

CKAP2BCHECLDN10OVGP1




CharacteristicsNo. of patients (%)High n (%)Low n (%)High P-valueLow n (%)High n (%)P-valueHigh n (%)Low n (%)P-valueHigh n (%)Low n (%)P-value
Primary therapy outcome395 (100) 0.451 0.005 0.521 0.012
Complete response276 (69.9)142 (51.4)134 (48.6) 127 (46.0)149 (54.0) 149 (54.0)127 (46.0) 148 (53.6)128 (46.4)
Partial response57 (14.4)23 (40.4)34 (59.6) 37 (64.9)20 (35.1) 25 (43.9)32 (56.1) 20 (35.1)37 (64.9)
Progressive disease37 (9.4)20 (54.1)17 (45.9) 26 (70.3)11 (29.7) 18 (48.6)19 (51.4) 12 (32.4)25 (67.6)
Stable disease25 (6.3)12 (48.0)13 (52.0) 12 (48.0)13 (52.0) 12 (48.0)13 (52.0) 11 (44.0)14 (56.0)
Tumor stage484 (100.0) 0.110 0.013 0.094 0.602
  II24 (5.0)13 (54.2)11 (45.8) 5 (20.8)19 (79.2) 17 (70.8)7 (29.2) 10 (41.7)14 (58.3)
  III381 (78.7)198 (52.0)183 (48.0) 197 (51.7)184 (48.3) 184 (48.3)197 (51.7) 193 (50.7)188 (49.3)
  IV79 (16.3)31 (39.2)48 (60.8) 38 (48.1)41 (51.9) 41 (51.9)38 (48.1) 37 (46.8)42 (53.2)
Platinum status287 (100) 0.377 0.079 0.017 0.153
  Resistance90 (31.4)45 (50.0)45 (50.0) 49 (54.4)41 (45.6) 36 (40.0)54 (60.0) 50 (55.6)40 (44.4)
  Sensitive197 (69.6)93 (47.2)104 (52.8) 88 (44.7)109 (55.3) 107 (54.3)90 (45.7) 95 (48.2)102 (51.8)

[i] P-value was evaluated using Pearsons χ2 test (2-sided).

Possible associations with drug resistance

The four genes, CKAP2, BCHE, CLDN10 and OVGP1, were dysregulated in platinum-resistant OC cells, as determined using microarray data retrieved from GEO Profiles and RT-qPCR analysis (Fig. 5). Compared with their expression in sensitive cells, expression of CKAP2 was increased in cisplatin-resistant cells, but decreased in carboplatin-resistant cells; BCHE, CLDN10 and OVGP1 all were decreased in cisplatin- and carboplatin-resistant cells, although the changes in CLDN10 and OVGP1 were not significant in one analysis (Fig. 5A).

Bioinformatics analyses were performed to further explain the associations of the four genes with drug resistance in OC. As shown in Fig. 6A, text mining indicated that BCHE had direct associations with drug resistance and OC, and all the four genes were indirectly associated with drug resistance and OC through six biological processes (P<0.05), including cell proliferation, apoptosis, mitotic nuclear division, phosphorylation, growth and gene expression. We further generated a protein/gene interaction network for the four genes and their products with 49 drug-resistance proteins/genes in OC, which including 15 tumor suppressors (19): BRCA1, BRCA2, CHEK2, FBXO32, MLH1, SULF1, IL24, CDKN2A, CDKN1A, TP53, TP73, PDCD4, PTEN, RASSF1 and WWOX, 25 oncoproteins (20): ACTN4, AKT1, AKT2, BAX, BCL2, MIEN1, CLU, CUZD1, DAXX, EGFR, ERBB2, FOS, JUN, IKBKE, KRAS, MET, MYC, NFKB1, NINL, NOTCH3, PIK3CA, RSF1, SRC, STAT3 and URI1, NEK2 (21), NEK11 (22), TRPC1 (12), CCL21 and SPARCL1 (23), GGNBP2 and RNASET2 (7), ALDH1A2 and ADH1B (17). As shown in Fig. 6B, CKAP2, BCHE, CLDN10 and OVGP1 had direct and indirect interactions with all of the 49 drug-resistant proteins/genes, and directly interacted with 11, 23, 19 and 13 of those 49 proteins/genes, respectively. The four genes and their products also directly interacted with each other, with OVGP1 and BCHE both directly interacting with CKAP2 and CLDN10.

Discussion

Biomarkers not only have prognostic implications, but are also helpful in monitoring treatment responses, surveillance for tumor recurrence and guidance of clinical decisions (24). Long-term survival of OC patients remains very poor as a result of recurrence and emergence of drug resistance, and the 5-year survival rates is only ~40% (2). Thus, prognostic biomarkers for OC patients are particularly necessary and crucial, and there is ongoing search for predictive biomarkers. In the present study, we identified a group of genes significantly associated with OS in 489 OC patients from a TCGA cohort (Table I), of which BCHE, CRISP3, LYVE1 and OVGP1 were identified as independent risk prognostic factors for OS (Table II). Four genes (CKAP2, BCHE, CLDN10 and OVGP1) were consistently associated with OS when the number of patients increased from 489 to 1583 (Fig. 1). Of these, upregulation of CKAP2, and downregulation of CLDN10 and OVGP1 were associated with shorter OS, whereas downregulation of BCHE was associated with longer OS. With the exception of CKAP2, the other three genes all were significantly related to PFS and PPS, with decreased BCHE expression associated with longer PFS and PPS in 1307 and 708 OC patients, respectively; and decreased expression of CLDN10 and OVGP1 were associated with shorter PFS and PPS (Fig. 4). In subgroups of patients with OC, we further verified that downregulation of OVGP1 was significantly associated with shorter OS in all OC subgroups, including the 752 patients treated with chemotherapy regimens that contained taxol, 763 treated with both platin and taxol, 1364 treated with platin, 371 patients with grade 1–2 disease, 968 with grade 3 disease, 1148 with stage III–IV disease, and 439 with TP53 mutations (Tables III and IV). All these results together suggest that CKAP2, BCHE, CLDN10 and OVGP1, (especially OVGP1), are potential biomarkers for predicting OS in OC.

The relationship between these genes and prognosis in OC is poorly studied. Upregulation of chromatin CKAP2 is an independent prognostic marker for shorter relapse-free survival in early-stage breast cancer (25), and for significantly worse prognoses in male patients with T1-2 gastric cancer (26). Downregulation of CLDN10 indicates shorter OS in some patients with lung adenocarcinoma (27). All these reports were consistent with our findings in this study. However, low levels of preoperative serum BCHE were validated as an independent negative prognostic factor for prostate cancer patients who undergo radical prostatectomy (28), and a predictor of shorter OS and DFS in patients with muscle-invasive bladder cancer who undergo radical cystectomy (29). These results were contrary to our findings that downregulation of BCHE was associated with longer OS. The associations of the four genes with OS in OC is rare, with only one study reporting that the downregulation of BCHE predicts a longer OS in OC patients, which was also based on the analysis of data in the TCGA cohort (30).

Serum CA125 is widely used to distinguish malignant from benign pelvic masses, to monitor response to treatment, to assess prognosis and to detect disease recurrence. Low levels of CA125 are normally an optimistic indicator of OC progression and prognosis. For example, preoperative CA125 <65 U/ml correlates with longer survival in OC patients (31); CA125 <50 U/ml correlates with significantly longer DFS and OS in patient with borderline ovarian tumors (32). However, no data are available that further predicts OS in subgroups of patients with low levels of CA125. CKAP2 expression was significantly associated with shorter OS in 515 OC patients who had low CA125 levels (the average CA125 below lower quartile) (Fig. 3), which suggests that this gene could be used for prognosis prediction with CA125. However, this study measured tissue levels of CA125, whereas the FDA-approved test for OC is serum-based (8). Associations between CKAP2 and OS in patients with low CA125 levels should be further studied.

The relationships of the four genes with drug resistance in cancer has not been widely reported. In the present study, we found that CKAP2, BCHE, CLDN10 and OVGP1 all were significantly dysregulated in cisplatin- and carboplatin-resistant cells (Fig. 5). Text mining indicated that the four genes were directly and indirectly associated with drug resistance through several biological processes. A protein/gene interaction network indicated that the four genes interacted with 49 drug-resistant proteins/genes in OC. In particular, BCHE and CLDN10 directly interacted with 23 and 19 of those 49 proteins/genes (Fig. 6). Downregulation of CLDN10 was also an indicator of drug resistance, and downregulation of OVGP1 was significantly associated with poor complete response (Table IV). All these results imply that dysregulation of the four genes affects drug resistance in OC. Furthermore, among patients treated with different regimens, high expression of CKAP2 was associated with shorter OS in patients treated with taxol-containing regimens, whereas downregulation of BCHE was associated with OS in patient treated by regimens without taxol; downregulation of CLDN10 was associated with shorter OS in patients treated with platin, and downregulation of OVGP1 predicted a shorter OS in patients treated with either platin or taxol (Table III). These results indicated that CKAP2 is associated with taxol resistance, BCHE and CLDN10 with platin resistance, and OVGP1 with both platin and taxol resistance.

The roles of CKAP2, BCHE, CLDN10 and OVGP1 in OC have rarely been studied. We previously reported that CKAP2 was upregulated, and BCHE, CLDN10 and OVGP1 were downregulated in OC, by a 4.26, −8.39, −7.04 and −102.95-fold changes, respectively (17). Downregulation of OVGP1 in OC was recently confirmed in a study that showed the protein to be less abundant in high-grade serous ovarian tumor fluids (malignant) than in benign serous cystadenoma tumor fluids (33). CLDN10 mRNA was specifically detected in five cancerous ovaries compared with three normal controls in experimental hens, which suggests that CLDN10 is a novel biomarker for detecting OC in chickens (34); however, its association with human OC was not discussed.

In summary, we identified 11 genes significantly related to OS in 489 OC patients. Of these, BCHE, CRISP3, LYVE1 and OVGP1 were identified as independent risk prognostic factors. Four genes CKAP2, BCHE, CLDN10, and particularly the OVGP1, were remarkably associated with OS in total of 1583 OC patients, and combination of these genes had much better prediction potential. Besides, comprehensive analysis indicated that CKAP2, BCHE, CLDN10 and OVGP1 might contribute to drug resistance in OC. This study has implicated genes that might be both prognostic markers and potential therapeutic targets to pursue in the treatment of ovarian cancer.

Acknowledgements

This study was supported by the National Natural Science Foundation of China (grant nos. 81302283, 81560424, 81660606 and 81460397), the China Postdoctoral Science Foundation (nos. 2014M552535XB and 2014M552291), the Natural Science Foundation of Guangxi (nos. 2014GXNSFCA118010, 2015GXNSFBA139115, 2015GXNSFAA139151 and 2014GXNSFBA118155).

References

1 

Coleman RL, Monk BJ, Sood AK and Herzog TJ: Latest research and treatment of advanced-stage epithelial ovarian cancer. Nat Rev Clin Oncol. 10:211–224. 2013. View Article : Google Scholar : PubMed/NCBI

2 

Siegel RL, Miller KD and Jemal A: Cancer statistics, 2016. CA Cancer J Clin. 66:7–30. 2016. View Article : Google Scholar : PubMed/NCBI

3 

Vaughan S, Coward JI, Bast RC Jr, Berchuck A, Berek JS, Brenton JD, Coukos G, Crum CC, Drapkin R, Etemadmoghadam D, et al: Rethinking ovarian cancer: Recommendations for improving outcomes. Nat Rev Cancer. 11:719–725. 2011. View Article : Google Scholar : PubMed/NCBI

4 

Pliarchopoulou K and Pectasides D: Epithelial ovarian cancer: Focus on targeted therapy. Crit Rev Oncol Hematol. 79:17–23. 2011. View Article : Google Scholar : PubMed/NCBI

5 

Cannistra SA: Cancer of the ovary. N Engl J Med. 351:2519–2529. 2004. View Article : Google Scholar : PubMed/NCBI

6 

Holohan C, Van Schaeybroeck S, Longley DB and Johnston PG: Cancer drug resistance: An evolving paradigm. Nat Rev Cancer. 13:714–726. 2013. View Article : Google Scholar : PubMed/NCBI

7 

Yin F, Liu L, Liu X, Li G, Zheng L, Li D, Wang Q, Zhang W and Li L: Downregulation of tumor suppressor gene ribonuclease T2 and gametogenetin binding protein 2 is associated with drug resistance in ovarian cancer. Oncol Rep. 32:362–372. 2014.PubMed/NCBI

8 

Gyorffy B, Lánczky A and Szállási Z: Implementing an online tool for genome-wide validation of survival-associated biomarkers in ovarian-cancer using microarray data from 1287 patients. Endocr Relat Cancer. 19:197–208. 2012. View Article : Google Scholar : PubMed/NCBI

9 

Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, et al: The cBio cancer genomics portal: An open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2:401–404. 2012. View Article : Google Scholar : PubMed/NCBI

10 

Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, et al: Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 6:pl12013. View Article : Google Scholar : PubMed/NCBI

11 

Edgar R, Domrachev M and Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 30:207–210. 2002. View Article : Google Scholar : PubMed/NCBI

12 

Liu X, Zou J, Su J, Lu Y, Zhang J, Li L and Yin F: Downregulation of transient receptor potential cation channel, subfamily C, member 1 contributes to drug resistance and high histological grade in ovarian cancer. Int J Oncol. 48:243–252. 2016.PubMed/NCBI

13 

Mostafavi S, Ray D, Warde-Farley D, Grouios C and Morris Q: GeneMANIA: A real-time multiple association network integration algorithm for predicting gene function. Genome Biol. 9:(Suppl 1). S42008. View Article : Google Scholar : PubMed/NCBI

14 

Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT, et al: The GeneMANIA prediction server: Biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 38:(Web Server). W214–W220. 2010. View Article : Google Scholar : PubMed/NCBI

15 

Jenssen TK, Laegreid A, Komorowski J and Hovig E: A literature network of human genes for high-throughput analysis of gene expression. Nat Genet. 28:21–28. 2001. View Article : Google Scholar : PubMed/NCBI

16 

Hedditch EL, Gao B, Russell AJ, Lu Y, Emmanuel C, Beesley J, Johnatty SE, Chen X, Harnett P, George J, et al: Australian Ovarian Cancer Study Group: ABCA transporter gene expression and poor outcome in epithelial ovarian cancer. J Natl Cancer Inst. 106:1062014. View Article : Google Scholar

17 

Liu X, Gao Y, Zhao B, Li X, Lu Y, Zhang J, Li D, Li L and Yin F: Discovery of microarray-identified genes associated with ovarian cancer progression. Int J Oncol. 46:2467–2478. 2015.PubMed/NCBI

18 

Li M, Balch C, Montgomery JS, Jeong M, Chung JH, Yan P, Huang TH, Kim S and Nephew KP: Integrated analysis of DNA methylation and gene expression reveals specific signaling pathways associated with platinum resistance in ovarian cancer. BMC Med Genomics. 2:342009. View Article : Google Scholar : PubMed/NCBI

19 

Yin F, Liu X, Li D, Wang Q, Zhang W and Li L: Tumor suppressor genes associated with drug resistance in ovarian cancer (Review). Oncol Rep. 30:3–10. 2013.PubMed/NCBI

20 

Liu X, Gao Y, Lu Y, Zhang J, Li L and Yin F: Oncogenes associated with drug resistance in ovarian cancer. J Cancer Res Clin Oncol. 141:381–395. 2015. View Article : Google Scholar : PubMed/NCBI

21 

Liu X, Gao Y, Lu Y, Zhang J, Li L and Yin F: Upregulation of NEK2 is associated with drug resistance in ovarian cancer. Oncol Rep. 31:745–754. 2014.PubMed/NCBI

22 

Liu X, Gao Y, Lu Y, Zhang J, Li L and Yin F: Downregulation of NEK11 is associated with drug resistance in ovarian cancer. Int J Oncol. 45:1266–1274. 2014.PubMed/NCBI

23 

Yin F, Liu X, Li D, Wang Q, Zhang W and Li L: Bioinformatic analysis of chemokine (C-C motif) ligand 21 and SPARC-like protein 1 revealing their associations with drug resistance in ovarian cancer. Int J Oncol. 42:1305–1316. 2013.PubMed/NCBI

24 

Wong KF, Xu Z, Chen J, Lee NP and Luk JM: Circulating markers for prognosis of hepatocellular carcinoma. Expert Opin Med Diagn. 7:319–329. 2013. View Article : Google Scholar : PubMed/NCBI

25 

Kim HS, Koh JS, Choi YB, Ro J, Kim HK, Kim MK, Nam BH, Kim KT, Chandra V, Seol HS, et al: Chromatin CKAP2, a new proliferation marker, as independent prognostic indicator in breast cancer. PLoS One. 9:e981602014. View Article : Google Scholar : PubMed/NCBI

26 

Kim YW, Eom BW, Kook MC, Kim HS, Kim MK, Hwang HL, Chandra V, Poojan S, Song Y, Koh JS, et al: Clinical implications of proliferation activity in T1 or T2 male gastric cancer patients. Exp Mol Med. 47:e1932015. View Article : Google Scholar : PubMed/NCBI

27 

Zhang Z, Wang A, Sun B, Zhan Z, Chen K and Wang C: Expression of CLDN1 and CLDN10 in lung adenocarcinoma in situ and invasive lepidic predominant adenocarcinoma. J Cardiothorac Surg. 8:952013. View Article : Google Scholar : PubMed/NCBI

28 

Koie T, Ohyama C, Hatakeyama S, Imai A, Yoneyama T, Hashimoto Y, Yoneyama T, Tobisawa Y, Hosogoe S, Yamamoto H, et al: Significance of preoperative butyrylcholinesterase as an independent predictor of biochemical recurrence-free survival in patients with prostate cancer treated with radical prostatectomy. Int J Clin Oncol. 21:379–383. 2016. View Article : Google Scholar : PubMed/NCBI

29 

Koie T, Ohyama C, Yamamoto H, Hatakeyama S, Imai A, Yoneyama T, Hashimoto Y, Kitayam M and Hirota K: Significance of preoperative butyrylcholinesterase as an independent predictor of survival in patients with muscle-invasive bladder cancer treated with radical cystectomy. Urol Oncol. 32:820–825. 2014. View Article : Google Scholar : PubMed/NCBI

30 

Willis S, Villalobos VM, Gevaert O, Abramovitz M, Williams C, Sikic BI and Leyland-Jones B: Single gene prognostic biomarkers in ovarian cancer: A meta-analysis. PLoS One. 11:e01491832016. View Article : Google Scholar : PubMed/NCBI

31 

Nagele F, Petru E, Medl M, Kainz C, Graf AH and Sevelda P: Preoperative CA 125: An independent prognostic factor in patients with stage I epithelial ovarian cancer. Obstet Gynecol. 86:259–264. 1995. View Article : Google Scholar : PubMed/NCBI

32 

Tang A, Kondalsamy-Chennakesavan S, Ngan H, Zusterzeel P, Quinn M, Carter J, Leung Y and Obermair A: Prognostic value of elevated preoperative serum CA125 in ovarian tumors of low malignant potential: A multinational collaborative study (ANZGOG0801). Gynecol Oncol. 126:36–40. 2012. View Article : Google Scholar : PubMed/NCBI

33 

Poersch A, Grassi ML, Carvalho VP, Lanfredi GP, Palma CS, Greene LJ, de Sousa CB, Carrara HH, Dos Candido Reis FJ and Faça VM: A proteomic signature of ovarian cancer tumor fluid identified by highthroughput and verified by targeted proteomics. J Proteomics. 145:226–236. 2016. View Article : Google Scholar : PubMed/NCBI

34 

Seo HW, Rengaraj D, Choi JW, Ahn SE, Song YS, Song G and Han JY: Claudin 10 is a glandular epithelial marker in the chicken model as human epithelial ovarian cancer. Int J Gynecol Cancer. 20:1465–1473. 2010.PubMed/NCBI

Related Articles

Journal Cover

May-2017
Volume 37 Issue 5

Print ISSN: 1021-335X
Online ISSN:1791-2431

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Gao Y, Liu X, Li T, Wei L, Yang A, Lu Y, Zhang J, Li L, Wang S, Yin F, Yin F, et al: Cross-validation of genes potentially associated with overall survival and drug resistance in ovarian cancer. Oncol Rep 37: 3084-3092, 2017
APA
Gao, Y., Liu, X., Li, T., Wei, L., Yang, A., Lu, Y. ... Yin, F. (2017). Cross-validation of genes potentially associated with overall survival and drug resistance in ovarian cancer. Oncology Reports, 37, 3084-3092. https://doi.org/10.3892/or.2017.5534
MLA
Gao, Y., Liu, X., Li, T., Wei, L., Yang, A., Lu, Y., Zhang, J., Li, L., Wang, S., Yin, F."Cross-validation of genes potentially associated with overall survival and drug resistance in ovarian cancer". Oncology Reports 37.5 (2017): 3084-3092.
Chicago
Gao, Y., Liu, X., Li, T., Wei, L., Yang, A., Lu, Y., Zhang, J., Li, L., Wang, S., Yin, F."Cross-validation of genes potentially associated with overall survival and drug resistance in ovarian cancer". Oncology Reports 37, no. 5 (2017): 3084-3092. https://doi.org/10.3892/or.2017.5534