王海崗，溫琪汾，穆志新，喬治軍

山西谷子核心資源群體結(jié)構(gòu)及主要農(nóng)藝性狀關(guān)聯(lián)分析

王海崗，溫琪汾，穆志新，喬治軍

（山西省農(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）

【】分析山西谷子地方品種遺傳多樣性和群體遺傳結(jié)構(gòu)，篩選與谷子農(nóng)藝性狀相關(guān)聯(lián)的分子標(biāo)記，為谷子雜交組合親本選配及分子標(biāo)記輔助育種提供依據(jù)。利用96對(duì)SSR標(biāo)記對(duì)595份山西谷子核心資源進(jìn)行全基因組掃描，采用PowerMarker 3.25軟件分析群體遺傳多樣性，利用STRUCTURE 2.3.4軟件分析群體遺傳結(jié)構(gòu)，使用TASSEL 2.1軟件中GLM（general linear model，Q）和MLM（mixed linear model，Q+K）2種方法，進(jìn)行表型和標(biāo)記關(guān)聯(lián)分析。96對(duì)SSR引物共擴(kuò)增出828個(gè)等位變異，平均每對(duì)引物擴(kuò)增到8.6個(gè)，變化范圍為2—26；基因多樣性指數(shù)變化范圍為0.005—0.941，平均為0.610；多態(tài)信息量變化范圍為0.005—0.938，平均為0.577；各位點(diǎn)雜合度變化范圍為0—0.050，平均位點(diǎn)雜合度僅為0.016。群體結(jié)構(gòu)分析將595份核心資源分為3個(gè)亞群。4 560個(gè)SSR位點(diǎn)成對(duì)組合中，共線性組合和非共線性組合之間都存在一定的連鎖不平衡。統(tǒng)計(jì)概率（<0.01）支持的LD成對(duì)位點(diǎn)1 955個(gè)，占全部位點(diǎn)組合的42.9%，平均值為0.23。通過(guò)GLM方法共檢測(cè)到12個(gè)極顯著性位點(diǎn)（＜0.01），表型變異解釋率為2.34%—13.94%，平均為6.33%，貢獻(xiàn)率較高的等位變異位點(diǎn)是CAAS2050（2=13.94%）和B153（2=11.36%）；通過(guò)MLM方法共檢測(cè)到9個(gè)極顯著性位點(diǎn)（＜0.01），表型變異解釋率為2.80%—9.22%，平均為5.16%，貢獻(xiàn)率較高的等位變異位點(diǎn)是P89（2=9.22%）和P3*（2=8.28%）；2種方法共同檢測(cè)到的極顯著性位點(diǎn)有7個(gè)。利用SSR標(biāo)記分析了595份山西谷子核心資源的遺傳多樣性和群體遺傳結(jié)構(gòu)。2種關(guān)聯(lián)分析模型中，GLM方法關(guān)聯(lián)到12個(gè)標(biāo)記與節(jié)數(shù)、株高、頸長(zhǎng)、莖粗、穗長(zhǎng)、穗粗、碼數(shù)、碼粒數(shù)、蛋白質(zhì)含量9個(gè)性狀相關(guān)；MLM方法關(guān)聯(lián)到9個(gè)標(biāo)記與節(jié)數(shù)、頸長(zhǎng)、葉寬、莖粗、穗粗、碼數(shù)、碼粒數(shù)、千粒重8個(gè)性狀相關(guān)。

谷子；山西；地方品種；SSR；遺傳多樣性；連鎖不平衡；關(guān)聯(lián)分析

0 引言

【研究意義】谷子（（L.）Beauv.）是起源于中國(guó)黃河流域的糧飼兼用作物，又稱粟，去殼后為小米，在中國(guó)的栽培歷史可追溯到11 500年前[1]，是世界范圍內(nèi)最古老的禾谷類作物之一。山西谷子種質(zhì)資源豐富，占全國(guó)谷子種質(zhì)資源的1/5[2]，山西傳統(tǒng)的名優(yōu)谷子品種較多，如東方亮、沁州黃、晉谷21號(hào)、晉谷10號(hào)，這些品種蘊(yùn)藏著豐富的優(yōu)異基因資源。優(yōu)異基因發(fā)掘是作物種質(zhì)資源研究的重要部分，對(duì)作物分子育種具有重要的實(shí)踐意義。關(guān)聯(lián)分析是基因發(fā)掘的有效途徑之一，以連鎖不平衡為基礎(chǔ)，通過(guò)對(duì)自然群體表型數(shù)據(jù)和基因型數(shù)據(jù)的統(tǒng)計(jì)檢測(cè)，尋找性狀可遺傳變異和標(biāo)記多態(tài)性之間的關(guān)聯(lián)[3-5]。【前人研究進(jìn)展】目前，關(guān)聯(lián)分析已在水稻庫(kù)源相關(guān)性狀[6]、紋枯病抗性[7]等；小麥株高[8]、旗葉葉綠素[9]、穗發(fā)芽抗性[10]等性狀；大豆農(nóng)藝性狀[9]、加工和品質(zhì)性狀[12-14]、幼苗期耐淹性[15]及棉花農(nóng)藝和纖維品質(zhì)性狀[16]、葉綠素含量[17]、適宜機(jī)采相關(guān)性狀[18]等已進(jìn)行了大量研究。且在大麥[19-21]、谷子[22-25]、小豆[28]、豇豆[29]等小宗作物也進(jìn)行了一定研究。山西谷子年播種面積在22.46×104hm2左右，占全國(guó)的1/4[30]，但是作為山西傳統(tǒng)雜糧作物的谷子分子輔助育種相對(duì)緩慢，充分挖掘山西本土谷子資源中的優(yōu)良等位變異，利用分子標(biāo)記輔助選擇，創(chuàng)制谷子新種質(zhì)是谷子資源研究急需開展的工作。【本研究切入點(diǎn)】本研究中所用材料主要來(lái)自20世紀(jì)50—80年代征集的山西谷子地方品種[31-32]，這些材料未經(jīng)過(guò)育種的選擇壓力，多數(shù)資源是通過(guò)自然選擇或突變后篩選而保留下來(lái)的種質(zhì)，是環(huán)境和歷史的載體材料。傳統(tǒng)育種主要根據(jù)作物的表型和產(chǎn)量性狀進(jìn)行雜交親本組合配制，對(duì)遺傳組成的差異考慮較少，導(dǎo)致親本來(lái)源單一，育成品種間遺傳基礎(chǔ)狹窄，很難培育出具有突破性的優(yōu)良品種?！緮M解決的關(guān)鍵問(wèn)題】本研究對(duì)595份山西谷子地方品種初選核心資源進(jìn)行遺傳多樣性分析，明確山西谷子地方品種的遺傳結(jié)構(gòu)，尋找與農(nóng)藝性狀顯著相關(guān)聯(lián)的分子標(biāo)記，為谷子種質(zhì)創(chuàng)新中的優(yōu)異等位變異發(fā)掘、親本組配和標(biāo)記輔助選擇等提供參考。

1 材料與方法

1.1 試驗(yàn)材料及表型鑒定

所用595份山西谷子地方品種來(lái)源于5 627份山西谷子資源構(gòu)建的初選核心資源[30]及近年新收集的谷子資源。2017—2018年在山西省農(nóng)業(yè)科學(xué)院東陽(yáng)試驗(yàn)基地種植，對(duì)抽穗期、節(jié)數(shù)、株高、頸長(zhǎng)、葉長(zhǎng)、葉寬、穗長(zhǎng)、穗粗、莖粗、單穗重、穗粒重、碼數(shù)、碼粒數(shù)、千粒重、蛋白質(zhì)含量等15個(gè)性狀進(jìn)行鑒定。

1.2 基因組DNA提取

谷子抽穗后取倒二葉，液氮速凍，-80℃保存?zhèn)溆谩２捎弥参锘蚪MDNA提取試劑盒（DP305-03）提取DNA，利用全功能酶標(biāo)儀（美國(guó)伯騰SynergyH1）測(cè)定DNA濃度，-20℃保存?zhèn)溆谩?/p>

