骨髓瘤和免疫区室的纵向分析揭示失调模式及与多发性骨髓瘤应答的关联

Longitudinal Profiling of Tumor and Immune Compartments Uncovers Patterns of Dysregulation and Associations with Response in Multiple Myeloma

作者信息Denis J Ohlstrom, William C Pilcher Jr, Marina E Michaud, Chaitanya Acharya, Sarthak Satpathy, Edgar Gonzalez-Kozlova, Reyka G Jayasinghe, Katherine Ferguson, Hope L Mumme, Shivani Nanda, Yizhe Song, Sowmitri Mantrala, Dimitra Karagkouni, Jessica Schulman, Nick Pabustan, Junia Vieira Dos Santos, Daniel W Sherbenou, Jonathan J Keats, Alexander M Gout, Steven Foltz, Alessandro Lagana, Taxiarchis Kourelis, Ravi Vij, Madhav V Dhodapkar, David Avigan, Hearn Jay Cho, Linda B Baughn, Ajay K Nooka, Sagar Lonial, Shaji Kumar, Mehmet K Samur, Ioannis S Vlachos, Li Ding, Sacha Gnjatic, George Mulligan, Manoj K Bhasin
PMID41364805
发布时间2026-03-04
DOI10.1158/2643-3230.BCD-25-0205

实验完整度

包含多中心纵向队列、CD138阴性单细胞测序、CD138阳性批量RNA测序和全基因组测序,并进行差异丰度、表达、通路和生存分析等多层级验证。

主要模型

多发性骨髓瘤患者骨髓样本 CD138阴性骨髓单个核细胞 CD138阳性骨髓瘤细胞

重点核对

样本来自MMRF CoMMpass研究(NCT01454297)的102例患者,共243份CD138阴性骨髓样本 CD138阴性单细胞RNA测序,10x Genomics Chromium平台,目标捕获5000个细胞/样本 CD138阳性批量RNA测序和全基因组测序,用于分析肿瘤区室 PFS分组:以中位PFS(1491天)为界,分为高于和低于中位PFS组 Shannon多样性指数用于评估免疫球蛋白多样性,并基于血清Ig测量值在推导和验证队列中验证

摘要

多发性骨髓瘤是一种由克隆性扩增的浆细胞构成的恶性肿瘤,其形成受到与免疫微环境(IME)复杂相互作用的影响。为了研究治疗应答和疾病进展的免疫相关因素,我们对来自102例患者的243份骨髓样本(631,226个细胞)进行了包括CD138阴性单细胞RNA测序在内的多组学分析,并对209份样本进行了CD138阳性批量RNA和全基因组测序。在纵向分析中,干扰素-γ信号与自体干细胞移植后T细胞记忆受损的标志物相关,而初始B细胞丰度和免疫球蛋白多样性与改善的无进展生存期(HR = 0.48;P = 2.3e-4)相关。在疾病进展时,多发性骨髓瘤细胞上调了癌症-睾丸抗原(CTAg)和免疫效应基因,同时伴有B细胞耗竭、骨髓源性抑制细胞表达富集和表型T细胞耗竭。这些发现突出了动态的免疫-肿瘤相互作用,将初始B细胞重建确定为持久应答的生物标志物,并将CTAg确定为高危疾病进展的潜在靶点。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
在多发性骨髓瘤的治疗和疾病进展过程中,骨髓免疫微环境的动态变化及其与治疗应答和预后的关联是什么?
核心机制
ASCT后,IFN-γ和TNF信号通路驱动CD8+ T细胞向凋亡而非记忆表型偏移;初始B细胞重建与免疫球蛋白多样性改善和无进展生存期延长相关。在疾病进展时,骨髓瘤细胞上调癌症-睾丸抗原(CTAg),并伴随B细胞耗竭、单核细胞髓系来源抑制细胞样表型和CD8+ T细胞耗竭,这些可能由CTAg表达的细胞因子和趋化因子(如IL10、IL32)驱动。
主要证据
基于CD138阴性单细胞RNA测序和CD138阳性批量RNA测序/全基因组测序的纵向分析,显示ASCT后CD8+ T细胞中IFN-γ和TNF信号富集,初始B细胞比例与Ig多样性正相关(R=0.45, P=0.012),CTAg富集样本中B细胞减少、MDSC评分升高、CD8+ T细胞耗竭标志物增加。
研究意义
该研究揭示了多发性骨髓瘤中免疫-肿瘤相互作用的动态变化,确定了初始B细胞重建和Ig多样性作为持久应答的生物标志物,并指出CTAg表达与免疫抑制相关,可能作为高危疾病进展的潜在治疗靶点。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

