题目：Single-cell transcriptomic profiling of zebrafish primary motor neurons reveals subtype-specific gene expression signatures

原文链接:https://pubmed.ncbi.nlm.nih.gov/41438093/

期刊：iScience

亮点

l 单细胞转录组学定义斑马鱼初级运动神经元亚型

l 七个pMN簇跨越早期、中期和末期分化状态

l 关键转录因子定义RoP、MiP和CaP神经元的身份

l 受体-配体图谱揭示驱动运动轴突导向的信号

摘要

运动神经元（MNs）将中枢神经系统与骨骼肌相连，实现运动功能。在斑马鱼中，RoP、MiP和CaP这三种运动神经元亚型已通过解剖学明确界定，但其分子特征仍未得到充分研究。我们对发育中的Tg(mnx1:gfp)斑马鱼胚胎进行单细胞转录组分析（该胚胎中的运动神经元被特异性标记），共鉴定出14个细胞簇，其中包括分化中的初级运动神经元。亚簇分析和拟时序分析进一步解析出7种初级运动神经元亚型，这些亚型被划分为早期、中间和终末分化状态，每种状态均具有独特的基因表达谱，涉及与轴突导向及细胞命运相关的转录因子和膜受体。单分子荧光原位杂交（smifish）的空间验证证实了脊髓内亚型特异性标志物的表达。蛋白-蛋白相互作用分析揭示了亚型特异性受体与运动轴突环境中存在的细胞外信号分子之间的潜在相互作用。本研究明确了初级运动神经元亚型的分子特征，为深入理解斑马鱼发育过程中运动神经元的身份特征及其在运动轴突决策中的作用提供了新见解。

图形摘要

引言

运动神经元（MNs）是连接中枢神经系统与肌肉并产生运动的关键结构。运动轴突与其肌肉靶标之间的精准连接对于实现协调运动至关重要。斑马鱼相较于其他脊椎动物结构相对简单，已成为研究运动神经系统发育的强大模型。 运动神经元的特化依赖于底板来源的音猬因子（Shh）信号通路，该通路诱导转录因子（TF）olig2的表达； 而初级运动神经元（pMNs）的数量则受Delta/Notch信号通路调控。在每个半体节中，根据在脊髓中的位置可区分出三种初级运动神经元亚型：嘴侧初级神经元（RoP）、中间初级神经元（MiP）和尾侧初级神经元（CaP）。受精后17小时（hpf），运动轴突沿共同路径向腹侧延伸，直至抵达对应水平肌隔的选择点。在此处，每条运动轴突会分支并支配特定区域：RoP支配外侧肌节，MiP延伸至背侧肌节，CaP则投射至腹侧肌节。在某些情况下，还存在第四种亚型，称为可变初级神经元（VaP），其发育特征与CaP相似，但会停留在选择点，并最终在受精后36小时退化。次级运动神经元（sMNs）在初级运动神经元生成之后产生，并于初级运动神经元形成8小时后开始支配肌节。

运动轴突在选择点的发散对于完整的肌神经支配至关重要。pMN的轴突自主导航，彼此独立，但它们的引导取决于与局部细胞外信号的相互作用。发育中的运动轴突由吸引和排斥的引导线索引导，这些线索可以是来自局部微环境的化学或机械信号。化学信号包括沿着运动轴突路径固定的众所周知的扩散分子或蛋白质，也就是接触介导的线索。在生长锥处，受体使MN能够检测和响应细胞外信号。在斑马鱼中，由于它们在未来生长的运动轴突轨迹中的位置，已知慢肌肉前体（也称为近轴细胞）是细胞外引导分子的主要生产者。

在脊椎动物和无脊椎动物中进行的几项研究表明，MN被逐步确定为不同的亚型，每个亚型对轴突引导线索的反应不同。这些亚型由特定的基因表达模式决定，特别是TF，包括LIM同源域蛋白的差异表达。在斑马鱼中，pMN身份的规范还涉及编码LIM同源域TF的基因的调节表达。每个pMN亚型暂时表达LIM基因的特定组合，包括islet1和islet2。 LIM基因表达的启动先于pMN分化，这可以通过几个小时后发生的轴突发生来证明。所有pMN在时间他们退出细胞周期，而轴突发生开始，CaP轴突主要表达islet2而MiP和RoP轴突主要表达islet1。有趣的是，Hutchinson和Eisen 证明，尽管islet1和islet2的表达模式受到严格调控，但这两个基因具有冗余功能。这表明额外的因素可能平行或上游作用，以交换意见的特异性到pMN亚型。Nkx6.1调节islet1表达的第二阶段，而旁轴中胚层派生信号也塑造pMN亚型身份。然而，我们对pMN身份和驱动选择点的电机轴突分歧的基础机制的理解在斑马鱼中仍然有限。

