基于SSR標記的255個棗品種親緣關系和群體遺傳結構分析

劉秀云，李 慧，2，劉志國，趙 錦，劉孟軍，4

（1河北農業大學/中國棗研究中心，河北保定 071001；2寧夏農林科學院種質資源研究所，寧夏銀川 750002；3河北農業大學生命科學學院，河北保定 071000；4北京林果業生態環境功能提升協同創新中心，北京 102206）

基于SSR標記的255個棗品種親緣關系和群體遺傳結構分析

劉秀云1，李 慧1，2，劉志國1，趙 錦3，劉孟軍1，4

（1河北農業大學/中國棗研究中心，河北保定 071001；2寧夏農林科學院種質資源研究所，寧夏銀川 750002；3河北農業大學生命科學學院，河北保定 071000；4北京林果業生態環境功能提升協同創新中心，北京 102206）

【目的】棗原產中國，種質資源豐富。對來自中國22個省區不同用途的255個棗品種進行SSR分析，揭示這些不同產地來源的棗種質資源之間的親緣關系和群體遺傳結構，為棗種質資源的科學管理和分子標記輔助育種提供參考。【方法】利用改良CTAB提取供試棗種質的基因組DNA，以前期棗基因組測序挖掘出的SSR引物為基礎，進行高效率引物篩選，并利用SSR分子標記技術對255份棗種質資源的基因組DNA進行PCR擴增，然后利用8%聚丙烯酰胺凝膠電泳分離，銀染后顯色。根據條帶有無統計數據，計算出多態性位點百分率（PIC），用NTSYS軟件進行UPGMA聚類分析；利用Structure軟件分析群體遺傳結構，計算出最適群體數目，構建遺傳結構圖。【結果】從64對SSR引物中篩選出23對高效率SSR引物，在供試材料中共檢測出117個多態性位點，各引物擴增的多態性位點數為2—10條，每對引物平均擴增多態位點數為5.09個，PIC值變幅為0.359—0.727，平均為0.548，這些多態性引物可應用到其他棗種質資源的研究中；建立了只需1—2個標記就可鑒別出來的部分棗品種的SSR指紋，可用于這些品種的快速分子鑒定；255個棗品種的聚類分析將所有棗品種分為15個亞類，包括4個大類和11個小類，不同品種間的相似系數范圍0.71—1.00，其中北京花生棗單獨聚為一類，與其他棗品種關系較遠；‘奉節雞蛋棗’和‘溆浦雞蛋棗’、‘陜西奶棗’和‘天津大馬牙棗’的相似系數均為 1.00；結合聚類圖、供試品種的用途和原產地分析，不同棗品種間的親緣關系與品種原產地有一定相關性，但和品種用途沒有顯著相關性。群體結構分析中，通過繪制K與ΔK的關系圖，K=15時，ΔK最大，據此將255個棗品種也同樣劃分為15個群體，與聚類分析結果基本一致；進一步分析表明，各群體中大部分品種血緣關系比較單一，較少品種含有其他類群的遺傳成分。總體看，山西或陜西的棗品種出現在絕大部分居群中，說明這兩個省的資源在不同群體間的基因交流中發揮了重要作用；南方栽培區域中來自湖南的棗品種形成了相對獨立的居群，可能是其起源相對單一，且在長期栽培過程中和其他產地棗品種間基因交流較少所致。上述結果表明，供試棗品種中與來源區域相關性明顯的品種由相同地域內棗品種演化而來，而另一部分與來源產地相關性不明顯的品種則是由不同區域間品種經過頻繁的基因交流和重組選育而來，融合了不同區域品種的特點，從而沒有了明顯的區域特征。【結論】不同地理環境在棗品種的群體進化中發揮了較重要的作用，影響了不同產地間棗種質資源的遺傳結構組成。

棗；SSR分子標記；產地；親緣關系；遺傳結構

0　引言

【研究意義】棗（Ziziphus jujuba Mill.）原產中國，栽培歷史悠久，遺傳變異豐富，“同物異名”和“同名異物”現象較為突出。采用傳統的形態學鑒定，費時費力，并且難以準確的鑒定性狀相似的品種，而利用分子標記構建棗品種指紋圖譜并進行遺傳結構分析，有助于棗品種鑒定、遺傳進化分析和分子輔助育種等研究。【前人研究進展】RAPD（random amplified polymorphic DNA）、AFLP（amplified fragment length polymorphism）等分子標記曾被成功用于棗的指紋圖譜構建和親緣關系分析。趙錦等［1］對棗樹品種、品系及其近緣種進行了RAPD分析，喬勇等［2］利用AFLP技術對21個棗品種進行了分析，上述研究中供試材料的聚類結果與傳統觀點基本一致。近年來，SSR（simple sequence repeat）和SNP（single nucleotide polymorphisms）在國際植物新品種權保護聯盟（UPOV）的BMT分子測試指南中被確定為構建DNA指紋數據庫的標準標記方法。由于SSR標記操作簡單，已成為當前各種作物建庫的首選標記［3］。SSR是由1—6個堿基對組成的簡單序列重復，又稱微衛星DNA，已經被應用在葡萄、獼猴桃、杏、蘋果等果樹上進行遺傳多樣性分析［4-7］。2012年，麻麗穎等［8］利用12對SSR引物也構建了 36份棗品種的指紋圖譜。遺傳結構分析研究已在蘋果［9-10］、梨［11］、杏［12-16］、毛櫻桃［17］、楊梅［18］、荔枝［19］、芒果［20］、獼猴桃［21］、核桃［22-23］、石榴［24］、桃［25］中開展。在棗樹研究中，李海濤［26］利用ISSR標記進行研究的結果表明，河南棗主栽品種與灰棗的總群體間分化程度較高，群體間基因交流較少，基因流可防止由遺傳漂變引起的群體間的遺傳分化；殷曉［27］利用ISSR標記對237份棗種質的遺傳多樣性進行了分析，結果表明黃河沿岸棗群體的雜合程度為重度雜合，遺傳多樣性豐富，變異程度較高，趨向遠緣雜交。【本研究切入點】目前尚未有利用SSR技術對大規模棗種質材料進行遺傳結構分析的報道，因此，在前人研究基礎上，利用前期基于棗基因組測序開發的SSR引物［28］，對來自全國22個省區主產地的255份棗種質進行更廣范圍的遺傳多樣性及群體結構研究。【擬解決的關鍵問題】分析中國棗種質資源的起源，探討其遺傳多樣性水平和遺傳結構，為指紋圖譜構建和分子輔助育種提供理論依據。

