王瑞云 季 煦 陸 平 劉敏軒 許 月 王 綸王海崗 喬治軍,*

?

利用熒光SSR分析中國(guó)糜子遺傳多樣性

王瑞云1,2,*,**季 煦1,**陸 平3劉敏軒3許 月4王 綸2王海崗2喬治軍2,*

1山西農(nóng)業(yè)大學(xué)農(nóng)學(xué)院, 山西太谷 030801;2山西省農(nóng)業(yè)科學(xué)院農(nóng)作物品種資源研究所 / 農(nóng)業(yè)部黃土高原作物基因資源與種質(zhì)創(chuàng)制重點(diǎn)實(shí)驗(yàn)室 / 雜糧種質(zhì)資源發(fā)掘與遺傳改良山西省重點(diǎn)實(shí)驗(yàn)室, 山西太原030031;3中國(guó)農(nóng)業(yè)科學(xué)院作物科學(xué)研究所, 北京100081;4吉林大學(xué)生命科學(xué)學(xué)院, 吉林長(zhǎng)春130012

分析糜子種質(zhì)資源的遺傳多樣性, 有助于了解糜子起源與進(jìn)化, 可為糜子優(yōu)異種質(zhì)發(fā)掘及資源高效利用提供理論基礎(chǔ)。本研究利用15個(gè)糜子特異性熒光SSR標(biāo)記檢測(cè)來(lái)源于中國(guó)11個(gè)省(區(qū))的132份糜子種質(zhì)資源, 檢測(cè)到107個(gè)等位變異, 每個(gè)位點(diǎn)等位變異數(shù)為2~14個(gè), 平均7個(gè); 基因多樣性指數(shù)為0.0936~0.8676, 平均0.5298; 多態(tài)性信息含量為0.0893~0.8538, 平均0.4864。采用遺傳距離的聚類將試驗(yàn)材料分為4類, 類群I來(lái)自東北春糜子區(qū), 類群II來(lái)自黃土高原春、夏糜子區(qū), 類群III來(lái)自于北方春糜子區(qū), 類群IV來(lái)自北方春糜子區(qū)和黃土高原春、夏糜子區(qū)。分析模型的遺傳結(jié)構(gòu)表明, 中國(guó)糜子資源來(lái)自4個(gè)(東北地區(qū)、黃土高原、北方地區(qū)和西北地區(qū))基因庫(kù), 與基于遺傳距離的聚類結(jié)果基本一致, 均與材料的地理起源相關(guān)。糜子遺傳變異豐富, 主要存在于糜子材料間。該結(jié)果從分子水平上準(zhǔn)確揭示了中國(guó)糜子的遺傳多樣性。

糜子; 熒光SSR; 遺傳多樣性; 聚類分析; 遺傳結(jié)構(gòu)

糜子(), 禾本科黍?qū)僖荒晟静荼局参? 主要分布于歐亞大陸干旱、半干旱地區(qū), 在俄羅斯、中國(guó)、印度、中東和中歐大面積種植[1-2]。目前俄羅斯、中國(guó)、烏克蘭等國(guó)家種質(zhì)資源庫(kù)中至少保存有24 680余份糜子資源[3-5]。糜子抗旱、耐瘠、耐鹽堿, 需水量少, 作為備荒作物及畜禽飼料被廣泛栽培[6-7]; 因其含有中和酸的堿性物質(zhì)及抗性淀粉, 作為保健食品在治療乳糜瀉、預(yù)防糖尿病和心血管病等方面已得到廣泛運(yùn)用[8-9]。

中國(guó)國(guó)家資源庫(kù)有8700余份糜子資源, 黃土高原和內(nèi)蒙古高原是糜子主產(chǎn)區(qū), 年均栽培面積50~60萬(wàn)公頃。目前, 糜子遺傳多樣性研究主要集中在表型性狀[10-15]、品質(zhì)特性[4]、抗逆性[4, 16-19]、抗病蟲(chóng)性[4, 20]等方面。劉峰等[10]、董孔軍等[11]、胡興雨等[12-13]和董俊麗等[14]研究了中國(guó)8643份糜子生物學(xué)特性差異, 發(fā)現(xiàn)穗重、生育期等在糜子種質(zhì)間存在豐富的遺傳變異。王瑞云等[15]調(diào)查了糜子葉片性狀特征(葉脈和氣孔密度)的解剖學(xué)差異, 發(fā)現(xiàn)氣孔密度遺傳變異豐富。王綸等[4]分析了中國(guó)6020份糜子資源的營(yíng)養(yǎng)品質(zhì)(蛋白質(zhì)、脂肪和賴氨酸), 篩選出優(yōu)異品質(zhì)種質(zhì)342份。王綸等[16]和王海崗等[17]評(píng)價(jià)了山西省520份糜子資源的抗旱性, 篩選出I級(jí)高度抗旱種質(zhì)和抗旱性強(qiáng)的地方品種“黃糜子”, 從生理、農(nóng)藝指標(biāo)方面揭示了糜子的抗旱機(jī)制。王綸等[18]和劉敏軒等[19]鑒定了6534份糜子資源的耐鹽性, 篩選出22份耐鹽種質(zhì)和衛(wèi)大黃糜等3份耐鹽性較強(qiáng)的糜子, 從芽苗期指標(biāo)揭示了糜子的耐鹽機(jī)制。王綸等[20]對(duì)中國(guó)6526份糜子資源進(jìn)行了抗黑穗病鑒定, 篩選出11份高抗種質(zhì)和574份抗病種質(zhì), 部分材料可直接用于黑穗病高發(fā)區(qū)。

