Refametinib

Osteosarcoma is the most common primary bone malignancy, with a very poor prognosis. Aberrant glycosylation is close involvement in osteosarcoma. Accordingly, this study aimed at investigating the role of glycosylation genes in the prognosis and therapy options of osteosarcoma. The microenvironment of osteosarcoma was assessed using estimate algorithm. A total of 20 immune-related glycosylation genes (IRGGs) was identified using Pearson correlation analysis. Accordingly, osteosarcoma patients were divided into C1 and C2 type using consensus clustering. Multiple algorithms (Xcell, MCP-counter, ssGSEA, epic, quantiseq), cancer immune cycle analysis, and GSVA were applied to estimate the immune, molecule and metabolism characteristics of osteosarcoma, indicating that C1 type was featured with high immune infiltration, high glycosylation, enriched MEK signaling, and good prognosis, while C2 type was characterized by more metastasis, enriched immunotherapy-positive gene signatures, high tumor mutation burden, and poor prognosis. Results from TIDE algorithm and immunotherapy datasets suggested the C2 type's preference of immune checkpoint inhibitors (ICIs), while data of GDSC, CMap analysis and cell experiments indicated that C1 type was sensitivity to MEK inhibitor PD0325901. In addition, univariate Cox and Lasso analysis was combined to establish an IRGGs' risk score containing 6 genes (B3GNT8, FUT7, GAL3ST4, GALNT14, HS3ST2, and MFNG). The data of DCA and ROC indicated its well prediction of prognosis in osteosarcoma. Finally, cellular location analysis showed that the 6 genes not only distributed in tumor cells but also in immune cells. In summary, the classification and risk score based on IRGGs effectively predicted the prognosis and therapy options of osteosarcoma. Further studies on IRGGs may contribute to the understanding of cancer immunity in osteosarcoma.

N6-methyladenosine (m6A) is the most predominant RNA epigenetic regulation in eukaryotic cells. Numerous evidence revealed that m6A modification exerts a crucial role in the regulation of tumor microenvironment (TME) cell infiltration in several tumors. Nevertheless, the potential role and mechanism of m6A modification in nasopharyngeal carcinoma (NPC) remains unknown. mRNA expression data and clinical information from GSE102349, and GSE53819 datasets obtained from Gene Expression Omnibus (GEO) was used for differential gene expression and subsequent analysis. Consensus clustering was used to identify m6A-related molecular patterns of 88 NPC samples based on prognostic m6A regulators using Univariate Cox analysis. The TME cell-infiltrating characteristics of each m6A-related subclass were explored using single-sample gene set enrichment (ssGSEA) algorithm and CIBERSORT algotithm. DEGs between two m6A-related subclasses were screened using edgeR package. The prognostic signature and predicated nomogram were constructed based on the m6A-related DEGs. The cell infiltration and expression of prognostic signature in NPC was determined using immunohistochemistry (IHC) analysis. Chi-square test was used to analysis the significance of difference of the categorical variables. And survival analysis was performed using Kaplan-Meier plots and log-rank tests. The NPC samples were divided into two m6A-related subclasses. The TME cell-infiltrating characteristics analyses indicated that cluster 1 is characterized by immune-related and metabolism pathways activation, better response to anit-PD1 and anti-CTLA4 treatment and chemotherapy. And cluster 2 is characterized by stromal activation, low expression of HLA family and immune checkpoints, and a worse response to anti-PD1 and anti-CTLA4 treatment and chemotherapy. Furthermore, we identified 1558 DEGs between two m6A-related subclasses and constructed prognostic signatures to predicate the progression-free survival (PFS) for NPC patients. Compared to non-tumor samples, REEP2, TMSB15A, DSEL, and ID4 were upregulated in NPC samples. High expression of REEP2 and TMSB15A showed poor survival in NPC patients. The interaction between REEP2, TMSB15A, DSEL, ID4, and m6A regulators was detected. Our finding indicated that m6A modification plays an important role in the regulation of TME heterogeneity and complexity.