Prognostic Matrisomal Gene Panel and Its Association with Immune Cell Infiltration in Head and Neck Carcinomas

Unsupervised Analysis Reveals the Involvement of Key Immune Response Genes and the Matrisome in Resistance to BRAF and MEK Inhibitors in Melanoma

Melanoma tumors exhibit a wide range of heterogeneity in genomics even with shared mutations in the MAPK pathway, including BRAF mutations. Consistently, adaptive drug resistance to BRAF inhibitors and/or BRAF plus MEK inhibitors also exhibits a wide range of heterogeneous responses, which poses an obstacle for discovering common genes and pathways that can be used in clinic for overcoming drug resistance. This study objectively analyzed two sets of previously published tumor genomics data comparing pre-treated melanoma tumors and BRAFi- and/or MEKi-resistant tumors. Heterogeneity in response to BRAFi and BRAFi/MEKi was evident because the pre-treated tumors and resistant tumors did not exhibit a tendency of clustering together. Differentially expressed gene (DEG) analysis revealed eight genes and two related enriched signature gene sets (matrisome and matrisome-associated signature gene sets) shared by both sets of data. The matrisome was closely related to the tumor microenvironment and immune response, and five out of the eight shared genes were also related to immune response. The PLXNC1 gene links the shared gene set and the enriched signature gene sets as it presented in all analysis results. As the PLXNC1 gene was up-regulated in the resistant tumors, we validated the up-regulation of this gene in a laboratory using vemurafenib-resistant cell lines. Given its role in promoting inflammation, this study suggests that resistant tumors exhibit an inflammatory tumor microenvironment. The involvement of the matrisome and the specific set of immune genes identified in this study may provide new opportunities for developing future therapeutic methods.

Integrative analysis of lysine acetylation-related genes and identification of a novel prognostic model for oral squamous cell carcinoma

Introduction: Oral squamous cell carcinoma (OSCC), which accounts for a high proportion of oral cancers, is characterized by high aggressiveness and rising incidence. Lysine acetylation is associated with cancer pathogenesis. Lysine acetylation-related genes (LARGs) are therapeutic targets and potential prognostic indicators in various tumors, including oral squamous cell carcinoma. However, systematic bioinformatics analysis of the Lysine acetylation-related genes in Oral squamous cell carcinoma is still unexplored. Methods: We analyzed the expression of 33 Lysine acetylation-related genes in oral squamous cell carcinoma and the effects of their somatic mutations on oral squamous cell carcinoma prognosis. Consistent clustering analysis identified two lysine acetylation patterns and the differences between the two patterns were further evaluated. Least absolute shrinkage and selection operator (LASSO) regression analysis was used to develop a lysine acetylation-related prognostic model using TCGA oral squamous cell carcinoma datasets, which was then validated using gene expression omnibus (GEO) dataset GSE41613. Results: Patients with lower risk scores had better prognoses, in both the overall cohort and within the subgroups These patients also had "hot" immune microenvironments and were more sensitive to immunotherapy. Disscussion: Our findings offer a new model for classifying oral squamous cell carcinoma and determining its prognosis and offer novel insights into oral squamous cell carcinoma diagnosis and treatment.

Identification of an individualized therapy prognostic signature for head and neck squamous cell carcinoma

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) are the most common cancers in the head and neck. Therapeutic response-related genes (TRRGs) are closely associated with carcinogenesis and prognosis in HNSCC. However, the clinical value and prognostic significance of TRRGs are still unclear. We aimed to construct a prognostic risk model to predict therapy response and prognosis in TRRGs-defined subgroups of HNSCC.

METHODS

The multiomics data and clinical information of HNSCC patients were downloaded from The Cancer Genome Atlas (TCGA). The profile data GSE65858 and GSE67614 chip was downloaded from public functional genomics data Gene Expression Omnibus (GEO). Based on TCGA-HNSC database, patients were divided into a remission group and a non-remission group according to therapy response, and differentially expressed TRRGs between those two groups were screened. Using Cox regression analysis and Least absolute shrinkage and selection operator (LASSO) analysis, candidate TRRGs that can predict the prognosis of HNSCC were identified and used to construct a TRRGs-based signature and a prognostic nomogram.

