胰腺导管细胞癌组织CSF3R基因表达对患者预后的影响及其机制

胰腺导管细胞癌组织CSF3R基因表达对患者预后的影响及其
机制
方天逸;王宇夫;林轩;王浩;崔云甫
【摘要】目的探讨粒细胞集落刺激因子3受体( CSF3R)基因表达对胰腺导管细胞癌患者预后的影响及其机制.方法从GEO数据库中下载包含CSF3R表达数据的人类胰腺导管细胞癌表达谱,包括132例份胰腺导管细胞癌样本.通过R2分析平台,以Kaplan-Meier法绘制生存曲线,分析CSF3R表达与患者预后的关系,筛选出与
CSF3R表达显著相关的基因并进行KEGG通路分析与GO分析.结果 CSF3R表达
与患者总体生存率呈负相关(P＜0.05),CSF3R高表达患者的总体生存率低于低表达患者(P＜0.05).检测出与CSF3R表达呈显著相关的基因共6 521个,其中正相关基
因3 040 个、负相关基因3 481 个. KEGG分析发现,有27 个高表达基因( C1QA、C1QB、C1QC、 C2、 C3AR1、 C5、 C5AR1、 FCAR、 FCGR2A、 FCGR2B、FCGR3A、FCGR3B、FGG、FPR1、FPR2、FPR3、HLA-DMB、HLA-DQA1、HLA-DQB1、ICAM1、ITGAM、ITGB2、MASP1、MASP2、MBL2、PLG、SELPLG)仅涉及1个信号通路,即金黄色葡萄球菌感染. GO分析发现,CSF3R及其
相关基因主要与细胞对刺激的反应、免疫反应、RAS调控、胰腺发育等相关.结论CSF3R表达与胰腺导管细胞癌患者总体生存率相关,其可能是预测患者预后的潜在
分子标记物.%Objective To investigate the relationship between the expression of granulocyte colony-stimulating factor-3 receptor (CSF3R) gene and the prognosis of patients with pancreatic ductal adenocarcinoma , and to analyze the correlated genes and cellular pathways.Methods The corresponding cancer datasets included 132 samples were extracted from the public database GEO to obtain the
patients’ information and mRNA expression data.The survival curve was drawn by Kap-lan-Meier method using the R2 bioinformatics analysis platform, and the relationship between the expression of CSF 3R and the prognosis of the patients was analyzed.Genes significantly associated with CSF 3R expression were enriched and analyzed for KEGG pathways and GO pathways.Results The datasets showed a significant negative correlation between CSF 3R mRNA expression level and the overall survival rate (P＜0.05).The overall survival rate was higher in patients with high CSF 3R expression than in patients with low CSF3R expression (P＜0.05).In addition, there were 6 521 genes associated with CSF3R expression, during which, 3 040 genes were positively correlated and 3 481 genes were negatively correlated.KEGG analysis showed that 27 gene were highly expressed
(C1QA, C1QB, C1QC, C2, C3AR1, C5, C5AR1, FCAR, FCGR2A, FCGR2B, FCGR3A, FCGR3B, FGG, FPR1, FPR2, FPR3, HLA-DMB, HLA-DQA1, HLA-DQB1, ICAM1, ITGAM, ITGB2, MASP1, MASP2, MBL2, PLG, SELPLG), and they involved only one signaling pathway, staphylococcus aureus infection.GO analysis found that CSF3R and its related genes were mainly related to cell response to stimulation, immune response, RAS regu-lation, and pancreas development.Conclusion CSF3R expression is negatively correlated with the overall survival rate of pa-tients with pancreatic ductal adenocarcinoma, and it may be a potential molecular marker for the study of prognosis of patients.
【期刊名称】《山东医药》
【年(卷),期】2018(058)024
【总页数】4页(P18-21)
【关键词】胰腺导管细胞癌;粒细胞集落刺激因子3受体;基因芯片;总体生存率【作者】方天逸;王宇夫;林轩;王浩;崔云甫
【作者单位】哈尔滨医科大学附属第二医院,哈尔滨150001;哈尔滨医科大学附属第二医院,哈尔滨150001;哈尔滨医科大学附属第二医院,哈尔滨150001;哈尔滨医科大学附属第二医院,哈尔滨150001;哈尔滨医科大学附属第二医院,哈尔滨150001
【正文语种】中文
【中图分类】R435.9
近年来，表观遗传学关于胰腺癌的研究取得了巨大进步，“炎癌转化”假说已经得到广泛认可[1]。

肿瘤组织中的炎性细胞，特别是嗜中性粒细胞是目前研究的热点[2]。

