- Short Report
- Open Access
Global gene expression profiling of oral cavity cancers suggests molecular heterogeneity within anatomic subsites
BMC Research Notesvolume 1, Article number: 113 (2008)
Oral squamous cell carcinoma (OSCC) is a frequent neoplasm, which is usually aggressive and has unpredictable biological behavior and unfavorable prognosis. The comprehension of the molecular basis of this variability should lead to the development of targeted therapies as well as to improvements in specificity and sensitivity of diagnosis.
Samples of primary OSCCs and their corresponding surgical margins were obtained from male patients during surgery and their gene expression profiles were screened using whole-genome microarray technology. Hierarchical clustering and Principal Components Analysis were used for data visualization and One-way Analysis of Variance was used to identify differentially expressed genes. Samples clustered mostly according to disease subsite, suggesting molecular heterogeneity within tumor stages. In order to corroborate our results, two publicly available datasets of microarray experiments were assessed. We found significant molecular differences between OSCC anatomic subsites concerning groups of genes presently or potentially important for drug development, including mRNA processing, cytoskeleton organization and biogenesis, metabolic process, cell cycle and apoptosis.
Our results corroborate literature data on molecular heterogeneity of OSCCs. Differences between disease subsites and among samples belonging to the same TNM class highlight the importance of gene expression-based classification and challenge the development of targeted therapies.
Head and neck squamous cell carcinomas (HNSCC) rank among the top ten most common cancers worldwide, with increasing rates in certain areas of the world . They can occur at different subsites, often invading more than one, each one with their own particular problems regarding management. In this study, we focused on oral squamous cell carcinoma (OSCC), a frequent neoplasm, which is usually aggressive and has unpredictable biological behavior and unfavorable prognosis. The accumulation of genetic alterations during oral tumorigenesis, leading to qualitative and quantitative changes in gene expression, is currently known but still largely unexplored [2, 3]. Unlike estrogen receptor or HER2 in breast cancer, no biomarkers are currently available for HNSCC prognosis, which depends largely on the stage at presentation, with the most important prognostic factor being the presence of neck node metastases . Improvements in specificity and sensitivity of diagnosis and in disease monitoring depend on the elucidation of the biological and molecular mechanisms underlying tumor development. In this context, large-scale transcriptome analysis may be used to assess patterns of gene expression, providing an opportunity to investigate the complex cascade of molecular events leading to tumor development and progression [5, 6].
Two recent publications have examined the global gene expression profiles of OSCC, both addressing the importance of understanding the molecular complexity of this malignancy [7, 8]. Unlike those authors, we did not focus on molecular differences between OSCC and normal mucosa, but we compared gene expression profiles of OSCC samples from different anatomic subsites. The biologic behavior of OSCC cannot be predicted and clinical reports suggest that molecular heterogeneity of anatomic subsites could play an important role in this scenario. For instance, cancer of the tongue seems to grow faster than other oral cavity cancers, with cervical metastases occurring more frequently in such cancers. Moreover, significant differences between floor of the mouth and tongue cancers in response to combined surgery and radiotherapy have been reported .
In order to evaluate molecular differences between OSCC anatomic subsites, we investigated samples of primary tongue and floor of the mouth tumors, and their corresponding surgical margins, in respect to their gene expression profile by means of DNA microarrays.
We selected samples of primary tongue and floor of the mouth tumors, and their corresponding surgical margins, from 9 male patients during surgery for squamous cell carcinomas. Samples consisted of either small tumors but already with evidences of metastasis (T2N+), or larger, non-metastatic samples (T3N0). Patients had a history of tobacco and alcohol abuse and no prior cancer treatment. A full description of the clinical data, including tumor stage is provided in Additional File 1. All participants provided written consent and the research protocol was approved by review boards of all involved institutions. Analysis of haematoxylin and eosin-stained sections indicated that each OSCC sample contained at least 70% tumor cells and the corresponding surgical margins were reported "tumor-free". Total RNA was extracted with TRIZOL Reagent (Invitrogen). Tissues adjacent to resected tumors were grouped into two pools of RNA for microarray analysis, corresponding to T2N+ and T3N0 samples, due to the limited amounts of tumor-free margin available in most cases.
Microarray experiments were carried out using CodeLink Whole Genome Bioarrays (GE Healthcare) and arrays were scanned on a GenePix 4000B Array Scanner (Axon Instruments), according to the recommended scanning procedures and settings. The data were treated with CodeLink feature extraction software v.4.0. A normalized signal for each transcript was obtained through quantile normalization . The array design and raw data files are available at the Gene Expression Omnibus database (GEO) under the accession number GSE9792. For global gene expression visualization, we used hierarchical clustering and Principal Components Analysis (PCA). Hierarchical clustering was performed using the Euclidean distance and the average linkage algorithm. One-way ANOVA was used to identify significant differences in the dataset. All the above-mentioned analyses were generated using Partek® software version 6.3 Copyright© 2007 Partek Inc., St. Louis, MO, USA. Functional annotations of differentially expressed genes were performed using Database for Annotation, Visualization and Integrated Discovery  using the parameters Gene Ontology (GO) Molecular Process term level 5 and KEGG Pathways (Kyoto Encyclopedia of Genes and Genomes).
In order to corroborate our findings, two publicly available microarray datasets on OSCC were used in addition to our data (GEO accession number GSE3524 and ArrayExpress accession number TABM302). Both experiments were carried out using Affymetrix HG-U133 GeneChips; GSE3524 corresponds to the profiling of 7 tongue and 9 floor of the mouth samples, while TABM302 corresponds to the profiling of 18 tongue samples and 10 samples reported as oral cavity. Statistical analysis and functional annotation of the differentially expressed genes followed the procedures previously described for our dataset (Additional Files 2 and 3).
Results and Discussion
By analyzing the dendogram resulting from hierarchical clustering and the PCA, we observed that OSCC samples did not group according to their pathological stages (Figures 1 and 2, respectively). In agreement with previous reports [14, 15], two major clusters reflected differences in global gene expression between non-tumoral tissues (cluster 1) and tumor samples (cluster 2). Tumor samples clustered mostly, but not exclusively, according to disease anatomic subsites. This result suggests molecular heterogeneity within TNM classes and corroborates previous data published by Ziober et al.  and Chung et al. In the later study, there were no consistent differences in gene expression among HNSCC subsites, but samples from oral cavity seemed more heterogeneous than samples from other sites when expression patterns were evaluated. Few studies have aimed to understand the molecular background of OSCC. In the context of molecular target therapies, drugs are designed to interact with specific molecules present in certain types or subtypes of tumors. Due to the heterogeneity within HNSCC, a thoroughly comprehension of the molecular characteristics underlying its subsites is needed before efficient therapy can be achieved.
Differential gene expression profiles
Using hierarchical clustering and PCA, two main groups of tumor samples (2a and 2b) were observed (Figure 1). Each group contained, mostly, either tongue or floor of the mouth samples. Differentially expressed genes between the two clusters were identified using ANOVA (Additional File 4). These genes (total number 1579) were functionally annotated using Gene Ontology (GO) and KEGG terms in an attempt to understand major molecular differences between the groups (Tables 1 and 2). Due to our small sample number, our results were compared with two additional OSCC microarray datasets, GSE3524 and TABM302. Despite broad differences in experimental design, probe sets, and possibly genetic background, some findings were consistent and are pointed out below.
Extensive differences in gene expression between clusters 2a and 2b were observed in respect to translation and mitosis-related GO terms (GO:0006412:translation, GO:0000087:M phase of mitotic cell cycle; GO:0007067:mitosis; GO:0000279:M phase). Translation initiation is regulated in response to mitogenic stimulation, thus coupled with cell cycle progression and cell growth. Several alterations in translational control occur in cancer, but there is still much to be discovered about their role in cancer development and progression . Changes in the expression or availability of components of the translational machinery can lead to global changes, such as an increase in the overall rate of protein synthesis and translational activation of the mRNA molecules involved in cell growth and proliferation. In agreement, differences between the two clusters were enriched in GO terms Ribonucleotide Biosynthetic Process, RNA Processing and mRNA Metabolic Process as well as in KEGG Proteasome, Cell Cycle, TGFB Signaling and Ribosome Pathways. Components of the translation machinery and related pathways represent good targets for cancer therapy [18, 19]. The differential expression of a number of genes involved in such processes suggests that the efficacy of oral cancer therapy should vary depending on the tumor subsite. Molecular mechanisms leading to differences in gene expression (i.e variations in mRNA sequences, changes in the availability of the translational machinery components, and activation of translation through abnormally activated signal transduction pathways) deserve to be investigated.
Differences in the expression of cell cycle progression-related genes were also observed in both GEO (GSE3524) and ArrayExpress (TABM302) datasets. In respect to the former, several genes associated with the KEGG MAPK Signaling Pathway presented differential expression when tongue samples were compared with floor of the mouth samples. Similar results were observed for GO terms concerning transcriptional regulation and apoptosis (Additional File 5). In agreement with our results, TABM302 dataset showed differential expression in genes associated with the GO terms Maintenance of Chromatin Architecture and DNA Packaging, both essential in cell cycle-related processes, and with the KEGG Ribosome Pathway (Additional File 6). Proliferation and cell death are also regulated by genes associated with the GO terms Small GTPase Mediated Signal Transduction and Regulation of Ras Protein Signal Transduction, both represented several times in the list of differentially expressed genes of the TABM302 dataset.
With respect to cell cycle and apoptosis, we also observed the over-expression in cluster 2b of three genes involved in the regulation of p27 phosphorylation during cell cycle progression (CKS1B, RBX1 and SKP1A) and genes from the MCM complex (minichromosomal maintenance genes), including MCM7 (Additional File 4). MCM7 protein is a marker for proliferation and it is upregulated in different tumors including neuroblastoma, prostate, cervical and hypopharyngeal carcinomas . SKP1A, MCM3, MCM6 and MCM7 were also differentially expressed in TABM302 dataset (Additional File 3), however, their over-expression was observed in tongue samples.
Widespread differences were also observed in oxidative phosphorylation pathways (Tables 2 and 3). Genes belonging to KEGG Oxidative Phosphorylation and Citrate Cycle (TCA cycle) Pathways were over-expressed in cluster 2b. Oxidative phosphorylation of ADP to ATP accompanies the oxidation of a metabolite through the operation of the respiratory chain. According to Warburg , alterations of respiratory machinery should result in compensatory increase in glycolytic ATP production and lead to carcinogenesis. In fact, malignant cells meet their energy (ATP) needs by glycolysis rather than through oxydative phosphorylation, probably an adaptation to hypoxia that develops as tumor grows . In accordance with our results, some authors have recently found a bioenergetic signature in several human tumors . Our observations suggest that distinct OSCC subsites differ on their glycolytic phenotype, a matter that deserves further investigation and may impact tumor progression management as well as therapy outcome.
Genes related to KEGG Oxidative Phosphorylation Pathway also showed differences in expression in GSE3524 dataset when subsites were compared; the same was observed for genes related to the GO term Organelle ATP Synthesis Coupled Electron Transport (Additional File 5). The expression of several cytochrome oxidases and ATPases was also distinct between subsites in TABM302 (Additional File 3).
KEGG Pathway analysis revealed genes over-expressed in cluster 2a which are associated with gap junctions and adherens junctions, as well as cytoskeleton organization and biogenesis (Table 4), including ACTN, DOCK1, GSN, ITG, LIMK, RAS, SOS, SSH and TIAM1. In agreement with our data, genes related to KEGG Regulation of Actin Cytoskeleton and Adherens Junction Pathways were also found differentially expressed in TABM302 (Additional File 6). Cancer therapy has focused on such pathways. For example, microtubule-targeted drugs interfere with the ability of the cancer cells to divide and multiply by disrupting microtubules . Gene-silencing may reduce expression of proteins involved in regulating the actin cytoskeleton, as LIMK, rendering cells more responsive to anticancer agents . Therefore, molecular differences between OSCC subsites may affect sensitivity to chemotherapy.
In addition to the data presented above, several genes over-expressed in cluster 2b have been previously associated with oral cancer, such as GSTM1, IL8, MGMT, NFKB1 and NUDT1; as well as with pancreatic cancer (APEX1, CDK4, FANCC, GSTM1, MGMT and SOD2), lung cancer (APEX1, CASP9, GSTM1, HSPA8, IL8, LIG1, MBD4, MGMT, MLH1, NUDT1 and SOD2), colorectal cancer (APEX1, BLM, GSTM1, IL8, MGMT, MLH1, NFKB1, NUDT1, RFC1, SOD2 and UMPS) and bladder cancer (APEX1, GSTM1, IL8, LIG1, MGMT and SOD2) (Additional File 4). We believe the differential expression between anatomic subsites of genes known to be involved in tumorigenesis deserves further investigation.
Gene expression profiling of lymph node positive and lymph node negative HNSCC has been addressed in the literature [14, 15, 26, 27]. Besides the differences in study designs, the lack of consistency concerning HNSCC gene expression signatures could be due to gene expression variation related to subsites. This is a preliminary study and as such presents results on a small number of samples. However, the idea of molecular heterogeneity within HNSCC is not new. We corroborate previously reported literature data and we suggest that differences between anatomic OSCC subsites could impact drug response and should be considered during the development of targeted therapies.
Shibuya K, Mathers CD, Boschi-Pinto C, Lopez AD, Murray CJ: Global and regional estimates of cancer mortality and incidence by site: II. Results for the global burden of disease 2000. BMC Cancer. 2002, 2: 37-10.1186/1471-2407-2-37.
Scully C, Field JK, Tanzawa H: Genetic aberrations in oral or head and neck squamous cell carcinoma (SCCHN): 1. Carcinogen metabolism, DNA repair and cell cycle control. Oral Oncol. 2000, 36 (3): 256-263. 10.1016/S1368-8375(00)00007-5.
Scully C, Field JK, Tanzawa H: Genetic aberrations in oral or head and neck squamous cell carcinoma 3: clinico-pathological applications. Oral Oncol. 2000, 36 (5): 404-413. 10.1016/S1368-8375(00)00023-3.
Greenberg JS, Fowler R, Gomez J, Mo V, Roberts D, El Naggar AK, Myers JN: Extent of extracapsular spread: a critical prognosticator in oral tongue cancer. Cancer. 2003, 97 (6): 1464-1470. 10.1002/cncr.11202.
Huang E, Cheng SH, Dressman H, Pittman J, Tsou MH, Horng CF, Bild A, Iversen ES, Liao M, Chen CM, et al: Gene expression predictors of breast cancer outcomes. Lancet. 2003, 361 (9369): 1590-1596. 10.1016/S0140-6736(03)13308-9.
van 't Veer LJ, Dai H, Vijver van de MJ, He YD, Hart AA, Mao M, Peterse HL, Kooy van der K, Marton MJ, Witteveen AT, et al: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415 (6871): 530-536. 10.1038/415530a.
Sticht C, Freier K, Knopfle K, Flechtenmacher C, Pungs S, Hofele C, Hahn M, Joos S, Lichter P: Activation of MAP kinase signaling through ERK5 but not ERK1 expression is associated with lymph node metastases in oral squamous cell carcinoma (OSCC). Neoplasia. 2008, 10 (5): 462-470.
Kondoh N, Ishikawa T, Ohkura S, Arai M, Hada A, Yamazaki Y, Kitagawa Y, Shindoh M, Takahashi M, Ando T, et al: Gene expression signatures that classify the mode of invasion of primary oral squamous cell carcinomas. Mol Carcinog. 2008
Ridge JA, Glisson BS, Horwitz EM, MN L: Cancers of the Head and Neck Region. Textbook of cancer management: a multidisciplinary approach. 2007, Melville, NY: PPR, Inc., 31-72. 10
Zelefsky MJ, Harrison LB, Fass DE, Armstrong J, Spiro RH, Shah JP, Strong EW: Postoperative radiotherapy for oral cavity cancers: impact of anatomic subsite on treatment outcome. Head Neck. 1990, 12 (6): 470-475. 10.1002/hed.2880120604.
Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19 (2): 185-193. 10.1093/bioinformatics/19.2.185.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4 (5): P3-10.1186/gb-2003-4-5-p3.
The Cancer Genome Anatomy Project. [http://cgap.nci.nih.gov/Pathways/]
Mendez E, Cheng C, Farwell DG, Ricks S, Agoff SN, Futran ND, Weymuller EA, Maronian NC, Zhao LP, Chen C: Transcriptional expression profiles of oral squamous cell carcinomas. Cancer. 2002, 95 (7): 1482-1494. 10.1002/cncr.10875.
Ziober AF, Patel KR, Alawi F, Gimotty P, Weber RS, Feldman MM, Chalian AA, Weinstein GS, Hunt J, Ziober BL: Identification of a gene signature for rapid screening of oral squamous cell carcinoma. Clin Cancer Res. 2006, 12 (20 Pt 1): 5960-5971. 10.1158/1078-0432.CCR-06-0535.
Chung CH, Parker JS, Karaca G, Wu J, Funkhouser WK, Moore D, Butterfoss D, Xiang D, Zanation A, Yin X, et al: Molecular classification of head and neck squamous cell carcinomas using patterns of gene expression. Cancer Cell. 2004, 5 (5): 489-500. 10.1016/S1535-6108(04)00112-6.
Meric F, Hunt KK: Translation initiation in cancer: a novel target for therapy. Mol Cancer Ther. 2002, 1 (11): 971-979.
Hidalgo M, Rowinsky EK: The rapamycin-sensitive signal transduction pathway as a target for cancer therapy. Oncogene. 2000, 19 (56): 6680-6686. 10.1038/sj.onc.1204091.
Palakurthi SS, Fluckiger R, Aktas H, Changolkar AK, Shahsafaei A, Harneit S, Kilic E, Halperin JA: Inhibition of translation initiation mediates the anticancer effect of the n-3 polyunsaturated fatty acid eicosapentaenoic acid. Cancer Res. 2000, 60 (11): 2919-2925.
Honeycutt KA, Chen Z, Koster MI, Miers M, Nuchtern J, Hicks J, Roop DR, Shohet JM: Deregulated minichromosomal maintenance protein MCM7 contributes to oncogene driven tumorigenesis. Oncogene. 2006, 25 (29): 4027-4032. 10.1038/sj.onc.1209435.
Warburg O: On the origin of cancer cells. Science. 1956, 123 (3191): 309-314. 10.1126/science.123.3191.309.
Assaily W, Benchimol S: Differential utilization of two ATP-generating pathways is regulated by p53. Cancer Cell. 2006, 10 (1): 4-6. 10.1016/j.ccr.2006.06.014.
Cuezva JM, Krajewska M, de Heredia ML, Krajewski S, Santamaria G, Kim H, Zapata JM, Marusawa H, Chamorro M, Reed JC: The bioenergetic signature of cancer: a marker of tumor progression. Cancer Res. 2002, 62 (22): 6674-6681.
Jordan MA: Mechanism of action of antitumor drugs that interact with microtubules and tubulin. Curr Med Chem Anticancer Agents. 2002, 2 (1): 1-17. 10.2174/1568011023354290.
Mendez E, Fan W, Choi P, Agoff SN, Whipple M, Farwell DG, Futran ND, Weymuller EA, Zhao LP, Chen C: Tumor-specific genetic expression profile of metastatic oral squamous cell carcinoma. Head Neck. 2007, 29 (9): 803-814. 10.1002/hed.20598.
Roepman P, Wessels LF, Kettelarij N, Kemmeren P, Miles AJ, Lijnzaad P, Tilanus MG, Koole R, Hordijk GJ, Vliet van der PC, et al: An expression profile for diagnosis of lymph node metastases from primary head and neck squamous cell carcinomas. Nat Genet. 2005, 37 (2): 182-186. 10.1038/ng1502.
We are grateful to the authors of microarray datasets GEO GSE3524 and ArrayExpress TABM302. The GENCAPO group acknowledges the financial support from Fundação de Amparo à Pesquisa do Estado de São Paulo/FAPESP (Grants 04/12054-9 and 04/14029-1) and Associação Beneficente Alzira Denise Hertzog Silva/ABADHS, and the researcher fellowships from Conselho Nacional de Pesquisas (CNPq), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Instituto Israelita de Ensino e Pesquisa Albert Einstein, and The Ludwig Institute for Cancer Research.
We also thank to the Head and Neck Genome Project GENCAPO members for data collection, study design, comments on the draft: Cury PM, Frizzera APZ (Departamento de Patologia, Faculdade de Medicina de São José do Rio Preto); de Carvalho MB, Silva AMA, Amar A, Barbieri RB, Bastos AU, Carvalho-Neto PB, Casemiro AF, Chedid H, Chiappini PBO, Correia LA, Costa ACW, Curioni OA, Franzi SA, Gazito D, Gutierres AP, Lehn CN, Martins AE, Mercante AMC, Porsani AF, Rapoport A, Rossi L, Santos M, Souza TB, Takamori JT (Hospital Heliópolis, São Paulo); Dias-Neto E, Ojopi EPB, Dias THG (Departamento e Instituto de Psiquiatria, Faculdade de Medicina, USP, São Paulo); Figueiredo DLA, Mamede RCM (Serviço de Cirurgia de Cabeça e Pescoço, Faculdade de Medicina de Ribeirão Preto, USP); Fukuyama EE, Góis-Filho JF, Cerione M, Cicco R, Settani F, Valentim PJ, Yamagushi F, Cominato ML, Mendes GS, Paiva R, Silva MJ (Serviço de Cirurgia de Cabeça e Pescoço, Instituto do Câncer Arnaldo Vieira de Carvalho, São Paulo); Leopoldino AM, Silva FAM (Departamento de Análises Clínicas, Toxicológicas e Bromatológicas, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, USP); Michaluart-Junior P, Moyses RA, Arap SS, Araújo NSS, Araújo-Filho V, Brandão LG, Cernea CR, Durazzo M, Ferraz AR, Gallo J, Guimarães PEM, Magalhães RP, Montenegro FLM, Silva-Filho GB, Smith RB, Stabenow E, Tavares MR, Turcano R, Volpi EM, Ramos O, Silva C (Departamento de Cirurgia de Cabeça e Pescoço, Faculdade de Medicina, USP, São Paulo); Moreira-Filho CA (Departamento de Pediatria, Faculdade de Medicina, USP, São Paulo); Nóbrega FG, Nóbrega MP, Canto AL, Macarenco R, Meneses C, Correa PMS, Bogossian AP (Departamento de Biociências e Diagnóstico Bucal, Faculdade de Odontologia, UNESP, São José dos Campos); Nunes FD, Souza SCOM, Rodini CO, Xavier FCA (Departamento de Estomatologia, Faculdade de Odontologia, USP, São Paulo); Okamoto OK (Departamento de Neurologia e Neurocirurgia, UNIFESP, São Paulo); Serafini LN (Departamento de Patologia, Faculdade de Medicina de Ribeirão Preto, USP); Severino P (Instituto de Ensino e Pesquisa Albert Einstein, São Paulo); Silva Jr WA, Brandão RM, Kaneto CM, Pinheiro DG, Santos ARD, Silva IT, Tarlá MVC (Departamento de Genética, Faculdade de Medicina de Ribeirão Preto, USP); Silveira NJF (Ciências da Computação, UNIVAP, São José dos Campos); Tajara EH, Rodrigues-Lisoni FC, Rodrigues RV, Polachini GM, Vidotto A, Cunha BR, Carmona-Raphe J (Departamento de Biologia Molecular, Faculdade de Medicina de São José do Rio Preto;); Wünsch-Filho V, Costa A, Figueiredo RO, Fortes CS, Inamine R, López RVM, Rodrigues AN (Departamento de Epidemiologia, Faculdade de Saúde Pública, USP, São Paulo); Zago MA (Departamento de Clínica Médica, Faculdade de Medicina de Ribeirão Preto, USP); SP, Brazil.
The authors declare that they have no competing interests.
PS carried out hybridization experiments, data analysis and drafted the manuscript; AMA coordinated sample collection and carried out clinical data analysis; PM carried out clinical data analysis for sample selection; OKO helped with data analysis and the study design; FDM participated in the study design and carried out clinical data analysis; GENCAPO team members were responsible for sample collection, initial on-site sample processing, provided the pathological analysis of the cases, obtained the informed consent and discussed the findings. CAMF participated in the study design and coordination; EHT conceived the study, carried out data analysis and coordination. All authors read and approved the final manuscript.