• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2017-08-22 05:59:47耿慶河王蘭芬王述民
    作物學報 2017年8期
    關(guān)鍵詞:粒長貢獻率表型

    耿慶河 王蘭芬 武 晶 王述民

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

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

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

    中國農(nóng)業(yè)科學院作物科學研究所,北京 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ū)。籽粒是普通菜豆主要的消費對象,在長期飲食習慣和文化的影響下,不同地區(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ǔ)相對復雜,籽粒大小又可進一步劃分為粒長、粒寬、粒厚等遺傳組分[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ū))由于長期地理隔離導致生殖隔離而形成了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è)科學院育成品系,來自安第斯基因庫,生長習性為矮生有限型,生育期100 d左右,籽粒較大,百粒重約70 g;F5910為哥倫比亞品種,來自中美基因庫,生長習性為矮生無限型,生育期90 d左右,籽粒較小,百粒重約20 g。親本在粒長、粒寬、粒厚等籽粒大小相關(guān)性狀間均存在較大差異。2015年5月中旬,在黑龍江省農(nóng)業(yè)科學院國家現(xiàn)代農(nóng)業(yè)示范區(qū)種植400株F2單株,每行10株,株距44 cm,行距50 cm,田間管理與當?shù)卮筇锕芾硐嗤?346株單株最終用于本試驗;同年11月底,在北京中國農(nóng)業(yè)科學院作物科學研究所昌平試驗基地溫室種植的包括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)境互作可能是導致這類情況出現(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緊密連鎖也可能是導致這種結(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)境下被重復檢測到,4個位點增效等位基因均來自安第斯基因庫的大粒親本龍270709,可作為精細定位的重要候選區(qū)段以及標記輔助選育的位點。

    [1]張赤紅,曹永生,宗緒曉,王志剛,王述民.普通菜豆種質(zhì)資源形態(tài)多樣性鑒定與分類研究.中國農(nóng)業(yè)科學,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)的分析.作物學報,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è)科技導報,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è)大學學報, 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定位.作物學報,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è)科學,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分析.作物學報, 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定位.作物學報, 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

    本研究由國家自然科學基金項目(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

    猜你喜歡
    粒長貢獻率表型
    水稻粒長遺傳及其功能基因研究進展
    一種通用的裝備體系貢獻率評估框架
    關(guān)于裝備體系貢獻率研究的幾點思考
    建蘭、寒蘭花表型分析
    秈稻粒長與稻米品質(zhì)的相關(guān)性及其育種應(yīng)用
    日本晴/R1126水稻重組自交系群體粒形性狀QTL定位
    В первой половине 2016 года вклад потребления в рост китайской экономики достиг 73,4 процента
    中亞信息(2016年10期)2016-02-13 02:32:45
    GABABR2基因遺傳變異與肥胖及代謝相關(guān)表型的關(guān)系
    慢性乙型肝炎患者HBV基因表型與血清學測定的臨床意義
    小麥粒型相關(guān)性狀的QTL定位分析
    久久精品国产综合久久久| www日本在线高清视频| 老司机午夜十八禁免费视频| 日韩制服骚丝袜av| 中文字幕色久视频| 日韩制服骚丝袜av| 亚洲精品乱久久久久久| 国产97色在线日韩免费| 日本色播在线视频| 久久久久久久久久久久大奶| 国产成人精品无人区| 精品免费久久久久久久清纯 | av网站免费在线观看视频| 亚洲av欧美aⅴ国产| 国产伦人伦偷精品视频| 国产免费福利视频在线观看| 成年美女黄网站色视频大全免费| 欧美xxⅹ黑人| 亚洲av欧美aⅴ国产| 只有这里有精品99| 操出白浆在线播放| 久久精品久久精品一区二区三区| 免费看av在线观看网站| 尾随美女入室| 色精品久久人妻99蜜桃| 国产一区二区在线观看av| 狂野欧美激情性xxxx| 国产精品 欧美亚洲| 激情五月婷婷亚洲| 9色porny在线观看| 黄色毛片三级朝国网站| 天堂8中文在线网| 人人澡人人妻人| av在线播放精品| 午夜福利视频精品| 女性被躁到高潮视频| 一级a爱视频在线免费观看| 另类精品久久| 亚洲精品国产av蜜桃| 午夜精品国产一区二区电影| 老熟女久久久| 美国免费a级毛片| 丁香六月欧美| 人人妻人人添人人爽欧美一区卜| 男女高潮啪啪啪动态图| 9色porny在线观看| 久久 成人 亚洲| av国产久精品久网站免费入址| 午夜福利影视在线免费观看| 国语对白做爰xxxⅹ性视频网站| 2021少妇久久久久久久久久久| 一级毛片我不卡| www日本在线高清视频| 亚洲精品第二区| 少妇粗大呻吟视频| 好男人电影高清在线观看| 国产精品麻豆人妻色哟哟久久| 考比视频在线观看| 亚洲av日韩精品久久久久久密 | 熟女少妇亚洲综合色aaa.| av在线播放精品| 国产精品久久久久久精品古装| 91麻豆精品激情在线观看国产 | 大片电影免费在线观看免费| 久久久久网色| av视频免费观看在线观看| 亚洲av欧美aⅴ国产| 午夜激情久久久久久久| 97精品久久久久久久久久精品| h视频一区二区三区| 日本a在线网址| av片东京热男人的天堂| 9热在线视频观看99| 国产成人系列免费观看| 亚洲一区二区三区欧美精品| 成人影院久久| 在线亚洲精品国产二区图片欧美| 少妇裸体淫交视频免费看高清 | av天堂久久9| 秋霞在线观看毛片| 女性生殖器流出的白浆| 91麻豆av在线| 丰满饥渴人妻一区二区三| a级片在线免费高清观看视频| 汤姆久久久久久久影院中文字幕| av在线播放精品| 18禁观看日本| 欧美日韩亚洲综合一区二区三区_| 免费不卡黄色视频| 我要看黄色一级片免费的| 久久国产精品男人的天堂亚洲| 最近最新中文字幕大全免费视频 | 国产成人一区二区在线| 少妇人妻 视频| 亚洲国产最新在线播放| 日本91视频免费播放| 免费观看人在逋| 可以免费在线观看a视频的电影网站| av不卡在线播放| 国产免费又黄又爽又色| 视频在线观看一区二区三区| 国产av国产精品国产| 好男人电影高清在线观看| 国产高清不卡午夜福利| 可以免费在线观看a视频的电影网站| 成年人午夜在线观看视频| 国产精品香港三级国产av潘金莲 | 狠狠精品人妻久久久久久综合| av有码第一页| 亚洲综合色网址| 下体分泌物呈黄色| 中文字幕人妻丝袜一区二区| 国产极品粉嫩免费观看在线| 精品久久久精品久久久| 精品一区二区三卡| 国产精品成人在线| 91精品伊人久久大香线蕉| 晚上一个人看的免费电影| 亚洲精品美女久久av网站| 久久精品久久久久久噜噜老黄| 热re99久久国产66热| 三上悠亚av全集在线观看| 欧美激情极品国产一区二区三区| 国产成人精品久久二区二区91| 99热网站在线观看| 亚洲精品久久午夜乱码| 国产精品一二三区在线看| 在线精品无人区一区二区三| 色婷婷av一区二区三区视频| 亚洲精品一卡2卡三卡4卡5卡 | 国产亚洲欧美在线一区二区| 国产男人的电影天堂91| 久久精品久久久久久久性| 99国产精品免费福利视频| 美女主播在线视频| 国产精品久久久久成人av| 欧美97在线视频| 麻豆av在线久日| 国产日韩一区二区三区精品不卡| 亚洲 国产 在线| 精品国产一区二区三区久久久樱花| 首页视频小说图片口味搜索 | 欧美精品一区二区免费开放| 一级毛片我不卡| 久久精品熟女亚洲av麻豆精品| 久久天堂一区二区三区四区| 在线观看免费日韩欧美大片| 一本久久精品| 熟女av电影| 69精品国产乱码久久久| 欧美性长视频在线观看| 一本一本久久a久久精品综合妖精| 操出白浆在线播放| 男人爽女人下面视频在线观看| 亚洲人成网站在线观看播放| 亚洲人成电影免费在线| a 毛片基地| 亚洲国产精品999| 99热国产这里只有精品6| 日韩免费高清中文字幕av| av有码第一页| 亚洲精品国产一区二区精华液| 日韩中文字幕欧美一区二区 | 曰老女人黄片| 91精品三级在线观看| 精品熟女少妇八av免费久了| 老司机影院成人| 久久久欧美国产精品| 香蕉国产在线看| 国产免费一区二区三区四区乱码| 男女无遮挡免费网站观看| 欧美日韩成人在线一区二区| 一级毛片我不卡| 一区福利在线观看| 亚洲欧洲国产日韩| 国产黄色视频一区二区在线观看| 成人免费观看视频高清| 国产精品免费大片| 天天躁夜夜躁狠狠躁躁| 男女之事视频高清在线观看 | 丁香六月天网| 久久久久久久久久久久大奶| av国产精品久久久久影院| 精品亚洲乱码少妇综合久久| 亚洲精品自拍成人| av国产久精品久网站免费入址| 久久狼人影院| 免费在线观看黄色视频的| 又黄又粗又硬又大视频| 亚洲精品美女久久久久99蜜臀 | 欧美黑人精品巨大| 国产成人av教育| 中文乱码字字幕精品一区二区三区| 亚洲精品一二三| 精品熟女少妇八av免费久了| 日本欧美国产在线视频| 欧美日本中文国产一区发布| av网站在线播放免费| 性色av一级| 黑丝袜美女国产一区| www.精华液| 精品国产超薄肉色丝袜足j| 亚洲成色77777| 高清av免费在线| 国产成人欧美| 亚洲av日韩在线播放| 久久综合国产亚洲精品| 精品一品国产午夜福利视频| 国产麻豆69| 水蜜桃什么品种好| 国产三级黄色录像| 亚洲精品国产色婷婷电影| 91字幕亚洲| 男的添女的下面高潮视频| av国产精品久久久久影院| 亚洲国产毛片av蜜桃av| 日本av免费视频播放| 一区二区三区乱码不卡18| 午夜免费男女啪啪视频观看| 少妇人妻久久综合中文| 蜜桃国产av成人99| 亚洲黑人精品在线| 亚洲精品av麻豆狂野| 亚洲精品久久成人aⅴ小说| 免费观看a级毛片全部| 亚洲自偷自拍图片 自拍| 欧美黑人欧美精品刺激| www.精华液| 自线自在国产av| 1024视频免费在线观看| 免费日韩欧美在线观看| 日韩熟女老妇一区二区性免费视频| √禁漫天堂资源中文www| 成年av动漫网址| 久久天躁狠狠躁夜夜2o2o | 交换朋友夫妻互换小说| 亚洲欧美一区二区三区久久| 黄色视频不卡| 成人18禁高潮啪啪吃奶动态图| 一个人免费看片子| 一区二区av电影网| 老汉色∧v一级毛片| 亚洲av日韩在线播放| 亚洲,一卡二卡三卡| 成年动漫av网址| bbb黄色大片| 99热全是精品| 天天操日日干夜夜撸| 啦啦啦中文免费视频观看日本| 国产又爽黄色视频| 久久精品久久久久久噜噜老黄| 欧美精品啪啪一区二区三区 | 亚洲,欧美,日韩| 成人国语在线视频| 久久久久久人人人人人| 十八禁人妻一区二区| 欧美日韩国产mv在线观看视频| 精品人妻在线不人妻| 国产精品秋霞免费鲁丝片| 亚洲欧美中文字幕日韩二区| 黄色视频不卡| 亚洲综合色网址| 啦啦啦啦在线视频资源| 亚洲av欧美aⅴ国产| 亚洲成人手机| 国产三级黄色录像| 午夜免费鲁丝| 亚洲,欧美精品.| 少妇的丰满在线观看| 女人被躁到高潮嗷嗷叫费观| 汤姆久久久久久久影院中文字幕| 中文乱码字字幕精品一区二区三区| 成人免费观看视频高清| 亚洲欧美日韩高清在线视频 | av线在线观看网站| 两性夫妻黄色片| 国产片特级美女逼逼视频| 色94色欧美一区二区| 青春草亚洲视频在线观看| av在线老鸭窝| 亚洲欧洲国产日韩| 亚洲国产成人一精品久久久| 亚洲人成电影观看| 黄色视频在线播放观看不卡| 久久99一区二区三区| 美女福利国产在线| 黄色一级大片看看| 久久99热这里只频精品6学生| 欧美变态另类bdsm刘玥| 日本猛色少妇xxxxx猛交久久| 极品少妇高潮喷水抽搐| 一区二区av电影网| 好男人电影高清在线观看| 国产精品成人在线| 亚洲熟女毛片儿| 大型av网站在线播放| netflix在线观看网站| 99国产精品免费福利视频| 欧美另类一区| 国产免费又黄又爽又色| 可以免费在线观看a视频的电影网站| 国产爽快片一区二区三区| 亚洲人成电影观看| 国产亚洲av高清不卡| 久久99一区二区三区| 91九色精品人成在线观看| 亚洲精品美女久久久久99蜜臀 | 人人妻人人澡人人爽人人夜夜| 天堂8中文在线网| 久久久欧美国产精品| 大片电影免费在线观看免费| 国产xxxxx性猛交| 精品久久久久久电影网| 国产在线一区二区三区精| 国产熟女午夜一区二区三区| 另类精品久久| 国产一区二区三区av在线| 狂野欧美激情性bbbbbb| 51午夜福利影视在线观看| 亚洲欧美一区二区三区久久| 成在线人永久免费视频| 午夜91福利影院| 美女高潮到喷水免费观看| a 毛片基地| 51午夜福利影视在线观看| 一级毛片我不卡| 久久精品人人爽人人爽视色| 99国产精品99久久久久| videosex国产| 色播在线永久视频| 国产精品一区二区在线观看99| xxxhd国产人妻xxx| 亚洲精品美女久久久久99蜜臀 | 国产一区二区三区综合在线观看| 亚洲欧美日韩高清在线视频 | h视频一区二区三区| 久久久久国产一级毛片高清牌| 97在线人人人人妻| 亚洲免费av在线视频| 亚洲av在线观看美女高潮| 这个男人来自地球电影免费观看| 国产亚洲av高清不卡| 母亲3免费完整高清在线观看| 我要看黄色一级片免费的| 大陆偷拍与自拍| 国产成人免费无遮挡视频| 欧美黑人欧美精品刺激| 精品少妇久久久久久888优播| 欧美激情 高清一区二区三区| 精品一区在线观看国产| 精品一区二区三卡| 午夜久久久在线观看| 国产精品久久久av美女十八| 精品久久蜜臀av无| 国产真人三级小视频在线观看| 少妇的丰满在线观看| 无遮挡黄片免费观看| 七月丁香在线播放| 无遮挡黄片免费观看| 性少妇av在线| 人人妻,人人澡人人爽秒播 | 啦啦啦啦在线视频资源| 麻豆国产av国片精品| 国产精品久久久久久精品电影小说| 国产欧美亚洲国产| 一级毛片女人18水好多 | 一本综合久久免费| 免费女性裸体啪啪无遮挡网站| 亚洲国产最新在线播放| 精品国产一区二区三区久久久樱花| 亚洲欧美日韩高清在线视频 | 久热这里只有精品99| 午夜激情久久久久久久| 久久精品熟女亚洲av麻豆精品| 爱豆传媒免费全集在线观看| 国产又爽黄色视频| 亚洲欧洲日产国产| 国产精品香港三级国产av潘金莲 | 亚洲精品美女久久久久99蜜臀 | 成年美女黄网站色视频大全免费| 中文乱码字字幕精品一区二区三区| 久久精品国产亚洲av高清一级| 我要看黄色一级片免费的| 免费女性裸体啪啪无遮挡网站| 亚洲自偷自拍图片 自拍| 中文字幕另类日韩欧美亚洲嫩草| 超碰成人久久| 国产成人精品无人区| 手机成人av网站| a级毛片在线看网站| 97精品久久久久久久久久精品| 中国美女看黄片| 99热全是精品| 麻豆国产av国片精品| 国产高清视频在线播放一区 | 99国产精品99久久久久| 大陆偷拍与自拍| 午夜福利影视在线免费观看| 国产成人精品久久二区二区91| 9191精品国产免费久久| 欧美变态另类bdsm刘玥| 亚洲av日韩在线播放| 日韩人妻精品一区2区三区| 亚洲,一卡二卡三卡| 女人高潮潮喷娇喘18禁视频| 国产精品偷伦视频观看了| 2021少妇久久久久久久久久久| 久久亚洲精品不卡| 欧美成狂野欧美在线观看| 无限看片的www在线观看| 亚洲男人天堂网一区| 女人精品久久久久毛片| 又紧又爽又黄一区二区| 亚洲av电影在线进入| 国产片内射在线| 日本五十路高清| 久久精品成人免费网站| 精品久久久精品久久久| 高清视频免费观看一区二区| 国产一区亚洲一区在线观看| 亚洲精品中文字幕在线视频| 亚洲一码二码三码区别大吗| 啦啦啦视频在线资源免费观看| 黄色毛片三级朝国网站| 青青草视频在线视频观看| 嫩草影视91久久| 九草在线视频观看| 欧美黄色淫秽网站| 亚洲精品国产av蜜桃| 亚洲欧美成人综合另类久久久| 一区二区日韩欧美中文字幕| 免费在线观看黄色视频的| 国产亚洲一区二区精品| 人成视频在线观看免费观看| 久久99热这里只频精品6学生| 热99久久久久精品小说推荐| av网站在线播放免费| 精品人妻一区二区三区麻豆| av福利片在线| 热re99久久精品国产66热6| 51午夜福利影视在线观看| 欧美日韩国产mv在线观看视频| 岛国毛片在线播放| 国产成人av激情在线播放| 蜜桃在线观看..| 多毛熟女@视频| 国产日韩欧美在线精品| 国产成人免费观看mmmm| 久久人人97超碰香蕉20202| 中文欧美无线码| 国产视频一区二区在线看| 交换朋友夫妻互换小说| 日韩 欧美 亚洲 中文字幕| av福利片在线| 欧美日韩黄片免| 国产精品国产三级专区第一集| 国产精品久久久久成人av| 国产欧美亚洲国产| av网站在线播放免费| 精品亚洲成a人片在线观看| 又大又黄又爽视频免费| 老汉色∧v一级毛片| 国产黄频视频在线观看| 一级黄色大片毛片| 日韩人妻精品一区2区三区| 成年人午夜在线观看视频| 欧美xxⅹ黑人| 亚洲精品国产av成人精品| 久久女婷五月综合色啪小说| 黄色毛片三级朝国网站| 国产在线视频一区二区| 欧美日韩视频高清一区二区三区二| 蜜桃在线观看..| 男女国产视频网站| 80岁老熟妇乱子伦牲交| 国产激情久久老熟女| 亚洲欧美色中文字幕在线| 国产精品亚洲av一区麻豆| 女人高潮潮喷娇喘18禁视频| 亚洲国产av新网站| 伊人久久大香线蕉亚洲五| 女警被强在线播放| 纯流量卡能插随身wifi吗| 别揉我奶头~嗯~啊~动态视频 | 一二三四在线观看免费中文在| 人成视频在线观看免费观看| 亚洲精品日本国产第一区| 久久人妻熟女aⅴ| 日韩免费高清中文字幕av| 国产av精品麻豆| 精品国产国语对白av| 免费高清在线观看日韩| 亚洲成人手机| 亚洲人成77777在线视频| 丝袜美足系列| 国产男女超爽视频在线观看| 99久久99久久久精品蜜桃| 看十八女毛片水多多多| 少妇猛男粗大的猛烈进出视频| 国产一区二区激情短视频 | 女性被躁到高潮视频| 高清黄色对白视频在线免费看| 在现免费观看毛片| 成年美女黄网站色视频大全免费| 超碰成人久久| 亚洲一码二码三码区别大吗| 手机成人av网站| 天天躁夜夜躁狠狠久久av| 黄色视频在线播放观看不卡| 国产精品 国内视频| 亚洲精品第二区| 十八禁人妻一区二区| 建设人人有责人人尽责人人享有的| 亚洲av国产av综合av卡| 日韩中文字幕欧美一区二区 | 婷婷成人精品国产| 69精品国产乱码久久久| 午夜视频精品福利| 精品国产乱码久久久久久小说| 成年av动漫网址| 少妇人妻久久综合中文| 国产亚洲av片在线观看秒播厂| 亚洲av男天堂| 国产亚洲精品第一综合不卡| 中文字幕色久视频| 日本欧美视频一区| 国产一区亚洲一区在线观看| 一个人免费看片子| 成人三级做爰电影| 国产av一区二区精品久久| 国产片内射在线| 久久av网站| 欧美日韩视频高清一区二区三区二| 在线精品无人区一区二区三| 一区二区三区四区激情视频| 日韩熟女老妇一区二区性免费视频| 欧美日韩国产mv在线观看视频| 高清黄色对白视频在线免费看| 亚洲欧美中文字幕日韩二区| 亚洲精品自拍成人| 热re99久久精品国产66热6| 美女主播在线视频| 久久精品国产a三级三级三级| 久久久久久亚洲精品国产蜜桃av| 巨乳人妻的诱惑在线观看| 国产精品亚洲av一区麻豆| 男人添女人高潮全过程视频| 丰满人妻熟妇乱又伦精品不卡| 免费在线观看完整版高清| 国产99久久九九免费精品| 久久女婷五月综合色啪小说| 无限看片的www在线观看| 纵有疾风起免费观看全集完整版| 久久久久国产精品人妻一区二区| 青草久久国产| 成年av动漫网址| 在线av久久热| h视频一区二区三区| 青春草亚洲视频在线观看| 99国产精品免费福利视频| 国产成人a∨麻豆精品| 亚洲人成77777在线视频| 2021少妇久久久久久久久久久| 男女免费视频国产| 亚洲欧美一区二区三区国产| 十八禁人妻一区二区| 免费在线观看影片大全网站 | 91九色精品人成在线观看| 在线精品无人区一区二区三| 啦啦啦在线免费观看视频4| 黄色视频不卡| 这个男人来自地球电影免费观看| 亚洲av在线观看美女高潮| 日韩视频在线欧美| 亚洲欧洲国产日韩| 一本综合久久免费| 国产成人一区二区在线| 国产视频首页在线观看| 国产精品久久久人人做人人爽| 婷婷色综合大香蕉| 久久精品久久久久久久性| 母亲3免费完整高清在线观看| 日韩免费高清中文字幕av| 久久久久网色| 手机成人av网站| 最近最新中文字幕大全免费视频 | 午夜两性在线视频| 精品少妇黑人巨大在线播放| 国产亚洲欧美精品永久| 老汉色∧v一级毛片| 日韩av在线免费看完整版不卡| 欧美日韩国产mv在线观看视频| 国产在线视频一区二区| 免费久久久久久久精品成人欧美视频| 在线观看人妻少妇| 后天国语完整版免费观看| 另类精品久久| 欧美精品人与动牲交sv欧美| 色婷婷久久久亚洲欧美| 色94色欧美一区二区| 久久国产亚洲av麻豆专区| 国产午夜精品一二区理论片| 多毛熟女@视频| 国产精品久久久av美女十八| 亚洲欧洲国产日韩| 91字幕亚洲| 视频在线观看一区二区三区| 欧美人与性动交α欧美软件| 国产高清不卡午夜福利| 热re99久久国产66热| 波野结衣二区三区在线| 亚洲精品日本国产第一区| 国产欧美亚洲国产| 精品少妇内射三级| 免费高清在线观看日韩|