探索青少年自闭症谱系障碍中的非典型空间-功能耦合：来自神经发育和转录组学结构的见解

《Frontiers in Neuroscience》：Exploring atypical spatial-functional coupling in adolescent autism spectrum disorder: insights from neurodevelopment and transcriptomic architecture

　　

自闭症谱系障碍（ASD）与非典型大规模脑网络组织相关，但空间-功能依赖关系如何与临床特征和分子参考图相关联仍不完全清楚。为了量化空间功能异质性（Sill）和相干持久性（Range），研究人员分析了来自162名ASD和175名典型发育（TD）青少年的静息态功能磁共振成像（fMRI）数据，年龄均为12–18岁。与TD组相比，ASD青少年在高级联合网络内（包括左侧语言网络和右侧后多模态网络）的Sill显著增加，而Range的组间差异未通过多重比较校正。在ASD组内，升高的Sill选择性地与更严重的社会情感症状相关，但与受限和重复行为无关。为了探索潜在的生物学相关性，研究人员整合了皮质基因表达参考数据，并识别出与区域Sill差异相关的转录组学模式。这些基因在突触信号、线粒体过程和胶质相关功能中富集，突显了空间-功能组织与分子参考图之间的多尺度对应关系。综上所述，这些结果展示了ASD中改变的空间-功能特性、临床严重程度以及与突触信号、线粒体过程和胶质相关功能相关的转录组学特征之间的统计关联，为大规模功能组织提供了补充性的空间视角。

**自闭症谱系障碍（ASD）中空间-功能耦合的异常：神经发育与转录组学结构的整合视角**

**研究背景与问题**
自闭症谱系障碍（ASD）是一种以社交沟通差异和受限重复行为为特征的神经发育疾病，其行为表现和神经生物学组织存在显著异质性。既往神经影像学研究已揭示ASD涉及大规模脑系统（特别是高级联合皮层）的分布式改变，但传统功能连接分析主要关注区域间耦合强度，而对功能信号在皮质表面的空间依赖结构——即功能相似性如何随空间距离变化——知之甚少。变程图（variogram）方法可量化空间异质性（Sill，反映最大功能差异幅度）和相干持久性（Range，反映功能耦合持续的空间距离），为描述皮质组织提供无先验图谱的补充指标。目前尚不清楚：ASD是否伴随空间-功能依赖特性的系统性改变？这些改变是否优先定位于特定功能网络？它们与临床症状严重程度及分子参考图（如神经递质受体分布和基因表达）有何统计对应关系？为回答这些问题，研究人员利用ABIDE-I大型公开数据集，对青少年ASD（12–18岁）进行了多层级分析。

**研究设计与主要结论**
研究人员分析了来自10个扫描站点的162名ASD青少年和175名典型发育（TD）青少年的静息态功能磁共振成像（rs-fMRI）数据，采用变程图建模计算每个顶点的Sill和Range，并进行组间比较、症状关联、网络分布、功能解码及与神经递质受体和转录组学参考图的整合分析。主要发现：ASD组在双侧颞叶和右侧运动皮层表现出显著升高的Sill（P_FWE < 0.001），而Range无显著组间差异；升高的Sill与社会情感症状严重程度呈正相关（r=0.30，P_FDR=0.037），但与受限重复行为无关；网络水平分析显示Sill改变集中于语言网络和后多模态网络；Sill改变的空间模式与去甲肾上腺素受体密度及突触信号、线粒体过程和胶质功能相关基因表达显著对应。该研究发表于《Frontiers in Neuroscience》，为ASD的功能组织提供了空间统计学视角，并建立了宏观影像与分子参考图之间的多尺度统计联系。

**关键技术方法**（不超过250字）
本研究使用来自ABIDE-I联盟的162名ASD和175名TD青少年（12–18岁）的rs-fMRI数据，样本选自10个扫描站点（包括NYU、SDSU、TRINITY等）。主要技术方法包括：（1）变程图建模：利用BrainSMASH工具箱计算基于测地距离的顶点间经验变程，通过非线性最小二乘法拟合指数模型，提取Sill（渐近方差）和Range（相关距离）两个参数；（2）顶点水平组比较：采用SurfStat进行双样本t检验，经随机场理论进行家族误差（FWE）校正；（3）转录组学整合：从艾伦人脑图谱（AHBA，6名成人神经典型供体）获取微阵列基因表达数据，通过Abagen工具箱映射至Schaefer100图谱，采用偏最小二乘（PLS）回归识别与Sill差异相关的基因模块，并利用Lake单细胞数据集进行细胞类型富集分析；（4）神经递质受体建模：整合来自Hansen等（2022）的9种神经递质受体/转运体密度图，通过多元线性回归和岭回归评估其对Sill差异空间模式的解释能力。