1.3 SSR分析

選取分布于9條染色體上的96對(duì)引物進(jìn)行多態(tài)性掃描（表1）[22-27]。SSR反應(yīng)體系包括模板DNA 1.0 μL、2×Tap PCR Mastermix 5.0 μL、10 mmol·L-1Forward primer 0.4 μL、10 mmol·L-1Reverse primer 0.4 μL和ddH2O 3.2 μL。PCR反應(yīng)程序?yàn)?5℃ 5 min；95℃ 30 s，58℃（不同引物退火溫度根據(jù)溫度梯度試驗(yàn)確定）30 s，72℃ 60 s，32個(gè)循環(huán)；72℃ 10 min；4℃保存。擴(kuò)增產(chǎn)物用8%聚丙烯酰胺凝膠電泳分離，銀染顯色。

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

根據(jù)SSR擴(kuò)增產(chǎn)物的電泳結(jié)果，參照DNA Ladder記錄擴(kuò)增產(chǎn)物條帶的分子量大小。使用Powermarker 3.25軟件計(jì)算每個(gè)位點(diǎn)等位基因數(shù)、基因型數(shù)、基因多樣性指數(shù)、多態(tài)性信息含量指數(shù)和雜合度，并計(jì)算Nei遺傳距離采用MEGA 5.0軟件構(gòu)建NJ聚類圖。

采用Structure2.3.4軟件分析群體遺傳結(jié)構(gòu)，K取值為1—10，3次重復(fù)；將MCMC（markov chain monto carlo）開始時(shí)的不作數(shù)迭代設(shè)為100 000次，再將不作數(shù)迭代后的MCMC設(shè)為100 000次，其余參數(shù)采用軟件默認(rèn)的設(shè)置。根據(jù)lnP(D)計(jì)算ΔK，并依據(jù)ΔK值選擇一個(gè)合適的K值。使用Tassel 2.1軟件計(jì)算SSR位點(diǎn)組合間的連鎖不平衡（linkage disequilibrium，LD）水平及支持概率，繪制LD配對(duì)檢測(cè)矩陣圖，使用一般線性模型（general linear model，GLM）和混合線性模型（mixed liner model，MLM）結(jié)合基因型數(shù)據(jù)、表型值和Q值，對(duì)SSR標(biāo)記與表型性狀進(jìn)行關(guān)聯(lián)分析，計(jì)算<0.01時(shí)標(biāo)記位點(diǎn)對(duì)表型變異的貢獻(xiàn)率（2）。采用STRUCTURE軟件計(jì)算Q值，用軟件SPAGeDi1.5軟件計(jì)算Kinship值。

2 結(jié)果

2.1 SSR多樣性

96對(duì)SSR引物在595份山西谷子地方品種核心資源中共檢測(cè)到828個(gè)等位變異，變化范圍2—26個(gè)，平均為8.6個(gè)，其中，標(biāo)記B109檢測(cè)到的等位變異最多為26個(gè)（表1）?；蛐蛿?shù)變化范圍在2—47，平均每個(gè)位點(diǎn)為13?；蚨鄻有灾笖?shù)變化范圍為0.0050—0.9411，平均為0.6097；多態(tài)信息量變化范圍為0.0050—0.9379，平均為0.5773；各位點(diǎn)雜合度變化范圍為0—0.050，平均位點(diǎn)雜合度僅為0.0159，說(shuō)明谷子異交率很小接近于純系。每條染色體上標(biāo)記平均數(shù)有10對(duì)，平均等位變異為92個(gè)；第3和第8染色體檢測(cè)到等位變異最多，均為107個(gè)，第7染色體最少為79個(gè)。第3染色體基因多樣性指數(shù)和多態(tài)性信息含量最高，分別為0.720和0.6912，第7染色體最低，分為0.5322和0.4940。第7染色體雜合度最小，為0.0114，第1染色體雜合度最大，為0.0212。

