• <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)
    午夜精品一区二区三区免费看| 韩国高清视频一区二区三区| 97热精品久久久久久| 少妇人妻精品综合一区二区| 亚洲成人一二三区av| av线在线观看网站| 中文字幕制服av| 日韩强制内射视频| 国产伦精品一区二区三区视频9| 国产三级在线视频| 国产精品国产三级专区第一集| 插逼视频在线观看| 一级毛片黄色毛片免费观看视频| 亚洲欧美日韩东京热| 国产精品99久久久久久久久| h日本视频在线播放| 精品一区二区三卡| 天天一区二区日本电影三级| 亚洲精品,欧美精品| 一边亲一边摸免费视频| 日本一二三区视频观看| 男女那种视频在线观看| 黄色日韩在线| 搡女人真爽免费视频火全软件| 校园人妻丝袜中文字幕| 国产精品美女特级片免费视频播放器| 老师上课跳d突然被开到最大视频| 99re6热这里在线精品视频| 日韩av不卡免费在线播放| 黑人高潮一二区| 欧美成人一区二区免费高清观看| 九九爱精品视频在线观看| 久久久久久久久久成人| 2021天堂中文幕一二区在线观| 亚洲欧洲日产国产| 亚洲欧美中文字幕日韩二区| 秋霞伦理黄片| 一区二区三区高清视频在线| 天堂√8在线中文| 一区二区三区乱码不卡18| 精品久久久久久成人av| 久久鲁丝午夜福利片| 免费看a级黄色片| 久久久久久久午夜电影| 少妇人妻精品综合一区二区| 日韩成人伦理影院| 中文欧美无线码| 国产高清国产精品国产三级 | 人人妻人人澡欧美一区二区| 国产成人精品婷婷| 美女cb高潮喷水在线观看| 天堂影院成人在线观看| 精品国产一区二区三区久久久樱花 | 亚洲av免费在线观看| 97热精品久久久久久| 亚洲av在线观看美女高潮| 日韩,欧美,国产一区二区三区| 黄片无遮挡物在线观看| 国产黄色小视频在线观看| 99热这里只有是精品在线观看| 天堂av国产一区二区熟女人妻| 午夜福利网站1000一区二区三区| xxx大片免费视频| 午夜精品在线福利| 精品久久久噜噜| 全区人妻精品视频| av又黄又爽大尺度在线免费看| 久久久久九九精品影院| 夜夜爽夜夜爽视频| 日韩欧美 国产精品| 亚洲18禁久久av| 在线 av 中文字幕| 国产探花在线观看一区二区| 久久久久久九九精品二区国产| 国产av不卡久久| 日韩亚洲欧美综合| 黄色欧美视频在线观看| 欧美精品一区二区大全| 久久久精品94久久精品| 亚洲人与动物交配视频| 波多野结衣巨乳人妻| 在线观看av片永久免费下载| 91久久精品国产一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 91在线精品国自产拍蜜月| 久热久热在线精品观看| 国产精品一区二区三区四区免费观看| 久久这里只有精品中国| 六月丁香七月| 乱码一卡2卡4卡精品| 美女脱内裤让男人舔精品视频| av播播在线观看一区| 免费大片18禁| 欧美 日韩 精品 国产| 国产成人freesex在线| 亚洲成人中文字幕在线播放| 久久久精品欧美日韩精品| 美女cb高潮喷水在线观看| 午夜免费激情av| 亚洲av一区综合| 欧美成人a在线观看| 99热6这里只有精品| 亚洲不卡免费看| 黄色一级大片看看| 久久久午夜欧美精品| 亚洲精品乱久久久久久| 一区二区三区四区激情视频| 97精品久久久久久久久久精品| 老女人水多毛片| 免费av不卡在线播放| 免费播放大片免费观看视频在线观看| 高清av免费在线| 国产老妇伦熟女老妇高清| 国产视频内射| 免费观看精品视频网站| 大陆偷拍与自拍| 女人久久www免费人成看片| 国产一区二区在线观看日韩| 偷拍熟女少妇极品色| 国产淫片久久久久久久久| 国产精品一二三区在线看| 欧美成人午夜免费资源| 成人一区二区视频在线观看| 一级爰片在线观看| 精品99又大又爽又粗少妇毛片| 精品久久久噜噜| 日韩av不卡免费在线播放| 一级毛片久久久久久久久女| 蜜桃久久精品国产亚洲av| 美女脱内裤让男人舔精品视频| 中文精品一卡2卡3卡4更新| 最近的中文字幕免费完整| 高清视频免费观看一区二区 | 久久久久久久国产电影| 亚洲成人av在线免费| 亚州av有码| 乱系列少妇在线播放| kizo精华| 久久精品国产自在天天线| 亚洲高清免费不卡视频| av在线老鸭窝| 久久精品人妻少妇| 日日摸夜夜添夜夜爱| 久久精品国产亚洲av天美| 精品国产露脸久久av麻豆 | 日韩欧美国产在线观看| 久久久亚洲精品成人影院| 亚洲av电影在线观看一区二区三区 | 国产一级毛片在线| 亚洲成人av在线免费| 少妇的逼好多水| 成年免费大片在线观看| h日本视频在线播放| 亚洲一级一片aⅴ在线观看| 水蜜桃什么品种好| 岛国毛片在线播放| 久久久久久久久久黄片| 亚洲精品成人av观看孕妇| 精品少妇黑人巨大在线播放| 七月丁香在线播放| 99热全是精品| 国产成人精品福利久久| av又黄又爽大尺度在线免费看| 99热全是精品| 夫妻性生交免费视频一级片| 丰满乱子伦码专区| 久久久久久久大尺度免费视频| 精品酒店卫生间| 亚洲国产欧美在线一区| 国产一级毛片在线| 国产精品一区二区性色av| 日韩强制内射视频| 精品久久久久久久久久久久久| 一区二区三区高清视频在线| 久久97久久精品| 久久久a久久爽久久v久久| 亚洲av国产av综合av卡| 亚洲婷婷狠狠爱综合网| 看黄色毛片网站| 精品久久久久久久久亚洲| 精品久久久精品久久久| 18+在线观看网站| 男人舔女人下体高潮全视频| 乱系列少妇在线播放| 亚洲欧美日韩东京热| 亚洲av免费在线观看| 麻豆精品久久久久久蜜桃| 色综合色国产| 床上黄色一级片| 嘟嘟电影网在线观看| 九九在线视频观看精品| 丰满乱子伦码专区| 天堂影院成人在线观看| 亚洲精品,欧美精品| 欧美成人一区二区免费高清观看| 成人一区二区视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 久久国内精品自在自线图片| 伊人久久精品亚洲午夜| 国产精品一区二区三区四区久久| 中文字幕制服av| 1000部很黄的大片| 69人妻影院| 最近最新中文字幕大全电影3| 汤姆久久久久久久影院中文字幕 | 国产伦在线观看视频一区| 国产成人免费观看mmmm| 建设人人有责人人尽责人人享有的 | 亚洲av免费高清在线观看| 久久亚洲国产成人精品v| 精品久久久久久电影网| 天堂俺去俺来也www色官网 | 免费观看无遮挡的男女| 赤兔流量卡办理| 九九久久精品国产亚洲av麻豆| 国产精品一区www在线观看| 日韩,欧美,国产一区二区三区| 777米奇影视久久| 少妇高潮的动态图| 国产精品人妻久久久影院| 国产中年淑女户外野战色| 夫妻性生交免费视频一级片| 一个人免费在线观看电影| 午夜精品一区二区三区免费看| 97人妻精品一区二区三区麻豆| 99久久精品国产国产毛片| 中国美白少妇内射xxxbb| 国产精品国产三级国产av玫瑰| 秋霞在线观看毛片| 全区人妻精品视频| 啦啦啦中文免费视频观看日本| 国产成人a区在线观看| 欧美日韩在线观看h| 午夜精品在线福利| 国产精品1区2区在线观看.| 国产精品日韩av在线免费观看| 免费黄网站久久成人精品| 亚洲av中文av极速乱| 日韩欧美精品免费久久| 国产探花极品一区二区| 精品久久国产蜜桃| 精品人妻一区二区三区麻豆| 亚洲欧美日韩无卡精品| 狂野欧美白嫩少妇大欣赏| 白带黄色成豆腐渣| 美女xxoo啪啪120秒动态图| 国产综合懂色| 免费播放大片免费观看视频在线观看| 国产视频内射| 久久久久久九九精品二区国产| 亚洲精品影视一区二区三区av| 久久久精品免费免费高清| 亚洲av.av天堂| 成年免费大片在线观看| 久久热精品热| 男女视频在线观看网站免费| 国产视频内射| 免费大片18禁| 男女边摸边吃奶| 亚洲av成人精品一二三区| 观看美女的网站| 能在线免费观看的黄片| 真实男女啪啪啪动态图| 精品酒店卫生间| 免费播放大片免费观看视频在线观看| 2021少妇久久久久久久久久久| av在线亚洲专区| 97人妻精品一区二区三区麻豆| 国产综合精华液| 亚洲人与动物交配视频| 成人午夜精彩视频在线观看| 寂寞人妻少妇视频99o| 中文欧美无线码| 日韩强制内射视频| 在线免费十八禁| 国产在视频线在精品| 成年版毛片免费区| 国产不卡一卡二| 精品人妻熟女av久视频| 亚洲电影在线观看av| 91久久精品国产一区二区成人| 看黄色毛片网站| 五月玫瑰六月丁香| 直男gayav资源| 美女大奶头视频| 成年人午夜在线观看视频 | 菩萨蛮人人尽说江南好唐韦庄| 欧美丝袜亚洲另类| 白带黄色成豆腐渣| 日韩成人av中文字幕在线观看| 亚洲在线观看片| 美女大奶头视频| 一级黄片播放器| 免费观看在线日韩| 久久久久九九精品影院| 欧美最新免费一区二区三区| 91久久精品国产一区二区成人| 国产伦一二天堂av在线观看| 久久久精品欧美日韩精品| 国产探花在线观看一区二区| 五月玫瑰六月丁香| 亚洲人成网站高清观看| 色综合亚洲欧美另类图片| 久久草成人影院| 亚洲精品乱码久久久久久按摩| 午夜免费激情av| 成人亚洲精品一区在线观看 | 美女大奶头视频| 熟女人妻精品中文字幕| 91狼人影院| 日产精品乱码卡一卡2卡三| 国产黄色视频一区二区在线观看| 免费黄色在线免费观看| av线在线观看网站| 免费大片黄手机在线观看| 国内少妇人妻偷人精品xxx网站| 国产成年人精品一区二区| 狂野欧美白嫩少妇大欣赏| 国产伦理片在线播放av一区| av免费观看日本| 免费看美女性在线毛片视频| 秋霞在线观看毛片| 99re6热这里在线精品视频| 久久99蜜桃精品久久| 日本爱情动作片www.在线观看| 一级av片app| 国产精品人妻久久久久久| 自拍偷自拍亚洲精品老妇| 国产av国产精品国产| 色视频www国产| 久久久久久久久久成人| freevideosex欧美| 在线观看美女被高潮喷水网站| 欧美丝袜亚洲另类| 国产探花极品一区二区| 国产 一区精品| 五月天丁香电影| 欧美人与善性xxx| 国产视频首页在线观看| 赤兔流量卡办理| 男女边吃奶边做爰视频| 中文字幕av在线有码专区| 一级片'在线观看视频| 国产大屁股一区二区在线视频| 男女下面进入的视频免费午夜| 欧美3d第一页| 国产三级在线视频| 国产精品伦人一区二区| 亚洲成人中文字幕在线播放| 亚洲精品国产av蜜桃| 亚洲欧美日韩无卡精品| 久久久久久久大尺度免费视频| 亚洲熟女精品中文字幕| 91精品国产九色| 亚洲熟女精品中文字幕| 嘟嘟电影网在线观看| 亚洲成人精品中文字幕电影| 欧美成人精品欧美一级黄| 91av网一区二区| 久久精品久久久久久久性| 亚洲av中文av极速乱| 午夜激情欧美在线| 日本与韩国留学比较| 久久人人爽人人片av| 午夜福利在线在线| 婷婷色av中文字幕| 又大又黄又爽视频免费| 亚洲精品一区蜜桃| 久久久久久伊人网av| 天堂影院成人在线观看| 国产精品人妻久久久影院| 国产成人a区在线观看| 国产一区亚洲一区在线观看| 日韩电影二区| 国产在视频线精品| 亚洲欧美精品自产自拍| 精品一区二区三区视频在线| 91午夜精品亚洲一区二区三区| 成年人午夜在线观看视频 | 淫秽高清视频在线观看| 国产伦理片在线播放av一区| 国产黄片视频在线免费观看| 国产毛片a区久久久久| 啦啦啦中文免费视频观看日本| 久久久久国产网址| 欧美成人一区二区免费高清观看| 精品熟女少妇av免费看| 亚洲精品日本国产第一区| 日韩欧美 国产精品| 亚洲欧美一区二区三区国产| 国产三级在线视频| 男人狂女人下面高潮的视频| 亚洲欧洲国产日韩| 精品亚洲乱码少妇综合久久| 美女黄网站色视频| 美女被艹到高潮喷水动态| 少妇的逼好多水| 天堂av国产一区二区熟女人妻| 日本猛色少妇xxxxx猛交久久| 一级毛片aaaaaa免费看小| 欧美另类一区| 伦精品一区二区三区| 国产高清有码在线观看视频| 亚洲成色77777| 亚洲国产最新在线播放| 1000部很黄的大片| 国产乱来视频区| 日本与韩国留学比较| 嘟嘟电影网在线观看| 久久精品久久久久久噜噜老黄| 欧美成人午夜免费资源| 亚洲av免费在线观看| 日韩电影二区| 99九九线精品视频在线观看视频| www.av在线官网国产| 91久久精品国产一区二区成人| 搞女人的毛片| 综合色av麻豆| av黄色大香蕉| 国产免费一级a男人的天堂| 插逼视频在线观看| 精品午夜福利在线看| av女优亚洲男人天堂| 免费无遮挡裸体视频| 五月玫瑰六月丁香| 久久久久久久国产电影| 久久久亚洲精品成人影院| 亚洲欧美日韩卡通动漫| 免费观看精品视频网站| 日韩精品有码人妻一区| 禁无遮挡网站| 人人妻人人看人人澡| 美女黄网站色视频| 亚洲成人一二三区av| 国产免费一级a男人的天堂| 丝袜美腿在线中文| 小蜜桃在线观看免费完整版高清| 欧美bdsm另类| 最近2019中文字幕mv第一页| 特级一级黄色大片| 国产亚洲91精品色在线| 亚洲图色成人| 亚洲自偷自拍三级| kizo精华| 午夜福利视频精品| 欧美三级亚洲精品| 少妇熟女欧美另类| 欧美日韩精品成人综合77777| 特大巨黑吊av在线直播| 97超视频在线观看视频| 久久这里有精品视频免费| av在线天堂中文字幕| 久久精品国产亚洲av涩爱| 一区二区三区四区激情视频| 日本免费在线观看一区| 国产成人aa在线观看| 亚洲欧美日韩无卡精品| 天堂影院成人在线观看| 欧美另类一区| 免费大片黄手机在线观看| 18禁裸乳无遮挡免费网站照片| 中文资源天堂在线| 亚洲精品亚洲一区二区| 一本一本综合久久| 欧美日韩一区二区视频在线观看视频在线 | 又黄又爽又刺激的免费视频.| 最新中文字幕久久久久| 亚洲av电影不卡..在线观看| 日本午夜av视频| 97精品久久久久久久久久精品| 三级国产精品欧美在线观看| av在线播放精品| 99久久精品国产国产毛片| 国产精品一区二区三区四区免费观看| 久久久久久久午夜电影| 91精品一卡2卡3卡4卡| 久久久久久国产a免费观看| 国产乱人偷精品视频| 精品人妻视频免费看| 午夜激情福利司机影院| 美女被艹到高潮喷水动态| 日韩三级伦理在线观看| 国产探花在线观看一区二区| 亚洲熟妇中文字幕五十中出| 成人综合一区亚洲| 亚洲精品色激情综合| 天堂网av新在线| 精品久久久久久久久av| 天堂影院成人在线观看| 亚洲伊人久久精品综合| 女人被狂操c到高潮| 国产爱豆传媒在线观看| 午夜精品在线福利| 一级毛片黄色毛片免费观看视频| 欧美成人一区二区免费高清观看| 看免费成人av毛片| 综合色丁香网| 一个人看视频在线观看www免费| 偷拍熟女少妇极品色| 色5月婷婷丁香| or卡值多少钱| 91精品一卡2卡3卡4卡| or卡值多少钱| 欧美高清成人免费视频www| 一级a做视频免费观看| 成人综合一区亚洲| 天堂√8在线中文| 最近最新中文字幕大全电影3| 国产精品国产三级国产av玫瑰| 国产高清三级在线| 黑人高潮一二区| 夜夜爽夜夜爽视频| 国产一区二区亚洲精品在线观看| 免费观看精品视频网站| 嫩草影院新地址| 日本-黄色视频高清免费观看| 日韩一区二区三区影片| 久久久久久久久大av| 少妇的逼水好多| 亚洲精品久久午夜乱码| 国产成人精品久久久久久| 国产69精品久久久久777片| 国产成人免费观看mmmm| 久久亚洲国产成人精品v| 97热精品久久久久久| 日韩精品青青久久久久久| 亚洲不卡免费看| 欧美另类一区| 国产男女超爽视频在线观看| 99热全是精品| 我的女老师完整版在线观看| 国产精品人妻久久久久久| 18禁在线无遮挡免费观看视频| 91午夜精品亚洲一区二区三区| 免费黄网站久久成人精品| 男女边吃奶边做爰视频| 一个人看的www免费观看视频| 国产片特级美女逼逼视频| 日韩中字成人| 一级毛片aaaaaa免费看小| 日日摸夜夜添夜夜爱| 国产爱豆传媒在线观看| 欧美高清性xxxxhd video| 一边亲一边摸免费视频| 国产探花极品一区二区| 男人爽女人下面视频在线观看| 日韩中字成人| a级一级毛片免费在线观看| 男女视频在线观看网站免费| 日韩精品青青久久久久久| 亚洲精品一区蜜桃| 听说在线观看完整版免费高清| 国产av码专区亚洲av| 久久久久性生活片| 在线免费观看的www视频| 日韩成人伦理影院| 美女cb高潮喷水在线观看| 少妇人妻一区二区三区视频| 国产成人a区在线观看| 日日干狠狠操夜夜爽| 两个人的视频大全免费| 九九爱精品视频在线观看| av一本久久久久| 欧美区成人在线视频| 国产麻豆成人av免费视频| 日本午夜av视频| 国产欧美日韩精品一区二区| 日韩av在线大香蕉| 精品人妻偷拍中文字幕| 高清在线视频一区二区三区| 欧美一区二区亚洲| 亚洲成人av在线免费| 麻豆av噜噜一区二区三区| 国产精品一区二区三区四区久久| 国产精品一区二区三区四区免费观看| 国内精品宾馆在线| 26uuu在线亚洲综合色| 成人欧美大片| 一区二区三区免费毛片| 汤姆久久久久久久影院中文字幕 | 国产毛片a区久久久久| 51国产日韩欧美| 九九久久精品国产亚洲av麻豆| 欧美日本视频| 春色校园在线视频观看| 3wmmmm亚洲av在线观看| 亚洲人成网站在线观看播放| 午夜免费激情av| 亚洲精品成人久久久久久| 日韩电影二区| 国模一区二区三区四区视频| 只有这里有精品99| 免费观看精品视频网站| 大片免费播放器 马上看| 久久久午夜欧美精品| 一级爰片在线观看| 午夜福利视频精品| 午夜福利在线在线| 欧美成人午夜免费资源| 久久人人爽人人爽人人片va| 欧美日本视频| 久久久精品94久久精品| 精品久久久精品久久久| 国产精品久久久久久久电影| 亚洲国产色片| 国产精品福利在线免费观看| 99热这里只有精品一区| 久久久成人免费电影| 亚洲在线自拍视频| 嫩草影院精品99| 99久久九九国产精品国产免费| 禁无遮挡网站| 蜜臀久久99精品久久宅男| 免费黄网站久久成人精品| 亚洲18禁久久av|