Name: 跟踪乳腺建筑特色与检测乳腺癌的磁共振扩散张量成像 (视频) | JoVE Journal
Uploaded: 2014-12-15T15:00:00.000+00:00
Duration: 15 min 48 s
Description: 学习使用磁共振扩散张量成像技术追踪乳腺微结构特征并检测乳腺癌的非侵入性方法。该技术通过分析水分子扩散系数和各项异性指数生成彩色参数图，无需造影剂即可识别恶性病变区域，为乳腺癌早期诊断和疗效监测提供安全可靠的影像学工具。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

我们感谢Nachum Stern先生和Fanny Attar女士提供的出色技术协助。H.D. 担任乳腺癌研究的Fred和Andrea Fallek讲席教授。

注意：所有实验方案均经以色列克法萨巴梅尔医疗中心机构审查委员会批准，并已从所有受试者获得签署的知情同意书。
1. 患者在磁共振成像扫描仪中的准备与体位摆放
<ol>
	<li>在左臂或右臂插入静脉导管，用于注射对比剂。</li>
	<li>确保患者取俯卧位，双侧乳房自然下垂，分别置于乳腺线圈的双侧开口中，如图5所示。</li>
	<li>确保患者头部置于枕头上，以保证头颈部舒适，双臂上举置于头侧，如图5所示。</li>
	<li>确保每侧乳房均处于中心位置，松散下垂并尽可能深入各自线圈开口内，如图5所示。</li>
</ol>
<img alt="用于成像研究的带有线圈和轮廓图的 MRI 装置示意图。" src="/files/ftp_upload/52048/52048fig5highres.jpg" data-altupdatedat="1787630408000" data-alt="图 5"> 
图 5：女性乳腺 MRI 扫描时的体位。注意俯卧位以及双侧乳房自然悬垂于乳腺线圈的双侧开口中。
<ol start="5">
	<li>将自动注射泵连接至静脉导管。</li>
</ol>
2. MRI 扫描
<ol>
	<li>使用MRI扫描仪获取定位图像，以定位双侧乳房，并确定视野（FOV）、扫描层数和层厚，确保扫描范围完全覆盖双侧乳房，并延伸至腋窝及胸壁。</li>
	<li>在后续使用MRI扫描仪进行的所有扫描序列中，每个病例均需保持相同的FOV、层数和层厚。例如，采集60个轴向切片，读取方向和相位方向的FOV为360 mm，层厚根据乳房大小在2至3 mm之间调整，相位编码方向为L &gt;&gt; R或R &gt;&gt; L。根据各序列所用的采集矩阵尺寸，调整后续序列中的面内分辨率，具体参见表1。</li>
</ol>
<table border="0" cellpadding="0" cellspacing="0">
	<tbody>
		<tr>
			<td>实验参数</td>
			<td>2D T2 加权</td>
			<td>2D DTI</td>
			<td>2D GE 场图</td>
			<td>3D DCE</td>
		</tr>
		<tr>
			<td>重复时间 (TR) msec</td>
			<td>5,500</td>
			<td>10,800</td>
			<td>669</td>
			<td>6.8</td>
		</tr>
		<tr>
			<td>回波时间 (TE) msec</td>
			<td>122</td>
			<td>120</td>
			<td>4.92</td>
			<td>2.49</td>
		</tr>
		<tr>
			<td>TE(2)</td>
			<td></td>
			<td></td>
			<td>7.38</td>
			<td></td>
		</tr>
		<tr>
			<td>翻转角，度</td>
			<td></td>
			<td></td>
			<td>60</td>
			<td>18</td>
		</tr>
		<tr>
			<td>采集矩阵</td>
			<td>640 x 428</td>
			<td>192 x 192</td>
			<td>64 x 192</td>
			<td>428 x 428</td>
		</tr>
		<tr>
			<td>平均次数</td>
			<td>1</td>
			<td>1</td>
			<td>1</td>
			<td>1</td>
		</tr>
		<tr>
			<td>拼接次数</td>
			<td>4</td>
			<td>1</td>
			<td>1</td>
			<td>1</td>
		</tr>
		<tr>
			<td>涡轮因子</td>
			<td>21</td>
			<td></td>
			<td></td>
			<td></td>
		</tr>
		<tr>
			<td>每层的回波链数量</td>
			<td>11</td>
			<td></td>
			<td></td>
			<td></td>
		</tr>
		<tr>
			<td>回波间隔，msec</td>
			<td>10</td>
			<td>0.79</td>
			<td></td>
			<td></td>
		</tr>
		<tr>
			<td>带宽，Hz/像素</td>
			<td>300</td>
			<td>1860</td>
			<td>1532</td>
			<td>560</td>
		</tr>
		<tr>
			<td>平面内扫描分辨率</td>
			<td>0.56 x 0.84</td>
			<td>1.9 x 1.9</td>
			<td></td>
			<td>0.8 x 0.8</td>
		</tr>
		<tr>
			<td>平面内图像分辨率</td>
			<td>0.56 x 0.56</td>
			<td>1.9 x 1.9</td>
			<td>1.9 x 1.9</td>
			<td>0.8 x 0.8</td>
		</tr>
		<tr>
			<td>采集时间，分:秒</td>
			<td>4:26</td>
			<td>6:09</td>
			<td>1:28</td>
			<td>1:06 (x9)</td>
		</tr>
		<tr>
			<td>扩散梯度数量</td>
			<td></td>
			<td>30</td>
			<td></td>
			<td></td>
		</tr>
		<tr>
			<td>脂肪抑制</td>
			<td>无</td>
			<td>FAT-SAT / SPAIR1</td>
			<td></td>
			<td>无</td>
		</tr>
		<tr>
			<td>b 值，sec/mm2</td>
			<td></td>
			<td>0, 700</td>
			<td></td>
			<td></td>
		</tr>
	</tbody>
</table>
表1：本研究中所用序列的实验参数。 1FAT-SAT：脂肪饱和。SPAIR：频谱衰减反转恢复。
<ol start="3">
	<li>使用 MRI 扫描仪的用户界面计算机定位包含双侧乳腺及腋窝的区域，并定义匀场区域。采用迭代匀场策略以优化磁场29。
	<ol>
		<li>在用户界面计算机屏幕上显示质子谱，通过将频率中心对准水共振频率，然后对准脂肪共振频率来调整匀场。重复匀场过程，直至在质子谱中脂肪与水信号达到最佳分离，并使信号强度和峰形最优化。确认照射频率已对准水共振频率。</li>
	</ol>
	</li>
	<li>通过用户界面计算机应用一种高空间分辨率、无脂肪饱和的二维横断面T2加权快速自旋回波多层序列，采用广义自校准部分并行采集技术（GRAPPA），加速因子为2。根据表1中“2D T2加权”列的数值设定屏幕上弹出的各项参数。</li>
	<li>通过用户界面计算机应用一种二维横断面、脂肪抑制的自旋回波扩散张量成像（DTI）序列，采用两次重聚焦的回波平面成像序列30，并使用GRAPPA技术，加速因子为2。根据表1中“2D DTI”列的数值设定屏幕上弹出的各项参数。</li>
	<li>通过用户界面计算机应用一个场图映射序列，用于校正回波平面成像（EPI）中的几何畸变，并按照Jezzard和Balaban所述方法获取相位差图像31。确保该序列包含采集两个不同同相位回波时间的二维横断面梯度回波（GE）图像，且相位编码方向与DTI序列（步骤2.5）一致。根据表1中“2D GE场图映射”列的数值设定屏幕上弹出的各项参数。</li>
	<li>通过用户界面计算机应用一种动态对比增强扫描方案，采用无脂肪抑制的三维快速梯度回波序列，参数根据三点时间法（3TP）进行优化32。根据表1中“3D DCE”列的数值设定屏幕上弹出的各项参数。
	<ol>
		<li>按照步骤2.7中的方案记录两次初始的对比前图像。</li>
	</ol>
	</li>
	<li>在第二次对比前扫描结束前15秒，以0.1 mmol/kg体重剂量，通过自动注射泵以2 ml/sec速率注射对比剂钆喷酸葡胺，随后以2 ml/sec速率推注20 ml生理盐水冲洗。
	<ol>
		<li>在对比剂注射后54、120、186、252、318、384和450秒，按照步骤2.7中的方案连续采集七组三维数据集。</li>
	</ol>
	</li>
	<li>通过用户界面计算机应用一种高空间分辨率、脂肪抑制的二维横断面T2加权快速自旋回波多层序列，采用GRAPPA技术，加速因子为2。实验参数应与步骤2.4中所述T2加权成像所用参数相似。 
	注：MRI扫描方案各步骤的时间流程总结见图6。</li>
