研究文章

利用相位对比和四维流磁共振成像对阿尔茨海默病中脑脊液-血管关联的定量评估

56 次观看

DOI:

10.3791/71768

2026年9月3日

 ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  , 

通讯作者: Dingtai Wei <wdtai83.@163.com>

* These authors contributed equally

本文内容

摘要

相位对比磁共振成像(PC-MRI)和四维血流磁共振成像(4D-Flow)揭示了阿尔茨海默病中脑脊液(CSF)动力学与脑血管系统之间关系的改变,反映了潜在的类淋巴系统功能变化,并提供了一种无创、定量的方法,可用于潜在的辅助生物标志物及疾病监测。

摘要

从神经病理学角度来看,迄今为止阿尔茨海默病(AD)最广为人知的病理改变主要是tau蛋白的过度磷酸化,这会导致细胞内神经原纤维缠结(NFTs)的形成以及淀粉样蛋白沉积——即细胞外淀粉样蛋白(Aβ)及其他蛋白聚集体的累积,从而引发神经原纤维缠结。类淋巴系统由动脉搏动驱动,通过血管周围间隙(PVS)清除包括淀粉样蛋白-β在内的代谢废物。在AD中,已分别观察到血管功能障碍和慢波睡眠异常,二者均被认为可能损害类淋巴系统的清除功能。然而,这些因素之间的因果关系尚未得到证实,而所提出的自我强化循环——即血管功能衰竭可能加剧蛋白沉积并进一步导致血管损伤——仍属推测。本研究采用二维电影相位对比磁共振成像(PC-MRI)和四维流磁共振成像(4D-Flow),定量分析了脑脊液(CSF)在中脑导水管(CA)中的流动以及颈内动脉(ICA)中的血流等血流动力学参数。在AD患者中观察到脑脊液-血管关联的改变,提示类淋巴系统可能参与了疾病进程。脑脊液与血管流动参数的整合为AD的早期检测与监测提供了一个有前景的生物标志物框架。PC-MRI与4D-Flow揭示了阿尔茨海默病中脑脊液-血管关联的改变,反映了潜在的类淋巴功能变化,为无创、定量的辅助性生物标志物开发及疾病监测提供了可能方法。

引言

阿尔茨海默病(AD)是最常见的痴呆类型,是一种以记忆丧失和认知功能障碍为特征的进行性、不可逆的神经系统疾病1,2。随着全球发病率迅速上升3,早期检测与干预已在全球范围内受到广泛关注。

由于阿尔茨海默病(AD)具有复杂的病理生理机制,其确切病因至今尚不明确。主要的神经病理学特征包括由tau蛋白过度磷酸化导致的细胞内神经原纤维缠结以及细胞外的β-淀粉样蛋白(Aβ)斑块1。越来越多的证据表明,Aβ清除障碍而非其过度产生,是驱动淀粉样蛋白沉积、促进疾病发生和进展的主要因素4。大脑中Aβ清除的一个关键通路是类淋巴系统(glymphatic system),这是一种由水通道蛋白-4(AQP4)介导的血管周围废物清除网络5,6。类淋巴功能障碍会损害脑脊液(CSF)通过血管周围间隙向脑实质的流入,从而阻碍脑脊液与组织间液(ISF)之间的交换,并影响包括Aβ在内的代谢废物的清除7。通过正电子发射断层扫描(PET)或脑脊液分析可检测到AD生物标志物的改变,这些改变可能在临床症状出现前数年即已存在8。在AD模型小鼠中的研究表明,类淋巴系统的异常出现在Aβ沉积之前9,10,凸显了其通过血管周围通路在脑内废物清除中的重要作用。

脑脊液动力学与脑血流动力学对于理解衰老和认知功能障碍至关重要。脑脊液维持中枢神经系统的稳态,促进代谢废物清除,并提供机械保护作用。脑血流的流动模式与多种因素相关,包括脑组织的顺应性、血管完整性、颅内压的变化11以及皮层tau蛋白沉积12,13。脑血流量、动脉搏动和血管舒缩运动是脑血管健康的关键指标,尤其在老年人群和阿尔茨海默病(AD)患者中具有重要意义14。异常的脑脊液动力学会导致多种神经系统疾病,包括阿尔茨海默病和正常压力脑积水(NPH)15,16。近期研究发现,阿尔茨海默病的严重程度可能与脉络丛通透性的改变有关,提示脑脊液循环障碍在AD的发生发展中起一定作用17

4D 流动 MRI 和 2D 电影相位对比(PC)MRI 等先进成像技术显著提升了评估这些动态变化的能力。PC-MRI 可无创量化一个心动周期内颅内动脉、静脉及脑脊液(CSF)流动的速度、方向和信号强度18。4D 流动 MRI 能够捕捉随时间变化的三维血流,从而实现对血流动力学的全面评估,包括管壁剪切应力和脉搏波速度19。早期基于 2D 电影 PC-MRI 的基础研究已证实其在量化脑脊液和血液流动动力学方面的价值1,20,21,22。例如,通过 2D PC MRI 测量到的老年人群中升高的脑动脉搏动已被证实与微血管损伤和认知功能下降相关23。近期,采用 4D 血流 MRI 的研究进一步揭示了脑血管动力学与认知障碍之间的关系,部分研究表明动脉搏动增强与认知功能下降相关,尤其在阿尔茨海默病(AD)和轻度认知障碍(MCI)人群中更为显著24。同时,实时 PC-MRI 的最新进展为捕捉瞬时流动现象提供了更深入的视角25

本研究旨在利用相位对比磁共振成像(PC-MRI)量化脑导水管中脑脊液的流体动力学标志物以及颈部大血管的血流动力学指标,探讨阿尔茨海默病(AD)、轻度认知障碍(MCI)与正常对照(NC)组之间的差异,并分析这些参数与认知评分的相关性。整合这些影像学方法可能有助于提升神经系统疾病的诊断能力,并为靶向治疗策略提供依据15

方案

参与者
2022年11月至2023年11月期间,参与者从宁德市医院神经内科及宁德师范学院招募。根据纳入/排除标准,共纳入35例患者(16例阿尔茨海默病患者,19例轻度认知障碍患者)和25例健康老年对照者(NC)。所有参与者及其认知障碍患者的家属均签署了知情同意书。本研究遵循赫尔辛基宣言的原则进行。宁德市医院伦理委员会批准了本研究(批准号:20200901)。

