留言板

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

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

闽楠PLR基因家族鉴定及响应激素的表达分析

王倩清 张毓婷 张俊红 刘慧 童再康

王倩清, 张毓婷, 张俊红, 刘慧, 童再康. 闽楠PLR基因家族鉴定及响应激素的表达分析[J]. 浙江农林大学学报, 2022, 39(6): 1173-1182. doi: 10.11833/j.issn.2095-0756.20220351
引用本文: 王倩清, 张毓婷, 张俊红, 刘慧, 童再康. 闽楠PLR基因家族鉴定及响应激素的表达分析[J]. 浙江农林大学学报, 2022, 39(6): 1173-1182. doi: 10.11833/j.issn.2095-0756.20220351
WANG Qianqing, ZHANG Yuting, ZHANG Junhong, LIU Hui, TONG Zaikang. Identification of PLR gene family and expression of responsive hormones in Phoebe bournei[J]. Journal of Zhejiang A&F University, 2022, 39(6): 1173-1182. doi: 10.11833/j.issn.2095-0756.20220351
Citation: WANG Qianqing, ZHANG Yuting, ZHANG Junhong, LIU Hui, TONG Zaikang. Identification of PLR gene family and expression of responsive hormones in Phoebe bournei[J]. Journal of Zhejiang A&F University, 2022, 39(6): 1173-1182. doi: 10.11833/j.issn.2095-0756.20220351

闽楠PLR基因家族鉴定及响应激素的表达分析

doi: 10.11833/j.issn.2095-0756.20220351
基金项目: 浙江省农业新品种选育重大专项(2021C2070-10);国家自然科学基金资助项目(32171828)
详细信息
    作者简介: 王倩清(ORCID: 0000-0003-0116-0551),从事闽楠遗传学研究。E-mail: 974317348@qq.com
    通信作者: 童再康(ORCID: 0000-0002-0202-6496),教授,博士,从事林木新品种选育与定向培育研究工作。E-mail: zktong@zafu.edu.cn
  • 中图分类号: S722.3

