驢遺傳育種現(xiàn)狀與發(fā)展趨勢

侯浩賓，李海靜，張 莉*

（1.中國農(nóng)業(yè)科學院北京畜牧獸醫(yī)研究所，北京 100193；

2.國家膠類中藥工程技術(shù)研究中心，東阿阿膠股份有限公司，山東 東阿 252201）

驢（Equus asinus L.）是重要的役用家畜，在世界農(nóng)耕文明中扮演了重要角色。據(jù)史料記載，中國北方和西部少數(shù)民族在3 500年前已經(jīng)馴養(yǎng)了驢，并通過馬和驢雜交獲得了騾，由于數(shù)量稀少，異常珍貴，當時被稱為奇畜[1]。隨著農(nóng)業(yè)機械和交通運輸工具的發(fā)展普及，驢的役用功能逐漸降低，而肉用、乳用和藥用等功能逐漸增加。分子遺傳學的發(fā)展，使得動物遺傳育種由常規(guī)的表型選擇發(fā)展到現(xiàn)在的標記輔助選擇，大大提高了動物遺傳改良的進程，尤其是基因組學的發(fā)展為挖掘驢重要經(jīng)濟性狀的功能基因提供了可能。驢遺傳育種的研究主要集中在生產(chǎn)性能、mtDNA序列多態(tài)性、微衛(wèi)星標記和功能基因等方面；繁育技術(shù)的研究主要集中在同期發(fā)情、人工授精和體外胚胎生產(chǎn)等方面。本文對驢遺傳育種研究進行了綜述，并對未來驢遺傳育種研究進行展望。

1 我國驢遺傳資源概況

1.1 驢的類型與分布

我國是養(yǎng)驢大國，驢品種資源豐富，其中24個具有地方特色的品種收錄于《中國畜禽遺傳資源志-馬驢駝志》[2]。按其生態(tài)類型可將驢分為干旱沙漠型（新疆驢）、干旱半荒漠型（涼州驢）、黃土丘陵溝壑型（佳米驢）、高原山地型（西南驢）、中部平原丘陵型（德州驢）[3]。按體型特征可分為大型驢、中型驢和小型驢三種類型[4]。大型驢體高超過130 cm，體重在260 kg左右，主要包括平原地區(qū)（山東、陜西、山西、河北等）的德州驢、關(guān)中驢、廣靈驢和晉南驢等；中型驢體高在110～130 cm，體重180 kg左右，主要包括高原地區(qū)（陜西、甘肅、山西、河北）與平原（河南）地區(qū)的慶陽驢、淮陽驢、佳米驢和沁陽驢等；小型驢體高低于110 cm，體重在130 kg左右，其分布范圍較廣，主要分布在新疆、內(nèi)蒙古、寧夏、甘肅、陜西以及云南、四川、西藏等省、自治區(qū)，包括新疆驢、云南驢、川驢和青海毛驢等。

1.2 驢的存欄量

2016年我國驢存欄量456.9萬頭，比2015年下降了15.7%，創(chuàng)歷史新低。從1990年以來，我國驢的存欄量連續(xù)26年呈下降趨勢，平均每年約下降25.5萬頭[5]。驢存欄量連續(xù)下滑的原因主要是役用功能被農(nóng)業(yè)機械和現(xiàn)代交通工具所取代以及驢本身的繁殖效率低、繁育速度慢的生物學特性[6]。

2 驢遺傳育種研究進展

通過對驢遺傳育種發(fā)表論文數(shù)量的統(tǒng)計，排在前五位的國家由高到低排序分別是意大利、中國、巴西、西班牙和美國，其中意大利、中國和巴西近十年來有所增加[7]。意大利驢存欄量雖比中國少得多，也不及巴西，但意大利的卡梅里諾大學（Univ Camerino）、巴里大學（Univ Bari）和卡塔尼亞大學（Univ Catania）均在驢上進行了深入研究，涉及驢肉品質(zhì)與胴體性狀[8]、理化特性[9]，驢奶蛋白組分[10]、及驢種群多樣性研究[11]等。國內(nèi)對于驢研究主要集中在不同品種的遺傳多樣性及起源[12]、品種資源的保護與利用[13]、體尺與體重的相關(guān)性[14]、驢肉營養(yǎng)成分[15]等方面。

2.1 驢生產(chǎn)性能研究

2.1.1 產(chǎn)肉性能