</ol>
<img alt="MRI序列时间线；示意图；包含定位像、T2、DTI、场图、DCE等过程；共30.5分钟。" src="/files/ftp_upload/52048/52048fig6highres.jpg" data-altupdatedat="1787630408000" data-alt="图6"> 
图6：临床MRI检查流程的时间进程。
3. 图像处理
<ol>
	<li>将通过所有扫描序列（见步骤 2.4 - 2.9）获得的完整数据集传输至远程工作站或个人计算机，并使用专门用于分析乳腺扩散张量成像（DTI）和动态对比增强磁共振成像（DCE-MRI）的软件程序处理全部图像数据集。
	<ol>
		<li>此处，使用基于 C++ 编写的自研软件包执行下述步骤，该软件包已在 MATLAB19 编程环境中完成验证。通常情况下，商用扫描仪提供用于脑部 DTI 实验的图像处理工具，这些工具可被调整用于乳腺成像，但可能不具备本自研软件包中的全部处理功能</li>
	</ol>
	</li>
	<li>按照图 7所示流程图，使用专用的 DTI 图像处理软件对 DTI 数据集进行图像处理。</li>
</ol>
<img alt="DTI的MRI流程图；过程包括数据加载、拟合、PCA、参数图、追踪。" src="/files/ftp_upload/52048/52048fig7highres.jpg" data-altupdatedat="1787630408000" data-alt="Figure 7"> 
图7：DTI处理流程的步骤示意图。
<ol style="margin-left: 40px;">
	<li>使用 DTI 软件程序评估乳腺外部及剩余组织中的噪声水平，在 3 至 4 个感兴趣区域（每个区域约 1 cm2）内测量，并确定最大噪声水平。</li>
	<li>使用 DTI 软件程序在所有高于最大噪声水平的像素中，对全部 60 个切片逐像素计算对称张量 D 的六个扩散系数 Dxx、Dyy、Dzz、Dxy、Dxz、Dyz（如下所述33，见图8） 
	<img alt="扩散张量矩阵公式；描述三维空间中的扩散；示意图；用于 MRI 分析。" fo:content-width="2in" src="/files/ftp_upload/52048/52048eq1.jpg" data-altupdatedat="1787630408000" data-alt="方程 1"> 
	其中 Dxy=Dyx、Dzy=Dyz、Dzx=Dxz。
	<ol>
		<li>为计算张量分量，使用基于 Stejskal-Tanner 方程的非线性回归拟合程序：Si,j(b)=Si(0)exp(-bD)，其中 Si,j(b) 表示在方向 j 施加扩散梯度时像素 i 的信号强度，Si(0) 表示同一像素在 b = 0 mm2·s 时的信号强度。该拟合通过改进的 Levenberg-Marquardt 算法34，最小化每个扩散梯度方向上六个扩散张量变量分量对应的 30 个非线性函数的平方和。此步骤由免费软件 cminpack 完成（<a href="http://devernay.free.fr/hacks/cminpack/index.html" target="_blank">http://devernay.free.fr/hacks/cminpack/index.html</a>）。</li>
	</ol>
	</li>
</ol>
<img alt="弥散张量成像示意图、非线性回归、对称张量方程、数据分析。" src="/files/ftp_upload/52048/52048fig8highres.jpg" data-altupdatedat="1787630408000" data-alt="图8"> 
图8：计算对称扩散张量的步骤示意图。
<ol start="3" style="margin-left: 40px;">
	<li>使用 DTI 软件通过对每个像素应用主成分分析（PCA）35 来对称对角化扩散张量。 
	注意：PCA 是一种常用的方法，用于降低数据的维度，并找到一种线性变换，将相关的像素参数映射到一组新的不相关且正交的坐标轴系统中。该过程为每个像素生成三个特征向量（ν1, ν2, ν3），定义了与组织扩散框架一致的椭球形三个正交轴上的扩散方向，以及相应的三个按从高到低排列的扩散特征值，用于确定三个方向上的扩散系数（λ1, λ2, λ3）（见 图 9 和图 10）。</li>
</ol>
<img alt="使用矩阵 D、特征值、特征向量进行主成分分析并建模为椭球体的示意图。" src="/files/ftp_upload/52048/52048fig9highres.jpg" data-altupdatedat="1787630408000" data-alt="图 9"> 
图 9：计算步骤的示意图，展示获得与组织扩散框架一致的特征值和特征向量，并将其建模为椭球体的过程。
<img alt="显示主成分分析、特征向量和特征值的扩散张量成像示意图；显示癌细胞的影响。" src="/files/ftp_upload/52048/52048fig10highres.jpg" data-altupdatedat="1787630408000" data-alt="图 10"> 
图 10：计算正常导管组织和癌组织中扩散特征向量与特征值的主要处理步骤示意图。
<ol start="4" style="margin-left: 40px;">
	<li>使用 DTI 软件程序计算每个像素的表观扩散系数（ADC），其定义为三个特征值的平均值 &lt; λ &gt; = (λ1 + λ2 +λ3) / 3.</li>
	<li>使用 DTI 软件程序计算每个像素的最大绝对各向异性指数，该指数定义为差值（λ1 – λ3).</li>
	<li>使用 DTI 软件程序，根据以下公式计算每个像素的分数各向异性（FA）指数，其取值范围为 0（各向同性扩散）至 1（单方向自由扩散）： 
	<img alt="Fractional Anisotropy equation; diagram; diffusion tensor imaging; tensor calculation." fo:content-width="3in" src="/files/ftp_upload/52048/52048eq2.jpg" width="300" data-altupdatedat="1787630408000" data-alt="Equation 2"></li>
	<li>使用 DTI 软件程序在双侧乳腺的每个切片中逐像素显示矢量图，以展示主特征向量 ν 的方向1，以及一张用三种颜色表示 ν 主要方向的彩色编码图1 并将其叠加在同一层面的T2加权图像上（图9）。此步骤使用免费软件 gtk 完成（<a href="http://gtk.org" target="_blank">http://gtk.org</a>用于图形制作。</li>
	<li>使用DTI软件程序构建扩散张量成像参数图，以逐像素方式显示双侧乳腺各层面中所有扩散张量参数λ的数值1，λ2, λ3, ADC, λ1-λ3，以及FA，并将其叠加到同一层面的T2加权图像上（图11).</li>
</ol>
<img alt="乳腺MRI扩散张量成像、T2 MRI、v1向量图、λ1-λ3、ADC、FA图用于分析。" src="/files/ftp_upload/52048/52048fig11highres.jpg" width="700" data-altupdatedat="1787630408000" data-alt="图11"> 
图11：一名24岁健康志愿者双侧乳腺同一中央切面的参数化扩散张量成像（DTI）图谱叠加于该切面T2加权图像上的结果。 方向图（Dir.）中，红色表示从左到右，绿色表示从前到后，蓝色表示从头到脚的方向。注意，向量图中的向量以白色线条标示。所有扩散系数及最大各向异性指数的单位均为10-3 mm2/sec。FA（分数各向异性）无量纲。向量图中的比例尺代表20 mm长度，T2图像中的比例尺在其余所有图像中也代表20 mm长度。
<ol start="3">
	<li>如有必要，使用 DTI 软件程序通过加载如步骤 2.6所述获得的相位差图来校正几何畸变。将梯度回波图像中所有高于噪声水平像素的相位差图范围调整为 -π 至 +π，以表示以弧度为单位的相位偏移 Δɸ (x,y,z)（图 12）。
	<ol>
		<li>使用 DTI 软件程序根据以下公式计算像素的空间位移： 
		<img alt="磁共振成像公式 Δx=Δϕ(x,y,z)/(2π·ΔTE·BWp)，用于空间编码的关键公式。" fo:content-width="1.5in" src="/files/ftp_upload/52048/52048eq3.jpg" width="150" data-altupdatedat="1787630408000" data-alt="公式 3"> 
		其中 ΔTE 为两次回波之间的时间差，等于 2.46 × 10-3 毫秒，BWp 为实验中的相位带宽，等于 13.2 Hz/像素。根据场图对像素进行位移，并重新计算扩散张量参数（图 12）。</li>
	</ol>
	</li>
