吳 磊, 覃媛鈺, 譚光輝, 李杰章, 張依裕, 蔣會梅, 楊遠清, 向程舉

(1.貴州大學(xué)高原山地動物遺傳育種與繁殖教育部重點實驗室，貴州 貴陽 550025；2.貴州大學(xué)貴州省動物遺傳育種與繁殖重點實驗室，貴州 貴陽 550025；3.貴州大學(xué)動物科學(xué)學(xué)院，貴州 貴陽 550025；4.畢節(jié)市畜禽遺傳資源管理站，貴州 畢節(jié) 551700)

IL-8屬于C-X-C亞家族，是一種非常重要的趨化因子，具有介導(dǎo)細胞在炎癥部位聚集、活化以及修復(fù)損傷的組織，刺激中性粒細胞和淋巴細胞發(fā)生趨化游走等功能[1-2].趨化因子是一種普遍存在的小分子量趨化蛋白，通過免疫調(diào)節(jié)細胞在感染期間遷移到感染部位進而參與宿主防御[3].IL8R 屬于視紫紅質(zhì) G-蛋白偶聯(lián)受體超家族，有兩種類型的基因，分別為IL8RA(只與IL-8有高度親和力)和IL8RB(能高度結(jié)合其他趨化因子)[4-5].研究表明，2種受體的氨基和羧基顯著不同，氨基的不同導(dǎo)致2種受體功能和作用有大的特異性[6].雞的白細胞介素8受體(interleukin-8 receptor bata, IL8RB)基因定位于7號染色體，全長2 981 bp，由2個外顯子和1個內(nèi)含子組成.IL8RB基因具有介導(dǎo)活化中性粒細胞的作用，可促進中性粒細胞脫顆粒，釋放貯存酶，增強中性粒細胞的吞筮功能，啟動超氧離子釋放，導(dǎo)致機體局部發(fā)生炎癥反應(yīng)，達到殺滅病原菌的目的[7].Matsushima et al[8]檢測到IL8RB基因與許多實體腫瘤的耐藥性、侵襲和轉(zhuǎn)移相關(guān)，可能影響家禽抗逆性.Luppi et al[9]證實IL8RB基因通過受體酪氨酸激酶磷酸化激活生長因子信號反應(yīng)，能導(dǎo)致細胞增殖、遷移和轉(zhuǎn)移的增強.綜合前人研究結(jié)果可知，IL8RB基因是與某些疾病相關(guān)的基因，可能與家禽抗逆性有很大聯(lián)系，但關(guān)于畜禽IL8RB基因多態(tài)性的研究報道極為少見.

本研究以貴州地方雞種長順綠殼蛋雞、黔畫烏雞和威寧雞為研究對象，采用DNA直接單個測序篩選IL8RB基因的單核苷酸多態(tài)性(single nucleotide polymorphism, SNP)位點，計算其單倍型、雙倍型與蛋殼品質(zhì)的關(guān)聯(lián)性，探討IL8RB基因突變與黔畫烏雞肉品質(zhì)的聯(lián)系，旨在為今后進一步研究該基因?qū)﹄u肉品質(zhì)及生理功能的影響提供參考.

1 材料與方法

1.1 試驗動物及樣品采集

長順綠殼蛋雞110只(♂38，♀72)、黔畫烏雞55只(♂15，♀40)和威寧雞76只(♂26，♀50)均隨機采自貴州大學(xué)家禽研究所科研雞場.雞只均健康無病，通過翅靜脈采血和抗凝肝素鈉抗凝的方式采集血樣0.5～1.0 mL，按照DNeasy Blood & Tissue Kit(250)試劑盒說明書的方法逐只提取血樣基因組DNA進行1%瓊脂糖凝膠電泳檢測，采用NanoDrop 2000型核酸濃度檢測儀(美國NanoDrop公司)測定其濃度，并稀釋成終濃度為50～100 ng·μL-1，置-20 ℃保存?zhèn)溆?屠宰全部黔畫烏雞，取胸肌組織，參照李利等[10]的方法測定雞肉品質(zhì).

1.2 引物設(shè)計

