1. 领域背景与文献
文献英文标题：Single-cell RNA sequencing reveals the expression characteristics and potential associated pathways of hemoglobin in chondrocytes of osteoarthritis；发表期刊：BMC Molecular and Cell Biology；影响因子：5.3（2023年）；研究领域：骨关节炎发病机制与软骨代谢调控。
领域共识：骨关节炎是全球最高发的退行性关节疾病，随着人口老龄化、肥胖率上升发病率逐年升高，目前临床缺乏早期特异性诊断标志物，现有治疗手段仅能缓解症状无法阻断病程进展，终末期患者只能接受关节置换手术，造成巨大的医疗与社会负担。正常软骨组织无血管分布，软骨细胞长期处于1%-6%氧浓度的低氧微环境中，既往研究认为软骨细胞主要依赖缺氧诱导因子（HIF）调控的糖酵解途径供能以适应低氧环境。近年研究首次发现软骨细胞可表达血红蛋白形成无膜复合体“hedy”，实现局部氧存储与供给，是软骨细胞低氧适应的新机制，但该机制的具体调控模式、在骨关节炎进展中的作用尚未明确，且不同功能亚群软骨细胞的血红蛋白表达谱存在研究空白。
本研究针对上述领域未解决问题，基于单细胞转录组测序技术解析正常人与骨关节炎患者软骨细胞的血红蛋白表达特征，明确其调控通路及相关分子靶点，为骨关节炎的发病机制解析与治疗新策略开发提供学术依据。
2. 文献综述解析
本研究的文献综述按领域认知迭代逻辑展开，首先梳理哺乳动物氧转运的经典研究脉络，明确红细胞中的血红蛋白负责全身氧运输，肌细胞、神经元可分别表达肌红蛋白、神经红蛋白实现局部氧储备，满足高氧消耗生理活动的需求，为软骨细胞血红蛋白的功能解读奠定理论基础。
随后作者梳理软骨低氧适应的研究进展，明确既往研究已证实软骨细胞主要通过HIF通路调控糖酵解供能，近年新发现的软骨细胞血红蛋白低氧适应机制尚未解析表达模式；同时作者系统总结软骨细胞异质性的研究成果，明确目前已鉴定出增殖性软骨细胞、前肥大软骨细胞、纤维软骨细胞、稳态软骨细胞等7种功能亚群，各亚群在骨关节炎进展中分别承担增殖、基质合成、纤维化调控等不同功能，但不同亚群的血红蛋白表达差异尚无相关报道。
最后作者梳理SPP1编码的骨桥蛋白的研究现状，明确骨桥蛋白在骨关节炎组织中表达显著升高，高表达SPP1的软骨细胞具有更强的血管生成能力与衰老表型，参与骨关节炎进展，但SPP1与软骨细胞血红蛋白表达的关联尚未被揭示。现有研究的优势在于已初步证实软骨细胞存在血红蛋白表达、明确了软骨细胞亚群的功能异质性，局限性在于未解析血红蛋白在不同软骨细胞亚群的分布特征、未明确其在骨关节炎中的调控机制，本研究的创新点在于首次在单细胞分辨率下解析骨关节炎与正常软骨中血红蛋白的表达谱，明确了其相关调控通路及与SPP1的关联，填补了领域研究空白。
3. 研究思路总结与详细解析
本研究的整体目标是明确软骨细胞中血红蛋白的表达特征及在骨关节炎进展中的潜在作用，核心科学问题为血红蛋白在不同软骨细胞亚群的分布模式、高表达血红蛋白软骨细胞的功能特征及上游调控机制，技术路线形成“公共数据挖掘→生物信息学分析→体内实验验证”的完整闭环：首先下载公开单细胞转录组数据集完成质控与细胞注释，随后分析血红蛋白在不同软骨亚群的表达特征，再对比血红蛋白高低表达细胞的差异基因与通路、细胞间通讯特征，明确潜在调控分子SPP1，最后通过动物模型验证两者的表达关联，推导功能机制。
3.1 单细胞数据获取与预处理
实验目的是获取合格的正常人与骨关节炎患者软骨单细胞转录组数据，为后续分析奠定基础。方法细节为从基因表达综合数据库（GEO）下载GSE169454数据集，包含3例正常膝关节软骨样本、4例骨关节炎膝关节软骨样本的10× Genomics单细胞测序数据，使用R语言Seurat包进行质控，过滤线粒体基因占比≥15%、基因表达数<200的低质量细胞，对合格细胞进行对数标准化、主成分分析，使用harmony包去除不同样本的批次效应。结果解读为质控后共获得17个细胞亚群，去除红细胞、免疫细胞、成纤维细胞等非软骨细胞群，基于SOX9、COL2A1、ACAN等软骨细胞标记基因，鉴定得到对照组13个软骨细胞亚群、骨关节炎组13个软骨细胞亚群，初步可视化显示HBB、HBA1、HBA2等血红蛋白编码基因在骨关节炎组软骨细胞中的表达占比显著高于对照组，对应结果见图1。
文献未提及该环节具体实验产品，领域常规使用R语言相关生物信息学分析包、公共基因表达数据库完成该环节分析。
3.2 软骨细胞亚群注释与血红蛋白表达分布分析
实验目的是明确血红蛋白在不同软骨细胞亚群中的表达特征。方法细节为使用t分布随机邻域嵌入（t-SNE）对软骨细胞进行无监督分群，基于已报道的标记基因注释得到调控性软骨细胞、稳态软骨细胞、增殖性软骨细胞、前肥大软骨细胞、效应软骨细胞、纤维软骨细胞、肥大软骨细胞共7个亚群，统计各亚群中HBB高表达细胞的占比。结果解读为对照组中调控性软骨细胞占比最高，骨关节炎组中纤维软骨细胞、前肥大软骨细胞占比显著升高；骨关节炎组软骨细胞中HBB、HBA1、HBA2的表达水平显著高于对照组；HBB在增殖性软骨细胞、纤维软骨细胞、前肥大软骨细胞中的表达水平最高，分别为0.44、0.42、0.36，提示这三类亚群是骨关节炎中血红蛋白的主要表达细胞，对应结果见图2。
文献未提及该环节具体实验产品，领域常规使用R语言相关生物信息学分析包完成该环节分析。
3.3 血红蛋白高低表达软骨细胞的差异分析
实验目的是解析高表达血红蛋白的软骨细胞的功能特征。方法细节为将HBB表达log2(FC)>1的细胞定义为HBB+软骨细胞，其余为HBB-软骨细胞，使用FindMarkers函数筛选两组差异表达基因，分别进行京都基因与基因组百科全书（KEGG）通路富集分析、基因集富集分析（GSEA）。结果解读为共筛选得到5236个差异表达基因，其中HBB+细胞中2694个基因上调、2542个基因下调；HBB+细胞中血红蛋白组成、核糖体构建、蛋白合成相关基因显著上调，凋亡、程序性细胞死亡、蛋白泛素化、细胞连接相关基因显著下调；KEGG富集显示HBB+细胞中富集HIF-1、TNF、mTOR、凋亡通路，HBB-细胞中富集氧化磷酸化、活性氧、TGF-β、铁死亡通路，提示两类细胞的生理功能存在显著差异，对应结果见图3。
文献未提及该环节具体实验产品，领域常规使用R语言相关生物信息学分析包完成该环节分析。
3.4 细胞间通讯分析
实验目的是解析高表达血红蛋白的软骨细胞的细胞间交互特征。方法细节为使用CellChat工具分析HBB+、HBB-细胞群中各软骨亚群的通讯数量、强度及相关信号通路。结果解读为HBB+细胞群的细胞通讯数量与HBB-群相近，但通讯强度略低；两组中胶原通路均为最主要的细胞互作通路，但HBB+群中SPP1通路的表达强度显著高于HBB-群；SPP1相关信号主要由调控性软骨细胞分泌，被包括前肥大软骨细胞、增殖性软骨细胞在内的所有亚群接收，对应结果见图4。
文献未提及该环节具体实验产品，领域常规使用R语言CellChat包完成该环节分析。
3.5 SPP1与血红蛋白的关联验证及动物实验验证
实验目的是明确SPP1与软骨细胞血红蛋白表达的关联及在体内的表达特征。方法细节为筛选SPP1高表达（SPP1+）与低表达（SPP1-）软骨细胞的差异基因，进行通路富集分析；构建内侧半月板不稳定（DMM）诱导的8周龄SPF级C57BL/6雄性小鼠骨关节炎模型，每组n=6，造模8周后取关节组织进行冰冻切片，开展免疫荧光染色检测血红蛋白、SPP1的表达。结果解读为SPP1+细胞主要分布于骨关节炎组，在调控性软骨细胞、增殖性软骨细胞、前肥大软骨细胞中的占比最高，分别为0.78、0.44、0.43；SPP1+细胞中HBB、HBA1、HBA2等血红蛋白编码基因的表达水平显著升高，mTOR信号通路、胰岛素信号通路显著富集，mTOR通路相关分子MTOR、IGFBP1、EIF4EBP2的表达显著上调；动物实验显示骨关节炎组关节软骨中血红蛋白阳性细胞占比为42.3±3.6%（n=6，P<0.0001），SPP1阳性细胞占比为38.7±4.1%（n=6，P<0.01），均显著高于假手术组，对应结果见图5。
实验所用关键产品：Proteintech的SPP1抗体、Abcam的抗兔Alexa荧光偶联二抗、Invitrogen的4",6-二脒基-2-苯基吲哚（DAPI）、碧云天的免疫染色封闭液、徕卡的冰冻切片机、蔡司的荧光显微镜。
4. Biomarker研究及发现成果
本研究涉及的潜在Biomarker包括软骨细胞表达的血红蛋白（HBB、HBA1、HBA2）及SPP1编码的骨桥蛋白，筛选验证逻辑完整：首先基于公共单细胞数据集筛选骨关节炎与正常软骨的差异表达分子，随后验证其在不同软骨功能亚群的表达特异性，再通过差异分析、细胞间通讯分析明确两者的调控关联，最后通过动物模型体内验证表达一致性。
两类Biomarker的检测样本包括人膝关节软骨单细胞测序样本、小鼠膝关节软骨组织切片，验证方法涵盖单细胞转录组测序定量、免疫荧光染色半定量。统计学数据显示骨关节炎组软骨细胞HBB表达水平较对照组显著升高（P<0.01），SPP1+细胞中血红蛋白编码基因的表达较SPP1-细胞显著升高（P<0.001），动物实验中骨关节炎组血红蛋白阳性细胞占比为42.3±3.6%（n=6，P<0.0001），SPP1阳性细胞占比为38.7±4.1%（n=6，P<0.01），均显著高于假手术组，文献未明确报道两类分子作为诊断标志物的ROC曲线、敏感性特异性数据。
核心成果显示，软骨细胞表达的血红蛋白可作为骨关节炎进展中软骨细胞低氧适应状态的潜在标志物，高表达血红蛋白的软骨细胞主要为增殖性软骨细胞、纤维软骨细胞、前肥大软骨细胞，富集HIF-1、mTOR通路，提示其与骨关节炎进展中软骨细胞的增殖、基质合成过程相关；SPP1可作为高表达血红蛋白软骨细胞的上游调控标志物，其通过调控mTOR通路参与血红蛋白表达的调控。本研究的创新性在于首次在单细胞分辨率下明确了软骨细胞血红蛋白的表达谱，首次发现SPP1与软骨细胞血红蛋白表达的关联，推测两者的调控轴可能成为骨关节炎治疗的新靶点，相关机制需后续细胞功能实验进一步验证。