在动物模型中监测抗肿瘤治疗效果的新方法：结合功能MRI和纹理分析

背景

功能性磁共振成像（fMRI）发展非常迅速，因为它提供了无创性和精确的成像，特别是其区分组织特征的能力。此外，利用病变的特征，fMRI可实时、无损地测量体内的病理过程，以便早期诊断和治疗评估。多b值扩散加权成像（DWI）和动态增强MRI（DCE-MRI）这两种新型fMRI扫描技术可以潜在地检测乳腺癌等主要疾病。总的来说，DCE-MRI在检测乳腺癌方面显示出很高的敏感性（89-100%），DWI在预测适当的治疗方案和监测治疗反应方面显示出了实用性[1]. 肿瘤内血管异质性对肿瘤治疗至关重要。因此，抗血管生成治疗被认为是一种非常有前途的预防肿瘤生长和转移的新策略。这两种功能性MRI技术能够测量微血管结构并反映其通透性[2]. DCE-MRI的几个定性和半定量参数，从时间-强度曲线的简单半定量检查到更复杂的示踪动力学模型，可以提供肿瘤内血管通透性的信息[三]. 此外，表观扩散系数（ADC）的值基于肿瘤组织的相对信号强度变化，在多b值DWI中随着b值的增加而增加，可以提供细胞水平的显微结构信息。ADC值的变化与组织和细胞密度呈负相关[4,5]. 因此，这两种成像方法有可能用于监测和评估治疗早期抗血管生成治疗的疗效。

最近的临床研究表明，贝伐单抗（bevacizumab）是一种基因工程人源化单克隆抗体，由于其抗血管内皮生长因子（VEGF）的活性，在治疗各种肿瘤方面非常有效。贝伐单抗可以特异性地与VEGF结合，阻止VEGF与VEGFR的结合，从而抑制新血管的形成，并以低毒性抑制肿瘤生长[6]. 作为对照，另一种常用的化疗药物紫杉醇可以与β-微管蛋白结合并稳定微管，从而抑制细胞有丝分裂和细胞增殖[7]. 如上所述，一种有前途的方法是联合使用多b值DWI和DCE-MRI来评估贝伐单抗与紫杉醇相比的抗血管生成活性。

为了确保我们研究的准确性，我们采用了另一种新技术——纹理分析，来分析和验证成像结果。纹理分析作为一种新的肿瘤成像生物标志物，可以量化肿瘤的区域异质性，这是公认的恶性肿瘤特征，与侵袭性生物学、不良预后和治疗耐药性有关[8]. 因此，该图像处理算法可以通过评估纹理粗糙度的分布来扫描细微的肿瘤内异常。重要的纹理参数，包括平均强度、灰度直方图分布的标准差、熵（灰度分布的不规则性）、偏度（直方图的不对称性）和峰度（直方图的平坦性），可以反映从解剖结构到生物功能的各种信息[9]. 先前的研究表明，与其他成像和生物学参数相比，粗糙纹理特征可能反映CD34定义的潜在血管系统[10]. 根据这项研究，对功能性MRI检查结果进行纹理分析并评估结果之间的相关性是有价值的。

方法

结果

讨论

在本研究中，我们旨在开发一种实用的方法来评估早期抗肿瘤治疗的疗效。以前的研究表明，血管生成可以为肿瘤提供营养和氧气，因此在肿瘤进展中起着至关重要的作用[16]. 肿瘤在血管受累时呈指数级增长，但在无血管环境中缓慢线性增长[17]. 因此，抗血管生成在肿瘤治疗中具有不可替代的作用，以肿瘤血管生成为靶点的抗肿瘤药物已成为近年来的研究热点。贝伐单抗是FDA批准的第一种抑制肿瘤血管生成的药物，其在阻断VEGF诱导的血管生成方面具有很高的亲和力，可以诱导内皮细胞的增殖和迁移，增加微血管的通透性[18]. 通常，评估药物是否成功抑制肿瘤血管生成的金标准是MVD计数。然而，在临床实践中，通过计算微血管密度，几乎不可能连续地从患者体内清除肿瘤组织来实时观察抗血管生成治疗的疗效。令人鼓舞的是，我们的研究证实，这个问题可以通过一种新的多参数融合分析来解决。

