肝细胞癌(HCC)是全球第6大常见恶性肿瘤,同时也是导致癌症相关死亡的第3大原因
[1]。超过90%的HCC发生在慢性肝病背景下
[2],HBV和HCV感染仍是全球HCC的主要诱因
[3, 4];此外,代谢功能障碍相关脂肪性肝病(MASLD),已被公认为肝细胞癌发生的重要驱动因素
[5]。肝细胞癌早期诊断率不足30%,70%~80%的患者初诊即为中晚期
[6, 7]。
肿瘤相关巨噬细胞(TAMs)是肝细胞癌(HCC)肿瘤微环境中最丰富的免疫细胞群之一
[8, 9]。TAMs以M2样表型为主,该表型与肝细胞癌进展和耐药密切相关。TAM通过促进免疫抑制、肿瘤血管新生和肿瘤干细胞形成,进而促进肝细胞癌进展以及耐药
[9-14]。近年来研究表明,膜联蛋白(ANXAs)家族成员在调控肿瘤相关巨噬细胞表型以及免疫应答中发挥重要作用。如在肝细胞癌中,巨噬细胞ANXA1表达水平上调,其表达水平与PD-L1呈正相关;下调ANXA1的表达可增加M1/M2巨噬细胞比例和促进T细胞活化,从而抑制肝癌细胞增殖和迁移
[15]。在胶质母细胞瘤中,ANXA2和ANXA4主要在M2型巨噬细胞中表达,参与免疫抑制微环境形成
[16]。在非酒精性脂肪性肝炎模型中,ANXA5诱导肝脏巨噬细胞由M1向M2表型极化
[17]。但是,ANXAs家族对肝细胞癌预后以及肿瘤相关巨噬细胞影响尚不清楚。
为探索ANXAs家族对肝细胞癌预后以及肿瘤相关巨噬细胞的影响,对HCC多组学数据进行生物信息学分析,并且通过体外实验使用肝癌细胞培养上清与巨噬细胞THP-1共培养,探究对ANXA2和ANXA5对巨噬细胞极化的影响,为肝细胞癌肿瘤免疫以及靶向治疗,提供新的思路和分子靶点。
1 材料和方法
1.1 生物信息学数据集来源
空间转录组数据集:GSE203612;TCGA-LIHC数据集:TCGA数据库中TCGA-LIHC项目STAR流程的RNAseq数据并提取TPM格式的数据以及临床数据。
1.2 生物信息学数据集信息数据分析
肿瘤与正常组织的基因表达差异分析,使用Wilcoxon Rank Sum Tests分析;使用R语言的survival包单因素Cox生存分析,计算风险比(HR)及其95%置信区间(CI)。采用Sparkle数据库(
https://www.grswsci.top)肝细胞癌组织切片执行空间转录组分析。采用Kaplan-Meier生存分析方法,使用R语言的survival包进行数据分析,利用survfit函数进行Log-rank检验,评估患者生存曲线的差异是否具有统计学意义。使用R包cluster Profiler进行GSEA富集分析。分析数据为TCGA数据库中TCGA-LIHC数据集,以调整后的dp<0.05、错误发现率<0.25和标准化富集评分|NES|>1为显著富集;利用单样本基因集富集分析(ssGSEA)方法,Spearman算法计算基因与免疫细胞浸润水平的相关性。
1.3 主要试剂以及仪器、细胞培养
1.3.1 主要试剂以及仪器
DMEM 培养基、1640培养基;胎牛血清(Viva Cell);胰蛋白酶(Solarbio);PBS(Biosharp);佛波酯PMA(LIANKE);ANXA2蛋白、ANXA5蛋白、肯帕罗酮、KLF4抗体(MCE);抗ANXA2抗体、抗ANXA5抗体、抗ANXA11抗体、抗CD206抗体(ProteinTech);RNA提取试剂盒(Omega);逆转录试剂(ThermoFisher);荧光定量PCR试剂盒(Takara);CD206-FITC、 PD-L1-PE(Biolegend);流式细胞仪(Agilent);荧光显微镜(Olympus)。
1.3.2 细胞培养
MHCC-97H细胞(普诺赛),使用含10%胎牛血清的DMEM培养,并收集肿瘤培养细胞上清(TCM)。THP-1细胞培养:THP-1细胞(普诺赛),使用含有10%胎牛血清的1640培养基培养;使用100 ng/mL佛波酯PMA诱导贴壁后,加入TCM共培养48 h。
1.4 免疫荧光检测
本研究所用病理组织切片取自石河子大学第一附属医院,且研究方案已通过伦理委员会审批(伦理批号:KJ2025-290-01)。石蜡切片经脱蜡复水,抗原修复等步骤后加入5%山羊血清封闭;分别加抗CD206抗体(1∶200)、抗ANXA2抗体(1∶100)、抗ANXA5抗体(1∶100)以及抗ANXA11抗体(1∶100),4 ℃孵育过夜;PBST洗脱后,二抗室温避光孵育1 h;再次使用PBST洗脱后,加入含DAPI的封片剂,5 min后置于荧光显微镜观察。
1.5 实时荧光定量PCR
提取总RNA,逆转录为cDNA;取2 μL cDNA,加10 μL SYBR Green Mix、上下游引物各1 μL,补水至20 μL;以2
-ΔΔCt计算相对表达量,使用β-actin作为内参。引物序列如
表1(塞维尔生物合成)。
1.6 流式细胞术
使用细胞刮,收集THP-1细胞,使用PBS洗涤;加入表面抗体PD-L1-PE,4 ℃避光30 min;洗涤,固定破膜后,加入CD206-FITC室温避光30 min;PBS洗两遍后,加入300 μL PBS重悬,上机检测。设同型对照与单染管。
1.7 统计学分析
所有分析均使用R语言(4.3.3版本)和GraphPad Prism9进行。两组间比较采用t检验,多组间的比较则采用方差分析;P<0.05认为差异具有统计学意义。
2 结果
2.1 ANXA2、ANXA5和ANXA11在肝细胞癌中表达上调,与肝细胞癌预后不良相关
对TCGA数据库收录的肝细胞癌(LIHC)数据集进行分析,结果显示ANXAs家族表达特征:ANXA1/2/4/5/6/7/9/11/13在肿瘤组织中上调,而ANXA3、ANXA8和 ANXA10 mRNA水平降低(
P<0.05,
图1A、B)。单因素Cox回归进一步筛选出4个与总生存(OS)密切相关的成员:ANXA2、ANXA5、ANXA10和ANXA11。其中,ANXA2、ANXA5、ANXA11的风险比(HR)均大于1,提示其高表达是预后不良的危险因素;而ANXA10的HR小于1,表明其上调可能降低死亡风险,发挥潜在的保护作用(
P<0.05,
图1C)。PPI网络分析证实,ANXA家族成员在分子功能层面存在相互作用(
图1D)。Kaplan-Meier生存分析显示,ANXA2、ANXA5和ANXA11高表达的HCC患者总体生存率降低于低表达患者(
P<0.05,
图1E)。
2.2 肝细胞癌组织恶性区域的ANXA2、ANXA5和ANXA11转录水平上调
对GSE203612数据集中LIHC8肝细胞癌组织切片执行空间转录组分析,对形态学表达矩阵进行聚类,结合肿瘤细胞的高拷贝数变特征,鉴定出恶性区域和非恶性组织区域(
图2A、B)。采用了Seurat包中的SpatialFeaturePlot函数,对ANXA2、ANXA5和ANXA11在肝细胞癌组织的表达水平分析,结果显示ANXA2、ANXA5和ANXA11在肝细胞癌中表达上调(
图2C~H,
P<0.05)。
2.3 ANXA2、ANXA5和ANXA11 mRNA水平与肝细胞癌巨噬细胞浸润水平呈正相关
对TCGA-LIHC转录组数据进行ssGSEA分析,结果显示,ANXA2、ANXA5和ANXA11与巨噬细胞、NKCD56bright、Th2细胞等浸润呈正相关,与Th17细胞、细胞毒T细胞的浸润呈负相关(
P<0.05,
图3A、C、E)。使用GSE203612数据集中LIHC8肝细胞癌组织切片,对ANXA2、ANXA5和ANXA11在肝细胞癌组织的表达水平与肝细胞癌组织免疫细胞浸润情况分析,分析显示ANXA2、ANXA5和ANXA11与巨噬细胞浸润水平呈正相关(
P<0.05,
图3B、D、F)。
2.4 ANXA2、ANXA5和ANXA11与免疫抑制及黏附相关基因的mRNA水平呈正相关
正相关基因TOP30显示,ANXA2、ANXA5和ANXA11与ITGAM、LIF和FSCN1等参与黏附、胞外基质降解以及TGF-β1、CXCL8等基因的mRNA水平呈正相关(
P<0.05,
图4A~C)。PPI互作网络显示,ANXA2、ANXA5与上述关键基因(ITGAM、FSCN1、MMP2、LIF、TGFβ1、CXCL8、JAK3等)形成紧密互作,而ANXA11仅与ANXA2存在直接关联(
图4D)。
2.5 ANXA2、ANXA5和ANXA11与IL-4/IL-13及PD-1信号通路显著正向富集
生物学功能(BP)分析表明,ANXA2、ANXA5和ANXA11与巨噬细胞对集落刺激因子的响应、白细胞介导的炎症反应等生物学过程密切相关(
P<0.05,
图5A)。基因富集分析(GSEA)显示,ANXA2、ANXA5和ANXA11与IL-4/IL-13、IL-10和PD-1等信号通路正向富集(
P<0.05,
图5B)。
2.6 ANXA2、ANXA5与M2型巨噬细胞标志CD206在肝细胞癌组织中共表达
免疫荧光显示,ANXA2、ANXA5与M2型巨噬细胞标志CD206在肝细胞癌组织中表达上调,且与CD206共定位;而ANXA11在肝细胞癌组织以及癌旁组织表达水平无差异(
图6A~C)。使用人肝癌细胞MHCC-97H细胞系的条件培养基TCM,处理巨噬细胞THP-1,ANXA2、ANXA5表达均上调,而ANXA11表达无差异(
P<0.05,
图6D)。
2.7 ANXA2与ANXA5诱导巨噬细胞M2型标志CD206表达水平上调
与对照组相比,巨噬细胞CD206、PD-L1和IL-10表达上调(
P<0.001,
图7A~C)。为探索巨噬细胞ANXA2、ANXA5上调对巨噬细胞表型的影响,在巨噬细胞培养时分别加入ANXA2、ANXA5纯化蛋白,ANXA2和ANXA5可诱导巨噬细胞IL-10 mRNA水平升高(
P<0.001,
图7D),而ANXA5处理组巨噬细胞IL-10 mRNA有升高趋势,但差异无统计学意义图(
图7D)。ANXA2和ANXA5均能诱导巨噬细胞M2型标志物CD206表达,且ANXA2和ANXA5联合使用诱导的M2型标志物CD206表达高于单独使用(
P<0.001,
图7E),但ANXA2和ANXA5对PD-L1表达均无影响(
图7F)。
2.8 肝癌细胞培养上清通过促进巨噬细胞KLF4表达上调ANXA2、ANXA5水平
ANXA2、ANXA5 mRNA水平与KLF4 mRNA呈正相关(
图8A)。使用肝癌细胞TCM处理THP-1细胞,KLF4表达上调,而经KLF4抑制剂肯帕罗酮处理后,KLF4表达明显受到抑制(
P<0.05,
图8B)。进一步研究发现,使用KLF4抑制剂可有效逆转TCM诱导的ANXA2、ANXA5表达上调(
P<0.05,
图8C),同时显著降低巨噬细胞CD206和PD-L1的表达(
P<0.001,
图8D)。
3 讨论
本研究通过生物信息学分析,揭示ANXA2、ANXA5和ANXA11在肝细胞癌中表达上调,与肝细胞癌不良预后相关。免疫细胞浸润和空间转录组学分析显示ANXA2、ANXA5和ANXA11与巨噬细胞浸润呈正相关,且与PD-1、IL-4/IL-13信号通路显著正向富集。免疫荧光结果显示,ANXA2、ANXA5与M2型巨噬细胞标志CD206在肝细胞癌组织中共表达,而ANXA11在肝细胞癌组织与癌旁组织表达水平无差异;ANXA2、ANXA5纯化蛋白可诱导巨噬细胞M2型标志CD206表达上调,但是不影响PD-L1的表达。Krüppel样因子4(KLF4)作为重要的转录因子,是诱导和维持巨噬细胞M2样表型的关键分子
[18, 19]。此外,本研究发现肝癌细胞TCM可通过上调KLF4表达,促进巨噬细胞ANXA2、ANXA5表达。但肝癌细胞TCM未能诱导巨噬细胞ANXA11表达上调;尽管ANXA11与PD-1、IL-4/IL-13信号通路相关显著正向富集;但PPI互作网络显示,ANXA11与参与巨噬细胞黏附、极化及免疫抑制功能的基因(如ITGAM、MMP9/MMP2、LIF、TGF-β1、CXCL8、JAK3等)不存在紧密互作。以上结果提示,肝癌细胞TCM可能通过“KLF4-ANXA2/5-CD206”轴驱动肿瘤相关巨噬细胞向M2表型极化,从而塑造免疫抑制微环境;而非通过ANXA11诱导巨噬细胞向M2型极化。
基于免疫细胞浸润和空间转录组学的生物信息学分析显示,ANXA2、ANXA5和ANXA11与巨噬细胞浸润呈正相关。肿瘤微环境中的TAMs与肝细胞癌进展和耐药密切相关
[8, 9];免疫荧光结果显示ANXA2、ANXA5均与CD206共定位,这表明ANXA2、ANXA5可能与TAMs趋化以及免疫抑制性表型相关。体外补充ANXA2、ANXA5蛋白可促进CD206表达上调,而不影响PD-L1,说明ANXA2和ANXA5可能通过IL-4/IL-13通路驱动M2型巨噬细胞极化,这与IL-4/IL-13信号正向富集的结果一致。有研究显示,基质金属蛋白酶28(MMP28)可与ANXA2相互作用,驱动TAMs的浸润和M2型极化
[20];在糖尿病模型中,重组ANXA5蛋白促进巨噬细胞从M1向M2表型转变
[21]。多项证据表明,ANXA2与ANXA5可能通过诱导巨噬细胞向M2型极化而发挥相关生物学功能。
KLF4在部分肿瘤中被认为是抑癌基因
[22, 23],但在HBV相关HCC中,HBV X蛋白可通过激活KLF4转录,进而上调HMGB1的表达促进肝细胞癌进展
[24]。此外有研究证实,KLF4是诱导和维持巨噬细胞M2样表型的关键分子
[18]。本研究发现肝癌细胞通过促进巨噬细胞KLF4表达,进而上调巨噬细胞M2型标志CD206以及PD-L1的表达。此外,抑制KLF4,可明显逆转肝癌细胞培养上清诱导的巨噬细胞ANXA2、ANXA5和M2型标志CD206的表达。相关性分析显示,ANXA2、ANXA5 mRNA水平与KLF4 mRNA呈正相关。据此,推测肝癌细胞培养上清可能通过KLF4,驱动巨噬细胞ANXA2、ANXA5的表达上调,进而巨噬细胞M2型标志CD206表达增加。
尽管分析显示,ANXA2、ANXA5与PD-1信号通路正向富集,但体外实验未发现ANXA2、ANXA5对巨噬细胞PD-L1的诱导作用;巨噬细胞上PD-L1的表达主要受IFN-γ-STAT1信号轴调控
[25, 26],据此推测,ANXA2与ANXA5不影响巨噬细胞对IFN-γ的应答,而是通过IL-4/IL-13途径促进巨噬细胞向M2型极化,进而抑制T细胞的功能。尽管本实验证实,KLF4促进巨噬细胞上调ANXA2、ANXA5的表达,但其机制尚未完全阐明。此外,后续研究将通过筛选靶向ANXA2、ANXA5的小分子药物以及体内外实验验证,为肝细胞癌联合免疫治疗提供新方案。
综上,ANXA2与ANXA5通过驱动肿瘤相关巨噬细胞向M2型极化,参与塑造肝细胞癌免疫抑制性微环境,有望成为HCC预后评估及联合免疫治疗的新型分子靶点。
国家自然科学基金(82573336)