驢的品種不同，其屠宰性能、肉質(zhì)理化特征及加工特性有很大差異。在品種方面，德州驢、廣靈驢、泌陽驢和疆岳驢可以向肉用方向開發(fā)；關(guān)中驢和廣靈驢肌間脂肪含量較高，可以向高檔肉方向開發(fā)；而云南驢產(chǎn)肉率較低，不適合作為肉用驢開發(fā)[16]。在年齡方面，德州驢的屠宰率、凈肉率和胴體產(chǎn)肉率均隨年齡的增加而逐漸下降，2.5～3.5歲改良德州母驢的屠宰性能最佳，食用品質(zhì)好，為最適宜屠宰和食用的年齡段[17]。意大利地區(qū)的Martina Franca驢在不同屠宰年齡以及不同哺乳方式下驢肉品質(zhì)具有明顯差異，結(jié)果表明18月齡的自然哺乳驢要比12月齡的驢具有更高的肌間脂肪和蛋白質(zhì)含量[18]。研究者還對不同年齡育肥驢的驢肉進行了營養(yǎng)成分分析，發(fā)現(xiàn)隨年齡增加，驢肉中干物質(zhì)、脂肪、能量上升，粗蛋白、灰分下降[19]。在育肥方面，斷奶驢不宜作為肥育肉用；成年退役驢以63 d強度肥育較好；1.5歲驢可短期（63 d，83 d）強度肥育生產(chǎn)優(yōu)質(zhì)驢肉，或肥育數(shù)月生產(chǎn)高中檔驢肉[20]。不同性別的德州驢其生產(chǎn)性能和屠宰性能存在顯著差異，母驢比公驢日增重高9.44%，肥育性也要高于公驢，而出皮率低于公驢[21]。在雜交改良方面，利用新疆驢和關(guān)中驢雜交改良，可以提高驢的產(chǎn)肉性能[22]。在體型選擇方面，體重和體尺相關(guān)性研究表明疆岳驢的胸圍對體重的影響最大，體高次之，體長最??；疆岳驢在選育、提純復(fù)壯和保種時應(yīng)加大對胸圍、體高和體長的選擇力[23]。

2.1.2 產(chǎn)奶性能

與牛、羊相比，驢的產(chǎn)奶特性非常特殊。國外報道了驢奶主要成分的晝夜節(jié)律性，觀察到體細胞數(shù)和pH沒有晝夜節(jié)律性，乳脂和乳糖含量在夜間達到高峰，蛋白質(zhì)含量在白天達到高峰[24]。國內(nèi)的疆岳驢在一個泌乳期內(nèi)，最高擠乳量可達893 kg，具有選育為乳用型驢品種的潛力[25]。意大利Ragusana驢的日產(chǎn)奶量、泌乳周期和驢奶的特性變化均取決于產(chǎn)駒季節(jié)，春天驢的產(chǎn)奶量和質(zhì)量都是最好的；在泌乳期，驢奶產(chǎn)量持續(xù)下降，蛋白質(zhì)和脂肪含量也呈下降趨勢，而乳糖含量呈上升趨勢[26]。

2.2 驢遺傳多樣性研究

2.2.1 線粒體DNA（mtDNA）研究

細胞色素b（Cytb）是線粒體的一個氧化磷酸化系統(tǒng)III蛋白復(fù)合體組成部分，由線粒體基因組Cytb基因編碼。mtDNA序列多態(tài)性分析揭示了埃塞俄比亞驢的遺傳多樣性和母系遺傳特征，結(jié)果表明埃塞俄比亞可能是非洲家驢多樣性的中心之一[27]。國內(nèi)首次對德州驢的全線粒體基因組進行了測序和組裝，系統(tǒng)發(fā)育分析表明德州驢可能起源于非洲[28]。在我國四個家驢品種群體中發(fā)現(xiàn)了豐富的mtDNA多樣性和兩個截然不同的線粒體母性世系，即索馬里和努比亞種間mtDNA D-loop區(qū)的遺傳變異及遺傳指數(shù)分析表明該群體具有較高的群體遺傳多樣世系；家驢Cytb基因的遺傳變異格局表明這兩個世系沒有發(fā)生過群體擴張事件[29]。意大利六個驢品種間mtDNA D-loop區(qū)的遺傳變異及遺傳指數(shù)分析表明該群體具有較高的群體遺傳多樣性[30]。中國四個地方驢品種mtDNA D-Loop區(qū)部分序列分析揭示了我國的家驢可能起源于非洲野驢的一個分支—索馬里野驢[31]。同時，我國4個家驢品種21個個體的mtDNA D-loop區(qū)399 bp序列的研究表明我國的驢遺傳多樣性正逐步喪失，需要加強驢種質(zhì)資源保護[32]。