表1 96對(duì)SSR引物遺傳多樣性

續(xù)表1 Continued table 1

續(xù)表1 Continued table 1

2.2 群體結(jié)構(gòu)分析

采用Structure 2.3.4軟件對(duì)595份山西谷子核心資源進(jìn)行群體遺傳結(jié)構(gòu)分析，計(jì)算每份材料的Q值（每份種質(zhì)被分配到對(duì)應(yīng)亞群中的概率）。隨著K值的增大，LnP(D)遞增，未出現(xiàn)拐點(diǎn)，因此無(wú)法準(zhǔn)確判斷K值（圖1）。參照EVANNO等[33]方法，當(dāng)K=3時(shí)，ΔK出現(xiàn)最大值，因此該群體材料被劃分為3個(gè)亞群（圖2）。按照亞群中個(gè)體的概率≥0.6時(shí)，個(gè)體被分配到相應(yīng)聚類的亞群；亞群中個(gè)體的概率≤0.6時(shí)，個(gè)體被分配到一個(gè)混合類群的分配原則，分析每份核心資源對(duì)應(yīng)的Q值。結(jié)果表明，595份材料歸屬為3個(gè)亞群，分別包含77、222和187份材料，3個(gè)亞群各占比例為12.9%、37.3%和31.4%；另有109份材料在任意一個(gè)亞群內(nèi)Q值均＜0.6，單獨(dú)劃為一個(gè)混合群，占比例為18.3%。

基于SSR遺傳距離的11個(gè)地市聚類分析可以看出（圖3），運(yùn)城、晉城、臨汾、長(zhǎng)治、陽(yáng)泉聚為第一類；大同、朔州、忻州、太原聚為第二類；晉中、呂梁聚為第3類。劃分的3個(gè)類群按照緯度高低自北向南依次分開，這與山西復(fù)雜多樣的地理生態(tài)類型有一定關(guān)系，說(shuō)明山西谷子種質(zhì)資源存在地理生態(tài)型的分化。

A：ln P(D)值與K值的變化圖；B：ΔK值隨K值的變化圖 A: Line chart of K and ln P(D); B: magnitude of ΔK as a function of K

圖2 595份山西谷子地方品種的群體遺傳結(jié)構(gòu)

2.3 連鎖不平衡分析

關(guān)聯(lián)分析以不同位點(diǎn)等位基因間的連鎖不平衡為基礎(chǔ)，描述群體的重組史。對(duì)分布于谷子全基因組96對(duì)SSR標(biāo)記進(jìn)行連鎖不平衡分析，4 560個(gè)SSR位點(diǎn)成對(duì)組合中，共線性組合（同一染色體）和非共線性組合（不同染色體）之間都存在一定的連鎖不平衡（圖4）。統(tǒng)計(jì)概率（<0.01）支持的LD成對(duì)位點(diǎn)1 955個(gè)，占全部位點(diǎn)組合的42.9%，平均值為0.23。LD成對(duì)位點(diǎn)的值范圍主要集中在0—0.2和0.2—0.4（表2），在第2和第4染色體上未發(fā)現(xiàn)在0.5以上的LD位點(diǎn)，說(shuō)明這兩條染色體較其他7條染色體受到的選擇壓力??；分析共線性位點(diǎn)組合的LD，發(fā)現(xiàn)第4、第6和第9染色體上的LD成對(duì)位點(diǎn)分別有31、35和36個(gè)，高于其他染色體，且第6染色體上平均值最高為0.28，說(shuō)明群體材料中，該染色體上發(fā)生重組選擇要大于其他染色體。

2.4 SSR標(biāo)記與性狀關(guān)聯(lián)分析

采用Tassel 2.1軟件分析了595份山西谷子核心資源在一點(diǎn)2年2種模型的15個(gè)性狀表型和基因型，結(jié)果顯示，在<0.01下，共檢測(cè)到14個(gè)SSR標(biāo)記與11個(gè)性狀相關(guān)聯(lián)（表3）。GLM方法檢測(cè)到12個(gè)標(biāo)記與節(jié)數(shù)、株高、頸長(zhǎng)、莖粗、穗長(zhǎng)、穗粗、碼數(shù)、碼粒數(shù)和蛋白質(zhì)9個(gè)性狀相關(guān)聯(lián)。表型變異解釋率為2.34%—13.94%，平均值為6.33%。表型變異解釋率較大的標(biāo)記是CAAS2050和B153，解釋率分別為13.94%和11.36%；MLM方法檢測(cè)到9個(gè)標(biāo)記與節(jié)數(shù)、頸長(zhǎng)、葉寬、莖粗、穗粗、碼數(shù)、碼粒數(shù)和千粒重8個(gè)性狀相關(guān)聯(lián)。表型變異解釋率為2.80%—9.22%，平均為5.16%。貢獻(xiàn)率較高的位點(diǎn)是P89和P3*，解釋率分別為9.22%和8.28%。比較GLM和MLM結(jié)果，CAAS6023、P3*、In6-2、CAAS2050和B117 5個(gè)標(biāo)記在2種模型條件下同時(shí)被檢測(cè)到。B153、P58、B117、In6-2、CAAS6023和B142標(biāo)記同時(shí)與2個(gè)或多個(gè)性狀相關(guān)聯(lián)，這可能與性狀相關(guān)有一定關(guān)系，說(shuō)明控制該性狀的等位基因可能連鎖或一因多效。