</ol>
<img alt="使用相位差图进行MRI分析、λ1校正、扩散成像、诊断放射学示意图。" src="/files/ftp_upload/52048/52048fig12highres.jpg" data-altupdatedat="1787630408000" data-alt="图12"> 
图12：使用相位差图进行几何畸变校正的示例。 图像来自一名47岁的侵袭性小叶癌志愿者。左侧第一列显示了左乳癌变层面的T2加权图像和相位差图。第二列和第3rd列的第一行显示了校正前的λ1图，第二行显示了校正后的λ1图。λ1图叠加在相应的T2加权图像上。在第2nd列中，彩色编码标尺使用水的扩散系数全范围（0.8–3.0）×10-3 mm2/sec，而在第3rd列中使用了缩减范围，阈值为1.7×10-3 mm2/sec，以突出左乳癌变组织与正常乳腺组织之间的对比度。T2图像中的比例尺在所有图像中均代表20 mm的长度。
<ol start="4">
	<li>使用能够提供增强动力学参数逐像素图谱的软件（如 3TP 图像处理软件17）对 DCE-MRI 数据集进行图像处理。</li>
</ol>

<ol><li>Turnbull, L. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Dynamic+contrast-enhanced+MRI+in+the+diagnosis+and+management+of+breast+cancer%5BTitle%5D)+AND+%22NMR+Biomed%22%5BJournal%5D)">Dynamic contrast-enhanced MRI in the diagnosis and management of breast cancer.</a> NMR Biomed. 22 (1), 28-39 (2009).</li>
<li>Sardanelli, F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Magnetic+resonance+imaging+of+the+breast%3A+recommendations+from+the+EUSOMA+working+group%5BTitle%5D)+AND+%22Eur+J+Cancer%22%5BJournal%5D)">Magnetic resonance imaging of the breast: recommendations from the EUSOMA working group.</a> Eur J Cancer. 46 (8), 1296-1316 (2010).</li>
<li>Lehman, C. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Accuracy+and+interpretation+time+of+computer-aided+detection+among+novice+and+experienced+breast+MRI+readers%5BTitle%5D)+AND+%22AJR+Am+J+Roentgenol%22%5BJournal%5D)">Accuracy and interpretation time of computer-aided detection among novice and experienced breast MRI readers.</a> AJR Am J Roentgenol. 200 (6), 683-689 (2013).</li>
<li>Riedl, C. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Magnetic+resonance+imaging+of+the+breast+improves+detection+of+invasive+cancer%2C+preinvasive+cancer%2C+and+premalignant+lesions+during+surveillance+of+women+at+high+risk+for+breast+cancer%5BTitle%5D)+AND+%22Clin+Cancer+Res%22%5BJournal%5D)">Magnetic resonance imaging of the breast improves detection of invasive cancer, preinvasive cancer, and premalignant lesions during surveillance of women at high risk for breast cancer.</a> Clin Cancer Res. 13 (20), 6144-6152 (2007).</li>
<li>Heywang-Köbrunner, S. H., Hacker, A., Sedlacek, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Magnetic+resonance+imaging%3A+the+evolution+of+breast+imaging%5BTitle%5D)+AND+%22Breast%22%5BJournal%5D)">Magnetic resonance imaging: the evolution of breast imaging.</a> Breast. 22 (2), 77-82 (2013).</li>
<li>Peters, N. H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Meta-analysis+of+MR+imaging+in+the+diagnosis+of+breast+lesions%5BTitle%5D)+AND+%22Radiology%22%5BJournal%5D)">Meta-analysis of MR imaging in the diagnosis of breast lesions.</a> Radiology. 246 (1), 116-124 (2008).</li>
<li>Warner, E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Systematic+review%3A+using+magnetic+resonance+imaging+to+screen+women+at+high+risk+for+breast+cancer%5BTitle%5D)+AND+%22Ann+Intern+Med%22%5BJournal%5D)">Systematic review: using magnetic resonance imaging to screen women at high risk for breast cancer.</a> Ann Intern Med. 148 (9), 671-679 (2008).</li>
<li>Benndorf, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Breast+MRI+as+an+adjunct+to+mammography%3A+Does+it+really+suffer+from+low+specificity%3F+A+retrospective+analysis+stratified+by+mammographic+BI-RADS+classes%5BTitle%5D)+AND+%22Acta+Radiol%22%5BJournal%5D)">Breast MRI as an adjunct to mammography: Does it really suffer from low specificity? A retrospective analysis stratified by mammographic BI-RADS classes.</a> Acta Radiol. 51 (7), 715-721 (2010).</li>
<li>Thomassin-Naggara, I., De Bazelaire, C., Chopier, J., Bazot, M., Marsault, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Trop+I+Diffusion-weighted+MR+imaging+of+the+breast%3A+advantages+and+pitfalls%5BTitle%5D)+AND+%22Eur+J+Radiol%22%5BJournal%5D)">Trop I Diffusion-weighted MR imaging of the breast: advantages and pitfalls.</a> Eur J Radiol. 82 (3), 435-443 (2013).</li>
<li>Bogner, W., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Diffusion-weighted+MR+for+differentiation+of+breast+lesions+at+3.0+T%3A+how+does+selection+of+diffusion+protocols+affect+diagnosis%5BTitle%5D)+AND+%22Radiology%22%5BJournal%5D)">Diffusion-weighted MR for differentiation of breast lesions at 3.0 T: how does selection of diffusion protocols affect diagnosis.</a> Radiology. 253 (2), 341-351 (2009).</li>
<li>Diffusion tensor imaging of the breast: preliminary clinical findings [abstr]. Proceedings of the Fourteenth Meeting of the International Society for Magnetic Resonance in Medicine. Partridge, S. C., et al. 2006 May 6-12, Seattle, Washington, , International Society for Magnetic Resonance in Medicine. Berkeley, CA. (2006).</li>
<li>Partridge, S. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Diffusion+tensor+MRI%3A+preliminary+anisotropy+measures+and+mapping+of+breast+tumors%5BTitle%5D)+AND+%22J+Magn+Reson+Imaging%22%5BJournal%5D)">Diffusion tensor MRI: preliminary anisotropy measures and mapping of breast tumors.</a> J Magn Reson Imaging. 31 (2), 339-347 (2010).</li>
<li>Partridge, S. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Diffusion+tensor+magnetic+resonance+imaging+of+the+normal+breast%5BTitle%5D)+AND+%22Magn+Reson+Imaging%22%5BJournal%5D)">Diffusion tensor magnetic resonance imaging of the normal breast.</a> Magn Reson Imaging. 28 (3), 320-328 (2010).</li>
<li>Baltzer, P. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Diffusion+tensor+magnetic+resonance+imaging+of+the+breast%3A+a+pilot+study%5BTitle%5D)+AND+%22Eur+Radiol%22%5BJournal%5D)">Diffusion tensor magnetic resonance imaging of the breast: a pilot study.</a> Eur Radiol. 21 (1), 1-10 (2011).</li>
<li>Wang, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Optimization+of+the+parameters+for+diffusion+tensor+magnetic+resonance+imaging+data+acquisition+for+breast+fiber+tractography+at+1.5%5BTitle%5D)+AND+%22T.+Clin+Breast+Cancer%22%5BJournal%5D)">Optimization of the parameters for diffusion tensor magnetic resonance imaging data acquisition for breast fiber tractography at 1.5.</a> T. Clin Breast Cancer. 14 (1), 61-67 (2014).</li>
<li>Method and apparatus for ductal tube tracking imaging for breast cancer and diagnosis and product. US Patent. , US8526698 B2 (2008).</li>
<li>Novel MRI method for breast cancer detection based on diffusion tensor tracking of the ductal trees [abstr]. Eyal, E., et al. Eighteenth Meeting of the International Society for Magnetic Resonance in Medicine, 2010 May 1-7, Stockholm, Sweden, , International Society for Magnetic Resonance in Medicine. Berkeley, CA. 362(2010).</li>
<li>Breast cancer detection and diagnosis based on diffusion tensor imaging [abstr]. Furman-Haran, E., et al. Nineteenth Meeting of the International Society for Magnetic Resonance in Medicine, 2011 May 9-13, Montreal, Quebec, , International Society for Magnetic Resonance in Medicine. Berkeley, CA. 515(2011).</li>
<li>Eyal, E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Parametric+diffusion+tensor+imaging+of+the+breast%5BTitle%5D)+AND+%22Invest+Radiol%22%5BJournal%5D)">Parametric diffusion tensor imaging of the breast.</a> Invest Radiol. 47 (5), 284-291 (2012).</li>
<li>Tagliafico, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Diffusion+tensor+magnetic+resonance+imaging+of+the+normal+breast%3A+reproducibility+of+DTI-derived+fractional+anisotropy+and+apparent+diffusion+coefficient+at+3.0+T%5BTitle%5D)+AND+%22Radiol+Med%22%5BJournal%5D)">Diffusion tensor magnetic resonance imaging of the normal breast: reproducibility of DTI-derived fractional anisotropy and apparent diffusion coefficient at 3.0 T.</a> Radiol Med. 117 (6), 992-1003 (2012).</li>
<li>Cakir, O., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Comparison+of+the+diagnostic+performances+of+diffusion+parameters+in+diffusion+weighted+imaging+and+diffusion+tensor+imaging+of+breast+lesions.Eur%5BTitle%5D)+AND+%22J+Radiol%22%5BJournal%5D)">Comparison of the diagnostic performances of diffusion parameters in diffusion weighted imaging and diffusion tensor imaging of breast lesions.Eur.</a> J Radiol. 82 (12), e801-e806 (2013).</li>
<li>Tsougos, I., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+contribution+of+diffusion+tensor+imaging+and+magnetic+resonance+spectroscopy+for+the+differentiation+of+breast+lesions+at+3T%5BTitle%5D)+AND+%22Acta+Radiol%22%5BJournal%5D)">The contribution of diffusion tensor imaging and magnetic resonance spectroscopy for the differentiation of breast lesions at 3T.</a> Acta Radiol. 55 (1), 14-23 (2014).</li>
<li>Wiederer, P. azahrS., Leo, C., Nanz, D., Boss, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Quantitative+breast+MRI%3A+2D+histogram+analysis+of+diffusion+tensor+parameters+in+normal+tissue%5BTitle%5D)+AND+%22Magn+Reson+Mater+Phy%22%5BJournal%5D)">Quantitative breast MRI: 2D histogram analysis of diffusion tensor parameters in normal tissue.</a> Magn Reson Mater Phy. 27, 185-193 (2014).</li>
<li>Nissan, N., Furman-Haran, E., Shapiro-Feinberg, M., Grobgeld, D., Degani, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Diffusion-tensor+MR+imaging+of+the+breast%3A+hormonal+regulation%5BTitle%5D)+AND+%22Radiology%22%5BJournal%5D)">Diffusion-tensor MR imaging of the breast: hormonal regulation.</a> Radiology. 271 (3), 672-680 (2014).</li>
<li>Cooper, A. P. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aCooper%2C+AP+%22On+the+Anatomy+of+the+breast%22">On the Anatomy of the breast</a>. , Longman, Orme, Green, Brown, and Longmans. London. (1840).</li>
<li>Ohtake, T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Computer-assisted+complete+three-dimensional+reconstruction+of+the+mammary+ductal%2Flobular+systems%3A+implications+of+ductal+anastomoses+for+breast-conserving+surgery%29%5BTitle%5D)+AND+%22Cancer%22%5BJournal%5D)">Computer-assisted complete three-dimensional reconstruction of the mammary ductal/lobular systems: implications of ductal anastomoses for breast-conserving surgery).</a> Cancer. 91 (12), 2263-2272 (2001).</li>
<li>Going, J. J., Moffat, D. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Escaping+from+Flatland%3A+clinical+and+biological+aspects+of+human+mammary+duct+anatomy+in+three+dimensions%5BTitle%5D)+AND+%22J+Pathol%22%5BJournal%5D)">Escaping from Flatland: clinical and biological aspects of human mammary duct anatomy in three dimensions.</a> J Pathol. 203 (1), 538-544 (2004).</li>
<li>Hancu, I., Govenkar, A., Lenkinski, R. E., Lee, S. K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((On+shimming+approaches+in+3T+breast+MRI%5BTitle%5D)+AND+%22Magn+Reson+Med%22%5BJournal%5D)">On shimming approaches in 3T breast MRI.</a> Magn Reson Med. 69 (3), 862-867 (2013).</li>
<li>Basser, P. J., Jones, D. K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Diffusion-tensor+MRI%3A+theory%2C+experimental+design+and+data+analysis+-+a+technical+review%5BTitle%5D)+AND+%22NMR+Biomed%22%5BJournal%5D)">Diffusion-tensor MRI: theory, experimental design and data analysis - a technical review.</a> NMR Biomed. 15 (7-8), 465-467 (2002).</li>
<li>Reese, T. G., Heid, O., Weisskoff, R. M., Wedeen, V. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Reduction+of+eddy-current-induced+distortion+in+diffusion+MRI+using+a+twice-refocused+spin+echo%5BTitle%5D)+AND+%22Magn+Reson+Med%22%5BJournal%5D)">Reduction of eddy-current-induced distortion in diffusion MRI using a twice-refocused spin echo.</a> Magn Reson Med. 49 (1), 1771-1782 (2003).</li>
<li>Jezzard, P., Balaban, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Correction+for+geometric+distortion+in+echo+planar+images+from+B0+field+variations%5BTitle%5D)+AND+%22Magn+Reson+Med%22%5BJournal%5D)">Correction for geometric distortion in echo planar images from B0 field variations.</a> Magn Reson Med. 34 (1), 65-73 (1995).</li>
<li>Kelcz, F., Furman-Haran, E., Grobgeld, D., Degani, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Clinical+testing+of+high-spatial+resolution+parametric+contrast-enhanced+MR+imaging+of+the+breast%5BTitle%5D)+AND+%22AJR+Am+J+Roentgenol%22%5BJournal%5D)">Clinical testing of high-spatial resolution parametric contrast-enhanced MR imaging of the breast.</a> AJR Am J Roentgenol. 179 (6), 1485-1492 (2002).</li>
<li>Le Bihan, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Diffusion+tensor+imaging%3A+concepts+and+applications%5BTitle%5D)+AND+%22J+Magn+Reson+Imaging%22%5BJournal%5D)">Diffusion tensor imaging: concepts and applications.</a> J Magn Reson Imaging. 13 (4), 534-546 (2001).</li>
<li>Marquardt, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((An+Algorithm+for+Least-Squares+Estimation+of+Nonlinear+Parameter%5BTitle%5D)+AND+%22SIAM+Journal+on+Applied+Mathematics%22%5BJournal%5D)">An Algorithm for Least-Squares Estimation of Nonlinear Parameter.</a> SIAM Journal on Applied Mathematics. 11 (2), 431-441 (1963).</li>
<li>Jolliffe, I. T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Principal+Component+Analysis%3A+A+Beginner%27s+Guide+%E2%80%94+I.+Introduction+and+application%5BTitle%5D)+AND+%22Weather%22%5BJournal%5D)">Principal Component Analysis: A Beginner's Guide — I. Introduction and application.</a> Weather. 45 (10), 375-382 (1990).</li>
<li>Mayr, N. A., Staples, J. J., Robinson, R. A., Vanmetre, J. E., Hussey, D. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Morphometric+studies+in+intraductal+breast+carcinoma+using+computerized+image+analysis%5BTitle%5D)+AND+%22Cancer%22%5BJournal%5D)">Morphometric studies in intraductal breast carcinoma using computerized image analysis.</a> Cancer. 67 (11), 2805-2812 (1991).</li>
<li>Ramsay, D. T., Kent, J. C., Hartmann, R. A., Hartmann, P. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Anatomy+of+the+lactating+human+breast+redefined+with+ultrasound+imaging%5BTitle%5D)+AND+%22J+Anat%22%5BJournal%5D)">Anatomy of the lactating human breast redefined with ultrasound imaging.</a> J Anat. 206 (6), 525-534 (2005).</li>
<li>Gooding, M. J., Mellor, M., Shipley, J. A., Broadbent, K. A., Goddard, D. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Automatic+mammary+duct+detection+in+3D+ultrasound%5BTitle%5D)+AND+%22Med+Image+Comput+Comput+Assist+Interv%22%5BJournal%5D)">Automatic mammary duct detection in 3D ultrasound.</a> Med Image Comput Comput Assist Interv. 8 (1), 434-441 (2005).</li>
<li>Diffusion tensor based reconstruction of the ductal tree [abstr]. Reisert, M. M., Eyal, E., Grobgeld, D., Degani, H., Hennig, J. Nineteenth Meeting of the International Society for Magnetic Resonance in Medicine, 2011 May 9-13, Montreal, Quebec, , International Society for Magnetic Resonance in Medicine. Berkeley, CA. 3649(2011).</li>
<li>Jezzard, P., Clare, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Sources+of+distortion+in+functional+MRI+data%5BTitle%5D)+AND+%22Hum+Brain+Mapp%22%5BJournal%5D)">Sources of distortion in functional MRI data.</a> Hum Brain Mapp. 8 (2-3), 80-85 (1999).</li>
<li>Jones, D. K., Cercignani, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Twenty-five+pitfalls+in+the+analysis+of+diffusion+MRI+data%5BTitle%5D)+AND+%22NMR+Biomed%22%5BJournal%5D)">Twenty-five pitfalls in the analysis of diffusion MRI data.</a> NMR Biomed. 23 (7), 803-820 (2010).</li>
</ol>

