基于SNP芯片信息分析西藏四個(gè)本地牛種的血統(tǒng)組成

摘 要： 旨在探究西藏4個(gè)本地牛種的遺傳結(jié)構(gòu)和血統(tǒng)組成。本研究以22頭高原荷斯坦牛與20頭高原娟姍牛以及45頭美國(guó)婆羅門牛、46頭云南瘤牛和11頭婆羅門牛為參照標(biāo)準(zhǔn)品種，待測(cè)樣品為4個(gè)西藏本地品種的177頭牛。通過(guò)Illumina GGP100K芯片雜交獲取所有待測(cè)個(gè)體的SNPs數(shù)據(jù)，利用MEGA X和Admixture等軟件對(duì)SNPs進(jìn)行遺傳結(jié)構(gòu)與血統(tǒng)組成分析。結(jié)果表明，西藏牛較其他西藏本地牛種與高原荷斯坦牛和高原娟姍牛有更近的遺傳關(guān)系；部分樟木牛和阿沛甲咂牛將西藏牛分為兩個(gè)大群；日喀則駝峰牛與云南瘤牛遺傳關(guān)系更近。在日喀則駝峰牛群體中，高原牛血統(tǒng)比例超過(guò)90%，在阿沛甲咂牛與西藏牛群體中，高原牛血統(tǒng)也超過(guò)了50%，僅在樟木牛群體中，高原牛血統(tǒng)比例為33.15%。其次，在樟木牛群體中，瘤牛血統(tǒng)占比高達(dá)18.27%，在阿沛甲咂牛、西藏牛與日喀則駝峰牛群體中，瘤牛血統(tǒng)較少（5%）。在阿沛甲咂牛群體內(nèi)有高達(dá)35.94%的普通牛血統(tǒng)，在樟木牛群體內(nèi)，奶牛血統(tǒng)高達(dá)48.58%。本研究對(duì)西藏本地牛種今后的雜交擴(kuò)繁和橫交固定提供了重要指導(dǎo)作用，也為阿沛甲咂牛與樟木牛兩個(gè)瀕危品種的資源保護(hù)和利用提供了理論參考。

關(guān)鍵詞： 阿沛甲咂牛；樟木牛；西藏牛；日喀則駝峰牛；SNPs數(shù)據(jù)；遺傳結(jié)構(gòu)；血統(tǒng)組成

中圖分類號(hào)： S823.2

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：0366-6964（2024）10-4377-14

收稿日期：2024-04-07

基金項(xiàng)目：西藏自治區(qū)重點(diǎn)研發(fā)計(jì)劃（CGZH2023000257）；國(guó)家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系（CARS-36）；西藏自治區(qū)區(qū)域協(xié)同創(chuàng)新專項(xiàng)（QYXTZX-LS2021-01）

作者簡(jiǎn)介：馬龍剛（1996-），男，甘肅天水人，碩士，主要從事分子及數(shù)量遺傳學(xué)研究，E-mail： 13369452684@163.com

*通信作者：王雅春，主要從事分子及數(shù)量遺傳學(xué)研究，E-mail：wangyachun@cau.edu.cn

Study on the Genetic Composition of Four Local Tibetan Cattle Breeds Based on SNP Chip

Analysis

MA" Longgang1， LIU" Nan2， NI MA Qunzong2， WANG" Xin3， LU" Jian4， MAO" Huaming5， CHEN" Gong

6， DAN ZENG Ouzhu2， LA BA Ciren2， ZHANG" Nan2， YU" Fuqing4， WANG" Yachun1*

（1.College of Animal Science and Technology， China Agricultural University， Beijing 100193，

China;

2.Animal Disease Prevention and Control Center of Tibet Autonomous Region， Lhasa 850009，

China;

3.Beijing Dairy Cattle Center， Beijing 100192，" China;

4.National Animal Husbandry

Station， Beijing 100125，" China;

5.College of Animal Science and Technology ，Yunnan Agricultural

University， Kunming 650201， China;

6.Institute of Life Science and Bioengineering， Beijing JiaoTong

University， Beijing 100044， China）

Abstract：" The study aimed to explore the genetic structure and bloodline composition of 4 local cattle breedings in Tibet. In this study， 22 plateau Holstein cattle， 20 plateau Jersey cattle， 45 American Brahmin cattle， 46 Yunnan Zebu and 11 Brahmen cattle as reference breeds， the samples to be tested were 177 cattle of 4 local breeds in Tibet. SNPs data were obtained by Illumina GGP100K chip hybridization， and the genetic structure and bloodline composition of SNPs were analyzed by MEGA X and Admixture. The results showed that Tibetan cattle had a closer genetic relationship with plateau Holstein cattle and plateau Jersey cattle than other local cattle breeds in Tibet. Some Zhangmu cattle and Apeijia cattle divided Tibetan cattle into two large subpopulations; There was a close genetic relationship between Shigatse hump cattle and Yunnan Zebu. In the Shigatse hump cattle population， the proportion of plateau cattle bloodline was more than 90%， in the Apaijia cattle population and Tibet cattle population， the proportion of plateau cattle bloodline was more than 50%， only in the Zhangmu cattle population， the proportion of plateau cattle bloodline was 33.15%. Secondly， bos indicus bloodline accounted for 18.27% in Zhangmu cattle population， while bos indicus bloodline was less （ 5%） in Apaijiza cattle， Tibet cattle and Shigatse hump cattle population. There were 35.94% dairy cattle bloodline in Apaijia population and 48.58% dairy cattle bloodline in Zhangmu population. This study provided important guidance for the future cross-breeding and cross-fixing of the native cattle species in Tibet， and also provided theoretical reference for the conservation and utilization of the resources of the two endangered breeds of Apaijiza cattle and Zhangmu cattle.

