Identification of key genes in osteosarcoma by meta-analysis of gene expression microarray

Osteosarcoma (OS) is one of the most malignant tumors in children and young adults. To better understand the underlying mechanism, five related datasets deposited in the Gene Expression Omnibus were included in the present study. The Bioconductor ‘limma’ package was used to identify differentially expressed genes (DEGs) and the ‘Weighted Gene Co-expression Network Analysis’ package was used to construct a weighted gene co-expression network to identify key modules and hub genes, associated with OS. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes overrepresentation analyses were used for functional annotation. The results indicated that 1,405 genes were dysregulated in OS, including 927 upregulated and 478 downregulated genes, when the cut off value was set at a ≥2 fold-change and an adjusted P-value of P<0.01 was used. Functional annotation of DEGs indicated that these genes were involved in the extracellular matrix (ECM) and that they function in several processes, including biological adhesion, ECM organization, cell migration and leukocyte migration. These findings suggested that dysregulation of the ECM shaped the tumor microenvironment and modulated the OS hallmark. Genes assigned to the yellow module were positively associated with OS and could contribute to the development of OS. In conclusion, the present study has identified several key genes that are potentially druggable genes or therapeutics targets in OS. Functional annotations revealed that the dysregulation of the ECM may contribute to OS development and, therefore, provided new insights to improve our understanding of the mechanisms underlying OS.