首页|基于细胞焦亡基因前列腺癌预后模型的构建及肿瘤微环境分析

基于细胞焦亡基因前列腺癌预后模型的构建及肿瘤微环境分析

扫码查看
目的 基于细胞焦亡基因构建前列腺癌(PCa)预后模型并分析其肿瘤微环境.方法 从癌症基因组图谱(TCGA)数据库中获取502例PCa患者和51例健康者的基因转录组数据、基因突变数据,提取差异表达的细胞焦亡基因.从基因表达综合(GEO)数据库中获取96例PCa患者的临床特征及细胞焦亡基因表达数据,筛选出与预后相关的细胞焦亡基因.从ImmPort数据库中下载免疫细胞数据.采用聚类分析确定PCa患者的最佳亚型分组并采用主成分分析(PCA)进行验证.对不同亚型PCa患者的差异交集基因进行基因本位(GO)功能注释和京都基因与基因组百科全书(KEGG)富集分析.将单因素Cox分析获得的显著差异表达基因纳入Lasso回归模型,获取与预后相关的基因.将PCa患者随机均等分为训练集和测试集,分别进行Lasso回归分析,根据训练集风险评分的中位值将训练集、测试集分为高风险组和低风险组.采用多因素Cox回归模型分析PCa患者预后的独立影响因素并构建预后模型.应用R语言包分析细胞焦亡基因与免疫细胞含量、干细胞含量与风险评分的相关性及高低风险组患者的肿瘤突变负荷差异.结果 差异分析获得37个显著差异表达的细胞焦亡基因.根据细胞焦亡基因的表达量及临床数据分析,获得10个与PCa预后相关的细胞焦亡基因.根据累积分布函数图的K值将PCa分为A、B、C三个亚型,且经PCA分析验证分型结果可靠.GO及KEGG富集分析结果显示,109个差异交集基因主要通过细胞外基质组织(EMO)、细胞外结构组织(ESO)、肽基酪氨酸磷酸化(PTP)、肽基酪氨酸修饰(PTM)等通路激活蛋白酪氨酸激酶(PTK)、跨膜受体蛋白激酶(TRPK)等分子功能.Lasso回归模型筛选出B细胞淋巴瘤/白血病3(BCL3)、骨形成蛋白2(BMP2)、C-C趋化因子受体5(CCR5)、几丁质酶3样蛋白2(CHI3L2)、肝癌缺失蛋白1(DLC1)、主要组织相容性复合体Ⅱ类DQ alpha2(HLA-DQA2)、分泌球蛋白家族3A成员1(SC-GB3Al)、serpin家族E成员1(SERPINE1)、溶质载体家族14成员1(SLC14A1)9个与预后相关的基因.差异分析显示,训练集与测试集预后情况大致相同,且随着风险评分的升高,患者死亡人数增加.多因素Cox分析结果显示,年龄、肿瘤分期(T分期、N分期)、风险程度均是PCa患者预后的独立影响因素.预后模型预测PCa患者第6、9、12年的生存情况与患者的实际生存情况比较相符.相关性分析显示,免疫细胞与细胞焦亡基因存在相关性(P<0.01).对高风险组(n=376)与低风险组(n=126)进行差异分析,结果显示,高风险组患者的基质细胞评分和肿瘤突变负荷均明显高于低风险组,差异均有统计学意义(P<0.01).干细胞含量与风险评分呈正相关(P<0.01).结论 本文成功构建了基于细胞焦亡基因的PCa预后模型,并进行了肿瘤微环境分析,发现了 BCL3、BMP2、CCR5、CHI3L2、DLC1、HLA-DQA2、SCGB3A1、SERPINE1、SLC14A 1 9 个细胞焦亡基因,可预测 PCa 患者的预后和免疫治疗反应,有利于进一步指导临床实践.
Construction of a prognostic model of prostate cancer based on pyroptosis gene and analysis of tumor microenvironment
Objective To construct a prognostic model of prostate cancer(PCa)based on pyroptosis gene and ana-lyze the tumor microenvironment.Method The transcriptome and mutation data of 502 PCa patients and 51 healthy sub-jects were obtained from The Cancer Genome Atlas(TCGA)database,and the differentially expressed pyroptosis genes were extracted.The clinical features and expression data of pyroptosis gene in 96 PCa patients were obtained from Gene Expression Omnibus(GEO)database,and the prognostic-related pyroptosis genes were screened out.The immune cell da-ta were downloaded from the ImmPort database.The optimal subtypes of PCa samples were determined by cluster analy-sis and verified by principal component analysis(PCA).The differential intersection genes of different subtypes of PCa patients were annotated by Gene Ontology(GO)function and enriched by Kyoto Encyclopedia of Genes and Genomes(KEGG).The significantly differentially expressed genes obtained by univariate Cox analysis were incorporated into Las-so regression model to obtain the genes associated with prognosis.PCa patients were randomly and equally divided into training set and test set and Lasso regression analysis was performed respectively.The training set and test set were divid-ed into high risk group and low risk group according to the median value of the risk score of the training set.Multivariate Cox regression model was used to analyze the independent prognostic factors of PCa patients and to construct a prognos-tic model.R language package was used to analyze the correlation between pyroptosis genes and immune cell content,stem cell content and risk score and the difference of tumor mutation load in high and low risk groups.Result The differ-ential analysis obtained 37 significantly differentially expressed pyroptosis genes.According to the expression of pyropto-sis genes and clinical data analysis,10 genes related to PCa prognosis were obtained.According to the K value of cumula-tive distribution function(CDF)diagram,PCa were divided into three subtypes A,B and C,and the classification results were verified by PCA analysis to be reliable.GO and KEGG enrichment analysis showed that the 109 differentially inter-secting genes mainly activated protein tyrosine kinase(PTK)and transmembrane receptor protein kinase(TRPK)and oth-er molecular functions through extracellular matrix organization(EMO),extracellular structure organization(ESO),pepti-dyl-tyrosine phosphorylation(PTP),peptidyl-tyrosine modification(PTM)and other pathways.Lasso regression model screened out 9 prognostic related genes including B-cell leukemia/lymphoma 3(BCL3),bone morphogenetic protein 2(BMP2),C-C motif chemokine receptor 5(CCR5),chitinase 3 like 2(CHI3L2),deleted in liver cancer 1(DLC1),major histocompatibility complex,class Ⅱ,DQ alpha 2(HLA-DQA2),secretoglobin family 3A member 1(SCGB3A1),serpin family E member 1(SERPINE1),solute carrier family 14 member 1(SLC14A1).Analysis of differences showed that out-comes in the training and test sets were roughly the same,and that the number of patients who died increased with the in-crease in risk scores.Multivariate Cox analysis showed that age,tumor stage(T stage,N stage)and risk degree were inde-pendent factors influencing the prognosis of PCa patients.The prognosis model predicted the survival of PCa patients at the 6th,9th and 12th year,which was consistent with the actual survival of patients.Correlation analysis showed that there was a correlation between immune cells and pyroptosis genes(P<0.01).The difference between high-risk group(n=376)and low-risk group(n=126)was analyzed,and the results showed that the stromal cell score and tumor mutation load of high-risk group were significantly higher than those of low-risk group,and the differences were statistically significant(P<0.01).Stem cell content was positively correlated with risk score(P<0.01).Conclusion The PCa prognostic model based on pyroptosis genes is successfully constructed in this study,which finds that 9 pyroptosis genes,including BCL3,BMP2,CCR5,CHI3L2,DLC1,HLA-DQA2,SCGB3A1,SERPINE1,SLC14A1,can be used to predict the prognosis and immunotherapy response of PCa patients.It is helpful to further guide clinical practice.

pyroptosis geneprostate cancerprognostic modeltumor microenvironment

邵波、万水、田申、马园园、陈丹霞、杜沂宸

展开 >

昭通市中医医院泌尿外科,云南 昭通 657000

芜湖市中医医院泌尿外科,安徽 芜湖 241000

芜湖市中医医院皮肤科,安徽 芜湖 241000

芜湖市中医医院肛肠科,安徽 芜湖 241000

展开 >

细胞焦亡基因 前列腺癌 预后模型 肿瘤微环境

芜湖市科技计划芜湖市科技计划云南中医药大学校院联合基金

2021cg072022cg37XYLH202352

2024

癌症进展
中国医学科学院,北京协和医学院

癌症进展

影响因子:1.004
ISSN:1672-1535
年,卷(期):2024.22(2)