引言
一个受精卵究竟是如何在空间和时间的三维坐标系中，精准无误地折叠出一个包含无数精细器官的复杂生命的？这是发育生物学中最令人着迷的核心问题。长期以来，我们观察生命发育的视角往往是割裂的：要么通过传统的显微镜捕捉细胞群的形态变化，却看不清内部基因表达的暗流涌动；要么将胚胎打碎，通过单细胞测序获取海量的基因数据，却彻底丢失了这些细胞在胚胎中的三维空间坐标。
当我们试图理解生命这张极其复杂的工程蓝图时，如果仅仅知道有哪些建筑材料，却不知道这些材料被放置在了大楼的哪一层、哪一个房间，显然是远远不够的。
3月12日，《Science》的研究报道“Whole-embryo spatial transcriptomics at subcellular resolution from gastrulation to organogenesis”，研究人员开发了一种全新的全胚胎成像平台。他们以前所未有的亚细胞分辨率，系统性地追踪了斑马鱼胚胎从原肠胚形成到早期器官发生阶段的基因表达动态。
突破维度的桎梏：看见基因表达的绝对坐标
为了在完整胚胎中实现高通量、高分辨率的基因表达分析，研究团队面临着巨大的技术挑战。传统的空间转录组学技术（Spatial transcriptomics）大多依赖于对组织进行二维切片，这在处理具有复杂三维结构的胚胎时，不仅容易丢失层级信息，而且往往只能捕获丰度较高的基因。此外，诸如原肠胚形成期这种伴随着剧烈细胞迁移和组织变形的关键阶段，任何二维平面的观测都显得力不从心。
为了打破这一僵局，研究人员巧妙地开发了一种名为 weMERFISH（全胚胎多重纠错荧光原位杂交）的新型成像平台。这一平台并没有简单地照搬已有的组织切片技术，而是针对数百微米厚的完整胚胎进行了三项关键的底层技术重构。
首先，他们引入了探针锚定技术。研究人员在合成靶向RNA的一级DNA探针时，加入了一种丙烯酰胺修饰（acrydite modification）。这一看似细微的化学修饰，使得探针能够与覆盖在胚胎样本上的聚丙烯酰胺凝胶发生共价交联。这就相当于在柔软脆弱的胚胎内部搭建了一层坚固的分子钢筋网。这种处理不仅使得样本在清除脂质和蛋白质后变得完全透明，还能在长达一个月的时间内保持极高的稳定性，彻底摆脱了RNA易降解的困扰。更重要的是，它允许使用极高浓度的甲酰胺反复洗脱荧光探针，从而支持多达60轮以上的成像循环。
其次，为了提高观测的灵活性，研究团队对编码本（Codebook）进行了创新设计。他们没有将用于多重成像的组合条形码直接整合到一级探针中，而是采用了一种基于接头探针（linker probes）的混合策略。这种设计使得研究人员可以根据需要，随时在单基因单分子荧光原位杂交（smFISH）和高通量MERFISH模式之间无缝切换。
最后，为了弥补共聚焦显微镜在对厚组织进行单分子成像时必然存在的信号衰减问题，研究人员引入了双层分支树信号放大系统（two-level branching signal amplification），成功将荧光信号强度提升了三倍。
通过这一系列深度的技术改造，研究团队成功在三个关键发育阶段（50%外包期的原肠胚起始阶段、75%外包期的原肠胚中期、6体节期的早期器官发生阶段）的斑马鱼胚胎中，对 495个关键基因 的转录本进行了亚细胞分辨率的精准定位。数据清晰地显示，在包含 103,646个细胞 的海量数据集中，weMERFISH的数据表现出极高的稳定性。不同生物学重复胚胎之间的基因转录本数量呈现出极高的相关性，相关系数达到了惊人的 R=0.996。与传统的smFISH相比，其转录本回收率达到了 50%至90%，而假阳性率被严格控制在每个细胞每个基因少于1个转录本。
重新定义细胞版图：隐藏在微米级边界中的秘密
当我们在亚细胞分辨率下审视整个胚胎的三维空间时，生命的复杂与有序远超想象。基于weMERFISH获取的单细胞内mRNA分子的绝对数量和种类，研究人员对胚胎内的细胞进行了重新聚类和空间映射。
以发育到6体节期的胚胎为例，研究团队在三维空间中将其精准地划分为了 51种 具有特定空间位置的细胞类型。值得注意的是，许多在以往单细胞RNA测序（scRNA-seq）研究中被粗略归为一类的细胞群体，在融入了空间坐标后，展现出了截然不同的异质性。例如，轴旁中胚层中的近轴细胞（Adaxial cells）曾被认为是一个单一的群体，但在weMERFISH的图谱中，它们根据所处的空间位置，被清晰地拆分为了三个亚群：位于已形成体节的亚群、位于正在形成体节的亚群，以及位于尾芽区域的亚群。这种空间上的细分，直接反映了细胞在不同发育阶段的成熟度差异。
空间转录组的威力不仅在于细化已知的细胞类群，更在于发现此前未知的、具有特定生物学功能的微型空间结构。在6体节期的胚胎中，研究团队在脊索后端的尖部发现了一个特殊的空间区域。这是一个狭窄的过渡带，介于前方正在分化成熟的脊索细胞和后方不断增殖的脊索祖细胞之间。数据表明，这个区域的细胞显著下调了与G1/S期转换相关的基因（如 mcm2、mcm3、mcm4 等），同时高度上调了细胞周期负调控因子 cdkn1ca。这提示我们，在剧烈的胚胎伸长过程中，胚胎巧妙地设置了一个 “细胞周期缓冲带”。这个缓冲带可能起到了协调前后两端不同发育节奏的作用：它允许前方的细胞有序地退出细胞周期并走向终末分化，同时保证后方的干细胞池能够持续提供细胞来源，以支撑尾部的不断延伸。
另一个令人意外的发现出现在原肠胚形成的起始阶段。传统观点通常认为，原肠胚的边缘区域（gastrula margin）是一个高度活跃的区域，大量的中内胚层基因在这里被激活。然而，weMERFISH的数据却揭示，在这个边缘区域，存在一组多达 83个基因 的表达水平显著低于胚胎的其他部分。基因本体（GO）富集分析表明，这些被大规模抑制的基因主要涉及有丝分裂DNA复制起始、染色质组织和细胞周期进程。这一发现挑战了既有认知，表明在细胞发生剧烈命运转变和形态学重排的原肠胚边缘，细胞周期和基础代谢活动可能会经历一次系统性的“踩刹车”过程，从而为随后的大规模分化腾出资源和时间。
此外，这些基因表达的空间边界具有极高的刻板性。在50%和75%外包期的重复胚胎之间，对应空间结构域的物理边界重合度达到了 20至30微米，这仅仅相当于一到两个细胞的直径。生命的底层秩序在微米尺度上展现出了惊人的精确度。
RNA的亚细胞漂流：空间位置决定命运
如果我们把视角进一步拉近，深入到单个细胞的内部，会发现RNA分子并非在细胞质中漫无目的地随机游走。它们的空间分布本身就蕴含着深刻的生物学功能。
weMERFISH的高分辨率使得研究人员能够量化每个基因的转录本到细胞膜的平均距离。数据表明，RNA的亚细胞定位呈现出强烈的基因特异性。例如，在50%外包期的深层细胞中，akap12b、kank1a 和 lima1a 等基因的转录本呈现出明显的近膜分布特征。特别是 akap12b，它在整个原肠胚形成期的所有细胞中，以及在6体节期的脊索和近轴细胞中，都稳定地锚定在细胞膜附近。相反，像 her1、cdx4 和 aspm 等基因的转录本则始终远离细胞膜，更倾向于在细胞核周围聚集。
更有意思的是，这种亚细胞层面的定位并非一成不变，而是跟随着胚胎的发育进程发生动态重组。研究团队系统性地比较了50%和75%外包期之间转录本距离细胞膜的变化。他们发现，znf185 基因的转录本在原肠胚形成的起始阶段，只是弥散地分布在包被层（Enveloping layer, EVL）细胞的细胞质中；但是，当胚胎发育到原肠胚中期时，这些转录本突然开始向EVL细胞的细胞间连接处高度富集。这种在特定时间节点发生的RNA定向迁移，暗示着 znf185 可能在原肠胚中期参与了维持或重塑上皮细胞屏障的复杂过程。
RNA分子在细胞内的微米级漂流，不仅反映了基因表达后的空间调控，更可能直接决定了所翻译蛋白质的局部浓度和作用靶点，这是理解细胞内部精细化管理的重要线索。
从数百到数万：以算法跨越观测的极限
尽管weMERFISH技术已经能够同时观测数百个基因，但这与构成生命蓝图的数万个基因相比，依然只占一小部分。为了突破这一观测极限，研究团队引入了单细胞多组学（scMultiome）数据，这包含在同一单细胞核层面同时测定的RNA测序（scRNA-seq）和染色质可及性测序（ATAC-seq）数据。
通过使用先进的Tangram计算方法，研究人员以weMERFISH测得的数百个基因为锚点，将scMultiome数据中的 25,872个基因 的空间表达模式成功推算（impute）到了完整的三维胚胎模型中。为了验证这一庞大推算结果的可靠性，研究团队将推算数据与斑马鱼信息网络（ZFIN）数据库中已有的105个高质量显色原位杂交结果进行了逐一比对。结果显示，无论是直接测量的基因还是通过算法推算的基因，其空间分布特征都与历史验证数据高度吻合。
这种维度上的拓展，使得这本生命蓝图瞬间变得丰满起来。通过这个被命名为MERFISHEYES的在线图谱系统，研究人员得以发现数百个以前从未被精确标注过空间位置的基因。例如，他们明确了在6体节期，cacna2d1a 和 tnni2b.1 特异性表达于近轴细胞，而 itm2cb 则富集于尾芽区域。
更令人振奋的是，既然可以通过RNA推算RNA，那么同样可以通过RNA作为桥梁，推算染色质的可及性。研究团队利用相同的方法，将 294,954个染色质开放区域（ATAC peaks）的空间位置映射到了6体节期的胚胎中。这意味着，我们不仅能看到哪些基因在哪些细胞里表达，还能看到为了促成这些表达，基因组上哪些特定的调控开关被物理性地“打开”了。在推算出的ATAC峰中，有73,036个显著富集于胚胎本部的特定细胞类型中（调整后的P值 < 0.05，log2倍数变化 > 0.5）。
从数百个观测数据点扩展到数十万个推算数据点，这绝非单纯的数量堆砌，而是为我们理解基因调控网络提供了前所未有的全基因组空间视角。
基因调控的折叠逻辑：定制与共享的智慧
有了全基因组的表达和染色质可及性数据，一个长期困扰发育生物学的核心问题终于迎来了大量数据的直接回应：那些在多个不同组织中同时表达的基因，究竟是如何被调控的？
在分析6体节期的胚胎时，研究团队鉴定出了2275个在特定组织中富集表达的基因。其中，超过一半（1316个基因）同时在两个或两个以上的不同组织或细胞类型中表达。这些基因是如何在不同时间、不同地点被精准唤醒的？
数据向我们展示了两种截然不同的底层调控逻辑。第一种逻辑可以被称为 “独立定制”。当一个基因在两个转录组特征差异较大的组织中表达时，它的顺式调控元件（cis-regulatory elements）往往是高度组织特异性的。以 hdlbpa 基因为例，这个基因同时在脊索（Notochord）和脊索前板（Prechordal plate）这两个组织中表达。但是，当你观察其基因组附近的ATAC峰时会发现，在这两个组织中打开的调控开关几乎完全不同。在其附近的8个主要ATAC峰中，有7个仅仅在脊索或脊索前板中的一个组织里处于开放状态。研究人员分析了所有48个在这两个组织中共同表达的基因，发现这种单一组织富集的调控模式是普遍存在的。这意味着，对于差异较大的组织，生命倾向于为同一个基因配备两套完全独立的“物理开关”，以确保互不干扰。
第二种逻辑则是 “资源共享”。当一个基因在两个功能相近或具有谱系联系的组织中表达时，调控逻辑发生了显著改变。以 palm1a 基因为例，它在近轴细胞和体节中共同表达，这两个组织在发育谱系上具有极强的连续性。数据表明，palm1a 基因两侧的许多染色质区域在近轴细胞和体节中是同时开放的。在这种情况下，相似的细胞状态允许它们共用一套广谱的调控元件。
这两种调控逻辑的揭示，充分说明了在漫长的进化过程中，基因组并不是随意堆砌的，而是形成了一套极其经济且安全的调控构架。差异大的组织用定制开关避免信号串扰，关系近的组织用共享开关提高响应效率。
突破预期的通讯网络：胚胎外层的隐秘对话
发育不仅仅是单个细胞按图索骥的过程，更是数以万计细胞之间复杂的社交网络。细胞与细胞之间通过分泌配体（Ligands）和表达受体（Receptors）进行着持续的通讯。由于weMERFISH保留了完整的三维组织结构，这为解析空间层面的细胞间通讯（Cell-cell communication）提供了完美的素材。
通过运用CellChat算法分析原肠胚期的空间转录组数据，研究团队不仅验证了经典的形态发生素（Morphogen）信号传导模式（如腹侧的BMP信号和边缘区域的Nodal/Wnt信号），还意外地在胚胎的最外层——包被层（EVL）中发现了未曾预料的隐秘对话。
长期以来，EVL被认为主要起着物理保护的作用，是一层相对处于生物学静息状态的胚外表皮结构。然而，数据清晰地显示，神经调节蛋白（Neuregulin） 和 颗粒蛋白（Granulin） 信号通路在EVL中表现出了极高的活性。具体而言，神经调节蛋白信号的配体和受体都严格限制在EVL内部表达，这意味着EVL细胞之间在进行极其密集的内部沟通；而颗粒蛋白的配体虽然在胚胎中广泛表达，但其受体却特异性地局限在EVL中，这暗示着深层细胞可能正在持续向外层发送某种未知的指令。
此外，研究人员甚至在EVL中检测到了经典的BMP、Wnt和FGF信号的配体和受体。这不仅颠覆了对EVL单纯作为“保护壳”的认知，更提示我们，在早期图式形成（Patterning）的过程中，这些经典的信号通路可能在胚胎最外层的细胞群中扮演着独立且复杂的调控角色。
捕捉时间的形状：RNA转录速率与胚胎的塑形
在生物学图像中，时间往往是被凝固的。我们得到的常常只是一张静止的切片或数据点。然而，weMERFISH的高分辨率带来了一个意外的惊喜：它不仅能够记录RNA在空间上的分布，还能通过区分细胞核内的未剪接RNA和细胞质内的成熟RNA，来推断基因表达的动态趋势，这就是所谓的 RNA转录速率（RNA Velocity）。
如果一个基因的转录本大量滞留在细胞核内，说明该基因正在被积极地转录激活（正速率）；反之，如果细胞核内空空如也，而细胞质中充满了该基因的转录本，则意味着该基因的转录已经停止，正在走向降解（负速率）。通过这种方式，研究人员可以在一张静止的荧光图片中，算出每一个细胞在转录层面正在朝着哪个方向“运动”。
当研究人员将这种转录层面的“时间箭头”投射到6体节期的胚胎空间中时，他们看到了一幅令人震撼的画面。在轴旁中胚层中，RNA转录速率的矢量场清晰地描绘出了一条从后方尾芽处的祖细胞，向前方的成熟体节平滑过渡的成熟轨迹。这表明，细胞的成熟并非是跳跃式的状态切换，而是在物理空间上呈现出一种连续的、渐进的转录状态演化。
更令人震惊的是，在原肠胚形成期，这些代表着基因转录变化方向的箭头，竟然与胚胎中细胞群的大规模形态发生运动（Morphogenetic movements）发生了不可思议的重合。在75%外包期，RNA转录速率的矢量图完美地勾勒出了边缘细胞向植物极的运动轨迹，以及腹侧-外侧细胞向背侧中轴发生会聚延伸（Convergent extension）的趋势。
这不仅仅是一个巧合。它暗示着在胚胎发育的最早期，全基因组水平的转录重编程方向，与细胞在宏观物理尺度上的运动方向是深度绑定的。基因表达的内在动力学，在很大程度上映射了胚胎外观形态的塑形过程。
命运轨迹的重构：动态重编程的基因启示录
虽然RNA转录速率能够预测短暂的未来，但由于原肠胚形成期的细胞发生着极其剧烈的物理位移，单纯的空间测序仍然无法告诉我们：现在位于这里的细胞，过去究竟经历了什么基因表达的变化？
为了彻底打通时间与空间的壁垒，研究团队开发了 MERFISH-FATE 分析流程。这一流程的核心思想是：将实时活体细胞追踪成像（Live-cell tracking）的数据，与固定时间点的weMERFISH空间转录组数据进行强制性的空间配准和融合。
通过这种强大的跨模态整合，研究团队以前所未有的清晰度，重构了细胞在原肠胚期剧烈运动过程中的基因表达演变史。他们发现，看似相同的最终基因表达模式，其背后的历史轨迹可能截然不同。
让我们聚焦中胚层特化因子 tbxta（也称T或Brachyury）。如果在50%外包期和75%外包期分别观察胚胎，你会发现 tbxta 在原肠胚边缘区域似乎一直都在表达。传统的推论自然是：边缘细胞在整个原肠胚期一直“维持”着 tbxta 的表达。
然而，MERFISH-FATE的数据狠狠地纠正了这一错觉。活细胞追踪结果表明，在50%外包期位于边缘的细胞，在发育到75%外包期时，大部分已经发生了内卷（Involution）并离开了边缘区域，在这个过程中，它们实际上“关闭”了 tbxta 的表达；而在75%外包期占据边缘位置的细胞，其实是来自于之前位于更偏动物极位置的细胞，它们在向边缘运动的过程中，全新“激活”了 tbxta 的表达。
因此，tbxta 在边缘的持续存在，并不是因为单批细胞的坚持，而是由于两批细胞的接力。它的调控模式实际上是一个复杂的复合体：既包含了像 mixl1 基因那样在早期边缘表达然后被抑制的过程，也包含了像 lef1 基因那样在后期边缘被全新激活的过程。这种维持、激活、抑制及混合调控的多样化策略，充分展示了细胞在面对复杂形态发生时的动态适应能力。我们看到的空间图式，不过是四维时空变幻下的一个切面。
边界的诞生：秩序源于原位重塑而非物理分拣
生命的巧妙之处在于它能够从混沌中建立起极其明确的物理边界，将不同的器官和组织严格分隔开来。在原肠胚发育过程中，脊索和体节之间是如何形成一道不可逾越的鸿沟的？这是发育生物学的一个经典议题。
研究人员首先量化了75%外包期胚胎中每个细胞周围50微米范围内的转录异质性。结果显示，在脊索祖细胞和体节中胚层祖细胞之间，存在着整个胚胎中最锐利的基因表达边界。在这个边界上，大量基因的表达水平在短短30微米的物理距离内发生了从极低到极高的剧烈跳变。
这种极其锐利的边界是如何诞生的？一种传统的观点认为（特别是在脊髓等组织中观察到的现象），由于原肠胚期细胞运动剧烈，不同命运的祖细胞最初可能是混合在一起的，随后它们通过差异性的细胞粘附力，像油和水一样发生了物理上的“分拣”（Cell sorting），最终形成了清晰的边界。
为了验证这一假说，研究团队利用光片显微镜（Light-sheet microscopy），从囊胚晚期开始，持续追踪了脊索和体节祖细胞的长距离运动轨迹。然而，追踪结果与物理分拣模型大相径庭。
数据清楚地表明，谱系限制在发育的极早期就已经完成。在囊胚期（Sphere stage）追踪的56个谱系中，只有5个包含了混合命运的细胞。当发育到75%外包期时，这些祖细胞完全且绝对地只产生脊索细胞，或者只产生体节细胞。更具决定性的是，从囊胚期到原肠胚中期的整个剧烈运动过程中，脊索祖细胞和体节祖细胞的运动轨迹几乎是完全平行的。统计显示，仅有不到4%的细胞跨越了这两种祖细胞之间的隐形边界。
这一有力的证据链表明，在50%外包期时，许多基因在脊索和体节的前体区域依然是广泛且模糊地表达的；而在接下来的几个小时里，这些细胞并没有发生大规模的位置交换。相反，它们是在保持相对空间位置不变的前提下，通过一种高度协调的基因 原位重塑（In situ remodeling），在原地集体改变了自身的转录状态。
它们不是通过“搬家”来划定界限，而是通过“原地换装”来构建秩序。细胞的转录程序与物理位置的协同进化，在此刻展现出了令人敬畏的精确性。
结语
从一个肉眼难以察觉的微小受精卵，发育为一个包含心跳、血液流淌和神经反射的完整生命体，是一场历经亿万年进化的奇迹。以往，我们只能盲人摸象般地去触碰这个奇迹的某个侧面。
然而，通过weMERFISH技术与单细胞多组学及活体成像的深度融合，我们终于拥有了一副可以看透生命折叠过程的三维透视镜。在这个多维度的坐标系中，我们不仅看到了数万个基因如何在微米级的空间内精准定位，看到了染色质的隐秘开关如何根据组织间的关系灵活配置，更看到了基因转录的动态如何驱动并重塑了细胞的命运与物理边界。
器官发育并非是在一维时间轴上机械播放的代码，而是一场在三维空间中实时响应、动态折叠的复杂交响乐。这项研究不仅为发育生物学提供了一份无价的资源图谱，更在底层逻辑上重塑了我们对生命空间秩序的认知。随着这些底层规律的不断破译，也许在不久的将来，我们不仅能够理解生命是如何被构建的，还能真正掌握干预甚至修复这本复杂蓝图的钥匙。
参考文献
Wan Y, El Kholtei J, Jenie I, Colomer-Rosell M, Liu J, Zhang Q, Acedo JN, Du LY, Codina-Tobias M, Wang M, Zheng W, Lin E, Chuang TH, Mayseless O, Sawh A, Mango SE, Yu G, Bintu B, Schier AF. Whole-embryo spatial transcriptomics at subcellular resolution from gastrulation to organogenesis. Science. 2026 Mar 12;391(6790):eadt3439. doi: 10.1126/science.adt3439. Epub 2026 Mar 12. PMID: 41818362.
声明：本文仅用于分享，不代表平台立场，如涉及版权等问题，请尽快联系我们，我们第一时间更正，谢谢！
往期热文：
