《Scientific Reports》：Identification of biomarkers associated with endoplasmic reticulum stress-related cell death in osteoporosis based on bulk and single-cell transcriptomic analyses and experimental validation
编辑推荐：
骨质疏松(OP)是一种复杂的骨骼疾病，其发病机制尚未完全阐明。本研究通过整合转录组、单细胞测序与实验验证，聚焦内质网应激(ER Stress)与程序性细胞死亡(PCD)的关联，致力于发现OP的新型诊断生物标志物。研究人员通过机器学习筛选出CAMKK2与DAPK3这两个关键基因，并证实其在OP中显著下调。该研究为理解OP的细胞死亡机制、开发新的诊断工具和潜在治疗靶点提供了新思路。
随着年龄增长，骨骼就像一栋年久失修的建筑，内部支撑结构变得脆弱，这正是骨质疏松症带来的困扰。骨质疏松是一种以骨量减少、骨微结构破坏为特征的代谢性骨病，患者骨折风险显著增高，严重影响生活质量。尽管医学界对其已有认识，但其确切的发病机制仍如一团迷雾，尤其是细胞层面的精细调控过程尚不清晰。近年来，科学家们将目光投向了细胞内一个名为“内质网”的细胞器。内质网是蛋白质合成与加工的重要车间，当它功能紊乱时，会引发“内质网应激”，这可能是触发骨细胞（如成骨细胞、破骨细胞）走向程序性细胞死亡(PCD)的关键开关。这种异常的细胞死亡，被认为是导致骨形成减少、骨吸收失衡，最终发展为骨质疏松的重要推手。那么，能否找到与内质网应激相关的细胞死亡密切关联的关键分子，作为诊断骨质疏松的“信号灯”和治疗的新靶点呢？这正是发表在《Scientific Reports》上的一项研究所要回答的核心问题。
为了揭开这个谜团，研究人员采取了一套从大数据挖掘到实验验证的组合拳。他们首先从公共数据库获取了骨质疏松患者和健康对照者的转录组数据作为训练集(GSE56815)和验证集(GSE56814)，以及一份单细胞RNA测序(scRNA-seq)数据(GSE147287)。通过生物信息学分析，他们筛选出差异表达基因(DEGs)，并将其与已知的内质网应激相关基因、程序性细胞死亡相关基因取交集，得到了候选基因池。随后，利用机器学习方法进一步筛选，并通过表达水平分析确定了核心的生物标志物。为了深入探究这些标志物的功能，研究团队还进行了基因集富集分析(GSEA)、免疫细胞浸润分析、药物预测与分子对接。更重要的是，他们利用单细胞测序数据，深入骨髓微环境，识别了关键细胞类型，分析了细胞间的通信网络，并追踪了标志物基因在细胞分化过程中的动态表达变化。最后，通过逆转录定量聚合酶链反应(RT-qPCR)实验，在独立样本中对生物信息学的发现进行了验证。
1. 差异表达分析与候选基因鉴定
研究人员通过比较骨质疏松患者与对照组的基因表达谱，识别出大量的差异表达基因(DEGs)。将这些DEGs分别与内质网应激相关基因集、程序性细胞死亡相关基因集取交集后，最终得到了28个共同的候选基因。这28个基因位于内质网应激与细胞死亡通路的交叉点，是潜在的枢纽分子。
2. 机器学习筛选与生物标志物确认
为了从28个候选基因中找出最具有诊断价值的核心基因，研究者运用了机器学习算法进行筛选。结合表达水平分析，最终将目标锁定在两个基因上：钙/钙调素依赖的蛋白激酶激酶2(CAMKK2)和死亡相关蛋白激酶3(DAPK3)。分析显示，这两个基因在骨质疏松患者的骨组织样本中表达水平 consistently（一致地）下调。
3. 生物标志物的功能与通路富集分析
通过基因集富集分析(GSEA)，研究发现CAMKK2和DAPK3的表达与“细胞因子-细胞因子受体相互作用”等信号通路显著相关。这提示这两个生物标志物可能通过调节免疫和炎症微环境参与骨质疏松的病理过程。进一步的免疫细胞浸润分析发现，这两个基因的表达水平与骨髓中活化的树突状细胞的比例存在显著相关性，暗示它们可能影响骨免疫平衡。
4. 潜在药物预测
基于计算模拟的药物预测和分子对接分析表明，一些已知药物，如骨化三醇（calcitriol，一种活性维生素D）和达那唑（danazol），可能与CAMKK2或DAPK3蛋白产生稳定的相互作用。这为老药新用或基于新靶点的药物研发提供了初步的计算依据。
5. 单细胞层面的深入解析
研究并未止步于组织水平的分析。利用单细胞RNA测序(scRNA-seq)数据，研究人员深入剖析了骨髓微环境的细胞组成。他们发现，骨髓间充质干细胞(BM-MSCs)是骨髓中一类非常重要的细胞群体，并且其基因表达特征与骨质疏松状态密切相关，因此被确定为潜在的关键细胞。细胞通信分析进一步揭示，BM-MSCs之间存在复杂的相互作用网络，例如由血管生成素样蛋白4(ANGPTL4)和钙粘蛋白11(CDH11)介导的BM-MSC自我通信。这表明BM-MSCs群体内部通过特定信号分子进行自我调控，可能影响其成骨分化和功能。
6. 生物标志物在细胞分化中的动态表达
为了了解CAMKK2和DAPK3在骨形成过程中的作用，研究者利用拟时间序列分析追踪了这两个基因在BM-MSCs分化过程中的表达变化。结果发现，CAMKK2的表达水平随着分化阶段而波动；而DAPK3的表达则呈现出“高-低-高”的动态模式。这种动态变化提示它们可能精确调控着骨髓间充质干细胞向成骨细胞分化的不同阶段。
7. 实验验证
所有计算分析的结果最终都需要实验的证实。研究人员通过逆转录定量聚合酶链反应(RT-qPCR)实验，在独立的样本集中检测了CAMKK2和DAPK3的表达。实验结果与生物信息学分析结果一致，确认了这两个基因在骨质疏松样本中确实存在表达下调，从而在实验层面验证了它们作为生物诊断标志物的可靠性。
综上所述，本研究通过整合多组学数据与实验验证，系统性地探索了骨质疏松症中内质网应激与程序性细胞死亡的关联。研究成功鉴定出CAMKK2和DAPK3是两个与内质网应激相关细胞死亡密切关联的诊断性生物标志物，它们在骨质疏松中表达下调。功能分析表明它们可能通过影响细胞因子通路和免疫微环境（特别是活化的树突状细胞）发挥作用。计算药物预测提示了潜在的干预方向。在单细胞层面，研究将骨髓间充质干细胞(BM-MSCs)锁定为关键细胞，并揭示了其内部复杂的通信网络，以及生物标志物在其分化过程中的动态表达模式。这项研究的意义在于，它不仅为骨质疏松提供了新的、具有潜在应用价值的诊断标志物（CAMKK2和DAPK3），更重要的是，它将疾病的机制研究从组织层面深入到单细胞与细胞间通信层面，揭示了内质网应激通过调控特定细胞死亡通路影响骨髓间充质干细胞功能，从而参与骨质疏松发生发展的新机制。这为未来开发针对该通路的新型治疗策略（如调节内质网应激、靶向CAMKK2/DAPK3或干预BM-MSCs通信）奠定了重要的理论基础。
