IF50.5｜转录组相似的神经元在功能与形态上存在地形差异

斑马鱼视顶盖中，同一转录组类型的神经元因空间位置不同而呈现功能与形态多样性。

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

研究背景

神经元可根据基因表达组合分为不同细胞类型，其转录组被认为编码了形态、连接和功能的遗传指令。在秀丽隐杆线虫和视网膜神经节细胞等系统中，转录组类型与形态功能高度一致。然而，在哺乳动物视觉皮层等区域，同一转录组类型的神经元在功能调谐和形态上表现出显著异质性，提示外在因素可能影响神经元表型。斑马鱼视顶盖是一个层状结构，接收视网膜的拓扑投射，并沿前后、背腹和浅深轴存在功能特化，是研究转录组、位置与表型关系的理想模型。

科学问题

卡在哪里： 此前研究多假设转录组类型、形态类型和功能类型等价，但缺乏在同一系统中同时获取转录组、功能和形态数据并系统比较的证据。斑马鱼视顶盖虽有单细胞转录组研究，但覆盖不全，且未将细胞类型与空间位置、功能反应和形态特征关联起来。

本文要回答： 系统绘制斑马鱼视顶盖神经元转录组类型图谱，并检验转录组类型是否决定功能反应和形态，以及空间位置在其中扮演的角色。

技术路线

作者主要用了：scRNA-seq、HCR原位杂交、双光子钙成像、转基因稀疏标记、形态追踪、SVM分类。

作者首先对6-7天斑马鱼视顶盖进行单细胞RNA测序，聚类得到66个神经元转录组类型。随后利用多重HCR原位杂交将标记基因空间定位到标准脑图谱，揭示转录组类型的层状分布。接着通过双光子钙成像记录视觉刺激下的神经元活动，再对同一批细胞进行HCR标记以匹配转录组身份，从而关联功能与类型。此外，利用转基因品系稀疏标记特定转录组类型的神经元，追踪其形态和投射。最后整合空间坐标，分析位置对功能和形态的影响。

图 1　Fig. 1。图1展示单细胞转录组聚类结果。a图显示45,766个细胞分为25个主要簇；b图将神经元亚群进一步细分；c和d分别展示兴奋性和抑制性神经元的再聚类；e图展示兴奋性神经元各簇的差异表达基因点图。整体呈现视顶盖神经元的高度分子多样性。

核心亮点

亮点 1｜视顶盖神经元转录组多样性

对6-7天斑马鱼视顶盖进行单细胞RNA测序，获得45,766个细胞，聚类为25个主要细胞群，包括3个祖细胞群、1个放射状胶质细胞群、14个神经元群和6个非神经元群。进一步细分神经元得到33个兴奋性类型（e1-e33）和33个抑制性类型（i1-i33），共66个神经元转录组类型。通过差异表达分析，为20个类型找到互斥的单一标记基因，其余类型需稀疏组合标记。该图谱比此前研究更全面，覆盖了更多细胞类型。

图 2　Fig. 2。图2展示标记基因的空间表达模式。a图显示在标准脑图谱的横切面和冠状面上测量像素强度；b图展示多个标记基因在SPV浅深轴上的z-score强度分布，揭示兴奋性和抑制性标记的层状分离；c-f图通过细胞质心最近邻距离聚类将SPV分为浅、中、深三层。

亮点 2｜转录组类型沿浅深轴分层排列

通过多重HCR原位杂交将标记基因定位到标准脑图谱，发现兴奋性神经元主要位于SPV浅层，抑制性神经元位于深层，胆碱能神经元介于两者之间。个别类型打破此规律，如兴奋性标记bhlhe23表达于深层。基于细胞质心最近邻距离的层次聚类将SPV分为浅、中、深三个分子层，表明转录组类型在空间上呈层状组织，且与神经元出生时间相关。

图 3　Fig. 3。图3展示组合标记的空间分布。a-c图针对atf5b+细胞，显示其亚型（如atf5b+ etv1+和atf5b+ cebpa+）在浅深轴和前后轴上的差异分布；d-f图针对sp5l+细胞，显示sp5l+ chata+和sp5l+ neurod1+等亚型的空间分离。

亮点 3｜组合标记揭示前后轴亚型分布

对表达atf5b和sp5l的神经元进行迭代HCR，检测亚标记基因组合。发现atf5b+ etv1+细胞位于SPV深层，而atf5b+ cebpa+细胞富集于浅层后部；sp5l+ chata+细胞位于中间层，sp5l+ neurod1+细胞集中于后部。表明组合标记可进一步细分转录组类型在前后轴上的空间分布，且与神经发生和分化进程相关。

图 4　Fig. 4。图4展示功能与转录组的关联。a图示意实验流程：先进行双光子钙成像，再对同一批细胞进行HCR标记；b图显示功能数据与HCR数据配准后，1,304个ROI被分配到9个标记基因之一；c图展示各标记基因类型的平均反应谱；d图显示所有标记阳性神经元间的反应相关性矩阵，呈现功能多样性；e图比较类型内与类型间的平均相关性。

亮点 4｜同一转录组类型内视觉反应异质

对1,304个同时具有功能记录和HCR标记的神经元分析，发现同一转录组类型内神经元对视觉刺激的反应存在多样性，但平均而言，同一类型内神经元间的反应相关性高于不同类型间（9个标记基因中6个显著）。例如atf5b类型对运动点刺激反应强，而itpr1b类型对 looming 和 OFF ramp 反应强。表明转录组类型对功能有一定约束，但并非完全决定。

图 5　Fig. 5。图5展示功能类型与空间位置的关系。a图对7,094个高反应神经元进行层次聚类，得到17个功能类型；b图UMAP嵌入显示功能类型分布；c图展示各功能类型的反应向量；d图将HCR标记的神经元映射到UMAP上，显示特定转录组类型在功能空间中的局部富集；e图展示转基因品系的功能映射；f图显示itpr1b+神经元中不同功能类型的空间分离；g图量化t/f-clusters的空间重叠；h和i图显示SVM分类结果，表明位置比功能更能预测转录组类型。

亮点 5｜空间位置塑造功能表型

将功能反应聚类为17个功能类型（f-types），发现没有f-type特异性地富集于某一转录组类型。将同一转录组类型和同一功能类型的神经元（t/f-clusters）按空间位置分析，发现它们常沿前后轴分离，如itpr1b+神经元中f-type 4和9几乎不重叠。量化空间重叠显示t/f-clusters的空间分离显著高于随机对照。SVM分类表明，细胞体位置比功能特征更能预测转录组类型。

亮点 6｜形态多样性同样受位置影响

利用转基因品系稀疏标记atf5b、cort、itpr1b和sp5l神经元，追踪形态发现同一转录组类型包含多种投射和分层模式。例如cort神经元包括单层窄分层、单层宽分层和双层分层等类型。atf5b和sp5l神经元在视顶盖前部倾向于分化为中间神经元，后部倾向于投射神经元。表明形态表型也受空间位置调节。

生信可带走

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

• 方法栈： scRNA-seq、HCR原位杂交、双光子钙成像、转基因稀疏标记、形态追踪、SVM分类
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事单细胞转录组、空间组学或神经环路研究，关注细胞类型与功能映射的研究者。

带走一句

做细胞类型分析时，别把转录组簇当成功能或形态的同质单元；空间坐标和发育历史是必须纳入的协变量。