为了填补这一空白，我们使用单细胞RNA测序（scRNAseq）方法来表征(GFP+)分离细胞中MNs亚型的基因表达谱，这些细胞来自发育中的Tg（mnx1： gfp）斑马鱼胚胎，其中mnx1启动子驱动MNs中的绿色荧光蛋白（GFP）表达。数据的生物信息学分析揭示了在pMN分化过程中获得的亚型特异性基因表达特征，并富集在编码TFs和细胞外信号分子受体的基因中。我们的发现为这些特异性表达特征提供了深入而全面的理解，为驱动pMN分化和运动轴突引导的机制提供了新的线索。

结果

单细胞RNA测序识别转录不同的簇，这些簇概括了斑马鱼胚胎中MN分化阶段

在选择点，pMNs发散以沿着追踪良好的亚型特异性定型轨迹导航，并支配亚型特异性肌瘤区域。为了破译pMN亚型的基因表达特征，我们决定对分离的pMNs进行单细胞RNA测序（scRNAseq）。由于pMNs仅代表所有主干细胞（180 pMNs/胚胎）的一小部分，因此我们在进行scRNAseq之前，使用来自24 hpf Tg（mnx1：gfp）胚胎主干的异质细胞悬浮液，对(GFP+)MN群体进行了富集排序。对scRNAseq数据的无偏见聚类分析显示，总共包含4,040个细胞的14个集群，除了MNs之外，还包括已知表达内源mnx1的细胞类型，例如中间神经元、胰腺细胞和脊索细胞。集群1,2,3和5在寡2+ /sox19a+ /pcna+细胞中富集，表明它们对应于神经元祖细胞。预计24 hpf斑马鱼胚胎中细胞分化不同阶段的存在。每30分钟添加一次体细胞，以轮齿-尾部进展，直到在24 hpf完成体细胞发生，导致在同一胚胎中同时存在口端成熟体细胞和尾端不成熟体细胞。在斑马鱼中，pMNs和腹侧纵向中间神经元（VeLD）都出现在寡2+祖细胞中，但后者在有丝分裂后不表达islet1。使用elavl4作为有丝分裂后神经元的标记，我们确定了集群5和8，都富集在elavl4+/islet1−细胞中，这很可能对应于VeLD中间神经元。根据elavl4， islet1的基因表达水平，和mnx1，这是24 hpf pMNs的已知特征，集群3被分配给pMNs。因此，我们的数据涵盖了从寡2+祖细胞到成熟中间神经元和pMNs的分化阶段。

图1. 斑马鱼运动神经元的单细胞RNA测序分析（A）24小时受精后斑马鱼胚胎单个半体节中pMN发育的示意图。虚线代表水平肌隔。标注出了RoP（吻侧）、MiP（中间）和CaP（尾侧）初级运动神经元的位置。(B) 单细胞RNA测序分析的分步实验方案。红色箭头表示解剖方案。散点图展示绿色荧光蛋白荧光强度（x轴）与粒度（y轴）的关系。框选区域显示了分选时筛选出的GFP阳性pMN群体（标识为((GFP+))）。(C) 将测序读段比对至斑马鱼基因组并对荧光激活细胞分选分离的细胞（标识为(GFP+)）进行聚类后得到的均匀流形逼近与投影（UMAP）图。每个聚类以不同颜色突出显示。(D) 展示分离的GFP阳性细胞中olig2、sox19a、pcna、elavl4、islet1和mnx1基因表达情况的UMAP特征图。另参见表S1。

不同的基因表达特征描述pMN亚型和轨迹

