Abstract

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 [2-6]. 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 [7-12]. 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,13-16]. 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

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 ​(Table1).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.

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 ​(Figure1),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).

Open in a separate window

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 ​(Table22 and Figure ​Figure22).

Open in a separate window

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 [23-25] 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,28-30]. 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,34-36]. 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 [37-41]. 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 [43-45]. 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.

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.

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1471-2407/8/98/prepub

Supplementary Material

Acknowledgements

References