青海駱駝線粒體DNA D-loop區遺傳多樣性研究

高雪 晁生玉 王啟菊 孟克達拉 烏蘭巴特爾 賈功雪

摘 要 旨在研究青海駱駝遺傳多樣性，揭示青海駱駝的母系起源進化。采集柴達木地區3個青海駱駝類群耳組織樣品88份并提取基因組DNA，利用PCR擴增和基因測序分析線粒體DNA D-loop序列，并與蒙古、內蒙古雙峰駝進行比較。結果表明：D-loop序列中共檢測到96個多態位點，定義了27種單倍型。3個青海駱駝類群占有16個單倍型，而蒙古雙峰駝獨有7個單倍型，內蒙古雙峰駝獨有4個單倍型。系統發育分析顯示出兩個獨立的分支：第一支為青海駱駝3個類群與蒙古雙峰駝，且德令哈雙峰駝與其他兩個類群間存在明顯界限；第二支為內蒙古雙峰駝。青海駱駝與內蒙古雙峰駝的遺傳距離均較遠，但與蒙古雙峰駝的遺傳距離較近。因此，青海駱駝母系起源更傾向于蒙古雙峰駝而非內蒙古雙峰駝。

關鍵詞 柴達木雙峰駝；線粒體DNA；系統發育；遺傳多樣性；單倍型

青海駱駝是中國優良的地方雙峰駝（Bactrian camel，Camelus bactrianus）品種之一，主要分布于青海省柴達木盆地境內的荒漠、半荒漠或半荒漠灌叢草場，是世界分布海拔最高的駱駝品種。作為當地重要的生產生活資料，青海駱駝的產肉、產絨、產乳、役用性能對于社會經濟發展做出了重要貢獻[1]。然而，隨著飼養模式的落后、養殖收益的下降與經濟轉化效率的不足，青海駱駝數量迅速由1983年的28 100峰降至2015年的6 668峰[2]。此后隨著品種保護工作的重視，青海駱駝數量得到了逐步的恢復，2018年已恢復至10 304峰。

關于青海駱駝的起源一直存有爭議，目前較為明確的是青海土著駱駝自13世紀以來長期經歷蒙古駱駝的血液滲透，形成了以蒙古駱駝為主要成分的混血種[3]。線粒體DNA（Mitochondrion DNA，mtDNA）具備嚴格的母系遺傳特征并且在世代傳遞中幾乎不發生重組，因此常被用于研究哺乳動物品種起源與系統發育。mtDNA D-loop區由于進化速度快、變異水平高，是最常用的種群遺傳多樣性研究標記之一[4-5]。多項研究利用mtDNA探討了國內外駱駝品種的遺傳多樣性與種群進化關系[6-7]，然而青海駱駝的品種起源與遺傳特征仍未得到充分闡明。因此，本研究借助mtDNA D-loop序列對青海駱駝的3個典型類群（都蘭、德令哈、莫河）進行分析，并與蒙古、內蒙古雙峰駝D-loop序列進行比較，以闡明青海駱駝不同類群及其與蒙古、內蒙古雙峰駝間的進化關系，為青海駱駝資源的保護與利用提供理論依據。

1 材料與方法

1.1 試驗動物及樣品

本試驗共采集88份青海駱駝耳組織樣品，覆蓋范圍包括都蘭、德令哈和莫河3個典型類群，分別采自青海省海西蒙古族藏族自治州德令哈市蓄集鄉、都蘭縣香日德鎮與烏蘭縣茶卡鎮。樣品保存于裝有75%酒精的離心管中，運回實驗室后于-80 ℃條件下保存。采樣點及樣本信息見表1。

1.2 儀器及試劑

StarSpin柱式動物DNA提取試劑盒（D111）、StarPrep快速膠回收試劑盒（D205）與 ?2×Taq PCR StarMix試劑盒（A012）購自北京康潤誠業生物科技有限公司；Nanodrop超微量分光光度計（ND1000）購自賽默飛世爾科技（中國）有限公司；PCR儀（A200）購自杭州朗基科學儀器有限公司。

