集成机器学习构建透明细胞肾细胞癌免疫影像分型及预测患者生存的价值.pdf

集成机器学习构建透明细胞肾细胞癌免疫影像分型及预测患者生存的价值李巧1，王宇昊2，夏一凡1，张玉东1*1南京医科大学第一附属医院放射科，2泌尿外科，江苏南京210029摘要 目的：基于集成机器学习（machine learning，ML）及增强CT构建透明细胞肾细胞癌（clear cell renal cell carcinoma，ccRCC）免疫影像（immunoradiomics，ImRad）分型并探讨其对患者预后生存的预测价值。方法：收集癌症基因图谱（the cancergenome atlas，TCGA）数据库113例ccRCC患者的临床、影像及基因表达信息，提取全肿瘤影像组学特征，基于集成ML进行特征筛选并构建免疫浸润、肿瘤突变负荷、免疫耗竭相关基因的ImRad分型。进一步多因素Cox回归分析ImRad对患者总生存（overall survival，OS）的预测效能。结果：经ML构建30个ImRad分型，经五折法验证，基于朴素贝叶斯算法对肿瘤免疫微环境的预测效能最佳（曲线下面积0.7170.956）。与基于临床、病理及ImRad单模态指标对比，融合临床病理及ImRad的模型预测OS的效能最佳。ImRad特征中，Rad激活态肥大细胞等8个特征是OS的独立预测因子。结论：基于集成ML及CT组学分析可预测ccRCC免疫微环境并提高预测患者术后生存的效能。关键词 透明细胞肾细胞癌；CT影像组学；机器学习；免疫浸润中图分类号 R737.11文献标志码 A文章编号 10074368（2023）09126508doi：10.7655/NYDXBNS20230913Immuneradiomic phenotype based on ensemble machine learning in predicting survivalof clear cell renal cell carcinomaLI Qiao1，WANG Yuhao2，XIA Yifan1，ZHANG Yudong1*1Department of Radiology，2Department of Urology，the First Affiliated Hospital of Nanjing Medical University，Nanjing 210029，ChinaAbstract Objective：To construct hybrid immuneradiomic（ImRad）phenotypes of clear cell renal cell carcinoma（ccRCC）basedon ensemble machine learning（ML）and contrastenhanced CT，and to investigate its predictive value for survival.Methods：Theclinical，CT imaging and geneexpression information of 113 ccRCC patients were collected from TCGA database.Radiomic featureswere extracted from whole tumor.ImRad predictors were constructed on tumor immune infiltration，tumor mutational burden，immuneexhaustion gene expression after feature selection based on ensemble ML.Further，the predictive value of ImRad for overall survival（OS）was assessed using multivariate Cox regression analysis.Results：Among the 30 ImRad contructed by ensemble ML andvalidated by 5 folds cross validation，Naive Bayes algorithm achieved the generally best performance（area under the curve：0717-0.956）.Clinicopathologic immune comprehensive model on predicting OS achieved the best performance，surpassing single modalityindicators based on clinical，pathological and ImRad.Among ImRad features，RadMast_cells_activated was one of eight independentpredictive factors for patientsprognosis.Conclusion：CT radiomics based on ensemble ML can predict immune microenvironmentand improve the prediction efficiency of postoperative survival of ccRCC.Key words clear cell renal cell carcinoma；CT radiomics；machine learning；immune infiltrationJ Nanjing Med Univ，2023，43（09）：12651272基金项目国家自然科学基金（82272082）通信作者（Correspondingauthor），Email：影像医学研究南京医科大学学报（自然科学版）Journal of Nanjing Medical University（Natural Sciences）第43卷第9期2023年9月1265南京医科大学学报第43卷第9期2023年9月透明细胞肾细胞癌（clear cell renal cell carcinoma，ccRCC）是肾细胞癌最常见的类型 1-2，相比乳头状和嫌色细胞亚型，更容易进展、转移，预后较差 3-4。患者状态、肿瘤分期、细胞核分级等是临床实践中常用的预后指标。然而，随着ccRCC个体化治疗的不断发展，特别是晚期肾癌，免疫治疗仅对20%30%的患者有效，目前尚无有效靶标对患者进行筛选和疗效预测5。肿瘤微环境包括细胞外基质、可溶性分子和肿瘤基质细胞，对肿瘤患者的免疫治疗疗效及预后影响重大 6-7。肿瘤突变负荷（tumor mutational burden，TMB）即每个肿瘤基因组编码区域的突变总数，是免疫治疗的生物标志物，TMB越高的肿瘤患者免疫治疗效果往往越好 8-9。微卫星不稳定性（microsatelliteinstability，MSI）即重复 DNA 束中核苷酸的自发丢失或获得，可以辅助诊断胃肠道、子宫内膜和结直肠肿瘤的表型，为各种癌症提供治疗决策信息 10-11。基于机器学习（machine learning，ML）的CT影像组学分析可捕获肉眼无法识别的图像特征，并揭示病理生理相关的潜在生物医学图像变化，预测治疗反应和生存期 11-13。因此，本研究利用集成ML算法，将ccRCC影像组学特征与肿瘤免疫微环境表达相关联，建立多模通道的免疫影像（immunoradiomics，ImRad）分型，并进一步研究该表型对患者术后生存的预测价值。1对象和方法1.1对象纳入南京医科大学第一附属医院经手术病理证实为ccRCC、术前行肾脏CT增强扫描且具有完整临床病理及免疫基因资料的患者113例。其中，男77例，女36例，年龄（59.412.9）岁（2688岁）；AJCC分期68.1%为T12期，31.9%为T34期；41.6%为低级别细胞核分级，58.4%为高级别。113例患者均有标准随访结果，随访内容包括治疗后有无复发及总生存时间（overall survival，OS），随访时间为186个月，删失病例OS按最长随访时间计算。中位OS为31（1054）个月。患者临床、影像数据从公开数据库癌症影像档案（the cancer imaging archive，TCIA）下载（http：/），对应免疫、基因数据从癌症基因图谱（the cancer genome atlas，TCGA）库中获得（https：/portal.gdc.cancer.gov）14。利用癌症样本转录谱的特性来推断肿瘤细胞以及浸润的各种正常细胞含量，获得3个CD8+T浸润评分、TMB、MSI、4种免疫耗竭相关基因（LAG3、CD244、PDCD1及TIGIT）15以及21种免疫细胞得分（使用R语言estimate包中的ESTIMATE算法），根据得分中位数分为高、低表达两组。1.2方法1.2.1CT图像勾画和特征提取使用国内开发的软件 SRhythm Multi Label 和Oncology Imaging Analysis（OCIA，上海师范大学上海磁共振重点实验室），分别由1名具有3年腹部影像诊断经验的放射科医师和1名具有6年临床手术经验的泌尿外科医生在CT皮质期轴位图像上对病灶边缘进行逐层勾画，然后融合成三维感兴趣区（region of interest，ROI）；每一次勾画均由另1名具有 15 年腹部影像诊断经验的放射科医师进行校对。校准完毕后利用Pyradiomics软件包提取肿瘤影像组学特征，包括14个形状特征、18个一阶直方图特征以及75个纹理特征、24个灰度共生矩阵（graylevel cooccurrence matrix，GLCM）特征、16个灰度游程矩阵（gray level runlength matrix，GLRLM）特征、16 个灰度区域大小矩阵（gray level size zone matrix，GLSZM）特征、5个邻域灰度差分矩阵（neighborhoodgraytone difference matrix，NGTDM）特征、14个灰度相关矩阵（gray leveldependencematrix，GLDM）特征 16。1.2.2模型构建构建肿瘤ImRad：降维及特征筛选。考虑高维数据ML建模产生的过拟合风险，首先采用4种特征筛选算法，即信息增益（information gain，Info.gain）、增益率（gain ratio）、基尼系数下降（gini decrease，Gini）、快速相关滤波器（fast correlation basedfilter，FCBF），分别对影像组学特征降维筛选。集成ML建模。利用随机森林（random forest，RF）、梯度增强（gradient boosting，GDBT）、逻辑回归（logisticregression，LR）、支 持 向 量 机（support vector machines，SVM）、朴素贝叶斯（naive bayes，NB）、k最近邻（knearest neighbor，kNN）、神经网络（neural network，NN）及堆叠式集成学习（stacking learning）8种算法在经筛选后的20个组学特征上建立二分类预测模型 17，预测目标包含3个CD8+T浸润评分、1个TMB表达、1个MSI表达、4个免疫耗竭相关基因及21种免疫细胞浸润表达结果，共建立30个ImRad表型。多种分类器算法的消融实验以获得最优算法，五折法用于小样本数据下模型的准确性验证。生存分析：纳入指标包括患者性别、年龄、1266AJCC TNM分期、细胞核分化程度、ImRad分型、3个CD8+T浸润评分、TMB、MSI、4种免疫耗竭关键基因及21种免疫细胞表达。对于维度较多的临床病理免疫影像（clinicopathologicimmunoradiomics，ClinicPtImRad）融合模型，利用LassoCox回归将不相关特征的回归系数缩小到零，得到系数不为零的特征。使用最大相关最小冗余法（mRMRe）基于有无肿瘤复发对851个组学特征进行降维，获取30个关键特征构建影像组学（radiomics，Rad）预后模型。进一步多因素Cox回归分析，最终构建如下预后模型：Rad 及 ImRad 模型；临床病理免疫（clinicopathologic immune，ClinicPtIm）模型；ClinicPtImRad模型。1.3统计学方法采用受试者工作特征（receiver operating characteristic，ROC）曲线和 ROC 曲线下面积（area undercurve，AUC）评价分类器效能。多因素Cox回归分析计算各模型中预测因子高、低表达组间OS的风险比（hazards ratio，HR）及其 95%可信区间（confidenceinterval，CI）。Cindex 评价模型预测效能。采用KaplanMeier生存曲线分析独立预测因子的生存影响并进行LogRank 检验。采用SPSS 26、R3.6.1和MeCalc19.9.4软件进行统计学分析。双尾P 0.05为差异有统计学意义。2结果2.1ImRad免疫影像表型构建在两步骤算法的 32 种组合中（4 种特征筛选8 种ML分类器算法），NB算法较其他分类器建模有更好的分类效能（AUC：0.7170.956，表1）。基于NB算法构建的ImRad表型能够有效监测肿瘤组织记忆B细胞（AUC=0.954）、静息态记忆CD4+T细胞（AUC=0.942）、激 活 态 自 然 杀 伤 细 胞（AUC=0.931）、嗜酸性粒细胞（AUC=0.956）和调节性T细胞（AUC=0.908）浸润（P 0.05）。2.2多模组学融合预测OS效能分别对比分析Rad、ImRad、ClinicPtIm、ClinicPtImRad对OS的预测效能。采用梯度提升树算法（gradient boosting machine，GBMCox）对 Rad 预测因子建模，获取 1 年、3 年、5年OS的校准曲线图（图1）。30个经mRMRe降维后的关键Rad特征中，waveletLLL_glcm_Imc1、wavelet LHL_ngtdm_Strength、wavelet_HHH_glrlm_Long Run Low Gray Level Emphasis为独立预测因子（表2），Rad模型预测OS的Cindex为0.756（95%CI：0.6360.876）。30个基于独立预测因子建立ImRad预后指数（表2，图2A），其预测OS的Cindex为0.857（95%CI：0.7870.927），优于Rad模型。ClinicPtIm指标中，AJCC分期、TMB、初始B细胞等11个指标是患者OS的独立预测因子，ClinicPtIm模型预测OS的Cindex为 0.924（95%CI：0.8720.976），优于 Rad 和 ImRad（表2，图2B）。融合模型ClinicPtImRad预测OS的Cindex为0.938（95%CI：0.9020.974），优 于 Rad、ImRad 和ClinicPtIm模型（表3，图2C）。各独立预测因子的生存分析中，Rad静息态树突状细胞、Rad激活态肥大细胞及RadM1巨噬细胞3个免疫影像特征的KaplanMeier曲线在高、低表达组间存在显著差异（P 0.05）。AJCC分期越早预后越好；激活态肥大细胞及MSI高表达均提示预后不佳，而静息态树突状细胞高浸润提示预后良好，其高、低表达组间生存曲线均存在显著差异（P 0.05，图3）。3讨论ccRCC多位于肾皮质，单发多见，通常大小不一，圆形或椭圆形，和周围肾实质分界较清晰，可有假包膜，常有钙化、出血、坏死及囊变。显微镜下肿瘤细胞体积较大，圆形或多边形，胞质丰富，透明或颗粒状，间质具有丰富的毛细血管和血窦。CT平扫表现主要取决于肿瘤本身的成分和血供的差异，以及肿瘤内有无出血、坏死、囊变、钙化等。通常CT平扫时密度与周围肾实质呈等或略低密度，体积较小的肿瘤密度多均匀，而体积较大者密度多不均匀，与瘤体内容易出血有关。绝大多数ccRCC为富血供肿瘤，增强扫描动脉期多明显不均匀或条纹状强化，强化程度多与邻近肾皮质相仿，甚至略高，静脉期及延迟期强化程度较周围肾实质减低，多呈“快进快出”的特点。病理研究显示其早期高强化与透明细胞的小泡状结构密切相关。本研究基于CT增强图像及集成ML构建ccRCC的免疫影像表型，并探讨其对患者预后的预测价值。目前基于全肿瘤影像组学特征分析主要用于肿瘤细胞异质性研究12，而本研究通过集成ML及全肿瘤ROI影像组学特征分析，构建ccRCC的多重免疫影像分型，在技术方法及临床应用方面具有一定探索性。本研究通过影像免疫关联，构建了30个第43卷第9期2023年9月李巧，王宇昊，夏一凡，等.集成机器学习构建透明细胞肾细胞癌免疫影像分型及预测患者生存的价值 J.南京医科大学学报（自然科学版），2023，43（9）：1265-12721267南京医科大学学报第43卷第9期2023年9月ImRad免疫影像分型，其中8个ImRad分型与患者OS独立相关，基于ImRad分型构建的Cox模型对OS的预测效能优于传统方法。而且发现Rad静息态树突状细胞、RadM1巨噬细胞及Rad激活态肥大细胞3个ImRad标志物表达水平不同，其OS存在显著差异，提示免疫影像分型对患者临床预后具有较好的预测价值，从而为优化临床治疗方案提供重要参考依据。与既往研究对比，本研究有以下创新：首先，通过泌尿专业及影像专业医生双盲参与肿瘤ROI勾画，对于形态不规则的肿瘤，结合临床和影像的经验，如强化环、肿瘤包膜隆突等，从而提高勾画精准性，并由经验更丰富的高年资医师进行校准。其次，本研究在技术上进行了多重算法对比及集成，并通过消融实验获得最优参数配置，提升了模型的鲁棒性，五折法用于小样本数据模型的准确性验证，所构建的30个ImRad分型中，基于贝叶斯算法的模型表现最佳，在多重验证中均取得比较可靠的结果。另外，在Cox回归分析中，通过Rad、ImRad、ClinicPtIm以及ClinicPtImRad多重对比，证实基于影像免疫关联的ImRad模型优于单纯Rad模型，ClinicPtImRad优于ClinicPtIm模型。该结果提示，多重免疫与高通量影像特征具有微观宏观上的关联性，通过构建免疫影像分型，有望为临床提供类别肿瘤基质评分肿瘤免疫细胞浸润评分肿瘤纯度评分TMBMSI初始B细胞记忆B细胞浆细胞CD8+T细胞静息态记忆CD4+T细胞激活态记忆CD4+T细胞辅助性滤泡T细胞调节性T细胞 T 细胞静息态自然杀伤细胞激活态自然杀伤细胞单核细胞M0巨噬细胞M1巨噬细胞M2巨噬细胞静息态树突状细胞激活态树突状细胞静息态肥大细胞激活态肥大细胞嗜酸性粒细胞中性粒细胞LAG3CD244TIGITPDCD1NB0.759*0.806*0.780*0.772*0.724*0.813*0.954*0.847*0.744*0.942*0.729*0.774*0.908*0.765*0.820*0.931*0.739*0.717*0.767*0.742*0.724*0.869*0.869*0.889*0.956#0.767*0.806*0.752*0.835*0.787*Stack0.6420.7600.7920.6920.6850.7640.7870.6700.6920.7910.6660.6720.7670.5990.6640.6910.6320.6060.6070.6980.6210.7310.6620.8050.8340.6700.7420.7390.7620.718GB0.6130.8130.7770.5690.6660.6480.570.6380.5130.6440.6050.6590.7020.6680.7280.5580.7210.5950.6740.7240.6330.5460.5900.7860.6260.6540.7550.6650.7620.817*RF0.6510.821*0.7950.6110.6460.6410.7880.5840.6090.8010.6170.6980.6750.6410.7370.5130.6110.5440.6690.6890.6420.5490.6610.901*0.7570.6170.7140.6450.7770.757NN0.6150.7850.832*0.4110.6150.5300.7250.5340.4100.7070.5500.6480.7230.5700.5760.5470.5810.5260.6140.5950.5720.3770.6880.8030.8010.6240.7290.6430.7580.702LR0.5030.5340.6370.5450.5560.4280.5670.5170.5230.6500.4100.5380.6950.5670.7180.3060.3800.5990.6420.6450.5970.4990.7070.6890.7150.6550.6050.4990.5930.611kNN0.5670.6120.5800.5970.5580.5450.6830.6070.5710.5730.3920.6240.5990.6220.6160.5470.6250.6550.6270.5900.6090.6340.4690.7590.5520.5730.5610.5610.6730.628SVM0.4360.7900.7730.4340.5450.5170.7650.4160.4300.5060.3860.6600.5750.6640.6630.4040.4070.4390.6780.6680.5830.5040.5320.7940.4320.5620.7390.6950.7530.731表1集成机器学习构建免疫影像ImRad表型及其AUCTable 1ImRad and AUC constructed by ensemble machine learning algorithms降维算法包括FCBF、Info.Gain、gain ratio、Gini，分类算法包括RF、GBDT、LR、SVM、NB、kNN、NN。*：为各种分类器中最高AUC；#：为显著高AUC值。1268一种可靠的无创性影像检查方法来预测肿瘤免疫微环境状态，从而为提高患者免疫治疗效果提供可靠依据。以往研究发现TMB可用于预测肺癌、膀胱癌、黑色素瘤及结直肠癌等实体肿瘤的免疫治疗疗效18，MSI决定了胃肠癌患者对免疫治疗的反应10-11，笔者首次将TMB和MSI应用于ccRCC中，其对患者术后预后的预测价值在 ClinicPtIm 和ClinicPtImRad模型中均得到验证，高MSI和TMB预示着复发死亡风险更高，集成ML构建RadTMB和RadMSI的AUC、精确度及准确度均0.70，揭示肿瘤TMB及MSI表型与影像组学特征具有一定关联性。同时本研究发现，高AJCC分期（和）是不良预后的独立风险因素，这与既往研究一致3-4。另外对CD8+T细胞浸润的3个关键评分分析，发现高Rad肿瘤基质评分及肿瘤纯度评分提示预后良好，可能因为免疫炎症型（具有密集的CD8+T细胞浸润）肿瘤免疫微环境具有较好的免疫原性和对免疫治疗的应答性，从而提高患者远期生存率19。21种免疫细胞浸润表达及其ImRad分型中，Rad静息态树突状细胞及静息态树突状细胞浸润均提示预后良好，其OS曲线在高、低表达组间均存在显著差异；RadM1巨噬细胞及M1巨噬细胞浸润提示预后良好，RadM1巨噬细胞的生存曲线在高、低表达组间存在显著差异，说明影像免疫分型有预测免疫细胞活性的能力；Rad浆细胞高表达提示预后良好，但浆细胞浸润预后不良，这可能与小样本数据分布不均有关；高Rad激活态肥大细胞及激活态肥大细胞患者预后不良，其生存曲线高、低表达组间均存在显著差异；高Rad嗜酸性粒细胞表达的患者预后不预测生存00.20.40.60.81.01年OS3年OS5年OS实际生存1.00.80.60.40.20基于终点状态有无复发使用最大相关最小冗余法对851个影像组学特征进行降维得到30个特征，采用最新梯度提升树算法（GBM）进行Rad生存分析，并绘制1年、3年、5年生存期的校准曲线图（R语言的mRMRe包、survival包、gbm包及rms包）。图1Rad模型预测1年、3年、5年生存期的校准曲线图Figure 1Calibration curve of Rad on predicting 1year，3year，5year survival period其上3个模型系数Omnibus检验的以2为底的对数似然值分别为176.065、115.517、191.548。预后模型RadImRadClinicPtIm独立预测因子waveletLLL_glcm_Imc1waveletHHH_glrlm_Long Run Low Gray Level EmphasiswaveletLHL_ngtdm_StrengthRadPDCD1Rad肿瘤基质评分Rad浆细胞RadM1巨噬细胞Rad静息态树突状细胞Rad激活态树突状细胞Rad激活态肥大细胞Rad嗜酸性粒细胞AJCC分期肿瘤免疫细胞浸润评分肿瘤纯度评分TMB初始B细胞静息态自然杀伤细胞M0巨噬细胞M1巨噬细胞M2巨噬细胞静息态树突状细胞激活态肥大细胞Cindex（95%CI）0.756（0.6360.876）0.857（0.7870.927）0.924（0.8720.976）HR（95%CI）00.01（0.000.05）00.50（0.250.99）28.23（24.7558.41）00.20（0.070.54）00.14（0.040.42）00.25（0.070.89）00.20（0.070.55）00.10（0.030.31）00.21（0.070.60）43.70（7.82244.9）03.75（1.0213.60）63.70（13.60298.70）07.79（1.5439.30）00.05（0.010.38）10.40（2.2049.30）10.10（1.7458.90）00.05（0.010.36）00.19（0.050.77）00.18（0.050.61）06.57（1.7824.20）00.11（0.030.44）151.00（26.00876.90）表2Rad、ImRad及ClincPtIm模型多因素Cox分析Table 2Multivariate Cox analysis of Rad，ImRad and ClincPtIm modelsP值0.0060.0490.0090.002 0.0100.0330.002 0.0100.004 0.0100.045 0.0100.0130.0030.0030.0100.0030.0200.0060.0050.002 0.010第43卷第9期2023年9月李巧，王宇昊，夏一凡，等.集成机器学习构建透明细胞肾细胞癌免疫影像分型及预测患者生存的价值 J.南京医科大学学报（自然科学版），2023，43（9）：1265-12721269南京医科大学学报第43卷第9期2023年9月0.010.1110100RadPDCD1Rad肿瘤基质评分Rad浆细胞Rad巨噬细胞Rad静息态树突状细胞Rad激活态树突状细胞Rad激活态肥大细胞Rad嗜酸性粒细胞Common effect modelHeterogeneity：I2=87%，2=3.539 6，P0.010.200.140.260.200.110.2143.793.750.0020.0100.0330.0020.0100.0040.0100.0450.080.540.050.430.070.890.070.560.040.320.070.617.83244.961.0313.690.220.48Hazard Ratio95%CIHRP值ImRadCindex=0.8570.01 0.1110100AJCC分期肿瘤免疫细胞浸润评分肿瘤纯度评分TMB初始B细胞激活态记忆CD4+T细胞静息态自然杀伤细胞M0巨噬细胞M1巨噬细胞M2巨噬细胞静息态树突状细胞激活态肥大细胞Common effect modelHeterogeneity：I2=91%，2=6.891 4，P0.0163.797.800.0610.4310.133.510.050.200.186.580.12151.00Cindex=0.9240.0100.0130.0030.0030.0100.0720.0030.0200.0060.0050.0020.01013.62298.731.5539.320.010.392.2149.321.7458.900.8913.770.0100.360.050.770.060.621.7824.280.030.4426.00876.911.102.62Hazard Ratio95%CIHRP值Clinic_Pt_ImCindex=0.938性别AJCC分期MSI初始B细胞浆细胞静息态树突状细胞激活态肥大细胞RadPDCD1Rad肿瘤基质评分Common effect modelHeterogeneity：I2=92%，2=5.404 3，P0.010.3263.564.0510.1412.640.3639.850.130.160.0400.0100.0350.0030.0050.0490.0100.0100.0010.110.9513.52298.831.1014.892.2445.952.1474.680.131.008.22193.280.040.400.060.470.731.68Hazard Ratio95%CIHRP值Clinic_Pt_lmRad0.010.1110100CAB图2ImRad（A）、ClinicPtIm（B）和ClinicPtImRad（C）的森林图Figure 2Forest maps of ImRad（A），ClinicPtIm（B），and ClinicPtImRad（C）models该模型系数的Omnibus检验，以2为底的对数似然值为111.313。ClinicPtImRad性别AJCC分期MSI初始B细胞肥大细胞静息态树突状细胞激活态肥大细胞RadPDCD1Rad肿瘤基质评分LASSO系数-0.20-2.22-0.27-0.11-0.33-0.06-1.57-0.88-0.71Cindex（95%CI）0.938（0.9020.974）HR（95%CI）00.32（0.110.95）63.50（13.50298.80）04.05（1.1014.80）10.10（2.2245.90）12.60（2.1474.60）00.36（0.130.99）39.80（8.22193.20）00.13（0.040.39）00.16（0.050.46）P值0.040 0.0100.0350.0030.0050.049 0.010 0.0100.001表3ClinicPtImRad的LASSOCox回归模型Table 3LASSOCox regression of ClinicPtImRad model1270良，高Rad激活态树突状细胞的患者预后良好而激活态树突状细胞浸润对预后无影响，既往研究表明嗜酸性粒细胞同时有促肿瘤和抗肿瘤功能20，但具体截断值有待进一步研究。树突状细胞介导的抗原特异性T细胞抗癌免疫治疗在临床上取得了巨大成功21。另外，静息态自然杀伤细胞、M0巨噬细胞、M2 巨噬细胞高表达是预后良好的独立预测因子。4 种 免疫耗竭关键基因中（LAG3、CD244、PDCD1、TIGIT），RadPDCD1 是 ccRCC 患者术后复发死亡的独立预测因子，高RadPDCD1表达预后良好，这与 Miao等22研究结果稍有不同，两组结果差异的原因有待进一步扩大样本量研究。本研究存在局限性，第一，来自公共数据集的图像质量存在较大差异，扫描参数不尽相同，这可能会影响分析结果的准确性。第二，本研究是回顾性研究，存在样本选择偏倚的风险。第三，研究样本量有限，将113例ccRCC进行五折交叉验证以保证预测模型的过拟合风险，准确性及可重复性需要在独立验证集中进一步评估。此外，阳性表达记忆B细胞、浆细胞、静息态记忆CD4+T细胞、调节性T细胞、激活态自然杀伤细胞和嗜酸性粒细胞等免疫细胞浸润的患者数量较少，可能会导致测量偏倚，这也是其他组学研究中常面临的技术问题 23。因此，在临床应用之前，基于有限数据集的研究结果需要通过多中心和大规模的研究进行前瞻性验证。总之，本研究通过影像免疫表达关联，证实增强CT影像组学特征在预测ccRCC患者的肿瘤免疫微环境和TMB方面具有良好的临床应用潜力，构建的ImRad免疫影像分型可用于ccRCC患者的术后生存预测，将有助于临床医生对ccRCC患者进行更好的风险分层和临床治疗。参考文献1 SUNG H，FERLAY J，SIEGEL R L，et al.Global cancerstatistics 2020：GLOBOCAN estimates of incidence andmortality worldwide for 36 cancers in 185 countriesJ.CA Cancer J Clin，2021，71（3）：209-2492 BUKAVINA L，BENSALAH K，BRAY F，et al.Epidemiology of renal cell carcinoma：2022 update J.Eur Urol，HR=2.4295%CI：1.065.56P0.05020406080100时间（月）Rad静息态树突状细胞高表达（n=69）低表达（n=44）100806040200生存率（%）HR=2.9995%CI：1.257.15P=0.01020406080100时间（月）静息态树突状细胞100908070605040生存率（%）高表达（n=73）低表达（n=40）020406080100时间（月）M1巨噬细胞100908070605040生存率（%）高表达（n=58）低表达（n=55）HR=39.8295%CI：9.51166.71P0.01020406080100时间（月）激活态肥大细胞100908070605040302010生存率（%）高表达（n=11）低表达（n=102）HR=59.8995%CI：9.72369.76P0.01020406080100时间（月）记忆B细胞100908070605040302010生存率（%）高表达（n=7）低表达（n=106）HR=7.8995%CI：3.0020.53P0.01020406080100时间（月）AJCC分期1009080706050403020生存率（%）高表达（n=36）低表达（n=77）HR=2.9395%CI：1.296.67P=0.01020406080100时间（月）MSI100908070605040生存率（%）高表达（n=58）低表达（n=55）HR=2.7295%CI：1.216.13P0.05020406080100时间（月）RadM1巨噬细胞1009080706050403020生存率（%）高表达（n=56）低表达（n=57）HR=58.7795%CI：9.55361.52P0.01020406080100时间（月）Rad激活态肥大细胞100806040200生存率（%）高表达（n=7）低表达（n=106）图3部分独立预测因子的KaplanMeier生存曲线Figure 3KaplanMeier survival curves of partial independent predict