利用線粒體DNA控制區鑒別狼與狗可行性分析

江泓慧 李 波，2* 馬 躍，2 徐艷春，2*

(1.東北林業大學野生動物資源學院，哈爾濱，150040;2.國家林業局野生動植物檢測中心，哈爾濱，150040)

灰狼(Canis lupus)發源于距今300多萬年前的更新世，由湯氏熊分化為現代犬屬動物［1］。狗(Canis familiaris)是狼的近親，通過對狼和狗表現型和基因型的研究發現狼是狗的野生祖先［2－3］，即狗由早期人類從灰狼馴化而來［4］，馴養始于亞洲，距今不超過16300年［5］。狼和狗具有相同的染色體組型，可以雜交并產生后代，有學者將二者描述為同一物種［6］。

狼作為數量有限的野生動物受到眾多國家法律保護，已被列入《世界自然保護聯盟瀕危物種紅色名錄》，在世界范圍內受到嚴格的保護和管理，禁止非法捕殺及販賣。狗則被視為人類的伴侶動物，不被野生動物相關法律保護。狼皮、牙、胃等在非法市場具有較大的需求，非法偷獵，走私狼產品在中國境內和周邊地區時有發生。當執法部門查獲相關的狼產品時，嫌疑人為逃避法律懲罰經常謊稱其產品來源于狗。為打擊非法偷獵、走私活動，保護野生狼的種群，急需建立一套有效的狼與狗鑒別方法。

傳統的物種鑒別方法可以分為宏觀和微觀兩種。宏觀鑒別中的形態學［7］及行為學［8－9］方法是快速且準確的。但通常需要涉案檢材有完整的特征。當所繳獲物品為動物的內臟或局部組織時，則需要采取微觀的分子生物學技術進行鑒別。其中，最常見的是利用線粒體DNA(mtDNA)作為分子標記的方法。多拷貝的mtDNA沒有基因重組現象，嚴格遵循母系遺傳，其進化速率是單拷貝核基因的5～10倍［10］，它是動物系統發育分析和物種識別常用的工具。通過mtDNA，包括線粒體細胞色素b基因(Cyt b)及控制區［11－13］等對狼和狗的研究始于20世紀。最近的犬科動物分子系統學研究發現犬屬的狼和狗等為系統樹中獨立的分支，其分支時間較赤狐(Vulpes vulpes)、藍狐(Alopex lagopus)和貉(Nyctereutes procyonoides)早［14］，驗證了狼和狗具有較其他犬科動物更為親近的關系。也有研究比較了狼和狗線粒體全基因組的結構特點，并進一步探索其物種起源與進化［15］。

在物種識別方面，有人利用mtDNA控制區對克羅地亞境內狼和狗進行鑒別，結果證明該段序列可以成功將狼和狗鑒別開來［16］。由此，本研究比對中國境內及周邊國家狼和狗mtDNA控制區序列，通過序列差異比較和系統發育分析探索利用控制區對狼和狗進行物種鑒別的可行性。該分子方法的建立可以豐富狼和狗鑒別手段，為動物法醫學相關鑒別提供助力［17］。

1 材料與方法

1.1 試驗材料

試驗材料來源于國家林業局野生動植物檢測中心保存的狼和狗肌肉組織樣本。全部樣本來源于中國境內。其中56份樣本為狗，18份樣本為狼。另有來源于日本、蒙古、阿富汗等國的24條狗mtDNA控制區序列，8條狼 mtDNA控制區序列及2條外群灰狐(Urocyon cinereoargenteus)mtDNA控制區序列下載自GenBank。經統計共108份試驗材料，80份來自于狗，26份來自于狼，2份來自于灰狐。本實驗所用全部試驗材料信息詳見表1。

1.2 DNA提取

肌肉組織采用AxyPrep基因組DNA提取試劑盒(康寧生命科學有限公司，吳江)提取，于－20℃條件下保存。

1.3 PCR擴增

使用引物 DL 14和 DH 7［18］擴增目標 mtDNA片段。PCR擴增體系為50 μl，包含10 ng模版DNA，1×EasyTaq Mix(北京全式金生物技術有限公司，北京)及上下游引物各200 nM。使用Eppendorf Mastercycler Pro PCR儀進行擴增，擴增條件為:94℃預變性5 min，隨后35個循環(94℃/30 s，50℃/30 s，72℃/60 s)，最終72℃保溫5 min。采用AxyPrep DNA凝膠回收試劑盒(康寧生命科學有限公司，吳江)，按照操作手冊進行純化回收。

