Expression of col6a1 in sarcoma tissue and its significance
HU Weiquan,ZHAO Qinfei,ZHU Longyu,HU Suping,CHEN Fangfang
DOI: https://doi.org/10.13362/j.jpmed.202306014
2023-01-01
Abstract:Objective To investigate the expression of COL6A1 in sarcoma tissue and its significance based on bioinformatics analysis. Methods The GEO database was used to download the GSE16088, GSE21122, and GSE49972 datasets to analyze the expression level of COL6A1-3. The Kaplan-Meier survival curves were used to analyze the impact of the expression level of COL6A1-3 on the prognosis of patients with sarcoma. A Spearman correlation analysis was performed using GSVA package to investigate the correlation between the expression level of COL6A1 and immune cell infiltration. The LinkedOmics database was used to analyze the genes co-expressed with COL6A1, and the ggstatsplot package of R software was used to investigate the correlation between COL6A1 and COL6A2. The ClusterProfiler package of R software was used to perform GO functional enrichment analysis and KEGG pathway analysis of the genes co-expressed with COL6A1. The GSEA analysis was used to investigate the important kinase, miRNA, and transcription factor targets of the genes co-expressed with COL6A1. Results Compared with normal tissue, there was a significant increase in the expression level of COL6A1 in osteosarcoma, soft tissue sarcoma, and clear cell sarcoma of the kidney. The survival analysis showed that the high expression of COL6A1 was associated with low overall survival rate and di-sease-specific survival rate (χ2=6.37,4.57,P<0.05). The immune cell infiltration analysis showed a significantly positive correlation between COL6A1 and the level of neutrophil infiltration (r=0.280,P<0.05). The co-expression analysis showed that 2 444 genes were significantly positively correlated with the expression level of COL6A1, while 3 711 genes were negatively correlated with the expression level of COL6A1 (FDR<0.05). The correlation analysis also showed a significantly positive correlation between COL6A1 and COL6A2 (r=0.830,P<0.05). The GO functional enrichment analysis showed that the genes co-expressed with COL6A1 in sarcoma were mainly involved in the biological processes including immune cell infiltration. The GSEA analysis showed that there were five most important kinase targets that interact with the genes co-expressed with COL6A1, i.e., PRKCA, SRC, LYN, PRKAA1, and PAK1; the important miRNA targets were miR-506, miR-200B, miR-200C, miR-429, and miR-330, and the important transcription factor targets were V\$ETS2_B, V\$AP1_C, V\$PU1_Q6, V\$RORA1_01, and V\$AP1_Q2. Conclusion In conclusion, the COL6A1 gene is highly expressed in sarcoma tissue, which is closely associated with the poor prognosis of patients with sarcoma. The COL6A1 gene might serve as a valuable prognostic indicator and a potential therapeutic target for sarcoma.