王颖，赵玉哲，韩青青，李思佳，张静

1河北医科大学研究生学院，石家庄 050000；2河北省人民医院腺体外科

甲状腺癌是内分泌系统中常见的恶性肿瘤，近年来发病人数持续增长[1]。2019年国家癌症中心的统计数据显示，甲状腺癌在全人群中发病率已位居第七位，甲状腺癌逐渐引起关注，但其发病机制还不明确[2]。甲状腺癌中乳头状癌最为常见，高分辨率超声结合细针穿刺活检一度被认为是术前诊断甲状腺癌的黄金搭档，但部分患者的甲状腺肿瘤相对较小，从而影响肿瘤性质的精准判断。

基因芯片又称DNA微阵列，是一种生物芯片，已成为大规模高效获取癌症基因表达谱信息的重要手段。微阵列分析是研究肿瘤基因、寻找肿瘤药物治疗的分子靶点、监测预后的一种新方法。然而，由于实验样本的异质性，使用不同的检测平台和数据处理方法会导致结果不一致。基于此，本研究以基因芯片研究为基础，应用整合的生物信息学技术进行生物信息数据挖掘，对PTC进行全面系统的生物信息分析，以期探索与PTC发生发展相关的关键基因。

1 资料与方法

1.1 临床资料 从GEO数据库(Https://www.ncbi.nlm.nih.gov/geo/)下载GSE3678、GSE60542、GSE6004和GSE53157数据集。这4个数据集均属于GPL 570平台，共包含53个组织样本，其中包括24个正常甲状腺组织样本和29个甲状腺乳头状癌样本。

1.2 差异基因的筛选 使用R语言Limma R软件包筛选每个数据集的DEGs。以|logFC|≥1和P<0.05为差异有统计学意义，并将各数据集中的基因按logFC排序，使用RoburRankAggreg(RRA)R包(Https://cran.rstudio.com/bin/windows/contrib/3.5/RobustRankAggreg_1.1.zip)对4组基因文件进行集成，并将上调和下调的DEGs列表保存起来，供后续分析。

1.3 功能注释和通路富集分析 应用David 6.8数据库(Https://david.ncifcrf.gov/)对整合的DEGs进行GO注释和KEGG途径富集分析，以了解基因参与的生物过程和信号传导途径。……