基于重測(cè)序數(shù)據(jù)評(píng)估南陽(yáng)牛保種效果

摘 要： 旨在通過遺傳多樣性和群體結(jié)構(gòu)分析，評(píng)估南陽(yáng)牛的保種效果。本研究基于30頭健康南陽(yáng)牛的全基因組重測(cè)序數(shù)據(jù)，通過變異檢測(cè)獲得單核苷酸多態(tài)性（single nucleotide polymorphism， SNP）信息，綜合分析群體遺傳多樣性、群體結(jié)構(gòu)、親緣關(guān)系以及連續(xù)純和片段（ROH）分布特征，對(duì)南陽(yáng)牛保種群的保種效果進(jìn)行全面評(píng)估。結(jié)果顯示：1）共檢測(cè)到高質(zhì)量SNPs位點(diǎn)數(shù)目是25 929 389個(gè)；2）保種群遺傳多樣性豐富，核苷酸多樣性為0.002 9，多態(tài)性標(biāo)記比例為0.888 7，最小等位基因頻率為0.186 1，期望雜合度為0.274 9，觀測(cè)雜合度為0.255 7；3）世代有效群體含量在1 000代前為2 834頭，在20代前為149頭，呈逐年下降趨勢(shì)；4）主成分分析結(jié)果顯示南陽(yáng)牛保種群沒有明顯分層；5）個(gè)體間的遺傳距離介于0.80～0.89，平均遺傳距離為0.83±0.02；6）30頭南陽(yáng)牛個(gè)體共分為7個(gè)家系，公?？煞譃?個(gè)家系；7）共檢測(cè)到ROH片段數(shù)目是4 992個(gè)，平均每頭南陽(yáng)牛的ROH片段長(zhǎng)度為約74.41 Mb，總長(zhǎng)度為2.18 Gb，基于ROH的平均近交系數(shù)為0.031，0.5 Mb以下的ROH片段占比最多為75.72％，2～4 Mb的ROH片段占比最少為0.16％。綜上所述，南陽(yáng)牛保種群體遺傳多樣性豐富，沒有出現(xiàn)明顯分層，近交水平較低，但仍有個(gè)別個(gè)體近交較高。因此，未來的保種工作應(yīng)加強(qiáng)選種和選配管理，以促進(jìn)群體的可持續(xù)發(fā)展。

關(guān)鍵詞： 南陽(yáng)牛；全基因組重測(cè)序；遺傳多樣性；群體結(jié)構(gòu)；保種

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

文獻(xiàn)標(biāo)志碼：A 文章編號(hào)： 0366-6964（2024）09-3876-11

Evaluation of the Conservation Effect in Nanyang Cattle Based on Resequencing Data

LIU" Siyu1， ZHANG" Man1， ZHANG" Yan1， WEI" Zhitong1， QI" Xinglei2， GAO" Tengyun1， LIU" Xian3， LIANG" Dong1*， FU" Tong1*

（1.College of Animal Science and Technology， Henan Agricultural University，

Zhengzhou 450046，" China;2.Center of Animal Husbandry Technical Service in Biyang， Biyang 463700，" China; 3.Henan Provincial Animal Husbandry Station， Zhengzhou 450008，" China）

Abstract：" The study aimed to evaluate the conservation effect of Nanyang cattle by analyzing genetic diversity and population structure. In this study， based on the whole genome resequencing data of 30 healthy Nanyang cattle， single nucleotide polymorphism （SNP） information was obtained through mutation detection. The genetic diversity， population structure， phylogenetic relationship and runs of homozygosity （ROH） distribution characteristics were analyzed comprehensively to evaluate the conservation effect of Nanyang cattle. The results showed as follows： 1） The total number of high quality SNPs sites was 25 929 389; 2） The conservation population was rich in genetic diversity， nucleotide diversity was 0.002 9， polymorphism marker ratio was 0.888 7， minimum allele frequency was 0.186 1， expected heterozygosity was 0.274 9， observed heterozygosity was 0.255 7; 3） The effective population size was 2 834 before 1 000 generations and 149 before 20 generations， which showed a decreasing trend year by year; 4） The results of principal component analysis showed that there was no obvious stratification in Nanyang cattle conservation population; 5） The genetic distance between individuals ranged from 0.80 to 0.89， and the average genetic distance was 0.83±0.02; 6） The 30 Nanyang cattle were divided into 7 families， and the bulls were divided into 5 families; 7） A total of 4 992 ROH fragments were detected， the average length of ROH fragments per Nanyang cattle was about 74.41 Mb， the total length was 2.18 Gb， the average inbreeding coefficient based on ROH was 0.031， and the maximum proportion of ROH fragments below 0.5 Mb was 75.72%， the proportion of ROH fragments of 2-4 Mb was at least 0.16%. In summary， the genetic diversity of the Nanyang cattle conservation population is rich， there is no obvious stratification， the inbreeding level is low， but there are still some individuals with high inbreeding. Therefore， the future conservation work should strengthen the management of seed selection and mating， in order to promote the sustainable development of the population.

