IF31.7｜多族裔GWAS揭示甲状腺癌与良性结节遗传异质性

883个位点、570个新发现，PRS可识别侵袭性甲状腺癌高危人群。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

甲状腺疾病高度流行且遗传度较高，甲状腺癌是最常见的内分泌恶性肿瘤，家族研究估计遗传因素对甲状腺癌易感性贡献高达53%，自身免疫性甲状腺病遗传度约75%。甲状腺结节在一般人群中超声检出率可达65%，但细针穿刺活检大部分结果为良性或不确定，导致大量不必要的侵入性操作。既往GWAS已在欧洲人群中发现若干甲状腺癌相关位点，但样本量有限、多族裔覆盖不足，且对良性结节与恶性结节的遗传关联缺乏系统比较。

科学问题

卡在哪里： 此前GWAS样本量有限，主要基于欧洲人群，对非欧洲族裔的遗传结构了解不足；良性结节与甲状腺癌之间的遗传重叠和特异性未被系统解析；缺乏将遗传风险评分应用于区分良恶性结节及预测侵袭性甲状腺癌的临床验证。

本文要回答： 通过大规模多族裔GWAS元分析，全面鉴定甲状腺癌、良性结节性甲状腺肿、Graves病、淋巴细胞性甲状腺炎和原发性甲状腺功能减退的遗传位点，解析疾病间遗传相关性和生物学通路，并评估PRS在良恶性结节鉴别和侵袭性甲状腺癌筛查中的临床价值。

技术路线

作者主要用了：GWAS元分析、cov-LDSC、TWAS、KEGG/Reactome通路分析、PRS。

作者联合19个生物库开展五类甲状腺疾病的GWAS元分析，共纳入约290万个体的基因型数据。首先进行变异发现和严格质控，随后通过遗传相关性分析、TWAS和通路富集解析功能机制，最后在独立生物库中构建并验证PRS，重点评估其区分良恶性结节和识别侵袭性甲状腺癌的能力。

图 1　Fig. 1。图1展示了VTB联盟的整体研究框架：从19个生物库的GWAS数据收集与质控，到元分析鉴定遗传位点，再到功能推断（遗传相关性、TWAS、通路分析）和临床转化（PRS开发与验证）。该图概括了从变异发现到临床应用的三阶段设计逻辑。

核心亮点

亮点 1｜883个位点，570个新发现

元分析共鉴定出883个独立显著位点（P≤5×10-8），其中313个为此前已报道的甲状腺性状相关位点，570个为新发现位点。多数先导变异位于内含子区（n=407），46个外显子变异中43个为错义突变，可能影响蛋白功能。例如TERF1罕见错义变异（8:73046129:G:A）仅在EUR-like分层分析中显著（β=1.32，P=1.08×10-9），而DIO1变异与甲状腺功能减退显著相关。这些发现大幅扩展了甲状腺疾病的遗传图谱。

图 2　Fig. 2。图2的热图展示了五种甲状腺疾病之间的遗传相关性，星号表示FDR<0.05的显著相关。右侧圆形图突出与甲状腺癌和良性结节显著相关的位点，红蓝点分别代表主要与甲状腺癌或良性结节相关的位点，PTCSC2位点（黄色）对两种疾病效应方向相反。左侧显示Graves病特异位点。

亮点 2｜甲状腺疾病间遗传相关性

遗传相关性分析显示淋巴细胞性甲状腺炎与原发性甲状腺功能减退几乎完全相关（rg=0.97），与Graves病相关性为0.62。甲状腺癌与良性结节性甲状腺肿的遗传相关性为0.41，提示两者既有共享遗传基础，也存在疾病特异性成分。Graves病与甲状腺癌的遗传相关性为0.20，与良性结节为0.31。这些结果为理解甲状腺疾病谱系的遗传结构提供了定量证据。

图 3　Fig. 3。图3展示了PRS在CCPM生物库中的表现。a部分显示PRSThC versus All的AUC为0.692，PRSThC versus BNG的AUC为0.670。b部分显示按PRS十分位分层的甲状腺癌风险，最高十分位风险显著升高。c部分显示PRS与侵袭性特征（如肿瘤大小、多灶性、淋巴结转移等）的关联，星号表示Bonferroni校正后显著。

亮点 3｜细胞周期与DNA修复基因驱动甲状腺癌

甲状腺癌特异性位点富集于细胞周期检查点、中心体/动粒功能、微管附着和染色体分离相关基因（如CDCA7L、CENPE、CHEK2、TP53等），以及DNA修复和损伤应答基因（ATM、DCLRE1B、EXO1等）。而良性结节特异性位点则涉及IGF1和FGF信号通路。端粒维持基因（TERT、STN1等）同时与良恶性结节相关。通路分析证实细胞周期、衰老和凋亡是甲状腺癌风险的关键生物学过程。

图 4　Fig. 4。图4提出甲状腺结节形成的双过程模型：增生（多克隆滤泡细胞增殖，无恶性潜能）和肿瘤形成（体细胞遗传改变驱动的克隆性生长）。GWAS表型定义（良性与恶性结节）与生物学机制（增生与肿瘤）之间的错配导致了表观遗传多效性。图中标注了与BNG相关但与ThC无关的通路和基因（如胰岛素信号通路）。

亮点 4｜PRS区分良恶性结节并预测侵袭性

在CCPM生物库中，PRSThC versus All区分甲状腺癌与对照的AUC为0.692，最高十分位个体患甲状腺癌风险是最低十分位的10.7倍。针对临床鉴别良恶性结节的任务，PRSThC versus BNG的AUC为0.670，最高十分位结节为恶性的风险是最低十分位的7.8倍。PRSThC versus BNG与结构性复发风险、肿瘤大小、多灶性、淋巴结转移和结外侵犯显著相关。模拟筛查侵袭性甲状腺癌时，PRSThC versus All的AUC达0.741，敏感性0.803。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： GWAS元分析、cov-LDSC、TWAS、KEGG/Reactome通路分析、PRS
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事甲状腺疾病遗传学、癌症风险预测或PRS临床转化的研究者。

带走一句

大规模多族裔GWAS结合功能注释和PRS验证，可系统解析复杂疾病的遗传异质性，并为临床风险分层提供可转化的遗传工具。