IF50.5｜空间转录组揭示人皮层分层与分区特化

MERFISH 单细胞空间图谱显示皮层六层结构在孕22周已分子定型，V1/V2边界孕20周已现。

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

研究背景

人类大脑皮层由六层结构和数十个功能分区组成，不同区域在神经元亚型和细胞构筑上差异显著。单细胞转录组研究虽已推进了对皮层发育的分子刻画，但细胞解离过程丢失了空间位置信息，无法判断相邻层和区域中形态相似的细胞类型是空间分离还是混杂。已有研究提出原图谱假说和原皮层假说来解释皮层分区形成机制，但人类胎儿组织获取困难，导致人脑皮层分区过程仍不清楚。空间转录组技术为在完整组织结构中同时分析分子和结构特征提供了可能，但基于捕获点的原位测序方法缺乏单细胞分辨率，而基于成像的MERFISH方法可实现真正单细胞水平的转录本检测。

科学问题

卡在哪里： 此前单细胞转录组研究丢失了空间背景，无法确定细胞类型在相邻层和区域中的空间分布关系；原位测序方法受限于捕获点包含多个细胞，难以达到单细胞分辨率；人类胎儿皮层发育过程中层结构建立和区域特化的时间线及分子机制仍不明确。

本文要回答： 利用MERFISH技术构建人类胎儿皮层高分辨率空间单细胞图谱，揭示皮层六层结构建立的时间窗口和区域特化的分子模式。

技术路线

作者主要用了：MERFISH、snRNA-seq、10x Visium、深度学习核分割、免疫荧光。

作者收集了孕15至34周的人类胎儿皮层样本，覆盖8个主要皮层区域。首先设计包含300个基因的MERFISH探针面板，涵盖细胞类型标志物和已发表scRNA-seq数据中的簇特异性基因。针对胎儿脑组织细胞密度极高的挑战，开发了基于CellPose 2.0的深度学习核分割模型，实现单核自动分割并适度膨胀以富集转录本。对约1600万个细胞进行层次聚类，得到8个细胞大类、33个细胞类型和114个亚型。同时利用snRNA-seq和Visium对连续切片进行互补分析，并通过ENVI方法将snRNA-seq数据与MERFISH整合，扩展基因覆盖范围。

图 1　Fig. 1。图1展示研究流程和细胞图谱概览。a为采样与工作流程示意图；b、c为MERFISH单细胞UMAP图，分别按H1细胞类和H2细胞类型着色；d为GW15-GW34主要皮层区域H2细胞类型的空间分布图；e为H2细胞类型标志基因表达点图；f为MERFISH H2簇与已发表中期妊娠scRNA-seq数据集的细胞类型对应热图。该图建立了整个图谱的框架和细胞分类体系。

核心亮点

亮点 1｜六层结构分子定型早于形态分层

通过分析兴奋性神经元亚型的层分布，发现皮层六层结构在孕22周已可由分子标志清晰界定，而细胞构筑分层直到孕34周才出现。层4-6在孕15周可见，层3-6在孕20周可见，符合由内向外的皮层形成顺序。层定义亚型的分布模式在不同区域和孕周间高度稳定，且随时间逐渐细化。这一发现将皮层分层的时间线提前了约3个月，表明分子分层与层构建几乎同步发生。

图 2　Fig. 2。图2展示皮层层的分子定义和层分布复杂性。a左为GW34 PFC中层2-6由选定EN H3亚型组织定义的空间图，右为层定义EN亚型的CD分布脊线图，虚线表示基于CD分布计算的层边界；b为GW34 PFC中EN-ET、深层EN-IT和上层EN-IT的脊线图，揭示进一步的空间复杂性；c为空间图和脊线图显示GW22时六层结构已可分辨。该图核心是证明分子分层早于形态分层。

亮点 2｜V1/V2边界在孕20周突然建立

在孕20周的枕叶皮层中，MERFISH分析揭示V1和V2之间存在由四对兴奋性神经元亚型定义的清晰边界，跨越所有皮层层次和板下层，其中层3和层4的转变最为急剧。V1特异性亚型仅存在于V1，而V2富集亚型在其他区域也有分布。通过扩展的960基因面板在孕18、20和21周样本中验证，该边界在孕20和21周出现，但孕18周尚未形成。这一离散边界与大多数皮层区域沿前后轴呈现的连续梯度模式形成鲜明对比，表明V1特化涉及独特的局部机制。

