留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于外周血标志物初步探讨irAEs预测模型及价值

邓俊 王均 王茜 高嫦娥 陈晓 史明霞

邓俊, 王均, 王茜, 高嫦娥, 陈晓, 史明霞. 基于外周血标志物初步探讨irAEs预测模型及价值[J]. 昆明医科大学学报.
引用本文: 邓俊, 王均, 王茜, 高嫦娥, 陈晓, 史明霞. 基于外周血标志物初步探讨irAEs预测模型及价值[J]. 昆明医科大学学报.
Zhang Dong Xian , Hong Shi Jun , Zhao Li Ping , He Yu Ting , Zhang Yun Long , Hu Yun Bo , Du Qi , Guan Li Na , Li Li Hua . Related Analysis on Disability after Craniocerebral Injury Caused by Road Traffic Accident[J]. Journal of Kunming Medical University, 2016, 37(11): 134-138.
Citation: Jun DENG, Jun WANG, Xi WANG, Change GAO, Xiao CHEN, Mingxia SHI. Prediction Model and Its Value of IrAEs Based on Peripheral Blood Markers[J]. Journal of Kunming Medical University.

基于外周血标志物初步探讨irAEs预测模型及价值

基金项目: 202201AY070001-050(HSP110DE9在MSS型结直肠癌中抗PD-1免疫治疗的作用及机制研究)
详细信息
    作者简介:

    邓俊(1986~),女,四川苍溪人,医学硕士,主治医师,主要从事恶性肿瘤免疫治疗研究工作

    通讯作者:

    史明霞,E-mail:shmxia2002@sina.com

  • 中图分类号: R730.3

