操利超,张核子,余晨笛,卢晓萍,巴 颖
(深圳市核子基因科技有限公司,广东 深圳 518071)
随着科学的发展,结肠癌在诊断筛查、预后和治疗等方面已取得了长足的进步,但其仍是世界上死亡率最高的恶性肿瘤之一[1]。有研究表明[2-4],肿瘤微环境可能在肿瘤的发生、发展和治疗中发挥着重要的作用。近年来,一些研究从肿瘤免疫微环境中挖掘多种肿瘤分子标志物,用于肿瘤的预后、诊断和免疫治疗[5-8]。CLDN11 是claudin 家族的成员之一,是一种编码蛋白基因。有研究表明,claudin 家族成员在多种肿瘤组织中均出现异常表达情况,且与多种肿瘤的预后相关[9-12]。然而,尚缺少研究全面探讨CLDN11 与多种肿瘤的预后和肿瘤微环境特征之间的关联。基于此,本研究将以结肠癌为例,重点探讨CLDN11 作为潜在的标志物与结肠癌预后、肿瘤微环境特征和免疫治疗的相关性,现报道如下。
1.1 数据来源 结肠癌的mRNA 表达数据集和对应的临床信息、生存数据信息从UCSC Xena 平台(https://xenabrowser.net/datapages/)下载,选择队列为GDC TCGA Colon Cancer(TCGA-COAD)。总计471 个样本,其中肿瘤样本432 个,正常样本39 个,对应的临床信息包括性别、年龄、病理分期、TMN 分期、总生存时间和无病生存时间。
1.2 TIMER2.0 在线分析 使用TIMER2.0(http://timer.cistrome.org/)进行在线分析[13]。打开主页,选择“Exploration”菜单,点击“Gene_DE”选项卡,在文本框中输入“CLDN11”,点击“Submit”按钮,查看CLDN11 在各种肿瘤中表达情况。接着,选择“Immune”菜单,点击“Gene”选项卡,在“Gene Expression”文本框中输入“CLDN11”,在“Immune Infiltrates”下选项中选择所有的免疫细胞类型,点击“Submit”按钮,查看CLDN11 的表达与各种免疫细胞侵润程度的相关性。P<0.05 为判断显著性的标准。
1.3 GEPIA2 在线分析 使用GEPIA2(http://gepia2.cancer-pku.cn/)进行在线生存分析[14]。打开主页,在左侧菜单栏选择“Expression Analysis”,选择子菜单“Survival Analysis”进入生存分析页面。在该页面中,选择“Gene”选项并填入“CLDN11”,“Methods”分别进行“Overall Survival”和“Disease Free Survival(RFS)”分析,“Datasets”每次选择1 个数据集,其它参数默认,点击“Plot”提交分析,即可分析CLDN11的表达水平与不同肿瘤的总生存时间和无病生存时间的相关性。其中,P<0.05 作为CLDN11 是否可以作为独立的预后因子的判断标准。
1.4 CLDN11 的表达与微卫星状态、免疫治疗靶点、结肠癌肿瘤分期的关联 选择TCGA-COAD 队列中的肿瘤样本,根据其微卫星状态,将样本分为微卫星稳定组(MSS)、微卫星高不稳定组(MSI-H)和微卫星低不稳定组(MSI-L)。利用Kruskal-Wallis 检验分析比较这三组中CLDN11 表达水平的差异性。并利用R 软件中的cor.test 分析CLDN11 的表达与CD274、CTLA4、HAVCR2、LAG3、PDCD1 和 TIGIT的表达相关性,显著性判断标准为P<0.05。根据样本的肿瘤分期信息,利用Kruskal-Wallis 检验和Wilcoxon 检验分析CLDN11 的表达与肿瘤TMN 分析、病理分期的相关性。上述分析过程利用R 语言(版本3.6.1,https://www.r-project.org/)完成。
2.1 CLDN11 在泛癌种中的表达情况和预后分析CLDN11 在膀胱癌(BLCA)、乳腺癌(BRCA)、结肠癌(COAD)、头颈癌(HNSC)、肾嫌色细胞癌(KICH)、肾透明细胞癌(KIRC)、肾乳头细胞癌(KIRP)、肝癌(LIHC)、肺腺癌(LUAD)、肺鳞状细胞癌(LUSC)、前列腺癌(PRAD)、直肠癌(READ)、子宫内膜样癌(UCEC)等肿瘤中高表达(P<0.05),而在胆管癌(CHOL)和甲状腺癌(THCA)中低表达(P<0.05),见图1A。对于结肠癌(COAD)、肾乳头细胞癌(KIRP)、肾透明细胞癌(KIRC)和胃癌(STAD),低表达组(Low CLDN11 Group)的总生存时间和无病生存时间等预后指标均表现较好,而对于黑色素瘤(SCKM),高表达组(High CLDN11 Group)的预后表现较好(P<0.05),见图1B~图1K。
图1 CLDN11 在泛癌种中的表达与生存预后分析
图1 CLDN11 在泛癌种中的表达与生存预后分析(续)
2.2 结肠癌样本中CLDN11 的表达与免疫细胞侵润程度的相关性分析 CLDN11 的表达水平与CD4+T细胞(CD4+TCell)、巨噬细胞(Macrophage)、中性粒细胞(Neutrophil)和树突状细胞(Dendritic Cell)的侵润水平正相关(cor>0.1,P<0.05),见图2。
图2 CLDN11 的表达与免疫细胞侵润程度的相关性
2.3 结肠癌样本中CLDN11 的表达与微卫星状态、免疫治疗靶点的相关性分析 CLDN11 在微卫星稳定组和微卫星低不稳定组的表达没有明显差异(P=0.4),而微卫星状态越稳定,CLDN11 的表达水平越高(P<0.05),见图3A。此外,CLDN11 的表达与常见的免疫治疗靶点基因的表达显著正相关(CD274:cor=0.347,P=1.639e-15;CTLA4:cor=0.293,P=2.857e-11;HAVCR2:cor=0.489,P=3.308e-31;LAG3:cor=0.369,P=2.119e-17;PDCD1:cor=0.327,P=8.046e-14;TIGIT:cor=0.407,P=3.549e-21),见图3B~图3G。
图3 CLDN11 的表达与免疫治疗靶点的相关性
2.4 CLDN11 表达水平与在结肠癌肿瘤分期的关联CLDN11 表达量越高,预示着肿瘤的体积和侵润的范围越大(P=0.019),肿瘤发生转移的可能性越高(P=0.0076),与肿瘤相关的淋巴结转移越多(P=0.0043),病理分期程度越高(P=0.013),见图4。
图4 CLDN11 的表达与结肠癌肿瘤分期的相关性
图4 CLDN11 的表达与结肠癌肿瘤分期的相关性(续)
本研究结果表明,CLDN11 在大多数肿瘤中低表达。CLDN11 由于被高度甲基化而导致表达下调,并通过与微管蛋白相关作用,阻断微管蛋白的聚合和细胞迁移活动,从而导致肿瘤的发生[15]。有研究表明[15-17],CLDN11 在多种肿瘤中作为一种潜在的肿瘤抑制基因而发挥作用。此外,本研究发现CLDN11 可以作为多种肿瘤的独立的预后预测因子,且在大多数肿瘤CLDN11 低表达组的预后相对良好。
在结肠癌中,本研究和既往报道的文献[18]均表明CLDN11 低表达的患者预示着总生存期越长;而在乳腺癌中则相反,CLDN11 高表达的患者预后相对较差[19,20]。其次,在结肠癌中,CLDN11 的表达水平越高,CD4+T 细胞、巨噬细胞、中性粒细胞和树突状细胞的侵润水平侵润越明显,这几种免疫细胞在抗肿瘤免疫中有着重要的调控作用。Borst J 等[21]概述了介导CD4+T 细胞调控的细胞动力学和膜受体,可增强肿瘤内细胞毒性T 淋巴细胞(scytotoxic T lymphocytes,CTL)的抗肿瘤活性。有研究表明[22],巨噬细胞是具有高度可塑性的细胞,具有多种功能,包括组织发育和体内平衡、清除细胞碎片、清除病原体和调节炎症反应等,以巨噬细胞为基础,减少抗炎性巨噬细胞并增加促炎(抗肿瘤)巨噬细胞,达到肿瘤免疫治疗的目的。
在肿瘤微环境中,中性粒细胞具有不同的功能,能够对包括癌症在内的各种炎症信号做出应答反应,对肿瘤患者的预后会有一定的影响[23]。树突状细胞是一种异质性白细胞群,由驱动特定类型免疫反应的不同亚群组成,能启动和调节适应性免疫反应,从而达到抗肿瘤免疫反应的目的[24]。因此,通过调控CLDN11 的表达水平,可间接调控抗肿瘤免疫反应,从而达到免疫治疗的目的。此外,在结肠癌患者中,CLDN11 的表达水平越低,微卫星状态越不稳定。有研究表明,微卫星高不稳定患者对免疫检查点抑制剂越敏感,将从免疫治疗中获益[25-27],且CLDN11 的表达与多种常见的免疫治疗靶点基因的表达呈现显著的正相关,进一步表明CLDN11 可作为潜在的免疫治疗靶点。另外,在结肠癌中,CLDN11 的表达水平与肿瘤的增值、转移和侵袭呈显著正相关,表明CLDN11 可作为结肠癌分期等级的预测因子。
然而,本研究也存在一些不足之处:首先,本研究还需要更多的样本队列来进行交叉实验验证;其次,CLDN11 的异常表达对于肿瘤的发生、发展和预后的影响尚待进一步开展机制研究和探讨;最后,CLDN11 作为潜在的免疫治疗靶点是否具有临床价值尚待开展临床实验验证和作用机制研究。
综上所述,CLDN11 在评估结肠癌患者复发风险分层、肿瘤分期、免疫治疗等方面具有潜在的应用价值。