关于《transcriptomic landscape of the blastema niche in regenerating adult axolotl limbs at single-cell resolution》的学术研究报告

第一， 研究作者、机构与发表信息 本研究由Nicholas D. Leigh, Garrett S. Dunlap, Kimberly Johnson, Rachelle Mariano, Rachel Oshiro, Alan Y. Wong, Donald M. Bryant, Bess M. Miller, Alex Ratner, Andy Chen, William W. Ye, Brian J. Haas & Jessica L. Whited 共同完成。研究团队主要来自哈佛医学院骨科外科系、哈佛干细胞研究所、布莱根妇女医院、麻省理工学院和哈佛大学布罗德研究所以及哈佛大学干细胞与再生生物学系。该研究成果以题为《Transcriptomic landscape of the blastema niche in regenerating adult axolotl limbs at single-cell resolution》的论文形式，于2018年发表于学术期刊 Nature Communications。

第二， 学术背景与研究目的 本研究的科学领域聚焦于再生生物学，具体研究墨西哥钝口螈（Axolotl）这种具有终生再生完整四肢等复杂多组织结构能力的特殊两栖动物模型。尽管长期以来，科学家们通过组织学、移植实验和批量组织RNA测序（bulk RNA-seq）对再生过程中的伤口表皮（wound epidermis, WE）和芽基（blastema）进行了广泛研究，但构成再生肢体的各种细胞类型的分子身份仍不明确。批量测序技术只能提供细胞群体的平均基因表达谱，掩盖了关键细胞类型特异性转录本的信息，这阻碍了精确界定不同组织（如伤口表皮和芽基）在再生中的具体贡献。

因此，本研究旨在利用新兴的单细胞RNA测序（single-cell RNA-seq）技术，以前所未有的分辨率解析成年钝口螈稳态肢体及再生不同阶段肢体细胞的转录组景观。具体目标包括：1) 系统性地鉴定在肢体稳态和再生过程中存在的细胞类型及其多样性；2) 识别再生诱导的特异性基因；3) 推断芽基细胞可能的分化轨迹；4) 提出成纤维细胞样芽基祖细胞（fibroblast-like blastema progenitor cells）的分子身份。最终，本研究希望建立一个分子框架，为未来应用分子技术深入研究这些细胞群体在肢体再生中的作用奠定基础。

第三， 详细研究流程与方法 本研究是一个系统性的单细胞转录组学分析项目，主要流程包括样本准备、单细胞测序、生物信息学分析和实验验证四个核心部分。

1. 样本准备与单细胞捕获： * 研究对象与样本量： 研究使用了超过1年龄、体长17-21.5厘米的成年钝口螈，从未经历过截肢。共从13只不同的个体中收集了13个独立的样本，涵盖了稳态（homeostasis）前肢以及截肢后三个关键再生阶段：伤口愈合期（wound healing）、早期芽基期（early-bud blastema）和中期芽基期（medium-bud blastema）。每个时间点有2到6个生物学重复。 * 样本处理： 截肢后，收集整个再生体，包括伤口表皮和芽基。为了获得单细胞悬液，研究团队采用了一种优化的细胞分离方案，该方案基于并改进了先前发表的方法。特别指出，为了平衡细胞获取的偏差（该方案倾向于获取更多表皮细胞），在早期芽基期去除了全部伤口表皮，在中期芽基期仅将一半伤口表皮放入收集管，以确保获得足够的芽基细胞。 * 单细胞测序技术： 研究采用了基于微流控的高通量单细胞RNA测序平台——inDrops。该技术通过将单个细胞与带有独特条形码（barcode）的微凝胶珠共同包裹在油滴中，实现对数千个单细胞进行同时的转录本捕获和条形码标记。随后，构建文库并在Illumina NextSeq500平台上进行测序。

2. 数据处理与生物信息学分析： * 数据预处理： 使用修改版的inDrops处理流程对测序数据进行解复用，提取细胞条形码和唯一分子标识符（UMI）。将测序读数比对到已发表的钝口螈转录组参考序列上，并生成每个细胞的基因表达计数矩阵。 * 无偏聚类与细胞类型鉴定： 使用Seurat单细胞分析工具包对超过25,000个细胞进行分析。流程包括：数据质控（过滤低基因数或高线粒体RNA表达的细胞）、数据标准化、识别高变基因、进行线性降维（主成分分析，PCA）。随后，基于基因表达模式的相似性进行无监督聚类（即不预先使用已知的细胞标记基因），并通过t-分布随机邻域嵌入（t-SNE）方法在二维空间可视化聚类结果。根据每个细胞簇的特异性高表达基因（标记基因），研究者推断其可能的细胞身份，如表皮细胞、造血细胞、成纤维细胞样细胞、肌源性细胞、内皮细胞、周细胞、施万细胞等。 * 拟时序（Pseudotime）与轨迹分析： 为了研究细胞分化的动态过程，研究使用了两种算法： * Monocle： 用于分析表皮细胞的分化轨迹。通过将细胞沿假定的分化路径排序，揭示了从基底表皮/伤口表皮（起始点）到中间层，最终到外层小分泌细胞（SSCs）的向外分化过程。 * URD： 用于重建芽基内间充质细胞（如成纤维细胞样芽基细胞、施万细胞、肌源性芽基细胞、内皮细胞、周细胞）的发育轨迹。该算法通过计算细胞间的转移概率和伪时间，构建分支树状图，以推断不同细胞类型（如软骨样、关节样、成骨样细胞等）可能的分化路径和共同祖细胞。