Identification of PLR gene family and expression of responsive hormones in Phoebe bournei

  • 摘要:   目的  解析闽楠Phoebe bournei木脂素合成途径关键酶基因Pinoresinol-lariciresinol reductase (PLR)特征及不同激素处理下的表达模式,探索闽楠PLR基因家族的生物学功能。  方法  利用生物信息学方法鉴定闽楠基因组中的PLR基因家族成员,分析基因结构、保守基序、染色体定位、系统进化、启动子顺式作用元件和激素响应。  结果  从闽楠全基因组内共鉴定出8个闽楠PLR (PbPLR)成员,不均匀地分布于第1、2、5和10号染色体。系统进化分析表明:PbPLR功能可能与底物立体选择性有关。启动子元件分析发现:PbPLRs启动子区域存在脱落酸、水杨酸和茉莉酸甲酯等激素响应元件;基因表达分析表明:PbPLRs基因的表达具有组织特异性和激素响应特征,其中PbPLRs基因在根中表达量最高,其次茎,叶中表达量较低;3种激素处理时,PbPLR2表达诱导最强。  结论  从闽楠基因组中鉴定出的8个PbPLRs基因,其基因特征和不同成员可能具有不同催化酶立体选择性,成员表达模式多样化,可能受多种激素诱导,具有组织特异性。图6表2参25
  • 图  1  PbPLR基因结构及蛋白保守基序分析

    Figure  1  Analysis of PbPLR gene structure and protein conserved motifs

    图  2  闽楠和多个物种的PLR家族系统发育树

    Figure  2  PLR family phylogenetic tree of P. bournei and several species

    图  3  PbPLR家族基因的顺式作用元件分析

    Figure  3  Analysis of putative cis-elements in the promoter regions of PbPLRs

    图  4  PbPLRs基因在ABA处理下的表达模式

    Figure  4  Expression patterns of PbPLRs under ABA treatment

    图  5  PbPLRs基因在SA处理下的表达模式

    Figure  5  Expression patterns of PbPLRs under SA treatment

    图  6  PbPLRs基因在MeJA处理下的表达模式

    Figure  6  Expression patterns of PbPLRs under MeJA treatment

    表  1  实时荧光定量PCR (RT-qPCR)引物信息

    Table  1.   RT-qPCR primer information

    引物名称引物序列(5′→3′)
    PbPLR1 F: GGCACAAATCTCATCAGCAA; R: GCCAGCTTTTATGGCCTGTA
    PbPLR2 F: TCTGGATTACGTGCAACAGG; R: AGCGTTTGAGATAGTCCTTCG
    PbPLR3 F: AGGTCTGGATTACGCACAGC; R: CGTTTGAGATAGTCCTTCACTCG
    PbPLR4 F: GGAATTCCTTGCTCAAATGAAA; R: TTTCGAAATTTGCCAGACAAC
    PbPLR5 F: GTGGAACAAGAGCCCAAATG; R: CTGGCAAGTGGGTGAAAAAT
    PbPLR6 F: ATGGGGTGGATAGCAATGAC; R: CCATAACACGCTCTTTATATGAACC
    PbPLR7 F: AAGCAAGCAGTCAAACAACC; R: AAAGAAGCCTACGGCATGG
    PbPLR8 F: GTCGCCGGAAACCACTTG; R: AGTTATTTGGTAAACCAACCTTGA
    PbEF1α F: CATTCAAGTATGCGTGGGT; R: ACGGTGACCAGGAGCA
    下载: 导出CSV

    表  2  闽楠PLR基因家族的基本信息

    Table  2.   Basic information of PLR gene family in P. bournei

    基因名称平均分子量/kDa等电点起始位置染色体位置氨基酸数量/个
    PbPLR135 441.985.9082085409~82088341chr01336
    PbPLR235 320.246.2382091009~82098484chr01312
    PbPLR335 378.056.32111385794~111390362chr02312
    PbPLR437 844.407.7480076163~80082267chr05312
    PbPLR535 378.056.3281056124~81058143chr05284
    PbPLR634 970.246.4152893843~52901318chr10316
    PbPLR734 910.236.1052876836~52885817chr10316
    PbPLR813 320.638.7652801277~52802274chr10123
    下载: 导出CSV
  • [1] 金重为, 邰瓞生, 尤纪雪. 天然耐腐木材的抗腐力及其在腐朽过程中化学成分的变化[J]. 林业科学, 1989, 25(5): 447 − 452. doi:  10.3321/j.issn:1001-7488.1989.05.014

    JIN Zhongwei, TAI Diesheng, YOU Jixue. Decay resistance of natural durable species and the changes in amounts of major components during decay of wood by Coriolus versicolor and Gloeophyllum trabeum [J]. Sci Silv Sin, 1989, 25(5): 447 − 452. doi:  10.3321/j.issn:1001-7488.1989.05.014
    [2] SCHULTZ T P, NICHOLAS D D, HENRY W P, et al. Review of laboratory and outdoor exposure efficacy results of organic biocide: antioxidant combinations, an initial economic analysis and discussion of a proposed mechanism [J]. Wood Sci Technol, 2005, 37: 175 − 184.
    [3] 周妮, 齐锦, 王燕高, 等. 桢楠现代木和阴沉木精油化学成分的GC-MS分析[J]. 西北农林科技大学学报(自然科学版), 2015, 43(6): 136 − 140.

    ZHOU Ni, QI Jin, WANG Yangao, et al. GC-MS analysis of chemical components of essential oils from recent and ancient buried Phoebe zhennan woods [J]. J Northwest A&F Univ Nat Sci Ed, 2015, 43(6): 136 − 140.
    [4] 吴大荣, 朱政德. 福建省罗卜岩自然保护区闽楠种群结构和空间分布格局初步研究[J]. 林业科学, 2003, 39(1): 23 − 30. doi:  10.3321/j.issn:1001-7488.2003.01.004

    WU Darong, ZHU Zhengde. Preliminary study on structure and spatial distribution pattern of Phoebe bournei in Luoboyan Nature Reserve in Fujian Province [J]. Sci Silv Sin, 2003, 39(1): 23 − 30. doi:  10.3321/j.issn:1001-7488.2003.01.004
    [5] STIRLING R, MORRIS P I. Potential contributions of lignans to decay resistance in western red cedar [J]. Wood Sci Technol, 2016, 50: 399 − 412. doi:  10.1007/s00226-015-0784-y
    [6] 贾贤, 黄秋生, 刘光华, 等. 我国楠木资源的研究现状[J]. 中国园艺文摘, 2014, 30(10): 55 − 59. doi:  10.3969/j.issn.1672-0873.2014.10.023

    JIA Xian, HUANG Qiusheng, LIU Guanghua, et al. Research status of Phoebe resources in China [J]. Chin Hortic Abstr, 2014, 30(10): 55 − 59. doi:  10.3969/j.issn.1672-0873.2014.10.023
    [7] NAKATSUBO T, MIZUTANI M, SUZUKI S, et al. Characterization of Arabidopsis thaliana pinoresinol reductase, a new type of enzyme involved in lignan biosynthesis [J]. J Biol Chem, 2008, 283(23): 15550 − 15557. doi:  10.1074/jbc.M801131200
    [8] DINKOVA-KOSTOVA A T, GANG D R, DAVIN L B, et al. (+)-Pinoresinol/(+)-lariciresinol reductase from Forsythia intermedia [J]. J Biol Chem, 1996, 271(46): 29473 − 29482. doi:  10.1074/jbc.271.46.29473
    [9] CHU A, DINKOVA A, DAVIN L B, et al. Stereospecificity of (+)-pinoresinol and (+)-lariciresinol reductases from Forsythia intermedia [J]. J Biol Chem, 1993, 268(36): 27026 − 27033. doi:  10.1016/S0021-9258(19)74213-6
    [10] XIAO Ying, JI Qian, GAO Shouhong, et al. Combined transcriptome and metabolite profiling reveals that IiPLR1 plays an important role in lariciresinol accumulation in Isatis indigotica [J]. J Exp Bot, 2015, 66(20): 6259 − 6271. doi:  10.1093/jxb/erv333
    [11] von HEIMENDAHL C B I, SCHAFER K M, EKLUND P, et al. Pinoresinol-lariciresinol reductases with different stereo specificity from Linum album and Linum usitatissimum [J]. Phytochemistry, 2005, 66: 1254 − 1263. doi:  10.1016/j.phytochem.2005.04.026
    [12] RENOUARD S, CORBIN C, LOPEZ T, et al. Abscisic acid regulates pinoresinol-lariciresinol reductase gene expression and secoisolariciresinol accumulation in developing fax (Linum usitatissimum L. ) seeds [J]. Planta, 2012, 235: 85 − 98. doi:  10.1007/s00425-011-1492-y
    [13] CORBIN C, RENOUARD S, LOPEZ T, et al. Identifcation and characterization of cis-acting elements involved in the regulation of ABA- and/or GA-mediated LuPLR1 gene expression and lignan biosynthesis in fax (Linum usitatissimum L. ) cell cultures [J]. J Plant Physiol, 2013, 170(5): 516 − 522. doi:  10.1016/j.jplph.2012.11.003
    [14] MIN T, KASAHARA H, BEDGAR D L, et al. Crystal structures of pinoresinol-lariciresinol and phenylcoumaran benzylic ether reductases and their relationship to isofavone reductases [J]. J Biol Chem, 2003, 278(50): 50714 − 50723. doi:  10.1074/jbc.M308493200
    [15] SUN Yanni, BUHLER J. Designing patterns and profiles for faster HMM search [J]. IEEE/ACM Trans Comput Biol Bioinform, 2009, 6(2): 232 − 243. doi:  10.1109/TCBB.2008.14
    [16] LETUNIC I, BORK P. 20 years of the SMART protein domain annotation resource [J]. Nucleic Acids Res, 2018, 46: 493 − 496. doi:  10.1093/nar/gkx1177
    [17] HEMMATI S, SCHMIDT T J, FUSS E, et al. (+)-Pinoresinol/(−)-lariciresinol reductase from Linum perenne Himmelszelt involved in the biosynthesis of justicidin B [J]. FEBS Lett, 2007, 581: 603 − 610. doi:  10.1016/j.febslet.2007.01.018
    [18] RICKLEFS E, GIRHARD M, URLACHER V B. Three-steps in one-pot: whole-cell biocatalytic synthesis of enantiopure (+)- and (−)-pinoresinol via kinetic resolution[J/OL]. Microb Cell Fact. 2016, 15: 78[2022-01-20]. doi:  10.1186/s12934-016-0472-0.
    [19] HANO C, MARTIN I, FLINIAUX O, et al. Pinoresinol-lariciresinol reductase gene expression and secoisolariciresinol diglucoside accumulation in developing fax (Linum usitatissimum) seeds [J]. Planta, 2006, 224: 1291 − 1301. doi:  10.1007/s00425-006-0308-y
    [20] HEMMATI S, von HEIMENDAHL C B I, KLAES M, et al. Pinoresinol-lariciresinol reductases with opposite enantiospecifcity determine the enantiomeric composition of lignans in the diferent organs of Linum usitatissimum L. [J]. Planta Med, 2010, 76: 928 − 934. doi:  10.1055/s-0030-1250036
    [21] CHIANG N T, MA Liting, LEE Y R, et al. The gene expression and enzymatic activity of pinoresinol-lariciresinol reductase during wood formation in Taiwania cryptomerioides Hayata [J]. Holzforschung, 2019, 73(2): 197 − 208. doi:  10.1515/hf-2018-0026
    [22] MOUMMOU H, KALLBERG Y, TONFACK L B, et al. The plant short-chain dehydrogenase (SDR) superfamily: genome-wide inventory and diversification patterns[J/OL]. BMC Plant Biol, 2012, 12: 219[2022-02-25]. doi:  10.1186/1471-2229-12-219.
    [23] MARKULIN L, CORBIN C, RENOUARD S, et al. Pinoresinol-lariciresinol reductases, key to the lignan synthesis in plants [J]. Planta, 2019, 249(6): 1695 − 1714. doi:  10.1007/s00425-019-03137-y
    [24] FUJITA M, GANG D R, DAVIN L B, et al. Recombinant pinoresinol-lariciresinol reductases from western red cedar (Thuja plicata) catalyze opposite enantiospecific conversions [J]. J Biol Chem, 1999, 274: 618 − 627. doi:  10.1074/jbc.274.2.618
    [25] CORBIN C, DROUET S, MATELJAK I, et al. Functional characterization of the pinoresinol-lariciresinol reductase-2 gene reveals its roles in yatein biosynthesis and fax defense response [J]. Planta, 2017, 243: 405 − 420.
  • [1] 尚林雪, 王群, 张国哲, 赵雨, 顾翠花.  紫薇LiCMB1基因的克隆及表达特性分析 . 浙江农林大学学报, doi: 10.11833/j.issn.2095-0756.20220333
    [2] 卓娟, 侯丹, 林新春.  毛竹PhebHLH6基因克隆及表达分析 . 浙江农林大学学报, doi: 10.11833/j.issn.2095-0756.20220553
    [3] 周佩娜, 党静洁, 邵永芳, 石遵睿, 张琳, 刘潺潺, 吴啟南.  荆芥HD-Zip基因家族的全基因组鉴定及分析 . 浙江农林大学学报, 2023, 40(1): 12-21. doi: 10.11833/j.issn.2095-0756.20220390
    [4] 王绍良, 张雯宇, 高志民, 周明兵, 杨克彬, 宋新章.  毛竹磷转运蛋白Ⅰ家族基因鉴定及表达模式 . 浙江农林大学学报, 2022, 39(3): 486-494. doi: 10.11833/j.issn.2095-0756.20210471
    [5] 孟超敏, 耿翡翡, 卿桂霞, 周佳敏, 张富厚, 刘逢举.  陆地棉磷高效基因GhMGD3的克隆与表达分析 . 浙江农林大学学报, 2022, 39(6): 1203-1211. doi: 10.11833/j.issn.2095-0756.20220145
    [6] 娄永峰, 高志民.  毛竹早期光诱导蛋白基因克隆及功能分析 . 浙江农林大学学报, 2021, 38(1): 93-102. doi: 10.11833/j.issn.2095-0756.20200237
    [7] 王楠楠, 董彬, 杨丽媛, 赵宏波.  梅花2个PmWRKY2基因克隆及在逆境胁迫下的表达模式 . 浙江农林大学学报, 2021, 38(4): 812-819. doi: 10.11833/j.issn.2095-0756.20200706
    [8] 庞天虹, 钱婕妤, 付建新, 顾翠花, 张超.  桂花己糖激酶基因家族成员的序列及表达分析 . 浙江农林大学学报, 2021, 38(2): 225-234. doi: 10.11833/j.issn.2095-0756.20200370
    [9] 黄元城, 郭文磊, 王正加.  薄壳山核桃全基因组LBD基因家族的生物信息学分析 . 浙江农林大学学报, 2021, 38(3): 464-475. doi: 10.11833/j.issn.2095-0756.20200454
    [10] 阮诗雨, 张智俊, 陈家璐, 马瑞芳, 朱丰晓, 刘笑雨.  毛竹GRF基因家族全基因组鉴定与表达分析 . 浙江农林大学学报, 2021, 38(4): 792-801. doi: 10.11833/j.issn.2095-0756.20200544
    [11] 肖业玲, 童再康, 吴建辉.  闽楠微扦插繁育体系的建立 . 浙江农林大学学报, 2020, 37(1): 188-194. doi: 10.11833/j.issn.2095-0756.2020.01.025
    [12] 蒋琦妮, 付建新, 张超, 董彬, 赵宏波.  桂花OfAP1基因的克隆及表达分析 . 浙江农林大学学报, 2019, 36(4): 664-669. doi: 10.11833/j.issn.2095-0756.2019.04.005
    [13] 李冰冰, 刘国峰, 魏书, 黄龙全, 张剑韵.  烟草NtPLR1基因克隆与表达分析 . 浙江农林大学学报, 2017, 34(4): 581-588. doi: 10.11833/j.issn.2095-0756.2017.04.003
    [14] 姜顺邦, 张怡, 韦小丽, 范辉华, 徐小琴, 段如雁.  不同闽楠优树子代苗期生长及光合生理特性比较 . 浙江农林大学学报, 2016, 33(1): 51-59. doi: 10.11833/j.issn.2095-0756.2016.01.007
    [15] 赵婷, 韩小娇, 刘明英, 乔桂荣, 蒋晶, 姜彦成, 卓仁英.  东南景天耐镉相关基因SaFer的克隆与功能初步分析 . 浙江农林大学学报, 2015, 32(1): 25-32. doi: 10.11833/j.issn.2095-0756.2015.01.004
    [16] 赵传慧, 周厚君, 童再康, 林二培, 黄华宏, 牛明月.  光皮桦成花相关MADS-box基因BlMADS1的克隆与表达 . 浙江农林大学学报, 2015, 32(2): 221-228. doi: 10.11833/j.issn.2095-0756.2015.02.008
    [17] 庞景, 童再康, 黄华宏, 林二培, 刘琼瑶.  杉木纤维素合成酶基因CesA的克隆及表达分析 . 浙江农林大学学报, 2015, 32(1): 40-46. doi: 10.11833/j.issn.2095-0756.2015.01.006
    [18] 杨希宏, 黄有军, 陈芳芳, 黄坚钦.  山核桃FLOWERING LOCUS C同源基因鉴定与表达分析 . 浙江农林大学学报, 2013, 30(1): 1-8. doi: 10.11833/j.issn.2095-0756.2013.01.001
    [19] 范繁荣.  闽楠群落种间关联性的灰色系统理论分析 . 浙江农林大学学报, 2008, 25(1): 33-36.
    [20] 李铁华, 朱祥云.  闽楠种子活力测定方法的研究 . 浙江农林大学学报, 2003, 20(3): 321-324.
  • 加载中
  • 链接本文:

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

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2022/6/1173

