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

    番茄中Rpi-blb2同源基因的挖掘與鑒定

    2016-05-27 01:50:45王俊生胡利宗劉紅占王永立
    廣西植物 2016年4期

    王俊生, 胡利宗, 劉紅占, 王永立, 陳 龍

    ( 周口師范學(xué)院 生命科學(xué)與農(nóng)學(xué)學(xué)院, 河南 周口 466001 )

    ?

    番茄中Rpi-blb2同源基因的挖掘與鑒定

    王俊生, 胡利宗, 劉紅占, 王永立, 陳龍

    ( 周口師范學(xué)院 生命科學(xué)與農(nóng)學(xué)學(xué)院, 河南 周口 466001 )

    摘要:番茄晚疫病是番茄生產(chǎn)中的主要病害之一,經(jīng)常會造成較大的經(jīng)濟(jì)損失。晚疫病生理小種的變異和進(jìn)化常會導(dǎo)致番茄品種原有的遺傳抗性喪失,因此不斷挖掘新的抗性基因,改良番茄晚疫病抗性是番茄抗病育種的長期任務(wù)。該研究采用BLAST同源比對的方法,以馬鈴薯野生近緣種的晚疫病抗性蛋白序列Rpi-blb2為種子序列,在NCBI蛋白質(zhì)序列數(shù)據(jù)庫中檢索得到11條番茄蛋白質(zhì)序列,這些序列與種子序列相似性為78%~83%,屬于番茄疾病抗性蛋白家族,并對該家族成員進(jìn)行了基因結(jié)構(gòu)、基因定位、序列保守結(jié)構(gòu)域和進(jìn)化關(guān)系等分析。結(jié)果表明:該家族中10條序列分布在第Ⅵ條染色體上,1條分布在第Ⅴ染色體上;6號染色體上的10序列呈現(xiàn)2個抗病基因簇分布,在染色體上分別占據(jù)2個和3個基因位點;10條同源蛋白是Rpi-blb2的共同垂直同源蛋白,但不具有平行同源關(guān)系,大多數(shù)成員定位于細(xì)胞質(zhì)。按照蛋白質(zhì)保守結(jié)構(gòu)域和基因定位的不同可分為三類,第一類共4條系列,包含有DUF3542和NB-ARC兩個保守結(jié)構(gòu)域特征序列;第二類共6條序列,與馬鈴薯Rpi-blb2蛋白一樣,僅包含NB-ARC保守結(jié)構(gòu)域特征序列,在這2類蛋白序列的NB-ARC結(jié)構(gòu)域均位于序列中部;第三類(僅包含XP_004239406.1)雖然也具有與第一類蛋白相似的DUF3542和NB-ARC結(jié)構(gòu)域,但在結(jié)構(gòu)域兩端的非保守區(qū)序列較短,且位于5號染色體上,因此將其單獨歸為1類。前兩類蛋白成員相應(yīng)的基因具有1~2個內(nèi)含子,第3類蛋白不含內(nèi)含子。該研究結(jié)果為利用生物技術(shù)選育番茄抗性品種提供了理論基礎(chǔ)。

    關(guān)鍵詞:番茄, 晚疫病, 基因, 垂直同源, 進(jìn)化分析

    番茄(Solanumlycopersicum)是在世界范圍內(nèi)廣泛種植的主要蔬菜之一,被各國人民所喜愛。番茄晚疫病(Phytophthora infestans)是番茄生產(chǎn)中的主要病害之一,可導(dǎo)致植株莖葉死亡和果實變褐甚至腐爛,每年在番茄主產(chǎn)區(qū)均造成較大的經(jīng)濟(jì)損失(Kamoun ,2001 ),已成為番茄生產(chǎn)的主要障礙之一。由于病源生理小種的不斷變異,導(dǎo)致生產(chǎn)上番茄品種的抗性逐漸喪失,加重了晚疫病病害的發(fā)生(趙統(tǒng)敏等,2006)。利用抗菌劑類藥物進(jìn)行防治不僅會造成環(huán)境污染,還會使病原菌產(chǎn)生抗耐藥性,促進(jìn)生理小種變異(Daayf & Platt, 2012),不利于晚疫病的控制,而選育抗性品種不僅可以從根本上克服番茄晚疫病的發(fā)生,而且經(jīng)濟(jì)有效。鑒定和克隆晚疫病抗性基因是選育抗性品種的基礎(chǔ),尤其一些廣譜主效基因在育種上的應(yīng)用會產(chǎn)生較好的效果。近年來,在番茄抗晚疫病基因克隆方面取得了較好的進(jìn)展,如Chunwongse et al(2002)、Moreau et al(1998)、Brouwer et al(2004)和Kole et al(2006)分別克隆了抗番茄晚疫病的幾個主效基因,如Ph-1、Ph-2、Ph-3和Ph-4,但隨后發(fā)現(xiàn)番茄中出現(xiàn)了能克服以上主效基因的晚疫病菌株(賈芝琪等,2009);Zhang et al(2013)把野生醋栗番茄中獲得的抗性基因Ph-3轉(zhuǎn)入番茄栽培種中使其獲得了廣譜晚疫病菌抗性,并且對其進(jìn)行了染色體精細(xì)定位。然而目前已有的番茄野生資源中,尚沒有發(fā)現(xiàn)更好的其他抗性基因,尤其廣譜性晚疫病抗性基因。

    從植物野生近緣種中挖掘鑒定抗性基因并導(dǎo)入栽培物種中,選育抗性植物品種是常用育種策略之一。近年來,從馬鈴薯野生種中克隆鑒定出了20多個抗性基因(周軍會,2008)。其中,Rpi-blb1和Rpi-blb2是從馬鈴薯野生種(Solanumbulbocastanum)中分離出的晚疫病廣譜抗性基因(Van der Vossen et al,2003,2005),并且Rpi-blb2基因被轉(zhuǎn)入番茄中獲得了抗馬鈴薯晚疫病菌株的番茄植株(Van der Vossen et al, 2005),證明了這種抗性可以在番茄和馬鈴薯之間進(jìn)行轉(zhuǎn)移,但是否可以抵抗番茄晚疫病生理小種尚未可知。馬鈴薯和番茄同屬茄科植物二倍體物種,均具有12對染色體,是比較基因組學(xué)研究的理想材料。因此,本研究以野生馬鈴薯抗病基因Rpi-blb2的蛋白序列為種子序列,基于番茄蛋白質(zhì)序列數(shù)據(jù)庫、基因組數(shù)據(jù)庫、遺傳圖譜和物理圖譜等生物信息資源,利用生物信息學(xué)方法挖掘番茄中與馬鈴薯Rpi-blb2基因同源的抗性蛋白(或基因)家族,并對其基因結(jié)構(gòu)信息、定位信息、序列進(jìn)化關(guān)系等進(jìn)行了較為全面的分析,為揭示該家族蛋白質(zhì)的結(jié)構(gòu)與功能關(guān)系,并利用相關(guān)信息進(jìn)行番茄抗晚疫病的分子育種奠定基礎(chǔ)。

    1材料與方法

    1.1 番茄抗性蛋白家族序列來源

    以馬鈴薯抗晚疫病基因Rpi-blb2(核酸序列Accession number DQ122125.1:7967bp)相應(yīng)蛋白序列(Accession number為AAZ95005.1,1267個氨基酸)作為種子序列,在NCBI(http://www.cbi.nlm.nih.gov/) 番茄非冗余蛋白質(zhì)數(shù)據(jù)庫進(jìn)行BLAST同源搜索,蛋白序列相似性設(shè)定為≥40%,相應(yīng)的核酸序列從蛋白質(zhì)序列檢索結(jié)果中超鏈接進(jìn)入。

    1.2 蛋白基本信息和基因結(jié)構(gòu)、基因定位信息獲取

    通過在蛋白比對結(jié)果網(wǎng)頁中,獲得比對序列之間的相似性和蛋白序列基本信息,鏈接目標(biāo)蛋白序列,可以獲得其相應(yīng)的核酸序列信息;通過相應(yīng)的GeneID超鏈接或者利用序列網(wǎng)頁中提供的Map Viewer超鏈接可以獲得相應(yīng)基因的定位信息和基因結(jié)構(gòu)信息、基因符號和基因功能描述等;也可以從相應(yīng)的mRNA序列(XM004245875.1)注釋信息中獲得基因定位識別號進(jìn)行超鏈接,從而獲得相應(yīng)基因定位信息和基因結(jié)構(gòu)信息。

    1.3 蛋白分子量、等電點和保守結(jié)構(gòu)域分析和亞細(xì)胞定位

    利用ExPASy網(wǎng)站提供的在線Prot-Param軟件(http://web.expasy.org/protparam/)進(jìn)行相對分子質(zhì)量和理論等電點(pI)等分析;利用CDD(http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) 在線軟件分析蛋白序列的保守結(jié)構(gòu)域;利用Plant-mPloc 軟件(http://www.csbio. sjtu.edu.cn/bioinf/plant-multi/)進(jìn)行亞細(xì)胞定位預(yù)測。

    1.4 蛋白質(zhì)序列進(jìn)化樹構(gòu)建

    通過Conserved Domain(http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi)在線工具,搜索蛋白質(zhì)序列保守域;在BLAST-protein結(jié)果頁,選擇所有同源目標(biāo)序列后,利用Tree View功能顯示比對序列與同源序列的距離進(jìn)化樹(采用Neighbor Joining法,kimura protein 距離,序列最大差異值設(shè)定為0.7,有根樹顯示);下載比對結(jié)果中的目標(biāo)序列到本地硬盤,采用DNAman軟件進(jìn)行多序列比對(采用缺省參數(shù)),分析番茄SLy-Rpi-blb2蛋白家族成員之間進(jìn)化關(guān)系。

    2結(jié)果與分析

    2.1 番茄Sly-Rpi-blb2蛋白家族與馬鈴薯Rpi-blb2蛋白是垂直同源關(guān)系

    2.1.1 番茄同源抗性蛋白序列基本信息以野生馬鈴薯抗晚疫病基因Rpi-blb2(核酸序列Accession number 為DQ122125.1:7967 bp)相應(yīng)蛋白序列(Accession number為AAZ95005.1,1 267個氨基酸)作為種子序列,在番茄非冗余蛋白質(zhì)數(shù)據(jù)庫進(jìn)行BLASTp同源搜索,獲得了18條相似性在40%以上的同源序列,其中與11條序列的相似性為78%~83%,11條序列之間的相似性為85%~100%。進(jìn)一步比對查詢這11條序列的相關(guān)信息(表1),發(fā)現(xiàn)這些序列長度不等,最長為1 257 個氨基酸,最短566 個氨基酸,其中有10條定位于第6號染色體,1條定位于5號染色體,均與番茄抗病性(晚疫病和根結(jié)線蟲病抗性)密切相關(guān);位于6號染色體上的10條序列分子量接近,等電點分布在4.85~5.19之間,非常接近,表明這些蛋白質(zhì)以酸性氨基酸為主,作用環(huán)境偏酸性。因此,推測這些蛋白序列可能與野生馬鈴薯Rpi-blb2蛋白互為共同垂直同源蛋白(Coorthologs),他們之間互為平行同源蛋白。位于5號染色體上的XP_004239406.1序列長度僅566個氨基酸,推測該蛋白的祖先基因序列可能發(fā)生了部分復(fù)制、轉(zhuǎn)移和插入等進(jìn)化事件。因此,把這11條蛋白質(zhì)序列暫命名為SLy-Rpi-blb2蛋白亞家族。

    2.1.2 番茄SLy-Rpi-blb2蛋白家族基因定位和亞細(xì)胞定位基因定位結(jié)果(表2,圖1)表明,該蛋白家族分布在6號染色體上的10個成員定位在2.3 M~2.75 M之間,其中8個成員分布在5個基因位點,并且集中分布在2個區(qū)段,成簇分布,其中第一族包括3條序列(XP_004240499.1、NP_001234063.1和XP_004240498.1),并且XP_004240499.1和NP_001234063.1占據(jù)基因組上相同位點,這與兩條序列來自不同抗性品種相關(guān);第二簇包括5條序列,分別為ABI96216.1、XP_004240523.1、AAC97933.1、NP-001234622.1和ABI96218.1,其中后3條序列位于基因組上同一位點;另外2個成員(AAC32253.1和AAC32252.1)未能定位于染色體具體位置。該蛋白家族中的各成員在基因組成上比較簡單(表2),僅有1~3個外顯子數(shù)。采用Plant-mPloc在線軟件進(jìn)行亞細(xì)胞定位結(jié)果表明,11條蛋白序列均定位于細(xì)胞質(zhì)中,這與野生馬鈴薯Rpi-blb2蛋白的亞細(xì)胞定位結(jié)果相同。

    2.1.3 番茄SLy-Rpi-blb2蛋白家族的保守結(jié)構(gòu)域分析 對11條同源蛋白序列通過Conserved Domain在線工具,搜索其保守域,結(jié)果表明這些蛋白的精確保守域共有3類,第一類共4條系列,包括XP_004240523.1、AAC97933.1、NP_001234063.1和AAC32252.1,同時包含有DUF3542和NB-ARC兩個保守結(jié)構(gòu)域特征序列,第二類包括XP_004240499.1、NP_001234622.1、AAC32253.1、ABI96218.1、ABI96216.1和XP_004240498.1,僅包含NB-ARC保守結(jié)構(gòu)域特征序列,并且與馬鈴薯Rpi-blb2蛋白保守結(jié)構(gòu)一致(圖1:B),因此推測第二類蛋白與馬鈴薯Rpi-blb2蛋白(AAZ95005.1)的進(jìn)化距離最近,與這類蛋白的祖先序列最為接近,分化較早,在這2類蛋白序列中,NB-ARC結(jié)構(gòu)域均位于序列中段;第三類(僅包含一條序列XP_004239406.1)雖然也具有DUF3542和NB-ARC結(jié)構(gòu)域,并且2個結(jié)構(gòu)域間距與第一類非常接近,但在結(jié)構(gòu)域兩端的其他序列非常短,并且該序列位于5號染色體上,因此將其單獨歸為1類。

    表 1 馬鈴薯Rpi-blb2蛋白在番茄中的同源蛋白序列信息

    表 2 番茄SLy-Rpi-blb2蛋白家族基因定位與結(jié)構(gòu)信息

    注: “—”表示未能獲得相關(guān)信息。

    Note: “—”indicates failed to get relevant information.

    圖 1 番茄SLy-Rpi-blb2蛋白家族成員在基因組上的分布位置Fig. 1 Location for members of SLy-Rpi-blb2 protein family in tomato genome

    2.2 番茄SLy-Rpi-blb2蛋白家族成員的序列比對、基因結(jié)構(gòu)與進(jìn)化

    為了進(jìn)一步確定SLy-Rpi-blb2蛋白家族成員之間以及與馬鈴薯Rpi-blb2的同源關(guān)系,利用馬鈴薯Rpi-blb2蛋白與番茄蛋白數(shù)據(jù)庫中相似性達(dá)到40%以上的18條抗性蛋白(包括晚疫病抗性蛋白)序列構(gòu)建拓?fù)浣Y(jié)構(gòu)樹,結(jié)果(圖3)表明,野生馬鈴薯Rpi-blb2 蛋白(圖3中紅色箭頭所指)與表1中的11條番茄蛋白序列歸為一大類(圖3:A中蛋白家族I),其他7條序列分為蛋白家族II和蛋白家族Ⅲ。表明這11條序列與馬鈴薯晚疫病抗性蛋白序列Rpi-blb2可能為共同垂直同源蛋白。將11條序列與種子序列進(jìn)行多序列比對分析(圖4),結(jié)果表明,有4條序列(分別為ABI96218.1(996 個氨基酸)、AAC32253.1(1204 個氨基酸)、AAC32252.1(1206 個氨基酸)和XP_004239406.1(566 個氨基酸)的長度與種子序列差異較大,其他7條與番茄Rpi-blb2序列基本可以完全匹配,這7條序列與馬鈴薯Rpi-blb2蛋白序列相比,在N端存在2個明顯的較短缺失序列(圖4:Ⅰ,缺失3個氨基酸;圖4:Ⅱ,缺失6個氨基酸);ABI96218.1與ABI96216.1相比,最大差異在于C端缺失了249個氨基酸序列,這可能是兩者的祖先基因因為部分染色體片段缺失造成的,也可能是堿基突變造成編碼序列提前終止引起的。基因結(jié)構(gòu)(圖3:B)表明,該家族成員有3種基因結(jié)構(gòu)類型,第一類含有2個內(nèi)含子(5條序列),第二類含有一個較小的內(nèi)含子(2條序列),第三類不具有內(nèi)含子結(jié)構(gòu)(XP_004239406),比較這些基因內(nèi)含子序列發(fā)現(xiàn),第一個內(nèi)含子序列相對較長,且長度差異較大(556~2 199 bp),序列間相似性差異也較大;第二內(nèi)含子位于靠近蛋白質(zhì)N段的編碼區(qū)內(nèi),所有同源蛋白相應(yīng)核酸序列的第二內(nèi)含子長度均為75 bp,且序列高度保守,這表明第二內(nèi)含子是一個相當(dāng)保守的區(qū)域,可能對這類蛋白具有重要作用,因此進(jìn)化過程中不易發(fā)生改變。

    3討論

    通過對生物信息數(shù)據(jù)進(jìn)行加工獲取有用知識,指導(dǎo)科研方向,是生物信息學(xué)常用的方法。馬鈴薯和番茄同屬茄科,進(jìn)化關(guān)系密切,因此利用馬鈴薯晚疫病抗性基因Rpi-blb2序列,從番茄蛋白數(shù)據(jù)庫中獲得了一個包含11條蛋白序列的蛋白家族SLy-Rpi-blb2,通過鑒定該家族成員與Rpi-blb2的進(jìn)化關(guān)系、基因結(jié)構(gòu)、蛋白保守域特性等,確定該家族與種子序列是垂直同源關(guān)系。隨后對家族成員之間的進(jìn)化事件進(jìn)行了分析,表明該家族的共同祖先基因在漫長的進(jìn)化歷史中發(fā)生了基因復(fù)制、重排、可變剪接以及水平轉(zhuǎn)移事件,從而形成了較為豐富的后代同源基因(蛋白),但與祖先基因(蛋白)序列相比,可能由于保守域序列的改變,從而導(dǎo)致不同成員在抗病的能力上發(fā)生了改變。

    3.1 番茄SLy-Rp-Blb2蛋白家族與Rpi-blb2可能是垂直同源關(guān)系

    NBS是指核苷酸結(jié)合位點,植物疾病抗性蛋白中的NBS是一個大約320個氨基酸的高度保守序列,具有這種保守序列的蛋白可能是利用蛋白質(zhì)構(gòu)象變化調(diào)節(jié)NTP的水解,從而通過信號轉(zhuǎn)導(dǎo),使寄主對病原菌正確識別,促進(jìn)寄主抗性防御信號的啟動(Belkhadir et al, 2004; Leipe et al, 2004)。在植物抗性蛋白質(zhì)中NB-ARC結(jié)構(gòu)域就是NBS-LRR保守結(jié)構(gòu)域(Stuart Seah et al, 2007)。本研究中的種子序列Rpi-blb2(Van der Vossen et al, 2005) 具有顯著的NB-ARC保守結(jié)構(gòu)域序列特征,屬于AAA-ATPase 超級家族(圖2),可能與這類蛋白參與病毒產(chǎn)生的植物過敏反應(yīng)或程序性死亡有關(guān)(Van der Biezen1,1998)。DUF3542保守域序列的功能尚未可知,但具有DUF3542結(jié)構(gòu)域的蛋白序列長度一般在516~1 284個氨基酸,并且這類蛋白多與NB-ARC結(jié)構(gòu)域有關(guān)(http://www.ncbi.nlm.nih.gov/ Structure/cdd/wrpsb.cgi?INPUT_TYPE=live&SEQUENCE=XP_004240523.1)。在番茄SLy-Rp-Blb2蛋白家族中,有7條序列與馬鈴薯Rp-blb2蛋白具有相同的保守結(jié)構(gòu)域,僅具有NB-ARC結(jié)構(gòu)域序列特征,另4條序列同時具有NB-ARC和DUF3542結(jié)構(gòu)域序列特征。通過對馬鈴薯Rpi-blb2蛋白和番茄蛋白數(shù)據(jù)庫中相似性較高(40%以上)的同源蛋白構(gòu)建序列進(jìn)化樹(圖3)和亞細(xì)胞定位分析(表2),進(jìn)一步證明番茄SLy-Rp-Blb2蛋白家族與野生馬鈴薯Rp-blb2蛋白可能是垂直同源關(guān)系,由共同祖先種通過物種分化形成的。

    圖 2 番茄SLy-Rpi-blb2蛋白家族成員結(jié)構(gòu)域Fig. 2 Conserved domain for tomato SLy-Rpi-blb2 protein family

    圖 3 番茄SLy-Rpi-blb2的系統(tǒng)進(jìn)化樹(A)及SLy-Rpi-blb2家族成員基因結(jié)構(gòu)(B) 基因結(jié)構(gòu)圖中,上部粗線段表示mRNA序列,細(xì)線段表示基因組中內(nèi)含子序列,下部粗線表示蛋白編碼區(qū)。Fig. 3 Phylogenetic tree of tomato SLy-Rpi-blb2 protein family Above blod lines indicate mRNA, slim lines indicate intron in genome, below blod lines indicate coding sequence.

    圖 4 番茄SLy-Rpi-blb2蛋白家族成員序列比對Fig. 4 Aligned multi-sequence comparison for members of tomato SLy-Rpi-blb2 protein subfamily

    在SLy-Rpi-blb2蛋白家族中,序列AAC97933.1(Mi-1)及其抗性蛋白同源序列(NP_001234063.1、NP_001234622.1、AAC32252.1和AAC32253.1)所在的染色體區(qū)段是利用胚胎挽救技術(shù)從茄科物種Solanumperuvianum中漸滲入番茄栽培種中的(Smith, 1944),且導(dǎo)入的染色體片段位于番茄6號染色體短臂,而后通過育種實踐和長期自然進(jìn)化,把滲入的染色體片斷減少到了約650 kb的區(qū)間,這一區(qū)間的抗性基因和不良性狀的緊密連鎖導(dǎo)致抗性基因與其它品種優(yōu)良基因重組率嚴(yán)重降低(Liharska et al, 1996;Kaloshian et al,1998;Schulte et al,2006)。但雜交和回交等技術(shù)的頻繁利用,使?jié)B入了攜帶Mi-1.1基因的染色體區(qū)段不斷發(fā)生染色體內(nèi)部重組、復(fù)制、重排以及染色體之間的交換,從而導(dǎo)致現(xiàn)在的抗性番茄品種(攜帶Mi-1.1基因區(qū)段)在這一區(qū)段與不具有Mi-1.1基因位點的易感根結(jié)線蟲品種相應(yīng)區(qū)段不完全一樣,這一點已得到證明(Kuang et al,2005; Stuart et al, 2007)。因此,番茄SLy-Rpi-blb2蛋白家族中的序列不再是簡單的平行同源關(guān)系。從系統(tǒng)樹(圖3:A,有根樹,距離長短表示進(jìn)化速率快慢)的末端節(jié)點到共同祖先節(jié)點距離來看,番茄11條序列的共同祖先序列為圖3中節(jié)點B代表的序列,從節(jié)點B到11條末端序列進(jìn)化中,只有序列XP_004240523.1進(jìn)化距離最短,并且沒有產(chǎn)生序列分化,其他10條序列到節(jié)點B之間存在相對復(fù)雜的進(jìn)化過程(存在較多的節(jié)點),因此推測XP_004240523.1是番茄中這類蛋白的祖先序列(或最原始序列),其他序列可能是由該序列進(jìn)化而來。當(dāng)然這僅是根據(jù)構(gòu)建的距離進(jìn)化樹推斷,如果改變進(jìn)化樹構(gòu)建方法或考慮基因結(jié)構(gòu)及相應(yīng)蛋白序列的保守域特性,可能會有不同的推斷結(jié)論。

    3.2 番茄SLy-Rpi-blb2蛋白在基因組上成簇分布

    植物基因組中抗性基因一般具有成簇存在的特點(Hulbert et al, 2001; Cannon et al, 2002),這類基因編碼的蛋白的一個主要功能就是作為檢測分子可以偵測到特定病原物和害蟲的感染。Van der Vossen et al(2003)研究表明,在野生馬鈴薯Rpi-blb1基因(RGA2-blb)座位兩側(cè)分布有RGA3-blb、RGA1-blb和RGA4-blb同源基因(互為平行同源基因),這些基因編碼的平行同源蛋白的序列相似性很高(70%~81%),并且Rpi-blb1基因是祖先基因RGA1-blb和RGA3-blb通過重組事件進(jìn)化而來。已有證據(jù)表明(Aaron & Barbara, 2007)抗性基因簇可以通過重組過程中的錯配,產(chǎn)生新的抗性基因序列,改變抗性特異性(抵抗新的病毒生理小種)或改變抗性基因的表達(dá)模式,而這一特性促進(jìn)了抗性基因的進(jìn)化;另外,抗性基因成簇分布,利于這些平行基因之間的序列交流,從而促進(jìn)和加速了抗性基因的進(jìn)化。

    在茄科植物不同物種的6號染色體短臂區(qū)均分布著抗性基因集合,執(zhí)行對不同病原菌的抗性(Stuart et al,2007)。本研究獲得的番茄SLy-Rpi-blb2蛋白家族共11條序列,10條成簇分布在6號染色體的短臂序列 (2 327 335~2 750 696 bp,總長624 kb)中,在染色體上的位置較靠近,執(zhí)行不同的抗性功能(表1)。其中第一簇有3條序列,占2個基因位點(圖1,Cluster 1p),在Mi-1B基因上游的基因位點LOC101265033長度為3 191 bp,是一個假定的晚疫病抗性蛋白同源序列,但不具有編碼蛋白功能,屬于假基因;第二簇有5條序列(Cluster 2p),占3個基因位點,這一簇中,還分布兩個基因位點,分別是LOC1012267122和LOC101267412,其相應(yīng)的蛋白序列均與運輸抑制響應(yīng)蛋白高度相似,這可能與該簇基因?qū)膊】剐缘捻憫?yīng)密切相關(guān)。另外2條AAC32252.1和AAC32253.1未能定位到具體的位置,但這2條序列來源于抗根結(jié)線蟲栽培品種Motelle,這一品種中的抗性基因均來源于茄科物種Solanumperuvianum中(Smith, 1944),這一結(jié)論與Stuart et al(2007)基本一致。

    番茄SLy-Rpi-blb2蛋白家族中11個成員間的序列相似性在85%以上,來自于抗性品種(Motelle或VFNT Cherry)的成員有5個,包括AAC97933.1、NP_001234063.1、NP_001234622.1、AAC32253.1和AAC32252.1,其他6個成員來自于易感根結(jié)線蟲品種(Heinz 1706),因此部分成員在基因組被定位在相同位置(圖1),如NP_001234063.1和XP_004240499.1分享同一基因位點,AAC97933.1、NP_001234622.1和ABI96218.1分享同一基因位點。這些蛋白在番茄抗病反應(yīng)中的功能目前并不完全清楚。Mi-1蛋白序列(AAC97933.1)與馬鈴薯Rpi-blb2具有82%的相似性(van der Vossen et al, 2005),近來被證明是在ETR3蛋白參與下,通過感知荷爾蒙來介導(dǎo)番茄植株對根結(jié)線蟲的抗性(Mantelin et al,2013),其同源序列Mi-1.2被證明具有根結(jié)線蟲抗性的功能(Nombela et al, 2003),其他同源蛋白由于沒有相應(yīng)的獨立表現(xiàn)型,因此對其功能的鑒定比較困難(Stuart et al, 2007)。但由于可變剪接產(chǎn)生的較短的抗性蛋白質(zhì)同源序列,如本研究中的ABI96218.1(Mi-1G)編碼996 個氨基酸,具有完整N段區(qū)域,能反向調(diào)節(jié)Mi-1.2(Stuart et al, 2007),從而調(diào)節(jié)番茄抗病性。定位于5號染色體上的蛋白質(zhì)序列XP_004239406.1(566 個氨基酸),與其他10條序列相比,該序列在N段和C段分別缺失了大約500 個氨基酸和100個氨基酸的片段(圖2和圖4),但該序列仍保留了2個保守域的序列特征(圖2),從這一點來看,該序列與圖2中第一類序列應(yīng)歸為同一類;該基因編碼區(qū)序列是完整的,這意味著該基因(蛋白)可能具有新的功能,在后續(xù)功能研究中應(yīng)作為重點研究之一,然而該基因與其他同源基因的進(jìn)化關(guān)系目前尚不能完全確定。

    本研究獲得的同源蛋白序列(基因),大多數(shù)功能尚未被鑒定,因此在后續(xù)研究中可以重點鑒定其表達(dá)模式和性狀表現(xiàn),以分析和鑒定其功能,為番茄抗病分子輔助育種提供理論依據(jù)。

    參考文獻(xiàn):

    AARON RF,BARBARA JB, 2007. The evolution of resistance genes in multi-protein plant resistance systems [J]. Curr Opin Genet & Dev,17:493-499.

    BELKHADIR Y,SUBRAMANIAM R,DANGL JL, 2004. Plant disease resistance protein signaling: NBS-LRR proteins and their partners [J]. Curr Opin Plant Biol(4),7:391-399.BROUWER DJ,JONES ES,CLAIR DAS, 2004. QTL analysis of quantitative resistance toPhytophthorainfestans(late blight) in tomato and comparisons with potato [J]. Genome,47(3): 475-492.

    CANNON SB,ZHU HY,BAUMGARTEN AM,et al, 2002. Diversity,distribution and ancient taxanomic relationships within the TIR and non-TIR NBS-LRR resistance gene subfamilies [J]. J Mol Evol,54(4):548-562.

    Chunwongse J,Chunwongse C,Black L,et al, 1998. Molecular mapping of thePh-3 gene forlatebligh resistance in tomato [J]. J Hortic Sci Biotechnol,77(3): 281-286.

    DAAYF F,PLATT HW, 2012. Changes in metalaxyl resistance among glucose phosphate isomerase genotypes ofPhytophthorainfestansin Canada during 1997 and 1998 [J]. Am J Potato Res,77(5): 311-318.

    HULBERT SH,WEBB CA,SMITH SM,et al, 2001. Resistance gene: complexes: evolution and utilization [J]. Ann Rev Phytopathol,39:285-312.

    JIA ZQ,CUI YH,LI Y,et al, 2009. Expression the potato late blight resistant geneR3a,R1 andRBin tomato [J]. Acta Hortic Sin,36(8):1 153-1 160. [賈芝琪,崔艷紅,李穎,等, 2009. 馬鈴薯抗晚疫病基因R3a、R1和RB在番茄中的表達(dá) [J]. 園藝學(xué)報,36(8):1 153-1 160.]

    KAMOUN S, 2001. Non-host resistance toPhytophthora: novel prospects for a classical problem [J]. Curr Opin Plant Biol,4(4): 295-300.

    KOLE C,ASHRAFI H,LIN G,et al, 2006. Identification and molecular mapping of a new R gene,Ph-4,conferring resistance to late blight in tomato [C]. Solanaceae Conf,Univ of Wisconsin,Madison,Abstr44.

    KALOSHIAN I,YAGHOOBI J,LIHARSKA T,et al, 1998. Genetic and physical localization of the root-knot nematode resistance locus Mi in tomato [J]. Mol Gen Genet,257(3):376-385.KUANG H,WEI F,MARANO MR,et al, 2005. The R1 resistance gene cluster contains three groups of independently evolving,type R1 homologues and shows substantial structural variation among haplotypes ofSolanumdemissum[J]. Plant J,44(1): 37-51.

    LEIPE DD,KOONIN EV,ARAVIND L, 2004. STAND,a class of P-loop NTPases including animal and plant regulators of programmed cell death: multiple,complex domain architectures,unusual phyletic patterns,and evolution by horizontal gene transfer [J]. J Mol Biol,343(1):1-28.LIHARSKA TB,KOORNNEEF M,VAN WORDRAGEN M,et al, 1996. Tomato chromosome 6: effect of alien chromosomal segments on recombinant frequencies [J]. Genome,39(3):485-491.

    MANTELIN S,BHATTARAI KK,JHAVERI TZ,et al, 2013. Mi-1-mediated resistance toMeloidogyneincognitain tomato may not rely on ethylene but hormone perception through ETR3 participates in limiting nematode infection in a susceptible host [J]. PLoS ONE,8(5):63 281.

    MOREAU P,THOQUET P,OLIVER J,et al, 1998. Genetic mapping ofPh-2,a single locus controlling partial resistance toPhytophthorainfestansin tomato [J]. Mol Plant Micr Interact,11(4): 259-269.

    NOMBELA G,WILLIAMSON VM,MUNIZ M, 2003. The root-knot nematode resistance geneMi-1.2 of tomato is responsible for resistance against Whitefly Bemisia tabaci [J]. Mol Plant Microl Interact,16(7):645-649.

    SMITH PG, 1944. Embryo culture of a tomato species hybrid. Proc Am Soc Hortic Sci,44:413-416.

    SCHULTE D,CAI DG,KLEINE M,et al, 2006. A complete physical map of a wild beet (Beta procumbens ) translocation in sugar beet [J]. Mol Gen Genom,275(5):504-511.

    ZHANG C,LIU L,ZHENG Z,et al, 2013. Fine mapping of thePh-3 gene conferring resistance to late blight (Phytophthorainfestans) in tomato [J]. Theor Appl Gene,126(10):2 643-53.

    ZHOU JH, 2008. Advancement of late blight resistance genes in potato [J]. Biotechnol Bull,5: 13-17. [周軍會,2008. 馬鈴薯晚疫病抗病基因研究進(jìn)展 [J]. 生物技術(shù)通報,5: 13-17.]

    STUART S,ADAM C,TELLEEN,et al, 2007. Introgressed and endogenousMi-1 gene clusters in tomato diver by complex rearrangements in flanking sequences and show sequence exchange and diversifying selection among homologue [J]. Theor Appl Genet,114(7):1 289-1 302.

    VAN DER VOSSEN EA,GROS J,SIKKEMA A,et al, 2005. TheRpi-blb2 gene from solanum bulbocastanum is anMi-1 gene homolog conferring broad-spectrum late blight resistance in potato [J]. Plant J,44(2): 208-222.

    VAN DER BIEZEN1 EA,JONES JD, 1998. The NB-ARC domain: a novel signaling motif shared by plant resistance gene products and regulators of cell death in animals [J]. Curr Biol,8(7): 226-228.

    VAN DER VOSSEN EA,SIKKEMA A,HEKKERT B,et al, 2003. An ancient R gene from the wild potato species solanum bulbocastanum confers broad-spectrum resistance to phytophthora infestans in cultivated potato and tomato [J]. Plant J,36(6): 867-882.ZHAO TM,ZHOU CY,YU WG,et al, 2006. Research progress on tomato late blight and resistance breeding [J]. Jiangsu J Agric Sci,4(2): 175-180. [趙統(tǒng)敏,鄒茶英,余文貴,等, 2006. 番茄晚疫病及其抗病育種研究 [J]. 江蘇農(nóng)業(yè)學(xué)報,4(2): 175-180.]

    Discovery and identification of homology gene toRpi-blb2 in tomato (Solanumlycopersicum)

    WANG Jun-Sheng, HU Li-Zong, LIU Hong-Zhan, WANG Yong-Li, CHEN Long

    (CollegeofLifeSciencesandAgriculture,ZhoukouNormalUniversity, Zhoukou 466001, China )

    Abstract:Late blight is the most destructive disease in tomato (Solanum lycopersicum) cultivation worldwide and will result in large economic losses. As the variation and evolution of original races, mining new resisitance gene and improving the genetic resistance to late blight is a major issue in breeding new varieties of tomato , because of the new races will overcome the resistance gene-mediated genetic resistance from commercial varieties. Rpi-blb2, a potato late blight resistance protein as seed sequence, has 11 homology sequences from tomato protein sequences databank in NCBI and has a similarity to seed sequence range from 78%-83%. The above data were obtained by blast-protein way in this study. Then the gene structure, gene location, conserved domain feature of 11 protein sequences, and homology sequences of Rpi-blb2 belong to disease resistance protein sup-family of tomato and their evolution relationship were revealed and discussed. The results showed that ten of eleven distributed in the sixth chromosome, one the fifth chromosome; ten sequences distributed the sixth chromosome showed two gene clusters, containing two and three gene locus in the sixth chromosome, respectively. They were co-orthologous for Rpi-blb2 protein, but not paralogous each other, the most of them located in cytoplasm. According to the difference in structure domain and gene location, they were divided into three categories. The first category, a total of four members, contains DUF3542 and NB-ARC conser domain; The second category, a total of six members, like potatoes Rpi-blb2 protein, containing only NB-ARC conserved domain, the NB-ARC domain is located in the middle of these two kinds of protein sequences; The third category (only contains XP_004239406. 1) also has DUF3542 and NB-ARC domain and located on chromosome 5, but the non-conserved sequence at the ends of the domain is shorter. The members of the first two categories have one or two introns in gene structure, the third category does not contain introns. The study provides theory basis for breeding tomato disease resistance cultivar by molecular biology techniques.

    Key words:tomato (Solanum lycopersicum), late blight, gene, orthologous, phylogenetic analysis

    中圖分類號:Q943.2,S641.2

    文獻(xiàn)標(biāo)識碼:A

    文章編號:1000-3142(2016)04-0462-09

    作者簡介:王俊生(1970-),男,陜西蒲城人,博士,副教授,主要從事生物信息教學(xué)與植物分子育種,(E-mail)wjs0099@163.com。

    基金項目:河南省科技廳基礎(chǔ)與前沿技術(shù)研究項目 (122300410127); 河南省教育廳項目 (14B210022) [Supported by the Basic and Forefront Technology Research from Henan Office of Science and Technology (122300410127); the Project of Henan Education Office (14B210022)]。

    *收稿日期:2014-11-12修回日期: 2015-04-03

    DOI:10.11931/guihaia.gxzw201408016

    王俊生,胡利宗,劉紅占,等. 番茄中Rpi-blb2同源基因的挖掘與鑒定[J]. 廣西植物, 2016, 36(4):462-470

    WANG JS,HU LZ,LIU HZ,et al. Discovery and identification of homology gene toRpi-blb2 in tomato (Solanumlycopersicum) [J]. Guihaia, 2016, 36(4):462-470

    91精品伊人久久大香线蕉| 一级爰片在线观看| 蜜桃国产av成人99| 日韩免费高清中文字幕av| 男男h啪啪无遮挡| 少妇猛男粗大的猛烈进出视频| 国产成人a∨麻豆精品| 一边摸一边做爽爽视频免费| 狠狠精品人妻久久久久久综合| 午夜视频国产福利| 一级毛片 在线播放| av国产精品久久久久影院| 黄网站色视频无遮挡免费观看| 男女无遮挡免费网站观看| 成人亚洲欧美一区二区av| 99香蕉大伊视频| 青春草亚洲视频在线观看| 9热在线视频观看99| 新久久久久国产一级毛片| 看免费成人av毛片| 欧美日韩视频精品一区| 一区二区三区四区激情视频| 久久久久久久大尺度免费视频| 宅男免费午夜| 国产极品天堂在线| 黄色怎么调成土黄色| 99热全是精品| 亚洲精品美女久久久久99蜜臀 | 美女国产高潮福利片在线看| 午夜免费观看性视频| 久久人妻熟女aⅴ| 中国三级夫妇交换| 中文乱码字字幕精品一区二区三区| 三上悠亚av全集在线观看| 日韩,欧美,国产一区二区三区| 在线观看免费日韩欧美大片| 国产成人精品婷婷| 一本—道久久a久久精品蜜桃钙片| 精品少妇黑人巨大在线播放| 免费女性裸体啪啪无遮挡网站| 国产黄频视频在线观看| 精品久久蜜臀av无| √禁漫天堂资源中文www| 国产日韩一区二区三区精品不卡| 一本久久精品| 香蕉精品网在线| 久久毛片免费看一区二区三区| 亚洲av电影在线进入| 亚洲欧美一区二区三区黑人 | 欧美激情极品国产一区二区三区 | 高清不卡的av网站| 亚洲美女搞黄在线观看| 99视频精品全部免费 在线| 视频区图区小说| 啦啦啦视频在线资源免费观看| 欧美成人精品欧美一级黄| 日本黄色日本黄色录像| 久久久久精品久久久久真实原创| 国产在视频线精品| 久久精品aⅴ一区二区三区四区 | 亚洲av男天堂| 精品一区二区三区视频在线| kizo精华| 欧美日韩精品成人综合77777| 精品少妇久久久久久888优播| 国产精品久久久久久久久免| 国产高清不卡午夜福利| 高清av免费在线| a级毛片黄视频| 免费人成在线观看视频色| 欧美精品一区二区大全| 久久影院123| 超碰97精品在线观看| 亚洲国产日韩一区二区| 成人国产av品久久久| 精品久久蜜臀av无| 国产精品国产三级国产专区5o| 国产成人午夜福利电影在线观看| 欧美人与性动交α欧美精品济南到 | 在线精品无人区一区二区三| 亚洲国产av新网站| 亚洲国产色片| 欧美人与善性xxx| 如日韩欧美国产精品一区二区三区| 综合色丁香网| 日韩制服丝袜自拍偷拍| 制服丝袜香蕉在线| 国产男女内射视频| 最近2019中文字幕mv第一页| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 五月伊人婷婷丁香| 亚洲精品av麻豆狂野| 中文字幕免费在线视频6| 性高湖久久久久久久久免费观看| 国产一区亚洲一区在线观看| av女优亚洲男人天堂| 国产成人午夜福利电影在线观看| 99热国产这里只有精品6| 丝袜美足系列| 久久狼人影院| 国产精品一国产av| 午夜老司机福利剧场| 国产精品麻豆人妻色哟哟久久| 国产成人精品久久久久久| 国产国拍精品亚洲av在线观看| 国产成人aa在线观看| a级片在线免费高清观看视频| 亚洲美女黄色视频免费看| 亚洲 欧美一区二区三区| 天天躁夜夜躁狠狠躁躁| 欧美日韩综合久久久久久| 爱豆传媒免费全集在线观看| 99热这里只有是精品在线观看| 99国产精品免费福利视频| 26uuu在线亚洲综合色| 精品卡一卡二卡四卡免费| 乱码一卡2卡4卡精品| 精品国产露脸久久av麻豆| av又黄又爽大尺度在线免费看| 成年人午夜在线观看视频| 一区二区日韩欧美中文字幕 | 下体分泌物呈黄色| 国产在线视频一区二区| 最近最新中文字幕免费大全7| 日本欧美国产在线视频| 80岁老熟妇乱子伦牲交| 视频区图区小说| 我要看黄色一级片免费的| 秋霞伦理黄片| 国产国语露脸激情在线看| 日韩欧美精品免费久久| 18+在线观看网站| videosex国产| 精品久久蜜臀av无| 国产成人欧美| 久久国产精品大桥未久av| 日韩中文字幕视频在线看片| 激情视频va一区二区三区| 看十八女毛片水多多多| √禁漫天堂资源中文www| 51国产日韩欧美| tube8黄色片| 深夜精品福利| 最近最新中文字幕大全免费视频 | 亚洲一区二区三区欧美精品| 99香蕉大伊视频| 国产精品三级大全| 黄色毛片三级朝国网站| 国产成人aa在线观看| 91成人精品电影| 欧美激情国产日韩精品一区| 男女免费视频国产| 街头女战士在线观看网站| 大码成人一级视频| 亚洲欧美日韩卡通动漫| 国产无遮挡羞羞视频在线观看| 少妇被粗大的猛进出69影院 | 国产成人a∨麻豆精品| 一个人免费看片子| 亚洲精品日韩在线中文字幕| 丰满少妇做爰视频| 熟女av电影| 亚洲国产av新网站| 欧美成人午夜精品| 久久国产精品大桥未久av| 少妇的逼好多水| 天堂俺去俺来也www色官网| 精品酒店卫生间| 自拍欧美九色日韩亚洲蝌蚪91| 少妇熟女欧美另类| 狠狠婷婷综合久久久久久88av| 国国产精品蜜臀av免费| 韩国精品一区二区三区 | 国产精品久久久久久精品古装| 国产高清不卡午夜福利| 久久综合国产亚洲精品| 丁香六月天网| 日本色播在线视频| 国产片特级美女逼逼视频| 国产毛片在线视频| 欧美精品国产亚洲| 老熟女久久久| 亚洲激情五月婷婷啪啪| 嫩草影院入口| 国产不卡av网站在线观看| 国产一区二区在线观看日韩| 久久久久国产精品人妻一区二区| 久久久久久久久久久免费av| 成人18禁高潮啪啪吃奶动态图| 免费黄网站久久成人精品| 亚洲欧美精品自产自拍| 少妇猛男粗大的猛烈进出视频| 精品视频人人做人人爽| www.色视频.com| 熟女人妻精品中文字幕| 国产精品久久久久久久电影| 啦啦啦视频在线资源免费观看| 久久狼人影院| 一级毛片 在线播放| 麻豆精品久久久久久蜜桃| 波野结衣二区三区在线| 中国国产av一级| videos熟女内射| 最新的欧美精品一区二区| 18禁在线无遮挡免费观看视频| h视频一区二区三区| 少妇被粗大猛烈的视频| 国产片特级美女逼逼视频| 久久久久国产网址| 狂野欧美激情性xxxx在线观看| 亚洲精品一区蜜桃| 哪个播放器可以免费观看大片| 色婷婷久久久亚洲欧美| 久久久久久久久久人人人人人人| 视频中文字幕在线观看| 久热这里只有精品99| 国产欧美日韩综合在线一区二区| 黄网站色视频无遮挡免费观看| 香蕉丝袜av| 美女xxoo啪啪120秒动态图| 男女午夜视频在线观看 | 99久久人妻综合| 日韩av不卡免费在线播放| 亚洲婷婷狠狠爱综合网| 熟女av电影| 各种免费的搞黄视频| 国产69精品久久久久777片| 日本免费在线观看一区| 高清不卡的av网站| 欧美+日韩+精品| 亚洲精品色激情综合| 成人二区视频| 成人毛片60女人毛片免费| 卡戴珊不雅视频在线播放| 秋霞伦理黄片| 美女国产视频在线观看| 18禁动态无遮挡网站| 精品福利永久在线观看| 国产在视频线精品| 90打野战视频偷拍视频| 色视频在线一区二区三区| 男女无遮挡免费网站观看| 日韩精品有码人妻一区| 妹子高潮喷水视频| 亚洲精品视频女| 性色avwww在线观看| 欧美激情国产日韩精品一区| 欧美少妇被猛烈插入视频| 高清在线视频一区二区三区| 久久精品国产自在天天线| 狂野欧美激情性bbbbbb| 全区人妻精品视频| 国产成人a∨麻豆精品| 久久亚洲国产成人精品v| 爱豆传媒免费全集在线观看| 亚洲少妇的诱惑av| 最新中文字幕久久久久| 超碰97精品在线观看| 成人毛片a级毛片在线播放| 男男h啪啪无遮挡| 国产伦理片在线播放av一区| 国产极品粉嫩免费观看在线| 中文字幕av电影在线播放| 波多野结衣一区麻豆| 99久久综合免费| 天天躁夜夜躁狠狠躁躁| 熟女电影av网| 国产高清国产精品国产三级| 欧美成人精品欧美一级黄| 亚洲欧美精品自产自拍| 一级毛片 在线播放| 韩国av在线不卡| 美女脱内裤让男人舔精品视频| 亚洲精品乱久久久久久| 国内精品宾馆在线| 一个人免费看片子| 欧美日本中文国产一区发布| 伊人亚洲综合成人网| 久久这里只有精品19| 91久久精品国产一区二区三区| 伊人久久国产一区二区| 国产精品不卡视频一区二区| 亚洲成人av在线免费| 香蕉国产在线看| 一边亲一边摸免费视频| 亚洲婷婷狠狠爱综合网| 侵犯人妻中文字幕一二三四区| 亚洲少妇的诱惑av| 十八禁网站网址无遮挡| 成人毛片a级毛片在线播放| 亚洲国产av影院在线观看| 精品第一国产精品| 亚洲国产精品成人久久小说| 熟妇人妻不卡中文字幕| 国产又色又爽无遮挡免| 国产精品久久久久久久电影| 国产国拍精品亚洲av在线观看| 免费观看性生交大片5| 欧美精品亚洲一区二区| 新久久久久国产一级毛片| 亚洲,欧美,日韩| 91午夜精品亚洲一区二区三区| 精品一区二区免费观看| 久久青草综合色| 少妇 在线观看| 夜夜爽夜夜爽视频| 90打野战视频偷拍视频| 国产精品一国产av| 欧美精品av麻豆av| av不卡在线播放| 久久精品国产鲁丝片午夜精品| 99re6热这里在线精品视频| 日韩av在线免费看完整版不卡| 色吧在线观看| 亚洲精品av麻豆狂野| 韩国高清视频一区二区三区| 欧美97在线视频| 欧美日韩一区二区视频在线观看视频在线| 国产精品无大码| 国产片内射在线| 黄色视频在线播放观看不卡| 久久久久久久亚洲中文字幕| 免费高清在线观看日韩| 看非洲黑人一级黄片| 久久综合国产亚洲精品| 久久久久国产精品人妻一区二区| 亚洲三级黄色毛片| av.在线天堂| 亚洲av男天堂| 国产免费一区二区三区四区乱码| 久久久久久久大尺度免费视频| 国产成人一区二区在线| 亚洲一区二区三区欧美精品| 久久久国产一区二区| 国内精品宾馆在线| 精品亚洲成国产av| videosex国产| 视频区图区小说| 国产成人免费无遮挡视频| 久热久热在线精品观看| 宅男免费午夜| 亚洲一区二区三区欧美精品| 色哟哟·www| a 毛片基地| 午夜福利乱码中文字幕| 边亲边吃奶的免费视频| 国产男人的电影天堂91| 亚洲四区av| 欧美精品高潮呻吟av久久| 国产男女超爽视频在线观看| 搡老乐熟女国产| 欧美精品高潮呻吟av久久| 永久免费av网站大全| 亚洲,一卡二卡三卡| 在线观看一区二区三区激情| 亚洲欧美一区二区三区黑人 | 人人妻人人澡人人看| 久久久久精品性色| 免费av不卡在线播放| av天堂久久9| 国产高清三级在线| 大香蕉久久成人网| 国产老妇伦熟女老妇高清| 少妇高潮的动态图| 少妇的逼好多水| 黑人巨大精品欧美一区二区蜜桃 | 久久久久精品性色| 国产免费现黄频在线看| 成人18禁高潮啪啪吃奶动态图| 免费黄网站久久成人精品| 日产精品乱码卡一卡2卡三| 最近手机中文字幕大全| 日韩成人伦理影院| 国产欧美日韩综合在线一区二区| 水蜜桃什么品种好| 久久久久久久久久人人人人人人| 久久综合国产亚洲精品| 国产免费又黄又爽又色| 久久精品国产自在天天线| 国产精品99久久99久久久不卡 | 国产爽快片一区二区三区| 日韩av在线免费看完整版不卡| 春色校园在线视频观看| 91精品三级在线观看| 免费久久久久久久精品成人欧美视频 | 免费黄网站久久成人精品| 亚洲av国产av综合av卡| 最近手机中文字幕大全| 看十八女毛片水多多多| 久久久久人妻精品一区果冻| 国产日韩欧美在线精品| 一区二区av电影网| 黄网站色视频无遮挡免费观看| 午夜福利网站1000一区二区三区| 91精品伊人久久大香线蕉| 中国美白少妇内射xxxbb| 九色成人免费人妻av| 亚洲av在线观看美女高潮| 国产精品一国产av| 51国产日韩欧美| 亚洲第一区二区三区不卡| 卡戴珊不雅视频在线播放| 伦精品一区二区三区| 国产精品欧美亚洲77777| 新久久久久国产一级毛片| 看非洲黑人一级黄片| 黄网站色视频无遮挡免费观看| 边亲边吃奶的免费视频| 性色avwww在线观看| 黄色配什么色好看| 男女国产视频网站| 少妇精品久久久久久久| 少妇人妻 视频| 香蕉精品网在线| 蜜臀久久99精品久久宅男| 国产 精品1| 亚洲第一区二区三区不卡| 热re99久久国产66热| 男女下面插进去视频免费观看 | 国产高清不卡午夜福利| 久久精品久久精品一区二区三区| 久久久欧美国产精品| 欧美精品国产亚洲| 久久久a久久爽久久v久久| 99热国产这里只有精品6| 国产无遮挡羞羞视频在线观看| 久久久亚洲精品成人影院| 又黄又粗又硬又大视频| 免费高清在线观看日韩| 亚洲 欧美一区二区三区| 欧美激情 高清一区二区三区| 少妇的逼好多水| 国产精品麻豆人妻色哟哟久久| 亚洲美女视频黄频| 午夜精品国产一区二区电影| 久久99精品国语久久久| 水蜜桃什么品种好| 高清黄色对白视频在线免费看| 国语对白做爰xxxⅹ性视频网站| 亚洲中文av在线| 国产精品女同一区二区软件| 国产午夜精品一二区理论片| 亚洲欧洲日产国产| 精品一区二区三卡| 爱豆传媒免费全集在线观看| 高清在线视频一区二区三区| 视频区图区小说| 亚洲精品国产av蜜桃| 亚洲美女黄色视频免费看| 麻豆乱淫一区二区| 亚洲精品日韩在线中文字幕| 国产精品一区www在线观看| 满18在线观看网站| 婷婷色综合www| 精品国产露脸久久av麻豆| 亚洲高清免费不卡视频| 全区人妻精品视频| 伊人亚洲综合成人网| 欧美国产精品va在线观看不卡| 99国产综合亚洲精品| 成人无遮挡网站| 大香蕉久久网| 久久精品久久精品一区二区三区| 国产成人午夜福利电影在线观看| 多毛熟女@视频| 午夜免费观看性视频| 午夜激情久久久久久久| 精品亚洲成国产av| 精品久久久久久电影网| 亚洲图色成人| 国产欧美日韩综合在线一区二区| 久热久热在线精品观看| 欧美精品一区二区免费开放| 热re99久久精品国产66热6| 只有这里有精品99| 国产黄色视频一区二区在线观看| 国产福利在线免费观看视频| 亚洲精品国产av成人精品| 久久久久久久久久人人人人人人| 热re99久久国产66热| 国产一区二区激情短视频 | 精品久久国产蜜桃| 日韩视频在线欧美| 精品午夜福利在线看| 国产精品麻豆人妻色哟哟久久| 桃花免费在线播放| 亚洲少妇的诱惑av| 亚洲av.av天堂| 黑人巨大精品欧美一区二区蜜桃 | 极品少妇高潮喷水抽搐| 人人妻人人澡人人看| 久久精品国产鲁丝片午夜精品| 日韩不卡一区二区三区视频在线| 狠狠精品人妻久久久久久综合| 少妇人妻 视频| 乱人伦中国视频| 一级毛片 在线播放| 新久久久久国产一级毛片| 中国三级夫妇交换| 亚洲国产毛片av蜜桃av| 日韩大片免费观看网站| 午夜福利视频在线观看免费| 我的女老师完整版在线观看| 男人舔女人的私密视频| 精品国产一区二区久久| 久久精品国产a三级三级三级| 日韩制服骚丝袜av| 十八禁高潮呻吟视频| 国产精品秋霞免费鲁丝片| 欧美成人午夜免费资源| 免费播放大片免费观看视频在线观看| 国产爽快片一区二区三区| 精品国产一区二区久久| 亚洲国产精品国产精品| 极品人妻少妇av视频| 久久久国产欧美日韩av| 五月玫瑰六月丁香| 国产精品久久久av美女十八| 久久人人爽人人爽人人片va| 精品视频人人做人人爽| 国国产精品蜜臀av免费| 久久毛片免费看一区二区三区| 制服人妻中文乱码| 伦精品一区二区三区| 最近中文字幕高清免费大全6| 五月开心婷婷网| 亚洲 欧美一区二区三区| 伦理电影大哥的女人| 国语对白做爰xxxⅹ性视频网站| 天堂中文最新版在线下载| 曰老女人黄片| 男女啪啪激烈高潮av片| 建设人人有责人人尽责人人享有的| 欧美精品一区二区免费开放| 久久久久久人人人人人| 色婷婷av一区二区三区视频| 丝袜喷水一区| 国产亚洲午夜精品一区二区久久| 国产午夜精品一二区理论片| 天美传媒精品一区二区| 九草在线视频观看| 在线观看免费高清a一片| 大陆偷拍与自拍| 热99国产精品久久久久久7| 午夜福利网站1000一区二区三区| 熟女人妻精品中文字幕| 久久99精品国语久久久| 国产精品蜜桃在线观看| 一区二区av电影网| 国产成人精品一,二区| 九色成人免费人妻av| 免费看不卡的av| 综合色丁香网| 夫妻性生交免费视频一级片| 国产在视频线精品| 亚洲一码二码三码区别大吗| 欧美日韩av久久| 国产精品成人在线| 三上悠亚av全集在线观看| 国产成人免费观看mmmm| 免费人成在线观看视频色| 中国国产av一级| 在线天堂中文资源库| 精品国产国语对白av| 男女免费视频国产| 在线免费观看不下载黄p国产| 欧美日韩成人在线一区二区| 亚洲精品国产av蜜桃| 成人国产av品久久久| 日韩伦理黄色片| 久久久a久久爽久久v久久| 国产精品免费大片| 在线精品无人区一区二区三| 大片免费播放器 马上看| 男女下面插进去视频免费观看 | 国产精品一国产av| 亚洲伊人色综图| 久久久国产一区二区| 熟女av电影| 18禁在线无遮挡免费观看视频| 亚洲国产精品一区二区三区在线| 三级国产精品片| 日韩不卡一区二区三区视频在线| 亚洲国产精品国产精品| 日本wwww免费看| 青春草视频在线免费观看| 免费大片18禁| 美女脱内裤让男人舔精品视频| av女优亚洲男人天堂| 午夜av观看不卡| 老熟女久久久| 日韩av免费高清视频| 97超碰精品成人国产| 国产一级毛片在线| 26uuu在线亚洲综合色| 免费黄频网站在线观看国产| 男女国产视频网站| 国产在线免费精品| 免费日韩欧美在线观看| 一区二区三区精品91| av线在线观看网站| 啦啦啦在线观看免费高清www| 亚洲欧美成人精品一区二区| 99久久综合免费| 亚洲一码二码三码区别大吗| 啦啦啦啦在线视频资源| 日日撸夜夜添| 考比视频在线观看| 国精品久久久久久国模美| 成年美女黄网站色视频大全免费| 亚洲五月色婷婷综合| 久久精品国产自在天天线| 天天影视国产精品| 成人漫画全彩无遮挡| 夫妻性生交免费视频一级片|