Skip to main content
  • Research article
  • Open access
  • Published:

Prognostic Impact of Array-based Genomic Profiles in Esophageal Squamous Cell Cancer

Abstract

Background

Esophageal squamous cell carcinoma (ESCC) is a genetically complex tumor type and a major cause of cancer related mortality. Although distinct genetic alterations have been linked to ESCC development and prognosis, the genetic alterations have not gained clinical applicability. We applied array-based comparative genomic hybridization (aCGH) to obtain a whole genome copy number profile relevant for identifying deranged pathways and clinically applicable markers.

Methods

A 32 k aCGH platform was used for high resolution mapping of copy number changes in 30 stage I-IV ESCC. Potential interdependent alterations and deranged pathways were identified and copy number changes were correlated to stage, differentiation and survival.

Results

Copy number alterations affected median 19% of the genome and included recurrent gains of chromosome regions 5p, 7p, 7q, 8q, 10q, 11q, 12p, 14q, 16p, 17p, 19p, 19q, and 20q and losses of 3p, 5q, 8p, 9p and 11q. High-level amplifications were observed in 30 regions and recurrently involved 7p11 (EGFR), 11q13 (MYEOV, CCND1, FGF4, FGF3, PPFIA, FAD, TMEM16A, CTTS and SHANK2) and 11q22 (PDFG). Gain of 7p22.3 predicted nodal metastases and gains of 1p36.32 and 19p13.3 independently predicted poor survival in multivariate analysis.

Conclusion

aCGH profiling verified genetic complexity in ESCC and herein identified imbalances of multiple central tumorigenic pathways. Distinct gains correlate with clinicopathological variables and independently predict survival, suggesting clinical applicability of genomic profiling in ESCC.

Peer Review reports

Background

Esophageal squamous cell carcinoma (ESCC) is a major cause of cancer-related mortality worldwide. Despite advances in diagnostic methods and combined treatment modalities, the majority of the tumors are diagnosed at advanced stages and the overall 5-year survival rate remains 40%. ESCC develops through a multistep process from dysplasia, through carcinoma in situ to invasive carcinoma, and the acquisition of genetic alterations is tightly related to the dysplasia-carcinoma sequence [1]. The characterization of genetic alterations inherently linked to ESCC development and an in-depth understanding of the molecular mechanisms underlying carcinogenesis and growth control may therefore provide information relevant for early tumor detection, refined prognosis and development of novel targeted therapeutics.

Loss of heterozygozity (LOH) studies and conventional comparative genomic hybridization (CGH) analyses have demonstrated genetic complexity in ESCC and have identified multiple recurrent copy number alterations, namely gains of 3q, 5p, 7q, 8q, 11q, 12p, 20p and 20q [26]. Amplifications of regions harbouring oncogenes e.g. 7p12 (EGFR), 8q24 (MYC), 17q21 (FGFR) and 11q13 (e.g. CCND1, FGF4/3, and EMS1) have consistently been observed [712]. Losses, albeit at a lower frequency than gains, have recurrently involved 3p, 5q, 9p, 13q, 18q and 21q and include target genes such as FHIT, APC, RB1 and CDKN2A [3, 1316]. Moreover, some of the changes identified, e.g. gain of 8q24, 11q13, 12p, and 20q12 and loss of 3p have been associated with poor prognosis, but genetic alterations and biological characteristics have so far had a limited impact on clinical prognostication and treatment. Taken together, these findings suggest that genetic profiling can be a useful diagnostic tool in ESCC, whereas its prognostic role remains uncertain. We used array-based CGH (aCGH) for high resolution mapping in 30 ESCC and demonstrate that copy number changes detected by aCGH provide prognostic information beyond that of classical clinicopathological variables.

Methods

Tumor tissue

All 30 patients, 24 men and 6 women with a mean age of 64 (range 54–78) years, were recruited from the southern Sweden health care region and had undergone primary esophagectomy at the Department of Surgery, Lund University Hospital. None of the patients had received neoadjuvant radiotherapy or chemotherapy. Tumor tissue collected at surgery was stored at -80°C until DNA extraction. Stage (according to the International Union Against Cancer) was I in 3 cases, II in 9, III in 7, and IV in 11 cases (see Additional file 1). Presence of = 50% tumor cells in the tissue was verified by touch imprints, which were stained with Hematoxilin-Eosin and evaluated by a gastrointestinal pathologist (B.H.). All deaths were ESCC related, the median follow up was 21 (range 1–46) months for the survivors. Written informed consent was provided by all patients and the study was approved by the Lund University ethics committee.

BAC Array Platform

We used the 32 k human genome high-resolution BAC re-arrayed clone set, Version 1.0 from the BACPAC Resource Center at Children's Hospital Oakland Research Institute, Oakland (CA, US), produced at the Swegene DNA Microarray Resource Center (GEO platform repository accession GPL4723) [17], Department of Oncology, Lund University with a resolution >80 kb.

DNA isolation, labelling and hybridization

DNA was extracted, labelled and hybridized as previously described [17]. A commercial obtained DNA, derived from a pool of normal human males (Promega, Madison, WI, USA) was used as reference. Scanning was performed using an Agilent microarray scanner (Agilent Technologies, Palo Alto, CA, USA).