为探究所鉴定的运动神经元前体细胞（pMNs）能否基于其亚身份进行聚类，我们在整个数据集中选取了表达编号为(levels >2)的mnx1阳性细胞（对数归一化表达值），以特异性分离pMNs，并开展了新的聚类分析。分析结果显示共得到7个亚簇，总计包含390个细胞。对不同亚簇间重叠的基因标志物进行分析，发现其存在统计学上显著的相似性（这是不同细胞亚型的典型特征），同时也存在显著差异，表明这些亚型具有独特的基因表达特征。对亚簇进行拟时序分析后，揭示了pMNs的成熟进程：亚簇1代表低分化状态，亚簇3、6和7代表高分化状态；亚簇2、4和5则包含处于中间分化状态的细胞。

为进一步阐明运动神经元分化过程背后的分子机制，我们对参与拟时序轨迹构建的前1000个基因进行了基因本体论（GO）富集分析。富集程度最高的GO条目包括“轴突发育”“神经元投射形态发生”和“神经元分化”，凸显了与pMNs成熟相关的关键生物学过程。截至目前，已被鉴定的pMN分化相关基因标志物数量较少。 研究人员通过原位杂交技术观察到，在pMNs分化与特化过程中LIM基因的表达呈现连续波次变化，这使得在分化细胞中鉴定特定亚型成为一项具有挑战性的工作。

因此，我们首先检测了islet1和islet2a在高分化细胞富集的亚簇（亚簇3、6和7）中的表达情况。结合islet2a在亚簇6中的高表达，以及单分子荧光原位杂交（smiFISH）实验结果——该实验显示第三种运动神经元中islet2a表达水平极高——我们确定亚簇6为尾侧肌节前体细胞（CaPs）。腹侧肌节前体细胞（RoP）和内侧肌节前体细胞（MiP）亚型在islet1和islet2a的表达谱上表现出较高的相似性；不过，基于亚簇3的islet1表达量高于亚簇7，我们将亚簇3归为RoP亚型，亚簇7则为MiP亚型。

随后，我们对7个亚簇进行了轨迹分析，发现亚簇4处存在一个分支点。其中一个分支指向亚簇7（已鉴定为MiP亚型），另一个分支则指向亚簇5，在拟时序分析中，亚簇5是亚簇3（对应RoP亚型）的直接前体。相比之下，CaP亚型（亚簇6）呈现出一条独特的分化轨迹（图2E）。这些研究结果表明，存在一个共同的前体细胞群（亚簇1），其后续分化为两种中间前体细胞：一种是CaP亚型的前体（亚簇2），另一种是RoP和MiP亚型的共同前体（亚簇4）。

通过本次分析，我们得以追踪关键标志物（包括islet1、islet2a以及此前被鉴定为分化MiPs标志物的mnx2b）在各pMN亚型拟时序中的表达模式，如今已能清晰区分RoP与MiP两个谱系。然而，我们的研究数据并未识别出一个 VaP 亚簇，这可能是因为它们在脊髓中的丰度较低，也可能是因为其基因表达谱与另一个亚型（极有可能是 CaPs）过于相似。

图2. 分化中的pMNs呈现独特的基因表达特征（A）对初级运动神经元（pMNs，mnx1 表达量＞2 的细胞）进行亚聚类后得到的 UMAP 图。不同聚类以不同颜色表示。(B) 展示各亚聚类拟时间变化的 UMAP 图。(C) 对参与构建拟时间轨迹的前 1000 个基因开展的基因本体（GO）分析。(D) UMAP 特征图，显示 islet1 与 islet2a 在各亚聚类中的表达（左图）；针对 islet1、islet2a 基因，在 24 hpf Tg (mnx1:gfp) 胚胎上开展的整体标本单分子低成本荧光原位杂交（smiFISH）结果（右图）。RoP（嘴侧初级）、MiP（中间初级）以及 CaP（尾侧初级）运动神经元。比例尺为 15 μm。(E) 展示亚聚类轨迹分析结果的 UMAP 图。(F) 柱状图：拟时间轴上，各初级运动神经元亚型中 islet1、islet2a 和 mnx2b 的平均表达量（对数标准化表达值）。每个聚类均标注平均表达量。另参见附图 S1 与附表 S1

独特的转录因子组定义了运动神经元前体细胞的亚型

