劉文萍，呂偉，黎冬華，任果香，張艷欣，文飛，韓俊梅，張秀榮

（1山西省農(nóng)業(yè)科學(xué)院經(jīng)濟作物研究所，太原 030031；2中國農(nóng)業(yè)科學(xué)院油料作物研究所，武漢 430062）

芝麻種質(zhì)資源成株期抗旱性關(guān)聯(lián)分析

劉文萍1，呂偉1，黎冬華2，任果香1，張艷欣2，文飛1，韓俊梅1，張秀榮2

（1山西省農(nóng)業(yè)科學(xué)院經(jīng)濟作物研究所，太原 030031；2中國農(nóng)業(yè)科學(xué)院油料作物研究所，武漢 430062）

【目的】利用33個多態(tài)性SSR分子標記分別與18個不同的抗旱性狀進行關(guān)聯(lián)分析，發(fā)掘與抗旱相關(guān)的主要基因位點，為抗旱基因定位和功能標記開發(fā)提供基礎(chǔ)；通過對100份芝麻種質(zhì)資源進行抗旱性鑒定，發(fā)掘優(yōu)異的耐旱種質(zhì)，為芝麻抗旱育種提供指導(dǎo)?！痉椒ā坎捎门柙院头磸?fù)干旱法，對芝麻種質(zhì)資源群體進行成株期抗旱性鑒定獲得表型指標測定值，利用SAS、SPSS和隸屬函數(shù)等進行統(tǒng)計分析，綜合評價其抗旱性，利用GLM模型和MLM模型，將表型數(shù)據(jù)與分子標記進行關(guān)聯(lián)分析?！窘Y(jié)果】研究群體干旱脅迫處理后，材料間響應(yīng)差異明顯，考察的表型性狀測定值均小于對照；干旱脅迫條件下18個性狀值的變異系數(shù)平均為0.31，高于對照（平均為0.19）；處理與對照間各性狀指標經(jīng)配對t檢驗，均達極顯著水平；通過連續(xù)變數(shù)的次數(shù)分布統(tǒng)計方法、主成分分析和隸屬函數(shù)分析，篩選出10個與抗旱性響應(yīng)關(guān)系密切的指標，并篩選出12份高抗旱種質(zhì)；基于芝麻基因組篩選出的33個多態(tài)性SSR標記掃描供試材料，共檢測到170個等位變異，平均每個標記5.15個；利用structure數(shù)學(xué)模型對供試群體進行遺傳結(jié)構(gòu)分析，可分為2個亞群；利用GLM模型和MLM模型分別檢測到120個和63個標記位點與供試群體抗旱系數(shù)顯著關(guān)聯(lián)（P＜0.05），表型變異解釋率分別為3.85%—14.30%和4.00%—12.5%，解釋率大于10%的標記位點分別有12個和3個，其中，位點4033-3和4033-2均與第一主成分因子第2次復(fù)水前萎蔫葉片數(shù)顯著關(guān)聯(lián)，且變異解釋率均為最高，分別達14.3%和12.5%，2個模型共同檢測到的標記位點有5個。通過引物序列在基因組上的位置比對，發(fā)現(xiàn)3個可能存在芝麻抗旱相關(guān)基因的基因組區(qū)段?！窘Y(jié)論】利用綜合評價方法，篩選出柳林芝麻3號、g80、8602-2等12份高抗旱芝麻種質(zhì)，同時利用GLM和MLM 2個模型檢測到與第2次復(fù)水前萎蔫葉片數(shù)顯著關(guān)聯(lián)的標記位點（位點4033-3和位點4033-2），且變異解釋率最高，分別達14.3%和12.5%。

芝麻種質(zhì)資源；抗旱性；SSR標記；關(guān)聯(lián)分析

0 引言

【研究意義】芝麻隸屬胡麻科，一年生草本植物，是中國主要的油料作物之一，在中國有著悠久的種植歷史，多數(shù)省份皆有種植[1]。山西地處干旱半干旱地區(qū)，降水量少，干旱嚴重，導(dǎo)致芝麻出苗率低、生長緩慢、產(chǎn)量降低，成為限制芝麻增產(chǎn)的主要限制因子，因此，研究芝麻的抗旱性，挖掘抗旱芝麻種質(zhì)資源，同時，利用SSR分子標記與抗旱性狀進行關(guān)聯(lián)分析，發(fā)掘與抗旱相關(guān)的主要基因位點，開發(fā)抗旱分子標記，并將其應(yīng)用到芝麻育種中，對提高芝麻的抗旱性和產(chǎn)量水平具有很重要意義?！厩叭搜芯窟M展】國內(nèi)外已對玉米、水稻、大豆、小麥、花生、油菜[2-7]等作物進行了較系統(tǒng)的抗旱性研究。早在20世紀70年代末，CARTER等[8]通過植株萎蔫程度對300余份大豆種質(zhì)進行了抗旱性鑒定；SONGSRI等[9]、劉吉利等[10]、張智猛等[11]對花生的生長發(fā)育和生理生態(tài)等方面的干旱適應(yīng)機制進行了研究；NOROUZI等[12]對甘藍型油菜研究表明，干旱脅迫下葉片含水量[13]、萎蔫指數(shù)[14]、種子活力指數(shù)[7]、相對根體積[15]以及脯氨酸[16]等生理指標發(fā)生明顯變化，千粒重、生物量等下降[17]，可作為油菜抗旱性鑒定的指標。在芝麻抗旱性研究方面國內(nèi)外也有一些報道，MENSAH等[18]利用盆栽法對芝麻苗期進行間隔不同天數(shù)澆水處理研究，結(jié)果表明，