分子標(biāo)記是評(píng)估作物遺傳多樣性的有效工具, 其中, SSR標(biāo)記在染色體上分布廣泛、再現(xiàn)性好、多態(tài)性信息含量高, 在糜子遺傳多樣性研究中最為實(shí)用。Hu等[21]和Rajput等[22-23]用385個(gè)種間SSR標(biāo)記(來(lái)自柳枝稷、小麥、大麥、燕麥和水稻)分析了中國(guó)和美國(guó)糜子資源的遺傳多樣性。Hunt等[24]、董俊麗等[14]和Liu等[25]用90個(gè)糜子特異性SSR標(biāo)記研究了中國(guó)和歐亞大陸糜子的遺傳差異。常規(guī)SSR的聚丙烯酰胺凝膠電泳檢測(cè)由于不同批次的擴(kuò)增條件、電泳時(shí)間和染色條件不一致, 無(wú)法準(zhǔn)確識(shí)別等位變異, 基因分型結(jié)果準(zhǔn)確性差。自動(dòng)熒光SSR (正鏈引物5'端標(biāo)記有熒光染料)檢測(cè)系統(tǒng)能準(zhǔn)確讀出基因片段大小、定量擴(kuò)增產(chǎn)物、分析通量大、費(fèi)用低[26], 已經(jīng)在玉米、高粱等植物中運(yùn)用[27-28], 但在糜子上應(yīng)用很少[24]。本試驗(yàn)利用15個(gè)糜子特異性熒光SSR標(biāo)記研究糜子主產(chǎn)區(qū)的132份試材, 從分子水平上分析糜子資源的遺傳多樣性和群體結(jié)構(gòu), 旨在為糜子的種質(zhì)資源研究、品種遺傳改良、優(yōu)異種質(zhì)鑒定和基因挖掘提供分子依據(jù)。

1 材料與方法

1.1 試驗(yàn)材料

132份糜子資源(見(jiàn)附表1)來(lái)源于11個(gè)省(區(qū)), 一部分由中國(guó)農(nóng)業(yè)科學(xué)院作物科學(xué)研究所(60份)、山西省農(nóng)業(yè)科學(xué)院農(nóng)作物品種資源研究所(4份)和山西省農(nóng)業(yè)科學(xué)院高寒區(qū)作物研究所(11份)提供; 另一部分為本實(shí)驗(yàn)室近5年收集的當(dāng)?shù)仄贩N(57份)。根據(jù)糜子栽培生態(tài)區(qū)劃[21, 25, 29-30], 132份糜子資源所屬生態(tài)區(qū)包括東北春糜子區(qū)(17份), 黃土高原春、夏糜子區(qū)(48份), 北方春糜子區(qū)(57份: 包括西北地區(qū)31份、山西高寒區(qū)21份和蒙古高原5份), 北方夏糜子區(qū)(10份)。將材料種植在溫室的直徑10 cm塑料盆中。

1.2 DNA提取

取生長(zhǎng)15~20 d的幼苗葉片, 用改良CTAB法[31]提取糜子基因組DNA。從每份樣品取3~5株幼苗葉片混合提取(約2 g)。用1%瓊脂糖凝膠電泳檢測(cè)DNA質(zhì)量, 終濃度調(diào)至30 ng μL–1。

1.3 SSR引物、PCR擴(kuò)增和自動(dòng)熒光檢測(cè)

參照Cho等[32]的引物序列(見(jiàn)附表2), 由上海生工生物工程技術(shù)服務(wù)有限公司合成。25 μL PCR體系含2.5 μL 10×緩沖液、2 μL MgCl2(25 mmol L–1)、正向引物(1 mmol L–1) 0.5 μL [5￠端標(biāo)記熒光素FAM (藍(lán)色)、ROX (紅色)、HEX (綠色)或TAMRA (黃色)(見(jiàn)附圖1)]、反向引物(1 mmol L–1) 0.5 μL、0.5 μL dNTPs (10 mmol L–1)、0.2 μLDNA聚合酶(5 U L–1)、1 μL模板DNA和17.8 μL ddH2O。反應(yīng)程序?yàn)?5℃ 3 min; 95℃ 30 s, 60℃ 30 s, 72℃ 30 s, 10個(gè)循環(huán); 95℃ 30 s, 55℃ 30 s, 72℃ 30 s, 20個(gè)循環(huán); 72℃ 6 min。用15對(duì)引物(見(jiàn)附表2)擴(kuò)增糜子試材, 用ABI3730XL (美國(guó)應(yīng)用生物系統(tǒng)公司) DNA分析儀對(duì)PCR產(chǎn)物進(jìn)行自動(dòng)熒光檢測(cè)。

1.4 數(shù)據(jù)分析

用Genemapper 4.0 (美國(guó)應(yīng)用生物系統(tǒng)公司)軟件自動(dòng)生成每個(gè)位點(diǎn)的圖譜文件, 給出PCR擴(kuò)增片段長(zhǎng)度和峰值(熒光強(qiáng)度值), 以PDF和Microsoft Excel數(shù)據(jù)導(dǎo)出。在每個(gè)位點(diǎn)對(duì)所有材料的等位基因依據(jù)最大峰值打分。

每個(gè)標(biāo)記檢測(cè)一個(gè)位點(diǎn), 每一條多態(tài)性條帶視為一個(gè)等位基因。用PowerMarker 3.25[33]計(jì)算每對(duì)引物的多樣性參數(shù), 包括等位基因數(shù)(a)、等位基因頻率、基因多樣性指數(shù)()、多態(tài)性信息含量指數(shù)(PIC)和稀有等位基因數(shù)等。用PopGen1.32[34]計(jì)算不同樣品間的Nei’s遺傳距離。用MEGA 5.0[35]構(gòu)建Neighbour-Joining聚類圖。用Structure 2.2[36]鑒定132份糜子資源的遺傳類群。用GenALEx 6.5[37]分析分子變異、評(píng)價(jià)群組內(nèi)和群組間的遺傳分化。

