中后脑类器官单细胞多组学图谱与形态发生素筛选:从中脑到后脑身份的程序化诱导

Single-cell multi-omic atlas and morphogen screening informs midbrain and hindbrain organoid engineering

📣 摘要

Azbukina 等人(2026, Nature Neuroscience, JIF=25.0)针对中后脑类器官(midbrain/hindbrain organoid)身份特化的精确调控这一长期难题,基于 104,452 个单细胞、跨越 7/15/30/60/90/120 天 5 个时间点的 scMultiome(scRNA-seq + scATAC-seq 联合)数据,系统比较了两种主流的 WNT-CHIR/SHH/FGF8 模式化诱导方案,并通过 48 组形态发生素组合筛选(43 组纳入分析,177,718 个细胞 snRNA-seq)扩展了后脑细胞类型多样性。核心发现:两种中脑类器官协议均产生约 57% 后脑特征神经元(包括 glycinergic / Purkinje 细胞),Protocol 1 生成 89% 多巴胺能神经元,Protocol 2 生成 82% 谷氨酸能 + 91% 中脑 GABAergic 神经元;通过基于 Perturb-seq 的多重单细胞 TF 敲除筛选(EBF3 / OTX2 / MYT1L / ONECUT2 等 12 TF),识别出 PAX6/PAX7 为腹侧延髓关键命运调节因子。意义:为类器官脑区身份的精确程序化诱导提供了单细胞分辨率的调控图谱与可重复的形态发生素组合方案。

核心技术清单:

• scMultiome 单细胞多组学(10x Genomics Multiome,RNA + ATAC 同步)
• MERSCOPE 空间转录组(Protocol 1 Day 30 切片)
• 双重 SMAD 抑制 + WNT/SHH/FGF8 模式化诱导
• Pando 框架 + SCENIC 转录因子调控网络推断
• Perturb-seq 多重单细胞 TF 敲除筛选(dCas9 + 12 TF gRNA 文库)
• 组合形态发生素筛选(48 组条件,AP 轴 + DV 轴双轴)

🔬 研究背景

中后脑结构(中脑/后脑/小脑)是大脑协调运动、感觉整合、认知功能的核心区域,其发育异常涉及帕金森病、自闭症谱系障碍、共济失调等重大神经疾病。但现有脑类器官体系的后脑细胞类型多样性严重不足,主流的中脑类器官诱导方案(基于 Qian 等 / Jo 等)虽标榜"中脑身份",实际产生大量非预期后脑细胞。核心问题:

1. 能否在体外精确程序化诱导中脑 vs 后脑身份的细胞命运?
2. 现有中脑类器官方案的真实细胞组成如何(中脑 vs 后脑比例)?
3. 哪些形态发生素组合可扩展后脑细胞类型多样性?
4. 后脑区域化神经发生的关键转录因子调控层级是什么?

本研究采用"两套协议对照 + 跨时间点单细胞多组学 + 48 组形态发生素组合筛选 + 多重单细胞 TF 敲除"四步策略回答上述问题。

🧪 实验方法

类器官培养

使用 3 株人 iPSC 系(HOIK1、WIBJ2 来自 HipSci 资源;WTC 来自 Coriell)。mTeSR Plus 培养基 + Matrigel 包板,ROCK 抑制剂 Y-27632(5 μM)用于传代后第 1 天。两种诱导协议:

• Protocol 1(Qian 等):连续 WNT-CHIR + SHH + FGF8 模式化
• Protocol 2(Jo 等):序贯 WNT-CHIR → SHH + FGF8 模式化

单细胞多组学

10x Genomics Multiome(scRNA-seq + scATAC-seq 同步)。104,452 个高质量单细胞横跨 7/15/30/60/90/120 天 5 个时间点。3 株 iPSC 独立重复。snRNA-seq 单独用于形态发生素筛选(177,718 个细胞)。

空间转录组

MERSCOPE 平台对 Protocol 1 第 30 天类器官切片(180 个标记基因 panel)进行原位空间转录组分析。

形态发生素组合筛选

48 组形态发生素处理条件,覆盖 AP 轴(WNT-CHIR、RA、FGF8/2/17、R-spondin 2/3)+ DV 轴(BMP4、胰岛素、SHH/SHH 激动剂 purmorphamine)。第 8-14 天施加处理,第 5 周收样,单核 snRNA-seq。每条件 2 个生物学重复。无监督聚类得 25-26 种分子特征细胞类型。48 组中 5 组因类器官发育过小排除,43 组纳入分析。

