中国南方人群TGFβ3和ARNT2基因多态性与非综合征性唇腭裂的关联研究

[摘要] 目的 探讨TGFβ3和ARNT2基因与非综合征性唇腭裂发生的关系。 方法 收集了2006～2010年的119

>> SKI基因多态性与中国华东地区部分人群非综合征性唇腭裂的相关研究 宁夏回汉族人群Wnt3基因单核苷酸多态性与非综合征型唇腭裂相关性的研究 中国人群5，10―亚甲基四氢叶酸还原酶基因C677T多态性与非综合征性唇腭裂关系的Meta分析 中国东北人群中转化生长因子―α基因rs3771494单核苷酸多态性与非综合征唇腭裂的关系 中国南方人内脏最健康 TBX22基因多态性与非综合征性唇腭裂相关关系的研究 Msx1基因多态性与非综合征性唇腭裂 IRF6基因多态性与非综合征性唇腭裂的相关性研究进展 不同季节南方人工改良草地载畜量研究 “进步的总统”还是“狡猾的南方人”？ 《南方人物周刊》的报道特色 一个南方人眼中的哈尔滨 ORMDL3基因多态性与南方成人哮喘发生的相关研究 南方人VS北方人：我们是如何相爱相杀的 右江流域壮族人群酒精依赖程度与ADH和ALDH基因多态性的关联性研究 盛南方:南方人才市场的掌舵人 MTHFR基因与非综合征性唇腭裂的关系研究进展 《南方人物周刊》人物报道中叙述主体的特点 《南方人物周刊》封面人物报道的叙事策略 为家乡效力还是为联邦效力南方人的抉择 常见问题解答 当前所在位置：l）M行传递失衡（TDT）检验。统计时将患病同胞拆分为独立样本家系；采用以家庭为基础的关联检验（FBAT）模型，应用统计软件fbat 1.41，在不同的遗传模型下，以子代，父母基因型作为基因型数据，以子代是否受累作为表型数据，对SNP位点与非综合征性唇腭裂表型之间是否存在关联进行线性回归模型检验。

2 结果

2.1 各位点的基因型分析

TGFβ3基因rs2300607的多态性表现为纯合子TT（140、107、45 bp），纯合子AA（247、45 bp），杂合子TA（245、140、107、45 bp）。rs2268625的多态性表现为纯合子TT（245 bp），纯合子CC（65、180 bp），杂合子TC（245、180、65 bp）；ARNT2基因rs2228099的多态性表现为纯合子GG（264 bp），纯合子CC（236、29 bp），杂合子GC（264、236、29 bp）。rs3768017的多态性表现为纯合子TT（408 bp），纯合子GG（379、29 bp），杂合子TG（408bp、329、29 bp）。

2.2 Hardy-Weinberg（H-W）平衡检验

对于汉族CLO和CLP核心家系，患儿和父母各基因型的观察值与期望值间吻合度良好，分布未偏离H-W平衡（P > 0.05）（表2）；而对于汉族CPO核心家系，有三个位点的基因型不符合H-W平衡，不能进行后面的分析。其他少数民族家系数较少，且民族众多，分层后统计差异无统计学意义（故以下统计数据均来自汉族家系）。

表2 汉族唇裂、唇腭裂和单纯腭裂家系中两基因各SNP位点的

H-W平衡检验

2.3 传递不平衡检验（TDT）

进行基因分型的核心家系有119个，由于TDT要求运用杂合子父母及其患病子女做分析，故需排除纯合子父母家系，余下有效家系106个用于统计分析。结果显示TGFβ3的rs2268625等位基因与唇腭裂的发生有显著关联（P = 0.02），ARNT2的rs2228099等位基因与单纯腭裂的发生有显著关联（P = 0.04）。（TDT检验只针对符合H-W平衡的SNP位点进行，具体数据见表3～4）。

2.4 各SNP位点的FBAT分析结果

在FBAT加性遗传模式下，这四个SNP位点的单点分析结果与上面的TDT分析结果非常接近（表5）。在用FBAT分析TGFβ3基因与ARNT2基因的单体型数据时，将汉族CLO和CLP家系合并统计以获得更多的信息家系数，由于FBAT程序中默认设置最小的信息家系数为10，信息家系数低于10的主题（等位基因、基因型和单体型）将不被分析。结果显示TGFβ3基因2个SNP位点构建的4个单体型中，TT单体型的频率最高（0.529），TC单体型的频率最低（0.047），但对该单体型的分析获得一个有意义的结果（P = 0.045）（表6）；ARNT2基因2个SNP位点构建的4个单体型中，CT单体型和GG单体型都与CL/P的发生存在显著关联（P = 0.008、0.012）（表7）。

3 讨论

