GPRASP1 is a candidate anti-oncogene and correlates with immune microenvironment and immunotherapeutic efficiency in head and neck cancer

Background G-protein-coupled receptor-associated sorting protein 1 (GPRASP1) plays an important role in tumorigenesis. However, GPRASP1 specific role has not been clarified in head and neck cancer (HNC). Methods HNC RNA sequencing (RNA-seq) datasets, DNA methylation data, somatic mutation data, copy number variation (CNV) data, and corresponding clinicopathologic information were acquired from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. A comprehensive evaluation was performed to explore the relationship of GPRASP1 expression with clinicopathologic characteristics, CNV, and DNA methylation. Additionally, we employed HNC tissue microarray (TMA) to further confirm the relation between GPRASP1 expression and clinical features. Then, we systematically associated the GPRASP1 with immunological properties from numerous perspectives, such as immune cell infiltration, immune-related pathways, immune checkpoint inhibitors (ICIs), immunomodulators, immunogenicity, and immunotherapy. Results Analyzing TCGA, GEO, and TMA datasets, GPRASP1 is significantly down-regulated in HNC compared to normal tissues. The expression of GPRASP1 is significantly negatively correlated with clinical features (perineural invasion, histologic grade, T stage, and TNM stage), and is an independent predictor of favorable prognosis, regardless of other clinicopathological features (HR: 0.42, 95% CI 0.20-0.91, p = 0.028). The etiological investigation found that the abnormal expression of GPRASP1 was related to DNA methylation, not CMV. Subsequently, the high expression of GPRASP1 was significantly correlated with immune cell infiltration (CD8(+) T cell, tumor infiltrating lymphocyte), immune-related pathways (cytolytic activity, check-point, human leukocyte antigen), ICIs (CTLA4, HAVCR2, LAG3, PDCD1, and TIGIT), immunomodulators (CCR4/5, CXCL9, CXCR3/4/5), and immunogenicity (immune score, neoantigen, tumor mutation burden). Finally, immunophenoscore and tumor immune dysfunction and exclusion analysis demonstrated that GPRASP1 expression levels can accurately predict the immunotherapeutic response. Conclusion GPRASP1 is a promising candidate biomarker that plays a role in the occurrence, development, and prognosis of HNC. Evaluating GPRASP1 expression will aid in the characterization of tumor microenvironment infiltration and orient more efficient immunotherapy strategies.