與莖粗相關(guān)的標(biāo)記有5個(gè)，表型變異解釋率為4.45%—9.77%，平均值為6.72%，貢獻(xiàn)率最高的位點(diǎn)為P89；與碼數(shù)相關(guān)的標(biāo)記有4個(gè)，表型變異解釋率為2.58%—13.94%，平均值為5.76%，貢獻(xiàn)率最高的位點(diǎn)為CAAS2050；與蛋白質(zhì)相關(guān)的位點(diǎn)有3個(gè)，表型變異解釋率為4.25%—11.36%，平均值為8.61%，貢獻(xiàn)率最高的位點(diǎn)為B153；與節(jié)數(shù)相關(guān)的位點(diǎn)有2個(gè)，表型變異解釋率為5.20%—5.97%，平均值為5.63%，貢獻(xiàn)率最高的位點(diǎn)為P58；與穗長(zhǎng)相關(guān)的位點(diǎn)有2個(gè)，表型變異解釋率為2.34%—6.05%，平均值為4.16%，貢獻(xiàn)率最高的位點(diǎn)為B153；與穗粗相關(guān)的位點(diǎn)有2個(gè)，表型變異解釋率為4.57%—6.46%，平均值為5.94%，貢獻(xiàn)率最高的位點(diǎn)為B142。

表2 不同D'值范圍LD成對(duì)位點(diǎn)數(shù)

圖4 9條染色體SSR位點(diǎn)的連鎖不平衡分布

3 討論

3.1 山西谷子資源的遺傳多樣性

優(yōu)異基因發(fā)掘是種質(zhì)資源高效利用和作物遺傳改良的基礎(chǔ)，種質(zhì)資源的遺傳多樣性評(píng)價(jià)和群體遺傳結(jié)構(gòu)的分析是優(yōu)異基因發(fā)掘的前提條件。關(guān)聯(lián)分析將蘊(yùn)藏在不同種質(zhì)中的等位變異與表型性狀建立了聯(lián)系，通過(guò)尋找顯著關(guān)聯(lián)位點(diǎn)，可以發(fā)掘種質(zhì)資源中優(yōu)異等位基因。山西谷子資源豐富，栽培歷史悠久，名優(yōu)品種眾多，但是傳統(tǒng)谷子育種只是簡(jiǎn)單的品種雜交，通過(guò)應(yīng)用分子標(biāo)記輔助選擇聚合有利基因，改良谷子地方品種，是山西谷子育種和種質(zhì)創(chuàng)制的新途徑。本研究采用SSR標(biāo)記對(duì)山西谷子地方品種核心資源進(jìn)行遺傳多樣性評(píng)價(jià)，分析群體遺傳結(jié)構(gòu)，挖掘與農(nóng)藝性狀關(guān)聯(lián)的分子標(biāo)記，對(duì)指導(dǎo)山西谷子資源的利用和育種親本選擇具有指導(dǎo)意義。

通過(guò)對(duì)595份山西谷子核心資源進(jìn)行分子標(biāo)記掃描，發(fā)現(xiàn)平均等位基因8.6個(gè)，平均基因多樣性指數(shù)0.610，平均PIC為0.577，低于Wang等[22]和Jia等[23]對(duì)全國(guó)谷子地方品種和育成品種研究結(jié)果，同時(shí)也低于Vetriventhan等[24]的研究，說(shuō)明山西谷子資源的遺傳多樣性要低于全國(guó)和世界的谷子資源。但與GUPTA等[25]和Jia等[34]研究結(jié)果相近，這與等位變異統(tǒng)計(jì)方法和試驗(yàn)選擇引物有一定關(guān)系。山西行政區(qū)域從北緯34°34′到40°44′橫跨6個(gè)緯度，海拔從180 m到3 058 m，地貌類型復(fù)雜多樣，全省各縣均有谷子資源分布，形成了多類型谷子種植區(qū)，按照緯度和海拔分為春播早熟區(qū)、春播中熟區(qū)、春播晚熟區(qū)和夏播區(qū)[35]。山西谷子核心資源按照11個(gè)地市來(lái)源劃分成3個(gè)類群，基本按照緯度高低由北到南劃分開來(lái)[36]，說(shuō)明山西谷子種質(zhì)資源存在地理生態(tài)型分化，也從分子水平證明了山西省是春、夏谷2種類型同時(shí)存在的省份。同時(shí)，明確了各地市谷子資源的遺傳距離，為山西谷子生態(tài)區(qū)域劃分和育種親本選配提供了分子證據(jù)。

表3 GLM和MLM關(guān)聯(lián)分析

SN：節(jié)數(shù)；PH：株高；PL：頸長(zhǎng)；LW：葉寬；DMS：莖粗；PL：穗長(zhǎng)；PD：穗粗；PBNP：碼數(shù)；SNPB：碼粒數(shù)；GW：千粒重；PC：蛋白含量；GLM：一般線性模型；MLM：混合線性模型

SN: stem node number; PH: plant height; PL: peduncle length; LW: Leaf width; DMS: diameter of main stem; PL: panicle length; PD: panicle diameter; PBNP: primary branch number per panicle; SNPB: spikelet number per primary branch; GW: 1000-grain weight; PC: protein content; GLM: general linear model; MLM: mixed linear model

3.2 SSR位點(diǎn)連鎖不平衡及關(guān)聯(lián)分析

