Skip to main content

Functional identification of PsMYB57 involved in anthocyanin regulation of tree peony



R2R3 myeloblastosis (MYB) genes are widely distributed in plants and comprise one of the largest transcription factor gene families. They play important roles in the regulatory networks controlling development, metabolism, and stress responses. Researches on functional genes in tree peony are still in its infancy. To date, few MYB genes have thus far been reported.


In this study, we constructed a comprehensive reference gene set by transcriptome sequencing to obtain R2R3 MYB genes. The transcriptomes of eight different tissues were sequenced, and 92,837 unigenes were obtained with an N50 of 1662 nt. A total of 48,435 unigenes (77.98%) were functionally annotated in public databases. Based on the assembly, we identified 57 R2R3 MYB genes containing full-length open reading frames, which clustered into 35 clades by phylogenetic analysis. PsMYB57 clustered with anthocyanin regulation genes in Arabidopsis and was mainly transcribed in the buds and young leaves. The overexpression of PsMYB57 induced anthocyanin accumulation in tobacco, and four detected anthocyanin structural genes, including NtCHS, NtF3’H, NtDFR, and NtANS, were upregulated. The two endogenous bHLH genes NtAn1a and NtAn1b were also upregulated and may work in combination with PsMYB57 in regulating anthocyanin structural genes.


Our study offers a useful reference to the selection of candidate MYB genes for further functional studies in tree peony. Function analysis of PsMYB57 is helpful to understand the color accumulation in vegetative organs of tree peony. PsMYB57 is also a promising resource to improve plant color in molecular breeding.


Myeloblastosis (MYB) transcription factors (TFs) are widely distributed in all eukaryotes and are one of the largest families of TFs in plants. MYB TFs are characterized by a conserved DNA binding domain (DBD) known as the MYB domain near the N-terminus. The MYB domain is comprised of up to four incomplete repeats, each of which consists of about 50–53 amino acids and forms three α-helices [1, 2]. The second and third helices form a helix-turn-helix (HTH) structure, and the third α-helix of each repeat is a DNA-recognition helix that binds specific DNA sequences in the major groove [3, 4]. Typically, the MYB repeat contains three regularly spaced tryptophans, which stabilize the structure of the MYB domain by forming a hydrophobic core [5]. In plants, the first Trp of the third repeat is always substituted by phenylalanine (Phe) or isoleucine (Ile) [3, 6, 7]. Compared with the conservative MYB domain, the C-terminal regions of MYB proteins are more flexible. They typically function as transacting domains (TADs) and are responsible for the regulatory activity of the protein [8].