1.4 數據分析

純化產物由哈爾濱新海基因生物有限公司利用3730測序儀采用一代測序技術測定，以引物DL 14單向測序進行分析測定。用MEGA 5軟件對所測序列進行比對，截取并保留序列約676 bp的片段，分析單倍型成對序列遺傳距離;用Arlequin 3.5.2.2［19］分析單倍型頻率;用dnasp 5分析多態位點，單倍型多樣性和核酸多樣性;用IQ－TREE和MrBayes3.2.6軟件分別構建ML系統發育樹和貝葉斯樹，選擇兩個灰狐單倍型為外群。

2 結果

2.1 序列分析

分析106條狼和狗mtDNA控制區約676 bp部分序列，共得到47個單倍型結果，其中12個來源于狼，35個來源于狗，且狼和狗之間沒有共享單倍型。這些單倍型包含了62個(9.17%)多態位點(表2)。狼有51個多態位點(7.54%)，狗有33個多態位點(4.88%)。狼有29個特異多態位點，狗有11個特異多態位點，二者占全部多態位點的64.52%。分析序列核苷酸多樣性(Pi)和單倍型多樣性(Hd)，結果表明狼 Pi值為 0.015，Hd值為 0.914，狗 Pi值為0.012，Hd值為0.957。狼的單倍型頻率介于3.85%～15.40%之間，狗的單倍型頻率介于1.25% ～10.00%之間(表3)，狼和狗種內單倍型頻率主要集中于低頻區域。表明不同的mtDNA控制區單倍型在狼和狗各自世系間分布較為均勻，且存在優勢單倍型。

2.2 系統發育分析

ML和BI方法構建的單倍型系統發育樹拓撲結構一致，各分支的支持度存在差異(圖1)。系統發育分析表明，全部狼和狗的單倍型獨立于外群形成2大分支或世系(支持度為100%)。第一分支中包含3個狼單倍型、8個狗單倍型，其中3個狼單倍型(WN03、AB007373、AB007374，25%)與5個狗單倍型共享同一世系。第二分支中包含36個單倍型(9個狼單倍型、27個狗單倍型)，其中狼單倍型(75%)與狗單倍型分屬不同分支。且狼單倍型 AB007375，WE01和WB01，WB02和 WC05總是獨立成支(支持度 ＞80%)。

3 討論

mtDNA被廣泛用于動物遺傳進化，系統發育，群體遺傳結構等方面的研究［13，20－23］。分別有研究對青海湖附近18頭狼［13］及24個不同品種共73只狗［24］的mtDNA控制區序列進行研究。結果表明同一地區狼之間及不同品種狗之間mtDNA控制區均具有廣泛的多態性，證明二者遺傳多樣性較高。為比較狼和狗mtDNA控制區單倍型并分析得到更加廣泛的多態性提供了可能性依據。本研究分析狼和狗單倍型序列多態位點的差異，結果狼有29個特異性位點，狗有11個特異性位點，證明狼和狗之間具有更高的遺傳多樣性。二者的顯著差異是本方法鑒別狼和狗的基礎，說明在mtDNA控制區676 bp段序列范圍內，狼和狗具有各自獨立的進化，為本研究鑒別方法提供了可行性。

利用mtDNA進行動物系統發育分析，有助于物種的起源［25］及分類地位［26］等方面的研究，對于不同種動物則可進一步對進化多樣性及物種鑒別進行研究。已有研究利用mtDNA控制區對羊、滇金絲猴(Rhinopithecus bieti)、沙氏下鱵魚(Hyporhamphus sajori)等動物的進化多樣性及物種鑒別［27－31］展開研究。在犬科動物方向，有研究利用mtDNA控制區對濟州犬與其他犬種的種內系統發育關系進行了研究［32］;在利用mtDNA控制區對克羅地亞地區的狼和狗進行種間區分的研究中，系統發育分析結果表示狼和狗具有各自獨立的世系，可以利用此方法鑒別狼和狗［19］。本研究利用全部47個狼和狗單倍型構建系統發育樹，在來自中國境內及周邊國家的12個狼單倍型中有9個(75%)與狗單倍型均分屬不同分支。這說明系統發育樹中單倍型的世系分布也具有一定的鑒別力。而另一方面，長期進化發展及基因交流等原因為狼和狗的基因分化帶來一定的干擾因素，導致少部分狼(25%)和狗具有相似的mtDNA單倍型，使得本方法的鑒別準確度降低。因而根據狼和狗mtDNA控制區單倍型在系統發育樹中的分布結果鑒別狼和狗，準確率較高，但有無法鑒別的情況存在。

