本篇文献分享发表于Progress in Neuro-Psychopharmacology and Biological Psychiatry期刊，第一作者为Peiyi Wu，共同第一作者为Ziyi Wang，通讯作者为YanqingTang。本公众号所发布内容旨在与大家分享学术新知，促进交流学习，版权归原作者或原出处所有，感谢各位学者的辛勤成果。

1.引言

青少年重型抑郁症 (AMDD)是一个全球性健康问题，影响全球约13%的青少年，是该人群疾病负担的主要原因。与成人起病抑郁症不同，青少年抑郁症出现在神经生物学显著重组的时期，此时突触修剪、髓鞘形成和神经递质系统成熟共同塑造认知和情感功能。这些发育变化造成了精神疾病易感性增强的时期。事实上，青少年起病的抑郁症与成人起病抑郁症相比，临床轨迹更严重，治疗结局更差。这种特异性易感性凸显了阐明青少年抑郁症独特神经生物学机制的迫切需求。

尽管对成人抑郁症神经生物学的研究取得了进展，但青少年抑郁症的机制基础仍未得到充分阐明。传统的神经影像学研究通常孤立地检查结构或功能改变，对神经系统的整合特性提供的见解有限。近年来多模态神经影像的进展突显了结构-功能耦合在维持脑网络效率中的关键作用，结构-功能(SC-FC)耦合代表白质连接与功能共激活模式之间的一致性。结构-功能耦合为表征脑网络的整合架构提供了一个强大框架，揭示了结构支架如何约束和实现功能动态。在健康发育过程中，结构-功能耦合遵循层级轨迹，在初级感觉区域增强，而在高阶联合网络中解耦合以支持认知特化。这种平衡的破坏与神经精神疾病有关，包括精神分裂症和成人抑郁症。然而，这些耦合模式在青少年抑郁症中是否以及如何发生特异性改变仍不清楚。

在成人重型抑郁症中，多模态研究开始阐明神经化学底物与结构-功能耦合异常之间的关系。例如，正电子发射断层扫描一致记录了前扣带皮层和前额叶皮层等区域中血清素转运体结合减少和多巴胺受体可用性改变。功能磁共振成像研究进一步表明，这些神经化学缺陷与情绪和认知控制网络内的结构-功能耦合受损相关。然而，这些研究通常受限于小样本量和横断面设计，未能完全阐明精确的机制联系。此外，发育中大脑的神经化学与网络耦合之间的关系尤其不明确。单胺能功能障碍——涉及血清素能、多巴胺能和去甲肾上腺素能信号——是抑郁症发病机制的核心。在青春期，这些系统经历实质性成熟，其中多巴胺能通路显示出特别显著的变化，影响奖赏加工和情绪调节。关键的是，这些发育性神经化学轨迹与结构-功能耦合异常之间的关系是理解青少年抑郁症的重要空白。

近年来多模态整合的进展，如JuSpace工具箱，使得将神经递质系统映射到神经影像表型成为可能，为研究神经化学如何与网络水平异常相关提供了一个框架。据据目前所知所知，尚无研究应用JuSpace工具箱来探究青少年抑郁症中神经递质系统与结构-功能耦合改变之间的关系。应用这一方法，可能有助于确定青少年抑郁症是否涉及反映其发育背景的特异性神经递质参与模式。

除神经递质系统外，基因表达模式同样影响脑网络组织和精神疾病易感性。然而，虽然转录组学研究已开始识别成人抑郁症中脑连接改变的分子特征，但青少年抑郁症中结构-功能耦合异常的分子基础仍不明确。青少年抑郁症可能涉及反映其独特神经发育背景的不同转录组谱，可能涉及控制突触修剪、髓鞘形成和关键期可塑性的通路。此外，青少年抑郁症中细胞类型特异性表达模式可能与成人有根本不同。可能涉及发育敏感的胶质细胞群和仍在精炼中的未成熟神经回路。理解这些发育性转录组差异对于制定适龄干预措施至关重要。

