丹参酮生物合成分子调控机制的研究进展

李 珅，林爱真，杨 媛，沈亚芳，饶 盈，羊 健，刘云辉，王 洋，周 伟

（1.浙江农林大学 农业与食品科学学院 浙江省农产品品质改良技术研究重点实验室，浙江 杭州311300； 2.浙江省泰顺县西旸镇农业公共服务中心，浙江 泰顺 325509；3.临安天目山林场，浙江 杭州311311；4.浙江省农业科学院 农业部/浙江省植保生物技术重点实验室，浙江 杭州 310021）

丹参Salvia miltiorrhiza是传统的中药材，其主要药用成分包括脂溶性的丹参酮类化合物和水溶性的酚酸类化合物。其中丹参酮是一类松香烷二萜化合物。20世纪30年代末，有学者首次发现丹参酮ⅡA和丹参酮ⅡB是天然的抗氧化剂。随后几十年，又连续发现了丹参酮Ⅰ、隐丹参酮、二氢丹参酮Ⅰ、异丹参酮、去甲基丹参酮等脂溶性有效成分。其药理作用主要表现在抗菌消炎、抗氧化、抗肿瘤、抗心肌缺氧、激素样活性等方面，并能够改善学习记忆能力，在临床上主要用于心脑血管疾病的治疗［1-2］。丹参有效成分产量较低，可利用资源有限，使其应用推广受到很大的限制。因此，筛选获得高丹参酮产量的优良品种或是通过其他方法解决药源问题是很有必要的。笔者从丹参酮生物合成途径及其分子调控机制等方面展开综述，同时总结了以往研究中存在的问题，并结合自身对研究前景进行了分析，为后续深入研究提供参考。

1 丹参酮的生物合成途径

丹参酮是一类二萜类化合物，其代谢合成方式有2条途径（图1），分别为位于细胞质中的甲羟戊酸途径 （mevalonate,MVA）［3］和位于质体中的非甲羟戊酸途径（2-C-methyl-d-erythritol-4-phosphate,MEP）［4］。MVA合成途径的上游基因包括AACT，HMGS，HMGR，MK，PMK和MDC；MEP合成途径上游基因包括DXS，DXR，CMS，CMK，MCS，HDS和 IDS；中游合成酶基因包括IPPI，GPPS，FPPS和GGPPS；下游合成酶基因包括CPS，KSL，CYP76AH1，CYP76AH3，CYP76AK1以及未知的P450s基因；在下游合成途径中柯巴基焦磷酸（copalyl diphosphate,CPP）经类贝壳杉烯合酶（KSL）催化形成丹参酮二烯（miltiradiene），丹参酮二烯经细胞色素 P450家族蛋白CYP76AH1催化合成铁锈醇（ferruginol），铁锈醇经过CYP76AH3和CYP76AK1的多步催化生成11,20-二羟基铁杉醇和11,20-二羟基柳杉酚，最后在其他CYP家族或其他酶作用下产生次丹参酮（miltirone），进而合成丹参酮类成分［5］。目前，从铁杉醇合成丹参酮之间的多步催化反应步骤仍不清楚，需要进一步的研究。

图1 丹参酮代谢合成路径Figure 1 Biosynthetic pathway involved in the biosynthesis of tanshinone in Salvia miltiorrhiza

2 诱导子对丹参酮生物合成的调控作用

诱导子能诱发植物次生代谢过程中酶的活性，从而增加次生代谢物的产量，有时可以诱导出新的代谢化合物。依据来源可划分为生物诱导子和非生物诱导子两大类。生物诱导子是指植物体在防御过程中为对抗微生物入侵而产生的防御性物质，如分生孢子、分解细胞壁的酶、细胞壁的碎片。非生物诱导子是指所有非植物细胞的固有成分，但能触发植物细胞形成抗毒素信号的物质［6］。

2.1 生物诱导子对丹参酮合成的影响

