姜雪丽，赵 辛，赵铁铮

甲状腺癌是一种好发于中老年女性的内分泌系统常见恶性肿瘤，多数病理类型为甲状腺乳头状癌[1-4]。据统计，2017年全球新发甲状腺癌病例占所有恶性肿瘤新发病例3.4%，在女性恶性肿瘤新发病例中甲状腺癌居第5位[5]。我国每年甲状腺癌新发病例和死亡病例数占全世界15%左右[6]。近年甲状腺癌患者病死率逐年增加[7-9]。目前，临床上对甲状腺癌患者预后判断主要基于临床病理特征，但由于恶性肿瘤有明显异质性，且患者个体差异较大，故此预测手段具有明显局限性[10]。近年，肿瘤基因和蛋白表达预测患者预后正逐步受到重视。

本研究通过对GEO和TCGA数据库中甲状腺癌患者转录组数据进行联合分析，筛选出甲状腺癌中差异表达的基因;随后通过单因素Cox回归分析，筛选出与预后相关的差异基因;最后通过多因素Cox回归分析成功构建5个差异基因预测甲状腺癌患者预后的模型，并结合患者临床资料对模型进行验证，用于判断甲状腺癌患者预后，现报告如下。

1 资料与方法

1.1资料收集和差异分析 在GEO数据库中以“Thyroid cancer”为关键词，选取人甲状腺癌组织样本转录组数据集GSE138198和GSE50901，使用R软件及Sva功能包校正批次，利用Limma包分析差异基因，认为符合|log2(Fold Change)|>0.585且P<0.05的基因为差异基因。在TCGA数据库中下载人甲状腺癌RNA-Seq数据，共含有58份正常甲状腺样本和509份甲状腺癌样本，利用R软件及Limma包分析差异基因，认为符合|log2(Fold Change)|>2且P<0.01的基因为差异基因。甲状腺癌患者的临床资料利用R软件进行整理，包括性别、年龄、临床分期、风险值评分、生存时间和生存状态等，去除临床资料不完整患者，最终保留498例采用单因素和多因素Cox回归分析进行甲状腺癌患者临床病理特征与预后关系的研究。

1.2Cox比例风险回归模型建立和验证 利用R软件将TCGA数据库中甲状腺癌患者随机分为两组，即试验组和验证组，使用Survival包对在GEO和TCGA数据库中均发生改变基因与患者预后的关系进行单因素Cox回归分析；使用Glmnet包对单因……