2 結(jié)果與分析

2.1 SSR標(biāo)記的多態(tài)性位點(diǎn)擴(kuò)增分析

用15個(gè)標(biāo)記在132份糜子材料中擴(kuò)增出可賦分峰譜(見(jiàn)附圖1)。其中12個(gè)標(biāo)記(GB-PaM-004、GB-PaM-013、GB-PaM-014、GB-PaM-025、GB-PaM- 023、GB-PaM-061、GB-PaM-096、GB-PaM-107、GB-PaM-115、GB-PaM-121、GB-PaM-126、GB-PaM- 134)擴(kuò)增出一組峰譜, 為單位點(diǎn)擴(kuò)增, 3個(gè)標(biāo)記(GB-PaM-066、GB-PaM-094和GB-PaM-145)擴(kuò)增出2組峰譜, 為雙位點(diǎn)(用a和b表示)擴(kuò)增, 共擴(kuò)增出18個(gè)位點(diǎn)。其中標(biāo)記GB-PaM-094b在全部樣品中的擴(kuò)增片段相同, 不具多態(tài)性, 因此15個(gè)標(biāo)記共擴(kuò)增出17個(gè)多態(tài)性位點(diǎn)(表1)。

表1 15個(gè)SSR標(biāo)記的遺傳參數(shù)

主要等位基因和稀有等位基因分別為基因頻率大于或小于5%的等位基因。

Major allele or rare allele indicates its gene frequency was higher or less than 5%, respectively.

2.2 SSR標(biāo)記的遺傳多樣性分析

從表1可以看出, 17個(gè)位點(diǎn)共有107個(gè)等位變異; 每個(gè)位點(diǎn)等位變異數(shù)為2~14個(gè)(平均7個(gè)), 其中, GB-PaM-134最高, 其次是GB-PaM-023和GB-PaM-126, GB-PaM-013和GB-PaM-145最低。PIC為0.0893~0.8538, 平均0.4864; GB-PaM-126最高, GB-PaM-066b最低。8個(gè)位點(diǎn)的PIC均大于0.5, 為高度多態(tài)性位點(diǎn); 6個(gè)位點(diǎn)(GB-PaM-004、GB-PaM- 013、GB-PaM-014、GB-PaM-025、GB-PaM-061和GB-PaM-066a)的PIC為0.25~0.50, 為中度多態(tài)性位點(diǎn); 3個(gè)位點(diǎn)(GB-PaM-066b、GB-PaM-145a和GB- PaM-145b)的PIC小于0.25, 為低度多態(tài)性位點(diǎn)?；蚨鄻有灾笖?shù)為0.0936~0.8676, 平均0.5298。其中, GB-PaM-126最大, GB-PaM-066b最小。

2.3 基于遺傳距離的聚類分析

聚類分析(圖1)表明, 132份糜子材料聚為4個(gè)群組。群組I共6份材料, 青海(4份)、寧夏(1份)和內(nèi)蒙(1份), 均屬于北方春糜子區(qū)。群組II為31份材料, 黑龍江(2份)、吉林(2份), 河北(1份)、山東(3份), 甘肅(9份)、寧夏(7份), 山西(5份), 陜西(2份); 屬于黃土高原春、夏糜子區(qū)、東北春糜子區(qū)、北方春糜子區(qū)和北方夏糜子區(qū)。群組III包括8份材料, 山西(5份)、陜西(2份), 甘肅(1份); 屬于黃土高原春、夏糜子區(qū)和北方春糜子區(qū)。群組IV包括87份材料, 山西黃土高原生態(tài)區(qū)(36份)、陜西(1份), 黑龍江(10份)、吉林(1份)、遼寧(2份), 山西高寒區(qū)(18份)、內(nèi)蒙古(4份)、甘肅(5份)、寧夏(3份)、青海(1份), 河北(5份)、山東(1份), 主要來(lái)自黃土高原春、夏糜子區(qū)、北方春糜子區(qū)和東北春糜子區(qū)。從圖1還可以看出, 來(lái)自山西的地方品種有93%聚在群組IV中, 育成品種中有71.4%聚在群組III中。基于多態(tài)性SSR位點(diǎn)劃分的類群與地理來(lái)源關(guān)系密切, 說(shuō)明不同省區(qū)的糜子資源具有地域特色。

2.4 中國(guó)糜子種質(zhì)資源群體遺傳結(jié)構(gòu)

用Structure 2.2軟件評(píng)估中國(guó)不同省(區(qū))的糜子材料, 發(fā)現(xiàn)遺傳群體數(shù)在=2和=4處有明顯的峰(圖2), 因此在2種模式下分析糜子的遺傳結(jié)構(gòu)。=2時(shí), 劃分為2個(gè)類群(圖3)。類群I為紅色部分(80份), 主要來(lái)自北方春糜子區(qū)和東北春糜子區(qū), 包括山西(14份)和甘肅(13份), 其次是黑龍江(12份)和寧夏(10份); 類群II為綠色部分(52份), 主要來(lái)自黃土高原春、夏糜子區(qū), 山西最多(31份)(圖4)。從圖4可以看出, 紅色和綠色群組以長(zhǎng)城為界, 分別分布在長(zhǎng)城以北(北方春糜子區(qū)和東北春糜子區(qū))和長(zhǎng)城以南(黃土高原春、夏糜子區(qū))。=4時(shí), 劃分為4個(gè)類群(圖3)。類群I為紅色部分(25份), 主要來(lái)自東北春糜子區(qū), 黑龍江最多(12份), 代表東北糜子基因庫(kù)。類群II為綠色部分(37份), 主要來(lái)自黃土高原春、夏糜子區(qū), 山西最多(34份), 代表黃土高原基因庫(kù)。類群III為藍(lán)色部分(30份), 主要來(lái)自北方春糜子區(qū)。類群IV為黃色部分(39份), 主要來(lái)自北方春糜子區(qū)和黃土高原春、夏糜子區(qū), 山西最多(17份),甘肅(8份)和寧夏(7份)次之(圖4)。

