《NeuroImage》：From Cell to Circuit: Altered Functional Network Topology in iPSC-derived Neuronal Circuits in Major Depressive Disorder
编辑推荐：
患者来源的诱导神经元在细胞水平上模拟了脑网络的组织结构。研究人员利用钙成像技术记录了人类神经网络中的神经元活动。网络分析揭示了抑郁症中网络分离性和聚类性的降低。在抑郁症患者的网络中观察到功能连接减少。抑郁症中的信息处理损伤在微观尺度上即可被观察到。
患者来源的诱导神经元在细胞水平上模拟了脑网络的组织结构。研究人员利用钙成像技术记录了人类神经网络中的神经元活动。网络分析揭示了抑郁症中网络分离性和聚类性的降低。在抑郁症患者的网络中观察到功能连接减少。抑郁症中的信息处理损伤在微观尺度上即可被观察到。
论文解读
一、研究背景与科学问题
大脑是一个高度复杂的网络系统，具有独特的功能拓扑结构，能够在远距离脑区之间实现高效的信息整合，同时保持一定程度的分离性以支持局部神经通信和功能特化。这种整合与分离过程之间的平衡对于脑功能至关重要，它使大脑能够满足高度动态的认知需求并产生复杂、适应性的行为。近年来，人类脑网络拓扑结构的改变已被用于解释多种精神疾病，包括重度抑郁症（MDD）。
重度抑郁症是一种影响全球约2.8亿人的精神病理状态，其特征是持续的情绪低落和认知功能紊乱。尽管抑郁症十分普遍，但其确切的病理生理机制仍然难以捉摸。目前关于抑郁症功能网络拓扑改变的证据大多来源于宏观尺度的研究，即使用功能磁共振成像（fMRI）进行的全脑成像研究。然而，在细胞水平上，即抑郁症背景下神经网络的微观功能组织，在很大程度上仍未被探索。
图论为研究复杂网络拓扑提供了数学框架，其中神经元网络被建模为由节点和边组成的图。聚类系数衡量网络中相邻且功能相似的节点形成密集簇的倾向，反映网络的分离性；全局效率则量化信息在整个网络中并行传输的效率，反映网络的整合水平。既往宏观尺度研究在抑郁症网络拓扑改变方面结果并不一致，部分研究报告了全局效率增强，另一些则显示网络整合和聚类性降低，还有研究未发现显著拓扑差异。这些不一致被归因于样本异质性和数据预处理方法的差异。值得注意的是，目前尚缺乏针对抑郁症微观尺度网络拓扑的图论研究。
人类诱导多能干细胞（hiPSCs）技术的出现为在细胞水平上进行疾病建模和药物发现开辟了新途径。该技术可将人类体细胞重编程为多能干细胞，再分化为各种神经元命运，同时保留供体细胞的遗传组成。在抑郁症研究中，hiPSCs已被用于生成血清素能神经元以模拟对抗抑郁药物SSRIs的耐药性，并探索多巴胺能神经元中氯胺酮依赖的神经可塑性。然而，hiPSC来源神经元回路的功能网络架构及其与宏观水平改变的比较，在MDD背景下仍 largely 未被探索。
本研究直接针对这一知识空白，通过比较抑郁症患者hiPSCs来源的微观网络与匹配对照组的网络在功能组织属性上的差异，旨在提供与MDD相关的细胞水平网络改变的深入见解。
二、研究设计与主要技术方法
研究人员从德国雷根斯堡大学医院皮肤科获取了9名MDD患者和9名年龄、性别匹配的健康对照者的皮肤活检样本。MDD诊断基于ICD-10标准，抑郁严重程度使用21项汉密尔顿抑郁量表（HAMD21）评估。研究人员通过附加体质粒载体将人成纤维细胞重编程为hiPSCs，随后分化为神经祖细胞（NPCs），再诱导为成熟神经元（iNeurons），在21天体外培养（DIV）时进行成像。
主要技术方法包括：使用ratiometric钙指示剂Fura-2/AM进行活细胞钙成像，记录自发神经元活动的20分钟视频；采用非刚性运动校正（NoRMCORRE）进行运动伪影校正；通过相关图像和CellProfiler建立半自动分割流程识别活跃神经元；提取每个分割区域的时间序列；计算成对Pearson相关系数构建功能连接（FC）矩阵；应用多个绝对阈值（rthresh范围±0.2-0.8，步长0.1）对FC矩阵进行二值化；使用NetworkX包构建二值无向图；计算聚类系数、全局效率和平均节点度等图论指标；采用三因素重复测量方差分析（ANOVA）和蒙特卡洛置换检验进行统计比较。
三、研究结果
3.1 人口学特征与样本特征
研究人员检查了年龄、性别、体重指数（BMI）和HAMD评分的组间差异。对照组和MDD组在性别（每组各4男5女）和年龄（对照组33.22±3.44岁，MDD组32.33±3.82岁，t(16)=-0.19，p=0.85）方面无显著差异。BMI同样无组间差异（对照组22.67±0.72，MDD组24.58±1.29，t(16)=1.29，p=0.22）。然而，两组HAMD评分存在显著差异（对照组0.4±0.25，MDD组26±1.63，t(16)=12.19，p<0.001），符合预期。
3.2 神经元标记物表达
iPS来源神经元表达了典型的神经元标记物β-III微管蛋白（B3TUB）和微管相关蛋白2（MAP2），并对特异性突触蛋白囊泡谷氨酸转运体1（VGLUT1）呈免疫阳性，表明NPCs成功分化为具有谷氨酸介导信号传导功能的功能性神经元。
3.3 基本图属性
研究人员在二值化网络进行图分析前计算了绝对相关值的均值和方差。未观察到这些描述性统计量的组间差异（平均相关：F(1,8)=0.44，p=0.53；相关方差：F(1,8)=0.11，p=0.75），表明组间差异并非由整体连接强度或变异性的变化驱动。
3.4 分离性属性
研究人员通过计算聚类系数评估培养神经元网络的分离性属性。结果显示，患者来源神经元培养物的聚类系数较对照组显著降低（Group主效应：F(1,8)=5.97，p=0.04，ges=0.31），表明MDD患者体外网络的分离能力下降。该结果经5000次置换检验确认（pperm=0.04）。聚类系数随网络规模增大而增加（Network size主效应：F(37,296)=73.59，p<0.001），随rthresh降低（即更宽松）而增加（Threshold主效应：F(2.42,19.34)=1003.25，p<0.001）。事后t检验显示，大多数聚类系数组间差异出现在较大网络规模下。
3.5 整合性属性
研究人员评估了全局效率以衡量网络的整合属性。数据显示MDD患者来源神经元网络的全局效率较对照组呈总体降低趋势，但未达到统计学显著性（无Group主效应：F(1,8)=3.92，p=0.08，ges=0.23，pperm=0.1）。然而，检测到显著的Group×Network size×Threshold三向交互效应（F(222,1776)=1.57，p<0.001，ges=0.01）。全局效率随网络规模增大而显著增加（Network size主效应：F(37,296)=23.56，p<0.001），随rthresh增大而显著降低（Threshold主效应：F(1.24,9.89)=1157.36，p<0.001）。
3.6 整体网络连接性
平均节点度是图论分析中最基本的度量之一，代表网络中节点平均拥有的功能连接数。研究人员发现MDD患者来源神经元培养物的平均节点度较对照组显著降低（Group主效应：F(1,8)=5.44，p=0.048，ges=0.26），经置换检验确认（pperm=0.035）。平均节点度随网络规模增大而显著增加（Network size主效应：F(37,296)=167.82，p<0.001），在更宽松阈值下也显著增加（Threshold主效应：F(6,48)=622.24，p<0.001）。事后t检验显示，组间差异在较大网络规模（超过20个细胞）时更为明显。
3.7 零模型比较
研究人员将观察网络与保持网络规模和度分布的随机化零模型网络进行比较。对照组和MDD患者网络的零模型之间在聚类系数（无Group主效应：F(1,8)=4.58，p=0.06）和全局效率（无Group主效应：F(1,8)=3.92，p=0.08）方面均无显著差异。组内比较确认，对照组（network type主效应：F(1,4)=96.57，p<0.001）和患者组（F(1,4)=111.58，p<0.001）的聚类系数均显著高于其大小和度匹配的随机化对应网络。两组观察网络的全局效率均显著低于随机化网络。这些结果证实两组均表现出结构化、非随机的网络组织，且观察到的聚类系数差异反映真实的拓扑改变而非度分布差异的伪影。
四、讨论与结论
研究人员在讨论中指出，患者来源神经元培养物中聚类系数的显著降低标志着抑郁症细胞水平上局部组织改变和网络分离性及功能特化能力下降。聚类拓扑对神经功能具有重要影响，其破坏可能对网络动力学产生显著影响，进而影响神经活动，并可能促成既往在MDD宏观尺度上报道的动态储备受限。此外，分离性拓扑对于网络抵抗扰动的鲁棒性和韧性至关重要，缺乏强局部互连性的网络可能无法限制和阻止病理过程在系统内的传播。
平均节点度的降低与既往尸检研究一致，这些研究提示MDD症状与前额叶皮层和海马突触密度降低相关。研究人员推测观察到的节点度降低可能反映突触密度减少，阻碍了神经元间连接的形成。这一假设与既往体外和体内研究一致，这些研究强调了突触丢失对MDD网络连接性和动力学的影响。
关于全局效率未观察到显著组间差异，研究人员认为尽管分离性测量发生显著改变，细胞水平上全局整合能力可能得以维持，整合缺陷可能仅在更高组织水平上通过多个细胞群体或脑区之间的相互作用而显现。
研究人员承认研究存在若干局限性：对照组非抑郁状态的分类仅依赖自我报告；培养物为纯神经元单培养，缺乏胶质细胞等支持性细胞类型；样本量较小限制了更大网络的分析。此外，所有MDD患者在活检时均接受过抗抑郁治疗，无法确定观察到的网络差异反映稳定的特质脆弱性、短暂的状态/治疗相关效应，还是两者的组合。
总之，这是首个使用iPSCs研究抑郁症功能网络拓扑的研究。研究结果支持将抑郁症视为涉及异常信息传递和功能拓扑紊乱的网络障碍的观点。MDD患者来源神经网络的聚类系数和节点度降低与多项宏观神经影像学研究报告的类似减少一致。研究结果表明，与MDD相关的拓扑改变在基本细胞水平即可被检测到，但尚不清楚这些细胞改变是否会持续并扩展影响完全发育的大规模脑网络。未来研究应探索跨空间尺度的微-宏观网络拓扑关联，并利用患者来源iPSC模型进行个性化功能药物筛选。