水分不足會直接影響芝麻植株的生長，株高、葉面積、地下部干重和地上部干重等降低，導(dǎo)致單株產(chǎn)量顯著下降；孫建等[19]利用盆栽試驗對芝麻苗期進行干旱脅迫研究，結(jié)果表明，芝麻苗期進行干旱處理對千粒重和產(chǎn)量的影響較大。【本研究切入點】作物的抗旱性屬于較復(fù)雜的數(shù)量性狀，受多基因控制[20]，關(guān)聯(lián)分析是研究復(fù)雜數(shù)量性狀的有效方法[21]，中國芝麻抗旱性分子基礎(chǔ)研究報道較少，僅見黎冬華等[22]采用不同濃度PEG 6000脅迫發(fā)芽芝麻種子基于SSR、SRAP和AFLP分子標記的關(guān)聯(lián)分析研究報道，而國外對芝麻抗旱性分子基礎(chǔ)研究還尚未見報道。【擬解決的關(guān)鍵問題】本研究通過對芝麻自然研究群體進行抗旱表型鑒定，獲得抗旱相關(guān)表型數(shù)據(jù)，并與全基因組分子標記數(shù)據(jù)進行關(guān)聯(lián)分析，發(fā)掘與抗旱相關(guān)性狀基因位點，篩選得到抗旱種質(zhì)資源，為芝麻抗旱遺傳改良提供理論基礎(chǔ)。

1 材料與方法

1.1 試驗材料

選用來源于山西和陜西芝麻種質(zhì)資源共計100份，由國家芝麻種質(zhì)資源中期庫（中國農(nóng)業(yè)科學(xué)院油料作物研究所）提供，并進行了4代純化。

1.2 干旱處理

試驗于2014年5月至10月在山西省農(nóng)業(yè)科學(xué)院經(jīng)濟作物研究所人工遮雨棚內(nèi)進行，采用盆栽法[23]，每個材料種植6盆，其中3盆用于干旱處理，3盆為對照，共計600盆。盆缽口徑35 cm，盆栽土用芬蘭進口營養(yǎng)土與蛭石按1∶1等量混合，每盆裝3 kg，每盆澆等量的水1 L，待適墑時播種。出苗后間苗2次，每盆保留均勻生長的苗5株，在干旱處理前保持正常供水。在芝麻4對真葉展開時進行干旱處理，采用反復(fù)干旱法，處理3次，每次干旱處理后待50%左右的材料出現(xiàn)永久萎蔫時進行復(fù)水，之后全部材料按正常管理澆水直至成熟。

1.3 表型性狀調(diào)查

處理期調(diào)查：在干旱處理的每次復(fù)水前同一時間段內(nèi)（3次復(fù)水前均在上午9：00以前）調(diào)查每株總?cè)~片數(shù)（全展葉）、萎焉葉片數(shù)（全展葉）；每次復(fù)水后第2天早上調(diào)查每株總?cè)~片數(shù)（全展葉）、萎焉葉片數(shù)（全展葉）。成熟期調(diào)查：調(diào)查每個材料處理和對照每株的株高、始蒴高度、單株蒴果數(shù)、單株產(chǎn)量、千粒重、單株地上部分生物量、地下部分生物量，計算平均值。

1.4 抗旱能力的綜合評價

各指標抗旱系數(shù)Xj=處理測定值/對照測定值；負向指標的抗旱系數(shù)Xj=1-處理測定值/對照測定值（1）

每一個材料各綜合指標的隸屬函數(shù)值由公式（2）求得[24-25]，式中，Xj表示第j個綜合指標，Xmax和Xmin分別表示第j個綜合指標的最大值和最小值。

由公式（3）計算各綜合指標的權(quán)重[24-25]，式中，Wj表示第j個綜合指標在所有綜合指標中的重要程度；Pj為各材料第j個綜合指標的貢獻率。

由公式（4）計算各材料的綜合抗旱能力（D值）。

用SPSS和SAS 9.1軟件進行相關(guān)性分析、主成分分析。

1.5 基因組DNA的提取和SSR分子標記分析

每份材料取苗期嫩葉的混合樣，采用CTAB法[26]提取基因組DNA，經(jīng)紫外分光光度計法檢測其質(zhì)量和濃度，稀釋至20 ng·μL-1，置-20℃冰箱保存。PCR體系及擴增程序參照文獻[27]。PCR擴增產(chǎn)物經(jīng)6%聚丙烯酰胺凝膠電泳，銀染顯色后照相。

1.6 關(guān)聯(lián)分析方法數(shù)據(jù)統(tǒng)計分析

以Structure 2.3.1軟件進行群體遺傳結(jié)構(gòu)分析，估計最佳群體組群數(shù)K，其取值范圍為1—10，將MCMC（Markov Chain Monte Carlo）開始時的不作數(shù)迭代（length of burn-in period）設(shè)為100 000次，再將不作數(shù)迭代后的MCMC設(shè)為100 000次，迭代次數(shù)（number of iterations）設(shè)置為5，依據(jù)似然值最大原則，選取合適的K值為群體數(shù)目，并繪出基于模型的群體遺傳結(jié)構(gòu)圖。采用Tassel 2.1 軟件GLM（general linear model）和MLM（mixed linear model）2種模型程序，進行分子標記數(shù)據(jù)和表型數(shù)據(jù)的關(guān)聯(lián)分析。