Key words： Nanyang cattle; whole genome resequencing; genetic diversity; population structure; breed conservation

*Corresponding authors：LIANG Dong， E-mail：554230517@qq.com; FU Tong， E-mail： futong2004@126.com

南陽(yáng)牛產(chǎn)于河南省南陽(yáng)市，是我國(guó)五大黃牛品種之一，具有體軀高大、結(jié)構(gòu)勻稱、役肉兼用、耐粗飼、適應(yīng)性強(qiáng)、性情溫順的特征。1988年將南陽(yáng)牛收錄于《中國(guó)牛品種志》，2000年，農(nóng)業(yè)部將南陽(yáng)牛列入《國(guó)家畜禽品種保護(hù)名錄》，2006年列入《國(guó)家畜禽遺傳資源保護(hù)名錄》，標(biāo)志著對(duì)其重要性的認(rèn)可［1-2］。2002年，國(guó)家質(zhì)量技術(shù)監(jiān)督總局對(duì)南陽(yáng)牛進(jìn)行了原產(chǎn)地標(biāo)記域名注冊(cè)，為南陽(yáng)牛產(chǎn)業(yè)的快速發(fā)展奠定了基礎(chǔ)。南陽(yáng)牛作為中國(guó)著名的地方優(yōu)良品種之一，在畜牧業(yè)中占有重要地位［3］。但隨著社會(huì)上肉牛雜交改良力度的加大，南陽(yáng)牛血統(tǒng)愈來愈混雜，純種南陽(yáng)牛數(shù)量急劇下降，使得南陽(yáng)牛保種及開發(fā)利用工作顯得尤為迫切［4］。南陽(yáng)市采取了積極措施，設(shè)立了南陽(yáng)市黃牛良種繁育場(chǎng)和南陽(yáng)黃?？萍贾行模诵娜罕３?個(gè)以上家系，并完善系譜和技術(shù)管理檔案。同時(shí)，河南省畜牧部門還通過與夏洛萊、皮埃蒙特、德國(guó)黃牛等品種進(jìn)行改良，成功培育出了新品種夏南牛、皮南牛和德南牛，為南陽(yáng)牛的保護(hù)和開發(fā)利用做出了重要貢獻(xiàn)［5-7］。因此，為了加強(qiáng)對(duì)南陽(yáng)牛的保種和開發(fā)利用，不使地方肉牛資源出現(xiàn)缺失，南陽(yáng)牛的保種任務(wù)已經(jīng)迫在眉睫。

由于我國(guó)南陽(yáng)牛保種工作起步較晚，技術(shù)、資金投入有限，缺乏系統(tǒng)的數(shù)據(jù)記錄以及實(shí)際生產(chǎn)中易出現(xiàn)記錄錯(cuò)誤等原因，傳統(tǒng)基于系譜的遺傳結(jié)構(gòu)分析容易受到缺失記錄和錯(cuò)誤記錄的干擾，難以得出準(zhǔn)確結(jié)論。作為一種高密度、多態(tài)性強(qiáng)且遺傳上穩(wěn)定的技術(shù)，單核苷酸多態(tài)性（single nucleotide polymorphism， SNP）已成為研究畜禽遺傳變異的關(guān)鍵工具［8］。SNP分型主要通過SNP芯片、簡(jiǎn)化基因組測(cè)序及全基因組重測(cè)序（whole genome sequencing， WGS）等方法，廣泛用于分析群體結(jié)構(gòu)、進(jìn)行全基因組選擇和評(píng)估種質(zhì)資源［9-10］。WGS技術(shù)的成熟和成本降低促進(jìn)了近年來對(duì)牛［11］、豬［12］、馬［13］、羊［14］、禽［15］等物種群體結(jié)構(gòu)和遺傳多樣性研究的快速發(fā)展。近年來，已有多項(xiàng)研究利用WGS技術(shù)對(duì)南陽(yáng)牛進(jìn)行了全基因組重測(cè)序分析，為深入了解其血統(tǒng)組成和遺傳多樣性提供了重要見解［16-18］。但是，針對(duì)南陽(yáng)牛保種效果進(jìn)行全方面評(píng)估尚未見報(bào)道。

