留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

陆地棉低磷胁迫应答基因GhGDPD1的克隆与表达分析

孟超敏 耿翡翡 卿桂霞 张富厚 李雪林 刘逢举

孟超敏, 耿翡翡, 卿桂霞, 张富厚, 李雪林, 刘逢举. 陆地棉低磷胁迫应答基因GhGDPD1的克隆与表达分析[J]. 浙江农林大学学报, 2023, 40(4): 723-730. doi: 10.11833/j.issn.2095-0756.20220624
引用本文: 孟超敏, 耿翡翡, 卿桂霞, 张富厚, 李雪林, 刘逢举. 陆地棉低磷胁迫应答基因GhGDPD1的克隆与表达分析[J]. 浙江农林大学学报, 2023, 40(4): 723-730. doi: 10.11833/j.issn.2095-0756.20220624
MENG Chaomin, GENG Feifei, QING Guixia, ZHANG Fuhou, LI Xuelin, LIU Fengju. Cloning and expression analysis of low phosphorus stress response gene GhGDPD1 in Gossypium hirsutum[J]. Journal of Zhejiang A&F University, 2023, 40(4): 723-730. doi: 10.11833/j.issn.2095-0756.20220624
Citation: MENG Chaomin, GENG Feifei, QING Guixia, ZHANG Fuhou, LI Xuelin, LIU Fengju. Cloning and expression analysis of low phosphorus stress response gene GhGDPD1 in Gossypium hirsutum[J]. Journal of Zhejiang A&F University, 2023, 40(4): 723-730. doi: 10.11833/j.issn.2095-0756.20220624

陆地棉低磷胁迫应答基因GhGDPD1的克隆与表达分析

doi: 10.11833/j.issn.2095-0756.20220624
基金项目: 国家自然科学基金资助项目(31101184)
详细信息
    作者简介: 孟超敏(ORCID: 0000-0003-2783-3996),副教授,从事作物营养遗传育种研究。E-mail: chaominm@haust.edu.cn
  • 中图分类号: S330;Q781

