[论文解读] Modeling the propagation of tumor fronts with shortest path and diffusion models -- implications for the definition of the clinical target volume

Thomas Bortfeld, Gregory Buti|arXiv (Cornell University)|Feb 26, 2022

Mathematical Biology Tumor Growth参考文献 31被引用 9

一句话总结

本文提出一种反应-扩散模型，作为放射肿瘤学中定义临床靶体积（CTV）的标准几何GTV-CTV边缘的机制性替代方案。通过利用扩散和最短路径（Eikonal）方程对肿瘤前沿传播进行建模，结果表明两种模型生成的波前形状相似，但扩散模型可通过扩散系数与增殖率之比实现可控的平滑性，从而为复杂各向异性组织（如大脑）中的CTV勾画提供更具生物学基础且可调节的框架。

ABSTRACT

Objective: The overarching objective is to make the definition of the clinical target volume (CTV) in radiation oncology less subjective and more scientifically based. The specific objective of this study is to investigate similarities and differences between two methods that model tumor spread beyond the visible gross tumor volume (GTV): 1. The shortest path model, which is the standard method of adding a geometric GTV-CTV margin, and 2. The reaction-diffusion model. Approach: These two models to capture the invisible tumor "fire front" are defined and compared in mathematical terms. The models are applied to example cases that represent tumor spread in non-uniform and anisotropic media with anatomical barriers. Main Results: The two seemingly disparate models bring forth traveling waves that can be associated with the front of tumor growth outward from the GTV. The shape of the fronts is similar for both models. Differences are seen in cases where the diffusive flow is reduced due to anatomical barriers, and in complex spatially non-uniform cases. The diffusion model generally leads to smoother fronts. The smoothness can be controlled with a parameter defined by the ratio of the diffusion coefficient and the proliferation rate. Significance: Defining the CTV has been described as the weakest link of the radiotherapy chain. There are many similarities in the mathematical description and the behavior of the common geometric GTV-CTV expansion method, and the definition of the CTV tumor front via the reaction-diffusion model. Its mechanistic basis and the controllable smoothness make the diffusion model an attractive alternative to the standard GTV-CTV margin model.

研究动机与目标

• 通过用机制性模型替代启发式几何边缘，减少临床靶体积（CTV）定义中的主观性。
• 比较标准几何GTV-CTV扩展（最短路径模型）与反应-扩散模型在肿瘤扩散方面的数学和行为相似性。
• 评估解剖屏障和组织各向异性对两种模型的影响，特别是在大脑等非均匀介质中。
• 证明反应-扩散模型相较于当前CTV勾画实践，是一种更具科学基础且可调控的替代方案。

提出的方法

• 利用度量张量表示组织各向异性和阻力，通过Eikonal方程（最短路径模型）模拟肿瘤前沿传播。
• 采用快速推进法（FMM）求解Eikonal方程，以高效、精确地计算从GTV出发的距离。
• 应用具有空间变化扩散张量D(r)和增殖率ρ的反应-扩散方程，模拟肿瘤细胞密度的演化。
• 将肿瘤前沿定义为肿瘤细胞密度图在某一阈值处的等值面（iso-surface）。
• 通过黎曼度量变换，正式建立各向异性介质中Eikonal模型与反应-扩散模型之间的关联。
• 在包含解剖屏障和方向性组织结构的几何测试案例中验证模型行为。

实验结果

研究问题

• RQ1 在非均匀和各向异性介质中，最短路径（Eikonal）模型与反应-扩散模型生成的肿瘤前沿形状有何异同？
• RQ2 解剖屏障在多大程度上改变了两种模型的传播模式？
• RQ3 反应-扩散模型能否生成比几何GTV-CTV边缘更平滑、更符合生物学特性的前沿？
• RQ4 扩散系数与增殖率之比如何影响反应-扩散模型中前沿的平滑性？
• RQ5 反应-扩散模型能否作为当前标准几何CTV扩展的机制性更优替代方案？

主要发现

• 在均匀介质中，最短路径（Eikonal）模型与反应-扩散模型均生成从GTV向外传播的行进波前，整体形状相似。
• 在存在解剖屏障的情况下，反应-扩散模型产生的波前更平滑、渐变更缓和，而Eikonal模型则表现出更尖锐的过渡。
• 反应-扩散模型波前的平滑性可直接通过扩散系数与增殖率之比进行控制，从而实现可调的波前形态。
• 在白质纤维束等复杂各向异性环境中，反应-扩散模型更能准确捕捉方向性扩散，与已知生物学通路一致。
• 数学形式化证明了通过黎曼度量变换，Eikonal方程与反应-扩散模型在特定条件下具有直接等价性。
• 反应-扩散模型为CTV定义提供了更具生理学合理性且可调控的框架，尤其在神经肿瘤学领域具有优势。