根據(jù)GenBank收錄的雞IL8RB基因序列(NC-006094.5)，采用網(wǎng)上引物設(shè)計軟件Primer3 Input (version 0.4.0)設(shè)計3對引物，引物信息見表1.F1/R1擴增區(qū)域包括外顯子1、F2/R2和F3/R3擴增區(qū)域包括外顯子2，引物由昆泰銳(武漢)生物技術(shù)有限責(zé)任公司合成.

1.3 PCR擴增

3對引物采用相同的PCR反應(yīng)體系和擴增程序.PCR反應(yīng)體系：1 μL基因組DNA模板、上游引物(10 μmol·L-1)和下游引物(10 μmol·L-1)各1 μL、10 μL 2×Es Taq MasterMix、7 μL雙蒸水，總體積為20 μL.PCR擴增程序：94 ℃預(yù)變性5 min， 94 ℃變性40 s，59 ℃退火40 s，72 ℃延伸40 s，共35個循環(huán)，72 ℃終延伸6 min.待反應(yīng)結(jié)束后，把PCR反應(yīng)產(chǎn)物用移液槍加在1.5%瓊脂糖槽中進行凝膠電泳，電泳儀電壓設(shè)置為135 V，電泳時間設(shè)定為30 min，電泳完成后立即放到凝膠成像儀上拍照，查看結(jié)果并截圖保存.選取特異性較好的PCR擴增產(chǎn)物送往昆泰銳(武漢)生物技術(shù)有限責(zé)任公司測序.

表1 IL8RB基因的引物序列

1.4 mRNA二級結(jié)構(gòu)和蛋白質(zhì)高級結(jié)構(gòu)的預(yù)測

