Abstract
Purpose: We aimed to find key molecules associated with chemoresistance in ovarian cancer using gene expression profiling as a screening tool.
Experimental Design: Using two newly established paclitaxel-resistant ovarian cancer cell lines from an original paclitaxel-sensitive cell line and four supersensitive and four refractory surgical ovarian cancer specimens from paclitaxel-based chemotherapy, molecules associated with chemoresistance were screened with gene expression profiling arrays containing 39,000 genes. We further analyzed 44 genes that showed significantly different expressions between paclitaxel-sensitive samples and paclitaxel-resistant samples with permutation tests, which were common in cell lines and patients' tumors.
Results: Eight of these genes showed reproducible results with real-time reverse transcription-PCR, of which indoleamine 2,3-dioxygenase gene expression was the most prominent and consistent. Moreover, by immunohistochemical analysis using a total of 24 serous-type ovarian cancer surgical specimens (stage III, n = 21; stage IV, n = 7), excluding samples used for GeneChip analysis, the Kaplan-Meier survival curve showed a clear relationship between indoleamine 2,3-dioxygenase staining patterns and overall survival (log-rank test, P = 0.0001). All patients classified as negative survived without relapse. The 50% survival of patients classified as sporadic, focal, and diffuse was 41, 17, and 11 months, respectively.
Conclusion: The indoleamine 2,3-dioxygenase screened with the GeneChip was positively associated with paclitaxel resistance and with impaired survival in patients with serous-type ovarian cancer.
Ovarian cancer is one of the primary causes of death related to gynecologic malignancies (1). Nearly 65% of ovarian cancer patients die from their disease within 5 years (2). Although ovarian cancer is considered highly responsive to combination therapy with paclitaxel and carboplatin (3), cancer recurs rapidly in >50% of responsive patients, and in many cases, the recurring cancer cells develop chemoresistance (4). Therefore, countering chemoresistance is essential for ovarian cancer management.
Properties within tumor cells that may lead to drug resistance in ovarian cancer include multidrug resistance proteins and mismatched repair processes (e.g., alterations in the p53 pathway; refs. 5–7). In addition, various molecules have been documented as candidates for chemoresistance in ovarian cancer (8–12). However, molecular targeting to overcome chemoresistance has not yet been delineated in ovarian cancer.
The development of microarray methods for large-scale analysis of gene expression makes it possible to search systematically for key molecules that may be involved in chemoresistance (13). We have already applied this approach to ovarian cancer (14) as well as to other cancers (15, 16). In previous works on ovarian cancer, gene expression profiling was used to distinguish types of ovarian cancer (17), malignant transformation from normal tissue (18, 19), serous uterine from ovarian cancers (20), or metastatic from nonmetastatic disease (21). Although some advances have been seen in chemoresistance of childhood acute lymphoblastic leukemia as well as other types of cancers (22–24), the technology has not elucidated a set of genes associated with chemoresistance, a critical factor for improving prognosis in most cancers.
In this experiment, GeneChip was applied to screen molecules expressed differentially between chemoresistant and chemosensitive cell lines as well as cancer cells derived from patients who were either clinically sensitive or resistant to chemotherapy. The clinical significance of a prominent molecule was further confirmed with immunohistochemical analysis to predict recurrence after chemotherapy.
Materials and Methods
Tumor specimens. The Jikei University School of Medicine Ethics Review Committee approved the study protocol with informed consent from all patients. A total of 32 ovarian cancer surgical specimens were obtained at the Jikei University Hospitals. Tumors were histologically classified according to the WHO international system and staged according to the International Federation of Gynecology and Obstetrics (25). All of the 32 cases underwent debulking surgery, and the sizes of the residual tumors were <2 cm in all cases. All cases were serous cystadenocarcinomas. There were 25 stage III cases and 7 stage IV cases.
Among the 32 cases, 4 patients with stage IIIc were diagnosed as having achieved a pathologic complete response according to a second look operation after six courses of chemotherapy, including paclitaxel; cancer did not recur in these patients for >1 year. These cases were termed “supersensitive.” In addition, we also selected four patients with stage IIIc who showed progressive disease during chemotherapy, including paclitaxel; these cases were termed “refractory.” Three of four supersensitive cases completed six courses of paclitaxel (180 mg/m2)-carboplatin (AUC 5), and one supersensitive case underwent six courses of paclitaxel only due to an allergic reaction to carboplatin in the first course. On the contrary, refractory cases underwent two to four courses of paclitaxel (180 mg/m2)-carboplatin (AUC 5) and could not completed six courses of paclitaxel-carboplatin due to progression of the disease. These four supersensitive and four refractory specimens were used for RNA extraction, Affymetrix GeneChip analysis (Santa Clara, CA), and real-time reverse transcription-PCR (RT-PCR).
Excluding the cases used for the GeneChip, the residual 24 surgical specimens were used for immunohistochemical analysis.
Establishment of paclitaxel-resistant ovarian cancer cell lines. Using a human serous ovarian cancer cell line, 2008 (provided by Dr. S.B. Howell, Department of Medicine and the Rebecca and John Moores Cancer Center, University of California-San Diego, La Jolla, CA), we developed two kinds of novel clones resistant to paclitaxel after 40 weeks as follows: 2008/PX2 cells were obtained by biweekly medium changes with 800 ng/mL paclitaxel followed by a 2-hour exposure to paclitaxel, where doses of paclitaxel were escalated stepwise to 6,200 ng/mL; 2008/PX24 cells were obtained by biweekly medium changes with chronic exposure to 2 ng/mL paclitaxel, where doses of paclitaxel were escalated stepwise to 29 ng/mL.
The resistances of these original paclitaxel-sensitive clones and newly developed paclitaxel-resistant 2008 clones were evaluated according to established methods: in vitro 2-(2-methoxy-4-nitrophenyl)-3-(4-nitrophenyl)-5-(2,4-disulfonyl)-2H-tetrazolium assay (26) and murine model in vivo (27). Briefly, for the in vitro experiments, a single-cell suspension of 2008, 2008/PX2, or 2008/PX24 in DMEM supplemented with 10% fetal bovine serum was seeded to a 96-well plate at 3,000 cells per well. Then, the cells were treated with a range of concentrations of paclitaxel and cisplatin (carboplatin is a derivative of cisplatin) from 0.00019 to 50 μmol/L with a 2-fold serial dilution. After 4 days of incubation at 37°C in a humidified incubator containing 5% CO2, 2-(2-methoxy-4-nitrophenyl)-3-(4-nitrophenyl)-5-(2,4-disulfonyl)-2H-tetrazolium reagent (Cell Counting Kit-8, Dojindo Laboratories, Tokyo, Japan) was added to each well, and the plates were further incubated for a few hours at 37°C. Finally, the absorbance at 450 nm was measured, and the antiproliferating activity of each drug was calculated using the formula: (1 − T / C) × 100 (%), where T and C represent the mean difference in absorbance at 450 nm of the cells treated with drugs (T) and that of the untreated control cells (C). The IC50 was obtained from three independent experiments (Table 1).
Cell line . | IC50 for paclitaxel* (ratio)† (nmol/L) . | IC50 for cisplatin* (ratio)† (μmol/L) . |
---|---|---|
2008 | 2.2 ± 0.54 (1.0) | 0.50 ± 0.14 (1.0) |
2008/PX2 | 200 ± 58 (92) | 0.50 ± 0.42 (1.0) |
2008/PX24 | 120 ± 22 (57) | 0.24 ± 0.07 (0.48) |
Cell line . | IC50 for paclitaxel* (ratio)† (nmol/L) . | IC50 for cisplatin* (ratio)† (μmol/L) . |
---|---|---|
2008 | 2.2 ± 0.54 (1.0) | 0.50 ± 0.14 (1.0) |
2008/PX2 | 200 ± 58 (92) | 0.50 ± 0.42 (1.0) |
2008/PX24 | 120 ± 22 (57) | 0.24 ± 0.07 (0.48) |
IC50 presents mean ± SD obtained from three independent experiments.
The numbers in parentheses indicate the ratio of IC50.
For in vivo experiments, a single-cell suspension of 2008, 2008/PX2, or 2008/PX24 (1 × 107 cells per mouse) was s.c. inoculated into the right flank of five female mice (BALB/c nu/nu). The tumor volume was estimated by two-dimensional measurements using the equation: ab2 / 2, where a and b represent tumor length and width, respectively. When the tumor volume reached 200 to 300 mm3, 40 mg/kg paclitaxel, 80 mg/kg paclitaxel, or vehicle was given i.v. once weekly for 3 weeks (vehicle: 10% Cremophore/0% ethanol/80% saline).
RNA extraction. Cryostat sections containing >80% cancer cells were microdissected and prepared as tumor specimens. Total RNA from ovarian tumors and cell lines were isolated using the hot phenol method (28). Total RNA was isolated from three different cultures of each cell line. We also the scraped ovarian surface epithelium from three menopausal patients with leiomyoma of the uterus who underwent total hysterectomy and bilateral salpingo-oophorectomy with informed consent, and the ovarian surface epithelium was immortalized by SV40 T antigen alone and with SV40 T antigen/human telomerase reverse transcriptase. All six immortalized cell lines are nontumorigenic, and immunocytochemical analysis showed a similar staining pattern to normal ovarian surface epithelium.8
In preparation.
Microarray. Human genome-wide gene expression was examined using the Human Genome U133 Array (HG-U133 Set: GeneChip), which contains ∼45,000 probe sets, representing >39,000 transcripts derived from ∼33,000 well-substantiated human genes (http://www.affymetrix.com/products/arrays/specific/hgu133.affx). Double-stranded cDNA was synthesized, and the cDNA was subjected to in vitro transcription in the presence of biotinylated nucleotide triphosphates. The biotinylated cRNA (10 μg) was hybridized with a probe array for 16 hours at 45°C, and the hybridized biotinylated cRNA was stained with streptavidin-phycoerythrin and then scanned with a Gene Array Scanner. The fluorescence intensity of each probe was quantified using a computer program, GeneChip Analysis Suite 5.0 (Affymetrix). The expression level of a single RNA was determined as the average fluorescence intensity among the intensities obtained by 11-paired (perfect-matched and single nucleotide–mismatched) probes consisting of 25-mer oligonucleotides. If the intensities of mismatched probes were very high, gene expression was judged to be absent even if a high average fluorescence was obtained with the Microarray Analysis Suite 5.0 program. The data were processed with Affymetrix's default variables, except for scaling (target intensity, 1,000), without normalization procedures to calculate the level of gene expression as the signal.
Quantitative real-time reverse transcription-PCR. Taqman reverse transcription reagents (Applied Biosystems, Foster City, CA) were applied for cDNA synthesis. The SYBR Green reagents kit (Applied Biosystems) was used for quantitative real-time RT-PCR analysis and done according to the manufacturer's recommendations. During RT-PCR, reactions were continuously monitored with an ABI Prism 7700 Sequence Detector (Applied Biosystems). Glyceraldehyde-3-phosphate dehydrogenase messages were used as the internal control. Primers for indoleamine 2,3-dioxygenase (IDO) and glyceraldehyde-3-phosphate dehydrogenase were purchased from Applied Biosystems.
Immunohistochemical analysis. For the immunohistochemical study, formalin-fixed, paraffin-embedded sections were used. Immunostaining was done using the labeled streptavidin-biotin peroxidase complex method with the Ventana auto-immunostaining system (Ventana Japan, Yokoyama, Japan). Murine monoclonal antibody against human IDO (1:1,000; ref. 29) was used. The antigen retrieval procedure was done with a microwave oven in DAKO antigen retrieval solution for 10 minutes at 95°C to efficiently stain the sample. The sections (DAKO Cytomation, Glostrup, Denmark) were developed with 3,3′-diaminobenzidine with 0.3% H2O2 and counterstained with hematoxylin. We used surgical specimens that were analyzed with the GeneChip and real-time RT-PCR as positive and negative controls. All of them showed consistent expression of IDO as the results of mRNA expression by real-time RT-PCR. Positive and negative controls were run in parallel for every stain.
Statistics. Hierarchical clustering was analyzed with Spotfire software version 8.0 (Spotfire, Somerville, MA). The Z-score (i.e., the SD from the normal mean value of raw data transformed by log2 in each gene) was used for normalization. First, all genes were included for hierarchical clustering. Second, to adjust the significant level to account for multiple testing in the data sets, permutation tests were applied for gene screening to detect differential expression between chemoresistant and chemosensitive cell lines and patients' tumors. The distribution of maximum t statistics based on 10,000 random permutations was compared with the observed values to determine the P and its 95% confidence interval for each gene using Stata 8.0 (Stata Corp., College Station, TX). Finally, these screened genes were recomputed with hierarchical clustering under sample sets of cell lines and patients' tumors, cell lines alone, and patients' tumors alone.
The association between the stage of cancer and the staining pattern was analyzed with the χ2 test. Survival curves of the patients were compared using the Kaplan-Meier method. These analyses were done by the log-rank test using Stata 8.0.
Results
Establishment of paclitaxel-resistant ovarian cancer cell lines. After 40 weeks of exposure to paclitaxel, the 2-(2-methoxy-4-nitrophenyl)-3-(4-nitrophenyl)-5-(2,4-disulfonyl)-2H-tetrazolium assay confirmed the development of two cell lines resistant to paclitaxel but still sensitive to cisplatin as follows: the ratio of IC50 for paclitaxel between 2008 and 2008/PX2 increased to 92, whereas that for cisplatin remained at 1.0; the ratio of IC50 for paclitaxel between 2008 and 2008/PX24 was 57, whereas that for cisplatin was 0.48 (Table 1). Thus, the degree of resistance against paclitaxel was greater in 2008/PX2 than in 2008/PX24, whereas the sensitivity against cisplatin remained the same.
Next, the resistance to paclitaxel of these new cell lines was examined using a murine in vivo model and compared with that of the parental cell line, 2008 (Fig. 1). The growth of 2008 in mice was almost completely suppressed by treatment with paclitaxel at 40 and 80 mg/kg (left), whereas at the same doses of paclitaxel the growth of 2008/PX2 and 2008/PX24 was only partially suppressed (middle and right). Thus, the two new cell lines were more resistant to paclitaxel than 2008 both in vitro and in vivo.
Screening with gene expression profiling. All cell lines (2008, 2008/PX2, and 2008/PX24) and eight surgical tumors from patients (four supersensitive and four refractory) were simultaneously analyzed under hierarchical clustering using all of gene expression data (Fig. 2). Although the cell lines and surgical tumors were clearly differentiated, the nature of the chemosensitivity or chemoresistance was independent of the clusters created by the analysis.
Then, the permutation tests were applied at a cutoff point of 0.05 to screen genes that differentially expressed chemosensitivity and chemoresistance, including both cell lines and surgical tumors. As a result, 44 genes (P < 0.05) were selected as candidates associated with chemoresistance or chemosensitivity and reanalyzed with hierarchical clustering (Fig. 3A). The 44 genes were classified into major two clusters: 17 kinds of genes were down-regulated (Table 2A) representing cluster A in Fig. 3A, whereas 27 genes were up-regulated representing cluster B in Fig. 3A in both chemoresistant cell lines and resistant surgical tumors (Table 2B). Furthermore, we repeated hierarchical clustering restricted to either cell lines alone (Fig. 3B) or surgical tumors alone (Fig. 3C).
ID . | Gene . | P* . | 95% Confidence interval* . | R/S ratio both† . | R/S ratio cell line‡ . | R/S ratio patients' tumor§ . | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Permutation analyses of 17 genes down-regulated in both chemoresistant cell lines and resistant tumor specimens | ||||||||||||
200963_x_at | Ribosomal protein L31 | 0.0076 | 0.0060-0.0095 | 0.83 | 0.93 | 0.75 | ||||||
200008_s_at | GDP dissociation inhibitor 2 | 0.0015 | 0.0008-0.0025 | 0.68 | 0.91 | 0.50 | ||||||
200052_s_at | Interleukin enhancer binding factor 2, 45 kDa | 0.0080 | 0.0063-0.0099 | 0.68 | 0.96 | 0.46 | ||||||
200008_s_at | GDP dissociation inhibitor 2 | 0.0086 | 0.0069-0.0106 | 0.67 | 0.89 | 0.51 | ||||||
200036_s_at | Ribosomal protein L10a | 0.0326 | 0.0292-0.0363 | 0.75 | 0.75 | 0.77 | ||||||
201105_at | Lectin, galactoside binding, soluble, 1 (galectin 1) | 0.0077 | 0.0061-0.0096 | 0.53 | 0.41 | 0.66 | ||||||
200741_s_at | Ribosomal protein S27 (metallopanstimulin 1) | 0.0318 | 0.0284-0.0354 | 0.78 | 0.92 | 0.69 | ||||||
221729_at | Collagen, type V, α2 | 0.0040 | 0.0029-0.0054 | 0.37 | 0.50 | 0.36 | ||||||
218350_s_at | Geminin, DNA replication inhibitor | 0.0081 | 0.0064-0.0101 | 0.43 | 0.93 | 0.28 | ||||||
204031_s_at | DKFZp564J157 protein | 0.0213 | 0.0186-0.0243 | 0.69 | 0.79 | 0.60 | ||||||
201524_x_at | Ubiquitin-conjugating enzyme E2N (UBC13 homologue, yeast)∥ | 0.0219 | 0.0191-0.0250 | 0.83 | 0.98 | 0.67 | ||||||
200934_at | DEK oncogene (DNA binding) | 0.0295 | 0.0263-0.0330 | 0.63 | 0.75 | 0.52 | ||||||
200081_s_at | Ribosomal protein S6 | 0.0329 | 0.0295-0.0366 | 0.69 | 0.60 | 0.80 | ||||||
200081_s_at | Ribosomal protein S6 | 0.0304 | 0.0271-0.0340 | 0.69 | 0.61 | 0.80 | ||||||
208782_at | Follistatin-like 1∥ | 0.0314 | 0.0281-0.0350 | 0.59 | 0.79 | 0.42 | ||||||
201054_at | Heterogeneous nuclear ribonucleoprotein A0 | 0.0229 | 0.0201-0.0260 | 0.69 | 0.56 | 0.86 | ||||||
201518_at | Chromobox homologue 1 (HP1β homologue Drosophila) | 0.0390 | 0.0353-0.0430 | 0.63 | 0.88 | 0.43 | ||||||
Permutation analyses of 17 genes up-regulated in both chemoresistant cell lines and resistant tumor specimens | ||||||||||||
205812_s_at | Sulfotransferase family, cytosolic, 1C, member 2 | 0.0016 | 0.001-0.003 | 1.58 | 1.59 | 1.63 | ||||||
201358_s_at | Coatomer protein complex, subunit β | 0.0028 | 0.0019-0.0040 | 1.66 | 1.50 | 1.90 | ||||||
200063_s_at | Nucleophosmin (nucleolar phosphoprotein B23, numatrin) | 0.0052 | 0.0039-0.0068 | 1.24 | 1.22 | 1.24 | ||||||
204386_s_at | Mitochondrial ribosomal protein 63 | 0.0111 | 0.0091-0.0134 | 1.30 | 1.21 | 1.41 | ||||||
224778_s_at | Homo sapiens, clone IMAGE:5259584, mRNA | 0.0103 | 0.0084-0.0125 | 1.56 | 1.70 | 1.41 | ||||||
200039_s_at | Proteasome (prosome, macropain) subunit, βtype, 2¶ | 0.0056 | 0.0042-0.0073 | 1.45 | 1.33 | 1.52 | ||||||
201288_at | Rho-GDP dissociation inhibitor β | 0.0078 | 0.0062-0.010 | 3.13 | 4.45 | 2.23 | ||||||
225301_s_at | Myosin VB | 0.0074 | 0.0058-0.0093 | 1.42 | 1.32 | 1.53 | ||||||
200055_at | TAF10 RNA polymerase II, TATA box binding protein-associated factor, 30 kDa¶ | 0.0169 | 0.0145-0.0196 | 1.55 | 1.61 | 1.52 | ||||||
200055_at | TAF10 RNA polymerase II, TATA box binding protein-associated factor, 30 kDa | 0.0148 | 0.0125-0.0174 | 1.59 | 1.73 | 1.52 | ||||||
202442_at | Adaptor-related protein complex 3, σ1 subunit | 0.0022 | 0.0014-0.0033 | 1.37 | 1.43 | 1.32 | ||||||
208709_s_at | Nardilysin (N-arginine dibasic convertase) | 0.0409 | 0.0371-0.0450 | 1.43 | 1.20 | 1.66 | ||||||
201046_s_at | RAD23 homologue A (Saccharomyces cerevisiae) | 0.0409 | 0.0371-0.0450 | 1.57 | 1.08 | 2.32 | ||||||
207988_s_at | Actin-related protein 2/3 complex, subunit 2, 34 kDa | 0.0395 | 0.0358-0.0435 | 1.42 | 1.29 | 1.54 | ||||||
206445_s_at | HMT1 hnRNP methyltransferase-like 2 (S. cerevisiae) | 0.0421 | 0.0382-0.0462 | 1.44 | 1.28 | 1.61 | ||||||
217939_s_at | Hypothetical protein FLJ20080 | 0.0419 | 0.0381-0.0460 | 1.41 | 1.28 | 1.58 | ||||||
201470_at | Glutathione S-transferase ω1 | 0.0012 | 0.0006-0.0021 | 1.45 | 1.09 | 1.91 | ||||||
202060_at | SH2 domain binding protein 1 (tetratricopeptide repeat containing) | 0.0313 | 0.0280-0.0349 | 1.75 | 1.64 | 1.87 | ||||||
200005_at | Eukaryotic translation initiation factor 3, subunit 7ζ, 66/67 kDa | 0.0112 | 0.0092-0.0135 | 1.66 | 1.15 | 2.35 | ||||||
226622_at | Hypothetical protein FLJ14408 | 0.0078 | 0.0062-0.0097 | 3.39 | 1.77 | 4.85 | ||||||
202090_s_at | Ubiquinol-cytochrome c reductase (6.4 kDa) subunit¶ | 0.0189 | 0.0163-0.0218 | 1.33 | 1.10 | 1.61 | ||||||
203190_at | NADH dehydrogenase (ubiquinone) Fe-S protein 8, 23 kDa (NADH-coenzyme Q reductase) | 0.0289 | 0.0257-0.0324 | 1.59 | 1.41 | 1.74 | ||||||
201077_s_at | NHP2 nonhistone chromosome protein 2-like 1 (S. cerevisiae)¶ | 0.0418 | 0.0380-0.0459 | 1.34 | 1.20 | 1.43 | ||||||
210029_at | IDO¶ | 0.0413 | 0.0375-0.0454 | 2.13 | 1.34 | 3.20 | ||||||
209374_s_at | Immunoglobulin heavy constant μ¶ | 0.0290 | 0.0258-0.0325 | 4.92 | 2.17 | 23.07 | ||||||
218261_at | Adaptor-related protein complex 1, μ2 subunit | 0.0416 | 0.0378-0.0457 | 2.03 | 1.47 | 2.62 | ||||||
204070_at | Retinoic acid receptor responder (tazarotene induced) 3 | 0.0185 | 0.0160-0.0213 | 2.47 | 1.63 | 3.40 |
ID . | Gene . | P* . | 95% Confidence interval* . | R/S ratio both† . | R/S ratio cell line‡ . | R/S ratio patients' tumor§ . | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Permutation analyses of 17 genes down-regulated in both chemoresistant cell lines and resistant tumor specimens | ||||||||||||
200963_x_at | Ribosomal protein L31 | 0.0076 | 0.0060-0.0095 | 0.83 | 0.93 | 0.75 | ||||||
200008_s_at | GDP dissociation inhibitor 2 | 0.0015 | 0.0008-0.0025 | 0.68 | 0.91 | 0.50 | ||||||
200052_s_at | Interleukin enhancer binding factor 2, 45 kDa | 0.0080 | 0.0063-0.0099 | 0.68 | 0.96 | 0.46 | ||||||
200008_s_at | GDP dissociation inhibitor 2 | 0.0086 | 0.0069-0.0106 | 0.67 | 0.89 | 0.51 | ||||||
200036_s_at | Ribosomal protein L10a | 0.0326 | 0.0292-0.0363 | 0.75 | 0.75 | 0.77 | ||||||
201105_at | Lectin, galactoside binding, soluble, 1 (galectin 1) | 0.0077 | 0.0061-0.0096 | 0.53 | 0.41 | 0.66 | ||||||
200741_s_at | Ribosomal protein S27 (metallopanstimulin 1) | 0.0318 | 0.0284-0.0354 | 0.78 | 0.92 | 0.69 | ||||||
221729_at | Collagen, type V, α2 | 0.0040 | 0.0029-0.0054 | 0.37 | 0.50 | 0.36 | ||||||
218350_s_at | Geminin, DNA replication inhibitor | 0.0081 | 0.0064-0.0101 | 0.43 | 0.93 | 0.28 | ||||||
204031_s_at | DKFZp564J157 protein | 0.0213 | 0.0186-0.0243 | 0.69 | 0.79 | 0.60 | ||||||
201524_x_at | Ubiquitin-conjugating enzyme E2N (UBC13 homologue, yeast)∥ | 0.0219 | 0.0191-0.0250 | 0.83 | 0.98 | 0.67 | ||||||
200934_at | DEK oncogene (DNA binding) | 0.0295 | 0.0263-0.0330 | 0.63 | 0.75 | 0.52 | ||||||
200081_s_at | Ribosomal protein S6 | 0.0329 | 0.0295-0.0366 | 0.69 | 0.60 | 0.80 | ||||||
200081_s_at | Ribosomal protein S6 | 0.0304 | 0.0271-0.0340 | 0.69 | 0.61 | 0.80 | ||||||
208782_at | Follistatin-like 1∥ | 0.0314 | 0.0281-0.0350 | 0.59 | 0.79 | 0.42 | ||||||
201054_at | Heterogeneous nuclear ribonucleoprotein A0 | 0.0229 | 0.0201-0.0260 | 0.69 | 0.56 | 0.86 | ||||||
201518_at | Chromobox homologue 1 (HP1β homologue Drosophila) | 0.0390 | 0.0353-0.0430 | 0.63 | 0.88 | 0.43 | ||||||
Permutation analyses of 17 genes up-regulated in both chemoresistant cell lines and resistant tumor specimens | ||||||||||||
205812_s_at | Sulfotransferase family, cytosolic, 1C, member 2 | 0.0016 | 0.001-0.003 | 1.58 | 1.59 | 1.63 | ||||||
201358_s_at | Coatomer protein complex, subunit β | 0.0028 | 0.0019-0.0040 | 1.66 | 1.50 | 1.90 | ||||||
200063_s_at | Nucleophosmin (nucleolar phosphoprotein B23, numatrin) | 0.0052 | 0.0039-0.0068 | 1.24 | 1.22 | 1.24 | ||||||
204386_s_at | Mitochondrial ribosomal protein 63 | 0.0111 | 0.0091-0.0134 | 1.30 | 1.21 | 1.41 | ||||||
224778_s_at | Homo sapiens, clone IMAGE:5259584, mRNA | 0.0103 | 0.0084-0.0125 | 1.56 | 1.70 | 1.41 | ||||||
200039_s_at | Proteasome (prosome, macropain) subunit, βtype, 2¶ | 0.0056 | 0.0042-0.0073 | 1.45 | 1.33 | 1.52 | ||||||
201288_at | Rho-GDP dissociation inhibitor β | 0.0078 | 0.0062-0.010 | 3.13 | 4.45 | 2.23 | ||||||
225301_s_at | Myosin VB | 0.0074 | 0.0058-0.0093 | 1.42 | 1.32 | 1.53 | ||||||
200055_at | TAF10 RNA polymerase II, TATA box binding protein-associated factor, 30 kDa¶ | 0.0169 | 0.0145-0.0196 | 1.55 | 1.61 | 1.52 | ||||||
200055_at | TAF10 RNA polymerase II, TATA box binding protein-associated factor, 30 kDa | 0.0148 | 0.0125-0.0174 | 1.59 | 1.73 | 1.52 | ||||||
202442_at | Adaptor-related protein complex 3, σ1 subunit | 0.0022 | 0.0014-0.0033 | 1.37 | 1.43 | 1.32 | ||||||
208709_s_at | Nardilysin (N-arginine dibasic convertase) | 0.0409 | 0.0371-0.0450 | 1.43 | 1.20 | 1.66 | ||||||
201046_s_at | RAD23 homologue A (Saccharomyces cerevisiae) | 0.0409 | 0.0371-0.0450 | 1.57 | 1.08 | 2.32 | ||||||
207988_s_at | Actin-related protein 2/3 complex, subunit 2, 34 kDa | 0.0395 | 0.0358-0.0435 | 1.42 | 1.29 | 1.54 | ||||||
206445_s_at | HMT1 hnRNP methyltransferase-like 2 (S. cerevisiae) | 0.0421 | 0.0382-0.0462 | 1.44 | 1.28 | 1.61 | ||||||
217939_s_at | Hypothetical protein FLJ20080 | 0.0419 | 0.0381-0.0460 | 1.41 | 1.28 | 1.58 | ||||||
201470_at | Glutathione S-transferase ω1 | 0.0012 | 0.0006-0.0021 | 1.45 | 1.09 | 1.91 | ||||||
202060_at | SH2 domain binding protein 1 (tetratricopeptide repeat containing) | 0.0313 | 0.0280-0.0349 | 1.75 | 1.64 | 1.87 | ||||||
200005_at | Eukaryotic translation initiation factor 3, subunit 7ζ, 66/67 kDa | 0.0112 | 0.0092-0.0135 | 1.66 | 1.15 | 2.35 | ||||||
226622_at | Hypothetical protein FLJ14408 | 0.0078 | 0.0062-0.0097 | 3.39 | 1.77 | 4.85 | ||||||
202090_s_at | Ubiquinol-cytochrome c reductase (6.4 kDa) subunit¶ | 0.0189 | 0.0163-0.0218 | 1.33 | 1.10 | 1.61 | ||||||
203190_at | NADH dehydrogenase (ubiquinone) Fe-S protein 8, 23 kDa (NADH-coenzyme Q reductase) | 0.0289 | 0.0257-0.0324 | 1.59 | 1.41 | 1.74 | ||||||
201077_s_at | NHP2 nonhistone chromosome protein 2-like 1 (S. cerevisiae)¶ | 0.0418 | 0.0380-0.0459 | 1.34 | 1.20 | 1.43 | ||||||
210029_at | IDO¶ | 0.0413 | 0.0375-0.0454 | 2.13 | 1.34 | 3.20 | ||||||
209374_s_at | Immunoglobulin heavy constant μ¶ | 0.0290 | 0.0258-0.0325 | 4.92 | 2.17 | 23.07 | ||||||
218261_at | Adaptor-related protein complex 1, μ2 subunit | 0.0416 | 0.0378-0.0457 | 2.03 | 1.47 | 2.62 | ||||||
204070_at | Retinoic acid receptor responder (tazarotene induced) 3 | 0.0185 | 0.0160-0.0213 | 2.47 | 1.63 | 3.40 |
P value and 95% confidence interval obtained by permutation t tests.
Mean raw data of chemoresistant cell lines plus resistant patients' tumors/mean raw data of chemosensitive cell lines plus supersensitive patients' tumors.
Mean raw data of chemoresistant cell lines/mean raw data of chemosensitive cell lines.
Mean raw data of resistant patients' tumors/mean raw data of supersensitive patients' tumors.
Reduced expression of the gene in chemoresistant cell lines and resistant patients' tumors was confirmed with RT-PCR.
Increased expression of the gene in chemoresistant cell lines and resistant patients' tumors was confirmed with RT-PCR.
Firstly, we assigned priority to 27 genes among 44 genes by either being reported as genes associated with carcinogenesis or being associated with notable pathways. Then, we selected eight genes that showed reproducible results by real-time RT-PCR comparing with the results of GeneChip analysis (Table 3). In particular, IDO was highly and consistently expressed in both chemoresistant cell lines and tumors from refractory patients but not in chemosensitive cell lines and tumors (Fig. 4). This finding was most prominent among these eight genes.
Gene . | Function . |
---|---|
IDO | Tryptophan degradation |
Immunoglobulin heavy constant μ | Immunity |
Proteasome | Cleaving peptides in an ATP/ubiquitin–dependent process |
Ubiquitin-conjugating enzyme E2N | Ubiquitin-conjugating enzyme |
Ubiquinol-cytochrome c reductase | Mitochondrial respiratory chain |
TAF10 | RNA polymerase II |
NHP2 nonhistone chromosome protein 2-like 1 | RNA binding protein |
Follistatin-like1 | Cell growth |
Gene . | Function . |
---|---|
IDO | Tryptophan degradation |
Immunoglobulin heavy constant μ | Immunity |
Proteasome | Cleaving peptides in an ATP/ubiquitin–dependent process |
Ubiquitin-conjugating enzyme E2N | Ubiquitin-conjugating enzyme |
Ubiquinol-cytochrome c reductase | Mitochondrial respiratory chain |
TAF10 | RNA polymerase II |
NHP2 nonhistone chromosome protein 2-like 1 | RNA binding protein |
Follistatin-like1 | Cell growth |
Expression of indoleamine 2,3-dioxygenase protein in pathologic specimens. Expression of IDO protein was further confirmed using pathologic specimens obtained from 24 patients with stage III or IV serous ovarian cancer, excluding samples used for GeneChip analysis. The staining patterns were classified as negative (n = 7; Fig. 5F), sporadic (n = 12; Fig. 5E), focal (n = 3; Fig. 5D), or diffuse (n = 2; Fig. 5C). There was no association between stage of cancer and staining pattern using the χ2 test.
Indoleamine 2,3-dioxygenase protein expressions and relapse-free survival. First, overall survival was compared between patients with stage III disease (n = 17) and stage IV disease (n = 7) using the log-rank test; no significant difference was noted. Next, Kaplan-Meier survival curves were generated based on the IDO staining pattern (Fig. 6). In contrast to clinical stages, staining patterns of IDO impaired survival (log-rank test, P = 0.0001). All patients classified as negative survived without relapse. The 50% survival of patients classified as sporadic, focal, and diffuse was 41, 17, and 11 months, respectively. We also established a scoring system considering both pattern and intensity, and statistical analysis showed significant differences among every score (data not shown).
Discussion
We screened and identified IDO from 39,000 transcripts as a strong prognostic factor expressed in serous ovarian cancer. Most previous works using gene expression profiling were able to identify a bulk of genes that were highly expressed or suppressed in clinical subgroups of patients, such as those with a differential prognosis (30) or a pathologic type (31). However, fewer studies have shown a single molecule that can be used to clinically distinguish specific subgroups of disease (32, 33). Although microarray technology may be powerful enough to enhance the predictive ability of the prognosis (34), the cost of this technology is still high. In this study, we used microarray technology as a screening tool to identify key molecules associated with chemoresistance in serous ovarian cancer.
Gene expression profiling of novel chemoresistant cell lines was compared with an original chemosensitive cell line to exclude individual differences. However, this approach may pick up genes associated not only with chemoresistant-specific molecules but also with the concurrent changes obtained during the 40 weeks of culture. In contrast, using differential expressions of genes using patients' cells derived from a small sample size, it may be difficult to detect chemoresistant genes, although we carefully selected eight patients who were in the same clinical stage but had a clear contrast between chemosensitive disease and chemoresistant disease in clinical settings. Few previous articles attempted to validate the results obtained from cell lines in patients' cells (35). In this study, a hierarchical clustering of gene expression profiling showed a prominent difference between cell lines and surgically resected patients' tumors but not between chemosensitivity and chemoresistance. Therefore, the permutation tests were applied to abstract chemoresistance-associated genes common to both cell lines and patients' cells. In the selected 27 genes, only 8 were confirmed with real-time RT-PCR, suggesting that the results of the GeneChip and permutation tests cutoff at 0.05 may include some false-positive information. Levels of up-regulation in IDO expression were more prominent in results of real-time RT-PCR than with the GeneChip, which may be due to differences in the methods used to quantify the amounts of RNA expression.
We were able to validate the clinical importance of IDO expression retrospectively using 24 clinical paraffin-embedded specimens, excluding cases used for GeneChip analyses. For patients with advanced serous ovarian cancer, staining patterns of IDO protein expression clearly differentiated between those with a good prognosis and those with a poor prognosis; these prognoses were not predicted by standard clinical staging. This evidence may provide credence to the strategy of starting with genome-wide screening with gene expression profiling using microarray technology, narrowing the number of genes, and ending up with a single gene to link to clinical end points.
IDO, which is a rate-limiting enzyme that catabolizes tryptophan to kynurenine, first attracted a great deal of attention because it could protect against fetal rejection due to immune surveillance (36–38). Recently, tumor cells were also shown to express IDO and to escape the immune surveillance of the host (39, 40) by degrading local tryptophan, which suppresses T cells (41, 42) and natural killer cell proliferation (43, 44). All patients who were negative for IDO survived without relapse, although the duration of survival was impaired depending on the pattern of IDO expression. This finding may be explained by the suppression of antitumor immune activities via IDO expression. On the other hand, the recurrence-free survival rate of IDO-positive patients with hepatocellular carcinoma was shown to be significantly higher than that of IDO-negative patients (45). According to their report, IDO-positive cells were identified to be tumor-infiltrating cells, not tumor cells, by immunohistochemical analysis. Although we also examined the staining pattern of tumor-infiltrating cells in the ovarian cancer portion, few cells showed positive staining. On the contrary, positive staining of tumor cells was much more prominent than that of noncancerous cells in all sporadic, focal, and diffuse patterns. Thus, the clinical significance of IDO expression being associated with prognosis in patients with serous-type ovarian cancer may not be universal to all types of cancer.
In this study, greater expression of IDO was confirmed not only in tumors from chemoresistant patients but also in chemoresistant cell lines, suggesting that IDO may affect chemosensitivity through intracellular mechanisms. Recently, IDO expression was shown to be suppressed by nitric oxide, which is known to mediate chemosensitivity in tumor cells via scavenging the production of large quantities of cytosolic superoxide anions (46). On the contrary, hypoxia-induced drug resistance seems to result, in part, from the downstream suppression of endogenous nitric oxide production (47–49). Therefore, the expression of IDO may be a parallel phenomenon to other mechanisms for chemoresistance, such as nitric oxide production, and may not cause chemoresistance directly. Just recently, Muller et al. reported that IDO inhibition cooperated with diverse chemotherapeutic agents to effectively promote the regression of established breast tumors that are refractory to chemotherapy (50). They used MMTV-Neu mice, a well-accepted transgenic mouse model of breast cancer, and showed that combining the IDO inhibitor 1-methyl-dl-tryptophan with paclitaxel resulted in a significant tumor decrease compared with paclitaxel alone (P = 0.0010). Their report supports our data, indicating that IDO is positively associated with paclitaxel resistance and impaired survival. They also indicated that Bin1 loss elevated the signal transducers and activators of transcription 1– and nuclear factor-κB–dependent expression of IDO. Nuclear factor-κB activation suppresses the apoptotic potential of chemotherapeutic agents (51). We speculate that IDO might be positively associated with paclitaxel resistance through the suppression of the apoptotic potential of paclitaxel.
Three of four supersensitive cases underwent paclitaxel-carboplatin, and one case underwent paclitaxel alone due to the hypersensitivity reaction to carboplatin. Moreover, neither 2008/PX2 nor 2008/PX24 showed cross-resistance for cisplatin. Carboplatin is a derivative of cisplatin and is a platinum compound. Using both surgical specimens and cell lines, we purified genes associated with paclitaxel resistance to prevent genes associated with platinum resistance. However, we speculate that IDO not only plays a role in paclitaxel resistance but also has an indirect effect on platinum in vivo. The latter speculation is supported by Muller et al. (50). They also showed that 1-methyl-dl-tryptophan with cisplatin also resulted in a significant tumor decrease compared with cisplatin alone. Future study should focus on the functional insights regarding the IDO gene for chemoresistance to paclitaxel by gene knockdown, such as the RNA interference technique.
In conclusion, IDO screened with the GeneChip was positively associated with paclitaxel resistance and impaired survival in patients with serous-type ovarian cancer.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Acknowledgments
We thank Mika Endo (Pharmaceutical Research Department 2, Research Division, Chugai Pharmaceutical Co., Ltd.) for assistance with in vivo experiments and Misako Arima (Department of Pathology, Jikei University School of Medicine) for technical assistance with immunohistochemistry.