跨生命周期颈脊髓形态计量学图表：基于大规模健康队列的规范化参考框架

《Human Brain Mapping》：Charting Cervical Spinal Cord Morphometry Across the Lifespan

　　

脊髓形态计量学（spinal cord morphometry）提供了反映神经健康的重要生物标志物，但临床解读常因个体间变异以及缺乏覆盖整个人类生命周期的规范参考而受到干扰。为填补这一空白，研究人员生成了首个全面的颈脊髓形态计量学生命周期图表。研究人员从78,269次T¹加权脑MRI扫描中，从33个数据集（年龄0–100岁）筛选出最终横断面队列，包含24,251名健康或认知正常对照者，且均具有颈髓覆盖。为克服对比度变异性，研究人员采用最先进的对比度无关深度学习分割方法，提取C₁至C₇的横截面积（cross-sectional area, CSA）、前后径（anteroposterior, AP）和左右径（right–left, RL/横径）以及形状指标（AP/RL比值、离心率和实心度）。研究人员使用位置、尺度和形状广义可加模型（Generalized Additive Models for Location, Scale, and Shape, GAMLSS）对规范轨迹进行建模。所得图表揭示了明显的非线性生命周期变化：儿童期至青春期快速生长，成年早期至中期达到成熟高峰（例如CSA在35岁左右），随后逐渐下降。研究量化了颈髓沿头尾方向显著的区域差异以及一致的性别差异（男性尺寸指标大于女性）。脊髓轨迹与脑白质和脑干体积表现出强烈的时间耦合，提示中枢神经系统发育与老化的整合性。这些生命周期图表提供了一个稳健的规范框架，能够对个体脊髓形态计量学进行年龄和性别特异性的百分位评分。该资源为区分典型变异与病理改变提供了关键工具，增强了脊髓MRI在发育和神经退行性疾病研究中的临床实用性。

论文解读：跨生命周期的颈脊髓形态计量学图表

脊髓形态计量学（spinal cord morphometry）提供反映脊髓宏观结构完整性的定量指标，如横截面积（cross-sectional area, CSA）、前后径（anteroposterior, AP）、左右径（right–left, RL）等，在多种神经系统疾病中具有重要临床价值。然而，个体间生物学差异、不同扫描仪与采集协议的方法学差异以及椎体节段定位差异导致测量变异性大，长期缺乏覆盖全生命周期的规范参考。既往研究多受限于窄年龄范围、较小样本量、单一椎体水平或简化统计模型。为克服这些限制，研究人员利用大规模脑MRI数据集，构建了首个颈脊髓形态计量学生命周期图表。

研究人员从78,269次T¹加权MRI会话中，经健康对照筛选、视觉质量控制和每名参与者单个合格扫描选择，最终纳入来自33个数据集的24,251名健康或认知正常对照，年龄0–100岁。采用对比无关深度学习分割提取C₁–C₇共六种形态特征，并用GAMLSS进行规范建模，揭示了非线性生命周期轨迹、区域差异、性别差异以及与脑体积的关系。该资源使个体测量可转化为年龄和性别调整的百分位评分，有助于区分典型变异与病理变化，增强脊髓MRI在发育和神经退行性疾病研究中的临床实用性。

关键方法方面，样本来自已聚合的大规模公开和本地数据集，最终横断面队列包含24,251名0–100岁健康或认知正常对照。图像处理使用Spinal Cord Toolbox v7.1中的对比无关深度学习模型自动分割脊髓并标记椎体；在原生空间中提取CSA、AP直径、RL（横径）直径、AP/RL比值、离心率（eccentricity）和实心度（solidity）。脑体积使用FreeSurfer v7.2.0提取。统计建模采用GAMLSS，以广义伽马分布拟合每个椎体水平和每种形态特征，年龄效应用分数多项式刻画，性别作为位置和尺度的固定效应，数据集作为随机效应。通过重复随机扫描选择和留一数据集法评估模型稳定性，并使用最大似然估计数据集特定偏移以进行样本外对齐。

研究结果如下：

