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

留言板

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

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

金丝楸幼苗响应盐碱胁迫的生理和转录组分析

郜新强 王小艳 焦伟 王静 郑丽月 王丹宁 王兴云 侯立江 尚增振 刘彦珍

引用本文:
Citation:

金丝楸幼苗响应盐碱胁迫的生理和转录组分析

    通讯作者: 刘彦珍, 1135953615@qq.com
  • 中图分类号: S722

Physiological and transcriptomic analysis of Catalpa bungei seedlings in response to saline-alkali stresses

    Corresponding author: Liu Yan-zhen, 1135953615@qq.com
  • CLC number: S722

  • 摘要: 目的 研究不同盐碱胁迫对金丝楸幼苗生长、光合和生理指标的影响,并结合转录组测序分析,探究楸树耐盐碱的生理机制和分子机制。 方法 采用盆栽法对金丝楸幼苗进行不同盐碱胁迫处理,分析其生物量、光合及生理指标对不同盐碱响应的差异,采用Illumina高通量测序技术进行转录组测序,通过生物信息学分析盐碱胁迫对转录水平的影响。 结果 不同盐碱胁迫下,金丝楸幼苗叶片受伤害程度为Na2CO3>混合盐碱>NaCl;新增株高和地径、地上部和根的干质量和鲜质量、生物量、根冠比均受到明显抑制,并随盐碱浓度增加而抑制加强,但生长胁迫指数均随浓度的增加而降低;丙二醛(MDA)含量和相对电导率都随胁迫浓度增加而不同程度上升,超氧化物歧化酶(SOD)活性、可溶性糖含量、脯氨酸(Pro)含量、叶绿素总量和光合速率都呈现先上升后下降的趋势。转录组测序共产生约60.4 Gb原始数据,组装得到55 793个Unigenes,其中29 534(52.93%)个Unigenes获得了注释;通过差异表达基因(DEGs)分析,3个比较组(CK vs NaCl,CK vs Na2CO3 和 CK vs 混合盐碱)分别筛选出1 779、2 835和4 059 个DEGs;DEGs GO富集分析表明,膜的整体成分、膜的内在成分、催化活性、类异戊二烯代谢和合成过程、氧化还原酶活性等条目被显著富集;DEGs KEGG分析表明,苯丙素生物合成、淀粉和蔗糖代谢、植物激素信号转导、萜类主干生物合成和精氨酸代谢等通路被显著富集;此外,在DEGs中鉴定的bHLH、ERF、MYB-related、NAC、C2H2、WRKY、MYB和bZIP转录因子家族成员最多。 结论 金丝楸主要通过积累可溶性糖和Pro,提高SOD酶活和光合作用来抵御盐碱胁迫,但都呈现“低促高抑”的现象,说明其具有一定阈值。金丝楸通过调节膜成分、催化活性、类异戊二烯代谢和生物合成过程、苯丙素生物合成、淀粉和蔗糖代谢、植物激素信号转导等生物过程和代谢途径,并结合有关转录因子共同响应盐碱胁迫。本研究为深入研究楸树耐盐碱生理机制和分子机制提供科学的理论依据。
  • 图 1  金丝楸不同伤害等级的叶片盐碱胁迫症状

    Figure 1.  The symptom of different damage level on leaves of C. bungei under different saline-alkali stress

    图 2  特异和共有差异表达基因的维恩图分析

    Figure 2.  Venn diagrams analysis of specific and shared DEGs

    图 3  DEGs GO富集分析(A Vs C)

    Figure 3.  GO enrichment analysis of the DEGs (A Vs C)

    图 4  DEGs GO富集分析(A Vs D)

    Figure 4.  GO enrichment analysis of the DEGs (A Vs D)

    图 5  DEGs GO富集分析(A Vs E)

    Figure 5.  GO enrichment analysis of the DEGs (A Vs E)

    图 6  DEGs KEGG富集分析(A Vs C)

    Figure 6.  KEGG enrichment analysis of the DEGs (A Vs C)

    图 7  DEGs KEGG富集分析(A Vs D)

    Figure 7.  KEGG enrichment analysis of the DEGs (A Vs D)

    图 8  DEGs KEGG富集分析(A Vs E)

    Figure 8.  KEGG enrichment analysis of the DEGs (A Vs E)

    图 9  结论示意图

    Figure 9.  Schematic diagram of conclusion

    表 1  不同盐碱胁迫对金丝楸幼苗生长的影响

    Table 1.  Effects of different saline-alkali stress on the growth of C. bungei seedlings

    处理
    Treatment
    新增株高
    Net growth of plant height/
    cm
    新增地径
    Net growth of stem diameter/
    mm
    地上部overground part
    鲜质量
    Fresh weight/g
    干质量
    Dry weight/g
    CK018.68 + 0.621.23 + 0.0924.60 + 1.145.70 + 0.19
    NaCl胁迫
    NaCl stress
    5012.58 + 0.62**1.50 + 0.13*27.06 + 1.086.75 + 0.25*
    10010.00 + 0.53**0.98 + 0.09*15.94 + 0.07**3.31 + 0.05**
    1508.14 + 0.36**0.66 + 0.10**10.13 + 0.14**2.14 + 0.09**
    2005.92 + 0.79**0.40 + 0.07**6.92 + 0.18**1.54 + 0.08**
    Na2CO3胁迫
    Na2CO3 stress
    5017.64 + 0.98*1.26 + 0.1018.87 + 0.62**4.68 + 0.18**
    10012.68 + 0.58**1.03 + 0.0911.20 + 0.84**3.14 + 0.22**
    15010.48 + 1.36**0.57 + 0.09**6.87 + 0.24**1.73 + 0.09**
    2006.58 + 0.74**0.37 + 0.05**4.80 + 0.10**1.12 + 0.05**
    混合盐碱胁迫
    Mixed saline-alkali
    stress
    5016.74 + 0.77**1.18 + 0.0835.18 + 0.42**9.13 + 0.20**
    10013.58 + 1.06**0.84 + 0.1315.56 + 0.42**3.35 + 0.29**
    15010.10 + 0.62**0.59 + 0.08**13.09 + 0.47**2.67 + 0.30**
    2009.50 + 1.22**0.37 + 0.08** 6.76 + 0.33** 1.40 + 0.05**
    处理
    Treatment
    根 root生物量
    Biomass/g
    根冠比
    Root-shoot ratio
    生长胁迫指数
    Growth stress index
    鲜质量
    Fresh weight/g
    干质量
    Dry weight/g
    CK08.99 + 0.151.28 + 0.056.97 + 0.240.22 + 0.01
    NaCl胁迫
    NaCl stress
    508.00 + 0.16*0.90 + 0.03*7.65 + 0.270.13 + 0.01**1.10 + 0.04
    1004.08 + 0.09**0.48 + 0.03**3.97 + 0.07**0.15 + 0.01**0.54 + 0.01
    1503.29 + 0.08**0.33 + 0.03**2.47 + 0.12**0.15 + 0.01**0.35 + 0.02
    2002.57 + 0.07**0.26 + 0.04**1.79 + 0.11**0.17 + 0.02*0.26 + 0.02
    Na2CO3胁迫
    Na2CO3 stress
    504.90 + 0.27**0.74 + 0.09*5.42 + 0.20**0.16 + 0.02*0.78 + 0.03
    1002.55 + 0.43**0.31 + 0.07**3.44 + 0.29**0.10 + 0.01**0.49 + 0.04
    1502.51 + 0.45**0.27 + 0.04**2.00 + 0.12**0.16 + 0.01*0.29 + 0.02
    2001.31 + 0.19**0.13 + 0.02**1.25 + 0.07**0.12 + 0.01**0.18 + 0.01
    混合盐碱胁迫
    Mixed saline-alkali
    stress
    5012.44 + 0.90*2.02 + 0.18*11.15 + 0.37**0.22 + 0.021.60 + 0.05
    1006.09 + 0.92*0.65 + 0.06**4.00 + 0.34**0.19 + 0.02**0.57 + 0.05
    1502.62 + 0.56**0.28 + 0.04**2.95 + 0.33**0.10 + 0.01**0.42 + 0.05
    2001.86 + 0.27**0.22 + 0.04**1.62 + 0.08**0.15 + 0.02*0.23 + 0.01
    注:“*”表示处理间差异显著(P<0.05),“**”表示处理间差异极显著(P<0.01),(x±SD,n=3),下同。
      Note: “*” indicated significant difference among different treatments (P<0.05), “**” indicated extremely significant difference among different treatments (P<0.01), (x±SD,n=3).The same as below.
    下载: 导出CSV

    表 2  不同盐碱胁迫对金丝楸幼苗生理和光合指标的影响

    Table 2.  Effects of different saline-alkali stress on physiological and photosynthetic indicators of C. bungei seedlings

    处理
    Treatment
    相对电导率/%
    Relative conductivity
    可溶性糖含量/
    (ug·g−1)
    Contents of
    soluble sugar
    脯氨酸含量/
    (ug·g−1)
    Contents of Pro
    SOD酶活性/
    (mg·g−1)
    SOD enzymatic activity
    丙二醛含量/
    (ug·g−1)
    Contents of MDA
    叶绿素总量
    Contents of total chlorophyll
    光合速率
    photosynthetic rate
    CK06.35 + 0.2314.68 + 0.640.028 + 0.0024.88 + 0.137.23 + 0.3437.40 + 2.039.63 + 0.94
    NaCl胁迫
    NaCl stress
    5018.04 + 1.42**16.36 + 0.920.034 + 0.0035.27 + 0.788.61 + 0.5438.04 + 2.2211.80 + 1.08*
    10022.27 + 1.78**18.92 + 1.05*0.049 + 0.002*7.42 + 0.63*9.11 + 0.75*40.20 + 3.2713.07 + 0.42*
    15034.59 + 2.67**19.24 + 0.75*0.041 + 0.001*6.51 + 0.52*9.43 + 0.66*36.30 + 2.8910.40 + 1.14
    20048.82 + 3.94**19.03 + 0.84*0.037 + 0.002*4.33 + 0.359.06 + 0.49*35.40 + 0.977.40 + 0.46*
    Na2CO3胁迫
    Na2CO3 stress
    5025.33 + 3.65**15.93 + 0.560.039 + 0.001*5.92 + 0.577.82 + 0.4838.50 + 1.2713.40 + 0.44*
    10036.25 + 2.61**19.87 + 0.90*0.052 + 0.003*8.64 + 0.79*8.49 + 0.5141.08 + 3.07*13.67 + 1.10
    15043.73 + 4.82**17.54 + 0.750.044 + 0.002*7.58 + 0.63*9.82 + 0.75*37.88 + 1.5811.70 + 0.95
    20052.57 + 4.21**17.06 + 0.790.035 + 0.0035.58 + 0.759.55 + 0.81*33.70 + 0.44*9.43 + 0.49
    混合盐碱胁迫
    Mixed saline-
    alkali stress
    5017.68 + 1.23**16.03 + 0.670.036 + 0.0015.23 + 0.378.05 + 0.6134.28 + 1.2914.67 + 0.78**
    10024.54 + 2.01**17.85 + 0.830.047 + 0.003*7.94 + 0.51*8.74 + 0.5641.02 + 2.7314.37 + 1.99*
    15033.78 + 3.24**19.08 + 0.91*0.043 + 0.002*7.16 + 0.69*9.27 + 0.94*36.72 + 1.6010.03 + 0.57
    20050.36 + 4.27**18.91 + 0.76*0.032 + 0.0015.19 + 0.579.46 + 0.83*35.84 + 1.468.57 + 0.76
    下载: 导出CSV
  • [1] 张建锋, 张旭东, 周金星, 等. 世界盐碱地资源及其改良利用的基本措施[J]. 水土保持研究, 2005, 12(6):32-34.

    [2] 范国强, 翟晓巧, 毕会涛, 等. 金丝楸的离体培养和植株再生[J]. 植物生理学通讯, 2002, 38(4):359.

    [3] 辛建华, 许庆标, 刘建军, 等. 不同树龄楸树种质资源的嫁接成活影响因素[J]. 山东林业科技, 2021, 51(4):71-74. doi: 10.3969/j.issn.1002-2724.2021.04.016

    [4] 张 蕊, 石 娜, 孙玲凌, 等. 楸树无糖组织培养快繁技术初探[J]. 江苏农业科学, 2022, 50(1):45-50.

    [5] 王改萍, 王良桂, 王晓聪, 等. 楸树嫩枝扦插生根发育及根系特征分析[J]. 南京林业大学学报:自然科学版, 2020, 44(6):94-102.

    [6] 岑云昕, 刘 佳, 陈发菊, 等. 农杆菌介导的楸树遗传转化体系[J]. 林业科学, 2021, 57(8):195-204. doi: 10.11707/j.1001-7488.20210820

    [7] 杨丹丹, 马玲玲, 李 亚, 等. 楸树EST-SSR标记开发及种质资源的遗传多样性分析[J]. 分子植物育种, 2020, 18(4):1216-1223.

    [8]

    Grabherr M G, Haas B J, Yassour M, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome[J]. Nat Biotechn, 2011, 29(7): 644-652. doi: 10.1038/nbt.1883
    [9]

    Bahieldin A, Atef A, Sabir J S M. RNA-Seq analysis of the wild barley (H. spontaneum) leaf transcriptome under salt stress[J]. Comptes Rendus Biologies, 2015, 338(5): 285-297. doi: 10.1016/j.crvi.2015.03.010
    [10]

    Wang W S, Zhao X Q, Li M, et al. Complex molecular mechanisms underlying seedling salt tolerance in rice revealed by comparative transcriptome and metabolomic profiling[J]. Journal of Experimental Botany, 2016, 67(1): 405-419. doi: 10.1093/jxb/erv476
    [11] 董 明, 再吐尼古丽·库尔班, 吕 芃, 等. 高粱苗期耐盐性转录组分析和基因挖掘[J]. 中国农业科学, 2019, 52(22):3987-4001. doi: 10.3864/j.issn.0578-1752.2019.22.005

    [12] 毛伟兵, 陈发菊, 王长兰, 等. 楸树雄性不育花芽转录组测序及分析[J]. 林业科学, 2017, 53(6):141-150. doi: 10.11707/j.1001-7488.20170617

    [13] 张恩亮, 马玲玲, 杨如同, 等. IBA诱导楸树嫩枝扦插不定根发育的转录组分析[J]. 林业科学, 2018, 54(5):48-61. doi: 10.11707/j.1001-7488.20180506

    [14] 张 丽. 平欧杂种榛抗盐碱生理机制研究及其耐盐性评价[D]. 中国林业科学研究院, 2015.

    [15]

    Li B, Dewey C N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome[J]. BMC Bioinformatics, 2011, 12(1): 323-338. doi: 10.1186/1471-2105-12-323
    [16]

    Marioni J C, Mason C E, Mane S M, et al. RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays[J]. Genome research, 2008, 18(9): 1509-1517.
    [17]

    Jin J, Tian F, Yang D C, et al. PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants[J]. Nucleic acids research, 2017, 45(D1): D1040-D1045. doi: 10.1093/nar/gkw982
    [18]

    Vicente O, Boscaiu M, Naranjo M A, et al. Responses to salt stress in the halophyte Plantago cyassifolia (Plantaginaceae)[J]. Journal of Arid Environments, 2004, 58(4): 463-481. doi: 10.1016/j.jaridenv.2003.12.003
    [19]

    Yang C W, Chong J N, Li C Y, et al. Osmotic adjustment and ion balance traits of an alkali resistant halophyte Kochia sieversiana during adaptation to salt and alkali conditions[J]. Plant Soil, 2007, 294(12): 263-276.
    [20]

    Nasr S M H, Parsakhoo A, Naghavi H, et al. Effect of salt stress on germination and seedling growth of Prosopis juliflora (Sw.)[J]. New Forests, 2012, 43(1): 45-55. doi: 10.1007/s11056-011-9265-9
    [21]

    Megdiche W, Amor N B, Debez A, et al. Salt tolerance of the annual halophyte Cakile maritima as affected by the provenance and the developmental stage[J]. Acta Physiol Plant, 2007, 29(4): 375-384. doi: 10.1007/s11738-007-0047-0
    [22] 李子英, 丛日春, 杨庆山, 等. 盐碱胁迫对柳树幼苗生长和渗透调节物质含量的影响[J]. 生态学报, 2017, 37(24):8511-8517.

    [23]

    Anower M R, Peel M D, Mott I W, et al. Physiological processes associated with salinity tolerance in an alfalfa half-sib family[J]. Journal of Agronomy and Crop Science, 2017, 203(6): 506-518. doi: 10.1111/jac.12221
    [24]

    Li R L, Shi F C, Fukuda K, et al. Effects of salt and alkali stresses on germination, growth, photosynthesis and ion accumulation in alfalfa (Medicago sativa L.)[J]. Soil Science and Plant Nutrition, 2010, 56(5): 725-733. doi: 10.1111/j.1747-0765.2010.00506.x
    [25]

    Yang C W, Xu H H, Wang L L, et al. Comparative effects of salt-stress and alkali-stress on the growth, photosynthesis, solute accumulation, and ion balance of barley plants[J]. Photosynthetica, 2009, 47(1): 79-86. doi: 10.1007/s11099-009-0013-8
    [26]

    Meloni D A, Oliva M A, Martinez C A, et al. Photosynthesis and activity of superoxide dismutase, peroxidase and glutathione reductase in cotton under salt stress[J]. Environmental and Experimental Botany, 2003, 49(1): 69-76. doi: 10.1016/S0098-8472(02)00058-8
    [27] 孙海菁, 王树凤, 陈益泰. 盐胁迫对6个树种的生长及生理指标的影响[J]. 林业科学研究, 2009, 22(3):315-324. doi: 10.3321/j.issn:1001-1498.2009.03.002

    [28] 邓敏捷, 张晓申, 范国强, 等. 四倍体泡桐对盐胁迫生理响应的差异[J]. 中南林业科技大学学报, 2013, 33(11):42-46. doi: 10.3969/j.issn.1673-923X.2013.11.008

    [29]

    Santos C V. Regulation of chlorophyll biosynthesis and degradation by salt stress in sunflower leaves[J]. Scientia Horticulturae, 2004, 103(1): 93-99. doi: 10.1016/j.scienta.2004.04.009
    [30]

    Sicilia A, Santoro D F, Testa G, et al. Transcriptional response of giant reed (Arundo donax L.) low ecotype to long-term salt stress by unigene-based RNAseq[J]. Phytochemistry, 2020, 177: 112436. doi: 10.1016/j.phytochem.2020.112436
    [31] 李 红, 李 波, 邬婷婷, 等. 紫花苜蓿耐苏打盐碱相关基因的转录组学分析[J]. 草地学报, 2019, 27(4):848-859.

    [32] 张 瑞, 张夏燚, 赵 婷, 等. 基于转录组分析垂丝海棠响应盐碱胁迫的分子机制[J]. 园艺学报, 2022, 49(2):237-251.

    [33]

    Khan SA, Li MZ, Wang SM, et al. Revisiting the role of plant transcription factors in the battle against abiotic stress[J]. International Journal of Molecular Sciences, 2018, 19(6): 1634-1663. doi: 10.3390/ijms19061634
    [34]

    Erpen L, Devi HS, Grosser JW, et al. Potential use of the DREB/ERF, MYB, NAC and WRKY transcription factors to improve abiotic and biotic stress in transgenic plants[J]. Plant Cell Tiss Organ Cult, 2018, 132(1): 1-25. doi: 10.1007/s11240-017-1320-6
    [35]

    Li W, Pang S, Lu Z, et al. Function and mechanism of WRKY transcription factors in abiotic stress responses of plants[J]. Plants, 2020, 9(11): 1515-1530. doi: 10.3390/plants9111515
    [36]

    Li H, Gao Y, Xu H, et al. ZmWRKY33, a WRKY maize transcription factor conferring enhanced salt stress tolerances in Arabidopsis[J]. Plant Growth Regul, 2013, 70(3): 207-216. doi: 10.1007/s10725-013-9792-9
    [37]

    Song Y S, Li J L, Sui Y, et al. The sweet Sorghum SbWRKY50 is negatively involved in salt response by regulating ion homeostasis[J]. Plant Molecular Biology, 2020, 102(6): 603-614. doi: 10.1007/s11103-020-00966-4
    [38]

    Zhu H, Jiang Y, Guo Y, et al. A novel salt inducible WRKY transcription factor gene, AhWRKY75, confers salt tolerance in transgenic peanut[J]. Plant Physiology and Biochemistry, 2021, 160: 175-183. doi: 10.1016/j.plaphy.2021.01.014
    [39]

    Puranik S, Sahu P P, Srivastava P S, et al. NAC proteins: Regulation and role in stress tolerance[J]. Trends in Plant Science, 2012, 17(6): 369-381. doi: 10.1016/j.tplants.2012.02.004
    [40]

    Hong Y, Zhang H, Huang L, et al. Overexpression of a stress-responsive NAC transcription factor gene ONAC022 improves drought and salt tolerance in rice[J]. Frontiers in Plant Science, 2016, 7(4): 1-19.
    [41]

    Mao X and Chen S. Novel NAC transcription factor TaNAC67 confers enhanced multi-abiotic stress tolerances in Arabidopsis[J]. PLoS ONE, 2014, 9(1): e84359-84374. doi: 10.1371/journal.pone.0084359
    [42]

    Li S, Wang N, Ji D, et al. A GmSIN1/GmNCED3s/GmRbohBs Feed-forward loop acts as a signal amplifier that regulates root growth in soybean exposed to salt stress[J]. The Plant Cell, 2019, 31: 2107-2130. doi: 10.1105/tpc.18.00662
    [43]

    Banerjee A, Roychoudhury A. Abscisic-acid-dependent basic leucine zipper (bZIP) transcription factors in plant abiotic stress[J]. Protoplasma, 2015, 254(1): 1-14.
    [44]

    Liang C, Meng Z, Meng Z, et al. GhABF2, a bZIP transcription factor, confers drought and salinity tolerance in cotton (Gossypium hirsutum L.)[J]. Scientific Reports, 2016, 6: 35040-35054.
    [45]

    Wang C, Lu G, Hao Y, et al. ABP9, a maize bZIP transcription factor, enhances tolerance to salt and drought in transgenic cotton[J]. Planta, 2017, 246(6): 1-19.
    [46]

    Yang Y, Yu T F, Ma J, et al. The soybean bZIP transcription factor gene GmbZIP2 confers drought and salt resistances in transgenic plants[J]. International Journal of Molecular Sciences, 2020, 21(2): 670-689. doi: 10.3390/ijms21020670
    [47]

    Li J L, Han G L, Sun C F, et al. Research advances of MYB transcription factors in plant stress resistance and breeding[J]. Plant Signaling & Behavior, 2019, 14(8): 1-9.
    [48]

    Li B, Fan R, Guo S, et al. The Arabidopsis MYB transcription factor, MYB111 modulates salt responses by regulating flavonoid biosynthesis[J]. Environmental and Experimental Botany, 2019, 166(2): 807-834.
    [49]

    Wang F, Ren X, Zhang F, et al. A R2R3-type MYB transcription factor gene from soybean, GmMYB12, is involved in flavonoids accumulation and abiotic stress tolerance in transgenic Arabidopsis[J]. Plant Biotechnology Reports, 2019, 13: 219-233. doi: 10.1007/s11816-019-00530-7
    [50]

    Chen H H, Lai L Y, Li L X, et al. AcoMYB4, an Ananas comosus L. MYB transcription factor, functions in osmotic stress through negative regulation of ABA signaling[J]. International Journal of Molecular Sciences, 2020, 21(16): 5727-5742. doi: 10.3390/ijms21165727
    [51]

    Yang T R, Yao S F, Hao L, et al. Wheat bHLH-type transcription factor gene TabHLH1 is crucial in mediating osmotic stresses tolerance through modulating largely the ABA-associated pathway[J]. Plant Cell Reports, 2016, 35(11): 2309-2323. doi: 10.1007/s00299-016-2036-5
    [52]

    Chen H C, Cheng W H, Hong C Y, et al. The transcription factor OsbHLH035 mediates seed germination and enables seedling recovery from salt stress through ABA-dependent and ABA-independent pathways, respectively[J]. Rice, 2018, 11(1): 50-67.
    [53]

    Qiu J R, Huang Z, Xiang X Y, et al. MfbHLH38, a Myrothamnus flabellifolia bHLH transcription factor, confers tolerance to drought and salinity stresses in Arabidopsis[J]. BMC Plant Biology, 2020, 20(1): 542-556. doi: 10.1186/s12870-020-02732-6
  • [1] 田林尹丹丹成铁龙夏新莉尹伟伦 . 盐胁迫下比拉底白刺差异表达基因的转录组分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2020.01.001
    [2] 梁锐涛韩维栋杨少瑕陈蓓蓓 . 无瓣海桑根响应盐胁迫的转录组分析. 林业科学研究,
    [3] 宗建伟张志龙薛可欣王爽爽杨雨华 . 盐碱胁迫下文冠果生理响应机制研究. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2021.03.018
    [4] 周婧李卓蓉吴江婷 . 不同氮形态处理条件下杨树根尖差异表达基因的特征分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2022.02.006
    [5] 张丽贾志国马庆华王贵禧 . 盐碱胁迫对平欧杂种榛生长及叶片内源激素含量的影响. 林业科学研究,
    [6] 李露双董文慧丁兴萃章志远孙春娃蔡函江 . 麻竹笋转录组测序及苦涩味物质合成基因差异表达分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2018.04.006
    [7] 袁长春莫健新袁柳娇余如凤刘倩钟军弟刘锴栋 . 重金属胁迫下白骨壤数字基因表达谱分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2017.02.004
    [8] 张振张含国周宇张磊于宏影张莉 . 红松4个组织的转录组数据分析与次生代谢产物的表达差异初探. 林业科学研究,
    [9] 邓楠史胜青常二梅刘建锋兰倩江泽平 . 基于中麻黄萌发种子转录组的黄酮类化合物合成途径基因的挖掘. 林业科学研究,
    [10] 马庆国王治军徐慧敏李明昆李慧民张俊佩张建国 . 金丝楸优良无性系遗传多样性和亲缘关系的AFLP分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2020.05.018
    [11] 刘宇徐焕文刘桂丰姜静 . 赤霉素GA4+7处理下白桦无性系生长及差异基因表达分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2017.01.025
    [12] 王帅邵芬娟李论芦强邱德有 . 穗花杉的转录组测序及其转录组特性分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2017.05.008
    [13] 饶龙兵杨汉波郭洪英段红平陈益泰 . 基于桤木属转录组测序的SSR分子标记的开发. 林业科学研究,
    [14] 马婧邓楠褚建民纪敬史胜青江泽平成铁龙 . 泡泡刺高通量转录组鉴定及其黄酮类代谢途径初步分析. 林业科学研究,
    [15] 杜明凤丁贵杰 . 基于马尾松干旱转录组的抗旱功能SSR位点分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2018.05.002
    [16] 李太强刘雄芳万友名李正红李钰莹刘秀贤马宏 . 滇东南濒危植物长梗杜鹃转录组微卫星特征分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2017.04.001
    [17] 李珊珊曾艳飞何彩云张建国 . 基于沙棘转录组序列开发EST-SSR分子标记. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2017.01.0010
    [18] 刘雄芳李太强李正红万友名刘秀贤张序安静马宏 . 云南干热河谷地区余甘子转录组分析. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2018.05.001
    [19] 刘丽婷温强黄小春刘琪璟 . 蕨类植物芒萁幼孢子体转录组高通量测序及特征分析. 林业科学研究,
    [20] 封润霞赵婕张苏芳王建军魏建荣刘建凤 . 绒毛白蜡(Fraxinus velutina Torr)韧皮部响应白蜡窄吉丁(Agrilus planipennis Fairmaire)危害的转录组变化. 林业科学研究, doi: 10.13275/j.cnki.lykxyj.2021.01.006
  • 加载中