1.3 總DNA的提取

取青海駱駝耳組織樣品0.1 g，使用DNA提取試劑盒提取總DNA，并使用超微量分光光度計測定DNA濃度與質量，隨后將樣品濃度稀釋至100 ng/μL。

1.4 D-loop序列的擴增與測序

使用引物D-loop-F（5′-ATACTGATTACG- CTGGTCTTG-3′）和D-loop-R（5′-ACTAATAGAAAGGCTGGGATC-3′）擴增mtDNA D-loop區序列片段。PCR擴增體系如下：2×Taq PCR StarMix 25.0 μL，上、下游引物各1.0 μL，DNA樣品2.0 μL，最后用ddH2O補齊至50.0 μL。PCR擴增程序如下：94 ℃預變性2 min； ?94 ℃變性30 s，59 ℃退火30 s，72 ℃延伸1 min，共30個循環；72 ℃終延伸10 min。PCR擴增產物使用20 g/L瓊脂糖凝膠電泳，觀察確認目的片段后切下膠塊，使用膠回收試劑盒進行回收純化，純化樣品送至生工生物工程（上海）股份有限公司進行Sanger測序。

1.5 統計分析

測序結果使用Chromas 2.6軟件進行原始序列峰圖檢測?；谇嗪ｑ橊勯L期存在蒙古駱駝血液滲透的觀點，自NCBI下載了蒙古與內蒙古雙峰駝D-loop序列（表2），與青海駱駝測序結果匯總分析。在MEGA 7.0軟件中使用MUSCLE比對序列同源性并校正，計算序列的核苷酸組成與遺傳距離；使用DnaSP 5.10軟件對多態位點數、單倍型數、核苷酸多樣性、單倍型多樣性、單倍型數、平均核苷酸差異、中性檢驗及遺傳分化指數等進行統計。以野生雙峰駝（wild Bactrian camel，Camelus ferus）D-loop序列NC_009628.2為外群，構建鄰接法（Neighbor-Joining，NJ）系統發育樹，并用iTOL網站（https：//itol.embl.de/）美化；使用Network 4.2軟件繪制單倍型網絡圖。

2 結果與分析

2.1 青海駱駝D-loop序列DNA多態性分析

經PCR擴增，3個青海駱駝類群獲得的D-loop序列片段與預期目的片段大?。? 300～ ?1 500 bp）一致，且條帶單一（圖1），保證了后續分析的可靠性。以野生雙峰駝D-loop序列為參考，分別比對3個青海駱駝類群測序結果與蒙古、內蒙古雙峰駝序列。經校正與修剪后，共獲得118條長度為759 bp的D-loop序列（圖2）。

序列中4種堿基（T、C、G和A）的平均含量分別為24.7%、27.9%、28.9%和18.5%，A+T含量（53.6%）高于C+G含量（46.4%），符合mtDNA D-loop區富含A/T堿基的序列特征。5個類群共檢測到96個多態位點，整個群體的核苷酸多樣性、核苷酸平均差異數和單倍型多樣性分別為0.014 93、10.851和0.845。在3個青海駱駝類群中，都蘭雙峰駝具有最高的核苷酸多樣性，莫河雙峰駝次之，而德令哈雙峰駝最低。與之對應，中性檢驗結果顯示，都蘭和莫河雙峰駝TajimaD結果均為負值且顯著偏離中性（P < ?0.05），表明二者在形成進程中均發生過明顯的種群擴張，而德令哈雙峰駝群體則較為保守（表3）。

2.2 雙峰駝群體歷史動態分析