在这项临床前研究中，我们发现许多重要的成像参数对不同的治疗很敏感。在添加贝伐单抗后，联合治疗组的功能MRI和纹理分析的变化非常显著，导致肿瘤体积与其他组相比存在差异。DWI在反映组织的微观结构（高b值）和血液灌注状态（低b值）方面具有很大的优势，尤其是其关键参数ADC[19,20]. 因此，如果治疗有效，细胞完整性将被破坏，那么ADC_缓慢的从高b值DWI中得出的值将因水扩散的增强而上升[21]这得到了我们的研究结果的支持。随着肿瘤中央部位坏死的发生，ADC_缓慢的未经任何治疗，对照组的数值略有增加。然而，当贝伐单抗阻断血管生成时，肿瘤生长所需的营养素将不足，由此导致的细胞密度降低将导致ADC显著增加_缓慢的值。同时，实验数据表明，紫杉醇对细胞有丝分裂的抑制导致的细胞密度降低不如贝伐单抗，但ADC增加_缓慢的与对照组相似。此外，f值直接评估了血液灌注，并显示出各组之间低b值DWI的显著差异。结果与ADC对比_缓慢的抗血管生成治疗在治疗开始后15天导致f值显著降低，但其他两组则表现出相反的趋势。此外，f值的变化与MVD密切相关，但ADC的变化_缓慢的与微血管计数呈显著负相关。DWI参数和组织学结果的非常有意义的相关性与早期研究完全一致，研究表明DWI可用于监测血管靶向药物的早期治疗效果[22].

DCE-MRI是组织血液灌注无创评估的最常用技术，是通过跟踪Gd-DTPA的药代动力学来监测各种治疗效果的有效方法[23]. 反映血管通透性、血流量和血容量的最常用参数是K^反式.与其他参数相结合，如K_{电动自行车}，K^反式可以在一定程度上反映肿瘤血管生成的程度[24,25]. 我们的研究表明，高钾^反式对照组随着肿瘤的生长而出现数值。这一发现与两个治疗组（如K^反式价值不断下降。K的增加^反式数值表明肿瘤血液灌注和高毛细血管通透性的增加为肿瘤生长提供了更多的营养，并最终加速了肿瘤细胞的增殖。在实验后期，对照组的皮下肿瘤体积明显大于其他两组，这为K提供了很好的验证^反式此外，紫杉醇和贝伐单抗的不同机制导致两个治疗组的下降趋势显著不同。紫杉醇通过抑制微管系统具有明确的抗肿瘤作用。然而，一些学者证实贝伐单抗可以改善紫杉醇的给药和疗效[26]. 贝伐单抗抑制血管生成和血管通透性可确保紫杉醇的浓度。体积的显著变化，K^反式与单纯紫杉醇组和对照组相比，联合组的其他影像学参数可能发生变化，因为治疗持续时间不够长，导致紫杉醇单药组与对照组之间存在明显差异。令人鼓舞的是，组织学结果与DCE-MRI一致。MVD计数与K^反式通过贝伐单抗治疗，联合用药组VEGF的表达降低。近年来，K受到了越来越多的关注_{电动自行车}之前的研究表明，K的高基线值_{电动自行车}对应于血浆和血管外细胞间隙（EES）之间药物的高交换分数，表明潜在的优越治疗效果[27]. 最有可能的是，个体差异、肿瘤细胞坏死等因素导致造影剂残留在间质空间并导致血管外细胞外渗透量的误差，最终导致K缺乏显著变化_{电动自行车}在我们的研究中。另一方面，K_{电动自行车}也受到V的显著影响_{e（电子）}根据Tofts的说法，这可能由细胞密度、囊性变性和组织反应等因素决定[28]，V_{e（电子）}不是一个非常稳定的因素，因为它经常受到病变周围水肿的影响。然而，当我们分析DCE-MRI和DWI之间的相关性时，我们发现K_{电动自行车}与ADC呈负相关_快速的从低b值DWI中提取。由于双指数模型中的ADC值主要反映肿瘤密度特征，肿瘤密度的增加必然会影响造影剂从EES返回血浆的速率。因此，从上述分析可以得出结论，多b值DWI和DCE在评估血管生成功能和肿瘤灌注方面是互补的。

