期刊俱乐部：单细胞转录组学解码再生与疾病中的命运决定

· Research Updates

期刊俱乐部：单细胞转录组学解码再生与疾病中的命运决定

Journal Club: Single-Cell Transcriptomics Decodes Fate Decisions in Regeneration and Disease

Journal Club: 单细胞转录组学解码再生与疾病中的命运决定

引言

在生物学中，健康与疾病之间的差异往往取决于单个细胞的命运决定。为什么一个受损的肾细胞能够干净利落地再生，而其邻近细胞却会引发一系列导致纤维化的级联反应？植物根细胞如何精确地计时从分裂到伸长的转变？数十年来，这些问题都是在组织或群体层面进行研究的，掩盖了单细胞尺度上发生的关键异质性和动态转变。单细胞RNA测序（scRNA-seq）的出现，彻底改变了我们解析这些命运决定的能力，以前所未有的分辨率揭示了基因调控网络和细胞间通讯。

本月的期刊俱乐部汇集了三篇论文，它们共同展示了这种方法的强大力量。第一篇发表在《科学》杂志上的里程碑式研究，确定了一个动态的SOX9转录开关，它是决定受损肾上皮细胞是愈合还是驱动纤维化的主调控因子。第二篇同样发表在《科学》上的论文，利用拟南芥根中的时间序列scRNA-seq，绘制了激素油菜素内酯如何协调从细胞增殖到伸长的发育转变。第三篇发表在《细胞》杂志上的全面综述，提供了必不可少的方法学基础，综述了用于scRNA-seq注释、整合和细胞间通讯推断的计算工具，这些工具使得前两篇论文中的发现成为可能。

这些工作共同突出了一个共同主题：动态的、依赖于环境的转录程序充当了细胞传感器，整合内部状态和外部信号，从而决定功能结果。 通过超越静态快照，研究人员现在可以追踪单个细胞在再生与病理之间，或从一个发育阶段到下一个发育阶段的关键节点上的轨迹。

论文1：SOX9开关——从肾脏修复到纤维化

主要发现

  1. 单细胞分辨率下的二元命运决定： 急性肾损伤（AKI）后，激活转录因子SOX9的近端肾小管上皮细胞（PTECs）会遵循两条路径之一。成功恢复顶端-基底极性和上皮结构的细胞会沉默SOX9（SOX9开-关）并在无纤维化的情况下愈合。形成鲜明对比的是，未能恢复极性的细胞会维持SOX9表达（SOX9开-开），标志着一条病理谱系。
  2. SOX9开-开细胞被重编程并具有促纤维化特性： 持续的SOX9开-开状态以钙粘蛋白6（CDH6） 的共表达为标志。这些SOX9+CDH6+细胞经历表观遗传重编程，进入祖细胞样状态，并且至关重要的是，成为WNT配体（包括WNT4和WNT7b）的有效分泌者。
  3. 通过WNT信号枢纽直接联系纤维化： 单细胞分析显示，这些SOX9+CDH6+细胞产生"显著的单细胞WNT活性"，作用于邻近的损伤诱导的WNT响应细胞（WRCs）——主要是成纤维细胞。这种旁分泌信号驱动成纤维细胞活化为肌成纤维细胞，即纤维化中关键的胶原蛋白产生细胞。从SOX9开-开细胞中基因敲除WNT分泌装置减少了肌成纤维细胞的积累。
  4. 在人类疾病中的临床相关性： 对移植的人类肾脏的分析显示，包含SOX9+CDH6+WNT2B+单细胞的簇与纤维化程度之间存在强相关性，从而在小鼠模型中发现的基础上验证了人类病理学。

方法与研究设计

该研究采用了以scRNA-seq为中心的多模态谱系追踪方法。作者使用Sox9-CreERT2小鼠模型永久标记了损伤早期激活SOX9的所有细胞，使他们能够追踪这些谱系随时间的命运（从AKI到慢性肾病）。这整合了：

  • 单核ATAC测序（snATAC-seq）： 用于分析SOX9开-关与SOX9开-开细胞核的表观遗传景观，显示后者恢复了成熟的PTEC特征，而前者则保持在祖细胞样状态。
  • 谱系特异性基因组占据分析（CUT&RUN）： 以时间分辨的方式直接将SOX9结合与靶基因激活联系起来，确定Wnt2为直接的SOX9开-开靶点。
  • 功能验证： 使用条件性敲除模型（例如，成纤维细胞中的Ctnnb1缺失）和体外上皮培养系统来测试所识别通路的必要性。

