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

留言板

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

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

杨树MYC基因家族成员表达模式研究

胡梦璇 宋学勤 刘颖丽 赵树堂

引用本文:
Citation:

杨树MYC基因家族成员表达模式研究

    通讯作者: 赵树堂, shutang@caf.ac.cn
  • 中图分类号: S718.46

Analysis of Expression Patterns of MYC Family Members in Populus

    Corresponding author: ZHAO Shu-tang, shutang@caf.ac.cn
  • CLC number: S718.46

  • 摘要: 目的 探究bHLH亚家族成员MYC基因在杨树生长发育及对环境响应中的表达模式,为解析茉莉酸信号通路关键基因MYC调控杨树生长发育及抗逆响应提供参考。 方法 利用生物信息学方法鉴定杨树基因组中的MYC基因家族成员并对各成员的基因结构、保守基序进行了系统分析;在此基础上,采用实时荧光定量PCR技术检测了MYC成员在不同组织、不同植物激素响应及逆境胁迫下的表达模式。 结果 毛果杨基因组中含有10个MYC成员,这些成员在进化上相对保守,均具有典型的bHLH结构域,根据进化关系将其分为3个分支。组织表达分析发现,大部分MYC家族成员在根中均具有较高表达量,分枝Ⅱ两对基因在茎中呈相反表达模式。不同激素处理及非生物胁迫条件下MYC基因存在显著表达差异,但基因对之间多存在类似的表达模式。 结论 杨树MYC基因家族成员可能参与不同的生物学过程,研究结果为深入解析杨树MYC基因的功能奠定了基础。
  • 图 1  杨树和拟南芥MYC基因家族系统进化树

    Figure 1.  Phylogenetic tree of MYC family members between Populus and Arabidopsis

    图 2  杨树MYC家族成员基因结构(A)和蛋白保守结构域(B)分布

    Figure 2.  Gene structure (A) and distribution of conversed motifs (B) of PtrMYC family members in Populus

    图 3  杨树MYC基因组织表达分析

    Figure 3.  Expression analyses of MYCs in Populus

    图 4  生长素、细胞分裂素、茉莉酸处理后杨树MYC基因的表达模式

    Figure 4.  Expression patterns of PtrMYCs under NAA、6-BA、MeJA treatment

    图 5  胁迫条件下MYC基因的表达模式

    Figure 5.  Expression patterns of PtrMYCs under 4 ℃, NaCl and PEG6000 treatment

    表 1  实时荧光定量PCR引物序列

    Table 1.  The primers used for quantitative real-time PCR

    引物
    Primer
    序列(5′−3′)
    Primer sequence(5′−3′)
    引物
    Primer
    序列(5′−3′)
    Primer sequence(5′−3′)
    PtrMYC2a-F TATGCCAATGTGTCTGTGATGAATG PtrMYC2a-R TCTTGCGTGTAAAACCGACTCC
    PtrMYC2b-F AATGCTGGGATTCCGTCAACTAC PtrMYC2b-R TGATGATTGTTATTGTTCGCAGTCT
    PtrbHLH14a-F GTGGTCCCGAATGTGTCAAAGATG PtrbHLH14a-R ATGTTCGCATTACCAGAGATTTTTG
    PtrbHLH14b-F TACTCAAGCAAGGTTAGCAACAAGC PtrbHLH14b-R ATCAACCAATCTCTCCAGGTCCAT
    PtrbHLH14c-F GCTCGTGAAACAGAAGAAGCAAAT PtrbHLH14c-R TAGAGGAAAGTCCGAGTCTGAATGC
    PtrbHLH14d-F TGACACTATGGACAATCAAAGCACC PtrbHLH14d-R ATTCACATTATCGGATTGGACCCTA
    PtrbHLH3a-F CAACAAGCACGAATCAAGTTTATGG PtrbHLH3a-R GGACAAGTCATCATTTGGCTGC
    PtrbHLH3b-F GGAGTCTCACCCCGTTTCTACC PtrbHLH3b-R CAGTCTTGAGCAGTGATTTGATGTT
    PtrbHLH17a-F CAAGGGACACCAGGTGAGGTTTAT PtrbHLH17a-R CCTGCTTTTGACCTTGGTAGTGATT
    PtrbHLH17b-F AGGCTTCATTGTTGGGAGATGCTAT PtrbHLH17b-R AATGCTCCCCGACTTTTCCCTT
    PtrActin-F AAACTGTAATGGTCCTCCCTCCG PtrActin-R GCATCATCACAATCACTCTCCGA
    下载: 导出CSV

    表 2  杨树MYC基因家族成员基本信息

    Table 2.  Basic information of MYC family members in Populus

    基因名称
    Gene
    基因ID
    Gene number
    基因在染色体位置
    Location
    蛋白质长度
    Protein length/aa
    等电点
    PI
    分子量
    MW/kDa
    亚细胞定位
    Subcellular localization
    PtrMYC2a Potri.001G142200 Chr01( + )11649287-11651987 648 5.38 71.11 细胞核
    PtrMYC2b Potri.003G092200 Chr03( − )11839709-11842763 660 5.49 72.09 细胞核
    PtrBHLH14a Potri.001G083500 Chr01( − )6613720-6615721 472 7.17 52.73 细胞核/细胞质
    PtrBHLH14b Potri.003G147300 Chr03( + )16158005-16160025 466 7.72 52.34 细胞质/细胞核
    PtrBHLH14c Potri.002G176900 Chr02( + )13775139-13776876 492 5.40 54.70 细胞核
    PtrBHLH14d Potri.014G103700 Chr14( + )6950999-6952650 494 5.91 54.87 细胞核/叶绿体
    PtrBHLH3a Potri.002G042000 Chr02( − )2756159-2758458 504 6.59 54.92 细胞核/叶绿体
    PtrBHLH3b Potri.005G221100 Chr05( + )22329285-22332096 505 6.13 55.44 细胞核
    PtrBHLH17a Potri.002G172101 Chr02( − )13241064-13244032 619 6.14 68.25 细胞核
    PtrBHLH17b Potri.014G099700 Chr14( − )6539405-6542974 620 6.16 68.02 细胞核
    下载: 导出CSV
  • [1]

    BROWSE J. Jasmonate passes muster: A receptor and targets for the defense hormone[J]. Annual Review of Plant Biology, 2009, 60(1): 183-205. doi: 10.1146/annurev.arplant.043008.092007
    [2]

    ACOSTA I F, FARMER E E. Jasmonates[J]. The Arabidopsis Book, 2010, 8: e0129. doi: 10.1199/tab.0129
    [3]

    CHINI A, BOTER M, SOLANO R. Plant oxylipins: COI1/JAZs/MYC2 as the core jasmonic acid-signalling module[J]. The FEBS Journal, 2009, 276(17): 4682-4692. doi: 10.1111/j.1742-4658.2009.07194.x
    [4]

    SONG C, CAO Y P, DAI J, et al. The Multifaceted Roles of MYC2 in plants: toward transcriptional reprogramming and stress tolerance by jasmonate signaling[J]. Frontiers in Plant Science, 2022, 13: 868874. doi: 10.3389/fpls.2022.868874
    [5] 李 罡, 李文龙, 许雪梅, 等. MYC2转录因子参与植物发育调控的研究进展[J]. 植物生理学报, 2019, 55(2):125-132. doi: 10.13592/j.cnki.ppj.2018.0479

    [6]

    PIRES N, DOLAN L. Origin and diversification of basic-helix-loop-helix proteins in plants[J]. Molecular Biology Evolution, 2010, 27(4): 862-874. doi: 10.1093/molbev/msp288
    [7]

    LORENZO O, CHICO J M, SANCHEZ-SERRANO J J, et al. JASMONATE-INSENSITIVE1 encodes a MYC transcription factor essential to discriminate between different jasmonate-regulated defense responses in Arabidopsis[J]. The Plant Cell, 2004, 16(7): 1938-1950. doi: 10.1105/tpc.022319
    [8]

    DOMBRECHT B, XUE G P, SPRAGUE S J, et al. MYC2 differentially modulates diverse jasmonate-dependent functions in Arabidopsis[J]. The Plant Cell, 2007, 19(7): 2225-2245. doi: 10.1105/tpc.106.048017
    [9]

    CHEN Q, SUN J, ZHAI Q, et al. The basic helix-loop-helix transcription factor MYC2 directly represses PLETHORA expression during jasmonate-mediated modulation of the root stem cell niche in Arabidopsis[J]. The Plant Cell, 2011, 23(9): 3335-3352. doi: 10.1105/tpc.111.089870
    [10]

    ZHANG Q, XIE Z, ZHANG R, et al. Blue light regulates secondary cell wall thickening via MYC2/MYC4 activation of the NST1-directed transcriptional network in Arabidopsis[J]. The Plant Cell, 2018, 30(10): 2512-2528. doi: 10.1105/tpc.18.00315
    [11]

    LUO F, ZHANG Q, XIN H, et al. A Phytochrome B-PIF4-MYC2/MYC4 module inhibits secondary cell wall thickening in response to shaded light[J]. Plant Communications, 2022, 3(6): 100416. doi: 10.1016/j.xplc.2022.100416
    [12]

    WANG H P, LI Y, PAN J J, et al. The bHLH transcription factors MYC2, MYC3, and MYC4 are required for jasmonate-mediated inhibition of flowering in Arabidopsis[J]. Molecular Plant, 2017, 10(11): 1461-1464. doi: 10.1016/j.molp.2017.08.007
    [13]

    QI T, HUANG H, SONG S, et al. Regulation of jasmonate-mediated stamen development and seed production by a bHLH-MYB complex in Arabidopsis[J]. The Plant Cell, 2015, 27(6): 1620-1633. doi: 10.1105/tpc.15.00116
    [14]

    SONG S, HUANG H, GAO H, et al. Interaction between MYC2 and ETHYLENE INSENSITIVE3 modulates antagonism between jasmonate and ethylene signaling in Arabidopsis[J]. The Plant Cell, 2014, 26(1): 263-279. doi: 10.1105/tpc.113.120394
    [15]

    LI T, XU Y, ZHANG L, et al. The jasmonate activated transcription factor MdMYC2 regulates ETHYLENE RESPONSE FACTOR and ethylene biosynthetic genes to promote ethylene biosynthesis during apple fruit ripening[J]. The Plant Cell, 2017, 29(6): 1316-1334. doi: 10.1105/tpc.17.00349
    [16]

    DOB A, LAKEHAL A, NOVAK O, et al. Jasmonate inhibits adventitious root initiation through repression of CKX1 and activation of RAP2.6L transcription factor in Arabidopsis[J]. Journal of Experimental Botany, 2021, 72(20): 7107-7118. doi: 10.1093/jxb/erab358
    [17]

    KUMAR S, STECHER G, TAMURA K. MEGA7: Molecular evolutionary genetics analysis Version 7.0 for bigger datasets[J]. Molecular Biology Evolution, 2016, 33(7): 1870-1874. doi: 10.1093/molbev/msw054
    [18]

    JANSSON S, DOUGLAS C J. Populus: A Model System for Plant Biology[J]. Annual Review of Plant Biology, 2007, 58: 435-458. doi: 10.1146/annurev.arplant.58.032806.103956
    [19]

    SUN J, XU Y, YE S, et al. Arabidopsis ASA1 is important for jasmonate-mediated regulation of auxin biosynthesis and transport during lateral root formation[J]. The Plant Cell, 2009, 21(5): 1495-1511. doi: 10.1105/tpc.108.064303
    [20]

    HUANG C F, YU C P, WU Y H, et al. Elevated auxin biosynthesis and transport underlie high vein density in C4 leaves[J]. Proceedings National Academy of Sciences of USA, 2017, 114(33): E6884-E6891.
    [21]

    QI T, WANG J, HUANG H, et al. Regulation of jasmonate-induced leaf senescence by antagonism between bHLH subgroup IIIe and IIId factors in Arabidopsis[J]. The Plant Cell, 2015, 27(6): 1634-1649. doi: 10.1105/tpc.15.00110
    [22]

    ZHUO M N, SAKURABA Y, YANAGISAWA S. A jasmonate-activated MYC2–Dof2.1–MYC2 transcriptional loop promotes leaf senescence in Arabidopsis[J]. The Plant Cell, 2020, 32(1): 242-262. doi: 10.1105/tpc.19.00297
    [23]

    KAZAN K, MANNERS J M. The interplay between light and jasmonate signalling during defence and development[J]. Journal of Experimental Botany, 2011, 62(12): 4087-4100. doi: 10.1093/jxb/err142
    [24]

    LI Y, YANG X, LI X. Role of jasmonate signaling pathway in resistance to dehydration stress in Arabidopsis[J]. Acta Physiologiae Plantarum, 2019, 41(6): 1-12. doi: 10.1007/s11738-019-2897-7
    [25]

    AGRAWAL R, SHARMA M, DWIVEDI N, et al. MEDIATOR SUBUNIT17 integrates jasmonate and auxin signaling pathways to regulate thermomorphogenesis[J]. Plant Physiology, 2022, 189(4): 2259-2280. doi: 10.1093/plphys/kiac220
    [26]

    WANG Y, XU H, LIU W, et al. Methyl jasmonate enhances apple’ cold tolerance through the JAZ-MYC2 pathway[J]. The Plant Cell, 2019, 136(1): 75-84.
    [27]

    DING F, WANG C, XU N, et al. The ethylene response factor SlERF. B8 triggers jasmonate biosynthesis to promote cold tolerance in tomato[J]. Environmental and Experimental Botany, 2022, 203: 105073. doi: 10.1016/j.envexpbot.2022.105073
  • [1] 薛丽刘晓霞罗莹吕中睿张建国饶国栋 . 毛果杨MAP65基因家族的扩张与表达分析. 林业科学研究, 2021, 34(2): 92-101. doi: 10.13275/j.cnki.lykxyj.2021.02.010
    [2] 刘聪张洋夏德安陈雪冰魏志刚 . 毛果杨PLD基因家族全基因组水平鉴定及其盐胁迫下的表达分析. 林业科学研究, 2021, 34(3): 23-36. doi: 10.13275/j.cnki.lykxyj.2021.03.003
    [3] 单雪萌王思宁朱成磊高志民 . 毛竹 PeCPD 基因克隆与表达分析. 林业科学研究, 2019, 32(5): 58-66. doi: 10.13275/j.cnki.lykxyj.2019.05.008
    [4] 王菁李爱王春国宋文芹陈成彬 . 日本落叶松 UDPGDH基因的cDNA克隆和表达分析. 林业科学研究, 2013, 26(S1): 76-81.
    [5] 王江英范正琪殷恒福李辛雷吴斌李纪元 . 杜鹃红山茶CaAPX基因的克隆、表达及功能分析. 林业科学研究, 2016, 29(4): 471-479.
    [6] 张恺恺杨立莹丰美静张林凤陈段芬邱德有杨艳芳 . 曼地亚红豆杉NAC基因家族鉴定及表达分析. 林业科学研究, 2022, 35(2): 97-103. doi: 10.13275/j.cnki.lykxyj.2022.02.011
    [7] 黄国文管天球赵雨云陈莫林刘宏辉 . 油茶转录因子基因CoSOC1-like的克隆和表达分析. 林业科学研究, 2022, 35(2): 129-139. doi: 10.13275/j.cnki.lykxyj.2022.02.015
    [8] 安静万友名马宏刘雄芳张秀姣曹毓蓉李正红 . 地涌金莲MlCYP734A6基因的克隆与表达分析. 林业科学研究, 2021, 34(3): 37-45. doi: 10.13275/j.cnki.lykxyj.2021.03.004
    [9] 章晶晶郭英华赵树堂卢孟柱 . 杨树PtRRI基因的组织特异性表达模式分析. 林业科学研究, 2018, 31(2): 34-40. doi: 10.13275/j.cnki.lykxyj.2018.02.005
    [10] 赵岩秋周厚君魏凯丽江成宋学勤卢孟柱 . 杨树中Ⅰ类KNOX基因结构、表达与功能分析. 林业科学研究, 2018, 31(4): 118-125. doi: 10.13275/j.cnki.lykxyj.2018.04.017
    [11] 朱建峰李万峰杨文华韩素英齐力旺 . 中间锦鸡儿生长发育过程中5个miRNAs及其靶基因的表达模式分析. 林业科学研究, 2013, 26(S1): 45-51.
    [12] 李龙张立峰齐力旺韩素英 . 日本落叶松体细胞胚胎发生相关基因LaSERK1的克隆与表达分析. 林业科学研究, 2013, 26(6): 673-680.
    [13] 张婷丁贵杰文晓鹏 . 马尾松紫色酸性磷酸酶基因PmPAP1的克隆与表达模式分析. 林业科学研究, 2016, 29(6): 797-806.
    [14] 史倩倩周琳李奎王雁 . 云南野生黄牡丹PlbHLH3转录因子基因的克隆与表达. 林业科学研究, 2015, 28(4): 488-496.
    [15] 程华李琳玲许锋王燕程水源 . 银杏过氧化氢酶基因CAT1的克隆及表达分析. 林业科学研究, 2010, 23(4): 493-499.
    [16] 林琳李健李慧玉穆怀志姜静 . 逆境胁迫下柽柳脂质转运蛋白基因 (ThLTP)的克隆与功能初步分析. 林业科学研究, 2012, 25(4): 492-499.
    [17] 孙化雨娄永峰李利超赵韩生高志民 . 毛竹TIPs基因家族成员组织表达模式研究. 林业科学研究, 2016, 29(4): 521-528.
    [18] 徐向东任逸秋张利李煜王丽娟卢孟柱 . 杨树PIF基因家族成员表达模式研究. 林业科学研究, 2018, 31(2): 19-25. doi: 10.13275/j.cnki.lykxyj.2018.02.003
    [19] 徐浩杨克彬朱成磊李英高志民 . 毛竹肉桂酰辅酶A还原酶基因PeCCR功能初步研究. 林业科学研究, 2020, 33(2): 77-84. doi: 10.13275/j.cnki.lykxyj.2020.02.010
    [20] 张俊红吴涛韩素英齐力旺张守攻 . 落叶松体胚发育中12个针叶树特异microRNAs表达分析. 林业科学研究, 2012, 25(4): 411-418.
  • 加载中
