Cell | 26个时间点、8大器官：STED-EC图谱如何重新定义血管内皮细胞发育研究？

构筑时空基座：七万个细胞的生命史诗

血管的发育是一个高度动态且连续的过程，任何粗糙的时间切片都可能导致关键信息的遗漏。为了捕捉这一过程中的每一个细微变化，研究人员展现了极其严谨的实验设计思路。他们并没有满足于几个标志性的发育节点，而是建立了一个名为 sted-ec (single-cell and time-series embryonic development endothelial cell atlas) 的图谱。

这个图谱的宏大与精细程度令人震惊。研究人员从小鼠胚胎发育的第7.0天 (e7.0) 开始，直至出生前的第19.0天 (e19.0)，以0.5天为间隔，连续采集了26个发育时间点的样本。这种高密度的采样策略，实际上是将一部原本可能只有关键帧的“幻灯片”，变成了一部流畅的“高清电影”。

在样本处理上，研究人员对脑、肺、肝、肾、心脏等8个主要器官进行了精细的解剖和标记。通过单细胞转录组测序技术 (scrna-seq)，他们最终获得了 70,915 个高质量的内皮细胞数据。

通过无监督聚类和内皮评分分析，研究人员在e7.0的胚胎中就检测到了表达 cdh5、pecam1 和 kdr 等典型内皮标志物的细胞群，这与之前的谱系示踪结果高度吻合，确认了内皮细胞发生的起始点。随后，利用基于随机森林分类器 (random forest classifier) 的迭代注释策略，研究人员成功地在这些细胞的“海洋”中，识别出了非器官特异性内皮细胞 (unorganotypic ecs)、中枢神经系统 (cns) 内皮细胞、肺内皮细胞、淋巴内皮细胞、心内膜内皮细胞、肝内皮细胞、胃肠道内皮细胞、肾内皮细胞以及肌肉和皮肤内皮细胞等多个亚群。

这不仅仅是一个数据的堆砌，更是一个从混沌到有序的梳理过程。数据显示，随着胚胎发育的进行，器官特异性内皮细胞的多样性显著增加。例如，肝脏内皮细胞在 e11.5之后出现了显著扩增，这与肝脏在发育中期的快速膨胀在时间上完美契合。sted-ec图谱的建立，实际上为我们提供了一个标准化的参考坐标系，让我们能够定量地去衡量每一个内皮细胞在发育长河中的位置。

关键的时间窗口：身份认同的觉醒

拥有了如此高分辨率的图谱，首要解决的问题便是：内皮细胞究竟何时获得其“器官户口”？

为了回答这个问题，研究人员设计了一套巧妙的评分系统——“成熟度评分” (maturity score)。他们利用成体器官特异性内皮细胞的组合分子标记，对胚胎各阶段的内皮细胞进行打分。这种评分并非简单的是非判断，而是量化了细胞与其最终成熟状态的相似程度。

分析结果揭示了一个令人深思的现象：在发育早期，各器官来源的内皮细胞成熟度评分普遍较低且保持稳定，这意味着它们在此时尚未分化出明显的个性。然而，在特定的时间点之后，这种评分出现了不可逆转的攀升。

数据显示，肺内皮细胞的特异性评分在 e12.0 之后开始持续上升，而中枢神经系统内皮细胞则更早，在 e9.5 左右便开始了其特异性分化之路。尽管不同器官的启动时间略有差异，但绝大多数器官特异性内皮细胞的身份确立都落在 e9.0 至 e13.5 这一时间窗口内。这表明，小鼠胚胎发育的中期 (mid-gestational stage) 是内皮细胞获得器官特异性的关键时期。

为了验证这一计算预测，研究人员在相应的发育时间点进行了免疫荧光染色。实验结果清晰地展示了关键标记物的时空表达模式。例如，肺特异性标记物 foxf1 在 e12.0 的肺部血管中开始清晰可见，而在更早的阶段则难以检测。

更有趣的是，通过应用基于最优传输理论 (optimal transport) 的 moscot 算法进行轨迹推断，研究人员重建了内皮细胞的分化路径。结果显示，除了心内膜内皮细胞展现出较为独特的路径外，其他器官的内皮细胞在 e7.0 至 e9.0 之间共享着相似的分化轨迹，而在 e9.0 之后则开始分道扬镳，走向各自的命运终点。