作者无任何利益冲突需要披露。

本研究展示了在3 T高空间分辨率（约8 mm3）下进行的弥散张量成像（DTI）能够对双侧乳腺全部纤维腺组织中的水弥散张量参数进行活体测量。in vivo 在本研究过程中开发的一套算法与软件，实现了对大规模DTI数据集的分析，并可同时显示所有乳腺切片的各向弥散系数λ1、λ2、λ3、ADC以及各向异性指数λ1–λ3和FA的弥散参数图。本研究还强调了各类弥散系数及各向异性指数在乳腺癌检测与诊断中的潜在价值。研究结果表明，主弥散系数λ1是检测乳腺癌的主要参数，能够有效区分恶性乳腺组织与正常乳腺组织。另一个具有高灵敏度但特异性明显低于λ1的独立辅助参数是最大各向异性指数，可用于对主弥散系数检测结果的确认。
结果还展示了磁共振成像（MRI）在测量方面的应用能力 体内 水的定向扩散，并量化双侧乳房整个导管/腺体系统中扩散的各向异性。为确保各向异性测量具有均匀的精度，采用了30个扩散梯度方向。选择120毫秒的长回波时间，并优化空间分辨率，从而揭示了导管内水扩散受限的情况。正常导管的平均尺寸据报道为90 μm，约70%位于40–100 μm范围内。36根据爱因斯坦方程，自由水扩散的平均位移 x = (6Dt)1/2，其中 D 是自由水扩散系数，t 为扩散时间。在我们的实验中，x 约为 25 μm，因此只有一小部分导管内水分子受到乳腺导管壁的限制，导致乳腺组织的各向异性分数（FA）值相对较低（上限值约 0.3），相比之下，脑白质中的 FA 值较高（≥0.5）。在阿斯特利·库珀（Sir Astley Cooper）的研究中25 研究发现，人类乳腺组织由多个独立的叶组成，每叶均包含一个导管树。对哺乳期乳腺的初步超声研究37 尝试了使用二阶形状测量来检测部分导管结构38然而，迄今为止尚无任何成像方法能够成功显示双侧乳房的完整导管系统 体内从DTI结果中获取完整3D导管树系统的算法仍在开发中，但矢量图已清晰展现出揭示乳腺树精细且高度多样解剖结构的潜力。近期已有基于我们DTI数据集对乳腺进行完整3D追踪的初步尝试报道。39.
在本研究过程中，注意到存在一些技术局限性，主要源于脂肪抑制效率低下以及EPI失真。脂肪抑制通过脂肪饱和序列（FAT-SAT）实现，该方法在致密型乳腺中通常效果良好。然而，在脂肪含量较高的乳腺中，采用了对B1场不均匀性更不敏感且效率更高的谱选择衰减反转恢复技术（SPAIR）。基于EPI的扩散成像方案还受到梯度诱导的涡流、B0磁场不均匀性以及磁化率差异的额外限制40,41。这些局限性通过在水和脂肪信号上进行迭代式的自动与手动匀场，并选择尽可能低的回波间隔而得以最小化。此外，如方案中所述，偶尔还会应用几何失真的后处理校正。总体而言，对于水信号较强的致密型乳腺，能够克服上述技术局限；然而，本研究中5例脂肪含量极高的乳腺病例因上述限制而无法进行分析。
综上所述，已开发出一套用于乳腺扩散张量 MRI 的实验方案及图像处理工具。这种完全安全、快速且无创的方法能够精确解析乳腺结构，有助于在临床中检测乳腺癌。研究发现，主扩散系数 λ1 和最大各向异性指数 λ1-λ3 可作为两个独立的扩散参数用于癌症检测。对乳腺癌患者开展的临床研究表明，利用这两个参数的检测效能与动态对比增强 MRI（DCE MRI）相当。因此，该方法基于对乳腺结构和癌症发生机制的基础理解，依赖于揭示可量化的显著变化，且具备安全、快速的特点，使其成为研究导管增殖等发育过程相关问题的有力工具，并适用于在大规模前瞻性试验中评估其在乳腺癌筛查中的临床价值。

