菠萝MYB基因AcoMYB1的克隆与表达分析

李穆 蔡元保 杨祥燕 黄思婕 李季东 谭秦亮 邱文武 方位宽

摘  要：根据菠萝转录组的测序结果克隆到1个MYB转录因子基因，命名为AcoMYB1，GenBank登录号为XM_020230319。该基因cDNA全长1221 bp，开放阅读框（ORF）为747 bp，编码一个含有248个氨基酸的蛋白。序列分析表明，AcoMYB1氨基酸序列N端具有2个保守的SANT结构域，属于R2R3类MYB转录因子。生物信息学分析表明，AcoMYB1是不稳定的亲水蛋白，不具有跨膜结构和信号肽，可能定位于细胞质，二级结构以α-螺旋和无规则卷曲为主。实时荧光定量PCR分析表明，AcoMYB1在菠萝干旱、低温逆境胁迫处理下受诱导表达，整体上表现出“先升后降”的趋势;在早熟品种和晚熟品种的果实发育过程中也被诱导表达，表现为“升-降-升”的趋势，特别是在果實发育早期和果实成熟后期受诱导表达的强度较为突出。由此推测AcoMYB1作为正调控因子参与菠萝冷害、干旱逆境胁迫的响应过程，并在菠萝果实早期发育及后期成熟发挥调控作用。

关键词：菠萝;MYB基因;基因表达;果实发育;逆境胁迫

中图分类号：S668.3      文献标识码：A

Abstract： A MYB transcription factor gene， belonging to R2R3， was cloned by the results of the transcriptome se-quencing of pineapple （Ananas comosus）， which was named AcoMYB1 and the GenBank accession number was XM_020230319. The full length cDNA and coding region （ORF） was 1221 bp and 747 bp， encoding 248 amino acids with two conserved SANT domains at N-terminal. By bioinformatics analysis， AcoMYB1 was an unstable hydrophilic protein without transmembrane structure and signal peptide， which may be located in the cytoplasm， and the secondary structure mainly had alpha helixes and random coils. Real-time fluorescence quantitative PCR analysis showed that AcoMYB1 was induced by drought and low temperature stress， and showed a trend of “up first then down”， and the expression was also induced in the fruit development of early and late maturing varieties， which showed a tendency of “up-down-up”， especially in the early stage of fruit development and the later stage of fruit maturation. Therefore， it is speculated that AcoMYB1 as a positive regulatory factor is involved in the response process of chilling injury and drought stress， and plays a regulatory role in the early fruit development and late fruit maturation of pineapple （Ananas comosus）.

Keywords： Ananas comosus; MYB genes; gene expression; fruit development; adversity stress

DOI： 10.3969/j.issn.1000-2561.2021.05.002

MYB转录因子是真核生物中最大的一类转录因子家族，一般具有DNA结合结构域、转录激活结构域和不完全界定的负调节区3个保守功能域[1]。其中，DNA结合结构域最为保守，一般由1～4个不完全重复的R基序（R1、R2和R3）组成，每个R基序约有52个氨基酸残基[2]。根据结构域特征可将MYB转录因子分为1R-MYB、R2R3-MYB、R1R2R3-MYB和4R-MYB共4个亚家族，其中R2R3-MYB转录因子是植物中最丰富的转录因子类型之一[3-4]。MYB类转录因子在植物体内有着广泛的生物学功能，几乎参与了植物各种生长发育和代谢过程，并在响应逆境胁迫中发挥着关键的调控作用[5]，这对于阐明植物抗逆分子机制及其遗传改良具有重要意义。

在MYB转录因子中，关于1R-MYB亚家族成员的功能研究相对较少，1R-MYB在调节植物生物钟、维持染色体结构完整性，以及调控DNA结合蛋白上起重要作用[6-7]。R2R3-MYB亚家族可能是由R1R2R3-MYB亚家族基因丢失了一个R1基序演化形成[8]，且是数量最多的MYB转录因子，如拟南芥（Arabidopsis thaliala）有130多个，水稻（Oryza sativa）有109个，葡萄（Vitis vinifera）有117个，大豆（Glycinemax）有244个，玉米（Zea mays）有157个等[9]。R2R3-MYB亚家族的相关功能研究在MYB转录因子中也较多，R2R3-MYB是参与新陈代谢途径最广泛的一类，主要通过调控基因表达而发挥关键作用[2]。