已有研究表明[22,24]，谷子LD衰減距離在15—40 cM，說(shuō)明谷子的LD水平較高，有利于進(jìn)行全基因組關(guān)聯(lián)分析，進(jìn)而發(fā)掘種質(zhì)資源中的有利基因。本研究利用覆蓋全基因組的96對(duì)SSR引物進(jìn)行標(biāo)記掃描，4 560個(gè)位點(diǎn)成對(duì)組合中，不論是共線性組合還是非共線性組合，都有LD存在。統(tǒng)計(jì)概率（<0.01）支持的LD成對(duì)位點(diǎn)占42.9%，較高水平的LD成對(duì)位點(diǎn)（>0.5）主要分布在第1、第2、第6和第9染色體上及與其組合的位點(diǎn)，說(shuō)明這4條染色體面臨的選擇壓力大于其他染色體，這些材料多經(jīng)歷了較多農(nóng)藝性狀的選擇。本研究發(fā)現(xiàn)核心資源的整體連鎖不平衡水平較低（平均值為0.23），說(shuō)明山西谷子資源在經(jīng)歷人為重組和自然突變的歷程較少，且谷子是自花授粉作物，自然異交率低，地方品種中攜帶的大量?jī)?yōu)異等位變異仍單獨(dú)存在，所以山西谷子育種仍有進(jìn)一步提升空間。關(guān)聯(lián)分析在極顯著（<0.01）條件下共檢測(cè)到14個(gè)與農(nóng)藝性狀和蛋白質(zhì)含量相關(guān)的位點(diǎn)，對(duì)表型變異解釋率平均值為5.98%，在谷子育種中可利用這些標(biāo)記進(jìn)行材料的等位變異標(biāo)記輔助選擇研究；另外，有6個(gè)標(biāo)記位點(diǎn)同時(shí)與2個(gè)以上性狀相關(guān)聯(lián)，其中B142位點(diǎn)與4個(gè)性狀相關(guān)聯(lián)，這可能是性狀相關(guān)和一因多效的遺傳基礎(chǔ)，說(shuō)明這些標(biāo)記可以同時(shí)進(jìn)行改良多個(gè)目標(biāo)性狀，可用于多個(gè)性狀的聚合育種研究。

4 結(jié)論

595份谷子核心資源遺傳結(jié)構(gòu)存在3個(gè)類群，11個(gè)地市按照緯度由高到低劃分成3類，山西谷子資源存在地理生態(tài)型分化。共線性或非共線性位點(diǎn)組合都有LD存在，核心資源群體有利于關(guān)聯(lián)分析和優(yōu)異等位變異挖掘。2種關(guān)聯(lián)模型下分別找到12個(gè)與節(jié)數(shù)、株高、頸長(zhǎng)、莖粗、穗長(zhǎng)、穗粗、碼數(shù)、碼粒數(shù)、蛋白質(zhì)含量相關(guān)聯(lián)的標(biāo)記和9個(gè)與節(jié)數(shù)、頸長(zhǎng)、葉寬、莖粗、穗粗、碼數(shù)、碼粒數(shù)、千粒重相關(guān)聯(lián)的標(biāo)記。

[1] 賈冠清, 刁現(xiàn)民. 谷子((L.) P. Beauv.)作為功能基因組研究模式植物的發(fā)展現(xiàn)狀及趨勢(shì).生命科學(xué), 2017, 29(3): 292-301.

JIA G Q, DIAO X M. Current status and perspectives of researches on foxtail millet ((L.) P. Beauv.): A potential model of plant functional genomics studies., 2017, 29(3): 292-301. (in Chinese)

[2] 王海崗, 秦慧彬, 田翔, 呂建珍, 陳凌, 王君杰, 曹曉寧, 劉思辰, 王綸, 溫淇汾, 穆志新, 喬治軍. 山西谷子地方品種表型多樣性分析. 中國(guó)農(nóng)學(xué)通報(bào), 2018, 34(32): 19-24.

WANG H G, QIN H B, TIAN X, Lü J Z, CHEN L, WANG J J, CAO X N, LIU S C, WANG L, WEN Q F, MU Z X, QIAO Z J. Phenotypic diversity of local foxtail millet varieties in Shanxi., 2018, 34(32): 19-24. (in Chinese)

[3] 王榮煥, 王天宇, 黎裕. 關(guān)聯(lián)分析在作物種質(zhì)資源分子評(píng)價(jià)中的應(yīng)用. 植物遺傳資源學(xué)報(bào), 2007, 8(3): 366-372.

WANG R H, WANG T Y, LI Y. Application of association analysis in molecular evaluation of crop germplasm resources., 2007, 8(3): 366-372. (in Chinese)

[4] 楊小紅, 嚴(yán)建兵, 鄭艷萍, 余建明, 李建生. 植物數(shù)量性狀關(guān)聯(lián)分析研究進(jìn)展. 作物學(xué)報(bào), 2007, 33(4): 523-530.

YANG X H, YAN J B, ZHENG Y P, YU J M, LI J S. Reviews of association analysis for quantitative traits in plants., 2007, 33(4): 523-530. (in Chinese)

[5] 馮建英, 溫陽(yáng)俊, 張瑾, 章元明. 植物關(guān)聯(lián)分析方法的研究進(jìn)展. 作物學(xué)報(bào), 2016, 42(7): 945-956.

FENG J Y, WEN Y J, ZHANG J, ZHANG Y M. Advances on methodologies for genome-wide association studies in plants., 2016, 42(7): 945-956. (in Chinese)

[6] 趙宏亮, 陳凱, 張強(qiáng), 徐正進(jìn), 徐建龍, 黎志康. 基于連鎖不平衡水稻源庫(kù)相關(guān)性狀的關(guān)聯(lián)分析. 核農(nóng)學(xué)報(bào), 2015, 29(4) : 674-684.

ZHAO H L, CHEN K, ZHANG Q, XU Z J, XU J L, LI Z K. Association analysis on source sink-related traits in rice based on linkage disequilibrium., 2015, 29(4): 674-684. (in Chinese)

[7] 孫曉棠, 盧冬冬, 歐陽(yáng)林娟, 胡麗芳, 邊建民, 彭小松, 陳小榮, 傅軍如, 賀曉鵬, 賀浩華, 朱昌蘭. 水稻紋枯病抗性關(guān)聯(lián)分析及抗性等位變異發(fā)掘. 作物學(xué)報(bào), 2014, 40(5): 779-787.