意义

这项工作在理解器官疤痕形成方面，为上皮损伤和基质纤维化之间提供了一个机制桥梁，填补了长期存在的空白。它超越了SOX9在损伤后被重新激活的观察，表明其活性的持续时间是决定结果的关键因素。通过将SOX9+CDH6+细胞状态识别为病理信号枢纽，该研究提供了一个精确的治疗靶点。旨在促进SOX9失活或破坏这些特定细胞WNT信号的干预措施，可以在保留再生能力的同时预防纤维化——这是肾脏病学和再生医学的一个主要目标。

论文2：植物根发育中激素作用的细胞图谱

图：植物根发育中的油菜素内酯信号网络。来源：Science (开放获取)

主要发现

  1. 激素响应的时空特异性： 对超过210,856个拟南芥根细胞进行油菜素内酯（BR）或其抑制剂芸苔素唑处理的时间序列scRNA-seq，确定了伸长皮层是BR介导的基因表达的主要部位。这解决了BR在何处以及何时发挥作用以促进生长的模糊性。
  2. 发育转变的调控： 皮层细胞的拟时间轨迹分析表明，BRs促进从增殖状态到伸长状态的转变。这种转变与细胞壁组织和生物合成基因的特异性上调相关，为激素信号传导和细胞形态发生之间提供了直接的分子联系。
  3. 关键调控因子的识别： 从scRNA-seq数据推断基因调控网络，预测并随后通过实验验证了两个BR响应转录因子——拟南芥同源盒7（HAT7） 和GT-2样1（GTL1）——是皮层细胞伸长的关键调控因子。GTL1被证明与核心BR转录因子BES1相互作用以控制共同靶点。
  4. 细胞类型特异性功能： 作者使用组织特异性CRISPR技术，在皮层中特异性敲除BR受体BRI1，证明了该单细胞层中的BR信号传导对于伸长区正常的细胞扩张是必需的，而对分生组织细胞长度影响甚微。

方法与研究设计

核心创新在于将时间序列scRNA-seq应用于发育中植物器官的激素扰动。作者在处理BR或其拮抗剂后的多个时间点生成了根细胞的转录组，捕捉了响应的动态变化。然后他们：

  • 使用聚类和轨迹推断（例如，UMAP、拟时间）将细胞置于发育连续体上。
  • 为每种细胞类型和发育阶段推断上下文特异性基因调控网络，从相关性走向因果性。
  • 采用组织特异性基因操作（皮层特异性pCASP1::Cas9）进行精确的功能验证，从计算预测无缝过渡到生物学测试。

意义

这项研究展示了scRNA-seq如何能够揭示激素作用的上下文依赖性，这是植物和动物生物学中的一个主要挑战。它从一份BR响应基因列表，发展到一个空间和时间上解析的网络模型。通过识别HAT7和GTL1作为皮层中特异性的伸长调控因子，它为如何精确设计植物生长性状提供了蓝图。该数据集本身是植物学界的巨大资源，支持未来无数关于细胞类型特异性响应的分析。

论文3：单细胞生物学的计算引擎室

主要发现

