IF50.5｜人类胚胎骨骼发育的多组学图谱

单细胞多组学揭示人类胚胎骨与关节形成的细胞状态和调控网络。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

人类骨骼发育始于受孕后6-8周，颅骨通过膜内成骨，四肢滑膜关节通过软骨内成骨。此前对小鼠的研究已鉴定出GDF5+间带细胞、HIC1+成纤维祖细胞等关键群体，但人类胚胎中这些过程的细胞和分子基础仍不清楚。单细胞技术已用于胎儿肢体研究，但缺乏跨区域、多时间点的配对转录与表观数据，且颅骨发育的细胞图谱尚未建立。

科学问题

卡在哪里： 此前人类胚胎骨骼发育研究主要依赖单细胞转录组，缺乏染色质可及性信息，无法解析调控网络；空间定位也局限于少数基因。颅骨不同成骨模式的细胞轨迹和调控机制未被系统描述。此外，骨关节炎等复杂疾病的遗传信号与发育期细胞状态的联系尚未建立。

本文要回答： 构建5-11 PCW人类胚胎关节和颅骨的多组学图谱，解析骨与软骨形成的细胞轨迹、调控网络和空间组织，并连接发育细胞状态与骨关节炎、颅缝早闭等疾病。

技术路线

作者主要用了：snRNA-seq、snATAC-seq、空间转录组、ISS、SCENIC+、CellOracle。

作者对12个供体的髋、膝、肩关节及颅骨前后区域进行配对snRNA-seq和snATAC-seq，获得336,162个高质量液滴。同时利用155重ISS对完整胚胎肢体进行空间定位，并用Visium分析颅缝。通过整合转录与表观数据，使用SCENIC+构建基因调控网络，结合RNA速率和轨迹分析推断分化路径。开发ISS-Patcher工具将液滴数据标签映射到空间数据，并利用OrganAxis定义颅骨成熟轴。最后通过fGWAS和SNP2Cell连接细胞类型与疾病遗传信号。

图 1　Fig. 1。图1展示实验设计：五个解剖区域采样示意图，12个供体信息，以及336,162个液滴的UMAP，显示8个主要细胞区室。间充质细胞在所有区域占主导，肌源性细胞在颅骨缺失。该图建立了整个图谱的细胞分类基础。

核心亮点

亮点 1｜滑膜关节的带状组织

通过亚聚类InterzoneChon群体，鉴定出早期、关节、纤维和GDF5hi等间带亚群。空间定位显示早期间带细胞弥散分布，关节间带富集于未来关节软骨区域，纤维间带富集于肩关节。SCENIC+预测早期间带细胞表达PRRX1，低RUNX2但中等SOX5/6/9，提示向软骨分化倾斜。GDF5hi亚群低成骨和成软骨活性，可能维持未分化状态。这些结果揭示了关节形成初期的空间分区和分子梯度。

图 2　Fig. 2。图2聚焦关节间带分区：RNA速率UMAP显示七个间带亚群的轨迹，SCENIC+热图展示转录因子表达与靶基因可及性。ISS-Patcher工作流示意图及后肢空间图显示早期间带弥散分布，关节间带富集于未来关节软骨，纤维间带富集于肩关节。

亮点 2｜颅缝成骨的调控网络

鉴定出CranialMes、FacialMes、PArchMes等颅骨祖细胞，以及表达TWIST1、ZIC1的SutureMes1/2和HHIP+PreOB。空间映射显示SutureMes位于冠状缝，HHIP+PreOB位于骨化前沿。SCENIC+网络显示LMX1B、TWIST1等抗成骨转录因子通过TCF12、PRRX2间接抑制RUNX2和HHIP，维持缝未骨化状态。OrganAxis分析表明沿成熟轴抗成骨因子下调，促成骨因子上调，形成空间分子开关。

