《Frontiers in Pharmacology》：Selective behavioral and prefrontal transcriptomic correlates of acute clozapine exposure in socially isolated mice
编辑推荐：
摘要专业翻译
背景： 氯氮平（clozapine）是治疗难治性精神分裂症的有效抗精神病药物，但其急性分子应答尚未完全阐明。社交隔离饲养可产生精神分裂症相关的行为异常，并允许同时评估饲养条件与药物暴露的影响。
目的： 表征急性氯氮平刺激在社交隔离小鼠中的行为与前
摘要专业翻译
背景： 氯氮平（clozapine）是治疗难治性精神分裂症的有效抗精神病药物，但其急性分子应答尚未完全阐明。社交隔离饲养可产生精神分裂症相关的行为异常，并允许同时评估饲养条件与药物暴露的影响。
目的： 表征急性氯氮平刺激在社交隔离小鼠中的行为与前额叶皮质（prefrontal cortex, PFC）转录组相关性，重新评估先前优先关注的Ank3与Nrxn1，并为候选基因提供外部细胞情境假说。
方法： 研究人员采用四组饲养×药物设计，对雄性C57BL/6J小鼠进行旷场实验（open-field test）、新物体识别实验（novel-object recognition test）和前脉冲抑制实验（prepulse inhibition test, PPI）。为实现可重复的全转录组审计，研究人员利用包含错误发现率（false discovery rate, FDR）控制的因子DESeq2模型对11个公开可获取的前额叶皮质文库进行了重新定量。历史配对优先化、网络分析、机器学习筛选、分子对接和分子动力学模拟均被视为探索性分析。成年野生型前额叶皮质单细胞RNA测序数据用于细胞定位、CellChat分析和虚拟敲除，而胚胎皮质参考数据则单独用于CellOracle转录因子扰动分析。
结果： 急性氯氮平最显著的行为效应体现在对隔离相关高运动活性的抑制上，而识别记忆未获改善，前脉冲抑制结果因刺激强度而异。初步DESeq2分析在FDR <0.05水平鉴定出56个饲养×药物互作基因；其中Lamb2是FDR支持最强的转录组信号。Nrxn1在历史配对对比中表现出方向相反的改变，但在因子模型中未达到FDR显著性，而Ank3未显示相应的因子模型支持。在成年前额叶皮质单细胞数据中，Lamb2富集于内皮细胞，而Ank3与Nrxn1广泛分布于神经元群体。虚拟Ank3与Nrxn1敲除分别鉴定出16个和15个显著的非靶基因，其中有5个共享基因，且未发现显著的相互候选效应。胚胎CellOracle分析识别出Ank3与Nrxn1依赖于模型的上游转录因子敏感性。
结论： 这些发现界定了急性氯氮平刺激的选择性行为与前额叶转录组相关性，而非广义的行为挽救或临床治疗应答预测。Lamb2是初步模型中FDR支持最强的信号，而Ank3与Nrxn1仍属探索性计算候选。外部的单细胞与虚拟扰动分析提供了可检验的细胞和调控假说，而非直接结合、因果调控或氯氮平诱导效应的体内证据。
论文解读文章
一、研究背景与问题
精神分裂症是一种严重的精神疾病，表现为阳性、阴性和认知症状。氯氮平（clozapine）是目前对许多治疗抵抗性精神分裂症患者最有效的抗精神病药物，但其临床相关作用发生在持续治疗过程中，不能简化为单次给药后的转录反应。然而，急性实验范式有助于分离即时药物相关效应与饲养相关的分子差异，前提是不将其描述为临床应答预测模型。
研究氯氮平应答的机制面临一个核心设计问题：药物暴露后的分子差异本身无法区分治疗性应答、急性药理学效应、疾病状态或非特异性运动效应。社交隔离饲养（social isolation rearing, SIR）可产生可重复的精神分裂症相关行为异常，而前额叶皮质（prefrontal cortex, PFC）是研究精神分裂症相关行为与转录表型的相关组织。因此，该模型能够识别急性的、依赖饲养条件的药物相关性特征，但若无应答者异质性与独立验证，则无法建立临床应答生物标志物。
本研究旨在将行为测试与PFC RNA测序相结合，检测社交隔离小鼠中急性氯氮平暴露的即时相关性。主要问题是行为效应是否依赖于饲养条件，以及哪些转录组信号表现出饲养×药物互作或氯氮平后的反向变化。网络分析、特征选择、结构模拟、外部单细胞重注释、细胞通讯分析与虚拟扰动均被视为探索性假说生成层，而非独立验证。
二、研究方法概述
研究人员采用了四组饲养×药物设计（社会群居SOC vs. 社交隔离SIR，溶剂对照vs. 氯氮平处理），对雄性C57BL/6J小鼠依次进行旷场实验（OFT）、新物体识别实验（NOR）和前脉冲抑制实验（PPI），并在最后一次行为测试后约2.5–3小时收集PFC组织进行bulk RNA测序。关键分析方法包括：基于因子DESeq2模型的差异表达分析（使用11个公开可获取的PFC文库，来源于BioProject PRJNA1014397）；历史名义P值配对比较与机器学习的候选基因筛选（随机森林、SVM-RFE、LASSO、ElasticNet）；分子对接（AutoDock Vina）与分子动力学模拟（Amber 24）；外部成年野生型小鼠PFC单细胞RNA-seq数据（GSE181021）的细胞定位、CellChat细胞通讯分析与scTenifoldKnk虚拟敲除；以及胚胎皮质参考数据（GSE274763）的CellOracle上游转录因子扰动分析。
三、研究结果
3.1 行为学评估结果
旷场实验： 运动距离显示出显著的药物主效应（F(1,34)=31.166, P<0.001）和显著的饲养×药物交互作用（F(1,34)=9.898, P=0.003），而饲养主效应不显著。Tukey校正比较显示，SIR-溶剂组小鼠比SOC-溶剂组运动更远（P=0.005），氯氮平显著降低了SIR小鼠的运动距离（P<0.001）。结果表明急性氯氮平与SIR小鼠运动距离的显著减少相关，但可能同时反映了饲养依赖性效应和一般性药物相关运动抑制。
新物体识别实验： 辨别比率显示出显著的饲养主效应（F(1,34)=16.687, P<0.001），但药物主效应和交互作用均不显著。SIR-溶剂组小鼠的辨别比率低于SOC-溶剂组（P=0.001），SIR内部溶剂与氯氮平的比较不显著。结果表明存在饲养相关的识别记忆缺陷，但在此测试方案下无急性氯氮平挽救的证据。
前脉冲抑制实验： 在69 dB下观察到显著的饲养和药物主效应，但交互作用不显著；在78 dB下所有效应均不显著；在88 dB下仅有显著的饲养主效应，SIR-溶剂组小鼠的PPI低于SOC-溶剂组（P=0.038）。总体而言，PPI结果显示了选定强度下的饲养相关缺陷，但不支持一致的急性氯氮平效应。
行为学发现总结： 急性氯氮平对隔离相关高运动活性有最清晰的效果，而识别记忆未改善，PPI结果因刺激强度而异。行为数据支持选择性急性行为效应，而非社交隔离相关异常的全面逆转。
3.2 差异基因表达分析
在11个PFC文库中，每个样本处理了1910万至2270万个读段对，77.9%–90.4%成功伪比对到小鼠GRCm39转录组。主成分分析（PCA）显示PC1和PC2分别解释了29.2%和22.5%的变异。两个SIR-氯氮平样本与SOC-氯氮平文库聚类，而SIR2PFC在PC2上分离，SOC3PFC未与其他SOC-溶剂文库聚类。所有文库均保留在初步分析中。
在FDR<0.05水平下，溶剂条件下SIR vs. SOC对比包含124个基因，SOC内氯氮平vs.溶剂对比包含245个基因，SIR内氯氮平vs.溶剂对比包含18个基因。共有56个基因显示出显著的饲养×药物交互作用。其中，Lamb2显示出正向交互作用（log2 fold change=1.261, FDR=5.38×10-5），并且是唯一一个在三个预设DESeq2对比（溶剂条件下的隔离效应、SIR组内的氯氮平效应以及饲养×药物交互作用）中均显著的基因。然而，这一双对比方向性成分在排除SIR2PFC后减弱，表明其对小型不平衡的组结构敏感。
Nrxn1在溶剂条件下SIR小鼠中相对于SOC小鼠降低，而在SIR小鼠中氯氮平处理后升高，但在因子模型中未达到FDR显著性（交互作用FDR=0.485）。Ank3未显示显著的交互作用效应（FDR=0.931）。在edgeR准似然法和limma-voom敏感性分析中，无任何饲养×药物交互作用基因通过FDR校正。对56个FDR显著交互作用基因的富集分析未发现任何经过多重检验校正后显著的GO或KEGG条目。
3.3 历史286基因重叠集的探索性富集分析
历史配对差异表达基因（DEG）对比产生了286个基因的重叠。GO富集分析显示，这些基因在NMDA受体聚集、轴突延伸调节、突触囊泡聚集和突触传递等生物学过程中显著富集；细胞组分涉及细胞粘附复合物、突触前细胞骨架和突触后致密区；分子功能涉及信号素受体活性和转录共调节因子结合。KEGG分析突出了核糖体、胆碱能突触和GABA能突触通路的富集。
3.4 蛋白质-蛋白质相互作用网络
286个历史探索性重叠基因的蛋白互作网络包含了髓鞘相关基因、突触相关基因和精神分裂症相关基因。
3.5 Ank3与Nrxn1的探索性计算优先化
四种机器学习方法（随机森林、SVM-RFE、LASSO、ElasticNet）的交集优先化了Ank3和Nrxn1。由于特征选择和性能可视化使用了同一六样本矩阵，该分析不能提供分类或治疗应答预测的无偏估计，仅作为探索性特征优先化程序保留。在独立重新处理的公共文库中，Nrxn1在溶剂条件下SIR vs. SOC中降低且在SIR内氯氮平后方向性升高，但均未通过FDR校正；Ank3未显示任何显著的饲养、药物或交互作用效应。因此，因子RNA-seq模型未能验证任一候选基因。
3.6 描述性数据集内候选基因可视化
基因水平的ROC曲线和列线图展示了在该矩阵内的分离情况，但不估计外部预测、校准或临床效用。
3.7 探索性免疫特征评分
使用小鼠特异性免疫标记矩阵进行的探索性签名评分未发现Ank3或Nrxn1与估计签名之间的稳健关联。
3.8 候选氯氮平应答蛋白的分子对接分析
分子对接结果显示，氯氮平与代表Ank3和Nrxn1相关蛋白区域的结构模型之间存在预测的结合模式。排名最高的对接分数分别为-9.4 kcal/mol（Ank3）和-7.0 kcal/mol（Nrxn1）。
3.9 分子动力学模拟
氢键分析显示模拟过程中氢键数量波动在0-1范围内。RMSD分析表明Ank3-氯氮平复合物大致保持在约4 ?以内，而Nrxn1-氯氮平复合物在约6 ?附近波动。RMSF分析显示Ank3-氯氮平复合物的残基波动较低。SASA和回转半径分析描述了复合物的溶剂暴露和紧实度特征。MM/GBSA计算的结合自由能分别为-70.29 kcal/mol（Ank3-氯氮平）和-66.99 kcal/mol（Nrxn1-氯氮平）。自由能景观描述了模拟复合物的主导构象区域。
3.10 成年PFC细胞情境与细胞通讯分析
成年野生型小鼠PFC单细胞数据包含5303个细胞，涵盖兴奋性神经元、抑制性神经元、内皮细胞、星形胶质细胞、小胶质细胞、少突胶质细胞系细胞和少突胶质前体细胞。Lamb2表达集中于内皮细胞（884个细胞中263个阳性，29.8%），而Ank3和Nrxn1广泛分布于神经元群体。CellChat分析鉴定了以内皮细胞为中心的LAMININ信号通路，以及在神经元亚型间的Nrxn1相关接触通讯。
3.11 候选相关神经元程序的计算机模拟扰动
在成年PFC神经元参考中，虚拟Ank3敲除鉴定出16个FDR显著的非靶基因，虚拟Nrxn1敲除鉴定出15个，两者共享5个基因。相互候选效应在两个方向上均不显著（FDR=1.00）。在胚胎皮质投影神经元参考中，CellOracle分析显示Mef2c对Nrxn1有最清晰的正向模式效应，Tcf4显示最强逆模式效应；Bclaf1和Hivep2对Ank3显示正向模式效应。
四、讨论与结论总结
讨论部分总结： 本研究使用社交隔离饲养模型检测了急性氯氮平暴露在PFC中产生的精神分裂症相关行为表型和分子变化。行为学数据的四组因子分析解决了原始统计报告中的不一致之处。急性氯氮平产生了领域选择性行为谱，其最清晰的效应是运动距离上的饲养×药物交互作用，主要由SIR小鼠的显著减少驱动。识别记忆未获改善，PPI仅在选定强度下显示饲养相关差异。因子转录组再分析大幅收窄了转录组主张，56个基因在FDR<0.05水平显示交互作用，Lamb2是唯一满足双对比方向性标准的基因，但该结果对排除PCA异常样本SIR2PFC敏感。历史网络和机器学习工作流优先化的Ank3和Nrxn1在因子DESeq2模型中均未达到FDR阈值。外部单细胞分析为候选基因提供了细胞类型特异性假说。虚拟敲除未支持Ank3-Nrxn1间的定向调控关系。研究受限于急性给药方案、仅限雄性队列和仅限于PFC组织。回顾性qPCR验证因缺乏存档组织而不可行。
结论部分翻译： 在社交隔离小鼠中，急性氯氮平产生了选择性行为谱，伴随PFC转录组变化。对公开可获取文库的因子再分析鉴定出FDR支持的饲养×药物交互作用信号，其中Lamb2是最强的转录组候选。Nrxn1保留了方向一致但FDR不显著的支持，而Ank3缺乏因子模型支持。成年PFC单细胞分析将Lamb2主要定位于内皮细胞，将Ank3/Nrxn1定位于神经元群体，而虚拟扰动分析生成了外部调控假说。这些发现需要独立的分子和功能验证，并不能确立氯氮平应答的因果机制。