Cloning and expression analysis of low phosphorus stress response gene GhGDPD1 in Gossypium hirsutum

  • 摘要:   目的  基于前期陆地棉‘新陆早 19’Gossypium hirsutum‘Xinluzao 19’根部低磷胁迫基因表达谱芯片差异表达序列数据分析,挖掘相关基因,并对其克隆与表达分析。  方法  克隆陆地棉‘新陆早 19’GhGDPD1基因并进行基因组DNA与cDNA测序分析,借助生物信息学方法分析GhGDPD1的基因结构和进化关系;采用半定量RT-PCR技术与实时荧光定量PCR (RT-qPCR)的方法检测该基因于根、茎、叶、花等4个组织的基因表达量的变化以及低磷胁迫下的表达模式。  结果  成功克隆GhGDPD1基因,开放阅读框序列长度为1 149 bp,共编码382个氨基酸,属于GDPD家族,其中存在一个保守结构域,名为GDPD_GDE5_like_1_plant。半定量RT-PCR和RT-qPCR试验结果均显示:GhGDPD1基因主要表达于根,中量表达于花和茎,微量表达于叶。该基因在受到低磷胁迫的刺激后,立即会对低磷胁迫做出应答,胁迫4 h相对表达量达到最高值。  结论  首次成功克隆到陆地棉‘新陆早 19’GhGDPD1基因,获得了GhGDPD1基因的组织表达以及低磷胁迫下的表达模式,为深入解析棉花GhGDPD1基因的生物学功能以及培育棉花磷高效利用新种质提供科学参考。图7表1参18
  • 图  1  GhGDPD1基因编码序列的克隆

    Figure  1  Cloning of GhGDPD1 CDS

    图  2  GhGDPD1基因结构分析

    Figure  2  GhGDPD1 gene structure analysis

    图  3  GhGDPD1蛋白的三级结构预测分析结果

    Figure  3  Predicted three-dimension structure of GhGDPD1 protein

    图  4  GhGDPD1与其他物种GDPD家族同源蛋白的多序列比对分析结果

    Figure  4  Multi-alignment of GhGDPD1 amino acid sequence with other GDPD family in different plants

    图  5  GhGDPD1蛋白系统进化树分析

    Figure  5  GhGDPD1 phylogenetic tree analysis

    图  6  半定量RT-PCR分析GhGDPD1在棉花不同组织中的表达

    Figure  6  Expression of GhGDPD1 in different tissues of cotton by semi-quantitative RT-PCR results

    图  7  GhGDPD1基因在不同组织(A)及低磷胁迫处理下(B)的表达模式

    Figure  7  Expression analysis of GhMGD1 in different tissues (A) and low phosphorus stress (B)

    表  1  试验所用引物信息

    Table  1.   Primers used in the study

    引物引物序列(5′→3′)引物引物序列(5′→3′)
    Sense primer ATTTTTCCCTCTCTTACTCTATCCC GhActin-F ATCCTCCGTCTTGACCTTG
    Anti-sense primer GTAGGGACAAGTTAGTGGTGTATCA GhActin-R TGTCCGTCAGGCAACTCAT
    通用引物M13F TGTAAAACGACGGCCAGT GhGDPD1-F TTCTCTGTCTCTCTACTCGTCTCGT
    通用引物M13R CAGGAAACAGCTATGACC GhGDPD1-R TCTATGCCCTATTACCAAAAACTTC
    下载: 导出CSV
  • [1] DISSANAYAKA D M S B, PLAXTON W C, LAMBERS H, et al. Molecular mechanisms underpinning phosphorus-use efficiency in rice [J]. Plant,Cell &Environment, 2018, 41(7): 1483 − 1496.
    [2] HINSINGER P. Bioavailability of soil inorganic P in the rhizosphere as affected by root-induced chemical changes: a review [J]. Plant and Soil, 2001, 237(2): 173 − 195.
    [3] KYLA S, WILLIAM P. Molecular mechanisms of phosphorus metabolism and transport during leaf senescence [J]. Plants, 2015, 4(4): 773 − 798.
    [4] VENEKLAAS E J, LAMBERS H, BRAGG J, et al. Opportunities for improving phosphorus-use efficiency in crop plants [J]. New Phytologist, 2012, 195(2): 306 − 320.
    [5] PLAXTON W, TRAN H T. Metabolic adaptations of phosphate-starved plants [J]. Plant Physiology, 2011, 156(3): 1006 − 1015.
    [6] MEHRA P, GIRI J. Rice and chickpea GDPDs are preferentially influenced by low phosphate and CaGDPD1 encodes an active glycerophosphodiester phosphodiesterase enzyme [J]. Plant Cell Reports, 2016, 35(8): 1699 − 1717.
    [7] ZHENG Bin, CHEN Dan, FARQUHAR M G. MIR16, a putative membrane glycerophosphodiester phosphodiesterase, interacts with RGS16 [J]. Proceedings of the National Academy of Sciences of the United States of America, 2000, 97(8): 3999 − 4004.
    [8] CORDA D, MOSCA M G, OHSHIMA N, et al. The emerging physiological roles of the glycerophosphodiesterase family [J]. The FEBS Journal, 2014, 281(4): 998 − 1016.
    [9] LEMIEUX M J, HUANG Yafei, WANG Daneng. Glycerol-3-phosphate transporter of Escherichia coli: structure, function and regulation [J]. Research in Microbiology, 2004, 155(8): 623 − 629.
    [10] 王欣. 拟南芥AtGDPD5基因在低磷环境下功能的研究[D]. 哈尔滨: 东北农业大学, 2011.

    WANG Xin. Studies on the Function of Arabidopsis AtGDPD5 Gene in Low Phosphorus Environment [D]. Harbin: Northeast Agricultural University, 2011.
    [11] van der REST B, BOISSON A M, GOUT E, et al. Glycerophosphocholine metabolism in higher plant cells. evidence of a new glyceryl-phosphodiester phosphodiesterase [J]. Plant Physiology, 2002, 130(1): 244 − 255.
    [12] MEHRA P, PANDEY B K, VERMA L, et al. A novel glycerophosphodiester phosphodiesterase improves phosphate deficiency tolerance in rice [J]. Plant,Cell &Environment, 2019, 42(4): 1167 − 1179.
    [13] 王桂峰, 张杰, 王安琪. 全国棉花生产格局时景下山东省棉花生产保护区支撑体系构建[J]. 山东农业科学, 2020, 52(5): 130 − 135.

    WANG Guifeng, ZHANG Jie, WANG Anqi. Construction of support system for cotton production protection reserve in Shandong Province under national cotton production pattern [J]. Shandong Agricultural Sciences, 2020, 52(5): 130 − 135.
    [14] 王桂峰, 徐勤青, 王安琪. 近二十年全国及山东植棉变化与山东省棉花生产保护区支撑体系构建思考[J]. 棉花科学, 2020, 42(1): 3 − 13.

    WANG Guifeng, XU Qinqing, WANG Anqi. Changes of cotton planting in China and Shandong Province in recent 20 years and the construction of supporting system of cotton production protection area in Shandong Province [J]. Cotton Sciences, 2020, 42(1): 3 − 13.
    [15] 张敏, 盛建东, 白灯莎·买买提艾力, 等. 不同磷效率棉花根系形态和磷酸酶活性对供磷强度的响应[J]. 棉花学报, 2017, 29(3): 283 − 291.

    ZHANG Min, SHENG Jiandong, Baidengsha Maimaitiaili, et al. Response of root morphology and phosphatase activity of cotton to phosphorus supply [J]. Cotton Science, 2017, 29(3): 283 − 291.
    [16] PAZ-ARES J, PUGA M I, ROJAS-TRIANA M, et al. Plant adaptation to low phosphorus availability: core signaling, crosstalks, and applied implications [J]. Molecular Plant, 2022, 15(1): 104 − 124.
    [17] CHENG Yuxiang, ZHOU Wenbin, SHEERY N I E, et al. Characterization of the Arabidopsis glycerophosphodiester phosphodiesterase (GDPD) family reveals a role of the plastid-localized AtGDPD1 in maintaining cellular phosphate homeostasis under phosphate starvation [J]. The Plant Journal, 2011, 66(5): 781 − 795.
    [18] CHENG Lingyun, BUCCIARELLI B, LIU Junqi, et al. White lupin cluster root acclimation to phosphorus deficiency and root hair development involve unique glycerophosphodiester phosphodiesterases [J]. Plant Physiology, 2011, 156(3): 1131 − 1148.
  • [1] 杨勇, 张俊红, 韩潇, 张毓婷, 杨琪, 童再康.  闽楠bZIP基因家族鉴定和脱落酸处理下的表达分析 . 浙江农林大学学报, 2024, 41(2): 275-285. doi: 10.11833/j.issn.2095-0756.20230342
    [2] 王书伟, 周明兵.  毛竹ICE基因家族的全基因组鉴定及低温胁迫下的表达模式分析 . 浙江农林大学学报, doi: 10.11833/j.issn.2095-0756.20230445
    [3] 尚林雪, 王群, 张国哲, 赵雨, 顾翠花.  紫薇LiCMB1基因的克隆及表达特性分析 . 浙江农林大学学报, 2023, 40(2): 330-337. doi: 10.11833/j.issn.2095-0756.20220333
    [4] 洪方蕾, 陆瑶, 俞世姣, 胡芷诺, 缪云锋, 钟诗蔚, 赵宏波.  桂花OfABFs基因克隆和表达分析 . 浙江农林大学学报, 2023, 40(3): 481-491. doi: 10.11833/j.issn.2095-0756.20220264
    [5] 卓娟, 侯丹, 林新春.  毛竹PhebHLH6基因克隆及表达分析 . 浙江农林大学学报, 2023, 40(4): 731-737. doi: 10.11833/j.issn.2095-0756.20220553
    [6] 兰智鑫, 侯丹, 吴蔼民, 林新春.  毛竹PeCIGRs基因的克隆及表达分析 . 浙江农林大学学报, 2023, 40(5): 982-990. doi: 10.11833/j.issn.2095-0756.20220761
    [7] 王倩清, 张毓婷, 张俊红, 刘慧, 童再康.  闽楠PLR基因家族鉴定及响应激素的表达分析 . 浙江农林大学学报, 2022, 39(6): 1173-1182. doi: 10.11833/j.issn.2095-0756.20220351
    [8] 陈雪冰, 刘聪, 程赫, 姜廷波, 夏德安, 魏志刚.  毛果杨ZHD家族全基因组水平鉴定及在干旱胁迫下的表达分析 . 浙江农林大学学报, 2022, 39(3): 465-474. doi: 10.11833/j.issn.2095-0756.20210373
    [9] 孟超敏, 耿翡翡, 卿桂霞, 周佳敏, 张富厚, 刘逢举.  陆地棉磷高效基因GhMGD3的克隆与表达分析 . 浙江农林大学学报, 2022, 39(6): 1203-1211. doi: 10.11833/j.issn.2095-0756.20220145
    [10] 陆丹迎, 程少禹, 章颖佳, 刘志高, 金梦婷, 董彬, 张寿洲, 彭豪, 戴梦怡, 王卓为, 赵宏波, 申亚梅.  景宁木兰PIF转录因子的生物信息学分析及极端遮阴条件下的表达模式 . 浙江农林大学学报, 2021, 38(3): 445-454. doi: 10.11833/j.issn.2095-0756.20200488
    [11] 庞天虹, 钱婕妤, 付建新, 顾翠花, 张超.  桂花己糖激酶基因家族成员的序列及表达分析 . 浙江农林大学学报, 2021, 38(2): 225-234. doi: 10.11833/j.issn.2095-0756.20200370
    [12] 王楠楠, 董彬, 杨丽媛, 赵宏波.  梅花2个PmWRKY2基因克隆及在逆境胁迫下的表达模式 . 浙江农林大学学报, 2021, 38(4): 812-819. doi: 10.11833/j.issn.2095-0756.20200706
    [13] 蒋琦妮, 付建新, 张超, 董彬, 赵宏波.  桂花OfAP1基因的克隆及表达分析 . 浙江农林大学学报, 2019, 36(4): 664-669. doi: 10.11833/j.issn.2095-0756.2019.04.005
    [14] 武丹阳, 杨洋, 李慧玉.  3个白桦BpBEE基因的克隆与表达分析 . 浙江农林大学学报, 2017, 34(1): 137-144. doi: 10.11833/j.issn.2095-0756.2017.01.019
    [15] 李冰冰, 刘国峰, 魏书, 黄龙全, 张剑韵.  烟草NtPLR1基因克隆与表达分析 . 浙江农林大学学报, 2017, 34(4): 581-588. doi: 10.11833/j.issn.2095-0756.2017.04.003
    [16] 侯传明, 郑雅文, 王正加, 徐英武.  山核桃MADS-like基因的克隆与分析 . 浙江农林大学学报, 2015, 32(1): 33-39. doi: 10.11833/j.issn.2095-0756.2015.01.005
    [17] 庞景, 童再康, 黄华宏, 林二培, 刘琼瑶.  杉木纤维素合成酶基因CesA的克隆及表达分析 . 浙江农林大学学报, 2015, 32(1): 40-46. doi: 10.11833/j.issn.2095-0756.2015.01.006
    [18] 沈辰, 裘佳妮, 黄坚钦.  山核桃COP1 E3连接酶的全长克隆及表达分析 . 浙江农林大学学报, 2014, 31(6): 831-837. doi: 10.11833/j.issn.2095-0756.2014.06.002
    [19] 杨希宏, 黄有军, 陈芳芳, 黄坚钦.  山核桃FLOWERING LOCUS C同源基因鉴定与表达分析 . 浙江农林大学学报, 2013, 30(1): 1-8. doi: 10.11833/j.issn.2095-0756.2013.01.001
    [20] 王国立, 安华明, 秦巧平, 李孟娇, 刘真真, 陈佳莹, 周倩, 张岚岚.  柑橘果实成熟特异基因CsPMEI/InvI的克隆与序列分析 . 浙江农林大学学报, 2013, 30(3): 336-342. doi: 10.11833/j.issn.2095-0756.2013.03.005
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.20220624

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2023/4/723

