基于SNP 芯片對里岔黑豬進行遺傳多樣性與遺傳結構分析的研究

李 曉，崔 超，王 源，于光輝*

（1.青島農業大學動物科技學院，山東青島 266109；2.膠州市農業農村局，山東青島 266300）

里岔黑豬是膠州市里岔鎮群眾四千年選育逐步形成的具有特色的、極為罕見的優良地方豬種，是山東省的重要地方品種豬[1-2]，其毛色全黑，嘴筒長直，額有縱紋，具有耐粗飼、適應能力強、繁殖性能好等特色，椎骨數較其他豬種多1～2 節，是我國優良的地方品種，素有“南太湖，北里岔”之稱[3]。

自上世紀以來，國內引種熱潮高漲，大量國外豬種（長白、大白、杜洛克等）被引入，由于引入豬種具有瘦肉率高、生長周期短等特點[4]，對中國地方豬造成較大沖擊，且部分地方豬種被盲目雜交或者被淘汰，導致血統不純和群體規模銳減[5]，進行種質資源保護顯得尤為重要。

以人工記錄系譜的方式進行遺傳結構分析，會因記錄錯誤等原因，影響最終選種育種效果。通過分子水平進行群體遺傳結構分析能夠精準了解里岔黑豬群體現狀，對地方豬種保護和利用具有重要意義[6]。隨著保種技術的不斷發展，通過SNP 芯片保種逐漸成為主流[7-8]。SNP 芯片作為現階段重要的育種工具之一，被廣泛應用于各個領域，其工作原理是通過芯片上固定的DNA序列與目的基因發生堿基互補配對，繼而進行序列中的SNP 位點鑒定[9-11]。本實驗利用SNP 芯片對茂華里岔黑豬良種場的在繁群體進行群體遺傳結構和多樣性分析，為里岔黑豬的選育及品種利用率提高等方面提供分子水平的參考依據。目前對里岔黑豬群體遺傳結構的研究處于一個相對初級階段，通過系統評價該群體有利于更好地保護該遺傳資源。

1 材料與方法

1.1 實驗材料 選擇膠州市茂華里岔黑豬良種場共109頭能繁里岔黑豬為研究對象，其中公豬13 頭、母豬96 頭，采集其耳組織置于含無水乙醇的離心管中帶回實驗室。

1.2 實驗方法

1.2.1 DNA 提取與基因分型 采用酚-氯仿抽提法提取耳組織DNA，經紫外分光光度儀和凝膠電泳對DNA 質量進行檢測，檢測合格后的DNA 樣本使用“中芯一號”芯片（北京康普森農業科技有限公司）進行SNP 基因分型，然后利用Plink（V1.90）軟件[12]對基因型數據進行質控，保留常染色體上的位點，剔除最小等位基因頻率小于0.01 的位點，質控后剩余39 644 個SNP 位點。

1.2.2 親緣關系分析 將Plink（V1.90）軟件質控后的數據進行主成分分析并計算個體間的親緣關系[13]，采用Origin 2021 軟件繪制第一主成分與第二主成分的二維圖，分析品種內的群體結構，使用R 腳本構建熱圖。通過G 矩陣（V2）軟件和RStudio 來計算親緣關系值和熱圖。

1.2.3 群體結構分析 將基因型數據質控后的39 644 個SNP 位點進行群體遺傳結構分析，通過MEGA7.0 軟件[14]構建鄰接法進化發育樹（Neighbor-Joining tree,NJ tree）進行聚類分析，通過Plink（V1.90）軟件進行群體結構分析。

1.2.4 ROH 分析 ROH 的長度和頻率用來反應該群體的歷史，使用Plink（V1.90）軟件分析并統計每個樣本中ROH 長度、分布和數目，利用公式計算ROH 片段總長度與常染色體基因總長度比值，得到基于ROH 的近交系數[15]。

其中，k 是個體的ROH 數目，L 是該物種的常染色體基因組的長度（豬V10.2 版本的基因組，常染色體長度約為2 450 462.292 kb）。

2 結果與分析

2.1 里岔黑豬遺傳多樣性分析 為了探究現存里岔黑豬群體的遺傳結構，將測得的里岔黑豬芯片數據進行遺傳多樣性分析，發現膠州市茂華里岔黑豬良種場中里岔黑豬的有效群體含量為5.88，多態性遺傳標記為0.825，期望雜合度為0.346 8，觀察雜合度為0.334 9；里岔黑豬群體的觀察雜合度小于期望雜合度，可能是群體內近交造成的。