本研究聚焦南陽(yáng)牛保種群體，采用全基因組重測(cè)序技術(shù)深入分析其群體結(jié)構(gòu)和遺傳多樣性。通過評(píng)估保種效果，為南陽(yáng)牛遺傳資源的有效保護(hù)和可持續(xù)開發(fā)利用提供科學(xué)的理論依據(jù)。本研究將深化對(duì)南陽(yáng)牛遺傳特性和群體動(dòng)態(tài)的理解，為其遺傳管理和持續(xù)利用提供實(shí)質(zhì)性指導(dǎo)。

1 材料與方法

1.1 試驗(yàn)材料

本試驗(yàn)選取河南省南陽(yáng)市黃牛良種繁育場(chǎng)30頭南陽(yáng)牛（公牛12頭，母牛18頭）作為研究對(duì)象。所有試驗(yàn)牛均在相同的設(shè)施條件和環(huán)境下的圈舍中飼養(yǎng)，并執(zhí)行相同的飼養(yǎng)標(biāo)準(zhǔn)。在采樣時(shí)，這些牛均是健康狀態(tài)。所有樣本的耳組織被采集并存放在含95%酒精的2 mL凍存管中，-20 ℃保存。

1.2 試驗(yàn)方法

1.2.1 DNA提取及測(cè)序

使用南京諾唯贊生物科技股份有限公司的組織基因組DNA提取試劑盒（DC112）提取南陽(yáng)牛耳組織中的基因組DNA，并利用瓊脂糖凝膠電泳及紫外分光光度計(jì)評(píng)估其完整性、濃度與產(chǎn)量。使用DNBSEQ-T7平臺(tái)對(duì)每個(gè)個(gè)體的基因組DNA進(jìn)行150 bp雙末端reads的測(cè)序。Trimomatic［19］軟件用于過濾fastq數(shù)據(jù)，生成clean reads。Clean reads通過BWA-MEM［20］算法與參考基因組（ARS-UCD1.3）進(jìn)行比對(duì)。使用Picard［21］的MarkDuplicates模塊和GATK（v 3.8）［22］的IndelRealigner模塊分別去除PCR重復(fù)和重新比對(duì)indels。變異檢測(cè)采用UnifiedGenotyper模塊，過濾參數(shù)為：“QDlt;2.0， FSgt;60.0， 等。使用VCFtools（v 0.1.16）［23］將變異數(shù)據(jù)從VCF格式轉(zhuǎn)為Plink格式。質(zhì)量控制中，本研究使用plink移除了最小等位基因頻率小于0.05、Hardy-Weinberg平衡P值小于10-6的SNP位點(diǎn)和缺失基因型大于10%的個(gè)體。

1.2.2 遺傳多樣性分析

本研究對(duì)SNP數(shù)據(jù)進(jìn)行基因功能注釋所使用到的軟件是SnpEff［24］。主要是通過PLINK（v 1.90）［25］軟件計(jì)算南陽(yáng)牛的最小等位基因頻率（minor allele frequency， MAF）、期望雜合度（expected heterozygosity， He）、觀測(cè)雜合度（observed heterozygosity， Ho）和多態(tài)性標(biāo)記比例（proportion of polymorphic marker， PN），評(píng)估保種群的遺傳多樣性。VCFtools（v 0.1.16）軟件計(jì)算核苷酸多樣性（π），設(shè)置滑動(dòng)窗口50 kb和步長(zhǎng)20 kb。SNeP（v 1.1）［26］軟件基于連鎖不平衡計(jì)算歷史有效群體大?。∟e）。

1.2.3 南陽(yáng)牛群體結(jié)構(gòu)分析

本研究利用PLINK（v 1.90）軟件進(jìn)行主成分分析（PCA），并用R語(yǔ)言的ggplot2包對(duì)前兩主成分進(jìn)行可視化。通過PLINK計(jì)算遺傳距離，構(gòu)建IBS距離矩陣和tassel［27］軟件構(gòu)建G矩陣，R語(yǔ)言進(jìn)行結(jié)果可視化。再次使用PLINK計(jì)算距離矩陣，通過R包“ape”、“phangorn”、“seqinr”建立NJ系統(tǒng)發(fā)育樹，利用iTOL（https：//itol.embl.de/）網(wǎng)站進(jìn)行可視化。

