围产期人脑发育中丘脑拓扑的稳定性：基于扩散MRI和高斯混合模型（GMM）聚类的丘脑分区研究

《Brain Topography》：Stable Thalamic Topography Across Perinatal Development: a Diffusion MRI Study

　　

丘脑是脑连接的关键中继枢纽，但其在早期人脑发育中的拓扑组织仍未完全理解。该研究通过优化并应用聚类框架对新生儿丘脑进行分区并描绘其发育轨迹，从而刻画关键的围产期发育阶段。研究人员使用基于高斯混合模型（Gaussian Mixture Model, GMM）的聚类研究丘脑拓扑，该聚类应用于从扩散磁共振成像（diffusion MRI, dMRI）数据获得的体素级纤维取向分布函数。使用三个围产期数据集：来自发育人类连接体项目（Developing Human Connectome Project, dHCP）的足月出生婴儿（n = 68，MRI时年龄：37.4–44.7周）、dHCP早产儿数据（n = 48，29.3–37.4周）以及来自Zürich的新生儿对照数据集（n = 42，39.0–48.7周）。聚类使用成人重测扩散MRI数据集（n = 30）优化。研究人员将丘脑细分为七个聚类，其中GMM显示中等可重复性（重测重叠范围为0.65至0.73），并与数字化Morel图谱具有可比的解剖对齐。稳健回归分析显示，各丘脑聚类的绝对体积随年龄显著增加（FDR校正p < 0.001）。通过相对体积表征的拓扑保持稳定，提示妊娠晚期丘脑拓扑得以保留。两个聚类，分别与腹前核（ventral-anterior nucleus, VA）和枕核群（pulvinar group, PuG）密切对应，显示出中等程度的年龄相关拓扑变化（FDR校正p < 0.05，adjusted R² = 0.039和0.101）。对协调化体积的分组分析显示，早产儿与足月出生新生儿之间无显著组效应。未观察到性别差异。研究发现表明，丘脑拓扑在围产期大部分已建立；然而，早产儿缺乏系统性的足月等效扫描限制了研究人员就早产影响得出结论的能力。解释可能受到扩散MRI特征的间接性质以及聚类技术中等可重复性的限制。

研究背景与问题
丘脑（thalamus）是神经环路发育的关键结构，参与轴突引导和神经元迁移，并作为感觉中继调节皮层区域间信息交换，支持感觉运动、认知、情绪和注意功能。丘脑发育紊乱与双相障碍、自闭症、认知障碍、情绪功能障碍等发育或精神疾病相关。早产也与丘脑发育和组织改变相关。早期脑发育中，丘脑建立拓扑结构，为神经通路和与皮层的相互连接提供支架，显著影响皮层发育。人类脑发育具有前额叶连接丰富、背内侧核神经元数量多等独特方面，难以从其他物种外推，且病理如何影响丘脑发育仍需人类数据研究。新生儿丘脑环路发育尤其敏感；早产儿显示丘脑体积减少，可能反映早产对丘脑成熟和整体脑结构的影响。传统丘脑核团依靠组织学和手动分割，MRI尤其扩散磁共振成像（diffusion MRI, dMRI）和功能磁共振成像（functional MRI, fMRI）实现体内分区；fMRI显示功能特化但与组织学图谱对齐不佳。扩散张量成像（diffusion tensor imaging, DTI）广泛用于丘脑分区，但发育脑前髓鞘化轴突降低扩散各向异性，削弱长程连接预测准确性，且DTI主要方向可产生空间不连续区域。约束球面去卷积（constrained spherical deconvolution, CSD）尤其单壳三组织CSD（single-shell 3-tissue CSD, SS3T-CSD）估计白质纤维取向分布（white matter fiber orientation distribution, WM-FOD）并显式建模灰质和脑脊液，更适合高自由水新生儿脑。围产期和产后早期感觉刺激快速增加，影响丘脑皮层连接或丘脑结构发育，但丘脑拓扑在该期如何变化尚不完全清楚。因此，研究人员旨在以局部扩散属性作为丘脑连接和结构的代理，描述丘脑拓扑组织发育轨迹；基于既往验证，假设基于WM-FOD的数据驱动聚类足以反映早期产后发育中的丘脑拓扑，并进一步优化聚类以生成解剖有意义的分区。

