1、论著 收稿日期 基金项目 邯 郸 市 科 学 技 术 研 究 与 发 展 计 划 项 目 任 务（）作者简介 梁新科（），男，河北磁县人，河北省邯郸市邯钢医院主治医师，医学学士，从事甲状腺疾病诊治研究。通信作者。：基于生物信息学分析 在甲状腺乳头状癌中的表达及相关分析梁新科，王鹏飞，霍胜男，李建英，侯庆田（河北省邯郸市邯钢医院甲状腺科，河北 邯郸 ）摘要目的 通过生物信息学分析探讨甲状腺乳头状癌（，）的关键基因及其潜在调控分子靶标，探索其在 中的分子调控机制。方法 通过 数据库检索并下载甲状腺乳头状癌的基因表达谱，使用 中的“”软件包寻找甲状腺乳头状癌与癌旁组织之间的是否存在差异长链非编码（，
2、）。使用 构建及临床相关性分析。通过 进行富集分析和注释。结果 经 数据库筛选，共 例甲状腺乳头状癌的肿瘤组织与 例癌旁组织被纳入分析，共得到了 个差异表达 ，其中包括 个高表达基因与 个低表达的基因。使用动态树切割法识别基因模块，设定模块内最低基因个数为 ，最终得到 个相应的模块，结果显示，蓝色模块与甲状腺乳头状癌显著相关。蓝色模块内差异 共 有 个，对 模 块 内 个 差 异 进 行 单 因 素 和 多 因 素 回 归 分 析，得 到 。随后对 绘制受试者工作特征曲线，曲线下面积为 。通过和 分析，揭示 的潜在通路与功能，结果显示 可能通过脂肪细胞因子信号途径、信号通路、钙黏蛋白结合参与细
3、胞黏附、对凋亡信号通路的正调控等功能参与甲状腺乳头状癌的病理过程。结论 在甲状腺乳头状癌患者与健康人群的甲状腺组织的表达差异有统计学意义，基因有望成为甲状腺乳头状癌的生物标志基因，为后续的 的研究提供重要依据。关键词甲状腺癌，乳头状；生物信息学；：中图分类号 文献标志码 文章编号 （），（，）（），（），第 卷第期 年月河北医科大学学报 ，（），（），；甲状腺乳头状癌（，）是最常见的甲状腺癌类型，占甲状腺恶性肿瘤的 ，并且全球发病率逐年增加。虽然甲状腺乳头状癌普遍预后良好，但其肿瘤的复发和转移亦严重影响患者的生存率。近年来，基因检测得以广泛普及，可以研究基因表达与疾病之间的关系，为肿瘤机制的研
4、究提供了明确的方向。本研究通过生物信息学方法筛选分析 与健康人群的甲状腺组织的差异表达基因，寻找其发挥作用的相关信号通路，探索 的发病机制，为 诊断及治疗提供新的研究方向。资 料 与 方 法一般资料选取 数据库中 年月 年 月收录并具有完整随访信息的甲状腺乳头状癌样本 例与癌旁组织样本 例。纳入标准：肿瘤组织样本来自甲状腺乳头状癌患者，对照组为正常组织；样本须同时含有基因表达谱数据与临床信息。数 据 的 分 析使 用 中 的“”软件包对 数据进行差异分析，寻找甲状腺乳头状癌组织与正常组织之间是否存在差异长 链 非 编 码（，），并绘制火山图和热图。差异 筛选标准为：，。构 建 及 临 床 相
5、关 性 分 析使 用“”包将该数据中的探针 转换成基因符号以进行进一步分析。使用软件 中 的“”软 件 包 进 行 样 本 层 次 聚 类（），过 滤 变 异 最 大 的 基 因。利 用 软件包中的“”函数调节参数的权重。使用 软件包将数据中有相关性和相邻关系的 计算成为拓扑重叠矩阵（，），然后计算相应的相异度（）。使用作为距离度量，使用动态树切割法进行分层聚类以识别模块，将高度相似的模块通过聚类标记并合并。最小基因模块 ，合并相似模块的阈值用于搜索在心肌甲状腺乳头状癌中起重要作用的模块。获得基因模块后根据模块化特征将甲状腺乳头状癌患者的临床特征和模块内的基因联合分析，探究模块内基因的生物学意
6、义。使用相关性热图探究基因与模块的相关性，并进行聚类分析。此外，算法还对在同一模块中基因进行相互作用网络的预测。是一种分析软件，可以为生物学家提供生物网络分析和二维可视化。模块中差异 的筛选 是一种生物学分析软件。本研究中，使用维恩图，将 模块中的 与差 异 取 交 集，得 到 模块中差异 。统 计 学 方 法本 研 究 在 中 使 用“”包进行单因素和多因素 回归分析，进行模 型 变 量 筛 选。受 试 者 工 作 特 征（，）曲线用于评估基因的诊断价值。为差异有统计学意义。功能富集分析选取甲状腺乳头状癌样本，将其分为 低表达组和 高表河 北 医 科 大 学 学 报第 卷第期达组，随后导入
7、中进行基因本体论（，）分析和京都基因与基因组百科全书（，）分析。结果差异 的筛选经 数据库筛选，共获得差异表达 个，其中包括高表达基因 个与低表达基因 个（），见图。表达情况见图。加权基因共表达网络构建及临床相关性分析“”函数的结果显示，当权重参数 时，（）与 （）之间相关系数的平方。选取软阈值 构建 ，见图。使用动态树切割法识别基因模块，设定模块内最低基因个数为 ，最终得到 个相应的模块，见图。结果显示，模块 与甲状腺乳头状癌显著相关，见图。基因与模块的相关性分别展现在图 与图 中。预测出的 模块中基因之间的相互作用见图。图甲状腺乳头状癌基因表达差异甲状腺乳头状癌基因表达差异火山图；甲状腺乳
8、头状癌 表达情况 图加权基因共表达网络构建构建 ；基因模块；蓝色模块与甲状腺乳头状癌显著相关；基因与模块的相关性；基因与模块的相关性；蓝色模块中基因之间的相互作用 模块中差异基因的筛选通过维恩图将得到甲状腺乳头状癌组织与正常组织之间的 个差异河 北 医 科 大 学 学 报第 卷第期 与 模块中 个模块 取交集，得到蓝色模块内差异 共有 个，见图。图蓝色模块中的差异 预后模型的构建与验证对蓝色模块内 个差异 进行单因素和多因素 回归分析（图、），得到 。随后对 绘制 曲 线，曲 线 下 面 积（，）值为 （图）。基因富集分析通过 和 分析，揭示 的 潜 在 通 路 与 功 能，结 果 显 示 通
9、过脂肪细胞因子信号途径、丝裂原活化 蛋 白 激 酶（，）信号通路、钙黏蛋白结合参与细胞黏附、对凋亡信号通路的正调控等功能参与甲状腺乳头状癌的病理过程，见图。图差异 的回归分析单因素 回归分析；多因素 回归分析；曲线 图富集分析显示 的潜在通路与功能 低表达组的 分析；高表达组的 分析；低表达组的 分析；高表达组的 分析 讨论甲状腺癌是最常见的内分泌肿瘤，占人类所有恶性肿瘤的。被认为是复发率最高的甲状腺癌类型，但同时也被认为是具有良好预后 的 惰 性 癌 症，占 已 分 化 甲 状 腺 癌 的 。对于大多数 患者，预后良好，然而，在某些情况下会出现侵袭性行为，的患者表现 出 临 床 或 隐 匿
10、性 颈 部 淋 巴 结 转 移。已被证实参与 的细胞迁移和增殖过程。因 此，寻 找 肿 瘤 与 正 常 组 织 之 间 的 差 异 有助于进一步了解甲状腺乳头状癌的发病机制，为 的术前诊断提供生物标志物和治疗靶点。基于 生 物 信 息 学 分 析，研 究 显 示，在甲状腺乳头状癌中表达与癌旁组织差异有统计学意义（），同时 和 结果显示 可能通过脂肪细胞因子信号河 北 医 科 大 学 学 报第 卷第期途径、信号通路、钙黏蛋白结合参与细胞黏附、对凋亡信号通路的正调控等功能参与甲状腺乳头状癌的病理过程。近年来研究显示，肥胖患病率的增加与包括乳腺癌、结肠癌和前列腺癌在内的各种人类癌症的发病率、进展和预
11、后较差有关。最近，一些研究显示甲状腺癌的发病率随着肥胖率的显著上升而增加，并且越来越多的证据表明肥胖与甲状腺癌风险增加之间存在关联。脂肪因子是脂肪细胞分泌的各种酶、激素、细胞因子、生长因子、蛋白质和其他生物活性物质，包括脂联素、瘦素、抵抗素和白细胞介素等。据报道，体重指数越高的患者脂肪因子浓度越高。由于营养过剩或运动不足所致肥胖时，脂肪因子将出现异常分泌，并参与多种病态的发生。脂肪因子参与以下病理和生理过程，如胰岛素敏感性、炎症和增殖，这些在肿瘤发生和发展过程中均发挥重要作用。因此，脂肪因子可能是将肥胖与甲状腺癌联系起来的目标之一。本研究显示，可能通过脂肪因子参与了甲状腺乳头状癌的发生发展过程
12、。通路代表一种普遍存在的信号转导通路，在细胞活力、凋亡和侵袭的调控中发挥重要作用。信号通路能够被受体酪氨酸激酶、蛋白偶联受体和部分细胞因子受体激活，将细胞外的信号刺激传递给细胞，引起细胞增殖、凋亡、分化等一系列反应，进而促进肿瘤的发生发展。先前的研究表明，信号通路在多种癌症中被激活，包括胃癌、肺癌、卵巢癌 和肝癌。这一途径已被证明在肿瘤发生发展中发挥关键作用，因此，成为包括甲状腺癌在内的恶性肿瘤的重要治疗靶点。属 于 丝 氨 酸苏 氨 酸 蛋 白 激 酶，信号级联是肿瘤细胞增殖、分化、凋亡和耐药性的关键途径。此外，激活可以磷酸化核转录因子和蛋白激酶，调节相关基因的转录并参与各种生理过程。信号通
13、路由个子通路组成：（）、（）、和 （）。研究表明，信号通路的激活促进细胞生长、分化、迁移和存活。和 的激活涉及细胞分化、凋亡和细胞存活的变化。研究显示，通路中，的突变及 重排能够促进甲状腺滤泡细胞向乳头状甲状腺癌的转化，这种改变被认为是甲状腺癌发生和进展的早期表现。据 等 报道，根据甲状腺癌的遗传学变异，有 的患者伴随 基因突变。其中 突变占甲状腺乳头状癌的 。等 研究显示 基 因 激 活 的 非 编 码（，）通过激活 信号通路进而对乳头状甲状腺癌细胞增殖和迁移产生影响。还 可 通 过 调 节 上 皮 间 质 转 化（）信号通路发挥作用，上皮间质转化被认为是癌细胞从原发部位向远处组织或淋巴结转
14、移过程中的一个重要步骤，通过该过程上皮细胞逐渐获得间充质特征并增强迁移和侵袭的特性。甲状腺乳头状癌的多灶性病变显示出淋巴结转移倾向，的淋巴结转移导致预后不良的结果，研究显示参与了 细胞的局部转移过程。的分子特征包括将极化上皮表型标志物（如 ）转换为间充质标志物（如 、纤连蛋白等）。转录因子（如 、和 等）的 激 活 触 发的发生。据报道，的激活会增加 、和 钙黏蛋白的水平，并降低 钙黏蛋白的表达水平，从而导致 。通路已被广泛证明在肿瘤进展和转移中起关键作用，可作为甲 状 腺 癌 的 临 床 治 疗 靶 点，。本 研 究 显 示 可能通过 通路，改变细胞黏附状态，形成环境，促进 的淋巴结转移，导
15、致 的进展。综上所述，本研究基于 甲状腺乳头状癌及癌旁组织的基因数据，显示 在甲状腺乳头状癌和癌旁组织中的表达差异有统计学方法意义，而 且 通 过和 分 析 其 可 能 通 过 通路、脂肪因子通路，调节 及凋亡的发生，从而参与甲状腺乳头状癌的发生发展过程，为验证 的作用提供理论依据。参考文献 ，（），（）：，（），（）：河 北 医 科 大 学 学 报第 卷第期 ，（）：，（）：，（）：，：，：（）：，（）：，（）：，（）：，（）：，（）：，（）：，（）：，（）：，（）：，：，（）：，（）：，（）：，（）：，（）：，（）：，：，（）：，（）：，（）：，：，（）：（本文编辑：何祯）梁新科等基于生物信息学分析 在甲状腺乳头状癌中的表达及相关分析
