全 文 : Guihaia Apr. 2016ꎬ 36(4):471-478
http: / / journal.gxzw.gxib.cn
http: / / www.guihaia-journal.com
DOI: 10.11931 / guihaia.gxzw201406013
石鹏ꎬ曹红星ꎬ李东霞ꎬ等. 油棕等热带植物 DXS基因的生物信息学分析 [J]. 广西植物ꎬ 2016ꎬ 36(4):471-478
SHI PꎬCAO HXꎬLI DXꎬet al. Bioinformatics analysis of DXS gene from six tropical plants including oil palm (Elaeis guineensis Jacq.)[J]. Guihaiaꎬ 2016ꎬ
36(4):471-478
油棕等热带植物 DXS基因的生物信息学分析
石 鹏ꎬ 曹红星ꎬ 李东霞ꎬ 王 永ꎬ 雷新涛∗
( 中国热带农业科学院 椰子研究所 海南省热带油料作物生物学重点实验室ꎬ 海南 文昌 571339 )
摘 要: 油棕等热带植物含有丰富的胡萝卜素和维生素 E等类异戊二烯物质ꎬ类胡萝卜素和甾醇等类异戊二
烯物质在植物生命活动中扮演重要角色ꎬ并且对保护人类健康具有重要意义ꎬMEP 途径是合成类异戊二烯的
重要途径之一ꎮ DXS是 MEP 途径中的第一个限速酶ꎬ其功能在油棕等热带植物中极其保守ꎮ 为了弄清油棕
等热带植物 DXS的结构和功能特点ꎬ该研究利用生物信息学工具和软件对以油棕等热带植物类异戊二烯合
成关键基因 DXS为对象ꎬ进行核酸和氨基酸序列的理化性质、蛋白质结构以及功能结构域等分析ꎬ探讨了不
同物种间的亲缘关系ꎮ 结果表明:DXS基因起始密码子均为 ATGꎬ终止密码子则分为 TAG、TAA 和 TGAꎬDXS
蛋白质属于不具有信号肽的亲水性蛋白ꎬ可能作为转运蛋白在叶绿体基质中发挥作用ꎬ未发现明显的跨膜结
构域ꎬ磷酸化位点有 36个ꎬ其中丝氨酸、苏氨酸和酪氨酸位点分别为 17、11和 8个ꎬ无规则卷曲和 α ̄螺旋是蛋
白质二级结构主要的结构元件ꎬ三级结构预测具有 DXS 酶特征ꎬ硫胺素焦磷酸盐结合位点和 PLN02582 保守
结构域ꎬ不同植物 DXS功能结构域非常保守ꎬ可以作为判断不同物种间亲缘关系的重要依据ꎮ 该研究结果为
油棕等热带植物 DXS的结构、功能分析和利用提供了进一步的信息ꎬ为其品质性状分子机制研究及遗传改良
奠定了基础ꎮ
关键词: 油棕ꎬ DXS基因ꎬ 生物信息学ꎬ 系统进化树
中图分类号: Q811.4 文献标识码: A 文章编号: 1000 ̄3142(2016)04 ̄0471 ̄08
Bioinformatics analysis of DXS gene from six tropical
plants including oil palm (Elaeis guineensis Jacq.)
SHI Pengꎬ CAO Hong ̄Xingꎬ LI Dong ̄Xiaꎬ WANG Yongꎬ LEI Xin ̄Tao∗
( Hainan Key Biological Laboratory of Tropical Oil Cropsꎬ Coconut Research Instituteꎬ
Chinese Academy of Tropical Agricultural Sciencesꎬ 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 im ̄
portant 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 tropi ̄
cal plants. In order to understand the structure and function characteristics of DXS from tropical plants as oil palmꎬ DXS
收稿日期: 2014 ̄12 ̄19 修回日期: 2015 ̄02 ̄15
基金项目: 国家林业局“948”项目(No.2013 ̄4 ̄40)[Supported by the Plan of Introducing the Advanced Technology for Agriculture of State Forestry Ad ̄
ministration(No.2013 ̄4 ̄40)]ꎮ
作者简介: 石鹏(1987 ̄)ꎬ男ꎬ湖南安乡人ꎬ硕士ꎬ研究实习员ꎬ主要从事油棕种质资源评价等研究ꎬ(E ̄mail)624811593@ qq.comꎮ
∗通讯作者: 雷新涛ꎬ博士ꎬ研究员ꎬ主要从事热带油料作物遗传育种等研究ꎬ(E ̄mail)xtlei@ 263.comꎮ
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 tyro ̄
sine 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 impor ̄
tant 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
油棕(Elaeis guineensis Jacq.)属棕榈科多年生
乔木ꎬ是热带产油植物ꎬ果实中类胡萝卜素的含量是
决定油棕品质的一个重要指标ꎮ 粗棕榈油中含有
600 mgkg ̄1类胡萝卜素ꎬ其中 80%都是 β ̄胡萝卜
素和 α ̄胡萝卜素ꎬ这也是粗棕榈油呈黄色或红棕色
的原因(雷新涛和曹红星ꎬ2013)ꎮ 巴西橡胶、可可、
白木香、杜仲和麻疯树都是重要的热带植物ꎬ对油棕
等热带植物进行 DXS的生物信息学分析ꎬ对于了解
这些重要植物的 DXS蛋白质结构和功能具有重要意
义ꎬ同时也为分析这些热带植物间进化关系提供
依据ꎮ
类异戊二烯是最大的一类天然产物ꎬ目前为止
已知的有 3 万多种 ( Sacchettini & Poulterꎬ1997)ꎮ
一些类异戊二烯在细胞中有着重要作用ꎬ如甾醇是
细胞膜结构的重要组分(Erickꎬ2008)ꎬ长萜醇在蛋
白质糖基化中担任糖基载体 ( Vincent & Dirkꎬ
2011)ꎬ泛素和质体醌在光合作用和呼吸作用中扮
演电子转运体的角色(Keisuke et alꎬ2010)ꎬ脱落酸、
细胞分裂素、赤霉素和甾醇作为激素控制着生理和
发育过程( Jonathan & Neilꎬ2001)ꎮ 尽管各种类异
戊二烯的结构和功能存在差异ꎬ但是所有的都由异
戊烯焦磷酸(IPP)和它的同分异构体二甲基丙烯焦
磷酸(DMAPP)构成ꎮ IPP 和 DMAPP 能通过两个不
同生物合成路线:甲羟戊酸(MVA)途径ꎬ其主要存
在真菌、酵母和动物细胞质中ꎻ2 ̄C ̄甲基 ̄D ̄赤藻糖
醇 ̄4 ̄磷酸(MEP)途径ꎬ其存在原生生物ꎬ大多数细
菌ꎬ绿藻和高等植物质体中 ( Lombard & Moreiraꎬ
2011ꎻWanke et alꎬ2001ꎻVranova et alꎬ2013)ꎮ MEP
途径的第一步是丙酮酸和 3 ̄磷酸甘油醛在 1 ̄脱氧
木酮糖 ̄5 ̄磷酸合成酶(DXS)的催化下生成 1 ̄脱氧
木酮糖 ̄5 ̄磷酸(DXP)ꎮ DXS是一种新的转酮酶ꎬ由
DXS基因翻译而来(Elizabeth et alꎬ2009ꎻHan et alꎬ
2013)ꎮ 在植物和大肠杆菌中ꎬ由此反应产生的
DXP 被用来在质体中生成 IPPꎮ DXS基因最初在大
肠杆菌中被克隆和分析ꎬ其它细菌和高等植物中随
后也找到了同源基因(Cordoba et alꎬ2011ꎻZhang et
alꎬ2009ꎻGong et alꎬ2006ꎻLuisa et alꎬ1998ꎻKuzuyama
et alꎬ2000ꎻHahn et alꎬ2001)ꎮ 胡萝卜素和叶黄素统
称为类胡萝卜素ꎬ是所有进行光合作用生物中的色
素ꎮ 在高等植物中ꎬ类胡萝卜素在绿色组织光合作
用中扮演重要角色ꎮ 其中ꎬβ ̄胡萝卜素是维生素 A
的前体ꎬ具有抗癌的作用ꎬ其缺乏会导致失明和早逝
(Rao et alꎬ2007ꎻShan & Liꎬ2008ꎻAguila & Manuelꎬ
2012)ꎮ 由于类胡萝卜素的重要性ꎬ高等植物中编
码类胡萝卜素生物合成的酶基因成为了增加类胡萝
卜素生物合成的潜在位点ꎮ 有文献报道证实 DXS 是
途径中的限速酶ꎬMEP 途径中的基因高表达与植物
中更高含量的类胡萝卜素有直接关系(Pengꎬ2008)ꎮ
生物信息学兴起于 21世纪ꎬ利用信息学工具分
析核酸和蛋白质等大分子ꎬ大规模挖掘这些大分子
的生物学意义(Luscombe et alꎬ2001)ꎮ 目前ꎬ植物
中开展 DXS的生物信息学分析还不多ꎬ尤其在油棕
等热带植物中还很少报道ꎮ 本文主要以油棕等热带
植物的 DXS核酸和蛋白质序列为分析对象ꎬ利用大
量的生物信息学工具挖掘其序列的生物学意义ꎬ为
进一步研究油棕等热带植物 DXS 基因功能和类胡
萝卜素代谢过程奠定基础ꎮ
1 材料与方法
1.1 数据来源
核酸和蛋白质序列均来源于NCBI ( National
274 广 西 植 物 36卷
表 1 生物信息学分析使用的工具和软件
Table 1 Tools and softwares used by bioinformatics analysis
工具或软件
Tool or software 功能 Function 网址 Website
ProtParam 蛋白质理化性质分析
Physicochemical property analysis
http: / / web.expasy.org / protparam /
ORF Finder 查找开放阅读框
Open reading frame
http: / / www.ncbi.nlm.nih.gov / gorf / gorf.html
ProtScale 蛋白质亲水性和疏水性分析
Hydrophobicity or hydrophilicity analysis
http: / / web.expasy.org / protscale /
SOPMA 蛋白质二级结构预测
Secondary structure prediction
http: / / npsa ̄pbil.ibcp. fr / cgi ̄bin / npsa_automat. pl? page =
npsa_sopma.html
SWISS ̄MODEL 蛋白质三级结构预测
Tertiary structure prediction
http: / / swissmodel.expasy.org / interactive
TMHMM Server V.2.0 蛋白质跨膜结构域预测
Transmembrane structure prediction
http: / / www.cbs.dtu.dk / services / TMHMM ̄2.0 /
NetPhos 2.0 Server 蛋白质磷酸化位点预测
Phosphorylation sites prediction
http: / / www.cbs.dtu.dk / services / NetPhos /
PSORT 蛋白质亚细胞定位
Subcellular localization
http: / / psort.hgc.jp / form.html
SignalP 4.1 Server 蛋白质信号肽预测
Signal peptide prediction
http: / / www.cbs.dtu.dk / services / SignalP /
TargetP 1.1 Server 蛋白质导肽预测
Leading peptide prediction
http: / / www.cbs.dtu.dk / services / TargetP /
Conserved Domains
Database
蛋白质保守结构域分析
Conserved domain analysis
http: / / www.ncbi.nlm.nih.gov / Structure / cdd / wrpsb.cgi
Blastn / Blastp 核酸和蛋白质比对
Nucleic acid and protein alignment
http: / / blast.ncbi.nlm.nih.gov / Blast.cgi
MEGA 6.06 系统进化树构建
Phylogenetic tree construction
http: / / www.megasoftware.net /
Center for Biotechnology Information)数据库中已有的
油棕等热带植物 1 ̄脱氧木酮糖 ̄5 ̄磷酸合成酶基因
的核酸序列及其对应的蛋白质序列ꎬ如油棕(Elaeis
guineensis Jacq.ꎬAY583783、AAS99588)、白木香(Aq ̄
uilaria sinensisꎬ JX860326、 AFU75321 )、 杜 仲
(Eucommia ulmoidesꎬJX458817、AFU93068)、巴西橡
胶(Hevea brasiliensisꎬAB294698、BAF98288)、麻疯树
( Jatropha curcasꎬ KF955327、 AHL17583 )、 可 可
(Theobroma cacaoꎬXM007035370、XP007035432)ꎻ另
外ꎬ构建系统进化树还用到了以下几种植物的 DXS
蛋白质序列:高良姜(Alpinia officinarumꎬAEK69518)、
阳春砂(Amomum villosum ꎬACR02668)、拟南芥(Ara ̄
bidopsis thalianaꎬNP193291)、巴豆(Croton stellatopilo ̄
susꎬ BAF75640 )、 海 岛 棉 ( Gossypium barbadenseꎬ
ABN13970 )、 苜 蓿 ( Medicago truncatulaꎬ
XP003609567 )、 烟 草 ( Nicotiana tabacumꎬ
ACF60511)、粳稻(Oryza sativaꎬNP001055524)、葛麻
姆 ( Pueraria montanaꎬ AAQ84169 )、蓖麻 ( Ricinus
communisꎬ XP002516843 )、 罗 汉 果 ( Siraitia
grosvenoriiꎬAEM42997)、番茄(Solanum lycopersicumꎬ
NP001234672)、马铃薯(S. tuberosumꎬNP001275130)、
南非醉茄(Withania somniferaꎬAFI98878)、玉米(Zea
maysꎬNP001157805)ꎮ
1.2 方法
所用生物信息学工具和软件如表 1 所示ꎬ具体
参数设置见结果ꎮ
2 结果与分析
2.1 油棕等热带植物 DXS基因核酸及氨基酸序列分析
利用 ORF Finder 和 ProtParam 对油棕、巴西橡
胶、麻疯树、杜仲、可可和白木香 6 种热带植物 DXS
基因的核酸序列及其氨基酸序列进行分析ꎮ 分析结
果如表 2所示ꎬ不同植物 DXS 基因开放读码框的碱
基数及所编码的氨基酸残基数变化不大ꎬ起始密码
子均为 ATGꎬ终止密码子分为 3 种ꎬ即油棕、巴西橡
胶和白木香均为 TAGꎬ麻疯树和杜仲为 TAAꎬ可可
则为 TGA ꎻ而且不同植物 DXS 的分子量、理论等电
点和高含量氨基酸成分基本一致ꎻ Leu、Ala、Gly 是
这 6种热带植物 DXS中含量最丰富的氨基酸ꎻ根据
3744期 石鹏等: 油棕等热带植物 DXS基因的生物信息学分析
表 2 6种热带植物 DXS核酸及氨基酸序列分析
Table 2 Analysis of DXS nucleic acid and amino acid from 6 tropical plants
来源
Resource
起始
密码子
Initiation
codon
终止
密码子
Termination
codon
aa残基数
Number of
amino
acid
residue
相对分子
量 / KD
Relative
molecular
weight /
kilodalton
理论等
电点
Theoretical
isoelectric
point
高含量氨基酸
Amino acid of
high content
亲水性
总平均值
Grand
average of
hydropathicity
不稳定
指数
Instability
index
油棕
Elaeis guineensis Jacq.
ATG TAG 707 76.44 7.00 Leu(10.5) Ala(9.6) Gly(9.6) -0.120 40.42
巴西橡胶
Hevea brasiliensis
ATG TAG 720 77.44 7.34 Leu(10.3) Gly(10.0) Ala(9.7) -0.063 39.55
麻疯树
Jatropha curcas
ATG TAA 720 77.49 7.12 Leu(10.3) Ala(10.1) Gly(9.9) -0.096 41.76
杜仲
Eucommia ulmoides
ATG TAA 712 76.96 8.14 Leu(9.7) Gly(9.7) Ala(9.4) -0.076 39.54
可可
Theobroma cacao
ATG TGA 722 77.59 6.94 Leu(10.4) Ala(10.4) Gly(9.8) -0.069 39.75
白木香
Aquilaria sinensis
ATG TAG 717 76.91 6.65 Ala(10.9) Leu(9.9) Gly(9.6) -0.113 41.39
不稳定性指数小于 40为稳定性蛋白质判断ꎬ巴西橡
胶、杜仲和可可的 DXS 为稳定类蛋白质ꎬ而油棕、麻
疯树和白木香的 DXS 为不稳定类蛋白质ꎻ且它们的
亲水系数较低ꎬ均为亲水性蛋白质ꎮ
2.2 油棕等热带植物 DXS 亲水性 /疏水性的预测和
分析
利用 ProtScale对油棕氨基酸序列的疏水性 /亲
水性进行预测(图 1)ꎮ 从图 1发现ꎬ多肽链第 55 位
的 His 具有最低的分值-2.656ꎬ亲水性最强ꎻ而第
185位的 Gly 具有最高的分值 2.067ꎬ疏水性最强ꎮ
整体来看ꎬ亲水氨基酸分布比较均匀ꎬ且数量大于疏
水氨基酸ꎬ故推测油棕 DXS蛋白质为亲水性蛋白ꎮ
2.3 油棕等热带植物 DXS二级结构的预测和分析
利用 SOPMA对油棕 DXS氨基酸序列进行二级
结构预测(图 2)ꎮ 从图 2 可以发现ꎬ油棕 DXS 由
40.45%的无规则卷曲(Random coil)ꎬ38.19%的 α ̄
螺旋 ( Alpha helix )ꎬ 15. 28% 的延伸链 ( Extended
strand)和6.08%β ̄转角(Beta turn)组成ꎮ 在整个蛋
白质结构中ꎬ无规则卷曲和 α ̄螺旋是油棕 DXS最主
要的结构元件ꎬ而延伸链和 β ̄转角则次之ꎮ
2.4 油棕等热带植物 DXS蛋白三级结构的预测和分析
利用 SWISS ̄MODEL的同源建模方法预测油棕
DXS蛋白的三级结构ꎬ获得油棕 DXS 三维结构模式
(图 3)ꎮ 同源建模分析发现其与蛋白质模板
2o1s.1.A有 50.58%的相似性ꎬ且根据三级结构特征
推测为 1 ̄脱氧木酮糖 ̄5 ̄磷酸合成酶(1 ̄deoxy ̄D ̄xy ̄
lulose ̄5 ̄phosphate synthase)ꎮ 从图 3 可以看出ꎬ紫
色卷曲状结构代表的无规则卷曲遍布在油棕 DXS
三级结构中ꎮ
2.5 油棕等热带植物 DXS跨膜结构域的预测和分析
利用 TMHMM对不同热带植物 DXS 氨基酸序
列的跨膜结构域进行预测(图 4)ꎬ未发现明显的跨
膜结构域ꎮ
2.6 油棕等热带植物 DXS翻译后修饰的预测和分析
利用 NetPhos 2.0 Server 对油棕 DXS 的翻译后
修饰预测的结果表明(图 5)ꎬ整个多肽链中分值在
0.5(阈值)以上的氨基酸位点有 36个ꎬ可知油棕 DXS
的磷酸化位点有 36个ꎬ其中 Ser:17ꎬThr:11ꎬTyr:8ꎮ
2.7 油棕等热带植物 DXS的亚细胞定位分析
利用 PSORT 对油棕等热带植物的 DXS 基因的
编码产物进行亚细胞定位分析ꎬ结果表明ꎬ油棕和其
它热带植物的 DXS 蛋白质定位在叶绿体中基质的
可能性最大ꎬ该结果说明热带植物 DXS 蛋白可能主
要是在叶绿体基质中行使功能ꎮ
2.8 油棕等热带植物 DXS信号肽和导肽的预测及分析
利用 SignalP 4.1 Server 对油棕 DXS 氨基酸序
列的信号肽进行预测(图 6)ꎮ 图 6 结果表明ꎬ油棕
DXS不含有信号肽ꎮ
同时利用 TargetP1.1 Server 进行导肽分析也证
实了其作为叶绿体转运蛋白(分值为 0.409)的可能
性比较大(图 7)ꎬ该蛋白存在信号肽(分值为0.030)
的可能性很小ꎮ 此结果与亚细胞定位及信号肽预测
结果基本一致ꎮ
2.9 油棕等热带植物 DXS保守结构域的预测和分析
通过 NCBI中的保守结构域数据库 (Conserved
DomainDatabase ꎬCDD)工具对DXS基因的氨基酸
474 广 西 植 物 36卷
图 1 油棕 DXS的疏水性 /亲水性预测和分析
Fig. 1 Hydrophobicity or hydrophilicity prediction
and analysis of EgDXS
序列分析发现ꎬDXS 含有硫胺素焦磷酸盐 ( TPPꎬ
Thiamine pyrophosphate)和嘧啶 ( PYRꎬPyrimidine)
结合位点ꎬTPP _DXS、TPP _PYR_DXS_TK_like 和
Transketolase_C保守域ꎬTPP_enzymes、TPP_enzyme_
PYR 和 Transketolase _ C superfamily 家 族 和
PLN02582保守结构域ꎮ 其中 PLN02582 是 1 ̄脱氧
木酮糖 ̄5 ̄磷酸合成酶的特征序列(图 8)ꎮ
2.10 油棕等热带植物 DXS序列比对及系统进化分析
用 Blast程序比对油棕与其它植物 DXS 的核酸
及氨基酸序列同源性ꎬ结果表明油棕 DXS 的核酸序
列与其它植物如高良姜 ( HQ874656 )、麻疯树
( KF955327 )、 葛 麻 姆 ( FJ455512 )、 橡 胶
(AB294698)、番茄 ( NM001247743 ) 和南非醉茄
(JQ710678)的一致性分别是 80%、80%、79%、79%、
图 2 油棕 DXS二级结构预测 蓝色线条ꎬ红色线条ꎬ绿色线条和紫色线条分别表示 α ̄螺旋ꎬ延伸链ꎬβ ̄转角和无规则卷曲ꎮ 横坐标
和纵坐标分别表示二级结构在所处的位置和数量ꎮ
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.
79%和 79%ꎻ蛋白质的 Blastp 比对发现ꎬ油棕 DXS
的氨基酸序列与其他植物的氨基酸序列有更高的相
似 性ꎬ 如 与 长 春 花 ( AGL40532 )、 阳 春 砂
( ACR02668 )、 玉 米 ( ACG27905 )、 麻 疯 树
( AHL17583 )、 高 良 姜 ( AEK69518 ) 和 罗 汉 果
(AEM42997) 的一致性分别为 86%、 85%、 86%、
86%、86%和 85%ꎮ 以上结果表明ꎬ植物 DXS 氨基
酸序列之间的相似性高于核酸序列ꎮ
把 21 个来自不同植物的 DXS 氨基酸序列用
MEGA 6.06软件采用 Neighbor ̄Joining法ꎬ进行1 000
次 bootstrap统计学检验构建系统进化树(图 9)ꎮ 图
9结果表明ꎬ21种不同植物的 DXS 氨基酸序列可以
分为 3大类ꎬ巴西橡胶和蓖麻等 10种植物分为第Ⅰ
类ꎬ其中大戟科植物麻疯树、巴豆、蓖麻和巴西橡胶
在此类ꎬ而白木香、可可和杜仲等热带植物也聚在此
类ꎻ豆科植物和茄科植物分在第Ⅱ类ꎬ其中同为豆科
植物的葛麻姆和苜蓿聚成一簇的信心指数为
100%ꎬ马铃薯、番茄、南非醉茄和烟草为茄科ꎬ其中
马铃薯和番茄聚成一簇的信心指数为 100%ꎻ油棕、
姜科和禾本科植物分在第Ⅲ类ꎬ其中同为禾本科的
粳稻和玉米聚成一簇的信心指数为 100%ꎬ同为姜
科的高良姜和阳春砂聚成一簇的信心指数为
100%ꎮ 从聚类分析结果可以发现ꎬ用 DXS 氨基酸
序列构建的系统进化树可以比较真实地反映不同植
5744期 石鹏等: 油棕等热带植物 DXS基因的生物信息学分析
图 3 油棕 DXS三级结构预测
Fig. 3 Tertiary structure prediction of EgDXS
图 4 油棕 DXS的跨膜结构域预测
Fig. 4 Transmembrane structure prediction of EgDXS
图 5 油棕 DXS翻译后修饰预测 磷酸化位点
分为丝氨酸ꎬ苏氨酸和酪氨酸 3种ꎮ
Fig. 5 Posttranslational modification prediction of EgDXS
Phosphorylation sites are divided into serineꎬ threonine and tyrosine.
物之间的进化关系ꎬ对判断植物间的亲缘关系具有
很好的参考价值ꎮ
图 6 油棕 DXS 信号肽预测 C ̄score. 切割位点的可能
性ꎻ S ̄score. 氨基酸是信号肽的可能性ꎻ Y ̄score. 结合 C ̄score
和 S ̄score的导数ꎬ能更好地反映切割位点的可能性ꎮ
Fig. 6 Signal peptide prediction of EgDXS C ̄score. Possi ̄
bility 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导肽预测 Len. 蛋白质长度ꎻ cTP. 作为
叶绿体转运蛋白的可能性ꎻ mTP. 作为线粒体转运蛋白的可能
性ꎻ 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.
3 讨论与结论
类胡萝卜素对于植物来说是必需的ꎬ并且在为
人类提供营养和保护健康方面具有重要作用ꎮ 本文
利用生物信息学工具和软件对油棕等热带植物的
DXS基因的核酸及蛋白质理化性质、结构组成和功
能特征等进行了预测和分析ꎬ并构建了油棕等 21 种
植物的系统进化树ꎮ 预测结果表明ꎬ不同植物 DXS
的理化性质差异不大ꎬ蛋白质稳定性不一ꎬ都属于亲
水性蛋白质ꎮ 跨膜结构域、信号肽和导肽预测结果
表明ꎬ油棕等热带植物的 DXS 都没有发现明显的跨
膜结构域和信号肽ꎬ主要是在叶绿体基质中发挥作
用ꎬ可能作为叶绿体转运蛋白行使功能ꎮ DXS 的保
守结构域分析发现ꎬDXS 含有硫胺素焦磷酸盐和嘧
啶结合位点ꎬ TPP_DXS、TPP_PYR_DXS_TK_like和
674 广 西 植 物 36卷
图 8 油棕 DXS保守结构域预测
Fig. 8 Conserved domain prediction of EgDXS
图 9 油棕等 21种植物 DXS氨基酸序列系统进化树分析
Fig. 9 Phylogenetic tree analysis of Amino acid series from 21 plants including oil palm
Transketolase_C保守域ꎬTPP_enzymes、TPP_enzyme_
PYR 和 Transketolase _ C superfamily 家 族 以 及
PLN02582保守结构域ꎬ其中 PLN02582 保守结构域
是 1 ̄脱氧木酮糖 ̄5 ̄磷酸合成酶的特征结构域ꎮ
在 MEP 途径中ꎬ以丙酮酸和 3 ̄磷酸甘油醛为底
物ꎬ先后经过 8个酶的催化最终形成了合成类异戊
二烯的通用底物 IPP 和 DMAPPꎮ 实验证实作为
MEP 途径第一个酶 DXS 在 IPP 和 DMAPP 合成过
程中发挥重要作用ꎮ 拟南芥中对 MEP 途径中的酶
研究较多ꎬ目前 MEP 途径中的所有酶都已经确认ꎮ
MEP 途径的第一个反应就是在 DXP 合成酶 DXS 的
催化作用下甘油醛 ̄3 ̄磷酸硫胺素缩合生成 1 ̄脱氧
木酮糖 ̄5 ̄磷酸ꎮ 研究表明ꎬ通过 MEP 途径合成类
胡萝卜素主要受到限速酶 DXSꎬDXR 和 HDR 的控
制( Sawitri & Wallieꎬ 2005ꎻ Elizabeth et alꎬ 2009)ꎮ
Paruicia et al(2006)认为在充分了解植物类胡萝卜
素合成机制的基础上ꎬ可以通过生物技术手段提高
植物中的类胡萝卜素ꎬ为人类提供类胡萝卜素含量
更高的食物ꎮ 目前ꎬ番茄、马铃薯、油菜和水稻中都
已经开展了相关转基因实验ꎮ 植物中类胡萝卜素的
代谢和生物合成途径已经得到广泛研究ꎬ合成途径
中几乎所有编码合成酶的基因都被分离出来ꎮ
Sanae & Akemi(2006)研究菊花花瓣颜色发现ꎬ是类
胡萝卜素的合成和积累机制不同形成了黄色和白色
两种花瓣颜色ꎬ其中 LCYEꎬLCYB 和 DXS 的表达量
与花瓣中类胡萝卜素的含量息息相关ꎮ 在番茄果实
发育过程中ꎬDXS的表达量与果实中类胡萝卜素含
量呈正相关关系ꎬ说明 DXS在番茄果实发育和器官
形成过程中发挥重要作用(Luisa et alꎬ2000)ꎮ 莱茵
衣藻中研究发现ꎬDXS和 DXR等基因的表达量随着
光周期波动ꎬDXS等基因的表达量表现出典型的昼
夜节律模式ꎬ并且也会受到中间代谢产物的反馈调
7744期 石鹏等: 油棕等热带植物 DXS基因的生物信息学分析
节(Tian et alꎬ2010)ꎮ 因此ꎬ在植物类胡萝卜素合成
代谢过程中ꎬDXS 是调控类胡萝卜素合成的关键
酶ꎬ对其进行生物信息学分析为今后深入研究具有
重要意义ꎮ 油棕和巴西橡胶等热带植物 DXS 生物
信息学分析ꎬ可以寻找这些不同科的热带植物在
DXS上的共同特征和不同点ꎬ希望为热带植物 DXS
结构和功能研究奠定基础ꎮ 植物激素、类胡萝卜素
和甾醇等类异戊二烯物质在热带植物生命活动中具
有重要作用ꎬ而它们在热带植物中的具体功能作用
及其控制基因研究较少ꎬ需要进一步开展相关研究ꎬ
揭开类异戊二烯在热带植物中的生物合成途径及其
调控机制ꎮ 油棕等 21种植物 DXS 氨基酸序列比对
和系统进化树分析发现ꎬ不同植物的 DXS 序列相似
性较高ꎬ功能域非常保守ꎬ以 DXS 氨基酸序列构建
的进化树能准确反映不同植物间的亲缘关系ꎬ也可
作为评价油棕等热带植物种质资源的重要依据ꎮ
参考文献:
AGUILA RSꎬ MANUEL RCꎬ 2012. Carotenoid biosynthesis in Ara ̄
bidopsis: 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 Bi ̄
olꎬ 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 of Ginkgo biloba DXS gene enco ̄
ding 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 in Catharanthus
roseu [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 in Rhodobacter capsulatus
[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 of in vivo plastoquinone 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 from Streptomyces sp. 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]. 北京: 中国农业出版社: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 from Escherichia coli encoding a tran ̄
sketolase ̄like enzyme that catalyzes the synthesis of D ̄1 ̄deoxyx ̄
ylulose ̄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: regula ̄
tory 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 di ̄
versity [J]. Scienceꎬ 277:1 788-1 789.
SANAE Kꎬ AKEMI Oꎬ 2006. Regulation of carotenoid biosynthesis
in petals and leaves of chrysanthemum (Chrysanthemum morifo ̄
lium) [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 palm Elaeis guineensis Jacq. [J]. Plant Sciꎬ169:571-578.
SHAN Lꎬ LI Lꎬ 2008. Carotenoid metabolism: biosynthesisꎬ regula ̄
tionꎬ and beyond [J]. J Integr Plant Biolꎬ 50(7):778-785.
TIAN HSꎬ CHENG QLꎬ YUAN YHꎬ et alꎬ 2010. Coordinated reg ̄
ulation of gene expression for carotenoid metabolism in
Chlamydomonas reinhardtii [J]. J Integr Plant Biolꎬ 52(10):
868-878.
VINCENT Cꎬ DIRK JLꎬ 2011. From glycosylation disorders to doli ̄
chol 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. Isopre ̄
noid 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 char ̄
acterization 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.
874 广 西 植 物 36卷