
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Experimental Therapeutics, Preclinical Pharmacology |
1Laboratory of Cell Biology, Cancer Research Institute and Cancer Research Center, Seoul National University, Seoul, Korea2 Research Institute and Hospital, National Cancer Center, Goyang, Gyeonggi, Korea
| ABSTRACT |
|---|
|
|
|---|
Experimental Design: We performed global gene expression analysis in the acquired drug-resistant gastric cancer cell lines to the commonly used drugs 5-fluorouracil, doxorubicin, and cisplatin using Affymetrix HG-U133A microarray. The gene expression patterns of 10 chemoresistant gastric cancer cell lines were compared with those of four parent cell lines using fold-change and Wilcoxons test for data analysis.
Results: We identified over 250 genes differentially expressed in 5-fluorouracil-, cisplatin-, or doxorubicin-resistant gastric cancer cell lines. Our expression analysis also identified eight multidrug resistance candidate genes that were associated with resistance to two or more of the tested chemotherapeutic agents. Among these, midkine (MDK), a heparin-binding growth factor, was overexpressed in all drug-resistant cell lines, strongly suggesting that MDK might contribute to multidrug resistance in gastric cancer cells.
Conclusions: Our investigation provides comprehensive gene information associated with acquired resistance to anticancer drugs in gastric cancer cells and a basis for additional functional studies.
| INTRODUCTION |
|---|
|
|
|---|
Microarray technologies have been widely used for comprehensive gene expression analysis as well as mutation and single nucleotide polymorphism detection (7, 8, 9, 10, 11, 12, 13, 14) . In particular, large-scale microarray analysis of gene expression enables researchers to analyze simultaneous changes in thousands of genes and identify significant patterns. Because the genome-wide expression analysis in a variety of drug-resistant gastric cancer cells has not yet been performed, we have used the recently developed Affymetrix HG-U133A high-density oligonucleotide microarray for analysis of the global gene expression.
The drugs 5-fluorouracil (5-FU), doxorubicin, and cisplatin are widely used in the treatment of various malignancies: 5-FU is a well-known antimetabolite that acts as a thymidylate synthase inhibitor (3) ; doxorubicin targets topoisomerase II by interfering with the catalytic cycle (15) ; and cisplatin intercalates into DNA, leading to DNA damage in the cancer cells (16) . In this study, we have examined genes that are differentially expressed in 5-FU-, doxorubicin-, or cisplatin-induced chemoresistant gastric cancer cell lines, as compared with their drug-sensitive parent cell lines. We identified genes showing altered expression in resistant cell lines, as well as several potential multidrug resistance candidate genes that were associated with resistance to two or more of the chemotherapeutic agents.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
MTT Assay.
Sensitivities of the drug-resistant and parent cell lines to 5-FU, doxorubicin, and cisplatin were determined by MTT assay as described previously (18)
. Briefly, single-cell suspensions were obtained by pipette disaggregation of the floating cells or by trypsinization of monolayer cultures. The number of cells plated into 96 wells was determined after preliminary cell growth studies so that untreated cells were in exponential growth phase at the time of initial harvest and at the end of the 4-day incubation. Equal number of cells was inoculated into each well in RPMI 1640 supplemented with 10% fetal bovine serum. For each drug, 510 concentrations were used, covering a 35-log concentration range that was chosen to span the 50% inhibitory concentration determined by preliminary assays. After 4 days of culture, MTT (Sigma Chemical Co., St. Louis, MO) was added to each well and was incubated at 37°C for an additional 4 h. The medium was aspirated from plates leaving about 30 µl of medium in each well. Care was taken not to disturb the formazan crystals at the bottom of the wells. One hundred fifty µl of DMSO was added to each well, and the plates were placed on a shaker for 15 min to solubilize the formazan crystals. The plates were then read immediately at 540 nm on a scanning multiwell spectrophotometer (ELISA reader; Biotek Instruments Inc., Burlington, VT). All of the data points represent the mean value of a minimum of six wells. Table 1
shows the IC50 representing the drug concentration resulting in 50% growth inhibition and the relative resistance calculated from the ratio of the IC50 of the drug-resistant cell lines versus that of the parent cell lines.
RNA Preparation and Affymetrix GeneChip Hybridization.
Total RNA was extracted using the Trizol reagent (Life Technologies, Inc. Carlsbad, CA) according to the manufacturers instructions. Genes expressed in 10 drug-resistant cell lines and their parent cell lines were analyzed on a high-density oligonucleotide microarray (HG-U133A; Affymetrix, Santa Clara, CA) containing 22,282 transcripts. Target preparation and microarray processing procedures were performed as described in the Affymetrix GeneChip Expression Analysis Manual (Affymetrix, Santa Clara, CA). Briefly, the extracted total RNA was purified with an RNeasy kit (Qiagen, Valencia, CA). Twenty µg of total RNA was used to synthesize double-strand cDNA with SuperScript II reverse transcriptase (Life Technologies, Inc. Rockville, MD) and a T7-(dT)24 primer (Metabion, Germany). Then, biotinylated cRNA was synthesized from the double-stranded cDNA using the RNA Transcript Labeling kit (Enzo Life Sciences, Farmingdale, NY) and was purified and fragmented. The fragmented cRNA was hybridized to the oligonucleotide microarray, which was washed and stained with streptavidin-phycoerythrin. Scanning was performed with an Agilent Microarray Scanner (Agilent Technologies, Palo Alto, CA).
Data Analysis.
GeneChip analysis was performed based on the Affymetrix GeneChip Manual (Affymetrix Inc., Santa Clara, CA) with Microarray Analysis Suite (MAS) 5.0, Data Mining Tool (DMT) 2.0, and Microarray Database software. All of the genes represented on the GeneChip were globally normalized and scaled to a signal intensity of 500. Fold changes were calculated by comparing transcripts between parent and acquired drug-resistant cell lines. The Microarray Analysis Suite software used the Wilcoxons test to generate detected (present or absent) and changed (increased or decreased) calls, and used the calls to statistically determine whether a transcript was expressed or not, and whether it was relatively increased, decreased, or unchanged. After being filtered through a "present" call (p < 0.05), a transcript was considered differentially expressed when it satisfied one of the following two conditions: 1) by fold change, transcripts increased or decreased >1.5-fold; 2) by one-sided Wilcoxons rank test, transcripts average fold change exceeded 1.5, with an "increased" (p < 0.003) or "decreased" (p > 0.997) call. In the case of 5-FU, all transcripts meeting the above conditions in at least three of four cell lines were considered differentially expressed. Hierarchical clustering and dendrogram figures were generated using Cluster and TreeView software (http://rana.stanford.edu). A three-dimensional graph of the multidimensional scaling was created using SPSS (SPSS Inc., Chicago, IL) and SigmaPlot (SPSS Inc.). To evaluate the statistical significance of eight genes differentially expressed in two or three of the drug subsets, paired t test was performed using SPSS (SPSS Inc.).
Quantitative and Semiquantitative RT-PCR.
We selected two genes for real-time quantitative reverse transcription (RT)-PCR and seven genes for semiquantitiative RT-PCR for validation of the microarray data. Five µg of total RNA was used for creation of single-stranded cDNA using the SuperScript Preamplification System for First Strand cDNA Synthesis (Life Technologies, Inc., Rockville, MD). The cDNA was diluted and quantitatively equalized for PCR amplification. For real-time quantitative RT-PCR of MDK (midkine) and BIRC5, TaqMan PCR method using a 7900HT Sequence Detection System (Applied Biosystems, Foster City, CA) was performed according to the manufacturers instructions. We used the primers and probes provided as Assays-on-Demand Gene Expression Products (Applied Biosystems). The expression of seven genes (MDK, FOSB, NQO1, DDB2, ITGB4, ABCC1, and IGFBP2) was verified by semiquantitative RT-PCR. The primer sets for PCR amplification were as follows: MDK forward 5'-ATGCAGCACCGAGGCTTCCT-3', reverse 5'-ATCCAGGCTTGGCGTCTAGT-3'; FOSB forward 5'-GAGAGGGGAAGAGACAAAGT-3', reverse 5'-CTTCATCCTCACACAGGACT-3'; NQO1 forward 5'-TGGAGAATATTTGGGATGAG-3', reverse 5'-AATCCAGGCTAAGGAATCTC-3'; DDB2 forward 5'-GGAGATATCATGCTCTGGAA-3', reverse 5'-GGCTACTAGCAGACACATCC-3'; ITGB4 forward 5'-TTCCAAATCACAGAGGAGAC-3', reverse 5'-CTTGAGGTTGTCCAGATCAT-3'; ABCC1 forward 5'-CTGACAAGCTAGACCATGAATGT-3', reverse 5'-TCACACCAAGCCGGCGTCTTT-3'; and IGFBP2 forward 5'- TTCCAGTTCTGACACACGTA-3', reverse 5'-GACACAGGGGTTCAAAAATA -3'.
PCR was carried out with 1 µl of cDNA as follows: initial denaturation at 94°C for 5 min followed by 2530 cycles of 94°C for 30 s, 55°C for 30 s, 72°C for 1 min, followed by a final elongation at 72°C for 7 min.
Western Blot Analysis.
To investigate the correspondence between mRNA and protein of MDK, we performed Western blot analysis in SNU-620, SNU-620R-5-FU/1000, SNU-620R-CIS/2000, SNU-638, and SNU-638R-CIS/400. MDK protein was detected with an MK antibody (Santa Cruz Biotechnology, Inc., Santa Cruz, CA). Western blot analysis was performed as described previously (19)
.
| RESULTS |
|---|
|
|
|---|
Validation of Microarray Results.
To verify the expression of the genes identified in microarray experiments, either real-time quantitative RT-PCR or semiquantitative RT-PCR was performed using the same RNA as that used in the microarray analysis. We tested two genes (MDK and BIRC5) for real-time quantitative RT-PCR and seven genes (MDK, FOSB, NQO1, DDB2, ITGB4, ABCC1, and IGFBP2) for semiquantitative RT-PCR, and found that the results were in good agreement with those from the microarray data, in that the observed differences were not significant.
Differentially Expressed Genes in 5-FU-Resistant Gastric Cancer Cells.
In 5-FU-resistant gastric cancer cell lines, we identified 38 differentially expressed genes (Table 2)
, most of which are involved in cell proliferation, metabolic pathways, cell growth, and cell organization. In addition, a subset of differentially expressed genes was associated with signaling, responses to external stimuli, and cell adhesion. In particular, we observed up-regulation of cell growth regulators (NOV and IGFBP2) and genes involved in nucleobase, nucleotide, and nucleic acid metabolism (HOXB6, SF3B3, TYMS, DGUOK, and DDB2). Also, the overexpression of TYMS (thymidylate synthetase) was consistent with the action of 5-FU as a TYMS inhibitor. The signal transducers sorcin (SRI), MDK, and NQO1 were up-regulated, and the calcium-binding protein SRI was up-regulated in all of the 5-FU-resistant cell lines. A novel hypoxia-inducible factor 1 (HIF-1)-responsive gene, RTP801, recently identified as involved in both pro- and antiapoptotic activities (20)
, was down-regulated, as were several calcium ion-binding molecules (S100P, CALB2, and LDLR).
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
It was suggested that each cancer cell represents a different pattern of drug-resistance gene expression even within cells clonally derived from the same cancer, and may be expected to exhibit a considerable amount of heterogeneity with respect to drug resistance (4)
. Because our four gastric cancer parent cell lines were genetically different (17)
, the overall gene expression pattern of each resistant cell line tended to be more similar to the same-parent-originated resistant cell lines than to the same-drug-treated resistant cell lines, as shown in multidimensional scaling (Fig. 1A)
.
|
The doxorubicin-resistant cell lines provided the most differentially expressed genes. This suggests that the doxorubicin-resistant cells, although derived from different gastric cancer patients, contain more consistent molecular changes than the other drug-resistant cells. When we performed clustering analysis between all drug-resistant cell lines and the differentially expressed genes, we observed that the doxorubicin-resistant cell lines were clustered together (Fig. 1B)
. A previous study demonstrated that at different time points, etoposide-resistant melanoma cells showed stable gene deregulation, whereas cisplatin- and fotemustine-resistant cells showed substantial variation (11)
. Together with our results, this suggests that anticancer drugs targeting topoisomerase II, such as doxorubicin and etoposide, exhibit coherent resistance in different cell types and during different periods.
Cell cycle deregulation is an important molecular event in the acquisition of drug resistance (5) . Most of the genes we identified as overexpressed in doxorubicin-resistant gastric cancer cells were involved in the cell cycle, including the mitotic cell cycle-associated genes BUB1, BUB1B, CDKN3, RRM1, and RRM2. It was previously shown that cell cycle genes had an interesting expression pattern in microarray experiment after exposure to doxorubicin in MCF-7 cells and a subset of these genes was also constitutively overexpressed in MCF-7 doxorubicin-resistant cells (7) . Thus, these results suggest that cell cycle genes might significantly contribute to the doxorubicin resistance. In our microarray result, the up-regulations of BIRC5 were observed in doxorubicin-resistant cells. BIRC5, also termed survivin, is a protein responsible for inhibiting apoptosis and preventing cell death (24) . Negative correlation between BIRC5 (survivin) and 5-FU derivatives in untreated human cancer cell lines was suggested (9) , implying that the overexpression of BIRC5 might be associated with drug-resistance. RRM2, one of two ribonucleotide reductase subunits, is a rate-limiting enzyme in DNA synthesis and repair (25) . Previously, it was reported that the overexpression of RRM2 mRNA and protein was found in gemcitabine-resistant cells, implicating the gene in drug resistance (25) . Finally, the ATP-binding cassette (ABC) transporter ABCC1 (MRP1) was up-regulated in all doxorubicin-resistant gastric cancer cells, as compared with ABCC1 expressions in the parent cells. Drug resistance through drug efflux pumps has been well described, and the ABC transporters are also known as energy-dependent drug efflux pumps (4) . In both our doxorubicin-resistant cells and their parent gastric cancer cells, ABCB1 (MDR1) expressions were not detectable. It was previously demonstrated that the levels of MDR1 RNA was relatively low in gastric cancer cell lines, whereas intermediate or high levels were present in colorectal carcinoma cell lines (6 , 18) .
A total of 38 genes were differentially expressed in 5-FU-resistant gastric cancer cell lines. The relationship between the overexpression of TYMS and 5-FU resistance is well-characterized (3) , and was noted in this investigation except in one case (SNU-638R5-FU/50000). Although SNU-638R5-FU/50000 represented the highest degree of resistance to 5-FU (>12,500-fold), TYMS expression was slightly decreased in these cells per our microarray data. This may be explained by the proposition that high doses of 5-FU reduce the activity of TYMS (26) . Because the SNU-638R5-FU/50000 cells were cultured at high 5-FU concentrations (380 µM) over almost 2 years, it is possible that TYMS transcription levels had been down-regulated, and that the acquired 5-FU resistance was not dependent on the TYMS mechanism. A more possible explanation for the TYMS down-regulation in SNU-638R-5-FU/50000 cells is that TYMS may be mutated and has a decreased affinity for 5-fluoro-dUMP. In our microarray data, DYPD expression was not significant in 5-FU resistant cells, although a previous report suggested that DYPD expression was negatively correlated with 5-FU sensitivity in cancer cells (21) . Another up-regulated gene is that for SRI, a soluble resistance-related calcium-binding protein known to be overexpressed in various drug-resistant cancer cells (27 , 28) , although this is the first identification of its association with 5-FU resistance. All four of the 5-FU-resistant cell lines showed up-regulated SRI expression, and, although increased expression of SRI is believed to result from gene amplification (28) , it is not yet known how SRI confers multidrug resistance. CLU (clusterin), also overexpressed in all of the 5-FU-resistant cell lines, is a ubiquitous glycoprotein that is highly overexpressed in some normal and malignant tissues undergoing apoptosis (29) . The overexpression of CLU has been shown to inhibit apoptotic cell death in cisplatin-treated renal cell carcinoma, suggesting that CLU confers a chemoresistance phenotype through its antiapoptotic activity (29 , 30) . Moreover, it has been demonstrated that suppression of CLU expression by an antisense oligonucleotide had a chemosensitizing effect (29 , 30) . From these results, it can be suggested that CLU is a potent candidate gene for gastric cancer cell drug resistance.
In cisplatin-resistant gastric cancer cells, the differentially expressed genes were mainly associated with stress response, transport, cell cycle, and metabolism. ATP7A, a copper-transporting P-type adenosine triphosphatase, was overexpressed in cisplatin-resistant gastric cancer cells. The closely related ATP7B is well known for conferring resistance to cisplatin (31)
, but little is known about the association of ATP7A with cisplatin resistance. ATP7A is highly homologous with ATP7B and has a similar export function for copper. A recent report (32)
showed ATP7A overexpression in a cisplatin-resistant ovarian cancer cell line. These results suggest a potential role for ATP7A in cisplatin resistance, as well as ATP7B, although the mechanism has not yet been elucidated (32)
. The Fos gene family members, FOS and FOSB, showed up-regulation in two of three cisplatin-resistant cells (not listed in Table 4
), consistent with previous reports suggesting a relationship between Fos gene expression and drug resistance (33)
.
In this study, we have identified eight genes whose differential expression was associated with resistance to more than one of the three chemotherapeutic drugs tested. Of these, the MDK gene was consistently overexpressed in all 5-FU-, doxorubicin-, and cisplatin-resistant gastric cancer cells compared with the expressions in their parent cells, strongly suggesting that it may contribute to multidrug resistance. The correspondence between mRNA and protein of MDK was observed in SNU-620R-5-FU/1000 and SNU-638R-CIS/400 by Western blot analysis. MDK, a heparin-binding growth factor, is associated with promoting neuronal survival, inducing neurite outgrowth (34, 35, 36) , and has been implicated in carcinogenesis and angiogenesis. Found at low levels in normal adult tissue, MDK is frequently overexpressed in esophageal, gastric, colon, pancreatic, hepatocellular, lung, breast, urinary bladder carcinoma, neuroblastoma and Wilms tumors (35 , 36) . Transfection of an MDK antisense oligodeoxynucleotide into CMT-93 mouse rectal carcinoma cells markedly suppressed tumor growth, suggesting that MDK could be a potent target for cancer therapy (36 , 37) . MDK has cytoprotective activity and has rescued cisplatin-induced apoptotic cell death in both murine kidney and cultured G401 cells, implying that MDK confers a cellular growth advantage by functioning as an antiapoptotic factor (34) . Recently, the overexpression of HB-EGFGF (heparin-binding epidermal growth factor-like growth factor) was observed in the cisplatin-resistant and 5-FU-resistant gastric tumors by microarray analysis, suggesting that HB-EGFGF is a candidate chemoresistant-related gene of gastric cancer (38) . Interestingly, these results support that heparin-binding growth factors such as MDK or HB-EGFGF might be one of the important genes contributing to chemoresistance in gastric cancers. Accordingly, we hypothesize that the overexpressed MDK can contribute to the growth progression of gastric cancer cells with an acquired multidrug resistance and may be a potent target to restore chemosensitivity. To clarify the relationship between MDK overexpression and multidrug resistance, it will be necessary to further investigate whether the inhibition of MDK overexpression in drug-resistant cancer cells leads to recovered sensitivity.
In summary, we have identified genes that are differentially expressed in 5-FU-, doxorubicin-, and cisplatin-resistant gastric cancer cells. Some of the identified genes were previously known to be associated with drug-resistance. To date, most studies on drug resistance in gastric cancer have focused on candidate gene work targeting a limited number of drug resistance genes identified in other cancers. In contrast, our use of microarray technology for a genome-wide screen of 10 drug-resistant gastric cancer cell lines has provided useful information on differentially expressed genes and possible new candidate multidrug resistance genes in gastric cancer cells. Although further in vivo validations for the identified genes are required because the materials used in the study were induced drug-resistant cells in vitro, this important information may lead to the discovery of new drug resistance targets, and perhaps to the development of better cancer chemotherapy strategies.
| FOOTNOTES |
|---|
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.
Notes: Drs. H. C. Kang, I-J. Kim, and J-H. Park contributed equally to this work.
Requests for reprints: Jae-Gahb Park, National Cancer Center, 809 Madu-dong, Ilsan-gu, Goyang, Gyeonggi, 411-764, Korea. Phone: 82-31-920-1501; Fax: 82-31-920-1511; E-mail: park{at}ncc.re.kr
Received 7/11/03; revised 10/20/03; accepted 10/21/03.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. D. Minna, L. Girard, and Y. Xie Tumor mRNA Expression Profiles Predict Responses to Chemotherapy J. Clin. Oncol., October 1, 2007; 25(28): 4329 - 4336. [Full Text] [PDF] |
||||
![]() |
N. Keshelava, E. Davicioni, Z. Wan, L. Ji, R. Sposto, T. J. Triche, and C. P. Reynolds Histone Deacetylase 1 Gene Expression and Sensitization of Multidrug-Resistant Neuroblastoma Cell Lines to Cytotoxic Agents by Depsipeptide J Natl Cancer Inst, July 18, 2007; 99(14): 1107 - 1119. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Boyer, W. L. Allen, E. G. McLean, P. M. Wilson, A. McCulla, S. Moore, D. B. Longley, C. Caldas, and P. G. Johnston Pharmacogenomic identification of novel determinants of response to chemotherapy in colon cancer. Cancer Res., March 1, 2006; 66(5): 2765 - 2777. [Abstract] [Full Text] [PDF] |
||||
![]() |
I.-J. Kim, H. C. Kang, S.-G. Jang, K. Kim, S.-A Ahn, H.-J. Yoon, S. N. Yoon, and J.-G. Park Oligonucleotide microarray analysis of distinct gene expression patterns in colorectal cancer tissues harboring BRAF and K-ras mutations Carcinogenesis, March 1, 2006; 27(3): 392 - 404. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. C. Mallory, G. Crudden, A. Oliva, C. Saunders, A. Stromberg, and R. J. Craven A Novel Group of Genes Regulates Susceptibility to Antineoplastic Drugs in Highly Tumorigenic Breast Cancer Cells Mol. Pharmacol., December 1, 2005; 68(6): 1747 - 1756. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Nakatsu, Y. Yoshida, K. Yamazaki, T. Nakamura, S. Dan, Y. Fukui, and T. Yamori Chemosensitivity profile of cancer cell lines and identification of genes determining chemosensitivity by an integrated bioinformatical approach using cDNA arrays Mol. Cancer Ther., March 1, 2005; 4(3): 399 - 412. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Michael and M.M. Doherty Tumoral Drug Metabolism: Overview and Its Implications for Cancer Therapy J. Clin. Oncol., January 1, 2005; 23(1): 205 - 229. [Abstract] [Full Text] [PDF] |
||||
![]() |
D.-k. Lee, Y. H. Kim, J.-S. Kim, and W. Seol Induction and characterization of taxol-resistance phenotypes with a transiently expressed artificial transcriptional activator library Nucleic Acids Res., August 10, 2004; 32(14): e116 - e116. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |