Abstract
We identified the microRNA (miRNA) expression signature of head and neck squamous cell carcinoma (HNSCC) tissues by RNA sequencing, in which 168 miRNAs were significantly upregulated, including both strands of the miR-31 duplex (miR-31-5p and miR-31-3p). The aims of this study were to identify networks of tumor suppressor genes regulated by miR-31-5p and miR-31-3p in HNSCC cells. Our functional assays showed that inhibition of miR-31-5p and miR-31-3p attenuated cancer cell malignant phenotypes (cell proliferation, migration, and invasion), suggesting that they had oncogenic potential in HNSCC cells. Our in silico analysis revealed 146 genes regulated by miR-31 in HNSCC cells. Among these targets, the low expression of seven genes (miR-31-5p targets: CACNB2 and IL34; miR-31-3p targets: CGNL1, CNTN3, GAS7, HOPX, and PBX1) was closely associated with poor prognosis in HNSCC. According to multivariate Cox regression analyses, the expression levels of five of those genes (CACNB2: p = 0.0189; IL34: p = 0.0425; CGNL1: p = 0.0014; CNTN3: p = 0.0304; and GAS7: p = 0.0412) were independent prognostic factors in patients with HNSCC. Our miRNA signature and miRNA-based approach will provide new insights into the molecular pathogenesis of HNSCC.