免费文献传递   相关文献

Cloning and analysis of the cytosolic glyceraldehyde 3-phosphate dehydrogenase (GAPC) gene from Solanum tuberosum and Arabidopsis thaliana

马铃薯和拟南芥GAPC酶基因的克隆及分析



全 文 :书马铃薯和拟南芥犌犃犘犆酶基因的克隆及分析
邹雪,张烨,吴明阳,王西瑶
(四川农业大学农学院,四川 成都611130)
摘要:胞质三磷酸甘油醛脱氢酶(GAPC)是糖酵解中的关键酶,近年的研究表明其对盐、低温、高温、氧化、高渗、低
磷等多种逆境胁迫均有响应。本实验采用RT-PCR,分别从马铃薯和拟南芥中克隆GAPC的编码区序列(CDS),
并作生物信息学比较,同时构建2个基因的植物表达载体。结果表明,犛狋犌犃犘犆和犃狋犌犃犘犆2的CDS长度均为
1017bp,相似性84%,编码338个氨基酸,相似性为92%。蛋白分子量分别为36.65和36.91kDa,理论等电点为
6.34和6.67,均属稳定蛋白;三级结构预测表明两者和水稻(2e5rC)GAPC蛋白结构相似,具有 GAPC保守功能
域。多种植物犌犃犘犆基因的序列进化分析表明同科属植物的同源性较高可归于同一分支,与现有的分类系统相对
应。以pBI121为基础载体,构建了由CaMV35S启动犌犃犘犆基因的植物表达载体,导入农杆菌EHA105,经PCR
检测确定获得阳性克隆。本研究为获得转犌犃犘犆基因的材料并进一步研究GAPC在逆境中的功能奠定基础。
关键词:马铃薯;拟南芥;胞质3磷酸甘油醛脱氢酶;基因克隆;生物信息学分析
中图分类号:S532.03;Q943.2  文献标识码:A  文章编号:10045759(2014)01023909
犇犗犐:10.11686/cyxb20140129  
  3磷酸甘油醛脱氢酶(glyceraldehyde3phosphatedehydrogenase,GAPDH)是生物中普遍存在的酶。在高
