耿慶河 王蘭芬 武 晶 王述民

中國農(nóng)業(yè)科學(xué)院作物科學(xué)研究所,北京 100081

普通菜豆籽粒大小與形狀的QTL定位

耿慶河 王蘭芬 武 晶 王述民*

中國農(nóng)業(yè)科學(xué)院作物科學(xué)研究所,北京 100081

普通菜豆是世界上最重要的食用豆類作物之一,其籽粒大小和形狀與產(chǎn)量及外觀品質(zhì)密切相關(guān)。本研究以來自安第斯基因庫的大粒品種龍270709和來自中美基因庫的小粒品種F5910配置雜交組合,獲得的F2分離群體分別在哈爾濱大田與北京昌平溫室種植,對百粒重、粒長、粒寬、粒厚、長寬比和長厚比6個籽粒性狀進行了相關(guān)性分析和QTL定位。相關(guān)性分析表明,百粒重與粒長、粒寬、粒厚、長寬比、長厚比5個衡量籽粒大小和形狀的性狀均顯著正相關(guān)。利用基于完備區(qū)間作圖方法的IciMapping 4.1進行QTL定位,哈爾濱環(huán)境下定位到38個與百粒重、粒長、粒寬、粒厚、長寬比、長厚比相關(guān)的QTL,表型貢獻率介于2.39%~17.37%之間,分布在除第1染色體外的其余10條染色體上;北京昌平環(huán)境下定位到21個上述性狀的QTL,表型貢獻率介于5.92%~22.53%之間,分布在第1、第3、第6、第7、第8、第9和第11染色體上。其中,百粒重QTLSW7與SW7’,SW6.1與SW6’,粒長QTLSL6.1與SL6.1’,粒厚QTLSH11與SH11’在2個環(huán)境下的標記區(qū)間重疊或者重合,SW7、SW6.1、SL6.1、SW6’和SL6.1’的表型貢獻率在10%以上。

普通菜豆;籽粒大小;粒形;數(shù)量性狀位點(QTL)

普通菜豆(Phaseolus vulgaris L.)是世界上種植面積最大的食用豆類,在非洲和美洲的部分地區(qū)作為一種重要的食物來源被廣泛種植[1-3]。在我國,普通菜豆是重要的出口創(chuàng)匯作物,主要分布在黑龍江、內(nèi)蒙古、山西、陜西、四川、貴州、云南等省區(qū)。籽粒是普通菜豆主要的消費對象,在長期飲食習(xí)慣和文化的影響下,不同地區(qū)的人們對籽粒的大小和形狀形成了不同的偏好,籽粒大小和形狀作為重要的品質(zhì)性狀深刻影響著普通菜豆的商業(yè)價值[2]。同時,籽粒大小也是普通菜豆產(chǎn)量要素之一,在長期育種實踐中扮演著不可或缺的角色。因此,解析普通菜豆籽粒大小的遺傳基礎(chǔ)將有效地促進普通菜豆的育種進程。

籽粒性狀在水稻中的研究較為深入,已鑒定出大量相關(guān)基因,其中包括控制粒長和粒重的GS3[4]基因,控制粒寬和粒重的GW2[5]基因,控制粒寬和粒重的qSW5[6]基因等,與粒重相關(guān)的基因GS3、GS5、qSW5、GW5、GW8、GIF1在馴化和育種進程中與產(chǎn)量密切相關(guān)[4-10],這些基因的分離和鑒定極大地促進了水稻的育種進程。普通菜豆籽粒大小QTL定位的研究基礎(chǔ)相對薄弱,而且大部分集中在國外。普通菜豆的籽粒大小是多基因控制的數(shù)量性狀[11],遺傳基礎(chǔ)相對復(fù)雜,籽粒大小又可進一步劃分為粒長、粒寬、粒厚等遺傳組分[12]。1923年,Sax[11]報道了與百粒重連鎖的P基因。隨著分子標記的不斷發(fā)展,遺傳連鎖作圖已成為解析數(shù)量性狀遺傳基礎(chǔ)的重要方法。2006年,Blair等[13]利用哥倫比亞大粒栽培種與野生種雜交獲得的 BC2F3:5群體以及SSR、SCAR、朊蛋白標記構(gòu)建的連鎖圖譜定位到10個與百粒重相關(guān)的QTL,位于第2、第3、第6、第7、第8、第9、第10和第11染色體上;2010年,Pérez-Vega等[2]利用來自2個基因庫間的材料構(gòu)建的RIL群體在第2、第3、第6、第8和第10染色體上定位到5個與粒長相關(guān)的QTL,在第3、第6和第7染色體上定位到3個與粒寬相關(guān)的QTL,在第6、第8染色體上定位到2個與粒厚相關(guān)的QTL,以及在第6、第8染色體上定位到3個與百粒重相關(guān)的QTL; 2014年,Fernando等[14]利用2個來自安第斯基因庫的親本雜交獲得的185個RILs定位到6個與粒長相關(guān)的QTL,位于第1、第2、第6、第7和第10染色體上,3個與粒寬相關(guān)的QTL,位于第2、第7和第9染色體上,2個與粒厚相關(guān)的QTL,位于第2、第9染色體上,以及4個與百粒重相關(guān)的QTL,位于第1、第6和第9染色體上。之前的報道中所用的定位群體的親本多來自同一基因庫,親本間的多態(tài)性較低,在一定程度上影響了QTL定位的準確度與精確度。中美起源中心(主要包括墨西哥、危地馬拉、哥倫比亞等國家)和安第斯起源中心(主要包括秘魯南部、玻利維亞、阿根廷等國家和地區(qū))由于長期地理隔離導(dǎo)致生殖隔離而形成了2個獨立的基因庫,在形態(tài)特征方面,安第斯起源中心的栽培種籽粒大小(百粒重大于40 g)普遍大于中美起源中心栽培種的籽粒大小(百粒重25~40 g,或小于25 g),因此發(fā)掘不同起源中心種質(zhì)的優(yōu)異籽粒性狀基因有利于推進普通菜豆的育種進程。普通菜豆2個基因庫典型材料全基因組測序的完成為解析重要農(nóng)藝性狀的遺傳基礎(chǔ)提供了新的契機[3,15]。