我们预期，将我们的单细胞RNA测序数据集作为整体进行分析，能更好地反映特定运动神经元前体亚型的身份，因为转录因子会通过组合作用来抑制或增强靶基因的表达。因此，我们分析了每种运动神经元前体亚型相较于另外两种亚型的前50个差异表达基因，以此来刻画运动神经元前体亚型之间最显著的特征。这些基因中有许多编码的蛋白质属于转录因子家族、膜受体家族、细胞骨架蛋白家族和离子通道家族nels。我们鉴定出的大量差异表达转录因子基因与神经元细胞命运特化相关，包括Hox蛋白、 Mnx蛋白、 Islet蛋白和Sox蛋白。 此外，鉴定出的受体如神经纤毛蛋白和丛蛋白、Robo蛋白、 整合素、 以及神经营养因子受体均与轴突导向有关。离子通道还通过整合机械信号支持轴突生长和导向。例如，我们鉴定为CaP亚型标志物的TRPC3已被证实参与特化神经元网络的发育。此外，多种离子通道对于整合突触信号和产生动作电位至关重要。 离子通道的特异性

每个pMN亚型中的离子通道差异表明其电生理特性存在变化。值得注意的是，转录因子（TF）基因家族在所有pMN亚型中差异表达程度最高，在RoP、MiP和CaP亚型中，前50个基因里分别有10%、28%和14%呈现差异表达。部分转录因子在多个pMN亚型中共同表达。其中检测到了已被鉴定的转录因子基因标记物，例如CaP簇中的islet2a和MiP簇中的mnx2b。然而，差异表达转录因子的点图分析清晰显示，分化的RoP、MiP和CaP亚型表达组合性的亚型特异性转录因子基因库：RoP亚型富含lmo1、mafaa、maff、nfe2l1b、nfe2l3、ebf3b；MiP亚型以cdx4、ddit3、hoxd9a和Jun的表达为特征；而CaP亚型则特异性表达isl2a、mafa、pbx3b、tbx3b和pdrm14。我们基于lmo1（一种共转录因子）以及cdx4和mafa（转录因子）在pMN亚型中高度显著的差异表达，且在可能干扰pMN信号的邻近组织（如肌细胞）中无已知表达的特点，选择这三个基因进行单分子原位荧光杂交（smiFISH）验证。使用靶向所选基因的特异性探针进行单分子低成本荧光原位杂交（smiFISH，表S6）验证了单细胞RNA测序（scRNAseq）数据：lmo1的染色信号主要出现在RoP亚型，cdx4信号集中于MiP亚型，mafa信号则在CaP/VaP亚型中显著（定量结果见图3C和3D）。这一结果支持了我们的亚聚类分析和细胞轨迹分析结果（图2E），即脊髓内这两个相邻亚型即便基因表达谱并非完全一致，也具有高度相似性。对pMN阳性亚型（RoP、MiP、CaP）占比的定量分析进一步证实了上述观察（图3D）。此外，在与mafa探针杂交的pMN中，与CaP相邻的四个VaP均呈阳性。总体而言，这些研究结果表明，pMN亚型表达独特的亚型特异性转录因子组合，这对确立其细胞身份至关重要，同时也支持了VaP是CaP亚群发生进一步自发性退化和退变的一个亚群的假说。

图3. 亚型特异性基因的身份及差异表达（A）饼图展示了每种pMN亚型前50个差异表达基因（DEGs）中包含的转录因子（TF）、膜受体、离子通道和细胞骨架相关基因的数量。包括吻侧初级运动神经元（RoP）、中间初级运动神经元（MiP）和尾侧初级运动神经元（CaP）。(B) 点图展示了三种pMN亚型中差异最显著的转录因子。点的颜色代表基因的平均表达量（对数归一化表达值），点的大小对应每个簇中表达该基因的细胞百分比。(C) 24小时受精期（hpf）Tg(mnx1:gfp)胚胎的全 mounts 荧光原位杂交（FISH），靶向lmo1、cdx4或mafa基因（呈白色）。箭头指示某一pMN亚型内的染色信号富集。包括吻侧初级运动神经元（RoP）、中间初级运动神经元（MiP）、尾侧初级运动神经元（CaP）和可变型初级运动神经元（VaP）。比例尺为10微米。(D) 各亚型（RoP、MiP或CaP）中表达每种标记物（lmo1、cdx4和mafa）的pMN占比，相对于其他亚型的占比。lmo1：(n_pMN=) 6，胚胎数（n_embryos）=4；cdx4：(n_pMN = 9)，胚胎数（n_embryos）=6；mafa：(n_pMN =7)、(n_embryos = 5)，费希尔精确检验（Fisher’s exact test）；(^{* * *} p value <0.001)。(^{* * *} p value <0.0001)。另见表S2-S6。

