馬士龍，李小偉，李響，謝書瓊，劉益麗，唐嬌，江明鋒*

（1.西南民族大學(xué)畜牧獸醫(yī)學(xué)院，四川 成都 610041；2.青藏高原動物遺傳育種資源保護(hù)與利用國家教育部重點實驗室，四川 成都 610041；3.西南民族大學(xué)青藏高原研究院，四川 成都 610041；4.四川省龍日種畜場，四川 紅原 624401）

牦牛是非常寶貴的家畜遺傳資源基因庫。中國擁有極其豐富的牦牛資源，總量有1600多萬頭，占世界的90%以上，全國形成了18個地方品種和2個培育品種［1］，保護(hù)好牦牛品種資源對于牦牛產(chǎn)業(yè)可持續(xù)發(fā)展具有重要意義。麥洼牦牛，屬于乳肉兼用型地方優(yōu)良品種，因中心產(chǎn)區(qū)是四川省阿壩藏族羌族自治州的紅原縣原麥洼部落而得名；據(jù)調(diào)查，麥洼牦牛是由阿壩牦牛引入康北牦牛、果洛牦牛的血緣，經(jīng)長期選育形成［2］；1987年被《四川省家畜家禽品種志》［3］收錄，1988年載入《中國牛品種志》［3］。由于早期阿壩州草場超載過牧、飼養(yǎng)管理不足、牛群近親繁殖等導(dǎo)致麥洼牦牛群體有退化趨勢，國家于1953年建立了四川龍日種畜場，成立麥洼牦牛保種場［4］，并經(jīng)長期的生產(chǎn)實踐，在2014年組建了純黑（全身被毛為黑色）、粉嘴（全身被毛為黑色，嘴唇毛色為粉色）和弗洛（全身被毛出生時黑色，逐漸變灰色）保種群［5］；其中，弗洛群的產(chǎn)肉性能高于粉嘴群和全黑群，純黑群的產(chǎn)奶性能高于粉嘴群和弗洛群。近年來，基于性能特征劃分的麥洼牦牛群體的生產(chǎn)性能均得到了提高［6］，但3個保種群的遺傳多樣性現(xiàn)狀，是否已經(jīng)出現(xiàn)遺傳分化的問題一直未展開相關(guān)研究；因此，對麥洼牦牛進(jìn)行遺傳多樣性和遺傳進(jìn)化分析，評價不同保種群的保種效果顯得尤為重要。

自20世紀(jì)以來，國內(nèi)學(xué)者已經(jīng)通過微衛(wèi)星標(biāo)記［7-8］、血液蛋白基因座標(biāo)記［9］、線粒體DNA標(biāo)記［10-12］等方法對麥洼牦牛、九龍牦牛、斯布牦牛等品種進(jìn)行群體遺傳進(jìn)化分析，然而這些DNA標(biāo)記僅能覆蓋部分且少量的牦牛基因組，所取得的基因組遺傳變異數(shù)據(jù)也是非常有限的。而新一代測序技術(shù)實現(xiàn)了單核苷酸多態(tài)性（single nucleotide polymorphism，SNP）分子標(biāo)記的高通量檢測，可從整個基因組水平上準(zhǔn)確高效地研究牦牛的遺傳進(jìn)化關(guān)系［13］?；蚍中蜏y序（genotyping by sequencing，GBS）技術(shù)是在二代測序基礎(chǔ)上發(fā)展的簡化基因組測序技術(shù)，通過酶切加tag的方法進(jìn)行測序，獲取酶切位點附近基因序列信息，進(jìn)而檢測大量且準(zhǔn)確性高的變異SNP位點信息［14］，被廣泛應(yīng)用于遺傳圖譜構(gòu)建、全基因組關(guān)聯(lián)分析和群體遺傳分析等領(lǐng)域［15-18］。

本研究利用GBS技術(shù)研究了麥洼牦牛3個保種群的遺傳多樣性、群體結(jié)構(gòu)和親緣關(guān)系，評價保種效果，為今后麥洼牦牛遺傳資源的保護(hù)與利用提供了可借鑒方法和理論依據(jù)，為下一步進(jìn)行麥洼牦牛的良種選育和生產(chǎn)性能提高奠定基礎(chǔ)。