先天性唇腭裂发生在胚胎发育6～12周时，在致畸因素的作用下，面突、腭突的外胚间质细胞生长减慢或停止，上皮细胞与间充质细胞之间的相互作用受到干扰，腭中嵴上皮细胞不能正常增生、分化、凋零、退变，而是分化为口腔样的复层鳞状上皮，面突与腭突就只能部分融合或不能融合，出现不同程度的唇腭裂[8]。这个过程受到体内多种细胞因子及信息分子的调节，TGFβ族作为一类多肽生长因子也参与其中，动物实验表明，在胚腭突接触前，TGFβ3诱导腭胚中嵴上皮细胞膜外出现丝状伪足，丝状伪足表面有硫酸软骨蛋白多糖，介导腭突黏着。在腭融合时TGFβ3可能作为凋亡的T导因子调节腭中嵴上皮细胞的转归，还可诱导腭成纤维细胞中基质金属蛋白酶-13的高度表达，而基质金属蛋白酶降解细胞外基质是腭融合中的关键步骤[9-10]。ARNT2是近年来被发现与唇腭裂的发生相关的侯选基因，它位于人类15q23-q25染色体，Abbott等[11]证实在C57BL/6N小鼠早期腭发育的不同时期中，ARNT mRNA有不同程度的上调表达。Barrow等[12]、Muraoka等[13]对AXB和BXA近交系小鼠进行基因组扫描，证实小鼠7号染色体上的ARNT2等位基因是CPO致畸易感位点，ARNT2基因敲除小鼠会产生包括CPO在内的畸形。Tami等[14]发现，TGFβ3或其信号通路中的其他因子与ARNT2基因编码蛋白所在的芳香烃受体通路有关，在暴露于某些环境因素时，两者之间有相互补充影响的作用。

关于TGFβ3基因与唇腭裂发生关系的研究最早在高加索人群中进行，如Lidral等[15]分析美洲高加索CL/P和CPO患者人群中的TGFβ3基因，发现第5外显子区X5.1与CL/P之间存在传递不平衡，X5.1-5′UTR单体型与CL/P之间也存在显著的传递不平衡。Maestri等[16]分析了160个家系受累患者唇腭裂的类型、人种、家族史与TGFβ3基因标记之间的关系，也发现TGFβ3与CL/P患者和CPO患者呈显著性相关。但Moreno等[17]和Beaty等[18]分析TGFβ3基因上的短串联重复序列，却未发现TGFβ3基因与唇腭裂发生存在关联或连锁。上述研究都是用微卫星遗传标记筛查法来寻找易感基因，使用SNP这类新型遗传标记寻找某基因与唇腭裂发生关系的研究较少。

我国人群中，研究TGFβ3基因与非综合征性唇腭裂关系的实验进行最多。众多的实验结论或相同或相悖，但设计方法多采用的是病例对照法，这种设计很难保证病例组和对照组来自遗传背景完全相同的人群，可能会因为人群分层（病例和对照在遗传背景上的差异）而导致虚假的关联结果。而病例父母对照研究则可避免这种差异，因为其本身就是患儿与正常父母之间的对照，可以得到更为真实的结果。FBAT检验是利用双亲基因型，受累子女基因型及其表型的家系数据，进行双等位基因或多等位基因与疾病表型的关联分析，它是由传递不平衡检验（TDT）扩展到对任何表型，多基因或环境因素的研究。其主要优势在于能有力地克服人群混合、人群分层以及性状的不同选择等混杂因素对分析结果的影响，允许在一个家系中存在多个受累对象，并能在不同的遗传模式下进行相关性检验分析等。本文就用这两种方法分析我国南方人群中TGFβ3和ARNT2基因与非综合征性唇腭裂的关系。

本文中TDT和FBAT单点分析结果都提示，TGFβ3基因的rs2268625位点与CLP的发生存在显著关联（P = 0.02），从亲代到子代该位点T等位基因传递更多；但该位点与CLO的发生没有关系。在对TGFβ3的单体型分析中，虽然TC单体型可能会增加汉族人CL/P发生的风险（P = 0.047），但该单体型的频率低于0.05，需要扩大样本量来增加说服力，再做进一步的解释。这与Eisaburo等[19]对112个日本CL/P核心家系的研究结果略有不同，他们发现TGFβ3基因的rs2300607（A>T）在日本人中存在显著传递不平衡（P = 0.041），rs2300607与rs2268625组成的单体型在TDT检验中得到的P值更低（P = 0.025），在日本人群中，这两个SNP位点构成的单体型比单个SNP位点与CL/P的关联性更强。造成这种现象的原因需要进一步探讨，目前认为，可能与人群中基因频率的不同有关，该位点等位基因A在日本CL/P患者[20]中的频率是0.37，而在中国人群CL/P患者中的频率则是0.42。

ARNT2的SNP位点rs2228099与CLO的发生有显著关联（P = 0.04），rs3768017有增加CLP发生的趋势，但未达到显著性水平（P = 0.07），我们发现，在FBAT检验中rs2228099与rs3768017构建的CT单体型和GG单体型比单个SNP位点得到的P值更低（P = 0.008、0.012），与CL/P发生存在更强的关联性。这些结果说明ARNT2是中国汉族人发生CL/P的易感基因，而ARNT2上与CL/P的发生关联性更强的SNP位点，应该位于与rs2228099/rs3768017单体型相连锁的区域。