为解决这些关键知识空白，此研究进行了一项综合性的青少年抑郁症结构-功能耦合多模态研究，整合了弥散张量成像、静息态功能磁共振成像、空间神经递质相关分析和转录组学分析。此研究旨在：(1)表征青少年抑郁症中区域和网络水平的结构-功能耦合异常；(2)使用机器学习评估这些改变的诊断和临床效用；(3)通过空间相关分析与神经递质受体/转运体分布探究其神经化学相关性；(4)通过整合艾伦人脑图谱(AHBA)的全脑基因表达数据识别耦合异常的分子和细胞特征，并将这些特征映射到特定细胞群和发育时间窗口。

此研究假设青少年抑郁症与反映青春期持续神经发育过程的独特结构-功能耦合改变相关。具体而言，此研究预期这些改变与在此期间经历实质性成熟的神经递质系统相关，特别是与青少年神经发育中的基因表达谱相关。通过整合从宏观网络动力学到分子机制的多个分析层面，此研究为青少年抑郁症的神经生物学架构提供了见解。这种多层次表征可能揭示青少年抑郁症背后的神经发育机制，为更具发育针对性的干预提供信息。

2.方法

数据采集和分析流程展示在图1中。

图1分析流程概述。

2.1 参与者与临床评估

2.1.1 参与者

此研究共纳入205名12-18岁青少年，包括116名青少年重型抑郁症患者和89名健康对照。青少年抑郁症患者招募自两家医院的精神科：中国医科大学附属第一医院和盛京医院。健康对照通过广告从当地社区和学校招募。此研究遵循《赫尔辛基宣言》，并获得两家参与机构伦理委员会的批准。所有参与者和法定监护人在研究程序开始前均签署书面知情同意书。

抑郁症诊断由两名具有执业资格的精神科医生使用DSM-IV和DSM-V的结构化临床访谈执行。所有青少年抑郁症参与者均符合首次发作、未用药抑郁症的标准，且HAMD-17评分≥17。

2.1.2 临床评估

临床评估在入组后72小时内进行，包括一系列综合评估：

（1）症状严重程度评估：17项汉密尔顿抑郁评定量表、14项汉密尔顿焦虑评定量表、简明精神病评定量表、杨氏躁狂评定量表。

（2）认知和行为评估：MATRICS共识认知成套测验、自动思维问卷。

（3）情绪调节评估：认知情绪调节问卷。

2.2 MRI数据采集与预处理

静息态功能磁共振成像、T1加权图像和弥散张量成像使用两台3T扫描仪采集：GE Signa和Siemens PRISMA，在质量控制过程中因严重运动伪影（位移>2 mm）或信噪比低于协议规定阈值而排除20名参与者的数据。参与者被要求在扫描期间闭眼、保持清醒并清空思绪。使用泡沫垫和耳塞以减少头动和扫描仪噪音。每个站点的详细成像参数如下：

站点1：中国医科大学附属第一医院(n=170)，使用GE Signa 3.0 T扫描仪：

T1成像：3D快速扰相梯度回波序列（TR/TE = 7.1/3.2 ms，矩阵=240×240，视野=240×240 mm²，176层连续1 mm层厚，体素大小=1.0×1.0×1.0 mm³）

静息态功能磁共振成像：梯度回波平面成像（TR=2000 ms，TE=30 ms，翻转角=90°，层厚=3 mm，35层，无间隔，视野=240×240 mm²，矩阵=64×64，200个时间点，时长400 s）

弥散张量成像：单次自旋回波平面成像（TR/TE=17000/85.4 ms，矩阵=120×120，视野=240×240 mm²，65层连续2 mm层厚，25个方向，b=1000 s/mm²，1个b=0体积，体素大小=2.0×2.0×2.0 mm³）

站点2：中国医科大学盛京医院（n=35），使用Siemens PRISMA 3.0 T扫描仪：

T1成像：3D快速扰相梯度回波序列（TR/TE=2500/2.38 ms，矩阵=240×240，视野=256×240 mm²，192层连续0.8 mm层厚，体素大小=0.8×0.8×0.8 mm³）

静息态功能磁共振成像：梯度回波平面成像（TR=500 ms，TE=30 ms，翻转角=60°，层厚=3.5 mm，35层，无间隔，视野=240×240 mm²，矩阵=64×64，时长487 s）

弥散张量成像：单次自旋回波平面成像（TR/TE=2800/63 ms，矩阵=108×108，视野=240×240 mm²，75层连续2 mm层厚，64个方向，b=1000 s/mm²，1个b=0体积，体素大小=2.0×2.0×2.0 mm³）