构建纵向队列和单细胞图谱

表征多发性骨髓瘤患者骨髓免疫微环境在不同疾病阶段的细胞组成。

收集102例患者的243份CD138阴性骨髓样本,进行单细胞RNA测序,经质控和聚类后,注释出T细胞、NK细胞、B细胞、髓系细胞等主要免疫细胞类型。

2

比较诱导治疗与ASCT对免疫微环境的影响

评估不同一线治疗方案(单纯诱导 vs 诱导+ASCT)对免疫细胞组成和基因表达的影响。

比较基线与首次应答样本的免疫细胞亚群比例变化,并分析CD8+ T细胞和B细胞的转录程序。

3

分析CD8+ T细胞偏倚和B细胞重建的机制

揭示ASCT后CD8+ T细胞偏向凋亡而非记忆表型的分子机制,以及B细胞重建的促增殖特征。

对CD8+ T细胞进行GSEA、转录因子活性分析和细胞通讯分析,对B细胞进行GSEA和转录因子活性分析。

4

评估初始B细胞丰度与临床结局的相关性

确定ASCT后初始B细胞丰度和Ig多样性与无进展生存期的关系。

将31例ASCT后患者按中位PFS分为GMpfs和LMpfs组,比较初始B细胞比例,并计算基于Ig基因表达的Shannon多样性指数,同时分析血清Ig多样性与PFS的关联。

5

分析CD16+单核细胞与B细胞活化的潜在相互作用

探究CD16+单核细胞是否通过配体-受体相互作用支持B细胞活化和Ig多样性。

将非B细胞亚群丰度变化与Ig多样性指数相关,并进行细胞通讯分析,评估CD16+单核细胞和B细胞中相关基因表达。

6

鉴定进展期CTAg表达及其与免疫抑制的关联

识别疾病进展时骨髓瘤细胞的转录变化,并评估CTAg表达与免疫微环境改变的关系。

对67例患者的基线和首次进展CD138阳性样本进行差异表达分析,鉴定CTAg富集样本,并比较其与免疫细胞组成和基因表达的关系。

7

纵向追踪CTAg表达和克隆演化

评估CTAg表达在多次进展中的动态变化及其与染色质修饰基因突变的关系。

分析13例具有3个及以上纵向样本的患者,追踪CTAg表达和遗传改变,并进行亚克隆重建。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
Chromium控制器10x Genomics--
Chromium Next GEM单细胞3' GEM、文库与凝胶珠试剂盒v2.110x Genomics--
NovaSeq 6000测序仪Illumina--
Cell Ranger10x Genomicsv6.0.1
EasySepStemCell Technologies--
RoboSepStemCell Technologies--
死细胞去除试剂盒Miltenyi Biotec--
NIH/3T3细胞ATCCCRL-1658

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
样本采集和处理
样本来自MMRF CoMMpass研究,采集时间点(基线、应答、进展),CD138阴性细胞分离方法,冷冻保存条件(90%胎牛血清+10% DMSO,浓度5-30百万细胞/mL)
阅读提示:Methods: Sample Collection, CD138neg Cell Isolation and Cryopreservation
单细胞RNA测序
10x Genomics Chromium平台,3'文库试剂盒v2.1,目标捕获5000细胞/样本,测序深度25,000-50,000 reads/细胞
阅读提示:Methods: scRNA-seq Sample Preparation, Library Construction, and Sequencing
数据处理和质控
Cell Ranger版本,质控阈值(UMI≥1000,基因≥200,线粒体≤20%),批次校正方法(Harmony),双细胞去除工具
阅读提示:Methods: Summary of scRNA-seq Processing, Quality Control, and Clustering
差异丰度和表达分析
线性混合效应模型设置,差异表达显著阈值(adjusted P<0.05),GSEA方法,TF推断方法
阅读提示:Methods: Summary of Differential Abundance Analyses, Summary of Expression Analyses
Ig多样性分析
Shannon多样性指数计算基于Ig基因表达和血清Ig测量,患者纳入标准(接受单次ASCT,有PFS/OS数据,临床协变量完整)
阅读提示:Methods: Summary of Ig Diversity Analyses
CTAg表达分析
CTAg候选列表来源和筛选标准,CTAg富集定义,差异表达阈值
阅读提示:Methods: Transcriptomic Analyses of Longitudinal CD138pos Bulk RNA-seq Samples
遗传和亚克隆分析
WGS数据处理管道,SNV过滤标准,亚克隆重建参数(QuantumClone和clonevol设置)
阅读提示:Methods: Summary of Evaluation of Cytogenetic Alterations and Mutations in CD138pos Whole-Genome Sequencing, Summary of Subclonal Reconstruction