1　材料與方法

試驗于2015年在河北農業大學進行。

1.1材料

選取來自全國22個省區及主產地的255個棗樹品種，供試品種信息詳見表1。所有品種材料均采自山西太谷國家棗種質資源圃，于2013年6月采集供試材料的新鮮幼嫩葉片，置于冰盒中帶回實驗室，轉入-20℃冰箱貯藏，以備提取基因組DNA。

1.2方法

基因組DNA提取方法按照肖京［29］的方法進行。

SSR引物為筆者實驗室前期開發出的引物［28］，SSR擴增在BIO-RAD公司生產BS97 MyCycler型的PCR儀上進行，反應體系（12.5 μL）中含有0.5 μL基因組DNA（含量10—25 ng），6.3 μL Taq MasterMix（天根），正向引物和反向引物各0.5 μL，4.7 μL ddH2O。SSR擴增程序：94℃預變性3 min，然后94℃變性30 s，57℃退火30 s，72℃延伸30 s，27個循環，最后72℃延伸7 min。擴增產物用10%聚丙烯酰胺凝膠進行分離，銀染技術檢測。

1.3統計分析

依據SSR擴增產物在電泳凝膠上的相對位置，對擴增條帶進行記錄，分別以“1”和“0”表示同一位置譜帶的有無，建立原始數據陣，用NTSYS-pc 2.10e軟件［30］進行后續分析。

計算多態性位點百分率P（%）=（k/n）×100，其中k為多態位點數，n為所測位點總數。SSR位點的多態性信息量（poly-morphism information content，PIC）按如下計算公式進行計算：PIC=1-ΣPi2，式中Pi表示第i個等位位點出現的頻率［31］。

利用NTSYS中的Qualitative date模塊計算任意兩個品種間的相似系數（GS）。其計算公式為GS=2Nij/ （Ni+Nj），其中Nij為材料i和j共有的擴增片段總數，Ni為材料i中出現的擴增片段數，Nj為材料j中出現的擴增片段數。以Clustering程序中SHAN進行UPGMA（非加權組平均法）聚類分析，并通過Treeplot模塊生成聚類圖。

將數據轉換成Structure默認的數據格式，并通過Structure軟件進行群體遺傳結構分析。將K值設為2 —20，采用馬爾可夫鏈的蒙特卡洛模擬算法（markov chain monte carlo，MCMC），不作數迭代（length of burn-in period）設為10 000次，再將不作數迭代后的MCMC設為100 000次。但依此得到的數據所繪制的折線圖常常不能提供真正的 K值。因此，在此基礎上推算對數變化率ΔK來更為準確的反應真正的K值［32］。ΔK的計算公式為：ΔK=m（|L（K+1）-2L（K）+L（K-1）|）/s［L（K）］。L（K）為每個K對應的對數值，s為標準差，m為平均值。根據ΔK的變化規律來確定最適的群體數目，構建遺傳結構圖。計算群體中品種基因組變異源于該群體中的概率Q，用于分析其遺傳成分［33］。

2　結果

2.1高效率 SSR引物的篩選與多態性位點百分率 P的計算

利用供試材料中30份棗品種的基因組DNA為模板，從課題組前期設計并篩選出的SSR引物［28］中進行二次引物篩選。根據擴增效果，從中選取23對多態性和重復性好的SSR引物（表2，圖1），用于255份棗種質的SSR分析。

23對引物在255份供試棗種質中共擴增出117條多態性條帶，多態位點百分率達 100%。各引物擴增的多態性條帶數為2—10條，平均每對引物擴增5.09個多態性條帶，擴增片段范圍90—290 bp（圖2）。BOTSTEIN等［31］提出了衡量基因變異程度高低的多態性信息量指標PIC；PIC＞0.5時，該引物為高度多態性信息引物；0.25＜PIC＜0.5，為中度多態性信息引物；PIC＜0.25，為低度多態性信息引物。本研究中篩選出的引物 PIC值變幅為 0.359—0.727（平均值0.548），均為高度或重度多態性信息引物（表2）。

表1　供試品種材料Table 1 Cultivars used in the study