本次实验发现，TGFβ3基因与我国南方汉族人群中非综合征性唇裂合并腭裂的发生有一定的相关性，但它是否是CL/P和CPO发生的主要易感基因仍有待研究。非综合征性唇腭裂是一种多基因疾病，可能还有其他在我国人群唇腭裂发生中起主导或辅助作用的致病基因。我们的后续实验将继续讨论其他相关基因与唇腭裂发生的关系，为探讨唇腭裂的发病机制奠定基础。

[参考文献]

[1] 梁娟，王艳萍.我国南方非综合征性总唇裂的流行病学[J].临床口腔医学杂志，2000，16（2）：122-123.

[2] Murray JC. Face facts：genes，environment，and cleft [J]. Am J Hum Genet，1995，57（2）：227.

[3] Stanier P，Moore GE. Genetics of cleft lip and palate：syndromic genes contribute to the incidence of non-syndromic clefts [J]. Hum Mol Genet，2004，13 Spec No 1：R73-R81.

[4] Loffredo LC，Souza JM，Freitas JA. Oral clefts and vitamin in supplementation cleft palate [J]. Craniofac J，2001，38（1）：76.

[5] Mitchell LE，Risch N. Mode of inheritance of nonsyndromic cleft lip with or without cleft palate：a reanalysis [J]. Am J Hum Genet，2007，51（2）：323-332.

[6] Prescott NJ，Winter RM，Malcolm S. Nonsyndromic cleft lip and palate：complex genetics and environmental effects [J]. Ann Hum Genet，2001，65（Pt 6）：505-515.

[7] Vieira AR，Orioli IM. Candidate genes for nonsyndromic cleft lip and palate [J]. ASDC J Dent Child，2001，68（4）：272-279，229.

[8] 东英，张建波，金岩.胚胎发育期TGF-β超家族成员与腭裂发生的关系[J].国外医学：口腔医学分册，2004，31（2）：132-134.

[9] Britto JA，Evans RD，Hayward RD，et al. Toward pathogenesis of Apert cleft palate：FGF，FGFR，and TGF beta genes are differentially expressed in sequential stages of human palatal shelf fusion [J]. Cleft Palate Craniofac J，2002，39（3）：332-340.

[10] Nugent P，Ma L，Greene J. Differential expression and biological activity of retinoic acid-induced TGFB isoforms in embryomic palatemesenchymal cells [J]. J Cell Physiol，1998，177（1）：36-46.

[11] Abbott BD，Schmid JE，Brown JG. RT-PCR quantification of AHR，ARNT，GR，and CYP1A1 mRNA in craniofacial tissues of embryonic mice exposed to 2，3，7，8-tetrachlorodibenzo-p-dioxin and hydrocortisone [J].Toxicol Sci，1999，47（1）：76-85.

[12] Barrow LL，Wines ME，Romitti PA，et al. ryl hydrocarbon receptor nuclear translocator 2 （ARNT2）：structure，gene mapping，polymorphisms，and candidate evaluation for human orofacial clefts [J]. Teratology，2002，66（2）：85-90.

[13] Muraoka N，Shum L，Fukumoto S，et al. Transforming growth factor-beta3 promotes mesenchymal cell proliferation and angiogenesis mediated by the enhancement of cyclin D1，Flk-1，and CD31 gene expression during CL/Fr mouse lip fusion [J]. Birth Defects Res A Clin Mol Teratol，2005，73（12）：956-965.

[14] Tami LT，Emily AS，Christopher AB. Transforming growth factor-β3 restores fusion in palatal shelves exposed to 2，3，7，8-tetrachlorodibenzo-p-dioxin [J]. J Biol Chem，2005，280（13）：12742-12746.

[15] Lidral AC，Romitti PA. Association of MSX1 and TGFB3 with nonsyndromic clefting in humans [J]. Am J Huma Genets，1998，63（2）：557-568.

[16] Maestri NE，Beaty TH，Hetmanski J. Application of trans-mission disequilibrium tests to nonsyndromic oral clefts：including candidate genes and environmental exposures in the models [J]. Am J Med Genet，1997，73（3）：337-344.

[17] Moreno LM，Burgos MA，Marazita ML. Genetic analysis of candidate loci in nonsyndromic cleft lip families from Antioquia-Colombia and Ohio [J]. Am J Med Genet A，2004，125A（2）：134-144.

[18] Beaty TH，Wang H，Hetmanski JB. A case-control study of nonsyndromic oral clefts in Maryland [J]. Ann Epidemiol，2001，11（6）：434-442.

[19] Eisaburo I，Akira W，Yoko N. PAX9 and TGFB3 are linked to susceptibility to nonsyndromic cleft lip with or without cleft palate in the Japanese：population-based and family-based candidate gene analyses [J]. J Hum Genet，2006，51（1）：38-46.

[20] Tanabe A，Taketani S，Endo-Ichikawa Y，et al. Analysis of the candidate genes responsible for non-syndromic cleft lip and palate in Japanese people [J]. Clin Sci（Lond），2000，99（2）：105-111.

（收稿日期：2016-12-23 本文辑：李亚聪）