• <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级毛片在线看网站| 老司机福利观看| 亚洲avbb在线观看| 白带黄色成豆腐渣| 久久精品91蜜桃| 国内揄拍国产精品人妻在线| 亚洲欧美一区二区三区黑人| aaaaa片日本免费| 天堂av国产一区二区熟女人妻| 夜夜看夜夜爽夜夜摸| 在线观看一区二区三区| 岛国在线免费视频观看| 校园春色视频在线观看| 免费av不卡在线播放| 一级a爱片免费观看的视频| 国产精品自产拍在线观看55亚洲| 国产97色在线日韩免费| 亚洲一区二区三区不卡视频| 精华霜和精华液先用哪个| 亚洲九九香蕉| 色视频www国产| 人人妻人人看人人澡| 国产高潮美女av| 中文字幕人妻丝袜一区二区| 欧美xxxx黑人xx丫x性爽| 色综合欧美亚洲国产小说| 一本综合久久免费| 五月伊人婷婷丁香| 亚洲黑人精品在线| av黄色大香蕉| 小说图片视频综合网站| 国产精品一区二区三区四区久久| 国产淫片久久久久久久久 | 亚洲精品粉嫩美女一区| 后天国语完整版免费观看| 两个人看的免费小视频| 999久久久精品免费观看国产| 美女午夜性视频免费| 婷婷六月久久综合丁香| 俺也久久电影网| 麻豆成人午夜福利视频| av黄色大香蕉| 欧美三级亚洲精品| 亚洲精品粉嫩美女一区| 美女高潮的动态| 日日摸夜夜添夜夜添小说| 午夜两性在线视频| 国产午夜精品久久久久久| av女优亚洲男人天堂 | 国产精品精品国产色婷婷| 最新中文字幕久久久久 | 91老司机精品| 999久久久国产精品视频| 国产精品久久久人人做人人爽| 夜夜爽天天搞| 男女那种视频在线观看| 亚洲专区中文字幕在线| 国内毛片毛片毛片毛片毛片| 午夜福利成人在线免费观看| 99国产精品一区二区三区| 97碰自拍视频| 免费大片18禁| 日韩欧美在线二视频| 国产精品1区2区在线观看.| 精品国产三级普通话版| 国产精华一区二区三区| 国产1区2区3区精品| 国产乱人视频| 嫩草影院精品99| 亚洲欧美激情综合另类| 免费无遮挡裸体视频| 久久这里只有精品中国| 欧美乱妇无乱码| 亚洲av成人精品一区久久| 动漫黄色视频在线观看| 久久久国产精品麻豆| 中文字幕精品亚洲无线码一区| 亚洲中文日韩欧美视频| 欧美又色又爽又黄视频| 亚洲aⅴ乱码一区二区在线播放| 熟女少妇亚洲综合色aaa.| 激情在线观看视频在线高清| 欧美黄色片欧美黄色片| 亚洲人成电影免费在线| 无限看片的www在线观看| 少妇人妻一区二区三区视频| 这个男人来自地球电影免费观看| 国产av麻豆久久久久久久| 91字幕亚洲| a级毛片a级免费在线| 免费在线观看成人毛片| 搞女人的毛片| 午夜免费观看网址| 99久久国产精品久久久| 亚洲最大成人中文| 丰满人妻一区二区三区视频av | 国产日本99.免费观看| 国产伦在线观看视频一区| 午夜影院日韩av| а√天堂www在线а√下载| 成人三级做爰电影| 国产精品1区2区在线观看.| 美女黄网站色视频| 精品一区二区三区四区五区乱码| 少妇熟女aⅴ在线视频| 日本 av在线| 最近最新免费中文字幕在线| 级片在线观看| 伊人久久大香线蕉亚洲五| 久久久精品欧美日韩精品| 99久久精品国产亚洲精品| 国产精品亚洲美女久久久| 18禁黄网站禁片免费观看直播| 国产成人精品久久二区二区91| 99久久成人亚洲精品观看| 亚洲电影在线观看av| 亚洲熟妇熟女久久| 国产精品亚洲av一区麻豆| 天天一区二区日本电影三级| 99热只有精品国产| 午夜激情欧美在线| 白带黄色成豆腐渣| 99热只有精品国产| 欧美黄色淫秽网站| 国产亚洲精品久久久久久毛片| 亚洲熟妇中文字幕五十中出| 免费在线观看亚洲国产| 婷婷亚洲欧美| 日韩三级视频一区二区三区| 成年免费大片在线观看| 亚洲人成网站在线播放欧美日韩| 最近视频中文字幕2019在线8| 国产v大片淫在线免费观看| 一个人看的www免费观看视频| 97超视频在线观看视频| 国内精品久久久久久久电影| av福利片在线观看| 久久精品影院6| xxxwww97欧美| 欧美日韩综合久久久久久 | 精品久久久久久久久久久久久| 精品国产亚洲在线| 90打野战视频偷拍视频| 校园春色视频在线观看| 国产乱人视频| 91在线精品国自产拍蜜月 | 欧美三级亚洲精品| 国产精品爽爽va在线观看网站| 久久精品国产99精品国产亚洲性色| 欧美日韩精品网址| 国产成人福利小说| 两人在一起打扑克的视频| 婷婷亚洲欧美| 久久久国产精品麻豆| 老司机福利观看| av天堂中文字幕网| 人人妻人人澡欧美一区二区| 国产一区在线观看成人免费| 我的老师免费观看完整版| 亚洲国产欧洲综合997久久,| 午夜激情福利司机影院| 久久精品国产综合久久久| 国产男靠女视频免费网站| 中文字幕人成人乱码亚洲影| 我要搜黄色片| 别揉我奶头~嗯~啊~动态视频| 精品国内亚洲2022精品成人| 小蜜桃在线观看免费完整版高清| 亚洲美女黄片视频| 欧美zozozo另类| 欧美一级毛片孕妇| 免费无遮挡裸体视频| 国产私拍福利视频在线观看| 狠狠狠狠99中文字幕| 国产精品亚洲av一区麻豆| 老司机深夜福利视频在线观看| 精品人妻1区二区| 露出奶头的视频| 狂野欧美白嫩少妇大欣赏| 人妻丰满熟妇av一区二区三区| 国产精品久久久久久久电影 | 日本与韩国留学比较| 成人特级av手机在线观看| 最近在线观看免费完整版| 亚洲av成人精品一区久久| 精品国产美女av久久久久小说| 中文字幕熟女人妻在线| 亚洲人成电影免费在线| 亚洲 欧美一区二区三区| 精品一区二区三区视频在线观看免费| 欧美午夜高清在线| 亚洲精华国产精华精| 亚洲avbb在线观看| 曰老女人黄片| 日韩 欧美 亚洲 中文字幕| 国产黄色小视频在线观看| e午夜精品久久久久久久| 欧美高清成人免费视频www| 亚洲 国产 在线| 国产午夜精品久久久久久| 国产真实乱freesex| 国内精品久久久久精免费| 久久久久久人人人人人| 1000部很黄的大片| 黑人巨大精品欧美一区二区mp4| 婷婷六月久久综合丁香| 欧美日韩中文字幕国产精品一区二区三区| 亚洲乱码一区二区免费版| 欧美黄色片欧美黄色片| 女警被强在线播放| 中文字幕av在线有码专区| 国产三级中文精品| 亚洲国产欧美人成| 99在线视频只有这里精品首页| 亚洲欧洲精品一区二区精品久久久| 三级男女做爰猛烈吃奶摸视频| 色在线成人网| 亚洲av日韩精品久久久久久密| or卡值多少钱| 亚洲成人中文字幕在线播放| 久久久国产精品麻豆| 午夜福利视频1000在线观看| 婷婷丁香在线五月| 不卡av一区二区三区| 久久国产乱子伦精品免费另类| 99国产综合亚洲精品| 国产视频内射| 超碰成人久久| 男人舔女人下体高潮全视频| 欧美激情在线99| 色av中文字幕| 欧美日韩综合久久久久久 | 国产一区二区在线观看日韩 | 叶爱在线成人免费视频播放| 国产v大片淫在线免费观看| 亚洲在线自拍视频| 日韩 欧美 亚洲 中文字幕| 国产精品美女特级片免费视频播放器 | 精品国产三级普通话版| 首页视频小说图片口味搜索| 欧美+亚洲+日韩+国产| 一边摸一边抽搐一进一小说| 亚洲成av人片在线播放无| 亚洲,欧美精品.| 亚洲欧美精品综合久久99| 成人av一区二区三区在线看| 毛片女人毛片| 一个人看的www免费观看视频| 国产一区二区三区在线臀色熟女| 日韩大尺度精品在线看网址| 久久中文看片网| 欧美一级毛片孕妇| 亚洲自拍偷在线| 国产淫片久久久久久久久 | xxx96com| 欧美av亚洲av综合av国产av| av片东京热男人的天堂| 18禁观看日本| 麻豆成人午夜福利视频| 国产99白浆流出| 夜夜爽天天搞| 欧美不卡视频在线免费观看| 91av网站免费观看| 真人做人爱边吃奶动态| 日韩三级视频一区二区三区| 午夜精品一区二区三区免费看| 国产高清激情床上av| 老司机午夜十八禁免费视频| 伊人久久大香线蕉亚洲五| 美女cb高潮喷水在线观看 | tocl精华| 亚洲专区国产一区二区| 一个人观看的视频www高清免费观看 | 欧美又色又爽又黄视频| 久久九九热精品免费| 老司机在亚洲福利影院| 99国产精品一区二区蜜桃av| 特级一级黄色大片| 我要搜黄色片| 欧美日韩国产亚洲二区| 国产精品自产拍在线观看55亚洲| 一区二区三区国产精品乱码| 久久精品影院6| 国产人伦9x9x在线观看| 久久这里只有精品中国| 久久精品国产99精品国产亚洲性色| 18禁观看日本| 一进一出抽搐动态| 一区二区三区高清视频在线| 1024手机看黄色片| 麻豆久久精品国产亚洲av| 欧美成人性av电影在线观看| 欧美绝顶高潮抽搐喷水| 91av网站免费观看| 欧美乱色亚洲激情| 天天添夜夜摸| 欧美乱码精品一区二区三区| 噜噜噜噜噜久久久久久91| 欧美一级a爱片免费观看看| 少妇丰满av| e午夜精品久久久久久久| 亚洲专区中文字幕在线| 女警被强在线播放| 91字幕亚洲| 最近在线观看免费完整版| 熟女电影av网| 午夜免费激情av| 成年版毛片免费区| 99久久无色码亚洲精品果冻| 免费无遮挡裸体视频| 搡老妇女老女人老熟妇| 亚洲va日本ⅴa欧美va伊人久久| 综合色av麻豆| 看黄色毛片网站| 久久久水蜜桃国产精品网| 国产综合懂色| 久久精品国产清高在天天线| 亚洲天堂国产精品一区在线| 免费电影在线观看免费观看| 中文在线观看免费www的网站| 亚洲美女视频黄频| 99热这里只有精品一区 | 美女 人体艺术 gogo| 黄色 视频免费看| 午夜福利视频1000在线观看| 999久久久精品免费观看国产| 18禁裸乳无遮挡免费网站照片| 国产精品99久久99久久久不卡| 两个人视频免费观看高清| 国产午夜精品论理片| 精品久久久久久久毛片微露脸| 免费在线观看日本一区| 国产三级黄色录像| 可以在线观看毛片的网站| 在线免费观看的www视频| 亚洲国产精品久久男人天堂| 午夜免费成人在线视频| 国产视频内射| 少妇人妻一区二区三区视频| 亚洲欧美精品综合久久99| 在线视频色国产色| 成人亚洲精品av一区二区| 日韩欧美精品v在线| 很黄的视频免费| 日韩欧美国产一区二区入口| xxxwww97欧美| 久久久国产欧美日韩av| 免费在线观看成人毛片| 村上凉子中文字幕在线| 少妇的丰满在线观看| 欧美日本亚洲视频在线播放| 成年女人毛片免费观看观看9| 女人高潮潮喷娇喘18禁视频| 国产成人欧美在线观看| 国产不卡一卡二| 校园春色视频在线观看| 男女下面进入的视频免费午夜| 99久久综合精品五月天人人| 国产一区在线观看成人免费| 午夜成年电影在线免费观看| 国产精品一及| 久久久精品大字幕| 床上黄色一级片| 美女被艹到高潮喷水动态| 变态另类丝袜制服| 在线国产一区二区在线| 欧美性猛交╳xxx乱大交人| 国产亚洲精品综合一区在线观看| 极品教师在线免费播放| 天天一区二区日本电影三级| 国产精品久久久久久人妻精品电影| 亚洲国产精品合色在线| 中亚洲国语对白在线视频| 成年版毛片免费区| 18禁裸乳无遮挡免费网站照片| 午夜两性在线视频| 国产淫片久久久久久久久 | 亚洲成av人片在线播放无| 欧美黄色片欧美黄色片| 成熟少妇高潮喷水视频| 日本黄大片高清| 欧美性猛交黑人性爽| 日本在线视频免费播放| 麻豆久久精品国产亚洲av| 91久久精品国产一区二区成人 | 亚洲av电影不卡..在线观看| 日韩精品青青久久久久久| 午夜福利在线观看吧| 99国产精品99久久久久| 亚洲国产精品sss在线观看| 久9热在线精品视频| 欧美不卡视频在线免费观看| 综合色av麻豆| 国产在线精品亚洲第一网站| 亚洲欧美日韩无卡精品| 亚洲国产色片| 黄色成人免费大全| 俄罗斯特黄特色一大片| 亚洲av成人精品一区久久| 国产一区二区在线观看日韩 | 又大又爽又粗| 在线十欧美十亚洲十日本专区| 欧美日韩中文字幕国产精品一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 99国产极品粉嫩在线观看| 国产亚洲精品久久久com| 日韩国内少妇激情av| 琪琪午夜伦伦电影理论片6080| 午夜久久久久精精品| 在线视频色国产色| 国产成人啪精品午夜网站| 啦啦啦免费观看视频1| 午夜福利成人在线免费观看| 国产私拍福利视频在线观看| 日日夜夜操网爽| or卡值多少钱| 欧美激情久久久久久爽电影| 性色avwww在线观看| 国产一区二区三区视频了| 亚洲欧美日韩高清专用| av天堂在线播放| 国产成人欧美在线观看| 久久久成人免费电影| 真实男女啪啪啪动态图| 19禁男女啪啪无遮挡网站| 一级毛片高清免费大全| 精品免费久久久久久久清纯| 久久国产精品影院| 99精品久久久久人妻精品| 日韩精品青青久久久久久| 欧美黄色淫秽网站| 国产蜜桃级精品一区二区三区| 黑人欧美特级aaaaaa片| 综合色av麻豆| 成人国产综合亚洲| 日本一二三区视频观看| ponron亚洲| e午夜精品久久久久久久| 99re在线观看精品视频| 两个人看的免费小视频| a级毛片a级免费在线| 搞女人的毛片| 这个男人来自地球电影免费观看| 久久久久久九九精品二区国产| 亚洲成a人片在线一区二区| 午夜福利欧美成人| 亚洲av成人精品一区久久| 国产亚洲欧美在线一区二区| 欧美乱色亚洲激情| 在线观看午夜福利视频| 悠悠久久av| 桃色一区二区三区在线观看| 欧美成人免费av一区二区三区| 12—13女人毛片做爰片一| 国产精华一区二区三区| 手机成人av网站| 免费大片18禁| 亚洲av片天天在线观看| 天堂√8在线中文| 黑人欧美特级aaaaaa片| 美女午夜性视频免费| 亚洲成av人片免费观看| 曰老女人黄片| av在线天堂中文字幕| www.熟女人妻精品国产| 美女免费视频网站| 99热6这里只有精品| 亚洲专区国产一区二区| 国产精品女同一区二区软件 | 午夜视频精品福利| 国产精品久久视频播放| 免费在线观看视频国产中文字幕亚洲| 亚洲精品国产精品久久久不卡| 制服人妻中文乱码| 一区福利在线观看| 国产亚洲欧美98| 少妇丰满av| 精品久久久久久久久久免费视频| 午夜福利免费观看在线| 网址你懂的国产日韩在线| 真人一进一出gif抽搐免费| 久久欧美精品欧美久久欧美| 日韩欧美在线二视频| 天堂√8在线中文| 国产伦在线观看视频一区| 国产精品久久视频播放| 久久精品人妻少妇| 91麻豆av在线| 十八禁网站免费在线| 一本综合久久免费| 高清毛片免费观看视频网站| 嫩草影视91久久| 人妻久久中文字幕网| 婷婷六月久久综合丁香| 一级毛片高清免费大全| 男女视频在线观看网站免费| 久久久久久九九精品二区国产| 精品久久蜜臀av无| 在线免费观看的www视频| netflix在线观看网站| 日韩av在线大香蕉| 国产精品综合久久久久久久免费| www国产在线视频色| 久久中文看片网| 欧美日本视频| 欧美xxxx黑人xx丫x性爽| 999久久久国产精品视频| 成人亚洲精品av一区二区| 天堂动漫精品| 亚洲专区字幕在线| 给我免费播放毛片高清在线观看| 午夜福利18| 久久国产精品人妻蜜桃| 免费看美女性在线毛片视频| av天堂中文字幕网| 亚洲人成伊人成综合网2020| 日韩成人在线观看一区二区三区| 制服人妻中文乱码| 俄罗斯特黄特色一大片| 亚洲国产高清在线一区二区三| 99久久精品国产亚洲精品| 一进一出好大好爽视频| 国产精品电影一区二区三区| av视频在线观看入口| 国产一级毛片七仙女欲春2| 国产视频内射| 99热只有精品国产| 欧美绝顶高潮抽搐喷水| 好男人在线观看高清免费视频| 色在线成人网| 亚洲在线观看片| 日本免费a在线| 日韩中文字幕欧美一区二区| 美女高潮的动态| 岛国视频午夜一区免费看| 国产成人精品久久二区二区免费| 国产精品 国内视频| 男人舔女人的私密视频| 午夜成年电影在线免费观看| 成人特级av手机在线观看| 免费av不卡在线播放| 精品一区二区三区四区五区乱码| 香蕉av资源在线| 国产亚洲精品av在线| 久久久久九九精品影院| 亚洲国产精品sss在线观看| 午夜精品一区二区三区免费看| 日本a在线网址| 不卡一级毛片| 色av中文字幕| 亚洲av日韩精品久久久久久密| av福利片在线观看| 桃色一区二区三区在线观看| 88av欧美| 久久久久国产精品人妻aⅴ院| 国产真实乱freesex| 精品国产三级普通话版| 亚洲av成人一区二区三| 亚洲五月婷婷丁香| 搡老岳熟女国产| 脱女人内裤的视频| 亚洲国产欧美一区二区综合| av天堂在线播放| 国产在线精品亚洲第一网站| 日日干狠狠操夜夜爽| 999精品在线视频| 美女扒开内裤让男人捅视频| 国产精品爽爽va在线观看网站| 国产精品亚洲一级av第二区| 久久久国产成人精品二区| 久久精品影院6| h日本视频在线播放| 在线免费观看不下载黄p国产 | 国产亚洲精品av在线| 国产亚洲精品久久久com| 国产真人三级小视频在线观看| 免费一级毛片在线播放高清视频| 在线观看午夜福利视频| 欧美av亚洲av综合av国产av| 日日夜夜操网爽| 欧美xxxx黑人xx丫x性爽| 村上凉子中文字幕在线| 啦啦啦韩国在线观看视频| 91在线精品国自产拍蜜月 | 成在线人永久免费视频| 成年女人看的毛片在线观看| 欧美成人免费av一区二区三区| 男女床上黄色一级片免费看| 手机成人av网站| 欧美色欧美亚洲另类二区| 亚洲国产欧洲综合997久久,| 老司机在亚洲福利影院| 成年免费大片在线观看| 亚洲熟妇熟女久久| 又紧又爽又黄一区二区| 成年女人毛片免费观看观看9| 国产精品九九99| 午夜两性在线视频| 精品日产1卡2卡| 亚洲av电影在线进入| 欧美日韩黄片免| 久久伊人香网站| 观看美女的网站| 精品电影一区二区在线| 国产 一区 欧美 日韩| 亚洲av熟女| 成人国产一区最新在线观看| 日本一本二区三区精品| 免费搜索国产男女视频| 村上凉子中文字幕在线| 丁香六月欧美| 日韩高清综合在线| 国产又色又爽无遮挡免费看| 看黄色毛片网站| 真人一进一出gif抽搐免费|