研究人员围绕生物诱导子对丹参酮合成的影响做了许多研究。晏琼等［6-7］利用酵母提取物（YE）、寡聚半乳糖醛酸和真菌诱导子等3种生物诱导子分别处理丹参毛状根，发现都能显著提高丹参酮的产量。施加100 g·L-1的真菌诱导子，总丹参酮的产量提升4.70倍；施加100 g·L-1的YE，总丹参酮的产量提升3.20倍；施加100 g·L-1的寡半乳糖醛酸，总丹参酮的产量提升4.20倍。GE等［8］考察了不同浓度的β-氨基丁酸（BABA）和YE单独或是协同对丹参酮代谢合成的影响，与对照相比最多可提高4.50倍。BABA+YE协同作用，在BABA处理3 d后再加入YE，对丹参酮的积累更加有效。YAN等［9］发现密旋链霉菌Streptomyces pactum处理丹参毛状根后，GGPPS，HMGR，DXS，DXR等4个关键酶基因的表达量显著提升，总丹参酮产量提高12.61倍［9］。研究发现，丹参毛状根中丹参酮ⅡA和隐丹参酮的积累和产量受YE诱导。丹参酮的产量经诱导处理6 d后测得的最大值为0.64 mg·g-1，其中HMGS，DXR，DXS2，CMK，IPPI，CPS等6个关键酶基因的表达量显著上调。YE和银离子（Ag+）共同诱导处理表明，丹参酮ⅡA和隐丹参酮的产量显著增加。丹参酮诱导处理9 d后测得的丹参酮的最大值为2.08 mg·g-1，其中AACT，HMGS，DXR，IPPI，CMK，FPPS，GGPPS，CPS等8个关键酶基因表达显著上调［10］。 上述研究表明，生物诱导子对于丹参酮的合成具有促进作用。因此，新型高效生物诱导子的挖掘和开发利用将是后续重要研究方向之一。

2.2 非生物诱导子对丹参酮合成的影响

非生物诱导子能引起或者促进植物细胞与防御机制有关的次生代谢产物的生物合成。目前，已证实的非生物诱导子约有100多种，主要包括茉莉酸甲酯（methyljasmonate，MeJA），茉莉酸（jasmonate，JA），水杨酸（alicylic acid，SA），重金属盐类和稀土元素。郭肖红等［11］研究发现，低浓度的铜离子（Cu2+）， 镁离子（Mg2+）， 锌离子（Zn2+）， 亚铁离子（Fe2+）， 锰离子（Mn2+）以及稀土元素均能促进丹参酮ⅡA的合成。KAI等［10］研究了0.10 mmol·L-1MeJA诱导处理丹参毛状根后丹参酮ⅡA和隐丹参酮的积累变化。处理9 d后，HMGR，DXR，DXS2，GGPPS，CPS和KSL等基因的表达显著上调，丹参酮的产量达到最大值 0.93 mg·g-1干质量（DW）， 相比于对照组 0.16 mg·g-1（DW）增加了约 5.78 倍。 HAO 等［12］发现SA处理丹参毛状根36 h后，丹参酮产量上升1.63倍。银离子处理毛状根后，HMGR，DXS2，IPPI，GGPPS等基因的表达显著上调；丹参酮ⅡA和隐丹参酮的积累逐渐增加，处理9 d后测得的丹参酮的最大产量为 0.35 mg·g-1， 相比于对照组的 0.16 mg·g-1（DW）增加约 2.20 倍。 晏琼等［6-7］利用银离子（Ag+），钴离子（Co2+）和α-氨基异丁酸（AIB）等3种非生物诱导子分别处理丹参毛状根，通过检测丹参酮的产量发现，此3种非生物诱导子都能提高丹参酮的积累，但诱导能力比生物诱导子弱；其中50.00 mmol·L-1钴离子获得了最高的丹参酮产量（1.75 mg·g-1），是对照组产量的4.10倍；30.00 mmol·L-1银离子处理后总丹参酮的产量提高了2.10倍；50.00 mmol·L-1AIB处理后总丹参酮产量提高了2.30倍。ZHOU等［13］利用金属元素0.01 mmol·L-1镧（La）对丹参毛状根进行诱导处理，结果表明丹参酮代谢途径上关键酶基因包括HMGR,HMGS,KSL,DXR,CMK的表达量均上调，丹参酮Ⅰ，丹参酮Ⅱ和隐丹参酮的产量均升高。由此可推测镧可能通过促进丹参酮代谢途径上游关键酶基因的表达而增加丹参酮的积累。房翠萍［14］研究了外源激素对丹参毛状根生长和丹参酮ⅡA生物合成的影响，结果显示在1/2MS培养基中单独施加2.00 mg·L-16-苄氨基腺嘌呤（6-BA）能使丹参酮ⅡA产量增加3.01倍；同时添加3.00 mg·L-16-BA与0.20 mg·L-1NAA能使丹参酮ⅡA产量增加0.93倍；同时添加1.00 mg·L-16-BA与2 mg·L-1赤霉素（GA）能使丹参酮ⅡA浓度增加1.69倍；同时添加1.00 mg·L-1激动素（KT）与1.00 mg·L-1吲哚丁酸（IBA）能使丹参酮ⅡA产量增加0.21倍；此外还研究了pH值对丹参毛状根生长和丹参酮ⅡA生物合成的影响，结果发现培养基pH值为pH 9时，丹参酮产量最高。YANG等［15］利用聚乙二醇（PEG）和脱落酸（ABA）对丹参毛状根进行诱导处理，发现PEG和ABA通过增加毛状根中活性氧含量促进丹参酮的代谢合成，推测PEG和ABA可能通过ABA途径诱导丹参毛状根中内源茉莉酸（JA）的活性，从而刺激丹参酮的代谢积累［15-16］（表1）。

