The Oncologist, Vol. 12, No. 4, 406-417, April 2007; doi:10.1634/theoncologist.12-4-406 © 2007 AlphaMed Press
BAC Clones Related to Prognosis in Patients with Esophageal Squamous Carcinoma: An Array Comparative Genomic Hybridization StudyaDepartment of Surgery and Molecular Oncology, Medical Institute of Bioregulation, Kyushu University, Beppu, Japan; bDepartment of Surgical Oncology, Graduate School of Medicine, Tokyo Medical and Dental University, Tokyo, Japan; cResearch Institute for Cell Engineering, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Japan; dDepartment of Oncological Science (Surgery II), Faculty of Medicine, Oita University, Yufu, Japan Key Words. Laser microdissection • Genome • Microarray • Clinical samples • Survival curve Correspondence: Masaki Mori, M.D., Ph.D., F.A.C.S., Department of Surgery and Molecular Oncology, Medical Institute of Bioregulation, Kyushu University, 4546 Tsurumibaru, Beppu 874-0838, Japan. Telephone: 81-977-27-1650; Fax: 81-977-27-1651; e-mail: mmori{at}beppu.kyushu-u.ac.jp or Takashi Hirano, Ph.D., Research Institute for Cell Engineering, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Higashi, Tsukuba, Ibaraki 305-8566, Japan. Telephone: 81-29-861-6152; Fax: 81-29-861-6144; e-mail: hirano-takashi{at}aist.go.jp Received December 4, 2006; accepted for publication February 1, 2007.
Purpose. The prognosis of patients with esophageal carcinoma is poor. To identify genomic alterations associated with poor patient prognosis, we analyzed whole DNA copy number profiles of esophageal squamous carcinomas (ESCs) using array-based comparative genomic hybridization (aCGH). Materials and Methods. Twenty-one operated and two biopsied cases of esophageal squamous cancer were examined for study. Each sample was laser microdissected to obtain pure cancer cell populations. The extracted DNA was analyzed using aCGH. Results. One of the most representative alterations was a previously reported amplification at 11q13.3. In addition, some novel alterations, such as deletion of 16p13.3, were identified. Of the 19 patients who were reassessed more than 5 years after the operation, nine were still living and 10 had died from disease recurrence. When aCGH profiles from the surviving group and the deceased group were compared, significant differences were recognized in 68 of 4,030 bacterial artificial chromosome (BAC) clones. Almost half of these clones were present at nine limiting regions in 4q, 13q, 20q, and Xq. For 22 of these 68 BAC clones, there also was a significant difference in the Kaplan-Meier survival curve, using the log-rank test, when comparing patients who had an alteration in a particular clone with those who did not. Conclusions. aCGH study of esophageal squamous cancer clearly identified BAC clones that are related to the prognosis of patients. These clones give us the opportunity to determine specific genes that are associated with cancer progression. Disclosure of potential conflicts of interest is found at the end of this article.
Comparative genomic hybridization (CGH) has been developed to analyze whole chromosomal aberrations of an entire genome [1]. This method has identified chromosomal aberrations in several kinds of cancer. Recently, DNA microarrays have been applied to CGH (array CGH), and this combined technique has a great potential for comprehensive analysis of both the relative DNA copy number and altered chromosomal regions in cancer [23]. In addition, the completion of the human genome sequence has enabled identification of candidate genes within regions of DNA associated with cancers [4]. Esophageal cancer is one of the most prevalent cancers worldwide. Most esophageal cancers seen in Japanese patients show squamous cancer histologically, while adenocarcinomas are frequently seen in Caucasians [5]. The 5-year survival rate of esophageal cancer is <40% despite the development of treatment modalities such as surgery, chemotherapy, and radiotherapy [6]. Thus, in order to develop better treatments, it is desirable to identify genomic alterations associated with esophageal squamous cancers (ESCs). Several previous studies have used chromosomal CGH (cCGH) to identify chromosomal aberrations associated with ESC. Noguchi et al. [7] studied the relationship between a gain in the 3p chromosomal region and tumor progression. Yen et al. [8] studied the relationship between a gain in 5p and 7q and a loss in 4p, 9p, and 11q and prognosis. Shinomiya et al. [9] reported the possible involvement of the DPI gene in the 13q34 amplicon. Adding to these reports, array-based CGH (aCGH) has been used recently in some reports. Ishizuka et al. [10] detected the amplification of eight oncogene loci using aCGH containing 51 oncogenes. Arai et al. [11] also detected the amplification of eight oncogene loci using the same aCGH. However, there are few reports studying whole genome profiles in ESC using aCGH. One of the most important factors in studies such as ours is the preparation of samples. It is necessary to obtain pure samples consisting only of cancer cells to evaluate the genomic changes associated with cancer. Samples obtained from cancer tissue usually contain both cancer cells and stromal cells. The recent development of laser microdissection (LMD) enables us to obtain pure cancer cells from the samples. In this study, we used a human bacterial artificial chromosome (BAC) array containing 4,030 human BAC clones and performed aCGH for 23 cases of ESC and two cases of normal esophageal epithelium. The BAC array does not require moving averages and data manipulation, in comparison with that of an oligoarray, because large clone inserts provide a strong signal as a result of the length of DNA available. In addition, the same clone on the array can be used for other applications, for example, fluorescence in situ hybridization (FISH) [12]. DNA samples were prepared through laser microdissected sections. Consequently, we identified 68 BAC clones that were well associated with the prognosis of postoperative patients. For 22 of these 68 BAC clones, there was a significant difference in the Kaplan-Meier survival curve, using the log-rank test, when comparing patients who had an alteration in a particular clone with those who did not.
Patients Cancer tissue samples were taken from 27 patients with ESC who underwent esophagectomy with lymph node dissection or endoscopic biopsy between 1999 and 2003 at our hospital. The age of the patients ranged from 36 to 82 years (mean, 63.5 years). Written informed consent was obtained from each patient. We evaluated 23 of the 27 cases that passed a strict DNA quality examination. Twenty-one samples were operated cases and two samples were biopsy cases. Two of the 21 operated patients died of noncancer related causes 6 months and 7 months after the operation. The other 19 patients were re-evaluated >5 years after the initial operation. Ten of the patients died of cancer recurrence within 334 months (mean survival time, 10.1 months). One of the nine surviving patients had local recurrence. Using the tumor node metastasis (TNM) classification system of the International Union Against Cancer [13], the 19 patients who underwent an operation were classified as follows: one with T1 tumors, three with T2 tumors, 13 with T3 tumors, and one with T4 tumors. Pathologically, all the tumors were squamous cancer. Lymph node metastases were present in 16 of the 19 patients (84.2%).
Tissue Samples and Preparation of the Cancer Cell Population by LMD LMD was performed as shown in Figure 1. Ten to 15 slides were dissected for each case. Sections were microdissected using the LMD system with a 337 nm nitrogen UV laser (Leica Laser Microdissection System, Leica Microsystems). Target cells were dropped immediately into a microcentrifuge tube cap filled with 30 µl of lysis buffer (Qiagen, Hiden, Germany). At least 10,000 cancer cells were collected from each sample, and then genomic DNA was extracted with the QIAamp DNA Micro Kit (Qiagen) according to the manufacturer's instructions. The DNA concentration was calculated using a NanoDrop (ND-1000; NanoDrop Technologies, Wilmington, DE) according to the manufacturer's instructions. We were able to extract 5003,000 ng genomic DNA from each sample.
Microarray Hybridization and Scanning The MACarray Karyo 4000 array slide (Macrogen Inc., Seoul, Korea, http://www.macrogen.com) was used for this study. It contains 8,060 elements, representing 4,030 human BAC clones selected from the original BAC library (96,768 clones) established by Macrogen Inc. in duplicate and randomly distributed throughout one microarray slide. The BACs, 4,030 clones, employed in this array were selected as fulfilling the following criteria: (a) the end-sequencing of the clones was complete and (b) the chromosomal locations of the clones were determined by FISH physical mapping. The entire genome was covered by BAC clones with a resolution of <1 M base pairs. Test DNA and human female DNA (0.5 µg each) were labeled with Cyanine3- or Cyanine5-dCTP (Perkin Elmer, Wellesley, MA) using the random primer method (BioPrime DNA labeling kit, Invitrogen, Carlsbad, CA). Labeled probes were purified by spin column (QiaQuick PCR Purification Column, Qiagen) and dissolved in 140 µl of hybridization solution (Macrogen) containing 100 µl of Cot-1 DNA solution (Macrogen) and 4 µl of yeast tRNA solution (Macrogen). The probe solution was heated to 70°C for 10 minutes to denature the DNA, then incubated for 60 minutes at 37°C to block repetitive sequences. The array slide was pretreated with salmon sperm DNA for 30 minutes at room temperature and then mounted on the slide processor (GeneMachines HybStation, Genomic Solutions, Ann Arbor, MI). After injection of 120 µl of probe, hybridization was performed for 72 hours on the HybStation with continuous agitation, followed by posthybridization washes as follows: 50% formamide/2x standard saline citrate (SSC) for 15 minutes at 46°C, 0.1% SDS/2x SSC for 30 minutes at 46°C, PN buffer (0.1M Na2PO4/0.1% NonDiet P-40, Nacalai Tesque, Kyoto, Japan) for 15 minutes at 46°C, and 2x SSC for 5 minutes at 46°C. After washing, the slides were dehydrated in ethanol. Following hybridization, the slides were scanned at 532 and 635 nm by a GenePix 4000A (Axon Instruments, Sunnyvale, CA). Next, MacViewer software (http://www.macrogen.com/eng/biochip/macviewer.jsp/), an analytical software program developed by Macrogen, was used to locate spots automatically on the Cy3 and Cy5 image acquisitions and to calculate fluorescence ratios. The MacViewer software automatically analyzed and summarized the results as follows: (a) averaged the ratios of the replicates and calculated the standard deviation, (b) rejected individual spot data based on several criteria (including weak fluorescent signals), (c) adjusted the Cy5/Cy3 ratios such that ratios of the normal genomic regions were always equal to 0, despite variations in dye labeling efficiency, and (d) plotted data relative to the position of the clones in the human genome (according to July 2003 University of CaliforniaSanta Cruz cartography). Data from the aCGH are summarized in online supplementary Table 1.
Statistical Analysis
aCGH of Normal Esophageal Squamous Epithelium: Quality Controls To confirm the profile of normal epithelium, we performed aCGH of normal epithelial DNA for two cases (case 14 and 24). As shown in Figure 2 (cases 14N and 14T), few alterations were seen in genomic DNA extracted from normal esophageal squamous cells (case 14N) compared with those in the cancer cell sample (case 14T). In cases 14N and 24N, only 18 (0.4%) and 31 (1%) loci of the 4,030 BAC clones deviated from the ±0.5 range of the log2 ratio, respectively. These regions were omitted for further analyses in this study.
BAC Clones Showing a Nonrandom Gain or Loss in Cancer Samples Figure 2 shows log2 ratio plots for the 4,030 BAC clones from four (cases 14T, 6T, 12T, and 28T) of the 23 ESC cases. Gains in 3q, 8q, and 11q and losses in 3p, 4p, 9p, and 13q have been reported in previous cCGH studies, and we also verified those alterations using our methods. Table 1A and 1B list the BAC clones showing a gain (>0.5 of the log2 ratio) or a loss (<0.5 of the log2 ratio), according to the number of cases with a gain or a loss. The clones containing known oncogenes, such as CCND1, MYC, and PIK3CA, showed gains in 15, 14, and 12 cases, respectively. In addition, the clones containing the known tumor suppressor gene FHIT showed losses in 15 cases. The BAC clones showing a higher gain (>1.0 of the log2 ratio) are indicated with an asterisk in Table 1
Correlation Between a Copy Number Gain or Loss and Patient Prognosis In 19 patients who underwent a curative operation and were followed for >5 years after the operation, nine survived and 10 died from recurrent disease (nonsurvivors). To elucidate the clones that correlate with postoperative prognosis, we counted the number of cases representing a gain or loss (deviation from the log2 ratio of ±0.5) in each of the 4,030 BAC clones, and with this number, a 2 test was performed between the nine survivors and the 10 nonsurvivors. A significant difference (p < .05) was detected in 92 of the 4,030 clones. Twenty-four of these 92 clones were omitted because they contain a region known to have genomic variants, according to the Database of Genomic Variants, so 68 clones were used. Two clones showed a normal ratio in survivors but a gain in nonsurvivors, 25 showed a normal ratio in survivors but a loss in nonsurvivors, nine showed a gain in survivors but a normal ratio in nonsurvivors, and 32 showed a loss in survivors but a normal ratio in nonsurvivors (Table 2). Almost half of these clones (31 of the 68 BAC clones) were located at nine limiting regions: 5.3 Mb of 4q13.34q21.21 (76.481.7 Mb, 5 clones), 2.5 Mb of 4q25 (109.4111.9 Mb, 3 clones), 4.1 Mb of 4q264q27 (119.9124.0 Mb, 3 clones), 5.1 Mb of 13q14.1113 (40.045.3 Mb, 3 clones), 0.5 Mb of 13q14.23 (49.6 Mb-50.2 Mb, 3 clones), 1.6 Mb of 20q12 (37.138.7 Mb, 3 clones), 1.9 Mb of 20q1213.12 (40.242.1 Mb, 3 clones), 3.3 Mb of Xq1213.1 (64.968.2 Mb, 4 clones), and 1.4 Mb of Xq21.33 (93.795.1 Mb, 4 clones). Furthermore, we plotted Kaplan-Meier survival curves for each of these 68 BAC clones, and used a log-rank test to distinguish between patients who had an alteration in a specific clone and those who did not. Consequently, 22 of these 68 BAC clones showed a significant difference (p < .05) using the log-rank test: one clone showed a normal ratio in survivors and a gain in nonsurvivors (Table 2A), nine showed a normal ratio in survivors and a loss in nonsurvivors (Table 2B), two showed a gain in survivors and a normal ratio in nonsurvivors (Table 2C), and 10 showed a loss in survivors and a normal ratio in nonsurvivors (Table 2D). Figure 4 shows the survival curves of six representative BAC clones of the 22 total. Although the number of cases in this study is small, we performed univariate and multivariate analyses with copy number changes in these six representative BAC clones and with the clinicopathological factors of each case. Univariate analysis showed that the following factors were significantly related to postoperative survival: lymphatic invasion, venous invasion, and copy number changes in each of the six BAC clones (online supplementary Table 2A). Multivariate regression analysis, using lymphatic invasion, venous invasion, and alteration in each of the six BAC clones, indicated that the copy number changes of four of the clones (BAC no. 164, 2622, 5282, and 5622) were independent prognostic factors (online supplementary Table 2
Previous cCGH studies of esophageal cancer identified several characteristic chromosomal aberrations, such as gains in 3q, 8q, and 11q and losses in 3p, 4p, 9p, and 13q [79, 18, 19]. These changes were similarly identified in our current study. The two previous aCGH studies disclosed that amplification of the regions of 11q13 containing the Cyclin D1, FGF3/4, and EMS1 genes were each identified in 9 of 32 cases [10] and in 5, 4, and 7 of 20 cases [11], respectively. In our study, we detected this amplification as a 1.1-Mb amplicon corresponding to six BAC clones in more cases (Table 1A and Fig. 3). The higher detectability in our study may be partly a result of strict sample preparation with LMD [20]. On the other hand, the remaining 8 of the 23 cases did not show a gain in this region, suggesting that this change may not always occur in ESC. Amplification at 11q13 has been reported previously not only in ESC but also in oral cancer, head and neck cancer, nasopharyngeal cancer, and several other types of cancer [2125]. For example, Hui et al. [21] reported this amplification as a 5.3-Mb amplicon at 11q13.113.3 as the most frequently detected gain in 26 nasopharyngeal cancers. They found that CCND1 is a target oncogene of nasopharyngeal cancer by analyzing six genesMEN1, CCND1, FGF3, EMS1, GARP, and PAK1with high-resolution aCGH. Other alterations in the clones at 17q12, containing ERBB2, or 8p11, containing FGFR1, were also identified in two and four cases, respectively, similar to previous reports [10]. However, the alterations in these clones were demonstrated as a part of widely gained regions, and thus are much different from the 11q13.3 amplifications, with a clear-cut, narrow change. Also, one clone (BAC no. 5359) with a novel loss (indicating homozygous deletion) was recognized at 16p13.3, containing the genes ZNF434, ZNF174, ZNF597, FLJ14154, LOC390671, and CLUAP1, which has not been reported previously in ESC (online supplementary Table 1). In addition, one clone (BAC no. 5308) with a novel loss was recognized at 11p15.4 (online supplementary Table 1) containing the gene CDKN1C, known as a tumor suppressor gene [2629]. However this clone was omitted from our analysis because it contains a region known to have genomic variants according to the Database of Genomic Variants.
With respect to the relationship between patient prognosis and chromosomal gain or loss, cCGH analysis has previously demonstrated that gains in 5p, 7q, and 12p and losses in 4p, 9p, and 11q are associated with poor prognosis in ESC patients [8, 18]. However, in our study, no such correlations were demonstrated using aCGH analysis. A chromosomal gain in 7q and losses in 4p and 11q were observed in both survivors and nonsurvivors. However, there were no significant differences. There were no chromosomal gains or losses in 5p, 12p, or 9p. Next, we studied the correlation between patient prognosis and DNA copy number change for each clone. Using a Many reports have described that DNA copy number alterations are usual phenomena in poor prognosis patients with many kinds of cancer, including esophageal cancer [8] [18]. However, this study clarified that alterations are frequently seen in patients with a good prognosis, compared with those with a poor prognosis, in a large number of clones. This is a novel finding, and needs further study to disclose its significance. Amplifications and deletions in genomic DNA are associated with high and low expression levels of the corresponding mRNA, respectively, for the most part [37, 38]. Prediction of the prognosis of cancer patients can be performed with DNA examination, as well as with RNA expression examination, protein expression examination, and so on. This study clearly demonstrates that some selected loci of BAC clones are very useful to predict the prognosis of patients with ESC. Studies using RNA expression require more strict sample preparation than those using DNA, thus, we consider DNA examination to be more practical than the use of RNA expression. Indeed Kyomoto et al. [39] reported that genomic amplification of the region that contains the Cyclin D1 gene is a more potent prognostic factor than its protein overexpression in human head and neck squamous cancer. It is desirable to develop methods that can easily identify gains or losses in specific genomic regions for practical examination. To our knowledge, our study is the first to demonstrate the correlation between alterations in specific loci of BAC clones and patient prognosis using aCGH. This study is very precise, because we used LMD samples obtained from frozen tissue. Of course, LMD is difficult to use in routine examination in its present form; however, if much easier methods were developed, this would be very useful in practice. In addition, the BAC clones identified in this study give us the opportunity to determine which specific genes are associated with cancer progression or metastasis.
The authors indicate no potential conflicts of interest.
We thank Y. Nakagawa for her technical assistance. Financial support came from Grant-in-Aid for Scientific Research, Japan Society for the Promotion of Science, No. 17109013 (S), and Ministry of Education, Culture, Sports, Science and Technology, No. 18790964(B). The authors Masaki Mori and Takashi Hirano contributed equally.
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||