《The FASEB Journal》：Integrative Transcriptomic Profiling Delineates LTBP4-Driven Fibroblast Reprogramming Across the Ossification Continuum of the Ligamentum Flavum
编辑推荐：
为了描绘黄韧带骨化（OLF）在正常、未成熟和成熟阶段的动态转录组图谱，鉴定核心基因，并筛选阶段特异性候选药物，研究人员收集了3例多节段OLF患者的正常、未成熟和成熟OLF组织（共9个样本）进行RNA-seq。研究人员进行了差异表达、功能富集和动态趋势分析，并使
为了描绘黄韧带骨化（OLF）在正常、未成熟和成熟阶段的动态转录组图谱，鉴定核心基因，并筛选阶段特异性候选药物，研究人员收集了3例多节段OLF患者的正常、未成熟和成熟OLF组织（共9个样本）进行RNA-seq。研究人员进行了差异表达、功能富集和动态趋势分析，并使用外部GEO数据集验证关键基因，同时整合单细胞和空间转录组学以确定细胞起源和空间定位，最后开展了基于AI的药物预测和分子对接。在正常与未成熟组织之间鉴定出1026个差异表达基因（DEGs）；在未成熟与成熟组织之间鉴定出713个DEGs；在正常与成熟组织之间鉴定出1285个DEGs。上调基因富集于RNA加工和核糖体生物发生，而下调基因富集于细胞外基质（ECM）-受体相互作用。LTBP4在骨化过程中呈进行性下调（AUC=0.88）。公共多组学数据集（GSE255942）的单细胞转录组学和空间转录组学分析显示，LTBP4表达于成纤维细胞，且在骨化组织中的表达低于正常组织。在成纤维细胞中虚拟敲除Ltbp4可重现OLF转录组特征。基于AI的药物预测鉴定出阶段特异性候选药物，其中BRD-K59831625与LTBP4的结合能最低（?8.9 kcal/mol）。持续的LTBP4下调驱动OLF进展，成纤维细胞是其核心效应细胞。研究人员还鉴定出阶段特异性候选药物。这些发现为OLF的分子诊断和精准治疗提供了新靶点。
**LTBP4驱动的黄韧带骨化成纤维细胞重编程：整合转录组学研究的突破**
黄韧带骨化（ossification of the ligamentum flavum，OLF）是导致胸椎管狭窄和脊髓压迫的重要原因，在东亚人群中患病率显著高于西方国家。该病常累及多节段，并可能合并后纵韧带骨化，临床异质性强、误诊率高。目前后路椎板切除减压是主要手术方式，但术后脊髓损伤、脑脊液漏、症状复发或邻近节段骨化进展等并发症发生率较高，关键难点在于术前难以准确评估不同节段骨化病变的成熟度与进展倾向。既往研究多聚焦于少数候选基因或经典信号通路（如TGF-β/BMP、Wnt/β-catenin等），缺乏对OLF从起始到成熟全过程转录组动态变化的系统描绘。为此，研究人员依托影像-病理分类，利用同一患者术中获取的正常、未成熟和成熟黄韧带组织，结合多组学验证与AI药物预测，系统绘制OLF的分子进化图谱，并筛选核心基因和阶段特异性干预药物。该成果发表于《The FASEB Journal》。
研究人员从宁夏医科大学总医院纳入3例多节段胸椎OLF患者（2男1女），术中依据CT值和MRI信号特征（影像-病理分类）精确区分并采集正常（N）、未成熟（I）和成熟（M）黄韧带组织，共9个样本行RNA测序（RNA-seq）。通过差异表达基因筛选、基因本体（GO）和京都基因与基因组百科全书（KEGG）富集分析、动态趋势分析，并结合外部数据集GSE181716验证，鉴定出核心基因LTBP4。随后整合公共单细胞转录组数据集GSE126060和空间转录组数据集GSE255942，解析LTBP4的细胞来源与空间分布；利用scTenifoldKnk进行虚拟敲除；基于CMap数据库进行药物重定位，通过分子对接评估候选药物与LTBP4的结合亲和力。
**3.1 患者人口统计学和样本概览**：3例患者均经术前影像和术中病理证实为胸椎OLF，依据骨化过程分为N、I、M三组，每组3个样本，共9份标本。典型影像显示，未成熟骨化灶密度不均、边界模糊，成熟骨化灶密度均匀、边界清楚；手术标本中不同阶段骨化灶肉眼可辨。
**3.2 差异表达基因**：在N vs I、I vs M、N vs M三组比较中分别鉴定出1026、713和1285个差异表达基因（DEGs），且不同阶段DEG的数目和上下调模式存在明显差异，提示疾病进展过程中基因表达发生大规模动态改变。
**3.3 功能富集分析**：GO和KEGG富集分析显示，上调DEGs主要富集于RNA加工、核糖体生物发生、DNA复制和细胞周期等细胞增殖与蛋白合成相关通路；下调DEGs主要富集于细胞外基质（ECM）组织、黏着斑和ECM-受体相互作用等微环境重塑通路。RNA代谢相关通路的显著富集是新的发现，提示OLF进展中可能存在广泛的转录后调控重编程。
**3.4 关键基因LTBP4的筛选和验证**：通过韦恩图和热图筛选出16个持续上调和34个持续下调的基因。其中LTBP4在正常、未成熟和成熟三个阶段中呈最显著的进行性下调。GeneMANIA调控网络显示LTBP4与TGF-β信号通路及ECM相关基因（如FBN1、FBN2、LTBP1）密切互作。外部数据集GSE181716验证LTBP4在OLF组织中显著低表达，ROC曲线显示其区分骨化与正常组织的曲线下面积（AUC）为0.88。
**3.5 单细胞转录组特征**：对小鼠异位骨化模型（GSE126060）的单细胞数据分析显示，损伤后第7天成纤维细胞与其他细胞类型（如巨噬细胞、内皮细胞）的通信强度达峰值，第21天显著下降；Ltbp4在成纤维细胞中清晰表达，并随损伤时间进行性降低，这与人类OLF标本的观察高度一致。
**3.6 空间转录组定位**：对多组学数据集GSE255942的空间转录组分析显示，正常组织中成纤维细胞散在分布，而骨化组织中成纤维细胞聚集于骨化灶周围；LTBP4表达与成纤维细胞定位高度重叠，且骨化组表达强度低于正常组，进一步证实成纤维细胞是LTBP4的主要来源。
**3.7 Ltbp4虚拟敲除后的差异表达基因**：利用scTenifoldKnk在正常组单细胞数据中虚拟敲除Ltbp4后，共有46个基因显著上调，包括ECM重塑相关基因（Col3a1、Col5a1、Col6a1、Dcn、Nid1）、免疫炎症相关基因（C1qa/b/c、Cd68、Adgre1、Csf1r）和溶酶体蛋白酶基因（Ctss、Ctsc、Lgmn），其中Timp2上调最显著。该转录特征与人类OLF组织的差异表达谱相似，提示LTBP4是ECM产生和免疫激活的负调控因子。
**3.8 药物预测和分子对接**：基于三组DEGs进行CMap药物预测，分别选取各阶段评分最高的前3种候选药物，并以LTBP4为靶蛋白进行分子对接。早期（N vs I）的BRD-A41450521和BRD-K93023739结合能均为?7.5 kcal/mol；进展期（I vs M）的BRD-K14880289结合能最低（?7.7 kcal/mol）；终末期（N vs M）的BRD-K59831625结合能最低（?8.9 kcal/mol），显示最强结合亲和力。
讨论部分指出，LTBP4的持续下调可能解除其对TGF-β信号的负调控，导致局部TGF-β活性异常升高，驱动成纤维细胞向成骨样细胞转分化。与腰椎管狭窄中LTBP4表达上调的研究不一致，可能反映疾病特异性和阶段特异性。虚拟敲除结果支持LTBP4作为ECM和免疫基因主负调控因子的作用。BRD-K59831625可能通过抑制ALK2发挥抗骨化效应。研究局限性包括样本量较小、外部队列来自急性损伤模型等。研究结论可翻译为：本研究首次系统描绘了黄韧带骨化从正常到成熟阶段的动态转录组图谱，揭示了以LTBP4持续下调和RNA代谢通路激活为特征的分子机制，鉴定成纤维细胞为骨化过程的核心效应细胞，并通过虚拟敲除实验证明Ltbp4缺失可直接诱导ECM和免疫基因激活，模拟人类OLF转录组。此外，基于不同病理阶段预测了具有阶段特异性干预潜力的候选药物，为OLF的分子诊断、进展预测和精准治疗提供了新理论依据。未来需在更大队列中验证核心基因表达变化，并通过体内外功能实验阐明LTBP4的具体作用，为开发靶向治疗策略奠定基础。