SUN X T, LU D D, OUYANG L J, HU L F, BIAN J M, PENG X S, CHEN X R, FU J R, HE X P, HE H H, ZHU C L. Association mapping and resistant alleles analysis for sheath blight resistance in rice., 2014, 40(5): 779-787. (in Chinese)

[8] 魏添梅, 昌小平, 閔東紅, 景蕊蓮. 小麥抗旱品種的遺傳多樣性分析及株高優(yōu)異等位變異挖掘. 作物學(xué)報(bào), 2010, 36(6): 895-904.

WEI T M, CHANG X P, MIN D H, JING R L. Analysis of genetic diversity and tapping elite alleles for plant height in drought-tolerant wheat varieties., 2010, 36(6): 895-904. (in Chinese).

[9] 張冬玲, 張洪娜, 郝晨陽(yáng), 王蘭芬, 李甜, 張學(xué)勇. 花后冠層溫度對(duì)小麥產(chǎn)量的影響及幾個(gè)關(guān)聯(lián)SSR位點(diǎn)的效應(yīng)分析. 作物學(xué)報(bào), 2015, 41(4): 548-556.

ZHANG D L, ZHANG H N, HAO C Y, WANG L F, LI T, ZHANG X Y. Influence of canopy temperature (ct) during grain-filling period on yield and effects of several ct-associated ssr loci., 2015, 41(4): 548-556. (in Chinese)

[10] 朱玉磊, 王升星, 趙良俠, 張德新, 胡建幫, 曹雪連, 楊亞杰, 常成, 馬傳喜, 張海萍. 以關(guān)聯(lián)分析發(fā)掘小麥整穗發(fā)芽抗性基因分子標(biāo)記. 作物學(xué)報(bào), 2014, 40(10): 1725-1732.

ZHU Y L, WANG S X, ZHAO L X, ZHANG D X, HU J B, CAO X L, YANG Y J, CHANG C, MA C X, ZHANG H P. Exploring molecular markers of preharvest sprouting resistance gene using wheat intact spikes by association analysis., 2014, 40(10): 1725-1732. (in Chinese)

[11] 楊勝先, 牛遠(yuǎn), 李夢(mèng), 魏世平, 劉曉芬, 呂海燕, 章元明. 栽培大豆農(nóng)藝性狀的關(guān)聯(lián)分析及優(yōu)異等位變異挖掘. 中國(guó)農(nóng)業(yè)科學(xué), 2014, 47(20): 3941-3952.

YANG S X, NIU Y, LI M, WEI S P, LIU X F, Lü H Y, ZHANG Y M. Association mapping of agronomic traits in soybean (L. Merr.) and mining of novel alleles., 2014,47(20): 3941-3952. (in Chinese)

[12] 文自翔, 趙團(tuán)結(jié), 鄭永戰(zhàn), 劉順湖, 王春娥, 王芳, 蓋鈞鎰. 中國(guó)栽培和野生大豆農(nóng)藝品質(zhì)性狀與SSR標(biāo)記的關(guān)聯(lián)分析: I.群體結(jié)構(gòu)及關(guān)聯(lián)標(biāo)記. 作物學(xué)報(bào), 2008, 34(7): 1169-1178.

WEN Z X, ZHAO T J, ZHENG Y Z, LIU S H, WANG C E, WANG F, GAI J Y. Association analysis of agronomic and quality traits with ssr markers inandin china: i. population structure and associated markers., 2008, 34(7): 1169-1178. (in Chinese)

[13] 文自翔, 趙團(tuán)結(jié), 鄭永戰(zhàn), 劉順湖, 王春娥, 王芳, 蓋鈞鎰. 中國(guó)栽培和野生大豆農(nóng)藝及品質(zhì)性狀與SSR標(biāo)記的關(guān)聯(lián)分析: Ⅱ.優(yōu)異等位變異的發(fā)掘. 作物學(xué)報(bào), 2008, 34(8): 1339-1349.

WEN Z X, ZHAO T J, ZHENG Y Z, LIU S H, WANG C E, WANG F, GAI J Y. Association analysis of agronomic and quality traits with ssr markers inandin china: ii. exploration of elite alleles., 2008, 34(8): 1339-1349. (in Chinese)

[14] 范虎, 文自翔, 王春娥, 王芳, 邢光南, 趙團(tuán)結(jié), 蓋鈞鎰. 中國(guó)野生大豆群體農(nóng)藝加工性狀與SSR關(guān)聯(lián)分析和特異材料的遺傳構(gòu)成. 作物學(xué)報(bào), 2013, 39(5): 775-788.

FAN H, WEN Z X, WANG C E, WANG F, XING G N, ZHAO T J, GAI J Y. Association analysis between agronomic-processing traits and SSR markers and genetic dissection of specific accessions in Chinese wild soybean population., 2013, 39(5): 775-788. (in Chinese)

[15] 孫慧敏, 張軍, 趙團(tuán)結(jié), 蓋鈞鎰. 亞洲地區(qū)中、外大豆品種幼苗期耐淹性與SSR標(biāo)記的關(guān)聯(lián)分析. 作物學(xué)報(bào), 2010, 36(10): 1615-1623.

SUN H M, Zhang J, ZHAO T J, GAI J Y. Association analysis between submergence tolerance and SSR markers in domestic and foreign soybean cultivars in asia., 2010, 36(10): 1615-1623. (in Chinese)

[16] 聶新輝, 尤春源, 鮑健, 李曉方, 惠慧, 劉洪亮, 秦江鴻, 林忠旭. 基于關(guān)聯(lián)分析的新陸早棉花品種農(nóng)藝和纖維品質(zhì)性狀優(yōu)異等位基因挖掘. 中國(guó)農(nóng)業(yè)科學(xué), 2015, 48(15): 2891-2910.