遺傳距離與遺傳分化指數分析顯示，3個青海駱駝類群間的遺傳距離較為接近，而都蘭雙峰駝與莫河雙峰駝間的遺傳分化指數最小，表明二者之間可能發生一定程度的基因交流（表4）。內蒙古雙峰駝與其他4個類群的遺傳差異均遠大于蒙古雙峰駝，進一步證實了青海駱駝可能與蒙古雙峰駝而非內蒙古雙峰駝存在更為密切的親緣關系。

2.3 系統發育樹構建與單倍型網絡關系分析

為了進一步分析青海駱駝3個類群的系統發育關系，以野生雙峰駝mtDNA D-loop序列為外群構建了NJ系統發育樹。結果顯示，所有樣本明顯分為兩個獨立的分支：第一支分支是青海駱駝3個類群與蒙古雙峰駝，另一支則是內蒙古雙峰駝，兩支在進化上存在明顯界限（圖3）。

由表5可知，單倍型分析共界定了27種單倍型Hap1～Hap27。其中：都蘭雙峰駝有7種，德令哈雙峰駝有6種，莫河雙峰駝有9種，蒙古雙峰駝有7種，內蒙古雙峰駝有4種。27個單倍型中，Hap5和Hap8在不同雙峰駝類群中出現的頻率最高，在3個青海駱駝類群中均有分布。

單倍型網絡分析呈現的不同類群親緣關系與系統發育樹結果基本一致。內蒙古雙峰駝與其他類群之間明顯存在差異，相比之下蒙古雙峰駝與青海駱駝則更為接近。值得關注的是，青海駱駝3個類群間存在單倍型共享：Hap6為都蘭和莫河雙峰駝共享；Hap9為都蘭和德令哈雙峰駝共享；Hap5與Hap8則在3個類群中均有發現（圖4）。

3 討? 論

厘清品種的系統進化關系與遺傳多樣性是家畜遺傳資源保護的必要前提。圍繞雙峰駝的研究早期多采用細胞染色體核型分析與血液蛋白多態性分析的手段，獲得的信息量有限[8-9]。近年來，借助基因組DNA測序手段與多種分子標記的輔助分析，更多的遺傳信息被揭示[10]。程佳等[11]使用線粒體細胞色素b部分序列檢測了5個雙峰駝類群，分析表明中國家養雙峰駝起源于同一母系，但與現存的野生雙峰駝祖先并不相同。這一結論隨后也被野生和家養雙峰駝的全基因組測序結果所證實[12]。何曉紅等[13]使用微衛星引物分析了10個雙峰駝類群的遺傳關系，證實新疆的雙峰駝類群聚為一類，而青海、甘肅、內蒙古與蒙古雙峰駝具有較近的親緣關系。而另一項針對中國、蒙古與俄羅斯等多國的11個雙峰駝類群進行比較的結果則證實，不同地域間的家養雙峰駝類群分化特征并不明顯，存在一定的基因交流[14]。

本研究基于遺傳距離與遺傳分化指數分析進一步表明，在家養雙峰駝內部依然存在一定程度的遺傳分化，青海駱駝、蒙古雙峰駝與內蒙古雙峰駝處于不同的進化分支。據記載，青海駱駝的起源至少有3個：早期駱駝可能隨吐谷渾人、蒙古人征戰由內蒙古帶入青海；解放前部分哈薩克人由新疆帶駝進入柴達木盆地；而解放初期青海由甘肅購進數千峰駱駝屯牧繁育于盆地[15]。本研究進一步比較青海駱駝的3個類群發現，都蘭雙峰駝與莫河雙峰駝存在相對較多的基因交流，而德令哈雙峰駝顯示出與其他兩個類群明顯不同的分化方向。研究結果與先前研究存在一定的差異[13，16]，但與白俊艷等[17]基于表型測定結果的聚類分析相符，這可能與之前試驗采用的樣品數量、地理區域與分子標記選擇相關。因此，還需要通過進一步研究分析青海駱駝不同起源的可能性。

4 結? 論