2.2.2 微衛(wèi)星標記研究

國內(nèi)外對利用微衛(wèi)星標記揭示驢不同品種間的遺傳變異性進行了諸多報道。國外學者用13個馬微衛(wèi)星位點對5個瀕危的西班牙驢品種遺傳多樣性進行了比較，根據(jù)DA距離矩陣構(gòu)建樹狀圖表現(xiàn)出西班牙驢品種間的差異不大[33,34]。Matassino等[35]利用微衛(wèi)星分子標記評價了Amiata驢和拉齊奧區(qū)（意大利行政區(qū)名）Viterbese驢種群的遺傳多樣性以及相似性。研究表明Viterbese驢遺傳變異的丟失主要是來源于被忽視的亞群體結(jié)構(gòu)（Wahlund效應(yīng)），而不是近親繁殖、無效等位基因或者選擇的影響。國內(nèi)學者利用24對微衛(wèi)星標記對中國8個大、中型驢品種遺傳多樣性進行了檢測分析，采用鄰接法構(gòu)建系統(tǒng)發(fā)生樹，并進行了系統(tǒng)發(fā)生分析，結(jié)果表明該群體基因多態(tài)性和遺傳多樣性相對較高；其中大型品種關(guān)中驢、晉南驢、廣靈驢、德州驢和中型品種慶陽驢聚為一類，中型品種泌陽驢、淮陽驢和佳米驢聚為另一類，各驢種的分子系統(tǒng)發(fā)生關(guān)系與其育成史和地理分布基本一致[36]。新疆3個地方驢群體的微衛(wèi)星標記研究表明該群體均表現(xiàn)出高度多態(tài)性[37]。烏頭驢和三粉驢是德州驢的兩個品系，通過檢測兩個品系13對微衛(wèi)星位點揭示了兩個驢品系的遺傳雜合度均偏低，遺傳變異程度較小，其中三粉驢的遺傳多態(tài)性高于烏頭驢[38]。此外，我國5個驢品種（關(guān)中驢，德州驢，華北驢，晉南驢，廣靈驢）的微衛(wèi)星標記分析還揭示了這5個品種的遺傳多樣性相對較高，但廣靈驢和晉南驢出現(xiàn)雜合子赤字，需要加強遺傳資源的保護[39]。

2.3 驢重要經(jīng)濟性狀功能基因研究

2.3.1 生長和皮用性狀相關(guān)基因

Myostatin（MSTN）基因?qū)趋兰∩L發(fā)育具有負調(diào)控作用，該基因的變異對調(diào)控哺乳動物骨骼肌生長起重要作用。LIU等[40]分析了中國13個驢品種MSTN基因的單核苷酸多態(tài)性（SNPs），對4個新的單核苷酸多態(tài)性進行了檢測和基因分型，構(gòu)建了不同物種間MSTN基因進化樹，9個聚類結(jié)果與物種分化的事實一致。國內(nèi)科研人員還對Cytb（細胞色素b）基因[41]、生長素基因[42]、DGAT2基因[43]與生長性狀進行了關(guān)聯(lián)分析，基因分型結(jié)果表明Cytb基因與云南驢的尻寬和德州驢的體高顯著相關(guān)；生長素基因序列在1267位的C→G可能影響到驢的生長發(fā)育；DGAT2基因存在多態(tài)性，其第3內(nèi)含子有兩種等位基因A和B，存在AA和AB兩種基因型，其中AB型存在堿基發(fā)生A→G的突變，使賴氨酸突變?yōu)榫彼幔梢宰鳛轶H體尺性狀的候選分子遺傳標記。驢皮含有豐富的膠原蛋白，是阿膠的重要原料。新疆驢驢皮中I型膠原蛋白2個亞基基因（COL1A1、COL1A2）mRNA的表達量及其蛋白含量的差異分析表明6～7歲年齡段的驢皮中COL1A1、COL1A2基因的表達量及其蛋白含量顯著高于其他各年齡段，確定了6～7歲年齡段為宰驢取皮的最佳年齡段[44]。

2.3.2 乳用性狀相關(guān)基因