1.2.4 全基因組ROH統(tǒng)計(jì)分析

ROH檢測(cè)采用PLINK（v 1.90）［25］的-homozyg選項(xiàng)，以0.5～1 Mb、1～2 Mb、2～4 Mb、gt;4 Mb將其ROH分為4個(gè)等級(jí)并計(jì)算近交系數(shù)（FROH），通過R語(yǔ)言軟件包ggplot2可視化最終結(jié)果。

2 結(jié) 果

2.1 南陽(yáng)牛群體遺傳多樣性分析

本研究對(duì)30頭南陽(yáng)牛進(jìn)行了全基因組重測(cè)序，平均測(cè)序深度達(dá)到12.65×。產(chǎn)生的clean reads有6 972 566 582個(gè)，其中平均99.87％的clean reads可以比對(duì)到參考基因組上。經(jīng)過變異檢測(cè)和質(zhì)量控制，30頭南陽(yáng)牛共獲得25 929 389個(gè)高質(zhì)量SNPs。在常染色體中，1號(hào)染色體上的SNP數(shù)目最多，包含1 897 811個(gè)SNPs；25號(hào)染色體上的SNP數(shù)目最少，包含502 704個(gè)SNPs（圖1）?；蚬δ茏⑨尳Y(jié)果顯示（圖2），南陽(yáng)牛全基因組中SNPs主要分布在基因間區(qū)（73 677 282個(gè)，占48.89％）以及內(nèi)含子區(qū)（57 520 323個(gè)，占38.17％）。同義突變與錯(cuò)義突變所占比例較少，分別包括462 843與209 938個(gè)SNPs。

南陽(yáng)牛保種群遺傳多樣性的參數(shù)，包括最小等位基因頻率0.186 1，多態(tài)性標(biāo)記比例0.888 7，期望雜合度0.274 9，觀測(cè)雜合度0.255 7，以及核苷酸多樣性為0.002 9，詳細(xì)信息見表1。觀測(cè)雜合度略低于期望雜合度，差異不顯著。

南陽(yáng)牛保種群歷史有效群體大小的結(jié)果顯示，在1 000代前時(shí)，Ne為2 834頭；在150代前時(shí)，Ne為689頭；在20代前時(shí)，Ne為149頭（圖3）。這種變化趨勢(shì)說明南陽(yáng)牛保種群的有效群體大小逐代減少。

2.2 南陽(yáng)牛群體結(jié)構(gòu)分析

主成分分析結(jié)果顯示，主成分1解釋了16.8%的遺傳變異，主成分2解釋了13.96%的遺傳變異。保種群內(nèi)部，12頭公牛（NY-M）分布均勻，未表現(xiàn)出明顯的聚集趨勢(shì)，見圖4所示。IBS的遺傳距離矩陣顯示，30頭南陽(yáng)牛個(gè)體間的IBS遺傳距離范圍為0.80～0.89，平均值為0.83±0.02，見圖5。G矩陣揭示的親緣關(guān)系如圖6，通過顏色深淺表示，親緣關(guān)系系數(shù)接近1時(shí)顏色較淺，表明個(gè)體間關(guān)系較近。大多數(shù)南陽(yáng)牛個(gè)體顯示中到低等級(jí)的親緣關(guān)系，僅少數(shù)表現(xiàn)出較高的親緣度。南陽(yáng)牛保種群的"" 遺傳距離和親緣關(guān)系分析表明，個(gè)體間遺傳距離較大，近交程度低，但不排除潛在的近交風(fēng)險(xiǎn)。

使用NJ法［28］為南陽(yáng)牛保種群體構(gòu)建了遺傳距離矩陣，并將親緣關(guān)系系數(shù)超過0.1的母牛歸入相應(yīng)的公牛家系中，結(jié)果顯示30頭南陽(yáng)牛可分為7個(gè)家系（圖7），而12頭公牛分為5個(gè)家系（圖8）。家系劃分如下：家系A(chǔ)包含3頭公牛；家系B包含2頭公牛和5頭母牛；家系C包含6頭母牛；家系D包含1頭公牛和3頭母牛；家系E包含2頭母牛；家系F包含2頭公牛；家系G包含4頭公牛和2頭母牛（表2）。

2.3 南陽(yáng)牛群體基因組ROH及近交系數(shù)分析結(jié)果