乳腺癌是全球女性癌症中最常见的类型。早期发现乳腺癌对于提高患者的生活质量和生存率具有关键作用。目前乳腺癌检测的标准方法主要基于X射线乳腺摄影和超声扫描。这些技术在检测致密乳腺组织中的病灶时敏感性和特异性不足，从而推动了其他技术的发展，包括乳腺磁共振成像（MRI）。动态对比增强（DCE）MRI已被确立为检测和诊断乳腺癌的有力工具1,2 ，并常借助计算机辅助诊断手段3。目前该技术主要用于特定情况，例如高风险患者4，而非常规筛查，可能原因在于其成本较高、需要注射对比剂、缺乏标准化，以及在区分良恶性病灶时特异性存在差异，其数值范围从较低/中等水平5,6 到结合乳腺X线摄影与DCE-MRI所获得的高水平不等7,8 。近年来，扩散加权MRI及其生成的表观扩散系数（ADC）图被评估为DCE-MRI的补充方法，研究表明ADC值有助于区分癌症、良性病灶和正常乳腺组织9,10。此外，针对健康志愿者和乳腺病灶患者开展的乳腺扩散张量成像（DTI）研究已在1.5 T11-15和3 T16-24磁场强度下启动。大多数此类研究报道了ADC和分数各向异性（FA）值11,12,14,15,20-23，并发现这两个参数具有可重复性，其中ADC值的可重复性高于FA13,20。这些研究结果表明，与正常组织和良性病灶相比，恶性病灶表现出较低的ADC值，然而关于FA值及其诊断能力的研究结果存在矛盾11,12,14,20-23。在一系列3 T DTI研究中，还报告了乳腺组织坐标系中三个张量特征值和特征向量的数值，结果以主特征向量的矢量图以及特征值、ADC、FA和最大各向异性指数的参数图形式呈现16-19,24。在这些研究中，主扩散特征值和最大各向异性被证明是检测和诊断癌性病灶最敏感的独立参数
乳房由纤维腺体组织和脂肪组织构成。纤维腺体组织进一步分为多个小叶，其大小和形状具有高度变异性。每个小叶的显微结构包括发挥功能的乳腺导管树及形成腺体组织的相应腺泡，以及周围的结缔纤维组织。大多数乳腺恶性肿瘤起源于导管或腺泡上皮细胞的异常增殖，进而发展 原位 癌变组织通过浸润周围组织转化为浸润性癌。因此，乳腺导管/小叶结构是研究乳腺恶性转化的关键区域。
1840年，阿斯特利·库珀爵士首次通过向乳房切除标本的导管内注射彩色石蜡，对导管树的结构特征进行了离体（ex vivo）研究25。近年来，已有研究利用乳房切除标本，通过对少数人类全乳腺导管树进行计算机辅助追踪实现了相关成像26,27。本文所展示的工作表明，通过活体（in vivo）扩散张量成像获得的参数可提供与乳腺组织独特微结构特征相关的信息，从而实现无创性乳腺癌检测。
乳腺扩散张量成像的物理原理基于磁共振成像（MRI）测量和量化受限环境中水分子各向异性扩散的能力28。通常情况下，均质溶液中的水分子扩散是自由且各向同性的；然而，当水分子运动因受到不可渗透壁的限制而受阻时，扩散则变为各向异性：沿壁方向为快速的自由扩散，而垂直于壁的方向则为较慢的受限扩散（图1）。组织中的水分子扩散过程较为复杂，依赖于细胞内外腔室的结构和生理特征，包括细胞大小、细胞密度、细胞外间隙的迂曲度以及跨膜水交换，同时也受血管和淋巴网络存在的影响（图2）。
<img alt="扩散过程示意图：各向同性自由扩散与各向异性受限扩散；运动方式比较。" src="/files/ftp_upload/52048/52048fig1highres.jpg" data-altupdatedat="1787630408000" data-alt="图1"> 
图1：自由扩散与受限扩散。 水分子自由扩散（左）和受不可穿透壁限制的扩散（右）的示意图。
<img alt="显示分子运动受限和受阻的细胞内和细胞外扩散示意图。" src="/files/ftp_upload/52048/52048fig2highres.jpg" data-altupdatedat="1787630408000" data-alt="Figure 2"> 
图2：组织中的复杂扩散。 细胞系统中水扩散的示意图，显示水分子在细胞外和细胞内区室中的运动，以及这两个区室之间的水交换（箭头）。
由于乳腺具有特殊的结构特征，水分子在乳腺导管和小叶中的扩散表现为受限且各向异性的运动：沿导管和小叶管壁方向的扩散接近自由扩散，而在垂直于管壁的方向上，扩散则受到由双层细胞和基底膜构成的管壁限制。因此，导管/腺体系统中的扩散相对较快且具有各向异性。另一方面，导管周围结缔纤维组织中的扩散则快速且各向同性，这是由于该组织含水量高、细胞密度低所致（图3和图4）。当发生恶性病变时，癌细胞阻塞导管和小叶，增加了水分子运动的迂曲度和限制性，导致各个方向的扩散系数以及各向异性运动均下降（图3和图4）。
<img alt="癌细胞扩散对比；各向异性与各向同性水扩散；示意图与图像" src="/files/ftp_upload/52048/52048fig3highres.jpg" data-altupdatedat="1787630408000" data-alt="图 3"> 
图 3：乳腺小叶中的扩散。 穿过小叶的剖面示意图及其中单个小叶内水分子的扩散情况。左侧：水分子扩散受到小叶壁的限制，表现为沿小叶壁方向扩散较快，而垂直于小叶壁方向的扩散受限。右侧：含癌细胞的小叶中的扩散情况。细胞外间隙中的扩散受到显著阻碍，但在各个方向上相似，因此几乎呈各向同性。
<img alt="分形生长模式示意图，包含分支模拟，展示微生物扩散动态。" src="/files/ftp_upload/52048/52048fig4highres.jpg" data-altupdatedat="1787630408000" data-alt="图4"> 
图4：导管树系统中的水扩散。 左图：注入彩色石蜡的乳腺导管，显示其辐射方向及分支间连接25。中图：正常导管树的示意图，箭头表示导管内（黑色箭头）和结缔组织中（绿色箭头）的扩散方向。右图：含有两个癌细胞病灶（紫色）的导管树示意图，红色箭头表示癌组织中的扩散。 
本文详细描述了扩散张量扫描方法，以及用于检测乳腺恶性肿瘤的DTI数据集的处理算法和软件分析。所有癌症均通过乳腺活检和/或手术标本的组织病理学结果确诊。我们还描述了获取乳腺解剖特征的T2加权扫描方案，以及作为评估DTI检测敏感性参考方法的DCE扫描方案。<a href="https://www.jove.com/files/ftp_upload/52048/52048fig4large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>