图(6) / 表(2)
计量
  • 文章访问数:  160
  • HTML全文浏览量:  24
  • PDF下载量:  85
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-05-16
  • 录用日期:  2022-08-08
  • 修回日期:  2022-07-15
  • 网络出版日期:  2022-11-21
  • 刊出日期:  2022-11-21

闽楠PLR基因家族鉴定及响应激素的表达分析

doi: 10.11833/j.issn.2095-0756.20220351
    基金项目:  浙江省农业新品种选育重大专项(2021C2070-10);国家自然科学基金资助项目(32171828)
    作者简介:

    王倩清(ORCID: 0000-0003-0116-0551),从事闽楠遗传学研究。E-mail: 974317348@qq.com

    通信作者: 童再康(ORCID: 0000-0002-0202-6496),教授,博士,从事林木新品种选育与定向培育研究工作。E-mail: zktong@zafu.edu.cn
  • 中图分类号: S722.3

摘要:   目的  解析闽楠Phoebe bournei木脂素合成途径关键酶基因Pinoresinol-lariciresinol reductase (PLR)特征及不同激素处理下的表达模式,探索闽楠PLR基因家族的生物学功能。  方法  利用生物信息学方法鉴定闽楠基因组中的PLR基因家族成员,分析基因结构、保守基序、染色体定位、系统进化、启动子顺式作用元件和激素响应。  结果  从闽楠全基因组内共鉴定出8个闽楠PLR (PbPLR)成员,不均匀地分布于第1、2、5和10号染色体。系统进化分析表明:PbPLR功能可能与底物立体选择性有关。启动子元件分析发现:PbPLRs启动子区域存在脱落酸、水杨酸和茉莉酸甲酯等激素响应元件;基因表达分析表明:PbPLRs基因的表达具有组织特异性和激素响应特征,其中PbPLRs基因在根中表达量最高,其次茎,叶中表达量较低;3种激素处理时,PbPLR2表达诱导最强。  结论  从闽楠基因组中鉴定出的8个PbPLRs基因,其基因特征和不同成员可能具有不同催化酶立体选择性,成员表达模式多样化,可能受多种激素诱导,具有组织特异性。图6表2参25