细胞外信号与运动神经元受体的配对决定了选择点处亚型的分化

我们接下来的研究重点是分析膜受体的特异性表达，以及这些受体在向肌肉靶标移动的过程中，与所处微环境中可能遇到的细胞外信号的潜在相互作用。pMN亚型细胞会沿着共同路径行进，直至抵达分支点，随后分化并支配整个肌节。目前，关于这一决策过程的相关研究仍较为匮乏。因此，我们推测对RNA测序数据进行生物信息学分析，以此预测配体-受体相互作用，或许能为阐明这一问题提供线索。在我们的在单细胞RNA测序（scRNAseq）数据集中，我们鉴定出了几种已知参与运动轴突导向并能与细胞外信号分子相互作用的受体，包括 Ephrin、神经纤毛蛋白（neuropilin）、丛蛋白（plexin）、Roundabout（Robo）、骨形态发生蛋白（BMP）受体、钙粘蛋白（cadherin）、整合素（integrin）、成纤维细胞生长因子（FGF）以及胶质细胞源性神经营养因子（GDNF）受体。其中，部分受体在所有运动神经元前体（pMN）亚型中均有表达，而另一些则表现出亚型特异性表达。我们近期研究发现，分化中的慢肌前体细胞（SMPs）是发育中运动轴突的细胞外导向分子的主要来源。为进一步探究pMN表达的受体与SMP产生的细胞外信号分子之间潜在的相互作用（二者可共同调控运动轴突导向），我们利用此前从分化中的SMPs中获得的scRNAseq数据。我们借助STRING数据库对每种pMN亚型进行了蛋白-蛋白相互作用网络预测分析，将本研究中鉴定的pMN受体表达情况（青色标注）与SMPs表达的已知结构性导向信号（品红色标注）进行整合。内侧运动神经元前体细胞（MiPs）和尾侧运动神经元前体细胞（CaPs）均表达itgb1a和itga6b，二者可形成整合素α6β1受体，该受体能与基底膜蛋白相互作用：层粘连蛋白γ1（lamc1）、层粘连蛋白α5（lama5）、Ⅳ型胶原α1（col4a1）以及Ⅳ型胶原α2（col4a2）（图4C和4D）。MiPs特异性表达ephb2b和robo1，其已知结合伴侣分别为ⅩⅨ型胶原α1（col19a1）、Netrin 5（Ntn5）和Netrin 1a（Ntn1a）（红色圆圈，图4C），这表明它们在引导MiP轴突背向生长并越过选择点抵达肌肉靶标的过程中发挥核心作用。我们还分析了运动轴突路径中可扩散性导向信号与受体的预测配体-受体相互作用。例如，尾侧运动神经元前体细胞（RoPs）特异性表达bmpr2a，其蛋白产物可与可扩散性蛋白Inhbb和Fsta相互作用（图S2A）。CaPs表达ntrk2a，该基因编码Ntrk2a蛋白，已知该蛋白可与Fgf17、神经生长因子B（Ngfb）和血小板衍生生长因子AB（Pdgfab）相互作用（图S2B），而MiPs则未表现出亚型特异性的相互作用（图S2C）。这些预测分析表明，亚型特异性受体与导向信号共同塑造了pMN轴突的导向过程，并为未来的实验验证提供了理论框架。

图4. 运动神经元亚型特异性的pMN表达膜受体与细胞外信号之间的相互作用网络（A）点图展示了三种pMN亚型中差异表达最显著的膜受体和类受体蛋白。点的颜色代表基因平均表达量（对数标准化表达值），点的大小对应每个簇中表达该基因的细胞百分比。（B–D）利用STRING数据库构建的RoP（B）、MiP（C）和CaP（D）型pMN表达的受体（青色）（ArrayExpress数据库编号：E-MTAB-15160）与慢肌前体细胞表达的结构型细胞外蛋白（品红色）（SRA数据库编号：PRJNA973246）的蛋白质-蛋白质相互作用网络。网络通过富集分析(p-value <10^{-16})构建。亚型特异性相互作用以红色圆圈标注，共有相互作用以黑色圆圈标注。RoP（吻侧）、MiP（中间）和CaP（尾侧）初级运动神经元。另见图S2。