30頭南陽(yáng)牛中發(fā)現(xiàn)了4 992個(gè)ROH片段，平均每頭南陽(yáng)牛約有166.4個(gè)ROH片段。所有ROH片段的總長(zhǎng)度達(dá)2.18 Gb，平均每頭南陽(yáng)牛的ROH片段長(zhǎng)度約74.41 Mb。在所有染色體中，第7號(hào)染色體的ROH片段最多，共336個(gè)；而第25號(hào)染色體的最少，僅54個(gè)（圖9A）。對(duì)南陽(yáng)牛ROH片段進(jìn)行分類統(tǒng)計(jì)發(fā)現(xiàn)，0～0.5 Mb長(zhǎng)度的片段最多，共3 780個(gè)，占總數(shù)的75.72%；而2～4 Mb長(zhǎng)度的片段最少，共8個(gè)，僅占0.16%（圖9B）。南陽(yáng)牛保種群的ROH長(zhǎng)度和數(shù)量表明了其豐富的遺傳多樣性。30頭南陽(yáng)牛ROH的近交系數(shù)介于0.012～0.058之間，平均值為0.031（圖9C）。

3 討 論

畜禽遺傳資源保護(hù)致力于維護(hù)遺傳多樣性，對(duì)國(guó)家畜產(chǎn)品供應(yīng)和種業(yè)振興至關(guān)重要。中國(guó)是畜種業(yè)強(qiáng)國(guó)，培育了眾多牛品種，其中南陽(yáng)牛作為杰出的地方品種，因其體型龐大、耐力強(qiáng)、適應(yīng)粗飼、肉質(zhì)鮮美及其特有的大理石紋理，在畜牧業(yè)發(fā)揮著關(guān)鍵作用，對(duì)肉牛產(chǎn)業(yè)和新品種培養(yǎng)極為重要。2000年，南陽(yáng)牛被納入國(guó)家畜禽品種保護(hù)名錄，政府從此加強(qiáng)保護(hù)并開展了開發(fā)與選育。自1986年起，南陽(yáng)牛數(shù)量不斷增加，一度達(dá)到約240萬頭。但隨著農(nóng)業(yè)機(jī)械化的提高和外來品種的引入，自1997年起其數(shù)量開始逐年減少。為保護(hù)這一品質(zhì)資源，南陽(yáng)市建立了南陽(yáng)市黃牛良種繁育場(chǎng)（原南陽(yáng)黃牛場(chǎng)）和南陽(yáng)黃牛科技中心（原南陽(yáng)黃牛研究所），對(duì)南陽(yáng)黃牛進(jìn)行品種保護(hù)和改良利用。本研究通過對(duì)30頭南陽(yáng)牛進(jìn)行全基因組重測(cè)序，共獲取了25 929 389個(gè)高質(zhì)量SNPs位點(diǎn)。利用這些數(shù)據(jù)評(píng)估了南陽(yáng)牛的遺傳多樣性和群體結(jié)構(gòu)，為保種工作提供了科學(xué)依據(jù)。

群體遺傳多樣性既映射了物種的進(jìn)化軌跡，也展示了其對(duì)環(huán)境的適應(yīng)性［29-32］。較高的遺傳多樣性通常預(yù)示著更強(qiáng)的環(huán)境適應(yīng)力，對(duì)于促進(jìn)遺傳資源的持續(xù)利用具有重要意義。Mei等［16］的研究表明，南陽(yáng)牛群體的遺傳結(jié)構(gòu)呈現(xiàn)出多樣性，并且在全球牛種中展現(xiàn)出獨(dú)特的遺傳特征；另一項(xiàng)研究由Zhang等［17］進(jìn)行，發(fā)現(xiàn)南陽(yáng)牛中存在著豐富的遺傳多樣性，并且呈現(xiàn)出一定程度的品種特異性；此外，Chen等［18］的研究進(jìn)一步揭示了南陽(yáng)牛的遺傳背景，指出其可能受到來自不同地區(qū)和品種的遺傳影響，表現(xiàn)出復(fù)雜的遺傳結(jié)構(gòu)。這些研究結(jié)果表明，南陽(yáng)牛的遺傳多樣性和血統(tǒng)組成受到多種因素的影響，包括地理環(huán)境、人工選擇和品種間交流等。根據(jù)已有研究，南陽(yáng)牛是黃牛和瘤牛的雜交品種，其具有較高的遺傳多樣性［17］。未來的研究可以進(jìn)一步探索南陽(yáng)牛群體的遺傳背景和演化歷史，以更好地指導(dǎo)其保種工作和遺傳資源管理。