CSN1S2基因編碼產(chǎn)生的αS2酪蛋白是反芻動物和馬科動物奶中三個鈣敏感酪蛋白之一。CSN1S2 cDNA的克隆測序和序列分析，鑒定了CSN1S2 I基因的兩個突變（外顯子12第119 nt T→C；外顯子14第12 nt G→A）和CSN1S2 II基因的一個突變（內(nèi)含子1第197 nt A→G）[45]。PRLR基因是一類重要的參與乳蛋白合成的基因，能夠調(diào)控乳蛋白基因的表達[46]。該基因的突變與疆岳驢泌乳性狀存在關(guān)聯(lián)性，可以作為疆岳驢乳用型選育的分子遺傳標記[47]。驢β-LG II基因座具有高度多態(tài)性，目前國外學者已經(jīng)發(fā)現(xiàn)了 6 個變異分別命名為 A[48]、B、C[49]、D[50]、E[51]和 F[52]，其中 F 位點與驢乳 β-LG II蛋白不表達或低表達有關(guān)。

2.3.3 毛色相關(guān)基因

KIT基因的變異在馬中能夠引起一系列白色表型，可以選擇KIT基因作為影響驢白斑和白色兩個毛色表型的候選基因。研究人員通過對KIT基因21個外顯子的突變分析發(fā)現(xiàn)白色驢存在錯義突變（c.662A＞C；p.Tyr221Ser）；而白斑驢存在第二種突變 （c.1978+2T＞A）[53]。許多動物的白化癥是由酪氨酸酶（tyrosinase，TYR）基因的突變引起的，酪氨酸酶是合成真黑素和棕黑素兩種黑色素的關(guān)鍵酶。意大利的Asinara白化驢TYR基因突變分析發(fā)現(xiàn)了一個高度保守的氨基酸位置有一個錯義突變 （c.604C＞G；p.His202Asp）與該驢的白化表型顯著相關(guān)[54]。另外，還有研究發(fā)現(xiàn)ASIP基因的一個突變（c.34902T65＞65C）與法國驢品種的黑色表型相關(guān)[55]。

3 驢繁殖技術(shù)研究

3.1 同期發(fā)情

同期發(fā)情在很多家畜和家禽中已經(jīng)得到了很好的應(yīng)用，目前該技術(shù)在驢上的應(yīng)用和研究尚未成熟。何良軍等[56]運用PGF2α（前列腺素F2α）和氯前列烯醇處理,誘導(dǎo)驢同期發(fā)情,對比兩種藥物的同期發(fā)情發(fā)現(xiàn)PGF2α比氯前列烯醇對驢同期發(fā)情的作用效果好。而人為給同一頭母驢注射PMSG（孕馬血清促性腺激素）與自然發(fā)情相比,驢的發(fā)情和排卵周期無顯著差異，并且PMSG處理不能使母驢超數(shù)排卵[57]。

3.2 人工授精技術(shù)

驢的繁殖周期很長，且為單胎動物，群體擴繁慢。運用人工授精技術(shù)可以提高優(yōu)秀種公驢的利用價值、提高母驢繁殖力，同時還有助于恢復(fù)驢瀕危品種的群體。首先，要選擇體型高大勻稱，睪丸發(fā)育良好的種公驢；對于選好的種公驢可全年進行精液的采集和冷凍保存，而冬季為發(fā)情淡季，適當調(diào)教種公驢后才可采精；隔天采一次或每周采三次是最理想的采精頻率[58]。國內(nèi)科研人員通過對5種冷凍保存液的精液保存效果比較發(fā)現(xiàn)：乳糖-葡糖-甘油-卵黃液為最佳種冷凍保存液（P＜0.05）；甘油為最佳冷凍保護劑,其適宜濃度為5%，精子凍后活率平均為0.36，畸形率≤27.4%，頂體完整率82.2%，解凍后常溫下精子存活時間為8.2 h[59]。同時，國外有研究表明利用甘油或者乙二醇作為冷凍保護液對精液活力和受孕率無顯著影響[60]。國外還報道了自動化冷凍精液的方法相比傳統(tǒng)的方法，精液表現(xiàn)出較高的活力；運用驢新鮮精液增加每一次劑量的精子細胞數(shù)以及增加人工授精次數(shù)能夠提高母驢的受孕率[61]。直腸引導(dǎo)深角授精技術(shù)（TRG）和宮腔鏡深角授精技術(shù)（HYS）是目前人工授精最常用的2種方式[62]，研究表明TRG和HYS兩種授精技術(shù)對驢受孕率無顯著影響[63,64]。同時，Hayden等[65]發(fā)現(xiàn)當輸精量足夠時，利用這兩種授精技術(shù)處理的母驢受孕率無顯著差異。

3.3 體外胚胎生產(chǎn)