1 材料與方法

1.1 試驗動物和方法

試驗于2020年6月在四川省阿壩藏族羌族自治州紅原縣國家級肉牛核心育種場—龍日種畜場進(jìn)行，該地區(qū)處于查真梁子山北部區(qū)域，海拔3500～4000 m，地勢平坦；只有冷、暖季之分，屬大陸性高原寒溫帶季風(fēng)氣候，牧草豐茂，是全國著名的地方優(yōu)良品種-麥洼牦牛、藏綿羊、河曲馬主產(chǎn)區(qū)［4］。從3個麥洼牦牛保種群全黑群（QH）、粉嘴群（FZ）、弗洛群（FL）中隨機選取406頭麥洼牦牛（公∶母=1∶3.5）（表1），群體內(nèi)個體間血緣關(guān)系未知，采集血樣樣本。具體步驟為用一次性注射器從牦牛頸外靜脈采集5 mL血液于乙二胺四乙酸（ethylene diamine tetraacetic acid，EDTA）抗凝管，放入冰盒中帶回，放置-80℃冰箱中保存?zhèn)溆谩?/p>

表1 3個保種群的采集信息Table 1 The sample information of 3 breeding populations

1.1.1 主要試劑及儀器 血液基因組DNA提取試劑盒（D3392-01）購自O(shè)mega Bio-Tek公司；異丙醇、無水乙醇、瓊脂糖、50×TAE、6×Loading Buffer、DNA Marker（BM5000）等購自上海生工生物工程技術(shù)服務(wù)有限公司和愛必信（上海）生物科技有限公司。

主要儀器有離心機（Eppendorf，5804R，德國）、電子天平（Denver，TP-214，美國）、電泳儀（Bio-rad，美國）、凝膠成像系統(tǒng)（Invitrogen，美國）、紫外分光光度計（Thermo Scientific，NanoDrop 2000，美國）；Qubit?2.0熒光計（Invitrogen，美國）。

1.1.2 DNA提取 按照血液基因組DNA提取試劑盒的說明書提取406個血液樣本的基因組DNA，利用1%瓊脂糖凝膠電泳（100 V，40 min）和紫外分光光度計檢測DNA的純度（A260nm/A280nm值）和完整性，使用Qubit?2.0進(jìn)行DNA濃度測定，樣品純度要求A260nm/A280nm值為1.8～2.0。

1.2 GBS試驗流程

1.2.1 GBS文庫構(gòu)建 文庫構(gòu)建參照董世武等［18］的方法，檢測合格的基因組DNA用Mse I限制性核酸內(nèi)切酶進(jìn)行酶切，酶切后基因組片段兩端加上帶有barcode的接頭（120 bp）后，利用PCR對每個樣品進(jìn)行擴增，然后對樣品進(jìn)行混合，選擇需要的片段（265～315 bp）進(jìn)行文庫構(gòu)建。

1.2.2 文庫庫檢及上機測序 構(gòu)建好的文庫用Qubit?2.0進(jìn)行初步定量，具體操作為稀釋文庫至1 ng·μL-1，再用Agilent 2100檢測文庫的插入片段大小是否符合265～315 bp，然后為保證文庫質(zhì)量，使用Q-PCR方法對文庫的有效濃度進(jìn)行準(zhǔn)確定量（文庫有效濃度＞2 nmol·L-1）。庫檢合格后，把不同文庫按照有效濃度及每個樣本可產(chǎn)生10萬個標(biāo)簽的目標(biāo)下機數(shù)據(jù)量的需求進(jìn)行樣本混合，即混池，然后進(jìn)行Illumina Hiseq PE150雙末端測序。

1.3 數(shù)據(jù)統(tǒng)計與分析

