Transcriptomic analyses reveal the underlying pro-malignant functions of PTHR1 for osteosarcoma via activation of Wnt and angiogenesis pathways

Background: Increasing evidence has indicated parathyroid hormone type 1 receptor (PTHR1) plays important roles for the development and progression of osteosarcoma (OS). However, its function mechanisms remain unclear. The goal of this study was to further illuminate the roles of PTHR1 in OS using microarray data. Methods: Microarray data were available from the Gene Expression Omnibus database under the accession number GSE46861, including six tumors from mice with PTHR1 knockdown (PTHR1.358) and six tumors from mice with control knockdown (Ren. 1309). Differentially expressed genes (DEGs) between PTHR1.358 and Ren. 1309 were identified using the LIMMA method, and then, protein-protein interaction (PPI) network was constructed using data from STRING database to screen crucial genes associated with PTHR1. KEGG pathway enrichment analysis was performed to investigate the underlying functions of DEGs using DAVID tool. Results: A total of 1163 genes were identified as DEGs, including 617 downregulated (Lef1, lymphoid enhancer-binding factor 1) and 546 upregulated genes (Dkk1, Dickkopf-related protein 1). KEGG enrichment analysis indicated upregulated DEGs were involved in Renin-angiotensin system (e.g., Agt, angiotensinogen) and Wnt signaling pathway (e.g., Dkk1), while downregulated DEGs participated in Basal cell carcinoma (e.g., Lef1). A PPI network (534 nodes and 2830 edges) was constructed, in which Agt gene was demonstrated to be the hub gene and its interactive genes (e.g., CCR3, CC chemokine receptor 3; and CCL9, chemokine CC chemokine ligand 9) were inflammation related. Conclusions: Our present study preliminarily reveals the pro-malignant effects of PTHR1 in OS cells may be mediated by activating Wnt, angiogenesis, and inflammation pathways via changing the expressions of the crucial enriched genes (Dkk1, Lef1, Agt-CCR3, and Agt-CCL9).