虽然多b值DWI和DCE-MRI为监测肿瘤生长和肿瘤治疗效果提供了大量信息，但这两种成像技术可能受到许多因素的影响，例如肿瘤组织的不均匀性、，皮下肿瘤模型和动物在成像过程中的运动产生的伪影[29]. 此外，临床图像在反映病变的细胞和分子特征方面也有一些局限性，如细胞增殖和代谢、坏死和缺氧[30]. 最近，越来越多的研究试图通过纹理分析来阐明医学图像异质性的测量，纹理分析是一种二阶统计技术，其参数来源于局部特征的分布，可以更好地描述组织特征、图像分割、，并预测治疗反应和生存率[31,32]. 因此，该潜在工具的主要优点是可以最大限度地利用临床图像中的信息，而无需额外采集[9]. 在我们的研究中必须充分利用这一优势。通过测量未增强的T2加权MRI，我们发现所有移植瘤荷瘤小鼠在治疗前都处于相同的状态，但随着治疗和各种处理，三组在所有SSF下的熵值显著增加。熵表示ROI中像素的无序程度，其值越高，组织的无序度越大。之前的一份出版物显示，严重程度与纹理粗糙度有关，纹理粗糙度与葡萄糖摄取测量值相关（从FDG-PET、，第页 = 0.51,对 = 0.03) [33]. 因此，很明显，葡萄糖代谢的增加使这4个T1同种异体移植瘤的生长速度增加，这与所有小鼠肿瘤的增大相一致。根据Ng等人[34]，肿瘤组织的异质性随着生长而增加。根据Ganeshan等人[10]和Henriksson等人[30]肿瘤内部图像异质性的增加可能与区域肿瘤细胞数、增殖、缺氧、血管生成和坏死的差异有关。因此，通过贝伐单抗和紫杉醇联合治疗的抗血管生成和抑制细胞有丝分裂的作用，肿瘤的微观结构，包括细胞、细胞外基质和微血管会受到干扰，在细胞和分子水平上产生一系列变化，这些变化太细微，无法使用传统的成像诊断技术进行检测。持续的变化最终导致病灶内像素平均值的显著差异（平均强度，对 < 0.05），平均值（SD，对 < 0.03). 由于缺乏强有力的有效化疗，其他两组治疗后均未出现明显变化。在一项核医学研究中，学者们发现，水分分布更不均匀的肿瘤（即SD和阳性像素平均值更高，MPP）更容易发生糖酵解[35]. 我们的经验证据也支持这一结论。当贝伐单抗阻断血管生成时，组织灌注的减少限制了肿瘤的氧气供应，导致与治疗前相比，对糖酵解能量的依赖性显著（P_MP公司对 < 0.05). 支持这一说法的另一项发现是，平均值、SD和MPP的变化都发生在中、粗纹理尺度上，这更倾向于反映基于Chowdhury等人的调查的基因组学分析中的生物特征[35]. 此外，上述分析适用于不同组之间的比较。细胞和分子水平上的差异，如单一治疗和联合治疗引起的抗增殖、缺氧、血管生成和坏死，最终导致解剖结构的差异（在精细和中等纹理尺度下），体现了像素平均值（平均值、，对 < 0.05）和正像素（MPP，对 < 0.05）。如上所述，在传统成像参数中可以观察到这些主要结构变化。在我们的研究中，纹理分析不是独立的；它与功能磁共振成像密切相关。熵与ADC显著负相关_缓慢的(第页 = − 0.503,对 < 0.05). 较高的熵表示非均质性增加，这在一定程度上表示水扩散受到限制（较低的表观扩散系数）。令人惊讶的是，MPP值的增加与ADC显著正相关_缓慢的(第页 = 0.603,对 = 0.013），可能是因为糖酵解环境（较高的MPP）产生的代谢物增加了细胞膜的渗透性并促进了水分子的扩散。然而，需要进一步确认。

诚然，我们的研究存在一些局限性。6周龄裸鼠的脆弱性和其他因素导致了实验期间的高死亡率；因此，在有限数量的小鼠中建立了动物肿瘤模型。此外，皮下肿瘤模型中气-软组织边界处DWI的敏感性伪影[29]动物在成像过程中的运动，以及植入的肿瘤比原发肿瘤更均质的事实导致了不可避免的系统误差。在进一步的研究中，我们将努力克服这些局限性，探索更加多样化的多模态融合成像方法。

结论

本研究显示，在临床前乳腺癌模型中，使用多b值DWI、DCE-MRI和纹理分析成功监测抗血管生成治疗的早期阶段。通过多参数的集成，DWI、DCE和TexRad提供了从生物学特征到解剖结构的全面而有价值的信息。更令人鼓舞的是，治疗前后关键参数的显著变化都与组织学结果有良好的相关性和一致性。这三种成像和分析技术相辅相成，可能作为非侵入性生物标记物，用于指导治疗算法和监测未来临床试验中抗血管生成治疗的早期反应。

作者的贡献

MM、JL和HX设计了实验；MM进行了实验；MM分析数据并撰写手稿；QW、JL、YL和TS提供了材料和分析工具，并进行了统计分析；HX和ZJ设计并监督了整个研究；ZJ参与了讨论、数据解释和手稿修订。所有作者审查并批准了最终稿。