图(5) / 表(2)
计量
  • 文章访问数:  5586
  • HTML全文浏览量:  2660
  • PDF下载量:  119
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-10-26
  • 录用日期:  2022-12-27
  • 网络出版日期:  2023-01-04
  • 刊出日期:  2023-06-20

杨树MYC基因家族成员表达模式研究

    通讯作者: 赵树堂, shutang@caf.ac.cn
  • 林木遗传育种全国重点实验室,中国林业科学研究院林业研究所,北京 100091

摘要:  目的 探究bHLH亚家族成员MYC基因在杨树生长发育及对环境响应中的表达模式,为解析茉莉酸信号通路关键基因MYC调控杨树生长发育及抗逆响应提供参考。 方法 利用生物信息学方法鉴定杨树基因组中的MYC基因家族成员并对各成员的基因结构、保守基序进行了系统分析;在此基础上,采用实时荧光定量PCR技术检测了MYC成员在不同组织、不同植物激素响应及逆境胁迫下的表达模式。 结果 毛果杨基因组中含有10个MYC成员,这些成员在进化上相对保守,均具有典型的bHLH结构域,根据进化关系将其分为3个分支。组织表达分析发现,大部分MYC家族成员在根中均具有较高表达量,分枝Ⅱ两对基因在茎中呈相反表达模式。不同激素处理及非生物胁迫条件下MYC基因存在显著表达差异,但基因对之间多存在类似的表达模式。 结论 杨树MYC基因家族成员可能参与不同的生物学过程,研究结果为深入解析杨树MYC基因的功能奠定了基础。