Key words： Apejiaza cattle; Zhangmu cattle; Tibetan cattle; Shigatse hump cattle; SNPs data; genetic structure; bloodline composition

*Corresponding author： WANG Yachun， E-mail：wangyachun@cau.edu.cn

西藏地處高原地區(qū)，其獨(dú)特的高寒生態(tài)環(huán)境和嚴(yán)苛的氣候條件為畜牧業(yè)的發(fā)展提出了獨(dú)特的挑戰(zhàn)［1］。在這獨(dú)特的生存條件下，西藏本地黃牛品種（阿沛甲咂牛、樟木牛、日喀則駝峰牛和西藏牛）成為研究的焦點(diǎn)之一［2］。這些品種以其獨(dú)特且出色的高原適應(yīng)性、抗寒性和粗放飼養(yǎng)條件下的生存能力而備受關(guān)注。在世界珍稀地方黃牛品種中，西藏本地黃牛因其獨(dú)特的遺傳特征和生態(tài)適應(yīng)性（高寒、耐寒及耐粗飼）而引起了廣泛的科研興趣［2］。特別是含有瘤牛血統(tǒng)的西藏本地黃牛，其獨(dú)特的遺傳基因使其在高寒地區(qū)表現(xiàn)出卓越的生存和繁衍能力。這些品種在高原地區(qū)簡(jiǎn)易的飼養(yǎng)條件下能夠保持較高的產(chǎn)乳產(chǎn)肉役用性能，為當(dāng)?shù)匦竽翗I(yè)的可持續(xù)發(fā)展提供了有力的支持。西藏本地黃牛品種不僅是當(dāng)?shù)匦竽翗I(yè)的重要資源，更是全球研究高原適應(yīng)性的重要材料［3］。通過(guò)深入研究西藏本地黃牛的遺傳背景、生理機(jī)制以及在不同環(huán)境條件下的表現(xiàn)，可以為培育更適應(yīng)高原環(huán)境的新型牛品種提供理論依據(jù)和實(shí)踐經(jīng)驗(yàn)。

西藏本地牛作為珍稀的地方品種，其血統(tǒng)組成包含了豐富的遺傳信息［4］。研究西藏本地牛的血統(tǒng)組成有助于保護(hù)該地區(qū)畜牧資源的遺傳多樣性。其適應(yīng)高原的特殊基因可用于培育新品種，提高生產(chǎn)性能和經(jīng)濟(jì)效益。此外，有益基因的了解可制定更有效的育種策略，推動(dòng)全球高原適應(yīng)性研究與應(yīng)用，促進(jìn)農(nóng)牧民的經(jīng)濟(jì)增長(zhǎng)。有研究報(bào)道，血統(tǒng)組成對(duì)于保護(hù)遺傳多樣性［5-6］、挖掘高原適應(yīng)性基因［6］、推動(dòng)畜牧業(yè)可持續(xù)發(fā)展［7］以及為全球高寒地區(qū)養(yǎng)殖提供經(jīng)驗(yàn)教訓(xùn)都具有深遠(yuǎn)的意義。

本研究根據(jù)阿沛甲咂牛、樟木牛、日喀則駝峰牛和西藏牛的基因芯片數(shù)據(jù)，以適應(yīng)高海拔地區(qū)環(huán)境的荷斯坦牛、娟姍牛以及云南地區(qū)的云南瘤牛和婆羅門牛和下載數(shù)據(jù)美國(guó)婆羅門牛為參照，分析177頭西藏本地牛的基因組結(jié)構(gòu)和血統(tǒng)組成，為今后保護(hù)遺傳多樣性、挖掘高原適應(yīng)性基因、推動(dòng)畜牧業(yè)可持續(xù)發(fā)展提供基礎(chǔ)。

1 材料與方法

1.1 試驗(yàn)材料

在西藏拉薩市、林芝市工布江達(dá)縣、日喀則市和聶拉木縣采集西藏本地牛種177頭，其中包括西藏牛61頭、阿沛甲咂牛39頭、日喀則駝峰牛57頭和樟木牛20頭。其中阿沛甲咂牛與樟木牛為該群體的全部個(gè)體。同時(shí)在拉薩市采集高原荷斯坦牛22頭［8］，高原娟姍牛20頭［8］和云南地區(qū)采集的云南瘤牛46頭和婆羅門牛11頭的血液樣品。同時(shí)在WIDDE［9］上獲取了45頭美國(guó)婆羅門牛的770K芯片數(shù)據(jù)（參考基因組UMD3.1）并使用liftOver （https：∥genome. ucsc.edu/cgi-bin/hgLiftOver） 進(jìn)行基因組版本轉(zhuǎn)換最終獲得參考基因組版本為ARS_UCD1.2的SNP數(shù)據(jù)［9］，本研究后續(xù)分析使用9個(gè)品種的321頭牛，其中荷斯坦牛與娟姍牛均在西藏拉薩市克服高原反應(yīng)存活并繁殖超過(guò)2代，云南瘤牛、美國(guó)婆羅門牛與婆羅門牛作為瘤牛血統(tǒng)參考群體（樣品信息詳情見表1）。在后續(xù)的分析中，荷斯坦牛與娟姍牛分組為奶牛血統(tǒng)（NN）；云南瘤牛、美國(guó)婆羅門牛與婆羅門牛分組為瘤牛血統(tǒng)（LN）；西藏牛、阿沛甲咂牛、日喀則駝峰牛和樟木牛分組為高原牛血統(tǒng)（XZ）。

