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

    木薯NOX基因家族成員的生物信息學(xué)及其表達(dá)分析

    2019-09-10 07:22:44梁大成
    關(guān)鍵詞:木薯外顯子結(jié)構(gòu)域

    胡 廣 梁大成

    摘要:【目的】搜索鑒定木薯還原型煙酰胺腺嘌呤二核苷酸磷酸氧化酶(NOX)基因家族成員,分析其生物學(xué)信息及在不同組織和非生物脅迫下的表達(dá)模式,為深入研究該家族基因?qū)δ臼韷K根發(fā)育的調(diào)控機(jī)理及提高木薯抗逆性提供理論參考?!痉椒ā繌腏GI數(shù)據(jù)庫(kù)下載木薯的基因組序列,由pfam數(shù)據(jù)庫(kù)下載4個(gè)NOX蛋白結(jié)構(gòu)域相關(guān)的HMM文件,通過(guò)HMMER程序篩選NOX基因家族成員,利用生物信息學(xué)分析軟件對(duì)該家族蛋白的理化性質(zhì)、氨基酸序列和結(jié)構(gòu)特征進(jìn)行分析。基于GEO數(shù)據(jù)庫(kù)的轉(zhuǎn)錄組測(cè)序(RNA-Seq)數(shù)據(jù),對(duì)木薯NOX基因家族成員在不同組織和干旱脅迫下的表達(dá)模式進(jìn)行分析?!窘Y(jié)果】從木薯基因組中共鑒定到15個(gè)NOX基因家族成員,包括9個(gè)NOX基因和6個(gè)FRO基因(NOX祖先基因)。FRO基因的外顯子數(shù)(6~9個(gè))明顯少于NOX基因的外顯子數(shù)(11~14個(gè)),且FRO基因的內(nèi)含子序列較NOX基因短,NOX基因的第一個(gè)外顯子相對(duì)保守。15個(gè)NOX基因家族成員分布在木薯的9條染色體和1個(gè)結(jié)構(gòu)支架上,且均被定位在染色體末端,除6號(hào)染色體有4個(gè)基因、8號(hào)和14號(hào)染色體有2個(gè)基因外,其他染色體均只有1個(gè)基因。木薯NOX蛋白的分子量和氨基酸數(shù)明顯高于FRO蛋白,而理論等電點(diǎn)(pI)差異不明顯。FRO蛋白含6~9個(gè)保守跨膜結(jié)構(gòu)域,而NOX蛋白均只含4個(gè)跨膜螺旋區(qū),且均位于第350~750位氨基酸殘基。木薯NOX家族蛋白有10個(gè)保守結(jié)構(gòu)域,其中Motif1和Motif8存在于所有的NOX和FRO蛋白中,Motif6為NOX蛋白所特有,Motif3只獨(dú)立存在于FRO蛋白中。NOX蛋白均具有典型的保守結(jié)構(gòu)域和相似的蛋白質(zhì)結(jié)構(gòu),其脫氫酶結(jié)構(gòu)域明顯折疊成催化域裂縫結(jié)構(gòu),且C末端氨基酸極度保守。由RNA-Seq數(shù)據(jù)的基因表達(dá)譜可知,NOX基因家族成員在不同木薯葉片、中靜脈、葉柄、莖、側(cè)芽、莖尖分生組織、根尖分生組織、易碎胚性愈傷組織、胚性愈傷組織(OES)和貯藏根中均有表達(dá),但呈不同的表達(dá)模式,其中,Manes.14G042600和Manes.S089900在儲(chǔ)藏根中高表達(dá),Manes.06G128700和Manes.09G172500在須根和根尖分生組織中高表達(dá)。在木薯品種KU50和新選048的葉片中,Manes.14G042600基因在干旱脅迫下呈上調(diào)表達(dá),推測(cè)NOX基因在特定組織及響應(yīng)非生物脅迫過(guò)程中扮演重要角色?!窘Y(jié)論】木薯NOX基因家族成員具有保守基因結(jié)構(gòu),其編碼蛋白具有保守的功能域,尤其是C端氨基酸序列;NOX基因表達(dá)具有組織特異性,部分成員的表達(dá)受非生物脅迫的影響。FRO蛋白和NOX蛋白具有較明顯的共進(jìn)化關(guān)系。

    關(guān)鍵詞: 木薯;還原型煙酰胺腺嘌呤二核苷酸磷酸氧化酶(NOX);基因家族成員;生物信息學(xué);干旱脅迫

    中圖分類號(hào): S533.035.3? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)志碼: A 文章編號(hào):2095-1191(2019)10-2178-10

    Bioinformatics of cassava NOX gene family and

    its expression analysis

    HU Guang1,2, LIANG Da-cheng1,2*

    (1Engineering Research Center of Ecology and Agricultural Use of Wetland, Ministry of Education/Hubei? Key Laboratory of Waterlogging Disaster and Wetland Agriculture, Jingzhou, Hubei? 434025, China;

    2Hubei Collaborative Innovation Center for Grain Industry,Jingzhou, Hubei? 434025, China)

    Abstract:【Objective】Gene family members of reduced nicotinamide adenine dinucleotide phosphate oxidase(NOX) in cassava were identified through bioinformatic mining of publicly available data,and gene expression profiles in di-fferent tissues and under various abiotic stresses were also analyzed,which laid? theoretical references for further studies on the regulation mechanism of NOX genes on cassava tuberization and improvement of stress resistance of cassava. 【Method】Cassava genome sequences were downloaded from GJI database and four NOX domain-related HMM files were downloaded from the pfam protein database,then the HMMER program was performed to characterize the NOX family members. The physical and chemical properties,amino acid sequence and structural characteristics of the family protein were analyzed by bioinformatics software. Based on transcriptome sequencing(RNA-Seq) data at GEO database,the expression profiles of cassava NOX gene family members was further analyzed in various tissues and under drought stress. 【Result】A total of 15 members of the NOX gene family were identified from the cassava genome,including 9 NOX genes,6 FRO genes (or the ancestral NOX genes). The exons of NOX(11-14) gene greatly outnumbered that of the FRO family(6-9). The intron sequence of FRO gene was shorter than that of NOX gene and the first exon of NOX gene was relatively conserved. Fifteen members of the NOX gene family were distributed on 9 different chromosomes and 1 scaffold of cassava,all of them were anchored near the endpoints of chromosomes. Except for chromosome 6(4 genes) and chromosome 14(2 genes),there was only one NOX gene in each chromosome. The molecular weight and amino acid number of NOX protein in cassava were? lower than FRO protein, but the theoretical isoelectric point(pI) difference was not obvious. FRO proteins harbored 6 to 9 conserved transmembrane domains,whereas all the NOX proteins had only 4 transmembrane domains. There are 10 conserved domains of NOX family proteins in cassava, with Motif1 and Motif8 in all NOX and FRO proteins, Motif6 only in NOX proteins, and Motif3 only in FRO proteins. Three-dimension homologous modeling further revealed that all NOX proteins had typical conserved domains and similar protein structures. The dehydrogenase domain folded obviously into the catalytic domain fissure structure and the C-terminal sequence appeared as an extremely conserved motif. The expression profile of the NOX genes using RNA-Seq data displayed different expression patterns across different tissues,including leaves,mid vein,petiole,stem,lateral bud, shoot apical meristem(SAM), root apical meristem(RAM), friable embryogenic callus(FEC), organized embryogenic structures(OES) and storage roots. Some of them were highly enriched in storage root,for example? Manes.14G042600 and Manes.S089900,however,Manes.06G128700 and Manes.09G172500 were highly expressed in fibrous root and RAM. Notably,in the leaves of varieties cv. KU50 and Xinxuan 048,Manes.14G042600 was upregulated under drought stress,implying a potential role of NOX gene in respon-ding to abiotic stress. 【Conclusion】Cassava NOX gene family members contain conserved gene structure and their encoded proteins also contain conserved functional domains,particularly in their C-terminal regions,suggesting the underlying co-evolutionary relationship between FRO protein and NOX protein. The NOX genes are tissue specific,and some of them are responsive to abiotic stress.

    Key words: cassava; reduced nicotinamide adenine dinucleotide phosphate oxidase(NOX); gene family members;bioinformatics; drought stress

    0 引言

    【研究意義】木薯(Manihot esculenta Crantz)是一種重要的塊根糧食作物,主要種植于熱帶和亞熱帶地區(qū),其塊莖中淀粉含量較高(70%~90%),約是水稻光合速率的2倍(Ihemere et al.,2006;Santa Brígida et al.,2014;Oyelakin et al.,2015)。研究發(fā)現(xiàn),塊莖的發(fā)育和儲(chǔ)藏根采后生理退化與活性氧(Reactive oxygen species,ROS)的產(chǎn)生密切相關(guān),活性氧能誘導(dǎo)超氧化物離子引起非生物脅迫,在響應(yīng)多種與植物氧化還原反應(yīng)相關(guān)的生物學(xué)過(guò)程中發(fā)揮重要作用(Xu et al.,2013;Liu et al.,2017);而還原型煙酰胺腺嘌呤二核苷酸磷酸氧化酶(NOX)是生成活性氧的主要催化酶(Bedard and Krause,2007;Lambeth and Neish,2014)。因此,利用生物信息學(xué)方法研究木薯NOX基因家族,揭示其成員數(shù)量、理化特性、染色體位置分布、基因結(jié)構(gòu)、進(jìn)化關(guān)系、蛋白三維結(jié)構(gòu)和時(shí)空表達(dá)特征等,對(duì)研究木薯塊莖的發(fā)育生物學(xué)特性具有重要意義?!厩叭搜芯窟M(jìn)展】NOX普遍存在于植物、動(dòng)物和真菌中(Bedard and Krause,2007)。在動(dòng)物細(xì)胞中,NOX是一種面向細(xì)胞外空間的酶復(fù)合物,包含膜結(jié)合異質(zhì)二聚體即黃細(xì)胞色素b558(由gp91phox和p22phox組成)及吞噬細(xì)胞的4種胞質(zhì)成分即p47phox、p67phox、p40phox和小分子蛋白GTPase Rac(Vignais,2002;Quinn and Gauss,2004)。當(dāng)細(xì)胞受到刺激后,胞質(zhì)成分與黃酮細(xì)胞色素b558相互作用,激活NOX,催化發(fā)生化學(xué)反應(yīng),最終產(chǎn)生過(guò)氧化氫(H2O2)和羥基自由基(OH)等活性氧(Vignais,2002;Quinn and Gauss,2004;Mizrahi et al.,2006)。在植物細(xì)胞中,NOX被稱為呼吸爆發(fā)氧化酶同源物(Rbohs),也是哺乳動(dòng)物吞噬細(xì)胞gp91phox的同源物,調(diào)控活性氧的生成(Sagi and Fluhr,2006)?;钚匝踉谥参锷L(zhǎng)調(diào)控和次生代謝反應(yīng)中發(fā)揮重要作用,不僅破壞由氧化代謝產(chǎn)生有害因子的細(xì)胞并導(dǎo)致其程序性細(xì)胞死亡,還作為植物的信號(hào)分子響應(yīng)非生物脅迫(Mittler et al.,2011;Dietz et al.,2016)。近年來(lái),已有大量研究表明,NOX基因參與植物的發(fā)育調(diào)控和逆境脅迫反應(yīng)。在水稻中,已有9個(gè)NOX基因家族成員被分離鑒定出來(lái),其中,OsNOX1、OsNOX2和OsNOX3基因在干旱和鹽脅迫下過(guò)表達(dá)能明顯提高水稻的抗逆性和耐受性(Nagano et al.,2016);OsNOX9基因高效表達(dá)能促進(jìn)根的形成(Yamauchi et al.,2017)。在擬南芥中,NOX基因家族成員AtRbohD在低氧脅迫下通過(guò)H2O2介導(dǎo)的乙烯信號(hào)釋放途徑而抑制根系伸長(zhǎng)(Kaya et al.,2014),在葉片中,AtRbohD基因通過(guò)響應(yīng)脫落酸激素的信號(hào)傳導(dǎo)而控制氣孔開(kāi)合(Chaouch et al.,2012)。在小麥中,NOX基因家族成員在不同組織中均有特異表達(dá),其中,TaNOX3和TaNOX13基因在小麥葉片中高效表達(dá),TaNOX11基因主要在莖中表達(dá),TaNOX12基因則主要在幼穗中表達(dá),TaNOX8基因在冷害和熱害脅迫下均呈下調(diào)表達(dá)(李日和陳坤明,2015;Hu et al.,2018),說(shuō)明這些基因在調(diào)控小麥的生長(zhǎng)發(fā)育和抗逆反應(yīng)方面發(fā)揮重要作用。在番茄中,SlRboh1基因高表達(dá)促使活性氧積累,以適應(yīng)更高的二氧化碳濃度和油菜素類固醇含量(Xia et al.,2014)。在草莓中,NOX基因過(guò)表達(dá)可明顯提高植株的冷害耐受性,抑制NOX基因表達(dá)會(huì)減少活性氧生成量導(dǎo)致細(xì)胞程序性死亡,耐受性降低(Zhang et al.,2018)。【本研究切入點(diǎn)】迄今,有關(guān)NOX基因家族的研究主要集中于溫帶禾本科植物或模式植物等,而鮮見(jiàn)針對(duì)木薯NOX基因家族的生物學(xué)信息及時(shí)空表達(dá)特征的研究報(bào)道,對(duì)不同逆境條件下該基因家族的分子調(diào)控機(jī)制尚不清楚。【擬解決的關(guān)鍵問(wèn)題】從木薯基因組中鑒定出NOX基因家族成員,分析其基因結(jié)構(gòu)、系統(tǒng)發(fā)育進(jìn)化關(guān)系、染色體位置分布、編碼蛋白的保守結(jié)構(gòu)域和三維結(jié)構(gòu)及不同組織和非生物脅迫下的表達(dá)情況,為深入研究該家族基因?qū)δ臼韷K根發(fā)育的調(diào)控機(jī)理及提高木薯抗逆性提供理論參考。

    1 材料與方法

    1. 1 木薯NOX基因家族成員檢索及鑒定

    從JGI數(shù)據(jù)庫(kù)(https://genome.jgi.doe.gov/)下載木薯的基因組序列,包括編碼區(qū)(CDS)序列和基因組GFF3注釋文件。利用pfam數(shù)據(jù)庫(kù)下載NOX蛋白結(jié)構(gòu)域相關(guān)的4個(gè)HMM文件,即NADPH_Ox、Ferric_reduct、NAD_Binding6和FAD_Binding8,并使用HMMER程序篩選NOX基因家族成員,E閾值設(shè)為10-E5,并在CDD數(shù)據(jù)庫(kù)進(jìn)行驗(yàn)證(Eddy,2011;Marchler-Bauer et al.,2017)。

    1. 2 木薯NOX基因家族成員基因結(jié)構(gòu)分析及染色體定位

    將NOX基因cDNA序列與木薯基因組數(shù)據(jù)庫(kù)(http://www.phyzome.net/cassava)中相應(yīng)的基因組DNA序列進(jìn)行比對(duì),分析NOX基因中內(nèi)含子和外顯子的結(jié)構(gòu)特點(diǎn),并利用GSDS對(duì)基因結(jié)構(gòu)進(jìn)行可視化(Guo et al.,2007)。從木薯基因組數(shù)據(jù)中獲得NOX基因的染色體位置和每條染色體的總長(zhǎng)度,利用MapGene2Chrom(http://mg2c.iask.in/mg2c_v2.0/)在線工具繪制NOX基因的染色體位置分布圖。

    1. 3 蛋白理化性質(zhì)及保守結(jié)構(gòu)域預(yù)測(cè)

    利用ExPASy預(yù)測(cè)NOX基因家族成員編碼蛋白的理化性質(zhì)(Gasteiger et al.,2003);使用HMMER在線預(yù)測(cè)蛋白的跨膜結(jié)構(gòu)域;利用MEME在線預(yù)測(cè)蛋白的保守基序,基序最大數(shù)目設(shè)置為10,基序長(zhǎng)度設(shè)為6~50個(gè)氨基酸(Bailey et al.,2009),并將預(yù)測(cè)的保守基序和蛋白結(jié)構(gòu)域進(jìn)行對(duì)應(yīng)分析。

    1. 4 蛋白三維結(jié)構(gòu)預(yù)測(cè)

    利用SWISS-MODEL預(yù)測(cè)木薯NOX蛋白的三維結(jié)構(gòu)。同源建模模板為NOX5蛋白(Protein Databank ID 5o0x.1.A),并利用PyMOL顯示NOX蛋白脫氫酶結(jié)構(gòu)域的保守末端。

    1. 5 系統(tǒng)發(fā)育進(jìn)化樹(shù)構(gòu)建

    利用ClusterW進(jìn)行多重序列比對(duì),采用MEGA 7.0以最大似然法(Maximum-Likelihood)構(gòu)建系統(tǒng)發(fā)育進(jìn)化樹(shù),并利用自舉重采樣法,設(shè)置1000次自舉重復(fù)對(duì)分支的可靠性進(jìn)行評(píng)估(Kumar et al.,2009)。

    1. 6 木薯轉(zhuǎn)錄組測(cè)序(RNA-Seq)數(shù)據(jù)獲取及表達(dá)特性分析

    從GEO數(shù)據(jù)庫(kù)下載11個(gè)木薯組織的RNA-Seq數(shù)據(jù)(GSE82279)及干旱脅迫下葉片組織的RNA-Seq數(shù)據(jù)(GSE98537)(Wislon et al.,2017)。利用R語(yǔ)言,基于FPKM值繪制基因表達(dá)譜,以Log2Fold Change(Log2FC)值繪制差異表達(dá)熱圖。

    2 結(jié)果與分析

    2. 1 木薯NOX基因家族成員鑒定及基因結(jié)構(gòu)分析結(jié)果

    從木薯基因組中共鑒定獲得15個(gè)NOX基因家族成員,包括9個(gè)NOX基因和6個(gè)FRO基因(NOX祖先基因),根據(jù)基因外顯子位置分布可將其明顯劃分為兩個(gè)類群(NOX基因和FRO基因),如圖1所示。FRO基因的外顯子個(gè)數(shù)明顯少于NOX基因,其中Manes.06G16800基因最少,僅含有6個(gè)外顯子,其他FRO基因含7~9個(gè),且相對(duì)于NOX基因,F(xiàn)RO基因的內(nèi)含子序列較短。NOX基因的外顯子數(shù)為11~14個(gè),平均為13個(gè),但Manes.02G016200、Manes.08G057200和Manes.S089900基因的內(nèi)含子區(qū)域較其他基因長(zhǎng),且NOX基因的第一個(gè)外顯子相對(duì)保守,推測(cè)是其編碼蛋白重要的功能域。

    2. 2 木薯NOX基因家族成員的染色體位置分布情況

    利用MapGene2Chrom繪制木薯NOX基因家族成員在木薯染色體上的位置,結(jié)果如圖2所示。這些基因分布在木薯的9條染色體和1個(gè)結(jié)構(gòu)支架上,且均被定位在染色體末端,除6號(hào)染色體、8號(hào)染色體和14號(hào)染色體外,其他染色體均只有1個(gè)基因。其中,6號(hào)染色體上分布的NOX基因和FRO基因最多,且NOX基因Manes.06G128700與2個(gè)FRO基因Manes.06G109600和Manes.06G151600的距離較近;14號(hào)染色體上分布的NOX基因Manes.14G042600與FRO基因Manes. 14G023200相距也較近(圖2),說(shuō)明上述FRO基因和NOX基因可能存在連鎖關(guān)系。

    2. 3 木薯NOX家族蛋白的理化性質(zhì)及跨膜結(jié)構(gòu)域預(yù)測(cè)結(jié)果

    木薯NOX家族蛋白(NOX蛋白和FRO蛋白,下同)的理化性質(zhì)如表1所示。NOX蛋白的氨基酸殘基數(shù)757~950個(gè),平均為897個(gè);相對(duì)分子量87.2~107.9 kD;平均為101.9 kD;理論等電點(diǎn)(pI)8.90~9.74,平均為9.49。FRO蛋白的氨基酸殘基數(shù)556~722個(gè),平均為649;相對(duì)分子量63.3~81.6 kD,平均為73.1;pI 7.72~10.06,平均為9.07??梢?jiàn),木薯NOX蛋白的分子量和氨基酸數(shù)明顯高于FRO蛋白。FRO蛋白包含有6~9個(gè)保守跨膜結(jié)構(gòu)域,而NOX蛋白均只含4個(gè)跨膜螺旋區(qū),且均位于第350~750位氨基酸殘基,說(shuō)明NOX蛋白氨基酸序列極度保守,以保證NOX蛋白的跨膜功能。

    2. 4 木薯NOX家族蛋白的結(jié)構(gòu)域分析結(jié)果

    利用MEME共鑒定出木薯NOX家族蛋白的10個(gè)保守結(jié)構(gòu)域(Motif1~Motif10),如圖3和圖4所示。其中,Motif6屬于EF_hand結(jié)構(gòu)域,是NOX蛋白特有的保守序列;Motif5、Motif9和Motif10屬于Ferric_reductase結(jié)構(gòu)域;Motif2屬于FAD_Binding8結(jié)構(gòu)域;Motif1、Motif4、Motif7和Motif8屬于NAD_Binding6結(jié)構(gòu)域。Motif1和Motif8存在于所有的NOX和FRO蛋白中,表明這些保守結(jié)構(gòu)域是連接FRO和NOX蛋白的橋梁,是進(jìn)化過(guò)程中保留下來(lái)的重要結(jié)構(gòu)。在FRO蛋白中,Manes.06G151600.3.p和Manes.06G151600.4.p缺失Motif2,可能是Manes. 06G151600基因選擇性剪接的結(jié)果,且只有Manes. 01G254300.1.p含有Motif9,而Motif9是NOX蛋白特有的保守序列,說(shuō)明FRO和NOX基因間存在進(jìn)化關(guān)系;Motif3只獨(dú)立存在于FRO蛋白中,說(shuō)明其是FRO蛋白獨(dú)立行使功能的必要位點(diǎn)。在NOX蛋白中,各結(jié)構(gòu)域相對(duì)穩(wěn)定,只有Manes.04G28100.1.p的N端重復(fù)1個(gè)Motif10。

    2. 5 系統(tǒng)發(fā)育進(jìn)化分析結(jié)果

    由圖5可知,不同物種的NOX家族蛋白可分為兩大分支(NOXs和FROs),其中,NOXs分支是含NOX蛋白結(jié)構(gòu)域的蛋白,F(xiàn)ROs分支則只含有FRO蛋白。NOXs分支包括5個(gè)苔蘚植物小立碗蘚(Physcomitrella patens)NOX蛋白、9個(gè)木薯NOX蛋白和10個(gè)擬南芥NOX蛋白,說(shuō)明隨著物種的進(jìn)化,NOX蛋白逐漸增加,其中,Manes.09G172500.1.P與擬南芥AT5G47910.1、Manes.08G118500.1.P與擬南芥AT1G09090.1、Manes.16G092100.1.P與擬南芥AT1G19230.1親緣關(guān)系較近,推測(cè)其具有相似的功能;5個(gè)苔蘚植物小立碗蘚NOX蛋白獨(dú)立聚為一個(gè)分支,僅存在于低等苔蘚植物中,可能是由于環(huán)境和生物的選擇壓力,植物在從水生到陸地進(jìn)化過(guò)程中NOX基因家族經(jīng)歷了較劇烈的演化過(guò)程。FROs分支又可被分為2個(gè)小分支(FROs1和FROs2),可能對(duì)應(yīng)2個(gè)不同的FRO蛋白亞型。

    2. 6 木薯NOX蛋白的三維結(jié)構(gòu)預(yù)測(cè)結(jié)果

    利用SWISS-MODEL構(gòu)建木薯NOX蛋白的三維結(jié)構(gòu)模型(圖6)。木薯NOX蛋白三維模型與模板的相似度為0.32~0.37,且不同NOX蛋白的三維結(jié)構(gòu)也較相似,說(shuō)明NOX蛋白的三維結(jié)構(gòu)預(yù)測(cè)結(jié)果較可靠。木薯NOX蛋白的脫氫酶結(jié)構(gòu)域明顯折疊成催化域裂縫結(jié)構(gòu),且C端保守結(jié)構(gòu)可能參與其中,這種結(jié)構(gòu)多以FDFHKENF出現(xiàn),推測(cè)其對(duì)FAD和NADPH的結(jié)合和反應(yīng)發(fā)揮重要作用(Magnani et al.,2017)。

    2. 7 木薯NOX基因家族成員的表達(dá)譜分析結(jié)果

    由圖7知,NOX基因家族成員在不同木薯組織中均有差異表達(dá)。其中,Manes.06G109600和Manes. 06G151600在葉和中脈中存在特異性表達(dá),Manes. 06G128700和Manes.09G172500在須根和根尖分生組織中均呈高表達(dá),Manes.08G118500和Manes.S089900在脆性胚性愈傷組織中也存在特異性表達(dá)。值得注意的是,Manes.14G042600、Manes.S089900、Manes.01G254300和Manes.06G168000在葉、中脈、葉柄和莖中均呈高表達(dá),而在木薯的重要器官儲(chǔ)藏根中,僅Manes.14G042600和Manes.S089900高表達(dá),推測(cè)二者是活性氧信號(hào)通路中的負(fù)向調(diào)控基因,其高表達(dá)有利于延緩木薯儲(chǔ)藏根的劣化(Zidenga et al.,2012;Xu et al.,2013)。

    2. 8 干旱脅迫條件下木薯葉中NOX基因家族成員差異表達(dá)分析結(jié)果

    由于木薯NOX和FRO基因在葉中有明顯表達(dá),選取兩個(gè)代表性的木薯品種KU50和新選048,探究木薯葉中NOX基因在干旱脅迫下的差異表達(dá),結(jié)果如圖8所示。在KU50的9個(gè)NOX基因中有3個(gè)在干旱脅迫后存在明顯的差異表達(dá)(基于log2FC>1折疊變化),即Manes.06G128700下調(diào)表達(dá),Manes. 08G057200和Manes.14G042600則明顯上調(diào)表達(dá)。在新選048的9個(gè)NOX基因,僅Manes.02G016200、Manes.09G172500、Manes.14G023200、Manes.14G0-42600和Manes.16G092100在葉中上調(diào)表達(dá),表明其在響應(yīng)干旱脅迫下發(fā)揮作用。另外,Manes.14G 042600在KU50和新選048中均呈明顯上調(diào)表達(dá)。

    3 討論

    至今,已報(bào)道了若干個(gè)不同類型的NOX家族蛋白,包括早期的NOX蛋白雛形鐵還原酶(FREs)或鐵還原氧化酶(FROs)、NOX1~NOX5、DUOX1和DUOX2,其中,NOX5為植物所特有(Bedard et al.,2007)。在植物中,NOX被稱為Rbohs,主要結(jié)構(gòu)包括6個(gè)保守的跨膜結(jié)構(gòu)域,即C端的FAD和NADPH結(jié)合結(jié)構(gòu)域、2個(gè)血紅素輔基基團(tuán)和2個(gè)N端Ca2+結(jié)合的EF手性基序,根據(jù)這些結(jié)構(gòu)基序可知Ca2+的調(diào)控方式(Geiszt,2006;Sagi and Fluhr,2006)。本研究發(fā)現(xiàn),木薯NOX蛋白C端具有保守的氨基酸殘基,與上述前人結(jié)果一致,但其氨基酸殘基基序FDFHKENF不同于其他物種,推測(cè)其為木薯所特有。

    前人研究表明,早期藻類植物基因組中只含極少數(shù)的FRO基因,不含NOX基因,如團(tuán)藻和衣藻中分別含有1個(gè)和3個(gè)FRO基因(Chang et al.,2016);苔蘚植物小立碗蘚基因組含有5個(gè)NOX基因和6個(gè)FRO基因,卷柏中含有9個(gè)NOX基因和8個(gè)FRO基因(Chang et al.,2016);擬南芥基因組中含有10個(gè)NOX基因和8個(gè)FRO基因;水稻基因組含有9個(gè)NOX基因和2個(gè)FRO基因,玉米基因組含有14個(gè)NOX基因和2個(gè)FRO基因,小麥基因組中含有36個(gè)NOX基因和10個(gè)FRO基因(Hu et al.,2018)。本研究發(fā)現(xiàn)木薯基因組含有15個(gè)NOX基因家族成員,包括9個(gè)NOX基因和6個(gè)FRO基因??梢?jiàn),從數(shù)量上來(lái)看,隨著生物的進(jìn)化,NOX蛋白從無(wú)到有,數(shù)量逐漸增多,所占比例也逐漸增大。本研究還發(fā)現(xiàn),木薯NOX蛋白的平均相對(duì)分子量和氨基酸數(shù)明顯高于FRO蛋白,與Hu等(2018)對(duì)小麥NOX基因家族成員的研究結(jié)果一致,推測(cè)在進(jìn)化過(guò)程中FRO基因發(fā)生基因擴(kuò)張和重排,逐漸進(jìn)化為NOX基因。在基因結(jié)構(gòu)方面,擬南芥、水稻和大麥的多數(shù)NOX基因含有10~14個(gè)外顯子,僅有少數(shù)NOX基因如擬南芥的AtRbohD含有8個(gè)外顯子(Kaya et al.,2014)。本研究中,木薯NOX基因的外顯子數(shù)為11~14個(gè),平均為13個(gè),與水稻、擬南芥和大麥基本一致,表明NOX基因結(jié)構(gòu)在進(jìn)化過(guò)程中保持相對(duì)穩(wěn)定,其原因可能是木薯與擬南芥、水稻和大麥等物種處在同一進(jìn)化梯度。NOX家族蛋白均屬于跨膜蛋白,大多數(shù)物種如水稻和草莓等含6個(gè)跨膜結(jié)構(gòu)域(Sagi and Fluhr,2006;Zhang et al.,2018)。但本研究發(fā)現(xiàn)木薯NOX蛋白均只含4個(gè)跨膜結(jié)構(gòu)域,而FRO蛋白含6~9個(gè)跨膜結(jié)構(gòu)域,少于水稻和草莓。另外,F(xiàn)RO作為NOX的祖先蛋白,還保留其自有的結(jié)構(gòu)功能域,如Motif3是木薯FRO蛋白所特有,不存在于NOX蛋白中,由此推測(cè)木薯FRO蛋白除進(jìn)化為NOX蛋白外還具有其他功能,有待進(jìn)一步挖掘。

    在水稻中過(guò)表達(dá)OsNOX9基因能介導(dǎo)根的形成(Yamauchi et al.,2017),在玉米中ZmNOX13和ZmNOX7基因參與調(diào)控根發(fā)育(Rajhi et al.,2011)。因此,推測(cè)木薯儲(chǔ)藏根中高表達(dá)的基因Manes.14G04- 2600和Manes.S089900具有類似的功能,參與調(diào)控木薯儲(chǔ)藏根的發(fā)育。在擬南芥葉片中,AtRbohD基因通過(guò)響應(yīng)脫落酸激素的信號(hào)傳導(dǎo)而控制氣孔開(kāi)合(Chaouch et al.,2012),推測(cè)木薯葉片中差異表達(dá)的NOX基因Manes.02G016200、Manes.06G128700、Manes.08G057200、Manes.09G172500、Manes.14G0-42600和Manes.16G092100與AtRbohD基因具有類似功能,參與到活性氧的信號(hào)通路中。由于植物在生長(zhǎng)發(fā)育過(guò)程中面臨各種環(huán)境脅迫,包括生物性脅迫和非生物脅迫(Mittler et al.,2011;Dietz et al.,2016),而活性氧作為重要的信號(hào)分子參與組織間及細(xì)胞間的信息運(yùn)輸,是響應(yīng)脅迫激發(fā)生理變化的重要物質(zhì)(Liang et al.,2018)。至今尚未清楚木薯NOX基因家族成員的功能及調(diào)控機(jī)制,為此,本研究采用生物信息學(xué)方法分析木薯NOX基因家族成員的特征,發(fā)現(xiàn)木薯NOX基因家族成員均具有其特定的表達(dá)模式,有利于今后深入研究該家族成員對(duì)木薯塊根發(fā)育的調(diào)控功能及提高植株抗非生物脅迫的能力。此外,在后續(xù)研究中,將對(duì)NOX基因家族成員進(jìn)行轉(zhuǎn)基因試驗(yàn),分析其對(duì)木薯塊根發(fā)育的具體表型效應(yīng)影響。

    4 結(jié)論

    木薯NOX基因家族成員具有保守基因結(jié)構(gòu),其編碼蛋白具有保守的功能域,尤其是C末端氨基酸序列;NOX基因的表達(dá)具有組織特異性,且部分成員的表達(dá)受非生物脅迫影響。FRO蛋白和NOX蛋白具有較明顯的共進(jìn)化關(guān)系。

    參考文獻(xiàn):

    李日,陳坤明. 2015. 小麥NOX家族基因鑒定及TaNOX12基因的克隆[D]. 楊凌:西北農(nóng)林科技大學(xué). [Li R,Chen K M. 2015. Identification of NOX gene family and cloning of TaNOX12 gene in wheat[D]. Yangling:Northwest A & F University.]

    Bailey T L,Boden M,Buske F A,F(xiàn)rith M,Grant C E,Cle-menti L,Ren J,Li W W,Noble W S. 2009. Meme suite:Tools for motif discovery and searching[J]. Nucleic Acids Research,37(2):202-208.

    Bedard K,Krause K H. 2007. The NOX family of ROS-gene-rating NADPH oxidases:Physiology and pathophysiology[J]. Physiological Reviews,87(1):245-313.

    Bedard K,Lardy B,Krause K H. 2007. NOX family NADPH oxidases:Not just in mammals[J]. Biochimie,89(9):1107-1112.

    Chang Y L,Li W Y,Miao H,Yang S Q,Li R,Wang X,Li W Q,Chen K M. 2016. Comprehensive genomic analysis and expression profiling of the NOX gene families under abiotic stresses and hormones in plants[J]. Genome Bio-logy and Evolution,8(3):791-810.

    Chaouch S,Queval G,Noctor G. 2012. AtRbohF is a crucial modulator of defence-associated metabolism and a key actor in the interplay between intracellular oxidative stress and pathogenesis responses in Arabidopsis[J]. The Plant Journal,69(4):613-627.

    Dietz K J,Mittler R,Noctor G. 2016. Recent progress in understanding the role of reactive oxygen species in plant cell signaling[J]. Plant Physiology,171(3):1535-1539.

    Eddy S R. 2011. Accelerated profile HMM searches[J]. PLoS Computational Biology,7(10):e1002195.

    Gasteiger E,Gattiker A,Hoogland C,Ivanyi I,Appel R D,Bairoch A. 2003. Expasy:The proteomics server for in-depth protein knowledge and analysis[J]. Nucleic Acids Research,31(13):3784-3788.

    Geiszt M. 2006. NADPH oxidases:New kids on the block[J]. Cardiovascular Research,71(2):289-299.

    Guo A Y,Zhu Q H,Chen X,Luo J C. 2007. GSDS:A gene structure display server[J]. Hereditas,29(8):1023-1026.

    Hu C H,Wei X Y,Yuan B,Yao L B,Ma T T,Zhang P P,Wang X,Wang P Q,Liu W T,Li W Q,Meng L S,Chen K M. 2018. Genome-wide identification and functional analysis of NADPH oxidase family genes in wheat during development and environmental stress responses[J]. Frontiers in Plant Science,9:906.

    Ihemere U,Arias G D,Lawrence S,Sayre R. 2006. Genetic modification of cassava for enhanced starch production[J]. Plant Biotechnology,4(4):453-465.

    Kaya H,Nakajima R,Iwano M,Kanaoka M M,Kimura S,Takeda S,Kawarazaki T,Senzaki E,Hamamura Y,Higashiyama T,Takayama S,Abe M,Kuchitsu K. 2014. Ca2+-activated reactive oxygen species production by Arabidopsis RbohH and RbohJ is essential for proper pollen tube tip growth[J]. Plant Cell,26(3):1069-1080.

    Kumar S,Stecher G,Tamura K. 2009. Mega7:Molecular evolutionary genetics analysis version 7.0 for bigger datasets[J]. Molecular. Biology and Evolution,33(7):1870-1874.

    Lambeth J D,Neish A S. 2014. Nox enzymes and new thin-king on reactive oxygen:A double-edged sword revisited[J]. Annual Review of Pathology,9:119-145.

    Liang D C. 2018. A salutary role of reactive oxygen species in intercellular tunnel-mediated communication[J]. Frontiers in Cell and Development Biology,6(2):1-24.

    Liu B L,Zhao S,Tan F,Zhao H,Wang D D,Si H J,Chen Q. 2017. Changes in ROS production and antioxidant capa-city during tuber sprouting in potato[J]. Food Chemistry,237(15):205-213.

    Magnani F,Nenci S,Millana Fananas E,Ceccon M,Romero E,F(xiàn)raaije M W,Mattevi A. 2017. Crystal structures and atomic model of NADPH oxidase[J]. Proceeding of the National Academic of Science,114(26):6764-6769.

    Marchler-Bauer A,Bo Y,Han L,He J,Lanczycki C J,Lu S,Chitsaz F,Derbyshire M K,Geer R C,Gonzales N R,Gwadz M,Hurwitz D I,Lu F,Marchler G H,Song J S,Thanki N,Wang Z,Yamashita R A,Zhang D,Zheng C,Geer L Y,Bryant S H. 2017. CDD/SPARCLE:Functional classification of proteins via subfamily domain architectures[J]. Nucleic Acids Research,45(1):200-203.

    Mittler R,Vanderauwera S,Suzuki N,Miller G,Tognetti V B,Vandepoele,Klaas V,Gollery M,Shulaev V,Breusegem F V. 2011. ROS signaling:The new wave?[J]. Trends in Plant Science,16(6):300-309.

    Mizrahi A,Berdichevsky Y,Ugolev Y,Molshanski-Mor S,Nakash Y,Dahan I,Alloul N,Gorzalczany Y,Sarfstein R,Hirshberg M,Pick E. 2006. Assembly of the phagocyte NADPH oxidase complex chimeric constructs derived from the cytosolic components as tools for exploring structure-function relationships[J]. Journal of Leukocyte Biology,79(5):881-895.

    Nagano M,Ishikawa T,F(xiàn)ujiwara M,F(xiàn)ukao Y,Kawano Y,Kawai-Yamada M,Shimamoto K. 2016. Plasma membrane microdomains are essential for Rac1-RbohB/H-mediated immunity in rice[J]. Plant Cell,28(8):1966-1983.

    Oyelakin O O,Opabode J T,Raji A A,Ingelbrecht I L. 2015. A cassava vein mosaic virus promoter cassette induces high and stable gene expression in clonally propagated transgenic cassava(Manihot esculenta crantz)[J]. South African Journal of Botany,97(5):184-190.

    Quinn M T,Gauss K A. 2004. Structure and regulation of the neutrophil respiratory burst oxidase:Comparison with nonphagocyte oxidases[J]. Journal of Leukocyte Biology,76(4):760-781.

    Rajhi I,Yamauchi T,Takahashi H,Nishiuchi S,Shiono K,Watanabe R,Mliki A,Nagamura Y,Tsutsumi N,Nishizawa N K,Nakazono M. 2011. Identification of genes expressed in maize root cortical cells during lysigenous aere-nchyma formation using laser microdissection and microa-rray analyses[J]. New Phytologist,190(2),351-368.

    Sagi M,F(xiàn)luhr R. 2006. Production of reactive oxygen species by plant NADPH oxidases[J]. Plant Physiology,141(2),336-340.

    Santa Brígida A B,dos Reis S P,Costa C N,Cardoso C M,Lima A M,de Souza C R. 2014. Molecular cloning and characterization of a cassava translationally controlled tumor protein gene potentially related to salt stress response[J]. Molecular Biology Reports,41(3):1787-1797.

    Vignais P V. 2002. The superoxide-generating NADPH oxidase:Structural aspects and activation mechanism[J]. Ce-llular Molecular Life Sciences,59(9):1428-1459.

    Wilson M C,Mutka A M,Hummel A W,Berry J,Chauhan R D,Vijayaraghavan A,Taylor N J,Voytas D F,Chitwood D H,Bart R S. 2017. Gene expression atlas for the food security crop cassava[J]. New Phytologist,213(4):1632-1641.

    Xia X J,Gao C J,Song L X,Zhou Y H,Shi K,Yu J Q. 2014. Role of H2O2 dynamics in brassinosteroid-induced stomatal closure and opening in Solanum lycopersicum[J]. Plant Cell Environment,37(9):2036-2050.

    Xu J,Duan X,Yang J,Beeching J R,Zhang P. 2013. Enhanced reactive oxygen species scavenging by overproduction of superoxide dismutase and catalase delays postharvest physiological deterioration of cassava storage roots[J]. Plant Physiology,161(3):1517-1528.

    Yamauchi T,Yoshioka M,F(xiàn)ukazawa A,Mori H,Nishizawa N K,Tsutsumi N,Yoshioka H,Nakazono M. 2017. An NADPH oxidase RBOH functions in rice roots during lysigenous aerenchyma formation under oxygen-deficient conditions[J]. Plant Cell,29(4):775-790.

    Zhang Y T,Li Y L,He Y W,Hu W J,Zhang Y,Wang X R,Tang H R. 2018. Identification of NADPH oxidase family members associated with cold stress in strawberry[J]. Fe-deration of European Biochemical Societies Open Bio,8(4):593-605.

    Zidenga T,Leyva-Guerrero E,Moon H,Siritunga D,Sayre R. 2012. Extending cassava root shelf life via reduction of reactive oxygen species production[J]. Plant Physiology,159(4):1396-1407.

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

    猜你喜歡
    木薯外顯子結(jié)構(gòu)域
    外顯子跳躍模式中組蛋白修飾的組合模式分析
    刮木薯
    外顯子組測(cè)序助力產(chǎn)前診斷胎兒骨骼發(fā)育不良
    蛋白質(zhì)結(jié)構(gòu)域劃分方法及在線服務(wù)綜述
    外顯子組測(cè)序助力產(chǎn)前診斷胎兒骨骼發(fā)育不良
    柬埔寨拜靈木薯喜獲大豐收,市場(chǎng)價(jià)格保持穩(wěn)定
    挖木薯
    重組綠豆BBI(6-33)結(jié)構(gòu)域的抗腫瘤作用分析
    組蛋白甲基化酶Set2片段調(diào)控SET結(jié)構(gòu)域催化活性的探討
    泛素結(jié)合結(jié)構(gòu)域與泛素化信號(hào)的識(shí)別
    精品国产三级普通话版| 久久精品国产清高在天天线| 日日摸夜夜添夜夜添小说| 日本撒尿小便嘘嘘汇集6| 亚洲av电影在线进入| 搡老岳熟女国产| 欧美最新免费一区二区三区 | 日本三级黄在线观看| 18美女黄网站色大片免费观看| 久久人人爽人人爽人人片va | 亚洲av日韩精品久久久久久密| 精品人妻视频免费看| 黄色一级大片看看| 观看免费一级毛片| xxxwww97欧美| 男女视频在线观看网站免费| 久久中文看片网| 中文字幕av成人在线电影| 老司机深夜福利视频在线观看| 日韩欧美国产在线观看| 午夜激情福利司机影院| 亚洲aⅴ乱码一区二区在线播放| 亚州av有码| 国产三级在线视频| 波野结衣二区三区在线| 丁香六月欧美| 精品日产1卡2卡| 久久久国产成人免费| 夜夜爽天天搞| 99在线人妻在线中文字幕| 亚洲精品成人久久久久久| 国产黄片美女视频| 中文字幕久久专区| 麻豆av噜噜一区二区三区| 久久久久九九精品影院| 长腿黑丝高跟| 精品久久久久久,| 久久久久国产精品人妻aⅴ院| 欧美日韩瑟瑟在线播放| 亚洲第一电影网av| 亚洲成人久久爱视频| 性插视频无遮挡在线免费观看| 国产精品久久久久久精品电影| 精品一区二区三区人妻视频| 少妇熟女aⅴ在线视频| 亚洲欧美日韩高清在线视频| 欧美黄色淫秽网站| 成年人黄色毛片网站| 69人妻影院| 国产国拍精品亚洲av在线观看| 美女免费视频网站| 国产精品一区二区性色av| 日韩国内少妇激情av| 十八禁国产超污无遮挡网站| 免费看日本二区| 午夜精品久久久久久毛片777| 国产精品久久久久久人妻精品电影| 精品乱码久久久久久99久播| 丁香欧美五月| 国产精品久久久久久久电影| 亚洲一区二区三区不卡视频| 成人亚洲精品av一区二区| 狠狠狠狠99中文字幕| 国产视频一区二区在线看| netflix在线观看网站| 亚洲av成人精品一区久久| 久久精品国产99精品国产亚洲性色| 日韩av在线大香蕉| 乱人视频在线观看| 国产欧美日韩精品亚洲av| 色噜噜av男人的天堂激情| 国产精品久久久久久亚洲av鲁大| 日韩中文字幕欧美一区二区| 黄色配什么色好看| 99久久精品一区二区三区| 久久伊人香网站| 日本黄色片子视频| 性插视频无遮挡在线免费观看| 给我免费播放毛片高清在线观看| aaaaa片日本免费| 成年女人毛片免费观看观看9| 午夜精品一区二区三区免费看| 精品久久久久久久久久免费视频| 国产探花极品一区二区| 中文字幕av在线有码专区| 97超级碰碰碰精品色视频在线观看| 成人特级av手机在线观看| 偷拍熟女少妇极品色| 制服丝袜大香蕉在线| 性插视频无遮挡在线免费观看| 毛片一级片免费看久久久久 | 在线观看av片永久免费下载| 久久精品国产自在天天线| 99热只有精品国产| 日本精品一区二区三区蜜桃| 日本精品一区二区三区蜜桃| 在线免费观看的www视频| 亚洲av成人不卡在线观看播放网| 成年女人看的毛片在线观看| 好男人电影高清在线观看| 欧美日韩中文字幕国产精品一区二区三区| 啪啪无遮挡十八禁网站| 美女高潮的动态| 在线免费观看不下载黄p国产 | 国产精品久久久久久人妻精品电影| 亚洲片人在线观看| 又紧又爽又黄一区二区| 精品久久久久久久久久免费视频| 一级作爱视频免费观看| 嫩草影院精品99| 日韩欧美国产一区二区入口| 久久久久国产精品人妻aⅴ院| 精品久久久久久久久久久久久| 真人一进一出gif抽搐免费| 色综合亚洲欧美另类图片| 麻豆久久精品国产亚洲av| 脱女人内裤的视频| 国产精品伦人一区二区| 午夜免费男女啪啪视频观看 | 久久久久免费精品人妻一区二区| 校园春色视频在线观看| 桃红色精品国产亚洲av| 综合色av麻豆| 午夜福利欧美成人| 国内精品美女久久久久久| 中文字幕av在线有码专区| 级片在线观看| 国产精品久久电影中文字幕| 91麻豆av在线| 久久久久性生活片| 欧美成狂野欧美在线观看| 免费人成视频x8x8入口观看| 国产精品爽爽va在线观看网站| 亚洲,欧美精品.| 国产蜜桃级精品一区二区三区| 搡老岳熟女国产| 欧美黑人欧美精品刺激| 最近视频中文字幕2019在线8| 亚洲aⅴ乱码一区二区在线播放| 久久久国产成人免费| 91麻豆av在线| 又爽又黄无遮挡网站| 国产成人av教育| 日韩精品青青久久久久久| 国语自产精品视频在线第100页| 国产在线男女| 亚洲三级黄色毛片| 国产精品电影一区二区三区| 91九色精品人成在线观看| 亚洲最大成人中文| 国产精品国产高清国产av| 亚洲av二区三区四区| av天堂中文字幕网| 少妇的逼水好多| 黄色丝袜av网址大全| 99在线视频只有这里精品首页| 国产精品自产拍在线观看55亚洲| 国产私拍福利视频在线观看| 国产黄片美女视频| 内射极品少妇av片p| 男女那种视频在线观看| 国产真实乱freesex| 18美女黄网站色大片免费观看| 一级毛片久久久久久久久女| 免费高清视频大片| 国产精品久久久久久精品电影| 97人妻精品一区二区三区麻豆| 久久久久久九九精品二区国产| 日本成人三级电影网站| 欧美绝顶高潮抽搐喷水| 赤兔流量卡办理| 国产av一区在线观看免费| 最新在线观看一区二区三区| 最近最新免费中文字幕在线| 少妇裸体淫交视频免费看高清| 欧美性猛交黑人性爽| 我的女老师完整版在线观看| 亚洲欧美日韩无卡精品| 91字幕亚洲| 日本成人三级电影网站| 亚洲欧美日韩东京热| 真人做人爱边吃奶动态| 久久人妻av系列| 内地一区二区视频在线| 精品午夜福利视频在线观看一区| 午夜激情欧美在线| 久久久国产成人免费| 免费观看的影片在线观看| 色播亚洲综合网| 两个人视频免费观看高清| 欧美日韩乱码在线| or卡值多少钱| 欧美性猛交黑人性爽| 亚洲内射少妇av| 国产黄a三级三级三级人| АⅤ资源中文在线天堂| 亚洲综合色惰| 久久亚洲精品不卡| 国产欧美日韩一区二区三| 我的女老师完整版在线观看| 深夜精品福利| 精品久久国产蜜桃| 久久久久久大精品| 国产精品自产拍在线观看55亚洲| 深爱激情五月婷婷| 亚洲av熟女| 99国产精品一区二区蜜桃av| 日本 欧美在线| 男人的好看免费观看在线视频| 久99久视频精品免费| a在线观看视频网站| 久久午夜福利片| 精品久久久久久久末码| 国产精品嫩草影院av在线观看 | 亚洲av日韩精品久久久久久密| 免费看美女性在线毛片视频| 一个人观看的视频www高清免费观看| 日韩高清综合在线| 日韩 亚洲 欧美在线| 老司机午夜福利在线观看视频| 日本在线视频免费播放| 亚洲 欧美 日韩 在线 免费| 日韩成人在线观看一区二区三区| 久久精品综合一区二区三区| 国产欧美日韩一区二区精品| 午夜精品一区二区三区免费看| 亚洲 欧美 日韩 在线 免费| 日韩欧美一区二区三区在线观看| 久久国产乱子免费精品| 男女做爰动态图高潮gif福利片| 国产国拍精品亚洲av在线观看| 俄罗斯特黄特色一大片| 一本一本综合久久| 一本精品99久久精品77| 亚洲va日本ⅴa欧美va伊人久久| 午夜免费成人在线视频| 一二三四社区在线视频社区8| 97热精品久久久久久| 99热这里只有是精品在线观看 | av天堂在线播放| 变态另类丝袜制服| 国产精品嫩草影院av在线观看 | 一进一出抽搐动态| 级片在线观看| 亚洲色图av天堂| av在线观看视频网站免费| 美女黄网站色视频| 噜噜噜噜噜久久久久久91| 欧美激情在线99| 久久性视频一级片| 欧美丝袜亚洲另类 | 中出人妻视频一区二区| 尤物成人国产欧美一区二区三区| 中文字幕av在线有码专区| 欧美日韩国产亚洲二区| 欧美日韩综合久久久久久 | 色尼玛亚洲综合影院| 看片在线看免费视频| 亚洲成人免费电影在线观看| 美女大奶头视频| 舔av片在线| 可以在线观看的亚洲视频| 午夜久久久久精精品| 国产精品野战在线观看| 在线观看免费视频日本深夜| 亚洲成人中文字幕在线播放| 91久久精品国产一区二区成人| 亚洲第一欧美日韩一区二区三区| 久久亚洲真实| 中文字幕精品亚洲无线码一区| 免费看日本二区| 亚洲精品乱码久久久v下载方式| 日本成人三级电影网站| 久久久久久久久大av| 美女cb高潮喷水在线观看| 成人美女网站在线观看视频| 亚洲av成人精品一区久久| 日韩欧美国产一区二区入口| 亚洲精华国产精华精| 国产精品美女特级片免费视频播放器| 免费看日本二区| 成人精品一区二区免费| 成人特级黄色片久久久久久久| 亚洲精品成人久久久久久| 午夜福利免费观看在线| 亚洲成av人片免费观看| 亚洲男人的天堂狠狠| 亚洲在线观看片| 日本熟妇午夜| 神马国产精品三级电影在线观看| 欧美日韩乱码在线| 亚洲天堂国产精品一区在线| 激情在线观看视频在线高清| 99久久成人亚洲精品观看| 亚洲人成电影免费在线| 久久精品综合一区二区三区| 欧美精品国产亚洲| 亚州av有码| 免费搜索国产男女视频| 亚洲一区高清亚洲精品| 窝窝影院91人妻| 免费无遮挡裸体视频| 亚洲欧美日韩无卡精品| 午夜福利免费观看在线| 真人一进一出gif抽搐免费| 高清在线国产一区| 欧美日韩瑟瑟在线播放| 97碰自拍视频| 亚洲熟妇中文字幕五十中出| 欧美成人一区二区免费高清观看| 91麻豆av在线| 91狼人影院| 欧美又色又爽又黄视频| 久久人人爽人人爽人人片va | 直男gayav资源| 国产精品亚洲一级av第二区| 精品人妻偷拍中文字幕| 黄片小视频在线播放| 国产精品久久久久久人妻精品电影| 亚洲欧美激情综合另类| 成人欧美大片| 三级男女做爰猛烈吃奶摸视频| 成人三级黄色视频| 欧美一区二区国产精品久久精品| 一本久久中文字幕| 亚洲第一欧美日韩一区二区三区| av专区在线播放| h日本视频在线播放| 日韩欧美国产一区二区入口| 欧美精品国产亚洲| 性色avwww在线观看| 亚洲avbb在线观看| 好看av亚洲va欧美ⅴa在| 中文字幕av成人在线电影| 久久99热这里只有精品18| 嫩草影院新地址| 国产精品伦人一区二区| 51午夜福利影视在线观看| 搡老熟女国产l中国老女人| 中文字幕人成人乱码亚洲影| 高清毛片免费观看视频网站| 如何舔出高潮| 日韩中文字幕欧美一区二区| 久久久久久国产a免费观看| 欧美高清性xxxxhd video| 九九在线视频观看精品| 久久精品国产自在天天线| 亚洲黑人精品在线| 日韩欧美在线二视频| 日韩中文字幕欧美一区二区| 一级毛片久久久久久久久女| 尤物成人国产欧美一区二区三区| 亚洲一区二区三区色噜噜| 又黄又爽又免费观看的视频| 午夜免费男女啪啪视频观看 | 午夜激情欧美在线| 欧美色欧美亚洲另类二区| 热99re8久久精品国产| 亚洲欧美日韩高清在线视频| 国产一级毛片七仙女欲春2| 色在线成人网| 国产av不卡久久| 黄色视频,在线免费观看| 国产精品亚洲av一区麻豆| 午夜精品一区二区三区免费看| 亚洲 欧美 日韩 在线 免费| 精品国产三级普通话版| 欧美zozozo另类| 久久久色成人| 亚洲成av人片在线播放无| 亚洲色图av天堂| 丰满的人妻完整版| 亚洲精品影视一区二区三区av| 日日夜夜操网爽| 成人国产综合亚洲| 中文字幕av成人在线电影| 变态另类丝袜制服| 99在线视频只有这里精品首页| 日本在线视频免费播放| 全区人妻精品视频| 日本精品一区二区三区蜜桃| 亚洲人成电影免费在线| 欧美日韩瑟瑟在线播放| 一进一出好大好爽视频| 国产日本99.免费观看| netflix在线观看网站| 麻豆成人av在线观看| 人人妻,人人澡人人爽秒播| 色播亚洲综合网| 91九色精品人成在线观看| 亚洲美女视频黄频| 日本 欧美在线| 欧美潮喷喷水| av中文乱码字幕在线| 51午夜福利影视在线观看| 欧美国产日韩亚洲一区| 欧美最新免费一区二区三区 | 成年人黄色毛片网站| 18+在线观看网站| 一a级毛片在线观看| 男女之事视频高清在线观看| 欧美国产日韩亚洲一区| 日本五十路高清| 亚洲五月天丁香| 精品一区二区三区av网在线观看| 亚洲 欧美 日韩 在线 免费| 床上黄色一级片| 丰满人妻熟妇乱又伦精品不卡| 国产乱人视频| 一本一本综合久久| 在线天堂最新版资源| 老鸭窝网址在线观看| 可以在线观看毛片的网站| 美女免费视频网站| 1024手机看黄色片| 99久久九九国产精品国产免费| 欧美性猛交黑人性爽| 国产综合懂色| 亚洲成人中文字幕在线播放| 美女 人体艺术 gogo| 国产精品亚洲一级av第二区| 亚洲av熟女| 99国产精品一区二区三区| 黄片小视频在线播放| 国产成人欧美在线观看| 中文亚洲av片在线观看爽| 免费大片18禁| 国内精品久久久久精免费| 九九久久精品国产亚洲av麻豆| aaaaa片日本免费| 三级国产精品欧美在线观看| 搡老妇女老女人老熟妇| 国产精品野战在线观看| 丁香欧美五月| 亚洲人与动物交配视频| 亚洲国产高清在线一区二区三| 天美传媒精品一区二区| 亚洲一区二区三区色噜噜| 男人和女人高潮做爰伦理| 每晚都被弄得嗷嗷叫到高潮| 美女cb高潮喷水在线观看| 精品国产亚洲在线| 美女 人体艺术 gogo| 少妇人妻一区二区三区视频| 精品久久久久久久久久久久久| www.999成人在线观看| 赤兔流量卡办理| 国产黄色小视频在线观看| 黄色一级大片看看| 成年人黄色毛片网站| 日韩精品青青久久久久久| 免费在线观看亚洲国产| 一区福利在线观看| 啦啦啦观看免费观看视频高清| 丰满乱子伦码专区| 日日夜夜操网爽| 欧美性感艳星| 变态另类成人亚洲欧美熟女| 色精品久久人妻99蜜桃| ponron亚洲| 两人在一起打扑克的视频| avwww免费| 亚洲在线自拍视频| 美女xxoo啪啪120秒动态图 | 男女视频在线观看网站免费| 身体一侧抽搐| 国产精品亚洲av一区麻豆| 婷婷六月久久综合丁香| 国产精品电影一区二区三区| 性色av乱码一区二区三区2| 久久久久国产精品人妻aⅴ院| 日韩高清综合在线| 男人舔奶头视频| 国产成人aa在线观看| 国产一区二区三区视频了| aaaaa片日本免费| 一区福利在线观看| 成人一区二区视频在线观看| 一夜夜www| 午夜两性在线视频| 久久午夜亚洲精品久久| 日韩欧美国产一区二区入口| 美女高潮的动态| 99久久精品热视频| 人妻丰满熟妇av一区二区三区| 最近视频中文字幕2019在线8| 免费无遮挡裸体视频| 好男人在线观看高清免费视频| 亚洲精品一区av在线观看| 成熟少妇高潮喷水视频| 亚洲av免费在线观看| 真人做人爱边吃奶动态| 色播亚洲综合网| 亚洲色图av天堂| 中亚洲国语对白在线视频| 看十八女毛片水多多多| 国内毛片毛片毛片毛片毛片| 天堂av国产一区二区熟女人妻| 亚洲激情在线av| 亚洲第一欧美日韩一区二区三区| 免费看a级黄色片| x7x7x7水蜜桃| 日韩免费av在线播放| 熟女人妻精品中文字幕| 色哟哟哟哟哟哟| 成年人黄色毛片网站| 国产一区二区三区视频了| 不卡一级毛片| 久久久国产成人免费| 久久热精品热| 免费在线观看成人毛片| 自拍偷自拍亚洲精品老妇| 性色avwww在线观看| 精品无人区乱码1区二区| 欧美区成人在线视频| 性色avwww在线观看| 在线天堂最新版资源| 最后的刺客免费高清国语| 亚洲第一区二区三区不卡| 好男人在线观看高清免费视频| 小说图片视频综合网站| 日韩av在线大香蕉| 白带黄色成豆腐渣| 夜夜爽天天搞| 亚洲精品一卡2卡三卡4卡5卡| 日本 欧美在线| av天堂中文字幕网| 亚洲精品影视一区二区三区av| 变态另类成人亚洲欧美熟女| a在线观看视频网站| 亚洲专区国产一区二区| 国产成人影院久久av| 欧美一区二区精品小视频在线| 成人av在线播放网站| 国产免费一级a男人的天堂| 国产人妻一区二区三区在| 91午夜精品亚洲一区二区三区 | 亚洲精品久久国产高清桃花| 变态另类成人亚洲欧美熟女| 高清毛片免费观看视频网站| 亚洲一区二区三区色噜噜| 成人av在线播放网站| 波多野结衣高清无吗| 午夜激情福利司机影院| 亚洲精品亚洲一区二区| 非洲黑人性xxxx精品又粗又长| 午夜影院日韩av| 首页视频小说图片口味搜索| 国产精品一及| 中文资源天堂在线| 午夜免费激情av| 深夜a级毛片| 啦啦啦观看免费观看视频高清| 欧美成人性av电影在线观看| 精品一区二区三区人妻视频| 欧美日韩亚洲国产一区二区在线观看| 国产一区二区亚洲精品在线观看| 国产精品国产高清国产av| 可以在线观看的亚洲视频| 草草在线视频免费看| 亚洲午夜理论影院| www.www免费av| 村上凉子中文字幕在线| 高潮久久久久久久久久久不卡| 国产伦人伦偷精品视频| 99久久久亚洲精品蜜臀av| 特大巨黑吊av在线直播| 国产一区二区亚洲精品在线观看| 变态另类成人亚洲欧美熟女| 日本黄色片子视频| 久久中文看片网| 99国产精品一区二区三区| 成人av在线播放网站| 久久精品国产亚洲av香蕉五月| 9191精品国产免费久久| 18禁裸乳无遮挡免费网站照片| 日韩欧美精品v在线| 淫妇啪啪啪对白视频| 亚洲aⅴ乱码一区二区在线播放| 91午夜精品亚洲一区二区三区 | 亚洲国产高清在线一区二区三| 成年女人毛片免费观看观看9| 床上黄色一级片| 亚洲专区国产一区二区| 亚洲一区高清亚洲精品| 国产黄a三级三级三级人| 在线观看66精品国产| 国产69精品久久久久777片| 亚洲av免费在线观看| 成人国产综合亚洲| 免费在线观看成人毛片| 国产高清视频在线观看网站| 欧美激情久久久久久爽电影| 非洲黑人性xxxx精品又粗又长| 黄片小视频在线播放| 又黄又爽又免费观看的视频| 中文字幕精品亚洲无线码一区| 欧美黄色片欧美黄色片| 欧美不卡视频在线免费观看| aaaaa片日本免费| 欧美精品国产亚洲| 亚洲专区国产一区二区| 国产精品伦人一区二区| 国产亚洲欧美在线一区二区| av天堂中文字幕网| 久久草成人影院| 亚洲内射少妇av| 女同久久另类99精品国产91| 国产精品亚洲av一区麻豆| 99在线人妻在线中文字幕| 三级男女做爰猛烈吃奶摸视频| 男人狂女人下面高潮的视频| 床上黄色一级片| 亚洲中文字幕日韩| 亚洲久久久久久中文字幕|