静息态功能磁共振成像数据预处理包括：时间层校正和头动校正、归一化、去线性趋势、时间带通滤波、噪声回归和平滑。T1加权图像预处理包括：偏置校正、颅骨剥离、归一化至蒙特利尔神经研究所标准空间。弥散张量成像数据预处理包括：去噪、B1场校正、涡流校正、弥散张量图像与T1加权图像配准、张量模型拟合、确定性纤维追踪。

2.3 结构与功能连接的构建

网络节点使用人脑连接组图谱定义（210个皮层亚区，36个皮层下亚区）。结构连接通过确定性纤维追踪重建，使用连续追踪纤维分配算法，终止条件为转角>45°或各向异性分数<0.1。采集过程中调整弥散梯度方向以优化张量计算，生成246×246的结构连接矩阵，连接权重基于流线计数。通过双向纤维追踪验证矩阵对称性。对于功能连接，在空间配准后将相同的Brainnetome分区应用于功能空间。所有246个区域间的区域BOLD时间序列进行相关分析，皮尔逊相关系数使用Fisher变换转换为z分数。这种方法在保留个体拓扑变异的同时，能够直接进行结构-功能连接比较。

2.4 区域结构-功能耦合的计算

对于每位参与者，通过计算每个脑区与其他脑区结构连接与功能连接之间的相关性来获得区域结构-功能耦合。首先，排除自连接，将非零结构连接边归一化到高斯分布以减轻偏态连接权重的影响。区域耦合定义为每个节点的结构连接强度向量与功能连接强度向量之间的皮尔逊相关系数，仅计算非零结构连接的边。这生成每个受试者的246维向量，其中每个元素代表区域特异性结构-功能耦合。全局耦合计算为所有区域值的平均值。为评估网络层面的组织结构，脑区被划分为Yeo的7个功能网络。

2.5 结构-功能耦合的病例-对照分析

为评估青少年抑郁症组和健康对照组之间结构-功能耦合的组间差异，采用一般线性模型，以诊断组为自变量，同时控制性别、年龄、扫描站点和平均帧间位移作为协变量。随后进行青少年抑郁症组与健康对照组比较的双侧t检验。进一步分析Yeo7功能网络内的子网络，以表征网络特异性结构-功能耦合改变。采用错误发现率方法进行多重比较校正，显著性阈值设为p < 0.05。

2.6 诊断与症状预测

为评估结构-功能耦合改变的诊断价值，以显著组间差异的皮层区域的结构-功能耦合值为输入特征，构建机器学习模型。采用两种分类算法：极限梯度提升法（XGBoost)）和支持向量机（SVM)。模型性能通过敏感性、特异性和受试者工作特征曲线下面积指标进行评估，以表征青少年抑郁症诊断的判别能力。特征重要性通过跨交叉验证迭代的平均权重确定。

对于症状预测，采用与分类模型相同的区域结构-功能耦合特征集，构建随机森林回归模型。预测准确性通过预测临床评分与实际临床评分之间的斯皮尔曼相关系数量化。

2.7 结构-功能耦合与神经递质系统的相关性

为探究结构-功能耦合改变背后的神经化学基础，使用JuSpace（版本1.5）分析区域耦合差异图与神经递质受体/转运体分布之间的空间相关性，该工具箱是经验证的多模态神经影像整合工具。神经递质图谱来自正电子发射断层扫描和放射自显影研究，包括血清素（变体5HT1a、5HT1b、5HT2a、血清素转运体、SERT）、多巴胺（D1和D2受体、多巴胺转运体、DAT、6-氟-(18F)-L-3,4-二羟基苯丙氨酸，PDOPA）、γ-氨基丁酸（GABAa型，GABAa）、阿片类（μ-阿片受体, MU）、去甲肾上腺素（去甲肾上腺素转运体, NAT）、乙酰胆碱（囊泡性乙酰胆碱转运体, VAChT）和谷氨酸（代谢型谷氨酸受体5, mGluR5）。

在Brainnetome图谱定义的脑区之间，计算受体/转运体密度值与结构-功能耦合差异t值之间的斯皮尔曼秩相关系数。为确定统计显著性同时考虑神经影像数据中的空间依赖性，采用严格的置换检验方法（10,000次置换）并进行空间自相关校正。所得p值使用FDR方法进行多重比较校正，显著性阈值设为p < 0.05。