1.2 芯片數(shù)據(jù)測(cè)定

從血樣中提取基因組DNA，用 Illumina GGP100K 進(jìn)行測(cè)序，獲得原始數(shù)據(jù)95 256個(gè)SNPs位點(diǎn)信息。使用plink［10］軟件過(guò)濾掉最小等位基因頻率小于0.05、檢測(cè)率低于90% 的SNPs位點(diǎn)。得到的184頭西藏本地牛種的基因型數(shù)據(jù)與參考群體的基因型數(shù)據(jù)進(jìn)行合并，取共有的SNPs位點(diǎn)進(jìn)行后續(xù)分析。

1.3 遺傳結(jié)構(gòu)與血統(tǒng)組成分析

合并后的276個(gè)樣品的基因型數(shù)據(jù)使用plink［10］進(jìn)行主成分分析（PCA，--pca），同時(shí)利用MEGA X［11］軟件，使用鄰接法（Neighbor Joining，NJ）［12］且基于最大復(fù)合似然法［13］得到的遺傳距離矩陣對(duì)樣本進(jìn)行進(jìn)化樹分析。最后使用Admixture［14］軟件對(duì)這58頭待測(cè)樣本的血統(tǒng)成分進(jìn)行推斷。所有結(jié)果的可視化均通過(guò)R語(yǔ)言實(shí)現(xiàn)。

2 結(jié) 果

2.1 遺傳結(jié)構(gòu)分析

將SNP基因型數(shù)據(jù)使用plink軟件質(zhì)控后進(jìn)行主成分分析（PCA），8個(gè)品種276個(gè)樣品經(jīng)質(zhì)控后獲得81 784個(gè)SNPs用于后續(xù)分析。PCA結(jié)果（圖1A）發(fā)現(xiàn)，高原荷斯坦牛（HO）、高原娟姍牛（JS）在第一主成分（PC1）、第二主成分（PC2）和第三主成分（PC3）明顯與西藏牛（XZ）、阿沛甲咂牛牛（AP）、日喀則駝峰牛（TF）和樟木牛（ZM）區(qū)分開，云南瘤牛（LN）和婆羅門牛（PLM）在第一主成分（PC1）、第二主成分（PC2）與西藏4個(gè)本地牛種區(qū)分較為明顯。而西藏4個(gè)本地牛品種在第三主成分（PC3）與西藏4個(gè)本地牛種并未形成較好的差異。揭示了西藏4個(gè)本地牛種與LN和PLM之間存在的關(guān)系較近。在西藏4個(gè)本地牛種中，日喀則駝峰牛（TF）品種一致性較強(qiáng)，整體呈現(xiàn)較為緊湊的聚集狀態(tài)，阿沛甲咂牛（AP）、西藏牛（XZ）和樟木牛（ZM）之間關(guān)系較為復(fù)雜，尤其是ZM與XZ之間存在著較大程度的重疊。表明現(xiàn)階段存在的ZM群體中，XZ群體參與了較為頻繁的雜交過(guò)程。AP在一定程度上表現(xiàn)出了與TF之間較大的差異性以及與XZ之間明顯的相似性與差異性。ZM暫時(shí)未表現(xiàn)出較大的差異性，但在一定程度上與HO和JS表現(xiàn)出了一定程度的相似性。

結(jié)合圖1D可見奶牛組、瘤牛組和高原牛組之間的關(guān)系，PC1解釋了3個(gè)群組的42.51%的變異，PC2解釋了10.2%的變異，表明不同的群組之間存在一定的差異，結(jié)合圖1D、E和F可見，高原群體更偏向于瘤牛群體。

由圖1A、B、C可見，西藏4個(gè)本地品種在一定程度上受LN與PLM的影響。表現(xiàn)出一定的相似性。這在外形上也有一定表現(xiàn)，例如：西藏4個(gè)本地品種或多或少存在一定的肩峰與垂皮，同時(shí)ZM生存的環(huán)境也是高濕熱地區(qū)，與LN和PLM生存環(huán)境存在一定的相似性。

進(jìn)化樹的結(jié)果表明（圖2），西藏牛、樟木牛與荷斯坦牛和娟姍牛有著更近的遺傳關(guān)系。同時(shí)阿沛甲咂牛與樟木牛將西藏牛分為兩個(gè)大群，表明阿沛甲咂牛、樟木牛和西藏牛有著緊密的遺傳關(guān)系，日喀則駝峰牛與瘤牛有著更近的遺傳關(guān)系，表現(xiàn)在外貌上都具有垂皮與肩峰。根據(jù)進(jìn)化樹的結(jié)果而言，阿沛甲咂牛與日喀則駝峰牛有著更好的群體一致性，樟木牛與西藏牛群體一致性較差，群體內(nèi)個(gè)體差異較大。

2.2 血統(tǒng)組成分析

根據(jù)養(yǎng)殖人員的介紹以及西藏本地牛種的記錄，本研究使用了在高原地區(qū)存活并繁殖超過(guò)2代的荷斯坦牛與娟姍牛以及云南地區(qū)的瘤牛和婆羅門牛4個(gè)品種作為參考群體，分析西藏本地牛種的血統(tǒng)組成。

