编者按：四联疗法已成为新诊断多发性骨髓瘤（NDMM）的标准诱导治疗，但治疗反应异质性显著，亟需分子预测因素指导早期治疗决策。近期，Nature Medicine发表MIDAS试验辅助研究，对423例NDMM患者进行整合基因组、转录组和表观遗传学分析，以刻画分子异质性及其与Isa-KRd四联诱导后可测量残留病（MRD）的关联。本文特整理核心要点，以飨读者。
为什么需要多组学MRD预测：四联时代风险分层需重新评估
MM是遗传和临床高度异质的浆细胞恶性肿瘤。基因组、转录组和甲基化组共同构成互补的分子维度，联合指导MM发病机制和患者分层。然而，在大型前瞻性临床试验队列中整合三个分子层面的分析仍缺乏。鉴于四联方案已成为NDMM新标准诱导治疗，临床结局的分子决定因素需在新标准治疗背景下重新评估。
MIDAS试验是一项随机、开放标签III期研究，根据统一Isa-KRd四联诱导后MRD状态评估巩固和维持策略。在本辅助研究中，研究者对423例NDMM进行整合多组学分析，以识别四联诱导后MRD状态的分子预测因素。
研究设计：423例MIDAS患者，WGS+RNA-seq+全基因组酶促甲基化测序
研究采用全基因组测序（WGS）、RNA测序和全基因组酶促甲基化测序。经严格质控，获得421例WGS、418例RNA-seq和422例甲基化数据，416例样本三种数据齐全。368/394例入组患者可评估诱导后MRD，44.8%在10⁻⁶灵敏度下为MRD阴性。
转录组亚组反应异质：HP1/HP3和MS组MRD阴性率最高，CD2和MF组最低
基于基因表达谱的无监督分类识别出八个分子组：CD1、CD2、MS、MF、LB和HP1-HP3。诱导后MRD阴性率在不同分子组间差异显著（P=1.9×10⁻¹⁰），校正临床协变量后仍显著。CD2和MF组MRD阴性率最低，HP1、HP3和MS组最高。
差异表达分析显示，MRD阳性样本中170个基因上调，包括CCND1、MAFA和B细胞分化相关基因PAX5；MRD阴性病例中217个基因上调，富集于表观遗传修饰因子、囊泡运输、细胞增殖调控和MYC靶点。B细胞分化程序在CD1/CD2和HP2组高表达，这些组MRD阴性率较低；MYC通路在HP1/HP3和MF组激活。分化状态和增殖能力定义了增殖vs B细胞样转录轴，该轴同时 underlie 转录组亚组和亚组内MRD反应。
在CASSIOPEIA试验257例接受四联诱导患者中验证：CCND1过表达和CD1/CD2亚组与MRD阳性相关，HP1亚组与MRD阴性相关；MRD相关差异表达基因和通路与MIDAS高度一致。
图1. MIDAS试验423例患者的分子分析
基因组不稳定表型与MRD：高FAG与MRD阴性相关，APOBEC高肿瘤与MRD阳性相关
WGS显示每例MM中位9,360个体细胞突变和26个结构变异（SV）。突变负荷在分子组间差异显著，MF组最高，因APOBEC特征高活性。SV数量也在组间差异，MF和MS组最高，CD1/CD2最低。染色体碎裂见于16%的MM，在MF和MS组更常见。染色体不稳定性以异常基因组分数（FAG）衡量，HP组最高，CD1/CD2最低。
较高FAG与MRD阴性强相关（FDR校正P=2.9×10⁻⁵），与HP1/HP3组高MRD阴性率一致。总体突变负荷与MRD无关，但SBS9和SBS18与MRD阴性相关。相反，APOBEC高肿瘤与MRD阳性相关（P=0.015），且与MF组高度重叠。
研究还发现罕见DNA修复缺陷病例：2例超突变肿瘤携带COSMIC ID2和SBS6特征，均检出MSH2致病性变异，提示MMRd。SV3特征与BRCA2改变状态显著相关，在MMRF队列中得到验证，提示MM中存在同源重组缺陷特征。
图2. 基因组不稳定表型及其与诱导后MRD的关联
WGS揭示MM新驱动基因：IRAK1和FOXR1
dN/dS分析揭示25个基因具有显著多于偶然预期的有害突变。罕见但显著复发性突变涉及组蛋白连接蛋白H1.4（14例，3.3%）和核糖体蛋白RPL10（10例，2.4%）。通过系统筛选与显著过表达相关的复发性SV断点基因，识别出2个潜在新癌基因：IRAK1和FOXR1。
IRAK1在15例MM（3.6%）中显示断点位于1.5 Mb内，伴显著过表达；在MMRF队列（779例）中验证，57例（7.3%）有IRAK1近端SV和一致过表达。FOXR1在6例MM（1.4%）中因染色体内倒位导致大量过表达，产生与DDX6和RPS25的融合转录本；在MMRF、MANHATTAN和CASSIOPEIA队列中识别出类似FOXR1改变。
驱动改变全景与MRD：整合突变、SV和局灶性CNA，每例样本中位识别出4个驱动改变。驱动改变数量在达到MRD阴性患者中显著更高（P=0.024）。就个体驱动基因而言，CCND1改变与MRD阴性率较低相关（P=5.6×10⁻⁷）；相反，MYC（P=2.5×10⁻⁴）、TXNDC5（P=0.0064）和BRAF（P=0.0098）改变与更高MRD阴性率相关。
图3. 新驱动事件及MRD关联
DNA甲基化分类：九个MRD相关亚组，表观遗传特征提供最强单组学预测
WGEM-seq数据分类识别出九个不同甲基化组（M1-M9），与基因组和转录组特征强相关。M1和M4分别与MF和MS转录组组对应；t(11;14)肿瘤聚集于M2和M3；超二倍体肿瘤分布于M6-M9。与正常浆细胞相比，所有MM亚组均存在广泛低甲基化，主要影响静息染色质、Polycomb抑制域和异染色质。
MRD分析显示甲基化组间差异显著（FDR校正P=6.7×10⁻¹¹）。M1-M3和MS与较低MRD阴性率相关；M4和超二倍体亚组M7、M9与较高MRD阴性相关。表观遗传特征提供了最强的单组学MRD预测。
驱动易位诱导局部和全基因组甲基化改变：所有易位癌基因因与Ig位点并置而发生广泛低甲基化。基因体甲基化与基因表达正相关（6,677个基因），M4（t(4;14)）组中KRAS信号基因同时过表达和基因体高甲基化。
多模态MRD预测模型：AUC 0.77，优于临床和基因组模型
全队列中20个特征与MRD状态显著相关，涵盖全局基因组特征、驱动改变、突变特征及转录组和甲基化亚组。将361例具有完整多组学数据和MRD状态的患者分为训练集（2/3，n=241）和留出验证集（1/3，n=120），采用嵌套交叉验证框架开发多模态随机森林模型。
在留出验证集中，多组学模型AUC为0.77（95% CI 0.67-0.85），预测MRD状态准确率71%，校准截距和斜率接近理想值。多组学模型区分度高于基于ISS、R-ISS和高危细胞遗传学异常的临床预测模型及基因组模型（P<0.015）；与转录组模型（AUC 0.73）和表观遗传模型（AUC 0.75）的AUC差异不显著，但校准指标更优。表达和甲基化簇、SBS9、SBS18、FAG以及CCND1和TXNDC5改变在重采样中一致被选中。
讨论与局限性
本研究证明整合多组学数据与MRD状态可阐明诱导治疗早期反应的分子决定因素。基因表达亚组显示显著不同的诱导后MRD率，B细胞样分化特征与较低MRD阴性率相关；表观遗传特征提供最强的单组学MRD预测。诱导后MRD阴性并不完全符合既定风险类别，MF和MS亚组传统上均被认为高危但MRD阴性率相反，提示既定风险类别可能不适用于当代四联诱导，或MRD主要捕捉早期治疗反应而非长期疾病控制。
局限性包括：单臂研究无对照组；多组学模型需外部验证以评估泛化性；MIDAS纵向设计将允许未来与PFS和OS关联。
总结
本研究建立了MIDAS多组学队列作为理解NDMM治疗反应分子决定因素的重要资源。基因组、转录组和表观遗传特征共同定义了与Isa-KRd诱导后早期MRD反应相关的分子景观，多模态整合模型优于临床和单一组学模型。这些发现为MM精准治疗提供了持久框架。
参考文献
Magrangeas F, Cherkaoui M, Bessonneau-Gaborit V, et al. Molecular predictors of minimal residual disease in newly diagnosed multiple myeloma treated with quadruplet therapy. Nat Med. 2026. doi:10.1038/s41591-026-04692-0.