黃勛和 趙雪敏 柯振華 陳龍星 李紀泰 段巖麗 鐘福生

摘要：應用微衛(wèi)星位點LEI0258研究我國烏骨雞的遺傳多態(tài)性與分子進化。以8種烏骨雞和4種黑羽雞為研究材料，通過基因分型與測序，構建中介網絡圖，分析烏骨雞的遺傳多樣性并探討其起源。240份樣品共檢測到67個等位基因，片段長度為181～507 bp。烏骨雞保持著與黑羽雞相當?shù)倪z傳變異水平，觀察雜合度、期望雜合度和多態(tài)信息含量分別為0.813、0.962、0.957。經測序鑒定出38個等位基因，部分等位基因為個別烏骨雞品種所特有。10個等位基因與25種已知MHC血清型相匹配?；趥纫韰^(qū)變異信息的中介網絡圖將38個等位基因分為6條進化枝，未發(fā)現(xiàn)烏骨雞特有的進化枝。研究結果為烏骨雞的保種選育與開發(fā)利用及起源進化研究提供了理論基礎。

關鍵詞：烏骨雞;微衛(wèi)星位點;LEI0258;多態(tài)性;進化

中圖分類號： S831.2? 文獻標志碼： A? 文章編號：1002-1302（2019）05-0040-06

收稿日期：2017-11-07

基金項目：廣東省自然科學基金項目（編號：2014A030307018）;廣東省公益研究與能力建設項目（編號：2015A020208020、2016A030303068）;嘉應學院省市共建重點建設項目[編號：嘉院（2017）27號]。

作者簡介：黃勛和（1982—），男，廣東河源人，博士，副教授，主要從事家雞遺傳資源與進化研究。E-mail：hxh826@jyu.edu.cn。

通信作者：鐘福生，博士，教授，主要從事動物生產與畜牧工程。E-mail：zfs@jyu.edu.cn。

微衛(wèi)星位點LEI0258是近年來逐漸應用于家雞遺傳多樣性與分子進化的遺傳標記，與常規(guī)的微衛(wèi)星不同，LEI0258由R13（CTATGTCTTCTTT）和R12（CTTTCCTTCTTT）2個重復單元組成，同時側翼區(qū)轉換、顛換或缺失等變異較多，因而該位點具有較高的多態(tài)性[1-6]。微衛(wèi)星位點LEI0258定位于雞16號染色體的主要組織相容性復合體B區(qū)域（major histocompatibility complex B region，MHC-B），其遺傳多態(tài)性通常與宿主的抗病能力呈正相關[7-8]，并且與MHC-B血清型有較好對應關系[9]。因此，微衛(wèi)星位點LEI0258對于評估保育群體的遺傳多樣性水平與保護潛力具有一定的參考價值。

烏骨雞是獨具中國特色的畜禽品種之一，因其特有的藥用價值和我國的飲食文化，在近千年的環(huán)境選擇及人工培養(yǎng)下，形成了近20個體型外貌、生產性能各異的地方品種[10]。微衛(wèi)星與線粒體DNA研究表明烏骨雞具有較高的遺傳多樣性水平[11-13]，而烏骨雞LEI0258多態(tài)性的研究尚屬空白。本研究在應用微衛(wèi)星位點LEI0258對我國主要烏骨雞品種進行遺傳多態(tài)性與分子進化分析，評估保育雞種的遺傳變異水平，探討烏骨雞的遺傳起源，為烏骨雞的保種選育與利用提供科學理論依據(jù)。

1 材料與方法

1.1 材料與試劑

烏骨雞（烏雞）及黑羽雞（非烏骨雞）樣品采自各自保種場，每個品種20羽（表1），翅下靜脈取血，終濃度70%的乙醇于-80 ℃保存，采用傳統(tǒng)的酚氯仿法提取基因組DNA，PCR相關試劑購自寶生物工程（大連）有限公司，引物由廣州艾基生物技術有限公司合成。

1.2? PCR擴增與基因分型

擴增引物序列參照McConnell等的序列[1]，正向引物5′端添加FAM熒光標記。PCR擴增體系為30 μL：3.0 μL 10×PCR buffer，2.4 μL dNTP mixture（2.5 mmol/L），上、下游引物（20.0 μmol/L）各0.3 μL，0.3 μL rTaq酶（5 U/μL），50～100 ng 模板DNA，滅菌ddH2O補齊至30 μL。PCR擴增條件：94 ℃預變性4 min;94 ℃ 30 s，63 ℃ 1 min，72 ℃延伸 1 min，38個循環(huán);72 ℃延伸7 min。經1.5%瓊脂糖凝膠電泳檢測后送上海翼和應用生物技術有限公司ABI 3730平臺進行基因分型。

1.3 等位基因測序