2 結(jié)果

2.1 供試材料對干旱脅迫的響應(yīng)分析

對供試材料考察了18個性狀，包括總?cè)~片數(shù)、萎焉葉片數(shù)、生物量、單株蒴果數(shù)、單株產(chǎn)量、株高等，處理與對照間存在明顯差異，各項指標平均值均小于對照，表明各個材料均受到干旱脅迫的影響。干旱脅迫下，18個性狀值的變異系數(shù)在0.07—0.65，平均為0.31，對照的變異系數(shù)在0.06—0.40，平均為0.19，除株高外，其他性狀干旱脅迫的變異系數(shù)均大于對照。處理與對照間各性狀指標經(jīng)配對t檢驗，均為極顯著，表明干旱處理對各性狀產(chǎn)生了明顯的抑制作用（表1）。處理組與對照組相比18個性狀的平均抗旱系數(shù)在0.326—0.951，均小于1.000，也表明測定的各項指標均受到干旱脅迫不同程度的影響（表2）。以上結(jié)果表明，研究材料間和調(diào)查的18個性狀均受到干旱脅迫，對干旱脅迫的響應(yīng)存在明顯差異。

2.2 各項指標抗旱性綜合評價

參照連續(xù)變數(shù)的次數(shù)分布統(tǒng)計方法，將供試材料18個性狀的抗旱系數(shù)以組距為0.2分成7個組區(qū)間，制作成次數(shù)分布表（表3）。結(jié)果表明，在同一組區(qū)間各性狀的抗旱系數(shù)分布次數(shù)相差較大，其中，第1次復(fù)水前總?cè)~片數(shù)、第2次復(fù)水前總?cè)~片數(shù)、千粒重這三個指標的抗旱系數(shù)在0.9＜Xj≤1.1區(qū)間的分布頻率分別為82%、50%和53%，表明這三個性狀對干旱脅迫反應(yīng)遲鈍；而第1次復(fù)水前萎蔫葉片數(shù)、第2次復(fù)水前萎蔫葉片數(shù)、第2次復(fù)水后萎蔫葉片數(shù)、第3次復(fù)水前總?cè)~片數(shù)、第3次復(fù)水前萎蔫葉片數(shù)、第3次復(fù)水后萎蔫葉片數(shù)、地上部鮮重、地上部干重、成熟期株高、單株產(chǎn)量這10個指標的抗旱系數(shù)在Xj＜0.9區(qū)間分布頻率分別為92%、77%、67%、86%、100%、83%、62%、60%、68%和62%，表明這10個性狀對干旱脅迫的反應(yīng)較敏感。由抗旱系數(shù)相關(guān)系數(shù)矩陣可看出（表4），各單項指標間存在一定的相關(guān)性，說明芝麻抗旱性是一個復(fù)雜性狀。

表1 各性測定值的差異分析Table 1 Variance analysis of mred value between characters

表2 各單指標的抗旱系數(shù)Table 2 Drought resistance coefficients of single index

續(xù)表2 Continued table 2

續(xù)表2 Continued table 2

續(xù)表2 Continued table 2

表3 供試材料各性狀指標的抗旱系數(shù)在不同區(qū)間的分布Table 3 Different distributions of drought-resistance coefficients of tested materials traits indexes

對干旱較敏感的10個性狀抗旱系數(shù)進行主成分分析，結(jié)果表明，前5個因子的累積貢獻率達到82%，因而取前5個因子，將以上10個性狀轉(zhuǎn)換為5個新的獨立的綜合指標，分別用CI（1）、CI（2）、CI（3）、CI（4）、CI（5）表示。第一主成分中第2次復(fù)水前萎蔫葉片數(shù)、第3次復(fù)水后萎蔫葉片數(shù)有較強的載荷量；第二主成分中地上部干重有較強的載荷量；第三主成分中，第1次復(fù)水前萎蔫葉片數(shù)、第3次復(fù)水前總?cè)~片數(shù)、單株產(chǎn)量有較強的載荷量；第四主成分中，地上部鮮重、株高有較強的載荷量；第五主成分中，第2次復(fù)水后萎蔫葉片數(shù)、第3次復(fù)水前萎蔫葉片數(shù)有較強的載荷量（表5）。

根據(jù)公式（2）、（3）、（4）求得供試材料抗旱性綜合評價D值（表2）。根據(jù)D值大小對供試材料進行抗旱性排序，篩選出12份高抗旱芝麻種質(zhì)資源（D值＞0.65），占供試材料的12%，其中，包括柳林芝麻3號、吳堡縣岔上鄉(xiāng)步墕村芝麻、g80、8602-2、柳林芝麻1號、臨縣芝麻2號、臨縣芝麻5號、g62、g45、吳堡縣岔上鄉(xiāng)郭家墕村芝麻、2012-48-01、四棱芝麻。

2.3 關(guān)聯(lián)分析

