| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Biology and Genetics |
Division of Molecular Pharmacology, Cancer Chemotherapy Center, Japanese Foundation for Cancer Research, Toshima-ku, Tokyo 170-8455 [S. D., K. Y., T. Y.]; Laboratory of Molecular Medicine, Human Genome Center, Institute of Medical Science, The University of Tokyo, Minato-ku, Tokyo 108-8639 [O. K., R. Y., H. Z., T. K., Y. N.]; and SNP Research Center, RIKEN (Institute of Physical and Chemical Research), Minato-ku, Tokyo 108-8639 [T. T.], Japan
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Two major factors need to be considered when regarding the efficacy of particular anticancer drugs: (a) local blood or tissue drug concentrations can be influenced by the activities of metabolic enzymes, which may activate or inactivate the drug, and actions of drug transporters; and (b) differences in individual drug efficacies may also reflect the intrinsic susceptibility of cancer cells to those anticancer drugs. In the latter case, a number of genes has been reported to influence chemosensitivity, e.g., cancer cells expressing high levels of P-glycoprotein encoded by the ABCB1 (formerly MDR-1) gene are often resistant to a large subset of anticancer drugs, including vincristine, etoposide, and paclitaxel (taxol; Ref. 1 ). However, it has become obvious that the susceptibility of cancer cells to particular anticancer drugs cannot be predicted by a single factor but is determined by many factors that influence overall sensitivity. Therefore, to establish an appropriate protocol for the prediction of chemosensitivity, we need to accumulate information on the sets of genes that pharmacologically characterize cancer cells on treatment with particular anticancer drugs.
Taking advantage of a panel of 39 well-characterized human cancer cell lines (2) and a cDNA microarray system consisting of 9216 genes (3) , we attempted to reveal genes associated with cancer cell chemosensitivity. Here, we report the construction of an integrated database of gene expression profiles and drug sensitivity patterns for 39 human cancer cell lines and the identification of gene sets that are likely to be involved in chemosensitivity.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Growth Inhibition Assay and Data Processing.
Growth inhibition was assessed as changes in total cellular protein after 48 h of drug treatment using a sulforhodamine B assay. The GI50 was calculated as described previously (2
, 4) .
Identification of Gene Expression Profiles by cDNA Microarray.
Cell lines grown as monolayers to log phase were washed twice with PBS, and total RNA was extracted with TRIzol reagent (Life Technologies, Inc.). After treatment with DNase I (Boehringer Mannheim) to remove contaminating genomic DNA, poly(A)+ RNA was extracted using an mRNA Purification Kit (Amersham Pharmacia Biotech). Probe cDNA-labeling reactions for each sample were performed as described previously (5)
, except that an oligodeoxythymidylate primer was used instead of random hexamers. Each sample was reverse transcribed in the presence of Cy5-labeled dCTP. A mixture of mRNA from all 39 cell lines was prepared as the control probe. The mRNA mixture was amplified by T7-based amplification and reverse transcribed in the presence of Cy3-labeled dCTP, as described previously (5, 6, 7)
. Sample and control-labeled probes were mixed together and hybridized to cDNA microarray slides that contained 9216 human genes (3)
. Hybridized slides were scanned, and the fluorescence intensities of Cy5 (sample) and Cy3 (control) for each gene spot quantified by using a GenIII microarray scanner (Amersham Pharmacia Biotech) with Array Vision software (Imaging Research, Inc.). Each slide contained 52 housekeeping genes to normalize the signal intensities of the fluorescent dyes. The intensities of Cy5 and Cy3 were adjusted so that the mean Cy5 and Cy3 intensities of probe cDNA binding to the housekeeping genes were equivalent.
Classification of Anticancer Drugs According to Drug Activity Pattern Against Human Cancer Cell Lines.
A hierarchical clustering method was applied to the chemosensitivity data. Before applying the clustering algorithm, the GI50 values were log transformed, and the absolute values (|log10GI50|) were used for correlation analysis. Correlations were assessed by the Pearson correlation coefficient as described previously (2
, 8)
. Analyses were performed using "Gene Spring" software (Silicon Genetics).
Classification of Human Cancer Cell Lines According to Gene Expression Profiles.
A hierarchical clustering method was applied to the gene expression data. To obtain reproducible clusters, we selected only the 991 genes that passed the cutoff filter [control Cy3 signals were >50,000 relative fluorescent units in >80% of cases, i.e., 32 of 39 cell lines examined; the fluorescence ratio (Cy5:Cy3) was greater than the mean ratio by at least 5-fold in at least 3 of 39 cell lines]. Analysis was performed using Gene Spring. Before applying the clustering algorithm, the fluorescence ratio for each spot was log transformed (log2Cy5:Cy3), and the data for each sample centered to remove experimental biases. Correlations were assessed by the Pearson correlation coefficient as described (9
, 10)
.
Correlation Analysis between Gene Expression Profiles and Chemosensitivity Profiles.
We calculated the degree of similarity between drug activity and gene expression pattern using the Pearson correlation coefficient by the following formula:
![]() |
where xi represented the log expression ratio (log2Cy5:Cy3) of gene x in cell i, whereas yi was the log sensitivity (|log10GI50|) of cell i to drug y. xm represented the mean of the log expression ratio of gene x, and ym represented the mean sensitivity (|log10GI50|) of the drug. In this analysis, we selected genes that passed the cutoff filter (signal intensities were >25,000 relative fluorescent units or signal:noise ratios were >3 in either Cy3 or Cy5 in >80% of cases, i.e., 32 of 39 cell lines examined). We then selected genes with expression patterns that showed significant correlation to drug activity patterns. A significant correlation was defined as a P < 0.05 and a slope of the regression line >1.5, where the difference of the |log10GI50| values between the most and the least sensitive cell lines was fixed as 1.
Clustering Analysis of Drug and Gene Profiles on the Basis of Pearson Correlation Coefficients between Drugs and Genes.
We performed clustering analysis on the drugs and genes on the basis of the Pearson correlation coefficients, as described previously (10)
. We used all 1071 genes that passed the selection criteria described in the previous section for at least 1 of 55 drugs. The algorithm used to calculate the correlation between drug a and drug b was the standard correlation coefficient by the following formula:
![]() |
where ai represented the Pearson correlation coefficient between drug a and gene i, whereas bi was the coefficient between drug b and gene i.
Similarly, the algorithm used to calculate the correlation between gene c and gene d was the standard correlation coefficient by the following formula:
![]() |
where cj represented the Pearson correlation coefficient between gene c and drug j, whereas dj was the coefficient between gene d and drug j.
| RESULTS |
|---|
|
|
|---|
|
|
|
10 of the 55 drugs and found 50 genes that were significantly correlated to sensitivity to a relatively wide range of anticancer drugs (Fig. 3)
|
|
|
(SFN), LIM domain kinase 2 (LIMK2), and cathepsin H (CTSH) commonly showed negative correlations, whereas genes encoding DDB2 (DDB2) and the ALL-fused gene from chromosome1q (AF1Q) revealed positive correlations (Fig. 5, AC)| DISCUSSION |
|---|
|
|
|---|
Our database, comprising both gene expression data and drug activity data for the same set of cell lines, was first developed by Ross et al. and Scherf et al. (9
, 10)
using 60 human cancer cell lines (NCI60). These studies examined gene expression profiles in 60 cell lines that were pharmacologically characterized on treatment with various kinds of anticancer drugs. The number of cell lines used in the present study was 39. Of these, 9 cell lines (6 derived from gastric cancers and 3 from breast cancers) were not included in the NCI60 panel (4)
. Using this different set of 39 cell lines, we performed hierarchical clustering with respect to 55 anticancer drugs on the basis of growth inhibition and Pearson correlation coefficients between drugs and genes. As shown in Figs. 1
and 4
, drugs with similar mechanisms of action clustered most of the time, in agreement with results of an earlier study that used a different set of cell lines (10)
. Moreover, we have shown previously that the mode of action of a novel compound can be predicted by examining its growth inhibitory activities against the panel of 39 cell lines (2
, 8) . This indicated that our cell line panel had sufficient power to classify drugs according to mode of action. The previous study by Scherf et al. (10)
showed expected correlations, such as that between 5-fluorouracil and dihydropyrimidine dehydrogenasae (DPYD). A moderately negative correlation between 5-fluorouracil and TYMS was demonstrated in our study (r = -0.29, P < 0.1). In addition, a stronger negative correlation between 5-fluorouracil derivative carmofur and TYMS was also observed (r = -0.38, P < 0.01; data not shown). Furthermore, positive correlations between the DNA topo II inhibitor etoposide and DNA topo II
(TOP2A, r = 0.39, P < 0.05) and ß (TOP2B, r = 0.29, P < 0.1) were also found (data not shown). Although the Pearson correlation coefficients were not high, together these moderate correlations might be sufficient as the sensitivity of cancer cells to anticancer drugs appears not to be determined by the expression of a single gene.
Clustering analysis of cell lines revealed that overall profiles of gene expression in cancer cells reflected their tissues of origin (Fig. 2)
. Therefore, some of the retrieved genes preferentially expressed in cells derived from specific tissues, e.g., keratin genes related to epithelial cell type (data not shown). Similarly, some of the drugs used in our study tended to show tissue-oriented efficacy (Fig. 1)
. Therefore, some of the genes reflecting cells tissue of origin correlating with chemosensitivity might be associated with sensitivity to such drugs. On that occasion, the correlation should be validated within each set of cell lines derived from the same tissues. On the other hand, there was a substantial population of genes that was not strongly associated with tissue of origin, including AKR1B1, LIMK2, and DDB2 (data not shown). These genes are conceivably good candidates for predictive markers of drug efficacy.
Clustering analysis of the drugs on the basis of Pearson correlation coefficients between drug activity and gene expression levels revealed that drugs with similar modes of action clustered into the same groups and that the sets of chemosensitivity-associated genes were similar between each drug within a group (Fig. 4)
. Thus, we were able to retrieve genes that commonly correlated within each of the groups of drugs with similar functions (Fig. 5, AD
, data not shown), e.g., we identified the genes encoding survivin (BIRC5) and C-IAP1 of apoptosis1 (BIRC2) that revealed negative correlations with 5-fluorouracil derivatives, although survivin showed either no significant negative correlations or moderately positive correlations to the other genotoxic anticancer drugs (Fig. 5D)
. The inhibitor of apoptosis family of proteins, including survivin, has been shown to play important roles in the suppression of apoptosis and is conceivably a chemoresistance factor (15
, 16)
. Thus, the most important aspect of our analysis is that we were able to retrieve different gene sets that may be used as putative diagnostic markers for chemosensitivity prediction with respect to each class of drugs with similar mechanisms of action. The next step will be to confirm the predictive power of our findings. In a recently published study, Staunton et al. (17)
identified putative predictive markers of chemosensitivity and showed the feasibility of chemosensitivity prediction by transcriptional profiling. In our current study, using a different methodology, we could narrow down the number of candidate genes from thousands to dozens that predict chemosensitivity, although we have not completed to evaluate our findings. The validation studies using dozens of cell lines outside the panel are under way.
In summary, we have identified different sets of genes that may act as predictive markers for chemosensitivity to drugs with similar mechanisms of action. Data-mining methodologies described here and elsewhere (17) may be useful to develop systems to select suitable chemotherapies in the future. However, statistical correlations are not the proof of causal relationships between gene expression and chemosensitivity. If these causal relationships can be demonstrated experimentally, the resultant gene products may explain the mechanisms of anticancer drugs. In addition, to elucidate the mechanisms of anticancer drugs, it is also necessary to determine gene expression changes after administration of drugs (18) . These studies will help us to discover new targets of cancer chemotherapy, as well as predictive markers for conventional drugs.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
1 Supported by Grant-in-Aid for Scientific Research on Priority Areas (C) from the Ministry of Education, Culture, Sports, Science and Technology, Japan. ![]()
2 To whom requests for reprints should be addressed, at Division of Molecular Pharmacology, Cancer Chemotherapy Center, Japanese Foundation for Cancer Research, 1-37-1 Kami-Ikebukuro, Toshima-ku Tokyo 170-8455, Japan. Phone: 81-3-5394-4068; Fax: 81-3-3918-3716; E-mail: yamori{at}ims.u-tokyo.ac.jp ![]()
3 The abbreviations used are: NCI, National Cancer Institute; GI50, drug concentration required for 50% growth inhibition; topo, DNA topoisomerase; TYMS, thymidylate synthetase. ![]()
Received 7/ 2/01. Accepted 12/10/01.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. Lynch, Y. Fukuda, P. Krishnamurthy, G. Du, and J. D. Schuetz Cell Survival under Stress Is Enhanced by a Mitochondrial ATP-Binding Cassette Transporter That Regulates Hemoproteins Cancer Res., July 1, 2009; 69(13): 5560 - 5567. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Shichiri, N. Fukai, Y. Kono, and Y. Tanaka Rifampicin as an Oral Angiogenesis Inhibitor Targeting Hepatic Cancers Cancer Res., June 1, 2009; 69(11): 4760 - 4768. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Takahata, Y. Inoue, H. Tsuda, I. Imoto, D. Koinuma, M. Hayashi, T. Ichikura, T. Yamori, K. Nagasaki, M. Yoshida, et al. SKI and MEL1 Cooperate to Inhibit Transforming Growth Factor-{beta} Signal in Gastric Cancer Cells J. Biol. Chem., January 30, 2009; 284(5): 3334 - 3344. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Nakatsu, T. Nakamura, K. Yamazaki, S. Sadahiro, H. Makuuchi, J. Kanno, and T. Yamori Evaluation of Action Mechanisms of Toxic Chemicals Using JFCR39, a Panel of Human Cancer Cell Lines Mol. Pharmacol., November 1, 2007; 72(5): 1171 - 1180. [Abstract] [Full Text] [PDF] |
||||
![]() |
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] |
||||
![]() |
A. Jimeno, B. Rubio-Viqueira, M. L. Amador, V. Grunwald, A. Maitra, C. Iacobuzio-Donahue, and M. Hidalgo Dual mitogen-activated protein kinase and epidermal growth factor receptor inhibition in biliary and pancreatic cancer Mol. Cancer Ther., March 1, 2007; 6(3): 1079 - 1088. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Fryknas, L. Rickardson, M. Wickstrom, S. Dhar, H. Lovborg, J. Gullbo, P. Nygren, M. G. Gustafsson, A. Isaksson, and R. Larsson Phenotype-Based Screening of Mechanistically Annotated Compounds in Combination with Gene Expression and Pathway Analysis Identifies Candidate Drug Targets in a Human Squamous Carcinoma Cell Model J Biomol Screen, August 1, 2006; 11(5): 457 - 468. [Abstract] [PDF] |
||||
![]() |
S.-i. Yaguchi, Y. Fukui, I. Koshimizu, H. Yoshimi, T. Matsuno, H. Gouda, S. Hirono, K. Yamazaki, and T. Yamori Antitumor activity of ZSTK474, a new phosphatidylinositol 3-kinase inhibitor. J Natl Cancer Inst, April 19, 2006; 98(8): 545 - 556. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. D. Petty, K. M. Kerr, G. I. Murray, M. C. Nicolson, P. H. Rooney, D. Bissett, and E. S.R. Collie-Duguid Tumor Transcriptome Reveals the Predictive and Prognostic Impact of Lysosomal Protease Inhibitors in Non-Small-Cell Lung Cancer J. Clin. Oncol., April 10, 2006; 24(11): 1729 - 1744. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Wallqvist, J. Connelly, E. A. Sausville, D. G. Covell, and A. Monks Differential Gene Expression as a Potential Classifier of 2-(4-Amino-3-methylphenyl)-5-fluorobenzothiazole-Sensitive and -Insensitive Cell Lines Mol. Pharmacol., March 1, 2006; 69(3): 737 - 748. [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] |
||||
![]() |
K. Iwao-Koizumi, R. Matoba, N. Ueno, S. J. Kim, A. Ando, Y. Miyoshi, E. Maeda, S. Noguchi, and K. Kato Prediction of Docetaxel Response in Human Breast Cancer by Gene Expression Profiling J. Clin. Oncol., January 20, 2005; 23(3): 422 - 431. [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] |
||||
![]() |
H.-R. Park, A. Tomida, S. Sato, Y. Tsukumo, J. Yun, T. Yamori, Y. Hayakawa, T. Tsuruo, and K. Shin-ya Effect on Tumor Cells of Blocking Survival Response to Glucose Deprivation J Natl Cancer Inst, September 1, 2004; 96(17): 1300 - 1310. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Wu, P. de Kievit, L. Vahlkamp, D. Pijnenburg, M. Smit, M. Dankers, D. Melchers, M. Stax, P. J. Boender, C. Ingham, et al. Quantitative assessment of a novel flow-through porous microarray for the rapid analysis of gene expression profiles Nucleic Acids Res., August 27, 2004; 32(15): e123 - e123. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. D. Petty, M. C. Nicolson, K. M. Kerr, E. Collie-Duguid, and G. I. Murray Gene Expression Profiling in Non-Small Cell Lung Cancer: From Molecular Mechanisms to Clinical Application Clin. Cancer Res., May 15, 2004; 10(10): 3237 - 3248. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. R. Bani, M. I. Nicoletti, N. W. Alkharouf, C. Ghilardi, D. Petersen, E. Erba, E. A. Sausville, E. T. Liu, and R. Giavazzi Gene expression correlating with response to paclitaxel in ovarian carcinoma xenografts Mol. Cancer Ther., February 1, 2004; 3(2): 111 - 121. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Vekris, D. Meynard, M.-C. Haaz, M. Bayssas, J. Bonnet, and J. Robert Molecular Determinants of the Cytotoxicity of Platinum Compounds: The Contribution of in Silico Research Cancer Res., January 1, 2004; 64(1): 356 - 362. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Mariadason, D. Arango, Q. Shi, A. J. Wilson, G. A. Corner, C. Nicholas, M. J. Aranes, M. Lesser, E. L. Schwartz, and L. H. Augenlicht Gene Expression Profiling-Based Prediction of Response of Colon Carcinoma Cells to 5-Fluorouracil and Camptothecin Cancer Res., December 15, 2003; 63(24): 8791 - 8812. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Sugiyama, S. Dan, Y. Yoshida, F. Akiyama, K. Sugiyama, Y. Hirai, M. Matsuura, S. Miyata, M. Ushijima, K. Hasumi, et al. A Large-Scale Gene Expression Comparison of Microdissected, Small-Sized Endometrial Cancers with or without Hyperplasia Matched to Same-Patient Normal Tissue Clin. Cancer Res., November 15, 2003; 9(15): 5589 - 5600. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. E. Lamendola, Z. Duan, R. Z. Yusuf, and M. V. Seiden Molecular Description of Evolving Paclitaxel Resistance in the SKOV-3 Human Ovarian Carcinoma Cell Line Cancer Res., May 1, 2003; 63(9): 2200 - 2205. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Ding and C. R. Cantor A high-throughput gene expression analysis technique using competitive PCR and matrix-assisted laser desorption ionization time-of-flight MS PNAS, March 18, 2003; 100(6): 3059 - 3064. [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 |