English Abstract

王倩清, 张毓婷, 张俊红, 刘慧, 童再康. 闽楠PLR基因家族鉴定及响应激素的表达分析[J]. 浙江农林大学学报, 2022, 39(6): 1173-1182. doi: 10.11833/j.issn.2095-0756.20220351
引用本文: 王倩清, 张毓婷, 张俊红, 刘慧, 童再康. 闽楠PLR基因家族鉴定及响应激素的表达分析[J]. 浙江农林大学学报, 2022, 39(6): 1173-1182. doi: 10.11833/j.issn.2095-0756.20220351
WANG Qianqing, ZHANG Yuting, ZHANG Junhong, LIU Hui, TONG Zaikang. Identification of PLR gene family and expression of responsive hormones in Phoebe bournei[J]. Journal of Zhejiang A&F University, 2022, 39(6): 1173-1182. doi: 10.11833/j.issn.2095-0756.20220351
Citation: WANG Qianqing, ZHANG Yuting, ZHANG Junhong, LIU Hui, TONG Zaikang. Identification of PLR gene family and expression of responsive hormones in Phoebe bournei[J]. Journal of Zhejiang A&F University, 2022, 39(6): 1173-1182. doi: 10.11833/j.issn.2095-0756.20220351
  • 不同树种木材的天然耐腐性差异大,据统计每年全球因木材腐朽造成的损失占采伐利用木材的10%以上,严重影响木材的加工利用[1]。木材腐朽主要由腐朽菌和变色菌等真菌侵蚀木材细胞所引起,而木材的天然耐腐性则依靠自身代谢物的生物杀菌、抑菌活性来实现[2]。闽楠Phoebe bournei为樟科Lauraceae楠属Phoebe植物,是中国特有的珍稀用材与绿化树种。楠木具有极强耐腐性,阴沉木历经千年不腐,被誉为“木中之玉”,市场俗称“金丝楠乌木”[3]。闽楠和桢楠Phoebe zhennan作为“金丝楠木”最主要的产材树种,树干挺拔,木材颜色呈淡黄色,坚硬细密,纹理美观具有芳香香气,属上等家具良材[4]。因此,闽楠可作为研究木材天然耐腐性的良好材料,闽楠全基因组测序的完成为基因挖掘和功能研究提供了重要基础。

    已有研究发现:从木材分离的木脂素(lignan)类化合物具有抗菌活性,参与植物防御系统,增强木材耐久性和抗真菌能力。例如北美乔柏Thuja plicata心材中提取的木脂素通过清除自由基和螯合亚铁来干扰真菌氧化还原循环提高木材耐腐性[5],南洋杉Araucaria araucana中提取的松脂素(pinoresinol)和开环落叶松脂素(secoisolariciresinol)对木材腐朽菌具有较好抗性[6]。然而,关于木脂素生物合成调控途径尚不明确,从而限制其深入研究与生产应用。

    松脂素还原酶Pinoresinol-lariciresinol reductase (PLR)是一种编码松脂醇的落叶松脂素还原酶基因,属木脂素生物合成的关键酶,一分子松脂素在PLR松脂素-落叶松脂素还原酶作用下经两步还原成落叶松脂素和异落叶松脂素,也是木脂素生物合成途径中重要限速步骤之一。目前已从模式植物拟南芥Arabidopsis thaliana[7]和富含木脂素的植物(连翘Forsythia suspense[8-9]、菘蓝Isatis indigotica[10]、亚麻属Linum[11-13])中克隆得到PLR基因。迄今为止,所有以重组纯化蛋白为特征的PLR均可以将松脂素还原为落叶松树脂素以及将落叶松树脂素还原为异落叶松树脂素,但AtPrR1对落叶松树脂素的活性比松脂素低约35倍,AtPrR2对作为底物的落叶松树脂素完全没有催化活性。因此,MIN等[14]将它们命名为松脂醇还原酶,而不是松脂醇-落叶松树脂醇还原酶。PLR的表达受激素的调控,不同物种在不同激素处理下松脂醇-落叶松树脂醇还原酶的表达模式有所不同。在脱落酸(ABA)处理下,LuPLR1的表达增强,LuPLR1的表达似乎是由在种子成熟和休眠中起关键作用的植物激素ABA触发的[12-13]。所以赤霉酸(GA)对LuPLR1的表达具有相反的调节作用[13]。在茉莉酸甲酯(MeJA)的处理下,LuPLR2的表达量有明显提高,因为亚麻在MeJA处理和伤害后,LuPLR2基因表达和叶酸蛋白产量增加,这与这些木脂素参与植物防御反应一致。本研究选择耐腐性强的闽楠解析木脂素合成重要限速酶PLR基因,利用闽楠基因组数据,对PLR基因家族成员进行鉴定和基因结构、保守基序及染色体定位分析,构建系统发育进化树预测闽楠PLR (PbPLR)催化酶的立体选择性,为进一步研究PbPLR基因在闽楠木脂素合成过程中的生物学功能提供基础。

    • 采用BLAST和hmmsearch鉴定闽楠PLR家族成员。首先在美国国家生物技术信息中心(NCBI)数据库(https://www.ncbi.nlm.nih.gov/)下载拟南芥PLR家族蛋白序列,通过BLASTP同源比对检索,用HMMER 3.1软件[15]搜索闽楠PLR基因家族的候选序列。搜索获得的基因成员在人工筛选后提交至SMART数据库[16]确认保守结构域,确定PLR基因成员,并依据染色体位置进行命名。

    • 采用MapChart软件绘制PbPLR基因在闽楠基因组上的位置信息,并利用MCScanX软件分析PbPLR基因的共线性区域和复制情况。

    • 根据PbPLR基因编码区序列和基因组序列,采用在线GSDS (http://gsds.gao-lab.org)分析PbPLR基因结构。并通过在线工具MEME Suite (http://meme-suite.org/)分析PbPLR氨基酸序列中保守基序,设置保守基序个数为20。最后利用SnapGene软件进行PbPLR蛋白序列比对。

    • 利用植物蛋白数据库下载的拟南芥、亚麻、台湾杉Taiwania cryptomerioides等物种PLR蛋白序列,与PbPLR家族成员进行系统进化分析。利用MEGA X软件构建进化树,采用邻接法(neighbor-joining method,NJ),重复1 000次(Bootstrap为1 000),利用在线工具ITOL (https://itol.embl.de/)美化MEGA X软件输出的进化树。

    • 根据闽楠基因组数据,利用TBtools软件提取PbPLR基因序列转录起始位点上游2 000 bp序列,使用在线软件PlantCARE (http://intra.psb.ugent.be:8080/PlantCARE/)预测PbPLR基因启动子序列中的顺式调控元件。

    • 选取1.5年生长势一致的闽楠半同胞家系苗进行激素处理,分别采集叶、茎和根进行PbPLRs基因表达分析。基于前期预实验结果,本研究分别使用1 mmol·L−1 ABA、1 mmol·L−1水杨酸(SA)和2 mmol·L−1 MeJA处理闽楠实生苗,对土壤上部植株均匀喷洒激素,使土壤表面被激素浸湿,同时,以喷洒水为对照(ck)。激素喷施时间为2021年10月1日9:30、10月2日11:00、10月3日11:30、10月3日22:00、10月4日9:00、10月4日9:30。在处理后1、3、12、24、48和72 h取样并液氮速冻。采用CTAB法提取RNA,使用Takara公司的PrimeScriptTMRT reagent Kit (Perfect Real Time)试剂盒反转录合成cDNA。利用在线网站(http://primer3.ut.ee/)设计PbPLRs特异性引物(表1)用于实时荧光定量PCR实验。定量PCR反应体系(10 μL)为:SYBR Green 5.0 μL,上下游引物(10 μmol·L−1)各0.2 μL,模板0.8 μL,双蒸水3.8 μL。反应程序为:95 ℃预变性3 min,95 ℃ 5 s,60 ℃ 30 s (此步完成时采集荧光信号),共40个循环。利用CFX96™ Real-Time PCR检测系统(Bio-Rad,美国)。内参基因为PbEF1α,按照$2^{-\Delta \Delta {{C}}_{t}}$法计算目的基因相对表达量。运用SPSS软件对定量数据进行分析。

      表 1  实时荧光定量PCR (RT-qPCR)引物信息

      Table 1.  RT-qPCR primer information

      引物名称引物序列(5′→3′)
      PbPLR1 F: GGCACAAATCTCATCAGCAA; R: GCCAGCTTTTATGGCCTGTA
      PbPLR2 F: TCTGGATTACGTGCAACAGG; R: AGCGTTTGAGATAGTCCTTCG
      PbPLR3 F: AGGTCTGGATTACGCACAGC; R: CGTTTGAGATAGTCCTTCACTCG
      PbPLR4 F: GGAATTCCTTGCTCAAATGAAA; R: TTTCGAAATTTGCCAGACAAC
      PbPLR5 F: GTGGAACAAGAGCCCAAATG; R: CTGGCAAGTGGGTGAAAAAT
      PbPLR6 F: ATGGGGTGGATAGCAATGAC; R: CCATAACACGCTCTTTATATGAACC
      PbPLR7 F: AAGCAAGCAGTCAAACAACC; R: AAAGAAGCCTACGGCATGG
      PbPLR8 F: GTCGCCGGAAACCACTTG; R: AGTTATTTGGTAAACCAACCTTGA
      PbEF1α F: CATTCAAGTATGCGTGGGT; R: ACGGTGACCAGGAGCA
    • 基于拟南芥AtPrRs蛋白序列,应用BLAST和HMMER 3.1软件在闽楠基因组中搜索PLR家族成员。利用BLAST在闽楠基因组中鉴定出12个候选PLR基因,通过HMMER 3.1鉴定出10个成员,去除重复和结构域鉴定后,共鉴定出8个PLR基因家族成员(表2),按照其在染色体上位置命名为PbPLR1~PbPLR8。

      表 2  闽楠PLR基因家族的基本信息

      Table 2.  Basic information of PLR gene family in P. bournei

      基因名称平均分子量/kDa等电点起始位置染色体位置氨基酸数量/个
      PbPLR135 441.985.9082085409~82088341chr01336
      PbPLR235 320.246.2382091009~82098484chr01312
      PbPLR335 378.056.32111385794~111390362chr02312
      PbPLR437 844.407.7480076163~80082267chr05312
      PbPLR535 378.056.3281056124~81058143chr05284
      PbPLR634 970.246.4152893843~52901318chr10316
      PbPLR734 910.236.1052876836~52885817chr10316
      PbPLR813 320.638.7652801277~52802274chr10123
    • PbPLRs基因不均匀地分布在闽楠第1、2、5和10号染色体,其中10号染色体上数目最多,含3个PbPLRs成员,2号染色体仅1个成员。闽楠基因复制事件检测结果表明:PbPLR1和PbPLR3经历了片段复制,而PbPLR1和PbPLR2以及PbPLR6、PbPLR7和PbPLR8经历了串联复制。

    • 利用GSDS软件绘制PbPLRs基因结构图(图1A)可见:8个PbPLRs均含有至少2个外显子,其中PbPLR2、PbPLR3、PbPLR4和PbPLR7各含4个外显子,每个成员均有内含子。

      图  1  PbPLR基因结构及蛋白保守基序分析

      Figure 1.  Analysis of PbPLR gene structure and protein conserved motifs

      对PbPLR蛋白序列进行保守基序(Motif)分析,共检测到14个保守基序(图1B)。多序列比对分析发现:PbPLRs都具有保守的NADPH结合位点和赖氨酸(Lys)活性位点(图1C),PbPLR1含有特异Motif 11、Motif 13和Motif 14;PbPLR 6和PbPLR 8含有特有的Motif 10。

    • 为研究PbPLRs家族的进化关系并预测其催化酶功能,选取裸子植物台湾杉、北美乔柏,以及被子植物拟南芥、连翘和亚麻等物种的PLR蛋白序列构建系统进化树(图2)。系统进化分析表明:裸子植物的台湾杉TcPLRs和北美乔柏TpPLRs属同一分支,即裸子植物和被子植物PLRs不论其底物立体选择性如何,均可分为2分支群[17-20]。然而,被子植物PLRs根据其底物立体选择性而分属不同支群,根据位于同一分支的LcPLR1、LaPLR1和LuPLR2的底物偏好,推测PbPLR1、PbPLR6、PbPLR7和PbPLR8偏好(+)-松脂醇[(+)-pinoresinol]和(+)-落叶松脂醇[(+)-lariciresinol];而PbPLR2、PbPLR3和PbPLR4同属于一个小分支,与偏好(−)-松脂醇和(−)-落叶松脂醇的LuPLR1和偏好(+)-松脂醇和(−)-落叶松脂醇的LpPLR1亲缘关系较近,推测PbPLR2、PbPLR3和PbPLR4可将(−)-松脂醇还原为(−)-落叶松脂醇。PbPLR5与其他PLR成员亲缘关系较远,不能预测其选择性。

      图  2  闽楠和多个物种的PLR家族系统发育树

      Figure 2.  PLR family phylogenetic tree of P. bournei and several species

    • 启动子区域的顺式调控元件可能影响基因表达模式。分别提取PbPLR基因上游2 000 bp序列进行顺式作用元件分析,发现启动子区域具多种类型响应元件(图3A),主要包括激素响应元件如脱落酸(ABA)响应元件、赤霉素(GA)响应元件、茉莉酸甲酯(MeJA)响应元件和水杨酸(SA)响应元件,与环境胁迫响应元件等,对PbPLR基因启动子区的激素响应元件统计分析表明:PbPLRs含有5种激素响应元件:脱落酸、水杨酸、赤霉素、茉莉酸甲酯和生长素(Auxin),每个PbPLRs至少含有2个激素响应元件,其中PbPLR2激素响应元件最多(6个),PbPLR4仅2个(图3B)。在全部PbPLR启动子区域的响应元件中,激素响应元件最多的是脱落酸(9个),其次为赤霉素(7个),生长素最少仅2个(图3C)。

      图  3  PbPLR家族基因的顺式作用元件分析

      Figure 3.  Analysis of putative cis-elements in the promoter regions of PbPLRs

    • 应用实时荧光定量PCR (RT-qPCR)反应检测了ABA、SA和MeJA处理后根、茎和叶中的PbPLRs基因的表达水平(图4~6)。结果表明:8个PbPLRs基因对ABA、SA和MeJA均有响应,处理1~48 h,基因表达水平显著上调,48 h时达最高值,之后上调幅度降低。其中,ABA处理48 h后根部PbPLR2表达量上调倍数最大,为ck的400倍,茎部表达水平上调为ck的300倍,但叶片中表达上调不明显。相较于在根部和叶部,ABA处理后的PbPLR6在茎部表达量上调倍数最大,48 h时的表达量上调倍数为ck的35倍;ABA处理48 h后的PbPLR4在叶中上调倍数最大,为ck的30倍(图4)。

      图  4  PbPLRs基因在ABA处理下的表达模式

      Figure 4.  Expression patterns of PbPLRs under ABA treatment

      SA处理后,根中PbPLR2表达量上调倍数最大,为ck的400倍,PbPLR6在根、茎、叶中均上调表达,而PbPLR1、PbPLR4和PbPLR8上调倍数较低(图5)。

      图  5  PbPLRs基因在SA处理下的表达模式

      Figure 5.  Expression patterns of PbPLRs under SA treatment

      MeJA处理48 h后,根中PbPLR2表达量上调倍数为ck的500倍,且茎中也达400倍;PbPLR3在根中表达上调为ck的140倍(图6)。PbPLR5对MeJA处理响应不明显,茎和根中表达量上调幅度均较低,但ABA(图4)和SA(图5)处理后PbPLR5的表达量上调显著(P<0.05),该结果与PbPLR5启动子区域含有的激素响应元件一致。

      图  6  PbPLRs基因在MeJA处理下的表达模式

      Figure 6.  Expression patterns of PbPLRs under MeJA treatment

    • 木脂素类化合物具有较为广泛的应用前景,PLR作为植物木脂素生物合成的重要限速酶基因之一,引起了广泛的关注。例如,在亚麻中鉴定到2个PLR基因(LuPLR1和LuPLR2)[12],模式植物拟南芥中共有2个PrR基因(AtPrR1和AtPrR2)[7],台湾杉中发现了3个PLR基因(TcPLR1、TcPLR2和TcPLR3)等[21]。本研究从闽楠基因组中共鉴定到8个PbPLRs基因,不均等分布于闽楠的4条染色体上,且均含有PLR家族保守的NADPH结合位点和赖氨酸活化位点。PLR属于short-chain dehydrogenase/reductase (SDR)超级家族,该家族广泛参与植物初级代谢和次级代谢过程,且在高等植物中具有明显功能分化[22]。闽楠中数量较多的PLR成员可能参与不同代谢通路,并发生功能分化。

      利用不同物种PLR蛋白构建系统进化树,发现裸子植物PLR蛋白单独聚为一类,而被子植物PLR系统聚类与催化酶的立体选择性相关,该结果与前人研究结果相一致。NAKATSUBO等[7]发现被子植物PLR (FiPLR1、LaPLR1、LuPLR1、LpPLR1、AtPrR1和AtPrR2)和裸子植物PLR (TpPLR1和TpPLR2)分别聚类。PLR催化酶的底物立体选择性直接决定了下游木脂素的结构骨架,是形成木脂素类化合物结构多样性的关键节点。PLR立体选择性最早是在金钟连翘Forsythia intermedia中发现,功能研究发现2个FiPLR都偏好选择催化(+)-松脂醇和(+)-落叶松脂醇为(−)-异落叶松脂醇[9]。随后,在其他物种中鉴定的PLRs同样具有明显立体选择性[23]。拟南芥AtPrR1可催化2种松脂醇异构体,而AtPrR2仅能催化(−)-松脂醇形成(−)-落叶松脂醇[7]。根据已有催化酶的功能研究,推测闽楠PbPLR1、PbPLR6、PbPLR7和PbPLR8偏好选择(+)-松脂醇和(+)-落叶松脂醇,而PbPLR2、PbPLR3、PbPLR4和PbPLR5聚类的PLR蛋白的立体选择性不同,无法通过系统进化树推测其底物偏好性。FUJITA等[24]研究发现:北美乔柏的TpPLRs聚类为一个分支却具有相反的立体选择性。PLR家族成员的蛋白序列差异与底物立体选择性并不完全等同,对底物的选择偏好性可能随物种共同经历多次进化;或形成于物种进化的早期,随后的碱基突变并不影响其催化功能[25]

      通过启动子顺式作用元件分析发现:PbPLR基因的启动子区域具有多种响应元件,推测PbPLR表达模式可能受生物及非生物因素调控。例如光响应、厌氧诱导、耐旱和防御等响应元件,推测PbPLR表达可能受光信号及非生物胁迫等因素影响,并通过次生代谢物的催化合成而参与植物防御系统。此外,PbPLRs基因启动子区域具有多种激素响应调控元件,其中ABA、MeJA和SA等3种响应元件数目最多。不同物种PLR基因表达受不同激素调控。例如,ABA处理下亚麻LuPLR1和菘蓝LiPLR1表达显著上调[10, 12-13];MeJA诱导LuPLR2、IiPLR1和鬼臼草Podophyllum hexandrum PhPLR表达量上调[10, 25];而GA对LuPLR1表达起抑制作用[21]。对PbPLR表达分析表明:PbPLR对ABA、MeJA和SA均有响应,这与启动子元件分析相吻合。MeJA处理下PbPLR2响应最显著,在处理48 h根部表达量最高;ABA处理下PbPLR4响应最为显著,在处理48 h时叶片表达量最高;SA处理下PbPLR6响应最为显著,在处理48 h时茎表达量最高。本研究表明:PbPLR调控模式较为复杂,可能受不同植物激素诱导,同时PbPLRs组织特异性表达可能与其结构功能分化相关但与系统进化分支无关。推测闽楠木脂素合成具有较为复杂的调控网络,需进一步研究其生物合成调控机制。

参考文献 (25)

目录

    /

    返回文章
    返回