多尺度DCE-BOLD MRI揭示胶质瘤IDH基因分型的血管解耦联模式：一项前瞻性基于生境的研究

《Neurotherapeutics》：Multiscale DCE-BOLD MRI reveals vascular uncoupling patterns for IDH genotyping in adult glioma: A prospective habitat-based study

　　

异柠檬酸脱氢酶（IDH）野生型胶质瘤表现出异质性血管功能障碍，但全肿瘤分析可能掩盖血管通透性和功能储备的区域变异。研究人员开发了一种多尺度生境框架，整合了动态对比增强磁共振成像（DCE-MRI）衍生的容积转移常数（Ktrans）和静息态血氧水平依赖功能磁共振成像（rs-BOLD-fMRI）衍生的相对脑血管反应性（rCVR），并纳入10毫米瘤周边缘区域，用于无创IDH基因分型。这项前瞻性单中心队列研究纳入132例患者（52例IDH突变型和80例IDH野生型），按时间顺序分为训练队列（n=92）和内部测试队列（n=40）。高斯混合模型将连续的联合Ktrans-rCVR分布划分为三个重叠的影像学定义的瘤内亚区；解剖学定义的瘤周区域被单独分析。Habitat-Peri10mm模型在内部测试中取得了最高的受试者工作特征曲线下面积（AUC）点估计值0.914（95%置信区间，0.831–0.997），而Habitat模型和Whole-tumor模型的AUC分别为0.854和0.753。成对DeLong检验未达到统计学显著性（p值分别为0.173和0.054）。来自亚区2的Ktrans一阶熵特征在SHapley Additive exPlanations（SHAP）分析中排名最高，并与Ki-67（r=0.57）和微血管密度（r=0.27）相关。生境衍生的瘤内异质性评分在IDH野生型肿瘤中高于IDH突变型肿瘤（0.79 vs. 0.40，p<0.001；AUC=0.824），尽管级别不平衡限制了IDH特异性解释。这种DCE-BOLD多尺度生境框架表征了与IDH相关的血管异质性模式，并显示出无创IDH基因分型的潜力。

研究背景、现存问题与研究目的

胶质瘤是成人中最常见且最具侵袭性的原发性恶性脑肿瘤，其治疗抵抗和不良预后根源于高度异质且动态的肿瘤微环境（TME），其中异常的血管结构和区域性缺氧是恶性进展的关键驱动因素。异柠檬酸脱氢酶（IDH）突变状态与微血管增殖和坏死等血管病理特征密切相关。然而，当前的IDH基因分型依赖于具有内在空间偏倚的有创手术取样，可能遗漏最具侵袭性的亚克隆，而常规的结构性磁共振成像（MRI）对此功能性异质性的信息提供有限。因此，亟需一种能够解析肿瘤血管功能空间复杂性的无创成像策略。

为应对这一挑战，研究人员开展了一项前瞻性研究，旨在开发并评估一种结合动态对比增强磁共振成像（DCE-MRI）衍生的容积转移常数（K^trans）和静息态血氧水平依赖功能磁共振成像（rs-BOLD-fMRI）衍生的相对脑血管反应性（rCVR）的多尺度生境成像框架。该框架通过无监督聚类方法划分瘤内亚区，并扩展至解剖学定义的10毫米瘤周区域（Peri10mm），以无创方式表征与IDH状态相关的血管通透性与功能储备的空间耦合模式，并探索其与组织病理学标志物及预后的关联。

主要关键技术方法

1. 1.
   前瞻性单中心队列：研究纳入132例弥漫性胶质瘤患者（52例IDH突变型，80例IDH野生型），来源于中国医科大学附属第一医院（2024年6月至2026年1月），并按时间顺序分为训练队列（n=92）和内部测试队列（n=40）。
2. 2.
   多参数MRI采集与处理：使用3.0T扫描仪获取数据。通过标准流程从rs-BOLD-fMRI数据生成rCVR图，并使用Tofts模型从DCE-MRI数据定量计算药代动力学参数图（包括K^trans）。