驢的胚胎體外生產(chǎn)技術(shù)遠遠落后于其他家畜，主要是由于對驢排卵機制的研究基礎(chǔ)薄弱，而且驢活體采卵、卵母細胞體外保存與培養(yǎng)、體外受精等技術(shù)的應(yīng)用尚處在起步階段。驢卵泡直徑大約每天增加1.5～3 mm，優(yōu)勢卵泡平均每天增長2.1 mm，發(fā)情周期在20～24 d之間，發(fā)情一般持續(xù)5～7 d，約在發(fā)情停止前1 d排卵[66]。驢卵母細胞的體外成熟（IVM）比馬卵母細胞的IVM慢，并且驢卵母細胞的IVM的最佳時間可能是34 h[67]。杜興雨等[68]通過分析不同濃度冷凍保護液及冷凍載體對驢卵母細胞成熟率的影響，可以篩選出較優(yōu)的冷凍保護液，研究表明開放式拉長塑料細管作為冷凍載體的效果顯著優(yōu)于普通玻璃細管。這些條件為建立驢的體外受精（IVF）技術(shù)及體外胚胎生產(chǎn)技術(shù)提供了前提條件。

4 未來發(fā)展趨勢

近年來，隨著分子遺傳育種的發(fā)展，家畜重要經(jīng)濟性狀形成的遺傳機制不斷被揭示。在功能基因研究方面，篩選和鑒定了諸多調(diào)控家畜繁殖、抗病、產(chǎn)肉、產(chǎn)奶等重要經(jīng)濟性狀的功能基因，部分重要經(jīng)濟性狀的QTL位點被定位。目前，我國驢品種的培育主要是通過傳統(tǒng)的育種方式實現(xiàn)的。如何將常規(guī)育種與分子育種結(jié)合起來培育新品種是驢遺傳育種發(fā)展的重要研究方向。

4.1 驢品種資源評價與保護利用

建立驢品種表型數(shù)據(jù)庫和遺傳資源基因庫，建設(shè)種質(zhì)資源共享平臺。通過對驢品種資源系統(tǒng)化的產(chǎn)肉、產(chǎn)奶、產(chǎn)皮、抗逆性等重要性狀的評價，為保護種質(zhì)資源，培育專門化的驢品種提供理論依據(jù)。一些瀕危的驢品種可能在產(chǎn)肉、產(chǎn)奶等生產(chǎn)性能方面不具有優(yōu)勢，但對于環(huán)境的適應(yīng)性、抗逆性等方面具有重要的科研價值，是驢產(chǎn)業(yè)可持續(xù)發(fā)展的寶貴資源，應(yīng)給予重點保護。在保護的同時也要加強驢品種的選育與利用，而運用現(xiàn)代分子生物學技術(shù)和繁殖技術(shù)可以加快驢遺傳育種的進程。

4.2 基因組測序在驢育種領(lǐng)域的應(yīng)用

伴隨著全基因組測序技術(shù)的發(fā)展，牛、馬、牦牛、家豬、山羊等重要家畜的基因組已被陸續(xù)解析，大量影響重要經(jīng)濟性狀的候選基因被鑒定出來。最近，丹麥科學家在Science上公布了最新的家驢近染色體級別的參考基因組，為未來驢重要經(jīng)濟性狀功能基因的挖掘與鑒定提供了重要條件。利用全基因組關(guān)聯(lián)分析（GWAS）可以定位調(diào)控重要經(jīng)濟性狀的位點，運用功能基因組數(shù)據(jù)分析揭示基因功能、表達調(diào)控機制和調(diào)控網(wǎng)絡(luò)，進而通過定位、分離、克隆出與驢產(chǎn)肉、產(chǎn)奶、產(chǎn)皮、抗病等相關(guān)的功能基因。同時，可以運用在其他動物上已經(jīng)鑒定出來影響重要經(jīng)濟性狀的功能基因作為候選基因，在驢基因組上進行驗證。

參考文獻：

[1]吳錦淑,楊從軍,段玉蘭.淺談中國驢種起源與馴化[J].山東畜牧獸醫(yī),2008,29（3）:16-17.

[2]中國畜禽遺傳資源委員會.中國畜禽遺傳資源志—馬驢駝志[M].北京:中國農(nóng)業(yè)出版社,2011.

[3]孫玉江,徐紀尊,蔣濤,等.中國驢種遺傳資源保護利用研究[J].中國草食動物科學,2006,26（6）:32-34.

