摘 要： 旨在分析子午嶺黑山羊的群體遺傳多樣性和親緣關(guān)系以及家系結(jié)構(gòu)，為子午嶺黑山羊的保護和利用提供依據(jù)。本研究通過簡化基因組測序（super-genotyping by sequencing，Super-GBS）技術(shù)對99只（10公/89母）成年子午嶺黑山羊進行全基因組SNP檢測。利用Plink軟件計算觀測雜合度（Ho）、期望雜合度（He）、多態(tài)信息含量（PIC）、核苷酸多樣性（Pi）、有效等位基因數(shù)（Ne）及次要等位基因頻率（MAF）等6項遺傳多樣性指標(biāo)；GCTA軟件進行主成分分析及基因組親緣關(guān)系G矩陣構(gòu)建；Plink軟件構(gòu)建IBS遺傳距離矩陣，R語言繪制熱圖；PHYLP構(gòu)建系統(tǒng)發(fā)育樹；detect RUNS工具檢測ROH。結(jié)果表明，99只子午嶺黑山羊個體共檢測到996 042個SNPs。子午嶺黑山羊群體的PIC、Pi、Ne及MAF值分別是0.161、0.193、1.295、0.130，且Ho（0.167）低于He（0.192），說明子午嶺黑山羊群體遺傳多樣性較低。G矩陣和IBS遺傳距離結(jié)果均表明子午嶺黑山羊群體間大部分個體間親緣關(guān)系較遠。主成分分析結(jié)果揭示子午嶺黑山羊群體內(nèi)部不存在明顯分化，系統(tǒng)發(fā)育樹結(jié)果說明子午嶺黑山羊群體公羊可大致分為6個家系，所有家系公羊數(shù)量較少。子午嶺黑山羊群體的近交系數(shù)FROH為0.049 6，說明子午嶺黑山羊群體內(nèi)部近交程度相對較低。綜上所述，子午嶺黑山羊群體遺傳多樣性較低，大部分個體間親緣關(guān)系較遠，群體內(nèi)近交程度較低，后期應(yīng)注意后代的選育，避免血統(tǒng)流失。

關(guān)鍵詞： 子午嶺黑山羊；簡化基因組測序；群體結(jié)構(gòu)；遺傳多樣性；家系結(jié)構(gòu)；親緣關(guān)系

中圖分類號：S827.2

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

文章編號：0366-6964（2024）10-4334-12

收稿日期：2024-03-25

基金項目：甘肅省自然基金重點項目（22JR5RA759）；甘肅省重點人才項目（2023RCXM67）

作者簡介：茍想珍（1978-），女，甘肅通渭人，碩士生，正高級畜牧師，主要從事羊遺傳育種與繁殖的研究，E-mail： gxz1621@sina.com

*通信作者：王 珂，主要從事動物遺傳育種與繁殖，E-mail：592374483@qq.com；馬友記，主要從事羊生產(chǎn)研究，E-mail：yjma@gsau.edu.cn

Evaluation of the Population Structure of the Ziwuling Black Goat Based on Super-GBS

GOU" Xiangzhen1， YANG" Junxiang1， ZHAO" Zihui1， FENG" Lingxia1， CHEN" Wanhui1， LI" Yujie2，

ZHANG" Zhongyu2， MA" Keyan3， JIANG" Dongping4， CHANG" Rong5， WEN" Yazhou6， WANG" Ke1*， MA

Youji3*

（1.Key Laboratory of Livestock and Poultry Resources （Goat） Evaluation and Utilization of Ministry

of Agriculture and Rural Affairs， Animal Husbandry and Veterinary Research Institute of Gansu， Pingliang

744000，" China;

2.Station of Huachi County Animal Science and Veterinary， Qingyang 745000，

China;

3.College of Animal Science and Technology， Gansu Agricultural University， Lanzhou 730070，

China;

4.Pingchuan District Animal Husbandry and Veterinary Technical Service Center， Baiyin 730900，

China;

5.Animal Health Supervision Institute of Ganzhou District， Zhangye 734000，" China;

6.

Gansu Sheep Breeding Technology Extension Station， Zhangye 734000，" China）

