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

    長江下游粳稻稻瘟病廣譜抗性基因組合模式分析

    2021-05-12 15:00:44吳云雨肖寧余玲蔡躍潘存紅李育紅張小祥黃年生季紅娟戴正元李愛宏
    中國農(nóng)業(yè)科學(xué) 2021年9期
    關(guān)鍵詞:廣譜抗粳稻稻瘟病

    吳云雨,肖寧,4,余玲,蔡躍,潘存紅,李育紅,張小祥,黃年生,季紅娟,戴正元,李愛宏,3

    長江下游粳稻稻瘟病廣譜抗性基因組合模式分析

    吳云雨1,2,肖寧1,2,4,余玲1,2,蔡躍1,2,潘存紅1,3,李育紅1,2,張小祥1,2,黃年生1,2,季紅娟1,2,戴正元1,3,李愛宏1,2,3

    1江蘇里下河地區(qū)農(nóng)業(yè)科學(xué)研究所,江蘇揚州 225007;2揚州大學(xué)江蘇省作物基因組學(xué)和分子育種重點實驗室,江蘇揚州 225009;3揚州大學(xué)江蘇省糧食作物現(xiàn)代產(chǎn)業(yè)技術(shù)協(xié)同創(chuàng)新中心,江蘇揚州 225009;4中國農(nóng)業(yè)科學(xué)院植物保護研究所植物病蟲害生物學(xué)國家重點實驗室,北京 100193

    【】基因聚合是實現(xiàn)水稻稻瘟病廣譜抗性的有效途徑之一。通過構(gòu)建粳稻背景下不同雙基因聚合系,利用長江下游粳型稻瘟病菌()菌株評價其抗性效應(yīng)并解析其抗性效應(yīng)產(chǎn)生的構(gòu)成因子,為長江下游粳稻抗稻瘟病育種提供廣譜抗性基因組合模式和種質(zhì)資源。以粳稻07GY31為背景的Piz基因座不同復(fù)等位基因(、、、、和)單基因系為核心,利用不完全NCII交配設(shè)計,分別與其他廣譜抗性基因(、和)單基因系雜交,經(jīng)分子標記輔助選擇和農(nóng)藝性狀篩選,共構(gòu)建18種不同基因組合的雙基因聚合系。2019年利用長江下游粳稻種植區(qū)采集、分離的109個稻瘟病代表性菌株進行苗瘟、穗瘟人工接種鑒定及不同病圃的自然誘發(fā)鑒定,評價不同雙基因聚合系的抗性效應(yīng),并分析雙基因聚合系抗性效應(yīng)的構(gòu)成因子。genotyping by sequencing(GBS)分析表明所構(gòu)建的雙基因聚合系均具有較高的背景恢復(fù)率,分布于97.08%(PPL)—99.08%(PPL)。表明除了目標基因區(qū)域不同外,所有雙基因聚合系的遺傳背景幾乎完全與受體親本07GY31一致。同時人工接菌鑒定表明絕大部分雙基因聚合系苗瘟和穗瘟抗性水平都優(yōu)于單基因系。其中苗瘟抗性效應(yīng)較好的聚合系分別為PPL、PPL、PPL、PPL、PPL、PPL、PPL、PPL、PPL, 而穗瘟抗性效應(yīng)較好的聚合系分別為PPL、PPL、PPL、PPL、PPL、PPL、PPL。不同抗性基因聚合后產(chǎn)生不同的效應(yīng),其中互補效應(yīng)高且能有效表達是提高雙基因聚合系苗瘟和穗瘟抗性的關(guān)鍵因子。雙基因聚合系PPL、PPL和PPL在苗瘟和穗瘟的人工接種,以及在不同病圃的自然誘發(fā)鑒定中均表現(xiàn)穩(wěn)定的廣譜抗性,同時,農(nóng)藝性狀調(diào)查結(jié)果也表明這3個雙基因聚合系的基本農(nóng)藝性狀與輪回親本07GY31基本一致,因此,基因組合/、/和/是適于長江下游粳稻的廣譜抗性基因組合模式??剐曰虻慕M合方式影響聚合系的抗性水平,互補效應(yīng)高且能有效表達是粳型雙基因聚合系抗性效應(yīng)提高的關(guān)鍵因子。本研究構(gòu)建的雙基因聚合系及其抗性效應(yīng)分析為長江下游廣譜稻瘟病抗性粳稻品種的精準培育提供了種質(zhì)資源和理論支撐。

    長江下游;粳稻;稻瘟病菌;稻瘟?。换蚓酆?;效應(yīng)分析

    0 引言

    【研究意義】由稻瘟病菌()引起的稻瘟病是威脅世界水稻安全生產(chǎn)的主要真菌類病害之一[1],控制該病害最為經(jīng)濟有效、綠色安全的策略是利用抗病基因培育抗病品種[2],但目前已鑒定的絕大多數(shù)稻瘟病抗性基因為垂直抗性基因,其中部分抗性基因存在不能兼顧苗瘟和穗瘟抗性、抗譜窄等缺點,同時由于稻瘟病菌的致病性分化快、生理小種多、變異頻繁,導(dǎo)致抗病品種應(yīng)用3—5年后便“喪失”抗性變?yōu)楦胁∑贩N??剐曰蚓酆媳蛔C明是培育廣譜抗病品種的主要途徑之一[3-4],然而,不同的抗性基因聚合后會產(chǎn)生不同的抗性效應(yīng)。因此,合理利用不同抗性基因進行聚合并評價其抗性效應(yīng),進而發(fā)掘有效的廣譜抗性基因組合模式對培育廣譜抗病品種、保障水稻安全生產(chǎn)具有重要現(xiàn)實意義和實踐應(yīng)用價值?!厩叭搜芯窟M展】宿主對病原菌的抗性反應(yīng)可分為完全抗性和部分抗性兩種類型,完全抗性通常由單個主效抗病基因控制,而部分抗性通常由多個微效抗性基因控制[5-6]。水稻稻瘟病主效抗病基因大多為來自STAND(signal transduction ATPases with numerous domains)亞家族中的NBS-LRR蛋白家族成員[7],NBS-LRR蛋白以直接或間接方式識別病原菌無毒蛋白后激活下游基因表達從而激發(fā)水稻抗病免疫反應(yīng)[8-11]。目前,已鑒定的稻瘟病主效抗性基因超過100個,部分抗性基因約347個[12-13]。其中有28個主效抗性基因(、、、、、、、、、、、、、、、/、/、、/、、、、、和)[14-16]及5個部分抗性基因(、、、和)[17-21]被成功克隆,為稻瘟病抗性育種提供了有效基因資源,但由于稻瘟病菌群體的高度變異及新的毒性種群的出現(xiàn),許多攜帶單一抗病基因(、、和)的抗性品種在大規(guī)模種植一段時間后失去了抗性效應(yīng)變?yōu)楦胁∑贩N[22-24]。為降低抗性喪失風(fēng)險,利用不同抗病基因進行基因聚合是培育廣譜、持久抗病品種的一個有效途徑。Jiang等[25]和Xiao等[26]分別報道聚合//和/的改良系對中國南方稻區(qū)稻瘟病菌菌株具有廣譜抗性,而聚合則對印度稻瘟病菌群體表現(xiàn)出廣譜抗性[4]。但抗病基因聚合后并不簡單表現(xiàn)出抗性效應(yīng)累加,有的還存在正向甚至負向互作。國際水稻研究所以Co39為背景構(gòu)建攜帶不同抗性基因的聚合系,鑒定結(jié)果顯示Piz/的基因聚合系抗性不如單基因抗性品系[27]。同樣,Xiao等[28]研究發(fā)現(xiàn)在粳稻07GY31背景下雙基因聚合系穗瘟抗性顯著低于單基因系;Wu等[29]以秈稻揚稻6號為背景,構(gòu)建了15種不同的抗性基因聚合系,發(fā)現(xiàn)/、/和/為秈稻背景下的廣譜抗性基因組合模式,同時也發(fā)現(xiàn)部分基因組合產(chǎn)生了負向互作效應(yīng)。此外,由于稻瘟病菌生理小種依賴于水稻宿主生存,生理小種基因型分布與水稻宿主種植區(qū)域有關(guān),來自秈稻、粳稻種植區(qū)域的生理小種被明顯分為“秈稻群”和“粳稻群”[30-31]。因此,在不同種植區(qū)需選擇與當(dāng)?shù)氐疚敛【悍怯H和性的抗性基因進行基因聚合,確保基因聚合后可有效提高其抗性水平。【本研究切入點】長江下游地區(qū)為我國主要粳稻種植區(qū)之一,約占全國粳稻總種植面積的40%[32-33],該地區(qū)粳稻主要含有、、、和等抗性基因[34-35]。近年來,隨著耕作制度和氣候條件的變化,稻瘟病在長江下游粳稻種植區(qū)常態(tài)化發(fā)生。如2014年江蘇水稻稻瘟病發(fā)生面積98萬公頃,2016年上升至128.7萬公頃,且發(fā)病面積呈逐年增加趨勢,水稻生產(chǎn)存在嚴重安全隱患。表明長江下游應(yīng)用品種現(xiàn)有攜帶的抗性基因已遠不能控制稻瘟病的發(fā)生,亟需導(dǎo)入新的抗性基因/基因組合?!緮M解決的關(guān)鍵問題】以粳稻07GY31為背景的Piz基因座不同復(fù)等位基因(、、、、和)單基因系為核心,采用不完全NCII交配設(shè)計,分別與其他廣譜抗性基因(、和)單基因系雜交,經(jīng)過分子標記輔助選擇和農(nóng)藝性狀篩選,構(gòu)建18種不同基因組合的雙基因聚合系。從長江下游不同粳稻種植區(qū)采集稻瘟病菌菌株進行苗瘟和穗瘟的人工接種鑒定,以及開展病圃的自然誘發(fā)鑒定,評價不同雙基因聚合系的抗性效應(yīng),并分析雙基因聚合系抗性效應(yīng)的構(gòu)成因子,從而明確適宜長江下游粳稻的稻瘟病廣譜抗性基因組合模式,以期為該地區(qū)粳稻抗稻瘟病育種提供種質(zhì)資源。

    1 材料與方法

    1.1 供試材料及菌株

    輪回親本07GY31由江蘇省農(nóng)業(yè)科學(xué)院糧食作物研究所育成,該品系農(nóng)藝性狀優(yōu)良,但高感稻瘟病。NIL、NIL、NIL、NIL、NIL、NIL、NIL、NIL和NIL9個廣譜抗性基因近等基因系均以07GY31為輪回親本構(gòu)建而成[36]。

    2013—2019年分別從江蘇、安徽、浙江、江西、湖北等長江下游粳稻種植區(qū)采集并分離109個稻瘟病菌菌株應(yīng)用于本研究。稻瘟病菌的單孢分離、菌株培養(yǎng)及菌液制備按照Puri等[37]報道的方法進行。

    1.2 分子標記檢測

    1.2.1 基于分子標記的前景選擇 在雙基因聚合系構(gòu)建過程中,采用簡易TPS快速提取法提取水稻基因組DNA[38]。利用表1所示引物進行基因擴增,以便對各目標基因進行前景選擇。20 μL PCR反應(yīng)體系含MgCl2(25 mmol·L-1)2.0 μL,10×PCR buffer 2.0 μL,引物(上下游引物,10 μmol·L-1)1.5 μL,dNTP(10 mmol·L-1)0.4 μL,DNA聚合酶(5 U·μL-1)0.2 μL,模板DNA 2.0 μL;ddH2O 11.9 μL。PCR反應(yīng)程序為94℃預(yù)變性5 min;94℃變性45 s;55—58℃(視不同引物而定)退火45 s;72℃延伸1 min,進行35個循環(huán),最后72℃充分延伸10 min。根據(jù)擴增產(chǎn)物的大小,分別采用8%的變性聚丙烯酰胺凝膠及4%瓊脂糖凝膠電泳對目標擴增產(chǎn)物進行電泳分離。

    1.2.2 基于genotyping by sequencing(GBS)的背景分析 利用DNAsecure Plant kit試劑盒(Qiagen,USA)提取所有試驗材料的基因組DNA。利用BioPhotometer plus核酸測定儀(Eppendorf,Germany)檢測核酸濃度及質(zhì)量。高質(zhì)量的DNA用于GBS文庫構(gòu)建和測序,首先應(yīng)用限制性內(nèi)切酶H I和I對基因組DNA進行酶切,酶切后的片段參考Poland等[39]方法進行建庫并用第二代測序儀Illumina HiSeq 2000 Sequencer(Illumina,USA)進行雙末端(pair end)測序。測序深度為~10×coverage,去除原始測序序列中的接頭序列和低質(zhì)量的短讀序列,剩下的高質(zhì)量序列采用BWA(Burrows-Wheeler Aligner)程序?qū)⑵淦ヅ涞饺毡厩鐓⒖蓟蚪M上[40]。利用GATK(V4.1)軟件[41]提取SNP,每個聚合系以輪回親本對照,利用Tassel(V5.0)去除聚合系與輪回親本間相同的SNP,不同SNP在染色體上的密度做圖參照R語言包(Cmplot V3.6.0),隨機選擇3萬個SNP標記用于背景恢復(fù)率的計算,背景恢復(fù)率=(相同標記數(shù)/總標記數(shù))×100%。

    表1 用于檢測目的基因的分子標記詳細信息

    1.3 抗性鑒定

    1.3.1 人工接種鑒定 2019年5月在江蘇里下河地區(qū)農(nóng)業(yè)科學(xué)研究所萬?;夭捎?09個稻瘟病菌菌株對9個單基因系、18個雙基因聚合系進行苗瘟抗性鑒定。供試水稻材料播種于塑料盆中,加入營養(yǎng)土于人工氣候室中育苗(27—30℃,正常光照)。所有供試水稻材料都隨機設(shè)置3個重復(fù),并以輪回親本07GY31及普感品種LTH作為陰性對照。待秧苗長至3—4葉期時移至玻璃接種箱內(nèi)。將孢子懸浮液配制成濃度為5×104個/ml,并加入0.02% Tween-20作為表面活性劑。用連接空氣壓縮機的噴槍在接種箱內(nèi)噴霧接種,每個秧盤內(nèi)接種40—50 ml的孢子懸浮液。26℃黑暗保濕24 h,然后轉(zhuǎn)移至可控溫室內(nèi),在溫度為26—28℃、相對濕度>95%的溫濕環(huán)境下生長。7—10 d后調(diào)查、記載各家系植株葉片上的病斑數(shù)目及反應(yīng)類型,病情調(diào)查按Mackill等[42]分級標準進行。

    2019年正季在自然條件下,采用人工注射接種法對以上材料進行穗瘟抗性鑒定。參考Wu等[43]的選擇標準,從109個稻瘟病菌菌株中挑選出28個代表性菌株用于穗瘟接種鑒定。每個試驗材料種植成120株的小區(qū),每區(qū)10行,每行12株,株行距13.3 cm×25 cm,采用完全隨機區(qū)組設(shè)計(RCBD),每個試驗材料重復(fù)種植兩次,常規(guī)水肥管理。8月初選擇孕穗期劍葉葉枕與倒二葉距離為3—4 cm左右的穗苞進行接種,每個菌株接種每個株系10穗,每穗注射1 mL,孢子接種濃度為5×104個/ml,并對接種后的稻穗作好標記。接種21 d后調(diào)查發(fā)病情況,穗瘟鑒定評價參照Wu等[36]標準??剐灶l率(resistance frequency,RF)=(接種后表現(xiàn)抗的菌株數(shù)/接種總菌株數(shù))×100%。

    1.3.2 病圃的多點自然誘發(fā)鑒定 2019年正季將以上單基因系及雙基因聚合系分別種植于江蘇金壇、安徽廬江和浙江長興的常年稻瘟病重發(fā)區(qū)。每個試驗材料栽10行,每行10株,株行距13.3 cm×25 cm,采用完全隨機區(qū)組設(shè)計(RCBD),每個試驗材料重復(fù)種植3次。為了使每個鑒定材料都能與誘發(fā)品種接觸,每隔兩個鑒定小區(qū)的材料播一行感病誘發(fā)品種。病圃水層保持10 cm,進行田間自然誘發(fā)。試驗田不使用任何殺菌劑,其他按常規(guī)肥水管理。齊穗30 d后調(diào)查每個鑒定材料的穗發(fā)病情況。調(diào)查、評價方法參照Wu等[29],健康稻穗率(health panicle proportion,HPP)=(小區(qū)未感染稻瘟病穗數(shù)/小區(qū)總穗數(shù))×100%。

    1.4 基本農(nóng)藝性狀調(diào)查

    除用于接種鑒定外,各株系材料同時種植一份用于農(nóng)藝性狀調(diào)查。5月18日播種,6月13日進行移栽,單本栽插,每區(qū)10行,每行12株,株行距13.3 cm×25 cm,常規(guī)水肥管理。適期取每株系中間第2—9行之間的5株材料進行抽穗期、株高、單株有效穗數(shù)、每穗總粒數(shù)、結(jié)實率、千粒重、單株產(chǎn)量等基本農(nóng)藝性狀的調(diào)查。

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

    利用Microsoft Excel 2016進行數(shù)據(jù)整理和表格繪制,多重比較采用IBM SPSS Statistics 22進行數(shù)據(jù)分析,雙基因聚合系抗性效應(yīng)的構(gòu)成因子分析的?;鶊D及雙基因聚合系在不同病圃抗性表現(xiàn)的三維散點圖均使用Python語言包(v 3.7.7)進行作圖。

    2 結(jié)果

    2.1 構(gòu)建的雙基因聚合系具有較高的背景恢復(fù)率

    雙基因聚合系構(gòu)建方法參考Wu等[29]的方法,采用不完全雙列雜交法(NCII)配制18個雜交F1組合,針對每個組合,開展目標基因的特異性標記檢測,篩選目標雙基因均為雜合帶型的F1單株進行自交,獲得1 000個左右F2單株的群體。結(jié)合標記輔助選擇和農(nóng)藝性狀篩選,在不同F(xiàn)2群體中分別篩選出4—22個目標雙基因純合且基本農(nóng)藝性狀類似輪回親本07GY31的雙基因聚合單株,將這些單株加代到F3種植成株系。進一步選擇2—4個農(nóng)藝性狀與07GY31基本一致的F3株系,通過GBS分析進行背景回復(fù)率檢測。選擇背景恢復(fù)率最高的一個株系用于后續(xù)的抗性鑒定評價。GBS分析結(jié)果表明不同雙基因聚合系均具有較高的背景恢復(fù)率,普遍在97.00%以上,分布于97.08%(PPL)—99.08%(PPL)(圖1)。表明除了目標基因區(qū)域不同外,所有雙基因聚合系的遺傳背景幾乎完全與受體親本07GY31一致。

    2.2 雙基因聚合系抗性水平普遍優(yōu)于單基因系

    利用采集到的稻瘟病菌菌株對9個單基因系、18個雙基因聚合系及輪回親本分別在苗期及抽穗期進行人工接種鑒定,結(jié)果表明不同抗性基因聚合后,其苗瘟及穗瘟的抗性水平顯著高于單基因系(圖2-A、2-B)。就苗瘟抗性而言,不同基因聚合后產(chǎn)生了不同的聚合效應(yīng),導(dǎo)致不同聚合系之間的抗性水平存在顯著差異(圖2-C)。根據(jù)抗性水平的不同,可將不同聚合系劃分為兩個抗性等級,第一級為PPL、PPL、PPL、PPL、PPL、PPL、PPL、PPL、PPL9個雙基因聚合系,抗性頻率在88.99%—99.08%;第二級則為以PPL為代表的另外9個雙基因聚合系,抗性頻率分布在28.44%—73.39%。針對高抗性水平基因而言,盡管NIL的抗性頻率高達94.50%,但是與和聚合后其抗性頻率均提高了4.58%。對于低抗性水平基因而言,其單基因系的抗性頻率為27.52%,與、和聚合后其抗性頻率則分別提高了20.19%、13.76%和0.92%。

    針對穗瘟而言,不同基因聚合后同樣產(chǎn)生不同聚合效應(yīng),導(dǎo)致不同聚合系間抗性水平存在顯著差異(圖2-D)??傮w而言,穗瘟的抗性頻率普遍要低于苗瘟。根據(jù)抗性水平的不同,同樣將聚合系抗性水平劃分為兩個等級,第一級為PPL、PPL、PPL、PPL、PPL、PPL、PPL7個雙基因聚合系,抗性頻率在64.29%—85.71%。第二級則為以PPL為代表的11個雙基因聚合系,抗性頻率分布在21.43%—50.00%。盡管NIL的穗瘟抗性頻率為82.14%,但是與和聚合后其抗性頻率均提高了3.57%。值得注意的是,單基因系抗性頻率僅為35.71%,與(RF=28.57%)聚合后,其雙基因聚合系PPL的抗性頻率達到64.29%;而與(RF=39.29%)和(RF=35.71%)聚合后,其雙基因聚合系PPL和PPL的抗性頻率同樣達到46.43%,相對NIL的抗性水平得到顯著提高。由此可見,對于部分抗性效應(yīng)不好的單基因,選擇合適的基因組合方式,可達到顯著提高抗性水平的目的。

    2.3 互補效應(yīng)為提升雙基因聚合系抗性水平的關(guān)鍵因子

    為進一步分析聚合系抗性提高的原因,將聚合系對稻瘟病菌菌株的抗性效應(yīng)歸為重疊效應(yīng)、互補效應(yīng)和互作效應(yīng)3種類型。分析各效應(yīng)與苗瘟及穗瘟的抗性關(guān)系,發(fā)現(xiàn)重疊效應(yīng)和互補效應(yīng)與聚合系的抗性顯著相關(guān),其中互補效應(yīng)及其能否有效表達是影響聚合系抗性水平的關(guān)鍵因子。在苗瘟抗性中,具有較高抗性頻率的第一級的9個聚合系均有較高的互補效應(yīng),其平均互補效應(yīng)為50.87%,并且46.48%的互補效應(yīng)能有效表達(即聚合的兩個目的基因一方為抗病表型,可以互補另一方的感病表型,導(dǎo)致聚合系為抗病表型),提高了聚合系的抗性效應(yīng),只有4.39%的互補效應(yīng)不能有效表達(即聚合的兩個目的基因一方為抗病表型不能互補另一方的感病表型,導(dǎo)致聚合系為感病表型),降低了聚合系的抗性效應(yīng)。而在以PPL為代表的抗性頻率較低的第二級的另外9個聚合系中,其平均互補效應(yīng)為43.94%,且只有27.32%的互補效應(yīng)得到有效表達,提高了聚合系的抗性效應(yīng),而有16.62%的互補效應(yīng)不能有效表達,導(dǎo)致聚合系出現(xiàn)感病表型(圖3-A)。同樣地,在穗瘟抗性反應(yīng)中,具有較高抗性頻率的第一級的7個聚合系均有較高的互補效應(yīng),其平均互補效應(yīng)為50.00%,且有39.80%互補效應(yīng)提高了聚合系的抗性效應(yīng),如穗瘟抗性表現(xiàn)最好的PPL、PPL和 PPL的互補效應(yīng)分別為46.43%、57.14%和60.71%。而以PPL為代表的抗性頻率較低的第二級的另外11個聚合系中,其平均互補效應(yīng)為43.08%,且僅有22.64%的互補效應(yīng)得以有效表達,提高了聚合系的抗性水平,而有20.44%的互補效應(yīng)不能表達,從而降低了聚合系的抗性效應(yīng)??剐孕?yīng)最低的PPL和PPL其互補效應(yīng)均僅為32.14%,且均只有14.28%的互補效應(yīng)提高了聚合系的抗性效應(yīng)(圖3-B)。此外,所有基因組合中都存在一定的正向互作效應(yīng)(即聚合的兩個目的基因雙方都為感病表型,聚合后表現(xiàn)為抗病表型),導(dǎo)致基因聚合后抗性效應(yīng)得到進一步的提高,如單基因系穗瘟抗性頻率僅為35.71%,與聚合后,分別產(chǎn)生了10.72%的重疊效應(yīng),25.00%的互補效應(yīng)及28.57%的互作效應(yīng),導(dǎo)致其雙基因聚合系PPL的穗瘟抗性頻率達到了64.29%。因此,抗病基因的組合方式直接決定了聚合系的抗性水平,選擇互補效應(yīng)高且能有效表達的抗病基因組合模式對培育廣譜抗性品種具有重要作用。

    圓圈代表R基因所在的基因座位置The circle represents the locus of the R gene

    A、B:雙基因聚合系、單基因系和輪回親本間苗瘟和穗瘟的綜合比較分析Comprehensive comparative analysis on seedling blast and panicle blast RF of PPLs, NILs and the recurrent parent;C、D:雙基因聚合系、單基因系和輪回親本的苗瘟、穗瘟抗性表現(xiàn)Resistance performances of PPLs, NILs and the recurrent parent for seedling blast and panicle blast resistance

    2.4 含有Pigm的基因組合在不同病圃間抗性表現(xiàn)最為穩(wěn)定

    為進一步明確不同聚合系在自然種植、病原菌強脅迫條件下的發(fā)病情況,各聚合材料及輪回親本07GY31分別種植于江蘇金壇、安徽廬江和浙江長興3個稻瘟病重發(fā)區(qū),進行田間自然誘發(fā)鑒定。結(jié)果顯示,輪回親本07GY31在這3個地區(qū)均屬于高感類型(表2),表明當(dāng)?shù)氐疚敛“l(fā)病嚴重,是理想的稻瘟病自然誘發(fā)鑒定病圃。同時對比人工接種鑒定結(jié)果與自然鑒定結(jié)果,發(fā)現(xiàn)人工鑒定的各個材料的穗瘟抗性表現(xiàn)與在金壇、廬江和長興的自然鑒定抗性表現(xiàn)基本一致,其決定系數(shù)(2)分別為0.7828、0.6788和0.8236(圖4-A)。此外,PPL、PPL和PPL3個雙基因聚合系在不同鑒定點抗性效應(yīng)最好也最為穩(wěn)定,小區(qū)健康稻穗率均在94.00%以上,與人工鑒定結(jié)果基本一致(圖4-B、表2)。其他雙基因聚合系的田間穗瘟抗性水平存在顯著差異,且在不同鑒定點間抗性水平表現(xiàn)明顯的特異性。如PPL在廬江的小區(qū)健康稻穗率為91.00%,表現(xiàn)為3級抗性水平;而在金壇和長興的小區(qū)健康稻穗率分別為77.00%和69.00%,表現(xiàn)為7級感病水平(圖4-B、表2)。表明基因組合的抗性表現(xiàn)存在地區(qū)特異性,對金壇和長興的稻瘟病菌群體表現(xiàn)一定程度的特異親和性。

    A:雙基因聚合系苗瘟抗性效應(yīng)分析桑基圖Analysis of resistance effect of PPLs against seedling blast with sankey diagram;B:雙基因聚合系穗瘟抗性效應(yīng)分析?;鶊DAnalysis of resistance effect of PPLs against panicle blast with sankey diagram。OE:重疊效應(yīng)overlapping effect;CE:互補效應(yīng)complementary effect;IE:互作效應(yīng)interaction effect

    2.5 雙基因聚合系的基本農(nóng)藝性狀表現(xiàn)

    考察了雙基因聚合系和輪回親本07GY31的基本農(nóng)藝性狀(表2)。結(jié)果表明,大多數(shù)雙基因聚合系的農(nóng)藝表現(xiàn),如株高、抽穗期以及單株穗數(shù)、每穗粒數(shù)、結(jié)實率、千粒重及單株產(chǎn)量等均與輪回親本07GY31無明顯差異,而攜帶有的雙基因聚合系與輪回親本相比則表現(xiàn)為生育期提前,從產(chǎn)量構(gòu)成性狀看,盡管其千粒重與輪回親本無差異,但每穗總粒數(shù)和結(jié)實率顯著降低,最終導(dǎo)致單株產(chǎn)量顯著降低。

    3 討論

    3.1 含有Pigm的基因組合為有效的抗性基因組合模式

    稻瘟病菌生理小種的組成及頻率在不同地區(qū)、不同年份間存在很大的差異和變化。研究結(jié)果表明,我國稻瘟病菌生理小種的地理分布存在明顯的地域差異,且秈、粳稻上病菌的致病性也有明顯分化,秈型小種主要分布于南方秈稻區(qū),而粳型小種則以北方粳稻區(qū)和長江下游秈、粳稻混栽區(qū)居多[30,44]。在大田種植過程中水稻與稻瘟病菌之間會發(fā)生協(xié)同進化,因此評價一個抗性基因的效應(yīng)不能僅用少數(shù)幾個菌株進行鑒定;同樣在抗性品種推廣應(yīng)用之前,也需用采用當(dāng)?shù)厣鷳B(tài)區(qū)有代表性的菌株加以鑒定。在本研究中,利用長江下游109個粳型代表性菌株對構(gòu)建的18個基因組合進行抗性評價,并分別在江蘇金壇、安徽廬江和浙江長興進行自然誘發(fā)鑒定,結(jié)果發(fā)現(xiàn)雙基因聚合系PPL、PPL、PPL在苗瘟和穗瘟人工接種中表現(xiàn)出廣譜抗性,其苗、穗瘟抗性頻率分別在85.71%、85.71%和82.14%以上,并在不同病圃也表現(xiàn)出穩(wěn)定的廣譜抗性,小區(qū)健康稻穗率分布在94.00%—100.00%。前期筆者采集了158個來自南方稻區(qū)的秈型菌株對秈稻背景的不同基因組合進行抗性鑒定,/、/和/3個基因組合同樣表現(xiàn)出廣譜抗性,其苗瘟和穗瘟抗性頻率分別在92.41%和86.67%以上,并且在3個南方秈稻病圃里面表現(xiàn)出穩(wěn)定而廣譜的抗性水平[29]。這些結(jié)果表明基因組合/、/和/不僅對長江下游粳型稻瘟病菌菌群表現(xiàn)出廣譜抗性,同樣對南方秈型菌群表現(xiàn)廣譜抗性,是具有重要應(yīng)用價值的廣譜稻瘟病抗性基因組合模式。

    A:3個病圃的不同雙基因聚合系的小區(qū)健康穗率與人工接種下穗瘟抗性頻率之間的關(guān)系Relationship between HPP of PPLs at three disease nurseries and the RF of PPLs in artificial inoculation evaluation;B:雙基因聚合系在不同病圃的穗瘟抗性表現(xiàn)Resistance performance of PPLs against panicle blast among different disease nurseries。JT:金壇Jintan;LJ:廬江Lujiang;CX:長興Changxing

    3.2 不同基因組合方式影響聚合系的抗性效應(yīng)

    在優(yōu)良受體親本中聚合非等位抗性基因可有效提高目標材料對稻瘟病廣譜抗性,但是對所要聚合的目標基因的選擇尤為重要。因為兩個非等位抗性基因聚合后,并不簡單地表現(xiàn)為單個抗病基因的抗譜之間的簡單累加,而是抗性基因之間表現(xiàn)為極顯著的基因相互作用,產(chǎn)生多種聚合效應(yīng)[45-47],在本研究中將其歸類為重疊效應(yīng)、互補效應(yīng)和互作效應(yīng)。基因聚合效應(yīng)的構(gòu)成因子極為復(fù)雜,不同基因組合方式會產(chǎn)生不同的聚合效應(yīng);同一個基因與不同基因聚合也會產(chǎn)生不同的聚合效應(yīng)。以含有的雙基因聚合系穗瘟抗性效應(yīng)為例,單基因系穗瘟抗性效應(yīng)僅為35.71%,與(RF=39.29%)聚合后分別產(chǎn)生了14.29%的重疊效應(yīng)、28.57%的互補效應(yīng)和3.57%的互作效應(yīng),導(dǎo)致雙基因聚合系PPL的抗性效應(yīng)為46.43%;與(RF=35.71%)聚合后分別產(chǎn)生了10.71%的重疊效應(yīng)、21.43%的互補效應(yīng)和14.29%的互作效應(yīng),導(dǎo)致雙基因聚合系PPL的抗性效應(yīng)為46.43%;而與(RF=28.57%)聚合后則分別產(chǎn)生了7.14%的重疊效應(yīng)、28.57%的互補效應(yīng)和28.57%的互作效應(yīng),導(dǎo)致雙基因聚合系PPL的抗性效應(yīng)為64.28%(圖2-D)。兩個非等位基因聚合后相互作用的分子機制十分復(fù)雜,目前還不清楚,然而聚合能夠識別不同菌株的抗性基因仍是提高目標材料廣譜抗性的有效途徑之一,但是在聚合過程中必須對單個目標基因的抗性水平及聚合后的抗性效應(yīng)進行比較分析,選擇互補效應(yīng)高且能產(chǎn)生較強互作效應(yīng)的抗病基因組合模式,才能有效提高目標材料的廣譜抗性。

    表2 18個雙基因聚合系及輪回親本07GY31的基本農(nóng)藝性狀表現(xiàn)

    同列數(shù)據(jù)后不同小寫字母表示差異顯著(<0.05)Different lowercase letters after the data in the same column indicate significant difference (<0.05)

    3.3 不同受體遺傳背景影響基因組合的抗性效應(yīng)

    受體的遺傳背景影響抗性基因/組合的抗性效應(yīng)[48]。研究表明在3種不同遺傳背景下對來自日本的稻瘟病菌菌株表現(xiàn)為部分抗性,而在TP309遺傳背景下則對來自中國的稻瘟病菌菌株表現(xiàn)完全抗性[22,49]??剐曰蚝驮诙i、粳基因組中均有分布,但僅在秈稻背景下有較好的抗性效應(yīng),而僅在粳稻背景下檢測到對抗性表型有較大的貢獻[34]。同樣在水稻白葉枯病抗性中也發(fā)現(xiàn)受體品種的遺傳背景可能會影響抗性基因的抗性表型,如基因組合在5個不同水稻受體品種中的抗性表現(xiàn)差異顯著,表明在某些水稻品種背景中可能存在抗性基因組合的抑制因子[50]。在本研究中,筆者發(fā)現(xiàn)基因組合在粳稻07GY31背景下的苗瘟及穗瘟抗性頻率分別為69.72%和46.43%,而在秈稻揚稻6號背景下的苗瘟及穗瘟抗性頻率則分別達到了95.57%和83.33%[29];同樣地,基因組合在粳稻07GY31背景下的苗瘟及穗瘟抗性頻率分別為37.61%和32.14%,而在秈稻揚稻6號背景下的苗瘟及穗瘟抗性頻率則分別達到了91.14%和70.00%[29];有趣的是基因組合/在粳稻背景下的苗瘟及穗瘟抗性頻率分別為91.74%和75.00%,高于秈稻揚稻6號背景下的苗瘟及穗瘟抗性頻率(86.71%和66.67%)[29]。這些結(jié)果表明不同遺傳背景中可能存在不同的抗性基因組合促進/抑制因子,這就暗示在不同生態(tài)區(qū),不同遺傳背景下選擇合適的抗性基因組合,做到抗性基因組合的合理布局,才能有針對性地培育廣譜持久抗性品種。

    4 結(jié)論

    抗性基因的組合方式直接影響了聚合系的抗性水平,而互補效應(yīng)高且能有效表達是粳型雙基因聚合系抗性效應(yīng)提高的關(guān)鍵因子。雙基因聚合系PPL、PPL和PPL在苗瘟和穗瘟人工接種及病圃自然誘發(fā)鑒定中均表現(xiàn)出穩(wěn)定的抗性,為長江下游粳稻的廣譜抗性基因組合模式。研究結(jié)果可為長江下游廣譜稻瘟病抗性粳稻品種選育提供種質(zhì)資源和理論支撐。

    [1] DEAN R, VAN KAN J A, PRETORIUS Z A, HAMMOND-KOSACK K E, DI PIETRO A, SPANU P D, RUDD J J, DICKMAN M, KAHMANN R, ELLIS J, FOSTER G D. The top 10 fungal pathogens in molecular plant pathology. Molecular Plant Pathology,2012, 13(4): 414-430.

    [2] KHUSH G S, JENE K. Current status and future prospects for research on blast resistance in rice (L)//Advances in Genetics, Genomics and Control of Rice Blast DiseaseSpringer, 2009: 1-10.

    [3] GOUDA P K, SAIKUMAR S, VARMA C M, NAGESH K, THIPPESWARMY S, SHENOY V, RAMESHA M S, SHASHIDHAR H E. Marker-assisted breeding ofandgenes imparting resistance to rice blast in PRR78, restorer line of Pusa RH-10 Basmati rice hybrid. Plant Breeding,2013, 132(1): 61-69.

    [4] ELLUR R K, KHANNA A, YADAV A, PATHANIA S, RAJASHEKARA H, SINGH V K, GOPALA KRISHNAN S, BHOWMICK P K, NAGARAJAN M, VINOD K K, PRAKASH G, MONDAL K K, SINGH N K, VINOD PRABHU K, SINGH A K. Improvement of Basmati rice varieties for resistance to blast and bacterial blight diseases using marker assisted backcross breeding. Plant Science,2016, 242: 330-341.

    [5] NIKS R E, QI X, MARCEL T C. Quantitative resistance to biotrophic filamentous plant pathogens: concepts, misconceptions, and mechanisms. Annual Review of Phytopathology,2015, 53: 445-470.

    [6] WANG G L, MACKILL D J, BONMAN J M, MCCOUCH S R, CHAMPOUX M C, NELSON R J. RFLP mapping of genes conferring complete and partial resistance to blast in a durably resistant rice cultivar. Genetics,1994, 136(4): 1421-1434.

    [7] LUKASIK E, TAKKEN F L. STANDing strong, resistance proteins instigators of plant defence. Current Opinion in Plant Biology,2009, 12(4): 427-436.

    [8] RAY S, SINGH P K, GUPTA D K, MAHATO A K, SARKAR C, RATHOUR R, SINGH N K, SHARMA T R. Analysis ofgenome reveals a fungal effector, which is able to induce resistance response in transgenic rice line containing resistance gene,. Frontiers in Plant Science,2016, 7: 1140.

    [9] Jia Y, McAdams S A, Bryan G T, Hershey H P, Valent B. Direct interaction of resistance gene and avirulence gene products confers rice blast resistance. The EMBO Journal,2000, 19(15): 4004-4014.

    [10] Park C H, Shirsekar G, Bellizzi M, Chen S, Songkumarn P, Xie X, Shi X, Ning Y, Zhou B, Suttiviriya P, Wang M, Umemura K, Wang G L. The E3 ligase APIP10 connects the effector AvrPiz-t to the NLR receptor Piz-t in rice. PLoS Pathogens,2016, 12(3): e1005529.

    [11] Fujisaki K, Abe Y, Ito A, Saitoh H, Yoshida K, Kanzaki H, Kanzaki E, Utsushi H, Yamashita T, Kamoun S, Terauchi R. Rice Exo70 interacts with a fungal effector, AVR-Pii, and is required for AVR-Pii-triggered immunity. The Plant Journal,2015, 83(5): 875-887.

    [12] Ashkani S, Rafii M, Shabanimofrad M, Ghasemzadeh A, Ravanfar S A, Latif M. Molecular progress on the mapping and cloning of functional genes for blast disease in rice (L.): current status and future considerations. Critical reviews in biotechnology,2016, 36(2): 353-367.

    [13] Sharma T, Rai A, Gupta S, Vijayan J, Devanna B, Ray S. Rice blast management through host-plant resistance: retrospect and prospects. Agricultural Research,2012, 1(1): 37-52.

    [14] Xie Z, Yan B X, Shou J Y, Tang J, Wang X, Zhai K R, Liu J Y, Li Q, Luo M Z, Deng Y W, He Z H.A nucleotide-binding site-leucine-rich repeat receptor pair confers broad-spectrum disease resistance through physical association in rice. Philosophical Transactions of the Royal Society B, 2019, 374(1767): 20180308.

    [15] WANG B H, Ebbole D J, WANG Z H. The arms race betweenand rice: diversity and interaction ofandgenes. Journal of Integrative Agriculture, 2017, 16(12): 2746-2760.

    [16] Zhao H j, Wang X y, Jia Y l, Minkenberg B, Wheatley M, Fan J b, Jia M H, Famoso A, Edwards J D, Wamishe Y, Valent B, Wang G L, Yang Y n. The rice blast resistance geneencodes an atypical protein required for broad-spectrum disease resistance. Nature communications, 2018, 9(1): 2039.

    [17] Fukuoka S, Yamamoto S I, Mizobuchi R, Yamanouchi U, Ono K, Kitazawa N, Yasuda N, Fujita Y, Nguyen T T, Koizumi S, SUGIMOTO K, MATSUMOTO T, YANO M. Multiple functional polymorphisms in a single disease resistance gene in rice enhance durable resistance to blast. Scientific Reports,2014, 4: 4550.

    [18] Fukuoka S, Saka N, Koga H, Ono K, Shimizu T, Ebana K, Hayashi N, Takahashi A, Hirochika H, Okuno K, Yano M. Loss of function of a proline-containing protein confers durable disease resistance in rice. Science,2009, 325(5943): 998-1001.

    [19] Hayashi N, Inoue H, Kato T, Funao T, Shirota M, Shimizu T, Kanamori H, Yamane H, Hayano-Saito Y, Matsumoto T, Yano M, Takatsuji H. Durable panicle blast-resistance geneencodes an atypical CC-NBS-LRR protein and was generated by acquiring a promoter through local genome duplication. The Plant Journal,2010, 64(3): 498-510.

    [20] Xu X, Hayashi N, Wang C T, Fukuoka S, Kawasaki S, Takatsuji H, Jiang C J. Rice blast resistance gene, a member of a resistance gene cluster on chromosome 4, encodes a nucleotide-binding site and leucine-rich repeat protein. Molecular Breeding,2014, 34(2): 691-700.

    [21] Inukai T, Nagashima S, Kato M.is a race-specific partial-resistance allele at the Pid3 blast resistance locus in rice. Theoretical and Applied Genetics,2019, 132(2): 395-404.

    [22] Variar M, Cruz CV, Carrillo M, Bhatt J, Sangar R. Rice blast in India and strategies to develop durably resistant cultivars//Advances in Genetics, Genomics and Control of Rice Blast DiseaseSpringer, 2009: 359-373.

    [23] 宛柏杰, 劉凱, 趙紹路, 朱靜雯, 劉艷艷, 張桂云, 朱國永, 王愛民, 唐紅生, 孫明法, 嚴國紅. 水稻抗稻瘟病基因、、和在骨干親本中的分布以及對穗頸瘟抗性的作用. 西南農(nóng)業(yè)學(xué)報, 2020, 33(1): 1-6.

    WAN B j, LIU K, ZHAO S l, ZHU J w, LIU Y y, ZHANG G y, ZHU g y, WANG A m, TANG H s, SUN M f, YAN g h. Distribution of rice blast resistance genes,,and

    [24] 朱勇良, 范方軍, 謝裕林, 伍應(yīng)保, 喬中英, 張建棟. 江蘇省遲熟中粳新品系稻瘟病抗病基因檢測與抗性評價. 江蘇農(nóng)業(yè)科學(xué), 2018, 46(19): 106-109.

    ZHU Y L, FAN F J, XIE Y L, WU Y B, QIAO Z Y, ZHANG J D. Detection and evaluation of blast resistance genes in late mature medium japonica lines in Jiangsu Province. Jiangsu Agricultural Sciences, 2018, 46(19): 106-109. (in Chinese)

    [25] Jiang H c, Feng Y t, Bao L, Li X, Gao G j, Zhang Q l, Xiao J h, Xu C g, He Y q. Improving blast resistance of Jin 23B and its hybrid rice by marker-assisted gene pyramiding. Molecular Breeding,2012, 30(4): 1679-1688.

    [26] Xiao W m, Yang Q y, Huang M, Guo T, Liu Y z, Wang J f, Yang G l, Zhou J y, Yang J y, Zhu X y, Chen Z q, Wang H. Improvement of rice blast resistance by developing monogenic lines, two-gene pyramids and three-gene pyramid through MAS. Rice,2019, 12(1): 78.

    [27] Hittalmani S, Parco A, Mew T, Zeigler R, Huang N. Fine mapping and DNA marker-assisted pyramiding of the three major genes for blast resistance in rice. Theoretical and Applied Genetics,2000, 100(7): 1121-1128.

    [28] Xiao N, Wu Y y, Pan C h, Yu L, Chen Y, Liu G q, Li Y H, Zhang X x, Wang Z p, Dai Z y, Liang C Z, Li A h. Improving of rice blast resistances inby pyramiding major R genes. Frontiers in Plant Science,2017, 7: 1918.

    [29] Wu Y y, Xiao N, Chen Y, Yu L, Pan C h, Li Y h, Zhang X x, Huang N s, Ji H j, Dai Z y, Chen X j, Li A h. Comprehensive evaluation of resistance effects of pyramiding lines with different broad-spectrum resistance genes againstin rice (L.). Rice,2019, 12(1): 11.

    [30] Zhong Z H, Chen M L, Lin L Y, Han Y J, Bao J D, Tang W, Lin L L, Lin Y H, Somai R, Lu L,. Population genomic analysis of the rice blast fungus reveals specific events associated with expansion of three main clades. The ISME Journal,2018, 12(8): 1867-1878.

    [31] Xiao N, Wu Y Y, Wang Z P, Li Y H, Pan C H, Zhang X X, Yu L, Liu G Q, Zhou C H, Ji H J, Huang N S, Jiang M, Dai Z Y, Li A H. Improvement of seedling and panicle blast resistance inrice varieties followingintrogression. Molecular Breeding, 2018, 38: 142.

    [32] 陳波, 周年兵, 郭保衛(wèi), 黃大山, 陳忠平, 花勁, 霍中洋, 張洪程. 南方稻區(qū)“秈改粳”研究進展. 揚州大學(xué)學(xué)報 (農(nóng)業(yè)與生命科學(xué)版), 2017, 38(1): 67-72, 88.

    CHEN B, ZHOU N B, GUO B W, HUANG D S, CHEN Z P, HUA J, HUO Z Y, ZHANG H C. Progress of “rice torice” in southern China. Journal of Yangzhou University (Agricultural and Life Science Edition), 2017, 38(1): 67-72, 88. (in Chinese)

    [33] 殷敏, 劉少文, 褚光, 徐春梅, 王丹英, 章秀福, 陳松. 長江下游稻區(qū)不同類型雙季晚粳稻產(chǎn)量與生育特性差異. 中國農(nóng)業(yè)科學(xué), 2020, 53(5): 890-903.

    YIN M, LIU S W, CHU G, XU C M, WANG D Y, ZHANG X F, CHEN S. Differences in yield and growth traits of differentvarieties in the double cropping late season in the lower reaches of the Yangtze River.Scientia Agricultura Sinica, 2020, 53(5): 890-903. (in Chinese)

    [34] Wu Y y, Xiao N, Yu L, Pan C h, Li Y h, Zhang X x, Liu G q, Dai Z y, Pan X b, Li A h. Combination patterns of major R genes determine the level of resistance to thein rice (L.). PLoS One,2015, 10(6): e0126130.

    [35] 王小秋, 杜海波, 陳夕軍, 李明友, 王嘉楠, 許志文, 馮志明, 陳宗祥, 左示敏. 江蘇近年育成粳稻新品種/系的稻瘟病抗性基因及穗頸瘟抗性分析. 中國水稻科學(xué), 2020, 34(5): 413-424.

    WANG X Q, DU H B, CHEN X J, LI M Y, WANG J N, XU Z W, FENG Z M, CHEN Z X, ZUO S M. Analysis of blast resistant genes and neck blast resistance ofrice varieties/lines recently developed in Jiangsu Province. Chinese Journal of Rice Science, 2020, 34(5): 413-424. (in Chinese)

    [36] Wu Y y, Chen Y, Pan C h, Xiao N, Yu L, Li Y h, Zhang X x, Pan X b, Chen X j, Liang C z, Dai Z y, Li A h. Development and evaluation of near-isogenic lines with different blast resistance alleles at thelocus inrice from the lower region of the Yangtze River, China. Plant Disease,2017, 101(7): 1283-1291.

    [37] Puri K D, Shrestha S M, Chhetri G B K, Joshi K D. Leaf and neck blast resistance reaction in tropical rice lines under green house condition. Euphytica,2009, 165(3): 523-532.

    [38] 盧揚江, 鄭康樂. 提取水稻DNA的一種簡易方法. 中國水稻科學(xué), 1992, 6(1): 47-48.

    Lu Y J, Zheng K L. A simple method for isolation of rice DNA. Chinese Journal of Rice Science, 1992, 6(1): 47-48. (in Chinese)

    [39] Poland J A, Brown P J, Sorrells M E, Jannink J L. Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach. PLoS One,2012, 7(2): e32253.

    [40] Kawahara Y, de la Bastide M, Hamilton J P, Kanamori H, McCombie W R, Ouyang S, Schwartz D C, Tanaka T, Wu J, Zhou S,. Improvement of theNipponbare reference genome using next generation sequence and optical map data. Rice, 2013, 6(1): 4.

    [41] Walker M A, Pedamallu C S, Ojesina A I, Bullman S, Sharpe T, Whelan C W, Meyerson M. GATK PathSeq: a customizable computational tool for the discovery and identification of microbial sequences in libraries from eukaryotic hosts. Bioinformatics, 2018, 34(24): 4287-4289.

    [42] Mackill D, Bonman J. Inheritance of blast resistance in near-isogenic lines of rice. Phytopathology,1992, 82(7): 746-749.

    [43] Wu Y y, Yu L, Pan C h, Dai Z y, Li Y h, Xiao N, Zhang X x, Ji H j, Huang N s, Zhao B h,. Development of near-isogenic lines with different alleles of Piz locus and analysis of their breeding effect under Yangdao 6 background. Molecular Breeding,2016, 36(2): 12.

    [44] 彭洪江, 張杰, 饒宗文, 彭士鐘, 吳先麗. 不同生態(tài)區(qū)的稻瘟病調(diào)查. 西南農(nóng)業(yè)學(xué)報, 1995, 8(1): 59-64.

    PENG H J, ZHANG J, RAO Z W, PENG S Z, WU X L. Investigation on rice blast in different ecological areas.Southwest China Journal of Agricultural Sciences, 1995, 8(1): 59-64. (in Chinese)

    [45] Divya B, Biswas A, Robin S, Rabindran R, Joel A J. Gene interactions and genetics of blast resistance and yield attributes in rice (L.). Journal of Genetics,2014, 93(2): 415-424.

    [46] Chaipanya C, Telebanco-Yanoria M J, Quime B, Longya A, Korinsak S, Korinsak S, Toojinda T, Vanavichit A, Jantasuriyarat C, Zhou B. Dissection of broad-spectrum resistance of the Thai rice variety Jao Hom Nin conferred by two resistance genes against rice blast. Rice,2017, 10(1): 18.

    [47] Chen X l, Jia Y l, Jia M H, Pinson S R M, Wang X y, Wu B M. Functional interactions between major rice blast resistance genes,and, and minor blast resistance quantitative trait loci. Phytopathology,2018, 108(9): 1095-1103.

    [48] Li W, Deng Y W, Ning Y S, He Z H, Wang G L. Exploiting broad-spectrum disease resistance in crops: from molecular dissection to breeding. Annual Review of Plant Biology,2020, 71: 575-603.

    [49] Xu X, Lv Q M, Shang J J, Pang Z Q, Zhou Z Z, Wang J, Jiang G H, Tao Y, Xu Q, Li X B, Zhao X F, Li S G, Xu J C, Zhu L H. Excavation oforthologs with differential resistance spectra toin rice resource. PLoS One,2014, 9(3): e93275.

    [50] Pradhan S K, Nayak D K, Mohanty S, Behera L, Barik S R, Pandit E, Lenka S, Anandan A. Pyramiding of three bacterial blight resistance genes for broad-spectrum resistance in deepwater rice variety, Jalmagna. Rice,2015, 8(1): 51.

    Construction and Analysis of Broad-spectrum Resistance Gene Combination Pattern forRice in Lower Region of the Yangtze River, China

    Wu YunYu1,2, Xiao Ning1,2,4, Yu Ling1,2, Cai Yue1,2, Pan CunHong1,3, Li YuHong1,2, Zhang XiaoXiang1,2, Huang NianSheng1,2, Ji HongJuan1,2, Dai ZhengYuan1,3, Li AiHong1,2,3

    1Lixiahe Institute of Agricultural Sciences of Jiangsu, Yangzhou 225007, Jiangsu;2Jiangsu Key Laboratory of Crop Genomics and Molecular Breeding, Yangzhou University, Yangzhou 225009, Jiangsu;3Jiangsu Co-Innovation Center for Modern Production Technology of Grain Crops, Yangzhou University, Yangzhou 225009, Jiangsu;4State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193

    【】Gene pyramiding is one of the most effective ways to achieve broad-spectrum resistance against. The objective of this study is to construct a set of polygene pyramiding lines (PPLs) under the background ofrice, to evaluate their resistance performances and analysis the components of their resistance effects usingstrains collected from lower region of the Yangtze River, China, thus providing broad-spectrum resistance gene combination pattern and germplasm resources forrice resistance breeding in lower region of the Yangtze River, China.【】Monogenic lines with multiple alleles of the Piz locus (,,,,and) with the background ofrice 07GY31 as the backbone, crossed with other broad-spectrum resistance gene,and, respectively using the incomplete NCII mating design. A total of 18 different PPLs were constructed using marker-assisted selection (MAS) and agronomic traits screening. Artificial inoculation assays at seedling and heading stage with 109 representativestrains collected from lower region of the Yangtze River, combined with natural induction identification under multiple field environments were conducted to evaluate the resistance performances of different PPLs, and to analyze the component factors of the resistance effects of the PPLs.【】Genotyping by sequencing (GBS) analysis shows that the constructed PPLs all have a high background recovery rate, which was ranging from 97.08% (PPL) to 99.08% (PPL), indicated that the genetic background of all PPLs was almost fully identical to that of the recurrent parent 07GY31. The seedling blast and panicle blast resistance levels of most PPLs were significantly higher than those of monogenic lines under artificial inoculation conditions, the PPLs with better resistance to seedling blast are PPL, PPL, PPL, PPL, PPL, PPL, PPL, PPLand PPL, respectively, and the PPLs with outstanding performance in panicle blast are PPL, PPL, PPL, PPL, PL, PPLand PPL, respectively. Different resistance gene combinations produced different effects after pyramiding. High complementary effect and which could be fully expressed is the key factor for the improvement of the resistance level of seedling blast and panicle blast of the PPLs. In addition, PPL, PPLand PPLdisplayed broad-spectrum resistance in artificial inoculation at seedling and heading stage, and showed stable broad-spectrum resistance under different disease nurseries. Besides, agronomic traits evaluation also showed PPLs with these three gene combinations were at par to the recurrent parent. Therefore,,andare broad-spectrum resistance gene combination patterns suitable forrice resistance breeding in lower region of the Yangtze River, China.【】The combination pattern of resistance genes affects the resistance level of the PPLs, and high complementary effect and which could be fully expressed is the key factor for the improvement of the resistance level of the PPLs inbackground. In addition, the development of PPLs and component factors analysis in this study provides valuable theoretical support and innovative germplasm resources for the precise breeding broad-spectrumvarieties in lower region of the Yangtze River, China.

    lower region of the Yangtze River;rice;; rice blast; gene pyramiding; effect analysis

    10.3864/j.issn.0578-1752.2021.09.006

    2020-07-13;

    2020-08-24

    國家自然科學(xué)基金(31801342,31971868)、江蘇省農(nóng)業(yè)科技自主創(chuàng)新基金(CX(18)2022)、江蘇省重點研發(fā)計劃(BE2019339)、江蘇省農(nóng)業(yè)重大新品種創(chuàng)制項目(PZCZ201702)、江蘇省作物基因組學(xué)和分子育種重點實驗室(BM2018003)、江蘇省自然科學(xué)基金(BK20181216)、中國農(nóng)業(yè)科學(xué)院植物病蟲害生物學(xué)國家重點實驗室開放基金(SKLOF201909)

    吳云雨,E-mail:wuyunyuyu@163.com。通信作者李愛宏,E-mail:yzlah@126.com

    (責(zé)任編輯 岳梅)

    猜你喜歡
    廣譜抗粳稻稻瘟病
    基于Meta-QTL和RNA-seq的整合分析挖掘水稻抗稻瘟病候選基因
    我國雙季早粳稻實現(xiàn)“零的突破”
    解讀粳稻品種南方種植播期
    挖掘這些基因,讓作物抗各種?。?020.3.27 科技日報)
    豫南粳稻機械化適應(yīng)品種篩選及利用研究
    不同藥劑防治苗稻瘟病、葉稻瘟病效果試驗研究
    播期播量對晚粳稻寧84農(nóng)藝性狀及產(chǎn)量的影響
    生物綠肥在稻瘟病防治中的應(yīng)用與示范
    幾種藥劑防治水稻稻瘟病穗瘟效果試驗
    關(guān)于頭孢唑肟合成工藝的研究
    久久精品影院6| 日韩欧美三级三区| 嫩草影院新地址| 亚洲av免费高清在线观看| 中出人妻视频一区二区| 久久久欧美国产精品| 免费观看在线日韩| 高清日韩中文字幕在线| 亚洲婷婷狠狠爱综合网| 亚洲国产精品国产精品| 亚洲色图av天堂| 精品人妻视频免费看| 亚洲av成人av| 国产成人a∨麻豆精品| 国产单亲对白刺激| 久久久久免费精品人妻一区二区| 国产精品久久久久久久电影| 亚洲精品日韩av片在线观看| 欧美一区二区国产精品久久精品| 亚洲aⅴ乱码一区二区在线播放| 乱系列少妇在线播放| 亚洲在线自拍视频| 国产一级毛片七仙女欲春2| 久久久久久伊人网av| 欧美在线一区亚洲| 国产精品永久免费网站| 熟女电影av网| 波多野结衣高清作品| 99riav亚洲国产免费| 国内久久婷婷六月综合欲色啪| 岛国毛片在线播放| 丰满的人妻完整版| 亚洲精品国产av成人精品| 欧洲精品卡2卡3卡4卡5卡区| 成人亚洲欧美一区二区av| 色噜噜av男人的天堂激情| 狠狠狠狠99中文字幕| 亚洲一区高清亚洲精品| 国产免费男女视频| 少妇人妻一区二区三区视频| 欧美zozozo另类| 免费看日本二区| 欧美日韩国产亚洲二区| 少妇被粗大猛烈的视频| 亚洲精品久久国产高清桃花| 亚洲欧美日韩高清在线视频| 亚洲精品乱码久久久久久按摩| 九色成人免费人妻av| 日本欧美国产在线视频| 天堂av国产一区二区熟女人妻| 悠悠久久av| 国产男人的电影天堂91| 最近的中文字幕免费完整| 99精品在免费线老司机午夜| 亚洲精品自拍成人| 校园人妻丝袜中文字幕| 日本三级黄在线观看| 久久久久久久久久久免费av| 久久99热这里只有精品18| 免费搜索国产男女视频| 国产 一区 欧美 日韩| 国产老妇女一区| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲欧美清纯卡通| 免费人成在线观看视频色| 国产精品一区二区性色av| 长腿黑丝高跟| 国产日本99.免费观看| 人人妻人人看人人澡| 大型黄色视频在线免费观看| 偷拍熟女少妇极品色| 成年免费大片在线观看| 久久久久久久久久久免费av| av又黄又爽大尺度在线免费看 | 久久久久性生活片| 久久鲁丝午夜福利片| 亚洲精品成人久久久久久| 一区二区三区四区激情视频 | 亚洲国产欧美在线一区| 最近的中文字幕免费完整| 在线国产一区二区在线| 人人妻人人澡欧美一区二区| 欧美激情在线99| 国产91av在线免费观看| 成年女人看的毛片在线观看| 成年女人永久免费观看视频| 91精品国产九色| 日韩,欧美,国产一区二区三区 | 欧美zozozo另类| 亚洲不卡免费看| 成人亚洲欧美一区二区av| 亚洲一级一片aⅴ在线观看| 毛片女人毛片| 亚洲欧美成人精品一区二区| 免费观看a级毛片全部| 狠狠狠狠99中文字幕| 亚洲美女搞黄在线观看| 亚洲一区二区三区色噜噜| 久久久久久国产a免费观看| 韩国av在线不卡| av在线亚洲专区| 精华霜和精华液先用哪个| 丝袜美腿在线中文| 亚洲av成人精品一区久久| 国产精品综合久久久久久久免费| 一夜夜www| 成人亚洲精品av一区二区| 三级经典国产精品| www.色视频.com| 丝袜美腿在线中文| 青春草国产在线视频 | 免费av不卡在线播放| 12—13女人毛片做爰片一| 国产精品免费一区二区三区在线| 别揉我奶头 嗯啊视频| 国产精品国产三级国产av玫瑰| 国产成人91sexporn| 男女啪啪激烈高潮av片| 91麻豆精品激情在线观看国产| 两个人视频免费观看高清| 又爽又黄无遮挡网站| 麻豆一二三区av精品| 欧美区成人在线视频| 五月伊人婷婷丁香| 黄色视频,在线免费观看| 国产黄色视频一区二区在线观看 | 国产美女午夜福利| 人妻久久中文字幕网| 少妇的逼水好多| 三级男女做爰猛烈吃奶摸视频| 三级国产精品欧美在线观看| 国产成人影院久久av| 最近的中文字幕免费完整| 中文在线观看免费www的网站| 久久久久免费精品人妻一区二区| 久久99精品国语久久久| 国产一区二区亚洲精品在线观看| 狂野欧美白嫩少妇大欣赏| 十八禁国产超污无遮挡网站| 亚洲成人av在线免费| 高清日韩中文字幕在线| 日韩亚洲欧美综合| 亚洲av成人精品一区久久| 国内精品美女久久久久久| 99在线人妻在线中文字幕| 国产精品久久视频播放| 久久草成人影院| av在线老鸭窝| 伦理电影大哥的女人| 欧美+亚洲+日韩+国产| 日韩三级伦理在线观看| 欧美激情在线99| 成人性生交大片免费视频hd| 久久人妻av系列| 看免费成人av毛片| 久久精品国产自在天天线| 深爱激情五月婷婷| 国产v大片淫在线免费观看| 久久久久免费精品人妻一区二区| 国产爱豆传媒在线观看| 日韩三级伦理在线观看| 能在线免费看毛片的网站| 成人亚洲精品av一区二区| 赤兔流量卡办理| 两个人的视频大全免费| 女同久久另类99精品国产91| 中文字幕av成人在线电影| 亚州av有码| 久久午夜福利片| 一卡2卡三卡四卡精品乱码亚洲| 小蜜桃在线观看免费完整版高清| 夜夜看夜夜爽夜夜摸| 亚洲国产精品sss在线观看| 一级毛片aaaaaa免费看小| 国产午夜福利久久久久久| 亚洲av男天堂| 久久精品国产亚洲av香蕉五月| 国内少妇人妻偷人精品xxx网站| 国产成人91sexporn| 国产精品久久视频播放| 直男gayav资源| 亚洲成人av在线免费| 丰满的人妻完整版| 亚洲五月天丁香| 尤物成人国产欧美一区二区三区| 欧美激情国产日韩精品一区| 久久人人精品亚洲av| 国产精品精品国产色婷婷| 麻豆av噜噜一区二区三区| 久久久久久久久久成人| 欧美+亚洲+日韩+国产| 黑人高潮一二区| 免费搜索国产男女视频| 极品教师在线视频| 波多野结衣巨乳人妻| 亚洲人成网站在线播放欧美日韩| 久久99热这里只有精品18| 亚洲高清免费不卡视频| 亚洲第一电影网av| 欧美最新免费一区二区三区| 国内揄拍国产精品人妻在线| 两个人视频免费观看高清| 精品日产1卡2卡| 搞女人的毛片| 国产精品久久久久久亚洲av鲁大| 国产精品日韩av在线免费观看| 午夜爱爱视频在线播放| 天堂av国产一区二区熟女人妻| 亚洲精品自拍成人| 青春草亚洲视频在线观看| 久久久国产成人精品二区| 少妇熟女aⅴ在线视频| 少妇人妻一区二区三区视频| 69人妻影院| 美女 人体艺术 gogo| 国产激情偷乱视频一区二区| 国产高潮美女av| 亚洲精品粉嫩美女一区| 在线免费十八禁| av专区在线播放| 在线观看av片永久免费下载| 精品久久国产蜜桃| 热99re8久久精品国产| 久久亚洲国产成人精品v| 亚洲自拍偷在线| 综合色丁香网| 中文字幕制服av| 久久久久久久久久久免费av| 国产蜜桃级精品一区二区三区| 一区二区三区免费毛片| 最近最新中文字幕大全电影3| 欧美精品一区二区大全| 99九九线精品视频在线观看视频| 91精品国产九色| 中国国产av一级| 又粗又硬又长又爽又黄的视频 | 麻豆av噜噜一区二区三区| 国产午夜精品一二区理论片| 夜夜看夜夜爽夜夜摸| 一级毛片aaaaaa免费看小| 国产在线男女| 日本黄色视频三级网站网址| 蜜臀久久99精品久久宅男| 亚洲人成网站在线观看播放| 黄色一级大片看看| 亚洲欧洲国产日韩| 18禁在线无遮挡免费观看视频| 九九爱精品视频在线观看| 大香蕉久久网| 成人国产麻豆网| 日韩欧美精品免费久久| 亚洲av成人精品一区久久| 12—13女人毛片做爰片一| 丰满乱子伦码专区| 国产69精品久久久久777片| 人妻制服诱惑在线中文字幕| 热99在线观看视频| 亚洲精品久久国产高清桃花| 欧美潮喷喷水| 日韩av不卡免费在线播放| 最后的刺客免费高清国语| 国产探花极品一区二区| 99九九线精品视频在线观看视频| 中文字幕av在线有码专区| 久久精品国产亚洲av香蕉五月| 日本黄色片子视频| 波多野结衣高清无吗| 亚洲av免费高清在线观看| 一个人看的www免费观看视频| 成人漫画全彩无遮挡| 精品国内亚洲2022精品成人| 男人舔女人下体高潮全视频| 如何舔出高潮| 国产 一区精品| 午夜激情福利司机影院| 欧美三级亚洲精品| 久久久久九九精品影院| 欧美区成人在线视频| 岛国在线免费视频观看| 日本熟妇午夜| 亚洲av一区综合| 美女被艹到高潮喷水动态| 久久中文看片网| 老司机福利观看| 国产精品久久久久久亚洲av鲁大| 久久精品夜色国产| 天天躁夜夜躁狠狠久久av| 99riav亚洲国产免费| 婷婷亚洲欧美| 成人午夜高清在线视频| 搞女人的毛片| 日本av手机在线免费观看| a级一级毛片免费在线观看| av免费在线看不卡| 国产久久久一区二区三区| 国产精品无大码| 三级毛片av免费| 免费观看在线日韩| 亚洲最大成人手机在线| 高清午夜精品一区二区三区 | 人人妻人人看人人澡| 亚洲av成人av| 青春草视频在线免费观看| 99久久中文字幕三级久久日本| 精品日产1卡2卡| 亚洲乱码一区二区免费版| 三级毛片av免费| 边亲边吃奶的免费视频| 人妻少妇偷人精品九色| 天天躁夜夜躁狠狠久久av| 亚洲av免费高清在线观看| 亚洲人成网站在线播放欧美日韩| 一本久久精品| 亚洲电影在线观看av| 美女内射精品一级片tv| 亚洲色图av天堂| 91在线精品国自产拍蜜月| 欧美日本视频| 偷拍熟女少妇极品色| 国产成人午夜福利电影在线观看| 久久久久久久午夜电影| 啦啦啦韩国在线观看视频| 99热只有精品国产| 国产亚洲av片在线观看秒播厂 | 人体艺术视频欧美日本| 国产精品久久视频播放| 日本成人三级电影网站| 午夜免费激情av| 国语自产精品视频在线第100页| 91在线精品国自产拍蜜月| 成人永久免费在线观看视频| 免费看美女性在线毛片视频| 欧美日韩国产亚洲二区| 日本av手机在线免费观看| 色5月婷婷丁香| 综合色av麻豆| 欧美日本视频| 日本免费a在线| 尾随美女入室| 一夜夜www| 精品少妇黑人巨大在线播放 | 国产单亲对白刺激| 免费观看a级毛片全部| 免费看av在线观看网站| 91av网一区二区| 99热这里只有是精品在线观看| 国产乱人视频| 亚洲欧美日韩高清专用| 我的老师免费观看完整版| 春色校园在线视频观看| 久久精品91蜜桃| 久久久久九九精品影院| 国产成人精品久久久久久| 毛片女人毛片| 欧美一区二区精品小视频在线| 精品一区二区三区视频在线| 国产爱豆传媒在线观看| 少妇熟女欧美另类| 国产精品无大码| 成人鲁丝片一二三区免费| 人人妻人人澡欧美一区二区| 干丝袜人妻中文字幕| 国产午夜精品久久久久久一区二区三区| 亚洲成a人片在线一区二区| 亚洲真实伦在线观看| 美女 人体艺术 gogo| 国产精品久久久久久亚洲av鲁大| 美女高潮的动态| 亚洲自拍偷在线| 久久精品国产99精品国产亚洲性色| 日韩欧美精品v在线| 亚洲成a人片在线一区二区| 只有这里有精品99| 国产单亲对白刺激| 一区福利在线观看| 亚洲国产精品国产精品| 国产成年人精品一区二区| 成人美女网站在线观看视频| 国产精品三级大全| 日韩在线高清观看一区二区三区| 一个人免费在线观看电影| 国产成人freesex在线| 岛国在线免费视频观看| 波野结衣二区三区在线| 99riav亚洲国产免费| 欧美日韩一区二区视频在线观看视频在线 | videossex国产| 欧美区成人在线视频| 性插视频无遮挡在线免费观看| 中国美白少妇内射xxxbb| 亚洲自偷自拍三级| 久久九九热精品免费| 国产老妇伦熟女老妇高清| 91狼人影院| 日韩欧美一区二区三区在线观看| 51国产日韩欧美| 国产精品一区二区在线观看99 | 大又大粗又爽又黄少妇毛片口| 国产一区二区在线观看日韩| 国产黄色视频一区二区在线观看 | 国产亚洲5aaaaa淫片| 毛片女人毛片| 三级经典国产精品| 超碰av人人做人人爽久久| 中文字幕久久专区| 免费av毛片视频| 国产69精品久久久久777片| 91久久精品电影网| 精品久久久久久久末码| 欧美人与善性xxx| 亚洲精品粉嫩美女一区| 又爽又黄a免费视频| 国产一区二区三区av在线 | 国产私拍福利视频在线观看| 啦啦啦观看免费观看视频高清| 亚洲国产欧美人成| 亚洲,欧美,日韩| 99国产精品一区二区蜜桃av| 亚洲中文字幕日韩| 亚洲人与动物交配视频| 深夜a级毛片| 成人亚洲精品av一区二区| 久久99热这里只有精品18| 真实男女啪啪啪动态图| 一级毛片电影观看 | 夜夜看夜夜爽夜夜摸| 久久久国产成人免费| 亚洲av中文av极速乱| 女人被狂操c到高潮| 午夜a级毛片| 中文在线观看免费www的网站| 中文精品一卡2卡3卡4更新| 亚洲av中文字字幕乱码综合| 淫秽高清视频在线观看| 亚洲国产精品合色在线| 日韩高清综合在线| 日韩人妻高清精品专区| 免费看av在线观看网站| 国产成人91sexporn| 国产色婷婷99| 欧美一区二区亚洲| 日韩在线高清观看一区二区三区| 久久久久久久午夜电影| 亚洲国产日韩欧美精品在线观看| 丝袜美腿在线中文| 欧美不卡视频在线免费观看| 日韩成人av中文字幕在线观看| 搞女人的毛片| 日韩高清综合在线| 亚洲精品亚洲一区二区| 免费看av在线观看网站| 久久久久久伊人网av| 99精品在免费线老司机午夜| 天美传媒精品一区二区| 在线天堂最新版资源| 99在线视频只有这里精品首页| 此物有八面人人有两片| 欧美极品一区二区三区四区| 在线免费观看的www视频| 麻豆一二三区av精品| 国产黄色小视频在线观看| 国产精品1区2区在线观看.| 日韩亚洲欧美综合| 日韩欧美精品v在线| 亚洲国产精品国产精品| 国产精品三级大全| 国产探花在线观看一区二区| 国产精品爽爽va在线观看网站| 亚洲最大成人中文| 亚洲最大成人手机在线| 久久99蜜桃精品久久| 欧美另类亚洲清纯唯美| 国产成年人精品一区二区| 免费av观看视频| 亚洲国产精品成人综合色| 日本-黄色视频高清免费观看| 最近2019中文字幕mv第一页| 国产不卡一卡二| 精品久久久久久久久久免费视频| a级毛片a级免费在线| 久久久久国产网址| 久久国内精品自在自线图片| 2021天堂中文幕一二区在线观| 国产91av在线免费观看| 哪里可以看免费的av片| 日本熟妇午夜| 亚洲av免费高清在线观看| av在线观看视频网站免费| 日本三级黄在线观看| 亚洲内射少妇av| 国产探花极品一区二区| 亚洲四区av| 神马国产精品三级电影在线观看| 国产久久久一区二区三区| 国产国拍精品亚洲av在线观看| 中文资源天堂在线| 99riav亚洲国产免费| 亚洲精品自拍成人| 日本在线视频免费播放| 国产一区二区亚洲精品在线观看| 国产精品一区二区三区四区免费观看| 久久久久网色| 日韩中字成人| 免费观看精品视频网站| 一区二区三区四区激情视频 | 性欧美人与动物交配| 久久精品国产清高在天天线| 91精品一卡2卡3卡4卡| 国产精品久久久久久av不卡| 久久热精品热| 九九爱精品视频在线观看| 精品久久久久久久久av| 亚洲最大成人中文| 久久久国产成人免费| 五月伊人婷婷丁香| 亚洲三级黄色毛片| 99国产精品一区二区蜜桃av| 一区二区三区高清视频在线| 国产91av在线免费观看| 国产精品日韩av在线免费观看| 少妇高潮的动态图| 能在线免费观看的黄片| 搡女人真爽免费视频火全软件| 青春草视频在线免费观看| 亚洲av第一区精品v没综合| 欧美三级亚洲精品| 一级av片app| 国产高清视频在线观看网站| 日本熟妇午夜| 在线免费十八禁| 久久精品国产自在天天线| 99精品在免费线老司机午夜| 色播亚洲综合网| 久久热精品热| 久久久久久久久久久丰满| 综合色av麻豆| 在线观看免费视频日本深夜| 一级黄片播放器| 九九久久精品国产亚洲av麻豆| 小蜜桃在线观看免费完整版高清| 色播亚洲综合网| 天天躁夜夜躁狠狠久久av| 高清毛片免费看| 淫秽高清视频在线观看| 插逼视频在线观看| 日韩一区二区三区影片| 久久精品国产清高在天天线| 国产伦在线观看视频一区| 国内精品宾馆在线| 精品久久久久久久久亚洲| 国产精品电影一区二区三区| 毛片一级片免费看久久久久| 日韩av不卡免费在线播放| 女人被狂操c到高潮| avwww免费| 国产又黄又爽又无遮挡在线| 床上黄色一级片| 少妇的逼水好多| 国产精品女同一区二区软件| 不卡视频在线观看欧美| 亚洲高清免费不卡视频| 色哟哟哟哟哟哟| 久久久久性生活片| 精品人妻视频免费看| 国产精品久久电影中文字幕| 日日撸夜夜添| 午夜精品一区二区三区免费看| 欧美性猛交黑人性爽| 国产免费男女视频| 六月丁香七月| 国产一级毛片七仙女欲春2| 亚洲成人久久爱视频| 亚洲av中文av极速乱| 成人亚洲精品av一区二区| 亚州av有码| 国产av麻豆久久久久久久| 亚洲国产欧美人成| 91狼人影院| 嫩草影院入口| 身体一侧抽搐| 久久精品国产亚洲网站| 久久国产乱子免费精品| 色尼玛亚洲综合影院| 长腿黑丝高跟| 国产伦一二天堂av在线观看| 夜夜爽天天搞| 精品一区二区三区人妻视频| 非洲黑人性xxxx精品又粗又长| 18禁黄网站禁片免费观看直播| 国产久久久一区二区三区| 国产精品一区二区三区四区久久| 午夜免费男女啪啪视频观看| 两个人视频免费观看高清| 中文亚洲av片在线观看爽| 熟女人妻精品中文字幕| 少妇人妻一区二区三区视频| 在线播放国产精品三级| 国产av不卡久久| 尾随美女入室| 精品久久久久久久末码| 边亲边吃奶的免费视频| 又爽又黄无遮挡网站| 婷婷六月久久综合丁香| 中文字幕久久专区| 成人亚洲欧美一区二区av| av女优亚洲男人天堂| 免费看a级黄色片| 直男gayav资源| 麻豆国产97在线/欧美| 特级一级黄色大片| 午夜老司机福利剧场| 久久久久久九九精品二区国产| 久久精品久久久久久久性| 男女做爰动态图高潮gif福利片| 国产色婷婷99|