根據(jù)基因分型的結果，參照黃勛和等的方法[5]挑選測序樣品，送廣州艾基生物技術有限公司進行雙向測序。

1.4 統(tǒng)計分析

觀察雜合度（observed heterozygosity，Ho）、期望雜合度（expected heterozygosity，HE）、多態(tài)信息含量（polymorphic information content，PIC）、等位基因數(shù)（NA）和等位基因頻率等由軟件Cervus 3.0.3[14]統(tǒng)計。為與黃羽雞[5]進行比較，應用軟件ADZE 1.0[15]對等位基因多態(tài)性進行稀釋處理。使用軟件MEGA 6.0[16]執(zhí)行序列比對。去掉重復序列后，使用軟件SplitsTree 4.10[17]構建基于側翼序列的中介網絡圖（median-joining network），進化枝的命名參考Chazara等[2]和黃勛和等[5-6]的方法。

2 結果與分析

2.1 基因型多態(tài)性

240個樣品基因分型共檢測出67個等位基因，片段大小范圍為181～507 bp，頻率最高的是等位基因205，為0.104 2，其次是193（0.068 8）和310（0.064 6）（表2）。烏骨雞觀察雜合度、期望雜合度和多態(tài)信息含量分別為0.813、0.962和0.957，與黑羽雞和黃羽雞相當，表明中國烏骨雞保持著較高的遺傳變異水平（表3）。其中，德化黑雞的多態(tài)信息含量最高（0.929），余干烏骨雞最低（0.809）。德化黑雞的等位基因豐度最高（14.389），黃羽烏雞最低（9.554）;江山烏骨雞的私有等位基因豐度最高（4.957），黃羽烏雞最低（1.014）。烏骨雞的等位基因豐度與私有等位基因豐度均高于黑羽雞和黃羽雞（表3）。

2.2 核苷酸序列分析

抽檢67個樣品進行雙向測序，共檢測出38個等位基因，片段大小在193～501 bp，其中10個為首次發(fā)現(xiàn)（表4）。3個等位基因與毛細管電泳分型結果一致，分別是193、205和237。微衛(wèi)星位點LEI0258由2個重復單元R13（CTATGTCTTCTTT）和R12（CTTTCCTTCTTT）組成，重復次數(shù)分別1～17次和3～28次。R13上游有55 bp（不含引物序列），含有8個變異位點，以轉換為主，但在-29～-30出現(xiàn)了“TT”缺失的現(xiàn)象。R12下游有53 bp（不含引物序列），含有5個變異位點，以顛換為主，5堿基位置出現(xiàn)了T/C轉換的現(xiàn)象。部分等位基因出現(xiàn)下游11～18位置“ATTTTGAG”缺失的現(xiàn)象，但僅限于小于241 bp的等位基因。等位基因261在11～18突變?yōu)锳TTTTGGG，不同于Fulton等報道的“ATTTGAG”[9]，但與Han等[3]、黃勛和等[5-6]的報道相同。

2.3 側翼區(qū)中介網絡圖分析

去掉R13、R12重復序列后，構建側翼區(qū)變異信息的中介網絡圖（圖1）。38個等位基因可分為6條進化枝（A-F），每條進化枝有1～15個等位基因。其中，A枝的等位基因最多，為15個;E枝只含有1個等位基因。有些等位基因只出現(xiàn)在特定的烏骨雞品種，如419（德化黑雞）、421（金湖烏鳳雞）;而有些等位基因則分布在多個烏骨雞品種，如193、205、249和333。每條進化枝均有黃羽雞與烏骨雞共享的等位基因，其中A、B、C同樣存在于紅原雞中[2，9]。10個等位基因的血清型分別與25種MHC單倍型相對應，分布在A、B、D和E進化枝中，其中2個等位基因為黑羽雞特有（表4）。

3 討論與結論

烏骨雞是我國寶貴的家禽遺傳資源，具有重要的經濟價值和科學研究價值。微衛(wèi)星位點是家雞品種保育效果動態(tài)監(jiān)測的重要遺傳標記[18-19]。微衛(wèi)星LEI0258遺傳變異分析顯示，烏骨雞保持著與黑羽雞相當?shù)倪z傳多樣性水平（表3），表明烏骨雞保育狀態(tài)較好，有利于后續(xù)的保種選育與開發(fā)利用。從期望雜合度、觀察雜合度和多態(tài)信息含量可以看出，與常規(guī)微衛(wèi)星相比，微衛(wèi)星LEI0258位點普遍具有較高的遺傳變異性[5-6]，這可能與所選研究樣本以及遺傳標記不同有關。

