留言板

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

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

牡丹响应高温胁迫的转录组分析及PsHSP基因表达

郝力慧 董彬 朱绍华 马进

郝力慧, 董彬, 朱绍华, 马进. 牡丹响应高温胁迫的转录组分析及PsHSP基因表达[J]. 浙江农林大学学报, 2021, 38(4): 802-811. doi: 10.11833/j.issn.2095-0756.20200529
引用本文: 郝力慧, 董彬, 朱绍华, 马进. 牡丹响应高温胁迫的转录组分析及PsHSP基因表达[J]. 浙江农林大学学报, 2021, 38(4): 802-811. doi: 10.11833/j.issn.2095-0756.20200529
HAO Lihui, DONG Bin, ZHU Shaohua, MA Jin. Transcriptome analysis and PsHSP gene expression of Paeonia suffruticosa in response to high temperature stress[J]. Journal of Zhejiang A&F University, 2021, 38(4): 802-811. doi: 10.11833/j.issn.2095-0756.20200529
Citation: HAO Lihui, DONG Bin, ZHU Shaohua, MA Jin. Transcriptome analysis and PsHSP gene expression of Paeonia suffruticosa in response to high temperature stress[J]. Journal of Zhejiang A&F University, 2021, 38(4): 802-811. doi: 10.11833/j.issn.2095-0756.20200529

牡丹响应高温胁迫的转录组分析及PsHSP基因表达

doi: 10.11833/j.issn.2095-0756.20200529
基金项目: 浙江农林大学科研发展基金资助项目(2017FR016);元成环境股份企业研发项目(201801)
详细信息
    作者简介: 郝力慧(ORCID: 0000-0002-7687-5933),从事园林植物遗传育种研究。E-mail: 2018105252002@stu.zafu.edu.cn
    通信作者: 马进(ORCID: 0000-0003-0534-5836),副教授,从事园林植物遗传育种及抗性研究。E-mail: majinzjl@163.com
  • 中图分类号: Q943

