TYK2基因与家兔肠炎的易感性研究

付陆 赵明德 吴小燕 贾先波 赖松家

摘要 [目的] 研究TYK2基因在家兔肠炎发生中的遗传效应，探讨其是否可以作为家兔抗病育种的辅助选择标记。[方法] 通過构建病例组-对照组肠炎兔群和低纤维诱导的肠炎兔群，PCR产物纯化后直接测序，以及qRT-PCR法检测TYK2基因在回肠和结肠组织中的mRNA表达水平。[结果] 发现TYK2基因外显子区域有5个cSNP，只有c.1477（c.1477，C>T）位点发生非同义突变，导致氨基酸p.Leu 404 Phe（L>F）的改变，C等位基因增加了肠炎的易感性（OR：1.36；95% CI 1.269～2.151；P=0.019），而等位基因T对肠炎的发生起到了一定的保护作用（OR：0.81，95% CI 0.352～0.960；P=0.018）。TYK2基因不同基因型和不同肠炎程度的mRNA水平均呈现差异表达（P<0.05）。[结论] TYK2（c.1477，C>T）是家兔肠炎易感风险基因，为家兔的抗病育种提供了一个可靠的辅助选择标记。

关键词 TYK2；肠炎；遗传易感性

Investigation of Genetic Susceptibility of TYK2 Gene in Rabbit Enteritis

FU Lu1， 2， ZHAO Ming-de1， WU Xiao-yan1， LAI Song-jia2* et al

（1.Laboratory Animal Centre， Southwest Medical University， Luzhou， Sichuan 646000；2.Institute of Animal Genetics， Sichuan Agricultural University，Chengdu， Sichuan 611130）

Abstract [Objective] To study the genetic effects of TYK2 gene in rabbit enteritis and to explore whether it could be used as a Marker-Assisted Selection （MAS） for resistance breeding in rabbits. [Method] To construct case-control rabbit enteritis groups and low fiber induced rabbit enteritis groups， the PCR products were purified and sequenced directly， and the TYK2 gene mRNA levels in ileum and colon tissues were detected by qRT-PCR. [Result] Five cSNPs were detected in exon region of TYK2 gene， and only c.1477 （c.1477， C>T） loci was non synonymous mutation， which caused to the change of amino acid p.Leu 404 Phe （L>F）. The C allele increased the susceptibility of enteritis （OR： 1.36；95% CI 1.269～2.151；P=0.019）， and the T allele was a protective effect on the occurrence of the enteritis （OR： 0.81， 95% CI 0.352～0.960；P=0.018）. The TYK2 mRNA levels in different genotypes and in different degree of enteritis both presented significant difference （P<0.05）. [Conclusion] TYK2 （c.1477， C>T） is a susceptible risk gene of rabbit enteritis， and provides a MAS for the resistance breeding.

Key words TYK2；Enteritis；Genetic susceptibility

家兔为人类提供大量的优质肉、毛和皮等产品，随着消费的迅速增长，规模化养兔成为发展的必然趋势。然而，薛家兵等[1]、谷子林[2]报道患黏液性肠病的兔绝大部分死亡，发病率50%～70%，死亡率90%以上，一些兔场发病死亡率甚至高达100%。家兔腹泻死亡的致病机理十分复杂，防治难度很大，被列入近几年三大兔病之一（兔瘟、兔肠道疾病和体表真菌病）[3]。因此，研究家兔腹泻和肠炎的发病机制，从而控制死亡率已经成为当前养兔业中抗病育种的热点。

TYK2（Tyrosine kinase 2）作为酪氨酸蛋白激酶家族（Janus kinases，JAKs）的一个重要成员，参与机体细胞的增殖、分化、凋亡[4]，介导机体免疫调控[5-6]和肿瘤生成等一系列生命活动过程[7-8]。Minegishi等[5]研究发现TYK2基因纯合子突变的病人被临床诊断为高IgE综合征，表明其在多种细胞因子的天然免疫和获得性免疫中扮演重要作用。研究表明TYK2可以通过一种棕色脂肪组织（BAT）的分化来调控小鼠肥胖[9]。既然TYK2基因广泛参与到上述生理生化反应中，它是否在家兔肠炎发生中也扮演一定的作用，值得探索。

