《Journal of the Neurological Sciences》：Data-driven progression trajectories across clinical subtypes of amyotrophic lateral sclerosis
编辑推荐：
肌萎缩侧索硬化症（amyotrophic lateral sclerosis, ALS）是一种进行性且致命的神经退行性疾病，以上、下运动神经元退行性变为特征。临床上，ALS表现为异质性的运动功能障碍模式，包括上肢为主型、下肢为主型和延髓为主型，各亚型与不同的疾
肌萎缩侧索硬化症（amyotrophic lateral sclerosis, ALS）是一种进行性且致命的神经退行性疾病，以上、下运动神经元退行性变为特征。临床上，ALS表现为异质性的运动功能障碍模式，包括上肢为主型、下肢为主型和延髓为主型，各亚型与不同的疾病轨迹和预后相关。然而，传统的分类方法主要基于横断面数据，未能充分捕捉疾病进展的时间动态特征。为弥补这一不足，研究人员采用亚型与分期推断（Subtype and Stage Inference, SuStaIn）这一概率机器学习框架，对来自大型ALS队列的ALSFRS-R（ALS功能评定量表修订版）条目水平数据进行分析，重建了各数据驱动亚型随病程分期的功能衰退轨迹，并推断个体所处的疾病分期。研究识别出三种具有临床可解释性的功能受累模式：上肢为主型（S1）、下肢为主型（S2）和延髓为主型（S3）。在180天和365天随访队列中，延髓为主型亚型的功能衰退速度显著快于上肢和下肢为主型亚型（Kruskal–Wallis检验，p < 0.001），该差异在调整年龄和性别后仍然显著。纵向稳定性分析显示，180天队列中63.0%的参与者保持亚型不变（Cohen's κ = 0.42），365天队列中为53.8%（Cohen's κ = 0.26）；分期进展的单调性（稳定或进展）在180天和365天队列中分别达85.0%和90.8%。生存分析表明，SuStaIn分期（每增加10期）与死亡风险显著相关，而亚型本身在整体队列中未显示显著的独立预后价值；但在低分期亚组中，亚型间存在显著生存差异。该研究证实，基于ALSFRS-R条目水平的SuStaIn建模能够重建具有临床意义的ALS功能进展轨迹，为ALS患者分层提供了一种可扩展的框架。
论文解读文章
一、研究背景与目的
肌萎缩侧索硬化症（amyotrophic lateral sclerosis, ALS）是一种进行性且致命的神经退行性疾病，以上、下运动神经元同时退行性变为核心病理特征。临床上，ALS呈现高度异质性的运动功能障碍模式，传统上依据主要受累部位分为上肢为主型、下肢为主型和延髓为主型三大类。这种基于发病部位的三分法框架在临床实践中被广泛使用，不同亚型对应不同的疾病轨迹和预后——例如延髓为主型通常较肢体为主型预后更差。然而，这一传统分类本质上属于横断面描述，仅反映某一时间点的功能受累格局，无法充分刻画疾病进展的时间动态特征。ALS的病程演变具有明显的时序性——从局灶性起病逐步向邻近解剖区域扩散，最终累及多系统运动功能——因此，仅依赖静态分类难以全面反映疾病进展的全貌。
为突破这一局限，需要能够同时建模表型异质性和疾病进展动态的数据驱动方法。亚型与分期推断（Subtype and Stage Inference, SuStaIn）是一种概率机器学习框架，可利用横断面数据同时识别疾病亚型并重建各亚型的时序进展模式。该方法已在多种神经退行性疾病中成功应用，揭示了具有生物学意义的亚型和分期依赖性轨迹。与此同时，ALS功能评定量表修订版（ALSFRS-R）作为广泛采用的临床量表，无需特殊设备即可完成亚型和分期分配，为这一框架在不同临床场景中的规模化应用提供了可行性基础。
基于上述背景，本研究旨在将SuStaIn应用于大型ALS队列的ALSFRS-R条目水平数据，重建各数据驱动亚型随病程分期的功能衰退轨迹，并检验这些亚型是否对应临床公认的功能受累模式，进而评估其在纵向稳定性、功能衰退速度和生存预后方面的临床相关性。研究人员假设，SuStaIn能够识别出具有临床可解释性的ALSFRS-R条目水平功能受累模式，并重建与纵向功能衰退和生存相关的分期依赖性进展轨迹。
二、主要技术方法
本研究使用来自PRO-ACT（Pooled Resource Open-Access ALS Clinical Trials）公共数据库的纵向临床数据，该数据库汇集了多项ALS临床试验的去标识化数据。研究人员从2025年11月12日下载的基线记录中筛选出3067名符合条件的ALS患者纳入分析。核心分析采用SuStaIn算法，该算法是一种概率机器学习框架，能够同时识别疾病亚型并重建各亚型的时序进展轨迹。由于ALSFRS-R条目为有序分类变量（0–4分），分析采用pySuStaIn中的有序模型，以12个ALSFRS-R条目（Q1–Q12）作为输入特征，并在预处理阶段对年龄和性别进行校正。模型通过200,000次马尔可夫链蒙特卡洛迭代和100次随机初始化进行评估，基于贝叶斯信息准则进行模型选择，并通过最大似然估计为每位参与者分配最可能的亚型和分期。主要模型设定最大亚型数为3（N_Smax = 3），以ALSFRS-R总分为48分的参与者作为功能参考对照组。此外，研究人员还进行了复杂度分析（允许最多4个亚型）和敏感性分析（扩展对照组定义为ALSFRS-R ≥ 47分）以检验模型的稳健性。纵向稳定性通过Cohen's kappa系数和转移矩阵评估，功能衰退率以每30天的ALSFRS-R总分变化量量化，生存分析采用Kaplan–Meier曲线和Cox比例风险回归模型。
三、研究结果
3.1 基线特征
研究共纳入3067名ALS患者。基线分析显示，三个SuStaIn亚型在年龄和性别分布上存在显著差异：S1（上肢为主型）参与者显著年轻于S2（下肢为主型）和S3（延髓为主型）。基线SuStaIn分期在亚型间存在显著差异（P < 0.001），其中延髓为主型亚型的平均分期（8.2 ± 6.3）低于上肢为主型（9.2 ± 6.5）和下肢为主型（8.8 ± 5.7）。尽管分期存在差异，三个亚型的ALSFRS-R总分并无显著差异，表明基线时总体功能损害程度相当。ALSFRS-R条目水平分析揭示了亚型特异性的功能受累模式：S1在上肢条目（Q4–Q7）得分显著较低，S2在下肢条目（Q8–Q9）得分显著较低，S3在延髓条目（Q1–Q3）和呼吸条目（Q10–Q11）得分显著较低。这些模式与临床上公认的上肢为主型、下肢为主型和延髓为主型表型高度对应，支持了数据驱动亚型的临床可解释性。值得注意的是，Q12（呼吸功能不全）在亚型间无显著差异。
3.2 亚型与分期推断
SuStaIn模型识别出三种具有不同进展模式的亚型。在上肢为主型（S1）中，早期受累集中于上肢条目（Q4–Q6），随后较快进展至邻近运动域（Q7–Q9），延髓条目（Q1–Q3）和呼吸条目（Q10–Q12）在较晚阶段才受累。在下肢为主型（S2）中，下肢条目（Q8–Q9）最早受累，随后进展至上肢条目（Q4–Q7）。在延髓为主型（S3）中，延髓条目（Q1–Q3）在最早阶段即受累，随后出现多域广泛受累；值得注意的是，呼吸困难（Q10）在该亚型中相对较早出现。复杂度分析（允许4个亚型）显示，下肢和延髓亚型得以保留，而上肢亚型被细分为两个密切相关的亚型，两者均以上肢条目早期受累为特征；由于未出现明确的新临床表型，为保持可解释性并避免过拟合，研究保留三亚型模型作为主要模型。敏感性分析（扩展对照组为ALSFRS-R ≥ 47分）显示亚型结构、分期进展模式和亚型比例与主要分析一致，表明结果在不同模型设定下具有内部一致性。
3.3 纵向稳定性
在180天队列（n = 2637）中，63.0%的参与者保持亚型不变，对应中等程度的一致性（Cohen's κ = 0.42）；分期进展的单调性（稳定或进展）达85.0%。在365天队列（n = 1506）中，亚型稳定性较低，53.8%的参与者保持亚型不变，对应一般程度的一致性（Cohen's κ = 0.26）；分期进展的单调性达90.8%。亚型稳定性随时间下降可能反映了疾病进展过程中临床表型的趋同——随着功能损害向多域扩散，晚期阶段亚型间的区分度降低。
3.4 各亚型ALSFRS-R衰退速率
在180天和365天两个队列中，亚型间的ALSFRS-R衰退速率均存在显著差异（Kruskal–Wallis检验，p < 0.001）。事后分析显示，延髓为主型（S3）在两个队列中的衰退速率均显著快于上肢为主型（S1）和下肢为主型（S2）。在180天队列中，S1的衰退速率显著快于S2，但在365天队列中S1与S2之间无显著差异。这些发现在调整基线年龄和性别后仍然成立。
3.5 生存分析
Kaplan–Meier生存曲线提示S3较S1有更快的生存下降趋势，但差异未达统计学显著（log-rank检验，p = 0.851）。在调整年龄和性别的Cox比例风险模型中，亚型与死亡率无显著关联；而较高的SuStaIn分期与死亡风险增加显著相关。在包含亚型、分期、年龄和性别的多变量模型中，较高的SuStaIn分期仍与死亡风险显著相关，而亚型未显示显著关联。当按中位SuStaIn分期（分期=8）分层后，低分期队列中观察到亚型间的显著生存差异（log-rank检验，p = 0.008），其中S1与S2之间以及S1与S3之间存在显著差异；高分期队列中未观察到显著亚型差异（log-rank检验，p = 0.271）。研究人员指出，这些事后分层和配对分析应视为探索性发现，因为记录死亡数有限且多重比较可能增加假阳性风险。
四、讨论与结论
本研究的核心贡献不在于发现新的ALS亚型或替代传统的发病部位分类，而在于利用广泛可及的功能量表重建了有序的功能损害模式，并推导出与死亡率相关的疾病分期。延髓为主型亚型更快的功能衰退与既往临床观察一致，提示延髓受累与更具侵袭性的疾病进程相关。该亚型在基线时呼吸困难（Q10）和端坐呼吸（Q11）损害显著更重，且这些呼吸症状在疾病进展的中期阶段（分期15–25）即出现；尽管肺功能数据仅覆盖部分参与者，延髓为主型亚型在基线时仍呈现呼吸功能下降的趋势（p = 0.054），提示应重视该亚型患者的早期呼吸功能评估。SuStaIn分期与死亡率的关联不应被视为完全独立于功能损害，其价值在于基于条目水平异常的分布和顺序提供疾病累及范围的结构化估计，从而补充ALSFRS-R总分等传统汇总指标。亚型相关的生存差异仅在低分期队列中显现，可能反映了基线疾病分期异质性对整体队列中亚型差异的掩盖效应；延髓为主型患者因症状更易识别而较早进入临床评估，这一领先时间偏倚也可能削弱整体队列中的亚型生存差异。由于死亡率数据仅覆盖队列的15%，整体队列中未观察到亚型与生存的显著关联也可能受限于统计效能。
研究存在若干局限性：PRO-ACT数据库的回顾性性质可能引入选择偏倚；2688名参与者因ALSFRS-R条目数据不完整而被排除，可能影响亚型分布和进展轨迹的估计；疾病持续时间数据不可得，呼吸和死亡率数据不完整；可靠的发病部位信息不一致，因此SuStaIn衍生亚型应被解释为ALSFRS-R条目水平的功能受累模式而非传统发病部位分类的直接替代；功能参考组规模小且由ALS患者而非健康对照组成；ALSFRS-R作为半定量临床量表可能受主观性和评估者间变异影响。此外，所识别的亚型可能部分反映ALSFRS-R量表固有的解剖-功能结构——该量表条目本身即评估延髓、上肢、下肢和呼吸域——因此解剖学可识别模式的出现并不应被单独解释为存在不同生物学亚型的证据，未来需整合影像学、神经电生理、遗传学或生物标志物数据以进一步细化亚型特征。
研究结论：SuStaIn在ALS中重建了三种具有临床可解释性的功能进展轨迹，分别以上肢、下肢和延髓受累为主。延髓为主型亚型与更快的功能衰退相关，较高的SuStaIn分期与增加的死亡风险相关。未来需与既定预后指标进行前瞻性比较，以确定该框架能否提供额外的预后信息并支持临床决策。