R2R3-MYB转录因子在植物细胞形态建成、生长发育等系列生理过程中发挥重要作用[10]，如棉花GhMYB9和GhMYB146能够调控棉花纤维的发育[11-12]，拟南芥AtMYB37、AtMYB38和AtMYB84能够调控其侧分生组织的形态建成[13]，杨树PtrMYB3和PtrMYB20参与了次生木质部的生物合成[14]，大豆GmMYB181调控其生殖器官的发育[15]。R2R3-MYB转录因子也能调控初级、次级代谢产物的生物合成，如拟南芥、葡萄等植物中多个R2R3-MYB基因调控类黄酮生物合成[16]、苯丙烷代谢途径[17-18]以及花青素的生物合成等[19-20]。此外，R2R3-MYB转录因子作为调控蛋白，能够响应生物和非生物胁迫以提高植物的抗逆性[5， 21]，如AtMYB44通过调控EIN2增强拟南芥对菜蛾和桃蚜的抗虫性[22]，野生大豆GsMYB15异源过表达能够增强转基因擬南芥对棉铃虫和盐胁迫的抗性[23]，OsMYB91通过调控SLR1的表达提高水稻的耐盐性[24]，OsMYB3R-2转基因水稻的抗旱、抗寒和耐盐性有明显提高[25]，烟草NtMYB44b能够响应MeJA、GA、ABA激素处理及干旱胁迫[26]，马铃薯StMYB44通过负调控StPHO1应对低磷胁迫[27]。

菠萝（Ananas comosus）又称凤梨，是世界四大热带水果之一，其风味独特，营养丰富，深受国内外消费者的青睐。目前，陈哲等[28]从菠萝基因组中筛选鉴定出103个R2R3-MYB转录因子，并分析其基因结构、编码蛋白和系统进化，Liu等[29]对菠萝全基因组的R2R3-MYB特征和逆境下的表达情况做了分析。而基于转录组信息分析菠萝R2R3-MYB转录因子的研究鲜为报道。因此，本研究从菠萝转录组中筛选和克隆获得1个诱导表达显著上调的R2R3-MYB基因，并通过生物信息学方法和RT-qPCR分析对其生物学特性与功能进行研究，以期为阐明菠萝R2R3-MYB基因在响应冷害、干旱逆境胁迫及其在果实发育中的分子机制提供理论基础。

1  材料与方法

1.1  植物材料和处理

试验材料为广西农业科学院广西亚热带作物研究所菠萝种质资源圃提供的‘金菠萝（MD-2）品种、‘巴厘（Comte de Paris）早熟品种和‘pérola晚熟品种。分别对苗龄为60 d的‘金菠萝组培苗进行20%浓度的PEG6000干旱胁迫和4 ℃低温胁迫处理，处理时间分别设为0、1、3、6、12、24 h，其中0 h处理的植株作为对照，每个处理设置3个重复，采集相应的叶片。在早熟品种‘巴厘谢花（从第1朵自然开花的小花算起，下同）后20、35、50、60、65、70 d分别采集果肉，在晚熟品种‘pérola谢花后20、40、60、80、100、120 d分别采集果肉。以上所有材料取样后进行液氮速冻，保存于?80 ℃超低温冰箱中，并迅速提取总RNA。

1.2  菠萝AcoMYB1基因的cDNA克隆

参照蔡元保等[30]的CTAB改良法提取菠萝叶片和果肉的总RNA，并按照Fermentas（MBI）反转录试剂盒的使用说明书进行cDNA第1链合成。利用本课题组从菠萝转录组中筛选获得1个MYB基因，根据其EST序列设计开放阅读框ORF的扩增引物AcoMYB1-F和AcoMYB1-R（表1），并获得菠萝AcoMYB1基因的最大开放阅读框ORF序列，所用引物全部由生工生物工程（上海）股份有限公司合成。

1.3  生物信息学分析

利用NCBI网站的Blast程序进行蛋白质序列的同源检索，用ORF Finder查找最大开放阅读框ORF，用SMART软件预测功能结构域，用ProtParam软件分析蛋白质的基本理化性质，用PSORT软件预测蛋白质的亚细胞定位，用SignalP软件预测信号肽，用TMpred软件预测蛋白质的跨膜区段，用SOPMA软件预测蛋白质的二级结构，用Swiss-model软件预测蛋白质的三维结构（表2）。在NCBI网站进行蛋白同源比对，并用MEGA 5.05软件构建系统进化树。

1.4  菠萝AcoMYB1基因表达分析

