基于全基因組重測序SNP標記的148份馬鈴薯種質遺傳多樣性分析

韓志剛，郝文勝，謝 銳，郭景山，伊六喜，侯建華*

(1 內蒙古農業大學 農學院，呼和浩特 010019； 2 內蒙古自治區農牧業科學院，呼和浩特 010031)

中國的馬鈴薯研究始于20世紀30年代，在40年代分批次從美國、日本等國家引種，50年代利用這些引進的資源開始了有規模的自主育種工作[1]。經歷半個多世紀，中國的馬鈴薯育種工作從無到有，從傳統育種到現在的與分子標記輔助育種相結合的發展經歷，取得了長足的進步。截止2020年利用前人引進的資源中國的馬鈴薯育種家們選育審定登記的新品種873份，而這些自主選育的新品種在中國的不同生態地區發揮著各自的優勢，為當地老百姓帶來了實惠，同時，有些國外種質資源至今仍在大面積種植，例如，美國的‘布爾班克’、‘大西洋’，加拿大的‘夏波蒂’，荷蘭的‘費烏瑞它’等品種?！最^翁’、‘卡它丁’等老品種更是育種家們沿用已久的經典親本[2]。因此，對國外引進的馬鈴薯種質資源及國內選育新品種的充分利用就顯得十分必要，而要有效合理地利用這些資源，就必須了解其遺傳背景、親緣關系，對其遺傳多樣性進行研究[3]。分子標記無疑是對各種生物群體的遺傳多樣性研究中最常用的方法。邸宏[4]利用AFLP和RAPD分子標記對國內主栽品種、地方品種、新型栽培種及抗PVY品種等4類資源種植共157份品種進行了遺傳多樣性分析，表明新型栽培種和地方品種的遺傳基礎豐富于普通栽培種，普通栽培種與新型栽培種在遺傳組成上有明顯區別，而地方品種與普通栽培種無明顯差異。徐敏[1]用SSR分子標記對國內審定的203份品種遺傳多樣性分析發現，中國馬鈴薯審定品種遺傳多樣性較低。程永芳[5]基于SRAP標記通過UPGMA法可將54份種質資源分為6大類群。單核苷酸多態性(single nucleotide polymorphism，SNP)指相同或不同的物種個體間的基因組DNA序列同一區間上的單個核苷酸存在差別的現象，是由美國學者Lander提出的DNA遺傳分子標記，起初應用于醫學[6]。隨著測序技術的不斷發展SNP標記以廣泛應用于各類植物的遺傳學研究，而馬鈴薯SNP標記的研究還處于初級階段[7-10]。王艦[11]基于簡化基因組SNP標記對288份種質資源的遺傳多樣性研究發現，馬鈴薯種質資源地域間差異不明顯。Sofía等[12]用SNP標記對231個不同基因型馬鈴薯品種做了群體結構分析和主成分分析。Li等[13]用 6 487 006 個高質量 SNPs 構建了有根鄰接系統發育樹，以S.chomatophilumBitter 作為外群，對‘clade 4’成員和栽培種進行了主成分分析，這兩個分析得到類似的結果，將野生馬鈴薯分成兩個亞群，其中一個亞群的物種位于秘魯，為北方野生亞組；而另一亞群的物種位于阿根廷，玻利維亞和智利，為南方野生亞組，這些分析得到了一致的結果，相對于南方野生亞組，栽培馬鈴薯與北方野生群體的親緣關系更近。以上研究報道表明利用分子標記技術分析馬鈴薯不同品種間的遺傳多樣性是可行的。

本研究以內蒙古農牧業科學院馬鈴薯課題組多年收集的國內外500余份馬鈴薯種質資源為基礎，篩選出性狀穩定、表現優良、系譜明確的148份核心種質資源，通過全基因組重測序，以全基因組范圍內的SNP標記，全面挖掘遺傳變異信息，進一步了解148份種質資源的親緣關系和遺傳多樣性，旨在提高馬鈴薯種質資源的利用效率，為馬鈴薯定向育種提供依據。

1 材料和方法

1.1 供試材料

材料為四倍體馬鈴薯普通栽培種148個(表1)，國內品種108個，其中華北地區48個、西北地區22個、東北地區16個、華中地區4個、西南地區15個以及華東地區3個品種；國外品種40個，包括歐洲22個、美洲17個及亞洲日本的1個品種。這些材料為系譜明確且田間表型穩定的馬鈴薯優良品種。