表1 诱导子对丹参酮代谢合成的影响Table 1 Effect on the biosynthesis of tanshinone induced by elicitors in Salvia miltiorrhiza

3 转录因子对丹参酮代谢合成的调控作用

转录因子又称为反式作用元件，是能够特异结合真核基因启动子区域顺式作用元件的脱氧核糖核酸（DNA）结合蛋白，通过蛋白与基因之间、蛋白与蛋白之间的相互作用激活或抑制转录［17］。近年来，随着对植物转录因子研究和认识的不断深入，转录调控被认为是一种能有效调控植物次生代谢途径的新手段。超表达转录因子能够有效地克服多基因共转化存在的弊端，全面调节相关基因的表达，对整个代谢途径的调控效果要优于单纯提高单个或多个结构基因的表达。目前，研究人员对转录因子在丹参次生代谢调控方面进行了深入研究，研究较多的转录因子有基本螺旋-环-螺旋（basic helix-loop-helix，bHLH）、植物中髓细胞组织增生蛋白（myelocytomatosis，MYCs），WRKY和JAZ等几个类型（表2）。

3.1 bHLH类转录因子

植物bHLH转录因子家族因其成员含有高度保守的bHLH区域而被命名为bHLH转录因子。bHLH结构域由50～60个左右氨基酸组成，包含10～15个氨基酸组成的碱性氨基酸区和40个左右氨基酸构成的HLH区［18］。50%以上的植物碱性氨基酸区含有高度保守的His5-Glu9-Arg13序列，这使得碱性氨基酸区能够识别 DNA 上 E-box（5′-CANNTG-3′）和 G-box（5′-CACGTG-3′）位点， 并与之结合［19-20］。 由于 HLH结构域中疏水氨基酸之间常常互作从而促进蛋白二聚体的形成，因此bHLH蛋白常以二聚体的形式发挥作用［21］。目前已报道的植物bHLH家族转录因子已有很多，其中拟南芥Arabidopsis thaliana中报道的bHLH家族成员有164个，水稻Oryza sativa中bHLH家族成员有180个［22］，而烟草Nicotiana tabacum和葡萄Vitis vinifera中bHLH转录因子分别有190个和191个以上［23-24］。通过转录组测序，127个丹参bHLH类转录因子已被分离，其中7个成员被发现能响应甲基茉莉酸 （MeJA）信号诱导。SmbHLH37，SmbHLH74和SmbHLH92被证实与丹参萜类物质的代谢合成密切相关［25］。由此推测，bHLH类转录调控因子在JA信号途径中调控丹参酮和丹酚酸的生物合成具有重要的作用。此外，丹参转录因子SmbHLH1和SmbHLH93的全长基因被相继克隆，分析表明此2个基因的表达与丹参酮的代谢合成密切相关［26-27］，这些研究为从分子水平全面解析bHLH类转录因子调控丹参酮的代谢合成机制奠定了扎实的工作基础。

3.2 MYC类转录因子