1.3.1 測序數(shù)據(jù)及SNP的質(zhì)控 Illumina Hiseq測序平臺得到的原始圖像文件經(jīng)堿基識別分析轉(zhuǎn)化為測序序列數(shù)據(jù)（raw reads），這些raw reads需進(jìn)行質(zhì)量控制，剔除含有接頭序列的reads和低質(zhì)量reads之后得到的數(shù)據(jù)即為有效數(shù)據(jù)（clean data）。

以野牦牛（Bos mutus）（生物項目數(shù)據(jù)庫：PRJNA221623）為參考基因組，通過BWA軟件（參數(shù)：mem-t 4-k 32-M）［19］將雙端有效數(shù)據(jù)比對到參考基因組，統(tǒng)計平均比對率、平均測序深度和平均覆蓋度。采用SAMTOOLS［20］軟件進(jìn)行群體SNP的檢測，得到位點變異文件，以測序深度（sequencing depth）為3，缺失（miss）為0.2，次要等位基因頻率（minor allele frequency，MAF）大于0.05條件過濾SNP位點，最后獲得高質(zhì)量的SNP位點用于后續(xù)分析研究。

1.3.2 群體結(jié)構(gòu)與遺傳多樣性分析 基于個體間SNP差異，通過GCTA［21］軟件計算特征向量和特征值，用R語言繪制第一主成分與第二主成分的主成分分析（principal component analysis，PCA）二維圖。使用TreeBest［22］軟件中的P-distances模型計算遺傳距離矩陣并構(gòu)建系統(tǒng)進(jìn)化樹，引導(dǎo)值（bootstrap values）設(shè)置為1000，檢驗可靠性。采用ADMIXTURE軟件構(gòu)建群體遺傳結(jié)構(gòu)和群體世系信息，假定的祖先群體個數(shù)為2～20，根據(jù)交叉驗證錯誤值（cross validation error，CV）確定最優(yōu)K值（一般CV最小的為最優(yōu)K值）。使用PopLDdecay（v 3.40）［23］進(jìn)行連鎖不平衡分析（linkage disequilibrium，LD），并用R軟件繪制LD衰減圖。

利用PLINK v 1.90［24］計算觀測雜合度（observed heterozygosity，Ho）和期望雜合度（expected heterozygosity，He）等遺傳多樣性指標(biāo)。使用R包genepop統(tǒng)計分析每個SNP位點在所有個體中的近交系數(shù)（inbreeding coefficients，F(xiàn)is）和群體遺傳分化系數(shù)（genetic differentiation index，F(xiàn)st）?；诤Y選后高質(zhì)量SNP信息，進(jìn)行選擇消除分析，采用PLINK軟件進(jìn)行選擇信號檢測，以100 Kb為窗口、10 Kb為step進(jìn)行滑動的核苷酸多樣度（π）比率和Fst計算，分別以π比率和Fst值top 10%作為閾值，對群體選擇信號進(jìn)行分析，其top 10%的交集為候選位點，并對受選擇位點所在基因進(jìn)行GO富集分析。

2 結(jié)果與分析

2.1 血液基因組DNA提取

提取406頭牦牛樣本血液基因組DNA，經(jīng)1%瓊脂糖凝膠電泳檢測發(fā)現(xiàn)，條帶清晰明亮（圖1），各項參數(shù)均符合測序所需標(biāo)準(zhǔn)。

圖1 部分樣本血液基因組DNA電泳圖Fig.1 Genome DNA electrophoresis of some blood samples

2.2 測序數(shù)據(jù)產(chǎn)出與質(zhì)控

對406份牦牛血液基因組進(jìn)行GBS測序和質(zhì)量控制（詳見1.3.1），獲得有效堿基數(shù)279.211～2456.350 Mb，共產(chǎn)生327.517 Gb原始數(shù)據(jù)，平均為806.690 Mb；過濾后共產(chǎn)生327.512 Gb有效數(shù)據(jù)，有效率達(dá)99.990%，平均為806.680 Mb。堿基檢測錯誤率為0.040%～0.060%，平均僅有0.044%（表2）；達(dá)到Q20質(zhì)量的堿基占比為93.83%～97.58%，平均為95.73%；GC堿基含量為36.56%～40.28%，平均為38.88%；獲得有效讀長數(shù)為1.94～17.08 Mb，平均為5.60 Mb；有效讀長比對到參考基因上的占比為92.60%～99.68%，平均比對率為98.75%；測序堿基總量與基因組的比值，即測序深度最小為5.65×，測序深度最大為31.64×，平均測序深度為14.33×，平均覆蓋度為3.46%（表3）。表明測序質(zhì)量高，GC分布正常，樣本均沒有被污染，此次建庫測序成功。