基于芝麻基因組開發(fā)多態(tài)性較好的33個SSR標記在100份芝麻供試材料中共檢測到170個等位變異，平均每個標記5.15個。利用structure數(shù)學(xué)模型對供試群體進行遺傳結(jié)構(gòu)分析，由此判斷供試群體可被分為2個亞群，并繪制供試材料群體結(jié)構(gòu)圖（圖1）。

利用GLM模型檢測到120個位點與以上對干旱脅迫響應(yīng)敏感的10個性狀顯著關(guān)聯(lián)（P＜0.05），每個性狀的位點數(shù)在4—31個，變異解釋率為3.85%—14.30%（表6），有28個位點同時與2個或多個性狀關(guān)聯(lián)；位點2177-3關(guān)聯(lián)的頻次最高，同時與5個指標顯著關(guān)聯(lián)，分別為第1次復(fù)水前萎蔫葉片數(shù)、第2、3次復(fù)水前、后萎蔫葉片數(shù)，解釋率為6.43%—12.76%，其中第2次復(fù)水前萎蔫葉片數(shù)、第3次復(fù)水后萎蔫葉片數(shù)均屬于第一主成分中的因子；解釋率大于10%的標記位點有12個，其中1816-2、1664-3、1550-5、2177-3、1816-4、4033-3等6個位點與2個第一主成分因子（第2次復(fù)水前萎蔫葉片數(shù)和第3次復(fù)水后萎蔫葉片數(shù)）顯著關(guān)聯(lián)，4033-3是所有位點中變異解釋率最高的，對第2次復(fù)水前萎蔫葉片數(shù)解釋率達14.30%。

表4 各指標抗旱系數(shù)的相關(guān)系數(shù)矩陣Table 4 Correlation coefficient matrix between drought resistance indexes

表5 主成分載荷矩陣及貢獻率Table 5 Component matrix and contribution rate

圖1 基于SSR標記的100份芝麻種質(zhì)遺傳結(jié)構(gòu)Fig. 1 Population structure of 100 sesame materials based on SSR marker

利用MLM模型檢測到63個與10個性狀顯著關(guān)聯(lián)（P＜0.05），每個性狀的位點數(shù)分別在3—11個，變異解釋率為4.00%—12.50%，有13個位點同時與2個或多個性狀關(guān)聯(lián)；位點1816-3同時與4個指標顯著關(guān)聯(lián)，分別為第1、2次復(fù)水前萎蔫葉片數(shù)和第2、3次復(fù)水后萎蔫葉片數(shù)，解釋率為4.32%—6.32%，其中第2次復(fù)水前萎蔫葉片數(shù)、第3次復(fù)水后萎蔫葉片數(shù)均屬于第一主成分中的因子；解釋率大于10%的標記位點有3個，分別是2209-7、1935-6和4033-2，其中，變異解釋率最高的位點4033-2所關(guān)聯(lián)的性狀也為第2次復(fù)水前萎蔫葉片數(shù)，解釋率達12.5%。

利用2種模型均能檢測到的標記位點有5個，其中位點4033-6、1664-1同時與第1次復(fù)水前萎蔫葉片數(shù)顯著關(guān)聯(lián)，解釋率分別為4.46%—5.88%和4.07%—4.64%，該指標屬于第三主成分中的因子；位點2218-1與第2次復(fù)水后萎蔫葉片數(shù)關(guān)聯(lián)顯著，解釋率為9.9%—13.66%，該指標屬于第五主成分中的因子；位點1550-3同時與地上部鮮重、地上部干重顯著關(guān)聯(lián)，解釋率分別為5.36%—6.08%和4.09%—7.88%，這兩個指標分別屬于第四、二主成分中的因子。

表6 與抗旱性顯著相關(guān)（P＜0.05）的標記位點數(shù)和變異解釋率Table 6 Number of loci associated with drought resistance（P＜0.05）and phenotypic variation explanation rates

對第一主成分因子性狀關(guān)聯(lián)到的位點進一步分析，其中貢獻率較大且重復(fù)率高的標記位點主要有：4033-2和4033-3對第2次復(fù)水前萎蔫葉片數(shù)變異解釋率分別為12.5%和14.3%，2177-3對第2次復(fù)水前萎蔫葉片數(shù)、第3次復(fù)水后萎蔫葉片數(shù)變異解釋率分別為11.86%和12.21%，1816-2和1826-4對第3次復(fù)水后萎蔫葉片數(shù)變異解釋率分別為10.44%和13.99%，這些標記位點分別來自4033、2177和1816 3個引物，通過比對這三個引物序列在芝麻基因組上的位置，結(jié)果表明，4033位于LG4的1 981 965 bp附近，2177位于LG6的11 802 028 bp附近，1816位于LG11的10 438 634 bp附近，推測以上3個基因組區(qū)段可能存在芝麻抗旱性相關(guān)基因。

3 討論

