TIMP1与EFEMP1在直肠癌中的表达水平及预后意义

2020-05-09谭莹练国达陈少杰李佳佳陈尚祥黄开红陈茵婷

岭南现代临床外科 2020年2期

谭莹，练国达，陈少杰，李佳佳，陈尚祥，黄开红，陈茵婷

直肠癌是常见的恶性实体肿瘤，肿瘤的转移常使患者失去完全切除肿瘤的手术机会，降低了直肠癌的预后［1］。针对肿瘤转移机制的研究十分重要。已知组织金属蛋白酶抑制剂-1（tissue inhibitor of metalloproteinases-1，TIMP1）存在两种作用模式，一方面，TIMP1通过结合膜受体CD63激活胞内通路，从而发挥生物学功能［2-4］；另一方面，TIMP1作为金属蛋白酶（metalloproteinases，MMPs）的抑制剂被人们熟知，TIMP1与MMPs之间的动态平衡，维持着细胞外基质的稳定，与肿瘤转移密切相关［5，6］。有文献报道，胃癌中，TIMP1 与 SERPINE1相关，而SERPINE1与EMT密切相关，TIMP1影响SERPINE1从而影响肿瘤转移［7］。TIMP1在多种肿瘤中高表达，且因其与预后密切相关，还可作为肿瘤预后的预测因子［8-11］。然而，TIMP1在直肠癌中的表达水平及是否能预测直肠癌转移与预后研究较少。

EFEMP1是纤维蛋白糖蛋白家族中的一员，其主要功能是维护基膜稳定性及细胞外基质完整性。有文献报道，在肺癌、肝癌和鼻咽癌中，EFEMP1可抑制肿瘤的侵袭和转移，其表达水平与肿瘤预后相关，可作为肿瘤的预测因子［12-14］。在宫颈癌和骨肉瘤中，EFEMP1可通过诱导EMT从而促进肿瘤转移［15，16］。然而 EFEMP1在直肠癌中的研究相对较少。

本研究旨在探讨TIMP1在直肠癌中的表达水平，以及TIMP1表达量在不同临床分期间的差异，探讨TIMP1与直肠癌预后的关系。探讨与TIMP1相关的差异基因的分子功能及生物学功能。从而探讨TIMP1与EFEMP1的相关性，以及EFEMP1在直肠癌中的表达水平及与预后的关系。

1 材料和方法

1.1 数据来源

直肠癌癌旁组织与癌组织各个基因表达量及临床资料来自于癌症基因组图谱（the cancer genome atlas，TCGA）。

1.2 数据分析方法

分析正常组织与直肠癌组织不同临床分期间的TIMP1的表达差异，分析TIMP1表达水平与直肠癌患者总体生存时间与无瘤生存时间之间的关系。利用 Decenter，参数设置为 log2FC＞1，FDR＜0.05，筛选正常组织与直肠癌组织的差异基因，绘制火山图。根据TIMP1的表达量，将直肠癌患者分为TIMP1高表达组与TIMP1低表达组，利用Decenter，参数设置为log2FC＞1，FDR＜0.05，筛选两组的差异基因，绘制火山图与热图。将差异基因进行基因富集分析，利用FUNRICH进行细胞成分（cellular component，CC），分子功能（molecular function，MF），生物过程（biological process，BP），生物通路（biological pathway）分析。寻找与TIMP1功能可能密切相关的基因，利用Origin 2018分析TIMP1与EFEMP1的相关性。分析EFEMP1表达水平与直肠癌患者总体生存时间与无瘤生存时间之间的关系。

1.3 统计学处理

TIMP1的表达水平经四分位数分成四等份并绘制箱式图进行比较，数据以中位数（四分位间距）表示，采用T检验与秩和检验。使用Kaplan-Meier生存分析来绘制OS及PFS曲线，并采用Logrank检验比较不同表达水平的生存差异从而得出P值。

2 结果

2.1 TIMP1在正常组织与直肠癌中的表达情况及与预后的关系

如图1A所示，TCGA数据库中10例直肠正常组织TIMP1表达量为51.03（42.58，59.40），165例直肠癌组织中TIMP1表达量为174.39（130.71，242.09）。直肠癌组织中TIMP1表达量高于正常组织（P＜0.05）。如图1B所示，30例Ⅰ期患者TIMP1表达量 152.48（107.13，203.03），50例Ⅱ期患者TIMP1表达量 176.84（130.73，235.22），51例Ⅲ期患者 TIMP1表达量 195.35（133.00，296.44），24例Ⅳ期患者 TIMP1表达量 185.65（144.92，223.03）。Ⅰ、Ⅱ、Ⅲ、Ⅳ期直肠癌组织中TIMP1水平均明显高于正常组织（P＜0.05），而各分期之间TIMP1表达水平无统计学意义（P＞0.05）。生存曲线分析显示，TIMP1高表达患者OS与PFS均明显低于低表达患者（P＜0.05）（图1C、D）。

