Osteosarcoma is a common and highly malignant tumour in children and teenagers that is characterized by drug resistance and high metastatic potential. Patients often develop pulmonary metastasis and have a low survival rate. However, the mechanistic basis for pulmonary metastasis remains unclear. To identify key gene and pathways associated with pulmonary metastasis of osteosarcoma, the authors downloaded the gene expression dataset GSE85537 and obtained the differentially expressed genes (DEGs) by analyzing high-throughput gene expression in primary tumours and lung metastases. Subsequently, the authors performed gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) enrichment analyses and a protein-protein interaction (PPI) network was constructed and analyzed by Cytoscape software. In total, 2,493 genes were identified as DEGs. Of these, 485 genes (19.45%) were upregulated, and the remaining 2,008 genes (80.55%) were downregulated. The authors identified the predominant GO categories and KEGG pathways that were significantly over-represented in the metastatic OS samples compared with the non-metastatic OS samples. A PPI network was constructed, and the results indicated that ALB, EGFR, INS, IL6, CDH1, FYN, ERBB2, IL8, CXCL12 and RAC2 were the top 10 core genes. The enrichment analyses of the genes involved in the top three significant modules demonstrated that the DEGs were principally related to neuroactive ligand-receptor interaction, the Rap1 signaling pathway, and protein digestion and absorption. Together, these data elucidated the molecular mechanisms of OS patients with pulmonary metastasis and provide potential therapeutic targets. However, further experimental studies are needed to confirm these results.