采用实时荧光定量RT-qPCR方法分析菠萝AcoMYB1基因的表达情况，AcoMYB1基因的特异引物为QP-F和QP-R，菠萝内参基因为18S rRNA基因，引物为Q18S-F和Q18S-R（表1）。通过RT-qPCR分析菠萝AcoMYB1基因在品种‘金菠萝幼苗期分别经干旱（20% PEG6000）和低温（4 ℃）处理0、1、3、6、12、24 h，以及早熟品种‘巴厘谢花后20、35、50、60、65、70 d和晚熟品种‘pérola谢花后20、40、60、80、100、120 d的表达情况。RT-qPCR反应程序为：94 ℃预变性3 min;94 ℃预变性30 s，58.4 ℃退火30 s，72 ℃延伸30 s;72 ℃延伸5 min，18S rRNA和AcoMYB1循环数分别为20和30。每个试验设3次重复，采用2?ΔΔCT方法分析处理数据，确定AcoMYB1基因的相对表达量。

2  结果与分析

2.1  菠萝AcoMYB1基因的克隆

以菠萝品种‘金菠萝cDNA第1链为模板，AcoMYB1-F和AcoMYB1-R为开放阅读框ORF的扩增引物。通过RT-PCR扩增和序列分析表明，菠萝AcoMYB1基因的cDNA全长序列为1221 bp，ORF长为747 bp，编码一个含248个氨基酸的蛋白质。利用NCBI网站的BLAST程序进行蛋白质序列同源检索显示，AcoMYB1蛋白含有典型的MYB结构域，属于MYB蛋白。因此，将该基因命名为AcoMYB1，GenBank登录号为XM_020230319。

2.2  AcoMYB1蛋白的多序列比对

经过多序列比对分析表明（图1），AcoMYB1蛋白与海藻、油棕、石刁柏等其他植物来源的MYB家族蛋白具有较高的相似性，即氨基酸序列相似性均在60%以上。蛋白质结构域预测表明，AcoMYB1蛋白N端含有2个保守的SANT结构域，属于MYB转录因子的R2R3亚家族蛋白，其中第1个高度保守的SANT结构域含有51个氨基酸残基（第11～61位），第2个高度保守的SANT结构域含有49个氨基酸残基（第64～112位），而非保守序列主要集中于C端。

AcoMYB1：菠萝（XP_020085908.1）;MYBAS1-like isoform X1：海藻（XP_008795797.1）;MYBAS1：

油棕（XP_029123914.1）;MYBAS2-like isoform X2：石刁柏（XP_020276648.1）;MYBAS1-like：

铁皮石斛（XP_020682892.2）;MYB59-like isoform X2：蝴蝶兰（XP_020583391.1）。

AcoMYB1： Ananas comosus （XP_020085908.1）; MYBAS1-like isoform X1： Sargassum （XP_008795797.1）; MYBAS1： Elaeis guineensis （XP_029123914.1）; MYBAS2-like isoform X2： Asparagus officinalis （XP_020276648.1）; MYBAS1-like： Dendrobium officinale （XP_020682892.2）; MYB59-like isoform X2： Phalaenopsis （XP_020583391.1）.

2.3  AcoMYB1蛋白特性及结构分析

通过ProtParam软件分析AcoMYB1蛋白的理化性质显示，AcoMYB1蛋白的相对分子量为28.25 kDa，理论等电点（pI）为5.93，不稳定系数为72.77，亲水平均数为?0.892，因此，推测AcoMYB1蛋白为不稳定亲水酸性蛋白。PSORT软件分析表明，AcoMYB1蛋白可能定位于细胞质。SignalP软件和TMpred软件预测表明，AcoMYB1蛋白不存在信号肽和跨膜结构。SOPMA软件分析结果显示（图2），AcoMYB1蛋白二级结构中主要含有α-螺旋（47.58%）和无规则卷曲（40.32%），其次是折叠延伸链（7.26%）和β-转角（4.84%），且AcoMYB1蛋白2个保守的SANT结构域（R2和R3）主要含有α-螺旋和无规则卷曲。

利用Swiss-model软件基于同源建模预测AcoMYB1蛋白的三维结构，以明确其蛋白结构域的结构特征。以WER转录因子2B的晶体结构（PDB ID：6kks.1A）为同源模板，构建AcoMYB1蛋白的三维预测模型。AcoMYB1蛋白三维结构模型的结果表明（图3），AcoMYB1蛋白的氨基酸序列与模板序列一致性为54.55%，覆盖率为44%，GMQE值为0.36，QMEAN4值为0.87，主要以单体形式存在，无配体;而且获得的R2、R3：2个保守的SANT结构域。

