- Research article
- Open Access
Clinical relevance of breast cancer-related genes as potential biomarkers for oral squamous cell carcinoma
BMC Cancer volume 14, Article number: 324 (2014)
Squamous cell carcinoma of the oral cavity (OSCC) is a common cancer form with relatively low 5-year survival rates, due partially to late detection and lack of complementary molecular markers as targets for treatment. Molecular profiling of head and neck cancer has revealed biological similarities with basal-like breast and lung carcinoma. Recently, we showed that 16 genes were consistently altered in invasive breast tumors displaying varying degrees of aggressiveness.
To extend our findings from breast cancer to another cancer type with similar characteristics, we performed an integrative analysis of transcriptomic and proteomic data to evaluate the prognostic significance of the 16 putative breast cancer-related biomarkers in OSCC using independent microarray datasets and immunohistochemistry. Predictive models for disease-specific (DSS) and/or overall survival (OS) were calculated for each marker using Cox proportional hazards models.
We found that CBX2, SCUBE2, and STK32B protein expression were associated with important clinicopathological features for OSCC (peritumoral inflammatory infiltration, metastatic spread to the cervical lymph nodes, and tumor size). Consequently, SCUBE2 and STK32B are involved in the hedgehog signaling pathway which plays a pivotal role in metastasis and angiogenesis in cancer. In addition, CNTNAP2 and S100A8 protein expression were correlated with DSS and OS, respectively.
Taken together, these candidates and the hedgehog signaling pathway may be putative targets for drug development and clinical management of OSCC patients.
Oral squamous cell carcinoma (OSCC) is the most common malignancy form in the head and neck region, accounting for about 260,000 new cases and 124,000 OSCC-related deaths worldwide annually [1, 2]. In western countries, the etiology of some newly diagnosed primary SCCs of the head and neck has shifted from tobacco and alcohol abuse to human papillomavirus (HPV) infections, possibly as a result of changes in sexual practices [3–8]. Despite aggressive treatment modalities, 5-year survival rates for advanced head and neck cancers have remained low and relatively unchanged (about 50-60%) for several decades, partially due to early locoregional recurrences within 2 years of initial treatment . There is therefore a pressing need for molecular predictors that enable earlier detection of the disease, describe tumor behavior, and improve risk assessment to identify patients at risk for recurrence and OSCC-related death.
Molecular profiling has become a common and effective method for cancer gene discovery and classification of cancer. Almost a decade ago, Chung et al. identified four intrinsic subtypes for head and neck squamous cell carcinoma with clinical and biological implications . Consequently, one of the subtypes with the most unfavorable prognosis also displayed strikingly similar transcriptional patterns with the breast carcinoma basal-like phenotype and lung squamous cell carcinoma. Recently, additional evidence of shared cellular processes was found between breast carcinoma and oral squamous cell carcinoma, i.e. mechanisms for tumor lymphangiogenesis and metastasis to the regional lymph nodes as well as HER2/neu polymorphisms [11, 12]. These findings suggest that cancers derived from different sites of origin may perturb common signaling pathways and thereby display similar tumor characteristics . To test this hypothesis, we evaluated the prognostic potential of 16 putative prognostic biomarkers (AZGP1, BTG2, CBX2, CNTNAP2, DNALI1, LOC389033, NME5, PIP, S100A8, SCUBE2, SERPINA11, STC2, STK32B, SUSD3, UBE2C, and WHSC1L1) previously identified in breast carcinoma using oral squamous cell carcinoma [14, 15]. Interestingly, several of the putative biomarkers have been implicated in the carcinogenesis of more than one cancer form. We, and others, have also been able to show that AZGP1, S100A8, and STK32B, as well as CNTNAP2 are associated with the basal-like phenotype and lymph node metastasis, respectively [15–17].
Here, we investigated the prognostic potential of the gene expression signature in relation to clinical outcome, disease-specific survival (DSS) and/or overall survival (OS), in two steps. First, transcriptional levels for each gene were evaluated with respect to the clinical endpoints using publicly available Affymetrix one-channel microarray (n = 168) and Illumina RNASeq datasets (n = 198) for OSCC from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) repositories, respectively. Second, because correlation between mRNA/protein levels is frequently low, Cox proportional hazards models for DSS and OS were also calculated using immunohistochemical protein expression patterns from 43 OSCC cases together with established clinicopathological features (tumor size and node status or tumor size and age, respectively).
To evaluate the prognostic potential of the AZGP1, BTG2, CBX2, CNTNAP2, DNALI1, LOC389033, NME5, PIP, S100A8, SCUBE2, SERPINA11, STC2, SUSD3, STK32B, UBE2C, and WHSC1L1 genes in OSCC specimens, three patient cohorts were compiled primarily from squamous cell carcinomas of the oral cavity. External gene expression datasets and corresponding clinical information for Cohorts I-II were compiled from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) repositories, respectively. Cohort I included two Affymetrix U133 Plus 2.0 GeneChip datasets (GEO accession numbers GSE41613 and GSE42743) containing 168 OSCC samples (oropharynx samples were excluded from the analysis) . Cohort II consisted of normalized RNAseq by Expectation-Maximization (RSEM) gene datasets from 198 OSCC patients (oral cavity: buccal mucosa, floor of mouth, tongue), which were downloaded from The Broad Institute TCGA GDAC (http://gdac.broadinstitute.org/runs/stddata__2014_01_15/). Cohort III consisted of 43 OSCC cases originating from the oral cavity (buccal gingiva, floor of mouth, tongue), which had been diagnosed between 1997-2004 at Sahlgrenska University Hospital in Gothenburg, Sweden. All patients underwent diagnostic battery inclusive biopsy of the primary tumor, palpation of the neck, radiological examination with MRT and/or CT, and TNM classified according to the American Joint Committee on Cancer (AJCC) staging system. Surgical excision of the primary tumor and supraomohyoid neck dissection (SOHND) were performed. In total, 16 patients had cervical lymph node metastases (pN1) of which 5/16 patients had micrometastases (pNmic) as assessed using anti-human monoclonal cytokeratin AE1/AE3, and 27 patients were lymph node-negative (pN0). Lymph node-positive patients received post-operative radiotherapy to the neck, whereas pN0 and pNmic patients were followed up clinically. All patients were followed up for at least five years during which seven patients (16%) developed local and/or regional recurrence, including two patients with pN1, three patients with pNmic, and two patients with pN0 disease. Three of the five patients with micrometastases developed recurrence (60%), of which two (40%) died within three years due to OSCC-related causes. The clinicopathological features for Cohorts I-III are summarized in Table 1.
For Cohort III, 45 FFPE samples corresponding to the 43 patients were obtained from the Department of Pathology at Sahlgrenska University Hospital and used in immunohistochemistry experiments in accordance with the Declaration of Helsinki and approved by the Medical Faculty Research Ethics Committee (Gothenburg, Sweden). The ethics committee approved a waiver of written consent to use the tumor specimens in the study. Histological classification and TNM staging of the tumor specimens were performed according to the WHO classification and International Union Against Cancer (UICC), respectively [19, 20]. Optimal antibody dilutions and assay conditions were achieved for immunohistochemistry using OSCC as positive controls. Four micrometer full-face FFPE sections were pretreated using the Dako PTLink system (Dako, Carpinteria, CA, USA) and processed using the Dako Envision™ FLEX High pH Link Kit (pH 9) for p16, AZGP1, BTG2, CBX2, CNTNAP2, NME5, S100A8, SCUBE2, SERPINA11, STC2, SUSD3, STK32B, SUSD3, UBE2C, and WHSC1L1 as listed in Additional file 1: Table S1. Peroxidase-catalyzed diaminobenzidine was used as the chromogen, followed by hematoxylin counterstain. The slides were then rinsed with deionized water, dehydrated in absolute alcohol, followed by 95% alcohol, cleared in xylene, and mounted. H & E staining was performed on one FFPE section to facilitate histological assessment. The degree of lymphoplasmacytic infiltration (inflammatory infiltration) was classified as minimal (few inflammatory cells), moderate (1-2 mm margin), and strong (>2 mm margin) according to the density of inflammatory cells. Immunostaining was evaluated by a head and neck pathologist, blinded to patient clinical outcome, and scored as previously described using the semi-quantitative H-score method to calculate the sum of the percentage and intensity of positively stained tumor cells within the invasive tissue component (negative staining = 0; weak staining = 1+; moderate staining = 2+; strong staining = 3+). The H-score ranged from 0 to 300, where H-score = (1 ×%1+) + (2 ×%2+) + (3 ×%3+) . The X-tile software (version 3.6.1) was used to determine an H-score cut-off for positive staining by dichotomizing patients according to H-score value and clinical outcome, as listed in Additional file 1: Table S1 . FFPE specimens lacking an invasive tissue component were removed from the analysis. Each tumor specimen was scored once, where multiple FFPE sections representing the same tumor were averaged. Staining was evaluated in the invasive and peritumoral stromal/normal tissue components.
Statistical analyses were performed using a 0.05 P-value cutoff (two-sided) in R/Bioconductor (version 2.15.0). Putative prognostic biomarkers for OSCC were identified in two steps. First, the prognostic potential of aberrant biomarker gene expression was evaluated in external microarray and RNASeq datasets (Cohorts I-II). Then, predictive models for DSS and OS were developed using biomarker protein expression (Cohort III).
Evaluation of gene expression patterns for the 16-marker signature in external microarray and RNASeq datasets
Univariate Cox proportional hazard models were calculated for each gene using the endpoints disease specific-survival (DSS) and/or overall survival (OS). OSCC survival rates were defined as a) the period from initial diagnosis to OSCC-related death for DSS and b) period from initial diagnosis to death from any cause for OS. Data processing and Cox regression analysis of the Affymetrix one-channel microarray datasets (Cohort I) and normalized RNASeq RSEM values (Cohort II) were performed using Nexus Expression 3.0 (BioDiscovery).
Development of a predictive model for DSS and OS using protein expression
Survival rates (DSS and OS) at different protein expression levels were depicted with Kaplan-Meier curves and tested with log-rank test. The relationship between clinicopathological features and protein expression was evaluated using two-tailed Fisher’s exact test. Multivariate analysis was conducted using the Cox proportional hazard model for DSS or OS with stepwise selection to assess the predictive strength and additive accuracy of protein expression after adjusting for established clinicopathological features (tumor size and node status or tumor size and age, respectively). A concordance index (C-index) for the time-dependent area under the ROC curve (AUC (t)) was calculated to assess model predictive performance, varying from C-index = 0.5 (no predictive power) to C-index = 1 (perfect prediction).
Prognostic potential of the molecular biomarkers in external gene expression microarray and RNASeq datasets
In previous work, we showed the clinical significance of 16 candidate molecular biomarkers (AZGP1, BTG2, CBX2, CNTNAP2, DNALI1, LOC389033, NME5, PIP, S100A8, SCUBE2, SERPINA11, STC2, STK32B, SUSD3, UBE2C, and WHSC1L1) in invasive breast carcinoma [14, 15, 23]. To investigate whether these putative prognostic biomarkers may also play a pivotal role in the aggressive nature of OSCCs, the effect of altered gene expression patterns on clinical outcome was evaluated using two external OSCC patient cohorts (Cohorts I-II; Table 1). Cox proportional hazard models were calculated for each gene with relation to clinical endpoints (DSS and/or OS; Table 2).
In Cohort I, two genes (LOC389033 and SERPINA11) were not found on the Affymetrix platform and therefore excluded from the analysis. Univariate Cox regression analysis showed that low levels of AZGP1 (P = 0.001), BTG2 (P = 0.020), PIP (P = 0.010), and SCUBE2 (P = 0.026) were indicative of a more unfavorable prognosis, whereas elevated levels of STC2 (P = 0.001) and UBE2C (P = 0.042) had an adverse effect on DSS. For OS, low levels of BTG2 (P = 0.005) and elevated levels of CBX2 (P = 0.005), STC2 (P = 0.001), and STK32B (P = 0.015) were predictive of outcome. In addition, low S100A8 mRNA levels (P = 0.053) were borderline significant for OS.
For Cohort II, the 16-gene signature was evaluated in RNASeq expression profiling data for 198 OSCC patients. Univariate Cox regression analysis showed that elevated levels of UBE2C mRNA levels (P = 0.005) were indicative of OS. On the other hand, elevated levels of NME5 (P = 0.050), low S100A8 levels (P = 0.072), and low SCUBE2 levels (P = 0.067) were borderline significant for OS. Furthermore, low S100A8 levels (P < 0.001; log2ratio = -2.61) and elevated UBE2C levels (P = 0.007; log2ratio = 0.713) were significantly associated with high histological grade.
Protein expression levels of the molecular biomarkers in OSCC specimens
Protein expression levels for the candidate biomarkers were evaluated using immunohistochemistry with 45 full-face FFPE specimens representing 43 OSCC patients (Cohort III; Table 1). PIP and DNALI1 were excluded from further analysis due to low expression levels in OSCC samples, whereas LOC389033 was excluded because the gene is not expressed at the protein level. Immunopositivity was shown for all of the examined proteins in peritumoral normal mucous membrane, the salivary glands, and dysplasia, with the exception of AZGP1, SUSD3, UBE2C, and WHSC1L1. AZGP1 and UBE2C were strongly positive in the basal cell layer; SUSD3 was positive in the mucous membrane but negative in the salivary glands, whereas WHSC1L1-positivity was shown in the layers of muscle tissue. In addition, no S100A8 staining was observed in the basal cell layer. In invasive tissue, immunopositivity for the 14 analyzed proteins ranged from 14-86% with CNTNAP2 and WHSC1L1 having the lowest and highest incidence rates, respectively (Table 3). Interestingly, there was only one reported case (6%) of SCUBE2-positivity in lymph node-positive tumors, compared with SCUBE2-positivity in 41% of lymph node-negative tumors. In addition, p16 immunopositivity was observed in 8/43 tumor specimens (19%).
Correlation of the molecular biomarkers with clinicopathological features
To investigate whether heterogeneous protein expression of the analyzed antigens is clinically relevant, a correlation analysis was performed with established clinicopathological features (Additional file 2: Table S2). S100A8 was strongly associated with tumor differentiation (P = 0.009), e.g. tumors with enhanced S100A8 expression levels were frequently well differentiated (64%) compared with 17% in S100A8-negative tumors. In addition, SCUBE2 was significantly associated with lymph node status (P = 0.01), CBX2 with tumor inflammatory infiltration (P = 0.03), and STK32B with tumor size (P = 0.04). Interestingly, a high proportion of SCUBE2-positive tumors (11/12) were lymph node-negative and all STK32B-positive tumors (15/15) were smaller in size (T1-T2 tumors); minimal peritumoral inflammatory infiltration was found in tumors with reduced CBX2 levels. Additionally, we also found a slight indication that UBE2C and SCUBE2, SERPINA11, and NME5 were associated with tumor differentiation (P = 0.06 and P = 0.09, respectively), inflammatory infiltration (P = 0.09), and p16 expression (P = 0.08), respectively.
Prognostic significance of the molecular biomarkers
Next, we examined the prognostic significance of the proposed biomarkers using disease-specific survival and overall survival. OSCC patients with tumors displaying enhanced CNTNAP2 levels had significantly shorter DSS (P = 0.010; HR (95% CI) = 5.70 (1.27-25.57)), whereas patients with S100A8-negative tumors had significantly shorter OS (P = 0.0063; HR (95% CI) = 0.10 (0.014-0.76); Figure 1). Our data suggest a slight association between SCUBE2 expression and DSS (P = 0.090), as well as UBE2C expression and OS (P = 0.074). CNTNAP2 expression had no significant effect on DSS after adjusting for tumor size and lymph node status (P = 0.10 - HR (95% CI) = 3.56 (0.78-16.17)). Furthermore, outcome prediction was not improved using a predictive model for DSS including CNTNAP2 expression, lymph node status, and tumor size (C-index = 0.949) compared with a model containing lymph node status and tumor size (C-index = 0.941; Figure 2). Following multivariate analysis adjusting for tumor size, lymph node status, differentiation and age, S100A8 was still statistically significant (P = 0.013 - HR (95% CI) = 0.11 (0.013-0.92; Table 4)). Combining S100A8 in a predictive model for OS with tumor size, lymph node status, differentiation and age improved outcome prediction significantly from 0.605 to 0.833 (Figure 2).
Oral squamous cell carcinoma is a heterogeneous disease with diverse clinical, pathological, and biological behavior . Nevertheless, the strongest determinants of prognosis still include tumor stage and the presence of cervical metastases at the time of diagnosis, as well as the time to locoregional recurrences [9, 24–26]. Unfortunately, up to 50% of OSCCs are diagnosed at an advanced stage with 5-year survival rates at approximately 60%, e.g. delayed diagnosis [27–30]. Therefore, many patients could benefit greatly from complementary molecular markers, which may help guide treatment decisions and be of value in the development of new therapeutic agents. Extensive efforts are currently being made to identify and validate biomarkers based on the biology of oral cancers that can complement established clinicopathological features and improve clinical management of the disease. Recent work to characterize OSCC using transcriptomic profiling has mainly focused on the identification of biomarkers for disease progression and lymph node metastasis prediction [31–38]. Surprisingly, few gene expression signatures have been developed to improve patient risk assessment [18, 34, 39, 40]. Although transcriptome analyses can give an indication of biological activity within a tissue, mRNA and protein levels may not correlate because gene expression is controlled by a multistage system. We have therefore performed immunohistochemistry using readily available FFPE samples to evaluate the clinical significance of the 16 putative biomarkers in OSCC. Here, we show the prognostic potential of the S100A8 and CNTNAP2 proteins, as well as the relationship between SCUBE2, CBX2, and STK32B protein levels and important clinicopathological features for OSCC, i.e. regional metastasis to the cervical lymph nodes, tumor inflammatory infiltration, and tumor size, respectively.
In the present work, we evaluated the applicability of breast cancer prognostic biomarkers for OSCC, given the biological similiarities between the two cancer types . In breast carcinoma, we showed the recurrent upregulation of the CBX2, CNTNAP2, S100A8, UBE2C, and WHSC1L1 genes as well as downregulation of the AZGP1, BTG2, DNALI1, LOC389033, NME5, PIP, SCUBE2, SERPINA11, STC2, STK32B, and SUSD3 genes in more aggressive tumors [14, 15, 23]. The prognostic potential of individual biomarkers in the 16-gene signature were evaluated in three OSCC patient cohorts containing clinical information, two of which (Cohorts I-II) were compiled from publicly available Affymetrix and Illumina RNASeq datasets. The Affymetrix one-channel system and Illumina RNAseq platform have recently been shown to correlate well, in particular for high abundance genes . Although the majority of the markers in the signature were in agreement for both breast carcinoma and OSCC, five proteins (S100A8, STC2, STK32B, SUSD3, and WHSC1L1) were inversely regulated in aggressive OSCC samples. In addition, several of the biomarkers which showed promise at the mRNA level (CBX2, S100A8, SCUBE2, and STK32B) were also either predictive of clinical outcome or associated with clinicopathological features at the protein level. The differences in the prognostic potential of specific biomarkers at the mRNA or protein levels may possibly be the result of small sample sizes, discordant mRNA-protein expression patterns or differences in sample preparation for microarray analyses and immunohistochemistry, i.e. microarray is frequently performed using a tumor mass, which consists of both malignant and nonmalignant cells, resulting in an over- or underestimation of expression levels, whereas cell type-specific protein expression patterns can be easily interpreted using immunohistochemistry. In breast carcinoma, elevated levels of the S100A8 protein were associated with moderate/strong tumor inflammatory infiltration and significantly shorter DSS rates. In contrast, lower S100A8 protein levels were associated with significantly shorter OS rates in OSCC, but there was no association with tumor infiltration. These findings suggest differences in cell type-specific gene expression and thereby activation and/or inhibition of diverse cellular mechanisms. Furthermore, CBX2, CNTNAP2, S100A8, SCUBE2, and STK32B protein levels were significantly associated with established clinicopathological features for OSCC and/or patient clinical outcome. We propose here two potential prognostic biomarkers for OSCC (the CNTNAP2 protein for DSS and the S100A8 protein for OS). There may be several reasons why these candidate biomarkers were associated with different clinical endpoints: differences in gene function and whether the biomarkers correlated with strong determinants of OSCC-related prognosis, e.g. tumor stage, lymph node metastasis or time to locoregional recurrences.
Chen et al. identified upregulation of CNTNAP2 mRNA levels as a distinctive characteristic of OSCC in comparison with normal mucosa . In a genome-wide association study for oral cancers, CNTNAP2 was found to be associated with cell migration . The CNTNAP2 gene also functions as a cell adhesion molecule and has been found to be either methylated or deleted in several different cancer types, e.g. glioma, myleoid leukemia, and pancreatic adenocarcinoma [43–45]. In the present study, elevated CNTNAP2 protein levels were prevalent in tumors from pN1 patients, but this association was not statistically significant likely due to the low number of events. These findings and the known function of the CNTNAP2 protein in tumorigenesis suggest that it may be a determinant of metastatic spread to the cervical lymph nodes and thereby a more aggressive phenotype with an adverse effect on DSS. Although elevated CNTNAP2 levels were predictive of a more unfavorable prognosis in univariate Cox regression models, disease-specific survival rates were not significantly improved by including the protein in a model containing lymph node status and tumor size. These results show the profound effect node metastases have on clinical outcome. In future studies, it will be necessary to alter the study design to test the prognostic potential of CNTNAP2 further, e.g. using a larger cohort containing node-negative patients. Furthermore, the S100A8 gene, also known as MRP-8 or calgranulin A, has been studied extensively in head and neck cancer. Inclusion of S100A8 protein expression patterns in multivariate models together with tumor size, node status, and age significantly improved OS prediction. In addition, we found an association between low S100A8 levels (mRNA levels in Cohort II and protein levels in Cohort III) and tumor differentiation. S100A8 belongs to the S100 gene family, several of which form a gene cluster on chromosome 1q21 and are commonly induced by chronic and acute inflammation [46–48]. In the present study, positive S100A8 immunostaining was also observed in a peritumoral inflammatory cell population, albeit similar in both S100A8-positive and negative tumors. Overexpression of the S100A8 protein has been observed in breast, colorectal, gastric, lung, pancreatic, and prostate cancer, wheras underexpression has been shown in various squamous cell carcinomas of the head and neck . Recently, the S100A8/A9 heterocomplex was identified as a regulator of cell cycle progression and cell proliferation in cancer cell lines originating from the head and neck region .
OSCC outcome prediction may also be improved by identifying biomarkers that can predict which tumors will inevitabily spread to the cervical nodes. When we initiated this study, we believed that the CNTNAP2 gene would be a strong candidate for lymph node status because of its association with metastatic spread to the axillary lymph nodes in breast carcinoma . However in the present study, loss of SCUBE2 protein expression in the tumor epithelial component was the only marker significantly associated with lymph node status and borderline significant for DSS. Interestingly, all 12 patients with SCUBE2-positive tumors were long-term disease-specific survivors and only 1/12 tumors had metastasized to the cervical lymph nodes.These findings are consistent with reports in breast carcinoma where elevated levels of the SCUBE2 protein not only inhibited proliferation of breast cancer cell lines, but was indicative of a more favorable prognosis [23, 51]. SCUBE2 is a multidomain, secreted glycoprotein that was first implicated in functioning upstream of hedgehog signaling in zebrafish, but has since been shown to also play a pivotal role in human tumorigenesis [52, 53]. Furthermore, SCUBE2 serves two different functions by using either the CUB domain at the C-terminal end or the EGF-like repeats at the N-terminus, i.e. control of cell proliferation by antagonizing bone morphogenetic protein (BMP) activity and regulation of cell-cell contacts, and possibly tumor metastasis, by forming complexes with E-cadherin in adherens junctions, respectively .
The role CBX2 and STK32B play in OSCC tumorigenesis is largely unknown. In a correlation analysis between protein expression patterns and established clinicopathological features, a clear link was found between CBX2-positivity and tumor inflammatory infiltration as well as STK32B-positivity and small tumor size (T1-T2). CBX2 plays a role in epigenetic regulation and hematopoietic stem cell differentiation, whereas the STK32B gene is one of several serine/threonine kinases (STK32B, STK36, and STK39) with similar gene expression patterns in basal-like breast cancers [16, 55, 56]. Consequently, STK36 is a component of the hedgehog signaling pathway. The function of the STK32B gene is not known, but single nucleotide polymorphisms in the STK32B gene have been associated with oral clefts .
Hedgehog signaling has also been implicated in OSCC and basal cell carcinoma, a skin cancer frequently found in the head and neck region [58, 59]. Leovic and colleagues illustrated the clinical relevance and aberrant expression of key hedgehog signaling components in SCCs of the oral cavity and oropharynx, i.e. PTCH1, GLI1, SMO, and SHH. In Cohort II, several genes which play a role in hedgehog signaling were also found to be associated with high histological grade (the KIF7, TCTN1, and TCTN2 genes) and metastatic spread to the cervical lymph nodes (the FGF4 and PTCHD1 genes). These findings further emphasize the importance of hedgehog signaling in OSCC.
In summary, our results illustrate that cancer is a complex, heterogeneous disease. However, certain clinically relevant biomarkers may be useful therapeutic targets for several different cancer types. Despite the limited number of patients included in the study, the integrative analysis identified at least five of the putative biomarkers originally developed in breast carcinoma (CBX2, CNTNAP2, S100A8, SCUBE2, and STK32B) that also improve risk assessment for OSCC patients and may play a crucial role in cancer-related processes, e.g. hedgehog signaling, regulation of cell cycle progression, cell proliferation, metastasis. These targets should be studied in a larger cohort to further evaluate their clinical significance in OSCC.
squamous cell carcinoma of the oral cavity
cervical lymph node metastases
- H & E:
Hematoxylin and eosin
Formalin-fixed, paraffin-embedded tumor tissues
- AUC (t):
time-dependent area under the ROC curve.
Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D: Global cancer statistics. CA Cancer J Clin. 2011, 61 (2): 69-90. 10.3322/caac.20107.
Lozano R, Naghavi M, Foreman K, Lim S, Shibuya K, Aboyans V, Abraham J, Adair T, Aggarwal R, Ahn SY, Alvarado M, Anderson HR, Anderson LM, Andrews KG, Atkinson C, Baddour LM, Barker-Collo S, Bartels DH, Bell ML, Benjamin EJ, Bennett D, Bhalla K, Bikbov B, Bin Abdulhak A, Birbeck G, Blyth F, Bolliger I, Boufous S, Bucello C, Burch M, et al: Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012, 380 (9859): 2095-2128. 10.1016/S0140-6736(12)61728-0.
Dahlstrand H, Dahlgren L, Lindquist D, Munck-Wikland E, Dalianis T: Presence of human papillomavirus in tonsillar cancer is a favourable prognostic factor for clinical outcome. Anticancer Res. 2004, 24 (3b): 1829-1835.
Hammarstedt L, Dahlstrand H, Lindquist D, Onelov L, Ryott M, Luo J, Dalianis T, Ye W, Munck-Wikland E: The incidence of tonsillar cancer in Sweden is increasing. Acta Otolaryngol. 2007, 127 (9): 988-992. 10.1080/00016480601110170.
Hammarstedt L, Lindquist D, Dahlstrand H, Romanitan M, Dahlgren LO, Joneberg J, Creson N, Lindholm J, Ye W, Dalianis T, Munck-Wikland E: Human papillomavirus as a risk factor for the increase in incidence of tonsillar cancer. Int J Cancer. 2006, 119 (11): 2620-2623. 10.1002/ijc.22177.
Hashibe M, Brennan P, Chuang SC, Boccia S, Castellsague X, Chen C, Curado MP, Dal Maso L, Daudt AW, Fabianova E, Fernandez L, Wunsch-Filho V, Franceschi S, Hayes RB, Herrero R, Kelsey K, Koifman S, La Vecchia C, Lazarus P, Levi F, Lence JJ, Mates D, Matos E, Menezes A, McClean MD, Muscat J, Eluf-Neto J, Olshan AF, Purdue M, Rudnai P, et al: Interaction between tobacco and alcohol use and the risk of head and neck cancer: pooled analysis in the International Head and Neck Cancer Epidemiology Consortium. Cancer Epidemiol Biomarkers Prev. 2009, 18 (2): 541-550. 10.1158/1055-9965.EPI-08-0347.
Oguejiofor KK, Hall JS, Mani N, Douglas C, Slevin NJ, Homer J, Hall G, West CM: The Prognostic Significance of the Biomarker p16 in Oropharyngeal Squamous Cell Carcinoma. Clin Oncol (R Coll Radiol). 2013, 25 (11): 630-638. 10.1016/j.clon.2013.07.003.
Syrjanen K, Syrjanen S, Lamberg M, Pyrhonen S, Nuutinen J: Morphological and immunohistochemical evidence suggesting human papillomavirus (HPV) involvement in oral squamous cell carcinogenesis. Int J Oral Surg. 1983, 12 (6): 418-424. 10.1016/S0300-9785(83)80033-7.
Kernohan MD, Clark JR, Gao K, Ebrahimi A, Milross CG: Predicting the prognosis of oral squamous cell carcinoma after first recurrence. Arch Otolaryngol Head Neck Surg. 2010, 136 (12): 1235-1239. 10.1001/archoto.2010.214.
Chung CH, Parker JS, Karaca G, Wu J, Funkhouser WK, Moore D, Butterfoss D, Xiang D, Zanation A, Yin X, Shockley WW, Weissler MC, Dressler LG, Shores CG, Yarbrough WG, Perou CM: 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.
Cueni LN, Hegyi I, Shin JW, Albinger-Hegyi A, Gruber S, Kunstfeld R, Moch H, Detmar M: Tumor lymphangiogenesis and metastasis to lymph nodes induced by cancer cell expression of podoplanin. Am J Pathol. 2010, 177 (2): 1004-1016. 10.2353/ajpath.2010.090703.
Wu CH, Tu HF, Gong NR, Chen YW, Liu CJ, Kao SY: The Val allele of HER-2 codon 655 predicts the progression of oral squamous cell carcinoma. Oral Oncol. 2009, 45 (7): 579-583. 10.1016/j.oraloncology.2008.07.018.
Futreal PA, Coin L, Marshall M, Down T, Hubbard T, Wooster R, Rahman N, Stratton MR: A census of human cancer genes. Nat Rev Cancer. 2004, 4 (3): 177-183. 10.1038/nrc1299.
Mollerstrom E, Kovacs A, Lovgren K, Nemes S, Delle U, Danielsson A, Parris T, Brennan DJ, Jirstrom K, Karlsson P, Helou K: Up-regulation of cell cycle arrest protein BTG2 correlates with increased overall survival in breast cancer, as detected by immunohistochemistry using tissue microarray. BMC Cancer. 2010, 10: 296-10.1186/1471-2407-10-296.
Parris TZ, Kovacs A, Aziz L, Hajizadeh S, Nemes S, Semaan M, Forssell-Aronsson E, Karlsson P, Helou K: Additive effect of the AZGP1, PIP, S100A8 and UBE2C molecular biomarkers improves outcome prediction in breast carcinoma. Int J Cancer. 2014, 134 (7): 1617-1629. 10.1002/ijc.28497.
Finetti P, Cervera N, Charafe-Jauffret E, Chabannon C, Charpin C, Chaffanet M, Jacquemier J, Viens P, Birnbaum D, Bertucci F: Sixteen-kinase gene expression identifies luminal breast cancers with poor prognosis. Cancer Res. 2008, 68 (3): 767-776. 10.1158/0008-5472.CAN-07-5516.
McKiernan E, McDermott EW, Evoy D, Crown J, Duffy MJ: The role of S100 genes in breast cancer progression. Tumour Biol: j Int Soc Oncodev Biol Med. 2011, 32 (3): 441-450.
Lohavanichbutr P, Mendez E, Holsinger FC, Rue TC, Zhang Y, Houck J, Upton MP, Futran N, Schwartz SM, Wang P, Chen C: A 13-gene signature prognostic of HPV-negative OSCC: discovery and external validation. Clin Cancer Res. 2013, 19 (5): 1197-1203. 10.1158/1078-0432.CCR-12-2647.
World Health Organization Classification of Tumours. Pathology and Genetics of Head and Neck Tumours. Edited by: Barnes LEJW, Reichart P, Sidransky D. 2005, Lyon: IARC Press
The TNM Classification of Malignant Tumours. Edited by: Sobin LH, Gospodarowicz MK, Wittekind C. 2009, Oxford: Wiley-Blackwell, 7
McCarty KS, Miller LS, Cox EB, Konrath J, McCarty KS: Estrogen receptor analyses. Correlation of biochemical and immunohistochemical methods using monoclonal antireceptor antibodies. Arch Pathol Lab Med. 1985, 109 (8): 716-721.
Camp RL, Dolled-Filhart M, Rimm DL: X-tile: a new bio-informatics tool for biomarker assessment and outcome-based cut-point optimization. Clin Cancer Res. 2004, 10 (21): 7252-7259. 10.1158/1078-0432.CCR-04-0713.
Parris TZ, Danielsson A, Nemes S, Kovacs A, Delle U, Fallenius G, Mollerstrom E, Karlsson P, Helou K: Clinical implications of gene dosage and gene expression patterns in diploid breast carcinoma. Clin Cancer Res. 2010, 16 (15): 3860-3874. 10.1158/1078-0432.CCR-10-0889.
De Araujo RF, Jr BCA, Clebis NK, De Moura SA, Lopes Costa Ade L: Prognostic significance of the anatomical location and TNM clinical classification in oral squamous cell carcinoma. Med Oral Patol Oral Cir Bucal. 2008, 13 (6): E344-E347.
Grandi C, Alloisio M, Moglia D, Podrecca S, Sala L, Salvatori P, Molinari R: Prognostic significance of lymphatic spread in head and neck carcinomas: therapeutic implications. Head Neck Surg. 1985, 8 (2): 67-73. 10.1002/hed.2890080202.
Kademani D, Bell RB, Bagheri S, Holmgren E, Dierks E, Potter B, Homer L: Prognostic factors in intraoral squamous cell carcinoma: the influence of histologic grade. J Oral Maxillofac Surg: official j Am Asso Oral and Maxillofacial Surg. 2005, 63 (11): 1599-1605. 10.1016/j.joms.2005.07.011.
Brouha XD, Tromp DM, Hordijk GJ, Winnubst JA, de Leeuw JR: Oral and pharyngeal cancer: analysis of patient delay at different tumor stages. Head Neck. 2005, 27 (11): 939-945. 10.1002/hed.20270.
Howlader NNA, Krapcho M, Neyman N, Aminou R, Waldron W, Altekruse SF, Kosary CL, Ruhl J, Tatalovich Z, Cho H, Mariotto A, Eisner MP, Lewis DR, Chen HS, Feuer EJ, Cronin KA, Edwards BK: SEER Cancer Statistics Review, 1975-2008. 2011, Bethesda, MD: National Cancer Institute
Onizawa K, Nishihara K, Yamagata K, Yusa H, Yanagawa T, Yoshida H: Factors associated with diagnostic delay of oral squamous cell carcinoma. Oral Oncol. 2003, 39 (8): 781-788. 10.1016/S1368-8375(03)00075-7.
Seoane-Romero JM, Vazquez-Mahia I, Seoane J, Varela-Centelles P, Tomas I, Lopez-Cedrun JL: Factors related to late stage diagnosis of oral squamous cell carcinoma. Med Oral Patol Oral Cir Bucal. 2012, 17 (1): e35-e40.
Brinkmann O, Kastratovic DA, Dimitrijevic MV, Konstantinovic VS, Jelovac DB, Antic J, Nesic VS, Markovic SZ, Martinovic ZR, Akin D, Spielmann N, Zhou H, Wong DT: Oral squamous cell carcinoma detection by salivary biomarkers in a Serbian population. Oral Oncol. 2011, 47 (1): 51-55. 10.1016/j.oraloncology.2010.10.009.
Chen C, Mendez E, Houck J, Fan W, Lohavanichbutr P, Doody D, Yueh B, Futran ND, Upton M, Farwell DG, Schwartz SM, Zhao LP: Gene expression profiling identifies genes predictive of oral squamous cell carcinoma. Cancer Epidemiol Biomarkers Prev. 2008, 17 (8): 2152-2162. 10.1158/1055-9965.EPI-07-2893.
Ginos MA, Page GP, Michalowicz BS, Patel KJ, Volker SE, Pambuccian SE, Ondrey FG, Adams GL, Gaffney PM: Identification of a gene expression signature associated with recurrent disease in squamous cell carcinoma of the head and neck. Cancer Res. 2004, 64 (1): 55-63. 10.1158/0008-5472.CAN-03-2144.
Lohavanichbutr P, Houck J, Doody DR, Wang P, Mendez E, Futran N, Upton MP, Holsinger FC, Schwartz SM, Chen C: Gene expression in uninvolved oral mucosa of OSCC patients facilitates identification of markers predictive of OSCC outcomes. PLoS One. 2012, 7 (9): e46575-10.1371/journal.pone.0046575.
Nguyen ST, Hasegawa S, Tsuda H, Tomioka H, Ushijima M, Noda M, Omura K, Miki Y: Identification of a predictive gene expression signature of cervical lymph node metastasis in oral squamous cell carcinoma. Cancer Sci. 2007, 98 (5): 740-746. 10.1111/j.1349-7006.2007.00454.x.
O'Donnell RK, Kupferman M, Wei SJ, Singhal S, Weber R, O'Malley B, Cheng Y, Putt M, Feldman M, Ziober B, Muschel RJ: Gene expression signature predicts lymphatic metastasis in squamous cell carcinoma of the oral cavity. Oncogene. 2005, 24 (7): 1244-1251. 10.1038/sj.onc.1208285.
Reis PP, Waldron L, Perez-Ordonez B, Pintilie M, Galloni NN, Xuan Y, Cervigne NK, Warner GC, Makitie AA, Simpson C, Goldstein D, Brown D, Gilbert R, Gullane P, Irish J, Jurisica I, Kamel-Reid S: A gene signature in histologically normal surgical margins is predictive of oral carcinoma recurrence. BMC Cancer. 2011, 11: 437-10.1186/1471-2407-11-437.
Ye H, Yu T, Temam S, Ziober BL, Wang J, Schwartz JL, Mao L, Wong DT, Zhou X: Transcriptomic dissection of tongue squamous cell carcinoma. BMC Genomics. 2008, 9: 69-10.1186/1471-2164-9-69.
Mendez E, Houck JR, Doody DR, Fan W, Lohavanichbutr P, Rue TC, Yueh B, Futran ND, Upton MP, Farwell DG, Heagerty PJ, Zhao LP, Schwartz SM, Chu C: A genetic expression profile associated with oral cancer identifies a group of patients at high risk of poor survival. Clin Cancer Res. 2009, 15 (4): 1353-1361. 10.1158/1078-0432.CCR-08-1816.
Pramana J, Van den Brekel MW, van Velthuysen ML, Wessels LF, Nuyten DS, Hofland I, Atsma D, Pimentel N, Hoebers FJ, Rasch CR, Begg AC: Gene expression profiling to predict outcome after chemoradiation in head and neck cancer. Int J Radiat Oncol Biol Phys. 2007, 69 (5): 1544-1552. 10.1016/j.ijrobp.2007.08.032.
Guo Y, Sheng Q, Li J, Ye F, Samuels DC, Shyr Y: Large scale comparison of gene expression levels by microarrays and RNAseq using TCGA data. PLoS One. 2013, 8 (8): e71462-10.1371/journal.pone.0071462.
Bhatnagar R, Dabholkar J, Saranath D: Genome-wide disease association study in chewing tobacco associated oral cancers. Oral Oncol. 2012, 48 (9): 831-835. 10.1016/j.oraloncology.2012.03.007.
Bralten LB, Gravendeel AM, Kloosterhof NK, Sacchetti A, Vrijenhoek T, Veltman JA, van den Bent MJ, Kros JM, Hoogenraad CC, Sillevis Smitt PA, French PJ: The CASPR2 cell adhesion molecule functions as a tumor suppressor gene in glioma. Oncogene. 2010, 29 (46): 6138-6148. 10.1038/onc.2010.342.
De Weer A, Poppe B, Vergult S, Van Vlierberghe P, Petrick M, De Bock R, Benoit Y, Noens L, De Paepe A, Van Roy N, Menten B, Speleman F: Identification of two critically deleted regions within chromosome segment 7q35-q36 in EVI1 deregulated myeloid leukemia cell lines. PLoS One. 2010, 5 (1): e8676-10.1371/journal.pone.0008676.
Omura N, Li CP, Li A, Hong SM, Walter K, Jimeno A, Hidalgo M, Goggins M: Genome-wide profiling of methylated promoters in pancreatic adenocarcinoma. Cancer Biol Ther. 2008, 7 (7): 1146-1156. 10.4161/cbt.7.7.6208.
Cross SS, Hamdy FC, Deloulme JC, Rehman I: Expression of S100 proteins in normal human tissues and common cancers using tissue microarrays: S100A6, S100A8, S100A9 and S100A11 are all overexpressed in common cancers. Histopathology. 2005, 46 (3): 256-269. 10.1111/j.1365-2559.2005.02097.x.
Odink K, Cerletti N, Bruggen J, Clerc RG, Tarcsay L, Zwadlo G, Gerhards G, Schlegel R, Sorg C: Two calcium-binding proteins in infiltrate macrophages of rheumatoid arthritis. Nature. 1987, 330 (6143): 80-82. 10.1038/330080a0.
Thorey IS, Roth J, Regenbogen J, Halle JP, Bittner M, Vogl T, Kaesler S, Bugnon P, Reitmaier B, Durka S, Graf A, Wockner M, Rieger N, Konstantinow A, Wolf E, Goppelt A, Werner S: The Ca2 + -binding proteins S100A8 and S100A9 are encoded by novel injury-regulated genes. J Biol Chem. 2001, 276 (38): 35818-35825. 10.1074/jbc.M104871200.
Gebhardt C, Nemeth J, Angel P, Hess J: S100A8 and S100A9 in inflammation and cancer. Biochem Pharmacol. 2006, 72 (11): 1622-1631. 10.1016/j.bcp.2006.05.017.
Khammanivong A, Wang C, Sorenson BS, Ross KF, Herzberg MC: S100A8/A9 (calprotectin) negatively regulates G2/M cell cycle progression and growth of squamous cell carcinoma. PLoS One. 2013, 8 (7): e69395-10.1371/journal.pone.0069395.
Cheng CJ, Lin YC, Tsai MT, Chen CS, Hsieh MC, Chen CL, Yang RB: SCUBE2 suppresses breast tumor cell proliferation and confers a favorable prognosis in invasive breast cancer. Cancer Res. 2009, 69 (8): 3634-3641. 10.1158/0008-5472.CAN-08-3615.
Kawakami A, Nojima Y, Toyoda A, Takahoko M, Satoh M, Tanaka H, Wada H, Masai I, Terasaki H, Sakaki Y, Takeda H, Okamoto H: The zebrafish-secreted matrix protein you/scube2 is implicated in long-range regulation of hedgehog signaling. Current biology: CB. 2005, 15 (5): 480-488. 10.1016/j.cub.2005.02.018.
Tsai MT, Cheng CJ, Lin YC, Chen CC, Wu AR, Wu MT, Hsu CC, Yang RB: Isolation and characterization of a secreted, cell-surface glycoprotein SCUBE2 from humans. Biochem J. 2009, 422 (1): 119-128. 10.1042/BJ20090341.
Lin YC, Chen CC, Cheng CJ, Yang RB: Domain and functional analysis of a novel breast tumor suppressor protein, SCUBE2. J Biol Chem. 2011, 286 (30): 27039-27047. 10.1074/jbc.M111.244418.
Hatano A, Matsumoto M, Higashinakagawa T, Nakayama KI: Phosphorylation of the chromodomain changes the binding specificity of Cbx2 for methylated histone H3. Biochem Biophys Res Commun. 2010, 397 (1): 93-99. 10.1016/j.bbrc.2010.05.074.
Klauke K, Radulovic V, Broekhuis M, Weersing E, Zwart E, Olthof S, Ritsema M, Bruggeman S, Wu X, Helin K, Bystrykh L, de Haan G: Polycomb Cbx family members mediate the balance between haematopoietic stem cell self-renewal and differentiation. Nat Cell Biol. 2013, 15 (4): 353-362. 10.1038/ncb2701.
Ingersoll RG, Hetmanski J, Park JW, Fallin MD, McIntosh I, Wu-Chou YH, Chen PK, Yeow V, Chong SS, Cheah F, Sull JW, Jee SH, Wang H, Wu T, Murray T, Huang S, Ye X, Jabs EW, Redett R, Raymond G, Scott AF, Beaty TH: Association between genes on chromosome 4p16 and non-syndromic oral clefts in four populations. Eur j hum genet : EJHG. 2010, 18 (6): 726-732. 10.1038/ejhg.2009.228.
Leovic D, Sabol M, Ozretic P, Musani V, Car D, Marjanovic K, Zubcic V, Sabol I, Sikora M, Grce M, Glavas-Obrovac L, Levanat S: Hh-Gli signaling pathway activity in oral and oropharyngeal squamous cell carcinoma. Head Neck. 2012, 34 (1): 104-112. 10.1002/hed.21696.
Wong CS, Strange RC, Lear JT: Basal cell carcinoma. BMJ (Clinical res). 2003, 327 (7418): 794-798. 10.1136/bmj.327.7418.794.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/14/324/prepub
The authors would like to thank Ulric Pedersen for his technical assistance. This work was supported by grants from the Swedish Cancer Society (KH), the King Gustav V Jubilee Clinic Cancer Research Foundation (KH), Assar Gabrielsson Research Foundation (LA), and the Wilhelm and Martina Lundgren Research Foundation (TZP and LA).
The authors declare that they have no competing interests.
TZP, LA, AK, CYC, PK, KH contributed to the study design. TZP, SH, MS, AK performed the immunohistochemical analyses. TZP performed the data and statistical analyses, and drafted the manuscript. LA collected clinical data, performed the evaluation of lymph nodes, data analyses, and drafted the manuscript. AK performed the pathologic review of the FFPE samples. SN performed the statistical analyses. AK, SH, SN, MS, CYC, PK, KH critically reviewed the manuscript. All authors have read and approved the final manuscript.
Toshima Z Parris, Luaay Aziz contributed equally to this work.
Electronic supplementary material
Additional file 2: Table S2: Relationship between protein expression and clinicopathological features in OSCC (Cohort III). (XLSX 28 KB)
About this article
Cite this article
Parris, T.Z., Aziz, L., Kovács, A. et al. Clinical relevance of breast cancer-related genes as potential biomarkers for oral squamous cell carcinoma. BMC Cancer 14, 324 (2014). https://doi.org/10.1186/1471-2407-14-324
- Oral squamous cell carcinoma
- Outcome prediction
- Molecular biomarker
- Model validation