[4]謝成俠.中國馬驢品種志[M].上海：上海科學技術(shù)出版社,1987.

[5]農(nóng)業(yè)部畜牧業(yè)司.2017年中國畜牧業(yè)統(tǒng)計[M].北京:中國農(nóng)業(yè)出版社,2017.

[6]張莉,杜立新.對我國驢產(chǎn)業(yè)發(fā)展的思考與建議[J].草食家畜,2015，（5）:1-5.

[7]魏珣,賈敬敦,孫康泰,等.基于文獻計量的世界家畜種業(yè)科技創(chuàng)新研究態(tài)勢分析[J].中國農(nóng)業(yè)科學,2015,48（13）:2622-2634.

[8]Polidori P,Vincenzetti S,Cavallucci C,et al.Quality of donkey meat and carcass characteristics[J].Meat Science,2008,80（4）:1222.

[9]Polidori P,Cavallucci C,Beghelli D,et al.Physical and chemical characteristics of donkey meat from Martina Franca breed[J].Meat Science,2009,82（4）:469-471.

[10]Vincenzetti S,Polidori P,Mariani P,et al.Donkey's milk protein fractions characterization[J].Food Chemistry,2008,106（2）:640-649.

[11]Rizzi R,Tullo E,Cito A M,et al.Monitoring of genetic diversity in the endangered Martina Franca donkey population[J].Journal of Animal Science,2011,89（5）:1304-1311.

[12]王全喜,陳建興,閔令江,等.中國家驢品種遺傳多樣性及母系起源研究[J].內(nèi)蒙古農(nóng)業(yè)大學學報（自然科學版）,2012,33（2）:7-11.

[13]封赟,張金川,江曉軍.關(guān)中驢的保種與開發(fā)利用[J].農(nóng)業(yè)科技與信息,2009，（5）:49-50.

[14]肖國亮,呂長鵬,伊力哈木,等.新疆驢成年母驢體重及體尺性狀的相關(guān)關(guān)系[J].當代畜禽養(yǎng)殖業(yè),2007,28（6）:55-56.

[15]尤娟,羅永康,張巖春.驢肉脂肪和脂肪酸組成特點及與其他畜禽肉的分析比較[J].食品科技,2009,34（2）:118-120.

[16]周楠,韓國才,柴曉峰,等.驢的產(chǎn)肉、理化指標及加工特性比較研究[J].畜牧獸醫(yī)學報,2015,46（12）:2314-2321.

[17]周楠,謝鵬,鄭世學,等.不同年齡改良德州母驢屠宰性能、肉品理化指標及加工特性研究[J].黑龍江畜牧獸醫(yī),2014（9）.

[18]De P P,Tateo A,Maggiolino A,et al.Martina Franca donkey meat quality:Influence of slaughter age and suckling technique[J].Meat Science,2017，（134）:128.

[19]雷亞非,侯文通.不同年齡肥育驢肉的營養(yǎng)成分分析[C].2015中國驢業(yè)發(fā)展大會.2015，1-9.

[20]侯文通.不同年齡和營養(yǎng)水平下驢肥育性能測定[J].草食家畜,2016，（2）:21-29.

[21]劉桂芹,曲洪磊，種肖玉,等.性別對生長驢生產(chǎn)性能、屠宰性能及器官指數(shù)的影響[J].飼料工業(yè),2017，（3）:61-64.

[22]王培基,焦多成,高景輝,等.關(guān)新雜交驢部分體尺和產(chǎn)肉性能測定[J].家畜生態(tài)學報,2007,28（2）:35-36.

[23]肖海霞,托乎提·阿及德,石國慶,等.疆岳驢體重和體尺性狀的相關(guān)和回歸分析[J].中國奶牛,2012（23）:27-31.

[24]Piccione G,Fazio F,Caola G,et al.Daily rhythmicity in nutrient content of asinine milk[J].Livestock Science,2008,116（1-3）:323-327.

[25]肖國亮,周小玲,陳根元,等.疆岳驢泌乳特征和產(chǎn)乳性能分析[J].中國奶牛,2015，（1）:16-20.

[26]Giosuè C,Alabiso M,Russo G,et al.Jennet milk production during the lactation in a Sicilian farming system[J].Animal An International Journal of Animal Bioscience,2008,2（10）:1491-1495.

[27]Kefena E,Dessie T,Tegegne A,et al.Genetic diversity and matrilineal genetic signature of native Ethiopian donkeys （Equus asinus）inferred from mitochondrial DNA sequence polymorphism[J].Livestock Science,2014,167（9）:73-79.