結(jié)合圖3與表2-表6結(jié)果表明，在阿沛甲咂牛群體的血統(tǒng)中，約59.25%的高原牛血統(tǒng)、4.81%的瘤牛血統(tǒng)以及35.94%的普通牛血統(tǒng)。在樟木牛群體的血統(tǒng)中，約33.15%的高原牛血統(tǒng)、約48.58%的普通牛血統(tǒng)及約18.27%的瘤牛血統(tǒng)。在西藏牛群體的血統(tǒng)中，約53.41%的高原牛血統(tǒng)、43.88%的普通牛血統(tǒng)及2.71%的瘤牛血統(tǒng)。在日喀則駝峰牛群體的血統(tǒng)中，約90.53%的高原牛血統(tǒng)、1.87%瘤牛血統(tǒng)及7.60%的普通牛血統(tǒng)。每個(gè)品種的個(gè)體具體的血統(tǒng)組成見表6。該分析與PCA和進(jìn)化樹的結(jié)果基本一致，既西藏4個(gè)本地品種之間具有一定的相似性，具有不同比例高原牛血統(tǒng)且占比較大；除樟木牛外，其他3個(gè)品種的瘤牛血統(tǒng)占比幾乎可以忽略不計(jì)（低于5%）；在阿沛甲咂牛與阿沛甲咂牛群體內(nèi)，有一定比例的普通牛血統(tǒng)。

3 討 論

西藏本地牛種在高原地區(qū)表現(xiàn)出出色的高寒適應(yīng)性和抗寒性，能夠在極端的氣候條件下生存，并適應(yīng)粗放飼養(yǎng)條件［15］。研究其生理機(jī)制和遺傳特征有助于培育更適應(yīng)高寒環(huán)境的新型牛品種，提高牛的生產(chǎn)性能，為高原地區(qū)的畜牧業(yè)發(fā)展提供重要支持［16］。西藏本地牛種具有豐富的遺傳多樣性和獨(dú)特的生態(tài)適應(yīng)性［17］。通過(guò)研究西藏本地牛種，可以更好地保護(hù)和合理利用地方畜種資源，防止基因資源的喪失和退化［18］。作為高原地區(qū)的典型代表，牦牛和西藏本地牛種研究成果對(duì)其他高寒地區(qū)的畜牧業(yè)具有啟示作用［19］。通過(guò)了解西藏本地牛種在高寒環(huán)境下的適應(yīng)性和生態(tài)影響，可以為其他類似地區(qū)提供經(jīng)驗(yàn)教訓(xùn)，促進(jìn)全球高寒地區(qū)畜牧業(yè)的可持續(xù)發(fā)展［20］。Zhang等［21］利用拷貝數(shù)變異證明了分布在5個(gè)不同氣候帶的中國(guó)本地牛具有表型多樣性和復(fù)雜的環(huán)境適應(yīng)性。例如，西藏牛、阿沛甲咂牛和日喀則駝峰牛棲息在青藏高原，這些牛對(duì)惡劣的氣候條件（如低氣壓和低氧）表現(xiàn)出極強(qiáng)的耐受性［17，21］。中國(guó)牛品種具有豐富的遺傳資源，為了解拷貝數(shù)變異在適應(yīng)特定環(huán)境，特別是高海拔地區(qū)、生存、維持健康和繁殖的作用［21-22］。

認(rèn)識(shí)和了解西藏本地牛種的血統(tǒng)組成是進(jìn)一步開展瀕危牛種保護(hù)及提高畜牧水平進(jìn)行雜交改良［23］與橫交固定［24］選育的基礎(chǔ)，但由于牧民缺乏專業(yè)指導(dǎo)且多數(shù)為散養(yǎng)，維持著隨機(jī)交配的擴(kuò)繁模式導(dǎo)致絕大多數(shù)的西藏本地牛種無(wú)配種和系譜記錄，最終導(dǎo)致樟木牛全群存活數(shù)量不足，如阿沛甲咂牛公牛僅為1頭，無(wú)法從繁殖記錄中推導(dǎo)計(jì)算其具體的血統(tǒng)組成比例。因此，本研究選擇使用基因組范圍內(nèi)的SNP數(shù)據(jù)來(lái)對(duì)各品種全群個(gè)體進(jìn)行血緣溯源分析。