近年來，中國國內(nèi)外對芝麻抗旱性研究多局限于大田干旱脅迫和利用PEG模擬干旱對芝麻種子芽期進行抗旱性鑒定，本研究對芝麻成株期進行反復(fù)干旱試驗，綜合利用與抗旱響應(yīng)關(guān)系密切的多項指標進行芝麻抗旱性綜合評價，為芝麻抗旱性精準鑒定提供了方法參考和理論基礎(chǔ)。截至目前，國內(nèi)外對芝麻抗旱性研究的供試材料數(shù)量和指標數(shù)均少于本研究。KADKHODAIE等[28]評價了10份芝麻品種在不同土壤水分條件下的抗旱性差異；孫建等[19]對5個芝麻品種進行苗期干旱脅迫研究，測定指標均少于10個。本研究對100份芝麻資源測定了18個相關(guān)表型性狀，其中，總?cè)~片數(shù)和萎蔫葉片數(shù)屬首次采用。采用連續(xù)變數(shù)的次數(shù)分布統(tǒng)計方法對供試材料進行分析發(fā)現(xiàn)，千粒重對干旱脅迫反應(yīng)較遲鈍，這與孫建等[19]等的研究結(jié)果不一致，而與祁旭升等[25]在胡麻抗旱性研究所得結(jié)果相一致；利用綜合分析方法得到第2次復(fù)水前萎蔫葉片數(shù)和第3次復(fù)水后萎蔫葉片數(shù)2個與抗旱響應(yīng)最為密切，為芝麻抗旱鑒定提供了可靠的新的指標。

植物的抗旱性是受多因素影響的、且非常復(fù)雜的數(shù)量性狀，不同指標的選擇不同，研究得到的結(jié)論也不一樣。因此，單一指標或少數(shù)指標的簡單分析很難且可靠地反應(yīng)不同資源品種的抗旱性。近年來，雖然國內(nèi)外學(xué)者針對不同作物從不同角度研究和提出了更加全面的抗旱評價鑒定方法[29-30]，然而現(xiàn)今對抗旱性鑒定評價的指標和方法均沒有準確可靠的統(tǒng)一規(guī)定，因此，建立科學(xué)的抗旱性綜合評價方法是芝麻抗旱性鑒定的基礎(chǔ)。目前普遍認為多指標多方法相結(jié)合的抗旱性綜合評價比較可靠[31]，朱宗河等[32]對油菜抗旱性的研究表明，利用隸屬函數(shù)等分析對甘藍型油菜耐旱性進行綜合評價，可篩選甘藍型油菜耐旱種質(zhì)。因此，本研究采用綜合評價方法，可以較好地揭示各指標與抗旱性的關(guān)系，這樣既考慮了各指標間的相互關(guān)系，又考慮到各指標的重要性。該方法已在大豆[33]、小麥[34]等作物上應(yīng)用，取得了比較理想的結(jié)果。不僅可以縮短芝麻抗旱性鑒定周期，同時為下一步研究耐旱分子機制提供理論依據(jù)。

芝麻的抗旱性是復(fù)雜的數(shù)量性狀，關(guān)聯(lián)分析是研究復(fù)雜數(shù)量性狀的有效方法，國內(nèi)外在芝麻抗旱分子基礎(chǔ)研究方面，僅見黎冬華等[22]采用PEG 6000脅迫發(fā)芽芝麻種子基于SSR、SRAP和AFLP分子標記的關(guān)聯(lián)分析研究報道，利用GLM模型獲得2個解釋率（4.965%和4.18%）較高的SRAP標記和1個解釋率（4.38%）較高SSR標記。本研究對芝麻成株期進行抗旱鑒定，利用全基因組開發(fā)的33個多態(tài)性SSR標記與抗旱相關(guān)性狀進行關(guān)聯(lián)分析，利用GLM和MLM 2種模型分別檢測出120個和63個與抗旱性有顯著關(guān)聯(lián)（P＜0.05）的分子標記位點，變異解釋率分別為3.85%—14.30%和4.00%—12.5%，其中，2種模型均檢測到與第一主成分因子第2次復(fù)水前萎蔫葉片數(shù)顯著關(guān)聯(lián)的標記位點，且變異解釋率均為最高，分別達14.3%（位點4033-3）和12.5%（位點4033-2），通過比對4033引物序列在芝麻基因組上的位置，進一步發(fā)現(xiàn)基因SIN_1006195在其附近，該基因?qū)儆谝褕蟮琅c抗旱緊密相關(guān)的NAC基因家族，可以作為芝麻抗旱后續(xù)研究的候選基因之一。研究結(jié)果為芝麻抗旱功能分子標記開發(fā)和抗旱基因發(fā)掘奠定了重要基礎(chǔ)。

本研究通過綜合評價，還篩選出一批高抗旱種質(zhì)材料，為抗旱育種研究提供了關(guān)鍵材料。

4 結(jié)論

采用抗旱系數(shù)進行綜合評價，獲得柳林芝麻3號、g80、8602-2等高抗旱芝麻種質(zhì)12份，可用于抗旱性育種；利用GLM和MLM 2個模型均檢測到與第一主成分因子第2次復(fù)水前萎蔫葉片數(shù)顯著關(guān)聯(lián)且變異解釋率最高的標記位點（位點4033-3和位點4033-2），此標記可應(yīng)用于分子標記輔助選擇。

[1] 黃鳳洪. 花生芝麻加工技術(shù). 北京: 金盾出版社, 1995.

HUANG F H.The Processing Technology of Peanut and Sesame. Beijing: Shield Press, 1995. (in Chinese)

[2] NEMETH M, JANDA T, HORVATH E, PALDI E, SZALAI G. Exogenous salicylic acid increases polyamine content but may decrease drought tolerance in maize.Plant Science, 2002, 162: 569-574.