纳入标准为年龄在50至85岁之间,愿意接受磁共振成像(MRI)检查并签署知情同意书,以及能够在家人协助下或通过多个检查阶段完成影像学检查。排除标准包括存在脑肿瘤;严重心血管疾病;精神病史;既往创伤性脑损伤;主要神经系统疾病,如中枢神经系统感染、帕金森病、亨廷顿病、癫痫或多发性硬化;严重全身性疾病;不可移除的金属植入物;严重影响认知评估的重度感觉障碍;图像伪影;以及左利手。

阿尔茨海默病(AD)、轻度认知障碍(MCI)和正常认知(NC)的诊断由两名神经科医生根据美国国家衰老研究所-阿尔茨海默病协会(NIA-AA)标准确定。可能的AD诊断需符合认知和功能进行性下降的标准;MCI需存在客观认知损害但日常功能保持完好;NC则要求无认知或神经系统异常。所有受试者的诊断均经共识达成。生物标志物确认在本研究中并非普遍可及,这是本研究的一项局限性。本研究使用的所有工具和材料均列于材料表中。

临床数据
收集了人口统计学和临床数据,包括年龄、性别、教育程度、病史(高血压、糖尿病、高脂血症)以及体格检查结果。采用简易精神状态检查量表(MMSE)和蒙特利尔认知评估量表(MoCA)评估认知功能。

图像采集
使用专用头颈线圈在3.0T磁共振成像(MRI)扫描仪上进行成像。成像方案包括T2加权液体衰减反转恢复(T2-FLAIR)、扩散加权成像(DWI)、矢状位T1加权成像以及相位对比MRI(包含2D电影相位对比和4D流序列)。所有受试者均采取仰卧位,头先进入扫描仪,使用泡沫垫固定头部,并佩戴耳部保护装置以减少运动伪影。

二维电影相位对比 MRI
根据矢状面 T1 加权定位图像确定,在中脑导水管垂直方向上采集二维电影相位对比序列(快速 PC 电影切片)。成像参数如下:视野,24 × 24 cm;层厚,5.0 mm;重复时间/回波时间,20.2/6.1 ms;矩阵,288 × 256;像素大小,0.83 × 0.94 mm;带宽,31.25 kHz;激发次数,1.0;翻转角,20°;编码速度,12 cm/s;外周门控;每个心动周期采集 20–30 个时相。选择该编码速度以充分捕捉导水管脑脊液流动,避免混叠伪影。采集时间根据受试者心率在 1 至 3 分钟之间变化。

四维血流磁共振成像
采用一种4D血流序列(Sag 4D Flow)对颈部颈动脉区域进行扫描。成像参数如下:视野范围,24 × 24 cm;层厚,1.0 mm;层间距,2.0 mm;层数,40–50层;重复时间/回波时间,6.6/2.8 ms;矩阵,180 × 180;像素大小,1.33 × 1.33 mm;带宽,62.50 kHz;激发次数,1.0;翻转角,15°;编码速度,50 cm/s;每个心动周期采集20–30个时相。选择50 cm/s的编码速度可在充分捕捉颈动脉血流的同时,最大限度地减少混叠伪影。采集时间根据心率不同,介于1至3分钟之间。所有图像均使用专用的后处理软件进行分析(材料表)。

图像处理与分析
脑脊液相关参数的分析
使用专用的后处理软件对二维电影相位对比(cine PC)的幅值图和相位图进行分析(图1A–D)。两名放射科医师(每人具有三年神经影像经验)在中脑导水管最大横截面积处,沿其外缘独立勾画感兴趣区(ROIs),该位置首先在矢状位T1加权定位像上确定,并在幅值图像上确认。每位放射科医师对每个ROI进行两次勾画,两次间隔至少两周。采用软件内置的静止组织校正算法进行背景相位校正,该方法通过从速度编码的相位图像中减去静止组织区域(置于导水管邻近的中脑实质内)的平均相位,以最小化涡流和麦克斯韦项的影响。为每位受试者生成以下参数:最大收缩期流速(收缩期峰值流速,PSV)、最大舒张期流速(舒张期峰值流速,PDV)、净搏动容积以及流量曲线。记录导水管的横截面积(图2)。采用双因素随机效应模型的组内相关系数(ICC)评估评分者间和评分者内可靠性,所有参数的ICC值均超过0.90。

脑部MRI扫描对比;矢状位、轴向视图、神经影像学分析。
图1:脑脊液流动成像的采集与处理。A)矢状位T1加权定位图像,显示2D电影相位对比切片的定位,该切片垂直于中脑导水管(虚线)。(B)电影2D序列的幅度图像,显示中脑导水管的解剖结构。(C)电影2D序列的相位图像,显示脑脊液(CSF)流动信号强度。缩写:CA = 中脑导水管;CSF = 脑脊液。请点击此处查看此图的放大版本。

速度与时间关系图;数据点显示运动分析;物理数据分析。
图2:典型的脑脊液流动速度曲线。 时间-流速曲线显示在一个心动周期内脑脊液呈双相性、准正弦型的流动模式。该曲线展示了收缩期和舒张期流速的变化情况。缩写:CSF = cerebrospinal fluid。请点击此处查看此图的放大版本。

与颈内动脉血流相关参数的分析
采用专用后处理软件中的 4D Flow 模块对四维血流数据进行重建和分析。通过一阶多项式拟合静态组织区域,进行相位偏移校正。使用基于阈值的半自动分割工具并结合手动优化,分离出左侧颈内动脉(LICA),生成三维血管造影容积。系统自动计算血管中心线,并在 C1 段(约位于颈动脉分叉上方 2–3 cm 处)沿中心线垂直方向设置三个横截面平面,相邻平面间距为 5 mm。选择左侧颈内动脉进行分析,因其解剖走行和血流特征较为稳定;未评估右侧颈内动脉,以缩短扫描时间并减轻受试者负担。尽管本项探索性研究采用单侧评估方案以确保可行性,但未来研究中应进行双侧评估。对于每个平面,在每个心动时相(每个心动周期约 20 个时相)均采用半自动边缘检测算法勾画血管管腔,并进行手动调整(图 3A–D)。

根据时间分辨的管腔轮廓,计算了以下参数:最大血流速度、每心动周期的血流量以及平均壁面剪切应力(WSS)。WSS 使用以下公式推导得出:

WSS = μ × (∂v/∂r)|_壁面 (1)