Abstract：" This study aimed to assess the genetic diversity and population structure of the Ziwuling black goat， laying the groundwork for the preservation and utilization of this breed. The whole-genome SNP analysis of a total of 99 adult Ziwuling black goats （10 males and 89 females） was performed by super-genotyping by sequencing（super-GBS）. Genetic diversity parameters including observed heterozygosity （Ho）， expected heterozygosity （He）， polymorphic information content （PIC）， nucleotide diversity （Pi）， effective allele number （Ne）， and minor allele frequency （MAF） were calculated using Plink software. Principal component analysis and construction of genetic relationship G matrix were conducted using GCTA software. IBS genetic distance matrix was built using Plink software， and a heatmap was plotted using R language. A phylogenetic tree was constructed using PHYLP， while runs of homozygosity （ROH） were detected using the detect RUNS tool. A total of 996 042 SNPs were detected in the 99 individuals of Ziwuling black goats. The PIC， Pi， Ne， and MAF values of Ziwuling black goats were 0.161， 0.193， 1.295， 0.130， respectively， indicating a relatively low genetic diversity with Ho（0.167） lower than He（0.192）. Both the G matrix and IBS genetic distance results showed that most individuals in the Ziwuling black goat population had distant genetic relationships. The principal component analysis revealed no obvious differentiation within the Ziwuling black goat population， while the phylogenetic tree indicated that the rams could be roughly divided into 6 lineages， with few rams in each lineage. The inbreeding coefficient FROH of the Ziwuling black goat population was 0.049 6， suggesting a relatively low level of inbreeding within the population. In conclusion， the genetic diversity of the Ziwuling black goat population is relatively low， with most individuals having distant kinship relations and a low level of inbreeding within the group. It is recommended to pay attention to the breeding of offspring in order to avoid loss of genetic diversity in the future.

Key words： Ziwuling black goat; Super-GBS；population structure; genetic diversity; family structure；kinship

*Corresponding authors： WANG Ke， E-mail：592374483@qq.com; MA Youji， E-mail：yjma@gsau.edu.cn

畜禽遺傳資源的多樣性是培育優(yōu)良畜禽品種的物質(zhì)基礎(chǔ)，可以提供應(yīng)對未知需求的適應(yīng)性［1-2］。因此，保護畜禽遺傳資源對于提升畜牧業(yè)綜合實力、可持續(xù)發(fā)展能力、種源自給率以及市場競爭力具有重要意義。

子午嶺黑山羊于1982年被列入《甘肅省畜禽品種志》，命名為隴東黑山羊。2011年收錄于《中國畜禽遺傳資源志 羊志》。在2021年公布的《國家畜禽遺傳資源名錄（2021版）》中被正式命名為子午嶺黑山羊［3］。子午嶺黑山羊原產(chǎn)于子午嶺兩側(cè)的甘肅省華池縣、環(huán)縣、合水縣和陜西省陜北及渭北等地，現(xiàn)主要集中于甘肅省東部（甘肅省寧縣、鎮(zhèn)原、慶城、定西等地區(qū)），截至2021年，現(xiàn)存欄子午嶺黑山羊126 143只。然而，由于畜牧業(yè)結(jié)構(gòu)和市場結(jié)構(gòu)調(diào)整等因素，為提高產(chǎn)絨量，當(dāng)?shù)匾脒|寧絨山羊及內(nèi)蒙古絨山羊與子午嶺黑山羊進行雜交，導(dǎo)致純種子午嶺黑山羊的數(shù)量急劇下降，面臨滅絕的威脅。子午嶺黑山羊現(xiàn)有一個保種場——隴東黑山羊保種育種場，其中現(xiàn)存欄子午嶺黑山羊1 000余只，出欄500只左右。因此，迫切需要開展保護工作。目前的文獻報道中，關(guān)于子午嶺黑山羊的遺傳多樣性和遺傳結(jié)構(gòu)的研究尚屬空白，這限制了對其資源的有效保護和利用。

單核苷酸多態(tài)性（single nucleotide polymorphism， SNP）是廣泛應(yīng)用于基因組水平的遺傳標(biāo)記基因分型方法，通過對群體內(nèi)的個體進行全基因組范圍內(nèi)的SNPs比較分析，可以揭示群體的遺傳多樣性、個體間的親緣關(guān)系以及群體結(jié)構(gòu)［4-7］?；蚍中蜏y序（genotyping by sequencing，GBS）技術(shù)是在二代測序基礎(chǔ)上發(fā)展的簡化基因組測序技術(shù)［8］，能夠檢測大量且準(zhǔn)確性高的變異SNP位點信息，廣泛應(yīng)用于高分辨率的遺傳圖譜構(gòu)建、群體遺傳學(xué)、物種鑒定、基因組演化等研究領(lǐng)域［9-11］。例如，馬士龍等［12］基于GBS技術(shù)研究了3個麥洼牦牛群體的遺傳結(jié)構(gòu)，為其后續(xù)保種計劃的制定提供了依據(jù)。劉海林等［13］利用GBS技術(shù)對60只湘東黑山羊進行測序，確立了湘東黑山羊的家系結(jié)構(gòu)分析，為湘東黑山羊的保種和選育提供了參考依據(jù)。Super-GBS技術(shù)作為GBS技術(shù)的進階，在SNP標(biāo)記的重現(xiàn)性、穩(wěn)定性、分型準(zhǔn)確性等特點上優(yōu)于GBS，提供了更高的位點密度和更多的SNPs，已在遺傳圖譜構(gòu)建、QTL定位和全基因組關(guān)聯(lián)分析等領(lǐng)域展現(xiàn)出其應(yīng)用潛力［14］。

