Gene expression profiling analysis of osteosarcoma cell lines

Osteosarcoma (OS) is the most common type of primary bone malignancy and has a poor prognosis. To investigate the mechanisms of osteosarcoma, the present analyzed the GSE28424 microarray. GSE28424 was downloaded from the Gene Expression Omnibus, and included a collective of 19 OS cell lines and four normal bone cell lines, which were used as controls. Subsequently, the differentially expressed genes (DEGs) were screened using the Limma package in Bioconductor. Gene Ontology (GO) and pathway enrichment analysis of the DEGs was performed using the Database for Annotation, Visualization and Integrated Discovery, interactions between the proteins encoded by the DEGs were identified using STRING, and the protein-protein interaction (PPI) network was visualized using Cytoscape. In addition, modular analysis of the PPI network was performed using the Clique Percolation Method (CPM) in CFinder. A total of 1,170 DEGs were screened, including 530 upreguated and 640 downregulated genes. The enriched functions included organelle fission, immune response and response to wounding. In addition, RPL8 was observed to be involved with the ribosomal pathway in module A of the PPI network of the DEGs. PLCG1, SYK and PLCG2 were also involved in the B-cell receptor signaling pathway in module B and the Fc-epsilon RI signaling pathway in module C. In addition, AURKA (degree=39), MAD2L1 (degree=38), CDCA8 (degree=38), BUB1 (degree=37) and MELK (degree=37) exhibited higher degrees of connectivity in module F. The results of the present study suggested that the RPL8, PLCG1, PLCG2, SYK, MAD2L1, AURKA, CDCA8, BUB1 and MELK genes may be involved in OS.