Genome-wide analysis of three-way interplay among gene expression, cancer cell invasion and anti-cancer compound sensitivity
- Yi-Chiung Hsu1, 2,
- Hsuan-Yu Chen1,
- Shinsheng Yuan1,
- Sung-Liang Yu3,
- Chia-Hung Lin1,
- Guani Wu1,
- Pan-Chyr Yang†2, 4, 5 and
- Ker-Chau Li†1, 6Email author
© Hsu et al.; licensee BioMed Central Ltd. 2013
Received: 22 May 2012
Accepted: 11 December 2012
Published: 16 April 2013
Chemosensitivity and tumor metastasis are two primary issues in cancer management. Cancer cells often exhibit a wide range of sensitivity to anti-cancer compounds. To gain insight on the genetic mechanism of drug sensitivity, one powerful approach is to employ the panel of 60 human cancer cell lines developed by the National Cancer Institute (NCI). Cancer cells also show a broad range of invasion ability. However, a genome-wide portrait on the contributing molecular factors to invasion heterogeneity is lacking.
Our lab performed an invasion assay on the NCI-60 panel. We identified invasion-associated (IA) genes by correlating our invasion profiling data with the Affymetrix gene expression data on NCI-60. We then employed the recently released chemosensitivity data of 99 anti-cancer drugs of known mechanism to investigate the gene-drug correlation, focusing on the IA genes. Afterwards, we collected data from four independent drug-testing experiments to validate our findings on compound response prediction. Finally, we obtained published clinical and molecular data from two recent adjuvant chemotherapy cohorts, one on lung cancer and one on breast cancer, to test the performance of our gene signature for patient outcome prediction.
First, we found 633 IA genes from the invasion-gene expression correlation study. Then, for each of the 99 drugs, we obtained a subset of IA genes whose expression levels correlated with drug-sensitivity profiles. We identified a set of eight genes (EGFR, ITGA3, MYLK, RAI14, AHNAK, GLS, IL32 and NNMT) showing significant gene-drug correlation with paclitaxel, docetaxel, erlotinib, everolimus and dasatinib. This eight-gene signature (derived from NCI-60) for chemosensitivity prediction was validated by a total of 107 independent drug tests on 78 tumor cell lines, most of which were outside of the NCI-60 panel. The eight-gene signature predicted relapse-free survival for the lung and breast cancer patients (log-rank P = 0.0263; 0.00021). Multivariate Cox regression yielded a hazard ratio of our signature of 5.33 (95% CI = 1.76 to 16.1) and 1.81 (95% CI = 1.19 to 2.76) respectively. The eight-gene signature features the cancer hallmark epidermal growth factor receptor (EGFR) and genes involved in cell adhesion, migration, invasion, tumor growth and progression.
Our study sheds light on the intricate three-way interplay among gene expression, invasion and compound-sensitivity. We report the finding of a unique signature that predicts chemotherapy survival for both lung and breast cancer. Augmenting the NCI-60 model with in vitro characterization of important phenotype-like invasion potential is a cost-effective approach to power the genomic chemosensitivity analysis.
KeywordsNCI-60 Invasion Metastasis Microarray Chemotherapy
Gene-expression profiling of tumors from patient cohorts has been used to develop gene signatures for clinical outcome prediction. Recently, a signature combining estrogen receptor (ER) status and predicted chemo/endocrine responsiveness succeeded in identifying patients with high probability of survival following taxane and anthracycline chemotherapy . However, the biological mechanism of the genes involved in such cohort-initiated genomic predictors may not be easy to elucidate. On the other hand, in vitro chemosensitivity experiments on cancer cell lines, such as the NCI-60 cell line panel, are helpful in elucidating the complex relationship between drug responsiveness and gene expression. Despite this, the more challenging problem of how to translate the elucidated relationship for clinical outcome prediction still awaits more studies.
In addition to chemosensitivity, metastasis is another major issue in studying treatment efficacy for many cancers, including invasive breast and lung cancer. Like the drug-sensitivity heterogeneity, tumor cells often exhibit a wide range of invasion ability. Such invasion heterogeneity may exist not only between different cancer types, but also among the individual cells from the same malignant neoplasm of a patient. More subtly, it is possible that the clinical outcome of chemotherapy may hinge on the growth-inhibition of the more invasive cells rather than the less invasive cells in the neoplasm. However, a characterization of genes associated with both the invasion potential and drug-sensitivity heterogeneity is lacking.
Molecular markers of tumor invasion potential: to identify the set of IA genes whose expression levels are likely to be indicative of the invasion potential of a tumor;
Drug sensitivity prediction by tumor-invasion markers: to evaluate how the expression levels of IA genes in a tumor are likely to be indicative of the tumor’s resistance or responsiveness to an anti-cancer drug;
Drug discovery with predictable sensitivity: to find anti-cancer drugs whose efficacies correlate with tumor-invasion potential and can be predicted by tumor-invasion markers;
Clinical validation: to demonstrate the use of the IA gene signature for predicting clinical outcome.
NCI-60 is a diverse panel of 60 cell lines used by the Development Therapeutics Program (DTP) of the National Cancer Institute (NCI) to screen more than 100,000 compounds since 1990 [2–5]. These human cancer cell lines are derived from patients with leukemia, melanoma, lung cancer, colon cancer, central nervous system cancer, ovarian cancer, renal cancer, breast cancer and prostate cancers. The molecular characteristics of these cell lines have been subjected to various DNA microarray studies using both Affymetrix (Santa Clara, CA, USA) and spotted cDNA array technology .
Both drug sensitivity and gene expression profiles of the NCI-60 panel are available from the public domain; however, there is no public invasion phenotype data for NCI-60. So we conducted the invasion assay for 53 solid tumor cell lines from the NCI-60 panel in our lab.
Matrigel invasion assay
We purchased the panel of NCI human cancer cell lines (NCI-60) to conduct tumor invasion assays. The suspension leukemia cancer cell lines were excluded. All cells were grown in tissue culture flasks at 37°C in 5% CO2 in RPMI 1640 with 2 mM L-glutamine, and 10% fetal bovine serum, all from Invitrogen, Eugene, OR. The invasion capacities for cell lines were examined by using membrane invasion culture system. The HTS FluoroBlok inserts containing 8-μm pores (Falcon, Becton Dickinson, Franklin Lakes, NJ) were coated 30 μg matrigel (BD Biosciences, San Jose, CA, USA). The cells were suspended in RPMI containing 10% FBS and seeded into the upper wells of the chamber (2.5 × 104 cells/well). After incubating for 24 hours at 37°C, the membrane of the transwell was fixed for 10 minutes at room temperature with methanol and stained for 30 minutes with 50 μg/ml propidium iodide (Sigma, St. Louis, MO). The number of cells in each blot was counted under a microscope with Analytical Imaging Station system (Imaging Research Inc., St. Catharines, ON, Canada). For each cell line, we reported the average invasion cell count after three repeats (n = 3).
Cell line gene expression data
Two gene expression datasets for NCI-60 cell lines, produced by Gene Logic (Gaithersburg, MD, USA) using Affymetrix U95, and U133plus2, respectively, were downloaded from the DTP . Probe mapping between U133plus2 and U95 was provided by Affymetrix.
The independent gene expression data of 78 cell lines used in the drug response validation were obtained either from the GEO website, GSE6569 (n = 23, breast cancer) , GSE9633 (n = 16, prostate cancer) , and GSE4127 (n = 29, lung cancer)  or by e-mail request (n = 10, lung, Balko et al.) . All but nine cell lines were outside of the NCI-60 panel (Additional file 1: Table S1).
We downloaded the NCI-60 chemosensitivity data for anti-cancer drugs with known molecular mechanisms, including those used in targeted therapy . Drug sensitivity was measured by the negative of log10 GI50 (−logGI50). Seven compounds, inactive in all cell lines, were excluded. A total of 99 drugs were analyzed.
The real-time reverse transcription quantitative PCR analysis
The expressions of eight signature genes and a control gene TBP in nine lung cancer cell lines of NCI-60 were measured by reverse transcription qPCR(RT-qPCR) with specific Taqman probes and primer sets (Additional file 1: Table S2). The transcripts were amplified with Taqman One-Step RT-PCR Master Mix Reagent (Applied Biosystems) and a detection system (ABI Prism 7900HT, Applied Biosystems). Gene expression was quantified in relation to the expression of the control gene with the use of sequence detector software and the relative quantification method (Applied Biosystems).
The comparison of the tumor invasion ability between tissue types was performed by ANOVA. We applied the normal score transformation to preprocess the gene expression data before computing the Pearson correlations. The invasion-associated genes were obtained by computing the Pearson's correlation coefficients between gene expression profiles and invasion ability profile in NCI-60 cell lines. Clustering was used to order the IA probes. Student’s t-test was used in validating the drug sensitivity prediction of gene signatures. All statistical analyses were performed in the R language environment .
Cancer chemotherapy cohorts
Two chemotherapy cohorts were used in this study. The lung cancer cohort came from JBR.10, a randomized controlled trial of adjuvant vinorelbine/cisplatin versus observation alone . We used the adjuvant cisplatin/vinorelbine therapy arm where after surgery (n = 71). The breast cancer cohort of 508 patients came from a prospective multicenter study, conducted at the M. D. Anderson Cancer Center, of which the patients were those with newly diagnosed Human Epidermal Growth Factor Receptor 2 (HER2)–negative breast cancer under chemotherapy containing sequential taxane and anthracycline–based regimens (followed by endocrine therapy if ER positive) .
Four cohorts where the patients were systemically untreated after surgery were used. Three of them were the breast cancer cohorts (GSE2034 n = 286 , GSE7390 n = 198  and GSE11121 n = 200 ). The fourth, a lung cancer cohort, was the OBS arm where patients were under observation after surgery (n = 62).
Kaplan-Meier survival curves were obtained and compared by log-rank tests. Multivariate Cox proportional hazard regression analysis was used to evaluate independent prognostic factors, such as age, gender, tumor stage nodal status and histological grade.
Identification of 633 IA genes
We identified the IA genes from Affymetrix U133plus2 NCI-60, and U95 NCI-60 microarray gene expression datasets by a two-stage procedure which required that (a) the gene expression profile must be significantly correlated with the invasion profile in both datasets, and (b) the sign of correlation must be consistent in both datasets (see the Additional file 1, Supplementary information). A total of 744 probes, which represented 633 distinct genes, were obtained. We estimated the false discovery rate (FDR) at the confirmation stage to be 0.08 (FDR = (2,417*0.025)/744; where 2,417 = total number of probes considered at the confirmation stage and 0.025 = the P-value cutoff with the sign consistency criterion).
Functional enrichment analysis of IA genes
We use the MetaCore™ (Thomson Reuters, New York, NY), a web-based computational platform designed for system biology and drug discovery, to conduct functional enrichment analysis of IA genes. We input the set of 633 IA genes (744 probes) and used the Functional Ontology Enrichment tool with the default settings. The results show that our IA genes were enriched in cell adhesion and cytoskeleton remodeling pathways (Additional file 1: Table S3).
Additional file 1: Table S3 list the significantly enriched pathways and networks (P <0.0001). The results indicate the involvement of cell adhesion and cytoskeleton remolding actin cytoskeleton network (ACTN1, ACTN4, ATCB, ACTG1, ZYX, VCL and CFL), integrin signaling (ITGA3, ITGB1, CAV1 and CAV2), matrix metalloproteinase signaling (CDH2, CD44, TGFB2, TGFBR2, JAK1, SMAD1 and SMAD3), microtubule cytoskeleton (TUBB and TUBB6), and myosin signaling (MYLK1, LIMK1, MYL6 and MYH9). We further examined the Integrin-mediated cell adhesion and migration pathway, (P = 5.49E-8, input/total nodes = 11/48), which has a much higher ratio of root to total nodes and significant P-value. The Integrin-mediated cell adhesion and migration pathway shows that ITGA3 promotes focal adhesion kinase (FAK) autophosphorylation and creates a binding site for c-Src. EGFR signaling also activates the FAK/Src pathway [19, 20]. FAK activation regulates the ERK-mediated phosphorylation and activation of Myosin light chain kinase (MYLK) contributes to cell-matrix adhesion dynamics . Integrin recruits FAK and a cytoskeletal protein vinculin and alpha-actin to focal contacts. c-Src and ERK2-mediated phosphorylation of FAK1 promotes its release from focal contacts and ERK2-mediated phosphorylation of paxillin promotes the association of non-phosphorylated FAK1 with paxillin at new or growing focal contact sites . Activation of MYLK together with inactivation of PAK1 contributes to cell-matrix adhesion dynamics. PAK1 phosphorylation leads to the activation of LIM-kinase 1 (LIMK1) , inhibition of MYLK, activation of myosin regulatory light chains (MRLC) . FAK is a tyrosine kinase which interacts with the important oncogene c-Src. FAK signaling is important for integrin regulated cell adhesion and migration. Notably, this pathway is the target of dasatinib , a drug we address further in this paper. A simplified version of the Integrin-mediated cell adhesion and migration pathway is shown Additional file 1: Figure S2, which features three IA genes, EGFR, ITGA3, MYLK.
Association of drug sensitivity with IA gene expression
To evaluate how the abundance of the IA gene transcripts in a cancer cell correlates with the cell’s response to each of the 99 anti-cancer drugs for which the drug-sensitivity profiles on NCI-60 were available, we computed Pearson's correlation between the gene expression profile of each IA gene and the chemosensitivity profile of each compound. A positive correlation means that cell lines with higher gene expression are more sensitive to the drug, while a negative correlation indicates the opposite. The gene-drug correlation results are shown in a heat map (Additional file 1: Figure S3A). We identified all the significant correlations (P <0.05), and displayed the findings in Additional file 1: Figure S3B. An IA probe showing a significant correlation with a drug would be called a drug-sensitivity-correlated probe for that drug. We counted the number of drug-sensitivity-correlated IA probes for each compound.
A higher count indicates the availability of more IA genes for the sensitivity prediction of that drug. Such a drug would be more likely to have differential anti-cancer effects among the tumors showing differential invasion potential. Drugs with higher counts are this study’s primary interest.
Individually speaking, the compound with the highest count is zoledronic acid, a bisphosphonate drug used to prevent skeletal fractures in cancer patients and to treat osteoporosis. Interestingly, recent studies indicate zoledronic acid can prevent skeletal metastases through inhibition of invasion and angiogenesis of cancer cells [26, 27]. However, we were unable to find the drug sensitivity validation and gene expression data outside the NCI panel for zoledronic acid.
We find that on average, the group of tubulin-binding agents ranks the highest, followed by the targeted therapy. The count for each compound in these two groups is given in Figure 4B. We continued our study on the 17 drugs from these two groups.
Gene-drug heat map for 17 drugs
Dasatinib is a clinically studied SRC inhibitor for cancer therapy . Our functional enrichment analysis of IA genes shows that the Src signaling pathway played an important role in integrin-mediated cell adhesion and migration pathway (Additional file 1: Figure S2). Erlotinib inhibits the intracellular phosphorylation of tyrosine kinase associated with the EGFR . Everolimus is an mTOR inhibitor . In the anti-tubulin group, docetaxel is a semisynthetic side chain analog of paclitaxel.
Selection of drug-sensitivity-correlated IA genes for five compounds
In addition to the mechanism sharing and the clustering pattern observed, we also noticed that the numbers of sensitivity-correlated IA genes that these five compounds had are higher than other drugs in the respective mechanism groups (Figure 4A).
This prompted us to ask if these five compounds might share any drug-sensitivity-correlated IA genes that could be used for drug-sensitivity prediction. By inspecting the drug-sensitivity-correlated IA gene list for these five compounds, 26 common probes are identified (Additional file 1: Table S5). Among them, 19 probes are in the Affy HG U-133A chip and 7 probes are in the HG U-133B chip. Because most microarray data are available in the public domain used only the U-133A chip, we exclude the seven U133-B probes. We use the standard deviation (SD) of the expression across all 53 cell lines to further exclude probes showing low expression variation. Only those probes with SD ranking within the top 10% are retained. Finally, we check the sign of correlation between gene expression and drug sensitivity to ensure the consistency in predicting sensitivity or resistance.
A final set of eight IA genes (EGFR, AHNAK, GLS, IL32, ITGA3, MYLK, NNMT, and RAI14) are obtained. We confirm the quality of the microarray gene expression data by performing the qPCR assay for these eight genes on nine lung cancer cell lines in NCI-60 (Additional file 1: Figure S4). A total of 72 data points are plotted, each dot representing the expression level of a gene measured by qPCR (horizontal axis) and the microarray (vertical axis). A significant positive trend (Pearson correlation = 0.69, 2.129e-11) is observed, confirming the consistency of the two assays.
Biological functions of eight invasion-associated genes
The eight IA genes we identified are associated with cellular cytoskeleton, cell invasion and oncogenetic signaling, including the well-known cancer driver EGFR, and FAK-Src signaling (ITGA3, MYLK); ITGA3 (integrin, alpha 3) mediates cell survival and invasion through FAK-Src signaling [31–33]. MYLK is a key target of FAK–Src signaling . Glutaminase (GLS) was shown to be up-regulated in MYC induction cancer cells and inhibition of GLS decreased the cancer cell growth [35, 36]. Retinoic acid induced 14 (RAI14) is an actin cytoskeleton protein regulated by retinoic acid . High expression of interleukin 32 (IL32) was shown to cause a worse clinical outcome in lung cancer patients . AHNAK nucleoprotein was reported as pseudopod-specific proteins in different metastatic human tumor cell lines . The activity of nicotinamide N-methyltransferase (NNMT) in catalyzing the N-methylation of nicotinamide is important for biotransformation of drug and xenobiotic compounds . NNMT was identified as a tumor biomarker [41–43], promoting cell migration  and cell invasion .
Validating gene signature with independent cell lines
In these validation experiments, we found nine cell lines from the NCI-60 panels. We removed these cell lines from the validation data and conducted the comparison again (Additional file 1: Figure S5). All significant differences in the first three experiments were still valid. For the fourth experiment, the P-value increased to 0.0788, mainly because the sample size (n = 8) was too small.
Clinical outcome prediction in adjuvant chemotherapy lung and breast cancer patients
To test if the eight genes have the predictive power in clinical outcome evaluation, we searched the public domain for chemotherapy cohort studies involving any of the two anti-microtubule compounds we used in obtaining the eight genes. Two recently published studies on the vinorelbine containing regimen for lung cancer  and the taxane-containing regimen for treating invasive breast cancer  were obtained. We use the simple average expression of the eight-genes as the risk score to dichotomize patients into two groups.
Multivariate Cox regression analysis of the eight-gene signature for predicting relapse-free survival in cancer patients
Lung cancer (n = 71)
1.76 to 16.1
1.00 to 1.11
Gender (Male vs Female)
0.48 to 3.68
Stage (1 vs 2)
1.12 to 6.45
0.56 to 3.73
Breast cancer (n = 462)
1.19 to 2.76
Age (>50 vs ≦50)
0.70 to 1.56
Clinical nodal status (positive vs negative)
1.45 to 4.18
Clinical tumor stage (T3 or T4 vs T1 or T2)
1.20 to 2.70
ER status (IHC positive vs negative)
0.29 to 0.67
The invasive or metastatic potential of a malignant neoplasm and the growth-inhibition of tumor cells by a therapeutic agent are two common denominators of patient survival in cancer systemic therapy. While cell line models have been used to predict treatment response or patient survival by genes associated with drug sensitivity in the cancer cell lines [46–49], these studies have not considered the varying metastasis potential of a tumor. To weigh in the interplay of both factors directly, our pre-clinical gene signature discovery method features the co-integration of invasion phenotypes and compound-sensitity profiles with gene expression at the full genome scale.
Our approach is based on the observation that invasion ability and drug sensitivity are both phenotypes of the cell lines available for study. Each phenotype is naturally associated with its own set of molecular determinants. We hypothesize the potential overlap between the set for invasion ability phenotype and the set for drug responsiveness phenotype. By identifying these common determinants, then we may use these shared determinants to estimate the overall invasion potential of the cancer cells in a tumor and also use it to predict the drug response at the same time. However, because the tumor microenvironment in a patient is different from the growth environment of cancer cell lines monitored in a lab, the robustness of an invasion molecular marker becomes an important factor for increasing the chance of success in clinical applications.
An alternative strategy of analyzing the three-way interaction of invasion, gene expression and drug response would be to correlate invasion with drug response first. Once the most correlated drugs are identified, then genes correlated with response to these drugs can be used to predict the drug sensitivity. However, we did not pursue this line of analysis further because the phenotye-phenotype correlation is often weaker than phenotype-genomic determinant correlation. As a matter of fact, our data show that most drug-invasion correlations appear weak; only four drugs pass the statistical significance and the best two correlations are only 0.39 and −0.35 (Additional file 1: Table S6). Our approach overcomes the limitations of weak phenotype-phenotype correlation by looking for statistical evidence of correlations directly from the genomic determinant. This helps improve the robustness of the genetic marker thus obtained.
Previously, without considering drug sensitivity, our team performed invasion profiling for the nine lung cancer cell lines of NCI-60 to obtain a four gene signature for clinical outcome prediction . We find that among the four genes ANKRD49 and LPHN1 are in the IA gene list and only ANKRD49 has a significant correlation with paclitaxel and docetaxel. The four-gene signature failed to predict the survival outcome for the two validation cohorts receiving adjuvant chemotherapy (Additional file 1: Figure S7).
To gain robustness of our gene signature, instead of using different panels of tissue origins in NCI-60 to obtain different sets of IA genes for different types of cancer, we used the invasion data from all 53 NCI60 solid tumor cell lines and obtained 633 IA genes. Then a series of statistical analyses were designed to increase the robustness of the final eight-gene signature in predicting drug sensitivity for the selected compounds. The eight-gene score differed between drug-resistant and drug-sensitive cell lines (Figure 6). We succeeded in applying our gene signature to one lung cancer cohort and one breast cancer cohort, of which the patients received a regimen containing an anti-microtubule agent. The success in using the same signature to predict patient outcome for different types of cancer showed the robustness of this gene signature.
The eight-gene signature showed a positive correlation with the sensitivity of targeted therapy compounds and a negative correlation with the sensitivity of anti-microtubule compounds (Figure 5B). Because high values of the eight-gene signature correlate with high invasion potential in cancer cells, this suggests that the direct correlation between the invasion profile and the sensitivity profiles of anti-MT drugs may be negative. This is indeed the case, but the correlation is weak (Additional file 1: Table S6, correlation = −0.16, -0.28 for paclitaxel and docetaxel, respectively). On the other hand, the correlation between the eight-gene score and the drug sensitivity is stronger (−0.41, -0.54, respectively). Similarly, the correlation between the eight-gene score and the sensitivity for erlotinib, dasatinib and everolimus is 0.46, 0.52, 0.44, respectively, which is again stronger than the correlation between invasion and drug sensitivity (0.26, 0.24, 0.09, respectively). Therefore, despite the weak correlation between the invasion phenotype and the drug sensitivity phenotype, the eight-gene signature is an effective genomic marker for invasion potential and it can be used to predict the drug’s differential growth-inhibition efficacy that varies between cancer cells of higher invasion potential and those of lower invasion potential.
When applying to the patient’s tumor specimen, the eight-gene signature provides an averaged profile of the gene expression by individual cells with varying invasion potential. A low eight-gene score indicates that the overall invasion potential of the tumor is low and the chance of the patient’s favorable response to regimens containing anti-microtubule compounds increases. On the other hand, a high eight-gene score predicts the abundance of the cells of higher invasion potential, which are harder to eradicate by anti-microtubule compounds, but may be more likely to succumb to the said targeted therapy. This suggests the combined use of targeted therapy like dasatinib or erlotinib with anti-microtubule agents to increase the regimen efficacy of chemotherapy alone. There have been several studies on augmenting the anticancer effect of chemotherapy with targeted therapy. Erlotinib was shown to be more sensitive in the doxorubicin-resistant human breast cancer cell lines and paclitaxel-resistant human ovarian cancer cell lines  and the sensitivity was positively correlated with EGFR expression. More references were provided in Additional file 1, Supplementary information Text II.
There is room to improve our eight-gene signature for drug-sensitivity prediction. The overlap in distribution between the drug-sensitivity group and the drug-resistant group (Figure 6) suggests that drug response in cell lines is a very complex phenotype which is not fully characterized by our gene signature. Other genomic components, such as DNA copy number, single-nucleotide polymorphisms, methylation and microRNA, have not been considered in our study. In addition, differences in lab environment may also contribute to the variations observed in the data.
We have shown that augmenting the NCI-60 model with in vitro characterization of important phenotypes like invasion potential is a cost-effective approach to power the genomic chemosensitivity analysis. Our analysis delineates the complex three-way interplay of gene expression, cancer cell’s invasion potential and cancer cell’s responsiveness to an anti-cancer compound. We report the identification of a unique eight-gene signature for both lung and breast cancer, which predicted the relapse-free survival of adjuvant chemotherapy patients. The signature features the cancer hallmark EGFR and genes involved in cell adhesion cell migration, cell invasion, tumor growth and tumor progression. The discovery of prognostic biomarkers for chemotherapy patients remains critical toward improving the efficacy of cancer treatment. The eight-gene signature obtained here may be useful for the development of individualized cancer therapy. Our method of gene discovery may be applicable in studying other cancers.
Central nervous system
Development therapeutics program
Epidermal growth factor receptor
Focal adhesion kinase
False discovery rate
Invaded cell counts
Integrin, alpha 3
Myosin regulatory light chains
Myosin light chain kinase
National cancer institute
This research is supported by NSC 98-2314-B-001-001-MY3 MIB, Institute of Statistical Science, Academia Sinica, Taiwan, and AS-100-TP-AB2. We are thankful to three reviewers and Editor for constructive suggestions.
- Hatzis C, Pusztai L, Valero V, Booser DJ, Esserman L, Lluch A, Vidaurre T, Holmes F, Souchon E, Wang H, Martin M, Cotrina J, Gomez H, Hubbard R, Chacon JI, Ferrer-Lozano J, Dyer R, Buxton M, Gong Y, Wu Y, Ibrahim N, Andreopoulou E, Ueno NT, Hunt K, Yang W, Nazario A, DeMichele A, O'Shaughnessy J, Hortobagyi GN, Symmans WF: A genomic predictor of response and survival following taxane-anthracycline chemotherapy for invasive breast cancer. JAMA. 2011, 305: 1873-1881. 10.1001/jama.2011.593.View ArticlePubMedGoogle Scholar
- Shoemaker RH: The NCI60 human tumour cell line anticancer drug screen. Nat Rev Cancer. 2006, 6: 813-823.View ArticlePubMedGoogle Scholar
- Paull KD, Shoemaker RH, Hodes L, Monks A, Scudiero DA, Rubinstein L, Plowman J, Boyd MR: Display and analysis of patterns of differential activity of drugs against human tumor cell lines: development of mean graph and COMPARE algorithm. J Natl Cancer Inst. 1989, 81: 1088-1092. 10.1093/jnci/81.14.1088.View ArticlePubMedGoogle Scholar
- Grever MR, Schepartz SA, Chabner BA: The National Cancer Institute: cancer drug discovery and development program. Semin Oncol. 1992, 19: 622-638.PubMedGoogle Scholar
- Holbeck SL, Collins JM, Doroshow JH: Analysis of Food and Drug Administration-approved anticancer agents in the NCI60 panel of human tumor cell lines. Mol Cancer Ther. 2010, 9: 1451-1460. 10.1158/1535-7163.MCT-10-0106.View ArticlePubMedPubMed CentralGoogle Scholar
- Development Therapeutics Program NCI/NIH. http://www.dtp.nci.nih.gov/.
- Huang F, Reeves K, Han X, Fairchild C, Platero S, Wong TW, Lee F, Shaw P, Clark E: Identification of candidate molecular markers predicting sensitivity in solid tumors to dasatinib: rationale for patient selection. Cancer Res. 2007, 67: 2226-2238. 10.1158/0008-5472.CAN-06-3633.View ArticlePubMedGoogle Scholar
- Wang XD, Reeves K, Luo FR, Xu LA, Lee F, Clark E, Huang F: Identification of candidate predictive and surrogate molecular markers for dasatinib in prostate cancer: rationale for patient selection and efficacy monitoring. Genome Biol. 2007, 8: R255-10.1186/gb-2007-8-11-r255.View ArticlePubMedPubMed CentralGoogle Scholar
- Gemma A, Li C, Sugiyama Y, Matsuda K, Seike Y, Kosaihira S, Minegishi Y, Noro R, Nara M, Seike M, Yoshimura A, Shionoya A, Kawakami A, Ogawa N, Uesaka H, Kudoh S: Anticancer drug clustering in lung cancer based on gene expression profiles and sensitivity database. BMC Cancer. 2006, 6: 174-10.1186/1471-2407-6-174.View ArticlePubMedPubMed CentralGoogle Scholar
- Balko JM, Potti A, Saunders C, Stromberg A, Haura EB, Black EP: Gene expression patterns that predict sensitivity to epidermal growth factor receptor tyrosine kinase inhibitors in lung cancer cell lines and human lung tumors. BMC Genomics. 2006, 7: 289-10.1186/1471-2164-7-289.View ArticlePubMedPubMed CentralGoogle Scholar
- R language. http://www.r-project.org/.
- Zhu CQ, Ding K, Strumpf D, Weir BA, Meyerson M, Pennell N, Thomas RK, Naoki K, Ladd-Acosta C, Liu N, Pintilie M, Der S, Seymour L, Jurisica I, Shepherd FA, Tsao MS: Prognostic and predictive gene signature for adjuvant chemotherapy in resected non-small-cell lung cancer. J Clin Oncol. 2010, 28: 4417-4424. 10.1200/JCO.2009.26.4325.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Y, Klijn JG, Zhang Y, Sieuwerts AM, Look MP, Yang F, Talantov D, Timmermans M, Meijer-van Gelder ME, Yu J, Jatkoe T, Berns EM, Atkins D, Foekens JA: Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet. 2005, 365: 671-679.View ArticlePubMedGoogle Scholar
- Desmedt C, Piette F, Loi S, Wang Y, Lallemand F, Haibe-Kains B, Viale G, Delorenzi M, Zhang Y, d'Assignies MS, Bergh J, Lidereau R, Ellis P, Harris AL, Klijn JG, Foekens JA, Cardoso F, Piccart MJ, Buyse M, Sotiriou C: Strong time dependence of the 76-gene prognostic signature for node-negative breast cancer patients in the TRANSBIG multicenter independent validation series. Clin Cancer Res. 2007, 13: 3207-3214. 10.1158/1078-0432.CCR-06-2765.View ArticlePubMedGoogle Scholar
- Schmidt M, Bohm D, von Torne C, Steiner E, Puhl A, Pilch H, Lehr HA, Hengstler JG, Kolbl H, Gehrmann M: The humoral immune system has a key prognostic impact in node-negative breast cancer. Cancer Res. 2008, 68: 5405-5413. 10.1158/0008-5472.CAN-07-5206.View ArticlePubMedGoogle Scholar
- Ross DT, Scherf U, Eisen MB, Perou CM, Rees C, Spellman P, Iyer V, Jeffrey SS, Van de Rijn M, Waltham M, Pergamenschikov A, Lee JC, Lashkari D, Shalon D, Myers TG, Weinstein JN, Botstein D, Brown PO: Systematic variation in gene expression patterns in human cancer cell lines. Nat Genet. 2000, 24: 227-235. 10.1038/73432.View ArticlePubMedGoogle Scholar
- Roschke AV, Tonon G, Gehlhaus KS, McTyre N, Bussey KJ, Lababidi S, Scudiero DA, Weinstein JN, Kirsch IR: Karyotypic complexity of the NCI-60 drug-screening panel. Cancer Res. 2003, 63: 8634-8647.PubMedGoogle Scholar
- Bussey KJ, Chin K, Lababidi S, Reimers M, Reinhold WC, Kuo WL, Gwadry F, Kouros-Mehr H, Fridlyand J, Jain A, Collins C, Nishizuka S, Tonon G, Roschke A, Gehlhaus K, Kirsch I, Scudiero DA, Gray JW, Weinstein JN, Jay A: Integrating data on DNA copy number with gene expression levels and drug sensitivities in the NCI-60 cell line panel. Mol Cancer Ther. 2006, 5: 853-867. 10.1158/1535-7163.MCT-05-0155.View ArticlePubMedPubMed CentralGoogle Scholar
- Sieg DJ, Hauck CR, Ilic D, Klingbeil CK, Schaefer E, Damsky CH, Schlaepfer DD: FAK integrates growth-factor and integrin signals to promote cell migration. Nat Cell Biol. 2000, 2: 249-256. 10.1038/35010517.View ArticlePubMedGoogle Scholar
- Goi T, Shipitsin M, Lu Z, Foster DA, Klinz SG, Feig LA: An EGF receptor/Ral-GTPase signaling cascade regulates c-Src activity and substrate specificity. EMBO J. 2000, 19: 623-630. 10.1093/emboj/19.4.623.View ArticlePubMedPubMed CentralGoogle Scholar
- Danen EH, van Rheenen J, Franken W, Huveneers S, Sonneveld P, Jalink K, Sonnenberg A: Integrins control motile strategy through a Rho-cofilin pathway. J Cell Biol. 2005, 169: 515-526. 10.1083/jcb.200412081.View ArticlePubMedPubMed CentralGoogle Scholar
- Mitra SK, Hanson DA, Schlaepfer DD: Focal adhesion kinase: in command and control of cell motility. Nat Rev Mol Cell Biol. 2005, 6: 56-68. 10.1038/nrm1549.View ArticlePubMedGoogle Scholar
- Edwards DC, Sanders LC, Bokoch GM, Gill GN: Activation of LIM-kinase by Pak1 couples Rac/Cdc42 GTPase signalling to actin cytoskeletal dynamics. Nat Cell Biol. 1999, 1: 253-259. 10.1038/12963.View ArticlePubMedGoogle Scholar
- Sanders LC, Matsumura F, Bokoch GM, de Lanerolle P: Inhibition of myosin light chain kinase by p21-activated kinase. Science. 1999, 283: 2083-2085. 10.1126/science.283.5410.2083.View ArticlePubMedGoogle Scholar
- Araujo J, Logothetis C: Dasatinib: a potent SRC inhibitor in clinical development for the treatment of solid tumors. Cancer Treat Rev. 2010, 36: 492-500. 10.1016/j.ctrv.2010.02.015.View ArticlePubMedPubMed CentralGoogle Scholar
- Peng H, Sohara Y, Moats RA, Nelson MD, Groshen SG, Ye W, Reynolds CP, DeClerck YA: The activity of zoledronic Acid on neuroblastoma bone metastasis involves inhibition of osteoclasts and tumor cell survival and proliferation. Cancer Res. 2007, 67: 9346-9355. 10.1158/0008-5472.CAN-06-4508.View ArticlePubMedGoogle Scholar
- Santini D, Virzi V, Fratto ME, Bertoldo F, Sabbatini R, Berardi R, Calipari N, Ottaviani D, Ibrahim T: Can we consider zoledronic acid a new antitumor agent? Recent evidence in clinical setting. Curr Cancer Drug Targets. 2010, 10: 46-54. 10.2174/156800910790980223.View ArticlePubMedGoogle Scholar
- Nautiyal J, Majumder P, Patel BB, Lee FY, Majumdar AP: Src inhibitor dasatinib inhibits growth of breast cancer cells by modulating EGFR signaling. Cancer Lett. 2009, 283: 143-151. 10.1016/j.canlet.2009.03.035.View ArticlePubMedGoogle Scholar
- Bulgaru AM, Mani S, Goel S, Perez-Soler R: Erlotinib (Tarceva): a promising drug targeting epidermal growth factor receptor tyrosine kinase. Expert Rev Anticancer Ther. 2003, 3: 269-279. 10.1586/1473718.104.22.1689.View ArticlePubMedGoogle Scholar
- Houghton PJ: Everolimus. Clin Cancer Res. 2010, 16: 1368-1372. 10.1158/1078-0432.CCR-09-1314.View ArticlePubMedPubMed CentralGoogle Scholar
- Edick MJ, Tesfay L, Lamb LE, Knudsen BS, Miranti CK: Inhibition of integrin-mediated crosstalk with epidermal growth factor receptor/Erk or Src signaling pathways in autophagic prostate epithelial cells induces caspase-independent death. Mol Biol Cell. 2007, 18: 2481-2490. 10.1091/mbc.E06-04-0261.View ArticlePubMedPubMed CentralGoogle Scholar
- Kline CL, Olson TL, Irby RB: Src activity alters alpha3 integrin expression in colon tumor cells. Clin Exp Metastasis. 2009, 26: 77-87. 10.1007/s10585-008-9215-x.View ArticlePubMedGoogle Scholar
- Mitra SK, Schlaepfer DD: Integrin-regulated FAK-Src signaling in normal and cancer cells. Curr Opin Cell Biol. 2006, 18: 516-523. 10.1016/j.ceb.2006.08.011.View ArticlePubMedGoogle Scholar
- Webb DJ, Donais K, Whitmore LA, Thomas SM, Turner CE, Parsons JT, Horwitz AF: FAK-Src signalling through paxillin, ERK and MLCK regulates adhesion disassembly. Nat Cell Biol. 2004, 6: 154-161. 10.1038/ncb1094.View ArticlePubMedGoogle Scholar
- Gao P, Tchernyshyov I, Chang TC, Lee YS, Kita K, Ochi T, Zeller KI, De Marzo AM, Van Eyk JE, Mendell JT, Dang CV: c-Myc suppression of miR-23a/b enhances mitochondrial glutaminase expression and glutamine metabolism. Nature. 2009, 458: 762-765. 10.1038/nature07823.View ArticlePubMedPubMed CentralGoogle Scholar
- Seltzer MJ, Bennett BD, Joshi AD, Gao P, Thomas AG, Ferraris DV, Tsukamoto T, Rojas CJ, Slusher BS, Rabinowitz JD, Dang CV, Riggins GJ: Inhibition of glutaminase preferentially slows growth of glioma cells with mutant IDH1. Cancer Res. 2010, 70: 8981-8987. 10.1158/0008-5472.CAN-10-1666.View ArticlePubMedPubMed CentralGoogle Scholar
- Kutty RK, Kutty G, Samuel W, Duncan T, Bridges CC, El-Sherbeeny A, Nagineni CN, Smith SB, Wiggert B: Molecular characterization and developmental expression of NORPEG, a novel gene induced by retinoic acid. J Biol Chem. 2001, 276: 2831-2840. 10.1074/jbc.M007421200.View ArticlePubMedGoogle Scholar
- Sorrentino C, Di Carlo E: Expression of IL-32 in human lung cancer is related to the histotype and metastatic phenotype. Am J Respir Crit Care Med. 2009, 180: 769-779. 10.1164/rccm.200903-0400OC.View ArticlePubMedGoogle Scholar
- Shankar J, Messenberg A, Chan J, Underhill TM, Foster LJ, Nabi IR: Pseudopodial actin dynamics control epithelial-mesenchymal transition in metastatic cancer cells. Cancer Res. 2010, 70: 3780-3790. 10.1158/0008-5472.CAN-09-4439.View ArticlePubMedGoogle Scholar
- Aksoy S, Brandriff BF, Ward A, Little PF, Weinshilboum RM: Human nicotinamide N-methyltransferase gene: molecular cloning, structural characterization and chromosomal localization. Genomics. 1995, 29: 555-561. 10.1006/geno.1995.9966.View ArticlePubMedGoogle Scholar
- Xu J, Moatamed F, Caldwell JS, Walker JR, Kraiem Z, Taki K, Brent GA, Hershman JM: Enhanced expression of nicotinamide N-methyltransferase in human papillary thyroid carcinoma cells. J Clin Endocrinol Metab. 2003, 88: 4990-4996. 10.1210/jc.2002-021843.View ArticlePubMedGoogle Scholar
- Roessler M, Rollinger W, Palme S, Hagmann ML, Berndt P, Engel AM, Schneidinger B, Pfeffer M, Andres H, Karl J, Bodenmuller H, Ruschoff J, Henkel T, Rohr G, Rossol S, Rosch W, Langen H, Zolg W, Tacke M: Identification of nicotinamide N-methyltransferase as a novel serum tumor marker for colorectal cancer. Clin Cancer Res. 2005, 11: 6550-6557. 10.1158/1078-0432.CCR-05-0983.View ArticlePubMedGoogle Scholar
- Tomida M, Mikami I, Takeuchi S, Nishimura H, Akiyama H: Serum levels of nicotinamide N-methyltransferase in patients with lung cancer. J Cancer Res Clin Oncol. 2009, 135: 1223-1229. 10.1007/s00432-009-0563-y.View ArticlePubMedGoogle Scholar
- Wu Y, Siadaty MS, Berens ME, Hampton GM, Theodorescu D: Overlapping gene expression profiles of cell migration and tumor invasion in human bladder cancer identify metallothionein 1E and nicotinamide N-methyltransferase as novel regulators of cell migration. Oncogene. 2008, 27: 6679-6689. 10.1038/onc.2008.264.View ArticlePubMedGoogle Scholar
- Tang SW, Yang TC, Lin WC, Chang WH, Wang CC, Lai MK, Lin JY: Nicotinamide N-methyltransferase induces cellular invasion through activating matrix metalloproteinase-2 expression in clear cell renal cell carcinoma cells. Carcinogenesis. 2011, 32: 138-145. 10.1093/carcin/bgq225.View ArticlePubMedGoogle Scholar
- Staunton JE, Slonim DK, Coller HA, Tamayo P, Angelo MJ, Park J, Scherf U, Lee JK, Reinhold WO, Weinstein JN, Mesirov JP, Lander ES, Golub TR: Chemosensitivity prediction by transcriptional profiling. Proc Natl Acad Sci USA. 2001, 98: 10787-10792. 10.1073/pnas.191368598.View ArticlePubMedPubMed CentralGoogle Scholar
- Bild AH, Yao G, Chang JT, Wang Q, Potti A, Chasse D, Joshi MB, Harpole D, Lancaster JM, Berchuck A, Olson JA, Marks JR, Dressman HK, West M, Nevins JR: Oncogenic pathway signatures in human cancers as a guide to targeted therapies. Nature. 2006, 439: 353-357. 10.1038/nature04296.View ArticlePubMedGoogle Scholar
- Neve RM, Chin K, Fridlyand J, Yeh J, Baehner FL, Fevr T, Clark L, Bayani N, Coppe JP, Tong F, Speed T, Spellman PT, DeVries S, Lapuk A, Wang NJ, Kuo WL, Stilwell JL, Pinkel D, Albertson DG, Waldman FM, McCormick F, Dickson RB, Johnson MD, Lippman M, Ethier S, Gazdar A, Gray JW: A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes. Cancer Cell. 2006, 10: 515-527. 10.1016/j.ccr.2006.10.008.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee JK, Havaleshko DM, Cho H, Weinstein JN, Kaldjian EP, Karpovich J, Grimshaw A, Theodorescu D: A strategy for predicting the chemosensitivity of human cancers and its application to drug discovery. Proc Natl Acad Sci USA. 2007, 104: 13086-13091. 10.1073/pnas.0610292104.View ArticlePubMedPubMed CentralGoogle Scholar
- Hsu YC, Yuan S, Chen HY, Yu SL, Liu CH, Hsu PY, Wu G, Lin CH, Chang GC, Li KC, Yang PC: A four-gene signature from NCI-60 cell line for survival prediction in non-small cell lung cancer. Clin Cancer Res. 2009, 15: 7309-7315. 10.1158/1078-0432.CCR-09-1572.View ArticlePubMedGoogle Scholar
- Dai Q, Ling YH, Lia M, Zou YY, Kroog G, Iwata KK, Perez-Soler R: Enhanced sensitivity to the HER1/epidermal growth factor receptor tyrosine kinase inhibitor erlotinib hydrochloride in chemotherapy-resistant tumor cell lines. Clin Cancer Res. 2005, 11: 1572-1578. 10.1158/1078-0432.CCR-04-0993.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1741-7015/11/106/prepub