3. 实验验证： * RNA原位杂交（RNA in situ hybridization）： 为了验证单细胞测序数据中识别的细胞类型特异性或再生特异性基因的表达模式，研究对多个候选基因（如KRT12, OTOG, FREM2, TRAC, APOEB, DPT, LUM, KRT17等）在稳态和不同再生阶段的肢体切片上进行了RNA原位杂交实验。这些实验直观地显示了目标基因在特定细胞群体或组织结构中的空间定位，与单细胞聚类分析结果相互印证。 * 组织化学染色： 使用α-萘酚乙酸酯（NSE）染色来检测芽基中的髓系细胞，作为对单细胞数据中发现的免疫细胞群体的补充验证。 * 苏木精-伊红（H&E）染色： 用于展示中期芽基期再生肢体的基本组织结构。

第四， 主要研究结果 本研究通过单细胞转录组分析，系统性地绘制了钝口螈肢体再生过程中芽基生态位的细胞图谱，取得了多项重要发现：

1. 揭示了再生肢体的高度细胞异质性： 在稳态和再生的不同阶段，共鉴定出11至19个转录组学上不同的细胞群体。这些群体涵盖了表皮谱系、间充质谱系和造血谱系。重要的是，几乎所有在再生过程中发现的细胞群体在稳态肢体中也存在，包括一群与推测的成纤维细胞样芽基（FLB）细胞相似的成纤维细胞群体。

2. 解析了伤口表皮的复杂构成与分化轨迹： * 在稳态表皮中，鉴定出了基底表皮、增殖表皮细胞、中间表皮、小分泌细胞（SSCs），以及推测的离子细胞和朗格汉斯细胞。 * 在再生伤口表皮中，离子细胞和朗格汉斯细胞缺失，但出现了Leydig细胞（表达CHS1和AGR2A）。研究还发现了多个中间伤口表皮群体。 * 拟时序分析证实，无论是稳态还是再生状态，表皮细胞都遵循一个从基底（祖细胞库）向外（SSCs）的分化轨迹。基因KRT17的表达沿此轨迹递减，其动态表达模式通过RNA原位杂交在真实的再生时间轴上得到验证。 * 识别了在基底伤口表皮中特异性表达的再生诱导基因，如FREM2，其表达模式通过原位杂交得到确认。

3. 鉴定了芽基生态位中丰富的免疫细胞： * 除了已知对再生至关重要的巨噬细胞，研究还发现适应性免疫细胞（T细胞和B细胞） 存在于所有采样时间点，包括中期芽基期。RNA原位杂交证实T细胞（表达TRAC）定位于芽基内部。 * 在先天免疫细胞中，鉴定了树突状细胞、新募集的巨噬细胞（表达TREM2）和中性粒细胞（表达CAMP）等群体。髓系细胞标记物APOEB的表达通过RNA原位杂交和NSE染色在不同再生阶段得到验证，显示其在伤口愈合期达到峰值。

4. 阐明了神经相关细胞和血管细胞的来源： 在数据集中发现了施万细胞、周细胞和内皮细胞群体。特别值得注意的是，在伤口愈合期，施万细胞主要不表达SOX2，而在中期芽基期，许多施万细胞获得了SOX2表达，这可能意味着它们去分化到了一个前体样状态，这与小鼠指尖再生中的观察类似。