因此，本研究采用Super-GBS技術(shù)，目的是深入探究子午嶺黑山羊保種群體的遺傳多樣性、群體結(jié)構(gòu)和親緣關(guān)系，旨在為該品種的保護和有效利用提供科學(xué)依據(jù)。

1 材料與方法

1.1 試驗動物及樣品采集

本研究選擇來自隴東黑山羊保種育種場的99只（10公/89母）無親緣關(guān)系子午嶺黑山羊為研究對象，所有個體的飼養(yǎng)標(biāo)準(zhǔn)與環(huán)境條件一致。其中HC001、HC016、HC036、HC068、HC071、HC072、HC076、HC077、HC087和HC093為公羊，其余均為母羊。

頸靜脈采血5 mL，用含抗凝劑的采血管收集，顛倒混勻8～10次，油性記號筆標(biāo)記樣品信息并編號，24 h內(nèi)帶回實驗室-20 ℃～-80 ℃冷凍后，干冰運輸。

1.2 試驗方法

1.2.1 基因組DNA提取與質(zhì)量鑒定

按照血液基因組DNA提取試劑盒的說明書提取99個血液樣本的基因組DNA，利用NanoDrop 2000檢測DNA的純度，用1%瓊脂糖凝膠電泳檢測DNA的完整性，使用Qubit 測定每個樣品的最終濃度。

1.2.2 簡化基因組測序

所有樣品委托上海歐易生物科技有限公司進行建庫測序。選擇GCF_001704415.2作為參考基因組進行酶切預(yù)測，最終使用PstI-HF/MspI對DNA進行酶切，長度為300～700 bp 的片段被回收進行PCR擴增后，使用Qubit測定PCR產(chǎn)物濃度，文庫質(zhì)檢合格后使用Illumina Novaseq PE150測序平臺進行測序。

1.3 數(shù)據(jù)處理

1.3.1 數(shù)據(jù)質(zhì)控與基因組比對及SNP變異鑒定

對Raw Reads進行質(zhì)控后得到 Clean Reads（質(zhì)控標(biāo)準(zhǔn)：去除接頭序列；去除 N（非AGCT）堿基大于等于5的Reads；去除 5＇端和 3＇端低質(zhì)量堿基；去除長度小于75 bp 的 Reads；去除質(zhì)量低于15的堿基數(shù)超過整條Reads堿基數(shù)40%的低質(zhì)量Reads）。利用bwa［15］軟件將Clean Reads比對到參考基因組上并使用GATK［16］軟件進行SNP檢測。為了降低SNP檢測的錯誤率，選用QD≥2.0的標(biāo)準(zhǔn)進行過濾，只保留同時滿足該條件的突變位點。使用SnpEff［17］軟件對得到的SNP進行注釋，以確定SNP在基因元件的位置、對氨基酸的變化影響等。

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

利用Plink（V1.90）軟件計算觀測雜合度（Ho）、期望雜合度（He）、多態(tài)信息含量（PIC）、核苷酸多樣性（Pi）、有效等位基因數(shù)（Ne）及次要等位基因頻率（MAF）等遺傳多樣性指標(biāo)。

本研究采用鄰接法（neighbor-joining method）［18］構(gòu)建所有個體的進化樹。通過treebest［19］軟件計算距離矩陣，進化樹的可靠性通過bootstrap法進行檢驗（重復(fù)1 000 次）。鑒于公羊?qū)Ρ７N群的重要性，采用PHYLIP中的neighbor對10只公羊個體單獨進行了系統(tǒng)進化樹的構(gòu)建。使用GCTA（V1.92.1）軟件［20］對獲得的SNPs標(biāo)記進行PCA分析，得到影響最大的兩個特征向量。

1.3.3 親緣關(guān)系分析

利用Plink（V1.90）軟件計算群體個體間的狀態(tài)同源（identical by state， IBS）和個體間的遺傳距離。通過GCTA（V1.92.1）軟件分析群體親緣關(guān)系，通過R軟件繪制G矩陣和IBS距離矩陣熱圖。

1.3.4 ROH分析