Image processing and data analysis

Scanned arrays were analyzed using Gene Pix Pro 4.1 (Axon Instruments, MDS Analytical Technologies, Ontario, Canada) and bad spots were "flagged" during manual inspection. The quantified data matrix was loaded into Bio Array Software Environment (BASE)[18] and filtering and normalization were performed herein [19]. Correction of background intensities of Cy3 and Cy5 were calculated using median-feature and median-local background intensities of the uploaded file. Intensity ratios were calculated from the background tumor channel (ch1) divided with the reference channel (ch2). Spots flagged as bad were filtered out from further analysis and were regarded as missing values. A signal to noise ratio (SNR) = 5 was set for both channels. Data were normalized using an implementation of a pin-based Lowess algorithm [20] in BASE excluding the X chromosome. A moving average smoothing algorithm with a 200 kbp sliding window was used and a BASE-adapted CGH-plotter software was used to identify regions of gains and losses excluding the X and the Y chromosomes [21]. A region of gain or loss was defined as two or more consecutive clones showing an absolute log2 ratio ≥ 0.2 and high-level amplifications as a log2 ratio ≥ 1.5.

Statistical Analysis

A Chi2 test was used to identify differences between copy number changes and clinicopathological characteristics. For survival analysis, the Kaplan Meier method was used to estimate relevant event variables, and the log-rank test was used to compare survival between two strata. The Cox proportional hazards model was used for univariate and multivariate survival analyses. Cox analysis was first used for all clones, thereafter two or more consecutive significant (P < 0.05) clones (interrupted only by one or two clones) were re-analysed as regions affected by gains and losses. Regions significantly (P < 0.05) linked to stage/outcome in univariate analysis were included in the multivariate analysis. All covariates were evaluated for adherence to the assumption of proportional hazards by calculating Schoenfeld residuals. Correlations between changes within a same pathway were analysed using pairwise correlation. CGH profiles were classified as gains or losses according to the dominating variable in the region corresponding to the gene locus. After Bonferroni correction correlations with P < 0.05 were considered significant. Stata 9.2 (StataCorp LP, College Station, TX 77845, USA) was used for the statistical calculations.

Immunohistochemical Staining

Serial 4-μm sections from one representative paraffin-embedded tumor block were used for immunostaining using a monoclonal antibody against human EGFR at a dilution of 1:50 (DAKO, Glostrup, Denmark). The slides were evaluated as no staining; 1+ (cytoplasmic staining or discontinuous membrane staining); 2+ (membrane staining with moderate intensity), and 3+ (intense staining with retained membranous staining). Interpretation of the staining was performed by two of the authors (AC and MN), who were blinded to the copy number changes and the clinical data.

Results

The median number of losses per tumor was 6.9% clones (range 0.1–24%) and the median number of gained clones was 9.8% (range 1.1–21%). Recurrent gains identified in at least 60% of the tumors involved chromosome regions 5p, 7p, 7q, 8q, 10q, 11q, 12p, 14q, 16p, 17p, 19p, 19q, and 20q. Likely target genes include TERT, EGFR, MYC, MYEOV, CCND1, FGF4, FGF3, CTTN and AKT1. Loss of genetic material in at least 40% of the samples affected 3p, 5q, 8p, 9p and 11q (Table 1). A homozygous deletion of 9p21.3, corresponding to the CDKN2A locus, was identified in one case and was verified with PCR using CDKN2A-specific primers (sequences available from the authors upon request, data not shown). High-level amplifications were observed in 33 regions and recurrently involved 11q22 (harbouring PDFG in 3 tumors), 11q13 (in 11 tumors, encompassing MYEOV, CCND1, FGF4, FGF3, PPFIA, FAD, TMEM16A, CTTS and SHANK2) and 7p11 (including EGFR in 4 tumors). Overexpression of EGFR was validated by immunohistochemistry, which revaled a highly positive (3+) staining in 19 tumors, 12 of which had copy number gain of 7p11. All 4 tumors with HLA of the EGFR locus showed 2+ or 3+ EGFR staining.

Table 1 Summary of the most frequent copy number gains and losses. Table listing the most frequent copy number gains and losses, sorted in decreasing order, in relation to tumor stage. Columns give the cytoband, start reporter, end reporter, size (Mbp), and candidate genes in the region.

To evaluate whether genes from a same signaling pathway were correlated, a complementary analysis of pairs of genomic loci harboring genes (CDKN2A, MDM2, RB, MDM4, CDKN1A, PIK3CA, PTEN, AKT1, TP53, MYC, CCND1, CCNE1, BCL2, CDK4, E2F3) from central tumorigenic pathways was performed. Significant correlations (P < 0.05) were identified between gain of the MDM4 locus and loss of the RB locus, gain of MDM2 and gain of TP53, gain of BCL2 and loss of CDKN1A, gain of BCL2 and loss of CDKN2A, gain of CCND1 and loss of CDKN2A. Significant negative correlations were found between gain of MDM2 and gain of PIK3CA, gain of PIK3CA and loss of E2F3, loss of PTEN and gain AKT1, and concordant gain of AKT1 and gain of CDK4. Pairwise correlation analysis between members of the EGFR pathway identified significant (P < 0.05) correlations between gains of AKT1/HER3, MAPK1/PIK3CA, and MAPK1/AKT1 and between gain of EGFR and loss of PTEN.