Based on the repeat number in the DBD domain, MYB proteins are classified into four subfamilies, including R2R3-MYB, R1R2R3-MYB, 4R-MYB, and MYB-related. R2R3-MYB proteins contain two repeats in the DBD domain and constitute the largest group of the MYB family. Currently, 8746 R2R3 MYB sequences are available in the plant transcription factor database ( In some plants with available complete genome sequences, the R2R3 MYB family gene has been systematically studied. 126 R2R3 MYBs were reported in Arabidopsis [9], 123 in Jatropha curcas [10], 94 in pineapple [11], and 128 in peach [12]. Since the first plant MYB gene C1 was cloned from maize, functional characterization of the R2R3 MYB gene has been conducted in various plants. A summary of MYB gene function in Arabidopsis showed that these genes play important roles in multiple plant-specific processes, such as primary and secondary metabolism, cell fate and identity, developmental processes, and responses to biotic and abiotic stress [9].

The anthocyanin pathway is a branch of the flavonoid pathway and includes two clusters of coregulated structural genes: early biosynthetic genes (EBG) and late biosynthetic genes (LBG). In the R2R3 MYB family, a small group of members are key factors in regulating anthocyanin biosynthesis and usually interact with bHLH factors to form complexes that regulate structural genes by binding to their promoters. In Arabidopsis, three MYBs, including MYB11, MYB12, and MYB111, mainly regulate EBGs and control the biosynthesis of flavonols. Four MYBs, including MYB75, MYB90, MYB113, and MYB114, mainly regulate LBGs. In nature, some tissue colors in plants are caused by the mutation of MYB genes. For example, a mutation in the promoters of Pr in purple cauliflower [13], MdMYB10 in red flesh apple [14], and Ruby in blood oranges [15] leads to the heterotopic accumulation of anthocyanins. Recently, the overexpression of members in non-model plants have been reported to lead to high levels of anthocyanin accumulation, such as PtrMYB119 in Populus trichocarpa [16], DcMYB6 in purple carrot [17], MaAN2 in grape hyacinth [18], and OjMYB1 in Oenanthe javanica [19].

Tree peony are among the most valued ornamental flower with a variety of colors in China. Researches on functional genes in tree peony are still in its infancy. To date, few studies in the MYB gene family have been reported. Due to the unpublished and highly complex genome of tree peony, transcriptome sequencing has become a suitable method for obtaining gene sequences. Several transcriptomic studies of tree peony have been made [20,21,22]. However, such studies have not provided a global genetic overview. In this study, we performed global transcriptome sequencing and constructed a reference gene sequence library of tree peony. The R2R3 MYB genes were screened and systematically analyzed in order to provide a foundation for further functional gene research. The anthocyanin regulatory MYB gene was functionally characterized, which will facilitate improvements in plant color in tree peony.


Tissue collection and total RNA extraction

Paeonia suffruticosa Andr. cv. Er Qiao is a traditional variety in China. Erqiao plants were kindly identified and provided by Luoyang Research Institute of Peony (Luoyang, China), and were grown under field conditions. From five healthy plants, we collected samples of bud, young leaf, red petal, pink petal, petal spot, stamen, pistil, and seed. The samples of each tissue type were pooled. The samples were rapidly frozen in liquid nitrogen and stored at − 80 °C. Total RNA extraction was performed using the CTAB-LiCl method [23], and the contamination of genomic DNA was eliminated using DNase I (TaKaRa, Dalian, China).

Sequencing and de novo assembly

The cDNA libraries were constructed and sequenced at BGI Co., Ltd. (Shenzhen, China), and sequencing was performed on a BGISEQ-500 platform. The raw sequence data reported in this paper have been deposited in the Genome Sequence Archive [24] in BIG Data Center [25] Beijing Institute of Genomics (BIG), Chinese Academy of Sciences, under accession numbers CRA001327 that are publicly accessible at The raw reads were filtered by removing the adapter sequences, low quality reads, and reads with more than 20% Q < 20 bases. The remaining clean reads generated from each sample were de novo assembled using Trinity [26], and all assembled datasets were further clustered using TGICL to generate a Nr unigene set [27]. The ORFs of all unigenes were predicted using Transdecoder v3.0.1.

Gene annotation

To predict the functions of the genes, the assembled unigenes were used to query seven public databases using BLAST with an E-value of <1e− 5, including the Nr, nucleotide (Nt), SWISS-PROT, InterPro, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and eukaryotic orthologous group (KOG). Based on the Nr annotation, Blast2GO was used to obtain GO annotations [28], and WEGO software was used to perform GO functional classification for all unigenes [29]. Pathway assignments were also carried out based on the KEGG database.

Identification of R2R3-MYB genes in tree peony

To identify the maximum number of MYB domain-containing sequences, de novo assembled unigenes were used to construct a local database using the BLAST program, and the HMM profile of the MYB DNA-binding domain (PF00249) downloaded from Pfam was used to isolate all possible homologs in tree peony with HMMER 3.0 [30]. The default parameters were adopted, and the cutoff value was set to 0.001. The redundant proteins were removed and sequences without complete ORFs were discarded. The remaining sequences were detected using PROSITE ( and the SMART online service ( to eliminate any sequences that do not contain the R2R3 MYB domain. Calculation of protein isoelectric point and molecular weight was performed using the ExPASy proteomics server ( Subcellular localization of proteins was predicted using the online service of Cell-PLoc 2.0 (

Gene sequence analysis

For conservation analysis of the R2R3 MYB domain in tree peony, all obtained R2R3 MYB proteins were aligned using ClustalW, and the sequence logos for the R2 and R3 repeats were generated using the Weblogo v3 online tool [31]. To examine the phylogenetic relationships and evolutionary history of the MYB gene family, full-length proteins of the MYB protein from tree peony and Arabidopsis were used to generate a phylogenetic tree based on the neighbor-joining method using MEGA5.0 software [32], and tree nodes were evaluated with 1000 bootstrap replicates. The phylogenetic tree was further modified using the iTOL online service ( Sequences and functional annotations of all Arabidopsis MYB proteins were obtained from the Arabidopsis Information Resource (TAIR) database.

Expression analysis by qRT-PCR

Total RNA was extracted as described above, and the first-strand cDNA was reversed-transcribed with 1 μg of total RNA using the PrimeScript™ 1st Strand cDNA Synthesis Kit (TaKaRa, Dalian, China). qRT-PCR experiments were performed using SYBR® Premix Ex Taq™ II kit (TaKaRa, Dalian, China) on an ABI7500 system according to the manufacturer’s instructions. The thermal cycling conditions were as follows: an initial heat denaturing step at 95 °C for 3 min; followed by 40 cycles of 95 °C for 20 s, 58 °C for 20 s, and 72 °C for 20 s. Each sample was amplified with three biological replicates and three technical replicates. Gene transcription levels were calculated using the 2−∆∆CT comparative cycle threshold method [33]. Ubiquitin and tubulinA1 were used as an internal control for Paeonia suffruticosa [34] and tobacco [18], respectively, to normalize the relative expression levels of the analyzed genes. Primers used for PCR amplification are shown in Table S1.

Plasmid construction and tobacco transformation

To verified assembly genes, ten paired primers were designed to amplifacte randomly selected MYB genes (Table S2). Each of the purified PCR products was ligated to pMD19-T vector and transferred into E.coil DH5α competent cells. Three clones of each gene were sequenced using sanger sequencing method. The coding region sequence of PsMYB57 was amplified by PCR using the following primers: forward primer: 5′-ACGCGTCGACATGGAGGGAATGTTAGGATTGAGAAAAG-3′, reverse primer: 5′-ATTGGCTGCAGTTAATTCACCGCTTGCCCTTCAGCACTTAATAG-3′. Following agarose gel electrophoresis, the PCR product was digested with the SalI and PstI enzymes. The target fragment was linked to the pCAMBIA1300–35 s vector. The recombinant plasmid was transferred into Agrobacterium strain Gv3101 by electrical transformation. Tobacco plants (Nicotiana tabacum L. cv. SR1) were transformed with Agrobacterium using a previously described protocol [35], and transgenic lines showing color changes in the leaves were used for further analysis.


Transcriptome sequencing and assembly

As the genomic information was very limited for tree peony, we performed a global transcriptome sequencing to identify MYB gene family members. Eight libraries were constructed from the total RNA of the tissues, including the bud, young leaf, pink petal, red petal, petal spot, pistil, stamen, and seed. After filtering, each library generated over 6.5 Gb data with Q20 scores over 96.5%. A total of 52.82 Gb data were obtained (Table S3). The data were de novo assembled using Trinity. After further assembly with TGICL, a total of 92,837 unigenes were generated with an N50 of 1662 nt. Among the assembled unigenes, 49,376 were longer than 500 bp, 10,543 were longer than 1000 bp, and 1600 were longer than 1500 bp.

Functional annotation

Gene functions were predicted by querying seven public databases, and a total of 48,435 unigenes (77.98%) were functionally annotated (Table S4). Among them, 45,046 unigenes (72.52%) obtained hits in the non-redundant (Nr) database, 31,786 obtained hits in the SWISS-PROT database, and 33,713 unigenes obtained hits in the InterPro database. Based on the Nr annotation, the top sequence matches obtained from BLASTX are shown in Fig. S1. The sequences were most similar to Vitis vinifera (35.2%), followed by Nelumbo nucifera (6.14%), Theobroma cacao (6.12%), Jatropha curcas (3.99%), and Prunus mume (3.39%).

Identification of R2R3 MYB genes in tree peony

To identify R2R3 MYB genes in tree peony, the MYB DNA-binding domain (PF00249) was queried in the de novo assembly using HMMER 3.0, and hits with full-length open reading frames (ORFs) were selected. The redundant sequences were removed and the MYB domain of the remaining sequences was further verified using SMART. A total of 57 R2R3 MYB genes were ultimately identified. As two R2R3 MYB genes in our transcriptome have previously been deposited in Genbank as PsMYB1 and PsMYB2, the remaining 55 genes were provisionally named PsMYB3 to PsMYB57 (Genbank accession numbers: MK377190- MK377244). To evaluate the assembly genes, ten randomly selected MYB genes were cloned and sequenced using sanger sequencing method. Results showed that the MYB genes had 97 to 100% identity with assembly genes. In general, the sanger sequencing results were in rough accordance with the electronic data of assembly. The MYB proteins ranged from 170 (PsMYB33) to 560 (PsMYB9) amino acids in length (Table S5). The predicted isoelectric point of the R2R3-MYB proteins ranged from 4.86 (PsMYB43) to 10.22 (PsMYB3). Subcellular localization analysis revealed that all 57 MYB proteins were localized in the nucleus.

To investigate the MYB domain features and the conservation of amino acids, sequence logos were generated using the aligned R2 and R3 motifs (Fig. 1). The results indicated an even distribution of a series of highly conserved tryptophan residues (W), which are considered as a hallmark of the MYB domain. In the R2 domain, three highly-conserved W were present at positions 5, 26, and 48, but tyrosine (Y) replaced the W at position 48 in PsMYB38. While two highly conserved W exist at positions 81 and 100 in the R3 domain, the first W residue was generally replaced by phenylalanine (F), isoleucine (I), and leucine (L), and substitutions with the amino acid methionine (M) were also observed in PsMYB52. Furthermore, F and Y substituted the W at the W-100 position in PsMYB17 and PsMYB18, respectively. In addition to the highly conserved W residues, D-10, L-13, G-21, C-44, and R-47 in the R2 repeat, and E-66, G-78, R-91, and T-92 in the R3 repeat were completely conserved. These results reveal that R2R3-MYB proteins are highly conserved in tree peony.

Fig. 1

The R2 and R3 MYB repeats are highly conserved across all R2R3 MYB proteins in tree peony. The sequence logos of the R2 (a) and R3 (b) MYB repeats are based on full-length alignments of all R2R3 MYB proteins. The bit score exhibits the information content for each position in the sequence. Asterisks indicate the conserved tryptophan residues (Trp) in the MYB domain

Phylogenetic analysis of the MYB proteins

To explore the putative function of MYBs in tree peony, a phylogenetic tree was constructed using full-length R2R3 MYBs, including 57 proteins from tree peony and 124 proteins from Arabidopsis. As is shown in Fig. 2, a total of 181 MYBs were clustered into 35 clades, whereas the proteins PsMYB46 and AtMYB139 were not grouped into any clade. In 26 of the 35 clades, members from both tree peony and Arabidopsis were present. However, in some clades, MYBs were unequally represented. For example, C25, C32, and C33 included just one or two PsMYBs but at least seven AtMYBs, while PsMYB members were more abundant than that in Arabidopsis in the C2 and C4 clades. We also detected some species-specific proteins. There were no PsMYB grouped in 7 clade, including clade 9, 10, 14, 21, 22, 26 and 27, whereas clade C7 contained no members from Arabidopsis. Overall, our classification of the MYBs corroborates that of Arabidopsis.

Fig. 2

Phylogenetic analysis of R2R3 MYB family genes in tree peony and Arabidopsis. The full-length amino acid sequences of MYB proteins were aligned using CLUSTALW and the phylogenetic tree was constructed using the neighbor-joining method in MEGA5.0. Subgroup names are included next to each clade together with a short name, and the putative functions of MYB proteins in tree peony are annotated alongside the names. Green and orange shading are used to distinguish the different clades

According to the gene annotation summary in Arabidopsis [9], the functions of PsMYBs in different clades were predicted and were found to be involved in a broad range of physiological functions. The functions of the three clades were associated with the regulation of secondary metabolism, including anthocyanin biosynthesis (C1), proanthocyanidins (PAs) and tannis (C4), and flavonol biosynthesis (C5); three clades were associated with determining cell fate and identity, including trichome initiation (C3), trichome branching (C13), and stomatal differentiation (C35); five clades were associated with plant developmental processes, including axillary meristem formation (C23), anther development (C24), stamen development (C25), later organ separation (C32), and shoot morphogenesis (C34); and five clades were associated with responses to biotic and abiotic stress, including abscisic acid response (C15), hormone response (C28, C29), and stress response (C12, C30).

R2R3 MYB genes involved in anthocyanin biosynthesis and their expression profiles in different tissues

The alignment of R2R3-MYBs involving in anthocyanins regulation revealed a highly conserved R2R3 repeat domain in the N-terminal of PsMYB57, and a [DE]Lx2[RK]x3Lx6Lx3R motif, which is necessary for interaction with bHLH proteins, was present in the R3 repeat (Fig. S2a). While in the C-terminal variable region, we found a signature motif KPXPR(S/T) F, which is specific to MYB proteins that activate anthocyanin biosynthesis. The phylogenetic analysis of R2R3-MYBs from different plant species showed that PsMYB57 was clustered into the clade of anthocyanin, and was more closely to Ruby, a MYB transcription factor involved in the anthocyanin pathway in Citrus sinensis (Fig. S2b). The R2R3 domain and completed protein sequence of PsMYB57 shared 88 and 51% identity with Ruby, respectively. The sequence analysis indicated that PsMYB57 was potentially involved in the regulation of the anthocyanin synthetic pathway as well, which was consistent with the analysis of MYB gene family in tree peony. In ‘Er Qiao’, pigment was mainly accumulated in tissues of bud, young leaf, red petal and petal spot (Fig. 3a-h). As shown in Fig. 3i, PsMYB57 was mainly expressed in bud and young leaf, but its transcript was not detectable or expressed at a very low level in other tissues.

Fig. 3

Tissues color of tree peony and transcription level analysis of PsMYB57. a to h display colors of bud, young leaf, red petal, pink petal, petal spot, pistil, stamen and seed of tree peony. i The expression profile of PsMYB57 in different tissues of tree peony. Ubiquitin was the reference gene to normalize the expression of PsMYB57. Each column represents means ±SD from three independent experiments

Overexpression of PsMYB57 induces anthocyanin accumulation in tobacco

To verify the function of putative anthocyanin regulation genes, coding region sequence of PsMYB57 was transferred into tobacco under the control of the cauliflower mosaic virus 35S promoter. In the process of obtaining hygromycin-resistant callus, no colored callus was found. The transgenic lines usually display red patches on the leaves (Fig. 4 a–h); we found that three lines displayed significant red patches on the leaves, while eight lines displayed only light red patches. None of the 11 lines had red pigments in the young leaves. In the reproductive organs, the transgenic lines showed red coloring in the sepals and pericarps (Fig. 4 i–l), but no obvious change in flower color was observed. Anthocyanin content was determined in the different plants with dark red or light red patches, and the transgenic lines contained higher levels of anthocyanins, whereas anthocyanins were not detected in the control plants (Fig. 4 m).

Fig. 4

Phenotypic observation and anthocyanin content in PsMYB57 transgenic tobacco lines. a to c are the PsMYB57 transgenic tobacco lines; d is the wild-type tobacco line; e to g are the leaves of the transgenic tobacco lines, h is leaf of wild-type tobacco line; i is three transgenic flowers; j is the wild-type flower; k is the sepals of the transgenic lines; l is the sepals of the wild-type lines; and m is the total anthocyanin content of the leaves in the wild-type and three transgenic lines. Three biological replicates were performed for anthocyanins detection and the asterisk indicates statistical significance (P < 0.05)

PsMYB57 promotes the expression of anthocyanin pathway genes in tobacco

Quantitative real-time PCR (qRT-PCR) analysis verified that PsMYB57 was transcribed in the transgenic lines but was not detected in the control plants. In T1 generation of the three PsMYB57 transgenic tobacco lines, genes transcript levels was detected in tissues of leaves, sepals and pericarps (Fig. 5). Results showed that the transcript levels of one EBG gene (NtCHS) and two LBG genes (NtDFR and NtANS) in the anthocyanin pathway were remarkably increased. NtF3’H was highly expressed in detected tissues of PsMYB57 transgenic lines, but showed only a slight increase in sepals of transgenic line 3. We also detected the transcript levels of two endogenous bHLH genes, NtAn1a and NtAn1b, which facilitate the regulation of anthocyanin biosynthesis by MYB in tobacco [36]. The transcription level of NtAn1a was not detectable in tissues of wild-type plant. The two bHLHs were upregulated in the tissues of all three detected transgenic lines. Overall, structural anthocyanin genes and bHLH genes showed similar trends as PsMYB57 in the transgenic lines.

Fig. 5

Transcript levels of structural anthocyanin genes and endogenous bHLH genes in wild-type and transgenic tobacco. a to f are transcription levels of NtCHS, NtF3’H, NtDFR and NtANS, NtAn1a and NtAn1b in leaves, sepals and pericarps of wild-type and transgenic tobacco lines. Transcription level of NtAn1a was not detectable in wild-type lines, and transgenic line 3 used as control in (e). CK was wild-type tobacco lines, Line 1 to Line 3 were three PsMYB57 transgenic tobacco lines. Three T1 generation lines from each transgenic tobacco were collected as samples. TubulinA1 was used as an internal control to normalize the relative expression levels of the analyzed genes. qRT-PCR analysis was performed with three biological and technical replicates per experiment. An asterisk indicates statistical significance (P < 0.05)


With the rapid development of high-throughput sequencing technology, whole-genome sequences of more animals and plants have become available. However, the genomes of the majority of species have not yet been sequenced. Furthermore, the genomes of some animals and plants are large and complex and thus pose significant challenges to the sequencing work. To promote functional gene research in species without reference genomes, it is necessary to establish as complete a gene reference database as possible. For instance, in Chinese giant salamander [37], whitefly [38], and Spodoptera frugiperda [39], the transcriptomes of different tissues were used to construct a reference gene set. Tree peony has a large genome of about 13 Gb, but complete genome information has not yet been available. Furthermore, MYB genes in tree peony are highly limited in the Genbank database. To identify as many MYB genes as possible, eight different tree peony tissues were collected and their transcriptomes sequenced to construct a reference gene library, which included 92,837 assemble unigenes. Our data not only provide comprehensive gene sequences for identifying as many MYB genes as possible, but also offer a valuable resource to future biological research in tree peony.

Based on the reference transcriptome, a total of 57 MYBs were identified in tree peony. Our observed MYBs were lower than those reported in plants with reference genomes, as only MYBs with complete ORFs were analyzed for subsequent research on gene function. Furthermore, some members with low transcription levels or with spatiotemporal-specific transcription may not be recovered. We found a series of evenly distributed and highly conserved Trp residues in our 57 MYB members, which have been reported in pineapple [11], watermelon [40], and sweet orange [6]. A conserved Trp residue is regarded as a landmark of the MYB domain and plays key roles in sequence-specific DNA binding [5, 41]. In addition to the conserved Trp residues, there are nine amino acids that were completely conserved in different positions. We suggest that the MYB domain in tree peony is highly evolutionarily conserved.

The functions of many R2R3 MYBs have been characterized and several members are associated with the control of plant-specific processes, including primary and secondary metabolism, developmental processes, cell fate and identity, and stress responses [9]. Based on the phylogenetic analysis, the 57 MYBs in tree peony and 124 MYBs in Arabidopsis were classified into 35 clades. In general, members within the same subfamily are assumed to have recent common evolutionary origins and usually display similar sequences and functions [42]. Our results provide a reference for exploring the functions of MYBs. For example, PsMYB57 was grouped together with Arabidopsis AtMYB75 and AtMYB90 in clade 1, which is associated with anthocyanin regulation [43]. PsMYB54 and PsMYB55 were grouped into clade 3 with the Arabidopsis AtMYB0, AtMYB66, and AtMYB33, which represent the functional clade of trichome initiation [44,45,46]. PsMYB10, PsMYB50, and PsMYB51 were found in group C29 and are associated with stress responses [47]. Two clades (C6 and C9) do not possess Arabidopsis MYB members, suggesting that these proteins might have specialized roles that were either lost in Arabidopsis or gained after divergence from the last common ancestor.

Flower color is an important economic trait that attracts many breeders in tree peony. In numerous plants, MYB genes have been verified to regulate the flavonoid pathway and are responsible for tissue color. The heterologous expression of MYB genes in tobacco drastically elevates anthocyanin contents in the leaves, as does VlMYBA2 in grape [48] and MaAN2 in grape hyacinth [18]. Additionally, the stable tobacco transformants of MdMYBA, which control anthocyanin biosynthesis in apple, do not accumulate pigments in the leaves [49]. MYB usually interacts with a bHLH protein to form a complex to regulate anthocyanin structural genes, and a previous study [50] suggested that the absence of anthocyanins in MdMYBA transgenic lines is due to a lack of a bHLH partner. In our study, the overexpression of PsMYB57 upregulated structural anthocyanin genes and two endogenous bHLH genes, NtAn1a and NtAn1b, in tobacco. We suggested that NtAn1a and NtAn1b are the partners of PsMYB57 in the regulation of structural genes. We note that the leaf color of the PsMYB57 transgenic lines differed from the reported lines transformed with VlMYBA2 and MaAN2, in which the full leaves were dark red in color. One possible explanation is that the PsMYB57-NtbHLH complex may be unstable or have lower activity, and a tree peony bHLH partner is necessary for the full functioning of PsMYB57 in tobacco. Overall, the overexpression of PsMYB57 upregulated structural anthocyanin genes in tobacco, leading to anthocyanin accumulation. The results indicated that PsMYB57 is an important anthocyanin regulator in tree peony.

In the tree peony variety ‘Er Qiao’, the bud, young leaf, and flower tissues typically display a red pigment, and the petal spot is a dark red color. However, PsMYB57 was mainly expressed in the bud and young leaf. We suggested that PsMYB57 mainly regulates color in the bud and young leaves, and thus unknown MYB genes might exist to take charge of flower color.


We constructed a comprehensive reference transcriptome of tree peony, including 92,837 unigenes, and identified 57 R2R3-MYB genes containing full-length ORFs. Functional prediction of the R2R3 MYBs verified that PsMYB57 regulates anthocyanin biosynthesis in the buds and young leaves in tree peony. Our study provides a basis for the exploitation of candidate MYB genes for the genetic engineering of tree peony. Furthermore, PsMYB57 can contribute to further researches on the improvement of flower colorization.

Availability of data and materials

The datasets supporting the conclusions of this article are available in the Genome Sequence Archive in BIG Data Center, Beijing Institute of Genomics (BIG), Chinese Academy of Sciences, with the accession numbers CRA001327 (, and the GenBank data libraries with the accession numbers MK377190-MK377244.





Transcription factors;


Quantitative real-time PCR.


DNA binding domain






Early biosynthetic gene


Late biosynthetic gene


Transacting domains


Gene Ontology


Kyoto Encyclopedia of Genes and Genomes


Eukaryotic orthologous group




  1. 1.

    Rabinowicz PD, Braun EL, Wolfe AD, Bowen B, Grotewold E. Maize R2R3 Myb genes: sequence analysis reveals amplification in the higher plants. Genetics. 1999;153(1):427.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Stracke R, Werber M, Weisshaar B. The R2R3-MYB gene family in Arabidopsis thaliana. Curr Opin Plant Biol. 2001;4(5):447–56.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Ambawat S, Sharma P, Yadav NR, Yadav RC. MYB transcription factor genes as regulators for plant responses: an overview. Physiol Mol Biol Plants. 2013;19(3):307–21.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Williams CE, Grotewold E. Differences between plant and animal Myb domains are fundamental for DNA binding activity, and chimeric Myb domains have novel DNA binding specificities. J Biol Chem. 1997;272(1):563–71.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Ogata K, Hojo H, Aimoto S, Nakai T, Nakamura H, Sarai A, Ishii S, Nishimura Y. Solution structure of a DNA-binding unit of Myb: a helix-turn-helix-related motif with conserved tryptophans forming a hydrophobic core. P Natl Acad Sci. 1992;89(14):6428.

    CAS  Article  Google Scholar 

  6. 6.

    Hou XJ, Li SB, Liu SR, Hu CG, Zhang JZ. Genome-wide classification and evolutionary and expression analyses of Citrus MYB transcription factor families in sweet Orange. PLoS One. 2014;9(11):e112375.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    He Q, Jones DC, Li W, Xie F, Ma J, Sun R, Wang Q, Zhu S, Zhang B. Genome-wide identification of R2R3-MYB genes and expression analyses during abiotic stress in Gossypium raimondii. Sci Rep. 2016;6:22980.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Martin C, Paz-Ares J. MYB transcription factors in plants. Trends Genetics. 1997;13(2):67–73.

    CAS  Article  Google Scholar 

  9. 9.

    Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010;15(10):573–81.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Peng X, Liu H, Wang D, Shen S. Genome-wide identification of the Jatropha curcas MYB family and functional analysis of the abiotic stress responsive gene JcMYB2. BMC Genomics. 2016;17(1):251.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Liu C, Xie T, Chen C, Luan A, Long J, Li C, Ding Y, He Y. Genome-wide organization and expression profiling of the R2R3-MYB transcription factor family in pineapple (Ananas comosus). BMC Genomics. 2017;18(1):503.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Zhang C, Ma R, Xu J, Yan J, Guo L, Song J, Feng R, Yu M. Genome-wide identification and classification of MYB superfamily genes in peach. PLoS One. 2018;13(6):e0199192.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Chiu L-W, Zhou X, Burke S, Wu X, Prior RL, Li L. The purple cauliflower arises from activation of a MYB transcription factor. Plant Physiol. 2010;154(3):1470.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Espley RV, Hellens RP, Putterill J, Stevenson DE, Kutty-Amma S, Allan AC. Red colouration in apple fruit is due to the activity of the MYB transcription factor, MdMYB10. Plant J. 2007;49(3):414–27.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Butelli E, Licciardello C, Zhang Y, Liu J, Mackay S, Bailey P, Reforgiato-Recupero G, Martin C. Retrotransposons control fruit-specific, cold-dependent accumulation of anthocyanins in blood oranges. Plant Cell. 2012;24(3):1242–55.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Cho JS, Nguyen VP, Jeon H-W, Kim M-H, Eom SH, Lim YJ, Kim WC, Park EJ, Choi YI, Ko JH. Overexpression of PtrMYB119, a R2R3-MYB transcription factor from Populus trichocarpa, promotes anthocyanin production in hybrid poplar. Tree Physiol. 2016;36(9):1162–76.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Xu ZS. Feng K, Que F, Wang F, Xiong AS: a MYB transcription factor, DcMYB6, is involved in regulating anthocyanin biosynthesis in purple carrot taproots. Sci Rep. 2017;7:45324.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Chen K, Liu H, Lou Q, Liu Y. Ectopic Expression of the Grape Hyacinth (Muscari armeniacum) R2R3-MYB Transcription Factor Gene, MaAN2, Induces Anthocyanin Accumulation in Tobacco. Front Plant Sci. 2017;8(965).

  19. 19.

    Feng K, Xu ZS, Que F, Liu JX, Wang F, Xiong AS. An R2R3-MYB transcription factor, OjMYB1, functions in anthocyanin biosynthesis in Oenanthe javanica. Planta. 2018;247(2):301–15.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Zhou H, Cheng F-Y, Wang R, Zhong Y, He C. Transcriptome comparison reveals key candidate genes responsible for the unusual Reblooming trait in tree peonies. PLoS One. 2013;8(11):e79996.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Zhang C, Wang Y, Fu J, Dong L, Gao S, Du D. Transcriptomic analysis of cut tree peony with glucose supply using the RNA-Seq technique. Plant Cell Rep. 2014;33(1):111–29.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Zhang Y, Cheng Y, Ya H, Xu S, Han J. Transcriptome sequencing of purple petal spot region in tree peony reveals differentially expressed anthocyanin structural genes. Front Plant Sci. 2015;6:964.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Gambino G, Perrone I, Gribaudo I. A rapid and effective method for RNA extraction from different tissues of grapevine and other woody plants. Phytochem Anal. 2008;19(6):520–5.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Wang Y, Song F, Zhu J, Zhang S, Yang Y, Chen T, Tang B, Dong L, Ding N, Zhang Q. GSA: genome sequence archive. Genom Proteom Bioinf. 2017;15(1):14–8.

    Article  Google Scholar 

  25. 25.

    BIG Data Center Members: Database Resources of the BIG Data Center in 2018. Nucleic Acids Res. 2018, 46(D1):D14-D20. Doi:

  26. 26.

    Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, et al. TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003;19(5):651–2.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006;34(suppl 2):W293–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Mistry J, Finn RD, Eddy SR, Bateman A, Punta M. Challenges in homology search: HMMER3 and convergent evolution of coiled-coil regions. Nucleic Acids Res. 2013;41(12):e121.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Crooks GE, Hon G, Chandonia J-M, Brenner SE. WebLogo: a sequence logo generator. Genome Res. 2004;14(6):1188–90.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Livak KJ, Schmittgen TD. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods. 2001;25(4):402–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    YanJie W, Li D, Chao Z, XiaoQing W. Reference gene selection for real-time quantitative PCR normalization in tree peony (Paeonia suffruticosa Andr.). J Agr Biotech. 2012;20(5):521–8.

    Article  Google Scholar 

  35. 35.

    Horsch R, Fry JE, Hoffmann NL, Eichholtz DZ, Rogers SG, Fraley RT. A simple and general method for transferring genes into plants. Science. 1985;227(4691):1229–31.

    CAS  Article  Google Scholar 

  36. 36.

    Bai Y, Pattanaik S, Patra B, Werkman JR, Xie CH, Yuan L. Flavonoid-related basic helix-loop-helix regulators, NtAn1a and NtAn1b, of tobacco have originated from two ancestors and are functionally active. Planta. 2011;234(2):363–75.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Geng X, Li W, Shang H, Gou Q, Zhang F, Zang X, Zeng B, Li J, Wang Y, Ma J, et al. A reference gene set construction using RNA-seq of multiple tissues of Chinese giant salamander, Andrias davidianus. GigaScience. 2017;6(3):gix006.

    CAS  Article  Google Scholar 

  38. 38.

    Wang XW, Luan JB, Li JM, Bao YY, Zhang CX, Liu SS. De novo characterization of a whitefly transcriptome and analysis of its gene expression during development. BMC Genomics. 2010;11(1):400.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Legeai F, Gimenez S, Duvic B, Escoubas JM, Gosselin Grenet AS, Blanc F, Cousserans F, Séninet I, Bretaudeau A, Mutuel D, et al. Establishment and analysis of a reference transcriptome for Spodoptera frugiperda. BMC Genomics. 2014;15(1):704.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Xu Q, He J, Dong J, Hou X, Zhang X. Genomic survey and expression profiling of the MYB gene family in watermelon. Hortic Plant J. 2018;4(1):1–15.

    Article  Google Scholar 

  41. 41.

    Wang F, Suo Y, Wei H, Li M, Xie C, Wang L, Chen X, Zhang Z. Identification and Characterization of 40 Isolated Rehmannia glutinosa MYB Family Genes and Their Expression Profiles in Response to Shading and Continuous Cropping. Int J Mol Sci. 2015:16(7).

  42. 42.

    Zheng X, Yi D, Shao L, Li C. In silico genome-wide identification, phylogeny and expression analysis of the R2R3-MYB gene family in Medicago truncatula. J Integr Agr. 2017;16(7):1576–91.

    CAS  Article  Google Scholar 

  43. 43.

    Gonzalez A, Zhao M, Leavitt JM, Lloyd AM. Regulation of the anthocyanin biosynthetic pathway by the TTG1/bHLH/Myb transcriptional complex in Arabidopsis seedlings. Plant J. 2008;53(5):814–27.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Tominaga-Wada R, Nukumizu Y, Sato S, Kato T, Tabata S, Wada T. Functional divergence of MYB-related genes, WEREWOLF and AtMYB23 in Arabidopsis. Biosci Biotechnol Biochem. 2012;76(5):883–7.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Kirik V, Schnittger A, Radchuk V, Adler K, Hulskamp M, Baumlein H. Ectopic expression of the Arabidopsis AtMYB23 gene induces differentiation of trichome cells. Dev Biol. 2001;235(2):366–77.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Bloomer RH, Juenger TE, Symonds VV. Natural variation in GL1 and its effects on trichome density in Arabidopsis thaliana. Mol Ecol. 2012;21(14):3501–15.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Jung C, Seo JS, Han SW, Koo YJ, Kim CH, Song SI, Nahm BH, Choi YD, Cheong JJ. Overexpression of AtMYB44 enhances Stomatal closure to confer abiotic stress tolerance in transgenic Arabidopsis. Plant Physiol. 2008;146(2):623.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Geekiyanage S, Takase T, Ogura Y, Kiyosue T. Anthocyanin production by over-expression of grape transcription factor gene VlmybA2 in transgenic tobacco and Arabidopsis. Plant Biotechnol Rep. 2007;1(1):11–8.

    Article  Google Scholar 

  49. 49.

    Ban Y, Honda C, Hatsuyama Y, Igarashi M, Bessho H, Moriguchi T. Isolation and functional analysis of a MYB transcription factor gene that is a key regulator for the development of red coloration in apple skin. Plant Cell Physiol. 2007;48(7):958–70.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Allan AC, Hellens RP, Laing WA. MYB transcription factors that colour our fruit. Trends Plant Sci. 2008;13(3):99–102.

    CAS  Article  PubMed  Google Scholar 

Download references


We gratefully acknowledge Huiping Ma and Zhengfeng Peng for helpful discussions.


This work was supported by the National Natural Science Foundation of China (31400602), the Henan Province Science and Technology Breakthrough Project (172102310054), and the Applied Science and Technology Research Fund of Luoyang Normal University (2015-YYJJ-003). The funding body had no contribution in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information




YZ, SX and YC designed the research. YZ, JW, XW and RL performed the experiments. YZ, SX, YC and JH analyzed the data. YZ wrote the manuscript and SX and YC revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Yanzhao Zhang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Table S1. Primers used for qRT-PCR analysis of tree peony and transgenic tobacco samples.

Additional file 2:

Table S2. Primers used for PCR amplification of MYB genes.

Additional file 3:

Table S3. Clean reads quality metrics of sequencing project.

Additional file 4:

Table S4. Unigene annotation in public databases.

Additional file 5:

Table S5. MYB deduced amino acid sequence characteristics and predicted subcellular location of genes.

Additional file 6:

Figure S1. The top sequence matches based on the Nr annotation.

Additional file 7:

Figure S2. Sequence alignment and phylogenetic analysis of PsMYB57 and its orthologs. (a) Alignment of the deduced amino acid sequences of PsMYB57 and other R2R3 MYBs associated with anthocyanin from different plant species. The R2 and R3 domain are indicated above the alignment. The motif [DE]Lx2[RK]x3Lx6Lx3R in the R3 repeat is indicated by dark triangles. The motif KPXPR(S/T) F is shown by a red box. (b) Phylogenetic tree was built using the neighbor-joining method using MEGA 5 software, and bootstrap value was set to 1000. PsMYB57 was marked by a red box. Putative functions of all R2R3 MYBs are listed on the right.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhang, Y., Xu, S., Cheng, Y. et al. Functional identification of PsMYB57 involved in anthocyanin regulation of tree peony. BMC Genet 21, 124 (2020).

Download citation


  • Transcriptome
  • MYB gene family
  • Tree peony
  • PsMYB57
  • Anthocyanin