在4個(gè)品種的血統(tǒng)組成分析中，可以發(fā)現(xiàn)阿沛甲咂牛群體的普通牛血統(tǒng)比例最高，結(jié)合該群體的養(yǎng)殖模式，推測(cè)由于當(dāng)?shù)仫曫B(yǎng)模式——僅保留1頭公牛用于繁殖其余公牛出售，導(dǎo)致群體的選擇強(qiáng)度較高，同時(shí)該群體主要用于擠奶供給制作相關(guān)奶制品導(dǎo)致奶牛雜交的血統(tǒng)比例升高（奶牛參考荷斯坦牛與娟姍牛），高原牛血統(tǒng)的高比例影響該群體在高原地區(qū)的生存與繁殖；由于歷史原因，導(dǎo)致該品種血統(tǒng)內(nèi)存在一定的瘤牛血統(tǒng)。樟木牛群體生存環(huán)境不同于其他3個(gè)品種，當(dāng)?shù)貧夂驖駶?rùn)但仍屬于高原地區(qū)，其海拔較拉薩市與林芝市低。特殊的氣候環(huán)境導(dǎo)致樟木牛血統(tǒng)組成中高原牛血統(tǒng)比例在4個(gè)高原牛種中最低。日喀則駝峰牛屬于西藏地區(qū)的牛種，其起源受到印度、尼泊爾、不丹等地引入的瘤牛公牛的影響，擁有部分瘤牛血統(tǒng)［25］，這與本研究結(jié)果相似，同時(shí)本研究發(fā)現(xiàn)瘤牛的血統(tǒng)比例在下降。日喀則駝峰牛是4個(gè)品種中唯一仍然役用的品種，導(dǎo)致該群體中出現(xiàn)較多的閹牛。由于存在重體力勞動(dòng)，該品種個(gè)體需要更大比例的高原牛血統(tǒng)用于在高原地區(qū)從事體力勞動(dòng)；由于生存環(huán)境較為干燥及幾乎不用于擠奶導(dǎo)致普通牛血統(tǒng)與瘤牛血統(tǒng)較低。劉詩(shī)平［26］報(bào)道了阿沛甲咂牛與樟木牛的重要性及解決瀕危的重要措施。

卜李那等［27］基于基因組SNPs數(shù)據(jù)，利用Admixture和RFMix軟件分析了4個(gè)從市場(chǎng)選取的豬肌肉樣品的遺傳成分，可以估計(jì)每個(gè)個(gè)體的遺傳成分來(lái)源，為品牌豬肉鑒定和品牌保護(hù)提供參考。鐘海安等［28］利用基因芯片數(shù)據(jù)，用西門塔爾牛和荷斯坦牛作為參照群分析了58頭藏黃牛（西藏牛）的遺傳結(jié)構(gòu)與血統(tǒng)組成，證明了雜交藏黃牛（西藏牛）的血統(tǒng)中，荷斯坦奶牛占比最高為 50%，其次是本地藏黃牛（高原牛）約 40%，這與本研究的結(jié)果基本一致。宋科林等［29］也選擇使用SNP芯片分析徽縣青泥黑豬遺傳多樣性和遺傳結(jié)構(gòu)，同樣證明了SNP芯片分析遺傳結(jié)構(gòu)估計(jì)血統(tǒng)組成的準(zhǔn)確性，為本研究的結(jié)果提供了支持；Quan等［30］利用1 163頭中國(guó)本地牛研究遺傳結(jié)構(gòu)與馴化歷史，發(fā)現(xiàn)早勝牛、秦川牛和平?jīng)雠?duì)本土牛遺傳資源貢獻(xiàn)最大，建議優(yōu)先保護(hù)，這一觀點(diǎn)證明了探究西藏本地牛血統(tǒng)組成的必要性。Edea等［31］使用69 903個(gè)SNPs的數(shù)據(jù)，分析亞洲和非洲本地牛的遺傳結(jié)構(gòu)，促進(jìn)其改良與保護(hù)。馬浩然等［32］通過(guò)10個(gè)肉牛品種的SNP數(shù)據(jù)，將西藏牛歸入北方牛群體，表明西藏牛與延黃牛和蒙古牛的遺傳關(guān)系更近，與本研究PCA結(jié)果相似，進(jìn)一步支持了西藏牛與北方牛群體的密切遺傳聯(lián)系［33］。

本研究利用基因芯片獲得參考群體（荷斯坦牛、娟姍牛、美國(guó)婆羅門牛、瘤牛和云南婆羅門牛）和待測(cè)群體樣品（阿沛甲咂牛、樟木牛、西藏牛和日喀則駝峰牛）的基因組芯片數(shù)據(jù)，利用Admixture分析基因組結(jié)構(gòu)，估計(jì)出4個(gè)西藏本地牛品種177個(gè)待測(cè)個(gè)體的血統(tǒng)組成比例，發(fā)現(xiàn)這些西藏本地牛品種中，日喀則駝峰牛高原牛血統(tǒng)含量最高，瘤牛血統(tǒng)出現(xiàn)在除樟木牛外的3個(gè)品種中，普通牛血統(tǒng)也不同程度的出現(xiàn)在阿沛甲咂牛、樟木牛和日喀則駝峰牛群體中。本研究結(jié)果符合西藏4個(gè)本地牛種的培育歷史和雜交牛的外貌特征。研究西藏本地黃牛意義重大，主要體現(xiàn)在以下幾個(gè)方面：保護(hù)珍貴的基因資源［34］，防止遺傳多樣性喪失［35］；揭示其獨(dú)特的生態(tài)適應(yīng)機(jī)制［36］，為高原及極端環(huán)境畜牧業(yè)提供借鑒［36］；提升其肉質(zhì)、奶質(zhì)的營(yíng)養(yǎng)價(jià)值，促進(jìn)當(dāng)?shù)亟?jīng)濟(jì)發(fā)展［37］；傳承和保護(hù)西藏獨(dú)特的農(nóng)牧文化［36］；探索其天然疾病抵抗力，開發(fā)新的疫苗和疾病防治策略［38-39］。通過(guò)研究西藏本地黃?？梢詫?shí)現(xiàn)科學(xué)、經(jīng)濟(jì)、生態(tài)和文化的多重效益，為高原畜牧業(yè)的可持續(xù)發(fā)展提供支持［40-41］。

4 結(jié) 論