MYC类转录因子是植物茉莉酸类激素响应途径中的核心转录因子。MYC具有多种调节功能，广泛存在于动植物中。MYC类转录因子中MYC2研究最为深入。MYC2转录因子含有bHLH保守结构域，属于bHLH类转录因子家族成员。MYC2在N端含有1个JID结构域，与JAZ结合相关［28］。ZHOU等［29］研究发现，丹参SmMYC2a能和丹酚酸代谢途径的结构基因SmHCT6，SmCYP98A14的E-BOX互作；而SmMYC2b能和SmCYP98A14的E-BOX互作，直接调控丹酚酸而间接调控丹参酮的代谢合成。王晓荣［30］研究发现，SmMYC2-like能和SmJAZ3互作，在丹参毛状根中过表达AtMYC2后，丹参酮代谢合成关键酶基因SmGGPPS的表达上调，总丹参酮产量达16.00 mg·g-1（DW），高出对照5.40倍。而过表达SmMYC2-like后，总丹参酮产量增加3.50倍。

3.3 WRKY类转录因子

WRKY转录因子家族是高等植物十大转录因子家族成员之一。WRKY转录因子的保守结构域约由60个氨基酸残基组成，其中靠近氨基酸N端的7个保守氨基酸残基WRKYGQK为WRKY结构域的核心序列，它们的变异会导致DNA结合活性减弱或者丧失；羧基C末端的锌指类似结构在植物的进化中可能起到重要的作用［31-32］。

在植物进化过程中，WRKY转录因子家族成员不断扩大。 LI等［33］克隆了61个WRKYs基因，发现其中42个WRKYs成员能够响应MeJA和银离子（Ag+）胁迫，与丹参酮的代谢合成密切相关。郝小龙［34］研究发现，在丹参毛状根中过表达SmWRKY3后，丹参酮合成关键酶基因SmCPS表达上调，总丹参酮产量达3.98 mg·g-1（DW），相比于对照提高了1.83倍；而抑制表达SmWRKY3后，SmCPS的表达下调，总丹参酮产量仅为1.98 mg·g-1（DW）。因此，WRKY转录因子能够促进丹参酮的合成，是丹参酮合成的正向调控因子。

3.4 JAZ类转录因子

含有ZIM结构域的基因称为JASMONATE ZIM-DOMAIN（JAZ）基因。ZIM结构域一般位于JAZ蛋白的中间，主要由28个氨基酸组成；在N端含有（TIF［F/Y］XG）结构域，在C末端有2个不变的丙氨酸。Jas结构域（SLX2FX2KRX2RX5PY）位于JAZ蛋白的C端，在JAZ蛋白中极其保守。研究发现，拟南芥的SCFCOI1复合体依赖JA-Ile，并通过Jas结构域与JAZ1，JAZ3，JAZ9，JAZ10蛋白结合；Jas结构域是JAZ蛋白核定位信号的关键结构。JAZ参与JA信号介导的次生代谢物的合成在拟南芥中研究得比较透彻。JAZ在拟南芥中有12个成员，其中至少8个JAZ蛋白成员都与MYC2相互作用。拟南芥MYC2的同源基因MYC2-LIKE成员JAM1，JAM2和 JAM3也参与JA信号途径的次生代谢合成（表2）；其中JAM1和JAM2分别与6个不同的JAZ成员相互作用，而JAM3却不能与任何JAZ蛋白成员互作。近期丹参JAZ蛋白参与JA信号介导的次生代谢物的合成机制研究也取得一些进展。在丹参中超表达JAZ3和JAZ9基因能使毛状根中丹参酮的产量急剧降低；酵母双杂试验证实SmJAZ9蛋白能与bHLH类转录因子拟南芥的AtMYC2互作，而SmJAZ3不能与AtMYC2互作，推测SmJAZ9通过与bHLH类转录因子互作，SmJAZ3通过与其他bHLH类转录因子互作而调控丹参酮的生物合成［5］。研究证实将SmJAZ8遗传转化丹参，超表达SmJAZ8的转基因毛状根系表现出对JA信号不敏感的表型；SmJAZ8的超表达导致SmJAZ1，SmJAZ2和SmJAZ3的表达量急剧下调，从而提高丹参酮和花青素的代谢合成；究其原因发现，SmJAZ8蛋白跟拟南芥的AtJAZ8类似，缺少LPIARR的结构域，而此结构域是JA信号诱导JAZ蛋白降解释放正向调节因子的关键区域。可见拟南芥JAZ蛋白在JA信号介导的次生代谢合成作用机制的解析为在丹参中研究JAZ调控丹参酮和丹酚酸的合成提供了很好的参考。丹参JAZ家族其他成员的功能有待于进一步地解析。

