PTHR1 May Be Involved in Progression of Osteosarcoma by Regulating miR-124-3p-AR-Tgfb1i1, miR-27a-3p-PPARG-Abca1, and miR-103/590-3p-AXIN2 Axes.

Abstract

Our previous study has indicated that the parathyroid hormone type 1 receptor (PTHR1) may play important roles in development and progression of osteosarcoma (OS) by regulating Wnt, angiogenesis, and inflammation pathway genes. The goal of this study was to further illuminate the roles of PTHR1 in OS by investigating upstream regulation mechanisms (including microRNA [miRNA] and transcription factors [TFs]) of crucial genes. The microarray dataset GSE46861 was downloaded from the Gene Expression Omnibus database, in which six tumors with short hairpin RNA (shRNA) PTHR1 knockdown (PTHR1.358) and six tumors with shRNA control knockdown (Ren.1309) were collected from mice. Differentially expressed genes (DEGs) between PTHR1.358 and Ren.1309 were identified using the linear models for microarray data (LIMMA) method, and then the miRNA-TF-mRNA regulatory network was constructed using data from corresponding databases, followed by module analysis, to screen crucial regulatory relationships. OS-related human miRNAs were extracted from the curated Osteosarcoma Database. Gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) tool. As a result, the miRNA-TF-mRNA regulatory network, including 1049 nodes (516 miRNA, 25 TFs, and 508 DEGs) and 15942 edges (interaction relationships, such as Pparg-Abca1 and miR-590-3p-AXIN2), was constructed, from which three significant modules were extracted and modules 2 and 3 contained interactions between miRNAs/TFs and DEGs such as miR-103-3p-AXIN2, miR-124-3p-AR-Tgfb1i1, and miR-27a-3p-PPARG-Abca1. miR-27a-3p was a known miRNA associated with OS. Abca1, AR, and miR-124-3p were hub genes in the miRNA-TF-mRNA network. Tgfb1i1 was involved in cell proliferation, Abca1 participated in the cholesterol metabolic process, and AXIN2 was associated with the canonical Wnt signaling pathway. Furthermore, we also confirmed upregulation of miR-590-3p and downregulation of AXIN2 in the mouse OS cell line K7M2-WT transfected with PTHR1 shRNA. In conclusion, PTHR1 may play important roles in progression of OS by activating miR-124-3p-AR-Tgfb1i1, miR-27a-3p-PPARG-Abca1, and miR-103/590-3p-AXIN2 axes.