Transcriptome analysis and PsHSP gene expression of Paeonia suffruticosa in response to high temperature stress

  • 摘要:   目的  夏季高温可导致牡丹Paeonia suffruticosa生长受限,花期缩短,观赏品质降低。研究牡丹耐热基因对牡丹热胁迫的分子机制。  方法  以牡丹品种‘羽红’‘Yuhong’为材料,对高温(40 ℃)和对照(25 ℃)处理后的叶片进行转录组测序,分析响应高温表达的差异基因;同时利用实时荧光定量PCR(qRT-PCR)验证差异基因,并对热激蛋白基因HSPs进行时空表达分析。  结果  测序共获得45.97 Gb数据,与对照组(25 ℃)相比,高温处理后4 220个基因上调表达,3 453个基因下调表达。通过基因本体(GO)分析发现,这些差异基因主要富集于代谢等生物学过程、细胞和膜结构等细胞成分以及结合和催化活性等分子功能条目;另外,通过全基因组及代谢途径数据库(KEGG)分析发现,碳代谢途径的差异基因数量最多。对牡丹热激蛋白PsHSP基因定量分析发现,PsHSP基因随着高温处理时间呈上升表达趋势,24 h达到最高,之后便呈现下降趋势。  结论  高温显著影响牡丹的代谢和生物合成等过程,进而影响牡丹的生长和发育。PsHSP基因可在短期内迅速响应高温胁迫并参与牡丹的耐热调控。图 4表6参34
  • 图  1  单基因长度分布

    Figure  1  Demo unigene length distribution

    图  2  非冗余蛋白质注释分布

    Figure  2  Nr annotation distribution

    图  3  差异基因定量表达及转录组测序比较

    Figure  3  Real-time fluorescence quantitative results for differential genes

    图  4  高温处理后不同时间HSP的qRT-PCR表达变化

    Figure  4  qRT-PCR detection of HSP in different time periods after high temperature treatment

    表  1  实时荧光定量基因引物信息

    Table  1.   Real-time fluorescence quantitative gene primer information

    基因编号基因名称引物序列(5′→3′)退火温度/℃
    C65824.graph_c0PsHSP70GCAGCCGTTCAGGCAGCAA58.3
    ACACCGCCAGCAGTTTCG
    C114411.graph_c0Ps4CL1AGGCTATGCTACTCACGC53.5
    TCATCCTCGGCAATCTG
    C106006.graph_c0PsCAT5GAGTGTTGCGTTAATAGGCG56.5
    GCACTTGGACAGGGAAGGTAT
    C111201.graph_c0PsCHL1GAGTGCTGAGATTTTAAGTGGC56.7
    CAACCCTTTCTACCTGACGA
    C120725.graph_c0PsAMP1ATTGTACGGCAGAACAGGC55.3
    GTGTAGCTCAGAACCCATTAGAC
    C112526.graph_c0PsSAMDCGTCCAGGCGAATGTTCTAATG56.7
    GTCGAAGGTCTTCAATGCTG
    C119988.graph_c0PsRPS4GGCTTTCATCCGATTTGTTG55.4
    ATGCCACCTGTCCTTACCTG
    C97637.graph_c0PsCCD4GTGAGGCATTTGCTTTCCG55.3
    TGCTGGGTTCATTCCTATCTG
    C111758.graph_c0PsPPIL1TCCCCTCCAACTATTCCTCC52.2
    AACTCAACCCCTTCCGATT
    C120729.graph_c1Hypothetical proteinCACATGAAGAGGACGACCAA56.5
    CCACATACATCAACCCAGGAC
    UbiquitinUbiquitinGACCTATACCAAGCCGAAG55.1
    CGTTCCAGCACCACAATC
    下载: 导出CSV

    表  2  样品测序数据评估统计表

    Table  2.   Statistical table of sample sequencing data evaluation

    样品名称读取碱基数/个总碱基数/个GC碱基数
    比例/%
    Q30碱基
    百分比/%
    ck1 24 156 058 7 200 895 95844.9696.78
    ck2 24 256 024 7 236 093 26644.7696.89
    ck3 25 256 661 7 530 655 89244.5396.74
    TE1 30 430 627 9 077 213 03243.6896.79
    TE2 22 967 108 6 856 944 37843.9596.73
    TE3 27 110 006 8 071 678 01043.6896.71
    合计154 176 48445 973 480 536
    下载: 导出CSV

    表  3  单基因簇的注释统计

    Table  3.   Unigene annotation statistics

    数据库注释300≤长度<
    1 000 nt
    长度≥
    1 000 nt
    数量百分比/%
    COG6 91528.441 6985 217
    GO14 56859.935 0119 557
    KEGG8 41634.622 8135 603
    KOG13 66956.234 7618 908
    Pfam15 10862.154 09611 012
    Swissprot15 42663.464 82510 601
    Eggnog21 79889.677 89113 907
    Nr24 00998.769 50914 500
    总数24 309100.009 73114 578
    下载: 导出CSV

    表  4  差异表达基因本体注释

    Table  4.   Difference expression gene GO annotation statistics

    注释本体
    类别
    编号功能类型基因注释数量比例/%注释本体
    类别
    编号功能类型基因注释数量比例/%
    生物过程1代谢过程1 70226.42细胞组分26复杂大分子复合体2423.73
    2细胞进程1 48923.1127细胞外区域630.97
    3单一生物进程1 14817.8228细胞连接540.83
    4生物调节4577.0929共质体530.82
    5定位4366.7730腔上包膜480.74
    6应激反应3665.6831超分子复合物170.26
    7细胞组分及生物合成2133.3132细胞外区域部分80.12
    8信号转导1281.9933病毒体40.06
    9发育过程1231.9134病毒体部分40.06
    10多细胞机体进程1221.8935拟核20.03
    11繁殖741.15分子功能36催化活性1 71939.95
    12繁殖进程741.1537结合1 52335.09
    13有机体进程490.1638转运活性2917.63
    14生化解毒350.5439结构分子活性581.52
    15生长150.2340核酸结合转录因子活性461.21
    16免疫系统进程90.1441分子功能调节器461.21
    17细胞活动20.0342信号传感器活动350.92
    18节律进程10.0243抗氧化活性270.71
    19生物黏附0044分子传感器活性260.68
    细胞组分20细胞1 31620.2745电子载体活动210.55
    21细胞组分1 30220.0646转录因子活性140.37
    22细胞膜1 23218.9847养分储藏活性40.10
    23细胞膜组成93514.4048蛋白质标签10.03
    24细胞器86613.3449翻译调节器活动10.03
    25细胞器成分3455.3250金属伴侣蛋白活性00
      说明:生物过程编号为1~19;细胞组分编号为20~35;分子功能编号为36~50
    下载: 导出CSV

    表  5  高温处理下差异表达基因的代谢通路分析

    Table  5.   Analysis of metabolic pathways for gene expression of differences under high temperature treatment

    全基因组及代谢途径通路术语基因数目/个全基因组及代谢途径编号校正后P
    苯丙烷生物合成47ko009401.00E−07
    内质网蛋白质加工59ko041415.50E−02
    脂肪酸延伸率14ko000623.41E−05
    萜类骨架的生物合成26ko009006.79E−05
    光合生物中的碳固定38ko007107.39E−05
    氮代谢15ko009107.84E−05
    碳代谢99ko012001.23E−04
    光合作用25ko001953.69E−04
    类黄酮生物合成12ko009415.23E−04
    脂肪酸降解22ko000719.32E−04
    亚油酸代谢6ko005917.97E−03
    α-亚麻酸代谢18ko005929.71E−03
    鞘脂代谢12ko006009.81E−03
    丙酮酸代谢37ko006201.11E−02
    氰氨基酸代谢17ko004601.12E−02
    缬氨酸、亮氨酸和异亮氨酸的降解19ko002801.30E−02
    类固醇生物合成12ko001001.33E−02
    脂肪酸代谢25ko012121.61E−02
    谷胱甘肽代谢27ko004801.65E−02
    其他聚糖降解10ko005111.70E−02
    光合作用-天线蛋白质8ko001969.17E−03
    糖酵解/糖异生41ko000102.24E−02
    油菜素类固醇生物合成6ko009052.34E−02
    甘氨酸、丝氨酸和苏氨酸的代谢21ko002605.25E−02
    下载: 导出CSV

    表  6  牡丹叶片响应高温的差异热激蛋白基因

    Table  6.   Peony leaves respond to high temperature heat shock DEGs

    基因编号基因命名基因表达量log2A亚细胞位置预测
    对照组25 ℃高温组40 ℃
    c112621.graph_c0PsHSP12.61295.986.83内质网
    c107551.graph_c0PsHSP259.462 254.785.24细胞核
    c114303.graph_c1PsHSP3161.17440.531.45线粒体
    c112152.graph_c0PsHSP411.64162.013.80线粒体
    c117526.graph_c2PsHSP57.06172.154.61细胞核
    c103807.graph_c0PsHSP610.7926.401.29线粒体
    c116391.graph_c0PsHSP7150.152 659.724.15叶绿体/细胞质
    c100675.graph_c0PsHSP821.49215.133.32叶绿体/细胞质
    c117374.graph_c1PsHSP9100.54591.552.88叶绿体
    c114184.graph_c0PsHSP1011.60136.373.56叶绿体/细胞质
    c119657.graph_c2PsHSP1163.383 980.175.97细胞质
    c100778.graph_c0PsHSP12102.14614.302.59叶绿体/线粒体
    c120033.graph_c1PsHSP13127.21608.072.26细胞核
    c114516.graph_c1PsHSP145.7928.502.30细胞膜
    c118563.graph_c0PsHSP159.9358.462.56叶绿体
    c111652.graph_c0PsHSP1615.3434.671.18细胞核
    c65824.graph_c0PsHSP1733.61873.894.19线粒体
    下载: 导出CSV
  • [1] 郭建平. 气候变化对中国农业生产的影响研究进展[J]. 应用气象学报, 2015, 26(1): 1 − 11.

    GUO Jianping. Advances in impacts of climate change on agricultural production in China [J]. J Appl Meteorol Sci, 2015, 26(1): 1 − 11.
    [2] KOTAK S, LARKINDALE J, LEE U, et al. Complexity of the heat stress response in plants [J]. Curr Opin Plant Biol, 2007, 10(3): 310 − 316. doi:  10.1016/j.pbi.2007.04.011
    [3] 许蓝心, 田如男. 观赏植物耐热性研究进展[J]. 亚热带植物科学, 2019, 48(1): 96 − 102.

    XU Lanxin, TIAN Runan. Advances in research on ornamental plant heat tolerance [J]. Subtrop Plant Sci, 2019, 48(1): 96 − 102.
    [4] 郝召君, 周春华, 刘定, 等. 高温胁迫对芍药光合作用、叶绿素荧光特性及超微结构的影响[J]. 分子植物育种, 2017, 15(6): 2359 − 2367.

    HAO Zhaojun, ZHOU Chunhua, LIU Ding, et al. Effects of high temperature stress on photosynthesi, chlorophyll fluorescence and ultrastructure of herbaceous peony (Paeonia lactiflora Pall.) [J]. Mol Plant Breed, 2017, 15(6): 2359 − 2367.
    [5] 骆俊, 韩金蓉, 王艳, 等. 高温胁迫下牡丹的抗逆生理响应[J]. 长江大学学报(自然科学版), 2011, 8(2): 223 − 226, 287 − 288.

    LUO Jun, HAN Jinrong, WANG Yan, et al. Response of heat stress on the physiological biochemistry of Paeonia suffruticosa [J]. J Yangtze Univ Nat Sci Ed, 2011, 8(2): 223 − 226, 287 − 288.
    [6] 彭勇政, 刘智媛, 朱晓非, 等. 5个月季品种高温处理后生理指标变化及其耐热性评价[J]. 上海交通大学学报(农业科学版), 2019, 37(5): 53 − 58.

    PENG Yongzheng, LIU Zhiyuan, ZHU Xiaofei, et al. Physiological index changes and heat tolerance evaluation of five rose cultivars after high temperature treatment [J]. J Shanghai Jiaotong Univ Agric Sci, 2019, 37(5): 53 − 58.
    [7] 李永华, 孙丽丹, 苏志国, 等. 短期高温对牡丹叶片光合作用及相关生理指标的影响[J]. 河南科学, 2008, 26(3): 291 − 293.

    LI Yonghua, SUN Lidan, SU Zhiguo, et al. Effects of short-term high temperature on photosynthesis and related physiological indices in the leaves of Paeonia suffruticosa [J]. Henan Sci, 2008, 26(3): 291 − 293.
    [8] 骞光耀, 孔祥生, 张淑玲. 3个牡丹品种对高温胁迫的生理响应[J]. 江苏农业科学, 2017, 45(12): 103 − 105.

    QIAN Guangyao, KONG Xiangsheng, ZHANG Shuling. Physiological response of three Paeonia suffruticosa varieties to high temperature stress [J]. Jiangsu Agric Sci, 2017, 45(12): 103 − 105.
    [9] 贾琪, 吴名耀, 梁康迳, 等. 基因组学在作物抗逆性研究中的新进展[J]. 中国生态农业学报, 2014, 22(4): 375 − 385.

    JIA Qi, WU Mingyao, LIANG Kangjing, et al. Advances in applications of genomics in stress resistance studies of crops [J]. Chin J Eco-Agric, 2014, 22(4): 375 − 385.
    [10] 祝小云. 切花非洲菊苗期在高温胁迫中的生理生化响应和转录组分析[D]. 杭州: 浙江农林大学, 2016.

    ZHU Xiaoyun. Physiological and Biochemical Responses of Gerbera Cultivars to Heat Stress and Transcriptome Analysis[D]. Hangzhou: Zhejiang A&F University, 2016.
    [11] LI Yongfang, WANG Yixing, TANG Yuhong, et al. Transcriptome analysis of heat stress response in switchgrass (Panicum virgatum L.) [J]. BMC Plant Biol, 2013, 13(1): 153 − 164. doi:  10.1186/1471-2229-13-153
    [12] WANG Wangxia, VINOCUR B, SHOSEYOY O, et al. Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response [J]. Trends Plant Sci, 2004, 9(5): 244 − 252. doi:  10.1016/j.tplants.2004.03.006
    [13] OHAMA N, SATO H, SHINOZAKI K, et al. Transcriptional regulatory network of plant heat stress response [J]. Trends Plant Sci, 2017, 22(1): 53 − 65. doi:  10.1016/j.tplants.2016.08.015
    [14] BUCHANAN B B, GRUISSEM W, JONES R L. 等. 植物生物化学与分子生物学 (2004)[M]. 瞿礼嘉, 顾红雅, 白书农, 等译. 北京: 科学出版社, 2004.
    [15] LIU Liang, CHEN Jiyun, YANG Bo, et al. Crystal structure and function of an unusual dimericHsp20.1 provide insight into the thermal protection mechanism of small heat shock proteins [J]. Biochem Biophys Res Commun, 2015, 458(2): 429 − 434. doi:  10.1016/j.bbrc.2015.01.134
    [16] HU Xiuli, YANG Yanfang, GONG Fangping, et al. Protein sHSP26 improves chloroplast performance under heat stress by interacting with specific chloroplast proteins in maize (Zea mays) [J]. J Proteomics, 2014, 115: 81 − 92.
    [17] GAO Lulu, MA Yuzhu, WANG Peng, et al. Transcriptome profiling of clematis apiifolia: insights into heat-stress responses [J]. DNA Cell Biol, 2017, 36(11). doi:  10.1089/dna.2017.3850
    [18] 王莲英. 中国牡丹品种图志[M]. 北京: 中国林业出版社, 1997: 2−28.
    [19] 史倩倩. 基于转录组测序滇牡丹花色形成分子调控机理研究[D]. 北京: 中国林业科学研究院, 2015.

    SHI Qianqian. Studies on the Molecular Mechanism of Paeonia delavayi Flower Color Formation[D]. Beijing: Chinese Academy of Forestry, 2015.
    [20] 张庆雨, 于蕊, 谢力行, 等. 牡丹种子脂肪酸合成相关基因PrLPAAT4的克隆与表达分析[C]//中国园艺学会观赏园艺专业委员会、国家花卉工程技术研究中心. 中国观赏园艺研究进展(2018). 北京: 中国林业出版社, 2018: 502 − 511.
    [21] 郑艳伟. 江南牡丹品种资源调查与引种栽培研究[D]. 杭州: 浙江农林大学, 2009.

    ZHENG Yanwei. The Study on the Investigation of Tree Peony Varieties and Its Introduction Cultivation in the South Yangtse River of China[D]. Hangzhou: Zhejiang A&F University, 2009.
    [22] 张佳平, 李丹青, 聂晶晶, 等. 高温胁迫下芍药的生理生化响应和耐热性评价[J]. 核农学报, 2016, 30(9): 1848 − 1856.

    ZHANG Jiaping, LI Danqing, NIE Jingjng, et al. Physiological and biochemical responses to the high temperature stress and heat resistance evaluation of Paeonia lactiflora Pall. cultivars [J]. J Nucl Agric Sci, 2016, 30(9): 1848 − 1856.
    [23] MALONE J H, OLIVER B C. Microarrays, deep sequencing and the true measure of the transcriptome[J]. BMC Biol, 2011, 9(1): 34. doi: 10.1186/1741-7007-9-34.
    [24] PATEL R K, MUKESH J. NGS QC Toolkit: a toolkit for quality control of next generation sequencing data[J]. PLoS One, 2012, 7(2): e30619. doi: 10.1371/journal.pone.0030619.
    [25] TRAPNELL C, HENDRIEKSON D G, SAUVAGEAU M, et al. Differential analysis of gene regulation at transcript resolution with RNA-seq [J]. Nat Biotechnol, 2013, 3l (1): 46 − 53.
    [26] LIVAK K J, LIVAK SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method [J]. Methods, 2001, 25(4): 402 − 408. doi:  10.1006/meth.2001.1262
    [27] ARYA M, SHERGILL I S, WILLIAMSON M, et al. Basic principles of real-time quantitative PCR [J]. Expert Rev Mol Diagn, 2005, 5 (2): 209 − 219. doi:  10.1586/14737159.5.2.209
    [28] 马月萍, 戴思兰, 马艳蓉. 荧光定量PCR技术在植物研究中的应用[J]. 生物技术通报, 2011(7): 37 − 45.

    MA Yueping, DAI Silan, MA Yanrong. Application of technique of quantitative real-time PCR in research of lants [J]. Biotechnol Bull, 2011(7): 37 − 45.
    [29] 王彦杰, 董丽, 张超, 等. 牡丹实时定量PCR分析中内参基因的选择[J]. 农业生物技术学报, 2012, 20(5): 521 − 528.

    WANG Yanjie, DONG Li, ZHANG Chao, et al. Reference gene selection for real-time quantitative PCR normalization in tree peony (Paeonia suffruticosa Andr.) [J]. J Agric Biotechnol, 2012, 20(5): 521 − 528.
    [30] WAHID A, GELANI S, ASHRAF M, et al. Heat tolerance in plants: an overview [J]. Environ Exp Bot, 2007, 61(3): 199 − 223. doi:  10.1016/j.envexpbot.2007.05.011
    [31] YIN Hui, CHEN Qiuming, YI Mingfang. Effects of short-term heat stress on oxidative damage and responses of antioxidant system in Lilium longiflorum [J]. Plant Growth Regul, 2008, 54(1): 45 − 54.
    [32] ZHAO Daqiu, XIA Xing, SU Jianghong, et al. Overexpression of herbaceous peony HSP70 confers high temperature tolerance[J]. BMC Genomics, 2019, 20(1): 70. doi: 10.1186/s12864-019-5448-0.
    [33] 宫本贺. 百合热激转录因子L1HSFA1及其下游热激蛋白L1HSP70响应热胁迫的机制解析[D]. 北京: 中国农业大学, 2014.

    GONG Benhe. Mechanism Analysis of Response to Heat Stress of LIHSFA1 and its Downstream LIHSP70 from Lily(Lilium longiflorum)[D]. Beijing: China Agricultural University, 2014.
    [34] SOO M P, CHOO B H. Class Ⅰ small heat-shock protein gives thermotolerance in tobacco [J]. J Plant Physiol, 2002, 159(1): 25 − 30. doi:  10.1078/0176-1617-00660
  • [1] 李叶华, 陈爽, 赵冰.  铝胁迫对绣球组培苗生长及生理特性的影响 . 浙江农林大学学报, 2020, 37(6): 1064-1070. doi: 10.11833/j.issn.2095-0756.20200121
    [2] 王颖颖, 赵冰, 李莹.  丛枝菌根真菌对杜鹃花耐热性的影响 . 浙江农林大学学报, 2019, 36(4): 733-740. doi: 10.11833/j.issn.2095-0756.2019.04.013
    [3] 许蔷薇, 楼雄珍, 杨彬, 林二培, 童再康.  云锦杜鹃转录组分析 . 浙江农林大学学报, 2019, 36(6): 1190-1198. doi: 10.11833/j.issn.2095-0756.2019.06.017
    [4] 尹跃, 安巍, 赵建华, 王亚军, 樊云芳, 曹有龙.  黑果枸杞转录组SSR信息分析及分子标记开发 . 浙江农林大学学报, 2019, 36(2): 422-428. doi: 10.11833/j.issn.2095-0756.2019.02.025
    [5] 黄耀辉, 张超, 周莉花, 赵宏波.  基于转录组序列的夏蜡梅SSR位点特征与引物开发 . 浙江农林大学学报, 2017, 34(4): 589-596. doi: 10.11833/j.issn.2095-0756.2017.04.004
    [6] 安苗苗, 刘静, 郦元, 周明兵.  花叶矢竹转录组中的转座子表达分析 . 浙江农林大学学报, 2016, 33(6): 935-943. doi: 10.11833/j.issn.2095-0756.2016.06.003
    [7] 贾思振, 杨恒伟, 颜志明, 魏跃.  高温胁迫下水杨酸对菊花幼苗生理生化指标的影响 . 浙江农林大学学报, 2016, 33(3): 449-454. doi: 10.11833/j.issn.2095-0756.2016.03.011
    [8] 龚仲幸, 何勇, 朱祝军.  水杨酸对高温胁迫下一串红耐热性的影响 . 浙江农林大学学报, 2015, 32(5): 701-707. doi: 10.11833/j.issn.2095-0756.2015.05.007
    [9] 姜晓玲, 黄秋娴, 李虹, 赵嘉平.  植物类甜蛋白基因家族研究进展 . 浙江农林大学学报, 2012, 29(2): 279-287. doi: 10.11833/j.issn.2095-0756.2012.02.019
    [10] 何勇清, 方佳, 余敏芬, 方仲相, 江波, 潘寅辉, 郑炳松.  植物质膜内在水通道蛋白PIPs的分子生物学研究进展 . 浙江农林大学学报, 2012, 29(3): 446-452. doi: 10.11833/j.issn.2095-0756.2012.03.020
    [11] 徐小博, 贾文庆, 刘会超.  牡丹品种间花瓣营养成分差异 . 浙江农林大学学报, 2012, 29(5): 729-733. doi: 10.11833/j.issn.2095-0756.2012.05.015
    [12] 李宗艳, 肖娟, 蒙进芳, 王锦.  丽江牡丹和中原牡丹叶片结构微形态比较 . 浙江农林大学学报, 2011, 28(1): 115-120. doi: 10.11833/j.issn.2095-0756.2011.01.018
    [13] 马进, 刘志高, 郑钢.  差异蛋白质组学及其在植物盐胁迫响应研究中的应用 . 浙江农林大学学报, 2011, 28(1): 139-143. doi: 10.11833/j.issn.2095-0756.2011.01.022
    [14] 孔红, 成仿云.  滇牡丹分类处理的细胞学与分子生物学证据 . 浙江农林大学学报, 2010, 27(4): 601-605. doi: 10.11833/j.issn.2095-0756.2010.04.021
    [15] 朱小楼, 楼炉焕, 王慧.  土壤干旱胁迫对4种薹草植物生理生化特性的影响 . 浙江农林大学学报, 2009, 26(5): 656-662.
    [16] 郑艳伟, 范义荣, 郭晨瑛, 范庆君, 张庆宝.  浙皖两地牡丹的栽培及应用现状 . 浙江农林大学学报, 2009, 26(6): 835-841.
    [17] 朱澜, 李雪芹, 贾晓琳, 王斌, 金松恒.  高温胁迫对高羊茅光合作用的影响 . 浙江农林大学学报, 2009, 26(5): 652-655.
    [18] 曹福亮, 欧祖兰.  水杨酸对银杏幼苗抗高温胁迫能力的影响 . 浙江农林大学学报, 2008, 25(6): 756-759.
    [19] 向佐湘, 许桂芳, 蒋文君.  干旱胁迫对4种刺篱植物抗性生理生化指标的影响 . 浙江农林大学学报, 2007, 24(1): 7-11.
    [20] 李纪元, 李辛雷, 范妙华, 田敏, 范正琪.  高温胁迫下15 个茶花品种的耐热性 . 浙江农林大学学报, 2006, 23(6): 636-640.
  • 加载中
  • 链接本文:

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

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2021/4/802