編號Code名稱Cultivar原產地Origin編號Code編號Code名稱Cultivar原產地Origin名稱Cultivar原產地Origin 76　新鄭長雞心棗（D）Xinzhengchangjixinzao河南He’nan 102　喀什葛爾小棗（B）Kashigeerxiaozao山西Shanxi 128　朝陽齊脆棗（F）Chaoyangqicuizao遼寧Liaoning 77　南京棗（P）Nanjingzao浙江Zhejiang 103　吾庫扎克小棗（B）Wukuzhakexiaozao山西Shanxi 129　朝陽棉套棗（D）Chaoyangmiantaozao遼寧Liaoning 78　保定斑棗（F）Baidingbanzao河北Hebei 104　新疆小圓棗（D）Xinjiangxiaoyuanzao山東Shandong 130　鮑莊大鈴棗（D）Baozhuangdalingzao山東Shandong 79　合陽鈴鈴棗（F）Heyanglinglingzao陜西Shannxi 105　泰安靈棗（D）Tai’anlingzao河南He’nan 131　鮑莊小棗（B）Baozhuangxiaozao山東Shandong 80　延川狗頭棗（B）Yanchuangoutouzao陜西Shannxi 106　泰安疙瘩棗（F）Tai’angedazao山東Shandong 132　鮑莊小圓棗（B）Baozhuangxiaoyuanzao山東Shandong 81　雞心棗（D）Jixinzao河南He’nan 107　泰安大脆棗（F）Tai’andacuizao湖南Hu’nan 133　北京31（U）Beijing 31北京Beijing 82　北京纓絡棗（B）Beijingyingluozao北京Beijing 108　河津晉棗（F）Hejinjinzao湖南Hu’nan 134　鮑莊尖棗（U）Baozhuangjianzao山東Shandong 83　武鄉碳棗（D）Wuxiangyazao山西Shanxi 109　萬榮脆棗（B）Wanrongcuizao湖南Hu’nan 135　朝陽無核棗（F）Chaoyangwuhezao遼寧Liaoning 84　武鄉甜棗（D）Wuxiangtianzao山西Shanxi 110　萬榮耬疙瘩（D）Wanronglougeda湖南Hu’nan 136　北京大脆棗（F）Beijingdacuizao北京Beijing 85　太原十月紅（D）Taiyuanshiyuehong山西Shanxi 111　臨猗圓鈴棗（D）Linyiyuanlingzao湖南Hu’nan 137　太谷郎棗（B）Taigulangzao山西Shanxi 86　阜陽螞蟻棗（B）Fuyangmayizao安徽Anhui 112　定襄小棗（B）Dingxiangxiaozao山西Shanxi 138　太谷壺瓶酸（F）Taiguhupingsuan山西Shanxi 87　西雙版納小棗（D）Xishuangbannaxiaozao云南Yunnan 113　寧陽宣鈴棗（F，P）Ningyangxuanlingzao新疆Xinjiang 139　榆次面棗（U）Yucimianzao山西Shanxi 88　昆明棗（U）Kunmingzao云南Yunnan 114　鎮平九月寨（D）Zhenpingjiuyuezhai新疆Xinjiang 140　榆次奶頭棗（F）Yucinaitouzao山西Shanxi 89　云南2號（U）Yunnan 2云南Yunnan 115　成武冬棗（F）Chengwudongzao新疆Xinjiang 141　臨猗珍珠棗（U）Linyizhenzhuzao山西Shanxi 90　山西哈八棗（D）Shanxihabazao山西Shanxi 116　溆浦小果算盤棗（F）Xupuxiaoguosuanpan山東Shandong 142　榆次長木棗（D）Yucichangmuzao山西Shanxi 91　永濟雞蛋棗（F）Yongjijidanzao山西Shanxi 117　溆浦巖棗（D）Xupuyanzao山東Shandong 143　太谷沒心紅（B）Taigumeixinhong山西Shanxi 92　山西小令棗（U）Shanxixiaolingzao山西Shanxi 118　溆浦香棗（D）Xupuxiangzao山東Shandong 144　大饃棗（U）Damozao山西Shanxi 93　湖北鈴鐺棗（U）Hubeilingdangzao湖北Hubei 119　溆浦木棗（D）Xupumuzao山西Shanxi 145　交城牙棗（D）Jiaochengyazao山西Shanxi 94　河北圓棗（B）Hebeiyuanzao河北Hebei 120　溆浦米棗（F）Xupumizao山西Shanxi 146　太谷沒心核（U）Taigumeixinhe山西Shanxi 95　湖北雞心棗（U）Hubeijixinzao湖北Hubei 121　臨猗甜酸棗（F）Linyitiansuanzao山西Shanxi 147　晉棗1號（B）Jinzao 1山西Shanxi 96　西安羊奶棗（F）Xi’anyangnaizao陜西Shannxi 122　臨猗脖脖棗（D）Linyibobozao山西Shanxi 148　新蔡大圓豐棗（U）Xincaidayuanfengzao河南He’nan 97　大荔林檎棗（B）Dalilinqinzao陜西Shannxi 123　河津木棗（D）Hejinmuzao山西Shanxi 149　駐馬店羊角棗（F）Zhumadianyangjiaozao河南He’nan 98　佳縣細腰腰棗（D）Jiaxianxiyaoyaozao陜西Shannxi 124　北京馬牙棗（F）Beijingmayazao北京Beijing 150　保定月光（F）Baodingyueguang河北Hebei 99　溆浦薄皮棗（F）Xupubopizao湖南Hu’nan 125　泰安金絲（B）Tai’anjinsi山東Shandong 151　運城蛤蟆棗（F）Yunchenghamazao山西Shanxi 100　猴頭棗（B）Houtouzao山東Shandong 126　山東梨棗（F）Shandonglizao山東Shandong 152　運城脆棗（F）Yunchengcuizao山西Shanxi 101　滕州落地紅（D）Tengzhouluodihong山東Shandong 127　慶云小梨棗（F）Qingyunxiaolizao山東Shandong 153　榆次晚紅棗（U）Yuciwanhongzao山西Shanxi