從母源進化上看，青海駱駝與蒙古雙峰駝間遺傳分化程度較小，但與內蒙古雙峰駝間存在較大的遺傳差異，其起源可能更傾向于蒙古雙峰駝而非內蒙古雙峰駝。青海駱駝3個類群間存在一定的基因交流，但是德令哈雙峰駝表現出與其他兩個類群間明顯的遺傳分化。研究進一步闡明了青海駱駝的遺傳多樣性與系統發育關系，對于青海駱駝的品種保護與資源開發具有指導意義。

參考文獻 Reference：

［1］ 張玲勤.青海駱駝亟待保護[J].中國畜牧雜志，2007， 43（5）：63-64.

ZHANG L Q.Qinghai camels need urgent protection[J].Chinese Journal of Animal Science，2007，43（5）：63-64.

[2] 艾德強，閻 萍.青海駱駝種質資源保護現狀與開發利用對策[J].中國草食動物科學，2018，38（2）：68-70.

AI D Q，YAN P.Protection situation and utilization countermeasures of Qinghai camel germplasm resources[J].Chinese Herbivore Science，2018，38（2）：68-70.

[3] 賀新民.青海駱駝來源和發展的探討——兼論羌族來源及其畜牧簡史[J].農業考古，1995（1）：292-298.

HE X M.Discussion on the origin and development of Qinghai camel：also about the origin of Qiang people and their animal husbandry[J].Agricultural Archaeology，1995 （1）：292-298.

[4] 李新建，薛亞輝，王明宇，等.河南地方豬mtDNA D-loop序列遺傳多樣性及系統進化分析[J].河南農業大學學報，2020，54（3）：452-461.

LI X J，XUE Y H，WANG M Y，et al. Genetic diversity and phylogenetic relationships of Henan native pig breeds based on the sequence of mtDNA D-loop region[J].Journal of Henan Agricultural University，2020，54（3）：452-461.

[5] 董澤生，高 雪，吳燕華，等.藏系綿羊mtDNA和Y染色體遺傳多樣性的檢測與系統進化關系的分析[J].黑龍江畜牧獸醫，2019，568（4）：45-47，51，179.

DONG Z? SH，GAO X，WU Y H，et al.Genetic diversity and phylogenetic evolution of mitochondrial DNA and? ?Y-chromosome in Tibetan sheep[J].Heilongjiang Animal Science and Veterinary Medicine，2019，568（4）：45-47，51，179.

[6] MING L，YI L，SA R，et al.Genetic diversity and phylogeographic structure of Bactrian camels shown by mitochondrial sequence variations[J].Animal Genetics，2017，48（2）：217-220.

[7] 張成東，任戰軍，楊雪嬌.蘇尼特家養雙峰駝mtDNA Cyt b基因和D-loop序列的遺傳多樣性[J].西北農林科技大學學報（自然科學版），2015，43（3）：32-42.

ZHANG CH D，REN ZH J，YANG X J.Genetic diversity of mtDNA Cyt b gene and D-loop sequence of Sunite domestic bactrian camel[J].Journal of Northwest A&F University （Natural Science Edition），2015，43（3）：32-42.

[8] 張才駿.青海雙峰駝血液蛋白質多態性的研究[J].中國畜牧雜志，1996，32（5）：33-34.

ZHANG C J.Study on blood protein polymorphism of Qinghai Bactrian camels[J].Chinese Journal of Animal Science，1996，32（5）：33-34.

[9] BALMUS G，TRIFONOV V A，BILTUEVA L S，et al.Cross-species chromosome painting among camel，cattle，pig and human：further insights into the putative Cetartiodactyla ancestral karyotype[J].Chromosome Research，2007，15（4）：499-515.

