白质纤维束集合水平轴突半径分布与传导速度预测对沿轴口径变异的稳健性：二维组织学采样验证及人类样本量指导

《Communications Biology》：White matter fiber bundles enable robust conduction and time-critical signaling despite axon caliber variation

　　

轴突口径（axon caliber）决定信号传导速度（conduction velocity），并在构成脑内布线（wiring）的纤维束内部及纤维束之间变化。研究人员对轴突组织的认识主要来自二维（2D）组织学，其假设轴突为圆柱形且口径由半径描述。然而，近期三维（3D）组织学揭示半径沿轴存在显著变异，从而对基于2D的轴突半径与传导速度推断提出疑问。研究人员在来自雄性大鼠的4.5×10⁵个3D轴突重建中显示，尽管存在这种个体变异，轴突纤维束在集合水平（ensemble level）仍具有稳定的半径分布。这些分布可被2D横截面忠实捕获，并可产生稳健的传导速度预测，因为沿轴变异仅具有中等影响。大轴突中的传导尤其稳定，支持时间关键性信号传递。在验证2D采样后，研究人员利用来自2D组织学的4.6×10⁷个人类胼胝体（corpus callosum）轴突，确定神经科学应用中的样本量需求。研究结果表明，轴突组织如何产生稳定的集合统计量与传导，从而解释并强化数十年来基于2D组织学的研究。

研究背景：脑白质通过纤维束连接远隔脑区，纤维束由排列的髓鞘化轴突组成；轴突口径（axon caliber）是传导速度（conduction velocity）的关键决定因素，粗轴突传导快但薄轴突更多且代谢经济。轴突口径也具有神经疾病生物标志物意义，并原则上可通过磁共振成像（MRI）无创获取。既有认识主要来自二维（2D）组织学横截面，并以圆柱模型（cylindrical model）解释，即把轴突视为半径恒定的均匀管道；由此得到右偏半径分布、跨物种尾部差异及纤维束间系统差异，并与电生理和扩散磁共振成像（diffusion MRI）等测量一致。问题在于2D切片通常远窄于所采样纤维束，样本量可能不足以刻画长尾或验证参数模型；串行切片电子显微镜的三维（3D）重建又显示单个轴突偏离圆柱模型，沿轨迹存在半径变异和波动。此前仅有约50根大轴突的线索提示半径分布沿纤维束稳定，尚需大规模定量验证；沿轴变异对扩散MRI和传导速度预测的影响也不清楚，尤其不同大小轴突是否差异受影响。因此研究人员利用最大规模3D大鼠轴突重建和2D人类数据，检验沿轴半径变异的后果、2D与3D采样稳健性及参数描述。研究显示：尽管单个轴突半径沿轴变化，集合水平半径分布沿纤维束稳定，2D横截面能忠实代表3D纤维束性质并给出稳健传导速度预测；大轴突传导尤其稳定；人类长尾需要更大样本；参数模型难以准确刻画长尾。这为长期2D组织学结论提供直接验证，并为采样设计提供指导，意义在于解释并强化数十年2D组织学研究，同时提示扩散MRI解释需考虑沿轴变异。

关键方法：研究人员再分析两个公开数据集。大鼠白质3D电子显微镜数据来自5只成年雄性Sprague-Dawley大鼠（3只创伤性脑损伤（traumatic brain injury, TBI），2只假手术对照），共10个体积、19个感兴趣区（region of interest, ROI）、约4.5×10⁵个髓鞘化轴突；人类胼胝体2D光学显微镜数据来自2名死后供者（1男61岁，1女60岁），35个ROI、约4.6×10⁷个轴突。分析包括沿轴半径轮廓、模拟2D横截面、Wasserstein距离（Wasserstein distance, W）、算术平均半径r?、MRI可见半径r_MRI、子采样和Akaike信息准则（Akaike information criterion, AIC）参数拟合。

研究结果：（1）沿轴半径变异对传导速度影响中等（Along-axon radius variation has modest impact on conduction velocity）。通过3D半径轮廓计算系数变异（coefficient of variation, CoV），中位为0.27，最高可达0.6；CoV随平均半径增加而降低，粗轴突的绝对半径变异饱和。模型估计传导速度中位降低约4%，沿轴扩散降低约21%；大轴突变异较小、传导更稳定，支持时间关键性信号传递。（2）2D横截面可捕获轴突半径分布（2D cross-sections capture axon radius distributions）。通过沿纤维束追踪局部分布并比较2D横截面与3D汇集分布，W显示约0.04 μm系统偏移，2D偏向较小半径；ROI内2D与3D差异小，且小于跨ROI解剖差异；r?和r_MRI在2D与3D间强相关（R≥0.97），但2D值分别低约12%和10.5%；以随机抽样比较表明差异来自2D短轴近似而非横截面内空间结构。结论是2D分布稳定代表3D，但存在几何近似偏差。（3）样本量需求在分布主体与尾部指标间不同（Sample size requirements differ for bulk and tail metrics）。通过大规模参考分布子采样，n=10³时分布主体可较好捕获，但人类尾部需约10⁵；n=10²时采样误差可超过ROI间解剖变异；r?无偏且变异<5%，r_MRI在小样本显著低估（人类n=10²约?30%，大鼠约?10%），收敛所需样本大鼠约10³至10⁴、人类约10⁵。结论是主体指标样本需求较低，尾部指标需更大样本且物种依赖。（4）参数分布遗漏长尾（Parametric distributions miss long tails）。通过大样本拟合六种候选参数分布并以AIC排序，大鼠拟合紧密，人类拟合发散且多数不能同时表示众数和尾部；广义极值分布（generalized extreme value, GEV）在多数ROI胜出，但形状偏差和尾部指标误差仍大；r?误差<5%，r_MRI误差大，人类多数分布超过25%，GEV在人类尾部误差最大；Gamma分布亦表现差。结论是AIC选择不保证尾部准确，参数建模尤其尾部应用需谨慎。

讨论总结：2D代表性源于沿轴变异在轴突间不相关且局部波动无系统趋势，2D横截面采样与随机采样统计等价；残余低估来自短轴近似，圆形等效半径可减小偏差但精度较低。传导速度的圆柱近似在口径变异下仍成立，约4%中位降低不淹没纤维束或物种差异；大轴突CoV约0.15，小轴突可达约0.30，猴子有限数据支持大轴突稳定，功能上作为时间关键活动高速通路，代谢上直径平方成本促使优化。扩散MRI中沿轴变异影响更复杂：半径分布与r_MRI可在圆柱模型下用2D组织学验证，但沿轴扩散信号可能被误归因于半径；沿轴形态本身也包含可测信息，需3D组织学验证。样本量指导为：主体约10³，尾部人类约10⁵，用于最大传导速度或r_MRI参考数据需尾部采样。参数分布可紧凑描述但应谨慎，尤其尾部。局限包括3D主要来自大鼠、猴子仅50根大轴突、TBI与对照ROI合并但2D与3D一致性接近恒等线；传导模型为局部圆柱序列和慢度积分，未纳入盐跃、节间长度、节点几何和轴突波动；固定收缩和2D分辨率限制也影响观测分布。研究结论部分翻译：总体而言，本研究重新评估了白质组织学的一个长期前提。传统上，研究人员使用2D组织学切片研究轴突半径分布，隐含假设轴突为圆柱形。虽然3D组织学已证伪该假设，但研究人员表明，尽管轴突口径存在变异，来自2D的集合水平结论——包括半径分布和传导速度预测——仍然有效。除强化数十年既有研究外，研究人员为神经科学应用中的高效组织学采样提供指导。