?

F：鮮食 Fresh；D：制干 Dry；B：兼用 Fresh and dry；P：加工 Processing；O：觀賞 Ornamental；U：不詳 Unknown

圖1　篩選出的引物（JSSR97、JSSR128）對30份棗材料的擴增結果Fig.1 Amplified results of 30 jujube cultivars with selected primer pairs JSSR97 and JSSR128

圖2　引物JSSR183對255份棗材料的擴增結果Fig.2 Amplified results of 255 jujube cultivars with primer JSSR183

表2　篩選出的23對多態性SSR引物Table 2 The 23 SSR primer pairs selected for this study

2.2部分棗品種特有SSR標記的挖掘

供試棗品種的SSR分析結果以Excel表格形式保存，不同品種的SSR指紋以SSR引物名稱為前綴、以該標記（擴增帶）的分子量為后綴，得到每個品種在某個標記的帶型編號。在255個供試品種中，有40個棗品種只需1—2個特有的SSR標記即可鑒別出來（表3）。從表1和表3可以看出，很多同一區域來源的棗品種（如同樣來源于陜西省關中平原東北部蒲城縣的‘蒲城直社疙瘩棗’‘蒲城直社疙瘩棗’和‘蒲城晉棗’可以被1—2個特有SSR標記區分開，這為地方品種的準確鑒定和資源管理提供了有效手段。

表3　40個棗品種的特有SSR指紋Table 3 Special SSR fingerprints of 40 jujube cultivars

2.3供試棗品種的聚類分析

利用23對引物擴增出的117條多態性條帶，可將供試255個品種的251個品種分開，占98.43%，依據23對SSR引物擴增位點構建的0-1矩陣，計算出供試材料間的配對遺傳相似系數（GS），其變化范圍為0.71 —1.00。由聚類圖（圖3）可看出，在聚類閥值為0.71處，北京花生棗單獨聚到一類，說明其與其他棗品種親緣關系較遠，分析原因可能是北京花生棗是獨立起源的品種，也可能是其在長期栽培中與其他棗品種基因交流較少所致；其余254個品種分為14個亞類，包括4個大亞類（每類24—56個品種）和10個小亞類（1—10個品種）。在10個小亞類中，湖南的‘溆浦巖棗’和陜西的‘合陽鈴鈴棗’亦均單獨為一個小類；同時有少部分品種沒有分開，如‘奉節雞蛋棗’和‘溆浦雞蛋棗’、‘陜西奶棗’和‘天津大馬牙棗’之間的相似系數為1.00，表明‘奉節雞蛋棗’和‘溆浦雞蛋棗’、‘陜西奶棗’和‘天津大馬牙棗’的親緣關系很近。‘奉節雞蛋棗’和‘溆浦雞蛋棗’分別來自重慶、湖南，這兩個品種棗果的果形都是近圓形；‘陜西奶棗’和‘天津大馬牙棗’分別來自陜西和天津，果形均較尖，本試驗中檢測出的兩兩品種間較近的親緣關系也可能是不同地區間栽培引種造成的，也可能是兩兩品種是同物異名，或者是本研究中所用引物數量較少，尚不能將這些近緣品種完全分開。

從聚類圖中還可以看出：（1）一些類群中，棗品種間的遺傳關系與其原產地有明顯的相關性。如I和II組的棗品種基本都是來自山西，II組的棗品種基本都來自北京，而O組的基本來自湖南；來自陜西的‘大荔鈴鈴棗’‘大荔蜂蜜罐’和‘大荔疙瘩棗’聚在一起，來自山東的‘泰安大脆棗’‘酥脆棗’和‘泰安疙瘩棗’也都聚在一起。（2）另一些類群中（如C、E、K組），棗品種來自于不同的產地，遺傳關系與地理分布似乎沒有必然的聯系。這些類群相對比較分散，每一小類群中品種數量不多。如Ⅳ類群中7個棗品種來自于6個不同地區；Ⅴ類群中20個品種主要來自河北、山西和山東；Ⅵ類群中品種主要來自江蘇和湖南。（3）結合聚類圖和供試品種的用途和原產地分析，不同棗品種間的親緣關系與品種原產地有一定相關性，但和品種用途沒有顯著相關性，例如第一大類中就包括了制干、鮮食、加工等不同用途的品種。

2.4供試棗種質的遺傳結構分析

基于上述獲得的SSR分析數據，利用STRUCTURE軟件分析供試棗群體的遺傳結構，繪制K與ΔK的關系圖（圖4），K=15時，ΔK最大。據此，將255個棗品種分成15個群組，進一步以K=15繪制了基于數學模型模擬的255個棗品種遺傳結構的貝葉斯分配圖（圖5）。并計算出了15個居群中最大Q值分布（表4）。

表4　各群體Q值分布Table 4 Distribution of Q-value of fifteen groups

圖3　基于SSR分析的255個棗品種聚類圖Fig.3 Dendrogram of 255 jujube cultivars based on SSR markers

圖4　棗255個品種的K值與ΔK的關系Fig.4 Graphical relationship between K and ΔK for 255 jujube cultivars

圖5　255個棗品種的遺傳結構分析Fig.5 Population structure of 255 jujube cultivars

圖5中不同的色塊代表不同的群體，縱坐標表示各居群中的品種占某居群祖先成分的比例，橫坐標為樣品編號。從圖5中可看出，根據各個品種的混合比率，將255份供試材料分成15個居群，進一步分析表明，一些居群（如3、8、9、12、14居群）遺傳背景比較一致，絕大部分來源于相同的栽培區域，混雜程度相對較低；而有的居群（如2、4、5、13居群）分配率較高，群體間混合程度相對較高，來源于不同的栽培區域，群體間的遺傳物質交流相對較多。從遺傳結構分析也發現，不同群體的品種和其用途沒有顯著相關性，但與聚類分析（圖3）的結果相似。

Structure群體遺傳結構分析中，當某一材料在某類群中的Q≥0.6時，則認為該材料血緣關系相對比較單一，否則認為該材料血緣關系來源復雜［33］。從根據供試材料在各個群體中的Q值分布（表4），255個供試材料中有36個品種Q值小于0.6，占所有供試材料總數的14.1%，Q值大于0.8和0.9的品種分別占68.2%和56.5%，說明各群體中大部分品種血緣關系比較單一，較少品種含有其他類群的遺傳成分。整體看，來自棗主栽區域山西（居群 6、10、11）、陜西（居群14）和河北（居群1）的品種形成了以當地品種為主的居群，說明這些地區在長期的栽培演化過程中選育出了豐富的棗品種資源，其中來自山西或陜西的棗品種出現在絕大部分居群中，說明這兩個省的資源在不同群體間的基因交流中發揮了重要作用；南方區域只有來自湖南的棗品種形成了相對獨自的居群，可能是其起源相對單一，且在較長一段時間的栽培過程中，栽培環境比較封閉，近親繁殖，和其他產地棗品種間基因交流較少所致（居群8）。

3　討論

本研究選用23對SSR引物對255個棗品種進行了分析，共擴增出117條多態性譜帶，多態位點比率達到100%。在棗樹上，麻麗穎等［8］報道的SSR標記多態性比率也為 100%，而李莉等［34］報道的RAPD標記多態性為 99.78%，白瑞霞［35］用 AFLP標記的多態性則僅為 64.30%，謝永波［36］用 AFLP標記的多態性為78.7%，說明本研究的SSR標記的多態性效率較高。分析原因一方面是因為本研究中所用引物經過了二次篩選，另一方面也是因為供試材料多、遺傳多樣性豐富而使得多態位點比率較高。同時，本試驗也注意到，不同 SSR引物對棗品種的鑒別能力有較大差異，并不是所有引物的多態性都較高，進一步說明篩選高效率引物的必要性。

中國棗品種繁多，不同品種間復雜的遺傳關系一直是研究熱點。本研究利用 SSR技術進行了 22個省區主產地 255個棗品種的遺傳多樣性和群體結構分析，提供了大量的多態性數據信息，為將來表型數據和基因數據的關聯性分析奠定了良好的基礎。聚類和群體結構兩種分析結果均表明，在棗品種的群體結構中大部分品種與來源區域有一定的相關性。相關性不明顯的棗品種可能是區域間品種經過頻繁的基因交流或重組培育而成，從而融合了不同區域品種的特點，沒有了明顯的區域特征。前人利用SSR［17］、ISSR［18］、AFLP［2］進行的研究也表明，遺傳多樣性在一定程度上能夠反映品種的地理分布特點。李海濤［26］對河南主栽棗品種的 ISSR分析表明，棗的遺傳變異主要來自于居群內的變異，來自于居群間的變異貢獻率較低。殷曉［27］對黃河沿岸六省48個棗品種進行的遺傳多樣性研究表明，地理區域內變異（93%）遠大于地理區域間分化（7%）。本研究中，255個供試材料中Q值小于0.6的品種36個，僅占14.1%（表4），也進一步說明大部分品種血緣關系比較單一，以居群內或地理區域內變異為主。這些結果均表明，現有的棗種質資源可能是不同區域起源的，在物種演化過程中區域內品種間進行了頻繁的基因交流和重組，進而產生了更為豐富的資源結構組成。

4　結論

利用23對SSR引物在255個供試品種中共檢測出117個多態性位點，平均每對引物擴增多態位點數為5.09，PIC值變幅為0.359—0.727，平均為0.548；為其中40個棗品種找到了1—2個各自特有的SSR標記；255個棗品種的群體遺傳結構與聚類分析結果基本一致，不同地理環境在棗品種的群體進化中發揮了較重要的作用，但不同棗品種的遺傳關系與品種用途沒有顯著相關性。

References

［1］ 趙錦，劉孟軍. 棗樹品種，品系及其近緣種的RAPD分析. 中國農業科學，2003，36（5）: 590-594. ZHAO J，LIU M J. RAPD analysis on the cultivars，strains and relatedspecies of Chinese jujube （Ziziphus jujuba Mill.）. Scientia Agricultura Sinica，2003，36（5）: 590-594. （in Chinese）

［2］ 喬勇，趙錦，楊海旭，劉孟軍. 21個棗品種（系）的AFLP指紋分析. 植物遺傳資源學報，2009，10（2）: 205-210. QIAO Y，ZHAO J，YANG H X，LIU M J. AFLP analysis on 21 cultivars and strains of Chinese jujube. Journal of Plant Genetic Resources，2009，10（2）: 205-210. （in Chinese）

［3］ Union for the Protection of New Varieties of Plants. Guidelines for DNA-profiling: Molecular marker selection and database construction. Switzerland，2007: 3-4.

［4］ 張萌. 基于SSR分子標記的葡萄種質資源遺傳多樣性分析及品種鑒定［D］. 江蘇: 南京農業大學，2012. ZHANG M. Analysis of genetic diversity and cultivar identification of grapevine germplasm resources based on SSR molecular markers ［D］. Jiangsu: Nanjing Agricultural University，2012. （in Chinese）