图(9) / 表(2)
计量
  • 文章访问数:  1490
  • HTML全文浏览量:  390
  • PDF下载量:  10
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-04-14
  • 录用日期:  2022-05-25
  • 网络出版日期:  2022-06-17

金丝楸幼苗响应盐碱胁迫的生理和转录组分析

    通讯作者: 刘彦珍, 1135953615@qq.com
  • 1. 安阳工学院生物与食品工程学院,河南 安阳 455000
  • 2. 滑县园林绿化管理处,河南 安阳 455000
  • 3. 三峡大学生物与制药学院,湖北 宜昌 443000

摘要:  目的 研究不同盐碱胁迫对金丝楸幼苗生长、光合和生理指标的影响,并结合转录组测序分析,探究楸树耐盐碱的生理机制和分子机制。 方法 采用盆栽法对金丝楸幼苗进行不同盐碱胁迫处理,分析其生物量、光合及生理指标对不同盐碱响应的差异,采用Illumina高通量测序技术进行转录组测序,通过生物信息学分析盐碱胁迫对转录水平的影响。 结果 不同盐碱胁迫下,金丝楸幼苗叶片受伤害程度为Na2CO3>混合盐碱>NaCl;新增株高和地径、地上部和根的干质量和鲜质量、生物量、根冠比均受到明显抑制,并随盐碱浓度增加而抑制加强,但生长胁迫指数均随浓度的增加而降低;丙二醛(MDA)含量和相对电导率都随胁迫浓度增加而不同程度上升,超氧化物歧化酶(SOD)活性、可溶性糖含量、脯氨酸(Pro)含量、叶绿素总量和光合速率都呈现先上升后下降的趋势。转录组测序共产生约60.4 Gb原始数据,组装得到55 793个Unigenes,其中29 534(52.93%)个Unigenes获得了注释;通过差异表达基因(DEGs)分析,3个比较组(CK vs NaCl,CK vs Na2CO3 和 CK vs 混合盐碱)分别筛选出1 779、2 835和4 059 个DEGs;DEGs GO富集分析表明,膜的整体成分、膜的内在成分、催化活性、类异戊二烯代谢和合成过程、氧化还原酶活性等条目被显著富集;DEGs KEGG分析表明,苯丙素生物合成、淀粉和蔗糖代谢、植物激素信号转导、萜类主干生物合成和精氨酸代谢等通路被显著富集;此外,在DEGs中鉴定的bHLH、ERF、MYB-related、NAC、C2H2、WRKY、MYB和bZIP转录因子家族成员最多。 结论 金丝楸主要通过积累可溶性糖和Pro,提高SOD酶活和光合作用来抵御盐碱胁迫,但都呈现“低促高抑”的现象,说明其具有一定阈值。金丝楸通过调节膜成分、催化活性、类异戊二烯代谢和生物合成过程、苯丙素生物合成、淀粉和蔗糖代谢、植物激素信号转导等生物过程和代谢途径,并结合有关转录因子共同响应盐碱胁迫。本研究为深入研究楸树耐盐碱生理机制和分子机制提供科学的理论依据。