[3] FU B Y, XIONG J H, ZHU L H, ZHAO X Q, XU H X, GAO Y M, LI Y S, XU J L, LI Z K. Identification of functional candidate genes for drought tolerance in rice.Molecular Genetics and Genomics, 2007, 278: 599-609.

[4] 李燦東, 蔣洪蔚, 郭泰, 王志新, 吳秀紅, 鄭偉, 陳慶山, 胡國華.大豆耐旱選擇群體基因型分析與株高QTL定位. 大豆科學(xué), 2011, 30(1): 15-19.

LI C D, JIANG H W, GUO T, WANG Z X, WU X H, ZHENG W, CHEN Q S, HU G H. QTL identification of plant height and analysis of genotype to soybean in selection population.Soybean Science, 2011, 30(1): 15-19. (in Chinese)

[5] PAKNIYAT H, TAVAKOL E. RAPD markers associated with drought tolerance in bread wheat (Triticum aestivumL.).Pakistan Journal of Biological Sciences, 2007, 10: 3237-3239.

[6] 賀鴻雁, 孫存華, 杜偉, 李揚. PEG 6000脅迫對花生幼苗滲透調(diào)節(jié)物質(zhì)的影響. 中國油料作物學(xué)報, 2006, 28(1): 76-78.

HE H Y, SUN C H, DU W, LI Y. Effects of PEG 6000 osmotic stress on osmolytes of peanut seed-ling.Chinese Journal of Oil Crop Sciences, 2006, 28(1): 76-78. (in Chinese)

[7] 李震, 楊春杰, 張學(xué)昆, 鄒崇順, 程勇, 鄭普英, 李桂英. PEG脅迫下甘藍型油菜品種(系)種子發(fā)芽耐旱性鑒定. 中國油料作物學(xué)報, 2008, 30(4): 438-442.

LI Z, YANG C J, ZHANG X K, ZOU C S, CHENG Y, ZHENG P Y, LI G Y. Evaluation of drought tolerance in rapeseed (Brassica napusL.) during germination under PEG6000 stress.Chinese Journal of Oil Crop Sciences, 2008, 30(4): 438-442. (in Chinese)

[8] CARTER T E JR, RUFTY T W, KUO C G. Soybean plant introduction exhibiting drought and aluminum tolerance. Adaptation of food crops to temperature and water stress. Proceedings of an international symposium, Taiwan, 1993: 335-346.

[9] SONGSRI P, JOGLOY S, KESMALA T, VORASOOT N, AKKASAENG C, PATANOTHAI A. Responses of reproductive characters of drought resistant peanut genotypes to drought.Asian Journal of Plant Sciences, 2008, 7(5): 427-439.

[10] 劉吉利, 趙長星, 吳娜, 王月福, 王銘倫. 苗期干旱及復(fù)水對花生光合特性及水分利用效率的影響. 中國農(nóng)業(yè)科學(xué), 2011, 44(3): 469-476.

LIU J L, ZHAO C X, WU N, WANG Y F, WANG M L. Effects of drought and rewatering at seedling stage on photosynthetic characteristics and water use efficiency of peanut.Scientia Agricultura Sinica, 2011, 44(3): 469-476. (in Chinese)

[11] 張智猛, 戴良香, 丁紅, 陳殿緒, 楊偉強, 宋文武, 萬書波. 中國北方主栽花生品種抗旱性鑒定與評價. 作物學(xué)報, 2012, 38(3): 495-504.

ZHANG Z M, DAI L X, DING H, CHEN D X, YANG W Q, SONG W W, WAN S B. Identification and evaluation of drought resistance in different peanut varieties widely grown in Northern China.Acta Agronomica Sinica, 2012, 38(3): 495-504. (in Chinese)

[12] NOROUZI M, TOORCHI M, SALEKDEH G H, MOHAMMADI S A, NISHABOURI M R, AHARIZAD S. Effect of water deficit on growth, grain yield and osmotic adjustment in rapeseed.Journal of Food, Agriculture & Environment, 2008, 6(2): 132-138.

[13] ESLAM B P. Evaluation of physiological indices, yield and its components as screening techniques for water deficit tolerance in oilseed rape cultivars.Journal of Agricultural Science and Technology, 2009, 11: 413-422.

[14] 李震, 吳北京, 陸光遠, 程勇, 鄒崇順, 張學(xué)昆. 不同基因型油菜對苗期水分脅迫的生理響應(yīng). 中國油料作物學(xué)報, 2012, 34(1): 33-39.

LI Z, WU B J, LU G Y, CHENG Y, ZOU C S, ZHANG X K. Differences in physiological responses ofBrassica napusgenotypes under water stress during seedling stage.Chinese Journal of Oil Crop Science, 2012, 34(1): 33-39. (in Chinese)

[15] 張霞, 謝小玉. PEG脅迫下甘藍型油菜種子萌發(fā)期抗旱鑒定指標的研究. 西北農(nóng)業(yè)學(xué)報, 2012, 21(2): 72-77.

ZHANG X, XIE X Y. Studies on identification indexes of drought resistance by PEG during seed germination of rapeseed (Brassica napusL.).Acta Agriculturae Boreali-occidentalis Sinica, 2012, 21(2): 72-77. (in Chinese)

