- Research article
- Open Access
- Open Peer Review
Integrated mutation, copy number and expression profiling in resectable non-small cell lung cancer
BMC Cancervolume 11, Article number: 93 (2011)
The aim of this study was to identify critical genes involved in non-small cell lung cancer (NSCLC) pathogenesis that may lead to a more complete understanding of this disease and identify novel molecular targets for use in the development of more effective therapies.
Both transcriptional and genomic profiling were performed on 69 resected NSCLC specimens and results correlated with mutational analyses and clinical data to identify genetic alterations associated with groups of interest.
Combined analyses identified specific patterns of genetic alteration associated with adenocarcinoma vs. squamous differentiation; KRAS mutation; TP53 mutation, metastatic potential and disease recurrence and survival. Amplification of 3q was associated with mutations in TP53 in adenocarcinoma. A prognostic signature for disease recurrence, reflecting KRAS pathway activation, was validated in an independent test set.
These results may provide the first steps in identifying new predictive biomarkers and targets for novel therapies, thus improving outcomes for patients with this deadly disease.
Non-small cell lung cancer (NSCLC) is the commonest cause of cancer death in Western communities. Current treatments offer the potential of cure only to the small number of patients who present with early stage NSCLC, whilst outcomes for those with advanced disease remain poor. Recent advances including adjuvant chemotherapy and targeted biological therapies have lead to modest improvements in survival for small subgroups of patients. Clearly new treatment approaches are required to substantially improve outcome. As has been the case for other tumour types, molecular profiling techniques have the potential to provide benefit through improved understanding of disease pathogenesis, identification of subgroups in whom current therapies are most likely to be effective and in the development of novel therapies.
Genetic heterogeneity is a feature of NSCLC, with varying combinations of multiple molecular alterations contributing to tumour development . A key challenge for high-throughput molecular profiling techniques is to distinguish between genes whose expression is altered directly by heritable changes in gene function and those where changes are an inevitable down-stream consequence of primary changes to genes directly involved in disease pathogenesis. Correlation of transcriptional and genomic data allows more focussed analysis of the large number of genetic alterations identified by molecular profiling. This study incorporates the results of both transcriptional and genomic profiling for clinically relevant subgroups of NSCLC to identify genes of potential predictive or pathogenic importance in this deadly disease.
After obtaining institutional ethics approval, patients with stage I-IIIA NSCLC seen by the St Vincent's Hospital Combined Lung Service between February 2004 and July 2006 and planned for curative resection were invited to participate. Exclusion criteria included age <18 years, administration of neoadjuvant chemotherapy, and inability to provide informed consent. Integrated demographic, radiological, pathological and outcome data was collected for all consenting patients.
In addition, a small number of samples collected earlier and stored in the Peter MacCallum (PeterMac) tissue bank were utilised after approval by the PeterMac Tissue Management Committee.
Samples of tumour (≥1 cm 3) were selected from fresh specimens, and then stored whole at -180°C. Only those specimens containing >75% tumour cells and <25% necrosis were used in molecular studies. Both RNA and DNA were isolated from each sample for analysis using established protocols (see additional files 1 and 2). Transcriptional profiling was performed using 10,500 element cDNA microarrays (PeterMac, Melbourne, Australia) . Genomic profiling using 2400 element bacterial artificial chromosome (BAC) arrays was completed at the University of San Francisco, California, USA . Detailed description of transcriptional and genomic profiling is included in additional files 3 and 4.
All samples were tested for TP53 mutations and all adenocarcinoma (AC) and large cell carcinoma (LCC) samples were screened for KRAS mutations using high resolution melting analysis [4, 5] with or without DNA sequencing.
The effect of histology, presence or absence of KRAS or TP53 mutation, tumour size and metastasis status, recurrence within 12 months of surgery and survival on gene expression was explored. After removing control genes, analysis was conducted in CRAN, R Bioconductor using the LIMMA [6, 7] package to generate detailed lists of gene expression differences with significance p values between each subgroup of interest. To account for multiple testing, p values of <0.005 were considered statistically significant. Gene lists were then interrogated using publicly available programs (Intelligent Systems and Bioinformatics Laboratory, Wayne State University, Detroit, MI, USA, http://vortex.cs.wayne.edu/projects.htm) to identify gene ontology and molecular pathway patterns.
Changes in the normalized and smoothed genomic data  were assigned stepwise copy change levels from -2 to 3 (-2 = homozygous loss, -1 = heterozygous loss, 1 = single copy gain, 2 = gain of two copies, 3 = high-level gain, 0 = normal copy number). Using these standardised copy number values, it was possible to make comparisons of the frequency of each level of copy change at each BAC location for specific groups of interest.
To compare our data with that of other groups of patients with early stage NSCLC, we performed comparisons with publicly available external data sets obtained via the NCBI Gene Expression Omnibus (GEO) website. Different platforms were reconciled using HUGO approved gene symbols, and extracted gene expression data were log2 transformed, centred and scaled across samples in order to emphasise relative expressions as opposed to absolute values.
Integration of transcriptional and genomic profiles
Integration of transcriptional and genomic datasets was performed by investigating levels of expression for genes located in regions of copy number variation between groups of interest. Genes whose differences in expression varied in the same direction as differences in copy number between two groups (e.g. relative over-expression of genes in a region of increased copy number) were viewed as genes of interest.
Molecular and clinical data was available for 69 patients who underwent surgery for NSCLC between May 1999 and July 2006. Demographic and pathologic data are included in table 1. Median follow-up for surviving patients exceeds 35 months. After a median follow-up of 36 months (1 - 80) for all patients, 28/68 patients developed recurrent disease, and 23/68 patients died of NSCLC (one patient with disseminated disease at diagnosis excluded from analysis). Comparable to other series of early stage NSCLC , five year overall survival rates approximated 55%.
All samples demonstrated significant chromosomal instability with an average of 43.6 chromosomal breakpoints per sample (defined by a change in the stepwise copy number along a chromosome), with over 100 and 150 regions of high-level gain (+3) and loss (-3) respectively. There was also a very high rate of low-level genomic alteration (both gains and losses). On average, over 10 whole arm losses or duplications were seen per sample, with a rate of isochromosome formation of 1.6 per genome (duplication of one arm with loss of the opposing arm of the same chromosome). These results are consistent with the highly disordered nature of lung cancer genomes. Comparisons between clinical subgroups of interest revealed remarkably similar degrees of aneuploidy and chromosomal disorder in almost all groups (Table 2). Specifically, neither prognosis, degree of histologic differentiation, K-Ras or TP53 status was associated with evidence of greater aneuploidy.
Charts documenting the averaged standardised copy numbers at each BAC location were generated to enable visual comparisons of genomic profiles between AC (23) and SCC (17) (figure 1). Regions of shared change (amplification of 5p15.33-p13.13, loss of 3p26.3-p13) and several regions of difference (SCC: gains at 3q12.1-q28 and 12p13.33-p12.1, losses of chromosome 4 (4pter-4qter); AC: gains of 6p25.1-p21.31, losses of 6q13-27, and a greater magnitude of change in several regions of shared loss (9p, 13q, 17p, 18q, 19p, 19q and 22q) or gain (7p, 7q, 8q)) were identified. These findings are consistent with previous studies [9–16].
Associations with mutation status in TP53, KRAS and EGFR pathways
Clinical data demonstrated a trend to greater rate of TP53 mutation in SCC than AC (TP53 mutation in 9/12 (75%) SCC and 9/18 (50%) AC, p = 0.083). Amplification of 3q was also more frequent in SCC than AC samples (p = 0.004). When analysing all samples, no relationship was found between 3q amplification and TP53 mutation (p = 0.99). However, when analysing SCC and AC separately, a statistically significant relationship between TP53 mutation and 3q amplification was detected in AC samples, with amplification of 3q being significantly more common in TP53 mutant cancers (1/8 and 4/10 samples with 3q amplification in TP53 wild-type (wt) and mutant AC groups respectively, p = 0.027). Both TP53 mutant and wt samples more frequently demonstrated copy number loss at the TP53 locus (17p13) than gain (11/18 and 9/12 samples in TP53 mutant and wt groups respectively).
Adenocarcinomas were screened for EGFR and KRAS mutations. Small numbers of EGFR mutant tumours limited detailed analysis. Good quality genomic profiles were available for only 5 tumours with EGFR mutation and no significant differences were seen between the profiles of EGFR mutant and non-mutant tumours.
Genomic profiles were available for 8 KRAS mutant tumours, 31 wt and 16 tumours with unknown KRAS mutation status (untested SCC samples). In averaged copy number charts, KRAS mutant tumours showed a predilection for losses of 1p36.32-p13.2, 6q11.1-q27, 11p13-11q13.2 and 11q21-12p13.1 and gains of 1q21.1-q43 (Figure 2).
Associations with metastasis, tumour recurrence and NSCLC-specific survival
To investigate the notion of inherent metastatic potential, molecular profiles of large (>4 cm) non-metastatic tumours and small (<2 cm) metastatic (nodal or distant) tumours were compared. Genomic profiles of 3 'metastatic' and 8 'non-metastatic' tumours revealed some differences in the magnitude of copy number changes, without regions of clear difference between the two groups. There was no correlation between genomic changes and tumour recurrence or survival. There were differences in the magnitude of gene copy number changes at 7p, 8q, 9p, 15q and 17p in recurrent compared to non-recurrent tumours. Contained within these regions are the MYC oncogene (8q), as well as TP53 (17p), and the CDKN2A locus (containing p14(ARF) and p16 tumour suppressor genes (9p) (TSG's)).
Ranking the genes by moderated t-statistics and selecting a p value cut-off of <0.005, 310 genes with differential expression between 16 SCC and 25 AC samples were identified, representing the biological processes of cell adhesion, epidermis development, keratinisation and keratinocyte differentiation. A significant proportion of these genes had roles in antigen processing and presentation, and the phosphatidylinositol signalling pathway. Thirty of 310 genes in the differentiating gene list were located on chromosome 3 (p = 0.0098), implicating genomic changes at this locus in determining NSCLC phenotype. This is consistent with the genomic data, which indicates gain of 3q is associated with SCC histology.
Associations with mutation status in TP53, K-Ras and EGFR pathways
Expression levels of 67 genes differed significantly between TP53 mutant (17) and wt (21) tumours. Many of the biological functions represented by these genes were also strongly represented by the genes differentially expressed between AC and SCC. In addition, 20/67 discriminating genes were also included in the gene list differentiating SCC from AC. Hierarchical clustering based on the expression of these 67 differentially expressed genes not only segregated TP53 mutant from wild-type tumours, but also resulted in clustering of SCC samples with the TP53 mutant tumours. Our results suggest that the gene expression signature observed for TP53 mutant tumours may be at least in part related to SCC histology rather than TP53 biology.
Transcriptional profiles of AC and LCC tumours with (8) and without (31) KRAS mutation were compared. Biological processes represented by 108 differentiating genes included cell growth, second-messenger mediated signalling, chromosome organisation and biogenesis, and gene regulation (mediated via histones and their effects on biosynthesis and nucleosome assembly) (table 3). These findings are consistent with other published studies linking KRAS mutation to increased translation of cancer related proteins, and chromosome instability [17, 18]. The low frequency of EGFR mutant cancers precluded statistically meaningful analysis of transcriptional data according to EGFR genotype.
Associations with metastasis, tumour recurrence and NSCLC-specific survival
Transcriptional profiles identified 39 genes that differentiated between 19 'metastatic' and 35 'non-metastatic' tumours, with molecular pathways involved in protein translation most strongly represented (MRPL33, RPL12, RPL27A, RPS5, RPS9). Comparison of expression profiles of 14 tumours recurring within 12 months of surgery to remaining samples identified 60 genes with differential expression between the two groups, with a common theme of RAS activation represented in ontological and single gene analyses. Included in the differentiating gene list were MAPK1, DUSP11 and DUSP13, PTPN11, and PIK3CB, all having roles in signal transduction and the MAPK pathway. The phosphatidylinositol signalling pathway was also significantly over-represented in ontology analysis.
Expression levels of only 38 genes differed significantly between deceased and surviving patients. 18 of these genes were shared with gene lists of recurrent vs. non-recurrent tumours. Few biological processes were represented by more than one gene, and clear patterns of gene ontology were not apparent.
Correlation with External Data Sets
Comparison of our differential gene list for recurrence with the discriminating gene list for survival in GSE11117 (transcriptional and survival data for 41 NSCLC samples, using Novachip Human 34.5k microarray interrogating ~34,500 transcripts for each sample; http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE11117) identified 40 matched transcripts (additional file 5). Log-transformed expression values of the 40 transcripts were used to classify the samples from GSE11117 into two subgroups using a correlation, average-linkage hierarchical clustering R package. Kaplan-Meier curves for these external samples using the 40 transcripts matched to our recurrence gene list demonstrated statistically significant survival prediction (figure 3, p < 0.0153)), with 21 and 20 samples in each group.
Integration of Genomic and Transcriptional Profiles
To determine whether integration of genomic and expression data added to the predictive value of these datasets for histologic classification, we identified 34 genes whose copy number and expression varied between AC and SCC, 24 of which demonstrated concordant differences in copy number and expression (table 4). Notably, 17 of 24 differentiating genes were located on chromosome 3q.
Associations with mutation status
Examination of transcriptional data from 1p, 1q, 6q, 11p, 11q and 12p (regions of genomic difference) identified 25 genes with concordant differences in copy number and expression between K-Ras mutant and wt tumours (table 5). A number of genes demonstrated reduced copy number and expression in KRAS mutant tumours, including putative tumour suppressor genes (FOXO3, EXTL2, PPP2R1B), negative regulators of the receptor tyrosine kinase oncogenic pathways (PTPRK, DGKZ, NCAM1), and negative regulators of Ras (EPHB2). Several over-expressed genes located in regions of amplification play roles in constitutive KRAS activation (PTGS2/COX2), enhanced transactivation of the EGFR (RGS2), enhanced invasive potential (ECM1), and MAPK/ERK activation (PTGS2/COX2).
Associations with metastasis, tumour recurrence and NSCLC-specific survival
Investigation of genomic and transcriptional data identified only 2 genes (ARFGEF1 and PENK) whose copy number and expression differentiated 'metastatic' from 'non-metastatic' tumours, neither of which have been previously implicated in malignancy. Similarly, correlation of transcriptional and genomic data identified only 3 genes with concordant differences in recurrence and survival comparisons - SMARCA2, MINK and RECK.
The clinical, demographic and pathologic characteristics of this NSCLC cohort are consistent with the published literature. The transcriptional and genomic profiles identified in this study should therefore be generalisable to other patients with early-stage NSCLC. The tumour samples analysed demonstrated substantial genomic instability, with comparisons between subgroups failing to demonstrate any significant difference. Previous studies of copy number changes in NSCLC have found no association between age, gender, histology, stage or tumour grade and the degree of genomic instability [19–21]. The absence of difference in the degree of genomic abnormalities between KRAS mutant and wt tumours is interesting, as both our transcriptional and genomic data imply enhanced activity of genes involved in chromosome structure and organisation in KRAS mutant tumours. We recognise that there were a small number of KRAS mutant tumours available for comparison and this may have limited our analysis.
Consistent with previously reported studies [9–16], the major differences in copy number and gene expression profiles between AC and SCC of the lung involved chromosome 3q. The strong independent correlation with amplification and over-expression at this locus suggests a causal relationship in SCC for genes in this region which warrant further investigation. These include TP73L, a gene extensively implicated in SCC, whose expression was most strongly correlated with the SCC phenotype, and which has been previously reported to be a putative oncogene [22–32]. While the role of TP73L in squamous cell pathogenesis remains unclear, a recent study of SiRNA mediated TP73L inhibition in SCC resulted in reduced cell survival with maintenance of squamous characteristics . These results suggest that TP73L is important in SCC cell survival. Other genes previously shown to be over-expressed in SCC were included in our differentiating list (CSTA [33, 34], FGFBP1), and warrant functional validation.
The differential copy number and expression levels between AC and SCC of TNFSF10/TRAIL and ABCC5, which have roles in apoptosis and chemoresistance respectively, may have implications for treatment of NSCLC. Recently published clinical data  suggest there are histotype-specific differences in response to systemic therapies. Validation of the differential activity of the roles of these genes and sensitivity to conventional and novel chemotherapeutic agents may be an area for future research.
Published data on the relationship between TP53 mutations and histotype in NSCLC is conflicting [36, 37]. SCC were associated with more frequent TP53 mutations than AC in our dataset. Cigarette smoking is a causal factor for both SCC phenotype and TP53 mutation . However, we also observed a correlation between 3q amplification and TP53 mutations in AC samples. This suggests that the apparent association between TP53 mutations and SCC may be mediated by the relationship between TP53 mutations and amplification of regions of 3q. We caution that this study is underpowered to draw strong conclusions regarding the role of TP53 in NSCLC pathogenesis.
Tumours possessing mutations of KRAS express genes playing key roles in cell growth, chromosome organisation and gene regulation. As previously reported, we identified amplification and over-expression of COX2 in KRAS mutant tumours. KRAS mutant tumours did not demonstrate mutations in EGFR consistent with previous reports in both NSCLC and colon carcinoma which suggest that KRAS mutations predict resistance to EGFR antagonists [29, 38–41]. Several reports link NCAM1 to Ras-dependent activation of ERK MAPK's [42, 43]. Reduced copy number and expression of NCAM1 in tumours bearing KRAS mutations, as seen in our data, has not previously been reported. Further research into a KRAS mutation profile may yield simple and reliable immunohistochemical markers of KRAS mutation, thereby significantly reduce the cost of determining KRAS status in clinical practice.
The gene expression profile observed in 'metastatic' tumours is consistent with a growing body of literature implicating deregulated protein synthesis in the development and metastatic potential of human cancers [44, 45]. Increased mRNA translation is a critical downstream function of many cancer related genes, and many gene products with roles in metastasis are not mutated but inappropriately expressed in malignant cells (e.g. VEGF, c-Myc, fos, Her2Neu, PDGF) . Opportunities for therapeutic intervention currently in development include oncolytic viruses that require deregulated protein translation for their replication , or agents that inhibit mTOR, an integral factor in protein translation (eg. temsirolimus (CCI-779), everolimus (RAD001) and deforolimus (AP23573)).
While the small number of recurrences and deaths due to NSCLC in our tumour-set makes it difficult to draw strong conclusions, transcriptional profiles linked to tumour recurrence suggest KRAS pathway activation. This may be due to a higher proportion of AC and LCC vs. SCC in the 'recurrent' group. Other regions of copy number change demonstrate genomic gains in the region of c-Myc and losses in the region of p16 (INK4a, CDNK2A) in recurrent or non-survivor tumours, supporting a prognostic association of the Myc:CDNK2A ratio in NSCLC, as has been described in head and neck SCC . Specific genes linked to recurrence or survival include SMARCA2 (implicated in the regulation of gene expression cell cycle control and oncogenesis), MINK (linked to the JNK MAP kinase pathway)  and RECK, which has putative roles in the suppression of tumour growth, invasion, angiogenesis and metastasis . KRAS mutation has been associated with reduced expression of RECK in NSCLC , consistent with the clinical observation of poor outcome in patients with KRAS mutation bearing NSCLC. Activation of the Ras pathway may reduce RECK expression and thereby increase tumour recurrence. Importantly, our prognostic gene signature was validated in an independent test set, suggesting that these findings may eventually yield prognostic markers in resected early-stage NSCLC to better select patients for adjuvant treatments.
Several molecular alterations have been identified in association with NSCLC histotype, KRAS mutation, TP53 mutation, metastatic potential, disease recurrence and survival. Although the size of the current study is small, our findings are in many cases consistent with those of previous studies, and have been validated in the case of the prognostic classifier in an independent test set. In addition, several novel molecular changes associated with clinically relevant endpoints have been demonstrated. It is hoped that these results will contribute to identifying new predictive markers and targets for novel therapies to improve treatment selection and better outcomes for patients with this deadly disease.
Maione P, et al: Combining targeted therapies and drugs with multiple targets in the treatment of NSCLC. Oncologist. 2006, 11: 274-84. 10.1634/theoncologist.11-3-274. 3
Morgan T, et al: Molecular profiling of Giant Cell Tumour of Bone and Osteoclastic Localization of Ligand for Receptor Activator of Nuclear Factor kB. American Journal of Pathology. 2005, 167: 117-28. 10.1016/S0002-9440(10)62959-8.
Weiss MM, et al: Genomic alterations in primary gastric adenocarcinomas correlate with clinicopathological characteristics and survival. Cell Oncol. 2004, 26: 307-17. 5-6
Krypuy M, et al: High resolution melting for mutation scanning of TP53 exons 5-8. BMC Cancer. 2007, 7: 10.1186/1471-2407-7-168. 168 PMID 17764544
Krypuy M, et al: High resolution melting analysis for the rapid and sensitive detection of mutations in clinical samples: KRAS codon 12 and 13 mutations in non-small cell lung cancer. BMC Cancer. 2006, 6: 10.1186/1471-2407-6-295. 295 PMID 17184525
Smyth G: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Statistical applications in genetics and molecular biology. 2004, Article 3
Smyth G, et al: The use of with-in array replicate spots for assessing expression in microarray experiments. Bioinformatics. 2005, 2067
Olshen A, et al: Circular binary segmentation for the analysis of array-based DNA copy number data. Biostatistics. 2004, 5: 557-72. 10.1093/biostatistics/kxh008.
Petersen I, et al: Patterns of chromosomal imbalances in adenocarcinoma and squamous cell carcinoma of the lung. Cancer Res. 1997, 57: 2331-5. 12
Bjorkqvist A, et al: DNA gains in 3q occur frequently in squamous cell carcinoma of the lung, but not in adenocarcinoma. Genes Chromosomes Cancer. 1998, 22: 79-82. 10.1002/(SICI)1098-2264(199805)22:1<79::AID-GCC11>3.0.CO;2-D. 1
Pei J, et al: Genomic imbalances in human lung adenocarcinomas and squamous cell carcinomas. Genes Chromosomes Cancer. 2001, 31: 282-87. 10.1002/gcc.1145. 3
Massion P, et al: Genomic copy number analysis of non-small cell lung cancer using array comparative genomic hybridisation: implications of the phosphatidylinositol 3-kinase pathway. Cancer Res. 2002, 62: 3636-40. 13
Massion P, et al: Early involvement of the phosphatidylinositol 3-kinase/Akt pathway in lung cancer progression. Am J Respir Crit Care Med. 2004, 170: 1088-94. 10.1164/rccm.200404-487OC. 10
Kim T, et al: Genome-wide screening of genomic alterations and their clinicopathologic implications in non-small cell lung cancers. Clin Cancer Res. 2005, 11: 8235-42. 10.1158/1078-0432.CCR-05-1157. 23
Tonon G, et al: High-resolution genomic profiles of human lung cancer. Proc Nat Ac Science USA. 2005, 9625-30. 10.1073/pnas.0504126102.
Garnis C, et al: High resolution analysis of non-small cell lung cancer cell lines by whole genome tiling path array CGH. Int J Cancer. 2006, 118: 1556-64. 10.1002/ijc.21491.
Graff J, Zimmer S: Translational control and metastatic progression: enhanced activity of the mRNA cap-binding protein eIF-4E selectively enhances translation of metastasis-related mRNAs. Clinical and Experimental Metastasis. 2003, 20: 265-73. 10.1023/A:1022943419011.
Bilanges B, Stokoe D: Mechanisms of translational deregulation in human tumors and therapeutic intervention strategies. Oncogene. 2007, 5973-5990. 10.1038/sj.onc.1210431.
Sahin AA, et al: Flow cytometric analysis of the DNA content of non-small cell lung cancer. Ploidy as a significant prognostic indicator in squamous cell carcinoma of the lung. Cancer. 1990, 65: 530-7. 10.1002/1097-0142(19900201)65:3<530::AID-CNCR2820650325>3.0.CO;2-6. 3
Granone P, et al: DNA flow cytometric analysis in patients with operable non-small cell lung carcinoma. Eur J Cardiothorac Surg. 1993, 7: 351-5. 10.1016/1010-7940(93)90065-J. 7
Kasprzyk M, et al: Prognostic value of DNA ploidy: 5-year follow-up of patients with resectable squamous cell carcinoma (SCC) of the lung. Lung Cancer. 2006, 51: 201-6. 10.1016/j.lungcan.2005.10.014. 2
Choi H, et al: Differential expression of p53 gene family members p63 and p73 in head and neck squamous tumorigenesis. Human Pathology. 2002, 33: 158-64. 10.1053/hupa.2002.30722. 2
Barbieri C, Barton C, Pietenpol J: Delta Np63 alpha expression is regulated by the phosphoinositide 3-kinase pathway. The journal of biological chemistry. 2003, 51: 51408-14. 10.1074/jbc.M309943200.
Glickman J, et al: Expression of p53-related protein p63 in the gastrointestinal tract and in esophageal metaplastic and neoplastic disorders. Human Pathology. 2001, 32: 1157-65. 10.1053/hupa.2001.28951. 11
Di Como C, et al: p63 expression profiles in human normal and tumour tissues. Clin Cancer Res. 2002, 8: 494-501. 2
Pruneri G, et al: p63 in laryngeal squamous cell carcinoma: evidence for a role of TA-p63 down-regulation in tumorigenesis and lack of prognostic implications of p63 immunoreactivity. Lab Invest. 2002, 82: 1327-34. 10
Massion P, et al: Significance of p63 amplification and overexpression in lung cancer development and prognosis. Cancer Res. 2003, 63: 7113-21. 21
Au N, et al: P63 expression in lung carcinoma: a tissue microarray study of 408 cases. Appl Immunohistochem Mol Morphol. 2004, 12: 240-47. 3
Van Cutsem E, Nowacki M, Lang I: Randomized phase III study of irinotecan and 5FU/FA with or without cetuximab in the first-line treatment of patients with metastatic colorectal cancer: The CRYSTAL trial. J Clin Pathol. 2007, 25: 164s abstr 4000
Thurfjell N, et al: Complex p63 mRNA isoform expression patterns in squamous cell carcinoma of the head and neck. Int J Oncol. 2004, 25: 27-35. 1
Nylander K, et al: Differential expression of p63 isoforms in normal tissues and neoplastic cells. J Pathol. 2002, 198: 417-27. 10.1002/path.1231. 4
De Young MP, et al: Tumor-specific p73 up-regulation mediates p63 dependence in squamous cell carcinoma. Cancer Res. 2006, 66: 9362-8. 10.1158/0008-5472.CAN-06-1619. 19
Ebert E, et al: Expression of cysteine protease inhibitors stefin A, stefin B, and cystatin C in human lung tumor tissue. Adv Exp Med Biol. 1997, 421: 259-65.
Werle B, et al: Cystatins in non-small cell lung cancer: tissue levels, localization and relation to prognosis. Oncol Rep. 2006, 16: 647-55. 4
Scagliotti GV, et al: Phase III study comparing cisplatin plus gemcitabine with cisplatin plus pemetrexed in chemotherapy-naive patients with advanced-stage non-small-cell lung cancer. J Clin Oncol. 2008, 26: 3543-51. 10.1200/JCO.2007.15.0375. 21
Takahashi T, et al: p53: a frequent target for genetic abnormalities in lung cancer. Science. 1989, 246: 491-4. 10.1126/science.2554494. 4929
Sekido Y, Fong K, Minna J: Molecular biology of lung cancer. [book auth.] V DeVita, S Hellman and S Rosenberg. Cancer - principles and practice of oncology. s.l. 2001, Lippincott Williams & Wilkins, 1: 917-25.
Pao W, et al: KRAS mutations and primary resistance of lung adenocarcinomas to gefitinib or erlotinib. PloS Med. 2005, 2: c73-10.1371/journal.pmed.0020073. 3
Peeters M, et al: Association of progression-free survival, overall survival, and patient-reported outcomes by skin toxicity and KRAS status in patients receiving panitumumab monotherapy. Cancer. 2009, 115: 1544-54. 10.1002/cncr.24088. 7
Bokemeyer C, et al: Fluorouracil, leucovorin, and oxaliplatin with and without cetuximab in the first-line treatment of metastatic colorectal cancer. J Clin Oncol. 2009, 27: 663-71. 10.1200/JCO.2008.20.8397. 5
Amado RG, et al: Wild-type KRAS is required for panitumumab efficacy in patients with metastatic colorectal cancer. J Clin Oncol. 2008, 26: 1626-34. 10.1200/JCO.2007.14.7116. 10
Schmid RS, et al: NCAM stimulates the Ras-MAPK pathway and CREB phosphorylation in neuronal cells. J Neurobiol. 1999, 38: 542-58. 10.1002/(SICI)1097-4695(199903)38:4<542::AID-NEU9>3.0.CO;2-1. 4
Son H, et al: Reciprocal actions of NCAM and tPA via a Ras-dependent MAPK activation in rat hippocampal neurons. Biochem Biophys Res Commun. 2002, 298: 262-8. 10.1016/S0006-291X(02)02453-1. 2
Ramaswamy S, et al: A molecular signature of metastasis in primary solid tumours. Nature Genetics. 2003, 33: 49-54. 10.1038/ng1060.
Weigelt B, et al: Gene expression profiles of primary breast tumours maintained in distant metastases. Proc Nat Acad Sciences USA. 2003, 100: 15901-05. 10.1073/pnas.2634067100. 26
Akervall J, et al: The gene ratios c-MYC:cyclin-dependent kinase (CDK)N2A and CCND1:CDKN2A correlate with poor prognosis in squamous cell carcinoma of the head and neck. Clin Cancer Res. 2003, 9: 1750-5. 5
Hu Y, et al: Identification and functional characterization of a novel human Misshapen/Nck interacting kinase-related kinase, hMINK-beta. J Biol Chem. 2004, 54387-54397. 10.1074/jbc.M404497200.
Oh J, et al: The membrane-anchored MMP inhibitor RECK is a key regulator of extracellular matrix integrity and angiogenesis. Cell. 2001, 107: 789-800. 10.1016/S0092-8674(01)00597-9.
Chang HC, Cho CY, Hung WC: Downregulation of RECK by promoter methylation correlates with lymph node metastasis in non-small cell lung cancer. Cancer Sci. 2007, 98: 169-73. 10.1111/j.1349-7006.2006.00367.x. 2
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/11/93/prepub
The authors would like to acknowledge the assistance of Michael Krypuy from the Department of Pathology, Peter MacCallum Cancer Centre in undertaking high resolution melting analysis for KRAS mutations. Funding for this project was received from The University of Melbourne (scholarship for Dr G Newnham), the Australasian Society of Cardiac and Thoracic Surgeons Foundation, and The National Health and Medical Research Council of Australia.
The authors declare that they have no competing interests.
GN performed DNA and RNA extraction and transcriptional microarray studies, and drafted the manuscript. MC, SM and GW participated in the design of the study and co-ordinated collection of clinical data. AD performed and supervised mutation analyses. HD performed mutational analyses. JL and NT performed statistical analysis. KO performed pathology review. DT and MC conceived the idea of the study. DT supervised molecular studies and assisted in the preparation of the manuscript. All authors read and approved the final copy of the manuscript.