The pan-cancer multi-omics landscape of key genes of sialylation combined with RNA-sequencing validation

Sialylation, the process of salivary acid glycan synthesis, plays a pivotal function in tumor growth, immune escape, tumor metastasis, and resistance to drugs. However, the association between sialylation and prognosis, tumor microenvironment (TME), and treatment response in a variety of cancers remains unclear. A comprehensive survey of the expression profile, prognostic value, and genetic and epigenetic alterations of sialylation-related genes was performed in pan-cancer. Subsequently, the single-sample gene set enrichment analysis (ssGSEA) algorithm was used to compute sialylation pathway scores in pan-cancer. Correlations of sialylation pathway scores with clinical features, prognosis, and TME were evaluated using multiple algorithms. Finally, the efficacy of the sialylation pathway score in determining the effect of immunotherapy was evaluated. The expression of sialylation-related genes were verified by RNA-sequencing. Significant differences were observed in sialylation-related genes expression between tumors and adjacent normal tissues for most cancer types. Sialylation pathway scores differed according to the type of tumor, where the poor prognosis was correlated with high sialylation pathway scores in uveal melanoma (UVM) and pancreatic adenocarcinoma (PAAD). In addition, sialylation pathway scores were positively associated with the ImmuneScore, StromalScore and immune-related pathways. Moreover, the level of immune cells infiltration was higher in tumors with higher sialylation pathway scores. Finally, patients with high sialylation pathway scores were more sensitive to immunotherapy. Sialylation-related genes are essential in pan-cancer. The sialylation pathway score may be used as a biomarker in oncology patients.