讨论

传统上，三个pMN亚型是根据其在脊髓中的位置进行分类的，但其分子身份仍知之甚少。本研究发现，RoP、MiP和CaP既存在共同的基因表达模式，也具有独特的基因表达模式，这确定了它们独特的分子特征，从而明确证明它们是不同的pMN亚型。pMN亚型的基因表达特征主要包括转录因子（TFs）和受体。为细胞外信号分子提供了新的见解，不仅阐明了驱动pMN亚型多样化和运动轴突导向的机制。

我们的单细胞RNA测序（scRNAseq）数据揭示了运动神经元前体细胞（pMN）亚型内独特的基因表达模式。近期针对小鼠的单细胞RNA测序研究已发现三种分子特征不同的轴索运动神经元（MN）亚型。 将这些研究数据与我们的进行对比后，我们发现，在小鼠和斑马鱼的运动神经元亚型中，MAF、ETS、FOX和EBF家族转录因子（TFs），以及斑马鱼中作为小鼠Pbx3同源基因的Pbx3b，均存在差异表达（本研究）。同样，包括Nrp1a/Nrp1、Nrp2b/Nrp2、Epha4b/Epha4以及Robo家族在内的轴突导向受体，在不同运动神经元亚型间也表现出显著的差异表达。LIM基因已在多个物种中被证实，是通过调控与运动神经元亚型特化关键转录因子的差异表达来塑造运动神经元身份的核心因子。在斑马鱼中，对LIM基因的研究可追溯至1995年，当时通过常规原位杂交技术完成了相关验证。该研究证实，斑马鱼LIM同源框基因在运动神经元前体细胞中呈组合表达，不同的运动神经元前体细胞表达这些基因的不同组合。我们的单细胞RNA测序数据不仅证实了这些发现，还精准揭示了包括isl1、isl2a以及mnx2b在内的LIM基因的表达模式。

先前的移植实验表明，pMN 细胞在完成最后一次分裂后、轴突发生前就已确定了特定的命运。对我们单细胞RNA测序数据的生物信息学分析与这些发现相符，证明存在共同的pMN 祖细胞，因此pMN 亚型的身份并非预先确定，而是在其分化的后期阶段获得的。细胞移植实验显示，LIM 同源框基因的表达会影响pMN 在脊髓内的位置。 这些发现早期就表明，pMN 细胞表面表达的受体使其能够响应分化信号，而后续研究证实这些信号源自体节中胚层。

对pMN亚型基因表达谱的分析揭示了其在转录因子、膜受体和离子通道上的富集特征。与其他物种的相关研究报道一致，30种转录因子的组合表达指导着运动神经元（MN）的分化与特化：RoP亚型富含lmo1、mafaa、maff、nfe2l1b、nfe2l3、ebf3b基因，MiP亚型以cdx4、ddit3、hoxd9a和Jun基因的表达为特征，而CaP亚型则特异性表达isl2a、mafa、pbx3b、tbx3b和pdrm14基因。单分子荧光原位杂交（Smi-FISH）实验证实，新鉴定出的候选基因可可靠地作为各亚型的标记基因。

膜受体是pMN细胞在选择点进行轴突寻路和决策的重要参与者。我们之前的研究表明，分化中的脊髓前体细胞（SMPs）是基质组基因的主要来源，这些基因为生长中的运动轴突提供及时且空间调控的局部信号，引导其向肌肉靶标移动。 本研究中，我们的蛋白质-蛋白质相互作用网络分析进一步揭示了pMN细胞表达的受体与分化中SMPs产生、且存在于运动轴突通路中的细胞外信号分子之间潜在的体内相互作用。这些分子包括XV型、XVIII型和XIX型胶原蛋白、层粘连蛋白及腱生蛋白C，还包括可扩散的导向信号，凸显了这些相互作用在pMN细胞轴突导向过程中的复杂性与动态性。最重要的是，我们的配体-受体分析揭示了pMN亚型的特异性，这表明尽管所有pMN都遵循共同的发育路径，但其轨迹在后续阶段出现分化。

