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

Clinical correlation of opposing molecular signatures in head and neck squamous cell carcinoma



The concept of head and neck cancers (HNSCC) having unique molecular signatures is well accepted but relating this to clinical presentation and disease behaviour is essential for patient benefit. Currently the clinical significance of HNSCC molecular subtypes is uncertain therefore personalisation of HNSCC treatment is not yet possible.


We performed meta-analysis on 8 microarray studies and identified six significantly up- (PLAU, FN1, CDCA5) and down-regulated (CRNN, CLEC3B and DUOX1) genes which were subsequently quantified by RT-qPCR in 100 HNSCC patient margin and core tumour samples.


Retrospective correlation with sociodemographic and clinicopathological patient details identified two subgroups of opposing molecular signature (+q6 and -q6) that correlated to two recognised high-risk HNSCC populations in the UK. The +q6 group were older, male, and excessive alcohol users whilst the –q6 group were younger, female, paan-chewers and predominantly Bangladeshi. Additionally, all patients with tumour recurrence were in the latter subgroup.


We provide the first evidence linking distinct molecular signatures in HNSCC with clinical presentations. Prospective trials are required to determine the correlation between these distinct genotypes and disease progression or treatment response. This is an important step towards the ultimate goal of improving outcomes by utilising personalised molecular-signature-guided treatments for HNSCC patients.

Peer Review reports