AcoMYB1蛋白三维模型与MYB转录因子的典型结构域高度相似，主要由α-螺旋和无规卷曲构成2个保守的SANT结构域（即DNA结合结构域R2和R3），这与AcoMYB1蛋白的二级结构预测结果相符合。

2.4  AcoMYB1蛋白的系统发育分析

为了确定菠萝AcoMYB1蛋白在MYB转录因子中的系统进化地位，将来源于不同植物共16个MYB典型蛋白构建系统进化树。分析结果表明（图4），MYB蛋白可明显分为4个进化分支，其中AcoMYB1与月季MYB59 isoform X1、苹果和罂粟MYB59-like isoform X1和石刁柏MYBAS2-like isoform X2的亲缘关系较近，具有相同的进化起源，推测他们具有类似的结构和功能;与菠萝MYBAS1-like、蝴蝶兰MYB59-like isoform X2和扁桃MYB59的亲缘关系较远。

2.5  菠萝AcoMYB1基因在不同胁迫和果实不同发育时期的表达分析

以菠萝18S rRNA为参照基因，利用实时荧光定量PCR分析菠萝AcoMYB1基因在不同胁迫和果实不同发育时期的表达情况。不同胁迫条件下的表达结果表明（图5），AcoMYB1基因在幼苗期经4 ℃低温胁迫和PEG模拟的干旱胁迫后，叶片中的表达量都有明显提高，随着时间延长诱导表达的效果越明显，整体上表现出“先升后降”的趋势，推测该基因参与菠萝冷害、干旱逆境胁迫响应。在果实不同发育时期的表达结果表明（图6），早熟品种‘巴厘和晚熟品种‘pérola果实发育成熟过程中，AcoMYB1基因均被诱导表达， 都表现为“升-降-升”的表达趋势，尤其在果实发育早期和果实成熟后期受诱导表达的强度比较突出，推测该基因在菠萝果实早期发育及后期成熟发挥调控作用。

3  讨论

本研究克隆获得了1个菠萝MYB转录因子基因，命名为AcoMYB1。氨基酸序列分析表明，AcoMYB1具有R2R3类MYB转录因子特有的2个高度保守SANT结构域（R2和R3），这与陈哲等[28]研究结果一致。多序列比对分析显示，AcoMYB1与海藻、油棕、石刁柏等其他植物来源的MYB家族蛋白相似性均在60%以上，表明R2R3类MYB转录因子在进化过程中非常保守。蛋白质理化性质分析显示，AcoMYB1蛋白为不稳定亲水酸性蛋白，其结构不稳定，易溶解。这与R2R3类MYB转录因子其他成员的研究结果一致[23-24]。系统发育分析显示，AcoMYB1与MYB59成员的亲缘关系最近，具有相同的进化起源，推测他们具有类似的结构和功能，为进一步研究AcoMYB1的功能特征奠定基礎。

R2R3-MYB轉录因子作为植物重要的调控蛋白，能够响应非生物胁迫以提高植物的抗逆性[5， 21]。在非生物胁迫下，拟南芥AtMYB60和AtMYB96能够响应干旱胁迫[31];番茄SpMYB和小麦TaODORANT1能够响应多种非生物胁迫，都能提高转基因烟草的耐盐性和耐旱性[32-33];过表达SlMYB41的转基因番茄耐旱性显著高于野生型植株[34]。这些已报道的R2R3-MYB转录因子大多作为正调控因子，其表达受到多种非生物胁迫的诱导。本研究中实时荧光定量PCR分析表明，AcoMYB1基因在幼苗期经4 ℃低温胁迫和PEG模拟的干旱胁迫后，叶片中的表达量都有明显提高。由此推测AcoMYB1作为正调控因子，参与菠萝冷害、干旱逆境胁迫的响应过程。

此外，R2R3-MYB转录因子在植物生长发育等系列生理过程中同样发挥重要作用[10]。在果实发育成熟过程中，不同品种的草莓FaMYB5基因[35]、辣椒中茄科特有MYB基因MYB31[36]、樱桃pacMYBA基因[37]等都大量表达，尤其是后期的表达量明显高于前期，推测这些基因参与果实发育成熟及其代谢产物的生物合成。本研究中实时荧光定量PCR分析表明，在早熟品种和晚熟品种果实发育成熟过程中，AcoMYB1基因都大量表达，在果实发育早期和果实成熟后期受诱导表达的强度比较突出，推测该基因在菠萝果实早期发育及后期成熟发挥调控作用。至于AcoMYB1基因是否调控菠萝果实代谢产物的生物合成还有待于进一步研究。