［5］ 徐小彪，廖嬌，黃春輝，辜青青，曲雪艷，劉善軍，陳金印. 基于EST-SSR標記分析獼猴桃種質遺傳關系. 果樹學報，2012，29（2）: 212-216. XU X B，LIAO J，HUANG C H，GU Q Q，QU X Y，LIU S J，CHEN J Y. Genetic relationships from Kiwifruit germplasms based on EST-SSR markers. Journal of Fruit Science，2012，29（2）: 212-216. （in Chinese）

［6］ 王玉安，歐巧明，陳建軍，胡達，王發林. 甘肅地方杏品種資源的SSR遺傳多樣性分析. 西北農業學報，2013，22（3）: 98-102. WANG Y A，OU Q M，CHEN J J，HU D，WANG F L. Analysis of genetic diversities in Gansu local apricot varieties with simple sequence repeat （SSR） markers. Acta Agriculturae Borealioccidentalis Sinica，2013，22（3）:98-102. （in Chinese）

［7］ 巴巧瑞. 蘋果栽培品種親緣關系的SSR和SRAP分析研究［D］. 楊陵: 西北農林科技大學，2011. BA Q R. Genetic relationship analysis of apple cultivars using SSR and SRAP markers ［D］. Yangling: Northwest A & F University，2011. （in Chinese）

［8］ 麻麗穎，孔德倉，劉華波，王斯琪，李穎岳，龐曉明. 36份棗品種SSR指紋圖譜的構建. 園藝學報，2012，39（4）: 647-654. MA L Y，KONG D C，LIU H B，WANG S Q，LI Y Y，PANG X M. Construction of SSR fingerprint on 36 Chinese jujube cultivars. Acta Horticulturae Sinica，2012，39（4）: 647-654. （in Chinese）

［9］ ZHANG C Y，CHEN X S，HE T M，LIU X L，FENG T，YUAN Z H. Genetic structure of Malus sieversii population from Xinjiang，China，revealed by SSR markers. Journal of Genetics and Genomics，2007，34（10）: 947-955.

［10］ 張春雨，陳學森，林群，苑兆和，張紅，張小燕，劉崇祺，吳傳金.新疆野蘋果群體遺傳結構和遺傳多樣性的 SRAP分析. 園藝學報，2009，36（1）: 7-14. ZHANG C Y，CHEN X S，LIN Q，YUAN Z H，ZHANG H，ZHANG X Y，LIU C Q，WU C J. SRAP markers for population genetics structure and genetic diversity in Malus sieversii from Xinjiang，China. Acta Horticulturae Sinica，2009，36（1）: 7-14. （in Chinese）

［11］ 劉晶. 中國豆梨與川梨的遺傳多樣性和群體遺傳結構研究［D］. 杭州: 浙江大學，2013. LIU J. Studies on genetic diversity and structure of Pyrus calleryana and P. pashia in China ［D］. Hangzhou: Zhejiang University，2013. （in Chinese）

［12］ 章秋平，劉東成，劉威生，劉碩，張愛民，劉寧，張玉萍. 華北生態群普通杏遺傳多樣性與群體結構分析. 中國農業科學，2013，46（1）: 89-98. ZHANG Q P，LIU D C，LIU W S，LIU S，ZHANG A M，LIU N，ZHANG Y P. Genetic diversity and population structure of the north China populations of apricot （Prunus armeniaca L.）. Scientia Agricultura Sinica，2013，46（1）: 89-98. （in Chinese）

［13］ 張俊環，王玉柱，孫浩元，楊麗. 不同用途杏品種群的AFLP分析.果樹學報，2011，28（4）: 610-616. ZHANG J H，WANG Y Z，SUN H Y，YANG L. Genetic analysis of apricot （Armeniaca） by fluorescent-AFLP markers. Journal of Fruit Science，2011，28（4）: 610-616. （in Chinese）

［14］ YUAN Z H，CHEN X S，HE T M，FENG J R，FENG T，ZHANG C Y. Population genetic structure in apricot （Prunus armeniaca L.）cultivars revealed by fluorescent-AFLP markers in southern Xinjiang，China. Journal of Genetics and Genomics，2007，34（11）: 1037-1047.

［15］ 何天明，陳學森，高疆生，張大海，徐麟，吳燕. 新疆栽培杏群體遺傳結構的SSR分析. 園藝學報，2006，33（4）: 809-812. HE T M，CHEN X S，GAO J S，ZHANG D H，XU L，WU Y. Using SSR markers to study population genetic structure of cultivated apricots native to Xinjiang. Acta Horticulturae Sinica，2006，33（4）: 809-812. （in Chinese）

［16］ 苑兆和，陳學森，張春雨，何天明，馮建榮，馮濤. 普通杏群體遺傳結構的熒光AFLP分析. 園藝學報，2008，35（3）: 319-328. YUAN Z H，CHEN X S，ZHANG C Y，HE T M，FENG J R，FENG T. Population genetic structure in apricot （Armeniaca L.） cultivars revealed by fluorescent-AFLP markers. Acta Horticulturae Sinica，2008，35（3）: 319-328. （in Chinese）

高含硫天然氣凈化廠中控室暖通系統技術優化與應用評價…………………………………………………………(1):110

［17］ 宛甜，蔡宇良，馮瑛，張雪，何恒流. 野生毛櫻桃 SSR遺傳多樣性和遺傳結構分析. 西北植物學報，2013，33（8）: 1544-1550. WAN T，CAI Y L，FENG Y，ZHANG X，HE H L，WAN T. Genetic diversity and genetic structure of wild Prunus tomentosa Thub. based on simple sequence repeats markers. Acta Botanica Boreali-Occidentalia Sinica，2013，33（8）: 1544-1550. （in Chinese）

［18］ 賴恭梯，劉煒婳，張梓浩，馮新，林玉玲，劉生財，祁芳斌，賴鐘雄.戴云山野生楊梅自然群體遺傳結構的 ISSR分析. 熱帶作物學，2013，34（10）: 1863-1870. LAI G T，LIU W H，ZHANG Z H，FENG X，LIN Y L，LIU S C，QI F B，LAI Z X. Genetic structure of natural populations revealed by ISSR in Myrica rubra from daiyunshan mountains. Chinese Journal of Tropical Crops，2013，34（10）: 1863-1870. （in Chinese）

［19］ 昝逢剛，吳轉娣，曾淇，張惠云，李明芳，鄭學勤. 荔枝種質遺傳多樣性的SRAP分析. 分子植物育種，2009，7（3）: 562-568. ZAN F G，WU Z D，ZENG Q，ZHANG H Y，LI M F，ZHENG X Q. Genetic diversity analysis of litchi germplasm by SRAP markers. Molecular Plant Breeding，2009，7（3）: 562-568. （in Chinese）

［20］ 余賢美，艾呈祥. 杧果野生居群遺傳多樣性ISSR分析. 果樹學報，2007，24（3）: 329-333. YU X M，AI C X. Genetic diversity of wild Mangifera indica populations detected by ISSR. Journal of Fruit Science，2007，24（3）: 329-333. （in Chinese）

［21］ 劉亞令，李作洲，姜正旺，劉義飛，黃宏文. 中華獼猴桃和美味獼猴桃自然居群遺傳結構及其種間雜交漸滲. 植物生態學報，2008，32（3）: 704-718. LIU Y L，LI Z Z，JIANG Z W，LIU Y F，HUANG H W. Genetic structure and hybridization introgession in natural populations of two closely related Actinidia species，A. Chinensis and A. Deliciosa. Journal of Plant Ecology （Chinese Version），2008，32（3）: 704-718. （in Chinese）

［22］ 劉曉麗，陳學森，張美勇，陳曉流，何天明，張立杰，張春雨. 普通核桃（Juglans regia） 3個群體遺傳結構的SSR分析. 果樹學報，2008，25（4）: 526-530. LIU X L，CHEN X S，ZHANG M Y，CHEN X L，HE T M，ZHANG L J，ZHANG C Y. Population genetic structure analysis of Juglans regia using SSR markers. Journal of Fruit Science. 2008，25（4）: 526-530. （in Chinese）

［24］ YUAN Z H，YIN Y L，QU J L，ZHU L Q，LI Y. Population genetic diversity in Chinese pomegranate （Punica granatum L.） cultivars revealed by fluorescent-AFLP markers. Journal of Genetics and Genomics，2007，34（12）: 1061-1071.

［25］ 呂志江，李疆，吾買爾夏提·塔漢，曾斌，羅淑萍. 新疆野扁桃種質資源遺傳多樣性的ISSR分析. 果樹學報，2010，27（6）: 918-923. Lü Z J，LI J，OMIR S T，ZENG B，LUO S P. ISSR analysis for genetic diversity of Amygdalus ledebouriana germplasm from xinjiang，China. Journal of Fruit Science，2010，27（6）: 918-923. （in Chinese）

［26］ 李海濤. 河南棗主栽品種及灰棗群體遺傳變異分析［D］. 鄭州: 河南農業大學，2008. LI H T. Analysis of genetic variation in Henan main cultivated Ziziphus jujuba varieties and Z. ‘Huizao’ population ［D］. Zhengzhou: Agricultural University of Henan，2008. （in Chinese）

［27］ 殷曉. 基于 SSR標記的中國棗遺傳多樣性研究［D］. 楊凌: 西北農林科技大學，2013. YIN X. Genetic diversity and population structure of Chinese jujube analysed by SSR markers ［D］. Yangling: Northwst A & F University，2013. （in Chinese）

［28］ XIAO J，ZHAO J，LIU M J，LIU P，DAI L，ZHAO Z H. Genome-wide characterization of simple sequence repeat （SSR） loci in Chinese jujube and jujube SSR primer transferability. PLoS One，2015，10（5）: e0127812.

［29］ 肖京. 棗基因組 SSR位點特征分析及引物開發［D］. 保定: 河北農業大學，2014. XIAO J. Characterization of SSR loci in jujube genome and development of SSR primers ［D］. Baoding: Agricultural University of Hebei，2014. （in Chinese）

［30］ ROHLF F J. NTSYS-PC: Numerical taxonomy and multivariate analysis system version 1.80. Setauket，New York: Distribution by Exeter SoftWare，1994.

［31］ BOTSTEIN D，WHITE R L，SKOLNICK M，DAVIS R W. Construction of a genetic linkage map in man using restriction fragment length polymorphism. American Journal of Human Genetics，1980，32（3）: 314-331.

［32］ EVANNO G，REGNAUT S，GOUDET J. Detecting the number of clusters of individuals using the software structure: A simulation study. Molecular Ecology，2005，14（8）: 2611-2620.

［33］ 吳承來，張倩倩，董炳雪，李圣福，張春慶. 我國部分玉米自交系遺傳關系和遺傳結構解析. 作物學報，2010，36（11）: 1820-1831.WU C L，ZHANG Q Q，DONG B X，LI S F，ZHANG C Q. Analysis of genetic structure and genetic relationships of partial maize inbred lines in China. Acta Agronomica Sinica，2010，36（11）: 1820-1831. （in Chinese）

［34］ 李莉，彭建營，白瑞霞. 中國棗屬植物親緣關系的RAPD分析. 園藝學報，2009，36（4）: 475-480. LI L，PENG J Y，BAI R X. Studies on the phylogenetic relationship of Chinese Ziziphus by RAPD technique. Scientia Agricultura Sinica，2009，36（4）: 475-480. （in Chinese）

［35］ 白瑞霞. 棗種質資源遺傳多樣性的分子評價及其核心種質的構建［D］. 保定: 河北農業大學，2008. BAI R X. Studies on genetic diversity and core collection construction of Ziziphus jujuba germsplasm resources using AFLP and SRAP markers. Baoding: Agricultural University of Hebei，2008. （in Chinese）

［36］ 謝永波. 棗屬種質資源形態學評價及品種 AFLP鑒定［D］. 泰安:山東農業大學，2014. XIE Y B. Morphological evaluation on germplasm resources jujube and cultivar identification by AFLP ［D］. Tai’an: Shangdong Agricultural University，2014. （in Chinese）

（責任編輯 趙伶俐）

Genetic Diversity and Structure of 255 Cultivars of Ziziphus jujuba Mill.

LIU Xiu-yun1，LI Hui1，2，LIU Zhi-guo1，ZHAO Jin3，LIU Meng-jun1，4
（1Reserch Center of Chinese Jujube，Agricultural University of Hebei，Baoding 071001，Hebei；2Germplasm Resources Institute of Ningxia Academy of Agriculture and Forestry Sciences，Yinchuan 750002；3College of Life Science，Agricultural University of Hebei，Baoding 071001，Hebei；4Beijing Collaborative Innovation Center for Eco-Enviromental Improvement with Forestry and Fruit Trees，Beijing 102206）

【Objective】 There are abundant jujube germplasm resources in China. A total of 255 cultivars of Ziziphus jujuba Mill. from 22 provenances were used as materials to reveal their genetic diversity and phylogenetic relationship by SSR analysis，and the results of analysis would help us to manage jujube germplasm resources and offer references for molecular maker-assisted breeding. 【Method】 Good genomic DNA was extracted from young leaves of jujube germplasm resources following the improved CTAB method，and then were amplified by simple sequence repeat molecular markers to analyze genetic diversity and genetic structure with the selected high-efficiency primer pairs which were excavated based on the genome sequencing. Separation of the amplified fragments was performed on 8% denaturing polyacrylamide gels and the gels were stained with AgNO3for visualizing the SSR fragments. The data were counted by presence or absence of the band and the percentage of polymorphic loci （PIC） was calculated. UPGMA cluster analysis was carried by software NTSYS，the optimal number of groups and population genetic structure was analyzed by software Structure.【Result】Totally，117 polymorphic alleles were revealed with 23 primer pairs which was selected from 64 primer pairs，each primer amplified polymorphic loci ranged from 2 to 10，with an average of 5.09 for each primer pairs. Polymorphism information content （PIC） values for the primer pairs ranged from 0.359 to 0.727，with an average of 0.548，these polymorphisms primers could be further applied to other study. The fingerprint for some jujube cultivars was established with 1-2 markers，providing a reference for the management of jujube germplasm. Meanwhile，based on the UPGMA cluster analysis，255 cultivars were divided into fifteen subgroups，which included four big groups and eleven small groups. Similarity coefficients among the cultivars were between 0.71 to 1.00，‘Beijinghuashengzao’ was clustered into one separate group，which has a distant relationship with other cultivars. The similarity coefficients of ‘Fengjiejidanzao’ and ‘Xupujidanzao’，‘Shannxinaizao’ and ‘Tianjindamayazao’ were both 1.00. In some subgroups the genetic relationship between cultivars and their provenances has a significant positive correlation，but the cultivars and their uses has no significant correlation. Based on K and ΔK values，255 jujube cultivars were also divided into fifteen populations by the population genetic structure analysis. The kinship among cultivars in the same population was relatively simple，and a few cultivars contained genetic component of other groups. The cultivars from Shanxi or Shannxi were distributed in most populations，indicating jujube cultivars of the two provinces played important roles in the gene exchange among populations. The jujube cultivars from Hunan of the South region formed a relatively alone population，indicating that the cultivars might be from the same source，or in the long-term cultivation few times of gene exchange were happened in Hunan cultivars with other populations. Different geographical environment played a key role in the evolution of jujube germplasm populations，some cultivars were selected from the same geographical environment and the others were selected by genetic recombination among those cultivars from the various geographical environment. Meanwhile，the consistency of the two different methods was further verifed the accuracy of the results，which provide useful clues and reference for the genetic diversity and structure of jujube germplasm.【Conclusion】Geographical environment play significant roles in the population evolution of jujube cultivars，affecting the genetic structure composition between different habitats.

Chinses jujube； SSR markers； provenances； genetic diversity； population structure

2016-01-18；接受日期：2016-04-13

“十二五”國家科技支撐計劃（2013BAD14B03）、北京市2011協同創新中心項目（PXM2016_014207_000038）、河北農業大學青年學術帶頭人項目

聯系方式：劉秀云，E-mail：jingxianyu705@126.com。李慧，E-mail：1063979927@qq.com。劉秀云和李慧為同等貢獻作者。通信作者劉孟軍，E-mail：lmj1234567@aliyun.com。通信作者趙錦，E-mail：zhaojinbd@126.com