多重单细胞 TF 敲除(Perturb-seq)

诱导型 Cas9 iPSC + 12 TF gRNA 混合文库(EBF3 / OTX2 / MYT1L / ONECUT2 等),慢病毒转导 + 双 SMAD 抑制神经诱导构建嵌合型后脑类器官。Day 7-14 多西环素诱导 TF 敲除,Day 30/70 收样,单细胞 RNA + gRNA 身份同步捕获。

数据分析

• 调控网络推断:Pando 框架(TF-靶基因共表达 + 基序)
• TF 活性定量:SCENIC chromVAR
• 迁移学习映射:scPoli + 加权 k-近邻(wkNN)映射到原代人胎脑细胞图谱(Mannens et al)
• 形态发生素筛选聚类:Shannon 多样性指数 + 正则化线性模型

数据规模总览

📊 核心结果与配图

发现 1:两种中脑类器官协议均产生显著后脑细胞身份

*图 1 caption: 单细胞图谱揭示人类后脑类器官发育的两套协议对比。a, Protocol 1 / Protocol 2 实验设计(连续 vs 序贯 CHIR/SHH/FGF8 模式化)。b, 明场形态(7-120 天)。c, 免疫荧光 DAPI/FOXA2/MAP2/OTX2 染色(200μm 标尺)。d-e, UMAP 注释细胞类型 + 脑区身份 + 分化时间。f, 堆叠柱状图:两类协议在不同天数的细胞类型动态。g, 细胞类型标记基因热图。h-i, 类器官细胞映射到原代人胎脑参考(Mannens et al):Protocol 1 → 腹侧中脑,Protocol 2 → 背侧中脑。j, 与原代脑细胞的转录组相似度(中位数 0.6-0.8)。k-o, MERSCOPE 空间转录组:第 30 天 Protocol 1 类器官中祖细胞 / 神经元 / 多巴胺能细胞的有序空间分布。

发现 2:染色质可及性图谱揭示后脑身份建立的时序层级

*图 2 caption: a-d, UMAP 按细胞类型 / 脑区身份 / 分化时间 / 协议上色。e-f, t-SNE 染色质开放峰聚类 + DAR 高亮。g, DAR 在不同峰簇的占比。h, 26 簇 × 13 细胞类型的层次聚类热图(蓝-4 到红+4)。i, FSTL5 / CAMK4 基因位点染色质开放轨迹(多巴胺能簇特异性)。j, scPoli + wkNN 迁移学习流程图。k-l, 原代人胎脑参考映射 UMAP(端脑/中脑/后脑/小脑/延髓/脑桥)。m, 与原代脑细胞的分子主题相似度小提琴图。

发现 3:多重单细胞 TF 敲除识别后脑命运决定调控网络

*图 3 caption: a, Pando 框架推断的后脑 TF 调控网络(节点大小=PageRank 中心性)。b, 调控网络对细胞类群变异的解释度。c, 散点图:NPC vs 神经元阶段 TF 脑区偏好(后脑 vs 中脑)。d, NPC 与神经元脑区特征三维相关性(表达差异/基序/调控区证据)。e, 神经元 UMAP(谷氨酸/早晚期多巴胺/GABA/甘氨酸)。f, TF 主成分分析(脑区分离)。g, chromVAR 基序活性(LHX1/OTX2/FOXA2)。h, 双 gRNA TF-KO 慢病毒载体图谱。i, Perturb-seq 实验流程(嵌合型类器官 + 多西环素诱导)。j, Day 30/70 UMAP 细胞组成。k, EBF3/OTX2/MYT1L/ONECUT2 内源表达。l, SCENIC 调控网络活性热图。m, 对照 vs TF-KO 调控网络活性差异箱线图。

发现 4:组合形态发生素筛选显著扩展后脑细胞类型多样性

*图 4 caption: a, 35 处理条件组合索引(基础调控阶段 + 形态发生素组合)。b, 26 簇 UMAP + 对照/处理组对比。c, 主热图:每个细胞类型在各处理条件下的富集/耗竭显著性(绿=富集/灰=耗竭)。d, Condition 7/45 在 DV-AP 坐标的细胞富集密度。e, 所有处理相对于对照的 DV-AP 偏移散点图。f, 三组并列 UMAP:映射到原代人胎脑参考(最大存在得分 + 解剖分区 + 神经递质类型)。g, 对照 vs Condition 7 / 已发表类器官方案 vs Condition 45 的目标细胞群分布。h, 形态发生素筛选相对 HNOCA 人脑参考图谱的新增细胞状态增益 UMAP。i, 5 大脑区条件富集的小提琴图(脑区特异性神经元亚型)。j, 总结热图:不同形态发生素方案特异性富集的细胞类型。