Prediction Model and Its Value of IrAEs Based on Peripheral Blood Markers

  • 摘要:   目的  基于外周血标志物探讨irAEs预测模型及价值。  方法  回顾性收集2020年12月至2023年12月昆明医科大学第一附属医院就诊且使用PD-1/PD-L1抗体治疗的825例恶性肿瘤患者的基线临床资料、实验室检查、irAEs随访结果,根据是否存在irAEs分为irAEs组和non-irAEs组,组间及组内的差异性分析采用t检验、秩和检验、卡方检验、Fisher确切概率法;运用LASSO、Ridge、Elastic-net logistic回归筛选预测因子并建立irAEs风险预测模型。   结果  136例患者经历178次irAEs,其中主要为内分泌毒性占42.64%,肝炎35.29%,肺炎20.58%,≥ G3级占19.07%,累及两种以上器官占总irAEs人数的24.26%。单因素分析结果显示,基线CD4+ T细胞计数、IL-6、IL-17、TSH、GLB和ALB与irAEs存在一定关联;通过Ridge、LASSO和Elastic-Net Logistic回归模型筛选出GLB、ALB、IL-17、TSH为重要风险因素,结果显示三类算法AUC均超过0.800。内部验证集LASSO-Logistic AUC为0.800(95%CI = 0.739~0.862)。外部验证集AUC为0.800( 95%CI0.739~0.861),且DCA曲线结果提示该预测模型的净收益率最高。  结论  GLB、ALB、IL-17、TSH是irAEs的独立预测因子,以它们为基础的irAEs预测模型预测效能良好。
  • 肺癌是全球最常见的癌症,占所有肿瘤例数的11.6%[1]。根据最新Global cancer statistics统计,每年估计诊断出220万新肺癌病例,约180万例肺癌死亡,造成了重大的社会负担和经济损失[2]。肺癌的早期治疗策略优先考虑解剖性肺切除联合淋巴结清扫。然而,术后治愈率和生存率仍然很低,30%~70%的术后患者出现肿瘤复发或转移,治疗效果非常不理想[3]。因此,确定新型和有效的生物标志物和治疗靶点对于肺癌防治是至关重要。

    作为一种重要的细胞内矿物营养物质,铜在几个重要的细胞过程中发挥作用,包括线粒体呼吸、氧代谢、铁摄取以及一些生物途径的调节[45]。此前的研究已经观察到铜与疾病状态之间的联系,例如,在各种恶性肿瘤中发现了更高水平的铜[45]。最近,Tsvetkov等[6]揭示了一种称为“铜死亡”的新机制,即过高的细内铜浓度会导致细胞死亡。值得注意地。已有研究证实铜死亡参与肺癌的进程。例如,PSMD11被鉴定为铜死亡和免疫相关基因,并且促进肺癌细胞的增殖和转移[7]。此外,AC144450.1/miR-424-5p/CBX2轴被鉴定为铜死亡相关的竞争性内源RNA机制,并且该轴的激活有利于肺癌细胞的生长[8]。这为肺癌治疗干预组合提供了潜在的方向。然而,铜死亡在肺癌中的研究仍然处于初步阶段。

    本研究拟通过生物信息学鉴定与肺癌预后和免疫微环境相关的铜死亡基因。旨在为肺癌的治疗提供新型的生物标志物和铜死亡相关的理论基础。

    本研究使用的肺癌数据集来源于TCGA(https://portal.gdc.cancer.gov/),包括肺腺癌(lung adenocarcinom,TCGA-LUAD)和肺鳞细胞癌(lung squamous cell carcinomas,TCGA-LUSC)队列。从TCGA数据库提取队列的count和fpkm格式数据。为避免模型构建受到不同队列研究间批次效应的影响,log2转换的fpkm(log2-fpkm)值采用sva包进行批次效应校正,校正中将样本类型(肿瘤或正常组织)作为协变量,避免批次校正掩盖组织类型间生物学效应的差异。最后筛选患者生存信息完整的肿瘤样本数据及所有正常组织样本数据作为分析使用的数据集,其中包括616例肿瘤样本和108例正常组织样本。Count数据用于差异分析,log2-fpkm数据用于预后模型构建。从此前的文献中共检索到38个CRGs[9],包括ATP7B、CDKN2A、DLD、DPYD、FDX1、GLRX5、GLS、ISCA2、LIPT1、MTF1、NDUFA1、NDUFA8、NDUFB10、NDUFB2、NDUFB6、NDUFC1、NDUFC2、NDUFV2、PDHA1、PLAT、POLD1、PPAT、SLC31A1、SDHB、TIMMDC1、DLAT、GCSH、DBT、DLST、LIAS、LIPM、LIPA、LIPT2、PDHB、ACO2、NLRP3、ATP7A和NFE2L2。

    采用R软件Deseq2包[10]对TCGA肺癌队列count格式数据执行差异表达分析,以鉴定差异表达基因(differentially expressed genes,DEGs)。阈值为|log2-fold change (FC)|≥2且校准的P-value<0.05。采用R软件ggplot2包[11]可视化DEGs的火山图。采用R软件VennDiagram包[12]获取和可视化DEGs和CRGs的交集,以鉴定DE-CRGs。

    DE-CRGs的富集分析包括GO和KEGG富集分析,以识别DE-CRGs的生物学功能和可能涉及的信号通路。DE-CRGs采用R软件clusterProfiler包[13]进行GO和KEGG富集分析,并采用R软件ggplot2包[11]进行气泡图的可视化。

    采用STRING数据库预测DE-CRGs的互作关系。最低要求的互作分数设置为中等可行度(0.4)。DE-CRGs互作关系的TSV从STRING数据库导出,并且导入Cytoscape软件。通过Cytoscape软件的Network Analysis进行PPI网络的可视化,并统计每个节点的连接度、介数中心性和接近中心性。

    采用R软件Survival包[14]进行TCGA肺癌队列log2-fpkm格式数据(共616个肺癌样本,包括344例死亡和272例存活)的Log-rank检验和Cox回归分析,以识别DE-CRGs表达和临床特征(年龄、性别、是否吸烟、TNM分期、残留肿瘤情况)与肺癌患者预后的相关性,筛选可以独立用于预后风险预测的DE-CRGs。计算每个变量的风险比(hazard ratio,HR)和95%置信区间(confidence interval,CI)。P-value<0.05被认为具有统计学意义。采用R软件survminer包进行森林图和K-M曲线的可视化。将上述筛选出的独立临床因素用于预测列线图的构建。使用R软件rms包[15]构建预测列线图和相应的校准曲线。校准曲线越接近代表最佳预测的45°线,代表列线图的预后预测性能越好。

    本研究采用CIBERSORT和ESTIMATE算法评估LIPT1表达与肺癌患者免疫微环境的相关性。CIBERSORT算法是基于线性支持向量回归原理对人类免疫细胞亚型的表达矩阵进行去卷积的一种算法[16]。本研究通过CIBERSORT算法分析了肺癌队列中22种人类免疫细胞亚群的比例,及LIPT1表达与22种免疫细胞浸润的相关性。Estimate算法通过使用表达数据预测恶性肿瘤组织中的基质细胞和免疫细胞存在的工具[17]。本研究通过ESTIMATE算法评估肿瘤纯度,并分析LIPT1与免疫分数和基质分数之间的相关性。

    采用R软件(版本4.1.1)进行统计分析。Student's t检验和Mann-Whitney检验用于比较肺癌组织和正常组织之间的基因表达。采用Kaplan-Meier(K-M)曲线及Log-rank统计检验对高表达LIPT1和低表达LIPT1的患者生存进行差异分析。使用单变量和多变量Cox回归分析来确定肺癌的独立的预后因素。P<0.05为差异有统计学意义。

    R软件Deseq2包分析显示,相较于正常组织,TCGA-LUAD和TCGA-LUSC队列的肺癌组织中共存在5269个DEGs,包含1130个下调的DEGs和4139个上调的DEGs,见图1A。韦恩图显示,DEGs和CRGs共存在11个交集,包括POLD1、DLD、ACO2、DLST、NDUFA1、DBT、NDUFB10、LIPT1、DLAT、NDUFB6、GLRX5,见图1B。这些交集为DE-CRGs。

    图  1  肺癌中DE-CRGs的鉴定
    A:火山图展示了TCGA-LUAD和TCGA-LUSC队列中DEGs的表达和差异性。每个点代表一个基因;红色和蓝色的点分别代表上调和下调DEGs;灰色的点代表表达没有显著差异性的基因;B:韦恩图展示了DEGs和CRGs的交集和并集数量。
    Figure  1.  Identification of DE-CRGs in lung cancer

    GO富集分析显示,11个DE-CRGs主要参与有氧呼吸、硫化合物代谢过程、氨基酸分解代谢过程、三羧酸循环和ATP代谢过程的调控,见图2A。KEGG富集分析显示,11个DE-CRGs主要调控碳代谢、氧化磷酸化、氨基酸代谢、乙醛酸盐和二羧酸盐代谢和柠檬酸循环,见图2B

    图  2  肺癌中DE-CRGs的富集分析
    A~B:气泡图展示了DE-CRGs的(A)GO和(B)KEGG富集分析的结果。点的大小代表富集到对应术语中的DE-CRGs数量;点的颜色代表校准的P值的大小,颜色越蓝代表校准的P值越小。
    Figure  2.  Enrichment analysis of DE-CRGs in lung cancer

    DE-CRGs相关PPI网络包含10个节点和19个互作关系对,其中不包含POLD1,见图3A。该网络中,DLST、ACO2、DLD、DLAT、DBT是连接度排名前五的DE-CRGs,见图3B

    图  3  肺癌中DE-CRGs相关PPI网络的构建
    A:基于STING数据库和Cytoscape软件构建的DE-CRGs相关PPI网络。节点越大代表连接度越大;B:表格展示了每个节点的连接度、介数中心性和接近中心性。
    Figure  3.  Construction of DE-CRGS-related PPI networks in lung cancer

    单因素COX分析表明,11个DE-CRGs中仅有DBT和LIPT1与肺癌患者生存相关,见图4A。进一步的多因素COX分析表明,LIPT1、年龄、TNM分期和残留肿瘤是肺癌患者的独立风险因素,见图4A。Log-rank检验同样表明,相较于高表达LIPT1的肺癌患者,低表达LIPT1的肺癌患者的生存可能更低,见图4B。本研究构建了预测肺癌患者预后的独立风险因素相关列线图,见图4C。校准曲线显示,在第1年、第2年和第3年,列线图预测的患者生存可能与患者实际的生存/结局相近,见图4D

    图  4  与肺癌患者预后相关的DE-CRGs鉴定
    A:单因素和多因素COX分析识别与肺癌患者生存相关的DE-CRGs和临床特征;B:Log-rank检验鉴定LIPT1表达与患者生存可能的相关性;C:基于LIPT1、年龄、TNM分期和残留肿瘤构建的列线图,以预测肺癌患者的生存可能;D:列线图的校准曲线,以观察列线图的预后预测性能。
    Figure  4.  Identification of DE-CRGs associated with prognosis in patients with lung cancer

    ESTIMATE算法分析显示,LIPT1表达与肺癌患者的免疫分数、基质分数和ESTIMATE分数呈现显著的负相关,见图5A。CIBERSORT算法分析显示,LIPT1表达与T细胞滤泡辅助细胞、静息肥大细胞、单核细胞、嗜酸性粒细胞、树突细胞、M1和M2型巨噬细胞、γδ T细胞、活化的自然杀伤细胞细胞、CD8+ T细胞的浸润呈现正相关,与M0型巨噬细胞、激活的肥大细胞、浆细胞、中性粒细胞、静息的自然杀伤细胞和Treg细胞显现显著的负相关,见图5B

    图  5  LIPT1与肺癌患者免疫微环境的相关性分析
    A:ESTIMATE算法识别LIPT1表达与免疫分数、基质分数和Estimate分数的相关性;B:CIBERSORT算法鉴定LIPT1表达与22种免疫细胞浸润的相关性。
    Figure  5.  Correlation analysis of LIPT1 and immune microenvironment in patients with lung cancer

    本研究在肺癌中共鉴定了11个DE-CRGs,包括ACO2、POLD1、DLD、DLST、NDUFA1、DBT、NDUFB10、LIPT1、DLAT、NDUFB6、GLRX5。有趣地,部分DE-CRGs在肺癌进程中的功能已被证实。Shideh等[18]证实,ACO2在非小细胞肺癌组织中低表达,这预示这患者预后变差,并且ACO2通过抑制CISD1导致铁饥饿反应来降低FeS转运,这抑制了癌细胞的增殖。此外,PM2.5通过翻译和转录激活DLAT介导的糖酵解重编程,这加速了非小细胞肺癌的生长[19]。然而,这些研究均不涉及铜死亡的调控。结合本研究的发现提示,ACO2和DLAT对肺癌的抑癌和促癌作用可能与铜死亡有关。值得注意地,富集分析证实,11个DE-CRGs主要调控能量代谢。这是可以预见的。因为,过量的铜会损害线粒体呼吸中的三羧酸循环,并刺激铜中毒,导致脂酰化和FeS簇蛋白的明显丢失以及热休克蛋白70水平的增加,进而引起细胞的程序性死亡[6]

    本研究证实,LIPT1低表达是肺癌的独立风险因子,并且LIPT1相关的列线图对肺癌患者预后具有卓越的预测性能。LIPT1是编码线粒体硫辛酸代谢的重要基因,其编码的LIPT1蛋白将硫辛酰基从硫辛酰腺苷酸转移到甘氨酸裂解系统蛋白H和2-氧酸脱氢酶E2亚基,进而参与硫辛酸代谢[620]。目前,有关于LIPT1的研究主要集中于Leigh综合征[21]、硫辛酸代谢相关的致命疾病[22]、恶性肿瘤(膀胱癌、黑色素瘤和肝细胞癌)[2325]。在膀胱癌和黑色素瘤中,LIPT1的高表达预示着患者预后良好,这与本研究的结果是一致的。然而,在肝细胞癌中,LIPT1表现为高表达,并且是肝细胞癌预后不良的独立风险因素。这表明,LIPT1在不同细胞类型中以完全相反的方式发挥作用,这一点对于理解癌症中的细胞命运决定至关重要。然而,LIPT1在肺癌发生和发展中的功能是未知的,尤其是对铜死亡的调控。值得注意地,在黑色素瘤患者中,上调的LIPT1表达可能会抑制Treg的浸润,从而提高免疫治疗效果[24]。本研究发现,LIPT1表达与Treg的浸润呈现显著的负相关。这个发现提示,LIPT1可能也通过调控Treg的浸润,参与肺癌免疫治疗疗效的干预。此外,本研究还发现,LIPT1表达与T细胞滤泡辅助细胞、肥大细胞、树突细胞、巨噬、γδ T细胞、自然杀伤细胞细胞、CD8+ T细胞和中性粒细胞的浸润有关。这些发现表明,LIPT1可能是肺癌免疫治疗的有效靶点,为肿瘤患者的临床治疗提供了新的希望。然而,LIPT1表达与癌症患者免疫检查点之间的关联仍需在更多的临床前和临床试验中探索。

    综上所述,本研究利用生物信息学分析探讨了与肺癌患者预后和免疫微环境相关的铜死亡基因。本研究证实,LIPT1是一种有效预测肺癌患者预后,并且与肿瘤免疫微环境相关的铜死亡生物标志物。本研究为进一步研究LIPT1在肺癌发生和治疗中的具体机制奠定了基础。

  • 图  1  irAEs特征

    A:类型,B:分级,C:占比率%,D:中位时间。

    Figure  1.  Characteristics of irAEsA

    图  2  筛选流程图

    Figure  2.  Flow chart of screening

    图  3  irAEs特征(外部验证)

    A:类型,B:分级,C:占比率%,D:中位时间。

    Figure  3.  Characteristics of irAEs(External validation )

    图  4  路径筛选系数

    A:二项式偏差;B:路径系数。

    Figure  4.  Path filtering coefficient

    图  5  ROC曲线

    A:三种算法训练集的ROC曲线;B:三种算法验证集的ROC曲线。

    Figure  5.  ROC curves A:ROC curves of training sets by three algorithms,B:ROC curves of three algorithm for validation sets

    图  6  矫正曲线及决策曲线

    A:训练集DCA矫正曲线;B:内部测试集DCA矫正曲线;C:训练集DCA决策曲线;D:内部测试集DCA决策曲线。

    Figure  6.  DCA correction curve and DCA decision curve

    图  7  列线图及外部测试ROC曲线

    A:风险预测模型列线图;B:外部测试集ROC曲线。

    Figure  7.  Nomogram of risk prediction model and ROC curve of external test set

    表  1  训练集、内部测试集和外部测试集人口学及治疗特征的描述性与可比性分析

    Table  1.   Descriptive and comparable analysis of demographic and treatment characteristics for training set,internal test set and external test set

    变量 训练集 内部测试集 外部测试集
    非irAE(n = 199)/
    irAE(n = 88)
    统计量 P 非irAE(n = 74)/
    irAE(n = 48)
    统计量 P 非irAE(n = 224)/
    irAE(n = 192)
    统计量 P
    性别 χ2 = 0.21 0.642 χ2 = 2.14 0.144 χ2 = 5.19 0.023
     男 144(68.57)/66(31.43) 55(57.29)/41(42.71) 176(57.14)/132(42.86)
     女 55(71.43)/22(28.57) 19(73.08)/7(26.92) 48(44.44)/60(55.56)
    年龄 χ2 = 0.01 0.966 χ2 = 0.06 0.812 χ2 = 1.40 0.237
     < 60 48(30.77)/108(69.23) 27(40.30)/40(59.70) 132(51.56)/124(48.44)
     ≥ 60 40(30.53)/91(69.47) 21(38.18)/34(61.82) 92(57.50)/68(42.50)
    BMI χ2 = 5.38 0.069 χ2 = 0.35 0.839 χ2 = 8.03 0.018
     肥胖 52(59.77)/35(40.23) 22(64.71)/12(35.29) 64(44.44)/80(55.56)
     偏瘦 31(73.81)/11(26.19) 8(57.14)/6(42.86) 20(62.50)/12(37.50)
     正常 116(73.42)/42(26.58) 44(59.46)/30(40.54) 140(58.33)/100(41.67)
    ECOG评分 χ2 = 8.39 0.003 χ2 = 0.05 0.832 χ2 = 0.10 0.749
     1 164(66.13)/84(33.87) 67(60.36)/44(39.64) 184(53.49)/160(46.51)
     2 35(89.74)/4(10.26) 7(63.64)/4(36.36) 40(55.56)/32(44.44)
    瘤种 χ2 = 9.13 0.058 χ2 = 5.68 0.224 χ2 = 13.91 0.008
     肺癌 90(65.69)/47(34.31) 41(55.41)/33(44.59) 100(52.08)/92(47.92)
     肝癌 20(80.00)/5(20.00) 7(63.64)/4(36.36) 20(45.45)/24(54.55)
     其他a 41(80.39)/10(19.61) 12(80.00)/3(20.00) 40(52.63)/36(47.37)
     食管癌 23(76.67)/7(23.33) 4(44.44)/5(55.56) 40(76.92)/12(23.08)
     头颈癌 25(56.82)/19(43.18) 10(76.92)/3(23.08) 24(46.15)/28(53.85)
    肿瘤分期 χ2 = 5.37 0.146 - 0.177 - 0.001
     Ⅰ期 7(77.78)/2(22.22) 0(0.00)/1(100.00) 0(0.00)/4(100.00)
     Ⅱ期 51(79.69)/13(20.31) 25(71.43)/10(28.57) 8(33.33)/16(66.67)
     Ⅲ期 58(69.05)/26(30.95) 21(63.64)/12(36.36) 92(67.50)/68(42.50)
     Ⅳ期 83(63.85)/47(36.15) 28(52.83)/25(47.17) 124(54.38)/104(46.62)
    治疗线数 χ2 = 5.70 0.116 χ2 = 0.67 0.735 χ2 = 2.09 0.351
     一线 128(73.56)/46(26.44) 44(57.14)/33(42.86) 144(54.55)/120(45.45)
     二线 34(58.62)/24(41.38) 12(63.16)/7(36.84) 76(54.29)/64(45.71)
     三线 11(57.89)/8(42.11) 3(75.00)/1(25.00) 4(33.33)/8(66.67)
    用药模式 χ2 = 12.08 0.001 χ2 = 5.41 0.004 - 0.088
     单药 2(2.02)/30(10.45) 11(12.94)/1(2.70) 36(60.00)/24(40.00)
     联合化疗 45(45.45)/152(52.96) 46(54.12)/22(59.46) 128(52.46)/116(47.54)
     联合化疗+靶b 20(20.20)/36(12.54) 12(14.12)/5(13.51) 56(54.83)/48(45.17)
     联合化疗
     +放疗
    16(16.16)/36(12.54) 6(7.06)/6(16.22) 0(0.00)/4(100.00)
     其他c 16(16.16)/33(11.50) 10(11.76)/3(8.11) 4(100.00)/0(0.00)
      注:a:代表胃恶性肿瘤、胰腺恶性肿瘤、胸腺恶性肿瘤、恶性黑色素瘤、胆囊恶性肿瘤、结肠恶性肿瘤、妇科恶性肿瘤、恶性淋巴瘤、肉瘤、输尿管恶性肿瘤;b:靶向药物均为抗血管生成药物,如贝伐珠单抗、仑伐替尼、安罗替尼等;c.其他包括免疫联合化疗+放疗+靶向、免疫联合放疗、免疫联合介入、免疫联合免疫等;-:代表Fisher确切概率法。
    下载: 导出CSV

    表  2  训练集、内部测试集和外部测试集实验室数据与可比性分析

    变量 训练集 内部测试集 外部测试集
    非irAE(n = 199)/
    irAE(n = 88)
    统计量 P 非irAE(n = 74)
    /irAE(n = 48)
    统计量 P 非irAE(n = 224)/
    irAE(n = 192)
    P
    CD4+T 659.59±167.41/
    531.41±167.89
    t =
    8.22
    < 0.001 662.72±184.21/
    492.15±162.99
    t =
    6.35
    < 0.001 517.00[443.00,591.00]/
    518.00[387.00,591.00]
    Z = −0.91 0.363
    CD8+T 435.65±141.20/
    475.78±150.93
    t = −2.99 0.003 435.56±161.47/
    438.29±154.42
    t = −0.11 0.909 259.00[198.00,369.00]/
    316.00[226.50,391.00]
    Z =
    2.78
    0.006
    CRP 23.27[16.74,40.23]/
    23.27[13.56,44.63]
    Z = 0.38 0.707 23.27[19.99,44.88]/
    23.27[15.24,48.82]
    Z = 0.36 0.722 1.69[0.69,2.36]/
    3.18[0.45,5.61]
    Z =
    4.60
    < 0.001
    IL-5 3.87[2.56,5.54]/
    4.53[2.92,6.04]
    Z = 2.55 0.011 3.67[2.63,5.38]/
    4.88[3.67,6.42]
    Z = 3.16 0.002 2.62[1.78,4.24]/
    2.24[1.43,3.86]
    Z = −2.56 0.011
    IL-6 18.15±12.99/
    21.85±21.95
    Z = 1.28 0.202 16.17[15.95,19.52]/
    17.80[12.14,35.43]
    Z = 2.13 0.033 6.50[4.88,8.86]/
    3.22[2.07,4.78]
    Z = −11.72 < 0.001
    IL-10 2.43±3.10/
    2.65±1.14
    t = −0.89 0.371 2.00±0.95/
    2.88±1.10
    t = −5.92 < 0.001 2.57[1.63,4.59]/
    2.15[1.19,3.59]
    Z = −3.23 0.001
    IL-17 1.95[1.35,2.80]/
    4.27[3.36,5.03]
    Z = 12.92 < 0.001 1.88[1.31,2.75]/
    4.34[2.95,5.18]
    Z = 7.39 < 0.001 2.34[1.57,2.69]/
    6.05[4.44,7.47]
    Z = 12.61 < 0.001
    LDH 227.46±66.41/
    245.97±66.55
    t = −2.99 0.003 244.09±73.02/
    233.57±69.19
    t =
    0.97
    0.332 172.50[156.50,219.50]/
    184.50[161.50,246.50]
    Z =
    2.84
    0.005
    TSH 2.36[1.39,4.35]/
    4.40[3.07,6.26]
    Z = 9.09 < 0.001 2.29[1.58,3.94]/
    4.18[2.88,6.03]
    Z = 5.22 < 0.001 1.40[1.05,2.51]/
    2.18[1.28,3.01]
    Z =
    3.44
    < 0.001
    白细胞 6.37[5.20,8.02]/
    6.08[4.86,7.41]
    Z = −1.67 0.095 6.21[5.25,7.79]/
    6.06[4.91,8.29]
    Z = −0.69 0.559 5.76[4.61,7.37]/
    5.88[5.18,7.41]
    Z = −0.41 0.680
    中性粒
    细胞
    4.11[3.05,5.50]/
    3.82[2.75,4.81]
    Z = −1.60 0.11 3.89[3.04,5.34]/
    3.83[3.14,5.69]
    Z = −0.58 0.897 3.58[2.77,4.92]/
    3.92[2.73,4.63]
    Z =
    1.12
    0.263
    淋巴
    细胞
    1.48[1.11,1.99]/
    1.47[1.15,1.88]
    Z = −0.18 0.860 1.51[1.11,1.96]/
    1.42[1.03,1.78]
    Z = 0.13 0.110 1.60[1.05,2.09]/
    1.26[0.93,1.86]
    Z = −3.43 < 0.001
    血小板 220.00[164.00,289.00]/
    230.00[177.00,281.00]
    Z = 0.57 0.568 242.00[180.00,307.00]/
    220.50[178.50,278.50]
    Z = −0.56 0.374 234.00[187.00,265.50]/
    254.50[202.50,296.50]
    Z =
    3.21
    0.001
    NLR
    (N/L)
    2.73[1.77,4.27]/4
    .96[4.34,5.82]
    Z = 11.21 < 0.001 2.67[1.90,3.79]/
    5.60[4.46,6.97]
    Z = −0.89 < 0.001 2.51[1.81,3.20]/
    2.94[2.03,4.78]
    Z =
    4.20
    < 0.001
    PLR
    (P/L)
    154.64[106.49,210.79]/172.79[130.56,242.64] Z = 3.27 0.001 158.86[108.00,224.43]/
    190.03[146.23,274.81]
    Z = 8.56 0.005 148.16[103.15,197.10]/
    187.35[128.36,236.40]
    Z =
    4.52
    < 0.001
    球蛋白 34.52[29.71,39.45]/
    34.89[29.00,39.80]
    Z = −0.53 0.597 34.31±7.18/
    33.96±5.96
    Z = 2.81 0.738 28.40[26.25,32.35]/
    29.00[25.90,32.35]
    Z =
    1.64
    0.530
    白蛋白 39.80[36.50,43.80]/
    41.70[37.30,44.60]
    Z = 2.08 0.038 40.90[36.90,44.90]/
    41.40[36.90,43.85]
    Z = −0.02 0.981 42.05[38.65,43.95]/
    42.50[38.90,44.20]
    Z =
    0.29
    0.774
    PNI 46.35±15.67/
    41.72±6.60
    t =
    3.67
    < 0.001 47.10[42.05,51.45]/
    40.58[37.23,44.55]
    Z = −0.02 < 0.001 49.38[45.60,53.53]/
    48.18[44.80,51.40]
    Z = −1.98 0.048
    下载: 导出CSV

    表  3  惩罚项最小误差对应的log(λ)值

    Table  3.   The log (λ) value corresponding to the minimum error of the penalty term

    CategoryRidgeLASSOElastic-net
    lambda.min0.051551290.018100650.02071571
    lambda.1se0.20811340.031631390.06943083
    下载: 导出CSV

    表  4  LASSO-Logistic风险回归模型

    Table  4.   Risk regression model by LASSO-logistic

    变量 OR(95%CI Estimate SE Wald χ2 P
    截距 1.365(0.462,4.038) 0.311 0.553 0.560 0.574
    球蛋白 1.245(1.167,1.328) 0.219 0.033 6.630 < 0.001
    IL-17 1.090(1.061,1.120) 0.086 0.014 6.180 < 0.001
    白蛋白 1.399(1.249,1.566) 0.336 0.058 5.810 < 0.001
    TSH 1.546(1.305,1.831) 0.435 0.086 5.040 < 0.001
    下载: 导出CSV
  • [1] 张雨峤,梅伟健. 免疫检查点抑制剂治疗实体瘤的标志性成果[J]. 实用医学杂志,2024,40(2):272-277.
    [2] 吴皎,杨润祥. 免疫检查点抑制剂相关不良事件识别与管理[J]. 中国临床医生杂志,2023,51(10):1135-1137+1132. doi: 10.3969/j.issn.2095-8552.2023.10.001
    [3] 叶淑芳,刘苗兰,仇晓霞. 基于真实世界的研究—646例肿瘤患者使用免疫检查点抑制剂相关不良反应的临床特征[J]. 肿瘤药学,2024,14(1):72-79.
    [4] Les,I,Martínez,M,Pérez-Francisco,I,et al. Predictive Biomarkers for checkpoint inhibitor immune-related adverse events[J]. Cancers (Basel),2023,15(5):1629. doi: 10.3390/cancers15051629
    [5] 缪延栋,赵雪寒,张怡明,等. 免疫检查点抑制剂治疗恶性肿瘤所致免疫相关不良反应预测标志物的研究进展[J]. 中华肿瘤防治杂志,2024,31(6):383-390.
    [6] Egami,S,Kawazoe,H,Hashimoto,H,et al. Peripheral blood biomarkers predict immune-related adverse events in non-small cell lung cancer patients treated with pembrolizumab: A multicenter retrospective study[J]. J Cancer,2021,12(7):2105-2112. doi: 10.7150/jca.53242
    [7] 徐瑞华,李进等,免疫检查点抑制剂相关的毒性管理指南2023[M]. 北京: 人民卫生出版社,2023: 4-132.
    [8] Zheng,Y,Kim,R,Yu,T,et al. Real-World clinical and economic outcomes in selected immune-related adverse events among patients with cancer receiving immune checkpoint inhibitors[J]. Oncologist,2021,26(11):e2002-e2012. doi: 10.1002/onco.13918
    [9] Jing,Y,Yang,J,Johnson,DB,et al. Harnessing big data to characterize immune-related adverse events[J]. Natrev Clin Oncol,2022,19(4):269-280.
    [10] 徐瑞华,李进等,免疫检查点抑制剂临床应用指南2023[M]. 北京: 人民卫生出版社,2023: 4-132.
    [11] Shimozaki,K,Sukawa,Y,Sato,Y,et al. Analysis of risk factors for immune-related adverse events in various solid tumors using real-world data[J]. Future Oncol,2021,17(20):2593-2603. doi: 10.2217/fon-2020-0861
    [12] Luongo,C,Morra,R,Gambale,C,et al. Higher baseline TSH levels predict early hypothyroidism during cancer immunotherapy[J]. J Endocrinol Invest,2021,44(9):1927-1933. doi: 10.1007/s40618-021-01508-5
    [13] Johnson,D,Patel,AB,Uemura,MI,et al. IL17A blockade successfully treated psoriasiform dermatologic toxicity from immunotherapy[J]. Cancer Immunlo Res,2019,7(6):860-865.
    [14] Husain,B,Kirchberger,MC,Erdmann,M,et al. Inflammatory markers in autoimmunity induced by checkpoint inhibitors[J]. J CANCER RES CLIN,2021,147(6):1623-1630.
    [15] Esfahani,K,Miller,WH. Reversal of Autoimmune toxicity and loss of tumor response by interleukin-17 blockade[J]. New Engl J Med,2017,376(20):1989-1991.
    [16] KANG D H,PARK C K,CHUNG C,et al. Baseline serum interleukin-6 levels predict the response of patients with advanced non-small cell lung cancer to PD-1/PD-L1 inhibitors[J]. Immune Netw,2020,20(3):e27.
    [17] Dharmapuri,S,Ozbek,U,Saeed,A,et al. Relationship between systemic inflammatory response markers and immune treatment related toxicity (IrAEs) in hepatocellular carcinoma (HCC)[J]. J Clin Oncol,2022,40(16_suppl):e16204-e16204.
    [18] Nakamura,Y,Tanaka,R,Maruyama,H,et al. Correlation between blood cell count and outcome of melanoma patients treated with anti-PD-1 antibodies[J]. Jpn J Clin Oncol,2019,49(5):431-437.
    [19] Suazo-Zepeda,E,Bokern,M,Vinke,PC,et al. Risk factors for adverse events induced by immune checkpoint inhibitors in patients with non-small-cell lung cancer: A systematic review and meta-analysis[J]. CANCER IMMUNOL IMMUN,2021,70(11):3069-3080.
    [20] Zhang,Y,Zhao,J,Hu,W,et al. Construction of a near-term predictive model for irAEs induced by PD-1 inhibitors[J]. J Clin Oncol,2022,40(16_suppl):3147-3147.
    [21] 刘剑,李敏菁. 全身免疫炎症指数对非小细胞肺癌免疫检查点抑制剂疗效的预测价值[J]. 实用医学杂志,2022,38(7):904-908. doi: 10.3969/j.issn.1006-5725.2022.07.024
    [22] Egami,S,Kawazoe,H,Hashimoto,H,et al. Peripheral blood biomarkers predict 8immune-related adverse events in non-small cell lung cancer patients treated with pembrolizumab: A multicenter retrospective study[J]. Cancer,2021,12(7):2015-2112.
  • [1] 赵珊, 高赛, 李堂春, 赵钟鸣, 吴亚婷, 郑敏.  云南省孕产妇死亡率预测模型的构建与评估, 昆明医科大学学报. 2025, 46(2): 110-117. doi: 10.12259/j.issn.2095-610X.S20250216
    [2] 王远珍, 魏红艳, 田仁海, 陈永臻, 许丹青, 张映媛, 常丽仙, 刘春云, 刘立.  原发性肝癌干预前并发慢性肝衰竭的风险预测模型建立与评价, 昆明医科大学学报. 2025, 46(3): 1-9.
    [3] 刘东艳, 刘燕, 任治铭, 王枫, 汪勇.  基于Logistic回归分析的恶性肿瘤患者化疗相关骨髓抑制的影响因素分析及预测模型建立, 昆明医科大学学报. 2025, 46(4): 1-7.
    [4] 闫兰竹, 乔顺义, 张艳丽, 赵二强, 杨虎, 林静, 付海艳.  免疫球蛋白在肝癌TACE治疗过程中的变化及预测疗效的价值, 昆明医科大学学报. 2024, 45(4): 128-134. doi: 10.12259/j.issn.2095-610X.S20240418
    [5] 徐文秀, 莫小凤, 杨祥民, 杨丝露, 吴凡, 李特.  Logistic回归与决策树模型在碘造影剂不良反应预测中的应用, 昆明医科大学学报. 2024, 45(9): 70-75. doi: 10.12259/j.issn.2095-610X.S20240911
    [6] 姬燕梅, 李文俊, 李青芸, 郭妮, 蒙妮, 周丹, 李秋宇, 金醒昉.  急性缺血性脑卒中后认知障碍相关因素分析及列线图模型构建, 昆明医科大学学报. 2024, 45(5): 73-81. doi: 10.12259/j.issn.2095-610X.S20240511
    [7] 曾慧娟, 田波, 袁红伶, 何杰, 李冠羲, 茹国佳, 许敏, 詹东.  机器学习算法构建慢性肾脏病伴高血压或糖尿病的预测模型, 昆明医科大学学报. 2024, 45(3): 99-105. doi: 10.12259/j.issn.2095-610X.S20240315
    [8] 李梦, 牛奔, 胡瑾, 马燕琳.  免疫检查点抑制剂引起肿瘤患者甲状腺功能障碍的临床分析, 昆明医科大学学报. 2024, 45(1): 67-72. doi: 10.12259/j.issn.2095-610X.S20240111
    [9] 字正梅, 郝丽, 王南, 陈国兵.  联合预测模型对急性A型主动脉夹层破裂的诊断价值, 昆明医科大学学报. 2023, 44(1): 35-40. doi: 10.12259/j.issn.2095-610X.S20230128
    [10] 杨镜玉, 许宁, 张雨涛, 黄凤昌, 蒋元明, 殷亮.  基于基线T2WI图对进展期直肠癌经放化同期治疗后病理完全缓解的多模型预测比较, 昆明医科大学学报. 2023, 44(5): 117-124. doi: 10.12259/j.issn.2095-610X.S20230512
    [11] 周敏, 马智慧, 李加艳, 范建华, 林灵, 余亭颖, 张慧芳, 刘立.  肝硬化并胸水再发危险因素的预测, 昆明医科大学学报. 2022, 43(5): 149-154. doi: 10.12259/j.issn.2095-610X.S20220524
    [12] 贾吉, 陶四明.  基于血浆渗透压建立急性ST段抬高型心肌梗死重症患者发生院内死亡风险预测模型, 昆明医科大学学报. 2022, 43(12): 58-65. doi: 10.12259/j.issn.2095-610X.S20221212
    [13] 王虹, 杨德兴, 王强, 周维钰, 唐杰夫, 王振方, 付凯, 刘圣哲, 刘荣.  ICU脓毒症患者发生再喂养综合征的危险因素分析及预测模型建立, 昆明医科大学学报. 2022, 43(11): 44-51. doi: 10.12259/j.issn.2095-610X.S20221102
    [14] 张桂梅, 陈蜀, 宋云华, 吴阳, 周虹媛.  AECOPD患者再入院危险因素分析及预测模型的构建, 昆明医科大学学报. 2022, 43(8): 184-190. doi: 10.12259/j.issn.2095-610X.S20220830
    [15] 樊浩, 刘幸, 沈凌筠, 李海雯, 余春红, 李婧炜.  基于AAT与细胞因子构建活动性肺结核的Nomogram预测模型, 昆明医科大学学报. 2022, 43(8): 106-112. doi: 10.12259/j.issn.2095-610X.S20220816
    [16] 戴靖榕, 李婕, 何旭, 李杨, 李燕.  云南某医院老年医学科住院患者抑郁状态的危险因素分析及风险预测模型的构建, 昆明医科大学学报. 2021, 42(11): 20-26. doi: 10.12259/j.issn.2095-610X.S20211104
    [17] 刘幸, 刘媛, 刘向芳, 陈洁, 陈永刚, 孙辉, 马玲.  头孢哌酮钠舒巴坦钠不良反应分析与预测模型构建, 昆明医科大学学报. 2021, 42(6): 139-145. doi: 10.12259/j.issn.2095-610X.S20210622
    [18] 杨镜玉, 李文亮.  MSI在结直肠癌中的应用研究新进展, 昆明医科大学学报. 2020, 41(04): 150-156.
    [19] 胡万芹.  血清Th1/Th2细胞因子联合β-HCG预测先兆流产结局的价值, 昆明医科大学学报. 2012, 33(06): -.
    [20] 万蓉.  ARIMA乘积季节模型在食源性疾病月发病率预测中的应用, 昆明医科大学学报. 2012, 33(06): -.
  • 期刊类型引用(2)

    1. 赵琳,樊少军,于建云,王文杰. 大客车正面碰撞车内乘员损伤AIS评分的个案分析. 昆明医科大学学报. 2018(04): 125-129 . 本站查看
    2. 屠建权. 道路交通事故致颅脑损伤的法医学伤残评定分析. 世界最新医学信息文摘. 2017(19): 176-177 . 百度学术

    其他类型引用(0)

  • 加载中
图(7) / 表(4)
计量
  • 文章访问数:  16
  • HTML全文浏览量:  18
  • PDF下载量:  1
  • 被引次数: 2
出版历程
  • 收稿日期:  2024-10-19

目录

/

返回文章
返回