利用在線軟件RNAfold (http：//nibiru.tbi.univie.ac.at/cgi-bin/RNAWebSuite/RNAfold.cgi)對IL8RB基因突變前后編碼區(qū)序列mRNA的二級結(jié)構(gòu)進行預(yù)測；利用在線軟件SOPMA(https：//npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html)預(yù)測和分析IL8RB基因編碼區(qū)突變前后蛋白質(zhì)的二級結(jié)構(gòu)；利用在線軟件SWISS-MODEL(https：//swissmodel.expasy.org/interactive)預(yù)測IL8RB基因編碼區(qū)突變前后蛋白質(zhì)的三級結(jié)構(gòu).

1.5 統(tǒng)計分析

運用SHEsis網(wǎng)上在線軟件(http://analysis.bio-x.cn/)計算SNP位點的基因型頻率、等位基因頻率、單倍型頻率、基因型分布卡方值(χ2)、連鎖不平衡的D′和r2值.根據(jù)Nei et al[11]和Botstein et al[12]的報道計算觀察雜合度(He)、有效等位基因數(shù)(Ne)和多態(tài)信息含量(PIC).運用SPSS 19.0軟件中的一般線性模型分析SNP位點基因型或雙倍型與所測定性狀指標(biāo)的相關(guān)性，模型為Y=μ+G+e(Y為性狀觀測值，G為基因型效應(yīng)或雙倍型效應(yīng)，μ為群體均值，e為隨機殘差).采用最小顯著性差異法(LSD)進行多重比較，結(jié)果用平均值±標(biāo)準(zhǔn)差表示.

2 結(jié)果與分析

2.1 IL8RB基因SNP的篩查與鑒定

對IL8RB基因PCR產(chǎn)物進行直接測序，序列比對結(jié)果(圖1)表明，共發(fā)現(xiàn)2個SNP變異位點g.22589443 C>T和g.22589512 C>T，均處于IL8RB基因外顯子2的錯義突變(3個雞種變異位點一致).g.22589443 C>T突變共產(chǎn)生3種基因型，引起亮氨酸變成絲氨酸；g.22589512 C>T突變共產(chǎn)生2種基因型，引起蘇氨酸變成甲硫氨酸.

2.2 IL8RB基因SNP位點的遺傳特性

對IL8RB基因在長順綠殼蛋雞、黔畫烏雞和威寧雞鑒定到的3個SNP位點進行遺傳特性分析.結(jié)果(表2)顯示，3個雞品種g.22589443 C>T位點中的CC基因型和C等位基因分別為優(yōu)勢基因型和優(yōu)勢等位基因，雜合度和有效等位基因數(shù)較高，多態(tài)信息含量屬于中度多態(tài)位點，長順綠殼蛋雞g.22589443 C>T位點的基因型分布極顯著偏離了Hardy-Weinberg平衡(P<0.01)，而黔畫烏雞和威寧雞顯著偏離了Hardy-Weinberg平衡(P<0.05).3個雞品種g.22589512 C>T位點中的CC基因型和C等位基因分別為優(yōu)勢基因型和優(yōu)勢等位基因，雜合度和有效等位基因數(shù)較低，屬低度多態(tài)，基因型分布均未偏離平衡狀態(tài)(P>0.05).

表2 IL8RB基因SNP位點的群體遺傳特性1)

2.3 IL8RB基因SNP位點的連鎖不平衡分析

對長順綠殼蛋雞、黔畫烏雞和威寧雞IL8RB基因2個SNP位點g.22589443 C>T和g.22589512 C>T進行連鎖不平衡分析.結(jié)果(表3)顯示，2個SNP位點的D′值均為1，r2值小于0.33，根據(jù)Ardlie et al[13]和Slatkin[14]的報道，當(dāng)|D′|>0.8和r2>0.33時認為SNP位點間存在強連鎖不平衡，揭示了本研究發(fā)現(xiàn)的2個SNP位點間不存在強連鎖不平衡.

表3 IL8RB基因SNP位點的連鎖不平衡D′和r2值

2.4 IL8RB基因單倍型及雙倍型分析

對長順綠殼蛋雞、黔畫烏雞和威寧雞IL8RB基因的2個SNP位點進行單倍型和雙倍型分析.結(jié)果(表4)顯示，在檢測的3個雞群體IL8RB基因中，g.22589443 C>T和g.22589512 C>T位點共存在3種單倍型和4種雙倍型，單倍型H1(CC)和雙倍型H1H1(CCCC)出現(xiàn)的頻率最高，單倍型H3(CT)和雙倍型H2H2(TTCC)出現(xiàn)的頻率最低.

2.5 IL8RB基因生物信息學(xué)分析

就長順綠殼蛋雞、黔畫烏雞和威寧雞IL8RB基因外顯子2上的g.22589443 C>T和g.22589512 C>T位點組合對mRNA序列二級結(jié)構(gòu)的影響進行預(yù)測.結(jié)果(圖2)顯示，CC組合mRNA序列結(jié)構(gòu)的自由能為-2 415.74 kJ·mol-1，TC組合的自由能為-2 412.81 kJ·mol-1，CT組合的自由能為-2 431.23 kJ·mol-1，穩(wěn)定性為：CT>CC>TC，TC與GenBank上收錄的序列相同，揭示突變使自由能降低，穩(wěn)定性升高.對蛋白質(zhì)二級結(jié)構(gòu)和三級結(jié)構(gòu)分析發(fā)現(xiàn)，2個SNP位點引起氨基酸變化但并未引起蛋白質(zhì)結(jié)構(gòu)變化.

2.6 IL8RB基因變異位點對黔畫烏雞常規(guī)肉質(zhì)指標(biāo)的影響

IL8RB基因變異位點g.22589443 C>T和g.22589512 C>T與黔畫烏雞肉品質(zhì)的關(guān)聯(lián)性分析結(jié)果(表5)顯示， 2個SNP位點對黔畫烏雞肉品質(zhì)指標(biāo)的影響均不顯著(P>0.05),雙倍型H1H1(CCCC)個體的肉色指標(biāo)顯著高于H1H3(CCCT)個體(P<0.05)，其他雙倍型各指標(biāo)之間的差異均不顯著(P>0.05).

表4 IL8RB基因SNP位點的單倍型和雙倍型頻率

表5 IL8RB基因變異位點與黔畫烏雞常規(guī)肉質(zhì)指標(biāo)的關(guān)聯(lián)性1)

1)肉質(zhì)的各項指標(biāo)均在屠宰后40 min內(nèi)測定；同列數(shù)據(jù)后附不同字母者表示差異顯著(P<0.05)，附相同字母者或無字母者表示差異不顯著(P>0.05).

3 討論

