《Journal of Hepatocellular Carcinoma》:Investigating the Expression and Function of PPP2R1A in Hepatocellular Carcinoma: A Multi-Omics Analysis with in vitro Validation
编辑推荐:
背景:肝细胞癌(HCC)的早期诊断与有效治疗仍面临困难。蛋白磷酸酶2A调节亚基Aα(PPP2R1A)作为PP2A全酶的支架(Aα)亚基,在多种癌症类型中发生改变,但其在HCC中的表达模式、功能意义及转化价值尚未被充分阐明。方法:整合公共转录组、单细胞及空间转录
背景:肝细胞癌(HCC)的早期诊断与有效治疗仍面临困难。蛋白磷酸酶2A调节亚基Aα(PPP2R1A)作为PP2A全酶的支架(Aα)亚基,在多种癌症类型中发生改变,但其在HCC中的表达模式、功能意义及转化价值尚未被充分阐明。方法:整合公共转录组、单细胞及空间转录组数据集(TCGA-LIHC、GTEx、GEO、ICGC-LIHC、HPA、TISCH2)与药物敏感性数据库,刻画PPP2R1A在HCC中的表达及其与生存、生物学通路和肿瘤免疫浸润的关联。通过RT-qPCR、CCK-8、Transwell、流式细胞术细胞周期检测及β-半乳糖苷酶实验,在Huh7细胞中评估PPP2R1A敲低的功能效应。结果:PPP2R1A在HCC组织中相对正常肝组织显著过表达,并在若干其他癌症类型中亦升高。Kaplan–Meier分析显示高PPP2R1A表达与较短的总生存期和疾病特异性生存期相关,但多因素校正后并非独立预后因素(风险比[HR] = 1.051,p = 0.789)。PPP2R1A表达与细胞周期、MYC、EGFR、MAPK和WNT信号相关转录程序呈正相关,与p53通路活性呈负相关。在单细胞和空间分辨率上,PPP2R1A富集于恶性细胞和肿瘤区域,其表达与若干计算推断的免疫细胞群丰度相关。Huh7细胞中PPP2R1A敲低降低增殖与侵袭能力,诱导G2期阻滞,并降低MYC、CDK1和CCNB1的mRNA表达,而不诱导衰老。药物敏感性筛选提示AT-7519等候选药物值得在PPP2R1A高表达HCC中进一步研究。结论:PPP2R1A在HCC中持续过表达,并至少部分通过细胞周期调控促进恶性表型。上述发现具有假说生成性质,PPP2R1A在HCC中的潜在诊断、预后或治疗价值需专门的前瞻性与功能验证。
研究背景:肝细胞癌(HCC)是全球常见且侵袭性强的恶性肿瘤,多发生于肝硬化或慢性肝炎患者,早期缺乏明显症状,诊断常被延迟,多数患者确诊时已属晚期,且复发率和死亡率高。HCC治疗面临肿瘤异质性、耐药及个体间治疗反应差异等多重挑战。近年来免疫检查点抑制剂(ICBs)等免疫治疗成为晚期HCC的标准治疗选择,但有效应答率仍偏低,多数患者无法获益,因此提高免疫治疗应答率、准确筛选获益人群并减少治疗副作用是亟待解决的关键科学问题。蛋白磷酸酶2A(PP2A)作为一种关键丝氨酸/苏氨酸磷酸酶,参与细胞周期调控、凋亡和DNA损伤修复等过程,其表达异常或功能失调与肿瘤发生发展密切相关;已有研究提示抑制PP2A可诱导新抗原产生、促进细胞毒性T细胞浸润并增强免疫检查点抑制剂应答,可将"冷肿瘤"转化为"热肿瘤"。PPP2R1A基因编码的蛋白是PP2A复合物的关键调节亚基,其失调与多种肿瘤发生密切相关,在多种常见肿瘤中表达显著升高,但PPP2R1A调控HCC发生发展的机制研究仍有限,具体功能通路尚未完全阐明。因此研究人员采用生物信息学分析与实验验证相结合的策略,探究PPP2R1A在HCC发生发展中的作用,以期为HCC早期诊断与治疗提供新的理论依据,并为临床免疫治疗的精准化与个体化提供新靶点与策略。
主要关键技术方法:研究整合TCGA-LIHC(367例肿瘤、159例正常肝组织)、GTEx、GEO系列数据集(GSE146115、GSE112790、GSE10141、GSE104580、GSE109211、GSE113996、GSE116174、GSE144269、GSE14520、GSE39791、GSE146773、GSE54236、GSE76427)、ICGC-LIHC(n = 260)、HPA(版本23.0)、TISCH2(LIHC_GSE146115,4例HCC标本共3200个单细胞)及TCGA来源LIHC12切片的10x Visium空间转录组数据,并行PRISM、CTRP、GDSC1、GDSC2药物敏感性分析;生物信息学分析通过Guangreshengwu(GRSWSci)在线平台完成。体外实验采用Huh7人肝细胞癌细胞(购自中国科学院细胞库),以si-PPP2R1A敲低,行RT-qPCR、CCK-8、Transwell、流式细胞术细胞周期检测及β-半乳糖苷酶染色;统计采用单因素方差分析(ANOVA)与t检验。
研究结果:
PPP2R1A在HCC组织中上调并在ROC分析中可区分HCC与正常组织:泛癌分析显示PPP2R1A在多数肿瘤类型中相对正常组织升高,提示其上调并非HCC特异;在TCGA-LIHC与GTEx合并数据集(367例HCC肿瘤与159例正常肝样本)中,HCC组织PPP2R1A mRNA表达显著高于正常肝组织(p < 0.0001);PPP2R1A在HPA数据库的癌细胞系(包括HCC来源细胞系)中广泛表达。ROC分析显示PPP2R1A在TCGA-GTEx合并数据集具较强区分能力(AUC = 0.78),在单独TCGA队列中为中等区分能力(AUC = 0.947),但合并数据集的AUC可能部分反映TCGA与GTEx样本的组织来源与批次构成差异,结果需谨慎解读,本身不足以确立临床诊断价值。
PPP2R1A与HCC患者临床预后相关:在TCGA-LIHC队列(n = 367,其中365例具完整生存数据)中,Kaplan–Meier分析显示PPP2R1A高表达组(红曲线)预后显著差于低表达组(蓝曲线),包括总生存期(OS)(log-rank p = 0.009)、疾病特异性生存期(DSS)(p = 0.011)、无进展间期(PFI)(p = 0.071)和无病间期(DFI)(p = 0.141)。多因素Cox回归校正年龄、性别、肿瘤分级、TNM分期、T分期、N分期和M分期后,PPP2R1A表达并非HCC患者的独立预后因素(HR = 1.051,95% CI:0.730–1.514,p = 0.789);单因素Cox回归则显示TNM分期(HR = 1.671,95% CI:1.361–2.050,p < 0.001)、T分期(HR = 1.683,95% CI:1.404–2.018,p < 0.001)和M分期(HR = 4.009,95% CI:1.260–12.759,p = 0.019)为显著预后因素。
PPP2R1A表达与HCC关键致癌信号通路相关:基于TCGA-LIHC队列(高表达组184例、低表达组183例)的差异表达分析显示,高、低表达组间多个基因集存在显著差异,尤其是MYC和G2/M等细胞周期相关通路;部分代谢通路在PPP2R1A高表达组激活(如氧化磷酸化),另一些受抑制(如赖氨酸降解)。CancerSEA功能状态相关性分析(ICGC-LIHC队列,n = 260)显示,PPP2R1A表达与细胞周期(r = 0.38,p < 3.3e-09)、DNA损伤应答(r = 0.39,p < 1.9e-09)和侵袭(r = 0.31,p < 1.9e-06)功能状态呈强正相关。PROGENy通路验证(TCGA-LIHC,n = 374)显示PPP2R1A表达与EGFR(r = 0.455,p < 1.80e-20)、WNT(r = 0.290,p < 1.12e-08)、MAPK(r = 0.272,p < 7.65e-09)、缺氧(r = 0.163,p < 0.01)和TGFβ(r = 0.156,p < 0.01)通路活性评分正相关,与雄激素(r = ?0.287,p < 1.68e-08)、JAK-STAT(r = ?0.203,p < 7.80e-05)、p53(r = ?0.201,p < 8.83e-05)和雌激素(r = ?0.179,p < 0.001)通路活性评分负相关;这些关联源于转录通路活性推断,属相关性结果,不能证明PPP2R1A直接激活或抑制上述通路。
单细胞RNA测序揭示PPP2R1A在HCC中的表达与定位特征:基于GSE146115数据集(4例HCC标本共3200个单细胞)分析肿瘤微环境(TME),UMAP降维识别出恶性细胞、单核/巨噬细胞及多种免疫细胞等特征细胞群,PPP2R1A在恶性细胞中的表达高于其他细胞群;跨细胞类型比较显示恶性细胞PPP2R1A表达显著高于其他细胞类型,柱状图进一步表明PPP2R1A优先表达于恶性细胞,PPP2R1A阳性细胞比例在恶性细胞中显著更高。鉴于该数据集规模有限(4例标本3200个细胞),上述结果为描述性,各主要细胞群细胞数及患者间变异贡献未正式评估,需在更大单细胞队列中确认。
空间转录组揭示高PPP2R1A表达富集于肿瘤区域并与免疫细胞浸润减少相关:对HCC组织切片(TCGA数据库LIHC12,约5000个条形码spots,每个spot直径约55 μm、含1–10个细胞)的10x Visium空间转录组分析显示,PPP2R1A表达与肿瘤细胞分布高度相似,高表达区域与肿瘤细胞所在微区域高度重叠,正常组织区域表达相对较低;肿瘤微环境细胞组成分析确认PPP2R1A在肿瘤中表达。恶性与非恶性组织区域定量分析显示恶性区域PPP2R1A表达显著高于非恶性组织。相关性分析显示PPP2R1A与肿瘤细胞强正相关,与多种免疫及基质细胞类型负相关,即PPP2R1A表达与肿瘤细胞含量正相关,并在空间上与免疫细胞浸润相对较低的区域相关。由于每个Visium spot可能含多个细胞,上述空间相关性反映spot水平细胞组成而非单细胞测量,且源自单一组织切片。基于TIMER2.0的免疫浸润分析(TCGA-LIHC,n = 373)显示PPP2R1A与58种免疫细胞类型显著相关(FDR < 0.05),与Neutrophils_MCPcounter强负相关(rho = ?0.445,p < 0.001),与Th1_cells_xCell强正相关(rho = 0.421,p < 0.001);这些评分来自不同去卷积算法与不同标记基因集,不应直接相互比较,相关性仅描述关联,不能证明PPP2R1A主动调节免疫浸润。
Si-PPP2R1A抑制HCC增殖与侵袭并诱导G2期阻滞:在Huh7细胞中行敲低实验(每个实验三次生物学重复),RT-qPCR验证si-PPP2R1A组PPP2R1A mRNA表达显著低于对照组(p < 0.0001);CCK-8增殖实验显示si-PPP2R1A组细胞增殖在48 h和72 h显著降低(p < 0.0001);Transwell实验显示si-PPP2R1A组侵袭能力显著降低(p < 0.0001);细胞周期分析显示si-PPP2R1A组S期细胞显著减少、G2期细胞显著增多。RT-qPCR显示敲低PPP2R1A显著降低MYC、CDK1和CCNB1的mRNA表达(p < 0.01),提示PPP2R1A敲低诱导G2期阻滞与上述基因mRNA表达降低相关,但机制确立前尚需蛋白水平验证。β-半乳糖苷酶染色显示si-PPP2R1A组较对照组无显著增加,表明PPP2R1A敲低不诱导衰老。
化疗药物敏感性分析与小分子化合物预测:基于TCGA-LIHC队列(n = 367)分析PPP2R1A表达与药物敏感性的相关性。PRISM数据库中AS-703026(r = 0.263,p = 6.24 × 10?9)、trametinib(r = 0.258,p = 1.17 × 10?8)和Ro-4987655(r = 0.255,p = 1.87 × 10?8)与PPP2R1A表达显著正相关,erdafitinib(r = ?0.196,p = 0.0165)、narasin(r = ?0.195,p = 0.0184)和idronoxil(r = ?0.192,p = 0.0247)显著负相关;CTRP数据库中PD318088(r = 0.210,p = 5.80 × 10?9)、selumetinib(r = 0.195,p = 4.48 × 10?8)和saracatinib(r = 0.154,p = 7.29 × 10?8)显著正相关,CCT036477(r = ?0.220,p = 1.18 × 10?8)、PRIMA-1(r = ?0.198,p = 3.03 × 10?8)和ML162(r = ?0.190,p = 8.57 × 10?8)显著负相关;GDSC1数据库中AT-7519(r = 0.141,p = 3.12 × 10?5)、Alectinib(r = 0.125,p = 2.13 × 10?4)和CCT-018159(r = 0.123,p = 2.96 × 10?4)显著正相关;GDSC2数据库中PD0325901(r = 0.132,p = 6.15 × 10?4)、AZD5438(r = 0.109,p = 4.72 × 10?3)和Ulixertinib(r = 0.106,p = 5.67 × 10?3)显著正相关,Nutlin-3a(r = ?0.082,p = 0.034)和AZD4547(r = ?0.078,p = 0.043)显著负相关。通路注释提示部分候选药物(如AT-7519)参与细胞周期调控,与前述PPP2R1A与细胞周期调控的观察一致;上述分析基于公共药物基因组数据库中PPP2R1A表达与药物反应指标的关联,仅提示关联,不能确立疗效、耐药或协同作用。
讨论部分总结:研究人员系统地从表达特征、预后价值、作用机制、细胞功能和药物敏感性等多维度探究PPP2R1A基因,揭示其在HCC进展中的关键作用。需强调PPP2R1A编码PP2A全酶的Aα支架亚基,全文"PPP2R1A表达"指公共数据集中测得的转录本丰度,本研究未测定PP2A酶活性。表达层面,PPP2R1A在HCC组织mRNA水平显著高于正常肝组织,单细胞与空间水平均显示其优先表达于恶性细胞和肿瘤区域;但所分析公共队列为混合病因HCC,HBV、HCV、酒精及代谢性肝病相关HCC的分子与免疫特征差异较大,其对各病因亚组的适用程度尚待确定。预后层面,Kaplan–Meier分析显示高表达与较短总生存期和疾病特异性生存期相关,但多因素Cox回归校正临床病理变量后并非独立预后因素,观察到的生存关联可能主要反映潜在肿瘤分期和其他临床协变量,不应将PPP2R1A单独描述为具有确立的预后价值。机制层面,通路富集、功能实验和微环境分析聚焦细胞周期调控与肿瘤免疫微环境重塑两方面,这些分析识别的是关联而非机制;细胞周期上,高表达组显著富集于MYC信号通路和G2/M细胞周期检查点等增殖相关通路,敲低显著抑制HCC细胞增殖与侵袭并诱导G2期阻滞,β-半乳糖苷酶染色显示对衰老无显著影响,提示其调控肿瘤细胞主要依赖周期破坏而非诱导衰老;免疫微环境上,空间转录组显示PPP2R1A与免疫细胞浸润负相关,高表达区域与肿瘤细胞分布重叠且CD4+ T细胞、CD8+ T细胞和NK细胞等抗肿瘤免疫细胞浸润减少,此空间关联与bulk转录组相关性一致,但相关数据无法推断因果。治疗意义与药物敏感性方面,敲低对恶性表型的抑制作用为PPP2R1A作为候选靶点提供初步实验支持,但不能确立其可被选择性或安全地药物靶向,正常肝组织中功能及正常肝细胞靶向选择性与安全性未评估;药物敏感性分析为联合策略提供线索,如AT-7519与PPP2R1A表达呈强正相关且靶向细胞周期通路,与PPP2R1A周期调控功能高度一致,但这些相关性数据具有假说生成性质,协同作用须经剂量矩阵联合实验及正式协同评估(如联合指数或等效线图分析)直接检验,本研究数据不支持协同结论。耐药与免疫治疗反应方面,PPP2R1A表达与药物AUC值正相关提示高表达与这些药物基因组数据集中药物敏感性降低相关,但未经实验验证,不应视为确立的化疗耐药标志物;免疫浸润分析显示PPP2R1A与中性粒细胞(rho = ?0.445)和Th1细胞(rho = 0.421)等不同免疫细胞群显著相关,方向异质,不应解读为PPP2R1A具有统一免疫抑制作用,Th1细胞评分与中性粒细胞评分来自不同去卷积算法(分别为xCell和MCPcounter)及不同标记基因集,不可直接比较,且bulk转录组评分与肿瘤纯度及基质含量共变;HCC中肿瘤相关中性粒细胞可发挥促瘤与免疫抑制功能,与推断中性粒细胞负相关不等同于抗肿瘤免疫增强,故免疫发现被视为探索性和方法依赖性,本研究不主张PPP2R1A可预测免疫治疗结局。诊断方面,尽管ROC分析中PPP2R1A可区分HCC与正常肝组织,但基于公共组织队列的常规基因水平ROC曲线无法确立临床效用,本研究未在临床环境中评估其诊断或分层价值。局限性包括:发现主要源于公共数据库与体外实验,需独立临床队列验证;下游分子机制未完全明确,通路层面关联为相关性;缺乏动物模型验证,未行体内敲低与药物联合实验;MYC、CDK1、CCNB1仅在转录水平检测,需Western blotting蛋白验证及CDK1活性与G2/M调控因子磷酸化状态评估;需临床样本免疫组化(IHC)验证蛋白表达;免疫微环境关联为相关性,需共培养、细胞因子谱分析或体内免疫分析确立因果;空间分析基于单一LIHC12切片且免疫细胞丰度为计算推断;PPP2R1A与CD8+ T细胞浸润降低关联的分子基础未研究(如CXCL9、CXCL10、CCL5等趋化因子或IDO1等免疫调节酶);本研究测定转录本丰度而非PP2A功能,未评估PP2A复合物组成、磷酸酶活性及相关底物磷酸化状态;单细胞分析基于小型公共数据集(4例HCC标本3200个细胞)并经在线平台处理,细胞群细胞数、质控指标和标本水平效应未系统评估;分析未按HCC病因或病毒状态分层;PPP2R1A未与甲胎蛋白(AFP)等既有临床标志物比较,未评估其加入常规临床病理变量后能否改善判别或校准;功能实验仅在单一Huh7细胞系采用单一siRNA敲低,需在其他HCC细胞系以独立沉默策略(第二种siRNA、shRNA或CRISPR敲除)及PPP2R1A再表达挽救实验确认,且未测试PPP2R1A药理学抑制剂。未来研究将聚焦:在独立、更大HCC队列中验证PPP2R1A表达及临床相关性;在注释完善的HCC队列中行病因与病毒状态分层分析;利用蛋白-蛋白相互作用和转录组测序鉴定下游靶点;在动物模型中评估PPP2R1A靶向策略单用或联合免疫治疗的疗效;利用临床免疫治疗队列或体内模型研究PPP2R1A表达与免疫治疗反应的关系,并在更多HCC细胞系中以过表达或敲低验证其在耐药中的作用。
研究结论(译文):综上所述,该多维分析表明PPP2R1A在包括HCC在内的多种癌症类型中广泛过表达,并优先表达于恶性细胞。高PPP2R1A表达在Kaplan–Meier分析中与较短的总生存期和疾病特异性生存期相关,但并非独立预后因素。PPP2R1A表达与细胞周期及增殖相关转录程序相关,Huh7细胞中PPP2R1A敲低降低增殖与侵袭并诱导G2期阻滞,表明其对恶性表型的功能性贡献至少部分由细胞周期调控介导。上述发现具有探索性和假说生成性质;PPP2R1A在HCC中是否具有诊断、预后或治疗效用,需在独立临床队列和更确定的功能研究中验证。