• 中国中文核心期刊
  • 中国科学引文数据库(CSCD)核心库来源期刊
  • 中国科技论文统计源期刊(CJCR)
  • 第二届国家期刊奖提名奖

留言板

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

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

无瓣海桑根响应盐胁迫的转录组分析

梁锐涛 韩维栋 杨少瑕 陈蓓蓓

引用本文:
Citation:

无瓣海桑根响应盐胁迫的转录组分析

    通讯作者: 陈蓓蓓, beibeichenk@outlook.com
  • 中图分类号: S728.5

Transcriptome Analysis of Sonneratia apetala Root in Response to Salt Stress

    Corresponding author: CHEN Bei-bei, beibeichenk@outlook.com
  • CLC number: S728.5

  • 摘要: 摘要 : 目的 初步探究无瓣海桑的耐盐分子机制,筛选出无瓣海桑抗盐候选基因,为后续功能验证实验及林木抗盐性遗传育种奠定分子基础。 方法 以1年生无瓣海桑幼苗为材料,用500 mmol•L−1 NaCl分别处理0 d(对照组)和10 d(处理组),取不同条件下的根部组织进行转录组测序,并结合三代全长转录组数据进行后续生物信息学分析。 结果 (1)与对照组相比,NaCl处理10 d后,无瓣海桑幼苗根系中共有14401个差异表达基因,其中7153个上调,7248个下调。(2)GO分析发现,共有11068个差异基因在47个GO 条目得到注释。(3)在KEGG富集分析中,共有6189个差异基因富集到134条通路,其中,共有14条通路显著富集(P值<0.01,Q值<0.05)。(4)通过进一步对差异基因进行功能注释分析,共筛选出抗盐候选基因89个,其中,抗氧化基因24个,渗透调节物质基因22个,植物激素基因19个,蛋白激酶基因10个,转录因子基因14个。 结论 活性氧清除、渗透调节、植物激素、蛋白激酶及转录因子相关基因参与调控无瓣海桑盐逆境适应过程。
  • 图 1  6个样本表达量相关性热图

    Figure 1.  Correlation heat map of six samples based on gene expression levels

    图 2  组间差异基因火山图

    Figure 2.  Volcanic map of differential expression genes among different groups

    图 3  差异基因富集数目最多的20个GO条目

    Figure 3.  the 20 most abundant GO terms for differentially expressed genes

    图 4  差异基因最显著富集的20个KEGG通路

    Figure 4.  the first 20 significant enriched KEGG pathways for differentially expressed genes

    图 5  无瓣海桑根系响应盐胁迫关键基因

    Figure 5.  Key genes in response to salt stress in root of S. apetala

    图 6  89个无瓣海桑根系抗盐候选基因表达量聚类热图

    Figure 6.  Cluster heat map of the 89 salt-resistance candidate genes in roots of S. apetala

    表 1  高通量转录组测序数据统计表

    Table 1.  Statistical of high-throughput transcriptome sequencing data

    样本编号
    Sample ID
    总原始读数
    Total Raw Reads /M
    总有效读数
    Total Clean Reads /M
    比对率
    Mapping ratio/%
    Q20/%Q30/%有效读序百分率
    Clean Reads Ratio/%
    RCK_a43.8243.0787.4996.8991.9798.28
    RCK_b43.8243.2687.8996.9292.0398.71
    RCK_c43.8243.2387.9496.8991.9798.66
    RT_a43.8243.3588.0897.5393.2198.93
    RT_b43.8243.2387.6297.7893.8898.65
    RT_c43.8243.2587.6297.7793.8598.71
    注:RCK,空白对照组,RT,盐胁迫处理10天的处理组,a, b, c表示3个生物学重复;Q20、Q30,质量值大于或等于20、30的碱基数占总碱基数的比例。
    Notes: RCK represents blank control group, RT represents samples treated with salt for 10 days, and a, b, c represent the three biological repetitions; Q20, Q30, the proportion of bases with value ≥ 20, 30.
    下载: 导出CSV

    表 2  无瓣海桑根系部分抗盐候选功能基因

    Table 2.  List of partial salt-resistance candidate genes in root of S. apetala

    基因IDlog2FCKEGG功能注释参与通路
    活性氧清除 isoform_61079 6.46 超氧化物歧化酶,铜锌家族
    superoxide dismutase, Cu-Zn family
    过氧物酶体
    Peroxisome
    isoform_173361 1.47 超氧化物歧化酶,铁锰家族
    superoxide dismutase, Fe-Mn family
    isoform_287367 6.10 过氧化物酶
    peroxidase
    苯丙素生物合成
    Phenylpropanoid biosynthesis
    isoform_57012 −4.69
    isoform_25838 3.16 过氧化氢酶
    catalase
    MAPK信号通路-植物
    MAPK signaling pathway-plant
    isoform_58072 −1.83 谷胱甘肽过氧化物酶
    Glutathione peroxidase
    谷胱甘肽代谢
    Glutathione metabolism
    isoform_67198 5.39 黄酮醇合成酶
    flavonol synthase
    类黄酮生物合成
    Flavonoid biosynthesis
    渗透调节 isoform_17235 7.54 海藻糖6-磷酸磷酸酶
    trehalose 6-phosphate phosphatase
    淀粉和蔗糖代谢
    Starch and sucrose metabolism
    isoform_134633 −4.24 海藻糖6-磷酸合酶/磷酸酶
    trehalose 6-phosphate synthase/phosphatase
    isoform_9570 3.75 水苏糖合成酶
    stachyose synthetase
    谷胱甘肽代谢
    Galactose metabolism
    isoform_49176 7.98 棉子糖合成酶
    raffinose synthase
    isoform_16906 5.82 δ-1-吡咯啉-5-羧酸合成酶
    delta-1-pyrroline-5-carboxylate synthetase
    氨基酸的生物合成
    Biosynthesis of amino acids
    isoform_31419 −1.20
    isoform_272673 5.53 甜菜碱醛脱氢酶
    betaine-aldehyde dehydrogenase
    甘氨酸、丝氨酸和苏氨酸代谢
    Glycine, serine and threonine metabolism
    isoform_228012 2.21 脯氨酸脱氢酶
    proline dehydrogenase
    精氨酸和脯氨酸代谢
    Arginine and proline metabolism
    isoform_8893 −2.72
    植物激素 isoform_6363 8.35 生长素应答因子
    auxin response factor
    植物激素信号转导
    Plant hormone signal transduction
    isoform_119200 6.23 生长素应答GH3基因家族
    auxin responsive GH3 gene family
    isoform_126871 3.48 生长素响应蛋白
    auxin-responsive protein IAA
    isoform_278367 4.36 SAUR家族蛋白
    SAUR family protein
    isoform_50365 −7.94
    isoform_48396 6.73 脱落酸受体PYR/PYL家族
    abscisic acid receptor PYR/PYL family
    isoform_49285 −6.23
    isoform_221018 −2.92 蛋白磷酸酶2C
    protein phosphatase 2C
    isoform_267229 4.96 脱落酸应答元件结合因子
    ABA responsive element binding factor
    isoform_43536 −2.46
    蛋白激酶 isoform_216952 4.44 促分裂原活化蛋白激酶激酶激酶1
    mitogen-activated protein kinase kinase kinase1
    MAPK信号通路-植物
    MAPK signaling pathway-plant
    isoform_182097 4.28 促分裂原活化蛋白激酶4
    mitogen-activated protein kinase4
    isoform_140969 1.38 促分裂原活化蛋白激酶6
    mitogen-activated protein kinase 6
    isoform_161961 −7.52 钙依赖性蛋白激酶
    calcium-dependent protein kinase
    植物-病原体相互作用
    Plant-pathogen interaction
    下载: 导出CSV
  • [1] 胡 涛, 张鸽香, 郑福超, 等. 植物盐胁迫响应的研究进展[J]. 分子植物育种, 2018, 16(9):3006-3015.

    [2] 牛 恋. 植物耐盐生理机制及抗盐性[J]. 河北农机, 2018(2):56.

    [3] 欧阳少明. 无瓣海桑的造林技术和管理措施[J]. 农业科技与信息, 2016(15):139-139. doi: 10.3969/j.issn.1003-6997.2016.15.116

    [4] 廖 岩, 赵 肖, 陈桂珠. 盐胁迫对无瓣海桑幼苗根茎叶膜保护系统的影响[J]. 海洋环境科学, 2009, 28(02):154-158. doi: 10.3969/j.issn.1007-6336.2009.02.010

    [5] 高慧娟, 吕昕培, 王润娟, 等. 转录组测序在林草植物抗逆性研究中的应用[J]. 草业学报, 2019, 28(12):184. doi: 10.11686/cyxb2019352

    [6] 赵陆滟, 曹绍玉, 龙云树, 等. 全长转录组测序在植物中的应用研究进展[J]. 植物遗传资源学报, 2019, 20(6):1390-1398.

    [7] 李玉梅, 李书娴, 李向上, 等. 第三代测序技术在转录组学研究中的应用[J]. 生命科学仪器, 2018, 16(Z1):114-121 + 113.

    [8] 周卫星, 石海鹤. 高通量测序中序列拼接算法的研究进展[J]. 计算机科学, 2019, 46(5):36-43. doi: 10.11896/j.issn.1002-137X.2019.05.005

    [9] 马东娜, 张兴坦, 魏柳锋, 等. 基因组二代测序数据与三代测序数据的混合校正和组装[J]. 基因组学与应用生物学, 2018(04):1547-1555.

    [10] 梅 瑜, 李向荣, 蔡时可, 等. 药食同源植物甘葛藤的全长转录组分析[J]. 华北农学报, 2021, 36(5):10-17. doi: 10.7668/hbnxb.20192213

    [11] 郝赛琦. 红树植物白骨壤根系结构与转录组特征研究[D]. 厦门: 厦门大学, 2017: 11-57.

    [12]

    Su W, Ye C, Zhang Y, et al. Identification of putative key genes for coastal environments and cold adaptation in mangrove Kandelia obovata through transcriptome analysis[J]. Science of the total environment, 2019, 681: 191-201. doi: 10.1016/j.scitotenv.2019.05.127
    [13]

    Yang Y, Yang S, Li J, et al. De novo assembly of the transcriptomes of two yellow mangroves, Ceriops tagal and C. zippeliana, and one of their terrestrial relatives, Pellacalyx yunnanensis[J]. Marine genomics, 2015, 23: 33-36. doi: 10.1016/j.margen.2015.04.003
    [14]

    Guo W, Wu H, Zhang Z, et al. Comparative analysis of transcriptomes in Rhizophoraceae provides insights into the origin and adaptive evolution of mangrove plants in intertidal environments[J]. Frontiers in Plant Science, 2017, 8: 795. doi: 10.3389/fpls.2017.00795
    [15]

    Yang Y, Yang S, Li J, et al. Transcriptome analysis of the Holly mangrove Acanthus ilicifolius and its terrestrial relative, Acanthus leucostachyus, provides insights into adaptation to intertidal zones[J]. BMC genomics, 2015, 16(1): 1-12. doi: 10.1186/1471-2164-16-1
    [16] 刘婷婷, 莫玉剑, 欧成川, 陈蓓蓓. sss盐胁迫下无瓣海桑差异表达转录因子的转录组信息分析[J/OL]. 分子植物育种: 1-16[2022-05-22]. http://kns.cnki.net/kcms/detail/46.1068.S.20220407.1145.002.html.

    [17]

    Chen Y, Chen Y, Shi C, et al. SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data[J]. GigaScience, 2018, 7(1): gix120.
    [18]

    Langmead B, Salzberg S L. Fast gapped-read alignment with Bowtie 2[J]. Nature methods, 2012, 9(4): 357-359. doi: 10.1038/nmeth.1923
    [19]

    Vuong H, Truong T, Tran T, et al. A revisit of RSEM generative model and its EM algorithm for quantifying transcript abundances[J]. bioRxiv, 2018: 503672.
    [20]

    Love M I, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2[J]. Genome biology, 2014, 15(12): 1-21.
    [21]

    Kanehisa M, Araki M, Goto S, et al. KEGG for linking genomes to life and the environment[J]. Nucleic acids research, 2007, 36(suppl_1): D480-D484.
    [22] 吴 寒. 活性氧在植物体内的作用及其清除体制[J]. 广东蚕业, 2018, 52(3):18-18.

    [23] 郭明欣, 刘佳佳, 侯琳琳, 等. 植物体内活性氧的产生及清除机制研究进展[J]. 科技视界, 2021(8):104-106.

    [24] 黎 家, 李传友. 新中国成立70年来植物激素研究进展[J]. 中国科学:生命科学, 2019, 49(10):1227-1281.

    [25] 朱婷婷, 王彦霞, 裴丽丽, 等. 植物蛋白激酶与作物非生物胁迫抗性的研究[J]. 植物遗传资源学报, 2017, 18(4):763-770.

    [26]

    Wang Y, Xu W, Chen Z, et al. Gene structure, expression pattern and interaction of Nuclear Factor-Y family in castor bean (Ricinus communis)[J]. Planta, 2018, 247(3): 559-572. doi: 10.1007/s00425-017-2809-2
    [27]

    Zhu J K. Abiotic stress signaling and responses in plants[J]. Cell, 2016, 167(2): 313-324.
    [28] 付晴晴, 李商锐, 刘小瑞, 等. NaCl胁迫对不同耐盐性葡萄株系根系活性氧代谢的影响[J]. 中外葡萄与葡萄酒, 2017(3):12-15.

    [29] 鲁 琳, 杨尚谕, 刘维东, 等. 基于转录组测序花烟草响应盐胁迫活性氧清除相关基因的挖掘[J]. 生物技术通报, 2020, 36(12):42-53.

    [30]

    Kosar F, Akram N A, Sadiq M, et al. Trehalose: a key organic osmolyte effectively involved in plant abiotic stress tolerance[J]. Journal of Plant Growth Regulation, 2019, 38(2): 606-618. doi: 10.1007/s00344-018-9876-x
    [31] 李 辉, 李德芳, 邓 勇, 等. 红麻海藻糖生物合成关键酶基因HcTPPJ的克隆及响应逆境的表达分析[J]. 作物学报, 2020, 46(12):1914-1922.

    [32] 张莹莹, 吕 燕, 宣雯燕, 等. 龙须菜(Gracilariopsis lemaneiformis)海藻糖-6-磷酸合成酶(TPS)对逆境胁迫的响应[J]. 海洋与湖沼, 2021, 52(3):777-785. doi: 10.11693/hyhz20200700225

    [33]

    Zhang H, Dong J, Zhao X, et al. Research progress in membrane lipid metabolism and molecular mechanism in peanut cold tolerance[J]. Frontiers in Plant Science, 2019, 10: 838. doi: 10.3389/fpls.2019.00838
    [34]

    Funck D, Baumgarten L, Stift M, et al. Differential contribution of P5CS isoforms to stress tolerance in Arabidopsis[J]. Frontiers in plant science, 2020, 11: 1483.
    [35]

    Gong Z, Xiong L, Shi H, et al. Plant abiotic stress response and nutrient use efficiency[J]. Science China Life Sciences, 2020, 63(5): 635-674. doi: 10.1007/s11427-020-1683-x
    [36] 刘海洋, 陈玉珍. ABA受体结构及功能与ABA信号通路研究进展[J]. 中国农学通报, 2019, 35(21):75-81. doi: 10.11924/j.issn.1000-6850.casb19010079

    [37] 张 弛, 蔚静玲, 储谟立, 等. 拟南芥蛋白磷酸酶PP2C31的盐胁迫响应功能研究[J]. 中国科技论文, 2018, 13(18):2070-2075. doi: 10.3969/j.issn.2095-2783.2018.18.004

    [38]

    Zhang M, Su J, Zhang Y, et al. Conveying endogenous and exogenous signals: MAPK cascades in plant growth and defense[J]. Current Opinion in Plant Biology, 2018, 45(Pt A): 1-10.
    [39]

    Abulfaraj A A. Stepwise signal transduction cascades under salt stress in leaves of wild barley (Hordeum spontaneum)[J]. Biotechnology & Biotechnological Equipment, 2020, 34(1): 860-872.
    [40]

    Kumar K, Raina S K, Sultan S M. Arabidopsis MAPK signaling pathways and their cross talks in abiotic stress response[J]. Journal of Plant Biochemistry and Biotechnology, 2020, 29(4): 700-714. doi: 10.1007/s13562-020-00596-3
    [41] 汪芳珍, 杨成行, 何子华, 等. 盐处理下旱生植物沙芥蛋白激酶相关基因的差异表达分析[J]. 草业学报, 2021, 30(10):116-124. doi: 10.11686/cyxb2020368

    [42] 潘凌云, 马家冀, 李建民, 等. 植物盐胁迫应答转录因子的研究进展[J]. 生物工程学报, 2022, 38(1):50-65.

    [43]

    Du C, Ma B, Wu Z, et al. Reaumuria trigyna transcription factor RtWRKY23 enhances salt stress tolerance and delays flowering in plants[J]. Journal of Plant Physiology, 2019, 239: 38-51.
    [44]

    Lin J, Dang F, Chen Y, et al. CaWRKY27 negatively regulates salt and osmotic stress responses in pepper[J]. Plant Physiology and Biochemistry, 2019, 145: 43-51. doi: 10.1016/j.plaphy.2019.08.013
  • [1] 田林尹丹丹成铁龙夏新莉尹伟伦 . 盐胁迫下比拉底白刺差异表达基因的转录组分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2020.01.001
    [2] 郜新强王小艳焦伟王静郑丽月王丹宁王兴云侯立江尚增振刘彦珍 . 金丝楸幼苗响应盐碱胁迫的生理和转录组分析. 林业科学研究,
    [3] 周婧李卓蓉吴江婷 . 不同氮形态处理条件下杨树根尖差异表达基因的特征分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2022.02.006
    [4] 刘聪张洋夏德安陈雪冰魏志刚 . 毛果杨PLD基因家族全基因组水平鉴定及其盐胁迫下的表达分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2021.03.003
    [5] 李露双董文慧丁兴萃章志远孙春娃蔡函江 . 麻竹笋转录组测序及苦涩味物质合成基因差异表达分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2018.04.006
    [6] 袁长春莫健新袁柳娇余如凤刘倩钟军弟刘锴栋 . 重金属胁迫下白骨壤数字基因表达谱分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2017.02.004
    [7] 林雪锋颉洪涛虞木奎陈顺伟 . 盐胁迫下3种海滨植物形态和生理响应特征及耐盐性差异. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2018.03.013
    [8] 李云郑德璋陈焕雄廖宝文郑松发陈相如 . 红树植物无瓣海桑引种的初步研究. 林业科学研究,
    [9] 张振张含国周宇张磊于宏影张莉 . 红松4个组织的转录组数据分析与次生代谢产物的表达差异初探. 林业科学研究,
    [10] 昝启杰王勇军廖宝文郑德璋陈玉军 . 深圳福田无瓣海桑+海桑-秋茄人工林结构的研究. 林业科学研究,
    [11] 邓楠史胜青常二梅刘建锋兰倩江泽平 . 基于中麻黄萌发种子转录组的黄酮类化合物合成途径基因的挖掘. 林业科学研究,
    [12] 李云郑德璋廖宝文郑松发宋湘豫 . 盐度与温度对红树植物无瓣海桑种子发芽的影响. 林业科学研究,
    [13] 廖宝文李玫郑松发陈玉军郑馨仁 . 外来种无瓣海桑种内、种间竞争关系研究. 林业科学研究,
    [14] . 利用无瓣海桑控制入侵种互花米草的初步研究. 林业科学研究,
    [15] 朱可峰廖宝文章家恩 . 广州市南沙红树植物无瓣海桑、木榄人工林生物量的研究. 林业科学研究,
    [16] 李玫廖宝文郑松发陈玉军 . 无瓣海桑对乡土红树植物的化感作用. 林业科学研究,
    [17] 王臣虞木奎张翠吴统贵袁健军周成云 . 盐胁迫下3个楸树无性系光合特征研究. 林业科学研究,
    [18] 刘宇徐焕文刘桂丰姜静 . 赤霉素GA4+7处理下白桦无性系生长及差异基因表达分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2017.01.025
    [19] 王帅邵芬娟李论芦强邱德有 . 穗花杉的转录组测序及其转录组特性分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2017.05.008
    [20] 张华新宋丹刘正祥 . 盐胁迫下11个树种生理特性及其耐盐性研究. 林业科学研究,
  • 加载中
