IF36.1｜单细胞筛选揭示神经类器官对形态发生素的响应规律

系统评估形态发生素对神经类器官区域特化的影响，为优化类器官模型提供参考。

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

研究背景

人类中枢神经系统的复杂结构和功能源于早期发育中精确协调的模式化事件。神经板形成后，沿前后轴和背腹轴的形态发生素浓度梯度指导神经管区域化。目前对中枢神经系统区域化的认识主要来自非人动物模型，其细胞结构、行为和基因表达模式与人类存在差异。干细胞来源的神经类器官为在体外研究人类特异性神经发育提供了机会，但不同实验室的定向分化方案常产生不同结果，且方案开发效率低下，通常依赖少数标记基因的批量检测，并需数月培养后验证，忽视了早期神经上皮模式化这一关键阶段。因此，亟需对形态发生素在人类神经类器官早期区域化中的作用进行系统调查。

科学问题

卡在哪里： 已有多种生成区域化神经类器官的方案，但不同方案对同一脑区的产出存在差异，且缺乏对早期神经上皮模式化的系统研究。现有方法通常依赖少量标记基因的批量测量，通量低，验证周期长，难以整合各研究的发现并评估其对人类特异性生物学的相关性。

本文要回答： 系统评估形态发生素的时间、浓度和组合对人类神经类器官早期区域特化的影响，并考察细胞系和神经诱导方法对形态发生素响应的作用。

技术路线

作者主要用了：scRNA-seq、细胞哈希、SCENIC、微流控梯度、HCR原位杂交。

作者建立了一个实验流程，在96孔板中生成神经类器官，并在不同时间窗口、浓度和组合下施加形态发生素（如SHH、WNT、FGF8、RA、BMP4/BMP7等）。第21天收集类器官，解离为单细胞，使用细胞哈希技术标记不同条件，混合后进行单细胞转录组测序。通过参考相似性谱（RSS）与人类早期发育脑图谱比对，进行细胞类型注释。随后，利用SCENIC推断调控子活性，分析形态发生素与调控子的关联。此外，还比较了微流控梯度与静态浓度两种模式化方式对背腹轴模式的影响。

图 1　Fig. 1。图1展示了实验流程和整体数据概览。a部分显示类器官培养和形态发生素处理的时间线，b部分展示细胞哈希和单细胞测序流程，c-e部分为RSS-UMAP投影，显示不同细胞类型和区域身份。f-i部分显示不同细胞系在对照和处理条件下的细胞类型分布，j部分展示各条件的Shannon多样性指数，表明大多数条件产生多种细胞类型。

核心亮点

亮点 1｜形态发生素时间与浓度效应

研究发现形态发生素的作用具有时间和浓度依赖性。例如，RA在早期促进非神经外胚层，晚期则促进视网膜和 hindbrain 祖细胞；低浓度RA促进视网膜祖细胞，高浓度则富集 hindbrain 和后部底板祖细胞。WNT抑制剂XAV939在早期和中期处理可富集端脑祖细胞，而WNT激活剂CHIR则导致神经嵴和外周神经系统神经元产生。BMP4和BMP7在早期阻止神经特化，但BMP4在后期诱导非神经外胚层和视网膜命运，BMP7则维持多能性。SHH和FGF8均可诱导端脑命运，但SHH的作用窗口更早。这些结果通过细胞类型富集分析和HCR染色验证。

图 2　Fig. 2。图2聚焦于HES3细胞系中形态发生素时间和浓度对细胞类型的影响。a和b部分为热图，显示不同处理条件下细胞类型的富集/耗竭评分。c部分展示RA处理窗口和浓度对视网膜、hindbrain、非神经外胚层和后部底板祖细胞比例的影响。d部分为HCR染色，验证低浓度RA促进视网膜祖细胞（SFRP2+），高浓度RA促进腹侧后部命运（FOXA2+）。

亮点 2｜调控子活性揭示分子网络

