BACKGROUND: Osteosarcoma is the most common primary pediatric and adolescent bone malignancy. An imbalance in copper homeostasis caused by copper ion accumulation could increase intracellular toxicity and regulate cancer cell growth. This study aimed to identify long non-coding RNAs (lncRNAs) associated with cuproptosis to predict prognosis and target drug use to improve patient survival.
METHODS: RNA sequencing and relevant clinical information of ninety-three osteosarcoma patients were obtained from the TARGET database. We then identified thirteen prognostic cuproptosis-related lncRNAs(CRLncs) using coexpression and univariate Cox analyses. The prognostic risk model with three CRLncs was constructed using the least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analysis. Patients were divided into low-risk and high-risk subgroups using the median risk score. The tumor microenvironment (TME) and immune status of identified subgroups were analyzed using ESTIMATE, CIBERSORT, MCP-counter, xCELL, EPIC, and ssGSEA analyses. Functional analyses were conducted to elucidate the underlying mechanisms, including GO, KEGG, GSVA, and GSEA analyses. Also, the relationships between the model, tumor immunity, and drug sensitivity were explored. Lastly, the expression level of ZNF37BP, AL353759.1, and AC005034.5 was validated. in vitro
RESULTS: We constructed a model containing three CRLncs (ZNF37BP, AL353759.1, and AC005034.5) and validated its excellent prognostic and predictive power. The AUC curves for 1-year, 3-year, and 5-year survival probabilities were 0.76, 0.84, and 0.89, respectively. Patients in the high-risk group had a shorter overall survival (OS) time than those in the low-risk. The stroma score and cancer-associated fibroblasts (CAFs) were significantly higher in the low-risk group. Immune cells such as T cells CD4 naive, T cells gamma delta, NK cells resting, dendritic cells resting, and mast cells activated were significantly upregulated in the high-risk group. Based on functional analyses, the PI3K-Akt pathway was identified as a critical metabolic pathway in osteosarcoma. Additionally, we obtained three potentially effective drugs for OS: erlotinib, MP470, and WH-4-023 targeting the PI3K-Akt pathway. The expression level of ZNF37BP was significantly elevated in OS cell lines than in normal osteoblast hFOB1.19 cells, and that of ATP7A, LIPT1, AL353759.1, and AC005034.5 were decreased considerably in OS cell lines.
CONCLUSION: Cuproptosis-related lncRNAs are correlated with the CAFs of osteosarcoma, and this could serve as a foundation for OS survival prediction and treatment.