紅色群組(=2)分化成紅色、藍(lán)色和黃色(=4)。21份山西糜子資源中10份(47.6%)為藍(lán)色, 10份(47.6%)為黃色; 13份甘肅糜子資源中3份(23.1%)為藍(lán)色, 8份(61.5%)為黃色。

綠色群組(=2)在=4模式下絕大部分仍為綠色, 小部分為黃色和藍(lán)色, 極少部分為紅色。43份山西資源中33份(76.7%)仍為綠色; 7份(16.3%)為黃色; 2份(4.7%)為藍(lán)色; 1份(2.3%)為紅色。

=2和=4結(jié)構(gòu)圖中各群組的遺傳多樣性衡量指標(biāo)見(jiàn)表2。=2時(shí), 紅色群組資源的多樣性比綠色群組豐富, 說(shuō)明長(zhǎng)城以北(北方春糜子區(qū)和東北春糜子區(qū))資源的遺傳多樣性高于長(zhǎng)城以南(黃土高原春、夏糜子區(qū))的資源。=4時(shí), 各個(gè)分類群組的遺傳多樣性參數(shù)值差別較小。

Delta根據(jù)Evanno等[38]的方法計(jì)算得到, 針對(duì)基因庫(kù)數(shù)目()建模。

Deltawas calculated by the method of Evanno et al.[38]and model was established by the number of genepools ().

2.5 分子遺傳變異(AMOVA)分析

對(duì)Structure遺傳結(jié)構(gòu)劃分類群的方差分析(表3)表明, 遺傳結(jié)構(gòu)分類群的遺傳變異極顯著(< 0.001)。其中, 65%發(fā)生在樣品間, 26%發(fā)生在地區(qū)內(nèi)樣品間, 說(shuō)明糜子種質(zhì)資源的遺傳變異主要存在于糜子資源間。

2.6 遺傳距離和遺傳結(jié)構(gòu)聚類分析的比較

比較遺傳距離聚類和遺傳結(jié)構(gòu)(=4)分析結(jié)果表明(表4), 內(nèi)蒙古、寧夏、青海和山東4個(gè)省(區(qū)) 2種聚類結(jié)果完全一致; 甘肅省除1份材料外, 分為3個(gè)群組; 黑龍江省除2份材料外, 分為1個(gè)群組; 吉林省除1份材料外, 分為2個(gè)群組; 遼寧省除1份材料外, 分為1個(gè)群組; 山西省除1份材料外, 分為3個(gè)群組; 陜西省除1份材料外, 分為2個(gè)群組?？梢?jiàn)基于模型的遺傳結(jié)構(gòu)聚類與基于遺傳距離聚類結(jié)果基本一致, 均與糜子資源的地理起源相關(guān)。

3 討論

3.1 基于SSR標(biāo)記的多樣性分析

分子標(biāo)記是揭示糜子遺傳多樣性的有效手段, 已經(jīng)有不少標(biāo)記, 如ITS、ISJ、AFLP、ISSR、RAPD、SSR、SNP等[6, 21, 25, 30, 39-44]用于分析糜子遺傳差異。其中, SSR以共顯性好、多態(tài)性豐富成為糜子最實(shí)用的檢測(cè)標(biāo)記。2009年, Hu等[21]首先用其他作物的SSR標(biāo)記研究中國(guó)糜子資源的遺傳多樣性, 由于存在種間差異, 引物通用性低, 開(kāi)發(fā)糜子特異性SSR并篩選出多態(tài)性高的核心引物, 對(duì)準(zhǔn)確評(píng)價(jià)糜子遺傳差異具有重要意義。Cho等[32]首次構(gòu)建了25個(gè)糜子微衛(wèi)星標(biāo)記, 隨后, 董俊麗等[14]、季煦等[45]和連帥等[46]用其中的部分標(biāo)記分析了中國(guó)糜子的遺傳多樣性。然而, 由于傳統(tǒng)SSR檢測(cè)所用的聚丙烯酰胺凝膠電泳以分子量Marker作對(duì)照, 用肉眼粗略估計(jì)擴(kuò)增片段大小, 難以區(qū)分幾個(gè)堿基對(duì)的差異, 檢測(cè)到的變異數(shù)低[14, 25, 45-46]。2011年, Hunt等[24]首次利用SSR熒光標(biāo)記分析歐亞大陸糜子資源的遺傳多樣性, 但僅包括少量中國(guó)材料(33份), 不足以準(zhǔn)確反映中國(guó)糜子的遺傳變異情況。本試驗(yàn)15個(gè)SSR標(biāo)記中, 12個(gè)與董俊麗等[14]所用相同, 但由于引物標(biāo)記了熒光, 能檢測(cè)到聚丙烯酰胺凝膠電泳識(shí)別不出的等位變異。本研究檢測(cè)到每對(duì)引物的平均等位基因數(shù)為7個(gè), 而前人研究結(jié)果均小于6個(gè)[14, 24-25, 32]。本研究每對(duì)引物檢測(cè)到的PIC值為0.4864, 而前人研究結(jié)果最高為0.3926[14, 24-25, 32]。黃河流域的黃土高原是中國(guó)糜子的起源中心[21], 本研究所用材料多數(shù)來(lái)源于該區(qū)域, 故檢測(cè)到豐富的遺傳多樣性。

橫坐標(biāo)的數(shù)字代表糜子材料序號(hào)。

Numbers in the horizontal axis represent serial number of the accession.

表2 遺傳結(jié)構(gòu)圖中各分類群的多樣性統(tǒng)計(jì)

表3 遺傳結(jié)構(gòu)圖中糜子資源的AMOVA分析