R2R3-MYB转录因子是植物中最丰富的转录因子类型之一，是改良植物遗传特性的良好基因资源，可为植物抗逆育种提供一种新的技术途径。本研究表明，AcoMYB1基因可能参与菠萝逆境胁迫响应以及果实发育成熟的调控，但其具体的生物学功能还需采用亚细胞定位、酵母单杂交及遗传转化等技术进一步验证。

参考文献

[1] Frampton J. MYB transcription factors： their role in growth， differentiation and disease[M]. Dordrecht： Springer Nether-lands， 2004： 6-8.

[2] Dubos C， Stracke R， Grotewold E， et al. MYB transcription factors in Arabidopsis[J]. Trends in Plant Science， 2010， 15（10）： 573-581.

[3] 陈  娜， 迟晓元， 潘丽娟， 等. MYB转录因子在植物盐胁迫调控中的研究进展[J]. 植物生理学报， 2015， 51（9）： 1395-1399.

[4] Liu J， Osbourn A， Ma P. MYB transcription factors as regu-lators of phenylpropanoid metabolism in plants[J]. Molecular Plant， 2015， 8（5）： 689-708.

[5] 牛义岭， 姜秀明， 许向阳. 植物转录因子MYB基因家族的研究进展[J]. 分子植物育种， 2016， 14（8）： 2050-2059.

[6] Yu E Y， Kim S E， Kim J H， et al. Sequence-specific DNA recognition by the MYB-like domain of plant telomeric pro-tein RTBP1[J]. Journal of Biological Chemistry， 2000， 275（31）： 24208-24214.

[7] Schaffer R， Landgraf J， Accerbi M， et al. Microarray analy-sis of diurnal and circadian-regulated genes in Arabidopsis[J]. Plant Cell， 2001， 13（1）： 113-123.

[8] Jiang C， Gu J， Chopra S， et al. Ordered origin of the typical two-and three-repeat Myb genes[J]. Gene， 2004， 326（4）： 13-22.

[9] Stracke R， Holtgr?we D， Schneider J， et al. Genome-wide identification and characterisation of R2R3-MYB genes in sugar beet （Beta vulgaris）[J]. BMC Plant Biology， 2014， 14（1）： 249-265.

[10] Ambawat S， Sharma P， Yadav N R， et al. MYB transcription factor genes as regulators for plant responses： An over-view[J]. Physiology and Molecular Biology of Plants， 2013， 19（3）： 307-321.

[11] 倪志勇， 邱迎风， 加得拉·吐留汗， 等. 棉花GhMYB146基因克隆及亚细胞定位分析[J]. 核农学报， 2016， 30（11）： 2119-2126.

[12] 邱迎风， 倪志勇， 刘  娜， 等. 陆地棉GhMYB9基因克隆及结合特性分析[J]. 核农学报， 2017， 31（3）： 440-446.

[13] Müller D， Schmitz G， Theres K. Blind homologous R2R3 MYB genes control the pattern of lateral meristem initiation in Arabidopsis[J]. The Plant Cell， 2006， 18（3）： 586-597.

[14] McCarthy R L， Zhong R， Fowler S， et al. The poplar MYB transcription factors， PtrMYB3 and PtrMYB20， are involved in the regulation of secondary wall biosynthesis[J]. Plant and Cell Physiology， 2010， 51（6）： 1084-1090.

[15] Yang H， Xue Q， Zhang Z， et al. GmMYB181， a soybean R2R3-MYB protein， increases branch number in transgenic Arabidopsis[J]. Frontiers in Plant Science， 2018， 9： 1027.

[16] 邢  文， 金晓玲. 调控植物类黄酮生物合成的MYB转录因子研究进展[J]. 分子植物育种， 2015， 13（3）： 689-696.

[17] Hichri I， Barrieu F， Bogs J， et al. Recent advances in the transcriptional regulation of the flavonoid biosynthetic pathway[J]. Journal of Experimental Botany， 2011， 62（8）： 2465-2483.

[18] Ma D， Reichelt M， Yoshida K， et al. Two R2R3-MYB pro-teins are broad repressors of flavonoid and phenylpropanoid metabolism in poplar[J]. The Plant Journal， 2018， 96（5）： 949-965.