<table><tbody><tr><td>磁性钆喷酸葡胺</td><td>Soreq, Yavne, 以色列</td><td></td><td>0.5 M</td></tr><tr><td>3 特斯拉磁共振成像仪，MAGNETON Trio</td><td>Siemens, Erlangen, 德国</td><td></td><td>120</td></tr><tr><td>双侧乳腺阵列线圈</td><td>Siemens, Erlangen, 德国</td><td></td><td>4 通道</td></tr><tr><td>双侧乳腺阵列线圈</td><td>In-Vivo, Orlando, 佛罗里达州</td><td></td><td>7 通道</td></tr><tr><td>自动注射泵，Spectris Solaris MR 注射器</td><td>Medrad, Indianola, 宾夕法尼亚州</td><td></td><td></td></tr><tr><td>DTI 图像处理软件</td><td>自主开发</td><td></td><td>产权属于 &lt;a href=&quot;http://www.yedarnd.com/ContactUs.asp&quot; target=&quot;_blank">Yeda Research and Development Co. Ltd&lt;/a&gt;</td></tr><tr><td>3TP 图像处理软件</td><td>自主开发</td><td></td><td>产权属于 &lt;a href=&quot;http://www.yedarnd.com/ContactUs.asp&quot; target=&quot;_blank">Yeda Research and Development Co. Ltd&lt;/a&gt;</td></tr></tbody></table>

追踪乳腺结构特征并利用磁共振扩散张量成像检测乳腺癌

Tracking the Mammary Architectural Features and Detecting Breast Cancer with Magnetic Resonance Diffusion Tensor Imaging

乳腺癌是全球女性中最常见的癌症。早期发现乳腺癌对于提高乳腺癌患者的生活质量和生存率具有关键作用。本文介绍了一种基于扩散张量成像（DTI）追踪乳腺结构元素的新型乳腺癌检测方法。
本文重点介绍为适应乳腺纤维腺组织的扩散特性及其在恶性转化过程中的变化而设计的扫描方案、图像处理算法及相关软件。最终输出结果为逐像素的矢量图，可追踪整个乳腺导管腺体树的结构，同时提供扩散张量系数和各向异性指数的参数图。
通过检测包括68名经组织病理学检查确诊为乳腺癌患者在内的女性志愿者，对该方法检测乳腺癌的效率进行了测试。与正常乳腺组织相比，含有癌细胞的区域表现出扩散系数和最大各向异性指数显著降低，从而为勾画恶性肿瘤边界的内在对比提供了依据。总体而言，DTI参数检测乳腺癌的敏感性较高，特别是在致密型乳腺中，其性能可与目前需要注射对比剂的标准乳腺MRI方法相媲美。因此，该方法为乳腺癌检测提供了一种完全无创、安全且敏感的工具。