研究人员使用成人重测数据集优化聚类，在dHCP足月、dHCP早产和Zürich新生儿数据中应用，将丘脑分为七簇，GMM中等可重复，与Morel图谱可比。绝对体积随年龄显著增加，相对体积拓扑稳定，VA和PuG有中等年龄相关拓扑变化；早产与足月无显著组效应，无性别差异。结论：丘脑拓扑在围产期大部分已建立；缺乏早产儿系统足月等效扫描限制结论，解释受dMRI间接性和聚类中等重复性限制。意义：框架跨扫描仪、发育阶段和早产儿泛化，支持丘脑成熟和病理研究。

关键技术方法
数据来源：成人重测数据来自Boekel et al.的重测库，最终30例；dHCP neonatal release 2采样120例，排除4例后116例，含48例早产和68例足月；Zürich新生儿来自University Children’s Hospital Zurich的Heart and Brain Research Group Zurich，42例。方法：dMRI经SS3T-CSD估计WM-FOD，以球谐（spherical harmonics, SH）阶0、2、4、6、8表示；融合空间坐标与WM-FOD特征，采用K-means和GMM聚类，贝叶斯优化空间权重α与SH尺度s；Morel图谱重映射，Dice相似系数（Dice similarity coefficient, DSC）和Hausdorff距离评估；稳健线性回归、ComBat-GAM协调化、FDR校正及队列/性别效应分析。

研究结果
参数优化（Parameter Optimization）：通过成人重测数据优化聚类，K-means与GMM均可应用，并以DSC与Morel图谱评估。结论是高空间权重一般改善对齐，但中等SH贡献重要；SH阶4和6优于8。

基于K-means聚类与Morel图谱之间DSC的网格搜索（Grid Search Based on DSC Between K-means Cluster and the Morel Atlas）：网格搜索显示，高空间权重改善对齐；低SH阶4和6在较宽空间权重范围DSC较高，SH阶8无额外收益。

贝叶斯优化（Bayesian Optimization）：固定最大空间权重0.9、0.8、0.7、0.5；最优DSC可在空间权重显著低于1时实现；α约0.5、s约10^?2平衡空间邻近和扩散特征；空间和SH距离等贡献DSC约0.7。

扫描-重扫描可重复性（Scan-rescan Reproducibility）：GMM中度重复，簇平均DSC为0.65–0.73，变异较大；silhouette和gap统计倾向两簇，但最终七簇为先验的、基于Morel核团分组、生物可解释分区。

组织学图谱定义的丘脑核团与基于WM-FOD的聚类之间的相似性（Similarity Between Histology Atlas-based Thalamic Nuclei and the WM-FOD-based Clusters）：在重测数据，Cluster 2与内侧背侧与中央组（Mediodorsal and Central, MDC）一致强对齐，平均DSC为0.78±0.11；Cluster 6与枕与膝状体组（Pulvinar and Geniculate, PuG）中度重叠0.71±0.06；Cluster 4与腹外侧与腹内侧组（Ventrolateral and Ventromedial, VLM）部分对齐0.52±0.09；其他分布或最小重叠。GMM较K-means与解剖核团对齐更低更可变，尤其MDC和PuG外。

