Complement C1q (C1qA, C1qB, and C1qC) May Be a Potential Prognostic Factor and an Index of Tumor Microenvironment Remodeling in Osteosarcoma

The tumor microenvironment (TME) has important effects on the tumorigenesis and development of osteosarcoma (OS). However, the dynamic mechanism regulating TME immune and matrix components remains unclear. In this study, we collected quantitative data on the gene expression of 88 OS samples from The Cancer Genome Atlas (TCGA) database and downloaded relevant clinical cases of OS from the TARGET database. The proportions of tumor-infiltrating immune cells (TICs) and the numbers of immune and matrix components were determined by CIBERSORT and ESTIMATE calculation methods. Protein-protein interaction (PPI) network construction and Cox regression analysis were conducted to analyze differentially expressed genes (DEGs). The complement components C1qA, C1qB and C1qC were then determined to be predictive factors through univariate Cox analysis and PPI cross analysis. Further analysis found that the levels of C1qA, C1qB and C1qC expression were positively linked to OS patient survival time and negatively correlated with the clinicopathological feature percent necrosis at definitive surgery. The results of gene set enrichment analysis (GSEA) demonstrated that genes related to immune functions were significantly enriched in the high C1qA, C1qB and C1qC expression groups. Proportion analysis of TICs by CIBERSORT showed that the levels of C1qA, C1qB and C1qC expression were positively related to M1 and M2 macrophages and CD8+ cells and negatively correlated with M0 macrophages. These results further support the influence of the levels of C1qA, C1qB and C1qC expression on the immune activity of the TME. Therefore, C1qA, C1qB and C1qC may be potential indicators of remodeling in the OS TME, which is helpful to predict the prognosis of patients with OS and provide new ideas for immunotherapy for OS.