3.1 跨生命周期非线性动态：在代表性C₃水平，CSA、AP、RL直径在儿童期至青春期快速增长，成年早期平台，老年期逐渐下降；AP/RL比值随年龄下降，离心率升高，实心度晚期略降。该结论基于GAMLSS拟合的百分位曲线与原始数据的对照。

3.2 区域变异：沿C₁–C₇层面，尺寸指标一般尾侧减小，C₄–C₅出现明显颈膨大，主要由RL直径驱动；形状指标具有椎体水平特异性的复杂轨迹。该发现通过对各椎体水平分别建模获得。

3.3 生命周期里程碑：CSA峰值出现在33–42岁（C₃为35岁），AP直径峰值一般早于RL直径峰值。该结论由拟合中位曲线的峰值年龄得出。

3.4 年化变化率：变化率在生命早期最大，如CSA在婴儿期增加约21%–75%/年、幼儿期增加6.7%–9.5%/年；晚期CSA下降0.23%–0.37%/年，AP直径下降速度大于RL直径。该数值由拟合曲线的一阶导数推导。

3.5 群体变异性：CSA的变异系数在C₁–C₅为8%–12%，C₆–C₇较高；AP和RL直径变异系数均低于10%，且随年龄大致稳定。该结论通过按生命阶段计算变异系数获得。

3.6 性别差异：男性CSA和直径普遍大于女性，形状指标性别差异较小，提示需要性别分层图表。该结论通过将性别作为GAMLSS位置和尺度固定效应建模获得。

3.7 脑–脊髓关系：脊髓CSA轨迹与脑白质和脑干体积轨迹最相似；CSA与脑体积的相关性在儿童期最强，并从头侧到尾侧递减。该结论通过比较中位轨迹和分阶段Pearson相关分析获得。

3.8 样本外对齐：在未用于构建图表的多发性硬化（multiple sclerosis, MS）数据集VUIIS-MS中，对齐后对照平均百分位为50.4、中位为48.8，MS参与者平均百分位为35.4、中位为38.3；人为引入5%偏差的实验证实数据集偏移可被校正；阿尔茨海默病神经影像学倡议（ADNI）队列中，阿尔茨海默病患者14%低于第5百分位。该结论通过控制组估计偏移并应用于患者获得。

3.9 模型稳定性：重复随机扫描选择的中位变异系数为0.05%，留一数据集法为0.12%，C₆–C₇及稀疏采样年龄变异较大。该结论通过敏感性分析获得。

讨论部分总结：这些图表为临床和科研提供了将绝对测量转换为年龄、性别调整百分位的工具，可在临床中用于识别低于第5百分位的非典型萎缩，并在研究中减少方差、提供标准化终点。图表显示脊髓结构成熟持续至30多岁，与中枢神经系统髓鞘化延续至成年早期的文献一致；老年期AP直径下降快于RL，与脑衰老类似。脑–脊髓强相关性表明中枢神经系统以整合方式发育和老化，并支持按性别分层。局限性包括横断面设计、C₆/C₇样本较少、婴幼儿期采样稀疏、基于椎体而非功能节段的参考框架、分数多项式可能难以捕捉急剧转变、广义伽马形状参数设为常数、各特征独立建模等。未来需纳入纵向数据、胸腰段脊髓和临床人群验证。

研究结论翻译：通过使用大规模脑影像队列和对比无关分割方法，研究人员生成了首个描述人类生命周期中颈脊髓形态计量学非线性动态的全面规范性图表。这些图表提供了典型生长、成熟、老化、性别差异和群体变异性的参考数据，揭示了脊髓与脑白质及脑干平行的发育过程。类似于既定的生长图表，这些脊髓图表提供了一个将个体测量值转换为年龄和性别调整百分位数评分的框架，以及一个用以基准化脊髓健康并检测提示神经发育或神经退行过程的非典型变异的标准化工具。这项工作为在临床实践和研究中更精确、个体化地评估脊髓完整性奠定了基础，有望提高诊断准确性、增强临床试验统计功效，并加深对中枢神经系统发育和衰老的理解。未来纳入纵向数据、胸/腰段脊髓和临床人群的扩展将进一步完善这些重要的神经解剖资源。生命周期图表、模型拟合和样本外对齐代码已公开。