其中 µ = 0.004 Pa·s 为假设的血液黏度,∂v/∂r 为垂直于血管壁的速度梯度,通过在每个周向位置对近壁速度数据拟合抛物线分布,并在管腔周长上取平均获得。生成了血流速度、矢量及壁面剪切应力(WSS)图,并绘制了时间-血流曲线和时间-WSS曲线。三个平面的数值取平均,以获得每位受试者的单一代表性数值。所有图像处理均由两名放射科医师独立完成,统计分析采用其测量结果的均值。图4).

血管血流分析;图A-E:血流速度,MRI切片;血流动力学研究结果。
图3:左侧颈内动脉的三维重建与血流动力学分析。A)对血管进行分割后,左侧颈动脉的三维容积渲染重建图像。(B)沿C1段设置的三个横截面分析平面(Flow 1–3)的位置。(C)血流速度图。(D)显示血流方向的矢量图。(E)管壁剪切应力图。缩写:LICA = left internal carotid artery;WSS = wall shear stress。请点击此处查看该图的放大版本。

流速与轴向壁面剪切应力图;生物系统中流体动力学的比较数据分析。
图 4:左侧颈内动脉后处理曲线。 一个心动周期内的时–流速曲线(左图)以及时–壁面剪切应力曲线(右图)。不同颜色代表三个采样平面(红色、黄色、绿色)。缩写:LICA = 左侧颈内动脉;WSS = 壁面剪切应力。 请点击此处查看该图的放大版本。

统计学分析
数据采用标准统计学软件进行分析(材料表)。通过Shapiro–Wilk检验评估正态分布的变量以均值±标准差(SD)表示,并采用单因素方差分析(one-way ANOVA)进行比较。非正态分布变量以中位数(四分位距,IQR)表示,并采用Kruskal–Wallis检验进行比较。分类变量采用卡方检验进行比较。组间比较采用以年龄和教育程度为协变量的协方差分析(ANCOVA),并对多重比较进行Bonferroni校正。采用Spearman相关分析和多重线性回归分析评估影像学指标与认知测量指标之间的关系,并对血管危险因素(高血压、糖尿病、高脂血症)进行校正。采用受试者工作特征(ROC)曲线分析评估各影像学参数的诊断效能,计算曲线下面积(AUC)、敏感性和特异性。ROC分析未对协变量进行校正,以评估各参数本身的鉴别能力。然而,鉴于观察到的组间在年龄和教育程度上的差异,对这些协变量进行校正的补充分析仍得出一致的结果。以双侧p < 0.05为差异具有统计学意义。本研究中使用的所有材料和设备的完整清单见材料表

结果

临床基线信息和认知评分
该队列包括16名阿尔茨海默病(AD)患者、19名轻度认知障碍(MCI)患者和25名正常对照(NC)参与者。各组在性别、体重、身高、BMI、心血管风险因素和心率方面具有可比性(p > 0.05),但在年龄(p = 0.026)和受教育年限(p = 0.033)上存在差异(表1)。如预期所示,MMSE和MoCA的总分及各子项评分在各组间差异均有统计学意义(均 p < 0.001),表现为AD < MCI < NC(表2)。

脑脊液相关参数的组间差异比较
所有导水管均通畅。CA中的脑脊液流动在心动周期中呈现双相脉动(图2)。在相位图像上,高信号(白色图5A)表示从头侧向尾侧方向的流动(与编码梯度方向相同),而低信号(黑色图5B)表示从尾侧向头侧方向的流动(与编码梯度方向相反)。在校正协变量后,MCI组的净搏动容积高于NC组(p = 0.041)。与NC组相比,MCI组(p = 0.024)和AD组(p = 0.002)的每搏排出量(PDV)较低;AD组的峰值流速(PSV)也低于NC组(p = 0.001)和MCI组(p = 0.025)。各组之间的导水管横截面积差异无统计学意义(表3)。

样本的光学显微镜图像 A 和 B,比例尺为 5 mm,用于对比材料分析。
图 5:脑导水管中脑脊液流动的相差图像。A)相位图像显示脑脊液呈高信号(白色),表示流动方向为头侧至尾侧(与编码梯度方向相同)。(B)相位图像显示脑脊液呈低信号(黑色),表示流动方向为尾侧至头侧(与编码梯度方向相反)。比例尺:5 mm。缩写:CA = 脑导水管;CSF = 脑脊液。 请点击此处查看该图的放大版本。

颈内动脉相关参数的组间差异比较
AD组的左侧颈内动脉(LICA)最大血流速度低于NC组(p = 0.001)和MCI组(p = 0.023)。AD组的平均壁面剪切应力(WSS)低于NC组(p = 0.003)。每心动周期的血流量在各组之间无显著差异(表4)。

脑脊液和颈内动脉相关指标与认知评分的相关性分析结果
颈动脉峰值收缩期流速、峰值舒张期流速、左侧颈内动脉最大流速以及平均壁面剪切应力均与MMSE和MoCA总分呈正相关(均 p < 0.001)(表5)。这些指标与定向、计算、回忆和记忆等子评分的相关性更强。

多变量分析
在校正高血压、糖尿病和高脂血症后,LICA 的最大血流速度与 CA PSV 呈正相关(β = 0.155,p < 0.001),而 LICA 的平均壁面剪切应力与 CA PDV 呈负相关(β = −0.134,p = 0.005)(表 6)。

受试者工作特征分析
受试者工作特征分析显示,颈动脉(CA)收缩期峰值流速、舒张期峰值流速、左侧颈内动脉(LICA)最大血流速度以及平均壁面剪切应力能够区分阿尔茨海默病(AD)与正常对照(NC)(AUC 0.760–0.848,p < 0.05)(表7)(图6)。鉴于AD组与NC组在年龄和教育水平上存在差异,进一步对这些协变量进行校正后重新进行ROC分析,结果相似,提示年龄或教育水平并未对鉴别效能产生显著混杂影响。其中,LICA最大血流速度的AUC最高(0.848),在截断值为36.91 cm/s时,敏感性为72.0%,特异性为93.8%。在区分AD与轻度认知障碍(MCI)方面,CA收缩期峰值流速(PSV)、LICA最大血流速度和平均壁面剪切应力(WSS)显示出中等诊断价值(AUC 0.740–0.743,p < 0.05)(表8)(图7)。

