糖皮质激素受体通过Map7d2与Erk2信号协同调控成骨细胞分化

《Frontiers in Immunology》：Glucocorticoid receptor regulates osteoblast differentiation in association with Map7d2 and Erk2 signaling

　　

糖皮质激素受体（glucocorticoid receptor，GR）是骨骼稳态的重要调控因子，但其调控成骨分化的分子机制尚未完全阐明。本研究利用两种互补的原代小鼠成骨细胞模型——骨髓基质细胞（bone marrow stromal cells，BMSCs）和骨来源原代成骨细胞，结合比较转录组学分析，探究了GR的作用。在两种细胞中，GR缺失均显著损害成骨分化，表现为碱性磷酸酶（alkaline phosphatase，Alp）活性降低及成骨标志物表达改变。RNA测序（RNA sequencing，RNA-seq）显示，GR敲低后BMSCs中有3,530个差异表达基因（differentially expressed genes，DEGs），原代成骨细胞中有4,657个DEGs。尽管两种成骨细胞模型的基因表达谱存在显著差异，但通路富集分析揭示了共同的生物学过程，包括MAPK信号传导和细胞迁移。比较转录组学分析进一步鉴定出仅17个在两种细胞中同向一致调控的基因，其中Map7d2成为一个此前未被识别的GR响应候选基因。与单独GR缺失相比，GR与Map7d2联合缺失可部分恢复细胞Alp活性及部分成骨标志物的表达。蛋白质印迹分析进一步显示，联合缺失GR与Map7d2后，Erk2磷酸化及Alp蛋白表达得到恢复，表明部分成骨读出的恢复与Erk2信号的恢复相关联。综上，本研究鉴定出Map7d2是一个与Erk2信号及成骨分化相关的、此前未被识别的GR响应候选基因，揭示了两类原代小鼠成骨细胞模型对GR缺失的共同分子反应，并为GR调控成骨分化的转录及信号网络提供了新见解。

论文解读：糖皮质激素受体通过Map7d2与Erk2信号协同调控成骨细胞分化

一、研究背景与科学问题

骨骼是动态组织，终身通过骨形成与骨吸收的协调活动维持结构完整、矿物质稳态和力学强度。成骨细胞、破骨细胞和骨细胞之间的平衡一旦被打破，可导致骨质疏松、炎性骨丢失和骨折愈合障碍等多种骨骼疾病。成骨细胞来源于骨髓和骨膜等不同骨骼微环境中的异质性骨骼干/祖细胞（skeletal stem and progenitor cell，SSPC）群体。成骨分化由复杂的转录和信号网络调控，涉及Runx2、Sp7及AP-1家族等谱系定义转录因子，以及Wnt/β-catenin、BMP/TGF-β、Hedgehog、Notch、PI3K/AKT和MAPK/ERK等信号通路。

糖皮质激素受体（glucocorticoid receptor，GR）由Nr3c1基因编码，属于核受体超家族，是配体依赖性转录调控因子，参与细胞增殖、分化、代谢、细胞外基质组织和炎症反应等多种生物学过程。大量证据表明糖皮质激素和GR调控成骨细胞生成、成骨细胞增殖与分化、成骨细胞和骨细胞凋亡、Wnt及细胞因子信号以及骨形成。遗传学研究进一步揭示了GR在成骨谱系细胞中的直接和细胞特异性功能。然而，GR调控成骨分化的分子机制仍不完全清楚。为此，研究人员利用两种独立的小鼠原代成骨细胞模型——骨髓基质细胞（BMSCs）和长骨来源原代成骨细胞，结合功能分析与比较转录组学，旨在鉴定GR调控成骨分化中稳健、共享的调控通路和GR响应候选基因。

二、主要技术方法

研究人员从12周龄小鼠长骨中分离原代成骨细胞和BMSCs。采用siRNA SMARTpool技术分别敲低GR和Map7d2，进行反向转染后诱导成骨分化。RNA-seq在Illumina HiSeq 2500平台完成，使用DESeq2进行差异表达分析，DEGs定义为|log₂FC|≥1且FDR<0.01，通路富集分析使用Enrichr和Metascape。蛋白质印迹检测Erk1/2磷酸化及Alp蛋白表达。定量Alp活性以pNPP底物法测定并归一化至细胞活力。

三、研究结果