2.8 结构-功能耦合改变的转录组学相关性

2.8.1 区域基因表达数据获取

转录组谱来自艾伦人脑图谱，该图谱提供了来自6个死后成人大脑的微阵列表达数据，包含3702个空间注释的组织样本。区域基因表达数据使用abagen工具箱进行处理，该工具箱是专门用于处理AHBA微阵列数据的Python包。

鉴于AHBA数据处理的已知变异性，实施了标准化的预处理流程。由于AHBA数据集中右侧半球覆盖有限，该方法聚焦于左侧半球分析。处理流程生成了与此研究神经影像数据对齐的123个皮层区域中15,631个基因的标准化表达值。

2.8.2 空间协变分析

为揭示青少年抑郁症中结构-功能耦合异常的分子相关性，进行了空间转录组学分析，将全脑基因表达与耦合改变相关联。以艾伦人脑图谱中15,631个基因的表达水平为预测因子，以123个脑区中结构-功能耦合差异的病例-对照t值为响应变量，实施偏最小二乘回归。

第一个偏最小二乘分量被确定为主要空间协变模式。通过保留空间自相关结构的置换检验（10,000次迭代）评估统计显著性，确定观察到的区域统计t图与PLS1基因权重之间的协方差是否超过偶然预期。

2.8.3 抑郁症相关基因的分析

为探究抑郁症相关基因在结构-功能耦合异常中的作用，此研究采用了一组先前验证的24个重型抑郁症基因，该基因集来自AHBA的“皮层中通过原位杂交表征的1000个基因”。这些抑郁症相关基因包括ADRA2A、AVPR1B、CHRM2、CNR1、CREB1、CRH、CRHR1、CRHR2、CUX2、GAD2、GPR50、HTR1A、HTR1B、HTR1D、HTR3A、HTR5A、MAOA、PDE1A、SLC6A2、SLC6A4、SST、TAC1、TPH1和TPH2，涵盖抑郁病理生理中涉及的血清素能、多巴胺能和应激反应系统。

为探究这些抑郁症相关基因的表达模式与观察到的结构-功能耦合改变之间的空间关系。在识别出MDD基因集与处理后的AHBA转录组数据集之间的重叠基因后，进行了空间相关分析，以量化每个基因的区域表达与结构-功能耦合的病例-对照差异之间的关系。所有p值使用FDR方法进行多重比较校正，显著性阈值设为p < 0.05。

2.9 PLS1正负权重基因的功能富集分析

为阐明具有显著PLS1权重的基因的生物学基础，使用在线生物信息学平台进行综合功能富集分析。分别分析PLS1正权重(PLS1+)和负权重(PLS1-)的基因，以识别与青少年抑郁症中结构-功能耦合异常相关的独特功能特征。

在三个领域进行基因本体论富集分析(GO)：分子功能(MF)、生物过程(BP)和细胞组分(CC)。此外，进行京都基因与基因组百科全书通路富集分析(KEGG)以识别相关的生物学通路。所有富集分析采用超几何检验，并进行FDR多重比较校正。pFDR < 0.05的富集项被认为具有统计学显著性。

2.10 青少年抑郁症特异性基因特征的细胞类型定位与发育轨迹分析

2.10.1 区域结构-功能耦合改变的细胞表征

为描绘结构-功能耦合异常的细胞基础，将PLS分析识别的基因表达模式映射到特定皮层细胞群。利用5项基于死后皮层组织单细胞转录组研究的综合细胞类型参考数据集，该数据集提供了7种主要脑细胞类型的基因表达特征：内皮细胞、星形胶质细胞、小胶质细胞、少突胶质细胞、少突胶质细胞前体细胞(OPCs),、兴奋性神经元和抑制性神经元。

使用通过GSVA R包实现的单样本基因集富集分析(ssGSEA)，基于细胞类型特异性标志基因的表达谱量化每个脑区的细胞类型富集分数。采用此方法，将PLS1正权重与负权重基因列表与细胞类型特异性参考数据集进行系统比较，以识别结构-功能耦合改变的主要细胞。细胞类型富集的统计显著性使用置换检验结合FDR校正进行评估，显著性阈值设为p < 0.05。