Despite genetic complexity (Figure 1), the genomic profiles were found to correlate with tumor stage, differentiation, and development of metastases with a lower (mean 13% versus 25%) number of changes in highly differentiated ESCC than in poorly differentiated tumors (see Additional file 2). A Chi2 test identified ~400 clones that mapped to 6 genomic regions (2p11.1-2q11.2, 2q35, 3p21.31, 4q12, 4q21.3-4q28.3, 5q12.1) that were significantly (P < 0.05) more often found in poorly differentiated tumors. The chromosomal regions most frequently affected by copy number gains differed between stage I and stage II-IV tumors. Gains of 6p25.3, 12p13.33 and 17p13.1 were present in all stage I tumors. Stage II-IV tumors showed a higher number of changes, with the most common gains affecting 12p13.33 (100%), 5p15.33 (85%), 20q13.33 (85%) and 11q13.3 (77%). Overall, copy number losses were more common in stage II tumors and the most frequent losses were shared by different tumor stages without a significant difference in frequency (see Additional file 3 and 4).

Figure 1
figure 1

A) Genome-wide frequency plot of DNA copy number gains (red) and losses (green) for all 30 ESCC tumors. B)Genome-wide copy number profile for a Stage I tumor (ESCC 33). Characteristic alterations for a Stage I tumors are gains on 6p, 12p, and 17p and losses 21q. C)Genome-wide copy number profile for a Stage IV tumor (ESCC 51). Characteristic for Stage IV tumors are gains on 12p, 19q and 20q and frequent losses on 3p, 5q and 9p.

Smaller (pT1) tumors had the lowest number of copy number changes (mean 12% of the clones affected), whereas pT2-T4 tumors showed alterations of 27%, 16%, and 18% of the clones. When gains/losses were correlated to presence of lymph node metastases, 838 clones that mapped to 21 genomic regions were identified (see Additional file 5). Copy number gain in 7p22.3 was significantly (P = 0.01) associated with lymph node metastases, and correctly predicted nodal metastases in 63% of the patients, as represented by an area under ROC curve of 0.73. Changes in 19 regions, corresponding to 1074 clones, correlated to metastasis at diagnosis (see Additional file 6) with copy number gain of 8q21.3 being significantly associated with the presence of metastasis (P = 0.03). Gain of 8q21.3 could classify 60% of the patients with distant metastasis with 40% sensitivity and 70% specificity.

Univariate analysis verified that stage and tumor size were associated with prognosis; stage HR 1.6, P = 0.04 and pT HR 1.8, P = 0.05. When copy number gains and losses were correlated to prognosis, Cox proportional hazards analysis identified 1284 clones, with a p-value <0.05, mapping to 30 regions (see Additional file 7). When gains and losses were separately considered, 7 regions remained significantly associated to prognosis in univariate analysis. When these regions were entered with stage into multivariate analysis, gain of 1p36.32 and gain of 19p13.3 independently predicted poor prognosis (Table 2 and Figure 2).

Table 2 Regions correlated to prognosis in univariate and multivariate Cox analysis. Table listing the regions significantly correlated to prognosis in univariate and multivariate analysis. Columns give the type of aberration, cytoband, hazard ratio (HR), P-value and confidence interval (CI).
Figure 2
figure 2

Kaplan-Meier survival plots of the two prognostic regions in multivariate analysis. A)Highly significant difference in survival between patients without gain and with gain of 1p36.32 (P = 0.005). B)Difference in survival between patients without gain and with gain of 19p13.3 (P = 0.01).

Discussion

Array-based genomic profiling of ESCC confirms the genetic complexity suggested by earlier studies that have applied other means of genetic profiing, e.g. cytogenetics, conventional CGH, and LOH analysis. We found copy number gains and losses affecting median 19% of the genome, identified multiple high-level amplications, and demonstrated an association between copy number alterations and stage, differentiation and prognosis, suggesting clinical applicability of genomic profiling in ESCC.

The 11q13 region is central in ESCC development and alterations herein were identified in 20/31 tumors. CCND1 is a likely target, but several other candidate genes, e.g. FGF4, FGF3, CTTN and SHANK2, showed high-level amplification. This amplicon harboured MYEOV, which has previously been associated with ESCC and described to be co-amplified with CCND1 [22]. The RB pathway is frequently targeted in ESCC carcinogenesis [2325] and its activation seems to be dependent mainly on CCND1 amplification. In our sample set no significant correlations were observed between the gains/losses observed in the members of the RB pathway (CCND1, CCNE1, E2F3 and CDKN2A). Gain of 14q32.3, which includes the AKT1 oncogene, was identified in half of the samples. The PTEN-PIK3CA-AKT signalling cascade is frequently deregulated in several types of cancers and expression of PIK3CA has been strongly associated with elevated AKT activity. An increased copy number of PIK3CA is primarily detected in tumors with retained PTEN expression [26], and indeed, none of the 11 tumors with PIK3CA gain showed loss of the PTEN locus at 10q23.3, whereas 7 tumors showed PTEN loss without change at PIK3CA locus. The pairwise analysis showed a negative (P = 0.005) correlation of both copy number gains. Expression data from array-based oligonucleotide arrays were available from 8 samples (unpublished data) and verified overexpression of PIK3CA in 7 of these tumors, which further supports PIK3CA and PTEN acting as mutually exclusive tumorigenic events [27]. Gain of 7p11.2 was identified in half of the tumors and included high-level amplifications in 4 tumors. The most likely target gene herein, EGFR, is overexpressed in a multitude of malignancies and including ESCC [7, 2830]. Immunostaining for EGFR was highly positive (3+) in 12 out of 14 tumors with copy number gain of EGFR, thus suggesting, as previously reported [28, 31, 32], that copy number gain leads to high protein expression in a significant fraction of tumors (86% in our sample set).

