The head and neck squamous cell carcinoma (HNSCC) transcriptome has been profiled extensively, nevertheless, identifying biomarkers that are clinically relevant and thereby with translational benefit, has been a major challenge. The objective of this study was to use a meta-analysis based approach to catalog candidate biomarkers with high potential for clinical application in HNSCC. Data from publically available microarray series (N = 20) profiled using Agilent (4X44K G4112F) and Affymetrix (HGU133A, U133A_2, U133Plus 2) platforms was downloaded and analyzed in a platform/chip-specific manner (GeneSpring software v12.5, Agilent, USA). Principal Component Analysis (PCA) and clustering analysis was carried out iteratively for segregating outliers; 140 normal and 277 tumor samples from 15 series were included in the final analysis. The analyses identified 181 differentially expressed, concordant and statistically significant genes; STRING analysis revealed interactions between 122 of them, with two major gene clusters connected by multiple nodes (MYC, FOS and HSPA4). Validation in the HNSCC-specific database (N = 528) in The Cancer Genome Atlas (TCGA) identified a panel (ECT2, ANO1, TP63, FADD, EXT1, NCBP2) that was altered in 30% of the samples. Validation in treatment naïve (Group I; N = 12) and post treatment (Group II; N = 12) patients identified 8 genes significantly associated with the disease (Area under curve>0.6). Correlation with recurrence/re-recurrence showed ANO1 had highest efficacy (sensitivity: 0.8, specificity: 0.6) to predict failure in Group I. UBE2V2, PLAC8, FADD and TTK showed high sensitivity (1.00) in Group I while UBE2V2 and CRYM were highly sensitive (>0.8) in predicting re-recurrence in Group II. Further, TCGA analysis showed that ANO1 and FADD, located at 11q13, were co-expressed at transcript level and significantly associated with overall and disease-free survival (p<0.05). The meta-analysis approach adopted in this study has identified candidate markers correlated with disease outcome in HNSCC; further validation in a larger cohort of patients will establish their clinical relevance.

Head and neck squamous cell carcinomas (HNSCC) includes multiple subsites that exhibit differential treatment outcome, which is in turn reflective of tumor stage/histopathology and molecular profile. This study hypothesized that the molecular profile is an accurate prognostic adjunct in patients triaged based on clinico-pathological characteristics. Towards this effect, publically available micro-array datasets (n = 8), were downloaded, classified based on HPV association (n = 83) and site (tongue n = 88; laryngopharynx n = 53; oropharynx n = 51) and re-analyzed (Genespring; v13.1). The significant genes were validated in respective cohorts in The Cancer Genome Atlas (TCGA) for correlation with clinico-pathological parameters/survival. The gene entities (n = 3258) identified from HPV based analysis, when validated in TCGA identified the subset specifically altered in HPV+ HNSCC (n = 63), with three genes showing survival impact (RPP25, NUDCD2, NOVA1). Site-specific meta-analysis identified respective differentials (tongue: 3508, laryngopharynx: 4893, oropharynx: 2386); validation in TCGA revealed markers with high incidence (altered in >10% of patients) in tongue (n = 331), laryngopharynx (n = 701) and oropharynx (n = 404). Assessment of these genes in clinical sub-cohorts of TCGA indicated that early stage tongue (MTFR1, C8ORF33, OTUD6B) and laryngeal cancers (TWISTNB, KLHL13 and UBE2Q1) were defined by distinct prognosticators. Similarly, correlation with perineural/angiolymophatic invasion, identified discrete marker panels with survival impact (tongue: NUDCD1, PRKC1; laryngopharynx: SLC4A1AP, PIK3CA, AP2M1). Alterations in ANO1, NUDCD1, PIK3CA defined survival in tongue cancer patients with nodal metastasis (node+ECS-), while EPS8 is a significant differential in node+ECS-laryngopharyngeal cancers. In oropharynx, wherein HPV is a major etiological factor, distinct prognosticators were identified in HPV+ (ECHDC2, HERC5, GGT6) and HPV-(GRB10, EMILIN1, FNDC1). Meta-analysis in combination with TCGA validation carried out in this study emphasized on the molecular heterogeneity inherent within HNSCC; the feasibility of leveraging this information for

Laryngo-pharyngeal squamous cell carcinomas are one of the most common head and neck cancers. Despite the presence of a large body of information, molecular biomarkers are not currently used in the diagnosis, treatment and management of patients for this group of cancer. Here, we have profiled expression of genes and microRNAs of larynx and hypopharynx tumors using high-throughput sequencing experiments. We found that matrix metalloproteinases along with SCEL, CRNN, KRT4, SPINK5, and TGM3 among others have significantly altered expression in these tumors. Alongside gene expression, the microRNAs hsa-miR-139, hsa-miR-203 and the hsa-miR-424/503 cluster have aberrant expression in these cancers. Using target genes for these microRNAs, we found the involvement of pathways linked to cell cycle, p53 signaling, and viral carcinogenesis significant (P-values 10−13, 10−9 and 10−7 respectively). Finally, using an ensemble machine-learning tool, we discovered a unique 8-gene signature for this group of cancers that differentiates the group from the other tumor subsites of head and neck region. We investigated the role of promoter methylation in one of these genes, WIF1, and found no correlation between DNA methylation and down-regulation of WIF1. We validated our findings of gene expression, 8-gene signature and promoter methylation using q-PCR, data from TCGA and q-MSP respectively.Data presented in this manuscript has been submitted to the NCBI Geo database with the accession number GSE67994.

Sample processing protocols that enable compatible recovery of differentially expressed transcripts and proteins are necessary for integration of the multiomics data applied in the analysis of tumors. In this pilot study, we compared two different isolation methods for extracting RNA and protein from laryngopharyngeal tumor tissues and the corresponding adjacent normal sections. In Method 1, RNA and protein were isolated from a single tissue section sequentially and in Method 2, the extraction was carried out using two different sections and two independent and parallel protocols for RNA and protein. RNA and protein from both methods were subjected to RNA-seq and iTRAQ-based LC-MS/MS analysis, respectively. Analysis of data revealed that a higher number of differentially expressed transcripts and proteins were concordant in their regulation trends in Method 1 as compared to Method 2. Cross-method comparison of concordant entities revealed that RNA and protein extraction from the same tissue section (Method 1) recovered more concordant entities that are missed in the other extraction method (Method 2) indicating heterogeneity in distribution of these entities in different tissue sections. Method 1 could thus be the method of choice for integrated analysis of transcriptome and proteome data.