目的 采用全外显子组测序(WES)结合靶向分析识别中国人群甲状腺乳头状癌(PTC)新的易感基因变异。方法 选取2018—2020年北京医院行手术治疗的PTC患者286例,用于识别PTC候选基因的变异位点。选取2021—2022年北京医院行手术治疗的PTC患者214例,用于验证变异位点的突变频率。采用WES检测所有研究对象外周血白细胞基因组DNA,并采用靶向分析识别全基因组关联研究(GWAS)定位的28个PTC候选基因的蛋白编码区,及其侧翼选择性剪接序列的罕见变异。检测变异携带者癌组织和癌旁组织中变异位点等位基因的表达情况。通过病例-对照关联分析评估变异位点对PTC发生风险的影响。对照人群的基因分型数据分别来自gnomAD数据库、TOPMed数据库、ChinaMAP数据库和HUABIAO数据库。结果 WES结合靶向分析识别出红细胞膜蛋白4.1样4A(EPB41L4A)基因上的2个极罕见错义突变(rs1455421289位点c.163G>C:p.D55H和rs761977647位点c.855G>C:p.W285C),可能潜在影响编码蛋白的功能。在用于识别PTC候选基因变异位点的286例PTC患者中,有3例携带rs1455421289位点c.163G>C:p.D55H变异,1例携带rs761977647位点c.855G>C:p.W285C变异。这2个位点在4个对照人群中的携带频率均<0.05%。在用于计算变异位点突变频率的214例PTC患者中,有1例携带rs1455421289位点c.163G>C:p.D55H变异,未发现rs761977647位点c.855G>C:p.W285C变异。EPB41L4A rs1455421289的突变频率为0.8%(4/500),rs761977647的突变频率为0.2%(1/500)。EPB41L4A rs1455421289变异导致突变型mRNA的表达丰度显著高于野生型mRNA(P<0.01)。与4个对照人群比较,携带rs1455421289位点c.163G>C:p.D55H变异者PTC发生风险显著升高[比值比(OR)]值分别为40.3、85.4、213.5、541.1,95%可信区间(CI)分别为4.5~361.5、9.5~765.4、57.2~797.2、98.9~2 960.8)。结论 基于WES识别出1个新的PTC易感基因变异EPB41L4A rs1455421289,该变异位点携带者PTC的发生风险显著增高。