English Abstract

  • 当前,全球土地盐碱化是人类面临的主要生态危机之一,给农业生产带来了巨大损失,严重威胁人类生存。据不完全统计,全球盐碱地面积约为9.54 × 109 hm2,约占可耕地面积的10%,而我国盐碱地面积约有9.87 × 108 hm2,位居世界第三,大部分为土壤盐化和碱化同时存在的复合盐碱地[1]。目前,科学合理地开发和利用盐碱地已经成为我国农业和林业发展过程中面临的严峻挑战。

    楸树(Catalpa bungei C. A. Mey.)是我国特有的珍贵优质用材和园林观赏树种,属于紫葳科(Bignoniaceae)梓属(Catalpa)落叶乔木,其材质好、用途广、经济价值高,自古就有“木王”的美誉。金丝楸自然分布于河南、山东和安徽等省,是楸树的一个地理变种,具有特殊金色纹理、生长快和干型通直的优点,近年来得到大面积推广[2]。目前,楸树的研究主要集中在嫁接繁育[3]、组培快繁[4]、扦插繁殖[5]、遗传转化[6]、遗传多样性分析[7]等方面,关于楸树耐盐碱的研究鲜有报道。认识和了解楸树耐盐碱机制,对指导楸树抗盐碱遗传改良和利用盐碱地具有重要的理论价值和现实意义。

    近年来,转录组测序技术(RNA-Seq)发展迅速,该技术无须了解物种基因组信息,弥补了非模式植物转录组研究中缺乏基因组信息的不足[8]。利用转录组解析植物的抗逆机制已经成为逆境研究的热点之一。目前,利用RNA-Seq已进行了大麦[9]、水稻[10]、高粱[11]等盐胁迫的转录组分析,对这些植物响应盐胁迫机理的探讨,加深了人们对植物耐盐机制的认识。但是,通过转录组分析楸树耐盐碱分子机制的研究尚未见报道,仅有通过转录组分析楸树雄性不育和不定根发育分子机制的极少数报道[12-13]

    本研究以金丝楸幼苗为试验材料,分析不同盐碱胁迫对金丝楸生长情况、光合及生理特性的影响,探究其耐盐碱生理机制。进一步利用RNA-Seq,对处理组和对照组进行DEGs分析,在转录水平上对楸树耐盐碱机制进行探索,以期为楸树应答盐碱胁迫的分子机制研究提供理论参考,同时也为重要基因的克隆及其功能研究等奠定基础。

    • 试验于安阳工学院温室中进行,2020年6月,选择生长健壮、形态一致、根系完整的金丝楸幼苗种于盆中,每盆种3株,共种40盆。正常栽培管理,60 d后进行盐碱胁迫处理。分别设置中性盐A(NaCl)、碱性盐B(Na2CO3)、混合盐碱C(NaCl∶Na2CO3=1∶1)3种盐碱胁迫处理,总盐浓度设50、100、150和200 mmol·L−1 4个梯度,A1-A4、B1-B4、C1-C4分别代表从低到高的NaCl、Na2CO3、混合盐碱胁迫的4个浓度,共12个处理浓度。2020年8月,将试验材料随机分成13组,每组3盆,每盆浇2.0 L相应浓度的处理液,对照组浇等体积的清水。期间正常栽培管理,胁迫30 d后进行生长、光合及生理指标的测定。胁迫7 d后取长势一致且大小相近的4组(CK、A2、B2和C2)叶片用于转录组测序,生物学重复2次。

    • 各处理组在盐碱胁迫后定期(14 d、21 d和28 d)对植株和叶片进行危害调查,根据胁迫伤害症状,设置分级标准(图1)如下:

      图  1  金丝楸不同伤害等级的叶片盐碱胁迫症状

      Figure 1.  The symptom of different damage level on leaves of C. bungei under different saline-alkali stress

      0级:枝条和叶片未见盐碱胁迫危害症状;

      1级:有极少部分叶尖、叶缘和叶脉变黄,枝条外观正常;

      2级:有少部分的叶尖、叶缘焦枯,枝条外观正常;

      3级:约有1/4叶片有叶尖、叶缘焦枯和落叶现象;

      4级:约有1/2叶片枝枯、叶落直至死亡。

    • 对照组和处理组均选择3株进行编号,处理前用卷尺和游标卡尺测定株高和地径,胁迫30 d后,再次测定上述指标,计算净生长量。

    • 胁迫30 d后,取对照组和处理组编号植株的根和地上部称其鲜质量,80 ℃烘干至恒质量,称其干质量,计算生物量(根干质量 + 地上部干质量)、根冠比(根干质量/地上部干质量)及生长胁迫指数(处理组植株干质量/对照组植株干质量)[14]

    • 用SPAD-502Plus叶绿素仪测定叶绿素总量,EBA-PE1001高光效仪测定光合速率。

    • 用电导仪法测定相对电导率,用硫代巴比妥酸法测定MDA含量,用氮蓝四唑法测定SOD酶活,用蒽酮显色法测定可溶性糖含量,用酸性茚三酮法测定Pro含量[14]

    • 样品经过RNA抽提、纯化、建库后,基于Illumina HiSeq测序平台进行双末端测序。

    • 对原始数据进行过滤,将带接头、小于50 bp、平均质量在Q20以下的Reads去除,得到高质量序列(Clean Reads),使用Trinity对Clean Reads进行De novo拼接得到转录本序列,挑选最长的转录本作为Unigene[8]

    • 用Blast对Unigenes进行NR、GO、KEGG、eggNOG、SwissProt和Pfam注释;对所有Unigenes进行GO功能分析和KEGG代谢途径分析。用转录组表达定量软件RSEM[15],计算每个基因的FPKM值[16],来代表每个基因的表达水平。

    • 用DESeq对基因表达进行差异分析,筛选DEGs条件为:表达差异倍数|log2FoldChange|>1,显著性P-value<0.05。然后将得到的DEGs再进行GO和KEGG分析。

    • 将蛋白质序列与Plant TFDB数据库比较[17],分析得到转录因子及其所属家族信息。

    • 不同盐碱胁迫后,金丝楸均出现盐碱胁迫危害症状。Na2CO3胁迫对叶片伤害较大,随胁迫浓度和时间的增加,伤害等级增大,胁迫处理14 d,150 mmol·L−1时伤害等级达到3级,而200 mmol·L−1时达到4级;而NaCl胁迫对叶片伤害较小,胁迫处理21 d,150 mmol·L−1时伤害等级达到2级,200 mmol·L−1时达到4级;混合盐碱胁迫下叶片伤害介于Na2CO3和NaCl之间。

      不同盐碱胁迫下,金丝楸的生长情况受到不同程度的影响(表1)。NaCl胁迫下,50 mmol·L−1 时新增株高被极显著的抑制,且浓度越高抑制作用越显著;Na2CO3胁迫对新增株高的影响较小,在50 mmol·L−1时新增株高显著降低,且随浓度增加株高极显著降低;混合盐碱胁迫对新增株高抑制作用最轻,但也都达到极显著水平。此外,低浓度NaCl和Na2CO3胁迫对地径有增粗作用,在50 mmol·L−1NaCl胁迫下,新增地径显著增加,而在150和200 mmol·L−1不同盐碱胁迫下,新增地径都极显著减小。

      表 1  不同盐碱胁迫对金丝楸幼苗生长的影响

      Table 1.  Effects of different saline-alkali stress on the growth of C. bungei seedlings

      处理
      Treatment
      新增株高
      Net growth of plant height/
      cm
      新增地径
      Net growth of stem diameter/
      mm
      地上部overground part
      鲜质量
      Fresh weight/g
      干质量
      Dry weight/g
      CK018.68 + 0.621.23 + 0.0924.60 + 1.145.70 + 0.19
      NaCl胁迫
      NaCl stress
      5012.58 + 0.62**1.50 + 0.13*27.06 + 1.086.75 + 0.25*
      10010.00 + 0.53**0.98 + 0.09*15.94 + 0.07**3.31 + 0.05**
      1508.14 + 0.36**0.66 + 0.10**10.13 + 0.14**2.14 + 0.09**
      2005.92 + 0.79**0.40 + 0.07**6.92 + 0.18**1.54 + 0.08**
      Na2CO3胁迫
      Na2CO3 stress
      5017.64 + 0.98*1.26 + 0.1018.87 + 0.62**4.68 + 0.18**
      10012.68 + 0.58**1.03 + 0.0911.20 + 0.84**3.14 + 0.22**
      15010.48 + 1.36**0.57 + 0.09**6.87 + 0.24**1.73 + 0.09**
      2006.58 + 0.74**0.37 + 0.05**4.80 + 0.10**1.12 + 0.05**
      混合盐碱胁迫
      Mixed saline-alkali
      stress
      5016.74 + 0.77**1.18 + 0.0835.18 + 0.42**9.13 + 0.20**
      10013.58 + 1.06**0.84 + 0.1315.56 + 0.42**3.35 + 0.29**
      15010.10 + 0.62**0.59 + 0.08**13.09 + 0.47**2.67 + 0.30**
      2009.50 + 1.22**0.37 + 0.08** 6.76 + 0.33** 1.40 + 0.05**
      处理
      Treatment
      根 root生物量
      Biomass/g
      根冠比
      Root-shoot ratio
      生长胁迫指数
      Growth stress index
      鲜质量
      Fresh weight/g
      干质量
      Dry weight/g
      CK08.99 + 0.151.28 + 0.056.97 + 0.240.22 + 0.01
      NaCl胁迫
      NaCl stress
      508.00 + 0.16*0.90 + 0.03*7.65 + 0.270.13 + 0.01**1.10 + 0.04
      1004.08 + 0.09**0.48 + 0.03**3.97 + 0.07**0.15 + 0.01**0.54 + 0.01
      1503.29 + 0.08**0.33 + 0.03**2.47 + 0.12**0.15 + 0.01**0.35 + 0.02
      2002.57 + 0.07**0.26 + 0.04**1.79 + 0.11**0.17 + 0.02*0.26 + 0.02
      Na2CO3胁迫
      Na2CO3 stress
      504.90 + 0.27**0.74 + 0.09*5.42 + 0.20**0.16 + 0.02*0.78 + 0.03
      1002.55 + 0.43**0.31 + 0.07**3.44 + 0.29**0.10 + 0.01**0.49 + 0.04
      1502.51 + 0.45**0.27 + 0.04**2.00 + 0.12**0.16 + 0.01*0.29 + 0.02
      2001.31 + 0.19**0.13 + 0.02**1.25 + 0.07**0.12 + 0.01**0.18 + 0.01
      混合盐碱胁迫
      Mixed saline-alkali
      stress
      5012.44 + 0.90*2.02 + 0.18*11.15 + 0.37**0.22 + 0.021.60 + 0.05
      1006.09 + 0.92*0.65 + 0.06**4.00 + 0.34**0.19 + 0.02**0.57 + 0.05
      1502.62 + 0.56**0.28 + 0.04**2.95 + 0.33**0.10 + 0.01**0.42 + 0.05
      2001.86 + 0.27**0.22 + 0.04**1.62 + 0.08**0.15 + 0.02*0.23 + 0.01
      注:“*”表示处理间差异显著(P<0.05),“**”表示处理间差异极显著(P<0.01),(x±SD,n=3),下同。
        Note: “*” indicated significant difference among different treatments (P<0.05), “**” indicated extremely significant difference among different treatments (P<0.01), (x±SD,n=3).The same as below.

      植物在盐碱胁迫下生物量的变化可直观反映生长状况(表1)。与对照组相比,NaCl和混合盐碱处理组,50 mmol·L−1时地上部干质量和鲜质量都有增加,但随浓度增加而下降;Na2CO3处理组地上部干质量和鲜质量随着浓度增加而下降。根干质量和鲜质量只在50 mmol·L−1混合盐碱处理组有增加,其它处理组都随浓度增加而下降。生物量的变化趋势与地上部干质量的变化基本一致。各处理组的根冠比均显著降低。不同盐碱胁迫下生长胁迫指数均随胁迫浓度的增加而降低,200 mmol·L−1时各处理组生长胁迫指数分别为0.26、0.18、0.23,可以看出Na2CO3胁迫降幅较大,NaCl和混合盐碱胁迫对生长胁迫指数影响较小。这些结果表明,Na2CO3和混合盐碱胁迫对金丝楸的生长发育抑制更严重,与对叶片危害的影响一致。

    • 不同盐碱胁迫对金丝楸生理和光合指标的影响如表2所示。各处理组中,MDA含量和相对电导率都随处理浓度的增加而显著上升。NaCl和混合盐碱处理组中,可溶性糖含量在50~150 mmol·L−1时逐渐积累,而在200 mmol·L−1时小幅下降;而Na2CO3处理组中,在100 mmol·L−1时达到最大值,随浓度增加而下降。各处理组中,Pro含量都明显增加,都在100 mmol·L−1时达到最大值,随处理浓度的增加而逐渐下降,但仍明显高于对照组。各处理组中,SOD酶活性在中低浓度都有所增加,其中100~150 mmol·L−1时显著增加,200 mmol·L−1时下降。各处理组中,低浓度时叶绿素总量和光合速率增加,都在100 mmol·L−1时达到最大值,随处理浓度的增加而逐渐下降。

      表 2  不同盐碱胁迫对金丝楸幼苗生理和光合指标的影响

      Table 2.  Effects of different saline-alkali stress on physiological and photosynthetic indicators of C. bungei seedlings

      处理
      Treatment
      相对电导率/%
      Relative conductivity
      可溶性糖含量/
      (ug·g−1)
      Contents of
      soluble sugar
      脯氨酸含量/
      (ug·g−1)
      Contents of Pro
      SOD酶活性/
      (mg·g−1)
      SOD enzymatic activity
      丙二醛含量/
      (ug·g−1)
      Contents of MDA
      叶绿素总量
      Contents of total chlorophyll
      光合速率
      photosynthetic rate
      CK06.35 + 0.2314.68 + 0.640.028 + 0.0024.88 + 0.137.23 + 0.3437.40 + 2.039.63 + 0.94
      NaCl胁迫
      NaCl stress
      5018.04 + 1.42**16.36 + 0.920.034 + 0.0035.27 + 0.788.61 + 0.5438.04 + 2.2211.80 + 1.08*
      10022.27 + 1.78**18.92 + 1.05*0.049 + 0.002*7.42 + 0.63*9.11 + 0.75*40.20 + 3.2713.07 + 0.42*
      15034.59 + 2.67**19.24 + 0.75*0.041 + 0.001*6.51 + 0.52*9.43 + 0.66*36.30 + 2.8910.40 + 1.14
      20048.82 + 3.94**19.03 + 0.84*0.037 + 0.002*4.33 + 0.359.06 + 0.49*35.40 + 0.977.40 + 0.46*
      Na2CO3胁迫
      Na2CO3 stress
      5025.33 + 3.65**15.93 + 0.560.039 + 0.001*5.92 + 0.577.82 + 0.4838.50 + 1.2713.40 + 0.44*
      10036.25 + 2.61**19.87 + 0.90*0.052 + 0.003*8.64 + 0.79*8.49 + 0.5141.08 + 3.07*13.67 + 1.10
      15043.73 + 4.82**17.54 + 0.750.044 + 0.002*7.58 + 0.63*9.82 + 0.75*37.88 + 1.5811.70 + 0.95
      20052.57 + 4.21**17.06 + 0.790.035 + 0.0035.58 + 0.759.55 + 0.81*33.70 + 0.44*9.43 + 0.49
      混合盐碱胁迫
      Mixed saline-
      alkali stress
      5017.68 + 1.23**16.03 + 0.670.036 + 0.0015.23 + 0.378.05 + 0.6134.28 + 1.2914.67 + 0.78**
      10024.54 + 2.01**17.85 + 0.830.047 + 0.003*7.94 + 0.51*8.74 + 0.5641.02 + 2.7314.37 + 1.99*
      15033.78 + 3.24**19.08 + 0.91*0.043 + 0.002*7.16 + 0.69*9.27 + 0.94*36.72 + 1.6010.03 + 0.57
      20050.36 + 4.27**18.91 + 0.76*0.032 + 0.0015.19 + 0.579.46 + 0.83*35.84 + 1.468.57 + 0.76
    • 对试验中的8个样品测序后分别获得了约50 000 000条原始序列,共测得约60.4 Gb原始数据。过滤后Clean datas约占原始序列数据的93%;碱基数约占原始碱基数的93%;Q20和Q30值分别在97%和93%以上。混合2个样品的Clean datas,用Trinity软件进行De novo组装,共得到171 339个Transcripts,55 793个Unigenes。Transcript和Unigene序列平均长度分别为1 768.34 bp和1 275.65 bp,N50分别为2 597 bp和2 263 bp,GC含量分别为39.59%和38.72%。

    • 对55 793个Unigenes进行基因功能注释,在NR、GO、KEGG、Pfam、eggNOG和Swiss-Prot数据库中分别注释28 722个(51.48%)、13 185个(23.63%)、11967个(21.45%)、15 837个(28.39%)、27 242个(48.83%)和21 694个(38.88%),其中,29 534(52.93%)个Unigenes最少被1个数据库注释,5 124(9.18%)个Unigenes被6个数据库注释。28 722个Unigenes在NR数据库中获得注释,其中5.8% Unigenes与之完全匹配,15.94% Unigenes不完全匹配。其中,相似度在95%~100%、80%~95%、60%~80%、40%~60%和≤40%的Unigenes比例分别为5.98%、32.12%、35.43%、20.45%和6.02%。被注释到芝麻(Sesamum indicum L.)、地黄猴面花(Erythranthe guttata DC.)、施蒴苣苔(Dorcoceras hygrometricum B.)、葡萄(Vitis vinifera L.)、咖啡(Coffea canephora PF.)、大麦(Hordeum vulgare L.)、可可(Theobroma cacao L.)和其他物种Unigenes比例分别为55.95%、14.67%、3.35%、2.09%、0.97%、0.91%、0.76%和21.31%。

      基因GO注释后,将成功注释的基因依照GO的三个大类进行分类。Unigenes富集最多的是代谢过程和细胞过程、细胞膜和细胞、结合和催化活性。KEGG 代谢通路上,基因注释最多的分别为碳水化合物代谢、翻译、信号分解代谢、运输分解代谢和内分泌系统。eggNOG比对注释,基因注释最多目录分别为功能未知、仅通用功能预测、信号传导机制和翻译后修饰、蛋白质周转、伴侣。

    • 为了进一步了解楸树对不同盐碱胁迫反应的转录组特征,分析比较了对照组和3种处理组的DEGs。对照组与NaCl处理组(A vs C)产生1 799个DEGs(上调844个,下调935个),对照组与Na2CO3处理组(A vs D)产生2 835个DEGs(上调1 326个,下调1 509个),对照组与混合盐碱处理组(A vs E)产生4 059个DEGs(上调2 295个,下调1 764个),说明这些DEGs与不同盐碱胁迫应答密切相关。根据差异分析的结果,维恩图进一步展示了各比较组间特异DEGs的个数,以及比较组间的重叠关系。A vs C与A vs D之间有1 036个共有DEGs,A vs D与A vs E之间有1 530个共有DEGs,A vs C与A vs E之间有963个共有DEGs,三个比较组之间有717个共有DEGs。A vs C有497个特异DEGs,A vs D有986个特异DEGs,A vs E有2 283个特异DEGs(图2)。说明共有DEGs可能由不同的盐碱胁迫共同引起的,特异DEGs可能响应不同的盐碱胁迫反应。

      图  2  特异和共有差异表达基因的维恩图分析

      Figure 2.  Venn diagrams analysis of specific and shared DEGs

    • 对DEGs进一步GO富集分析发现,在细胞组分、分子功能和生物过程方面存在差异。在细胞组分方面,膜的整体成分、膜的内在成分、膜等在不同盐碱胁迫下被显著富集。在分子功能方面,催化活性、氧化还原酶活性、信号传感器活性、单加氧酶活性、水解酶活性、水解o-糖基化合物等在不同盐碱胁迫下被显著富集。在生物过程方面,类异戊二烯代谢和生物合成过程、防卫反应、碳水化合物代谢过程等在不同盐碱胁迫下被显著富集(图3~5)。

      图  3  DEGs GO富集分析(A Vs C)

      Figure 3.  GO enrichment analysis of the DEGs (A Vs C)

      图  4  DEGs GO富集分析(A Vs D)

      Figure 4.  GO enrichment analysis of the DEGs (A Vs D)

      图  5  DEGs GO富集分析(A Vs E)

      Figure 5.  GO enrichment analysis of the DEGs (A Vs E)

      通过对DEGs进一步KEGG分析发现,DEGs主要集中在代谢途径。NaCl胁迫下显著富集的通路为:萜类主干生物合成、苯丙素的生物合成、植物激素信号转导、光合生物中的碳固定(图6);Na2CO3胁迫下显著富集的通路为:淀粉和蔗糖代谢、苯丙素的生物合成、氰氨基酸代谢、精氨酸生物合成(图7);混合盐碱胁迫下显著富集的通路为:淀粉和蔗糖代谢、萜类主干生物合成、苯丙素的生物合成、植物激素信号转导(图8)。

      图  6  DEGs KEGG富集分析(A Vs C)

      Figure 6.  KEGG enrichment analysis of the DEGs (A Vs C)

      图  7  DEGs KEGG富集分析(A Vs D)

      Figure 7.  KEGG enrichment analysis of the DEGs (A Vs D)

      图  8  DEGs KEGG富集分析(A Vs E)

      Figure 8.  KEGG enrichment analysis of the DEGs (A Vs E)

    • 植物转录因子广泛参与多种生物过程,并在转录水平上调控基因的表达,是许多逆境胁迫响应基因的重要调控因子,在植物对逆境胁迫的应答过程中发挥重要作用。为了进一步研究楸树中可能参与盐碱胁迫的转录因子,我们共鉴定出6 725个差异表达转录因子,分属于57个不同的转录因子家族,其中bHLH、ERF、MYB-related、NAC、C2H2、WRKY、MYB和bZIP转录因子家族成员最多。

    • 盐碱胁迫对植物最直接的影响就是抑制其正常的生长发育[18]。本研究发现,Na2CO3胁迫对金丝楸叶片伤害程度最大,NaCl最小,与前人研究结果一致[19]。不同盐碱胁迫导致了金丝楸新增株高和地径、生物量、根冠比不同程度的降低,并且随着胁迫浓度的增加而加强(表1),表明其对不同盐碱胁迫的忍耐能力存在差异,这与其它植物中的研究结果相同[20]。分析对新增株高的影响程度依次为NaCl > Na2CO3 > 混合盐碱,三种胁迫的中高浓度都极显著抑制新增地径的增长(表1)。对地上部和根的干质量以及生物量积累的影响程度依次为Na2CO3 > NaCl > 混合盐碱,但50 mmol·L−1时NaCl和Na2CO3胁迫对地径有增粗作用,其中NaCl胁迫显著增加,50 mmol·L−1时混合盐碱胁迫对地上部和根的干质量和鲜质量都极显著增加(表1),其它植物的研究也得到类似结果[21]

    • 盐碱胁迫还会影响植物的生理代谢和光合作用。MDA是植物体内膜脂过氧化的产物,常作为对逆境反应强弱的指标。本研究发现,不同盐碱胁迫后MDA都不同程度增加,其中在中高浓度胁迫后显著增加(表2),该结果与前人研究结果一致[22]。Pro能调节渗透势和改善细胞内环境,可溶性糖能维持细胞渗透势和提供碳源,它们是植物通过渗透调节作用对逆境胁迫的一种生理响应[23]。本研究发现,不同盐碱胁迫后,Pro和可溶性糖含量都显著上升,在100或150 mmol·L−1时达到最大值,浓度继续升高后会小幅降低,但仍高于对照组(表2),这一结果与之前关于盐碱胁迫对Pro和可溶性糖含量影响的研究结果相吻合[24-25]。说明在一定浓度范围内楸树可以通过增加体内Pro和可溶性糖的含量,降低细胞渗透势,超过这个范围,Pro和可溶性糖合成受到抑制。盐碱胁迫会使植物产生大量的活性氧自由基,破坏细胞膜系统,抗氧化酶活性的升高可增强自由基清除能力,减轻伤害[26]。本研究发现,不同盐碱胁迫下SOD都随浓度增加呈先升后降的趋势(表2),这与孙海菁等之前发现的情况一致[27]。说明当盐碱胁迫超出阈值后,过量的自由基对植物造成损伤,SOD酶活降低。叶绿素是植物光合作用的物质基础,是反映光合能力的重要指标。本研究发现,在低浓度盐碱胁迫下叶绿素总量升高,光合速率显著提高(表2),说明低浓度胁迫能够促进叶绿素的合成和光合作用,而在高浓度胁迫下叶绿素总量和光合速率下降。该结果和地上部生物量的结果基本一致。该结果与前人研究结果不完全一致[28-29],可能与物种差异有关,有待我们进一步研究。

    • RNA-Seq是研究植物非生物胁迫分子机制的良好手段[8,30],楸树作为非模式植物,缺乏基因组背景,RNA-Seq为从转录水平研究楸树耐盐碱分子机制、基因挖掘及相关代谢通路的发现提供了可能。本研究分别对NaCl、Na2CO3和混合盐碱处理组及对照组进行转录组测序及分析。DEGs GO富集结果表明,在细胞组分方面,膜的整体成分、内在成分和膜在不同盐碱胁迫下均显著富集,说明膜及其成分在不同盐碱胁迫下都发挥积极作用。在分子功能方面,富集结果相差较大,催化活性和氧化还原酶活性主要响应NaCl胁迫;信号传感器活性和单加氧酶活性主要响应Na2CO3胁迫;催化活性和水解o-糖基化合物主要响应混合盐碱胁迫,说明不同分子功能响应不同的盐碱胁迫。在生物过程方面,类异戊二烯代谢和生物合成过程在不同盐碱胁迫下都显著富集,说明类异戊二烯对植物耐盐碱非常重要;二甲基烯丙基二磷酸生物合成过程、防卫反应和碳水化合物代谢过程分别单独显著富集(图3-5)。DEGs KEGG分析结果表明,苯丙素生物合成在不同盐碱胁迫中均显著富集,植物激素信号转导和萜类主干生物合成在NaCl和混合盐碱胁迫中都显著富集,淀粉和蔗糖代谢在Na2CO3和混合盐碱胁迫中都显著富集,还有光合生物碳固定、精氨酸代谢和生物合成分别在NaCl、Na2CO3胁迫中单独显著富集(图6-8)。说明苯丙素生物合成、植物激素信号转导、淀粉和蔗糖代谢、萜类主干生物合成、类异戊二烯通路在应答盐碱胁迫时至关重要,与前人转录组测序分析结果一致[11,31-32]。以上结果暗示这些功能与途径可能在楸树响应不同盐碱胁迫中起着重要作用,这些注释信息也为进一步挖掘楸树抗不同盐碱基因提供了数据资源。

    • 转录因子在植物对多种逆境胁迫的应答过程中发挥重要作用[33]。目前,已报道的响应盐碱胁迫的转录因子有WRKY、NAC、bZIP、MYB、bHLH等家族成员[34]。WRKY是一类锌指蛋白家族,对植物抗逆有着重要作用[35]。玉米ZmWRKY33、高粱SbWRKY50和花生AhWRKY75等WRKY转录因子对耐盐性起作用[36-38]。NAC是植物特有的转录因子家族,在植物发育、逆境应答等过程中发挥多重作用[39]。水稻ONAC022、小麦TaNAC67和大豆GmSIN1等NAC转录因子与耐盐能力有关[40-42]。bZIP转录因子参与依赖ABA信号的途径,调控相关抗盐碱基因的表达[43]。棉花GhABF2、玉米ABP9和大豆GmbZIP2等bZIP转录因子控制盐碱耐性[44-46]。MYB转录因子广泛参与植物的生长发育、逆境响应和激素信号转导[47]。拟南芥MYB111、大豆GmMYB12和菠萝AcoMYB4等MYB转录因子调控盐胁迫响应基因[48-50]。bHLH转录因子同样在非生物胁迫应答中起着重要作用。小麦TabHLH1、水稻OsbHLH035和杨梅MfbHLH38等bHLH转录因子应答盐胁迫反应[51-53]。本研究鉴定出bHLH、ERF、MYB-related、NAC、C2H2、WRKY、MYB和bZIP转录因子家族成员最多,这些转录因子与不同盐碱耐性相关,说明楸树耐盐碱机制是由多种转录因子共同调控的。

    • 本研究通过对不同盐碱胁迫下金丝楸幼苗生长情况、生理指标和光合特性进行测定及分析,新增株高和地径、生物量、根冠比等都明显受到抑制,生长情况受抑制程度为Na2CO3 > 混合盐碱 >NaCl。金丝楸幼苗主要通过积累渗透调节物质、提高抗氧化酶系统和光合作用来共同抵御盐碱胁迫,但都呈现“低促高抑”的现象,具有一定阈值。通过对对照组和3种处理组叶片进行转录组测序及分析,3个比较组分别筛选出1 779、2 835和4 059个DEGs。通过DEGs GO功能和KEGG通路分析,理清了DEGs富集的生物过程与代谢通路。金丝楸通过调节膜成分、催化活性、类异戊二烯代谢和生物合成过程、苯丙素生物合成、淀粉和蔗糖代谢、植物激素信号转导等生物过程和代谢途径,并结合相关转录因子共同响应盐碱胁迫(图9)。综上所述,本研究的首次尝试为全面解析楸树耐盐碱机制研究提供理论参考,同时也为后期重要基因克隆及其功能分析等奠定基础。

      图  9  结论示意图

      Figure 9.  Schematic diagram of conclusion

参考文献 (53)

目录

    /

    返回文章
    返回