本研究糜子試材來(lái)自11個(gè)省(區(qū)), 以每個(gè)地區(qū)為一個(gè)群體, 借助SSR標(biāo)記在DNA水平上分析了糜子種質(zhì)資源的遺傳差異。其中, 3對(duì)引物(GB-PaM- 066、GB-PaM-094和GB-PaM-145)均擴(kuò)增出2個(gè)變異位點(diǎn), 且GB-PaM-066的2個(gè)變異位點(diǎn)均有多態(tài)性, GB-PaM-094僅一個(gè)位點(diǎn)具有多態(tài)性, 這與Hunt等[24]研究結(jié)果一致。本研究發(fā)現(xiàn)GB-PaM-145有2個(gè)多態(tài)性位點(diǎn), 而Hunt等[24]僅檢測(cè)到1個(gè)(位點(diǎn)b), 這可能與研究材料不同有關(guān), 本研究用的中國(guó)資源(132份)是Hunt等[24](33份)的4倍。

王銀月等[43]用高通量測(cè)序結(jié)合PCR擴(kuò)增篩選出糜子特異性SSR引物116對(duì)(雙堿基重復(fù)102對(duì)、三堿基重復(fù)14對(duì))。前人用到的引物多為雙堿基和三堿基重復(fù), 四堿基重復(fù)僅一個(gè)(EF117731)[14, 32]。今后研究需要開(kāi)發(fā)五堿基和六堿基重復(fù)及復(fù)合堿基SSR, 為糜子品質(zhì)和抗性優(yōu)異基因挖掘、遺傳多樣性、群體遺傳結(jié)構(gòu)及遺傳連鎖圖譜構(gòu)建等研究提供更多可靠的分子檢測(cè)工具。

3.2 糜子種質(zhì)資源遺傳結(jié)構(gòu)分析

132份糜子資源聚類結(jié)果中, 山西的地方品種93%聚在群組IV中, 育成品種71.4%聚在群組III中, 說(shuō)明中國(guó)不同省的糜子資源具有地域特色, 這很可能與人們的食物偏好和文化差異有關(guān)。分析發(fā)現(xiàn), 山西的地方品種在4個(gè)群組中均有呈現(xiàn), 反映出山西糜子資源類型比較豐富, 這與董俊麗等[14]的研究結(jié)果一致。山西的育成品種大多集中在群組III中, 說(shuō)明當(dāng)代育成品種遺傳基礎(chǔ)狹窄, 可能與種質(zhì)資源的利用不充分有關(guān)。建議當(dāng)代的糜子育種家要加強(qiáng)育種材料的交換, 擴(kuò)大育成品種的遺傳基礎(chǔ), 提高品種的適應(yīng)性和育種水平。

表4 基于遺傳距離和遺傳結(jié)構(gòu)(K=4)各分類群組在不同省份的糜子分布

本研究基于遺傳距離的聚類將132份糜子分為4個(gè)群組, 與基于模型的遺傳結(jié)構(gòu)(=4)分類結(jié)果相符。不同地理來(lái)源糜子的遺傳多樣性分析表明, 山西資源的多項(xiàng)遺傳多樣性衡量指標(biāo)(等位基因數(shù)、H和PIC)均最高, 說(shuō)明山西的糜子資源多樣性相對(duì)較豐富, 這與前人研究結(jié)果[14, 47]相符, 可能也與本研究中山西糜子材料數(shù)多有關(guān)。分析山西64份糜子資源的遺傳結(jié)構(gòu)表明大部分資源屬于黃土高原春、夏糜子區(qū), 其中, 山西資源在不同類群中均有呈現(xiàn), 這與以往研究結(jié)果相符[14, 25]。

3.3 糜子群體遺傳結(jié)構(gòu)類群間的遺傳分化

群體間存在顯著遺傳分化的物種分化系數(shù)大于0.15[48], Hunt等[35]研究歐亞大陸糜子資源的遺傳多樣性發(fā)現(xiàn),=2和=6時(shí)糜子材料的分化系數(shù)分別為0.162和0.324, 均高于本研究結(jié)果, 原因在于其材料來(lái)源廣泛, 地理差異明顯。本研究AMOVA分析發(fā)現(xiàn)=2和=4時(shí)糜子材料的分化系數(shù)分別為0.088和0.084, 說(shuō)明中國(guó)糜子結(jié)構(gòu)類群間遺傳分化不明顯, 這與遺傳結(jié)構(gòu)類群間的多樣性豐富程度類似相符。今后研究需進(jìn)一步豐富育種材料, 以取得更加可靠的結(jié)果。

4 結(jié)論

中國(guó)糜子資源遺傳多樣性豐富, 利用模型的遺傳結(jié)構(gòu)聚類和基于遺傳距離的聚類結(jié)果一致, 均可以把材料劃分為4個(gè)群組。

[1] Lu H, Zhang J, Liu K, Wu N, Li Y, Zhou K, Ye M, Zhang T, Zhang H, Yang X, Shen L, Xu D, Li Q. Earliest domestication of common millet () in East Asia extended to 10000 years ago., 2009, 106: 7367–7372

[2] Martin J H, Leonard W H, Stamp D L. Principles of field crop production, 3rd edn. New York: MacMillan Publishing Co., Inc., 1976. pp 415–429

[3] Kurzeva A, Romanova O, Krylov A. Genetic resources of common millet (L.) at the Vavilov Research Institute of Plant Industry (VIR). In:Chai Y, Feng B, eds. Advances in Broomcorn Millet Research: Proceedings of the 1st International Symposium on Broomcorn Millet. Yangling: Northwest A&F University Press, 2012. pp 76–81