RESULT

A total of 1896 differentially expressed TRRGs were screened, including 1530 upregulated genes and 366 downregulated genes. Then, 206 differently expressed TRRGs that was significantly associated with the survival were chosen using univariate Cox regression analysis. Finally, a total of 20 candidate TRRGs genes were identified by LASSO analysis to establish a signature for risk prediction, and the risk score of each patient was calculated. Patients were divided into a high-risk group (Risk-H) and a low-risk group (Risk-L) based on the risk score. Results showed that the Risk-L patients had better overall survival (OS) than Risk-H patients. Receiver operating characteristic (ROC) curve analysis revealed great predictive performance for 1-, 3-, and 5-year OS in TCGA-HNSC and GEO databases. Moreover, for patients treated with post-operative radiotherapy, Risk-L patients had longer OS and lower recurrence than Risk-H patients. The nomogram involves risk score and other clinical factors had good performance in predicting survival probability.

CONCLUSIONS

The proposed risk prognostic signature and Nomogram based on TRRGs are novel promising tools for predicting therapy response and overall survival in HNSCC patients.

Prognostic Neurotransmitter Receptors Genes Are Associated with Immune Response, Inflammation and Cancer Hallmarks in Brain Tumors

Glioblastoma multiforme (GBM) is one of the most aggressive forms of cancer. Neurotransmitters (NTs) have recently been linked with the uncontrolled proliferation of cancer cells, but the role of NTs in the progression of human gliomas is still largely unexplored. Here, we investigate the genes encoding for neurotransmitter receptors (NTRs) by analyzing public transcriptomic data from GBM and LGG (low-grade glioma) samples. Our results showed that 50 out of the 98 tested NTR genes were dysregulated in brain cancer tissue. Next, we identified and validated NTR-associated prognostic gene signatures for both LGG and GBM. A subset of 10 NTR genes (DRD1, HTR1E, HTR3B, GABRA1, GABRA4, GABRB2, GABRG2, GRIN1, GRM7, and ADRA1B) predicted a positive prognosis in LGG and a negative prognosis in GBM. These genes were progressively downregulated across glioma grades and exhibited a strong negative correlation with genes associated with immune response, inflammasomes, and established cancer hallmarks genes in lower grade gliomas, suggesting a putative role in inhibiting cancer progression. This study might have implications for the development of novel therapeutics and preventive strategies that target regulatory networks associated with the link between the autonomic nervous system, cancer cells, and the tumor microenvironment.

The Role of the Extracellular Matrix and Tumor-Infiltrating Immune Cells in the Prognostication of High-Grade Serous Ovarian Cancer

Ovarian cancer is the eighth global leading cause of cancer-related death among women. The most common form is the high-grade serous ovarian carcinoma (HGSOC). No further improvements in the 5-year overall survival have been seen over the last 40 years since the adoption of platinum- and taxane-based chemotherapy. Hence, a better understanding of the mechanisms governing this aggressive phenotype would help identify better therapeutic strategies. Recent research linked onset, progression, and response to treatment with dysregulated components of the tumor microenvironment (TME) in many types of cancer. In this study, using bioinformatic approaches, we identified a 19-gene TME-related HGSOC prognostic genetic panel (19 prognostic genes (PLXNB2, HMCN2, NDNF, NTN1, TGFBI, CHAD, CLEC5A, PLXNA1, CST9, LOXL4, MMP17, PI3, PRSS1, SERPINA10, TLL1, CBLN2, IL26, NRG4, and WNT9A) by assessing the RNA sequencing data of 342 tumors available in the TCGA database. Using machine learning, we found that specific patterns of infiltrating immune cells characterized each risk group. Furthermore, we demonstrated the predictive potential of our risk score across different platforms and its improved prognostic performance compared with other gene panels.