目的:在中国人非综合征型唇裂伴或不伴腭裂(non-syndromic cleft lip with or without palate,NSCL/P)核心家系中,利用全外显子组测序探索与NSCL/P发病相关的新生突变位点.方法:对22个中国NSCL/P核心家系进行全外显子组测序,采用基因组分析工具包(Genome Analysis ToolKit,GATK)通过对比亲代与子代同一位点的等位基因识别新生突变位点,采用SnpEff软件对位点进行功能注释.对新生突变位点进行富集分析,检验全外显子区域内存在的新生突变数量是否高于预期值,以及是否存在包含新生突变数量显著高于预期值的基因.通过查阅文献总结既往研究提示与NSCL/P发病存在较强证据支持的基因,根据注释信息筛选能够引起蛋白质改变的新生突变位点,对该类位点所在基因编码的蛋白质与NSCL/P相关基因编码的蛋白质进行交互作用分析.利用R软件的denovolyzeR包进行富集分析(Bonferroni多重检验校正:P=0.05/n,n为基因个数).利用STRING数据库预测新生突变所在基因与已知NSCL/P致病基因编码的蛋白质间的交互作用.结果:全外显子组测序得到的位点中共有339 908个位点通过质量控制,经GATK软件比对共筛选出345个高置信度新生突变,其中错义突变44个,无义突变1个,经典剪接位点2个,同义突变20个,内含子区或基因间区位点278个.富集分析显示,全外显子组中引起蛋白质改变的新生突变数量显著高于预期值(P＜0.05),KRTCAP2、HMCN2、ANKRD36C、ADGRL2和DIPK2A 5个基因所含的新生突变位点高于预期(P＜0.05/(2×19 618)).蛋白质交互作用分析纳入46个包含能够引起蛋白质序列改变的新生突变所在的基因及13个既往研究提示与NSCL/P存在关联的基因,两类基因编码的蛋白质之间存在6组交互作用,其中RGPD4与SUMO1编码的蛋白质的交互作用证据可信度最高,STRING数据库交互作用评分为0.868.结论:研究为NSCL/P的发病提供了新的证据,对携带新生突变的基因进行功能分析有助于揭示复杂疾病的遗传结构.

新生突变;富集分析;蛋白质交互作用;非综合征型唇裂伴或不伴腭裂

陈曦;王斯悦;薛恩慈;王雪珩;彭和香;范梦;王梦莹;武轶群;秦雪英;李劲;吴涛;朱洪平;李静;周治波;陈大方;胡永华

北京大学公共卫生学院流行病与卫生统计学系,北京 100191;北京大学口腔医学院·口腔医院口腔颌面外科,国家口腔医学中心,国家口腔疾病临床医学研究中心,口腔生物材料和数字诊疗装备国家工程研究中心,口腔数字医学北京市重点实验室,国家卫生健康委员会口腔医学计算机应用工程技术研究中心,国家药品监督管理局口腔生物材料重点实验室,北京 100081;北京大学口腔医学院·口腔医院儿童口腔科,国家口腔医学中心,国家口腔疾病临床医学研究中心,口腔生物材料和数字诊疗装备国家工程研究中心,口腔数字医学北京市重点实验室,国家卫生健康委员会口腔医学计算机应用工程技术研究中心,国家药品监督管理局口腔生物材料重点实验室,北京 100081

北京大学学报（医学版）

[1]陈曦;王斯悦;薛恩慈;王雪珩;彭和香;范梦;王梦莹;武轶群;秦雪英;李劲;吴涛;朱洪平;李静;周治波;陈大方;胡永华-.基于核心家系全外显子组测序数据探索新生突变与非综合征型唇腭裂的关联)[J].北京大学学报（医学版）,2022(03):387-393

非综合征型唇裂,非综合征型唇裂伴或不伴腭裂,ToolKit,SnpEff,denovolyzeR,KRTCAP2,HMCN2,ANKRD36C,ADGRL2,DIPK2A,RGPD4

家系,全外显子组测序,序数,数据探索,新生突变,非综合征型唇腭裂,syndromic,cleft,lip,or,without,palate,NSCL,突变位点,基因组分析,工具包,Genome,Analysis,GATK,亲代,子代,等位基因,基因识别,功能注释,富集分析,子区域,变数,预期值,证据支持,基因编码,交互作用分析,Bonferroni,多重检验,STRING,致病基因,点中,个位,软件比对,置信度,错义突变,无义突变,剪接,同义突变,内含子,所含,蛋白质交互作用,蛋白质序列,序列改变,存在关联,SUMO1,可信度,功能分析,遗传结构

0.233395