Abstract:The present study aimed to understand the molecular mechanisms underlying osteosarcoma metastasis. Microarray dataset GSE49003 was downloaded from the Gene Expression Omnibus database and used for analysis. Raw expression data were preprocessed using the preprocessCore, impute and aggregate packages in R. Differentially expressed genes (DEGs) between metastatic and non‑metastatic osteosarcoma cell lines were screened using the limma package following exclusion of DEGs with a higher significance in intra‑groups compared with inter‑groups using the genefilter package. Enrichment analysis was performed on DEGs using TargetMine, followed by identification of transcription factors (TFs) and microRNAs (miRNAs). Regulatory networks were constructed using Cytoscape software. A total of 248 upregulated and 208 downregulated genes were obtained. The upregulated genes were significantly enriched in the following pathways: Downregulation of transforming growth factor β (TGF‑β) receptor signaling and TGF‑β receptor signaling activates SMADs; these upregulated genes included protein phosphatase 1, regulatory subunit 15A, transforming growth factor, β receptor II and ubiquitin carboxyl‑terminal hydrolase L5. In addition, some upregulated genes were enriched in lung cancer disease ontology, including epidermal growth factor receptor (EGFR), insulin‑like growth factor 2 mRNA binding protein 3 (IGF2BP3), runt‑related transcription factor 3 (RUNX3) and secreted frizzled‑related protein 1 (SFRP1). Conversely, the downregulated genes were significantly enriched in extracellular matrix‑associated pathways or functions, such as collagen, type XII, α 1; collagen, type I, α 1; collagen, type IV, α 1; and collagen, type V, α 1. In addition, some downregulated genes were significantly enriched in the TGF‑β signaling pathway, including bone morphogenetic protein 4, inhibitor of DNA binding 3 and SMAD family member 6. A total of 10 TFs and 84 miRNAs (e.g. miR-21-5p) were deemed to be associated with DEGs. In conclusion, DEGs enriched in the downregulation of TGF‑β receptor signaling, TGF‑β receptor signaling activates SMADs and TGF‑β signaling pathways, as well as the extracellular matrix may be implicated in the progression of osteosarcoma metastasis. Dysregulated EGFR, IGF2BP3, RUNX3 and SFRP1 may contribute to the metastasis of osteosarcoma to the lungs. In addition, screened TFs and miR-21-5p may be associated with metastasis via target genes.

An Integrative Analysis of DNA Methylation in Osteosarcoma

A 3-DNA Methylation Signature As a Novel Prognostic Biomarker in Patients with Sarcoma by Bioinformatics Analysis

Integrated Genome-Wide Methylation and Expression Analyses Reveal Key Regulators in Osteosarcoma

Integral Analyses of Competing Endogenous RNA Mechanisms and DNA Methylation Reveal Regulatory Mechanisms in Osteosarcoma

Identification and Analysis of Key Genes in Osteosarcoma Using Bioinformatics

Development of a Risk Score Model for Osteosarcoma Based on DNA Methylation-Driven Differentially Expressed Genes

Quantitative analysis of promoter hypermethylation in multiple genes in osteosarcoma.

Integrative analysis of immune‐related multi‐omics profiles identifies distinct prognosis and tumor microenvironment patterns in osteosarcoma

Epigenetic Regulation of CXCL12 Plays a Critical Role in Mediating Tumor Progression and the Immune Response in Osteosarcoma

Exploring the Molecular Mechanisms of Osteosarcoma by the Integrated Analysis of Mrnas and Mirna Microarrays.

The Regulatory Role of Aberrant Methylation of Microrna-34A Promoter CpGs in Osteosarcoma

Identification of Tumorigenicity-Associated Genes in Osteosarcoma Cell Lines Based on Bioinformatic Analysis and Experimental Validation.

Identification Of Differentially Expressed Genes Induced By Aberrant Methylation In Oral Squamous Cell Carcinomas Using Integrated Bioinformatic Analysis

Molecular characterization of metastatic osteosarcoma: Differentially expressed genes, transcription factors and microRNAs

Comprehensive investigation of tumor immune microenvironment and prognostic biomarkers in osteosarcoma through integrated bulk and single-cell transcriptomic analysis

Identification of Key Genes and Pathways in Osteosarcoma by Bioinformatics Analysis

Integrated Analysis of Gene Expression and Genomic Aberration Data in Osteosarcoma (OS)

Identification of CD20, ECM, and ITGA as Biomarkers for Osteosarcoma by Integrating Transcriptome Analysis

The role of SPI1-TYROBP-FCER1G network in oncogenesis and prognosis of osteosarcoma, and its association with immune infiltration

An Integrated Methylation and Gene Expression Microarray Analysis Reveals Significant Prognostic Biomarkers in Oral Squamous Cell Carcinoma

Integrated analysis of gene expression and DNA methylation profiles in ovarian cancer