图 3　Fig. 3。图3展示层依赖基因表达的动态变化和区域富集基因。a为经典皮层标志基因在成人PFC和中期妊娠PFC中的层表达模式比较表；b为GW22 PFC和V2中CP内层依赖基因的层表达小提琴图，虚线框突出区域间差异；c、d为CBLN2等基因的z-score空间表达图；e为额叶富集基因的时序表达热图；f为多通道原位杂交标志物组合建议；g为前后轴梯度分布的EN亚型示例；h为跨五个神经元类别共享的前后富集基因。该图展示了层和区域特化的分子基础。

亮点 3｜NPY标记V1/V2边界的二元开关

在V1/V2边界附近，大多数差异表达基因在1-2毫米范围内逐渐变化，但NPY表现出明确的二元开关模式，在V1的层3和层4中强烈富集。该模式在孕20和21周出现，孕18周不存在。在猕猴胎儿皮层中也观察到由NPY标记的V1/V2边界，且V1特异性兴奋性神经元亚型表达NPY、PDZRN4和ABI3BP，表明这一分子边界在灵长类中保守。

图 4　Fig. 4。图4展示V1/V2边界的发现和特征。a为GW20-21时V1和V2无形态差异的示意图和染色图；b、c为枕叶皮层选定EN亚型的空间图，显示V1/V2之间由EN亚型急剧转变标记的清晰边界；d为各层和SP中边界转变的量化；e为成对亚型的层分布相似性；f为V1特异性亚型仅存在于V1，而V2富集亚型分布更广；g为边界附近基因表达的空间转变，NPY显示二元开关；h为Visium数据验证；i为边界在CP和SP中直线贯穿，但IZ、oSVZ、iSVZ和VZ均质；j为Visium斑点的UMAP图显示V1和V2平行轨迹在SP处分叉。该图确立了V1/V2离散边界的时空特征。

亮点 4｜V1神经元特化涉及两套转录程序

对V1富集基因集的时间分析显示，EN-IT亚型中V1富集基因在孕20周、孕34周和成年期之间发生显著变化，NPY和IGFBPL1在孕20周强富集但成年期消失，而TRPM3和CUX1在孕34周和成年期富集但孕早期无差异。相比之下，EN-ET亚型中V1富集基因集随时间更稳定，如TSHZ3和GABRA5持续中度富集。这表明V1特化至少涉及两套不同的转录程序：一套在孕中期瞬时激活于EN-IT，另一套在EN-ET中持续维持。

图 5　Fig. 5。图5展示V1神经元特化的时间进展和突触发生上调。a为GW34时V1和V2含不同神经元亚型的MERFISH空间图；b为V1富集亚型仅存在于V1的直方图；c为成年期V1和V2也含不同神经元亚型；d为各时间点V1-V2差异基因集在EN-IT中的累积表达热图，显示基因集随时间变化；e为NPY和IGFBPL1在GW20强富集但成年期消失；f、g为EN-ET中V1富集基因集更稳定；h为V1富集基因的GO分析富集突触组装；i为CellChat显示V1层4神经元接收的NRXN信号在EN来源中富集；j为网络图；k为免疫染色显示孕20周V1突触点密度高于V2。该图将V1特化与突触发生联系起来。

亮点 5｜V1层4神经元突触发生早期上调

通过ENVI将snRNA-seq数据与MERFISH整合，对V1特异性和V2特异性层3/4兴奋性神经元亚型进行差异表达分析，发现V1富集基因的GO分析显著富集细胞粘附和突触组装相关通路，且多数相关细胞组分为突触相关膜成分。CellChat配体-受体分析显示V1特异性层4神经元接收的NRXN-NLGN信号在EN来源中显著富集，而其他上层EN-IT仅从IN来源接收该信号。免疫染色证实孕20周V1突触点密度高于V2，但孕18周无差异，提示V1层4神经元在孕20周左右发生突触发生的早期上调。

生信可带走

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

• 方法栈： MERFISH、snRNA-seq、10x Visium、深度学习核分割、免疫荧光
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事神经发育、空间转录组或单细胞多组学分析的研究者，尤其是关注皮层区域化和层特化机制的人员。

带走一句

空间单细胞数据能揭示传统方法无法捕捉的发育事件，分子分层和区域边界可早于形态变化出现，分析时应整合空间与单细胞多模态数据。