2.2 種群親緣關系分析 本實驗將109 個樣本的39 644個SNP 進行主成分分析。由圖1 可知，膠州市茂華里岔黑豬良種場所有能繁豬之間親緣關系適中，多樣性相對較高。但在第一主成分中，右側的個體之間親緣關系比較緊密，表明親緣關系較近，可能存在近交。進而構建了該群體的G 矩陣（圖2），也驗證了主成分分析的結果，大部分樣本處于中等親緣關系范圍，部分樣本間親緣關系較近，表明存在一定的近親繁殖趨勢，配種計劃中應避免該樣本相遇。

圖1 里岔黑豬主成分分析結果

圖2 里岔黑豬基于G 矩陣的基因組親緣關系分析

2.3 里岔黑豬群體遺傳結構分析 本實驗所用里岔黑豬的血統信息均是由公豬得來，將13 頭公豬結合G 矩陣和聚類分析。結果顯示，現存公豬被分為8 個家系（圖3），自上而下每個家系分別有1、1、1、2、1、2、2、3 頭公豬，具體家系信息見表1。結合群體中的母豬和公豬的關系，除了8 個家系外，另有28 頭母豬的血緣關系與任一公豬的親緣關系均低，因此將其單獨分類。

表1 里岔黑豬家系構建

圖3 里岔黑豬能繁公豬樣本聚類分析

2.4 近交程度分析 在109 頭里岔黑豬中，共檢測到4 250個ROH 片段，個體ROH 的總長度為0～1 008.33 Mb，平均ROH長度為343.56±18.251Mb，ROH長度在300～450Mb之間的個體數目最多，占39.45%（圖5）。里岔黑豬個體含有的ROH 數量為0～70 個，平均為38.99±1.461 個，其中40～50 個ROH 的個體最多，占34.86%（圖6）。通過對群體中各個體的ROH 統計，得到每個個體基于ROH 的近交系數值，由圖7 可以得出，膠州市茂華里岔黑豬良種場近交系數主要集中在0.1～0.2 之間，當前整個保種群體的平均近交系數值為0.14。

圖5 里岔黑豬ROH 長度的樣本數

圖6 里岔黑豬ROH 數量的樣本數

圖7 ROH 的近交系數分布情況

3 討 論

3.1 遺傳多樣性分析 里岔黑豬作為山東省重要的地方豬種之一，其系譜記錄一直采用傳統的人工記錄而后統計匯總的方法，記錄過程中可能會出現錯寫或者漏寫等不確定因素，導致系譜錯誤或者無系譜的結果?；赟NP 芯片數據對里岔黑豬現存群體進行整體評估，將芯片數據構建的家系與記錄系譜相對照，能夠準確了解里岔黑豬現階段的遺傳多樣性和保種現狀。隨著芯片技術的開發，Geneseek[16]、中芯一號、液相50 K[17]等的出現，降低了后續的育種成本，能夠促進我國育種發展進度，從分子層面對其性狀進行分析并指導育種。本研究所用里岔黑豬群體處于世代封閉模式，由于后期系譜丟失導致家系不明確，因此均作為同一世代進行測序。使用“中芯一號”芯片測得里岔黑豬有效群體含量為5.88，多態性遺傳位點為0.825，這可能是群體內近交導致的結果，又由于里岔黑豬群體的期望雜合度小于觀察雜合度，這進一步佐證了群體內近交的現狀。本次測序結果較已報道的其他中國地方豬種（金華豬等）的有效群體含量略低[18]，且明顯低于商品豬有效群體大小[19]；劉彬等[20]利用“中芯一號”芯片對青峪豬群體進行分析，其中有效群體含量最高的為F1世代有12 頭，最低為F3世代僅有3 頭；F1世代多態性遺傳位點最高為0.659 2，F3世代最低為0.536，發現3 個重疊世代中多態性遺傳位點與遺傳多樣性豐富度隨世代均有下降趨勢，這符合中國地方豬種現狀。里岔黑豬雜合度的值低于之前報道的中國其他豬的雜合度[21]，這意味著遺傳變異性相較于以前略有降低。根據基因芯片結果分析可得，可以在維持膠州市茂華里岔黑豬良種場現有家系不減少的情況下，制定合理的配種方案，減少近親交配，維持其遺傳多樣性。

