• <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)
    高清在线视频一区二区三区| 亚洲精品aⅴ在线观看| 国产精品 国内视频| 日韩大片免费观看网站| 精品国产国语对白av| 国产精品国产三级国产专区5o| 在线观看人妻少妇| 亚洲成色77777| 69精品国产乱码久久久| 亚洲色图综合在线观看| 丰满饥渴人妻一区二区三| 青春草视频在线免费观看| 亚洲人成77777在线视频| 少妇被粗大猛烈的视频| 少妇猛男粗大的猛烈进出视频| 国产精品国产三级国产专区5o| 久久久精品免费免费高清| 国产免费一级a男人的天堂| 日日撸夜夜添| av不卡在线播放| 高清在线视频一区二区三区| 夫妻性生交免费视频一级片| 又大又黄又爽视频免费| 日日爽夜夜爽网站| 国产一级毛片在线| 伦精品一区二区三区| 欧美精品人与动牲交sv欧美| 亚州av有码| 一二三四中文在线观看免费高清| 人妻系列 视频| 欧美97在线视频| 免费高清在线观看日韩| 国产精品秋霞免费鲁丝片| 在线观看美女被高潮喷水网站| 午夜激情av网站| 久久久国产欧美日韩av| 肉色欧美久久久久久久蜜桃| 我的女老师完整版在线观看| 亚洲国产最新在线播放| av黄色大香蕉| 亚洲av二区三区四区| 国产午夜精品久久久久久一区二区三区| 亚洲欧洲日产国产| 日韩免费高清中文字幕av| 国产片特级美女逼逼视频| 美女cb高潮喷水在线观看| 18禁在线无遮挡免费观看视频| 午夜激情av网站| 久久久精品区二区三区| 高清在线视频一区二区三区| 久久av网站| 日本-黄色视频高清免费观看| 最近中文字幕2019免费版| 女性生殖器流出的白浆| 你懂的网址亚洲精品在线观看| 天天躁夜夜躁狠狠久久av| 久久99精品国语久久久| 极品人妻少妇av视频| 欧美少妇被猛烈插入视频| 国产高清有码在线观看视频| 免费观看的影片在线观看| 国产视频内射| 少妇被粗大猛烈的视频| 国产欧美另类精品又又久久亚洲欧美| xxxhd国产人妻xxx| 午夜福利视频精品| 久久青草综合色| 高清黄色对白视频在线免费看| 成人毛片60女人毛片免费| 22中文网久久字幕| 99精国产麻豆久久婷婷| 亚洲精品乱码久久久久久按摩| 最新中文字幕久久久久| 午夜激情久久久久久久| 黑人猛操日本美女一级片| 成人18禁高潮啪啪吃奶动态图 | 青青草视频在线视频观看| 一区二区三区免费毛片| 久久久久久人妻| 亚洲精品第二区| 如何舔出高潮| 亚洲第一av免费看| 99re6热这里在线精品视频| 欧美激情国产日韩精品一区| 日日摸夜夜添夜夜添av毛片| 高清黄色对白视频在线免费看| 午夜福利在线观看免费完整高清在| 国产熟女欧美一区二区| 亚洲av成人精品一二三区| 爱豆传媒免费全集在线观看| 久久ye,这里只有精品| 久久精品久久精品一区二区三区| 亚洲高清免费不卡视频| 最近中文字幕高清免费大全6| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久精品电影小说| 麻豆精品久久久久久蜜桃| 丝瓜视频免费看黄片| 婷婷成人精品国产| 乱码一卡2卡4卡精品| 91精品一卡2卡3卡4卡| 黑人猛操日本美女一级片| 桃花免费在线播放| 国产成人精品一,二区| 亚洲第一av免费看| 亚洲av不卡在线观看| 国国产精品蜜臀av免费| 亚洲综合色网址| 亚洲国产精品一区三区| 99九九在线精品视频| 国产黄色视频一区二区在线观看| 精品99又大又爽又粗少妇毛片| a级片在线免费高清观看视频| 九九在线视频观看精品| 国产精品99久久久久久久久| 狂野欧美白嫩少妇大欣赏| 美女内射精品一级片tv| 午夜激情av网站| 内地一区二区视频在线| 精品久久久噜噜| 精品视频人人做人人爽| 一二三四中文在线观看免费高清| 国产片内射在线| 女性被躁到高潮视频| 黄色视频在线播放观看不卡| 在线观看www视频免费| 少妇精品久久久久久久| 精品一品国产午夜福利视频| 日本黄色日本黄色录像| 99国产综合亚洲精品| 中文天堂在线官网| 国产极品粉嫩免费观看在线 | 亚洲精品乱久久久久久| 国产高清国产精品国产三级| 日本av手机在线免费观看| 日韩免费高清中文字幕av| 99久国产av精品国产电影| 日韩精品免费视频一区二区三区 | 色婷婷久久久亚洲欧美| 亚洲国产精品成人久久小说| 水蜜桃什么品种好| 久久精品久久久久久久性| 久久狼人影院| 国产av精品麻豆| 欧美人与性动交α欧美精品济南到 | 色婷婷av一区二区三区视频| 久久综合国产亚洲精品| 精品人妻熟女av久视频| 欧美精品高潮呻吟av久久| 亚洲人成77777在线视频| 精品久久久久久久久亚洲| 亚州av有码| 久久久久国产精品人妻一区二区| 夫妻性生交免费视频一级片| 多毛熟女@视频| 热re99久久精品国产66热6| 女性生殖器流出的白浆| 国产精品欧美亚洲77777| 欧美3d第一页| 边亲边吃奶的免费视频| 亚洲成人手机| 熟女av电影| 一区二区av电影网| 久久人人爽人人片av| 看免费成人av毛片| 黑人猛操日本美女一级片| 欧美 日韩 精品 国产| 肉色欧美久久久久久久蜜桃| 精品少妇久久久久久888优播| 新久久久久国产一级毛片| 久久久亚洲精品成人影院| 亚洲,一卡二卡三卡| 久久国产精品男人的天堂亚洲 | 日本猛色少妇xxxxx猛交久久| 国产视频内射| 久久99一区二区三区| 久久综合国产亚洲精品| 日本欧美视频一区| 亚洲美女视频黄频| 青春草国产在线视频| 午夜免费男女啪啪视频观看| 女性生殖器流出的白浆| 日韩熟女老妇一区二区性免费视频| 亚洲av男天堂| 日韩av不卡免费在线播放| 熟女av电影| 在线天堂最新版资源| 少妇丰满av| 国产精品秋霞免费鲁丝片| 国产视频首页在线观看| 大片免费播放器 马上看| 亚洲av欧美aⅴ国产| 欧美精品一区二区免费开放| 成人国语在线视频| 18在线观看网站| av黄色大香蕉| 天天操日日干夜夜撸| 国产精品无大码| 中文精品一卡2卡3卡4更新| 日韩免费高清中文字幕av| 又黄又爽又刺激的免费视频.| 精品国产一区二区三区久久久樱花| 久久精品久久久久久噜噜老黄| 91久久精品国产一区二区成人| 日韩制服骚丝袜av| 黄色毛片三级朝国网站| 中文字幕亚洲精品专区| 一区二区日韩欧美中文字幕 | 高清欧美精品videossex| 亚洲国产欧美在线一区| 男女国产视频网站| 亚洲高清免费不卡视频| 飞空精品影院首页| 亚洲欧洲精品一区二区精品久久久 | 欧美日韩一区二区视频在线观看视频在线| 亚洲国产av影院在线观看| 午夜福利视频精品| 亚洲人成网站在线播| 欧美激情极品国产一区二区三区 | 成人国产麻豆网| 毛片一级片免费看久久久久| 内地一区二区视频在线| 久久久国产精品麻豆| 日本wwww免费看| 国产免费一区二区三区四区乱码| 18禁观看日本| 欧美精品一区二区大全| 插阴视频在线观看视频| 亚洲精品美女久久av网站| 一边亲一边摸免费视频| 精品亚洲成国产av| 少妇丰满av| 97精品久久久久久久久久精品| 亚洲成色77777| 中文字幕最新亚洲高清| 日韩电影二区| 国产一区二区在线观看av| 久久久久国产精品人妻一区二区| 午夜福利在线观看免费完整高清在| 日本av手机在线免费观看| 青春草视频在线免费观看| 视频区图区小说| 免费黄色在线免费观看| 满18在线观看网站| av国产久精品久网站免费入址| 人妻 亚洲 视频| 亚洲国产av影院在线观看| 丰满饥渴人妻一区二区三| 日本-黄色视频高清免费观看| 91在线精品国自产拍蜜月| 男女无遮挡免费网站观看| 搡女人真爽免费视频火全软件| 少妇熟女欧美另类| 伊人亚洲综合成人网| 午夜福利视频精品| 伊人久久精品亚洲午夜| 国产日韩一区二区三区精品不卡 | 亚洲欧美成人综合另类久久久| 99热6这里只有精品| 一级毛片黄色毛片免费观看视频| av卡一久久| 久久亚洲国产成人精品v| 中国三级夫妇交换| 国产精品国产三级国产专区5o| 国产精品国产av在线观看| 日韩一本色道免费dvd| 欧美日韩成人在线一区二区| 国产熟女午夜一区二区三区 | 国产不卡av网站在线观看| 色网站视频免费| 欧美精品亚洲一区二区| 日本91视频免费播放| 免费观看在线日韩| 少妇人妻久久综合中文| 亚洲婷婷狠狠爱综合网| 乱人伦中国视频| 亚洲性久久影院| 欧美精品一区二区大全| 中文精品一卡2卡3卡4更新| 2021少妇久久久久久久久久久| 国产色婷婷99| 天天躁夜夜躁狠狠久久av| 一级二级三级毛片免费看| 蜜桃国产av成人99| 女性生殖器流出的白浆| 99久久中文字幕三级久久日本| 男人爽女人下面视频在线观看| 80岁老熟妇乱子伦牲交| 国产精品偷伦视频观看了| 国精品久久久久久国模美| videos熟女内射| 欧美精品一区二区大全| 成人18禁高潮啪啪吃奶动态图 | av在线老鸭窝| av福利片在线| 韩国高清视频一区二区三区| 成人毛片60女人毛片免费| 成人无遮挡网站| 99re6热这里在线精品视频| 97在线视频观看| 国产精品麻豆人妻色哟哟久久| 欧美bdsm另类| 欧美 日韩 精品 国产| 汤姆久久久久久久影院中文字幕| 老司机影院成人| 欧美日韩综合久久久久久| 一区二区三区四区激情视频| 久久久欧美国产精品| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产精品国产精品| 精品国产一区二区久久| 免费黄色在线免费观看| 大陆偷拍与自拍| 啦啦啦中文免费视频观看日本| 久久久欧美国产精品| 日本猛色少妇xxxxx猛交久久| 一个人免费看片子| 一级毛片 在线播放| 三上悠亚av全集在线观看| 中文欧美无线码| 久久久久网色| 永久免费av网站大全| 99久久精品国产国产毛片| 人人妻人人澡人人爽人人夜夜| 色视频在线一区二区三区| 熟女av电影| 国产片内射在线| 秋霞伦理黄片| 国产精品三级大全| av视频免费观看在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲少妇的诱惑av| 七月丁香在线播放| tube8黄色片| 桃花免费在线播放| 少妇精品久久久久久久| 777米奇影视久久| 丁香六月天网| 久久久久久伊人网av| 亚洲美女搞黄在线观看| 亚洲伊人久久精品综合| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产av影院在线观看| 久久久久久久国产电影| 少妇人妻 视频| 亚洲人成网站在线播| 3wmmmm亚洲av在线观看| 18禁观看日本| 少妇人妻久久综合中文| 久久久久久久久久久丰满| 看非洲黑人一级黄片| 最近的中文字幕免费完整| 伊人久久精品亚洲午夜| 久久国产精品男人的天堂亚洲 | 国产精品秋霞免费鲁丝片| 老司机亚洲免费影院| 午夜免费观看性视频| 精品国产乱码久久久久久小说| 国产片特级美女逼逼视频| 不卡视频在线观看欧美| 国产精品国产三级国产av玫瑰| 国产精品一二三区在线看| 久久久久久久久久人人人人人人| 国产精品三级大全| 欧美成人午夜免费资源| 亚洲精品视频女| 大片免费播放器 马上看| 久久这里有精品视频免费| www.色视频.com| 麻豆精品久久久久久蜜桃| 日本-黄色视频高清免费观看| 最近中文字幕高清免费大全6| 日韩视频在线欧美| 国产欧美日韩综合在线一区二区| 三级国产精品片| av在线观看视频网站免费| 天堂俺去俺来也www色官网| 日本与韩国留学比较| 国产69精品久久久久777片| 黄片无遮挡物在线观看| 亚洲精品国产av蜜桃| 亚洲欧洲国产日韩| 亚洲精品乱码久久久v下载方式| 一区二区av电影网| 丝袜喷水一区| 在线观看三级黄色| 欧美变态另类bdsm刘玥| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品,欧美精品| 日产精品乱码卡一卡2卡三| 一本—道久久a久久精品蜜桃钙片| 亚洲欧美中文字幕日韩二区| 亚洲av二区三区四区| 精品午夜福利在线看| 亚洲欧美日韩另类电影网站| 欧美最新免费一区二区三区| 在线观看免费日韩欧美大片 | 中国国产av一级| 最近2019中文字幕mv第一页| 看十八女毛片水多多多| 国产一区二区在线观看av| 韩国av在线不卡| 最近中文字幕2019免费版| 9色porny在线观看| 一级爰片在线观看| 免费观看a级毛片全部| 精品一区二区三区视频在线| 久久久久久久久久久免费av| 亚洲精品国产色婷婷电影| 精品人妻一区二区三区麻豆| 欧美日韩在线观看h| 成年人午夜在线观看视频| 久久久久久久久久久免费av| 久久午夜综合久久蜜桃| 国产成人精品一,二区| 又大又黄又爽视频免费| 国产无遮挡羞羞视频在线观看| 亚洲精品aⅴ在线观看| 最后的刺客免费高清国语| 久久精品国产亚洲av涩爱| av在线播放精品| 波野结衣二区三区在线| 99久国产av精品国产电影| 91久久精品电影网| 性色av一级| 欧美bdsm另类| av在线老鸭窝| 男女高潮啪啪啪动态图| 精品人妻熟女毛片av久久网站| 少妇猛男粗大的猛烈进出视频| 最近中文字幕2019免费版| 免费日韩欧美在线观看| 这个男人来自地球电影免费观看 | 日韩制服骚丝袜av| 中文字幕久久专区| 中国国产av一级| 国产欧美亚洲国产| 美女cb高潮喷水在线观看| 人妻夜夜爽99麻豆av| 欧美xxxx性猛交bbbb| 中文欧美无线码| 熟妇人妻不卡中文字幕| 亚洲av福利一区| 少妇被粗大的猛进出69影院 | 我要看黄色一级片免费的| 91久久精品电影网| 成人二区视频| 99热网站在线观看| 80岁老熟妇乱子伦牲交| 久久精品久久精品一区二区三区| 免费av不卡在线播放| 热99国产精品久久久久久7| 午夜福利网站1000一区二区三区| 午夜激情av网站| 在线观看免费日韩欧美大片 | 国产精品免费大片| 在线精品无人区一区二区三| 丰满少妇做爰视频| 久久热精品热| 亚洲av不卡在线观看| 99久久中文字幕三级久久日本| 成年女人在线观看亚洲视频| 人人妻人人添人人爽欧美一区卜| 亚洲人成77777在线视频| 在线 av 中文字幕| 91aial.com中文字幕在线观看| 丝瓜视频免费看黄片| 亚洲av福利一区| 考比视频在线观看| 日日爽夜夜爽网站| 成年美女黄网站色视频大全免费 | 草草在线视频免费看| 亚洲第一av免费看| 久久综合国产亚洲精品| 大陆偷拍与自拍| 国产 精品1| 哪个播放器可以免费观看大片| 97在线人人人人妻| 久久久精品94久久精品| 黄色一级大片看看| 亚洲精品日本国产第一区| 最新的欧美精品一区二区| 久久久久久久大尺度免费视频| 我的女老师完整版在线观看| 国产不卡av网站在线观看| 久久国产亚洲av麻豆专区| 制服丝袜香蕉在线| 另类精品久久| 色吧在线观看| 国产成人午夜福利电影在线观看| 热re99久久国产66热| 永久免费av网站大全| 亚洲成人一二三区av| 日本vs欧美在线观看视频| 精品人妻在线不人妻| 亚洲,一卡二卡三卡| 欧美xxⅹ黑人| 亚洲成人一二三区av| 99久久中文字幕三级久久日本| 亚洲人成网站在线观看播放| 男男h啪啪无遮挡| 国产成人精品无人区| 999精品在线视频| 久久久久久久久大av| 美女中出高潮动态图| 黄色怎么调成土黄色| 美女脱内裤让男人舔精品视频| 免费黄色在线免费观看| 日本wwww免费看| 久久久精品免费免费高清| 日日爽夜夜爽网站| 欧美精品一区二区大全| 中文字幕亚洲精品专区| 日韩av免费高清视频| 日韩,欧美,国产一区二区三区| 我要看黄色一级片免费的| 大香蕉久久成人网| 老女人水多毛片| 一级二级三级毛片免费看| 精品国产一区二区三区久久久樱花| 少妇人妻精品综合一区二区| 国产精品秋霞免费鲁丝片| 考比视频在线观看| 18禁在线无遮挡免费观看视频| 男女啪啪激烈高潮av片| 国产深夜福利视频在线观看| 激情五月婷婷亚洲| 国产男女超爽视频在线观看| 99热6这里只有精品| 菩萨蛮人人尽说江南好唐韦庄| 精品亚洲成国产av| 精品午夜福利在线看| 成年女人在线观看亚洲视频| 亚洲熟女精品中文字幕| 久久午夜综合久久蜜桃| 一区二区日韩欧美中文字幕 | 99热这里只有精品一区| 91国产中文字幕| 性色av一级| 汤姆久久久久久久影院中文字幕| 亚洲国产色片| 两个人的视频大全免费| 狂野欧美激情性bbbbbb| 汤姆久久久久久久影院中文字幕| 妹子高潮喷水视频| 久久精品久久久久久噜噜老黄| 久久精品国产鲁丝片午夜精品| 我的女老师完整版在线观看| 日韩视频在线欧美| 欧美少妇被猛烈插入视频| 久久亚洲国产成人精品v| 久久久久久久精品精品| 九草在线视频观看| 国产免费一区二区三区四区乱码| 久久人妻熟女aⅴ| 精品人妻熟女av久视频| 久久99精品国语久久久| 久热久热在线精品观看| 美女国产视频在线观看| 丝袜喷水一区| 日本vs欧美在线观看视频| 在线观看免费视频网站a站| 91精品伊人久久大香线蕉| 午夜日本视频在线| 丁香六月天网| 日本爱情动作片www.在线观看| 人妻夜夜爽99麻豆av| 精品99又大又爽又粗少妇毛片| 男女国产视频网站| 国产高清不卡午夜福利| 久久精品久久久久久久性| 777米奇影视久久| 国产精品一区二区三区四区免费观看| 国产老妇伦熟女老妇高清| 成人国产麻豆网| 亚洲精品一二三| 2021少妇久久久久久久久久久| 美女主播在线视频| 久久婷婷青草| 一区二区三区免费毛片| 中文欧美无线码| 精品视频人人做人人爽| 亚洲高清免费不卡视频| 成人手机av| 男女无遮挡免费网站观看| 亚洲高清免费不卡视频| av播播在线观看一区| 久久狼人影院| 久久午夜福利片| 国产黄频视频在线观看| 久久狼人影院| 午夜91福利影院| 日本色播在线视频| av视频免费观看在线观看| 亚洲欧美清纯卡通| 日韩中字成人| 国产精品人妻久久久久久| 亚洲四区av| 搡女人真爽免费视频火全软件| 亚洲无线观看免费| 曰老女人黄片| 精品久久久久久电影网| 亚洲无线观看免费| 国产精品久久久久久精品电影小说| 菩萨蛮人人尽说江南好唐韦庄| 国产精品99久久久久久久久| 丝袜美足系列| 我的老师免费观看完整版| 亚洲不卡免费看| 一级二级三级毛片免费看| 天天操日日干夜夜撸| 91精品国产九色| 激情五月婷婷亚洲| 亚洲经典国产精华液单| 天堂俺去俺来也www色官网|