骨缺损是由创伤、感染、骨质疏松或骨肿瘤切除等原因引起的骨科常见疾病
[1-3]。自体骨移植是临床治疗的金标准
[3, 4],其具有良好的骨传导和骨诱导特性,但自体骨的来源有限,且骨的质量具有可变性,例如,部分骨质疏松患者的骨质量难以达到骨缺损填充的标准
[5, 6],此外,取骨可能延长手术时间,并且术后供体部位易出现严重疼痛和并发症
[4, 7]。因此,迫切需要寻求能达到自体骨移植效果的替代治疗方案。骨组织工程
[8-10]作为一种有前景的替代疗法,通过将具有良好骨诱导性及骨传导性的复合生物材料和干细胞植入骨缺损部位,利用支架材料支持细胞生长和细胞外基质沉积
[11],从而加速新骨再生。
间充质干细胞(MSCs)是一类中胚层来源的多能干细胞,主要存在于骨髓、脂肪等结缔组织和器官间质中。它们具有自我更新和多向分化潜能,即能够分化为多种间质组织,如骨、软骨、脂肪、骨髓造血组织等。因其来源广泛且成骨能力优异,MSCs成为骨组织工程研究的“宠儿”
[12]。脂肪间充质干细胞(ADSCs)作为MSCs的重要成员之一,相较于骨髓间充质干细胞,其具有易获取且数量丰富的特点,因此ADSCs被认为是骨组织工程中最适合使用的细胞
[13]。
根据既往研究报道,间充质干细胞成骨和成脂分化命运的决定取决于多种转录因子和信号通路的调控
[14, 15],而成脂分化和成骨分化之间存在负相关关系
[16],即“成脂-成骨”平衡。找到促进“成脂-成骨”平衡向成骨分化倾斜的关键节点因子,能够帮助ADSCs在适宜的培养环境中实现定向的成骨分化。找到推动ADSCs向成骨谱系分化的关键因素,并利用骨组织工程为ADSCs提供适宜的微环境,确保移植至骨缺损部位的ADSCs能向成骨谱系分化并促进骨缺损修复,是本研究的核心及关键。
高迁移率族蛋白A2(HMGA2)是具有109个氨基酸的非组蛋白染色体蛋白。由于包含AT-hook结构域,HMGA2能够结合富含腺嘌呤和胸腺嘧啶的核苷酸序列,从而调控靶基因的转录活性。此外,HMGA2可以通过别构效应促进转录因子与靶区域的结合,间接影响基因表达
[17, 18]。通常认为HMGA2只是调控干细胞成脂分化的重要因子
[19-22],然而,部分研究表明HMGA2也参与了MSCs的成骨分化
[23-25]。我们在前期实验中也发现,
Hmga2在mADSCs成骨分化进程中出现了差异性表达。然而,
Hmga2在ADSCs成骨分化进程中的作用机制尚不明确。
Hmga2的表达水平与干细胞分化命运决定息息相关,目前对于HMGA2在成骨分化中的作用仍存有较大争议。研究指出,LncRNA LOXL1-AS1/miR-196a-5p/
HMGA2信号轴通过调节C/EBPβ介导的PPARγ的表达促进间质/间充质干细胞的成骨分化
[23];miR-664a-5p
[24],let-7
[25],mirR-98
[26]均可直接靶向
Hmga2促进骨髓间充质干细胞或成骨前体细胞的成骨分化;而miR-497-5p
[27]则通过靶向
Hmga2并激活下游的JNK信号通路加速成骨前体细胞的成骨分化进程,以上研究结果均表明,
Hmga2的下调直接或间接地促进成骨分化。但也有研究得出相反的结论,即
Hmga2的下调抑制成骨细胞谱系的成骨分化
[28-30]。本研究通过生物信息学分析和文献调研,结合前期实验结果,提出如下科学假设:HMGA2不仅参与了ADSCs的成脂分化,也参与了其成骨分化,并在ADSCs“成脂-成骨”分化平衡中起关键作用,决定其分化命运。因此,阐明HMGA2对ADSCs的调控作用及其分子机制,有助于完善ADSCs的成骨分化机制,找到推动ADSCs向成骨谱系分化的关键因子,并有望为ADSCs治疗骨缺损的骨组织工程应用提供新的思路。
1 材料和方法
1.1 材料
I型胶原酶(Worthington Biochemical),间充质干细胞专用胎牛血清(Vivacell),β-甘油磷酸(Sigma),维生素C(Sigma),地塞米松(Sigma),DMEM低糖培养基及Opti-MEM减血清培养基(Gibco),青霉素-链霉素溶液(Invitrogen),碱性磷酸酶染液及定量检测试剂盒(碧云天),茜素红溶液(赛维尔生物),全细胞蛋白裂解液及BCA定量试剂盒(凯基生物),Masson三色染色液(索莱宝),苏木素伊红(HE)染色试剂盒(雷根生物),Epson Perfection V600 扫描仪(EPSON),化学发光(ECL)发光液(爱必胜生物),Lipofectamine 3000脂质体(Invitrogen),逆转录及实时定量PCR检测试剂盒(诺唯赞生物),Runt相关转录因子2(RUNX2)抗体(CST,12556S),骨桥蛋白(OPN)抗体(CST),骨钙素(OCN)抗体(Sigma),快速RNA提取试剂盒(EZB),siRNA合成及qPCR引物合成(北京擎科生物科技股份有限公司)。SPF级6周龄雄性Balb/c小鼠及4周龄雄性C57BL/6J小鼠均购自南方医科大学动物实验中心,动物的使用方案获得南方医科大学动物医学伦理委员会审查批准(审批号:L2020005),并严格遵循美国NIH有关实验动物照料及使用的指南。
1.2 方法
1.2.1 GEO数据库以及Rstudio软件分析
使用GEOquery程序包获取GEO(GSE175624,GSE125331,GSE63754)数据集,以P<0.05且|log2FC|>1为标准,筛选ADSCs成脂或成骨分化后显著上调和下调的基因。使用Venn图进行差异基因交叉富集,使用pheatmap程序包根据标准化后的表达矩阵exp绘制热图。使用clusterProfiler程序包对表达矩阵进行GO富集分析。差异基因表达谱数据DEG使用limma包t检验和经验贝叶斯方法进行统计学分析;GO富集分析使用clusterProfiler包超几何检验或Fisher精确检验进行统计学分析。P<0.05为差异具有统计学意义。
1.2.2 mADSCs提取、培养、成骨诱导分化
选用4周龄雄性C57BL/6J小鼠(n=3),剪开小鼠下腹部至颈部皮肤,钝性分离皮肤以暴露腹股沟及腋窝处的白色脂肪组织。用剪刀镊子取下完整的大块组织并置于PBS中备用,漂洗血污及毛发后,将组织转移至干净的细胞皿中,用洁净的眼科剪将脂肪组织剪碎至糜状。转移组织至50 mL无菌离心管中,添加约10倍组织体积量的0.1% I型胶原酶,置于37 ℃恒温摇床内,85 r/min震荡消化。待脂肪组织完全消化后,将消化液转移至15 mL离心管中,离心弃上清,用含10%胎牛血清的L-DMEM重悬细胞,并以2×106细胞接种于25 cm2培养瓶中,在培养箱中孵育(37 ℃,5% CO2),全换液2~3 d/次,细胞贴壁融合至80%左右进行传代培养。将第3代细胞接种于12孔板,1×104 /孔,进行成骨诱导分化培养(L-DMEM,50 μg/mL维生素C,10 nmol/L地塞米松和10 mmol/L β-甘油磷酸),3 d/次换液。
1.2.3 碱性磷酸酶(ALP)染色及定量
室温下用4%多聚甲醛固定细胞10 min,按说明书配制碱性磷酸酶染液(碧云天),弃固定液,PBS清洗细胞3遍,24孔板每孔加入染液300 μL覆盖细胞,避光染色1 h,弃染液,去离子水洗2遍终止染色,晾干,扫描仪扫描,光学显微镜下拍照比较ALP染色的程度。向24孔板中加入全细胞裂解溶液(不含磷酸酶抑制剂),以50 μL/孔收集蛋白质,并使用BCA蛋白定量检测试剂盒和碱性磷酸酶检测试剂盒定量总蛋白和碱性磷酸酶含量。
1.2.4 茜素红(ARS)染色及定量
将细胞固定在4%多聚甲醛中,并在2% ARS溶液中染色。沉积的钙被染成红色。扫描仪扫描,光学显微镜下拍照比较着色钙结节数量。并用10% 的乙酸萃取和10% 的氨水中和对沉积的钙进行半定量。
1.2.5 免疫印迹实验
收集细胞,提取总蛋白,采用BCA蛋白定量试剂盒测定总蛋白质浓度,每孔上样30 μg,进行SDS-PAGE电泳,转膜,5% 脱脂牛奶液封闭l h。以β-actin为内参,加抗RUNX2(CST,1∶1000)、抗OPN(CST,1∶1500)、抗OCN(Sigma,1∶1000),4 ℃摇床孵育过夜,TBST洗3次,10 min/次,加入相应二抗,室温孵育l h,化学发光液染色后进行观察和图像分析。
1.2.6 RT-qPCR法
通过RT-qPCR检测相关基因的表达水平。PrimerBank上查找目的基因引物(
表1),通过RNA快速提取试剂盒提取细胞mRNA,进行RNA纯度和完整性检测,最后进行逆转录反应,选择
Gapdh作为内参。
1.2.7 小鼠不可自愈合颅骨缺损模型
用戊巴比妥钠麻醉6周龄的雄性Balb/c小鼠(
n=16),在右侧颅骨上钻取全层骨缺损(直径3 mm)。Blank组缺损未做任何处理,MeGC
[31](甲基丙烯酰化壳聚糖水凝胶)组只加入水凝胶进行填充,si-NC组和si-
Hmga2组分别加入MeGC包裹了转染si-NC或si-
Hmga2的ADSCs。术后6周,断颈处死小鼠并取材其颅骨组织。然后进行μCT扫描、组织学HE和Masson三色染色,染色步骤参照试剂盒说明书,蔡司显微镜拍照采集染色结果并进行骨缺损修复评估。
1.2.8 μCT扫描
将取材的颅骨固定在4%多聚甲醛溶液中48 h。然后对所有样本均使用Viva CT40(Scanco Medical AG, Bassersdorf)进行高分辨率μCT图像采集。采用55 kVp的电压值和15 μm的分辨率曝光50 min扫描。使用 ImageJ 测定所有样品的新骨表面积,并将其归一化为原始缺损表面积(直径 3 mm)。使用 SkyScan CT-Analyzer 程序提取骨体积密度、骨小梁数、骨小梁厚度和骨小梁分离度。
1.3 统计学分析
本文实验均至少独立重复3次,实验结果统计学分析以及图表制作使用GraphPad Prism 8.0软件,统计结果以均数±标准差表示。两组单因素组间统计采用Student's t-test分析;两组以上的单因素组间统计采用One-way ANOVA分析,并使用Tukey进行多重比较检验。当P<0.05时,认为差异具有统计学意义。
2 结果
2.1 Hmga2基因在成脂分化中显著高表达
根据“成脂-成骨”平衡理论(
图1A),通过Venn图对两个数据集中所有与成脂分化相关的基因进行交叉富集分析,共有504个基因在两个数据集中表达显著上调(
图1B)。进一步对DEG进行GO分析(
P<0.05,
图1C)。本研究对“fat cell differentiation”下的26个基因进行了可视化处理,并发现在成脂分化过程中稳定上调表达的转录调控因子
HMGA2(
图1D,E)。
2.2 抑制Hmga2促进mADSCs中ALP表达和矿化结节形成
RT-qPCR结果显示,与NC组相比,转染si-
Hmga2后,mADSCs细胞中
Hmga2的mRNA水平显著降低,其中si-
Hmga2-2的敲低效率最高(
P<0.05,
图2A)。后续的功能实验选用si-
Hmga2-2进行。分别在成骨诱导的第3天和第7天使用ALP染色及定量评估mADSCs的成骨能力,结果显示,抑制
Hmga2的表达后mADSCs的ALP染色显著增强,ALP活性相较于对照组分别提高至1.2倍和1.3倍(
P<0.001,
图2C、E)。成骨诱导14 d进行ARS染色,si-
Hmga2组相较于对照组红染矿化结节显著增多,定量结果显示钙结节数量为对照组的1.4倍,差异具有统计学意义(
P<0.001,
图2G)。
2.3 抑制Hmga2促进成骨分化标志物的表达,降低成脂分化标志物的表达
在成骨诱导的第3天,si-
Hmga2组Osterix(
Osx)基因的表达水平显著高于对照组,同时,成脂分化转录因子过氧化物酶体增殖物激活受体γ(
Pparγ)表达下调,成脂分化抑制因子delta样非典型Notch配体1(
Dlk1)表达上调,差异具有统计学意义(
P<0.01,
图3A-C);成骨诱导的第7天,si-
Hmga2组Runt相关转录因子2(Runx2)和骨桥蛋白(Opn)的基因(
P<0.01,
图3D-E)和蛋白表达水平(
P<0.05,
图3F)相较于对照组均显著增加。到了成骨诱导的第14天,si-
Hmga2组的骨钙素(OCN)蛋白表达也显著高于对照组(
P<0.05,
图3G)。上述实验结果与ALP及ARS染色结果相互印证,说明在mADSCs的分化进程中,
Hmga2的敲低会抑制mADSCs的成脂分化,推动mADSCs向成骨谱系分化(
图4)。
2.4 抑制Hmga2的表达促进体内骨缺损修复
术后6周取材,μCT扫描三维重建图像显示si-
Hmga2组ADSCs的骨缺损修复效果最佳,可见大量新骨形成,si-NC组和MeGC组在骨缺损孔中央亦有少量新骨形成,而Blank组只观察到新生骨从边缘处生长(
图5A)。通过骨计量学分析新生骨面积、骨体积分数(BV/TV%)、骨小梁数目(Tb.N)等,抑制
Hmga2后ADSCs促进骨缺损修复效果(
P<0.01,
图5B)。组织切片HE和Masson染色显示si-
Hmga2组有更丰富的胶原覆盖,且有较大面积新生骨形成(
图6)。
2.5 HMGA2在成骨分化过程中的蛋白互作网络分析及靶点预测
HMGA2归属于其中一个拥有54份蛋白质的聚类类别中(
图7A),对该聚类下的蛋白网络进行了进一步分析,并围绕HMGA2的可能互作关系绘制了卫星图(
图7B)。在聚类分析后,各节点的相互作用蛋白数目决定了它们在网络图中的位置分布,节点之间连线颜色的深浅表示相互作用的可能性大小,而HMGA2则位于网络图的第2圈。网络图的中心节点是由CDH1编码的钙粘蛋白-1,该蛋白参与调节细胞间粘附、迁移和增殖功能,与HMGA2作为转录因子,调控细胞周期和增殖功能密切相关。HMGA2预测的可能互作蛋白包括SMAD7、CDH1、CDH2、SNAI1、SMAD9、IGF2BP3、ALDH1A1(
表2)。在预测结果中,IGF2BP3被认为是HMGA2下游靶点的可能性最高,同时,GO及KEGG分析提示HMGA2可能通过激活TGF-β信号通路来调节成骨分化。
3 讨论
治疗骨不连和大面积骨缺损,并有效增强骨组织修复仍然是临床实践中的重大挑战
[32]。传统手术方法如假体植入和骨移植虽然具有一定临床疗效,但存在感染、疼痛、高手术成本以及需要额外手术等缺点
[4],限制了其应用。因此,迫切需要新的策略或改进现有治疗方法,以提供个性化诊治并加快患者康复。ADSCs作为骨组织工程的“宠儿”,探索其成骨分化的机制并实现在骨组织工程应用中的定向分化,能够有力地推动ADSCs在骨缺损修复中的应用
[13]。本研究从公开数据库GEO收集并筛选了hADSCs成脂分化前后的RNA-seq数据集进行分析,结果表明hADSCs成脂分化前后存在大量差异表达基因。“成脂-成骨”分化存在此消彼长的平衡,因此明确促进成骨分化,同时抑制成脂分化的关键节点分子具有重要的意义。众多研究指出
Hmga2在成骨分化进程中尤其是在MSCs的成骨分化中发挥着十分重要的作用,然而,目前尚未有研究指出
Hmga2与“成脂-成骨”平衡存在直接相关关系,同时,不同研究报道的
Hmga2在成骨分化中的作用存在分歧。
本研究明确了敲低Hmga2能够促进ADSCs成脂分化并抑制成骨分化,具有调控其“成脂-成骨”分化平衡的能力。在体外细胞和体内骨缺损模型中均证实了抑制Hmga2的表达能够促进mADSCs向成骨谱系分化,促进骨再生修复,明确了Hmga2的抑制促进ADSCs成骨分化这一结论。结合上述实验结果,我们认为Hmga2在ADSCs成骨分化命运决定中发挥着关键节点分子的作用,调控“成脂-成骨”平衡,这也是本研究的核心及创新点。本研究通过数据库及PPI蛋白互作网络分析,预测出HMGA2的下游潜在靶点IGF2BP3及TGF-β信号通路。遗憾的是,本研究通过生物信息学预测找到的靶点和通路未能进一步探索和验证,目前尚停留在理论预测阶段。在后续研究中,将会重点完成HMGA2对成骨分化影响的调控网络构建和机制的研究,以阐明HMGA2影响“成脂-成骨”分化平衡过程中的信号通路。
本研究的另一个创新点在于将抑制
Hmga2促进ADSCs成骨分化的机制应用于动物模型。在免疫缺陷小鼠的颅骨缺损不可自愈合模型中
[33],利用MeGC光敏水凝胶包裹抑制
Hmga2表达的mADSCs制备的组织工程系统展现出了卓越的骨再生修复效果,这一实验结果有力地支持了
Hmga2作为间充质干细胞分化命运调节因子的证据,更重要的是,这为解决大面积骨缺损、衰老或肥胖伴随骨质疏松等骨骼疾病提供了新的基因治疗靶点。
综上所述,HMGA2作为ADSCs成骨分化命运决定的关键节点分子,通过抑制Hmga2可以改变“成脂-成骨”平衡,抑制ADSCs的成脂分化,促进其成骨分化。在此过程中,成骨标记物RUNX2、OPN和OCN表达显著上调,矿化钙结节显著增加;此外,抑制Hmga2还能够加速小鼠体内骨缺损修复,在骨组织工程修复骨缺损的应用中展现出巨大的潜力。