[10] [ZK（#]明 亮，伊 麗，何 靜，等.雙峰駝起源與進化的分子遺傳學研究進展[J].家畜生態學報，2017，38（3）：5-9.

MING L，YI L，HE J，et al.Research progress on molecular genetics of origin and domestication of Bactrian camel[J].Acta Ecologiae Animals Domastici，2017，38（3）： ?5-9.

[11] 程 佳，任戰軍，王 樂，等.基于Cytb基因的家養雙峰駝分子系統發育研究[J].西北農林科技大學學報（自然科學版），2009，37（12）：17-21.

CHENG J，REN ZH J，WANG L，et al.Molecular phylogeny of domestic bactrian camel revealed by partial cytochrome b gene[J].Journal of Northwest A&F University （Natural Science Edition），2009，37（12）：17-21.

[12] JIRIMUTU，WANG ZH，DING G H，et al.Genome sequences of wild and domestic bactrian camels[J].Nature Communication，2012，3：1202.

[13] 何曉紅，韓秀麗，關偉軍，等.采用微衛星標記分析10個雙峰駝群體的遺傳變異和群體間的遺傳關系[J].生物多樣性，2012，20（2）：199-206.

HE X H，HAN X L，GUAN W J，et al.Genetic variability and relationship of 10 Bactrian camel populations revealed by microsatellite markers[J].Biodiversity Science，2012，20（2）：199-206.

[14] MING L，YI L，GUO F C，et al.Molecular phylogeny of the Bactrian camel based on mitochondrial Cytochrome b gene sequences[J].Genetics and Molecular Research，2016，15（3）：15038983.

[15] 焦小鹿，寧金友，馬清德，等.青海省畜禽遺傳資源志[M].西寧：青海人民出版社，2013：104-108.

JIAO X L，NING J Y，MA Q D，et al.Genetic Resources of Livestock and Poultry in Qinghai? Province[M].Xining：Qinghai Peoples? Publishing House，2013：104-108.

[16] 柏 麗，周俊文，李新海，等.中國6個雙峰駝群體微衛星遺傳多樣性分析[J].中國畜牧獸醫，2017，44（6）：1734-1745.

BAI L，ZHOU J W，LI X H，et al.Analysis on the genetic diversity of six Bactrian camel populations in China? ?using microsatellite marker[J].China Animal Husbandry and Veterinary Medicine，2017，44（6）：1734-1745.

[17] 白俊艷，烏仁套迪，巴 拉，等.中國4個雙峰駝品種的聚類分析[J].黑龍江畜牧獸醫，2015，485（17）：236-238.

BAI J Y，WUREN T，BA L，et al.Clustering analysis of four Bactrian camel populations in China[J].Heilongjiang Animal Science and Veterinary Medicine，2015，485（17）：236-238.

Genetic Diversity of Mitochondrial DNA D-loop Region in Qinghai Camel

Abstract To evaluate the genetic diversity and maternal origin of Qinghai camel，a total of 88 ear tissue samples were collected from three distributed ecotypes of Qinghai camel.Genomic DNA extraction was performed to obtain the samples，after which the D-loop region on mitochondrial DNA was amplified and sequenced，and sequences were compared with those of Mongolian and Inner Mongolian Bactrian camels.The results showed that a total of 96 polymorphic loci were detected in D-loop sequences and 27 haplotypes were identified.There were 16 haplotypes shared by Qinghai camel.The Mongolian Bactrian camel had seven haplotypes，while the Inner Mongolian Bactrian camel had four haplotypes.Phylogenetic analysis revealed two distinct branches.The first branch consisted of three ecotypes of Qinghai camel and Mongolian Bactrian camel，with the Delingha Bactrian camel exhibiting a clear distinction from the other two ecotypes.The second branch was predominantly represented by Inner Mongolian Bactrian camel.The genetic distance between Qinghai camel and Inner Mongolia Bactrian camel was significantly greater than that between Qinghai camel and Mongolia Bactrian camel.In conclusion，the maternal origin of Qinghai camel is more closely related to Mongolian Bactrian camel than to the Inner Mongolian Bactrian camel.

Key words Qaidam Bactrian camel; Mitochondrial DNA; Phylogeny; Genetic diversity; Haplotype