表2 转录因子对丹参酮代谢合成的影响Table 2 Effect on the biosynthesis of tanshinone regulated by transcription factors in Salvia miltiorrhiza

4 可变剪接的分子调控作用

真核生物基因转录产生的mRNA前体大多数只按一种方式进行剪接，产生一种成熟的mRNA分子，只翻译成一种蛋白质。但有的基因mRNA前体通过不同的剪接方式（按不同的剪接位点）产生不同的mRNA剪接异构体，此过程称为可变剪接。RNA可变剪接不涉及遗传信息的永久性改变，是真核基因表达调控的重要手段，是调节基因表达和产生蛋白质多样性的重要机制，是导致人类基因和蛋白质数量变异的重要原因。

现如今可变剪接作为代谢和发育另一个重要的调节机制已在许多植物中被鉴定，如拟南芥［36-37］，大豆 Glycine max［38］， 蔓越橘 Brachypodium distachyon［39］和水稻［40］。 而在丹参中， 利用短读高通量测序技术（NGS）和单分子实时（SMRT）长读测序技术，结合对根部各组织进行测序发现，KSL1，HDR，HMGR，AACT3，MK和 PMK基因存在着选择性剪接调控的现象［41］。由此可见，丹参酮代谢合成途径中关键合成酶基因的可变剪接对丹参酮的代谢合成起到了重要的调节作用。某些转录调控因子是否也存在可变剪接的调控方式有待进一步研究。

5 总结与展望

丹参这种传统的中草药具有重要的药用价值。近年来野生丹参品质严重退化，濒临灭绝，导致丹参原药的供应极不稳定，解决药源问题已成为研究的热点。目前，对丹参酮次生代谢合成的研究重点主要包括以下几个方面：①丹参酮代谢合成途径的解析。分离丹参酮代谢合成酶或是修饰酶，通过底物与产物的鉴定，解析丹参酮次生代谢化合物的合成途径。②转录因子及其调控机制和调控网络的解析。转录因子可同时调控代谢途径上的1个或是多个基因的表达，在丹参酮次生代谢合成中起重要的调控作用。为阐明丹参酮代谢合成的调控机制，利用诱导子处理丹参，结合转录组学、蛋白质组学和代谢组学揭示候选转录因子的表达与产物代谢流的分配规律，以揭示转录因子间的表达与互作机制对丹参酮合成代谢的调控规律。③丹参酮基因工程的研究。在代谢途径和调控网络充分解析的基础上，实现代谢调控的可预测性。如采用时空特异性表达的启动子、多基因共转化技术，实现丹参有效成分在特异部位或是全株系大量积累和特定性状改良的目标。

随着研究的深入，进一步从转录调控水平全面解析丹参酮的代谢分子调控机制对于利用现代基因工程手段培育高品质丹参新品种/系具有重要的指导意义，也具有广阔的应用前景。因此，应从以下几个方面开展研究工作：①丹参酮代谢合成关键调控因子可变剪接的挖掘与分子机制解析。②非编码RNA在丹参酮代谢合成过程中的调控作用。③丹参酮合成生物学的研究。通过工程酵母生产丹参酮前体物，并利用有机半合成技术体外合成丹参酮，这将改变丹参药源生产方式，为丹参系列药物走向国际市场奠定工作基础，无疑将具有重要的应用前景。

［1］ CHEN Xiuping,GUO Jiajie,BAO Jiaolin,et al.The anticancer properties of Salvia miltiorrhiza Bunge （Danshen）：a systematic review ［J］.Med Res Rev,2014,34（4）：768-794.doi：10.1002/med.21304.

［2］ WANG Xin,LEE W,ZHOU Xuelin,et al.A pharmacodynamic-pharmacokinetic （PD-PK） study on the effects of danshen （Salvia miltiorrhiza） on midazolam,a model CYP3A probe substrate,in the rat ［J］.Phytomedicine,2010,17（11）：876-883.doi：10.1016/j.phymed.2010.05.007.