[16] 謝小玉, 張霞, 張兵. 油菜苗期抗旱性評價及抗旱相關(guān)指標變化分析. 中國農(nóng)業(yè)科學(xué), 2013, 46(3): 476-485.

XIE X Y, ZHANG X, ZHANG B. Evaluation of drought resistance and analysis of variation of relevant parameters at seedling stage of rapeseed (Brassica napusL.).Scientia Agricultura Sinica, 2013, 46(3): 476-485. (in Chinese)

[17] ABBASIAN A, RAD A H S. Investigation the response of rapeseed cultivars to moisture regimes in different growth stages.Journal of Central European Agriculture, 2011, 12(2): 353-366.

[18] MENSAH J K, OBADONI B O, ERUOTOR P G, ONOMEIRIEGUNA F. Simulated flooding and drought effects on germination, growth and yield parameters of sesame (Sesamum indicumL.).African Journal of Biotechnology, 2006, 5(13): 1249-1253.

[19] 孫建, 饒月亮, 樂美旺, 顏廷獻, 顏小文, 周紅英. 干旱脅迫對芝麻生長與產(chǎn)量性狀的影響及其抗旱性綜合評價. 中國油料作物學(xué)報, 2010, 32(4): 525-533.

SUN J, RAO Y L, LE M W, YAN T X, YAN X W, ZHOU H Y. Effects of drought stress on sesame growth and yield characteristics and comprehensive evaluation of drought tolerance.Chinese Journal of Oil Crop Sciences, 2010, 32(4): 525-533. (in Chinese)

[20] 唐富福. 水稻農(nóng)藝和品質(zhì)若干性狀的全基因組初步關(guān)聯(lián)分析[D].杭州: 浙江大學(xué), 2013.

TANG F F. Genome-wide association analysis in some agronomic and quality traits of rice[D]. Hangzhou: Zhejiang University, 2013. (in Chinese)

[21] 薛瑞, 周廣奇, 胡新文, 郭建春, 唐燕瓊. 柱花草種質(zhì)抗旱性綜合評價. 中國農(nóng)學(xué)通報, 2009, 25(11): 224-233.

XUE R, ZHOU G Q, HU X W, GUO J C, TANG Y Q. Comprehensive evaluation of drought resistance in stylosanthes germplasms.Chinese Agricultural Science Bulletin, 2009, 25(11): 224-233. (in Chinese)

[22] 黎冬華, 劉文萍, 張艷欣, 王林海, 危文亮, 高媛, 丁霞, 王蕾, 張秀榮. 芝麻耐旱性的鑒定方法及關(guān)聯(lián)分析. 作物學(xué)報, 2013, 39(8): 1425-1433.

LI D H, LIU W P, ZHANG Y X, WANG L H, WEI W L, GAO Y, DING X, WANG L, ZHANG X R. Identification method of drought tolerance and association mapping for sesame (Sesamum indicumL.).Acta Agronomica Sinica, 2013, 39(8): 1425-1433. (in Chinese)

[23] 孫建, 張秀榮, 張艷欣, 王林海, 黃波. 濕害處理對不同生育時期芝麻葉片保護酶活性和種子產(chǎn)量的影響. 應(yīng)用與環(huán)境生物學(xué)報, 2009, 15(6): 790-795.

SUN J, ZHANG X R, ZHANG Y X, WANG L H, HUANG B. Effects of waterlogging on leaf protective enzyme activities and seed yield of sesame at different growth stages.Chinese Journal of Applied & Environmental Biology, 2009, 15(6): 790-795. (in Chinese)

[24] 劉海卿, 孫萬倉, 劉自剛, 武軍艷. 北方寒旱區(qū)白菜型冬油菜抗寒性與抗旱性評價及其關(guān)系. 中國農(nóng)業(yè)科學(xué), 2015, 48(18): 3743-3756.

LIU H Q, SUN W C, LIU Z G, WU J Y. Evaluation of drought resistance and cold resistance and research of their relationship at seedling stage of winter rapeseed (Brassica campestrisL.) in cold and arid regions in North China.Scientia Agricultura Sinica, 2015, 48(18): 3743-3756. (in Chinese)

[25] 祁旭升, 王興榮, 許軍, 張建平, 米君. 胡麻種質(zhì)資源成株期抗旱性評價. 中國農(nóng)業(yè)科學(xué), 2010, 43(15): 3076-3087.

QI X S, WANG X R, XU J, ZHANG J P, MI J. Drought-resistance evaluation of flax germplasm at adult plant stage.Scientia Agricultura Sinica, 2010, 43(15): 3076-3087. (in Chinese)

[26] POREBSKI S, BAILEY L G, BAUM B R. Modification of a CTABDNA extraction protocol for plants containing high polysaccharide and polyphenol components.Plant Molecular Biology Reporter, 1997, 15: 8-15.

[27] ZHANG Y X, ZHANG X Y, HUA W, WANG L H, CHE Z. Analysis of genetic diversity among indigenous landraces from sesame (Sesamum indicumL.) core collection in China as revealed by SRAP and SSR markers.Genes and Genomics, 2010, 32: 207-215.

[28] KADKHODAIE A, ZAHEDI M, RAZMJOO J, PESSARAKLI M. Changes in some anti-oxidative enzymes and physiological indices among sesame genotypes (Sesamum indicumL.) in response to soil water deficits under field conditions.Acta Physiology Plant, 2014, 36: 641-650.