本研究對(duì)西藏4個(gè)本地品種（阿沛甲咂牛、樟木牛、西藏牛和日喀則駝峰牛）的177個(gè)待測(cè)個(gè)體的遺傳結(jié)構(gòu)與血統(tǒng)組成進(jìn)行了分析。結(jié)果發(fā)現(xiàn)，西藏牛較其他西藏本地牛種與荷斯坦牛和娟姍牛有更近的遺傳關(guān)系；日喀則駝峰牛和樟木牛與瘤牛具有更近的遺傳關(guān)系。但是在血統(tǒng)分析中，日喀則駝峰牛群體具有最高的高原牛血統(tǒng)比例（91%），在阿沛甲咂牛與西藏牛群體中高原牛血統(tǒng)超過(guò)了50%。在阿沛甲咂牛、西藏牛與日喀則駝峰牛群體中，瘤牛血統(tǒng)較少（5%左右）。本研究結(jié)果可以為西藏本地牛種今后的雜交擴(kuò)繁和橫交固定提供重要指導(dǎo)作用，也為阿沛甲咂牛與樟木牛兩個(gè)瀕危品種的資源保護(hù)和利用提供理論參考。

參考文獻(xiàn)（References）：

［1］ 張核真，路紅亞，洪健昌，等.藏西北地區(qū)氣候變化及其對(duì)草地畜牧業(yè)的影響［J］.干旱區(qū)研究，2013，30（2）：308-314.

ZHANG H Z，LU H Y，HONG J C，et al.Climate change and its effect on steppe animal husbandry in Northwest Tibet［J］.Arid Zone Research，2013，30（2）：308-314.（in Chinese）

［2］ 鄧 兵，晉美旺杰，任增幫，等.基于RAD-seq簡(jiǎn)化基因組測(cè)序的西藏黃牛遺傳多樣性研究［J］.江蘇農(nóng)業(yè)科學(xué)，2021， 49（16）：153-157.

DENG B，JINMEI W J，REN Z B，et al.Study on genetic diversity of Tibetan cattle based on RAD-seq simplified genome sequencing［J］.Jiangsu Agricultural Sciences，2021，49（16）：153-157.（in Chinese）

［3］ 鄭玉才.淺析西藏山南市地方土種高原黃牛資源保護(hù)［J］.畜禽業(yè)，2019，30（10）：51.

ZHENG Y C.Analysis on the protection of local soil species plateau yellow cattle resources in Shannan of Tibetan［J］.Livestock and Poultry Industry，2019，30（10）：51.（in Chinese）

［4］ 唐建華，陳曉英，宋天增，等.西藏黃牛種質(zhì)資源保護(hù)與利用研究［J］.中國(guó)牛業(yè)科學(xué)，2016，42（3）：48-51.

TANG J H，CHEN X Y，SONG T Z，et al.The research on conservation and utility of Tibet Cattle germplasm resources［J］.China Cattle Science，2016，42（3）：48-51.（in Chinese）

［5］ 黃 興，柴志欣，信金偉，等.西藏牛亞科部分群體線粒體DNA遺傳多樣性研究［J］.西南民族大學(xué)學(xué)報(bào)（自然科學(xué)版），2019，45（2）：117-124.

HUANG X，CHAI Z X，XIN J W，et al.Genetic diversity of mitochondrial DNA of some subfamilies of cattle in Tibet［J］.Journal of Southwest Minzu University （Natural Science Edition），2019，45（2）：117-124.（in Chinese）

［6］ 張桂香，鄭友民，王志剛，等.我國(guó)部分黃牛品種線粒體D-loop區(qū)遺傳多樣性與起源分化［J］.遺傳，2009，31（2）：160-168.

ZHANG G X，ZHENG Y M，WANG Z G，et al.Genetic diversity and origin of mitochondria DNA D-loop region of some Chinese indigenous cattle breeds［J］.HEREDITAS（Beijing），2009，31（2）：160-168.（in Chinese）

［7］ 劉錫武，劉 鵬，高玉君.西藏日喀則市黃牛改良現(xiàn)狀與對(duì)策［J］.中國(guó)牛業(yè)科學(xué)，2009，35（4）：63-65.

LIU X W，LIU P，GAO Y J.Current situation and measures to improve yellow cattle in Shigatse of Tibet［J］.Chinese Cattle Science，2009，35（4）：63-65.（in Chinese）

［8］ 黃上真，馬龍剛，婁文琦，等.高原地區(qū)奶牛血液指標(biāo)的影響因素分析［J］.畜牧獸醫(yī)學(xué)報(bào)，2023，54（5）：1964-1978.

HUANG S Z，MA L G，LOU W Q，et al.Analysis of influencing factors on blood indicators of dairy cows at high-altitude area［J］.Acta Veterinaria et Zootechnica Sinica，2023，54（5）：1964-1978.（in Chinese）

［9］ SEMPéRé G，MOAZAMI-GOUDARZI K，EGGEN A，et al.WIDDE：a Web-Interfaced next generation database for genetic diversity exploration，with a first application in cattle［J］.BMC Genomics，2015，16：940.

［10］ PURCELL S，NEALE B，TODD-BROWN K，et al.PLINK：a tool set for whole-genome association and population-based linkage analyses［J］.Am J Hum Genet，2007，81（3）：559-575.

［11］ KUMAR S，STECHER G，LI M，et al.MEGA X：molecular evolutionary genetics analysis across computing platforms［J］.Mol Biol Evol，2018，35（6）：1547-1549.

［12］ ZHANG W，SUN Z R.Random local neighbor joining：a new method for reconstructing phylogenetic trees［J］.Mol Phylogenet Evol，2008，47（1）：117-128.

