|
Advertisement |
|
|
||
 |
|
|||||
|
|||||
|
 |
||||||
|
||||||
|
||||||
|
Journal of Clinical Oncology, Vol 22, No 19 (October 1), 2004: pp. 3965-3972 © 2004 American Society of Clinical Oncology. DOI: 10.1200/JCO.2004.01.094 Identification of Novel Prognosticators of Outcome in Squamous Cell Carcinoma of the Head and NeckFrom the Laboratory of Epithelial Cancer Biology, Head and Neck Service, Departments of Surgery, Epidemiology and Biostatistics, and Medicine, Memorial Sloan-Kettering Cancer Center, New York, NY; and Department of Pediatrics, Baylor College of Medicine, Houston, TX Address reprint requests to Bhuvanesh Singh, MD, Laboratory of Epithelial Cancer Biology, Head and Neck Service, Memorial Sloan-Kettering Cancer Center, 1275 York Ave, New York, NY 10021; e-mail: singhb{at}mskcc.org
PURPOSE: The goal of this study was to identify chromosomal aberrations associated with poor outcome in patients with head and neck squamous cell carcinoma (HNSCC).
PATIENTS AND METHODS: We assessed the global genomic composition of 82 HNSCCs from previously untreated patients with comparative genomic hybridization (CGH). The CGH data were subcategorized into individual cytogenetic bands. Only genomic aberrations occurring in more than 5% of cases were analyzed, and redundancies were eliminated. Each aberration was submitted to univariate analysis to assess its relationship with disease-specific survival (DSS). We used Monte Carlo simulations (MCS) to adjust P values for the log-rank approximate RESULTS: CGH analysis revealed a recurrent pattern of chromosomal aberrations typical for HNSCC. Univariate analysis revealed 38 abnormalities that were correlated with DSS. After controlling for multiple comparisons and confounding effects of stage, five chromosomal aberrations were significantly associated with outcome, including amplification at 11q13, gain of 12q24, and losses at 5q11, 6q14, and 21q11 (MCS adjusted P = .0009 to P = .01). CONCLUSION: HNSCC contains a complex pattern of chromosomal aberrations. A sequential approach to control for multiple comparisons and effect of confounding variables allows the identification of clinically relevant aberrations. The significance of each individual abnormality merits further consideration.
Head and neck squamous cell carcinoma (HNSCC) can arise anywhere along the upper aerodigestive tract and is etiologically associated with tobacco and alcohol exposure.1 Overall, HNSCC portend a grave clinical course, with 5-year survival rates less than 50%.2 Selection of optimal treatment for individual patients requires accurate assessment of malignant potential at presentation. Traditionally, prognostic stratification of HNSCC has relied on analysis of clinicopathologic variables, which are inconsistent predictors of clinical behavior.3 Because HNSCC results from a sequential accumulation of genetic abnormalities, recent focus has shifted toward prognostication based on molecular profiling.4,5 The advent of global genomic screening based on comparative genomic hybridization (CGH) has contributed significantly to the molecular characterization of HNSCC.6-9 CGH analysis of HNSCC has identified gains of chromosomal regions 1q21, 3q26.3, 5p15, 7p12, 8q24, 9q34, 11q13, and 20q12 and deletions of 3p, 4q, 5q, 7q22, 8p23, 9p21, 13q12-24, 17p, 18q21, 21q11-21, and 22q.10 Studies investigating the prognostic utility of individual alterations have suggested that gains of 3q26 and 11q13 and deletions of 8p23 and 22q may be valuable markers of aggressive disease.11,12 However, the density of the information resulting from CGH analysis has limited attempts at defining individual abnormalities associated with outcome. The simultaneous analysis of multiple genetic and clinicopathologic factors is limited by several factors, especially issues pertaining to multiplicity of testing, with a significant possibility for inaccurate correlations.13-15 In the present study, we used a novel statistical approach to simultaneously assess the prognostic utility of CGH-detected DNA copy number changes in a large panel of HNSCC while controlling for multiple correlations. This study identified gains of 11q13 and 12q24 and deletions of 5q11-15, 6q14-21 and 21q11-21 as independent prognostic factors in patients with HNSCC. The approach presented here is applicable to the analysis of large-scale molecular-cytogenetic profiling of cancer in general.
Tissue Acquisition and Clinical Information Randomly selected tumor samples from patients undergoing curative treatment for HNSCC were obtained at the time of biopsy or surgical resection in the operating room following guidelines established by the institutional review board. Only previously untreated cases with biopsy-proven, primary squamous cell carcinomas involving the oral cavity or laryngopharyngeal complex were included in the study. Clinical information was collected independent of the laboratory data from medical records, and all identifier data were terminally coded to maintain patient anonymity. The median follow-up time was 22.5 months overall and 29 months (range, 2 to 249 months) for survivors. Table 1 lists the clinical information characterizing these cases.
Comparative Genomic Hybridization All tumor specimens were examined histologically and confirmed to contain greater than 70% tumor. Test DNA was extracted from tumor specimens and reference DNA was extracted from normal placenta from a healthy donor as described previously. CGH was performed as described previously.6,10,16 In brief, tumor DNA from each case was labeled with fluorescein-12-dUTP (NEN Dupont, Boston, MA) and reference DNA from normal human placenta was labeled with Texas red 5-dUTP using nick translation (Gibco BRL, Rockville, MD). After competitive hybridization of equal amounts (2 µg) of test and reference DNA to metaphase spreads from normal lymphocytes, seven to 10 separate metaphases were captured and processed using the Quantitative Image Processing System (Quips Pathvysion system; Applied Imaging, Santa Clara, CA). Red, green, and blue fluorescence intensities were analyzed for all metaphase spreads, normalized to a standard length, and statistically combined to show the red to green signal ratio and 95% CIs for the entire chromosome. Copy number changes were detected based on the variance of the red to green ratio profile from the standard of 1. Ratio values of 1.2 and 2.0 were defined as thresholds for gains and amplifications, respectively, and losses were defined as ratio of 0.8 or less.
Prestatistical Data Analysis
Statistical Analysis
In our analysis, we performed a sequential series of steps to control for issues present in the analysis of CGH data. First, we eliminated rare chromosomal abnormalities found in less than 5% of the sample because any abnormality with low prevalence is likely not to be useful as a prognostic factor. Second, because redundant data do not provide additional prognostic information, we selected a single representative of each set of closely linked bands whose abnormalities co-occurred across patients in the data set (corresponding to the minimal common region of involvement). Then, we submitted the reduced list of chromosomal abnormalities to univariate analysis to assess the strength of relationship with DSS time using the log-rank test. We adjusted the P values by doing a large number of Monte Carlo simulations because the
Further, we used Hochberg-Benjamini's13 false discovery rate (FDR) controlling procedure to adjust for multiple testing of the large number of abnormalities when selecting a subset of the most significant abnormalities to submit to stepwise Cox proportional hazards regression, adjusting for significant demographic/clinical variables. The FDR multiplies the i-th smallest P value by a factor m/i, where m is the total number of significant tests. The FDR procedure starts with the largest nominal P value and stops at the first adjusted P value. Hochherg and Benjamini's approach is more powerful than the very conservative Bonferroni method. All abnormalities whose univariate tests resulted in Hochberg-Benjamini's procedure adjusted P value of less than .15 were considered in the stepwise Cox proportional hazards regression model. We again adjusted the P values from the
CGH Analysis Application of CGH to the analysis of the genomic composition of 82 cases of HNSCC demonstrated a consistent pattern of DNA copy number changes (Fig 1). Gain of genomic loci most commonly involved 3q26-27, 5p14-15.3, 7p12-22, 8q24, 9q34, 11q13, 17q25, 19, and 22q13. Deletions of genomic material most commonly involved 3p, 4q31, 5q11-15, 9p13-21, 11q14-25, 13q31, and 18q21-23. In addition, recurrent high-level amplifications were identified at 2q32, 3q26, 4p15.3 to 16, 5p15, 7q11.2-p12, 7q21, 8p11, 8q24, 9p22-24, 11q13, 12p13, 18p, and 19p.
Prestatistical Analysis We extended the assessment of DNA copy number changes across our samples of HNSCC by analyzing each cytogenetic band individually. This analysis showed a total of 595 DNA copy number abnormalities (gains, losses, amplifications) of cytogenetic bands (273 losses, 253 gains/amplifications, and 69 amplifications) that were found in at least one patient's sample (Table 2). However, because this approach is associated with the introduction of redundancy that may influence the statistical power (see Patients and Methods), we chose to select one abnormality among a set of closely linked (neighboring) DNA copy number changes of chromosomal bands (minimal common region). In addition, to reduce noise and improve statistical power, we excluded all abnormalities that occurred in less than 5% of cases. After these modifications, we found 207 copy number abnormalities in at least four patients, including 104 losses, 95 gains/amplifications, and eight amplifications (Table 2).
Outcome Analysis Thirty-eight abnormalities were associated with DSS at nominal P < .05 (Fig 2). These included five at P < .0001 (gain at 12q24 and 14q11; losses at 5q11, 6q14, and 21q11). Gain at 8p11 and amplification of 11q13 were among seven nominally significant at .0001 < P < .001. All abnormalities significant at P < .05 by the log-rank statistic  2 approximation remained at P < .05 in the Monte Carlo permutation simulations. Twenty abnormalities were significant at P < .015 in the Monte Carlo simulations and remained significant at P < .15 after the Hochberg-Benjamini adjustment for multiple testing. The identified prognostic abnormalities and tumor-node-metastasis system stage (which was the only univariate prognostic demographic/clinical characteristics associated with DSS) were then submitted to stepwise Cox proportional hazards regression (Table 3). Five (12q24, 5q11, 6q14, 21q11, and 11q13) were jointly significant at P < .05 when adjusted for tumor-node-metastasis system stage in the final regression equation (Table 4). CIs for regression coefficients were not adjusted for multiple comparisons. We did further Monte Carlo simulations that adjusted the joint significance levels based on the 2 statistics for the regression coefficients (Table 4).
No significant differences in clinicopathologic features were observed based on the presence or absence of one or more of the identified prognostic abnormalities (Table 5), except sex, which was not a prognostic factor of DSS. Of the 23 patients who had one or more prognostic genetic abnormalities, 87% (20 of 23 patients) were male (P = .009). To graphically show the effect of the five selected genetic abnormalities, we divided patients post hoc into three cohorts: stage I to III patients with no abnormality, stage IV patients with no abnormality, and patients with one or more of the five prognostic chromosomal abnormalities. The survival plots for these three are shown in Fig 3. Our data confirm association between 11q13 overrepresentation and disease-specific survival of HNSCC. Moreover, prognostic significance of 12q24 gain and deletions of 5q11-15, 6q14-21, and 21q11-21 has not been previously described and may aid the research of treatment planning of HNSCC.
HNSCC development and progression is driven by an accumulation of genetic events.21 The advent of global genetic screening methodologies, such as CGH, has allowed an improved understanding of the genomic profile of HNSCC, identifying recurrent gains, losses, and amplifications on virtually all chromosomal arms.10 Although the overall genomic profile of HNSCC is highly consistent, individual cases differ significantly in the types of accumulated genetic abnormalities, likely reflecting the heterogeneous clinical phenotype of HNSCC. Accordingly, previous studies have associated individual genomic aberrations with clinicopathologic factors, including lymphatic metastasis, tumor grade, and presence or absence of etiologic factors, including tobacco, alcohol, and human papillomavirus exposure.22-25 Despite these efforts, few studies have assessed the global DNA copy number profile of HNSCC in relationship to its clinical outcome. To date, only one group has associated CGH-detected DNA copy number abnormalities with clinical outcome in a large number of patients with HNSCC.11 This study correlated gains of 3q26 and 11q13 and deletion of 8p21-22 with survival outcomes. However, because of limitations in the statistical approach, especially the lack of correction for multiple testing, the results are difficult to interpret. In this study, we have applied a rigorous statistical approach to simultaneously assess the impact of multiple CGH-detected events on clinical outcome of a large panel of randomly ascertained HNSCCs. Application of CGH to the genetic analysis of 82 HNSCCs confirmed previous studies showing the presence of a recurrent pattern of DNA copy number changes in these tumors. To objectify the comparison of DNA copy number abnormalities that may vary in size and location between individual tumor cases, we chose to analyze the genomic constitution of our cases for each cytogenetic band. Subsequently, we reduced redundancy introduced by this approach by selecting a single locus representative (minimal common region) of closely linked (neighboring) abnormalities for outcome analyses. In addition, we attempted to reduce noise by excluding abnormalities present in less than 5% of cases, because these were unlikely to have an impact on outcome.
We evaluated 207 individual copy number abnormalities for univariate association with outcome, revealing 38 abnormalities that correlated with outcome at P < .05. Given the large number of genomic abnormalities tested, several of the loci linked with outcome likely represent false-positive results.13 We performed Monte Carlo permutation simulations to adjust for inaccuracies of the The association between 11q13 overrepresentation and outcome supports the validity of the statistical approach presented here. 11q13 amplification is strongly linked with poor clinical outcome of squamous cell carcinomas of the upper aerodigestive tract.11,26,27 Cyclin D1, a cell cycle regulator whose overexpression affects cellular survival, is a candidate oncogene within 11q13 locus.28 Antisense cyclin D1 treatment inhibits HNSCC cell proliferation, induces apoptosis, and leads to tumor shrinkage, supporting its role in HNSCC pathogenesis.29 Additionally, recent studies suggest that EMS1 and TAOS1 may drive selection for 11q13 amplification.30,31 Similar to 11q13 amplification, 12q24 amplification and deletions of 5q11-15, 6q14-21 and 21q11-21 have been described in squamous cell carcinomas arising in various organs including the head and neck, anogenital tract, and lung.32,33 Recently, Fitzsimmons et al34 reported evidence for the presence of an unknown mortality gene driving selection for 6q14 deletion in advanced tongue carcinomas. However, prognostic significance of the abnormalities has not been assessed. Moreover, gene targets driving selection for these abnormalities remain ill-defined. It is of significance that our previous work using gene array analysis identified several genes within the identified regions whose dysregulation is correlated with outcome of HNSCC.35 These genes include DHFR (5q11-13), MAP3K7 (6q14), EST H75853 (6q14), SCAD (12q24), IFNAR1 (21q22), and NCAM2 (21q21). The role of these genes in HNSCC requires further delineation. Bockmuhl et al11,36 suggested 3q26.3 amplification and 8p23 deletion as independent prognostic markers of HNSCC. In the present study, these abnormalities were commonly detected, but no outcome associations were identified. To assess the role of 3q26.3 copy number in more detail, we defined precise 3q26.3 copy numbers in an independent panel of 50 HNSCCs by fluorescence in situ hybridization using a novel probe.37 Although outcome was linearly associated with 3q26.3 copy number, only high copy numbers of 3q26.3 were independently associated with survival. Because 3q amplification seems to occur exclusively in patients with a history of tobacco/alcohol use,25 the findings of our study may have been diluted by the inclusion of a significant number of nonsmokers. Indeed, high-level 3q amplification was a significant predictor of outcome in the subset of smokers with HNSCC. With regard to the generalizability of our findings, several limitations should be taken into account. First, although the present study includes one of the largest cohorts of HNSCC that have been assessed by CGH to date, it is clear that our findings are based on a relatively small number of cases. Second, it needs to be taken into account that our study population was skewed toward cases involving the laryngopharyngeal complex as a result of specimen availability. Nonetheless, this report represents the first study on the impact of CGH-detected DNA copy number changes in HNSCC analyzed using a rigorous statistical approach. The statistical approach outlined here may prove valuable in the analysis of multiple biomedical factors in cancer, such as microarray-based analyses. In addition to 11q13 amplification, we propose amplification of 12q24 and deletions of 5q11-15, 6q14-21, and 21q11-21 as predictors of HNSCC outcome. Akin to 11q13 amplification, these alterations may also have prognostic value in squamous cell carcinomas arising in other anatomic loci, and merit further study.
The authors indicated no potential conflicts of interest.
Supported in part by the Young Investigator Award from the American Society of Clinical Oncology and the George H.A. Clowes, Jr, MD, FACS, Memorial Research Career Development Award from the American College of Surgeons. V.B.W. and W.S. contributed equally to this work. Authors' disclosures of potential conflicts of interest are found at the end of this article.
1. Maier H, Dietz A, Gewelke U, et al: Tobacco and alcohol and the risk of head and neck cancer. Clin Investig 70:320-327, 1992[Medline]
2. Shah JP, Lydiatt W: Treatment of cancer of the head and neck. CA Cancer J Clin 45:352-368, 1995 3. Lydiatt WM, Shah JP, Hoffman HT: AJCC stage groupings for head and neck cancer: Should we look at alternatives? A report of the Head and Neck Sites Task Force. Head Neck 23:607-612, 2001[CrossRef][Medline] 4. Sidransky D: Molecular genetics of head and neck cancer. Curr Opin Oncol 7:229-233, 1995[Medline]
5. Estilo CL, O-Charoenrat P, Ngai I, et al: The role of novel oncogenes squamous cell carcinoma-related oncogene and phosphatidylinositol 3-kinase p110alpha in squamous cell carcinoma of the oral tongue. Clin Cancer Res 9:2300-2306, 2003
6. Kallioniemi A, Kallioniemi OP, Sudar D, et al: Comparative genomic hybridization for molecular cytogenetic analysis of solid tumors. Science 258:818-821, 1992
7. Speicher MR, Howe C, Crotty P, et al: Comparative genomic hybridization detects novel deletions and amplifications in head and neck squamous cell carcinomas. Cancer Res 55:1010-1013, 1995
8. Redon R, Muller D, Caulee K, et al: A simple specific pattern of chromosomal aberrations at early stages of head and neck squamous cell carcinomas: PIK3CA but not p63 gene as a likely target of 3q26-qter gains. Cancer Res 61:4122-4129, 2001
9. Singh B, Reddy PG, Goberdhan A, et al: p53 regulates cell survival by inhibiting PIK3CA in squamous cell carcinomas. Genes Dev 16:984-993, 2002
10. Singh B, Gogineni SK, Sacks PG, et al: Molecular cytogenetic characterization of head and neck squamous cell carcinoma and refinement of 3q amplification. Cancer Res 61:4506-4513, 2001
11. Bockmuhl U, Schluns K, Kuchler I, et al: Genetic imbalances with impact on survival in head and neck cancer patients. Am J Pathol 157:369-375, 2000 12. Ashman JN, Patmore HS, Condon LT, et al: Prognostic value of genomic alterations in head and neck squamous cell carcinoma detected by comparative genomic hybridization. Br J Cancer 89:864-869, 2003[CrossRef][Medline] 13. Benjamini Y, Hochburg Y: Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Stat Soc B 57:289-300, 1995
14. Reiner A, Yekutieli D, Benjamini Y: Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics 19:368-375, 2003
15. Jain AN, Chin K, Borresen-Dale AL, et al: Quantitative analysis of chromosomal CGH in human breast tumors associates copy number abnormalities with p53 status and patient survival. Proc Natl Acad Sci U S A 98:7952-7957, 2001
16. Wreesmann VB, Ghossein RA, Patel SG, et al: Genome-wide appraisal of thyroid cancer progression. Am J Pathol 161:1549-1556, 2002
17. Wessels LF, van Welsem T, Hart AA, et al: Molecular classification of breast carcinomas by comparative genomic hybridization: A specific somatic genetic profile for BRCA1 tumors. Cancer Res 62:7110-7117, 2002 18. An International System for Human Cytogenetic Nomenclature (1985) ISCN 1985: Report of the Standing Committee on Human Cytogenetic Nomenclature. Birth Defects Orig Artic Ser 21:1–117, 1985 19. Kallioniemi OP, Kallioniemi A, Piper J, et al: Optimizing comparative genomic hybridization for analysis of DNA sequence copy number changes in solid tumors. Genes Chromosomes Cancer 10:231-243, 1994[Medline] 20. Cahill DP, Kinzler KW, Vogelstein B, et al: Genetic instability and Darwinian selection in tumours. Trends Cell Biol 9:M57-60, 1999[CrossRef][Medline]
21. Califano J, van der Riet P, Westra W, et al: Genetic progression model for head and neck cancer: Implications for field cancerization. Cancer Res 56:2488-2492, 1996
22. Bockmuhl U, Petersen S, Schmidt S, et al: Patterns of chromosomal alterations in metastasizing and nonmetastasizing primary head and neck carcinomas. Cancer Res 57:5213-5216, 1997 23. Bockmuhl U, Wolf G, Schmidt S, et al: Genomic alterations associated with malignancy in head and neck cancer. Head Neck 20:145-151, 1998[CrossRef][Medline]
24. Kutler D, Wreesmann VB, Goberdhan A, et al: Human papilloma virus DNA and p53 polymorphisms in squamous cell carcinomas from Fanconi anemia patients. J Natl Cancer Inst 95:1718-1722, 2003 25. Singh B, Wreesmann VB, Pfister D, et al: Chromosomal aberrations in patients with head and neck squamous cell carcinoma do not vary based on severity of tobacco/alcohol exposure. BMC Genet 3:22, 2002[CrossRef][Medline] 26. Shinozaki H, Ozawa S, Ando N, et al: Cyclin D1 amplification as a new predictive classification for squamous cell carcinoma of the esophagus, adding gene information. Clin Cancer Res 2:1155-1161, 1996[Abstract]
27. Bellacosa A, Almadori G, Cavallo S, et al: Cyclin D1 gene amplification in human laryngeal squamous cell carcinomas: Prognostic significance and clinical implications. Clin Cancer Res 2:175-180, 1996 28. Donnellan R, Chetty R: Cyclin D1 and human neoplasia. Mol Pathol 51:1-7, 1998[Abstract]
29. Sauter ER, Nesbit M, Litwin S, et al: Antisense cyclin D1 induces apoptosis and tumor shrinkage in human squamous carcinomas. Cancer Res 59:4876-4881, 1999 30. Brookes S, Lammie GA, Schuuring E, et al: Amplified region of chromosome band 11q13 in breast and squamous cell carcinomas encompasses three CpG islands telomeric of FGF3, including the expressed gene EMS1. Genes Chromosomes Cancer 6:222-231, 1993[Medline]
31. Huang X, Gollin SM, Raja S, et al: High-resolution mapping of the 11q13 amplicon and identification of a gene, TAOS1, that is amplified and overexpressed in oral cancer cells. Proc Natl Acad Sci U S A 99:11369-11374, 2002 32. Harris CP, Lu XY, Narayan G, et al: Comprehensive molecular cytogenetic characterization of cervical cancer cell lines. Genes Chromosomes Cancer 36:233-241, 2003[CrossRef][Medline] 33. Balsara BR, Testa JR: Chromosomal imbalances in human lung cancer. Oncogene 21:6877-6883, 2002[CrossRef][Medline] 34. Fitzsimmons SA, Ireland H, Barr NI, et al: Human squamous cell carcinomas lose a mortality gene from chromosome 6q14.3 to q15. Oncogene 22:1737-1746, 2003[CrossRef][Medline]
35. Belbin TJ, Singh B, Barber I, et al: Molecular classification of head and neck squamous cell carcinoma using cDNA microarrays. Cancer Res 62:1184-1190, 2002 36. Bockmuhl U, Ishwad CS, Ferrell RE, et al: Association of 8p23 deletions with poor survival in head and neck cancer. Otolaryngol Head Neck Surg 124:451-455, 2001[CrossRef][Medline]
37. Singh B, Stoffel A, Gogineni S, et al: Amplification of the 3q26.3 locus is associated with progression to invasive cancer and is a negative prognostic factor in head and neck squamous cell carcinomas. Am J Pathol 161:365-371, 2002 Submitted January 14, 2004; accepted June 30, 2004. This article has been cited by other articles:
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|||||||||||
|
|||||||||||
|
|
Copyright © 2004 by the American Society of Clinical Oncology, Online ISSN: 1527-7755. Print ISSN: 0732-183X
|
TOP
ABSTRACT
INTRODUCTION
