整合批量与单细胞转录组学构建肺腺癌胞葬相关预后模型并验证关键基因LDHA

Integrated bulk and single-cell transcriptomics to develop an efferocytosis-related prognostic model for lung adenocarcinoma and validate the key gene LDHA.

作者信息Weipeng Zhang, Quanyue Chen, Yue Li, Li Zhang, Mingxi Wang, Boxu Sun, Shuo Yang, Fan Yang, Haiying Dong
PMID41345972
发布时间2025-12-04
DOI10.1186/s40001-025-03622-z

实验完整度

高

整合多组学数据(TCGA、GEO、scRNA-seq、stRNA-seq)构建并验证预后模型,并结合体外细胞实验(THP-1来源巨噬细胞极化、共培养功能实验)进行功能验证。

主要模型

THP-1细胞来源的巨噬细胞(M0/M2) A549细胞系 H1975细胞系 TCGA-LUAD和GEO(GSE31210, GSE26939, GSE3141)LUAD患者队列

重点核对

THP-1细胞向M2巨噬细胞极化的诱导条件:100 ng/mL PMA处理24小时,20 ng/mL IL-4和IL-13处理48小时 LDHA敲低慢病毒载体及对照载体,嘌呤霉素(1.5 μg/mL)筛选72小时 共培养体系:0.4 μm孔径Transwell小室,M2巨噬细胞与A549或H1975细胞共培养48小时 细胞增殖(CCK-8)、迁移和侵袭(Transwell)实验的检测时间点,迁移24小时、侵袭48小时 细胞实验独立重复≥3次,统计分析采用Student's t检验

摘要

背景:胞葬作用,即吞噬细胞对凋亡细胞的清除,发挥抗炎和促肿瘤的双重作用,而现有研究尚未充分阐明其在肺腺癌(LUAD)中的潜在机制。本研究旨在鉴定LUAD中与胞葬作用相关的预后基因及潜在治疗靶点。方法:从UCSC Xena和GEO数据库获取LUAD相关基因表达数据集,从GeneCards数据库获取胞葬相关基因(ERGs)。利用无监督共识聚类分析将LUAD患者分为两个不同的亚组,随后进行差异基因表达分析以鉴定差异表达基因(DEGs)。随后,通过LASSO-Cox回归构建了一个四基因预后模型,并通过Kaplan-Meier(KM)分析和受试者工作特征(ROC)曲线评估进一步验证。单细胞RNA测序(scRNA-seq)数据分析用于研究与胞葬作用相关的细胞亚群,并确定所涉及的关键基因。进行体外实验以确定所选关键基因在LUAD进展中的功能重要性。结果:成功为LUAD构建了一个包含四个基因的预后特征,并计算风险评分以将患者分为低风险和高风险组。值得注意的是,两个风险组呈现出不同的生存结局和免疫细胞浸润(ICI)。scRNA-seq数据分析揭示巨噬细胞是LUAD中主要的胞葬相关细胞类型,并确定LDHA为一个关键调控基因。拟时序轨迹分析显示,随着巨噬细胞分化为肿瘤相关巨噬细胞(TAMs),LDHA表达和胞葬相关活性逐渐增加。此外,体外实验数据显示,巨噬细胞中LDHA沉默促进了向M1表型的极化,同时抑制了M2极化。此外,将LUAD细胞与LDHA沉默的M2巨噬细胞共培养,显著降低了癌细胞的增殖、迁移和侵袭。细胞通讯分析提示,LDHA通过调节巨噬细胞与癌症相关成纤维细胞(CAFs)之间的SPP1通路,激活ECM受体及下游效应分子,从而促进M2型TAMs。结论:我们为LUAD构建了一个胞葬相关预后基因模型。此外,LDHA可能靶向M2巨噬细胞,突显了其作为LUAD治疗靶点的潜力。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
胞葬相关基因在肺腺癌中的预后意义及其潜在机制,特别是关键基因LDHA在巨噬细胞极化和肿瘤进展中的作用。
核心机制
LDHA可能通过调控巨噬细胞与癌症相关成纤维细胞(CAFs)之间的SPP1信号通路(包括SPP1-CD44和SPP1-integrin),激活ECM受体及下游效应分子,从而促进M2型肿瘤相关巨噬细胞(TAMs)极化。
主要证据
在肺腺癌中构建并验证了包含LDHA、PLK1、KRT8和FCRL1的预后模型;单细胞RNA测序和拟时序分析显示LDHA表达和胞葬活性随巨噬细胞向TAMs分化而增加;体外实验表明LDHA敲低促进M1极化并抑制M2极化,且LDHA敲低的M2巨噬细胞共培养抑制LUAD细胞增殖、迁移和侵袭。
研究意义
该研究为肺腺癌的预后评估和治疗干预提供了新的视角,并提示LDHA作为肺腺癌潜在治疗靶点的可能性。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

胞葬相关亚组鉴定

基于胞葬相关基因表达特征,对LUAD患者进行分子分型并评估其预后差异。