［13］ TAMURA K，NEI M，KUMAR S.Prospects for inferring very large phylogenies by using the neighbor-joining method［J］.Proc Natl Acad Sci U S A，2004，101（30）：11030-11035.

［14］ ALEXANDER D H，NOVEMBRE J，LANGE K.Fast model-based estimation of ancestry in unrelated individuals［J］.Genome Res，2009，19（9）：1655-1664.

［15］ KONG Z W，ZHOU C S，LI B，et al.Integrative plasma proteomic and microRNA analysis of Jersey cattle in response to high-altitude hypoxia［J］.J Dairy Sci，2019，102（5）：4606-4618.

［16］ 張劍搏，丁學(xué)智，Anum A A，等.高原土著動(dòng)物適應(yīng)性進(jìn)化的研究進(jìn)展［J］.畜牧獸醫(yī)學(xué)報(bào)，2019，50（9）：1723-1736.

ZHANG J B，DING X Z，AHMAD A A，et al.Advances in research on adaptive evolution of native animals of tibetan plateau［J］.Acta Veterinaria et Zootechnica Sinica，2019，50（9）：1723-1736.（in Chinese）

［17］ CHEN N B，XIA X T，HANIF Q，et al.Global genetic diversity，introgression，and evolutionary adaptation of indicine cattle revealed by whole genome sequencing［J］.Nat Commun，2023，14（1）：7803.

［18］ 賈功雪，丁路明，徐尚榮，等.青藏高原牦牛遺傳資源保護(hù)和利用：?jiǎn)栴}與展望［J］.生態(tài)學(xué)報(bào)，2020，40（18）：6314-6323.

JIA G X，DING L M，XU S R，et al.Conservation and utilization of yak genetic resources in Qinghai-Tibet plateau：problems and perspectives［J］.Acta Ecologica Sinica，2020，40（18）：6314-6323.（in Chinese）

［19］ 馬志杰.牦?；蚪M和轉(zhuǎn)錄組研究進(jìn)展［J］.中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào)，2019，24（7）：79-86.

MA Z J.Research progress of the genome and transcriptome of yak［J］.Journal of China Agricultural University，2019，24（7）：79-86.（in Chinese）

［20］ QIU Q，ZHANG G J，MA T，et al.The yak genome and adaptation to life at high altitude［J］.Nat Genet，2012，44（8）：946-949.

［21］ ZHANG Y R，HU Y，WANG X G，et al.Population structure，and selection signatures underlying high-altitude adaptation inferred from genome-wide copy number variations in Chinese indigenous cattle［J］.Front Genet，2020，10：1404.

［22］ 干珠扎布，胡國(guó)錚，高清竹，等.藏北高原草地生態(tài)治理與畜牧業(yè)協(xié)同發(fā)展模式研究［J］.中國(guó)工程科學(xué)，2019，21（5）：93-98.

HASBAGAN G，HU G Z，GAO Q Z，et al.Synergetic mode for grassland ecological management and animal husbandry development in Northern Tibet Plateau［J］.Strategic Study of CAE，2019，21（5）：93-98.（in Chinese）

［23］ 王 莉，旦增洛桑，馬金英，等.娟姍牛雜交改良西藏本地黃牛的效果調(diào)查［J］.中國(guó)奶牛，2015（2）：13-15.

WANG L，DANZENG L S，MA J Y，et al.Hybrid improvement of Jersey cattle crossbreeding with Tibet local yellow cattle［J］.China Dairy Cattle，2015（2）：13-15.（in Chinese）

［24］ 張國(guó)慶，郭妮妮，胡雄貴，等.山南市黃牛改良橫交固定研究及推廣應(yīng)用［J］.現(xiàn)代農(nóng)業(yè)科技，2016（24）：248，251.

ZHANG G Q，GUO N N，HU X G，et al.Study and application of improved crossbreed fixation of cattle in Shannan City［J］.Modern Agricultural Science and Technology，2016（24）：248，251.（in Chinese）

［25］ 楊章平，常 洪，李相運(yùn)，等.西藏南部高寒地區(qū)瘤牛群體形態(tài)及生態(tài)特征研究［J］.揚(yáng)州大學(xué)學(xué)報(bào)：農(nóng)業(yè)與生命科學(xué)版，2002，23（1）：23-26.

YANG Z P，CHANG H，LI X Y，et al.Study on morphology and ecology of the indicus populations from cooled region of South Tibet［J］.Journal of Yangzhou University： Agricultural and Life Science Edition，2002，23（1）：23-26.（in Chinese）

［26］ 劉詩(shī)平.絕不讓雪域高原瀕危種質(zhì)資源滅失——西藏阿沛甲咂牛與樟木牛搶救性保種全面展開［J］.中國(guó)畜牧業(yè)，2022（24）：10-11.

LIU S P.Never let the endangered germplasm resources be lost in the snowy plateau-the rescue and conservation of Jiaza cattle and Zhangmu cattle in Tibet are carried out in an all-round way［J］.China Animal Industry，2022（24）：10-11.（in Chinese）

［27］ 卜李那，彭業(yè)博，張春媛，等.基于祖源分析的豬遺傳成分鑒定［J］.中國(guó)畜牧獸醫(yī)，2020，47（2）：544-553.

BU L N，PENG Y B，ZHANG C Y，et al.Ancestry analysis of genetic components of pigs［J］.China Animal Husbandry amp; Veterinary Medicine，2020，47（2）：544-553.（in Chinese）