图 3　Fig. 3。图3展示颅骨成骨的空间轨迹：光片显微镜图像显示8.5和10 PCW胚胎颅骨中COL2A1和SP7染色，冠状缝和矢状缝区域清晰。OrganAxis定义的成熟轴显示TWIST1+ SutureMes在缝区富集，HHIP+PreOB在骨化前沿富集，成骨标志物沿轴递增。

亮点 3｜成骨轨迹与血管生成

推断出两条主要成骨轨迹：膜内成骨从CranialMes经SutureMes到HHIP+PreOB；软骨内成骨从LimbMes到前成骨细胞和成骨细胞。空间数据显示LimbMes在7.3 PCW肢体中定位于骨骺，前成骨细胞位于初级骨化中心。在颅骨中，SutureMes和HHIP+PreOB高表达VEGFA/B，与内皮细胞VEGF受体共定位，提示成骨前沿的血管生成耦合。CellPhoneDB预测内皮细胞通过NOTCH、FGF2等信号促进成骨分化。

图 4　Fig. 4。图4展示软骨细胞异质性和施万细胞谱系：UMAP显示软骨细胞亚型，点图展示标记基因表达。施万细胞UMAP显示SOX9+ enSC群体，点图显示其共表达软骨和施万标志。RNA-ISH图像显示MPZ和SOX9在髋软骨内共定位。

亮点 4｜施万细胞作为软骨细胞来源

发现SOX9+ endoneurial Schwann cells (SOX9+ enSC) 共表达软骨细胞标志（SOX9、COL2A1）和施万细胞标志（MPZ、SOX10）。轨迹分析显示其来源于施万细胞前体，表达间充质和HOX基因。RNA-ISH和ISS证实MPZ与SOX9在髋软骨内共定位。该群体可能代表人类软骨细胞的非经典来源，与小鼠中施万细胞分化为软骨细胞的报道一致。

图 5　Fig. 5。图5展示骨关节炎遗传信号富集：fGWAS结果显示膝骨关节炎信号富集于软骨细胞，髋骨关节炎信号富集于成骨细胞。SNP2Cell示意图展示网络传播方法。热图显示前成骨细胞在髋骨关节炎中富集最高，关节软骨细胞在膝骨关节炎中富集。

亮点 5｜骨关节炎遗传信号的细胞类型特异性

fGWAS分析显示膝骨关节炎信号富集于软骨细胞状态，而髋骨关节炎信号富集于前成骨细胞、成骨细胞和骨细胞。SNP2Cell进一步识别出前成骨细胞和关节软骨细胞中的疾病相关子网络，涉及NFATC1、FOXA3、RUNX2等转录因子。这些结果表明髋和膝骨关节炎可能分别通过成骨和成软骨途径在发育期奠定风险基础。

亮点 6｜颅缝早闭基因的扰动模拟

将颅缝早闭相关基因与成骨轨迹交叉分析，发现13/22个基因富集于膜内成骨祖细胞。CellOracle模拟TWIST1、MSX2、LMX1B敲除导致SutureMes2中高速度偏移，方向与成骨一致。网络分析显示这些转录因子相互调控，并与TCF12、ALX4等颅缝早闭基因形成紧密网络，提示冗余调控维持缝未闭合。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： snRNA-seq、snATAC-seq、空间转录组、ISS、SCENIC+、CellOracle
• 公开数据： E-MTAB-14385
• 代码： 公开仓库 http://github.com/emdann/milopy、https://github.com/Teichlab/HCA_Heart_ver2、https://github.com/Teichlab/drug2cell；本地已克隆：PHM、iss_patcher、skeletal_dev_atlas、snp2cell、milopy、microaligner、drug2cell
• 谁该点开原文： 从事骨骼发育、单细胞多组学、空间转录组或骨关节炎遗传学研究的科研人员。

带走一句

发育期单细胞多组学数据结合空间定位和调控网络分析，能揭示疾病遗传信号的细胞类型特异性，为复杂骨骼疾病提供机制线索。