通過ROH分析可以用于評估個體的近交程度，ROH越長，近交程度越高。本研究使用R包（detectRUNS）［21］工具，在window-free的情況下統(tǒng)計子午嶺黑山羊群體的ROH分布。通過Plink（V1.90）軟件計算得到每個樣本的ROH長度（參數(shù)：minSNP=20，ROHet=FALSE，maxGap=106，minLengthBps=250 000，maxOppRun=1，maxMissRun=1），進而分析得到子午嶺黑山羊的近交系數(shù)［22］。

2 結(jié) 果

2.1 DNA提取和質(zhì)量檢測

共提取99份血液DNA，質(zhì)檢合格，瓊脂糖凝膠電泳檢測結(jié)果條帶清晰明亮（圖1），可用于后續(xù)基因組測序。

2.2 數(shù)據(jù)產(chǎn)出、質(zhì)控及SNPs分布

本試驗共獲得1 795 029 254個raw reads，過濾后得到1 762 002 636個clean reads。由表1可知，樣本比對率范圍為99.66%～99.94%，平均測序深度為6.7×，測序平均Q30為90.98%，平均GC含量為51.18%，數(shù)據(jù)質(zhì)量達到后續(xù)分析要求。經(jīng)過突變位點檢測并過濾，共獲得996 042個SNPs位點，根據(jù)其在染色體上的位置，繪制染色體分布圖，發(fā)現(xiàn)SNPs非均勻分布在各染色體上，其中Chr1 SNPs數(shù)量最多，為51 223個，Chr27 SNPs數(shù)量最少，為19 145個（圖2）。由圖3可知，大部分SNPs位于基因間區(qū)（55.79%）和內(nèi)含子區(qū)域（38.81%）。

2.3 遺傳多樣性分析

根據(jù)表2可知，子午嶺黑山羊群體SNP位點的多態(tài)信息含量（PIC）平均值為0.161。平均觀測雜合度（Ho）為0.167，低于期望雜合度（He）0.192。值得注意的是，子午嶺黑山羊群體內(nèi)部部分個體發(fā)生了純合突變，則該位點觀測雜合度為0。有效等位基因數(shù)（Ne）平均值為1.295。核苷酸多樣性（Pi）平均值為0.193，次要等位基因頻率（MAF）平均值為0.130。

2.4 遺傳距離和親緣關(guān)系分析

通過G矩陣分析子午嶺黑山羊個體之間的親緣關(guān)系，根據(jù)其可視化結(jié)果（圖4）可知，子午嶺黑山羊群體大部分個體間親緣關(guān)系較遠（藍色部分），少數(shù)個體親緣關(guān)系較近（偏紅部分）。

IBS遺傳距離矩陣（圖5）結(jié)果表明，子午嶺黑山羊群體IBS遺傳距離為0.817 0～0.9862，平均遺傳距離為0.831 3。其中10只公羊之間的平均遺傳距離為0.833 6。IBS遺傳距離矩陣與G矩陣結(jié)果一致，大部分個體親緣關(guān)系較遠（藍色小格），少部分個體間親緣關(guān)系較近（偏黃及偏紅的方格）。

2.5 群體結(jié)構(gòu)分析

2.5.1 主成分分析

由PCA（圖6）可知，第一主成分和第二主成分分別解釋了2%和1%的遺傳變異。主成分分析結(jié)果顯示，子午嶺黑山羊大部分個體緊密聚集在一起，無明顯分層，但仍然存在一定程度的遺傳異質(zhì)性，表現(xiàn)為部分個體的分散分布。特別值得注意的是，10只公羊樣本中HC001與其他9只公羊樣本在遺傳上并未形成明顯的聚集，這一發(fā)現(xiàn)表明子午嶺黑山羊的保種群體可能由至少兩個不同的家系組成。

2.5.2 系統(tǒng)發(fā)育樹分析

系統(tǒng)發(fā)育樹（圖7A）結(jié)果表明，子午嶺黑山羊群體分布在兩個主要分支上。鑒于公羊?qū)φ麄€子午嶺黑山羊保種群體的重要性，本研究對10只公羊進行了基于IBS遺傳距離的聚類分析。分析結(jié)果（圖7B）揭示了這些公羊分屬于6個不同的家系，其中家系4和家系6各包含了3只公羊，而其他4個家系各由1只公羊單獨構(gòu)成。

2.6 ROH分析

99個子午嶺黑山羊個體共檢測到34 246個ROHs，其中33 960個（99.17%）ROHs長度位于0.25～1 Mb之間，284個（0.83%）ROHs長度＞1 Mb。由圖8A可知，個體水平上，個體ROH數(shù)量最多為1 023個，累計長度423.33 Mb，最少的僅有28個，累計長度8.31 Mb?；赗OH計算子午嶺黑山羊群體近交系數(shù)FROH，由圖8B可知，子午嶺黑山羊個體近交系數(shù)范圍為0.003 3～0.171 7，平均FROH值為0.049 6，其中公羊近交系數(shù)FROH值為0.051 1。