5. 深入解析了芽基细胞的分子身份与分化潜能： * 确认了Pax7+肌肉卫星细胞来源的肌源性芽基细胞的存在，它们在早期芽基期高表达PAX7，在中期芽基期高表达MYF5，而SIX1和EYA1在两个时期都标记该群体。 * 在所有时间点都发现了成纤维细胞样芽基（FLB）细胞，它们高表达LUM、DPT、POSTN等成纤维细胞标记。RNA原位杂交证实了这些标记在芽基中的表达。 * 轨迹分析（URD） 揭示了FLB细胞可能的多向分化潜能。分析结果表明，FLB细胞可能产生朝向软骨样、关节样、滑膜成纤维细胞样、成骨样细胞分化的轨迹，它们可能共享一个共同的未分化远端芽基祖细胞。而周细胞、内皮细胞、施万细胞、肌源性芽基细胞则显示出各自独立的分化轨迹。此外，还发现了一群罕见的、表达肌腱细胞标记物Scleraxis（SCX）但不表达成纤维细胞标记的细胞，提示肌腱可能由一种罕见的组织驻留祖细胞再生。 * 这些结果支持了两种假说：部分细胞（如肌肉卫星细胞、可能的肌腱祖细胞）来源于罕见的组织驻留祖细胞；而另一些细胞（如施万细胞、成纤维细胞）则可能通过去分化形成祖细胞池。特别是FLB细胞表现出最高的多能性，可能形成再生肢体的关节、软骨和骨骼。

第五， 研究结论与意义 本研究通过单细胞RNA测序技术，首次在单细胞分辨率上系统绘制了成年钝口螈肢体再生过程中芽基生态位的转录组图谱，为理解肢体再生这一复杂过程建立了一个详尽的分子框架。

科学价值： 1. 资源价值： 该研究提供了一个宝贵的分子资源库，详细定义了参与钝口螈肢体再生的细胞群体，并提供了区分这些群体的分子标记。这将极大地方便未来通过谱系追踪、基因功能操纵等分子技术来精确评估各细胞群体在再生中的贡献。 2. 机制见解： 研究揭示了再生过程的细胞复杂性，挑战了将芽基视为均质细胞团的传统观点。它表明再生是多种细胞类型（包括表皮细胞、多种间充质祖细胞、免疫细胞、神经血管相关细胞）协同作用的结果。对表皮分化轨迹、免疫细胞组成、以及间充质细胞分化路径的解析，为深入研究这些细胞间相互作用如何指导再生提供了具体靶点和方向。 3. 理论支持： 研究结果同时支持了肢体再生的“罕见祖细胞”和“去分化”两种来源假说，并指出成纤维细胞样细胞是形成再生骨骼结构的主要多能性来源，而其他谱系（肌肉、肌腱、神经、血管）的细胞可能更多受到其原始身份的谱系限制。研究还发现许多细胞群体重新表达了发育相关的转录因子（如SIX1, EYA1, SCX, SOX2, FOXD3），强烈提示肢体再生在某种程度上重演了发育程序。

应用价值： 对具有超强再生能力的动物模型进行深入研究，有助于识别在哺乳动物（包括人类）中可能被抑制或丢失的关键再生通路和细胞行为。本研究鉴定的关键细胞类型、特异性标记基因以及动态变化的信号通路，为在哺乳动物系统中尝试激活类似再生程序提供了潜在的干预靶点。

第六， 研究亮点 1. 技术先进性与系统性： 这是首次在成年钝口螈肢体再生这一经典模型中应用高通量单细胞RNA测序技术，对超过25,000个细胞进行无偏倚的系统性分析，覆盖了稳态到多个关键再生阶段，数据全面且具有时间维度。 2. 发现新颖性： 研究不仅验证了已知细胞类型（如基底伤口表皮、巨噬细胞、Pax7+肌肉卫星细胞），更发现了许多之前未被充分认识的细胞群体，如伤口表皮中的Leydig细胞、多种表皮亚型、芽基中丰富的T细胞和B细胞等适应性免疫细胞、以及具有不同分化潜能的成纤维细胞样细胞亚群。 3. 深入解析细胞状态与命运： 通过拟时序和轨迹分析，动态地描绘了表皮细胞的分化路径以及芽基中间充质细胞可能的分化轨迹，将静态的细胞类型分类提升到了动态的细胞状态转变研究，提出了FLB细胞多能性分化的具体假设路径。 4. 数据验证严谨： 研究不仅依赖生物信息学分析，还通过大量的RNA原位杂交实验对关键发现进行了空间定位验证，增强了结论的可靠性。 5. 为领域提供新视角与工具： 该研究彻底改变了人们对再生芽基细胞构成复杂性的认识，强调了非芽基细胞（如免疫细胞、伤口表皮）在构建再生微环境中的重要性。所提供的分子标记和细胞图谱是未来功能研究的基石。

第七， 其他有价值的内容 本研究还坦诚讨论了单细胞RNA测序技术的局限性，例如细胞捕获效率、转录本捕获深度可能导致的某些低表达基因（如KAZALDI1）未被检测到。作者指出，随着钝口螈基因组和转录组注释的不断完善，以及计算方法的进步，重新分析该数据集可能会获得更精细的见解。此外，研究将所有单细胞RNA-seq数据和细胞-基因矩阵公开于GEO数据库（GSE121737），并将分析代码开源，体现了研究的可重复性和对科学社区的贡献。