采用病例组-对照组试验设计拟在480只新西兰兔中，筛选TYK2基因的编码区SNP（cSNP），并剖析其中重要cSNP在家兔肠炎发生中的遗传效应；在低纤维日粮（CF=8.96%）诱导的60只新西兰肠炎患病兔中，荧光定量PCR法检测TYK2基因在回肠和结肠组织中的mRNA差异表达水平，为家兔抗病育种提供分子辅助选择标记。

1 材料与方法

1.1 试验兔群的构建及样本采取

选取四川农业大学教学试验兔场饲养的新西兰兔为试验材料，随机选择35日龄健康无病、三代内无亲缘关系的480只新西兰兔为试验动物，构建病例组-对照组肠炎易感性试验兔群。选取三代内无亲缘关系、胎次相近、健康无病、体重相当的35日龄新西兰断奶仔兔60只，作为低纤维日粮肠炎试验兔群，饲喂低纤维饲料（CF=8.96%）。上述2个群体构建的具体方法、分组情况及组织样品采取参考Zhang等[10]的方法。

1.2 DNA和RNA的提取及检测

试验采用DNA提取试剂盒AxyPrepTM Muitisource Genomic DNA Miniprep Kit（Axygen Biosciencs，USA），按照试剂盒说明书要求，从家兔耳组织中提取全基因组DNA。试验采用Trizol法（Invitrogen），严格按照其说明提取回肠和结肠组织总RNA。DNA和RNA的检测方法相同，先进行1.5%琼脂糖凝胶电泳检测，后用核酸蛋白仪进行OD260/OD280值测定。

1.3 PCR扩增及产物测序

根据NCBI提供的家兔TYK2基因参考序列（XM_017338771.1），设计3对引物（表1）进行编码区的扩增。采用10 μL的反应体系：2×Taq PCR Mix（5 μL）、Forward primer/Reverse primer（各0.3 μL）、Template cDNA（0.8 μL）和ddH2O（3.6 μL）。将上述溶液混合后，按以下条件进行PCR反应：94 ℃变性3 min，94 ℃变性30 s、退火（具体的退火温度由梯度PCR结果决定）30 s、72 ℃延伸60～90 s，扩增35～38个循环，72 ℃终延伸10 min，4 ℃保存。将PCR产物进行1.5%琼脂糖凝胶电泳检测，胶回收纯化后得到的PCR产物送至北京六合华大基因测序中心测序。

1.4 实时荧光定量PCR

先使用反转录试剂盒PrimeScriptTM RT Reagent Kit with gDNA Eraser（Cat.# RR047Q，TaKaRa），进行cDNA的合成，具体操作按试剂盒的说明书进行。然后采用2-△△CT法，以HPRT和GAPDH基因作为双内参，在CFX96TM Real-time System上检測TYK2基因在回肠和结肠组织中的mRNA表达水平。根据TYK2基因参考序列，设计荧光定量引物（表1），由北京六合华大基因有限公司合成。按照SYBR Premix Ex TaqTM II（TaKaRa，Dalian）试剂盒的使用说明书加样，使用标准品作为板间对照，并设NTC，每个待测样3个重复。

1.5 数据处理

基因分型的个体数直接以n表示，mRNA相对表达量表示为平均值±标准误。使用SAS 9.2软件分析TYK2基因与肠炎的遗传效应。运用SPSS 19.0分析处理，组间的差异比较采用单因素方差分析。P<0.05表示其差异具有统计学意义。

2 结果与分析

2.1 低纤维日粮诱导家兔肠炎炎症程度的分类

