Volume 38 Issue 4
Aug.  2021
Turn off MathJax
Article Contents
Article Navigation >  Journal of Zhejiang A&F University  > 2021  >  38(4): 802-811

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
shu

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

doi: 10.11833/j.issn.2095-0756.20200529
  • 1.

    Zhejiang Provincial Key Laboratory of Germplasm Innovation and Utilization for Garden Plants, Zhejiang A&F University, Hangzhou 311300, Zhejiang, China

  • 2.

    Key Laboratory of National Forestry and Grassland Administration on Germplasm Innovation and Utilization for Southern Garden Plants, Zhejiang A&F University, Hangzhou 311300, Zhejiang, China

  • 3.

    College of Landscape Architecture, Zhejiang A&F University, Hangzhou 311300, Zhejiang, China

  • Received Date: 2020-08-18
  • Rev Recd Date: 2021-03-11
    • Available Online: 2021-08-09
  • Publish Date: 2021-08-20
  •   Objective  High temperature in summer can lead to growth restriction, shorter flowering period and lower ornamental quality of Paeonia suffruticosa. This study aims to explore the molecular mechanism of heat stress induced by heat tolerance genes in P. suffruticosa.  Method  Using P. suffruticosa ‘Yuhong’ as material, transcriptome sequencing of leaves treated with high temperature (40 ℃) and room temperature (25 ℃) was conducted to analyze the differentially expressed genes in response to high temperature. Meanwhile, the differential genes were verified by real-time fluorescence quantitative PCR (qRT-PCR), and the spatio-temporal expression of heat shock protein gene HSPs was analyzed.  Result  A total of 45.97 Gb data were obtained by sequencing. Compared with the control group (25 ℃), 4 220 genes were up-regulated and 3 453 genes were down-regulated after high temperature treatment. Through gene ontology (GO) analysis, it was found that these differential genes were mainly concentrated in biological processes such as metabolism, cellular components such as cell and membrane structure, and molecular functional items such as binding and catalytic activity. In addition, genome-wide and metabolic pathway database (KEGG) analysis showed that the number of differential genes in carbon metabolic pathway was the largest. Quantitative analysis of heat shock protein PsHSP gene showed that the expression of PsHSP gene increased with the time of high temperature treatment, reached the peak at 24 h, and then showed a downward trend.  Conclusion  High temperature significantly affects the metabolism and synthesis of P. suffruticosa, and further affects its growth and development. PsHSP gene can rapidly respond to high temperature stress in a short time and participate in heat tolerance regulation in P. suffruticosa. [Ch, 4 fig. 6 tab. 34 ref.]
  • [1] DENG Xuan, CHEN Chunbing, DENG Jing, LIU Lianlian, LI Juan, ZHA Xingfu.  Transcriptome analysis of abscised and normal surviving young fruit of Morus alba . Journal of Zhejiang A&F University, 2023, 40(1): 45-54. doi: 10.11833/j.issn.2095-0756.20220205
    [2] WANG Jiaqi, WANG Xin, DENG Xiaomei, WU Aimin.  Analysis of candidate genes for terpene synthesis in Ormosia henryi based on metabolome and transcriptome . Journal of Zhejiang A&F University, 2023, 40(5): 970-981. doi: 10.11833/j.issn.2095-0756.20220737
    [3] LIU Yiling, LIU Qiyuan, DONG Shuai, ZHANG Min, LUO Yan, LI Zhenhua.  Analysis of molecular networks for dark germination of shallow photodormant Nicotiana tabacum seeds based on transcriptomic and proteomic data . Journal of Zhejiang A&F University, 2023, 40(2): 237-243. doi: 10.11833/j.issn.2095-0756.20220515
    [4] HUANG Yizi, QIAN Wang, QIU Shan, WANG Wenxin, HUANG Huahong, LIN Erpei.  Identification and expression analysis of AP2/ERF gene family in Betula luminifera . Journal of Zhejiang A&F University, 2022, 39(6): 1183-1193. doi: 10.11833/j.issn.2095-0756.20220331
    [5] LI Yehua, CHEN Shuang, ZHAO Bing.  Effect of aluminum stress on growth and physiological characteristics of Hydrangea tissue culture seedlings . Journal of Zhejiang A&F University, 2020, 37(6): 1064-1070. doi: 10.11833/j.issn.2095-0756.20200121
    [6] WANG Yingying, ZHAO Bing, LI Ying.  Heat resistance of Rhododendron with arbuscular mycorrhizal fungi . Journal of Zhejiang A&F University, 2019, 36(4): 733-740. doi: 10.11833/j.issn.2095-0756.2019.04.013
    [7] XU Qiangwei, LOU Xiongzhen, YANG Bin, LIN Erpei, TONG Zaikang.  Transcriptome sequencing and analysis of Rhododendron fortunei . Journal of Zhejiang A&F University, 2019, 36(6): 1190-1198. doi: 10.11833/j.issn.2095-0756.2019.06.017
    [8] YIN Yue, AN Wei, ZHAO Jianhua, WANG Yajun, FAN Yunfang, CAO Youlong.  SSR information in transcriptome and development of molecular markers in Lycium ruthenicum . Journal of Zhejiang A&F University, 2019, 36(2): 422-428. doi: 10.11833/j.issn.2095-0756.2019.02.025
    [9] JIA Sizhen, YANG Hengwei, YAN Zhiming, WEI Yue.  Physiological and biochemical indexes of exogenous salicylic acid on chrysanthemum seedlings with high temperature stress . Journal of Zhejiang A&F University, 2016, 33(3): 449-454. doi: 10.11833/j.issn.2095-0756.2016.03.011
    [10] GONG Zhongxing, HE Yong, ZHU Zhujun.  Heat resistance of Salvia splendens with salicylic acid added and high temperature stress . Journal of Zhejiang A&F University, 2015, 32(5): 701-707. doi: 10.11833/j.issn.2095-0756.2015.05.007
    [11] HE Yong-qing, FANG Jia, YU Min-fen, FANG Zhong-xiang, JIANG Bo, PAN Yin-hui, ZHENG Bing-song.  Advances in molecular biology for plasma membrane intrinsic proteins （PIPs）：a review . Journal of Zhejiang A&F University, 2012, 29(3): 446-452. doi: 10.11833/j.issn.2095-0756.2012.03.020
    [12] XU Xiao-bo, JIA Wen-qing, LIU Hui-chao.  Nutritional component differences in petal of Paeonia suffruticosa . Journal of Zhejiang A&F University, 2012, 29(5): 729-733. doi: 10.11833/j.issn.2095-0756.2012.05.015
    [13] MA Jin, LIU Zhi-gao, ZHENG Gang.  Differential proteomics and its application in plant salt stress . Journal of Zhejiang A&F University, 2011, 28(1): 139-143. doi: 10.11833/j.issn.2095-0756.2011.01.022
    [14] LI Zong-yan, XIAO Juan, MENG Jin-fang, WAN Jin.  Comparative leaf anatomy for six cultivars of Paeonia suffruticosa from Lijiang and Zhongyuan，China . Journal of Zhejiang A&F University, 2011, 28(1): 115-120. doi: 10.11833/j.issn.2095-0756.2011.01.018
    [15] KONG Hong, CHENG Fang-yun.  Cytological and molecular biological evidence for taxonomic revision of Paeonia delavayi . Journal of Zhejiang A&F University, 2010, 27(4): 601-605. doi: 10.11833/j.issn.2095-0756.2010.04.021
    [16] ZHENG Yan-wei, FAN Yi-rong, GUO Chen-ying, FAN Qing-jun, ZHANG Qing-bao.  Investigation on cultivation and application of tree peony in Zhejiang and Anhui . Journal of Zhejiang A&F University, 2009, 26(6): 835-841.
    [17] ZHU Lan, LI Xue-qin, JIA Xiao-lin, WANG Bin, JIN Song-heng.  High temperature stress with photosynthesis in Festuca arundinacea . Journal of Zhejiang A&F University, 2009, 26(5): 652-655.
    [18] CAO Fu-liang, OU Zu-lan.  Ginkgo biloba seedling tolerance to high-temperature stress using salicylic acid . Journal of Zhejiang A&F University, 2008, 25(6): 756-759.
    [19] XIANG Zuo-xiang, XUGui-fang, JIANG Wen-jun.  Drought stress on physiological and biochemical processes in four spiny plant species . Journal of Zhejiang A&F University, 2007, 24(1): 7-11.
    [20] LI Ji-yuan, LI Xin-lei, FAN Miao-hua, TIAN Min, FAN Zheng-qi.  Heat tolerance of 15 Camellia cultivars under heat stress . Journal of Zhejiang A&F University, 2006, 23(6): 636-640.
  • [1]
    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.
    [3]
    XU Lanxin, TIAN Runan. Advances in research on ornamental plant heat tolerance [J]. Subtrop Plant Sci, 2019, 48(1): 96 − 102.
    [4]
    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]
    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]
    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]
    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]
    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]
    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]
    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.
    [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.
    [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.
    [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.
    [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).
    [19]
    SHI Qianqian. Studies on the Molecular Mechanism of Paeonia delavayi Flower Color Formation[D]. Beijing: Chinese Academy of Forestry, 2015.
    [21]
    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]
    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.
    [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.
    [28]
    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]
    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.
    [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]
    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.
  • 加载中
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(4)  / Tables(6)

Get Citation
PDF
XML

Article views(1231) PDF downloads(62) Cited by()

Related
Proportional views

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

doi: 10.11833/j.issn.2095-0756.20200529
  • 1. Zhejiang Provincial Key Laboratory of Germplasm Innovation and Utilization for Garden Plants, Zhejiang A&F University, Hangzhou 311300, Zhejiang, China
  • 2. Key Laboratory of National Forestry and Grassland Administration on Germplasm Innovation and Utilization for Southern Garden Plants, Zhejiang A&F University, Hangzhou 311300, Zhejiang, China
  • 3. College of Landscape Architecture, Zhejiang A&F University, Hangzhou 311300, Zhejiang, China
  • Received Date: 2020-08-18
  • Rev Recd Date: 2021-03-11
  • Available Online: 2021-08-09
  • Publish Date: 2021-08-20

Abstract:   Objective  High temperature in summer can lead to growth restriction, shorter flowering period and lower ornamental quality of Paeonia suffruticosa. This study aims to explore the molecular mechanism of heat stress induced by heat tolerance genes in P. suffruticosa.  Method  Using P. suffruticosa ‘Yuhong’ as material, transcriptome sequencing of leaves treated with high temperature (40 ℃) and room temperature (25 ℃) was conducted to analyze the differentially expressed genes in response to high temperature. Meanwhile, the differential genes were verified by real-time fluorescence quantitative PCR (qRT-PCR), and the spatio-temporal expression of heat shock protein gene HSPs was analyzed.  Result  A total of 45.97 Gb data were obtained by sequencing. Compared with the control group (25 ℃), 4 220 genes were up-regulated and 3 453 genes were down-regulated after high temperature treatment. Through gene ontology (GO) analysis, it was found that these differential genes were mainly concentrated in biological processes such as metabolism, cellular components such as cell and membrane structure, and molecular functional items such as binding and catalytic activity. In addition, genome-wide and metabolic pathway database (KEGG) analysis showed that the number of differential genes in carbon metabolic pathway was the largest. Quantitative analysis of heat shock protein PsHSP gene showed that the expression of PsHSP gene increased with the time of high temperature treatment, reached the peak at 24 h, and then showed a downward trend.  Conclusion  High temperature significantly affects the metabolism and synthesis of P. suffruticosa, and further affects its growth and development. PsHSP gene can rapidly respond to high temperature stress in a short time and participate in heat tolerance regulation in P. suffruticosa. [Ch, 4 fig. 6 tab. 34 ref.]

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
shu

    HTML

  • 全球变暖导致植物在生长发育和环境适应性等方面面临更加多样化的挑战[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]，相关耐热研究报道较少。牡丹对高温较为敏感，夏季高温严重制约牡丹的生长和发育，致使花期缩短，观赏品质降低，因此，挖掘牡丹耐热基因对牡丹热胁迫分子机制探究，以及牡丹育种和推广具有重要的意义。

1.   材料与方法

    1.1.   植物材料

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

    • 1.2.   叶片总RNA提取及cDNA文库构建

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

    • 1.3.   转录组测序及数据分析

  • 将测序获得的原始序列读长(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)功能分析。

    • 1.4.   基因表达量及差异基因分析

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

    • 1.5.   实时荧光定量PCR

  • 使用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次生物学重复。

    基因编号基因名称引物序列(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

    Table 1.  Real-time fluorescence quantitative gene primer information

2.   结果与分析

    2.1.   测序数据组装结果

  • 对照组(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)。

    样品名称读取碱基数/个总碱基数/个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

    Table 2.  Statistical table of sample sequencing data evaluation

    Figure 1.  Demo unigene length distribution

    • 2.2.   转录组单基因簇的数据库功能注释

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

    数据库注释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

    Table 3.  Unigene annotation statistics

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

    Figure 2.  Nr annotation distribution

    • 2.3.   响应高温胁迫的差异基因功能富集分析

  • 以差异倍数≥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)等过程。

    注释本体
    类别    		编号功能类型基因注释数量比例/%注释本体
    类别    		编号功能类型基因注释数量比例/%
    生物过程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

    Table 4.  Difference expression gene GO annotation statistics

    将高温组(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)。表明在高温刺激后，牡丹中的多种代谢途径、生物合成途径和内质网蛋白加工途径等均发生变化，积极参与应激响应。

    全基因组及代谢途径通路术语基因数目/个全基因组及代谢途径编号校正后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

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

    • 2.4.   实时荧光定量分析

  • 为了进一步验证转录组数据，随机选取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)、叶绿素酶(PsCHL1，c111201.graph_c0)、类胡萝卜素裂解双加氧酶(PsCCD4，c97637.graph_c0)、腺苷甲硫氨酸脱羧酶样原酶(PsSAMDC，C112526.graph_c0)以及假定蛋白(c120729.graph_c1)，使用qRT-PCR的方法进行验证分析(图3)，两者在变化趋势和差异倍数总体相符，表明测序结果相对可信。

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

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

    • 2.5.   高温胁迫下牡丹‘羽红’热激蛋白HSP基因挖掘及表达分析

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

    基因编号基因命名基因表达量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线粒体

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

3.   讨论

  • 全球变暖使得世界上许多地区农业受到不同程度的热害影响。近年来，植物耐热性研究逐渐受到关注。短期或持续高温胁迫引起植物形态和生理生化方面的改变，影响植物生长发育的同时也导致农林经济产量下降[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基因可在短期内迅速响应高温胁迫并参与牡丹的耐热调控。本研究可为进一步探究牡丹耐高温胁迫的分子机制提供依据。

Reference (34)

Catalog

    浙ICP备11046845号

    Copyright by 2019 Editorial Department of Journal of Zhejiang A&F University

    Address: No. 666, Wu Su street, Ling'an District, Hangzhou, ZhejiangPost Code: 311300Tel: 0571-63732749

    E-Mail: :zlxb@zafu.edu.cn

    Supported by: Beijing Renhe Information Technology Co. Ltd  support: info@rhhz.net

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return