虽然膜受体和转录因子（TFs）通常分别与轴突导向和细胞命运特化相关，但在pMN亚型中差异表达的离子通道基因的发现却出乎意料。此前研究表明，在斑马鱼中，运动神经元（MN）亚型在轴突生长之前以及轴突延伸直至48小时胚胎期（hpf）期间表现出电生理差异，这提示pMN亚型内离子通道的亚型特异性表达或表达水平存在差异，这与我们的研究结果一致。此外，Kelly等人证实，与次级运动神经元相比，一组在受精后4天（dpf）初级运动神经元中特异性富集的电压门控离子通道基因，能够促进驱动游泳反应的快速放电。然而，从功能角度来看，目前仍不清楚pMN亚型为何表现出不同的电生理特性。这一问题仍需进一步研究，而本研究对pMN亚型特异性离子通道的鉴定将为相关研究提供助力。

近期的单细胞RNA测序研究揭示了斑马鱼运动神经元（MNs）中广泛的分子多样性。D'Elia及其团队利用标记运动神经元、中间神经元和少突胶质细胞的olig2转基因品系，对2天胚胎（dpf）的轴运动神经元进行了图谱分析，鉴定出15种运动神经元亚型，其中多种亚型的验证性标记基因（如nrp1a和gfra1b）与本研究中发现的标记基因存在重叠。他们还鉴定出mecom和prdm16作为背侧快速运动神经元的标记基因，这两个基因被证实对快速运动神经元的发育至关重要，不过本研究的数据集未检测到这两个基因。早期研究通过分析olig2:EGFP阳性细胞的转录谱，确定了24小时胚胎（hpf）时期运动神经元前体细胞（pMN）的转录特征；与之相呼应，本研究利用sox19a1将1、2、4和6号细胞簇归为神经元前体细胞。此外，针对成年斑马鱼的研究表明，运动神经元亚型的分子多样性会延续至成年期。总体而言，这些研究为本研究结果的对比验证提供了丰富的参考框架，同时也凸显了运动神经元亚型特化过程中既保守又具有发育阶段特异性的特征。

总之，我们的单细胞RNA测序分析全面呈现了斑马鱼运动神经系统发育过程中pMN亚型的全貌，研究表明pMN亚型是通过多步过程获得其身份的。我们已鉴定出亚型特异性的基因表达特征，包括转录因子和细胞外信号分子受体，这些特征可能在塑造斑马鱼pMN命运和引导运动轴突轨迹方面发挥关键作用。本研究为进一步探究斑马鱼中运动神经元发育、轴突导向及命运决定的精确机制提供了宝贵资源，也为其他脊椎动物的相关研究提供了参考。

关于木芮生物ABOUT MURUIBIO

苏州木芮生物科技有限公司（以下简称“木芮生物”）成立于2018年1月22日，是一家致力于以斑马鱼为模式生物提供生物医学基础科研服务、斑马鱼实验操作试剂盒、斑马鱼临床疾病模型、其他生物学实验试剂以及进行斑马鱼应用技术转化的高科技公司。木芮生物以斑马鱼为模式生物，建立起了上百种的临床疾病模型，深入探索相关疾病发生的深层次疾病机制， 进而为临床治疗疾病提供可行的治疗方案或者药物筛选机制。到目前为止，木芮生物已经建立起包括遗传、行为、 细胞、 生化分子等斑马鱼相关成熟的技术，并且依托这些技术成功将斑马鱼应用到基础科研、毒理测试、癌症药物筛选、 中药药效成分筛选、保健品开发和功效评价。木芮生物经过七年的快速发展，公司现已有核心创业团队成员20余人，其中博士研究生1人，硕士研究生6人，本科生10余人。木芮生物成立后，陆续与江苏省产业研究院和苏州纳米应用技术研究所等科研单位建立起了科研合作，建立了“模式生物技术与应用研发中心”；公司于2020年被认定为“国家高新技术企业”，于 2021年底获得苏州工业园区“领军成长企业”称号。