NIE X H, YOU C Y, BAO J, LI X F, HUI H, LIU H L, QIN J H, LIN Z X. Exploration of elite alleles of agronomic and fiber quality traits in xinluzao cotton varieties by association analysis., 2015, 48(15): 2891-2910. (in Chinese)

[17] 劉其寶, 李黎貝, 張馳, 宿俊吉, 魏恒玲, 王寒濤, 喻樹迅. 陸地棉葉片葉綠素含量與SSR標(biāo)記的關(guān)聯(lián)分析及優(yōu)異等位變異的挖掘. 中國(guó)農(nóng)業(yè)科學(xué), 2017, 50(18): 3439-3449.

LIU Q B, LI L B, ZHANG C, SU J J, WEI H L, WANG H T, YU S X. Association analysis of leaf chlorophyll content with SSR markers and exploration of superior alleles in upland cotton., 2017, 50(18): 3439-3449. (in Chinese)

[18] 王娟, 董承光, 劉麗, 孔憲輝, 王旭文, 余渝. 棉花適宜機(jī)采相關(guān)性狀的SSR標(biāo)記關(guān)聯(lián)分析及優(yōu)異等位基因挖掘. 作物學(xué)報(bào), 2010, 36(10): 1615-1623.

WANG J, DONG C G, LIU L, KONG X H, WANG X W, YU Y. Association analysis and exploration of elite alleles of mechanical harvest-related traits with SSR markers in upland cotton cultivars (l)., 2010, 36(10): 1615-1623. (in Chinese)

[19] 賴勇, 王鵬喜, 范貴強(qiáng), 司二靜, 王晉, 楊軻, 孟亞雄, 李葆春, 馬小樂, 尚勛武, 王化俊. 大麥SSR標(biāo)記遺傳多樣性及其與農(nóng)藝性狀關(guān)聯(lián)分析. 中國(guó)農(nóng)業(yè)科學(xué), 2013, 46(2): 233-242.

LAI Y, WANG P X, FAN G Q, SI E J, WANG J, YANG K, MENG Y X, LI B C, MA X L, SHANG X W, WANG H J. Genetic diversity and association analysis using SSR markers in barley., 2013, 46(2): 233-242. (in Chinese)

[20] 姜曉東, 郭剛剛, 張京. 中國(guó)大麥地方品種的遺傳多樣性及α-淀粉酶活性的全基因組關(guān)聯(lián)分析. 中國(guó)農(nóng)業(yè)科學(xué), 2013, 46(4): 668-677.

JIANG X D, GUO G G, ZHANG J. Analysis on the genetic diversity in Chinese barley landrace and genomic wide association of α-amylase activity genes., 2013, 46(4): 668-677. (in Chinese)

[21] 賴勇, 王鵬喜, 范貴強(qiáng), 司二靜, 王晉, 楊軻, 孟亞雄, 李葆春, 馬小樂, 尚勛武, 王化俊. 大麥遺傳多樣性及連鎖不平衡分析. 作物學(xué)報(bào), 2013, 39(12): 2154-2161.

LAI Y, WANG P X, FAN G Q, SI E J, WANG J, YANG K, MENG Y X, LI B C, MA X L, SHANG X W, WANG H J. Genetic diversity and linkage disequilibrium analysis in barley., 2013, 39(12): 2154-2161. (in Chinese)

[22] WANG C F, JIA G Q, ZHI H, NIU Z G, CHAI Y, LI W, WANG Y F, LI H Q, LU P, ZHAO B H, DIAO X M. Genetic diversity and population structure of chinese foxtail millet [(L.) Beauv.] landraces., 2012, 7(2): 769-777.

[23] JIA G Q, LIU X T, SCHNABLE J C, NIU Z G, WANG C F, LI Y H, WANG S J, WANG S Y, LIU J R, GUO E H, ZHI H, DIAO X M. Microsatellite variations of elite setaria varieties released during last six decades in China., 10(5): e0125688.

[24] Vetriventhan M, Upadhyaya Hari D, Anandakumar C R, Senthilvel S, Varshney R K, Parzies H K. Population structure and linkage disequilibrium of ICRISAT foxtail millet ((L.) P. Beauv.) core collection., 2014, 196: 423-435

[25] GUPTA S, KUMARI K, MUTHAMILARASAN M, PARIDA W K, PRASAD M. Population structure and association mapping of yield contributing agronomic traits in foxtail millet., 2014, 33: 881-893

[26] 王春芳. 利用微衛(wèi)星標(biāo)記分析中國(guó)谷子地方品種的群體結(jié)構(gòu)與遺傳多樣性[D]. 石家莊: 河北師范大學(xué), 2011.

WANG C F. Population structure and genetic diversity of foxtail millet landraces in China assessed by microsatellites[D]. Shijiazhuang: Hebei Normal University, 2011. (in Chinese)

[27] Zhang S, Tang C J, Zhao Q, Li J, Yang L F, Qie L F, Fan X K, Li L, Zhang N, Zhao M C, Liu X T, Chai Y, Zhang X, Wang H L, Li Y T, Li W, Zhi H, Jia G Q, Diao X M. Development of highly polymorphic simple sequence repeat markers using genome- wide microsatellite variant analysis in foxtail millet [(L.) P. Beauv.]., 2014, 15: 78.

[28] 白鵬, 程須珍, 王麗俠, 王素華, 陳紅霖. 小豆遺傳差異、群體結(jié)構(gòu)和連鎖不平衡水平的SSR分析. 作物學(xué)報(bào), 2014, 40(5): 788-797.

BAI P, CHENG X Z, WANG L X, WANG S H, CHEN H L. Genetic diversity, population structure and linkage disequilibrium in adzuki bean by using SSR markers., 2014, 40(5): 788-797. (in Chinese)

[29] 熊海錚. 豇豆遺傳多樣性及若干農(nóng)藝性狀關(guān)聯(lián)分析[D]. 杭州: 浙江大學(xué), 2016.

