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

    菜豆普通細(xì)菌性疫病抗性基因定位

    2017-02-05 13:57:36朱吉風(fēng)王蘭芬朱振東王述民
    作物學(xué)報(bào) 2017年1期
    關(guān)鍵詞:菜豆抗病細(xì)菌性

    朱吉風(fēng) 武 晶 王蘭芬 朱振東 王述民

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

    菜豆普通細(xì)菌性疫病抗性基因定位

    朱吉風(fēng) 武 晶 王蘭芬 朱振東 王述民*

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

    菜豆普通細(xì)菌性疫病是世界上危害普通菜豆生產(chǎn)的最嚴(yán)重的病害之一。龍蕓豆5號(hào)是我國黑龍江省的主栽品種, 對(duì)菜豆普通細(xì)菌性疫病表現(xiàn)出良好的抗性。為定位來源于龍蕓豆5號(hào)的抗性基因, 本研究構(gòu)建了包含785個(gè)單株的F2分離群體?;谠撊后w構(gòu)建了一張包含206個(gè)SSR標(biāo)記, 總長度1648.42 cM, 標(biāo)記間平均遺傳距離8.00 cM的遺傳圖譜。圖譜包含12個(gè)連鎖群, 各連鎖群平均長度137.37 cM, 連鎖群上標(biāo)記數(shù)量3~35個(gè)。結(jié)合溫室表型鑒定結(jié)果, 采用QTL IciMapping v4.0軟件的完備區(qū)間作圖法進(jìn)行QTL定位和效應(yīng)估計(jì)。接種14 d后在Pv06染色體上檢測(cè)到一個(gè)抗病QTL。該位點(diǎn)位于標(biāo)記p6s249與p6s183之間, 加性效應(yīng)值為0.44, 說明增效基因來源于龍蕓豆5號(hào), LOD值為5.93, 表型貢獻(xiàn)率為4.61%, 該抗病QTL的效應(yīng)值相對(duì)較低, 將在培育穩(wěn)定持久的抗菜豆普通細(xì)菌性疫病的品種中發(fā)揮作用。最后, 對(duì)抗性基因緊密連鎖的11對(duì)SSR引物與菜豆普通細(xì)菌性疫病抗性的關(guān)聯(lián)分析表明, SSR標(biāo)記p6s249與菜豆普通細(xì)菌性疫病抗性極顯著關(guān)聯(lián)(P<0.001), 該標(biāo)記可用于抗病分子育種。

    普通菜豆; 普通細(xì)菌性疫病; SSR; QTL; 關(guān)聯(lián)分析

    普通菜豆(Phaseolus vulgaris L.)是人類生活中重要的食用豆類之一, 富含蛋白質(zhì)、碳水化合物、維生素等人體所必需的營養(yǎng)成分[1]。2014年全世界普通菜豆籽粒產(chǎn)量為2500萬噸, 占食用豆總產(chǎn)量的50%左右(FAOSTAT 2015)。然而, 由地毯草黃單胞桿菌菜豆致病變種(Xanthomonas axonopodis pv. phaseoli, Xap)[2]所引起的菜豆普通細(xì)菌性疫病(common bacterial blight, CBB)嚴(yán)重威脅著普通菜豆的生產(chǎn), 已成為全世界普通菜豆生產(chǎn)中的主要病害[3]。菜豆普通細(xì)菌性疫病為種傳細(xì)菌性病害, 在普通菜豆全生育期內(nèi)皆可發(fā)生, 危害寄主植物的莖、葉、嫩莢和種子, 致使產(chǎn)量及品質(zhì)下降; 由該病導(dǎo)致的普通菜豆產(chǎn)量損失一般為 20%~60%, 嚴(yán)重時(shí)高達(dá)80%, 甚至絕產(chǎn)[4]。這就迫切要求我們鑒定抗菜豆普通細(xì)菌性疫病的優(yōu)異種質(zhì), 并從中挖掘優(yōu)異基因或有效功能標(biāo)記, 將其有效應(yīng)用于普通菜豆抗病分子育種。

    迄今為止, 國內(nèi)外已報(bào)道定位了24個(gè)抗菜豆普通細(xì)菌性疫病的QTL[3,5-9], 分布于普通菜豆的11條連鎖群上[3,6,9]。其中, 位于染色體Pv07上的菜豆普通細(xì)菌性疫病抗病位點(diǎn)最多, 有 5個(gè)[10-12]; 而在染色體Pv01[11]、Pv04[7]、Pv06[7]和Pv10[10-11]上各檢測(cè)到的一個(gè)抗病位點(diǎn)。在已定位的抗病位點(diǎn)中, 對(duì)SAP6[13]、BC420[7]和SU91[14]的研究較為深入。SAP6來源于普通菜豆(Montana 5)[13], 而BC420、SU91均來源于普通菜豆的近緣種寬葉菜豆(PI 319443)[15]。SAP6位于染色體Pv10[13], 可解釋35%的表型變異[13]。BC420位于染色體Pv06[7], 可解釋62%~63%的表型變異[12]。SU91位于染色體 Pv08[14], 可解釋 14%~17%的表型變異[16]。由于已報(bào)道的24個(gè)抗性基因的分子標(biāo)記多為顯性標(biāo)記如RAPD、SCAR標(biāo)記等, 不能鑒別雜合子和純合子, 限制了其在抗病分子育種中的應(yīng)用。在已定位的抗病基因中, 僅針對(duì)BC420和SU91基因開發(fā)了共顯性標(biāo)記且被廣泛應(yīng)用于普通菜豆的抗病分子育種[17]。為此, 需要定位克隆新的菜豆普通細(xì)菌性疫病抗性基因, 并開發(fā)共顯性標(biāo)記, 為抗病分子輔助選擇育種提供更多、更有效的分子標(biāo)記。

    普通菜豆龍蕓豆5號(hào)在溫室和田間對(duì)我國菜豆普通細(xì)菌性疫病病菌Xap菌株XS2均表現(xiàn)出較好的抗性。目前, 關(guān)于龍蕓豆5號(hào)的抗病性遺傳分析研究尚未報(bào)道。本研究利用龍蕓豆5號(hào)與感病材料龍蕓豆4號(hào)配置的 F2雜交群體, 以 SSR分子標(biāo)記定位來源于龍蕓豆5號(hào)的抗病基因, 并為抗病分子育種提供有效的分子標(biāo)記。

    1 材料與方法

    1.1 試驗(yàn)材料

    普通菜豆抗病品種龍蕓豆5號(hào)(龍5)和感病品種龍蕓豆4號(hào)(龍4)為黑龍江省農(nóng)業(yè)科學(xué)院的2個(gè)育成品種, 由中國農(nóng)業(yè)科學(xué)院作物科學(xué)研究所提供。以龍5為父本、龍4為母本, 通過雜交和自交獲得F2分離群體共785個(gè)單株, 用于表型鑒定、遺傳分析和抗病QTL定位。

    取保存于-80℃超低溫冰箱的Xap菌株XS2于牛奶吐溫培養(yǎng)基[18]上活化 4~5 d (28℃±2℃), 然后轉(zhuǎn)入牛肉膏蛋白胨液體培養(yǎng)基中振蕩培養(yǎng)48 h (200轉(zhuǎn)min-1, 28℃±2℃), 用滅菌蒸餾水稀釋至菌液濃度為1×108cfu mL-1, 用于接種。

    1.2 表型鑒定

    以親本材料龍5為抗病對(duì)照、龍4為感病對(duì)照,將親本與F2群體播種于口徑為23 cm × 18 cm × 18 cm的花盆中, 每盆5粒, 每親本材料4盆, 控制溫室白天溫度為28℃±2℃、晚上溫度為20℃±2℃, 待第一片三出復(fù)葉完全展開后采用針刺葉片法接種[19]。接種 14 d后調(diào)查接種植株葉片的發(fā)病情況, 參照Zapata等[20]評(píng)級(jí)標(biāo)準(zhǔn), 調(diào)查記錄每株接種葉片的發(fā)病嚴(yán)重度, 并計(jì)算親本龍5和龍4的平均發(fā)病級(jí)別,利用軟件SAS 9.1對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析。

    1.3 基因組DNA提取及PCR擴(kuò)增

    采用改良CTAB法[21]提取龍5、龍4及F2單株的基因組DNA。用于遺傳分析的SSR引物包括兩部分, 共計(jì)3186對(duì)。第1部分來源于已發(fā)表文獻(xiàn)[22-23],計(jì)381對(duì); 第2部分是本實(shí)驗(yàn)室基于普通菜豆全基因組序列[24]所開發(fā), 計(jì)2805對(duì)。上述SSR引物全部用于親本間的多態(tài)性分析, 選取在親本間有明顯多態(tài)性且?guī)颓宄囊锓治鯢2群體的遺傳連鎖性。

    PCR擴(kuò)增體系為15 μL, 含20 ng模板DNA、0.2 μmol L-1正反引物(Invitrogen, USA)、0.25 mmol L-1dNTPs (dATP、dCTP、dGTP與dTTP)、1.5 μL的10×Taq緩沖液(含1.5 mmol L-1Mg2+)以及1 U DNA聚合酶。PCR擴(kuò)增程序?yàn)?5℃預(yù)變性5 min; 95℃變性30 s, 53℃退火45 s, 72℃延伸45 s, 循環(huán)35次; 最后72℃延伸10 min。PCR擴(kuò)增產(chǎn)物經(jīng)8%非變性聚丙烯酰胺凝膠電泳檢測(cè), 銀染法[25]染色后觀察記錄結(jié)果。

    1.4 QTL定位與關(guān)聯(lián)性分析

    利用QTL IciMapping v4.0[26]構(gòu)建F2群體的遺傳連鎖圖譜, 設(shè)置LOD≥3.0、最大遺傳距離30 cM, 選用Kosambi作圖函數(shù)[27]作圖。參照普通菜豆全基因組序列[24], 將連鎖群與染色體一一對(duì)應(yīng)。

    結(jié)合構(gòu)建的遺傳連鎖圖譜和表型鑒定結(jié)果, 采用完備區(qū)間作圖法(inclusive composite interval mapping, ICIM)[28]檢測(cè)龍5中LOD值≥3.0的抗性基因位點(diǎn)。

    利用 TASSEL 2.1軟件[29]中一般線性模型(general linear model, GLM)對(duì)目標(biāo)區(qū)段所在染色體上的SSR標(biāo)記與CBB抗感結(jié)果進(jìn)行標(biāo)記-性狀關(guān)聯(lián)分析, 結(jié)合QTL檢測(cè)結(jié)果確定與抗病性緊密關(guān)聯(lián)的位點(diǎn)。

    2 結(jié)果與分析

    2.1 表型鑒定結(jié)果

    接種14 d后, 抗病親本龍5接菌葉片沒有出現(xiàn)感病癥狀, 表現(xiàn)出較強(qiáng)抗性(圖 1), 平均發(fā)病級(jí)別為3.88 (圖2); 感病親本龍4接菌葉片的接種區(qū)和對(duì)照區(qū)均枯萎變色(圖 1), 平均發(fā)病級(jí)別為 7.64 (圖 2),表現(xiàn)為感病。F2群體表現(xiàn)出明顯的抗感分離, 最大值超過了高值親本, 最小值低于低值親本, 表明性狀同時(shí)具有正向和負(fù)向超親優(yōu)勢(shì), 且呈連續(xù)分布(圖2), 其偏度(Skewness)值為 0.46, 峰度(Kurtosis)值為-0.30, 二者的絕對(duì)值均小于 0.5, 符合正態(tài)分布, 適合進(jìn)行QTL定位分析。

    圖1 接種Xap菌株XS214 d后龍5與龍4的發(fā)病情況Fig. 1 Phenotype of Long 5 and Long 4 inoculated with Xap strain, XS2

    圖2 接種Xap菌株XS214 d后F2群體的抗性等級(jí)分布Fig. 2 Frequency distribution of ratings in F2population for resistance to Xap strain, XS2

    2.2 遺傳連鎖分析

    利用3186對(duì)SSR引物擴(kuò)增親本龍5和龍4基因組 DNA表明, 有 228對(duì)引物在親本間具有多態(tài)性。利用這些多態(tài)性引物對(duì)785個(gè)F2單株進(jìn)行遺傳連鎖分析, 獲得含206個(gè)SSR標(biāo)記的遺傳連鎖圖譜,共分12個(gè)連鎖群, 其中最短連鎖群為Pv02, 僅含3個(gè)標(biāo)記、最長連鎖群為Pv08, 含35個(gè)標(biāo)記(圖3)。依據(jù)標(biāo)記的染色體信息, 確定每個(gè)連鎖群對(duì)應(yīng)的染色體, 染色體Pv03由于標(biāo)記密度不夠而斷開為2個(gè)連鎖群(Pv03a與Pv03b), 其余連鎖群與其染色體一一對(duì)應(yīng)(圖3)。該遺傳圖譜總長度為1648.42 cM, 平均遺傳距離為8.00 cM; 標(biāo)記間最短遺傳距離為0.60 cM, 位于Pv06上的p6s267與p6s123標(biāo)記之間; 標(biāo)記間最大遺傳距離為 36.92 cM, 位于 Pv06上的p6s126與 p6s277標(biāo)記之間(圖 3); 標(biāo)記密度最大的是 Pv04連鎖群, 其標(biāo)記間平均遺傳距離只有5.65 cM。

    2.3 抗病QTL定位與關(guān)聯(lián)性分析

    結(jié)合遺傳連鎖圖譜和 F2表型鑒定結(jié)果, 利用軟件QTL IciMapping v4.0進(jìn)行ICIM定位分析, 接種后14 d僅在Pv06染色體上檢測(cè)到一個(gè)抗病QTL, 該位點(diǎn)位于分子標(biāo)記 p6s249與 p6s183之間(圖 4-A),其LOD值為5.93 (圖4-B), 共解釋4.61%的表型變異。此外, 在 Pv06染色體上檢測(cè)到的抗病 QTL加性效應(yīng)值為0.44, 說明增效基因來源于龍蕓豆5號(hào),且該抗病QTL能增強(qiáng)龍蕓豆5號(hào)對(duì)菜豆普通細(xì)菌性疫病的抗性。

    利用GLM法對(duì)抗病QTL所在連鎖群(Pv06)上的所有標(biāo)記關(guān)聯(lián)性分析表明, 該連鎖群上有 11對(duì)SSR引物與 CBB發(fā)病等級(jí)顯著相關(guān)(P<0.05), 其中個(gè)SSR標(biāo)記p6s302、p6s281、p6s267、p6s249、p6s210、PVBR163與CBB抗性的相關(guān)性達(dá)到P<0.001的顯著水平(表 1)。結(jié)合本研究所檢測(cè)到的抗性 QTL的側(cè)翼標(biāo)記(p6s249與p6s183)(圖4-A), 發(fā)現(xiàn)與菜豆普通細(xì)菌性疫病抗性位點(diǎn)緊密關(guān)聯(lián)的標(biāo)記為 p6s249, 該引物在G19833基因組中預(yù)擴(kuò)增目標(biāo)片段為216 bp (表1)。p6s249標(biāo)記在龍5基因組中擴(kuò)增出的目標(biāo)片段位于201~238 bp之間, 在龍4擴(kuò)增出的片段位于238~242 bp之間(圖5)。

    3 討論

    由地毯草黃單胞桿菌菜豆致病變種[2]引起的菜豆普通細(xì)菌性疫病是影響全世界菜豆生產(chǎn)的一種主要病害[3]。研究表明菜豆普通細(xì)菌性疫病的抗性是受少數(shù)主效基因控制的數(shù)量性狀[8,12]。目前, 已報(bào)道的24個(gè)與CBB抗性相關(guān)的QTL中, 位于染色體Pv07上有 5個(gè)[10-12], Pv02上 4個(gè)[7,10-11], Pv08上 3個(gè)[5,14], Pv03[7,10]、Pv05[7]、Pv09[11]和 Pv11[10]上各 2 個(gè), Pv01[11]、Pv04[7]、Pv06[7]和Pv10[10-11]上各1個(gè)。這些抗病QTL中, 位于Pv06染色體上的(BC420)因其具有較高貢獻(xiàn)率而得到深入研究[6,17]。

    本研究在Pv06染色體上檢測(cè)到一個(gè)抗Xap菌株XS2的 QTL, 該抗病位點(diǎn)位于分子標(biāo)記 p6s249與p6s183之間(圖4-A), LOD值為5.93 (圖4-B), 貢獻(xiàn)率為 4.61%, 抗病效應(yīng)值較低, 可能將在培育穩(wěn)定持久的抗菜豆普通細(xì)菌性疫病的品種中發(fā)揮作用。依據(jù)SSR引物p6s249與p6s183在G19833物理圖譜上的位置可初步斷定本研究獲得的抗病QTL位于Pv06染色體28~30 Mb的物理距離處(圖4-C)。根據(jù)Shi等[17]所確定的位于 Pv06染色體上的抗性基因BC420的候選基因序列比對(duì)分析, 該抗病位點(diǎn)位于第6染色體4 Mb的物理距離處, 與本研究檢測(cè)到的抗病位點(diǎn)相距較遠(yuǎn)(圖4-C), 因此我們所發(fā)現(xiàn)的抗病位點(diǎn)可能是一個(gè)來源于普通菜豆的新抗病位點(diǎn)。

    圖3 利用F2群體構(gòu)建的SSR分子標(biāo)記遺傳連鎖圖譜Fig. 3 Linkage map of F2population with SSR markers

    圖4 F2群體中抗病QTL的定位Fig. 4 Mapping of QTL for CBB resistance in F2populationB: 紅線代表在Pv06上檢測(cè)到的抗病QTL的LOD值。B: The red line means the LOD score of the major QTL for CBB resistance on chromosome Pv06.

    圖5 p6s249在龍5、龍4和F2基因組中的擴(kuò)增結(jié)果Fig. 5 Amplification products of p6s249 in Long 5, Long 4, and F2populationA: p6s249在F2中擴(kuò)增出的與龍5一致的帶型; B: p6s249在F2中擴(kuò)增出的與龍4一致的帶型; H: p6s249在F2中擴(kuò)增出的雜合帶型; P1: p6s249在龍5中擴(kuò)增出的帶型; P2: p6s249在龍4中擴(kuò)增出的帶型。A: products of p6s249 in F2consistent with Long 5; B: products of p6s249 in F2consistent with Long 4; H: products of the hybrid in F2by p6s249; P1: products of p6s249 in Long 5; P2:products of p6s249 in Long 4.

    表1 SSR引物與CBB抗性間的關(guān)聯(lián)分析Table 1 Testing of association between SSR primers and CBB resistance

    利用與目的基因緊密連鎖的 DNA標(biāo)記進(jìn)行分子標(biāo)記輔助選擇有利于加快作物育種進(jìn)程、提高育種效率[30]。目前, DNA標(biāo)記已被廣泛應(yīng)用于水稻[31]、玉米[32]、小麥[33]等作物的分子育種。在眾多的DNA標(biāo)記中, SSR標(biāo)記因其具有多態(tài)性高、共顯性、數(shù)量豐富、信息量高、重復(fù)性好、DNA需求量低、易檢測(cè)、簡單經(jīng)濟(jì)等優(yōu)點(diǎn), 成為作物分子標(biāo)記輔助選擇育種中最好用的一種標(biāo)記類型[34]。此外, 由DNA的特殊序列標(biāo)記(如RFLP、RAPD等)轉(zhuǎn)換而來的STS、SCAR和SNP等標(biāo)記也同樣被廣泛應(yīng)用于分子標(biāo)記輔助選擇育種[34]。分子標(biāo)記輔助選擇育種技術(shù)在菜豆普通細(xì)菌性疫病抗病育種中也得到應(yīng)用, 特別是3個(gè)重要的SCAR標(biāo)記BC420、SU91和SAP6已被用于菜豆普通細(xì)菌性疫病的 MAS育種[6,35]。然而,由于SCAR標(biāo)記同其他共顯性標(biāo)記一樣不能區(qū)分雜合位點(diǎn), 在一定程度上限制了它的應(yīng)用。Shi等[17]開發(fā)了共顯性候選基因標(biāo)記 BC420-CG14和 SU91-CG11來代替SU91和BC420用于CBB抗病育種。本研究獲得了一個(gè)與CBB極顯著關(guān)聯(lián)的SSR標(biāo)記p6s249, 該標(biāo)記所在位點(diǎn)遠(yuǎn)離 BC420, 為一個(gè)新抗病位點(diǎn)連鎖的標(biāo)記; 此外, 該標(biāo)記來源于普通菜豆,更易于普通菜豆育種選擇, 所以p6s249可能有助于今后的MAS抗病育種。

    4 結(jié)論

    定位了 1個(gè)普通菜豆抗菜豆普通細(xì)菌性疫病QTL, 該位點(diǎn)位于標(biāo)記p6s249與p6s183之間, 共解釋 4.61%的表型變異率, 效應(yīng)值較小, 將在培育穩(wěn)定持久抗菜豆普通細(xì)菌性疫病品種中發(fā)揮作用。p6s249與菜豆普通細(xì)菌性疫病抗性極顯著關(guān)聯(lián)(P<0.001), 可被用于抗病分子育種。

    [1] Bitocchi E, Nanni L, Bellucci E, Rossi M, Giardini A, Zeuli P S, Logozzo G, Stougaard J, McClean P, Attene G, Papa R. Mesoamerican origin of the common bean (Phaseolus vulgaris L.) is revealed by sequence data. Proc Natl Acad Sci USA, 2012, 109: E788-E796

    [2] Schuster M L, Coyne D P. Biology, epidemiology, genetics and breeding for resistance to bacterial pathogens of Phaseolus vulgaris L. Hort Rev, 1981, 3: 28-58

    [3] Singh S P, Schwartz H F. Breeding common bean for resistance to diseases: a review. Crop Sci, 2010, 50: 2199-2223

    [4] Lema-Marquez M, Teran H, Singh S P. Selecting common bean with genes of different evolutionary origins for resistance to Xanthomonas campestris pv. phaseoli. Crop Sci, 2007, 47: 1367-1374

    [5] Jung G, Skroch P W, Nienhuis J, Coyne D P, Arnaud-Santana E, Ariyarathne H M, Marita J M. Confirmation of QTL associated with common bacterial blight resistance in four different genetic backgrounds in common bean. Crop Sci, 1999, 39: 1448-1455

    [6] Miklas P N, Kelly J D, Beebe S E, Blair M W. Common bean breeding for resistance against biotic and abiotic stresses: from classical to MAS breeding. Euphytica, 2006, 147: 105-131

    [7] Tar’an B, Michaels T E, Pauls K P, Mapping genetic factors affecting the reaction to Xanthomonas axonopodis pv. phaseoli in Phaseolus vulgaris L. under field conditions. Genome, 2001, 44: 1046-1056

    [8] Liu S, Yu K, Park S J. Development of STS markers and QTL validation for common bacterial blight resistance in common bean. Plant Breed, 2008, 127: 62-68

    [9] Shi C, Navabi A, Yu K. Association mapping of common bacterial blight resistance QTL in Ontario bean breeding populations. BMC Plant Biol, 2011, 11: 52

    [10] Jung G, Coyne D P, Skroch P W, Nienhuis J, Arnaud-Santana E, Bokosi J, Ariyarathne H M, Steadman J R, Beaver J S, Kaeppler S M. Molecular markers associated with plant architecture and resistance to common blight, web blight, and rust in common beans. J Am Soc Hort Sci, 1996, 121: 794-803

    [11] Ariyarathne H M, Coyne D P, Jung G, Skroch P W, Vidaver A K, Steadman J R, Miklas P N, Bassett M. Molecular mapping of disease resistance genes for halo blight, common bacterial blight, and bean common mosaic virus in a segregating population of common bean. J Am Soc Hort Sci, 1999, 124: 654-662

    [12] Yu K, Park S J, Zhang B, Haffner M, Poysa V. An SSR marker in the nitrate reductase gene of common bean is tightly linked to a major gene conferring resistance to common bacterial blight. Euphytica, 2004, 138: 89-95

    [13] Miklas P N, Coyne D, Grafton K F, Mutlu N, Reiser J, Lindgren D T, Singh S P. A major QTL for common bacterial blight resistance derives from the common bean great northern landrace cultivar Montana No.5. Euphytica, 2003, 131: 137-146

    [14] Pedraza F, Gallego G, Beebe S, Tohme J. Marcadores SCAR y RAPD para la resistencia a la bacteriosis comun (CBB). In: Singh S P, Voysest O, eds. Taller de Mejoramiento de frijol para el Siglo XXI: Bases Para Una Estrategia Para America Latina. CIAT Cali Colombia, 1997. pp 130-134

    [15] Thomas C V, Waines J G. Fertile backcross and allotetraploid plants from crosses between tepary beans and common beans. J Hered, 1984, 75: 93-98.

    [16] Vandemark G J, Fourie D, Miklas P N. Genotyping with real-time PCR reveals recessive epistasis between independent QTL conferring resistance to common bacterial blight in dry bean. Theor Appl Genet, 2008, 117: 513-522

    [17] Shi C, Yu K, Xie W, Perry G, Navabi A, Pauls K P, Miklas P N, Fourie D. Development of candidate gene markers associated to common bacterial blight resistance in common bean. Theor Appl Genet, 2012, 125: 1525-1537

    [18] Sheppard J W, Kurowski C, Remeeus P M. International Rules for seed testing, 7-021: Detection of Xanthomonas axonopodis pv. phaseoli and Xanthomonas axonopodis pv. phaseoli var. fuscans on Phaseolus vulgaris. International Seed Testing Association (ISTA), Bassersdorf, Switzerland. 2007, http://www.seedtest.org/ upload/cms/user/7-021.pdf

    [19] Zapata M. Proposed of a uniform screening procedure for the evaluation of variability of Xanthomonas axonopodis pv. phaseoli and resistance on leaves of Phaseolus vulgaris under greenhouse conditions. Annu Rep Bean Improv Coop, 2006, 49: 213-214

    [20] Zapata M, Beaver J S, Porch T G. Dominant gene for common bean resistance to common bacterial blight caused by Xanthomonas axonopodis pv. phaseoli. Euphytica, 2011, 179: 373-382

    [21] Afanador L K, Hadley S D, Kelly J D. Adoption of a “mini-prep”DNA extraction method for RAPD marker analysis in common bean (Phaseolus vulgaris L.). Annu Rep Bean Improv Coop, 1993, 36: 10-11

    [22] Hanai L R, Santini L, Camargo L E A, Fungaro M H P, Gepts P, Tsai S M, Vieira M L C. Extension of the core map of common bean with EST-SSR, RGA, AFLP, and putative functional markers. Mol Breed, 2010, 25: 25-45

    [23] Chen M, Wu J, Wang L, Zhang X, Blair M W, Jia J, Wang S. Development of mapped simple sequence repeat markers from common bean (Phaseolus vulgaris L.) based on genome sequences of a Chinese landrace and diversity evaluation. Mol Breed, 2014, 33: 489-496

    [24] Schmutz J, McClean P E, Mamidi S, Wu G A, Cannon S B, Grimwood J, Jenkins J, Shu S, Song Q, Chavarro C, Torres-Torres M, Geffroy V, Moghaddam S M, Gao D, Abernathy B, Barry K, Blair M, Brick M A, Chovatia M, Gepts P, Goodstein D M, Gonzales M, Hellsten U, Hyten D L, Jia G, Kelly J D, KudrnaD, Lee R, Richard M M, 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: 706-713

    [25] Zhang C L, Wang Y, Chen H, Lan X Y, Lei C Z. Enhance the efficiency of single-strand conformation polymorphism analysis by short polyacrylamide gel and modified silver staining. Anal Biochem, 2007, 365: 286-287

    [26] 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 populations. Crop J, 2015, 3: 269-283

    [27] Kosambi D D. The estimation of map distances from recombination values. Ann Eugen, 1943, 12: 172-175

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

    [29] Bradbury P, Zhang Z, Kroon D E, Casstevens T M, Ramdoss Y, Buckler E S. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics, 2007, 23: 2633-2635

    [30] Ribaut J M, Hoisington D. Marker-assisted selection: New tools and strategies. Trends Plant Sci, 1998, 3: 236-239

    [31] Toojinda T, Tragoonrung S, Vanavichit A, Siangliw J L, Pa-In N, Jantaboon J, Siangliw M, Fukai S. Molecular breeding for rainfed lowland rice in the Mekong Region. Plant Prod Sci, 2005, 8: 330-333

    [32] Abalo G, Tongoona P, Derera J, Edema R. A comparative analysis of conventional and marker-assisted selection methods in breeding maize streak virus resistance in maize. Crop Sci, 2009, 49: 509-520

    [33] Gupta P K, Langridge P, Mir R R. Marker-assisted wheat breeding: present status and future possibilities. Mol Breed, 2010, 26: 145-161

    [34] Collard B C, Mackill D J. Marker-assisted selection: an approach for precision plant breeding in the twenty-first century. Philos Trans R Soc Lond, 2008, 363: 557-572

    [35] O’Boyle P D, Kelly J D, Kirk W W. Use of marker-assisted selection to breed for resistance to common bacterial blight in common bean. J Am Soc Hort Sci, 2007, 132: 381-386

    Mapping of Common Bacterial Blight Resistance Gene in Common Bean

    ZHU Ji-Feng, WU Jing, WANG Lan-Fen, ZHU Zhen-Dong, and WANG Shu-Min*
    Institute of Crop Science, Chinese Academy of Agricultural Sciences, Beijing 100081, China

    Common bacterial blight disease is a serious disease affecting the production of common bean (Phaseolus vulgaris L.) worldwide. The common bean germplasm Longyundou 5 is the main cultivar carrying common bacterial blight resistance gene in Heilongjiang province. To study the genetic mechanisms behind this resistance, we constructed 785 F2plants from Longyundou 5. Linkage analysis was performed on this population by using SSR markers. A linkage map covering 1648.42 cM with an average marker distance of 8.00 cM among 206 SSR markers was constructed. This map contained 12 linkage groups with a mean group length of 137.37 cM and the number of loci ranging from three to thirty-five. Combined with the results of phenotypic evaluation, QTL analysis was performed by the inclusive composite interval mapping method with QTL IciMapping v4.0. A QTL was detected between p6s249 and p6s183 on chromosome Pv06, with an additive effect of 0.44, which means the favourable gene of this locus is from Longyundou 5. The LOD score of the QTL was 5.93 and the total phenotypic variation at 14 days after inoculation was 4.61%. These results showed the effect of the QTL was low, which will play an important role in breeding durable resistant varieties of common bean. Association analysis between 11 SSR primers linked with common bacterial blight resistance and the disease rating showed that SSR marker p6s249 is excellently associated (P<0.001) with common bacterial blight resistance and can be used for marker assisted selection.

    Phaseolus vulgaris L; Common bacterial blight; SSR; QTL; Association analysis

    10.3724/SP.J.1006.2017.00001

    本研究由國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)(CARS-09), 國家科技支撐計(jì)劃項(xiàng)目(2013BAD01B03-18a)和中國農(nóng)業(yè)科學(xué)院科技創(chuàng)新工程項(xiàng)目資助。

    This study was supported by the China Agriculture Research System (CARS-09), the National Key Technology R&D Program of China (2013BAD01B03-18a), and the Agricultural Science and Technology Innovation Program of the Chinese Academy of Agricultural Sciences.

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

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

    稿日期): 2016-05-04; Accepted(接受日期): 2016-07-11; Published online(

    日期): 2016-07-28.

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

    猜你喜歡
    菜豆抗病細(xì)菌性
    我國小麥基因組編輯抗病育種取得突破
    細(xì)菌性食物中毒事件調(diào)查研究
    基于羊細(xì)菌性疾病的預(yù)防控制分析
    預(yù)防細(xì)菌性食物中毒的主要方法
    bZIP轉(zhuǎn)錄因子在植物激素介導(dǎo)的抗病抗逆途徑中的作用
    葡萄新品種 優(yōu)質(zhì)又抗病
    細(xì)菌性肝膿腫64例診治分析
    菜豆采后生理及貯藏技術(shù)的研究進(jìn)展
    番茄果實(shí)感染灰霉病過程中H2O2的抗病作用
    春播菜豆如何實(shí)現(xiàn)一播全苗
    長江蔬菜(2013年23期)2013-03-11 15:09:39
    国产精品一二三区在线看| 亚洲无线观看免费| 免费大片18禁| 校园人妻丝袜中文字幕| 大片免费播放器 马上看| 菩萨蛮人人尽说江南好唐韦庄| freevideosex欧美| 色综合色国产| 国产成人freesex在线| 国产探花极品一区二区| 亚洲国产日韩欧美精品在线观看| 国产成人午夜福利电影在线观看| 亚洲人成网站高清观看| 国产在视频线在精品| 亚洲人成网站在线播| 亚洲人成网站在线播| 乱码一卡2卡4卡精品| 精品久久久久久成人av| 久热久热在线精品观看| 熟妇人妻不卡中文字幕| 美女高潮的动态| 别揉我奶头 嗯啊视频| 狠狠精品人妻久久久久久综合| 丝瓜视频免费看黄片| 99热这里只有精品一区| 国产 一区 欧美 日韩| 乱码一卡2卡4卡精品| 国产欧美日韩精品一区二区| 国产一区有黄有色的免费视频 | 一区二区三区免费毛片| 91精品伊人久久大香线蕉| 搡老妇女老女人老熟妇| 精品久久久久久久末码| 婷婷色综合www| 亚洲精品国产av成人精品| 久久久久久久久大av| 亚洲精品乱久久久久久| 国产精品爽爽va在线观看网站| 国产黄色小视频在线观看| 午夜老司机福利剧场| 一夜夜www| 久久午夜福利片| 99久国产av精品国产电影| 久99久视频精品免费| 人妻系列 视频| 97精品久久久久久久久久精品| 又爽又黄无遮挡网站| 美女高潮的动态| 亚洲久久久久久中文字幕| 久久久久久久久久久丰满| 亚洲人成网站在线播| 国产色婷婷99| 中文字幕人妻熟人妻熟丝袜美| 最近的中文字幕免费完整| 精品一区二区三区视频在线| 女人被狂操c到高潮| av黄色大香蕉| 国产综合懂色| 国产欧美日韩精品一区二区| 内射极品少妇av片p| 边亲边吃奶的免费视频| 亚洲怡红院男人天堂| 亚洲av成人av| 超碰av人人做人人爽久久| 一个人免费在线观看电影| 69av精品久久久久久| 欧美成人一区二区免费高清观看| 午夜精品在线福利| 国产69精品久久久久777片| 成人亚洲精品av一区二区| 99热6这里只有精品| 美女高潮的动态| 麻豆av噜噜一区二区三区| 久久精品国产鲁丝片午夜精品| 国产成人精品一,二区| 80岁老熟妇乱子伦牲交| 成人综合一区亚洲| 十八禁网站网址无遮挡 | 校园人妻丝袜中文字幕| 欧美性感艳星| 亚洲18禁久久av| 美女国产视频在线观看| 99热网站在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品人妻一区二区三区麻豆| 久久精品久久久久久久性| 亚洲成人一二三区av| 欧美成人精品欧美一级黄| 男女那种视频在线观看| 国产一区二区三区综合在线观看 | 人体艺术视频欧美日本| 精品久久久久久久久av| 午夜激情欧美在线| 人妻制服诱惑在线中文字幕| 最近最新中文字幕免费大全7| 日本欧美国产在线视频| 精品久久久久久久人妻蜜臀av| 国产在视频线在精品| 日本-黄色视频高清免费观看| 免费不卡的大黄色大毛片视频在线观看 | 中文乱码字字幕精品一区二区三区 | 国产不卡一卡二| 日韩av在线大香蕉| 一级爰片在线观看| 蜜臀久久99精品久久宅男| 精品一区在线观看国产| 亚洲国产成人一精品久久久| 成年女人看的毛片在线观看| 国产视频首页在线观看| 日本三级黄在线观看| 一级毛片黄色毛片免费观看视频| 2018国产大陆天天弄谢| 日本黄色片子视频| 夫妻性生交免费视频一级片| 久久久久久久久大av| 哪个播放器可以免费观看大片| 精品久久久久久久久av| 国产精品日韩av在线免费观看| 欧美性感艳星| 天天躁日日操中文字幕| 日本欧美国产在线视频| 秋霞在线观看毛片| 久热久热在线精品观看| 欧美一级a爱片免费观看看| 亚洲人成网站在线观看播放| 久久久久久久久久久免费av| 欧美区成人在线视频| 国产黄色免费在线视频| av免费观看日本| 男女边吃奶边做爰视频| 一边亲一边摸免费视频| 亚洲精品久久午夜乱码| 18禁在线无遮挡免费观看视频| 亚洲精品乱码久久久v下载方式| 十八禁网站网址无遮挡 | 国产男女超爽视频在线观看| 国产av码专区亚洲av| 亚洲精品久久久久久婷婷小说| 国产午夜精品久久久久久一区二区三区| 久久99精品国语久久久| 欧美日韩在线观看h| 在线免费观看的www视频| 国产淫片久久久久久久久| 看黄色毛片网站| 91av网一区二区| 能在线免费看毛片的网站| 国产免费又黄又爽又色| 3wmmmm亚洲av在线观看| 国内精品宾馆在线| 高清视频免费观看一区二区 | 免费人成在线观看视频色| 亚洲av中文av极速乱| 久久精品久久精品一区二区三区| 男女那种视频在线观看| 2021天堂中文幕一二区在线观| 亚洲欧美日韩卡通动漫| 99久久精品一区二区三区| 国产三级在线视频| 国产在视频线精品| 亚洲性久久影院| 亚洲av成人精品一二三区| 日韩亚洲欧美综合| 九草在线视频观看| 亚洲av成人av| 成年av动漫网址| 国产熟女欧美一区二区| 激情五月婷婷亚洲| 直男gayav资源| 国产高潮美女av| 欧美性感艳星| 免费黄频网站在线观看国产| 黄片无遮挡物在线观看| 国产老妇伦熟女老妇高清| 亚洲,欧美,日韩| 亚洲av.av天堂| 精品人妻熟女av久视频| 亚洲色图av天堂| 国产国拍精品亚洲av在线观看| 建设人人有责人人尽责人人享有的 | 午夜福利在线观看吧| 欧美97在线视频| 亚洲精品,欧美精品| 在线观看免费高清a一片| 男女那种视频在线观看| 久久精品国产自在天天线| 一本久久精品| 岛国毛片在线播放| 777米奇影视久久| 久久99热6这里只有精品| 三级男女做爰猛烈吃奶摸视频| 床上黄色一级片| 亚洲国产成人一精品久久久| 亚洲精品日韩在线中文字幕| 久久国内精品自在自线图片| 精品少妇黑人巨大在线播放| av免费在线看不卡| 午夜福利成人在线免费观看| 成人毛片a级毛片在线播放| 亚洲av一区综合| 高清在线视频一区二区三区| av播播在线观看一区| 婷婷色综合www| 尾随美女入室| 麻豆国产97在线/欧美| 久久久久久久大尺度免费视频| 免费观看在线日韩| 男女视频在线观看网站免费| 欧美+日韩+精品| 七月丁香在线播放| 国产精品伦人一区二区| 久久热精品热| 久久午夜福利片| 国产又色又爽无遮挡免| 国产亚洲91精品色在线| 女的被弄到高潮叫床怎么办| 日韩av不卡免费在线播放| 国产精品久久久久久av不卡| 特级一级黄色大片| 色哟哟·www| 少妇的逼水好多| 日韩伦理黄色片| 日本一二三区视频观看| 99re6热这里在线精品视频| 免费电影在线观看免费观看| 中文字幕免费在线视频6| 成人毛片a级毛片在线播放| 在线观看免费高清a一片| 精品久久久久久久人妻蜜臀av| 又大又黄又爽视频免费| 97精品久久久久久久久久精品| www.av在线官网国产| 国产精品国产三级国产av玫瑰| 青春草国产在线视频| 国产精品久久久久久久电影| 91精品伊人久久大香线蕉| 赤兔流量卡办理| 日日撸夜夜添| 亚洲av电影不卡..在线观看| av网站免费在线观看视频 | 久久久久久国产a免费观看| 免费av不卡在线播放| 成年版毛片免费区| 26uuu在线亚洲综合色| 一夜夜www| 五月天丁香电影| 亚洲成人一二三区av| 在线播放无遮挡| 97人妻精品一区二区三区麻豆| 啦啦啦中文免费视频观看日本| 亚洲av不卡在线观看| 最近中文字幕2019免费版| av免费在线看不卡| 久99久视频精品免费| 国产精品不卡视频一区二区| 美女大奶头视频| 我要看日韩黄色一级片| 99热全是精品| 国产激情偷乱视频一区二区| 亚洲欧美精品自产自拍| av又黄又爽大尺度在线免费看| 久久国内精品自在自线图片| 亚洲精品国产av蜜桃| 国产成人freesex在线| 午夜福利在线观看免费完整高清在| 日韩制服骚丝袜av| 亚洲欧美清纯卡通| 久久久久久久亚洲中文字幕| 黄色配什么色好看| 日日啪夜夜撸| 99热网站在线观看| 一个人观看的视频www高清免费观看| 亚洲成人久久爱视频| 久久久久九九精品影院| 美女高潮的动态| 亚洲一区高清亚洲精品| 可以在线观看毛片的网站| 精品人妻视频免费看| 18禁裸乳无遮挡免费网站照片| 青春草视频在线免费观看| 黄片wwwwww| 三级毛片av免费| 在线天堂最新版资源| 亚洲欧美日韩卡通动漫| 日韩国内少妇激情av| 成人特级av手机在线观看| 国产综合懂色| 老女人水多毛片| 精品人妻视频免费看| 欧美xxxx性猛交bbbb| 天堂中文最新版在线下载 | 最近视频中文字幕2019在线8| 免费少妇av软件| 成人国产麻豆网| 又爽又黄无遮挡网站| 亚洲激情五月婷婷啪啪| 亚洲欧美一区二区三区国产| 午夜福利高清视频| 午夜爱爱视频在线播放| 中文天堂在线官网| 简卡轻食公司| ponron亚洲| 99久久精品国产国产毛片| 亚洲av二区三区四区| 国产精品女同一区二区软件| 日韩在线高清观看一区二区三区| 亚洲国产欧美人成| 国产激情偷乱视频一区二区| 国产男人的电影天堂91| av在线老鸭窝| 黑人高潮一二区| 欧美成人午夜免费资源| 免费播放大片免费观看视频在线观看| 婷婷色av中文字幕| 亚洲精品日韩av片在线观看| 一级av片app| 亚洲精品自拍成人| 99热6这里只有精品| 日韩欧美 国产精品| 最后的刺客免费高清国语| 国产精品.久久久| 国产精品女同一区二区软件| 麻豆成人av视频| 夫妻午夜视频| 少妇的逼水好多| 有码 亚洲区| 精品久久久噜噜| 色吧在线观看| 久久久久久久亚洲中文字幕| 精品午夜福利在线看| 久久久久久久久久久丰满| 国产不卡一卡二| 久热久热在线精品观看| 欧美日韩精品成人综合77777| 欧美丝袜亚洲另类| 午夜日本视频在线| 一级黄片播放器| 国产高潮美女av| 男人和女人高潮做爰伦理| 全区人妻精品视频| 亚洲av免费高清在线观看| 日本免费在线观看一区| 男插女下体视频免费在线播放| av福利片在线观看| 欧美成人一区二区免费高清观看| 99视频精品全部免费 在线| 男女视频在线观看网站免费| av天堂中文字幕网| 欧美xxxx黑人xx丫x性爽| 精品酒店卫生间| 亚洲欧美日韩卡通动漫| 久久久久久久久久久免费av| videossex国产| 黄色日韩在线| videossex国产| .国产精品久久| 亚洲高清免费不卡视频| 精品国产露脸久久av麻豆 | 天天一区二区日本电影三级| 亚洲欧美日韩卡通动漫| 亚洲色图av天堂| 国产成人a∨麻豆精品| 欧美不卡视频在线免费观看| 好男人视频免费观看在线| 久久亚洲国产成人精品v| 色综合亚洲欧美另类图片| 国产人妻一区二区三区在| 亚洲精品一二三| 又大又黄又爽视频免费| 最近的中文字幕免费完整| 看十八女毛片水多多多| 亚洲欧洲日产国产| 精品久久久久久成人av| 日韩精品青青久久久久久| 午夜福利视频1000在线观看| 少妇熟女aⅴ在线视频| 国产精品人妻久久久影院| 色视频www国产| 永久免费av网站大全| 国产免费视频播放在线视频 | 国产男人的电影天堂91| 蜜桃久久精品国产亚洲av| 神马国产精品三级电影在线观看| 最后的刺客免费高清国语| 欧美激情国产日韩精品一区| 日韩av免费高清视频| av线在线观看网站| 天堂中文最新版在线下载 | 97人妻精品一区二区三区麻豆| 亚洲精品久久午夜乱码| 成人一区二区视频在线观看| 午夜福利在线观看免费完整高清在| 亚洲自偷自拍三级| 中文字幕亚洲精品专区| 亚洲欧美精品自产自拍| 亚洲熟妇中文字幕五十中出| 亚洲精品日本国产第一区| 精品国内亚洲2022精品成人| 我的老师免费观看完整版| 午夜福利视频精品| 别揉我奶头 嗯啊视频| 男女啪啪激烈高潮av片| 18禁在线无遮挡免费观看视频| 赤兔流量卡办理| 亚洲av.av天堂| 超碰97精品在线观看| 国产高潮美女av| 亚洲精品影视一区二区三区av| 97超碰精品成人国产| 国产 亚洲一区二区三区 | 51国产日韩欧美| h日本视频在线播放| 狂野欧美激情性xxxx在线观看| 最近最新中文字幕大全电影3| 日韩中字成人| 99热全是精品| 麻豆久久精品国产亚洲av| 大片免费播放器 马上看| 亚洲,欧美,日韩| 亚洲欧洲日产国产| 国产 一区 欧美 日韩| 少妇丰满av| a级毛色黄片| 人妻夜夜爽99麻豆av| 老司机影院成人| 久久精品人妻少妇| 久热久热在线精品观看| 美女国产视频在线观看| 久久精品熟女亚洲av麻豆精品 | av黄色大香蕉| 一边亲一边摸免费视频| 91精品一卡2卡3卡4卡| 天美传媒精品一区二区| 人妻制服诱惑在线中文字幕| 一区二区三区乱码不卡18| av在线蜜桃| 少妇人妻精品综合一区二区| 日韩欧美精品v在线| 国产91av在线免费观看| 黑人高潮一二区| 日韩欧美一区视频在线观看 | 亚洲成色77777| 最近手机中文字幕大全| 男女边摸边吃奶| 麻豆av噜噜一区二区三区| 日韩伦理黄色片| 国产大屁股一区二区在线视频| 高清av免费在线| 亚洲国产日韩欧美精品在线观看| 久久韩国三级中文字幕| 身体一侧抽搐| 日韩一区二区三区影片| 久久精品夜夜夜夜夜久久蜜豆| 最近最新中文字幕大全电影3| 少妇的逼好多水| 尾随美女入室| 大片免费播放器 马上看| 欧美xxⅹ黑人| 一级黄片播放器| 国产一区二区在线观看日韩| 男女那种视频在线观看| 日韩强制内射视频| 国产精品.久久久| 丝袜美腿在线中文| 看黄色毛片网站| 亚洲av免费在线观看| 欧美极品一区二区三区四区| 精华霜和精华液先用哪个| 三级国产精品片| 午夜福利高清视频| 亚洲精品影视一区二区三区av| 成人欧美大片| 在线免费观看的www视频| 成年av动漫网址| 日本午夜av视频| a级毛色黄片| 身体一侧抽搐| 欧美zozozo另类| 中文字幕av在线有码专区| 亚洲熟妇中文字幕五十中出| 亚洲天堂国产精品一区在线| 国产一区二区在线观看日韩| 亚洲成人一二三区av| 久久精品人妻少妇| 色视频www国产| 久久这里只有精品中国| 可以在线观看毛片的网站| 少妇裸体淫交视频免费看高清| 三级经典国产精品| 成人av在线播放网站| 国产午夜精品久久久久久一区二区三区| 久久精品国产亚洲av天美| 一级片'在线观看视频| 看黄色毛片网站| 亚洲最大成人av| 国产在线男女| 国产国拍精品亚洲av在线观看| 九九爱精品视频在线观看| 黄片wwwwww| 晚上一个人看的免费电影| 亚洲欧美日韩无卡精品| 国产精品熟女久久久久浪| 国产有黄有色有爽视频| 色播亚洲综合网| 白带黄色成豆腐渣| 午夜福利在线观看吧| 九草在线视频观看| 精品久久国产蜜桃| 国产一区有黄有色的免费视频 | 淫秽高清视频在线观看| 老师上课跳d突然被开到最大视频| 九九久久精品国产亚洲av麻豆| 能在线免费看毛片的网站| 亚洲性久久影院| 国产综合精华液| 久久精品国产亚洲av涩爱| 青春草亚洲视频在线观看| 尾随美女入室| 在线观看美女被高潮喷水网站| 久久久久久久久久成人| 校园人妻丝袜中文字幕| 日韩国内少妇激情av| 国精品久久久久久国模美| 女人久久www免费人成看片| 卡戴珊不雅视频在线播放| 免费大片黄手机在线观看| 欧美激情在线99| av天堂中文字幕网| 精品久久久精品久久久| 欧美三级亚洲精品| 麻豆成人午夜福利视频| 男的添女的下面高潮视频| 少妇人妻精品综合一区二区| 免费高清在线观看视频在线观看| 免费在线观看成人毛片| 久久这里有精品视频免费| 日韩av免费高清视频| 免费人成在线观看视频色| 亚洲精品成人久久久久久| 在线观看一区二区三区| 亚洲精品456在线播放app| 久久久午夜欧美精品| 国产成人aa在线观看| 汤姆久久久久久久影院中文字幕 | 欧美 日韩 精品 国产| 91av网一区二区| 中文字幕免费在线视频6| 婷婷六月久久综合丁香| 亚洲美女搞黄在线观看| www.色视频.com| 一级爰片在线观看| 有码 亚洲区| 亚洲欧美成人精品一区二区| 在线观看美女被高潮喷水网站| 日韩国内少妇激情av| 69av精品久久久久久| 国产老妇女一区| 精品欧美国产一区二区三| 1000部很黄的大片| 亚洲精品第二区| 日本与韩国留学比较| 久久草成人影院| 99久久精品热视频| 国产高清三级在线| 国国产精品蜜臀av免费| 91精品一卡2卡3卡4卡| 亚洲精品久久午夜乱码| 亚洲av二区三区四区| 禁无遮挡网站| 国产乱人视频| 国产v大片淫在线免费观看| 一边亲一边摸免费视频| 成人亚洲精品一区在线观看 | 夫妻午夜视频| 色网站视频免费| a级毛色黄片| 国产精品一区二区三区四区久久| 国产成人精品福利久久| 国产中年淑女户外野战色| 国产 亚洲一区二区三区 | 国产伦理片在线播放av一区| 大香蕉97超碰在线| 最新中文字幕久久久久| 日韩在线高清观看一区二区三区| 久久国内精品自在自线图片| 久久精品国产自在天天线| 中文字幕av在线有码专区| 男插女下体视频免费在线播放| 亚洲成人中文字幕在线播放| 久久久精品免费免费高清| 好男人视频免费观看在线| 街头女战士在线观看网站| 成人美女网站在线观看视频| 在线观看免费高清a一片| 日产精品乱码卡一卡2卡三| 狠狠精品人妻久久久久久综合| 国产成人免费观看mmmm| 国产人妻一区二区三区在| 国产一区有黄有色的免费视频 | 舔av片在线| 搞女人的毛片| 青春草亚洲视频在线观看| av在线亚洲专区| 校园人妻丝袜中文字幕| 搞女人的毛片| 一区二区三区免费毛片| 亚洲国产高清在线一区二区三| 男人爽女人下面视频在线观看| 国产av在哪里看| 欧美成人精品欧美一级黄| 亚洲怡红院男人天堂| 免费人成在线观看视频色| 一级毛片电影观看| 不卡视频在线观看欧美| 国产精品久久视频播放| 日本免费在线观看一区|