本研究采用分別來自安第斯基因庫和中美基因庫的龍270709和F5910為親本配制雜交組合,獲得F2分離群體作為定位群體,有利于發(fā)掘來自不同基因庫的新基因。根據(jù)普通菜豆全基因組序列設(shè)計SSR標記,有利于QTL信息的比對和綜合分析。多環(huán)境F2定位數(shù)據(jù)利于獲得穩(wěn)定的QTL。本研究為進一步精細定位籽粒性狀相關(guān)基因及分子標記輔助育種奠定了良好基礎(chǔ)。

1 材料與方法

1.1 材料及田間試驗

以龍270709為母本,F5910為父本構(gòu)建F2群體。龍270709為2005年黑龍江省農(nóng)業(yè)科學(xué)院育成品系,來自安第斯基因庫,生長習(xí)性為矮生有限型,生育期100 d左右,籽粒較大,百粒重約70 g;F5910為哥倫比亞品種,來自中美基因庫,生長習(xí)性為矮生無限型,生育期90 d左右,籽粒較小,百粒重約20 g。親本在粒長、粒寬、粒厚等籽粒大小相關(guān)性狀間均存在較大差異。2015年5月中旬,在黑龍江省農(nóng)業(yè)科學(xué)院國家現(xiàn)代農(nóng)業(yè)示范區(qū)種植400株F2單株,每行10株,株距44 cm,行距50 cm,田間管理與當(dāng)?shù)卮筇锕芾硐嗤?346株單株最終用于本試驗;同年11月底,在北京中國農(nóng)業(yè)科學(xué)院作物科學(xué)研究所昌平試驗基地溫室種植的包括188個F2單株的群體用于本試驗,種植密度及田間管理同上,溫室溫度設(shè)置為白天28℃±2℃,夜晚20℃±2℃。成熟后,單株收取干籽粒,用于表型鑒定。

1.2 性狀調(diào)查

收獲后調(diào)查百粒重(seed weight,SW)、粒長(seed length,SL)、粒寬(seed width,SWI)、粒厚(seed height, SH)、長寬比(seed length-to-width ratio,SLW)和長厚比(seed length-to-height ratio)。粒長為平行于種臍,上下2個頂點之間的距離,粒寬為垂直于種臍,上下2個面之間的距離,粒厚為種臍到它對立面的距離[14],長寬比=粒長/粒寬,長厚比=粒長/粒厚。從每個單株取健康籽粒20粒,采用SC-G自動考種分析儀測量粒長、粒厚,采用游標卡尺測量粒寬,利用電子分析天平測量百粒重。

1.3 統(tǒng)計分析

采用SAS9.2軟件和Microsoft Excel 2007進行各性狀間相關(guān)性分析,以及基本統(tǒng)計量的計算。

1.4 QTL定位