2.10.2 细胞类型特异性基因特征的功能通路分析

为理解与结构-功能耦合改变相关的细胞类型特异性基因表达模式的潜在生物学机制，进行了综合功能富集分析。使用在线生物信息学平台，对细胞类型特异性特征与PLS1加权基因之间的重叠基因进行GO和KEGG通路分析。pFDR < 0.05的富集项和通路被认为具有统计学显著性。

2.10.3 青少年抑郁症相关基因的发育表达映射

为表征结构-功能耦合改变相关基因的发育轨迹，使用细胞类型特异性表达分析(CSEA)工具将PLS1正负基因列表映射到BrainSpan发育转录组图谱。该分析为每个基因集生成跨越15个发育阶段和12个脑区的富集分数。通过整合来自BrainSpan数据集的人类数据，检查了抑郁症相关基因在整个神经发育过程中的时空表达模式。

在脑海科技云平台中，内置了结构-功能耦合分析与神经影像-基因关联分析模块，支持用户批量处理多模态数据，并确保每一步参数设置都有据可查。此外，平台的项目管理模块可清晰记录数据筛选标准、排除被试原因、分析版本、质量控制等信息，极大提升了多模态整合研究的透明度和可复现性。感兴趣可联系预约产品演示。

3.结果

3.1 人口学与临床特征

青少年重型抑郁症患者与健康对照组的人口学和临床特征呈现在表1中。两组在年龄(p = 0.735)、教育年限(p = 0.522)或性别分布(p = 0.198)方面无显著差异。

青少年抑郁症组在临床测量中表现出显著较高的评分，包括HAMD-17(p < 0.001)、HAMA-14(p < 0.001)和BPRS(p < 0.001)。青少年抑郁症组在MCCB总分上表现出显著较低的认知表现(p < 0.001)。自杀意念和自杀企图在青少年抑郁症参与者中的发生率分别为62.2%和42.06%。

关于CERQ，与健康对照组相比，青少年抑郁症参与者在自责(p = 0.001)、沉思(p < 0.001)和责怪他人(p = 0.008)上表现出显著较高的评分。相反，他们在积极重评上得分显著较低(p = 0.005)。青少年抑郁症组在ATQ测量的自动负性思维上也表现出显著较高的评分(p < 0.001)。

表1所有受试者的人口统计学与临床特征。

3.2 区域结构-功能耦合异常

青少年抑郁症和健康对照组在全脑区域中表现出不同的结构-功能耦合值。在健康对照组中，区域耦合值范围为-0.670至0.682（图2A），显示某些脑区具有强结构-功能对齐而其他区域则解耦合。尽管青少年抑郁症组有相似的变异范围（-0.667至0.621；图2B），但其最大正相关低于健康对照组(0.621 vs. 0.682)，表明抑郁症中结构-功能同步峰值的一致性降低。

区域分析揭示了多个脑区和网络中显著的结构-功能耦合异常（表2，图2C）。在默认模式网络中，观察到双侧颞中回区域的耦合显著降低。相反，在躯体运动网络内的左侧中央后回和视觉网络内的左侧外侧枕叶皮层观察到耦合增加。

皮层下区域也表现出双向耦合改变，右侧海马在相邻亚区中同时显示耦合增加和减少。基底节显示耦合显著降低，而腹侧注意网络内的右侧扣带回也显示耦合降低。

网络水平分析进一步证实了这些区域模式，DMN中结构-功能耦合显著降低，SMN和VN中耦合显著增加（图2D）。

图2病例组与对照组间区域SC-FC耦合差异显示大脑各区域存在显著的空间变异性。

表2与健康对照组(HC)相比，结构-功能连接耦合改变在AMDD中的区域差异。

3.3 机器学习分类性能

XGBoost(AUC=0.83)和SVM(AUC=0.85)均能稳健区分青少年抑郁症患者与健康对照，具体敏感性、特异性见表3。特征重要性分析显示，DMN内的区域对分类性能贡献最大。包括准确性、敏感性和特异性在内的综合评价指标呈现在表3中。

还探究了区域结构-功能耦合改变预测青少年抑郁症症状严重程度的潜力。随机森林回归模型对临床症状显示出显著预测能力，在预测HAMD-17评分方面呈显著负相关(r = -0.49，p = 0.01)，但在预测其他症状方面未观察到统计学显著差异。

