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
ATG16L
1 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.