微衛(wèi)星位點LEI0258特異等位基因可作為品種鑒定的參考依據(jù)[3，5-6]。烏骨雞具有較高的營養(yǎng)價值，為提高生長生產性能和營養(yǎng)價值，各地開展了烏骨雞品種保種選育和新品系培養(yǎng)等工作。但市面上出現(xiàn)了許多良莠不齊的商品雞冒充地方烏骨雞品種，給鑒定工作帶來了極大的困擾。本研究發(fā)現(xiàn)，少數(shù)等位基因為個別品種所特有，如260（黃山黑雞）、308（金湖烏鳳雞）。將來還需加大研究樣品量以鑒定更多的品種特異等位基因，同時聯(lián)合其他分子標記[20-23]以及外貌特征開展品種鑒定檢測工作。

烏骨雞具有較強的環(huán)境適應能力，在我國大部分省份均有分布。MHC遺傳變異性通常與家禽抗病能力密切相關[24-26]，而微衛(wèi)星位點LEI0258等位基因與MHC單倍型有較強的相關性[9]。本研究根據(jù)已有材料整理了烏骨雞微衛(wèi)星位點LEI0258等位基因血清型與MHC單倍型的對應關系，后續(xù)還需開展烏骨雞MHC單倍型相關的試驗鑒定工作。

家雞馴化的起源地和時間迄今仍存在較大爭議[27-33]，基于側翼區(qū)的中介網絡圖反映了等位基因之間的關系[34]，對系統(tǒng)進化研究亦有一定的參考價值[5-6]。在烏骨雞所有的進化枝中，均存在與黑羽雞和黃羽雞共享的等位基因（圖1），烏骨雞難以從家雞中完全分離出來，提示烏骨雞可能是從家雞選育而來，但具體從哪些品種選育而來、具體選育時間等信息則需進一步研究。紅原雞也有部分等位基因分布于進化枝A、B和C中[2，9]，表明紅原雞是這些家雞的共同祖先。將來需要加大烏骨雞品種覆蓋度和樣品數(shù)量，從黑色素等特色基因全基因組角度探討烏骨雞的起源與擴散歷史。

致謝：感謝華南農業(yè)大學張細權教授與何丹林高級實驗師、安徽農業(yè)大學耿照玉教授、廣西大學夏中生教授與楊秀榮教授、江西省農業(yè)科學院武艷平博士、山東農業(yè)大學李顯耀教授提供樣品，以及廣東海洋大學杜炳旺教授、福建農林大學李昂教授在采集樣品時提供的幫助。

參考文獻：

[1]McConnell S K，Dawson D A，Wardle A，et al. The isolation and mapping of 19 tetranucleotide microsatellite markers in the chicken[J]. Animal Genetics，1999，30（3）：183-189.

[2]Chazara O，Chang C S，Bruneau N A，et al. Diversity and evolution of the highly polymorphic tandem repeat LEI0258 in the chicken MHC-B region[J]. Immunogenetics，2013，65（6）：447-459.

[3]Han B，Lian L，Qu L J，et al. Abundant polymorphisms at the microsatellite locus LEI0258 in indigenous chickens[J]. Poultry Science，2013，92（12）：3113-3119.

[4]金麗娜，韓建林. 雞MHC區(qū)域內微衛(wèi)星LEI0258和MCW0312的遺傳多態(tài)性[J]. 畜牧與獸醫(yī)，2015，47（7）：50-55.

[5]黃勛和，李麗芝，張金楓，等. 華南家雞MHC-B區(qū)域復合微衛(wèi)星位點LEI0258的遺傳多樣性與進化研究[J]. 畜牧獸醫(yī)學報，2016，47（11）：2175-2183.

[6]黃勛和，張金楓，譚坤鳳，等. 五華三黃雞MHC-B區(qū)域復合微衛(wèi)星位點LEI0258遺傳多樣性與進化研究[J]. 廣東農業(yè)科學，2016，43（6）：163-168.

[7]Bacon L D，Hunt H D，Cheng H H. Genetic resistance to Mareks disease[J]. Current Top Microbiology Immunogenetics，2001，255：121-141.

[8]Lee L F，Bacon L D，Yoshida S，et al. The efficacy of recombinant fowlpox vaccine protection against Mareks disease：Its dependence on chicken line and B haplotype[J]. Avian Diseases，2004，48（1）：129-137.

[9]Fulton J E，Juul-Madsen H R，Ashwell C M，et al. Molecular genotype identification of the Gallus gallus major histocompatibility complex[J]. Immunogenetics，2006，58（5/6）：407-421.

[10]國家畜禽遺傳資源委員會. 中國畜禽遺傳資源志：家禽志[M]. 北京：中國農業(yè)出版社，2011.

[11]湯青萍，陳寬維，李慧芳，等. 應用微衛(wèi)星標記對12個中國地方烏骨雞品種遺傳多樣性的研究[J]. 畜牧獸醫(yī)學報，2005，36（8）：755-760.

[12]魏 麟，劉勝貴，史憲偉. 雪峰烏骨雞自然群體遺傳多樣性的微衛(wèi)星分析[J]. 生物多樣性，2008，16（5）：503-508.