表3基于人类脑网组图谱评估XGBoost与 SVM 方法的诊断性能。

3.4 结构-功能耦合改变与神经递质系统的关联

空间相关分析揭示了青少年抑郁症中结构-功能耦合异常与多种神经递质系统之间的显著关联。其中与多巴胺能标志物的关联性最为显著，表明多巴胺信号在观察到的连接紊乱中起核心作用（图3）。

D2受体系统与结构-功能耦合改变的相关性最强(r = 0.363，pFDR = 0.020)，紧随其后的是多巴胺转运体密度(r = 0.360，pFDR = 0.020)。在全面检查多巴胺能系统时，观察到多巴胺合成能力标志物FDOPA也与耦合异常表现出强正空间相关(r = 0.348，pFDR = 0.020)。多巴胺D1受体系统同样显示出与耦合紊乱的显著关联，尽管较弱(r = 0.293，pFDR = 0.042)。

血清素能信号也表现出与结构-功能耦合异常的关联。具体而言，SERT密度与区域耦合改变呈显著正相关(r = 0.330，pFDR = 0.034)。耦合异常的神经化学特征也延伸至胆碱能神经传递，VAcHt密度与结构-功能耦合紊乱呈正空间相关(r = 0.293，pFDR = 0.042)。值得注意的是，GABAa受体的空间相关未通过多重比较校正(pFDR > 0.05)，这与强调GABA功能障碍的成人抑郁症研究形成对比。

图3 AMDD 中异常区域SC-FC耦合与特定神经递质谱之间的关联性。

3.5 结构-功能耦合改变的转录组学相关性

PLS1分析揭示了结构-功能耦合的病例-对照差异与基因表达模式之间存在显著空间相关（r = 0.345，pperm = 0.0001；图4A-B）。PLS1解释了转录谱与结构-功能耦合改变之间关系总方差的14.51%。

进一步检查了PLS1基因集与先前识别的MDD相关基因之间的重叠。在24个已知MDD相关基因中，13个与PLS1基因集重叠。还计算了单个基因表达模式与结构-功能耦合病例-对照t值之间的空间相关性，随后使用自旋置换检验进行空间自相关校正。经FDR校正后，两个基因在青少年抑郁症中显示出与结构-功能耦合改变的显著空间相关：HTR5A(r = -0.26，pspin = 0.003，pFDR = 0.022)和MAOA（r = 0.29，pspin = 0.0009，pFDR = 0.014；图4C-D）。

图4基因表达与区域SC-FC耦合差异之间的关系。

3.6 功能富集分析

使用单变量Z检验（|Z| ≥ 3.65），识别出1631个PLS1+基因和2062个PLS1-基因（所有pFDR < 0.0001）。PLS1-基因在与线粒体功能和神经退行性过程相关的通路中富集。最显著的KEGG通路富集包括氧化磷酸化、产热以及多种神经退行性疾病。GO分析显示与细胞呼吸、线粒体基因表达和组蛋白修饰相关的生物过程显著富集。细胞组分术语主要与线粒体结构相关，而分子功能类别包括各种转运体和转移酶活性（图5B、D、F）。

PLS1+基因在独特的生物学通路中富集。显著的KEGG通路包括氨基酸代谢、PPAR信号和内质网蛋白加工。GO分析识别出与脂肪酸代谢和应激反应以及有机酸生物合成相关的生物过程富集。细胞组分术语主要与膜结构和囊泡相关，而分子功能包括各种结合活性和氧化还原酶活性（图5A、C、E）。

图5通过偏最小二乘法分析鉴定的基因集功能特征分析

3.7 细胞类型特异性分析

3.7.1 结构-功能耦合改变背后的细胞类型特异性分子特征

细胞类型定位分析揭示了与青少年抑郁症中结构-功能耦合改变相关的独特分子特征（图6A）。PLS1+基因在星形胶质细胞、小胶质细胞和少突胶质细胞前体细胞中显著富集（图6B）。这些基因的功能分析显示星形胶质细胞富集于应激反应通路、嘌呤能信号和代谢过程，而胶质细胞群参与炎症和细胞外基质修饰通路（图6C）。