进一步的皮尔逊相关系数 (pearson's correlation coefficient) 分析量化了这种差异的演变。在 e12.5 时，心内膜、肝脏和淋巴内皮细胞的转录组差异已经显现；到了 e15.5，随着血脑屏障的建立，中枢神经系统内皮细胞与其他器官的差异变得愈发显著；而到了 e18.5，肺和肾脏内皮细胞的独特性也完全确立，整体的转录组差异水平已经逐渐接近成体状态。这一系列数据有力地证明，内皮细胞的器官特异性分化始于胚胎中期，并在随后的发育过程中不断强化和巩固。

基因模块与功能匹配：形式追随功能

在确立了“时间”之后，研究人员将目光投向了“内容”。不同器官的内皮细胞在获得特异性身份后，究竟表达了哪些基因？这些基因又如何服务于器官的生理功能？

通过对基因表达趋势的系统分析，研究人员发现内皮细胞的转录组模式可以被划分为四个主要模块 (modules)，这些模块根据基因在发育过程中的激活和抑制时序进行定义。令人惊讶的是，这四个模块在不同器官内皮细胞中的占比差异巨大。例如，在心内膜、中枢神经系统和肝脏内皮细胞中，模块1的基因占据了35%到40%；而在肺、肾和淋巴内皮细胞中，模块4的基因则占据了主导地位，比例高达50%到64%。这种差异暗示了不同器官内皮细胞采用了截然不同的基因调控程序。

当我们将这些差异表达基因 (degs) 与具体的生物学通路 (pathways) 关联时，一幅“形式追随功能”的精美画卷徐徐展开：

心内膜内皮细胞与心脏腔室发育、房间隔发育以及右心室形态发生相关的信号通路（如go:0003215, 0060413）表现出极高的活跃度，这与其作为心脏腔室衬里的解剖位置完美呼应。

中枢神经系统内皮细胞slc家族介导的跨膜转运（reactome: r-mmu-425407）以及葡萄糖跨膜转运蛋白活性（go:0055056）显著富集。这正是血脑屏障的核心功能所在——严格控制物质进出大脑，同时保证大脑的高能量代谢需求。

肝脏内皮细胞作为代谢枢纽，富集了氨基糖分解代谢过程（go:0046348）和类固醇激素生物合成过程。此外，wnt5a信号通路的富集则与其调控胎儿造血和肝细胞增殖的功能密切相关。

淋巴内皮细胞表现出与淋巴结发育、淋巴细胞趋化（go:0140131）以及抗原加工和提呈相关的高评分，凸显了其在免疫监视中的关键角色，而不仅仅是体液引流的管道。

具体的基因表达分析进一步佐证了这一点。例如，转录因子 gata5 在早期心内膜内皮细胞中特异性高表达，已知它能与nfatc相互作用调控内皮-心内膜的分化；而在肝脏内皮细胞中富集的转录因子 maf，则是肝血窦区域化分化所不可或缺的。

这些发现不仅仅是确认了已知的功能，更重要的是，sted-ec图谱展示了众多尚未被深入研究的器官特异性基因。例如，虽然 cdh5 在所有内皮细胞中普遍表达，但其他紧密连接和粘附连接蛋白（如 cldn11, cdh13）则显示出明显的器官偏好性。这种偏好性同样体现在分泌型因子的表达上，不同的白介素家族成员 (il1a, il7, il33) 在不同器官中的特异性分布，暗示了内皮细胞与其周围微环境之间存在着高度定制化的“方言”交流。

跨物种的鸿沟：相似的路径，不同的引擎

作为生物医学研究的基础模型，小鼠的研究结果在多大程度上能推演至人类，始终是一个绕不开的问题。利用sted-ec图谱覆盖全胚胎时期的优势，研究人员以肺发育为例，进行了一场跨越物种的时空对话。

研究人员整合了已发表的人类胚胎肺发育单细胞测序数据，时间跨度从受孕后第5周至第22周。通过典型相关分析 (cca) 等整合算法，研究人员发现，尽管物种不同，但人和小鼠的肺内皮细胞在整合后的umap空间中展现出了惊人的分布相似性。轨迹分析显示，两者的发育路径高度重合，都经历了从早期前体细胞向静脉、动脉及毛细血管亚群的分化过程。特别值得注意的是，人类的 tbx2+ 肺毛细血管内皮细胞与小鼠的 car4+ 亚群表现出了最高的转录组相似性。

然而，当研究人员深入到基因表达的微观层面时，却看到了意想不到的差异。在对比人和小鼠肺内皮细胞的基因表达趋势时，研究人员发现，在定义的四个基因表达模块中，仅有28.9%的基因表现出了相似的表达趋势，而高达 71.07% 的基因则显示出物种特异性的表达模式。

这种差异在具体的基因列表上更为直观：在肺内皮细胞富集的前200个基因中，只有 35.7% 是人和小鼠共有的（如 scn7a）。剩下的大部分基因则是物种特异性的（人类特异性占32.9%，小鼠特异性占31.4%）。功能富集分析显示，那些共有的基因主要集中在血管发育和上皮细胞增殖等基础生物学过程，而物种特异性的基因则可能涉及更精细的调控网络或适应性功能。

此外，细胞间通讯的分析也揭示了物种间的差异。在人类肺部，内皮细胞作为信号发送者的比例似乎更高，且存在着物种特异性的配体-受体相互作用对。

这一发现具有深远的启示意义：虽然小鼠和人类在器官发育的宏观结构和细胞类型演化上保持了高度的保守性（走着相似的路），但在分子驱动机制上却存在显著的差异（开着不同的车）。这提示我们在利用小鼠模型研究人类疾病或药物筛选时，必须对那些物种特异性的靶点保持高度警惕，同时也凸显了通过比较基因组学挖掘保守核心基因的重要性。

关键调控因子的发现与验证：casz1的角色

sted-ec图谱的价值不仅在于描述现象，更在于预测和发现新的关键调控因子。作为概念验证，研究人员通过图谱分析，锁定了一个名为 casz1 的转录因子。

分析显示，casz1 在肺内皮细胞中特异性高表达，并且在人和小鼠中保守。为了验证这一点，研究人员利用 casz1-em1zyliu 敲入报告小鼠（在内源 casz1 序列中插入了ha标签），通过免疫荧光染色证实了 casz1 蛋白特异性地存在于肺内皮细胞中，而在肝脏、心脏或脑的血管内皮中则未检测到表达（尽管在心肌细胞中有表达）。

此前关于 casz1 的研究主要集中在其对心肌细胞和神经元的功能调控上，其在内皮细胞中的作用几乎是一片空白。为了探究其功能，研究人员构建了内皮特异性 casz1 条件性敲除小鼠 (cdh5-creert2; casz1 fl/fl)，并通过他莫昔芬诱导在胚胎期进行基因敲除。

表型分析显示，内皮 casz1 的缺失导致了严重的肺血管发育缺陷。在 e15.5 时，敲除小鼠肺部的血管化区域显著减少。到了 e18.5，car4+ 的肺毛细血管内皮细胞——这一对气体交换至关重要的亚群——在敲除小鼠中急剧减少。单细胞测序分析进一步证实了这一发现，敲除小鼠的肺内皮细胞转录组特征发生了显著的“退化”。当将敲除小鼠 (e18.5) 的内皮细胞投射到 sted-ec 图谱上时，它们的转录组特征竟然类似于正常 e16.0 之前的状态。这意味着，失去 casz1 后，肺内皮细胞的发育被“卡”在了假腺体期 (pseudoglandular stage) 向小管期 (canalicular stage) 过渡的关键门槛上，无法完成成熟分化。

对差异表达基因的深入挖掘发现，casz1 的缺失导致了一系列关键血管生成调节因子（如 vegfa, kdr, tek, sox17）以及肺特异性分化标志物（如 foxf1, tmem100, tbx3）的表达下调。cut&tag测序技术进一步揭示，casz1 直接结合在这些关键基因的启动子或增强子区域，直接调控它们的转录。这不仅解释了表型，更构建了一个以 casz1 为核心的肺血管发育转录调控网络。

体外实验也佐证了这一点：在人脐静脉内皮细胞 (huvecs) 中过表达 casz1 会抑制细胞迁移并促进成管，而敲低 casz1 则产生相反的效果。这种功能与体内表型的互补，有力地确立了 casz1 作为肺内皮细胞发育关键“指挥官”的地位。

血管分泌这一信号：不仅仅是管道

血管内皮细胞早已不再被视为单纯的管道内衬，它们通过分泌“血管分泌因子” (angiocrine factors) 主动调控周围组织的生长和稳态。在 casz1 敲除小鼠中，研究人员观察到了一个有趣的现象：不仅血管出了问题，肺上皮细胞的发育也受到了严重的连带影响。

数据显示，在 e15.5 和 e18.5，敲除小鼠肺中 nkx2.1+ 的上皮细胞数量显著减少，肺泡ii型 (at2) 和i型 (at1) 上皮细胞的分化也受到抑制。edu掺入实验表明，这是由于上皮细胞的增殖能力下降所致。

既然 casz1 是在内皮细胞中特异性敲除的，为什么上皮细胞会遭殃？答案在于被破坏的“细胞间串扰” (crosstalk)。

通过对细胞通讯网络的分析，研究人员发现，在 casz1 缺失的情况下，从 car4+ 内皮细胞流向肺上皮细胞的信号强度显著减弱。在众多信号分子中，成纤维细胞生长因子1 (fgf1) 引起了研究人员的注意。sted-ec图谱显示，fgf1 是胚胎期肺内皮细胞特异性富集的血管分泌因子。而在 casz1 敲除小鼠的内皮细胞中，fgf1 的表达量出现了断崖式下跌。

cut&tag 数据证实，casz1 能够直接结合在 fgf1 基因的染色质区域（第18号染色体特定位点），直接驱动其表达。为了验证这一机制，研究人员建立了一个精巧的体外共培养体系：将 huvec 与 a549 肺上皮细胞进行共培养。结果显示，当敲低 huvec 中的 casz1 后，共培养的上皮细胞增殖显著减慢；而如果在培养体系中回补重组 fgf1 蛋白，则可以部分挽救这一增殖缺陷。

这一系列严密的证据链条，从转录因子结合、基因表达、信号分泌到靶细胞响应，清晰地描绘了 endothelial casz1 - fgf1 - epithelial proliferation 这一调控轴。它生动地展示了器官发育中“血管-组织”协同发育的精细机制：血管不仅要长好自己，还要通过分泌特定的生长因子，像“园丁”一样呵护周围实质细胞的生长。

从图谱到认知的飞跃

通过构建横跨26个时间点、覆盖8大器官的sted-ec图谱，这项研究不仅绘制了内皮细胞发育的精细路线图，更深刻地改变了我们对血管生物学的认知。

我们看到，内皮细胞的命运决定并非一蹴而就，而是在胚胎中期的特定窗口期，通过精密的转录程序启动的。我们看到，尽管不同器官的血管看似相似，但其内在的基因调控网络却有着天壤之别，以适应各自独特的功能需求。我们更看到，即便在进化上高度保守的小鼠和人类，其血管发育的分子细节也存在着巨大的差异，这为转化医学研究敲响了警钟。

更为重要的是，以 casz1 为例的发现过程，展示了“数据驱动假说”的强大力量。sted-ec图谱不仅仅是一个数据库，它更像是一个蕴藏着无数宝藏的矿山，等待着研究人员去挖掘更多像 casz1 这样的关键因子，去解析更多像 fgf1 这样的交流密码。

随着单细胞技术、空间转录组技术与表观遗传学技术的进一步融合，未来的血管图谱将更加立体和多维。而这项研究，无疑为我们理解生命之河如何通过精密的蓝图，流淌进每一个器官的深处，滋养生命的繁荣，奠定了坚实的基石。

对于发育生物学、再生医学以及血管相关疾病的研究者而言，这份图谱 ( http://sted-ec.ccla.ac.cn ) 将是一个不可或缺的导航仪，指引我们在探索生命奥秘的航程中，驶向更深、更远的未知海域。