[4] 王綸, 王星玉, 溫琪汾, 武變娥. 中國(guó)黍稷種質(zhì)資源研究與利用. 植物遺傳資源學(xué)報(bào), 2005, 6: 474–477 Wang L, Wang X Y, Wen Q F, Wu B E. Research and utilization of proso millet germplasm resource in China., 2005, 6: 474–477 (in Chinese with English abstract)

[5] Kobyzeva L N, Gorbacheva S N, Grigorashchenko L V, Gorlacheva O V, Biryukova O V. Formation of National collection and millet breeding in Ukraine. In:Chai Y, Feng B, eds. Advances in Broomcorn Millet Research: Proceedings of the 1st International Symposium on Broomcorn Millet. Yangling: Northwest A&F University Press, 2012. pp 67–73

[6] Lagler R, Gyulai G, Humphreys M, Szabo Z, Horvath L, Bittsanszky A. Morphological and molecular analysis of common millet () cultivars compared to a DNA sample from the 15th century (Hungary)., 2005, 146: 77–85

[7] Bagdi A, Balázs G, Schmidt J, Szatmári M, Schoenlechner R, Berghofer E, T?m?sk?zia S. Protein characterization and nutrient composition of Hungarian proso millet varieties and the effect of decortication., 2011, 40: 128–141

[8] Nishizawa N, Fudamo Y. The elevation of plasma concentration of high-density lipoprotein chloesterol in mice fed with protein from proso millet ()., 1995, 59: 333–335

[9] Park M Y, Jang H H, Kim J B, Yoon H N, Lee J Y, Lee Y M, Kim J H, Park D S. Hog millet (L.): supplemented diet ameliorates hyperlipidemia and hepatic lipid accumulation in C57BL/6J-ob/obmice., 2011, 5: 511–519

媒體報(bào)道稱，一項(xiàng)報(bào)告顯示，共享經(jīng)濟(jì)已經(jīng)擴(kuò)展至中國(guó)全國(guó)，2017年的交易額達(dá)到4．92萬(wàn)億元，比上年增長(zhǎng)47%。正因?yàn)榇?，大公司和中?guó)政府正在力推這個(gè)經(jīng)濟(jì)潮流，認(rèn)為它是推動(dòng)經(jīng)濟(jì)發(fā)展的力量。據(jù)預(yù)測(cè)，到2020年，共享經(jīng)濟(jì)規(guī)模將占到中國(guó)國(guó)內(nèi)生產(chǎn)總值（GDP）的10%，而在2015年這一數(shù)字僅為3%。

[10] 劉峰. 糜子農(nóng)藝性狀的遺傳多樣性分析. 黑龍江農(nóng)業(yè)科學(xué), 2010, (3): 15–16 Liu F. Genetic diversity analysis of agronomic traits on proso millet.,2010, (3): 15–16 (in Chinese with English abstract)

[11] 董孔軍, 楊天育, 何繼紅. 甘肅省糜子地方資源農(nóng)藝性狀遺傳多樣性分析. 河北農(nóng)業(yè)科學(xué), 2012, 16(2): 1–3 Dong K J, Yang T Y, He J H. Analysis of genetic diversity of agronomic traits of local broomcorn millet germplasm resources in Gansu province., 2012, 16(2): 1–3 (in Chinese with English abstract)

[12] 胡興雨, 陸平, 賀建波, 王綸, 王星玉, 張紅生, 張宗文, 吳斌. 黍稷農(nóng)藝性狀的主成分分析與聚類分析. 植物遺傳資源學(xué)報(bào), 2008, 9: 492–496 Hu X Y, Lu P, He J B, Wang L, Wang X Y, Zhang H S, Zhang Z W, Wu B. Principal components and cluster analysis of agronomic traits of proso millet ()., 2008, 9: 492–496 (in Chinese with English abstract)

[13] 胡興雨, 王綸, 張宗文, 陸平, 張紅生. 中國(guó)黍稷核心種質(zhì)的構(gòu)建. 中國(guó)農(nóng)業(yè)科學(xué), 2008, 41: 3489–3502 Hu X Y, Wang L, Zhang Z W, Lu P, Zhang H S. Establishment of broomcorn millet core collection in China., 2008, 41: 3489–3502 (in Chinese with English abstract)

[14] 董俊麗, 王海崗, 陳凌, 王君杰, 曹曉寧, 王綸, 喬治軍. 糜子骨干種質(zhì)資源遺傳多樣性和遺傳結(jié)構(gòu)分析. 中國(guó)農(nóng)業(yè)科學(xué), 2015, 48: 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 structure of proso millet core germplasm., 2015, 48: 3121–3131 (in Chinese with English abstract)

[15] 王瑞云, 季煦, 劉曉歡, 楊陽(yáng), 徐鳳杰, 程冀南. 不同生態(tài)區(qū)糜子種質(zhì)間葉片及光合特性的多樣性分析. 山西農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2014, 34(2): 97–102 Wang R Y, Ji X, Liu X H, Yang Y, Xu F J, Cheng J N. Diversity of leaf traits and photosynthetic characteristics in broomcorn millet (L.) germplasms among different ecotype zone of China.(Nat Sci Edn), 2014, 34(2): 97–102 (in Chinese with English abstract)

[16] 王綸, 溫琪汾, 曹利萍, 王星玉. 黍稷抗旱種質(zhì)篩選及抗旱機(jī)理研究. 山西農(nóng)業(yè)科學(xué), 2007, 35(4): 31–34 Wang L, Wen Q F, Cao L P, Wang X Y. Drought-resistant germplasm screening and drought-resistant mechanism in proso millet., 2007, 35(4): 31–34 (in Chinese with English abstract)

[18] 王綸, 王星玉, 溫琪汾, 武變娥, 曹利萍. 中國(guó)黍稷種質(zhì)資源耐鹽性鑒定. 植物遺傳資源學(xué)報(bào), 2007, 8: 426–429 Wang L, Wang X Y, Wen Q F, Wu B E, Cao L P. Identification of salt tolerance in Chinese proso millet germplasm., 2007, 8: 426–429 (in Chinese with English abstract)