3.1 GR缺失损害BMSCs的成骨分化并改变转录程序

在BMSCs中，GR敲低显著增加细胞活力，但显著降低细胞Alp活性，定性Alp染色证实成骨分化受损。GR缺失显著降低Sp7和Alpl表达，而Runx2、Col1a1和Bglap无显著变化。RNA-seq鉴定出3,530个DEGs，其中2,110个上调、1,420个下调。通路富集分析显示细胞外结构组织、细胞外基质组织、细胞增殖正调控、细胞迁移调控、骨骼系统发育及MAPK级联正调控等通路显著富集。

3.2 GR缺失损害原代成骨细胞分化并改变转录程序

在原代成骨细胞中，GR敲低显著降低细胞活力、Alp活性，并显著降低Runx2、Sp7、Alpl和Col1a1表达，Bglap无显著变化。RNA-seq鉴定出4,657个DEGs，其中2,004个上调、2,653个下调。通路富集分析显示MAPK级联正调控、PI3K/PKB信号转导正调控、细胞迁移调控、ERK1和ERK2级联调控等通路显著富集。

3.3 比较转录组学鉴定BMSCs和原代成骨细胞中共同的GR调控基因

两个数据集中有2,728个基因重叠，但仅17个基因同向调控。其中Uchl1、Map7d2和Sema3e经qRT-PCR验证在两种细胞中均显著上调。Map7d2因在骨骼生物学中的作用此前未被研究而被优先选为功能研究对象。

3.4 联合缺失GR和Map7d2部分恢复成骨分化

在BMSCs中，GR敲低增加细胞活力并降低Alp活性；Map7d2敲低不影响细胞活力但同样降低Alp活性；联合敲低部分恢复Alp活性。基因表达层面，联合敲低显著恢复Sp7、Alpl和Bglap表达，但Runx2和Col1a1未显著恢复，提示部分且标志物特异性的挽救效应。

3.5 Map7d2与GR缺失BMSCs中Erk2信号改变相关

蛋白质印迹显示，GR敲低降低Erk2磷酸化而不影响总Erk2蛋白表达，Map7d2敲低增加Erk2磷酸化，联合敲低恢复Erk2磷酸化。Alp蛋白表达在GR敲低后显著降低，联合敲低后恢复。Erk1磷酸化变化不明显。

四、讨论与结论

本研究证实GR在调控成骨分化中的重要性，并鉴定出GR、Map7d2和Erk2信号之间此前未被识别的关联。GR缺失在两种细胞模型中均损害成骨分化，但细胞活力反应相反，提示成骨表型不能仅由细胞活力变化解释。超过2,700个基因在两种模型中重叠差异表达，但仅17个同向调控，反映模型间的生物学和技术差异。通路富集分析在MAPK信号和细胞迁移上收敛，支持聚焦共享通路和共同调控基因的研究策略。

Map7d2属于MAP7微管相关蛋白家族，其在骨骼生物学中的作用此前未被研究。联合缺失GR和Map7d2部分恢复Alp活性和部分成骨标志物表达，但挽救不完整且标志物特异性。Map7d2单独敲低降低Alp活性却增加Alpl和Bglap mRNA表达，提示转录本丰度与功能性Alp活性之间存在复杂关系。GR缺失与Erk2磷酸化降低相关，联合缺失恢复Erk2磷酸化及Alp蛋白表达，支持部分挽救与Erk2信号恢复的关联。

研究局限性包括：结论基于两种体外原代细胞模型，需在条件性GR缺陷小鼠和人类成骨细胞中验证；转录组学无法定义染色质水平机制或确立Map7d2为直接GR靶基因；Erk2磷酸化变化虽在两次独立实验中可重复，但样本量有限，未进行推断性统计分析。

研究结论：GR是两种独立原代小鼠细胞模型中高效成骨分化所必需的，研究鉴定出这些系统中共同的GR依赖性转录和信号反应。功能实验表明，与单独GR缺失相比，联合缺失GR和Map7d2部分恢复部分成骨读出，蛋白质印迹结果将该部分恢复与Erk2磷酸化和Alp蛋白表达的恢复相关联。这些发现将Map7d2鉴定为GR依赖性调控网络中参与成骨分化的潜在组分，为深入探究GR、Map7d2和Erk信号之间的分子机制奠定了基础。