这篇综述并未呈现新的生物学数据，而是综合了该领域的现状。其主要贡献是对方法学挑战和工具的结构化分析：

  1. 注释挑战： 该综述调查了细胞类型注释的方法，从基于标记基因的手动标记到自动分类器（SingleR、scPred）和参考图谱整合。它强调了错误注释对所有下游生物学解释的关键影响。
  2. 整合的必要性： 随着scRNA-seq数据集的激增，能够在保留生物学变异的同时消除技术批次效应的整合工具（Harmony、Seurat、SCVI）至关重要。该综述将此扩展到多模态整合，讨论了将scRNA-seq与空间转录组学（例如，用于Visium的Seurat、Tangram）、表观基因组学（scATAC-seq）和蛋白质组学（CITE-seq）相结合的方法。
  3. 推断细胞对话： scRNA-seq的一个主要前景是解码细胞间通讯。该综述整理了配体-受体相互作用工具（CellChat、NicheNet、CellPhoneDB），这些工具利用配体和受体对的表达来预测相互作用的细胞簇。它指出考虑新技术提供的空间背景对于验证这些推断的重要性。
  4. 持续存在的瓶颈： 作者指出了持续的挑战：处理稀疏和嘈杂的数据、某些平台的高成本和低通量，以及分析日益庞大的数据集（数百万细胞）的计算复杂性。

方法与研究设计

作为一篇全面的综述，该论文的方法是对已发表的计算工具和算法进行系统性的调查和比较。它将scRNA-seq分析流程结构化为逻辑步骤——预处理、聚类、注释、整合和通讯推断——并评估了每个类别中领先方法的优势、劣势和适用场景。

意义

这篇论文是任何进行或解释scRNA-seq实验的研究人员的必备手册。如果没有这里编目的复杂方法，论文1和2将不可能实现。通过对工具进行基准测试并澄清假设，该综述使研究人员能够选择稳健的分析策略，从而提高单细胞研究的可靠性和可重复性。它强调，生物学发现现在与计算创新密不可分。

综合与讨论

这三篇论文构成了一个强大的三部曲：一篇展示了通过先进单细胞方法实现的基础生物学发现，另一篇展示了类似方法在生命不同界中的应用，第三篇则提供了支撑前两者的方法论纲要。

解决核心问题的互补方法： 论文1和2都使用scRNA-seq来回答生物学中的一个核心问题：细胞状态转变是如何被控制的？ 在肾脏中，转变是从受损上皮到再生组织或纤维化组织。在根中，是从增殖到伸长的皮层细胞。两者都识别了主转录调控因子（SOX9；HAT7/GTL1），其动态活性决定了结果。两者也都强调了空间背景：导致SOX9活性锁定的失败的上皮极性，以及BRs促进伸长的特定皮层细胞层。
方法学创新的核心作用： 论文 3 明确指出了论文 1 和 2 中隐含的观点：它们的发现直接源于计算工具的开发。论文 1 对 SOX9+CDH6+ 枢纽的识别依赖于准确的聚类和注释。其与 WRCs 的关联则依赖于配体-受体推断工具。论文 2 中的轨迹分析和网络推断是综述中的核心主题。这些论文之间并无矛盾；相反，它们代表了从工具开发到生物学应用的无缝流程。

新兴的跨界模式： 一个有趣的模式是 “信号枢纽”细胞状态 的概念。在肾脏中，一个特定的、瞬态的细胞状态（SOX9+CDH6+）负责再生，但当其持续激活时，则成为一个病理性的信号中心（分泌 WNTs）。这反映了发育和癌症中的观点，即特定的祖细胞或过渡状态会分泌形态发生素来组织微环境。植物研究表明，一种激素如何定义一个具有独特转录输出的过渡状态（延伸皮层）。共同的主线是：细胞命运的决定不仅关乎其自身的命运，还关乎其指导邻近细胞的能力。

开放性问题与未来方向：

  • 治疗转化： SOX9 的开关能否通过药理学手段操控？小分子药物能否促进 SOX9 失活，或阻断来自 SOX9+CDH6+ 状态的特定 WNT 配体？
  • 保守原则： 其他修复过程（例如肝或肺纤维化）是否存在类似的“开关”转录调控因子？植物应激反应是否采用类似的持续性与瞬态转录状态来决定生存与死亡？
  • 多模态未来： 超越论文 1 中 scRNA-seq/snATAC-seq 组合的下一步，是在大规模水平上进行完全同步的多组学分析。论文 3 强调了这一方向。从同一个单细胞同时测量表观遗传学、转录组学和表面蛋白质，将如何改变我们对 SOX9+CDH6+ 等状态的理解？
  • 动态建模： 两篇生物学论文都捕捉了跨时间点的快照。未来的方向在于使用活体报告基因或通过计算整合时间序列数据到预测性动态模型中，实现对这类转变的真正实时追踪。

