Stemness and PD-L1 expression analysis of the risk score and enrichment analysis

We divided patients into four subgroups based on divergent immune cell infiltration, with the purpose to explore whether there is significant distinction among their risk scores. In addition, to explore the relationship between stemness and the risk score of our model, correlation analysis was preformed, which involved the relation between DNA methylation and risk score, and relation between RNA expression and risk score. Besides, other correlation analyses were also implemented to investigate the relationship of risk score with immune sore and stromal score. Due to the extensive participation of PD-L1 in cancer progression, we further compared the PD-L1 expression level in high and low-risk groups, which was created in box plots. Also, correlation analysis about risk score and PD-L1 expression was performed.

Hallmark enrichment analysis of DEGs (differentially expressed genes) was conducted to obtain the main signaling pathways in which these genes in the high-risk group enriched compared with those in the low-risk group. Additionally, we performed KEGG analysis to explore the biological processes in which prognostic genes in high-risk group involve in comparison with those in low-risk group.