3 討 論

動物遺傳多樣性對整個物種的適應(yīng)性和生存能力至關(guān)重要，它增強了物種對環(huán)境變化的適應(yīng)力和抗壓能力［23-25］。其中，雜合度是衡量群體遺傳多樣性的一個重要參數(shù)。隨著群體雜合度的提高，表示該群體所包含的遺傳信息也更豐富［26］。期望雜合度（He）是指群體中任意個體在特定位點上呈現(xiàn)雜合狀態(tài)的概率，而觀測雜合度（Ho）則是指群體中在某位點上呈現(xiàn)雜合狀態(tài)的個體數(shù)量占總數(shù)量的比例［27］。當(dāng)He＞Ho時，說明群體內(nèi)部可能發(fā)生近交或受到了人工選擇的影響；當(dāng)He＜Ho時則可能意味著新遺傳變異的引入［28］。本研究顯示，子午嶺黑山羊群體的觀測雜合度為0.167，低于期望雜合度（0.192），說明群體內(nèi)部發(fā)生了近交。此外，該群體的多態(tài)信息含量平均值為0.161，屬于低度多態(tài)性［29-30］，反映出較低的遺傳變異水平，低于濟寧青山羊（0.283）［31］、永登七山羊（0.221）［27］、隴南山羊（0.236）［32］及柯爾克孜羊（0.737）［33］，說明不同群體的變異與生存環(huán)境及受選擇程度相關(guān)。這揭示了遺傳變異與群體所處環(huán)境及選擇壓力的緊密聯(lián)系。研究結(jié)果突出了子午嶺黑山羊群體遺傳多樣性的降低可能與群體內(nèi)部的近交現(xiàn)象密切相關(guān)。

本研究基于測序獲得的大量SNPs分別進行了基因組親緣關(guān)系G矩陣和IBS遺傳距離矩陣分析，旨在評估子午嶺黑山羊群體的親緣關(guān)系。G矩陣能更真實的反映個體間的親緣關(guān)系，它比單純根據(jù)系譜信息計算的親緣關(guān)系更加準(zhǔn)確且接近實際情況［34］。IBS只考慮個體之間的遺傳標(biāo)記或等位基因的相似性，因此在不清楚群體系譜或沒有祖代樣本的情況下，仍能對群體的親緣關(guān)系進行分析［35］。本研究中，親緣關(guān)系G矩陣及IBS遺傳距離矩陣分析均表明子午嶺黑山羊大部分個體間親緣關(guān)系較遠，但仍有少部分個體親緣關(guān)系較近。劉彬等［36］指出，在實際的育種實踐中，由于1頭種公畜可能與多頭母畜交配，對種公畜的選擇標(biāo)準(zhǔn)通常較為嚴(yán)格。這一現(xiàn)象導(dǎo)致不同種公畜間的親緣關(guān)系相對較遠，從而使得種公畜群體的IBS遺傳距離略高于群體平均水平。本研究結(jié)果與此觀點相吻合，子午嶺黑山羊公羊的平均IBS遺傳距離（0.833 6）略高于整個群體的平均遺傳距離（0.831 3）。動物群體結(jié)構(gòu)研究在遺傳學(xué)領(lǐng)域占據(jù)著舉足輕重的地位。深入探究物種的群體結(jié)構(gòu)和親緣關(guān)系對于物種的有效保護、資源的可持續(xù)利用以及生物多樣性的維護具有重大意義［37］。本研究通過主成分分析發(fā)現(xiàn)子午嶺黑山羊群體內(nèi)部并未表現(xiàn)出明顯的層次分化。盡管如此，子午嶺黑山羊部分個體間仍存在顯著的遺傳背景差異，這與子午嶺黑山羊的種群動態(tài)相契合（引入外血與部分個體雜交）。系統(tǒng)發(fā)育樹的構(gòu)建進一步證實了這一點，發(fā)現(xiàn)子午嶺黑山羊群體主要分布在兩個顯著的分支上。此外，對10只公羊的系統(tǒng)發(fā)育分析揭示了它們被劃分為6個不同的家系，兩個家系各包括3只公羊，其余家系均只有1只公羊，后期應(yīng)注意后代的選育，避免血統(tǒng)流失。