PLS1-基因在兴奋性和抑制性神经元中显著富集（图6D）。这些神经元基因主要与突触功能、神经调节和激素相关信号通路相关（图6E）。

3.7.2 亚型特异性基因表达的发育轨迹分析

时空分析揭示了PLS1+和PLS1-基因在不同脑区中的独特发育特征（图6F-G）。PLS1-基因在青春期和成年早期表现出皮层和小脑的富集，同时在儿童中期至晚期有次级信号。

相反，PLS1+基因表现出分布的皮层下表达模式，在以下时期显著富集：（1）儿童中期至晚期和胎儿早期的丘脑；（2）婴儿晚期和儿童中期至晚期的纹状体；（3）婴儿晚期的海马；（4）胎儿晚期和胎儿早期至中期的杏仁核。

图6. PLS加权基因在特定细胞类型、外周炎症及皮质层中的转录特征评估。

4.讨论

此多模态研究识别了青少年抑郁症中的三个核心神经生物学特征：（1）网络特异性结构-功能耦合异常，表现为默认模式网络中耦合降低和躯体运动网络中耦合升高；（2）区别于成人抑郁症的多巴胺中心性神经化学特征；（3）涉及胶质-神经元交互作用的细胞类型特异性转录组特征。这些发现突显了青少年抑郁症如何源于网络整合、神经化学信号和分子通路的发育性紊乱，这些与青春期的神经发育背景相一致。

4.1 青少年抑郁症中的网络特异性结构-功能耦合异常

此研究观察到DMN内结构-功能耦合降低与感觉运动区域耦合升高相结合的模式，代表了青少年抑郁症中脑网络架构的根本重组。虽然DMN解耦合与成人抑郁症的发现相似，但SMN的超耦合似乎更特异于青少年表型。这种SMN超耦合可能反映了发育特异性过程，因为这些网络在青春期经历实质性精炼。此外，这些网络改变与青少年抑郁症队列中观察到的临床现象直接对应。DMN解耦合影响与自我参照加工和适应不良沉思相关的区域，与参与者CERQ上显著升高的沉思和自责评分一致。这种自我相关信息整合受损可能导致青少年抑郁症特征性的负性自我评价和认知僵化。同时，感觉运动网络中的超耦合可能代表一种试图通过增强身体意识来调节压倒性负性情绪的补偿机制。与认知重评相比，这可能构成一种较不有效的情感调节策略，认知重评在青少年抑郁症组中显著降低。

4.2 结构-功能耦合异常的神经化学相关性

结构-功能耦合异常的神经化学相关性揭示了青少年抑郁症中以多巴胺能系统参与为特征的独特特征。与成人抑郁症中血清素能和GABA能主导相比，研究结果指向青少年中多巴胺中心性模型。D2受体系统和多巴胺转运体显示出与耦合改变的最强空间相关。这与强调多巴胺在青春期关键重组（影响奖赏和情绪调节）的发育模型一致。与血清素能和胆碱能系统的较弱关联表明青少年抑郁症中神经递质参与的层级组织。这表明抑郁症起病的发育时间塑造了其神经化学特征，可能将青少年抑郁症定义为一个独特的神经生物学亚型。

4.3 结构-功能耦合改变背后的分子特征

转录组学分析识别了与结构-功能耦合改变空间相关的两个核心分子谱。与耦合差异正相关的基因在免疫和炎症反应以及脂肪酸代谢过程中显著富集，并映射到PPAR信号和内质网蛋白加工等KEGG通路。具体而言，在耦合增加区域内，PLS1+基因显示出对神经递质合成至关重要的氨基酸代谢通路和对神经炎症调节重要的PPAR信号的富集。

相反，与耦合负相关的基因在线粒体和能量代谢功能（包括氧化磷酸化和有氧呼吸）中强烈富集，并与神经退行性通路相关。这与DMN解耦合的观察一致，表明代谢功能障碍可能损害这一自我参照网络的整合能力。此外，PLS1-基因中组蛋白修饰通路的富集表明表观遗传机制也可能导致网络功能障碍。

