Abstract. The present study aimed to investigate the possible molecular mechanisms underlying the pathogenesis of metastatic osteosarcoma (OS), by examining the microarray expression profiles of normal samples, and metastatic and non-metastatic OS samples. The GSE9508 gene expression profile was downloaded from the Gene Expression Omnibus database, which included 11 human metastatic OS samples, seven non-metastatic OS samples and five normal samples. Pretreatment of the data was performed using the BioConductor package in R language, and the differentially expressed genes (DEGs) were identified by a t-test. Furthermore, function and pathway enrichment analyses of the DEGs were conducted using a molecule annotation system. A differential co-expression network was also constructed, and the submodules were screened using MCODE in Cytoscape. A total of 965 genes were identified as DEGs in metastatic OS. The DEGs were shown to participate in the regulation of DNA-dependent transcription, the composition of the nucleus, cytoplasm and membrane, and protein and nucleotide binding. Furthermore, the screened DEGs were significantly associated with the ribosome, axon guidance and the cytokine-cytokine receptor interaction pathway. Certain hub genes were identified in the constructed differential co-expression network, including matrix metalloproteinase 1 (MMP1), smoothened (SMO), ewing sarcoma breakpoint region 1 (EWSR1) and fasciculation and elongation protein ζ-1 (FEZ1). Brain selective kinase 2 (BRSK2) and aldo-keto reductase family 1 member B10 (AKRIB10) were present in the screened submodules. The results of the present study suggest that genes, including MMP1, SMO, EWSR1, FEZ1, BRSK2 and AKRIB10, may be potential targets for the diagnosis and treatment of metastatic OS.
IntroductionOsteosarcomas (OS) are among the most frequently occurring secondary malignancies in childhood cancer (1). OS most often originates in the metaphyses of long bones in adolescents and young adults (2). During the past 30 years, an optimal treatment strategy for OS has been developed, which consists of multi-agent chemotherapy and aggressive surgical resection of all sites of disease involvement (3). However, ~80% of patients with localized OS develop metastatic disease following surgical resection (4). Patients with primary metastatic OS are a heterogeneous group, and a five-year event-free survival rate of up to 75% is reported for patients presenting with unilateral lung metastases (5). Furthermore <20% of patients with high-grade osteosarcoma are clinically diagnosed with metastatic disease at the initial diagnosis, and long-term survival rates of patients with metastatic OS range between 10-40% (6,7). Therefore, due to the high rate of systemic spread, complete recovery following surgical treatment alone is rare (8).Recently, microarray analysis (9) has been widely used in screening for differentially expressed genes (DEGs), in order to identify potential genes that may be investigated for the treatment of various human diseases (10,11)...