MSC C1：肺腺癌骨转移微环境中连接基质可塑性、免疫调节与患者生存的潜在枢纽

《Chemical Biology & Drug Design》：MSC C1: A Potential Hub Linking Stromal Plasticity, Immune Regulation, and Patient Survival in Lung Adenocarcinoma

　　

肺腺癌（LUAD）骨转移直接导致不良预后。本研究旨在剖析伴骨转移LUAD的基质景观。对伴骨转移LUAD患者的单细胞RNA测序（scRNA-seq）数据进行无监督聚类，识别基质亚群。伪时间分析联合CytoTRACE干性评分重建间充质基质细胞（MSC）分化轨迹，同时CellChat刻画MSC亚群间的细胞–细胞通讯网络。对功能不同的MSC亚群特异性差异表达基因（DEGs）应用LASSO-Cox回归，构建预后风险评分（RiskScore）模型，并在独立批量RNA测序（RNA-seq）队列中就总生存（OS）、免疫细胞浸润和药物敏感性进行验证。一个独特的高干性MSC C1亚群定位于两条分化轨迹的分叉点，分别趋向基质重塑命运和分泌性、免疫调节表型。RiskScore模型将患者分为低风险和高风险组，两组总生存结局不同。高风险患者中观察到细胞毒性淋巴细胞浸润减少，且高风险亚组内药物敏感性更高。体外实验验证LDHA敲低损害A549细胞的增殖和迁移。本研究为风险适配治疗中靶向基质来源基因提供了依据。

论文主体解读

研究背景：肺腺癌（LUAD）是非小细胞肺癌的常见组织学亚型，也是癌症相关死亡的重要原因；骨是LUAD远处转移最常见部位之一，超过40%的晚期患者受累，骨转移直接导致不良预后，而针对骨转移龛的治疗仍有限。骨转移进展不仅依赖肿瘤细胞，还依赖肿瘤细胞与宿主微环境的动态串扰；骨髓生态位由基质、免疫、血管和造血成分构成，支持肿瘤定植、存活和耐药。间充质基质细胞（MSCs）是骨髓关键血管周围驻留细胞，参与维持骨骼稳态、调节免疫反应和细胞外基质（ECM）动态；在乳腺癌和前列腺癌等肿瘤中，肿瘤教育MSC可促进免疫抑制、血管生成、骨溶解和化疗耐药。然而，肺癌骨转移中MSC区室的功能异质性仍不清楚，主要受限于bulk转录组无法解析细胞多样性；人LUAD骨转移的scRNA-seq全面刻画有限，基质群体的转录状态、谱系轨迹和细胞间通讯网络认识不足，且基质来源分子模型能否作为预测预后和治疗反应的可靠生物标志物尚不明确。为此，研究人员整合LUAD骨转移scRNA-seq与多独立队列bulk RNA-seq及临床数据，刻画骨转移龛细胞组成、基质异质性和谱系轨迹，绘制细胞互作，鉴定基质特异基因构建预后模型，并评估免疫浸润和药物敏感性，体外验证LDHA。总体结论是：高干性MSC C1亚群位于两条分化轨迹分叉点并作为通讯枢纽；基于MSC C1标记的7基因RiskScore可分层总生存（OS）；高风险组细胞毒性淋巴细胞浸润减少且药物敏感性更高；LDHA敲低抑制A549细胞增殖和迁移。其意义在于为风险适配治疗中靶向基质来源基因提供依据，并提示将基质特征纳入LUAD骨转移风险评估。

关键技术方法（不超过250字）：研究人员使用GSE266330中10例肺癌骨转移患者scRNA-seq数据，以及TCGA-LUAD（490例原发肿瘤、58例癌旁）和GSE31210（226例LUAD）队列。关键方法包括DoubletFinder去除双细胞、decontX校正环境RNA、SCTransform归一化、PCA与Harmony去批次、UMAP降维及无监督聚类，依据CellMarker 2.0和既往研究注释；Monocle2伪时间分析联合CytoTRACE干性评分重建MSC轨迹；CellChat推断细胞通讯；limma筛选差异表达基因（DEGs），FindAllMarkers鉴定亚群标记，单因素Cox、最小绝对收缩和选择算子（LASSO）-Cox及赤池信息准则（AIC）逐步多因素Cox构建RiskScore，Kaplan–Meier（KM）/log-rank和timeROC评估；单样本基因集富集分析（ssGSEA）与TIMER分析免疫浸润；oncoPredict预测IC₅₀；体外采用qRT-PCR、CCK-8、划痕愈合、Transwell和Western blot，并以siRNA敲低LDHA。