图(7) / 表(1)
计量
  • 文章访问数:  183
  • HTML全文浏览量:  43
  • PDF下载量:  308
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-09-28
  • 修回日期:  2023-03-22
  • 录用日期:  2023-03-24
  • 网络出版日期:  2023-07-13
  • 刊出日期:  2023-08-20

陆地棉低磷胁迫应答基因GhGDPD1的克隆与表达分析

doi: 10.11833/j.issn.2095-0756.20220624
    基金项目:  国家自然科学基金资助项目(31101184)
    作者简介:

    孟超敏(ORCID: 0000-0003-2783-3996),副教授,从事作物营养遗传育种研究。E-mail: chaominm@haust.edu.cn

  • 中图分类号: S330;Q781

摘要:   目的  基于前期陆地棉‘新陆早 19’Gossypium hirsutum‘Xinluzao 19’根部低磷胁迫基因表达谱芯片差异表达序列数据分析,挖掘相关基因,并对其克隆与表达分析。  方法  克隆陆地棉‘新陆早 19’GhGDPD1基因并进行基因组DNA与cDNA测序分析,借助生物信息学方法分析GhGDPD1的基因结构和进化关系;采用半定量RT-PCR技术与实时荧光定量PCR (RT-qPCR)的方法检测该基因于根、茎、叶、花等4个组织的基因表达量的变化以及低磷胁迫下的表达模式。  结果  成功克隆GhGDPD1基因,开放阅读框序列长度为1 149 bp,共编码382个氨基酸,属于GDPD家族,其中存在一个保守结构域,名为GDPD_GDE5_like_1_plant。半定量RT-PCR和RT-qPCR试验结果均显示:GhGDPD1基因主要表达于根,中量表达于花和茎,微量表达于叶。该基因在受到低磷胁迫的刺激后,立即会对低磷胁迫做出应答,胁迫4 h相对表达量达到最高值。  结论  首次成功克隆到陆地棉‘新陆早 19’GhGDPD1基因,获得了GhGDPD1基因的组织表达以及低磷胁迫下的表达模式,为深入解析棉花GhGDPD1基因的生物学功能以及培育棉花磷高效利用新种质提供科学参考。图7表1参18

