Introduction
磁共振成像(MRI)是一种非侵入性工具,它提供了一个窗口的脑和脊髓在健康和疾病。 MRI已经彻底改变了临床诊断,但它也适用于实验室研究的重要工具。神经损伤或疾病的动物模型提供了一个平台,了解病理生理学和推动发现疗法。在这份报告中,我们展示了MRI的应用脊髓损伤的大鼠模型来研究微观损伤1采用弥散张量成像(DTI)的潜在生物标志物。成像生物标志物的潜力发现将有助于患者脊髓损伤的诊断和管理。这些标记物都可能在他们的翻译中发现的疗法在临床前模型中发挥作用,并能观察或预后的临床环境。
DTI MRI是一种特殊形式的测量微观运动水分子( 即扩散)。 DTI已经在神经系统是特别有利的,由于轴突的存在,其中扩散是不成比例更快沿轴突比垂直于它们,这提供了关于他们的取向和微结构的组合物的信息。从DTI衍生标量指数,其中包括在组织内的总体扩散的量度,平均弥散(MD),和定向扩散的依赖性的量度,分数各向异性(FA)2,3,已经看到了广泛的应用中表征的微观结构神经系统在健康和疾病4。这些指标揭示微观组织的功能,是无形的,通过其他大多数MRI方法。以前的努力证明,DTI检测颈髓大鼠1以下胸SCI中远程微观结构的变化。英国贸工部的变化远离病变可能反映了怎样的整个脊髓水库池塘损伤,并具有潜在的继发性损伤的一个标志物。
活体成像的大鼠脊髓提出了一些独特的挑战。最值得注意的是,脊髓受呼吸运动的影响,并且需要小心注意,以尽量减少使用多种方法的运动。在过去的研究中,固定设备扫描5中去除脊柱的运动。为颈髓成像,我们利用在一个头支架和耳棒,其中衰减的形式物理约束,但不消除由呼吸运动。此外,我们利用一个自定义的呼吸门控方案来以有效的方式进行同步的图像采集与呼吸周期。这些修饰使去除否则所造成的呼吸6大型整体运动引起的伪影。 DWI是微观运动,包括脑脊液流动和血液脉动,和运动这些小源高度敏感contamination的缓解也由呼吸门控方案。此外,脊髓具有小的横截面面积,并且表示视场的一小部分。用于颈椎成像,其中脊髓位于动物的身体内深处,圆柱形射频线圈有足够的信号穿透是需要图像颈髓具有高分辨率。在视场的降低是由外体积抑制(OVS),它也用于取消或弃土,从组织中的信号的脊髓外实现。这种方法被称为扰流梯度或外卷抑制,也有助于减少这些组织中残留的动物的运动,脑脊液流动,或血液脉动任何污染。
脊髓的排列也可以被利用来简化成像协议。在白质(WM)脊髓轴突几乎都取向为平行于脊髓的主轴线。日我们,而大脑的DWI需要至少沿6个方向,以确保该结果不依赖于磁铁(这个过程被称为扩散张量成像)内的位置的测量,在脊髓测量能够获得仅沿2个方向平行,并垂直于所述线7,8,以下简称为纵向和横向上。因此,扩散度及其他参数沿2个方向分别测定推论,并允许进入在健康和疾病或损伤的组织的微结构。
Protocol
注:伦理学声明:威斯康星医学院的机构护理和使用委员会(IACUC)和克莱门特J. Zablocki VA医学中心批准的所有程序。
1.动物准备和监控
- 麻醉大鼠在感应腔,用5%异氟烷在医用空气。当翻正反射是不存在,并且挤压后爪不产生撤回反射,减少麻醉至2%,在头部的第一俯卧动物转移到扫描仪床。保持2%异氟醚通过在整个过程的鼻锥设备,并保持医用空气以约1升/分钟的流速。取少量润滑软膏,以老鼠的眼睛,以避免损伤角膜,而在麻醉下。
- 放置一个呼吸监测皮带围绕牢固老鼠的躯干。皮带连接到一个呼吸门控系统。前推进老鼠放入扫描仪孔,CHEC k是呼吸监测计算机,以确保呼吸周期是清楚和一致。必要时调整皮带,因为这一步是必要的图像质量。
- 监视和通过直肠探头和温暖的空气加热系统保持动物的体温在37℃。通过调整麻醉介于1.2和2%的水平维持每分钟30-45次呼吸的呼吸速率。
- 定位在所述头部保持器的大鼠的咬入杆和螺钉在耳棒( 图1),和头滑入一正交线圈体积直至颈椎被定位在线圈的中心。
注:老鼠的肩膀可以防止进一步恶化到线圈。 - 推进大鼠和配套支架插入扫描仪孔。如果适用的话,调整调谐和线圈到适当频率的匹配电容器,并根据由所述线圈供应商提供的说明进行操作阻抗。
注:这里所描述的方法使用了9.4 t水平孔小动物的制度,但适用于小动物磁共振成像系统等领域的优势。
- 使用MRI系统的自动化程序来检测共振频率的,迭代地改善磁场的接收器增益的均匀性(匀场)中,射频功率的校准和调节。
- 使用该系统的软件界面,获得一个默认的三平面定位扫描,以确保正确定位。
- 点击“新建扫描”,选择tripilot,并点击“交通灯”,以获取图像。
- 确保颈椎的中心与所述磁体的两个中心和MRI线圈的中心对准。居中脊柱磁铁,在推或拉的摇篮,重新获取验证搜索扫描。
- 要调整positioÑ颈椎相对于所述MRI线圈的,从磁铁重新定位取出支架。如果有必要,重复这个过程,直到立场是一贯的。如果动物被重新定位,重复步骤2.1。
- 添加新的回波平面扩散加权自旋回波序列(DtiEpi)到当前摄像协议。
- 配置和使用,除了以下的默认设置与DWI序列获得弥散加权成像:
- 打开切片位置的图形界面,以规定12片的厚度为0.75毫米。定向片垂直于颈髓的主轴线。确保不同动物之间或跨越使用小脑作为内部参考的基部不同的成像会话一致切片定位。
- 设置饱和度乐队“开”。第4位的饱和带与脊髓外部的厚度为10毫米,以尽量减少从这些组织的信号,并减少其潜在的诱导构件( 图3)。设置呼吸门控'的'('触发模块“)来。
注:自定义呼吸门控要求的脉冲序列编程的知识和经验。如果这是不可用,一个解决办法是,减小条带的数量3-5和对TR至1秒,以确保所有的条带被在中间的动物的呼吸获得。重复该全序列与片的另一子集,以获得颈髓的全面覆盖。 - 单击工具箱图标,然后单击“编辑方法。”EPI的段数设置为4.更改相位编码方向左右。其他默认设置应该是:回波间隔= 0.3234毫秒,每EPI段总回声列车长度= 32。
注:相位编码中设置的左右方向,而不是前 - 后会减少运动的污染其他结构。 - 使用以下geometriCAL设置。矩阵大小= 128×128,并且在平面视场的= 25.6 X25.6毫米,导致平面内的空间分辨率= 0.200 X0.200毫米。确保层厚= 0.75毫米。切片为了='交错',层间距= 0毫米。
- 使用下面的扩散加权设置:DW测量模式='DW对比度',扩散梯度的持续时间(δ)= 7毫秒,扩散梯度分离(Δ)= 12毫秒,b值数= 8,所需的b值= 0 ,250,500,750,1000,1500,2500 3500毫米/秒2,扩散方向= 2,扩散加权方向= [1 0 0]号和[0 0 1](制成为在平面平行且正交于脊髓轴)。
注:使用这些设置,我们实现了b值高达3500秒/毫米2。硬件规格和其他系统性能特征可能限制b值,由于扩散梯度的持续时间(δ)和扩散梯度分离(Δ)是依赖于克radient的表现,这对我们的系统有:(最大梯度强度:440吨/米,最大压摆率:3440 T /米/秒)。用于测量峭度的,2 b值,具有至少2000秒/毫米2的较高b值,被推荐的。 - 使用下面的时间设置。回波时间(TE)= 27毫秒(通过输入0设置为最小),重复时间(TR)= 1800毫秒。
- 收购准备序列。上面列出的参数,总采集时间是约25分钟。
- 在所有的扫描,监测呼吸门控的软件,并调整“触发”之间的延迟时间(软件检测呼气的)和信号对MRI系统,使得仅出现在呼吸周期的静态(静止的)部分( 图收购2A,灰线的稳定部分)。 100-400毫秒之间的触发延时需要根据动物的呼吸模式。这将有助于减少ARTI与呼吸运动( 图3e)发生的事实。
- 如果可以,用定制的“反向光点”设置为“ON”,这就要求采集时间的25分钟重复上述步骤。
注意:如果自定义的“反向标志”序9(第3步骤期间所需的易感性伪影校正)是不可用的,只有一个单一的EPI相位编码方向是可能的,而相反的光点序列修饰允许相位编码方向的选择(右-to-向左或从左向右)。 - 当成像完成后,从支架上取下动物,并返回到它的笼子里。不要让动物无人看管,直到它恢复足够的意识,保持胸骨斜卧。
3.图像处理
- 从系统中DICOM格式导出数据直接从系统(首选)或使用自定义或锡尔将数据转换为NIFTI格式D-第三方软件。
- 进行敏感性伪影校正。
- 提取B =从每个0卷扫描到一个文件中,使用带有FSL或其他MRI软件包工具。是必需的一个文件的每个相位编码方向。
注意:例如,如果每次扫描由8扫描不同b值与扩散加权在横向方向上,其次是8扫描的长度方向扩散加权的,则该图像文件中包含B = 0的扫描中的1 次和第第 9卷,并且可以提取并组装用下面的壳代码:
fslroi $ {}起来_dwi_masked.nii.gz temp1中0 1
fslroi $ {}起来TEMP2 _dwi_masked.nii.gz 8 1
fslroi $ {}下降_dwi_masked.nii.gz TEMP3 0 1
fslroi $ {}下降_dwi_masked.nii.gz temp4 8 1
fslmerge -t blip_both temp1中TEMP2 TEMP3 temp4
(其中,在这种情况下,$向上和向下$是与正常的扫描和反向相位编码方向,分别地)。 使用FSL 10,11的'充值'命令创建降低图像失真文物修正文件。该修正适用于将用于创建参数映射的原始的DWI图像。
注:使用说明该命令可以在这里找到http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/TOPUP/TopupUsersGuide 。示例代码以在这种情况下使用的命令如下:
充值--imain = blip_both_nlmf_b0images_masked.nii --datain = .. / topup_data.txt --config =。/ b02b0_ratspine.cnf --out = topup_splines_nlmf --iout = $出去 - 详细--logout = topuplog.log
dwiup =`LS $ {}上升* dwi_nlmFilt.nii`
dwidown =`LS $ {}下来* dwi_nlmFilt.nii`
applytopup --imain = $ {} dwiup,$ {} dwidown = --datain .. / topup_data.txt --method =江淮--inindex = 1,$ IND --topup = topup_splines_nlmf --out = DWI _ $ {出} -v
复制和编辑$ {} FSLDIR的/etc/flirtsch/b02b0.cnf默认文件的大鼠旋人帘线由通过10倍减少的每个在--warpres和--fwhm线的数值。 - 提取B =从每个0卷扫描到一个文件中,使用带有FSL或其他MRI软件包工具。是必需的一个文件的每个相位编码方向。
注:从脊髓等程序段的投资回报已被记录在案14,15
fslstats -t DWI_corrected.nii.gz -k GM_mask.nii.gz -M
- 在一个数值计算程序复制的b值作为8 b值的向量。的b值在横向和纵向方向是相同的。如果可能的话,将有效b值,而不是标称b值,应该从扫描仪,它是从步骤2.3.5为“有效果B值”中列出的参数窗口获得。
- 使用数值计算程序的曲线拟合工具箱,以适应信号与b值数据为thË通过在命令提示符下键入cftools所需的模型。要做到这一点,请单击“数据...”,然后选择以Y型数据信号矢量和b值x数据。点击“配件...”和“拟合类型”下选择“自定义公式”,然后单击“新建”和“一般方程”,进入一个方程拟合。
S0 * EXP(-x。* D)“(1)
S0 * EXP(-x。* D +(1/6)*(X * D)^ 2 * K)“(2)
AI =( 平负-D T)/(平负+ D T)(3)
这是类似的分数各向异性(FA)从DTI模型计算。峰度的各向异性指数也可以使用横向和纵向峭度代替扩散率计算出来的。
注意此方法给出的模型参数如K T,D T,等等。另外,也可以使用曲线拟合工具箱的命令行操作的每个体素的脊柱内,以从模型创建的每个参数的值映射。替代拟合方法可以用于和其他地方详细描述。17
Representative Results
正确的步骤,以尽量减少运动伪影导致大鼠脊髓型颈椎病的高品质的弥散加权图像。使用自定义呼吸门控( 图2),饱和从组织不希望的信号的脊柱外侧( 图3B& C),以及磁场磁化率失真校正产生扩散加权图像,如在图4和5。不当或未选通的图像会导致重影( 图3E)的形式文物,而正确的门是免费的文物。
横跨12片的扩散加权图像的视觉检测揭示了涉及其微观脊髓的功能。具体而言,在组织导致在扩散加权图像更大的信号损耗,这是加剧具有更大扩散加权(b值)快扩散。随着扩散权重进行perpendicular脊髓轴,白质沿该线的周边显示明亮,因为扩散是缓慢的,限制在垂直于轴突。与此相反,该线的中心区域内的灰质出现暗的,因为它是由轴突和细胞体未全部沿单一方向排列的。相比较而言,扩散加权在并联方向导致白质与较暗的外观,因为扩散是快速沿轴突,而灰质是相对较亮。值得注意的是,单独的扩散加权图像被示出为不同的b值,由于平行和垂直的方向上具有不同b值白质和灰质之间的最佳对比度是重要的。
综合所有使用的数学形式主义扩散加权图像就可以显示扩散参数的地图。从白质和灰质的平均信号被暗算的DIF融合权重系数(b值)为平行和垂直的方向。这个量化数据强化了图4所示的扩散加权图像,具体而言,白质具有一个在扩散加权(纵向或横向)的方向上有很强的依赖性,而灰质是较少依赖于方向。同样地,使用等式扩散峭度的产量的扩散参数( 图6B),它突出这一相同的依赖定量映射拟合信号中的每个体素。白质具有高程度的各向异性的两个扩散(AID)和峰度测量(AIK)。因此,横向的扩散和峭度表明,从组织学研究已知脊髓底层的显微结构。这些扩散参数,这将在现场获得的,但麻醉动物,反映微观组织性质如轴突密度和直径。在改变所造成的伤害和疾病本身的措施将是无创性评估损害的后果,并承诺治疗的效果非常有用。因此大鼠脊髓型颈椎病的弥散加权成像可能成为脊髓损伤和脊髓疾病的临床前研究的一个工具。
图1:线圈和持有人的脊髓型颈椎病的MRI设计的自定义正交卷线圈(多蒂科学公司)用于图像颈椎具有高灵敏度和均匀性。麻醉和医用空气是由所指示的气体口递送到鼻锥,它适合舒适周围的大鼠的鼻子。呼出和过量的气体被轻微真空下排气管线捕获。鼠的头被固定在咬合杆周围放置门齿和耳杆放置德利奇ately内耳道。其他生理监测组件,包括呼吸监测器和温度传感器未示出。
图2:呼吸门控方案 。一个典型的呼吸迹线(灰色)和触发器(红色)从选通单元被示意性示出(A)中 。在典型的实现选通(B)的,单触发器用于获取所有切片(垂直线;在这里12示出)在时间均匀间隔的重复时间(TR)内。如果TR超过典型期间,可能会发生呼吸期间几片,并容易受到运动(红色)。在修改后的方案(C)中 ,片的一个子集被触发(此处6中示出),接着在延迟之后迅速获得与随后的触发后获得的其他条带。有效中,TR为通过重新排列序列内的延迟相同的两种方案之间。
图3:MRI片定位,饱和度带,和运动控制轴十二片被安排的侦察图像(A)最前片位于距脑干和小脑的交点一致的距离上。饱和带(B)的加入来消除感兴趣的区域之外的干扰信号。不扩散权重(C)和一个带有扩散加权(D)所采用的定制选通方案的图像清晰地显示了线的解剖,是免费的文物。与非优化的门控方案或不当呼吸门控,扩散加权图像显示伪影(E),为丢失电源线内的信号,或多个“鬼”的线,将损坏的后续分析,外请点击此处查看该图的放大版本。
图4:代表弥散加权像使用文本描述的优化,高品质的弥散加权图像与扩散加权横向应用(A)和纵向(B)脊髓主轴获得。不同b值示为每个方向提供用于说明性目的的白色和灰色物质之间的最佳对比度。对于每一个方向或b值,所有12个片约90秒被收购。加载/ 52390 / 52390fig4large.jpg“目标=”_空白“>点击此处查看该图的放大版本。
图5:反相编码校正过程中左栏显示了单片成像与DWI序列在本协议(“昙花一现达”的形象)概述。中间一列给出了序列获得第二次用“反向的blip”设置为'上'。请注意如何将显得捉襟见肘,在第一图像特征出现压缩在中间一栏。右列显示了使用充值纠正扩散加权图像。顶行是所述非扩散加权图像,中间行是与在横向方向上施加扩散加权的例子,且最底行是与在纵向方向上施加扩散加权的例子。//www.jove.com/files/ftp_upload/52390/52390fig5large.jpg“目标=”_空白“>点击此处查看该图的放大版本。
图6:计算扩散和峭度的地图的归一化信号(图像强度)绘制(A),因为扩散加权(b值),用于在横向(T)和纵向(L)的扩散编码方向的函数。高质量的地图扩散(D)的(B),峰度(K)和各向异性(AI)是从每个像素的信号,计算并显示脊髓组织的独特功能。具体而言,就是在白质和灰质之间的参数,以及白质区域的区域差异明显的区别。 请点击此处为viEW这个数字的放大版本。
Discussion
这里概述的方法可以在体内提供的大鼠脊髓高质量扩散加权图像。图像质量取决于许多因素,但脊髓具有很重要的几个独特的问题。
运动是一个重要的问题,即如果不纠正,将导致不能使用的图像。因此,它需要MRI会议期间仔细监测。如果图像伪影的初始扫描是与运动相一致观察,停止收购,并采取措施消除文物,因为这些是很难在后期处理来去除。确保呼吸道计算机从呼吸监测设备接收到一个强大的,定期的信号。呼吸带,可能需要进行调整以获得正确的张力,提供了一个一致的信号,但不限制动物的呼吸。维持麻醉在任何时候的适当水平; 1.5-2.0%,异氟醚已在我们的experien被使用CE。同样地,减少在动物和脊柱的整体运动是提供无伪影的图像的另一个重要方面。不像人类脊髓,其经历引起脑脊液脉动相关的心脏周期显著运动,在啮齿动物脑脊液脉动与呼吸周期18主要关联。虽然很难完全消除在帘线所有运动中,特别重要的是减少了运动的程度成为可能,这通常是通过试验和错误来实现。此外,大鼠各种神经损伤或病症可具有异常呼吸速率或其他生理并发症,可能需要的本文所概述的程序适应。
该修改用于呼吸门控脉冲序列,连同定制用于该目的的图像重建程序,使失真最小而引起不均匀磁场以致无法祛瘀的影响通过在MRI系统上执行的调整编
同样地,图像质量取决于成像时间的持续时间。在我们的例子中,限制扩散的加权的仅沿两个方向上的数目使在总成像时间的减少。这种方法的一个限制是,它不再是全张量分析(DTI),这是标准的许多其他研究兼容。或者,使用较少的平均值和多扩散方向或b值可以允许更好的表征,同时保持相同的采集时间。以前的研究已经表明,在2方向的方法提供了与6-方向(DTI)的方法19一致的信息,但必须小心,以确保切片(和扩散方向)的取向准确地沿着并垂直于该线。然而,获取多个b值允许更好表征和峭度的数学拟合,推荐在使用单一的B-VALU的即此外,全序列,重复用颠倒的相位编码方向这减少磁场敏感伪影的效果,并提高了通过平均化图像的整体质量。最后,在我们的协议中使用的图像的分辨率提供的白质和灰质明确分离。具有更高分辨率的图像是可能的,虽然这往往出现以牺牲更长的扫描时间或更多神器的潜力。
在射频线圈,脉冲序列,和后处理方法的改进都将有改善脊髓成像在该方法的将来的修改的效果。例如,表面线圈可能有利于类似于在小鼠中观察到的改善的图像质量。20,这些措施有被作为生物标志物用于脊髓损伤的临床诊断和管理有用的高可能性。
Disclosures
这篇文章出版费用由布鲁克公司被部分赞助。
Acknowledgments
我们感谢凯尔斯特赫利克,娜塔莎·威尔金斯和马特Runquist实验协助。通过研究和教育计划基金,推进更健康养老威斯康星州威斯康星医学院的一个组成部分,与克雷格·H·尼尔森基金会资助。
Materials
Name | Company | Catalog Number | Comments |
Small animal imaging RF coil | Doty | SAIP400-H-38-S | |
Respiratory gating system | SA Instruments | 1030 | |
MR scanner | Bruker | Biospec 94/30 USR |
References
- Jirjis, M. B., Kurpad, S. N., Schmit, B. D. Ex Vivo Diffusion Tensor Imaging of Spinal Cord Injury in Rats of Varying Degrees of Severity. J Neurotrauma. 30 (18), 1577-1586 (2013).
- Basser, P. J., Mattiello, J., Lebihan, D. Estimation of the Effective Self-Diffusion Tensor from the NMR Spin Echo. J Magn Reson. 103 (3), 247-254 (1994).
- Basser, P. J., Mattiello, J., LeBihan, D. MR diffusion tensor spectroscopy and imaging. Biophys J. 66 (1), 259-267 (1994).
- Song, S. -K., Sun, S. -W., Ju, W. -K., Lin, S. -J., Cross, A. H., Neufeld, A. H. Diffusion tensor imaging detects and differentiates axon and myelin degeneration in mouse optic nerve after retinal ischemia. NeuroImage. 20 (3), 1714-1722 (2003).
- Beckmann, N., Bruttel, K., Urban, L., Rudin, M. Signal changes in the spinal cord of the rat after injection of formalin into the hindpaw: characterization using functional magnetic resonance imaging. Proc Natl Acad Sci U S A. 94 (10), 5034-5039 (1997).
- Le Bihan, D., Poupon, C., Amadon, A., Lethimonnier, F. Artifacts and pitfalls in diffusion MRI. J Magn Reson Imaging. 24 (3), 478-488 (2006).
- Ford, J. C., Hackney, D. B., et al. MRI characterization of diffusion coefficients in a rat spinal cord injury model. Magn Reson Med. 31 (5), 488-494 (1994).
- Clark, C. A., Barker, G. J., Tofts, P. S. Magnetic resonance diffusion imaging of the human cervical spinal cord in vivo. Magn Reson Med. 41 (6), 1269-1273 (1999).
- Mohammadi, S., Nagy, Z., Hutton, C., Josephs, O., Weiskopf, N. Correction of vibration artifacts in DTI using phase-encoding reversal (COVIPER). Magn Reson Med. 68 (3), 882-889 (2012).
- Andersson, J. L. R., Skare, S., Ashburner, J. How to correct susceptibility distortions in spin-echo echo-planar images: application to diffusion tensor imaging. NeuroImage. 20 (2), 870-888 (2003).
- Smith, S. M., Jenkinson, M., et al. Advances in functional and structural MR image analysis and implementation as FSL. NeuroImage. 23, S208-S219 (2004).
- Behrens, T. E. J., Woolrich, M. W., et al. Characterization and propagation of uncertainty in diffusion-weighted MR imaging. Magn Reson Med. 50 (5), 1077-1088 (2003).
- Cook, P. A., Bai, Y., et al. Camino: Open-Source Diffusion-MRI Reconstruction and Processing. 14th Scientific Meeting of the International Society for Magnetic Resonance in Medicine. , 2759 (2006).
- Kim, J. H., Tu, T. -W., Bayly, P. V., Song, S. -K. Impact Speed Does Not Determine Severity of Spinal Cord Injury in Mice with Fixed Impact Displacement. J Neurotrauma. 26 (8), 1395-1404 (2009).
- Tu, T. -W., Kim, J. H., Yin, F. Q., Jakeman, L. B., Song, S. -K. The impact of myelination on axon sparing and locomotor function recovery in spinal cord injury assessed using diffusion tensor imaging. NMR Biomed. 26 (11), 1484-1495 (2013).
- Jensen, J. H., Helpern, J. A., Ramani, A., Lu, H., Kaczynski, K. Diffusional kurtosis imaging: The quantification of non-gaussian water diffusion by means of magnetic resonance imaging. Magn Reson Med. 53 (6), 1432-1440 (2005).
- Veraart, J., Sijbers, J., Sunaert, S., Leemans, A., Jeurissen, B. Weighted linear least squares estimation of diffusion MRI parameters: Strengths, limitations, and pitfalls. NeuroImage. 81, 335-346 (2013).
- Budgell, B. S., Bolton, P. S. Cerebrospinal Fluid Pressure in the Anesthetized Rat. J Manipulative Physiol Ther. 30 (5), 351-356 (2007).
- Tu, T. -W., Kim, J. H., Wang, J., Song, S. -K. Full Tensor Diffusion Imaging Is Not Required To Assess the White-Matter Integrity in Mouse Contusion Spinal Cord Injury. J Neurotrauma. 27 (1), 253-262 (2010).
- Kim, J. H., Song, S. -K. Diffusion tensor imaging of the mouse brainstem and cervical spinal cord. Nat Protoc. 8 (2), 409-417 (2013).