整合单细胞与Bulk转录组学揭示脂质代谢相关lncRNA特征及ELFN1-AS1/LYPLA1轴驱动骨肉瘤进展

《Frontiers in Immunology》：Integrative single-cell and bulk transcriptomics uncovers a lipid metabolism-related lncRNA signature alongside the ELFN1-AS1/LYPLA1 axis driving osteosarcoma progression

　　

背景：骨肉瘤是青少年中最常见的原发性恶性骨肿瘤，具有显著的瘤内异质性、早期肺转移以及转移性/化疗耐药患者生存改善停滞的特点。脂质代谢改变是恶性肿瘤的核心适应性标志，而脂质代谢相关长链非编码RNA（lipid metabolism-related long non-coding RNAs，LMRLs）在骨肉瘤中的预后意义和分子调控功能仍知之甚少。

方法：整合了来自GSE152048的单细胞RNA测序（single-cell RNA-seq）数据以及来自TARGET-OS的bulk转录组和临床数据。分别使用AUCell、Monocle和CellChat评估脂质代谢活性、恶性细胞状态转变和细胞间通讯。通过整合单细胞差异表达分析与加权基因共表达网络分析（weighted gene co-expression network analysis，WGCNA）鉴定脂质代谢相关基因。使用Pearson相关、Cox回归和LASSO分析构建预后lncRNA特征。进一步评估免疫浸润、肿瘤突变负荷、免疫检查点表达和药物敏感性。基于scTenifoldKnk的虚拟敲除筛选用于预测ELFN1-AS1的下游靶基因，随后对ELFN1-AS1/LYPLA1调控轴进行全面的体外和体内功能验证。

结果：单细胞分析揭示了骨肉瘤细胞群体中显著的脂质代谢异质性，伴随恶性细胞状态转变的动态重塑，以及高脂质代谢骨肉瘤细胞与基质或血管区室之间增强的通讯。通过整合单细胞和bulk转录组证据，鉴定了22个脂质代谢相关候选基因，并由此开发了包含AL133410.1、AL596247.1、ELFN1-AS1、IL10RB-DT和NECTIN3-AS1的五-LMRL特征。该特征有效将患者分为预后不同的风险组，并作为总生存期的独立预测因子。探索性计算分析表明，风险亚组表现出不同的免疫景观、突变谱和预测的药物反应性。虚拟敲除筛选将LYPLA1确定为ELFN1-AS1的关键下游效应分子。在特征lncRNA中，ELFN1-AS1在骨肉瘤组织和细胞中显著上调，并与不良预后相关。机制上，ELFN1-AS1通过触发LYPLA1依赖性脂质积累相关表型驱动骨肉瘤进展。使用奥利司他（orlistat，选择性FASN抑制剂）药理学阻断从头脂肪生成，可减弱ELFN1-AS1驱动的致癌表型。

结论：本研究建立了基于LMRL的骨肉瘤预后分层方法，反映了免疫、基因组和治疗异质性。结合虚拟敲除筛选和实验验证，研究人员将致癌性ELFN1-AS1/LYPLA1调控级联表征为骨肉瘤精准管理的潜在生物标志物和治疗靶点。

论文解读：整合单细胞与bulk转录组学揭示脂质代谢相关lncRNA特征及ELFN1-AS1/LYPLA1轴驱动骨肉瘤进展

骨肉瘤是青少年最常见的原发性恶性骨肿瘤，具有显著瘤内异质性、早期肺转移倾向，转移性/化疗耐药患者的生存改善长期停滞。脂质代谢重编程被认为是恶性肿瘤的核心适应性标志，但脂质代谢相关长链非编码RNA（LMRLs）在骨肉瘤中的预后价值与分子调控机制尚不明确。为此，研究人员整合单细胞与bulk转录组学，系统描绘骨肉瘤脂质代谢异质性，构建预后模型并深入解析关键lncRNA的功能机制。该研究发表于《Frontiers in Immunology》。

研究人员从GEO数据库获取GSE152048单细胞RNA测序数据（含11例骨肉瘤样本：7例原发、2例转移、2例复发），从TARGET数据库获取84例具有完整临床信息的骨肉瘤患者转录组与临床数据。采用Seurat和Harmony处理单细胞数据，AUCell量化脂质代谢活性，Monocle构建拟时轨迹，CellChat分析细胞通讯；通过单细胞差异表达基因与加权基因共表达网络分析（WGCNA）交集鉴定脂质代谢相关基因（LMRGs）；利用Pearson相关、LASSO和多因素Cox回归构建预后LMRLs特征；基于scTenifoldKnk进行虚拟敲除预测ELFN1-AS1下游靶基因。随后使用32对来自山东大学齐鲁医院的骨肉瘤及配对癌旁组织样本进行RT-qPCR验证，并通过脂质含量测定、油红O染色、CCK-8、克隆形成、EdU、Transwell及皮下异种移植模型等功能实验验证ELFN1-AS1/LYPLA1轴。

研究结果如下：

**单细胞图谱揭示骨肉瘤脂质代谢异质性与细胞间通讯**。对GSE152048中101,747个高质量细胞进行无监督聚类，鉴定出11种主要细胞类型。AUCell分析显示脂质代谢活性在不同细胞类型中呈异质性分布，高代谢活性主要富集于髓系和破骨样群体。骨肉瘤细胞再聚类得到7个亚群，高、低脂质代谢组间存在864个差异表达基因（DEGs），KEGG富集显示脂肪酸代谢、碳代谢等通路显著富集。CellChat分析发现骨肉瘤细胞与基质/血管成分存在广泛双向通讯。

