Systemic inflammatory events and localized disease, mediated by the microbiome, may be measured in saliva as head and neck squamous cell carcinoma (HNSCC) diagnostic and prognostic biomonitors. We used a 16S rRNA V3-V5 marker gene approach to compare the saliva microbiome in DNA isolated from Oropharyngeal (OPSCC), Oral Cavity Squamous Cell Carcinoma (OCSCC) patients and normal epithelium controls, to characterize the HNSCC saliva microbiota and examine their abundance before and after surgical resection.The analyses identified a predominance of Firmicutes, Proteobacteria and Bacteroidetes, with less frequent presence of Actinobacteria and Fusobacteria before surgery. At lower taxonomic levels, the most abundant genera were Streptococcus, Prevotella, Haemophilus, Lactobacillus and Veillonella, with lower numbers of Citrobacter and Neisseraceae genus Kingella. HNSCC patients had a significant loss in richness and diversity of microbiota species (p<0.05) compared to the controls. Overall, the Operational Taxonomic Units network shows that the relative abundance of OTU's within genus Streptococcus, Dialister, and Veillonella can be used to discriminate tumor from control samples (p<0.05). Tumor samples lost Neisseria, Aggregatibacter (Proteobacteria), Haemophillus (Firmicutes) and Leptotrichia (Fusobacteria). Paired taxa within family Enterobacteriaceae, together with genus Oribacterium, distinguish OCSCC samples from OPSCC and normal samples (p<0.05). Similarly, only HPV positive samples have an abundance of genus Gemellaceae and Leuconostoc (p<0.05). Longitudinal analyses of samples taken before and after surgery, revealed a reduction in the alpha diversity measure after surgery, together with an increase of this measure in patients that recurred (p<0.05). These results suggest that microbiota may be used as HNSCC diagnostic and prognostic biomonitors.

Tumor suppressor genes (TSGs) are commonly inactivated by somatic mutation and/or promoter methylation; yet, recent high-throughput genomic studies have not identified key TSGs inactivated by both mechanisms. We pursued an integrated molecular analysis based on methylation binding domain sequencing (MBD-seq), 450K Methylation arrays, whole exome sequencing, and whole genome gene expression arrays in primary head and neck squamous cell carcinoma (HNSCC) tumors and matched uvulopalatopharyngoplasty tissue samples (UPPPs). We uncovered 186 downregulated genes harboring cancer specific promoter methylation including PAX1 and PAX5 and we identified 10 key tumor suppressor genes (GABRB3, HOXC12, PARP15, SLCO4C1, CDKN2A, PAX1, PIK3AP1, HOXC6, PLCB1, and ZIC4) inactivated by both promoter methylation and/or somatic mutation. Among the novel tumor suppressor genes discovered with dual mechanisms of inactivation, we found a high frequency of genomic and epigenomic alterations in the PAX gene family of transcription factors, which selectively impact canonical NOTCH and TP53 pathways to determine cell fate, cell survival, and genome maintenance. Our results highlight the importance of assessing TSGs at the genomic and epigenomic level to identify key pathways in HNSCC, deregulated by simultaneous promoter methylation and somatic mutations.

While specific mutations in the tyrosine kinase domain of epidermal growth factor receptor (EGFR) identify tumors that are responsive to EGFR tyrosine kinase inhibitors (TKIs), these genetic alterations are present in only a minority of patients. Patients with tumors expressing wild-type (wt) EGFR lack reliable predictive markers of their clinical response to EGFR TKIs. Although epithelial-mesenchymal transition (EMT) has been inversely correlated with the response of cancers to EGFR-targeted therapy, the precise molecular mechanisms underlying this association have not been defined and no specific EMT-associated biomarker of clinical benefit has been identified. Here we show that during transforming growth factor-β (TGFβ)-mediated EMT, inhibition of the microRNAs 200 (miR200) family results in upregulated expression of mitogen-inducible gene 6 (Mig6), a negative regulator of EGFR. The Mig6-mediated reduction of EGFR occurs concomitantly with a TGFβ-induced EMT-associated kinase switch of tumor cells to an AKT-activated EGFR-independent state. In a panel of 25 cancer cell lines of different tissue origins, we find that the ratio of the expression levels of Mig6 and miR200c is highly correlated with EMT and resistance to erlotinib. Analyses of primary tumor xenografts of patient-derived lung and pancreatic cancers carrying wild type EGFR showed that the tumor Mig6(mRNA)/miR200 ratio was inversely correlated with response to erlotinib in vivo. Our data demonstrate that the TGFβ-miR200-Mig6 network orchestrates the EMT-associated kinase switch that induces resistance to EGFR inhibitors, and identify a low ratio of Mig6 to miR200 as a promising predictive biomarker of the response of tumors to EGFR TKIs.