Screening of candidate key genes associated with human osteosarcoma using bioinformatics analysis

The aim of the present study was to identify the key genes associated with osteosarcoma (OS) using a bioinformatics approach. Microarray data (GSE36004) was downloaded from the Gene Expression Omnibus database, including 19 OS cell lines and 6 normal controls. Differentially expressed genes (DEGs) in the OS cell lines were identified using the Limma package, and differentially methylated regions were screened with methyAnalysis in R. Copy number analysis was performed and genes with copy number gains/losses were further screened using DNAcopy and cghMCR packages. Functional enrichment analyses were performed using the Database for Annotation, Visualization and Integrated Discovery online tool, and protein-protein interactions were identified based on information obtained from the Search Tool for the Retrieval of Interacting Genes database. A total of 47 downregulated genes were screened in hyper-methylated regions, including the fragment crystallizable (Fc) region of immunoglobulin E, high affinity I, receptor for; γ polypeptide (FCER1G), leptin (LEP) and feline Gardner-Rasheed sarcoma viral oncogene homolog (FGR). In addition, a total of 17 upregulated genes, including the TPase family, AAA domain containing 2 (ATAD2) and cyclin-dependent kinase 4 (CDK4), exhibited copy number gains, while 5 downregulated genes, including Rho GTPase activating protein 9 (ARHGAP9) and major histocompatibility complex, class II, DO α (HLA-DOA), exhibited copy number losses. These results indicate that hyper-methylation of FCER1G, LEP, and FGR may serve a crucial function in the development of OS. In addition, copy number alterations of these DEGs, including ATAD2, CDK4, ARHGAP9 and HLA-DOA, may also contribute to OS progression. These DEGs may be candidate targets for the diagnosis and treatment of this disease.