谢鋆晖，余济春

（1.南昌大学第二附属医院，江西 南昌 330000；2.赣州市肿瘤医院头颈科，江西 赣州 341000）

甲状腺癌（Thyroid carcinoma，THYC）是由甲状腺滤泡细胞发展而来的一种内分泌系统的肿瘤，甲状腺乳头状癌（Papillary thyroid carcinoma，PTC）是最常见的类型（约占80%），它是一种分化型癌［1-2］。近年来，甲状腺癌的死亡率和乳头状甲状腺癌的发病率不断增加，给人们带来很大的困扰。

随着测序技术的发展和RNA-sequencing（RNAseq）测序技术的应用，基于mRNA水平的生物分子在各种疾病的诊断、预后和治疗中得到广泛关注［3］。RNA-seq使我们对基因的转录本及翻译机制有了更加深入的了解，将RNA序列检测基因融合分析应用于已知致病转录本的差异表达等，RNA-seq已成为科研工作中不可或缺的工具［4］。同时，它也是医学生物学（包括肿瘤领域）中重要的手段［5］。SINGH等通过RNA-seq揭示了甲状腺及甲状腺激素的调节机制［6］。QIU等通过转录组测序手段推测PTC进展的部分分子机制，并预测了治疗的生物标志物［7］。

本研究中我们通过GEO数据集的分析初步筛选出癌与正常组织的差异基因，并进一步结合肿瘤基因组图谱（The cancer genome atlas，TGCA）及人蛋白数据库数据发现CDH2、SeprinA1及FN1与甲状腺癌的关系最为密切。现报告如下。

1 材料与方法

1.1 基因表达数据从GEO数据库（https://www.ncbi.nlm.nih.gov/geo/）下载GSE64912数据集，共计22例RNA转录组及表达谱分析结果，其中包含18例乳头状甲状腺癌与4例健康对照［8］。

1.2 基因数据前处理通过在线网站iDEP.92（http：//bioinformatics.sdstate.edu/idep92/）对下载的基因表达谱进行数据分析。导入表达谱的基因读数（Read Counts）数据进行数据标准化，通过DESeq2软件包实现常规对数log2转换，减少均值相关方差，输出的数据用于后续分析。……