围产期丘脑拓扑发育轨迹的表征（Characterizing Developmental Trajectories of Thalamic Topography in the Perinatal Period）：协调化绝对体积稳健回归显示，早产dHCP中Cluster 2（MDC）、3（腹前组，Ventral Anterior, VA）、4（VLM）、6（PuG）显著正相关，VA调整R²=0.263。足月dHCP中2、3、4、5（腹后组，Ventral Posterior, VP）、6显著，VA最大R²=0.42。Zurich中4、5、6显著，效应较小；Cluster 5仅5例。合并三队列所有七簇显著增长（FDR校正p<0.001），VA的R²=0.610，PuG为0.521，MDC为0.367。相对体积大体稳定；合并仅Cluster 3（VA）、6（PuG）、7（外侧组，Lateral, L）随龄显著（FDR校正p<0.05，R²为0.039–0.101）；单队列除Zurich Cluster 5（VP）（n=5）外无显著。队列主效应、年龄×组交互不显著；性别主效应及年龄×性别交互FDR校正后不显著。结论：丘脑核团绝对体积增长，相对比例在围产期大体保留。

讨论总结与研究结论翻译
讨论：围产期和产后早期因感觉刺激快速增加及新生儿适应新环境而关键。绝对体积与扫描年龄正相关，符合脑生长、丘脑皮层（thalamocortical, TC）连接成熟和白质体积增加。相对体积稳定支持丘脑拓扑组织早期建立，在TC回路初始形成后可能稳定；与握手假说及人妊娠18–22周TC连接建立时间线一致，提示突触精化和髓鞘化均质影响核团，未导致围产期拓扑转移。早产与足月及Zürich新生儿协调化后无显著组效应，应解释为年龄和样本特异性观察，而非证明早产对丘脑拓扑无影响；缺乏早产儿系统足月等效扫描限制比较。虽然TC连接受早产影响，但该研究反映拓扑而非连接强度，方法可能不足以捕捉轻微早产相关变化。特定簇中VA相对体积随龄增加、L随龄降低，统计中等，可能受聚类不确定性混杂；与MD核发育性神经元过度生长间接证据一致，MD核主要前额叶连接且可能与VA邻近聚类，可部分解释随孕龄变化，作为丘脑-前额叶连接选择性重组的间接标志。研究未统计体素GMM分配确定性与年龄或队列关系，不能得出结论；所有体积分析基于硬最大概率标签。MRI聚类生物有效性有争议。该框架使用局部扩散和空间特征；与主特征向量或完整45维FOD相比，与Morel对齐和重复性相似；仅局部扩散无空间降低对应和重测可靠性；中阶SH特征加空间邻近改善生物合理性，高阶SH附加价值有限。生物有意义分区无需很高角分辨率，可适用于临床dMRI，如Zürich婴儿队列；框架跨扫描仪、发育阶段和早产儿泛化。与图谱比较仅部分核团良好对应，MDC、PuG实质重叠，Cluster 4和3混合对应；差异可能反映Morel基于细胞结构与聚类基于局部扩散属性，以及基于距离的标签分配可能偏好较大核团；总体支持生物效度。七簇选择结合生物合理性和探索统计验证；与既往七主要核团组一致，统计指标不一致：elbow法一般6或7簇，silhouette和gap一致倾向2簇，5簇中等但无一致支持；最终七簇为平衡稳健方案。局限性：排除零体积簇可能选择偏倚；合并dHCP和Zürich协议和招募差异，ComBat-GAM协调但未做留一队列敏感性分析；Morel为成人脑，配准误差和缺乏新生儿特异图谱；无专门足月等效早产队列，发育年龄效应和早产相关效应不能完全区分。早产不良影响可能数天或数周后出现，需与足月等效年龄足月对照比较。

研究结论翻译：该研究发现表明，丘脑拓扑在围产期大部分已建立；然而，早产儿缺乏系统性的足月等效扫描限制了研究人员就早产影响得出结论的能力。解释可能受到扩散MRI特征的间接性质以及聚类技术中等可重复性的限制。研究人员的主要发现是围产期各年龄丘脑体积的稳定性，提示神经连接早期建立，与发育理论一致，正在进行的突触发生等发育过程对丘脑拓扑影响相对较小。