Head and neck squamous cell carcinoma (HNSCC) is the 6th most common form of cancer worldwide. It is a multifactorial disease, with known risk factors including tobacco, alcohol, areca nut and human papilloma viruses (HPV). As with many cancers HNSCC occurs as a result of abnormal genetic alterations such as point mutations, amplifications, rearrangements and deletions of genes, paving the way for tumour progression [1]. However, no molecular testing technique has yet been developed to aid in early diagnosis and prognostic evaluation of HNSCCs. Whereas, other epithelial origin cancers such as breast and lung already have reliable diagnostic markers (mutant HER2 and EGFR, respectively [2, 3]) which are routinely used by oncologists to personalise treatments and improve outcomes for individual patients. Indeed in the UK, HNSCCs are one of the few cancers where incidence rates are still projected to rise in the future and mortality rates have not decreased despite the significant advances in oncological management (

Efforts are being made to find reliable HNSCC biomarkers that reflect the molecular make-up of the tumours. Through systemic reviews and meta-analysis, EGFR and cyclic D1 have been identified as potential serum diagnostic markers [4], and ANO1 and FADD reported as possible prognostic markers [5]. Genomic changes such as hypermethylation of RAS association domain family protein 1a (RASSF1A), a tumour suppressor gene, has been associated with a high risk of developing HNSCC [6]. However, currently none of these have translated into clinical application.

This study aims to explore the expression of potential HNSCC biomarkers for both diagnostic and prognostic purposes. Meta-analyses of eight independent HNSCC microarray studies was carried out to identify significantly up- and down-regulated genes in studies comparing HNSCC with normal oral mucosa [7,8,9,10,11,12,13,14]. The expression of a panel of likely genes was established in HNSCC patient samples and the molecular findings correlated with each patient’s clinical and histopathological features. We show, for the first time, that within HNSCC patients exist two sub-groups, which are molecularly and clinically distinct from each other. Knowledge of the existence of such heterogeneity will aid in developing personalised treatments to improve outcomes for HNSCC patients.


Clinical samples

The use of fresh clinical specimens collected in the UK was approved by the NHS Research Ethics Committee (06/MRE03/69). All tissue samples were previously collected according to local ethical committee-approved protocols and informed patient consent was obtained from all participants [15, 16]. Fresh tissue biopsies were preserved in RNALater (#AM7022, Ambion, Applied Biosystems, Warrington, UK) and stored short-term at 4 °C (1–7 days) prior to transportation and subsequent storage at − 20 °C until used. All frozen samples were digested with nuclease-free proteinase K at 60 °C prior to mRNA extraction (Dynabeads mRNA Direct kit, Invitrogen).

Transcriptome data mining

Transcriptome datasets were queried in the Oncomine ( database. The main inclusion criterion was that the studies must involve comparison between HNSCC tumour samples with normal tissues. Studies using HNSCC cell lines were excluded. At the time of analysis, there were 8 studies eligible for analysis (Fig. 1a). Differentially expressed genes were ranked according to their median P-values for over-expression and under-expression. Candidate genes were selected based on their top-ranking positions across the 8 studies resulting in a total of 20 upregulated and 20 downregulated genes were shortlisted (Fig. 1b) and were used for subsequent gene expression validation in our HNSCC cell line models using RT-qPCR.

Fig. 1
figure 1

Bioinformatics meta-analysis of 8 independent microarray studies on HNSCC tissues samples compared with normal oral tissues. a, Information for the 8 microarray studies: HNSCC anatomical sites, PubMed ID (PMID) referenced to published paper, microarray data archive (GEO or *Oncomine), number of patients’ tumour, normal and lymph-node metastastic (LNM) samples were as indicated. b, Based on statistical ranking of the most differentially expressed genes, top 10 upregulated and top 10 downregulated across the 8 studies were shortlisted for further validation on cell lines. c, Relative gene expression mRNA levels (Log2 Ratio) were measured using RT-qPCR and compared in a panel of 8 primary normal human oral keratinocytes (OK355, HOKG, OK113, NOK, NOK1, NOK3, NOK16 and NOK376) and 10 HNSCC cell lines (SCC4, SCC9, SCC15, SCC25, SqCC/Y1, UK1, VB6, CaLH2, CaDec12 and 5PT). We identified 3 most significantly upregulated (PLAU, FN1 and CDCA5) and 3 most significantly downregulated (CRNN, CLEC3B and DUOX1) in HNSCC cell lines and their gene putative functions (from NCBI’s Gene database) are listed in D

Reverse transcription quantitative PCR (RT-qPCR)

The RT-qPCR methodology was performed as described previously [15, 16]. Reverse transcription of purified mRNA were converted into cDNA using Transcriptor cDNA Synthesis kit (Roche Diagnostics Ltd., England, UK) and relative gene expression were performed using SYBR Green I Master (Roche) in the LightCycler 480 qPCR system (Roche) based on our published protocols [17,18,19] which are MIQE compliant [20]. Thermocycling begins with 95 °C for 5 min prior to 45 cycles of amplification at 95 °C for 10s, 60 °C for 6 s, 72 °C for 6 s, 76 °C for 1 s (data acquisition). A ‘touch-down’ annealing temperature intervention (66 °C starting temperature with a step-wise reduction of 0.6 °C/cycle; 8 cycles) was introduced prior to the amplification step to maximise primer specificity. Melting analysis (95 °C for 30s, 65 °C for 30s, 65–99 °C at a ramp rate of 0.11 °C/s with a continuous 5 acquisitions/°C) was performed at the end of qPCR amplification to validate single product amplification in each well. Relative quantification of mRNA transcripts was calculated based on an objective method using the second derivative maximum algorithm [21] (Roche). Sequences of the qPCR primers used in this study are provided in Table 1. All target genes were normalised to two stable reference genes (YAP1 and POLR2A) validated previously [17] to be amongst the most stable reference genes across a wide variety of primary human epithelial cells, dysplastic and squamous carcinoma cell lines, using the GeNorm algorithm [22]. No template controls (NTC) were prepared by omitting cells/tissue sample during RNA purification and eluates were used as NTCs for qPCR assays to monitor contamination.

Table 1 q6 Primer Sequences

Statistical analysis

Gene expression data were exported from Roche LightCycler LC480 Software as text files for subsequent analysis. Statistical analysis was carried out by the t-test on Graph Pad Prism software and Microsoft Excel and the Mann Whitney U test on SPSS software.


Microarray data mining and gene selection

The cancer microarray database Oncomine [23] ( was used to select 8 studies which analysed HNSCC cancer samples versus normal as shown in Fig. 1a. From this, top 20 differentially expressed genes were selected based on the reported P-value (> 0.001), out of which 10 were significantly upregulated and 10 significantly downregulated (Fig. 1b). Primers for each gene were designed using the Roche Applied Science Universal Probe Library Assay Design Centre for RT-qPCR assays. Initial testing of primer specificity and expression of selected biomarkers was established on cDNA from a panel of 8 normal primary oral keratinocytes and 10 HNSCC cell lines. Based on primer specificity and good reproducibility, three differentially expressed biomarkers were identified, of which, PLAU, FN1 and CDCA5 were found to be upregulated; whereas CRNN, CLEC3B and DUOX1 were downregulated when comparing HNSCC cell lines to normal oral epithelial cells in culture (Fig. 1c). For ease we called them q6 (quantification of selected six biomarkers). The q6 biomarkers were found to be involved in important cellular functions, listed in Fig. 1d.

Biomarkers identified molecularly distinct HNSCC subgroups

The selected candidate biomarkers were then validated on HNSCC patient tissue specimens by RT-qPCR which provided quantitative data on expression of these selected genes (relative to two reference genes YAP1 and POLR2A) on paired margin and tumour core tissue samples. Individual gene expressions were determined in each margin and core tumour tissue pairs. In order to obtain a clinically meaningful index value from the 6 genes in each patient, we derived an equation to summarise the degree of differential gene expression of the 6 genes between margin and core tissues in each patient:

$$ \mathrm{q}6\ \mathrm{Value}=\left(\mathrm{Sum}\ \mathrm{of}\ \mathrm{Log}2\ \mathrm{Ratios}\ \mathrm{of}\ \mathrm{the}\ 3\ \mathrm{upregulated}\ \mathrm{genes}\right)-\left(\mathrm{Sum}\ \mathrm{of}\ \mathrm{Log}2\ \mathrm{Ratios}\ \mathrm{of}\ \mathrm{the}\ 3\ \mathrm{downregulated}\ \mathrm{genes}\right) $$

We found that within the 100 patient samples analysed, two molecularly distinct populations could be identified (Fig. 2a). The majority (> 70%) of patients had positive q6 (+q6) values showing the predicted expression of q6 biomarkers from our HNSCC cell line data, with PLAU, FN1 and CDCA5 being upregulated and CRNN, CLEC3B and DUOX1 downregulated. An example of a patient with +q6 expression pattern is shown in Fig. 2b. On the opposite spectrum, about 20% of patients showed negative q6 (−q6) values indicating that these patients showed inverse expression of the q6 biomarkers whereby PLAU, FN1 and CDCA5 were downregulated and CRNN, CLEC3B and DUOX1 upregulated. An example of a patient with -q6 expression pattern is shown in Fig. 2c.

Fig. 2
figure 2

Validation of the 6 markers on 100 HNSCC patients with paired tumour core and margin tissue samples. a, The relative mRNA expression levels of each of the 6 markers were measured using RT-qPCR against two reference genes (YAP1 and POLR2A). Differential expression ratios (q6 values) were derived from Log2 Ratio of 3 upregulated markers (PLAU, FN1 and CDCA5) against 3 downregulated markers (CRNN, CLECB and DUOX1). Majority of HNSCC patients showed positive q6 values (indicating the expected expression pattern, as shown in panel b) whilst ~ 20% patients showed negative q6 values indicating an inversed expression patterns (as shown in panel c). d, Statistical analyses of sociodemographic and clinicopathological findings comparing +q6 and -q6 groups

Clinicopathological analysis of the two HNSCC (+q6 vs -q6) subgroups

In order to understand the different pattern of q6 biomarker expression, patient’s clinical reports were correlated with their molecular findings. Patient and tumour details were collated retrospectively for the 20 most positive and 20 most negative q6 values. In each group, data collection was incomplete for a number of patients as the clinical records were not traceable or the available clinical records were incomplete and therefore these individuals were removed from the analysis. Clinical data from 13 patients with +q6 values was compared to 8 patients with -q6 values. All patients were followed up for at least 2 years following surgical resection of their primary tumours and all tumours were histologically confirmed as HNSCC. The 18 primary HNSCC were excised and subsequently treated with post-operative radiotherapy (RT) with the tissue samples for this study being taken prior to RT whilst the 3 recurrences had originally been resected and all had post-op RT at that time.

We found statistically significant differences in age, sex, ethnicity, alcohol consumption and paan usage among the two groups. In the +q6 group the mean age was 63 with more males than females, while the -q6 group mean age was 56 with more females than males (P = 0.04). Additionally, more patients of Bangladeshi descent were found in the -q6 group (P = 0.04) (Fig. 2d).

Statistically significant difference was also found in the two groups with regards to associated risk factors. High levels of alcohol consumption in the +q6 group (P = 0.04), compared to the -q6 group who were often paan chewers (P = 0.02) (Fig. 2d). No difference was found in the smoking habits, tumour site and size among the two groups. Recurrent lesions were only found in the -q6 group although this was not statistically significant (Fig. 2d).

Prognostic values for q6 biomarkers on other cancer types

We further investigated the prognostic significance of the q6 markers on breast, ovarian, lung, gastric and liver cancers (Fig. 3) using the Kaplan Meier plotter transcriptome database [24] containing 54,675 genes on survival using 10,825 cancer samples (as of 16 Jan 2018). These include 5143 breast, 1816 ovarian, 2437 lung, 1065 gastric and 364 liver cancer patients with a mean follow-up of 69, 40, 49, 33 and 30 months, respectively ( We found that poor prognosis was associated with high expressions of PLAU1 and FN1 on ovarian, lung and gastric cancers. Similarly, high expression of CDCA5 were associated with poor prognosis of breast, lung and liver cancers. Low expression of CRNN, CLEC3B and DUOX1 were associated with poor prognosis of breast cancer. Downregulation of CLEC3B was also associated with poor prognosis in lung and liver cancers. Curiously, gastric cancer showed inverse relationship for CDCA5, CRNN and DUOX1 on prognosis compared to other cancer types (red asterisks; Fig. 3).

Fig. 3
figure 3

Prognostic significance analysis of the 6 markers on breast, ovarian, lung, gastric and liver cancers using the Kaplan Meier plotter transcriptome database containing 54,675 genes on survival using 10,825 cancer samples (as of 16 Jan 2018). *P < 0.05, **P < 0.01 and ***P < 0.001 showed expected prognostic patterns corresponding to each marker expression levels. Interestingly, * or *** (in red) showed inverse survival curve whereby high expression of e.g., CDCA5 was significantly associated with better prognosis in gastric cancer patients but poorer prognosis in breast, lung and liver cancer patients


In spite of increasing advancements in the management of HNSCC, the long-term survival rate remains unchanged over many decades at about 50%. Currently the mainstay of treatment for amenable cancers is surgical resection and reconstruction followed by adjunctive radiotherapy whilst other tumours can be managed using a combination of chemo and radiotherapy. Unlike breast and lung cancer patients, all HNSCC sufferers are subjected to the same combinations of treatment irrespective of the genetic makeup of their cancer. This is primarily because of the gap in our knowledge regarding molecular biomarkers that can be employed to stratify sub-populations and indicate the most suitable intervention based on the molecular profile of the individual tumour. Some progress towards stratification of HNSCC treatments has been made with the histological identification of HPV-driven HNSCCs as potentially separate clinical entities, but alternative treatment strategies such as deintensification protocols are still in the trial phase (e.g. the ComPARE trial in the UK: Our data suggests there are multiple other molecularly distinct subtypes of HNSCC, with distinct clinical presentations that may also require stratified treatment approaches.

Through bioinformatics data mining and validation using a panel of 8 normal oral keratinocyte lines and 10 HNSCC cell lines, we identified 6 candidate biomarkers that were differentially expressed, of which 3 (PLAU, FN1, CDCA5) were upregulated and 3 (CRNN, CLEC3B and DUOX1) downregulated in HNSCC. These 6 gene expression profiles in a cohort of 100 HNSCC patients’ tissues, were consistent with the cell line data in the majority of samples (+q6). From clinical correlation this +q6 cohort were found to be predominantly older males who consumed more than the recommended units of alcohol per week which is representative of the majority of HNSCCs found in the UK population as a whole. The alternative genetic profile was an unexpected inverse expression of q6 biomarkers (−q6) observed in about 20% of patients. This study utilised tissue samples collected from an area of East London with the highest concentration of Bangladeshi individuals in the Western world with specific cultural risk factors for HNSCC, such as ‘paan’ (areca nut) usage. This population is known to be at the higher risk of developing HNSCC compared with rest of the UK population and particularly at a younger age and mostly in women due to these cultural influences. Therefore, the finding that patients with -q6 values in this study are markedly different from the positive group in age (significantly younger) and gender (more females) as well as paan users is extremely interesting as it suggests that the q6 response confirms the epidemiological data, indicating that there are two distinct types of HNSCC in our study population. This is further supported by the statistically significant lower reported alcohol consumption in the -q6 group as alcohol usage is low in Muslim cultures while it is a significant issue in the wider UK population.

Of note also is that all recurrent HNSCC were in the -q6 group, although this finding was not statistically significant. This is potentially important as it suggests that the -q6 group could be more prone to recurrence. In addition to the small sample size, it must be noted that to be included in this study the recurrence must have been treated by surgical excision and therefore operable, so this finding must be interpreted with extreme caution but warrants further investigation.

Tissue samples for this study were collected at the time of tumour resection. The decision to treat surgically was made on presenting clinicopathological factors by a multi-disciplinary team of surgeons, oncologists and allied health professionals in conjunction with the patient’s wishes as per NHS best-practice policy. This would suggest the study sample represents a proportion of all head and neck tumours diagnosed in the study period with other tumours being either inoperable (e.g. due to size, position, metastatic spread, patient’s general health or patient wishes) or managed without a tissue sample being generated (e.g. chemo-radiotherapy or laser ablation). Patients in both groups were managed in the same manner, and our genetic analysis indicates that molecularly the tumours in each group were different and showed distinct expression of q6 biomarkers. This suggests that q6 biomarkers can be used to stratify HNSCC patients based on their molecular signatures.

Interpretation of clinical data from a retrospective study of this type must be done with caution particularly when assessing patient treatment modalities, which cannot ethically be influenced by the study design. The various forms of bias and presence of unknown confounders are a significant concern in this study as is the small sample size augmented by the inability to collect patient details on a proportion of each group. Nevertheless, our compelling data indicates the need for a prospective observational study of the correlation between patient factors and HNSCC treatment response.

A further study conducted in our laboratory on a gamma-irradiated resistant oral keratinocyte cell line demonstrated that the downregulation of PLAU, FN1 and CDCA5 appeared to be indicators of the tumour being resistant to radiation therapy (data not shown). Although this finding needs to be verified and validated by further study, the fact that these biomarkers were able to identify tumours that are potentially resistant or responsive to radiotherapy is potentially an important clinical finding as it identifies patients who, for example, may not respond to one treatment modality and therefore would benefit from personalised alternatives.

Further support for the prognostic significance of the q6 markers came from our analyses on breast, ovarian, lung, gastric and liver cancers using the Kaplan Meier plotter transcriptome database [24]. Overall (with exceptions), the 3 upregulated genes (PLAU1, FN1 and CDCA5) were generally associated with poor prognosis in these cancer types when gene expression levels were upregulated. Similarly, the 3 downregulated genes (CRNN, CLEC3B and DUOX1) were generally associated with poor prognosis when downregulated. These further confirms that the 3 upregulated genes tend to be oncogenes whilst the 3 downregulated genes tend to be tumour suppressor genes. Interestingly, gastric cancer showed inverse relationship for CDCA5, CRNN and DUOX1 on prognosis compared to other cancer types investigated. These results indicated that the q6 biomarkers may also have prognostic significance in many other cancer types.

We further looked into the literature to understand the functional significance of q6 biomarkers in HNSCC development and progression. PLAU (plasminogen activator, urokinase) encodes a serine protease involved in degradation of the extracellular matrix facilitating tumour cell migration and proliferation. PLAU has been shown to be a novel biomarker with high tumour expression levels in HNSCC and is linked to decreased survival rate, increased disease progression and relapse [25]. In prostate cancer and laryngeal squamous cell carcinoma, PLAU gene amplification was preferentially found in advanced stage, but not detected in benign lesions, suggesting PLAU may have a tumour stage-dependent expression pattern [26, 27].

FN1 (fibronectin-1) encodes two forms of fibronectin, soluble plasma fibronectin-1 and insoluble cellular fibronectin-1. The insoluble cellular fibronectin is involved in cell adhesion and migration processes including embryogenesis, wound healing, host defence and metastasis. FN1 is also a downstream target of SATB1 oncogene and is up regulated in salivary ductal carcinoma [28], oesophageal squamous cell carcinoma resulting in enhanced cell proliferation and migration [29]. FN1 also induces metalloproteinases, such as MMP9/MMP2 to promote invasion and metastasis [30,31,32].

The biomarker CDCA5 (cell division cycle associated 5 or human sororin gene) is involved in sister chromatid cohesion, separation and tumourigenesis [33]. A study on lung carcinoma has shown high levels of CDCA5 and its association with poor prognosis [34]. CDCA5 was found to be upregulated in 4 OSCC cell lines and its knockdown led to tumour cell growth inhibition in vitro and in vivo. The same study also found that high levels of CDCA5 immunostaining in OSCC tissues correlated significantly with poorer overall survival [35]. This suggests that CDCA5 has a significant role in OSCC progression, targeting CDCA5 may be a potentially useful diagnostic and therapeutic approach for OSCC patients.

CRNN (cornulin) also known as squamous epithelial heat shock protein 53 belongs to the “fused gene” family and is involved in epithelial immune response and differentiation [36]. In oesophageal squamous cell carcinoma, it is 5-fold downregulated in 89% of cases during transformation from normal to neoplastic cells [37]. Significant loss of CRNN expression is associated with advanced stage, invasiveness, lymph node metastasis and poor survival [37,38,39]. CRNN expression is reported to be downregulated in HNSCC [40] through loss of heterozygosity and microsatellite instability [41]. These findings highlight the role of CRNN in tumour progression and a possible prognostic marker to predict disease outcome.

CLEC3B (C-type lectin domain family 3, member B) encodes tetranectin protein which is a potential biomarker for metastatic oral cancer. Decreased levels of tetranectin have been assoiciated with cancer progression [42]. In ovarian and breast cancer, decreased serum levels of CLEC3B have been associated with poor treatment response [43, 44]. These findings support that CLEC3B may be used as a biomarker for metastasis.

DUOX1 (dual oxidase 1) encodes a glycoprotein and is a member of the NADPH oxidase family. This protein generates hydrogen peroxide and plays a role in antimicrobial defense at mucosal surfaces. It has been found that in 50% of lung cancers NADPH oxidase DUOX1 and DUOX2 go under epigenetic silencing via hypermethylation of CpG-rich promoter regions. Introducing normal levels of DUOX1 into lung cancer cell lines increased cell migration and wound repair without affecting cell growth [45]. The prognostic value of DUOX1 expression is highlighted in liver cancer with low levels of expression, while normal levels were indicative of disease-free survival [46]. To date the potential use of DUOX1 as a diagnostic or prognostic tool has not been explored in HNSCC.


We present the first reported correlation of distinct molecular signatures in HNSCC with the clinical presentation of the disease. Larger scale longitudinal studies are now warranted to establish the linkage between these different molecular subtypes and disease progression or treatment response. This is an important step towards the ultimate goal of improving outcomes by utilising personalised molecular-signature-guided treatments for HNSCC patients.

Availability of data and materials

All microarray datasets used in this study are available in the Gene Expression Omnibus and Oncomine databases. Accession numbers for each dataset have been listed within the manuscript. Any supporting data not included in this manuscript or reagents used in this study, which are not commercially available, will be provided to readers following a written request to the corresponding author.



complimentary DNA


Head and neck squamous cell carcinoma


human papilloma virus


lymph-node metastatic


no template control




reverse transcription quantitative polymerase chain reaction


  1. Williams HK. Molecular pathogenesis of oral squamous carcinoma. Mol Pathol. 2000;53(4):165–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Pao W, Chmielecki J. Rational, biologically based treatment of EGFR-mutant non-small-cell lung cancer. Nat Rev Cancer. 2010;10(11):760–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Pinto AC, Ades F, de Azambuja E, Piccart-Gebhart M. Trastuzumab for patients with HER2 positive breast cancer: delivery, duration and combination therapies. Breast. 2013;22(Suppl 2):S152–5.

    Article  PubMed  Google Scholar 

  4. Guerra EN, Rego DF, Elias ST, Coletta RD, Mezzomo LA, Gozal D, De Luca CG. Diagnostic accuracy of serum biomarkers for head and neck cancer: a systematic review and meta-analysis. Crit Rev Oncol Hematol. 2016;101:93–118.

    Article  PubMed  Google Scholar 

  5. Reddy RB, Bhat AR, James BL, Govindan SV, Mathew R, Ravindra DR, Hedne N, Illiayaraja J, Kekatpure V, Khora SS, et al. Meta-analyses of microarray datasets identifies ANO1 and FADD as prognostic markers of head and neck Cancer. PLoS One. 2016;11(1):e0147409.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Meng RW, Li YC, Chen X, Huang YX, Shi H, Du DD, Niu X, Lu C, Lu MX. Aberrant methylation of RASSF1A closely associated with HNSCC, a meta-analysis. Sci Rep. 2016;6:20756.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Cromer A, Carles A, Millon R, Ganguli G, Chalmel F, Lemaire F, Young J, Dembele D, Thibault C, Muller D, et al. Identification of genes associated with tumorigenesis and metastatic potential of hypopharyngeal cancer by microarray analysis. Oncogene. 2004;23(14):2484–98.

    Article  CAS  PubMed  Google Scholar 

  8. Estilo CL, O-charoenrat P, Talbot S, Socci ND, Carlson DL, Ghossein R, Williams T, Yonekawa Y, Ramanathan Y, Boyle JO, et al. Oral tongue cancer gene expression profiling: Identification of novel potential prognosticators by oligonucleotide microarray analysis. BMC Cancer. 2009;9:11.

    Article  PubMed  PubMed Central  Google Scholar 

  9. 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.

    Article  CAS  PubMed  Google Scholar 

  10. Kuriakose MA, Chen WT, He ZM, Sikora AG, Zhang P, Zhang ZY, Qiu WL, Hsu DF, McMunn-Coffran C, Brown SM, et al. Selection and validation of differentially expressed genes in head and neck cancer. Cell Mol Life Sci. 2004;61(11):1372–83.

    Article  CAS  PubMed  Google Scholar 

  11. Schlingemann J, Habtemichael N, Ittrich C, Toedt G, Kramer H, Hambek M, Knecht R, Lichter P, Stauber R, Hahn M. Patient-based cross-platform comparison of oligonucleotide microarray expression profiles. Lab Investig. 2005;85(8):1024–39.

    Article  CAS  PubMed  Google Scholar 

  12. Sengupta S, den Boon JA, Chen IH, Newton MA, Dahl DB, Chen M, Cheng YJ, Westra WH, Chen CJ, Hildesheim A, et al. Genome-wide expression profiling reveals EBV-associated inhibition of MHC class I expression in nasopharyngeal carcinoma. Cancer Res. 2006;66(16):7999–8006.

    Article  CAS  PubMed  Google Scholar 

  13. Toruner GA, Ulger C, Alkan M, Galante AT, Rinaggio J, Wilk R, Tian B, Soteropoulos P, Hameed MR, Schwalb MN, et al. Association between gene expression profile and tumor invasion in oral squamous cell carcinoma. Cancer Genet Cytogenet. 2004;154(1):27–35.

    Article  CAS  PubMed  Google Scholar 

  14. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Teh MT, Hutchison IL, Costea DE, Neppelberg E, Liavaag PG, Purdie K, Harwood C, Wan H, Odell EW, Hackshaw A, et al. Exploiting FOXM1-orchestrated molecular network for early squamous cell carcinoma diagnosis and prognosis. Int J Cancer. 2013;132(9):2095–106.

    Article  CAS  PubMed  Google Scholar 

  16. Ma H, Dai H, Duan X, Tang Z, Liu R, Sun K, Zhou K, Chen H, Xiang H, Wang J, et al. Independent evaluation of a FOXM1-based quantitative malignancy diagnostic system (qMIDS) on head and neck squamous cell carcinomas. Oncotarget. 2016;7(34):54555–63.

    PubMed  PubMed Central  Google Scholar 

  17. Gemenetzidis E, Bose A, Riaz AM, Chaplin T, Young BD, Ali M, Sugden D, Thurlow JK, Cheong SC, Teo SH, et al. FOXM1 upregulation is an early event in human squamous cell carcinoma and it is enhanced by nicotine during malignant transformation. PLoS One. 2009;4(3):e4849.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Teh MT, Gemenetzidis E, Chaplin T, Young BD, Philpott MP. Upregulation of FOXM1 induces genomic instability in human epidermal keratinocytes. Mol Cancer. 2010;9(1):45.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Waseem A, Ali M, Odell EW, Fortune F, Teh MT. Downstream targets of FOXM1: CEP55 and HELLS are cancer progression markers of head and neck squamous cell carcinoma. Oral Oncol. 2010;46(7):536–42.

    Article  CAS  PubMed  Google Scholar 

  20. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4):611–22.

    Article  CAS  PubMed  Google Scholar 

  21. Zhao S, Fernald RD. Comprehensive algorithm for quantitative real-time polymerase chain reaction. J Comput Biol. 2005;12(8):1047–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3(7):RESEARCH0034.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Rhodes DR, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, Barrette T, Pandey A, Chinnaiyan AM. ONCOMINE: a cancer microarray database and integrated data-mining platform. Neoplasia. 2004;6(1):1–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Lanczky A, Nagy A, Bottai G, Munkacsy G, Szabo A, Santarpia L, Gyorffy B. miRpower: a web-tool to validate survival-associated miRNAs utilizing expression data from 2178 breast cancer patients. Breast Cancer Res Treat. 2016;160(3):439–46.

    Article  CAS  PubMed  Google Scholar 

  25. Sepiashvili L, Hui A, Ignatchenko V, Shi W, Su S, Xu W, Huang SH, O'Sullivan B, Waldron J, Irish JC, et al. Potentially novel candidate biomarkers for head and neck squamous cell carcinoma identified using an integrated cell line-based discovery strategy. Mol Cell Proteomics. 2012;11(11):1404–15.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Bloch M, Ousingsawat J, Simon R, Schraml P, Gasser TC, Mihatsch MJ, Kunzelmann K, Bubendorf L. KCNMA1 gene amplification promotes tumor cell proliferation in human prostate cancer. Oncogene. 2007;26(17):2525–34.

    Article  CAS  PubMed  Google Scholar 

  27. Hui L, Yang N, Yang H, Guo X, Jang X. Identification of biomarkers with a tumor stage-dependent expression and exploration of the mechanism involved in laryngeal squamous cell carcinoma. Oncol Rep. 2015;34(5):2627–35.

    Article  CAS  PubMed  Google Scholar 

  28. Zhang Z, Pan J, Li L, Wang Z, Xiao W, Li N. Survey of risk factors contributed to lymphatic metastasis in patients with oral tongue cancer by immunohistochemistry. J Oral Pathol Med. 2011;40(2):127–34.

    Article  PubMed  Google Scholar 

  29. Song G, Liu K, Yang X, Mu B, Yang J, He L, Hu X, Li Q, Zhao Y, Cai X, et al. SATB1 plays an oncogenic role in esophageal cancer by up-regulation of FN1 and PDGFRB. Oncotarget. 2017;8(11):17771–84.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Qian P, Zuo Z, Wu Z, Meng X, Li G, Wu Z, Zhang W, Tan S, Pandey V, Yao Y, et al. Pivotal role of reduced let-7g expression in breast cancer invasion and metastasis. Cancer Res. 2011;71(20):6463–74.

    Article  CAS  PubMed  Google Scholar 

  31. Shibata K, Kikkawa F, Nawa A, Suganuma N, Hamaguchi M. Fibronectin secretion from human peritoneal tissue induces Mr 92,000 type IV collagenase expression and invasion in ovarian cancer cell lines. Cancer Res. 1997;57(23):5416–20.

    CAS  PubMed  Google Scholar 

  32. Moroz A, Delella FK, Lacorte LM, Deffune E, Felisbino SL. Fibronectin induces MMP2 expression in human prostate cancer cells. Biochem Biophys Res Commun. 2013;430(4):1319–21.

    Article  CAS  PubMed  Google Scholar 

  33. Zhang N, Pati D. Sororin is a master regulator of sister chromatid cohesion and separation. Cell Cycle. 2012;11(11):2073–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Nguyen MH, Koinuma J, Ueda K, Ito T, Tsuchiya E, Nakamura Y, Daigo Y. Phosphorylation and activation of cell division cycle associated 5 by mitogen-activated protein kinase play a crucial role in human lung carcinogenesis. Cancer Res. 2010;70(13):5337–47.

    Article  CAS  PubMed  Google Scholar 

  35. Tokuzen N, Nakashiro K, Tanaka H, Iwamoto K, Hamakawa H. Therapeutic potential of targeting cell division cycle associated 5 for oral squamous cell carcinoma. Oncotarget. 2016;7(3):2343–53.

    Article  PubMed  Google Scholar 

  36. Contzler R, Favre B, Huber M, Hohl D. Cornulin, a new member of the "fused gene" family, is expressed during epidermal differentiation. J Invest Dermatol. 2005;124(5):990–7.

    Article  CAS  PubMed  Google Scholar 

  37. Pawar H, Maharudraiah J, Kashyap MK, Sharma J, Srikanth SM, Choudhary R, Chavan S, Sathe G, Manju HC, Kumar KV, et al. Downregulation of cornulin in esophageal squamous cell carcinoma. Acta Histochem. 2013;115(2):89–99.

    Article  CAS  PubMed  Google Scholar 

  38. Chen K, Li Y, Dai Y, Li J, Qin Y, Zhu Y, Zeng T, Ban X, Fu L, Guan XY. Characterization of tumor suppressive function of cornulin in esophageal squamous cell carcinoma. PLoS One. 2013;8(7):e68838.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Hsu PK, Kao HL, Chen HY, Yen CC, Wu YC, Hsu WH, Chou TY. Loss of CRNN expression is associated with advanced tumor stage and poor survival in patients with esophageal squamous cell carcinoma. J Thorac Cardiovasc Surg. 2014;147(5):1612–8 e1614.

    Article  CAS  PubMed  Google Scholar 

  40. Imai FL, Uzawa K, Nimura Y, Moriya T, Imai MA, Shiiba M, Bukawa H, Yokoe H, Tanzawa H. Chromosome 1 open reading frame 10 (C1orf10) gene is frequently down-regulated and inhibits cell proliferation in oral squamous cell carcinoma. Int J Biochem Cell Biol. 2005;37(8):1641–55.

    Article  CAS  PubMed  Google Scholar 

  41. Salahshourifar I, Vincent-Chong VK, Chang HY, Ser HL, Ramanathan A, Kallarakkal TG, Rahman ZA, Ismail SM, Prepageran N, Mustafa WM, et al. Downregulation of CRNN gene and genomic instability at 1q21.3 in oral squamous cell carcinoma. Clin Oral Investig. 2015;19(9):2273–83.

    Article  PubMed  Google Scholar 

  42. Arellano-Garcia ME, Li R, Liu X, Xie Y, Yan X, Loo JA, Hu S. Identification of tetranectin as a potential biomarker for metastatic oral cancer. Int J Mol Sci. 2010;11(9):3106–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Hogdall CK, Soletormos G, Nielsen D, Norgaard-Pedersen B, Dombernowsky P, Clemmensen I. Prognostic value of serum tetranectin in patients with metastatic breast cancer. Acta Oncol. 1993;32(6):631–6.

    Article  CAS  PubMed  Google Scholar 

  44. Hogdall CK, Hogdall EV, Hording U, Clemmensen I, Norgaard-Pedersen B, Toftager-Larsen K. Pre-operative plasma tetranectin as a prognostic marker in ovarian cancer patients. Scand J Clin Lab Invest. 1993;53(7):741–6.

    Article  CAS  PubMed  Google Scholar 

  45. Luxen S, Belinsky SA, Knaus UG. Silencing of DUOX NADPH oxidases by promoter hypermethylation in lung cancer. Cancer Res. 2008;68(4):1037–45.

    Article  CAS  PubMed  Google Scholar 

  46. Chen S, Ling Q, Yu K, Huang C, Li N, Zheng J, Bao S, Cheng Q, Zhu M, Chen M. Dual oxidase 1: a predictive tool for the prognosis of hepatocellular carcinoma patients. Oncol Rep. 2016;35(6):3198–208.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We are thankful to The Facial Surgery Research Foundation - Saving Faces for funding this study as part of the PhD program for FQ (to MTT). We thank Queen Mary Innovations (QMI) and The Rosetrees Trust for providing financial support (to AW). The authors are thankful to the Centre for Immunobiology and Regenerative Medicine (COIRM) for supporting this study.


This study was co-funded by the following: The Facial Surgery Research Foundation - Saving Faces for funding this study as part of the PhD program for FQ (to MTT), Queen Mary Innovations (QMI) and The Rosetrees Trust (to AW). The funding bodies listed here do not have any roles in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations



Study concepts & Study design: MTT. Data acquisition: FQ, AL, HHD, SH, HA. Quality control of data and algorithms: FQ, AL, AW, MTT. Data analysis and interpretation: FQ, AL, AW, MTT. Statistical analysis: AL, MTT. Manuscript preparation: FQ, AL, AM, HM, HD, MTT. Manuscript editing: FQ, AL, AM, HM, HD, MTT. Manuscript review: FQ, AL, AM, HM, HD, MTT. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Muy-Teck Teh.

Ethics declarations

Ethics approval and consent to participate

The use of fresh clinical specimens collected in the UK was approved by the NHS Research Ethics Committee (06/MRE03/69). All tissue samples were previously collected according to local ethical committee-approved protocols and informed patient consent was obtained from all participants. All cell lines used in this study were previously derived from anonymised tissue samples which were consented and ethnically approved for used in research.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Qadir, F., Lalli, A., Dar, H.H. et al. Clinical correlation of opposing molecular signatures in head and neck squamous cell carcinoma. BMC Cancer 19, 830 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: