基于SSR標記的黍稷種質資源遺傳多樣性及親緣關系研究

薛延桃，陸平，喬治軍，劉敏軒，王瑞云1,

?

基于SSR標記的黍稷種質資源遺傳多樣性及親緣關系研究

薛延桃1,2，陸平2，喬治軍3，劉敏軒2，王瑞云1,3

（1山西農業大學農學院，山西太谷 030801；2中國農業科學院作物科學研究所，北京 100081；3山西省農業科學院農作物品種資源研究所/農業部黃土高原作物基因資源與種質創制重點實驗室/雜糧種質資源發掘與遺傳改良山西省重點實驗室，太原 030031）

【目的】利用SSR標記，分析黍稷種質資源（野生材料和地方品種）的遺傳多樣性水平，揭示不同來源黍稷種質資源的親緣關系和遺傳群體結構差異，為黍稷起源進化研究奠定基礎。【方法】用6份地理差異顯著的黍稷種質資源對137對小宗作物課題組開發的具有多態性的SSR引物進行初步篩選，最終篩選103對條帶清晰、擴增良好且多態性穩定的SSR引物，利用這103對多態性SSR標記對146份黍稷材料進行PCR擴增，通過遺傳參數、聚類、遺傳結構等分析，評估不同個體間及不同群體間的遺傳多樣性，探討遺傳結構差異。【結果】103對SSR標記共檢測出308個等位基因（Na），平均值為2.99，平均Shannon-Weaver指數（I）為0.8478，平均期望雜合度為0.3642，平均多態性信息含量指數（）為0.5544。103對SSR標記的分布區間為0—1、1—2、2—3、3—4和4—5，分辨率范圍為0.334—4.002，77.67%的標記分布于區間1—4，具有適度分辨力。國內資源的觀測等位基因數（2.9126）、多樣性指數（0.8302）、期望雜合度（0.5023）、多態性信息含量指數（0.5278）均高于國外資源，遺傳多樣性更豐富。12個群體的遺傳距離的變化范圍為0.0783—0.5762，均值為0.2938；遺傳一致度變化范圍為0.5620—0.9247，均值為0.75，遺傳相似性與地理分布具有一定相關性，地理分布越近，遺傳距離越小，遺傳一致度越高。聚類分析在遺傳距離為0.15處可以把12個群體分為4個組群，其中南美洲和山西資源各自獨立分為一支，與其他資源親緣關系較遠。個體間聚類中，國內外資源劃分非常顯著，在遺傳距離為0.63處，146份黍稷資源可分為3大組群，組群Ⅰ和組群Ⅱ為國外資源，組群Ⅲ為國內資源。組群Ⅱ在遺傳距離為0.39處又分為3個亞群，組群Ⅲ在遺傳距離為0.45處分為5個亞群，其中亞洲與歐洲資源、中國河北與中國山西、中國內蒙古資源的遺傳關系較近。遺傳結構分析結果顯示國內外群體間存在明顯的遺傳分化，其中5個組群（組群2、組群5、組群6、組群7和組群9）為國內野生資源特有基因型，分布較為分散；2個組群（組群1和組群4）為國外資源特有基因型，分布較為集中。中國寧夏、南美洲資源的群體結構趨向單一化，中國河北、中國黑龍江、亞洲資源的群體結構趨向多元化。UPGMA聚類結果與遺傳結構分析結果一致，且不同地區黍稷資源群體間遺傳關系遠近均與其地理分布相關。【結論】野生資源的遺傳多樣性高于國外資源，其中中國河北群體的遺傳多樣性最豐富，中國河北可能是黍稷的起源中心。

黍稷；野生糜子；國外品種；SSR標記；遺傳多樣性；群體結構

0 引言

【研究意義】作物種質資源是開展作物品種選育的物質基礎，稀有特效種質對推動育種進程具有決定性作用，新的育種目標能否實現決定于育種單位所擁有的種質資源，此外，種質資源也是開展生物學理論研究和分子遺傳的重要基礎材料，因此，開展種質資源的收集、評估、保存和利用對中國農業研究和發展具有十分重要的意義。黍稷（L.），屬于禾本科黍屬（L.）[1]，一年生草本植物[2]。它是最古老的農作物之一，有10 000多年的栽培歷史，在世界范圍內都有種植[3-4]。黍稷生長周期短，抗旱耐貧瘠，是干旱半干旱地區的重要作物[5-6]。中國目前保存黍稷資源8 800多份，準確評估黍稷資源在形態學和分子水平的遺傳多樣性，有利于其合理開發與利用，并促進其育種工作和基礎理論方面的研究[7]。【前人研究進展】遺傳多樣性是生物多樣性的重要組成部分，通過遺傳多樣性分析可以為各稀有物種的收集、育種應用以及起源進化研究提供有利數據。早期關于黍稷遺傳多樣性的研究多集中于形態學特征。胡興雨等[8]和董孔軍等[9]分別對國家種質庫的8 016份黍稷種質資源和甘肅省600份黍稷地方資源的11個農藝性狀進行主成分分析和聚類分析，評估了參試材料在形態學水平上的遺傳多樣性并以此為基礎構建了黍稷初級核心種質。近年來，隨著分子生物技術的快速發展，分子標記在黍稷遺傳多樣性研究方面也得到了較廣泛的應用。在不同類型的分子標記中，SSR標記具有多等位基因、共顯性和高度多態性等優點[10]，廣泛應用于作物的遺傳多樣性評價中，是評估遺傳多樣性的最佳選擇[11-13]。Cho等[14]、王銀月等[15]和Rajput等[16]采用SSR-PCR技術分別開發了25對、116對、254對多態性較好的黍稷特異性SSR標記。王瑞云等[17]開發了85對高基元SSR標記，其分辨能力強，擴增多態性好；朱建楚[18]、董俊麗等[19]將分子水平和表型結合對黍稷資源進行了遺傳多樣性分析；Hu等[20]、Hu等[21]、連帥等[22]，Liu等[23]和王瑞云等[24]分別利用SSR標記評估了不同地理來源的黍稷資源的遺傳多樣性，研究表明聚類結果與地理分布密切相關，王瑞云等[17]和Hu等[21]的研究結果表明黃土高原地區的遺傳多樣性最豐富，推測黃土高原可能是黍稷的起源地。連帥等[25]利用63對SSR多態性引物對192份黍稷資源進行遺傳多樣性分析，結果表明，國內外群體遺傳分化不明顯，群體間親緣關系較近，個體間存在相互滲透。【本研究切入點】目前，黍稷遺傳多樣性研究多數集中于栽培品種與地方品種的研究，有關黍稷野生資源和國外地方品種的遺傳多樣性，特別是分子水平群體遺傳結構的研究鮮有報道。【擬解決的關鍵問題】本研究利用本課題組開發的黍稷基因組SSR標記對近年來新引進的國內主栽區的黍稷野生資源以及國外資源共計146份材料進行SSR標記水平的遺傳多樣性分析，旨在揭示不同來源黍稷種質資源的遺傳多樣性和群體遺傳結構差異水平，為進一步開展黍稷起源地、傳播途徑以及性狀進化等研究奠定基礎。

1 材料與方法

1.1 材料

所用的146份黍稷資源（12個群體）均由中國農業科學院作物科學研究所小宗作物課題組近幾年新收集和引進，包括85份國內野生資源和新引進的61份國外地方品種（電子附表1）。

1.2 DNA的提取與檢測

選取50 mg室內花盆培養出苗20 d左右的幼葉，用DNA提取試劑盒（北京鼎國昌盛生物有限公司）提取基因組DNA。采用1%的瓊脂糖凝膠電泳檢測DNA的質量，用超微量分光光度計（P360）測定其濃度。調試其終濃度至30 ng·μL-1，-20℃保存備用。

1.3 SSR分析

選用6份地理差異顯著的黍稷種質資源對137對小宗作物課題組開發的具有多態性的SSR引物進行初步篩選，最終篩選103對條帶清晰、擴增良好且多態性穩定的SSR引物，用于后續遺傳多樣性分析。

PCR擴增反應在MY-CYCLER PCR儀上進行。PCR反應體系包括1.6 μL 10×PCR buffer、0.8 μL 2.5 mmol·L-1dNTP、0.1 μL Taq酶（5.0 U·μL-1）、0.5 μL引物、0.1 μL DNA模版和5.5 μL ddH2O。PCR擴增程序為94℃ 5 min；94℃ 45 s，55℃ 50 s，72℃ 1 min，39個循環；72℃ 10 min；4℃保存。用8%的聚丙烯酰胺凝膠電泳分離PCR產物，硝酸銀染色顯影。

1.4 數據分析

分析SSR標記擴增條帶，如果在至少2份資源中觀察到不同大小的2條DNA條帶，則SSR標記被認為是多態性的，有條帶記為“1”，無條帶記為“0”。利用以下公式計算每個SSR標記物的分辨率（resolving power，Rp值）：Rp=∑，=1-(2×︱0.5-︱)，其中為某等位基因信息量，為某等位基因在6份材料中出現的頻率。分別按照地理來源與國家來源，將146份黍稷資源分為12個群體和30個群體，利用PowerMarker3.25和MEGA7計算不同群體間和群體內各位點的值及遺傳距離，并繪制聚類圖。利用PopGen1.32計算不同群體間和群體內各位點的遺傳多樣性指標，包括觀測等位基因數（observed number of alleles，Na）、有效等位基因數（effective number of alleles，Ne）、觀測雜合度（observed heterozygosity，Ho）、期望雜合度（expected heterozygosity，He）及Shannon-Weaver指數（I）。利用Structure2.3.4對黍稷資源進行遺傳結構分析。

2 結果

2.1 SSR標記分析

使用103對SSR標記對146份糜子資源進行檢測（表1），共檢測出308個等位基因（Na），平均值為2.9903，其中有效等位基因（Ne）為236.55個，平均值為2.2966，所占比重為76.8%。檢測的等位基因數的范圍為1（BM994、BM5191和LMX635）—6（BM5038），其中，62對標記檢測到多個等位基因（3—6），38對標記檢測到2個等位基因，3對標記檢測到單個等位基因。

Shannon-Weaver指數范圍（I）為0（BM994、BM5191和LMX635）—1.5964（BM5181），平均值為0.8478.。觀測雜合度（Ho）范圍為0（BM114、BM136、BM994、BM5191和LMX635）—0.9920（BM4851），平均值為0.4871。期望雜合度（He）范圍為0（BM994、BM5191和LMX635）—0.7992（BM5181），平均值為0.3642。Nei's期望雜合度范圍為0（BM994、BM5191和LMX635）—0.7951（BM5181）。值范圍為0.0640（BM994）—0.8035（BM4956），平均值為0.5544。一些SSR標記（BM114、BM136、BM994、BM5191和LMX635）的Shannon-Weaver指數、觀測雜合度、期望雜合度、Nei's期望雜合度為0，表明個體之間存在普遍異交。

引物分辨率（Rp值）指不同標記對不同基因型的辨別能力，引物分辨率越高，其對不同基因的分辨能力越強，因此從各標記的高分辨能力也可以觀察到這些標記的高多態性。103對SSR標記的Rp值范圍為0.334（BM1533、BM4739、BM4741、BM4947、LMX635和LMX1429）—4.002（LMX1761）。Rp值的分布區間為0—1、1—2、2—3、3—4和4—5，分布頻率分別為21.36%、45.63%、26.21%、5.83%和0.97%（圖1），由以上數據可知，77.67%的標記分布區間1—4，具有適度分辨力，選擇有效的標記利于后續遺傳多樣性的分析。

表1 103對SSR標記的遺傳參數

圖1 103個SSR標記的分辨率（Rp值）

2.2 不同群體間的遺傳多樣性分析

利用PopGen1.3.2軟件分析12個不同群體間黍稷資源的遺傳多樣性（表2），結果表明，各群體間等位基因（Na）變化范圍為1.4947—2.6602，最大的為中國河北群體，最小的為南美洲群體。有效等位基因（Ne）變化范圍為1.4510—2.0961，均值為1.8037。就國內外兩大群體而言（表3），國內各群體的有效等位基因數（2.2620±0.7862）大于國外群體（2.0429±0.6601），國內群體間的有效等位基因范圍為1.6302—2.0961（河北＞黑龍江＞遼寧＞寧夏＞內蒙古＞山西），均值為1.8740，比全部參試群體的平均有效等位基因高0.0703。國外群體間的有效等位基因范圍為1.4510—1.9860（歐洲＞亞洲＞北美洲＞大洋洲＞非洲＞南美洲），均值為1.7334，比全部參試群體的平均有效等位基因低0.0703。由以上數據可知國內外群體間遺傳豐富度差異不大，國內野生糜子遺傳豐富度略大于國外地方品種。

Shannon-Weaver指數（I）范圍為0.3246—0.7608，最大的為中國河北群體，最小的為南美洲群體。國內群體Shannon-Weaver指數（0.8302）大于國外群體（0.7434），其不同群體的變化趨勢與Ne一致。從Nei’s基因雜合度分析，觀測雜合度（Ho）最大的是中國遼寧群體（0.3443），最小的是中國黑龍江群體（0.2727），均值為0.3123。期望雜合度（He）最大的是中國河北群體（0.4789），最小的是南美洲群體（0.3123），均值為0.4143。國內外不同群體的觀測雜合度和期望雜合度的變化趨勢不一致，國內群體的Ho值（0.2972±0.3030）小于國外群體（0.3172±0.3405），He值（0.5023±0.1909）大于國外群體（0.4581±0.1967）。各參試群體的值變化范圍為0.1383（南美洲）—0.4898（中國黑龍江），國內群體的值大于國外群體。

表2 12個群體的遺傳多樣性分析

表3 國內外群體的遺傳多樣性分析

綜合以上遺傳參數可知，國內資源的遺傳多樣性高于國外資源。國內群體中河北、黑龍江的遺傳多樣性較高，國外群體中亞洲、歐洲的遺傳多樣性較高。

2.3 不同群體間遺傳相似性分析

利用PopGen1.3.2軟件分析12個群體的遺傳相似性（表4），結果表明，遺傳距離的變化范圍為0.0783—0.5762，均值為0.2938。遺傳一致度變化范圍為0.5620—0.9247，均值為0.7500。歐洲群體和亞洲群體的遺傳距離最小（0.0783），遺傳一致度最大（0.9247），說明這兩個群體的親緣關系較近。而南美洲群體和中國山西群體的遺傳距離最大（0.5762），遺傳一致度最小（0.5620），兩者的親緣關系較遠。由以上分析可知，遺傳相似性與地理分布具有一定相關性，地理分布越近，遺傳距離越小，遺傳一致度越高，親緣關系越近。

2.4 UPMGA聚類分析

基于UPMGA對12個群體進行聚類（圖2），遺傳距離0.15處，12個群體分為4個組群，分別為組群Ⅰ、組群Ⅱ、組群Ⅲ、組群Ⅳ。組群Ⅰ和組群Ⅱ為國外群體，組群Ⅲ和組群Ⅳ為國內群體。遺傳距離為0.18處南美洲群體獨立分為一支，與其余國內外群體遺傳關系均較遠，從地理上看，它與別的群體相距均較遠，各種環境因素、氣候條件的不同可能是造成該結果的主要原因。遺傳距離為0.15處，中國山西群體獨立分為一支，與其他國內群體相差較遠。國內群體中，中國河北群體、黑龍江群體、遼寧群體遺傳關系較近，值得注意的是，河北群體和黑龍江群體盡管地理分布相距較遠，但兩者的遺傳距離最小，且聚在同一個亞群。國外群體中，亞洲群體、歐洲群體和北美洲群體遺傳關系較近，大洋洲群體次之。

表4 12個群體遺傳距離和遺傳一致度

將國外群體按不同的國家分為24個群體，然后與國內6個群體進行聚類（圖3），結果顯示，遺傳距離為0.20處，國內群體與國外群體有明顯的分離，與以上結果一致。但是，從圖中可以看出，遺傳距離為0.10處，國外群體被分為若干個小群，幾乎沒有分為大的聚類，總體而言，亞洲與歐洲各群體的遺傳關系較近，與12個群體間聚類的結果吻合。

基于UPGMA對146份個體進行聚類（圖4），結果顯示國內外資源劃分非常顯著。在遺傳距離為0.63處，146份黍稷資源分為3大組群，組群Ⅰ和組群Ⅱ為國外資源，組群Ⅲ為國內資源。組群Ⅰ包括2份資源，均為亞洲資源。在遺傳距離為0.39處，組群Ⅱ可分為3個亞群，分別為亞群A、亞群B和亞群C。亞群A包括29份資源，分別為亞洲6份、歐洲15份、大洋洲3份、北美洲3份和非洲2份。亞群B包括24份資源，分別為亞洲13份、歐洲5份、南美洲2份、大洋洲2份、非洲1份和北美洲1份。亞群C包括6份資源，分別為亞洲4份、歐洲1份和北美洲1份。在遺傳距離為0.45處，組群Ⅲ可分為5個亞群，分別為亞群D、亞群E、亞群F、亞群G和亞群H。亞群D包括4份資源，中國山西1份和中國河北3份。亞群E包括20份資源，分別為中國遼寧2份和中國黑龍江18份。亞群F包括10份資源，分別為中國河北2份、中國遼寧2份和中國黑龍江5份。亞群G包括21份資源，分別為中國寧夏15份、中國內蒙古5份和中國河北1份。亞群H包括30份資源，分別為中國河北20份、中國山西2份、中國內蒙古6份和中國遼寧2份。從各組群和亞群分布可以看出，亞洲與歐洲資源、中國河北與中國山西、中國內蒙古資源的遺傳關系較近，各組群分布與地理來源有一定關聯。

2.5 遺傳結構分析

利用Structure2.3.4軟件對146份黍稷資源進行群體結構分析（圖5），并根據對數似然方差Var[Lnp](D)]對值繪制折線圖，最佳分組應為9組。國內外群體間存在明顯的遺傳分化（表5），組群2、組群5、組群6、組群7和組群9是國內野生資源特有基因型，各省份群體分布較為分散，其中中國河北、中國黑龍江、中國遼寧資源集中分布在組群5、組群7和組群9，中國寧夏、中國內蒙古資源集中分布在組群2和組群6，各組群分布與地理來源基本相符。組群1和組群4為國外材料特有基因型，各群體分布更為集中。組群3中除2份亞洲資源外，其余均為國內野生資源（80%）組群8中除1份中國山西資源外，其余均為國外資源（97%），中國山西野生糜子和亞洲其他國家地方品種存在不同程度的遺傳混雜，遺傳關系相對較近。此外，中國寧夏資源全部分布在組群2，南美洲資源全部聚集在組群8，表明這兩個地區的黍稷資源的群體結構趨向單一化。相反，中國河北資源的群體結構最復雜，趨向多元化，分別被聚到6個組群中；其次是中國黑龍江和亞洲資源，分別被聚到4個組群中。這些結果很好地印證了UPGMA聚類結果，均于各群體的地理分布有一定的相關性。

表5 12個群體的遺傳結構分析

圖2 基于遺傳距離的12個群體黍稷資源聚類圖

圖3 基于遺傳距離的30個群體黍稷資源聚類圖

中國山西 Shanxi, China； 中國寧夏Ningxia, China； 中國內蒙古，Inner Mongolia, China； 中國河北 Hebei, China； 中國遼寧 Liaoning, China； 中國黑龍江 Heilongjiang, China； 亞洲Asia； 歐洲Europe； 南美洲South America； 北美洲North America； 大洋洲 Oceania； 非洲 Africa

圖5 146份黍稷種質資源群體遺傳結構圖

3 討論

3.1 黍稷SSR標記的遺傳多樣性

SSR標記是分析遺傳多樣的有效標記，已被廣泛應用于多種作物的遺傳多樣性研究，如玉米[11]、小麥[12]、燕麥[13]、水稻[26]等。近年來，不少報道中開發了黍稷特異性SSR標記，并分析了不同黍稷資源的遺傳多樣性。Cho等[14]開發了25對多態性SSR標記，共檢測到110個等位基因，平均每個位點4.4個。連帥等[25]采用63對SSR標記對192份黍稷地方品種和野生品種進行遺傳多樣性分析，共擴增出161個等位基因，平均每個位點2.56個，值為0.0877—0.8020，平均0.4855。Hou等[27]采用14個多態性SSR標記檢測到43個等位基因，平均每個位點3個。王瑞云等[24]利用15個糜子特異性熒光SSR標記分析了132份國內黍稷資源的遺傳多樣性，檢測到107個等位變異，平均每個位點變異數為7個，值范圍為0.0893—0.8538，平均0.4864。本研究中利用103對SSR標記對146份黍稷資源進行遺傳多樣性分析，擴增出308個等位基因，平均每個位點2.99個，值為0.064—0.8035，平均0.5544，Rp值為0.334—4.002，77.67%的標記具有適度分辨力。等位基因數低于王瑞云等[24]的7個，Cho等[14]的4.4個，高于連帥等[25]的2.56個，與Hou等[27]的3個相當，這可能與引物篩選標準不同有關。平均值大于0.5，高于連帥等[25]和王瑞云等[24]的報道，表明所選引物為高度多態性信息引物，可用于黍稷遺傳多樣性和起源進化的進一步研究。

3.2 參試材料的遺傳多樣性以及與地理分布的關系

本研究中MEGA結果顯示，個體間聚類和群體間聚類均將國內外群體分為兩大組群，國內群體中，河北、黑龍江、遼寧資源親緣關系較近，各遺傳參數也表明這些地區黍稷資源的遺傳多樣性豐富。國外群體中亞洲、歐洲遺傳距離最小，遺傳一致度最大，親緣關系最近。中國山西群體與南美洲群體的遺傳距離最大，兩者均獨立分為一支，與其他材料遺傳差異較大，且兩者地理分布相距遠，環境氣候有很大差異，表明聚類分析與地理分布相關。STRUCTURE分析結果顯示國內外群體間存在明顯的基因型分型，國內外群體間地理位置相距較遠，地理距離的隔離對各群體的基因流具有一定阻礙。此外還發現中國山西資源與國外亞洲資源存在部分遺傳混雜，表明兩者親緣關系較近，且中國山西與亞洲其他國家地理相距較近，暗示中國山西與亞洲其他國家的黍稷資源間可能存在基因流。黍稷群體間可能存在遺傳距離與地理距離的相關性，該研究結果與前人[10,21,23-24]報道相似。本研究聚類分析和遺傳結構分析結果均表明國外群體遺傳分化不明顯，沒有大的聚類，個體間存在相互滲透，這與連帥等[25]研究結果相似，這可能與國外黍稷資源遺傳背景比較狹窄有關，因此，國外黍稷在育種上需加強種質資源的引進與創新。

3.3 黍稷起源進化研究現狀

黍稷是干旱半干旱地區的主要作物之一，具有重要的經濟價值。研究黍稷野生資源是開展黍稷起源進化研究的重要根據之一。自1753年，提出黍稷原生于印度起[29]，國內外學者對黍稷的分類和進化進行了一系列的研究。目前，中國是世界上黍稷考古發現最早、資料最豐富的國家，以黃河中下游為中心，西至新疆，東至黑龍江，多處發現黍稷的遺跡，是公認的黍稷起源中心[28-30]。高俊山等[31]對中國620份栽培黍稷品種酯酶同工酶譜類型地理分布的研究表明中國的栽培黍稷起源于黃土高原。Hu等[21]和王星玉等[32]也得出了相同的結論。中國黍稷專家魏仰浩將野黍定位為黍稷的亞種[29-30]。通過從酯酶同工酶譜和生態性狀分析栽培黍稷與野生黍稷的關系，表明栽培黍稷是以具有各種原始性狀的野生黍稷及其近緣種進化而來[29-31]。本研究利用SSR標記，從DNA分子水平分析了黍稷野生資源的遺傳多樣性，結果顯示黍稷野生資源中中國河北群體的遺傳多樣性最豐富，其群體結構趨向多元化，與黑龍江、遼寧群體遺傳關系較近，與內蒙古、山西資源個體間存在相互滲透。與栽培黍稷相比較，野生黍稷的遺傳多樣性更為豐富，在很大程度上保持著原先的遺傳狀態。一方面可能是因為來自于河北的參試材料較多；另一方面，河北一帶野生糜子分布廣泛，目前栽培黍稷種植區域也較大，因此，推測河北也有可能是栽培黍稷的初級或次級起源中心，關于黍稷起源中心尚需豐富野生資源進一步研究確定。

4 結論

國內材料野生資源的遺傳多樣性高于國外材料，其中河北、黑龍江的遺傳多樣性最豐富，河北可能是黍稷的起源中心。聚類結果與遺傳結構分析結果相似，均與地理分布有一定的相關性。

附表1 參試材料

Table S1 Broomcorn millet varieties and landraces used in this study

序號Number群體編號與名稱Group number and name材料編號Material number材料來源Material source 11中國山西 Shanxi, China大同-1中國山西大同 Datong Shanxi, China 2中國山西 Shanxi, China大同-2中國山西大同 Datong Shanxi, China 3中國山西 Shanxi, China山西崗縣野糜子中國山西崗縣 Gangxian Shanxi, China 42中國寧夏 Ningxia, China寧夏-1中國寧夏 Ningxia, China 5中國寧夏 Ningxia, China寧夏-2中國寧夏 Ningxia, China 6中國寧夏 Ningxia, China寧夏-3中國寧夏 Ningxia, China 7中國寧夏 Ningxia, China寧夏-4中國寧夏 Ningxia, China 8中國寧夏 Ningxia, China寧夏-5中國寧夏 Ningxia, China 9中國寧夏 Ningxia, China寧夏-6中國寧夏 Ningxia, China 10中國寧夏 Ningxia, China寧夏-7中國寧夏 Ningxia, China 11中國寧夏 Ningxia, China寧夏-8中國寧夏 Ningxia, China 12中國寧夏 Ningxia, China寧夏-9中國寧夏 Ningxia, China 13中國寧夏 Ningxia, China寧夏-10中國寧夏 Ningxia, China 14中國寧夏 Ningxia, China寧夏-11中國寧夏 Ningxia, China 15中國寧夏 Ningxia, China寧夏-12中國寧夏 Ningxia, China 16中國寧夏 Ningxia, China寧夏-13中國寧夏 Ningxia, China 17中國寧夏 Ningxia, China寧夏-14中國寧夏 Ningxia, China 18中國寧夏 Ningxia, China寧夏-15中國寧夏 Ningxia, China 193 中國內蒙古 Inner Mongolia, China赤峰-2中國內蒙古赤峰 Chifeng Inner Mongolia, China 20中國內蒙古 Inner Mongolia, China赤峰-4中國內蒙古赤峰 Chifeng Inner Mongolia, China 21中國內蒙古 Inner Mongolia, China赤峰-5中國內蒙古赤峰 Chifeng Inner Mongolia, China 22中國內蒙古 Inner Mongolia, China赤峰-6中國內蒙古赤峰 Chifeng Inner Mongolia, China 23中國內蒙古 Inner Mongolia, China赤峰-7中國內蒙古赤峰 Chifeng Inner Mongolia, China 24中國內蒙古 Inner Mongolia, China赤峰-8中國內蒙古赤峰 Chifeng Inner Mongolia, China 25中國內蒙古 Inner Mongolia, China內蒙-2中國內蒙古 Inner Mongolia, China 26中國內蒙古 Inner Mongolia, China內蒙-3中國內蒙古 Inner Mongolia, China 27中國內蒙古 Inner Mongolia, China內蒙-4中國內蒙古 Inner Mongolia, China 28中國內蒙古 Inner Mongolia, China內蒙-5中國內蒙古 Inner Mongolia, China 29中國內蒙古 Inner Mongolia, China內蒙-6中國內蒙古 Inner Mongolia, China 304 中國河北 Hebei, China承德-2中國河北承德 Chengde Hebei, China 31中國河北 Hebei, China承德-3中國河北承德 Chengde Hebei, China 32中國河北 Hebei, China承德-5中國河北承德 Chengde Hebei, China 33中國河北 Hebei, China承德-A中國河北承德 Chengde Hebei, China 34中國河北 Hebei, China承德-6中國河北承德 Chengde Hebei, China 35中國河北 Hebei, China承德-7中國河北承德 Chengde Hebei, China 36中國河北 Hebei, China承德-8中國河北承德 Chengde Hebei, China 37中國河北 Hebei, China承德-9中國河北承德 Chengde Hebei, China 38中國河北 Hebei, China承德B中國河北承德 Chengde Hebei, China 39中國河北 Hebei, China承德-10中國河北承德 Chengde Hebei, China 40中國河北 Hebei, China承德-11中國河北承德 Chengde Hebei, China 41中國河北 Hebei, China承德C中國河北承德 Chengde Hebei, China 42中國河北 Hebei, China承德-12中國河北承德 Chengde Hebei, China 43中國河北 Hebei, China承德D中國河北承德 Chengde Hebei, China 44中國河北 Hebei, China河北-4中國河北 Hebei, China 45中國河北 Hebei, China河北-5中國河北 Hebei, China 46中國河北 Hebei, China河北-6中國河北 Hebei, China 47中國河北 Hebei, China河北-9中國河北 Hebei, China 48中國河北 Hebei, China河北-10中國河北 Hebei, China 49中國河北 Hebei, China河北-14中國河北 Hebei, China 50中國河北 Hebei, China河北-15中國河北 Hebei, China 51中國河北 Hebei, China河北-16中國河北 Hebei, China 52中國河北 Hebei, China河北-17中國河北 Hebei, China 53中國河北 Hebei, China河北-19中國河北 Hebei, China 54中國河北 Hebei, China河北-20中國河北 Hebei, China 55中國河北 Hebei, China河北-7中國河北 Hebei, China 565 中國遼寧 Liaoning, China遼寧-1中國遼寧 Liaoning, China 57中國遼寧 Liaoning, China遼寧-2中國遼寧 Liaoning, China 58中國遼寧 Liaoning, China遼寧-3中國遼寧 Liaoning, China 59中國遼寧 Liaoning, China遼寧-6中國遼寧 Liaoning, China 60中國遼寧 Liaoning, China遼寧-7中國遼寧 Liaoning, China 61中國遼寧 Liaoning, China遼寧-8中國遼寧 Liaoning, China 62中國遼寧 Liaoning, China遼寧-9中國遼寧 Liaoning, China 636 中國黑龍江 Heilongjiang, China齊齊哈爾野-1中國黑龍江齊齊哈爾 Qiqihar Heilongjiang, China 64中國黑龍江 Heilongjiang, China齊齊哈爾野-2中國黑龍江齊齊哈爾 Qiqihar Heilongjiang, China 65中國黑龍江 Heilongjiang, China齊齊哈爾野-3中國黑龍江齊齊哈爾 Qiqihar Heilongjiang, China 66中國黑龍江 Heilongjiang, China齊齊哈爾野-4中國黑龍江齊齊哈爾 Qiqihar Heilongjiang, China 67中國黑龍江 Heilongjiang, China齊齊哈爾野-5中國黑龍江齊齊哈爾 Qiqihar Heilongjiang, China 68中國黑龍江 Heilongjiang, China齊齊哈爾野-6中國黑龍江齊齊哈爾 Qiqihar Heilongjiang, China 69中國黑龍江 Heilongjiang, China齊齊哈爾野-7中國黑龍江齊齊哈爾 Qiqihar Heilongjiang, China 70中國黑龍江 Heilongjiang, China黑龍江245-1中國黑龍江 Heilongjiang, China 71中國黑龍江 Heilongjiang, China黑龍江245-2中國黑龍江 Heilongjiang, China 72中國黑龍江 Heilongjiang, China黑龍江245-3中國黑龍江 Heilongjiang, China 73中國黑龍江 Heilongjiang, China黑龍江245-4中國黑龍江 Heilongjiang, China 74中國黑龍江 Heilongjiang, China黑龍江245-5中國黑龍江 Heilongjiang, China 75中國黑龍江 Heilongjiang, China黑龍江245-6中國黑龍江 Heilongjiang, China 76中國黑龍江 Heilongjiang, China黑龍江245-7中國黑龍江 Heilongjiang, China 77中國黑龍江 Heilongjiang, China黑龍江246-2中國黑龍江 Heilongjiang, China 78中國黑龍江 Heilongjiang, China黑龍江246-3中國黑龍江 Heilongjiang, China 79中國黑龍江 Heilongjiang, China黑龍江246-4中國黑龍江 Heilongjiang, China 80中國黑龍江 Heilongjiang, China黑龍江246-5中國黑龍江 Heilongjiang, China 81中國黑龍江 Heilongjiang, China黑龍江246-6中國黑龍江 Heilongjiang, China 82中國黑龍江 Heilongjiang, China黑龍江246-8中國黑龍江 Heilongjiang, China 83中國黑龍江 Heilongjiang, China黑龍江247-1中國黑龍江 Heilongjiang, China 84中國黑龍江 Heilongjiang, China黑龍江247-2中國黑龍江 Heilongjiang, China 85中國黑龍江 Heilongjiang, China黑龍江248-1中國黑龍江 Heilongjiang, China 867 亞洲 AsiaAMES32315格魯吉亞 Georgia 87亞洲 AsiaAMES32316格魯吉亞 Georgia 88亞洲 AsiaPI163298印度 India 89亞洲 AsiaPI163300印度 India 90亞洲 AsiaPI170592土耳其 Turkey 91亞洲 AsiaPI207501阿富汗 Afghanistan 92亞洲 AsiaPI211058阿富汗 Afghanistan 93亞洲 AsiaPI211059阿富汗 Afghanistan 94亞洲 AsiaPI211060阿富汗 Afghanistan 95亞洲 AsiaPI212108阿富汗 Afghanistan 96亞洲 AsiaPI202317阿富汗 Afghanistan 97亞洲 AsiaPI253953阿富汗 Afghanistan 98亞洲 AsiaPI253955阿富汗 Afghanistan 99亞洲 AsiaPI649381阿富汗 Afghanistan 100亞洲 AsiaPI269955巴基斯坦 Pakistan 101亞洲 AsiaPI269957巴基斯坦 Pakistan 102亞洲 AsiaPI269958巴基斯坦 Pakistan 103亞洲 AsiaPI269959巴基斯坦 Pakistan 104亞洲 AsiaPI346939哈薩克斯坦 Kazakhstan 105亞洲 AsiaPI427247尼泊爾 Nepal 106亞洲 AsiaPI427247尼泊爾 Nepal 107亞洲 AsiaPI427248尼泊爾 Nepal 108亞洲 AsiaPI427249尼泊爾 Nepal 109亞洲 AsiaPI427250尼泊爾 Nepal 110亞洲 AsiaPI433381中國臺灣 Taiwan, China 1118 歐洲 EuropePI649372法國 France 112歐洲 EuropePI477123德國 Germany 113歐洲 EuropePI209790德國 Germany 114歐洲 EuropePI531407德國 Germany 115歐洲 EuropePI232929匈牙利 Hungary 116歐洲 EuropePI289320匈牙利 Hungary 117歐洲 EuropePI289321匈牙利 Hungary 118歐洲 EuropePI289324匈牙利 Hungary 119歐洲 EuropePI290726英國 England 120歐洲 EuropePI346934烏克蘭 Ukraine 121歐洲 EuropePI346941烏克蘭 Ukraine 122歐洲 EuropePI346942烏克蘭 Ukraine 123歐洲 EuropePI346933蘇聯 Soviet Union 124歐洲 EuropePI346935蘇聯 Soviet Union 125歐洲 EuropePI346944蘇聯 Soviet Union 126歐洲 EuropePI346945蘇聯 Soviet Union 127歐洲 EuropePI442533比利時 Belgium 128歐洲 EuropePI516181羅馬尼亞 Romania 129歐洲 EuropePI531399保加利亞 Bulgaria 130歐洲 EuropePI531402捷克斯洛伐克 Czechoslovakia 131歐洲 EuropePI531406捷克斯洛伐克 Czechoslovakia 1329 南美洲 South AmericaPI202294阿根廷 Argentina 133南美洲 South AmericaPI202295阿根廷 Argentina 13410 北美洲 North AmericaPI296376加拿大 Canada 135北美洲 North AmericaPI677103美國 America 136北美洲 North AmericaPI649385美國 America 137北美洲 North AmericaPI649384美國 America 138北美洲 North AmericaPI649383美國 America 13911 大洋洲 OceaniaPI365040澳大利亞 Australia 140大洋洲 OceaniaPI365842澳大利亞 Australia 141大洋洲 OceaniaPI365845澳大利亞 Australia 142大洋洲 OceaniaPI367683澳大利亞 Australia 143大洋洲 OceaniaPI367684澳大利亞 Australia 14412 非洲 AfricaPI517016摩洛哥 Morocco 145非洲 AfricaPI517017摩洛哥 Morocco 146非洲 AfricaPI531419肯尼亞 Kenya

附表2 所用SSR引物的序列及退火溫度

Supplement table 2 Sequence of SSR primers and annealing temperature

位點 Locus正向引物 Forward primer (5'-3')反向引物 Reverse primer (5'-3')退火溫度 Tm (℃) BM5CCCTTCTTCCCTACTTTGCCGTGTGCGTGCATGGGTGT55 BM114ATCGTAGAAACCATTGGCCCTGACCCATGGACACTTTTCA55 BM136AATGTCACAGGTTTCCCTCGGCGAGAAAGAGGAGAGGGTT55 BM289TGGGACAATATGGCAAGGTTACAAATGCCTGATGGTAGGC55 BM295CACACAGATATTTGGCACCGTGAGGATCCGAAAAGATTGG55 BM306ATTTTCTGGGCAATTCAACGGTCCTCATCCCTTCCCTCTC55 BM344AGCACTGAGGCACAATTCCTGTGCTGGGGTTTGTGACTTT55 BM374CTACCGCTTCAAAACGAAGGTGTCCCACTCTCCTACCTACTACC55 BM378ATGGGATGCACAGGTACACATCCTTAGGTCATCGTCCTATTTG55 BM396TTGATTATGCTTTGGAGGGGCCTCTCCTTACACGGGGATT55 BM484GGACAGGATGAGGAGGATCACGGTTACCATCGCCTTCTTA55 BM552GCAGCAATAGATTTGCCTCCTCTAGCCTACCCGAACTTGC55 BM637CACACTTGTGCTGTTGGGTCTCTAGCCTACCCCAACTTGC55 BM787TCCGTGACATGATGCAGACTGCACTATTTCTGCAGGAGCC55 BM994AACAGAACCGCACATCAGTGTCAGGTGAGCGATTCTTCCT55 BM1303TTGACGCGTACAACAAGAGCACCGAAGAGAAGGTTGCTGA55 BM1332CGGTTCAGTATTCAGGGCATCACTACACGTACCCGGTCCT55 BM1419CGGTTCAGTATTCAGGGCATCACAAGGATCAGGCCAAAGT55 BM1477GCGAAAGATGAACTGCCAATGCTGCGACGGATACTGATCT55 BM1533CTCTTGTCGTCTTGGTCCGTCGTGCGTGTGTCGAGAGAG55 BM1745CAGCTGATCATTTGAAAAGTATGGTAGTTGTCGCGATCATGAGG55 BM4724TGCATGAAGAAGTGCAGGTGATGATGCCCAAGAGATCGAC55 BM4739GGATGAGAGGAAGGGGTTGTCTTCCCTCCAAACTTCCTCC55 BM4741GTCTGCGTACCTCGGAAGTCTTCCCTCCAAACTTCCTCCT55 BM4749TCTGCGTACCTCGGAAGTCTTTCCCTCCAAACTTCCTCCT55 BM4761GATCACAAAGGGGAACAGGAACATGAGCCAGCTCTCCACT55 BM4830GTTGAAGGCAGACACCTTCCCTCCTCCCTCTCCCTCCC55 BM4846TCTAGCCTACCCCAACTTGCACAAAACCGCTGATCTGCTT55 BM4847GGACGAGTCGGTGAAGAGAGAATTTCTGTAGCGGTGTCGC55 BM4851TTCAATCCAATCCAAGGAGCAACCTTGAAGACATGAGCTGG55 BM4858CCCAACTTGCTTGGGATAAAACCTGATGGCATATTGCTCC55 BM4871CCAGTCGGTGAAGAGAGACCCTCTGTACTCGAGGGCGGT55 BM4877CGCTCCTTTCCTTTCTGTTGAGGAGCCTACTGCTGGAACA55 BM4882TGGGTGTCAGCTTGACTCAGTGCTTACATCCACTGCTTGC55 BM4947GATTTGGGACGAGTCGGTTAAAAGAATTTTCGGGCGACTT55 BM4954GGCGGATCGATGATTGTAGTTAAACGAGGGCAGAAAGGAA55 BM4956GGCTGTAAATCCGCTCGTAGTCCCTCTAATGCAGCCAGTT55 BM4958TCCATAGCTTCCATAGCATTTTCCCCAAAGGATATGGTTCCAA55 BM4962TCCATTGCTAGCGTTCTCCTTAGCGCCAACATGGTGTAGA55 BM4965CCTCCATCGAGAAACCAACTTGCAGAGACATGCATCACAA55 BM4967CTTCACCCACCACACCTTCTAGAATCTGACACCTGCACCA55 BM4969CTTCACCCACCACACCTTCTGAATCTGACACCTGCACCAA55 BM4973TTATCTTCCTCCGCCTCATCTGCATCACAACTGAAGACGA55 BM4997TTCTATCGCTGTGGTGGTCACAACTGCTTCTACAACGGCA55 BM5037ATTCCACCCTTGACCTCCTTTGCATTGACTGCCTTGAAAC55 BM5038GGTCTCTTAAGGAGGGCTGGTGCATTGACTGCCTTGAAAC55 BM5043TGACGAGATCATCTGGCTTGGAAAAGACAGGAGCGCAAAC55 BM5181CCTCACGGTGCATCAGAGTCCAGAGAGCTAAGATATAGGGAATC55 BM5190ACGATCACGAATAAATCCGCTTGAGATGATTCCTGTGTGTGTC55 BM5191GCCCGATAGCCAACGTAATAATGTTTTGAAGGCCCTATGC55 BM5197CTCTGCAACACTGCCACCTATTTACACCGTCGAATTTTATGA55 BM5198CGGAAATAAAGCAACCATGCATTGTGTGCTAGGCCACCTC55 BM5199TCACGTCCCTACACTTCACGGAGGTCGACAGGAAATGCTC55 BM5200CGGAAATAAAGCAACCATGCGACTGACGTAGGGCGTTCTC55 BM5222TGGAGTGTGTGTGTATGTGTGTGCTTGTAGGGCTTGTGGCAAT55 BM5233CTGCTAGAAAAGGACGGGTGCTAGCTAGCCAGGCCCATC55 BM5236CATGCCTCATCCCCTTATTCCGAAACGGCAAACACAAATA55 BM5255GTAAATTCAGATGCTTGCCGCCATCCGAGTTCAAGTTGGT55 BM5259CATACTTCCACACGTGCCTGGGACTATTCTTTGCTACTCCGGT55 BM5260GCGTGCGATAATCATTCGTGCCTACAGAAATATCAGCTGCG55 BM5278CGCCGTTTCATTTTTCTTTCTCACACACACACTACACTCACACTC55 BM5279CGCCGTTTCATTTTCTTTCTCTCGCTCACACACACACTCA55 BM5292GCATGCATGTGTGTGTGTGTAACATATGTGGACGGTTGGC55 LMX258GTGTCTCTTTCGTCTTGCCCGGGACACTTCCACCATCATC55 LMX503GTGGTACAGCTGCTCGTTCAAGGAGGAACCAGGAAGCAAT55 LMX515CGTTTTCTCGCTACACACGATGGACAACGGAAAACGTACA55 LMX621ATGAATCACCCGATCCACATACGCCAACATCAGCATATCA55 LMX632GCTGTCGGTCAGTCCTGTTTACGCCAACATCAGCATATCA55 LMX635GCACACGCATCATCACAAGTGCTCATTCAACGACAGATGC55 LMX645GCTAGGTGTCGTTTTCTCGCGGGGCGTCTAAGCTTGATATT55 LMX653CGATGAACGAAAATTCACCCGTTCATTCGTCCAAATGCCT55 LMX684GGCCTCCTGTGGTCTTCTCTGACCCTCTCCTCCTACCACC55 LMX786CCTGGACACACACACACACATCTTGTCACTGTCGGCGTAG55 LMX691ACTCATGGTTACGGCAACTGGCGCGAGAGAGAGAGAGAGA55 LMX836GCGCAGTAATATATTTCAGTAATTCAGCATCATCGTCAAGACCTCA55 LMX1067ATCGACGACTAGGCCCTGTATGCGGAGTGTCTTTGTTCTG55 LMX1080ATCGACGACTAGGCCCTGTAGGCCGTCACTATATCTGTCACC55 LMX1220AGTTACTTTTCCACACCCGCTTTATTGGGAGTTCGGATGG55 LMX1251GCAGTTTTGTTTCGATCGGTTGGTGTCAAGACATGTGGATG55 LMX1380GCCTCCTGTCTTGTAGCGTCAGGGTAGGCTGAGAGCCTGT55 LMX1387TTTCAGGGACTGGACTGGACGTAGGGGGTAGCTGAGAGCC55 LMX1400AGCAACGGAGGTGAGAGAGATCGACACACACGACACACAC55 LMX1429AATATCCCTTTTGTCGCACGATGCATTGATGGGCTTGATT55 LMX1625ACCCAACCGTATATCCAACGTGTCACAGTTGTCCTGGCAT55 LMX1672ATTTGACCTGTGACCTCGCTCCTTTCTGTTTCTGCAAGCC55 LMX1760TGTGGGAGAGAAGTGGGCCAAGGAAGGAATAAACCGCA55 LMX1749CTCAACACCACACCTTCCCTCTTTTGATCCATGGGCGTAG55 LMX1761GAGATGGTGCGGATTCTGAGTCATTTCCACTGTCACTGCC55 LMX1959GTGGTGAGAAGAAGCCTTGCGATGGGTCGATCAGCAGAAT55 LMX2175GTGACCGACGACTACCGTTTCACCCCTTTCCTCTCGCTAT55 LMX2410CTCAACACCACACCTTCCCTCTTTTGATCCATGGGCGTAG55 LMX1940GCAGTGGGTCAGCTTATGGTTCTCTCTGTGTGTGTGCGTG55 LMX2019CCTCTCCTTACACGGGGATTTTGATTATGCTTTGGAGGGG55 LMX2068TGAGATTGGCATCAAGCAAGTTTCTGGTCAGTTCGGTCAG55 LMX2074GCCACACTAAATAAGCTTTGTGTCTGGTCGTCACTGATTACGGA55 LMX2185GAGTTAGAGGACAGCGTGGCTGCAGCAGAGAATGTGCTACT55 LMX2281TTCTCCGTCAGCTCACATTGTCCATTGTTCATTTAGTAGAAACCT55 LMX2305CGCACTAGCCCTTGTCTTTCCGCCCTACGAACAAATCACT55 LMX2370CCAAGTCCCAAGGTTTTTCAGAGAGAGAGAGGCGTGTGCT55 LMX2382CAACAAGGTTGGTTGGCTTTATGCTGCTGCAGATGTTTTG55 LMX2540TATTCGAGCCCCATTTCTTGGCGTTATCCGGATGATGAAG55 LMX2734CACACAGATATTTGGCACCGTGAGGATCCGAAAAGATTGG55 LMX2782GCCGGAGTATAGATCCGACAGTCAGGCCGTGAACGTTATT55

[1] WANG R Y, WANG H G, LIU X Y, LIAN S, CHEN L, QIAO Z J, MCINERNEY C E, WANG L. Drought-induced transcription of resistant and sensitive common millet varieties.2017, 27(4): 1303-1314.

[2] 王瑞云, 楊陽, 王海崗, 陳凌, 王綸, 陸平, 劉敏軒, 喬治軍. 糜子PmNCED1的克隆及其對 PEG 脅迫的響應. 核農學報, 2018, 32(2): 244-256.

WANG R Y, YANG Y, WANG H G, CHEN L, WANG L, LU P, LIU M X, QIAO Z J. Cloning of PmNCED1 and its response to PEG stress., 2018, 32(2): 244-256. (in Chinese)

[3] 王星玉, 王綸. 黍稷種質資源描述規范和黍稷標準. 北京: 中國農業出版社, 2006: 5-15.

WANG X Y, WANG L.. Beijing: China Agriculture Press, 2006: 5-15. (in Chinese)

[4] 王璐琳, 王瑞云, 何杰麗, 薛延桃, 陳凌, 王海崗, 喬治軍. 糜子特異性SSR標記的開發. 山西農業科學, 2018, 46(1):1-4,86

WANG L L, WANG R Y, HE J L, XUE Y T, CHEN L, WANG H G, QIAO Z J. Development of broomcorn millet specific SSR marker., 2018, 2018, 46(1):1-4,86. ( in Chinese)

[5] WANG R Y, HUNT H V, QIAO Z J, Wang L, HAN Y H. Diversity and cultivation of broomcorn millet (L.) in China: A review., 2016, 70(3): 332-342.

[6] WANG R Y, WANG H G, LIU X H, JI X, CHEN L, LU P, LIU M X, TENG B, QIAO Z J. Waxy allelic diversity in common millet (L.) in China. The Crop Journal, 2018, https://doi.org/10.1016/j.cj.2018.02.004.

[7] 王瑞云. 糜子遺傳多樣性及進化研究進展. 北京: 中國農業出版社, 2017.

WANG R Y.(L.). Beijing: China Agriculture Press, 2017 (in Chinese)

[8] 胡興雨, 陸平, 賀建波, 王綸, 王星玉, 張紅生, 張宗文, 吳斌. 黍稷農藝性狀的主成分分析與聚類分析. 植物遺傳資源學報, 2008, 9(4): 492-496.

HU X Y, LU P, HE J B, WANG L, WANG X Y, ZHANG H S, ZHANG Z W, WU B. Principal component analysis and cluster analysis of agronomic traits of broomcorn millet., 2008, 9(4): 492-496. (in Chinese)

[9] 董孔軍, 楊天育, 何繼紅. 甘肅省糜子地方資源農藝性狀遺傳多樣性分析. 河北農業科學, 2012, 16(2): 1-3.

DONG K J, YANG T Y, HE J H. Genetic diversity of agronomic traits in broomcorn millet local resources in Gansu Province., 2012, 16(2): 1-3. (in Chinese)

[10] RAJPUT S G, SANTRA D K. Evaluation of genetic diversity of proso millet germplasm available in the United States using simple-sequence repeat markers., 2016, 56(5): 2401-2409.

[11] LI L J, YANG K C, PAN G T, RONG T Z. Genetic diversity of maize populations developed by two kinds of recurrent selection methods investigated with SSR markers., 2008, 7(9): 1037-1045.

[12] MANIFESTO M M, SCHLATTER A R, HOPP H E, Suárez E Y, Dubcovsky J. Quantitative evaluation of genetic diversity in wheat germplasm using molecular markers., 2001, 41(3): 682-690.

[13] He X Y, Bj?rnstad ?. Diversity of North European oat analyzed by SSR, AFLP and DArT markers., 2012, 125(1): 57-70.

[14] CHo Y I, Chung J W, Lee G A, Ma K H, Dixit A, Gwag J G, Park Y J. Development and characterization of twenty-five new polymorphic microsatellite markers in proso millet (L.)., 2010, 32(3): 267-273.

[15] 王銀月, 劉敏軒, 陸平, 喬治軍, 楊天育, 李海, 崔喜艷. 構建黍稷分子遺傳圖譜SSR引物的篩選. 作物雜志, 2014(4): 32-38.

Wang Y Y, Liu M X, Lu P, QIAO Z J, Yang T y, Li H, Cui X y. Construction of molecular genetic map of broomcorn millet for SSR primer screening., 2014(4): 32-38. (in Chinese)

[16] Rajput S G, Plyler-Harveson T, Santra D K. Development and characterization of SSR markers in proso millet based on switchgrass genomics., 2014, 5(1): 175.

[17] 王瑞云, 劉笑瑜, 王海崗, 陸平, 劉敏軒, 陳凌, 喬治軍. 用高基元微衛星標記分析中國糜子遺傳多樣性. 中國農業科學, 2017, 50(20): 3848-3859.

Wang R y, Liu X y, Wang H g, Lu P, Liu M x, Chen L, qiao z j. Analysis of genetic diversity of Chinese broomcorn millet by using high elementary microsatellite markers., 2017, 50(20): 3848-3859. (in Chinese)

[18] 朱建楚. 基于PCR標記的糜子遺傳多樣性分析[D]. 楊凌: 西北農林科技大學, 2005.

Zhu J C. Analysis of genetic diversity of broomcorn millet based on PCR markers[D]. Yangling:Agriculture and Forestry , 2005. (in Chinese)

[19] 董俊麗, 王海崗, 陳凌, 王君杰, 曹曉寧, 王綸, 喬治軍. 糜子骨干種質遺傳多樣性和遺傳結構分析. 中國農業科學, 2015, 48(16): 3121-3131.

Dong J l, Wang H g, Chen L, Wang J j, Cao X n, Wang L, qiao z j. Analysis of genetic diversity and genetic structure of broomcorn millet of the Chinese., 2015, 48(16): 3121-3131. (in Chinese)

[20] Hu Y G, Zhu J, Liu F, Zhang Z, Chai Y, Weining S. Genetic diversity among Chinese landraces and cultivars of broomcorn millet () revealed by the polymerase chain reaction., 2008, 153(3): 357-364.

[21] Hu X, Wang J, Lu P, Zhang H. Assessment of genetic diversity in broomcorn millet (L.) using SSR markers., 2009, 36(8): 491-500.

[22] 連帥, 王瑞云, 馬躍敏, 劉笑瑜, 季煦. 不同生態區糜子種質資源的遺傳多樣性分析. 山西農業大學學報: 自然科學版, 2015, 35(3): 225-231.

Lian S, Wang R Y, Ma Y M, Liu X Y, Ji X. Analysis of genetic diversity of broomcorn millet germplasms in different ecotype zone.,2015, 35(3): 225-231. (in Chinese)

[23] Liu M x, Xu Y, He J h, Zhang S, Wang Y y, Lu P. Genetic diversity and population structure of broomcorn millet (L.) cultivars and landraces in china based on microsatellite markers., 2016, 17(3): 370.

[24] 王瑞云, 季煦, 陸平, 劉敏軒, 許月, 王綸, 王海崗, 喬治軍. 利用熒光SSR分析中國糜子遺傳多樣性. 作物學報, 2017, 43(4): 530-548.

Wang R y, Ji x, Lu P, Liu M x, Xu Y, Wang L, Wang H g, qiao z j. Analysis of genetic diversity of Chinese broomcorn millet by fluorescence SSR., 2017, 43(4): 530-548. (in Chinese)

[25] 連帥, 陸平, 喬治軍, 張琦, 張茜, 劉敏軒, 王瑞云. 利用SSR分子標記研究國內外黍稷地方品種和野生資源的遺傳多樣性. 中國農業科學, 2016, 49(17): 3264-3275.

LIAN S, LU P, QIAO Z J, ZHANG Q, ZHANG Q, LIU M X, WANG R Y. Genetic diversity in broomcorn millet (L.) from China and abroad by using SSR markers., 2016, 49(17): 3264-3275. (in Chinese)

[26] THOMSON M J, POLATO N R, PRASETIYONO J,TRIJATMIKO K R, SILITONGA T S, MCCOUCH S R. Genetic diversity of isolated populations of Indonesian landraces of rice (L.) collected in east Kalimantan on the island of Borneo, 2009, 2(1): 80-92.

[27] Hou S Y, Sun Z X, Li Y S, Wang Y, Ling H B, Xing G F, HAN Y H, Li H Y. Transcriptomic analysis, genic SSR development, and genetic diversity of proso millet (; Poaceae)., 2017, 5(7):1-11

[28] 王星玉. 黍稷史話. 種子通訊, 1985, 4: 21.

Wang X y. History of broomcorn millet., 1985, 4: 21. (in Chinese)

[29] 魏仰浩. 試論黍的起源. 農業考古, 1986, 2: 248-251.

Wei Y h. On the Origin of broomcorn millet., 1986, 2: 248-251. (in Chinese)

[30] 董玉琛, 鄭殿升. 中國作物及其野生近緣植物: 糧食作物卷. 北京:中國農業出版社, 2006.

Dong Y C, Zheng D s. Crops and Their Wild Relatives in China: Food Crops Rolls. Beijing: China Agriculture Press, 2006. (in Chinese)

[31] 高俊山, 徐新宇, 胡榮海, 魏仰浩. 糜黍酯酶同工酶的研究. 作物品種資源, 1990, 2: 2.

Gao J s, Xu X y, Hu R h, Wei Y h. Studies on esterase isoenzymes of broomcorn millet.s, 1990, 2: 2. (in Chinese)

[32] 王星玉, 王綸, 溫琪汾, 師穎. 山西是黍稷的起源和遺傳多樣性中心. 植物遺傳資源學報, 2009, 10(3): 465-470.

Wang X y, Wang L, Wen Q F, Shi Y. Shanxi is the origin and genetic diversity center of broomcorn millet., 2009, 10(3): 465-470. (in Chinese)

（責任編輯 李莉）

Genetic Diversity and genetic relationship of Broomcorn Millet (L.) germplasm based on SSR markers

XUE YanTao1,2, LU Ping2, QIAO ZhiJun3, LIU MinXuan2, WANG RuiYun1,3

(1College of Agriculture, Shanxi Agricultural University, Taigu 030801, Shanxi;2Institute of Crop Science, Chinese Academy of Agricultural Sciences, Beijing 100081;3Institute of Crop Germplasm Resources, Shanxi Academy of Agricultural Sciences/Key Laboratory of Crop Gene Resources and Germplasm Enhancement on Loess Plateau, Ministry of Agriculture/Shanxi Key Laboratory of Genetic Resources and Genetic Improvement of Minor Crops, Taiyuan 030031)

【Objective】The objective of this study is to analyze the genetic diversity and relationship of broomcorn millet landraces and wild materials by SSR markers, to provide available data for further evolutionary study of broomcorn millet.【Method】137 SSR primers are used to identify polymorphisms in six representatives which selected randomly from the total of accessions. A total of 103 primers produce clear and reproducible polymorphic fragments among the six accessions and then are used to amplify 146 broomcorn millet accessions. Genetic diversity and relationship between different individuals and populations is evaluated by analyzing genetic parameter, clustering, and genetic structure.【Result】 103 SSR markers detect a total of 308 alleles (Na) with an average of 2.99 for each SSR and the mean values of Shannon-Weaver index (I), Nei andwere 0.8478, 0.3642 and 0.5544, respectively. Their resolution range was 0.334-4.002 and more than 60% distribution at intervals of 1-4, indicated the moderate resolving power of these SSR. The observed number of alleles (2.9126), Shannon-Weaver index (0.8302), expected heterozygosity (0.5023), andvalue (0.5278) of broomcorn millet accessions in China were all higher than those in abroad, indicated more abundant genetic diversity in Chinese samples. The genetic distance of the 12 populations ranged from 0.0783 to 0.5762 with a mean of 0.2938. The genetic identity ranged from 0.5620 to 0.9247 with a mean of 0.75. We found that the genetic similarity had a certain correlation with geographical distribution. The closer geographical distribution, the smaller genetic distance, the higher genetic identity. Cluster analysis divided 12 populations into 4 groups at a genetic distance of 0.15. Among them, resources in South America and Shanxi were each independently divided into one group, which had a far-distance relationship with other resources. In the inter-individual clustering, the division of resources at home and abroad was very significant. At a genetic distance of 0.63, 146 broomcorn millet accessions could be divided into three groups. Group Ⅰand group Ⅱwere foreign accessions, and group Ⅲwas domestic accessions. Further, groupⅡwas divided into three subpopulations at a genetic distance of 0.39, and group Ⅲwas divided into five subpopulations at a genetic distance of 0.45. There had closer genetic relationship between Asia and European resources, as well as Hebei, Shanxi and Inner Mongolia in China resources than other populations. The result of genetic structure analysis showed that there is obvious genetic differentiation between the domestic populations and foreign populations. Five groups (Group 2, Group 5, Group 6, Group 7 and Group 9) were unique genotypes which owned by Chinese wild resources and distributed more widely, 2 groups (Group 1 and Group 4) were unique genotypes of foreign resources and have a relative narrow distribution. The population structures of Ningxia and South America tend to be independent, and the population structures of Hebei, Heilongjiang and Asia tend to be diversified. The UPGMA clustering results were consistent with the results of genetic structure analysis, and the genetic relationships were related to their geographical distribution.【Conclusion】The genetic diversity of wild accessions is higher than that of landraces, of which Hebei population has the most abundant genetic diversity, so we suppose Hebei province may be the sub-origin center of broomcorn millet.

L; wild broomcorn millet; foreign germplasm; SSR markers; genetic diversity; population structure

2018-05-07；

2018-06-05

國家自然科學基金（31271791）、山西省回國留學人員科研資助項目（2016-066）、國家現代農業產業技術體系建設專項（CARS-06-13.5-A16）、農業部谷子高粱產業體系（CARS-07-12[1].5-A1）、農業部作物種質資源保護項目（NB2012-2130135-25-06-1）

薛延桃，Tel：18234839761；E-mail：yantaoxue305@163.com。通信作者劉敏軒，Tel：010-62159962；E-mail：liuminxuan@caas.cn。通信作者王瑞云，Tel：15234420135；E-mail：wry925@126.com

10.3864/j.issn.0578-1752.2018.15.002