The CAV family, especially CAV1 and CAV2, is significantly associated with tumor development. In this study, we aimed to explore the pathogenic and prognostic roles of CAV1 and CAV2 in head and neck squamous cell carcinoma (HNSCC) through bioinformatic analysis and verified in human tissue. We analyzed expression profiles of CAV1 and CAV2 in HNSCC and in normal tissues via data from The Cancer Genome Altas. Prognostic significance was examined by Kaplan-Meier survival curve obtained from the Xena browser together with Cox regression analysis. Co-expressed genes were uploaded to GeneMANIA and applied to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses, showing interaction networks. Signaling pathways of CAV1 and CAV2 in HNSCC were analyzed by Gene Set Enrichment Analysis to elucidate potential regulatory mechanisms. Gene-drug interaction network was explored via Comparative Toxicogenomics Database. Immunohistochemistry was performed to verify theoretical results. Compared with normal tissues, expression levels of CAV1 and CAV2 were remarkably higher in HNSCC (p < 0.0001), which independently implies poor OS (CAV1: HR: 1.146, p = 0.027; CAV2: HR: 1.408, p = 0.002). Co-expressed genes (PXN, ITGA3, TES, and MET) were identified and analyzed by FunRich with CAV1 and CAV2, revealing a significant correlation with focal adhesion (p < 0.001), which has a vital influence on cancer progression. GSEA also showed cellular protein catabolic process (ES = 0.42) and proteasome complex (ES = 0.72), which is a key degradation system for proteins involved in oxidatively damaging and cell cycle and transcription, closely correlated with high expression of CAV2 in HNSCC. More importantly, we found the relationship between different immune cell infiltration degrees in the immune micro-environment in HNSCC and expression levels of CAV1/CAV2 (p < 0.0001). Gene-drug interaction network was checked via CTD. Moreover, tissue microarrays verified higher expression levels of CAV1/CAV2 in HNSCC (p < 0.0001), and the high expression subgroup indicated significantly poorer clinical outcomes (p < 0.05). The results revealed that CAV1 and CAV2 are typically upregulated in HNSCC and might predict poor prognosis.
Read full abstract