Figures
Abstract
R2R3-MYB proteins play role in plant development, response to biotic and abiotic stress, and regulation of primary and secondary metabolism. Little is known about the R2R3-MYB proteins in Scutellaria baicalensis which is an important Chinese medical plant. In this paper, nineteen putative SbMYB genes were identified from a S. baicalensis cDNA library, and eleven R2R3-MYBs were clustered into 5 subgroups according to phylogenetic reconstruction. In the S. baicalensis leaves which were sprayed with GA3, SbMYB2 and SbMYB7 had similar expression pattern with SbPALs, indicating that SbMYB2 and SbMYB7 might be involved in the flavonoid metabolism. Transactivation assay results showed that SbMYB2 and SbMYB7 can function as transcriptional activator. The expression of several flavonoid biosynthesis-related genes were induced or suppressed by overexpression of SbMYB2 or SbMYB7 in transgenic tobacco plants. Consistent with the change of the expression of NtDH29 and NtCHI, the contents of dicaffeoylspermidine and quercetin-3,7-O-diglucoside in SbMYB2-overexpressing or SbMYB7-overexpressing transgenic tobacco plants were decreased. The transcriptional level of NtUFGT in transgenic tobacco overexpressing SbMYB7 and the transcriptional level of NtHCT in SbMYB2-overexpressing tobacco plants were increased; however the application of GA3 inhibited the transcriptional level of these two genes. These results suggest that SbMYB2 and SbMYB7 might regulate the flavonoid biosynthesis through GA metabolism.
Citation: Yuan Y, Wu C, Liu Y, Yang J, Huang L (2013) The Scutellaria baicalensis R2R3-MYB Transcription Factors Modulates Flavonoid Biosynthesis by Regulating GA Metabolism in Transgenic Tobacco Plants. PLoS ONE 8(10): e77275. https://doi.org/10.1371/journal.pone.0077275
Editor: Ji-Hong Liu, Key Laboratory of Horticultural Plant Biology (MOE), China
Received: May 9, 2013; Accepted: August 30, 2013; Published: October 15, 2013
Copyright: © 2013 Yuan et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by the Scientific and technological research of Chinese medicine Program (06-07ZP51), Beijing nova program (2008B82) and the Self-selection program of China Academy of Traditional Chinese Medicine (ZZ2008062). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
MYB proteins present in all eukaryotes and play roles in a variety of plant-specific processes, as evidenced by their extensive functional characterization in Arabidopsis (Arabidopsis thaliana) [1], maize (Zea mays) [2], rice (Oryza sativa) [3], petunia (Petunia hybrida) [4], grapevine (Vitis vinifera L.) [5], poplar (Populus tremuloides) [6] and apple (Malus domestica) [7]. The increasing availability of plant genome sequence information has allowed comparisons and a better understanding of the evolution of this large family of transcription factors.
Most plant MYB proteins belong to the R2R3-MYB subfamily [8], and the Arabidopsis R2R3-type MYB factors encoded by the AtMYB genes have been categorized into 22 subgroups on the basis of the conserved amino acid sequence motifs [9]. Arabidopsis R2R3-MYB proteins have been found to be involved in primary and secondary metabolism, cell fate and identity, developmental processes and responses to biotic and abiotic stresses [10]. Some R2R3-MYB proteins are also involved in the regulation of the flavonoid biosynthetic pathway [11]. Overexpression of AtMYB75/PAP1 and AtMYB90/PAP2 resulted in a massive activation of phenylpropanoid biosynthetic genes and enhanced the accumulation of lignin, hydroxycinnamic acid esters, and purple anthocyanins [12]. AtMYB4 was shown to negatively regulate the expression of cinnamate 4-hydroxylase gene, then repress the synthesis of sinapoyl malate.
The roots of Scutellaria baicalensis Georgi are used to treat various diseases in Chinese traditional medicine. The active compounds of S. baicalensis include baicalin, baicalein, wogonoside, wogonin, neobaicalein, visidulin I, and oroxylin A, and these compounds exhibit anti-inflammatory, anti-tumor, and anti-HIV activities [13]. Baicalin is synthesized via the phenylpropanoid pathway by the activities of several enzymes (Figure 1), including phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), 4-coumarate:CoA ligase (4CL), chalcone synthase (CHS) and chalcone isomerase (CHI) [14]. β-glucuronidase (GUS) catalyze baicalin to baicalein [15,16]. Baicalein can be catalyzed back to baicalin by UDP-glucuronate: baicalein 7-O-glucuronosyltransferase (UBGAT) [17]. In tobacco, coumaroyl- phenylpropanoids were firstly synthesized by PAL, 4CL and C4H, and formed to caffeoyl- and feruloyl- phenylpropanoids by p-coumarate 3-hydroxylase and caffeic acid 3-O-methyltransferase [18]. CHS and CHI were also important genes in biosynthesis of anthocyanidins pathway [19]. Hydroxycinnamoyl-coenzyme A: putrescine acyltransferase (AT1) responsible for caffeoylputrescine biosynthesis in tobacco, and another acyltransferase DH29 was specific for spermidine conjugation to mediate the initial acylation step in dicaffeoylspermidine formation [20]. Hydroxycinnamoyl transferase (HCT) and caffeoyl-CoA O-methyltransferase (CCoAOMT) catalyzed the synthesis of shikimate and quinate esters in phenylpropanoid biosynthesis [21]. And both glucosyltransferase (GT) and UDP-glucose:flavonoid 3-O-glucosyltransferase (UFGT) were flavonoid-glucosyltransferases in tobacco (Figure 1) [22].
In our previous work, proteomics analysis showed that the protein level of a putative R2R3-MYB transcription factor in S. baicalensis roots was increased under water deficit condition. This R2R3-MYB has high identity with AtMYB113 which is involved in the regulation of anthocyanin biosynthesis, indicating that the R2R3-MYB transcription factor is also involved in flavonoid biosynthesis in S. baicalensis. The expression levels of several proteins related to GA metabolism were also affected. These prompted us to consider the possible linkage between flavonoid accumulation and hormone metabolism [23].
In this paper, we obtained nineteen S. baicalensis MYB transcription factors from a cDNA library, and performed a phylogeny and expression patterns analysis to yield an overview of the R2R3-MYB gene family in S. baicalensis. We further confirmed that SbMYB2 and SbMYB7 play roles in flavonoids biosynthesis.
Results
Identification of R2R3-MYB genes from S. baicalensis
We have developed a S. baicalensis full-length cDNA library (unpublished work). To identify R2R3 type MYB genes in S. baicalensis, a preliminary BLASTX search was performed using NR database in full-length cDNA library. Only hits with E values below e-30 were considered as members of this gene family. Eleven SbMYB genes have R2R3-MYB conserved domains and motifs, and their deduced proteins showed different lengths, isoelectric points, and molecular weights (Table S1; Table S2). The sequences of these nineteen SbMYB genes have been submitted to the GenBank with the accession number KC990835, KC990836, KF008651-KF008667.
Based on sequence similarity, the identified S. baicalensis R2R3-MYB proteins were clustered into 5 subgroups, according to clades with at least 50% bootstrap support (Figure 2). During the subfamily classification of the MYB genes, we also took into account the results of Stracke et al. [8] and Dubos et al. [26] for AtMYBs. The validity of our phylogenetic reconstruction is confirmed by the fact that it shows the same subgroups as those observed in previously constructed phylogenetic trees. SbMYB2, SbMYB7 and SbMYB11 belong to subgroup S14. SbMYB13 and SbMYB19 were clustered with OsMYB4 and ATMYB5, and SbMYB15 was clustered with AtMYB20, AtMYB43, AtMYB85, AtMYB42, AtMYB40 and AtMYB99. Only SbMYB8 belongs to subgroup S6, and SbMYB16 belongs to subgroup S18. In general, the gene functions of a clade appear highly but not absolutely conserved across plant species. Thus, knowledge of the gene functions of certain members will facilitate confirmation of paralogous and orthologous relationships.
The proteins are clustered into 23 subgroups, which are designated with a subgroup number (e.g., S1).
The expression pattern of S. baicalensis R2R3-MYB genes and flavonoid biosynthesis-related genes
The flavonoid accumulation in S. baicalensis might be related with GA hormone metabolism and some R2R3-MYB proteins might be involved in the flavonoid accumulation [23]. The expression of some R2R3-MYB genes and the flavonoid biosynthesis-related genes were investigated in the S. baicalensis leaves which were sprayed with GA3. The results showed that exogenous GA3 decreased the expression of Sb4CL, SbUBGAT, SbPAL1, SbPAL2 and SbPAL3, whereas the expression of SbC4H and SbCHS were increased by GA3 treatment (Figure 3). The expression of SbMYB2, SbMYB5, SbMYB7 and SbMYB12 was decreased after GA3 treatment, however GA3 treatment increased the expression of SbMYB8 (Figure 4). SbMYB2 and SbMYB7 had similar expression pattern with SbPALs, indicating that SbMYB2 and SbMYB7 might be involved in the flavonoid metabolism. The functions of these two genes were further investigated.
RT-PCR analysis of expression of SbPAL1, SbPAL2, SbPAL3, SbC4H, Sb4CL, SbCHS, SbGUS, and SbUBGAT in leaves of S. baicalensis after spraying GA3. Vertical bars indicate the standard deviation of three biological replicates. Asterisks indicate a significant difference at the P<0.05 level.
RT-PCR analysis of transcriptional level of SbMYB2, SbMYB5, SbMYB7, SbMYB8, SbMYB12, SbMYB13, and SbMYB19 in leaves of S. baicalensis at 0, 1, 2, 3 h after spraying GA3. Vertical bars indicate the standard deviation of three biological replicates. Asterisks indicate a significant difference at the P<0.05 level.
Subcellular localization of SbMYB2 and SbMYB7
Firstly, we determine the subcellular localization of SbMYB2 and SbMYB7. The full-length cDNA sequence of SbMYB2 and SbMYB7 were fused in front of the 5’ terminus of GFP reporter gene under the control of CaMV 35S promoter with the correct reading frame, respectively. The recombinant constructs of the SbMYB2-GFP and SbMYB7-GFP fusion gene and GFP alone were transformed into onion (Allium cepa) epidermal cells by particle bombardment. SbMYB7-GFP fusion protein accumulated mainly in the nucleus, suggesting that SbMYB7 is a nucleus-localized protein. Whereas SbMYB2-GFP fusion protein is located not only in nucleus but also in some other plastids，and GFP alone was present throughout the whole cell (Figure 5). These results are consistent with the predicted localization results (Table S3).
The recombinant constructs of the SbMYB2-GFP and SbMYB7-GFP fusion gene and GFP alone were transformed into onion (Allium cepa) epidermal cells by particle bombardment. A, D, G, pGEM-SbMYB2; B, E, H, pGEM-SbMYB7; C, F, I, empty vector pE3025.
Transactivation assay of SbMYB2 and SbMYB7
A yeast GAL4 system was used to determine the transcription activity of SbMYB2 and SbMYB7. The full-length cDNA of SbMYB2 and SbMYB7 was fused to the GAL4 DNA-binding domain of the pGBKT7 vector, and the fusion plasmid pBD-SbMYB2 and pBD-SbMYB7 was transformed into the yeast strain YGR2. Figure 6 showed that yeasts transformed with pBD-SbMYB2 or pBD-SbMYB7 could grow on the selection synthetic dextrose mediums lacking tryptophan and adenine (SD/-Trp/-Ade), and on the medium lacking tryptophan, adenine, and histidine (SD/-Trp/-Ade/-His). A healthy growth of yeast on both media were detected in the transformants containing the full-length cDNA of SbMYB2 and SbMYB7 compared with the control yeast transformed with empty vector, These results suggests that the SbMYB2 and SbMYB7 protein can function as transcriptional activator.
A,C, pBD-SbMYB2; B,D, pBD-SbMYB7. pGAL4 and pBD-GAL4 was used as a positive control and negative control, respectively.
Molecular characterization of transgenic tobacco lines overexpressing SbMYB2 and SbMYB7
To investigate the function of SbMYBs in plants, SbMYB2 and SbMYB7 were transformed into tobacco plants, respectively. The integration of SbMYB2 and SbMYB7 was confirmed using PCR analysis (Figure S1). The real-time RT-PCR analysis results showed that the expression of SbMYB2 and SbMYB7 was markedly increased in the transgenic plants (Table S4). Three independent transgenic lines (e10-25, e10-26 and e10-29) overexpressing SbMYB2 and three independent transgenic lines (e18-53, e18-b and e18-d) overexpressing SbMYB7 were selected for further analysis.
SbMYB2 and SbMYB7 regulates the expression of flavonoid biosynthesis-related genes
To investigate whether SbMYB2 and SbMYB7 regulates the flavonoid biosynthesis, the expression of several flavonoid biosynthesis-related genes including NtPAL1, NtPAL2, NtC4H, NtCHS, NtCHI, NtUFGT, NtGT4, NtAT1, NtDH29, NtHCT and NtCCoAMT1 in transgenic tobacco plants were measured using real-time RT-PCR with the specific primers. The expression of NtPAL1, NtPAL2, NtC4H and NtUFGT were markedly induced in transgenic tobacco plants expressing SbMYB2 or SbMYB7, indicating that SbMYB2 or SbMYB7 positively regulates the expression of these genes. The expression of NtCHI was decreased in transgenic plants overexpressing SbMYB2 or SbMYB7 and overexpression of SbMYB2 decreased the expression of NtGT4 (Figure 7), indicating that SbMYB2 negatively regulate the expression of CHI and GT4. The expression of NtAT1 was increased and the transcriptional level of NtDH29 was slightly decreased in both SbMYB2-overexpressing and SbMYB7-overexpressing transgenic tobacco plants (Figure 8). Because NtDH29 is an important gene which is involved in the biosynthesis of dicaffeoylspermidine [20], above results indicates that SbMYB2 and SbMYB7 could be related to the dicaffeoylspermidine formation. The expression of flavonoid related genes in wild type tobacco plants were not affected by GA treatment (Table S5).
Quantitative real-time PCR analysis of transcription level of NtPAL1, NtPAL2, NtC4H, NtCHS, NtCHI, NtGT4 and NtUFGT in leaves of SbMYB2-overexpressing and SbMYB7-overexpressing tobacco plants.
Quantitative real-time PCR analysis of transcript level of NtPAL1, NtPAL2, NtC4H, NtCHS, NtCHI, NtGT4, NtUFGT, NtAT1, NtCCoAMT1, NtDH29 and NtHCT in leaves of SbMYB2-overexpressing and SbMYB7-overexpressing tobacco plants at 0, 1, 2, 3 h after spraying GA3, respectively. Vertical bars indicate the standard deviation of three biological replicates. Asterisks indicate a significant difference at the P<0.05 level.
SbMYB2 and SbMYB7 negatively regulated the synthesis of dicaffeoylspermidine and flavone in transgenic tobacco
The overexpression of SbMYB2 or SbMYB7 led to the change of the expression of flavonoid biosynthesis-related genes, indicating that SbMYB2 or SbMYB7 might regulate the accumulation of flavonoid. To investigate whether the overexpression of SbMYB2 or SbSMB7 in transgenic tobacco plants affected the accumulation of flavonoid, we performed HPLC analysis of the leaf samples. Two peaks having significant areas were inhibited in both transgenic plants overexpressing SbMYB2 or SbMYB7, compared with wild type plants, and were identified as dicaffeoylspermidine and quercetin-3,7-O-diglucoside by LC-MS analysis (Figure 9). Overexpression of SbMYB2 or SbMYB7 decreased the accumulation of dicaffeoylspermidine and quercetin-3,7-O-diglucoside (Figure 9D), indicating that SbMYB2 and SbMYB7 negatively regulated flavonoid synthesis in transgenic tobacco plants.
A, UPLC analysis. B,C Dicaffeoylspermidine and quercetin-3,7-O-diglucoside were identified using LC-MS and LC-MS/MS. D, content of dicaffeoylspermidine and quercetin-3,7-O-diglucoside in transgenic tobacco.
Exogenous GA3 affected the expression of flavonoid pathway genes in transgenic tobacco plants over-expressing SbMYB2 or SbMYB7
To further analyze the function of SbMYB2 and SbMYB7 in possible linkage among flavonoid accumulation and GA metabolism, exogenous GA3 were sprayed on the leaves of transgenic plants overexpressing SbMYB2 or SbMYB7, and the expression of several flavonoid pathway genes were measured by real-time RT-PCR. Expression of NtPAL1 was decreased in both SbMYB2-overexpressing and SbMYB7-overexpressing transgenic plants at 2 and 3h after spay exogenous GA3. The expression of NtPAL2 was markedly increased in SbMYB2-overexpressing transgenic plants at 1h after spaying exogenous GA3, whereas was increased in SbMYB7-overexpressing transgenic plants at 2h after exogenous GA3 treatment. Transcriptional level of NtCHS was increased in SbMYB2 and SbMYB7-overexpressing transgenic plants at 3h after spaying exogenous GA3. The transcriptional levels of NtCHI and NtUFGT were decreased at 1h and increased at 2 and 3h in SbMYB2-overexpressing transgenic plants after spaying exogenous GA3, and decreased at 1h in SbMYB7-overexpressing transgenic plants after spaying exogenous GA3. The expression of NtHCT was decreased in SbMYB2-overexpressing transgenic plants at 3h after spaying exogenous GA3 (Figure 8). Gene expression pattern without GA application was also investigated and no difference was observed (Table S6).
SbMYBs bind with the box L sequence of the PAL promoter
Several reports have shown that MYB proteins regulate the expression of PAL by combining the box L [24], and box L is present in the NtPAL gene [25]. The box L sequence in the promoter sequence of NtPAL (GenBank:AB008199) was predicted as ACTTTG using Softberry (linux1.softberry.com). The sequence contains ACTTTG, which has been identified in several gene promoters as a component of binding sites for transcription factor. The interaction between SbMYB2 and SbMYB7 with NtPAL promoter sequence was assayed with electrophoretic mobility shift assay (EMSA) experiment. SbMYB2 and SbMYB7 were expressed in E. coli, respectively, for the use of EMSA analysis. No binding bands were detected with crude proteins of E. coli without or with empty vector (Figure 10 lane 1 and lane 2). SbMYB2 and SbMBY7 specifically bind with the box L sequence, and unlabeled probes inhibit the binding (Figure 10). These results confirmed that SbMYB proteins could combine to the box L sequence of NtPAL gene which is the target gene of MYB protein.
1, E.coli; 2, pGEX-4T-1; 3,4, SbMYB2; 5,6, SbMYB7; 3,5, biotin labeled box L probe; 4,6, biotin labeled and unlabeled box L probes.
Discussion
MYB proteins are key factors in regulatory networks controlling development, metabolism and responses to biotic and abiotic stresses. The duplication of R2R3-MYBs plays a key role in generating diversity of gene function [26]. AtMYB75/PAP1, AtMYB90/PAP2, AtMYB113 and AtMYB114 (subgroup 6) control anthocyanin biosynthesis in vegetative tissues [27]. Chalcone synthese (CHS), chalcone isomerase (CHI), flavonol 3-hydroxylase (F3H), flavonol 3’-hydroxylase (F3’H) and flavonol synthase (FLS) are positively regulated by subgroup 7 R2R3-MYBs, and dihydroflavonol-4-reductase (DFR), leucoanthocyanidin dioxygenase (LDOX) and anthocyanidin reductase (ANR) are activated by subgroup 6 R2R3-MYBs [28]. When ectopically expressed in PA-accumulating cells, all R2R3-MYB members in subgroup 4 were able to inhibit PA biosynthesis, suggesting that these proteins have the ability to inhibit the expression of DFR, LDOX and ANR [26].
In this study, we isolated nineteen full-length SbMYB genes from a S. baicalensis cDNA library, and eleven R2R3-MYBs with conserved R2R3 domain were divided into 5 subgroups based on the conservation of the DNA binding domain of Arabidopsis MYB proteins. Analysis of new plant genomes suggests that some MYB genes have evolved to fulfill lineage-specific functions [29]. Despite the divergence of the amino-acid sequence outside of the MYB domain, there are some conserved motifs that may contribute to function. Within subgroups that are conserved between divergent species, primary protein structures and biological functions are correlated, such as phenylpropanoid metabolism regulation by R2R3-MYB subgroups 6 [9]. Therefore, protein structure and gene expression patterns can help deduce the functions of new MYB proteins in plants. Using this way, some Arabidopsis MYB proteins were predicted to have the function of controlling flavonoids [30] and flavonol [31] biosynthesis. In this study, the transcriptional levels of SbMYB2 and SbMYB7 which belong to subgroup 14 were decreased by GA3 treatment, whereas the expression of SbMYB8 which belongs to subgroups 6 was increased after spraying exogenous GA. The transcriptional levels of SbCHS and SbC4H, two key genes which are involved in baicalein biosynthesis, were also increased. SbCHS and SbC4H have the similar expression pattern with SbMYB8, indicating that SbMYB8 was involved in the flavonoid biosynthesis in S. baicalensis based on subgroup classification and co-expressed analysis.
MYB proteins in subgroup 14 were believed to have functions on plant development [32]. Here, SbMYB2 and SbMYB7 which belong to subgroup 14 co-expressed with SbPALs, indicating that SbMYB2 and SbMYB7 might be involved in the flavonoid metabolism. To confirm this hypothesis, transgenic tobacco plants overexpressing SbMYB2 or SbMYB7 were developed. In transgenic plants, the transcriptional level of some flavonoid biosynthesis-related genes (NtPAL1, NtPAL2, NtC4H and NtUFGT) were increased, whereas the transcription levels of NtCHI and NtGT4 were decreased, suggesting that SbMYB2 and SbMYB7 could up-regulate the first step and down-regulate the last step of flavonoid biosynthesis. Consistent with the decreased expression of NtDH29 and NtCHI, the content of dicaffeoylspermidine and quercetin-3,7-O-diglucoside in transgenic tobacco plants was significantly decreased by overexpression of SbMYB2 or SbMYB7, suggesting phenylpropanoid-polyamine conjugates was negatively regulated by above SbMYBs. Dicaffeoylspermidine was a phenylpropanoid-polyamine conjugates and it has been shown that the N-coupling reaction of polyamines to phenolic acids (such as cinnamic, p-coumaric, caffeic, ferulic and sinapic acids) in plants is catalyzed by a specific class of acyltransferase enzymes. NaMYB8 silencing induces specific alterations in the accumulation of coumaroyl-containing metabolites and suppresses caffeoyl- and feruloyl- containing metabolites, and resulted in a strong suppression of dicaffeoylspermidine in Nicotiana [33].
Plant hormones affect the accumulation of secondary metabolites and R2R3-type MYB proteins also participate in mediating hormone actions [34]. It has been observed that ABA and GA3 treatment decreased CsPAL expression level and catechin content [35]. GA3 may inhibit the phenylpropanoid pathway through affecting PAL in Myrica rubra, pea and carrot [36-38]. Devaiah et al. [39] reported that AtMYB62 regulated phosphate starvation responses via changes in GA metabolism and signaling. Gibberellin acts through jasmonate to control the expression of AtMYB21, AtMYB24, and AtMYB57 to promote stamen filament growth in Arabidopsis [40]. Rice GaMYB is an important component of GA signaling in cereal aleurone cells and anther development [41]. In our previous study, the levels of total flavonoids and baicalin and the ratio of baicalin to baicalein in roots of S. baicalensis were decreased under water deficit condition after application of GAs, and these decreases were recovered after application of paclobutrazol [23]. The results in this paper also insist that GAs affected flavonoid metabolism in S. baicalensis. Over-expression of Arabidopsis MYR1 or MYR2 produced GA-deficient symptoms that were rescued by application of GA3 [42]. Our result also showed that the transcription level of NtUFGT in SbMYB7-overexpressing tobacco and NtHCT transcripts in SbMYB2-overexpressing tobacco was increased. These changes could be rescued by application of GA3. These results suggest that SbMYB2 and SbMYB7 might regulate the flavonoid biosynthesis through the negative effect on levels of bioactive GA.
ZmC1 regulates anthocyanin production together with ZmR in maize, suggesting that R2R3 MYB proteins are often involved in the combinatorial interaction of transcription factors for the generation of highly specific expression patterns [43]. The transcriptional levels of NtPAL1 and NtPAL2 were increased in transgenic tobacco plants overexpressing SbMYB2 or SbMYB7. NtPAL2 transcription level was also increased in transgenic plants overexpressing SbMYB2 and SbMYB7 after GA3 treatment. It has been reported that R2R3-MYB factors regulate the transcriptional activation of Pinus pinaster PAL by interaction with the promoter sequence containing AC elements [44], and MYB proteins regulate the expression of PAL by combining with the box L [24]. EMSA analysis clearly showed that SbMYB2 and SbMYB7 could combine to the box L sequence of NtPAL gene which is the target gene of MYB protein (Figure 10).
Within a subgroup, paralogs can control the same metabolic pathway in different cell types as a result of differences in expression patterns [27,30]. AtMYB66/WER and GL1, both clustering together in subgroup 15, can functionally complement each other and display different biological functions only because of their different spatial expression patterns [45]. Although both SbMYB2 and SbMYB7 belong to subgroup 14, they have 62.1% identity at nucleotide acid level and 45.8% identity at amino acid level, and they have different subcellular localization. The transcriptional level of flavonoid biosynthesis-related genes in SbMYB2- and SbMYB7-overexpressing transgenic tobacco after spray GA3 have different expression patterns, indicate that SbMYBs play redundant, but divergent roles in flavonoid biosynthesis and GA response. The results in this study suggested that SbMYB2 and SbMYB7 could affect phenylpropanoid biosynthesis, and SbMYB2 affects flavonoid accumulation through regulating the gibberellin (GA) signaling pathways.
Materials and Methods
Plant Materials and Growth Condition
The seeds of S. baicalensis were obtained from Institute of Chinese Materia Medica, Academy of Chinese Medical Sciences, Beijing, China), sterilized in 0.5% NaOCl for 5 min, then washed 3 times with sterile water, and placed in petri dishes to germinate. The seedlings 2 weeks after germination were transferred to individual pots (10 seedlings per pot) containing 500 g dried soil in climate chamber at 25°C with 16 h-light photoperiod under well-water condition. GA3 (100 uM) were sprayed on leaves of plants one week after transplant of S. baicalensis and transgenic tobacco. The leaves were sampled three times at 1, 2 and 3 h after spraying. The sample were rinsed three times in distilled water, and then stored at -80°C for further experiments
Identification of R2R3-MYB Protein in S. baicalensis
To identify R2R3-MYB genes, we performed a BLASTX algorithm [46] at the S. baicalensis full-length cDNA library (Yuan et al, unpublished) against the amino acid sequences in NR database (http://www.ncbi.nlm.nih.gov). The functional and structural domains were predicted by InterProScan [47] and Blast2GO [48] analysis, respectively.
Sequence Analysis
The deduced amino acid sequences were adjusted manually using BioEdit (version 7.0.0) with the default parameters (Pittsburgh Supercomputing Center; http:// www.psc.edu/biomed/genedoc/). Open read frame of 19 R2R3-MYB proteins were performed by BioEdit [49]. Theoretical isoelectric points and molecular weights were predicted using the Compute pI/MW tool on the ExPASy server (http://web.expasy.org/compute_pi/) [50]. The localizations of the deduced proteins were predicted on the ProtComp Version 9.0 (http://linux1.softberry.com/berry.phtml?topic=protcompan&group=programs &subgroup=proloc) as well as SubLoc v1.0 (http://www.bioinfo.tsinghua.edu.cn/SubLoc/). The conserved amino acids were analyzed by protein alignment using such tools as ClustalW and checked manually [51].
Construction of the Phylogenetic Trees
Phylogenetic analysis of the alignments was performed using ClustalW (Thompson, 1994) and MEGA 4.0 [52] for neighbor-joining analysis. The reliability of these tree topologies was evaluated using bootstrap support with 1000 replicates [53]. The sequences of 126 Arabidopsis R2R3-MYB proteins were downloaded from the TAIR Arabidopsis genome (http://www.arabidopsis.org/). The predicted proteins of 52 well-known plant R2R3-MYB genes were collected from the National Center for Biotechnology Information (NCBI, http://www.ncbi.nlm.nih.gov/).
Gene expression analysis in S. baicalensis
Total RNA was extracted from plant tissues using Trizol reagent (Invitrogen, USA) and pretreated with RNase-Free DNase (Promega, USA) to eliminate genomic DNA contamination. RNA integrity was analyzed on 1% agarose gel. RNA quantity was determined using a NanoDrop 2000C spectrophotometer (Thermo Scientific, USA).
Semiquantitative RT-PCR was carried out for SbPAL1 (HM062775), SbPAL2 (HM062776), SbPAL3 (HM062777), SBC4H (HM062778), Sb4CL (HM166767), SbCHS (AB008748), SbUBGAT (EF512580), SbGUS (AB040072), SbMYBs and Sb18S (FJ527609) using the One-Step RT-PCR kit (TakaRa) with specific primers (Table S7). The Sb18S gene was chosen as a loading control. The one-step RT-PCR was done as follows: 94°C for 3 min, 31 cycles of 94°C for 30 s, annealing temperature for 40 s, and 72°C for 40 s, and 72°C for 10 min.
Subcellular localization
The whole coding sequence of SbMYB2 and SbMYB7 was ligated into pE3025 vector [54] digested with EcoRI and KpnI to generate plasmids pGEM-SbMYB2 and pGEM-SbMYB7，respectively. In both plasmids, SbMYB-GFP fusion genes are under the control of CaMV 35S promoter. The construct was confirmed by sequencing and used for transient transformation of onion (Allium cepa) epidermis via a gene gun (Bio-Rad). After 24 h of incubation, GFP fluorescence in transformed onion cells was observed under a confocal microscope (Zeiss, Germany).
Transactivation assay
To determine the transactivation activity, the open reading frames of SbMYB2 and SbMYB7 were generated by PCR amplification, cloned into vector pBD-GAL4 which was digested with EcoRI and SalI, to construct pBD-SbMYB2, and pBD-SbMYB7, respectively. The constructs were transformed into YGR2 cells by the lithium acetate-mediated method. The transformed yeast strains were placed on SD/–Trp medium at 28 °C for 2 days. Yeast transformants from SD medium lacking Trp were then transferred and streaked onto solid SD agar lacking Trp/His/Ade (SD/–Trp/–His/–Ade) to score the growth response after 3 days. For the colony-lift filter assay (X-gal assay), the yeast was transferred to Whatman filter paper plus X-gal for transcription activation activity analysis within 8 h. pGAL4 and pBD-GAL4 was used as a positive control and negative control, respectively.
Tobacco transformation
SbMYB2 and SbMYB7 fragments were inserted into binary vectors, pCambia1305 to produce p35Spro-SbMYB2 and p35Spro-SbMYB7, respectively. The constructs were then transformed into Agrobacterium tumefaciens EHA105. Tobacco (Nicotiana tabacum) leaf discs were transformed via an A. tumefaciens mediated leaf disc procedure [55] and selected using 50 mg L-1 Hygromycin B and 200 mg L-1 carbenicillin. After rooting and acclimatization, regenerated plants were grown in a greenhouse to set seeds by self-pollination. T1 transgenic plants were used for further analyses.
Chemical analysis
To determine flavonoid content, 100 mg powdered tobacco leaf was extracted for 1 h in 1 mL ethyl alcohol. The solution was filtered through a membrane filter (0.2 μm), and flavonoid concentrations were determined using an UPLC-Q-Tof system with a 1.0 mL/min flow rate. UPLC was performed on a diamonsil C18 column (4.6 mm×250 mm, 5 μm). The detection wavelength was set at 354 nm and the column temperature was maintained at 30°C. The mobile phase consisted of acetonitrile-methanoic acid (A; 99.9:0.1, v/v) and deionized water-trifluoroacetic acid (B; 99.9:0.1, v/v). The initial condition was A–B (5:95, v/v) for 20 min, and this was linearly changed to A–B (10:90, v/v) at 20 min, to A–B (20:80, v/v) at 40 min, and to A–B (40:60, v/v) at 60 min. UPLC grade acetonitrile (E. Merck, Darmstadt, Germany) was used for the UPLC analysis. Dicaffeoylspermidine and quercetin-3,7-O-diglucoside were identified using LC-MS and LC-MS/MS. The injection volume of the sample solution was 20 μl, and the experiment was repeated six times.
Quantitative real-time PCR
Total RNA was reverse-transcribed using Reverse Transcriptase MMLV (Takara, China). PCRs were performed using SYBR Premix Ex Taq kits (TaKaRa, China) following the manufacturer’s instructions and conducted in triplicate using an ABI 7500 Real-Time PCR System (ABI, USA). Gene-specific primers of NtPAL1 (M84466), NtPAL2 (D17467), NtC4H (AJ937847), NtCHI (AB213651), NtCHS (AF311783), NtUFGT (GQ395697), NtGT4 (AB176522), NtAT1 (JN390826), NtDH29 (JN390824), NtCCoAMT (NTU62736), NtHCT (NTU62736), SbMYB2 and SbYB7 were designed using Primer3 (http://frodo.wi.mit.edu/primer3/). The primer sequences are listed in Table S7. The lengths of PCR products ranged from 100 to 250 bp. Ntactin was chosen as an endogenous control in studying gene expressions in various samples of transgenic tobacco. The specificity of amplification was assessed by melting curve analysis, and the relative abundance of genes was determined using the comparative Ct method as suggested in ABI 7500 Software v2.0.1 (ABI).
Expression of SbMYBs protein in E.coli
The open reading frame (ORF) of SbMYB2 and SbMYB7 was cloned into the expression vector pGEX-4T-1 and transformed into Transetta (DE3) chemically competent cells (Beijing TransGen Biotech Co., Ltd, China), respectively. The vector pGEX-4T-1 (+) allows in-frame cloning of PCR products resulting in a GST-tag attached at the N-terminal end of the recombinant protein. Expression of the recombinant protein was induced by adding isopropyl-β-D-1-thiogalactopyranoside (IPTG) and cells were harvested at 9h.
Electrophoretic Mobility Shift Assay
The boxes L sequence in the promoter sequence of NtPAL (GenBank:AB008199) was as ACTTTG using Softberry (linux1.softberry.com). Oligonucleotides of boxes L sequence were synthesized and labeled with biotin (Sangon Biotech (Shanghai) Co., Ltd., China) for chemiluminescence using a lightshift chemiluminescent electrophoretic mobility shift assay kit (Pierce). After labeling, complementary labeled strands were mixed together in an equimolar ratio and annealed at room temperature after denaturation at 90°C. Gel mobility shift assays were performed by incubating 0.5 ng of labeled probe with SbMYBs protein and competing oligonucleotides in binding buffer (10 mM Tris-HCl, pH 7.5, 50 mM NaCl, 1 mm dithiothreitol, 1 mm EDTA, 5% glycerol, and 1 μg/μl poly(dIdC) at room temperature for 30 min. Mixtures were size-fractionated on a non-denaturing 46% polyacrylamide gel followed by drying and transfer to nitrocellulose membranes and detection by streptavidin-HRP/chemiluminescence for biotin-labeled probes.
Supporting Information
Figure S1.
PCR analysis of transgenic tobacco. M, 2000bp DNA ladder; CK+, gene; CK-, wild-type tobacco.
https://doi.org/10.1371/journal.pone.0077275.s001
(TIF)
Table S1.
Blast results of MYB transcription factors in Scutellaria baicalensis.
https://doi.org/10.1371/journal.pone.0077275.s002
(DOC)
Table S2.
Characteristics of R2R3-MYB proteins in Scutellaria baicalensis.
https://doi.org/10.1371/journal.pone.0077275.s003
(DOC)
Table S3.
Subcellular localization predicted by ProtComp and SubLoc.
https://doi.org/10.1371/journal.pone.0077275.s004
(DOC)
Table S4.
Real-time RT-PCR analysis of transgenic tobacco.
https://doi.org/10.1371/journal.pone.0077275.s005
(DOC)
Table S5.
Transcriptional level of flavonoid biosynthesis genes in wild-type tobacco.
https://doi.org/10.1371/journal.pone.0077275.s006
(DOC)
Table S6.
Expression patterns of the genes in transgenic tobacco without GA application.
https://doi.org/10.1371/journal.pone.0077275.s007
(DOC)
Author Contributions
Conceived and designed the experiments: YJL LQH. Performed the experiments: YY CW JY. Analyzed the data: YY CW. Contributed reagents/materials/analysis tools: YY. Wrote the manuscript: YY CW YJL.
References
- 1. Shinozaki K, Yamaguchi-Shinozaki K, Urao T, Koizumi M (1992) Nucleotide sequence of a gene from Arabidopsis thaliana encoding a myb homologue. Plant Mol Biol 19: 493-499. doi:https://doi.org/10.1007/BF00023398. PubMed: 1623193.
- 2. Du H, Feng BR, Yang SS, Huang YB, Tang YX (2012) The R2R3-MYB transcription factor gene family in maize. PLOS ONE 7: e37463. doi:https://doi.org/10.1371/journal.pone.0037463. PubMed: 22719841.
- 3. El-Kereamy A, Bi YM, Ranathunge K, Beatty PH, Good AG et al. (2012) The rice R2R3-MYB transcription factor OsMYB55 is involved in the tolerance to high temperature and modulates amino acid metabolism. PLOS ONE 7: e52030. doi:https://doi.org/10.1371/journal.pone.0052030. PubMed: 23251677.
- 4. Spitzer-Rimon B, Farhi M, Albo B, Cna'ani A, Ben Zvi MM et al. (2012) The R2R3-MYB-Like Regulatory Factor EOBI, Acting downstream of EOBII, regulates scent production by activating ODO1 and structural scent-related genes in Petunia. Plant Cell 24: 5089-5105. doi:https://doi.org/10.1105/tpc.112.105247. PubMed: 23275577.
- 5. Czemmel S, Heppel SC, Bogs J (2012) R2R3 MYB transcription factors: key regulators of the flavonoid biosynthetic pathway in grapevine. Protoplasma 249 S2: 109-118. doi:https://doi.org/10.1007/s00709-012-0380-z. PubMed: 22307206.
- 6. Legay S, Sivadon P, Blervacq AS, Pavy N, Baghdady A et al. (2010) EgMYB1, an R2R3 MYB transcription factor from eucalyptus negatively regulates secondary cell wall formation in Arabidopsis and poplar. New Phytol 188: 774-786. doi:https://doi.org/10.1111/j.1469-8137.2010.03432.x. PubMed: 20955415.
- 7. Chagné D, Lin-Wang K, Espley RV, Volz RK, How NM et al. (2013) An ancient duplication of apple MYB transcription factors is responsible for novel red fruit-flesh phenotypes. Plant Physiol 161: 225-239. doi:https://doi.org/10.1104/pp.112.206771. PubMed: 23096157.
- 8. Stracke R, Werber M, Weisshaar B (2001) The R2R3-MYB gene family in Arabidopsis thaliana. Curr Opin Plant Biol 4: 447-456. doi:https://doi.org/10.1016/S1369-5266(00)00199-0. PubMed: 11597504.
- 9. Kranz HD, Denekamp M, Greco R, Jin H, Leyva A et al. (1998) Towards functional characterisation of the members of the R2R3-MYB gene family from Arabidopsis thaliana. Plant J 16: 263-276. doi:https://doi.org/10.1046/j.1365-313x.1998.00278.x. PubMed: 9839469.
- 10. Meissner RC, Jin H, Cominelli E, Denekamp M, Fuertes A et al. (1999) Function search in a large transcription factor gene family in Arabidopsis: assessing the potential of reverse genetics to identify insertional mutations in R2R3 MYB genes. Plant Cell 11: 1827-1840. doi:https://doi.org/10.2307/3871080. PubMed: 10521515.
- 11. Dubos C, Le Gourrierec J, Baudry A, Huep G, Lanet E et al. (2008) MYBL2 is a new regulator of flavonoid biosynthesis in Arabidopsis thaliana. Plant J 55: 940-953. doi:https://doi.org/10.1111/j.1365-313X.2008.03564.x. PubMed: 18532978.
- 12. Borevitz JO, Xia YJ, Blount J, Dixon RA, Lamb C et al. (2000) Activation tagging identifies a conserved MYB regulator of phenylpropanoid biosynthesis. Plant Cell 12: 2383-2393. doi:https://doi.org/10.2307/3871236. PubMed: 11148285.
- 13. Błach-Olszewska Z, Jatczak B, Rak A, Lorenc M, Gulanowski B et al. (2008) Production of cytokines and stimulation of resistance to viral infection in human leukocytes by Scutellaria baicalensis flavones. J Interf Cytok Res 28: 571-581. doi:https://doi.org/10.1089/jir.2008.0125. PubMed: 18771341.
- 14. Xu H, Park NI, Li X, Kim YK, Lee SY et al. (2010) Molecular cloning and characterization of phenylalanine ammonia-lyase, cinnamate 4-hydroxylase and genes involved in flavone biosynthesis in Scutellaria baicalensis. Bioresour Technol 101: 9715-9722. doi:https://doi.org/10.1016/j.biortech.2010.07.083. PubMed: 20708926.
- 15. Morimoto TN, Matsuda T, Tanaka H, Taura F, Furuya N et al. (1998) Novel hydrogen peroxide metabolism in suspension cells of Scutellaria baicalensis Georgi. J Biol Chem 273: 12606-12611. doi:https://doi.org/10.1074/jbc.273.20.12606. PubMed: 9575222.
- 16. Sasaki K, Taura F, Shoyama Y, Morimoto S (2000) Molecular characterization of a novel b-Glucuronidase from Scutellaria baicalensis Georgi. J Biol Chem 275: 27466-27472. PubMed: 10858442.
- 17. Nagashima S, Hirotani M, Yoshikawa T (2000) Purification and characterization of UDP-glucuronate: baicalein 7-O-glucuronosyltransferase from Scutellaria baicalensis Georgi. cell suspension cultures. Phytochem 53: 533-538. doi:https://doi.org/10.1016/S0031-9422(99)00593-2. PubMed: 10724177.
- 18. Cook CM, Daudi A, Millar DJ, Bindschedler LV, Khan S et al. (2012) Transcriptional changes related to secondary wall formation in xylem of transgenic lines of tobacco altered for lignin or xylan content which show improved saccharification. Phytochemistry 74: 79-89. doi:https://doi.org/10.1016/j.phytochem.2011.10.009. PubMed: 22119077.
- 19. Kumar V, Yadav SK (2013) Overexpression of CsANR increased flavan-3-ols and decreased anthocyanins in transgenic Tobacco. Mol Biotechnol, 54: 426–35. doi:https://doi.org/10.1007/s12033-012-9580-1. PubMed: 22872496.
- 20. Onkokesung N, Gaquerel E, Kotkar H, Kaur H, Baldwin IT et al. (2012) MYB8 controls inducible phenolamide levels by activating three novel hydroxycinnamoyl-coenzyme A: polyamine transferases in Nicotiana attenuate. Plant Physiol 158: 389-407. doi:https://doi.org/10.1104/pp.111.187229. PubMed: 22082505.
- 21. Hoffmann L, Maury S, Martz F, Geoffroy P, Legrand M (2003) Purification, cloning and properties of an acyltransferase controlling shikimate and quinate ester intermediates in phenylpropanoid metabolism. J Biol Chem 278: 95-103. PubMed: 12381722.
- 22. Senthil-Kumar M, Hema R, Suryachandra TR, Ramegowda HV, Gopalakrishna R et al. (2010) Functional characterization of three water deficit stress-induced genes in tobacco and Arabidopsis: an approach based on gene down regulation. Plant Physiol Biochem 48: 35-44. doi:https://doi.org/10.1016/j.plaphy.2009.09.005. PubMed: 19811926.
- 23. Yuan Y, Liu Y, Wu C, Chen S, Wang Z et al. (2012) Water deficit affected flavonoid accumulation by regulating hormone metabolism in Scutellaria baicalensis Georgi roots. PLOS ONE 7: e42946. doi:https://doi.org/10.1371/journal.pone.0042946. PubMed: 23077481.
- 24. Maeda K, Kimura S, Demura T, Takeda J, Ozeki Y (2005) DcMYB1 acts as a transcriptional activator of the carrot phenylalanine ammonia-lyase gene (DcPAL1) in response to elicitor treatment, UV-B irradiation and the dilution effect. Plant Mol Biol 59: 739-752. doi:https://doi.org/10.1007/s11103-005-0910-6. PubMed: 16270227.
- 25. Taguchi G, Sharan M, Gonda K, Yanagisawa K, Shlmosaka M et al. (1998) Effect of methyl jasmonate and elicitor on PAL gene expression in tobacco cultured cells. J Plant Biochem Biotechnol 7: 79-84. doi:https://doi.org/10.1007/BF03263040.
- 26. Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C et al. (2010) MYB transcription factors in Arabidopsis. Trends Plant Sci 15: 573-581. doi:https://doi.org/10.1016/j.tplants.2010.06.005. PubMed: 20674465.
- 27. Gonzalez A, Zhao M, Leavitt JM, Lloyd AM (2008) Regulation of the anthocyanin biosynthetic pathway by the TTG1/bHLH/Myb transcriptional complex in Arabidopsis seedlings. Plant J 53: 814-827. doi:https://doi.org/10.1111/j.1365-313X.2007.03373.x. PubMed: 18036197.
- 28. Preston J, Wheeler J, Heazlewood J, Li SF, Parish RW (2004) AtMYB32 is required for normal pollen development in Arabidopsis thaliana. Plant J 40: 979-995. doi:https://doi.org/10.1111/j.1365-313X.2004.02280.x. PubMed: 15584962.
- 29. Bailey PC, Dicks J, Wang TL, Martin C (2008) IT3F: a web-based tool for functional analysis of transcription factors in plants. Phytochemistry 69: 2417-2425. doi:https://doi.org/10.1016/j.phytochem.2008.08.003. PubMed: 18771782.
- 30. Lepiniec L, Debeaujon I, Routaboul JM, Baudry A, Pourcel L et al. (2006) Genetics and biochemistry of seed flavonoids. Annu Rev Plant Biol 57: 405-430. doi:https://doi.org/10.1146/annurev.arplant.57.032905.105252. PubMed: 16669768.
- 31. Stracke R, Ishihara H, Huep G, Barsch A, Mehrtens F et al. (2007) Differential regulation of closely related R2R3-MYB transcription factors controls flavonol accumulation in different parts of the Arabidopsis thaliana seedling. Plant J 50: 660-677. doi:https://doi.org/10.1111/j.1365-313X.2007.03078.x. PubMed: 17419845.
- 32. Liu R, Lü B, Wang X, Zhang C, Zhang S et al. (2010) Thirty-seven transcription factor genes differentially respond to a harpin protein and affect resistance to the green peach aphid in Arabidopsis. J Biosci 35: 435-450. doi:https://doi.org/10.1007/s12038-010-0049-8. PubMed: 20826953.
- 33. Kaur H, Heinzel N, Schöttner M, Baldwin IT, Gális I (2010) R2R3-NaMYB8 regulates the accumulation of phenylpropanoid-polyamine conjugates, which are essential for local and systemic defense against insect herbivores in Nicotiana attenuata. Plant Physiol 152: 1731-1747. PubMed: 20089770.
- 34. Jin H, Martin C (1999) Multifunctionality and diversity within the plant MYB-gene family. Plant Mol Biol 41: 577-585. doi:https://doi.org/10.1023/A:1006319732410. PubMed: 10645718.
- 35. Singh K, Kumar S, Rani A, Gulati A, Singh AP (2009) Phenylalanine ammonia-lyase (PAL) and cinnamate-4-hydroxylase (C4H) and catechins(flavan-3-ols) accumulation in tea. Funct Integr Genomics 9: 125–134. doi:https://doi.org/10.1007/s10142-008-0092-9. PubMed: 18679731.
- 36. Li X, Li S, Lin J (2003) Effect of GA3 spraying on lignin and auxin contents and the correlated enzyme activities in bayberry (Myrica rubra Bieb.) during flower bud induction. Plant Sci 164: 549-556. doi:https://doi.org/10.1016/S0168-9452(03)00004-9.
- 37. Russell DW, Galston AW (1969) Blockange by gibberellic acid on phytochrome effects on growth, auxin responses, and flavonoid biosynthesis in etiolated pea internodes. Plant Physiol 44: 1211–1216. doi:https://doi.org/10.1104/pp.44.9.1211. PubMed: 16657193.
- 38. Hinderer W, Peterson M, Seitz HU (1983) Inhibition of flavonoid biosynthesis by gibberellic acid in cell suspension cultures of Daucus carota. Planta 160: 544-549.
- 39. Devaiah BN, Madhuvanthi R, Karthikeyan AS, Raghothama KG (2009) Phosphate starvation responses and gibberellic acid biosynthesis are regulated by the MYB62 transcription factor in Arabidopsis. Mol Plants 2: 43–58. doi:https://doi.org/10.1093/mp/ssn081. PubMed: 19529828.
- 40. Cheng H, Song S, Xiao L, Soo HM, Cheng Z et al. (2009) Gibberellin acts through jasmonate to control the expression of MYB21, MYB24, and MYB57 to promote stamen filament growth in Arabidopsis. PLOS Genet 5, e1000440. PubMed: 19325888.
- 41. Aya K, Ueguchi-Tanaka M, Kondo M, Hamada K, Yano K et al. (2009) Gibberellin modulates anther development in rice via the transcriptional regulation of GAMYB. Plant Cell 21: 1453-1472. doi:https://doi.org/10.1105/tpc.108.062935. PubMed: 19454733.
- 42. Zhao C, Hanada A, Yamaguchi S, Kamiya Y, Beers EP (2011) The Arabidopsis Myb genes MYR1 and MYR2 are redundant negative regulators of flowering time under decreased light intensity. Plant J 66: 502-515. doi:https://doi.org/10.1111/j.1365-313X.2011.04508.x. PubMed: 21255164.
- 43. Goff SA, Cone KC, Chandler VL (1992) Functional analysis of the transcriptional activator encoded by the maize B-gene: evidence for a direct functional interaction between two classes of regulatory proteins. Genes Dev 6: 864-875. doi:https://doi.org/10.1101/gad.6.5.864. PubMed: 1577278.
- 44. Craven-Bartle B, Pascual MB, Cánovas FM, Avila CA (2013) A Myb transcription factor regulates genes of the phenylalanine pathway in maritime pine. Plant J, 74: 755–66. doi:https://doi.org/10.1111/tpj.12158. PubMed: 23451763.
- 45. Lee MM, Schiefelbein J (2001) Developmentally distinct MYB genes encode functionally equivalent proteins in Arabidopsis. Development 128: 1539-1546. PubMed: 11290293.
- 46. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389-3402. doi:https://doi.org/10.1093/nar/25.17.3389. PubMed: 9254694.
- 47. Mulder N, Apweiler R (2007) InterPro and InterProScan: tools for protein sequence classification and comparison. Methods Mol Biol 396: 59-70. doi:https://doi.org/10.1007/978-1-59745-515-2_5. PubMed: 18025686.
- 48. Conesa A, Götz S (2008) Blast2GO: A comprehensive suite for functional analysis in plant genomics. Int J Plants Genomics, 2008: 2008: 619832. PubMed: 18483572.
- 49. Hall TA (1999) BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp 41: 95-98.
- 50. Bjellqvist B, Basse B, Olsen E, Celis JE (1994) Reference points for comparisons of two-dimensional maps of proteins from different human cell types defined in a pH scale where isoelectric points correlate with polypeptide compositions. Electrophoresis 15: 529-539. doi:https://doi.org/10.1002/elps.1150150171. PubMed: 8055880.
- 51. Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22: 4673-4680. doi:https://doi.org/10.1093/nar/22.22.4673. PubMed: 7984417.
- 52. Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol, 24: 1596-1599. doi:https://doi.org/10.1093/molbev/msm092. PubMed: 17488738.
- 53. Kumar S, Nei M, Dudley J, Tamura K (2008) MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform 9: 299-306. doi:https://doi.org/10.1093/bib/bbn017. PubMed: 18417537.
- 54. Li DF, Zhang Y, Hu X, Shen X, Ma L et al. (2011) Transcriptional profiling of Medicago truncatula under salt stress identified a novel CBF transcription factor MtCBF4 that plays an important role in abiotic stress responses. BMC Plant Biol 11: 109. doi:https://doi.org/10.1186/1471-2229-11-109. PubMed: 21718548.
- 55. Horsch RB, Rogers SG, Fraley RT (1985) Transgenic plants. Cold Spring Harb Symp Quant Biol 50: 433-437. doi:https://doi.org/10.1101/SQB.1985.050.01.054. PubMed: 3868487.