等植物中,根据其存在位置、催化的反应可分为3类:1)磷酸化NAD依赖型和非磷酸化NADP依赖型,由4个相
同GapC亚基组成,存在于胞质中参与糖酵解,前者催化的反应产生含高能磷酸键的产物并在下一步反应中生成
3磷酸甘油酸和ATP,后者催化反应直接形成3磷酸甘油酸,没有ATP的产生;2)磷酸化NAD依赖型,存在于
质体中;3)磷酸化NADPH依赖型,由GapA和GapB两种亚基组成,存在于叶绿体中,参与卡尔文循环[12]。由
于糖酵解和卡尔文循环都是植物细胞能量代谢的主要途径,因此GAPDH在植物碳代谢和能量代谢中起关键作
用。Brinkmann等[3]、Russel 和Sachs[4]最初从玉米(犣犲犪犿犪狔狊)中克隆到3个胞质三磷酸甘油醛脱氢酶基因
(犌犃犘犆)[34],迄今已从多种植物中克隆到该基因,如拟南芥(犃狉犪犫犻犱狅狆狊犻狊狋犺犪犾犻犪狀犪)、水稻(犗狉狔狕犪狊犪狋犻狏犪)、西伯利
亚蓼(犘狅犾狔犵狅狀狌犿狊犻犫犻狉犻犮狌犿)[57]。最初人们认为GAPC蛋白仅参与糖酵解其mRNA在细胞内表达稳定,因此将
犌犃犘犆作为RTPCR(reversetranscriptionpolymerasechainreaction,反转录聚合酶链式反应)等检测中的内参
基因,但近年来越来越多的研究表明植物GAPC对盐、干旱、高渗、低温、高温、活性氧、低磷、病菌侵染等多种逆
境胁迫均有响应,并通过过量表达犌犃犘犆基因获得了抗盐胁迫的水稻和马铃薯(犛狅犾犪狀狌犿狋狌犫犲狉狅狊狌犿)[814]。
磷对重要粮食作物马铃薯的块茎形成以及淀粉积累都有良好的促进作用,并交叉影响抗逆表现,而我国农田
有2/3缺磷,因此研究马铃薯的磷吸收利用机制并通过生物技术增强马铃薯对磷的吸收和利用效率具有重要意
义[1518]。王西瑶等[10,19]从拟南芥耐低磷突变体中首先鉴定到犃狋GAPC2蛋白与耐低磷性状相关,进一步分析表
明低磷诱导该基因(犃狋犌犃犘犆2)转录增强,TDNA插入突变体表现出对低磷胁迫的敏感性。Penaloza等[9]研究
发现白羽扇豆(犔狌狆犻狀狌狊犪犾犫狌狊)在低磷胁迫的早期及晚期与NAD依赖型GAPDH同源的基因表达增强,推测糖
酵解在白羽扇豆响应低磷胁迫的特殊机制中起作用。但 Hajirezaei等[20]通过反义技术沉默马铃薯犌犃犘犆基因
(犛狋犌犃犘犆),表明该基因在马铃薯块茎形成的代谢调节上仅发挥次要作用。
犛狋犌犃犘犆和犃狋犌犃犘犆2所编码的蛋白在功能多样性上是否具有差异,犌犃犘犆的过量表达能否提高马铃薯对
第23卷 第1期
Vol.23,No.1
草 业 学 报
ACTAPRATACULTURAESINICA   
239-247
2014年2月
收稿日期:20130226;改回日期:20130402
作者简介:邹雪(1984),女,四川彭山人,在读博士。Email:zou_xue_2008@aliyun.com
通讯作者。Email:wxyrtl@163.com
低磷、盐等多种逆境胁迫的抗性,为此本实验分别克隆了这2个基因并对两者的序列和蛋白结构作对比分析;同
时利用多种植物的犌犃犘犆序列建立进化树,分析环境胁迫对GAPC进化及功能多样性产生的影响;最后分别构
建了这2个基因的植物表达载体,为进一步验证GAPC的功能多样性并获得耐逆境胁迫的马铃薯材料奠定基
础。
1 材料与方法
1.1 试验材料
野生型拟南芥种子(Columbia生态型),由四川农业大学玉米研究所张素芝教授提供。马铃薯栽培型品种米
拉,由四川农业大学农学院马铃薯研究开发中心提供。植物表达载体pBI121由四川农业大学李立芹副教授惠
赠,克隆载体pMDTM19TSimpleVector(TaKaRa)购自宝生物公司。
1.2 方法
1.2.1 犃狋犌犃犘犆2和犛狋犌犃犘犆基因的克隆 实验于2010年10-12月在四川农业大学农学院马铃薯研究开发
中心完成,分别取拟南芥和马铃薯幼苗叶片100mg,采用Trizol法提取总RNA,以40μLDEPC(diethylpyro
carbonate,焦碳酸二乙酯)处理水溶解。反转录合成cDNA第一链,体系为20μL,在用DEPC处理过的无菌0.2
mLPCR管中依次加入:oligodT1μL,RNA3μL,DEPC处理水8μL,轻微混匀并离心30s,65℃变性5min,冰
上急冷2min,按序加入5×ReactionBuffer4μL,RNase抑制剂 (20U/μL)1μL,10mmol/LdNTPMix2μL,
反转录酶 (200U/μL)1μL。轻微混匀稍离心,42℃,60min,72℃10min终止反应。以第一链为模板,PCR扩
增目的片段,根据NCBI上公布的犃狋犌犃犘犆2(拟南芥)和犛狋犌犃犘犆(马铃薯)的CDS全长序列(AccessionNo.
AY090275,AF527779),用Oligo6软件分别设计引物并引入犡犫犪Ⅰ、犛犿犪Ⅰ酶切位点。
犃狋犌犃犘犆2P1:5′TCTAGAATGGCTGACAAGAAGATCAGAAT3′,P2:5′CCCGGGTTAGGCCTTTGA
CATGTGAACG3′。
犛狋犌犃犘犆P1:5′TCTAGAATGGCCAATGGCAAGATCAAAAT3′,P2:5′CCCGGGTCAAGCCTTGGC
CATAT3′。
PCR扩增体系为25μL,反应体系中加入:ddH2O16μL,10×PCRbuffer2.5μL,P1(10μmol/L)0.5μL,
P2(10μmol/L)0.5μL,dNTPMixture(2.5mmol/L)2.0μL,cDNA,1.0μL(拟南芥或马铃薯),犜犪犓犪犚犪犈狓
犜犪狇0.5μL,MgCl22.0μL。PCR扩增反应条件,拟南芥:94℃预变性5min;94℃变性30s,58℃退火30s,72℃
延长60s,30个循环;72℃ 延长10min。马铃薯:95℃预变性3min;95℃变性30s,68℃退火延长90s,30个循
环;68℃ 延长10min。产物于1%琼脂糖凝胶电泳检测。
1.2.2 基因克隆及测序 电泳检测后回收纯化目的片段,与克隆载体pMDTM19TSimpleVector在4℃连接24
h,10μL连接体系如下:pMD19TVector0.5μL,InsertDNA3.5μL,SolutionI5μL,ddH2O1μL。将连接产
物转化大肠杆菌DH5α感受态,在涂有Xgal(20mg/mL)40μL 和IPTG(200mg/mL)4μL的含 Amp(100
μg/mL)的LB平板上作蓝白斑筛选。将筛选的白色单克隆提质粒作PCR检测和质粒酶切鉴定,酶切反应体系
20μL:犡犫犪Ⅰ和犛犿犪Ⅰ0.5/0.5μL,0.1%BSA0.2μL,10×TBuffer2.0μL,ddH2O1.8μL,质粒15.0μL。
30℃水浴6h,1%琼脂糖凝胶电泳检测。最后各随机选取3个克隆,3次重复测序(上海生工生物工程股份有限
公司)。
1.2.3 序列分析与比较 将测得的序列拼接后通过NCBI的BLAST寻找核苷酸同源序列(http://blast.nc
bi.nlm.nih.gov/)。利用 ProtParam 程序预测氨基酸序列分子量和理论等电点(pI)、不稳定系数等,用
ProScale程序预测蛋白疏水性,通过ConservedDomains程序进行保守功能区分析(http://au.expasy.org/)。
将氨基酸序列提交SwissModel预测三级结构,用Rasmol2.6查看图像。用ClustalX和DNAMAN对克隆的
两种基因作核苷酸和氨基酸序列比对,并将NCBI中不同植物的犌犃犘犆序列数据作进化树分析。
1.2.4 植物表达载体的构建 用内切酶犡犫犪I和犛犿犪I对表达载体pBI121和克隆载体pMDAtGAPC2、pMD
042 ACTAPRATACULTURAESINICA(2014) Vol.23,No.1
StGAPC进行酶切,酶切体系50μL:10×TBuffer5.0μL,0.1%BSA5.0μL,犡犫犪I和犛犿犪I1.0μL,质粒30.0
μL,ddH2O8.0μL,30℃进行4~6h。电泳后切胶回收相应片段连接,连接体系20.0μL:T4DNA连接酶(3
U/μL)2.0μL,10×Buffer2.0μL,pBI121Vector4.0μL,目的片段12.0μL,4℃24h。连接产物转化大肠杆菌
DH5α感受态,在含卡那霉素(Kan,50μg/mL)的平板上筛选阳性克隆,经质粒PCR、酶切鉴定后将正确的重组
质粒用冻融法转化到根癌农杆菌EHA105的感受态中,在含Kan+(50μg/mL)和利福平(25μg/mL)YEB平板
上挑选单菌落扩增培养后提质粒作PCR鉴定。
2 结果与分析
2.1 犃狋犌犃犘犆2和犛狋犌犃犘犆基因的克隆
从野生型拟南芥Col和栽培型马铃薯品种米拉中分别提取RNA,经凝胶电泳检测,如图1A所示,有28S,
18S以及5S三个条带,说明RNA的完整性较好。由Oligod(T)作引物转录合成cDNA第一链,以此为模板,根
据Genbank公布的序列设计引物,用高保真DNA聚合酶犜犪犓犪犚犪犈狓犜犪狇扩增目的条带,取2μL扩增产物进
行电泳检测,表明得到了预期大小的片段(约1000bp,图1B)。回收目的片段,与克隆载体pMDTM19TSimple
Vector连接后转化大肠杆菌DH5α,挑选白斑进行滞后质粒筛选(图2A),进行PCR及犡犫犪Ⅰ、犛犿犪Ⅰ双酶切检
测(图2B、图2C)。检测结果表明无论是扩增还是酶切均能获得1000bp左右片段。
图1 总犚犖犃提取及犚犜-犘犆犚扩增犌犃犘犆基因
犉犻犵.1 犜狅狋犪犾犚犖犃犲狓狋狉犪犮狋犻狅狀犪狀犱犌犃犘犆犵犲狀犲犪犿狆犾犻犳犻犮犪狋犻狅狀
 A:总RNA提取 RNAextraction;1~2:拟南芥犃狉犪犫犻犱狅狆狊犻狊;3~6:马铃薯Potato.B:犌犃犘犆基因的RT-PCR扩增产物 Amplificationprod
uctsof犌犃犘犆genebyRT-PCR;M:MarkerⅢ;1~2:拟南芥犃狉犪犫犻犱狅狆狊犻狊;3~4:马铃薯Potato. 
图2 含目的基因的克隆检测
犉犻犵.2 犐犱犲狀狋犻犳犻犮犪狋犻狅狀狅犳犮犾狅狀犲狊
  A:滞后质粒筛选Identificationofdelayedplasmid;B:目的片段PCR扩增
犌犃犘犆genePCRamplification;M:DNA MarkerⅢ;1:蓝斑质粒 Bluecolor
plasmid;2~4:pMDStGAPC;5~7:pMDAtGAPC2.C:犡犫犪I和犛犿犪I酶
切质粒Plasmiddigestionby犡犫犪Iand犛犿犪I;M:DNA MarkerⅢ;1~2:
pMDStGAPC;3~4:pMDAtGAPC2.
142第23卷第1期 草业学报2014年
  挑选重组子pMDAtGAPC2和pMDStGAPC各3个,每个3次重复作测序,其中pMDAtGAPC2有2个
重组子的序列与公布序列完全一致,另有1个发生A到G的变异。而pMDStGAPC的3个重组子与公布序列
均有差异:其中1个在404bp处由T变为C;1个在143bp处由T变为C,且871bp处由G变为A;第3个在
976bp处由T变为 A,由于3个重组子与公布序列的差异在不同位置,所以很可能是PCR扩增时碱基错配造
成,另选3个pMDStGAPC重组子测序,其中2个与公布序列完全一致,另1个在158bp处C变为T,601处A
变为G。将所获得的测序正确克隆扩繁后再作测序,仍与公布的序列完全相同,至此获得了与公布序列完全一致
的犃狋犌犃犘犆2和犛狋犌犃犘犆基因的克隆。
2.2 犛狋犌犃犘犆和犃狋犌犃犘犆2序列分析
2.2.1 犛狋犌犃犘犆和犃狋犌犃犘犆2氨基酸序列分析和结构域预测 测序结果表明,两者的开放阅读框均为1017bp,
编码338个氨基酸和1个终止密码子。两者的比对表明核苷酸序列相似性为84%,氨基酸序列相似性达到
92%。犛狋GAPC氨基酸序列分子量为36.646kDa,理论pI值6.34,整个氨基酸组成中缬氨酸(Val)含量最高达
10.4%,其次是丙氨酸(Ala)占9.8%,蛋白不稳定系数为20.21,属稳定蛋白;犃狋GAPC2氨基酸序列分子量为
36.913kDa,理论pI值6.67,整个氨基酸组成中缬氨酸(Val)含量最高达11.5%,其次是赖氨酸(Lys)占9.2%,
蛋白不稳定系数为22.6,属稳定蛋白。以上可以看出两者的基本性质相似。从图3A可以看出两者在氨基酸序
列上共有26处存在差异,有些是同类氨基酸间的转变如:KR、AV、ST、DE,有些则是非极性与极性间的转变
如:GK、AP、MT、AE。
图3 犛狋犌犃犘犆和犃狋犌犃犘犆2氨基酸序列与单体三维结构比较
犉犻犵.3 犃犾犻犵犿犲狀狋狅犳犪犿犻狀狅犪犮犻犱狊犲狇狌犲狀犮犲犪狀犱犿狅狀狅犿犲狉3犇狊狋狉狌犮狋狌狉犲狅犳犛狋犌犃犘犆犪狀犱犃狋犌犃犘犆2
 A:氨基酸序列比对,相同区域以阴影表示,下划线标记的8个氨基酸为GAPDH活性位点的特征序列Aligmentofaminoacidsequence,thesha
dedpartshowedidenticalsequences,theunderlineshowedtotalyconservedcharacteristicsequenceinactivesiteofGAPDH;B:单体三维结构 Mono
mer3Dstructure;1,折叠差异Differenceofstrands;2~3,α螺旋差异 Differenceofαhelices.
242 ACTAPRATACULTURAESINICA(2014) Vol.23,No.1
  利用NCBI网站的ConservedDomains程序作保守功能区分析表明犛狋GAPC和犃狋GAPC2具有GAPDH典
型的2个保守区域:NADB_Rossmann(登录号cl09931),即 N 端 NAD(P)结合域;Gp_dh_C (登录号
pfam02800),即C端行使糖结合功能和催化功能的区域。两者均含所有GAPDH催化活性位点固有的特征序列
ASCTTNCL(两多肽链上的154~161位)(图3A中下划线所示)。
2.2.2 犛狋GAPC和犃狋GAPC2蛋白结构预测 利用SwissModel在线分析功能预测犛狋GAPC和犃狋GAPC2蛋
白的三级结构,两者均与模板3e5rC(犗狉狔狕犪狊犪狋犻狏犪 GAPDH Cchain)的相似度最高,分别为85.629%和
85.586%,并得到了犛狋GAPC和犃狋GAPC2蛋白单体的3D预测结构(如图3B所示),三者的3D结构非常相似,
其中左边区域有一个由7个β折叠构成的马鞍型超二级结构位于催化结构域中,起着固定活性中心的作用,而右
边区域的βαβαβ结构形成Rossmann卷曲为结合NAD的区域。犛狋GAPC蛋白单体包含16个α螺旋,22个β折
叠,36个转角;犃狋GAPC2蛋白单体包含16个螺旋,22个β折叠,37个转角;而3e5rC有15个α螺旋,20个β折
叠,40个转角。图3B中1号和2号箭头分别表示这两种蛋白都有而3e5rC没有的2个β折叠和1个α螺旋,3
号箭头显示出两者在参与形成该螺旋的氨基酸数目上有差异,而三者在N末端和C末端均有差异。
2.2.3 犌犃犘犆基因序列的同源性比较 将不同科、属植物犌犃犘犆的核苷酸序列(CDS)以及氨基酸序列用DNA
MAN作同源比对。结果显示与马铃薯同科的植物中,普通蕃茄(犔狔犮狅狆犲狉狊犻犮狅狀犲狊犮狌犾犲狀狋狌犿)核苷酸序列同源性与
其最高为96.3%,其次是大叶烟草×花烟草(犖.犾犪狀犵狊犱狅狉犳犳犻犻×犖.狊犪狀犱犲狉犪犲)90.8%,野生型马铃薯(犛狅犾犪狀狌犿
犮犺犪犮狅犲狀狊犲)和普通烟草(犖犻犮狅狋犻犪狀犪狋犪犫犪犮狌犿)与之同源性较低为80%左右;而氨基酸序列间的同源性高于核苷酸
依次为98.6%,94.4%,89.3%,91.5%。犃狋犌犃犘犆2与犃狋犌犃犘犆1以及同科的芥菜(犅狉犪狊狊犻犮犪犼狌狀犮犲犪)、油菜(犅.
狀犪狆狌狊GAPC1和GAPC2)、白菜(犅.狉犪狆犪)等核苷酸序列同源性均为90%左右,而氨基酸同源性则高达98%~
96%。犛狋犌犃犘犆和犃狋犌犃犘犆2都与禾本科植物,如玉米、水稻、小麦(犜狉犻狋犻犮狌犿犪犲狊狋犻狏狌犿)等的同源性较低,核苷酸
为76%~78%,氨基酸为85%~87%。
根据多种植物犌犃犘犆基因的CDS序列建立如图4所示进化树,蕨类植物复活卷柏与其他被子植物处于两
大分支中,被子植物中的单子叶植物水稻、玉米、小麦、毛竹与其他双子叶植物处于不同分支,而双子叶植物中拟
南芥、油菜、白菜等十字花科在一个分支中,同属藜科的大洋洲滨藜和甜菜在同一分支,同属豆科的大豆和豌豆在
同一分支,其他科、属的植物也按此聚类,基本与现有的分类系统相符。但同为茄科的马铃薯栽培型、番茄、大叶
烟草×花烟草为一个分支,而马铃薯野生型、矮牵牛、普通烟草则出现在较远的另一分支。
2.3 犛狋犌犃犘犆和犃狋犌犃犘犆2基因植物表达载体构建
将植物表达载体pBI121和克隆载体分别酶切,并回收相应片段,纯化连接,构建了CaMV35S驱动犌犃犘犆基
因表达的载体。以重组子为模板PCR扩增犌犃犘犆基因,可扩增出约1000bp左右的片段(图5A),用犡犫犪Ⅰ和
犛犿犪Ⅰ双酶切重组质粒,也可获得1000bp左右的片段(图5B)。检测结果说明犌犃犘犆基因片段已插入表达载体
pBI121上,将获得的重组质粒分别命名为pBI121StGAPC和pBI121AtGAPC2。
采用冻融法将两种表达载体分别转入感受态农杆菌EHA105中,在含Kan的平板上获得单菌落,各选2个
单菌落提质粒,PCR扩增Kan抗性基因狀狆狋Ⅱ和目的基因犌犃犘犆,如图5所示,能扩增出预期的650和1000bp
左右的片段,说明获得了分别含pBI121StGAPC和pBI121AtGAPC2质粒的两种菌株。
3 讨论
胞质3磷酸甘油醛脱氢酶是由4个同型亚基组成的四聚体蛋白,Marri等[21]将拟南芥磷酸化GAPDH 基因
作了分类和表达研究,其中GapC为胞质内NAD依赖型的GAPDH,在糖酵解、糖质新生代谢中行使功能,由2
个基因犌犪狆犆1和犌犪狆犆2编码,两者cDNA相似性为89%,氨基酸序列相似性达97.9%,本实验克隆了其中的
犃狋犌犃犘犆2基因。目前认为造成克隆序列间差异的原因主要有品种差异、PCR扩增时的碱基错配、测序错误或者
发生中性突变[22]。本实验根据测序,筛选到了与GenBank公布的序列完全一致的克隆。而同时测序的几个克
隆中,与公布序列在不同位置有1个或2个碱基的变异,这可能是PCR扩增过程中的碱基错配造成。所克隆的
342第23卷第1期 草业学报2014年
图4 犌犃犘犆基因的进化树分析
犉犻犵.4 犘犺狔犾狅犵犲狀犲狋犻犮狋狉犲犲犪狀犪犾狔狊犻狊犫犪狊犲犱狅狀狀狌犮犾犲狅狋犻犱犲狅犳犌犃犘犆
 
图5 表达载体犘犆犚与酶切检测
犉犻犵.5 犐犱犲狀狋犻犳犻犮犪狋犻狅狀狅犳犲狓狆狉犲狊狊犻狅狀狏犲犮狋狅狉犫狔犲狀狕狔犿犲犱犻犵犲狊狋犻狅狀
 A:质粒PCR扩增目的片段GAPCgenePCRamplification;M:DNAMarkerⅢ;1:pBI121StGAPC;2:pBI121AtGAPC2;B:质粒酶切Enzyme
digestion;M:DNAMarkerDL15000TM;1,3:pBI121StGAPC和pBI121AtGAPC2质粒犡犫犪Ⅰ和犛犿犪Ⅰ双酶切 EnzymedigestionofpBI121
StGAPCandpBI121AtGAPC2by犡犫犪Ⅰand犛犿犪Ⅰ;2,4:pBI121StGAPC和pBI121AtGAPC2质粒犛犿犪Ⅰ单酶切EnzymedigestionofpBI121
StGAPCandpBI121AtGAPC2by犛犿犪Ⅰ.
犛狋犌犃犘犆和犃狋犌犃犘犆2核苷酸序列的相似性为84%,氨基酸序列相似性达92%,氨基酸转变会有疏水性变化,最
终会使3D结构产生差异,这些差异是否会造成犛狋犌犃犘犆和犃狋犌犃犘犆2功能间的差异还需进一步研究。
由于植物GAPC氨基酸序列间的同源性很高,遗传密码的兼并性可能会掩盖进化差别,降低进化树的可信
度,因此本实验利用犌犃犘犆的核苷酸序列构建进化树,发现与现有的分类系统基本相符,这与肖川等[6]对多个物
种犌犃犘犇犎 基因的分子进化研究结果相同,其研究表明叶绿体的犌犃犘犇犎 基因更接近原核生物,低等植物苔藓
442 ACTAPRATACULTURAESINICA(2014) Vol.23,No.1
较早分支,裸子植物和被子植物处于两大分支,被子植
图6 质粒转化农杆菌的犘犆犚验证
犉犻犵.6 犐犱犲狀狋犻犳犻犮犪狋犻狅狀狅犳犈犎犃105狋狉犪狀狊犳狅狉犿犲犱
犻狀狋狅狆犾犪狊犿犻犱犫狔犘犆犚
   M:DNAMarkerⅢ;1~2:pBI121StGAPC扩增犌犃犘犆基因 Ampli
fication犌犃犘犆geneofpBI121StGAPC;3~4:pBI121AtGAPC扩增
犌犃犘犆 基因 Amplification犌犃犘犆geneofpBI121AtGAPC2;5~6:
pBI121StGAPC扩增狀狆狋Ⅱ基因 Amplification狀狆狋Ⅱ geneofpBI121
StGAPC;7~8:pBI121AtGAPC2扩增狀狆狋Ⅱ基因 Amplification狀狆狋Ⅱ
geneofpBI121AtGAPC2.
物中单子叶植物和众多双子叶植物也处于不同分支,
认为犌犃犘犇犎 作为分子进化研究材料有较高可信度。
犌犃犘犆在同科、属间的同源性较高,说明在进化上较
为保守,这与其参与糖酵解这一基本生理活动有关。
对其他不同基因的研究所构建的进化树分析也表明各
基因的属内同源性较高可归于同一分支[2225]。研究表
明西伯利亚蓼的GAPDH 可提高酵母的抗盐碱胁迫
能力,并且该基因与盐生植物冰叶日中花的遗传距离
最近,推测逆境胁迫的作用使植物GAPDH基因的进
化方向和功能趋于相同[7]。本实验所构建的进化树
中,茄科植物出现了同属但不在一个分支,不同属反在
同一分支的现象,如马铃薯栽培型(茄属)、番茄栽培型
(番茄属)、大叶烟草×花烟草(烟草属)为一个分支,而
马铃薯野生型(茄属)、矮牵牛(矮牵牛属)、普通烟草
(烟草属)出现在另一分支。马铃薯栽培型在遗传上是四倍体而野生型是二倍体,可能对不同生长环境的适应及
染色体倍性的变化影响着犌犃犘犆基因的进化。进化树的构建是一个统计学问题,所构建出来的进化树只是对真
实进化关系的模拟,所以单凭犌犃犘犆序列建立的进化树还不能充分反映出马铃薯栽培型和野生型间的进化关
系。
GAPC是糖酵解中的关键酶之一,催化3磷酸甘油醛形成含高能磷酸键的1,3二磷酸甘油酸和NADH,前
者在下一步反应中生成3磷酸甘油酸和ATP,过去认为GAPC在细胞内表达稳定但近年来的研究表明GAPC
具有功能多样性。动物中的研究表明GAPC除参与糖酵解外,还参与细胞凋亡、调节 mRNA稳定性、氧化胁迫
应答、细胞自我吞噬、介导细胞信号转导、组蛋白基因调节等,并对GAPC在这些过程中的作用机理都有深入研
究和解释[26]。植物中发现GAPC能够响应多种逆境胁迫,但研究的深度滞后于动物。关于GAPC响应高温、低
温、活性氧、低磷、盐等多种逆境胁迫已被不同研究者所证实[812],而通过转基因技术也表明犌犃犘犆基因的过量表
达能提高植物对盐胁迫的抗性[1314],抑制表达则影响代谢和生长[27]。但目前关于GAPC在植物抗逆胁迫中如何
起作用的研究还不够深入,Zhang等[14]对水稻的研究表明犗狊犌犃犘犆3过量表达的转基因材料中,与水稻抗逆相关
的基因犆犪狋犃、犇犚犈犅2犃、犔犻狆9表达量高于对照,并且H2O2 含量低,植株的耐盐胁迫能力提高。这从一定程度上
解释了过量表达犌犃犘犆提高耐盐胁迫的原因。Gou等[28]利用双分子荧光互补充技术,发现H2O2 能促进拟南芥
GAPCs和磷脂酶D的相互作用从而增强磷脂酶D的活性,证明GAPC参与ABA诱导的逆境胁迫信号转导过
程,这一研究揭示了 GAPC的调节功能和作为逆境信号转导中的节点作用。本研究克隆了马铃薯和拟南芥
犌犃犘犆,并构建植物表达载体,这将为进一步通过转基因和酵母双杂等途径研究两者在马铃薯抗逆中的功能机制
奠定基础。
4 结论
本实验通过RTPCR技术,克隆了马铃薯和拟南芥的胞质3磷酸甘油醛脱氢酶基因的开放阅读框序列
(CDS),对其核苷酸和氨基酸序列以及蛋白结构作对比分析,并构建了植物表达载体,为下一步基因功能的研究
以及获得耐逆境胁迫的材料奠定基础。
参考文献:
[1] CerffR,ChambersSE.Subunitstructureofhigherplantglyceraldehyde3phosphatedehydrogenases(EC1.2.1.12andEC1.
2.1.13)[J].JournalofBiologicalChemistry,1979,254:60946098.
542第23卷第1期 草业学报2014年
[2] PlaxtonWC.Theorganizationandregulationofplantglycolysis[J].AnnualReviewofPlantPhysiologyandPlantMolecular
Biology,1996,47:185214.
[3] BrinkmannH,MartinezP,QuigleyF,犲狋犪犾.Endosymbioticoriginandcodonbiasofthenucleargeneforchloroplastglyceral
dehyde3phosphatedehydrogenasefrommaize[J].JournalofMolecularEvolution,1987,26(4):320328.
[4] RusselDA,SachsM M.Differentialexpressionandsequenceanalysisofthemaizeglyceraldehyde3phosphatedehydrogen
asegenefamily[J].ThePlantCel,1989,1(8):793803.
[5] ShihMC,HeinrichP,GoodmanHM.Cloningandchromosomalmappingofnucleargenesencodingchloroplastandcytosolic
glyceraldehyde3phosphatedehydrogenasefrom犃狉犪犫犻犱狅狆狊犻狊狋犺犪犾犻犪狀犪[J].Gene,1991,104(2):133138.
[6] 肖川,刘晓斐,詹树萱,等.水稻甘油醛3磷酸脱氢酶cDNA结构分析和分子进化[J].自然科学进展,1998,8(4):411
419.
[7] 李晓泽,刘关君,杨传平.西伯利亚蓼甘油醛3磷酸脱氢酶基因的cDNA克隆与序列分析[J].植物生理学通讯,2007,
43(1):4148.
[8] YangYJ,HawkBK,PengHP,犲狋犪犾.Stressresponsesandmetabolicregulationofclyceraldehyde3phosphatedehydrogen
asegenesin犃狉犪犫犻犱狅狆狊犻狊[J].PlantPhysiology,1993,101:209216.
[9] PenalozaE,GutierrezA,MartinezJ,犲狋犪犾.Differentialgeneexpressioninproteoidrootclustersofwhitelupin(犔狌狆犻狀狌狊犪犾
犫狌狊)[J].PlantPhysiology,2002,116:2836.
[10] WangXY,ChenYF,ZouJJ,犲狋犪犾.Involvementofacytoplasmicglyceraldehyde3phosphatedehydrogenaseGapC2in
lowphosphateinducedanthocyaninaccumulationin犃狉犪犫犻犱狅狆狊犻狊[J].ChineseScienceBuletin,2007,52(13):17641770.
[11] LaxaltAM,CassiaRO,SanlorentiPM,犲狋犪犾.Accumulationofcytosolicglyceraldehydes3phosphatedehydrogenaseRNA
underbiologicalstressconditionsandelicitortreatmentsinpotato[J].PlantMolecularBiology,1996,30(5):961972.
[12] NiuX,WangH,BressanRA,犲狋犪犾.Molecularcloningandexpressionofaglyceraldehyde3phosphatedehydrogenasegene
inadeserthalophyte,犃狋狉犻狆犾犲狓狀狌犿犿狌犾犪狉犻犪L[J].PlantPhysiology,1994,104(3):11051106.
[13] JeongMJ,ParkSC,ByunMO.Improvementofsalttoleranceintransgenicpotatoplantsbyglyceraldehydes3phosphate
dehydrogenasegenetransfer[J].MoleculesandCels,2001,12:185189.
[14] ZhangXH,RaoXL,ShiHT,犲狋犪犾.Overexpressionofacytosolicglyceraldehydes3phosphatedehydrogenasegene犗狊
犌犃犘犆3conferssalttoleranceinrice[J].PlantCelTissueandOrganCulture,2011,107:111.
[15] 门福义,刘梦云.马铃薯栽培生理[M].北京:中国农业出版社,1993.
[16] 王西瑶,朱涛,邹雪,等.缺磷胁迫增强了马铃薯的耐旱能力[J].作物学报,2009,35(5):875883.
[17] 杨乾,杨宏羽,刘玉汇,等.马铃薯适应低磷胁迫的生理生化效应[J].分子植物育种,2011,9(2):224229.
[18] 吴平,印莉萍,张立平,等.植物营养分子生理学[M].北京:科学出版社,2001.
[19] 王西瑶.与拟南芥突变体lpt1、lks1性状相关蛋白质的分离、鉴定与功能分析[D].北京:中国农业大学,2005.
[20] HajirezaeiMR,BiemeltS,PeiskerM,犲狋犪犾.Theinfluenceofcytosolicphosphorylatingglyceraldehydes3phosphatedehy
drogenase(GAPC)onpotatotubermetabolism[J].JournalofExperimentalBotany,2006,(23):115.
[21] MarriL,SparlaF,PupiloP,犲狋犪犾.Coordinatedgeneexpressionofphotosyntheticglyceraldehyde3phosphatedehydrogen
ase,phosphoribulokinase,andCP12in犃狉犪犫犻犱狅狆狊犻狊狋犺犪犾犻犪狀犪[J].JournalofExperimentalBotany,2005,56:7380.
[22] 牛继平,张金文,王旺田,等.马铃薯SGAs合成代谢途径末端SGT酶基因克隆及序列分析[J].草业学报,2012,21(3):
106116.
[23] 李玉坤,王学敏,王赞,等.东方山羊豆脱水蛋白基因的克隆及初步分析[J].草业学报,2012,21(1):176183.
[24] 李立芹,黄玉碧,王西瑶.马铃薯犠犚犓犢6基因的克隆\序列分析与原核表达研究[J].草业学报,2011,20(2):177183.
[25] 董洁,王学敏,王赞,等.紫花苜蓿二氢黄酮醇还原酶基因(MsDFR)的克隆与分析[J].草业学报,2012,21(2):123132.
[26] SiroverMA.Onthefunctionaldiversityofglyceraldehyde3phosphatedehydrogenase:Biochemicalmechanismsandregula
torycontrol[J].BiochimicaetBiophysicaActaGeneralSubjects,2011,1810:741751.
[27] RiusSP,CasatiP,IglesiasAA,犲狋犪犾.Characterizationof犃狉犪犫犻犱狅狆狊犻狊linesdeficientinGAPC1,acytosolicNADdepend
entglyceraldehyde3phosphatedehydrogenase[J].PlantPhysiology,2008,148:16551667.
[28] GouL,DevaiahSP,NarasimhanR,犲狋犪犾.Cytosolicglyceraldehyde3phosphatedehydrogensaesinteractwithphospholipase
642 ACTAPRATACULTURAESINICA(2014) Vol.23,No.1
dδtotransducehydrogenperoxidesignalsinthe犃狉犪犫犻犱狅狆狊犻狊responsetostress[J].ThePlantCel,2012,24:22002212.
犆犾狅狀犻狀犵犪狀犱犪狀犪犾狔狊犻狊狅犳狋犺犲犮狔狋狅狊狅犾犻犮犵犾狔犮犲狉犪犾犱犲犺狔犱犲3狆犺狅狊狆犺犪狋犲犱犲犺狔犱狉狅犵犲狀犪狊犲
(犌犃犘犆)犵犲狀犲犳狉狅犿犛狅犾犪狀狌犿狋狌犫犲狉狅狊狌犿犪狀犱犃狉犪犫犻犱狅狆狊犻狊狋犺犪犾犻犪狀犪
ZOUXue,ZHANGYe,WUMingyang,WANGXiyao
(ColegeofAgronomy,SichuanAgriculturalUniversity,Chengdu611130,China)
犃犫狊狋狉犪犮狋:Cytosolicglyceraldehyde3phosphatedehydrogenase(GAPC)takespartinglycolysisandisaclassi
calhousekeepingenzymeinbothprokaryotesandeukaryotes.IncreasingevidenceindicatesthatGAPCisin
volvedinvariousplantabioticandbioticstressresponsessuchassalt,lowphosphate,reactiveoxygenspecies
(ROS)andheatinphytophthorainfestans.The犃狋犌犃犘犆2genehasbeenfoundinresponsestolowphosphate
andosmoticstressin犃狉犪犫犻犱狅狆狊犻狊.AnalysisofsequencesandproteinstructureofGAPCcouldhelpinstudies
ofthefunctionandmechanismsofGAPCactinginplantstressresistance.ThefragmentsencodingGAPCwere
obtainedfrom犛狅犾犪狀狌犿狋狌犫犲狉狅狊狌犿 (StGAPC)and犃狉犪犫犻犱狅狆狊犻狊狋犺犻犪犾犻犪狀犪(犃狋犌犃犘犆2)seedlingsbyreversetran
scriptionpolymerasechainreaction(RT-PCR).TheCDSlengthsof犛狋犌犃犘犆and犃狋犌犃犘犆2were1017bp,
encodingaproteinsubunitof338aminoacids.Bioinformaticsanalysisdemonstratedthatthe犛狋犌犃犘犆and犃狋
犌犃犘犆2had84% similarityinnucleotidesequenceand92% similarityinaminoacidsequence.Molecular
weightsoftheproteinswere36.65and36.91kDa,thePI(theoreticalisoelectricpoint)were6.34and6.67,
andbothofthemhadstableproteinstructures.The3Dstructuresof犛狋犌犃犘犆and犃狋犌犃犘犆2werepredictedby
homologycomparativemodelingintheSwissModel,andshowedthatthe3Dstructureswerehighlysimilarto
犗狉狔狕犪狊犪狋犻狏犪GAPDHCchain(2e5rC).PredictedproteinscontainedconservedGAPCdomainswhichincluded
anNADbindingdomainconstructedbyβαβαβRossmanncoilandacatalyzeddomaininthesaddlesupersecond
arystructurecomposedofsevenβfoldings.CDSsequencesofGAPCshowedthatthesamefamilyofplantshas
ahighhomologyattributedtothesamebranchintheevolutionarytreeandthisevolutionaryrelationshipcorre
spondedtotheexistingclassificationsystem.pBI121wasconstructedasabasevectorandthegeneexpression
vectorintheplantwasdrivenbyaCaMV35Spromoter.Therecombinantvectorwasintroducedinto犃犵狉狅犫犪犮
狋犲狉犻狌犿strainEHA105andpositiveclonesweredeterminedbyPCR.Thisresearchisabasistoobtaintransgen
icmaterialsandforfurtherstudyonthefunctionofGAPC.
犓犲狔狑狅狉犱狊:犛狅犾犪狀狌犿狋狌犫犲狉狅狊狌犿;犃狉犪犫犻犱狅狆狊犻狊狋犺犪犾犻犪狀犪;cytosolicglyceraldehyde3phosphatedehydrogenase;
genecloning;bioinformaticanalysis
742第23卷第1期 草业学报2014年