[29] CONDON A G, RICHARDS R A, REBETZKE G J, FARQUHAR G D. Breeding for high water use efficiency.Journal of Experimental Botany, 2004, 55: 2447-2460.

[30] SUBRAHMANYAM D, SUBASH N, HARIS A, SIKKA A K. Influence of water stress on leaf photosynthetic characteristics in wheat cultivars differing in their susceptibility to drought.Photosynthetica, 2006, 44: 125-129.

[31] SUBRAHMANYAM D, SUBASH N, HARIS A, SIKKA A K. Influence of water stress on leaf photosynthetic characteristics in wheat cultivars differing in their susceptibility to drought.Photosynthetica, 2006, 44: 125-129.

[32] 朱宗河, 鄭文寅, 張學(xué)昆. 甘藍型油菜耐旱相關(guān)性狀的主成分分析及綜合評價. 中國農(nóng)業(yè)科學(xué), 2011, 44(9): 1775-1787.

ZHU Z H, ZHENG W Y, ZHANG X K. Principal component analysis and comprehensive evaluation on morphological and agronomic traits of drought tolerance in rapeseed (Brassica napusL.).Scientia Agricultura Sinica, 2011, 44(9): 1775-1787. (in Chinese)

[33] 楊萬明, 王敏, 李貴全, 杜維俊. PEG脅迫下大豆BIL群體芽期性狀與耐旱性評價. 中國油料作物學(xué)報, 2013, 35(5): 564-571.

YANG W M, WANG M, LI G Q, DU W J. Germination traits of soybean BIL population under PEG stress and assessment of their drought toleranceChinese Journal of Oil Crop Sciences, 2013, 35(5): 564-571. (in Chinese)

[34] 賈壽山, 朱俊剛, 王曙光, 史雨剛, 孫黛珍. 山西小麥地方品種萌發(fā)期的抗旱性. 華北農(nóng)學(xué)報, 2011, 26(2): 213-217.

JIA S S, ZHU J G, WANG S G, SHI Y G, SUN D Z. Drought resistance of Shanxi wheat landraces at sprouting stage.Acta Agriculturae Boreali-Sinica, 2011, 26(2): 213-217. (in Chinese)

（責任編輯 李莉）

Drought Resistance of Sesame Germplasm Resources and Association Analysis at Adult Stage

LIU WenPing1, Lü Wei1, LI DongHua2, REN GuoXiang1, ZHANG YanXin2, WEN Fei1, HAN JunMei1, ZHANG XiuRong2
(1Institute of Economic Crops Research, Shanxi Academy of Agricultural Sciences, Taiyuan 030031;2Institute of Oil Crops Research, Chinese Academy of Agricultural Sciences, Wuhan 430062)

【Objective】By association mapping between 33 SSR markers and 18 sesame drought-resistance phenotypes, identifying the associated loci, providing a basis for drought-resistance gene mapping and functional marker development, identifying drought-resistant sesame germplasm resources and providing parental materials and varieties for drought-resistance breeding.【Method】Sesame germplasm populations were cultivated using pots and stressed by drought-water-drought duplication treatment, morphology and growth index during the drought-stress period and mature period were analyzed, then the combinationmethod of SAS, SPSS, and membership function were used for drought-resistance evaluation. At the same time, association analysis was carried out using GLM and MLM models based on molecular markers and drought-resistance phenotype data.【Result】After drought stress, obviously different responses were observed among materials in the population, each index that related to drought resistance averaged smaller than the control; the averaged coefficient of variation of 18 traits was 0.31, which was higher than the control. T test of all traits between treatment and control showed a significant difference. A total of 10 indices related to drought-resistance were screened out, and 12 germplasms with high drought-resistance were identified by functional and clustering analysis. A total of 170 loci were detected by 33 SSR makers in the tested materials, the average loci number was 5.15 per marker. Genetic structure analysis indicated that the population consists of 2 sub-groups. There were 120 loci and 63 loci that were detected by association analysis based on GLM model and MLM model, respectively, the explanation rate ranged from 3.85% to 14.3%, and from 4.00% to 12.5%. There were 12 loci and 3 loci were detected with higher explanation rate than 10% by GLM and MLM, respectively, 5 loci were detected simultaneously by GLM and MLM. By comparison of the primers sequences with sesame genome, 3 genome regions were presumed located genes related with drought-resistance.【Conclusion】This research selected 12 sesame germplasms with high drought-resistance by comprehensive evaluation method, at the same time, 4033-3 and 4033-2 loci that were detected with the highest explanation rate by association analysis based on GLM model and MLM model, respectively, and the explanation rate were 14.3% and 12.5%.

sesame germplasm; drought resistance; SSR markers; association analysis

10.3864/j.issn.0578-1752.2017.04.003

2016-08-25；接受日期：2016-10-24

國家現(xiàn)代農(nóng)業(yè)芝麻產(chǎn)業(yè)技術(shù)體系（CARS-15）、山西省科技攻關(guān)項目（20140311010-3）、山西省農(nóng)業(yè)科學(xué)院科技自主創(chuàng)新能力提升工程項目（2016zzcx-07）、芝麻高產(chǎn)、優(yōu)質(zhì)、抗逆新品種選育（16 yzgc 009）

聯(lián)系方式：劉文萍，Tel：13593405471；E-mail：wenggeping@163.com