A Pan-Cancer Analysis Revealing the Dual Roles of Lysine (K)-Specific Demethylase 6B in Tumorigenesis and Immunity

Introduction: Epigenetic-targeted therapy has been increasingly applied in the treatment of cancers. Lysine (K)-specific demethylase 6B (KDM6B) is an epigenetic enzyme involved in the coordinated control between cellular intrinsic regulators and the tissue microenvironment whereas the pan-cancer analysis of KDM6B remains unavailable. Methods: The dual role of KDM6B in 33 cancers was investigated based on the GEO (Gene Expression Omnibus) and TCGA (The Cancer Genome Atlas) databases. TIMER2 and GEPIA2 were applied to investigate the KDM6B levels in different subtypes or stages of tumors. Besides, the Human Protein Atlas database allowed us to conduct a pan-cancer study of the KDM6B protein levels. GEPIA2 and Kaplan–Meier plotter were used for the prognosis analysis in different cancers. Characterization of genetic modifications of the KDM6B gene was analyzed by the cBioPortal. DNA methylation levels of different KDM6B probes in different TCGA tumors were analyzed by MEXPRESS. TIMER2 was applied to determine the association of the KDM6B expression and immune infiltration and DNA methyltransferases. Spearman correlation analysis was used to assess the association of the KDM6B expression with TMB (tumor mutation burden) and MSI (microsatellite instability). The KEGG (Kyoto encyclopedia of genes and genomes) pathway analysis and GO (Gene ontology) enrichment analysis were used to further investigate the potential mechanism of KDM6B in tumor pathophysiology. Results: KDM6B was downregulated in 11 cancer types and upregulated across five types. In KIRC (kidney renal clear cell carcinoma) and OV (ovarian serous cystadenocarcinoma), the KDM6B level was significantly associated with the pathological stage. A high level of KDM6B was related to poor OS (overall survival) outcomes for THCA (thyroid carcinoma), while a low level was correlated with poor OS and DFS (disease-free survival) prognosis of KIRC. The KDM6B expression level was associated with TMB, MSI, and immune cell infiltration, particularly cancer-associated fibroblasts, across various cancer types with different correlations. Furthermore, the enrichment analysis revealed the relationship between H3K4 and H3K27 methylation and KDM6B function. Conclusion: Dysregulation of the DNA methyltransferase activity and methylation levels of H3K4 and H3K27 may involve in the dual role of KDM6B in tumorigenesis and development. Our study offered a relatively comprehensive understanding of KDM6B’s dual role in cancer development and response to immunotherapy.