In our cohort, loss of PTEN was observed in 23% of the samples and was significantly correlated (P = 0.04) to EGFR gain, which may be relevant for resistance to EGFR inhibitors, since PTEN loss correlates with treatment resistency. Furthermore, gain of 17q12, harbouring ERBB2, was observed in 9 tumors and 6 of these showed concomitant gain of 7p11.2 (EGFR), which suggests that co-overexpression of ERBB2 and EGFR may apply also to ESCC [33]. High level amplification of ERBB2 correlated to overexpression (data not shown). Copy number gain of 5p15 was among the most frequent changes and the minimal region of overlap harbour some 20 identified genes, among which the telomerase regulator TERT, which has previously been shown to be overexpressed in ESCC and has been associated with prognosis in other tumor types [12, 3436]. Gains of 7p22.3, 8q22.3-qter and 20q11.21 were also frequently found and include the target genes MAD1L1 involved in TERT transcription, LRP12 and WISP1 linked to cell survival and p53-mediated apoptosis and TPX2 known to activate Aurora-A kinase [3741]. High-level amplifications affected 33 loci, among which recurrent high-level amplification peaks were detected at 7p11 (EGFR), 11q22 (cIAP1, MMP3 and PDGF), 11q13 (that harbours e.g. CCND1, FGF4, FGF3, CTTN and SHANK2), and 10q21 with unknown targets.

The most frequent recurrent copy number losses affected 3p, 5p, 8p, 9p, and 11q, which is consistent with other studies and these loci also contain several tumor suppressors linked to ESCC [2, 3, 8, 42]. Losses affecting the 9p21-p24 region, which contains CDKN2A and CDKN2B, were identified in 13/30 tumors. CDKN2A deletions have been associated with an invasive and metastatic phenotype and a homozygous CDKN2A deletion was identified in one sample [4345]. Frequent losses were also observed at 3p26-p14 which harbours THRB, RARB, TOP2B and FHIT. Pairwise correlations between frequently observed gains and losses identified 5 regions that were significantly more often affected by concurrent aberrations. Four of these were located on the same chromosome, whereas loss of 3p24 and 5q12 occurred at an increased incidence. These regions contain targets such as TOP2B, RARB, and TGFBR2 on 3p and PIK3R1 and RAD17 on 5q, and the association identified may indicate cross-talk between genes in these regions.

Conclusion

The accumulation of genetic changes is central in ESCC development and progression. Our study is the first to apply high-resolution aCGH to clinical prognostication in ESCC. Studies that have applied traditional CGH have suggested a prognostic independent role for genes located on 8q, 11q, 12p, 14q, and 20q [2, 16, 46], and a recent study using a 4k BAC array platform (with 1 Mb resolution) identified 4 clones, mapping to 3q29, 4q21.21, 8q24 and 8q24.3, linked to survival [47]. The high-resolution data here presented demonstrates extensive genetic complexity already in early stage tumors, supports the involvement of several key genes in ESCC, links gain of 7p22.3 to presence of nodal metastases and demonstrates that gains of 1p36.32 and 19p13.3 provide independent prognostic information (HR = 19.6 and 7.0 respectively). The different prognostic regions identified may be related to the inherent genetic complexity of ESCC, to differences in materials (e.g. study populations from different geographic areas with disparities in dietary and environmental exposures) and use of different genetic profiling technologies [48, 49]. Nevertheless, these results hold promise for the application as genetic classifiers and refined prognostic markers. Moreover, the recognition of recurrent rearrangements in central signaling pathways provides a basis for the development of selected and individualized targeted therapeutics in ESCC.