XIONG H Z. Genetic diversity and association analysis of agronomic traits in Cowpea[D]. Hangzhou: Zhejiang University, 2016. (in Chinese)

[30] 王海崗, 溫琪汾, 喬治軍, 穆志新. 山西谷子地方品種初選核心種質(zhì)構(gòu)建. 農(nóng)學(xué)學(xué)報(bào), 2019, 9(4): 26-31.

WANG H G, WEN Q F, QIAO Z J, MU Z X. Construction of candidate core collection local varieties of foxtail millet in Shanxi., 2019, 9(4): 26-31. (in Chinese)

[31] 溫琪汾, 王綸, 趙衛(wèi)紅, 暢建武. 山西省谷子種質(zhì)資源的繁種入庫(kù). 山西農(nóng)業(yè)科學(xué), 2002 , 30(4): 32-34.

WEN Q F, WANG L, ZHAO W H, CHANG J W. Reproduction and storage of foxtail millet germplasm in Shanxi., 2002, 30(4): 32-34. (in Chinese)

[32] 溫琪汾, 王綸, 王星玉. 山西省谷子種質(zhì)資源及抗旱種質(zhì)的篩選利用. 山西農(nóng)業(yè)科學(xué), 2005, 33(4): 32-33.

WEN Q F, WANG L, WANG X Y. The foxtail millet germplasm resources and screening and utilization of drought resistance germplasm in Shanxi., 2005, 33(4): 32-33. (in Chinese)

[33] EVANNO G, REGNAUTS S, GOUDET J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study., 2005, 14: 2611-2620.

[34] JIA X P, ZHANG Z B, LIU Y H, ZHANG C W, SHI Y S, SONG Y C, WANG T Y, LI Y. Development and genetic mapping of SSR markers in foxtail millet [(L.) P. Beauv.]., 2009, 118(4): 821-829.

[35] 古兆明, 古世祿. 山西谷子起源與發(fā)展研究. 北京: 中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社, 2007: 166-168.

GU Z M, GU S L.. Beijing: Agricultural Science and Technology Press, 2007: 166-168. (in Chinese)

[36] 任繼海, 溫琪汾, 王星玉, 王綸. 山西省谷子品種類型及其分布. 華北農(nóng)學(xué)報(bào), 1999, 14(3): 25-30.

REN J H, WEN Q F, WANG X Y, WANG L. Varieties of foxtail millet and their distribution in Shanxi province., 1999, 14(3): 25-30. (in Chinese)

Population Structure and Association Analysis of Main Agronomic Traits of Shanxi Core Collection in Foxtail Millet

WANG HaiGang, WEN QiFen, MU ZhiXin, QIAO ZhiJun

(Institute of Crop Germplasm Resources of 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)

The objective of this study is to detect the SSR markers associated with agronomic trait and analyze genetic diversity and genetic structure of foxtail millet landrace in Shanxi province. The results will be helpful for hybridization combination of parent materials and molecular marker assisted breeding.96 SSR markers on 9 chromosomes were genome-wide screened for polymorphism in core collection of 595 accessions. PowerMarker 3.25 software was used to estimate the polymorphism information of population. Population structure was analyzed using STRUCTURE 2.3.4 software. Then the data were associated with 96 SSR markers by GLM (general linear model, Q) and MLM (mixed linear model, Q+K).Totally 828 alleles were found with 96 SSR markers and 8.6 alleles were revealed with each marker in average ranged from 2-26. The gene diversity was from 0.005 to 0.941, averagely 0.610. The polymorphism information content (PIC) value ranged from 0.005 to 0.938 with the mean of 0.577. Heterozygosity per locus on average was 0.016, ranging from 0 to 0.050.All the 595 accessions were divided into three subgroups by analysis of population genetic structure. There was linkage disequilibrium (LD) among linked loci and unlinked loci pairs, and 1 955 out of 4 560 loci pairs (42.9%) had significant LD (< 0.01) with average′ value of 0.23. A total of 12 locus found by GLM method significantly at the level of<0.01 which explained 2.34%-13.94% of the phenotypic variance and the mean value was 6.33%. CAAS2050 (2=13.94%) and B153(2=11.36%) kept the max value. Meanwhile, 9 loci were found by MLM method significantly at the level of<0.01 which explained 2.80%-9.22% of the phenotypic variance and the mean value was 5.16%. P89(2=9.22%) and P3*(2=8.28%) kept the max value. A total of 7 loci were detected in common by GLM and MLM.Genetic diversity and population structure of 595 accessions were analyzed through SSR markers. In the two association analysis models, 12 markers were associated with nine traits including stem node number, plant height, peduncle length, diameter of main stem, panicle length, panicle diameter, primary branch number per panicle, spikelet number per primary branch, protein content by GLM. Nine markers were associated with eight traits including stem node number, peduncle length, leaf width, diameter of main stem, panicle diameter, primary branch number per panicle, spikelet number per primary branch, 1000-grain weight by MLM.

foxtail millet; Shanxi; landrace; SSR; genetic diversity; linkage disequilibrium; association analysis

10.3864/j.issn.0578-1752.2019.22.013

2019-04-22；

2019-06-20

山西省平臺(tái)基地和人才專項(xiàng)優(yōu)秀人才科技創(chuàng)新項(xiàng)目（201705D211026）、山西省農(nóng)業(yè)科學(xué)院特色農(nóng)業(yè)技術(shù)攻關(guān)（YGG17057）、國(guó)家谷子高粱產(chǎn)業(yè)技術(shù)體系（CARS-06-13.5-A16）

王海崗，E-mail：nkywhg@126.com。通信作者穆志新，E-mail：muzx2008@sina.com。通信作者喬治軍，E-mail：nkypzs@126.com

（責(zé)任編輯 李莉）