[19] 朱婷婷， 梁  東， 夏  惠. R2R3-MYB调控果实花色苷合成的研究进展[J]. 基因组学与应用生物学， 2016， 35（4）： 985-991.

[20] Chen K， Liu H， Lou Q， et al. Ectopic expression of the grape hyacinth （Muscari armeniacum） R2R3-MYB transcription factor gene， MaAN2， induces anthocyanin accumulation in tobacco[J]. Frontiers in Plant Science， 2017， 8（7）： 965-978.

[21] Pireyre M， Burow M. Regulation of MYB and bHLH tran-scription factors： a glance at the protein level[J]. Molecular Plant， 2015， 8（3）： 378-388.

[22] 吕贝贝. MYB44对HrpNEa诱导拟南芥抗虫防卫信号传导的调控作用[D]. 南京： 南京农业大学， 2012.

[23] Shen X J， Wang Y Y， Zhang Y X， et al. Overexpression of the wild soybean R2R3-MYB transcription factor GsMYB15 enhances resistance to salt stress and Helicoverpa armigera in transgenic Arabidopsis[J]. International Journal Molecular Sciences， 2018， 19（12）： 3958.

[24] Zhu N， Cheng S， Liu X， et al. The R2R3-type MYB gene OsMYB91 has a function in coordinating plant growth and salt stress tolerance in rice[J]. Plant Science， 2015， 236： 146-156.

[25] Yang A， Dai X， Zhang W H. A R2R3-type MYB gene Os-MYB2 is involved in salt cold and dehydration tolerance in rice[J]. Journal of Experimental Botany， 2012， 63（7）： 2541-2556.

[26] 余  婧， 郭玉双， 雷  波， 等. 烟草NtMYB44b转录因子基因克隆以及生物信息学和表达分析[J]. 植物生理学报， 2020， 56（2）： 200-208.

[27] Zhou X， Zha M， Huang J， et al. StMYB44 negatively regu-lates phosphate transport by suppressing expression of PHOSPHATE1 in potato[J]. Journal of Experimental Botany， 2017， 68（5）： 1265-1281.

[28] 陈  哲， 胡福初， 阮城城， 等. 菠萝R2R3-MYB基因家族鉴定与表达分析[J]. 热带作物学报， 2019， 40（10）： 1958-1971.

[29] Liu C， Xie T， Chen C， et al. Genome-wide organization and expression profiling of the R2R3-MYB transcription factor family in pineapple （Ananas comosus）[J]. BMC Genomics， 2017， 18（1）： 503.

[30] 蔡元保， 杨祥燕， 孫光明， 等. 菠萝花发育相关基因AcMADS1的克隆与组织表达特性分析[J]. 植物学报， 2014， 49（6）： 692-703.

[31] Seo P J， Xiang F， Qiao M， et al. The MYB96 transcription factor mediates abscisic acid signaling during drought stress response in Arabidopsis[J]. Plant Physiology， 2009， 151（1）： 275-289.

[32] Li J B， Luan Y S， Yin Y L. SpMYB overexpression in tobacco plants leads to altered abiotic and biotic stress res-ponses[J]. Gene， 2014， 547（1）： 145-151.

[33] Wei Q， Luo Q， Wang R， et al. A wheat R2R3-type MYB transcription factor TaODORANT1 positively regulates drought and salt stress responses in transgenic tobacco plants[J]. Frontiers in Plant Science， 2017， 8： 1374-1389.

[34] 赵盼盼. 番茄R2R3MYB转录因子家族鉴定及SlMYB41和SlMYB64基因功能研究[D]. 泰安： 山东农业大学， 2017.

[35] Schaart J G， Dubos C， Romero D F I， et al. Identification and characterization of MYB-bHLH-WD40 regulatory complexes controlling proanthocyanidin biosynthesis in strawberry （Fragaria × ananassa） fruits[J]. New Phytologist， 2013， 197（2）： 454-467.

[36] 朱张生， 孙彬妹， 蔡  文， 等. MYB31转录因子的天然变异导致极辣辣椒的形成[J]. 辣椒杂志， 2019， 17（4）： 33-50.

[37] Shen X J， Zhao K， Liu L L， et al. A role for PacMYBA in ABA-regulated anthocyanin biosynthesis in red-colored sweet cherry cv. Hong Deng （Prunus avium L.）[J]. Plant and Cell Physiology， 2014， 55（5）： 862-880.

责任编辑：黄东杰