基于多模态磁共振成像(Multimodal Magnetic Resonance Imaging)聚类的初发帕金森病两种影像定义脑表型鉴定

《NeuroImage: Clinical》：Multimodal magnetic resonance imaging-based clustering identifies two imaging-defined brain phenotypes in De novo Parkinson's disease

　　

摘要：帕金森病(Parkinson's disease, PD)存在显著的临床与神经生物学异质性，为其预后判断与治疗研发带来困难。研究人员旨在整合弥散磁共振成像(dMRI)与结构磁共振成像(sMRI)特征，采用数据驱动聚类方法鉴定影像定义的PD亚型。研究获取帕金森病进展标志物倡议(Parkinson's Progression Markers Initiative, PPMI)数据库中103例初发PD受试者与40例正常对照(Normal Control, NC)的基线神经影像学数据，提取209项影像特征，包括弥散MRI得出的各向异性分数(Fractional Anisotropy, FA)与平均扩散率(Mean Diffusivity, MD)，以及sMRI得出的皮层与皮层下体积。研究人员采用基于正常对照的、考虑交互作用的回归模型对影像特征进行人口学、解剖学与扫描仪相关混杂因素校正。对校正后特征的主成分(Principal Component, PC)实施Ward联接层次聚类(Hierarchical Clustering with Ward's Linkage)。最优聚类数由轮廓系数(Silhouette Index)与Calinski–Harabasz指数确定，聚类稳定性通过跨算法验证与自举法(Bootstrapping)评估。鉴定出两种影像定义的PD亚型：A亚型表现为多白质纤维束FA广泛升高、MD广泛降低，无显著体积萎缩，符号数字模态测验(Symbol Digit Modalities Test, SDMT)得分较高，发病年龄较早；B亚型表现为MD广泛升高、区域特异性FA降低、接近显著的皮层体积缩小、脑室扩大及横断面认知表现较低。48个月随访的线性混合效应模型(Linear Mixed-Effects Model, LME)显示认知轨迹分离：A亚型认知评分稳定，B亚型出现显著认知衰退。这些发现鉴定出早期初发PD中存在两种具不同认知轨迹的影像定义结构表型，为理解疾病最早阶段PD异质性提供了多模态影像框架。

论文解读：《Multimodal magnetic resonance imaging-based clustering identifies two imaging-defined brain phenotypes in De novo Parkinson's disease》发表于《NeuroImage: Clinical》

一、研究背景与立项依据

帕金森病(Parkinson's disease, PD)虽共享相同临床诊断，但其临床进程、症状谱、进展速率及治疗反应差异显著，传统"一刀切"的疾病模型面临挑战。既往PD亚型划分多由假设驱动转向多维临床与生物数据驱动，其中神经影像学是阐明PD异质性神经生物学基础的关键工具。结构磁共振成像(structural Magnetic Resonance Imaging, sMRI)可评估宏观结构脑退化，已有研究利用皮层厚度揭示与不同神经心理、运动特征及发病年龄相关的PD模式；但宏观结构改变可能是神经退行较晚期的表现，未必能捕捉早期微结构病理。弥散MRI(Diffusion MRI, dMRI/DTI)通过各向异性分数(Fractional Anisotropy, FA)与平均扩散率(Mean Diffusivity, MD)敏感探测白质微结构完整性，可在体积萎缩前检测到轴索退变与神经炎症改变。然而单一DTI用于显式患者分型仍不成熟，近期研究倾向将DTI与sMRI互补整合以跨宏观—微尺度全面表征PD相关神经退变。目前仍缺乏在无监督数据驱动聚类框架内专门整合sMRI与DTI的研究，且聚类前常未充分校正年龄、性别等混杂因素，导致亚型可能反映人口学而非真实生物学变异。本研究以初发(de novo)PD为对象，采用多模态神经影像无监督聚类并结合基于正常对照的、含交互项混杂校正框架，鉴定影像定义的PD亚型并探讨其临床与纵向预后意义。

二、主要关键技术方法概述

研究纳入帕金森病进展标志物倡议(Parkinson's Progression Markers Initiative, PPMI)数据库的103例初发PD患者与40例正常对照(NC)。采集3.0 T Siemens扫描仪的三维T₁加权结构像(TR=2300 ms, TE=3.0 ms, TI=900 ms, 1 mm各向同性体素)与单次激发平面回波弥散加权像(b=1000 s/mm2, 64方向, 2 mm各向同性体素)。sMRI经CAT12/SPM12预处理并行Harvard-Oxford图谱分割获取48个皮层区与21个皮层下区体积；DTI经MRtrix3去噪、FSL涡流校正与脑提取、DTIFIT拟合FA/MD图，非线性配准至FMRIB58_FA模板，应用JHU-ICBM标签(50区)与纤维束(20束)图谱提取70个白质区FA与MD均值。共209个影像特征。基于NC拟合含年龄、性别、受教育年限、总颅内体积(Total Intracranial Volume, TIV——DTI特征不含TIV)、扫描仪型号及组别与年龄×组别交互项的线性回归，取残差作为校正后特征。RobustScaler标准化后按方差阈值筛选保留177特征，PCA降维保留贡献>1%的主成分(18个,累积方差78.2%)，在此空间进行Ward联接层次聚类，以Silhouette与Calinski–Harabasz指数确定最佳聚类数k=2，并通过k-means交叉验证与自举法评估稳定性。组间人口学与临床差异用t检验/卡方检验与广义线性模型(General Linear Model, GLM)并False Discovery Rate(FDR)校正；亚型vs NC用GLM检验；纵向认知轨迹用以MoCA为因变量、含亚型×时间交互项的线性混合效应模型(Linear Mixed-Effects Model, LME)分析。

三、研究结果

3.1. Demographic and clinical characteristics(人口统计学与临床特征)

层次聚类确定两亚型：A亚型59例，B亚型44例。两亚型性别与受教育年限无差异；Subtype A发病年龄更早(58.7±5.95岁 vs 65.73±6.45岁, p<0.001)，基线年龄更低(60.73±6.38岁 vs 67.73±6.08岁, p<0.001)。运动症状(UPDRS III、震颤、强直、运动迟缓分项)及非运动量表(SCOPA-AUT、GDS、RBDQ等)在两亚型间无显著差异。Subtype A的SDMT(Symbol Digit Modalities Test)得分显著高于Subtype B(42.41±8.3 vs 37.41±9.46, adj. p=0.02)，提示A亚型信息处理速度相对保留。

3.2. Biological characterization(生物学特征刻画)

亚型间比较显示Subtype B比Subtype A多个皮层(缘上回、额岛盖、中央岛盖等)与皮层下区(杏仁核、丘脑等)体积减小(Cohen's d>0.8或>0.5)，伴左侧脑室扩大；Subtype B在多白质纤维束(胼胝体各部、额枕束、丘脑前辐射、扣带束、皮质脊髓束等)FA显著降低、MD显著升高。与NC比较，Subtype A无显著体积差异，但全脑多联络/连合/投射/脑干纤维束FA广泛升高(FDR p<0.001–0.021)、部分纤维束MD降低(FDR p≤0.023)；Subtype B与NC比未见达显著性体积萎缩(边缘性减小p≈0.065)，FA降低范围较局限(胼胝体膝/体、右丘脑前辐射等, FDR p=0.039–0.046)，而MD广泛升高涉及多联络/投射/连合纤维及中小脑脚(all FDR p=0.004–0.049)。

3.3. Correlations between biological and clinical parameters(生物学与临床参数相关性)

控制协变量后全PD样本中双侧杏仁核与颞中回体积及右侧上放射冠/上纵束MD与发病年龄显著相关(FDR校正)。按亚型分层未达FDR显著性，但Subtype A内呈中度相关趋势(r=0.30–0.42, raw p<0.01)，Subtype B内无此关联。SDMT与影像特征无FDR显著相关。

3.4. Validation of PD clusters(PD聚类验证)

100次自举重采样后，各亚型特征z值剖面与原聚类特征剖面高度相关(r=0.99)，表明生物学签名可重复。

3.5. Longitudinal clinical relevance(纵向临床相关性)

48个月LME模型显示Subtype A的MoCA变化率极小(β₁=?0.03分/年, p=0.70)，认知轨迹稳定；Subtype B显著下降(β₂=?0.47分/年, p<0.001)。亚型×时间交互显著(β₃=?0.43分/年, 95%CI ?0.69~?0.18, p<0.001)，即B亚型比A亚型每年多衰退约0.43 MoCA分，4年累计差异约1.7分，且该差异不因替换发病年龄为协变量而改变。

四、讨论总结与结论翻译

讨论指出两亚型非单纯严重程度梯度，而是具质异性的宏观—微结构组合模式：Subtype A的白质FA广泛升高伴MD降低及体积保留，可能为早期多巴胺能缺失后的代偿性微结构重组或交叉纤维选择性丢失所致(需NODDI等进一步验证)，对应较好横断面认知与稳定4年轨迹；Subtype B呈弥漫MD升高(部分伴FA降低)及趋势性皮层萎缩、脑室扩大，提示弥漫微结构改变表型，与较快认知衰退相符。该方法学上采用含年龄×组别交互项的规范化回归去除疾病修饰的正常老化效应，纳入扫描仪协变量，并经交叉算法与自举验证确保聚类稳健。局限性含样本量适中、DTI受交叉纤维影响、特征筛选与聚类距离假设、初发队列限制外推及未测遗传等生物因素。

结论(翻译)：

本研究通过多模态无监督聚类框架鉴定出初发PD的两种影像定义亚型，表明早期PD异质性不能用单一严重轴描述，而反映宏观与微结构组织质的差异组合。Subtype A以白质FA广泛升高、皮层与皮层下体积保留、较轻横断面认知损害及48个月随访中稳定认知轨迹为特征；Subtype B以MD弥漫升高、局限区域性FA降低、趋势性体积缩小、脑室扩大及显著更快认知衰退为特征。结果表明多模态特征融合对生物学有意义的PD分型不可或缺，宏观测量单独不足以捕获早期疾病异质性，关键病理信号主要为微结构性；FA与MD解离强调弥散指标互补价值，MD是早期PD弥散组织破坏的敏感标记。基于正常对照的交互感知混杂校正框架、扫描仪方差吸收及自举验证共同支持亚型有效性。Subtype A中FA升高与其保留认知轨迹相符，支持代偿或前退变解读，但最终需补充成像对比确认。这些发现为理解最早阶段PD异质性提供生物学框架，并为预后分层奠基。