《Human Brain Mapping》：Major Depressive Disorder With Non-Suicidal Self-Injury in Adolescents: Altered Morphological Inverse Divergence Gradient in Brain Structure and Gene Expression Patterns
编辑推荐：
伴有非自杀性自伤（NSSI）的重性抑郁障碍（MDD）青少年面临严重的身心健康挑战。该状况的发展受遗传、环境和脑发育因素影响。本研究首次采用形态学逆散度梯度网络（MIND-Gd）探讨伴NSSI的MDD青少年的脑结构网络异常及其转录机制。共241名MDD青少年参与
伴有非自杀性自伤（NSSI）的重性抑郁障碍（MDD）青少年面临严重的身心健康挑战。该状况的发展受遗传、环境和脑发育因素影响。本研究首次采用形态学逆散度梯度网络（MIND-Gd）探讨伴NSSI的MDD青少年的脑结构网络异常及其转录机制。共241名MDD青少年参与研究，其中134名伴NSSI，107名不伴NSSI。利用高分辨率T1加权磁共振成像（MRI），研究人员构建主要MIND梯度网络。该网络反映不同脑区皮质形态的相似性。结果显示，伴NSSI青少年左侧尾侧中额回和双侧颞下回MIND梯度值较高；相反，左侧舌回和右侧顶上小叶值较低。这些结构改变与背侧注意网络（DAN）的功能改变存在空间相关。进一步分析显示，这些差异与九类神经递质受体分布相关，包括5-HT1A、GABAa、谷氨酸和多巴胺。研究人员还进行偏最小二乘回归分析（PLS）。第一成分（PLS1）解释MIND-Gd空间方差的39.53%。该分析识别与NSSI相关关键基因，如GABRA5和MGST1。富集分析显示，PLS1+基因参与突触信号、神经发育和心境障碍通路。这些基因富集于星形胶质细胞和皮质I、II、V层。另一方面，PLS1?基因与染色质调控和细胞骨架代谢相关。这些基因主要表达于皮质III和IV层。总之，伴NSSI的MDD青少年显示脑结构梯度的特定异常。这些异常与神经递质受体分布和基因表达模式密切相关。研究结果为理解NSSI神经生物学机制提供多层面证据。
研究背景与问题：MDD是慢性、复发性心理障碍，受遗传、环境和神经生物因素影响，是青少年疾病负担的核心来源之一。MDD青少年中NSSI发生率较高，NSSI亦是自杀行为的重要预测因素。临床上NSSI多依赖自伤行为和可见瘢痕识别，难以揭示其神经生物学机制。既有研究提示MDD和NSSI与皮质结构及脑网络异常相关，尤其涉及结构连接组，但伴NSSI的MDD患者特定结构连接组异常仍有限，其与不伴NSSI的MDD患者差异亦不清楚。既往DWI纤维束成像、结构协方差网络（SCN）、形态学相似网络（MSN）和形态学逆离散网络（MIND）等方法各有局限；MIND仅需T1加权MRI即可构建较稳定网络，且与遗传关系更强。皮质梯度可刻画基因表达、细胞结构、髓鞘分布和进化扩展等连续轴组织特征。既往研究多集中于MSN，NSSI相关异常及其对疾病相关基因表达的依赖仍待探索。因此，本研究比较伴与不伴NSSI的MDD青少年，构建MIND-Gd并联合转录组数据，以探索NSSI相关脑结构梯度异常、神经递质受体分布及基因表达模式。
研究开展与意义：研究人员测量区域和全局MIND皮质形态相似性，建立主要MIND梯度，比较两组梯度层级结构，并将梯度映射至功能网络和皮质细胞网络；进一步分析异常梯度与活体受体分布图的空间关系，采用PLS关联解剖基因表达模式，识别相关基因，并行GO/KEGG、细胞类型和皮质层富集分析。研究显示伴NSSI的MDD青少年存在特定MIND-Gd改变，且与神经递质受体分布及突触信号、神经发育、染色质调控和细胞骨架相关基因表达模式相关，为理解NSSI神经生物学机制提供多层面证据。
关键技术方法（不超过250字）：本研究为回顾性分析，样本来自山东精神卫生中心，纳入241名符合ICD-10的青少年MDD患者，其中NSSI组134例、非NSSI组107例；采集3.0 T T1WI，经FreeSurfer预处理，在D-K 308图谱提取灰质体积（GMV）、皮质厚度（CT）、表面积（SA）、平均曲率（MC）和沟深（SD），以Kullback–Leibler（KL）散度构建MIND相似矩阵，经BrainSpace扩散嵌入获得主要MIND梯度；以一般线性模型（GLM）比较组间差异；采用JuSpace进行神经递质受体空间相关；AHBA数据经abagen映射至左半球152个脑区×15632个基因，行PLS、GO/KEGG、皮质层及细胞类型富集，并用spin置换和稳健性分析验证。
研究结果：
3.1 数据样本与特征：通过质量控制纳入241名青少年，NSSI组134名，非NSSI组107名；临床与人口学资料及全局数据残差分析显示平均残差接近零，提示模型无偏。
3.2 MIND-Gd组间差异：通过区域主要MIND梯度比较发现，NSSI组左侧尾侧中额回、左侧颞下回、左侧中颞回和右侧颞下回梯度升高，左侧舌回和右侧顶上小叶降低；非NSSI组平均区域MIND梯度与组差异t图呈负相关（r=?0.708）。Yeo-7网络中DAN呈边缘显著降低，von Economo分类未见显著组间差异。
3.3 神经递质受体系统与MIND-Gd改变关联：通过JuSpace中30个受体图与组差异t图的空间相关分析，发现9类神经递质受体图显著相关，包括5-HT
1A、GABA
a、谷氨酸和多巴胺等。
3.4 转录-神经影像关联：通过PLS分析，PLS1解释39.53%的宏观结构方差，PLS1贡献73.12%方差；PLS1评分与组差异t图呈正空间相关（r=0.557）。研究识别2209个PLS1+基因和2110个PLS1?基因，并从PLS1+中识别与NSSI密切相关的GABRA5和MGST1；GABRA5表达与组差异t图正相关（r=0.545），MGST1呈负相关（r=?0.448）。
3.5 PLS加权基因注释：通过Metascape功能富集，PLS1+基因富集于突触信号和神经发育相关通路，PLS1?基因富集于染色质调控、细胞骨架组织和核酸代谢；疾病富集显示PLS1+与抑郁、癫痫等神经精神障碍相关，PLS1?与自闭症和发育迟缓等神经发育障碍相关。
3.6 皮质层富集：通过已发表皮质层标记基因分析，PLS1+基因显著富集于皮质I、II、V层，PLS1?基因显著富集于III、IV层，两组基因富集层无重叠。
3.7 特定细胞类型富集：通过七类中枢神经系统细胞分析，PLS1+基因富集于兴奋性神经元、星形胶质细胞和抑制性神经元，PLS1?基因富集于兴奋性和抑制性神经元；仅PLS1+基因与星形胶质细胞显著相关。
3.8 敏感性与稳健性分析：通过回归TIV、Spearman相关构建MIND-Gd、不同簇大小阈值、随机拆分NSSI亚组及基因类别富集分析等，组间差异和富集结果保持稳健。
讨论总结：本研究首次系统比较伴与不伴NSSI的青少年MDD患者脑结构和网络改变。MIND较可靠、与细胞结构匹配较好，对个体脑连接差异更敏感；区域上额颞梯度升高、顶枕梯度降低，提示结构分化与分离改变；视觉系统及DAN相关注意机制可能与NSSI相关，但DAN结果为边缘显著。受体分析显示组差异图与九类受体图空间相关，涉及血清素、GABA、谷氨酸、多巴胺和阿片等系统；但该分析基于成人规范性PET模板，应谨慎解释为空间关联而非患者受体活动直接证据。PLS1+基因富集于突触信号、神经发育和心境障碍相关通路，PLS1?基因富集于染色质调控、细胞骨架和核酸代谢；GABRA5和MGST1与NSSI相关，后者与炎症假说相关。皮质层和细胞类型分析提示兴奋性/抑制性神经元及星形胶质细胞参与。局限包括仅用T1和五种形态特征、样本量有限、回顾性设计缺少自杀意念、冲动性、童年创伤、焦虑和用药等信息、女性占多数、AHBA为6名成人供者且仅左半球数据，基因结果应视为探索性空间关联。
研究结论翻译：总之，研究结果揭示伴与不伴NSSI的青少年MDD患者主要MIND梯度存在显著改变，表现为额叶和颞叶区域梯度值升高、顶叶和枕叶区域梯度值降低。这些改变与神经递质受体分布存在空间相关，并富集于神经生物学相关基因通路。然而，鉴于横断面和相关设计，这些发现应视为初步和假设生成性；在推断任何临床意义前，需在具有全面临床表型的大样本前瞻性队列中独立重复。
