张 婧,陈 洁,关贵文,张 婷,鲁凤民,陈香梅
(北京大学基础医学院病原生物学系暨北京大学感染病研究中心,北京 100191)
肝细胞癌(hepatocellular carcinoma,HCC)是世界范围内最常见的恶性肿瘤之一,其术后5年肿瘤复发转移率高达40%~70%[1]。HCC的致病因素主要包括慢性乙型肝炎病毒(hepatitis B virus,HBV)感染、慢性丙型肝炎病毒(hepatitis C virus,HCV)感染、食物黄曲霉毒素暴露、嗜酒以及非酒精性脂肪肝等[2]。我国的HCC发生多与慢性HBV感染有关,且大约70%~90%的乙肝相关HCC都有肝硬化背景[3],因此,阐明HCC特别是乙肝相关HCC的致病机制将对我国HCC防治具有重要意义。
趋化因子及其受体在肿瘤发生发展中发挥重要作用,趋化因子CXCL10是由干扰素-γ(interferon-γ, IFN-γ)诱导产生的一类趋化因子,CXCR3是CXCL10基因的唯一受体。近年的研究显示,CXCL10/CXCR3信号通路与炎症和肿瘤发生密切相关[4]。本研究通过对癌症基因组图谱(the cancer genome atlas,TCGA)和来自国际肝癌研究所(Liver Cancer Institute,LCI)数据库的HCC表达数据库挖掘,以及采用实时荧光定量PCR(quantitative real-time PCR,qPCR)方法检测45例乙肝相关HCC组织中CXCL10基因和CXCR3基因剪接异构体B(CXCR3B)的表达,分析CXCL10和CXCR3基因的异常表达与HCC患者预后的相关性,探索CXCL10基因和CXCR3基因在HCC发生发展中的作用。
45例原发性HCC患者的癌组织和配对癌旁组织样本取自郑州大学附属肿瘤医院接受肝癌手术切除的患者,HCC诊断标准参照中国2011年《原发性肝细胞癌诊疗规范》[5]。所有HCC患者均为男性,均有慢性HBV感染及肝硬化背景,并经病理学确诊。8例正常肝组织来源于肝血管瘤患者。本研究已获得北京大学生物医学伦理委员会审批(IRB00001052-12088), 所有研究对象均签署知情同意书。
TCGA肝细胞癌基因表达数据下载自美国Broad 研究所基因数据分析中心(http://gdac.broadinstitute.org/),采用FPKM(fragments per kilobase million)法对TCGA转录组测序数据进行标准化,并去除了部分重复检测的病例,最终纳入370例不同病因的HCC组织及50例癌旁组织的数据。
LCI肝细胞癌基因表达数据下载自高通量功能基因组数据GSE14520(http://www.ncbi.nlm.nih.gov/geo),最终纳入204例有随访信息的HCC患者,所有患者均有HBV感染背景。
RNA提取试剂Trizol购自美国Invitrogen公司,逆转录试剂盒购自立陶宛Fermentas公司,实时荧光定量PCR试剂购自美国Super Array Biosciences公司,实时荧光定量PCR仪购自瑞士Roche 公司(Light Cycler 480Ⅱ)。
将约100 mg的HCC、癌旁或正常肝组织分别置于装有液氮的研钵中,充分研磨至细粉状,移至无RNA酶的1.5 mL EP管中。加入1 mL Trizol试剂,按照说明书提取肝组织的总RNA。用NanoDrop2000测定RNA浓度,分装冻存于-80 ℃冰箱中备用。
按照逆转录试剂盒说明书,在20 μL的反应体系内,以5 μg总RNA为模板,用随机引物合成cDNA链。取1 μL逆转录产物用于qPCR实验,PCR扩增条件:95 ℃ 10 min,95 ℃ 30 s,60 ℃ 30 s,72 ℃ 30 s,共40个循环。PCR反应运行结束后,采用Roche软件分析每孔反应的Ct值,以管家基因C-TBP1的表达为内参照,计算CXCL10基因和CXCR3B基因的相对表达量,qPCR引物见表1。
表1 临床标本qPCR引物Table 1 qPCR primer sequences for clinical samples
F, forward primer; R, reverse primer.
使用SPSS 22.0及GraphPad Prism 5软件进行统计学分析。正态计量资料采用均数±标准差进行描述,其中配对资料用配对t检验,非配对资料用非配对t检验,相关性分析用Pearson相关;非正态资料采用中位数和上下四分位数[M(P25,P75)]描述,其中配对资料用Wilcoxon配对检验,非配对资料比较用Mann-WhitneyU检验,相关性分析用Spearman 秩相关。使用Kaplan-Meier法进行生存分析,Log-rank检验用于比较组间生存差异。r≥0.3即认为具有相关性,多重比较时采用Bonferroni法对P值进行校正,P<0.05即认为差异具有统计学意义。
为探索CXCL10及其受体CXCR3基因在HCC组织中的表达情况,我们首先分析TCGA数据库中370例HCC组织和50例癌旁组织中的基因表达,结果显示HCC组织中的CXCL10基因表达水平显著高于癌旁组织(P<0.001,图1A)。在TCGA数据库中,有50例HCC患者同时有HCC和配对癌旁组织的表达数据,分析CXCL10基因在配对HCC及癌旁组织中的表达,结果显示50例HCC组织中CXCL10基因的表达水平显著高于配对癌旁组织(P=0.018,图1B)。进一步采用Kaplan-Meier法分析了CXCL10基因表达水平与HCC患者预后的相关性,结果显示与CXCL10基因低表达(n=185)的HCC患者相比,CXCL10基因高表达(n=185)的HCC患者5年总体生存率升高,但差异无统计学意义(P=0.107,图1C)。
我们采用同样的方法也分析了TCGA数据库HCC组织中CXCR3基因的表达,结果表明370例HCC组织中CXCR3基因表达水平显著高于50例癌旁组织(P<0.001,图1D),且CXCR3基因在50例配对HCC组织中的表达亦显著高于癌旁组织(P=0.037,图1E)。Kaplan-Meier分析结果表明CXCR3基因高表达(n=185)的HCC患者5年总体生存率显著高于CXCR3基因低表达(n=185)的HCC患者(P=0.004,图1F)。
A, the expression of CXCL10 in nonpaired HCC tumor tissues (n=370) and non-tumor tissues (n=50); B, the expression of CXCL10 in paired HCC tumor tissues and non-tumor tissues (n=50); C, the prognostic of CXCL10 in HCC patients; D, the expression of CXCR3 in nonpaired HCC tumor tissues (n=370) and non-tumor tissues (n=50); E, the expression of CXCR3 in paired HCC tumor tissues and non-tumor tissues (n=50); F, the prognostic of CXCR3 in HCC patients. Results were showed as 图1 TCGA数据库HCC与癌旁组织中CXCL10基因和CXCR3基因的表达及预后分析Figure 1 Expression and prognostic of CXCL10 and CXCR3 in HCC tumor tissues and non-tumor tissues in TCGA database
我们进一步分析了LCI数据库中CXCL10基因和CXCR3基因在乙肝相关HCC样本中的表达情况,结果表明HCC组织中CXCL10基因表达水平显著高于癌旁组织(P=0.009,图2A)。Kaplan-Meier分析结果显示,CXCL10基因高表达(n=102)HCC患者的5年总体生存率显著高于CXCL10基因低表达(n=102)HCC患者(P=0.002,图2B),但CXCR3基因在HCC组织中的表达显著低于癌旁组织(P=0.003,图2C)。Kaplan-Meier分析结果所示,CXCR3基因高表达(n=102)HCC患者的5年总生存率高于CXCR3基因低表达(n=102)HCC患者(P=0.014,图2D)。
A, the expression of CXCL10 in HCC tumor tissues and non-tumor tissues (n=204); B, the prognostic of CXCL10 in HCC patients (n=204); C, the expression of CXCR3 in HCC tumor tissues and non-tumor tissues (n=204); D, the prognostic of CXCR3 in HCC patients (n=204). Results were showed as 图2 LCI数据库HCC与癌旁组织中CXCL10基因和CXCR3基因的表达及预后分析Figure 2 The expression and prognostic of CXCL10 and CXCR3 in HCC tumor tissues and non-tumor tissues in LCI database
为了验证TCGA和LCI数据库的结果,我们采用qPCR方法检测了45例乙肝相关HCC中CXCL10和CXCR3基因的表达水平,并以8例正常肝组织作为对照,结果显示CXCL10基因在乙肝相关HCC和配对癌旁组织里的表达均显著高于正常肝组织(P=0.010,P<0.001),但CXCL10基因表达水平在HCC和配对癌旁组织中的差异无统计学意义(P=1.000,图3A)。同样,CXCR3基因在HCC和配对癌旁组织里的表达水平也显著高于正常肝组织(P=0.004,P<0.001),但其在HCC组织中的表达水平与配对癌旁组织相比差异无统计学意义(P=0.374,图3B)。
A, mRNA expression of CXCL10 in HCC (n=45), non-tumor (n=45) and nomal (n=8) tissues; B, mRNA expression of CXCR3 in HCC (n=45), non-tumor (n=45) and normal (n=8) tissues. Results were showed as M (P25,P75).图3 qPCR方法检测45例乙肝相关HCC样本中CXCL10和CXCR3基因的表达Figure 3 mRNA expression of CXCL10 and CXCR3 of 45 HCC clinical samples with HBV infection background detected by qPCR
我们进一步分析了来自TCGA、LCI数据库以及我们检测的临床样本中CXCL10基因与CXCR3基因表达水平的相关性,结果如图4所示,在HCC组织中,来自TCGA数据库的CXCL10基因与CXCR3基因的表达水平间呈显著正相关,但在LCI数据库以及我们的临床样本中未见明显相关性(图4A),这一现象同样存在于癌旁组织中(图4B)。
A, correlation between the expression level of CXCL10 and CXCR3 in tumor tissues in TCGA datebase, LCI database and clinical samples; B, correlation between the expression level of CXCL10 and CXCR3 in non-tumor tissues in TCGA datebase, LCI database and clinical samples.图4 CXCL10与CXCR3基因表达水平在HCC和癌旁组织中的相关性分析Figure 4 Correlation between the expression level of CXCL10 and CXCR3 in tumor and non-tumor tissues
作为肿瘤细胞赖以生存的场所,肿瘤微环境可直接影响肿瘤细胞的生长、转移及预后[6],因此研究肿瘤微环境的特点及功能可为肿瘤的早期诊断和治疗提供新依据。趋化因子及其受体之间的相互作用是肿瘤微环境的关键环节,在促进肿瘤细胞增殖、招募免疫细胞及重塑细胞外基质等方面发挥重要作用。已有研究表明,CXCL10基因可通过招募表达CXCR3基因的T细胞、B细胞、单核/巨噬细胞和自然杀伤细胞等,参与多种炎症性疾病和肿瘤的发生与发展[7]。CXCL10/CXCR3在胃癌、乳腺癌、肺癌和多发性骨肉瘤组织中呈高表达,且与上述肿瘤的恶性表型密切相关[7],但CXCL10/CXCR3在HCC中的表达及作用鲜见报道。
CXCL10基因主要由成纤维细胞、内皮细胞、肝实质细胞、角化细胞等多种细胞表达。在细胞处于稳态条件下,CXCL10基因呈低表达,但在炎症细胞因子刺激下CXCL10基因表达上调。CXCL10基因的表达主要由IFN-γ诱导,其主要作用是趋化T淋巴细胞、单核细胞、自然杀伤细胞到达炎症部位,发挥抗炎免疫作用。CXCL10基因也是一种血管抑制因子,可招募抗肿瘤T淋巴细胞,发挥潜在的抗肿瘤效应,但肿瘤细胞的自分泌CXCL10信号也能促进肿瘤细胞的增殖、血管生成和转移[7]。Ren等[8]的研究发现,HCC组织中CXCL10基因表达上调,且CXCL10基因高表达与HCC患者的总体生存率呈负相关,但也有研究表明[9-10],HCC组织中CXCL10基因的表达水平与Th1、CD8+T细胞和自然杀伤样T细胞的浸润呈正相关,提示在某些HCC患者中CXCL10基因可能是通过增强抗肿瘤免疫抑制肿瘤的转移和复发,从而利于HCC患者的预后。为了进一步验证组织中CXCL10基因的表达水平与患者预后的关系,我们分析了TCGA和LCI数据库中CXCL10基因的表达情况及其与预后的关系。与Ren等[8]的发现一致,本研究结果显示CXCL10基因在HCC组织中的表达显著高于癌旁组织,但不同的是,CXCL10基因高表达HCC患者的总体生存率显著高于CXCL10基因低表达HCC患者,提示CXCL10基因高表达利于HCC患者的预后。此外,我们也在45例乙肝相关HCC样本中验证了CXCL10基因在HCC组织的高表达,但HCC和癌旁组织中CXCL10基因水平差异无统计学意义,可能与我们所检测癌旁组织均有肝硬化背景有关。
已有研究表明,CXCR3基因在乳腺癌、骨髓瘤、结肠癌等多种肿瘤中异常表达,且与肿瘤的转移和预后相关[9]。Ding等[11]和任颖等[12]的研究表明,HCC组织中CXCR3基因表达水平高于癌旁组织,HCC组织中CXCR3蛋白的阳性率也高于癌旁组织,且CXCR3蛋白阳性患者术后总体生存率低于阴性患者。本研究来自TCGA数据库370例HCC患者的转录组测序数据显示,CXCR3基因在HCC组织中的表达水平显著高于癌旁组织,且CXCR3基因高表达利于HCC患者的预后,但来自LCI数据库乙肝相关肝癌的分析结果与TCGA相反,CXCR3基因在HCC组织中的表达水平低于配对癌旁组织,而CXCR3基因高表达的HCC患者同样预后更好。我们检测了CXCR3基因在45例乙肝相关HCC中的表达,CXCR3基因在HCC组织和癌旁组织中的表达差异不具有统计学意义,但趋势与LCI数据库一致。本研究检测的45例HCC样本和LCI数据库里纳入的204例HCC均有HBV感染背景,TCGA数据库中仅有102例HCC患者有HBV感染,而Ding等[11]和任颖等[12]的研究没有分析HCC患者的慢性乙肝背景,提示不同研究中CXCR3基因表达不一致的原因可能与HCC患者的病因不同有关。与此结果一致的是,CXCL10基因和CXCR3基因表达水平呈正相关的现象仅存在于TCGA数据库,但在LCI数据库和我们检测的样本中这种相关性消失。由于HCC中浸润的淋巴细胞也可以表达CXCR3基因,淋巴细胞浸润程度会影响CXCR3基因mRNA水平的检测。此外,CXCR3基因有3个剪接异构体,分别为CXCR3A、CXCR3B和CXCR3alt。研究表明,不同CXCR3基因亚型的功能及分布不同。CXCR3A激活可诱发细胞的趋化和增殖,而CXCR3B激活则抑制细胞的迁移和增殖,促进凋亡[13]。本研究定量检测的是CXCR3B亚型,其在乙肝相关肝癌中的低表达可能与肿瘤的发生有关,但TCGA、LCI数据库以及另外两项研究中均未注明所检测的CXCR3基因亚型,因此,不同研究中的CXCR3基因表达结果不一致也可能与所检测的亚型不同有关。
综上所述,本研究结果证实CXCL10基因在不同病因的HCC组织中均呈高表达,而CXCR3基因在乙肝相关HCC组织中呈低表达,且CXCL10基因和CXCR3基因的高表达利于HCC患者的术后生存,提示CXCL10基因和CXCR3基因在HCC的发生、发展过程中发挥重要的作用。
(志谢:郑州大学附属肿瘤医院肝胆胰外科张玲教授提供肝组织样本)