IF50.5｜自闭症干细胞模型揭示发育收敛与分化机制

大规模类器官转录组揭示不同ASD突变在发育早期共享核心调控网络。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

自闭症谱系障碍（ASD）遗传异质性极高，已发现上百个高风险基因，但死后脑组织转录组分析显示不同遗传背景的ASD患者存在共同的分子失调模式。然而，这些收敛性变化如何在发育过程中产生，以及不同罕见突变是否共享机制，仍不清楚。由于大多数ASD风险基因在胎儿期表达高峰，而人脑研究多在发育后期进行，早期神经发育的作用难以直接评估。干细胞来源的类器官模型为在体外模拟人脑早期发育提供了可能，但既往研究多局限于单一突变或小样本，缺乏跨多种遗传形式的系统比较。

科学问题

卡在哪里： 此前研究未能回答：不同ASD相关突变是否在发育过程中收敛到共同的分子通路？收敛发生在哪个阶段？驱动收敛的核心调控因子是什么？此外，特发性ASD是否与已知突变形式共享表型也缺乏证据。

本文要回答： 利用大规模hiPS细胞来源的皮层类器官，系统比较8种ASD相关突变和特发性ASD在多个发育时间点的转录组动态，鉴定共享的调控网络及其驱动因子，并验证其对下游基因表达的影响。

技术路线

作者主要用了：hiPS细胞、皮层类器官、RNA-seq、WGCNA、RcisTarget、CRISPRi、IP-MS。

作者从96个hiPS细胞系出发，经全基因组测序质控后保留70个系（55个个体），涵盖8种ASD相关突变、特发性ASD和对照。每个系分化为皮层类器官（hCO），在分化第25、50、75、100天收集样本进行RNA-seq。通过差异表达、共表达网络（WGCNA）、转录因子结合基序富集（RcisTarget）等分析，识别跨突变共享的调控模块，并利用CRISPRi筛选和免疫沉淀-质谱（IP-MS）验证候选调控因子的功能和蛋白互作。

图 1　Fig. 1。图1展示了实验设计和数据质控：a为从hiPS细胞到皮层类器官再到测序的流程，标注了各ASD形式和对照的细胞系数目；b示意16p11.2缺失和22q13.3缺失的类器官分化方案；c显示同一样本内和不同样本间的基因表达Spearman相关性，表明高重复性。

核心亮点

亮点 1｜早期突变特异性，后期跨形式收敛

差异表达分析显示，各突变形式在分化早期（第25天）变化最大，且不同突变间共享的差异基因较少。但随着分化推进，不同突变形式的log2FC相关性显著增强（第100天与第25天相比P=2.7×10^-5），元分析在第100天鉴定出最多显著基因（6,680个），表明转录组变化随发育收敛。

图 2　Fig. 2。图2呈现转录组关系：a为基于logFC相关性的层次聚类热图，显示不同ASD形式随分化时间的聚类模式；b为各形式在不同时间点的DEG数量，早期变化最大；c为第25天各形式间DEG重叠；d为重叠基因的GO富集；e显示DEG中ASD风险基因的富集；f-h显示跨形式logFC相关性随时间增强，元分析在第100天基因数最多。

亮点 2｜M5模块：ASD风险基因富集且跨突变下调

WGCNA共表达网络分析发现模块M5显著富集ASD风险基因（OR=3.6, FDR=9.4×10^-12），且在16p11.2缺失、16p11.2重复、PCDH19突变等多个形式中显著下调。M5与胎儿脑中期共表达模块重叠，功能富集于组蛋白修饰和转录因子结合，提示其可能作为上游调控枢纽。

图 3　Fig. 3。图3展示生物学过程和细胞类型变化：a为各形式富集的GO术语（GSEA），显示早期神经发生相关通路下调，后期突触相关通路下调；b为去卷积预测的细胞比例变化，显示特定细胞类型（如CycPro、ExNeu-2）在部分形式中改变；c-f为WGCNA模块与ASD形式的关联，M5和M19显著富集ASD风险基因且跨形式下调。

亮点 3｜M5调控网络预测下游ASD基因

通过RcisTarget预测转录因子结合基序，发现M5和M1含有最多的上游调控因子。M5调控的基因在多个ASD相关模块中富集，且这些靶基因中261个为SFARI高风险基因（如CNTNAP2、NLGN1、FOXG1、PAX6等），其中68%与M5调控因子负相关，表明M5可能抑制这些基因表达。

图 4　Fig. 4。图4揭示M5作为上游调控模块：a为RcisTarget工作流程；b显示M5调控的模块及其权重；c为模块间调控关系热图；d显示M5与下游模块的负相关；e为M5内转录因子的下调情况；f列出受M5调控的ASD风险基因；g显示调控因子与靶基因的表达时序；h-i为M5蛋白互作网络及其与ASD风险基因的富集。

亮点 4｜蛋白互作网络验证M5核心复合物

IP-MS实验证实M5调控因子（如SMARCB1、SMARCA4、EP300、TBP等）在神经祖细胞中形成显著蛋白互作网络（P=0.001），且该网络直接富集ASD和NDD风险基因（SFARI OR=15, FDR=2.0×10^-4），支持M5作为多蛋白复合物在早期神经发育中发挥调控作用。

图 5　Fig. 5。图5验证M5蛋白互作网络：a为抗体特异性验证的western blot；b为IP-MS火山图示例，显示显著结合的蛋白；c-f为M5蛋白互作网络的连接度分析，表明其与SFARI基因、NDD基因等显著连接，验证了数据库预测的互作。

亮点 5｜CRISPRi筛选验证M5调控因子的下游效应

在NGN2诱导的神经祖细胞中对26个M5候选调控因子进行CRISPRi敲低，发现11个因子的敲低导致其预测靶基因的显著差异表达（如RBBP5敲低影响超过1,250个基因），且差异基因富集ASD风险基因，证实了M5调控网络的功能相关性。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： hiPS细胞、皮层类器官、RNA-seq、WGCNA、RcisTarget、CRISPRi、IP-MS
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事神经发育、ASD遗传学、类器官建模或转录组网络分析的研究者。

带走一句

多突变并行类器官模型结合网络分析，能揭示ASD遗传异质性背后的收敛机制，为寻找共同治疗靶点提供思路。