本研究對(duì)南陽(yáng)牛群體的遺傳多樣性進(jìn)行了深入分析。結(jié)果顯示，南陽(yáng)牛群體的平均最小等位基因頻率為0.186 1，指示其遺傳多樣性較高。同時(shí)，58%的SNPs位點(diǎn)多態(tài)性進(jìn)一步證實(shí)了這一點(diǎn)。南陽(yáng)牛的平均觀測(cè)雜合度為0.255 7，略低于期望值，暗示可能存在輕微的選擇壓力和近交現(xiàn)象。通過連鎖不平衡估計(jì)，20代前南陽(yáng)牛的有效群體規(guī)模約為149頭，顯示下降趨勢(shì)，這說明南陽(yáng)牛群體仍保持較高的遺傳多樣性，但仍然還存在遺傳多樣性喪失的風(fēng)險(xiǎn)。根據(jù)本研究結(jié)果，未來遺傳管理和持續(xù)利用工作可以考慮以下改進(jìn)措施：制定科學(xué)的繁殖計(jì)劃，避免近親繁殖，選擇遺傳多樣性高的個(gè)體進(jìn)行配種。收集南陽(yáng)牛的遺傳信息，為種群管理提供數(shù)據(jù)支持。在保持遺傳多樣性的同時(shí)，通過人工選擇引導(dǎo)種群進(jìn)化。提高養(yǎng)殖戶對(duì)遺傳多樣性保護(hù)和繁殖管理的意識(shí)和能力，確保管理措施的有效實(shí)施。

為防止外來品種干擾，保種場(chǎng)多采用內(nèi)群繁育，使得研究保種群體結(jié)構(gòu)對(duì)于其可持續(xù)發(fā)展極為重要［33-34］。主成分分析顯示，30頭南陽(yáng)牛沒有明顯分層，遺傳距離分析表明個(gè)體間親緣關(guān)系較遠(yuǎn)，暗示近交程度較低。聚類分析將30頭南陽(yáng)牛分為7個(gè)家系，顯示了家系間的個(gè)體數(shù)量差異。未來保種應(yīng)關(guān)注公牛分布，優(yōu)化配種策略以保證家系平衡。ROH產(chǎn)生于親代將同源單倍型完整傳遞給后代。長(zhǎng)ROH片段表明近期的近交歷史，而短ROH源自遠(yuǎn)祖。研究發(fā)現(xiàn)，30頭南陽(yáng)牛中共有4" 992個(gè)ROH片段，總長(zhǎng)2.18 Gb。其中，0～0.5 Mb的ROH片段最多，而2～4 Mb的最少。這說明南陽(yáng)牛近期無顯著近交。平均近交系數(shù)為0.031，反映了南陽(yáng)牛整體近交水平較低。為保持低近交水平，選配時(shí)應(yīng)注意高近交個(gè)體，引入新血統(tǒng)或遠(yuǎn)親家系的公牛進(jìn)行配種可以有效降低其后代的近交系數(shù)。也可以通過分子標(biāo)記技術(shù)識(shí)別近交個(gè)體，從而在選種和配種時(shí)予以排除。

4 結(jié) 論

本研究利用全基因組重測(cè)序數(shù)據(jù)分析了30頭南陽(yáng)牛的遺傳多樣性和群體結(jié)構(gòu)。發(fā)現(xiàn)南陽(yáng)牛保種群體遺傳多樣性較高，群體結(jié)構(gòu)未顯著分層。通過聚類分析，將這些牛分為7個(gè)家系，揭示了它們之間的親緣關(guān)系。雖然大部分南陽(yáng)牛近交水平較低，但仍有個(gè)別個(gè)體近交較高。因此，未來的保種工作應(yīng)加強(qiáng)選種和配種管理，以促進(jìn)群體的可持續(xù)發(fā)展。

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

［1］ XIA X，QU K，ZHANG G，et al.Comprehensive analysis of the mitochondrial DNA diversity in Chinese cattle［J］.Anim Genet，2019，50（1）：70-73.

［2］ ZHANG Y，WEI Z T，ZHANG M， et al.Population structure and selection signal analysis of Nanyang cattle based on whole-genome sequencing data［J］.Genes，2024，15（3）： 351.

［3］ XIA X，YAO Y，LI C，et al.Genetic diversity of Chinese cattle revealed by Y-SNP and Y-STR markers［J］.Anim Genet，2019，50（1）： 64-69.

