Establishment and Analysis of MicroRNA Prognostic Risk Model of Papillary Thyroid Carcinoma Based on TCGA Database
-
摘要:目的
构建基于微小RNA(miRNA)表达的预测乳头状甲状腺癌(papillary thyroid carcinoma, PTC)患者预后的生存模型。
方法从TCGA数据库官方网站上下载PTC miRNA测序数据和患者的临床资料,利用R3.6.0软件中的edgeR包筛选表达失调的miRNA。利用单因素Cox及Lasso回归分析筛选出与患者预后相关的miRNA(P < 0.05),进一步使用多因素Cox回归分析建立预后模型的风险评分方程risk score,构建生存预后模型,使用受试者工作特征曲线(ROC)来评价模型的敏感度和特异性。
结果与正常甲状腺组织相比,PTC组织中失调表达的miRNA共有75个(|log foldchange|≥2, FDR < 0.05),多因素Cox回归分析最终得到基于8个miRNA(hsa-mir-6730、hsa-mir-4709、hsa-mir-196a-2、hsa-mir-146b、hsa-mir-6860、hsa-mir-509-3、hsa-mir-513c、hsa-mir-515-1)的预测患者预后的风险模型。ROC曲线下面积(AUC)分析显示,该模型具有较好的敏感度和特异性(AUC > 0.8)。
结论成功构建了基于miRNA表达的风险预测模型,该模型可有效预测PTC患者的预后。
Abstract:ObjectiveTo construct a risk model for predicting the prognosis of papillary thyroid carcinoma (PTC) patients based on the expression of microRNAs (miRNA).
MethodsThe original microRNAs sequencing data and corresponding clinical data of PTC patients were downloaded from the official website of TCGA database. The aberrant expression of microRNAs in PTC was screened by edgeR package in R 3.6.0 software. Univariate and multivariate Cox regression analyses as well as the Lasso regression analysis were utilized to screen out the microRNAs which were related to the prognosis of patients. Then, the risk score equation was established by multivariate Cox regression analysis and the risk prognosis model was constructed. The receiver operating characteristic (ROC) curve was utilized to evaluate the sensitivity and specificity of the model.
ResultsCompared with normal tissues, 75 microRNAs (|log foldchange|≥2, FDR < 0.05) were dysregulated in PTC tissues. Multivariate Cox regression analysis revealed eight microRNAs (hsa-mir-6730, hsa-mir-4709, hsa-mir-196a-2, hsa-mir-146b, hsa-mir-6860, hsa-mir-509-3, hsa-mir-513c, hsa-mir-515-1) which could be used to construct risk model for predicting the prognosis of patients. The model had good sensitivity and specificity (AUC > 0.8).
ConclusionA risk prediction model based on the expression of microRNAs has been successfully constructed, and the model could effectively predict the prognosis of thyroid cancer patients.
-
Key words:
- Thyroid cancer /
- MicroRNAs /
- TCGA /
- Prognostic model
-
0 引言
鼻咽癌是最常见的头颈部肿瘤之一。我国为鼻咽癌高发地区,每年的发病率约为20/10万[1],由于鼻咽解剖结构及生物学行为的特殊性,很难行手术治疗,目前鼻咽癌公认和有效的治疗手段为放射治疗或以放疗为主的综合治疗。虽然放疗技术不断进步与放疗设备的不断更新,鼻咽癌的生存率有了较大的提高,但5年生存率仍为60%~80%[2],部分患者仍未能获得长期生存。TNM分期系统是鼻咽癌预后判断和指导治疗的重要依据,但临床发现同一分期患者即使接受相同的治疗方案,预后却不同[3-4],这提示鼻咽癌生物学差异的存在,仅基于解剖学信息的TNM临床分期系统还不能准确地预测鼻咽癌患者的预后。虽然EB病毒滴度、表皮生长因子受体、microRNA也可提示鼻咽癌的预后[5-7],但检测成本高,需要多中心合作,临床上可行性差。所以,亟需检测方便、价格低廉可预测鼻咽癌预后的标志物。
流行病学研究证实,约25%的肿瘤由炎性反应发展而来,其与肿瘤的发生发展密切相关并且影响肿瘤患者的预后[8]。炎性反应指标,如白细胞计数[9]、血小板计数[10-11]、中性粒淋巴细胞比(neutrophil-lymphocyte ratio, NLR)[12-13]、血小板淋巴细胞比(platelet-lymphocyte ratio, PLR)[14-15]被发现可作为肿瘤的独立预后因素。这些血液指标检测方便,价格低廉,可广泛应用于临床,评估患者预后。本研究通过对91例鼻咽癌患者临床资料进行回顾性分析,评价治疗前PLR和NLR与鼻咽癌患者预后的相关性,为评估预后提供客观依据。
1 资料与方法
1.1 临床资料
回顾性收集2009年1月至2013年9月期间于西安交通大学第一附属医院和陕西省人民医院初治的91例鼻咽癌患者,所有病例均经病理证实。临床资料完整。排除标准:(1)合并有免疫性疾病以及其他恶性肿瘤的患者;(2)治疗前合并有急性或慢性感染;(3)合并有血液系统疾病、血栓或出血性疾病;(4)合并有严重的肝、肾疾病;(5)治疗前曾接受过放疗或化疗;(6)无远处转移。记录患者治疗前的中性粒细胞计数、淋巴细胞计数及血小板计数结果。
1.2 治疗及随访方法
入选患者采用3D-CRT或IMRT根治性放疗(有或无化疗),Ⅰ期患者仅接受单纯放射治疗,Ⅱ、Ⅲ、Ⅳ期患者接受以顺铂和5-氟尿嘧啶为主的辅助或同步放化疗。鼻咽原发灶和颈部转移淋巴结剂量为(70~76)Gy/(7~8)w/(35~38)f,颈部预防区域剂量为(50~60)Gy/(5~6)w/(25~30)f。根据患者的临床分期及不良反应给予2~6周期的全身化疗,化疗方案为:顺铂25 mg/m2,第1~3天静脉滴注;5-氟尿嘧啶500 mg/m2,第1~5天静脉滴注,每21天重复1周期。患者治疗结束后均定期随访,治疗后前2年,每3月检查一次,2年后半年复查一次,5年后1年复查1次。随访截止时间为2016年9月。
1.3 统计学方法
采用SPSS19.0软件对数据进行统计学分析。绘制ROC曲线确定PLR和NLR与总生存期(overall survival, OS)及无进展生存期(progression-free survival, PFS)的相关性,选取截断值。应用Kaplan-Meier法进行生存分析并采用Log rank检验来检验。采用Cox比例风险回归模型分析多种因素对生存时间的影响。以P < 0.05为差异有统计学意义。
2 结果
2.1 鼻咽癌患者临床病理资料
91例患者的基本特征资料见表 1,中位年龄53岁(12~72)岁,女30例,男61例,男女比例2:1,Ⅰ、Ⅱ、Ⅲ、Ⅳ期患者分别为2、27、42、20例。单纯放疗患者9例,82例患者接受辅助或同步放化疗,所有患者均按期完成放化疗。中位随访时间为44月(6~87)月,其中44例出现复发或转移,39例患者死亡。患者的1、3、5年总生存率分别为92.3%、72.1%、56.8%,1、3、5年无进展生存率分别为82.4%、60.9%、53.3%。
表 1 91例鼻咽癌患者临床基本特征资料(n(%))Table 1 Basic clinical features of 91 nasopharyngeal carcinoma patients (n(%))2.2 ROC曲线选取PLR和NLR预后相关截断值
以OS作为终点,PLR、NLR为检测变量,绘制ROC曲线选取截断值分别为143.3、2.6,两者的曲线下面积分别为0.640、0.739,见图 1。
以PFS作为终点,PLR、NLR为检测变量,绘制ROC曲线选取截断值分别为143.3、2.6,两者的曲线下面积分别为0.657、0.694,见图 2。说明治疗前PLR、NLR与患者的预后存在相关性,利用ROC曲线选取的截断值进行进一步生存分析。
2.3 Kaplan-Meier生存分析、Cox单因素和多因素分析
PLR≥143.3组和PLR < 143.3组患者生存曲线比较,差异有统计学意义(P=0.022),见图 3~4。NLR≥2.6组和NLR < 2.6组患者生存曲线比较,差异有统计学意义(P=0.044),见图 5~6。
Cox单因素分析显示除性别、年龄以外,TNM分期、治疗前PLR≥143.3、NLR≥2.6均是影响鼻咽癌患者OS和PFS的不良预后因素(P < 0.05),见表 2。Cox多因素分析显示治疗前PLR≥143.3(RR=2.491, 95%CI=1.139~5.451, P=0.022)、NLR≥2.6(RR=2.186, 95%CI=1.021~4.682,P=0.044)是鼻咽癌患者OS的独立危险因素,而治疗前PLR≥143.3(RR=2.461,95%CI=1.242~4.874, P=0.010)是鼻咽癌患者PFS的独立危险因素,见表 3。
表 2 影响鼻咽癌患者生存预后的Cox单因素分析Table 2 Cox univariate analysis of prognostic factors for nasopharyngeal carcinoma patients表 3 影响鼻咽癌患者生存预后的Cox多因素分析Table 3 Cox multivariate analysis of prognostic factors for nasopharyngeal carcinoma patients3 讨论
鼻咽癌对放射线高度敏感,因此放疗成为主要治疗手段。随着三维适形放疗和调强放射治疗的临床应用,鼻咽癌的生存率较前明显提高,但5年生存率仍仅为60%~80%。多项研究表明鼻咽癌患者预后与众多因素有关,包括患者年龄、临床分期、EB病毒感染及贫血等。此外,肿瘤的预后还与机体本身的炎性反应有关。炎性反应包含中性粒细胞、淋巴细胞、血小板、C反应蛋白等多种指标,其中PLR、NLR已受到越来越多专家的关注。本研究发现治疗前PLR和NLR可能成为鼻咽癌的独立预后因素。
恶性肿瘤患者常伴随血小板的升高,实验研究表明血小板参与肿瘤细胞生长、转移及血管生成[16]。临床研究表明血小板数目升高与肿瘤患者较差预后相关[11, 17]。此外研究表明中性粒细胞可促使机体产生多种促肿瘤生长因子和蛋白酶,促进肿瘤的发生、发展[18]。而淋巴细胞参与机体的免疫反应是抗肿瘤免疫的重要组成部分,淋巴细胞减少说明机体免疫机制异常,抗肿瘤免疫力下降,为肿瘤生长、浸润和转移提供条件。随着肿瘤进展,机体内炎性反应与肿瘤失去平衡,体内淋巴细胞降低,而血小板、中性粒细胞升高,相应的PLR和NLR比值也增高,机体内促进肿瘤炎性反应与抗肿瘤炎性反应的平衡状态被打破。因此PLR和NLR是反应机体免疫情况的重要指标,两者的升高能促进肿瘤进展,导致肿瘤患者预后不良。既往研究结果显示高PLR和NLR可影响宫颈癌、乳腺癌、结直肠癌等恶性肿瘤的预后[19-21]。而目前关于PLR、NLR与鼻咽癌患者预后相关性的研究较少,Sun等[21]分析了251例鼻咽癌患者治疗前PLR和NLR,结果证明治疗前两者水平是影响鼻咽癌患者生存独立预后因素。本研究结果显示治疗前PLR、NLR与鼻咽癌患者的总生存期和无进展生存期具有相关性。Cox多因素分析提示PLR≥143.3、NLR≥2.6和TNM分期是影响鼻咽癌患者治疗后的独立危险因素。PLR≥143.3组患者有较短OS和PFS,而NLR≥2.6组患者有较差的OS,和本研究结果相一致。因此,高PLR、NLR的鼻咽癌患者总生存率要低于低PLR、NLR的患者,且高PLR的患者复发或转移风险明显增加。据此,临床上或许可以通过提高鼻咽癌患者免疫功能及降低机体炎性反应,改善患者的预后。
但由于本研究是一个相对小样本的回顾性研究,不能代表大部分的鼻咽癌患者,且随访时间较短,存在一定的局限性,因此需要进行多中心、大样本的前瞻性研究来进一步证实。
本研究结果表明,治疗前PLR和NLR水平与鼻咽癌患者预后具有相关性,可能是影响鼻咽癌患者预后的独立危险因素,NLR和PLR的获取具有简便、经济的优点,可以作为鼻咽癌患者病情评估的一个有益补充,值得推广。目前鼻咽癌相关有效预后指标较多,笔者将在今后的临床研究工作中继续探索,将本研究指标与已有的有效预后指标进行比较,从而提高治疗前PLR和NLR水平这一预后指标应用于临床的合理性及可靠性。
作者贡献焦赵爽:研究设计、实施、分析数据及论文撰写张怀念:数据整理、分析数据及修改论文 -
表 1 PTC组织中前十位差异表达的miRNA
Table 1 Top 10 differentially-expressed miRNA in papillary thyroid carcinoma (PTC) tissues
-
[1] 董芬, 张彪, 单广良.中国甲状腺癌的流行现状和影响因素[J].中国癌症杂志, 2016, 26(1): 47-52. http://d.old.wanfangdata.com.cn/Periodical/zgazzz201601008 Dong F, Zhang B, Shang GL. Distribution and risk factors of thyroid cancer in China[J]. Zhongguo Ai Zheng Za Zhi, 2016, 26(1): 47-52. http://d.old.wanfangdata.com.cn/Periodical/zgazzz201601008
[2] Raue F, Frank-Raue K. Thyroid Cancer: Risk-Stratified Management and Individualized Therapy[J]. Clinical Cancer Res, 2016, 22(20): 5012-5021. doi: 10.1158/1078-0432.CCR-16-0484
[3] Glaser SM, Mandish SF, Gill BS, et al. Anaplastic thyroid cancer: Prognostic factors, patterns of care, and overall survival[J]. Head Neck, 2016, 38(Suppl 1): E2083-E2090. https://www.ncbi.nlm.nih.gov/pubmed/26894506
[4] 林琳, 刘建华, 吴丽娜.甲状腺癌RNA分子研究进展[J].肿瘤防治研究, 2019, 46(11): 1031-1035. doi: 10.3971/j.issn.1000-8578.2019.19.0431 Lin L, Liu Jh, Wu LN. Research Progress of RNA Molecules in Thyroid Cancer[J]. Zhong Liu Fang Zhi Yan Jiu, 2019, 46(11): 1031-1035. doi: 10.3971/j.issn.1000-8578.2019.19.0431
[5] Li N, Liu X, Han L, et al. Expression of miRNA-146b-5p in patients with thyroid cancer in combination with Hashimoto's disease and its clinical significance[J]. Oncol Lett, 2019, 17(6): 4871-4876. doi: 10.3892/ol.2019.10173/download
[6] 刘延晴, 林岩松.碘难治性分化型甲状腺癌诊治策略及预后[J].中国实用外科杂志, 2019, 39(3): 216-220. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgsywkzz201903007 Liu YQ, Lin YS. Emerging management and its impact on radioiodine refractory differentiated thyroid cancer[J]. Zhongguo Shi Yong Wai Ke Za Zhi, 2019, 39(3): 216-220. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgsywkzz201903007
[7] 杨进宝, 李小毅, 商中华.甲状腺乳头状癌各亚型临床病理特点的研究进展[J].癌症进展, 2015, 13(1): 55-60. http://d.old.wanfangdata.com.cn/Periodical/azjz201501012 Yang JB, Li XY, Shang ZH. Advances in clinicopathological characteristics of various subtypes of papillary thyroid cancer[J]. Ai Zheng Jin Zhan, 2015, 13(1): 55-60. http://d.old.wanfangdata.com.cn/Periodical/azjz201501012
[8] Abdullah MI, Junit SM, Ng KL, et al. Papillary Thyroid Cancer: Genetic Alterations and Molecular Biomarker Investigations[J]. Int J Med Sci, 2019, 16(3): 450-460. doi: 10.7150/ijms.29935
[9] Raposo L, Morais S, Oliveira MJ, et al. Trends in thyroid cancer incidence and mortality in Portugal[J]. Eur J Cancer Prev, 2017, 26(2): 135-143. doi: 10.1097/CEJ.0000000000000229
[10] Bray F, Ferlay J, Laversanne M, et al. Cancer Incidence in Five Continents: Inclusion criteria, highlights from Volume X and the global status of cancer registration[J]. Int J Cancer, 2015, 137(9): 2060-2071. doi: 10.1002/ijc.29670
[11] Grogan RH, Kaplan SP, Cao H, et al. A study of recurrence and death from papillary thyroid cancer with 27 years of median follow-up[J]. Surgery, 2013, 154(6): 1436-1446, discussion 1446-1437. doi: 10.1016/j.surg.2013.07.008
[12] 马珊, 尹乐平.甲状腺乳头状癌手术方式的选择[J].实用医院临床杂志, 2017, 14(3): 124-127. doi: 10.3969/j.issn.1672-6170.2017.03.044 Ma S, Yin LP. Choice of surgical method of thyroid papillary carcinoma[J]. Yi Yong Yi Yuan Lin Chuang Za Zhi, 2017, 14(3): 124-127. doi: 10.3969/j.issn.1672-6170.2017.03.044
[13] Trzybulska D, Bobjer J, Giwercman A, et al. Serum microRNAs in male subfertility-biomarkers and a potential pathogenetic link to metabolic syndrome[J]. J Assist Reprod Genet, 2017, 34(10): 1277-1282. doi: 10.1007/s10815-017-0989-0
[14] Hill M, Tran N. MicroRNAs Regulating MicroRNAs in Cancer[J]. Trends Cancer, 2018, 4(7): 465-468. doi: 10.1016/j.trecan.2018.05.002
[15] Yu L, Xiang L, Feng J, et al. miRNA-21 and miRNA-223 expression signature as a predictor for lymph node metastasis, distant metastasis and survival in kidney renal clear cell carcinoma[J]. J Cancer, 2018, 9(20): 3651-3659. doi: 10.7150/jca.27117
[16] Liu Q, Diao R, Feng G, et al. Risk score based on three mRNA expression predicts the survival of bladder cancer[J]. Oncotarget, 2017, 8(37): 61583-61591. doi: 10.18632/oncotarget.18642
[17] Butz H, Patocs A. MicroRNAs in endocrine tumors[J]. EJIFCC, 2019, 30(2): 146-164. https://www.ncbi.nlm.nih.gov/pubmed/31263390
[18] Yue K, Wang X, Wu Y, et al. microRNA-7 regulates cell growth, migration and invasion via direct targeting of PAK1 in thyroid cancer[J]. Mol Med Rep, 2016, 14(3): 2127-2134. doi: 10.3892/mmr.2016.5477
[19] Wang J, Yang H, Si Y, et al. Iodine Promotes Tumorigenesis of Thyroid Cancer by Suppressing Mir-422a and Up-Regulating MAPK1[J]. Cell Physiol Biochem, 2017, 43(4): 1325-1336. doi: 10.1159/000481844
[20] Dong S, Meng X, Xue S, et al. microRNA-141 inhibits thyroid cancer cell growth and metastasis by targeting insulin receptor substrate 2[J]. Am J Transl Res, 2016, 8(3): 1471-1481. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4859632/
[21] Wang P, Meng X, Huang Y, et al. MicroRNA-497 inhibits thyroid cancer tumor growth and invasion by suppressing BDNF[J]. Oncotarget, 2017, 8(2): 2825-2834. doi: 10.18632/oncotarget.13747
[22] Chou CK, Liu RT, Kang HY. MicroRNA-146b: A Novel Biomarker and Therapeutic Target for Human Papillary Thyroid Cancer[J]. Int J Mol Sci, 2017, 18(3). pii: E636. doi: 10.3390/ijms18030636
[23] Qiu Z, Li H, Wang J, et al. miR-146a and miR-146b in the diagnosis and prognosis of papillary thyroid carcinoma[J]. Oncol Rep, 2017, 38(5): 2735-2740. doi: 10.3892/or.2017.5994
[24] Fu YT, Zhang DQ, Zhou L, et al. Has-MiR-196a-2 is up-regulated and acts as an independent unfavorable prognostic factor in thyroid carcinoma[J]. Eur Rev Med Pharmacol Sci, 2018, 22(9): 2707-2714. https://pubmed.ncbi.nlm.nih.gov/29771422/
[25] Zhang G, Liu Z, Han Y, et al. Overexpression of miR-509 Increases Apoptosis and Inhibits Invasion via Suppression of Tumor Necrosis Factor-alpha in Triple-Negative Breast Cancer Hs578T Cells[J]. Oncol Res, 2016, 24(4): 233-238. doi: 10.3727/096504016X14648701447977
[26] Tang J, Kong D, Cui Q, et al. Bioinformatic analysis and identification of potential prognostic microRNAs and mRNAs in thyroid cancer[J]. Peer J, 2018, 6: e4674. doi: 10.7717/peerj.4674
[27] Xia HL, Lv Y, Xu CW, et al. MiR-513c suppresses neuroblastoma cell migration, invasion, and proliferation through direct targeting glutaminase (GLS)[J]. Cancer Biomark, 2017, 20(4): 589-596. doi: 10.3233/CBM-170577
[28] Zhang K, Zhao Z, Yu J, et al. LncRNA FLVCR1-AS1 acts as miR-513c sponge to modulate cancer cell proliferation, migration, and invasion in hepatocellular carcinoma[J]. J Cell Biochem, 2018, 119(7): 6045-6056. doi: 10.1002/jcb.26802
[29] Wei HT, Guo EN, Liao XW, et al. Genomescale analysis to identify potential prognostic microRNA biomarkers for predicting overall survival in patients with colon adenocarcinoma[J]. Oncol Rep, 2018, 40(4): 1947-1958. https://www.spandidos-publications.com/or/40/4/1947
[30] Liao X, Wang X, Huang K, et al. Genome-scale analysis to identify prognostic microRNA biomarkers in patients with early stage pancreatic ductal adenocarcinoma after pancreaticoduodenectomy[J]. Cancer Manag Res, 2018, 10: 2537-2551. doi: 10.2147/CMAR.S168351
[31] 贾晓晨, 贾勇圣, 孟文静, 等.基于TCGA数据库建立的八基因预后模型在乳腺癌中的应用[J].天津医药, 2018, 46(8): 856-861. http://d.old.wanfangdata.com.cn/Periodical/tjyy201808017 Jia XC, Jia YS, Meng WJ, et al. Identification of prognostic eight-gene signature model in breast cancer using integrated TCGA database[J]. Tianjin Yi Yao, 2018, 46(8): 856-861. http://d.old.wanfangdata.com.cn/Periodical/tjyy201808017