2.2 TIMP1相关的差异基因筛选及基因富集分析

将患者分为TIMP1高表达组与低表达组，筛选出差异基因，随后绘制火山图（图2A）。将差异基因进行基因富集分析。细胞成分分析结果显示，差异基因大多表达细胞外蛋白，与细胞外基质密切相关（图3A）。分子功能分析结果显示，差异基因表达蛋白主要参与构成细胞外基质的组成结构，细胞黏附分子活动，细胞生长活动，蛋白抑制活动，钙离子结合等（图3B）。生物过程分析结果显示，差异基因表达蛋白主要参与细胞生长过程，信号传导过程等（图3C）。生物通路分析结果显示，差异基因主要与EMT通路相关（图3D）。在TIMP1高表达组与低表达组的差异基因中，选取与EMT相关的23个基因，绘制热图（图2B）。筛选出差异基因筛选在正常组织与直肠癌组织中存在表达差异的基因，绘制火山图（图2C）。综合筛选出9个，在正常组织与直肠癌组织中存在表达差异，且与TIMP1相关的EMT相关基因：EFEMP1、Fibulin-1（FBLN1）、生长停滞特异性1（growth arrest specific 1，GAS1）、双糖链蛋白多糖（biglycan，BGN）、视黄酸受体应答蛋白 2（Retinoic acid receptor responder protein 2，RARRES2）、泛素羧基末端水解酶L1（Ubiquitin carboxy-terminal hydrolase L1，UCHL1）、前列环素合酶（prostacyclin synthase，PTGIS）、耳蜗背核（the dorsal cochlear nucleus，DCN）、细胞色素P450 1B1（Cytochrome P450 1B1，CYP1B1）。

2.3 EFEMP1与TIMP1的相关性及EFEMP1与直肠癌预后的关系

如图4A所示，EFEMP1随着TIMP1表达量升高而降低，EFEMP1与TIMP1呈正相关相关（r=0.2795，P＜0.05）。如图4B所示，EFEMP1高表达患者总体生存时间低于低表达患者（P＜0.05）。如图4C所示，EFEMP1高表达患者无瘤生存时间低于低表达患者（P＜0.05）。

3 讨论

本研究结果显示直肠癌组织中TIMP1表达水平高于正常组织，不同临床分期直肠癌组织均高于正常组织。TIMP1高表达组预后差于TIMP1低表达组，验证了此前报道的TIMP1为结直肠癌的标记物。基因富集分析显示，根据TIMP1表达量筛选出的差异基因，其功能与细胞外基质的稳定性及EMT相关。EMT指上皮来源的肿瘤细胞发生基因或表观遗传水平改变，上皮标志物表达下降，间质标志物表达上调，失去细胞间紧密连接和细胞极性，获得间质细胞特性，形态转变为梭形，迁移和侵袭能力增强。一方面，TIMP1可阻断MMPs的功能，从而维持细胞外基质的稳定性，进而抑制肿瘤转移。另一方面，TIMP1还可与SERPINE1等EMT相关分子相互作用，从而影响肿瘤细胞EMT，促进肿瘤的转移。

在根据TIMP1表达量筛选出的差异基因中，进一步筛选出在正常组织与直肠癌组织中存在表达差异的基因。最后得到9个基因：EFEMP1、FBLN1、GAS1、BGN、RARRES2、UCHL1、PTGIS、DCN、CYP1B1。其中，EFEMP1与EMT密切相关，一方面，在骨肉瘤等肿瘤中，敲低EFEMP1可上调E-钙黏素（E-cadherin）表达，抑制N-钙黏素（N-cad-herin）、波形蛋白（Vimentin）、Snail、Slug和Twist的表达，从而促进肿瘤细胞EMT，促进肿瘤的转移。另一方面，在子宫内膜癌等肿瘤中，EFEMP1与E-cadherin呈正相关，与Vimentin，Snail and β-catenin呈负相关［17］。本研究结果显示，在直肠癌组织中，EFEMP1高表达患者OS与PFS均明显低于低表达患者。说明EFEPM1在直肠癌中可能通过促进EMT起到促进肿瘤转移的功能。EFEPM1与直肠癌预后相关，可将EFEPM1作为直肠癌标志物，预测直肠癌患者预后情况。结果显示TIMP1在直肠癌中高表达，且TIMP1与EFEMP1存在正相关，说明可能存在TIMP1激活EFEMP1的表达从而促进EMT，进而促进直肠癌的转移，此假设有待进一步基础实验研究证实。