[19] 劉敏軒, 張宗文, 吳斌, 陸平. 黍稷種質(zhì)資源芽、苗期耐中性混合鹽脅迫評(píng)價(jià)與耐鹽生理機(jī)制研究. 中國(guó)農(nóng)業(yè)科學(xué), 2012, 45: 3733–3743 Liu M X, Zhang Z W, Wu B, Lu P. Evaluation of mixed salt-tolerance at germination stage and seedling stage and the related physiological characteristics ofL., 2012, 45: 3733–3743 (in Chinese with English abstract)

[20] 王綸, 王星玉, 溫琪汾, 趙衛(wèi)紅, 劉金玉. 中國(guó)黍稷種質(zhì)資源抗黑穗病鑒定評(píng)價(jià). 植物遺傳資源學(xué)報(bào), 2008, 9: 497–501 Wang L, Wang X Y, Wen Q F, Zhao W H, Liu J Y. Identification and evaluation of resistance to dustbrand in Chinese proso millet germplasm resources., 2008, 9: 497–501 (in Chinese with English abstract)

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

[22] Rajput S G, Harveson T, Santra D K. Development and characterization of SSR markers in proso millet (L.) based on switch grass genomics., 2014, 5:175–186

[23] Rajput S G, Santra D K. Evaluation of genetic diversity of proso millet () germplasm available in the USA using SSR markers., 2016, doi: 10.2135/cropsci2015. 10.0644

[24] Hunt H V, Campana M G, Lawes M C, Park Y J, Bower M A, Howe C J, Jones M K. Genetic diversity and phylogeography of broomcorn millet (L.) across Eurasia., 2011, 20: 4756–4771

[25] 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). pii: E370. doi: 10.3390/ijms17030370

[26] 劉曉鑫, 謝傳曉, 趙琦, 路明, 曲延英, 張濱, 梁業(yè)紅, 張世煌. 基于SSR熒光標(biāo)記技術(shù)的玉米群體混合樣本基因頻率分析方法. 中國(guó)農(nóng)業(yè)科學(xué), 2008, 41: 3991–3998 Liu X X, Xie C X, Zhao Q, Lu M, Qu Y Y, Zhang B, Liang Y H, Zhang S H. Establishment of fluorescent SSR technique on detecting allelic frequency in maize (L.) populations with bulk sampling strategy., 2008, 41: 3991–3998 (in Chinese with English abstract)

[27] 郝晨陽(yáng), 王蘭芬, 賈繼增, 董玉琛, 張學(xué)勇. SSR熒光標(biāo)記和銀染技術(shù)的比較分析. 作物學(xué)報(bào), 2005, 31: 144–149 Hao C Y, Wang L F, Jia J Z, Dong Y C, Zhang X Y. Comparison of fluorescence and silver-staining detection systems of microsatellite markers., 2005, 31: 144–149 (in Chinese with English abstract)

[28] 李會(huì)勇, 王天宇, 黎裕, 石云素, 宋燕春, 陸平. TP-M13自動(dòng)熒光檢測(cè)法在高梁SSR基因型鑒定中的應(yīng)用. 植物遺傳資源學(xué)報(bào), 2005, 6: 68–70 Li H Y, Wang T Y, Li Y, Shi Y S, Song Y C, Lu P. Application of the TP-M13 automated fluorescent-lablled system of SSR genotyping in sorghum., 2005, 6: 68–70 (in Chinese with English abstract)

[29] 柴巖. 糜子. 北京: 中國(guó)農(nóng)業(yè)出版社, 1999. pp 165–168 Chai Y. Proso Millet (L.). Beijing: China Agriculture Press, 1999. pp 165–168 (in Chinese)

[30] 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: 357–364

[31] Edwards K, Johnstone C, Thompson C. A simple and rapid method for the preparation of plant genomic DNA for PCR analysis., 1991, 19: 1349

[32] 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: 267–273

[33] Liu K, Muse S V. PowerMarker: integrated analysis environment for genetic marker data., 2005, 21: 2128–2129

[34] Yeh F C, Boyle T J. Population genetic analysis of co-dominant and dominant markers and quantitative traits., 1997, 129: 157

[35] Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance and maximum parsimony methods., 2011, 28: 2731–2739

[36] Falush D, Stephens M, Pritchard J K. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies., 2003, 164: 1567–1587

[37] Peakall R, Smouse P E. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update., 2012, 8:2537–2539

[38] Evanno G, Regnaut S, Goudet J. Detecting the number of cluster of individuals using the software structure: a simulation study., 2005, 14: 2611–2620

[39] Li C, Dong Y, Liu M, Lu P, Li W, Wang Y, Cui X, Zhou H, Xu Y. Ancient DNA analysis of(broomcorn millet) from a Bronze Age cemetery in Xinjiang, China., 2016, 25: 469–477

[40] Karam D, Westra P, Nissen S J, Ward S M, Figueiredo J E F. Assessment of silver-stained AFLP markers for studying DNA polymorphism in proso millet (L.)., 2006, 29: 609–615

[41] Karam D, Westra P, Nissen S J, Ward S M, Figueiredo J E F. Genetic diversity among proso millet () biotypes assessed by AFLP technique., 2004, 22: 167–174

[42] M’Ribu H K, Hilu K W. Detection of interspecific and intraspecific variation inmillets through random amplified polymorphic DNA., 1994, 88: 412–416