该方法最初通过扫描处于不同激素水平的健康志愿者进行测试和验证。图11展示了使用我们软件获得的一名年轻健康志愿者乳腺中央切片的参数图和矢量图，该志愿者的纤维腺体组织比例相对较高，这在T2加权图像上可以清晰看出（所有灰色区域为纤维腺体组织，亮区为脂肪）。主扩散系数λ1的方向在矢量图v1中显示，大部分像素指向乳头方向。如预期所示，扩散张量系数的数值从λ1到λ2再到λ3依次降低。利用这三个扩散系数可计算平均扩散率ADC、分数各向异性（FA）以及最大各向异性λ1 - λ3。请注意，在FA图和λ1 - λ3图中，这两个指数的空间分布具有高度一致性。
在绝经后女性中，乳腺通常密度较低，如图13中一位63岁健康志愿者的示例所示。老年志愿者的扩散系数较低，但各向异性指数较高，这可能是因为导管和小叶的直径较小，从而对垂直于导管壁方向的扩散施加了更强的限制。另一个体现张量对导管尺寸敏感性的示例如图14中一位哺乳期志愿者所示。由于乳汁是脂肪球分散于以水为基础的液体中形成的胶体，其中含有溶解的碳水化合物、蛋白质聚集体及矿物质，其黏度高于导管内正常的水性液体，因此哺乳期乳腺的水分子扩散系数低于绝经前女性24。此外，靠近乳头区域的导管较大，导致垂直于导管方向的扩散限制低于非哺乳期乳腺，因而各向异性指数也较低。然而，在哺乳期乳腺的小叶—后部区域，各向异性仍然较高。
在恶性肿瘤患者中，弥散张量参数的主要变化表现为三个弥散系数显著降低。λ 的变化1 具有最高的信噪比19本研究纳入了68例经病理确诊的患者，其中包括33例浸润性导管癌（IDC）患者和19例导管原位癌患者 原位 （导管原位癌），13例为浸润性小叶癌（ILC），3例为其他恶性肿瘤。多名患者患有多灶性或多中心性乳腺癌。肿瘤大小范围为3至95毫米，中位大小为14毫米，四分位距为10至30毫米。在所有这些患者中，DCE与DTI的检出率相当。然而，有5例在DCE中显示强化并经病理证实为癌症的病例，由于脂肪乳腺中的技术问题（主要与磁场不均匀性和脂肪抑制不充分导致的畸变和伪影有关）未被纳入DTI分析。
图15、16和17展示了用于检测乳腺恶性肿瘤的两个主要扩散参数 λ1 和 λ1 - λ3 的典型参数图。该图除了显示T2加权图像外，还展示了采用3TP彩色编码方法进行DCE分析的结果。如引言中所述，导管或小叶及其周围存在的恶性细胞会阻碍细胞外间隙的扩散，显著降低扩散系数。此外，由于癌细胞以混乱的方式向各个方向扩散，不再具有明确的方向性，导管结构所导致的各向异性也随之消失。各向异性分数（FA）不是乳腺癌检测的合适参数，因为该参数经平均扩散率归一化后，在癌组织中的数值与正常组织相似，均较高19。然而，最大各向异性确实为癌症检测提供了依据（图15-17）。尽管如此，由于结缔纤维组织也接近各向同性，并表现出较低的最大各向异性值，因此该参数的特异性低于扩散系数 λ1，仅作为 λ1 的辅助参数使用。
图17 还展示了DTI在表征新辅助化疗反应方面的能力。在此例中，患者对治疗完全响应（4个周期的阿霉素+环磷酰胺，随后4个周期的紫杉醇）。事实上，治疗反应导致扩散系数显著升高至接近正常乳腺组织的水平，提示存在修复性结缔组织取代了癌细胞。在其他4例对治疗产生响应的患者中也获得了类似结果。
<img alt="MRI脑部分析；包括T2加权图像、矢量图和扩散指标（ADC、FA）" src="/files/ftp_upload/52048/52048fig13highres.jpg" data-altupdatedat="1787630408000" data-alt="图13"> 
图13：一名63岁健康志愿者乳腺中央切片的参数化DTI图叠加于同一层面T2加权图像上的输出结果。 Dir. — 方向图，红色表示从左到右，绿色表示从前到后，蓝色表示从头到脚方向。注意矢量以白色线条标示。所有扩散系数和最大各向异性指数的单位均为1×10-3 mm2/sec。FA无单位。矢量图中的比例尺代表20 mm长度，T2图像中的比例尺在所有其余图像中也代表20 mm长度。
<img alt="肩部MRI矢量图与热图；弥散张量成像；组织微结构的测量" src="/files/ftp_upload/52048/52048fig14highres.jpg" data-altupdatedat="1787630408000" data-alt="图14"> 
图14：叠加于同一层面T2加权图像上的40岁哺乳期志愿者乳腺中央切片的参数化DTI图输出结果。 Dir. — 方向图，其中红色表示从左到右，绿色表示从前到后，蓝色表示从头到脚的方向。所有弥散系数和最大各向异性指数的单位均为1×10-3 mm2/sec。FA无单位。请注意，在方向图和矢量图中，大多数像素均朝向乳头方向排列。同时请注意，与健康志愿者的图7和图10中的数值相比，弥散系数较低。矢量图中的比例尺代表20 mm的长度，T2图像中的比例尺在其余所有图像中也代表20 mm的长度。
<img alt="Breast tumor MRI analysis; diffusion map comparison; λ1-λ3 eigenvalues; imaging results." src="/files/ftp_upload/52048/52048fig15highres.jpg" data-altupdatedat="1787630408000" data-alt="Figure 15"> 
图15：λ的参数图1 和 λ1 - λ3 一名38岁患者患有弥漫性浸润性导管癌。 在第1个st DTI 参数以 1.7 x 10 为阈值呈现-3 毫米2/sec 对应 λ1 和 0.6 mm2/秒，对应λ1-λ3 （所有高于阈值的数值均以紫色标出）。在2nd 原始图像中，高于阈值的数值未着色，显示其下方的T2加权图像。该图还包含在1st 左侧一列为同一层面的T2加权图像及通过3TP方法获得的DCE参数图，右侧为DTI参数图。T2图像中的标尺在所有图像中均代表20 mm长度。请注意，DCE图像的面内空间分辨率约为DTI的两倍，但在3TP与λ的病灶位置和大小上具有高度的视觉一致性。1 参数图
<img alt="Breast MRI diagrams comparing T2, 3TP, λ1, and λ1-λ3 for perfusion and diffusion analysis." src="/files/ftp_upload/52048/52048fig16highres.jpg" data-altupdatedat="1787630408000" data-alt="Figure 16"> 
图16：λ的参数图1 和 λ1-λ3 一名60岁患者的低级别导管原位癌。 该图还包含了1中的内容st 左侧一列为同一层面的T2加权图像和通过3TP方法获得的DCE参数图，右侧为DTI参数图。T2图像中的标尺在所有图像中均代表20 mm长度。请注意DTI在富含脂肪的乳腺中检测癌症的能力。
<img alt="Breast MRI analysis of diffusion tensor imaging: T2, MIP, 3TP sequences, λ1, λ1-λ3 maps." src="/files/ftp_upload/52048/52048fig17highres.jpg" data-altupdatedat="1787630408000" data-alt="Figure 17"> 
图17：λ的参数图1 和 λ1 - λ3 一名39岁侵袭性小叶癌患者在新辅助化疗前后的状况。 该患者接受了两次扫描，分别在治疗前和手术前进行，治疗方案为4个周期的阿霉素（Adryamicin）联合环磷酰胺（Cycloxan）以及4个周期的紫杉醇（Taxol）。注意观察λ的增加1 和在 λ1-λ3 在对治疗有反应的癌症区域。1st 原始图像显示治疗前获得的图像以及第2天nd raw 显示了与图 1 中大致相同区域的图像st 行，术前获得。治疗前，癌变区域表现出 λ1 和 λ1-λ3 低于其阈值，治疗后这两项参数升高至超过其阈值。该图还包含在1中的st 左侧列为 T2 加权图像；在 2nd 柱状多图像投影（MIP），通过将对比前图像从对比后2分钟的图像中减去而获得；在3rd 通过3TP方法获得的DCE参数图。T2图像中的比例尺在所有图像中均表示20 mm的长度

我们介绍了如何利用磁共振成像获取乳腺扩散张量的参数图和矢量图。成像后处理的方案及最终输出结果经过优化，可用于追踪乳腺结构特征并检测乳腺恶性病变。

乳腺癌是全球女性中最常见的癌症病因。早期发现乳腺癌对于提高治疗质量和改善预后具有关键作用。

