《GLIA》：Single-Cell Integration of Mouse Models of Autism Reveals Conserved Transcriptomic Signatures of Excitatory-Inhibitory Imbalance and Glial Dysfunction
编辑推荐：
自闭症谱系障碍（autism spectrum disorder, ASD）的病因异质性使得识别核心分子机制变得复杂。本研究系统地整合了来自11种不同遗传和环境小鼠（Mus musculus）模型的单细胞RNA测序（single-cell RNA sequen
自闭症谱系障碍（autism spectrum disorder, ASD）的病因异质性使得识别核心分子机制变得复杂。本研究系统地整合了来自11种不同遗传和环境小鼠（Mus musculus）模型的单细胞RNA测序（single-cell RNA sequencing, scRNA-seq）数据，涵盖了超过300,000个细胞，涉及多个脑区和发育阶段。研究人员利用先进的计算协调方法（reciprocal PCA, RPCA），识别出超越模型特异性变异的稳健收敛转录组特征。研究结果揭示了神经元通讯通路的广泛改变，特别是涉及PTN、SLIT和PTPR通路的兴奋-抑制信号增强。此外，通过NicheNet进行的细胞间通讯分析优先筛选出与回路兴奋性改变相关的配体-受体对，为失调的细胞间交互作用提供了预测蓝图。至关重要的是，研究人员发现了显著的胶质细胞参与，其特征为星形胶质细胞功能抑制，表现为与水/钾稳态受损相关的Aqp4和Kcnj10下调，以及成熟少突胶质细胞代谢失调，表现为Olig家族上调并伴随氧化磷酸化增强。差异表达基因（differentially expressed genes, DEGs）与SFARI数据库所列的高置信度风险基因显著重叠，尤其是在兴奋性神经元亚型中。与人类死后前额叶和顶叶皮质的单核数据的比较表明，与突触蛋白合成受损（GRIN2B, CAMK2A）和少突胶质细胞介导的神经发育缺陷（NRXN, DLGAP）相关的保守特征。这项研究证实了E/I失衡和胶质功能障碍是ASD中普遍的细胞病理学特征，为以单细胞分辨率识别新的治疗靶点提供了一个稳健的大规模参考。
研究背景与目的
自闭症谱系障碍（ASD）是一种复杂的神经发育疾病，全球约1%的儿童受到影响。尽管全基因组关联研究（GWAS）已鉴定出如CHD7、SHANK3等高风险基因，但大多数病例源于多基因与环境因素的复杂相互作用，导致其生物学基础难以阐明。现有治疗手段未能针对根本机制。单细胞RNA测序（scRNA-seq）技术的发展为在细胞层面解析ASD的转录组变化提供了新机遇。然而，单一动物模型的研究常受限于特定实验条件，且方法学差异阻碍了跨研究的数据整合。因此，本研究的目的是通过整合多种ASD小鼠模型的scRNA-seq数据，利用先进的生物信息学算法消除批次效应并保留生物学信号，从而识别出超越模型特异性的、保守的ASD核心分子通路和细胞病理学特征。
研究方法概述
研究人员从11项独立研究中获取了涵盖多种遗传（如Shank3, Ube3a, Fmr1敲除等）和环境（丙戊酸VPA暴露）因素的小鼠ASD模型scRNA-seq数据，包括大脑皮层、海马、纹状体等多个脑区及胚胎期至成年等多个发育阶段，共计超过30万个细胞。关键技术方法包括：采用互作主成分分析（RPCA）进行数据整合与批次校正；通过层次聚类和经典标志基因进行细胞类型鉴定；运用Wilcoxon秩和检验和伪批量（pseudobulk）DESeq2框架进行差异表达基因（DEG）分析；利用CellChat和NicheNet推断细胞间通讯网络及配体-受体相互作用；并与SFARI数据库及人类死后脑组织单核RNA-seq数据（Gandal et al., 2022）进行比较验证。
研究结果
3.1 整合数据集的特征描述
研究人员成功整合了11项研究的数据，最终包含291,916个高质量单细胞。通过RPCA整合，UMAP可视化显示细胞按生物状态而非技术批次聚集。层次分类和标志基因表达（如Aqp4标记星形胶质细胞，Gad2标记抑制性神经元）精确鉴定了主要细胞类型。空间和时间分布分析表明，整合后数据保留了原始发育阶段和脑区的特异性。局部逆辛普森指数（iLISI）评估证实了良好的数据集混合效果。
3.2 ASD小鼠模型中细胞类型的转录组改变
对比ASD与对照组，研究发现大量DEGs。其中，编码甲状腺素转运蛋白的Ttr基因在多个细胞类型中显著上调；而Olig基因家族（Olig1, Olig2）在非神经元细胞（尤其是少突胶质细胞）中普遍上调。在星形胶质细胞中，研究人员发现其功能相关基因（如Aqp4, Kcnj10, Gja1, Glul）广泛下调，提示水/钾稳态、缝隙连接通讯和谷氨酸-谷氨酰胺循环等功能受损。在成熟少突胶质细胞中，多数DEGs上调，富集于氧化磷酸化通路。伪批量验证分析确认了这些转录组改变的稳健性。
3.3 整合模型与单个数据集的比较
整合数据集中的DEGs与各个单独研究子集的结果显示出高度相关性，尤其在抑制性和兴奋性神经元中，证实了跨模型保守的转录组失调。在VPA模型子集中，多个ASD相关基因（如Fmr1, Foxg1）表现出广泛的细胞类型下调。与Donnard et al.（2022）的独立数据集进行超几何富集检验，发现在星形胶质细胞、少突胶质细胞前体细胞等多种细胞类型中存在显著的方向一致性重叠，进一步验证了结果的可靠性。
3.4 ASD小鼠模型中细胞间通讯紊乱
CellChat分析显示，与对照组相比，ASD模型中抑制性与兴奋性神经元之间的相互作用强度显著增加，暗示了神经元信号的普遍失调。具体而言，涉及神经发育布线过程的多个关键信号通路，如PTN、SLIT、PTPR和SEMA3通路，在ASD条件下的相互作用数量增多。NicheNet分析进一步预测了由抑制性神经元（发送者）到兴奋性神经元（接收者）的关键配体-受体对（如Pdgfa, Adamts3, Reln）及其下游靶基因。
3.5 整合模型DEGs与临床相关基因的汇聚
将本研究DEGs与SFARI数据库比对，发现大量重叠基因，其中高置信度风险基因（评分1或2）占比较高，且在兴奋性神经元亚型（如L4/5 IT）和星形胶质细胞中富集显著。超几何检验证实了这种富集并非偶然。与人类死后皮质单核RNA-seq数据的跨物种比较显示，在多个细胞类型中存在显著的非随机方向一致的重叠。例如，在顶叶皮层的兴奋性神经元中，共享的DEGs在SynGO分析中富集于“突触前/后翻译”等生物学过程，涉及GRIN2B, CAMK2A等基因；在成熟少突胶质细胞中，共享的DEGs则富集于神经发育障碍相关通路，并形成一个包含DLGAP1, DLGAP2, NRXN1, NRXN3的功能模块。
讨论与结论
本研究通过对11种ASD小鼠模型的系统性单细胞转录组整合，首次构建了一个统一的ASD相关转录失调参考图谱。研究证实，兴奋-抑制（E/I）失衡和胶质细胞功能障碍并非模型特异的假象，而是跨越不同遗传和环境病因的、稳健且保守的核心分子特征。研究发现强调了星形胶质细胞功能抑制（如Aqp4, Kcnj10下调导致的水/钾稳态失衡）和成熟少突胶质细胞代谢失调（Olig家族上调及氧化磷酸化增强）在ASD病理中的重要性。与SFARI风险基因库及人类死后脑组织数据的显著交叉验证，进一步证明了这些发现具有临床转化意义。该研究为理解ASD的通用细胞病理学机制提供了高分辨率的参考，并为识别潜在的新型治疗靶点奠定了坚实基础。