[43] 王銀月, 劉敏軒, 陸平, 喬治軍, 楊天育, 李海, 崔喜艷. 構(gòu)建黍稷分子遺傳圖譜SSR引物的篩選. 作物雜志, 2014, (4): 32–38 Wang Y Y, Liu M X, Lu P, Qiao Z J, Yang T Y, Li H, Cui X Y. SSR primer selection to construct molecular genetic map for proso millet., 2014, (4): 32–38 (in Chinese with English abstract)

[44] Rajput S G, Santra D K, Schnable J. Mapping QTLs for morpho-agronomic traits in proso millet (L.)., 2016, 36: 37, DOI: 10.1007/s11032-016-0460-4

[45] 季煦, 王瑞云, 曳水瑛, 劉曉歡, 楊陽(yáng), 劉笑瑜, 連帥. 基于SSR的中國(guó)糜子遺傳多樣性分析. 山西農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2015, 35: 137–142Ji X, Wang R Y, Ye S Y, Liu X H, Yang Y, Liu X Y, Lian S. The analysis of genetic diversity in broomcorn millet (L.) based on SSR markers.(Nat Sci Edn), 2015, 35: 137–142 (in Chinese with English abstract)

[46] 連帥, 王瑞云, 馬躍敏, 劉笑瑜, 季煦. 不同生態(tài)區(qū)糜子種質(zhì)資源的遺傳多樣性分析. 山西農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2015, 35: 225–231 Lian S, Wang R Y, Ma Y M, Liu X Y, Ji X. Genetic diversity of broomcorn millet (L.) germplasms of different ecotype zone of China.(Nat Sci Edn), 2015, 35: 225–231 (in Chinese with English abstract)

[47] 王綸, 王星玉, 溫琪汾, 師穎. 山西是黍稷的起源和遺傳多樣性中心. 植物遺傳資源學(xué)報(bào), 2009, 10: 465–470 Wang L, Wang X Y, Wen Q F, Shi Y. Shanxi province is the center of the origin and genetic diversity of proso millet., 2009, 10: 465–470 (in Chinese with English abstract)

[48] Frankham R, Ballou J D, Briscoe D. Introduction to Conservation Genetics. Cambridge: Cambridge University Press, 2010. pp 309–336

(附圖1)

(附圖1)

附圖1 64號(hào)糜子樣品15個(gè)引物擴(kuò)增產(chǎn)物熒光標(biāo)記峰型圖

Supplementary fig. 1 SSR amplicon graph of common millet No. 64 sample by 15 primers

Analysis of Genetic Diversity in Common Millet () Using Fluorescent SSR in China

Wang Rui-Yun1,2,*,**, Ji Xu1,**, Lu Ping3, Liu Min-Xuan3, XU Yue4, Wang Lun2, Wang Hai-Gang2, and Qiao Zhi-Jun2,*

1College of Agriculture, Shanxi Agricultural University, Taigu 030801, China;2Institute 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, China;3Institute of Crop Science, Chinese Academy of Agricultural Sciences, Beijing 100081, China;4School of Life Sciences, Jilin University, Changchun 130012, China

Understanding of the genetic diversity of common millet germplasm resources is a basis of research on its origin and evolution, as well as the effective utilization of elite germplasm. In this study, we analyzed the genetic diversity of 132 common millet accessions collected from 11 provinces of China using 15 millet-specific fluorescent-labelled simple sequence repeat (SSR) markers. A total of 107 alleles were detected with 7 alleles per locus ranging from 2 to 14. The expected heterozygosity ranged from 0.0936 to 0.8676 with the mean value of 0.5298. The polymorphism information content (PIC) ranged from 0.0893 to 0.8538 with the mean value of 0.4864. The 132 accessions were clustered into four groups according to genetic distance between accessions, which showed geographic and eco-regional distributions. Group I is a collection of germplasm from spring-sowing area of Northeast China, group II is composed of partial accessions from the spring and summer-sowing areas of Loess Plateau, group III is composed of partial accessions from the spring-sowing area of Northern China, and the remaining accessions comprises Group IV with geographic overlaps to groups II and III. Model-based genetic structure analysis indicated that the genetic resources of Chinese common millet were from four gene pools in Northeast, Loess Plateau, Northern, and Northwest China. Obviously, this result was very similar to that of clustering analysis, showing the importance of geographical relationship. Rich genetic variations were also found in Chinese common millet resources, mainly from different accessions. The above results provide an overall evaluation of genetic diversity of common millet in China at a molecular level.

Common millet (L.); Fluorescent SSR; Genetic diversity; Cluster analysis; Genetic structure

10.3724/SP.J.1006.2017.00530

本研究由國(guó)家自然科學(xué)基金項(xiàng)目(31271791), 山西省回國(guó)留學(xué)人員科研資助項(xiàng)目(2016-066), 國(guó)家自然科學(xué)基金青年科學(xué)基金(31301386), 國(guó)家自然科學(xué)基金青年科學(xué)基金(31300279), 國(guó)家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)(CARS-07-13.5-A12), 山西省重點(diǎn)研發(fā)計(jì)劃(一般項(xiàng)目)(農(nóng)業(yè))項(xiàng)目(201603D221003-5)資助。

This study was supported by the National Natural Science Foundation of China (31271791), the Scientific Research Foundation for Returned Scholars in Shanxi province (2016-066), the National Natural Science Foundation of China for Young Scientists (31301386, 31300279), the China Agriculture Research System (CARS-07-13.5-A12), and General Project of the Shanxi Provincial Key Research and Development Program for Agriculture (201603D221003-5).

王瑞云, E-mail: wry925@126.com, Tel: 15234420135; 喬治軍, E-mail: nkypzs@126.com, Tel: 0351-7065530

**同等貢獻(xiàn)(Contributed equally to this work)

2016-08-12;

Accepted(接受日期): 2016-11-02;

Published online(網(wǎng)絡(luò)出版日期): 2016-11-15.

URL: http://www.cnki.net/kcms/detail/11.1809.S.20161115.1618.006.html