表2 麥洼牦牛測序數(shù)據(jù)統(tǒng)計Table 2 Summary of sequencing data of Maiwa yak

表3 測序數(shù)據(jù)質(zhì)量評估Table 3 Quality assessment of sequencing

2.3 SNP檢測

應(yīng)用SAMTOOLS等軟件對測序有效數(shù)據(jù)進(jìn)行SNP檢測，共鑒定出1628805個SNP位點，進(jìn)行質(zhì)控后獲得高質(zhì)量SNP位點126122個。質(zhì)控后的SNP經(jīng)過與牦牛參考基因組（生物項目數(shù)據(jù)庫：PRJNA 221623）比對分析，基因組的SNPs數(shù)目及基因組序列長度如圖2所示，總體上來看SNPs的數(shù)目與基因組序列長度總體上呈正相關(guān)；序列號為NW_005393843.1的片段最大（8785736 bp），所含SNP數(shù)目最多，為463個SNP位點；而由于基因組序列長度差異，不同基因組序列的SNP數(shù)目變化較大，NW_005393917.1的 片 段 較 大（542197 bp），卻 僅 有1個SNP位點。

圖2 質(zhì)控SNPs位點數(shù)與基因組序列片段長度變化趨勢Fig.2 Variation trend of SNPs number and genome sequence length

2.4 遺傳多樣性和群體結(jié)構(gòu)分析

2.4.1 遺傳多樣性分析 根據(jù)變異位點分析QH、FZ、FL群的雜合度（Ho、He）、核苷酸多樣度（π）、近交系數(shù)（Fis），結(jié)果見表4。3個群的觀測雜合度分別為0.3029、0.3042、0.3044，觀測雜合度均與期望雜合度接近，只有FZ的觀測雜合度稍高于期望雜合度，總體平均觀測雜合度與期望雜合度接近。3個群的核苷酸多樣度分別為0.3030、0.3020、0.3009，F(xiàn)L群最低；而與QH、FZ群的近交系數(shù)相比，F(xiàn)L群是最高的，為0.0209。

表4 麥洼牦牛3個育種群遺傳多樣性參數(shù)Table 4 Genetic diversity parameters of three preserved populations of Maiwa yak

連鎖不平衡（LD）模式見圖3，3個群體LD均隨距離增長逐漸平滑下降，各等位基因不存在強烈的連鎖不平衡水平（r2≤0.5），SNP距離 在0～100 kb的LD衰減 速度較快，LD系數(shù)r2約從0.5降至0.2以 下，SNP距離在100～300 kb的LD水平較低。當(dāng)r2=0.3時，各群體對應(yīng)LD衰減距離基本相等，說明各群體LD衰減速度差異不大。

圖3 LD隨距離增長的衰減Fig.3 Attenuation of LD with distance

由表5可以看出，3個群體遺傳分化系數(shù)Fst為0.02504～0.03513（介于0～0.05），其中QH和FZ群之間的遺傳分化程度最大（0.03513），兩群的遺傳距離同樣是最遠(yuǎn)（0.0358）。

表5 群體間的遺傳分化系數(shù)（Fst）（對角線下）和遺傳距離（DR）（對角線上）Table 5 Genetic differentiation index（below diagonal）and genetic distance（above diagonal）among three populations of Maiwa yak