图(4) / 表(6)
计量
  • 文章访问数:  138
  • HTML全文浏览量:  22
  • PDF下载量:  13
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-08-18
  • 修回日期:  2021-03-11
  • 网络出版日期:  2021-08-09
  • 刊出日期:  2021-08-20

牡丹响应高温胁迫的转录组分析及PsHSP基因表达

doi: 10.11833/j.issn.2095-0756.20200529
    基金项目:  浙江农林大学科研发展基金资助项目(2017FR016);元成环境股份企业研发项目(201801)
    作者简介:

    郝力慧(ORCID: 0000-0002-7687-5933),从事园林植物遗传育种研究。E-mail: 2018105252002@stu.zafu.edu.cn

    通信作者: 马进(ORCID: 0000-0003-0534-5836),副教授,从事园林植物遗传育种及抗性研究。E-mail: majinzjl@163.com
  • 中图分类号: Q943

摘要:   目的  夏季高温可导致牡丹Paeonia suffruticosa生长受限,花期缩短,观赏品质降低。研究牡丹耐热基因对牡丹热胁迫的分子机制。  方法  以牡丹品种‘羽红’‘Yuhong’为材料,对高温(40 ℃)和对照(25 ℃)处理后的叶片进行转录组测序,分析响应高温表达的差异基因;同时利用实时荧光定量PCR(qRT-PCR)验证差异基因,并对热激蛋白基因HSPs进行时空表达分析。  结果  测序共获得45.97 Gb数据,与对照组(25 ℃)相比,高温处理后4 220个基因上调表达,3 453个基因下调表达。通过基因本体(GO)分析发现,这些差异基因主要富集于代谢等生物学过程、细胞和膜结构等细胞成分以及结合和催化活性等分子功能条目;另外,通过全基因组及代谢途径数据库(KEGG)分析发现,碳代谢途径的差异基因数量最多。对牡丹热激蛋白PsHSP基因定量分析发现,PsHSP基因随着高温处理时间呈上升表达趋势,24 h达到最高,之后便呈现下降趋势。  结论  高温显著影响牡丹的代谢和生物合成等过程,进而影响牡丹的生长和发育。PsHSP基因可在短期内迅速响应高温胁迫并参与牡丹的耐热调控。图 4表6参34

