《Cancer Immunology, Immunotherapy》：RARRES1 marks an immune-cold, chemoresistance-associated malignant epithelial subpopulation enriched in pancreatic ductal adenocarcinoma
编辑推荐：
本研究针对胰腺导管腺癌（PDAC）治疗抵抗和免疫抑制微环境的临床难题，通过单细胞RNA测序（scRNA-seq）分析，鉴定出一个在疾病进展和肝转移样本中富集的恶性上皮亚群（c05）。该亚群具有免疫“冷”表型、多药抵抗特征，并高表达核心基因RARRES1。功能实验证实RARRES1促进吉西他滨抵抗，泛癌分析确立其作为不良预后的新型生物标志物。这项工作揭示了PDAC治疗抵抗的细胞机制，为靶向干预提供了新方向。
胰腺癌，尤其是其最主要的亚型——胰腺导管腺癌（PDAC），素有“癌王”之称，其凶险之处在于极高的死亡率与极差的治疗效果。尽管手术切除是潜在的治愈手段，但大多数患者确诊时已处晚期或发生转移，失去了手术机会。对于这些患者，化疗（如FOLFIRINOX或吉西他滨/白蛋白紫杉醇方案）成为一线选择。然而，疗效在患者间差异巨大，且常在治疗过程中出现耐药，最终导致疾病进展。另一个令人沮丧的现实是，PDAC的肿瘤微环境（TME）通常呈现免疫“冷”的特征，即缺乏能够攻击肿瘤的免疫细胞浸润，这使得在其他癌种中取得革命性成功的免疫检查点抑制剂疗法在PDAC中几乎无效。那么，驱动PDAC化疗抵抗和免疫逃逸的深层细胞机制究竟是什么？是否存在一个关键的、具有共同特征的恶性细胞亚群，同时介导了这两种致命的特性？这成为了研究者亟待攻克的科学谜题。
为了回答这些问题，研究团队将目光投向了能够在单细胞分辨率下解析肿瘤异质性的强大工具——单细胞RNA测序（scRNA-seq）。他们系统分析了来自27个PDAC样本（包括未经治疗和经化疗治疗的原发肿瘤及肝转移灶）的公开scRNA-seq数据，旨在跨越“治疗反应”和“转移状态”等多个维度，寻找那个与不良预后紧密相关的“罪魁祸首”细胞亚群。这项研究的结果最终发表于《Cancer Immunology, Immunotherapy》期刊。
研究人员开展研究主要应用了以下几个关键技术方法：首先，利用来自基因表达综合（GEO）数据库的27例PDAC样本（含17例原发、10例肝转移）的单细胞RNA测序数据进行细胞图谱构建和异质性分析。其次，通过推断拷贝数变异（CNV）从上皮细胞中鉴定恶性亚群。再者，运用基因集富集分析（GSEA）、细胞间通讯分析（CellChat）、非负矩阵分解（NMF）和拟时序轨迹分析（monocle）等多种生物信息学手段，深入刻画目标亚群的功能特征和演化关系。此外，还利用癌症基因组图谱（TCGA）和国际癌症基因组联盟（ICGC）的批量RNA测序数据进行生存分析和突变景观分析，并通过吉西他滨耐药细胞株的构建、蛋白质印迹法（Western blot）和细胞活力（MTT） assay等功能实验对核心基因进行验证。
研究结果
单细胞图谱揭示PDAC中与治疗反应和转移相关的细胞异质性
研究人员对包含不同治疗反应（部分缓解PR、疾病稳定SD、疾病进展PD）和部位（原发、肝转移）的样本进行单细胞分析，构建了包含8个主要细胞类型的图谱。他们发现，在疾病进展（PD）组中，上皮细胞的比例显著更高，而在部分缓解（PR）组中，T/NK/固有淋巴细胞（ILCs）、B/浆细胞等免疫细胞比例更高。与原发肿瘤相比，肝转移灶中上皮细胞比例也显著升高。这些结果表明，上皮细胞的富集与疾病进展和转移密切相关。
鉴定并表征一个化疗抵抗和转移性的恶性上皮亚群c05
通过对上皮细胞进行亚聚类和CNV分析，研究者从20个上皮亚群中鉴定出14个恶性亚群。其中，亚群c05表现出独特的特征：它在肝转移样本和化疗暴露样本中显著富集，尤其在疾病进展（PD）组中富集程度最高。生存分析显示，c05的特征基因高表达与TCGA及两个独立ICGC队列（PACA-CA, PACA-AU）中患者的总生存期缩短显著相关。通路富集分析表明，c05特异性激活了RAF/MAPK信号级联、MYC靶点等致癌通路以及铂类药物耐药等化疗抵抗通路。
亚群c05呈现免疫“冷”表型且细胞间通讯受损
细胞通讯分析显示，c05与其他恶性上皮细胞（尤其是经典型和基底样型）相比，发出的信号强度显著更弱。更重要的是，c05接收来自肥大细胞、髓系细胞、T/NK/ILCs和B/浆细胞等免疫细胞的信号输入非常有限，呈现出一种“绝缘”状态。基于批量RNA数据的免疫浸润分析（CIBERSORT）进一步支持了这一发现：c05特征高表达的样本表现出免疫抑制性调节性T细胞（Tregs）的增加，但缺乏抗肿瘤效应免疫细胞的浸润。NMF分析也鉴定出一个与c05富集相关的NMF5程序，该程序表现为免疫“冷”表型，且主要来源于PD样本；而免疫“热”的NMF1程序则主要来源于PR和SD样本。
亚群c05富集于进化轨迹终点并表现出多药抵抗表型
拟时序分析将c05和来源于PD样本的细胞定位在轨迹的终末分支，提示其处于一种分化的、治疗抵抗的状态。药物敏感性预测分析显示，与c05特征低表达样本相比，高表达样本对吉西他滨、5-氟尿嘧啶、奥沙利铂、伊立替康、紫杉醇等一线PDAC化疗药物的敏感性显著降低，即表现出多药抵抗。相反，这些样本对达沙替尼（酪氨酸激酶抑制剂）和曲美替尼（MEK抑制剂）表现出了更高的敏感性。突变景观分析还发现，c05高表达组具有更高的肿瘤突变负荷（TMB）和更高的KRAS突变频率（81% vs. 45%）。
RARRES1标记了PDAC中一个化疗抵抗和免疫排斥的表型
为了寻找驱动c05表型的核心基因，研究者交叉分析了四个对比组（c05 vs. 其他恶性亚群、PD vs. 其他反应组、转移 vs. 原发、化疗暴露 vs. 未治疗）的差异表达基因，发现RARRES1 (Retinoic Acid Receptor Responder 1) 是四个组中 consistently 上调的四个基因之一。功能实验证实，在吉西他滨耐药细胞株（L3.6pl-Res）中，RARRES1在mRNA和蛋白水平均显著上调。而利用短发夹RNA（shRNA）敲低RARRES1，能够有效恢复耐药细胞对吉西他滨的敏感性。拟时序分析显示RARRES1在轨迹终末的耐药阶段表达升高。肿瘤免疫功能障碍和排除（TIDE）分析表明，RARRES1高表达与更高的免疫排斥（exclusion）分数相关，提示其与免疫逃逸的空间排斥机制有关。此外，泛癌生存分析确立RARRES1高表达是肾上腺皮质癌、肾透明细胞癌、肺腺癌、卵巢癌等多种癌症不良预后的生物标志物。
研究结论与意义
本研究通过整合多组学分析，在PDAC中鉴定出一个独特的恶性上皮亚群c05。该亚群富集于化疗抵抗和肝转移样本中，与患者不良预后紧密相关，其特征是激活了促生存和耐药通路，处于肿瘤进化轨迹的终点，并塑造了一个免疫“冷”的微环境。研究者进一步将RARRES1确立为该亚群的核心上调基因，其表达不仅与化疗抵抗和免疫排斥表型相关，功能实验也证实其直接参与促进吉西他滨耐药。此外，RARRES1还被证明是一个潜在的泛癌不良预后生物标志物。
这项研究的意义在于，它超越了传统的批量转录组分型，在单细胞水平揭示了一个驱动PDAC治疗失败和转移的具体细胞实体（c05）及其关键分子标签（RARRES1）。这为理解PDAC的高度异质性和治疗抵抗提供了新的细胞机制视角。从临床转化角度看，c05特征或RARRES1表达有望成为预测化疗反应、识别高危患者的新型生物标志物。同时，该研究提示靶向c05/RARRES1轴（例如通过其分泌特性开发抗体偶联药物ADC，或利用其对达沙替尼等药物的敏感性设计联合疗法）可能成为克服PDAC耐药的新策略。尽管研究存在样本量、缺乏配对样本等限制，但它无疑为开发针对这一顽固癌症的精准治疗策略奠定了重要的理论基础。