以下实验的总体目标是通过一种完全无创且安全的成像方法，追踪乳腺的微观结构，以检测乳腺癌。该方法通过对双侧乳腺进行磁共振成像，并在整个乳腺纤维腺组织中测量水分子扩散张量，从而揭示正常及病变乳腺的微观结构和细胞特征。第二步，采用对称张量模型处理扩散张量成像数据集，以逐像素计算方向性扩散系数和各向异性指数。随后生成方向性扩散系数的矢量图和彩色编码参数图，同时生成平均表观扩散系数、分数各向异性指数和最大各向异性指数，用于确定检测恶性病变和评估病理特征的最佳参数。结果表明，扩散张量参数能够揭示乳腺的结构特征，并通过发现扩散系数值和最大各向异性值较正常组织降低的图像区域来检测恶性病变。这种新的扩散张量成像方法在乳腺癌诊断准确性方面与标准的动态对比增强三时间点方法相当。该技术相较于现有的磁共振成像方法（如动态对比增强磁共振）的主要优势在于，DTI 无需注射外源性对比剂，因此完全安全且无创。该方法还具有相对较快的检测速度，且其临床表现不受激素变化的影响。该方法通过检测因癌细胞增殖导致的细胞密度增加，有助于回答乳腺影像学中的关键问题，包括乳腺癌的检测、诊断与治疗、疗效监测，以及对其他具有指状腺体结构的器官（如前列腺和肾脏）的成像。通常，初学者在应用该方法时会遇到困难，因为乳腺扩散张量成像（DTI）需要通过对所有实验参数进行精细调整来优化扫描方案，同时还需要先进的原始数据集处理软件，以生成清晰的扩散张量参数的参数图。我们产生这一方法的想法，源于在一本1840年出版的乳腺解剖学书籍中看到的乳腺导管系统图示，以及近年来阿特金（Atkin）及其同事在2001年和2004年发表的关于乳腺导管系统的三维重建计算机模拟研究。该方法的可视化演示至关重要，因为该方法的应用和意义依赖于彩色编码的参数图，这些参数图可用于评估乳腺组织并检测癌症。该操作将由Edna Fran博士、Mara Shapiro Feinberg博士、Noam Nisan博士、Dove Gogel先生以及我们两位优秀的MRI技术人员Fania Tao女士和Holster先生进行演示。新的乳腺扩散张量成像方法将与当前标准的动态对比增强乳腺MRI方法进行比较评估，因此首先需在患者任一上肢置入静脉导管以注射对比剂，随后将患者摆位以进行乳腺扫描。患者必须采取俯卧位，双侧乳房自然悬垂于乳腺线圈的双侧开口中。同时，患者头部和颈部应枕于枕头上以确保舒适。现在让患者将双臂伸过头顶，并检查乳房的位置是否正确。再次确保每侧乳房均位于线圈开口的中心位置，自然下垂并尽可能深入线圈内。现在，将自动注射泵连接至静脉导管，移动患者进入扫描仪，并开始扫描。利用定位像，确定双侧乳房的位置，设定扫描视野，并确定完全覆盖双侧乳房（延伸至腋窝及胸壁）所需的层数和层厚。在所有扫描序列中保持这三个实验参数一致。接下来，使用扫描仪的接口计算机，定位包含双侧乳房和腋窝的区域，以定义匀场盒。然后采用迭代匀场策略，在观察质子谱信号的同时优化磁场。通过将频率分别对准水共振和脂肪共振进行匀场调节。为分离脂肪和水信号，需将照射频率中心对准水的共振频率，以优化信号强度和形态。为此，现应用一个二维横向T2加权快速自旋回波多层面序列，采用高空间分辨率，不使用脂肪饱和，结合GRAPPA技术，加速因子为2。根据表一设置可变参数。接下来，进行脂肪抑制自旋回波弥散张量成像（DTI）。采用两次重聚焦的回波平面成像序列，再次应用梯度反转相位编码（grpa），加速因子为2。对于自旋回波DTI，请使用以下参数，其中最关键的参数以紫色标出。这些参数包括回波时间te、空间分辨率、扩散梯度方向数以及梯度b值。接下来应用场图序列以校正平面回波成像中的几何畸变。按照Desert和Balaban所述方法获取相位差分图像，确保成像序列包含具有两个不同同相回波时间的二维横断面梯度回波图像，并确保其相位和编码方向与DTI序列一致。表一也列出了该成像序列应使用的参数值。下一步是采用动态对比增强方案，使用无脂肪抑制的3D快速梯度回波序列，并根据三点时间法优化参数。再次查阅表一以获取参数值。在第二次扫描前图像采集结束前15秒，使用计算机控制的自动注射器，以每秒2毫升的速度注射对比剂。随后以相同速度自动推注20毫升生理盐水作为冲洗。现在，在注射后七个固定时间点连续采集七组三维数据集。接下来，采用与T2加权成像类似的实验参数，创建一个二维横断位T2加权脂肪抑制快速自旋回波多层序列。整个检查时间为30.5分钟，其中DTI序列持续时间为6分钟。首先将全部数据集传输至远程工作站或个人计算机，然后使用专门用于分析乳腺扩散、张量成像和D-C-E-M-R-I的程序打开数据。此处使用一个自主开发的软件包来完成DTI图像处理的全部步骤。首先，在乳腺外部及剩余组织的三到四个感兴趣区域评估噪声水平，每个区域约100个像素。找出高于噪声水平的像素中的最大噪声值。利用30个梯度方向计算所有60个切片中对称张量D的六个扩散系数。使用Staal-Tanner方程和非线性回归拟合程序，然后通过主成分分析对每个像素中的对称扩散张量进行对角化。该过程为每个像素生成三个 iGen 载体，用于定义椭球体形状三个正交主轴（AE）上的扩散方向。这些值与组织的扩散坐标系一致，并对应决定各向异性扩散系数的扩散特征值。接下来，针对每个像素计算表观扩散系数，即三个特征值的平均值。然后计算最大绝对各向异性指数，即 Lambda 1 和 Lambda 3 的特征值之间的差值。现在计算每个像素的分数各向异性指数。数值范围从零到一，零表示各向同性扩散，一表示单方向自由扩散。完成所有这些计算后，利用软件的绘图功能生成主特征向量（即第一特征向量 V1）的矢量图。同时创建一幅用三种颜色显示 V1 主要方向的图谱，并将其叠加到同一层面的 T2 加权图像上。在 DTI 软件中，构建扩散张量成像，生成参数图，以显示每个层面中每个像素的全部扩散张量参数值。然后将这些值叠加到同一层面的T2加权图像上。该方法在一名年轻的健康志愿者身上进行了测试和验证。在T2加权图像上，纤维腺组织的相对较高比例可清晰地显示为灰色区域。明亮区域为脂肪组织。主扩散系数 Lambda one 的方向显示，大量像素指向乳头方向。与预期一致，扩散张量系数的数值从 Lambda one 到 Lambda two 再到 Lambda three 逐渐降低。利用这三个扩散系数可计算平均扩散率（A DC）、分数各向异性以及最大各向异性，其中分数各向异性与最大各向异性高度一致。对68例乳腺恶性肿瘤患者进行了DTI扫描。结果显示，癌组织中的最大扩散系数（Lambda one）和最大各向异性（Lambda one 减去 Lambda three）的参数图显著低于正常乳腺组织，且在乳腺癌检测中具有最高的效能。在这些患者中，Lambda one 和 Lambda one 减去 Lambda three 的低值区域与通过三点时间法彩色编码的动态对比增强 MRI 所识别出的癌变区域相一致。第二组 Lambda one 与 Lambda one 减去 Lambda three 图像显示了一位随后接受新辅助化疗的患者的癌变情况。在手术前成功完成新辅助化疗后，弥散张量成像（DTI）能够表征组织的治疗反应。扩散系数显著升高，反映出修复性结缔组织替代了癌细胞。一旦掌握序列采集与图像处理方法，扩散张量成像数据集可在数分钟内完成处理，生成覆盖整个乳腺组织的扩散参数图，从而实现对恶性病变的高灵敏度检测。在实施该操作过程中，需注意优化磁共振成像（MRI）扫描仪性能、磁场均匀性，并严格遵循扩散张量成像（DTI）的扫描协议及图像处理规范。完成该操作后，可进一步采用其他MRI方法（如动态对比增强成像）来验证DTI结果，从而提高乳腺癌诊断的准确性。该技术有望为研究人员探索正常乳腺组织的发育特征及激素调控机制，以及表征乳腺恶性病变发生前的早期改变提供新的途径。观看本视频后，您应能够充分理解如何将我们介绍的采集与图像处理工具适配并应用于您先进的 MRI 扫描仪，以及如何评估乳腺扩散张量参数图并识别乳腺恶性病变。

tracking-mammary-architectural-features-detecting-breast-cancer-with

Research

科研

JoVE Journal

Medicine

医学

跟踪乳腺建筑特色与检测乳腺癌的磁共振扩散张量成像

测量心肌组织应变等参数的准确且实用的方法具有重要的临床价值，因为已有研究表明，与常用指标射血分数（EF）相比，应变是检测收缩功能障碍更敏感且更早期的指标。...

基于特征追踪的磁共振衍生心肌应变评估

用于评估纳米颗粒在细胞内生物分布的生物成像方法，适用于纳米制剂化合物的治疗和诊断监测。通过组织学共定位分析，本文所述方法具有高灵敏度和特异性。这些方法为从啮齿类动物到人类的应用提供了转化路径。...

用于诊断与治疗监测的纳米材料注册生物成像

在骨转移的发病机制中，血管生成是一个关键过程，因此成为成像与治疗的重要靶点。本文介绍一种部位特异性的乳腺癌骨转移大鼠模型，并描述了用于无创成像血管生成的策略 体内...

利用磁共振成像、容积计算机断层扫描和超声在乳腺癌骨转移裸鼠模型中进行血管生成的多模态成像

近年来，利用磁共振成像（MRI）进行细胞示踪受到了广泛关注。本实验方案描述了使用富含氟（19F）的颗粒对树突状细胞进行标记、将标记后的细胞体内应用，以及通过19F/1H...

监测树突状细胞迁移的19F / 1H磁共振成像方法

我们展示了在早期乳腺X线摄影中检测结构扭曲的方法。通过使用Gabor滤波器和相位图分析有方向性的结构，以检测放射状组织模式的病灶。对每个病灶通过表征放射状结构的度量指标进行特征提取和分类。这些方法有助于乳腺癌的检测。...

在既往乳腺X线摄影中检测结构扭曲 通过 定向模式分析

弥散张量成像（DTI）基本上是一种基于磁共振成像（MRI）的技术，用于识别 体内...

神经退行性疾病分析中的弥散张量磁共振成像

采用弥散张量成像（DTI）技术以尝试显示视觉通路的主要部分。目标是使用经美国食品药品监督管理局（FDA）批准的标准商用工作站，该工作站可用于日常常规操作，以期减少患者术后对视觉通路的损伤。...

视觉通路的弥散张量成像——白质纤维束与脑部病变

神经肌肉疾病通常表现出时间上动态变化、空间上异质性以及多方面的病理特征。本方案的目的是利用无创磁共振成像方法对这种病理特征进行表征。...

骨骼肌疾病的定量磁共振成像

本文介绍了一种应用扩散张量成像参数评估脊髓受压的实验方案。...

慢性脊髓受压中的弥散张量磁共振成像

活体显微镜的发现向临床转化面临的主要挑战是其对组织的穿透深度较浅。本文介绍了一种小鼠背部窗口腔模型，该模型可实现活体显微镜与临床适用的成像方式（如CT、MRI）的共同配准，从而进行直接的空间关联分析，有望加速活体显微镜发现成果的临床转化。...

用于转化癌症研究中活体显微成像与磁共振成像联合应用的背侧皮肤折叠窗室肿瘤小鼠模型

本方案描述了用于乳腺成像的微波成像设备的成像操作流程。...

微波乳腺成像的临床影像

学习使用磁共振扩散张量成像技术追踪乳腺微结构特征并检测乳腺癌的非侵入性方法。该技术通过分析水分子扩散系数和各项异性指数生成彩色参数图，无需造影剂即可识别恶性病变区域，为乳腺癌早期诊断和疗效监测提供安全可靠的影像学工具。

追踪乳腺结构特征并利用磁共振扩散张量成像检测乳腺癌

探索更多视频