在60只低纤维日粮诱导家兔肠炎中（试验期间死亡和试验结束后屠宰），临床症状评分（CSS）为2和3分的有40只，其中35只的剖检症状分析评分（GLS）同样较高，主要表现为食欲下降、腹部胀气、盲肠内容物结块，大部分个体回肠内有胶冻状物质（22/35），部分个体结肠伴随着回肠出现胶冻状物质（17/35）。观察组织病理切片发现各肠断在发生严重肠炎时，都存在组织损伤和炎性损伤，并且回肠和结肠段病变得更为严重（图1）。患病家兔的回肠和结肠段，肠腔中未见红细胞和淋巴细胞内容物；肠绒毛多处毛细血管充血；部分肠上皮细胞坏死、溶解，可见少量中性粒细胞浸润（图1b和d中①标注）；肠绒毛中肠腺细胞发生坏死性溶解，可见嗜碱性粒细胞浸润（图1b和d中②标注）；黏膜固有层、黏膜下层结缔组织间的间质有炎性细胞浸润（图1d中③标注）。最终，将低纤维诱导试验兔群分为健康组（12只）、轻微肠炎组（13只）和严重肠炎组（35只）。

2.2 TYK2基因编码区扩增和直接测序结果

对新西兰兔TYK2基因的外显子区域，分别设计引物进行PCR扩增，扩增产物特异性好，用1.5%琼脂糖凝胶电泳检测PCR产物，电泳条带单一，明亮致密，大小与目的条带一致。将目的基因PCR扩增产物胶回收后，送至北京六合华大基因公司进行直接测序，共发现5个cSNP，分别为c.1095（exon7，G>A），c.1477（exon9，C>T），c.2013（exon13，C>T），c.2133（exon14，T>C），c.2166（exon14，C>T），其中c.1477位点发生非同义突变，导致氨基酸p.Leu 404 Phe（L>F）的改变，其余的cSNP均为同义突变。

2.3 c.1477与家兔肠炎易感性分析

试验共选取了480只新西兰兔（病例组253，对照组227），对c.1477位点采用直接测序的方法进行基因分型，结果发现纯合CC型、TT型、杂合CT型在病例组和对照组中个体数分别为156/118，20/58，77/51。

群体遗传结构分析发现在病例组和对照组中，C等位基因均处于优势地位，其频率分别为76.88%和63.22%（P<0.05），而T等位基因分别为23.12%和36.78%（P<0.05）。2群体中基因型频率CC（61.66% vs.51.98%），TT（7.91% vs.25.55%），CT（30.43% vs.22.47%）差异均显著（P<0.05）。通过SAS 9.2软件计算等位基因C增加了肠炎的易感性（OR：1.36；95% CI 1.269～2.151；P=0.019），而等位基因T对肠炎的发生起到了一定的保护作用（OR：0.81，95% CI 0.352～0.960；P=0.018）。HWE评估检验发现病例组和对照组都遵循HWE平衡定律（P>0.05）。使用SNPStasts软件评估5种遗传模型对该突变的可靠性。根据AIC和BIC值的大小，确定隐性模型为这个位点的最适模型。利用这个模型，关联性分析表明C/C基因型增加了肠炎的易感性（OR：1.73；95% CI 1.870～2.800；P<0.001）（表2）。

[10] ZHANG G W，GAO L，CHEN S Y，et al.Single nucleotide polymorphisms in the FTO gene and their association with growth and meat quality traits in rabbits[J].Gene，2013，527（2）：553-557.

[11] SCHLOLAUT W，HUDSON R，ROEDEL H G.Impact of rearing management on health in domestic rabbits：A review[J].World rabbit science，2013，21（3）：145-159.

[12] CARABAO R，BADIOLA I，CHAMORRO S，et al.Review.New trends in rabbit feeding：Influence of nutrition on intestinal health[J].Spanish journal of agricultural research，2008，6（S1）：15-25.

[13] BENNEGADI N，GIDENNE T，LICOIS D.Impact of fibre deficiency and sanitary status on non-specific enteropathy of the growing rabbit[J].Animal research，2001，50（5）：401-414.

[14] HALIGUR M，OZMEN O，DEMIR N.Pathological and ultrastructural studies on mucoid enteropathy in New Zealand rabbits[J].Journal of exotic pet medicine，2009，18（3）：224-228.

[15] VIGNAL A，MILAN D，SANCRISTOBAL M，et al.A review on SNP and other types of molecular markers and their use in animal genetics[J].Genetics selection evolution，2002，34（3）：275-306.

[16] ZHANG G W，WANG H Z，CHEN S Y，et al.A reduced incidence of digestive disorders in rabbits is associated with allelic diversity at the TLR4 locus[J].Vet immunol immunopathol，2011，144（3/4）：482-486.