References

  1. Metzger R, Schneider PM, Warnecke-Eberz U, Brabender J, Holscher AH: Molecular biology of esophageal cancer. Onkologie. 2004, 27 (2): 200-206. 10.1159/000076913.

    CAS  PubMed  Google Scholar 

  2. Yen CC, Chen YJ, Chen JT, Hsia JY, Chen PM, Liu JH, Fan FS, Chiou TJ, Wang WS, Lin CH: Comparative genomic hybridization of esophageal squamous cell carcinoma: correlations between chromosomal aberrations and disease progression/prognosis. Cancer. 2001, 92 (11): 2769-2777. 10.1002/1097-0142(20011201)92:11<2769::AID-CNCR10118>3.0.CO;2-M.

    Article  CAS  PubMed  Google Scholar 

  3. Pack SD, Karkera JD, Zhuang Z, Pak ED, Balan KV, Hwu P, Park WS, Pham T, Ault DO, Glaser M, et al: Molecular cytogenetic fingerprinting of esophageal squamous cell carcinoma by comparative genomic hybridization reveals a consistent pattern of chromosomal alterations. Genes, chromosomes & cancer. 1999, 25 (2): 160-168. 10.1002/(SICI)1098-2264(199906)25:2<160::AID-GCC12>3.0.CO;2-U.

    Article  CAS  Google Scholar 

  4. Du Plessis L, Dietzsch E, Van Gele M, Van Roy N, Van Helden P, Parker MI, Mugwanya DK, De Groot M, Marx MP, Kotze MJ, et al: Mapping of novel regions of DNA gain and loss by comparative genomic hybridization in esophageal carcinoma in the Black and Colored populations of South Africa. Cancer research. 1999, 59 (8): 1877-1883.

    CAS  PubMed  Google Scholar 

  5. Tada K, Oka M, Hayashi H, Tangoku A, Oga A, Sasaki K: Cytogenetic analysis of esophageal squamous cell carcinoma cell lines by comparative genomic hybridization: relationship of cytogenetic aberrations to in vitro cell growth. Cancer Genet Cytogenet. 2000, 117 (2): 108-112. 10.1016/S0165-4608(99)00160-0.

    Article  CAS  PubMed  Google Scholar 

  6. Shinomiya T, Mori T, Ariyama Y, Sakabe T, Fukuda Y, Murakami Y, Nakamura Y, Inazawa J: Comparative genomic hybridization of squamous cell carcinoma of the esophagus: the possible involvement of the DPI gene in the 13q34 amplicon. Genes, chromosomes & cancer. 1999, 24 (4): 337-344. 10.1002/(SICI)1098-2264(199904)24:4<337::AID-GCC7>3.0.CO;2-O.

    Article  CAS  Google Scholar 

  7. Ishizuka T, Tanabe C, Sakamoto H, Aoyagi K, Maekawa M, Matsukura N, Tokunaga A, Tajiri T, Yoshida T, Terada M, et al: Gene amplification profiling of esophageal squamous cell carcinomas by DNA array CGH. Biochem Biophys Res Commun. 2002, 296 (1): 152-155. 10.1016/S0006-291X(02)00836-7.

    Article  CAS  PubMed  Google Scholar 

  8. Arai H, Ueno T, Tangoku A, Yoshino S, Abe T, Kawauchi S, Oga A, Furuya T, Oka M, Sasaki K: Detection of amplified oncogenes by genome DNA microarrays in human primary esophageal squamous cell carcinoma: comparison with conventional comparative genomic hybridization analysis. Cancer Genet Cytogenet. 2003, 146 (1): 16-21. 10.1016/S0165-4608(03)00106-7.

    Article  CAS  PubMed  Google Scholar 

  9. Jiang W, Kahn SM, Tomita N, Zhang YJ, Lu SH, Weinstein IB: Amplification and expression of the human cyclin D gene in esophageal cancer. Cancer research. 1992, 52 (10): 2980-2983.

    CAS  PubMed  Google Scholar 

  10. Shinozaki H, Ozawa S, Ando N, Tsuruta H, Terada M, Ueda M, Kitajima M: Cyclin D1 amplification as a new predictive classification for squamous cell carcinoma of the esophagus, adding gene information. Clin Cancer Res. 1996, 2 (7): 1155-1161.

    CAS  PubMed  Google Scholar 

  11. Hollstein MC, Smits AM, Galiana C, Yamasaki H, Bos JL, Mandard A, Partensky C, Montesano R: Amplification of epidermal growth factor receptor gene but no evidence of ras mutations in primary human esophageal cancers. Cancer research. 1988, 48 (18): 5119-5123.

    CAS  PubMed  Google Scholar 

  12. Lu SH, Hsieh LL, Luo FC, Weinstein IB: Amplification of the EGF receptor and c-myc genes in human esophageal cancers. International journal of cancer. 1988, 42 (4): 502-505. 10.1002/ijc.2910420406.

    Article  CAS  PubMed  Google Scholar 

  13. Shiomi H, Sugihara H, Kamitani S, Tokugawa T, Tsubosa Y, Okada K, Tamura H, Tani T, Kodama M, Hattori T: Cytogenetic heterogeneity and progression of esophageal squamous cell carcinoma. Cancer Genet Cytogenet. 2003, 147 (1): 50-61. 10.1016/S0165-4608(03)00159-6.

    Article  CAS  PubMed  Google Scholar 

  14. Hu N, Roth MJ, Polymeropolous M, Tang ZZ, Emmert-Buck MR, Wang QH, Goldstein AM, Feng SS, Dawsey SM, Ding T, et al: Identification of novel regions of allelic loss from a genomewide scan of esophageal squamous-cell carcinoma in a high-risk Chinese population. Genes, chromosomes & cancer. 2000, 27 (3): 217-228. 10.1002/(SICI)1098-2264(200003)27:3<217::AID-GCC1>3.0.CO;2-A.

    Article  CAS  Google Scholar 

  15. Huang Y, Boynton RF, Blount PL, Silverstein RJ, Yin J, Tong Y, McDaniel TK, Newkirk C, Resau JH, Sridhara R, et al: Loss of heterozygosity involves multiple tumor suppressor genes in human esophageal cancers. Cancer research. 1992, 52 (23): 6525-6530.

    CAS  PubMed  Google Scholar 

  16. Kwong D, Lam A, Guan X, Law S, Tai A, Wong J, Sham J: Chromosomal aberrations in esophageal squamous cell carcinoma among Chinese: gain of 12p predicts poor prognosis after surgery. Hum Pathol. 2004, 35 (3): 309-316. 10.1016/j.humpath.2003.10.020.

    Article  CAS  PubMed  Google Scholar 

  17. Jonsson G, Staaf J, Olsson E, Heidenblad M, Vallon-Christersson J, Osoegawa K, de Jong P, Oredsson S, Ringner M, Hoglund M, et al: High-resolution genomic profiles of breast cancer cell lines assessed by tiling BAC array comparative genomic hybridization. Genes, chromosomes & cancer. 2007, 46 (6): 543-558. 10.1002/gcc.20438.

    Article  Google Scholar 

  18. The BioArray Software Environment (BASE) website. [http://base2.thep.lu.se/onk/]

  19. Saal LH, Troein C, Vallon-Christersson J, Gruvberger S, Borg A, Peterson C: BioArray Software Environment (BASE): a platform for comprehensive management and analysis of microarray data. Genome Biol. 2002, 3 (8): SOFTWARE0003-10.1186/gb-2002-3-8-software0003.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002, 30 (4): e15-10.1093/nar/30.4.e15.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Autio R, Hautaniemi S, Kauraniemi P, Yli-Harja O, Astola J, Wolf M, Kallioniemi A: CGH-Plotter: MATLAB toolbox for CGH-data analysis. Bioinformatics. 2003, 19 (13): 1714-1715. 10.1093/bioinformatics/btg230.

    Article  CAS  PubMed  Google Scholar 

  22. Janssen JW, Imoto I, Inoue J, Shimada Y, Ueda M, Imamura M, Bartram CR, Inazawa J: MYEOV, a gene at 11q13, is coamplified with CCND1, but epigenetically inactivated in a subset of esophageal squamous cell carcinomas. J Hum Genet. 2002, 47 (9): 460-464. 10.1007/s100380200065.

    Article  CAS  PubMed  Google Scholar 

  23. Kawakubo H, Ozawa S, Ando N, Kitagawa Y, Mukai M, Ueda M, Kitajima M: Alterations of p53, cyclin D1 and pRB expression in the carcinogenesis of esophageal squamous cell carcinoma. Oncology reports. 2005, 14 (6): 1453-1459.

    CAS  PubMed  Google Scholar 

  24. Nakajima Y, Miyake S, Tanaka K, Ogiya K, Toukairin Y, Kawada K, Nishikage T, Nagai K, Kawano T: The expressions of p21 and pRB may be good indicators for the sensitivity of esophageal squamous cell cancers to CPT-11: Cell proliferation activity correlates with the effect of CPT-11. Cancer science. 2004, 95 (5): 464-468. 10.1111/j.1349-7006.2004.tb03233.x.

    Article  CAS  PubMed  Google Scholar 

  25. Mathew R, Arora S, Khanna R, Mathur M, Shukla NK, Ralhan R: Alterations in p53 and pRb pathways and their prognostic significance in oesophageal cancer. Eur J Cancer. 2002, 38 (6): 832-841. 10.1016/S0959-8049(02)00007-2.

    Article  CAS  PubMed  Google Scholar 

  26. Altomare DA, Testa JR: Perturbations of the AKT signaling pathway in human cancer. Oncogene. 2005, 24 (50): 7455-7464. 10.1038/sj.onc.1209085.

    Article  CAS  PubMed  Google Scholar 

  27. Pedrero JM, Carracedo DG, Pinto CM, Zapatero AH, Rodrigo JP, Nieto CS, Gonzalez MV: Frequent genetic and biochemical alterations of the PI 3-K/AKT/PTEN pathway in head and neck squamous cell carcinoma. International journal of cancer. 2005, 114 (2): 242-248. 10.1002/ijc.20711.

    Article  CAS  PubMed  Google Scholar 

  28. Hanawa M, Suzuki S, Dobashi Y, Yamane T, Kono K, Enomoto N, Ooi A: EGFR protein overexpression and gene amplification in squamous cell carcinomas of the esophagus. International journal of cancer. 2006, 118 (5): 1173-1180. 10.1002/ijc.21454.

    Article  CAS  PubMed  Google Scholar 

  29. Sunpaweravong P, Sunpaweravong S, Puttawibul P, Mitarnun W, Zeng C, Baron AE, Franklin W, Said S, Varella-Garcia M: Epidermal growth factor receptor and cyclin D1 are independently amplified and overexpressed in esophageal squamous cell carcinoma. J Cancer Res Clin Oncol. 2005, 131 (2): 111-119. 10.1007/s00432-004-0610-7.

    Article  CAS  PubMed  Google Scholar 

  30. Gibault L, Metges JP, Conan-Charlet V, Lozac'h P, Robaszkiewicz M, Bessaguet C, Lagarde N, Volant A: Diffuse EGFR staining is associated with reduced overall survival in locally advanced oesophageal squamous cell cancer. Br J Cancer. 2005, 93 (1): 107-115. 10.1038/sj.bjc.6602625.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Janmaat ML, Gallegos-Ruiz MI, Rodriguez JA, Meijer GA, Vervenne WL, Richel DJ, Van Groeningen C, Giaccone G: Predictive factors for outcome in a phase II study of gefitinib in second-line treatment of advanced esophageal cancer patients. J Clin Oncol. 2006, 24 (10): 1612-1619. 10.1200/JCO.2005.03.4900.

    Article  CAS  PubMed  Google Scholar 

  32. Personeni N: Epidermal growth factor receptor gene copy number in esophageal cancer and outcome prediction to gefitinib: does intratumoral heterogeneity matter?. J Clin Oncol. 2006, 24 (34): 5465-10.1200/JCO.2006.07.5622. author reply 5466–5467

    Article  PubMed  Google Scholar 

  33. DiGiovanna MP, Stern DF, Edgerton SM, Whalen SG, Moore D, Thor AD: Relationship of epidermal growth factor receptor expression to ErbB-2 signaling activity and prognosis in breast cancer patients. J Clin Oncol. 2005, 23 (6): 1152-1160. 10.1200/JCO.2005.09.055.

    Article  CAS  PubMed  Google Scholar 

  34. Hsu CP, Lee LW, Shai SE, Chen CY: Clinical significance of telomerase and its associate genes expression in the maintenance of telomere length in squamous cell carcinoma of the esophagus. World J Gastroenterol. 2005, 11 (44): 6941-6947.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Barclay JY, Morris A, Nwokolo CU: Telomerase, hTERT and splice variants in Barrett's oesophagus and oesophageal adenocarcinoma. European journal of gastroenterology & hepatology. 2005, 17 (2): 221-227. 10.1097/00042737-200502000-00014.

    Article  CAS  Google Scholar 

  36. Zhang H, Jin Y, Chen X, Jin C, Law S, Tsao SW, Kwong YL: Cytogenetic aberrations in immortalization of esophageal epithelial cells. Cancer Genet Cytogenet. 2006, 165 (1): 25-35. 10.1016/j.cancergencyto.2005.07.016.

    Article  CAS  PubMed  Google Scholar 

  37. Ohta Y, Hamada Y, Saitoh N, Katsuoka K: Effect of the transcriptional repressor Mad1 on proliferation of human melanoma cells. Exp Dermatol. 2002, 11 (5): 439-447. 10.1034/j.1600-0625.2002.110507.x.

    Article  CAS  PubMed  Google Scholar 

  38. Garnis C, Coe BP, Zhang L, Rosin MP, Lam WL: Overexpression of LRP12, a gene contained within an 8q22 amplicon identified by high-resolution array CGH analysis of oral squamous cell carcinomas. Oncogene. 2004, 23 (14): 2582-2586. 10.1038/sj.onc.1207367.

    Article  CAS  PubMed  Google Scholar 

  39. Xu L, Corcoran RB, Welsh JW, Pennica D, Levine AJ: WISP-1 is a Wnt-1- and beta-catenin-responsive oncogene. Genes Dev. 2000, 14 (5): 585-595.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Wang XX, Liu R, Jin SQ, Fan FY, Zhan QM: Overexpression of Aurora-A kinase promotes tumor cell proliferation and inhibits apoptosis in esophageal squamous cell carcinoma cell line. Cell Res. 2006, 16 (4): 356-366. 10.1038/sj.cr.7310046.

    Article  CAS  PubMed  Google Scholar 

  41. Tong T, Zhong Y, Kong J, Dong L, Song Y, Fu M, Liu Z, Wang M, Guo L, Lu S, et al: Overexpression of Aurora-A contributes to malignant development of human esophageal squamous cell carcinoma. Clin Cancer Res. 2004, 10 (21): 7304-7310. 10.1158/1078-0432.CCR-04-0806.

    Article  CAS  PubMed  Google Scholar 

  42. Noguchi T, Takeno S, Kimura Y, Uchida Y, Daa T, Yokoyama S, Gabbert HE, Mueller W: FHIT expression and hypermethylation in esophageal squamous cell carcinoma. Int J Mol Med. 2003, 11 (4): 441-447.

    CAS  PubMed  Google Scholar 

  43. Takeuchi H, Ozawa S, Ando N, Kitagawa Y, Mukai M, Ueda M, Kitajima M: Further evidence that altered p16/CDKN2 gene expression is associated with lymph node metastasis in squamous cell carcinoma of the esophagus. Oncology reports. 2001, 8 (3): 627-632.

    CAS  PubMed  Google Scholar 

  44. Tokugawa T, Sugihara H, Tani T, Hattori T: Modes of silencing of p16 in development of esophageal squamous cell carcinoma. Cancer research. 2002, 62 (17): 4938-4944.

    CAS  PubMed  Google Scholar 

  45. Takeuchi H, Ozawa S, Ando N, Shih CH, Koyanagi K, Ueda M, Kitajima M: Altered p16/MTS1/CDKN2 and cyclin D1/PRAD-1 gene expression is associated with the prognosis of squamous cell carcinoma of the esophagus. Clin Cancer Res. 1997, 3 (12 Pt 1): 2229-2236.

    CAS  PubMed  Google Scholar 

  46. Ueno T, Tangoku A, Yoshino S, Abe T, Hayashi H, Toshimitsu H, Hashimoto K, Satoh T, Oga A, Furuya T, et al: Prediction of nodal metastasis by comparative genomic hybridization in biopsy specimens from patients with superficial esophageal squamous cell carcinoma. Clin Cancer Res. 2003, 9 (14): 5137-5141.

    CAS  PubMed  Google Scholar 

  47. Hirasaki S, Noguchi T, Mimori K, Onuki J, Morita K, Inoue H, Sugihara K, Mori M, Hirano T: BAC clones related to prognosis in patients with esophageal squamous carcinoma: an array comparative genomic hybridization study. The oncologist. 2007, 12 (4): 406-417. 10.1634/theoncologist.12-4-406.

    Article  CAS  PubMed  Google Scholar 

  48. Wang Z, Tang L, Sun G, Tang Y, Xie Y, Wang S, Hu X, Gao W, Cox SB, Wang JS: Etiological study of esophageal squamous cell carcinoma in an endemic region: a population-based case control study in Huaian, China. BMC cancer. 2006, 6: 287-10.1186/1471-2407-6-287.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Wang JM, Xu B, Rao JY, Shen HB, Xue HC, Jiang QW: Diet habits, alcohol drinking, tobacco smoking, green tea drinking, and the risk of esophageal squamous cell carcinoma in the Chinese population. European journal of gastroenterology & hepatology. 2007, 19 (2): 171-176. 10.1097/MEG.0b013e32800ff77a.

    Article  CAS  Google Scholar 

Pre-publication history

Download references

Acknowledgements

The study was financially supported by grants from the Knut and Alice Wallenberg Foundation via the Swegene program, the Swedish Cancer Society, the Swedish Research Council, the Nilsson Cancer Research Fund, the Kamprad Research Fund, the Region Skåne Research Funds and the American Cancer Society (GJ). AC was supported by a grant from the Portuguese Institute of Oncology (Lisbon) and from the Portuguese League Against Cancer (South Region Centre).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ana Carneiro.

Additional information

Competing interests

The author(s) declare that they have no competing interests.

Authors' contributions

AC designed the study, performed experiments, analyzed the data and drafted the manuscript. AI contributed to the design of the study, performed experiments and helped to draft the manuscript. AK carried out experiments and analyzed data. JJ provided cancer samples and clinical data. GJ helped the data analysis. POB performed and supervised statistical analysis. DF provided cancer samples. BH performed pathological analysis of tumor samples. MN conceived, designed and coordinated the study, and helped to draft the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12885_2007_1036_MOESM1_ESM.xls

Additional file 1: Clinicopathologic variables of tumors and survival of patients of the sample set. Table listing the clinicopathologic variables of tumors and survival of patients of the sample set analyzed in this study. (XLS 20 KB)

12885_2007_1036_MOESM2_ESM.tiff

Additional file 2: Number of gains and losses by differentiation. Boxplot of the number of gained and lost clones according to tumor differentiation. (TIFF 15 MB)

12885_2007_1036_MOESM3_ESM.tiff

Additional file 3: Number of gains and losses by tumor stage. Boxplot of the number of gained and lost clones according to tumor stage. (TIFF 16 MB)

12885_2007_1036_MOESM4_ESM.xls

Additional file 4: Copy number alterations, listed in decrescent prevalence order, in relation to tumor stage. Table listing the most frequent copy number alterations, sorted in decreasing order, in relation to tumor stage. Columns give the cytoband, start reporter, end reporter, size (Mbp), frequency of the change and candidate genes in the region. (XLS 29 KB)

12885_2007_1036_MOESM5_ESM.xls

Additional file 5: Chi2 test significant regions associated with positive lymph nodes (P < 0.05). Table listing the regions of aberrations that showed a significant statistical (P < 0.05) association with positive lymph nodes. Given are cytoband, start reporter, end reporter, size (Mbp) and the frequency across patients with positive lymph nodes and patients with negative lymph nodes. P-values were computed using Chi2 test. (XLS 24 KB)

12885_2007_1036_MOESM6_ESM.xls

Additional file 6: Chi2 test significant regions associated with metastasis at diagnosis (P < 0.05). Table listing the regions of aberrations that showed a significant statistical (P < 0.05) association with metastasis at diagnosis. Given are cytoband, start reporter, end reporter, size (Mbp) and the frequency across patients with metastasis at diagnosis and patients without metastasis at diagnosis. P-values were computed using Chi2 test. (XLS 21 KB)

12885_2007_1036_MOESM7_ESM.xls

Additional file 7: Regions correlated to prognosis identified by Cox analysis. Table listing the regions of aberrations significantly correlated (P < 0.05) to prognosis in univariate Cox analysis. Cytoband, start reporter, end reporter and size (Mbp) of the aberrant region are given. (XLS 22 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Carneiro, A., Isinger, A., Karlsson, A. et al. Prognostic Impact of Array-based Genomic Profiles in Esophageal Squamous Cell Cancer. BMC Cancer 8, 98 (2008). https://doi.org/10.1186/1471-2407-8-98

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2407-8-98

Keywords