在開花前,取親本和 F2單株幼嫩葉片,以CTAB法提取葉片DNA[16]。本實驗室根據(jù)2014年公布的普通菜豆全基因組序列(https://phytozome. jgi.doe.gov/)[3],開發(fā)了2000余對SSR標記。篩選獲得的親本間多態(tài)性標記用于本試驗基因型的鑒定。利用IciMapping 4.1計算標記間的遺傳距離,設(shè)置LOD≥3,利用 Kosambi作圖函數(shù)作圖,利用JoinMap 4根據(jù)遺傳距離繪制連鎖圖譜[17-19]。采用高用明等[20]推薦的 QTL命名方法命名。利用IciMapping 4.1檢測F2群體籽粒大小和形狀QTL,采用完備區(qū)間作圖方法,設(shè)置檢測步長1.0 cM和1000次排列檢驗確定QTL的臨界值(P<0.05)[17-18]。

2 結(jié)果與分析

2.1 表型變異及相關(guān)分析

從表1可以看出,母本龍270709的百粒重、粒長、粒寬、粒厚、長寬比、長厚比均高于父本F5910,其中百粒重的差異最為明顯。母本270709百粒重在黑龍江哈爾濱和北京昌平溫室2個環(huán)境下的均值差別較大,可能是品種特性與環(huán)境因素共同作用的結(jié)果。F2分離群體在2個環(huán)境下不同性狀間均存在較大程度的分離,變異系數(shù)介于5.17%~20.35%之間,其中百粒重的變異系數(shù)最大,2個地點分別為16.62%和20.35%,長厚比的變異系數(shù)最小,分別為5.17%和6.62%。各性狀的表型值為連續(xù)變異,接近正態(tài)分布(圖1)。

相關(guān)性分析結(jié)果(表2)顯示,百粒重分別與粒長、粒寬、粒厚、長寬比及長厚比顯著正相關(guān)。粒長分別與粒寬、粒厚、長寬比和長厚比顯著正相關(guān)。粒寬與粒厚顯著正相關(guān),與長寬比顯著負相關(guān),與長厚比在哈爾濱環(huán)境(環(huán)境1,下同)下顯著正相關(guān),在北京昌平環(huán)境下(環(huán)境2,下同)相關(guān)性不顯著。粒厚與長寬比顯著正相關(guān),與長厚比相關(guān)性不顯著。長寬比與長厚比顯著正相關(guān)。

2.2 遺傳圖譜構(gòu)建

利用哈爾濱種植的346株F2和北京昌平種植的188株F2分別構(gòu)建了連鎖圖譜(圖2和圖3)。346株F2構(gòu)建的遺傳圖譜包括100個SSR標記,分布于11條染色體上,每條染色體上的標記個數(shù)5~17個,連鎖圖譜總長439.88 cM,標記間平均遺傳距離4.40 cM,有2對標記間的遺傳距離大于20.00 cM,分別位于第10染色體P10S167~P10S186(31.91 cM)和第4染色體P4S2~P4S22(22.28 cM),其余標記遺傳距離均小于20 cM。188株F2構(gòu)建的遺傳圖譜包括88個SSR標記,與前者相比,第5、第6、第7、第9和第11染色體分別缺少了4、4、2、1、1個標記,為F2基因型缺失較嚴重的標記,其余標記均與前者相同,每條染色體上的標記個數(shù)3~13個,連鎖圖譜總長421.94 cM,標記間平均遺傳距離4.80 cM,標記間遺傳距離大于20 cM的2對標記同樣位于第10染色體P10S167~P10S186(31.18 cM)和第4染色體P4S2~ P4S22(26.57 cM),其余標記的遺傳距離均小于20 cM。

圖1 群體籽粒大小與形狀相關(guān)性狀頻率分布圖Fig.1 Frequency distribution of seed size and shape traits of F2population in Harbin and Beijing

表2 籽粒性狀間的相關(guān)性分析Table 2 Correlation coefficients for seed traits in common bean

2.3 QTL定位

利用IciMapping 4.1進行1000次排列獲得QTL顯著的LOD閾值均為3.0。哈爾濱環(huán)境下定位到的QTL在11條染色體上的分布如圖2所示,北京昌平環(huán)境下定位到的QTL在11條染色體上的分布如圖3所示。2個環(huán)境下QTL的定位結(jié)果見表3,其中哈爾濱環(huán)境下共定位到38個與百粒重、粒長、粒寬、粒厚、長寬比和長厚比相關(guān)的QTL,北京昌平環(huán)境下共定位到21個與上述6個性狀相關(guān)的QTL。

2.3.1 百粒重 環(huán)境1中共檢測到7個與百粒重相關(guān)的QTL(SW7、SW6.1、SW6.2、SW9、SW4、SW2和SW10),分布在第2、第4、第6、第7、第9和第10染色體上,表型貢獻率介于3.88%~10.97%之間,加性效應(yīng)均為正,表明增效等位基因均來自大粒親本龍270709。第7染色體上的SW7,表型貢獻率最大,為10.97%。第6染色體上的SW6.1和SW6.2,表型貢獻率分別為10.79%和7.76%,遺傳效應(yīng)較大。SW9、SW4、SW2、SW10分別位于第9、第4、第2和第10染色體上,表型貢獻率分別為 5.56%、5.27%、4.42%和3.88%,遺傳效應(yīng)較小。

環(huán)境2中共檢測到4個與百粒重相關(guān)的QTL (SW6’、SW7’、SW11’、SW9’),分布在第6、第7、第11和第9染色體上,表型貢獻率分別為15.88%、8.22%、8.16%和7.85%,加性效應(yīng)均為正。

比較2個環(huán)境下的定位結(jié)果比較,第6、第7和第9染色體上均檢測到百粒重QTL,SW7(P7S85~ P7S129)與SW7’(P7S85~P7S140)標記區(qū)間存在重合,與標記P7S85的遺傳距離均為0 cM。SW6.1與SW6’均位于標記P6S108~P6S266之間,與標記P6S108的遺傳距離分別為 4.50 cM和 0.32 cM。SW9 (P9S147~P9S155)與 SW9’(P9S136~P9S147)標記區(qū)間相鄰,與標記P9S147的遺傳距離分別為2.17 cM和0.47 cM。SW4、SW2、SW10環(huán)境貢獻率均小于SW9,未在環(huán)境2中檢測到,SW6.2貢獻率比SW7和SW6.1小,比SW9大,在環(huán)境2中未被檢測到。

2.3.2 粒長 環(huán)境1中共檢測到8個與粒長相關(guān)的QTL(SL6.1、SL7、SL6.2、SL9、SL5、SL10、SL2和SL4),分布在第2、第4、第5、第6、第7、第9和第10染色體上,除SL5外,其余位點加性效應(yīng)均為正。SL6.1和SL7遺傳效應(yīng)較大,表型貢獻率分別為17.37%和12.53%。其余6個QTL(SL6.2、SL9、SL5、SL10、SL2、SL4)表型貢獻率分別為8.57%、4.46%、4.25%、3.96%、2.94%和2.64%。

在環(huán)境2定位到5個粒長QTL(SL6.1’、SL7’、SL6.2’、SL1’和SL9’),分布于第1、第6、第7和第9染色體上,表型貢獻率介于6.01%~22.53%之間, SL1’加性效應(yīng)為負值,其余均為正值。SL6.1’與SL7’遺傳效應(yīng)較大,表型貢獻率分別為 22.53%和10.14%。SL6.2’、SL1’、SL9’表型貢獻率分別為9.10%、7.14%和6.01%。

2個環(huán)境下定位結(jié)果比較分析,均在第6、第7、第9染色體上定位到粒長QTL。SL6.1與SL6.1’均位于標記P6S108~P6S266之間,與標記P6S108之間的遺傳距離分別為 4.50 cM 和 0.32 cM。SL7 (P7S85~P7S129)和 SL7’(P7S140~P7S153)與標記P7S85的遺傳距離分別為 3 cM 和 4 cM。SL9 (P9S86~P9S129)和 SL9’(P9S136~P9S147)與標記P9S147間的遺傳距離分別為4.83 cM和0.47 cM。SL6.2與SL6.2’均位于第6染色體,但是距離較遠,與標記P6S185之間的遺傳距離分別是7.00 cM和28.82 cM。另外,在環(huán)境2定位到位于第1染色體的SL1’并沒有在環(huán)境1中檢測到,環(huán)境1中檢測到的的SL5、SL10、SL2、SL4表型貢獻率相對較小,沒有在環(huán)境2中檢測到。

2.3.3 粒寬 在環(huán)境1定位到3個粒寬QTL (SWI2、SWI6、SWI5),分布于第2、第5和第6條染色體上,表型貢獻率介于4.36%~7.68%之間,除SWI5加性效應(yīng)為負值外,其余位點加性效應(yīng)均為正值。

圖2 哈爾濱環(huán)境下(環(huán)境1)定位到的QTL在連鎖群上的分布Fig.2 Distribution of QTL for seed size and shape in linkage groups in Harbin,Heilongjiang(environment 1)

在環(huán)境2定位到3個粒寬QTL(SWI11’、SWI8’和SWI9’),分布于第8、第9、第11染色體上,表型貢獻率介于8.22%~10.72%之間,加性效應(yīng)均為正值。2個群體中定位到的粒寬QTL表型貢獻率均沒有超過10%,同時也沒有位于相同染色體上的位點,表明粒寬受環(huán)境影響較大。

圖3 北京昌平環(huán)境下(環(huán)境2)定位到的QTL在連鎖群上的分布Fig.3 Distribution of QTL for seed size and shape in linkage groups in Beijing(environment 2)

表3 哈爾濱環(huán)境(環(huán)境1)和北京昌平(環(huán)境2)的QTL定位結(jié)果Table 3 QTL detected in Harbin,Heilongjiang(environment 1)and Changping,Beijing(environment 2)

(續(xù)表3)

2.3.4 粒厚 在環(huán)境1定位到9個粒厚QTL (SH8、SH6.1、SH6.2、SH7、SH2、SH3.1、SH11、SH4和SH3.2),分布于第2、第3、第4、第6、第7、第8和第11染色體上,表型貢獻率介于3.29%~ 10.91%之間,在第3和第6染色體上分別檢測到2個QTL,在其余染色體上均檢測到1個QTL。SH8和SH3加性效應(yīng)為負值,其余位點遺傳效應(yīng)均為正值。SH8表型貢獻率最大,為10.91%,SH6.1表型貢獻率次之,為9.62%。其余位點(SH6.2、SH7、SH2、SH3.1、SH11、SH4、SH3.2)表型貢獻率分別為7.11%、6.84%、4.80%、4.26%、4.06%、3.29%和2.39%。在環(huán)境2定位到4個粒厚QTL(SH7’、SH1’、SH6’、SH11’),分布于第1、第6、第7、第11染色體上,SH1’加性效應(yīng)為負值,其余位點加性效應(yīng)均為正值。表型貢獻率在10%以上的位點有2個,分別是SH7’(17.78%)和SH1’(10.23%)。SH6’和SH11’表型貢獻率分別為8.92%和7.20%。

(續(xù)表3)

環(huán)境1中檢測到的SH8表型貢獻率最大,但在環(huán)境2中并未被檢測到。SH6.1(P6S108~P6S266)和SH6’(P6S275~P6S283)與標記P6S108間的遺傳距離分別為4.50 cM和23.68 cM。SH6.2(P6S316~P6S289)和SH6’(P6S275~P6S283)與標記P6S289間的遺傳距離分別為0.08 cM和21.22 cM?？梢奡H6’與SH6.1和SH6.2距離都比較遠。SH7(P7S85~P7S129)和SH7’(P7S140~P7S153)與標記P7S85間的遺傳距離分別為 0 cM和 4 cM。SH11與 SH11’均位于標記P11S21-P11S18之間,與標記P11S18間的遺傳距離分別為0.32 cM和0.90 cM。

2.3.5 長寬比 環(huán)境1中定位到4個長寬比QTL (SLW6、SLW7、SLW8、SLW10),分布于第6、第7、第8、第9染色體上,SLW8加性效應(yīng)值為負,其余均為正,表型貢獻率最大的QTL是SLW6,為9.27%,最小的是SLW10,為3.48%。環(huán)境2中定位到3個長寬比QTL(SLW6’、SLW8’、SLW7’),分布于第6、第7、第8染色體,SLW8’加性效應(yīng)值為負,其余為正,3個長寬比QTL(SLW6’、SLW8’、SLW7’)表型貢獻率分別為20.29%、9.88%和7.31%。

2個環(huán)境下均在第6、第7、第8染色體上定位到長寬比QTL。在第6染色體上定位到的SLW6 (P6S283~P6S316)和SLW6’(P6S266~P6S309)標記區(qū)間不同且不相鄰,與標記P6S316間的遺傳距離為8.08 cM和36.22 cM,相距較遠。在第7染色體上定位到的SLW7(P7S171~P7S182)與SLW7’(P7S153~ P7S171)標記區(qū)間相鄰,與標記P7S182間的遺傳距離分別為5.41 cM和11.80 cM。在第8染色體上定位到的SLW8(P8S141~P8S152)和SLW8’(P8S106~ P8S141)標記區(qū)間相鄰,與標記P8S152間的連鎖距離分別是4.11 cM和4.21 cM。

2.3.6 長厚比 環(huán)境1中定位到7個與長厚比相關(guān)的QTL(SLH3、SLH7、SLH6、SLH5、SLH8、SLH9和SLH10),分布在第3、第5、第6、第7、第8、第9和第10染色體上,加性效應(yīng)值均為正,表型貢獻率介于2.83%~10.99%之間。環(huán)境2中定位到2個與長厚比相關(guān)的QTL(SLH6’、SLH3’),位于第6和第3染色體上,表型貢獻率分別為31.59%和5.92%,加性效應(yīng)值均為正。

2個環(huán)境下均在第3和第6染色體上定位到長厚比 QTL,位于第 3染色體上的 SLH3(P3S79~ P3S137)與 SLH3’(P3S169~P3S171)標記區(qū)間不同,與標記P3S171間的遺傳距離分別為20.98 cM和0.76 cM,相距較遠。位于第 6染色體上的SLH6(P6S309~P6S275)與SLH6’(P6S108~P6S266)標記區(qū)間不同,與P6S108間的遺傳距離分別是7.90 cM和0.32 cM。

3 討論

3.1 QTL定位結(jié)果分析

普通菜豆的籽粒特性是重要的農(nóng)藝性狀,直接影響著普通菜豆的產(chǎn)量和商品性。一直被研究者關(guān)注,截至目前在菜豆中已經(jīng)定位到大量的與籽粒特性相關(guān)的QTL[2,13-14,21]。本研究在哈爾濱環(huán)境下共定位到38個與百粒重、粒長、粒寬、粒厚、長寬比和長厚比相關(guān)的QTL,在昌平環(huán)境下共定位到21個與上述性狀相關(guān)的QTL。其中,位于第2染色體的百粒重QTL SW2,標記區(qū)間P2S11~P2S13,在染色體上的物理位置為2.48~3.11 Mb,與Blair等[13]定位到的百粒重QTL sw2.2(2.819 Mb)位置重疊;位于第6染色體標記 P6S289~P6S185之間的百粒重 QTL SW6.2,物理位置在29.57~30.77 Mb之間,與Blair等[13]檢測到的百粒重QTL sw6.1(27.44 Mb)位置重疊;位于第1染色體上的粒長QTLSL1’,標記區(qū)間P1S77~P1S82,物理位置在39.77~42.35 Mb之間,與Fernando等[14]檢測到的SL1.2PP(32.41 Mb)位置相近。定位結(jié)果的一致性一方面說明定位結(jié)果的可靠性,另一方面也說明這些位點的重要性。我們還注意到在第4染色體上未見前人報道與籽粒大小相關(guān)的位點,而本研究在這一染色體上檢測到與百粒重、粒長、粒厚相關(guān)的QTL,位于標記P4S2~P4S22之間。

同一群體多環(huán)境表型鑒定及QTL分析,可以檢測到相對穩(wěn)定可靠的位點。特別是貢獻率(加性效應(yīng))較大的QTL更容易在不同環(huán)境條件下被連續(xù)檢測到[22]。本研究中2個環(huán)境下共檢測到4對標記區(qū)間重疊的QTL:百粒重QTL SW7和SW7’分別位于標記區(qū)間P7S85~P7S129和 P7S85~P7S140,百 粒 重 QTL SW6.1(SW6’)位于標記區(qū)間P6S108~P6S266,粒長QTL SL6.1(SL6.1’)位于標記區(qū)間 P6S108~P6S266,粒厚 QTL SH11(SH11’)位于標記區(qū)間 P11S18~ P11S21。其中3對QTL至少在一個環(huán)境下表型貢獻率在 10%以上,SW7(SW7’)表型貢獻率 10.91% (8.22%),SW6.1(SW6’)表型貢獻率10.79%(15.88%), SL6.1(SL6.1’)表型貢獻率17.37%(22.53%)。這種具有較高貢獻率且在不同的環(huán)境下可連續(xù)檢測到的QTL可以更好地用于分子標記輔助選擇育種[23]。4對QTL加性效應(yīng)均為正,表明增效等位基因均來自母本龍270709,這將有利于發(fā)掘和利用來自安第斯基因庫的普通菜豆優(yōu)異種質(zhì)的籽粒大小和形狀基因。除此之外,有3對QTL標記區(qū)間雖然不重疊但是相鄰,分別位于第7、第8、第9染色體上,位于第7染色體上的長寬比QTL SLW7和SLW7’分別位于標記區(qū)間P7S171~P7S182和P7S153~P7S171,位于第8染色體上的長寬比QTL SLW8和SLW8’分別位于標記區(qū)間P8S141~P8S152和P8S106~P8S141,位于第9染色體上的百粒重QTL SW9和SW9’分別位于標記區(qū)間P9S147~P9S155和P9S136~P9S147。SLW8(SLW8’)、SW9(SW9’)與SLW7(SLW7’)表型貢獻率均介于4%~10%之間,分別為4.09%(9.88%)、5.56%(7.85%)和6.71%(7.31%)?；蛐团c環(huán)境互作可能是導(dǎo)致這類情況出現(xiàn)的主要原因,尤其是對于表型貢獻率相對較小的位點,環(huán)境的影響更為明顯。精細定位或者標記選擇時可考慮將其標記區(qū)間擴大來利用這些位點。

另外,由于環(huán)境2所用群體比環(huán)境1所用群體相對較小,環(huán)境1中檢測到的表型貢獻率較小的位點在環(huán)境2中并未被檢測到,如SW2、SW4、SW10、SL2、SL4、SL5、SL10、SH2、SH3和SH4,包括與百粒重、粒長和粒厚相關(guān)的QTL,集中分布在第2、第3、第4、第5染色體上。群體大小影響檢測的效率,群體擴大,效應(yīng)值較小的位點也能被檢測到[24]。

3.2 染色體區(qū)段多效性分析

本研究中出現(xiàn)不同性狀定位在相同或相近染色體區(qū)段的現(xiàn)象。這種現(xiàn)象較為普遍,在小麥、玉米、大豆等作物中均有類似報道[25-29]。Tuberosa等[30]認為“一因多效”可能是這種現(xiàn)象出現(xiàn)的主要原因,例如水稻中的GS3基因既控制粒長,又與百粒重有關(guān), GW2同時與粒寬和百粒重相關(guān)[4-5]。另外,QTL緊密連鎖也可能是導(dǎo)致這種結(jié)果的一個原因。本研究中,第 6染色體 P6S108~P6S266區(qū)段、第 7染色體P7S85~P7S129區(qū)段內(nèi)均同時檢測到與百粒重、粒長、粒厚和長厚比均相關(guān)的QTL;第4染色體P4S2~P4S22區(qū)段檢測到與百粒重、粒長和粒厚相關(guān)的QTL;第2染色體P2S11~P2S13區(qū)段內(nèi)檢測到與百粒重、粒寬相關(guān)的QTL,P2S56~P2S68區(qū)段內(nèi)檢測到與粒長、粒厚相關(guān)的QTL;第10染色體P10S3~ P10S86區(qū)段內(nèi)檢測到與百粒重、粒長、長厚比相關(guān)的QTL;第6染色體P6S275~P6S283與第7染色體P7S140~P7S153區(qū)段內(nèi)均同時檢測到與粒長、粒厚相關(guān)的QTL;第9染色體P9S136~P9S147區(qū)段內(nèi)同時檢測到與百粒重、粒長、粒寬相關(guān)的QTL。百粒重與粒長、粒厚、粒寬、長厚比相關(guān)的QTL在同一染色體區(qū)段被同時檢測到,這為解釋百粒重與粒長、粒厚、粒寬及長厚比之間顯著的相關(guān)性提供了一定的遺傳基礎(chǔ)。

4 結(jié)論

黑龍江哈爾濱環(huán)境中定位到38個與百粒重、粒長、粒寬、粒厚、長寬比和長厚比相關(guān)的QTL,北京昌平環(huán)境下定位到上述性狀的QTL共21個,其中百粒重QTLSW7(SW7’)和SW6.1(SW6’),粒長QTL SL6.1(SL6.1’),粒厚QTLSH11(SH11’)在2個環(huán)境下被重復(fù)檢測到,4個位點增效等位基因均來自安第斯基因庫的大粒親本龍270709,可作為精細定位的重要候選區(qū)段以及標記輔助選育的位點。

[1]張赤紅,曹永生,宗緒曉,王志剛,王述民.普通菜豆種質(zhì)資源形態(tài)多樣性鑒定與分類研究.中國農(nóng)業(yè)科學(xué),2005,38: 27~32

Zhang C H,Cao Y S,Zong X X,Wang Z G,Wang S M. Morphological diversity and classification of common bean (Phaseolus vulgaris L.)germplasm resource in China.Sci Agric Sin,2005,38:27–32(in Chinese with English abstract)

[2]Pérez-Vega E,Pa?eda A,Rodríguez-SuárezC,Campa A, Giraldez R, Ferreira J J. Mapping of QTLs for morpho-agronomic and seed quality traits in a RIL population of common bean(Phaseolus vulgaris L.).Theor Appl Genet,2010, 120:1367–1380

[3]Schmutz J,McClean P E,Mamidi S,Wu G A,Cannon S B, Grimwood J,Jenkins J,Shu S Q,Song Q J,Chavarro C, Torres-Torres M,Geffroy V,oghaddam S M,Gao D Y,Abernathy B,Barry K,Blair M,Brick M A,Chovatia M,Gepts P,Goodstein D M,Gonzales M,Hellsten U,Hyte D L,Jia G F,Kelly J D, Kudrna D,Lee R,Richard M M S,Miklas P N,Osorno J M, Rodrigues J,Thareau V,Urrea C A,Wang M,Yu Y,Zhang M, Wing R A,Cregan P B,Rokhsar D S,Jackson S A.A reference genome for common bean and genome-wide analysis of dual domestications.Nat Genet,2014,46:707–713

[4]Fan C C,Xing Y Z,Mao H L,Lu T T,Han B,Xu C G,Li X H, Zhang Q F.GS3,a major QTL for grain length and weight and minor QTL for grain width and thickness in rice,encodes a putative transmembrane protein.Theor Appl Genet,2006,112: 1164–1171

[5]Song X J,Huang W,Shi M,Zhu M Z,Lin H X.A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase.Nat Genet,2007,39:623–630

[6]Shomura A,Izawa T,Ebana K,Ebitani T,Kanegae H,Konishi S, Yano M.Deletion in a gene associated with grain size increased yields during rice domestication.Nat Genet,2008,40:1023–1028

[7]Weng J F,Gu S H,Wan X Y,Gao H,Guo T,Su N,Lei C L, Zhang X,Cheng Z J,Guo X P,Wang J L,Jiang L,Zhai H Q,Wan J M.Isolation and initial characterization of GW5,a major QTL associated with rice grain width and weight.Cell Res,2008,18: 1199–1209

[8]Wang E T,Wang J J,Zhu X D,Hao W,Wang L Y,Li Q,Zhang L X,He W,Lu B R,Lin H X,Ma H,Zhang G Q,He Z H.Control of rice grain-filling and yield by a gene with a potential signature of domestication.Nat Genet,2008,40:1370–1374

[9]Li Y B,Fan C C,Xing Y Z,JiangY H,Luo L J,Sun L,Shao D,Xu C J.Natural variation in GS5 plays an important role in regulating grainsize and yield in rice.Nat Genet,2011,43:1266–1269

[10]Wang S K,Wu K,Yuan Q B,Liu X Y,Liu Z B,Lin X Y,Zeng R Z,Zhu H T.Control of grain size,shape and quality by OsSPL16 in rice.2012,Nat Genet,44:950–954

[11]Sax K.The association of size differences with seed-coat pattern and pigmentation in PHASEOLUS VULGARIS.Genetics,1923,8: 552–560

[12]Vallejos C E,Chase C D.Linkage between isozyme markers and a locus affecting seed size in Phaseolus vulgaris L.Theor Appl Genet,1991,81:413–419

[13]Blair M W,Iriarte G,Beebe S.QTL analysis of yield traits in an advanced backcross population derived from a cultivated Andean ×wild common bean(Phaseolus vulgaris L.)cross.Theor Appl Genet,2006,112:1149–1163

[14]Yuste-Lisbona F J,González A N,Capel C,Garcia-Alcazar M, Capel J,Ron A M D,Lozano R,Santalla M.Genetic analysis of single-locus and epistatic QTLs for seed traits in an adapted× nu?a RIL population of common bean(Phaseolus vulgaris L.). Theor Appl Genet,2014,127:897–912

[15]Vlasova A,Capella-Gutiérrez S,Rendón-Anaya M,Hernández-O?ate M E.Minoche A,Erb I,Camara F,Prieto-Barja P,Corvelo A,Sanseverino W,Westergaard G,Dohm J C,Pappas Jr G J, Saburido-Alvarez S,Kedra D,Gonzalez I,Cozzuto L, Gómez-Garrido J,Aguilar-Morón M A,Andreu N,Aguilar O, Garcia-Mas J,Zehnsdorf M,Vázquez M P,Delgado-Salinas A, Delaye L,Lowy E,Mentaberry A,Vianello-Brondani R P,García J,Alioto T,Sánchez F,Himmelbauer H,Santalla M,Notredame C, Gabaldón T,Herrera-Estrella A,Guigó R.Genome and transcriptome analysis of the Mesoamerican common bean and the role of gene duplications in establishing tissue and temporal specialization of genes.Genome Biol,2016,17:1–18

[16]Murray M G,Thompson W F.Rapid isolation of high molecular weight plant DNA.Nucl Acids Res,1980,8:4321–4325

[17]Li H H,Ye G Y,Wang J K.A modified algorithm for the improvement of composite interval mapping.Genetics,2007,175: 361–374

[18]Meng L,Li H H,Zhang L Y,Wang J K.QTL IciMapping: integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental population.Crop J, 2015,3:269–283

[19]Wan X Y,Wan J M,Jiang L,Wang J K,Zhai H Q,Weng J F, Wang H L,Lei C L,Wang J K,Zhang X.QTL analysis for rice grain length and fine mapping of an identified QTL with stable and major effects.Theor Appl Genet,2006,112:1258–1270

[20]高用明,朱軍,宋佑勝,何慈信,石春海,邢永忠.水稻永久F2群體抽穗期QTL的上位性及其與環(huán)境互作效應(yīng)的分析.作物學(xué)報,2004,30:849–854

Gao Y M,Zhu J,Song Y S,He C X,Shi C H,Xing Y Z.Use of permanent F2population to analyze epistasis and their interaction effects with environments for QTLs controlling heading date in rice.Acta Agron Sin,2004,30:849–854(in Chinese with English abstract)

[21]Hoyos-Villegas V,Song Q,Wright E M,Beebe S E,Kelly D J. Joint linkage QTL mapping for yield and agronomic traits in a composite map of three common bean RIL populations.Crop Sci, 2016,56:1–18

[22]Tanksley S D.Mapping ploygenes.Annu Rev Genet,1993,27: 205–233

[23]陳秀華,于麗娟,羅黎明,陳洪梅,劉麗.玉米分子標記輔助育種及標記開發(fā)研究進展.中國農(nóng)業(yè)科技導(dǎo)報,2016,18(1): 26–31

Chen X H,Yu L J,Luo L M,Chen H M,Liu L.Research progress on maize molecular marker-assisted breeding and marker development.J Agric Sci Technol,2016,18(1):26–31(in Chinese with English abstract)

[24]趙洪波,李明麗,魯紹雄,連林生,李國治.群體規(guī)模和性狀遺傳力對F2設(shè)計下QTL定位效果的影響.云南農(nóng)業(yè)大學(xué)學(xué)報, 2007,22(2):159–163

Zhao H B,Li M L,Lu S X,Lian L S,Li G Z.Study on the effects of population size and trait heritability on the accuracy of QTL mapping under F2design.J Yunnan Agric Univ,2007,22(2): 159–163(in Chinese with English abstract)

[25]Groos C,Robert N,Bervas E,Charmet G.Genetic analysis of grain protein-content,grain yield and thousand-kernel weight in bread wheat.Theor Appl Genet,2003,106:1032–1040

[26]張坤普,徐憲斌,田紀春.小麥籽粒產(chǎn)量及穗部相關(guān)性狀的QTL定位.作物學(xué)報,2009,35:270–278

Zhang K P,Xu X B,Tian J C.QTL mapping for grain yield and spike related traits in common wheat.Acta Agron Sin,2009,35: 270–278(in Chinese with English abstract)

[27]李文福,劉賓,彭濤,袁倩倩,韓淑曉,田紀春.利用DH和IF2兩個群體進行小麥粒重、粒型和硬度的QTL分析.中國農(nóng)業(yè)科學(xué),2012,45:3453–3462

Li W F,Liu B,Peng T,Yuan Q Q,Han S X,Tian J C.Detection of QTL for kernel weight,grain size,and grain hardness in wheat using DH and immortalized F2population.Sci Agric Sin,2012, 45:3453–3462(in Chinese with English abstract)

[28]覃鴻妮,晏萌,王召輝,郭瑩,王輝,孫海燕,劉志齋,蔡一林.玉米籽粒中花色苷和黑色素含量的QTL分析.作物學(xué)報, 2012,38:275–284

Qin H N,Yan M,Wang Z H,Guo Y,Wang H,Sun H Y,Liu Z Z, Cai Y L.QTL mapping for anthocyanin and melanin contents in maize kernel.Acta Agron Sin,2012,38:275–284(in Chinese with English abstract)

[29]鄭德波,楊小紅,李建生,嚴建兵,張士龍,賀正華,黃益勤.基于SNP標記的玉米株高及穗位高QTL定位.作物學(xué)報, 2013:549–556

Zheng D B,Yang X H,Li J S,Yan J B,Zhang S L,He Z H, Huang Y Q.QTL identification for plant height and ear height based on SNP mapping in maize(Zea mays L.).Acta Agron Sin, 2013,39:549–556(in Chinese with English abstract)

[30]Tuberosa R,Salvi S,Sanguineti M C,Landi P,Maccaferri M, Conti S.Mapping QTLs regulating morpho-physiological traits and yield:casestudies,shortcomings and perspectivesin drought-stressed maize.Ann Bot,2002,89:941–963

QTL Mapping for Seed Size and Shape in Common Bean

GENG Qing-He,WANG Lan-Fen,WU Jing,and WANG Shu-Min*

Institute of Crop Science,Chinese Academy of Agricultural Sciences,Beijing 100081,China

Common bean is one of the most important food legumes worldwide.Seed size and shape have a great effect on yield and seed quality.This research used an F2population derived from a cross of Andean cultivar Long 270709′Mesoamerican cultivar F5910 planted in Harbin,Heilongjiang and Beijing respectively to analyze the correlation among seed weight,seed length, seed width,seed height,seed length-to-width ratio and seed length-to-height ratio,and to detect QTL related to these traits.The correlation analysis showed that seed weight had a significant correlation with other traits related to seed size and shape.Software IciMapping 4.1 based on inclusive composite interval mapping was used to identify the additive QTL for seed size and shape.A total of 38 QTL were detected on 10 chromosomes except the first chromosome in Harbin environment with the PVE between 2.39%and 17.37%.A total of 21 QTL were detected on seven chromosomes(Chr.1,Chr.3,Chr.6,Chr.7,Chr.8,Chr.9,Chr.11)in Beijing environment with the PVE between 5.92%and 22.53%.Among them,four pairs of QTL were detected on the same marker interval,including seed weight QTL SW7 and SW7’,seed weight QTL SW6.1 and SW6’,seed length QTL SL6.1 and SL6.1’,and seed height QTL SH11 and SH11’.The PVE of SW7,SW6.1,SL6.1,SW6’,and SL6.1’was larger than 10%.

Common bean;Seed size;Seed shape;Quantitative trait locus(QTL)

(

):2016-12-21;Accepted(接受日期):2017-04-20;Published online(網(wǎng)絡(luò)出版日期):2017-04-27.

10.3724/SP.J.1006.2017.01149

本研究由國家自然科學(xué)基金項目(31471559),國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(CARS-09)和國家科技支撐計劃項目(2013BAD01B05-2-4)資助。

This study was supported by the National Natural Science Foundation of China(31471559),the China Agriculture Research System (CARS-09),and the National Key Technology R&D Program of China(2013BAD01B05-2-4).

*通訊作者(Corresponding author):王述民,E-mail:wangshumin@caas.cn

聯(lián)系方式:E-mail:17710428642@163.com

URL:http://kns.cnki.net/kcms/detail/11.1809.S.20170427.0948.016.html