2.4.2 群體遺傳結(jié)構(gòu)和親緣關(guān)系分析 基于個體基因組SNP差異程度，通過PCA二維圖的形式來展示個體的聚類情況，可以大致了解群體的遺傳結(jié)構(gòu)。由圖4可知，麥洼牦牛群體的PC1、PC2的貢獻(xiàn)率分別為4.75%、2.72%，F(xiàn)L群緊密聚在一起，與QH、FZ群部分個體關(guān)系近；FZ群個體較分散，QH群中有個體偏離群體與FZ群聚在一起，兩群存在明顯的分層。

圖4 PCA分析二維圖Fig.4 PCA analysis of two-dimensional graph

群體遺傳結(jié)構(gòu)指遺傳變異在物種或群體中的一種非隨機分布。按照地理分布或親緣關(guān)系可將一個群體分為若干亞群，處于同一亞群內(nèi)的不同個體親緣關(guān)系較近，而亞群與亞群之間則親緣關(guān)系稍遠(yuǎn)。從圖5可以看出，3個保種群的血統(tǒng)構(gòu)成，當(dāng)假設(shè)祖先群為2（K=2）時，QH、FZ、FL群都攜帶兩個祖代血緣；當(dāng)最 佳K值 為12時 ，F(xiàn)L群 可 分 為 兩 個群 ，F(xiàn)Z群 部 分 個體（棕色）可單獨劃分為一個群，QH群部分個體（粉色）可單獨劃分一個群，其他需根據(jù)系譜記錄進(jìn)行劃分。

圖5 麥洼牦牛structure群體遺傳結(jié)構(gòu)Fig.5 Genetic structure of Maiwa yak structure populations

406頭麥洼牦牛聚類的系統(tǒng)進(jìn)化樹（圖6）直觀地以圈圖的形式展示它們之間的親緣關(guān)系遠(yuǎn)近，結(jié)果可以看出，QH群聚成一個大支（綠色），群體存在離群個體；FL和FZ群組成另一支，F(xiàn)Z群逐漸分離出來，最后形成一小支（藍(lán)色），F(xiàn)L群有一支最后被分離出來（紅色），但整體有一半的個體與QH、FZ群聚在一起。

圖6 系統(tǒng)進(jìn)化樹Fig.6 Phylogenetic tree

2.4.3 群體選擇信號分析 Fst和π分析是兩種有效地檢測選擇消除區(qū)域的方法，特別是在挖掘與生存環(huán)境密切相關(guān)的基因功能區(qū)時可以得到較強的選擇信號，便于目標(biāo)基因的篩選。根據(jù)麥洼牦牛群體間Fst和π進(jìn)行選擇信號分析，篩選出3513個差異較大的SNP位點，共有104個差異基因受到選擇，進(jìn)行GO功能注釋分析，結(jié)果表明（圖7），這些受選擇基因主要富集在有絲分裂、氯離子運輸、節(jié)律過程、刺激反應(yīng)等生物學(xué)過程；主要涉及的細(xì)胞組成：氯離子通道復(fù)合體、細(xì)胞間連接、膜的錨定組分等；主要涉及的分子功能包括：氯離子通道活性、酶結(jié)合、蛋白質(zhì)二聚化活性等。KEGG結(jié)果富集到了28條通路，主要集中在生殖激素信號、內(nèi)/外分泌、鈣離子信號、cGMP-PKG信號、血管平滑肌收縮等調(diào)控通路上（表6）。

表6 受選擇基因的KEGG通路Table 6 KEGG pathway of differential expression genes

圖7 受選擇基因的GO富集條目Fig.7 GO enrichment items of selected gene

3 討論

牦牛是我國青藏高原珍稀而寶貴的遺傳資源，在高原畜牧業(yè)發(fā)展和生態(tài)系統(tǒng)平衡方面起著至關(guān)重要的作用。牦牛遺傳多樣性水平和遺傳結(jié)構(gòu)的研究不僅是保護(hù)牦牛遺傳資源的主要內(nèi)容，也是實施科學(xué)有效的品種保護(hù)措施的前提條件。隨著分子生物技術(shù)的快速發(fā)展，基于SNP的新一代測序技術(shù)成為主流方法，使得從基因組層面上研究牦牛遺傳結(jié)構(gòu)和遺傳多樣性成為可能，更全面精確地監(jiān)測牦牛保種效果。自麥洼牦牛品種被挖掘以來，其遺傳多樣性研究是保護(hù)這個地方牦牛品種資源的核心問題之一；因此，本研究對406頭麥洼牦牛進(jìn)行GBS測序探究其群體的遺傳多樣性及遺傳結(jié)構(gòu)，測序結(jié)果得到了806.68 Mb的有效數(shù)據(jù)量，共鑒定出126122個SNP位點，這些數(shù)據(jù)將為未來制備麥洼牦牛育種SNP芯片提供基礎(chǔ)。

