Integrated analysis of anti-tumor roles of BAP1 in osteosarcoma

BACKGROUND: This study aims to screen out differentially expressed genes (DEGs) regulated by BRCA1-associated protein 1 (BAP1) in osteosarcoma cells, and to analyze their biological functions. METHODS: The microarray dataset GSE23035 of BAP1-knockdown osteosarcoma cells was obtained from Gene Expression Omnibus (GEO) database, consisting of shControl, shBAP1#1 and shBAP1#2 samples. The DEGs between the BAP1-knockdown osteosarcoma cells and the untreated osteosarcoma cells were screened with limma package, and then subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Gene Set Enrichment Analysis (GSEA) was also performed for the three groups of samples. Hub genes in a protein-protein interaction (PPI) network of DEGs was filtered, and then subjected to prognostic analysis and correlation analysis with BAP1 in Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database. Besides, the correlation between BAP1 and biological processes/pathways was analyzed by Gene Set Variation Analysis (GSVA) method and the correlation between BAP1 and immune infiltration by CIBERSORT and ESTIMATE methods. The roles of BAP1 in regulating proliferation and epithelial-mesenchymal transition (EMT) were validated by CCK-8 and western blot. RESULTS: 58 upregulated DEGs and 81 downregulated DEGs were obtained with |logFC| ≥ 1 and adj.p < 0.05. Cell cycle, DNA repair, and focal adhesion were associated with BAP1 in datasets. Further, BAP1 was negatively correlated with naïve CD4 T cells infiltration. In vitro, BAP1 inhibited proliferation and EMT. CONCLUSION: BAP1 might be a tumor suppressor in osteosarcoma and a promising therapeutic target.