
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Oncology, Markers, Clinical Correlates |
Combination Chemotherapy in Advanced Hepatocellular Carcinoma
1 Department of Surgery and Clinical Oncology, Graduate School of Medicine, Osaka University, Osaka, Japan; and 2 Taisho Laboratory of Functional Genomics and 3 Laboratory of Theoretical Life Science, Nara Institute of Science and Technology, Nara, Japan
ABSTRACT
Purpose: The prognosis of hepatocellular carcinoma (HCC) is very poor, particularly in patients with tumors that have invaded the major branches of the portal vein. Combination chemotherapy with intra-arterial 5-fluorouracil and subcutaneous interferon-
has shown promising results for such advanced HCC, but it is important to develop the ability to accurately predict chemotherapeutic responses.
Experimental Design: We analyzed the expression of 3,080 genes using a polymerase chain reaction-based array in 20 HCC patients who were treated with combination chemotherapy after reduction surgery. After unsupervised analyses, a supervised classification method for predicting chemotherapeutic responses was constructed. To minimize the number of predictive genes, we used a random permutation test to select only significant (P < 0.01) genes. A leave-one-out cross-validation confirmed the gene selection. We also prepared an additional 11 cases for validation of predictive performance.
Results: Hierarchical clustering analysis and principal component analysis with all 3,080 genes revealed distinct gene expression patterns in responders (those with complete response or partial response) and nonresponders (those with stable disease or progressive disease) to the combination chemotherapy. Using a weighted-voting classification method with either all genes or only significant genes as assessed by permutation testing, the objective responses to treatment were correctly predicted in 17 of 20 cases (accuracy, 85%; positive predictive value, 100%; negative predictive value, 80%). Moreover, patients in the validation dataset could be classified into two distinct prognostic groups using 63 predictive genes.
Conclusions: Molecular analysis of 63 genes can predict the response of patients with advanced HCC and major portal vein tumor thrombi to combination chemotherapy with 5-fluorouracil and interferon-
.
INTRODUCTION
Hepatocellular carcinoma (HCC), the predominant histologic subtype of primary liver cancer, is one of the major causes of death from malignancy worldwide. Although recent progress in both diagnostic and surgical techniques has resulted in considerable improvement in the morbidity and mortality rates, the overall outcome remains far from satisfactory. The prognosis is miserable, particularly in patients with tumors that have invaded the major branch of the portal vein, and the median survival time is several months (1, 2, 3)
. Indeed, the 1-year survival rate was <50% even in resectable cases (3
, 4)
. We and others have recently developed a new treatment regimen for such advanced HCC patients: combination chemotherapy of 5-fluorouracil (5-FU) and interferon (IFN)-
(5
, 6)
. In follow-up studies, about 40% to 50% of patients showed promising responses to the therapy and lived >2 years. However, for patients who did not respond to this treatment, survival was too short to receive another type of treatment. Furthermore, this chemotherapy may have multiple adverse effects, including leukopenia, thrombocytopenia, and depression. Therefore, accurate prediction of chemosensitivity is desirable, not only so that nonresponding patients do not lose a limited chance to take advantage of other possible treatments, but also to eliminate suffering of these patients due to debilitating side effects. Unfortunately, there are currently no useful indicators to distinguish between patients who are likely to respond to this combination chemotherapy and patients who are not.
In this study, genes possessing the ability to predict patient responses to 5-FU and IFN-
combination chemotherapy were selected by gene expression profile analysis using adaptor-tagged competitive polymerase chain reaction (ATAC-PCR) technology, a polymerase chain reaction (PCR)-based array system (7)
. The performance of this prediction method was estimated with both a complete cross-validation and an external validation dataset.
MATERIALS AND METHODS
Tissues and Patients.
Between August 1998 and January 2003, 20 HCC patients with multiple tumors spreading to bilateral lobes with tumor thrombi in the major branches of the portal vein underwent palliative surgery to reopen the portal flow and recover liver function. After obtaining informed consent, we collected HCC tissue from the main resected tumors and isolated the total RNA from these samples for PCR-based array experiments. All 20 patients had visible tumors in the remaining liver. The patients were then subjected to a treatment regimen of combination chemotherapy with 5-FU and IFN-
as described previously (5)
. Briefly, patients received continuous arterial infusion of 5-FU (450500 mg/d) for the initial 2 weeks and subcutaneous IFN-
injection (5 million IU) 3 times per week for 4 weeks. The sizes of the remaining tumors in the liver, which are thought to be intrahepatic metastases from the resected primary tumor, were measured before and after treatment by either abdominal computed tomography scan or ultrasonography. The chemotherapeutic response was clinically evaluated according to the Eastern Cooperative Oncology Group (ECOG) criteria (8)
. In this study, responders were defined as patients with complete response (CR) or partial response (PR); nonresponders were defined as patients with stable disease or progressive disease (PD).
Furthermore, as an independent validation dataset, we collected liver tissue specimens and clinical data from 11 HCC patients who underwent hepatic resection of all visible tumors between October 1998 and October 2002. Although these patients had tumor(s) with major portal vein tumor thrombi, the area of spreading tumor(s) was limited to two segments of the liver. The treatment regimen of combination chemotherapy and the method of follow-up in the validation dataset were the same as the procedures used in the original dataset. All aspects of our study protocol were approved by the ethics committee of Osaka University Medical School.
Polymerase Chain Reaction-Based Array System.
To select genes expressed in liver tissues, we constructed three cDNA libraries: one from a mixture of HCC and nontumorous livers, one from normal livers, and one from metastatic liver cancers, as described previously (9)
. We designed PCR primers for ATAC-PCR reactions for a total of 2,666 genes from these expressed sequence tag collections. In total, we prepared 3,080 primers for ATAC-PCR; this total includes an additional 414 genes established in previous literature. The specificity of this gene selection provides an advantage over more universal gene sets, such as those selected from the UniGene database, which include genes not expressed in liver tissues. The ATAC-PCR experimental procedure was performed as described previously (10)
. The complete list of genes and detailed protocols for the ATAC-PCR experimental procedure are available on our web site.4
Analysis of Polymerase Chain Reaction-Based Array Data.
The relative expression level of each gene was calculated by calibration using a standard mixture of normal liver tissues, as described previously (11, 12, 13, 14, 15)
. After conversion to a logarithmic scale (base 2), the data matrix was normalized to a median of 0 by standardizing each sample.
As a supervised classification method, we adopted a weighted-voting (WV) algorithm generally used in gene expression profiling (16, 17, 18, 19)
. Briefly, we calculated the signal to noise (S2N) ratio, Si = (µR µN)/(
R +
N), where µ and
represent the mean and SD of expression for responders and nonresponders, respectively. The magnitude of the gene vote (vi) reflects the deviation of the test sample Xi value from the average of the two groups: vi = Si x (Xi (µR + µN)/2). We summed the vi values to obtain the total votes for a "good signature" group (VG) and a "poor signature" group (VP). The prediction strength is (VG |VP|)/(VG + |VP|), and we adopted a threshold of 0. If the strength is a positive number, the test sample belongs to the good signature group. If the strength is a negative number, the test sample belongs to the poor signature group. This model was evaluated by leave-one-out cross-validation (16, 17, 18, 19, 20, 21)
, wherein one sample was randomly withheld, and the model was regenerated with the remaining samples and then used to predict the group of the withheld sample. This process was repeated until every sample was tested, and the cumulative accuracy was recorded. When an external validation dataset was used to evaluate this WV model, we used the S2N ratio (Si) obtained from the original dataset. In calculating the weighted vote (vi), the mean expression level of the entire validation dataset was substituted for the average of the two groups ((µR + µN)/2) in the original dataset.
Permutation testing, which involves randomly permuting group labels to determine gene-group correlations, was used to determine statistical significance (16, 17, 18, 19) . The original score, So = |µR µN|, of each gene was calculated without permuting the labels (responder or nonresponder). The labels of all of the samples were then randomly permuted, and the scores were recalculated between two groups consisting of the new members. Repetition of this permutation 50,000 times provides a score number (No) larger than the original score (So). For each gene, the permutation P value was determined by P = No/50,000. To estimate the false discovery rate, which is the percentage of genes identified by chance (22) , we calculated the Q value of every gene using downloaded software.5 The Q value provides a measure of each features significance, automatically taking into account the problem of multiple testing (23) .
Hierarchical cluster analysis using an unweighted pair group method using arithmetic mean method with Pearsons correlation and principal component analysis (PCA) was performed using GeneMaths 2.0 software. Other statistical analyses were performed using StatView 5.0J software. Correlations between the responses to our combination chemotherapy and clinicopathological parameters were evaluated by the
2 test with the Yates correction. Overall and disease-free survival rates were calculated using the Kaplan-Meier method. Differences in survival curves were estimated by the log-rank test.
RESULTS
Patient Characteristics.
The characteristics of the 20 HCC patients in the original dataset are shown in Table 1
. All patients had multiple tumors spreading to bilateral lobes with tumor thrombi in the major branches of the portal vein and were treated with 5-FU and IFN-
combination chemotherapy after reduction surgery. Fifteen of the 20 patients were treated with three cycles of the above-mentioned regimen. The remaining five patients received only one or two cycles of 5-FU due to extensive progression of the tumors. According to the ECOG criteria of objective response, 8 patients (40%) were classified as responders, demonstrating either CR or PR for at least 4 weeks. The remaining 12 patients (60%) were classified as nonresponders, exhibiting either stable disease or PD. The clinicopathological characteristics of responders (CR and PR) and nonresponders (stable disease and PD) were compared by the
2 test (Table 2)
. There were no significant differences (P < 0.05) in any factor between the two response groups.
|
|
combination chemotherapy, i.e., responders and nonresponders, with only a few exceptions (Fig. 1)
|
|
combination chemotherapy was evaluated by leave-one-out cross-validation. In this first cross-validation analysis, we used all 3,080 genes to avoid any selection bias (Fig. 3A)
|
|
|
|
HCC remains a major problem worldwide and appears to be increasing in developed Western countries (25)
. HCC is notorious for its poor prognosis because the overall 5-year survival rate after resection has remained as poor as 35% to 50%, even if curative surgery is performed (26, 27, 28)
. Moreover, once the tumor invades the major branch of the portal vein, the prognosis is extremely poor (1, 2, 3, 4)
. Various complications including gastrointestinal bleeding, ascites, and hemorrhage from esophageal varices are provoked by stenosis of the portal vein. Until now, there has been no standard treatment for advanced HCC with portal tumor thrombi. However, we and others have recently developed a new combination chemotherapy regimen using intra-arterial 5-FU and subcutaneous IFN-
(5
, 6)
, which has proven to be effective in about half of such advanced HCC patients. In this study, we have demonstrated that a molecular prediction method using expression data of 63 genes was able to distinguish patients who responded to this therapy from those who did not respond.
DNA microarray technology allows parallel expression analysis of thousands of genes to address complex questions in tumor biology. Many trials predicting the prognosis of various human malignancies have been reported using DNA microarrays (29, 30, 31) . Moreover, a novel prediction method using gene expression profiling has recently been reported for the treatment of breast cancer patients with the taxanes (docetaxel and paclitaxel; ref. 21 ). Although DNA microarrays have contributed to such expression profiling studies to some degree, DNA microarrays can detect only a fraction of the changes in gene expression detectable by reverse transcription-PCR (32) . We therefore performed a high-throughput quantitative PCR based on ATAC-PCR (7) to analyze the genetic differences in HCC. This assay requires smaller amounts of RNA than DNA microarray analysis. PCR-based analysis of selected genes costs much less than DNA microarrays, which are likely to require at least several hundred spotted genes for diagnosis. The aforementioned benefits and the strength of this system for cancer research, established in previous reports on HCC, breast, and colorectal cancers make this technique a powerful means to obtain a better understanding of the molecular characteristics of cancers (11, 12, 13, 14, 15) .
In the unsupervised analyses with all examined genes, the difference of expression patterns between responders (CR and PR) and nonresponders (stable disease and PD) was clear. Of course, it was controversial how to divide two response groups. However, the response to this regimen with 5-FU and IFN-
was usually typical; that is, PR or stable disease cases were rare, as shown in Table 1
. Clinically, radiologically, and biochemically, the PR case in this study revealed a near CR, and the stable disease case was similar to PD cases. We therefore classified the PR patient as a responder and the stable disease patient as a nonresponder due to the conventional grouping method (33
, 34)
. In supervised analyses, we selected 63 significant genes by a random permutation test to minimize the number of predictive genes. These 63 predictive genes included some genes that have previously been reported to be associated with sensitivity to 5-FU or IFN. For example, a gene related to 5-FU sensitivity is tumor necrosis factor ligand superfamily member 6 (TNFSF6), also called FAS (CD95) ligand, which plays a major role in the induction of apoptosis in response to a variety of extracellular signals, including anticancer drugs. In hepatoma cells or human breast cancer, up-regulation of TNFSF6 was associated with poor sensitivity to 5-FU chemotherapy (35
, 36)
. Matrix metalloproteinase 9 (MMP9), a member of the secreted zinc metalloprotease group, is also inhibited by 5-FU treatment (37)
. Hexokinase 2 (HK2) plays an important role in intracellular glucose metabolism and is strongly associated with the metabolism of 2'-fluoro-2'-deoxy-D-mannose, which has been correlated to 5-FU response (38)
. A gene involved in sensitivity to IFN-
is OASL, a member of the 2'-5'-oligoadenylate synthetase gene families. OASL plays an important role in the antiviral effects of IFN and is involved in apoptosis and control of cellular growth (39
, 40)
. Baculoviral IAP repeat-containing 4 (BIRC4; XIAP), encodes a protein that significantly inhibits apoptosis. IFN inhibited BIRC4 function by inducing XIAP associated factor-1 (XAF-1; ref. 41
). Glutamyl aminopeptidase is a differentiation-related kidney glycoprotein of 160 kDa (gp160), and the absence of glutamyl aminopeptidase expression was predictive of IFN-
sensitivity in renal cell carcinoma (42)
. Calmodulin-1 (CALM1) plays a role in growth and the cell cycle and in signal transduction. Calmodulin inhibitors affect not only the uptake and processing of IFN inducers but also the release of IFN from induced cells (43)
. Our selection approach to these predictive genes was confirmed with both a complete cross-validation and an external validation dataset.
In our univariate analysis with clinicopathological factors, there was no significant indicator that distinguished responders and nonresponders. Patt et al. (6)
reported that low levels (
50 ng/mL) of serum
-fetoprotein (AFP) were a predictor of response to 5-FU and IFN-
combination chemotherapy, but the serum AFP level was not a significant indicator of chemotherapeutic response in our study. Because we treated only extremely advanced HCC patients with major portal vein tumor thrombi with this combination therapy, only 3 of the 20 patients showed low levels of serum AFP. It may be too early to make conclusions regarding this factor, but we believe that another appropriate method will be required to more accurately predict the response to this treatment.
Indeed, this combination chemotherapy of 5-FU and IFN-
is promising for advanced HCC, but we need to consider the side effects. IFN-
and/or 5-FU sometimes induce severe adverse effects, including myelosuppression, fever, and depression (44)
. In particular, adverse myelosuppression is an important factor in HCC cases not only because thrombocytopenia and/or leukopenia are frequently present before chemotherapy, but also because treatment often has to be discontinued due to these side effects. Moreover, the survival of patients with such extremely advanced HCC was too short to allow a secondary treatment (1, 2, 3, 4)
. Therefore, accurate prediction of sensitivity to the first-line chemotherapy is necessary for advanced HCC patients so that they do not lose their limited chance to receive treatment and the patients who will not respond to the therapy can be protected from debilitating side effects.
These results illustrate the potential of biological technology to advance diagnostic methods, allowing physicians to plan beyond empirical results toward a more molecularly welldefined, personalized therapy. Particularly for patients with poor prognosis, such as advanced HCC patients, accurate prediction of chemotherapeutic responses before treatment is necessary. Because our molecular prediction system involves only a small number of predictive genes and a simplified algorithm that does not require either statistical software or specialists, this system should easily lead to future clinical application.
ACKNOWLEDGMENTS
We thank Chiyuri Maruyama, Keiko Ikenaga, Mihoko Yoshino, and Satoko Kinjo for expert technical assistance as well as Dr. Hideo Ota for helpful suggestions.
FOOTNOTES
Grant support: Grant-in-Aid for the Development of Innovative Technology from the Ministry of Education, Culture, Sports, Science and Technology, Japan; Grants-in-Aid for the Second Term Comprehensive 10-Year Strategy for Cancer Control and Cancer Research from the Ministry of Health and Welfare, Japan; and Grants-in-Aid for Scientific Research on Priority Areas, and Basic Research from the Ministry of Education, Science, Sports, and Culture, Japan.
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.
Requests for reprints: Morito Monden, Department of Surgery and Clinical Oncology, Graduate School of Medicine, Osaka University, 2-2, Yamadaoka, Suita city, Osaka 565-0871, Japan. Phone: 81-6-6879-3251; Fax: 81-6-6879-3259; E-mail: monden{at}surg2.med.osaka-u.ac.jp
4 http://love2.aist-nara.ac.jp/laboratory/index_frame.html. ![]()
5 http://faculty.washington.edu/
jstorey/qvalue. ![]()
Received 2/ 6/04; revised 6/ 9/04; accepted 6/17/04.
REFERENCES
This article has been cited by other articles:
![]() |
R. W. C. Pang, J. W. Joh, P. J. Johnson, M. Monden, T. M. Pawlik, and R. T. P. Poon Biology of Hepatocellular Carcinoma Ann. Surg. Oncol., April 1, 2008; 15(4): 962 - 971. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Shirahata, K. Iwao-Koizumi, S. Saito, N. Ueno, M. Oda, N. Hashimoto, J. A. Takahashi, and K. Kato Gene Expression-Based Molecular Diagnostic System for Malignant Gliomas Is Superior to Histological Diagnosis Clin. Cancer Res., December 15, 2007; 13(24): 7341 - 7356. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. K. Abou-Alfa, L. Schwartz, S. Ricci, D. Amadori, A. Santoro, A. Figer, J. De Greve, J.-Y. Douillard, C. Lathia, B. Schwartz, et al. Phase II Study of Sorafenib in Patients With Advanced Hepatocellular Carcinoma J. Clin. Oncol., September 10, 2006; 24(26): 4293 - 4300. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. N. Vemuri and A. A. Aristidou Metabolic Engineering in the -omics Era: Elucidating and Modulating Regulatory Networks Microbiol. Mol. Biol. Rev., June 1, 2005; 69(2): 197 - 216. [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 |