[28]Sun Y,Jiang Q,Yang C,et al.Characterization of complete mitochondrial genome of Dezhou donkey（Equus asinus ）and evolutionary analysis[J].Current Genetics,2016,62（2）:383-390.

[29]Chen J,Sun Y,Manglai D,et al.Maternal genetic diversity and population structure of four Chinese donkey breeds[J].Livestock Science,2010，131（2-3）:272-280.

[30]Cozzi M C,Valiati P,Cherchi R,et al.Mitochondrial DNA genetic diversity in six Italian donkey breeds（Equus asinus）[J].Mitochondrial Dna Part A Dna Mapping Sequencing&Analysis,2017，1.

[31]孫偉麗,楊博輝,曹學亮,等.中國四個地方驢品種mtDNA D-Loop部分序列分析與系統(tǒng)進化研究[J].中國草食動物科學,2007,27（2）:7-10.

[32]雷初朝.中國四個畜種（黃牛、水牛、牦牛、家驢）線粒體DNA遺傳多樣性研究[D].楊凌：西北農(nóng)林科技大學，2002.

[33]Arangurenméndez J,Jordana J,Gomez M.Genetic diversity in Spanish donkey breeds using microsatellite DNA markers[J].Genetics Selection Evolution,2001,33（4）:1-10.

[34]Arangurenméndez J,Gómez M,Jordana J.Hierarchical analysis of genetic structure in Spanish donkey breeds using microsatellite markers[J].Heredity,2002,89（3）:207-211.

[35]Matassino D,Cecchi F,Ciani F,et al.Genetic Diversity and Variability in two Italian Autochthonous Donkey Genetic Types Assessed by Microsatellite Markers[J].Italian Journal of Animal Science,2014,13（1）:53-60.

[36]朱文進,張美俊,葛慕湘,等.中國8個地方驢種遺傳多樣性和系統(tǒng)發(fā)生關(guān)系的微衛(wèi)星分析[J].中國農(nóng)業(yè)科學,2006,39（2）:398-405.

[37]阿吉,王金富,路立里,等.新疆驢地方類群微衛(wèi)星標記與體尺性狀的相關(guān)分析[J].石河子大學學報（自科版）,2007,25（4）:449-52.

[38]劉艷艷,譚晴晴,范陽陽,等.烏頭驢和三粉驢微衛(wèi)星遺傳多態(tài)性分析[J].中國畜牧雜志,2017,53（1）:44-50.

[39]Di R,Liu Q Y,Xie F,et al.Evaluation of genetic diversity and population structure of five Chinese indigenous donkey breeds using microsatellite markers[J].Czech Journal of Animal Science,2017,62（5）:219-225.

[40]Liu D H,Han H Y,Zhang X,et al.The genetic diversity analysis in the donkey myostatin gene[J].Journal of Integrative Agriculture,2017,16（3）:656-663.

[41]Chen J,Song Z,Rong M,et al.The association analysis between Cytb,polymorphism and growth traits in three Chinese donkey breeds[J].Livestock Science,2009，126（1-3）:306-309.

[42]朱文進,陳娟,關(guān)學敏,等.驢生長激素基因的克隆與分析[J].中國農(nóng)業(yè)科學,2011,44（3）:579-586.

[43]王顏顏,托乎提·阿及德,肖海霞,等.新疆良種驢DGAT2基因第3內(nèi)含子PCR-SSCP多態(tài)性與體尺性狀的相關(guān)性分析[J].石河子大學學報（自科版）,2011,29（1）:40-44.

[44]王艷萍,高帥,劉迥,等.年齡、部位對新疆驢皮中Ⅰ型膠原蛋白表達影響白研究[J].中國畜牧雜志,2017（12）:41-44.

[45]Cosenza G,Pauciullo A,Annunziata A L,et al.Identification and characterization of the donkey CSN1S2 I and II cDNAs[J].Italian Journal of Animal Science,2010,9（2）:206-211.

[46]Viitala S,Szyda J,Blott S,et al.The role of the bovine growth hormone receptor and prolactin receptor genes in milk,fat and protein production in Finnish Ayrshire dairy cattle[J].Genetics,2006,173（4）:2151-2164.

[47]畢蘭舒,肖海霞,臧長江,等.疆岳驢催乳素受體基因多態(tài)性及其與泌乳性狀間的相關(guān)性分析[J].中國畜牧獸醫(yī),2017,44（1）:180-185.