隨著人們生活水平的快速提高，對禽品質(zhì)的要求也越來越高，通過提高家禽抗逆性，培育出安全優(yōu)質(zhì)的禽產(chǎn)品，可以滿足人們對高品質(zhì)禽產(chǎn)品的需求[15].通過多個數(shù)據(jù)庫檢索發(fā)現(xiàn)，IL8RB基因是與畜禽抗逆性有關(guān)的基因，但有關(guān)畜禽IL8RB基因的多態(tài)對生產(chǎn)性能的影響極少報道.IL8RB基因主要分布于中性粒細胞和髓樣前體細胞系中，具有介導(dǎo)活化中性粒細胞的作用，可促進脫顆粒，釋放貯存酶，增強吞噬功能，啟動超氧離子釋放，導(dǎo)致機體局部炎癥反應(yīng)，達到殺滅病原菌的目的[16-17]，IL-8真核表達質(zhì)粒作為免疫增強劑能夠顯著提高家禽的抗體效價水平[18].目前，已經(jīng)在人體IL8RA基因中鑒定出了SNP，并發(fā)現(xiàn)這些SNP與炎性疾病相關(guān)，并且已經(jīng)顯示其在不同的地理種群中有所不同[19-20].徐敏等[21]報道，CXCR2基因可以作為影響牛乳房炎和奶品質(zhì)性狀的可能主要候選基因，并將SNP+777和SNP+861兩個多態(tài)位點用于抗乳房炎性狀的選擇.本研究以貴州地方雞種長順綠殼蛋雞、黔畫烏雞和威寧雞為材料，對IL8RB基因進行SNP檢測，在IL8RB基因外顯子2上共篩選出2個SNP變異位點g.22589443 C>T和g.22589512 C>T.g.22589443 C>T突變共產(chǎn)生3種基因型，引起亮氨酸變成絲氨酸；g.22589512 C>T突變共產(chǎn)生2種基因型，引起蘇氨酸變成甲硫氨酸.g.22589443 C>T突變位點的雜合度和有效等位基因數(shù)較高，為中度多態(tài)，更利于遺傳選擇，在今后的遺傳選育工作中獲得更大遺傳進展的潛力更大[22].2個SNP位點聯(lián)合產(chǎn)生3種單倍型和4種雙倍型，理論上應(yīng)有4種單倍型和10種雙倍型，實際觀察值與理論值存在差異，可能由長期人工選育、樣本量不夠或者根本不存在另外一種單倍型引起的.基因在配子中的隨機分離及在合子中的隨機重組都會導(dǎo)致一定的誤差，引起基因頻率的變化，即隨機漂變引起基因頻率偏離Hardy-Weinberg平衡[23].本研究中,g.22589443 C>T位點的基因型分布在長順綠殼蛋雞中極顯著偏離了Hardy-Weinberg平衡(P<0.01)，在黔畫烏雞和威寧雞中顯著偏離了Hardy-Weinberg平衡(P<0.05)；而g.22589512 C>T位點的基因型分布在3個雞種群中均未偏離Hardy-Weinberg平衡(P>0.05).表明g.22589443 C>T位點可能受到遺傳漂變、突變或者人工選育等因素的影響[24].r2和D′值反映了連鎖不平衡的不同方面，D′值僅包括重組史，而r2值包括了重組史和突變史以及諸如自然或人工選擇、基因轉(zhuǎn)換、遺傳漂變等可以影響連鎖不平衡的歷史事件，D′值不適宜小樣本的研究應(yīng)用，r2值小于0.33時，位點之間不存在強連鎖不平衡[25]，表明該基因2個SNP位點趨向獨立遺傳.基因突變使IL8RB基因mRNA最小自由能降低，穩(wěn)定性升高，2個SNP位點雖然引起氨基酸變化，但蛋白質(zhì)的二級結(jié)構(gòu)和三級結(jié)構(gòu)并未改變，蛋白質(zhì)的功能與其空間結(jié)構(gòu)有密切關(guān)系，因此，突變沒有引起蛋白質(zhì)功能和空間結(jié)構(gòu)發(fā)生改變.對IL8RB基因變異位點與黔畫烏雞肉品質(zhì)關(guān)聯(lián)性的分析發(fā)現(xiàn)，H1H1(CCCC)個體的肉色指標(biāo)顯著高于H1H3(CCCT)個體(P<0.05)，且雙倍型H1H1各項指標(biāo)均在正常范圍，雙倍型H1H1是改善肉色的有利雙倍型.今后可以加大對IL8RB基因SNP的研究并擴大篩查范圍，對其多態(tài)性進行更深層次的探討與分析，為雞的保種選育提供參考.