图(6) / 表(2)
计量
  • 文章访问数:  435
  • HTML全文浏览量:  193
  • PDF下载量:  5
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-03-29
  • 录用日期:  2022-05-31
  • 网络出版日期:  2022-08-04

无瓣海桑根响应盐胁迫的转录组分析

    通讯作者: 陈蓓蓓, beibeichenk@outlook.com
  • 广东海洋大学滨海农业学院,广东 湛江, 524088

摘要: 摘要 : 目的 初步探究无瓣海桑的耐盐分子机制,筛选出无瓣海桑抗盐候选基因,为后续功能验证实验及林木抗盐性遗传育种奠定分子基础。 方法 以1年生无瓣海桑幼苗为材料,用500 mmol•L−1 NaCl分别处理0 d(对照组)和10 d(处理组),取不同条件下的根部组织进行转录组测序,并结合三代全长转录组数据进行后续生物信息学分析。 结果 (1)与对照组相比,NaCl处理10 d后,无瓣海桑幼苗根系中共有14401个差异表达基因,其中7153个上调,7248个下调。(2)GO分析发现,共有11068个差异基因在47个GO 条目得到注释。(3)在KEGG富集分析中,共有6189个差异基因富集到134条通路,其中,共有14条通路显著富集(P值<0.01,Q值<0.05)。(4)通过进一步对差异基因进行功能注释分析,共筛选出抗盐候选基因89个,其中,抗氧化基因24个,渗透调节物质基因22个,植物激素基因19个,蛋白激酶基因10个,转录因子基因14个。 结论 活性氧清除、渗透调节、植物激素、蛋白激酶及转录因子相关基因参与调控无瓣海桑盐逆境适应过程。