学者们致力于寻找与嗜中性粒细胞高度相关，甚至可在一定程度上代表肿瘤组织内嗜中性粒细胞募集程度的基因。

慢性中性粒细胞白血病是一种以血液中嗜中性粒细胞恶性增多为主的疾病，主要由粒细胞集落刺激因子3受体(CSF3R)基因突变导致[3]。

CSF3R不仅在骨髓中高表达，在胆囊、肺脏、膀胱、肾脏等组织也被检测到不同程度的表达[4]。

目前CSF3R基因在实体肿瘤中作用的研究刚刚起步。

2017年12月～2018年2月，本研究通过公开数据库GEO中的芯片数据，利用R2平台，采用生物信息学技术方法分析了CSF3R与胰腺导管细胞癌患者总体生存率的关系，并分析其相关的作用机制和信号通路，旨在为进一步深入研究
CSF3R基因与胰腺导管细胞癌的关系奠定基础。

1 材料与方法
1.1 材料从GEO数据库(https:///geo)中下载基因表达谱公共数据集GSE21501，此数据集包含CSF3R表达数据的人类胰腺导管细胞癌表达谱，包括132例份胰腺导管细胞癌患者的癌组织标本(30例份标本无生存资料)，分别来自北卡罗来纳大学教堂山分校、西北大学纪念医院、北岸大学、内布拉斯加医学中心和霍普金斯医疗机构等医疗机构。

其中原发性胰腺癌117例，转移性胰
腺癌15例。

T分期：T1期2例，T2期18例，T3期90例，T4期2例；N分期：N0期34例，N1期81例；M分期：M0期116例，M1期15例；分化程度：
低分化5例，中分化55例，高分化51例；切缘阳性16例、阴性63例；新辅助治疗113例。

注：内布拉斯加医学中心的15例份标本缺失T分期、N分期、分
化程度、切缘情况、新辅助治疗等资料；霍普金斯医疗机构的34例份标本缺失切缘情况资料。

1.2 实验方法采用Allprep试剂盒(Qiagen)从快速冰冻肿瘤组织中提取癌细胞RNA，纳升分光光度计定量Bioanalyzer2100评估RNA质量。

利用R2分析平台(http://r
2.amc.nl)进行基因表达情况与生存数据的相关性运算，通过Kyoto Encyclopedia of Genes and Genomes(KEGG)通路分析及Gene Ontology(GO)分析方法对CSF3R的相关基因通路进行初步分析。

生存数据的处理采用Kaplan-Meier法。

P<0.05为差异有统计学意义。

2 结果
2.1 CSF3R表达与患者总体生存率的关系 CSF3R高表达患者62例，低表达患者
40例；CSF3R表达水平与患者的总体生存率呈负相关(P=0.015)，高表达患者的
总体生存率低于低表达患者(P<0.05)。

见图1。

图1 CSF3R表达与患者总体生存率的关系
2.2 与CSF3R相关的基因利用R2分析平台，对数据集GSE21501中132个胰腺导管细胞癌患者的全部基因进行统计学分析，检测出与CSF3R表达呈显著相关的
基因共6 521个，其中正相关基因3 040个、负相关基因3 481个，前20位正
负相关基因见表1。

表1 与胰腺导管细胞癌患者CSF3R表达呈正负相关的前20位基因基因
rPFPR20.730<0.001LILRA30.711<0.001CXCR20.691<0.001AQP90.681<0.00
1C19orf590.681<0.001LILRA50.680<0.001SIRPD0.666<0.001CMTM20.657< 0.001LILRA50.643<0.001TSPAN160.640<0.001ANK3-
0.609<0.001PLA2G12A-0.607<0.001PCP4-0.581<0.001LNX2-
0.581<0.001PDIK1L-0.577<0.001DLG3-0.568<0.001MAP7-
0.568<0.001RPRD1A-0.566<0.001ATP9A-0.565<0.001BEX2-0.557<0.001 2.3 基因通路 KEGG分析发现，与CSF3R相关的6 521个基因中，有1 943个基因涉及不同的信号通路。

当取P<0.001时，27个CSF3R高表达基因(C1QA、
C1QB、C1QC、C2、C3AR1、C5、C5AR1、FCAR、FCGR2A、FCGR2B、FCGR3A、FCGR3B、FGG、FPR1、FPR2、FPR3、HLA-DMB、HLA-DQA1、HLA-DQB1、ICAM1、ITGAM、ITGB2、MASP1、MASP2、MBL2、PLG、SELPLG)仅涉及1个信号通路，即金黄色葡萄球菌感染。

GO分析发现，CSF3R
及其相关基因主要与细胞对刺激的反应、免疫反应、RAS调控、胰腺发育等相关。

见表2。

表2 与胰腺导管细胞癌患者CSF3R相关基因的GO分析GO Path编号生物学过
程基因数P70887细胞对化学刺激的反应692<0.00171310细胞对有机物质的反应962<0.00148583调节对刺激的反应560<0.0016955免疫反应
413<0.0016968细胞防御反应26<0.00146578RAS蛋白信号传导的调控
170<0.0016909吞噬作用71<0.00145321白细胞激活202<0.00138127ErbB
信号通路85<0.00131016胰腺发育34<0.001
3 讨论
胰腺导管细胞癌的恶性程度远高于其他实体肿瘤[5]。

胰腺癌具有特殊的解剖结构以及复杂的肿瘤微环境，肿瘤组织长期处于缺血和低氧环境，肿瘤细胞具有高度的异质性和异常的代谢模式[6]，但这些尚不能充分解释胰腺癌的发病机理。

随着研究的不断深入，胰腺癌组织全基因组测序技术的成熟和分子生物学的长足发展，为更深层次了解胰腺癌的生物学特征提供了条件，推动了胰腺癌早期诊断和特异性靶向治疗的发展。

CSF3R是一种蛋白质编码基因，又名粒细胞集落刺激因子受体，位于1号染色体的p34区，由17 272个碱基构成，有19个外显子。

该基因可以调控粒细胞的生成、分化和信号传导[7]。

CSF3R基因起初被认为在粒细胞的生长和发育中发挥重要作用[8]。

Maxson等[9]发现血液系统肿瘤中CSF3R突变，揭开了其在血液系统肿瘤研究的序幕。

目前认为，CSF3R突变与重度先天性中性粒细胞减少症关系密切，可诱发其发展为急性髓细胞白血病(AML)[10]。

CSF3R突变包括CSF3R
T618I和CSF3R T615A两种类型，常合并SETBP1基因突变[11]。

本研究对选取的芯片数据进行Kaplan-Meier生存曲线分析证实，CSF3R对于胰腺导管细胞癌患者的总体生存期有明显的预测作用，CSF3R高表达患者的总体生存期明显短于低表达患者。

尽管数据库中有30例患者生存资料缺失，但仍提示CSF3R是一个明显的促癌因素。

本研究通过R2分析平台检测出与CSF3R表达显著相关的基因共6 521个，经计算与其高度正相关的部分基因均发挥不同程度的促癌作用。

Hou等[12]报道，FPR2表达水平与胃癌恶性程度、淋巴结转移情况呈正相关；Wang等[13]报道，CXCR2在胃癌细胞的增殖和侵袭中起重要作用。

对这6 521个基因的KEGG分析结果显示，CSF3R高表达时仅涉及金黄葡萄球菌感染信号通路。

Wang等[14]认为，胰腺癌组织中炎性细胞浸润程度越高，患者预后
越差；Bayne等[15]报道，肿瘤来源的粒细胞-巨噬细胞集落刺激因子(GM-CSF)
是胰腺癌组织内炎症和免疫抑制的重要调节剂。

葡萄球菌感染灶内主要的炎性细胞为中性粒细胞。

我们推测CSF3R高表达导致患者总体生存率降低可能与中性粒细
胞的三种杀伤途径有关，即吞噬、脱颗粒和中性粒细胞外网(NETs)。

在吞噬过程中，为了消除和降解微生物，活性氧(ROS)、防御素和丝氨酸蛋白酶在吞噬体中释放。

ROS的持续释放导致由中性粒细胞局部浸润的组织微环境持续缺氧，而慢性
持续性低氧环境在胰腺导管细胞癌的发生、侵袭和转移中发挥重要作用。

Abdol Razak等[16]研究表明，胰腺癌细胞会刺激周围中性粒细胞产生NETs，而NETs
可促进ROS对肿瘤组织及其间质的杀伤作用，进而诱导多种致癌因子的产生。

Geller等[17]研究发现，胰腺癌对吉西他滨产生耐药是由于肿瘤组织中存在可以产生胞苷脱氨酶的特异性细菌；这些细菌可以代谢化疗药物，并募集中性粒细胞，使肿瘤产生耐药性，这也从一个新的角度阐释了细菌、炎症和肿瘤相互促进的关系。

本研究 GO分析结果显示，CSF3R与其相关基因在细胞炎症反应中发挥了一系列
功能，如免疫反应、白细胞激活等。

此外，其激活的ErbB信号通路也可以调控直肠癌细胞的增殖、浸润和转移[18]。

在本研究数据集中CSF3R高表达导致胰腺癌
患者不良预后，该结果为今后研究CSF3R在胰腺癌发展中作用的相关机制提供了
思路。

转录组学的发展使临床对胰腺癌有了新的认识。

转录组学与目前的细胞学研究不同，其虽可在短时间内检测到人类细胞或组织中所有基因的转录表达，但这种通过公共平台的分析方法无法描述细胞间化学相互作用的真实情况。

其次，其统计结果只适用于群体患者，而不适用于单一患者，因此需要更大的样本量进行进一步的研究。

通过公共基因芯片数据进行分析，具有快捷、简便、样本量大的优势。

本研究中CSF3R表达与胰腺癌患者的总体生存率呈负相关，可以为下一步CSF3R与胰腺癌相关分子机制的研究提供参考。

【相关文献】
[1] Zambirinis CP, Pushalkar S, Saxena D, et al. Pancreatic cancer, inflammation, and microbiome[J]. Cancer J, 2014,20(3):195-202.
[2] Wang Y, Fang T, Huang L, et al. Neutrophils infiltrating pancreatic ductal adenocarcinoma indicate higher malignancy and worse prognosis [J]. Biochem Biophys Res Commun, 2018,501(1):313-319.
[3] Zhang H, Means S, Schultz AR, et al. Unpaired extracellular cysteine mutations of CSF3R mediate gain or loss of function[J]. Cancer Res, 2017,77(16):4258-4267.
[4] Fagerberg L, Hallstrom BM, Oksvold P, et al. Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics[J]. Mol Cell Proteomics, 2014,13(2):397-406.
[5] Lavu H, Yeo CJ. Pancreatic ductal adenocarcinoma treatment--the past, present, and future[J]. Semin Oncol, 2015,42(1):4-7.
[6] Nagathihalli NS, Castellanos JA, Shi C, et al. Signal transducer and activator of transcription 3, mediated remodeling of the tumor microenvironment results in enhanced tumor drug delivery in a mouse model of pancreatic cancer[J]. Gastroenterology, 2015,149(7):1932-1943 e1939.
[7] Sheng G, Zhang J, Zeng Z, et al. Identification of a novel CSF3R-SPTAN1 fusion gene in an atypical chronic myeloid leukemia patient with t(1;9)(p34;q34) by RNA-Seq[J]. Cancer Genet, 2017(216-217):16-19.
[8] Wright CR, Brown EL, Della Gatta PA, et al. Regulation of granulocyte colony-stimulating factor and its receptor in skeletal muscle is dependent upon the type of inflammatory stimulus[J]. J Interferon Cytokine Res, 2015,35(9):710-719.
[9] Maxson JE, Gotlib J, Pollyea DA, et al. Oncogenic CSF3R mutations in chronic neutrophilic leukemia and atypical CML[J]. N Engl J Med, 2013,368(19):1781-1790. [10] Szuber N, Tefferi A. Chronic neutrophilic leukemia: new science and new diagnostic criteria[J]. Blood Cancer J, 2018,8(2):19.
[11] Fleischman AG, Maxson JE, Luty SB, et al. The CSF3R T618I mutation causes a lethal neutrophilic neoplasia in mice that is responsive to therapeutic JAK inhibition[J]. Blood, 2013,122(22):3628-3631.
[12] Hou XL, Ji CD, Tang J, et al. FPR2 promotes invasion and metastasis of gastric cancer cells and predicts the prognosis of patients[J]. Sci Rep, 2017,7(1):3153.
[13] Wang Z, Liu H, Shen Z, et al. The prognostic value of CXC-chemokine receptor 2
(CXCR2) in gastric cancer patients[J]. BMC Cancer, 2015(15):766.
[14] Wang WQ, Liu L, Xu HX, et al. Infiltrating immune cells and gene mutations in pancreatic ductal adenocarcinoma[J]. Br J Surg, 2016,103(9):1189-1199.
[15] Bayne LJ, Beatty GL, Jhala N, et al. Tumor-derived granulocyte-macrophage colony-stimulating factor regulates myeloid inflammation and T cell immunity in pancreatic cancer[J]. Cancer Cell, 2012,21(6):822-835.
[16] Abdol Razak N, Elaskalani O, Metharom P. Pancreatic cancer-induced neutrophil extracellular traps: a potential contributor to cancer-associated thrombosis[J]. Int J Mol Sci, 2017,18(3):E487.
[17] Geller LT, Barzily-Rokni M, Danino T, et al. Potential role of intratumor bacteria in mediating tumor resistance to the chemotherapeutic drug gemcitabine[J]. Science, 2017,357(6356):1156-1160.
[18] Szmida E, Karpinski P, Leszczynski P, et al. Aberrant methylation of ERBB pathway genes in sporadic colorectal cancer[J]. J Appl Genet, 2015,56(2):185-192.。