主要研究结果：
3.1 Mesenchymal Stromal Cells Were Enriched and Robustly Identified in the Bone Metastatic Niche：通过严格质控、归一化、降维和无监督聚类，从GSE266330获得26,822个高质量细胞，分为10个簇，并注释为T细胞、单核细胞、MSCs、B细胞、内皮细胞、造血干/祖细胞（HSPCs）、浆细胞和浆细胞样树突状细胞（pDCs）；MSC因跨样本一致存在、丰度较高且稳健表达THY1和LEPR而被选进一步分析。
3.2 Transcriptional Heterogeneity Underlied Functional Diversification of MSCs in Bone Metastasis：通过对MSC再聚类，t-SNE降维识别C1、C2、C3三个亚群，各具独特标记基因；KEGG富集显示C1富集于核糖体、ECM-受体相互作用和HIF-1信号通路，C2富集于PI3K-Akt和TNF信号通路，C3富集于cGMP–PKG和cAMP信号通路，提示MSC异质性促成骨转移中的功能多样化。
3.3 MSC Heterogeneity Reflected Divergent Differentiation Trajectories Toward Matrix-Remodeling and Secretory Fates：通过Monocle2轨迹推断，MSC沿连续轨迹发生状态转换并在末端形成两个分支；CytoTRACE评分最高的State 2被定为轨迹根，多数C1细胞位于初始状态附近；一支富集细胞连接组装调控、骨化、细胞外结构组织、ECM组织和细胞黏附负调控，趋向基质重塑/成骨命运，另一支富集内质网蛋白定位和蛋白靶向，提示分泌能力增强。
3.4 MSC C1 Acted as a Central Communication Hub：通过CellChat分析，MSC与多种细胞广泛互作，C1与邻近细胞通讯最活跃并接收更多传入信号；C1主要通过TGFβ1/3、VEGFA、PDGFs和WNT5B及TGFBR1/2、VEGFR1/2、PDGFRs、FZD1/4向单核细胞和内皮细胞传递信号，关联血管生成、免疫抑制和ECM重塑；C1接收B细胞、单核细胞和T细胞经MIF-ACKR3、MDK-SDC1、PTN-NCL和IGF1-IGF1R轴的信号，其中MIF-ACKR3富集，提示C1为双向通讯枢纽。
3.5 Construction and Validation of a Clinical Prognostic Model：通过将MSC C1标记与TCGA转录组整合，2508个DEGs中102个与C1标记重叠，单因素Cox筛选34个与OS显著相关基因，LASSO-Cox在λ=0.0172降至16个，AIC逐步选择最终建立7基因模型：RiskScore=0.297*PPFIBP1?0.343*DOCK4+0.474*LDHA+0.078*NT5E+0.152*C4orf48?0.205*NUPR1+0.129*ID1；TCGA中高风险组OS显著缩短，1–5年OS的受试者工作特征（ROC）曲线下面积（AUC）为0.65–0.75，GSE31210验证队列中模型仍可有效分层且AUC较高。
3.6 Risk Stratification Was Correlated With Immune Activation and Differential Chemotherapeutic Response：通过ssGSEA和TIMER分析，低风险组活化B细胞、CD8⁺ T细胞、效应记忆CD8⁺ T细胞、巨噬细胞和未成熟B细胞浸润显著更高；通过oncoPredict预测，13种药物与RiskScore显著相关，ERK_6604_1714与RiskScore显著负相关，提示高风险组对该化合物更敏感。
3.7 LDHA Promoted Tumor Cell Proliferation and Migration in LUAD：通过qRT-PCR，A549细胞中LDHA、NT5E、C4orf48、NUPR1和ID1 mRNA表达显著高于BEAS-2B，LDHA最高；si-LDHA#1敲低LDHA后，CCK-8显示A549增殖受抑，Transwell和划痕愈合显示迁移受损；Western blot显示上述基因蛋白在A549中显著更高，提示LDHA促进LUAD细胞增殖和迁移。

讨论与结论：讨论部分指出，MSC C1具干性、广泛细胞间通讯和临床相关性，肿瘤教育MSC是转移龛的主动调节者而非仅结构成分；C1近轨迹根、高CytoTRACE评分，代表早期低分化状态，并与免疫和血管区室双向通讯，TGFβ、VEGFA和WNT信号可能参与血管生成和免疫抑制，MIF-ACKR3等炎症信号提示其整合微环境信号。基于C1标记的模型在TCGA-LUAD和GSE31210中稳健，高风险与活化B细胞、CD8⁺ T细胞和效应记忆亚群浸润减少相关，提示免疫“冷”表型，支持基质驱动免疫抑制与不良结局相关。模型基因涉及代谢激活、免疫调节、黏附运动及应激适应；LDHA在RiskScore中贡献最强并经A549验证，PPFIBP1和NT5E（CD73）高表达与风险增加相关，DOCK4和NUPR1高表达与较低RiskScore和较好预后相关，ID1连接上皮–间质转化（EMT）和干性。ERK_6604_1714为研究性ERK通路抑制剂，与RiskScore负相关，提示高风险组可能对MAPK信号靶向抑制更敏感，并为KRAS突变常见背景下生物标志物富集临床试验提供理由。局限性包括scRNA-seq仅10例患者、训练和验证队列均为未选择原发LUAD且缺乏骨转移患者独立验证队列、ERK_6604_1714仅基于计算分析、C1促转移因果作用仍多为推断，需更大样本、额外骨转移数据集、体外体内实验及共培养或条件敲除研究验证。

研究结论翻译：总之，本研究揭示MSC区室内的转录异质性促成LUAD骨转移龛的功能多样化。C1亚群成为基质–免疫串扰的关键介导者及临床相关生物标志物的来源。通过整合单细胞生物学与临床结局，这些发现凸显将基质特征纳入伴骨转移LUAD患者风险评估和治疗决策的潜在价值。