• <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期

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

    ( 周口師范學院 生命科學與農學學院, 河南 周口 466001 )

    ?

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

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

    ( 周口師范學院 生命科學與農學學院, 河南 周口 466001 )

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

    關鍵詞:番茄, 晚疫病, 基因, 垂直同源, 進化分析

    番茄(Solanumlycopersicum)是在世界范圍內廣泛種植的主要蔬菜之一,被各國人民所喜愛。番茄晚疫病(Phytophthora infestans)是番茄生產中的主要病害之一,可導致植株莖葉死亡和果實變褐甚至腐爛,每年在番茄主產區(qū)均造成較大的經濟損失(Kamoun ,2001 ),已成為番茄生產的主要障礙之一。由于病源生理小種的不斷變異,導致生產上番茄品種的抗性逐漸喪失,加重了晚疫病病害的發(fā)生(趙統(tǒng)敏等,2006)。利用抗菌劑類藥物進行防治不僅會造成環(huán)境污染,還會使病原菌產生抗耐藥性,促進生理小種變異(Daayf & Platt, 2012),不利于晚疫病的控制,而選育抗性品種不僅可以從根本上克服番茄晚疫病的發(fā)生,而且經濟有效。鑒定和克隆晚疫病抗性基因是選育抗性品種的基礎,尤其一些廣譜主效基因在育種上的應用會產生較好的效果。近年來,在番茄抗晚疫病基因克隆方面取得了較好的進展,如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轉入番茄栽培種中使其獲得了廣譜晚疫病菌抗性,并且對其進行了染色體精細定位。然而目前已有的番茄野生資源中,尚沒有發(fā)現(xiàn)更好的其他抗性基因,尤其廣譜性晚疫病抗性基因。

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

    1材料與方法

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

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

    1.2 蛋白基本信息和基因結構、基因定位信息獲取

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

    1.3 蛋白分子量、等電點和保守結構域分析和亞細胞定位

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

    1.4 蛋白質序列進化樹構建

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

    2結果與分析

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

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

    2.1.2 番茄SLy-Rpi-blb2蛋白家族基因定位和亞細胞定位基因定位結果(表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ù)基因組上相同位點,這與兩條序列來自不同抗性品種相關;第二簇包括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在線軟件進行亞細胞定位結果表明,11條蛋白序列均定位于細胞質中,這與野生馬鈴薯Rpi-blb2蛋白的亞細胞定位結果相同。

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

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

    表 2 番茄SLy-Rpi-blb2蛋白家族基因定位與結構信息

    注: “—”表示未能獲得相關信息。

    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蛋白家族成員的序列比對、基因結構與進化

    為了進一步確定SLy-Rpi-blb2蛋白家族成員之間以及與馬鈴薯Rpi-blb2的同源關系,利用馬鈴薯Rpi-blb2蛋白與番茄蛋白數(shù)據(jù)庫中相似性達到40%以上的18條抗性蛋白(包括晚疫病抗性蛋白)序列構建拓撲結構樹,結果(圖3)表明,野生馬鈴薯Rpi-blb2 蛋白(圖3中紅色箭頭所指)與表1中的11條番茄蛋白序列歸為一大類(圖3:A中蛋白家族I),其他7條序列分為蛋白家族II和蛋白家族Ⅲ。表明這11條序列與馬鈴薯晚疫病抗性蛋白序列Rpi-blb2可能為共同垂直同源蛋白。將11條序列與種子序列進行多序列比對分析(圖4),結果表明,有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個氨基酸序列,這可能是兩者的祖先基因因為部分染色體片段缺失造成的,也可能是堿基突變造成編碼序列提前終止引起的?;蚪Y構(圖3:B)表明,該家族成員有3種基因結構類型,第一類含有2個內含子(5條序列),第二類含有一個較小的內含子(2條序列),第三類不具有內含子結構(XP_004239406),比較這些基因內含子序列發(fā)現(xiàn),第一個內含子序列相對較長,且長度差異較大(556~2 199 bp),序列間相似性差異也較大;第二內含子位于靠近蛋白質N段的編碼區(qū)內,所有同源蛋白相應核酸序列的第二內含子長度均為75 bp,且序列高度保守,這表明第二內含子是一個相當保守的區(qū)域,可能對這類蛋白具有重要作用,因此進化過程中不易發(fā)生改變。

    3討論

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

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

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

    圖 2 番茄SLy-Rpi-blb2蛋白家族成員結構域Fig. 2 Conserved domain for tomato SLy-Rpi-blb2 protein family

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

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

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

    在茄科植物不同物種的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,其相應的蛋白序列均與運輸抑制響應蛋白高度相似,這可能與該簇基因對疾病抗性的響應密切相關。另外2條AAC32252.1和AAC32253.1未能定位到具體的位置,但這2條序列來源于抗根結線蟲栽培品種Motelle,這一品種中的抗性基因均來源于茄科物種Solanumperuvianum中(Smith, 1944),這一結論與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個成員來自于易感根結線蟲品種(Heinz 1706),因此部分成員在基因組被定位在相同位置(圖1),如NP_001234063.1和XP_004240499.1分享同一基因位點,AAC97933.1、NP_001234622.1和ABI96218.1分享同一基因位點。這些蛋白在番茄抗病反應中的功能目前并不完全清楚。Mi-1蛋白序列(AAC97933.1)與馬鈴薯Rpi-blb2具有82%的相似性(van der Vossen et al, 2005),近來被證明是在ETR3蛋白參與下,通過感知荷爾蒙來介導番茄植株對根結線蟲的抗性(Mantelin et al,2013),其同源序列Mi-1.2被證明具有根結線蟲抗性的功能(Nombela et al, 2003),其他同源蛋白由于沒有相應的獨立表現(xiàn)型,因此對其功能的鑒定比較困難(Stuart et al, 2007)。但由于可變剪接產生的較短的抗性蛋白質同源序列,如本研究中的ABI96218.1(Mi-1G)編碼996 個氨基酸,具有完整N段區(qū)域,能反向調節(jié)Mi-1.2(Stuart et al, 2007),從而調節(jié)番茄抗病性。定位于5號染色體上的蛋白質序列XP_004239406.1(566 個氨基酸),與其他10條序列相比,該序列在N段和C段分別缺失了大約500 個氨基酸和100個氨基酸的片段(圖2和圖4),但該序列仍保留了2個保守域的序列特征(圖2),從這一點來看,該序列與圖2中第一類序列應歸為同一類;該基因編碼區(qū)序列是完整的,這意味著該基因(蛋白)可能具有新的功能,在后續(xù)功能研究中應作為重點研究之一,然而該基因與其他同源基因的進化關系目前尚不能完全確定。

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

    參考文獻:

    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在番茄中的表達 [J]. 園藝學報,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]. 生物技術通報,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]. 江蘇農業(yè)學報,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

    文獻標識碼:A

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

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

    基金項目:河南省科技廳基礎與前沿技術研究項目 (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

    黄色 视频免费看| 十八禁网站免费在线| 精品国内亚洲2022精品成人| xxxwww97欧美| 国产精品av久久久久免费| 亚洲国产精品999在线| 亚洲午夜精品一区,二区,三区| 99热这里只有精品一区 | 国产精品一区二区免费欧美| 午夜福利视频1000在线观看| 搡老岳熟女国产| 变态另类成人亚洲欧美熟女| 黄色丝袜av网址大全| 黄片大片在线免费观看| 久久精品影院6| 男人舔女人的私密视频| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品国产亚洲av高清一级| 日韩欧美 国产精品| 国产亚洲精品av在线| 精品久久久久久久久久久久久| 久99久视频精品免费| 国产精品日韩av在线免费观看| 两个人的视频大全免费| 看黄色毛片网站| 伦理电影免费视频| 一个人观看的视频www高清免费观看 | 香蕉国产在线看| 国产欧美日韩精品亚洲av| 国产精品久久久久久亚洲av鲁大| 天天一区二区日本电影三级| 国产高清videossex| а√天堂www在线а√下载| 少妇被粗大的猛进出69影院| 身体一侧抽搐| 国产精品乱码一区二三区的特点| 久久精品国产99精品国产亚洲性色| 久久久久久久精品吃奶| 国产成年人精品一区二区| 国产av在哪里看| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美一级毛片孕妇| 日本一二三区视频观看| 欧美黑人精品巨大| 悠悠久久av| www日本在线高清视频| 99热这里只有精品一区 | 99国产极品粉嫩在线观看| 精品国产亚洲在线| 国产精品亚洲av一区麻豆| 亚洲五月婷婷丁香| 十八禁网站免费在线| 嫩草影视91久久| 国产熟女xx| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品久久国产高清桃花| 欧美日韩中文字幕国产精品一区二区三区| 精品人妻1区二区| 久热爱精品视频在线9| 九色国产91popny在线| 久久这里只有精品19| or卡值多少钱| 黄色女人牲交| 日韩精品青青久久久久久| 日本成人三级电影网站| 日日摸夜夜添夜夜添小说| 国产成人欧美在线观看| 毛片女人毛片| 久久性视频一级片| 日本免费a在线| 色综合亚洲欧美另类图片| 午夜两性在线视频| 日日摸夜夜添夜夜添小说| 午夜成年电影在线免费观看| 两人在一起打扑克的视频| a在线观看视频网站| 久久人妻福利社区极品人妻图片| 他把我摸到了高潮在线观看| 亚洲精品在线观看二区| 日本黄色视频三级网站网址| 成人永久免费在线观看视频| 亚洲 欧美一区二区三区| 看免费av毛片| 午夜福利成人在线免费观看| 在线观看日韩欧美| 免费搜索国产男女视频| 欧美在线一区亚洲| 久久久精品欧美日韩精品| www.999成人在线观看| 久久亚洲精品不卡| 日韩成人在线观看一区二区三区| 一区二区三区激情视频| 日韩精品中文字幕看吧| 日韩欧美在线乱码| 国产99久久九九免费精品| 精品欧美一区二区三区在线| 国产成人一区二区三区免费视频网站| 久久久久九九精品影院| 国产激情偷乱视频一区二区| 中出人妻视频一区二区| 91国产中文字幕| 亚洲aⅴ乱码一区二区在线播放 | 亚洲国产精品久久男人天堂| 桃红色精品国产亚洲av| 亚洲无线在线观看| 亚洲精品美女久久久久99蜜臀| 女人高潮潮喷娇喘18禁视频| 国产亚洲精品av在线| 欧美中文综合在线视频| avwww免费| 亚洲国产精品999在线| 热99re8久久精品国产| a级毛片a级免费在线| 久久精品国产亚洲av香蕉五月| 亚洲色图av天堂| 高清在线国产一区| 国产欧美日韩一区二区精品| 亚洲精品色激情综合| 999精品在线视频| 国产熟女xx| 成人一区二区视频在线观看| 欧美黄色淫秽网站| 制服丝袜大香蕉在线| 男女做爰动态图高潮gif福利片| 精品福利观看| 怎么达到女性高潮| 欧美另类亚洲清纯唯美| 免费人成视频x8x8入口观看| 亚洲成人国产一区在线观看| 久久婷婷成人综合色麻豆| 国产成+人综合+亚洲专区| 男人的好看免费观看在线视频 | 亚洲精华国产精华精| 91麻豆精品激情在线观看国产| 精华霜和精华液先用哪个| 亚洲电影在线观看av| 一区二区三区国产精品乱码| 成年女人毛片免费观看观看9| 在线观看美女被高潮喷水网站 | 午夜两性在线视频| 精品第一国产精品| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品综合一区二区三区| 国产精品av视频在线免费观看| 老司机深夜福利视频在线观看| 成人18禁在线播放| bbb黄色大片| a级毛片在线看网站| 日韩欧美一区二区三区在线观看| 精品国产超薄肉色丝袜足j| 亚洲中文日韩欧美视频| 啦啦啦观看免费观看视频高清| 亚洲aⅴ乱码一区二区在线播放 | 99久久久亚洲精品蜜臀av| 亚洲国产日韩欧美精品在线观看 | 12—13女人毛片做爰片一| 国产av又大| 亚洲性夜色夜夜综合| 日韩大码丰满熟妇| 极品教师在线免费播放| 天天添夜夜摸| 国产亚洲av高清不卡| 老司机福利观看| 国产熟女午夜一区二区三区| 欧美日韩乱码在线| 亚洲精品美女久久av网站| 久久这里只有精品中国| 欧美+亚洲+日韩+国产| 亚洲欧美精品综合一区二区三区| 成人一区二区视频在线观看| 欧美极品一区二区三区四区| 亚洲色图 男人天堂 中文字幕| 国产伦人伦偷精品视频| 亚洲精品粉嫩美女一区| 黄色女人牲交| 国产一区二区三区视频了| 又紧又爽又黄一区二区| 亚洲最大成人中文| 神马国产精品三级电影在线观看 | av免费在线观看网站| 精品一区二区三区四区五区乱码| 桃红色精品国产亚洲av| 亚洲性夜色夜夜综合| 久久久久久久午夜电影| 久久99热这里只有精品18| av中文乱码字幕在线| 亚洲精品一卡2卡三卡4卡5卡| 黑人操中国人逼视频| 久久久久久久久免费视频了| 97超级碰碰碰精品色视频在线观看| 久9热在线精品视频| 一个人观看的视频www高清免费观看 | 日本精品一区二区三区蜜桃| 香蕉av资源在线| 精品电影一区二区在线| 变态另类成人亚洲欧美熟女| 丁香六月欧美| 丰满人妻熟妇乱又伦精品不卡| 精品久久久久久久末码| 亚洲国产精品成人综合色| www.999成人在线观看| 精品久久久久久,| 一进一出抽搐gif免费好疼| 久热爱精品视频在线9| 国产午夜精品论理片| 国产激情偷乱视频一区二区| 高潮久久久久久久久久久不卡| 免费看日本二区| 97超级碰碰碰精品色视频在线观看| 成人av在线播放网站| 免费看十八禁软件| 日韩av在线大香蕉| 51午夜福利影视在线观看| 欧美一区二区国产精品久久精品 | 久久天堂一区二区三区四区| av中文乱码字幕在线| 久久久久久人人人人人| 欧美三级亚洲精品| 亚洲全国av大片| 老汉色∧v一级毛片| 全区人妻精品视频| 变态另类成人亚洲欧美熟女| 久久精品亚洲精品国产色婷小说| 精品久久久久久久人妻蜜臀av| 亚洲欧美日韩高清在线视频| 国产激情偷乱视频一区二区| 一级毛片精品| 成人精品一区二区免费| 高潮久久久久久久久久久不卡| 99热这里只有是精品50| 91老司机精品| 日日夜夜操网爽| 亚洲国产精品久久男人天堂| 欧美极品一区二区三区四区| 亚洲成a人片在线一区二区| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产精品久久男人天堂| 老司机午夜福利在线观看视频| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久久精品吃奶| 欧美日韩亚洲综合一区二区三区_| 亚洲中文av在线| 午夜福利免费观看在线| a级毛片a级免费在线| 亚洲精品色激情综合| 宅男免费午夜| 宅男免费午夜| 一进一出好大好爽视频| 亚洲成人中文字幕在线播放| 国产男靠女视频免费网站| 老熟妇仑乱视频hdxx| 老司机靠b影院| 欧美成人午夜精品| 国产免费男女视频| 特级一级黄色大片| 亚洲一区二区三区不卡视频| 中文字幕高清在线视频| 亚洲欧美日韩高清在线视频| 中文在线观看免费www的网站 | 午夜精品久久久久久毛片777| 国产亚洲av嫩草精品影院| 欧美日韩一级在线毛片| 亚洲成av人片免费观看| 在线观看免费日韩欧美大片| 国产成人精品无人区| 岛国视频午夜一区免费看| 国产一级毛片七仙女欲春2| 色哟哟哟哟哟哟| 在线视频色国产色| 久久中文看片网| 成人午夜高清在线视频| 麻豆国产av国片精品| 法律面前人人平等表现在哪些方面| 黑人欧美特级aaaaaa片| 99国产精品99久久久久| 国产一区二区三区视频了| 久久久国产欧美日韩av| 88av欧美| 国产精品久久久久久人妻精品电影| 亚洲成av人片免费观看| 成人三级做爰电影| 99久久99久久久精品蜜桃| 国产精品1区2区在线观看.| 波多野结衣高清无吗| 91麻豆精品激情在线观看国产| 在线观看美女被高潮喷水网站 | 国产一区二区在线观看日韩 | 五月玫瑰六月丁香| 欧美成人性av电影在线观看| 18禁黄网站禁片午夜丰满| 成人18禁在线播放| 高潮久久久久久久久久久不卡| 国产精品一区二区三区四区免费观看 | 色老头精品视频在线观看| 白带黄色成豆腐渣| 在线观看免费视频日本深夜| 欧美日韩黄片免| av天堂在线播放| 免费一级毛片在线播放高清视频| 国产成人一区二区三区免费视频网站| 黄片大片在线免费观看| 久久香蕉激情| 很黄的视频免费| 夜夜夜夜夜久久久久| 99热这里只有精品一区 | 757午夜福利合集在线观看| 一二三四在线观看免费中文在| 999久久久精品免费观看国产| 我的老师免费观看完整版| 波多野结衣高清无吗| 精华霜和精华液先用哪个| 免费电影在线观看免费观看| av有码第一页| cao死你这个sao货| 国产精华一区二区三区| 日本一区二区免费在线视频| 禁无遮挡网站| 禁无遮挡网站| 欧美3d第一页| 国产成人啪精品午夜网站| √禁漫天堂资源中文www| svipshipincom国产片| 99热这里只有是精品50| 黑人操中国人逼视频| 成人国产综合亚洲| 国产爱豆传媒在线观看 | 国产精品亚洲一级av第二区| 精品国产超薄肉色丝袜足j| www.自偷自拍.com| 亚洲成人中文字幕在线播放| 日韩欧美三级三区| 久久久久久久午夜电影| 久久久久久大精品| 亚洲精品粉嫩美女一区| 夜夜爽天天搞| 久久精品成人免费网站| 精品熟女少妇八av免费久了| 午夜福利成人在线免费观看| 中文字幕久久专区| 又黄又爽又免费观看的视频| 亚洲一区二区三区色噜噜| 国产在线观看jvid| 日本 av在线| 美女午夜性视频免费| 变态另类丝袜制服| 狠狠狠狠99中文字幕| 毛片女人毛片| 国产亚洲欧美98| 午夜精品一区二区三区免费看| 曰老女人黄片| 91字幕亚洲| 在线永久观看黄色视频| 99国产综合亚洲精品| 99久久无色码亚洲精品果冻| 亚洲欧美日韩高清专用| 亚洲免费av在线视频| 在线国产一区二区在线| 亚洲欧美精品综合一区二区三区| 亚洲专区中文字幕在线| 国产精品免费一区二区三区在线| 一个人免费在线观看电影 | cao死你这个sao货| 久久久国产成人免费| 老司机午夜福利在线观看视频| 午夜精品一区二区三区免费看| 欧美性长视频在线观看| 少妇粗大呻吟视频| 久久亚洲精品不卡| 亚洲中文av在线| 一边摸一边抽搐一进一小说| 国产成人av教育| 色尼玛亚洲综合影院| 午夜精品一区二区三区免费看| 99久久99久久久精品蜜桃| 亚洲色图av天堂| 国产精品综合久久久久久久免费| 国产成人系列免费观看| 国产午夜精品久久久久久| 99久久无色码亚洲精品果冻| 国产区一区二久久| 亚洲在线自拍视频| 国产精品一区二区三区四区免费观看 | 悠悠久久av| 午夜福利欧美成人| 午夜a级毛片| or卡值多少钱| 国产亚洲av嫩草精品影院| 亚洲一区高清亚洲精品| 久久久精品大字幕| 久热爱精品视频在线9| 午夜免费观看网址| av视频在线观看入口| 国产精品美女特级片免费视频播放器 | 日本在线视频免费播放| 搡老熟女国产l中国老女人| 国产视频一区二区在线看| 此物有八面人人有两片| 国产精品一及| 床上黄色一级片| 国产成人影院久久av| 久久久久精品国产欧美久久久| 成人精品一区二区免费| av免费在线观看网站| 亚洲欧美日韩无卡精品| 哪里可以看免费的av片| 国产精品 国内视频| 国产av又大| 欧美成人免费av一区二区三区| 精品国产乱子伦一区二区三区| 变态另类成人亚洲欧美熟女| 精品人妻1区二区| av有码第一页| 久久九九热精品免费| 日日摸夜夜添夜夜添小说| 午夜免费成人在线视频| 亚洲精品美女久久av网站| cao死你这个sao货| 99久久精品国产亚洲精品| 91老司机精品| 免费在线观看视频国产中文字幕亚洲| 欧美成人午夜精品| 精品人妻1区二区| 狠狠狠狠99中文字幕| 后天国语完整版免费观看| 精品久久久久久成人av| 欧美最黄视频在线播放免费| 在线免费观看的www视频| 欧美成狂野欧美在线观看| 国产伦在线观看视频一区| 午夜免费成人在线视频| 男女之事视频高清在线观看| 国产高清有码在线观看视频 | 亚洲人与动物交配视频| 午夜精品一区二区三区免费看| 国产精品一区二区精品视频观看| 嫩草影视91久久| 久久精品aⅴ一区二区三区四区| 亚洲国产高清在线一区二区三| 亚洲真实伦在线观看| 中文字幕人成人乱码亚洲影| 一二三四社区在线视频社区8| 一个人免费在线观看电影 | 色综合亚洲欧美另类图片| 亚洲无线在线观看| 色av中文字幕| 五月玫瑰六月丁香| 深夜精品福利| 国产日本99.免费观看| 欧美一级a爱片免费观看看 | 国产精品久久久av美女十八| 国内毛片毛片毛片毛片毛片| 制服诱惑二区| 老司机福利观看| 老司机深夜福利视频在线观看| 久久精品人妻少妇| 露出奶头的视频| 1024视频免费在线观看| 午夜老司机福利片| 中文字幕久久专区| 我的老师免费观看完整版| av在线播放免费不卡| 中文字幕人成人乱码亚洲影| 亚洲av美国av| 成在线人永久免费视频| 久久久久性生活片| 欧美一级a爱片免费观看看 | 18美女黄网站色大片免费观看| 国产探花在线观看一区二区| www.自偷自拍.com| 少妇被粗大的猛进出69影院| av超薄肉色丝袜交足视频| 国产精品永久免费网站| 国产黄片美女视频| 久久婷婷人人爽人人干人人爱| 叶爱在线成人免费视频播放| 午夜福利视频1000在线观看| 久久久精品国产亚洲av高清涩受| 成人欧美大片| 在线观看www视频免费| 午夜福利欧美成人| 大型av网站在线播放| 国产激情久久老熟女| 亚洲精品一卡2卡三卡4卡5卡| 香蕉av资源在线| 国产精品 国内视频| 女人被狂操c到高潮| 看片在线看免费视频| 变态另类成人亚洲欧美熟女| 国产av在哪里看| 亚洲狠狠婷婷综合久久图片| 黄色片一级片一级黄色片| 精品第一国产精品| 一区二区三区激情视频| 十八禁网站免费在线| 午夜成年电影在线免费观看| 一进一出抽搐gif免费好疼| 51午夜福利影视在线观看| 91老司机精品| 欧美在线黄色| 午夜老司机福利片| 久久香蕉精品热| 黄频高清免费视频| 久久国产精品人妻蜜桃| 日韩欧美精品v在线| 久9热在线精品视频| 女人爽到高潮嗷嗷叫在线视频| 久久久久久久久久黄片| 91麻豆av在线| 亚洲专区中文字幕在线| 久久精品91无色码中文字幕| 国产精品国产高清国产av| 久久久水蜜桃国产精品网| 校园春色视频在线观看| 久久婷婷成人综合色麻豆| 黄色视频,在线免费观看| 亚洲天堂国产精品一区在线| 欧美av亚洲av综合av国产av| 这个男人来自地球电影免费观看| 在线观看舔阴道视频| 1024香蕉在线观看| 少妇被粗大的猛进出69影院| av在线天堂中文字幕| 亚洲人与动物交配视频| 日韩大尺度精品在线看网址| 国产一级毛片七仙女欲春2| 午夜免费观看网址| 1024手机看黄色片| 中文字幕高清在线视频| 人人妻,人人澡人人爽秒播| 午夜亚洲福利在线播放| 欧美3d第一页| 老熟妇乱子伦视频在线观看| 久久亚洲精品不卡| www.www免费av| 无限看片的www在线观看| 人人妻人人澡欧美一区二区| 日韩三级视频一区二区三区| 很黄的视频免费| 少妇熟女aⅴ在线视频| av天堂在线播放| 国产精品精品国产色婷婷| 99热这里只有精品一区 | 午夜两性在线视频| 亚洲狠狠婷婷综合久久图片| 免费观看精品视频网站| 神马国产精品三级电影在线观看 | www.www免费av| 国产激情偷乱视频一区二区| 久久久久久九九精品二区国产 | 欧美丝袜亚洲另类 | 一边摸一边抽搐一进一小说| 国产精品1区2区在线观看.| av在线播放免费不卡| 亚洲第一欧美日韩一区二区三区| 人人妻,人人澡人人爽秒播| 国产成人aa在线观看| 两个人的视频大全免费| 午夜免费激情av| 999精品在线视频| 国产真实乱freesex| 18禁观看日本| 美女 人体艺术 gogo| 男女午夜视频在线观看| 97超级碰碰碰精品色视频在线观看| 亚洲全国av大片| 亚洲人与动物交配视频| 欧美又色又爽又黄视频| 久99久视频精品免费| 每晚都被弄得嗷嗷叫到高潮| 在线观看免费日韩欧美大片| 欧美成人午夜精品| 久久久久国内视频| 黄色成人免费大全| 夜夜爽天天搞| 国产精品1区2区在线观看.| 一区二区三区国产精品乱码| 露出奶头的视频| 搡老熟女国产l中国老女人| 我要搜黄色片| 国语自产精品视频在线第100页| 色播亚洲综合网| 欧美不卡视频在线免费观看 | 日韩精品中文字幕看吧| 99久久无色码亚洲精品果冻| 黑人巨大精品欧美一区二区mp4| a在线观看视频网站| 老司机靠b影院| 午夜老司机福利片| 国产男靠女视频免费网站| 国产v大片淫在线免费观看| 在线观看一区二区三区| 国产97色在线日韩免费| 亚洲 欧美 日韩 在线 免费| 久久九九热精品免费| 午夜福利高清视频| 日本撒尿小便嘘嘘汇集6| 亚洲美女黄片视频| 精品久久蜜臀av无| 国产主播在线观看一区二区| 日本在线视频免费播放| 久久中文看片网| 看免费av毛片| 90打野战视频偷拍视频| 欧美日韩国产亚洲二区| 国产成人精品无人区| 国产亚洲av高清不卡| 亚洲成av人片免费观看| 亚洲专区中文字幕在线| 国内精品久久久久久久电影| 久久国产乱子伦精品免费另类| 成人亚洲精品av一区二区| 国产精品久久久人人做人人爽| 精品久久久久久久久久免费视频| 黄色视频不卡| 国产亚洲精品综合一区在线观看 |