English Abstract

  • 植物激素茉莉酸(JA)是脂肪酸衍生的羟脂素,在植物生长发育和逆境胁迫响应等过程中发挥重要作用[1]。JA通过其信号转导途径调控与多种激素互作或调控下游响应基因参与多种生物学过程。JA信号途径是一个由多基因参与的复杂调控过程[2]。 Coronatine insensitive 1 (COI1)是JA信号的受体,与JA结合后促进转录抑制因子jasmonate ZIM-domain (JAZ)家族蛋白泛素化并被26S蛋白酶体降解,从而解除JAZ对转录因子髓细胞组织增生蛋白(Myelocytomatosis proteins,MYCs)MYC2的抑制作用,最终激活下游响应基因的表达[3]。植物中MYC是JA信号途径中关键转录因子,通过形成COI/JAZs/MYC复合物发挥调控作用,参与多种激素信号转导过程[4]

    MYC转录因子含有保守的bHLH结构域,属于bHLH类转录因子家族。bHLH结构域由碱性域与HLH域组成, N端碱性域具有DNA结合功能,C端HLH域参与bHLH蛋白二聚化形成同源或异源二聚复合体[5]。bHLH类转录因子是转录因子中的超级家族,在拟南芥(Arabidopsis thaliana (L.) Heynh.)中共有158个家族成员,分为26个亚家族,其中,MYC转录因子分布于亚家族III [6],包含8个成员。已有研究证明,MYC转录因子在植物生长发育、次生代谢物的合成、逆境响应、植物激素信号转导及不同激素信号通路之间的相互作用中发挥关键作用[4,7-8]。MYC2是MYC类转录因子中研究最透彻的一个,在JA信号通路中发挥关键调控作用。在拟南芥根系发育过程中,MYC2可以直接结合到维持根尖细胞活性的关键转录因子PLT1PLT2的启动子上抑制其表达,进而抑制拟南芥主根分生区细胞分裂[9]。在蓝光诱导下,拟南芥MYC2/MYC4可以结合NST1启动子激活其表达,进而调控次生壁合成,促进花序轴纤维细胞次生壁加厚[10]。改变红光与远红光比例,光敏色素互作因子PIF4能够与MYC2/MYC4直接相互作用影响MYC蛋白核定位,抑制MYC2/MYC4对NST1的激活进而影响细胞次生壁加厚[11]。在花发育过程中,myc2/3/4三突变体及myc2/3、myc2/4双突变体表现出早开花表型[12],同时伴有花萼发育延迟、花药不能正常开裂或开裂延迟以及花粉粒存活率低等表型[13]。此外,MYC2参与多种植物激素间的调控,包括生长素、细胞分裂素、乙烯等。MYC2可以与乙烯不敏感转录因子EIN3相互作用,抑制EIN3的转录活性,从而抑制HLS1的表达,最终影响拟南芥黄化苗顶端弯钩的形成[14]。在果实成熟方面,MYC2转录因子能通过结合苹果乙烯合成途径中的ACC合成酶基因和ACC氧化酶基因启动子调控乙烯合成进而影响果实成熟[15]。最近研究表明,MYC2可以通过抑制细胞分裂素氧化酶/脱氢酶(CKX)基因来提高细胞分裂素含量,从而抑制不定根的起始[16]。综上所述,MYC类转录因子对植物生长发育及生物和非生物胁迫响应具有重要调控作用。

    尽管MYC基因家族已在拟南芥、水稻(Oryza sativa L.)等模式植物中研究得比较深入,但在林木中其功能研究还鲜有报道。杨树(Populus L.)是重要的人工林树种,而且具有基因组小、适应性强、遗传转化体系成熟等特性,已成为林木分子育种研究的模式树种。因此,本文系统分析了杨树MYC家族成员基因结构、蛋白保守结构域、组织表达特性及不同激素和胁迫响应模式,研究结果为进一步解析杨树MYC基因的生物学功能奠定了基础。

    • 在拟南芥基因组网站TAIR(TAIR - Home Page (arabidopsis.org))数据库中获取MYC基因家族8个成员的氨基酸序列,并以此为靶序列,在Phytozome 网站(https://phytozome-next.jgi.doe.gov/)中通过BLASTP方法比对检索毛果杨基因组数据库V4.1,获得MYC直系同源氨基酸序列。利用ExPASy(http://www.expasy.org/)在线工具分析杨树MYC蛋白基本理化性质。使用WoLEPSORT(WoLF PSORT: Protein Subcellular Localization Prediction Tool (genscript.com))在线工具分析杨树MYC蛋白的亚细胞定位情况。

    • 参考MEGA7.0[17],利用MEGA11中ClustalW程序,将毛果杨、拟南芥MYC蛋白序列进行多重序列比对;采用Neighbor-Joining(NJ)法,模型为p-distance,成对删除,校验参数为1000,构建系统发育进化树。iTOL(https://itol.embl.de/)在线网站进一步完善进化树。

    • 杨树MYC基因结构用GSDS2.0 (http://gsds.cbi.pku.edu.cn/)软件分析。杨树MYC蛋白结构域使用MEME(MEME - Submission form (meme-suite.org))在线分析,最大motif检索数为10,其余为默认参数。

    • 实验所用材料84K杨(Populus alba × P. glandulosa)由林木遗传育种全国重点实验室保存。RNA提取使用(RNA Easy Fast Plant Tissue Kit)试剂盒。分别提取温室培养2个月 84K杨的根、茎、叶以及茎的不同部位的总RNA,反转录合成cDNA第一条链,并以cDNA为模板进行实时荧光定量PCR分析(qRT-PCR)。反应体系为TB Green® Premix EX Taq ™ Ⅱ 5 μL,正向引物(10 μmol·L−1)0.5 μL,反向引物(10 μmol·L−1)0.5 μL,cDNA 1 μL,超纯水3 μL;qRT-PCR程序为:95 ℃ 30 s;95 ℃ 5 s;60 ℃ 30 s;50 ℃ 30 s,45个循环。引物序列见表1

      表 1  实时荧光定量PCR引物序列

      Table 1.  The primers used for quantitative real-time PCR

      引物
      Primer
      序列(5′−3′)
      Primer sequence(5′−3′)
      引物
      Primer
      序列(5′−3′)
      Primer sequence(5′−3′)
      PtrMYC2a-F TATGCCAATGTGTCTGTGATGAATG PtrMYC2a-R TCTTGCGTGTAAAACCGACTCC
      PtrMYC2b-F AATGCTGGGATTCCGTCAACTAC PtrMYC2b-R TGATGATTGTTATTGTTCGCAGTCT
      PtrbHLH14a-F GTGGTCCCGAATGTGTCAAAGATG PtrbHLH14a-R ATGTTCGCATTACCAGAGATTTTTG
      PtrbHLH14b-F TACTCAAGCAAGGTTAGCAACAAGC PtrbHLH14b-R ATCAACCAATCTCTCCAGGTCCAT
      PtrbHLH14c-F GCTCGTGAAACAGAAGAAGCAAAT PtrbHLH14c-R TAGAGGAAAGTCCGAGTCTGAATGC
      PtrbHLH14d-F TGACACTATGGACAATCAAAGCACC PtrbHLH14d-R ATTCACATTATCGGATTGGACCCTA
      PtrbHLH3a-F CAACAAGCACGAATCAAGTTTATGG PtrbHLH3a-R GGACAAGTCATCATTTGGCTGC
      PtrbHLH3b-F GGAGTCTCACCCCGTTTCTACC PtrbHLH3b-R CAGTCTTGAGCAGTGATTTGATGTT
      PtrbHLH17a-F CAAGGGACACCAGGTGAGGTTTAT PtrbHLH17a-R CCTGCTTTTGACCTTGGTAGTGATT
      PtrbHLH17b-F AGGCTTCATTGTTGGGAGATGCTAT PtrbHLH17b-R AATGCTCCCCGACTTTTCCCTT
      PtrActin-F AAACTGTAATGGTCCTCCCTCCG PtrActin-R GCATCATCACAATCACTCTCCGA
    • 激素处理:选择生长4周的84K组培苗,分别进行生长素(10 μmol·L−1 NAA)、细胞分裂素(10 μmol·L−1 6-BA)和茉莉酸(100 μmol·L−1 MeJA)处理。处理时间为:生长素(0.5、1、2、3、6 h)、细胞分裂素(0.5、1、2、3 h)、茉莉酸(1、3、6、9、12 h)。每个处理包含3次生物学重复,每次重复包含3株独立植株,取样位置为顶端向下第4片叶,以未处理的84K杨作为对照。

      非生物胁迫处理:84K杨温室培养4周后,分别进行4 ℃、20% PEG6000及200 mmol·L−1 NaCl处理,处理3、6、12、24 h后取顶端第4片叶,每处理重复3次,每重复3株植株混合取样,以未处理84K杨作为对照。

      使用RNA Easy Fast Plant Tissue Kit试剂盒提取处理后杨树叶片RNA,并通过反转录试剂盒合成单链cDNA。使用Primer Premier 5软件,在杨树MYC基因序列的非保守区设计qRT-PCR引物(表1),引物序列的Tm值为60~62 ℃,并由生工生物公司合成。qRT-PCR使用TB Green® Premix EX Taq ™ Ⅱ(TaKaRa Dalian,China)试剂盒,并根据Roche480实时荧光定量PCR仪说明书进行操作。反应体系为:TB Green® Premix EX Taq ™ Ⅱ 5 μL,上下游引物各0.5 μL,cDNA 1 μL,超纯水3 μL。qRT-PCR所用内参基因为杨树Actin基因(基因序列号:Potri.001G309500)。每个模板设置4个重复,所有实验均进行3次生物学重复,采用2−ΔΔct方法进行定量数据分析。

    • 利用Phytozome网站BLAST工具查找拟南芥MYCs在毛果杨基因组中同源基因,共筛选获得10个杨树MYC基因成员(表2),编码的蛋白质序列长度为466~660 aa,蛋白质分子量介于52.34~72.09 kDa,等电点为5.38~7.72。亚细胞定位预测结果显示:10个MYC蛋白亚细胞定位不同,其中,6个定位在细胞核,2个定位在细胞核和细胞质, 2个定位在细胞核与叶绿体。基于鉴定出的毛果杨PtrMYCs以及拟南芥MYC基因家族编码的氨基酸序列,使用NJ法构建系统发育进化树,结果(图1)显示:与拟南芥MYC基因家族类似,毛果杨MYC基因家族分为3个分枝,但不同分枝成员数量与拟南芥存在较大差异。分枝Ⅰ包含拟南芥MYC2/3/4/5共4个基因,但杨树中只有2个同源基因,分别命名为PtrMYC2aPtrMYC2b;分枝Ⅱ包含拟南芥AtbHLH14及杨树4个同源基因,分别命名为PtrbHLH14a14d;分枝Ⅲ包括拟南芥AtbHLH3AtbHLH17,每个基因在杨树中均有2个同源基因,命名为PtrbHLH3aPtrbHLH3bPtrbHLH17aPtrbHLH17b

      表 2  杨树MYC基因家族成员基本信息

      Table 2.  Basic information of MYC family members in Populus

      基因名称
      Gene
      基因ID
      Gene number
      基因在染色体位置
      Location
      蛋白质长度
      Protein length/aa
      等电点
      PI
      分子量
      MW/kDa
      亚细胞定位
      Subcellular localization
      PtrMYC2a Potri.001G142200 Chr01( + )11649287-11651987 648 5.38 71.11 细胞核
      PtrMYC2b Potri.003G092200 Chr03( − )11839709-11842763 660 5.49 72.09 细胞核
      PtrBHLH14a Potri.001G083500 Chr01( − )6613720-6615721 472 7.17 52.73 细胞核/细胞质
      PtrBHLH14b Potri.003G147300 Chr03( + )16158005-16160025 466 7.72 52.34 细胞质/细胞核
      PtrBHLH14c Potri.002G176900 Chr02( + )13775139-13776876 492 5.40 54.70 细胞核
      PtrBHLH14d Potri.014G103700 Chr14( + )6950999-6952650 494 5.91 54.87 细胞核/叶绿体
      PtrBHLH3a Potri.002G042000 Chr02( − )2756159-2758458 504 6.59 54.92 细胞核/叶绿体
      PtrBHLH3b Potri.005G221100 Chr05( + )22329285-22332096 505 6.13 55.44 细胞核
      PtrBHLH17a Potri.002G172101 Chr02( − )13241064-13244032 619 6.14 68.25 细胞核
      PtrBHLH17b Potri.014G099700 Chr14( − )6539405-6542974 620 6.16 68.02 细胞核

      图  1  杨树和拟南芥MYC基因家族系统进化树

      Figure 1.  Phylogenetic tree of MYC family members between Populus and Arabidopsis

    • 为了了解PtrMYC基因的结构特征,对其基因结构及保守基序进行分析,发现杨树所有MYC基因均不含内含子(图2A)。使用MEME软件对杨树MYC蛋白进行保守结构域分析,结果表明:所有MYC蛋白的结构域在进化上是高度保守的,Motif1、2、3、4、5、6、7、8存在于所有家族成员中,其中,Motif1、3、7和2、4、5、6、8为该家族基因保守的特征结构域即bHLH和bHLH-MYC_N。另外,有的家族分枝具有或缺失某些蛋白基序,如Motif9仅存在于分枝Ⅱ的家族成员中,Motif10仅存在于分枝Ⅰ和Ⅲ的成员中,大多数同源配对成员有高度相似的结构域组成,提示可能在基因功能上存在冗余(图2B)。

      图  2  杨树MYC家族成员基因结构(A)和蛋白保守结构域(B)分布

      Figure 2.  Gene structure (A) and distribution of conversed motifs (B) of PtrMYC family members in Populus

    • 基因组织特异性表达模式研究对预测基因功能具有重要的作用。为了探讨PtrMYC基因家族在杨树生长发育中的功能,本研究利用qRT-PCR方法分析了MYC基因各成员在不同器官和组织的表达特异性,结果(图3)显示:不同MYC家族成员具有不同的组织表达特异性,但基因对之间多存在类似的表达模式。除PtrbHLH3a外,其余9个PtrMYC基因在根中均有较高表达量,说明大部分MYC家族成员可能在根系发育过程中发挥作用;在从顶端到基部的不同茎段中,PtrbHLH14aPtrbHLH14b这一组基因对表达量呈现先上升后下降的趋势,说明它们可能在茎由初生生长向次生生长发育过程中具有重要调节作用;PtrbHLH14cPtrbHLH14d这组基因对表达量则呈现下降的趋势,说明可能在顶端分生组织发育中具有重要的作用。杨树MYC家族成员在不同器官和组织存在显著表达差异,说明不同成员可能在杨树生长发育不同生物学过程中发挥调控作用。

      图  3  杨树MYC基因组织表达分析

      Figure 3.  Expression analyses of MYCs in Populus

    • 为进一步了解杨树MYC家族基因对不同植物激素的响应情况,本研究利用qRT-PCR方法检测了其对生长素( NAA)、细胞分裂素(6-BA)、茉莉酸( MeJA)响应的表达模式(图4)。在NAA、6-BA、MeJA处理后,杨树MYC基因家族各成员表达模式各不相同。PtrbHLH14aPtrbHLH14b基因表达量受生长素诱导0.5 h后开始下调,表明生长素可以快速抑制特异MYC基因的表达。细胞分裂素处理0.5 h后, PtrbHLH14aPtrbHLH14bPtrbHLH14d基因表达量显著上调,说明它们可以快速响应细胞分裂素处理。大部分杨树MYC基因都能响应茉莉酸信号,茉莉酸处理后, PtrMYC2aPtrMYC2bPtrbHLH14cPtrbHLH14dPtrbHLH17aPtrbHLH17b基因表达量显著上调。值得注意的是,PtrbHLH14aPtrbHLH14b基因对能够响应生长素和细胞分裂素处理,但表达模式呈相反趋势,说明它们可能同时参与生长素与细胞分裂素信号通路。综上所述,MYC家族成员可能广泛参与不同激素调控植物生长发育的多个过程。

      图  4  生长素、细胞分裂素、茉莉酸处理后杨树MYC基因的表达模式

      Figure 4.  Expression patterns of PtrMYCs under NAA、6-BA、MeJA treatment

    • 为揭示杨树MYC基因对非生物胁迫的响应机制,本研究利用qRT-PCR方法检测了该家族基因在4 ℃、PEG6000和NaCl胁迫下的表达模式(图5)。在200 mmol·L−1 NaCl、20% PEG6000及4 ℃低温处理后,杨树MYC基因家族各成员表达都受到不同程度的诱导,且大部分基因表达量在处理后12 h达到高峰。低温处理可以诱导6个杨树MYC基因上调表达,其中PtrbHLH17aPtrbHLH17b两个基因受低温诱导上调表达最多,且在12 h达到峰值;PtrbHLH14aPtrbHLH14b下调表达,说明不同MYC基因响应低温表达方式不同。在NaCl处理之后,PtrMYC2bPtrbHLH14bPtrbHLH3aPtrbHLH3bPtrbHLH14aPtrbHLH14d基因的表达量先上调再下调,而PtrMYC2aPtrbHLH14cPtrbHLH17b基因则呈现相反表达模式。干旱处理后,PtrbHLH14aPtrbHLH14bPtrbHLH14dPtrbHLH3a基因表达量显著上调,其中,PtrbHLH3a表达量在处理后12 h时达到高峰,其它3个基因在9 h时达到高峰,说明不同MYC基因响应干旱时间不同。另外,同一分支基因响应相同非生物胁迫时具有相似表达模式,说明基因功能存在冗余。综上所述,MYC家族成员在不同非生物胁迫中表达情况存在差异,推测在植物响应逆境胁迫中发挥不同作用。

      图  5  胁迫条件下MYC基因的表达模式

      Figure 5.  Expression patterns of PtrMYCs under 4 ℃, NaCl and PEG6000 treatment

    • 林木在长期进化过程中形成了自己独特的生物学特征和环境适应性,而基因组学和基因功能研究发现,不同物种同源基因可能在功能上存在较大差异,因此,在模式树种杨树中开展基因功能研究对揭示林木特有的生物学过程具有重要意义[18]MYC基因家族是bHLH基因家族一个亚家族,对植物的生长发育及环境适应性具有重要作用。目前,在拟南芥中对MYC2基因进行了较为深入的研究,发现该基因不仅能参与调控根、叶、种子等器官的发育活动,并且还以调控开关的角色存在于多种激素通路互作过程中[5]。然而,模式树木毛果杨中,MYC家族基因基本特征及其逆境胁迫响应特性尚无报导。本研究根据拟南芥MYC基因家族氨基酸序列,在毛果杨基因组中通过序列比对获得了杨树MYC基因家族成员,并对其基因结构、组织表达特异性及对激素和逆境响应表达模式等进行了系统分析。在毛果杨基因组中共鉴定出10个MYC成员,系统进化树分析发现该基因家族分为3个分枝,其中,分枝Ⅰ包含2个成员,分枝Ⅱ有4个成员,分枝Ⅲ有4个成员,这与拟南芥中各分枝成员数量相差较大,说明MYC基因家族不同分枝在进化过程中可能存在基因扩张或丢失现象。蛋白结构域分析发现,PtrMYC蛋白含有多个基序,10个Motif中有8个Motif存在于所有MYC成员中(图2B),表明PtrMYC成员之间在蛋白结构上高度保守。

      基因组织表达特异性对进一步探究其功能具有重要作用。本研究中,杨树PtrMYC在不同器官和组织存在明显表达差异,表明它们可能在不同生长发育过程中发挥作用。拟南芥中研究发现,MYC2抑制主根分生区细胞分裂[9],同时促进侧根形成[19]。另外,MYC2可以通过调控细胞分裂素含量抑制不定根起始[16],说明MYC2转录因子对根系发育的调控是多方面的。组织表达分析表明,PtrMYC家族9个成员在根部表达水平较高,暗示它们可能参与杨树根系不同发育过程。在叶片发育过程中,MYC2能通过抑制生长素合成进而抑制叶脉发育[20]。Qi等[13]研究表明,MYC2/MYC3/MYC4在激活叶片衰老相关基因SAG29诱导叶片衰老中存在功能冗余,而bHLH03bHLH13bHLH14bHLH17SAG29的启动子结合抑制其表达,拮抗MYC2/MYC3/MYC4诱导的叶片衰老[21]。另外,MYC2-Dof2.1形成前馈机制促进叶片衰老[22]。本研究中发现,PtrbHLH3aPtrbHLH3bPtrbHLH17aPtrbHLH17b等在叶片中高丰度表达,可能参与杨树叶片发育和衰老过程的调控。

      植物激素之间的相互作用组成一个非常复杂的信号调控网络,在植物生长发育和逆境响应过程中起着非常重要的调控作用[23]MYC基因在不同激素之间互作过程中具有重要作用。在不定根和叶片发育过程中,MYC家族基因能够抑制细胞分裂素代谢和生长素合成影响不定根起始和叶脉发育 [16, 20]。拟南芥中,MYC2可以与其靶基因ERD1的启动子结合,激活其表达,进而激活JA诱导的脱水应激反应[24]。此外,茉莉酸和生长素协同作用调控拟南芥热形态建成[25]。在苹果愈伤组织中过表达MdMYC2增加了MdCIbHLH1MdCBF1MdCBF2MdCBF3的表达水平,导致苹果抗冻耐受性显著提高[26]。在番茄中MYC2激活乙烯响应因子SlERF.B8的转录,进而增强JA信号,增加了番茄耐寒性[27]。本研究发现,不同PtrMYC成员响应激素和非生物胁迫方式不同,一些MYC基因既受激素诱导,又响应某种非生物胁迫,如PtrMYC14aPtbHLH14b表达可以被生长素抑制,也可以被细胞分裂素诱导,同时还可以响应盐和干旱胁迫;而PtrMYC2b表达能够被茉莉酸诱导,同时还可以响应低温和盐胁迫,提示MYC基因家族可能参与不同植物激素相互作用,进而广泛参与到植物生长发育和非生物胁迫响应中,且不同成员可能参与不同生物学过程,但其具体调控机制尚需进一步研究。

    • 本文系统分析了杨树MYC基因家族的基因和蛋白结构、组织表达特异性以及响应激素和逆境的表达模式。杨树基因组中共有10个MYC成员,均含有bHLH转录因子典型结构域;组织表达分析发现,杨树MYC基因家族成员表达模式各异,可能参与了不同的生物学过程的调控。在不同的植物激素处理下,杨树MYC家族成员表达模式多样,表明不同的MYC基因可能参与到不同的植物激素信号通路。在不同的非生物胁迫下,杨树MYC基因家族各成员的表达情况也不同,表明杨树MYC基因在不同的非生物胁迫下发挥不同的功能。本研究为揭示MYC基因在杨树生长发育和胁迫响应过程中的基因功能奠定了基础。

参考文献 (27)

目录

    /

    返回文章
    返回