[48]Godovaczimmermann J,Conti A,Sheil M.Covalent structure of the minor monomeric β-lactoglobulin II component from donkey milk[J].Biol Chem Hoppe Seyler,1990,371（2）:871-879.

[49]Herrouin M,Mollé D,Fauquant J,et al.New Genetic Variants Identified in Donkey's Milk Whey Proteins[J].Journal of Protein Chemistry,2000,19（2）:105.

[50]Cunsolo V,Costa A,Saletti R,et al.Detection and sequence determination of a new variant beta-lactoglobulin II from donkey[J].Rapid Communications in Mass Spectrometry,2007,21（8）:1438–1446.

[51]Chianese L,Simone C D,Ferranti P,et al.Occurrence of qualitative and quantitative polymorphism at donkey beta-Lactoglobulin II locus[J].Food Research International,2013,54（1）:1273-1279.

[52]Criscione A,Cunsolo V,Tumino S,et al.Polymorphism at donkey β-lactoglobulin II locus:identification and characterization of a new genetic variant with a very low expression[J].Amino Acids,2018，1-12.

[53]Haase B,Rieder S,Leeb T.Two variants in the KIT gene as candidate causative mutations for a dominant white and a white spotting phenotype in the donkey[J].Animal Genetics,2015,46（3）:321.

[54]Utzeri V J,Bertolini F,Ribani A,et al.The albinism of the feral Asinara white donkeys （Equus asinus）is determined by a missense mutation in a highly conserved position of the tyrosinase （TYR）gene deduced protein[J].Animal Genetics,2016,47（1）:120-124.

[55]Abitbol M,Legrand R,Tiret L.A missense mutation in the agouti signaling protein gene （ASIP）is associated with the no light points coat phenotype in donkeys[J].Genetics Selection Evolution,2015,47（1）:28.

[56]何良軍,肖國亮,庫爾班·阿木提,等.驢同期發(fā)情初探[J].中國草食動物,2008,28（1）:34-34.

[57]何良軍,張建俊,王時偉.PMSG對驢卵泡發(fā)育的影響[J].中國草食動物,2009,29（6）:66.

[58]張瑞濤.德州驢冷凍精液及其主要影響因素研究[D].濟南：山東師范大學,2014.

[59]王慧偉,陳靜波,霍飛,等.驢細管精液冷凍效果研究[J].草食家畜,2008，（3）:18-20.

[60]Rota A,Panzani D,Sabatini C,et al.Donkey jack （Equus asinus）semen cryopreservation:studies of seminal parameters,post breeding inflammatory response,and fertility in donkey jennies[J].Theriogenology,2012,78（8）:1846-1854.

[61]De Oliveira J V,Freire Oliveira P V,Cm M E O,et al.Strategies to improve the fertility of fresh and frozen donkey semen[J].Theriogenology,2016,85（7）:1267-1273.

[62]張瑞濤,張燕,高運東,等.驢精液采集與冷凍相關(guān)技術(shù)研究[J].安徽農(nóng)學通報,2013,19（3）:121-124.

[63]Brinsko S P,Rigby S L,Lindsey A C,et al.Pregnancy rates in mares following hysteroscopic or transrectally-guided insemination with low sperm numbers at the utero-tubal papilla[J].Theriogenology,2003,59（3-4）:1001.

[64]Lindsey A C,Varner D D,Jr S G,et al.Hysteroscopic or rectally guided,deep-uterine insemination of mares with spermatozoa stored 18 h at either 5 degrees C or 15 degrees C prior to flow-cytometric sorting[J].Animal Reproduction Science,2005,85（1-2）:125.

[65]Hayden S S,Blanchard T L,Brinsko S P,et al.Pregnancy rates in mares inseminated with 0.5 or 1 million sperm using hysteroscopic or transrectally guided deep-horn insemination techniques[J].Theriogenology,2012,78（4）:914-920.

[66]何良軍,張建俊,肖國亮.驢發(fā)情期卵泡變化初探[J].經(jīng)濟動物學報,2008,12（3）:153-155.

[67]Goudet G,Douet C,Kaabouba-Escurier A,et al.Establishment of conditions for ovum pick up and IVM of jennies oocytes toward the setting up of efficient IVF and in vitro embryos culture procedures in donkey （Equus asinus）[J].Theriogenology,2016,86（2）:528-535.

[68]杜興雨,王彥妮,張焱如,等.不同濃度玻璃化冷凍保護液及冷凍載體對驢卵母細胞成熟率的影響[J].畜牧與飼料科學,2016,37（3）:4-6.