ROC 曲线图;分析 CA、LICA 变量的速度和血流量;数据的敏感性与特异性。
图 6 区分阿尔茨海默病与正常对照的受试者工作特征曲线。 脑导水管(CA)收缩期峰值速度、脑导水管舒张期峰值速度、左侧颈内动脉(LICA)最大血流速度以及 LICA 平均壁面剪切应力的 ROC 曲线。缩写:AD = 阿尔茨海默病;AUC = 曲线下面积;CA = 脑导水管;LICA = 左侧颈内动脉;NC = 正常对照;ROC = 受试者工作特征。 请点击此处查看该图的放大版本。

ROC曲线图,显示临床数据分析中的灵敏度与1-特异度的关系。
图7:用于区分阿尔茨海默病与轻度认知障碍的受试者工作特征曲线。 脑导水管峰值收缩速度、脑导水管峰值舒张速度、左侧颈内动脉最大血流速度及左侧颈内动脉平均壁面剪切力的受试者工作特征曲线。缩写:AD = 阿尔茨海默病;AUC = 曲线下面积;CA = 脑导水管;LICA = 左侧颈内动脉;MCI = 轻度认知障碍;ROC = 受试者工作特征。 请点击此处以查看此图的放大版本。

数据可用性:
支持本研究中统计分析的数值数据集——包括脑脊液流动参数、LICA血流动力学参数、诊断分组和心率——作为补充表1提供。该数据集被直接用于文中报告的所有组间比较、相关性分析、回归模型和ROC分析。研究参与者的认知评估量表见补充表2。为保护参与者隐私,人口统计学和临床数据不公开提供;其汇总统计信息已在表1表2中给出。由于数据量和隐私限制,原始DICOM影像数据不可获取。所有数据均可在合理申请并经机构审查委员会批准后,向通讯作者索取。本研究未生成任何自定义代码;分析使用的是商业软件(IBM SPSS Statistics)。

正常对照组 (n = 25)轻度认知障碍组 (n = 19)阿尔茨海默病组 (n = 16)p
性别 (男/女)13/124/155/110.094
年龄 (年)62.24 ± 4.8260.63 ± 8.1769.44 ± 8.490.026*
体重 (kg)59.00 ± 7.5961.00 ± 7.9755.06 ± 9.580.141
身高 (m)1.62 ± 0.081.59 ± 0.071.60 ± 0.080.413
BMI (kg/m²)22.39 ± 2.5123.69 ± 2.1821.54 ± 3.130.054
高血压 (是/否)5/207/126/100.349
糖尿病 (是/否)2/234/150/160.063
高脂血症 (是/否)6/191/182/140.194
睡眠障碍 (是/否)1/244/154/120.092
受教育年限 (年)5.00 (10.50)5.00 (12.00)0.00 (2.00)0.033*
心率69.52 ± 11.0267.63 ± 9.5970.75 ± 13.260.709

表1: 研究参与者的临床基线特征。 数据以均值±标准差、中位数(四分位距)或数量(百分比)表示。缩写:AD = 阿尔茨海默病;BMI = 体重指数;F = 女性;M = 男性;MCI = 轻度认知障碍;NC = 正常对照。*p < 0.05 表示差异具有统计学显著性。

正常对照组 (n = 25)轻度认知障碍组 (n = 19)阿尔茨海默病组 (n = 16)p
MMSE 评分
总分27.80 ± 1.2923.42 ± 2.3616.50 ± 2.000.000*
定向力10 (0)8 (3)5.5 (1)0.000*
记忆3 (0)3 (0)3 (0)0.525
计算力5 (0)4 (2)1 (3)0.000*
回忆2 (2)0 (2)0 (0)0.000*
语言8 (0)8 (0)6.5 (2)0.000*
结构模仿与视空间能力0 (1)0 (1)0 (0)0.494
MoCA 评分
总分24.80 ± 2.9618.79 ± 3.7410.63 ± 4.180.000*
视空间与执行功能4 (2)2 (2)1 (1)0.000*
命名3 (0)3 (1)2 (2)0.001*
注意力6 (1)5 (1)1.5 (5)0.000*
语言3 (2)2 (2)1 (3)0.067
抽象能力2 (1)1 (2)0 (2)0.001*
记忆与延迟回忆3 (2)0 (3)0 (0)0.000*
定向力6 (0)4 (2)3 (1)0.000*

表2: 三组间认知评分的比较 数据以均值 ± 标准差或中位数(四分位距)表示。缩写:AD = Alzheimer's disease;MCI = mild cognitive impairment;MMSE = Mini-Mental State Examination;MoCA = Montreal Cognitive Assessment;NC = normal control。 p < 0.05 表示具有统计学显著性差异。

NCMCIADNC 与 MCI 比较NC 与 AD 比较MCI 与 AD 比较
每搏总容积 (mL)0.47 ± 0.220.37 ± 0.020.37 ± 0.020.1880.3250.828
净每搏容积 (mL)0.02 ± 0.010.03 ± 0.020.03 ± 0.020.041*0.160.633
收缩期峰值流速 (cm/s)6.72 ± 1.795.96 ± 1.764.67 ± 1.170.1360.001*0.025*
舒张期峰值流速 (cm/s)5.60 ± 1.984.35 ± 1.833.68 ± 1.620.024*0.002*0.253
正向血流流速 (cm/s)0.17 ± 0.070.19 ± 0.110.17 ± 0.790.7160.6390.426
负向血流流速 (cm/s)0.08 ± 0.050.11 ± 0.080.11 ± 0.090.2230.410.774
正向血流量 (mL)0.04 ± 0.020.04 ± 0.020.03 ± 0.020.7930.6550.494
负向血流量 (mL)0.02 ± 0.010.03 ± 0.020.03 ± 0.040.2160.3720.814

表3: 各组间大脑导水管脑脊液流动参数的比较 数据以均值 ± 标准差表示。缩写:AD = Alzheimer's disease;CA = cerebral aqueduct;MCI = mild cognitive impairment;NC = normal control。 p < 0.05 表示具有统计学显著性差异。

正常对照组(NC)轻度认知障碍组(MCI)阿尔茨海默病组(AD)正常对照组 vs 轻度认知障碍组正常对照组 vs 阿尔茨海默病组轻度认知障碍组 vs 阿尔茨海默病组
血流量(mL)2.58 ± 0.762.49 ± 0.782.44 ± 0.740.840.6290.764
最大流速(cm/s)39.73 ± 7.0736.89 ± 7.9531.26 ± 4.250.2140.001*0.023*
平均壁面剪切应力(Pa)0.45 ± 0.160.38 ± 0.150.29 ± 0.830.1720.003*0.063