綜合特異性位點和單倍型世系分布分析結果，推測利用mtDNA控制區鑒別狼和狗是可行的。但mtDNA母系遺傳的特點使得本方法在雜合子鑒別方面存在障礙，如本研究全部(75%部分和25%部分)狼單倍型個體中是否存在雌狼和雄狗雜交的后代，或者雌狼和雜種后代產生的后代，需要進一步綜合常染色體或Y染色體等遺傳標記［33］進行父系鑒別。綜上，本研究方法對狼和狗的鑒別具有一定的可行性及參考意義，但仍需要更廣泛地選取樣本進行深入研究。

［1］ 高中信，袁力，付松華．狗與狼 ［J］．野生動物，1997，18(6):29－30．

［2］ Vilà C，Savolainen P，Maldonado J E，et al．Multiple and ancient origins of the domestic dog［J］．Science，1997，276(5319):1687－1689．

［3］ Wayne R K，Ostrander E A．Origin，genetic diversity，and genome structure of the domestic dog ［J］．BioEssays，1999，21(3):247－257．

［4］ Leonard J A，Wayne R K，Wheeler J，et al．Ancient DNA evidence for Old World origin of New World dogs［J］．Science，2002，298(5598):1613－1616．

［5］ Parker H G，Kim L V，Sutter N B，et al．Genetic structure of the purebred domestic dog［J］．Science，2004，304(5674):1160 －1164．

［6］ Vilà C，Wayne R K．Hybridization between wolves and dogs［J］．Conservation Biology，1999，13(1):195 －198．

［7］ Janssens L，Spanoghe I，Miller R，et al．Can orbital angle morphology distinguish dogs from wolves［J］．Zoomorphology，2016，135(1):149－158．

［8］ Range F，Virányi Z．Social learning from humans or conspecifics:differences and similarities between wolves and dogs［J］．Frontiers Psychology，2013，4:868．

［9］ Udell M A R．When dogs look back:inhibition of independent problem-solving behaviour in domestic dogs(Canis lupus familiaris)compared with wolves(Canis lupus) ［J］．Biology Letters，2015，11(9):20150489．

［10］ Brown W M，George M，Wilson A C．Rapid evolution of animal mitochondrial DNA ［J］．Proceedings of the National Academy of Sciences of the United States of America，1979，76(4):1967 －1971．

［11］ Pereira L，Van Asch B，Amorim A．Standardisation of nomenclature for dog mtDNA D-loop:a prerequisite for launching a Canis familiaris database ［J］．Forensic Science International，2004，141(2/3):99－108．

［12］ 竇海龍，沙未來，陳磊，等．基于線粒體D-Loop區的青海狼種群遺傳學研究［C］∥中國動物學會，中國生態學學會，中國野生動物保護協會．第八屆全國野生動物生態與資源保護學術研討會論文集，2012:152－153．

［13］ Eichmann C，Parson W．Molecular characterization of the canine mitochondrial DNA control region for forensic applications［J］．International Journal of Legal Medicine，2007，121(5):411－416．

［14］ 陳永久，張亞平，鄒希明．犬科的線粒體細胞色素b DNA序列及其分子系統學研究［J］．遺傳學報，2000，27(1):7－11．

［15］ 孟超．中國狼、犬科及部分食肉目動物線粒體全基因組研究 ［D］．曲阜:曲阜師范大學，2009．

［16］Sindiˇcic'M，Gomerˇcic'T，Galov A，et al．Mitochondrial DNA control region as a tool for species identification and distinction between wolves and dogs from Croatia ［J］．Veterinarski Arhiv，2011，81(2):249－258．

［17］ Budowle B，Garofano P，Hellman A，et al．Recommendations for animal DNA forensic and identity testing［J］．International Journal of Legal Medicine，2005，119(5):295－302．