通过SCENIC推断调控子活性，发现形态发生素与特定调控子存在关联。例如，SHH和RA共同关联的调控子包括GATA3和GLI1。SHH诱导的NKX2-1调控子激活需要较高浓度，而TCF7L2调控子仅在晚期和高浓度下激活。RA诱导FOXA2调控子活性在早期处理窗口，而HOX基因调控子呈浓度依赖性。CHIR相关调控子多与神经发生负相关，而SOX10等神经嵴调控子正相关。这些结果表明调节形态发生素的时间和浓度可以选择性激活不同的发育程序。

图 3　Fig. 3。图3展示组合筛选结果。a部分为HES3细胞系中各条件下细胞类型组成的饼图，显示不同形态发生素组合对细胞命运的影响。b部分为加入FGF8后的细胞类型组成变化。c和d部分为MMD距离分析，显示FGF8对特定条件的影响程度。e-g部分展示XAV939和SHH组合的非加性效应，通过非负最小二乘残差量化。

亮点 3｜组合筛选揭示形态发生素互作

设计组合实验，测试XAV939、RA、CHIR、SHH和FGF8等组合。在HES3细胞系中，XAV939与SHH组合产生协同效应，增加非神经外胚层细胞比例。在WTC细胞系中，RA单独处理后后部底板祖细胞占53%，加入SHH后比例升至87.8%，而加入FGF8则恢复hindbrain祖细胞。通过MMD距离和非负最小二乘残差分析，发现FGF8在HES3中对XAV939条件影响最大，而在WTC中对RA条件影响最大。这些结果揭示了形态发生素组合的非加性效应。

图 4　Fig. 4。图4展示可重复性筛选结果。a部分为实验时间线，b和c部分为RSS-UMAP投影，显示四种细胞系和两种诱导方法下的细胞类型分布。d部分为批次间可重复性分析，e部分比较不同诱导方法的细胞类型组成，f部分显示细胞系特异性细胞类型。g-l部分为AP-DV轴评分和变异性分析，显示高浓度形态发生素在各细胞系中效应一致。

亮点 4｜细胞系与诱导方法影响模式化

在四种hPS细胞系（H9、H1、WIBJ2、WTC）和两种神经诱导方法（MNIM和dual SMAD inhibition）中测试12种形态发生素条件。发现dual SMAD抑制产生更多神经元，但细胞系间变异性更大；MNIM则产生更多非神经外胚层组织，但可重复性更高。通过AP-DV轴评分，发现高浓度SHH和RA在各细胞系中效应一致，而低浓度条件则不一致。调控子分析显示，形态发生素对调控子活性的影响方向在细胞系间一致，但幅度受神经诱导方法影响。

图 5　Fig. 5。图5展示调控子活性分析。a部分为形态发生素-调控子网络，显示不同形态发生素关联的调控子。b部分展示调控子活性在细胞系间的一致性和变异性。c部分比较MNIM和dual SMAD条件下调控子对形态发生素的响应差异，如HOX调控子在CHIR处理下在两种条件下表现相反。

亮点 5｜微流控梯度与静态浓度比较

比较微流控梯度（MiSTR）和静态浓度（类器官）对SHH介导的腹侧化影响。发现类器官在低SHH浓度下即产生腹侧命运（如MGE祖细胞），而MiSTR组织在低浓度下仍保持背侧命运（如皮层/LGE祖细胞）。高浓度下，类器官产生下丘脑和底板命运，而MiSTR组织主要产生MGE命运。通过活体成像，发现类器官中NKX2-1表达呈局灶性并逐渐扩散，而二维培养中表达均匀。这些差异可能源于组织结构和细胞间信号放大。

生信可带走

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

• 方法栈： scRNA-seq、细胞哈希、SCENIC、微流控梯度、HCR原位杂交
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事神经类器官、发育生物学或单细胞组学分析的研究人员。

带走一句

形态发生素的时间、浓度和组合对神经类器官命运具有精细调控作用，但细胞系和诱导方法引入的变异性不可忽视。分析时应考虑这些因素，并利用调控子活性作为更稳定的分子读数。