［4］ 陳希鵑，肖喜東，吳明安.南陽(yáng)黃牛的保種、改良與種公牛飼養(yǎng)管理［J］.養(yǎng)殖與飼料，2020（3）：47-48.

CHEN X J，XIAO X D，WU M A.Breeding，improvement and breeding management of Nanyang yellow cattle［J］.Breeding and Feed，2020（3）：47-48.（in Chinese）

［5］ 杜書增，王冠立，王玉海，等.南陽(yáng)牛肉用選育改良及其產(chǎn)業(yè)化開發(fā)［J］.中國(guó)牛業(yè)科學(xué)，2016，42（2）：63-66.

DU S Z，WANG G L，WANG Y H，et al.Breeding and improvement of Nanyang cattle and its beed industrialization［J］.China Cattle Science，2016，42（2）：63-66.（in Chinese）

［6］ 歲豐軍.南陽(yáng)牛保種育種現(xiàn)狀［C］//《第八屆中國(guó)牛業(yè)發(fā)展大會(huì)》論文集.昌吉：中國(guó)畜牧業(yè)協(xié)會(huì)，2013：248-251.

SUI F J.Current situation of Nanyang cattle breeding［C］//Proceedings of the 8th China Cattle Industry Development Conference. Changji：China Animal Husbandry Association，2013：248-251.（in Chinese）

［7］ 茹寶瑞，高騰云.關(guān)于南陽(yáng)牛育種與生產(chǎn)的一些建議［C］//《第六屆中國(guó)牛業(yè)發(fā)展大會(huì)》論文集.重慶：中國(guó)畜牧業(yè)協(xié)會(huì)牛業(yè)分會(huì)，2011：310-313.

RU B R，GAO T Y.Some suggestions on breeding and production of Nanyang cattle［C］//Proceedings of the 6th China Cattle Industry Development Conference.Chongqing：China Animal Husbandry Association Cattle Industry Branch，2011：310-313.（in Chinese）

［8］ LIANG D，ZHAO P J，SI J F，et al.Genomic analysis revealed a convergent evolution of LINE-1 in coat color：a case study in water buffaloes （Bubalus bubalis）［J］.Mol Biol Evol，2021，38（3）：1122-1136.

［9］ KARIMI K，KOSHKOIYEH A E，F(xiàn)OZI M A，et al.Prioritization for conservation of Iranian native cattle breeds based on genome-wide SNP data［J］.Conserv Genet，2016，17（1）：77-89.

［10］ 劉繼強(qiáng)，郝曉東，武麗娜，等.全基因組SNP分型技術(shù)在畜禽遺傳育種研究中的應(yīng)用［J］.畜牧獸醫(yī)學(xué)報(bào)，2022， 53（12）： 4123-4137.

LIU J Q，HAO X D，WU L N，et al.Application of whole genome SNP genotyping technology in livestock and poultry genetics and breeding［J］.Acta Veterinaria et Zootechnica Sinica，2022，53（12）：4123-4137.（in Chinese）

［11］ HU M Y，SHI L L，YI W F，et al.Identification of genomic diversity and selection signatures in Luxi cattle using whole-genome sequencing data［J］.Anim Biosci，2024，37（3）：461-470.

［12］ 龍 熙，柴 捷，潘紅梅，等.基于SNP芯片的盆周山地豬群體選擇信號(hào)分析［J］.中國(guó)畜牧雜志，2023，59（12）：136-140.

LONG X，CHAI J，PAN H M，et al.Signal analysis of population selection of mountainous pigs in pelvic region based on SNP chip［J］.Chinese Journal of Animal Science，2023，59（12）：136-140.（in Chinese）

［13］ JAGANNATHAN V，GERBER V，RIEDER S，et al.Comprehensive characterization of horse genome variation by whole-genome sequencing of 88 horses［J］.Anim Genet，2019，50（1）：74-77.

［14］ 史露露，胡明月，賴偉寧，等.基于50K SNP芯片的夏洛來羊遺傳結(jié)構(gòu)及選擇信號(hào)分析［J］.黑龍江畜牧獸醫(yī)，2023（6）： 49-53，59，139-140.

SHI L L，HU M Y，LAI W N，et al.Genetic structure and selection signal analysis of Charolais sheep based on 50K SNP chip［J］.Heilongjiang Animal Science and Veterinary Medicine，2023（6）：49-53，59，139-140.（in Chinese）