English Abstract

  • 土壤盐渍化现象几乎遍布所有的气候区域以及地理环境,是影响农业和林业生产的主要问题。土壤盐渍化使植物受到渗透胁迫、离子毒害、膜透性改变及生理代谢紊乱等危害[1],严重制约了植物的生长,是在植物生产与栽培上一个非常严峻的问题。因此,通过对植物抗盐的分子调控机制进行探析,挖掘耐盐相关基因,并应用于新种质创制和新品种选育,对于缓解土壤盐渍化问题具有重要的现实意义。

    长期生长在滨海区域潮间带的红树林,由于其特殊的、有别于陆地和淡水植物的生长环境,演化出了一套特殊的盐逆境适应机制。植物的耐盐机制涉及渗透物质积累、蛋白积累、抗氧化酶的诱导、植物激素的诱导和光合作用路径的变化5个方面[2]。对于红树植物,不同的红树物种具体的耐盐机制会有所不同。无瓣海桑(Sonneratia apetala Buch.-Ham.)是海桑科海桑属红树乔木,具有生长快、耐淹浸、耐贫瘠、耐寒强等特点,是滨海滩涂的先锋造林树种,在维持海岸带生态平衡方面起着重要作用[3]。先前有关其耐盐特性的研究主要集中在形态、生理生态层面,如廖岩等[4]研究了不同盐度处理3个月下无瓣海桑根、茎、叶中的蛋白质、可溶性糖、丙二醛含量、抗氧化酶活性的变化情况,并探讨了它们的相互关系,为耐盐树种的选育提供了参考。然而,目前鲜见关于无瓣海桑耐盐分子机制方面的研究报道。

    转录组测序是转录组学研究的核心技术,被广泛应用于新基因的挖掘等[5]。目前,最新的测序技术已经发展到了第三代。三代测序技术具有单分子、超长读长、无需PCR扩增等优点,已在基因组学、转录组学、表观遗传学等领域上得到了应用[6]。在转录组学上,三代全长转录组测序(Iso-Seq)的应用,可弥补二代测序技术读段较短、无法覆盖整个转录本的缺点,从而做到从全长转录本水平上研究转录组[7]。但三代测序技术存在碱基错误率高、成本高等缺点,目前应用并不广泛[8]。二代测序有通量高、成本低的特点,且发展至今已取得了较好的工作基础,因此,二代和三代测序结合使用成为很多研究者的策略[9]。如梅瑜等[10]利用基于SMRT技术原理的PacBio平台三代测序技术,对甘葛藤(Pueraria thomsonii Benth)进行全长转录组测序及分析,从转录数据中鉴定了黄酮类生物合成基因、转录因子、R基因和SSR标记,为甘葛藤的选育和利用提供了基础。然而,对于红树植物,当前研究仅限于对白骨壤(Avicennia marina (Forsk.) Vierh.)[11]、秋茄(Kandelia candel (Linn.) Druce / Kandelia obovata)[12]、角果木(Ceriops tagal (perr.) C. B. Rob)[13]、竹节树(Carallia brachiate (Lour.) Merr.)[14]和老鼠簕(Acanthus ilicifolius L.)[15]等个别物种开展了转录组学相关研究。特别是,至今对无瓣海桑仅开展了以三代全长转录组数据集为基础的盐胁迫应答转录因子挖掘工作[16]。因此,本研究利用二代与三代结合的高通量测序技术,对盐胁迫处理前后的无瓣海桑根系进行转录组测序分析,从中鉴定并初步筛选出无瓣海桑根系响应盐胁迫的候选基因,可为深入研究无瓣海桑的耐盐基因组学、分子生物学和进一步揭示红树植物的耐盐机制提供理论基础。

    • 本实验原材料采自于广东省湛江市特呈岛红树林自然保护区(21°09′ ~ 21°10′ N, 110°25′ ~ 110°27′ E)。2019年10月采集无瓣海桑的种子,播种于人工配制的营养土中。一年后,选取长势一致的无瓣海桑幼苗,分别设置对照组(RCK)和处理组(RT)。处理组用500 mmol•L−1的NaCl溶液每3 d进行一次盐胁迫处理,对照组使用等量清水。对照组(RCK_a、RCK_b、RCK_c)和处理组(RT_a、RT_b、RT_c)均进行了3次生物学重复。10 d后分别取对照组和处理组根部组织液氮速冻,于−80 ℃冰箱保存。

    • 提取6个样本RNA,质检合格后构建三代全长转录组测序文库和二代RNA-seq测序文库,并进行测序及分析。二代测序文库构建方法如下:(1)用mRNA富集法或rRNA去除法对total RNA进行处理;(2)用打断buffer把获得的RNA片段化,随机的N6引物进行反转录,合成cDNA二链形成双链DNA,再将其末端补平并5′端磷酸化;(3)连接一个3'端有凸出“T”的鼓泡状的接头;(4)通过特异的引物进行PCR扩增得到PCR产物,将其热变性成单链,再用一段桥式引物将单链DNA环化得到单链环状DNA文库;(5)利用BGISEQ-500平台测序。

      PacBio测序文库构建方法如下:取6个样本RNA,等量混合,使用SMARTer™ PCR cDNA Synthesis试剂盒将混合的RNA样品反转录合成第一链cDNA,经PCR扩增合成第二链cDNA,然后对双链DNA经过二次PCR扩增后的片段进行损伤修复、末端修复,并连接SMRT适配体后构建全长转录组测序文库,利用PacBio平台进行测序。

      文库构建及测序工作由华大基因公司完成。

    • 为保证数据质量,需对测序得到的原始数据进行过滤、参考基因比对等前期处理。本研究使用华大自主研发的过滤软件SOAPnuke[17]对6个样本进行过滤,去除包含接头污染、未知碱基N含量大于10%、低质量的reads。使用Bowtie2[18]将clean reads比对到三代全长转录组数据库;利用RSEM计算各个样品的基因表达水平[19]。样本间的相关性分析通过皮尔逊相关系数(Pearson correlation coefficient)进行。

    • 本研究基于负二项分布原理的Deseq2模型分析方法[20],对两组样本进行差异检测分析,并以差异倍数值(fold change, FC)≥ 2或 ≤ 0.5,P-value < 0.01,Q-value < 0.05作为筛选差异基因的标准,并将该筛选得到的差异基因用于后续的GO和KEGG Pathway富集分析,前者通过agriGO网站(http://systemsbiology.cau.edu.cn/agriGOv2/)进行,后者借助基因数据库网站(https://www.kegg.jp/[21]完成。

    • 通过SOAPnuke软件对原始数据进行过滤后,如表1所示,有效读序百分率均在98%以上,Q20(碱基量 ≥ 20%)均大于96%且Q30(碱基量 ≥ 30%)均 > 91%,比对率均在87%以上。6个样品转录本reads覆盖度的峰值在90 ~ 100;空白对照RCK组和盐处理RT组的三次生物学重复间的相关系数极高( > 0.98,图1)。以上数据均说明6个样品的测序质量较高,测序数据可靠,可用于后续的差异表达分析。

      表 1  高通量转录组测序数据统计表

      Table 1.  Statistical of high-throughput transcriptome sequencing data

      样本编号
      Sample ID
      总原始读数
      Total Raw Reads /M
      总有效读数
      Total Clean Reads /M
      比对率
      Mapping ratio/%
      Q20/%Q30/%有效读序百分率
      Clean Reads Ratio/%
      RCK_a43.8243.0787.4996.8991.9798.28
      RCK_b43.8243.2687.8996.9292.0398.71
      RCK_c43.8243.2387.9496.8991.9798.66
      RT_a43.8243.3588.0897.5393.2198.93
      RT_b43.8243.2387.6297.7893.8898.65
      RT_c43.8243.2587.6297.7793.8598.71
      注:RCK,空白对照组,RT,盐胁迫处理10天的处理组,a, b, c表示3个生物学重复;Q20、Q30,质量值大于或等于20、30的碱基数占总碱基数的比例。
      Notes: RCK represents blank control group, RT represents samples treated with salt for 10 days, and a, b, c represent the three biological repetitions; Q20, Q30, the proportion of bases with value ≥ 20, 30.

      图  1  6个样本表达量相关性热图

      Figure 1.  Correlation heat map of six samples based on gene expression levels

    • 无瓣海桑根系在500 mmol•L−1的NaCl胁迫10 d后,共检测到20289个基因表达量发生了改变(图2)。两组样本间共有14401个基因差异表达,其中7153个差异基因上调表达,log2FC最大值为11.32(isoform_13051);7248个差异基因下调表达,log2FC最低值为−12.69(isoform_11580)。结果表明,下调表达的差异基因数量多于上调的数量。

      图  2  组间差异基因火山图

      Figure 2.  Volcanic map of differential expression genes among different groups

    • 将差异基因比对到Gene Ontology (GO)数据库,共有11068个差异基因注释到47条GO条目,注释率为76.9%,涉及3大主要功能分类:生物过程、细胞组分和分子功能,占比分别为25.02%,43.76%,31.21%,由此可知大多数差异基因的功能与细胞组分具有较大的相关性。根据GO功能注释分类结果(图3),这些差异基因在参与生物过程中,主要集中在细胞过程和代谢过程;在细胞组分中,主要集中在细胞、膜、膜部分;在分子功能中,主要集中在连接、催化活性。表明在高盐胁迫下,上述的亚类中涉及到的基因可能在无瓣海桑根系中响应强烈。

      图  3  差异基因富集数目最多的20个GO条目

      Figure 3.  the 20 most abundant GO terms for differentially expressed genes

    • 在KEGG分析结果中,共有6 189个差异基因富集到134条通路,以Q-value<0.05作为显著富集标准,共有14条显著富集的通路(图4),其中,富集差异基因数量最多的是植物激素信号转导(ko04075),涉及512个基因,最显著富集的植物-病原体相互作用(ko04626)涉及367个基因。富集分析结果表明,无瓣海桑根系基因主要通过14条显著通路参与盐胁迫的响应,并且其中涉及的基因可能对无瓣海桑根系耐盐性的调控发挥作用,关注这些显著通路,有助于对无瓣海桑根系耐盐候选基因的筛选。

      图  4  差异基因最显著富集的20个KEGG通路

      Figure 4.  the first 20 significant enriched KEGG pathways for differentially expressed genes

    • 通过进一步对差异基因进行功能注释分析,本研究共筛选出75个无瓣海桑根系抗盐候选功能基因,涉及12条KEGG通路(表2图5)。

      表 2  无瓣海桑根系部分抗盐候选功能基因

      Table 2.  List of partial salt-resistance candidate genes in root of S. apetala

      基因IDlog2FCKEGG功能注释参与通路
      活性氧清除 isoform_61079 6.46 超氧化物歧化酶,铜锌家族
      superoxide dismutase, Cu-Zn family
      过氧物酶体
      Peroxisome
      isoform_173361 1.47 超氧化物歧化酶,铁锰家族
      superoxide dismutase, Fe-Mn family
      isoform_287367 6.10 过氧化物酶
      peroxidase
      苯丙素生物合成
      Phenylpropanoid biosynthesis
      isoform_57012 −4.69
      isoform_25838 3.16 过氧化氢酶
      catalase
      MAPK信号通路-植物
      MAPK signaling pathway-plant
      isoform_58072 −1.83 谷胱甘肽过氧化物酶
      Glutathione peroxidase
      谷胱甘肽代谢
      Glutathione metabolism
      isoform_67198 5.39 黄酮醇合成酶
      flavonol synthase
      类黄酮生物合成
      Flavonoid biosynthesis
      渗透调节 isoform_17235 7.54 海藻糖6-磷酸磷酸酶
      trehalose 6-phosphate phosphatase
      淀粉和蔗糖代谢
      Starch and sucrose metabolism
      isoform_134633 −4.24 海藻糖6-磷酸合酶/磷酸酶
      trehalose 6-phosphate synthase/phosphatase
      isoform_9570 3.75 水苏糖合成酶
      stachyose synthetase
      谷胱甘肽代谢
      Galactose metabolism
      isoform_49176 7.98 棉子糖合成酶
      raffinose synthase
      isoform_16906 5.82 δ-1-吡咯啉-5-羧酸合成酶
      delta-1-pyrroline-5-carboxylate synthetase
      氨基酸的生物合成
      Biosynthesis of amino acids
      isoform_31419 −1.20
      isoform_272673 5.53 甜菜碱醛脱氢酶
      betaine-aldehyde dehydrogenase
      甘氨酸、丝氨酸和苏氨酸代谢
      Glycine, serine and threonine metabolism
      isoform_228012 2.21 脯氨酸脱氢酶
      proline dehydrogenase
      精氨酸和脯氨酸代谢
      Arginine and proline metabolism
      isoform_8893 −2.72
      植物激素 isoform_6363 8.35 生长素应答因子
      auxin response factor
      植物激素信号转导
      Plant hormone signal transduction
      isoform_119200 6.23 生长素应答GH3基因家族
      auxin responsive GH3 gene family
      isoform_126871 3.48 生长素响应蛋白
      auxin-responsive protein IAA
      isoform_278367 4.36 SAUR家族蛋白
      SAUR family protein
      isoform_50365 −7.94
      isoform_48396 6.73 脱落酸受体PYR/PYL家族
      abscisic acid receptor PYR/PYL family
      isoform_49285 −6.23
      isoform_221018 −2.92 蛋白磷酸酶2C
      protein phosphatase 2C
      isoform_267229 4.96 脱落酸应答元件结合因子
      ABA responsive element binding factor
      isoform_43536 −2.46
      蛋白激酶 isoform_216952 4.44 促分裂原活化蛋白激酶激酶激酶1
      mitogen-activated protein kinase kinase kinase1
      MAPK信号通路-植物
      MAPK signaling pathway-plant
      isoform_182097 4.28 促分裂原活化蛋白激酶4
      mitogen-activated protein kinase4
      isoform_140969 1.38 促分裂原活化蛋白激酶6
      mitogen-activated protein kinase 6
      isoform_161961 −7.52 钙依赖性蛋白激酶
      calcium-dependent protein kinase
      植物-病原体相互作用
      Plant-pathogen interaction

      图  5  无瓣海桑根系响应盐胁迫关键基因

      Figure 5.  Key genes in response to salt stress in root of S. apetala

    • 正常情况下,植物体内存在着活性氧(reactive oxygen species, ROS)产生和清除的动态平衡,但外界的不良环境会打破这个平衡,导致ROS的积累,产生毒害[22]。植物体内存在两种ROS清除机制——酶促清除和非酶促清除,前者主要涉及超氧化物歧化酶(superoxide dismutase, SOD)、过氧化氢酶(catalase, CAT)等,后者主要涉及抗坏血酸、类黄酮物质等[23]。本研究发现6类共24个活性氧清除基因(图6),其中,上调表达基因数量最多的是过氧化物酶(7),倍数最高为6.46倍(超氧化物歧化酶,铜锌家族超氧化物歧化酶)。这些上调表达的基因有利于无瓣海桑在盐胁迫前期,合成多种抗氧化酶类物质,提高体内抗氧化和ROS清除的能力,增强耐受性。

      图  6  89个无瓣海桑根系抗盐候选基因表达量聚类热图

      Figure 6.  Cluster heat map of the 89 salt-resistance candidate genes in roots of S. apetala

    • 植物在对盐胁迫的适应过程中,细胞会积累一定量的可溶性有机物质,作为渗透调节物质参与渗透调节,以适应外界的低水势[2]。在本研究中,发现10种共22个编码生成有机渗透调节物质的基因(图6),涉及水苏糖、棉子糖、蔗糖、海藻糖、脯氨酸、甜菜碱、热激蛋白等渗透调节物质。在这些基因中,部分基因存在上调表达,表明在盐胁迫下,无瓣海桑可通过合成一些渗透调节物质,以调节盐渍环境下植物细胞内的渗透势。

    • 植物体内的多种激素在其生长发育、代谢等生理过程中发挥着重要的调控作用[24]。本研究,发现7种共19个涉及2种植物激素(IAA与ABA)的基因(图6),其中,IAA途径相关的有生长素应答因子(auxin response factor, ARF)、生长素应答GH3家族基因(auxin responsive GH3 gene family, GH3)、SAUR家族蛋白(SAUR family protein, SAUR)和生长素响应蛋白(auxin-responsive protein IAA, IAA)共11个基因;涉及ABA途径的有脱落酸应答元件结合因子(ABA responsive element binding factor, ABF)、脱落酸受体PYR/PYL家族(abscisic acid receptor PYR/PYL family, PYL)和蛋白磷酸酶2C(protein phosphatase 2C, PP2C)共8个基因。IAA相关基因除GH3 基因均上调表达外, 其他同时存在上下调表达,而ABA相关基因大部分下调表达。上述结果表明,植物激素能够在无瓣海桑的信号转导途径中发挥作用,不同激素共同调控盐胁迫响应过程。

    • 蛋白激酶可调控植物的生长,在逆境中起到调控感受胁迫信号、启动各种非生物逆境响应的作用[25]。 本研究发现6种共10个蛋白激酶相关基因(图6),其中2种促分裂原活化蛋白激酶基因(mitogen-activated protein kinases, MAPKMPK4MPK6和2种丝氨酸/苏氨酸蛋白激酶基因PBS1CTR1均上调表达,1种促分裂原活化蛋白激酶基因(mitogen activated protein kinase, MEKKMEKK1和1种钙依赖性蛋白激酶基因(calcium-dependent protein kinase, CDPK)中均有2个上调表达,1个下调表达。可见这些蛋白激酶在无瓣海桑根系响应高盐胁迫时,发挥着尤为重要的调控作用。

    • 转录因子是一类功能蛋白质,具有调控基因表达的作用。当植物受到逆境胁迫时,转录因子通过与相应的顺式作用元件相互作用,调控并减轻逆境胁迫带来的伤害,以适应各种逆境[26]。本研究发现8种共14个差异表达的转录因子(图6),涉及5个家族,分别为WRKY家族、MYB家族、GRAS家族、bHLH家族和EREBP家族,其中WRKY家族发现的成员数量最多,暗示WRKY家族在无瓣海桑响应盐胁迫中起着重要的调控作用。

    • 植物对盐胁迫的响应是一个多基因参与的复杂过程,从微观上的基因响应到宏观上生理的变化,经历了一系列的信号传递与调节过程,最终以适应盐逆境[27]。通过转录组测序技术对植物盐胁迫下的转录调控机制进行研究分析,能够快速准确地初步筛选出其中的耐盐相关基因,并为进一步利用分子生物学手段深入探究其耐盐机制奠定数据基础。然而,对于红树植物,当前研究对无瓣海桑仅限于利用转录组学技术开展了盐胁迫相关转录因子的挖掘[16],而对转录因子以外的其他可能参与调控盐逆境适应过程的基因并未关注。本研究联合二代和三代转录组测序技术,获得了无瓣海桑对照组及500 mmol•L−1盐胁迫处理10 d后的根系转录组数据,在此基础上进行生物信息学分析,最终筛选出89个无瓣海桑根抗盐候选基因,其中一些基因对盐胁迫的响应起着十分重要的作用。

      在植物的ROS酶促清除机制中,SOD、过氧化物酶(peroxidase, POD)、CAT是常见的酶。研究表明,NaCl胁迫后耐盐性强的A34株系葡萄根系中的SOD、CAT、POD活性表现出快速升高并在较高水平后小幅度降低的趋势,这有利于维持活性氧代谢的平衡,缓解细胞膜脂过氧化[28]。鲁琳等[29]通过高通量测序技术对盐胁迫下花烟草(Nicotiana alata Link et Otto)的活性氧清除基因进行挖掘,发现编码SOD、POD、CAT的基因均发生了显著的上调表达,同时存在下调表达的基因。本研究也发现了类似的结果:在SOD中,SOD1 2个上调,1个下调,SOD2 1个上调表达;CAT 5个基因全部上调表达;在POD中,存在较多下调表达的基因。上述结果推测可能是由于植物体内的抗氧化机制存在着正负两种调控方式,且不同抗氧化酶在不同植物中的调控方式可能有所不同;对于无瓣海桑来说,CAT可能是以正向调控为主导,而POD可能主要是以负向调控为主。

      在逆境下,植物可通过调节海藻糖的含量,稳定蛋白质的结构与功能,提高在逆境中的生存能力[30]。李辉等[31]克隆了海藻糖合成关键酶基因海藻糖-6-磷酸磷酸酶基因(trehalose-6-phosphate phosphatase, TPP),并研究其在非生物胁迫下的作用,结果表明该基因表达显著上调。研究表明,高盐胁迫主要促进了龙须菜(Asparagus schoberioides Kunth)海藻糖6-磷酸合酶/磷酸酶1基因(trehalose 6-phosphate synthase/phosphatase 1,TPS1)、海藻糖6-磷酸合酶/磷酸酶2基因(trehalose 6-phosphate synthase/phosphatase 2, TPS2)和海藻糖6-磷酸合酶/磷酸酶4基因(trehalose 6-phosphate synthase/phosphatase 4, TPS4)的表达,而渗透胁迫抑制了TPS1TPS2和海藻糖6-磷酸合酶/磷酸酶3基因(trehalose 6-phosphate synthase/phosphatase 3, TPS3)的表达[32]。本研究发现2种共7个合成海藻糖的关键基因:海藻糖6-磷酸磷酸酶基因(otsB)和海藻糖6-磷酸合酶基因(TPS)。在这些基因中,3个表现出下调表达,可能是由于盐胁迫产生的渗透胁迫抑制了这些基因的表达。

      脯氨酸是植物中常见的有机渗透调节物质,而脯氨酸的积累是植物应对盐胁迫的普遍反应。在大多数植物物种中,δ-1-吡咯啉-5-羧酸合成酶(P5CS)催化谷氨酸生物合成脯氨酸的第一步[33]。在Funck等[34]的研究中,发现在拟南芥(Arabidopsis thaliana (L.) Heynh)NaCl处理的p5cs2突变体中,几乎没有Na+的积累,植株既没有褪绿,也没有光合作用减弱,可见P5CS2介导的脯氨酸合成在调节叶片Na+积累和抗盐胁迫方面具有重要调控功能。本研究中发现2个P5CS基因,1个上调,1个下调,表明P5CS基因在渗透调节中发挥了作用。

      脱落酸在盐、干旱、低温等非生物胁迫下会在植物体内被诱导产生,同时启动相关的信号途径,以响应非生物胁迫[35]。ABA信号途径由ABA受体、PP2C、SnRK2和ABF转录因子四部分组成[36]。研究表明,盐处理抑制了拟南芥中PP2C31基因的表达,且在不同时期盐处理下PP2C31基因出现了不同程度的下调表达,说明PP2C31蛋白在拟南芥响应盐胁迫中起负调控作用[37]。本研究发现4个PP2C基因均下调表达,推测PP2C蛋白负调控无瓣海桑的盐胁迫响应。

      丝裂原活化蛋白激酶(MAPK)途径在盐、高温、病原菌等逆境的应答中起着重要作用,它通过MAPKKK→MAPKK→MAPK逐级传递放大胁迫信号,由下游的MAPK激活调控相应基因的表达,引起植物对逆境的应答[38]。在野生大麦(Hordeum spontaneum L.)中,MEKK1-MKK2-MPK4/6[39]联级系统参与了盐胁迫的响应。拟南芥的MAPK联级中多个MAPK基因在多种非生物胁迫的响应中起着重要作用[40]。汪芳珍等[41]通过对沙芥根进行6 h 50 mmol•L−1的盐处理实验,结果表明,共有11个MAPK/MAPKK/MAPKKK蛋白激酶基因的表达显著上调。本研究中发现MEKK1MPK4MPK6基因大都上调表达,推测这些基因在无瓣海桑的MAPK联级途径响应盐胁迫中起着重要的作用,但无瓣海桑具体响应盐胁迫的MAPK联级系统中的基因及机制,还有待进一步研究。

      植物体转录因子在盐胁迫调控下游机制中发挥着重要作用。参与植物盐胁迫响应的转录因子家族有WRKY、MYB、bHLH、NAC等[42]。在盐胁迫的响应中,WRKY家族转录因子存在着正负2种调控机制。在拟南芥中,RtWRKY23可提高自身过氧化物酶水平,以增强盐胁迫耐受性[43]。转基因烟草中过表达的CaWRKY27抑制了抗氧化酶合成途径相关基因的表达,在响应盐胁迫中起负调控作用[44]。本研究发现WRKY22WRKY25均呈上调表达,说明这些转录因子可能正向调控无瓣海桑的耐盐性。此外,WRKY2WRKY33呈现上下调表达,这表明了WRKY家族转录因子参与无瓣海桑耐盐机制方式的复杂性。

    • 本研究通过二代RNA-seq和三代全长转录组测序技术的联合使用对无瓣海桑响应盐逆境的转录调控机制进行了解析,获得了盐逆境下无瓣海桑根系的差异表达基因及功能注释信息,并鉴定出一批与活性氧清除、渗透调节、植物激素、蛋白激酶、转录因子等相关的潜在耐盐候选基因,研究结果对揭示红树植物的抗盐机制、耐盐型红树新型品种选育及沿海生态环境保护具有重要意义。

参考文献 (44)

目录

    /

    返回文章
    返回