表4:各组左侧颈内动脉血流动力学参数比较。 数据以均值±标准差表示。缩写:AD = 阿尔茨海默病;LICA = 左侧颈内动脉;MCI = 轻度认知障碍;NC = 正常对照;Pa = 帕斯卡。p < 0.05 表示差异具有统计学意义。

CALICA
净搏出量(ml)收缩期峰值流速(cm/s)舒张期峰值流速(cm/s)最大流速(cm/s)平均壁面剪切应力(Pa)
rprprprprp
MMSE 评分
总分0.1720.190.4730.000*0.4430.000*0.4670.000*0.4360.000*
定向力0.1270.3340.4280.001*0.3850.002*0.4220.001*0.370.004*
记忆力0.1190.3640.0490.7090.1970.1310.0640.6290.1620.216
计算力0.1360.2990.2830.028*0.2610.044*0.4620.000*0.3850.002*
回忆0.1460.2670.4290.001*0.3480.006*0.4320.001*0.3550.005*
语言0.150.2520.3050.018*0.2230.0860.3040.018*0.3310.010*
结构模仿与视空间能力0.0430.7420.0870.5080.0680.6040.250.0540.1240.346
MoCA 评分
总分0.2090.110.4980.000*0.4030.001*0.4390.000*0.4070.001*
视空间与执行功能0.190.1460.3920.002*0.1810.1670.1670.2020.1250.34
命名0.160.2210.2350.0710.1720.1880.380.003*0.3640.004*
注意力0.1470.2640.3210.013*0.1930.140.3050.018*0.390.002*
语言0.1490.2570.2930.023*0.1890.1480.190.1460.1520.247
抽象能力0.3390.008*0.2150.0990.10.4470.110.4030.1580.227
记忆与延迟回忆0.1390.290.4390.000*0.3570.005*0.4960.000*0.4010.001*
定向力0.1360.2980.3860.002*0.380.003*0.3820.003*0.3070.017*

表5:脑脊液和颈内动脉参数与认知评分的相关性分析。 数据以Spearman相关系数(r)及相应的p值表示。缩写:CA = 大脑导水管;LICA = 左侧颈内动脉;MMSE = 简易精神状态检查;MoCA = 蒙特利尔认知评估量表;Pa = 帕斯卡。p < 0.05 表示具有统计学意义的相关性。

