
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Human Cancer Biology |
Authors' Affiliations: 1 Department of Microbiology, College of Medicine, Catholic University of Korea, Socho-gu, Seoul; 2 Korea National Cancer Center, Research Institute, Division of Cancer Control and Epidemiology, Gyeonggi-do; Departments of 3 Pathology and 4 Thoracic and Cardiovascular Surgery, College of Medicine, Dankook University, Cheonan, Chungnam, Republic of Korea; and 5 The Wellcome Trust Sanger Institute, Hinxton, Cambridge, United Kingdom
Requests for reprints: Yeun-Jun Chung, Department of Microbiology, College of Medicine, Catholic University of Korea, 505 Banpo-dong, Socho-gu, Seoul 137-701, Republic of Korea. Phone: 82-2590-1214; Fax: 82-2596-8969; E-mail: yejun{at}catholic.ac.kr.
| Abstract |
|---|
|
|
|---|
Experimental Design: One-megabase resolution array comparative genomic hybridization was applied to 29 squamous cell carcinomas and 21 adenocarcinomas of the lung. Tumor and normal tissues were microdissected and the extracted DNA was used directly for hybridization without genomic amplification. The recurrent genomic alterations were analyzed for their association with the clinicopathologic features of lung cancer.
Results: Overall, 36 amplicons, 3 homozygous deletions, and 17 minimally altered regions common to many lung cancers were identified. Among them, genomic changes on 13q21, 1p32, Xq, and Yp were found to be significantly associated with clinical features such as age, stage, and disease recurrence. Kaplan-Meier survival analysis revealed that genomic changes on 10p, 16q, 9p, 13q, 6p21, and 19q13 were associated with poor survival. Multivariate analysis showed that alterations on 6p21, 7p, 9q, and 9p remained as independent predictors of poor outcome. In addition, significant correlations were observed for three pairs of minimally altered regions (19q13 and 6p21, 19p13 and 19q13, and 8p12 and 8q11), which indicated their possible collaborative roles.
Conclusions: These results show that our approach is robust for high-resolution mapping of genomic alterations. The novel genomic alterations identified in this study, along with their clinicopathologic implications, would be useful to elucidate the molecular mechanisms of lung cancer and to identify reliable biomarkers for clinical application.
Some genomic aberrations in tumors have been suggested to be prognostic markers or can be used to identify the target genes for treatment or prevention (3, 4). Likewise, in other solid tumors, chromosomal aberrations are thought to be critical molecular events in the pathogenesis of lung cancer (5, 6). However, clinically applicable screening tools or prognostic markers are still underdeveloped. Because the lack of efficient screening methods and therapy accounts for the poor outcome of lung cancer, genome-wide assessment of aberrations could help in developing more accurate diagnostic and therapeutic strategies.
For this reason, previous cytogenetic studies using conventional comparative genomic hybridization (CGH) or fluorescence in situ hybridization have focused on identifying the chromosomal aberrations associated with NSCLC. Recurrent genomic alterations have been observed in NSCLC, including the gains of partial or whole chromosomal arms on 1q, 3q, 5p, and 8q along with the losses on 3p, 6q, 8p, 9p, 13q, and 17p (711). However, the
10 Mb resolution of conventional CGH is insufficient for the precise identification of submicroscopic changes (12). As accumulating evidence suggests that changes in the genomic dosage contribute to tumorigenesis by altering the expression levels of the cancer-related genes (13, 14), more detailed analyses with sufficient resolution are required.
For enhancing the resolution, array CGH using mapped bacterial or P1 artificial chromosomes (BAC/PAC) rather than metaphase chromosomes, has been recently developed (1517). This technique provides a high resolution that is directly related to the genomic density and insert size of the arrayed clones. Array CGH has emerged as a useful tool for detecting and mapping the genomic aberrations, which may contain putative oncogenes or tumor suppressor genes and for performing a molecular classification of tumors (18).
To see genomic alterations and their clinicopathologic implications in NSCLC, we applied genome-wide array CGH to the genomic DNA extracted from the microdissected tissues of 29 squamous cell carcinoma and 21 adenocarcinoma cases, on which the association study was done. Using this strategy, the genomic copy number changes specific to NSCLC including novel minimally altered regions (MAR) were identified. Those genomic alterations are likely to be related to tumorigenesis or the clinical outcomes of lung cancer.
| Materials and Methods |
|---|
|
|
|---|
Tissue preparation. After surgical resection, tumor and adjacent normal tissues from the same patient were collected separately and snap-frozen in a deep freezer. Frozen sections were prepared of 10 µm thickness on a gelatin-coated slide using 2800 Frigocut (Reighert-Jung, Germany). After H&E staining, tumor cellrich area (>60% of tumor cells) and histologically normal cell area were selected under the microscope and dissected manually. Microdissected tissues were transferred into the cell lysis buffer (1% proteinase-K in TE buffer) and DNA was extracted. DNA from normal tissue was used as reference DNA for array CGH. Extracted DNA was purified using a DNA purification Kit (Solgent, Daejeon, Korea) and used for dye labeling reactions.
Array comparative genomic hybridization and image analysis. We used human large insert clone arrays with 1 Mb resolution across the whole genome printed by the Sanger Institute Microarray Facility (19). DNA labeling, prehybridization, hybridization, and posthybridization processes were done as described previously (19, 20). Arrays were scanned using GenePix 4100A scanner (Axon Instruments, Union City, CA) and the image was processed using GenePix Pro 6.0.
Data processing, normalization, and mapping of BAC clones. Normalization and re-aligning raw array CGH data were done using the web-based array CGH analysis interface, ArrayCyGHt (http://genomics.catholic.ac.kr/arrayCGH/; ref. 21). Mapping of large insert clones was done according to the genomic location in the UCSC genome browser (May 2004 freeze). In total, 2,987 successfully mapped BAC clones out of initial 3,014 clones were processed subsequently. All the genomic coordinates such as cytogenetic bands or gene positions described in this study are based on the same version of the human genome available on the UCSC genome browser.
Data analysis for chromosomal alterations. To set the cutoff value for chromosomal alterations of individual large insert clones, we did a series of four independent normal hybridizations (three sex-matched and one male versus female hybridizations) as a control. The average SD value of the control batch was 0.081. Adopting the criteria of a previous study (22), the cutoff value for the copy number aberrations was set to be ±0.2 in log2 ratio in this study, >2-fold of control SD. The entire chromosome arm gain or loss was determined as previously described (23). Regional copy number change was defined as DNA copy number alteration limited to part of a chromosome. High-level amplification of clones was defined when their intensity ratios were >1.0 in log2 scale, and vice versa for homozygous deletion. The boundary of copy number change was assigned to be halfway between the two neighboring clones.
Definition of minimally altered regions. To define MARs of chromosomal gain or loss, we used CGH-Miner (http://www-stat.stanford.edu/
wp57/CGH-Miner/) to smooth the raw intensity ratio and to identify the breakpoints of chromosomal alterations (24). A series of four normal hybridizations were combined as a control and the analysis was done with recommended program variables. The significant gains or losses reported by the program were directly used for subsequent aligning procedures. Minimal regions of chromosomal gains and losses were determined by altered segments recurring for at least seven samples.
Statistical analysis. The significance of the differences in chromosomal arm changes between squamous cell carcinomas and adenocarcinomas was tested by two-sided Fisher's exact test. The correlations between recurrent genetic changes on minimally altered regions were assessed using univariate pairwise Pearson's correlation. For multiple comparisons, the step-down Sidak method was used to adjust the overall level of significance. In this case, the pairs of genetic changes on the same chromosomal arm were excluded for the concordance analysis. The correlations between genetic alterations and clinical variables were analyzed by two-sided Fisher's exact test. All the MARs as well as chromosomal arm changes were included in the analysis. For comparison, four kinds of clinical variables were treated as categorical variables such as age (<60 versus
60 years), stage (stages I and II as early versus stages III and IV as advanced), lymph node status (negative versus positive), and the disease recurrence (presence versus absence). Kaplan-Meier method was used for survival analysis and the difference between survival curves was compared using the log-rank test in univariate model. To identify independent prognostic factors after adjusting clinical variables such as age, sex, stage, treatment, metastasis, and recurrence, Cox regression was done. In all statistical analyses, P < 0.05 was considered significant.
| Results |
|---|
|
|
|---|
|
|
|
|
|
Association between genomic aberrations and clinical characteristics. Four types of clinical variables (age, stage, lymph node, and recurrence) were analyzed for their association with the genomic alterations identified (see Supplementary Table S4). Significant associations were observed for the MAR-L on 13q21 with cancers from those aged <60 and an advanced stage (stages III and IV). The chromosomal gain of Xq was also found to be associated with an advanced stage and disease recurrence. The MAR-G on 1p32 and a chromosomal gain of Yp were associated with being lymph nodenegative.
Survival analysis was done to assess the prognostic values of the genetic aberrations identified. Using Kaplan-Meier methods, we identified that six genetic aberrations associated with a relatively poor survival (Fig. 4); gain of 10p (P = 0.0091) and 16q (P = 0.0262), loss of 9p (P = 0.0082) and 13q (P = 0.0019), and the MAR-Gs on 6p21 and 19q13 (P = 0.0265 and 0.0295, respectively). Multivariate analyses using all the genetic alterations identified, as well as the clinical variables such as age, gender, stage, treatment, metastasis, and recurrence showed that four genetic alterations and three clinical variables remained independent factors to be significantly associated with a poor survival outcome (Table 3). One of the four genetic alterations was a MAR-G on 6p21 and the other three were a loss of 9p, and gains of 7p and 9p.
|
|
| Discussion |
|---|
|
|
|---|
The frequent chromosomal changes in this study are largely consistent with previous cytogenetic analysis (711), including a loss of the Y chromosome in male patients (5, 6). It is notable that the copy number alterations on the small chromosomes such as 19, 20, and 22 were much more frequent in our study. This might be due to the differences in the analytic methods. However, it is more likely to reflect the potential of array CGH to improve the low resolution of conventional CGH as described elsewhere (25). The genomic size of the high copy number changes ranged from 0.31 to 14.78 Mb, and most of them were <5 Mb. Genomic changes <5 Mb are likely to be novel because they cannot be detected using conventional CGH (12).
To date, two array-based CGH analyses of lung cancer have been published (26, 27). Our results are in general agreement with the previous array CGH results. For example, the 30-Mb sized recurrently gained region around 3q26 reported by Massion et al. (26) overlapped with the MAR-G on 3q26-q28 in this study. Other metaphase CGH studies have also shown frequent amplifications around the same region in squamous cell carcinomas (911). However, some of the MARs identified were not consistent with them. Gain of 3q26 was the only recurrent alteration in Massion et al.'s study, whereas we identified not only 3q26 but also 16 more MARs across various chromosomes. The difference is thought to be due to the resolution. This study placed 2,987 large insert clones in 1 Mb intervals (average of 125 clones per chromosome), whereas they used 348 BAC clones (average of 15 clones per chromosome, at most 78 clones on chromosome 3). That might explain why more MARs could be defined. In Jiang et al.'s study (27), the recurrent genomic alterations were only partly compatible with this result, and there was no clear minimal recurrent gain on 3q. They used a cDNA microarray rather than a BAC/PAC array for CGH analysis. Despite its advantages, the sensitivity and reliability for detecting copy number changes, particularly single copy changes, is known to be limited (18).
Several interesting cancer-related genes are located in the genomic alteration regions identified in this study. For example, PIK3CA on 3q26-q28, which is one of the genes in the most common MAR-G in this study, is believed to contribute to the tumorigenesis of squamous cell carcinomas by involving in the phosphoinositide-3-OH kinase signaling pathway (26). The ECT2 oncogene, which is located in the same locus, is known to activate the Rho signaling pathway leading to a malignant transformation (28). In our unpublished data,6 ECT2 is frequently overexpressed in primary lung cancers. Although, there has been no report of ECT2 overexpression in lung cancer, this suggests the involvement of ECT2 in malignant transformation or the progression of lung cancer. Most high-level amplifications were usually observed in one or two cases with the exceptions of the amplifications on 3q. They might reflect the individual nature of the genomic evolution for the respective NSCLC cases. Among the putative cancer-related genes in the high-level amplification region (Table 1), the expression of the AF1Q, TPM3, REL, SKIL, ECT2, BCL6, MLLT6, YES1, and HKR genes have not been reported in lung cancer.
A homozygous deletion on 10q23.31 observed in one squamous cell carcinoma case contains the well-known tumor suppressor gene, PTEN. PTEN is known to encode lipid phosphatase, which negatively controls the signaling proteins activated in the phosphoinositide-3-OH kinase pathway (29). This suggests a potential role of the phosphoinositide-3-OH kinase signaling pathway in NSCLC pathogenesis.
In contrast to the high copy number changes that were largely limited to a few samples, single copy changes were found in many more samples, which is indicative of a shared mechanism common to the earlier stage of NSCLC. Minimal recurrent gains and losses were successfully identified using high-resolution array CGH. Seventeen MARs of various sizes were defined. The MAR-Gs on 1p, 2p, 6p, 8p, 19p, and 20p along with MAR-Ls on 5q and 20q are believed to be novel features in lung cancer, which shows the advantage of genome-wide, high-resolution mapping of the genomic alterations. Interestingly, three pairs of MAR-Gs (19q13.1 and 6p21, 19p13 and 19q13.1, and 8p12 and 8q11-12) showed significant correlations among themselves, suggesting a possible collaborative role in the tumorigenesis of NSCLC. Further investigations will be needed to confirm the functional consequences of the associations between the MARs. Some of the MARs showed significant correlations with the clinical features. This suggests that the common single copy changes identified by high-resolution analysis can be useful biomarkers for the clinical characteristics of lung cancer.
Survival analysis revealed that six genetic alterations were associated with a poor survival outcome in the univariate model (Fig. 4). Among those six alterations, a loss of 9p was reported to be associated with a poor survival outcome (30). However, there has been no report about the association between the other five genomic alterations and survival outcomes in lung cancer. These genomic alterations might be a novel genetic indicator of the prognosis of NSCLC after the appropriate validation. In particular, two of these alterations are MARs, which appeared concordantly (P = 0.0482). These two MARs, MAR-Gs on 6p21 and 19q13, contain cancer-related genes such as PIM1, CCND3 (both in 6p21), and HKR1 (19q13). The high expression level of HKR1 after administering platinum drugs has been reported to be associated with the acquisition of resistance to chemotherapy (31). There is no report demonstrating an alteration of CCND3 and PIM1 proto-oncogene in lung cancer. However, both genes are well known to be involved in the tumorigenesis pathways of various tumors. Therefore, further investigations will be needed to evaluate their specific implications in lung cancer.
Subsequent Cox regression analysis identified seven factors, including four genomic alterations such as MAR-G on 6p21, 9p loss, 7q gain, and 9q gain, to be independent indicators of a poor survival outcome. This indicates that in addition to the clinical factors, precisely defined recurrent genetic alterations can be useful biomarkers for the prognosis of NSCLC. However, due to the limited number of samples in this study, further studies with a larger sample size will be needed to confirm the prognostic implication of these genomic alterations and to identify further reliable prognostic markers.
This study showed that a well-designed high-resolution array CGH could define more novel regions possibly associated with the tumorigenesis of lung cancer. Therefore, these results will give a clue for further studies to elucidate lung cancer pathogenesis or to develop biomarkers for predicting the prognosis or treatment response of lung cancer.
| Acknowledgments |
|---|
| Footnotes |
|---|
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.
Note: T-M. Kim and S-H. Yim contributed equally to this paper.
Supplementary data for this article are available at Clinical Cancer Research Online (http://clincancerres.aacrjournals.org/).
6 M.S. Kwon, H.M. Kang, and Y.J. Chung, manuscript in preparation. ![]()
Received 5/31/05; revised 8/ 3/05; accepted 9/ 1/05.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
W. Wang, M. R. Spitz, H. Yang, C. Lu, D. J. Stewart, and X. Wu Genetic Variants in Cell Cycle Control Pathway Confer Susceptibility to Lung Cancer Clin. Cancer Res., October 1, 2007; 13(19): 5974 - 5981. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Li, N. W. Todd, Q. Qiu, T. Fan, R. Y. Zhao, W. H. Rodgers, H.-B. Fang, R. L. Katz, S. A. Stass, and F. Jiang Genetic Deletions in Sputum as Diagnostic Markers for Early Detection of Stage I Non-Small Cell Lung Cancer Clin. Cancer Res., January 15, 2007; 13(2): 482 - 487. [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 | Cell Growth & Differentiation |