综合而言，这种转录组二分法表明，感觉运动网络中的超耦合可能由胶质细胞介导的神经免疫和代谢过程驱动，而高阶网络（如DMN）中的解耦合可能源于神经元能量代谢缺陷。这一解释通过已知MDD风险基因与耦合图的空间相关性得到进一步强化：与HTR5A表达的负相关表明，血清素能抑制控制降低可能导致青少年抑郁症的情绪调节缺陷，而与MAOA的正相关表明，单胺分解增加可能是神经化学分析中观察到的多巴胺能功能障碍的基础。

4.4 网络功能障碍的细胞类型特异性贡献

细胞类型富集分析揭示了一个引人注目的胶质-神经元二分法，为结构-功能耦合改变提供了机制性见解。PLS1+基因在胶质细胞中显著富集，特别是星形胶质细胞和小胶质细胞。这表明胶质细胞驱动的过程，如神经炎症信号和脂质代谢，可能促进SMN中观察到的超耦合。星形胶质细胞调节神经血管耦合和突触可塑性，而小胶质细胞通过修剪和细胞因子释放塑造神经回路。这些基因在应激反应和代谢通路中的同时富集将胶质细胞定位为整合环境应激与青少年抑郁症网络水平变化的关键角色。

相比之下，PLS1-基因主要在神经元群体中表达，并从根本上与氧化磷酸化和线粒体功能相关。这表明跨模态联合网络（如DMN）中的解耦合源于神经元生物能量能力的缺陷。鉴于其高能量需求，神经元特别容易受到线粒体功能障碍的影响，这可能损害突触传递和网络同步。因此，青少年抑郁症中的结构-功能解耦合可能反映了能量受损的神经元无法维持DMN所需的高成本功能整合。

4.5 发育轨迹与关键窗口

发育轨迹分析揭示了遗传脆弱性与环境因素相互作用以塑造青少年抑郁症风险的关键窗口。PLS1-基因在青春期和成年早期皮层和小脑中的富集指向这些区域在发育期的特定脆弱性——恰好在抑郁症通常出现的时期。相反，PLS1+基因在多个发育阶段（从胎儿早期到儿童中期至晚期）的分布皮层下表达表明，皮层下脆弱性可能在发育早期就已存在。这种时间特异性对干预时机具有重要影响，表明虽然某些风险因素可能在产前就已建立，但青春期代表了皮层网络整合变得易受破坏的第二个关键窗口。这些发现支持青少年发育的双重打击模型，即早期皮层下脆弱性与青春期快速神经重组期间后期皮层网络脆弱性相互作用。

4.6 临床转化与未来方向

研究结果强调了发育针对性干预的必要性。多巴胺能相关性的显著地位表明，针对该系统的药理学方法（如多巴胺受体调节剂）可能对青少年特别有前景。此外，细胞类型特异性转录组特征提示精准治疗的潜力。例如，星形胶质细胞基因在甘氨酸代谢通路中的富集可为针对胶质细胞介导的代谢支持或谷氨酸能调节的策略提供信息。类似地，靶向神经炎症的药物可能解决观察到的胶质细胞驱动的病理生理学。

除药理学外，研究结果支持神经调控干预。聚焦于加强DMN整合或正常化SMN活动的经颅磁刺激方案可能直接针对识别出的耦合异常。心理治疗可被增强以靶向网络特异性功能障碍：认知行为疗法针对DMN相关的沉思，基于正念的方法针对SMN相关的躯体过度激活。基因表达的时间特异性进一步表明，儿童期的预防性干预可能有助于减轻早期皮层下脆弱性。

此研究存在一定局限性。横断面设计无法推断发育轨迹的因果联系。此外，纳入任务态功能磁共振成像将允许直接评估结构-功能耦合改变如何调节情绪调节。最后，使用成人死后转录组数据是一个公认的局限；未来使用青少年特异性转录组数据的研究将极具价值。

5.结论

综上所述，此研究的多层次分析表明，青少年抑郁症涉及脑网络组织的发育特异性破坏，特征为DMN解耦合和SMN超耦合。这些改变以独特的以多巴胺为中心的神经化学特征为基础，反映了互补的转录组和细胞特征：胶质细胞主导的免疫代谢过程可能驱动网络超耦合，而神经元能量代谢缺陷可能是网络解耦合的基础。识别出这种胶质-神经元二分法为理解青少年抑郁症病理生理学提供了一个机制框架，并为发育针对性干预提供了多重靶点。