颈动脉
收缩期峰值流速(p 值舒张期峰值流速(p 值
左侧颈总动脉最大血流速度0.000*0.255
平均壁面剪切应力0.2230.005*

表6:脑脊液流速与颈内动脉搏动性的多变量回归分析。 数据以标准化回归系数(β)及相应的 p 值表示。缩写:CA = 大脑导水管;LICA = 左侧颈内动脉。p < 0.05 表示具有统计学显著性关联。

AUC灵敏度特异性诊断阈值p
CA
收缩期峰值流速0.8230.640.8755.940.001*
舒张期峰值流速0.760.920.53.2550.005*
LICA
最大流速0.8480.720.93836.910.000*
平均壁面剪切应力0.8260.80.750.3150.000*

表7:用于区分阿尔茨海默病与正常对照的受试者工作特征分析。AUC = 曲线下面积;AD = 阿尔茨海默病;CA = 脑导水管;LICA = 左颈内动脉;NC = 正常对照。p < 0.05 表示诊断效能具有统计学显著性。

AUC敏感性特异性诊断阈值P
CA
收缩期峰值流速0.740.7890.6884.90.016*
舒张期峰值流速0.6070.4210.8134.8750.282
LICA
最大流速0.7430.5790.93836.650.014*
平均壁面剪切应力0.740.4740.9380.3970.016*

表8:用于区分阿尔茨海默病与轻度认知障碍的受试者工作特征分析。AUC = 曲线下面积;AD = 阿尔茨海默病;CA = 脑导水管;LICA = 左侧颈内动脉;MCI = 轻度认知障碍。p < 0.05 表示诊断效能具有统计学显著性。

补充表 1:统计分析的完整数值数据集。 本表格包含用于所有统计分析的去标识化数值数据,包括组间比较、相关性分析、回归模型以及受试者工作特征(ROC)分析。请点击此处下载该文件。

补充表 2:研究参与者的去标识化认知评分。 本表包含所有参与者的去标识化认知评估数据,包括简易精神状态检查(MMSE)和蒙特利尔认知评估(MoCA)的总分及各子项得分。参与者标识信息已被移除,并替换为匿名的受试者编号,以保护患者隐私。缩写:MMSE = 简易精神状态检查;MoCA = 蒙特利尔认知评估。请点击此处下载该文件。

讨论

通过整合二维电影相位对比磁共振成像(2D cine PC-MRI)和四维流磁共振成像(4D flow MRI),本研究对阿尔茨海默病(AD)谱系中脑脊液动力学及颈段大血管血流动力学的改变及其相互关系进行了定量评估。主要研究结果总结如下:(1)AD和轻度认知障碍(MCI)患者的中脑导水管内脑脊液收缩期峰值和舒张期峰值流速均显著降低,且这些降低程度与认知评分呈相关性。(2)AD组患者左侧颈内动脉(LICA)的最大流速和壁面剪切应力(WSS)显著降低。(3)LICA最大流速与导水管脑脊液收缩期峰值流速呈正相关,而LICA平均WSS与脑脊液舒张期峰值流速呈负相关。(4)这些血流动力学参数,尤其是LICA最大流速,在区分AD患者与认知正常(NC)个体及MCI个体方面表现出良好的诊断效能。上述结果提示,在AD进展过程中,脑脊液-血管关联单元可能发生功能紊乱,且颈段大血管血流动力学可能作为颅内类淋巴功能的间接影像学指标。然而,该解释仍属推测性,需在未来研究中通过直接证据加以验证。

在导水管中观察到的双向、准正弦样26,27,28脑脊液(CSF)流动模式,与早期二维相位对比磁共振成像(2D PC-MRI)研究建立的经典描述一致29,21。在一个心动周期中,脑脊液在中央动脉(CA)内的最大收缩期流速始终高于最大舒张期流速,导致每搏净流量持续朝向头-足方向,这与既往研究结果相符30。然而,阿尔茨海默病(AD)和轻度认知障碍(MCI)患者脑脊液峰值流速的普遍降低,提示驱动这一生理振荡的主要动力受损。根据Monro-Kellie学说,在固定的颅内体积内,动脉血的脉动性流入是产生脑脊液往复运动的主要驱动力31,32。因此,本研究中颈动脉观察到的动脉脉动性减弱(表现为最大流速和壁面剪切应力WSS降低)很可能直接导致脑脊液振荡减弱。这一点得到了多元回归分析结果的有力支持,结果显示左侧颈总动脉(LICA)最大流速与导水管脑脊液峰值收缩期流速呈正相关。值得注意的是,尽管MCI组的脑脊液峰值流速呈下降趋势,但其每搏净流量却与正常对照(NC)组相比出现反常性升高。这可能提示一种 早期代偿机制:当神经血管调节功能开始受损时,脑灌注可能通过增加每搏流量来维持总体脑脊液流量,类似于脑血管储备研究中观察到的代偿性高灌注现象。随着AD病情进展,这种代偿最终失效,表现为AD组中脑脊液动力学全面下降。

本研究发现,阿尔茨海默病(AD)患者颈内动脉的血流速度和壁面剪切应力(WSS)显著降低。WSS 是指流动血液对血管壁施加的切向摩擦力,是维持内皮细胞健康的关键生物力学信号。WSS 降低与内皮功能障碍、动脉粥样硬化及血管重塑密切相关。在 AD 的背景下,颈动脉 WSS 减弱可能通过两条主要途径加剧疾病进展: 首先,它直接导致脑灌注压和脑血流量下降,这与 AD 中已有充分记载的脑低灌注现象一致。 大量研究表明,脑白质病变33,34,35和脑血流量减少36,37是 AD 的重要病理特征。其次,异常的血流动力学信号可能通过损害内皮功能,加重血脑屏障(BBB)破坏和血管炎症,从而促进 Aβ 沉积和 tau 蛋白病理扩散38,39。这进一步引发微循环障碍,促进 AD 的发生与发展40。因此,颈动脉大血管血流动力学的改变不仅是脑灌注减少的 原因 ,也可能是 AD 相关血管病理的 结果 和 放大器 。颈动脉血流动力学参数与导水管脑脊液(CSF)流速之间观察到显著相关性。越来越多的证据表明,动脉搏动是类淋巴系统(glymphatic system)——即沿血管周围间隙促进脑脊液与组织间液交换的脑部清除通路——的主要驱动力。观察到较高的颈动脉流速/壁面剪切应力与较高的收缩期 CSF 流速相关,可解释为较强的动脉搏动通过血管周围通路有效地将脑脊液驱动至脑实质内。近期研究进一步证实了颈动脉和颅内动脉搏动指数与脑淋巴系统之间的关联,后者在 AD 的发病机制和疾病进展中起关键作用41。因此,AD 患者因颈动脉搏动减弱所导致的脑脊液振荡减弱,很可能对应着类淋巴清除功能受损,从而造成 Aβ 等代谢废物清除效率下降。这一关联背后的机制可能涉及多种因素:AD 和轻度认知障碍(MCI)患者脑脊液中 Aβ、总 tau 蛋白和磷酸化 tau 蛋白水平的改变,以及局部脑低灌注和白质损伤,这些因素可能间接导致脑脊液流速异常42。该解释与 AD 动物模型中的观察结果一致,在这些模型中,类淋巴功能障碍被发现早于斑块沉积。然而需要指出的是,这些解释仍是间接的,因为本研究并未直接测量类淋巴功能,其他可能的解释——如动脉顺应性降低、局部脑萎缩或颅内顺应性改变——尚不能排除。

基于相位对比磁共振成像(PC-MRI)获得的流体动力学参数与认知评分(MMSE/MoCA)显著相关,尤其在定向力、计算力和延迟回忆等核心认知领域表现明显。这种关联可能反映了这些参数对与阿尔茨海默病(AD)相关认知功能障碍密切相关的神经病理改变具有敏感性,包括嗅皮层和颞叶皮层萎缩39、内侧颞叶和内侧顶叶皮层的tau蛋白沉积,以及前脑、海马、基底节和颞顶皮层区域脑灌注减少43。研究结果提示,脑脊液和颈动脉的流体力学状态可能与AD相关的神经功能障碍有关。受试者工作特征(ROC)分析显示,左侧颈总动脉(LICA)最大血流速度(截断值≈36.9 cm/s)在识别AD时具有高特异性(>93%)。该发现提示其潜在的转化价值:作为一种相对快速且标准化的MRI测量指标,颈动脉血流速度或可作为早期AD筛查或其他类型痴呆鉴别诊断的实用辅助影像学指标。尽管单独使用时敏感性可能中等,但将该指标与结构MRI、Aβ-PET或认知量表联合应用,有望增强多模态诊断模型的效能。然而,上述解释基于间接测量结果,仍需在更大规模、独立的队列中加以验证。

本研究存在若干局限性,在解读研究结果时应予以考虑。首先,横断面研究设计无法进行因果推断。其次,样本量较小,限制了统计效能和结果的普适性。第三,临床诊断未通过生物标志物确认(如脑脊液 Aβ42/p-tau 或淀粉样蛋白 PET),可能引入诊断不确定性。第四,仅分析了左侧颈内动脉。这种单侧评估方法可能无法捕捉脑血流动力学潜在的半球间差异,研究结果可能无法完全推广至右侧颈动脉供血区域。未来需要开展双侧评估的研究,以明确所观察到的关联是否具有对称性或侧别依赖性。第五,4D 流量 MRI 采用的编码速度(Venc)为 50 cm/s,尽管该参数选择旨在兼顾颈动脉血流捕捉与避免混叠伪影,但对于血流速度异常高的受试者,可能限制了速度测量的准确性。第六,缺乏 T1 加权的脑萎缩测量指标,因此无法判断血流速度降低是由实质萎缩引起,还是源于原发性血流动力学改变。第七,将所观察到的变化与类淋巴功能障碍相关联的解释仍是间接的,因本研究未进行直接的类淋巴系统成像或脑脊液清除率测量。第八,尽管 ROC 分析提示了一定的诊断价值,但这些结果来自单一队列,且主要分析中未对年龄和教育程度进行校正;尽管补充的校正分析结果一致,仍需在更大规模、独立的队列中加以验证。最后,ROC 分析结果需在独立队列中进一步验证。未来研究应采用纵向设计,纳入双侧颈动脉评估和结构成像,并整合多组学数据,以验证这些初步发现。

综上所述,本研究结合二维电影相位对比磁共振成像(2D cine PC-MRI)和四维流磁共振成像(4D flow MRI),揭示了在整个阿尔茨海默病(AD)谱系中,颅内脑脊液振荡与颈段大血管血流动力学的耦合衰减及其相互关系。这些改变与认知功能障碍的程度相关,并具有中等程度的诊断价值。研究结果支持血管因素在AD发病机制中的作用,并将血流动力学异常的范围扩展至颈动脉入颅部位。本研究提出一种假设:颈动脉搏动性血流可能影响颅内类淋巴功能,其搏动减弱可能作为类淋巴功能障碍的间接影像学标志。然而,这一解释仍属推测性,需在未来研究中通过直接证据加以验证。这些发现为AD的发病机制提供了流体动力学视角,并可能为中枢神经系统液体与屏障生理及病理相关的无创性动态监测与评估工具的开发提供依据。

披露

作者声明无相关财务或非财务利益需要披露。

致谢

感谢所有参与本研究的受试者及法定监护人,以及宁德市医院放射科技术人员提供的技术支持。本研究由宁德师范学院医学院融合发展专项研究项目[项目编号 2024ZX70]和中国福建省卫生技术项目[项目编号 2022CXA060]资助。

材料

本文使用的材料清单
姓名公司目录编号评论
28通道头颈线圈GE Healthcare, Waukesha, WI, USAGE SIGNA Architect系统组成部分;https://www.gehealthcare.com专用于头颈部联合成像的接收专用相控阵线圈;用于所有MRI序列中的信号接收。
2D电影相位对比MRI序列(Fast PC Cine Slice)GE Healthcare, Waukesha, WI, USAGE SIGNA Architect软件包组成部分二维回顾性门控相位对比序列;用于量化中脑导水管内脑脊液流速和每搏输出量。参数:Venc 12 cm/s,每个心动周期采集20–30个时相,扫描平面垂直于导水管。
4D血流MRI序列(Sag 4D Flow)GE Healthcare, Waukesha, WI, USAGE SIGNA Architect软件包组成部分时间分辨的三维相位对比序列;用于获取颈段颈动脉内空间分辨的血流动力学数据。参数:Venc 50 cm/s,每个心动周期20–30个时相,40–50层。
CVI42后处理工作站(版本5.14)Circle Cardiovascular Imaging, Calgary, Canadahttps://www.circlecvi.com/cvi42/专用的心血管及血流分析软件;用于血管分割、感兴趣区域勾画、流量定量(脑脊液和血液)以及管壁剪切应力计算。版本5.14。
GE SIGNA Architect 3.0T MRI扫描仪GE Healthcare, Waukesha, WI, USAhttps://www.gehealthcare.com/products/magnetic-resonance-imaging/signa-architect配备28通道头颈正交线圈的全身3.0T超导MRI系统;用于所有结构成像、2D电影相位对比及4D血流数据采集。
IBM SPSS Statistics(版本25)IBM Corp., Armonk, NY, USASCR_002865; https://www.ibm.com/products/spss-statistics统计分析软件;用于所有描述性统计、组间比较、相关性分析、回归分析及ROC曲线分析。版本25。

参考文献

  1. Graff-Radford J, et al. New insights into atypical Alzheimer's disease in the era of biomarkers. The Lancet Neurology. 2021;20(3):222-34.
  2. Zhang XW, Zhu XX, Tang DS, Lu JH. Targeting autophagy in Alzheimer’s disease: Animal models and mechanisms. Zoological Research. 2023;44(6):1132-45.
  3. Muroni A, et al. Point prevalence of epilepsy in dementia: A "real-world" estimate. Epileptic Disorders. 2024.
  4. Shokri-Kojori E, et al. β-Amyloid accumulation in the human brain after one night of sleep deprivation. Proceedings of the National Academy of Sciences. 2018;115(17):4483-8.
  5. Chandra A., et al. Aquaporin-4 polymorphisms predict amyloid burden and clinical outcome in the Alzheimer's disease spectrum. Neurobiology of Aging. 2021;97:1-9.
  6. Iliff JJ, et al. A Paravascular Pathway Facilitates CSF Flow Through the Brain Parenchyma and the Clearance of Interstitial Solutes, Including Amyloid β. Science Translational Medicine. 2012;4(147).
  7. Wu C, et al. Continuous Theta-Burst Stimulation Promotes Paravascular CSF-Interstitial Fluid Exchange through Regulation of Aquaporin-4 Polarization in APP/PS1 Mice. Mediators of Inflammation. 2022;2022:1-12.
  8. Buchhave P. Cerebrospinal Fluid Levels of β-Amyloid 1-42, but Not of Tau, Are Fully Changed Already 5 to 10 Years Before the Onset of Alzheimer Dementia. Archives of General Psychiatry. 2012;69(1).
  9. Peng W, et al. Suppression of glymphatic fluid transport in a mouse model of Alzheimer's disease. Neurobiology of Disease. 2016;93:215-25.
  10. Ben-Nejma IRH, et al. Altered dynamics of glymphatic flow in a mature-onset Tet-off APP mouse model of amyloidosis. Alzheimer's Research & Therapy. 2023;15(1):23.
  11. Algin O, Hakyemez B, Parlak M. The Efficiency of PC-MRI in Diagnosis of Normal Pressure Hydrocephalus and Prediction of Shunt Response. Academic Radiology. 2010;17(2):181-7.
  12. Cavallin L, et al. Voxel-based correlation between coregistered single-photon emission computed tomography and dynamic susceptibility contrast magnetic resonance imaging in subjects with suspected Alzheimer's disease. Acta Radiologica. 2008;49(10):1154-61.
  13. Eskildsen SF, et al. Increased cortical capillary transit time heterogeneity in Alzheimer's disease: a DSC-MRI perfusion study. Neurobiology of Aging. 2017;50:107-18.
  14. Rivera-Rivera LA, et al. Changes in intracranial venous blood flow and pulsatility in Alzheimer’s disease: A 4D flow MRI study. Journal of Cerebral Blood Flow & Metabolism. 2016;37(6):2149-58.
  15. Korbecki A, et al. Imaging of cerebrospinal fluid flow: fundamentals, techniques, and clinical applications of phase-contrast magnetic resonance imaging. Polish Journal of Radiology. 2019;84:240-50.
  16. Vikner T, et al. CSF dynamics throughout the ventricular system using 4D flow MRI: associations to arterial pulsatility, ventricular volumes, and age. Fluids and Barriers of the CNS. 2024;21(1).
  17. Duck CJ, et al. Choroid Plexus Volume and Permeability at Brain MRI within the Alzheimer's Disease Clinical Spectrum. Radiology. 2022;304(3):212400.
  18. Tucker O, Samuel M, Bruce G, Scott F. Noninvasive Calculation of Cerebrospinal Fluid Production Using Phase-Contrast Magnetic Resonance Imaging: First Implementation in Augusta, Georgia. Cureus. 2023;15(5):e39686.
  19. Zhuang B, Sirajuddin A, Zhao S, Lu M. The role of 4D flow MRI for clinical applications in cardiovascular disease: current status and future perspectives. Quantitative Imaging in Medicine and Surgery. 2021;11(9):4193-210.
  20. Bradley WG, et al. Normal-pressure hydrocephalus: evaluation with cerebrospinal fluid flow measurements at MR imaging. Radiology. 1996;198(2):523-9.
  21. Bateman GA. Pulse-wave encephalopathy: a comparative study of the hydrodynamics of leukoaraiosis and normal-pressure hydrocephalus. Neuroradiology. 2002;44(9):740-8.
  22. El Sankari S, et al. Cerebrospinal fluid and blood flow in mild cognitive impairment and Alzheimer’s disease: a differential diagnosis from idiopathic normal pressure hydrocephalus. Fluids and Barriers of the CNS. 2011;8(1):1-12.
  23. Pahlavian SH, et al. Cerebroarterial pulsatility and resistivity indices are associated with cognitive impairment and white matter hyperintensity in elderly subjects: A phase-contrast MRI study. Journal of Cerebral Blood Flow & Metabolism. 2020;41(3):670-83.
  24. Rivera-Rivera LA, et al. Assessment of vascular stiffness in the internal carotid artery proximal to the carotid canal in Alzheimer’s disease using pulse wave velocity from low rank reconstructed 4D flow MRI. Journal of Cerebral Blood Flow & Metabolism. 2020;41(2):298-311.
  25. Karki P, et al. Real-time phase-contrast MRI to monitor cervical blood flow and CSF flow: a feasibility study. Magnetic Resonance in Medicine. 2022;87(2):966-77.
  26. Lee J, et al. The impact of breathing patterns on CSF flow and global brain BOLD signal while awake: A study using real-time phase-contrast MRI and fast rsfMRI. Alzheimer's & Dementia. 2023;19(S19).
  27. Lotz J, et al. In vitro validation of phase-contrast flow measurements at 3 T in comparison to 1.5 T: Precision, accuracy, and signal-to-noise ratios. Journal of Magnetic Resonance Imaging. 2005;21(5):604-10.
  28. Pan L, et al. Validating the accuracy of real-time phase-contrast MRI and quantifying the effects of free breathing on cerebrospinal fluid dynamics. Fluids and Barriers of the CNS. 2024;21(1):25.
  29. Greitz D, et al. Pulsatile brain movement and associated hydrodynamics studied by magnetic resonance phase imaging. Neuroradiology. 1993;35(5):350-5.
  30. Virhammar J, Laurell K, Cesarini KG, Larsson EM. Aqueductal CSF Stroke Volume Is Increased in Patients with Idiopathic Normal Pressure Hydrocephalus and Decreases after Shunt Surgery. AJNR Am J Neuroradiol. 2019;40(3):453-9.
  31. Balédent O, et al. Relationship Between Cerebrospinal Fluid and Blood Dynamics in Healthy Volunteers and Patients with Communicating Hydrocephalus. Investigative Radiology. 2004;39(1):45-55.
  32. Jennifer H, Asha S, Stephen L, Sumit P. Systematic Approach to Pediatric Macrocephaly. Radiographics. 2023;43(5):e220159.
  33. Takase H, et al. Transcriptomic Profiling Reveals Neuroinflammation in the Corpus Callosum of a Transgenic Mouse Model of Alzheimer's Disease. Journal of Alzheimer's Disease. 2024.
  34. Yamamoto R, et al. Effects of hippocampal atrophy and white matter lesions on cognitive function in Alzheimer’s disease. Alzheimer's & Dementia. 2023;19(S17).
  35. Yaqing L, Jiaxin Z, Tian L, Junjian Z. White Matter and Alzheimer's Disease: A Bidirectional Mendelian Randomization Study. Neurology and Therapy. 2022;11(2):881-92.
  36. Fang X, et al. Amyloid beta accumulation in TgF344-AD rats is associated with reduced cerebral capillary endothelial Kir2.1 expression and neurovascular uncoupling. GeroScience. 2023;45(5):2909-26.
  37. Pyne J, et al. Total and posterior cerebral blood flow is associated with white matter hyperintensity volume and Alzheimer ’s-like patterns of brain atrophy in older adults. Alzheimer's & Dementia. 2023;19(S17).
  38. Axel M, et al. Undetectable gadolinium brain retention in individuals with an age-dependent blood-brain barrier breakdown in the hippocampus and mild cognitive impairment. Alzheimer's & Dementia. 2019;15(12):1568-75.
  39. Sofia M, et al. Effects of Interventions on Cerebral Perfusion in the Alzheimer's Disease Spectrum: A Systematic Review. Aging Research Reviews. 2022;79:101661.
  40. Costantino I. The overlap between neurodegenerative and vascular factors in the pathogenesis of dementia. Acta Neuropathologica. 2010;120(3):287-96.
  41. Xie L, et al. Higher intracranial arterial pulsatility is associated with presumed imaging markers of the glymphatic system: An explorative study. NeuroImage. 2024;288.
  42. Erdélyi-Bótor S, et al. Serum L-arginine and dimethylarginine levels in migraine patients with brain white matter lesions. Cephalalgia. 2016;37(6):571-80.
  43. Odano I, et al. Diagnostic approach with Z-score mapping to reduce artifacts caused by cerebral atrophy in regional CBF assessment of mild cognitive impairment (MCI) and Alzheimer's disease by [99mTc]-ECD and SPECT. Japanese Journal of Radiology. 2024.

重印与许可

标签