表1 基于SNP分析的四倍體馬鈴薯品種名稱來源及系譜

續表1 Continued Table 1

續表1 Continued Table 1

續表1 Continued Table 1

1.2 方 法

1.2.1 試驗材料的種植各試驗材料種植于呼和浩特市武川縣國家馬鈴薯產業技術體系呼和浩特試驗站試驗田內，各品種種植5株，行長1.5 m，株距0.3 m，行距0.9 m，機器開溝人工點種覆土，前作葵花。

1.2.2 樣品DNA的提取供試材料148個種植出苗后，選擇健壯無退化表現的植株取新鮮嫩葉2.0 g，參照Stewart等[14]提出的CTAB法提取基因組DNA，用NanoDrop核酸蛋白測定儀檢測DNA的純度，Qubit對DNA濃度進行精確定量。

1.2.3 SNP標記檢測148個樣品基因組DNA送至諾禾致源公司進行高通量重測序。高通量IlluminaHiSeqTM測序得到的原始圖像數據經轉化為序列數據，去除未測出的堿基等干擾信息后得到有效數據。通過BWA軟件，將這些有效數據與已公布的馬鈴薯參考基因組(https://www.ncbi.nlm.nih.gov/genome/?term=Solanum+tuberosum+L.)進行比對。SAMTOOLS軟件進行群體SNP的檢測，貝葉斯模型檢測群體中的多態性位點，通過有效過濾(測序深度>5X；最小等位基因頻率MAF>0.1；Call rate>90%)，獲得高質量的SNPs用于遺傳多樣性分析。

1.2.4 群體結構分析和聚類分析本研究以群體結構分析和聚類分析2種方法來探究148個馬鈴薯種質資源各品種間群體結構和親緣關系，兩者的結果可以相互驗證。運用TreeBestv1.9.2軟件計算距離矩陣，以此為基礎換算成遺傳相似系數，通過鄰接法(neighbor-joining method)構建群體聚類圖，引導值(bootstrap values)經過達1 000次計算獲得。采用PLINK進行群體結構分析，首先創建PLINK的輸入文件—Ped文件，然后利用admixture軟件構建群體遺傳結構和群體世系信息。PLINK (http://pngu.mgh. harvard.edu/-purcell/plink/)。

1.2.5 部分骨干親本雜交后實生種子統計2018和2019年間按照育種目標對148份種質資源內篩選出的部分骨干親本于馬鈴薯盛花期進行雜交授粉，待坐果成熟后對漿果內實生種子(True Potato Seed TPS)有無、數量等進行統計分析，并結合各品種間遺傳相似系數即親緣關系遠近，為不同親本間雜交組合選配提供參考。

2 結果與分析

2.1 DNA樣本檢測與重測序質量評估

148份馬鈴薯種質資源提取到的DNA濃度在20.6～166.0 μg/μL之間，平均為92.88 μg/μL，總量在2.47～17.42 μg范圍內，平均為7.44 μg，表明DNA質量符合后續重測序的要求。經每樣本重測序數據與參考基因組比對Q20值(錯誤率1%以下)在95.77%～98.04%，平均為96.92%，錯誤率在Q30值(錯誤率0.1%以下)在89.27%～94.16%之間，平均值為91.77%，GC含量在35.44%～36.34%之間，平均為35.79%。與參考基因組的比對率最高96.51%，最低93.53%，平均95.66%，測序深度最高21.91X，最低11.93X，平均深度為15.12X，平均覆蓋度為95.55%。以上數據表明各樣本與參考基因組的相似度達到重測序標準，并且足以滿足該群體的遺傳多樣性分析。

2.2 多態性SNP位點的分布

經重測序148個馬鈴薯樣本共獲得2 315 860個SNP位點，過濾篩選后獲得1 209 969個高質量SNP位點，其中被明確定位在染色體水平上的SNP位點為1 192 472個，比例為98.55%，這些SNP位點將用于148份種質資源的遺傳多樣性分析。

SNP位點從染色體分布情況來看差距較大，11號染色體沒有SNP位點，而最多的是5號染色體有258 909個位點占總數的21.71%，其次為1號染色體的194531個位點占比為16.08%，之后依次為2號、10號、6號、4號、9號、8號以及3號和12號染色體，SNP位點數量分別為153 224、149 793、133 432、123 134、58 965、54 202、30 106、27 671個，占比依次為12.17%、12.38%、11.03%、10.17%、4.87%、4.48%、2.49%和2.29%，除了11號染色體外的7號染色體SNP分布數量最少為8 485個，僅占0.70%(圖1)。

圖1 SNP位點在染色體上的分布Fig.1 Distribution of polymorphic SNPs on chromosomes

2.3 各品種間的遺傳相似系數

利用國內外148個馬鈴薯品種全基因組重測序獲得的SNP標記，通過TreeBestv1.9.2軟件計算出各品種間遺傳相似系數，共獲得10 878個遺傳相似值，范圍為0.784～0.958，平均值為0.842。遺傳相似系數最低的是華北地區山西省農科院高寒區作物所育成的‘系薯1號’和同為華北地區內蒙古農業大學育成的‘內農薯1號’之間的0.784，遺傳相似系數最高的是華北地區河北省張家口農科院育成的2個品種‘冀張薯-59’和‘冀張薯20號’之間的0.958。各品種間遺傳相似系數主要集中于0.800～0.880之間，有10 604個，占97.5%，小于0.800的有55個，占0.5%，大于0.88的有219個，占2.0%(圖2)。

圖2 品種間遺傳相似系數的次數分布Fig.2 Distribution of genetic similarity between varieties

按照各品種育成地區將中國馬鈴薯品種劃分為華北、西北、東北、華中華東、西南等5個區域，國外品種劃分為歐洲、美洲、亞洲(日本)等3個區域。對各區域內育成品種間遺傳相似性進行分析發現(表2)，中國華北地區平均遺傳相似系數為0.841，西北地區為平均系數為0.843，東北地區為0.842、華中華東地區是0.845，西南地區各品種間平均遺傳相似系數最高為0.848，國外歐洲品種間平均遺傳相似系數為0.844，美洲品種間平均遺傳相似系數為0.850。數據結果表明無論是國內還是國外各區域品種間平均遺傳相似系數差距不明顯。

表2 不同地區育成品種間遺傳相似系數

2.4 群體結構分析

將148個馬鈴薯品種通過PLINK和admixture軟件構建群體遺傳結構和群體世系信息。進行群體結構分析時將分群數K值設置為2-8。根據各品種間的交叉驗證錯誤率ΔK最小時K=6，說明最佳分群K值為6(圖3)，因此148個品種分為6個群最為接近于群體的真實情況(圖4)，圖中每一個直方柱子代表1個品種，每根柱子的顏色代表遺傳組分的構成。利用admixture軟件對每個個體的遺傳背景進行分組并計算相應的Q值，即某一材料基因組變異源于某個群體的概率，Q值>0.6視為來源較單一、遺傳背景較純，Q值<0.6視為來源混合，遺傳背景較為雜[15-16]。

圖3 交叉驗證錯誤率ΔK值隨K值的變化Fig.3 The variation of cross validation error rate ΔK with K value

淺藍色代表群組1；藍色代表群組2；淺綠色代表群組3；綠色代表群組4；粉色代表群組5；紅色代表群組6圖4 148份四倍體馬鈴薯種質SNPs群體結構分析(K=6)Light blue blocks represent Pop1; Blue blocks represent Pop2; Light green blocks represent Pop3; Green blocks represent Pop4; Pink blocks represent Pop5; Red blocks represent Pop6Fig.4 Population structure analysis of 148 tetraploid potato varieties (K=6) by using identified SNPs

本次分析的6個群組主色段分別為淺藍、藍、淺綠、綠、粉以及紅色。第一群組(Pop1)包括21個品種，主色段為淺藍色，其中5個品種的遺傳背景較純，分別為中國農科院蔬菜花卉研究所的‘中薯19號’、‘中蔬9518-15’，內蒙古烏蘭察布市農科院的‘烏蘭一號’以及國外美洲地區加拿大的‘亨特’、美國的‘Atlantic’，其他品種含有多種遺傳組分，屬于混合型。第二群組(Pop2)藍色為主色段，包括30個品種，其中6個品種的遺傳背景較純，分別是華北區中國農科院蔬菜花卉研究所的‘中薯10號’、‘中薯4號’及‘中薯21號’，河北省張家口農科院的‘壩薯9號’，西南地區云南省農科院經濟作物所的‘云薯506’，美洲地區加拿大的‘阿旺’。其他品種均為混合型。第三群組(Pop3)有33個品種，主色段是淺綠，其中華北區山西省高寒區作物所的‘晉薯24號’、內蒙古農牧業科學院的‘紅美’、河北省張家口市農業科學院的‘冀張薯5號’、中國農科院蔬菜花卉研究所的‘中蔬613’，西南地區貴州省威寧地區農科所的‘威176選2’以及歐洲地區俄羅斯的‘Epicure’和‘Cибиpяк’等7個品種的遺傳背景較純，其他品種為混合型。第四群組(Pop4)共18個品種，主色段為綠色，其中包括華北地區的‘中薯12號’、‘冀張薯-59’、‘冀張薯20號’，東北地區黑龍江省農科院克山分院的‘克新1號’以及美洲地區加拿大的‘貝勒斯利’和美國的‘Kennebec’等6個品種的遺傳背景較純。第五群組(Pop5)包括29個品種，主色段是粉色，其中6個品種遺傳背景較純分別是華北地區的‘冀張薯12號’和‘系薯1號’，西北地區甘肅省農業科學院馬鈴薯研究所的‘隴薯3號’，華東地區山東省希森集團的‘希森3號’以及歐洲地區荷蘭的‘荷蘭薯’和‘Favorita’，其余品種為混合型。第六群組(Pop6)包括17個品種，主色段為紅色，其中6個品種較遺傳背景較純，分別為西北地區甘肅省農業科學院馬鈴薯研究所的‘渭會四號’，西南地區四川省農科院作物所的‘川芋39號’，東北地區吉林省蔬菜花卉科學研究所的‘春薯3號’，歐洲地區德國的‘強壯’、‘北斗星’,以及美洲地區美國的‘Katahdin’，其他品種為混合型。

2.5 聚類分析

為探明148個馬鈴薯品種間的親緣關系，本研究以遺傳相似系數為基礎，通過鄰接法構建NJ聚類圖。結果表明(圖5)，148個馬鈴薯品種被劃分為3個類群。

圖5 基于SNP標記的148份四倍體馬鈴薯種質的進化樹Fig.5 Neighbor-joining tree for 148 tetraploid potato varieties based on SNP markers

第Ⅰ類有6個品種，也是3個類群中品種數最少的一個群，遺傳相似系數在0.852～0.956之間。包含華北地區的3個品種，西北及東北地區各1個品種，國外美洲的1個品種?！墒?8號’的父本‘3-4-4’和‘烏盟851’的母本‘474-7-9’有野生種的親緣，據此推測僅這6個品種在148個群體中形成獨立分支很可能與野生種親緣有關。

第Ⅱ類有80個品種，分為3個亞類，遺傳相似系數在0.784～0.956之間，平均為0.843，其中華北地區內蒙古華頌公司育成的‘華頌7號’與歐洲德國的‘Agria’遺傳相似系數最高為0.956，華北區內蒙古農業大學育成的‘內農薯1號’與山西省高寒區作物所育成的‘系薯1號’相似系數最低為0.784。Ⅱ-1亞類有18個品種，主要包括華北地區5個品種，西北地區5個品種，國外美洲的5個品種，其中西北地區甘肅農科院育成的‘隴薯系列’5個品種緊相鄰，美洲5個品種中除了‘亨特’為加拿大品種外，其他4個均為美國品種。從系譜上看，此亞類18個品種親本大多為國外引進品種，尤其是‘Shepody’、‘疫不加’、‘燕子’、‘菲多利’等都是國外的老品種?！甃enape’是‘Atlantic’的母本也是‘Belchip’的父本，因此‘Lenape’與‘Atlantic’相似系數較高為0.895，與‘Belchip’的相似系數為0.889?！畠绒r薯1號’與‘隴薯8號’的母本同為‘Atlantic’，兩者相似系數較高為0.896，‘隴薯11號’和‘隴薯14號’母本為‘L9712-2’，兩者遺傳相似系數也為較高的0.881。Ⅱ-2亞類共有30個品種，此亞類較為典型的是中國農科院蔬菜花卉研究所育成的‘中薯系列’7個品種聚在此亞群，其中‘中薯19’號與‘中薯21’號父母本相同，為同胞關系，兩者相似系數為0.899?！惺?9號’、‘中薯21號’、‘中薯18號’的父本同為‘C93.154’，此類半同胞關系的遺傳相似值也較高，‘中薯18號’與‘中薯19號’之間為0.854，與‘中薯21號’之間為0.876。‘中薯20號’的父本與‘中薯19號’和‘中薯21號’的母本同為品系‘92.187’，因此‘中薯20號’與‘中薯19號’的遺傳相似系數為0.858，與‘中薯21號’的相似值為0.868。歐洲德國品種‘Agria’為荷蘭品種‘Ramos’的母本，兩者遺傳相似系數為0.890。Ⅱ-3亞類有32個品種，為第Ⅱ類群中品種數最多的亞類，華北區品種有14個，其他品種無明顯的地域性，育成地區比較混雜。從系譜角度觀察，‘冀張薯14號’和‘冀張薯-88’為同胞關系，具有相同的父母本，兩者遺傳相似值是此亞群最高的0.955，‘隴薯12號’和‘隴薯13號’的父本同為‘L0202-2’屬半同胞關系相似值為0.942，‘冀張薯8’號是‘冀張薯24號’的母本，其相似值是0.905。

第Ⅲ類共有62個品種，分為2個亞類，遺傳相似值在0.789～0.958之間，平均值是0.847，其中華北區中國農科院蔬菜花卉研究所育成的‘中蔬613’和西北區甘肅農科院育成的‘隴薯5號’之間遺傳相似系數最低為0.789，最高的是‘冀張薯-59’和‘冀張薯20號’之間的0.958，也是148群體各品種間遺傳相似值最高的2個品種。Ⅲ-1亞類有22個品種，從聚類圖上看華北區中國農科院蔬菜花卉研究所育成的 ‘中薯系列’3個品種相鄰，河北張家口農科院育成的‘冀張薯-59’和‘冀張薯20’號相鄰，其他品種無明顯育成地區的一致性。由各品種系譜關系看，‘冀張薯-59’和‘冀張薯20號’為父母本相同的同胞關系，‘中薯3號’是‘中薯9號’的父本，兩者之間的相似值為0.949。‘涼薯3號’是‘川芋10號’的父本，相似值為0.856?！惺?號’和‘XS-06YH-03’的母本同為‘Shepody’，兩者的相似值為0.873。Ⅲ-2亞類共有40個品種，包括華北地區9個品種、西北地區6個品種、東北地區5個品種、西南地區5個品種、華中地區1個品種、歐洲8個品種及美洲的6個品種，品種育成區域較為混雜。東北地區5個品種均為黑龍江省農業科學院克山分院育成的‘克新系列’品種，其中‘克新23號’與‘克新25號’為同胞關系，在聚類圖上緊相鄰，遺傳相似值為0.909。‘Desiree’是‘青薯168’的父本，兩者聚類位置也是緊相鄰，相似值為0.884，另外‘Shepody’是‘Innorvator’的母本，兩者相似值為0.886，其他品種在系譜上無明顯親緣關系。

2.6 TPS情況統計分析

2018至2019年間按照“鮮食、高產(3 000 kg以上)”育種目標，篩選出骨干親本38份，配制349個雜交組合，坐果并有TPS的組合(視為成功組合)32個，占9.17%。雜交成功組合中‘中薯12號’(♀)與‘早大白’(♂)之間的遺傳相似系數最低為0.807，‘費烏瑞它’(♀)與‘冀張薯12號’(♂)的遺傳相似系數最高為0.895，32個組合平均為0.843(表3)。

表3 2018和2019年部分骨干親本雜交選配TPS統計

3 討 論

種植資源是作物基礎研究和遺傳改良的基礎物質，其數量和質量是定向育種的兩大制約因素[17]。而種質資源的遺傳多樣性是育種工作中各品種間的雜交選配、親緣關系遠近、雜種優勢強弱的重要指標，因此，對種質資源的遺傳背景，各品種間的遺傳相似系數的了解顯得格外重要[18]。隨著測序技術的發展從起初的第一代測序到現在第三代高通量測序使得各作物的起源、馴化、傳播，各群體間的遺傳多樣性研究更加深入明了。相對于簡化基因組測序，全基因組重測序理論上應該覆蓋整個基因組，測序深度10X以上可覆蓋整個基因組的90%～99%，足以達到群體遺傳多樣性分析。分子標記從最初的限制性片段長度多態性(RFLP)經歷了隨機擴增多態性DNA(RAPD)及簡單序列重復(SSR)直到現在的SNP，可謂層出不窮[19]。SNP作為第三代分子標記是基于第二代SSR等分子標記發展起來的新型標記技術，近年來已廣泛應用于水稻、玉米、小麥、棉花等各類作物的品種鑒定、遺傳圖譜構建以及全基因組關聯分析等研究工作[20-22]。相對于其他分子標記，SNP在基因組上做到了單個核苷酸水平，且具有分布廣泛、數量多、遺傳穩定性強等優點；其缺點是做到全基因組水平的SNP標記的話測序成本較高，易發生假陽性、錯誤率高等[23]。因此，為了確保獲得高質量的SNPs一定要嚴格把控測序深度、最小等位基因頻率、比對率等過濾標準。本研究中，為了獲得高質量的SNPs，148份樣品的平均測序深度為15.12、Q30值為91.77%、最小等位基因頻率高于10%，正是因為高標準的過濾參數使得低質量的SNPs被去除掉，這也許就是導致11號染色體未能檢出高質量SNP的原因。

群體結構分析主要是能夠明確該群體的分群合理性及每個個體遺傳組分的多少。本研究中，148份馬鈴薯品種被劃分為6個群組，Q值>0.6的品種有36個，占比24.3%，說明遺傳背景純、來源較單一的品種占比不大，大部分品種遺傳組分混雜程度較高，基因交流較多。這36個品種中來源為華北地區品種有16個占44.4%，國外品種有13個占36.1%，二者共占80.5%，說明華北區育成品種及國外品種的遺傳背景純正比例相對高于其他區域品種。這可能是由于中國馬鈴薯北方一季作區氣溫較低、降雨量少、無霜期短適宜生長的品種較少，供育種者們選擇的骨干親本較少有關。聚類分析呈現的是群體樣本DNA序列間的差異，也就是各品種間的親緣遠近關系。148份種質資源被劃分為3個群，遺傳相似系數0.784～0.958之間，平均0.842。按照品種來源即育成單位地理位置劃分為8個區域，各區域相似值在0.841～0.850之間，區域間平均相似系數差距不明顯。從聚類圖上看，西北地區‘隴薯系列’的7個品種集中在Ⅱ-1亞群，華北地區的‘中薯系列’6個品種則集中在Ⅱ-2亞群，‘冀張薯系列’的5個品種集中于Ⅱ-3亞群，而‘中薯系列’的3個品種和‘冀張薯系列’的2個品種相鄰于Ⅲ-1亞群，由此看出北方一季作區育成的部分品種聚在一起有著一定的地理區域性與段紹光等[24]利用SSR標記分析馬鈴薯品種遺傳多樣性結果相一致。華北區品種在聚類圖上顯示相鄰的結果與群體結構分析遺傳組分單一，遺傳背景狹窄相互驗證。其他區域品種總體上看來并無明顯的地理區域性，不同地理來源的馬鈴薯品種相互交錯分布，說明品種來源的地理差別與親緣關系并無必然聯系，此結果與王艦[11]分析結果一致。究其原因：首先，各育種單位相互引種頻繁，在長期的選育過程中存在基因交流的情況；其次，不同區域引種品種為適應當地環境，逐漸減少了遺傳差異所致。

群體結構分析與聚類分析對比發現，群體結構分析的Pop1(21個品種)與聚類分析的Ⅱ-1亞群(18個品種)有12個品種相同，Pop2(30個品種)和Ⅱ-2亞群(30個品種)有13個品種相同，Pop3(33個品種)與Ⅱ-3亞群(32個品種)的相同品種有13個，Pop4(18個品種)的8個品種聚在Ⅲ-1亞群(22個品種)，Pop5(29個品種)與Pop6(17個品種)共46個品種和Ⅲ-2亞群(40個品種)與Ⅰ群(6個品種)的46個品種相對應有15個品種相同，以上結果說明群體結構分析與聚類分析分群結果基本吻合，也相互驗證了148份種質資源分群的合理性。

系譜分析與全基因組SNP標記的聚類分析相結合發現，同胞關系的品種基本聚在一起。例如：‘冀張薯14號’、‘冀張薯20號’、‘冀張薯-59’及‘冀張薯-88’為父母本完全相同的同胞關系，聚類結果顯示‘冀張薯14號’和‘冀張薯-88’相鄰于Ⅱ-3亞群，‘冀張薯20號’和‘冀張薯-59’相鄰在Ⅲ-1亞群?！诵?3號’與‘克新25號’也是同胞關系，二者相鄰在Ⅲ-2亞群，‘中薯19號’和‘中薯21號’聚在Ⅱ-2亞群。半同胞關系的品種也基本聚在同一亞群中，在聚類分析結果中已舉例說明。由此可見，本次篩選到的SNPs可以揭示馬鈴薯品種資源的遺傳差異以及品種間的親緣關系。當然，系譜信息與聚類結果不一致的現象是存在的，例如：‘燕子’是‘烏蘭一號’的母本，‘Favorita’是‘蒙薯18號’的母本，‘Innovator’和‘XS-06YH-03’母本相同，‘隴薯8號’與‘冀張薯12號’母本相同，系譜信息顯示以上品種間確實存在親緣關系，但并未聚在同一個亞群。出現此類現象可能與以下幾點有關：一是馬鈴薯普通栽培種為同源四倍體，有48條染色體，形成配子多樣性高，遺傳分離較大；二是育種者往往注重高產、抗病等性狀，長期的定向育種造成某個性狀位點過于集中，從而造成其他性狀位點逐漸丟失；三是往往育種單位種質資源品種較多，在引種、繁種、保種時難免會出現系譜或品種名稱記載偏差的可能，導致同物異名、同名異物的結果。

馬鈴薯有別于其他主糧作物的最大特點是無性繁殖。優點在于只要獲得綜合性狀優良的個體后，通過無性繁殖的方式可以將其穩定下來，從而獲得穩定的品種品系；但是，由于無性繁殖的特征也使得現有馬鈴薯普通栽培種處于遺傳雜合狀態，而這對于新品種選育來說是不利的[25]。馬鈴薯普通栽培種均為同源四倍體材料，四體遺傳和高度雜合的特點提高了性狀分離的復雜程度，使得馬鈴薯雜交組合必須要有相當大的群體，從而提高優選率。育種工作中親本的選配對于是否能得到優良后代是至關重要的。父母本的遺傳背景差異大、親緣關系較遠的組合其雜種優勢更強，變異類型更多，優選幾率也會更高一些。本研究中通過2年部分骨干親本的雜交情況看出，馬鈴薯普通栽培種的雜交成功率很低，349個雜交組合中僅有32個組合雜交成功，占比僅為9.17%，從遺傳學角度考慮這與近緣雜交形成的雜種優勢不強、變異類型不多、遺傳背景狹窄等因素有著密切的關系。32個成功組合的平均父母本遺傳相似系數為0.843，因此推斷制定雜交組合選配時可以考慮0.843為父母本的遺傳相似參考值，育種目標相同情況下，應該考慮遺傳相似值在0.843以下的2個品種作為優選，以便提高雜交成功率。

本研究以國內外148份四倍體馬鈴薯普通栽培種為群體，通過第三代高通量測序手段進行全基因組重測序，以SNP為分子標記，遺傳相似系數為指標，利用群體結構分析與聚類分析相結合的方法，參照系譜信息，初步明確了各品種間的遺傳背景和親緣關系及群體的遺傳多樣性。旨在現有條件下盡可能避免親緣關系過近品種間雜交，有助于豐富育成品種的遺傳多樣性，為馬鈴薯普通栽培種育種工作的親本選配提供參考。本研究相對于同類研究，首次采用全基因組重測序技術，獲得的高質量SNP數量較大、各品種系譜信息較全，為分析結果的可靠性提供了保障。

研究發現，以SNP為分子標記的遺傳多樣性分析結果顯示，部分馬鈴薯品種間親緣關系具有一定的地理區域性，但并非必然關系。148份種質資源中大部分品種間遺傳相似系數較高，遺傳背景不夠豐富，雜交成功率低，因此，在育種中亟需引入新型種質資源，繼而為種質創新提供更豐富的遺傳基礎。