供实验室会议讨论

以下是 5 个引发讨论的问题：

  1. 方法与设计： 论文 1 使用了 Sox9-CreERT2 谱系追踪。这种基于“历史”的标记方法在定义 SOX9on-on 状态时存在哪些潜在局限性？新兴的转录组“时间戳”技术（例如，利用内含子读数或代谢标记）如何能补充这种方法？
  2. 普适性： SOX9 开关在肾脏近端小管中得到了优雅的证明。鉴于 SOX9 是软骨生成和其他过程中的关键转录因子，你认为其作为修复状态的“动态传感器”这一角色是肾脏特有的功能，还是一个可能适用于其在损伤后表达的其他组织的更广泛原则？
  3. 跨界联系： 论文 2 展示了激素（油菜素内酯）如何触发协调的基因网络以执行发育转变。在论文 1 中，触发或维持 SOX9on-on 状态的“激素”等价物是什么？是顶端-基底极性的丧失本身，还是由这种丧失产生的特定信号？
  4. 未来方向： 如果你要设计论文 1 的后续实验，你会优先考虑 (a) 高通量筛选选择性失活 CDH6+ 细胞中 SOX9 的药物，(b) 亚细胞分辨率空间转录组学以绘制 SOX9on-on 细胞周围精确的 WNT 配体梯度，还是 (c) 一项验证 SOX9+CDH6+ 细胞作为 CKD 进展预后生物标志物的人类队列研究？请证明你的选择。
  5. 更广泛的趋势： 论文 3 讨论了整合多模态数据的挑战。假设你拥有来自同一纤维化肾脏样本的 scRNA-seq、snATAC-seq 和空间蛋白质组学数据。利用这个整合数据集，你可以回答关于 SOX9+CDH6+ 状态的哪个具体生物学问题，而这是任何单一模态数据都无法单独回答的？

关键术语

  • 伪时间： 一种计算方法，利用 scRNA-seq 数据将单个细胞沿动态过程（如分化或对刺激的反应）的假设连续体排序，推断它们在该过程中的相对“进展”。
  • 谱系追踪： 一种遗传技术（通常使用 Cre-Lox 系统），用可遗传的标记（如荧光蛋白）永久标记一个细胞及其所有后代，使研究人员能够追踪初始细胞群随时间的命运。
  • 基因调控网络： 一组分子调控因子（如转录因子），它们彼此相互作用并与细胞中的其他元件相互作用，共同调控基因的表达水平，最终控制细胞状态和功能。论文 1 和 2 都从 scRNA-seq 数据中推断 GRN。
  • WNT 响应细胞： 论文 1 中对接收并响应 WNT 配体信号的细胞（主要是成纤维细胞）的功能分类，通常以靶基因（如 Axin2）的表达为标志。它们的激活是纤维化的关键步骤。
  • 多模态整合： 将同一生物样本中测量的不同类型数据（模态）——如基因表达（RNA）、染色质可及性（ATAC）和蛋白质丰度——进行组合的计算过程，以获得对细胞状态更全面的视图。

参考文献

  1. Shikhar Aggarwal, Zhanxiang Wang, David Rincon Fernandez Pacheco, Anna Rinaldi, A. Rajewski, J. Callemeyn, et al. (2024). SOX9 switch links regeneration to fibrosis at the single-cell level in mammalian kidneys. Science. https://doi.org/10.1126/science.add6371

  2. Trevor M. Nolan, Nemanja Vukašinović, Che-Wei Hsu, Jingyuan Zhang, I. Vanhoutte, R. Shahan, et al. (2023). Brassinosteroid gene regulatory networks at cellular resolution in the Arabidopsis root. Science. https://doi.org/10.1126/science.adf4721

  3. Changde Cheng, Wenan Chen, Hongjian Jin, & Xiang Chen (2023). A Review of Single-Cell RNA-Seq Annotation, Integration, and Cell–Cell Communication. Cells. https://doi.org/10.3390/cells12151970

Single-Cell Genomics scRNA-seq