PENK inhibits osteosarcoma cell migration by activating the PI3K/Akt signaling pathway
Journal of Orthopaedic Surgery and Research volume 15, Article number: 162 (2020)
This article reports the effects of proenkephalin (PENK) on osteosarcoma (OS) cell migration.
A Gene Expression Omnibus (GEO) dataset was used to identify differentially expressed genes (DEGs) in OS tumor samples and normal human osteoblasts. Tumor tissue and adjacent normal tissue were collected from 40 OS patients. MG63 cells were transfected with si-PENK. Transwell migration assays and wound healing assays were performed to compare the effect of PENK on migration. Moreover, LY294002 was used to identify the potential mechanism. Gene expression was examined via qRT-PCR and Western blotting.
Bioinformatic analysis revealed that PENK was downregulated in OS tumor samples compared with normal human osteoblasts. Moreover, PENK was identified as the hub gene of the DEGs. The PI3K/Akt signaling pathway was significantly enriched in the DEGs. Moreover, PENK was downregulated in OS and MG63 cells compared with the corresponding control cells. Silencing PENK promoted MG63 cell migration; however, treatment with LY294002 partially attenuated PENK silencing-induced OS cell migration.
PENK inhibits OS cell migration by activating the PI3K/Akt signaling pathway.
Osteosarcoma (OS), a tumor of mesenchymal origin, is a common bone malignancy [1, 2]. OS is characterized by direct formation of osteoid tissue or immature bone by tumor cells  and accounts for approximately 60% of pediatric bone tumors . Traditional treatment strategies for OS mainly include medical approaches such as multidrug chemotherapy and wide surgical excision . However, the inhibitory effect of these strategies on the progression of OS is limited. The median survival time of OS patients is 13 months, and the 2-year and 5-year survival rates of these patients are 44% and 20%, respectively .
Recent advances in molecular genetic studies have improved our understanding of OS pathogenesis. Exploring the hub genes involved in OS could help identify therapeutic targets for OS. Previous studies identified FGF5 as an early diagnostic marker for OS . In addition, Yang et al.  revealed that CCND1 expression was associated with poor prognosis of OS. However, no core gene that could be a therapeutic target for OS has been identified.
At the molecular level, multiple genes are differentially expressed during OS progression . Proenkephalin (PENK) was originally shown to be expressed in the mature nervous and neuroendocrine systems through the opioid pathway and to be involved in the regulation of cell death and survival . PENK has been identified as a key gene in multiple types of cancers [11,12,13].
However, the molecular mechanisms of PENK in OS remain elusive. In this study, the microarray dataset GSE12865 was downloaded from the Gene Expression Omnibus (GEO) repository, and bioinformatic analysis was performed. PENK was identified as a hub gene involved in the progression of OS. Moreover, our results suggested that PENK promotes OS cell migration by regulating the phosphatidylinositol 3-kinase (PI3K)/Akt signaling pathway.
Materials and methods
Differential gene analysis
The gene expression profile GSE12865, which is based on the Affymetrix Human Gene 1.0 ST Array platform, was downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/). The GSE12865 dataset contains 14 samples—12 OS tumor samples and 2 normal human osteoblast (HOB) samples. The limma R/Bioconductor software package was used in R (version 3.5.3) to identify the differentially expressed genes (DEGs) between the OS tumor samples and HOBs. The cutoff criteria were set as follows: |log fold change| > 2 and adjusted P value < 0.05. Volcano plots and heat maps for visualization of the DEGs were exported from R (version 3.5.3).
Gene Ontology (GO) term and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses
To explore the functions of significantly differentially expressed mRNAs and the corresponding signaling pathways, GO term and KEGG pathway enrichment analyses were conducted. GO term enrichment analysis was used to annotate the enrichment of the DEGs in gene functional terms in the molecular function (MF), biological process (BP), and cellular component (CC) categories. KEGG pathway enrichment analysis was used to identify signaling pathways enriched in the DEGs . The clusterProfiler package was used to perform GO and KEGG enrichment analyses of DEGs . P < 0.05 was considered to indicate statistical significance.
The Search Tool for the Retrieval of Interacting Genes (STRING) database (http://string-db.org/) was used to generate the protein-protein interaction (PPI) network dataset. A minimum required interaction score of ≥ 0.4 was selected as the cutoff value. The Molecular Complex Detection (MCODE) Cytoscape software plugin (version 3.7.1) was used to generate subclusters in the PPI network . The advanced options were set as follows: degree cutoff = 2, node score cutoff = 0.2, and k-core = 2.
Patients and tissues
A total of 40 pairs of OS tissue and matched adjacent normal tissue samples were obtained during surgery between January 2013 and January 2015. The inclusion criteria were as follows: (1) primary diagnosis of OS without any history of therapy, (2) no presence of other cancer-related diseases, (3) availability of data for final pathological diagnosis, and (4) no presence of systemic diseases, such as cardiovascular diseases. In addition, no restrictions were placed on the patients’ age or sex or on the anatomic site of the tumor. The patients’ information is summarized in Table 1. The patients did not differ significantly in age, sex ratio, tumor size, tumor anatomic location, serum level of lactate dehydrogenase, serum level of alkaline phosphatase, pathology, or metastasis status. However, a statistically significant difference in the clinical stage was found between patients with high expression of PENK and those with low expression of PENK (P < 0.05). Informed written consent forms were signed by all included patients, and the study protocol was approved by the Institutional Research Ethics Committee of Huangshi Central Hospital.
Human OS cell lines (MG63, U2OS, 143B, and SAOS2) and an immortalized human fetal osteoblast cell line (hFOB1.19) were obtained from the American Type Culture Collection (ATCC, Manassas, VA, USA). All of these cells were maintained in RPMI-1640 medium containing 10% fetal bovine serum (FBS, Gibco; Thermo Fisher Scientific) at 37 °C in a humidified atmosphere of 5% CO2. The culture medium was replaced every 2 days. Confluent monolayers were split by treatment with phosphate-buffered saline (PBS) and 0.05% trypsin solution containing 0.05% ethylenediaminetetraacetic acid (EDTA) (Nanjing Jiancheng Biotechnology Institute). When the cell monolayer was between 70% and 80% confluent, the cells were treated with small interfering RNA (siRNA) negative control (si-con), si-PENK-1, and si-PENK-2.
Human OS cell lines and hFOB1.19 suspensions (1 × 104 cells in 150 μL) were plated on coverslips in a confocal petri dish (NEST, Hong Kong, China) for transient transfection. DMEM (1 mL) containing 10% FBS was added to the dish, and the cells were cultured for 24 h prior to transfection. Transfection was performed using Lipofectamine 3000 (Thermo Fisher Scientific) and Opti-MEM reduced-serum medium (Life Technologies, Waltham, Massachusetts, USA) according to the manufacturer’s instructions.
siRNA sequences targeting PENK (si-PENK-1 and si-PENK-2) were purchased from GenePharma, Inc. The siRNA sequences were as follows: si-PENK-1 forward, 5′-GCAATCGAGATGGAACCAT-3′ and si-PENK-1 reverse, 5′-ATGGTTCCATCTCGATTGC-3′; si-PENK-2 forward, 5′-CCATCGGTCTACCTATTAT-3′ and si-PENK-2 reverse, 5′-TACGAACGACCAGCTTACC-3′; and si-con forward, 5′-UUCUCCGAACGUGUCACGUTT-3′ and si-NC reverse, 5′-ACGUGACACGUUCGGAGAATT-3′. LY294002 was purchased from MedChemExpress (MCE, Shanghai, China) and used at a concentration of 10 μM as previously described . MG63 cells were pretreated with LY294002 (10 μM) and were then stimulated with si-PENK.
Cell Counting Kit-8 (CCK-8) assay
The proliferation ability of the MG63 OS cell line was assessed using a CCK-8 assay (Dojindo, Japan). In brief, cells were plated in 96-well plates, and after incubation for 24 h, 48 h, and 96 h, CCK-8 solution was added. Then, the absorbance at 450 nm was measured with a spectrophotometer (Thermo Fisher Scientific).
Cell migration assay
Cells were seeded in 6-well plates and grown to confluence. The cell monolayers were scratched with a micropipette tip and cultured in FBS-free medium. Photomicrographs were acquired at 0 and 24 h after wounding. Each wound healing assay was independently repeated three times.
Transwell migration chambers (Corning Costar; pore size 8 mm) were used to assess the cell migration ability. A total of 5 × 104 cells were suspended in 100 μL of FBS-free medium and added to the top chamber. The bottom chamber was filled with standard medium. Cells were incubated at 37 °C for 24 h. Cells that failed to migrate or invade were removed from the top surface of the membranes by wiping. Migrated and invaded cells on the bottom surface of the membranes were fixed and stained with 0.2% crystal violet. The migrated or invaded cells from five random areas of the membrane were counted using light microscopy. All experiments were independently performed three times.
Reverse transcription PCR (RT-PCR)
Cells were extracted using TRIzol reagent (Invitrogen; Thermo Fisher Scientific). The concentration and purity of total RNA was determined by measuring the absorbance at 260 and 280 nm (A260/A280). Single-strand cDNA was synthesized by reverse transcription with a Superscript II Reverse Transcriptase Kit and oligo dT (12–18) primers (Invitrogen, Carlsbad, CA, USA) and amplified by RT-PCR with iQ SYBR Green Supermix (Bio-Rad, Hercules, CA, USA). Quantitative PCR (qPCR) was performed using SYBR Green PCR Master Mix (Applied Biosystems) following the manufacturer’s instructions. The expression levels of the target genes were calculated relative to that of the housekeeping gene GAPDH using the Stratagene Mx3000P software (Applied Biosystems) with the 2−ΔΔCt method. The primer sequences used for PCR are shown in Table 2.
Western blotting (WB)
Total protein was extracted with RIPA buffer from cells and was then separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). The separated proteins were transferred to polyvinylidene difluoride membranes, which were blocked with 5% nonfat milk dissolved in Tris-buffered saline containing 0.05% Tween 20 (TBS-T) for 1 h at room temperature. Membranes were then incubated overnight at 4 °C with primary antibodies specific for PENK (Abcam, 1:2000), PI3K (CST 8486, 1:1000), Akt (Abcam ab8448, 1:1000), p-Akt (Abcam ab8564, 1:1000), and GAPDH (Proteintech, 1:8000). Finally, membranes were incubated for 1 h with horseradish peroxidase (HRP)-conjugated secondary antibodies (1:500) and visualized using an enhanced chemiluminescence system according to the manufacturer’s instructions.
All statistical analyses were performed using the SPSS version 17.0 software. GraphPad Prism 7 was used to generate the graphs in this article. The data are expressed as the mean ± standard deviation (SD) values. Student’s t test or one-way analysis of variance was conducted to compare differences between groups. *P < 0.05 was considered to indicate a significant difference.
Bioinformatic analysis results
To decrease the fraction of false-positive findings, 14 sample values were log2 transformed and normalized by quantile normalization. The normalized data are shown in Fig. 1a. The volcano plots and the heatmap of the top 50 DEGs identified from GSE12865 are shown in Fig. 1b, c, respectively. A total of 568 DEGs were identified, among which 245 mRNAs were downregulated and 223 were upregulated. To explore the potential biological functions of the consensus genes, GO term analysis (Fig. 1d) and KEGG pathway analysis (Fig. 1e) were performed. The BP enrichment analysis results showed that the DEGs were significantly involved in extracellular matrix organization. The CC enrichment analysis results showed that the DEGs were significantly involved in the extracellular region. The MF enrichment analysis results showed that the DEGs were significantly involved in integrin binding. The KEGG pathway analysis results revealed that the PI3K/Akt signaling pathway was enriched in the DEGs. The PPI results showed that PENK exhibited the maximum binding score for interactions with other proteins (Fig. 1f).
PENK is downregulated in OS and MG63 cells
The qPCR results revealed that the relative expression of PENK was decreased in OS tumors compared with normal tissues (Fig. 2a). Moreover, downregulation of PENK gene expression was more obvious in MG63 cells than in the other cell lines examined. Therefore, MG63 cells were used for subsequent experiments (Fig. 2b).
The overall survival analysis results suggested that the overall survival of patients with high PENK expression was statistically significantly better than that of patients with low expression (Fig. 2c, P = 0.041).
Silencing of PENK does not affect MG63 cell proliferation
Both the mRNA (Fig. 3a) and protein expressions (Fig. 3b) of PENK were appreciably downregulated after transfection with si-PENK-1 or si-PENK-2. The silencing efficiency was nearly 85%. We selected si-PENK-1 for the subsequent CCK-8 assay to determine whether si-PENK affects OS cancer cell proliferation. The CCK-8 assay results revealed that the number of OS cancer cells in all different treatment groups increased gradually. However, after 24, 48, and 72 h of treatment, the average number of viable cells did not differ significantly between the si-PENK-1- and si-PENK-2-treated groups and the si-con-treated group (Fig. 3c).
Silencing of PENK promotes MG63 cell migration
The migration ability of OS cancer cells was evaluated by a transwell migration assay and a wound healing assay. Compared with si-con treatment, si-PENK treatment significantly enhanced cell migration (Fig. 4a), as further confirmed by the wound healing assay (Fig. 4b). These results revealed that silencing PENK enhanced the migratory activity of OS cancer cells. Moreover, the si-PENK-promoted migration of MG63 cells was significantly blocked by LY294002 (Fig. 4a, b) compared with that of cells treated with si-PENK alone.
Silencing of PENK promotes MG63 cell migration through the PI3K/Akt signaling pathway
Next, we sought to determine whether silencing PENK affected PI3K/Akt pathway activity in MG63 cells. Silencing PENK significantly increased the phosphorylation of Akt and PI3K compared with that in si-con-treated cells (Fig. 5). Compared with silencing PENK alone, treatment with LY294002, a PI3K inhibitor, decreased the phosphorylation of Akt and the expression of PI3K in MG63 cells (Fig. 5).
OS cell migration plays a crucial role in OS progression, and a better understanding of the potential mechanism of OS will help reduce the number of patients with OS. Here, we provided the first evidence that PENK might inhibit OS migration by downregulating the PI3K/Akt signaling pathway.
In this study, we first used a public dataset (GSE12865) to assess the DEGs between HOBs and OS cells. PENK was identified as the hub gene involved in the progression of OS, and the PI3K/Akt signaling pathway was enriched in the DEGs. Further experiments validated that PENK expression was lower in human OS cells and tissues than in control cells and tissues. Moreover, we revealed that silencing PENK increased the migration of MG63 cells, and that this effect was mediated by the PI3K/Akt signaling pathway.
This study is the first to explore the expression and function of PENK in OS, specifically in human OS cells. PENK is expressed in several cancers and performs diverse biological functions. PENK and PENK-derived peptides are associated with the progression of head and neck cancer  and pancreatic cancer . Tang et al.  found that PENK expression is closely associated with the prognosis of gastrointestinal stromal tumors, indicating that PENK may be a therapeutic target for gastrointestinal stromal tumors. Moreover, further studies conclusively showed aberrant methylation of PENK.
However, the role and potential mechanism of PENK in OS have not been comprehensively explained. First, analysis with the Bioconductor software package revealed that PENK was downregulated in OS tumor cells compared with HOBs. In addition, bioinformatic analysis further revealed that PENK was the hub gene involved in OS progression. Thus, we selected PENK for further evaluation.
The PCR and WB results further confirmed that PENK was downregulated in OS tumors. OS patients with high PENK expression had better overall survival than those with low expression. PENK, which is highly upregulated in the STA-ET-7.2 cell line, has been shown to promote stress-activated apoptosis through transcriptional repression of NF-kappaB- and p53-regulated gene targets. In addition, Aryee et al.  found that PENK was differentially expressed in Ewing sarcoma patients, and Rosen et al.  identified that PENK was expressed in nondifferentiated cells of diverse mesodermal lineages, and thus may play different roles in tumor development.
Next, siRNA was used to characterize downstream gene expression. Bioinformatic analysis revealed that PI3K was the most likely potential target of PENK.
The PI3K/Akt pathway is thought to be one of the most important oncogenic pathways in human cancers . Migration of OS cells plays a crucial role in the progression of malignant tumors, and accumulating evidence indicates that the PI3K/Akt pathway promotes these aggressive behaviors. Musculoskeletal sarcomas, including OS, commonly develop aberrant activation of the PI3K/Akt signaling pathway . Moreover, targeting the PI3K/Akt signaling pathway can be an approach for the treatment of sarcomas . Zhang et al.  also found that inhibiting the PI3K/Akt pathway suppressed the progression of Ewing sarcoma. Here, LY294002, a PI3K inhibitor, was used to block the PI3K-Akt signaling pathway, and the si-PENK-promoted migration of MG63 cells was significantly blocked by LY294002. Taken together, these results indicate that PENK inhibits MG63 cell migration through the PI3K/Akt signaling pathway.
In conclusion, we explored PENK expression in OS and found that PENK might inhibit OS migration by downregulating the PI3K/Akt signaling pathway. Therefore, upregulation of PENK may be a promising treatment approach for OS.
Availability of data and materials
We declare that the materials described in the manuscript will be freely available to all scientists for noncommercial purposes.
Gene Expression Omnibus
Differentially expressed genes
Kyoto Encyclopedia of Genes and Genomes
Search Tool for the Retrieval of Interacting Genes
Molecular Complex Detection
Sodium dodecyl sulfate-polyacrylamide gel electrophoresis
Tris-buffered saline with 0.05% Tween 20
Rickel K, Fang F, Tao J. Molecular genetics of osteosarcoma. Bone. 2017;102:69–79. https://doi.org/10.1016/j.bone.2016.10.017.
Zha Z, Han Q, Liu W, et al. lncRNA GAS8-AS1 downregulates lncRNA UCA1 to inhibit osteosarcoma cell migration and invasion. J Orthop Surg Res. 2020;15:38. https://doi.org/10.1186/s13018-020-1550-x.
Mirabello L, Troisi RJ, Savage SA. Osteosarcoma incidence and survival rates from 1973 to 2004: data from the surveillance, epidemiology, and end Results program. Cancer. 2009;115:1531–43. https://doi.org/10.1002/cncr.24121.
Wu B, Xing C, Tao J. Upregulation of microRNA-23b-3p induced by farnesoid X receptor regulates the proliferation and apoptosis of osteosarcoma cells. 2019, 14: 398 DOI: 10.1186/s13018-019-1404-6.
Goryn T, Pienkowski A, Szostakowski B, et al. Functional outcome of surgical treatment of adults with extremity osteosarcoma after megaprosthetic reconstruction-single-center experience. J Orthop Surg Res. 2019;14:346. https://doi.org/10.1186/s13018-019-1379-3.
Faisham WI, Mat Saad AZ, Alsaigh LN, et al. Prognostic factors and survival rate of osteosarcoma: a single-institution study. Asia-Pacific journal of clinical oncology. 2017;13:e104–10. https://doi.org/10.1111/ajco.12346.
Han D, Wang M, Yu Z, et al. FGF5 promotes osteosarcoma cells proliferation via activating MAPK signaling pathway. Cancer Manag Res. 2019;11:6457–66. https://doi.org/10.2147/cmar.s200234.
Yang G, He F, Duan H, et al. lncRNA FLVCR-AS1 promotes osteosarcoma growth by targeting miR381-3p/CCND1. Onco Targets Ther. 2020;13:163–72. https://doi.org/10.2147/ott.s214813.
Wu H, Zhang J, Dai R, et al. Transferrin receptor-1 and VEGF are prognostic factors for osteosarcoma. J Orthop Surg Res. 2019;14:296. https://doi.org/10.1186/s13018-019-1301-z.
Tegeder I, Geisslinger G. Opioids as modulators of cell death and survival--unraveling mechanisms and revealing new indications. Pharmacol Rev. 2004;56:351–69. https://doi.org/10.1124/pr.56.3.2.
Roperch JP, Incitti R, Forbin S, et al. Aberrant methylation of NPY, PENK, and WIF1 as a promising marker for blood-based diagnosis of colorectal cancer. BMC Cancer. 2013;13:566. https://doi.org/10.1186/1471-2407-13-566.
Salhia B, Kiefer J, Ross JT, et al. Integrated genomic and epigenomic analysis of breast cancer brain metastasis. PLoS One. 2014;9:e85448. https://doi.org/10.1371/journal.pone.0085448.
Zhang Y, Fang L, Zang Y, et al. Identification of core genes and key pathways via integrated analysis of gene expression and DNA methylation profiles in bladder cancer. Med Sci Monit, 2018, 24: 3024-3033 DOI: 10.12659/msm.909514.
Wang X, Yin Z, Zhao Y, et al. Identifying potential prognostic biomarkers in head and neck cancer based on the analysis of microRNA expression profiles in TCGA database. Mol Med Rep. 2020. https://doi.org/10.3892/mmr.2020.10964.
Zhou ZX, Chen XM, Zhang YQ, et al. Comprehensive analysis of long noncoding RNA and mRNA in five colorectal cancer tissues and five normal tissues. Biosci Rep. 2020. https://doi.org/10.1042/bsr20191139.
Cao L, Chen Y, Zhang M, et al. Identification of hub genes and potential molecular mechanisms in gastric cancer by integrated bioinformatics analysis. PeerJ. 2018;6:e5180. https://doi.org/10.7717/peerj.5180.
Zhao Z, Ma X, Ma J, et al. Naringin enhances endothelial progenitor cell (EPC) proliferation and tube formation capacity through the CXCL12/CXCR4/PI3K/Akt signaling pathway. Chem Biol Interact. 2018;286:45–51. https://doi.org/10.1016/j.cbi.2018.03.002.
Zagon IS, Smith JP, McLaughlin PJ. Human pancreatic cancer cell proliferation in tissue culture is tonically inhibited by opioid growth factor. Int J Oncol. 1999;14:577–84. https://doi.org/10.3892/ijo.14.3.577.
Suzuki M, Chiwaki F, Sawada Y, et al. Peripheral opioid antagonist enhances the effect of anti-tumor drug by blocking a cell growth-suppressive pathway in vivo. PLoS One. 2015;10:e0123407. https://doi.org/10.1371/journal.pone.0123407.
Tang D, Lin T, Wang Y, et al. High expression of proenkephalin is associated with favorable outcomes in patients with gastrointestinal stromal tumors. Cancer Manag Res. 2019;11:6681–90. https://doi.org/10.2147/cmar.s202044.
Aryee DN, Niedan S, Ban J, et al. Variability in functional p53 reactivation by PRIMA-1(Met)/APR-246 in Ewing sarcoma. Br J Cancer. 2013;109:2696–704. https://doi.org/10.1038/bjc.2013.635.
Rosen H, Polakiewicz RD, Benzakine S, et al. Proenkephalin A in bone-derived cells. Proc Natl Acad Sci U S A. 1991;88:3705–9. https://doi.org/10.1073/pnas.88.9.3705.
Zhang J, Yu XH, Yan YG, et al. PI3K/Akt signaling in osteosarcoma. Clin Chim Acta. 2015;444:182–92. https://doi.org/10.1016/j.cca.2014.12.041.
Scotlandi K, Picci P. Targeting insulin-like growth factor 1 receptor in sarcomas. Curr Opin Oncol. 2008;20:419–27. https://doi.org/10.1097/CCO.0b013e328302edab.
Lim HJ, Wang X, Crowe P, et al. Targeting the PI3K/PTEN/AKT/mTOR pathway in treatment of sarcoma cell lines. Anticancer Res, 2016, 36: 5765-5771 DOI: 10.21873/anticanres.11160.
Zhang S, Li D, Jiao GJ, et al. miR-185 suppresses progression of Ewing’s sarcoma via inhibiting the PI3K/AKT and Wnt/β-catenin pathways. Onco Targets Ther. 2018;11:7967–77. https://doi.org/10.2147/ott.s167771.
Ethics approval and consent to participate
Ethics approval was provided by The Huangshi Central Hospital, Affiliated Hospital of Hubei Polytechnic University, Edong Healthcare Group, and all patients provided consent for participation in this study.
Consent for publication
All the authors in this manuscript have read and approved the final version submitted, and no author has a conflict of interest in this submission.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhang, Hp., Yu, Zl., Wu, Bb. et al. PENK inhibits osteosarcoma cell migration by activating the PI3K/Akt signaling pathway. J Orthop Surg Res 15, 162 (2020). https://doi.org/10.1186/s13018-020-01679-6