［18］ Tsuda K，Kikkawa Y，Yonekawa H，et al．Extensive interbreeding occurred among multiple matriarchal ancestors during the domestication of dogs:evidence from inter-and intraspecies polymorphisms in theD-loop region ofmitochondrialDNA between dogsand wolves［J］．Genes ＆ Genetic Systems，1997，72(4):229 －238．

［19］ Excoffier L，Laval G，Schneider S．Arlequin(version 3.0):anintegrated software package for population genetics data analysis［J］．Evolutionary Bioinformatics Online，2005，1(4A):47－50．

［20］ Kim K S，Lee S E，Jeong H W，et al．The complete nucleotide sequence ofthe domestic dog(Canisfamiliaris)mitochondrial genome［J］．Molecular Phylogenetics and Evolution，1998，10(2):210－220．

［21］ Bjornerfeldt S，Webster M T，Vilà C．Relaxation of selective constraint on dog mitochondrial DNA following domestication［J］．Genome Research，2006，16(8):990－994．

［22］ Ishiguro N，Inoshima Y，Shigehara N．Mitochondrial DNA analysis of the Japanese wolf(Canis lupus hodophilax Temminck，1839)andcomparison with representative wolf and domestic dog haplotypes ［J］．Zoological Science，2009，26(11):765－770．

［23］ Verscheure S，Backeljau T，Desmyter S．Coding region SNP analysis to enhance dog mtDNA discrimination power in forensic casework ［J］．Forensic Science International:Genetics，2015，14:86－95．

［24］ Okumura N，Ishiguro N，Nakano M，et al．Intra-and interbreed genetic variations of mitochondrial DNA major non-coding regions in Japanese native dog breeds(Canis familiaris) ［J］．Animal Genetics，1996，27(6):397 －405．

［25］ Savolainen P，Zhang Yaping，Luo Jing，et al．Genetic evidence for an East Asian origin of domestic dogs［J］．Science，2002，298(5598):1610－1613．

［26］ 李隱俠，李齊發，劉振山，等．藏獒線粒體DNA D-loop區序列測定及其分類地位研究［J］．南京農業大學學報，2009，32(2):107－111．

［27］ 陳睿賾，張英杰，劉月琴．mtDNA分析在羊進化多樣性及種屬鑒別中的應用［J］．中國草食動物科學，2012(Z1):75－78．

［28］ 柳燕，李莉．線粒體 DNA cytb和 HVI的分析用于種屬鑒定［J］．中國司法鑒定，2006(4):14－17．

［29］ 卞曉東，張秀梅，肖永雙，等．線粒體DNA序列在沙氏下鱵魚魚卵鑒別上的應用 ［J］．中國海洋大學學報:自然科學版，2007，37(S1):111－116．

［30］ 趙亮，劉志瑾，任寶平，等．滇金絲猴線粒體D-loop片段SNPs位點的鑒別及初步分析［J］．獸類學報，2010，30(3):239－246．

［31］ Lalitha R，Chandavar V R．Intraspecific variations in Cyt b and D-loop sequences of Testudine species，Lissemys punctata from south Karnataka［J］．Journal of Advanced Research，2018，9:87 －95．

［32］ Kim M G，Kim N Y，Lee S S，et al．Phylogenetic relationships of jeju dogs to other domestic and foreign dog breeds determined by using mtDNA D-loop sequences［J］．Journal of Animal Science and Technology，2011，53(4):303－310．

［33］ Vilà C，Walker C，Sundqvist A K，et al．Combined use of maternal，paternal and bi-parental genetic markers for the identification of wolf-dog hybrids［J］．Heredity，2003，90(1):17－24．

圖1 基于mtDNA控制區單倍型構建的狼和狗系統發育樹Fig.1 Phylogenetic tree of mtDNA control region haplotypes in dogs and wolves

表1 狼和狗線粒體DNA控制區序列Tab.1 Mitochondrial DNA control region sequences of dogs and wolves

續表1

表2 狼和狗mtDNA控制區結構和序列多態性Tab.2 Structure and sequence polymorphisms in the mtDNA control region of dogs and wolves

表3 狼和狗單倍型頻率Tab.3 Haplotype frequencies of dogs and wolves