发现 5:腹侧延髓特异 TF 调控层级 PAX6/PAX7

*图 5 caption: a, NPC 分化三种神经元的对照 vs Treatment 1 实验设计 + 比例矩阵 + 相关性热图。b, 神经类群转录组相关性热图(8 类细胞)。c, 类群空间映射(背侧延髓/背侧中脑小脑/运动神经元/腹侧延髓/腹侧中脑)。d, 形态发生素-转录调控网络(节点=TF / 椭圆=形态发生素)。e, 形态发生素调控谱层次聚类(1-相关系数距离)。f, 腹侧延髓 NPC 特异因子散点图(PAX6/PAX7/SOX6/FOXP1/LMX1A)。g, PAX6/PAX7 表达 vs 调控子活性空间验证。h, Braun et al. 外部数据集匹配验证。i, 腹侧延髓亚类(NPC / RELN+ / POU6F2+ / PAX2+)空间分层。

💡 与产品关联

与现有产品线的关联(按代理关系与论文方法学匹配度):

9-16 实战铁律提醒:代理关系按飞轮 Agent_Vendors 表;本文方法学核心 = scMultiome + MERSCOPE + Perturb-seq,所有代理产品均为方法学互补 / 后续验证延伸,非论文直接使用。

🏷️ 结论

1. 后脑细胞身份由 NPC 阶段已建立的 TF 调控层级驱动:scATAC-seq + Perturb-seq 双证据链(发现 2 + 3)证明脑区身份特化不是终末分化期被动响应,而是 NPC 期就已确定的程序——这从根本上改变了类器官诱导的策略(从"分化期调整"转向"NPC 期锁定")
2. 组合形态发生素筛选是扩展后脑细胞类型多样性的有效路径:43 组条件中,RA + BMP4 是最强调控变量;Condition 7(中剂量 BMP4)与 Condition 45(R-spondin + 低剂量 BMP4)分别代表中脑 vs 后脑身份特化极性——这是体外工程化后脑类器官的可重复配方
3. PAX6/PAX7 是腹侧延髓身份的关键命运调节因子:与 RA 信号下游耦合,且拮抗 SOX6/FOXP1/LMX1A——为腹侧延髓命运的精确诱导提供了明确的 TF 干预靶点
4. "中脑类器官"标签需重新审视:两种主流中脑类器官协议均产生 57% 后脑特征神经元——体外"中脑"标签与体内身份存在系统性偏差,后续研究需要更严格的区域身份映射(MERSCOPE 空间转录组 + Mannens / Braun et al. 原代脑参考)

📖 文献信息

关键参考文献 (5 篇已 CrossRef 核验):

1. Mannens, C. C. et al. An integrated transcriptomic cell atlas of human neural organoids. Nature (2024) — DOI: 10.1038/s41586-024-08172-4
2. Braun, E. et al. Comprehensive cell atlas of the first-trimester developing human brain. Science (2023) — DOI: 10.1126/science.adf1226
3. Qian, X. et al. Brain-Region-Specific Organoids Using Mini-bioreactors for Modeling ZIKV Exposure. Cell (2016) — DOI: 10.1016/j.cell.2016.04.032
4. Jo, J. et al. Midbrain-like Organoids from Human Pluripotent Stem Cells Contain Functional Dopaminergic Neurons. Cell Stem Cell (2016) — DOI: 10.1016/j.stem.2016.07.005
5. Treutlein, B. et al. Dissecting direct reprogramming from fibroblast to neuron using single-cell RNA-seq. Nature (2016) — DOI: 10.1038/nature18323

📥 原文 PDF 下载

• 联系人:李东伟 ｜ 📞 18045130660 ｜ 🌐 https://www.biokingdom.top

转发语: 这篇 Nat Neurosci 2026 把中后脑类器官这件事讲透了:Protocol 1/2 都不是真正的"中脑"协议(都产生 57% 后脑神经元),真正的脑区身份锁定窗口在 NPC 期不是分化末期。43 组形态发生素筛选 + 12 TF Perturb-seq 给出了可重复的配方和调控靶点(PAX6/PAX7 = 腹侧延髓命运开关)。🧠 神经类器官客户可以拿去当"为什么需要 Nanolive Eve + 3Brain HD-MEA 做功能验证"的最强销售证据 ⚡。