［3］ LANGE B,RUJAN T,MARTIN W,et al.Isoprenoid biosynthesis：the evolution of two ancient and distinct pathways across genomes ［J］.Proc Natl Acad Sci USA,2000,97（24）：13172-13177.doi：10.1073/pnas.240454797.

［4］ ROHMER M,KNANI M,SIMONIN P,et al.Isoprenoid biosynthesis in bacteria：a novel pathway for the early steps leading toisopentenyldiphosphate ［J］.Biochem J,1993,295（2）：517-524.doi：10.1042/bj2950517.

［5］ SHI Min,ZHOU Wei,ZHANG Jianlin,et al.Methyl jasmonate induction of tanshinone biosynthesis in Salvia miltiorrhiza hairy roots is mediated by JASMONATE ZIM-DOMAIN repressor proteinsb ［J］.Sci Rep,2016,6（1）：20919.doi：10.1038/srep20919.

［6］ 晏琼,胡宗定,吴建勇.生物与非生物诱导子协同作用对丹参毛状根培养生产丹参酮的影响［J］.中国中药杂志， 2006， 31（3）： 188-191.YAN Qiong,HU Zongding,WU Jianyong.Synergistic effects of biotic and abiotic elicitors on the production of tanshinones in Salvia miltiorrhiza hairy root culture ［J］.China J Chin Mater Med,2006,31（3）：188-191.

［7］ 晏琼，胡宗定，吴建勇.生物和非生物诱导子对丹参毛状根培养生产丹参酮的影响［J］.中草药，2006，37（2）： 262-265.YAN Qiong,HU Zongding,WU Jianyong.Influence of biotic and abiotic elicitors on production of tanshinones in Salvia miltiorrhiza hairy root culture ［J］.China Tradit Herb Drug,2006,37（2）：262-265.

［8］ GE Xiuchun,WU Jianyong.Induction and potentiation of diterpenoid tanshinone accumulation in Salvia miltiorrhiza hairy roots by β-aminobutyric acid ［J］.Appl Microbiol Biot,2005,68（2）：183-188.doi：10.1007/s00253-004-1873-2.

［9］ YAN Yan,ZHANG Shuncang,YANG Dongfeng,et al.Effects of Streptomyces pactum Act12 on Salvia miltiorrhizahairy root growth and tanshinone synthesis and its mechanisms ［J］.Appl Biochem Biotechnol,2014,173（4）：883-893.doi：10.1007/s12010-014-0876-4.

［10］ KAI Guoyin,LIAO Pan,XU Hui,et al.Molecular mechanism of elicitor-induced tanshinone accumulation in Salvia miltiorrhiza hairy rootcultures ［J］.Acta Physiol Plant,2012,34（4）：1421-1433.doi：10.1007/s11738-012-0940-z.

［11］ 郭肖红，高文远，陈海霞，等.金属离子对丹参酮ⅡA和原儿茶醛生物合成的影响［J］.中国中药杂志，2005， 30（12）： 885-888.GUO Xiaohong,GAO Wenyuan,CHEN Haixia,et al.Effects of mineral cations on the accumulation of tanshinoneⅡA and protocatechuic aldehyde in the adventitious root culture of Salvia miltiorrhiza ［J］.China J Chin Mater Med,2005,30（12）：885-888.

［12］ HAO Xiaolong,SHI Min,CUI Lijie,et al.Effects of methyl jasmonate and salicylic acid on tanshinone production and biosynthetic gene expression in transgenic Salvia miltiorrhiza hairy roots ［J］.Biotechnol Appl Biochem,2015,62（1）：24-31.doi：10.1002/bab.1236.

［13］ ZHOU Jie,FANG Lei,WANG Xiao,et al.La dramaticaly enhances the accumulation of tanshinones in Salvia miltiorrhiza hairy root cultures ［J］.Earth Sci Res,2012,2（1）：187-192.doi：http：//dx.doi.org/10.5539/esr.v2n1p187.

［14］ 房翠萍.丹参多倍体诱导、毛状根培养及其丹参酮产量提高的研究［D］.合肥：安徽农业大学，2011.FANG Cuiping.Polyploidy Induction,Hairy Root Culture and Improve the Yield of Tanshinone in Salvia miltiorrhiza［D］.Hefei：Anhui Agriculture University,2011.