[13]徐文娟，朱文奇，束婧婷，等. 我國主要烏骨雞品種遺傳多樣性和系統(tǒng)進化研究[J]. 中國畜牧雜志，2014，50（23）：10-14.

[14]Kalinowski S T，Taper M L，Marshall T C. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment[J]. Molecular Ecology，2007，16（5）：1099-1106.

[15]Szpiech Z A，Jakobsson M，Rosenberg N A. ADZE：a rarefaction approach for counting alleles private to combinations of populations[J]. Bioinformatics，2008，24（21）：2498-2504.

[16]Tamura K，Stecher G，Peterson D，et al. MEGA6：molecular evolutionary genetics analysis version 6.0[J]. Molecular Biology and Evolution，2013，30（12）：2725-2729.

[17]Huson D H，Bryant D. Application of phylogenetic networks in evolutionary studies[J]. Molecular Biology and Evolution，2006，23（2）：254-267.

[18]Groeneveld L F，Lenstra J A，Eding H，et al. Genetic diversity in farm animals - a review[J]. Animal Genetics，2010，41（1）：6-31.

[19]張跟喜，丁馥香，王金玉，等. 利用微衛(wèi)星標記分析邊雞遺傳多樣性及保種效果[J]. 農業(yè)生物技術學報，2010，18（5）：944-950.

[20]Yacoub H A，F(xiàn)athi M M，Sadek M A. Using cytochrome b gene of mtDNA as a DNA barcoding marker in chicken strains[J]. Mitochondrial DNA，2015，26（2）：217-223.

[21]黃勛和，陳潔波，何丹林，等. DNA條形碼技術鑒定中國地方雞品種的重新評估[J]. 中國農業(yè)科學，2016，49（13）：2622-2633.

[22]秦玉梅，任 嵩，李佳玉，等. 雞FSHβ、ESRα基因多態(tài)性及其合并基因型與產蛋性能的關聯(lián)性分析[J]. 江蘇農業(yè)學報，2017，33（4）：854-862.

[23]陶志云，徐文娟，朱春紅，等. 高郵鴨GH基因SNP位點及其與早期體質量的相關性[J]. 江蘇農業(yè)學報，2017，33（1）：146-150.

[24]Rogers S L，Kaufman J. High allelic polymorphism，moderate sequence diversity and diversifying selection for B-NK but not B-lec，the pair of lectin-like receptor genes in the chicken MHC[J]. Immunogenetics，2008，60（8）：461-475.

[25]李福偉，逯 巖，雷秋霞，等. 山東地方雞MHCB-F基因遺傳變異與免疫性狀相關性研究[J]. 農業(yè)生物技術學報，2013，21（2）：179-184.

[26]Wang H Z，Ma T，Chang G B，et al. Molecular genotype identification of different chickens：major histocompatibility complex[J]. Journal of Science and Technology，2014，2：1-7.

[27]Liu Y P，Wu G S，Yao Y G，et al. Multiple maternal origins of chickens：out of the Asian jungles[J]. Molecular Phylogenetics and Evolution，2006，38（1）：12-19.

[28]Miao Y W，Peng M S，Wu G S，et al. Chicken domestication：an updated perspective based on mitochondrial genomes[J]. Heredity，2013，110（3）：277-282.

[29]Xiang H，Gao J Q，Yu B Q，et al. Early holocene chicken domestication in Northern China[J]. Proceedings of the National Academy of Sciences of the United States of America，2014，111（49）：17564-17569.

[30]Peng M S，Shi N N，Yao Y G，et al. Caveats about interpretation of ancient chicken mtDNAs from Northern China[J]. Proceedings of the National Academy of Sciences of the United States of America，2015，112（16）：E1970-E1971.

[31]Peters J，Lebrasseur O，Best J，et al. Questioning new answers regarding Holocene chicken domestication in China[J]. Proceedings of the National Academy of Sciences of the United States of America，2015，112（19）：E2415.

[32]Eda M，Lu P，Kikuchi H，et al. Reevaluation of early Holocene chicken domestication in Northern China[J]. Journal of Archaeological Science，2016，67：25-31.

[33]Peters J，Lebrasseur O，Deng H，et al. Holocene cultural history of red jungle fowl （Gallus gallus） and its domestic descendant in East Asia[J]. Quaternary Science Reviews，2016，142：102-119.

[34]E G X，Sha R，Zeng S C，et al. Genetic variability，evidence of potential recombinational event and selection of LEI0258 in chicken[J]. Gene，2014，537（1）：126-131.劉雪嬌，賈田惠，高同國，等. 貝萊斯芽孢桿菌3A3-15電擊轉化條件優(yōu)化[J]. 江蘇農業(yè)科學，2019，47（5）：46-49.