Abstract
Background: Osteosarcoma is a common solid malignancy of the bone in children and adolescents, and its metastasis and recurrence are the principal causes of poor treatment outcomes. Methods: Autophagy-related genes were used to cluster osteosarcoma patients by consensus clustering analysis using the GSE21257 database. Differentially expressed genes (DEGs) were identified by limma package. Multiple-gene risk signature was constructed using least absolute shrinkage and selection operator (LASSO) analysis and Cox regression analyses. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was used to determine gene expression levels. Then, single-cell RNA-sequencing dataset GSE152048 were used to identify the correlation between the DEGs and effector molecules expressed in specific tumor-infiltrating immune cells. Results: Two clusters were identified in the consensus clustering analysis, which were confirmed by principal component analysis. Limma analysis revealed that 15 genes were related, and 9 genes were screened using protein-protein interaction network and LASSO regression analysis. Cox regression analyses identified 5 genes. Combined with survival analysis, only the autophagy related 16 like 1 gene (ATG16L1) was significant. The results of qRT-PCR showed low expression levels of ATG16L1 in tumor cells group. Immune infiltration analysis revealed significantly lower expression of CD8+ T cells in the high ATG16L1 gene expression group. ScRNA-seq revealed that in the ATG16L1+ CD8 + T cell group, the expression of GZMB was lower, whereas the expression of ITGA1 was higher. These results showed that ATG16L1 is an immune-related gene, which is associated with poor prognosis in patients with osteosarcoma. Conclusion: ATG16L1 is a potential prognostic biomarker and immune signature and may be a therapeutic target for osteosarcoma.