**研究结果**
以下按照原文小标题总结，并说明通过什么研究得出什么结论：

**ASD-TD组间空间-功能依赖指标差异**
通过顶点水平双样本t检验（FWE校正），发现ASD青少年在双侧颞叶（左侧P_FWE < 0.001，Cohen’s d=0.53；右侧P_FWE < 0.001，d=0.50）和右侧运动皮层（P_FWE < 0.001，d=0.58）的Sill显著高于TD组，而Range无显著差异。敏感性分析（线性混合效应模型以站点为随机效应）确认结果稳健，头动指标（平均帧位移）与Sill无显著相关，排除了运动伪影的混淆。

**与常规rs-fMRI指标的对比**
将Sill和Range的组差异t图与局部一致性（ReHo）、低频振幅（ALFF）和度中心性（DC）的t图进行空间相关分析，发现弱相关，表明变程图指标捕捉了与这些传统指标部分不同的皮质功能组织特征。

**Sill改变区域与症状严重程度的关联**
在ASD组内，提取显著簇的平均Sill，控制年龄、性别、站点和平均FD后进行偏相关分析。结果显示Sill与社会情感评分（ADOS Social Affect）呈显著正相关（r=0.30，P_FDR=0.037），而与总分（r=0.25，P_FDR=0.051）和校准严重度评分（r=0.23，P_FDR=0.051）仅达到名义显著性，与受限重复行为（RRB）评分无显著关联（r=-0.06，P=0.611）。

**网络水平分布与功能特征**
将顶点Sill值平均至12个规范功能网络，经FDR校正后，ASD青少年在语言网络和后多模态网络中表现出显著更高的Sill（P_FDR < 0.05），提示改变集中于高级联合系统而非初级感觉网络。

**认知与功能关联**
将ASD–TD Sill差异图与来自Neurosynth的123个认知术语的第一主成分（PC1）进行空间相关分析（spin检验），发现显著正相关（r=0.33，P_spin < 0.001）。功能解码（基于24个代表主题术语）显示这些改变与语言、工作记忆、社会认知等认知域重叠。

**神经递质模型**
构建多元线性回归模型，以9种神经递质受体密度预测ASD–TD Sill差异图，模型解释了17.2%的方差（调整R2=0.172，r=0.30，P=0.002）。去甲肾上腺素是唯一通过FDR校正的显著预测因子（相对贡献18.6%，P_FDR=0.009），而乙酰胆碱、血清素和组胺仅达到名义显著性。岭回归分析进一步确认了去甲肾上腺素的稳健贡献。

**转录组学对应**
PLS回归发现PLS1解释35%的协方差。PLS+基因在细胞类型富集分析中显著富集于内皮细胞（End），PLS-基因富集于少突胶质细胞（Oli）。基因本体（GO）富集分析显示PLS+基因涉及线粒体内膜、蛋白质转运、突触结构等生物学过程。疾病富集分析显示PLS+基因与智力障碍、精神分裂症等神经发育/精神疾病相关，PLS-基因与阿尔茨海默病相关。所有结果均经FDR校正和spin置换检验确认。

**稳健性分析**
控制利手和智商等协变量后，Sill组差异模式基本不变。不同基因表达阈值（5%、15%）下GO富集结果高度一致，髓鞘相关类别始终显著。

**讨论与结论**
**研究讨论总结**
研究人员发现，ASD相关改变主要体现于Sill而非Range，反映了空间功能异质性而非广义空间范围的变化。这种改变定位于双侧颞叶和右侧运动皮层，集中于语言和后多模态网络，与社会情感症状而非RRB相关，提示空间-功能耦合异常与社交缺陷的潜在统计关联。神经递质模型显示去甲肾上腺素受体密度是Sill差异空间模式的最稳健预测因子，转录组学富集分析揭示突触信号、线粒体过程和胶质功能相关的基因表达模式，与已知ASD分子特征一致。然而，这些关联均为统计对应而非因果机制，且受限于横断面设计、多站点变异性、成人神经典型转录组数据等局限。研究强调，Sill作为反映功能相似性随空间距离衰减的指标，提供了与常规连接性分析互补的皮质组织视角。

**研究结论部分翻译**
综上所述，这些结果展示了ASD中改变的空间-功能特性、临床严重程度以及与突触信号、线粒体过程和胶质相关功能相关的转录组学特征之间的统计关联，为大规模功能组织提供了补充性的空间视角。