对来自TCGA-LUAD和GSE31210的742例LUAD样本,基于13个预后相关ERGs进行无监督共识聚类,分为ERGclusterA和ERGclusterB两个亚组,并比较其总生存期。

2

预后模型构建与验证

利用差异表达基因构建并验证LUAD患者的预后风险评分模型。

基于ERG亚组间的差异表达基因,通过LASSO-Cox回归和多因素Cox回归构建包含LDHA、PLK1、KRT8、FCRL1的四基因风险评分模型,并在训练集(TCGA-LUAD和GSE31210)和独立验证集(GSE26939、GSE3141)中通过KM曲线和ROC曲线评估其预测效能。

3

肿瘤微环境与治疗反应分析

探索不同风险组间免疫细胞浸润、肿瘤突变负荷及药物敏感性的差异。

利用多种算法(ESTIMATE、CIBERSORT、IOBR)评估风险组间的免疫细胞浸润;通过TIDE和IPS预测免疫治疗反应;使用maftools分析突变谱;利用oncoPredict预测化疗药物敏感性。

4

单细胞RNA测序分析

在单细胞水平鉴定LUAD中与胞葬作用相关的关键细胞类型和关键基因。

对scRNA-seq数据GSE200972进行质量控制、Harmony整合、Seurat聚类和SingleR注释;通过AddModuleScore计算胞葬评分;通过monocle进行拟时序分析揭示巨噬细胞分化轨迹;通过CellChat分析细胞通讯。

5

LDHA对巨噬细胞极化的功能验证

验证LDHA对巨噬细胞M1/M2极化的调控作用。

通过慢病毒介导的shRNA在THP-1细胞中敲低LDHA,经PMA诱导分化为M0,再用IL-4和IL-13诱导为M2,通过qRT-PCR和Western blot验证敲低效率,并检测M1/M2标志物表达。

6

LDHA对LUAD细胞恶性行为的影响

评估LDHA敲低的M2巨噬细胞对LUAD细胞增殖、迁移和侵袭能力的影响。

将LDHA敲低或对照的M2巨噬细胞与A549或H1975细胞共培养48小时,或使用其条件培养基(CM),通过CCK-8、集落形成、划痕愈合和Transwell迁移/侵袭实验检测LUAD细胞行为变化。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
RPMI 1640培养基----
青霉素-链霉素----
胎牛血清----
β-巯基乙醇----
佛波醇12-十四酸酯13-乙酸酯Beyotime--
重组人白细胞介素-4Novoprotein--
重组人白细胞介素-13Novoprotein--
Transwell小室(0.4 μm孔径)Corning--
LDHA靶向慢病毒载体GenePharma--
嘌呤霉素GenePharma--
TRNzol试剂Tiangen--
PrimeScript逆转录预混液TAKARA--
TB Green预混Ex Taq IITAKARA--
RIPA裂解液Boster--
蛋白酶抑制剂----
PVDF膜Boster--
Multi-rAb HRP-山羊抗兔二抗Proteintech--
增强化学发光试剂Boster--
细胞计数试剂盒-8Biosharp--
结晶紫----
Transwell小室(8 μm孔径)Corning--
Matrigel基质胶Corning--
抗LDHA抗体Wanleibio--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
数据获取与预后模型构建
TCGA-LUAD和GSE31210数据集的下载来源,GeneCards数据库中ERGs的筛选标准(relevance score > 0.4),患者分组标准(ERGs聚类),LASSO回归参数,多因素Cox回归纳入基因,风险评分公式中基因系数
阅读提示:Methods: Data sources, Consensus clustering based on ERGs, Development of prognostic model
细胞培养与分化
THP-1、A549、H1975细胞系的来源,培养基组成,PMA诱导浓度和时间,IL-4和IL-13浓度和时间,M2极化验证(标志物)
阅读提示:Methods: Cell culture, Differentiation of macrophages
慢病毒转染与敲低
LDHA shRNA序列信息,慢病毒载体构建,感染条件,嘌呤霉素筛选浓度和时间,敲低效率验证方法(qRT-PCR、WB)
阅读提示:Methods: Cell transfection; 表S2
共培养与条件培养基制备
Transwell小室孔径,共培养时间,CM制备步骤(血清饥饿时间、离心速度和时间、过滤膜孔径),CM与新鲜培养基混合比例
阅读提示:Methods: Co-culture system and preparation of conditioned medium
功能实验
CCK-8实验种植密度和检测时间,集落形成实验种植密度和培养时间,划痕愈合实验细胞密度,迁移实验细胞密度和孵育时间,侵袭实验Matrigel稀释比例和孵育时间
阅读提示:Methods: Cell proliferation assay, Colony formation assay, Wound healing assay, Cell migration and matrigel invasion assays
统计学分析
独立重复次数,统计检验方法(Student's t检验、Wilcoxon秩和检验等),显著性阈值
阅读提示:Methods: Statistical analysis