癌症是一类复杂且多样的疾病,其共同特征是细胞在体内异常增殖,通常形成肿瘤,并可能扩散到身体其他部位
[1]。癌症的形成涉及基因突变、基因表达调控失常、细胞增殖失控和细胞分化异常等多种因素,虽然存在许多不同类型的癌症
[2],如乳腺癌、肺癌、肝癌、结直肠癌等,但它们通常都表现出细胞的生长和分裂不受正常调控机制的约束。
AHCY(Adenosylhomocysteinase)基因成为癌症研究的一个重要焦点。AHCY基因编码的腺苷同型半胱氨酸酶(S‑adenosylhomocysteine hydrolase,SAHH)在一碳代谢途径中扮演关键角色,负责腺苷同型半胱氨酸(S‑adenosylhomocysteine,SAH)到同型半胱氨酸和腺苷的水解
[3],这个过程直接影响到细胞内甲基化反应的平衡。研究
[4‑5]发现,AHCY在某些癌症类型中表现异常,例如在肝癌、胃癌和结直肠癌等中,AHCY的表达水平可能会改变。
肿瘤微环境(tumor microenvironment,TME)由多种细胞组成,其中浸润性免疫细胞占很大比例
[6],TME在人类癌症的发生和发展中发挥着关键作用,然而,探讨TME中基质和免疫成分的动态调节机制仍然是该领域面临的重要挑战。随着免疫检查点阻断和其他策略(例如治疗性疫苗和工程T细胞)的免疫疗法的发展,肿瘤与免疫细胞的相互作用成为人们关注的焦点
[7]。作为经典抗癌疗法的替代方案,免疫疗法已证明对多种癌症类型有效,并被开发用于重新激活适应性和先天免疫系统,其目标是免疫细胞和肿瘤细胞之间的相互作用
[8],这表明AHCY在癌症中的角色值得进一步研究,以探寻其作为癌症治疗靶点的可能性。
在本研究中,我们通过基于各种数据库的数据挖掘分析,全面分析了AHCY的表达及其与肿瘤浸润免疫细胞(tumor‑infiltrating immune cells,TIIC)和相关免疫标志物的关系,并进一步可视化其在泛癌中的预后情况。结果表明,AHCY可能通过与TIIC的相互作用影响癌症患者的预后。AHCY对泛癌具有致癌作用,AHCY表达升高不利于人类癌症患者的生存时间。综上所述,AHCY不仅可作为评估肿瘤免疫浸润及患者不良预后的生物标志物,同时也是一个具有潜力的癌症治疗靶点
1 资料与方法
1.1 原始数据采集和处理
通过TCGA研究网络(
https://portal.gdc.cancer.gov/)对33种不同肿瘤类型的10 000多名肿瘤患者的大量临床和分子数据进行了概要分析。33种癌症的转录组RNA‑seq 数据提取自TCGA数据库。
1.2 癌症中AHCY的基因组改变
癌症患者AHCY状态的改变是从癌症基因组学在线cBioPortal数据库(
http://www.cbioportal.org/)获取的。AHCY的基因组改变包括拷贝数扩增、深度缺失、具有未知意义的错义突变和mRNA上调。
1.3 癌症中AHCY表达分析
肿瘤和匹配的正常组织之间AHCY差异表达的信息来自TCGA和基因型组织表达(GTEx)项目。GTEx(
http://gtexportal.org)是由美国国立卫生研究院(NIH)共同基金建立的组织库和数据资源。通过整合TCGA及GTEx数据资源,我们探索了AHCY在人类泛癌中的表达情况,选择log
2(TPM+1)转换的表达数据进行绘图。
1.4 AHCY表达与生存分析关系
GEPIA(
http://gepia.cancer-pku.cn)数据库是一个在线平台,用于使用标准处理方法剖析来自TCGA和GTEx项目的RNA测序表达数据。利用GEPIA的“生存”模块评估AHCY表达与癌症预后之间的相关性。GEPIA还提供交互功能,例如根据病理阶段进行分析。
1.5 AHCY表达与免疫的关系
TIMER数据库(
https://cistrome.shinyapps.io/timer/)是一个集成的网络服务器,用于评估不同癌症类型中TIIC的丰度,我们用TIMER数据库研究AHCY表达水平与TIIC丰度之间的关系,TIIC包括B细胞、CD4
+ T细胞、CD8
+ T细胞、中性粒细胞、树突状细胞(DC)和巨噬细胞。
1.6 AHCY表达与免疫细胞浸润之间的关系
越来越多的证据表明肿瘤免疫微环境在癌症的发生发展中发挥着重要作用。为了建立免疫和基质成分的估计比例与AHCY表达的关联,我们使用Sangerbox在线平台来估计TME中基质成分的比例。此外,结果以两种分数的形式展示:ImmuneScore和StromalScore。ImmuneScore或StromalScore估计的分数越高,与免疫或基质的比例呈正相关,即相应的分数越高,TME中相应成分比例越大。
1.7 AHCY表达与肿瘤免疫调节基因的相关性分析
肿瘤‑免疫系统相互作用综合资源库门户网站TISIDB(
http://cis.hku.hk/TISIDB/index.php)用于探索不同肿瘤中特定基因与免疫的相互作用。通过TISIDB数据库探讨了AHCY表达与免疫调节剂之间的关系(以
P<0.05为分界值)。Spearman检验用于分析其相关性。
1.8 AHCY表达与基因组异质性分析
肿瘤突变负荷(tumor mutation burden,TMB)正在成为一种新的、意义深远的生物标志物,用于预测免疫治疗效果,并以每个DNA兆碱基的突变总量计算,其中检测到的变异被定义为跨碱基的插入、碱基替换或删除。微卫星不稳定性(microsatellite instability,MSI)是一种分子肿瘤表型,是指短串联重复DNA束中核苷酸的自发丢失或增加。利用Sangerbox (
http://www.sangerbox.com/tool)数据库分析了AHCY基因表达与TMB和MSI之间的关系。
1.9 蛋白质‑蛋白质相互作用网络构建
GeneMANIA(
http://www.genemania.org)是一个交互式、直观的网站,用于构建蛋白质‑蛋白质相互作用(PPI)网络,生成有关基因功能预测的假设并检测具有相似功能的基因。该网络集成算法具有以下生物信息学方法:物理相互作用、共表达、共定位、基因富集分析、遗传相互作用和网站预测。在本研究中,GeneMANIA用于AHCY的PPI分析。
1.10 统计学分析
使用t检验分析来自TCGA和GTEx数据库的基因表达数据。使用Spearman相关分析在TIMER数据库中评估相关分析。通过Spearman相关性评估AHCY表达与免疫细胞丰度评分之间的相关性。所有分析均使用加载R软件包(ggplot2、circlize、clusterProfiler、DOSE和richplot)的R软件(版本3.5.1,www.r‑project.org)进行,以可视化结果。P<0.05被认为差异具有统计学意义,为数据分析提供了可信度。
2 结果
2.1 AHCY在癌症中的mRNA表达和遗传改变差异
由于AHCY在一碳代谢途径及细胞内甲基化中的潜在作用,它可能代表癌症诊断的重要新靶点或生物标志物。为了探究AHCY表达是否与癌症相关,我们评估了不同肿瘤和邻近正常组织中的AHCY表达。TCGA和GTEx数据库的数据显示,与在正常组织中相比,AHCY的mRNA在肾上腺皮质癌(ACC)、膀胱尿路上皮癌(BLCA)、乳腺浸润癌(BRCA)、宫颈鳞癌和腺癌(CESC)、胆管癌(CHOL)、结肠癌(COAD)、弥漫性大B细胞淋巴瘤(DLBC)、食管癌(ESCA)、多形成性胶质细胞瘤(GBM)、头颈癌(HNSC)、肝细胞癌(LIHC)、肺腺癌(LUAD)、肺鳞癌(LUSC)、卵巢癌(OV)、胰腺癌(PAAD)、嗜铬细胞瘤和副神经节瘤(PCPG)、前列腺腺癌(PRAD)、皮肤黑色素瘤(SKCM)、胃腺癌(STAD)、睾丸生殖细胞肿瘤(TGCT)、甲状腺癌(THCA)、胸腺瘤(THYM)、子宫内膜样癌(UCEC)、子宫肉瘤(UCS)肿瘤组织表达显著升高。这表明它可能在多种肿瘤的发展中充当致癌分子(
图1A)。人们普遍认为基因组突变与肿瘤发生密切相关,为了找出癌症中AHCY的基因组突变,对AHCY进行了比较分析。我们首先使用cBioPortal数据库检查了癌症患者AHCY基因的遗传改变。AHCY的基因改变谱表明,扩增是结肠直肠癌、胃食管癌症、非小细胞肺癌、子宫内膜癌、膀胱癌、卵巢癌、乳腺癌、子宫颈癌、头颈癌、子宫内膜样癌、软组织肉瘤、骨癌、成熟B细胞淋巴瘤、神经胶质瘤、前列腺癌改变的最重要的单一因素之一(
图1B)。
2.2 AHCY在癌症中的预后价值
接下来,我们进一步评估了GEPIA中AHCY对泛癌的预后价值[总体生存(OS)和无疾病生存(DFS)]。AHCY表达升高与OS和DFS较差显著相关,特别是与低表达水平相比,AHCY的高表达水平与ACC(
HR=2.3,
P=0.036)、BLCA(
HR=1.4,
P=0.029)和KIRP(
HR=1.5,
P=0.005 3)、LUAD(
HR=1.8,
P=0.004 9)和SARC(
HR=2.5,
P=0.004 5)较差的OS相关(
图2A),有趣的是,我们还发现了在KIRC(
HR=0.6,
P=0.001)中AHCY的高表达水平与较好的OS相关。与ACC(
HR=2.4,
P=0.014)、BLCA(
HR=1.5,
P=0.015)和KIRP(
HR=2,
P=0.025)中较差的DFS相关(
图2B)。上述数据表明,AHCY的表达水平是影响癌症患者预后的重要因素,但其相关性可能因肿瘤类型的不同而存在差异。
2.3 癌症中AHCY表达与免疫浸润水平的相关性
TIIC是调节多种癌症发生和进展的复杂微环境的重要组成部分
[9]。肿瘤浸润淋巴细胞的数量和活性状态是癌症患者生存时间的重要预测指标
[10]。因此,我们探讨了免疫浸润与AHCY表达之间的相关性。通过TIMER数据库探索AHCY表达与免疫浸润水平的系数来确定AHCY表达是否与各种癌症中的免疫浸润水平相关。结果表明,AHCY表达与17种癌症类型中的B细胞、21种癌症类型中的CD4
+ T细胞、17种癌症类型中的CD8
+ T细胞、18种癌症类型中的中性粒细胞、11种癌症类型中的巨噬细胞、15种癌症类型中的树突状细胞浸润水平显著相关(
图3)。
免疫评分(ImmuneScore)与基质评分(StromalScore)是TME分析中的两个重要量化指标,基于基因表达数据(如RNA‑seq或芯片数据)计算,用于评估肿瘤组织中免疫细胞和基质细胞的相对浸润程度。ImmuneScore是通过分析肿瘤转录组数据中免疫相关基因(如淋巴细胞、巨噬细胞、NK细胞等标记基因)的表达水平,量化肿瘤内免疫细胞的总体浸润程度。高评分提示肿瘤微环境中免疫细胞活跃。StromalScore是基于基质细胞相关基因(如成纤维细胞、血管内皮细胞、ECM成分基因)的表达,评估肿瘤间质(非肿瘤细胞成分)的占比。高评分反映显著的基质重塑或纤维化。我们发现AHCY水平在LUAD、SARC、KIRP、KIRC、BLCA和ACC与ImmuneScore呈负相关(
图4A),在LUAD、SARC、KIRP、KIRC、BLCA和ACC与StromalScore也呈负相关(
图4B)。
2.4 癌症中AHCY表达与免疫标记物集、TMB和MSI之间的相关性
免疫监视在确定各种类型癌症的预后中的重要性已被广泛接受。肿瘤可以通过利用免疫检查点基因(包括PD‑1和CTLA‑4)来逃避免疫反应。为了进一步研究AHCY与不同癌症免疫浸润程度之间的关联,我们分析了AHCY与免疫检查点基因表达之间的相关性。例如在GBMLG中,AHCY表达与ENTPD1、ICAM1、BTN3A1、BTN3A2、CD40、TNFSF4、CD70、IFNG、PRF1、IL2RA、CXCL10、CXCL9、CD80、TNFRSF9、CD27、ITGB2、CD40LG、CD28、ICOS、L4、LAG3、CTLA4、HAVCR2、KIR2DL1、KIR2DL3、PDCD1、SLAMF7、CD274、IDO1、TGFB1、CD276、VEGFA和IL10的表达呈正相关(
图5A)。这些结果表明AHCY的高表达可能在介导免疫逃避中发挥重要作用。此外,AHCY在KIPAN、THYM和KICH中与TMB呈正相关,在COAD和COADREAD中与TMB呈负相关(
图5B)。AHCY与SARC、THYM和UVM中的MSI呈正相关,与GBMLGG、COAD、PRAD、COADREAD、LGG、THCA、KIPAN和LUAD中的MSI呈负相关(
图5C)。AHCY表达与TMB之间的关联在不同癌症类型中存在显著差异,KIPAN、THYM和KICH中AHCY的较高表达与较高的TMB相关。较高的体细胞TMB与更好的OS和癌症患者ICI治疗的最佳亚组相关
[11‑12],所有这些数据共同表明AHCY高表达与癌症免疫广泛相关。
2.5 AHCY在癌症中的PPI网络分析
为了探讨AHCY参与癌症发生发展的潜在作用机制,我们使用GeneMANIA在线工具构建了AHCY的PPI网络,结果如
图6所示。如图所示,AHCY与BHMT(Betaine‑Homocysteine Methyltransferase)具有很强的物理相互作用,BHMT基因编码一种关键酶,即甜菜碱‑同型半胱氨酸甲基转移酶,该酶在人体代谢过程中扮演重要角色,特别是在同型半胱氨酸代谢中
[13]。BHMT主要作用于甜菜碱‑甲基化途径,催化甜菜碱将甲基转移给同型半胱氨酸,生成甲硫氨酸(methionine)
[14]。这个过程对维持体内甲基化反应的平衡非常重要,因为甲基化在DNA、RNA、蛋白质和脂质的代谢中起着关键作用,同时甲基化还在基因表达调控和基因组稳定性中扮演重要角色,BHMT的活性或表达的变化可能影响这些过程,从而与癌症相关联
[15],这恰好与共表达的结果一致。此外,预测AHCY与BCKDK、YOD1、BOP1和FBXL20具有显著相关性。
3 讨论
在本研究中,我们从多方面探索了AHCY在泛癌中的表达、其与肿瘤免疫浸润的关系、以及其作为预后因子的潜在价值。通过综合分析TCGA和GTEx数据库中的数据,我们发现AHCY在多种癌症中显示出异常的表达模式,这与某些癌症的预后和免疫特征存在显著相关性。这一发现为AHCY在癌症生物学及免疫治疗中的潜在作用提供了新的视角。
首先,我们的分析表明,与正常组织相比,AHCY的mRNA表达在多种癌症中均有所升高。这一现象表明AHCY在肿瘤发生和进展中可能扮演积极角色。具体而言,AHCY在膀胱癌、前列腺癌、结肠癌、肺腺癌等多种癌症中显著过表达,提示其可能参与这些癌症的致癌机制。这一发现也引发了关于AHCY功能的进一步思考,AHCY在一碳代谢途径中起着至关重要的作用
[16],调控甲基化反应,如果AHCY过度活跃,可能导致DNA和RNA的甲基化模式异常
[17‑18],从而影响基因表达和细胞命运。这种基因表达的失调可能是癌症发生的重要驱动因素。因此,AHCY在多种癌症中的过表达可能与癌细胞的异常增殖、转移和存活能力有关
[19‑20]。
除了观察到AHCY在多种癌症中的过表达外,我们还发现其表达水平与癌症患者的预后密切相关,在一些癌症中AHCY被确认为预后不良的危险因素。这一结果表明,AHCY的高表达可能预示着癌症的更具侵袭性、更高的转移风险以及更差的患者生存率,这也为AHCY作为潜在生物标志物的研究提供了基础
[20]。它可以用于识别高风险的患者群体,从而为临床治疗和管理策略提供指导。此外,AHCY的表达可能作为预测患者预后的参考指标,帮助临床医生制定更有效的治疗计划。
在这项研究中,我们还通过TIMER分析了AHCY表达与肿瘤免疫浸润之间的关系,发现二者在某些癌症中存在显著关联。免疫浸润在癌症的发生和发展中具有关键作用,因为免疫系统负责识别并摧毁异常细胞
[21]。因此,了解AHCY表达与免疫细胞浸润之间的关系,对于理解AHCY在癌症中的作用至关重要。分析结果显示,在某些癌症中,AHCY表达与多种免疫细胞的浸润呈负相关。这意味着AHCY的高表达可能与免疫系统对癌症的抑制作用减弱有关。特别是在肺腺癌和膀胱癌中,AHCY的高表达与较低的免疫浸润水平有关,这可能导致肿瘤更容易逃避免疫监视。同时,AHCY的表达也与免疫调节基因、TMB和MSI呈负相关。这些结果表明,AHCY的过表达可能影响癌症的免疫微环境
[22],从而影响患者对免疫治疗的反应。
最后,我们利用GeneMANIA数据库分析了AHCY与其他蛋白质之间的相互作用。这些分析揭示了AHCY可能参与的代谢途径和生物学过程,进一步证实了其在一碳代谢、甲基化反应等中的作用
[4,16]。这些结果为我们理解AHCY在癌症中的作用提供了更多线索,并可能有助于识别新的治疗靶点。例如,鉴于AHCY在一碳代谢中的关键角色,药物靶向AHCY或相关途径可能成为癌症治疗的新策略。
综合而言,这项研究揭示了AHCY在多种癌症中的表达异常,以及其与免疫浸润和患者预后之间的显著相关性。AHCY的高表达可能是某些癌症的不良预后指标,同时也可能与免疫逃逸和治疗耐药性相关。这些发现为进一步研究AHCY在癌症中的作用以及其作为治疗靶点的潜力奠定了基础。未来的研究可以探索针对AHCY的靶向治疗,并进一步研究其与免疫治疗反应的关系,以期为癌症治疗提供新的方法和策略。
湖北省自然科学基金资助项目(2019CFB306)