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

    栽培番茄與野生番茄NBS-LRR類抗病基因家族的全基因組鑒定及表達(dá)分析

    2021-09-08 00:00史建磊熊自立蘇世聞張海利宰文珊
    關(guān)鍵詞:結(jié)構(gòu)域抗病氨基酸

    史建磊 熊自立 蘇世聞 張海利 宰文珊

    摘要:【目的】對栽培番茄和野生番茄NBS-LRR類抗病基因家族進(jìn)行全基因組鑒定及表達(dá)分析,為NBS-LRR類抗病基因功能研究及其在植物抗病育種中的應(yīng)用提供理論參考?!痉椒ā繎?yīng)用生物信息學(xué)和比較基因組學(xué)方法,從栽培番茄和野生番茄[潘那利番茄(Solanum pennellii)和醋栗番茄(Solanum pimpinellifolium)]基因組中鑒定出NBS-LRR類抗病基因家族成員,明確其類型,并進(jìn)行序列特征、染色體定位、系統(tǒng)進(jìn)化及生物脅迫下表達(dá)模式分析?!窘Y(jié)果】從栽培番茄、潘那利番茄和醋栗番茄基因組中分別鑒定獲得238、202和217個(gè)NBS-LRR類抗病基因家族成員,根據(jù)這些抗病基因編碼的結(jié)構(gòu)域差異,將其分為TNL(TIR-NBS-LRR)、CNL(CC-NBS-LRR)、RNL(RPW8-NBS-LRR)、TN(TIR-NBS)、CN(CC-NBS)、RN(RPW8-NBS)、NL(NBS-LRR)和N(NBS)8種類型,其中TNL、CNL和RNL型基因分別為58、87和13個(gè)。番茄NBS-LRR類抗病蛋白整體呈弱酸性,具有不穩(wěn)定性和親水性,以絲氨酸(Ser)磷酸化為主,二級結(jié)構(gòu)以α-螺旋和無規(guī)則卷曲為主,定位于細(xì)胞核、細(xì)胞質(zhì)和葉綠體。栽培番茄和野生潘那利番茄的NBS-LRR類抗病基因不均勻地散布在13條染色體(0號虛擬染色體~12號染色體)上,分別有58.47%和52.48%的基因形成多個(gè)基因簇,有35.59%和22.77%的基因形成串聯(lián)重復(fù)。栽培番茄和2個(gè)野生番茄共477個(gè)NBS-LRR類抗病蛋白可聚為十大類群,其中N端為TIR的蛋白基本聚在一起,但N端為CC和RPW8結(jié)構(gòu)域的蛋白未能很好分開,聚在多個(gè)類群中;NL和N型基因散布于各類群中。在477個(gè)番茄NBS-LRR類抗病基因中共發(fā)現(xiàn)115對直系同源基因和8對旁系同源基因,主要分布在4號和5號染色體上,其中51.57%的NBS-LRR類抗病基因以同源基因形式存在,大多數(shù)Ka/Ks均小于1,說明其進(jìn)化中主要受到純化選擇?;谵D(zhuǎn)錄組測序數(shù)據(jù)(RNA-Seq)分析結(jié)果表明,多數(shù)NBS-LRR類抗病基因能響應(yīng)不同細(xì)菌脅迫?!窘Y(jié)論】番茄NBS-LRR類抗病基因具有成簇存在的特點(diǎn),在進(jìn)化過程中經(jīng)重復(fù)擴(kuò)張,已形成大量同源基因,且能響應(yīng)多種細(xì)菌脅迫。

    關(guān)鍵詞: 番茄;NBS-LRR類抗病基因;序列特征;染色體定位;系統(tǒng)進(jìn)化;生物脅迫;基因簇;基因重復(fù)

    中圖分類號: S641.203.6? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼: A 文章編號:2095-1191(2021)05-1158-09

    Abstract:【Objective】Genome-wide identification and expression analysis of NBS-LRR disease-resistantgenes in cultivated tomatoand its wild relatives were performed,to provide a theoretical basis for the research on the function of NBS-LRR disease-resistant genes and the application in plant disease resistance breeding. 【Method】NBS-LRR disease-resistant gene family members were identified and analysed on classification, sequence characteristics, chromosome location, phylogeny, and expression patterns under biotic stresses from cultivated(Solanum lycopersicum) and two wild(S. pennellii and S. pimpinellifolium) tomato species genomes using bioinformatics and comparative genomics methods. 【Result】In total,238,202 and 217 NBS-LRR disease-resistant genes were identified in cultivated tomato and S. pennellii and S. pim-pinellifolium genomes,respectively. Based on the differences in the domains encoded by these disease-resistant genes, they were divided into eight types,which were TNL(TIR-NBS-LRR), CNL(CC-NBS-LRR), RNL(RPW8-NBS-LRR), TN(TIR-NBS), CN(CC-NBS), RN(RPW8-NBS), NL(NBS-LRR) and N(NBS),with 58 TNL,87 CNL and 13 RNL genes,respectively. NBS-LRR disease-resistant genes were weakly acidic in general, and were unstable and hydrophilic,with serine(Ser) phosphorylation and secondary structure of alpha-helix and random coil,which were mainly located in the nucleus,cytoplasm and chloroplast. Chromosome localization revealed that the NBS-LRR disease-resistant genes in cultiva-ted tomato and S. pennellii genomes were unevenly distributed on 13 different chromosomes(virtual chromosome 0-chromosome12),with 58.47% and 52.48% genes forming multiple gene clusters,35.59% and 22.77% genes forming tandem duplicates,respectively. Systematic clustering divided 477 NBS-LRR disease-resistant genes in the three genomes into ten groups,the N-terminal TIR-type genes were basically clustered together,but CC-and RPW8-type genes were not well separated and clustered in multiple groups; NL- and N-type genes were scattered in various groups. Gene homology analysis revealed 115 pairs of orthologous genes and 8 pairs of paralogous genes in 477 NBS-LRR disease-resistantgenes, equivalent to 51.57% of NBS-LRR genes had a homologous gene,and that they were mainly distributed on chromosomes 4 and 5, with most Ka/Ks values being less than 1,indicating that purification selection had occurred in evolution. It revealed that the most NBS-LRR disease-resistant genes could respond to different bacterial stresses by transcriptome sequencing(RNA-seq)data analysis. 【Conclusion】The tomato NBS-LRR disease-resistant genes exist in cluster. During the evolution,a large number of homologous genes have appeared through gene expansion. Also,they can respond to a variety of bacterial stresses.

    Key words:tomato; NBS-LRR disease-resistant gene;sequence characteristics; chromosomal localization;phyletic evolution;biotic stress; gene cluster;gene duplication

    Foundation item: Major Science and Technology Project for Breeding New Agricultural(Vegetable) Varieties in Zhejiang(2016C02051-1-2); Wenzhou Basic Scientific Research Project(N20180003); General Scientific Research Pro-ject of Education Department in Zhejiang(Y201840413)

    0 引言

    【研究意義】植物在自然進(jìn)化過程中面對病原微生物的侵染已形成一整套識別與防御機(jī)制,在此過程中抗病基因發(fā)揮著非常重要的作用。NBS-LRR是植物中最大抗病基因家族之一,目前已克隆的抗病基因多為NBS-LRR類抗病基因家族成員(Dangl and Jones,2001)。NBS-LRR類抗病蛋白具有核苷酸結(jié)合位點(diǎn)(Nucleotide-binding site,NBS)和富亮氨酸重復(fù)(Leucine-rich repeat,LRR)典型特征結(jié)構(gòu)域(劉云飛等,2014a,2014b),通過與病原效應(yīng)子結(jié)合,觸發(fā)ETI(Effector-triggered immunity)免疫反應(yīng),從而抑制病原在植物細(xì)胞中的定殖與擴(kuò)展。番茄(Solanum lycopersicum)是世界主要蔬菜作物和科研模式植物之一,同時(shí)其近緣野生種中含豐富的抗病蟲和抗逆資源,鑒定與評價(jià)番茄NBS-LRR類抗病基因,對研究其數(shù)目特征、遺傳進(jìn)化和生物功能具有重要理論意義;同時(shí)在相關(guān)抗病基因深入挖掘的基礎(chǔ)上,對番茄栽培品種性狀改良和優(yōu)良新品種選育也具有積極實(shí)踐意義。【前人研究進(jìn)展】隨著越來越多植物基因組測序的完成,水稻(Oryza sativa)(Bai et al.,2002)、擬南芥(Arabidopsis thaliana)(Meyers et al.,2003)、苜蓿(Medicago sativa)(Ameline-Torregrosa et al.,2008)、馬鈴薯(S. tuberosum)(Lozano et al.,2012)、玉米(Zea mays)(Cheng et al.,2012)、番茄(劉云飛等,2014b)、獼猴桃(Actinidia chinensis)(文歡和劉永立,2015)、大麥(Hordeum vulgare)(Habachi-Houimli et al.,2018)、蝴蝶蘭(Phalaenopsis aphrodite)(蔣卉等,2018)、桑樹(Morus alba)(劉潮等,2019)、谷子(李任建等,2020)等多種植物NBS-LRR類抗病基因已被鑒定。研究結(jié)果顯示,NBS-LRR類抗病蛋白主要包含由氨基端(N端)到羧基端(C端)的TIR/CC、NBS和LRR 3個(gè)結(jié)構(gòu)域(Meyers et al.,2003),有的還含有RPW8等結(jié)構(gòu)域(Xiao et al.,2001)。其中,NBS又稱NB-ARC,是最保守的結(jié)構(gòu)域,由8個(gè)保守基序組成(Meyers et al.,2003),通過結(jié)合水解ATP和GTP并進(jìn)行信號轉(zhuǎn)導(dǎo),調(diào)節(jié)隨后的病原物識別(Fenyk et al,2012);LRR由多個(gè)亮氨酸或脯氨酸和天冬氨酸組成(Lozano et al.,2012),是最主要的病原物特異識別區(qū),通過蛋白互作識別病原物并決定抗病基因的?;院涂共∽V(Bai et al.,2002)。根據(jù)N端是否存在TIR(Toll/interleukin-1 receptor like)結(jié)構(gòu)域,可將NBS-LRR類抗病基因分為TNL(TIR-NBS-LRR)和nTNL(non-TIR-NBS-LRR)兩大類,后者通常會編碼卷曲螺旋(Coiled-coil)結(jié)構(gòu),也稱為CNL(CC-NBS-LRR)(Meyers et al.,2003)。但并非所有NBS-LRR類抗病蛋白均具有3個(gè)結(jié)構(gòu)域,通常會缺少N端或C端結(jié)構(gòu)域?!颈狙芯壳腥朦c(diǎn)】目前番茄全基因組測序已完成(Sato et al.,2012),為在全基因組水平上鑒定NBS-LRR類抗病基因提供了基礎(chǔ)條件。目前盡管已有番茄NBS-LRR類抗病基因家族的全基因組鑒定,但鮮見有關(guān)野生番茄NBS-LRR類抗病基因家族的研究報(bào)道?!緮M解決的關(guān)鍵問題】應(yīng)用生物信息學(xué)和比較基因組學(xué)方法,從栽培番茄和野生番茄[潘那利番茄(S. pennellii)和醋栗番茄(S. pimpinellifolium)]基因組鑒定出NBS-LRR類抗病基因家族成員,明確其類型,并進(jìn)行序列特征、染色體定位、系統(tǒng)進(jìn)化及生物脅迫下表達(dá)模式分析,為番茄相關(guān)基因功能驗(yàn)證及其在抗病育種中的應(yīng)用提供理論依據(jù)。

    1 材料與方法

    1. 1 番茄NBS-LRR類抗病基因家族成員鑒定

    從茄科基因組網(wǎng)站SGN(http://solgenomics.net/)下載番茄基因組序列及其注釋信息,并利用BioEdit 7.2.5構(gòu)建本地?cái)?shù)據(jù)庫。在Pfam(http://pfam.xfam.org/)中下載NB-ARC(PF00931)氨基酸序列,并以此為靶序列進(jìn)行BLASTp(E-value<1)比對檢索,獲得候選基因。然后通過Pfam、SMART和COILS Server鑒定TIR、CC、RPW8、NBS和LRR等結(jié)構(gòu)域,將不含NBS結(jié)構(gòu)域編碼序列的候選基因排除。

    1. 2 番茄NBS-LRR類抗病基因序列特征分析

    利用ExPASy在線網(wǎng)站的Compute pI/MW工具計(jì)算蛋白分子量和理論等電點(diǎn)(pI);以ProtParam在線網(wǎng)站進(jìn)行蛋白理化性質(zhì)預(yù)測;利用ProtScale分析蛋白質(zhì)親/疏水性;運(yùn)用NetPhos 3.1 Server在線程序分析蛋白磷酸化位點(diǎn);采用ExPASy在線網(wǎng)站的SOPMA預(yù)測蛋白二級結(jié)構(gòu);利用WoLF PSORT在線工具進(jìn)行蛋白亞細(xì)胞定位預(yù)測;并以WebLogo 3在線工具繪制保守序列Logo圖。

    1. 3 番茄NBS-LRR類抗病基因定位

    結(jié)合SGN網(wǎng)站中的基因注釋信息,利用MapDraw 2.1進(jìn)行基因染色體定位并繪制其分布圖。隨后對番茄NBS-LRR類抗病基因進(jìn)行基因簇和串聯(lián)重復(fù)分析,劃分基因簇的原則:(1)2個(gè)相鄰基因間距小于200 kb;(2)2個(gè)相鄰基因間其他基因不多于8個(gè)(Meyers et al.,2003;Yang et al.,2008);判斷串聯(lián)重復(fù)的原則:(1)相鄰基因距離小于100 kb;(2)基因間相似度大于70%(Huang et al.,2012)。采用在線程序PAL2NAL(http://www.bork.embl.de/pal2nal/index.cgi?example=Yes#RunP2N)計(jì)算基因間的選擇壓力。

    1. 4 番茄NBS-LRR類抗病基因系統(tǒng)進(jìn)化分析

    首先提取出NBS結(jié)構(gòu)域的氨基酸序列,利用ClustalX 1.83進(jìn)行多序列比對,再以MEGA 7.0的鄰接法(Neighbor-joining,NJ)(參數(shù)設(shè)置:Bootstrap method 1000,Poisson model,and Pairwise deletion)構(gòu)建系統(tǒng)發(fā)育進(jìn)化樹,去除Bootstrap支持率低于50%的節(jié)點(diǎn)。采用DNAMAN 7.0計(jì)算序列間的相似性。

    1. 5 番茄NBS-LRR類抗病基因表達(dá)分析

    從番茄功能基因組數(shù)據(jù)庫TFGD(http://ted.bti.cornell.edu/cgi-bin/TFGD/digital/home.cgi)下載番茄轉(zhuǎn)錄組測序(RNA-Seq)數(shù)據(jù)(D007),利用MeV 4.9.0繪制基因的表達(dá)熱圖,并分析基因的表達(dá)模式。

    2 結(jié)果與分析

    2. 1 番茄NBS-LRR類抗病基因家族成員鑒定結(jié)果

    從栽培番茄及野生潘那利和醋栗番茄的基因組序列中分別檢索獲得238、202和217個(gè)NBS-LRR類抗病基因。根據(jù)這些抗病基因編碼蛋白的結(jié)構(gòu)域差異,將其分為TNL(TIR-NBS-LRR)、CNL(CC-NBS-LRR)、RNL(RPW8-NBS-LRR)、TN(TIR-NBS)、CN(CC-NBS)、RN(RPW8-NBS)、NL(NBS-LRR)和N(NBS)8種類型(表1)。其中,N型基因最多,為226個(gè),占NBS-LRR類抗病基因總數(shù)(657個(gè))的34.4%;TNL、CNL和RNL型基因分別為58、87和13個(gè),共計(jì)158個(gè)。

    2. 2 番茄NBS-LRR類抗病蛋白理化性質(zhì)、二級結(jié)構(gòu)及亞細(xì)胞定位結(jié)果

    158個(gè)番茄TNL、CNL和RNL型NBS-LRR類抗病蛋白的理化性質(zhì)、二級結(jié)構(gòu)及亞細(xì)胞定位如表2所示。栽培番茄TNL型蛋白的氨基酸數(shù)為856~2871個(gè),分子量為98.40~325.87 kD,pI為5.43~8.77,不穩(wěn)定系數(shù)為37.88~51.45,磷酸化氨基酸平均為136個(gè),其中絲氨酸(Ser)、蘇氨酸(Thr)和酪氨酸(Tyr)分別占比67.11%、22.58%和10.30%;nTNL型蛋白的氨基酸數(shù)為696~2613個(gè),分子量為79.05~298.53 kD,pI為4.85~9.03,不穩(wěn)定系數(shù)為37.93~50.34,磷酸化氨基酸數(shù)平均為100個(gè),其中Ser、Thr和Tyr分別占比64.69%、27.07%和8.24%。潘那利番茄TNL型蛋白的氨基酸數(shù)為673~2871個(gè),分子量為76.33~326.23 kD,pI為5.48~8.96,不穩(wěn)定系數(shù)為37.25~52.73,磷酸化氨基酸數(shù)平均為118個(gè),其中Ser、Thr和Tyr分別占比66.87%、22.78%和10.34%;nTNL型蛋白的氨基酸數(shù)為617~1687個(gè),分子量為69.79~194.21 kD,pI為5.34~8.62,不穩(wěn)定系數(shù)為39.82~52.41,磷酸化氨基酸數(shù)平均為101個(gè),其中Ser、Thr和Tyr分別占比67.42%、24.80%和7.78%。醋栗番茄TNL型蛋白的氨基酸數(shù)為973~2871個(gè),分子量為111.40~325.79 kD,pI為5.42~9.23,不穩(wěn)定系數(shù)為40.45~51.45,磷酸化氨基酸數(shù)平均為128個(gè),其中Ser、Thr和Tyr分別占比67.02%、22.86%和10.11%;nTNL型蛋白的氨基酸數(shù)為645~2595個(gè),分子量為73.57~296.37 kD,pI為4.85~8.91,不穩(wěn)定系數(shù)為38.96~49.79,磷酸化氨基酸數(shù)平均為95個(gè),其中Ser、Thr和Tyr分別占比64.29%、27.14%和8.56%。番茄NBS-LRR類抗病蛋白整體上呈弱酸性,且具有不穩(wěn)定性和親水性,以Ser磷酸化為主,其二級結(jié)構(gòu)以α-螺旋和無規(guī)則卷曲為主,主要定位于細(xì)胞核、細(xì)胞質(zhì)和葉綠體。

    對番茄TNL、CNL和RNL型NBS-LRR類抗病蛋白的保守結(jié)構(gòu)域進(jìn)行分析,繪制得到相應(yīng)氨基酸序列的Logo圖(圖1-A、圖1-B和圖1-C)。保守氨基酸組成的NBS結(jié)構(gòu)域P-loop、Kinase2、RNBS-B和GLPL,TNL和nTNL型NBS-LRR類抗病基因不完全一致,如TNL型基因NBS保守結(jié)構(gòu)域的Kinase2一般以天冬氨酸(Asp)結(jié)尾,而nTNL型一般以色氨酸(Trp)結(jié)尾。說明番茄NBS-LRR類抗病基因在進(jìn)化過程中具有較高的保守性。

    2. 3 番茄NBS-LRR類抗病基因染色體定位結(jié)果

    通過獲取番茄NBS-LRR類抗病基因的注釋信息(醋栗番茄基因組注釋不完整,故未作分析),分別將栽培番茄的236個(gè)(另外2個(gè)基因不能確定染色體位置)和野生潘那利番茄的202個(gè)NBS-LRR類抗病基因分別不均勻地散布于12條染色體和0號虛擬染色體上(圖2-A和圖2-B)。番茄NBS-LRR類抗病基因集中在染色體端部,且具有成簇分布的特點(diǎn)。栽培番茄各染色體的基因數(shù)為1(0號虛擬染色體)~45(11號染色體)個(gè),平均為18.15個(gè),其中138個(gè)(58.47%)基因形成47個(gè)基因簇,平均每簇有2.94個(gè)基因,4號染色體10簇37個(gè)基因,最大簇7個(gè)基因,6號和4號染色體成簇基因均超過80.00%,其他染色體成簇基因數(shù)為0~16個(gè);同時(shí),84個(gè)(35.59%)基因?yàn)榇?lián)重復(fù)基因,以4號染色體的串聯(lián)重復(fù)基因最多,為29個(gè),6號染色體上串聯(lián)重復(fù)基因比例最大,為76.92%。潘那利番茄各染色體的基因數(shù)為3(0號虛擬染色體)~31(11號染色體)個(gè),平均為15.54個(gè),其中,106個(gè)(52.48%)基因形成39個(gè)基因簇,平均每簇2.72個(gè)基因,11號染色體11簇24個(gè)基因,4號染色體7簇23個(gè)基因,6號染色體最大簇8個(gè)基因,11號、4號和6號染色體成簇基因均超過70.00%,其他染色體成簇基因數(shù)為0~14個(gè);同時(shí),46個(gè)(22.77%)基因?yàn)榇?lián)重復(fù)基因,以4號染色體最多,為17個(gè),6號和4號染色體上串聯(lián)重復(fù)基因均超過50.00%。整體來看,番茄4號、6號和11號染色體NBS-LRR類抗病基因成簇和串聯(lián)重復(fù)特點(diǎn)最明顯。

    2. 4 番茄NBS-LRR類抗病蛋白系統(tǒng)進(jìn)化分析結(jié)果

    提取NBS-LRR類抗病蛋白NBS結(jié)構(gòu)域的氨基酸序列,去掉氨基酸序列過短的180個(gè)蛋白,對477個(gè)NBS-LRR類抗病蛋白進(jìn)行多重序列比對分析,結(jié)果發(fā)現(xiàn)這些抗病蛋白含有多個(gè)相對保守的氨基酸位點(diǎn),如24-I、26/29/31/127/188-G、32-K、33/34/133/134-T、36/192-A、47/166-F、102/105/152/160/189/191/200/234-L、106/107-D、128/230-S、135-R、190-P和208-W等。由系統(tǒng)發(fā)育進(jìn)化樹(圖3)可知,栽培番茄和2個(gè)野生番茄共477個(gè)NBS-LRR類抗病蛋白可聚為十大類群,其中,N端為TIR結(jié)構(gòu)域的蛋白中僅有1個(gè)蛋白聚在Ⅰ類群,其他65個(gè)蛋白聚在Ⅴ類群;N端為CC結(jié)構(gòu)域的蛋白聚在Ⅱ、Ⅲ、Ⅳ、Ⅶ、Ⅷ和Ⅹ等6個(gè)類群;N端為RPW8結(jié)構(gòu)域的蛋白聚在Ⅱ、Ⅲ、Ⅵ和Ⅸ等4個(gè)類群。因此,N端為TIR的蛋白基本聚在一起,但N端為CC和RPW8結(jié)構(gòu)域的蛋白未能很好分開,聚在多個(gè)類群中;NL和N型基因散布于各類群中。

    結(jié)合477個(gè)番茄NBS-LRR抗病基因的聚類分析結(jié)果及其核苷酸序列相似性,結(jié)果共發(fā)現(xiàn)115對直系同源基因和8對旁系同源基因,占比51.57%,主要分布在4號和5號染色體上。其中,栽培番茄分別與野生醋栗番茄和潘那利番茄形成97(84.35%)和11(9.57%)對直源基因,2個(gè)野生番茄間形成7(6.09%)對直源基因,且同源基因?qū)λ谌旧w對應(yīng),4號和5號染色體分別為27和16對,其他染色體為5~9對。野生潘那利番茄和栽培番茄各自分別形成7和1對旁源基因。計(jì)算同源基因間非同義(Ka)和同義(Ks)替換率及Ka/Ks值,結(jié)果發(fā)現(xiàn)有4對直源基因(Solyc04 g009120.2和Sopim04g009120、Solyc04g009130.3和Sopim04g009130、Solyc11g069020.2和Sopim11g06 9020、Solyc12g044190.2和Sopim12g044190)及2對旁源基因(Sopen10g029640和Sopen10g029660、Sopen11g028600和Sopen11g028610)的Ka/Ks值大于1,說明這6對同源基因受到自然正選擇作用,其他同源基因?qū)a/Ks值均小于1,為純化選擇。[Ⅰ][Ⅹ][Ⅸ][Ⅶ][Ⅵ][Ⅴ][Ⅳ][Ⅲ][Ⅱ][Ⅷ]

    2. 5 番茄NBS-LRR類抗病基因在細(xì)菌脅迫下的表達(dá)情況

    從54個(gè)栽培番茄TNL、CNL和RNL型NBS-LRR類抗病基因(其他5個(gè)基因數(shù)據(jù)缺失,2個(gè)野生番茄對應(yīng)注釋不可用,故未作分析)的表達(dá)聚類熱圖(圖4)可看出,不同基因可明顯分為兩大類群,其中,Ⅰ類群的16個(gè)基因表達(dá)量較低,Ⅱ類群的38個(gè)基因表達(dá)量較高。其中,31個(gè)基因在4或5種生物脅迫下有較高表達(dá)量,10個(gè)基因在所有生物脅迫下表達(dá)量均較低。Solyc02g082050.3和Solyc01g087200.3在2種假單胞菌和根癌農(nóng)桿菌脅迫下有較高表達(dá)量,Solyc08g0055 10.2在根癌農(nóng)桿菌脅迫下表達(dá)量最高,Agrobacterium tumefaciens和DC3000hrcQ-UfliC脅迫下基因表達(dá)量整體偏低。[Ⅰ][Ⅱ]

    3 討論

    植物NBS-LRR類抗病基因的數(shù)量和結(jié)構(gòu),在某種程度上代表了其識別病原物的能力。部分基因NBS結(jié)構(gòu)域不完整,可能與基因進(jìn)化過程中的丟失、假基因化或功能異化有關(guān)。本研究從栽培番茄及野生潘那利番茄和醋栗番茄中分別鑒定獲得238、202和217個(gè)NBS-LRR類抗病基因,共計(jì)657個(gè),根據(jù)這些基因結(jié)構(gòu)的差異進(jìn)一步將其劃分為8種類型(TNL、CNL、RNL、TN、CN、RN、NL和N),各類型中的分布數(shù)目各不相同,其中,TNL、CNL和RNL型基因數(shù)分別為58、87和13個(gè)。前人在與番茄同屬的馬鈴薯中發(fā)現(xiàn)419個(gè)NBS-LRR類抗病基因,其中TNL型基因37個(gè)、CNL型基因65個(gè),二者比值為0.57(Lozano et al.,2012),與本研究中TNL和CNL型基因數(shù)比值(0.67)基本一致,但明顯高于擬南芥的TNL和CNL型基因數(shù)比值1.80(Meyers et al.,2003)。這種不均勻分布可能是基因進(jìn)化和植物環(huán)境適應(yīng)的自然選擇結(jié)果。RNL型蛋白包含RPW8結(jié)構(gòu)域,與植物廣譜抗性密切相關(guān)(Xiao et al.,2003)。磷酸化作為重要的蛋白表觀修飾,參與基因的轉(zhuǎn)錄調(diào)控(唐文武等,2017)。而本研究發(fā)現(xiàn)番茄NBS-LRR類抗病蛋白以Ser磷酸化為主,主要是通過改變蛋白結(jié)構(gòu)激活酶活性。此外,本研究對番茄NBS-LRR類抗病基因的NBS結(jié)構(gòu)域進(jìn)行保守性分析,發(fā)現(xiàn)P-loop、Kinase2、RNBS-B和GLPL序列具有較高的相似性,與擬南芥(Meyers et al.,2003)和蝴蝶蘭(蔣卉等,2018)相關(guān)研究結(jié)果一致,表明不同物種NBS-LRR類抗病基因在進(jìn)化過程中具有較高的保守性。進(jìn)一步研究發(fā)現(xiàn),NBS保守結(jié)構(gòu)域的Kinase2末端的氨基酸(Asp或Trp)是區(qū)分CNL和TNL型的關(guān)鍵因素(Lozano et al.,2015)。本研究的染色體定位結(jié)果顯示,栽培番茄和野生潘那利番茄的NBS-LRR類抗病基因不均勻地散布在12條染色體和0號虛擬染色體上,其中以4號、5號和11號染色體較多,且主要位于染色體兩端,因此這些基因在進(jìn)化過程中可能更易變異;栽培番茄和野生潘那利番茄均有超過50.00%的基因成簇分布,且分別有35.59%和22.77%的基因?yàn)榇?lián)重復(fù)基因。可見,番茄NBS-LRR類抗病基因具有在染色體上成簇分布的特征,以4號、6號和11號染色體最突出,與擬南芥(Meyers et al.,2003)和馬鈴薯(Lozano et al.,2012)研究結(jié)果相似。基因簇形式有利于抗病基因間的序列交換,且串聯(lián)重復(fù)通常出現(xiàn)在基因簇中,表明串聯(lián)重復(fù)在NBS-LRR抗病基因家族的擴(kuò)張中具有重要作用,也是基因簇形成的主要原因(Liu et al.,2007;劉云飛等,2014a)。

    本研究通過構(gòu)建番茄NBS-LRR類抗病蛋白的系統(tǒng)發(fā)育進(jìn)化樹發(fā)現(xiàn),大部分N端為TIR結(jié)構(gòu)域的NBS-LRR類抗病蛋白聚類在一起,但N端為CC和RPW8結(jié)構(gòu)域的蛋白未能很好分開,其中以N端為CC結(jié)構(gòu)域的蛋白擴(kuò)張程度最大,散布于6個(gè)類群,其原因可能是含TIR結(jié)構(gòu)域的編碼基因在進(jìn)化過程中發(fā)生丟失現(xiàn)象,或含CC結(jié)構(gòu)域的編碼基因起源更早(Tarr and Alexander,2009)。一般而言,同一進(jìn)化分支上的基因具有較高相似性,系統(tǒng)發(fā)育進(jìn)化樹包含不同染色體上基因組成的混合進(jìn)化分支,反映出這些基因在染色體上的擴(kuò)張情況,可能是基因產(chǎn)生新功能的源泉;同一進(jìn)化分支包含不同物種基因,說明NBS-LRR類抗病基因可能按照物種特異性的方式進(jìn)行擴(kuò)張。本研究對番茄NBS-LRR類抗病基因與已克隆的NBS-LRR類抗病基因(Liu et al.,2007)進(jìn)行比對分析,結(jié)果顯示,NBS-LRR類抗病基因Solyc04g008130.2與Hero、Solyc06g008450.3與Mi-1、Solyc09g098130.2與Sw-5、Solyc09g018220.2與Tm-2、Solyc05g007850.2與Bs4、Solyc11g071995.1與I-2的相似性均在80%以上,說明其具有類似抗病蟲害功能。115對直系同源基因主要源自栽培番茄與野生醋栗番茄之間,且主要分布在4號和5號染色體上,證明其具有較高的遺傳相似性;8對旁系同源基因中有5對是串聯(lián)重復(fù)形成,其中1對來自栽培番茄、4對來自野生潘那利番茄。串聯(lián)重復(fù)基因中同源性偏低,可能與篩選條件過嚴(yán)有關(guān),也可能是基因重復(fù)后在長期演化中發(fā)生了功能異化。除了6對同源基因外,其他同源基因(95.12%)的Ka/Ks均小于1,表明其受純化選擇作用,基因進(jìn)化動力可能是遺傳漂變(Edwards,2000)。

    基于RNA-Seq的番茄NBS-LRR類抗病基因表達(dá)分析結(jié)果顯示,大多數(shù)基因在4~5種生物脅迫下表達(dá)量較高。雖然抗病基因表達(dá)促使植株產(chǎn)生抗性,但大量抗病基因的高表達(dá)會導(dǎo)致植物細(xì)胞死亡,因此為減少植物負(fù)擔(dān),抗病基因表達(dá)量通常較其他基因低(Zhai et al.,2011)。同時(shí),抗病基因的轉(zhuǎn)錄可能通過小RNA調(diào)控(Fei et al.,2013),或進(jìn)化為假基因甚至在基因組中丟失??梢姡参镌谂c環(huán)境長期互作中形成了完善的自我調(diào)節(jié)機(jī)制,以保證個(gè)體生存與群體繁衍。對于番茄NBS-LRR類抗病基因家族仍需進(jìn)行針對性地生物學(xué)功能鑒定及挖掘,以推動相關(guān)抗病基因的農(nóng)業(yè)生產(chǎn)應(yīng)用。

    4 結(jié)論

    番茄NBS-LRR類抗病基因具有成簇存在的特點(diǎn),在進(jìn)化過程中經(jīng)重復(fù)擴(kuò)張,已形成大量同源基因,且能響應(yīng)多種生物脅迫。

    參考文獻(xiàn):

    蔣卉,張晶,符真珠,董曉宇,王慧娟,李艷敏,高杰,王利民,張和臣. 2018. 蝴蝶蘭NBS-LRR家族基因挖掘和生物信息學(xué)分析[J]. 分子植物育種,16(9):2786-2794. doi:10. 13271/j.mpb.016.002786. [Jiang H,Zhang J,F(xiàn)u Z Z,Dong X Y,Wang H J,Li Y M,Gao J,Wang L M,Zhang H C. 2018. Mining and bioinformatics analysis of NBS-LRR gene family in Phalaenopsis[J]. Molecular Plant Breeding,16(9):2786-2794.]

    李任建,申哲源,李旭凱,韓淵懷,張寶俊. 2020. 谷子NBS-LRR類基因家族全基因組鑒定及表達(dá)分析[J]. 河南農(nóng)業(yè)科學(xué),49(2):34-43. doi:10.15933/j.cnki.1004-3268.2020. 02.005. [Li R J,Shen Z Y,Li X K,Han Y H,Zhang B J. 2020. Genome-wide identification and expression analysis of NBS-LRR gene family in Setaria italica[J]. Journal of Henan Agricultural Sciences,49(2):34-43.]

    劉潮,褚洪龍,韓利紅,楊云錦,高永,唐利洲. 2019. 桑樹NBS-LRR類基因家族的全基因組鑒定及其調(diào)控micro-RNAs分析[J]. 江蘇農(nóng)業(yè)學(xué)報(bào),35(3):544-553. doi:10.3969/ j.issn.1000-4440.2019.03.007. [Liu C,Chu H L,Han L H,Yang Y J,Gao Y,Tang L Z. 2019. Genome-wide identification of NBS-LRR genes and regulation analysis by microRNAs in mulberry[J]. Jiangsu Journal of Agricultural Sciences,35(3):544-553.]

    劉云飛,萬紅建,李志邈,葉青靜,王榮青,阮美穎,姚祝平,周國治,韋艷萍,楊悅儉. 2014a. 植物NBS-LRR抗病基因的結(jié)構(gòu)、功能、進(jìn)化起源及其應(yīng)用[J]. 分子植物育種,12(2):377-389. doi:10.13271/j.mpb.012.000377. [Liu Y F,Wan H J,Li Z M,Ye Q J,Wang R Q,Ruan M Y,Yao Z P,Zhou G Z,Wei Y P,Yang Y J. 2014. Analysis of plant NBS-LRR resistance gene:structure,function,origin,evolution and their application[J]. Molecular Plant Breeding,12(2):377-389.]

    劉云飛,萬紅建,韋艷萍,李志邈,葉青靜,王榮青,阮美穎,姚祝平,周國治,楊悅儉. 2014b. 番茄NBS-LRR抗病基因家族全基因組分析[J]. 核農(nóng)學(xué)報(bào),28(5):790-799. [Liu Y F,Wan H J,Wei Y P,Li Z M,Ye Q J,Wang R Q,Ruan M Y,Yao Z P,Zhou G Z,Yang Y J. 2014. Genome-wide analysis of NBS-LRR resistance genes in tomato[J]. Journal of Nuclear Agricultural Sciences,28(5):790-799.]

    唐文武,吳秀蘭,李桂花. 2017. 蕓薹屬作物花粉發(fā)育相關(guān)基因MS1的生物信息學(xué)分析[J]. 南方農(nóng)業(yè)學(xué)報(bào),48(12):2122-2128. doi:10.3969/j.issn.2095-1191.2017.12.02. [Tang W W,Wu X L,Li G H. 2017. Bioinformatics analysis for pollen development related gene MS1 of Brassica crops[J]. Journal of Southern Agriculture,48(12):2122-2128.]

    文歡,劉永立. 2015. ‘紅陽獼猴桃全基因組NBS-LRR類基因家族的生物信息學(xué)分析[J]. 果樹學(xué)報(bào),32(2):413-419. doi:10.13925/j.cnki.gsxb.20140408. [Wen H,Liu Y L. 2015. Bioinformatic analysis of the NBS-LRR gene family in kiwifruit ‘Hongyang[J]. Journal of Fruit Science,32(2):413-419.]

    Ameline-Torregrosa C,Wang B B,O'Bleness M S,Deshpande S,Zhu H,Roe B,Young N D,Cannon S B. 2008. Identification and characterization of nucleotide-binding site-leucine-rich repeat genes in the model plant Medicago truncatula[J]. Plant Physiology,146(1):5-21. doi:10. 1104/pp.107.104588.

    Bai J F,Pennill L A,Ning J C,Lee S W,Ramalingam J,Webb C A,Zhao B Y,Sun Q,Nelson J C,Leach J E,Hulbert S H. 2002. Diversity in nucleotide binding site-leucine-rich repeat genes in cereals[J]. Genome Research,12(12):1871-1884. doi:10.1101/gr.454902.

    Cheng Y,Li X,Jiang H,Ma W,Miao W,Yamada T,Zhang M. 2012. Systematic analysis and comparison of nucleotide-binding site disease resistance genes in maize[J]. The FEBS Journal,279(13):2431-2443. doi:10.1111/j.1742-4658. 2012.08621.x.

    Dangl J L,Jones J D. 2001. Plant pathogens and integrated defence responses to infection[J]. Nature,411(6839):826-833. doi:10.1038/35081161.

    Edwards A W F. 2000. The genetical theory of natural selection[J]. Genetics,154:1419-1426. doi:10.1038/126595a0.

    Fei Q L,Xia R,Meyers B C. 2013. Phased,secondary,small interfering RNAs in posttranscriptional regulatory networks[J]. The Plant Cell,25(7):2400-2415. doi:10.1105/ tpc.113.114652.

    Fenyk S,Campillo A D S E,Pohl E,Hussey P J,Cann M J. 2012. A nucleotide phosphatase activity in the nucleotide binding domain of an orphan resistance protein from rice[J]. Journal of Biological Chemistry,287(6):402-4032. doi:10.1074/jbc.M111.314450.

    Habachi-Houimli Y,Khalfallah Y,Mezghani-Khemakhem M,Makni H,Makni M,Bouktila D. 2018. Genome-wide identification,characterization,and evolutionary analysis of NBS-encoding resistance genes in barley[J]. 3 Biotech,8(11):453. doi:10.1007/s13205-018-1478-6.

    Huang S,Gao Y,Liu J,Peng X,Niu X,F(xiàn)ei Z,Cao S,Liu Y. 2012. Genome-wide analysis of WRKY transcription factors in Solanum lycopersicum[J]. Molecular Genetics and Genomics,287(6):495-513. doi:10.1007/s00438-012-0696-6.

    Liu J L,Liu X L,Dai L Y,Wang G L. 2007. Recent progress in elucidating the structure,function and evolution of di-sease resistance genes in plants[J]. Journal of Genetics and Genomics,34(9):765-776. doi:10.1016/S1673-8527(07)60087-3.

    Lozano R,Ponce O,Ramirez M,Mostajo N,Orjeda G. 2012. Genomewide identification and mapping of NBS-enco-ding resistance genes in Solanum tuberosum group phureja[J]. PLoS One,7(4):e34775. doi:10.1186/s12864-015-1554-9.

    Lozano R,Hamblin M T,Prochnik S,Jannink J L. 2015. Identification and distribution of the NBS-LRR gene family in the cassava genome[J]. BMC Genomics,16:360. doi:10.1371/journal.pone.0034775.

    Meyers B C,Kozik A,Griego A,Kuang H,Michelmore R W. 2003. Genome-wide analysis of NBS-LRR-encoding genes in Arabidopsis[J]. The Plant Cell,15(4):809-834. doi:10.1105/tpc.009308.

    Sato S,Tabata S,Hirakawa H,Asamizu E,Shirasawa K,Isobe S,Kaneko T,Nakamura Y,Shibata D,Aoki K. 2012. The tomato genome sequence provides insights into fleshy fruit evolution[J]. Nature,485(7400):635-641. doi:10. 1038/nature11119.

    Tarr D E K,Alexander H M. 2009. TIR-NBS-LRR genes are rare in monocots:Evidence from diverse monocot orders[J]. BMC Research Notes,2(197):1-10. doi:10.1186/1756-0500-2-197.

    Xiao S Y,Brown S,Patrick E. 2003. Enhanced transcription of the Arabidopsis disease resistance genes RPW8.1 and RPW8.2 via a salicylic acid-dependent amplification circuit is required for hypersensitive cell death[J]. The Plant Cell,15(1):33-45. doi:10.1105/tpc.006940.

    Xiao S,Ellwood S,Calis O,Patrick E,Li T,Coleman M,Turner J G. 2001. Broad-spectrum mildew resistance in Arabidopsis thaliana mediated by RPW8[J]. Science,291(5501):118-120. doi:10.1126/science.291.5501.118.

    Yang S,Zhang X,Yue J X,Tian D,Chen J Q. 2008. Recent duplications dominate NBS-encoding gene expansion in two woody species[J]. Molecular Genetics and Genomics,280(3):187-198. doi:10.1007/s00438-008-0355-0.

    Zhai J X,Jeong D H,De Paoli E,Park S,Rosen B D,Li Y P,González A J,Yan Z,Kitto S L,Grusak M A,Jackson S A,Stacey G,Cook D R,Green P J,Sherrier D J,Meyers B C. 2011. Micrornas as master regulators of the plant NBS-LRR defense gene family via the production of phased,transacting siRNAs[J]. Genes and Development,25(23):2540-2553. doi:10.1101/gad.177527.111.

    (責(zé)任編輯 陳 燕)

    猜你喜歡
    結(jié)構(gòu)域抗病氨基酸
    科學(xué)解讀食物中的蛋白質(zhì)如何分“優(yōu)劣”
    植物細(xì)胞內(nèi)存在“自殺神器”
    補(bǔ)氨基酸不如吃雞蛋
    結(jié)核分枝桿菌蛋白Rv0089的生物信息學(xué)分析
    黃星天牛中腸中內(nèi)切葡聚糖酶的鑒定與酶活性測定
    蛋白質(zhì)功能預(yù)測方法研究進(jìn)展
    優(yōu)質(zhì)高產(chǎn)小麥新品種寧麥24的選育及栽培技術(shù)
    復(fù)方氨基酸注射液不合理應(yīng)用分析
    RGA法標(biāo)記植物抗病基因的研究進(jìn)展
    2007年上海生物高考卷部分遺傳題解密
    亚洲av欧美aⅴ国产| 国产1区2区3区精品| 欧美精品啪啪一区二区三区 | 丝瓜视频免费看黄片| 午夜免费观看性视频| 国产精品av久久久久免费| 日本一区二区免费在线视频| 国产激情久久老熟女| 精品少妇久久久久久888优播| a级毛片黄视频| 久久久久久人人人人人| 十八禁网站网址无遮挡| 最黄视频免费看| 两个人看的免费小视频| 在线观看www视频免费| 成年动漫av网址| 狠狠婷婷综合久久久久久88av| 人人妻人人澡人人看| 黑人欧美特级aaaaaa片| 精品少妇一区二区三区视频日本电影| 热99久久久久精品小说推荐| 老鸭窝网址在线观看| 一本色道久久久久久精品综合| av国产精品久久久久影院| 一区二区三区激情视频| 高清视频免费观看一区二区| 免费人妻精品一区二区三区视频| 成人影院久久| 欧美日韩黄片免| 国产精品国产av在线观看| 99国产精品一区二区蜜桃av | 国产午夜精品一二区理论片| 视频在线观看一区二区三区| svipshipincom国产片| 欧美国产精品va在线观看不卡| 新久久久久国产一级毛片| 90打野战视频偷拍视频| 精品国产一区二区久久| 少妇猛男粗大的猛烈进出视频| 欧美亚洲日本最大视频资源| 十八禁网站网址无遮挡| 国产麻豆69| 国产成人精品久久二区二区免费| 自线自在国产av| 欧美激情极品国产一区二区三区| 麻豆国产av国片精品| 国产一区二区 视频在线| 国产精品免费视频内射| 男女床上黄色一级片免费看| 亚洲色图 男人天堂 中文字幕| 我的亚洲天堂| 成年av动漫网址| 19禁男女啪啪无遮挡网站| 69精品国产乱码久久久| 亚洲成人手机| 亚洲av在线观看美女高潮| 日韩熟女老妇一区二区性免费视频| 欧美人与性动交α欧美精品济南到| 校园人妻丝袜中文字幕| 亚洲专区中文字幕在线| 欧美日韩亚洲综合一区二区三区_| 国产精品秋霞免费鲁丝片| 中文字幕av电影在线播放| 国产高清视频在线播放一区 | 十八禁人妻一区二区| svipshipincom国产片| 丝袜脚勾引网站| 搡老岳熟女国产| 欧美激情高清一区二区三区| 免费在线观看黄色视频的| 91老司机精品| 国产一级毛片在线| 国产伦理片在线播放av一区| 日本av免费视频播放| 久久人人爽av亚洲精品天堂| 黄色 视频免费看| 深夜精品福利| 精品国产一区二区久久| 亚洲精品国产av成人精品| 亚洲欧美色中文字幕在线| 免费观看人在逋| 一本色道久久久久久精品综合| 夜夜骑夜夜射夜夜干| 亚洲成人免费av在线播放| 天天添夜夜摸| 国产一区亚洲一区在线观看| 中文乱码字字幕精品一区二区三区| 天天影视国产精品| 永久免费av网站大全| 可以免费在线观看a视频的电影网站| 老司机深夜福利视频在线观看 | 性色av乱码一区二区三区2| 国产精品亚洲av一区麻豆| 91老司机精品| 热re99久久精品国产66热6| 亚洲,欧美,日韩| 欧美人与性动交α欧美精品济南到| 日日摸夜夜添夜夜爱| 丝瓜视频免费看黄片| 99精国产麻豆久久婷婷| 国产成人欧美在线观看 | 久久精品亚洲av国产电影网| 国产av国产精品国产| 天天躁夜夜躁狠狠躁躁| 亚洲欧美一区二区三区黑人| 国产精品三级大全| 黄色 视频免费看| 看十八女毛片水多多多| 色综合欧美亚洲国产小说| 国产黄色视频一区二区在线观看| 国产伦理片在线播放av一区| 又大又黄又爽视频免费| 国产老妇伦熟女老妇高清| 久久亚洲精品不卡| 亚洲人成电影免费在线| 精品少妇内射三级| 国产淫语在线视频| 别揉我奶头~嗯~啊~动态视频 | 国产在线一区二区三区精| 99国产综合亚洲精品| 久久久久视频综合| 国产精品一二三区在线看| 九草在线视频观看| xxx大片免费视频| 97精品久久久久久久久久精品| 国产伦理片在线播放av一区| 亚洲欧美成人综合另类久久久| av电影中文网址| 青春草视频在线免费观看| www.999成人在线观看| 久久99精品国语久久久| 9色porny在线观看| 精品少妇一区二区三区视频日本电影| 午夜福利在线免费观看网站| www.精华液| 日韩一区二区三区影片| 国产欧美日韩精品亚洲av| 黄色片一级片一级黄色片| 日韩免费高清中文字幕av| 国产淫语在线视频| 国产极品粉嫩免费观看在线| 不卡av一区二区三区| 午夜福利乱码中文字幕| 国产极品粉嫩免费观看在线| 成年人午夜在线观看视频| 久久久欧美国产精品| 成年人黄色毛片网站| 国产免费视频播放在线视频| 少妇精品久久久久久久| 母亲3免费完整高清在线观看| 制服人妻中文乱码| 极品少妇高潮喷水抽搐| 免费人妻精品一区二区三区视频| 好男人视频免费观看在线| 搡老岳熟女国产| 婷婷色麻豆天堂久久| 亚洲精品久久午夜乱码| 日韩视频在线欧美| 亚洲精品久久久久久婷婷小说| 久久99热这里只频精品6学生| 亚洲图色成人| 男女边摸边吃奶| 50天的宝宝边吃奶边哭怎么回事| 精品人妻1区二区| 欧美精品人与动牲交sv欧美| 久久99热这里只频精品6学生| 91成人精品电影| 久久天躁狠狠躁夜夜2o2o | 你懂的网址亚洲精品在线观看| 国产亚洲午夜精品一区二区久久| 国产一级毛片在线| 伊人亚洲综合成人网| 人人澡人人妻人| 亚洲精品一区蜜桃| 欧美日韩成人在线一区二区| 国语对白做爰xxxⅹ性视频网站| 一区福利在线观看| 日韩制服骚丝袜av| 狂野欧美激情性bbbbbb| 国产av精品麻豆| 国产成人啪精品午夜网站| 免费看av在线观看网站| 国产福利在线免费观看视频| 看免费av毛片| 欧美变态另类bdsm刘玥| 性少妇av在线| 性色av一级| 亚洲人成77777在线视频| 国产激情久久老熟女| 欧美 日韩 精品 国产| 免费观看av网站的网址| 两个人看的免费小视频| avwww免费| 国产熟女午夜一区二区三区| av不卡在线播放| 18禁国产床啪视频网站| 国产成人欧美在线观看 | 久久久精品国产亚洲av高清涩受| 欧美日韩国产mv在线观看视频| 午夜福利乱码中文字幕| 国产国语露脸激情在线看| 久久久久网色| 90打野战视频偷拍视频| www.熟女人妻精品国产| 午夜福利视频在线观看免费| 成年美女黄网站色视频大全免费| 可以免费在线观看a视频的电影网站| 亚洲精品第二区| 久久女婷五月综合色啪小说| 国产高清国产精品国产三级| 午夜两性在线视频| a级片在线免费高清观看视频| 国产国语露脸激情在线看| 日本av手机在线免费观看| 晚上一个人看的免费电影| 波野结衣二区三区在线| 91老司机精品| 国产免费又黄又爽又色| 少妇人妻 视频| 日韩 亚洲 欧美在线| 国产亚洲欧美在线一区二区| 国产亚洲午夜精品一区二区久久| av天堂在线播放| 视频在线观看一区二区三区| 婷婷色综合www| 中文字幕人妻丝袜一区二区| 成人三级做爰电影| 新久久久久国产一级毛片| 嫁个100分男人电影在线观看 | 视频在线观看一区二区三区| 美国免费a级毛片| 宅男免费午夜| 大片电影免费在线观看免费| 91字幕亚洲| 亚洲国产欧美一区二区综合| av国产久精品久网站免费入址| 另类亚洲欧美激情| 精品久久蜜臀av无| 国产在线一区二区三区精| 亚洲专区中文字幕在线| 一区二区三区精品91| 国产精品久久久久成人av| 一二三四社区在线视频社区8| 国产成人欧美在线观看 | 日韩中文字幕视频在线看片| 中文字幕亚洲精品专区| 日韩大码丰满熟妇| 五月天丁香电影| 亚洲国产日韩一区二区| 国产亚洲午夜精品一区二区久久| 女人高潮潮喷娇喘18禁视频| 十八禁人妻一区二区| 一级毛片我不卡| 亚洲国产精品999| 国产精品亚洲av一区麻豆| 国产成人欧美| 亚洲激情五月婷婷啪啪| 夜夜骑夜夜射夜夜干| 国产午夜精品一二区理论片| 久久久久国产一级毛片高清牌| 每晚都被弄得嗷嗷叫到高潮| 亚洲av国产av综合av卡| 亚洲av片天天在线观看| 亚洲精品美女久久av网站| www.精华液| 新久久久久国产一级毛片| 后天国语完整版免费观看| 在线 av 中文字幕| 欧美日韩综合久久久久久| 久久热在线av| 人妻 亚洲 视频| 黄色片一级片一级黄色片| 亚洲成人免费电影在线观看 | 午夜激情av网站| 国产av精品麻豆| 免费在线观看日本一区| 深夜精品福利| 久久国产精品人妻蜜桃| 久久精品aⅴ一区二区三区四区| 极品人妻少妇av视频| 一区二区三区激情视频| 国产黄色免费在线视频| 国产精品久久久久久精品电影小说| 老司机影院毛片| 黄色a级毛片大全视频| 一区二区日韩欧美中文字幕| 国产一区亚洲一区在线观看| videos熟女内射| 69精品国产乱码久久久| 婷婷色av中文字幕| 亚洲精品久久久久久婷婷小说| 老司机影院成人| 90打野战视频偷拍视频| 校园人妻丝袜中文字幕| 久久久欧美国产精品| 可以免费在线观看a视频的电影网站| 欧美 日韩 精品 国产| 日韩一本色道免费dvd| 不卡av一区二区三区| 波野结衣二区三区在线| 精品少妇内射三级| 久久久国产精品麻豆| 国产一区二区 视频在线| 一区福利在线观看| 美女扒开内裤让男人捅视频| 日韩视频在线欧美| 国产男人的电影天堂91| 国产真人三级小视频在线观看| 少妇人妻久久综合中文| 精品人妻1区二区| 99香蕉大伊视频| 国产日韩欧美视频二区| 亚洲精品久久久久久婷婷小说| 美女国产高潮福利片在线看| 在线看a的网站| 亚洲国产最新在线播放| 99国产综合亚洲精品| 99国产精品99久久久久| 国产一区二区激情短视频 | 国产精品久久久久久精品古装| 欧美精品人与动牲交sv欧美| 亚洲欧美中文字幕日韩二区| 永久免费av网站大全| 中文字幕精品免费在线观看视频| 99精国产麻豆久久婷婷| 一区二区三区乱码不卡18| 亚洲av日韩精品久久久久久密 | 亚洲视频免费观看视频| 一级毛片电影观看| 久久久欧美国产精品| 国产精品三级大全| 精品人妻在线不人妻| 亚洲综合色网址| 新久久久久国产一级毛片| 亚洲av成人精品一二三区| 欧美 日韩 精品 国产| 亚洲美女黄色视频免费看| 99热全是精品| 久久精品亚洲av国产电影网| 在线观看www视频免费| 18在线观看网站| 亚洲欧美中文字幕日韩二区| 国产精品久久久久久人妻精品电影 | 久久久欧美国产精品| www.精华液| 国产成人系列免费观看| 又粗又硬又长又爽又黄的视频| 欧美av亚洲av综合av国产av| 大码成人一级视频| 国产xxxxx性猛交| 亚洲av男天堂| 涩涩av久久男人的天堂| 极品人妻少妇av视频| 亚洲国产成人一精品久久久| 99香蕉大伊视频| 另类精品久久| 最新的欧美精品一区二区| 国产1区2区3区精品| 国产不卡av网站在线观看| 精品一区在线观看国产| a级片在线免费高清观看视频| 国产人伦9x9x在线观看| 亚洲男人天堂网一区| 国产一区二区 视频在线| 18禁国产床啪视频网站| 50天的宝宝边吃奶边哭怎么回事| 成年人免费黄色播放视频| 黄色a级毛片大全视频| 男女床上黄色一级片免费看| 日日爽夜夜爽网站| www.av在线官网国产| 看免费av毛片| 人人妻人人澡人人爽人人夜夜| 真人做人爱边吃奶动态| 久久女婷五月综合色啪小说| 亚洲精品国产一区二区精华液| 国产精品一二三区在线看| 国产伦人伦偷精品视频| 国产精品久久久久成人av| 亚洲欧美一区二区三区黑人| 国产片特级美女逼逼视频| 国产成人精品久久久久久| 亚洲精品日本国产第一区| 久久精品国产亚洲av高清一级| 2021少妇久久久久久久久久久| 中文字幕另类日韩欧美亚洲嫩草| a级毛片黄视频| 一级片免费观看大全| 黄色视频在线播放观看不卡| av网站免费在线观看视频| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品国产av蜜桃| 午夜老司机福利片| 日本91视频免费播放| av又黄又爽大尺度在线免费看| 天天影视国产精品| 丰满人妻熟妇乱又伦精品不卡| 欧美变态另类bdsm刘玥| 久久久久久久久免费视频了| 五月开心婷婷网| 91精品伊人久久大香线蕉| 天天添夜夜摸| a级毛片黄视频| 久久天躁狠狠躁夜夜2o2o | 少妇人妻久久综合中文| 国产又爽黄色视频| 精品一区在线观看国产| 日本黄色日本黄色录像| 热re99久久国产66热| 亚洲九九香蕉| 亚洲欧美中文字幕日韩二区| 免费人妻精品一区二区三区视频| 在线亚洲精品国产二区图片欧美| 桃花免费在线播放| 国产真人三级小视频在线观看| 亚洲欧洲日产国产| 男女床上黄色一级片免费看| 99久久综合免费| 啦啦啦在线免费观看视频4| 捣出白浆h1v1| 午夜福利视频精品| 久久人妻熟女aⅴ| 下体分泌物呈黄色| 国产成人影院久久av| 日韩中文字幕欧美一区二区 | 国产精品久久久av美女十八| www.熟女人妻精品国产| 久久久久精品人妻al黑| 咕卡用的链子| 大型av网站在线播放| 久热爱精品视频在线9| 国产黄频视频在线观看| 亚洲国产成人一精品久久久| 成人国语在线视频| 免费少妇av软件| 美女扒开内裤让男人捅视频| 欧美黑人欧美精品刺激| 王馨瑶露胸无遮挡在线观看| 国产免费福利视频在线观看| 国产视频一区二区在线看| 两性夫妻黄色片| 亚洲 国产 在线| 爱豆传媒免费全集在线观看| 考比视频在线观看| av国产精品久久久久影院| 成人亚洲欧美一区二区av| 一区二区三区乱码不卡18| 国产女主播在线喷水免费视频网站| 亚洲欧美精品自产自拍| 亚洲精品国产av成人精品| 老鸭窝网址在线观看| 久久久久久免费高清国产稀缺| 大香蕉久久成人网| 成年av动漫网址| 可以免费在线观看a视频的电影网站| 人妻 亚洲 视频| 在线观看国产h片| 久久久久精品人妻al黑| 久久精品国产综合久久久| 精品免费久久久久久久清纯 | 国产欧美日韩综合在线一区二区| 新久久久久国产一级毛片| 又粗又硬又长又爽又黄的视频| 国产91精品成人一区二区三区 | 精品一区在线观看国产| 国产在线视频一区二区| 高清欧美精品videossex| 欧美精品人与动牲交sv欧美| 老司机靠b影院| 国产成人91sexporn| 国产xxxxx性猛交| 久久ye,这里只有精品| 日本黄色日本黄色录像| 欧美 亚洲 国产 日韩一| 亚洲黑人精品在线| 悠悠久久av| 操美女的视频在线观看| 国产av一区二区精品久久| 国产一区二区 视频在线| 午夜福利免费观看在线| www日本在线高清视频| 亚洲激情五月婷婷啪啪| 天天躁夜夜躁狠狠躁躁| 亚洲国产欧美在线一区| 亚洲,欧美,日韩| 亚洲av片天天在线观看| 久久免费观看电影| 黄网站色视频无遮挡免费观看| 一区二区av电影网| 日本欧美视频一区| 亚洲欧美日韩高清在线视频 | 亚洲国产精品999| 亚洲,一卡二卡三卡| 欧美日韩亚洲综合一区二区三区_| 下体分泌物呈黄色| 欧美日韩一级在线毛片| 精品亚洲成国产av| 国产精品 欧美亚洲| 捣出白浆h1v1| 老司机深夜福利视频在线观看 | 亚洲国产精品一区三区| 久久女婷五月综合色啪小说| 亚洲欧美精品综合一区二区三区| av福利片在线| 久久狼人影院| 91老司机精品| 成年美女黄网站色视频大全免费| 亚洲中文字幕日韩| 极品人妻少妇av视频| 精品一品国产午夜福利视频| 久久 成人 亚洲| 久久99一区二区三区| 亚洲精品在线美女| 操美女的视频在线观看| 亚洲成人国产一区在线观看 | 母亲3免费完整高清在线观看| 婷婷色av中文字幕| 一本—道久久a久久精品蜜桃钙片| 国产一卡二卡三卡精品| 大片免费播放器 马上看| 日韩一本色道免费dvd| 91字幕亚洲| 高清av免费在线| 久久鲁丝午夜福利片| 欧美 亚洲 国产 日韩一| 欧美精品一区二区免费开放| 日本91视频免费播放| 欧美久久黑人一区二区| 国产男女超爽视频在线观看| 一级黄色大片毛片| 精品一区二区三区四区五区乱码 | 婷婷丁香在线五月| 国产成人一区二区在线| 亚洲欧美清纯卡通| 一区二区三区精品91| 国产亚洲av片在线观看秒播厂| 欧美另类一区| 成人亚洲精品一区在线观看| 欧美日韩精品网址| 国产麻豆69| av视频免费观看在线观看| 黄色毛片三级朝国网站| av网站免费在线观看视频| 91字幕亚洲| 国产亚洲午夜精品一区二区久久| 国产亚洲欧美在线一区二区| 老司机在亚洲福利影院| 日本一区二区免费在线视频| h视频一区二区三区| av在线播放精品| 亚洲黑人精品在线| 狠狠精品人妻久久久久久综合| 欧美日韩亚洲国产一区二区在线观看 | 又大又爽又粗| 成年av动漫网址| 久久久久精品国产欧美久久久 | 欧美精品人与动牲交sv欧美| 999久久久国产精品视频| 欧美另类一区| 老司机影院毛片| 国产亚洲av高清不卡| av片东京热男人的天堂| 熟女av电影| 免费在线观看视频国产中文字幕亚洲 | 日本色播在线视频| 午夜影院在线不卡| 日日爽夜夜爽网站| 欧美在线一区亚洲| 亚洲av成人不卡在线观看播放网 | 欧美日韩精品网址| 午夜老司机福利片| 午夜福利一区二区在线看| 中文字幕高清在线视频| 交换朋友夫妻互换小说| 男人舔女人的私密视频| 国产精品三级大全| 老司机影院毛片| 秋霞在线观看毛片| 久久久久久亚洲精品国产蜜桃av| 熟女av电影| 色婷婷久久久亚洲欧美| 后天国语完整版免费观看| 在线观看免费视频网站a站| 丁香六月欧美| 日韩大码丰满熟妇| 亚洲少妇的诱惑av| 黄色a级毛片大全视频| 大话2 男鬼变身卡| 自线自在国产av| 中文字幕最新亚洲高清| 亚洲欧美一区二区三区黑人| 中文字幕亚洲精品专区| 国产高清videossex| 久久中文字幕一级| 老汉色∧v一级毛片| 国产成人啪精品午夜网站| videosex国产| 麻豆乱淫一区二区| 成人国产av品久久久| 亚洲伊人色综图| 亚洲精品国产色婷婷电影| 黑丝袜美女国产一区| 美女视频免费永久观看网站| 午夜久久久在线观看| 亚洲人成77777在线视频| 麻豆乱淫一区二区| 亚洲综合色网址| 大香蕉久久网| 日韩中文字幕欧美一区二区 | 欧美人与善性xxx| 国产片特级美女逼逼视频| 日本vs欧美在线观看视频| 在线精品无人区一区二区三| 亚洲欧美精品综合一区二区三区|