3.2 親緣關系與遺傳結構分析 為了分析茂華里岔黑豬良種場現存群體的遺傳關系，通過NJ tree 并進行聚類分析[22]。整個里岔黑豬群體共109 頭，可以分為8 個公豬的家系和1 個非公豬家系，共有3 個父系起源。目前里岔黑豬群體家系結構可以得出，不同家系公母豬數量和近交系數均不平衡，其中有4 個家系公豬數量最少僅有1 頭，這與蔡春波等[23]利用Illumina CAUPorince 50 K 芯片測得的結果相似，馬身豬被分成3 個家系，共20 頭種公豬，其中部分母豬同時被分到了幾個家系中并出現交叉現象，各家系公母豬數差異較大，家系結構不平衡；王晨等[24]采用Illumina CAUPorince 50 K芯片也得出了相似結果，139 頭藍塘豬被劃分為4 個公豬家系和1 個非公豬家系，平均每個家系4 頭公豬，但最少的家系僅2 頭公豬，家系結構不平衡。這可能與人工選種時僅按照體況健康、符合地方品種特征、性欲是否旺盛等留種標準有關，這也符合前期技術不發達時的選種現狀。在后續保種工作中，各家系需要多加留種，使里岔黑豬群體中各家系血脈維持在一個相對穩定的數量，以防血脈丟失。通過SNP 芯片數據檢測出的各公豬家系，可以通過冷凍精液的方式進行家系保護，有效降低保種支出。

3.3 近交程度分析 ROH 通常被用來估計基因組的近交系數[25]，ROH 長度和ROH 覆蓋的基因組比例可以很好地反映近交系的年齡和來源，從而反映近交的水平。ROH 的長度與個體間親緣成正比，ROH 片段越長其近交可能性越大，反之亦然[26]；ROH 片段數目越多一般認為該群體存在一定的近親繁殖[27]。在本實驗檢測的109 頭里岔黑豬群體中，共檢測到4 250 個ROH 片段，ROH 長度在500 Mb 以下占84.4%，最短ROH 為0 Mb，少部分個體ROH 長度達到900 Mb 以上，里岔黑豬個體含有的ROH 數量在40～50 個ROH 的個體最多，最少的僅有1 個ROH，平均近交系數值為0.14。相較于Liu 等[28]對梁山豬進行全基因組分析，得出涼山豬的大多數（80%）ROH 在100 Mb 以內，ROH 平均長度為63.24 Mb，群體的平均近交系數為0.026，里岔黑豬ROH 長度長于梁山豬，這可能是里岔黑豬在長期交配過程中，無意識的近交導致的結果，從兩者的平均近交系數也能反映出近交的情況；Xu 等[29]通過基因組還原和測序技術（GGRS）進行基因分型，202 頭金華豬的ROH 數量大多數（77%）由1～5 Mb 的片段構成，ROH 長度小于10 Mb 的占總數的92%。古代與近代祖先對金華豬群體基因組產生了一定的影響，盡管群體內大部分個體近交系數較低，但仍有較高水平的近交系數，在后續留種過程中逐步淘汰該種個體，保證其遺傳多樣性；Wang 等[30]利用Illumina 60 K 芯片將江西省的3個地方品種豬在全球豬群中進行比較，東鄉花豬的平均ROH 最高達到140.5 Mb、萍鄉兩頭烏平均ROH 最高達到38.7 Mb，稍高于野豬（20.6 Mb），玉山黑豬平均ROH 稍高為121.7 Mb。里岔黑豬相比較野豬和中國部分地方豬種（梁山豬、金華豬等）近交系數相對較大，可能是因為從小群體發展至如今規模，一直以本交作為配種方式，所以存在一定劣勢，后期可以通過SNP 數據合理安排配種計劃，必要時可以改變配種方式，采用人工授精[31]等方式保證群體結構的穩定。

3.4 制定適合的育種方案 通過分析里岔黑豬NJ 數據，從8 個家系中各選1 頭性狀優良的種公豬，將96 頭能繁母豬隨機分為8 組，每組12 頭。參考Windig 等[32]通過計算機模擬的幾種保種方案，按照輪換交配的方式，制定適合于膠州市茂華里岔黑豬良種場的育種方案：第1 輪按照家系A 與第1 組母豬交配；家系B 與第2 組母豬交配；家系C 與第3 組母豬交配……家系H 與第8 組母豬交配。第2 輪按照家系A 與第2 組母豬交配；家系B 與第3 組母豬交配；家系C 與第4 組母豬交配……家系H 與第1 組母豬交配。第3 輪家系A 與第3 組母豬交配；家系B 與第4 組母豬交配；家系C 與第5 組母豬交配……家系H 與第2 組母豬交配，按照步驟依次輪轉交配，能有效減緩近交系數增長速度。

4 結 論

本研究從分子水平基于里岔黑豬群體進行分析，發現里岔黑豬群體遺傳多樣性較豐富；各家系中個體數相對較少且差異明顯，需注意家系內留種；為防止膠州市茂華里岔黑豬良種場群體內近交系數增長過快，需制定合理配種方案，必要時可以從其他保種場引入新鮮血液。