遺傳雜合度是度量群體遺傳多樣性的最適參數(shù)，其觀測雜合度和期望雜合度的接近程度可反映群體遺傳多樣性的高低。本研究中麥洼牦牛3個保種群的觀測雜合度為0.3029～0.3044，屬于中度雜合群體，整個群的平均觀測雜合度與平均期望雜合度接近，這表明麥洼牦牛品種具有較高純度和豐富的遺傳多樣性，這與趙芳芳［25］的微衛(wèi)星法和師方等［11］的隨機擴增多態(tài)性脫氧核糖核酸（random amplified polymorphic DNA，RAPD）法結(jié)果相似；龍日種畜場為高原淺丘地貌，可利用草地面積1.667萬hm2，在這樣的相對隔離地理環(huán)境中麥洼牦牛得到了有效保護(hù)，較低的近交水平也說明保種群并沒有受外來品種及近交的影響。另外，粉嘴群的觀測雜合度高于全黑群，結(jié)果與侯孟典等［26］和柴志欣等［27］的報道一致；這表明粉嘴群所經(jīng)過的長期定向人工選擇強度高于全黑群，有研究［5］指出粉嘴個體相較于全黑個體的乳脂率較高，而擠奶量較低，推測粉嘴群已經(jīng)在乳質(zhì)等生產(chǎn)性狀上受到了強烈的選擇。群體間的遺傳分化與選育的方式和世代長短有關(guān)，四川省龍日種畜場的牦牛原種場采取的開放式群體繼代選育法，在世代選育過程中，按需地從麥洼牦牛主產(chǎn)區(qū)的選育群引進(jìn)優(yōu)良非近親種公牦牛進(jìn)行血緣更新［3］，使得場內(nèi)的麥洼牦牛群體內(nèi)個體間遺傳變異度變高；因此粉嘴和全黑群的遺傳分化指數(shù)和遺傳距離較遠(yuǎn)，分化呈現(xiàn)出不穩(wěn)定的混亂狀態(tài)，有部分血緣交叉和個體離群，這提示在保種過程中，應(yīng)該加強對種公牛的嚴(yán)格選留工作，避免近親交配的發(fā)生。另外，早期的弗洛個體其實并沒有從總?cè)褐袉为殑澐殖鰜?，而是后來依?jù)其毛色特征和產(chǎn)肉性能才組建的［5］，這印證了弗洛群是粉嘴群和全黑群的祖代個體雜交所形成的群落的可能；進(jìn)一步說明弗洛群含有部分全黑群和粉嘴群血統(tǒng)，有助于解釋弗洛群的出生毛色變化現(xiàn)象。龍日種畜場內(nèi)麥洼牦牛群體的遺傳結(jié)構(gòu)隨著時間的推移也發(fā)生了較大變化，推測是由較小的保種群體規(guī)模限制導(dǎo)致，群體結(jié)構(gòu)圖和系統(tǒng)進(jìn)化樹都表明全黑群大部分個體聚為一類，弗洛群小部分個體聚為一類，這說明全黑、弗洛群的保種效果相對較好，全黑群和弗洛群出現(xiàn)血統(tǒng)較純正的個體群，這些個體和粉嘴群一樣，它們的某些經(jīng)濟性狀可能已經(jīng)得到了固定。