［28］ 鐘海安，張 健，旦增洛桑，等.基于基因芯片分析西藏山南雜交藏黃牛血統(tǒng)組成［J］.中國(guó)畜牧雜志，2023，59（11）：166-169，174.

ZHONG H A，ZHANG J，DANZENG L S，et al.Analysis of pedigree composition of Shannan hybrid Tibetan yellow cattle in Tibet based on gene chip［J］.Chinese Journal of Animal Science，2023，59（11）：166-169，174.（in Chinese）

［29］ 宋科林，閆尊強(qiáng)，王鵬飛，等.基于SNP芯片分析徽縣青泥黑豬遺傳多樣性和遺傳結(jié)構(gòu)［J］.畜牧獸醫(yī)學(xué)報(bào)，2024，55（3）：995-1006.

SONG K L，YAN Z Q，WANG P F，et al.Analysis on genetic diversity and genetic structure based on SNP chips of Huixian Qingni black pig［J］.Acta Veterinaria et Zootechnica Sinica，2024，55（3）：995-1006.（in Chinese）

［30］ QUAN J Q，LI Y Q，YANG Y H，et al.Population genetic diversity and genetic evaluation models reveal the maternal genetic structure and conservation priority characteristics of indigenous cattle in China［J］.Global Ecol Conserv，2021，32：e1903.

［31］ EDEA Z，BHUIYAN M S A，DESSIE T，et al.Genome-wide genetic diversity，population structure and admixture analysis in African and Asian cattle breeds［J］.Animal，2015，9（2）：218-226.

［32］ 馬浩然，張路培，金生云，等.利用高密度SNP芯片評(píng)估中國(guó)地方肉牛品種基因組親緣關(guān)系［J］.畜牧獸醫(yī)學(xué)報(bào)，2023，54（10）： 4174-4185.

MA H R，ZHANG L P，JIN S Y，et al.Assessment of the genomic relationships for chinese indigenous beef cattle using high-density SNP chip［J］.Acta Veterinaria et Zootechnica Sinica，2023，54（10）：4174-4185.（in Chinese）

［33］ LYU Y，WANG F W，CHENG H J，et al.Recent selection and introgression facilitated high-altitude adaptation in cattle［J］.Sci Bull （Beijing），2024：S2095-9273（24）00380-3.

［34］ CHEN N B，ZHANG Z W，HOU J W，et al.Evidence for early domestic yak，taurine cattle，and their hybrids on the Tibetan Plateau［J］.Sci Adv，2023，9（50）：eadi6857.

［35］ CHEN N B，CAI Y D，CHEN Q M，et al.Whole-genome resequencing reveals world-wide ancestry and adaptive introgression events of domesticated cattle in East Asia［J］.Nat Commun，2018，9（1）：2337.

［36］ GAO X，WANG S，WANG Y F，et al.Long read genome assemblies complemented by single cell RNA-sequencing reveal genetic and cellular mechanisms underlying the adaptive evolution of yak［J］.Nat Commun，2022，13（1）：4887.

［37］ 增旦旺姆.西藏地區(qū)草地畜牧業(yè)發(fā)展現(xiàn)狀與重點(diǎn)任務(wù)［J］.畜禽業(yè)，2022，33（5）：82-83，86.

ZENGDAN W M.Development status and key tasks of grassland animal husbandry in Tibet［J］.Livestock and Poultry Industry， 2022，33（5）：82-83，86.（in Chinese）

［38］ 阿旺澤珍.西藏畜牧獸醫(yī)工作中動(dòng)物檢疫的思考［J］.畜禽業(yè)，2022，33（1）：77-78.

AWANG Z Z.Current situation and key tasks of grassland animal husbandry in Tibet：reflections on animal quarantine in animal husbandry and veterinary work in Tibet［J］.Livestock and Poultry Industry，2022，33（1）：77-78.（in Chinese）

［39］ 張春來(lái)，馬偉元.西藏畜牧獸醫(yī)工作中動(dòng)物檢疫現(xiàn)狀淺析［J］.吉林畜牧獸醫(yī)，2021，42（4）：128.

ZHANG C L，MA W Y.Analysis of animal quarantine status in animal husbandry and veterinary work in Tibet［J］.Jilin Animal Husbandry and Veterinary，2021，42（4）：128.（in Chinese）

［40］ 梁偉岸，王 平.農(nóng)業(yè)強(qiáng)國(guó)視域下中國(guó)畜牧業(yè)循環(huán)發(fā)展路徑探究［J/OL］.飼料研究，2024（12）：193-196［2024-07-18］.https：∥doi.org/10.13557/j.cnki.issn1002-2813.2024.12.034.

LIANG W A，WANG P.Exploration of circular development path of Chinese livestock industry from the perspective of an agricultural powerhouse［J/OL］.Feed Research，2024（12）：193-196［2024-07-18］.https：∥doi.org/10.13557/j.cnki.issn1002-2813.2024.12. 034.（in Chinese）

［41］ 鄒 驊，劉 沈.畜牧業(yè)數(shù)智化發(fā)展的現(xiàn)狀及對(duì)策建議［J］.畜牧產(chǎn)業(yè)，2024（6）：59-63.

ZOU H，LIU S.Current situation and countermeasures of the development of numerical intelligence in animal husbandry［J］.Animal Agriculture，2024（6）：59-63.（in Chinese）

（編輯 郭云雁）