傳統(tǒng)的近交系數(shù)計算依賴于詳盡的家族譜系信息，其精確度和可信度很大程度上取決于譜系數(shù)據(jù)的質(zhì)量和完整性［34］。隨著測序技術(shù)的發(fā)展，基于全基因組SNPs數(shù)據(jù)直接計算近交系數(shù)的方法層出不窮，其中FROH可以在不依賴系譜數(shù)據(jù)的條件下準(zhǔn)確地計算每個個體基因組中同源的百分比，從而對群體的近交程度進行量化分析［21］。Liu等［38］對我國5個本地綿羊品種進行了ROH分析，發(fā)現(xiàn)大尾寒羊近交程度最高（FROH=0.080 8），而呼倫貝爾羊近交系數(shù)最低（FROH=0.014 8）。張任豹等［31］利用SNP芯片對我國濟寧青山羊進行基因分型，通過ROH分析發(fā)現(xiàn)濟寧青山羊群體的FROH是0.047 9。本研究中，子午嶺黑山羊群體的近交系數(shù)FROH為0.049 6，其中公羊近交系數(shù)FROH值為0.051 1，這一數(shù)值高于柯爾克孜羊（0.008）［33］，低于川中黑山羊（0.12）［39］，接近于濟寧青山羊（0.047 9）［31］，這表明子午嶺黑山羊的近交程度相對較低。但是后期隨著保種時間延長，加上群體規(guī)模局限性及相對封閉的管理模式，近交系數(shù)存在上升的趨勢。因此，在以后的的育種工作中應(yīng)持續(xù)加強譜系檔案的管理，以避免近親繁殖的發(fā)生。

4 結(jié) 論

本研究通過Super-GBS測序技術(shù)分析了子午嶺黑山羊群體遺傳多樣性、親緣關(guān)系及群體結(jié)構(gòu)，發(fā)現(xiàn)子午嶺黑山羊群體遺傳多樣性較低；大部分個體間親緣關(guān)系較遠，但仍有部分個體間親緣關(guān)系較近。子午嶺黑山羊群體未出現(xiàn)明顯分層，群體整體近交水平相對較低。為保種場進一步提純選育和優(yōu)化配種提供了理論依據(jù)。

參考文獻（References）：

［1］ NOTTER D R.The importance of genetic diversity in livestock populations of the future［J］.J Anim Sci，1999，77（1）：61-69.

［2］ EYDIVANDI S，ROUDBAR M A，KARIMI M O，et al.Genomic scans for selective sweeps through haplotype homozygosity and allelic fixation in 14 indigenous sheep breeds from Middle East and South Asia［J］.Sci Rep，2021，11（1）： 2834.

［3］ 茍想珍，常 嶸，文亞洲，等.子午嶺黑山羊種質(zhì)資源現(xiàn)狀及保護措施探討［J］.畜牧獸醫(yī)雜志，2023，42（4）：44-47，49.

GOU X Z，CHANG R，WEN Y Z，et al.Discussion on the current situation and protection measures of Ziwuling black goat germplasm resources［J］.Journal of Animal Science and Veterinary Medicine，2023，42（4）：44-47，49.（in Chinese）

［4］ MA K Y，SONG J J，LI D P，et al.Genomic structure analysis and construction of DNA fingerprint for four sheep populations［J］.Animal，2024，18（4）：101116.

［5］ MOKHBER M，MORADI-SHAHRBABAK M，SADEGHI M，et al.A genome-wide scan for signatures of selection in Azeri and khuzestani buffalo breeds［J］.BMC Genomics，2018，19（1）：449.

［6］ FU Y R，LI S K，GUO Q Q，et al.Genetic diversity and population structure of two endemic Cupressus （cupressaceae） species on the Qinghai-Tibetan plateau［J］.J Genet，2019，98（1）：14.

［7］ RELLSTAB C，ZOLLER S，WALTHERT L，et al.Signatures of local adaptation in candidate genes of oaks （Quercus spp.） with respect to present and future climatic conditions［J］.Mol Ecol，2016，25（23）：5907-5924.

［8］ 胡亞亞，劉蘭服，冀紅柳，等.簡化基因組測序技術(shù)研究進展［J］.江蘇師范大學(xué)學(xué)報（自然科學(xué)版），2018，36（4）：63-68.

HU Y Y，LIU L F，JI H L，et al.Research progress on the reduced-representation genome sequencing technique［J］.Journal of Jiangsu Normal University （Natural Science Edition），2018，36（4）：63-68.（in Chinese）

［9］ EMERSON K J，MERZ C R，CATCHEN J M，et al.Resolving postglacial phylogeography using high-throughput sequencing［J］.Proc Natl Acad Sci U S A，2010，107（37）：16196-16200.