為鑒別3個保種群之間在基因組特定位點上存在的分化，對麥洼牦牛群體104個受到強烈選擇的基因進(jìn)行GO和KEGG分析，發(fā)現(xiàn)這些基因廣泛參與生殖機能調(diào)控、生長發(fā)育、刺激反應(yīng)等生物學(xué)相關(guān)過程。其中，特異型鈣調(diào)磷酸酶PPP3CC基因參與精子發(fā)生和睪丸成熟，缺失該基因的小鼠精子尾部中段會發(fā)生僵直、超活化能力降低［28］，該基因也已經(jīng)在牛體內(nèi)被檢測到［29］；鉀離子通道宿主因子KCNMA1與牛的繁殖性能相關(guān)，可增強動物性行為［30］。神經(jīng)生長因子NGF能有效促進(jìn)卵泡形成，通過調(diào)控原始卵泡中扁平細(xì)胞分化成次級卵泡中包裹卵母細(xì)胞的立方形顆粒細(xì)胞過程，在牛孤雌胚胎的早期發(fā)育起關(guān)鍵作用［31］；鳥嘌呤核苷酸結(jié)合蛋白GNAQ基因是通過影響小腦分泌激素從而影響動物毛色形成的重要調(diào)控基因［32］，也參與動物性腺的發(fā)育和生殖細(xì)胞的成熟［33］。這表明3個保種群在龍日種畜場的人工選擇條件下，公母牛的繁殖力均得到了增強，尤其公牛的性行為方面較為明顯。肌細(xì)胞增強因子MEF2C主要參與調(diào)控肌細(xì)胞分化，特別是在從胎牛到成年牛發(fā)育過程中，介導(dǎo)骨骼肌、心肌和平滑肌的細(xì)胞分化［34］；小GTP結(jié)合蛋白Rho下游效應(yīng)因子ROCK2對低氧環(huán)境誘導(dǎo)的肺動脈高血壓發(fā)揮著重要作用，并在成肌細(xì)胞分化中持續(xù)上調(diào)表達(dá)，參與肌肉形成過程［35］；另外，KCNJ3基因與心律發(fā)生有關(guān)［36］，這3個基因可能與麥洼牦牛的高原適應(yīng)性和肌肉發(fā)育調(diào)節(jié)相關(guān)，說明保種群在強烈的人工選擇中有向肉質(zhì)性狀方面發(fā)生了定向選擇的趨勢。谷氨酸突觸通路中的ITPR2基因是鈣離子跨膜轉(zhuǎn)運活性的關(guān)鍵調(diào)控因子，內(nèi)含子區(qū)存在動物抗環(huán)境壓力的潛在靶點［37］，這表明保種群在龍日種畜場的天然草場放牧和一定的補飼管理下，群體的抗應(yīng)激能力也得到了強化。酪氨酸蛋白激酶受體KIT基因是調(diào)控哺乳動物白色皮膚及被毛的基因，KIT基因突變會導(dǎo)致動物由花斑到全白不同程度的白色被毛，而牦牛的白色斑點就是由KIT基因在6號染色體上的易位形成的［38］，提示該基因可作為弗洛群的灰色表型發(fā)生的候選基因；另外，仍有許多受選擇功能基因有待進(jìn)一步研究注釋。本研究在基因組水平上進(jìn)一步說明龍日種畜場內(nèi)麥洼牦牛不同保種群經(jīng)歷了定向選擇，篩選的受選擇功能基因?qū)α私恹溚蓐笈５倪m應(yīng)性變化并獲得其優(yōu)良的經(jīng)濟性狀具有非常深遠(yuǎn)的意義。

4 結(jié)論

本研究利用GBS技術(shù)揭示了麥洼牦牛3個保種群的遺傳多樣性和遺傳結(jié)構(gòu)，粉嘴群和全黑群有遺傳分化趨勢，分群的效果得到了體現(xiàn)，3個保種群的保種效果較好，可在世代之間進(jìn)行比較，進(jìn)一步分析保種效果。測序出3513個顯著SNP位點和104個差異基因；3個保種群在毛色、肉質(zhì)、繁殖性狀上經(jīng)歷了人工選擇，這可為麥洼牦牛的保種選育工作提供建議。