English Abstract

郝力慧, 董彬, 朱绍华, 马进. 牡丹响应高温胁迫的转录组分析及PsHSP基因表达[J]. 浙江农林大学学报, 2021, 38(4): 802-811. doi: 10.11833/j.issn.2095-0756.20200529
引用本文: 郝力慧, 董彬, 朱绍华, 马进. 牡丹响应高温胁迫的转录组分析及PsHSP基因表达[J]. 浙江农林大学学报, 2021, 38(4): 802-811. doi: 10.11833/j.issn.2095-0756.20200529
HAO Lihui, DONG Bin, ZHU Shaohua, MA Jin. Transcriptome analysis and PsHSP gene expression of Paeonia suffruticosa in response to high temperature stress[J]. Journal of Zhejiang A&F University, 2021, 38(4): 802-811. doi: 10.11833/j.issn.2095-0756.20200529
Citation: HAO Lihui, DONG Bin, ZHU Shaohua, MA Jin. Transcriptome analysis and PsHSP gene expression of Paeonia suffruticosa in response to high temperature stress[J]. Journal of Zhejiang A&F University, 2021, 38(4): 802-811. doi: 10.11833/j.issn.2095-0756.20200529
  • 全球变暖导致植物在生长发育和环境适应性等方面面临更加多样化的挑战[1]。高温不仅影响植物表型变化,而且会引起植物体内细胞稳态失衡,生长和发育受到抑制,严重影响观赏植物的品质和产量[2-3]。如芍药Paeonia lactiflora在遭受高温胁迫后,植株出现萎蔫、畸形、变色等热害现象,甚至全株死亡[4]。在牡丹P. suffruticosa中也出现类似的表型现象[5]。另外,经过高温处理后,牡丹和月季Rosa chinensis植株内的脯氨酸、丙二醛、可溶性糖以及相对细胞膜透性显著升高,而净光合速率、超氧化物歧化酶和可溶性蛋白质含量等显著降低[6-8]。为了深入了解高温胁迫对观赏植物品质的影响,利用现代分子生物学以及转录组学等手段[9],挖掘关键基因并解析其调控耐热的分子机制,是培育和开发耐热优良品种的重要途径。对非洲菊Gerbera jamesonii切花高温处理后,转录组分析发现大部分单基因簇(unigene)可注释到细胞过程、代谢过程和结合功能等基因本体(GO)过程中[10];全基因组及代谢途径数据库(KEGG)中代谢功能途径中的单基因簇分布最多。在柳枝稷Panicum virgatum中发现[11]:长期高温胁迫条件下大量的基因发生差异表达,其中与ATP酶调节因子、伴侣结合和蛋白质折叠相关的基因显著上调,与蛋白质修饰、转录、磷和氮代谢途径相关的基因在热胁迫后显著下调。热激蛋白(heat shock protein, HSPs)是一类高温胁迫诱导的主要功能蛋白,其分子量为10~200 kDa,包括HSP60、HSP70、HSP90、HSP100和小分子量热激蛋白[12](small heat shock protein,sHSP)等5类。HSPs充当蛋白质质量控制的分子伴侣,作为分子伴侣和折叠催化剂细胞网络的中心成分,在热应激中被诱导并参与热应激过程中的信号转导[13];sHSP存在多数植物叶肉细胞中,在外界热诱导条件下的sHSP表达呈现上升的趋势[14-16]。在女萎Clematis apiifolia[17]中,高温处理后热激蛋白热胁迫后明显上调,而且sHSPs是对热胁迫反应最强烈的伴侣基因家族。HSPs不仅响应高温变化,还对水分、盐渗透和氧化胁迫等[12, 14]有响应。牡丹是芍药科Paeoniaceae芍药属Paeonia多年生落叶灌木,中国传统名花[18]。目前,关于牡丹基因的研究多集中于生长发育、花色基因挖掘以及抗氧化活性物质生物合成等方面[19-22],相关耐热研究报道较少。牡丹对高温较为敏感,夏季高温严重制约牡丹的生长和发育,致使花期缩短,观赏品质降低,因此,挖掘牡丹耐热基因对牡丹热胁迫分子机制探究,以及牡丹育种和推广具有重要的意义。

    • 实验材料牡丹‘羽红’P. suffruticosa ‘Yuhong’保存在浙江农林大学智能实验大棚。取长势一致的1年生幼苗移栽至花盆缓苗,1株·盆−1,正常管理。后转移至人工气候箱中,在温度25 ℃,光12 h/暗12 h,光强4 000 lx,湿度80%条件下预培养2周。处理组和对照组分别置于40和25 ℃人工气候箱中处理24 h,选取相同位置的叶片,使用液氮速冻并保存于−80 ℃冰箱于后期备用,设置3个生物学重复。

    • 使用Trizol试剂盒(天根,北京)提取牡丹叶片总RNA,按照说明书操作。用紫外分光光度计和质量分数1.5%琼脂糖凝胶电泳检查RNA质量,质量检测合格后,经百迈客生物科技有限公司(北京)使用Illumina HiSeq2500对6个独立样品进行双端文库构建以及RNA-Seq测序。

    • 将测序获得的原始序列读长(raw data)进行过滤组装,得到高质量的有效序列读长(clean reads),然后使用Trinity进行序列组装[23-24]。所得单基因簇序列用基于BLAST算法的搜索和注释,对照非冗余蛋白质(Nr)数据库(www.ncbi.nlm.nih.gov/)、全基因组及代谢途径(KEGG)数据库(www.genome.jp/kegg/kaas/)、蛋白直系同源簇(COG)数据库(www.ncbi.nlm.nih.gov/COG/)以及瑞士蛋白质序列及注解(Swissprot)数据库(www.expasy.ch/SPORT)等数据库进行注释。使用Blast2Go软件(www.blast2go.com/)进行基因本体论(GO)功能分析。

    • 使用Bowtie软件将测序所得的片段与非冗余基因序列库进行对比,结合RSEM(RNA-Seq by Expectation-Maximization)预估表达量水平[24-25]。在此基础上分析高温胁迫处理前后6组样本的差异表达基因,使用RPKM[每百万序列读长(Reads)中某基因每千长度序列读长值]统计基因的表达量;单基因簇表达丰度用FPKM值表示,以差异倍数|log2A|≥2。其中,A表示差异倍数或变化倍数(flod change),且错误发现率(false discovery rate)≤0.01为标准筛选差异表达基因。

    • 使用Primer Premier 5.0软件设计荧光定量PCR引物(表1),Prime Script RT reagent Kit With gDNA Eraser试剂盒(TaKaRa,北京)进行实时荧光定量PCR反应。反应体系:上下游引物(10 μmol·L−1)各0.4 μL,cDNA 1.0 μL,ddH2O 3.2 μL,SYBR Premix Ex Taq 5.0 μL。总反应程序:95 ℃预变性30 s,95 ℃变性5 s,60 ℃复性30 s,共38个循环。以$ {2^{ - \Delta \Delta {C_{\rm{t}}}}} $[26-28]计算目的基因相对表达量,其中以牡丹泛素延伸蛋白基因(Ubiquitin)为内参基因[29],设3次生物学重复。

      表 1  实时荧光定量基因引物信息

      Table 1.  Real-time fluorescence quantitative gene primer information

      基因编号基因名称引物序列(5′→3′)退火温度/℃
      C65824.graph_c0PsHSP70GCAGCCGTTCAGGCAGCAA58.3
      ACACCGCCAGCAGTTTCG
      C114411.graph_c0Ps4CL1AGGCTATGCTACTCACGC53.5
      TCATCCTCGGCAATCTG
      C106006.graph_c0PsCAT5GAGTGTTGCGTTAATAGGCG56.5
      GCACTTGGACAGGGAAGGTAT
      C111201.graph_c0PsCHL1GAGTGCTGAGATTTTAAGTGGC56.7
      CAACCCTTTCTACCTGACGA
      C120725.graph_c0PsAMP1ATTGTACGGCAGAACAGGC55.3
      GTGTAGCTCAGAACCCATTAGAC
      C112526.graph_c0PsSAMDCGTCCAGGCGAATGTTCTAATG56.7
      GTCGAAGGTCTTCAATGCTG
      C119988.graph_c0PsRPS4GGCTTTCATCCGATTTGTTG55.4
      ATGCCACCTGTCCTTACCTG
      C97637.graph_c0PsCCD4GTGAGGCATTTGCTTTCCG55.3
      TGCTGGGTTCATTCCTATCTG
      C111758.graph_c0PsPPIL1TCCCCTCCAACTATTCCTCC52.2
      AACTCAACCCCTTCCGATT
      C120729.graph_c1Hypothetical proteinCACATGAAGAGGACGACCAA56.5
      CCACATACATCAACCCAGGAC
      UbiquitinUbiquitinGACCTATACCAAGCCGAAG55.1
      CGTTCCAGCACCACAATC
    • 对照组(ck1、ck2和ck3)与高温(40 ℃)处理组(TE1、TE2和TE3)样本经过测序后,获得45.97 Gb的有效序列读长,且6个样品Q30碱基百分比值均大于等于96.71%(表2)。转录本序列数量(transcript number)总计151 939条,非冗余基因序列数量(unigene number)54 935条,有较高的组装完整性。经拼接后片段重叠群N50长度为1 561 nt,平均长度为1 042 nt,长度为300~500 nt的单基因簇数量占组装总量的36.1%,有19 834条。长度500~1000 nt的数量次之,占总量的28.94%。大于2 000 nt的占总量的13.33%(图1)。

      表 2  样品测序数据评估统计表

      Table 2.  Statistical table of sample sequencing data evaluation

      样品名称读取碱基数/个总碱基数/个GC碱基数
      比例/%
      Q30碱基
      百分比/%
      ck1 24 156 058 7 200 895 95844.9696.78
      ck2 24 256 024 7 236 093 26644.7696.89
      ck3 25 256 661 7 530 655 89244.5396.74
      TE1 30 430 627 9 077 213 03243.6896.79
      TE2 22 967 108 6 856 944 37843.9596.73
      TE3 27 110 006 8 071 678 01043.6896.71
      合计154 176 48445 973 480 536

      图  1  单基因长度分布

      Figure 1.  Demo unigene length distribution

    • 将单基因簇基因序列分别在各数据库注释,其中以非冗余蛋白质数据库(Nr)获得注释的比例最高,达到98.76%(表3)。

      表 3  单基因簇的注释统计

      Table 3.  Unigene annotation statistics

      数据库注释300≤长度<
      1 000 nt
      长度≥
      1 000 nt
      数量百分比/%
      COG6 91528.441 6985 217
      GO14 56859.935 0119 557
      KEGG8 41634.622 8135 603
      KOG13 66956.234 7618 908
      Pfam15 10862.154 09611 012
      Swissprot15 42663.464 82510 601
      Eggnog21 79889.677 89113 907
      Nr24 00998.769 50914 500
      总数24 309100.009 73114 578

      根据非冗余蛋白质数据库中的物种注释对比,做出对应注释图和E值分布图,大部分序列(29.13%)E值为0~10−150 (图2A),25.78%的序列相似度大于80.00%(图2B);另外,22.43%的单基因簇基因可注释到葡萄Vitis vinifera,6.56%的单基因簇序列可以注释到欧洲栓皮栎Quercus suber(图2C)。

      图  2  非冗余蛋白质注释分布

      Figure 2.  Nr annotation distribution

    • 以差异倍数≥2和错误发现率≤0.01为标准,共筛选获得7 673个差异基因,其中4 220个基因上调表达,3 453个基因下调表达。利用Web基因本体论注释软件进行基因本体富集分析(表4),有6 443个差异基因富集于生物过程,主要集中于代谢过程(1 702个,占26.42%,GO:0008152)、细胞进程(1 489个,占23.11%,GO:0009987)和单一生物进程(1 148个,占17.82%,GO:0044699);有6 491个差异基因富集于细胞组分,主要集中于细胞(1 316个,占20.27%,GO:0005623)、细胞组分(1 302个,占20.06%,GO:0044464)和细胞膜(1 232个,占18.98%,GO:0016020);有3 812个差异基因富集于分子功能,主要分布在催化活性(1 719个,占39.95%,GO:0003824)、结合活性(1 523个,占35.09%,GO:0005488)和转运活性(291个,占7.63%,GO:0005215)等过程。

      表 4  差异表达基因本体注释

      Table 4.  Difference expression gene GO annotation statistics

      注释本体
      类别
      编号功能类型基因注释数量比例/%注释本体
      类别
      编号功能类型基因注释数量比例/%
      生物过程1代谢过程1 70226.42细胞组分26复杂大分子复合体2423.73
      2细胞进程1 48923.1127细胞外区域630.97
      3单一生物进程1 14817.8228细胞连接540.83
      4生物调节4577.0929共质体530.82
      5定位4366.7730腔上包膜480.74
      6应激反应3665.6831超分子复合物170.26
      7细胞组分及生物合成2133.3132细胞外区域部分80.12
      8信号转导1281.9933病毒体40.06
      9发育过程1231.9134病毒体部分40.06
      10多细胞机体进程1221.8935拟核20.03
      11繁殖741.15分子功能36催化活性1 71939.95
      12繁殖进程741.1537结合1 52335.09
      13有机体进程490.1638转运活性2917.63
      14生化解毒350.5439结构分子活性581.52
      15生长150.2340核酸结合转录因子活性461.21
      16免疫系统进程90.1441分子功能调节器461.21
      17细胞活动20.0342信号传感器活动350.92
      18节律进程10.0243抗氧化活性270.71
      19生物黏附0044分子传感器活性260.68
      细胞组分20细胞1 31620.2745电子载体活动210.55
      21细胞组分1 30220.0646转录因子活性140.37
      22细胞膜1 23218.9847养分储藏活性40.10
      23细胞膜组成93514.4048蛋白质标签10.03
      24细胞器86613.3449翻译调节器活动10.03
      25细胞器成分3455.3250金属伴侣蛋白活性00
        说明:生物过程编号为1~19;细胞组分编号为20~35;分子功能编号为36~50

      将高温组(TE1、TE2、TE3)与对照组(ck1、ck2、ck3)相比,进行全基因组及代谢途径调控网络富集分析,共有949个差异基因(different expression gene set,DEGs)可注释到全基因组及代谢途径通路中(表5)。集中在碳代谢途径(carbon metabolism)和内质网蛋白质加工(protein processing in endoplasmic reticulum)的差异基因数量较多(99个和59个),其次是苯丙烷类生物合成(phenylpropanoid biosynthesis)和糖酵解/糖异生途径(glycolysis/gluconeogenesis)。表明在高温刺激后,牡丹中的多种代谢途径、生物合成途径和内质网蛋白加工途径等均发生变化,积极参与应激响应。

      表 5  高温处理下差异表达基因的代谢通路分析

      Table 5.  Analysis of metabolic pathways for gene expression of differences under high temperature treatment

      全基因组及代谢途径通路术语基因数目/个全基因组及代谢途径编号校正后P
      苯丙烷生物合成47ko009401.00E−07
      内质网蛋白质加工59ko041415.50E−02
      脂肪酸延伸率14ko000623.41E−05
      萜类骨架的生物合成26ko009006.79E−05
      光合生物中的碳固定38ko007107.39E−05
      氮代谢15ko009107.84E−05
      碳代谢99ko012001.23E−04
      光合作用25ko001953.69E−04
      类黄酮生物合成12ko009415.23E−04
      脂肪酸降解22ko000719.32E−04
      亚油酸代谢6ko005917.97E−03
      α-亚麻酸代谢18ko005929.71E−03
      鞘脂代谢12ko006009.81E−03
      丙酮酸代谢37ko006201.11E−02
      氰氨基酸代谢17ko004601.12E−02
      缬氨酸、亮氨酸和异亮氨酸的降解19ko002801.30E−02
      类固醇生物合成12ko001001.33E−02
      脂肪酸代谢25ko012121.61E−02
      谷胱甘肽代谢27ko004801.65E−02
      其他聚糖降解10ko005111.70E−02
      光合作用-天线蛋白质8ko001969.17E−03
      糖酵解/糖异生41ko000102.24E−02
      油菜素类固醇生物合成6ko009052.34E−02
      甘氨酸、丝氨酸和苏氨酸的代谢21ko002605.25E−02
    • 为了进一步验证转录组数据,随机选取10个上调和下调的差异基因,分别是热激蛋白70(PsHSP70,c65824.graph_c0)、核糖体蛋白S4(PsRPS4,c119988.graph_c0)、酸性黏多糖(PsAMP1,c120725.graph_c0)、肽基脯氨酰顺反异构酶(PsPPIL1,c111758.graph_c0)、4-香豆酸辅酶A连接酶(Ps4CL1,c114411.graph_c0)、阳离子氨基酸转运蛋白5(PsCAT5,c106006.graph_c0)、叶绿素酶(PsCHL1c111201.graph_c0)、类胡萝卜素裂解双加氧酶(PsCCD4,c97637.graph_c0)、腺苷甲硫氨酸脱羧酶样原酶(PsSAMDC,C112526.graph_c0)以及假定蛋白(c120729.graph_c1),使用qRT-PCR的方法进行验证分析(图3),两者在变化趋势和差异倍数总体相符,表明测序结果相对可信。

      图  3  差异基因定量表达及转录组测序比较

      Figure 3.  Real-time fluorescence quantitative results for differential genes

      图  4  高温处理后不同时间HSP的qRT-PCR表达变化

      Figure 4.  qRT-PCR detection of HSP in different time periods after high temperature treatment

    • 分析测序发现:17个热激蛋白基因在高温胁迫后显著上调表达(表6),这些热激蛋白多定位在细胞核、线粒体、叶绿体和细胞质等部位。选取6个PsHSP (PsHSP1、PsHSP2、PsHSP4、PsHSP5、PsHSP7、PsHSP11)分别在25和40 ℃处理下进行时空表达分析(图4)。在40 ℃高温处理下,6个PsHSP表达呈上升趋势,且全部于24 h达到最大值,之后便呈现下降趋势;25 ℃条件下,6个PsHSP表达变化相对稳定。

      表 6  牡丹叶片响应高温的差异热激蛋白基因

      Table 6.  Peony leaves respond to high temperature heat shock DEGs

      基因编号基因命名基因表达量log2A亚细胞位置预测
      对照组25 ℃高温组40 ℃
      c112621.graph_c0PsHSP12.61295.986.83内质网
      c107551.graph_c0PsHSP259.462 254.785.24细胞核
      c114303.graph_c1PsHSP3161.17440.531.45线粒体
      c112152.graph_c0PsHSP411.64162.013.80线粒体
      c117526.graph_c2PsHSP57.06172.154.61细胞核
      c103807.graph_c0PsHSP610.7926.401.29线粒体
      c116391.graph_c0PsHSP7150.152 659.724.15叶绿体/细胞质
      c100675.graph_c0PsHSP821.49215.133.32叶绿体/细胞质
      c117374.graph_c1PsHSP9100.54591.552.88叶绿体
      c114184.graph_c0PsHSP1011.60136.373.56叶绿体/细胞质
      c119657.graph_c2PsHSP1163.383 980.175.97细胞质
      c100778.graph_c0PsHSP12102.14614.302.59叶绿体/线粒体
      c120033.graph_c1PsHSP13127.21608.072.26细胞核
      c114516.graph_c1PsHSP145.7928.502.30细胞膜
      c118563.graph_c0PsHSP159.9358.462.56叶绿体
      c111652.graph_c0PsHSP1615.3434.671.18细胞核
      c65824.graph_c0PsHSP1733.61873.894.19线粒体
    • 全球变暖使得世界上许多地区农业受到不同程度的热害影响。近年来,植物耐热性研究逐渐受到关注。短期或持续高温胁迫引起植物形态和生理生化方面的改变,影响植物生长发育的同时也导致农林经济产量下降[1, 30-31]。高温是限制牡丹生长的重要因素且严重影响其观赏性。通过RNA-Seq测序发现:牡丹‘羽红’高温胁迫后共有7 673个基因差异表达,这些差异基因在基因本体注释主要集中在代谢过程、细胞过程和单一生物过程;全基因组及代谢途径数据库富集分析发现:大部分DEGs富集在碳代谢、内质网蛋白质的加工途径;同时,苯丙烷生物合成、光合生物中的碳固定、糖酵解/糖异生和丙酮酸代谢等途径均受到高温胁迫的影响。

      极端温度可以造成植物体内蛋白质功能障碍导致生长发育受限,热激蛋白HSPs参与植物细胞膜稳定及相关代谢过程[12, 30],在响应和调节高温胁迫的过程中发挥重要作用[32-34]。芍药PlHSP70基因过表达拟南芥Arabidopsis thaliana发现,PlHSP70对拟南芥高温的耐受性有正向调节作用,高温处理后转基因株系的叶绿素荧光值明显高于野生型植株,此外还具有更加完整的细胞膜、叶绿体和淀粉颗粒。百合Lilium brownie LlHSP70基因异源转化拟南芥,LlHSP70可以增强转基因植株的高温耐性;烟草Nicotiana tabacum Class Ⅰ 类热激蛋白TLHS1基因在烟草幼苗中的表达水平与耐热性有较强的相关性,TLHS1转至烟草幼苗经过高温胁迫1~4 h后,子叶开张率比转反义转基因烟草幼苗高出1倍以上。本研究发现:高温胁迫后,共有17个PsHSP基因在牡丹中上调表达,且PsHSP1、PsHSP2、PsHSP4、PsHSP5、PsHSP7和PsHSP11在高温处理后表达上调3倍以上;同时PsHSPs时空表达也发现,PsHSP随热胁迫时间增加而表达量上升,24 h后达到最大值后便逐渐下降,说明PsHSP基因可能在短期内迅速响应高温胁迫并参与牡丹的耐热调控。综上,通过转录组测序发现,高温显著影响牡丹的代谢和生物合成等过程,从而抑制牡丹的生长和发育。PsHSP基因可在短期内迅速响应高温胁迫并参与牡丹的耐热调控。本研究可为进一步探究牡丹耐高温胁迫的分子机制提供依据。

参考文献 (34)

目录

    /

    返回文章
    返回