3. 3.
   生境分析与影像组学建模：应用高斯混合模型（GMM）对联合的K^trans-rCVR图谱进行无监督聚类，划分出三个影像学定义的瘤内亚区。随后，利用PyRadiomics从全肿瘤、各亚区和瘤周区域提取影像组学特征，并通过递归特征消除（RFE）等方法进行特征选择和模型构建，采用支持向量机（SVM）、线性判别分析（LDA）和逻辑回归（LR）等分类器建立预测模型。

研究结果

• •
  人口统计学特征：IDH野生型患者年龄更大，更常伴有瘤周水肿，且Ki-67增殖指数和微血管密度（MVD）显著高于IDH突变型患者。所有IDH野生型肿瘤均为WHO 4级（胶质母细胞瘤），而IDH突变型肿瘤涵盖WHO 2-4级。
• •
  空间生境分析：三组分GMM解决方案具有最低的赤池信息量准则（AIC）值，平均轮廓得分为0.115，表明亚区在联合特征空间中存在重叠。生境分析揭示了三种不同的瘤内亚区，并根据其相对通透性-反应性特征被描述性地标记为：低通透性/低储备（亚区1）、中等通透性/高储备（亚区2）和高通透性/低储备（亚区3）。
• •
  影像组学模型性能：在内部测试队列中，整合了瘤内亚区和瘤周区域的Habitat-Peri10mm模型取得了最高的受试者工作特征曲线下面积（AUC）点估计值0.914，优于单独的Habitat模型（0.854）、全肿瘤模型（0.753）和瘤周模型（0.776）。然而，成对DeLong检验未显示统计显著性差异。决策曲线分析（DCA）表明所有模型均具有临床净获益。
• •
  模型可解释性与生物学关联：SHAP分析显示，来自亚区2的K^trans一阶熵特征对模型预测贡献最大。该特征与Ki-67（r=0.57）和MVD（r=0.27）呈正相关。整合的影像组学评分（Radscore）与Ki-67的相关性最强（r=0.61）。此外，亚区3的形状特征LeastAxisLength在有和无瘤周水肿的肿瘤间存在显著差异。
• •
  瘤内异质性（ITH）评分：基于主成分分析（PCA）计算的ITH评分在IDH野生型肿瘤中显著高于IDH突变型肿瘤（0.79 vs. 0.40，p<0.001），其对IDH基因分型的AUC为0.824。进一步的探索性分析显示，无论表皮生长因子受体（EGFR）扩增状态如何，IDH野生型患者的ITH评分始终显著更高。

讨论与结论

讨论部分指出，该研究的独特贡献在于联合表征了与通透性和反应性相关的成像变异。与既往结合灌注与弥散或其他参数的生境研究相比，K^trans-rCVR配对纳入了源自静息态BOLD成像的血管功能储备评估，从而能评估高通透性区域是否同时保留或丧失血管反应性，为描述通透性-反应性耦合模式提供了生理学基础。研究还讨论了各影像学定义亚区的潜在生物学意义，例如亚区2的异质性与增殖活性和血管特征相关，而亚区3的高通透性/低储备模式则符合血管解耦联的概念。同时，研究者也坦诚地指出了研究的局限性，包括单中心设计导致的过拟合风险、IDH状态与WHO级别的混杂效应、亚区重叠度高、组织病理学测量缺乏空间匹配以及随访数据不完整等问题。

研究结论：总之，本研究提出了一种基于DCE-BOLD的多尺度生境成像框架，能够无创地表征胶质瘤的瘤内和瘤周功能性异质性。通过整合K^trans和rCVR，该框架捕捉了三个GMM定义的瘤内亚区和解剖学定义的瘤周区域内通透性-反应性耦合的空间变异，并为IDH基因分型提供了一个较高的内部测试AUC点估计值。ITH评分进一步表征了与IDH状态相关的成像异质性。总体而言，这些发现为将影像学定义的血管模式与肿瘤病理生理学联系起来提供了一个基于生物学的框架。未来的前瞻性多中心研究，结合外部验证和成像-组织病理学配准，对于确认这些发现并推进临床转化是必要的。