近日，美国宾夕法尼亚大学Tan Kai、秦聆团队合作在Cell上发表了题为“Mapping the cellular biogeography of human bone marrow niches using single-cell transcriptomics and proteomic imaging”的文章。研究团队对29,325个非造血细胞进行单细胞RNA测序（scRNA-seq）分析，发现了9种不同转录亚型；同时分析了53,417个造血细胞，预测其与非造血细胞亚型的相互作用；并通过索引共检测（CODEX）对超120万个细胞进行空间分析，整合CODEX数据与scRNA-seq数据，将预测的细胞信号与空间邻近性相关联。

该研究发现了多种非造血细胞亚型，揭示了人类骨髓中非造血细胞和造血细胞的重要性和相互作用，以及早期骨髓造血过程中特定细胞类型的空间分布。利用CODEX图谱对新图像进行注释，发现了急性髓系白血病（AML）患者样本中骨髓MSC扩增，以及白血病母细胞与MSC共同富集的空间邻域。综上，这一空间分辨的骨髓多组学图谱为研究人类骨髓细胞相互作用提供了参考。

为揭示人类骨髓的细胞组成，研究团队开发了利用酶从股骨头组织中释放细胞的实验流程，并设计了一种富集策略。通过从16名接受髋关节置换手术的个体中收集新鲜股骨头的骨髓，研究人员利用scRNA-seq和CODEX共分析了53,417个造血细胞和29,325个非造血细胞，构建了骨髓综合scRNA-seq图谱。经质量过滤后，每个细胞的中位数为16,903个唯一分子标识符（UMIs）、3,117个独特基因和约3%线粒体reads。

造血细胞分析发现了三系造血和HSPC亚群，并改进了其他以造血为中心的骨髓scRNA-seq数据，捕获了完整的粒细胞分化轨迹。非造血细胞分析鉴定了三种主要细胞系——血管平滑肌细胞 （VSMC）、EC和间充质细胞（包括MSCs和成骨细胞）。骨髓中这些细胞内的细胞多样性远高于现有scRNA-seq参考数据集。

MSCs细胞存在显著异质性，研究团队进一步确定了其可分为骨细胞、脂肪细胞和成纤维细胞等不同类型。此外，研究还发现了多种在基因表达、细胞功能和分化能力上存在显著差异的独特MSCs亚群（图3），具有最高CXCL12水平和高表达脂肪系基因的Adipo-MSCs亚群；高表达脂肪系基因的THY1+ MSCs亚群，其具有独特的表达谱；以及Fibro-MSCs亚群，具有较高的成纤维细胞集落形成单位（CFU-F）形成能力和分化能力，可表达MSCs标志物（即THY1、DPT、NES和NT5E等）。

研究发现了两类骨髓内皮细胞，即动脉内皮细胞（AEC）和静脉内皮细胞（SEC）。将多个已发表的健康人类MSCs scRNA-seq数据集映射到骨髓综合scRNA-seq图谱中，发现该图谱捕获了更多非Adipo-MSC亚群，以及THY1+ MSCs细胞和Fibro-MSCs细胞，这强调了捕获非造血细胞异质性的重要性。

综上，该研究发现了人类骨髓MSCs中存在先前未被描述、功能上具有影响的多样性，并证明了骨髓综合scRNA-seq图谱的实用性，可用于表征其他scRNA-seq数据集中MSCs的相似性和差异性。

研究团队分析了MSC的异质性能否导致造血支持因子的转录差异（图4），发现不同MSC类型专门产生特定的支持因子，如CXCL12在Adipo-和THY1+ MSCs中表达最高，而IGF1主要在Fibro-MSCs中表达。

CellChat系统分析显示，成骨细胞和MSCs都可能维持HSPC，但大多数典型支持因子由MSCs亚群产生，尤其是Adipo-MSCs和THY1+MSCs；CD4+T细胞是骨髓中FLT3LG的主要生产者，其对骨髓形成十分重要。

为对scRNA-seq图谱中骨髓细胞类型进行单细胞定义，研究团队通过全玻片CODEX多重成像技术对12个样本进行成像，并设计验证了一个54-plex panel（53个抗体+ DAPI）（图5）。共识别出32种细胞类型，包括免疫表型HSC（Lin- CD34+ CD38- CD45RA- CD90+）和Schwann细胞等罕见细胞；个体之间的总体细胞类型分布相似。

CODEX细胞类型的蛋白质表达谱与scRNA-seq数据高度匹配，匹配的细胞类型得分非常高，最终细胞类型注释可用于“细胞表型图”中细胞拓扑结构可视化。该注释参考图是造血细胞在其微环境背景下的拓扑综合图谱，凸显了这两个图谱在了解人类骨髓细胞解剖学方面的互补性。

随后，研究团队通过细胞邻域分析和CellChat预测了骨髓内的细胞组织和相互作用。结果显示，共识别出15个细胞邻域（CNs），混合谱系区域（CN7 和 CN8）HSPCs最多，CN4和CN6两个邻域富含GMP/成髓细胞和早期髓系祖细胞（EMP），其沿骨内表面和小动脉周围强烈定位。CellChat预测显示，成骨细胞-VSMC相互作用并共定位，形成动脉-骨内膜生态位。

研究团队分析了6种关键微环境结构，包括脂肪细胞、小动脉、骨骼、巨噬细胞及其突起、窦状体和CXCL12+基质（图6）。结果显示，小动脉最接近骨骼；中心生态位以窦状体、脂肪细胞、基质和巨噬细胞为特征。有趣的是，所有MSCs细胞都靠近动脉-骨内膜生态位，但Adipo-MSC位于窦状体附近。

CODEX图谱显示，HSPCs接近巨噬细胞；GMP/髓母细胞以及EMPs接近脂肪细胞。相比之下，成熟的髓样细胞最接近血窦，而非小动脉或骨骼，揭示成熟过程发生在中央、血窦周围生态位，上述数据表明髓系发育受到空间限制。

研究团队以健康图谱为参考，使用骨髓CODEX panel探究了肿瘤进化和微环境的变化，分析了AML患者样本，发现AML患者的髓样细胞比例显著增加，出现非常明显的结构变化。邻域分析显示，AML样本中有4个blast富集的邻域——Dx- CN3 /CN5/CN9/CN12，其中两个邻域富集了MSC细胞群。进一步分析发现动脉-骨内膜早期髓系邻域可能在治疗抵抗中发挥作用。

综上所述，研究团队报告了迄今为止最全面的健康人类骨髓参考图谱，涵盖了造血和非造血细胞类型及其空间组织，揭示了人类骨髓中显著的非造血多样性与原位离散的空间表型，展示了CODEX图谱在AML患者中的应用。该研究构建的人类骨髓多组学图谱将为未来人类骨髓生态位的研究提供重要参考。