［10］ HOHENLOHE P A，AMISH S J，CATCHEN J M，et al.Next-generation RAD sequencing identifies thousands of SNPs for assessing hybridization between rainbow and westslope cutthroat trout［J］.Mol Ecol Resour，2011，11（S1）：117-122.

［11］ HOHENLOHE P A，BASSHAM S，ETTER P D，et al.Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags［J］.PLoS Genet，2010，6（2）：e1000862.

［12］ 馬士龍，李小偉，李 響，等.基于GBS簡化基因組測序評估3個麥洼牦牛保種群的遺傳結(jié)構(gòu)研究［J］.草業(yè)學(xué)報，2022，31（9）：183-194.

MA S L，LI X W，LI X，et al.Assessment of genetic structure of 3 Maiwa yak preserved populations based on genotyping-by-sequencing technology［J］.Acta Prataculturae Sinica，2022，31（9）：183-194.（in Chinese）

［13］ 劉海林，江為民，段洪峰，等.基于GBS測序的湘東黑山羊的親緣關(guān)系及近交系數(shù)［J］.湖南農(nóng)業(yè)大學(xué)學(xué)報（自然科學(xué)版），2024，50（2）：78-85.

LIU H L，JIANG W M，DUAN H F，et al.Genetic relationship and inbreeding coefficient among Xiangdong black goats based on GBS sequencing［J］.Journal of Hunan Agricultural University （Natural Sciences），2024，50（2）：78-85.（in Chinese）

［14］ 徐建霞，丁延慶，馮 周，等.基于Super-GBS的高粱株高和節(jié)間數(shù)QTL定位［J］.生物技術(shù)通報，2023，39（7）：185-194.

XU J X，DING Y Q，F(xiàn)ENG Z，et al.QTL mapping of sorghum plant height and internode numbers based on Super-GBS technique［J］.Biotechnology Bulletin，2023，39（7）：185-194.（in Chinese）

［15］ LI H，DURBIN R.Fast and accurate long-read alignment with burrows-wheeler transform［J］.Bioinformatics，2010， 26（5）：589-595.

［16］ MCKENNA A，HANNA M，BANKS E，et al.The genome analysis toolkit：a mapreduce framework for analyzing next-generation DNA sequencing data［J］.Genome Res，2010，20（9）：1297-1303.

［17］ CINGOLANI P，PLATTS A，WANG L L，et al.A program for annotating and predicting the effects of single nucleotide polymorphisms，SnpEff：SNPs in the genome of drosophila melanogaster strain w1118;iso-2;iso-3［J］.Fly，2012，6（2）：80-92.

［18］ SAITOU N，NEI M.The neighbor-joining method：a new method for reconstructing phylogenetic trees［J］.Mol Biol Evol，1987，4（4）：406-425.

［19］ RUAN J，LI H，CHEN Z Z，et al.TreeFam：2008 update［J］.Nucleic Acids Res，2008，36（S1）：D735-D740.

［20］ VANRADEN P M.Efficient methods to compute genomic predictions［J］.J Dairy Sci，2008，91（11）：4414-4423.

［21］ MARRAS G，GASPA G，SORBOLINI S，et al.Analysis of runs of homozygosity and their relationship with inbreeding in five cattle breeds farmed in italy［J］.Anim Genet，2015，46（2）：110-121.

［22］ PERIPOLLI E，MUNARI D P，SILVA M V G B，et al.Runs of homozygosity：current knowledge and applications in livestock［J］.Anim Genet，2017，48（3）：255-271.

［23］ TU Z C，ZHANG Y P，QIU H.Genetic diversity and divergence in Chinese yak （Bos grunniens） populations inferred from blood protein electrophoresis［J］.Biochem Genet，1997，35（1/2）：13-16.

［24］ BAAZAOUI I，BEDHIAF-ROMDHANI S，MASTRANGELO S，et al.Genome-wide analyses reveal population structure and identify candidate genes associated with tail fatness in local sheep from a semi-arid area［J］.Animal，2021，15（4）：100193.

［25］ LVAREZ I，F(xiàn)ERNNDEZ I，TRAORé A，et al.Genomic scan of selective sweeps in Djallonk （West African Dwarf） sheep shed light on adaptation to harsh environments［J］.Sci Rep，2020，10（1）：2824.

［26］ MASTRANGELO S，CIANI E，SARDINA M T，et al.Runs of homozygosity reveal genome-wide autozygosity in Italian sheep breeds［J］.Anim Genet，2018，49（1）：71-81.

［27］ 馬克巖，韓金濤，白雅琴，等.基于簡化基因組測序的永登七山羊遺傳多樣性分析［J］.畜牧獸醫(yī)學(xué)報，2023， 54（5）：1939-1950.

