
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Department of Gastroenterology, Graduate School of Medicine, University of Tokyo, Tokyo 113-8655 [M. M., Y. H., M. O., N. K., T. G., H. T., Y. S., M. O.]; Cellular Informatics Team, Computational Biology Research Center, Tokyo 135-0064 [S. N.]; and Department of Functional Genomics, Graduate School of Medicine, Chiba University, Chiba 260-8670, [S. N.], Japan
2 To whom requests for reprints should be addressed, at Department of Gastroenterology, Graduate School of Medicine, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8655, Japan. Phone: 81-3-3815-5411, extension 33056; Fax: 81-3-3814-0021; E-mail: moriyamam-int{at}h.u-tokyo.ac.jp
| Abstract |
|---|
|
|
|---|
In eight hepatoma cell lines (HLE, HLF, Huh7, Hep3B, PLC/PRF/5, SK-Hep1, Huh6, and HepG2), the baseline expression levels of 2300 genes were measured by cDNA microarray. The concentrations of eight anticancer drugs (nimustine, mitomycin C, cisplatin, carboplatin, doxorubicin, epirubicin, mitoxantrone, and 5-fluorouracil) needed for 50% growth inhibition were examined and used as a measure of chemosensitivity. These data were combined and comprehensive pair-wise correlations between gene expression levels and the 50% growth inhibition values were calculated. Significant correlations with significance were used to construct networks of similarity.
Fifty-two relations, including 42 genes, were selected. Among them, nearly 20% were various types of transporters, and most of them negatively correlated with chemosensitivity. Transporter associated with antigen processing 1 was associated with resistance to mitoxantrone, consistent with previous reports. Other transporters were not reported previously to associate with chemosensitivity. Resistance to doxorubicin and its analogue, epirubicin, were positively correlated with topoisomerase II ß expression, whereas it negatively correlated with expression of carboxypeptidases A3 and Z. Response to nimustine was associated with expression of superoxide dismutase 2.
Relevance networks identified several negative correlations between gene expression and resistance, which were missed by hierarchical clustering. Our results suggested the necessity of systematically evaluating the transporting systems that may play a major role in resistance in hepatoma. This may provide useful information to modify anticancer drug action in hepatoma.
| Introduction |
|---|
|
|
|---|
The development of hepatoma is a complicated process, associated with numerous genetic factors. Large-scale transcriptional profiling allows us to obtain large amounts of data about complicated biological processes. However, to interpret the result of DNA microarray experiments, we have to extract the biologically significant information from the many changes that are random or not related to the phenomena studied.
Relevance networks offer a method to construct networks of similarity with various type of information, including mutual information and the correlation status (6). In this study, we applied this methodology to provide new information on chemosensitivity in hepatoma.
| Materials and Methods |
|---|
|
|
|---|
Anticancer Drugs.
We tested eight drugs used currently for the treatment of hepatoma in clinical practice: an antimetabolite (5-FU 3), two platinum anticancer agents (CDDP and CBDCA), two alkylating agents (ACNU and MMC), and three (topo) II inhibitors (ADM, EpiADM, and MIT). CDDP, CBDCA, ADM, MMC, MIT, and 5-FU were purchased from Sigma Co. (St. Louis, MO). ACNU was obtained from Sankyo Co. (Tokyo, Japan). EpiADM was obtained from Pharmacia Co. (Peapack, NJ).
Measurement of Response to Drugs.
The response of each cell line to the anticancer drugs was determined as follows. Briefly, cells were grown in 24-well microplates (IWAKI Glass, Chiba, Japan) and exposed to the drugs for 48 h. Cell growth was assessed by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide assay. The drug concentration required to inhibit cell growth by 50% in comparison with untreated controls was designated as the growth inhibition activity of the drug (GI50). Data were averages of three independent experiments.
RNA Extraction.
Total cellular RNA was extracted using an acid guanidine isothiocyanate phenol-chloroform method according to the manufacturers instructions (ISOGEN Reagent; Nippon Gene Co., Tokyo, Japan). RNA was treated with RNase-free DNase to remove contaminating genomic DNA, and RNA integrity was confirmed by gel electrophoresis and ethidium bromide staining. Polyadenylic acid mRNA was obtained using the Oligotex-dT30 mRNA purification kit (TaKaRa Co., Tokyo, Japan).
cDNA Microarray.
We assessed gene expression using in-house cDNA microarray developed by ourselves in the Helix Research Institute (Chiba, Japan; Ref. 7). A total of 2300 named, sequence-validated human cDNAs (Research Genetics) were spotted onto carbodiimide-coated glass slides using a robot (SPBIO-2000; Hitachi Software Engineering Co., Tokyo, Japan).
Microarray experiments were performed as described previously (8), and for each gene, expression was compared with expression in normal human liver RNA (Clontech, Palo Alto, CA) using the Cy5:Cy3 ratios. To control for labeling differences, reactions were carried out at least in duplicate, and the fluorescent dyes were switched. Each pair of probes was hybridized to a separate microarray. We have confirmed previously the reproducibility of such experiments (710).
Microarray Data Analysis.
Spots were only included if their raw fluorescence intensities were at least 1.5 times the local background. We adjusted raw data of fluorescence ratios (Cy5:Cy3) by log transformation (base 2), median centering, and normalization. Subsequently, we filtered data on various genes for additional analyses. Genes were included if one or more cell lines had fluorescence ratios
2.0 or
0.5, and maximum minus minimum values of log fluorescent ratios
1.0.
To ascertain that the range of values of each gene expression ratio was sufficient for subsequent analysis, an entropy value, H, was calculated using the following formula:
![]() |
Relevance Networks.
Using log-transformed expression ratios of selected genes and the GI50 values of the anticancer drugs, similarity metrics, R, were calculated in a comprehensive pair-wise manner as a square of the Pearson correlation coefficient with original sign (positive or negative). The correlation coefficient, r, was calculated by the following formula:
![]() |
| Results |
|---|
|
|
|---|
|
6% of genes on our microarray encode transporters. Approximately 10% of selected genes were associated with cell cycle control and cell proliferation.
|
|
| Discussion |
|---|
|
|
|---|
Previous microarray research on this disease showed that some growth factors and detoxification-related genes might be associated with its cancer development (1218). In comparison, our current study suggested that low levels of expression of specific intra- and intercellular transporters might be associated with its resistance to drugs. Most of the transporters showed negative correlation with the GI50 values, and these correlations were not identified by hierarchical clustering. The gene expression profiles of these transport systems may be an efficient and practical way to understand chemosensitivity, because most of the genes code for membrane proteins, and it is difficult to analyze their function individually (in fact, none of them is included in the Protein Data Bank entry).
Our analysis also picked up several genes reported previously to be associated with response to anticancer drugs. TOP2B was associated with resistance to two anthracyclines, ADM and EpiADM. The
isoform of TOP2 is a target for anthracycline and anthracenedione, and its transcription was positively correlated with resistance to topo II inhibitors. The association between TOP2B expression and the resistance of cancers to topo II inhibitors correlated positively with susceptibility to topo II inhibitors, whereas the association between TOP2B expression and the resistance to cancers to topo II inhibitors varies among different organs and different drugs (11). In our experiments, TOP2B expression was positively correlated with resistance to ADM and EpiADM.
Expression of two carboxypeptidases, CPA3 and CPZ, was negatively correlated with resistance to ADM and EpiADM. A previous study showed that carboxypeptidase G2, an enzyme that hydrolyzes the COOH-terminal glutamate residue, improves response to another anticancer drug, methotrexate (19). It may be possible that in hepatoma, CPA3 and CPZcould decrease resistance to ADM and EpiADM.
CA1, which was reportedly associated with colon cancer development (20), showed negative correlation with resistance to EpiADM. TRG also showed negative correlation with resistance to EpiADM. This is consistent with prior reports that polymorphisms in TRG are associated with early onset colorectal cancer (21).
Superoxide, a toxic free radical, reportedly plays an important role in anticancer effects of anthracycline (11). SOD2 inactivates this free radical, but in this study its expression was not associated with resistance to the anthracycline derivates tested. On the other hand, SOD2 expression is negatively correlated with resistance to nitrosourea in various cancer cells (22). Our data also indicated a negative correlation between SOD2 expression and resistance to ACNU, also a nitrosourea.
TAP1, a member of multidrug resistance family, is associated with resistance to MIT (23). Consistent with this, we found that TAP1 expression showed a specific positive correlation with MIT GI50.
Cluster analysis does not always provide a complete interpretation of microarray data, and the application of more than one analysis technique can illuminate different relationships (24). Relevance network analysis can capture negative correlations that are missed in hierarchical clustering and simultaneously display them with positive correlations. On the other hand, relevance network analysis has several limitations. First, large-scale networks may not be detected if individual links are weak, because during analysis all of the links that are subthreshold are discarded. Second, the analysis may not necessarily assign equal statistical weight to all of the observations. Variances of measurements are inhomogeneous; that is, there may be very little variation in most of the samples and great variation in only a few of the samples, and the analysis may give greater weight to the latter. Third, this is not a globally optimal network. Fourth, comprehensive pair-wise calculation of the dissimilarity measure requires extensive computing resources (25).
In summary, our results highlighted the importance of various types of transporters that may affect chemosensitivity in hepatoma. However, most of them are ubiquitous molecules, and, thus, a study of these alone will not allow us to understand whole mechanism of chemosensitivity. Establishment of systems biology models, e.g., focused on transport systems rather than one or two individual transporters, will help to understand anticancer drug action in hepatoma (26).
| Footnotes |
|---|
1 Supported in part by the Health Science Research Grants for Medical Frontier Strategy Research from the Ministry of Health, Labor, and Welfare of Japan, and by grants-in-aid for scientific research from the Ministry of Education, Culture, Sports, Science and Technology, Japan. ![]()
3 The abbreviations used are: 5-FU, 5-fluorouracil; ACNU, nimustine; ADM, doxorubicin; CBDCA, carboplatin; CDDP, cisplatin; CA1, carbonic anhydrase 1; CPA3, carboxypeptidase A3; CPZ, carboxypeptidase Z; EpiADM, epirubicin; GI50, 50% growth inhibition activity of the drug; topo, topoisomerase; MIT, mitoxantrone; MMC, mitomycin C; SOD, superoxide dismutases; TAP, transporter associated with antigen processing; TRG, T-cell receptor
locus; TOP2B, DNA topoisomerase II ß. ![]()
4 Internet address: http://www.bioperl.org/. ![]()
5 Internet address: http://www.ncbi.nih.gov/. ![]()
6 Internet address: http://www.geneontology.org/. ![]()
7 Internet address: http://www.genome.ad.jp/kegg/. ![]()
Received 12/ 2/02; accepted 12/ 2/02.
| References |
|---|
|
|
|---|
: a microsatellite marker for colorectal cancer.Ann. Surg. Oncol.
,
9:88
93,2002
.This article has been cited by other articles:
![]() |
M. Tamura and P. D'haeseleer Microbial genotype-phenotype mapping by class association rule mining Bioinformatics, July 1, 2008; 24(13): 1523 - 1529. [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] |
||||
![]() |
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] |
||||
![]() |
L. I. Gomes, G. H. Esteves, A. F. Carvalho, E. B. Cristo, R. Hirata Jr., W. K. Martins, S. M. Marques, L. P. Camargo, H. Brentani, A. Pelosof, et al. Expression Profile of Malignant and Nonmalignant Lesions of Esophagus and Stomach: Differential Activity of Functional Modules Related to Inflammation and Lipid Metabolism Cancer Res., August 15, 2005; 65(16): 7127 - 7136. [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 |