**拟时轨迹与高/低脂质代谢细胞的通讯差异**。Monocle拟时分析显示脂质代谢活性沿恶性细胞状态转变发生动态重塑，软骨母细胞样及部分成骨细胞样亚群在拟时后期脂质代谢评分升高。高、低脂质代谢组均存在复杂通讯网络，但高脂质代谢组中骨肉瘤细胞与微环境细胞的相互作用更显著，主要涉及胶原-整合素、趋化因子、VEGF、MIF、FN1和SPP1等信号。

**WGCNA与单细胞整合鉴定骨肉瘤LMRGs**。在TARGET-OS队列中构建无尺度共表达网络，绿色模块与脂质代谢评分正相关最强（Cor=0.41，P=9e-05）。将该模块511个基因与单细胞DEGs取交集，获得22个LMRGs，富集于鞘脂代谢、糖鞘脂生物合成等通路。

**训练队列中LMRLs预后特征的构建**。以22个基因为枢纽，通过Pearson共表达分析获得958个LMRLs。单因素Cox和LASSO筛选后，多因素Cox建立包含AL133410.1、AL596247.1、ELFN1-AS1、IL10RB-DT和NECTIN3-AS1的五-LMRL特征。风险评分=AL133410.1×1.559 + AL596247.1×0.478 + ELFN1-AS1×0.426 + IL10RB-DT×(-1.237) + NECTIN3-AS1×(-0.405)。高风险组总生存显著较差，1/3/5年AUC分别为0.848/0.877/0.895，风险评分是独立预后因素。

**LMRLs风险模型的验证与列线图转化**。内部测试队列中，高风险组预后仍显著较差（P=0.0057），AUC为0.815/0.705/0.710；合并队列AUC为0.833/0.804/0.809。基于风险评分和临床变量构建的列线图在训练队列中预测1/3/5年总生存的AUC达0.940/0.748/0.774，校准曲线显示良好一致性。

**LMRLs风险亚组的功能、免疫、基因组和治疗异质性**。高、低风险组间有679个DEGs，富集于ECM受体互作、PI3K-Akt、MAPK等通路。低风险组免疫细胞浸润评分（自然杀伤T细胞、活化CD8⁺T细胞、巨噬细胞等）及免疫功能评分较高，而高风险组肿瘤纯度更高，且CD276、HAVCR2、LAG3等免疫抑制分子上调。高风险组肿瘤突变负荷（TMB）更高，且对伏立诺他、阿西替尼、索拉非尼等预测更敏感，低风险组对5-氟尿嘧啶等更敏感。

**ELFN1-AS1被确定为关键LMRL**。单细胞分析显示ELFN1-AS1主要在恶性骨肉瘤细胞中表达，且在成骨样OS-2亚群富集。scTenifoldKnk虚拟敲除预测下游基因中LYPLA1变化最显著。KEGG富集显示与氧化磷酸化、甘油磷脂代谢等有关。

**ELFN1-AS1高表达与骨肉瘤不良预后相关**。RT-qPCR证实ELFN1-AS1在32对骨肉瘤组织中显著上调，且转移组高于非转移组，与远处转移和晚期Enneking分期相关。高表达组无进展生存期显著缩短。

**ELFN1-AS1在体外和体内作为骨肉瘤进展的致癌驱动因子**。敲低ELFN1-AS1显著抑制U2OS和143B细胞增殖、克隆形成、迁移和侵袭，并延缓皮下移植瘤生长，减少肺转移结节数。

**ELFN1-AS1通过激活LYPLA1介导的脂质积累相关表型驱动骨肉瘤进展**。ELFN1-AS1敲低降低细胞内甘油三酯、总胆固醇和脂滴含量，并抑制FASN、ACLY等脂肪生成基因。ELFN1-AS1主要定位于细胞质，通过稳定LYPLA1 mRNA上调其表达。回补实验证实LYPLA1敲低或奥利司他（FASN抑制剂）可逆转ELFN1-AS1过表达引起的脂质积累和恶性表型；而LYPLA1过表达可挽救ELFN1-AS1敲低导致的表型变化。研究提出ELFN1-AS1通过稳定LYPLA1 mRNA触发脂质积累相关表型，从而促进骨肉瘤增殖和转移。

讨论部分指出，该研究揭示了脂质代谢状态与恶性细胞状态转变及微环境通讯的关联，支持LMRLs作为预后分层工具。低风险组可能呈现免疫活化微环境，高风险组则显示免疫抑制特征。研究局限性包括缺乏独立外部队列验证、ELFN1-AS1调控LYPLA1 mRNA的中间分子尚未明确、缺乏脂质组学和同位素示踪、皮下异种移植模型不能完全模拟骨微环境、免疫相关分析仅为计算推断等。

结论：该研究整合单细胞转录组、bulk RNA-seq、预后建模、虚拟敲除筛选和实验验证，系统表征了骨肉瘤中的LMRLs。研究揭示了恶性细胞状态间的脂质代谢异质性，建立了五-LMRL特征用于风险分层，并鉴定ELFN1-AS1为通过LYPLA1依赖性脂质积累驱动骨肉瘤进展的关键致癌lncRNA。ELFN1-AS1/LYPLA1轴有望成为预后生物标志物和治疗靶点，选择性FASN抑制剂如奥利司他靶向抑制从头脂肪生成可能为高危骨肉瘤提供代谢干预策略。