MA K Y，HAN J T，BAI Y Q，et al.Genetic diversity analysis of Yongdeng Qishan sheep based on specific-locus amplified fragment sequencing［J］.Acta Veterinaria et Zootechnica Sinica，2023，54（5）：1939-1950.（in Chinese）

［28］ XU S S，REN X，YANG G L，et al.Genome-wide association analysis identifies the genetic basis of fat deposition in the tails of sheep （Ovis aries）［J］.Anim Genet，2017，48（5）：560-569.

［29］ ZHANG W G，GAO X，ZHANG Y，et al.Genome-wide assessment of genetic diversity and population structure insights into admixture and introgression in Chinese indigenous cattle［J］.BMC Genet，2018，19（1）：114.

［30］ AZIMU W，MANATBAY B，LI Y，et al.Genetic diversity and population structure analysis of eight local chicken breeds of southern xinjiang［J］.Br Poult Sci，2018，59（6）：629-635.

［31］ 張任豹，周東輝，周李生，等.基于70K SNP芯片分析濟寧青山羊保種群體的遺傳結(jié)構(gòu)［J］.畜牧獸醫(yī)學(xué)報，2023， 54（7）：2836-2847.

ZHANG R B，ZHOU D H，ZHOU L S，et al.Analysis of genetic structure of conservation population in Jining gray goats based on 70 K SNP chip［J］.Acta Veterinaria et Zootechnica Sinica，2023，54（7）：2836-2847.（in Chinese）

［32］ 馬克巖，劉占經(jīng)，白雅琴，等.基于簡化基因組測序評估隴南山羊群體遺傳多樣性和群體結(jié)構(gòu)［J］.中國畜禽種業(yè)，2024，20（4）：8-17.

MA K Y，LIU Z J，BAI Y Q，et al.Evaluation of genetic diversity and population structure of Longnan goats using SLAF-seq［J］.The Chinese Livestock and Poultry Breeding，2024，20（4）：8-17.（in Chinese）

［33］ 李隱俠，牙生江·那斯?fàn)?，賽里克·都曼，?SNP芯片評估柯爾克孜羊群體遺傳多樣性和遺傳結(jié)構(gòu)［J］.畜牧獸醫(yī)學(xué)報，2023， 54（2）：572-583.

LI Y X，YA S N，SAI L D，et al.Evaluation of genetic diversity and genetic structure in Kirgiz sheep population based on SNPs chip［J］.Acta Veterinaria et Zootechnica Sinica，2023，54（2）：572-583.（in Chinese）

［34］ SHI L Y，WANG L G，LIU J X，et al.Estimation of inbreeding and identification of regions under heavy selection based on runs of homozygosity in a large white pig population［J］.J Anim Sci Biotechnol，2020，11：46.

［35］ 劉晨龍，盧 丹，周泉勇，等.利用高密度SNP芯片分析杭豬的群體遺傳結(jié)構(gòu)［J］.畜牧獸醫(yī)學(xué)報，2022，53（8）：2502-2513.

LIU C L，LU D，ZHOU Q Y，et al.Analysis of population genetic structure of Hang pigs by high density SNP chip［J］.Acta Veterinaria et Zootechnica Sinica，2022，53（8）：2502-2513.（in Chinese）

［36］ 劉 彬，沈林園，陳 映，等.基于SNP芯片分析青峪豬保種群體的遺傳結(jié)構(gòu)［J］.畜牧獸醫(yī)學(xué)報，2020，51（2）：260-269.

LIU B，SHEN L Y，CHEN Y，et al.Analysis of genetic structure of conservation population in Qingyu pig based on SNP chip［J］.Acta Veterinaria et Zootechnica Sinica，2020，51（2）：260-269.（in Chinese）

［37］ GUO J Z，TAO H X，LI P F，et al.Whole-genome sequencing reveals selection signatures associated with important traits in six goat breeds［J］.Sci Rep，2018，8（1）：10405.

［38］ LIU J X，SHI L Y，LI Y，et al.Estimates of genomic inbreeding and identification of candidate regions that"" differ between Chinese indigenous sheep breeds［J］.J Anim Sci Biotechnol，2021，12（1）：95.

［39］ 郭家中，孫學(xué)良，向秋楠，等.川中黑山羊基因組近交系數(shù)和選擇信號分析［J］.四川農(nóng)業(yè)大學(xué)學(xué)報，2022，40（5）：656-663.

GUO J Z，SUN X L，XIANG Q N，et al.Analyses of genomic inbreeding coefficients and selection signals in Chuanzhong black goats［J］.Journal of Sichuan Agricultural University，2022，40（5）：656-663.（in Chinese）

（編輯 郭云雁）