［15］ YANG Dongfeng,SHENG Dongfeng,DUAN Qimei,et al.PEG and ABA trigger the burst of reactive oxygen species to increase tanshinone production in Salvia miltiorrhiza hairy roots ［J］.J Plant Growth Regul,2012,31（4）：579-587.doi：10.1007/s00344-012-9268-6.

［16］ YANG Dongfeng,MA Pengda,LIANG Xiao,et al.PEG and ABA trigger methyl jasmonate accumulation to induce the MEP pathway and increase tanshinone production in Salvia miltiorrhiza hairy roots ［J］.Physiol Plantarum,2012,146（2）：173-183.doi：10.1111/j.1399-3054.2012.01603.x.

［17］ 李洁.植物转录因子与基因调控［J］.生物学通报，2004， 39（3）： 9-11.LI Jie.Plant transcription factors and gene regulation ［J］.Bull Biol,2004,39（3）：9-11.

［18］ 张鑫，宋经元，胡鸢雷，等.bHLH转录因子调控植物活性成分生物合成的研究进展［J］.药学学报，2014， 49（4）： 435-442.ZHANG Xin,SONG Jingyuan,HU Yuanlei,et al.Research progress of the regulation on active compound biosynthesis by the bHLH transcription factors in plants ［J］.Acta Pharm Sin,2014,49（4）：435-442.

［19］ WANG Jinyan,HU Zhongze,ZHAO Tongmin,et al.Genome-wide analysis of bHLH transcription factor and involvement in the infection by yellow leaf curl virus in tomato （Solanum lycopersicum）［J］.BMC Genom,2015,16（1）：1-14.doi：10.1186/s12864-015-1249-2.

［20］ PIRES N,DOLAN L.Origin and diversification of basic-helix-loop-helix proteins in plants ［J］.Mol Biol Evol,2010,27（4）：862-874.doi：10.1093/molbev/msp288.

［21］ FERRÉ-D’AMARÉ A R,PRENDERGAST G C,ZIFF E B,et al.Recognition by max of its cognate DNA through a dimeric b/HLH/Z domain ［J］.Nature,1993,363（6424）：38-45.doi：10.1038/363038a0.

［22］ XIONG Yuqing,LIU Tieyan,TIAN Chaoguang,et al.Transcription factors in rice：a genome-wide comparative analysis between monocots and eudicots ［J］.Plant Mol Biol,2005,59（1）：191-203.doi：10.1007/s11103-005-6503-6.

［23］ RUSHTON P J,BOKOWIEC M T,HAN Shengcheng,et al.Tobacco transcription factors：novel insights into transcriptional regulation in the Solanaceae ［J］.Plant Physiol,2008,147（1）：280-295.doi：https：//doi.org/10.1104/pp.107.114041.

［24］ JAILLON O,AURY J M,NOEL B,et al.The grapevine genome sequence suggests ancestral hexaploidization inmajor angiosperm phyla ［J］.Nature,2007,449（7161）：463-467.doi：10.1038/nature06148.

［25］ ZHANG Xin,LUO Hongmei,XU Zhichao,et al.Genome-wide characterisation and analysis of bHLH transcription factors related to tanshinone biosynthesis in Salvia miltiorrhiza ［J］.Sci Rep,2015,5：11244.doi：10.1038/srep11244.

［26］ 汪琬宜，蒋喜红，张利华，等.丹参转录因子 SmbHLH1基因的克隆和表达分析［J］.中国中药杂志，2011， 36（24）： 3416-3420.WANG Wanyi,JIANG Xihong,ZHANG Lihua,et al.Isolation and characteristics of SmbHLH1 gene in Salvia miltiorrhiza ［J］.China J Chin Mater Med,2011,36（24）：3416-3420.

［27］ 周宏骏，武玉翠，晋鑫鑫，等.丹参转录因子SmbHLH93的克隆及表达模式分析［J］.中草药，2014，45（23）： 3449-3455.ZHOU Hongjun,WU Yucui,JIN Xinxin,et al.Cloning and expression pattern analysis of transcription factor SmbHLH93 from Salvia miltiorrhiza ［J］.Chin Tradit Herb Drug,2014,45（23）：3449-3455.