［15］ GUO Y，RUBIN G J，R NNEBURG T，et al.Whole-genome selective sweep analyses identifies the region and candidate gene associated with white earlobe color in Mediterranean chickens［J］.Poult Sci，2024，103（1）：103232.

［16］ MEI C G，WANG H C，LIAO Q J，et al.Genetic architecture and selection of Chinese cattle revealed by whole genome resequencing［J］.Mol Biol Evol，2018，35（3）：688-699.

［17］ ZHANG S J，YAO Z，LI X M，et al.Assessing genomic diversity and signatures of selection in Pinan cattle using whole-genome sequencing data［J］.BMC Genomics，2022，23（1）：460.

［18］ 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.

［19］ BOLGER A M，LOHSE M，USADEL B.Trimmomatic：a flexible trimmer for Illumina sequence data［J］.Bioinformatics， 2014，30（15）：2114-2120.

［20］ LI H. Aligning sequence reads， clone sequences and assembly contigs with BWA-MEM［J］. arXiv preprint，2013，arXiv： 1303.3997.

［21］ EBBERT M T W，WADSWORTH M E，STALEY L A，et al.Evaluating the necessity of PCR duplicate removal from next-generation sequencing data and a comparison of approaches［J］.BMC Bioinformatics，2016，7：239.

［22］ 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.

［23］ DANECEK P，AUTON A，ABECASIS G，et al.The variant call format and VCFtools［J］.Bioinformatics，2011，27（15）：2156-2158.

［24］ CINGOLANI P.Variant annotation and functional prediction：SnpEff［J］.Methods Mol Biol，2022，2493：289-314.

［25］ 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.

［26］ BARBATO M，OROZCO-TERWENGEL P，TAPIO M，et al.SNeP：a tool to estimate trends in recent effective population size trajectories using genome-wide SNP data［J］.Front Genet，2015，6：109.

［27］ BRADBURY P J，ZHANG Z W，KROON D E，et al.TASSEL：software for association mapping of complex traits in diverse samples［J］.Bioinformatics，2007，23（19）：2633-2635.

［28］ WANG Y，DONG R L，LI X，et al.Analysis of the genetic diversity and family structure of the Licha Black pig population on Jiaodong Peninsula，Shandong Province，China［J］.Animals （Basel），2022，12（8）：1045.

［29］ EDWARDS C E，TESSIER B C，SWIFT J F，et al.Conservation genetics of the threatened plant species Physaria filiformis （Missouri bladderpod） reveals strong genetic structure and a possible cryptic species［J］.PLoS One，2021，16（3）：e0247586.

［30］ 劉晨龍，盧 丹，周泉勇，等.利用高密度SNP芯片分析杭豬的群體遺傳結(jié)構(gòu)［J］.畜牧獸醫(yī)學(xué)報(bào)，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）

［31］ 陳國(guó)宏，季從亮，王敏強(qiáng).12個(gè)中國(guó)地方雞種群體遺傳結(jié)構(gòu)及遺傳多樣性分析［J］.中國(guó)畜牧獸醫(yī)文摘，2006（5）：32.

CHEN G H，JI C L，WANG M Q.Analysis of population genetic structure and genetic diversity of 12 local chicken breeds in China［J］.Chinese Journal of Animal Science and Veterinary Abstracts，2006（5）：32.（in Chinese）

［32］ 李隱俠，牙生江·那斯?fàn)枺惱锟恕ざ悸?，?SNP芯片評(píng)估柯爾克孜羊群體遺傳多樣性和遺傳結(jié)構(gòu)［J］.畜牧獸醫(yī)學(xué)報(bào)，2023，54（2）：572-583.

LI Y X，YA SHENGJIANG·NASIER，SAI LIKE·DUMAN，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）

［33］ 馬克巖，韓金濤，白雅琴，等.基于簡(jiǎn)化基因組測(cè)序的永登七山羊遺傳多樣性分析［J］.畜牧獸醫(yī)學(xué)報(bào)，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）

［34］ 馬 鈞，樊安平，王武生，等.全基因組重測(cè)序解析秦川牛保種群遺傳多樣性和遺傳結(jié)構(gòu)［J］.遺傳，2023，45（7）：602-616.

MA J，F(xiàn)AN A P，WANG W S，et al.Analysis of genetic diversity and genetic structure of Qinchuan cattle conservation population using whole-genome resequencing［J］.Hereditas，2023，45（7）：602-616.（in Chinese）

（編輯 郭云雁）