English Abstract

孟超敏, 耿翡翡, 卿桂霞, 张富厚, 李雪林, 刘逢举. 陆地棉低磷胁迫应答基因GhGDPD1的克隆与表达分析[J]. 浙江农林大学学报, 2023, 40(4): 723-730. doi: 10.11833/j.issn.2095-0756.20220624
引用本文: 孟超敏, 耿翡翡, 卿桂霞, 张富厚, 李雪林, 刘逢举. 陆地棉低磷胁迫应答基因GhGDPD1的克隆与表达分析[J]. 浙江农林大学学报, 2023, 40(4): 723-730. doi: 10.11833/j.issn.2095-0756.20220624
MENG Chaomin, GENG Feifei, QING Guixia, ZHANG Fuhou, LI Xuelin, LIU Fengju. Cloning and expression analysis of low phosphorus stress response gene GhGDPD1 in Gossypium hirsutum[J]. Journal of Zhejiang A&F University, 2023, 40(4): 723-730. doi: 10.11833/j.issn.2095-0756.20220624
Citation: MENG Chaomin, GENG Feifei, QING Guixia, ZHANG Fuhou, LI Xuelin, LIU Fengju. Cloning and expression analysis of low phosphorus stress response gene GhGDPD1 in Gossypium hirsutum[J]. Journal of Zhejiang A&F University, 2023, 40(4): 723-730. doi: 10.11833/j.issn.2095-0756.20220624
  • 磷(P)是植物生长发育所必需的大量营养元素,对细胞组成、能量生成、新陈代谢和信号转导等关键细胞功能都具有重要作用[1]。由于土壤中磷的质量分数及流动性非常低,无机磷成为全球作物生产力的常见限制因素[2],缺磷是作物生产力损失的主要原因之一,因此,提高磷素利用效率对可持续粮食生产至关重要。为了提高磷的利用率,植物通过复杂的反应体系以改善从环境中获取磷素的途径,并在植物内循环中加以利用[3]。单纯地增加磷的吸收只能使磷从环境中输出,因此提高磷从衰老器官向幼嫩器官和发育器官的再利用效率是提高作物磷素利用效率的一个重要发展方向[4]。在许多植物中,磷脂降解是磷再活化的主要驱动因素[5]

    甘油磷酸二酯磷酸二酯酶(GDPD)将甘油磷酸二酯水解为甘油-3-磷酸和相应的醇,是甘油磷酸二酯代谢途径中的重要成员。这类基因被认为在磷酸盐稳态中发挥重要作用[6]。在原核生物和真核生物中都存在GDPD,它在各种生理过程中都发挥着重要作用[7],因此该蛋白具有生化和生物功能多样性。在哺乳动物中,GDPD与G蛋白信号转导、磷酸肌醇代谢、渗透保护、运动神经元分化、细胞骨架组织、成骨细胞分化和成肌分化有关[8]。第1个GDPD编码基因在大肠埃希菌Escherichia coli中被鉴定,它积极参与甘油和甘油-3-磷酸的摄取和代谢,因此在碳水化合物代谢和磷脂生物合成中发挥重要作用[9]。有关植物中GDPD的内容相对较少[10]。GDPD首先在磷饥饿胡萝卜Daucus carota var. sativa、美国梧桐Platanus occidentalis和拟南芥Arabidopsis thaliana细胞悬浮培养物的液泡和细胞壁中被发现[1011]。细胞膜以磷脂的形式占细胞有机磷的三分之一,因此,在缺磷条件下从膜磷脂中再活化磷可以成为提高磷利用效率的重要策略[12]

    中国的棉花Gossypium栽种始于宋元时期,通过引进新的品种打开了一种区域化、大面积的棉花种植模式[13],此后不断引入种质资源、培育棉花品种,逐步改善了棉花的纤维品质,提高了中国棉花的生产率及效益。1978—2019年中国棉花种植面积逐年减少,中国农业部要求致力提高单产、品质,增加效益[14]。因此必须不断优化棉花品种,提升棉花品质,而缺磷是作物生产力损失的主要原因之一。本研究以陆地棉‘新陆早19’ Gossypium hirsutum ‘Xinluzao 19’为试验材料,克隆了GhGDPD1基因。采用生物信息学工具分析其编码蛋白理化性质,并对其进行表达特性分析,旨在为深入解析棉花GhGDPD1基因的生物学功能提供科学参考,并为培育磷高效利用的棉花新品种提供基因资源。

    • 采用前期研究筛选出的磷高效陆地棉品种‘新陆早19’为试验材料,在河南科技大学农学院(34°59′N,112°42′E)进行大田种植,选取长势一致的‘新陆早19’植株采集根、茎、叶、花组织,用于分析GhGDPD1基因在4个组织中的表达状况。

      将‘新陆早19’种子播种于含干净湿润细沙的塑料盆内,28 ℃恒温培养箱(光照14 h、黑暗10 h)培养至棉苗三叶期,选择长势良好且一致的幼苗移至水培盆中,用 1/2 Hoagland 营养液培养1周后分为2个处理水平,即适磷处理(SP: 1.00 mmol·L−1)和低磷处理(LP: 0.01 mmol·L−1),磷源为磷酸二氢钾(KH2PO4)。为保证K+的浓度一致,以1.0 mmol·L−1为标准,低磷营养液中以氯化钾(KCl)补齐K+,其他营养成分含有2.0 mmol·L−1 Ca(NO3)2·4H2O、2.5 mmol·L−1KNO3、0.5 mmol·L−1NH4NO3、1.4 mmol·L−1MgSO4·7H2O、1.0×10−3 mmol·L−1 ZnSO4·7H2O、1.0×10−3 mmol·L−1MnSO4·H2O、1.0×10−4 mmol·L−1CuSO4·5H2O、 1.0×10−2 mmol·L−1H3BO3、0.5×10−5 mmol·L−1 (NH4)6MO7O24·4H2O、0.1 mmol·L−1EDTA-FeNa[15]。分别处理0、4、12、24、72 h后,选取3株生长一致的‘新陆早19’植株,并混合取其根部组织,然后迅速置于液氮中,−80 ℃保存备用。每日调节营养液 ,使pH保持在 6.5 左右,每3 d更换1次营养液。

    • 根据前期陆地棉根部低磷胁迫基因表达谱芯片差异表达基因序列进行分析,在美国国家生物技术信息中心(NCBI)网站的EST数据库中检索差异表达基因序列,得到该基因的相似序列。将所得相似性序列(覆盖率>50%,相似度>90%)使用DNASTAR的Seqman进行拼接得到重叠群(conting),将所得序列继续检索与拼接直至没有新的相似序列出现,所得即为结果序列重叠群,利用ORFfinder在线平台查找开放阅读框,进行目标基因GhGDPD1的克隆与分析。

      根据ORFfinder在线平台查找的开放阅读框,使用Primer 5.0软件设计引物,由生工生物工程有限公司合成(表1)。

      表 1  试验所用引物信息

      Table 1.  Primers used in the study

      引物引物序列(5′→3′)引物引物序列(5′→3′)
      Sense primer ATTTTTCCCTCTCTTACTCTATCCC GhActin-F ATCCTCCGTCTTGACCTTG
      Anti-sense primer GTAGGGACAAGTTAGTGGTGTATCA GhActin-R TGTCCGTCAGGCAACTCAT
      通用引物M13F TGTAAAACGACGGCCAGT GhGDPD1-F TTCTCTGTCTCTCTACTCGTCTCGT
      通用引物M13R CAGGAAACAGCTATGACC GhGDPD1-R TCTATGCCCTATTACCAAAAACTTC
    • 取液氮速冻整株‘新陆早19’植株在研钵中迅速研磨成粉末转移至离心管,利用CTAB法提取‘新陆早19’的基因组DNA,向其中添加200 μL的TE缓冲液溶解后置于−20 ℃下保存备用。

      以‘新陆早19’的基因组DNA为模板设计如下PCR扩增反应体系(20.0 μL)。冰上操作:2×M5 HiPer plus Taq HiFi PCR mix (with blue dye) 10.0 μL,Sense primer (10.0 μmol·L−1) 0.5 μL,Anti-sense primer (10.0 μmol·L−1) 0.5 μL,Template DNA 0.5 μL,Nuclease-free ddH2O 8.5 μL。反应程序为:95 ℃ 3 min;94 ℃ 25 s,55 ℃ 25 s,72 ℃ 50 s,35个循环;72 ℃ 5 min;−4 ℃低温保存。120 V,25 min,质量分数为1%琼脂糖凝胶电泳检测。

      将目的DNA片段进行胶回收纯化后连接pTOPO-T载体,具体体系为(5.0 μL):M5 HiPer pTOPO-TA Vector 0.5 μL,10×Enhancer 0.5 μL,纯化后的PCR产物 1.8 μL,灭菌水 2.2 μL。取5.0 μL连接产物转化50.0 μL大肠埃希菌Escherichia coli DH5ɑ感受态细胞,将菌液涂布在含氨苄青霉素的LB固体培养基上,37 ℃培养过夜。随机挑取单菌落,扩大培养4 h。以所得菌液作模板进行菌液PCR反应,体系如下(20.0 μL):2× M5 HiPer plus Taq HiFi PCR mix (with blue dye) 10.0 μL,通用引物M13F 0.5 μL,通用引物M13R 0.5 μL,Template DNA 0.5 μL,Nuclease-free ddH2O 8.5 μL。反应程序为:95 ℃ 3 min;94 ℃ 25 s,55 ℃ 25 s,72 ℃ 50 s,35个循环;72 ℃ 5 min;−4 ℃低温保存。检测符合后送至生工生物工程有限公司测序。

    • 取液氮速冻整株‘新陆早19’植株在研钵中迅速充分研磨成粉末,后续步骤依照所用试剂盒说明书进行,150 V,15 min,质量分数为1%琼脂糖凝胶电泳检测。参考试剂盒说明书完成cDNA的合成,−20 ℃保存备用。

      以‘新陆早19’的cDNA为模板设计如下PCR扩增反应体系(20.0 μL)。冰上操作:2×M5 HiPer plus Taq HiFi PCR mix (with blue dye) 10.0 μL,Sense primer (10.0 μmol·L−1) 0.5 μL,Anti-sense primer (10.0 μmol·L−1) 0.5 μL,cDNA 1.0 μL,Nuclease-free ddH2O 8.0 μL。PCR反应程序为:95 ℃ 3 min;94 ℃ 25 s,56 ℃ 25 s,72 ℃ 45 s,35个循环;72 ℃ 5 min;−4 ℃低温保存。120 V,25 min,质量分数为1%琼脂糖凝胶电泳检测。

      将目的DNA片段进行胶回收纯化后连接pTOPO-T载体,体系如下(5.0 μL):M5 HiPer pTOPO-TA Vector 0.5 μL,10×Enhancer 0.5 μL,纯化后的PCR产物 1.5 μL,灭菌水 2.5 μL。加入50.0 μL大肠埃希菌DH5ɑ感受态细胞转化,挑菌进行菌液PCR反应。电泳检测符合后送至生工生物工程有限公司进行测序。

    • 使用在线平台及软件对基因进行生物信息学分析,预测分析基因的结构、性质等。

    • 利用半定量RT-PCR技术,分析GhGDPD1基因在根、茎、叶、花4个组织中的表达状况。以‘新陆早19’各组织的cDNA为模板,GhActin作内参基因,所用引物见表1。设计如下PCR扩增反应体系(20.0 μL)。冰上操作:2×M5 HiPer plus Taq HiFi PCR mix (with blue dye) 10.0 μL,Sense primer (10.0 μmol·L−1) 0.5 μL,Anti-sense primer (10.0 μmol·L−1) 0.5 μL,cDNA 1.0 μL,Nuclease-free ddH2O 8.0 μL。PCR反应程序为:95 ℃ 3 min;94 ℃ 25 s,60 ℃ 25 s,72 ℃ 10 s,35个循环;72 ℃ 5 min;−4 ℃低温保存。120 V,25 min,质量分数为1%琼脂糖凝胶电泳检测。

    • 采用RT-qPCR技术检测GhGDPD1基因在不同组织中及低磷胁迫处理下的表达模式。以GhActin作为内参基因,所用引物见表1。采用SYBR® Green Pro Taq HS预混型qPCR试剂盒进行荧光定量,使用的荧光定量PCR仪器型号是CFX96,数据分析采用2−ΔΔCt法。

    • 以实验室前期低磷胁迫差异表达序列为探针,检索发现并下载了10条相似序列,利用DNASTAR将所得序列全部拼接得到新的重叠群,对新得到的序列再次进行检索,在第1次检索序列的基础上又增加了1条新的序列,将11条基因序列全部拼接,得到了1个新的重叠群,序列长度为1 462 bp。

    • 根据得到的克隆重叠群,横跨GhGDPD1基因的编码序列(coding sequence,CDS)区设计特异性引物,以‘新陆早19’植株的DNA和cDNA为模板,克隆得到GhGDPD1基因的CDS (图1),其条带符合目的条带大小。通过无缝克隆连接到pTOPO-T载体上,挑选阳性克隆测序。获得的DNA序列大于2 000 bp (图1,泳道2),与重叠群序列的一致性较高,均大于89%。获得cDNA序列1 200 bp (图1,泳道1),与重叠群序列比对一致性为99%。开放阅读框为长度为1 149 bp,共编码382个氨基酸。该编码区具有GDPD_GDE5_like_1_plant结构域,因此GhGDPD1基因是GDPD家族成员。

      图  1  GhGDPD1基因编码序列的克隆

      Figure 1.  Cloning of GhGDPD1 CDS

    • 将DNA测序结果序列与cDNA测序结果的最大开放阅读框序列导入Gene Structure Display Server 2.0,对该基因组序列与编码序列进行分析,表明该基因编码序列包含5′与3′非编码序列,共有5个内含子,6个外显子(图2)。

      图  2  GhGDPD1基因结构分析

      Figure 2.  GhGDPD1 gene structure analysis

    • 将基因编码蛋白导入ExPASy-ProtParam tool在线软件,分析可知:该基因共编码了382个氨基酸,所编码蛋白的分子式为C1930H3030N506O580S12,脂肪族氨基酸指数为94.84,分子量为42 987.05 Da,等电点为5.17,属于酸性蛋白。其中亮氨酸的相对含量占比最大,为10.2%,共39个。在该蛋白中带正电荷的氨基酸有38个,带负电荷的有54个。氨基酸组分会直接影响蛋白的亲疏水性,总亲水性平均系数为−0.179,不稳定系数为41.79,是亲水性不稳定蛋白。ProtScale analysis在线软件显示:此蛋白中精氨酸亲水性最强,亲水指数为−4.500 0,异亮氨酸疏水性最强,亲水指数为4.500 0。

    • 在NPS@:SOPMA secondary structure prediction网站预测该蛋白二级结构,其中占比最大的α-螺旋占40.84%,包含156个氨基酸,其次是无规则卷曲占比为38.22%,包含146个氨基酸。此外,此蛋白还含有延伸链(extended strand)和β-转角(beta turn),其中延伸链占比较多,为15.18%,有58个氨基酸,β-转角仅有22个氨基酸,占比为5.76%,是最少的。表明此蛋白的二级结构由α-螺旋和无规则卷曲占据主体地位所组成,各结构在各个区间分布较为均匀。由于无规则卷曲占比较多,预测其结构比较复杂。

      运用TMHMM-2.0分析发现:该基因所编码蛋白不含有跨膜结构域,属于非跨膜蛋白,382个氨基酸均位于细胞膜表面,并未形成跨膜螺旋区。通过SignalP-5.0进行信号肽分析,可知:GhGDPD1基因编码蛋白有信号肽的概率是0.000 3,其他的可能性则高达0.999 7。382个氨基酸中并未出现典型的信号肽趋势,即该蛋白不存在信号肽。

      在Plant-mPLoc在线网站内输入基因编码序列进行基因编码蛋白的亚细胞定位预测,预测定位在细胞膜(Cell membrane)。预测结果与跨膜结构预测分析结果相同,均为该蛋白位于膜外,推测该基因可能直接在细胞膜发挥作用。

    • 通过SWISS-MODEL Interactive Workspace进行同源建模,预测蛋白三级结构(图3)。评价同源建模 (QMEAN)结果为0.57,全球性模型质量估测 (GMQE)值为0.44,预测其含有甘油磷酸二酯磷酸二酯酶,与腾冲嗜热菌Thermoanaerobacter tengcongensis MB4的甘油磷酸二酯磷酸二酯酶的晶体结构类似。蛋白的三级结构与二级结构分析结果一致,均以α螺旋和无规则卷曲为主体。

      图  3  GhGDPD1蛋白的三级结构预测分析结果

      Figure 3.  Predicted three-dimension structure of GhGDPD1 protein

      将获得的基因编码序列提交到STRING 11.0在线网站,预测该基因可能的互作网络。分析结果显示:GhGDPD1基因编码蛋白能够预测到蛋白质的相互作用网络,因此预测GhGDPD1基因可能是通过多个蛋白互作进行调控。

      N-糖基化对蛋白稳定及功用均会产生极其重要的影响。将获得的基因编码序列提交到NetNGlyc 1.0在线软件,可知该蛋白存在3个潜在的N-糖基化位点(分别位于7~10、21~24、116~119氨基酸区段)。磷酸化位点与蛋白的功能、结构都有很大的关系。通过NetPhos 3.1在线网站分析可知:此蛋白中占比最大的是丝氨酸和苏氨酸磷酸化位点。GhGDPD蛋白除了可以调控蛋白的活性以外,在细胞的信号转导中也占据了重要的位置。

    • 使用NCBI的保守结构域数据库(CDD-search)查找GhGDPD1氨基酸序列的保守结构域。分析结果显示:GhGDPD1基因所编码蛋白含有1个特异匹配,为GDPD_GDE5_like_1_plant,匹配所属的超家族是PI-PLCc_GDPD_SF基因家族。GhGDPD1基因与其他物种的GDPD1氨基酸序列进行多重比对分析,同源性较高,物种间GDPD蛋白氨基酸长度基本一致,同源性达84.12%(图4)。

      图  4  GhGDPD1与其他物种GDPD家族同源蛋白的多序列比对分析结果

      Figure 4.  Multi-alignment of GhGDPD1 amino acid sequence with other GDPD family in different plants

      在NCBI数据库利用BLASTp检索氨基酸序列的同源序列,下载了榴莲Durio zibethinues、木槿Hibiscus syriacus、可可树Theobroma cacao、哥伦比亚锦葵Herrania umbratica、长蒴黄麻Corchorus olitorius、麻疯树Jatropha curcas、芒果Mangifera indica、胡杨Populus euphratica、开心果Pistacia vera、白梨Pyrus bretschneideri、金钱橘Citrus clementina等11条不同物种的蛋白序列。通过软件MEGA 5构建系统进化树(图5),‘新陆早19’蛋白序列与木槿的亲缘关系最近,其次与可可树、哥伦比亚锦葵、榴莲、长蒴黄麻的亲缘关系较近,与其余6种物种的亲缘关系较远,不在一个大支上。

      图  5  GhGDPD1蛋白系统进化树分析

      Figure 5.  GhGDPD1 phylogenetic tree analysis

    • 以棉花的GhActin为内参基因,通过半定量RT-PCR技术对GhGDPD1基因在‘新陆早19’的4个组织中的表达情况进行分析。从图6可知:GhGDPD1高表达于根,微量表达于叶,中量表达于茎和花。利用RT-qPCR检测GhGDPD1基因在‘新陆早19’植株的根、茎、叶、花4个组织中的表达结果(图7A),结果与半定量RT-PCR分析一致。由图7B可以看出:在低磷处理4~24 h时相对表达量逐渐下降,但相对表达量均比适磷情况下高,在72 h时其相对表达量又上升,仅略低于4 h时的相对表达量。其中,4、12、72 h时低磷处理比适磷处理显著表达(P<0.05),说明GhGDPD1基因在受到低磷胁迫刺激后,会立即对低磷胁迫做出应答,且具有持续性。

      图  6  半定量RT-PCR分析GhGDPD1在棉花不同组织中的表达

      Figure 6.  Expression of GhGDPD1 in different tissues of cotton by semi-quantitative RT-PCR results

      图  7  GhGDPD1基因在不同组织(A)及低磷胁迫处理下(B)的表达模式

      Figure 7.  Expression analysis of GhMGD1 in different tissues (A) and low phosphorus stress (B)

    • 磷是植物最重要的必需营养素之一,因为它是核酸、核苷酸和磷脂等关键大分子的组成部分,并且以磷酸盐或磷酯的形式参与新陈代谢和蛋白质调节。但磷的特性也导致植物对磷肥的低效利用。施加的磷肥只有15%~25%被植物吸收,其余的则渗入环境,造成土壤退化和水体富营养化等有害后果[16]。GDPD已被鉴定为拟南芥中的多基因家族[17]。所有已鉴定的拟南芥GDPD被分为2组:具有1个GDPD结构域的A型GDPD和具有2个推定的GDPD结构域的B型GDPD。其中B型GDPD具有植物特异性,磷酸二酯酶活性相对较低。据报道,GDPD基因能水解甘油磷酸二酯,并在磷饥饿时显著诱导,表明它们在低磷胁迫中存在潜在作用[18]

      本研究以陆地棉‘新陆早19’植株的cDNA为模板克隆到1个GDPD1基因。该基因序列开放阅读框的长度为1 149 bp,共编码382个氨基酸;该基因的结构比较复杂,有5个内含子,1个GDPD_GDE5_like_1_plant结构域,命名为GhGDPD1。生物信息学分析结果显示:GhGDPD1蛋白相对分子质量是42 987.05 Da,理论等电点为5.17,是酸性蛋白,具备不稳定性与亲水性。GhGDPD1蛋白的二级结构由α-螺旋和无规则卷曲占据主体地位所组成。通过同源建模可知:GhGDPD1基因三级结构预测其含有甘油磷酸二酯磷酸二酯酶。GhGDPD1基因编码蛋白不存在信号肽和跨膜结构域。GhGDPD1蛋白存在3个潜在的N-糖基化位点和多个酸化位点。亚细胞定位结果显示:GhGDPD1基因编码蛋白在细胞膜中,预测其能够通过蛋白质的相互作用参与调控。RT-qPCR与半定量RT-PCR结果表明:GhGDPD1基因在‘新陆早19’根中的表达量最大,在叶片中表达量最小,中量表达于花和茎。在受到低磷胁迫的刺激后,GhGDPD1基因会立即对低磷胁迫做出应答。GhGDPD1基因可能在根系发挥作用,在‘新陆早19’的磷高效利用信号调控过程中具有重要的作用。

参考文献 (18)

目录

    /

    返回文章
    返回