［28］ HONG Gaojie,XUE Xueyi,MAO Yingbo,et al.Arabidopsis MYC2 interacts with DELLA proteins in regulating sesquiterpene synthase gene expression ［J］.Plant Cell,2012,24（6）：2635-2648.doi：10.1105/tpc.112.098749.

［29］ ZHOU Yangyun,SUN Wei,CHEN Junfeng,et al.SmMYC2a and SmMYC2b played similar but irreplaceable roles in regulating the biosynthesis of tanshinones and phenolic acids in Salvia miltiorrhiza ［J］.Sci Rep,2016,6：22852.doi：10.1038/srep22852.

［30］ 王晓荣.丹参bHLH类转录因子MYCs的基因克隆及功能初步研究［D］.上海：上海师范大学，2015.WANG Xiaorong.Isolation and Function Analysis of bHLH Translation Factor in Salvia miltiorrhiza ［D］.Shanghai：Shanghai Normal University,2015.

［31］ RUSHTON P J,TORRES J T,PARNISKE M,et al.Interaction of elicitor-induced DNA-binding proteins with elicitor response elements in the promoters of parsley PR1 genes ［J］.Embo J,1996,15（20）：5690-5700.

［32］ XIE Zhen,ZHANG Zhonglin,ZOU Xiaolu,et al.Annotations and functional analyses of the rice WRKY gene superfamily reveal positive and negative regulators of abscisic acid signaling in aleurone cells ［J］.Plant Physiol,2005,137（1）：176.doi：https：//doi.org/10.1104/pp.104.054312.

［33］ LI Caili,LI Dongqiao,SHAO Fenjuan,et al.Molecular cloning and expression analysis of WRKY transcription factor genes in Salvia miltiorrhiza ［J］.BMC Genom,2015,16（1）：200.doi：10.1186/s12864-015-1411-x.

［34］ 郝小龙.丹参SmWRKY3和SmWRKY70转录因子的功能研究［D］.上海：上海师范大学，2014.HAO Xiaolong.The Functional Study of SmWRKY3 and SmWRKY70 Transcription Factor in Salvia miltiorrhiza ［D］.Shanghai：Shanghai Normal University,2014.

［35］ 李磊磊.丹参WRKY转录因子SmWRKY54的功能初步研究［D］.上海：上海师范大学，2016.LI Leilei.The Functional Study of SmWRKY54 Transcription Factor in Salvia miltiorrhiza ［D］.Shanghai：Shanghai Normal University,2016.

［36］ FILICHKIN S A,PRIEST H D,GIVAN S A,et al.Genome-wide mapping of alternative splicing in Arabidopsis thaliana ［J］.Genome Res,2010,20（1）：45-58.doi：10.1101/gr.093302.109.

［37］ MARQUEZ Y,BROWN J W,SIMPSON C,et al.Transcriptome survey reveals increased complexity of the alternative splicing landscape in Arabidopsis ［J］.Genome Res,2012,22（6）：1184-1195.doi：10.1101/gr.134106.111.

［38］ SHEN Yangting,ZHOU Zhengkui,WANG Zheng,et al.Global dissection of alternative splicing in paleopolyploid soybean ［J］.Plant Cell,2014,26（3）：996-1008.doi：https：//doi.org/10.1105/tpc.114.122739.

［39］ BRADEN W,GENGKON L,GAURAV S,et al.Genome-wide landscape of alternative splicing events in brachypodium distachyon ［J］.DNA Res,2013,20（2）：163-171.doi：10.1093/dnares/dss041.

［40］ ZHANG Guojie,GUO Guangwu,HU Xueda,et al.Deep RNA sequencing at single base-pair resolution reveals high complexity of the rice transcriptome ［J］.Genome Res,2010,20（5）：646-654.doi：10.1101/gr.100677.109.

［41］ XU Zhicao,PETERS R J,WEIRATHER J,et al.Full-length transcriptome sequences and splice variants obtained by a combination of sequencing platforms applied to different root tissues of Salvia miltiorrhiza and tanshinone biosynthesis ［J］.Plant J,2015,82（6）：951-961.doi：10.1111/tpj.12865.