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

    油棕等熱帶植物DXS基因的生物信息學(xué)分析

    2016-05-27 01:50:46曹紅星李東霞雷新濤
    廣西植物 2016年4期
    關(guān)鍵詞:油棕生物信息學(xué)

    石 鵬, 曹紅星, 李東霞, 王 永, 雷新濤

    ( 中國熱帶農(nóng)業(yè)科學(xué)院 椰子研究所 海南省熱帶油料作物生物學(xué)重點(diǎn)實(shí)驗(yàn)室, 海南 文昌 571339 )

    ?

    油棕等熱帶植物DXS基因的生物信息學(xué)分析

    石鵬, 曹紅星, 李東霞, 王永, 雷新濤*

    ( 中國熱帶農(nóng)業(yè)科學(xué)院 椰子研究所 海南省熱帶油料作物生物學(xué)重點(diǎn)實(shí)驗(yàn)室, 海南 文昌 571339 )

    摘要:油棕等熱帶植物含有豐富的胡蘿卜素和維生素E等類異戊二烯物質(zhì),類胡蘿卜素和甾醇等類異戊二烯物質(zhì)在植物生命活動中扮演重要角色,并且對保護(hù)人類健康具有重要意義,MEP途徑是合成類異戊二烯的重要途徑之一。DXS是MEP途徑中的第一個限速酶,其功能在油棕等熱帶植物中極其保守。為了弄清油棕等熱帶植物DXS的結(jié)構(gòu)和功能特點(diǎn),該研究利用生物信息學(xué)工具和軟件對以油棕等熱帶植物類異戊二烯合成關(guān)鍵基因DXS為對象,進(jìn)行核酸和氨基酸序列的理化性質(zhì)、蛋白質(zhì)結(jié)構(gòu)以及功能結(jié)構(gòu)域等分析,探討了不同物種間的親緣關(guān)系。結(jié)果表明:DXS基因起始密碼子均為ATG,終止密碼子則分為TAG、TAA和TGA,DXS蛋白質(zhì)屬于不具有信號肽的親水性蛋白,可能作為轉(zhuǎn)運(yùn)蛋白在葉綠體基質(zhì)中發(fā)揮作用,未發(fā)現(xiàn)明顯的跨膜結(jié)構(gòu)域,磷酸化位點(diǎn)有36個,其中絲氨酸、蘇氨酸和酪氨酸位點(diǎn)分別為17、11和8個,無規(guī)則卷曲和α-螺旋是蛋白質(zhì)二級結(jié)構(gòu)主要的結(jié)構(gòu)元件,三級結(jié)構(gòu)預(yù)測具有DXS酶特征,硫胺素焦磷酸鹽結(jié)合位點(diǎn)和PLN02582保守結(jié)構(gòu)域,不同植物DXS功能結(jié)構(gòu)域非常保守,可以作為判斷不同物種間親緣關(guān)系的重要依據(jù)。該研究結(jié)果為油棕等熱帶植物DXS的結(jié)構(gòu)、功能分析和利用提供了進(jìn)一步的信息,為其品質(zhì)性狀分子機(jī)制研究及遺傳改良奠定了基礎(chǔ)。

    關(guān)鍵詞:油棕, DXS基因, 生物信息學(xué), 系統(tǒng)進(jìn)化樹

    油棕(ElaeisguineensisJacq.)屬棕櫚科多年生喬木,是熱帶產(chǎn)油植物,果實(shí)中類胡蘿卜素的含量是決定油棕品質(zhì)的一個重要指標(biāo)。粗棕櫚油中含有600 mg·kg-1類胡蘿卜素,其中80%都是β-胡蘿卜素和α-胡蘿卜素,這也是粗棕櫚油呈黃色或紅棕色的原因(雷新濤和曹紅星,2013)。巴西橡膠、可可、白木香、杜仲和麻瘋樹都是重要的熱帶植物,對油棕等熱帶植物進(jìn)行DXS的生物信息學(xué)分析,對于了解這些重要植物的DXS蛋白質(zhì)結(jié)構(gòu)和功能具有重要意義,同時(shí)也為分析這些熱帶植物間進(jìn)化關(guān)系提供依據(jù)。

    類異戊二烯是最大的一類天然產(chǎn)物,目前為止已知的有3萬多種(Sacchettini & Poulter,1997)。一些類異戊二烯在細(xì)胞中有著重要作用,如甾醇是細(xì)胞膜結(jié)構(gòu)的重要組分(Erick,2008),長萜醇在蛋白質(zhì)糖基化中擔(dān)任糖基載體(Vincent & Dirk,2011),泛素和質(zhì)體醌在光合作用和呼吸作用中扮演電子轉(zhuǎn)運(yùn)體的角色(Keisuke et al,2010),脫落酸、細(xì)胞分裂素、赤霉素和甾醇作為激素控制著生理和發(fā)育過程(Jonathan & Neil,2001)。盡管各種類異戊二烯的結(jié)構(gòu)和功能存在差異,但是所有的都由異戊烯焦磷酸(IPP)和它的同分異構(gòu)體二甲基丙烯焦磷酸(DMAPP)構(gòu)成。IPP和DMAPP能通過兩個不同生物合成路線:甲羥戊酸(MVA)途徑,其主要存在真菌、酵母和動物細(xì)胞質(zhì)中;2-C-甲基-D-赤藻糖醇-4-磷酸(MEP)途徑,其存在原生生物,大多數(shù)細(xì)菌,綠藻和高等植物質(zhì)體中(Lombard & Moreira,2011;Wanke et al,2001;Vranova et al,2013)。MEP途徑的第一步是丙酮酸和3-磷酸甘油醛在1-脫氧木酮糖-5-磷酸合成酶(DXS)的催化下生成1-脫氧木酮糖-5-磷酸(DXP)。DXS是一種新的轉(zhuǎn)酮酶,由DXS基因翻譯而來(Elizabeth et al,2009;Han et al,2013)。在植物和大腸桿菌中,由此反應(yīng)產(chǎn)生的DXP被用來在質(zhì)體中生成IPP。DXS基因最初在大腸桿菌中被克隆和分析,其它細(xì)菌和高等植物中隨后也找到了同源基因(Cordoba et al,2011;Zhang et al,2009;Gong et al,2006;Luisa et al,1998;Kuzuyama et al,2000;Hahn et al,2001)。胡蘿卜素和葉黃素統(tǒng)稱為類胡蘿卜素,是所有進(jìn)行光合作用生物中的色素。在高等植物中,類胡蘿卜素在綠色組織光合作用中扮演重要角色。其中,β-胡蘿卜素是維生素A的前體,具有抗癌的作用,其缺乏會導(dǎo)致失明和早逝(Rao et al,2007;Shan & Li,2008;Aguila & Manuel,2012)。由于類胡蘿卜素的重要性,高等植物中編碼類胡蘿卜素生物合成的酶基因成為了增加類胡蘿卜素生物合成的潛在位點(diǎn)。有文獻(xiàn)報(bào)道證實(shí)DXS是途徑中的限速酶,MEP途徑中的基因高表達(dá)與植物中更高含量的類胡蘿卜素有直接關(guān)系(Peng,2008)。

    生物信息學(xué)興起于21世紀(jì),利用信息學(xué)工具分析核酸和蛋白質(zhì)等大分子,大規(guī)模挖掘這些大分子的生物學(xué)意義(Luscombe et al,2001)。目前,植物中開展DXS的生物信息學(xué)分析還不多,尤其在油棕等熱帶植物中還很少報(bào)道。本文主要以油棕等熱帶植物的DXS核酸和蛋白質(zhì)序列為分析對象,利用大量的生物信息學(xué)工具挖掘其序列的生物學(xué)意義,為進(jìn)一步研究油棕等熱帶植物DXS基因功能和類胡蘿卜素代謝過程奠定基礎(chǔ)。

    1材料與方法

    1.1 數(shù)據(jù)來源

    核酸和蛋白質(zhì)序列均來源于NCBI(National Center for Biotechnology Information)數(shù)據(jù)庫中已有的油棕等熱帶植物1-脫氧木酮糖-5-磷酸合成酶基因的核酸序列及其對應(yīng)的蛋白質(zhì)序列,如油棕(ElaeisguineensisJacq.,AY583783、AAS99588)、白木香(Aquilariasinensis,JX860326、AFU75321)、杜仲(Eucommiaulmoides,JX458817、AFU93068)、巴西橡膠(Heveabrasiliensis,AB294698、BAF98288)、麻瘋樹(Jatrophacurcas,KF955327、AHL17583)、可可(Theobromacacao,XM007035370、XP007035432);另外,構(gòu)建系統(tǒng)進(jìn)化樹還用到了以下幾種植物的DXS蛋白質(zhì)序列:高良姜(Alpiniaofficinarum,AEK69518)、陽春砂(Amomumvillosum,ACR02668)、擬南芥(Arabidopsisthaliana,NP193291)、巴豆(Crotonstellatopilosus,BAF75640)、海島棉(Gossypiumbarbadense,ABN13970)、苜蓿(Medicagotruncatula,XP003609567)、煙草(Nicotianatabacum,ACF60511)、粳稻(Oryzasativa,NP001055524)、葛麻姆(Puerariamontana,AAQ84169)、蓖麻(Ricinuscommunis,XP002516843)、羅漢果(Siraitiagrosvenorii,AEM42997)、番茄(Solanumlycopersicum,NP001234672)、馬鈴薯(S.tuberosum,NP001275130)、南非醉茄(Withaniasomnifera,AFI98878)、玉米(Zeamays,NP001157805)。

    表 1 生物信息學(xué)分析使用的工具和軟件

    1.2 方法

    所用生物信息學(xué)工具和軟件如表1所示,具體參數(shù)設(shè)置見結(jié)果。

    2結(jié)果與分析

    2.1 油棕等熱帶植物DXS基因核酸及氨基酸序列分析

    利用ORF Finder和ProtParam對油棕、巴西橡膠、麻瘋樹、杜仲、可可和白木香6種熱帶植物DXS基因的核酸序列及其氨基酸序列進(jìn)行分析。分析結(jié)果如表2所示,不同植物DXS基因開放讀碼框的堿基數(shù)及所編碼的氨基酸殘基數(shù)變化不大,起始密碼子均為ATG,終止密碼子分為3種,即油棕、巴西橡膠和白木香均為TAG,麻瘋樹和杜仲為TAA,可可則為TGA ;而且不同植物DXS的分子量、理論等電點(diǎn)和高含量氨基酸成分基本一致; Leu、Ala、Gly是這6種熱帶植物DXS中含量最豐富的氨基酸;根據(jù)不穩(wěn)定性指數(shù)小于40為穩(wěn)定性蛋白質(zhì)判斷,巴西橡膠、杜仲和可可的DXS為穩(wěn)定類蛋白質(zhì),而油棕、麻瘋樹和白木香的DXS為不穩(wěn)定類蛋白質(zhì);且它們的親水系數(shù)較低,均為親水性蛋白質(zhì)。

    表 2 6種熱帶植物DXS核酸及氨基酸序列分析

    2.2 油棕等熱帶植物DXS親水性/疏水性的預(yù)測和分析

    利用ProtScale對油棕氨基酸序列的疏水性/親水性進(jìn)行預(yù)測(圖1)。從圖1發(fā)現(xiàn),多肽鏈第55位的His具有最低的分值-2.656,親水性最強(qiáng);而第185位的Gly具有最高的分值2.067,疏水性最強(qiáng)。整體來看,親水氨基酸分布比較均勻,且數(shù)量大于疏水氨基酸,故推測油棕DXS蛋白質(zhì)為親水性蛋白。

    2.3 油棕等熱帶植物DXS二級結(jié)構(gòu)的預(yù)測和分析

    利用SOPMA對油棕DXS氨基酸序列進(jìn)行二級結(jié)構(gòu)預(yù)測(圖2)。從圖2可以發(fā)現(xiàn),油棕DXS由40.45%的無規(guī)則卷曲(Random coil),38.19%的α-螺旋(Alpha helix),15.28%的延伸鏈(Extended strand)和6.08%β-轉(zhuǎn)角(Beta turn)組成。在整個蛋白質(zhì)結(jié)構(gòu)中,無規(guī)則卷曲和α-螺旋是油棕DXS最主要的結(jié)構(gòu)元件,而延伸鏈和β-轉(zhuǎn)角則次之。

    2.4 油棕等熱帶植物DXS蛋白三級結(jié)構(gòu)的預(yù)測和分析

    利用SWISS-MODEL的同源建模方法預(yù)測油棕DXS蛋白的三級結(jié)構(gòu),獲得油棕DXS三維結(jié)構(gòu)模式(圖3)。同源建模分析發(fā)現(xiàn)其與蛋白質(zhì)模板2o1s.1.A有50.58%的相似性,且根據(jù)三級結(jié)構(gòu)特征推測為1-脫氧木酮糖-5-磷酸合成酶(1-deoxy-D-xylulose-5-phosphate synthase)。從圖3可以看出,紫色卷曲狀結(jié)構(gòu)代表的無規(guī)則卷曲遍布在油棕DXS三級結(jié)構(gòu)中。

    2.5 油棕等熱帶植物DXS跨膜結(jié)構(gòu)域的預(yù)測和分析

    利用TMHMM對不同熱帶植物DXS氨基酸序列的跨膜結(jié)構(gòu)域進(jìn)行預(yù)測(圖4),未發(fā)現(xiàn)明顯的跨膜結(jié)構(gòu)域。

    2.6 油棕等熱帶植物DXS翻譯后修飾的預(yù)測和分析

    利用NetPhos 2.0 Server對油棕DXS的翻譯后修飾預(yù)測的結(jié)果表明(圖5),整個多肽鏈中分值在0.5(閾值)以上的氨基酸位點(diǎn)有36個,可知油棕DXS的磷酸化位點(diǎn)有36個,其中Ser:17,Thr:11,Tyr:8。

    2.7 油棕等熱帶植物DXS的亞細(xì)胞定位分析

    利用PSORT對油棕等熱帶植物的DXS基因的編碼產(chǎn)物進(jìn)行亞細(xì)胞定位分析,結(jié)果表明,油棕和其它熱帶植物的DXS蛋白質(zhì)定位在葉綠體中基質(zhì)的可能性最大,該結(jié)果說明熱帶植物DXS蛋白可能主要是在葉綠體基質(zhì)中行使功能。

    2.8 油棕等熱帶植物DXS信號肽和導(dǎo)肽的預(yù)測及分析

    利用SignalP 4.1 Server對油棕DXS氨基酸序列的信號肽進(jìn)行預(yù)測(圖6)。圖6結(jié)果表明,油棕DXS不含有信號肽。

    同時(shí)利用TargetP1.1 Server進(jìn)行導(dǎo)肽分析也證實(shí)了其作為葉綠體轉(zhuǎn)運(yùn)蛋白(分值為0.409)的可能性比較大(圖7),該蛋白存在信號肽(分值為0.030)的可能性很小。此結(jié)果與亞細(xì)胞定位及信號肽預(yù)測結(jié)果基本一致。

    2.9 油棕等熱帶植物DXS保守結(jié)構(gòu)域的預(yù)測和分析

    通過NCBI中的保守結(jié)構(gòu)域數(shù)據(jù)庫 (Conserved

    Domain Database,CDD)工具對DXS基因的氨基酸序列分析發(fā)現(xiàn),DXS含有硫胺素焦磷酸鹽(TPP,Thiamine pyrophosphate)和嘧啶(PYR,Pyrimidine)結(jié)合位點(diǎn),TPP_DXS、TPP_PYR_DXS_TK_like和Transketolase_C保守域,TPP_enzymes、TPP_enzyme_PYR和Transketolase_C superfamily家族和PLN02582保守結(jié)構(gòu)域。其中PLN02582是1-脫氧木酮糖-5-磷酸合成酶的特征序列(圖8)。

    圖 1 油棕DXS的疏水性/親水性預(yù)測和分析Fig. 1 Hydrophobicity or hydrophilicity prediction and analysis of EgDXS

    2.10 油棕等熱帶植物DXS序列比對及系統(tǒng)進(jìn)化分析

    用Blast程序比對油棕與其它植物DXS的核酸及氨基酸序列同源性,結(jié)果表明油棕DXS的核酸序列與其它植物如高良姜(HQ874656)、麻瘋樹(KF955327)、葛麻姆(FJ455512)、橡膠(AB294698)、番茄(NM001247743)和南非醉茄(JQ710678)的一致性分別是80%、80%、79%、79%、

    圖 2 油棕DXS二級結(jié)構(gòu)預(yù)測 藍(lán)色線條,紅色線條,綠色線條和紫色線條分別表示α-螺旋,延伸鏈,β-轉(zhuǎn)角和無規(guī)則卷曲。橫坐標(biāo)和縱坐標(biāo)分別表示二級結(jié)構(gòu)在所處的位置和數(shù)量。Fig. 2 Secondary structure prediction of EgDXS Blue lines, red lines, green lines, purple lines represent alpha helix, extended strand, beta turn and random coil respectively. Horizontal and vertical axis represent position and quantity of secondary structure respectively.

    圖 3 油棕DXS三級結(jié)構(gòu)預(yù)測Fig. 3 Tertiary structure prediction of EgDXS

    圖 4 油棕DXS的跨膜結(jié)構(gòu)域預(yù)測Fig. 4 Transmembrane structure prediction of EgDXS

    圖 5 油棕DXS翻譯后修飾預(yù)測 磷酸化位點(diǎn)分為絲氨酸,蘇氨酸和酪氨酸3種。Fig. 5 Posttranslational modification prediction of EgDXSPhosphorylation sites are divided into serine, threonine and tyrosine.

    圖 6 油棕DXS信號肽預(yù)測 C-score. 切割位點(diǎn)的可能性; S-score. 氨基酸是信號肽的可能性; Y-score. 結(jié)合C-score和S-score的導(dǎo)數(shù),能更好地反映切割位點(diǎn)的可能性。Fig. 6 Signal peptide prediction of EgDXS C-score. Possibility of cleavage site; S-score. Possibility of signal peptide; Y-score. Derivative of the C-score combined with the S-score resulting in a better cleavage site prediction than the raw C-score alone.

    圖 7 油棕DXS導(dǎo)肽預(yù)測 Len. 蛋白質(zhì)長度; cTP. 作為葉綠體轉(zhuǎn)運(yùn)蛋白的可能性; mTP. 作為線粒體轉(zhuǎn)運(yùn)蛋白的可能性; SP. 信號肽的可能性。Fig. 7 Leading peptide prediction of EgDXS Len. Length of protein; cTP. Possibility of chloroplast transport protein; mTP. Possibility of mitochondira transport protein; SP. Possibility of signal pepetide.

    79%和79%;蛋白質(zhì)的Blastp比對發(fā)現(xiàn),油棕DXS的氨基酸序列與其他植物的氨基酸序列有更高的相似性,如與長春花(AGL40532)、陽春砂(ACR02668)、玉米(ACG27905)、麻瘋樹(AHL17583)、高良姜(AEK69518)和羅漢果(AEM42997)的一致性分別為86%、85%、86%、86%、86%和85%。以上結(jié)果表明,植物DXS氨基酸序列之間的相似性高于核酸序列。

    把21個來自不同植物的DXS氨基酸序列用MEGA 6.06軟件采用Neighbor-Joining法,進(jìn)行1 000次bootstrap統(tǒng)計(jì)學(xué)檢驗(yàn)構(gòu)建系統(tǒng)進(jìn)化樹(圖9)。圖9結(jié)果表明,21種不同植物的DXS氨基酸序列可以分為3大類,巴西橡膠和蓖麻等10種植物分為第Ⅰ類,其中大戟科植物麻瘋樹、巴豆、蓖麻和巴西橡膠在此類,而白木香、可可和杜仲等熱帶植物也聚在此類;豆科植物和茄科植物分在第Ⅱ類,其中同為豆科植物的葛麻姆和苜蓿聚成一簇的信心指數(shù)為100%,馬鈴薯、番茄、南非醉茄和煙草為茄科,其中馬鈴薯和番茄聚成一簇的信心指數(shù)為100%;油棕、姜科和禾本科植物分在第Ⅲ類,其中同為禾本科的粳稻和玉米聚成一簇的信心指數(shù)為100%,同為姜科的高良姜和陽春砂聚成一簇的信心指數(shù)為100%。從聚類分析結(jié)果可以發(fā)現(xiàn),用DXS氨基酸序列構(gòu)建的系統(tǒng)進(jìn)化樹可以比較真實(shí)地反映不同植物之間的進(jìn)化關(guān)系,對判斷植物間的親緣關(guān)系具有很好的參考價(jià)值。

    3討論與結(jié)論

    類胡蘿卜素對于植物來說是必需的,并且在為人類提供營養(yǎng)和保護(hù)健康方面具有重要作用。本文利用生物信息學(xué)工具和軟件對油棕等熱帶植物的DXS基因的核酸及蛋白質(zhì)理化性質(zhì)、結(jié)構(gòu)組成和功能特征等進(jìn)行了預(yù)測和分析,并構(gòu)建了油棕等21種植物的系統(tǒng)進(jìn)化樹。預(yù)測結(jié)果表明,不同植物DXS的理化性質(zhì)差異不大,蛋白質(zhì)穩(wěn)定性不一,都屬于親水性蛋白質(zhì)??缒そY(jié)構(gòu)域、信號肽和導(dǎo)肽預(yù)測結(jié)果表明,油棕等熱帶植物的DXS都沒有發(fā)現(xiàn)明顯的跨膜結(jié)構(gòu)域和信號肽,主要是在葉綠體基質(zhì)中發(fā)揮作用,可能作為葉綠體轉(zhuǎn)運(yùn)蛋白行使功能。DXS的保守結(jié)構(gòu)域分析發(fā)現(xiàn),DXS含有硫胺素焦磷酸鹽和嘧啶結(jié)合位點(diǎn), TPP_DXS、TPP_PYR_DXS_TK_like和Transketolase_C保守域,TPP_enzymes、TPP_enzyme_PYR和Transketolase_C superfamily家族以及PLN02582保守結(jié)構(gòu)域,其中PLN02582保守結(jié)構(gòu)域是1-脫氧木酮糖-5-磷酸合成酶的特征結(jié)構(gòu)域。

    圖 8 油棕DXS保守結(jié)構(gòu)域預(yù)測Fig. 8 Conserved domain prediction of EgDXS

    圖 9 油棕等21種植物DXS氨基酸序列系統(tǒng)進(jìn)化樹分析Fig. 9 Phylogenetic tree analysis of Amino acid series from 21 plants including oil palm

    在MEP途徑中,以丙酮酸和3-磷酸甘油醛為底物,先后經(jīng)過8個酶的催化最終形成了合成類異戊二烯的通用底物IPP和DMAPP。實(shí)驗(yàn)證實(shí)作為MEP途徑第一個酶DXS在IPP和DMAPP合成過程中發(fā)揮重要作用。擬南芥中對MEP途徑中的酶研究較多,目前MEP途徑中的所有酶都已經(jīng)確認(rèn)。MEP途徑的第一個反應(yīng)就是在DXP合成酶DXS的催化作用下甘油醛-3-磷酸硫胺素縮合生成1-脫氧木酮糖-5-磷酸。研究表明,通過MEP途徑合成類胡蘿卜素主要受到限速酶DXS,DXR和HDR的控制(Sawitri & Wallie,2005;Elizabeth et al,2009)。Paruicia et al(2006)認(rèn)為在充分了解植物類胡蘿卜素合成機(jī)制的基礎(chǔ)上,可以通過生物技術(shù)手段提高植物中的類胡蘿卜素,為人類提供類胡蘿卜素含量更高的食物。目前,番茄、馬鈴薯、油菜和水稻中都已經(jīng)開展了相關(guān)轉(zhuǎn)基因?qū)嶒?yàn)。植物中類胡蘿卜素的代謝和生物合成途徑已經(jīng)得到廣泛研究,合成途徑中幾乎所有編碼合成酶的基因都被分離出來。Sanae & Akemi(2006)研究菊花花瓣顏色發(fā)現(xiàn),是類胡蘿卜素的合成和積累機(jī)制不同形成了黃色和白色兩種花瓣顏色,其中LCYE,LCYB和DXS的表達(dá)量與花瓣中類胡蘿卜素的含量息息相關(guān)。在番茄果實(shí)發(fā)育過程中,DXS的表達(dá)量與果實(shí)中類胡蘿卜素含量呈正相關(guān)關(guān)系,說明DXS在番茄果實(shí)發(fā)育和器官形成過程中發(fā)揮重要作用(Luisa et al,2000)。萊茵衣藻中研究發(fā)現(xiàn),DXS和DXR等基因的表達(dá)量隨著光周期波動,DXS等基因的表達(dá)量表現(xiàn)出典型的晝夜節(jié)律模式,并且也會受到中間代謝產(chǎn)物的反饋調(diào)節(jié)(Tian et al,2010)。因此,在植物類胡蘿卜素合成代謝過程中,DXS是調(diào)控類胡蘿卜素合成的關(guān)鍵酶,對其進(jìn)行生物信息學(xué)分析為今后深入研究具有重要意義。油棕和巴西橡膠等熱帶植物DXS生物信息學(xué)分析,可以尋找這些不同科的熱帶植物在DXS上的共同特征和不同點(diǎn),希望為熱帶植物DXS結(jié)構(gòu)和功能研究奠定基礎(chǔ)。植物激素、類胡蘿卜素和甾醇等類異戊二烯物質(zhì)在熱帶植物生命活動中具有重要作用,而它們在熱帶植物中的具體功能作用及其控制基因研究較少,需要進(jìn)一步開展相關(guān)研究,揭開類異戊二烯在熱帶植物中的生物合成途徑及其調(diào)控機(jī)制。油棕等21種植物DXS氨基酸序列比對和系統(tǒng)進(jìn)化樹分析發(fā)現(xiàn),不同植物的DXS序列相似性較高,功能域非常保守,以DXS氨基酸序列構(gòu)建的進(jìn)化樹能準(zhǔn)確反映不同植物間的親緣關(guān)系,也可作為評價(jià)油棕等熱帶植物種質(zhì)資源的重要依據(jù)。

    參考文獻(xiàn):

    AGUILA RS, MANUEL RC, 2012. Carotenoid biosynthesis inArabidopsis: a colorful pathway [J]. Arabidopsis Book, 10:e0158.CORDOBA E, PORTA H, ARROYO A, et al, 2011. Functional characterization of the three genes encoding 1-deoxy-D-xylulose 5-phosphate synthase in maize [J]. J Exp Bot, 62(6):2 023-2 038.

    ERICK JD, 2008. Sterols and membrane dynamics [J]. J Chem Biol, 1(1-4):63-77.

    ELIZABETH C, MARI S, PATRICIA L, 2009. Unravelling the regulatory mechanisms that modulate the MEP pathway in higher plants [J]. J Exp Bot, 60(10):2 933-2 943.

    GONG YF, LIAO ZH, GUO BH, et al, 2006. Molecular cloning and expression profile analysis ofGinkgobilobaDXS gene encoding 1-deoxy-D-xylulose 5-phosphate synthase, the first committed enzyme of the 2-C-methyl-D-erythritol 4-phosphate pathway [J]. Plant Med, 72(4):329-335.

    HAN M, HEPPEL SC, SU T, et al, 2013. Enzyme inhibitor studies reveal complex control of meghyl-D-erythritolf-4-phosphate (MEP) pathway enzyme expression inCatharanthusroseu[J]. PLoS ONE, 8(5):e62467.

    HAHN FM, EUBANKS LM,TESTA CA, et al, 2001. 1-deoxy-D-xylulose-5-phosphate synthase, the gene product of open reading frame (ORF) 2816 and ORF 2895 inRhodobactercapsulatus[J]. J Bact, 183:1-11.JONATHAN DB, NEIL WP, 2001. Plant hormones and the control of physiological processes [J]. New Phytol, 152(3):375-407.

    KEISUKE Y, MASARU S, ICHIRO T,et al, 2010. Simultaneous determination ofinvivoplastoquinone and ubiquinone redox states by HPLC-based analysis [J]. Plant Cell Physiol, 51(5):836-841.

    KUZUYAMA T, TAKAGI M, TAKAHASHI S, et al, 2000. Cloning and characterization of 1-deoxy-D-xylulose-5-phosphate synthase fromStreptomycessp. Strain CL190, which uses both the mevalonate and nonmevalonate pathway for isopentenyl diphoshate biosynthesis [J]. J Bacteriol, 182(4):891-897.

    LEI XT,CAO HX, 2013. Oil Palm [M]. Beijing: China Agriculture Press:103-104. [雷新濤,曹紅星, 2013. 油棕 [M]. 北京: 中國農(nóng)業(yè)出版社:103-104.]

    LOMBARD J, MOREIRA D, 2011. Origins and early evolution of the mevalonate pathway of isoprenoid biosynthesis in the three domains of life [J]. Mol Biol Evol, 28(1):87-99.

    LUISA ML, NARCISO C, SURYA RP, et al, 1998. Cloning and characterization of a gene fromEscherichiacoliencoding a transketolase-like enzyme that catalyzes the synthesis of D-1-deoxyxylulose-5-phosphate, a common precuisor for isoprenoid, thiamin, and pyridoxol biosynthesis [J]. Proc Natl Acad Sci, 95(5):2 105-2 110.

    LUSCOMBE NM, GREENBAUM D, GERSTEIN M, et al, 2001. What is bioinformatics? A proposed definition and overview of the field [J]. Meth Inf Med, 40(4):346-358.

    LUISA ML, MANUEL RC, FRANCESCA G,et al, 2000. Carotenoid biosynthesis during tomato fruit development: regulatory role of 1-deoxy-D-sylulose 5-phosphate synthase [J]. Plant J, 22(6):503-513.

    PENG HC, 2008. Bioimage informatics: a new area of engineering biology [J]. Bioinformatics, 24(17):1 827-1 836.

    PATRICIA BP, MANUEL RC, 2006. Carotenoid biotechnology in plants for nutritionally improved foods [J]. Physiol Plant, 126(3):369-381.

    RAO AV, RAO LG, 2007. Carotenoids and human health [J]. Pharmacol Res, 55:207-216.

    SACCHETTINI JC, POULTER CD. 1997. Creating isoprenoid diversity [J]. Science, 277:1 788-1 789.

    SANAE K, AKEMI O, 2006. Regulation of carotenoid biosynthesis in petals and leaves of chrysanthemum (Chrysanthemummorifolium) [J]. Physiol Plant, 128(3):436-447.SAWITRI K, WALLIE S, 2005. Molecular cloning and expression of a cDNA encoding 1-deoxy-D-xylulose-5-phosphate synthase from oil palmElaeisguineensisJacq. [J]. Plant Sci,169:571-578.

    SHAN L, LI L, 2008. Carotenoid metabolism: biosynthesis, regulation, and beyond [J]. J Integr Plant Biol, 50(7):778-785.

    TIAN HS, CHENG QL, YUAN YH, et al, 2010. Coordinated regulation of gene expression for carotenoid metabolism inChlamydomonasreinhardtii[J]. J Integr Plant Biol, 52(10):868-878.

    VINCENT C, DIRK JL, 2011. From glycosylation disorders to dolichol biosynthesis defects: a new class of metabolic diseases [J]. J Inherit Metab Dis, 34(4):859-867.

    VRANOVA E, COMAN D, GRUISSEM W, 2013. Network analysis of the MVA and MEP pathway for isoprenoid synthesis [J]. Ann Rev Plant Biol,64:665-700.

    WANKE M, SKORUPINSKA TK, SWIEZEWSKA E, 2001. Isoprenoid biosynthesis via 1-deoxy-D-xylulose-5-phosphate/2-C-methyl-D-erythritol-4-phosphate (DOXP/MEP) pathway [J]. Acta Biochim Pol, 48(3):663-672.

    ZHANG M, LI K, ZHANG C, et al, 2009. Identification and characterization of class 1 DXS gene encoding 1-deoxy-D-xylulose-5-phosphate synthase, the first committed enzyme of the MEP pathway from soybean [J]. Mol Biol Rep, 36(5):879-887.

    Bioinformatics analysis ofDXSgene from six tropical plants including oil palm (ElaeisguineensisJacq.)

    SHI Peng, CAO Hong-Xing, LI Dong-Xia, WANG Yong, LEI Xin-Tao*

    (HainanKeyBiologicalLaboratoryofTropicalOilCrops,CoconutResearchInstitute,ChineseAcademyofTropicalAgriculturalSciences, Wenchang 571339, China )

    Abstract:The paper unraveled structure characteristics, function features and phylogenetic relationships among species through analyzing nucleic acid and protein sequences of key gene DXS controlling isoprenoids synthesis of tropical plants as oil palm (Elaeis guineensis Jacq.) by using bioinformatics tools. Isoprenoids like carotenoids and sterols play an important role in life of plant, and have important significance for the protection of human health. Tropical plants as oil palm contains abundant isoprenoids like carotenes and vitamin E, MEP is one of the important pathway for isoprenoids biosynthesis. DXS as the first rate limiting enzyme in MEP pathway, its functions were highly conserved in diverse tropical plants. In order to understand the structure and function characteristics of DXS from tropical plants as oil palm, DXS genes nucleic acid and protein sequences from different plants including oil palm were analyzed such as physicochemical property, protein structure, and conserved domain by the tools of bioinformatics in this study. The results showed that the start codon of DXS gene was ATG, but termination codon was divided into TAG, TAA and TGA; DXS protein belonged to a hydrophilic protein without a signal peptide; It might played a role in the chloroplast stroma as a transport protein; No transmembrane domain was found, there were 36 phosphorylation sites, including 17 serine, 11 threonine and 8 tyrosine sites respectively; Random coil and alpha helix were the major structural elements of secondary structure; It had DXS enzyme character, thiamin diphosphate (TPP) binding sites and PLN02582 conserved multi-domains with tertiary structure analysis; Domains of DXS protein from different plant, were very conservative, and could be used as an important basis for determining genetic relationships among species. The results would provide further information for structure and functional analysis of DXS, and lay the foundation for molecular mechanism research and genetic improvement of quality traits for tropical plants including oil palm.

    Key words:oil palm (Elaeis guineensis Jacq.), DXS gene, bioinformatics, phylogenetic tree

    中圖分類號:Q811.4

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

    文章編號:1000-3142(2016)04-0471-08

    作者簡介:石鵬(1987-),男,湖南安鄉(xiāng)人,碩士,研究實(shí)習(xí)員,主要從事油棕種質(zhì)資源評價(jià)等研究,(E-mail)624811593@qq.com。*通訊作者: 雷新濤,博士,研究員,主要從事熱帶油料作物遺傳育種等研究,(E-mail)xtlei@263.com。

    基金項(xiàng)目:國家林業(yè)局“948”項(xiàng)目(No.2013-4-40)[Supported by the Plan of Introducing the Advanced Technology for Agriculture of State Forestry Administration(No.2013-4-40)]。

    *收稿日期:2014-12-19修回日期: 2015-02-15

    DOI:10.11931/guihaia.gxzw201406013

    石鵬,曹紅星,李東霞,等. 油棕等熱帶植物DXS基因的生物信息學(xué)分析 [J]. 廣西植物, 2016, 36(4):471-478

    SHI P,CAO HX,LI DX,et al. Bioinformatics analysis ofDXSgene from six tropical plants including oil palm (ElaeisguineensisJacq.)[J]. Guihaia, 2016, 36(4):471-478

    猜你喜歡
    油棕生物信息學(xué)
    棕櫚栽培管理初探
    淺談醫(yī)學(xué)院校生物信息學(xué)專業(yè)青年教師規(guī)范培訓(xùn)模式的建立
    “PBL+E—learning”教學(xué)模式探索
    移動教學(xué)在生物信息學(xué)課程改革中的應(yīng)用
    今傳媒(2016年11期)2016-12-19 11:35:50
    中醫(yī)大數(shù)據(jù)下生物信息學(xué)的發(fā)展及教育模式淺析
    數(shù)據(jù)挖掘技術(shù)在生物信息學(xué)中的應(yīng)用
    生物信息學(xué)課堂危機(jī)及對策研究
    科技視界(2016年23期)2016-11-04 10:07:53
    泰國油棕果市價(jià)連跌2個月
    油棕
    今年生產(chǎn)40萬t棕油指標(biāo)
    亚洲成人中文字幕在线播放| 欧美精品人与动牲交sv欧美| 久久久久九九精品影院| 久久久亚洲精品成人影院| 涩涩av久久男人的天堂| 女人十人毛片免费观看3o分钟| 日韩av在线免费看完整版不卡| eeuss影院久久| 精品人妻一区二区三区麻豆| av卡一久久| 免费观看a级毛片全部| 高清视频免费观看一区二区| 国产精品三级大全| 国产在线男女| 国产成人aa在线观看| 国内少妇人妻偷人精品xxx网站| 成人美女网站在线观看视频| 亚洲欧美清纯卡通| 久久久久久久久大av| 国产老妇女一区| 成人无遮挡网站| 草草在线视频免费看| 国产成人freesex在线| 国产v大片淫在线免费观看| 国产综合精华液| 久久久亚洲精品成人影院| 精品视频人人做人人爽| 美女内射精品一级片tv| 国产伦精品一区二区三区视频9| 国产精品.久久久| 国产高清三级在线| 国产精品精品国产色婷婷| 日韩欧美 国产精品| 国产欧美日韩精品一区二区| 日韩一区二区三区影片| 搡老乐熟女国产| 在线看a的网站| 国产成人精品一,二区| 在线 av 中文字幕| 听说在线观看完整版免费高清| 免费看a级黄色片| 国产探花极品一区二区| 欧美日韩亚洲高清精品| av卡一久久| 丰满人妻一区二区三区视频av| 欧美精品一区二区大全| kizo精华| av在线app专区| 最近2019中文字幕mv第一页| 亚洲av一区综合| 久久久久久伊人网av| 国产免费福利视频在线观看| 一级毛片aaaaaa免费看小| 久久久久国产精品人妻一区二区| 在线观看国产h片| 国产伦精品一区二区三区视频9| 国产免费一级a男人的天堂| 最近最新中文字幕大全电影3| 亚洲精品日本国产第一区| 男女边吃奶边做爰视频| 一区二区三区精品91| 久久久久国产网址| 又爽又黄a免费视频| 青青草视频在线视频观看| 精华霜和精华液先用哪个| 久久久久久久久久人人人人人人| 免费少妇av软件| 欧美3d第一页| 欧美激情在线99| 亚洲av成人精品一二三区| 日韩一区二区视频免费看| 在线 av 中文字幕| 夫妻午夜视频| 视频区图区小说| 天天躁夜夜躁狠狠久久av| 日本熟妇午夜| 国产日韩欧美亚洲二区| 麻豆成人午夜福利视频| 国产精品一及| 国产免费福利视频在线观看| 人人妻人人爽人人添夜夜欢视频 | 国产精品福利在线免费观看| 白带黄色成豆腐渣| 99热这里只有是精品50| 视频区图区小说| 亚洲av二区三区四区| 高清毛片免费看| 亚洲精品影视一区二区三区av| 各种免费的搞黄视频| xxx大片免费视频| freevideosex欧美| 伦理电影大哥的女人| 国产欧美亚洲国产| 国产成人freesex在线| 国产成人freesex在线| 国产大屁股一区二区在线视频| 天堂中文最新版在线下载 | av国产久精品久网站免费入址| 69av精品久久久久久| 国产精品女同一区二区软件| 99久国产av精品国产电影| 麻豆精品久久久久久蜜桃| 亚洲久久久久久中文字幕| 男女无遮挡免费网站观看| kizo精华| 亚洲电影在线观看av| 日韩人妻高清精品专区| 免费电影在线观看免费观看| 一二三四中文在线观看免费高清| 国产乱来视频区| 高清在线视频一区二区三区| 小蜜桃在线观看免费完整版高清| 免费观看av网站的网址| 搡女人真爽免费视频火全软件| 亚洲国产欧美人成| 日韩亚洲欧美综合| 国产爽快片一区二区三区| 真实男女啪啪啪动态图| 国产一区二区三区综合在线观看 | 国产视频内射| 97精品久久久久久久久久精品| 成人黄色视频免费在线看| 免费av毛片视频| 色哟哟·www| videossex国产| 色视频www国产| 99热这里只有是精品50| 最近中文字幕2019免费版| 少妇人妻一区二区三区视频| 大片电影免费在线观看免费| 黄色怎么调成土黄色| 亚洲色图综合在线观看| 高清欧美精品videossex| a级毛色黄片| av黄色大香蕉| 国产亚洲av片在线观看秒播厂| 51国产日韩欧美| 亚洲精品国产色婷婷电影| 亚洲精品乱码久久久久久按摩| 国产人妻一区二区三区在| 嫩草影院入口| a级毛片免费高清观看在线播放| 欧美高清成人免费视频www| 亚洲婷婷狠狠爱综合网| 深爱激情五月婷婷| 久久久久久伊人网av| 精品久久久久久久久亚洲| 22中文网久久字幕| 国产一区二区三区综合在线观看 | 国产欧美另类精品又又久久亚洲欧美| 国产av国产精品国产| 国内精品宾馆在线| 新久久久久国产一级毛片| 国产精品99久久久久久久久| 欧美亚洲 丝袜 人妻 在线| 日韩av在线免费看完整版不卡| 一本久久精品| 大香蕉久久网| av专区在线播放| 国产成人精品一,二区| 成年免费大片在线观看| 青春草国产在线视频| 亚洲精品456在线播放app| 中文字幕制服av| 亚洲av成人精品一区久久| 日日摸夜夜添夜夜添av毛片| 欧美人与善性xxx| 精品熟女少妇av免费看| 亚洲欧美精品专区久久| 久久久久精品久久久久真实原创| 国产一区亚洲一区在线观看| 久久99热这里只频精品6学生| 亚洲欧美日韩卡通动漫| 国产视频首页在线观看| 亚州av有码| 一边亲一边摸免费视频| 久久人人爽人人片av| 伦精品一区二区三区| 少妇丰满av| 欧美成人一区二区免费高清观看| 午夜免费男女啪啪视频观看| 亚洲欧洲国产日韩| 三级国产精品片| 精品酒店卫生间| 黑人高潮一二区| 在线免费观看不下载黄p国产| 免费大片黄手机在线观看| 国产精品久久久久久精品电影| 精品一区二区三区视频在线| 黄色配什么色好看| 啦啦啦在线观看免费高清www| 亚洲精品成人av观看孕妇| 特级一级黄色大片| 国产男人的电影天堂91| 韩国高清视频一区二区三区| 免费观看在线日韩| 综合色av麻豆| 涩涩av久久男人的天堂| av播播在线观看一区| 欧美zozozo另类| 国产乱来视频区| 国产精品蜜桃在线观看| 国产成人aa在线观看| 校园人妻丝袜中文字幕| 一级二级三级毛片免费看| av国产久精品久网站免费入址| 高清毛片免费看| 边亲边吃奶的免费视频| 一区二区三区乱码不卡18| 亚洲人成网站在线观看播放| 最近最新中文字幕免费大全7| 国产成人aa在线观看| 亚洲色图av天堂| 欧美国产精品一级二级三级 | 国产精品爽爽va在线观看网站| 欧美极品一区二区三区四区| 中文字幕av成人在线电影| 校园人妻丝袜中文字幕| 人妻 亚洲 视频| 日韩伦理黄色片| 国产亚洲精品久久久com| 99re6热这里在线精品视频| 丰满人妻一区二区三区视频av| 免费在线观看成人毛片| 一个人看的www免费观看视频| 少妇的逼好多水| 国产精品久久久久久精品古装| 免费观看在线日韩| 大片免费播放器 马上看| kizo精华| 亚洲欧洲日产国产| 成人美女网站在线观看视频| 日韩中字成人| 人妻少妇偷人精品九色| 一级毛片黄色毛片免费观看视频| 18禁裸乳无遮挡免费网站照片| 亚洲在线观看片| 国产精品国产三级专区第一集| 天天躁日日操中文字幕| 欧美人与善性xxx| 大片电影免费在线观看免费| 乱码一卡2卡4卡精品| 国内少妇人妻偷人精品xxx网站| 午夜福利视频精品| 亚洲欧美日韩东京热| 中文精品一卡2卡3卡4更新| 91精品国产九色| 成人毛片60女人毛片免费| 午夜福利视频精品| 国产探花极品一区二区| 亚洲精品国产成人久久av| 在线观看人妻少妇| 日本wwww免费看| 秋霞在线观看毛片| 欧美精品人与动牲交sv欧美| 一级毛片我不卡| 色5月婷婷丁香| 国产成年人精品一区二区| 老司机影院成人| 国产片特级美女逼逼视频| 亚洲婷婷狠狠爱综合网| 免费电影在线观看免费观看| 国产精品久久久久久av不卡| 乱码一卡2卡4卡精品| 全区人妻精品视频| 国产一区二区亚洲精品在线观看| 一本久久精品| 中文字幕制服av| 九草在线视频观看| 男女边吃奶边做爰视频| 老师上课跳d突然被开到最大视频| 一区二区三区精品91| 热99国产精品久久久久久7| 久久久久久久久久成人| 99九九线精品视频在线观看视频| 人人妻人人看人人澡| 最近最新中文字幕免费大全7| 亚洲国产成人一精品久久久| av播播在线观看一区| 国产片特级美女逼逼视频| 国产男女超爽视频在线观看| 欧美zozozo另类| 亚洲无线观看免费| 听说在线观看完整版免费高清| 真实男女啪啪啪动态图| 久久99热这里只有精品18| 国产片特级美女逼逼视频| 色网站视频免费| 欧美性猛交╳xxx乱大交人| 观看免费一级毛片| 菩萨蛮人人尽说江南好唐韦庄| 国产亚洲一区二区精品| 欧美精品一区二区大全| 男人和女人高潮做爰伦理| .国产精品久久| 久久久久久久精品精品| 日韩在线高清观看一区二区三区| 国产亚洲午夜精品一区二区久久 | 日韩中字成人| 特大巨黑吊av在线直播| 日本午夜av视频| 狂野欧美白嫩少妇大欣赏| 日本免费在线观看一区| 国产成人a∨麻豆精品| 亚洲在久久综合| 激情 狠狠 欧美| 午夜精品一区二区三区免费看| 成人国产麻豆网| 干丝袜人妻中文字幕| 久久久亚洲精品成人影院| 免费观看的影片在线观看| 色综合色国产| 少妇裸体淫交视频免费看高清| 国产午夜精品一二区理论片| 色哟哟·www| 免费少妇av软件| 最近最新中文字幕免费大全7| 亚洲不卡免费看| 国产精品无大码| 51国产日韩欧美| 全区人妻精品视频| 精品久久久精品久久久| 国产白丝娇喘喷水9色精品| 亚洲av一区综合| 国产欧美日韩精品一区二区| 中国三级夫妇交换| 97在线人人人人妻| 亚洲av.av天堂| 亚洲久久久久久中文字幕| 亚洲欧美日韩无卡精品| 亚洲欧美日韩另类电影网站 | 大话2 男鬼变身卡| 亚洲精品日本国产第一区| 欧美三级亚洲精品| 99久久精品国产国产毛片| 成年版毛片免费区| 日本av手机在线免费观看| 小蜜桃在线观看免费完整版高清| 欧美国产精品一级二级三级 | 久久精品国产亚洲网站| 免费看a级黄色片| 国产黄色视频一区二区在线观看| 欧美精品人与动牲交sv欧美| 人妻少妇偷人精品九色| 午夜日本视频在线| 色视频www国产| 日产精品乱码卡一卡2卡三| 人妻少妇偷人精品九色| 男人舔奶头视频| 午夜爱爱视频在线播放| 亚洲一区二区三区欧美精品 | 国内揄拍国产精品人妻在线| 亚洲av不卡在线观看| 国产精品成人在线| 欧美 日韩 精品 国产| 国产午夜精品久久久久久一区二区三区| 久久久久国产精品人妻一区二区| 亚洲一区二区三区欧美精品 | 寂寞人妻少妇视频99o| 一区二区av电影网| 欧美三级亚洲精品| 777米奇影视久久| 免费黄色在线免费观看| 99精国产麻豆久久婷婷| 高清欧美精品videossex| 色网站视频免费| 能在线免费看毛片的网站| 国产白丝娇喘喷水9色精品| 国产免费视频播放在线视频| 国产毛片a区久久久久| 亚洲,一卡二卡三卡| 日韩国内少妇激情av| 免费大片18禁| 免费看光身美女| 嫩草影院新地址| 欧美激情国产日韩精品一区| 日日摸夜夜添夜夜爱| 国产欧美亚洲国产| 好男人视频免费观看在线| 人人妻人人澡人人爽人人夜夜| 卡戴珊不雅视频在线播放| 国产成人精品福利久久| 最后的刺客免费高清国语| 国产高清不卡午夜福利| 美女cb高潮喷水在线观看| 精品熟女少妇av免费看| 丝袜脚勾引网站| 中文字幕av成人在线电影| 一级黄片播放器| 日日啪夜夜撸| 婷婷色av中文字幕| 男女下面进入的视频免费午夜| 99热这里只有精品一区| 少妇人妻久久综合中文| 真实男女啪啪啪动态图| 日本猛色少妇xxxxx猛交久久| 国产淫语在线视频| 综合色av麻豆| 不卡视频在线观看欧美| 色播亚洲综合网| .国产精品久久| 亚洲欧美清纯卡通| 日韩人妻高清精品专区| 青春草国产在线视频| 国产 一区精品| 欧美成人精品欧美一级黄| 最后的刺客免费高清国语| 日本免费在线观看一区| 水蜜桃什么品种好| 亚洲成人精品中文字幕电影| 中文字幕亚洲精品专区| 国产一区二区在线观看日韩| 新久久久久国产一级毛片| www.色视频.com| 久久久久久久亚洲中文字幕| 黑人高潮一二区| 国产日韩欧美在线精品| 国产国拍精品亚洲av在线观看| 日韩欧美一区视频在线观看 | 美女主播在线视频| 国产久久久一区二区三区| 一级毛片 在线播放| 黄色欧美视频在线观看| 激情 狠狠 欧美| 国产伦理片在线播放av一区| 亚洲av电影在线观看一区二区三区 | 国产欧美亚洲国产| 一级毛片我不卡| 1000部很黄的大片| 一级黄片播放器| 久久国产乱子免费精品| av线在线观看网站| 啦啦啦在线观看免费高清www| 少妇的逼水好多| 成人漫画全彩无遮挡| 日韩,欧美,国产一区二区三区| 亚洲国产精品成人综合色| 波野结衣二区三区在线| 在线免费十八禁| 自拍偷自拍亚洲精品老妇| 蜜桃亚洲精品一区二区三区| av免费在线看不卡| 秋霞伦理黄片| 青青草视频在线视频观看| 久久久久久久久久久免费av| 高清欧美精品videossex| 久久热精品热| 欧美日韩亚洲高清精品| 亚州av有码| 精品久久久久久久人妻蜜臀av| 青春草亚洲视频在线观看| 不卡视频在线观看欧美| 日韩免费高清中文字幕av| 久久久久久九九精品二区国产| 久久精品国产亚洲网站| 亚洲美女视频黄频| 精品久久久噜噜| 国产精品爽爽va在线观看网站| 天天一区二区日本电影三级| 婷婷色综合www| 亚洲精品乱码久久久久久按摩| 久热久热在线精品观看| 国产精品av视频在线免费观看| 老师上课跳d突然被开到最大视频| 精品久久久久久久久av| 伊人久久精品亚洲午夜| 一级毛片我不卡| 精品人妻偷拍中文字幕| 国产精品久久久久久精品电影| 亚洲精品视频女| 中文乱码字字幕精品一区二区三区| 成人亚洲欧美一区二区av| 一个人观看的视频www高清免费观看| 日韩国内少妇激情av| 国产黄片视频在线免费观看| 久久久欧美国产精品| 色哟哟·www| 免费人成在线观看视频色| 精品午夜福利在线看| 在线观看人妻少妇| 日日撸夜夜添| 伊人久久国产一区二区| 菩萨蛮人人尽说江南好唐韦庄| 精品99又大又爽又粗少妇毛片| 99热国产这里只有精品6| 欧美高清性xxxxhd video| 日韩中字成人| 蜜桃久久精品国产亚洲av| 国产欧美另类精品又又久久亚洲欧美| 久久亚洲国产成人精品v| 国产av国产精品国产| 天堂中文最新版在线下载 | 免费看光身美女| 最近手机中文字幕大全| 午夜福利在线在线| 伊人久久国产一区二区| 综合色av麻豆| 老司机影院毛片| 51国产日韩欧美| 免费播放大片免费观看视频在线观看| 亚洲在线观看片| 超碰av人人做人人爽久久| 九九在线视频观看精品| 国产精品国产三级国产专区5o| 少妇丰满av| 国产淫语在线视频| 女人十人毛片免费观看3o分钟| 日本黄大片高清| 麻豆精品久久久久久蜜桃| 少妇人妻久久综合中文| 国产成人午夜福利电影在线观看| 卡戴珊不雅视频在线播放| 交换朋友夫妻互换小说| 久久久久久久国产电影| 国产黄a三级三级三级人| 久久鲁丝午夜福利片| 禁无遮挡网站| 日本与韩国留学比较| 涩涩av久久男人的天堂| 美女xxoo啪啪120秒动态图| 国产老妇女一区| 人妻系列 视频| 日本wwww免费看| 午夜视频国产福利| 18禁在线无遮挡免费观看视频| 91久久精品电影网| 欧美日韩亚洲高清精品| 国产成人福利小说| 成人免费观看视频高清| 国产精品一区www在线观看| 精华霜和精华液先用哪个| 夜夜爽夜夜爽视频| 免费观看无遮挡的男女| 久久久久久九九精品二区国产| 久久久久久久国产电影| 欧美另类一区| 亚洲经典国产精华液单| 午夜激情福利司机影院| 久久久久久伊人网av| 夜夜爽夜夜爽视频| 国产亚洲午夜精品一区二区久久 | 日本熟妇午夜| 日韩强制内射视频| 欧美高清成人免费视频www| 国产精品爽爽va在线观看网站| 国产午夜精品一二区理论片| 18禁在线播放成人免费| 99re6热这里在线精品视频| 午夜免费鲁丝| 看免费成人av毛片| 亚洲国产精品成人综合色| 三级经典国产精品| 亚洲国产最新在线播放| 永久网站在线| 黄色视频在线播放观看不卡| 91久久精品电影网| 免费在线观看成人毛片| av国产免费在线观看| av黄色大香蕉| 少妇的逼好多水| 人人妻人人爽人人添夜夜欢视频 | 国产亚洲av片在线观看秒播厂| 国产黄a三级三级三级人| 人妻 亚洲 视频| 男人爽女人下面视频在线观看| 2022亚洲国产成人精品| 国产一级毛片在线| 日本wwww免费看| 成人综合一区亚洲| 97在线视频观看| 十八禁网站网址无遮挡 | 白带黄色成豆腐渣| 亚洲高清免费不卡视频| 色网站视频免费| 一本一本综合久久| 男女边摸边吃奶| 国产精品一二三区在线看| 国产一区二区在线观看日韩| 老女人水多毛片| 高清av免费在线| 日日摸夜夜添夜夜添av毛片| 国产高清有码在线观看视频| 肉色欧美久久久久久久蜜桃 | 天天躁夜夜躁狠狠久久av| 美女高潮的动态| 久久久久精品久久久久真实原创| av在线老鸭窝| 男女国产视频网站| 国产伦在线观看视频一区| 91精品伊人久久大香线蕉| 自拍偷自拍亚洲精品老妇| 久久久久久久大尺度免费视频| 久久久久久久久大av| 精品久久久噜噜| 国产高清国产精品国产三级 | 日本-黄色视频高清免费观看| 国产精品一区www在线观看| 国产精品成人在线| 欧美潮喷喷水| 一级毛片我不卡| 婷婷色综合大香蕉| 国产伦精品一区二区三区四那| 精品久久久精品久久久| 日本午夜av视频| 国产午夜精品一二区理论片| 99久久精品国产国产毛片| 又大又黄又爽视频免费| 国产精品人妻久久久久久| 国产精品嫩草影院av在线观看| 一区二区av电影网| 久久久久久久久久成人| 欧美xxxx性猛交bbbb| 日韩大片免费观看网站| 伦理电影大哥的女人| 国产精品国产三级国产专区5o|