张天龙,赵继荣,陈祁青,马同,王闻奇,赵宁,薛旭,郭培尧,蔡毅
1.甘肃中医药大学中医临床学院,甘肃兰州730030;2.甘肃省中医院,甘肃兰州730050;3.甘肃省中医院医学影像中心,甘肃兰州730050
各种疾病导致的慢性腰腿痛已经影响了绝大多数人的生活,腰椎间盘突出症(Lumbar Disc Herniation,LDH)为骨科多发病、常见病,其所致的慢性腰腿痛亦是导致诸多治疗棘手的病痛之一[1]。LDH 是指椎间盘内的组织(髓核或纤维环)移位于椎体后间隙,挤压神经根硬膜囊引起慢性腰腿痛或下肢麻木无力的一种症候。在美国,各种原因所致的慢性腰腿痛已经影响了大约1亿人,每年造成的直接经济消耗高达6 000 多亿美元[2]。循证医学研究表明,对于LDH 所致的慢性腰腿痛,药物治疗并不能彻底缓解疼痛症状,且部分药物副作用大;而手术治疗虽然能明显改善早期患者腰腿痛,但随访时间加长,患者仍出现不同程度的慢性腰腿痛,疗效亦不满意。近年来,关于静息态功能磁共振成像(Resting State Functional Magnetic Resonance Imaging,RS-fMRI)技术在疼痛医学领域得到了广泛应用,LDH 神经损伤所致的慢性腰腿痛的病理机制已经不再受限于病灶局部的病理机制探索[3],而着眼于全身各重要器官。大量研究均证实持续疼痛与大脑功能的改变及中枢损害存在一定相关性。然而,大部分临床医生及学者对于该疾病的认识还停留于传统的局部骨科疾患,并未认识到LDH 所致的慢性腰腿痛引起的神经损害等病理特征与脑部功能活动间存在密切关系。RS-fMRI技术作为研究二者关系的媒介之一,为探索LDH 致慢性腰腿痛的机制与脑功能间关系提供了可能。RS-fMRI 技术分析复杂大脑功能网络的方法众多,没有统一的标准,其实验结果的可靠性和有效性也受到了质疑。因此,本研究通过分析有关RS-fMRI脑功能连接分析方法在LDH 致慢性腰腿痛机制研究的应用,对LDH 致慢性腰腿痛的脑机制病理做一综述,从而为LDH致慢性腰腿痛的诊疗提供新思路。
RS-fMRI作为一种非侵入式的脑功能检测技术,具有无创、无污染、易定位、图形清晰、可重复、活体观察等优势[4]。RS-fMRI 对不同脑区血氧水平依赖(Blood Oxygen Level Dependent,BOLD)信号波动序列之间的时间相关性进行记录,并检测神经受激惹状态下血流、血容量和血红蛋白氧合作用的变化,以此来观察相关脑功能区的变化情况[5]。LDH 的病理生理特点是持续的慢性腰腿痛或间歇性疼痛强度加剧,可引起大脑区域BOLD(粗体)信号的时间同步性、空间分布波动性等指标的显著变化[6],其相应脑默认模式网络中杏仁核中心与皮质和皮质下区域的功能连接发生增强或减少。大量的重叠研究发现,LDH 所致的慢性腰腿痛患者的初级感觉运动皮层连接信号相比健康对照组更弥散,初级感觉运动皮层与岛叶连接强度的降低与更高的自发疼痛水平相关[7]。解剖学和电生理神经学相关研究也证明与LDH 相关的临床症状会引起脑结构改变,其皮层区域之间功能连接强度也会发生改变,包括躯体感觉皮层、岛叶、前扣带回皮质、背外侧/背内侧前额叶皮质和皮质下结构[8-9]。LDH所致的慢性腰腿痛刺激相应神经根,异常的神经传导效应激发多个脑区BOLD信号及代谢改变,引起脑功能区中疼痛网络连接改变,使疼痛调节功能受损,这一过程涉及脊髓、脑干、大脑皮层等多个结构。基于此,LDH 所致的慢性腰腿痛不再是局部病变,而是一种牵涉感觉系统、注意系统、记忆系统、默认网络等多种属性脑功能改变的病理状态[10-11]。
当下关于RS-fMRI 分析计算方法复杂多样,但主要是从脑功能分化和脑功能连接方面来探索LDH所致的慢性腰腿痛与复杂脑功能网络存在的相关性。功能分化主要集中于单一的磁共振信号源的研究,包括低频振幅、比率低频和局域一致性等方法[12]。低频振幅本身噪声大,相关脑室出现能量增大,一定程度上对RS-fMRI 造成了干扰,影响实验结果的准确性[13]。比率低频对BOLD(粗体)信号的振幅进行分析,不再进行信号波的过滤,其优势在于避免了低频震荡中相关脑区BOLD(粗体)信号的丢失,可以为研究区域提供更为全面的信息,其缺点是仅关注频域内的信息,而人脑中远隔脑区间的协作主要依靠时间同步性实现,所以该方法只注重脑局部的功能活动,而忽略了协助研究脑区调控的其他脑区[14]。局域一致性法是将给定体素的时间序列与其相邻较紧密的体素时间序列进行肯德尔和谐系数度量时间的同步性,以此来说明复杂的脑功能网络在实验分析过程中有诸多影响因素,如空间平滑大小、频域滤波和均一化方法等,实验结果的准确性仍有待商榷[15]。功能连接的重点是分析不同脑功能区域之间的相互协调作用,利用多个脑区的磁共振信号的联系紧密性和强度来发现其脑部的非正常变化。RS-fMRI 数据构建和脑功能连接分析的方法主要有3 种,分别是独立成分分析法(Independent Components Analysis,ICA)、基于种子点的方法和图论分析方法。以上方法在分析LDH 所致的慢性腰腿痛的脑机制方面发挥了巨大作用,并初步揭示了LDH 大脑功能连接的失衡[16-18]。现就其中的脑功能连接分析方法在LDH 所致的慢性腰腿痛中的应用进行讨论。
ICA已被确定为分析多元数据的基本方法,它将数据的线性分解(变换),无需选定输出数据和考虑不同成分的条件,在因子分析和主成分分析方法失败的情况下成功找到数据中潜在的原始成分和信号源[19]。近年来ICA 与脑影像学的结合可谓是重大发现,尤其是针对LDH 所致的慢性腰腿痛,ICA 将整个大脑的异常功能信号分解为许多相关空间网络及其相关的时间序列,从而使各组成部分的空间独立性最大化,然后,对这些空间网络及时间过程采用体素的方式进行分析,显示群体间和内部的连结性差异,以此找到慢性腰腿痛患者的异常脑功能网络区域。Baliki 等[20]最早利用该方法分析LDH 所致的慢性腰腿痛与大脑功能网络的联系,发现在慢性腰腿痛的持续刺激性下,丘脑后部、纹状体腹侧区和杏仁延展区有大量的异常活动,前额叶皮质内也有较大的活动,而顶叶和岛叶区未见任何活动,暗示当慢性腰腿痛持续时,它会激活涉及情绪、认知和运动的大脑区域。Kobayashi 等[21]通过RS-fMRI 对8 名健康受试者与6 名LDH 慢性腰腿痛患者进行脑功能网络改变的ICA,发现LDH 患者的脑功能疼痛调节区较健康受试者异常活跃,主要集中在前额叶、岛叶、后扣带皮质等区域,进一步验证了慢性腰腿痛引起的脑功能区的异常改变。李霁等[22]利用ICA观察LDH 患者针刺镇痛前后静息态脑默认网络(Default Mode Network,DMN)的功能磁共振信息变化,发现针刺镇痛前后的持续效应与DMN 存在相关性,DMN 中参与疼痛调节的区域发生了从负向到正向的变化,疼痛病理损伤得到了修复,这为临床评价针刺对LDH所致的慢性腰腿痛的疗效提供了功能影像学方面的客观依据。冯秀梅等[23]在RS-fMRI 下对坐骨神经受损的大鼠进行针刺观察大脑疼痛加工区的信息传导,发现初级大脑皮层、次级躯体感觉区、岛叶、前扣带回、前额叶皮质、杏仁核、基底神经节、小脑、海马、顶叶和颞叶内等与疼痛相关的皮质区内灰质与白质密度发生了改变,进一步为临床针刺治疗LDH 所致的慢性腰腿痛的效果提供了病理脑机制改变基础。ICA是盲源分离的一种方法,将混合的信号分离成潜在的信息成分,结合RS-fMRI 应用于LDH 所致的慢性腰腿痛脑机制领域,成功地分离出了LDH 的病理改变脑功能区,明确了持久的慢性腰腿痛刺激下的脑机制,极大地促进了循证医学的发展。然而,ICA也存在不足:(1)在分析LDH 所致的慢性腰腿痛的脑功能影像中,受到仪器噪声、冥想、头动等多个独立成分的干扰;(2)ICA 没有按照任何特定的顺序输出数据,这需要更复杂的方法来推断正确的排序;(3)大多数ICA算法都是基于局部优化法,从一个随机的初始点开始,在每次迭变时都试图增加目标函数,这样的算法不一定能在全局网络中找到目标函数的最优解,而且增加了实验结果的随机性和差异化[24-25]。总之,ICA的不足使提取的脑功能网络的精确性和有效性仍然值得探究,掌握现有方法的原理将有助于更好地进行数据分析,清楚该方法的缺点则能帮助其改进,最终能更好地服务于临床。
基于种子点的相关分析方法在分析LDH 所致的慢性腰腿痛脑功能数据方面具有简单、灵敏和易于解释等优势,因此在RS-fMRI 研究中得到了广泛应用。利用该算法,在脑功能区域中手动选择感兴趣的区域作为种子,然后对该种子的时间轨迹与图像中其他种子像素的时间轨迹之间执行Pearson 相关,得到相关图,以其显示脑功能上相应的区域,相应的时间轨迹之间具有高相关值(强相关)的区域可能在功能上相似,而具有低相关值(弱相关)的区域可能在功能上不相关,从而得出所需脑功能激活的区域[26-27]。Tagliazucchi 等[28]沿Baliki 路径从12例LDH 所致的慢性腰腿痛患者和20 例健康受试者中建立了8 个成熟的DMN,经过严格的科学分析后共选出27 个与LDH 所致的慢性腰腿痛最紧密的种子,发现连接性最强的种子位于扣带皮层、前躯、辅助运动区和颞上回,这与之前在全脑功能和解剖连接网络中连接的种子点识别相一致,这一研究表明,在LDH 致慢性腰腿痛患者中,DMN 中枢之间的相关平衡被破坏,而脑区的其他中枢却没有受到干扰,这一结果为LDH 作为影响DMN 的条件提供了论证依据,并且有助于缩小LDH 脑机制的研究范围。Hemington 等[29]在基于互补种子的分析中证明了LDH 致慢性腰腿痛患者后扣带回皮质区是大脑疼痛网络改变的枢纽,其结论揭示了慢性腰腿痛脑网络功能连接异常与临床症状的相关性。Zhang等[30]在受试者注射高渗盐水的实验条件下研究慢性LDH 患者DMN 和岛叶的静息态网络信息传导变化,发现持续的慢性腰腿痛刺激破坏DMN 的皮层区域与岛叶之间的功能连接,经过种子点分析得出的实验性LDH 的脑机制结果进一步提高了我们对LDH 致慢性腰腿痛的脑区异常神经活动和疼痛机制的认识。谢洪武等[31]选取慢性腰腿痛活动最活跃的后扣带回区作为感兴趣的种子,观察慢性腰腿痛患者脑功能连接默认网络的特异性改变,证实了慢性腰腿痛能够引起疼痛中枢的病理改变,为临床开展LDH 致慢性腰腿痛的脑功能区的精细化治疗提供了靶向依据。基于种子点的分析方法相对于ICA 提高了其慢性腰腿痛脑机制的精确性,能最大化地找出与慢性腰腿痛最紧密的脑网络调节区域[32],但是该方法也有其局限性,受限于预先选定的区域,只能得出感兴趣的脑网络区域,忽视了疼痛加工区其他脑网络的参与。
图论分析方法最早被发现于数学,以作图为依据,图论中的图是由给定的节点与任何相关两节点所勾画出的图形,这种图形采用自身独特的语言和符号来精准而简洁地描述各种社会领域的复杂网络关系,为数学家、物理学家以及哲学家提供了描述各种复杂网络的共同语言和研究平台[33]。随着近几年影像技术的快速发展,图论分析方法已应用于描述复杂脑功能网络,尤其应用于在RS-fMRI技术的基础上对复杂脑网络的描述,以揭示多种神经疾病复杂的脑病理机制[34-35]。Salvador 等[36]在2005年首次引入图谱构建脑功能的复杂网络,将大脑分割成不同的区域,计算脑区之间信号的时间相关系数,通过统计检验确定较强的脑功能连接来获取被试者的脑网络形态,该网络形态显现了最短路径长度和相对大的聚类系数,证实了小世界网络特征、高效的网络连接以及稳定的拓扑性质等[37],为研究上百个复杂脑功能区的信息处理提供了新方案。Baliki 等[38]在2013年利用图论分析方法将坐骨神经受损大鼠的大脑按照标准的图谱分割成了96 个解剖区域(每个半球48 个区域),在保证脑功能网络的可评估性和小世界网络特性的基础上,经拓扑度量的计算与分析,证明了啮齿动物全脑网络中的小世界网络特性的改变与持续的神经病理性疼痛相关,这一发现与最近人类大脑对LDH 致慢性腰腿痛的成像结果有着高度的相似之处。Zhang 等[39]通过长时间针刺LDH致慢性腰腿痛患者的穴位,对RSfMRI 获得的数据进行图论分析,从而获得脑功能网络的拓扑结构,证实在复杂的脑网络结构中存在小世界网络,并且发现小世界网络受到明显的慢性腰腿痛刺激的影响,长时间的针刺改变了脑功能网络中小世界特性。图论分析方法具有其他分析方法不可比拟的优越性,但也存在不足[40-41]:(1)临床研究中图论分析在脑功能区的切分、功能连接量化和阈值定义图形等方面存在极大的差异,而这些元素是图论分析的基础;(2)在不同的研究中,所报道的脑功能网络指标有很大的差异,多个决策点上因缺乏标准化的方法会产生乘法结果,任何两个研究使用相同的分割方案、功能连接定义、阈值策略和网络度量的可能性都非常低,这与区分不同的病理生理机制或识别诊断共性的努力背道而驰;(3)图论分析方法的多样性削弱了数据共享的价值,数据只有共享才能促进各类疾病的病理机制的最大化发展。为此,必须要有标准统一的方法进行数据收集,然后进行其精确化的分析,才能得出更可靠和真实的成果[42],使图论分析方法在各个领域实现解决问题的能力。
RS-fMRI 技术在探索神经损伤所致的慢性腰腿痛脑病理机制方面已成为医学热点,受到了国内外学者的热捧,其对LDH 致慢性腰腿痛的脑机制的研究颠覆了传统LDH 病理机制。过去,大多数LDH 所致慢性腰腿痛的神经影像学研究都是横断面的,并不能从根本上推断出所观察到的脑功能和结构变化。RS-fMRI 从纵向研究追踪LDH 所致慢性腰腿痛的脑功能复杂网络区域信号的改变,得到了具体的脑功能区和相关连接的病理变化,未来将会改变LDH 所致的慢性腰腿痛诊疗格局,涌现出更多以脑网络改变为基础的治疗方法。但是RS-fMRI 技术自诞生以来,在分析神经损伤大脑病理机制方面附带了诸多问题,其不断开辟新的分析方法应用于其中,诸多方法各有其优势,不能形成统一标准,不同的分析方法得出的结果也各有差异;除此之外,LDH 所致的慢性腰腿痛脑病理机制复杂,涉及疼痛调节的功能网络较多,很少是单一网络的控制,临床治疗需要考虑诸多脑功能调节区域;最后,基于RS-fMRI 临床研究LDH 所致的慢性腰腿痛样本甚少,实验动物研究也甚少,没有大队列的研究,这阻碍了该技术成果在临床的推广。今后研究的重点应该是在RS-fMRI基础上探索更为全面的分析方法对LDH 慢性腰腿痛进行大样本的实验研究和临床研究,以此增加其结论的可靠性、重复性和有效性,寻求适当有效的疗法来减轻慢性腰腿痛,达到正常大脑结构和功能水平修复的目的。